Cours Complet RDM 2

RDM 5.
xx
de Yves DEBARD (IUT Le Mans)
Module
Elasticité Plane
RDM 5.xx Page 1
www.GenieCivilPDF.com
Sommaire
1. Rappel d'élasticité.......................................................................................................................... 5
1.1 Contraintes...................................................................................................................... 5
1.1.1 Etat des contraintes en un point ...................................................................... 5
→
1.1.2 Vecteur contrainte en un point pour une direction n .................................... 5
1.1.3 Contrainte normale et tangentielle.................................................................. 6
1.1.4 Contraintes principales et directions principales............................................ 6
1.1.5 Etats de contraintes particuliers...................................................................... 8
1.1.6 Cercles de Mohr.............................................................................................. 10
1.1.6.1 Cercles de Mohr des contraintes..................................................... 10
1.1.6.2 Cercles de Mohr en contraintes planes ........................................... 11
1.1.6.3 Construction de Mohr ..................................................................... 13
1.1.7 Etat des contraintes autour d'un point sur la surface d'un solide .................... 13
1.2 Déplacements - Déformations ........................................................................................ 15
1.2.1 Champ des déplacements................................................................................ 15
1.2.2 Etat des déformations au voisinage d'un point................................................
15
1.2.3 Allongement unitaire en A et pour une direction q ....................................... 15
1.2.4 Glissement. Distorsion.................................................................................... 16
1.2.5 Déformations principales et directions principales ........................................ 17
1.2.6 Etat de déformations planes............................................................................ 17
1.3 Loi de comportement ...................................................................................................... 18
1.4 Critères de limite élastique ............................................................................................. 19
1.4.1 Critère de Tresca ou du cisaillement maximal ............................................... 19
1.4.2 Critère du plus grand travail de distorsion. Critère de Von Mises ................. 23
1.4.3 Critères exprimés dans le cas de contraintes planes ....................................... 23
1.4.3.1 Critère de Tresca ............................................................................. 23
1.4.3.2 Critère de Von Mises ...................................................................... 24
1.5 Types particuliers de problème d’élasticité.................................................................... 24
1.5.1 Contraintes planes du plan (x1,x2) ................................................................ 25
1.5.2 Déformations planes ....................................................................................... 25
1.5.3 Problèmes axisymétriques méridiens.............................................................. 26
2. Fondement mécanique de la méthode des éléments finis.............................................................. 27
2.1 Notations......................................................................................................................... 27
2.2 Problème d'élasticité. Equations d'équilibre ................................................................... 28
2.3 Loi de comportement ...................................................................................................... 29
2.4 Energie de déformation élastique ................................................................................... 29
2.5 Théorème d'unicité.......................................................................................................... 30
2.6 Champ de déplacement virtuel admissible ..................................................................... 30
2.6.1 Définition ........................................................................................................ 30
2.6.2 Conséquence ................................................................................................... 30
2.7 Energie potentielle d'un système élastique ..................................................................... 30
2.7.1 Système à un degré de liberté ......................................................................... 31
2.7.2 Système à plusieurs degrés de liberté ............................................................. 32
2.7.3 Système continu .............................................................................................. 32
2.8 Approximation par éléments finis .................................................................................. 33
2.8.1 Définitions ...................................................................................................... 34
2.9 Méthode des éléments finis en élasticité, conduite à partir des déplacements............... 34
2.10 Application à l'élément triangulaire à trois noeuds ...................................................... 36
2.10.1 Construction de la loi d'interpolation............................................................ 36
2.10.2 Tenseur des déformations approchées .......................................................... 39
2.10.3 Loi de comportement .................................................................................... 39
2.10.4 Matrice de rigidité de l'élément .................................................................... 40
RDM 5.xx Page 2
2.10.5 Vecteur d'effort ............................................................................................. 40
2.11 Application à une poutre en flexion simple.................................................................. 41
2.11.1 Décomposition en éléments finis.................................................................. 41
2.11.2 Calcul de la matrice rigidité de chaque élément ........................................... 42
2.11.3 Assemblage de la matrice ............................................................................. 44
2.11.4 Introduction des conditions aux limites........................................................ 47
2.11.5 Résolution du système linéaire ..................................................................... 47
2.11.6 Calcul des déformations et des contraintes................................................... 48
2.11.6.1 Vecteur déplacement en un point quelconque .............................. 48
2.11.6.2 Tenseur des déformation en un point quelconque ........................ 48
2.11.6.3 Tenseur des contraintes en un point quelconque .......................... 48
2.11.7 Organigramme du processus de résolution ................................................... 49
2.12 Types d'éléments plus performants............................................................................... 50
3. Type d'éléments finis et indicateurs de choix................................................................................ 51
3.1 Problèmes pratiques posés à l'utilisateur de la méthode des éléments finis ................... 51
3.1.1 Choix des éléments ......................................................................................... 51
3.1.2 Influence du maillage. Etude sur un cas test................................................... 52
4. Description des possibilités du logiciel ......................................................................................... 56
4.1 Modélisation ................................................................................................................... 58
4.1.1 Modélisation de la géométrie.......................................................................... 59
4.1.2 Modélisation du maillage............................................................................... 64
4.1.2.1 Définition des paramètres de maillage............................................ 64
4.1.2.2 Choix de l’élément .......................................................................... 66
4.1.2.3 Vérification de la qualité du maillage............................................ 71
4.1.2.4 Sauvegarde du maillage ................................................................. 74
4.1.2.5 Maillage par blocs........................................................................... 76
4.2 Sauvegardes des études................................................................................................... 79
4.2.1 Modélisation mécanique ................................................................................. 80
4.2.1.1 Définition des épaisseurs ................................................................ 80
4.2.1.2 Définition des matériaux................................................................. 82
4.2.1.3 Définition des liaisons .................................................................... 84
4.2.1.4 Définition des cas de charges.......................................................... 85
4.2.1.5 Définition du modèle d’étude dynamique....................................... 88
4.3 Calculs ............................................................................................................................ 90
4.3.1 Calcul statique ................................................................................................ 90
4.3.2 Calcul dynamique ........................................................................................... 91
4.4 Exploitation des résultats................................................................................................ 93
4.4.1 Résultats du calcul théorique .......................................................................... 93
4.4.2 Résultats donnés par le logiciel ...................................................................... 94
5. Exemples de modélisation ............................................................................................................. 115
5.1 Couronne de pont 17x56................................................................................................. 115
5.1.1 Problème posé................................................................................................. 115
5.1.2 Données techniques ........................................................................................ 115
5.1.3 Traitement d’un modèle OSSATURE ............................................................ 115
5.1.3.1 Modélisation ................................................................................... 115
5.1.4 Traitement d’un modèle en élasticité plane (contraintes planes) avec le
module M. E. F. ....................................................................................................... 116
5.1.4.1 Modélisation ................................................................................... 116
5.1.4.2 Résultats de l’étude en élasticité plane ........................................... 116
5.1.5 Conclusion ...................................................................................................... 120
5.2 Capteur d'effort "PRECIA-precia".................................................................................. 122
5.2.1 Présentation..................................................................................................... 122
5.2.2 Réalisation ...................................................................................................... 123
5.2.2.1 Forme du corps d'épreuve ............................................................... 123
RDM 5.xx Page 3
5.2.2.2 Jauges de contraintes....................................................................... 123
5.2.3 Corps d'épreuve étudié.................................................................................... 123
5.2.3.1 Epaisseur ......................................................................................... 124
5.2.3.2 Matériau .......................................................................................... 124
5.2.4 Maillage .......................................................................................................... 124
5.2.5 Montage du corps d'épreuve ........................................................................... 124
5.2.6 Liaisons........................................................................................................... 125
5.2.7 Chargement ..................................................................................................... 125
5.2.8 Résultats.......................................................................................................... 125
5.3 Pièce support d'ATR 42 .................................................................................................. 128
5.3.1 L'avion ATR 42.............................................................................................. 128
5.3.2 Conditionnement d'air..................................................................................... 128
5.3.3 Données .......................................................................................................... 129
5.3.3.1 Caractéristiques techniques............................................................ 130
5.3.3.2 Appuis : ........................................................................................... 130
5.3.3.3 Charges extérieures......................................................................... 131
5.3.4 Travail demandé ........................................................................................... 131
5.3.4.1 En théorie des poutres..................................................................... 131
5.3.4.2 En élasticité plane ........................................................................... 131
5.3.5 Résultats.......................................................................................................... 132
5.3.5.1 Appuis de type pivots en C et D ( hyperstatique) ........................... 132
5.3.5.1.1 Maillage........................................................................... 132
5.3.5.1.2 Appuis et déplacements en A et B .................................. 133
5.3.5.1.3 Contraintes équivalente de Von Mises............................ 134
5.3.5.2 Avec appui ponctuel en D (isostatique) .......................................... 135
5.3.5.2.1 Déformée......................................................................... 135
5.3.5.2.2 Contraintes équivalente de Von Mises............................ 135
5.4 Support de galet freineur ................................................................................................ 136
5.4.1 Objectif de ce problème.................................................................................. 136
5.4.2 Présentation du système mécanique ............................................................... 136
5.4.2.1 Description du galet freineur .......................................................... 136
5.4.3 Etude statique.................................................................................................. 139
5.4.3.1 Etude statique du galet .................................................................... 139
5.4.3.2 Etude statique du support (15) ........................................................ 140
5.4.4 Etude en élasticité plane ................................................................................. 144
5.4.4.1 Description des liaisons .................................................................. 145
5.4.4.2 Description des charges .................................................................. 146
5.4.4.3 Conditions de l'étude....................................................................... 147
5.4.5 Modèle définitif .............................................................................................. 148
6. Annexe........................................................................................................................................... 154
6.1 Limites actuelles du logiciel M.E.F................................................................................ 155
6.1.1 Limites de calcul............................................................................................. 155
6.1.2 Sauvegarde des études .................................................................................... 155
6.2 Modélisation d’une articulation...................................................................................... 156
6.3 Courbes CETIM.............................................................................................................. 158
7. Bibliographie ................................................................................................................................. 161
RDM 5.xx Page 4
1. Rappel d'élasticité
1.1 Contraintes
1.1.1 Etat des contraintes en un point
On démontre dans le cours d'élasticité que, compte tenu des hypothèses physiques, l'état de
contrainte en un point A (figure 1. 1) est caractérisé par le tenseur des contraintes. C'est un
tenseur du second ordre symétrique. Dans une base orthonormée il est représenté par la
matrice des contraintes qui s'écrit :
X
3
σ33
σ32
σ23
σ31
σ13 σ22
≈
 σ 11 σ 12 σ 13 
A
σ21
≈
X
2
σ (A) =  σ 22 σ 23 
σ12
σ11
Sym σ 33 
X
1
Figure 1. 1 : Contraintes autour du point A
Dans cette notation du tenseur des contraintes, le premier indice indique la direction de la
normale à la facette, le deuxième : la direction de la contrainte
→
1.1.2 Vecteur contrainte en un point pour une direction n
→
Soit un point A d'un solide, et une direction repérée par un vecteur n (normale extérieure à la
matière) .
→
Soit une facette infiniment petite d'aire dS de normale n . Le vecteur contrainte au point A pour
→
la direction n s'écrit :
→ → ≈ →
T (A, n ) = σ (A) n
Soit :
RDM 5.xx Page 5
 T1  σ 11 σ 12 σ 13   n1  σ 11n1+σ 12n2+σ 13n3 
 T2  =  σ     
   12 σ 22 σ 23   n2  = σ 12n1+σ 22n2+σ 23n3 
 T3   σ 13 σ 23 σ 33   n3  σ 13n1+σ 23n2+σ 33n3 
1.1.3 Contrainte normale et tangentielle

Le vecteur contrainte en un point A et pour la direction n (figure 1. 2) peut être projeté :
- sur la normale, on obtient la contrainte normale :
→ ≈
σ = T (A ,n) . n = t n σ (A) n
σ 11n1+σ 12n2+σ 13n3   n1(σ 11n1+σ 12n2+σ 13n3) 

σ = [ n 1 n2 n3 ] σ n +σ n +σ n  =  n (σ n +σ n +σ n ) 
 12 1 22 2 23 3   2 12 1 22 2 23 3 
σ 13n1+σ 23n2+σ 33n3   n 3(σ 13n1+σ 23n2+σ 33n3) 
- sur le plan tangent, on obtient la contrainte tangentielle τ telle que:
τ ² = T(A , n) ² - σ ²
σ
T(A,n)
A
τ
Figure 1. 2 : Contrainte normale σ et tangentielle τ en un point A
1.1.4 Contraintes principales et directions principales

Mathématiquement on démontre : le tenseur des contraintes étant réel symétrique, il est
→
diagonalisable, c'est à dire qu'il existe un réel σi et une direction Xi telle que:
→ →
≈ (A) X = σ X
σ i i i
- les trois valeurs propres σi sont réelles (distinctes ou confondues) ;

→
- si les trois valeurs propres sont distinctes, les vecteurs propres correspondants X i
sont orthogonaux.
RDM 5.xx Page 6
X3
σ3
σ2
A
X2
σ1
X1
Figure 1. 3: Contraintes principales autour du point A
Traduction mécanique : si les trois contraintes principales σ1, σ2 et σ3 sont distinctes, il

→ → →
existe trois directions principales orthogonales correspondantes X 1, X 2 et X 3.
→
Ainsi pour une telle direction X i le vecteur contrainte :
→ → → →
T ( A , Xi ) = σ i Xi est colinéaire à la direction X i
En d'autres termes, la contrainte tangentielle pour cette direction est nulle. Dans cette direction
on a donc affaire, soit à une sollicitation de traction (σi > 0) , soit à une sollicitation de
compression (σi < 0).
Dans le repère principal la matrice des contraintes s'écrit alors :
≈ σ 1 

σ (A) =  σ2 

 σ 3  ( X , X ,X )
1 2 3
Détermination des contraintes principales et des directions principales (figure 1. 3)
- Les contraintes principales sont déterminées en écrivant que le déterminant

suivant est nul :
σ11 - σ σ 12 σ 13
σ 12 σ 22 - σ σ 23 = 0
σ 13 σ 23 σ 33 - σ
On aboutit à l'équation caractéristique : -σ3 + I1 σ² - I2 σ + I3 = 0
Dans cette équation I1, I2 et I3 sont les trois invariants du tenseur des contraintes (
quantités indépendantes de la base dans laquelle est exprimé le tenseur). Dans une
base quelconque ils ont pour expression :
I1 = σ11 + σ22 + σ33 = trace de σ(A)

I2 = (σ11 σ22 - σ12²) + ( σ22 σ33 - σ23²) + ( σ11 σ33 - σ13²)
≈
I3 = dét [ σ (A)]
RDM 5.xx Page 7
et, dans la base principale :
I1 = σ1 + σ2 + σ3 = trace de σ(A)
I2 = σ1 σ2 + σ2 σ3 + σ3 σ1
≈
I3 = dét [ σ (A)]
→
La direction principale X correspondant à la contrainte principale σ est déterminée en
écrivant :
σ 11 - σ σ 12 σ 13   X1  0 
 σ σ 22- σ σ 23   X2  = 0 
 12    
 σ 13 σ 23 σ 33- σ   X3  0 
≈
Le déterminant de la matrice σ (A) - σ δij étant nul, on ne dispose plus que de deux équations
indépendantes (par exemple les deux premières). On ne peut donc déterminer qu'une direction
et non un vecteur (l'une des composantes est arbitraire), c'est pourquoi l'on parle de direction
principale.
1.1.5 Etats de contraintes particuliers

a) Etat de contrainte uniaxial
Dans le repère principal, le tenseur des contraintes se réduit à :
σ 1 0 0 
≈
σ (A) =  0 0 0 
 0 0 0  (X ,X ,X )
1 2 3
Traction simple si σ1 > 0, compression simple si σ1 < 0
b) Etat de cisaillement simple

Soit un repère orthonormé ( A ; X1, X 2 , X 3 ), l'état de contraintes en A est un état de cisaillement

simple par rapport aux directions x 1 et x 2 si le tenseur des contraintes se réduit à :
0 τ 0 
≈
σ (A) = τ 0 0 
 0 0 0  ( x , x , x )
1 2 3
Les contraintes principales sont égales à :
σ1 = τ , σ2 = -τ , σ3 = 0
Dans le repère principal, le tenseur des contraintes admet la forme suivante :
RDM 5.xx Page 8
τ 0 0
≈
σ (A) = 0 −τ 0 
0 0 0  (X , X , X )
1 2 3
c) Etat de contraintes planes

On a affaire à un état plan de contraintes parallèlement au plan ( x 1, x 2) si :
σ13 = σ23 = σ33 = 0

σ11 = σ11 (x1,x2) , σ22 = σ22 (x1,x2) et σ12 = σ12 (x1,x2)
C'est le cas des plaques planes chargées dans leur plan.
Le tenseur des contraintes s'écrit alors dans une base quelconque :
σ 11 σ 12 0 
≈
σ (A) = σ 12 σ 22 0 
 0 0 0  ( x , x , x )
1 2 3

L'axe x 3 est donc direction principale et la contrainte principale correspondante est nulle.
Dans la base principale :
σ 1 0 0 
≈
σ (A) =  0 σ 2 0 
 0 0 0  (X , X , X )
1 2 3
d) Etat de contrainte dans une section droite de poutre
En tout point A d'une section droite, l'état de contrainte peut se représenter dans la base locale

classique de la théorie des poutres ( x 1, x 2, x 3 ) :

- x 1 tangent à la ligne moyenne ;

- x 2 et x 3 dans le plan de section droite et axes principaux de la section.
σ 11 σ 12 σ 13 
≈
σ (A) = σ 12 0 0 
σ 13 0 0  (x , x , x )
1 2 3
Traction (compression) simple :
RDM 5.xx Page 9
σ 0 0
≈
σ (A) =  0 0 0
 0 0 0 ( x , x , x )
1 2 3
Flexion pure :
σ 0 0
≈
σ (A) =  0 0 0
 0 0 0 ( x , x , x )
1 2 3
Flexion simple :
 σ σ 12 σ 13 
≈
σ (A) = σ 12 0 0 
σ 13 0 0  (x , x , x )
1 2 3
Torsion avec sections circulaires :
 0 σ 12 σ 13 
≈
σ (A) = σ 12 0 0 
σ 13 0 0  ( x , x , x )
1 2 3
1.1.6 Cercles de Mohr

1.1.6.1 Cercles de Mohr des contraintes
Supposons connues les trois contraintes principales σ1, σ2 et σ3 au point M. On peut montrer
→
(σ,τ) (appelé plan de Mohr), l'extrémité des vecteurs contraintes T (M, n )
que, dans le plan
admissibles, n tournant autour du point M , est la surface ombrée de la figure 1. 4.
RDM 5.xx Page 10
τ
τMax
τ
T(A,n)
σ3 σ2 σ σ1 σ
Figure 1. 4a : Tricercle de Mohr

Cas particulier important : Dans le cas où le vecteur n décrit un plan principal, par exemple
→ → →
( X 1, X 2 ), l'extrémité du vecteur T (M, n ) décrit, dans le plan ( σ, τ ) , un cercle centré sur
l'axe Oσ et ayant pour diamètre le segment σ1, σ2 . C'est notamment le cas en contraintes
planes.
1.1.6.2 Cercles de Mohr en contraintes planes
En contraintes plane σ3 = 0 ainsi la figure ci-dessus devient :
τ
τMax
τ
T(A,n)
σ3 = 0 σ2 σ σ1 σ
Figure 1. 4b : Tricercle de Mohr en contraintes planes
RDM 5.xx Page 11
→ →
Plaçons nous dans le plan principal (A ; X 1, X 2)
→
Considérons une facette dont la normale n fait un angle de θ par rapport à l'axe principal X 1.

Soit un vecteur unitaire t tangent à la facette, obtenu par rotation de - π/2 de n .

n = cosθ X1 + sinθ X2

t = sinθ X1 - cosθ X2
X2
τ
n
M
θ T(M,n)
X1 θ 2 θ σ2
facette σ1 B σ
C A
Figure 1. 5 : Construction de Mohr
Le vecteur contrainte admet pour expression :
→ → σ 1 0   cos θ  → →
T ( A, n ) =     = σ cos θ X + σ sin θ X
 0 σ 2   sin θ 
1 1 2 2
→
la contrainte normale est : σ = T ( A, n ) . n = σ1 cos²θ + σ2 sin²θ
et la contrainte tangentielle définie ici par :
→
τ = T ( A, n ) . t = σ1 cosθ sinθ - σ2 sinθ cosθ
en passant en arc double, il vient :
 σ1+σ 2 σ1-σ 2
 σ = + cos (2θ )
2 2

 τ = σ 1 - σ 2 sin (2θ )
 2
RDM 5.xx Page 12

Lorsque la facette orientée par la normale extérieure n tourne, l'angle θ varie. Il lui
correspond le point figuratif M (figure 1. 5)
Dans le plan de Mohr, le lieu de M lorsque θ varie est un cercle de centre :
σ1 σ2
σ=
+
2 σ1 - σ 2
C et de rayon
2
τ=0
1.1.6.3 Construction de Mohr
Convenons de classer les contraintes principales dans l'ordre suivant :

σ1 > σ2
Dans le plan de Mohr (Figure 1. 5) l'axe horizontal est gradué en contrainte normale σ et l'axe
vertical en contrainte tangentielle τ.
θ=0 θ= π/2
Le point A → σ = σ 1 et le point B → σ = σ2
τ=0 τ=0
Le point M est le point courant. Il correspond à une facette qui a tourné de θ par rapport à l'axe
→ → →
X 1. L'angle au centre ( CA ,CM ) est alors de 2θ.
Remarque : Ne pas oublier que σ3 = 0. Il y a donc trois cercles de Mohr.
1.1.7 Etat des contraintes autour d'un point sur la surface d'un
solide

Considérons un point A sur la surface d'un solide. Soit un repère (A; x 1, x 2, x 3) tel que :

- x 1 et x 2 soient situés dans le plan tangent ;

- x 3 est dirigé suivant la normale
n ,x 3
A
θ x2
x1
t
Figure 1. 6 : Point A sur la surface d'un solide
RDM 5.xx Page 13
Supposons qu'en A il n'y a pas de chargement. C'est donc que le vecteur chargement extérieur
→
est nul : Q (A) = 0
→ →
Les conditions aux limites sur le chargement imposent : T (A, x 3 ) = Q (A)
→
alors : T (A, x 3 ) = 0 (1)

