MATERIAUX 2006 13-17 Novembre 2006 Dijon, France

Homognisation numrique des proprits lastiques dun matriau sandwich

Moussa KARAMA , Lama ALBACHI
a,b a b

Ecole Nationale dIngnieurs de Tarbes/ LGP-Equipe CMAO, 47 avenue dAzereix BP 1629 65016 Tarbes cedex, moussa@[Link]

RESUME : Les matriaux sandwiches structure nid d'abeilles sont un bon compromis entre rsistance et lgret et leur coefficient de dilatation est proche de zro. Ces proprits sont trs recherches par exemple pour les rflecteurs et les antennes des satellites de communications et les structures aronautiques. Etant donn le cot de l'utilisation de ces matriaux pour des prototypes, il est intressant de connatre leurs proprits mcaniques pour prdire et calculer leur comportement dans des environnements spcifiques et varis. Ce papier est consacr au problme de lestimation correcte des caractristiques mcaniques des nids dabeilles (modules dYoung, coefficients de Poisson et modules de cisaillement). aprs des rappels thoriques sur la prvision de ces modules en fonction de la gomtrie des cellules, une estimation plus fine est effectue grce des calculs par lments finis.

MOTS-CLES : matriaux sandwich, nida, homognisation, modlisation analytique, modlisation numrique

Introduction

L'utilisation des composites dans l'arospatiale, l'aronautique et l'automobile est en croissance permanente sous toute forme de configurations. Ces domaines d'activit sont sans cesse la recherche de produits nouveaux plus rsistants, plus lgers ou insensibles aux changements de temprature. Les matriaux sandwiches sont composs de plusieurs parties interagissant entre elles. On trouve donc les peaux ou semelles ayant des caractristiques mcaniques propres et le cur ou me pouvant tre entre autres du nid d'abeilles (nida) ayant des proprits mcaniques spcifiques. L'assemblage de ces deux parties se ralise par collage ou soudage mettant en jeu un autre corps de comportement diffrent. L'intrt de l'tude est de pouvoir homogniser le nida indpendamment des peaux, c'est--dire remplacer le nida de gomtrie complique par un matriau homogne de proprits mcaniques similaires dans les modles analytiques et numriques. Cela permettra par la suite de prvoir le comportement mcanique des diffrentes combinaisons possibles nida/peaux en connaissant, d'une part les caractristiques du nida et, d'autre part celles des semelles. La connaissance du comportement mcanique des mes, en particulier en cisaillement, est videmment indispensable lors de la phase de conception des structures sandwiches. Lestimation de ce comportement consiste sappuyer sur la structure rptitive des nids dabeilles pour raliser des calculs dhomognisation. On peut alors prdire de nombreux coefficients lastiques [1]. Toutefois, dans le cas gnral, les modules de cisaillement peuvent seulement tre borns avec ce type dapproche [2].

MATERIAUX 2006 13-17 Novembre 2006 Dijon, France Le nid d'abeilles est assimilable un matriau orthotrope de part ses proprits directionnelles et ses plans de symtries. Son homognisation passe tout d'abord par la dtermination des constantes lastiques de l'ingnieur qui sont les modules d'lasticit, les coefficients de Poisson et les modules de cisaillement.

Caractristiques lastiques dun nid dabeille

Il est noter que suivant les littratures, les dnominations et indices peuvent varier. La gomtrie ainsi que les notations utilises par la suite sont dtailles ci-aprs. les modules d'lasticit : Ei, o Xi est la direction de traction. les coefficients de Poisson: ij , o Xi est la direction de traction et Xj est la direction de la contraction. les modules de cisaillement : Gij module de cisaillement dans le plan (Xi, Xj).

X2
t t'

X1

Figure 1. Gomtrie dun nid dabeilles Un nid dabeilles peut tre caractris par un angle de cellule et quatre paramtres dfinis partir des caractristiques suivantes: -

Angle de cellule Diamtre de cellule

b a t t h Longueur dune paroi centrale Longueur dune paroi incline Epaisseur dune paroi centrale Epaisseur dune paroi incline Hauteur du nid dabeilles

