Exos RDM

RDM – Ossatures
Manuel d’exercices
Yves Debard
Institut Universitaire de Technologie du Mans
Département Génie Mécanique et Productique
http://iut.univ-lemans.fr/ydlogi/index.html
26 juin 2006 – 29 mars 2011

Table des matières
1 Exemples 1
Exemple 1 : Portique plan . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
Exemple 2 : Treillis plan à nœuds articulés . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
Exemple 3 : Anneau plan . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
Exemple 4 : Plancher . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
Exemple 5 : Ossature spatiale . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
Exemple 6 : Modes propres d’un anneau plan . . . . . . . . . . . . . . . . . . . . . . . . . . 12
Exemple 7 : Ossature plane . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
2 Analyse statique 16
E1 : Treillis plan à noeuds articulés . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
E2 : Ossature plane . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
E3 : Ossature plane . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
E4 : Ossature plane . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
E5 : Ossature plane . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
E6 : Poutre droite . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
E7 : Poutre courbe . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
E8 : Ossature plane . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
E9 : Poutre à section droite variable soumise à son poids propre . . . . . . . . . . . . . . . . 26
E10 : Treillis spatial à nœuds articulés . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
E11 : Portique plan – poutre soumise à une variation de température . . . . . . . . . . . . . 29
E12 : Treillis plan – poutre soumise à une variation de température . . . . . . . . . . . . . . 30
E13 : Ossature plane – appui incliné . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31
3 Sections droites : caractéristiques et contraintes 32

S1 : Caractéristiques d’une section droite . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
S2 : Torsion d’une poutre rectangulaire . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
S10 : Contrainte normale dans une section droite : flexion déviée . . . . . . . . . . . . . . . 45
S11 : Contraintes dans une section droite : flexion-torsion . . . . . . . . . . . . . . . . . . . 46
S12 : Cisaillement du à l’effort tranchant . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 48
S13 : Contrainte normale dans une poutre à section droite variable . . . . . . . . . . . . . . 49
S14 : Contrainte normale dans une section droite : flexion déviée . . . . . . . . . . . . . . . 50
S15 : Section droite à parois minces . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51
S16 : Contraintes tangentielles dans un caisson multicellulaire . . . . . . . . . . . . . . . . . 53
3
S17 : Cisaillement dans un profil mince fermé et simplement cloisonné . . . . . . . . . . . . 55
S18 : Flexion - torsion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57
S19 : Contraintes normales dans une poutre à section droite variable . . . . . . . . . . . . . 59
4 Flambement eulérien 60
F1 : Ossature plane . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 60
F2 : Poutre droite . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 62
F3 : Poutre droite à section variable . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 63
F4 : Poutre console – flexion-torsion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 64
F5 : Lame équerre – flexion-torsion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 66
F6 : Lame équerre – flexion-torsion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 68
F7 : Flambement d’un mât vertical sous son poids propre . . . . . . . . . . . . . . . . . . . 71
F8 : Flambement d’une poutre droite . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 72
F9 : Flambement d’un cadre . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 73
5 Modes propres 75
D1 : Treillis plan à nœuds articulés . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 75
D2 : Poutre droite à section variable . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 76
D3 : Vibrations transversales d’une poutre droite bi-encastrée . . . . . . . . . . . . . . . . . 77
D4 : Portique plan . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 78
D5 : Ossature spatiale . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 79
D6 : Ossature plancher . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 80
D7 : Vibrations transversales d’une poutre droite libre . . . . . . . . . . . . . . . . . . . . . 81
D8 : Premier mode propre d’une poutre console avec masses . . . . . . . . . . . . . . . . . . 82
Références 83
Chapitre 1
Exemples
Exemple 1 : Portique plan
Référence : A. Giet, L. Géminard, Résistance des matériaux, tome 2, 1968, pages 148-156.
Données :
La structure plane représentée sur la figure est constituée de deux poutres de même section droite.
Soient A l’aire des sections droites et IZ leur moment quadratique par rapport à l’axe Z. L’ossature
est encastrée en 1 et articulée en 4. Les poutres sont en acier de module de Young E.
Le nœud 2 porte une force de composantes (P, 0, 0).
L’énergie de déformation due à l’effort tranchant est négligée (modèle de Bernoulli).
On donne :
L=2m
A = 16 cm2 , IZ = 135 cm4
E = 200000 MPa
P = 10000 N
2 RDM – Ossatures
Modélisation et calcul :
Les étapes de la modélisation sont :
Fichier
Nouvelle étude
Définir le type de l’ossature
Ossature plane
Entrer les coordonnées des 4 nœuds : (0,0) (0,1) (0,2) (2,2)
Poutres
Créer des poutres définies par leurs nœuds extrémités : 1 − 2 , 2 − 3 , 3 − 4
Sections droites
Section droite quelconque
A = 16 cm2 , IZ = 135 cm4
Liaisons
L’ossature est encastrée en 1 et articulée en 4
Cas de charges
Le nœud 2 porte une charge de composantes (10000, 0, 0) N.
Matériaux
Définir
Module de Young = 200000 MPa
Calculer
Paramètres
Modèle de Bernoulli
Calculer
Analyse statique
Enregistrer les données et lancer le calcul
Résultats
Exploiter les résultats du calcul
Résultats :
– Déplacements nodaux :
u2 = 2.2144 mm , v2 = −0.0017 mm , θ2z = −0.0388˚

u3 = 0.0245 mm , v3 = −0.0033 mm , θ3z = 0.1510˚
θ4z = −0.0754˚
– Actions de liaison :
R1x = −6077.4 N , R1y = 533.4 N , M1z = 3221.6 N.m

R4x = −3922.6 N , R4y = −533.4 N
Remarque : dans la référence, l’énergie de déformation due à l’effort normal est négligée.
Manuel d’exercices 3
Exemple 2 : Treillis plan à nœuds articulés
Référence : A. Giet, L. Géminard, Problèmes de résistance des matériaux, tome 1, 1973, page 52.
Problème :
La structure représentée sur la figure est composée de trois barres articulées entre elles. L’ensemble
est relié à l’extérieur par trois rotules en 2, 3 et 4.
Les trois barres ont la même section droite : carré plein de côté 10 mm.
Les poutres 1 − 2 et 1 − 4 sont en acier :
module de Young = 200000 MPa

coefficient de dilatation = 11 10−6 K−1
La poutre 1 − 3 est en laiton :

coefficient de dilatation = 18 10−6 K−1
Le nœud 1 porte une charge P~ de composantes (0, −10000, 0) N.
L’ossature subit une augmentation de température de 50 K.
Modélisation :
Nouvelle étude
Définir le type de l’ossature : Plane
Définir l’unité de longueur : m
Entrer les coordonnées des nœuds : (0, −0.8) , (−0.6, 0) , (0, 0) , (0.6, 0)
4 RDM – Ossatures
Poutres
Créer des poutres définies par leur nœud origine et leur nœud extrémité
Relaxations
Les trois poutres sont du type rotule-rotule (liaisons intérieures)
Sections droites
Section droite paramétrée
Carré plein de côté 10 mm
Matériaux
Modifier la couleur courante
Attribuer la couleur courante à la poutre 1 − 3 (bouton Élément)
Entrer les caractéristiques de la poutre en laiton (bouton Définir)
module de Young = 100000 MPa , coefficient de dilatation = 18E−6 K−1
Entrer les caractéristiques des poutres en acier ( bouton Définir)
module de Young = 200000 MPa , coefficient de dilatation = 11E−6 K−1
Liaisons
L’ossature est articulée en 2 , 3 et 4
Cas de charges
Le nœud 1 porte une force de composantes (0, −10000, 0) N
Variation de température = 50 K
Calculer
Analyse statique
Résultats :
– Déplacement du nœud 1 :
u1 = 0 , v1 = −0.96 mm
– Allongement des poutres :
∆1−2 = ∆1−4 = 0.768 mm , ∆1−3 = 0.960 mm
– Efforts normaux :
N1−2 = N1−4 = 4370 N , N1−3 = 3008 N
Remarque : pour extraire ces résultats, utiliser le bouton droit de la souris.

Exemple 3 : Anneau plan
Référence : solution analytique.
Données :
L’anneau de plan moyen {O, xy} et de section droite constante (carré plein de coté c) représenté sur
la figure est réalisé en acier de module de Young E et de coefficient de Poisson ν.
Le tronçon 6 − 2 porte une force uniformément répartie d’intensité linéique (0, p, 0).
Le tronçon 5 − 4 porte une force uniformément répartie d’intensité linéique (0, −p, 0).
L’énergie de déformation due à l’effort tranchant est prise en compte (modèle de Timoshenko).
On donne :
E = 200000 MPa , ν = 0.3
c = 10 mm , L = R = 50 mm
p = −10 N/mm
Modélisation :
Le problème présente une symétrie par rapport aux plans x = 0 et y = 0. Il suffit de modéliser le
quart de l’anneau.
Fichier
Bibliothèque
La géométrie existe dans la bibliothèque d’ossatures paramétrées
Ossature plane
Numéro 31 : R = 50 mm , L = 50 mm , l’arc est discrétisé en 20 éléments
6 RDM – Ossatures
Matériau
Définir
E = 200000 MPa , ν = 0.3
Sections droites
Section droite paramétrée
Carré plein de côté c = 10 mm
Liaisons/Symétries
La structure est symétrique par rapport au plan x = 0 : désigner le nœud 1
La structure est symétrique par rapport au plan y = 0 : désigner le nœud 3
Cas de charges
La poutre 1 − 2 une force uniformément répartie d’intensité (0, −10, 0) N/mm
Calculer
Paramètres
Modèle de Timoshenko
Calculer
Analyse statique
Résultats :
Référence :
– Déplacements :
(6 π 2 + 17π − 6) pR4 π pR2 (2 + π) pR2
v1 = + +
24 (2 + π) EIz 4 EA 4 GAky
= −0.324026 − 0.000982 − 0.005013 = −0.330021 mm
(π − 14) pR4 pR2 pR2
u3 = + −
6 (2 + π) EIz 2 EA 2 GAky
= 0.131992 − 0.000625 + 0.001950 = 0.133317 mm
– Actions de liaisons :
(14 + 3 π) pR2
F1x = 0 , M1z = = −18983 N.mm
6 (2 + π)
(2 + 3 π) pR2
F3y = −pR = 500 N , M3z = = −18567 N.mm
3 (2 + π)
– Moment fléchissant dans la section 2 :
4 pR2
Mfz2 = − = 6483 N.mm
3 (2 + π)
– Contraintes normales :
¾
σa (14 + 3 π) pR2
=∓ = ±113.90 MPa
σb (2 + π) c3
¾ ½
σc pR 2 (2 + 3 π) pR2 106.10
= 2 ∓ = MPa
σd c (2 + π) c 3 −116.10
Solution éléments finis :
– Déplacements :
v1 = −0.329765 mm , u3 = 0.133290 mm
F1x = 0 N , M1z = −18977 N.mm , F3y = 500 N , M3z = −18523 N.mm
– Moment fléchissant dans la section 2 : Mfz2 = 6477 N.mm
– Contraintes normales :
σa = 113.86 MPa , σb = −113.86 MPa , σc = 106.14 MPa , σd = −116.14 MPa
Remarque :
Avec le module RDM – Éléments finis (hypothèse contraintes planes, 600 triangles à 6 nœuds), on
obtient :
v1 = −0.328065 mm u3 = 0.133370 mm
σa = 113.96 MPa , σb = −113.96 MPa , σc = 99.66 MPa , σd = −124.20 MPa
La théorie des poutres courbes [3] donne :
σc = 99.10 MPa , σd = −124.00 MPa