Projetons la relation (1) sur la tangente t définie par θ et sur la normale.
→
T ( A , x3 ) . t = 0 (2)
→
T ( A , x3 ). x 3 = 0 (3)
Signification de la relation (2) :
σ 11 σ 12 σ 13  0 σ 13 
/ /
T ( A, x 3) =  σ 12 σ 22 σ 23  0 = σ 23 
 σ 13 σ 23 σ 33  1  σ 33 
/ / /
t = cos θ x 1 + sin θ x2
→
alors T (A, x3 ) . t = σ13 cos θ + σ 23 sin θ = 0
et ceci quel que soit θ , c'est donc que: σ13 = 0 et σ23 = 0
≈ σ 11 σ 12 0 

Ainsi le tenseur des contraintes s'écrit en A: σ ( A) = σ 12 σ 22 0

 
 0 0 σ 33 
Ainsi, lorsqu'il n'y a pas de chargement tangentiel à la surface au point A, la normale à la

≈
surface est une direction principale car: σ (A) n = σ 33 n

Le plan ( x 1, x 2) est alors plan principal.
Signification de la relation (3) :
→
T (A , x3 ) . x3 = 0 ⇒ σ33 = 0
≈ σ 11 σ 12 0

En définitive, le tenseur des contraintes s'écrit en A: σ ( A) = σ 12 σ 22 0

 
 0 0 0
Conclusion : Lorsque en un point A de la surface d'un solide il n'y a pas de charge extérieure,
l'état de contrainte est un état de contrainte plane dans le plan tangent à la surface.
Cette propriété est mise en oeuvre en extensométrie.
RDM 5.xx Page 14
1.2 Déplacements - Déformations
1.2.1 Champ des déplacements
Sous l'effet des efforts, la structure se déforme. Un point M de coordonnées (x1, x2, x3)
appartenant à la structure se déplace sous le chargement. Son déplacement est caractérisé par le
vecteur déplacement :
→
U (M) = u1 (x1,x2,x3) x 1 + u2 (x1,x2,x3) x 2 + u3 (x1,x2,x3) x 3
Comme on est en théorie des petites perturbations, les composantes u1, u2 et u3 sont "petites" .
1.2.2 Etat des déformations au voisinage d'un point

On se place ici dans le cas des petites déformations.
On démontre que l'état de déformation au voisinage d'un point A est caractérisé par le tenseur
des déformations.
C'est un tenseur du second ordre symétrique qui se déduit du champ des déplacements par la
relation :
1 ∂ ui ∂ uj
εij = ( + )
2 ∂ xj ∂ xi
Dans une base orthonormée il s'écrit en A :
≈  ε 11 ε 12 ε 13 
ε (A) =  ε 22 ε 23 
Sym ε 33 

1.2.3 Allongement unitaire en A et pour une direction q
→
Après déformation, la longueur ds1 du vecteur AA 1 de direction q 1 est devenue ds'1. On peut

alors définir l'allongement relatif en A et pour la direction q 1 (Figure 1. 7) :
ds'1 - ds1
C'est la quantité : e (A, q1 ) =
ds1
ds'1 et ds1 sont des longueurs infiniment petites.
On peut montrer que cette quantité s'exprime à partir du tenseur des déformations par :
≈
e ( A , q1 ) = t
q 1 ε (A) q 1
RDM 5.xx Page 15

Exemple: Calculons l'allongement unitaire dans la direction x1 (1,0,0)
ε 11 ε 12 ε 13  1 
/
e (A, x1) = [1 0 0] ε 12 ε 22 ε 23   0 = ε11
ε 13 ε 23 ε 33  0
On interprète ainsi les termes diagonaux du tenseur des déformations qui représentent les
allongements unitaires dans les trois directions orthonormées.
q2
q'2
A
2
ds 2 A'2
x2 x2
ds1 q1
A1
A
ds'2 θ
Après déformation q'1
A'
1
A ds'1
x1 x1
Figure 1. 7 : Déformation au voisinage d'un point A
1.2.4 Glissement. Distorsion

Après déformation, l'angle droit A1AA2 est devenu A'1A'A'2. Cette variation d'angle droit

s'appelle le glissement ou la distorsion. On peut montrer que le glissement g (A, q1 , q 2 ) se
calcule par :
≈
t
g (A,q1 , q 2 ) = - 2 q 1 q2 ε

Exemple : Calculons la distorsion de l'angle droit ( x 1, x 3)
ε 11 ε 12 ε 13  0 
/ /
g (A, x1, x3) = -2 [1 0 0] ε 12 ε 22 ε 23  0  = - 2 ε 13
ε 13 ε 23 ε 33  1 
On interprète ainsi les termes non diagonaux du tenseur des déformations qui représentent à un
facteur près les distorsions des trois angles droits d'un trièdre.
RDM 5.xx Page 16
1.2.5 Déformations principales et directions principales
De même que pour le tenseur des contraintes, le tenseur des déformations étant réel symétrique,
→
il est diagonalisable, c'est à dire qu'il existe un réel εi et une direction X i telle que :
≈ → →
ε (A) Xi = i Xi ε
→
Ainsi dans une telle direction X i il n'y a pas de glissement mais seulement un allongement.
Dans le repère principal en A, la matrice des déformations s'écrit alors :
≈ ε 1 
ε 
(A) =  ε2 

 ε 3  ( X , X ,X )
1 2 3
1.2.6 Etat de déformations planes

On a affaire à un état plan de déformations parallèlement au plan ( x 1, x 2) si le champ des
→
déplacements U (M) de tout point M peut se mettre sous la forme :
→
U (M) = u1 (x1, x2) x 1 + u2 (x1, x2) x 2 + C x3 x 3
Les composantes u1 et u2 ne sont fonctions que des deux seules variables x1 et x2 , C est une

constante : la composante suivant x 3 est une fonction affine en x3
Cette hypothèse est généralement admise lorsque l'on étudie des pièces cylindriques de

génératrice parallèle à l'axe x 3 , suffisamment longues pour que l'on puisse négliger les effets

aux extrémités, et chargées dans le plan ( x 1, x 2).
C'est le cas des canalisations de transport de fluide par exemple.
Dans ces conditions, le tenseur des déformations s'écrit :
1 ∂ ui ∂ uj
εij = ( + )
2 ∂ xj ∂ xi
Comme u1 et u2 ne dépendent pas de x3 , ε13 = ε23 = 0
et ε33 = C
∂ u1 ∂ u2 1 ∂ u1 ∂ u2
ε11 = , ε22 = et ε12 = ( + )
∂ x1 ∂ x2 2 ∂ x2 ∂ x1
Dans la suite on supposera C = 0.
RDM 5.xx Page 17
Le tenseur des déformations s'écrit alors dans une base quelconque :
≈ ε 11 ε 12 0 
ε (A) = ε 12 ε 22 0 
 
 0 0 0  ( x , x , x )
1 2 3

L'axe x 3 est donc direction principale et l'allongement unitaire correspondant est nul.
1.3 Loi de comportement

La linéarité de la loi de comportement de l'élasticité se traduit par la linéarité de la loi qui relie
tenseur des contraintes et tenseur des déformations.
Soit en notation indicielle :
σij = λ ε kk δij + 2 G ε ij
Dans cette relation : δij est le tenseur de Kronecker
δij = 1 si i = j , δij = 0 si i ≠ j
Dans certains manuels on note G : µ
λ et G sont les coefficients de Lamé, constants pour un matériau donné.
εkk = ε11 + ε22 + ε33 est le premier invariant du tenseur des déformations
Inversement, on peut exprimer le tenseur des déformations à partir de celui des contraintes :
1 + ν ν
ε ij =
E
σij -
E
σ kk δ ij
Dans cette relation :
• E est le module de Young et ν le coefficient de Poisson, constants pour un matériau

donné.
• εkk = ε11 + ε22 + ε33 est le premier invariant du tenseur des contraintes
Les relations entre les différents coefficients d'élasticité sont les suivantes :
Eν E
λ = , G=
(1- 2 ν) (1 + ν) 2 (1 + ν)
G est le module de Coulomb.
Si ν = 0,25 alors G = E / 2.5 , c'est le cas pour les matériaux ductiles.
Pour l'acier "doux" ( S235 par exemple) : E ≈ 200 000 MPa et G ≈ 80 000 MPa
RDM 5.xx Page 18
En décomposant sur les axes on obtient :
σ11 ν σ ν
ε11 = - ( σ 22 + σ 33) , ε 22 = 22 - (σ + σ11 )
E E E E 33
σ 33 ν σ12 σ13 σ
ε 33 = - (σ + σ11 ) , ε12 = , ε13 = , ε 23 = 23
E E 22 2G 2G 2G
1.4 Critères de limite élastique

On supposera dans la suite que la limite élastique en traction simple est égale à la limite
élastique en compression simple (matériaux ductiles). Soit σe cette limite.
On connaît bien le comportement d'un matériau dans le cas d'une sollicitation de traction
simple. Cette connaissance est liée à l'essai de traction simple statique.
Soit un état de contrainte complexe caractérisé en un point A par les trois contraintes
principales
σ1, σ2 et σ3.
Existe-t-il un moyen de savoir si, en ce point, la limite élastique est dépassée ?
On peut répondre à cette question par l'affirmative. On définit pour cela une contrainte de
traction simple σg présentant le même danger de dépassement de limite élastique que l'état de
contrainte complexe ( σ1, σ2 ,σ3).
σg est appelée contrainte équivalente. Bien entendu, dans le cas général, σg est une
contrainte fictive que l'on ne rencontre pas dans la pièce contrainte.
Il n'y a pas unicité du critère de limite élastique. Au cours de l'histoire de la mécanique des
milieux continus déformables, plusieurs critères ont été proposés. Certains sont plus ou moins
bien vérifiés en fonction du type de matériau sollicité et du type de sollicitation.
A l'heure actuelle les logiciels d'élasticité prennent en compte surtout les critère de Tresca et de
Von Mises que nous allons expliciter.
1.4.1 Critère de Tresca ou du cisaillement maximal

Pour ce critère, l'état limite est atteint lorsque la contrainte de cisaillement maximal admet la
valeur seuil τe déterminée par l'essai de torsion.
A l'aide de la représentation de Mohr, on sait déterminer le cisaillement maximal, il suffit de

tracer le plus grand des trois cercles de Mohr.
RDM 5.xx Page 19
τ
τMax
τ
T(A,n)
σ3 σ2 σ σ1 σ
Figure 1. 9 : Tricercle de Mohr et contrainte tangentielle maximale
a) Dans le cas particulier d'une sollicitation de traction simple :
σ1 seul est différent de zéro. Alors: τmax = σ1 / 2
τMax
σ
Ο σ1
Figure 1. 10: Traction uniaxiale
Le critère de Tresca impose : τmax = σ1/2 ≤ τe
Si σ1 atteint sa limite σe on obtient: τe = σe / 2
Ainsi le critère de Tresca impose cette dernière relation entre τe et σe
Rappel :
• σe est la limite élastique en traction simple.
• τe est la limite élastique en cisaillement simple.
RDM 5.xx Page 20
b) Dans le cas particulier d'une sollicitation en contraintes planes
Soit un état de contraintes planes caractérisé par les deux contraintes principales σ1 et σ2
(σ3 = 0).
On peut représenter cet état de contraintes dans le plan des contraintes (σ1, σ2). En effet à tout
état de contraintes (σ1 , σ2) on fait correspondre un point caractéristique M (σ1, σ2). Le point
M peut, à priori, décrire tout le plan des contraintes.
σ2
Ο σ1
Figure 1. 11 : Etat de contraintes planes
Distinguons plusieurs cas :
1°) σ1 et σ2 admettent le même signe
Dans ces conditions M appartient au premier ou au troisième quadrant.
a) Si σ1 > σ2
La plus grande scission est :
τmax = 0,5  σ1- σ3 or σ3 = 0 donc: τmax = 0,5  σ1 
et le critère de Tresca impose : 0,5  σ1  ≤ τe
Soit σ1 ≤ 2 τe = σe . Ainsi dans ces conditions:  σ1  ≤ σe
Si σ1 > 0 alors σ1 < σe

Si σ1 < 0 alors σ1 > - σe
b ) Si σ1 < σ2
La plus grande scission est : τmax = 0,5 σ2 

Ainsi dans ces conditions :  σ2  ≤ σe
Si σ2 > 0 alors σ2 ≤ σe
Si σ2 < 0 alors σ2 ≥ - σe
RDM 5.xx Page 21
τ
σ
σ3= 0 σ2 σ1
Figure 1. 12: σ1 et σ2 de même signe
2°) σ1 et σ2 admettent des signes opposés
Dans ces conditions M appartient au deuxième ou au quatrième quadrant.

τ max = 0,5 σ1 - σ2 
Si σ1 > σ2 alors τmax = 0,5 (σ1 - σ2) soit σ1 - σ2 ≤ σe
Si σ1 < σ2 alors τmax = 0,5 (σ2 - σ1) soit σ2 - σ1 ≤ σe
σ
σ2 σ3 = 0 σ1
Figure 1. 13 : σ1 et σ2 de signes différents
Conclusion :
Pour appliquer le critère de Tresca, il faut bien connaître les signes de σ1 et σ2 :
- si σ1 et σ2 ont un même signe σ1 < σe et σ2 < σe ;

- si σ1 et σ2 ont des signes différents: σ1 - σ2 < σe
La zone admissible pour le point M est l'intérieur du polygone de la figure 1. 14.
RDM 5.xx Page 22
σ2
σe
σ2 σ1 + σe
Polygone de Tresca
σe σ1
σe
σe σ2 σ1 σe
Figure 1. 14 : Polygone de Tresca
1.4.2 Critère du plus grand travail de distorsion. Critère de

Von Mises
Pour ce critère, l'état limite est atteint lorsque l'énergie de distorsion par unité de volume est
égale à l'énergie de distorsion unitaire limite du matériau.
L'énergie de distorsion par unité de volume s'exprime en fonction des contraintes principales
par:
dW(f ) 1 + ν
= [ ( σ1 - σ2 )2 + ( σ1 - σ3 )2 + ( σ2 - σ3 )2 ]
dv E
Dans le cas de la traction simple, seule σ1 ≠ 0 :
dW(f ) 1 + ν
2 σ1
2
=
dv E
Ce critère devant être valable quel que soit l'état de sollicitation, on doit donc avoir :
[ ( σ1 - σ2 )2 + ( σ1 - σ3 )2 + ( σ2 - σ3 )2 ] ≤ 2 σe ²
La contrainte de traction simple équivalente σg à l'état de contrainte complexe est alors en

fonction des contraintes principales telle que :
1
σg ² = ( σ1 - σ2 )2 + ( σ1 - σ3 )2 + ( σ2 - σ3 )2
2
et, en fonction des contraintes non principales :

1
σg ² = ( σ11 - σ22 )2 + ( σ11 - σ 33 )2 + ( σ22 - σ33 )2 + 6 ( σ12 ² + σ23 ² + σ13 ²)
2
1.4.3 Critères exprimés dans le cas de contraintes planes

Seules les contraintes σ11, σ22 et σ12 sont non nulles.
1.4.3.1 Critère de Tresca
σg = [(σ11 - σ22)2 + 4 σ122 ]1/2
Si de plus σ22 = 0 (flexion-torsion) alors : σg = ( σ112 + 4 σ122 ) 1/2
RDM 5.xx Page 23
1.4.3.2 Critère de Von Mises
En fonction des contraintes principales, il reste :
σg = ( σ12 + σ22 - σ1 σ2 ) 1/2
L'état limite est atteint pour: σ12 + σ22 - σ1 σ2 = σe2
C'est l'équation d'une ellipse dont le grand axe est incliné de 45° dans le plan ( σ1, σ2 )
En fonction des contraintes non principales, il vient :
σg = ( σ112 + σ222 - σ11 σ12 + 3 σ122 )1/2
Si de plus σ22 = 0 (flexion-torsion) alors :
σg = ( σ112 + 3 σ122 )1/2
σ2
σe
σ 2 σ1 + σe Ellipse de Von Mises
Polygone de Tresca
σe σ1
σe
σe σ2 σ1 σe
Figure 1. 15 : Représentation graphique des critères de limite élastique
σe est la limite élastique du matériau considéré.
Conclusion :
Il apparaît clairement que les critères de Tresca et de Von Mises ne donnent pas tout à fait les
mêmes résultats ( coefficients 3 ou 4 en flexion-torsion, notamment).
L'expérience montre que le critère de Von Mises est souvent plus satisfaisant que celui de
Tresca. De plus le critère de Von Mises est d'emploi plus facile que celui de Tresca car il se
traduit par une formule unique quadratique, donc sans problème de signe.
1.5 Types particuliers de problème d’élasticité

RDM 5.xx Page 24
1.5.1 Contraintes planes du plan (x1,x2)
C'est le cas de certaines plaques soumises à des forces parallèles à leur plan moyen.
Les contraintes et les déformations sont indépendantes de x3.
Le tenseur des contraintes admet la forme :
≈ σ 11 σ 12 0 
σ (A) = σ 
 12 σ 22 0 
 0 0 0  ( x , x , x )
1 2 3
Le tenseur des déformations admet la forme :
≈ ε 11 ε 12 0 
ε (A) = ε 12 ε 22 0  ε33 =
-ν
1 - ν
( ε11 + ε22 )
 0 0 ε 33  ( x , x , x )
1 2 3
Avec :
1
ε11 = ( σ11 - ν σ22 )
E
1
ε22 = ( σ22 - ν σ11 )
E
1 + ν
ε12 = σ12
E
1.5.2 Déformations planes

C'est le cas de certains solides cylindriques d'axe x 3, longs et soumis à des forces de surface

et de volume perpendiculaire à x 3. Ces forces sont indépendantes de la coordonnée x 3
≈ ε 11 ε 12 0
ε (A) = ε 12 ε 22 0
 0 0 0 ( x , x , x )
1 2 3
RDM 5.xx Page 25
≈ σ 11 σ 12 0 
σ (A) = σ 
 12 σ 22 0 
 0 0 σ 33  ( x , x , x )
1 2 3
σ33 = ν ( σ11 + σ22)
Avec :
1 + ν
ε11 = (1− ν ) σ11 - ν σ 22
E
1 + ν
ε22 = (1 − ν) σ22 - ν σ11
E
1 + ν
ε12 = σ12
E
1.5.3 Problèmes axisymétriques méridiens

Par définition, la structure, le chargement et les forces de liaisons admettent un axe de

révolution O x 3.
De plus toutes les forces sont situés dans des plans méridiens (pas de forces
circonférentielles).
Tous les paramètres sont indépendants de la coordonnées circonférentielle θ. On traite ce type

de problème en coordonnées cylindriques.
≈  ε r 0 ε rz 
ε (A) =  0 εθ 0 
 
ε rz 0 ε z  ( r ,θ ,z )
≈  σ r 0 σ rz 
σ (A) =  0 σθ 0 
 
σ rz 0 σ z  ( r ,θ , z )
RDM 5.xx Page 26
2. Fondement mécanique de la méthode des
éléments finis
Ce bref exposé a pour but de présenter brièvement à des utilisateurs mécaniciens potentiels
les idées essentielles de la méthode des éléments finis, appliquée à l'élasticité linéaire. Il ne
faut pas chercher dans ces quelques lignes une rigueur ni mécanique, ni mathématique (un
ouvrage entier serait alors nécessaire). On rappelle les points essentiels de la méthode, tant
mathématiques que mécaniques. C'est un premier "point d'entrée" relativement à cette
méthode. Les esprits curieux pourront se référer à la bibliographie fournie.
Afin de simplifier les notations, on se limite à un exposé en élasticité plane.
2.1 Notations
La structure occupe le domaine plan (D), délimitée par la surface (S ).
≈
σ (M) : tenseur des contraintes ;
≈
ε (M) : tenseur des déformations ;
Notations vectorielles :
→
U (M) : champ des déplacements réels ;
→
U *(M) : champ de déplacements virtuels ;
→
f (M) : champ des efforts volumiques ;
→ ≈
T (M,n) = σ (M) n : vecteur contrainte en M pour la direction n ;
→ →
→
Q (P): champ des efforts surfaciques sur (S), P∈(S) ;
Notations matricielles des vecteurs et des matrices :
{u} : vecteur déplacement associé à un point courant M ;

{d} : vecteur déplacement associé aux noeuds d'un élément ;
{D} : vecteur déplacement des noeuds de la structure entière ;
[E] : matrice relative à la loi de comportement ;
[B] : matrice associée à la relation {ε} = [B] {d} ;
[N] : matrice d'interpolation associée à la relation {u} = [N] {d} ;
[k] : matrice de rigidité associée à un élément ;
[K] : matrice de rigidité de la structure ;
{r} : vecteur d'effort appliqué aux noeuds d'un élément ;
{R} : vecteur d'effort appliqué aux noeuds de la structure ;
W : énergie de déformation de la structure ;

πP : énergie potentielle de la structure ;
m : nombre d'éléments de la structure discrétisée ;
RDM 5.xx Page 27
2.2 Problème d'élasticité. Equations d'équilibre
Dans tout ce qui suit, on considère que la structure est en équilibre (Figure 2. 1).
Les conditions aux frontières sont généralement de deux types :

→
- des forces surfaciques Q (P) sont imposées sur une partie SF de (S) ;
→
- des déplacements U O(P) sont imposés (appuis) sur une partie SU de (S) ;
On suppose que : S = SU ∪ SF et SU ∩ SF = Φ
Les conditions aux limites à vérifier sur S sont alors :
→ →
T (P, n ) = Q (P) pour tout point P de SF ; (1)
→ →
U (P) = U O(P) pour tout point P de SU ;
≈ →
σ (M) étant le tenseur des contraintes en M et f (M) le champ des efforts volumiques, les
équations d'équilibre s'écrivent :
σij,j + fi = 0 i et j ∈[1,3] (2)
Dans (2) on emploie la convention de l'indice muet.
La relation (2) conduit à trois équations aux dérivées partielles :
∂ σ 11 ∂ σ 12 ∂ σ 13
+ + + f1 = 0
∂ x1 ∂ x2 ∂ x3
∂ σ 21 ∂ σ 22 ∂ σ 23
+ + + f2 = 0
∂ x1 ∂ x2 ∂ x3
∂ σ 31 ∂ σ 32 ∂ σ 33
+ + + f3 = 0
∂ x1 ∂ x2 ∂ x3
RDM 5.xx Page 28
Q(A)
y
σi,j + f i =0 n
A SF
M
(D)
SU
x
Figure 2. 1 : Problème d'élasticité
2.3 Loi de comportement