Le nid d'abeilles peut tre assimil un matriau orthotrope. Pour le caractriser, il faut donc dterminer les quinze constantes de l'ingnieur: les modules d'lasticit E1, E2, E3, les coefficients de Poisson 12, 21, 23, 32, 13, 31, les modules de cisaillement G12, G21, G23, G32, G13, G31. La dtermination de ces caractristiques mcaniques peut tre effectue grce deux mthodes diffrentes : la modlisation analytique, la modlisation numrique.

MATERIAUX 2006 13-17 Novembre 2006 Dijon, France La mthode dhomognisation Il est ncessaire pour des raisons de temps de calcul et de puissance de calcul disponible de rduire les modles et donc d'utiliser les symtries du modle. Dans le cas du nida, le nombre de modles pouvant tre crs est abondant. Les modles utiliss par la suite sont guids par les tudes et les simulations effectues sur le nid d'abeilles par de nombreuses quipes de recherche [3], [2], afin de pouvoir au pralable comparer les rsultats des diffrentes simulations, amliorer et adapter si possible les modles au code de calcul Samcef [4]. La gomtrie est base de cellules hexagonales rgulires. Si l'on se rfre aux notations de la figure 1, les caractristiques sont les suivantes : = 6,35 mm, a = b = 3,66 mm, h = 9,525 mm, = 30 Les longueurs . lx, ly, lz du volume lmentaire reprsentatif sont (Figure 2) : lx = 54,986 mm, ly = 25,4 mm, lz = 9,525 mm. L'paisseur t est paramtre. Pour ce modle, toutes les simulations sont effectues avec quatre valeurs de t et la double paisseur la liaison entre deux cellules adjacentes est prise en compte. Le matriau considr est de l'aluminium 5056 avec les caractristiques suivantes : E0 = 68,940.103 MPa, G = 6,54.103 MPa, = 0,312. En comparant les rsultats obtenus partir de ces deux mthodes, on peut ainsi effectuer une approximation des constantes de l'ingnieur du nida. Les simulations ont t effectues sur le code de calcul Samcef. Le modle de nida prend en compte 40 cellules qui ont t mailles avec des lments plaques 4 nuds 6 degrs de libert par nud. Pour ce modle, le maillage choisi permet de ne pas trop alourdir le calcul tout en gardant une bonne prcision au niveau des rsultats. Un raffinement excessif ne donne pas de meilleurs rsultats donc le maillage retenu compte 1500 lments et 1725 nuds. Chaque face de l'hexagone (clinquant) est discrtise par 12 lments, 4 suivant l'axe Z et 3 suivant l'axe X.

lz

ly YZ 0 X

lx

Figure 2 : Volume lmentaire reprsentatif V.E.R.

Pour calculer les diffrents modules, un dplacement est impos sur une face du V.E.R. dans une direction donne, la face oppose tant fixe. Les symtries sont prises en compte au niveau des conditions limites. Neuf simulations sont ncessaires pour dterminer les neuf constantes lastiques du nida : Pour calculer les diffrents modules, un dplacement est impos sur une face du V.E.R. dans une direction donne, la face oppose tant fixe. Les symtries sont prises en compte au niveau des conditions limites. Neuf simulations sont ncessaires pour dterminer les neuf constantes lastiques du nida :

La simulation dune sollicitation de traction suivant une direction i (i = x, y, z) permet de dterminer les 3 modules de rigidit E1, E2 et E3 et les 6 coefficients de Poisson (Figure 3 a-c). V0 : volume avant dformation, Vf : volume aprs dformation.

MATERIAUX 2006 13-17 Novembre 2006 Dijon, France

V0

Vf Ux x

Vf

y V0 Uy x

V0

Vf Uz x

(a) E1, 12, 13

(b) E2, 21, 23

(c) E3, 31, 32

Figure 3. VER en dplacement impos suivant X (a) ; Y (b) ; Z (c)

La simulation des sollicitations de cisaillement suivant les plans xy, yz, xz (2 simulations par plan) permet de dterminer les 6 coefficients de cisaillement (Figure 4 d-i).

Ux

V0

y Vf

V0

Vf Uy

(d) G12 = 12 / 12

(e) G21 = 21 / 21

Ux

V0 Vf

Vf V0 Uz

x
(f) G13 = 13 / 13 (g) G31 = 31 / 31

Uy

Vf

Vf Uz

V0
(h) G32 = 32 / 32

V0
(i) G23 = 23 / 23