8 RDM – Ossatures
Exemple 4 : Plancher
Référence : W. Weawer, J. Gere, Matrix analysis of framed structures, Van Nostrand Reihnold,
1990, pages 342-345.
Problème :
L’ossature plancher représentée sur la figure est constituée de cinq poutres de même section droite.
Les sections 1 , 3 , 5 et 6 sont encastrées.
Le nœud 2 porte une force de composantes (0, 0, 50) kN et un couple de comosantes (0, 100, 0) kN.m.
La poutre 1 − 2 porte en son milieu une force ponctuelle de composantes (0, 0, −150) kN.
La poutre (5 − 4) porte sur toute sa longueur une charge uniformément répartie d’intensité linéique
(0, 0, −75) kN/m.
On donne :
L=2m
module de Young = 200000 MPa , coefficient de Poisson = 0.25
aire = 102 cm2 , constante de torsion de Saint Venant J = 2 105 cm4 , IZ = 105 cm4
P = 5000 daN
Nouvelle étude
Définir le type de l’ossature : Plancher

Entrer les coordonnées des nœuds
Poutres
Sections droites
Section quelconque
Aire = 100 cm2
Constante de torsion de Saint Venant : J = 2E5 cm4
Moment quadratique : IZ = 1E5 cm4
Liaisons
L’ossature est encastrée en 1 , 3 , 5 et 6
Cas de charges
Le nœud 2 porte une force Fz = 50 kN
Le nœud 2 porte un couple My = 100 kN.m
La poutre 1 − 2 porte une force ponctuelle Fz = −150 kN située à 3 m du nœud origine
La poutre 5 − 4 porte une force uniformément répartie fz = −75 kN/m
Matériau
Définir
Module de Young = 200000 MPa , coefficient de Poisson = 0.25
Calculer
Analyse statique
Résultats :
w2 = −1.2182 mm , θ2x = −0.35599 10−3 rad , θ2y = −0.14976 10−3 rad
w4 = −2.0993 mm , θ4x = 0.28856 10−3 rad , θ4y = 0.18376 10−3 rad
F1z = 93.528 kN , M1x = 9.493 kN.m , M1y = −163.092 kN.m
F3z = 34.452 kN , M3x = 14.240 kN.m , M3y = 76.393 kN.m

F5z = 214.940 kN , M5x = −11.543 kN.m , M5y = −239.068 kN.m
F6z = 57.080 kN , M6x = −128.588 kN.m , M6y = −7.351 kN.m
10 RDM – Ossatures
Exemple 5 : Ossature spatiale

Référence : J.-J. Barrau, S. Laroze, Calcul de structures par éléments finis, ENSAE, 1984.
Problème :
L’ossature spatiale représentée sur la figure est constituée de poutres dont les sections droites sont
des rectangles pleins.
Les coordonnées nodales sont :
nœud x (m) y (m) z (m)

1 0 0 0
2 0 0 4
3 0 8 4
4 0 11 4
5 3 8 4
6 3 8 0
Les caractéristiques élastiques du matériau sont : E = 100000 MPa et ν = 0.2987.
Les sections 1 et 6 sont encastrées.
Le nœud 4 porte une force F~ de composantes (0, 0, −1000) daN .
Nouvelle étude
Définir le type de l’ossature : Spatiale
Définir l’unité de longueur : m
Entrer les coordonnées des nœuds
Poutres
Les poutres sont définies par leur nœud origine et leur nœud extrémité
Matériaux
Module de Young = 100000 MPa , coefficient de Poisson = 0.2987
Sections droites
Changer les poutres 3 − 5 et 5 − 6 de groupe
Paramétrée
Désigner la poutre 2 − 3
Rectangle plein : 600 × 300 mm
Paramétrée
Paramétrée
Repère local
Modifier le repère local de la poutre 1 − 2 (angle = 90˚)
Liaisons
L’ossature est encastrée en 1 et 6
Cas de charges
Le nœud 4 porte une charge de composantes (0, 0, −1000) daN
Calculer
Paramètres du calcul
Modèle de Timoshenko
Calculer
Analyse statique
Résultats :
Moments aux extrémités des poutres (en daN.m)
élément Mto MfY o MfZo Mte MfY e MfZe

1−2 Référence -6 271.2 -389.6 -6 322 -104.7
RDM – Ossatures -5.6 271.5 -389.7 -5.64 322.8 -101.2
2−3 Référence 322.2 -6 -104.7 -322.2 96.6 -2513
RDM – Ossatures -322.8 -5.6 -101.2 -323.1 97.04 -2511
3−4 Référence 0 0 -3000 0 0 0
RDM – Ossatures 0 0 -3000 0 0 0
3−5 Référence 487.2 322.2 -96.6 487.2 -3581 117.1
RDM – Ossatures 488.6 322.8 -97.04 488.6 -3581 119.4
5−6 Référence 117.1 -3581 -487.2 117.1 -3632 -202
RDM – Ossatures 119.4 -3581 -488.6 119.5 -3632 -200.1
Exemple 6 : Modes propres d’un anneau plan
Référence : Guide de validation des progiciels de calcul de structures, AFNOR, 1990, page 208.
Problème :
L’ossature plane représentée sur la figure est constituée d’un anneau (centre O, rayon moyen R) et
d’une patte 1 − 2 de longueur L. L’ensemble est encastré en 1.
L’anneau et la patte ont des sections droites rectangulaires pleines.
Soient E le module de Young du matériau et ρ sa masse volumique.
On recherche les six premiers modes propres de cet anneau.
On donne :
R = 0.1 m , L = 0.0275 m
E = 72000 MPa , ρ = 2700 kg/m3
Section droite de l’anneau : Ha = 5 mm , Ba = 10 mm
Section droite de la patte : Hp = 3 mm , Bp = 10 mm
Modélisation :
Bibliothèque (une partie de la géométrie existe dans la bibliothèque d’ossatures paramétrées)

Définir le type d’ossature : Plane
Entrer le numéro de l’ossature paramétrée : 30
Rayon = 0.1 m , angles : 0 et 360 degrés , le cercle est discrétisé en 60 éléments
Poutres (création de la patte)
Ajouter une poutre verticale

Origine : nœud 1 , longueur = 0.0275 m
Matériau
Masse volumique = 2700 kg/m3
Sections droites
Changer la patte de groupe de section
Paramétrée
Désigner l’anneau
Rectangle plein : 5 x 10 mm
Paramétrée
Désigner la patte
Rectangle plein : 3 x 10 mm
Liaisons
La patte est encastrée en 1
Poutres
Discrétiser la patte en 6 éléments
Calculer
Modes propres
6 premiers modes propres
Résultats :
Fréquences en Hz :
Mode Référence RDM – Ossatures

1 28.8 28.81
2 189.3 189.30
3 268.8 268.60
4 641.0 640.52
5 682.0 681.65
6 1063.0 1062.70
Exemple 7 : Ossature plane

Référence : W. Weawer, J. Gere, Matrix analysis of framed structures, Van Nostrand Reihnold,
1990, page 283.
Données :
L’ossature plane représentée sur la figure est constituée de cinq poutres droites identiques articulées
entre elles.
Les caractéristiques de ces poutres sont :
Module de Young : E
Longueur : L
Aire de la section droite A
Les nœuds 1 et 2 sont articulés et le nœud 4 repose sur un appui simple (u4 = 0).
Le nœud 3 porte une force (P, −2 P ).

La poutre 1 − 2 porte en son milieu une force de composantes (0, 2 P ).
La poutre 2 − 4 porte en son milieu une force de composantes (0, −2 P ).
La poutre 3 − 4 porte sur toute sa longueur une charge triangulaire dont l’intensité à l’extrémité 4 a
pour composantes (0, −6 P/L).
On donne :
L = 1.5 m
section droite paramétrée : carré creux, côté extérieur c = 100 mm , t = 5 mm
P = 1000 daN
Modélisation :

Fichier
Nouvelle étude
Définir le type de l’ossature : Plane
Entrer les coordonnées des nœuds 1 et 2 : 0, 0 , 1.5, 0
Nœuds
Créer un nœud défini par un nœud de référence et ses coordonnées polaires :
nœud 3 : nœud de référence = 1 , coordonnées = (60˚, 1.5 m)
nœud 4 : nœud de référence = 2 , coordonnées = (60˚, 1.5 m)
Afficher ⇒ Échelle maximale
Poutres
Relaxations
Toutes les poutres sont du type rotule-rotule
Sections droites
Bibliothèque
Carré creux de côté 100 mm et d’épaisseur 5 mm
Liaisons
L’ossature est articulée en 1 et 2
L’ossature repose sur un appui simple (u = 0) en 4
Charges
Le nœud 3 porte une force de composantes (1000, −2000) daN
La poutre 1 − 2 porte en son milieu une force de composantes (0, 2000) daN
La poutre 3 − 4 porte sur toute sa longueur une charge triangulaire dont l’intensité en 4 est
égale à (0 − 4000) daN/m
La poutre 2 − 4 porte en son milieu une force de composantes (0, −2000) daN
Matériau
Définir
Calculer
Analyse statique
Résultats :
u3 = 0.02632 mm , v3 = −0.07895 mm , v4 = −0.15789 mm
R1x = 699 daN , R1y = 211 daN
R2x = 699 daN , R2y = 4789 daN
R4x = −2398 daN
– Efforts intérieurs sur la poutre 3 − 4 :
N3 = N4 = −667 daN , TY 3 = −1000 daN , TY 4 = 2000 daN

Chapitre 2
Analyse statique
E1 : Treillis plan à noeuds articulés
Référence : F. Frey – Analyse des structures et milieux continus, Presses Polytechniques et Uni-
versitaires Romandes, 1985, page 108.
Problème :
L’ossature plane représentée sur la figure est constituée de 9 poutres droites articulées entre elles.
L’ensemble est lié à l’extérieur par un appui simple en 4 et une rotule en 1.
La structure est en acier de module de Young E = 210000 MPa.
Les poutres sont des carrés creux de côté 70 mm et d’épaisseur 5 mm (bibliothèque).
Le nœud 1 porte une force : ½ ¾

0
{Q} = daN
−1800
Les nœuds 2 et 3 portent une force :
½ ¾
0
{P } = daN
−3600
Résultats :
Actions de liaison :
R1x = 0 , R1y = 5400 daN
R4y = 3600 daN
Efforts normaux :
N12 = N23 = −5143 daN
N34 = −3600 daN
N15 = 6278 daN
N56 = 3758 daN
N64 = 5091 daN
N25 = −3600 daN
N53 = 1883 daN
N63 = −4680 daN
E2 : Ossature plane
Référence : A. Jalil – Calcul pratique des structures, Eyrolles, 1985, page 55.
Problème :
L’ossature plane représentée sur la figure est constituée de 3 poutres droites soudées entre elles.
L’ensemble est lié à l’extérieur par un appui simple en 1 et une articulation en 4.
La structure est en acier.
Les trois poutres sont des HEA 600.
La poutre 1 − 2 porte en son milieu A une force : P~A = (0, −2000) daN.
La poutre 3 − 4 porte en son milieu C une force : P~C = (−1000, 0) daN.
La poutre 2 − 3 porte en son milieu B une force : P~B = (0, −2000) daN et sur le tronçon 2 − B une
charge uniformément répartie ~q = (0, −1000) daN/m.
Résultats :
Les actions de liaison sont : ½ ¾

~1 = 0
R daN
4679
½ ¾
~ 4 = 1000 daN
R
2321
Le moment fléchissant maximal est égal à 18301 daN.m et situé sur la poutre 2 − 3 à X = 2.66 m.
E3 : Ossature plane
Problème :
L’ossature plane représentée sur la figure est constituée de quatre poutres droites. L’ensemble est lié
à l’extérieur par deux rotules en 1 et 5. Les poutres 2 − 3 et 3 − 4 sont liées entre elles par une rotule.
La structure est en acier.
Les quatre poutres sont des HEA 600.

½ ¾
4000
Le noeud 2 porte une force P~ = daN.
0
½ ¾
1000
La poutre 1 − 2 porte une charge uniformément répartie ~q1 = daN/m.
0
½ ¾
0
Les poutres 2 − 3 et 3 − 4 portent une charge uniformément répartie ~q2 = daN/m.
−5000
Résultats :
Les actions de liaison sont : ½ ¾

