Mcanique des structures en matriaux composites par les mthodes 2013

analytiques et des lments finis - Statique

Chapitre2 lment de barre

Une barre est un lment soumis seulement deux forces axiales; il ny a ni flexion ni
cisaillement. Le seul champ considrer est le dplacement axial de la barre u(x). Ces
lments sont principalement utiliss dans les structures comportant des membrures
forces axiales comme treillis et colonnes.

2.1 Treillis
Structure comportant des membrures deux forces axiales. Les forces sappliquent en
gnral aux joints.

Charge

Barre deux forces

axiales

Figure 2.1 Treillis sous charge verticale

2.1.1 Formulation des lments de barre dans le systme de repre global

La force F agissant sur une barre de section A et dune longueur L, illustre la figure
2.3, est dfinie comme :

AE
=F ( )L
L
La barre peut se modlise comme un ressort dont la rigidit est :
AE
k=
L

1
Mcanique des structures en matriaux composites par les mthodes 2013
analytiques et des lments finis - Statique

Treillis simple: Treillis hyperstatique:

3 quations dquilibre 3 quations dquilibre
3 inconnus 4 inconnus

Figure 2.2 Exemple de treillis simple et hyperstatique

L
L

Figure 2.3 Barre soumise une charge axiale

Le treillis prsents la figure 2.4 est modlis par 6 lments de barre. Le modle
comporte 5 nuds et 6 lments.

2
Mcanique des structures en matriaux composites par les mthodes 2013
analytiques et des lments finis - Statique

Charge Charge

(3) (6)

4 5
3

(2) (4)
(5)

(1)

1 2

Figure 2.4 Treillis en charge

Le repre global (X,Y) est utilis dans la figure 2.5 afin de considrer la localisation des
joints, lorientation des membrures ainsi que les conditions aux limites et les charges
appliques. Le repre local ou lmentaire (x,y) se trouve sur la membrure pour la
description de son comportement sous chargement individuel. Les relations entre les deux
systmes de coordonns sont :

=
Uix u ix cos u iy sin
=
Uiy u ix sin + u iy cos
=
U jx u jx cos u jy sin
=
U jy u jx sin + u jy cos

ou sous forme matricielle :

{U} = [T ]{u}
o :
U ix cos sin 0 0 u ix
sin cos
U iy 0 0 u iy
={U} = , [T ] = , et {u}
U jx 0 0 cos sin u jx
U u
0 0 sin cos
jy jy

{U} et {u} reprsentent les dplacements aux nuds i et j dans les repres global et local
respectivement. [T] est la matrice permettant la transformation du repre local au repre
global.

3
Mcanique des structures en matriaux composites par les mthodes 2013
analytiques et des lments finis - Statique

Ujy
X

ujx
y
uj
Ujx Fjy
j

fjx
fjy
j
Fjx
x
y

x
y
Uiy
uix
y

ui

Uix Fiy

fix
i
fiy

i Fix

Figure 2.5 Relations entre les deux systmes de coordonn global et local

De faon analogue, on obtient pour les forces nodales (figure 2.6) :

=
Fix fix cos fiy sin
=
Fiy fix sin + fiy cos
=
Fjx f jx cos f jy sin
=
Fjy f jx sin + f jy cos
Forme matricielle :

{F} = [T ]{f }
o
Fix f ix

Fiy f iy
={F} = et {f }
Fjx f jx
F f
jy jy

4
Mcanique des structures en matriaux composites par les mthodes 2013
analytiques et des lments finis - Statique

ujx
Y
fjx=k(uix-ujx)

X ujx
j
fjx=k(ujx-uix)
(a)
uix ou j
y x

uix
i (b)
fix=k(uix-ujx)

i
fix=k(uix-ujx)
Figure 2.6 DCL dune membrure deux forces

Notons que la barre est une membrure 2 forces axiales seulement. Par consquent, les
forces (f iy , f jy ) et les dplacements nodaux (u iy , u jy ) sont nulles. quilibre des forces
internes dans la barre illutes la figure 2.6(b) donne :