En élasticité, en petites déformations et pour les matériaux isotropes linéaires elle s'écrit :
≈ ≈
σ (M) = [E] ε (M) (3)
En déformation plane :
 ν 
 1 1-ν
0 
σ 11    ε 11 
  ν  
σ 12  =  1 0  ε 12 
σ  1-ν   2ε 
 22  1 − 2ν   22 
 0 0
 2(1-ν ) 
et en contraintes planes :
 
σ 11  1 ν 0  ε 11 
  E    
σ 12  = ν 1 0  ε 12 
σ  (1+ν ) (1-ν )  1 −ν   2ε 
 22  0 0   22 
 2 
2.4 Energie de déformation élastique

En élasticité, en petites déformations et pour les matériaux isotropes linéaires, on appelle
énergie de déformation élastique par unité de volume en M, la fonction w(M) telle que :
RDM 5.xx Page 29
1
w(M) = σij(M) εij(M) (4)
2
Pour toute la structure, occupant un domaine D on obtient :
W= ∫ w(M) dV (5)
D
2.5 Théorème d'unicité

≈ →
La théorie de l'élasticité permet de démontrer que la solution ( σ (M), U (M)) vérifiant les
conditions aux limites définies par (1) est unique si l'énergie de déformation élastique (5) est
positive.
2.6 Champ de déplacement virtuel admissible
2.6.1 Définition
→
Un champ de déplacement U *(M) défini sur D est dit cinématiquement admissible, s'il est
continûment dérivable dans D et s'il vérifie les conditions de déplacement sur SU :
→ →
U *(P) = U O(P) quel que soit P appartenant à SU
2.6.2 Conséquence
→
U *(M) étant continûment dérivable, on peut à partir de ce champ calculer le tenseur des
≈* ≈
déformations ε (M) et par les lois de comportement, le tenseur des contraintes σ *(M). Mais
≈
ce tenseur des contraintes ne vérifie pas les équations d'équilibre et σ *(M). n ne vérifie pas la
→
→
condition aux frontière sur SF, sinon U *(M) serait le champ des déplacements réels en vertu
du
théorème d'unicité.
2.7 Energie potentielle d'un système élastique
On appelle énergie de potentielle de la structure occupant un domaine D, l'énergie définie

par :
RDM 5.xx Page 30
→ → → → →
πP (U) = W(U) - ∫D
f (M) . U (M) dV - ∫
S
T(P, n ) . U (P) dS

→
πP (U) est l'énergie potentielle de la structure associée au champ de déplacement U(M) ;
W(U) est l'énergie de déformation élastique de la structure associée au champ de

→
déplacement U(M) ;
→
la première intégrale représente le travail développé par les efforts volumiques f (M), et
moins cette intégrale représente l'énergie potentielle associé à ces efforts ;
→ →
la deuxième intégrale représente le travail développé par les efforts surfaciques T (P, n) , et
moins cette intégrale représente l'énergie potentielle associé à ces efforts ;
Considérons un solide déformable plus les charges supportées par celui-ci. Le principe des
travaux virtuels permet d'affirmer que :
Parmi tous les champs de déplacement virtuels admissibles, ceux qui satisfont les
équations d'équilibre rendent l'énergie potentielle extrémale.
Si cet extremum est un minimum, alors l'équilibre est stable.
2.7.1 Système à un degré de liberté

Exemple : Ressort de raideur k soumis à une charge P.
Le déplacement sous la charge P est x (un degré de liberté).

Dans ce cas l'énergie potentielle du système est :
1
πP = k x² - P x
2
Le déplacement virtuel δx qui conduit à l'équilibre statique est tel que :
δπP = 0
Soit : kx-P=0
P
La position d'équilibre est alors : xéq=
k
RDM 5.xx Page 31
W 2
1/2 k x
2
1/2 k x - Px
x
x éq
x
P
-Px
Figure 2. 2 : Energie potentielle d'un ressort
2.7.2 Système à plusieurs degrés de liberté

Si la configuration d'un système dépend de Nd degrés de liberté, l'énergie potentielle dépend
de ces Nd degrés de liberté et :
πP = πP (D1, D2,......, Dn)
∂πP ∂πP ∂πP

alors : δ π p = δD1 + δD2 + ................+ δDn
∂ D1 ∂ D2 ∂ Dn
t
∂ π P 
soit, en notation matricielle : δ π p =   {δ D }
 ∂ D1 
Si l'on cherche un extremum de πP, on doit avoir δπP = 0, quel que soit le champ de
déplacements infinitésimal {δD}, soit:
∂πP ∂ π P 
= 0 pour i ∈ [ 1, Nd ] , soit  = {0}
∂ D1 ∂ D 
On obtient ainsi un système algébrique de Nd équations avec Nd inconnues.
2.7.3 Système continu

Un système continu comporte une infinité de degrés de liberté. Par exemple en élasticité, à
tout point matériel M du solide on peut associer un vecteur déplacement :
→
U(M) = u1 x 1 + u2 x 2 + u3 x 3
soit trois degrés de liberté.
RDM 5.xx Page 32
Comme le nombre de "points" est infini , le nombre de degrés de liberté du système est infini.
La résolution des équations aux dérivées partielles (2) d'équilibre , soumises aux conditions
aux limites (1) n'est pas chose facile !
L'idée de la théorie des éléments finis est donc de réduire le système à un nombre Nd fini de
degrés de liberté, c'est à dire d'effectuer le calcul en un nombre fini de points appelés noeuds
du système, on obtient ainsi une solution approchée à partir de la résolution de Nd équations
algébriques. La solution est approchée en ce sens que les équations d'équilibre ne seront pas
rigoureusement satisfaites en tout point M du solide, le nombre de degrés de liberté étant fini.
Cette solution sera d'autant moins approchée que ce nombre de degrés de liberté sera
important.
2.8 Approximation par éléments finis

Raisonnons pour plus de facilité en déformation plane.
→
En élasticité, on considère généralement que le champ des déplacements U(M) est l'inconnue
première du problème.
→
Soit : U(M) = u1(x1,x2) x 1 + u2(x1,x2) x 2
x2
De m
O x1
Figure 2. 3 : Découpage en éléments finis
Appelons u l'une des deux composantes u1 ou u2. La méthode d'approximation par éléments
finis simplifie la construction de la fonction approchée u et s'adapte bien au calcul sur
ordinateur. Elle consiste à :
- définir un ensemble De de sous-domaines de D , sans recouvrements ni intersections

(l'exposant e signifie élément) (figure 2. 3) ;
- définir un champ de déplacements cinématiquement admissible (fonction approchée)

ue (x1, x2) différente sur chaque élément par la méthode d'approximation nodale.
RDM 5.xx Page 33
L'approximation nodale utilisée admet de plus les particularités suivantes :
- cette approximation ue (x1,x2) ne fait intervenir que les noeuds situés sur De et sur
sa frontière ;
- les fonctions ue (x1,x2) sont continues sur De et elles satisfont à des conditions de
continuité entre les différents sous-domaines De.
2.8.1 Définitions
Les sous domaines De sont appelés des éléments ;
Les points en lesquels la fonction u (x1, x2) coïncide avec la fonction exacte ue (x1, x2) sont
les noeuds d'interpolation ;
Les coordonnées (x1, x2) de ces noeuds sont les coordonnées nodales ;
Les valeurs des déplacements (d1, d2) au noeud d'interpolation considéré sont les variables
nodales.
La solution est obtenue quel que soit le point M en se servant de fonctions d'interpolation N,
soit :
{u} = [ N ] { d }

→
{u} est le vecteur comportant les deux composantes du vecteur déplacement U(M) ;
{d} est le vecteur qui comporte autant de composantes que l'élément possède de degrés de
liberté nd ;
enfin [N] est une matrice rectangulaire à 2 lignes et nd colonnes appelée matrice
d'interpolation. Chaque terme de cette matrice est fonction des coordonnées des noeuds de
l'élément et des coordonnées du point M.
2.9 Méthode des éléments finis en élasticité, conduite à

partir des déplacements
A partir du champ des déplacements, on peut construire le tenseur des déformations car :
ε ij =
u i, j + u j,i
2
On obtient sous forme condensée:
{ε} = [B] {d }
On obtient ainsi le potentiel d'un élément :
1 t  
{d}  ∫ [ B] [ E ][ B] dV  {d} - {d} ∫ [ N ] {f } dV- {d} ∫ [ N ] {Q} dS
t t t t t
πPe =
2  De  De Se
RDM 5.xx Page 34
L'énergie potentielle totale du système est la somme des énergies potentielles des éléments.
Faisons intervenir de plus le potentiel des charges P concentrées, directement appliquées aux
noeuds. On obtient alors pour les m éléments :
∑π - {D} {P}
t
πP = Pe
i =1
On appelle maintenant {D} le vecteur des déplacements de la structure totale. Il comporte 2

Nd composantes si la structures comporte Nd noeuds.
πPe=
 m
1 t  m    
{D} ∑ ∫ [ B] [ E ][ B] dV  {D} - {D} ∑  ∫ [ N ] {f } dV- ∫ [ N ] {Q} dS   - {D} {P}
t t t t t


2 
 i=1 De   De
 i=1
 Se  
La configuration d'équilibre est obtenue si δπP = 0, quel que soit le champ de déplacements
infinitésimal {δD}, soit :
∂πP ∂ π P 
= 0 pour i ∈ [ 1, Nd ] , soit  = {0}
∂ Di ∂ D 
On obtient ainsi un système algébrique de Nd équations avec Nd inconnues :
m   m   
∑ ∫ [ ] [ ][ ] { } ∑  ∫ [ N ] {f } dV+ ∫ [ N ] {Q} dS   + {P}
t t t
B E B dV  D =
 i=1 De   i=1  De Se  
∫ [ B] [ E ][ B] dV
t
Appelons : [k] =
De
chacune des m matrices de rigidité élémentaire du premier membre de la relation

précédente.
∫ [ N ] {f } dV + ∫ [ N ] {Q} dS
t t
De même notons : {r} =
De Se
le vecteur du second membre relatif à chaque élément.
On peut alors écrire de façon simplifiée :
m  m 
∑ [ k ] {D} = ∑ {r} + {P}
 i=1   i=1 
Cette relation peut encore s'écrire de façon plus concise :
[K] {D} = {R}
[K] s'appelle la matrice rigidité de la structure totale (connue après calcul) ;

{D} est le vecteur déplacement (inconnu) de la structure complète ;
{R} est le vecteur (connu) des efforts extérieurs appliqués à la structure.
RDM 5.xx Page 35
On obtient donc un système linéaire de Nd équations à Nd inconnues que l'on sait bien
résoudre numériquement par des méthodes de type Gauss ou par des méthodes itératives pour
les grands systèmes.
2.10 Application à l'élément triangulaire à trois noeuds

→
On suppose que les efforts f (M) de volume sont nuls. On détaille la méthode de calcul sur ce
type d'élément car les calculs engendrés sont simples, la méthode exposée n'en reste pas
moins générale. Cependant, les résultats issus d'un calcul avec ce type d'élément sont à
regarder de très près.....
vn
y
n
yn un vm
ym m
um
v
k
yk k
u
k
xk xn xm x
Figure 2. 4: Elément triangulaire à 3 noeuds
2.10.1 Construction de la loi d'interpolation

Comme on dispose de trois noeuds d'interpolation on peut écrire :
u (x,y) = α1 + α2 x + α3 y (20)
v (x,y) = α4 + α5 x + α6 y
αi sont des coefficients qu'il s'agit de déterminer en fonction des déplacements des noeuds.
Les relations ci-dessus montrent que l'approximation des déplacements est linéaire sur un
élément.
RDM 5.xx Page 36
 F1(x,y)  1 x y 0 0 0 
En notant [F] la matrice d'interpolation : [F] =   =  
 F2(x,y)  0 0 0 1 x y 
u (x,y) 
et {U} le vecteur des déplacements en un point M du domaine {U} =  
 v (x,y) 
On peut écrire matriciellement :
α 1 
α 
 2
u (x,y)  1 x y 0 0 0  α 3 
 =    
 v (x,y)   0 0 0 1 x y  α 4 
α 5 
 
α 6 
Il faut écrire ensuite que les déplacements en un point quelconque sont construits à partir de
ceux des noeuds (interpolation) :
uk = α1 + α2 xk + α3 yk
vk = α4 + α5 xk + α6 yk
um = α1 + α2 xm + α3 ym
vm = α4 + α5 xm + α6 ym
un= α1 + α2 xn + α3 yn
vn = α4 + α5 xn + α6 yn
Soit matriciellement :
 u k   1 xk yk 0 0 0  α 1 
 v  0 0 0 1 xk yk  α 2 
 k 
 um   1 xm ym 0 0 0  α 3 
 =    
 vm  0 0 0 1 xm ym  α 4 
un  1 xn yn 0 0 0  α 5 
     
 vn  0 0 0 1 xn yn  α 6 
et, sous forme condensée :
{d} = [A] {α} (21)
En inversant cette relation matricielle, on obtient :
{α} = [A]-1 {u} = [C] {d} ( car les déplacements des noeuds sont notés {d} (22)
RDM 5.xx Page 37
α 1   c11 0 c12 0 c13 0   uk 
α  c   
 2  21 0 c22 0 c23 0   vk 
α 3  1  c31 0 c32 0 c33 0   um 
Soit :  =    
α 4  2 ∆e  0 c11 0 c12 0 c13   vm 
α 5   0 c21 0 c22 0 c23   un 
     
α 6   0 c31 0 c32 0 c33   vn 
Les αi sont maintenant déterminés en fonction des coordonnées des noeuds k, m et n et des
déplacements de ceux-ci.
On note ∆e l'aire de l'élément qui est définie par :
1 xk yk
1
∆e = det 1 xm ym
2
1 xn yn
∆e est l'aire du triangle k, m, n. Il apparaît donc intéressant d'avoir un triangle d'aire maximale
car ∆e intervient au dénominateur ( en tout état de cause non nul ! ).
La matrice [C] contient les coordonnées des noeuds k, m et n ( décrit dans le sens
trigonométrique).
c11= xm yn - xn ym c12= xn yk - xk yn c13 = xk ym - xm yk
c21= ym - yn c22= yn- yk c23 = yk - ym
c31= xn - xm c32= xk- xn c33 = xm - xk (23)
uk 
v 
 k
 u (x,y)  1 1 x y 0 0 0   um 
 =   [C]  
 v (x,y)  2 ∆e  0 0 0 1 x y   vm 
 un 
 
 vn 
Ces relations définissent le champ de déplacement linéaire sur un élément.
RDM 5.xx Page 38
2.10.2 Tenseur des déformations approchées
En élasticité plane :
ε11 = u1,1 ε22 = u2,2 et ε12 = (u1,2 + u2,1)/2 (24)
Compte tenu des relations (20), il vient :
ε11 = α2 ε22 = α6 et ε12 = ( α3 + α5 )/2 (25)
Ainsi, à l'intérieur d'un élément les déformations et donc les contraintes sont constantes.
u k 
v 
ε 11  0  
k
c21 0 c22 0 c23
ε  = 1  0 c  um 
 22  2 ∆e  31 0 c32 0 c33   
vm
 2ε 12  c31 c21 c32 c22 c33 c23   
 un 
 
 vn 
∆e est l'aire du triangle k, m, n.
Soit en remplaçant les cij par leurs valeurs :

u k 
v 
ε 11   y m -y n 0 y n -yk 0 y k -y m 0   k
ε  =  0  um 
x m -x k 
1
 22  2 ∆e 
x n -x m 0 x k -x n 0  
 2ε 12   x n -x m y m -y n x k -x n y n -y k x m -x k yk -y m   vm 
 un 
 
 vn 
Le champ des déformations est donc constant à l'intérieur de l'élément.
Soit : {ε} = [B] {d} (27)
Notons que la matrice [B] est constante à l'intérieur d'un élément.
2.10.3 Loi de comportement

≈ ≈
Elle s'écrit : σ = [E] ε (28)
En contraintes planes :
RDM 5.xx Page 39
 
1 ν 0 
E  
[E] = ν 1 0  (29)
1− ν² 
1 −ν 
0 0 
 2 
et en déformations planes :
 
ν (1 −ν ) ν 0 
E  
[E] =  ν 1 −ν 0 
(1 + ν )(1 − 2 ν )
 1 −ν 
 0 0 
 2 
 E11 E12 0 

Dans les deux cas on peut écrire : [E] = E12 E 22 0

 
 0 0 E 33 
Remarquons que l'on passe de (29) à (30) en substituant :
E ν
à E et à ν
1− ν² 1− ν
2.10.4 Matrice de rigidité de l'élément

La matrice de rigidité de l'élément vaut :
∫ [ B] [ E ][ B] dV
t
[k] =
De
Dans ce cas, les matrices [B] et [E] étant constantes, l'intégration donne simplement :
[k] = [ B] [ E ][ B] dV = [B]t[E][B] V
t
∫
De
V est le volume de l'élément.
2.10.5 Vecteur d'effort
∫ [ N ] {Q} dS
t
Il s'écrit : {r} =
Se
Lorsque l'on prend en compte l'ensemble des éléments, ce vecteur s'annule sauf pour les
éléments dont les cotés sont des frontières chargées par un chargement extérieur.
Remarque :
RDM 5.xx Page 40
On a déjà fait remarquer que les contraintes à l'intérieur d'un élément sont uniformes. Ce type
d'élément est donc à proscrire dans les zones à fort gradient de contraintes. Par contre, en
traction, cet élément peut s'avérer intéressant.
2.11 Application à une poutre en flexion simple

Soit la poutre ci-dessous, en flexion simple, encastrée à son extrémité gauche et supportant
une charge verticale à droite.
Considérons la poutre d'épaisseur constante e.
Le processus de résolution est le suivant :
2.11.1 Décomposition en éléments finis

On définit le contour de la pièce et l'on décompose l'intérieur en un certain nombre d'éléments
(ici triangulaires).
F
2 4 6
(III) (IV)
H
(I) (II)
1 3 5
x
L L
Figure 2. 5
L'élément générique est défini géométriquement par la position de ses trois noeuds k, m et n :
vn
y
n
yn un vm
ym m
um
v
k
yk k
u
k
xk xn xm x
Figure 2. 6
RDM 5.xx Page 41
La poutre est ici séparée en 4 éléments ( I à IV ) et comporte 6 noeuds.
2.11.2 Calcul de la matrice rigidité de chaque élément

Si S est l'aire d'un élément, la matrice de rigidité d'un élément est :
[k] = [B]t [E] [B] e S
Faisons une étude en élasticité plane alors la matrice d'élasticité est :
 E11 E12 0 

[E] = E12 E 22 0

 
 0 0 E 33 
La matrice [B] est issue de la relation entre le vecteur déplacement des noeuds et le vecteur des
déformations : {ε} = [B] {d}
Cette relation s'explicite par :

 ul 
v 
ε 11   y m -yn 0 yn -yl 0 y l -y m 0   l
  1   um 
ε 22  =  0 x n -x m 0 x l -x n 0 x m -x l   
 2ε  2 ∆e  x -x y m -y n x l -x n y n -yl x m -x l yl -y m   vm 
 12   n m  un 
 
 vn 
∆e est l'aire du triangle k, m, n . Si l'on passe de k à m et de m à n par rotation dans le sens
trigonométrique. Elle se calcule par :
1 xk yk
1
∆e = det 1 xm ym
2
1 xn yn
Dans l'exemple, l'aire commune à tous les triangles est ∆e = H L / 2 soit 2 ∆e = H L
Pour l'élément choisi, la matrice [k] ne comporte que les coordonnées des noeuds et les
coefficients d'élasticité.
La matrice de rigidité traduit la relation qui existe entre les efforts appliqués par les noeuds
sur les éléments.
Xk  uk 
Y  v 
 k  k
Xm  um 
Cet effort sera noté : { Fe } =   et les déplacements associés : {d} =  
Ym   vm 
Xn  un 
   
Yn   vn 
RDM 5.xx Page 42
Ainsi :
 Xk    uk 
Y    v 
 k    k
 Xm   Ke  um 
  =    
 Ym     vm 
 Xn    un 
     
 Yn     vn 
Nous allons faire tout d'abord deux remarques essentielles pour la suite.
Faisons une remarque relative à la matrice élémentaire [Ke] associée à l'élément ( k, m, n).
On peut décomposer la matrice (6,6) en 9 blocs (2,2) dont la lecture est la suivante :
kmm : coefficients des forces exercées sur l'élément, au voisinage du noeud m, dues
aux déplacements du noeud de numéro m
kmn : coefficients des forces exercées sur l'élément, au voisinage du noeud m, dues
aux déplacements du noeud de numéro n
 Xk  uk 
Y  v 
 k  kkk kkm kkn   k
 Xm  k  um 
  =  mk kmm kmn   
 Ym   knk knm knn   vm 
 Xn  un 
   
 Yn   vn 
Matrice de rigidité de l'élément I :
 X3   u3 
 Y3   v3 
   k33,I k32,I k31,I   
 X2     u2 
  =  k23,I k22,I k21,I   
 Y2   k13,I k12,I k11,I   
v2
 X1   u1 
   
 Y1   v1 
Matrice de rigidité de l'élément II :
RDM 5.xx Page 43
 X3   u3 
 Y3   v3 
   k33,II k35,II k34,II   
 X5     u5 
  =  k53,II k55,II k54,II   
 Y5   k43,II k45,II k44,II   v5 
 X4   u4 
   
 Y4   v4 
Matrice de rigidité de l'élément III :
 X3   u3 
 Y3   v3 
   k33,III k34,III k32,III   
 X4     u4 
  =  k43,III k44,III k42,III   
 Y4
 v4 
 k23,III k24,III k22,III 
 X2   u2 
   
 Y2   v2 
Matrice de rigidité de l'élément IV :
 X5   u5 
 Y5   v5 
   k55,IV k56,IV k54,IV   
 X6     u6 
  =  k65,IV k66,IV k64,IV   
 Y6
 v6 
 k45,IV k46,IV k44,IV 
 X4   u4 
   
 Y4   v4 
2.11.3 Assemblage de la matrice

Cet assemblage consiste à traduire l'équilibre des noeuds et donc à expliciter la relation
matricielle :
 4   4 
- ∑ [ k ] {D} = ∑ {r} + {P}
 i=1   i=1 
Le signe moins apparaît devant [k] car il s'agit maintenant des actions des éléments sur le
noeud isolé
Quand on isole un noeud i, il intervient les actions de tous les éléments partageant ce
noeud.
RDM 5.xx Page 44
Exemple : Noeud numéro 3 :
III
II
I
Noeud 3
Figure 2. 7
Les éléments II, III et I agissent sur le noeud numéro 3 :
Il n'y a pas d'effort extérieur. Son équilibre se traduit par : FII→3 + FIII→3 + FI→3 = 0
2 4 2 4
(III)
(I) (II)
3 1 3 3 5
Figure 2. 8
FII→3 provient des déplacements des noeuds 3, 4 et 5

FI→3 provient des déplacements des noeuds 3,2 et 1
FIII→3 provient des déplacements des noeuds 3, 4 et 2
Ainsi : FII→3 = F33,II + F34,II + F35,II
Autrement dit, l'action de l'élément II sur le noeud 3 est due :
- aux déplacements du noeud 3 : → k33,II

et : FI→3 = F33,I + F32,I + F31,I
Autrement dit, l'action de l'élément I sur le noeud 3 est due :
- aux déplacements du noeud 3 : → k33,I

enfin : FIII→3 = F33,III + F34,III + F32,III
Autrement dit, l'action de l'élément III sur le noeud 3 est due :
- aux déplacements du noeud 3 : → k33,III

RDM 5.xx Page 45
Ces termes sont encadrés dans la matrice page suivante.
F
2 4 6
(III) (IV)
(I) (II)
1 3 5
Figure 2. 9
Il s'agit maintenant de constituer une matrice générale en "éclatant" les matrices de rigidité de
chaque élément en blocs (2,2) et en les plaçant convenablement dans la matrice générale.
On obtient ainsi la matrice associée à la numérotation choisie :
dép1 dép2 dép3 dép4 dép5 dép 6

u1 , v1 u2 , v2 u3 , v3 u4 , v4 u5 , v5 u6 , v6
force1 k11,I k12,I k13,I
X1 , Y1
k21,I k22,I k23,I k24,III
force2 + + +
X2 , Y2 k22,III k23,III k24,III
k31,I k32,I k33,I

force3 + +
X3 , Y3 k33,II k34,II k35,II
k32,III
+ +
k33,III k34,III
k42,II k43,IV k44,II k45,II k46,IV
force4 + + +
X4 , Y4 k43,III k44,III k45,IV
+
k43,IV
force5 k53,II k54,II k55,II k56,IV
X5 , Y5 + +
k54,IV k55,IV
force 6 k64,IV k65,IV k66,IV
X6 , Y6
Dans cette matrice, chaque sous-matrice [k] doit être précédée du signe moins car il s'agit des
actions des éléments sur le noeud isolé
On s'aperçoit que cette matrice s'organise autour de la diagonale principale (matrice bande) .
Ceci est dû à la numérotation choisie pour les noeuds.
RDM 5.xx Page 46
La demi-largeur de bande (1/2 LB) comptée en noeuds est égale à la plus grande différence
entre les indices des noeuds d'un même élément augmentée de un.
Le stockage de la matrice est alors celui de la bande.

On voit donc l'importance de la numérotation des noeuds. La plupart des logiciels comporte
une procédure de renumérotation des noeuds qui minimise la largeur de bande. Cette
procédure est totalement transparente pour l'utilisateur qui ne voit que sa propre numérotation
( si la numérotation a été faite manuellement).
On aboutit à un système linéaire de 12 équations à 12 inconnues, dont la matrice rigidité est

[Kg].
  u1  0 
   v  0 
   1  
  u2  0 
     
   v2  0 
 Kg  u3  0 
     
   v3  0 
   =  
   u4  0 
   v4  0 
     
   u5  0 
   v5  -F 
     
   u6  0 
   v6  0 
   
2.11.4 Introduction des conditions aux limites

Deux types de noeuds sont à considérer :
- les noeuds 3, 4, 5 et 6 sont à déplacements libres. Leurs déplacements sont

inconnus, par contre les efforts extérieurs appliqués en ces noeuds sont connus :
- nuls pour les noeuds 3, 4 et 5 ;
- (0, -F) pour le noeud 6
- les noeuds 1 et 2 sont encastrés sur le bâti, ainsi leurs déplacements sont nuls. Par
contre les efforts extérieurs ne sont pas connus.
Ainsi en un noeud, on connaît soit le déplacement, soit l'effort.
Le système linéaire précédent est donc restructuré pour ne garder dans ce système que les
noeuds dont les déplacements sont inconnus.
2.11.5 Résolution du système linéaire

Il s'effectue par la méthode de Gauss tant que les systèmes ne sont pas trop important.
Lorsque les systèmes sont très grands la méthode de Gauss peut conduire à des erreurs
importantes. On met alors en oeuvre des méthodes indirectes (itératives).
RDM 5.xx Page 47
Le résultat de cette résolution est l'obtention des déplacements des noeuds :
3, 4, 5 et 6
les déplacements des noeuds 1 et 2 sont nuls.
Ainsi, à ce stade, on connaît les déplacements de tous les noeuds du maillage.
2.11.6 Calcul des déformations et des contraintes

2.11.6.1 Vecteur déplacement en un point quelconque
Sur l'élément générique (l, m, n) les relations suivantes permettent de calculer les déplacements
en un point quelconque M de l'élément :

u

v
1

u
1
[C]
u (x,y)  1 1 x y 0 0 0 m
 =
 v (x,y)  2 ∆e
0 0 0 1 x
 y 
v

u
m

v
n
2.11.6.2 Tenseur des déformation en un point quelconque
Le tenseur des déformations se calcule par :
∂u ∂v 1 ∂u ∂v
εxx = , εyy = et εxy = ( + )
∂x ∂y 2 ∂y ∂x
uk 
v 
ε xx   y m -y n 0 y n -y k 0 yk-y m 0   k
  1  um 
Soit par : ε yy  =
 0 x n -x m 0 x k -x n 0 x m -x k   
 2ε  2 ∆e  x -x y m -y n x k -x n y n -y k x m -x k y k -y m   vm 
 xy  n m un 
 
 vn 
Le tenseur des déformations est constant sur l'élément.
2.11.6.3 Tenseur des contraintes en un point quelconque
Le tenseur des contraintes se calcule à partir du tenseur des déformations par la loi de
comportement :
≈ ≈
σ = [E] ε
RDM 5.xx Page 48
Le tenseur des contraintes est constant sur l'élément.
A partir de ces éléments fondamentaux, on peut calculer tous les éléments dérivés :
- critères de limite élastique ( Tresca, Von-Mises,......) ;
- les directions et les contraintes principales ;
- ....................
2.11.7 Organigramme du processus de résolution

Compte tenu de ce qui vient d'être dit, on peut résumer le processus à l'aide de l'organigramme
de principe suivant :
Entrée des données
Calcul des matrices élémentaires

dans le repère local
Assemblage de la matrice de rigidité globale
Introduction des conditions aux limites
Résolution du système d'équations

et
détermination des déplacements des noeuds
Calcul des effets élastiques

en fonction
des déplacements des noeuds
Edition des résultats
Figure 2. 10 : Algorithme de résolution
RDM 5.xx Page 49
2.12 Types d'éléments plus performants
On a vu que le triangle à 3 noeuds conduit à une approximation linéaire pour le champ des
déplacements et uniforme pour le tenseur des contraintes.
Des éléments plus performants sont utilisés lorsque l'on veut une approximation plus précise.
Eléments triangulaires
Linéaire : 3 noeuds Quadratique : 6 noeuds Cubique: 9 noeuds
Eléments quadrangulaires
Linéaire : 4 noeuds Quadratique incomplet : 8 noeuds
Quadratique complet : 9 noeuds Cubique : 12 noeuds
Figure 2.11 : Différents éléments utilisés en 2D
Exemple : Pour un domaine D à deux dimensions
Polynôme de degré 1 : 1, x, y → nd = 3
Polynôme de degré 2 : 1, x, y, x², y², xy → nd = 6
En d'autres termes :
- pour approximer linéairement u il faut disposer de 3 noeuds d'interpolation ;
RDM 5.xx Page 50
- pour approximer quadratiquement u il faut disposer de 6 noeuds d'interpolation ;
- pour approximer cubiquement u il faut disposer de 10 noeuds d'interpolation si l'on
veut disposer d'une base complète.
Dans certaines applications on se contente de base incomplète. Par exemple :
Base bilinéaire : 1, x, y, xy → nd = 4
3. Type d'éléments finis et indicateurs de choix

3.1 Problèmes pratiques posés à l'utilisateur de la
méthode des éléments finis
En pratique différents problèmes pratiques se posent à l'utilisateur de logiciel mettant en
oeuvre la méthode des éléments finis :
- choix des éléments du maillage: triangle, quadrangle, nombre de noeuds de ces

éléments...
- le choix des éléments étant fait, quelle doit être la densité de ce maillage ?
- validité de la solution approximative trouvée, la solution théorique n'étant
évidemment pas connue dans le cas général.
Soit un maillage M1 réalisé à l'aide de N éléments triangulaires à 6 noeuds.

Soit un autre maillage M2 réalisé à l'aide de 2N éléments triangulaires à 3 noeuds.
Bien qu'ayant le même nombre de noeuds, ces deux maillage ne sont pas équivalents. On a vu
que l'élément à 3 noeuds est à la base d'une approximation linéaire des déplacements, alors
que celui à 6 noeuds conduit à une approximation quadratique des déplacements et est donc
beaucoup plus performant.
Il est difficile à l'heure actuelle de répondre de façon très précise à ces différents points. Dans
la référence (11) les auteurs précisent d'ailleurs que :
"beaucoup d'études ont été consacrés aux aspects théoriques mais peu d'articles traitent de
la démarche à suivre pour analyser et concevoir une pièce mécanique dans le contexte
industriel".
De façon générale on peut cependant donner quelques recommandations de bon sens.

L'approximation de la solution d'un problème d'élasticité linéaire doit résulter de plusieurs
études :
- considérer plusieurs maillages avec le même type d'éléments ;

- utiliser des éléments de type différents.
Il semble que l'expérience joue un grand rôle dans la réponse aux questions posées. Les
auteurs précisent que les systèmes experts associés à des bases de connaissances devraient
jouer un rôle important dans les années à venir.
3.1.1 Choix des éléments
RDM 5.xx Page 51
- pour les problèmes plans ou axisymétriques où la zone à discrétiser ne comporte pas
d'accidents, les éléments triangulaires à 3 noeuds ou quadrangulaires à 4 noeuds sont
économiques ;
- si la flexion intervient, les éléments de type quadrangle à 8 noeuds, 9 noeuds ou

triangulaire à 6 noeuds sont performants.
3.1.2 Influence du maillage. Etude sur un cas test

Poutre encastrée :
Acier : E=210000 MPa , ν = 0,28

Epaisseur e = 10 mm
Longueur L = 90 mm
Hauteur H = 30 mm
Charge répartie p = 100 N/mm
On fait varier deux types de paramètres du maillage :
- le type d'éléments :
- triangles à 3 noeuds ;
- quadrangles à 4 noeuds ;
- triangles à 6 noeuds ;
- quadrangles à 8 noeuds ;
- la densité des éléments.
p = 100 N/mm
Acier 30
e = 10 mm
90
x
Figure 3. 1
Différents maillages (suivant x, suivant y) :
triangle à Quadrangle à Triangle à Quadrangle à Quadrangle à

3 noeuds 4 noeuds 6 noeuds 8 noeuds 9 noeuds
RDM 5.xx Page 52
a) maillage 3 x 1; g) maillage 3 x 1; m) maillage 3 x 1; p) maillage 3 x 1; s) maillage 3 x 1;
b) maillage 6 x 2; h) maillage 6 x 2; n) maillage 6 x 2; q) maillage 6 x 2; t) maillage 6 x 2;
c) maillage 9 x 3; i) maillage 9 x 3; o) maillage 9 x 3; r) maillage 9 x 3; u) maillage 9 x 3;
d) maillage 12 x 4; j) maillage 12 x
e) maillage 15 x 5; 4;
f) maillage 18 x 6; k) maillage 15 x 5;
l) maillage 18 x 6;
Nature des éléments Maillage Nb Nb Valeur absolue du