~1 = −1250
R daN
9583
½ ¾
~ −7750
R5 = daN
20417
Le moment fléchissant est maximal en 4 et Mfmax = 38750 daN.m
E4 : Ossature plane
Problème : l’ossature plane représentée sur la figure est constituée de quatre poutres droites. L’en-
semble est lié à l’extérieur par deux articulations en 1 et 5. Les poutres 2 − 3 et 3 − 4 sont liées entre
elles par une rotule.
La structure est en acier de module de Young 210000 MPa.
Les quatre poutres sont des HEA 600.
Le noeud 2 porte une force F~ = (2000, −5000, 0) daN et un couple C

~ = (0, 0, −3000) daN.m
Les poutres 2 − 3 et 3 − 4 portent une charge uniformément répartie ~q = (0, −1000, 0) daN/m projeté.
Résultats :
Actions de liaison :    
−250 −1750
~ 1 = 6000
R daN , ~5 =
R 5000 daN
   
0 0
E5 : Ossature plane
Référence : W. Weawer, J. Gere – Matrix analysis of framed structures, Van Nostrand Reihnold,
1990, page 228.
Problème :
L’ossature plane représentée sur la figure est constituée de six poutres droites articulées entre elles.
L’ensemble est lié à l’extérieur par deux articulations en 3 et 4.
L’ossature est en acier de module de Young E.
Les caractéristiques des poutres sont :
– poutres 1 − 4 et 3 − 2 : aire = A
– poutres 1 − 2 et 3 − 4 : aire = 0.6 A
– poutres 3 − 1 et 4 − 2 : aire = 0.8 A
La structure porte les charges suivantes :
– le noeud 2 porte une force P~1 de composantes (2P, P, 0).

– la poutre 2 − 4 porte en son milieu une force P~2 de composantes (P, −P, 0).
– la poutre 1 − 2 porte en son milieu un couple C ~ de composantes (0, 0, −1.2 P L).
– la poutre 3 − 1 porte sur toute sa longueur une charge uniformément répartie. La charge par
unité de longueur ~q a pour composantes : (2.5 P/L, 0, 0).
– la poutre 3 − 4 porte en son milieu une force P~3 de composantes (0, −2 P, 0).
Données numériques :
module de Young : E = 200000 MPa
L = 1.25 m
poutres 1 − 4 et 3 − 2 : carré plein de côté 50 mm

poutres 3 − 1 et 4 − 2 : rectangle plein de dimensions 40 × 50 mm
poutres 1 − 2 et 3 − 4 : rectangle plein de dimensions 30 × 50 mm
P = 1000 N
Résultats :
Déplacements :
   
25.00 10−3  26.53 10−3 
nœud 1 : 10.37 10−3 mm , nœud 2 : 10.05 10−3 mm
   
0 0
½ ¾ ½ ¾
~3 = −2890 ~4 = −2110
R N , R N
−5667 7667
Efforts aux extrémités des poutres (en N) :
poutre N à l’origine TY à l’origine N à l’extrémité TY à l’extrémité

1−2 610 2000 610 2000
3−4 0 −1000 0 1000
3−1 4147 - 1000 4147 1000
4−2 −4520 - 500 −3520 500
1−4 −2683 0 −2683 0
3−2 3150 0 3150 0
E6 : Poutre droite
Problème : la poutre droite d’axe x représentée sur la figure est encastrée à ses deux extrémités.
Les caractéristiques de la section droite sont :
aire = 10−3 m2 , IZ = 1.7 10−8 m4
Le module de Young est E = 2.1 1011 Pa.
Elle porte :
 
 0 
– sur toute sa longueur une force uniformément répartie p~ = −24000 N/m.
 
0
   
30000  0 
– au point d’abscisse x = 0.3 m une force P~ = 0 ~ =
N et un couple C 0 N.m.
   
0 −3000
 
 10000 
– au point d’abscisse x = 0.7 m une force Q~ = −20000 N.
 
0
Résultats :
Action de liaison : RAx = −24000 N
Déplacement : v (x = 0.5 m) = −4.90 10−2 m
Forces intérieures : TY (x = 0.5 m) = −540 N , MfZ (x = 0.5 m) = 2800 N.m

E7 : Poutre courbe
Référence : Solution analytique.
Problème :
L’ossature plane représentée sur la figure est constituée

d’une poutre courbe 1−2 de centre O et de rayon moyen R.
La section droite est un carré plein de côté c. La poutre
est encastrée en 1.
Elle porte en 1 une force de composante (0, P, 0).
Les caractéristiques élastiques du matériau sont E et ν.
On donne :
R = 60 mm , c = 30
E = 210000 MPa , ν = 0.28
P = 6000 N.
Modélisation :
Ossature paramétrée 30 (20 éléments, rayon = 60 mm , angle de départ = 0˚, angle de l’arc = 270˚).
Résultats :
Le déplacement vertical du point 2 est :

µ ¶
3 πP R 9π P R3 3 πP R
v2 = + +2 + = 0.0045 + 0.8291 + 0.0138 mm
|4 {z
EA}
|
4
{z
EIZ
}
4 GA k
| {z Z}
effort normal moment fléchissant effort tranchant
E
où A est l’aire de la section droite et IZ son moment quadratique par rapport à Z. G =
2(1 + ν)
est le module d’élasticité transversal. Le dernier terme représente l’influence du cisaillement transverse.
kZ = 5/6 est le coefficient d’aire cisaillée.
On obtient pour v2 (en mm) :
nombre d’éléments Timoshenko Bernoulli

10 0.8285 0.8148
20 0.8426 0.8289
40 0.8462 0.8324
référence 0.8474 0.8336
E8 : Ossature plane
Problème : l’ossature plane représentée sur la figure est constituée de trois poutres droites articulées
entre elles.
Elle est en acier de module de Young E L’ensemble est lié à l’extérieur par trois articulations en 1 et
2 et 3.
Les caractéristiques des poutres sont :
poutre 1 − 4 : rectangle plein 2.5 a × 2 a
poutre 2 − 4 : rectangle plein 1.5 a × 2 a
poutre 3 − 4 : carré plein de côté 2 a
La poutre 1 − 4 porte une charge d’intensité linéique ~q qui lui est perpendiculaire.
On donne :
L = 10 cm , E = 200000 MPa , a = 10 mm , q = 8 N/mm
Résultats :
Déplacements :
3 q L2 2 q L2
u4 = = 6 10−3 mm , v4 = − = −8 10−3 mm
2 E a2 E a2
Efforts normaux :
3
N1−4 = 0 , N2−4 = −2 q L = −1600 N , N4−3 = − q L = −1200 N
2
E9 : Poutre à section droite variable soumise à son poids propre
Problème : la poutre droite de longueur L représentée sur la figure est encastrée en 1.
Soient E et ρ respectivement le module de Young et la masse volumique du matériau. La section

droite est un rond plein dont le diamètre varie linéairement entre les sections 1 et 2. La poutre est
soumise son poids propre. Soit g l’accélération de la pesanteur.
On donne :
E = 200000 MPa , ρ = 8000 kg/m3
g = 10 m/s2
L = 1.2 m , D = 50 mm
Calculer le déplacement vertical et la rotation de la section 2.
Résultats :
– Flèche en 2 :
– Modèle de Bernoulli :
ρ g L4
v2 = − = −0.1105920 mm
3 E D2
– Modèle de Timoshenko :
ρ g L4 2 (1 + ν) L2 ρ g
v2 = − − = −0.1105920 − 0.0005824 = −0.1111744 mm
3 E D2 3 E ky
– Rotation de la section 2 :
ρ g L3
θ2z = − = −0.0079206˚
2 E D2
E10 : Treillis spatial à nœuds articulés