f ix k 0 k 0 u ix

f iy 0 0 0 0 u iy
=
f jx k 0 k 0 u jx

0 0 0 u jy
f jy 0
AE
o : k=
L
Forme matricielle :

5
Mcanique des structures en matriaux composites par les mthodes 2013
analytiques et des lments finis - Statique

{f } = [ K ]{u}
Or

{f } {
u}

[T ] {F} = [ K ][T ] {U}
1 1

o :
cos sin 0 0
sin cos 0 0
[T ] =
1

0 0 cos sin

0 0 sin cos

Multiplication des membres de lquation prcdente par [T] et simplification nous

donnent finalement la relation entre les forces et les dplacements nodaux dans le repre
globale :

{F} = [T ][ K ][ T ] {U}
1

Fix cos 2 sin cos cos 2 sin cos U ix

F
iy sin cos sin 2 sin cos sin 2 U iy
=k
Fjx cos 2 sin cos cos 2 sin cos U jx
F

jy sin cos sin 2
sin cos sin 2
U jy

La matrice de rigidit lmentaire dans le repre global est donc :

cos 2 sin cos cos 2 sin cos

sin cos sin
2
sin cos sin 2
K (e) =k
cos 2 sin cos cos 2 sin cos

sin cos sin 2 sin cos sin 2

Lexemple suivant prsente lapplication numrique du processus au complet en

commenant par le prtraitement puis la solution du systme dquations et en terminant
avec le post-traitement.

Exemple 2.1 Considrons le treillis prsent la figure 2.7 dont les caractristiques
sont :

6
Mcanique des structures en matriaux composites par les mthodes 2013
analytiques et des lments finis - Statique

A=5000mm2, E=13100MPa, L1 =L 3 =L 4 =L6 =900mm, L2 =L 5 =1 270.315mm

1 = 3 = 6 = 0o
2 = 135o
4 =90o
5 =45o

Calculez les dplacements aux joints, les contraintes dans les membrures et les ractions.

Solution

2000N 2000N

(3) (6)

4 5
3

(2) (4)
(5)

(1) 45o

1 2
Figure 2.7 Treillis de bois en charge

Prtraitement

Discrtisation

Tableau 2.1 Proprits des lments dans le modle du treillis la figure 2.7
lment Nud i Nud j o
1 1 2 0
2 2 3 135
3 3 4 0
4 2 4 90
5 2 5 45
6 4 5 0

Comportement des lments

k 1 = k 3 = k 4 = k 6 = (AE/L) = (5000x13100/900) = 72 778 N/mm

k 2 = k 5 = 51 562 N/mm

7
Mcanique des structures en matriaux composites par les mthodes 2013
analytiques et des lments finis - Statique

Matrices de rigidit lmentaires

X
y
i j
x

1 lment 1 2

3 lment 3 4
4 lment 6 5

Figure 2.8 lments (1), (3) et (6) dans les repres local et global

lment (1) 1 = 0 et k 1 = 72 778 N/mm :

La matrice de rigidit lmentaire dans le repre global est :

72778 0 -72778 0 U1X

0 0 0 0 U1Y
[ KEBG1] =
-72778 0 72778 0 U 2X

0 0 0 0 U 2Y

Matrice de rigidit lmentaire dans le repre global avec ses lments placs en position
pour lassemblage la matrice globale du treillis :

[KEBGA1] = [A1]T[KEG1][A1]

8
Mcanique des structures en matriaux composites par les mthodes 2013
analytiques et des lments finis - Statique

72778 0 -72778 0 0 0 0 0 0 0 U1X

0 0 0 0 0 0 0 0 0 0 U1Y

-72778 0 72778 0 0 0 0 0 0 0 U 2X

0 0 0 0 0 0 0 0 0 0 U 2Y
0 0 0 0 0 0 0 0 0 0 U 3X
[ KEBGA1] =
0 0 0 0 0 0 0 0 0 0 U 3Y
0 0 0 0 0 0 0 0 0 0 U 4X