noeuds éléments déplacement
suivant y en mm
Triangle à 3 noeuds a 8 6 0,063
b 21 24 0,110
c 40 54 0,142
d 65 96 0,160
e 96 150 0,170
f 133 216 0,176
Quadrangle à 4 noeuds g 8 3 0,133
h 21 12 0,171
i 40 27 0,182
j 65 48 0,186
k 96 75 0,188
l 133 108 0,190
Triangle à 6 noeuds m 21 6 0,185
n 65 24 0,192
o 133 54 0,193
Quadrangle à 8 noeuds p 18 3 0,186
q 53 12 0,192
r 106 27 0,193
Quadrangle à 9 noeuds s 21 3 0,189
t 65 12 0,193
u 133 27 0,1937
La figure suivante montre que le triangle à 3 noeuds est très peu performant. Même lorsque le
nombre de noeuds augmente, la solution ne tend que très lentement vers la solution exacte.
Ainsi, il ne sert à rien de se servir d'élément triangulaire à 3 noeuds en élevant

inconsidérément le nombre de noeuds.
L'élément quadrangle à 4 noeuds est légèrement plus performant que le triangle à 3 noeuds.
Par contre, triangle à 6 noeuds et quadrangle à 8 et 9 noeuds sont des éléments nettement plus
performants. La solution est rapidement approchée, même avec un faible nombre de noeuds.
On peut donc conclure de cette petite étude que :
- le choix du type d'élément est très important ;
RDM 5.xx Page 53
- en flexion, triangle à 6 noeuds et quadrangle à 8 et 9 noeuds sont des éléments
performants
Le test a porté sur la convergence des déplacements. En construction la convergence des

contraintes doit aussi être abordée. Sur l’exemple traité le maillage (s) de 3 éléments à 9
noeuds donne un résultat satisfaisant pour le déplacement. Il ne satisfait pas les contraintes. Il
sera donc nécessaire d’augmenter le nombre d’éléments pour satisfaire les deux conditions :
deplacements et contraintes
déplacement y
0,197
0,063 : triangle à 3 noeuds

: quadrangle à 4 noeuds
: triangle à 6 noeuds
nb de noeuds
40 100 150
Figure 3. 2
La valeur de référence en prenant en compte l’effort tranchant se calcule par la relation :

[Timoshenko : Résistance des matériaux Tome 1]
l4 l2
y=q + qα et vaut y = 0.197 mm
8 E IZ 2 AG
avec :
RDM 5.xx Page 54
3 4
α = section rectangulaire ;α = section circulaire
2 3
A = aire de la section droite
Iz = moment quadratique suivant z
E = module de YOUNG
E
G=
2 (1+ν )
ν = coefficient de Poisson
RDM 5.xx Page 55
4. Description des possibilités du logiciel
Etude de structures planes ou de révolution et calcul de caractéristiques de section droites
par la méthode des éléments finis en tenant compte des hypothèses suivantes :
- matériau homogène et isotrope
- comportement linéaire et élastique
- petits déplacements
Ecrans de RDM Version 5

Au lancement du programme M. E. F. l’écran ci dessous apparaît.
Il est composé :
- d’une barre titre RDM - Eléments finis
- d’une barre de menus des différents modules accessibles
- d’une zone graphique dans laquelle est indiqué :
• la version du logiciel et la date de mise à jour
• la signature de l’utilisateur
Figure 4.1
Après le choix du module, l’écran du programme apparaît (voir figures 4.2 et 4.3).
Il est composé :
- de la barre titre du module indiquant le type de problème suivi du nom de l’étude
traitée Elasticité : contraintes planes [ chemin\nom de l’étude ]
- de la barre de menus déroulants,
- de la zone graphique
- d’une zone de 2 lignes qui renseignent sur les actions en cours
RDM 5.xx Page 56
La saisie de données est effectuée dans des fenêtres flottantes qui apparaissent en fonction des
actions demandées.
Figure 4.2
Figure 4.3
RDM 5.xx Page 57
Description des possibilités offertes par le Module Eléménts Finis (M. E. F.) du logiciel.
Cette description est faite à partir de l’exemple suivant (figure 4.4):
Eprouvette rectangulaire d’épaisseur 5mm comportant une singularité de forme, encastrée à une
extrémité et
- soumise à un effort de 6000 N suivant son axe principal (sollicitation de traction).
- soumise à un couple de 250 Nm (sollicitation de flexion pure).
Zone d’épaisseur N°1 Zone d’épaisseur N°2
60 N/mm
Cas de charge numéro1
100 N/mm
Cas de charge numéro 2
-100 N/mm
Figure 4.4
La démarche d’étude d’une pièce mécanique passe par trois étapes :
1 - Modélisation géométrique et mécanique.
2 - Calcul.
3 - Exploitation des résultats.
Il est souvent nécessaire d’effectuer un retour à la modélisation pour modification après analyse
des résultats.
4.1 Modélisation
RDM 5.xx Page 58
4.1.1 Modélisation de la géométrie
La géométrie doit être aussi proche que possible de la géométrie réelle de la pièce à étudier. Ce
modèle peut être réalisé avec le module Dessin/Maillage du logiciel. Pour des géométries
complexes, il est souvent plus aisé d’utiliser un logiciel de D. A. O. permettant de la
sauvegarder dans un fichier au format standard I. G. E. S. ou D. X. F1.
Cette géométrie doit définir les contours extérieurs de la pièce (traits continus) ainsi que les
contours intérieurs (traits interrompus). Les contours intérieurs définissent des frontières
entre des zones de caractéristiques différentes (épaisseurs, matériaux, etc...). Un enlèvement de
matière (Zone sans matière type perçage débouchant) sera limité par un contour extérieur.
Le fichier au standard (I. G. E. S. ou D. X. F. ) est repris avec le menu Dessin/Maillage du
module M. E. F. .
Les écrans figure 4.5 et figure 4.6 apparaissent successivement :
Pointer sur Dessin/Maillage (figure 4.5)

Pointer sur Fichier (figure 4.6)
Glisser sur Importer
Glisser sur IGES
Renseigner le champ : Nom du fichier
Figure 4.5
1 I. G. E. S. est une norme d’échange graphique de données (Initial Graphics Exchange Spécification). Les entités
récupérées par RDM sont : le segment, le cercle, l’arc de cercle, (le point).
D. X. F. est un standard d’échange graphique de Autodesk Co.
RDM 5.xx Page 59
Figure 4.6
Figure 4.7
Après avoir précisé les unités employées, le logiciel récupère la géométrie, la traite et l’affiche
à l’écran. Si elle est incomplète, elle peut être modifiée. C’est dans cette phase que l’on
introduit les points libres qui permettent de définir un noeud imposé pour le maillage (points
particuliers qui peuvent servir pour positionner une charge, un appui, localement une
densification locale de maillage).
RDM 5.xx Page 60
Figure 4.8
Ce transfert étant terminé, on devra modéliser, si ce n’est déjà fait, les limites de zones
matérialisées par des traits interrompus et, si nécessaire, les points à mailler.
Les renseignements concernant un segment de la

géométrie sont obtenus en pointant l’élément et
en cliquant avec le bouton droit de la souris. Le
résultat apparaît dans une fenêtre identique à
celle ci-contre. Le segment pointé est le segment
bord gauche de l’éprouvette.
Les coordonnées d’un point sont obtenues de la
même manière.
Pour créer les points de contrainte de maillage il faut tout d’abord modéliser des points de
géométrie, et ensuite définir ces points comme points à mailler (noeuds caractéristiques de
maillage).
Définition du point géométrique (figure 4.9)
Pointer le menu déroulant Modéliser
Pointer sur Point
Sélectionner par pointé l’icône permettant de construire le point
Exécuter la commande. Une petite croix apparaît.
RDM 5.xx Page 61
Définition de ce point comme point a mailler (figures 4.10 et 4.11)
Pointer le menu déroulant Modéliser
Pointer sur Points à mailler
Pointer sur le point géométrique. La croix précédente est encapsulée.
Figure 4.9
Remarques :
Un point à mailler ne peut être créé que sur des points géométriques existants.
Les points à mailler peuvent être construits sur les contours, après construction du point
géométrique. Ils découpent la ligne du contour sur laquelle ils sont placés. Cela présente un
intérêt pour le traitement de problèmes avec maillage non uniforme de la frontière
(densification du maillage sur des lignes déterminées). Dans ce cas, la densification
automatique n’est pas (actuellement) possible autour de ces points. Il est évident que les points
extrémités de segments, peuvent êtres désignés directement (sans passer par le menu point)
comme points à mailler.
Par exemple sur la figure 4.10, six points sont définis « Points à mailler » :
- Le point milieu du segment vertical situé à droite (point servant à définir la zone d’application
des charges), construit comme point milieu d’un segment. Ce point découpe la ligne verticale
en deux segments (lignes 2 et 3).
- Les 5 points situés au droit de le singularité de forme (zone de densification du maillage),
construit comme points situés entre 2 points pour les 3 points isolés et directement pointés sur
RDM 5.xx Page 62
la géométrie pour les 2 points sur la frontière (ces points proviennent de la construction DAO
comme frontière de 2 arcs de cercle).
Noeud 2 Ligne 2
Ligne 1
Noeud 3
Noeud 4
Noeud 5 Noeud 1
Noeud 6
Ligne 3
Figure 4.10
Figure 4.11
Il est alors possible de passer au maillage en pointant sur le menu Mailler (Delaunay). Cette
action s’accompagne de la création d’un fichier temporaire de description de la géométrie
enregistré avec le nom $$$.GEO.
Pour enregistrer la géométrie avec un nom de fichier spécifique :
Pointer le menu déroulant Fichier

Pointer sur enregistrer
Renseigner le champ nom du fichier
RDM 5.xx Page 63
4.1.2 Modélisation du maillage
Il existe deux possibilités pour réaliser le maillage. Le choix proposé par défaut est l’utilisation
d’un mailleur automatique (Delaunay). Il est possible, en modifiant les options d’activer un
module de maillage par blocs offrant un plus grand choix d’éléments, mais d’un maniement
plus complexe.
Le maillage de la pièce sera réalisé à l’aide du mailleur automatique.
Pointer sur Mailler (Delaunay)
L’écran figure 4.12 est affiché. Sur le contour apparaît un découpage uniforme matérialisé par
des points. Ce découpage tient compte des paramètres par défaut qui sont sauvegardés dans le
fichier de configuration RDMDAT.INI après chaque modifications.
Figure 4.12
4.1.2.1 Définition des paramètres de maillage
L’accès à ces paramètres s’obtient en utilisant le menu Modéliser