Référence : W. Weawer, J. Gere – Matrix analysis of framed structures, Van Nostrand Reihnold,
1990, page 352.
Problème : la structure représentée sur la figure ci-dessous est constituée de 9 poutres articulées
entre elles.
Les coordonnées des nœuds sont (en m) :
noeud 1 2 3 4 5 6
x 3 0 -3 3 -3 0
y 0 0 0 5 5 5
z 0 3 0 0 0 -3
Soient E = 80000 MPa et ν = 0.3 les caractéristiques élastiques du matériau.
Les caractéristiques (section quelconque) des poutres sont :
4 − 5 , 4 − 6 et 5 − 6 : A = 100 cm2 , J = IY = IZ = 10000 cm4
1 − 4 , 1 − 6 , 3 − 5 , 3 − 6 , 2 − 4 et 2 − 5 : A = 200 cm2 , J = IY = IZ = 20000 cm4
L’ ensemble est fixé au mur par 3 rotules en 1, 2 et 3.
Le nœud 6 porte une force de composantes (48, 24, −24) kN. La poutre 1 − 4 porte en son milieu
une force de composantes (0, 0, −24) kN. La poutre 4 − 5 porte sur toute sa longueur une force
uniformément répartie d’intensité linéique (0, 0, 24) kN/m.
Résultats :
Déplacements nodaux :
nœud 4 : (0.44359, 0.30312, 2.08842) mm
nœud 5 : (0.02059, 0.33437, 1.79382) mm
nœud 6 : (0.41121, 1.34562, 2.10171) mm
nœud 1 : (−13.8, −74, −1.8) kN
nœud 2 : (−6, 204, −122.4) kN
nœud 3 : (−28.2, −154, 28.2 kN

E11 : Portique plan – poutre soumise à une variation de température
Problème : la structure plane représentée sur la figure est constituée de 3 poutres de même matériau
et de même section droite (rond creux de diamètre extérieur d et d’épaisseur t).
La poutre 2 − 3 est articulée en 2 et 3. L’ensemble est encastré en 1 et 4. Soient E et α respectivement

le module de Young et le coefficient de dilatation du matériau.
La poutre 2 − 3 subit une variation de température égale à ∆T .
On donne :
L = 1 m , H = 0.3 m , d = 80 mm , t = 5 mm
E = 210000 MPa , α = 13 10−6 K−1
∆T = 50 K
Résultats :
Soient A et Iz respectivement l’aire et le moment quadratique de la section droite.
L’allongement de la poutre 2 − 3 est égal à :

NL
δ = α ∆T L +
EA
où N est l’effort normal dans la poutre 2 − 3.
L’effort normal N est solution de l’équation :

µ ¶
1 N H3 L H3 1
δ=− soit N + = − α ∆T L
2 3 EIz 2 EA 3 EIz 2
On obtient :
N = −6071.3 N , δ = 0.62546 mm
E12 : Treillis plan – poutre soumise à une variation de température
Problème :
Le treillis plan à nœuds articulés représenté sur la figure

ci-contre est constituée de 5 poutres de même matériau
et de même section droite (carré creux de côté extérieur c
et d’épaisseur t) . Les poutres 1 − 2, 1 − 3 et 1 − 4 ont
la même longueur L. Le triangle 2 − 3 − 4 est équilaté-
ral.
L’ensemble est articulé en 2 et 4.
Soient E et α respectivement le module de Young et le

coefficient de dilatation du matériau.
La poutre 1 − 3 subit une variation de température égale à ∆T .
On donne :
L = 0.5 m , c = 40 mm , t = 5 mm
E = 200000 MPa , α = 12.5 10−6 K−1
∆T = 30 K
Résultats :
Soit A l’aire de la section droite.
L’effort normal dans les poutres 1 − 2, 1 − 3 et 1 − 4 est égal à :

√
− 3 α ∆T E A
N= √ = −12636 N
2+3 3
L’effort normal dans les poutres 2 − 3 et 3 − 4 est égal à :
−N
√ = 7296 N
3
Le déplacement vertical du point 1 est égal à :
2 NL
= −0.09026 mm
EA
L’allongement de la poutre 1 − 3 est égal à :
NL
δ = α ∆T L + = 0.14237 mm
EA
E13 : Ossature plane – appui incliné
Problème : la structure plane représentée sur la figure ci-dessous est constituée de 2 poutres de même
matériau et de même section droite (rond creux de diamètre extérieur d et d’épaisseur t).
Elle est articulée en 1 et repose en 3 sur un appui incliné à 45˚par rapport à l’axe x. Soit E le module
de Young du matériau.
La poutre (2 − 3) porte une charge uniformément répartie d’intensité (0, q, 0).
On donne :
L = 0.3 m , d = 30 mm , t = 5 mm , E = 210000 MPa , q = −1000 N/m
Modélisation : ajouter un changement de repère {x0 , y 0 } en 3, puis définir la liaison dans ce repère
local.
Résultats :
−2 qL
Posons : X = = 200 N
3
Les actions de liaison sont égales à (dans le repère {x, y}) :
½ ¾ ½ ¾ ½ ¾ ½ ¾
~1 = X 200 ~3 = −X −200
R = N , R = N
2X 400 X 200
Le moment fléchissant en 2 est égal à : −XL = −60 N.m
Soient A et Iz respectivement l’aire et le moment quadratique de la section droite. Le déplacement

horizontal du nœud 3 dans le repère {x, y} est égal à :
6 qL4 12 qL2
u3 = + = −0.27009 mm
27 EIz 9 EA
Chapitre 3
Sections droites : caractéristiques et

contraintes
S1 : Caractéristiques d’une section droite

Problème : considérons la section droite représentée sur la figure ci-dessous. Soient G le centre de
gravité et C le centre de torsion.
1. Première étude :
On donne : L = H = 100 mm , t = 20 mm.
Calculer les caractéristiques de la section droite pour plusieurs maillages.
2. Deuxième étude :
Pour t = 5, 10, 20, 30, 40 mm, calculer les caractéristiques de la section et comparer avec les
solutions analytiques valables pour les profils minces.
Modélisation :
Prendre une ossature spatiale quelconque, modéliser la section (section paramétrée) puis entrer dans
le menu Calculer section droite.
Résultats :
Pour éditer les caractéristiques, sélectionner la commande Caractéristiques du menu Fichier.

On obtient (la valeur en % représente l’écart avec la valeur obtenue avec le maillage le plus fin) :
Maillage J (cm4 ) Iω (cm6 ) YC ( mm ) kY kZ

100 TR3 72.97 6379 -64.78 0.6331 0.2408
6.77 % 2.95 % 0.89 % 2.81 % 2.99 %
400 TR3 70.13 6511 -65.10 0.6219 0.2366
2.62 % 0.94 % 0.40 % 0.99 % 1.20 %
20 TR6 71.34 6535 -65.01 0.6228 0.2378
4.39 % 0.58 % 0.54 % 1.14 % 1.71 %
100 TR6 68.80 6563 -65.26 0.6174 0.2348
0.67 % 0.15 % 0.15 % 0.26 % 0.43 %
150 TR6 68.73 6563 -65.26 0.6174 0.2347
0.57 % 0.15 % 0.15 % 0.26 % 0.38 %
400 TR6 68.45 6570 -65.32 0.6163 0.2341
0.16 % 0.05 % 0.06 % 0.08 % 0.13 %
2200 TR6 68.34 6573 -65.36 0.6158 0.2338
Les formules de résistance des matériaux (R.D.M.) valables pour les profils minces sont donnés
dans les références [1, 2, 4] :
t3 h2 b3 t (2 h + 3 b) 3 b2 b2
J= (h + 2 b) , Iω = , YC = − −
3 12 (h + 6 b) h + 6b h + 2b
où h = L − t et b = H − t/2
On obtient (M.E.F. = solution éléments finis obtenue avec ∼ 400 triangles à 6 nœuds) :
t (mm) J (cm4 ) Iω (cm6 ) YC (mm)

M.E.F. R.D.M. ∆ (%) M.E.F. R.D.M. ∆ (%) M.E.F. R.D.M. ∆ (%)
5 1.206 1.208 0.17 2500 2473 1.08 -74.45 -74.71 0.36
10 9.286 9.333 0.51 4269 4077 4.50 -72.12 -73.25 1.57
20 68.44 69.33 1.29 6569 5393 17.91 -65.32 -70.30 7.70
30 211.3 216.0 2.22 7835 5123 34.61 -55.17 -67.47 22.29
40 454.9 469.3 3.17 7653 4096 46.48 -41.40 -64.64 56.16
∆ représente l’écart entre la solution analytique et la solution éléments finis, cette dernière
servant de référence.
S2 : Torsion d’une poutre rectangulaire

Référence : S. Laroze, Mécanique des structures – Tome 2 : Poutres, Cépaduès, 2005, page 93.
Problème : la poutre console représentée sur la figure est en acier de caractéristiques élastiques E
et ν. Son extrémité libre est soumise à un couple de composantes (0, C, 0).
On donne :
E = 200000 MPa , ν = 0.3 , L = 1 m , a = 100 mm , C = 100 kN.m
Calculer la constante de torsion de la section droite, la rotation θ de l’extrémité libre de la poutre et

le cisaillement maximal τmax pour plusieurs maillages de la section.
Modélisation : activer le menu Calculer section droite du menu Modéliser.
Résultats :
Référence :  
∞
X
64 1 3nπ CL
J = a4 1 − tanh , θ=
π5 n 5 2 GJ
n=1,3,...
 
∞
X
aC  8  1
τmax = 1 − 2 
J π 3 n π
n=1,3,... n2 cosh
2
On obtient (activer le menu Contraintes sur section droite du menu Résultats) :
maillage J (cm4 ) θ (rad) τmax (MPa)

100 TR3 8040.34 0.0161685 114.77
400 TR3 7934.79 0.0163835 119.59
50 TR6 7913.64 0.0164273 124.55
100 TR6 7902.96 0.0164495 124.64
400 TR6 7899.86 0.0164560 124.73
référence 7899.51 0.0164567 124.75
Remarque : le cisaillement est maximal en A et B.

Problème :
Considérons la section droite représentée sur la figure ci-dessous :
On donne : H = 120 mm , L = 100 mm , t = 20 mm
Modélisation :
Prendre une ossature spatiale quelconque, modéliser la section (section paramétrée) puis activer le
menu Calculer section droite.
Résultats :
Maillage J (cm4 ) Iω (cm6 ) kY kZ

100 TR3 84.15 7766 0.3745 0.6768
8.78 % 2.61 % 3.14 % 4.33 %
400 TR3 80.20 7883 0.3676 0.6591
3.67 % 1.14 % 1.24 % 1.60 %
100 TR6 78.22 7946 0.3646 0.6515
1.11 % 0.35 % 0.41 % 0.43 %
200 TR6 77.94 7953 0.3642 0.6508
0.75 % 0.26 % 0.30 % 0.32 %
400 TR6 77.63 7963 0.3637 0.6498
0.35 % 0.14 % 0.17 % 0.17 %
800 TR6 77.51 7967 0.3634 0.6493
0.19 % 0.09 % 0.08 % 0.09 %
4552 TR6 77.36 7974 0.3631 0.6487
t3 (H − t)2 L3 t
J= (H − t + 2 L) Iω =
3 24
t (mm) J (cm4 ) Iω (cm6 )

M.E.F. R.D.M. ∆ (%) M.E.F. R.D.M. ∆ (%)
5 1.306 1.313 0.53 2749 2755 0.22
10 10.206 10.333 1.24 4991 5042 1.02
20 77.627 80.000 3.06 7963 8333 4.65
30 247.74 261.00 5.35 8990 10125 12.63
40 546.82 597.33 9.24 8285 10667 28.75
Problème :
Considérons la section droite représentée sur la figure ci-dessous :
On donne : H = 250 mm L = 100 mm , t = 20 mm.
On donne : H = 250 mm , L = 100 mm.
Modélisation :
Prendre une ossature spatiale quelconque, modéliser la section (section paramétrée puis activer le
Résultats :
Maillage J ( cm4 ) Iω ( cm6 ) YC (mm) kY kZ

100 TR3 473.1 125622 -59.24 0.3398 0.5143
19.71 % 5.08 % 0.49 % 2.81 % 3.86 %
400 TR3 416.7 130683 -59.36 0.3341 0.5018
5.44 % 1.25 % 0.29 % 1.18 % 1.33 %
100 TR6 401.5 131913 -59.48 0.3321 0.4973
1.59 % 0.32 % 0.08 % 0.58 % 0.42 %
200 TR6 398.6 132081 -59.48 0.3315 0.4969
0.86 % 0.20 % 0.08 % 0.39 % 0.34 %
400 TR6 396.6 132226 -59.51 0.3307 0.4958
0.35 % 0.09 % 0.03 % 0.15 % 0.12 %
800 TR6 395.9 132279 -59.52 0.3305 0.4955
0.18 % 0.05 % 0.02 % 0.09 % 0.06 %
2500 TR6 395.2 132342 -59.53 0.3302 0.4952
dans la référence [4] :
t3 h2 L3 t
h = H − 1.75 t J= (h + 14.75 L) Iω =
3 7

M.E.F. R.D.M. ∆ (%) M.E.F. R.D.M. ∆ (%)
5 6.96 7.25 4.17 41581 41573 0.02
10 53.7 56.9 5.96 77286 77223 0.08
20 397 451 13.60 132226 132071 0.12
30 1219 1505 23.46 166662 167170 0.30
40 2599 3531 35.86 182568 185143 1.41
Problème :
Considérons la section droite représentée sur la figure (UPN 400).
On donne (en mm) :
h = 400 , b = 110
tw = 14 , tf = 18
r = 18 , r1 = 9
Calculer les caractéristiques de la section droite

pour plusieurs maillages.
Modélisation :
Prendre une ossature spatiale quelconque, modéliser la section (bibliothèque) puis activer le menu
Calculer section droite.
Résultats :
maillage J (cm4 ) Iω (cm6 ) YC (mm) kY kZ

100 TR3 101.25 195071 -48.62 0.2353 0.5866
25.54 % 8.87 % 4.14 % 22.62 % 3.15 %
400 TR3 84.69 211369 -50.49 0.1975 0.5702
5.01 % 1.26 % 0.45 % 2.92 % 0.26 %
100 TR6 81.00 207173 -49.90 0.2058 0.5760
0.43 % 3.22 % 1.62 % 7.24 % 1.28 %
200 TR6 80.75 213908 -50.72 0.1922 0.5687
0.12 % 0.07 % 0.00 % 0.16 % 0.00 %
400 TR6 80.72 214076 -50.74 0.1919 0.5685
0.09 % 0.00 % 0.04 % 0.00 % 0.04 %
800 TR6 80.70 214109 -50.74 0.1918 0.5685
0.06 % 0.02 % 0.04 % 0.05 % 0.04 %
2600 TR6 80.65 214066 -50.72 0.1919 0.5687
Problème :
Considérons la section droite représentée sur la figure (IPN 500).
On donne (en mm) :
h = 500 , b = 185
tw = 18 , tf = 27
r = 18 , r1 = 10.8
Calculer les caractéristiques de la section droite pour plu-

sieurs maillages.
Modélisation :
Résultats :

100 TR3 511.26 1164816 0.4990 0.6546
36.72 % 11.49 % 2.72 % 26.13 %
400 TR3 398.11 1298520 0.4866 0.5338
6.46 % 1.33 % 0.16 % 2.85 %
100 TR6 402.15 1270004 0.4841 0.5556
7.54 % 3.5 % 0.35 % 7.05 %
200 TR6 376.44 1319475 0.4841 0.5171
0.67 % 0.26 % 0.35 % 0.37 %
400 TR6 375.95 1320964 0.4840 0.5161
0.54 % 0.37 % 0.37 % 0.56 %
800 TR6 375.44 1320812 0.4841 0.5160
0.40 % 0.36 % 0.35 % 0.58 %
3200 TR6 373.94 1316060 0.4858 0.5190
Problème : considérons la section droite représentée sur la figure : HEM 320 .
On donne (en mm) :
h = 359 , b = 309
tw = 21 , tf = 40
r = 27

sieurs maillages.
Modélisation :
Résultats :

100 TR3 1754.81 4813570 0.2356 0.7204
16.20 % 1.56 % 1.82 % 3.02 %
400 TR3 1607.62 4857678 0.2335 0.7082
6.46 % 0.66 % 0.91 % 1.27 %
100 TR6 1521.79 4888856 0.2315 0.6997
0.77 % 0.02 % 0.04 % 0.06 %
200 TR6 1514.09 4888907 0.2315 0.6996
0.26 % 0.02 % 0.04 % 0.04 %
400 TR6 1511.22 4889701 0.2314 0.6994
0.07 % 0.01 % 0.00 % 0.01 %
800 TR6 1510.42 4889938 0.2314 0.6993
0.02 % 0.00 % 0.00 % 0.00 %
3600 TR6 1510.13 4890017 0.2314 0.6993
Problème : considérons la section droite représentée sur la figure (cornière à ailes inégales et à coins
arrondis : [ 70 x 50 x 7 ]).
On donne (en mm) :
a = 70 , b = 50
t=7,r=7
r1 = 3.5

sieurs maillages.
Modélisation :
Résultats :
Maillage J (cm4 ) Iω (cm6 ) βω (cm6 ) YC (mm) ZC (mm) kY kZ

200 TR3 1.4592 2.9307 0.7301 -11.99 15.93 0.5171 0.3891
6.73 % 15.45% 21.05 % 1.64 % 1.30 % 2.34 % 3.26 %
400 TR3 1.4175 3.2488 0.8591 -12.07 16.02 0.5119 0.3836
3.68 % 6.27 % 7.10 % 0.98 % 0.74% 1.31 % 1.80 %
100 TR6 1.3686 3.4656 0.9125 -12.17 16.12 0.5059 0.3775
0.10 % 0.02 % 1.33 % 0.16 % 0.12% 0.12 % 0.19 %
200 TR6 1.3677 3.4659 0.9278 -12.19 16.14 0.5053 0.3768
0.04 % 0.01 % 0.32 % 0.00 % 0.00% 0.00 % 0.00 %
400 TR6 1.3673 3.4661 0.9243 -12.19 16.14 0.5053 0.3769
0.01 % 0.01 % 0.05 % 0.00 % 0.00% 0.00 % 0.03 %
800 TR6 1.3672 3.4663 0.9249 -12.19 16.14 0.5053 0.3768
0.00 % 0.00 % 0.01 % 0.00 % 0.00% 0.00 % 0.00 %
1360 TR6 1.3672 3.4663 0.9248 -12.19 16.14 0.5053 0.3768

Problème : considérons la section droite représentée sur la figure ci-dessous :
On donne : H = 200 mm , B = 120 mm , t = 20 mm.
Modélisation :
Prendre une ossature spatiale quelconque, modéliser la section (section paramétrée) puis activer le
Résultats :
Maillage J (cm4 ) Iω (cm6 ) βω (cm6 ) kY kZ

100 TR3 114.35 83705 131793 0.4772 0.4655
8.14 % 1.13 % 0.41 % 1.53 % 1.39 %
400 TR3 110.81 84186 132030 0.4732 0.4614
4.79 % 0.56 % 0.23 % 0.68 % 0.50 %
100 TR6 106.43 84566 132241 0.4710 0.4593
0.65 % 0.11 % 0.07 % 0.21 % 0.04 %
200 TR6 106.08 84590 132268 0.4707 0.4592
0.32 % 0.08 % 0.05 % 0.15 % 0.02 %
400 TR6 105.85 84627 132303 0.4703 0.4592
0.10 % 0.04 % 0.02 % 0.06 % 0.02 %
700 TR6 105.75 84655 132330 0.4700 0.4591
0.01 % 0.00 % 0.00 % 0.00 % 0.02 %
1240 TR6 105.74 84658 132333 0.4700 0.4591
t3 h2 b3 t (b + 2 h)
h=H −t b = B − 0.5 t J= (h + 2 b) Iω =
3 12 (2 b + h)

M.E.F. R.D.M. ∆ (%) M.E.F. R.D.M. ∆ (%)
5 1.791 1.792 0.06 30337 30335 0.01
10 13.96 14.00 0.29 53944 53923 0.04
20 105.8 106.7 0.85 84627 84452 0.21
30 337.5 342.0 1.33 98495 97945 0.56
40 753.0 768.0 1.99 100656 99556 1.09
50 1379 1417 2.76 95008 93381 1.71
60 2218 2304 3.88 84668 82605 2.44
S10 : Contrainte normale dans une section droite : flexion déviée
Problème : la poutre console représentée sur la figure ci-dessous est soumise en son extrémité libre
à une force de composantes (P, 0, 3 P ).
On donne :
L = 1 m , a = 100 mm , P = −10000 N
Étudier la contrainte normale dans la section encastrée.
Modélisation : sélectionner l’option Ossature spatiale.
Résultats :
Solution analytique : dans la section encastrée, la contrainte normale est égale à :

MfZ MfY
σ=− Y + Z
IZ IY
avec
a4 2 a4
IY = , IZ = , MfY = P L , MfZ = 3 P L
6 3
soit
−3 P L
σ= (3 Y − 4 Z)
2 a4
−15 P L
La contrainte de traction est maximale en A(a, −a/2) : σT = = 75 MPa .
2 a3
La contrainte de compression est maximale en B(−a, a/2) : σC = −75 MPa .
Méthode des éléments finis : pour extraire les résultats ci-dessus, activer le menu Contraintes
sur section du menu Résultats, désigner la poutre, puis entrer l’abscisse de la section encastrée
(commande Abscisse de la section du menu Modéliser).
S11 : Contraintes dans une section droite : flexion-torsion

Problème :
La structure représentée sur la figure ci-dessous est composée de deux poutres de section droite carrée
(côté c). Elle est encastrée en 1 et porte en 3 une force de composantes (0, 0, P ). Soient E et ν les
caractéristiques élastiques du matériau.
On donne :
E = 200000 MPa , ν = 0.3

L = 0.5 m , H = 0.4 m
c = 40 mm
P = −3000 N
Pour plusieurs maillages de la section droite, calculer :
– le déplacement vertical des nœuds 2 et 3.
– dans la section encastrée, la contrainte de cisaillement et la contrainte équivalente de Von Mises

au point M .
– dans la section encastrée, la position et la valeur du cisaillement maximal.
Modélisation :
Sélectionner l’option Ossature spatiale ou l’option Ossature plancher.
Résultats :
Solution analytique :
Les caractéristiques de la section sont :
c4
A = c2 , IZ = , J = 0.1405770 c4 , kY = 5/6
12
– le déplacement vertical du nœud 2 est égal à :
P L3 PL
w2 = + = −2.92969 − 0.01463 = −2.94431 mm
3 EIZ GA kY
– le déplacement vertical du nœud 3 est égal à :
µ ¶
P (L3 + H 3 ) P H 2 L P (L + H)
w3 = + + = −13.09931 − 0.02633 = −13.12564 mm
3 EIZ GJ GA kY
– en M et dans le repère {XY Z}, le tenseur des contraintes a pour expression :

 
σ 0 τ
6 LP 4.80387 HP
[σ(M )] =  0 0 0 avec σ = − 3 , τ=
c c3
τ 0 0
On obtient donc : σ = 140.63 MPa et τ = −90.07 MPa.

√
On en déduit la contrainte équivalente de Von Mises : σV M = σ 2 + 3 τ 2 = 210.03 MPa
– Le cisaillement est maximal en N . En ce point, le tenseur des contraintes a pour expression :

 
0 τ 0
 4.80387 HP 3P
[σ(N )] = τ 0 0 avec τ = + 2
c3 2c
0 0 0
Le cisaillement maximal est donc égal à : τmax = | τ | = 92.89 MPa.
Remarque : le deuxième terme de l’expression ci-dessus est dû à l’effort tranchant.
Méthode des éléments finis :
– Déplacements (utiliser le bouton droit de la souris) :
modèle de Bernoulli modèle de Timoshenko

maillage w2 (mm) w3 (mm) w2 (mm) w3 (mm)
50 TR3 −2.92969 −12.91785 −2.94308 −12.94196
400 TR3 = −13.05862 −2.94414 −13.08464
50 TR6 = −13.08275 −2.94431 −13.10906
400 TR6 = −13.09886 = −13.12519
référence −2.92969 −13.09931 −2.94431 −13.12564
– Contraintes : pour extraire les résultats demandés, activer le menu Contraintes sur section
du menu Résultats, désigner la poutre 1 − 2, puis entrer l’abscisse de la section encastrée
(commande Abscisse de la section du menu Modéliser).
maillage τmax (MPa) σVM (MPa)

100 TR3 86.27 205.20
400 TR3 89.93 206.59
100 TR6 92.49 209.66
400 TR6 92.77 209.89
référence 92.89 210.03
S12 : Cisaillement du à l’effort tranchant

Référence : théorie élémentaire du cisaillement.
Problème :
La poutre console représentée sur la figure est constituée d’une demi-poutrelle IPE 500.
L’extrémité de la poutre est soumise à une force de composantes (0, 0 − 100) kN.
Calculer le cisaillement au centre de gravité G et le cisaillement maximal pour plusieurs maillages.
Résultats :
Activer le menu Contraintes sur section du menu Résultats.
Solution analytique : au centre de gravité de la section, le cisaillement dû à l’effort tranchant TY

est égal à (théorie élémentaire du cisaillement) :
TY MZ
τG =
IZ tw
où MZ est le moment statique par rapport à l’axe Z de la surface de la section située au dessus de
l’axe Z.
Les caractéristiques de la section sont :

189.93
IZ = 3262.67 cm4 , MZ = 189.93 × 10.2 = 183.98 cm3 , tw = 10.2 mm
2
d’où τG = 55.28 MPa.
Méthode des éléments finis : on obtient (pour extraire la quantité τG , effectuer une coupe droite
parallèle à Y au voisinage de G) :
maillage τmax (MPa) τG (MPa)

100 TR3 55.53 55.53
400 TR3 55.35 55.34
100 TR6 57.02 55.58
400 TR6 56.90 55.28
S13 : Contrainte normale dans une poutre à section droite variable
Problème :
La poutre droite de longueur L représentée sur la figure est encastrée en 1. Soit E le module de
Young du matériau. La section droite est un rond plein dont le diamètre varie linéairement entre les
sections 1 et 2. La poutre porte en 2 une force de composantes (0, P, 0).
On donne :
L = 1 m , D2 = 50 mm , D1 = 2 D2 = 100 mm , P = −10000 N

Étudier la contrainte normale le long de la poutre.
Résultat :
Solution analytique : la contrainte normale maximale dans la section droite d’abscisse x est égale
à :
32 P (L − x)
σ=±
π D3
où le diamètre D de la poutre est égal à :
x ³ x´
D = D1 + (D2 − D1 ) = D2 2 −
L L
Cette contrainte normale est maximale dans la section d’abscisse L/2 :
128 P L
σmax = ± = ±120.72 MPa
27 π D23
Méthode des éléments finis : activer le menu Poutre du menu Résultats :
σmax = ±120.72 MPa à x = 0.5 m

S14 : Contrainte normale dans une section droite : flexion déviée

Problème :
La poutre droite de longueur L représentée sur la figure est encastrée en 1. La section droite est une
cornière à ailes inégales (grande aile : a, petite aile : b, épaisseur : t).
La poutre porte en 2 une force de composantes (0, 0, P ).
On donne : L = 0.8 m , a = 100 mm , b = 60 mm , t = 9 mm , P = −1000 N .
Étudier la contrainte normale dans la section encastrée.
Modélisation :
– Sélectionner l’option Ossature spatiale.
– Pour créer la poutre, désigner le point B puis le point A.
– Modifier l’orientation angulaire de la poutre.
Résultats :
Solution analytique : les composantes de la charge dans le repère {XY Z} sont (0, −P cos α, P sin α).
Dans la section encastrée, la contrainte normale est donc égale à :
MfZ MfY P L cos α P L sin α
σ=− Y + Z= Y − Z
IZ IY IZ IY
avec : IY = 22.9774 cm4 et IZ = 153.1764 cm4 .
La contrainte de traction est maximale en M (−27.629 mm , 25.498 mm) : σT = 43.65 MPa .
La contrainte de compression est maximale en N (63.412 mm , −16.841 mm) : σC = −51.02 MPa .
Méthode des éléments finis : pour extraire ces résultats, activer le menu Contraintes sur section
du menu Résultats, puis entrer l’abscisse de la section encastrée (commande Abscisse de la section
du menu Modéliser).
S15 : Section droite à parois minces
Référence : A. Bazergui, T. Bui-Quoc, A. Biron, G. McIntyre, C. Laberge, Résistance des

matériaux, Recueil de problèmes, tome 1, exercice 16.7, Éditions de l’École Polytechnique de Montréal.
Problème :
Les deux sections droites à parois minces représentées sur la figure ci-dessous ont une épaisseur
constante t.
On donne :
a = 100 mm , t = 10 mm.
Calculer pour chaque section droite :
– les caractéristiques en utilisant plusieurs maillages.
– la contrainte moyenne de cisaillement dans la paroi en A et B quand la section droite est soumise
à un moment de torsion MX = 10000 N.m
Modélisation :
Modéliser une poutre console spatiale soumise à un moment de torsion.
Définir la section (géométrie importée : fichier IGES ou .geo).

Résultats :
Pour évaluer les caractéristiques, activer le menu Calculer section droite du menu Modéliser.
Pour évaluer les contraintes, activer le menu Contraintes sur section droite du menu Résultats.
Section I (section ouverte) :
Référence :
J = 2 a t (a2 + t2 ) = 2020 cm4
MX MX
τA = = 24.75 MPa , τB = = 4.95 MPa
4 t (a2 + t2 ) 2 a (a2 + t2 )
On obtient :
– caractéristiques :
Maillage J (cm4 ) Iω (cm6 ) ZC (mm) kY kZ

270 TR3 2190 2015578 −248.69 0.1184 0.4238
1000 TR3 2153 2021152 −248.69 0.1175 0.4204
160 TR6 2140 2022219 −248.67 0.1172 0.4192
220 TR6 2137 2022294 −248.66 0.1171 0.4188
900 TR6 2131 2022509 −248.64 0.1169 0.4179
1800 TR6 2129 2022511 −248.63 0.1168 0.4176
– contraintes moyennes : τA = 24.54 MPa , τB = 4.70 MPa (∼ 900 triangles à 6 nœuds)
Section II (section fermée) :
Référence :
J = 10 a3 t = 10000 cm4
MX MX
τA = 2
= 5.00 MPa , τB = = 10.00 MPa
20 a t 10 a2 t
On obtient :
– caractéristiques :
Maillage J (cm4 ) Iω (cm6 ) ZC (mm) kY kZ

240 TR3 10577 578053 −66.06 0.3281 0.3735
850 TR3 10480 578457 −66.06 0.3249 0.3699
170 TR6 10427 579207 −66.00 0.3237 0.3684
240 TR6 10414 579239 −66.01 0.3231 0.3678
850 TR6 10388 579895 −66.07 0.3221 0.3667
1900 TR6 10379 579995 −66.09 0.3217 0.3663
– contraintes moyennes : τA = 5.10 MPa , τB = 9.90 MPa (∼ 900 triangles à 6 nœuds )

S16 : Contraintes tangentielles dans un caisson multicellulaire
Référence : Résistance des matériaux : cisaillement dans les poutres à parois minces.
Problème :
Considérons la poutre console dont la section droite (caisson rectangulaire à deux cloisons) est repré-
sentée ci-dessous. Les parois et les cloisons ont la même épaisseur t.
On donne : a = 500 mm , t = 20 mm.
La section droite est soumise à un moment de torsion MX = 1000 kN.m .
Évaluer le cisaillement moyen dans la paroi en A, B et C.
3. Troisième étude :
La section droite est soumise à un effort tranchant TY = 1000 kN .
Évaluer le cisaillement moyen dans la paroi en en A, B et C.
Modélisation :
Modéliser une poutre console spatiale.
La section droite est paramétrée : [ 1000 , 2000 , 1000 , 20 , 20 , 20 ] mm

Résultats :
1. Première étude : caractéristiques
Référence :
112 3
J= a t = 5600 103 cm4
5
On obtient :
Maillage J (cm4 ) kY kZ
100 TR3 5752 103 0.4385 0.4550
400 TR3 5683 103 0.4337 0.4502
100 TR6 5664 103 0.4322 0.4486
400 TR6 5650 103 0.4315 0.4475
1000 TR6 5648 103 0.4313 0.4473
2. Deuxième étude : cisaillement dû au couple de torsion Mx
Référence :
MX
τA = 4 τ τC = 3 τ , τB = τ avec τ =
56 a2 t
On obtient (∼ 1000 triangles à 6 nœuds) :
référence RDM-Ossatures
τA 14.29 MPa 14.33 MPa
τB 3.57 MPa 3.61 MPa
τC 10.71 MPa 10.76 MPa
3. Troisième étude : cisaillement dû à l’effort tranchant TY
Référence :
TY
τA = 0 τB = 5 τ τC = 4 τ avec τ=
32 at
référence RDM – Ossatures

τB 15.63 MPa 15.63 MPa
τC 12.50 MPa 12.52 MPa
S17 : Cisaillement dans un profil mince fermé et simplement cloisonné
Référence : S. Laroze, Mécanique des structures, tome 5, Cépaduès, 2005, pages 133,167.
Problème : considérons le gouvernail de profondeur dont la section droite est représentée ci-dessous.
On donne : L = 600 mm , R = 75 mm , t1 = 2 mm , t2 = 4 mm , t3 = 3 mm.
Évaluer les caractéristiques de la section : constante de torsion de Saint Venant J, constante de

gauchissement Iω , position du centre de cisaillement C, coefficients d’aire cisaillée (kY , kZ ).
Évaluer le cisaillement moyen dans les parois 1, 2 et 3 quand la section droite est soumise à un
couple de torsion Mt = 10000 N.m.
3. Troisième étude :
La section droite est soumise à un effort tranchant TY = 10000 N. Évaluer le cisaillement maxi-
mal.
Résultats :
1. Première étude : caractéristiques de la section
Référence : J = 1809 cm4 .
On obtient :
Maillage J (cm4 ) Iω (cm6 ) YC (mm) kY kZ

200 TR3 1828 531 -200.6 0.6166 0.0958
400 TR3 1820 297 -201.1 0.6146 0.0949
200 TR6 1819 370 -201.4 0.6140 0.0946
400 TR6 1816 343 -201.2 0.6136 0.0943
1200 TR6 1815 334 -201.2 0.6134 0.0942
2. Deuxième étude : cisaillement dû au couple de torsion
référence RDM - Ossatures

τ1 42.90 MPa 43.24 MPa
τ2 32.24 MPa 32.16 MPa
τ3 14.39 MPa 14.06 MPa
3. Troisième étude : cisaillement dû à l’effort tranchant TY
Référence : le cisaillement est maximal en A et B et vaut : τmax = 5.16 MPa.
On obtient (maillage : ∼ 1000 triangles à 6 nœuds) : τmax = 5.23 MPa.

S18 : Flexion - torsion
Problème : la structure représentée sur la figure est constituée de deux poutres identiques de lon-
gueur L. Soient E et ν les caractéristiques élastiques du matériau. L’ensemble est encastré en 1. Les
deux poutres portent une charge uniformément répartie d’intensité linéique (0, 0, p).
On donne :
E = 200000 MPa , ν = 0.3
L = 0.6 m , section droite : IPN 180
p = −1000 N/m
Dans chacun des cas suivants et pour plusieurs maillages de la section droite, évaluer le déplacement
vertical des sections 2 et 3.
Résultats :
Les caractéristiques J, kY . . . utilisées dans la solution de référence sont extraites de la bibliothèque

de profilés (maillage = 4 × 1993 triangles à 6 nœuds).
Cas 1 :
Référence :
11 pL4 3 pL2
w2 = + = −0.02057 − 0.00582 = −0.02638mm
24 EIZ 2 GA kY
7 pL4 2 pL2 pL4
w3 = + + = −0.02618 − 0.00775 − 9.28441 = −9.31834 mm
12 EIZ GA kY 2 GJ
Les déplacements obtenus sont (en mm) :

maillage w2 w3 w2 w3
400 TR3 −0.02057 −8.46166 −0.02635 −8.46937
800 TR3 = −8.83228 −0.02638 −8.84003
400 TR6 = −9.25606 −0.02640 −9.26384
800 TR6 = −9.26805 −0.02641 −9.27583
référence −0.02057 −9.31059 −0.02638 −9.31834
Cas 2 :
Référence :
11 pL4 3 pL2
w2 = + = −0.36576 − 0.00444 = −0.37020 mm
24 EIY 2 GA kZ
7 pL4 2 pL2 pL4

w3 = + + = −0.46551 − 0.00592 − 9.28441 = −9.75584 mm
12 EIY GA kZ 2 GJ
Les déplacements obtenus sont (en mm) :

maillage w2 w3 w2 w3
400 TR3 −0.36576 −8.90356 −0.37000 −8.90921
800 TR3 = −9.27178 −0.37071 −9.27759
400 TR6 = −9.69548 −0.37022 −9.70143
800 TR6 = −9.70751 −0.37023 −9.71347
référence −0.36576 −9.74992 −0.37020 −9.75584
S19 : Contraintes normales dans une poutre à section droite variable
Problème :
La poutre droite de longueur L représentée sur la figure est encastrée en 1. La section droite est un
carré plein dont le côté varie linéairement entre les sections 1 et 2.
La poutre est soumise en 2 à :
– une force de composantes (N, F, 0).

– un couple de composantes (0, 0, C).
On donne :
L = 1 m , c = 10 mm
N = 1000 N , F = 1 N , C = 1 N.m
Évaluer les contraintes normales dans les sections 1 et 2.
Résultats :
Solution analytique :
µ ¶ µ ¶
1 3 (C + F L) 1 3 (C + F L)
σ1,inf = 2 N + = 4 MPa , σ1,sup = 2 N − = 1 MPa
4c c 4c c
µ ¶ µ ¶
1 6C 1 6C
σ2,inf = 2 N + = 16 MPa , σ2,sup = 2 N − = 4 MPa
c c c c
Méthode des éléments finis : pour extraire ces résultats, activer le menu Poutre du menu Ré-
sultats, puis entrer l’abscisse de la section (commande Valeur en un point du menu Résultats).
Chapitre 4
Flambement eulérien
F1 : Ossature plane
Référence : S.P. Timoshenko, J.M. Gere, Théorie de la stabilité élastique, Dunod, 1966, page 69.
Problème :
L’ossature plane représentée sur la figure est constituée de trois poutres droites soudées entre elles.
L’ensemble est lié à l’extérieur par une rotule en 1 et 3, un appui simple en 2.
La structure est en acier de module de Young E.
Les poutres ont une section droite rectangulaire de dimensions (b, t).
Le nœud 2 porte une force de composante (F, 0).
On donne :
E = 200000 MPa
L=1m
t = 20 mm , b = 100 mm
F = −10 kN
Calculer le coefficient de charge critique en utilisant plusieurs maillages.

Résultats :
La charge critique est égale à :

13.9 EIZ
FC = = 18.53 kN
L2
Le coefficient de charge critique est donc égal à :
λC = 18.53
On obtient avec RDM – Ossatures :
Nombre d’éléments λC
3 24.82
5 18.70
10 18.53
20 18.52
F2 : Poutre droite
Référence : Z.P. Bažant,L. Cedolin, Stability of structures, Oxford, 1991, page 70.
Problème :
L’ossature plane représentée sur la figure est constituée de deux poutres droites de longueur L et de
section rectangulaire. Elle est liée à l’extérieur par une rotule en 1 et un appui simple en 2. Soit E le
module de Young du matériau. La poutre porte en 3 une force (P, 0).
On donne :
L = 0.8 m , b = 25 mm , t = 10 mm
E = 210000 MPa , P = −1000 N
Résultats :
b t3
Posons : IZ =
12
π 2 EIZ
PC = 0.1813 = 1223 N
L2
On en déduit λC = 1.223 .
2 1.227
4 1.223
10 1.223
F3 : Poutre droite à section variable
Référence : S.P. Timoshenko, J.M. Gere, Théorie de la stabilité élastique, Dunod, 1966, page 127.
Problème :
La poutre droite 1 − 2 de longueur L est encastrée en 1. Soit E le module de Young du matériau. La

section droite est un rond plein dont le diamètre varie linéairement entre les nœuds 1 et 2. La poutre
porte en 2 une force (P, 0).
On donne :
L = 1.2 m
d1 = 50 mm , d2 = 28.117 mm (Iz2 = 0.1 Iz1 )
E = 200000 MPa
P = −10000 N
Résultats :

EIz1
PC = 1.202 = 51218 N
L2
On en déduit : λC = 5.1218 .
1 5.322
2 5.154
10 5.127
F4 : Poutre console – flexion-torsion

Référence : Solution analytique
Problème :
La poutre droite (ossature spatiale), représentée sur la figure, a une longueur L et une section
constante (rectangle plein : b × t). Elle est en acier de constantes élastiques E et ν. Elle est encastrée
en 1.
Cas de charge 1 : le nœud 2 porte une force (0, 0, −P ).
Cas de charge 2 : la poutre porte une charge uniformément répartie sur toute sa longueur (0, 0, −q).
Cas de charge 3 : le nœud 2 porte un couple (M, 0, 0).
On donne :
L = 1.2 m , b = 100 mm , t = 6 mm
E = 200000 MPa , ν = 0.3
P = 100 N , q = 100 N/m , M = 100 Nm
Pour chaque cas de charge, calculer le coefficient de charge critique en utilisant plusieurs maillages
et plusieurs hypothèses de calcul (petites rotations ou rotations modérées).
Modélisation :
La poutre est une ossature spatiale.
Pour évaluer la constante de torsion de Saint Venant, activer le menu Calculer section droite (∼ 600
triangles à 6 nœuds). Les caractéristiques de la section sont : IY = 0.18 cm4 , J = 0.6928 cm4 .
Résultats :
Cas 1 :

4.0126 p
FC = EIY GJ = 1221 N (petites rotations et rotations modérées)
L2
Le coefficient de charge critique est donc égal à : λC1 = 12.21
Cas 2 :

12.85 p
qC = EIY GJ = 3257 N/m (petites rotations et rotations modérées)
L3
Cas 3 :
hypothèse petites rotations :

π √
La charge critique est égale à : MC = EIY GJ = 573.4 Nm
2L
hypothèse rotations modérées :

π√
La charge critique est égale à : MC = EIY GJ = 1146.7 Nm
L
Remarque : quand les rotations ne sont pas petites, le résultat dépend de la manière dont le
couple extérieur est appliqué. Le résultat ci-dessus est obtenu avec un couple semi-tangentiel [5].
nombre d’éléments λC1 λC2 λC3 λC4

1 18.25 57.48 6.322 12.644
2 13.05 38.61 5.882 12.644
3 12.54 34.74 5.799 11.995
10 12.24 32.75 5.739 11.514
20 12.21 32.62 5.735 11.479
solution analytique 12.21 32.57 5.734 11.467
F5 : Lame équerre – flexion-torsion

Référence : J.-H. Argyris, O. Hilpert, G.-A. Malejannakis, D.-W. Scharpf, On the geo-
metrical stiffness of a beam in space – a consistent v. w. approach, CMAME, vol 20, (1979), 105-131.
Problème :
La structure spatiale représentée sur la figure est composée de deux poutres droites de longueur L
et de section constante (rectangle plein : b × t ). Elle est encastrée en 1.
Soient E et ν les constantes élastiques du matériau.
Le nœud 3 porte une force (0, P, 0) où P peut être positif ou négatif.
On donne :
L = 240 mm , b = 30 mm , t = 0.6 mm
E = 71240 MPa , ν = 0.31
P = ±1 N
Calculer le coefficient de charge critique en utilisant plusieurs maillages et plusieurs hypothèses de

calcul (petites rotations ou rotations modérées).
Modélisation :
Modéliser la section droite comme une section quelconque :
Forme = 5
b t3 t b3
A = b t = 0.18 cm2 , IY = = 0.000054 cm4 , IZ = = 0.135 cm4
12 12
b t3
J= = 0.000216 cm4
3
Résultats :
Référence (avec 2 × 10 éléments) :
hypothèse petites rotations : λC (P > 0) = 0.5507 , λC (P < 0) = 0.4214
hypothèse rotations modérées : λC (P > 0) = 1.0880 , λC (P < 0) = 0.6804
petites rotations rotations modérées

Nombre d’éléments λC (P > 0) λC (P < 0) λC (P > 0) λC (P < 0)
2×1 0.5604 0.4269 1.1754 0.7085
2×2 0.5531 0.4227 1.1101 0.6873
2 × 10 0.5507 0.4214 1.0880 0.6804
2 × 20 0.5506 0.4213 1.0873 0.6802
F6 : Lame équerre – flexion-torsion

Référence : J.-H. Argyris, O. Hilpert, G.-A. Malejannakis, D.-W. Scharpf, On the geo-
metrical stiffness of a beam in space – a consistent v. w. approach, CMAME, vol 20, (1979), 105-131.
Problème :
La structure spatiale représentée sur la figure est composée de deux poutres droites de longueur L,
perpendiculaires entre elles et de section constante (rectangle plein : b × t).
Soient E et ν les constantes élastiques du matériau.
Les conditions aux limites sont :
nœud 1 : u = v = w = θy = θz = 0
nœud 3 : u = w = θy = θz = 0
Cas de charge 1 :
nœud 1 : un couple de composante (−M, 0, 0)

nœud 3 : un couple de composantes (M, 0, 0)
où M peut être positif ou négatif.
Cas de charge 2 :
nœud 2 : une force de composantes (0, 0, P )
où P peut être positif ou négatif.

On donne :
L = 240 mm , b = 30 mm , t = 0.6 mm
E = 71240 MPa , ν = 0.31
M = ±1 N.mm , P = ±1 N
Calculer le coefficient de charge critique en utilisant plusieurs maillages et plusieurs hypothèses de

calcul (petites rotations ou rotations modérées).
Modélisation :
Modéliser la section droite comme une section quelconque :
Forme = 5
b t3 t b3
A = b t = 0.18 cm2 , IY = = 0.000054 cm4 , IZ = = 0.135 cm4
12 12
b t3
J= = 0.000216 cm4
3
Résultats :
Cas de charge 1 :
hypothèse petites rotations : λC (M > 0) = 315.79 , λC (M < 0) = 937.84
hypothèse rotations modérées : λC ( M > 0 ) = 624.77 , λC ( M < 0 ) = 624.77
On obtient (4 modes demandés, précision sur le calcul des valeurs propres = 0.0001) :
nombre petites rotations rotations modérées

d’éléments λC (M > 0) λC (M < 0) λC (M > 0) λC (M < 0)
2×4 317.31 985.38 638.30 638.30
2 × 10 315.79 937.84 624.77 624.77
2 × 20 315.58 931.14 622.85 622.85
2 × 50 315.51 929.27 622.31 622.31
Remarque : la charge critique théorique (hypothèse rotations modérées) est égale à :

πp
MC = EIY GJ = 622.21 N.mm
L
pour M positif ou négatif. Cette valeur est indépendante de l’angle que font entre elles les deux
poutres.
Cas de charge 2 :
hypothèse petites rotations : λC (P > 0) = 19.326 , λC (P < 0) = 2.419
hypothèse rotations modérées : λC (P > 0) = 11.744 , λC (P < 0) = 3.947
On obtient (5 modes demandés, précision sur le calcul des valeurs propres = 0.0001) :
nombre petites rotations rotations modérées

d’éléments λC (P > 0) λC (P < 0) λC (P > 0) λC (P < 0)
2×4 15.419 2.420 12.265 3.951
2 × 10 14.908 2.419 11.744 3.947
2 × 20 14.836 2.419 11.672 3.946
Remarque : la valeur λC (P > 0, hypothèse petites rotations) donnée dans la référence corres-
pond au premier mode symétrique. On obtient (5e valeur propre) : 19.326 avec 20 éléments.
F7 : Flambement d’un mât vertical sous son poids propre
Référence : J. Courbon, Stabilité de l’équilibre élastique, Les Techniques de l’Ingénieur, C2040.
Problème :
Le mât représenté sur la figure est encastré à sa base et libre à son extré-
mité supérieure. Ce mât de hauteur H, de section droite constante : rond
plein de diamètre D est soumis à son poids propre. Soient E le module
de Young du matériau et ρ sa masse volumique. Soit g l’accélération de la
pesanteur.
On donne :
H = 4 m , D = 30 mm
E = 200000 MPa , ρ = 7800 kg/m3
g = 10 m/s2
Évaluer le coefficient de charge critique en utilisant plusieurs maillages.
Résultats :
La charge critique par unité de longueur est égale à :

E IZ
pC = 7.8373 = 973.804 N/m
H3
Le poids propre par unité de longueur étant égal à :
π D2
p= ρ g = 55.035 N/m ,
4
on en déduit :
pC
λC = = 17.662
p
1 17.779
2 17.707
3 17.673
4 17.666
10 17.662
F8 : Flambement d’une poutre droite
Problème :
La poutre droite représentée ci-dessous, de longueur L = 1.2 m et de section droite constante (rectangle
plein : 20 x 100 mm) est en acier de module Young E = 200000 MPa. Elle porte à son extrémité
supérieure une force de composantes (0, P = −1000) N.
Calculer le coefficient de charge critique pour les conditions aux limites suivantes :
Cas 1 2 3 4
Base encastrement rotule encastrement encastrement
Extrémité supérieure libre u=0 u=0 u = 0 , θz = 0
Résultats :
π 2 EIz
Référence : λC1 = 0.25 λ , λC2 = λ , λC3 = 2.04575 λ , λC4 = 4 λ avec λ =
|P | L2
On obtient :
nombre d’éléments λC1 λC2 λC3 λC4

1 23.018 111.110
2 22.858 92.073 191.750 370.37
3 22.849 91.530 188.100 373.550
4 22.847 91.432 187.340 368.300
5 22.847 91.405 187.110 366.720
20 22.846 91.385 186.950 365.550
solution analytique 22.846 91.385 186.951 365.541
F9 : Flambement d’un cadre

Référence : C. Massonnet, Résistance des matériaux, Dunod, 1968, page 410.
Problème :
Le cadre représenté sur la figure est constitué de quatre

poutres de longueur L et de section droite constante :
rectangle plein (cY × cZ ). Soit E le module de Young du
matériau. Le cadre est articulé en 1 et 4. Il porte en 2
et 3 deux forces égales de composantes (0, −P ).
On donne :
L = 0.6 m
cY = 10 mm , cZ = 50 mm
E = 200000 MPa
P = 1000 N
Calculer le coefficient de charge critique λC quand le déplacement horizontal du point 2 est libre et
quand celui-ci est nul.
Résultats :
– Première étude : le nœud 2 est libre

EIZ
PC = 5.68783 = 13166 N
L2
On en déduit :
PC
λC = = 13.166
P
4 13.194
8 13.181
12 13.168
16 13.165
référence 13.166
– Deuxième étude : le déplacement horizontal du nœud 2 est nul

EIZ
PC = 16.4634 = 38110 N
L2
On en déduit :
PC
λC = = 38.110
P
nombre d’éléments λC
4
8 38.468
12 38.209
16 38.144
40 38.111
référence 38.110
Chapitre 5
Modes propres
D1 : Treillis plan à nœuds articulés

Référence : M. Géradin, D. Rixen, Théorie des vibrations, Masson, 1996, page 265.
Problème : l’ossature plane représentée sur la figure est constituée de neuf poutres droites articulées
entre elles. Elle est liée à l’extérieur par une rotule en 1 et un appui simple en 2. Les poutres sont des
carrés creux de côté extérieur c et d’épaisseur t. Soient E le module de Young du matériau et ρ sa
masse volumique.
On donne :
L = 1 m , c = 40 mm , t = 5 mm , E = 200000 MPa , ρ = 8000 kg/m3
Calculer les 9 premières fréquences propres de membrane.
Modélisation : pour obtenir les vibrations de membrane, ne pas discrétiser les poutres.
Résultats : on obtient (fréquences en Hz) :
Mode Référence RDM – Ossatures

1 171.40 171.39
2 290.50 290.48
9 1663.5 1663.41
D2 : Poutre droite à section variable

Problème :
La poutre droite 1−2 de longueur L est encastrée en 1. Soient E le module de Young du matériau et ρ
sa masse volumique. La section droite est un rectangle plein dont les dimensions varient linéairement
entre les nœuds 1 et 2.
On donne :
L = 1 m , E = 200000 MPa , ρ = 7800 kg m−3
hY 1 = 40 mm , hZ1 = 50 mm
hY 2 = 10 mm , hZ2 = 10 mm
Calculer les 5 premières fréquences propres.
Modélisation :
Modéliser la poutre comme une ossature plane. Utiliser plusieurs maillages.
Résultats :
On obtient (fréquences en Hz) pour les modes de flexion :
Mode Référence 1 élément 2 éléments 5 éléments 10 éléments

1 56.55 56.81 56.59 56.55 56.55
2 175.79 180.81 175.83 175.73
3 389.01 404.07 390.21 388.65
4 702.36 903.32 714.02 701.25
5 1117.63 1186.12 1115.82
D3 : Vibrations transversales d’une poutre droite bi-encastrée

Référence : R.D. Blevins, Formula for natural frequency and mode shape, Krieger, 1993, page 108.
Problème :
L’ossature plane représentée sur la figure est constituée d’une poutre droite 1 − 2 de longueur L et de
section constante : carré plein de côté c. Elle est encastrée en 1 et 2. Soient E le module de Young du
matériau et ρ sa masse volumique.
On donne :
L = 1 m , E = 210000 MPa , ρ = 7800 kg m−3
c = 10 mm
Calculer les 5 premières fréquences propres en utilisant plusieurs maillages.
Résultats :
Référence :
s
1 h2i EIZ
fi = avec hi = 4.73004, 7.85320, 10.9956, 14.1372, 17.2788
2 π L2 ρA
Les fréquences en Hz obtenues avec RDM – Ossatures sont :
Mode Référence 2 éléments 3 éléments 10 éléments 20 éléments

1 53.34 54.20 53.55 53.34 53.33
2 147.02 195.38 149.93 147.03 147.00
3 288.22 348.62 288.39 288.12
4 476.45 692.49 477.37 476.19
5 711.73 715.02 711.22
D4 : Portique plan
Problème : l’ossature plane représentée sur la figure est constituée de 6 poutres droites de section
constante : rectangle plein (b, h). Elle est encastrée en 1 et 2. Soient E le module de Young du matériau
et ρ sa masse volumique.
On donne :
b = 29 mm , h = 4.8 mm
E = 210000 MPa , ρ = 7800 kg/m3
Mode Référence 6 éléments 20 éléments 60 éléments

1 8.8 8.79 8.78 8.78
2 29.4 29.52 29.44 29.44
3 43.8 52.93 43.87 43.85
4 56.3 86.77 56.35 56.30
5 96.2 118.64 96.41 96.18
13 335 343.36 335.48
D5 : Ossature spatiale
Référence : M. Petyt, Introduction to finite element vibration analysis, Cambridge University Press,
1990, page 108.
Problème : l’ossature spatiale représentée sur la figure est constituée de 16 poutres droites. Elle est
encastrée à sa base. Soient E et ν les caractéristiques élastiques du matériau et ρ sa masse volumique.
On donne :
E = 219900 MPa , ν = 0.3 , ρ = 7900 kg/m3
L = 1 m , c = 50 mm , b = 150 mm , h = 50 mm

1 11.8 11.81 11.81 11.81 11.81
3 15.38 15.38 15.38 15.38
4 34.1 34.13 34.11 34.11 34.11
6 43.28 43.25 43.24 43.24
7 134.76 122.05 121.59 121.56
10 178.04 153.70 152.81 152.75
D6 : Ossature plancher
Référence : J.P Rezette, F. Leleux, Calcul dynamique des structures par la méthode des éléments
finis, Les notes techniques du CETIM, 1974, page 58.
Problème : l’ossature plancher représentée sur la figure est constituée de 40 poutres droites ( ronds
pleins de diamètre d). Soient E et ν les caractéristiques élastiques du matériau et ρ sa masse volumique.
Les nœuds extérieurs reposent sur un appui simple.
On donne :
E = 200000 MPa , ν = 0.3 , ρ = 8000 kg/m3
d = 0.01 m
Calculer les 6 premières fréquences propres de flexion-torsion en utilisant plusieurs maillages.
Modélisation :
Ossature plancher paramétrée : numéro 50 : L = 0.8 m , H = 0.4 m , N = M = 4

1 96 96 96 96 96
2 165 165 165 165 165
3 278 278 276 275 275
4 306 306 301 300 300
5 369 370 361 361 361
6 468 469 453 452 452
Remarque : dans la référence, les calculs sont effectués avec 40 éléments.

D7 : Vibrations transversales d’une poutre droite libre

Problème : la poutre droite de longueur L représentée a une section constante : carré plein de côté c.
Soient E le module de Young du matériau et ρ sa masse volumique.
On donne :
E = 210000 MPa , ρ = 7800 kg/m3
L = 1.2 m , c = 20 mm
Étudier les 5 premiers modes propres élastiques de flexion en utilisant plusieurs maillages.
Calcul : introduire un décalage spectral égal à 20 Hz (il y a 3 modes rigides).
Résultats :
Référence :
s
h2i EIZ
fi =
2 π L2 ρA
avec
hi = 4.73004, 7.85320, 10.9956, 14.1372, 17.2788
On obtient (fréquences en Hz) :
Mode Référence 10 éléments 20 éléments 40 éléments

1 74.08 74.04 74.04 74.04
2 204.20 203.99 203.95 203.94
3 400.31 399.81 399.47 399.45
4 661.73 661.14 659.67 659.57
5 988.52 988.91 984.31 983.97
D8 : Premier mode propre d’une poutre console avec masses

Problème :
La poutre console de longueur L représentée sur la figure est un rectangle plein de base b et de
hauteur h. Soient E et ν les caractéristiques élastiques du matériau et ρ sa masse volumique. La
poutre porte une masse ponctuelle M à son extrémité et une masse uniformément répartie sur toute
sa longueur d’intensité m.
On donne :
L = 0.8 m , E = 200000 MPa , ν = 0.3 , ρ = 7800 kg m−3
b = 100 mm , h = 10 mm
M = 2 kg , m = 4 kg/m
Problème : étudier le premier mode propre en utilisant plusieurs maillages.
Résultats :
Référence :
s
1 3 E Iz bh3
f= 3
avec A = b h , Iz =
2π L (M + 0.24267 ( ρ A + m) L) 12
On obtient ( fréquences en Hz ) :
Maillage M =m=0 M 6= 0 , m = 0 M = 0 , m 6= 0 M 6= 0 , m 6= 0
1 élément 12.84 8.43 8.95 6.98
2 éléments 12.79 8.42 10.00 7.45
3 éléments 12.78 8.42 10.21 7.55
20 éléments 12.78 8.42 10.39 7.62
référence 12.78 8.39 10.39 7.59
Bibliographie
[1] J.-L. Batoz et G. Dhatt – Modélisation des structures par éléments finis, Volume 2. Poutres
et plaques, Hermès, 1990.
[2] W. D. Pilkey – Formulas for stress, strain and structural matrices, Wiley, 1994.
[3] S. P. Timoshenko – Résistance des matériaux, Tome 2. Théorie développée et problèmes, Dunod,
1968.
[4] W. C. Young et R. G. Budynas – Roarks formulas for stress and strain, McGraw-Hill, 2002.
[5] H. Ziegler – Principles of structural stability, 2 éd., Birkauser Verlag, 1977.

Exos RDM

Transféré par

Droits d'auteur :

Formats disponibles

Exos RDM

Transféré par

Informations du document

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Exos RDM

Transféré par

Droits d'auteur :

Formats disponibles

RDM – Ossatures

26 juin 2006 – 29 mars 2011

3 Sections droites : caractéristiques et contraintes 32

Exemple 1 : Portique plan

Le nœud 2 porte une force de composantes (P, 0, 0).

L’énergie de déformation due à l’effort tranchant est négligée (modèle de Bernoulli).

Les étapes de la modélisation sont :

u2 = 2.2144 mm , v2 = −0.0017 mm , θ2z = −0.0388˚

R1x = −6077.4 N , R1y = 533.4 N , M1z = 3221.6 N.m

Exemple 2 : Treillis plan à nœuds articulés

Les poutres 1 − 2 et 1 − 4 sont en acier :

module de Young = 200000 MPa

La poutre 1 − 3 est en laiton :

module de Young = 100000 MPa

Le nœud 1 porte une charge P~ de composantes (0, −10000, 0) N.

L’ossature subit une augmentation de température de 50 K.

Les étapes de la modélisation sont :

– Allongement des poutres :

∆1−2 = ∆1−4 = 0.768 mm , ∆1−3 = 0.960 mm

Remarque : pour extraire ces résultats, utiliser le bouton droit de la souris.

Exemple 3 : Anneau plan

Référence : solution analytique.

E = 200000 MPa , ν = 0.3

Les étapes de la modélisation sont :

Solution éléments finis :

F1x = 0 N , M1z = −18977 N.mm , F3y = 500 N , M3z = −18523 N.mm

– Moment fléchissant dans la section 2 : Mfz2 = 6477 N.mm

σa = 113.86 MPa , σb = −113.86 MPa , σc = 106.14 MPa , σd = −116.14 MPa

σc = 99.10 MPa , σd = −124.00 MPa

Les sections 1 , 3 , 5 et 6 sont encastrées.

L’énergie de déformation due à l’effort tranchant est négligée (modèle de Bernoulli).

Les étapes de la modélisation sont :

Définir le type de l’ossature : Plancher

w2 = −1.2182 mm , θ2x = −0.35599 10−3 rad , θ2y = −0.14976 10−3 rad

w4 = −2.0993 mm , θ4x = 0.28856 10−3 rad , θ4y = 0.18376 10−3 rad

F1z = 93.528 kN , M1x = 9.493 kN.m , M1y = −163.092 kN.m

F3z = 34.452 kN , M3x = 14.240 kN.m , M3y = 76.393 kN.m

Exemple 5 : Ossature spatiale

Les coordonnées nodales sont :

nœud x (m) y (m) z (m)

Les caractéristiques élastiques du matériau sont : E = 100000 MPa et ν = 0.2987.

Les sections 1 et 6 sont encastrées.

Le nœud 4 porte une force F~ de composantes (0, 0, −1000) daN .

Les étapes de la modélisation sont :

Moments aux extrémités des poutres (en daN.m)

élément Mto MfY o MfZo Mte MfY e MfZe

Exemple 6 : Modes propres d’un anneau plan

L’anneau et la patte ont des sections droites rectangulaires pleines.

Soient E le module de Young du matériau et ρ sa masse volumique.

On recherche les six premiers modes propres de cet anneau.

L’énergie de déformation due à l’effort tranchant est négligée (modèle de Bernoulli).

E = 72000 MPa , ρ = 2700 kg/m3

Section droite de l’anneau : Ha = 5 mm , Ba = 10 mm

Section droite de la patte : Hp = 3 mm , Bp = 10 mm

Les étapes de la modélisation sont :

Bibliothèque (une partie de la géométrie existe dans la bibliothèque d’ossatures paramétrées)

Ajouter une poutre verticale

Mode Référence RDM – Ossatures

Exemple 7 : Ossature plane