0 0 0 0 0 0 0 0 0 0 U 4Y
0 0 0 0 0 0 0 0 0 0 U 5X

0 0 0 0 0 0 0 0 0 0 U 5Y

o la matrice de slection est :

1 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0
[ A1] =
0 0 1 0 0 0 0 0 0 0

0 0 0 0 0 0 0 0 0 0

lment (3) 3 = 0 et k 3 = 72 778 N/mm :

72778 0 -72778 0 U 3X
0 0 0 0 U 3Y
[ KEBG3] =
-72778 0 72778 0 U 4X

0 0 0 0 U 4Y

0 0 0 0 1 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0
[ A3] =
0 0 0 0 0 0 1 0 0 0

0 0 0 0 0 0 0 0 0 0

9
Mcanique des structures en matriaux composites par les mthodes 2013
analytiques et des lments finis - Statique

0 0 0 0 0 0 0 0 0 0 U1X
0 0 0 0 0 0 0 0 0 0 U1Y

0 0 0 0 0 0 0 0 0 0 U 2X

0 0 0 0 0 0 0 0 0 0 U 2Y
0 0 0 0 72778 0 -72778 0 0 0 U 3X
[ KEBGA3] =
0 0 0 0 0 0 0 0 0 0 U 3Y
0 0 0 0 -72778 0 72778 0 0 0 U 4X

0 0 0 0 0 0 0 0 0 0 U 4Y
0 0 0 0 0 0 0 0 0 0 U 5X

0 0 0 0 0 0 0 0 0 0 U 5Y

lment (6) 6 = 0 et k 6 = 72 778 N/mm :

72778 0 -72778 0 U 4X
0 0 0 0 U 4Y
[ KEBG6] =
-72778 0 72778 0 U 5X

0 0 0 0 U 5Y

0 0 0 0 0 0 1 0 0 0
0 0 0 0 0 0 0 0 0 0
[ A6] =
0 0 0 0 0 0 0 0 1 0

0 0 0 0 0 0 0 0 0 0

10
Mcanique des structures en matriaux composites par les mthodes 2013
analytiques et des lments finis - Statique

0 0 0 0 0 0 0 0 0 0 U1X
0 0 0 0 0 0 0 0 0 0 U1Y

0 0 0 0 0 0 0 0 0 0 U 2X

0 0 0 0 0 0 0 0 0 0 U 2Y
0 0 0 0 0 0 0 0 0 0 U 3X
[ KEBGA6] =
0 0 0 0 0 0 0 0 0 0 U 3Y
0 0 0 0 0 0 72778 0 -72778 0 U 4X

0 0 0 0 0 0 0 0 0 0 U 4Y
0 0 0 0 0 0 -72778 0 72778 0 U 5X

0 0 0 0 0 0 0 0 0 0 U 5Y

lment 2 : 2 = 135o et k 2 = 51 562 N/mm :

j
x
Y
135o
x
X y 2

lment 2
i
Figure 2.9 lment (2) dans les repres local et global

0 0 1 0 0 0 0 0 0 0
0 0 0 1 0 0 0 0 0 0
[ A2] =
0 0 0 0 1 0 0 0 0 0

0 0 0 0 0 1 0 0 0 0

11
Mcanique des structures en matriaux composites par les mthodes 2013
analytiques et des lments finis - Statique

0 0 0 0 0 0 0 0 0 0 U1X
0 0 0 0 0 0 0 0 0 0 U1Y

0 0 2.578 -2.578 -2.578 2.578 0 0 0 0 U 2X

0 0 -2.578 2.578 2.578 -2.578 0 0 0 0 U 2Y
0 0 -2.578 2.578 2.578 -2.578 0 0 0 0 U 3X
[ KEBGA2] = 1.0e+004
0 0 2.578 -2.578 -2.578 2.578 0 0 0 0 U 3Y
0 0 0 0 0 0 0 0 0 0 U 4X

0 0 0 0 0 0 0 0 0 0 U 4Y
0 0 0 0 0 0 0 0 0 0 U 5X