Pointer sur Modéliser
Glisser sur Paramètres
Cette action est suivie de l’ouverture de la fenêtre Paramètres (voir figure 4.13)
La ligne Discrétisation uniforme de la frontière est cochée. C’est cette option qui
prédétermine le découpage uniforme de la frontière en n segments. La valeur de n, calculée
RDM 5.xx Page 64
automatiquement dépend directement de la valeur du paramètre Nombre d’éléments
souhaités.
Le Nombre d’éléments, triangles, ou quadrangles pour l’option maillage par blocs, peut
varier de 10 à 1700. Le nombre de noeuds autorisés est de 3000. Pour ce problème, 600
éléments sont souhaités. La valeur de ce champ sera modifiée en fonction des résultats obtenus.
Dilution (0.5..2) : ce champ définit une valeur qui permet de modifier en diminuant (< 1) ou
en augmentant (>1) de manière significative la taille des éléments lorsque l’on s’éloigne de la
frontière. En augmentant la valeur proposée par défaut (1.3) on va diminuer le nombre
d’éléments générés.
Amélioration : définit le nombre d’itérations de calculs lors de la phase d’amélioration du
maillage. Cette opération consiste à déplacer chaque noeud qui n’est pas situé sur la frontière
au centre de gravité des noeuds qui lui sont connectés. Répétée plusieurs fois elle
« uniformise » le maillage.
Qualité du Jacobien : valeur minimale de la qualité géométrique de l’élément créé selon le
critère du Jacobien. La valeur idéale qui peut être obtenue est 1 (triangle à bord droit). La
valeur par défaut 0.7 est une valeur très acceptable. Dans la pratique, il n’est pas utile de
modifier ce champ. Il est à noter que la qualité du Jacobien n’intervient que lorsqu’il existe des
éléments à bords curvilignes. Elle est donc liée au choix de l’élément de bibliothèque qui
servira au maillage (par exemple élément quadratique triangle à six noeuds à bords curvilignes).
Figure 4.13
Dans le cas où l’option Discrétisation uniforme de la frontière n’est pas cochée, chaque
ligne de la géométrie est discrétisée avec un élément. Il conviendra alors de définir
manuellement le nombre d’éléments à générer pour chaque ligne du domaine (contours
RDM 5.xx Page 65
extérieurs et intérieurs) en utilisant le sous menu Discrétiser la frontière du menu
Modéliser et pointer sur ces éléments comme le montre la figure 4.14.
La fonction est modale. Si plusieurs lignes doivent avoir la même discrétisation, il suffit de les
pointer les unes après les autres.
Pour modifier la valeur de discrétisation sur un élément frontière, on rappelle la fenêtre de
saisie de la valeur en appuyant sur le bouton droit de la souris.
Cette action, sur le bouton droit, quitte la modalité et rappelle la fonction en cours.
6 éléments 8 éléments 3 éléments

5 elements
11 éléments
Figure 4.14
4.1.2.2 Choix de l’élément
Après modification éventuelle des paramètres, discrétiser le domaine en sélectionnant l’élément

adéquat par pointé de cet élément dans la fenêtre d’icônes Mailler.
Si cette fenêtre n’apparaît pas à l’écran, on l’obtient en sélectionnant les menus :
RDM 5.xx Page 66
Modéliser
Discrétiser le domaine
La figure 4.15 montre un maillage avec les paramètres précédents pour un élément à 6 noeuds à
bords curvilignes en discrétisation uniforme de la frontière.
Figure 4.15
Un tel maillage peut ne pas paraître satisfaisant, notamment au droit des variations de section. Il
est alors possible de densifier localement ce maillage autour d’un noeud isolé en le
sélectionnant après avoir modifié la valeur du coefficient local (voir figure 4.16).
Ce paramètre local est modifié avec les menus :

Modéliser
Densifier
RDM 5.xx Page 67
Figure 4.16
La valeur par défaut est 1.
Cette valeur peut être modifiée et n’affectera que les noeuds isolés c’est à dire les noeuds qui
n’appartiennent pas à une frontière.
Il est donc important de prévoir cette éventualité de densification locale lors de la création de la
géométrie. Toutefois un retour à la géométrie est possible par le menu Dessin pour pallier un
oubli.
Les figures 4.17 à 4.21 montrent le maillage sur la même géométrie avec des paramètres différents :
Discrétisation uniforme de la frontière

Nombre d’éléments : 600
→ Dilution : 1.3
Amélioration : 3
→ Qualité du jacobien 0.9
→ Densification : 2 (autour du premier point isolé vertical. Noeud à
mailler 3)
Noeud à mailler 3
1173 Noeuds, 544 triangles, 2 sous domaines
Figure 4.17

RDM 5.xx Page 68
→ Dilution : 1.1
Amélioration : 3
→ Qualité du jacobien 0.7
→ Densification : 1
Figure 4.18

→ Dilution : 2
Amélioration : 3
Qualité du jacobien 0.7
Densification : 1
Figure 4.19
RDM 5.xx Page 69
Dilution : 2
Amélioration : 3
→ Densification : 3 (autour du deuxième point isolé vertical point à
mailler 4)
Noeud à mailler 4
Figure 4.20
Discrétisation non uniforme de la frontière (conforme à la figure 4.14)

Dilution : 2
Amélioration : 3
Densification : 1
Figure 4.21
Le maillage qui « à priori » donnera les meilleurs résultats est celui de la figure 4.21. La
discrétisation au droit de la zone perturbée étant la plus fine. Il ne présente pas d’écarts de
surface trop importants entre deux mailles adjacentes. Il pourrait néanmoins être amélioré en
augmentant le nombre d’éléments.
RDM 5.xx Page 70
4.1.2.3 Vérification de la qualité du maillage
Menu Afficher (figure 4.22)
Figure 4.22
Le maillage peut être visualisé de différentes manières :
Eléments séparés
Eléments colorés
Eléments séparés et colorés
La qualité géométrique peut être vérifiée à l’aide de 3 critères :

Distorsion : qualité idéale = 1. Cela correspond à un triangle équilatéral ou un carré parfait.
Noeud milieu : qualité idéale = 0. Utilisé pour vérifier la qualité d’éléments à noeud milieu à
bords curvilignes (exemple : triangle 6 noeuds).
Jacobien : Qualité idéale = 1. Utilisé pour vérifier la qualité d’éléments à noeud milieu à
bords curvilignes (exemple : triangle à 6 noeuds).
Pour une pièce qui ne possède pas de formes modélisées par des arcs de cercles ou cercles, il
est inutile de visualiser les critères Noeud milieu et Jacobien. On pourra remarquer que ces
deux critères sont très voisins, et donnent pratiquement le même indicateur.
Le critère Distorsion est celui qui permettra de visualiser de la manière la plus simple la
qualité.
En règle générale : [Référence : « Modélisation des ouvrages » HERMES éditeur.]
RDM 5.xx Page 71
Un maillage non uniforme sera correct si le rapport des surfaces de deux éléments adjacents est
inférieur à 2 (ou 0.5).
Les éléments « triangle » doivent ressembler à des triangles équilatéraux.
Les éléments ne doivent pas être écrasés. La plus petite dimension, la plus petite hauteur pour
un triangle, ne doit pas dépasser trois à cinq fois la plus grande dimension.
Un bon maillage résulte de la bonne association élément, sollicitation.
Qualité du maillage selon le critère Distorsion (figure 4.23)
Figure 4.23
Bien que la qualité du maillage apparaisse correcte dans ce problème, on peut remarquer qu’au
droit de la variation de section, la densité de maillage n’est pas très importante. Il est
souhaitable dans ce cas de prévoir d’autres modèles du maillage afin de voir l’incidence de
cette discrétisation dans la convergence des résultats. Si cela est possible, cette densification
n’entraînera pas d’augmentation importante du nombre des éléments. En effet, toute
augmentation du nombre d’éléments est pénalisante pour le temps de calcul.
A titre de comparaison les calculs vont être menés pour 2 modèles de maillage qui
correspondent aux figures 4.17 et 4.21.
Qualité du maillage selon le critère du Jacobien (figure 4.24)
RDM 5.xx Page 72
Figure 4.24
Qualité du maillage selon le critère du Noeud milieu (figure 4.25)
Figure 4.25
RDM 5.xx Page 73
La qualité géométrique d’un élément selon les trois critères peut être visualisée à l’écran en
pointant l’élément avec le bouton droit de la souris voir figure 4.26.
Seule la valeur minimale du jacobien est conservée et sert pour le critère de qualité.
Figure 4.26
Si les valeurs constatée pour un élément ne sont pas satisfaisantes, le seul moyen d’agir pour
tenter de les modifier est de changer les paramètres :
en augmentant le nombre d’éléments
en augmentant la valeur du Jacobien s’il s’agit d’un élément à bords curvilignes
en plaçant sur la frontière (ou ailleurs)avec des noeuds à mailler
en procédant éventuellement à un maillage non uniforme de la frontière.
Il faut, ensuite, relancer l’opération de maillage.

Il n’est pas possible d’agir directement sur le modèle maillé en déplacent « manuellement » un
noeud.
Il est possible, par ce menu, d’afficher d’autres informations :

Numéros
Un noeud
Noeuds
Un élément
Eléments
4.1.2.4 Sauvegarde du maillage
Toutes ces opérations étant terminées, la sauvegarde peut être effectuée.

Attention
La commande enregistrer n’effectue que la sauvegarde de la géométrie sans le maillage
(extension .GEO) mais, avec les derniers paramètres utilisés (densification, qualité Jacobien,
uniforme, ....).
Le maillage sera sauvegardé dans un fichier possédant l’extension .CAL après avoir
sélectionné le type de problème à traiter.
Dans le cas du problème traité le choix est Elasticité/thermique.
RDM 5.xx Page 74
A ce stade le logiciel propose de sauvegarder. Cette proposition de sauvegarde ne concerne que
la géométrie (extension de fichier .GEO) contour de la pièce et paramètres du maillage.
Accepter cette sauvegarde si cela n’a pas déjà été fait. Si cette sauvegarde n’est pas acceptée le
logiciel en crée une qui a pour nom $$$.GEO.
Le module Elasticité/thermique étant lancé, choisir le type de problème à traiter.
Dans le cas du problème traité notre choix se porte sur Contraintes planes.
Figure 4.27
Ce choix étant effectué il est alors possible de sauvegarder le maillage qui servira de base à la
résolution du problème. le fichier à pour extension .CAL.
Menus :
Fichier
Enregistrer
Figure 4.28
RDM 5.xx Page 75
4.1.2.5 Maillage par blocs
Le maillage par blocs correspond à une partition prédécoupée en « macro-éléments » destinés à

être subdivisés en éléments finis. Ce découpage correspond donc à un prémaillage. La forme de
ces éléments doit être simple (triangles ou quadrangles).
La figure 4.29 présente un découpage par blocs pour le problème traité. Ce découpage sera
réalisé à partir d’un logiciel de D. A. O. et importé dans le logiciel en utilisant un format
d’échange comme expliqué précédemment (voir modélisation de la géométrie).
Figure 4.29
L’accès au maillage par blocs en sélectionnant le menu déroulant Outils
Pointer sur Outils
Glisser sur Options
Cette action est suivie de l’ouverture de la fenêtre Options
Dans cette fenêtre cocher par pointé sur la zone pour rendre le mailleur par blocs actif.
Figure 4.30
Le menu Mailler par blocs apparaît et devient actif dans la barre de menu.
RDM 5.xx Page 76
Figure 4.31
Tous les blocs sont discrétisés à un élément. Pour réaliser le maillage il convient de discrétiser
les limites de zone.
Ponter sur Mailler par blocs
Pointer sur Mailler
Glisser sur Limites de zone
Figure 4.32
RDM 5.xx Page 77
Dans ce type de maillage le nombre d’éléments affecte les lignes opposées des différents blocs.
Pour des blocs de type quadrangle, deux lignes opposées possèdent le même nombre
d’éléments. Toutes les lignes d’un bloc de type triangle possèdent le même nombre d’éléments.
La figure 4.33 représente une discrétisation d’un élément quadrangle. Les éléments adjacents
sont affectés par la discrétisation.
Figure 4.33
La figure 4.34 représente une discrétisation d’un élément triangle. Les éléments adjacents
(quadrangles) sont affectés par la discrétisation.
Figure 4.34
La figure 4.35 représente une discrétisation de l’ensemble des zones.
Figure 4.35
RDM 5.xx Page 78
Le maillage par blocs offre une palette plus grande quand au choix des éléments de maillage.
Des éléments quadrangles à quatre, huit ou neuf noeuds sont disponibles.
Figure 4.36
Le modèle suivant est maillé avec un élément quadrangle à huit noeuds.

Il possède 1885 noeuds et 592 éléments
Figure 4.37
La sauvegarde est effectuée suivant la même démarche que précédemment.
Remarque :
Ce type de maillage est souvent difficile à réaliser. Le prémaillage étant une étape délicate. La
qualité des maillage obtenus est généralement moins bonne qu’avec le mailleur automatique.
Cependant dans certains cas précis où les éléments quadrangles sont utilisés il est nécessaire
d’y avoir recours.
4.2 Sauvegardes des études.

Une étude sera définie à l’aide de deux fichiers.
Un fichier possédant l’extension .GEO, contenant la définition géométrique de la structure

ainsi que les paramètres du maillage en cours lors de la dernière sauvegarde.
Un fichier possédant l’extension .CAL, contenant la définition géométrique du modèle maillé

ainsi que l’ensemble de toutes les données mécaniques de l’étude.
Ces données mécaniques sont :
RDM 5.xx Page 79
Le type du problème
Les matériaux
Les liaisons avec le milieu extérieur
Les cas de charges
Etc..
C’est le fichier avec l’extension .CAL qui servira aux différents calculs.
Une étude étant créée, fichier .CAL, il est IMPOSSIBLE de modifier la géométrie du maillage
sur ce fichier. Pour une autre géométrie du maillage, il faut faire appel au fichier .GEO et
générer un nouveau maillage et, par la suite une nouvelle étude.
Pour de la géométrie : EPROUVT.GEO
l’exemple qui illustre ce chapitre, les fichier ont pour nom :
Définition des études :

figure 4.17 : EPROUVT.CAL
figure 4.21 : EPMAN1.CAL
Il est ainsi possible de créer autant d’études que l’on souhaite.
4.2.1 Modélisation mécanique

4.2.1.1 Définition des épaisseurs
La définition des épaisseurs est obtenue à partir des menus :

Modéliser
Epaisseurs
Une fenêtre permet d’affecter des couleurs à des zones d’épaisseur différente, préalablement
définies dans la modélisation géométrique. Ces zones d’épaisseur différentes sont ici appelées
« groupe ».
Pour attribuer une couleur :

Pointer sur la palette la couleur choisie. Cette couleur devient couleur courante.
Attribuer la couleur à un Elément, un Groupe, tout le Domaine en pointant

un élément de la géométrie après avoir sélectionné la commande Elément,
Groupe, Domaine.
RDM 5.xx Page 80
Figure 4.38
Figure 4.39
La figure 4.39 montre le résultat pour les deux zones de l’éprouvette.
RDM 5.xx Page 81
Pour définir l’épaisseur de la zone :
Pointer sur le bouton Définir
Pointer un élément du groupe dans la fenêtre graphique, la fenêtre suivante est
ouverte, renseigner le champ Epaisseur.
Figure 4.40
Il est possible de traiter des études comportant au maximum 15 épaisseurs différentes.
Dans notre problème les deux zones auront la même épaisseur de 5 mm.
4.2.1.2 Définition des matériaux
La démarche est identique à la définition des épaisseurs.

La définition des matériaux est obtenue à partir des menus :
Modéliser
Matériaux
Une fenêtre permet d’affecter des couleurs à des zones de matériaux différents, préalablement
définies dans la modélisation géométrique. Ces zones de matériaux différents sont ici appelées
« groupe ».
Pour attribuer une couleur :
Pointer sur la palette la couleur choisie. Cette couleur devient couleur courante.
Figure 4.41
RDM 5.xx Page 82
Pour définir le matériau de la zone, pointer sur le bouton Définir
Pointer un élément du groupe dans la fenêtre graphique, la fenêtre figure 4.42 est
ouverte, renseigner les différents champs caractéristiques du matériau.
Les études traitées peuvent comporter 15 matériaux différents.
Figure 4.42
Les zones matériaux sont indépendantes des zones épaisseurs. Le découpage géométrique doit
tenir compte de ce paramètre.
Le bouton Bibliothèque permet de sélectionner un matériau défini dans une base de données
contenue dans le fichier MATERIAU.DAT
Le choix du matériau pour l’éprouvette est un alliage d’aluminium A-U4G
Figure 4.43
RDM 5.xx Page 83
La vérification des données Epaisseur et Matériau pour une zone ou un élément donné est
obtenue en pointant un élément de la zone (ou l’élément considéré) avec le bouton droit de la
souris. Les renseignements apparaissent dans une fenêtre indiquant les caractéristiques
mécaniques ainsi que la numérotation des noeuds.
Figure 4.44
4.2.1.3 Définition des liaisons
La définition des liaisons est obtenue à partir des menus :

Modéliser
Liaisons/Symétrie
Une fenêtre d’icônes indiquant les différents types de liaisons disponibles est ouverte.
Symétrie par un plan d’abscisse x = A

Symétrie par un plan d’ordonnée y = B
Liaisons réparties sur une frontière
Liaisons nodales
Appui élastique : Fy = -K dy
Appui élastique : Fx = -K dx
Déplacement imposé sur y
é Déplacement imposé sur x
Liaison rotule
Appui simple incliné

Déplacement sur y = 0
Déplacement sur x = 0
Figure 4.45
RDM 5.xx Page 84
Le choix NOD (liaisons NODales) ou REP (liaisons REParties sur une ligne frontière) affecte
la barre de titre de la fenêtre d’icônes.
La liaison encastrement est obtenue par un ensemble de liaisons rotule (dx = 0 ; dy = 0)

réparties sur tous les noeuds du bord gauche de l’éprouvette (ligne 1).
Dans la fenêtre d’icônes, pointer sur :

REP
Liaison rotule
Sélectionner la ligne 1
Les liaisons sont représentées dans la fenêtre graphique. La vérification est obtenue en utilisant
le bouton droit de la souris par pointé sur un noeud lié.
Figure 4.46
4.2.1.4 Définition des cas de charges
La définition des cas de charge est obtenue à partir des menus :

Modéliser
Cas de charges
RDM 5.xx Page 85
Une fenêtre d’icônes indiquant les différents chargements disponibles est ouverte.
Il est possible de créer 15 cas de charges différents de 300 charges chacun pour un problème
donné.
Gradient thermique
Pesanteur
Pression
Charge linéique
Charge nodale
Figure 4.47
Le chargement numéro 1 (éprouvette en traction) est constitué d’une charge répartie de 60

N/mm suivant l’axe x.
Sélectionner l’icône Charge linéique
Renseigner la boite de dialogue.
Figure 4.48
Sélectionner les 2 lignes frontière supportant ce chargement sur le bord droit de l’éprouvette
(lignes 2 et 3).
Le chargement est alors représenté dans la fenêtre graphique.
RDM 5.xx Page 86
Figure 4.49
Un deuxième cas de charge correspondant à un couple appliqué sur l’extrémité droite de

l’éprouvette est entré (éprouvette en flexion) en suivant la procédure suivante :
Modéliser
Ajouter un cas de charge
Sélectionner l’icône charge répartie

Fx = 100 N/mm ; Fy = 0 N/mm
Sélectionner la ligne 2 (bord droit de l’éprouvette situé au dessus
de l’axe de symétrie de direction x)
Action sur le bouton droit de la souris (arrêt de la modalité sur
la valeur de la charge). La fenêtre d’icônes revient en premier
plan.
Sélectionner l’icône charge répartie

Fx = -100 N/mm ; Fy = 0 N/mm
Sélectionner la ligne 3 (bord droit de l’éprouvette situé au dessous
de l’axe de symétrie de direction x)
Action sur le bouton droit de la souris (arrêt de la modalité sur
la valeur de la charge). La fenêtre d’icônes revient en premier
plan.
Action sur le bouton droit de la souris (arrêt de la modalité sur le
cas de charge).
RDM 5.xx Page 87
Le résultat apparaît dans la fenêtre graphique
Figure 4.50
L’utilisation du bouton droit de la souris permet de vérifier les valeurs de chargement en

sélectionnant la frontière (ou le noeud) sur laquelle est représenté ce chargement.
4.2.1.5 Définition du modèle d’étude dynamique
Cette définition est obtenue à partir des menus

Modéliser
Etude dynamique
Une fenêtre permettant d’indiquer les paramètres du calcul dynamique est ouverte.
Les paramètres à renseigner sont :
Nombre de fréquences (20 modes propres maximum)
Précision du calcul. Le logiciel arrête le calcul itératif sur la fréquence lorsque la
valeur indiquée est atteinte.
Calcul des énergies élémentaires.
Méthode de calcul souhaitée.
RDM 5.xx Page 88
Figure 4.51
La méthode Itération sur sous espace permet de calculer les n plus petites valeurs propres
d’un système de grande dimension.
La méthode Itération inverse est bien adaptée au calcul de la plus petite valeur propre.
La modélisation du problème étant terminée, la phase de calcul peut démarrer.

Avant de lancer le calcul, il est sage de sauvegarder les données du problème en utilisant les
menus :
Fichier
Enregistrer
Le fichier possède l’extension .CAL
Figure 4.52
RDM 5.xx Page 89
4.3 Calculs
Cette phase est automatique et ne demande aucune action particulière à l’opérateur.
4.3.1 Calcul statique
Le calcul est lancé en utilisant le menu :

Calculer
Analyse statique
La fenêtre de calcul affiche un certain nombre d’informations relatives au bon déroulement des
opérations. Le logiciel crée un fichier possédant l’extension .RES du déroulement de la phase
de calcul. Pour notre calcul, cela donne le fichier suivant :
+-------------------+
| Thermo-élasticité |
+-------------------+
Nom du projet : c:\rdmv5\livret\eprouvt
Date : 22 Décembre 1996

Début du calcul - Heure : 15/28/27
+------------+
| Assemblage |
+------------+
2312 inconnues
Mémorisation "Profil" : 138392 réels en 20 bloc(s) de 54 Ko
Table des actions de liaisons
Vecteur force global - Cas 1 2

