Identication dune loi de comportement biaxiale en grandes dformations pour des tricots lastomriques

Identication of elastomeric knitted fabrics

Issam GAIED* , Sylvain DRAPIER**
* Dpartement Mcanique et Matriaux

Centre Science des Matriaux et des Structure ** Dpartement Biomatriaux et Mcanique Centre Ingnierie et Sant 158, cours Fauriel 42023 Saint-tienne cedex 2 France drapier@emse.fr

RSUM. Les matriaux tricots ont intrinsquement un comportement orthotrope. Dans le cas des tricots lastomriques tudis ici, dont les dformations normales en fonctionnement sont comprises entre 80 et 120 %, le comportement biaxial en grandes dformations doit tre identi. Pour cela, une mthodologie spcique a t mise en place o les mesures exprimentales sont obtenues par corrlation dimages et les simulations reposent sur des modles lments nis o une formulation corrotationnelle a t implmente. Lidentication des surfaces de rponse est pilote par un algorithme de Levenberg-Marquardt qui conduit minimiser lcart, calcul au sens des moindre carrs linaires, entre les vecteurs des grandeurs mesures et simules. ABSTRACT. Knitted materials have intrinsically an orthotropic behaviour. In our elastomeric knitted fabrics, the use of which requires normal strains ranging from 80 to 120 %, the biaxial response must be characterized under large deformations. To achieve this, a specic methodology has been proposed in which experimental measures are obtained through a digital image correlation technique and simulated results rely on nite element models using a corrotationnal formulation implemented specically. The identication of the response surfaces is controlled by a Levenberg-Marquardt algorithm able to minimize the distance, calculated in the sense of non-linear least squares, between both vectors of measured and simulated results. MOTS-CLS : KEYWORDS:

Identication, grandes dformations, biaxial, corrlation dimages, modles EF Identication, nite strains, biaxial tension, digital image correlation, FE models

Nom de la Revue. Volume X - n X/1999, pages 1 X

Nom de la Revue. Volume X - n X/1999

1. Introduction Les techniques didentication par mthode inverse couvrent aujourdhui de nombreux domaines de la mcanique des matriaux et des structures. Ces mthodes permettent aussi bien didentier des proprits mcaniques que de caractriser la prsence de discontinuits par exemple [BON 98]. Pourtant, lapplication de ces mthodes au domaine des grandes dformations est encore trs limite. Dans le cas des tricots lastomriques que nous tudions plus particulirement, le domaine des grandes dformations constitue le domaine dutilisation permettant de gnrer des pressions de contention sur les membres humains. Ces matriaux structurs, possdant par nature deux directions privilgies, supportent des dformations normales comprises entre 80 et 120 %. Les efforts mis en jeu dans ces tricots sont 1 . Pourtant, la caractrisation ne du comportepar contre de lordre de 0,5 ment biaxial de ces tricots est essentielle, elle sinscrit dans un projet de prvision des distributions des pressions de contention dlivres sur un membre humain.

Dans cet article, nous prsentons une mthode spcique didentication de ce comportement biaxial en grandes dformations, base sur lutilisation dune mthode exprimentale sans contact et sur des simulations lments nis avec le code Zbulon c o une formulation corrotationnelle a t implmente spciquement. Lidentication de la loi de comportement postule utilise lalgorithme de Levenberg-Marquardt qui permet de minimiser la fonction objectif dnie comme la mesure, au sens des moindre carrs non-linaires, de lcart entre les vecteurs mesurs et simuls.

2. Stratgie didentication Par rapport lidentication dun comportement en petites dformations, le domaine des grandes dformations impose de traiter les mesures, autant que les simulations, de faon incrmentale. Comme indiqu sur la gure 1, lalgorithme didentication se compose dune partie exprimentale indpendante de la gnration successive des rsultats simuls.