0 0 0 0 0 0 0 0 0 0 U 5Y

lment 4 : 4 = 90o et k 4 = 72 778 N/mm :

x Y j
4
o
90

X x
y

lment 4

i 2

Figure 2.10 lment (4) dans les repres local et global

0.0000 0.0000 -0.0000 -0.0000 U 2X

0.0000 7.2778 -0.0000 -7.2778 U 2Y
[ KEBG4] = 1.0e+004
-0.0000 -0.0000 0.0000 0.0000 U 4X

-0.0000 -7.2778 0.0000 7.2778 U 4Y

12
Mcanique des structures en matriaux composites par les mthodes 2013
analytiques et des lments finis - Statique

0 0 0 0 0 0 0 0 0 0
0 0 0 1 0 0 0 0 0 0
[ A4] =
0 0 0 0 0 0 0 0 0 0

0 0 0 0 0 0 0 1 0 0

0 0 0 0 0 0 0 0 0 0 U1X
0 0 0 0 0 0 0 0 0 0 U1Y

0 0 0 0 0 0 0 0 0 0 U 2X

0 0 0 72778 0 0 0 -72778 0 0 U 2Y
0 0 0 0 0 0 0 0 0 0 U 3X
[ KEBGA4] =
0 0 0 0 0 0 0 0 0 0 U 3Y
0 0 0 0 0 0 0 0 0 0 U 4X

0 0 0 -72778 0 0 0 72778 0 0 U 4Y
0 0 0 0 0 0 0 0 0 0 U 5X

0 0 0 0 0 0 0 0 0 0 U 5Y

lment 5 : 5 = 45o et k 5 = 51 562 N/mm :

Y
x

45o
5

X
x
y
i

lment 5
2

Figure 2.11 lment (5) dans les repres local et global

13
Mcanique des structures en matriaux composites par les mthodes 2013
analytiques et des lments finis - Statique

2.5781 2.5781 -2.5781 -2.5781 U 2X

2.5781 2.5781 -2.5781 -2.5781 U 2Y
[ KEBG5] = 1.0e+004
-2.5781 -2.5781 2.5781 2.5781 U 5X

-2.5781 -2.5781 2.5781 2.5781 U 5Y

0 0 1 0 0 0 0 0 0 0
0 0 0 1 0 0 0 0 0 0
[ A5] =
0 0 0 0 0 0 0 0 1 0

0 0 0 0 0 0 0 0 0 1

0 0 0 0 0 0 0 0 0 0 U1X
0 0 0 0 0 0 0 0 0 0 U1Y

0 0 2.578 2.578 0 0 0 0 -2.578 -2.578 U 2X

0 0 2.578 2.578 0 0 0 0 -2.578 -2.578 U 2Y
0 0 0 0 0 0 0 0 0 0 U 3X
[ KEBGA5] = 1.0e+004
0 0 0 0 0 0 0 0 0 0 U 3Y
0 0 0 0 0 0 0 0 0 0 U 4X

0 0 0 0 0 0 0 0 0 0 U 4Y
0 0 -2.578 -2.578 0 0 0 0 2.578 2.578 U 5X

0 0 -2.578 -2.578 0 0 0 0 2.578 2.578 U 5Y

Assemblage de la matrice de rigidit total du modle

La matrice de rigidit du treillis dans le repre global est :

[KG]=[KEBGA1]+[KEBGA2]+[KEBGA3]+[KEBGA4]+[KEBGA5]+[KEBGA6]

14
Mcanique des structures en matriaux composites par les mthodes 2013
analytiques et des lments finis - Statique

0.7278 0 -0.7278 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0

-0.7278 0 1.2434 0 -0.2578 0.2578 0 -0.2578 -0.2578
0

0 0 0 1.2434 0.2578 -0.2578 0 -0.7278 -0.2578 -0.2578
0 0 -0.2578 0.2578 0.9856 -0.2578 -0.7278 0 0 0
[ KG ] = 1.0e+5
0 0 0.2578 -0.2578 -0.2578 0.2578 0 0 0 0
0 0 0 0 -0.7278 0 1.4556 0 -0.7278 0