Figure 4. VER sous sollicitation de cisaillement (d-i) Cette mthode base sur l'analyse des efforts aux liaisons, partir dun dplacement impos, permet par lutilisation des lois de comportement de dduire les constantes lastiques.

MATERIAUX 2006 13-17 Novembre 2006 Dijon, France

Rsultats Modle de Gibson : Le tableau 1 prsente les rsultats obtenus par simulations numriques pour une srie de conditions aux limites (dplacement impos suivant la direction X ) et pour 4 paisseurs (t) du feuillard et les compare aux valeurs analytiques obtenues par les formules de Gibson suivantes :

.
cos t E1 = E o a ( b a + sin ) sin 2
3

(1)

21 =

( b a + sin ) sin cos 2

cos 2 ( b a + sin ) sin

(2)

12 =

(3)

Samcef

Gibson Erreur %
Erreur %

Lx (mm) Ly (mm) Lz (mm) x y z Fx (N) x (MPa) E1 (MPa) 12 13 E1 (MPa) 12 13 E1 13

t (mm) 0,018 0,025 0,060 0,127 1 1 1 1 -4,62E-01 -4,62E-01 -4,62E-01 -4,61E-01 -1,43E-06 -2,83E-06 -1,63E-05 -7,29E-05 1,82E-02 1,82E-02 1,82E-02 1,82E-02 -1,82E-02 -1,82E-02 -1,82E-02 -1,81E-02 -1,50E-07 -2,97E-07 -1,71E-06 -7,66E-06 8,02E-02 2,23E-01 3,08E+00 2,92E+01 3,32E-04 9,22E-04 1,27E-02 1,21E-01 0,018 0,051 0,701 6,635 1 0,9999 0,9994 0,997 8,26E-06 1,63E-05 9,41E-05 4,21E-04 0,018 0,050 0,698 6,618 1 1 1 1 0 0 0 0 0,41 0,41 0,37 0,26
0 0,01 0,06 0,27

Tableau 1. Rsultats de la traction suivant X. On peut remarquer que l'cart des rsultats entre la simulation numrique et le modle de analytique est faible puisque l'erreur maximale pour E1 est de 0,41%. Lerreur sur les coefficients de Poisson est trs faible. Le tableau 2 prsente pour 4 paisseurs de feuillard les rsultats obtenus pour le module de cisaillement G23 et leur comparaison avec ceux obtenus par le modle de Gibson.

G23 = G

cos t b a + sin a

(4)

MATERIAUX 2006 13-17 Novembre 2006 Dijon, France

0,018 Ly (mm) 0,1 yz 0,0105 Fy (N) 1078,9 yz (MPa) 0,77 Samcef G23 (MPa) 73,57 Gibson G23 (MPa) 74,31 Erreur % G32 1

G23

t (mm) 0,025 0,060 0,1 0,1 0,0105 0,0105 1517 3641,2 1,09 2,61 103,45 248,30 104,49 250,77 1 0,99

0,127 0,1 0,0105 7707,6 5,52 525,59 530,80 0,98

Tableau 2. Rsultats en cisaillement pour un dplacement impos suivant la direction z.

Modle de Grdiac : Grdiac s'est intress particulirement aux modules de cisaillement hors-plan. Il s'est inspir du modle de Gibson, mais comparativement lui, il a pris en compte galement la double paisseur des clinquants centraux. Grdiac a dfini le nida par quatre grandeurs adimensionnelles :

t b h t , R2 = , R3 = , R4 = a a a t Deux cas sont diffrencis : R1 =

- 1 cas :
er

(5)

R 4 = 2 , G est le module de cisaillement du matriau utilis :

1 + R 2 sin 1 + R 2 sin 2 R G R1 G G xz R 2 (1 + R 2 ) cos (1 + R2 sin ) R2 cos 1

et

(6)

G yz =

cos R1 G 1 + R 2 sin

(7)

- 2 cas :

R4 = 1

1 + R 2 sin 1 + 2 R 2 sin 2 R1 G G xz R1 G R 2 (2 + R 2 ) cos 2 (1 + R 2 sin ) R 2 cos

et :

(8)

G yz =

cos R1 G 1 + R 2 sin

(9)

et lexpression analytique pour le module de cisaillement Gxz o llancement R3 est pris en compte. Cette expression est du type :
min G xz = G xz i +

