Elasticite Et RDM Part II
Elasticite Et RDM Part II
Elasticite Et RDM Part II
Remarques:
Mécanique générale, élasticité, résistance des matériaux sont des appellations
traditionnelles gardées par respect de l’histoire de la science et ne traduisent
pas le sens propre de ces trois modèles de formulation de science de la mécanique:
Appellations traditionnelles
Mécanique des solides indéformables : mécanique générale
Mécanique des solides élastiques avec faibles déformations : élasticité
Mécanique des milieux curvilignes élastiques : RDM
3
S. MEZLINI - T. DAAMI – A. DOGUI 3
Calculateurs Temps Réel
Résistance des Matériaux
II Milieux curvilignes: formulation torsorielle des solides déformables de type poutres
1. Hypothèses fondamentales de formulation
H1 : hypothèses sur la géométrie:
Poutre: un solide dont les dimensions di des sections S(s) perpendiculaires à la
fibre moyenne sont faibles par rapport à la dimension longitudinale.
S(s) Est tangent à la fibre moyenne s
G1
G1
G L x2 x( s )
•G
d i L x s y(s)
di x3 z (s)
G0
x s
G0
Pas de variation brusque des sections:
On définit la ligne moyenne: courbe passant par les centres de gravités G(s)des
sections S(s). Les sections droites S(s) du solide sont perpendiculaires à la ligne moyenne.
5
S. MEZLINI - T. DAAMI – A. DOGUI 5
Calculateurs Temps Réel
Chapitre 3: Problème de Saint Venant
II Milieux curvilignes: formulation torsorielle des solides déformables de type poutres
1. Hypothèses fondamentales
H4: hypothèse sur les efforts et le mouvement : hypothèse de Saint Venant
Suffisamment loin des points d’application des charges, des zones de
discontinuité géométrique et de variation brusque de section, la solution
(état de contrainte et de déformation) d’un problème de mécanique des
solides élastiques de type poutre ne dépend que du torseur équivalent
des efforts appliqués et non de la manière précise dont ces efforts sont
appliqués.
Variation brusque de la section
Zones non concernées par Saint Venant
6
S. MEZLINI - T. DAAMI – A. DOGUI 6
Calculateurs Temps Réel
Problème de Saint Venant
II Milieux curvilignes: formulation torsorielle des solides déformables de type poutres
1. Hypothèses fondamentales
Exemple : Simulation de l’état des contraintes selon le modèle de l’élasticité montrant
les zones d’application de la RDM
F F
2 2 F
Discontinuité
géométrique
100MPa
Suffisamment loin des points d’application des charges et des zones de discontinuité
géométrique, les deux répartitions donnent le même état de contraintes 100MPa
7
S. MEZLINI - T. DAAMI – A. DOGUI 7
Calculateurs Temps Réel
Résistance des Matériaux
II Milieux curvilignes: formulation torsorielle des solides déformables de type poutres
8
S. MEZLINI - T. DAAMI – A. DOGUI 8
Calculateurs Temps Réel
Résistance des Matériaux
II Milieux curvilignes: formulation torsorielle des solides déformables de type poutres
2. Conséquences des hypothèses : modèle milieu curviligne
F F F
20 000 m
T Construction métallique
250 100
A B Cr Couple résistant
t Tensions de la courroie
Cr t
T
T
S. MEZLINI - T. DAAMI – A. DOGUI 9 9
Calculateurs Temps Réel
Résistance des Matériaux
II Milieux curvilignes: formulation torsorielle des solides déformables de type poutres
2. Conséquences des hypothèses : Charge uniforme repartie
modèle milieu curviligne
10
S. MEZLINI - T. DAAMI – A. DOGUI 10
Calculateurs Temps Réel
Résistance des Matériaux
II Milieux curvilignes: formulation torsorielle des solides déformables de type poutres
2. Conséquences des hypothèses : modèle milieu curviligne
c. Problème non régulier:
Pertes d’informations sur les discontinuités géométriques dans le
solide et comment les efforts sont appliqués.
L’application qui à toute action on associe un torseur n’est pas bijective.
Ceci n’est pas contraignant pour exploiter cette modélisation simple qui donne
des résultats appliqués très économique en respectant les hypothèses .
Correction:
Coefficient de concentration des contraintes pour les discontinuité
géométrique
Modèle de Hertz : effet de contact des solides déformables
11
S. MEZLINI - T. DAAMI – A. DOGUI 11
Calculateurs Temps Réel
Résistance des Matériaux
III milieu curviligne et description des efforts :
1. Modélisation des efforts extérieurs
Cj s
sL
f s
s F i : Force concentrée
R : quelconque R : quelconque
M : quelconque M= 0
Rotule
Encastrement Sphérique
R : quelconque
y
x M= 0
Mx 0 R R Ry
pivot My 0
Mz 0 Appui Simple
M int
•
int ext/gauche ext/droite
int : traduit la réaction de la partie 1 (gauche)
sur la partie 2 (droite)
x3
int 1/2 2/1
R int
Les sollicitations
R s N T
int s T
M s M t M f
R(s )
int
N : Effort normal
T : Efforts tranchants
M t : Moment de torsion s
G(s)
M f : Moments de flexion N
Mt
M f (s)
M (s)int
R( A) R( B) Fi Rk f (s)ds 0 F
sA
M ( A) OA R ( A) M ( B ) OB R ( B ) (Ci OH i Fi )
sB
i
M k OP k Rk OG f (s) ds 0
sA
y
n j
Ces équations vont être exploitées pour l’étude par p
morceau de poutre (tronçons de poutre AB) entre les
points des discontinuités géométriques du milieu
curviligne et des charges (actions concentrées f ( y ) m x h
R( s) R( s ds) f ( s)ds 0
OG(s) f (s) ds 0
sA s2
M ( s) M ( s ds) OA R( s) OB R( s ds)
s1
sB sB
PB d
( OG f (s))ds M OG R P A
sA
s ds (M OG R)ds
A
sB
dM dR
( R OG )ds
sA
ds ds
sB
d M
sB
dR
s ()ds s ds R OG ( ds f (s))ds d ' où
dM
R ( s )
A A ds
0
R int ( s ) Rk R int ( s )
M int ( s ) M k M int ( s )
20
S. MEZLINI - T. DAAMI – A. DOGUI 20
Calculateurs Temps Réel
Résistance des Matériaux
III Milieux curvilignes et description des efforts
4. Equation d’équilibre locale
Application aux poutres rectilignes :les équations d’équilibre donnent
dTy dTz dM y dM z
dN
fx fy fz Tz T y
dx dx dx dx
dx
( L a) F
Ty T
y F L
a x
I + x
A B
aF ( L a) F
YB , YA dM z Mz aF
L L m T
dx L
dTy
0 Ty cte c M z mx b
dx
c, m,b constantes qu’on détermine facilement Grace
x
aux conditions aux limites x=0, x=a, x=L ( L a) F
a
L
21
S. MEZLINI - T. DAAMI – A. DOGUI 21
Calculateurs Temps Réel
Résistance des Matériaux
III Milieux curvilignes et description des efforts
4. Equation d’équilibre locale
Application aux structures en treillis (barres bi articulés, actions sur les nœuds).
Rj F1 F2 F3
L/4 L/4 L/4 L/4
j 4
6
2
M f3 x ds2 11
6
Découplage des effets
12 x2 13 x3 Contrainte de cisaillement, effet de M t et T
11 x1 Contrainte normale , effet de M f et N
11 x1 , x2 , x3 0 0
R int N x1
0 0 0 int G
M int M f 2 x 2 M f3 x3 G
0 0 0
σ=σ.n 0 avec
n 0, n2 , n3
Les équations d’équilibre:
σ11
divσ=0 =0 σ11 =σ11 x 2 , x 3
x1
2σ11 2σ11 2σ11
Les équations de Beltrami: = 2 = =0
x 2 x 3 x 2x 3
2
M f 2 M G 02 x1RG 03 11 x3ds
a x3ds b x2 x3ds c x32 ds
M f 3 M G 03 x1RG 02 11 x2 ds
a x2 ds b x22 ds c x2 x3ds
On choisit aussi x2,x3 selon les directions principales d’inertie: x2 x3ds 0
On aura donc:
N
N RG 01 RG11 a.S a
S
Mf2
M f 2 M G 02 x1RG 03 c x3 ds c.J 2 c
2
J2
Mf3
M f 3 M G 03 x1RG 02 b x2 ds b.J 3 b
2
J3
N Mf3 Mf2
11 a bx2 cx3 x2 x3
S J3 J2
31
S. MEZLINI - T. DAAMI – A. DOGUI 31
Calculateurs Temps Réel
Résistance des Matériaux
IV- Relations entre contraintes et efforts intérieurs
2. Problème de flexion simple
Même démarche, l’état des contraintes et le torseur des efforts intérieurs sont
de la forme :
11 12 13 R int T2 x 2 T3 x 3
12 0 0 int G
M int M f 2 x 2 M f3 x 3 G
13 0 0
32
S. MEZLINI - T. DAAMI – A. DOGUI 32
Calculateurs Temps Réel
Résistance des Matériaux
V Torseurs des petits déplacements et déformations
1 Torseur de petit déplacement ( torseur cinématique)
Dans la définition du travail Tf , le dual du torseur des efforts sera le
torseur cinématique U (s)
Sachant que T f R u M
( s )
R
effort s int
U ( s)
u ( s)
M
Hypothèse:
le mouvement rigidifiant la section S (s)
P
P S ( s) u( P) u( s) ( s) GP G
Rappel: En élasticité
1
A :1
ES
0 0 0 0 0
A: 1
0 0 0 0 0
n EJ 2 N
Equivalence en RDM: M
2 0 0
1
0 0 0 f 2
M f 3
3
EJ 3
Les coefficients sont identifiés n 1 M
0 0 0 0 0 t
par des essais simple de
2 GI T2
traction, flexion, torsion, en 1
s’appuyant sur le résulta des 3 0 0 0 0 0 T3
GS2
problèmes Saint Venant ou en
utilisant la méthode 0 0 0 0 0
1
GS3
énergétique.
Pour les section pleines, on montre que l’effet des efforts tranchants est
négligeable . Pour les poutre rectilignes et avec l’hypothèse de Bernoulli on
aura:
y
EIy ( x) M z ( x)
''
x
y (x)
Pour les section de profils minces les effets des efforts tranchant s sont très
importants .
37
S. MEZLINI - T. DAAMI – A. DOGUI 37
Calculateurs Temps Réel
Résistance des Matériaux
VI Loi de comportement et énergie de déformation
ds
T
T ds ( s) t 2 x 2
G
S La section n’est pas perpendiculaire à la
fibre moyenne après déformation
38
S. MEZLINI - T. DAAMI – A. DOGUI 38
Calculateurs Temps Réel
Résistance des Matériaux
VI méthode Energétique
Energie de déformation
1
En élasticité w
2 : dv
0
ext d ext l
N équations indépendantes
Cas 2: N-n=h
Dimensionnement
Torseur des efforts Torseur des efforts en terme d’efforts
extérieurs intérieurs
ext d s L.C.
U s D s
S. MEZLINI - T. DAAMI – A. DOGUI 43 43
Calculateurs Temps Réel
Méthode de résolution d’un problème RDM
2. Méthode énergétique
Théorème de Castigliano
Si une structure est soumise à une action A (une force ou un couple)
concentrée, alors la dérivée partielle de l’énergie de déformation par rapport à
cette action donne le déplacement de son point d’application P selon la
direction de A. w
P (deplacement suivant la direction de A) def
A( P)
Théorème de Menebrea: Pour un système hyperstatique, les inconnues
hyperstatiques prennent les valeurs qui minimisent l’énergie de déformation du système.
Ce théorème est un corollaire du théorème de Castigliano à savoir que les déplacements
relatifs aux actions de liaison sont nuls. wdef
liaison 0
Aaction liaison
On choisit h inconnues hyperstatiques / N-h inconnues s’écrivent en
fonction de h inconnues hyperstatiques
Démarche
s, h
hyperstatiques choisies hi
h
i
ext d i
w
0 Energie de déformation en
hi hi
h équations fonction des inconnues
supplémentaires hyperstatiques choisies
1 B
La structures est hyperstatique d’ordre 1, pour
1 soulever l’hyperstaticité on applique Manebréa
A
wdef
0
On introduit 1 , 2 , 3 Forces virtuelles M A
wdef wdef
1 2
1 0
1 0
1 , 2 , 3 caractérisant la déformée
1 1
47
S. MEZLINI - T. DAAMI – A. DOGUI 47
Calculateurs Temps Réel
Méthode de résolution d’un problème RDM
Dimensionnement élastique
Von Mises
éqV x2 3 2
Démonstration: 1
D
tr .1
3
x xy xz 2
3x xy xz
xy 0 0
xz 0 0 1
xy x 0
3
Tenseurs de contrainte dans une
poutre. Voir problème de Saint Venant
1
x
pour déterminer les termes de ce xz 0
tensueur 3
9 9
6 2
6 2
2 2
x 2 xy xz x 2 x 3 2
9
2 2
9
2
3
3 D D
éqV : x 3 2
2
Démonstration: x xy xz
Tenseurs de contrainte
xy 0 0
dans une poutre
xz 0 0
1 0
x xy xz
x x2 4 2
.1 0 xy 0 0 2
xz 0 2
x x2 4 2
3
2
éqT 4
2
x
2
53
S. MEZLINI - T. DAAMI – A. DOGUI 53
Calculateurs Temps Réel
Méthode de résolution d’un problème RDM
r/d
54
S. MEZLINI - T. DAAMI – A. DOGUI 54