0 0 0 -0.7278 0 0 0 0.7278 0 0
0 0 -0.2578 -0.2578 0 0 -0.7278 0 0.9856 0.2578

0 0 -0.2578 -0.2578 0 0 0 0 0.2578 0.2578

Application des conditions aux limites

U 1X =U 1Y =U 3X =U 3Y =0 et F 4Y =F 5Y =2000N impliquent le nouveau systme dquation

suivant :

1 0 0 0 0 0 0 0 0 0 U1X 0
0 1 0 0 0 0 0 0 0 U
0 1Y 0

-0.73 0 1.24 0 -0.26 0.26 0 0 -0.26 -0.26 U 2X 0

0 0 0 1.24 0.26 -0.26 0 -0.73 -0.26 -0.26 U 2Y 0
0 0 0 0 1 0 0 0 0 0 U 3X 0
1.0e+5 =
0 0 0 0 0 1 0 0 0 0 U 3Y 0
0 0 0 0 -0.73 0 1.46 0 -0.73 0 U 4X 0

0 0 0 -0.73 0 0 0 0.73 0 0 U 4Y 2000
0 0 -0.26 -0.26 0 0 -0.73 0 0.99 0.26 U 5X 0

0 0 -0.26 -0.26 0 0 0 0 0.26 0.26 U 5Y 2000

Sachant que U 1X =U 1Y =U 3X =U 3Y =0, on peut liminer les premire,deuxime, cinquime

et sixime lignes ainsi que les colonnes correspondant ces dplacements dans la matrice
de rigidit.

15
Mcanique des structures en matriaux composites par les mthodes 2013
analytiques et des lments finis - Statique

Rsolution du systme dquations :

1 0 0 0 0 0 0 0 0 0 U1X 0

0 1 0 0 0 0 0 0 0 0 U1Y 0

-0.73 0 1.24 0 -0.26 0.26 0 0 -0.26 -0.26 U 2X 0

0 0 0 1.24 0.26 -0.26 0 -0.73 -0.26 -0.26 U 2Y 0
0 0 0 0 1 0 0 0 0 0 U 3X 0
1.0e+5 =
0 0 0 0 0 1 0 0 0 0 U 3Y 0
0 0 0 0 -0.73 0 1.46 0 -0.73 0 U 4X 0

0 0 0 -0.73 0 0 0 0.73 0 0 U 4Y 2000
0 0 -0.26 -0.26 0 0 -0.73 0 0.99 0.26 U 5X 0

0 0 -0.26 -0.26 0 0 0 0 0.26 0.26 U 5Y 2000

Rsolution de ce systme dquations rduit nous donne les dplacements

correspondants :

U 2X -0.0868
U -0.2503
2y
U 4X 0.0274
{UREDUCED
= } = mm
U 4Y -0.2777
U 0.0548
5X
U 5Y -0.4688

La matrice des dplacements nodaux au complet est :

U1X 0
U 0
1Y
U 2X -0.0868

U 2Y -0.2503
U 3X 0
{UTOTAL
= } = mm
U 3Y 0
U 4X 0.0274

U 4Y -0.2777
U 0.0548
5X
U 5Y -0.4688

16
Mcanique des structures en matriaux composites par les mthodes 2013
analytiques et des lments finis - Statique

Post-Traitement

La matrice des ractions est :

{R}=[KG]*{UTOTAL}-{F}

R1X 6.3171
R 0
1Y
R 2X -0.1194

R 2Y -0.2385
R 3X
{R} =
= 3 -6.2093
10
R 3Y 4.2152
R 4X 0

R 4Y 0.0059
R 0.0116
5X
R 5Y 0.0174

Les dplacements nodaux et les forces internes dans chaque lment se calculent laide
des quations :

[u] = [T]-1 [U]

Et
f ix k 0 k 0 u ix
f
iy 0 0 0 0 u iy
=
f jx k0 k 0 u jx
f 0
0 0 0 u jx
jx
AE
o : k=
L
Les contraintes lmentaires se dterminent par lexpression :

