Chapitre 8 h2014 PDF
Chapitre 8 h2014 PDF
Chapitre 8 h2014 PDF
2014
x
z
O
Axe
neutre
La flche w dune poutre en flexion est obtenue partir de la relation avec le moment
appliqu M et la rigidit en flexion J=EI :
d2w
(8.1)
M
dx 2
o E est le module dlasticit du matriau et I dsigne le second moment de la section de
la poutre.
EI
1
d2w
2
x
dx
(8.2)
2014
M EI x J x
(8.3)
Cette relation est valide que pour les matriaux isotropes. Dans le cas des poutres en
stratifi de composites (Figure 8.2) le module dlasticit de chacun des plis varie selon
le matriau et lorientation des renforts. Par consquent, afin de calculer la flche dune
poutre en stratifi de matriaux composites, il est ncessaire dtablir sa rigidit
apparente J a =E a I.
L
q(x)
b
h
y
z
8.1
D
M xy 16
D12
D 22
D 26
0
D16 x
0
D 26 y
D66 0
xy
(8.4)
0x
2w 0
,
x 2
0y
2w 0
,
y 2
0xy 2
2w 0
xy
(8.5)
2014
0 D ' D '
'
D16
x 11
12
M x
0 '
'
'
(8.6)
y D12 D 22 D 26 M y
0 '
'
'
xy D16 D 26 D66 M xy
'
o Dij sont les lments de la matrice inverse de la matrice de rigidit en flexion-torsion
[D] du stratifi de composites.
Il y a deux cas distincts considrer : les poutres troites et les poutres larges. La
distinction est base sur le rapport b/h de la section (Voir figure 8.1). Le coefficient de
Poisson implique une distorsion transversale de la section (Figure 8.3). Cet effet est
localis seulement aux bords dune poutre large qui comporte comme une plaque.
Poutre troite
Poutre large
Figure 8.3 Distorsion de de la section : poutre aile troite et celle aile large
8.1.1
M
0 ),
b
alors :
M y M xy 0
(8.7)
(8.8)
2w 0
'
Mx
D16
xy
Remarque :
Les courbures y et xy sont fonctions du M x selon la relation (8.6). Elles dpendent alors
de la variable y. Cependant, leffet de la flexion et de la torsion induite est ngligeable
lorsque le rapport L/b est suffisamment lev. Cest pourquoi, la thorie des poutres
considre que la flche nest fonction que de x pour une poutre dont le rapport L/b est
lev :
2014
w 0 w 0 (x)
(8.9)
d2w 0
'
' M
0x D11
M x D11
2
dx
b
(8.10)
do :
bM x
1
M
'
D11
Mx
3
bh
x E fx I yy
E fx
12
3
bh
o:
I yy =
12
Le module de rigidit en flexion devient :
0x
1
x
et que
E fx
12
'
D11
h3
(8.11)
(8.12)
(8.13)
d2w 0
M
=
2
dx
E f I yy
(8.14)
=
J a E=
fx I
Sachant que :
2014
0x
x
z 0
y
k
y
xy
k
xy
(8.14)
=
{ (xy) }k Q {(xy) }k
k
(8.15)
Px
2
0x
d2 w 0
M
Px
=
2
E fx I
2E fx I
dx
Application des conditions aux limites :
x= 0 :
x= L/2 :
L
2
0x
L
2
(a)
M = w0 = 0
dw 0 /dx = 0
=
w0
x 3 ( ) 2
48E fx I
L
La flche au centre de la poutre o x = L/2 est :
PL3
PL3 '
=
w0 =
D11
48E fx I 48b
et les contraintes dans le ke pli sont :
0x
x
0
_
Q
z
y
y=
k
k
0
xy
xy k
2014
Exemple 8.2
Soit une poutre en appuis simples, faite de stratifi [0/60] S du composite carbone/poxy
AS3501 (v f =0.65), qui est soumise une charge rpartie q(x)= 0.2N/mm. Sachant que
les dimensions de la poutre sont L= 100mm et b = 5mm, et que les plis ont la mme
paisseur de 0.150mm, calculez :
1. La flche maximale;
2. Les contraintes au niveau suprieur du 2e pli.
Solution
Proprits du matriau : E 1 =138000MPa, E 2 =9000MPa, G 12 =6.900MPa, nu 12 =0.3
Stratifi symtrique :
Pli#
1
Direction
1 =0
2 =60
3 =-60
Hauteur
Z0=-0.450mm
Z1=-0.300mm
Z1=-0.300mm
Z2=-0.150mm
Z2=-0.150mm
Z3=0.000mm
Moment de flexion :
qL2 0.2 1002
=
M =
= 250N mm
8
8
M 200
M=
=
= 50N
x
b
5
1
'
D11
= 0.0001631393(
)
MPa mm3
b
5
=
J a E=
=
= 30648.66118MPa mm 4
fx I
'
D11 0.0001631393
5qL4
5qL4
5 0.2 1004
w
=
= =
= 8.4968mm
0max
384E fx I 384J a 384 30648.66118
Le vecteur de chargement est : [0;0;0;50;0;0]
et les dformations et les contraintes dans le ke pli sont respectivement :
2014
x
x
z 0
y
k
y
xy
k
xy
{ }
(xy) k
_
= Q { (xy) }
k
k
x 29.9613
y =
3.8738 MPa
6.3084
xy
8.1.2
Dans ce cas o le rapport b/h est lev, les courbures suivantes sont considres nulles :
0y =0xy =0
Lquation (8.4) devient :
M x D11
M D
y 12
D
M xy 16
D12
D 22
D 26
D16 0x
D 26 0
D66 0
do
=
M bM
=
bD110x
x
et la rigidit apparente de la poutre est :
=
E fx
12D11
=
et J a bD11
h3
(8.16)
(8.17)
Exemple 8.3
Soit une poutre en appuis simples, faite de stratifi [0/90/-30/30] S du composite
carbone/poxy AS4/3501-6 (E 1 =131GPa, E 2 =11.2GPa. G 12 =6.55GPa, 12 =0.28), qui
est soumise une charge rpartie q(x)= 0.38N/mm. Sachant que les dimensions de la
poutre sont L= 254mm et b = 12.7mm et que les plis ont la mme paisseur de 0.4mm,
calculez :
2014
La flche maximale;
Les contraintes au niveau suprieur du 2e pli.
Solution
Pli#
1
Direction
1 =0
2 =90
3 =-30
4 =30
Hauteur
Z0=-1.6mm
Z1=-1.2mm
Z1=-1.2mm
Z2=-0.8mm
Z2=-0.8mm
Z3=-0.4mm
Z3=-0.4mm
Z4=0mm
3
[D] MPa-mm = 1.0e + 005 * 0.1570 1.3177 0.0362
0.0975 0.0362 0.2497
=
M
0 0156
0.04170.0045.
Poutre troite :
=
J a E=
fx I
w
=
0max
b
12
=
= 3043183.8937948MPa mm 4
'
D11 4.17326078318695x10-6
5qL4
5qL4
5 0.38 2544
= =
=
384E fx I 384J a 384 3043183.8937948
6.76751218392041mm
2014
0.1014
0.0380
x -1.3303
y =
1.3631 MPa
-0.2987
xy
Poutre large :
J a bD
12.7x244984.810784553
=
=
= 3111307.09696382MPa-mm 4
11
w
=
0max
5qL4
5qL4
5 0.2 1004
= =
= 6.61933503743946mm
384E fx I 384J a 384 3111307.09696382
w 0(troite) w 0(larg e)
w 0(troite)
8.1.3
6.76751218392041mm 6.61933503743946mm
)% (=
)x100 2.19%
6.76751218392041mm
Poutre en sandwich
Exemple 8.4
Soit une poutre en sandwich dont les peaux se constituent de stratifi [0 2 /45 2 /-45 2] S des
plis dpaisseur t k = 0.15mm, du composite carbone /poxy AS-3501 (v f =0.65,
E 1 =138GPa, E 2 =9 GPa,G 12 =6.9GPa, nu 12 =0.3), qui est en appuis simples. Le noyau
dune paisseur de 6mm est en matriau isotrope (E=7MPa, G 12 = 2.9GPa, nu 12 =0.2).
Sachant que les dimensions de la poutre sont L=508mm, b=50.8mm et que la charge
rpartie applique transversalement la poutre est q=0.14N/mm, calculez :
1. La flche maximale;
2. Les contraintes au niveau suprieur du 2e pli.
Solution
Peau: AS3501 [02/452/-452]s
Noyau: matriau Isotrope
Peau: AS3501 [02/452/-452]s
x
z
2014
Pli#
1
Matriau
AS3501
Direction (o)
1 =0
2
3
4
5
6
7
8
9
AS3501
AS3501
AS3501
AS3501
AS3501
Matriau isotrope
Matriau isotrope
Matriau isotrope
2 =45
3 =-45
4 =-45
5 =45
6 =0
7 =0
8 =0
9 =0
Hauteur (mm)
Z 0 =-4.80
Z 1 =-4.50
Z 2 =-4.20
Z 3 =-3.90
Z 4 =-3.60
Z 5 =-3.30
Z 6 =-3.00
Z 7 =-2.00
Z 8 =-1.00
Z 9 =0.00
Programme dixhuit_plis_sym_v3
AS3501 = matriau #1
Matriau isotrope = matriau #2
qL2
8
M qL2 0.14 5082
M=
=
=
= 88.9N
x
b
8b
8 50.8
M=
'
D11
= 0.2859e 6
b
50.8
J a E=
177684505MPa mm 4
=
=
=
fx I
6
'
D11 0.2859e
w
=
0max
5qL4
5qL4
5 0.14 5084
= =
= 0.68324mm
384E fx I 384J a 384 177684505
xy
8.2
Il est noter que lme de la poutre en stratifi de composite illustre la figure 8.4 est
constitue de plis dont laxe x sont orientes perpendiculaire laxe neutre de la poutre.
10
2014
hme
tme
baile
Figure 8.4 Orientation des fibres dans les plis dune poutre en I faite de stratifi de
composite
xme
tme
(a)
(Nx)me
(b)
yme
Figure 8.5 (a) Repre global pour les ailes et lme; (b) repre local de lme
8.2.1 me
Il est remarquer que les plis de lme sont parallles avec laxe z du repre global de la
poutre. Par consquent, tel que montr la figure 8.5(a) la flexion de cette partie nest
pas cause par une charge transversale au PMG de lme comme dans la thorie classique
des plaques en stratifi de composites. La matrice [D] ne joue donc, aucun rle dans la
flexion de lme.
La figure 8.5(b) illustre la contrainte axiale x de lme au niveau partir de laxe
neutre de la poutre, qui est cause par le moment appliqu sur la poutre. Elle est
11
2014
(8.19)
do :
=
N x (me) (
1
'
11(me)
) x
(8.20)
h /2
h /2
N=
x d
h /2
h /2
1 2
) x d
'
A11
ou :
M me
h3
= (
)x = Ja x
'
12A11
(8.21)
x
Nx
0
1
[ A ] 0 =
et y
0
0
xy
0
[ D]
M x
0
0
(8.22)
do :
N=
(
x
1 0
)=
x ( 1' ) 0x
'
A11
A11
et M=
(
x
1 0
)x
'
D11
(8.23)
12
2014
12
1
(8.24)
+ ' ) 0x
'
A11 D11
Le moment de flexion total pour la poutre est la somme des moments appliqus lme et
aux deux ailes :
=
M aile b aile (N x 1 +=
M x ) baile (
M poutre =
M me + 2M aile =
J a 0x
do la rigidit apparente de la poutre est :
2
h3
1
Ja = '
)
+ 2b aile ( ' 1 + '
12A11(me)
A11(aile) D11(aile)
(8.25)
(8.26)
0y =0xy =0
La rigidit apparente de la poutre scrit :
12
h3
Ja = '
+ 2b aile ( '
+ D11(aile) )
12A11(me)
A11(aile)
(8.27)
Exemple 8.5
La poutre illustre la figure 8.6 est faite de stratifi de composite carbone/poxy
AS3501 (v f =.65). Sachant que la charge rpartie applique est 1.75N/mm et que
lpaisseur des plis est 0.1321mm, calculez :
1. La rigidit apparente de la poutre;
2. La flche maximale;
3. Les contraintes axiales au ke pli de lme et de laile.
13
2014
Aile [90/45/-45]s
19mm
me [90/45/-45]s
12.7mm
Figure 8.6
Solution
Pli#
1
Matriau
AS3501
2
3
Direction (o)
90
45
-45
Hauteur (mm)
Z0=-0.3963
Z1=-0.2642
Z2=-0.1321
Z3=0.0000
Programme sixplis_sym_v3
b aile =12.7mm
h=19mm
1 =9.5mm+3x0.1321mm=9.8963mm
L=508mm
qL2
M=
8
M qL2 1.75 5082
M=
=
=
= 4445N
x
b
8b
8 12.7
me:
'
=
A11
0.4685e 4mm / N
Aile:
'
=
A11
0.4685e 4mm / N
'
=
D11
0.0015N 1 mm 1
14
2014
2
h3
1
+ 2b aile ( ' 1 + '
Ja = '
)
12A11(me)
A11(aile) D11(aile)
=
Ja
9.896242
193
1
+
=
+
2
12.7(
) 65285691.9N mm
4
4
12 0.4687 10
0.4687 10
0.0015
La flche maximale :
=
w
0max
5qL4
5qL4
5 1.75 5084
= =
= 23.24mm
384E fx I 384J a 384 65285691.9
(x )k (Q11 )k 0x et (y )k (Q12 )k 0x
Aile troite :
Les dformations au ke pli de laile suprieure sont :
15
centre poutre
centre
2014
M
Ja
poutre k
'
D12
( x )au
'
D11
centre poutre
(x ) k (Q11 ) k ( x ) k (Q12 ) k ( y ) k ;
( y ) k (Q12 ) k ( x ) k (Q 22 ) k ( y ) k
Aile large :
( y ) 0 ( y ) k 0
do :
_
(x ) k (Q11 ) k ( x ) k ;
_
( y ) k (Q12 ) k ( x ) k
8.3
Dans le cas des poutres non symtriques il existe le couplage entre les contraintes et
moments rsultants. La relation charges-dformations est :
0(xy)
=
0
(xy)
A | B
B | D
N
=
M
A'
'
B
B' N
| D' M
|
(8.28)
Remarque : Les matrices [A], [B] et [D] ne sont pas les matrices inverses des matrices
[A], [B] et [D]. Seul dans le cas des stratifis symtriques que [A]=[A]-1 etc..
Lquation (8.28) avec N x =N y =N xy =M y =M xy = 0 donne les relations suivantes :
'
'
=
0x B11
M=
0x D11
Mx
x
(8.29)
Dans le cas des poutres non symtriques, laxe neutre ne passe pas au milieu de la
hauteur de la poutre mais se trouve un niveau o :
16
2014
'
'
0 = x =0x + z n 0x =B11
M x + z n D11
Mx
do :
zn =
'
B11
'
D11
(8.30)
Exemple 8.6
Soit une poutre en appuis simples faite de stratifi [0/90/-30/30] 2 de composite
unidirectionnel carbone/poxy AS3501 (v f =0.65). Les dimensions de la poutre sont L=
100mm et b= 5mm et lpaisseur des plis est de 0.125mm. Sachant que la charge
transverse q = 0.2N/mm, calculez;
1. La flche maximale;
2. Les contraintes au niveau infrieur du 2e pli.
Solution
Matriau AS3501 (v f =0.65) : E 1 =138GPa, E 2 =9GPa, G 12 =6.9GPa, Nu 12 =0.3.
q=0.2 N/mm, L=100mm, b=5mm
Pli#
1
Orientation (o)
1 =0
2
3
4
5
6
7
8
2 =90
3 =-30
4 =30
1 =0
2 =90
3 =-30
4 =30
Hauteur (mm)
Z 0 =-0.500
Z 1 =-0.375
Z 2 =-0.250
Z 3 =-0.125
Z 4 =0.000
Z 5 =0.125
Z 6 =0.250
Z 7 =0.375
Z 8 =0.500
Programme huit_plis_v3
Aile troite:
'
D11
= 0.000153583383205132(N-mm) 1
=
Ja
b
5
=
=
'
D11 0.000153583383205132
32555.6052722305N mm 2
Flche maximale :
17
2014
5qL4
5qL4
5 0.2 1004
= =
= 7.999mm
384E fx I 384J a 384 32555.6052722305
w
=
0max
0( x ,y )
=
0
( x ,y )
A'
B'
| B' N (x,y)
| D' M (x,y)
{ }=
(xy) k
{ }
{ }
_
_
0
0
Q
=
{
}
(xy)
Q ( ( xy ) + z k ( xy ) )
k
k
k
8.4
Contraintes de cisaillement xz
Dans le cas de flexion cylindrique ou flexion pure considre dans les sections
prcdentes le moment de flexion est considr constant et leffort tranchant est nul. En
pratique, dans une poutre soumise une charge quelconque, le moment flchissant passe
18
2014
par un extremum lorsque leffort tranchant est nul. Il est admis que lorsque la longueur de
la poutre est trs grande par rapport aux deux autres dimensions, la formule pour calculer
la contrainte axiale x est acceptable malgr une certaine faible erreur.
Lorsque leffort tranchant existe, la contrainte x varie par consquent do lexistence de
la contrainte de cisaillement xz . Cest le cas de flexion normale o leffort V0 qui est
abord ici.
Ligne de
centre
x
z
(a)
dx
dz
h/2
X + dX
X
xz
Axe
neutre
(b)
( x + x dx)bdz
x
h /2
(8.31)
do :
h /2
xz =
x
dz
x
19
2014
x= Q11 x + Q12 y
En remplaant lintgrale par rapport la hauteur du stratifi par la somme des intgrales
par rapport lpaisseur de chacun des plis et lexpression de x , lquation (8.31) peut
scrire comme :
k
k
_
_
=
xz
=
+
( (
)
dz)
(
[(Q
)
(Q
x k
11 k x
12 ) k y ]dz
x k zk 1
x k zk1
(8.32)
=
xz
k
_
D' M
d
M _
( [(Q11 ) k z + (Q12 ) k z 21
]dz
'
dx k zk 1
J
J
D11
_
_
D' k
dM 1
=
( ( [(Q11 ) k + (Q12 ) k 21
] zdz
'
dx J k
D11
k 1
(8.33)
dM
dx
=
xz
_
_
z 2 z k2 1
V
( [(Q11 ) k + (Q12 ) k ]( k
)
J k
2
(8.34)
o :
=
D'21
pour aile troite et =0 pour aile large.
'
D11
Exemple 8.7
Soit une poutre en appuis simples, faite de stratifi [0 4 /45 4 /-45 4] S du composite
unidirectionnel carbone/poxy AS4/3501-6( E 1 =131GPa, E 2 =11.2GPa, G 12 =6.55GPa,
nu12=0.28).
20
2014
Les dimensions de la poutre sont L= 254mm et b= 50.8mm et lpaisseur des plis est de
0.125mm. Sachant que la charge transverse q = 0.3N/mm, calculez la contrainte xz au
PMG de la poutre.
Solution
La poutre est considre comme poutre aile large cause du rapport : largeur/paisseur
b/h= 50.8/3 =16.9.
Programme sixplis_sym_v3 donne :
Groupe
de pli#
Orientation (o)
Hauteur
(mm)
1 =0
2
3
2 =45
3 =-45
z 0 =-1.5
z 1 =-1.0
z 2 =-0.5
z 3 =0
D 11
Q11
(N/mm)
13184.01
Q12
(N/mm)
3157.16
238095.33
43918,48
43918,48
30818,48
30818,48
238095.33
238095.33
J a ==
bD11 50.8 238095.33 =
12095242.60 N mm 2
V
=
qL 0.3 254
=
= 38.10 N
2
2
= -0.3836MPa
Commande Matlab :
Txz1=V/Ja*((qb1(1,1)+qb1(1,2))*0.5*(z1^2-z0^2))
Txz2=V/Ja*((qb1(1,1)+qb1(1,2))*0.5*(z1^2-z0^2)+(qb2(1,1)+qb2(1,2))*0.5*(z2^2z1^2))
Txz3=V/Ja*((qb1(1,1)+qb1(1,2))*0.5*(z1^2-z0^2)+(qb2(1,1)+qb2(1,2))*0.5*(z2^2z1^2)+(qb3(1,1)+qb3(1,2))*0.5*(z3^2-z2^2))
21
2014
8.4.2 Section en I
me :
Pli#k
tme
x
tk
(x)k
hme
Nx
=
(xz ) k
( x ) k
d
x
(8.35)
Considrant :
(x ) k (Q11 ) k x (Q12 ) k y ;
h me t me y 0 y 0;
x x
M
Jk
h /2
(Q11 ) k dM
d
M
=
(xz ) k (Q11 ) k =
( )d
dx J k
J k dx
or :
h /2
dM
dx
22
2014
do :
_
(Q ) V 1 2 h /2 V Q11 1 h 2
=
(xz ) k ( 11 k )( =
)(
( )( 2 )t k
Jk
tk 2
Jt k 2 4
La contrainte de cisaillement xz de lme est :
V
h2
1 N _
(xz )me
= ( xz =
)k
( 2 )( ) (Q11 ) k t k
Ja me t me 4
2
k
k
A
11
ou :
(xz )me ==
VA11
h2
( 2 )
2Ja me t me 4
(8.36)
Aile :
Dans le cas dune poutre symtrique dont aile est galement un stratifi symtrique :
_
_
z 2k z 2k 1
V
xz )aile
(=
( [(Q11 ) k + (Q12 ) k ](
)
J k
2
(8.37)
D'21
pour aile troite et =0 pour aile large.
D'11
23
2014
24