Matrice de rigidité
Coefficient maximal = 3.25806496314347E+0006
Temps de calcul : = 0 H 0 min 17 s
+-----------------------------------------+
| Factorisation de la matrice de rigidité |
+-----------------------------------------+
pivot minimal = 1.11250918959535E+0003

ddl = 2312
noeud = 171
+------------------+
| Analyse statique |
+------------------+
Problème 1 : Cas de charges 1

Déplacements nodaux
Action(s) de liaison
Contraintes
Création du fichier : c:\rdmv5\livret\eprouvt.E1
Problème 2 : Cas de charges 2

Déplacements nodaux
Action(s) de liaison
Contraintes
Création du fichier : c:\rdmv5\livret\eprouvt.E2
RDM 5.xx Page 90
Fin du calcul
4.3.2 Calcul dynamique
Le calcul est lancé en utilisant le menu :

Calculer
Analyse dynamique
La fenêtre de calcul affiche un certain nombre d’informations relatives au bon déroulement des
opérations. Le logiciel crée un fichier possédant l’extension .RES du déroulement de la phase
de calcul. Pour notre calcul, avec un mode propre, cela donne le fichier suivant :
+-------------------+
| Thermo-élasticité |
+-------------------+
Nom du projet : c:\rdmv5\livret\eprouvt
Date : 27 Janvier 1997

Début du calcul - Heure : 11/27/32
+------------+
| Assemblage |
+------------+
2312 inconnues
Mémorisation "Profil" : 138392 réels en 20 bloc(s) de 54 Ko
Table des actions de liaisons
Matrice de rigidité
Coefficient maximal = 3.25806496314347E+0006
Matrice masse
Coefficient maximal = 6.31675624043176E-0004
+-----------------------------------------+
| Factorisation de la matrice de rigidité |
+-----------------------------------------+
pivot minimal = 1.11250918959535E+0003

ddl = 2312
noeud = 171
+-------------------+
| Analyse dynamique |
+-------------------+
Résolution du système [k][v]=lamda[m][v]
Itération sur sous-espace
Dimension du sous espace = 2
Itération 1
Résolution du système [k][v]=lamda[m][v] par la méthode de Jacobi généralisée

RDM 5.xx Page 91
i k[i,i] m[i,i] Valeur propre
1 2.45205990021416E-0008 2.72542481970146E-0012 8996.9824

2 2.78499036938262E-0011 6.05259356245540E-0017 460131.7337
Cycle 1
1 2.39505272473878E-0008 2.72006535359444E-0012 8805.1293

2 2.78061965299288E-0011 6.03479587046894E-0017 460764.4919
Cycle 2
1 2.39505272473878E-0008 2.72006535359444E-0012 8805.1293

2 2.78061965299288E-0011 6.03479587046894E-0017 460764.4919
Termes non diagonaux maximaux
k [ 1 , 2 ] = 0.00000000000000E+0000 Facteur de couplage =

0.00000000000000E+0000
m [ 1 , 2 ] = 0.00000000000000E+0000 Facteur de couplage =
0.00000000000000E+0000
Convergence en 2 cycles et 2 transformations
Mode Pulsation Fréquence
1 2967.3438 472.2674
Itération 2
1 5.45738569300535E-0008 6.35546468506303E-0012 8586.9184

2 1.24933425439289E-0009 3.40548615674818E-0015 366859.2961
Cycle 1
1 5.45765193681493E-0008 6.35577835534668E-0012 8586.9136

2 1.24808811768233E-0009 3.25335122515961E-0015 383631.5329
Cycle 2
1 5.45765193681493E-0008 6.35577835534668E-0012 8586.9136

2 1.24808811768233E-0009 3.25335122515961E-0015 383631.5329
k [ 1 , 2 ] = 0.00000000000000E+0000 Facteur de couplage =

0.00000000000000E+0000
m [ 1 , 2 ] = 0.00000000000000E+0000 Facteur de couplage =
0.00000000000000E+0000
1 2930.3436 466.3787
Erreur maximale = 2.54125970865131E-0002
Itération 3
RDM 5.xx Page 92
1 5.45343757877119E-0008 6.35111934622840E-0012 8586.5771

2 1.24361069275234E-0009 3.24433237618964E-0015 383317.9060
Cycle 1
1 5.45343757877119E-0008 6.35111934622840E-0012 8586.5771

2 1.24361069275234E-0009 3.24433237618964E-0015 383317.9060
k [ 1 , 2 ] = 4.32982845158451E-0015 Facteur de couplage = 5.25766962413037E-

0007
m [ 1 , 2 ] = 7.34758436783227E-0019 Facteur de couplage = 5.11866579000563E-
0006
1 2930.2862 466.3695
Erreur maximale = 3.91876099209395E-0005
Convergence en 3 itérations
Fin du calcul
Si le problème n’est pas complètement défini (absence de liaisons, de chargement, d’épaisseur)

un message apparaît, demandant de revoir la modélisation. Le calcul n’est pas lancé.
Un message de type « Pivot négatif ou nul » pendant le calcul, indique que le système possède
des mobilités. Les liaisons avec le milieu extérieur sont incomplètement définies ou
insuffisantes.
Durant cette phase, il est bon de contrôler la valeur du pivot minimal. Si cette valeur est
inférieure à 1, cela peut donner des résultats qui ne sont pas conformes. La méthode de
résolution utilisant les pivots de Gauss, les « erreurs » de calcul numérique se trouvent
amplifiées. Il faut éventuellement, dans ce cas, revoir le modèle.
Le calcul étant terminé, le menu Résultat de la barre de menus déroulants devient actif.
Passer à la phase d’exploitation des résultats en sélectionnant ce menu.
4.4 Exploitation des résultats
4.4.1 Résultats du calcul théorique

Calculs à partir des courbes du CETIM (voir courbes en annexe).
Pour les deux cas de charges, les abaques du CETIM [2] donnent les coefficients de
concentration des contraintes pour l’accident géométrique.
RDM 5.xx Page 93
Les coefficients de formes valent :
r 10 d 60
= = 0.5 = = 0.6
t 20 D 100
Cas de charge 1 : traction
Le coefficient de concentration des contraintes est Kt = 2.2

F
Contrainte nominale σ nom = .
S
Avec F = 6000 N et S = 60×5 = 300 mm2.
σ nom = 20 MPa
La contrainte réelle est estimée à :

σ réelle = K t σ nom
σ réelle = 44 MPa
Cas de charge 2 : flexion
Le coefficient de concentration des contraintes est Kt = 2.02

M fz
Contrainte nominale σ nom = × ymax .
I GZ
Avec Mfz = 250000 Nmm ; IGZ = 5×603/12 = 90000 mm4 ; ymax = 30 mm
σ nom = 83.4 MPa
La contrainte réelle est estimée à :

σ réelle = K t σ nom
σ réelle = 168.4 MPa
4.4.2 Résultats donnés par le logiciel
L’accès aux résultats du calcul s’obtient en utilisant le menu Résultat qui, maintenant, est
marqué en caractères gras dans la barre de menu.
Pointer sur Résultats
RDM 5.xx Page 94
Si l’étude comporte plusieurs cas de charges, ou plusieurs modes dans le cas d’une étude
dynamique, une fenêtre s’ouvre proposant de choisir parmi ces différents calculs.
Figure 4.53
Pointer sur le cas choisi
Le tracé de la structure, chargée et déformée, apparaît à l’écran en traits forts. En traits

interrompus, il subsiste le tracé de la structure au repos.
RDM 5.xx Page 95
Figure 4.54
L’accès à l’ensemble des résultats qui peuvent être consultés s’obtient en utilisant le nouveau
menu déroulant Résultats.
Pointer sur Résultats
Figure 4.55
RDM 5.xx Page 96
La fenêtre déroulante ouverte indique les différents types de résultats consultables.
Déformée : Visualisée par défaut après avoir sélectionné l’étude à traiter.

Vecteur déplacement : Visualise le champ de déplacement aux différents noeuds de la
structure
y
CHAMP DE DEPLACEMENT
x
Figure 4.56
En pointant avec le bouton droit de la souris sur un noeud de la structure représentée, on accède
aux valeurs numériques des déplacements et des contraintes pour ce noeud ainsi que les actions
de liaison si elles existent.
Figure 4.57
Faces principales : Visualise les directions principales aux différents noeuds de la

structure par des doubles flèches orientées indiquant s’il s’agit de
traction ou de compression.
RDM 5.xx Page 97
En pointant avec le bouton droit de la souris sur un noeud de la
structure représentée, on accède aux valeurs numériques des
déplacements et des contraintes pour ce noeud ainsi que les actions
de liaison si elles existent.
Figure 4.58
Cercles de Mohr : Visualise en un noeud de la structure l’état de contrainte dans le
plan de Mohr.
Pointer un noeud de la structure
NOEUD 12 y COORDONNEES DU NOEUD :

Y
x = 200.00 , y = 80.00 mm
TAU TENSEUR DES CONTRAINTES :
Xx
SIGMA xx = 40.0 MPa
SIGMA yy = 3.4 MPa
SIGMA xy = -0.5 MPa
SIGMA zz = 0.0 MPa
CONTRAINTES PRINCIPALES :
zY y Xx SIGMA
SIGMA XX = 40.0 MPa
SIGMA YY = 3.4 MPa
SIGMA ZZ = 0.0 MPa
ANGLE X / x = -0.85 DEGRES
SIGMA TRESCA = 40.0 MPa
SIGMA VON MISES = 38.4 MPa
Figure 4.59
Actions de liaison : visualise les actions de liaisons aux noeuds liés et permet d’obtenir
la valeur sur un noeud.
Pointer un noeud de la structure
RDM 5.xx Page 98
Figure 4.60
Isovaleurs : un deuxième sous menu apparaît permettant de choisir d’autres
types de résultats. Ils sont ici visualisés sous forme de domaines
colorés, 1 couleur = 1 intervalle, pour lesquelles les grandeurs sont
sensiblement identiques. Le choix de ce menu donne accès à la
sélection de 10 nouvelles grandeurs représentables sous cette forme.
RDM 5.xx Page 99
Figure 4.61
Dans la figure 4.62, le choix est porté sur Sigma xx afin de pouvoir comparer avec les
résultats théoriques précédents.
Figure 4.62
RDM 5.xx Page 100
La zone de contrainte maximale est située au droit de l’accident de forme, comme le montre le
zoom localisé sur cette forme. La valeur maximale est voisine de 40 MPa.
Figure 4.63
Cette valeur nous permet de calculer un coefficient de concentration des contraintes de Kt = 2
inférieur à la valeur théorique calculée précédemment.
Nous avons donc une « erreur relative » de 9 %.
Le sous menu Traction (Compression) permet de représenter la contrainte principale

maximales positive (négative).
Dans le cas du menu Traction, seules les contraintes principales maximales positives sont
représentées. Les valeurs négatives sont occultées et mises à zéro. Cette représentation peut être
intéressante dans le cas d’études de formes géométriques simples, qui peuvent s’apparenter à
des poutres, soumises, à de la traction.
Dans le cas du menu Compression, seules les contraintes principales maximales négatives
sont représentées. Les valeurs positives sont occultées et mises à zéro.
RDM 5.xx Page 101
Figure 4.64
Figure 4.65
RDM 5.xx Page 102
Représentation suivant le critère de Von Mises.
Figure 4.66
Isovaleurs sur Déformée : Donne les mêmes résultats que le menu Isovaleurs sur la
structure déformée
Isolignes : Donne les mêmes résultats que le menu Isovaleurs. Ils sont
ici visualisés sous forme de lignes, frontières des domaines
précédents colorés, pour lesquelles les grandeurs sont
identiques. Le choix de ce menu donne accès à la même
sélection des 10 grandeurs précédentes.
Coupe suivant une droite : Permet de spécifier un plan de coupe définissant une section
(ou portion de section) par une droite définie par 2 points.
La coupe est ici définie au droit de l’accident de forme au voisinage de la section droite
minimale. Elle est caractérisée par les 2 points A et A’.
Les 2 points sont choisis hors matière. Le point A est l’origine de l’axe de coupe.
RDM 5.xx Page 103
Figure 4.67
La fenêtre coupe est ouverte.

Dans cette fenêtre choisir le menu Grandeurs.
Les grandeurs représentables sont les mêmes que celles du sous menu Isovaleurs.
RDM 5.xx Page 104
Figure 4.68
Représentation de la contrainte Sigma xx le long de la section droite.
MPa
38.35
33.63
28.90
24.18
19.45
14.73 mm
36.89 49.02 61.15 73.28 85.41 97.54
SIGMA xx
MIN = 14.73 MPa A 67.42 mm
MAX = 38.35 MPa A 36.89 mm
Figure 4.69
Sur la représentation graphique (figure 4.69) :L’axe horizontal est défini par la droite A A’, A
étant l’origine de l’axe. Cet axe représente l’abscisse des différents points calculés de la
section.
L’axe vertical représente la grandeur choisie.
Trois lignes de commentaires renseignent sur les valeurs caractéristiques.
Remarques :
Sur l’axe horizontal le début de la section est marqué à 36.89 mm. Cette valeur représente la
distance du point A, premier point de la définition de la droite de coupe, avec le point
intersection de la droite et de la matière. La fin est située à 97.54 mm soit à 60.65 mm ce qui
correspond à la hauteur de section à cet endroit.
La contrainte maximale est de 38.35 MPa (la section n’est pas exactement la section minimale)
prés de A et 34 MPa prés de A’.
On voit ici l’incidence du maillage sur les résultats obtenus. Utilisation d’une densification
locale prés de A, pas de densification prés de A’.
Coupe suivant une ligne : permet de définir une ligne quelconque à partir d’une portion du
contour de la structure.
La ligne est définie par 3 points :
Un point début B, un point fin B’, un point situé entre les deux premiers ce qui permet de
localiser la portion du contour étudiée.
Dans la représentation qui suit, le choix est porté sur l’entaille.
RDM 5.xx Page 105
Figure 4.70
Le résultat demandé est toujours la contrainte Sigma xx. On voit apparaître au milieu la
contrainte maximale de 40 MPa affichée par le menu Isovaleurs.
MPa
39.99
31.51
23.03
14.56
6.08
-2.40 mm
0.00 10.06 20.13 30.19 40.26 50.32
SIGMA xx
MIN = -2.40 MPa A 37.50 mm
MAX = 39.99 MPa A 25.09 mm
Figure 4.71
Les résultats du Cas de charges numéro 2 sont obtenus à partir du menu Etudes
Pointer sur Etudes
RDM 5.xx Page 106
Choix du Cas de charges numéro 2
La structure déformée apparaît.
Figure 4.72
Afin de comparer les résultats sur la contrainte Sigma xx, celle-ci est visualisée en
Isovaleurs.
RDM 5.xx Page 107
Figure 4.73
La contrainte maximale est de 120.05 MPa ce qui conduit à un coefficient de concentration des
contraintes de Kt = 1.44. Cette valeur est inférieure à la valeur théorique calculée
précédemment.
Nous avons donc une « erreur relative » voisine de 28 %.
Figure 4.74
On remarque que pour les deux sollicitations, et pour un même maillage, le pourcentage
d’erreur diffère. Il est plus important pour la flexion que pour la traction. Ceci doit conduire à
la réalisation de nouvelles études qui donneront des résultats de meilleure qualité.
Les deux études qui suivent donnent des résultats qui convergent sans pour autant augmenter
notablement le nombre d’éléments.
Etude epman.cal : 783 Noeuds 352 éléments
Cas de charges 1
RDM 5.xx Page 108
Figure 4.75
contraintes de Kt = 2.34. Cette valeur est supérieure à la valeur théorique calculée
précédemment.
RDM 5.xx Page 109
Figure 4.76
Cas de charges 2
Figure 4.77
La représentation est en Isovaleurs sur déformée.
RDM 5.xx Page 110
précédemment.
Nous avons donc une « erreur relative » de 14 %.
Figure 4.78
RDM 5.xx Page 111
Etude epman1.cal : 915 Noeuds 412 éléments
Cas de charges 1
Figure 4.79
contraintes de Kt = 2.38. Cette valeur est supérieure à la valeur théorique calculée
précédemment.
Remarque.
Cette valeur est très voisine de celle obtenue avec l’étude epman.cal. Il y a convergence des
résultats vers une valeur qui ce situe au voisinage de 48 à 49 Mpa, valeur différente de la valeur
théorique calculée.
RDM 5.xx Page 112
Figure 4.80
Cas de charges 2
Figure 4.81
RDM 5.xx Page 113
La représentation est en Isovaleurs sur déformée.
précédemment.
Il apparaît pour ce cas de charges une convergence des résultats vers la valeur théorique
calculée.
Figure 4.82
On voit donc ici la nécessité de plusieurs calculs afin d’assurer une convergence des résultats.
La pratique permet de minimiser le nombre de ces calculs.
Les bons résultats ne sont pas forcément obtenus avec un très grand nombre d’éléments.
Le maillage doit être densifié dans les zones qui présentent les accidents de forme afin
d’augmenter la précision du calcul.
RDM 5.xx Page 114
5. Exemples de modélisation
5.1 Couronne de pont 17x56

Dessin de définition partiel figure 5.1. 7.
Cette couronne est montée sur une boite pont de véhicule.
5.1.1 Problème posé

Tenir une tolérance de circularité de 0.01 sur un diamètre de 140 mm lors de l'usinage d'une
couronne de pont de transmission (couple 17x56) montée frettée sur le diamètre intérieur.
5.1.2 Données techniques

Matériau = 25CrMo4
Un serrage concentrique (3 mors) sur l’extérieur de la couronne est-il compatible avec les
contraintes imposées par le bureau d’étude ?
Compte tenu des efforts nécessaires au maintien correct de la couronne lors de l'usinage (calcul
du couple de serrage à partir du couple créé par les efforts de coupe), 3000N sur chaque mors,
la déformation engendrée doit être compatible avec la tolérance souhaitée.
5.1.3 Traitement d’un modèle OSSATURE
5.1.3.1 Modélisation
Structure composée de 12, 24, 36 éléments barre (discrétisation de la ligne moyenne de la

couronne).
Ossature plane. Eléments barre de section rectangulaire 23.3x22.1.
RDM 5.xx Page 115
dx = 0
D
3000 N 3000 N
B C
y
A
z x
Figure 5.1. 1
Liaisons externes : 1 noeud rotule en A, 1 noeud déplacement nul sur une direction x en D(la
rotule symbolise un mors)
Chargement : 2 charges symétriques qui symbolisent les 2 autres mors de chacune 3000 N.
Traiter le modèle avec les différentes discrétisations, comparer les résultats.
5.1.4 Traitement d’un modèle en élasticité plane (contraintes

planes) avec le module M. E. F.
5.1.4.1 Modélisation
Le modèle est une couronne d’épaisseur constante. L’étude est menée en contrainte plane (plan
de symétrie de la couronne).
L'acquisition de la géométrie est, dans ce cas, obtenue à partir d'un logiciel de dessin (DMT10)
par transfert via les formats IGES. Elle pouvait être obtenue de manière simple à l’aide du
module dessin/maillage en utilisant la bibliothèque d’éléments paramétrés du logiciel RDM
(élément de bibliothèque numéro 23).
1 - Modéliser la structure.
2 - Traiter le modèle.
3 - Analyser les résultats.
4 - Comparer avec les résultats obtenus avec le module OSSATURES.
5.1.4.2 Résultats de l’étude en élasticité plane

RDM 5.xx Page 116
Maillage de la structure (Delaunay en éléments 6 noeuds curvilignes)
Figure 5.1. 2
Epaisseur constante de 23.3 mm
Liaisons et actions mécaniques
C 3000 N
dx = 0
A B
Figure 5.1. 3
RDM 5.xx Page 117
En A et B appuis ponctuels inclinés (A : 25°, 30°, 35° ; B : 145° ,150° , 155°)
Déplacements
x Cas de charges
Figure 5.1. 4
Courbe des déplacements des génératrices du diamètre intérieur (d = 140.034 mm)

Point de départ = point D. Développement sur le périmètre.
RDM 5.xx Page 118
mm
2.05E-02
1.77E-02
1.50E-02
1.22E-02
9.46E-03
6.70E-03 mm
0.00 87.98 175.96 263.94 351.91 439.89
Module' du deplacement
MIN = 6.703E-03 mm A 311.59 mm
MAX = 2.046E-02 mm A 0.00 mm
Figure 5.1. 5
Contraintes de VON MISES
5 79 9 7 5 MPa
3 3
1 1
min= 0.23
1 1
max= 31.27
5 7 7 5 3 int= 3.10
1 3 1
1 1
3 3
1: 3.33
5 5 2: 6.44
3: 9.54
4: 12.65
5: 15.75
5 5 6: 18.86
7: 21.96
1 3 3 1
8: 25.06
1 1 9: 28.17
1 1 1 1
3 3
3 5 5 3
5 5
3 3
5 1 1 5
1 1
5 3 3 5
3 3
1 1
1 1
SIGMA VON MISES
Figure 5.1. 6
RDM 5.xx Page 119
5.1.5 Conclusion
Pour l'obtention de cet usinage, la prise en serrage concentrique n'est pas acceptable. La
déformation engendrée par cette prise de pièce est trop importante et dépasse les contraintes
imposées à la fabrication.
La solution adoptée est la réalisation d'un montage avec serrage sur flanc qui n’entraîne pas de
déformation radiale.
RDM 5.xx Page 120
Figure 5.1. 7
RDM 5.xx Page 121
5.2 Capteur d'effort "PRECIA-precia"
5.2.1 Présentation
On se propose d'étudier un capteur d'effort utilisé pour le pesage. Ce capteur est développé par
la société ATEX-PRECIA installée à PRIVAS ( 07) (figure 5.2. 1)
Figure 5.2. 1
Cette société développe toute une gamme de capteur :
Type Charge nominale Echelon Charge limite en Dimension plateau en

en kg minimum en g kg mm
CPA 3000/ 6 6 1 9 325 x 225
CPA 3000/ 10 10 1,5 15 325 x 225
CPA 3000/ 18 18 3,5 27 400 x 400
CPA 3000/ 40 40 7 60 600 x 600
CPA 3000/ 80 80 14 120 600 x 600
CPA 3000/ 200 200 28 300 800 x 800
CPA 3000/ 500 500 70 750 800 x 800
Le matériaux de ces capteurs sont :