f k(u ix u jx ) AE(u ix u jx ) u u jx
= = = = E( ix )
A A AL L

17
Mcanique des structures en matriaux composites par les mthodes 2013
analytiques et des lments finis - Statique

2.2 Colonne
Exemple 2.2 Soit une colonne en acier (E=199 948MPa) qui supporte 4 planchers telle
que prsente la figure 2.12. Sachant que la colonne a une section de 25613 mm2,
dterminez :

les dplacements axiaux au niveau de chaque plancher;

les contraintes dans la colonne dans chaque portion entre les planchers;
les dplacements aux niveaux A et B.

Solution

k= A*E/L = 25613 x 199948/4572= 1 120 137N/mm

K=[k -k;-k k]

1.1201 -1.1201
K 1.0e+6 *

-1.1201 1.1201

Figure 2.12 Colonne en charge. Adapte de rf [2].

18
Mcanique des structures en matriaux composites par les mthodes 2013
analytiques et des lments finis - Statique

KEG1=KEG2=KEG3=KEG4=K
lment (1) :

Matrice de slection :

A1=[1 0 0 0 0;0 1 0 0 0]

1 0 0 0 0
A1
0 1 0 0 0

KEGA1=A1'*KEG1*A1

1.1201 -1.1201 0 0 0

-1.1201 1.1201 0 0 0

KEGA1 1.0e+6 * 0 0 0 0 0
0 0 0 0 0

0
0 0 0 0

lment (2) :

A2=[0 1 0 0 0;0 0 1 0 0]

0 1 0 0 0
A2
0 0 1 0 0

KEGA2=A2'*KEG2*A2

0 0 0 0 0

0 1.1201 -1.1201 0 0

KEGA2 1.0e+6 * 0 -1.1201 1.1201 0 0
0 0 0 0 0

0
0 0 0 0

lment 3 :
A3=[0 0 1 0 0;0 0 0 1 0]

19
Mcanique des structures en matriaux composites par les mthodes 2013
analytiques et des lments finis - Statique

0 0 1 0 0
A3
0 0 0 1 0

0 0 0 0 0

0 0 0 0 0

KEGA3 1.0e+6 * 0 0 1.1201 -1.1201 0
0 0 -1.1201 1.1201 0

0 0
0 0 0

lment (4) :
A4=[0 0 0 1 0;0 0 0 0 1]

0 0 0 1 0
A4
0 0 0 0 1

KEGA4=A4'*KEG4*A4

0 0 0 0 0

0 0 0 0 0

KEGA4 1.0e+6 * 0 0 0 0
0
0 0 0 1.1201 -1.1201

0 0 -1.1201 1.1201
0

Matrice de rigidit globale :

KG=KEGA1+KEGA2+KEGA3+KEGA4

1.1201 -1.1201 0 0 0

-1.1201 2.2403 -1.1201 0 0

KG 1.0e+6 * 0 -1.1201 2.2403 -1.1201 0
0 0 -1.1201 2.2403 -1.1201

0 -1.1201 1.1201
0 0

Matrice des charges :

FG=-[0;222411;222411;222411;266893]

20
Mcanique des structures en matriaux composites par les mthodes 2013
analytiques et des lments finis - Statique

0

-222411

FG -222411 N

-222411

-266893

Application des conditions aux limites :

U 1 =0

1.1201 -1.1201 0 0 0 U1 0

-1.1201 2.2403 -1.1201 0 0 U 2 -222411

1.0e+6 * 0 -1.1201 2.2403 -1.1201 0 U 3 -222411


0 0 -1.1201 2.2403 -1.1201 U 4 -222411

0 0 0 -1.1201 1.1201 U 5 -266893

Rsolution de ce systme dquations rduit nous donne les dplacements

correspondants.
Le vecteur des dplacements est :

U1 u1 0


U 2

u 1 -0.8339


U 3 u1 -1.4693 mm

U 4 u -1.9061

1


U 5 u1 -2.1444