0.787 max i min G xz G xz i R3

(10)

On peut noter que pour les cellules hexagonales rgulires (=30 avec R4 = 1 , les modules sont gaux ) quelle que soit la direction. Le problme de ce modle est qu'il ne tient pas compte du paramtre R3. Quand R3 est petit, il ne donne pas une bonne valuation des modules.

MATERIAUX 2006 13-17 Novembre 2006 Dijon, France Le cas le plus courant correspond des nids dabeilles fabriqus partir de feuilles de mtal dcoupes, embouties et colles les unes sur les autres. Dans ce cas, t = 2 t . On a alors les relations qui changent sont les suivantes : I moment d'inertie suivant X3 du clinquant centrale BD : I =

h (2 t )3 12

b a + sin t G12 = G yx = E o a (b a )2 (1 + b 4 a ) cos

(11)

E3 = Eo

Section de la cellule 2 Surface hexagone

, Section de la cellule = 4 a t + 4 b t

(12)

b a +1 t E3 = E z = Eo a (b a + sin ) cos
et la contrainte de cisaillement 23.

(13)

b a + sin b a + sin 2 t t G G xz = G23 G cos (b a + sin ) a a (1 + b a ) cos

(14)

La synthse de lensemble des simulations effectues (avec les diffrents modles) pour la dtermination de lensemble des caractristiques du nida est ralise dans le tableau 3.

Nida6 E1 (MPa) E2 (MPa) E3 (MPa) 12 13 21 23 31 32 G12 (MPa) G21 (MPa) G23 (MPa) G23 (MPa) G13 (MPa) G13-mini (MPa) G13-maxi (MPa) G31 (MPa)

Samcef Gibson Grdiac 0,886 0,886 1954,97 0,999 0 0,999 0 0,33 0,33 0,475 0,17 238 57 389,2 0,828 0,828 1840 1 0 1 0 0,33 0,33 0,497 250 250 383 375 416,7

375 416,7 0,92

Tableau 3 : Caractristiques du nida obtenues par diffrents modles

Conclusion Les rsultats sont assez proches car l'cart entre Samcef et Gibson est d'environ 1%. Ce module est trs important. Il a le mme rle que le module G13 dans le comportement des structures sandwiches.

MATERIAUX 2006 13-17 Novembre 2006 Dijon, France

Le raffinement des conditions aux limites par couplage des nuds que nous avons dvelopp, permet de fortement recaler les simulations numriques par rapport au modle analytique de Gibson. On constate que la dtermination des modules dYoung E1 et E2 par les deux mthodes analytique et numrique est correcte car lcart entre le modle de Gibson et la simulation numrique atteint 0,41% pour E1 et E2. Par contre lerreur sur les coefficients de Poisson est toujours trs faible. Nous avons rencontr cependant un problme dans la dtermination de E3 qui savre assez lev par rapport au modle de Gibson (10%). Nous sommes cependant moins loigns que les donnes des fournisseurs tel que Hexcel [Hexcel, 1994], qui pour un clinquant dpaisseur t = 0,018 mm donnent un module E3 de 207 MPa, ce qui correspond un cart de 39%. Pour ce qui concerne les modules de cisaillement les carts entre les simulations numriques et les simulations analytiques sont relativement faibles pour les quatre paisseurs de clinquant que nous avons considres. Les rsultats que nous obtenons, pour un feuillard dpaisseur t = 0,018 mm, sont trs proches des donnes de Hexcel. En effet l'erreur sur le module de cisaillement G13 est de 5%, celle sur le module de G23 est de 12%.

References [1] L. J. Gibson & M. F Ashby, "Cellular solids, structure and properties", Pergamon Press., 1988. [2] M. Grdiac, "A finite element study of the transverse shear in honeycomb cores", Int. J. Solids structures, vol. 30, p. 1777-1788, 1993. [3] C. C. Chamis, R. A. Aiello, L. N. Murthy, Fiber composite sandwich thermostructural behaviour: Computational simulation, Journal of composites technology & research, vol. 10, p. 93-99,1988. [4] L. Al Bachi, Etude numrique et exprimentale du comportement d'une structure cellulaire priodique, application aux structures sandwiches, thse INPT deToulouse, juillet 2002. [5] Hexcel , Design of sandwich structures, Structural Materials Handbook, Section IV, 1994.