2.1. Mesures exprimentales Compte-tenu de la discontinuit des structures tricotes, et des grandes dformations, les mesures exprimentales sont ralises laide dune technique de corrlation dimages numriques, mthode de champ sans contact. Dans cette mthode, partant dune paire dimages, la position dune zone (fentre) de la premire image est recherche, dans la seconde image, par maximisation dune fonction de corrlation re-

. force / unit de longueur : mesure de force et rigidit utilise dans les tricots ou tissus dont
lpaisseur est difcilement contrlable

Identication tricots lastomriques

Figure 1. Principe de lidentication en grandes dformations.

prsentant la similitude des fentres au sens des niveaux de gris quelles contiennent. La granularit naturelle des tricots fait ici ofce de traceurs [GAI 05]. Le montage dvelopp permet de solliciter le tricot en traction sans le dtriorer (structure tricote fragile), mais surtout la camra CCD utilise pour lacquisition des images est xe sur un pantographe, ce qui permet de viser le mme point de rfrence pendant tout lessai. Ceci est ncessaire pour pouvoir, lors du traitement des images, extraire des dplacements absolus de points gomtriques partir des dplacements relatifs entre images. Enn, ce montage est quip dun systme de tension transverse permettant de maintenir la dimension transverse xe. Le coefcient de liage 2 des tri-

. assimilable un coefcient de Poisson pour les milieux continus, au sens des dformations
libres induites dans le sens transverse la sollicitation

Nom de la Revue. Volume X - n X/1999

cots tant trs lev, les dformations induites sont de lordre des dformations dans la direction de sollicitation [GAI 04].

2.2. Modlisation choisie Le modle lments nis de notre essai, qui nous permet dobtenir des rsultats simuls comparer aux rsultats exprimentaux, est ralis avec le code de calculs Zbulon c . Le maillage utilis est bidimensionnel, compos dlments quadratiques en contraintes planes.

La spcicit rside ici dans la formulation corrotationnelle quil est ncessaire de prendre en charge en grandes dformations. En utilisant le mta-language Zeb-front, nous avons implment cette formulation en retenant une mesure logarithmique pour les dformations. Cette formulation nous permet galement de vrier par simulation que dans notre cas de traction biaxiale, le tenseur de rotation rsultant de la dcomposition polaire de la dformation se rduit lidentit mme pour un matriau isotrope. Ceci correspond bien au comportement attendu pour nos tricots : deux directions privilgies, et absence de cisaillement.

2.3. Identication Finalement, le schma didentication peut-tre explicit. Les grandeurs exprimentales (forces aux bords et dplacements) sont mesures pour former le vecteur . La procdure didentication numrique consiste alors fournir un jeu de paramtres initial pour la loi de comportement choisie, puis gnrer successivement des jeux de rsultats simuls pour former la mesure de la distance entre les grandeurs mesures et simules . Il est noter que les rsultats des simulations non-linaires doivent tre au pralable extrapols aux points de mesures gomtriques pour tre compars aux mesures exprimentales. Si la fonction objectif ne satisfait pas le critre de convergence, un nouveau jeu de paramtres est alors gnr, travers lalgorithme de Levenberg-Marquardt de Z-optim [LER 99]. Le jeu de paramtres qui permet de satisfaire le critre de convergence correspond au meilleur jeu pour la loi de comportement postule.

 