Les contraintes dans les lments (sections entre les planchers) sont :

E(u j u i )
(e)
L

function y=SIGMA(E,L,uj,ui)
y=E*(uj-ui)/L;

21
Mcanique des structures en matriaux composites par les mthodes 2013
analytiques et des lments finis - Statique

SIGMA1=SIGMA(199948,4572,u2,u1)
SIGMA1 = -36.4708MPa
SIGMA2=SIGMA(199948,4572,u3,u2)
SIGMA2 = -27.7873MPa
SIGMA3=SIGMA(199948,4572,u4,u3)
SIGMA3 = -19.1026 MPa
SIGMA4=SIGMA(199948,4572,u5,u4)
SIGMA4 = -10.4202MPa

Dplacement au niveau A(u j =u 2 =-0.8339mm, u i = u 1 = 0mm)

Pour un lment linaire les fonctions de forme dans le repre local sont :
y y
N i 1 ; Nj
L L
(Voir figure 2.13 pour la relation entre les repres local et global dun lment de linaire
de barre.

Le point A appartient llment (1) do son dplacement est interpol avec les
dplacements aux nuds 1 et 2 :
y y
u (A) N1u1 N 2 u 2 (1 A )u1 ( A )u 2
L L
3048 3048
u (A) (1 )(0)+( )(-0.8339)
4572 4572
u (A) 0.5559mm

Dplacement au niveau B(u j = u 5 = -2.1444mm, u i = u 4 = -1.9061mm)

yB y
u (B) N 4 u 4 N 5 u 5 (1)u 4 ( B )u 5
L L
4572 2438 4572 2438
u (A) (1 )(-1.9061)+( )(-2.1444)
4572 4572
u (A) 2.01732mm

22
Mcanique des structures en matriaux composites par les mthodes 2013
analytiques et des lments finis - Statique

Figure 2.13 Relation entre les repres local et global dun lment de barre. Adapt de
rf. [2]
Rfrences
[1] Comprendre les lments finis - Principes, formulations et exercices
corrigs, Alain Chateauneuf, Ellipses, ISBN 978-2-7298-5430-0
[2] Finite Element Analysis Theory and Application with ANSYS, 3rd
Edition, Saeed Moaveni, Pearson Prentice Hall, ISBN 978-0-13-189080-0
[3] A first course in the Finite Element Method using Algor, Dany Logan,
PWS publishing Co, ISBN 0-534-94692-5

Exercices
2.1 Les membrures du treillis prsent la figure 1 ont une section de 1500mm2.
Sachant que leur module dlasticit est de 70GPa, calculez :

a. le dplacement du joint A;
b. les contraintes dans les membrures;
c. les ractions.

23
Mcanique des structures en matriaux composites par les mthodes 2013
analytiques et des lments finis - Statique

900N

mm
915 60o 60o

Figure 1

2.2 La figure 2 illustre un treillis dont les membrures ont une section de 1500mm2.
Sachant que leur module dlasticit E = 70GPa, calculez :

a. les dplacements des joints;

b. les contraintes dans les membrures;
c. les ractions.

1500N
3 4
1000mm

1 1500mm 2

24
Mcanique des structures en matriaux composites par les mthodes 2013
analytiques et des lments finis - Statique

Figure 2

2.3 Le systme illustr la figure 3 comporte des membrures en acier dont le module
dlasticit E = 200GPa. Dterminez :

a. les dplacements des points E et F;

b. les contraintes dans les membrures.

125mm 100mm

A= 322mm2 A= 515mm2
D E
F
15kN

300mm 250mm

Figure 3

2.4 Le poteau prsent la figure 4 est fait de tube en acier dont le module dlasticit
est de 210GPa. En ngligeant leffet du vent, calculez :

a. le dplacement du point A au sommet du poteau;

b. les contraintes dans les sections du poteau.

25
Mcanique des structures en matriaux composites par les mthodes 2013
analytiques et des lments finis - Statique

450N

1500mm
670N
A=484mm2

1500mm
900N
A=1390mm2

3000mm
A=1900mm2

26