- un Duralumin ( 2024 T4, AU4G ) pour les capteurs faibles charges ;
- un acier "à ressort" ( 45 SCD 6 ) pour les charges importantes ;
- un acier inoxydable pour les environnements agressifs.
Le montage de ces capteurs sur les balances s'effectue conformément au dessin (figure 5.2. 2).
RDM 5.xx Page 122
Figure 5.2. 2
5.2.2 Réalisation
5.2.2.1 Forme du corps d'épreuve
De type "à lunettes" : deux blocs de fixation massifs sont reliés par deux bras rigides affaiblis
en quatre points et réalisant un parallélogramme déformable élastiquement.
Le bras inférieur ne supportant pas de jauges pour des considérations de coût, est aminci pour
améliorer la linéarité des déformations du bras supérieur.
5.2.2.2 Jauges de contraintes
Ce sont des jauges pour capteurs, doubles pour faciliter les opérations de collage, en
constantan, auto-compensées en température pour le corps d'épreuve en aluminium. Elles ne
sont montées que sur les parties plates du bras supérieur pour réduire les coûts, améliorer la
précision et faciliter la protection.
5.2.3 Corps d'épreuve étudié

RDM 5.xx Page 123
75
R=11
40 30 20
192
Figure 5.2. 3
5.2.3.1 Epaisseur
e = 28 mm
5.2.3.2 Matériau
Duralumin AU4G : E = 74000 MPa, Rm = 430 MPa, Rp0,2 = 270 MPa
5.2.4 Maillage
1 3
2 4
Figure 5.2. 4 :Zonage du corps d'épreuve
Le zonage interne permet de densifier le maillage au voisinage des zones à fort gradient de
contraintes
(zones 1, 2, 3 et 4)
Le maillage s'effectue en éléments paramétriques à 6 noeuds. On discrétise la frontière en

resserrant les noeuds dans les zones 1 à 4.
figure 5.2. 5 :Maillage du corps d'épreuve
5.2.5 Montage du corps d'épreuve

RDM 5.xx Page 124
F = 120 N
S
40 mm
C
96 mm
Figure 5.2. 6
La pièce n'étant pas contrainte suivant la normale à son plan, on conduit une étude en
contraintes planes.
5.2.6 Liaisons
L'extrémité gauche est encastrée sur le bâti.
5.2.7 Chargement
Le torseur de chargement est un glisseur passant par le milieu C du plateau et de valeur
120 N.
Ce torseur réduit en C donne :

- une résultante de 120 N
- un moment de 120 x 96 = 11520 N mm
Ce moment est modélisé par deux glisseurs horizontaux en A et B de valeur
G telle que : 11520 = 40 x G
G = 288 N
5.2.8 Résultats
La déformée est telle que la zone massive droite se translate bien verticalement par rapport à la
zone encastrée.
RDM 5.xx Page 125
288 N
120 N
288 N
Figure 5.2. 7 :Corps d'épreuve déformé
mm
0.000 E+00
-2.362 E-02
-4.724 E-02
-7.087 E-02
-9.449 E-02
-1.181 E-01 mm
0.00 38.40 76.80 115.20 153.60 192.00
Deplacement
' dy
MIN = -1.181E-01 mm A 185.18 mm
MAX = 0.000E+00 mm A 0.00 mm
Figure 5.2. 8 :Déplacement vertical de la frontière supérieure
Cette figure montre clairement la déformation de type parallélogramme.

La contrainte équivalente de Von-Mises est maximale au voisinage de la zone amincie et vaut
19 MPa.
Ce résultat est à comparer à Rp0,2 = 270 MPa
RDM 5.xx Page 126
MPa
18.85
11.32
3.78
-3.75
-11.29
-18.82 mm
0.00 38.40 76.80 115.20 153.60 192.00
SIGMA xx
MIN = -18.82 MPa A 134.05 mm
MAX = 18.85 MPa A 57.79 mm
Figure 5.2. 9 :Contrainte σxx sur la frontière supérieure
Cette figure, ainsi que la suivante, montrent clairement le fort gradient de contraintes au
voisinage des zones de faible épaisseurs.
MPa
18.80
15.05
11.30
7.56
3.81
0.07 mm
0.00 38.40 76.80 115.20 153.60 192.00
SIGMA VON MISES
MIN = 0.07 MPa A 95.85 mm
MAX = 18.80 MPa A 57.79 mm
Figure 5.2. 10 :Contrainte équivalente de Von-Mises sur la frontière supérieure
RDM 5.xx Page 127
5.3 Pièce support d'ATR 42
5.3.1 L'avion ATR 42

L'avion ATR 42 peut transporter de 42 à 50 passagers avec un rayon d'action de 1700 km.
Motorisé par deux turbopropulseurs "PRATT et WHITNEY", chacun équipé d'une hélice
"HAMILTON" quadripale de diamètre 3,96 m. Cet avion de longueur 22,67 m, d'envergure 24,57
m, de hauteur hors tout 7,59 m, de masse maximale au décollage de 16700 kg a une capacité de
carburant de 4500 kg.
Pièce étudiée
Figure 5.3.1
5.3.2 Conditionnement d'air

Dans l'ATR 42, pour le confort des passagers et du personnel navigant, deux groupes de
climatisation maintiennent dans tout le fuselage un air à pression et température régulées ;
chacun des deux groupes de climatisation est alimenté en air chaud, sous pression, prélevé du
compresseur axial de l'un des deux turbopropulseurs correspondant.
RDM 5.xx Page 134
A faible régime, c'est à dire l'avion au sol ou en descente :
la veine d'air prépondérante en débit, prélève l'air "haute pression" du dernier étage du
compresseur axial et l'amène dans une électrovanne de régulation.
A régime de croisière ou de montée :
la veine d'air prépondérante en débit, prélève l'air "basse pression" d'un étage intermédiaire du
compresseur axial.
Dans tous les cas :

ces deux veines de "haute pression régulée" et de "haute pression" se raccordent sur une
tuyauterie commune située dans le longeron de voilure et la paroi intérieure du bord d'attaque de
l'aile et allant du turbopropulseur au fuselage . Cette tuyauterie commune rejoint le groupe de
climatisation placé dans la bossette de l'atterisseur principal. Cet air passe ensuite dans les
échangeurs ( régulation de température) puis dans un compresseur ( pressurisation) , ce
compresseur étant lui même actionné par une turbine utilisant une partie de cet air.
124
= = A 5,5
B 12
50
25 = = 25 58
33 2
y 2
21
23
140
36
A 30
A Section A-A
7,9
x 24
2 2
C 14 D
2 8 93
211,5
Figure 5.3. 2 : Dessin de la pièce support
RDM 5.xx Page 135
5.3.3 Données
La pièce étudiée, support de tuyauterie, mise en situation (figure 5.3. 1) et dessinée (figure 5.3.
2) fait partie d'un ensemble (figure 5.3. 3 et 4) qui maintient la tuyauterie commune le long du
longeron de voilure de l'avion.
Longeron der voilure
Support étudié
Figure 5.3. 3
Ref. Voilure
Vue suivant H
Figure 5.3. 4
5.3.3.1 Caractéristiques techniques
Nature du matériau : "Duralumin" ( AU4G ou 2024) ;

Mode d'obtention : usinée dans la masse ;
Protection anticorrosion : oxydation anodique chromique ;
Masse de la pièce : 189 g
Module de Young : E = 72000 MPa ;
RDM 5.xx Page 136
Module de Coulomb : G = 27000 MPa ;
5.3.3.2 Appuis :
Appuis pivots en C et D
5.3.3.3 Charges extérieures
? ? ?
En A : glisseur A = 600 x - 1600 y
? ? ?
En B : glisseur B = 600 x - 1600 y
5.3.4 Travail demandé
5.3.4.1 En théorie des poutres
124
= =
A
B
600 N 600 N
58
1600 N 1600 N
y
140
C D
x
93
211,5
Figure 5.3. 5
a) Elaborer un modèle d'étude qui permette de déterminer les déplacements des points
A et B ainsi que les contraintes maximales .
Faire les hypothèses nécessaires. Les préciser .
b) Traiter ce modèle à l'aide de "RDM Ossatures".
c) Analyser les résultats obtenus.
5.3.4.2 En élasticité plane
a) Elaborer un modèle d'étude qui permette de déterminer les déplacements des points
A et B ainsi que les contraintes maximales .
RDM 5.xx Page 137
Faire les hypothèses nécessaires. Les préciser .
b) Traiter ce modèle à l'aide de "RDM élasticité".
c) Analyser les résultats obtenus et les comparer à ceux obtenus à l'aide de la théorie
des poutres.
5.3.5 Résultats
Seule l'étude en élasticité plane a été conduite. Afin de montrer l'influence d'appuis
hyperstatiques sur la rigidité, on a mis la pièce successivement en appui isostatique et en appui
hyperstatique.
Bien que la pièce ne soit pas tout à fait symétrique par rapport au plan xOy ( section A-A) , on
conduit une étude en contraintes planes.
5.3.5.1 Appuis de type pivots en C et D ( hyperstatique)
5.3.5.1.1 Maillage
2 1
3
Figure 5.3. 6
Les zones 1,2 et 3 admettent une épaisseur de 2 mm

Les autres admettent une épaisseur de 8 mm
RDM 5.xx Page 138
Figure 5.3. 7 : Maillage de la structure
5.3.5.1.2 Appuis et déplacements en A et B
dx = 0,19
dx = 0,25 dy = - 0,01
dy = -0,13
A B
Figure 5.3. 8 : Déformée : pivots en C et D
RDM 5.xx Page 139
5.3.5.1.3 Contraintes équivalente de Von Mises
MPa
79.71
64.12
48.52
32.93
17.34
1.75 mm
0.00 29.99 59.98 89.97 119.96 149.95
SIGMA VON MISES

MIN = 1.75 MPa A 86.34 mm
MAX = 79.71 MPa A 40.86 mm
Figure 5.3. 9 : Contrainte de Von-Mises sur la frontière supérieure de la membrure

horizontale
Les maximums se situent en A et B dans les rayons de raccordement

MPa
98.37
78.98
59.60
40.22
20.83
1.45 mm
0.00 39.10 78.20 117.30 156.40 195.50
SIGMA VON MISES
MIN = 1.45 MPa A 13.64 mm
MAX = 98.37 MPa A 177.31 mm
Figure 5.3. 10 : Contrainte de Von-Mises sur la frontière inférieure de la membrure

horizontale
Le maximum se situe en C dans le rayon de raccordement
RDM 5.xx Page 140
5.3.5.2 Avec appui ponctuel en D (isostatique)
5.3.5.2.1 Déformée
Figure 5.3. 11 : Déformée : pivots en C et appui ponctuel en D
Déplacements du noeud B : dx = 2,28 mm ; dy = -1,06 mm

Déplacements du noeud A : dx = 0,81 mm ; dy = -0,54 mm
Ces déplacements sont beaucoup plus importants que pour la structure en appui hyperstatique.
5.3.5.2.2 Contraintes équivalente de Von Mises
RDM 5.xx Page 141
MPa
302.61
242.99
183.37
123.74
64.12
4.50 mm
0.00 29.99 59.98 89.97 119.96 149.95
SIGMA VON MISES
MIN = 4.50 MPa A 131.78 mm
MAX = 302.61 MPa A 63.60 mm
Figure 5.3. 12 : Contrainte de Von-Mises sur la frontière supérieure de la membrure

horizontale
Le maximum de contrainte équivalente se situe au voisinage du point A
La comparaison de ces résultats avec ceux obtenus avec appuis pivots en C et D montre tout
l'intérêt de mettre la structure en appuis hyperstatique. La rigidité de la structure s'en trouve
accrue.
5.4 Support de galet freineur

Ce problème est tiré de la première épreuve d’admissibilité du concours
d’Agrégation Interne de Mécanique 1996.
5.4.1 Objectif de ce problème

La façon de poser le problème à traiter est ici plus générale que la façon académique qui
consiste à se donner d’emblée les efforts appliqués à la pièce.
En effet, on se propose de vérifier la tenue aux efforts d’une pièce appartenant à un système
mécanique. On ne connaît donc pas directement les actions mécaniques s’exerçant sur les
pièces sur lesquelles on se propose de conduire une vérification.
Un calcul préalable, conduit en statique, est ici nécessaire pour déterminer les actions
mécaniques appliquées aux pièces à vérifier.
La donnée première du problème n’est d’ailleurs pas une action mécanique, mais un
déplacement de solide.
RDM 5.xx Page 142
On verra que mener à bien un tel problème nécessite un calcul itératif et l’élaboration d’un
algorithme de résolution.
5.4.2 Présentation du système mécanique

Référence au document 3 : « Fiche technique galet freineur 7302 »
Le galet freineur de type 7302 qui fait l'objet de cette étude est conçue par la société SIPA
ROLLER.
C'est un système mécanique qui se place dans les couloirs de stockage dynamique. Son but est
de maîtriser la vitesse gravitaire de déplacement des charges transportées (marchandises,
palettes, .... ). Ce galet est utilisable pour freiner des charges variant de 35 à 1000 kg. La
vitesse de déplacement de la charge est comprise entre 0,1 et 0,3 m/s.
5.4.2.1 Description du galet freineur
Références aux documents 1 : plan d’ensemble et 2 : nomenclature
Le galet freineur est supporté sur deux rails porteurs en tôle ( R) par deux colonnes support
(17) et (18) et poutres supports (15) et (16) déformables élastiquement. L'ensemble des
différents composants sont montés sur l'axe principal (11). Les ressorts (10) assurent un effort
normal déterminé entre le tambour (2) du galet et la palette de produit à transporter.
Le porte-satellite (1) est lié en rotation à la poutre support (15) par l'intermédiaire de son
embout hexagonal.
Le tambour (2) est entraîné en rotation par la palette à freiner.
Deux trains épicycloïdaux montés en série ont pour fonction d'augmenter la fréquence de
rotation du porte-mâchoires (6) d'un frein à inertie dont les mâchoires (7) freinent le tambour
(2) du galet freineur.
La liaison entre les deux colonnes support (17) et (18) et le tambour (2) du galet est assurée
par deux supports déformables élastiquement, repérés (15) et (16) .
L'objectif de ce travail consiste à effectuer l'étude en élasticité plane des supports (15) et (16)
du galet freineur, de façon à déterminer l'action mécanique entre la palette et le tambour pour
un dénivellement donné du galet par rapport à la voie de roulement.
Remarques :
- on désigne par galet, repéré (2), l'ensemble des éléments internes au galet
(repères 1 à 14, document 1) ;
- les deux supports, définis figure 5.3. 1, sont géométriquement identiques. Seul
celui repéré (15) est sollicité par le couple de freinage du galet au niveau de
l'accouplement hexagonal; le second support est repéré (16) ;
- toutes les liaisons des supports (15) et (16) avec le galet (2) et les colonnes
(17) et (18) présentent un jeu important.
RDM 5.xx Page 143
Les supports (15) et (16) sont réalisés dans un matériau synthétique nommé PA6. L'ensemble
galet est soumis :
?
- à l'action de la palette, supposée réductible à un glisseur faisant un angle ? avec y et
appliqué au point P, milieu de la génératrice supérieure du galet ;
- aux actions de deux ressorts de traction agissant entre les points H et J du support (15), et
H' et J' du support (16).
- aux actions des deux colonnes 17 et 18
Toutes les autres actions extérieures seront supposées négligeables.
y
u ?
A'
F
17
A 16
B'
O
18
15 L
C x
B 2
z
?
y
F Rayon R2
P
x
L1 L2
A K B
zone (G) zone (D) c
C
17 18
H J
FR -F R
Figure 5.4. 1 : Paramétrage du galet et des supports
Données numériques :
RDM 5.xx Page 144
Ce paragraphe regroupe les notations et valeurs numériques des paramètres du problème.
Certains d'entre eux seront définis ultérieurement.
Paramètres géométriques :
? rayon du galet : R2 = 42,5 mm;

? entraxe des colonnes (17) et (18): AB = A'B' = e = 110 mm
? distance des supports (15) et (16): AA' = BB' = L = 120 mm
? L1 = 68 mm ? L2 = 42 mm ?l = 6 mm
? d1 = 34 mm ? d2 = 10 mm ?c = 13 mm
? h = 32 mm ? a = 11 mm ? b = 3 mm
Contact galet-palette :
? action palette-galet (valeur initiale considérée pour le démarrage du calcul itératif) :
F = 450 N
? angle de frottement : ? = 36°
Ressort :
? L0 = 42 mm ?k = 22 N/mm
Matériau synthétique PA6 des supports :
? module de Young : E = 1850 Mpa (N/mm²)

? limite élastique en traction: Re = 49 Mpa (N/mm²)
5.4.3 Etude statique

L'intérêt de cette question repose sur la prise en compte des jeux dans les liaisons pour procéder
à l'étude statique du galet freineur, étude qui sans cette hypothèse présenterait un caractère
hyperstatique.
En plus des hypothèses exposées ci-dessus, on suppose pour cette étude statique (figure 5.4. B-
1) :
- que les supports (15) et (16) sont indéformables;
- que leurs liaisons avec les colonnes (17) et (18) sont réalisées avec un jeu j important
(jeu au rayon: j = R - r), qui permet de considérer ponctuels les contacts avec ces
pièces;
- que le frottement peut être négligé au contact des supports et des colonnes.
- que l'action de la palette sur le galet est modélisable par un glisseur noté :
??? ? ?
F ? ? F u?
?Fpalette ? galet? ? ? ? ?
?? 0 ??
P
5.4.3.1 Etude statique du galet
RDM 5.xx Page 145
a) L'ensemble constitué par les supports (15) et (16) et les axes (17) et (18) étant supposé
rigide, proposer des liaisons entre le galet (2) et cet ensemble conduisant à un modèle
isostatique et assurant l'équilibre du galet.
Réponse :
Pour assurer l’équilibre du galet, sa liaison avec l’ensemble du bâti supposé rigide doit être
globalement un encastrement.
Il est réalisé de façon pratique au moyen de deux liaisons avec les supports (15) et (16). Seul le
?
contact avec (15) est susceptible d’exercer la composante de moment suivant z assurant la
fonction de freinage du galet. On conviendra de représenter cette liaison par une liaison
?
sphérique à doigt de centre (C; z ) et la liaison avec (16) par une linéaire annulaire d’axe
?
( C'; z ) .
b) Procéder à l'étude statique du galet (2) et en déduire les actions exercées sur les
supports (15) et (16).
Application numérique: déterminer la valeur numérique de ces actions.
Réponse :
? ? ? ? ? ? ? ? ? ? ? ? ?
On note ( t , u , z ) la base déduite par rotation d’angle ? de ( x , y , z ) : ( x , y , z ) ? ( t , u , z )
Le galet est soumis à trois actions:
?? ? ?
? F ? ? F u?
- Action de la palette: ? Fpalette ? ( 2 ) ? ? ? ? ?
?? 0 ??
P
? ? ? ? ?
?X t ? Y u ? Z z ?
- Action du support (15): ? F15 ? 2? ? ? ? ?
?? N z ??
C
? ? ? ?
?X' t ? Y' u ?
- Action du support (16): ?F16 ? 2? ? ? ? ?
?? 0 ??
C'
Equilibre en O:
? ? ? ? ? ?
Résultante en projection sur ( t , u , z ) : Moment en projection sur ( t , u , z ) :
? L
? X ? X' = 0 ( X ? X' )
2
=0
?? L
? Y + Y'- F = 0 (Y - Y') = 0
? 2
?Z ? 0 N ? FR2 sin ? ? 0
??
RDM 5.xx Page 146
X ? X' ? 0
F X ? X' ? 0 N
Y ? Y' ? Y ? Y' ? 225 N
D'où: 2 AN :
Z? 0
Z? 0
N ? 11,24 Nm
N ? FR 2 sin ?
5.4.3.2 Etude statique du support (15)
La figure 5.4. 2 montre, avec une amplification importante des jeux et déplacements, la
géométrie de l'ensemble support (15) et colonnes (17) et (18) après application des efforts sur le
galet (2).
Les points A et B sont les centres des sections des arbres (17) et (18), de rayons r. Les points
A1 et B1 sont les centres des alésages du support, de rayons R. Les entraxes AB et A1 B1
présentent une différence de longueur notée ?, définie par: AB = e et A1 B1 = e + ? , avec ? < j =
R-r où j représente le jeu au rayon de la liaison.
? ? ? ? ? ?
Les contacts sont ponctuels, de normales n A et n B , telles que ( y , nA ) = ? A et ( y , nB) = ? B .
nA
?A
nB ?B
17 15 18
R
A B
A1 B1
r
y
C
x
Figure 5.4. 2 : Modélisation géométrique du montage du support (15) et des axes

(17) et (18)
1°)
a) Réaliser un schéma du support (15) et représenter graphiquement les actions mécaniques
qui lui sont appliquées.
b) Exprimer la relation géométrique liant les inclinaisons ? A et ? B des normales au contact
entre le support (15) et les axes (17) et (18) en fonction de j et ?.
c) Que devient cette relation si ? = 0 ?
Réponse :
RDM 5.xx Page 147
nB
nA
FA FB
?A ?B
17 15 18
A B
H
A1 B1
y
C
x -N z
?
F/2
Figure 5.4. 3
? ?
Soit AB ? e x l'entraxe des arbres (11) et (12) constituant le bâti; l'entraxe des alésages du
? ? ? ?
support a pour expression A1B1 = (e + ? ) x' avec ( x , x ') ? ? où ? est un angle de faible valeur.
D'autre part, le support (15) étant en contact avec les axes (17) et (18), on a les propriétés:
? ? ? ? ? ?
? A1A ? (R ? r ) nA ? j nA ? j( ? sin ? A x ? cos ? A y )
? ? ? ? ? ?
?? B1B ? ( R ? r ) nB ? j nB ? j(? sin ? B x ? cos ? B y )
En utilisant la propriété vectorielle suivante :