!#"%$'& 
)( *
!+&-,  
  

3. Tests et rsultats La procdure didentication dcrite ci-dessus (gure 1) a t teste, dabord en 1D sur une loi hyperlastique dont la forme est similaire la rponse macroscopique de nos tricots. Nous avons mis en vidence que la comparaison des champs de dplacement seuls ou des forces seules ne permettait pas de retrouver un loi connue. Cest en combinant les dplacements lintrieur du tricot et les forces au bord que les calculs convergent le plus vite et vers la solution connue. Ensuite, ces mmes essais ont t effectus en 2D sur des surfaces de rponse gnres.

Identication tricots lastomriques

Le cas des tricots lastomriques a ensuite t trait. Une rponse pragamatique a t recherche, elle a t pose par observation de la rponse uniaxiale en sens trame et chane des tricot. La rponse polynmiale est dordre 6 en sens trame, 4 en sens chane, et 2 pour le couplage trame-chane, soit un total de 18 coefcients identier. Sur la gure 2 sont reports les rsultats pour lidentication du comportement dun tricot Diaphane 2 commercialis par la socit Ganzoni France S.A. Dans ce cas prcis, 28 images ont t prises pour 4 longations diffrentes dans le sens transverse, soit 112 pas de calculs en termes de simulations EF. Le temps didentication est de 3,5 h sur un PC PIII. Les rsultats montrent un bon accord entre les surfaces de rponse identies et les mesures exprimentales. Lcart qui apparat aux faibles et grandes dformations, notamment en sens chane (gure 2-b), montre la difcult didentier simultanment la rponse complte du matriau. Dun point de vue pragmatique, ces tricots sont utiliss pour des dformations sens chane infrieures 50 %, ce qui est conrm par des essais o une dformation sens trame de 100 % en prsence dune dformation sens chane suprieure 70 % (gure 2-b) conduit une dgradation de la structure du tricot. On peut conclure que dans les zones dutilisation courantes, les surfaces identies reproduisent bien le comportement du tricot, compte-tenu des dformations leves, mais aussi et surtout du faible niveau deffort rencontr.

4. Conclusion Nous avons montr dans cet article, les possibilits didentication dune loi de comportement en grandes dformations, ds lors que les contraintes de mesures et de reprsentation sont traites avec des outils adapts. Les calculs didentication ne sont pas coteux si le problme inverse est bien conditionn, en combinant les grandeurs physiques qui constituent le vecteur des grandeurs comparer. Il faut maintenant afner la reprsentation des surfaces laide de rponses mieux adaptes, la difcult principale de mesure tant leve.

5. Bibliographie
[BON 98] B ONNET M., R EYNIER M., On the estimation of the geometrical support of modelling defects using the distributed error in constitutive equation , Proceedings of Inverse problems, control and shape optimisation, April 8-10, Carthage, Tunisia, 1998, p. 6570. [GAI 04] G AIED I., Contribution ltude de la contention : approche bidimensionnelle, identication du comportement biaxial en grandes dformations des tricots lastomriques , PhD thesis, cole Nationale Suprieure des Mines de Saint-tienne, 2004. [GAI 05] G AIED I., D RAPIER S., Identication de la loi de comportement dun tricot lastomrique , Revue des Matriaux Composites et Avancs, vol. paratre, 2005. [LER 99] L E R ICHE R., G UYON F., Least squares parameter estimation and the LevenbergMarquardt algorithm : deterministic analysis, sensitivities and numerical experiments , Rapport technique, 1999, Laboratoire de Mcanique de Rouen, INSA de Rouen.

Nom de la Revue. Volume X - n X/1999

(a)

(b) Figure 2. Rsultats de lidentication du comportement : forces en fonction des dformations sens trame et chane (a) force sens trame et (b) force sens chane.

Identication tricots lastomriques

Annexe pour le service de fabrication

Article pour les actes : Giens 2005 Auteurs : Issam GAIED* , Sylvain DRAPIER** Titre de larticle : Identication dune loi de comportement biaxiale en grandes dformations pour des tricots lastomriques Titre abrg : Identication tricots lastomriques Traduction du titre : Identication of elastomeric knitted fabrics Date de cette version : 10 fvrier 2005 Coordonnes des auteurs : tlphone : +33.64.77.42.00.79 tlcopie : +33.64.77.42.02.49 Email : drapier@emse.fr Logiciel utilis pour la prparation de cet article : A L TEX, avec le chier de style version 1.10 du 17/09/2001.

./1032547658:9 ;<87/>=38<?A@B476 ? ,

Formulaire de copyright : Joindre le formulaire de copyright sign, rcupr sur le web ladresse

;C0505DFEHG1GI1I1IF@J; 8/>=K8 ?195?1452>8>LK4>8M@N4O=