Modélisation du comportement mécanique lors du procédé de mise en forme et pyrolyse des interlocks CMC

Mechanical behavior modeling of CMC interlocks through the forming and pyrolysis processes

• Auteur / Autrice : Sylvain Mathieu
• Direction : Philippe Boisse
• Type : Thèse de doctorat
• Discipline(s) : Génie mécanique
• Date : Soutenance le 09/12/2014
• Etablissement(s) : Lyon, INSA
• Ecole(s) doctorale(s) : Ecole 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 / LaMCoS
• Jury : Président / Présidente : Jean-Michel Bergheau
• Jury : Examinateurs / Examinatrices : Philippe Boisse, Jean-Michel Bergheau, Christophe Binetruy, Laurent Orgéas, Florent Bouillon, Franck Dupé, Anthony Gravouil, Nahiène Hamila
• Jury : Rapporteurs / Rapporteuses : Christophe Binetruy, Laurent Orgéas

Mots clés

Mots clés contrôlés

Pyrolyse

Composites à matrice céramique

Mots clés libres

Mécanique

Matériaux

Matériau composite

Composite à matrice céramique

Renfort tissé

Pyrolyse

Comportement mécanique

Hyperélasticité

Mechanics

Materials

Composite Material

Ceramic matrix composite

Woven fabric

Pyrolysis

Mechanical behaviour

Hyperelasticity

Résumé

La simulation des procédés de production des composites à renforts tissés est un enjeu majeur pour les industries de pointe, où leur utilisation s’intensifie. La maitrise des procédés d’obtention des composites à matrice et fibres en céramique, notamment les étapes de mise en forme et de pyrolyse, s’avère primordiale. La connaissance et la simulation du comportement mécanique aux différentes étapes est nécessaire pour optimiser les performances des pièces finales. Deux approches de modélisation macroscopique des renforts tissés épais de composite sont détaillées : une approche continue classique et une approche semi-discrète. Pour cela, une loi de comportement hyperélastique initialement orthotrope est développée. Cette loi est basée sur l’observation phénoménologique des modes de déformation privilégiés, à partir desquels sont proposés des invariants physiques de la transformation. L’identification des paramètres matériaux nécessaires est décrite. Une version modifiée de cette loi, sans contribution en tension, est implémentée dans un élément semi-discret, où le travail en tension est alors pris en compte par des barres discrétisant le tissage réel. Les importantes différences de rigidités entre sollicitations en tension et en cisaillements font des renforts tissés épais des matériaux fortement anisotropes. Leur modélisation numérique met en évidence des phénomènes parasites ou des limitations liés à cette spécificité. Le phénomène de verrouillage en tension est tout d’abord mis en évidence. Une solution basée sur une formulation éléments finis enhanced assumed strain est proposée pour des éléments continus classiques ou semi-discrets. Puis des problèmes liés aux simulations numériques dominées par la flexion sont soulevés : l’hourglassing transverse et l’absence de résistance locale à la courbure. Dans le cas de l’hourglassing transverse, deux méthodes de rigidification de ces modes de déplacement sont proposées : par moyennage des dilatations dans l’élément ou par ajout d’une rigidité matérielle tangente supplémentaire. Pour l’introduction d’une résistance à la courbure, une méthode basée sur l’utilisation purement numérique de plaques rotation free est proposée. Celles-ci permettent le calcul de la courbure induisant, par l’intermédiaire d’un moment de flexion, des efforts internes supplémentaires. Finalement, la modélisation du retour élastique après pyrolyse de la matrice organique à précurseurs céramique est réalisée. Le comportement de la matrice pyrolysée est identifié expérimentalement à l’aide d’une loi hyperélastique isotrope transverse. L’addition de cette loi, qui prend comme référence la préforme déformée, à la loi de comportement initiale du renfort tissé permet de visualiser les déformations obtenues en fin de pyrolyse. Cette modélisation est comparée à des résultats expérimentaux.

Manufacture processes modeling of woven fabrics composites is a major stake for state-of-the-art industrial parts, where their usage is intensifying. Control of all the manufacturing stages of ceramic matrix composites, particularly the forming and pyrolysis steps, is essential. Understanding and simulation of the mechanical behavior at each stage is required to optimize the final product performances. Two macroscopic modeling approaches of thick woven fabric reinforcements are detailed: a continuous classical one and a semi-discrete one. An initially orthotropic hyperelastic constitutive law is thus established. This law is based on a phenomenological observation of the main fabric deformation modes, from where physical invariants of the deformation are suggested. The required material parameters identification is explained. A modified version of this law, without any tensile energetic contribution, is implemented in a semi-discrete element where the tensile work is taken into account by bars that discretize the real weaving. Thick woven reinforcements are highly anisotropic materials due to the large ratio between the tensile rigidity and the others. Their numerical modeling highlights spurious phenomena and limitations related to this specificity. The tension locking is firstly tackled. A remedy based on an enhanced assumed strain finite element formulation is suggested for classical continuum and semi-discrete elements. Problems linked to bending-dominated numerical simulations are brought to attention : transverse hourglassing and lack of local bending stiffness. For the transverse hourglassing situation, two stiffening technics are proposed : averaging the dilatation through the whole element or adding a supplementary tangent material rigidity in a specific direction. The local bending stiffness problem is solved by calculating the curvature inside the element by using rotation free plates. The induced bending moment leads to supplementary internal loads. Finally, the elastic springback following the pyrolysis of the polymer matrix with ceramic precursors is modeled. The constitutive behavior is experimentally identified with a transverse isotropic hyperelastic law. Added to the initial reinforcements’ hyperelastic law, with the preformed fabric as reference configuration, the pyrolysis induced deformations can be visualized. This final model is compared with experimental results.