? ? ? ? ? ? ?
( AB? BB1? B1A? A1A ) ? 0 , on obtient en projection sur ( x , y ) :
?e ? j(sin ? A ? sin ?B) ? (e ? ? ) ? 0
?
? j(cos ?A ? cos ? B) ? ? (e ? ?) ? 0
?
soit en particulier : sin ?B ? sin ?A ? (I)
j
La seconde relation confirme le caractère ‘petit’ de l’angle ? ? posé en hypothèse (ordre de

grandeur de j/e)...
Les dimensions e, R, r et ? étant données, la position du support (15) par rapport aux colonnes
est complètement définie par les deux angles ? A = ? B. La relation (I) liant ces paramètres met
en évidence la mobilité de degré 1 du support par rapport au bâti supposé fixe.
Dans le cas particulier ou ? = 0, on vérifie que ? A = ? B .
RDM 5.xx Page 148
2°) Exprimer les conditions d'équilibre du support (15). En déduire:
? ?
a) Les expressions des actions des axes en A et B, notées RA, et RB en fonction de
F, e, R2 , c, ? , ? A, et ? B .
b) Exprimer les conditions liant ? A et ? B à l'équilibre
3°)
a) Etudier le cas particulier où les entraxes des arbres et des alésages sont égaux (?
= 0).
b) Application numérique: en déduire les valeurs numériques des actions RA et R B .
? ?
L'hypothèse ? = 0 et les expressions correspondantes des actions RA et RB
sont retenues pour la suite de l'étude.
Réponse :
Dans cette question de statique , on suppose que les petits déplacements dus aux jeux entre le
support et les colonnes ne modifient pas sensiblement la géométrie de l’ensemble, donc les
relations traduisant l’équilibre. Cela revient à supposer que A et A1 sont confondus, de même
que B et B1.
Equilibre du support (15) :

? F? ?
? u ?
- Action du galet: ?F(2 ) ? (15)? ? ? 2
?
?
? ? FR sin ? z ?
C? 2 ?
?? ? ?
? RA ? RA n A?
- Action de l'axe (17): ? F(17 ) ? ( 15) ? ? ? ? ?
A ?? 0 ??
?? ? ?
? RB ? RB nB?
- Action de l'axe (18): ? F(18) ? ( 15) ? ? ? ? ?
?? 0 ??
B
? F
? ? ? ? RAsin? A ? RBsin? B + 2 sin? = 0 ( II)
Résultante en projection sur ( x , y ) : ?
F
? RAcos ? A ? RBcos? B ? cos ? = 0 ( III )
? 2
?
? ? ? ? ? ? F ? ?
Moment en C, en projection sur z : HA ? RA nA ? HB? RB nB ? HC? ?Nz? 0,
2
soit après développement : ? e (RA cos ? A ? RB cos ? B) ? F sin ? (2 R 2 ? c ) ? 0 (IV)
RDM 5.xx Page 149
Il est possible, en combinant les équations (III) et (IV) d’exprimer deux relations liant RA et ? A
d’une part, RB et ? B d’autre part :
? ? cos ? 2 R ? c ??
? RA cos ? A ? F ? ? sin ? ?? 2 ? (III')
? ? 4 ? 2 e ? ??
?
? RB cos ? B ? F ? cos ? ? sin ? ?? 2 R2 ? c ?? ? (IV' )
?? ? 4 ? 2 e ? ??
?
En replaçant RA et RB dans (II), on exprime une seconde condition sur ? A et ? B à l'équilibre :
? 1 2R 2 ? c ? ? 1 2 R2 ? c ? 1
tan ? A ? ? ? ? ? tan ?B ? ? ??? (II’)
?4 tan ? 2e ? ?4 tan ? 2e ? 2
à rapprocher de la condition géométrique (I) :
?
sin ?B ? sin ?A ?
j
Conclusion :
Les quatre relations (I), (II), (III) et (IV) lient les quatre paramètres RA, RB, ? A et ? B. On
déduit de la résolution de ce système la position du support à l’équilibre, ainsi que les actions
dans les liaisons en A et B.
La résolution dans le cas général n’est possible que numériquement. Seul le cas particulier ?=0
conduit à une solution analytique simple, hypothèse retenue dans la suite.
5.4.4 Etude en élasticité plane

Utilisation du logiciel d'élasticité plane «RDM Le Mans» pour étudier les contraintes et
déformations du support (15).
Les conditions de l'étude sont celles définies par la figure 5.3. 1, avec pour efforts :
? accouplement hexagonal :
R|F x ? F y U|
? ?
lF q? S|
galet ? support
x
V| avec F = 132 N
?
y
Fy = -182 N et Mz = -11.2 Nm
CT W
x
M z z
R U|
q? |S| X x? Y
? ?
? action de l'arbre (17) : lF V| avec X

A A y
(17) ? support ? = -55 N et YA = 75 N
AT W
A
0
R| X x ? Y ?
y U|
?
? action de l'arbre (18) : lF q S|

?
(18) ? support
B
?
B
V| avec X
B= -78 N et YB = 106 N
BT W
0
RDM 5.xx Page 150
? action du ressort hélicoïdal : FR = 135 N
Partant d’un dessin exécuté avec un logiciel graphique (Autocad ou DMT), on importe le
dessin sous forme de fichier d’échange DXF dans RDM.
On procède dans un premier temps à la définition géométrique du modèle, à son maillage en

éléments triangulaires curvilignes isoparamétriques à six noeuds. Le maillage a été resserré
dans la zone où l’on s’attend à un fort gradient de contraintes.
La figure du document 4 montre le résultat du maillage.
On procède ensuite à la définition :

- des épaisseurs ;
- du matériau (E = 1850 MPa)
Il s’agit ensuite de définir les conditions aux limites (liaisons avec les colonnes (17) et (18) et
actions extérieures).
5.4.4.1 Description des liaisons
La figure 5.4. 4 présente le groupe des boutons disponibles dans le logiciel pour la description des
liaisons des noeuds avec le bâti.
Noeud en appui élastique suivant y

Noeud en appui élastique suivant x
Déplacement imposé suivant y
Déplacement imposé suivant x
Noeud immobilisé en translation
Noeud en appui simple de normale quelconque
Noeud en appui simple de normale y
Noeud en appui simple de normale x
Figure 5.4. 4 : Descriptions du groupe des liaisons
Ces liaisons appartiennent à l'un des quatre types suivants :
- type 1: les trois premiers boutons correspondent à des appuis simples de

? ? ?
normale respectivement x (horizontale), y (verticale), et u de direction
quelconque, définie par l'utilisateur;
- type 2: le quatrième bouton correspond à une liaison du noeud interdisant son
déplacement en translation;
RDM 5.xx Page 151
- type 3: le cinquième et le sixième bouton correspondent à des liaisons imposant
un déplacement connu du noeud, défini par l'utilisateur, respectivement suivant
? ?
x (horizontal, noté dx) et y (vertical noté dy)
- type 4: les deux derniers boutons correspondent à des appuis élastiques, pour
lesquels l'action de liaison est proportionnelle au déplacement du noeud,
? ?
respectivement suivant x ( F x = - k? x) et suivant y ( F y = - k? y)
Question :
Sur la figure du document 4, définir et justifier le choix des liaisons retenues pour l'étude du
support de galet.
Indiquer avec précision quels noeuds font l'objet de liaisons, la nature de ces liaisons et les
propriétés géométriques associées.
Réponse :
Les liaisons au niveau de chaque colonne sont modélisées par des appuis simples sur les noeuds
?
situés à proximité immédiate de la direction u définie par l’angle ? .
Le fait de retenir plusieurs noeuds pour chacune des deux colonnes permet d’assurer l’équilibre
et de valider ultérieurement la position du support, par consultation des actions de contact
déterminées lors du calcul.
5.4.4.2 Description des charges
La figure 5.4. 5 présente le groupe des boutons disponibles dans le logiciel pour la description des
charges extérieures appliquées au support.
Action de pesanteur
Champ de pression
Champ d'effort
Glisseur appliqué à un noeud
Figure 5.4. 5 : Description du groupe des efforts
Ces charges appartiennent à l'un des quatre types suivants :
- type 1 : le premier bouton correspond à un glisseur appliqué à un noeud, de

composantes définies par l'utilisateur;
- type 2 : le second bouton correspond à une charge linéique, de composantes
définies par l'utilisateur et appliquée sur une partie du contour de la pièce;
- type 3 : le troisième bouton correspond à un champ de pression uniforme sur une
partie du contour extérieur de la pièce;
- type 4 : le dernier bouton correspond à l'action de pesanteur, définie par les données
relatives à la géométrie de la pièce (volume) et sa masse volumique
(matériau).
RDM 5.xx Page 152
Sur la figure du document 4 , définir et justifier les charges à appliquer au support pour procéder
à son étude. Indiquer avec précision quels sont les noeuds soumis à ces actions, leur type et leurs
propriétés géométriques et numériques.
Réponse :
Accouplement hexagonal avec le galet :
Les efforts effectivement appliqués au niveau de l’accouplement hexagonal ne peuvent pas
être connus avec précision ; on choisira par conséquent un mode de chargement dont le torseur
associé est conforme à celui résultant de l’action du galet. Par exemple:
? trois glisseurs appliqués aux sommets de l’hexagone H1, H2 et H3, de résultantes

?
F
identiques, égales à;
6
? deux glisseurs appliqués aux sommets H0 et H3 de l’hexagone, de résultantes
FR2 sin ?
identiques de normes égales à , disposés tangentiellement au cercle de centre
2?
C et de rayon ? .
Il ne faudra pas s’attacher aux résultats trop proches de ces points ( principe de Saint-Venant).
Action du ressort :
L’action du ressort sera modélisée par deux glisseurs opposés, appliqués aux noeuds H et J et
de norme FR=135 N pour démarrer le processus de calcul.
5.4.4.3 Conditions de l'étude
Préciser et justifier en quelques mots le type d'étude envisagée: contraintes planes ou

déformations planes.
Indiquer les paramètres à consulter au niveau des résultats pour vérifier l'aptitude de la pièce à
satisfaire aux exigences du fonctionnement.
Réponse :
Compte tenu de la faible épaisseur du support, de la symétrie plane de sa géométrie et des
efforts appliqués, une étude en contraintes planes est retenue.
En règle générale, une étude en éléments finis suppose que les conditions aux limites (liaisons
et efforts) sont connues.
Dans le cas particulier de ce problème, l’action FR du ressort dépend de la déformation du
support, qui dépend de l’action de la palette, elle -même fonction de l’action du ressort. Une
telle situation ne trouve sa solution qu’au prix d’une succession d’études approchant
progressivement de la solution effective.
La démarche utilisée est la suivante :

Les résultats de l’étude en statique constituent un point de départ ‘raisonnable’ de la recherche
itérative.
On suppose connues les valeurs des efforts F et FR ; l’application du calcul permet de
déterminer le déplacement du point C et la longueur effective du ressort. A partir de ces
résultats, on évalue de nouveau l’effort F pour atteindre un déplacement vertical de
l’accouplement hexagonal de 5 mm, et la valeur de FR à partir de l’allongement du ressort...
Ce processus est répété jusqu'à ce que le déplacement du point C soit voisin de 5 mm et
l’action du ressort en cohérence avec son allongement.
RDM 5.xx Page 153
Un algorithme possible de résolution est alors le suivant :
Donnée du déplacement
vertical ? du point P
Donnée de l'effort en P
Détermination par la statique

des actions mécaniques
appliquées à la pièce :
- efforts en A et B
(direction et norme)
- de l'action du ressort
Calcul par logiciel du déplacement

vertical ? i du point P
non ? i - ? ???
oui
Obtention des résultats
Figure 5.4. 6 : Algorithme de résolution
5.4.5 Modèle définitif

Le modèle d’étude en terme de :
- maillage ;
- appuis ;
- efforts
est conforme à celui schématisé sur la figure 5.4. 8.
Toutes ces données sont répertoriées dans le fichier de calcul : galet.cal, ne pas oublier que le
module de Young du matériau est de 1850 MPa.
Trois itérations sont nécessaires pour obtenir une convergence satisfaisante du processus.
RDM 5.xx Page 154
u
?
u
+A +B
H J
FR -F R
FR2 sin ?
y
2?
?
H0 C H3
F/6
H2
H1
H3
Figure 5.4. 7 : Modèle d’étude
RDM 5.xx Page 155
36° u 36°
u
+A +B
H J
185 N 185 N
Fx = 44 N
617 N Fy = -60 N
9 mm
R 2 = 42 mm
H0 C H3
F = 450 N
617 N
? = 36°
? = 9 mm H2
H1
Fr = 185 N
H3
Figure 5.4. 8 : Valeurs numériques
RDM 5.xx Page 156
Détail suivant E
du support 15
C D A-A
C-C D-D
15 14 1 3 4 5 6 7 2 12 13 16 11
B B
C D
B-B R
17
Vue F couvercle 13 enlevé
7
10
Sens de rotation
du tambour
A A
9
7'
embout hexagonal
18
Galet freineur 7302 Echelle #1 DOCUMENT 1

RDM 5.xx Page 157
18 1 Colonne support
17 1 Colonne support
16 1 Support droit PA 6
15 1 Support gauche PA 6
14 2 Roulement 42 x 10 PX C 9021 M
13 2 Couvercle PA 6.6 rouge
2 1 Embout PA 6
11 1 Axe principal
10 2 Ressort Inox
9 2 Ressort de rappel Z6 CN 18-08
8 2 Masselotte Pb 2
7 2 Mâchoire
6 1 Porte-mâchoire PA 6
5 2 Satellite PA 6
4 1 Porte satellite deux axes PA 6
3 3 Satellite PA 6
2 1 Tambour PA 6.6 rouge
1 1 Porte satellite trois axes PA 6
R 2 Rail porteur Tôle galvanisée
Repère Nbre Désignation Matériau Observations
Date:
Echelle 1
GALET
FREINEUR 7302
DOCUMENT 2 NOMENCLATURE PARTIELLE
RDM 5.xx Page 158
Figure 5.4. 9 : Dessin du support de galet
RDM 5.xx Page 159
6. Annexe
RDM 5.xx Page 160
6.1 Limites actuelles du logiciel M.E.F.
6.1.1 Limites de calcul
Le logiciel est capable de traiter, à l’heure actuelle, des études comportant au maximum :
3000 Noeuds
1700 Eléments triangles ou quadrangles
15 Matériaux différents
15 Epaisseurs différentes
300 Liaisons
15 Cas de charges
300 Charges par cas
10 Combinaisons de 5 cas de charges
20 Modes propres demandés
6.1.2 Sauvegarde des études
Les études sont sauvegardées dans deux fichiers :

Un fichier modèle de la géométrie avec l’extension .GEO.
Un fichier modèle de calcul avec l’extension .CAL
Il est impossible de modifier la géométrie à partir du fichier de calcul.
RDM 5.xx Page 161
6.2 Modélisation d’une articulation
La modélisation des conditions limites pour une liaison de type pivot peut être obtenue en utilisant
la liaison nodale appui incliné sur les différents noeuds de la frontière.
La difficulté réside dans l’orientation de l’ensemble des appuis. La position et l’orientation de ces
appuis dépendent de la discrétisation du cercle à partir duquel est construit le modèle maillé.
Pour un cercle discrétisé à quatre éléments (le découpage étant uniforme, les quatre secteurs
sont définis par un même angle), il est possible de placer au moins quatre appuis inclinés orientés
à 0°, 90°, 180° et 270° par rapport à la direction horizontale. Le premier noeud doit être dans ce
cas placé sur la direction horizontale passant par le centre du cercle discrétisé (cas général
lorsqu’il n’y a pas de point défini comme point à mailler sur le cercle).
Dans le cas ou il existe un point à mailler, ce point devient le point de départ de la discrétisation.
Dans ce cas, la position angulaire de ce point sur le cercle étant connue, il suffit de modifier les
valeurs précédentes en leur rajoutant la valeur de cette position.
Pour une discrétisation plus dense, il sera possible de placer autant de liaisons que nécessaire.
Point à mailler
situé à 30°
y y
y
x x x
Après calcul une étude des actions de liaisons est nécessaire. Dans la modélisation ces liaisons
sont bilatérales. De ce fait, certaines deviennent inopportunes car le sens de l’action de liaison
est incompatible avec le comportement d’une liaison pivot. Le sens doit être vers l’intérieur de la
matière du modèle étudié et non l’inverse.
y y
x x
RDM 5.xx Page 162
Il peut paraître alors utile, de supprimer les liaisons qui « perturbent » et refaire un nouveau
calcul jusqu’à obtenir un modèle de comportement cohérent avec la liaison modélisée.
y y
x x
Cette remarque qui concerne la vérification du comportement des liaisons peut être transféré à
toute autre liaison.
RDM 5.xx Page 163
6.3 Courbes CETIM
Courbes CETIM
d’après
Guide du dessinateur
Les concentrations de contraintes
Plaque avec deux entailles à fond semi-circulaire - TRACTION
r
P
D
d
t
r/t
0,03 0,04 0,05 0,06 0,07 0,08 0,09 0,10 0,12 0,14
6 0,16
Kt 0,18
0,20
0,22
5
0,26
0,30
0,40
4
0,50
0,60
0,80
3 1,00
1,20
1,40
2,00
3,00
2 4,00
5,00
6,00
10,00
1
0,4 0,5 0,6 0,7 0,8 0,9 1,0
d/D
RDM 5.xx Page 164
Courbes CETIM
Plaque avec deux entailles à fond semi-circulaire - FLEXION
M M
D
d
t
6 0,03 0,04 0,05 0,06 0,08 0,10 0,12 0,14 r/t

0,16
Kt
0,18
0,20
0,22
5
0,26
0,30
0,40
4
0,50
0,60
0,80
3
1,00
1,40
2,00
3,00
2 4,00
5,00
6,00
10,00
1
0,4 0,5 0,6 0,7 0,8 0,9 1,0
d/D
RDM 5.xx Page 165
RDM 5.xx Page 166
7. Bibliographie
(1) Yves DEBARD

Notice du logiciel "RDM"
(2) J.P. FAURIE et al.

Guide du dessinateur. Les concentrations de contraintes. CETIM
(3) J.P. HENRY et F.PARSY

Cours d'élasticité.
DUNOD Université 1982
(4) M. KERGUIGNAS et G. CAIGNAERT

Résistance des Matériaux.
DUNOD Université 1977
(5) G. SPINNLER
Cours polycopié de "Dimensionnement des organes de machine"
Ecole polytechnique fédérale de Lausanne 1985
(6) S. LAROZE et J.J. BARRAU: Mécanique des structures.

Tome 1. Solides élastiques plaques et coques
2e Edition
EYROLLES-MASSON 1988
(7) A. POTIRON . Cours de Mécanique des Milieux Continus .

Centre de l' ENSAM d'Angers
(8) G. DHATT et G. TOUZOT

Une présentation de la méthode des éléments finis
Collection Université de Compiègne
MALOINE (1981)
(9) The Finite Element Methode in Engineering Science

O.C. ZIENKIEWICZ
McGRAW-HILL (1977)
(10) Concept and Applications of Finite Element Analysis

Robert D. COOK
John WILEY & SONS, Inc. (1974)
(11) Modélisation des structures par éléments finis. Volume 1: Solides plastiques
Jean Louis BATOZ et Gouri DHATT
HERMES (1990)
(12) Computation methods for the solution of engineering problems.

C.A. BREBBIA et A.J. FERRANTE
RDM 5.xx Page 167
PENTECH PRESS CRANE, RUSSAK & CO., Inc New York
(13) Eléments finis et CAO.

J.C SABONNADIERE et J.L. COULOMB
HERMES (1986)
RDM 5.xx Page 168
RDM 5.xx Page 169

Cours Complet RDM 2

Transféré par

Droits d'auteur :

Formats disponibles

Cours Complet RDM 2

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

Cours Complet RDM 2

Transféré par

Droits d'auteur :

Formats disponibles

RDM 5.

RDM 5.xx Page 1

RDM 5.xx Page 2

RDM 5.xx Page 4

Figure 1. 1 : Contraintes autour du point A

RDM 5.xx Page 5

1.1.3 Contrainte normale et tangentielle

- sur la normale, on obtient la contrainte normale :

σ 11n1+σ 12n2+σ 13n3   n1(σ 11n1+σ 12n2+σ 13n3) 

- sur le plan tangent, on obtient la contrainte tangentielle τ telle que:

Figure 1. 2 : Contrainte normale σ et tangentielle τ en un point A

1.1.4 Contraintes principales et directions principales

- les trois valeurs propres σi sont réelles (distinctes ou confondues) ;

RDM 5.xx Page 6

Figure 1. 3: Contraintes principales autour du point A

Traduction mécanique : si les trois contraintes principales σ1, σ2 et σ3 sont distinctes, il

Dans le repère principal la matrice des contraintes s'écrit alors :

Détermination des contraintes principales et des directions principales (figure 1. 3)

- Les contraintes principales sont déterminées en écrivant que le déterminant

On aboutit à l'équation caractéristique : -σ3 + I1 σ² - I2 σ + I3 = 0

I1 = σ11 + σ22 + σ33 = trace de σ(A)

RDM 5.xx Page 7

1.1.5 Etats de contraintes particuliers

Dans le repère principal, le tenseur des contraintes se réduit à :

Traction simple si σ1 > 0, compression simple si σ1 < 0

b) Etat de cisaillement simple

Les contraintes principales sont égales à :

Dans le repère principal, le tenseur des contraintes admet la forme suivante :

RDM 5.xx Page 8

c) Etat de contraintes planes

σ13 = σ23 = σ33 = 0

C'est le cas des plaques planes chargées dans leur plan.

Le tenseur des contraintes s'écrit alors dans une base quelconque :

Dans la base principale :

d) Etat de contrainte dans une section droite de poutre

Traction (compression) simple :

RDM 5.xx Page 9

Torsion avec sections circulaires :

1.1.6 Cercles de Mohr

RDM 5.xx Page 10

Figure 1. 4a : Tricercle de Mohr

1.1.6.2 Cercles de Mohr en contraintes planes

En contraintes plane σ3 = 0 ainsi la figure ci-dessus devient :

Figure 1. 4b : Tricercle de Mohr en contraintes planes

RDM 5.xx Page 11

Figure 1. 5 : Construction de Mohr

Le vecteur contrainte admet pour expression :

et la contrainte tangentielle définie ici par :

en passant en arc double, il vient :

RDM 5.xx Page 12

Dans le plan de Mohr, le lieu de M lorsque θ varie est un cercle de centre :

1.1.6.3 Construction de Mohr

Convenons de classer les contraintes principales dans l'ordre suivant :

Remarque : Ne pas oublier que σ3 = 0. Il y a donc trois cercles de Mohr.

Figure 1. 6 : Point A sur la surface d'un solide

RDM 5.xx Page 13

Signification de la relation (2) :

Ainsi, lorsqu'il n'y a pas de chargement tangentiel à la surface au point A, la normale à la

Signification de la relation (3) :