Analyse mésoscopique par éléments finis de la déformation de renforts fibreux 2D et 3D à partir de microtomographies X

Meso-scale fe analysis of 2D and 3D fibrous reinforcements deformation based on X-ray computed tomography

• Auteur / Autrice : Naïm Naouar
• Direction : Philippe Boisse, Éric Maire
• Type : Thèse de doctorat
• Discipline(s) : Mécanique - Génie mécanique - Génie civil
• Date : Soutenance le 29/09/2015
• Etablissement(s) : Lyon, INSA
• Ecole(s) doctorale(s) : École doctorale Mécanique, Energétique, Génie Civil, Acoustique (Villeurbanne ; 2011-....)
• Partenaire(s) de recherche : Laboratoire : LaMCoS - Laboratoire de Mécanique des Contacts et des Structures (Lyon, INSA ; 2007-....) - Laboratoire de Mécanique des Contacts et des Structures [Villeurbanne] / LaMCoS
• Jury : Président / Présidente : Yves Rémond
• Jury : Examinateurs / Examinatrices : Philippe Boisse, Éric Maire, Yves Rémond, Christian Geindreau, Gilles Hivet, Stepan Lomov, Julien Schneider, Emmanuelle Vidal-Sallé
• Jury : Rapporteur / Rapporteuse : Christian Geindreau, Gilles Hivet

Mots clés

Mots clés contrôlés

Éléments finis, Méthode des

Armures (tissage)

Systèmes mésoscopiques

Composites à fibres

Mots clés libres

Matériau composite

Renfort tissé

Modélisation mésoscopique

Mise en forme

Hypoélasticité

Méthode des éléments finis

Microtomographie par rayon X

Woven reinforcement

Mesoscopic modelling

Shaping

Hypoelasticity

Finite element method

X-Ray microtomography

Résumé

La simulation à l'échelle mésoscopique de la déformation des renforts composites fournit des informations importantes. En particulier, elle donne la direction et la densité de fibres qui conditionne la perméabilité du renfort textile et les propriétés mécaniques du composite final. Ces analyses mésoscopiques par éléments finis dépendent fortement de la qualité de la géométrie initiale du modèle. Certains logiciels ont été développés pour décrire ces géométries de renforts composites. Mais, les géométries obtenues impliquent une simplification (notamment dans la section transversale de mèche) qui peut perturber le calcul de déformation du renfort. Le présent travail présente une méthode directe utilisant la microtomographie à rayon X pour générer des modèles éléments finis, basée sur la géométrie réelle de l'armure textile. Le modèle EF peut être obtenu pour tout type de renfort, plus ou moins complexe. Les problèmes d’interpénétrations de mèches sont évités. Ces modèles sont utilisés avec deux lois de comportement : une loi hypoélastique et une loi hyperélastique. Les propriétés de chacune d'entre elles, ainsi que les grandeurs caractéristiques nécessaires à leur implémentation dans le logiciel ABAQUS sont développées. Enfin, une identification des paramètres matériau à l'aide d'une méthode inverse est proposée. Les résultats obtenus pour les simulations de mise en forme sont comparés avec les résultats expérimentaux et montrent une bonne correspondance entre les deux.

The simulation at meso-scale of textile composite reinforcement deformation provides important information. In particular, it gives the direction and density of the fibres that condition the permeability of the textile reinforcement and the mechanical properties of the final composite. These meso FE analyses are highly dependent on the quality of the initial geometry of the model. Some software have been developed to describe composite reinforcement geometries. The obtained geometries imply simplification that can disrupt the reinforcement deformation computation. The present work presents a direct method using computed microtomography to determine finite element models based on the real geometry of the textile reinforcement. The FE model is obtained for any specificity or variability of the textile reinforcement, more or less complex. The yarns interpenetration problems are avoided. These models are used with two constitutive laws : a hypoelastic law and a hyperelastic one. An analysis of their properties is presented and their implementation in the software ABAQUS is detailed. Finally, an identification method is presented and the results of forming simulations are compared to experimental tests, which shows a good fit between the both.