Thèse Ramirez
Thèse Ramirez
*****************************************
!
*****************************************
Tableaux ........................................................................................................................................... 15
Introduction
1.1.1 Généralités................................................................................................................... 24
3
2.3 Conclusions ......................................................................................................................... 66
3.1 Introduction @ Comportement non linéaire des matériaux composites bobinés ...... 70
4
Chapitre IV: Caractérisation expérimentale de l'endommagement de structures en
composite bobiné 0 Essais sur éprouvettes entaillées
5.2.3 Simulation des éprouvettes entaillées [02, ±452, 902]s .......................................... 186
5
5.3 Application industrielle : réservoir bobiné hyperbare ................................................. 196
Annexes
Bibliographie
6
Table de figures
8
Figure 52.a) Dépendance de la réponse en fonction de la taille de maille. b) Dépendance de
la zone en fonction de la forme et taille de maille.................................................................... 109
Figure 53. Éléments finis en série............................................................................................... 109
Figure 54. Effet de l’endommagement d’un élément sur la structure. .................................. 110
Figure 55. Influence de la vitesse sur la réponse. ..................................................................... 110
Figure 56. test de l'algorithme Box @ Muller implementé dans l'Umat.................................. 115
Figure 57. Données générées avec un facteur de forme de 4049 et d'échelle de 6.02........ 116
Figure 58. Démarche de sélection de la contrainte a rupture ou de la fraction de fibres en
un point de gauss donne. ............................................................................................................. 117
Figure 59. a) Energie absolue vs temps. b) zoom dans la zone de basse et moyenne énergie.
......................................................................................................................................................... 118
Figure 60. Méthode choisie pour représenter des signaux énergétiques............................... 119
Figure 61. Dimensions des éprouvettes d'identification. ........................................................ 120
Figure 62. Eprouvette EF [±454]s et éprouvette réelle. ........................................................... 121
Figure 63. Identification des comportements à 0º et 90º........................................................ 122
Figure 64. Identification fissuration matricielle. ....................................................................... 123
Figure 65. Courbes de traction [±454]s utilisées pour l'identification.................................... 123
Figure 66. Identification de l'évolution de l'endommagement diffus et le glissement........ 124
Figure 67. Identification des paramètres visqueux................................................................... 125
Figure 68. Comparaison test @ simulation de la traction sur [±454]s bobiné. ...................... 125
Figure 69. Résultat du test de traction sur des fibres seules. .................................................. 126
Figure 70. Comparaison probabilité cumulée expérimentale – simulée. .............................. 127
Figure 71. Assemblage des micrographies................................................................................. 128
Figure 72. Taux volumique en fonction de la taille de fenêtre utilisée.................................. 129
Figure 73. Identification du taux volumique en utilisant des lois normales. ........................ 130
Figure 74. Identification des paramètres de la loi normale et dépendance par rapport à la
taille. ................................................................................................................................................ 130
Figure 75. Energie absolue détectée dans un essai de traction sur une éprouvette entaillée
[02, ±202]s. ...................................................................................................................................... 132
Figure 76. Dimensions des éprouvettes entaillées.................................................................... 138
Figure 77. Mandrin bobiné et ligne de coupure. ...................................................................... 139
Figure 78. a) Contour des éprouvettes. b) Types de motifs obtenus selon la largeur de
mèche.............................................................................................................................................. 139
Figure 79. Dispositif expérimental des essais de traction. ...................................................... 141
9
Figure 80. Capteurs EA et marqueurs........................................................................................ 141
Figure 81. Camera utilisée pour la corrélation d'image. .......................................................... 141
Figure 82. Courbe de traction à trois températures [02, ±202]s.............................................. 143
Figure 83. Courbe de traction à trois temperatures [02, ±302]s............................................... 143
Figure 84. Courbe de traction à trois températures [02, ±452, 902]s ....................................... 143
Figure 85. Courbes de traction [±454]s À trois températures. ................................................ 145
Figure 86. a) Rupture de l'éprouvette entaillée bobinée. b) rupture d'une éprouvette
stratifiée. ......................................................................................................................................... 146
Figure 87. Courbes de traction [±67.54]s à trois températures. .............................................. 147
Figure 88. DMTA dans une éprouvette à 0º (a) et 90º (b)...................................................... 147
Figure 89. Séquence des endommagements dans une éprouvette quasi @ isotrope............. 150
Figure 90. Micrographie (vue latérale) de l'entaille d'une éprouvette [02, ±202]s.................. 150
Figure 91. micrographie des entailles ([02, ±452,902]s). ............................................................ 151
Figure 92. Etat d’endommagement de l’éprouvette [02, ±302]s. ............................................. 151
Figure 93. Radiographie X post @ mortem d'une éprouvette [02,± 202]s................................ 152
Figure 94. Zone délaminée dans une éprouvette [02, ±202]s................................................... 153
Figure 95. Radiographie éprouvette [02, ±302]s......................................................................... 153
Figure 96. Fissures intra @ laminaires dans le pli à 90º. a) entaille b) hors entaille............... 154
Figure 97. Radiographie post @ mortem d'une éprouvette quasi – isotrope. ........................ 154
Figure 98. Comparaison des radiographies à trois températures [02, ±202]s........................ 155
Figure 99. Comparaison de radiographies à trois températures des éprouvettes quasi @
isotropes. ........................................................................................................................................ 156
Figure 100. Champ de déformation de cisaillement d'une éprouvette quasi – isotrope a
différents moments du chargement. .......................................................................................... 158
Figure 101. Observations de l'endommagement dans une éprouvette du Groupe II. ....... 159
Figure 102. Radiographies post @ mortem des éprouvettes [±454]s a trois températures. .. 160
Figure 103. Corrélation d’images sur des éprouvettes [±454]s à différentes étapes du
chargement exprimées en fraction de la charge à rupture....................................................... 161
Figure 104. Radiographies post @ mortem des éprouvettes [±67.54]s a trois températures.162
Figure 105. a) Cumul d'énergie absolue dans un essai de traction [02, 202]s. b) histogramme
3D des evenements@temps@position........................................................................................... 164
Figure 106. a) Vue de dessus de l'histogramme 3D. b) Augmentation de la zone
endommagée au cours du temps. ............................................................................................... 164
10
Figure 107. a) Cumul d'énergie (vert) et courbe de traction (bleu) et b) histogramme 3D
d'une éprouvette [02, ±302]s. ........................................................................................................ 165
Figure 108. Comparaison des surfaces émettrices. .................................................................. 166
Figure 109. a) Cumul d'énergie acoustique et b) histogramme 3d d'une éprouvette [02, ±452,
902]s.................................................................................................................................................. 166
Figure 110. Comparaison des cumuls d'énergie du drapage [02, ±202]s. ............................... 167
Figure 111. Comparaison cumul d'énergie du drapage quasi @ isotrope et zoom de la partie
basse énergie. ................................................................................................................................. 168
Figure 112. a) Comportement mécanique et acoustique et b) histogramme 3d d'une
éprouvette [±454]s......................................................................................................................... 168
Figure 113. Comparaison des cumuls d'énergie de la séquence [±454]s................................ 169
Figure 114. Vue latérale d'une éprouvette [02, ±452, 02]s......................................................... 172
Figure 115. Géométrie des éprouvettes entaillées.................................................................... 177
Figure 116. a) Modèle EF coque et b) 3D d’une éprouvette [±454]s.................................... 178
Figure 117. Maillage a) coque et b) 3D...................................................................................... 179
Figure 118. Vue du maillage dans l'épaisseur du modèle 3D quasi @ isotrope. .................... 179
Figure 119. Conditions aux limites et chargement imposé. .................................................... 179
Figure 120. Courbes de traction numérique et expérimentale de l'éprouvette entaillée
[±454]s. ............................................................................................................................................ 180
Figure 121. Comparaison des charges à rupture. ..................................................................... 181
Figure 122. Comparaison entre les mesures d'émission acoustique virtuelle (sim) et réelle
(exp). ............................................................................................................................................... 182
Figure 123. Comparaison du champ de déformation axiale experimental (gauche) et
numerique (droite). ....................................................................................................................... 183
Figure 124. Niveau d'endommagement diffus dans les trois phases (a, b, c) définies sur la
courbe d'EA (Figure 122). ........................................................................................................... 184
Figure 125. a) Rupture expérimentale, b) Modèle 3D (gauche : endommagement hors plan,
droite : endommagement matriciel), c) coque (endommagement matriciel). La couleur noire
représente la rupture totale de la matrice................................................................................... 184
Figure 126. Endommagement hors @ plan a l'intérieur de l'éprouvette juste avant la rupture
finale (Noir = Totalement endommagé). .................................................................................. 185
Figure 127. Comparaison expérience @ Simulation 2d @ 3D (déterministe). ......................... 186
Figure 128. Charges à rupture expérimentale et simulées de la séquence quasi @ isotrope. 187
Figure 129. Simulations probabilistes du drapage quasi @ isotrope........................................ 187
11
Figure 130. Charges à rupture (simulations probabilistes)...................................................... 188
Figure 131. Comparaison énergie cumulée expérimentale et simulée (séquence quasi @
isotrope).......................................................................................................................................... 188
Figure 132. Comparaison fissuration matricielle post @ mortem expérimentale et numérique
(noir = complètement endommagé). ......................................................................................... 189
Figure 133. Comparaison des cartographies numérique et expérimentale de
l'endommagement hors – plan a l’interface de différents plis. (noir = totalement
endommagé). ................................................................................................................................. 190
Figure 134. Déformée (grandissement x5) du centre de l'eprouvette entaillee quasi @
isotrope........................................................................................................................................... 191
Figure 135. Composite bobiné ondulé (a) modèle EF et b) micrographie) ......................... 192
Figure 136. Endommagement hors plan (en haut) pour différentes valeurs de contraintes
dans le sens des fibres (en bas).................................................................................................... 192
Figure 137. Comparaison de la cartographie de la fissuration matricielle avec une
micrographie post – mortem....................................................................................................... 193
Figure 138. Comparaison de la Cartographie de l'endommagement hors – plan avec une
micrographie post mortem. ......................................................................................................... 194
Figure 139. Comparaison de la cartographie de l'endommagement diffus avec une
micrographie post @ mortem........................................................................................................ 194
Figure 140. Champ de glissement plan γ12 a) expérimental et b) numérique, pour une valeur
de force axiale appliquée de 30% de la charge maximale dans les deux cas......................... 195
Figure 141. Réservoir bobiné de type IV................................................................................... 197
Figure 142. Fabrication du réservoir bobiné de type IV. ........................................................ 198
Figure 143. Géométrie du liner et de l’embase......................................................................... 199
Figure 144. a) Alternance d'orientations dans l'enroulement, b) Paramètres de calcul de
l'évolution des angles.................................................................................................................... 200
Figure 145. Comparaison des dômes EF et réels..................................................................... 202
Figure 146. Modèles 3D (droite) et axisymétrique (gauche)................................................... 202
Figure 147. Comparaison des masses des réservoirs numérique et réel. .............................. 203
Figure 148. Exemple de la variation des propriétés dans une structure bobinée, en fonction
de la variation d’angle dans le dôme........................................................................................... 204
Figure 149. Evolution des angles dans le réservoir.................................................................. 204
Figure 150. a) courbe de traction a rupture de l'acier 316l et b) identification de la loi
plastique.......................................................................................................................................... 205
12
Figure 151. Conditions aux limites et déplacements mesurés. ............................................... 206
Figure 152. Courbes de déplacement expérimentale et simulée du réservoir de type IV à
20ºC................................................................................................................................................. 207
Figure 153. Courbes de déplacement expérimentale et simulée du réservoir de type IV à
20ºC avec le déplacement axial corrigé...................................................................................... 208
Figure 154. a) Réservoir éclaté en mode sain et b) non sain. ................................................. 209
Figure 155. Déformation plastique équivalente dans l'embase à 600 bar............................. 209
Figure 156. Radiographie du réservoir. Détail du jeu existant. .............................................. 210
Figure 157. Simulation avec un jeu entre le composite et l'embase dans les dômes (le
comportement du composite est considéré ici comme élastique sans endommagement). 210
Figure 158. Endommagement des fibres (une valeur de 1 représente des fibres rompues)
(1580 bar, 1600 bar et 1620 bar). ................................................................................................ 211
Figure 159. Emission acoustique pendant l’essai d'éclatement, (temps en s en abscisse, 2000
seconds = 650 bar). ...................................................................................................................... 212
Figure 160. Rupture matricielle et hors @ plan a 600 bar. ........................................................ 212
Figure 161. a) Détail de l’endommagement hors @ plan et b) fuite détectée lors des essais.
......................................................................................................................................................... 214
Figure 162. Influence de la rupture matricielle dans le comportement du réservoir. ......... 215
Figure 163. Comparaison des déformations de certaines couches (en rouge, à droite,
couches circonférentielles)........................................................................................................... 216
Figure 164. Comparaison entre simulations probabiliste et deterministe de l'éclatement du
réservoir.......................................................................................................................................... 216
Figure 165. Comparaison entre deux simulations probabilistes 3D et axisymétrique. ....... 217
Figure 166. Endommagement des fibres de la première couche circonférentielle @simulation
probabiliste 3D (noir : fibre rompue). ....................................................................................... 218
Figure 167. a) Cellule élémentaire, b) Chargements élémentaires pour calculer l'influence
sur la rigidité (Huchette 2005)..................................................................................................... 223
Figure 168. Eprouvette testée. .................................................................................................... 224
Figure 169. Modèle éléments finis de l'essai de flexion trois points...................................... 224
Figure 170. Comparaison des courbes force @ déplacement expérimentale et numérique et
configuration de l’endommagement hors@plan (délaminage). ................................................ 224
Figure 171. Délaminages observés pendant l'essai de flexion. ............................................... 225
Figure 172. Evolution de l’endommagement par rapport au nombre de cycles pour trois
amplitudes différentes. ................................................................................................................. 228
13
Figure 173. Comparaison expérimentale et numérique de la courbe contrainte –
déformation au cours du cyclage. ............................................................................................... 228
Figure 174. Evolution du cumul dans un cycle de charge@décharge. .................................... 230
Figure 175. Influence du paramètre n dans l'évolution des variables d'endommagement (d).
......................................................................................................................................................... 231
Figure 176. Evolution expérimentale et simulée de la tenue résiduelle. ............................... 232
Figure 177. Chargements cYcliques appliqués. ........................................................................ 232
Figure 178. Evolution de l'endommagement selon le chargement appliqué........................ 233
Figure 179. Simulation cyclage et traction................................................................................. 234
Figure 180. Évolution de l’endommagement matriciel cyclique. ........................................... 235
Figure 181. a) géométrie du dôme généré et b) détail de l'assignation des angles............... 236
Figure 182. Modes d'éclatement a) sain et b) non sain............................................................ 249
Figure 183. a) Dôme du modèle 1 et b) modèle 2. .................................................................. 250
Figure 184. Eclatement des modèles 1 et 2............................................................................... 251
Figure 185. Rupture des fibres (rouge=totalement endommagé) pour a) modèle 1 et b)
modèle 2. ........................................................................................................................................ 251
Figure 186. Endommagement matriciel dans le moment de l'éclatement............................ 252
14
Tableaux
15
INTRODUCTION
Les matériaux composites à fibres longues sont de plus en plus une vraie alternative
aux matériaux métalliques grâce aux propriétés mécaniques spécifiques, l’adaptabilité aux
besoins mécaniques de la structure ou de la pièce à fabriquer et leur faible sensibilité aux
phénomènes de fatigue. Le domaine d’utilisation des composites est très vaste, puisqu’il
s’étend de l’industrie aérospatiale au matériel sportif. C’est pour ces raisons que l’utilisation
de ce type de matériau est en pleine croissance depuis nombreuses années.
Une des utilisations des composites à fibre longue en plein essor depuis plusieurs
années est la fabrication des structures cylindriques (dont les réservoirs). La fabrication de
ces structures de révolution est faite par bobinage de fibres (pré – imprégnées, voie
humide, fibres sèches…) sur un moule ou un mandrin. Cette procédure est qualifiée
d’enroulement filamentaire. Ce processus entraine une architecture finale du matériau
différente des stratifiées utilisés dans l’aérospatiale. Cette architecture particulière aura des
conséquences sur le comportement mécanique du matériau, notamment sur les modes et la
cinétique de dégradation.
Une des applications de ces structures bobinées est la fabrication des réservoirs pour le
stockage de gaz à haute pression. Plus particulièrement, dans le cadre de cette thèse, nous
étudierons des réservoirs conçus pour le stockage de l’hydrogène à très haute pression et
qui feront partie des voitures de l’avenir équipées de piles à combustible, en réponse aux
problématiques énergétiques et écologiques.
17
développement des technologies de l’hydrogène permettrait de réduire de 50% les
émissions de C02 et de 40% la consommation de pétrole d’ici 2050.
La filière hydrogène apparaît ainsi comme la voie la plus prometteuse pour répondre à
la problématique énergétique mondiale et à l’évolution du climat de la planète. Le stockage
hyperbare de l’hydrogène gazeux est aujourd’hui la voie la plus mature. Pour structurer les
réservoirs à des pressions élevées (en tenant compte des facteurs de sécurité) sans obtenir
des masses de système de stockage trop importantes, des matériaux composites hautes
performances sont utilisés. Alors que les structures composites minces représentent une
très grosse majorité des cas, les pressions mises en jeu dans les stockages d’hydrogène
conduisent à de fortes épaisseurs. Les connaissances associées au dimensionnement et à la
conception de structures épaisses sont plus limitées alors que, de fait, celles@ci sont
nettement plus sensibles aux variabilités de procédés et/ou de matériaux. Pour y faire face,
les industriels accroissent les quantités de fibres de carbone déposées sur les réservoirs au
détriment du prix et de la masse du réservoir. Cela constitue un véritable problème car la
fibre de carbone représente déjà 45% à 80% du prix du réservoir selon la fibre de carbone
utilisée. Le prix demeure un des verrous majeurs de cette technologie inhibant son
déploiement. La conception des réservoirs fait en général appel à l’utilisation de codes de
calcul plus ou moins évolués. Les codes de calcul utilisant les éléments finis pour le design
de l’architecture fibreuse comme ABAQUS, ANSYS, SAMCEF, SYSPLY, ZEBULON…
sont généralement employés en pré@dimensionnement. Le dimensionnement du réservoir
n’est généralement finalisé qu’à l’issue de phases itératives de fabrications suivies d’essais
destructifs (processus long et très coûteux). Les éléments finis utilisés pour les calculs sont
le plus souvent des éléments de type « coque ». Ceux@ci ne sont pas très bien adaptés pour
décrire le comportement en cas de fortes épaisseurs. Les éléments 3D sont très peu utilisés
car leur mise en œuvre est le plus souvent très laborieuse. En outre, pour la phase de
dimensionnement, les critères de rupture comme les critères énergétiques (TSAI@WU,
TSAI@HILL, HASHIN…) sont en général trop conservatifs et, en conséquence, non
utilisés car conduisant eux aussi à des surdimensionnements notables. Ils ne tiennent
18
notamment pas compte des processus le plus souvent cumulatif de dégradation de
l’amorçage jusqu’à la rupture de l’objet. Lorsque des modèles d’endommagement sont
utilisés, ils sont le plus souvent spécifiques à un type de matériau, à un type de sollicitation
ou à un type d’architecture. Par ailleurs, les structures simulées ne sont pas les structures
réelles car le procédé d’enroulement filamentaire provoque de nombreux défauts et les
entrecroisements de torons de fibres ou encore les variations réelles des angles et des
épaisseurs, notamment dans les dômes des réservoirs, ne sont pas ou mal connus. Il en
résulte qu’aujourd’hui, les outils de calcul par éléments finis pour la conception des
réservoirs et les méthodologies associées ne sont pas qualifiés et chaque utilisateur utilise
son code avec ses propres règles.
Le projet ANR auquel cette thèse a contribué, OSIRHYS IV, tente d’apporter des
réponses à ces verrous scientifiques. Il se penche sur le stockage haute pression de
l’hydrogène pour des réservoirs de type IV où le liner, assurant uniquement la fonction
d’étanchéité à l’hydrogène, est en matériau polymère et la fonction de résistance mécanique
est prise en charge par la coque composite (très épaisse). L’objet de ce projet est de faire un
état des lieux des codes de calculs EF pour le design et la simulation du comportement des
réservoirs de stockage hyperbare de l’hydrogène gazeux de type IV obtenus par
enroulement filamentaire. Le travail de thèse présenté ici trouve son application en
contribuant au projet OSIRHYS IV et en répondant plus particulièrement au problème de
la simulation de l’endommagement, et de la rupture, dans les structures épaisses en
composites bobinés.
Dans le premier chapitre, les matériaux composites seront présentés d’une façon très
générale en insistant sur les composites à fibre de carbone – matrice polymère bobiné. On
fera un bref rappel des particularités mécaniques de ce type de matériaux comme
l’anisotropie mécanique. Dans une deuxième partie, on présentera les différents outils
expérimentaux, comme l’émission acoustique et la corrélation numérique d’image qui
seront utilisés dans le cadre de cette thèse.
19
Le deuxième chapitre est consacré à la description des différentes approches de la
modélisation de l’endommagement des matériaux composites. On notera que la grande
majorité de ces outils sont dédiés aux stratifiés. On fera ainsi le point sur des modèles écrits
à l’échelle des microconstituants (« modèles micromécaniques »), à l’échelle du pli
unidirectionnel (« mésoscopiques ») comme le mésomodèle de Cachan ou le modèle
développé à l’ONERA. Des modèles de rupture progressive seront aussi présentés car ils
ont été les premières tentatives de simulation de la rupture de structures composites. On
finira ce chapitre en introduisant une nouvelle famille de modèles qui essaient d’incorporer
le caractère probabiliste de la rupture des fibres dans le comportement du composite.
Parmi tous ces modèles, seulement le dernier a été utilisé dans les composites bobinés, ce
qui nous donne une idée du manque d’études consacrées à ce type de matériau par rapport
aux stratifiés.
20
trois températures (@40ºC, 20ºC et 85ºC) pour nous permettre de mieux comprendre le
comportement de ce type de matériau en présence de concentration de contraintes et les
différents modes d’endommagement (amorçage et évolution). Tous ces essais ont été multi
– instrumentés (émission acoustique, micrographies, rayons X, DIC et suivi optique)
permettant l’étude des éprouvettes de différents points de vue, toujours complémentaires.
Cette campagne d’essais nous permettra aussi de construire une base de données pour
valider le comportement de notre modèle d’endommagement en simulant le comportement
des séquences les plus représentatives.
Cette validation du modèle ainsi que son application à la tenue des réservoirs bobinés
seront traités dans le chapitre 5. Ce chapitre est divisé en deux parties bien distinctes : dans
la première, on fera une comparaison entre les essais menés sur certaines éprouvettes
entaillées et les résultats fournis par le modèle. On essaiera de vérifier les capacités de
prédiction multi – échelle du modèle, c’est@à@dire de voir s’il est capable de reproduire
correctement le comportement mécanique macroscopique (courbes force – déplacement)
mais aussi le type d’endommagement présent dans le matériau ainsi que sa localisation et
son étendue. Une fois le modèle validé, on testera ses capacités dans une application
industrielle : la prédiction de la tenue d’un réservoir bobiné. La géométrie de cette structure
sera modélisée à l’aide du plugin Abaqus Wound Composite Modeler. On analysera les
prédictions d’éclatement, dans un premier temps déterministes. La richesse du modèle nous
permettra de prédire des phénomènes complexes comme la création des chemins de fuite
ou le mode d’éclatement (éjection ou non des extrémités). La prise en compte de la rupture
probabiliste apportera une amélioration significative dans la prédiction de l’éclatement
grâce à l’effet de redistribution des contraintes.
21
CHAPITRE I:
1.1.1 Généralités.
Les composites sont composés d’un ensemble de matériaux différents avec des
propriétés également différentes. La plupart du temps, il s’agit de renforts noyés dans une
matrice qui doit assurer la cohésion de l’ensemble, répartir les efforts et protéger les
renforts contre l’environnement. Les renforts ont une fonction purement structurelle, c'est@
à@dire qu’ils doivent résister aux contraintes imposées. L’avantage de ce type de matériaux
est la capacité de s’adapter à l’état de contraintes qu’il va supporter. De cette façon, on peut
obtenir un matériau optimisé pour la fonction qu’il va réaliser. En plus, ses propriétés
spécifiques sont notamment meilleures par rapport aux matériaux classiques comme les
métaux. Ces deux caractéristiques ont fait que l’utilisation des matériaux composites est en
pleine croissance dans le monde industriel en général et aéronautique en particulier (Figure
1 et Figure 2).
[1]
Figure 1. Evolution de l'utilisation des composites dans l'aéronautique.
24
(, ! -
[1]
Figure 2. Utilisation des composites CFRP dans le domaine industriel.
Tº Transition
E
Type Matériaux vitreuse
(g/cm3) (GPa) (GPa) (Cº)
. 1.2 2.8@3.5 0.05@0.08 60@100
Thermodurcissables
!. 1.1@1.5 3@5 0.1@0.15 290
. . 0.9 1.1@1.4 0.02@0.035 50@60
Thermoplastiques
. 1.14 1.2@2.5 0.06@0.08 65@100
/ 1.23 2.7@3.5 0.045 230
Thermostables 0.069@
. 1.43 3.2 360
0.089
Tableau 1. Propriétés mécaniques de quelques matrices organiques
Les renforts, quant à eux, sont souvent de nature filamentaire ou particulaire (Figure 3).
25
(, ! -
Les propriétés mécaniques de ces fibres diffèrent beaucoup selon la nature de la fibre.
On peut voir dans le graphique suivant l’ensemble de types des fibres selon leur contrainte
à rupture et leur module spécifique.
[1]
Figure 4. Contrainte à rupture et module spécifiques pour les différents types de fibres longues.
On peut voir dans la Figure 4 que les fibres de carbone ont les meilleures
performances. On ne trouve pas dans ce graphique les fibres d’origine végétale, comme le
lin ou le chanvre qui ont un faible coût mais aussi des propriétés mécaniques faibles. Son
utilisation est plutôt orientée vers l’industrie de l’automobile. Dans le cas des fibres longues,
qui sera traité dans la suite, les paramètres influant sur les propriétés sont notamment la
géométrie des fibres, la concentration des constituants, l’orientation du renfort.
26
(, ! -
Les multiples combinaisons existant entre les matrices et fibres qu’on vient de
présenter donnent lieu à une grande quantité de composites différents, par exemple :
• CFRP : fibres de carbone dans une matrice polymère. Très utilisés dans l’industrie
aérospatiale, l’énergie...
• GFRP : les fibres de carbone sont remplacées par des fibres de verre.
• CMC : composites à matrice céramique parmi lesquels on peut trouver des fibres de
carbone ou de SiC dans une matrice également SiC, C ou Al2O3. Ces composites
sont très utilisés dans les applications à très haute température.
• « Green Composites » : fibres végétales dans une matrice polymère.
La procédure d’élaboration liée à cette architecture est le « lay – up » qui peut être
manuel (« Hand Lay – up ») ou automatique. Cette procédure est assez couteuse,
notamment en raison de sa durée et le besoin d’un autoclave pour la cuisson de la pièce.
Néanmoins, les propriétés mécaniques obtenues sont très élevées, si bien que cette
procédure est la plus utilisée dans les applications visant les plus hautes performances
(aéronautique). Si, à la place des plis UD, on utilise des tissus, une nouvelle gamme de
processus de fabrication apparaît, comme le RTM ou l’infusion. Ces procédures permettent
de fabriquer des pièces à une vitesse supérieure ou de très grandes tailles (notamment
l’infusion).
27
(, ! -
Un autre procédé d’élaboration, très utilisé dans le cas des pièces possédant une
symétrie de révolution, et qui constituera l’objet d’étude de cette thèse, est l’enroulement
filamentaire : une machine vient déposer une mèche de fibres de carbone pré – imprégnée
sur un moule ou, par exemple, sur le liner polymère du réservoir de type IV. On trouve
plusieurs façons de déposer cette mèche [2] sur la structure comme on peut le voir sur la
figure suivante :
• le type d’enroulement.
• le taux d’imprégnation des fibres.
• la tension de la mèche.
• le type de trajectoire (dans le cas des réservoirs).
• la durée de la cuisson
• le chevauchement des fibres
L’architecture du composite est déterminée par l’angle formé par la mèche et l’axe de
symétrie. Un exemple de motif résultant de cette fabrication est indiqué sur la figure
suivante [3]. On notera l’entrecroisement des mèches qui sera responsable de certaines
particularités de comportement, comme on pourra le voir dans la suite.
28
(, ! -
Ce motif peut nous rappeler l’architecture d’un composite tissé. Le composite bobiné
est composé de cellules de forme et de taille variables selon les paramètres de fabrication.
Ces variations géométriques ont des conséquences sur le comportement mécanique. On
voit dans la figure suivante un exemple du bobinage à ±45º avec différents paramètres.
a) b)
Figure 9. a) Enroulement d'un tube, b) installation composée de structures bobinées.
29
(, ! -
Ces réservoirs peuvent être très différents par leur taille selon leur utilisation
(embarqué ou fixe). Une caractéristique spécifique qui aura une grande influence sur le
comportement de ces structures est la forte épaisseur de la coque composite. En effet, cette
forte épaisseur limitera la validité des techniques de conception et de simulation du
comportement mécanique (théorie des stratifiés, en particulier) comme on peut le voir dans
le graphique ci – dessous.
[4]
Figure 11. Evolution de la pression d'éclatement calculée pour un tube [90, ±20,90] selon un
modèle d'épaisseur fin ou épais (b=rayon extérieur, a=rayon intérieur).
La Figure 11 montre comment un modèle basé sur la théorie de stratifiés prédit une
évolution linéaire de la pression d’éclatement qui s’éloigne de la prédiction d’un modèle
épais. En plus, cette forte épaisseur provoquera des problèmes pendant la cuisson en raison
du gradient de température. Ce fait génère des contraintes résiduelles d’origine thermique
qui peuvent être significatives [5,6]. Cet effet d’épaisseur nous amènera dans le cadre de
cette thèse à formuler un modèle 3D pour essayer de mieux reproduire le comportement
mécanique de la coque composite.
Finalement on peut aussi trouver des structures de très grande taille fabriquées par
cette méthode dans le domaine aérospatiale, notamment les structures et les réservoirs des
fusées (Figure 12a) ou des parties d’aéronefs comme le fuselage (Figure 12b).
30
(, ! -
a) b)
Figure 12. a) Structure de fusée en composite bobiné, b) fuselage fabriqué par enroulement
filamentaire.
Une des caractéristiques mécaniques les plus importantes des composites est
l’anisotropie naturelle présente dans ce type de matériau. À la différence des matériaux
isotropes, lorsqu’on impose un champ de contraintes, les axes principaux du matériau ne
sont pas confondus avec les directions des contraintes principales.
Les composites à fibres unidirectionnelles (par exemple l’échelle du pli dans les
matériaux stratifiés) font partie la catégorie des matériaux à isotropie transverse, c'est@à@dire
que le matériau présente une isotropie dans le plan orthogonal à la direction privilégiée
(direction des fibres). La matrice de souplesse dans ce cas a la forme suivante:
31
(, ! -
1 -ν21 -ν21
E 0 0 0
E2 E2
1
-ν12 1 -ν
0 0 0
ε 11 E 1 E2 E2 σ 11
ε -ν -ν 1 σ
22 12 0 0 0 22
ε 33 E 1 E2 E2 σ 33
= σ 23
γ 23 0 2 (1+ ν )
0 0 0 0
γ 13 E2 σ 13
σ
γ 12 0 0 0 0
1
0 12
G 12
1
0 0 0 0 0
G 12
modules de Young associés à la direction des fibres (1) et transverse (2), le module de
Les travaux présentés dans cette thèse viseront à établir, à valider et à utiliser un
modèle d’endommagement pour prévoir la rupture d’une structure complexe (réservoir
bobiné). On suppose que les mécanismes d’endommagement présents dans le composite
bobiné seront très similaires (voir identiques) à ceux qu’on trouve dans les composites
stratifiés. C’est la raison pour laquelle on analysera dans cette section les différents modes
d’endommagement des stratifiés dans le but de montrer quel mécanisme nous devrons
prendre en compte dans la future modélisation du composite bobiné. On doit ajouter que,
même s’il existe une vaste base de données pour les composites stratifiés, les composites
bobinés n’ont pas été étudiés d’une manière si approfondie.
32
(, ! -
Interfaces
Fibre Matrice
Inter 0 pli Fibre0 matrice
Micro – fissuration
Décohésion fibre @
Rupture des fibres Micro – vides Délamination
matrice
Macro fissuration
Tableau 2. Types d’endommagement dans les composites stratifiés
Ces phénomènes sont schématisés sur la Figure 13. De nombreux facteurs ont une
influence sur l’apparition puis l’évolution des différents mécanismes d’endommagement.
Ces facteurs peuvent être le type de sollicitation (statique/fatigue), la direction du
chargement, le type d’empilement, la température, l’humidité…. La suite de cette section
passe en revue d’une manière plus détaillée les différents modes d’endommagement.
33
(, ! -
Figure 14. Fissure matricielle formee à partir de la coalescence des décohésions et microfissures
matricielles.
Les causes de ces décohésions sont variées. On peut citer, par exemple, la présence de
micro@défauts au niveau de l’interface fibre @ matrice (comme des micro@vides) qui entraîne
l’amorçage des décohésions, mais aussi les contraintes d’origine thermique (pendant la
cuisson) ou l’état de contrainte local induit par la répartition aléatoire des fibres [8].
34
(, ! -
1.2.3 Le délaminage.
35
(, ! -
[22] [23]
Figure 17. Micro délaminage en extrémité de fissure matricielle.
36
(, ! -
et de décohésions dans les zones riches en matrice présentes dans l’interface inter – pli. On
peut voir dans la Figure 18.
[26]
37
(, ! -
[29]
Figure 19. Délaminage observe dans une cornière a) sans fissuration matricielle et b) avec
fissuration.
Malgré tout, le délaminage est encore un des endommagements qui fait l’objet de
nombreuses études et qui pose des questions tant au niveau expérimental que numérique.
38
(, ! -
L’émission acoustique est une technique qui nous permet d’entendre des sons créés
par des événements se produisant dans les matériaux (notamment des endommagements).
Par exemple, lorsqu’une fissure est créée, une certaine énergie est libérée sous la forme
d’une onde élastique [38]. Cette technique nous permet d’enregistrer ces ondes en temps
réel. Pour capter ces ondes, on utilise deux (ou plus) capteurs piézoélectriques qui
convertissent l’énergie mécanique en signaux électriques. Un gel couplant est utilisé pour
assurer la transmission des ondes entre le matériau et le capteur piézoélectrique. Un
ordinateur pourvu d’une carte d’acquisition spécifique est chargé d’enregistrer et de traiter
les données via un logiciel déterminé (AEWin dans notre cas). Le schéma de la Figure 20
illustre le principe.
Cette technique est très utilisée pour le suivi d’endommagement via l’activité
acoustique. Certains auteurs l’ont utilisée pour étudier le mouvement de dislocations dans
les matériaux métalliques [39], la rupture des fibres [40], l’apparition des fissures dans l’acier
A316 [41] etc. Lorsqu’on étudie l’activité acoustique d’un matériau pendant un essai, deux
phénomènes différents peuvent apparaître. D’un côté, l’effet Kaiser désigne une émission
de signaux acoustiques une fois dépassée la charge atteinte précédemment (Figure 21
ABC), par exemple quand on applique des cycles croissants de charge. Sous certaines
conditions, on détecte une activité acoustique avant de dépasser la charge maximale (Figure
21 DEF). On parle alors d’effet Felicity, pour lequel on définit le rapport suivant :
avec
Équation 1
39
(, ! -
où est la force pour laquelle on détecte des échos et est la force maximale
L’utilisation d’un unique capteur permet d’entendre les signatures acoustiques. Avec
deux capteurs (ou plus), il est possible de localiser (dans une direction) la source de l’écho.
Des tests acoustiques préliminaires sont indispensables pour calculer la vitesse de
propagation d’onde dans le matériau. Une fois connue cette vitesse, chaque capteur
entendra l’écho à un instant différent, soit et . Ces différences de temps et la
connaissance de la vitesse de propagation des ondes permettent de repérer la source
d’émission acoustique. On peut utiliser plus de deux capteurs soit pour améliorer la
précision dans une localisation unidimensionnelle, soit pour être capable de localiser dans le
plan [44].
40
(, ! -
dans [48] et [49]. Les classificateurs non supervisés, contrairement aux premiers, n’ont pas
besoin d’entraînement pour classifier les échos. Une méthode non supervisé est la méthode
des k – moyennes utilisée récemment par Bertin [50].
La radiographie par rayons X est une technique qui nous permet d’obtenir des images
dont le contraste dépend de la densité du matériau et de son épaisseur. Grâce à cette
technique, on observe une projection de certains endommagements comme les fissures
matricielles (des lignes droites) ou le délaminage sous la forme de taches diffuses. Le
fonctionnement de cette technique est assez simple. D’un côté, nous avons le foyer de
rayons X, de l’autre une surface sensible aux rayons X et entre les deux on place
l’échantillon. Les différents endommagements subis par l’échantillon créent des espaces
« vides », mais ce vide n’ajoute pas assez de contraste pour pouvoir être observé. C’est
pourquoi on doit traiter l’échantillon avec un opacifiant qu’on vient déposer sur la surface
libre et qui remplit, par capillarité, ces espaces vides. Cet opacifiant est composé de 60g
d’iodure de zinc, 10 ml d’eau distillé, 10 ml d’alcool isopropylique et quelques gouttes de
mouillant photoflow.
Rayons X
Echantillon
Image
a) b)
Figure 22. a) Schema de projection des endommagements sur le plan. b) Image obtenu par rayons
X [51]
41
(, ! -
Image 1 Image 2
Une fois que les fenêtres sont identifiées, on mesure le déplacement existant et on
obtient le champ de déplacement complet. Il nous reste à dériver ce champ de
42
(, ! -
déplacement pour obtenir le champ de déformation. Pour corréler les deux images et
obtenir les champs de déplacements, on utilise un logiciel développé à l’Institut P’ depuis
quelques années appelé Granu. Une fois obtenus ces déplacements, on crée un modèle
éléments finis sous ABAQUS qui représente la même zone qui été photographiée. Le
champ de déplacement est donc imposé comme chargement et on obtient directement les
champs de déformation correspondante.
43
CHAPITRE II:
1% "1% #* "
*##'+ # $
#'$ &('%) *# * ($
(( , 1 2 - !
2.1 Introduction.
Nous présentons par la suite quelques exemples des différentes approches et leur
apport à l’étude de l’endommagement et de la prédiction de la rupture des structures
composites. Une des différences principales parmi les différents modèles réside dans
l’échelle à laquelle on se place pour construire le modèle. De cette façon, on trouve des
modèles qui envisagent l’échelle de la fibre et de la matrice, et qui prennent en compte le
comportement de la fibre d’une part, celui de la matrice d’autre part et leur interaction. De
cette manière, nous obtenons le comportement du pli à partir de ses éléments constituants.
Si on s’éloigne et on prend le pli unidirectionnel comme unité élémentaire, on parle de
modèles à l’échelle mésoscopique. Le modèle mésoscopique modélise le comportement du
pli UD, puis de structures de plus grande taille. Par contre, se placer à cette échelle oblige à
faire un effort conséquent de modélisation du comportement du pli. D’autres modèles sont
construits sur une approche hybride : c’est l’approche multi – échelle (changement d’échelle
selon les besoins).
47
(( , 1 2 - !
La prévision de la rupture peut être traitée par les modèles d’endommagement de deux
manières différentes : les approches de type « rupture progressive » [67,68] utilisent des
critères de rupture qui, une fois atteints, dégradent localement les propriétés dans certaines
directions, alors que d’autres modèles décrivent la dégradation progressive [69] des
propriétés et associent la rupture à la détection d’une divergence du comportement ou de la
simulation.
Ajoutons qu’il est possible d’enrichir les modèles cités précédemment de manière à
reproduire la nature statistique de certains phénomènes comme la rupture des fibres [70].
On dispose alors d’outils aptes à interpréter des effets de type « maillon faible », de tailles,
de gradients etc…
48
(( , 1 2 - !
Figure 26. Endommagements échelle micro 1)Macro fissure matricielle, 2)Micro délaminage en
pointe de fissure, 3)Endommagement diffus (microfissuration matricielle, décohésion fibre –
matrice) et 4) Macro délaminage d’après [71].
49
(( , 1 2 - !
Les critères d’initiation et de propagation dans ce type d’approche font intervenir les
forces thermodynamiques à l’échelle locale associées aux modes de rupture I, II et III,
l’épaisseur du pli et l’épaisseur seuil (voir Figure 27). En ce qui concerne le délaminage, le
modèle utilise aussi des surfaces minimales de rupture dont l’initiation et la propagation
dépendent des taux de restitution d’énergie locaux. La rupture de fibres est, elle, considérée
comme un phénomène fragile. Comme dans le cas des mécanismes « continus », on adopte
un volume cubique minimal élémentaire de hauteur pour lequel on définit une force
thermodynamique moyenne en traction, une force thermodynamique moyenne en
compression et une force thermodynamique de cisaillement. Avec ces trois éléments, on
peut construire deux critères différents, un pour une rupture en traction et l’autre pour une
rupture en compression.
50
(( , 1 2 - !
Équation 2
Équation 3
de Poisson liés aux trois directions du matériau, à savoir, 1 pour la fibre, 2 pour la direction
perpendiculaire aux fibres (dans le plan) et 3 pour la direction hors – plan. Les termes
et prennent la valeur 1 dans les cas de chargement positif et 0 autrement.
• : rupture de fibre
51
(( , 1 2 - !
• , , : fissuration transverse
Équation 4
Équation 5
Pour différencier les ruptures en traction et compression, deux valeurs seuils sont
utilisées : et . Si la force est supérieure à et alors la rupture en traction
est atteinte. Dans le cas où est supérieure à mais , la rupture atteint sera celle
Équation 6
Équation 7
52
(( , 1 2 - !
expérimentalement.
être mesurée, par exemple, par comptage de microfissures). La densité peut être liée
suivante :
Équation 8
les taux de restitution d’énergie critiques pour les modes I, II et III respectivement, on
Équation 9
• Si
53
(( , 1 2 - !
Équation 10
avec
Équation 11
Équation 12
• Sinon :
Équation 13
de ces variables est une combinaison des forces associées à chacune des variables. Le
schéma utilisé est très similaire à celui de la fissuration matricielle, en distinguant le cas de la
densité de fissures inférieure ou égale à la valeur de saturation. Les évolutions de ces
• Si
Équation 14
54
(( , 1 2 - !
Équation 15
Équation 16
• Autrement, si
Équation 17
Ce modèle a été utilisé largement pour le calcul de structures stratifiés, notamment des
éprouvettes trouées, traitées dans l’article [81], dans lequel le modèle d’endommagement
intra pli est implémenté dans un code industriel comme ABAQUS, ou dans la référence
[79] où le modèle est utilisé dans sa version la plus complète pour prédire le type de rupture
dans des éprouvettes trouées selon le type de séquence utilisé (groupé ou non groupé). Le
premier article aborde une simulation de la traction sur une éprouvette trouée [02, 902]s et
en montre les cartographies d’endommagement.
Le deuxième article poursuit l’analyse des éprouvettes trouées avec des drapages
différents ([45n/90n/@45n/0n]m s, les coefficients n et m prenant les valeurs 1,2 et 4). Un
exemple de résultat est donné sur la figure suivante :
b)
55
(( , 1 2 - !
a)
Figure 28. Courbes de traction pour différentes valeurs de n et m et cartographie pour un
drapage m=N=2. [79]
Ces courbes montrent l’effet de l’épaisseur des plis (groupé ou non groupé) sur le
comportement qui est en accord avec des données expérimentales montrant l’effet
bénéfique des plis fins et non groupés [50]. On peut voir sur la Figure 28b les
cartographies d’endommagement des fibres ( ), de densité de fissuration ( ) et délaminage
( ). Ce modèle, dans une version simplifiée sans endommagement d’interface, a été utilisé
par [82] pour modéliser des tissus de fibres de verre et époxy en statique et aussi en fatigue
comme on en parlera ultérieurement. Cependant, on n’a pas constaté d’utilisation de ce
modèle dans le cas spécifique des composites bobinés de forte épaisseur.
Cette approche est basée sur le calcul des variables à différentes échelles du matériau
composite. Une des approches multi@échelles les plus connues est celle construite par [85]
pour être appliquée aux stratifiés carbone – époxy. Ce modèle opère une transition
d’échelle entre l’échelle mésoscopique (celle correspondant au pli unidirectionnel) et
l’échelle dit macroscopique correspondant au stratifié complet. Le schéma général est
présenté ci@dessous.
56
(( , 1 2 - !
Cette approche peut être divisée en cinq étapes différentes : une méthode de
changement d’échelle, le choix du comportement mésoscopique, la définition d’un critère
de rupture mésoscopique, la dégradation des propriétés et la rupture du stratifié. La
méthode de changement d’échelle repose sur la Théorie classique des stratifiés. Le critère
de rupture est basé sur le multi – critère de Hashin [58]. Il envisage un mode de rupture de
fibre et un autre dit d’inter – fibre, tous les deux considérés soit en traction, soit en
compression. Le critère inter – fibre ( =1 en traction ou =1 en compression) a
l’expression suivante :
Équation 18
Autrement, si
Équation 19
Les termes qui apparaissent dans les deux expressions, , et , représentent les
transversale. L’introduction de ce terme permet d’avoir une plus grande résistance dans un
état de chargement en cisaillement et en légère compression. Les résistances effectives sont
calculées de la façon suivante :
57
(( , 1 2 - !
Équation 20
Une autre amélioration par rapport au critère classique de Hashin est l’intervention
d’une variable d’endommagement pour représenter la dégradation des résistances inter –
fibre provoquée par la rupture anticipée d’une seule fibre. Cette rupture, de caractère
statistique (par exemple selon une distribution de Weibull), provoque l’apparition de
décohésions aux interfaces fibre – matrice et de microfissures dans le voisinage. Ce
mécanisme peut être modélisé en utilisant une variable d’endommagement dont
l’évolution est basée sur des modèles classiques d’endommagement. Elle est calculée
comme indiqué ci@dessous.
Équation 21
définie par le coefficient . Une distinction entre traction et compression est faite aussi
Équation 22
Équation 23
que la résistance en traction est connectée avec une autre variable d’endommagement ( )
qui sera décrite ultérieurement. Cette variable représente l’endommagement inter fibre et sa
présence est due à l’influence que cet endommagement a sur les ruptures des fibres.
58
(( , 1 2 - !
On vient de décrire le critère mésoscopique utilisé pour déterminer la rupture d’un pli
unidirectionnel. A l’échelle macroscopique, le comportement thermo – viscoélastique est
donné par :
Équation 24
déformations d’origine visqueuse. Elles sont calculées via un modèle spectral développé au
sein de l’ONERA [86,87].
Équation 25
Équation 26
Équation 27
59
(( , 1 2 - !
Équation 28
Équation 29
Équation 30
L’effet de fermeture des fissures en fonction du signe des contraintes est introduit par
le biais des indices d’activation . Les constantes représentent l’effet de
Pour finir cette description du modèle ONERA pour des composites carbone/époxy,
il ne nous reste qu’à décrire le critère de rupture macroscopique (rupture du stratifié). En
effet, ces critères changent en fonction de l’application du stratifié. Le choix fait par les
modélisateurs est l’introduction d’un indicateur normalisé et directionnel défini par
l’équation suivante:
Équation 31
60
(( , 1 2 - !
plus fidèle et complet mais aussi plus complexe, eu égard au nombre de paramètres et à leur
identification. On peut citer comme exemple d’application de ce modèle les travaux décrits
dans les articles [89,90] où le modèle est utilisé pour prédire la rupture de joints composites.
Dans un premier temps, le joint a été modélisé par une éprouvette trouée et boulonnée
(Figure 30).
a) [89]
b) [89]
Figure 30. Quelques essais réalisés et modélisation EF.
Plusieurs essais ont été menés pour tester les capacités du modèle. Dans ces
simulations, un modèle de zone cohésive est introduit pour prendre en compte
l’endommagement des interfaces inter – pli. Le modèle donne de bons résultats en termes
de prédiction de rupture et de cartographies d’endommagement comme on peut le voir sur
les images suivantes. Ces bons résultats dans un cas de calcul de structures complexe
expliquent le succès de ce modèle et son utilisation pour dimensionner les structures
stratifiées.
a) @ [89] b) @ [89]
61
(( , 1 2 - !
Cette approche a été une des premières tentatives pour étudier l’effet de la rupture
locale dans le comportement d’une structure. Des exemples peuvent être trouvés dans des
articles comme [67,91] dans lesquels la méthode de rupture progressive (dite « Progressive
Failure Analysis ») est appliquée au calcul d’éléments finis.
L’avantage de cette méthode est la simplicité de son application. Elle ne définit pas de
forces thermodynamiques ni d’endommagement (au sens de l’évolution progressive d’une
variable indiquant l’état de dégradation du matériau). On décrira dans la suite la démarche
standard de ce type d’approche en donnant quelques exemples représentatifs. La méthode
peut être décrite synthétiquement de la façon suivante : on utilise un critère (par exemple de
type Hashin) pour déterminer en chaque point de Gauss le type de rupture et si l’élément
doit être rompu ou non. Si le critère est vérifié, les propriétés de l’élément concerné sont
dégradées. Cette dégradation peut être imposée différemment selon les auteurs. Par
exemple, dans [91], on considère trois types d’endommagement différents pour lesquels
trois types d’abattement de propriétés sont envisagés. Si la matrice de rigidité élastique du
pli a la forme :
σ 1 Q11 Q12 0 0 0 ε1
σ Q Q 0 0 0 ε 2
2 12 22
. τ 12 = 0 0 Q66 0 0 γ 12
τ 23 0 0 0 Q44 0 γ 23
τ 13 0 0 0 0 Q55 γ 13
Équation 32
les rigidités restantes après les différents modes de rupture (FF= rupture des fibres,
IFF=rupture interface fibre – matrice et Del= délaminage) sont :
62
(( , 1 2 - !
0 0 0 0 0
0 Q 0 0 0
22
Q FF = 0 0 0 0 0
0 0 0 Q44 0
0 0 0 0 0
Équation 33
Q11 0 0 0 0
0 0 0 0 0
Q IFF = 0 0 0 0 0
0 0 0 0 0
0 0 0 0 Q55
Équation 34
3) Délaminage
Q11 Q12 0 0 0
Q Q 0 0 0
12 22
Q DEL = 0 0 Q66 0 0
0 0 0 0 0
0 0 0 0 0
Équation 35
Une fois que les propriétés d’un élément sont abattues, cet élément restera toujours
avec ses nouvelles propriétés. Une fois que tous les éléments ont été cassés au cours d’un
incrément de chargement, la réponse de la structure doit être calculée à nouveau puisque la
rupture de certains éléments change la rigidité globale de la structure et donc sa réponse. Ce
calcul est itératif car il peut y avoir de nouveaux éléments qui vont rompre à cause du
nouvel état provoqué par la rupture d’autres éléments. On passe à l’incrément suivant dès
la convergence du calcul. Graphiquement, on peut résumer ce processus par la Figure 32
63
(( , 1 2 - !
a) [92] b) [92]
Figure 33. a) Maillages utilisés et b) courbes force – déplacement.
On peut voir sur ces figures l’influence du maillage dans la réponse de la structure,
phénomène très courant dans les simulations avec un comportement adoucissant (ce point
64
(( , 1 2 - !
sera traité dans le chapitre suivant). Les résultats obtenus sont cependant en très bon
accord avec les données expérimentales. Des simulations de panneaux avec raidisseurs ont
été aussi réalisées avec des résultats assez satisfaisants.
Les points forts de ces modèles sont principalement leur simplicité et leur facilité
d’implémentation dans des codes numériques. Ils ne nécessitent qu’un critère pour
déterminer le moment d’apparition de la rupture locale et choisir la forme de la matrice de
raideur dégradée. Néanmoins, ces modèles sont moins adaptés pour décrire finement les
endommagements de nature progressive, comme ceux induits par des chargements de
cisaillement, et qui consistent en une combinaison de phénomènes comme la décohésion
fibre/matrice en cisaillement, la microfissuration matricielle, le micro délaminage…
65
(( , 1 2 - !
Figure 34. Changement d’échelle micro – macro en ayant V.E.R et V.E.S. [96]
Ce modèle a été confronté à des essais sur une éprouvette entaillée avec de bons
résultats en ce qui concerne la prédiction des clusters de fibres rompues au cours de l’essai.
La description complète de cette validation se trouve dans [95]. Un des problèmes de ce
modèle est le coût de calcul associé qui rend assez difficile son utilisation à une échelle
industrielle.
2.3 Conclusions
66
(( , 1 2 - !
élevé qui oblige à les utiliser pour le calcul de petites pièces (ou de zones critiques de taille
restreinte), loin de la taille nécessaire pour le calcul de structure. On franchit une échelle
jusqu’à arriver à la taille du pli, l’échelle mésoscopique. A cette échelle, on trouve deux
approches différentes : le modèle mésoscopique de Cachan, placé dans le cadre
thermodynamique de la mécanique de l’endommagement et qui est un modèle de référence
à cette échelle, et le modèle ONERA qui, partant d’un critère construit sur la base de celui
de Hashin, est couplé avec un modèle d’endommagement qui fait évoluer la dégradation du
matériau. Le type de comportement non linéaire modélisé est différent pour chaque
approche : si le caractère endommageable est présent dans les deux cas, on note
l’adjonction du caractère viscoélastique de la matrice dans le cas du modèle ONERA. On
aura constaté que le premier peut simuler la déformation résiduelle mais pas les boucles de
décharge ou l’effet de la vitesse (absence de phénomènes visqueux) et le deuxième ne peut
pas simuler correctement les décharges [98]. Cependant, les deux modèles ont montré leur
performance dans le domaine de la simulation de structures stratifiés minces. En
augmentant davantage la simplicité, on arrive au modèle de rupture progressive qui dégrade
brutalement les propriétés du matériau dès qu’un certain critère est vérifié. Ce type de
modèle est plus simple mais moins fin dans la description des modes d’endommagement. Il
montre aussi des problèmes de convergence dus notamment à la diminution brusque des
propriétés. Une nouvelle génération de modèles d’endommagement se développe
actuellement en incorporant la nature stochastique de la rupture des fibres. Ces modèles
sont capables de bien reproduire la variabilité de la rupture des fibres et de cette façon de
mieux étudier l’effet de cette variabilité dans les simulations de structures.
Par contre, tous ces modèles (sauf le modèle probabiliste dans [94–96]) ont été
développés et utilisés pour simuler des composites stratifiés de faible épaisseur. On peut
voir que le cas des composites bobinés présente des particularités par rapport aux stratifiés.
En particulier, ces structures bobinées sont caractérisées par de fortes épaisseurs, ce qui
nécessite une fine description des phénomènes hors – plan. Afin de traiter le cas de ce type
de matériau peu abordé dans la littérature, nous proposerons une approche différente et
originale qui essaiera d’associer les divers aspects abordés jusqu’ici, à savoir la description
de différents types d’endommagement, la déformation résiduelle, la viscosité de la matrice,
le caractère probabiliste de la rupture des fibres et les fortes épaisseurs. On gardera à
l’esprit la nécessité de proposer une approche simple et de pouvoir l’implémenter dans un
code de calcul industriel.
67
CHAPITRE III:
MODELED’ENDOMMAGEMENT A
DIRECTIONS FIXES
((( , # - 4 5!
Les matériaux composites bobinés ne diffèrent pas des composites stratifiés classiques
en ce qui concerne les types d’endommagements qu’ils subissent. Ces endommagements
sont : la rupture des fibres, la rupture matricielle (macro – fissures), l’endommagement
diffus qui regroupe la décohésion fibre/matrice, la micro fissuration matricielle et le micro
délaminage et finalement, le macro délaminage. Il existe dans la littérature plusieurs études
sur les composites bobinés comme [50,99,100]. Dans ces études, on a constaté que dans la
direction des fibres ou la direction transverse, le comportement reste similaire à celui du
composite stratifié, c’est@à@dire un comportement linéaire avec une rupture très fragile.
Néanmoins, le comportement en cisaillement pur (drapages à ±45º) montre une différence
très claire. Les niveaux de déformation atteints sont plus élevés pour le composite bobiné
que pour les stratifiés. Malgré ces informations, il n’y a pas eu de comparaison de
comportement entre un composite stratifié et un composite bobiné, qu’on distingue par le
procédé de fabrication et donc l’architecture du matériau.
Dans le cadre du projet OSIRHYS IV, une campagne d’essais d’identification a été menée
d’une façon conjointe par Mines ParisTech et l’Institut Pprime. Le matériau bobiné utilisé
est composé de fibres T700SC et de matrice époxyde UF3369. Les drapages testés sont
les suivants : [08], [908], [02, 902]s et [±454]s. Les trois premiers drapages sont dédiés à l’étude
de la rupture des fibres et de la fissuration matricielle. Ces éprouvettes n’ont aucune
différence d’architecture par rapport à un composite stratifié courant. Le comportement en
traction simple des deux premiers drapages est montré sur la figure suivante :
70
((( , # - 4 5!
On peut voir sur la Figure 35 que le matériau suit un comportement attendu, très
linéaire avec une rupture brusque donc fragile. L’essai de traction d’une éprouvette à
drapage [908] nous donne l’information du seuil de rupture matricielle mais pas l’évolution
de la fissuration dans le cas où les plis à 90º sont placés entre des plis orientés dans la
direction de traction. Cette information est donnée par le troisième drapage ([02, 902]s). On
peut voir l’évolution de la densité des fissures lors d’un essai de traction quasi – statique
dans la Figure 36.
Fissuration matricielle
0,8
0,7
Test 1 Test 2
Densité fissures
0,6
0,5
0,4
0,3
0,2
0,1
0
0 20 40 60 80 100
% charge rupture
71
((( , # - 4 5!
l’architecture du matériau bobiné. Dans le cas d’un stratifié, cette notation correspond à
[+45, @45, +45, @45, +45, @45, +45, @45]s. Par contre, dans le cas du bobiné, on doit parler
de pseudo @ couches puisque, dans le même plan, on trouvera des orientations à +45° ou @
45° en fonction de la zone considérée, en raison de l’entrecroisement des mèches. Dans le
Chapitre 1, on a montré le motif créé par le procédé d’enroulement filamentaire. Les essais
présentés dans ce chapitre (et d’une manière générale dans tout ce manuscrit) ont été
réalisés sur des éprouvettes dont l’architecture est identique à celle des réservoirs (le
Chapitre 4 présentera la procédé d’élaboration) et prélevées de telle sorte qu’on trouve la
ligne d’entrecroisement sur l’axe central de l’éprouvette. On trouvera ainsi une alternance
des orientations des deux côtés de l’éprouvette comme le montre la Figure 37.
Figure 37. Schéma d'une éprouvette croisée. Chaque couleur représente une orientation
différente (positive +θ ou négative @θ)
La figure suivante montre la courbe contrainte – déformation axiales pour les deux
types de composites (stratifié vs. bobiné).
Bobiné 100
160
140 80
Stratifié
120
100 60
80 Bobiné
40
60
40 20 Stratifié
20
0 0
0 5 10 0 0,5 1 1,5
Déformation axiale (%) Déformation axiale (%)
72
((( , # - 4 5!
On peut voir ci – dessus que les fissures se sont arrêtées à proximité des
entrecroisements. Cet effet nous montre l’importance d’un paramètre supplémentaire dans
le design des structures bobinées : l’entrecroisement. Cette comparaison illustre la nécessité
de proposer un modèle capable de décrire une cinétique d’endommagement conforme à
ces essais et différente de celles des stratifiés. Pour résumer les types d’endommagement
présents dans ce composite, l’essai sur éprouvettes [±454]s nous montre trois phénomènes
qu’on doit représenter :
73
((( , # - 4 5!
Pour l’instant, on n’a analysé que les phénomènes dans le plan des plis mais les fortes
épaisseurs des structures à modéliser imposent de prendre en compte des phénomènes
hors@plan. La campagne d’essais réalisée dans le cadre du projet OSIRHYS IV n’a pas
explicitement abordé ces effets. En l’absence de données expérimentales hors − plan sur le
composite bobiné (par exemple des essais de flexion quatre point sur des cornières [101]),
nous ferons l’hypothèse de modéliser ce comportement de la même manière que le
comportement transverse [102]. Une caractérisation précise des propriétés hors – plan reste
une perspective de ce travail.
74
((( , # - 4 5!
(déformation) et d’un jeu de variables internes (scalaires ou tensorielles) dont le but est
75
((( , # - 4 5!
de décrire les phénomènes dissipatifs en jeu et agissant ainsi comme témoin de l’histoire du
solide. est donc de la forme :
Équation 36
Le choix a été fait ici d’un potentiel fonction de la déformation . D’autres théories
Équation 37
Équation 38
Équation 39
76
((( , # - 4 5!
Équation 40
Parallèlement aux lois d’état (Équation 38) et (Équation 39), il convient de préciser les
lois complémentaires qui décrivent l’évolution des variables internes , i.e. les processus
valeur scalaire, continue, convexe, non négative et nulle à l’origine dans l’espace de ses
variables. Dans l’hypothèse d’un système à dissipation normale, les lois complémentaires
deviennent :
Équation 41
Il est cependant souvent plus intéressant d’avoir des lois donnant l’évolution des en
fonction des variables duales . On définit alors, par une transformée de Legendre@
peut montrer que possède les mêmes propriétés que . La loi d’évolution des s’écrit
alors :
Équation 42
Équation 43
77
((( , # - 4 5!
Équation 44
Équation 45
Équation 46
avec
Équation 47
d’un potentiel thermodynamique (Équation 36) duquel découlent les lois d’état (Équation
38) et (Équation 39) et d’un pseudo@potentiel donnant l’évolution des variables
définit les lois constitutives du modèle dans le cadre précité des systèmes à dissipation
normale.
• une variable scalaire : cette représentation a été largement utilisée dans la littérature
(voir, par exemple, [105], en lien avec le concept de contrainte effective) en
78
((( , # - 4 5!
Équation 48
matériau lors d’un essai uniaxial. L’avantage d’une telle variable réside essentiellement dans
sa simplicité et la facilité de sa quantification (la comparaison des pentes de décharge d’une
courbe contrainte @ déformation avec la pente initiale fournit directement la valeur de ),
mais sa nature scalaire ne peut pas rendre compte du comportement fortement anisotrope
des matériaux considérés ici.
Leckie [108] apportent la preuve que seuls les tenseurs d’ordre pair en
• une variable tensorielle d’ordre deux : c’est cette dernière solution qui sera retenue
ici ; elle permet de rendre compte du caractère orienté des microfissures tout en
évitant les inconvénients évoqués au point précédent.
Une microfissure est caractérisée, d’une part, par son orientation (normale ) et,
d’autre part, par sa surface S. L’ordre deux impose une expression proportionnelle à ,
de la variable d’endommagement :
79
((( , # - 4 5!
Équation 49
Équation 50
Équation 51
80
((( , # - 4 5!
Cependant, Cormery et Welemane [111] ont montré que les décompositions spectrales
peuvent aboutir à la non unicité de l’énergie libre, et par conséquent à la perte de son
caractère de potentiel thermodynamique. Afin de pallier ces difficultés, [103] ont proposé
une alternative à la représentation de l’endommagement : l’approche à directions fixes.
l’étude, est fixe dans l’espace physique du matériau et chaque densité de microfissures
associée , est une variable interne du modèle amenée à évoluer. La contribution des
caractère fixe des directions d’endommagement. Ce caractère fixe aura une grande
importance dans la construction du modèle adapté aux composites. En effet, la variable
possède trois vecteurs propres qui évoluent avec le chargement ; la loi d’évolution doit
donc gérer simultanément la croissance des défauts et la rotation éventuelle des directions
propres en fonction du chargement. Dans l’approche alternative développée dans [103], les
directions sont fixes et seules les densités associées (scalaires) évoluent. Ces directions ne
81
((( , # - 4 5!
sont pas définies au travers d’une décomposition spectrale, ce qui en évite les inconvénients
précédemment mentionnés. Les tenseurs (et leur nombre) seront choisis de sorte à
Équation 52
est alors un invariant scalaire dans le groupe complet *, i.e. une fonction isotrope. On
supposera de plus que est un invariant polynomial. Il s’agit donc de déterminer une
tout invariant scalaire polynomial de ces mêmes arguments puisse s’exprimer par un
polynôme des invariants générateurs. Ici, et c'est un cas particulier, cet ensemble
d'invariants forme une base, qui est appelée base d'intégrité. On peut trouver une étude très
détaillée de la détermination de ces bases d’intégrité dans [112] ou [113]. De nombreux
travaux utilisent ce concept de représentation des fonctions tensorielles, notamment pour
formuler les lois d’état des solides anisotropes [114], [115].
• invariants de :
• invariants de :
• invariants mixtes :
82
((( , # - 4 5!
Cette base d’invariants est construite à partir de deux tenseurs uniquement. Dans le
cadre de cette thèse, le tenseur deviendra l’ensemble des tenseurs « directionnels »
Équation 53
Il est à noter qu’en choisissant des tenseurs directionnels unitaires (ce qui sera le cas
dans la suite), c’est@à@dire tels que , le nombre d’invariants de base s’en trouve
Dans cette section, on va s’appuyer sur les outils décrits précédemment pour
construire un modèle d’endommagement dédié aux composites bobinés. Le modèle décrit
dans la suite est placé dans la catégorie des modèles mésoscopiques, c’est@à@dire que notre
échelle de référence sera celle du pli homogénéisé. Cette échelle demande un effort du
83
((( , # - 4 5!
point de vue de la modélisation mais nous permettra plus tard de faire des simulations de
structures complexes.
L’utilisation de l’approche à directions fixes devient très avantageuse dans le cas des
matériaux composites car elle privilégie certaines directions d’endommagements d’une
façon très claire. Les modes d’endommagement décrits précédemment dictent le choix des
tenseurs les plus pertinents. Comme dans [116], ces tenseurs directionnels sont de la
dans le but d’être le plus précis possible dans la description de l’état d’endommagement et
pour prendre en compte les effets hors plan justifiés par la forte épaisseur de composite.
Puis, dans une version simplifiée, on réduira la dimension du modèle pour pouvoir l’utiliser
avec des éléments coques (hypothèse de contrainte planes) Dans la suite de ce manuscrit,
on présentera des comparaisons de ces deux versions. À la différence des matériaux pour
lesquels [116] avaient proposé leur modèle, les composites sont fortement anisotropes.
Cette particularité évite de traiter la difficulté de la représentation de l’isotropie initiale
[117]. Dans le cas des matériaux composites, on associe un tenseur directionnel à chaque
direction d’anisotropie, à laquelle correspondent un ou plusieurs modes
d’endommagement. Les tenseurs choisis sont :
ii. . Tenseur lié à la direction transverse (dans le plan du pli) qui représente la
iii. . Tenseur lié à la direction hors – plan par rapport au pli. On associe à ce
84
((( , # - 4 5!
tenseurs sont liés aux directions ±45º par rapport à deux axes donnés. Par exemple,
pour 7 et 6 7 8 , on obtient les tenseurs des directions +45º ( comme positif)
ou @45° ( comme négatif). Ces systèmes sont liés à l’endommagement associé aux
états de cisaillement (décohésion fibre/matrice progressive, microfissuration, micro
délaminage) et de délaminage provoqué par le cisaillement dans le plan du pli.
NHP
NHP NHP
N0-90n
N90-HPp N90-HPn N90-HPp N90-HPn
Ces tenseurs directionnels restent fixes et liés à l’orientation du pli dont la direction de
référence est l’orientation des fibres. Par exemple, lorsqu’on exerce une traction dans le
sens des fibres, on note une rupture des fibres dans le plan qui leur est perpendiculaire
(c’est@à@dire correspondant au tenseur ) On peut suivre un raisonnement identique pour
la rupture matricielle. Pour le cas de l’endommagement diffus (dont les tenseurs sont
orientés à ±45º par rapport aux fibres), les tenseurs associés sont orientés dans les
directions de cisaillement maximal.
85
((( , # - 4 5!
Cet ensemble de systèmes peut être réduit pour utiliser le modèle sous l’hypothèse de
contraintes planes. Le seul changement est la disparition du tenseur et les directions à
±45º associés. Cette version du modèle pourra être implémentée dans des calculs de coques
où l’hypothèse de contraintes planes est admise.
avec i = 0, 90, HP, pour représenter la rupture des fibres, la rupture transverse, la rupture
hors plan et l’endommagement diffus respectivement. Ce dernier endommagement englobe
l’ensemble des phénomènes qui apparaissent d’une façon - lors d’un état de
déformation de cisaillement. Ces phénomènes sont la décohésion fibre/matrice, la
microfissuration matricielle, le micro – délaminage et le macro – délaminage dans le plan du
pli. En ce qui concerne les autres phénomènes dissipatifs, on utilisera les notations pour
86
((( , # - 4 5!
Équation 54
Équation 55
i.
ii.
iii.
l’ensemble des variables et paramètres directionnels intervenant dans ces expressions. Ces
ensembles seront précisés dans la suite ainsi que chacun des termes du potentiel
thermodynamique pour montrer clairement leur construction et les choix faits pour
représenter les effets des différents phénomènes.
suivante
87
((( , # - 4 5!
Équation 56
Équation 57
Équation 58
88
((( , # - 4 5!
serait d’inclure dans tous les termes issus de la Théorie de Représentation des
On distingue trois termes différents qui ont des effets différents sur la rigidité du
cisaillement. Les constantes et agissent dans la sous matrice , respectivement sur les
termes hors diagonale et sur ceux de la diagonale. Les constantes restantes agissent sur
89
((( , # - 4 5!
β α α 0 0 0
α β α 0 0 0
α α β 0 0 0
C=
0 0 0 ζ 0 0
0 0 0 0 ζ 0
0 0 0 0 0 ζ
liant à la direction normale aux fissures. Pour placer les constantes hors diagonale,
obtient des termes qui représentent l’effet des déformations sur la contrainte . Enfin,
endommagements fragiles liés aux directions 0, 90 et HP. Pour illustrer l’effet induit par
l’endommagement fragile sur la rigidité du matériau par la forme de l’Équation 58, on
montre celui obtenu dans le cas de l’endommagement des fibres (lié à la direction 0). Le
potentiel particularisé pour la direction 0 a pour forme :
90
((( , # - 4 5!
Équation 59
Équation 60
Équation 61
Avec ce potentiel (Équation 58), on décrit l’influence de trois des quatre types
d’endommagement présents dans le composite bobiné. Il nous reste donc la description de
la partie du potentiel thermodynamique associée à l’endommagement diffus (progressif) qui
a une forme légèrement différente, comme le montre l’équation ci@dessous :
Équation 62
91
((( , # - 4 5!
Équation 63
totalement dégradé). Pour représenter ce fait, les constantes intervenant dans le potentiel
devront annuler le terme de rigidité correspondant pour traduire le fait qu’une valeur de
représente la rupture du matériau.
Équation 64
On peut voir que ces termes du potentiel thermodynamique ne sont pas construits
avec des fonctions quadratiques mais linéaires. De cette façon, on n’aura pas d’influence
sur la matrice de rigidité mais on obtiendra des d’origine visqueuse. La présence
de deux termes dans l’Équation 64indique que l’effet de la viscosité est principalement
observé sur les essais où le cisaillement joue un rôle important. C’est pourquoi on trouve
92
((( , # - 4 5!
Équation 65
Dans le cadre de la Thermodynamique des Processus Irréversibles, les lois d’état sont
obtenues par dérivation du potentiel thermodynamique. Ces lois d’état seront la contrainte
élastique (obtenue par dérivation par rapport à la déformation élastique), les forces
thermodynamiques associées à l’endommagement et la déformation résiduelle. Les
expressions obtenues sont détaillées dans la suite :
93
((( , # - 4 5!
Équation 66
Équation 67
94
((( , # - 4 5!
On voit dans ces expressions que toutes les déformations auront une influence sur la
valeur de la force thermodynamique. Les déformations associées aux termes diagonaux
ont cependant un rôle prépondérant par rapport aux autres. En ce qui concerne le
glissement, on suit la même démarche :
Équation 68
Équation 69
en couplant les deux phénomènes (glissement et endommagement). Une fois que les
forces thermodynamiques ont été décrites, on peut passer à l’étude des domaines de
réversibilités et des lois d’évolution.
Les domaines de réversibilité des variables internes sont construits à l’aide des forces
thermodynamiques et des fonctions dites « seuils ». L’expression générale des domaines de
réversibilité peut se mettre sous la forme :
Équation 70
95
((( , # - 4 5!
La forme de ces fonctions est dictée par le type d’évolution souhaité pour chaque
variable : au vu des résultats de l’étude du comportement élémentaire du composite bobiné,
on aura besoin de trois types d’évolution différents (voir Figure 44) :
Contrainte (MPa)
Contrainte (MPa)
700 600
1000
600
500
800 500
400
600 400
300
300
400
200 200
200 100
100
0 0 0
0 0,5 1 1,5 0 0,5 1 1,5 0 0,5 1 1,5
Déformation (%) Déformation (%) Déformation (%)
On commence par étudier le premier des trois types d’évolution, celui associé au
comportement fragile. La fonction seuil choisie pour représenter ce comportement a
l’expression suivante :
Équation 71
Équation 72
96
((( , # - 4 5!
Équation 73
Endommagement fragile
1,6
1,4
1,2
0,8
ρ
0,6
0,4
0,2
0
0 0,005 0,01
Déformation
L’endommagement évolue dès que atteint une valeur égal à . Quand la valeur
Si alors .
Si
Équation 74
97
((( , # - 4 5!
sont utilisées pour calibrer la cinétique d’évolution. De cette façon, on pourra adapter le
comportement aux constats expérimentaux. On montre sur la figure suivante les évolutions
de pour trois jeux différents de constantes et leur effet sur le comportement du matériau.
1,2 600
FI=1, FF=2
Contrainte (MPa)
FI=0.5, FF=3
0,8 400
FI=1, FF=2
0,6 300
ρ
FI=1.5, FF=2.5
FI=0.5, FF=3
0,4 200
0,2 100
0 0
0 1 2 3 4 0 0,002 0,004 0,006 0,008
F^0,5 Déformation
a) b)
Figure 46. a) Evolution de la variable d’endommagement pour trois combinaisons différentes des
constantes. b) Effets de ces évolutions sur le comportement du matériau.
Pour identifier ces paramètres, on utilisera les données issues des essais sur éprouvettes
élémentaires comme on le verra dans la partie « Identification ». On vient donc de
présenter la loi d’évolution utilisée dans le cas des endommagements entrainant une rupture
de type fragile. Dans la suite, on analysera la loi choisie pour l’endommagement diffus, qui
a une évolution beaucoup plus progressive. La fonction seuil correspondant à cet
endommagement a une forme plus complexe, comme on peut le voir ci – dessous :
98
((( , # - 4 5!
Équation 75
Équation 76
Équation 77
l’endommagement.
ii. , et . Constantes cinétiques. Elles ont la fonction de contrôler l’évolution
(vitesse et forme).
iii. . Cette constante nous permet de définir un état saturé différent de l’unité.
99
((( , # - 4 5!
0,6 0,6
0,5 0,5
0,4 0,4
ρ
ρ
0,3 C1v = 0.1 0,3 C2v = 0.01
C1v = 0.2 C2v = 0.05
0,2 0,2 C2v = 0.1
C1v = 0.4
0,1 0,1
0 0
0 0,005 0,01 0,015 0,02 0,025 0 0,005 0,01 0,015 0,02 0,025
Déformation Déformation
a) b)
Influence n Influence SAT
0,7 0,6
0,6 0,5
0,5
0,4
0,4
0,3
ρ
ρ
0,3 n = 0.5 SAT = 0.3
n=1 0,2 SAT = 0.4
0,2 n=2
SAT = 0.5
0,1 0,1
0 0
0 0,005 0,01 0,015 0,02 0,025 0 0,005 0,01 0,015 0,02 0,025
Déformation Déformation
c) d)
Figure 47. Influence des paramètres sur le comportement de l'endommagement diffus.
On peut voir que la loi d’évolution choisie est composée de deux parties : la partie
principale correspondant à une exponentielle et une deuxième partie linéaire. La première
partie de la loi (exponentielle) a été déjà utilisée dans d’autres modèles d’endommagement
comme le modèle ONERA [84] dédié aux composites à matrice céramique. Sur la Figure
47, on constate qu’on peut modifier à volonté la forme des courbes pour l’adapter au
comportement expérimental. On voit que affecte la cinétique de la partie exponentielle
(plus haute est sa valeur, plus lente est l’évolution de la variable). D’autre part, influence
la partie linéaire, comme on peut le voir sur la Figure 47b. La pente de la partie saturée
augmente quand on augmente la valeur de ce paramètre. Pour agir sur la forme de la
courbe, on utilise l’exposant (Figure 47c). Enfin, le paramètre change l’état saturé de
la partie exponentielle. L’identification de ces paramètres est faite par recalage de la courbe
de perte de rigidité. Cette courbe est obtenue à partir des pentes mesurées lors d’un essai de
CRP (chargement répété progressif) sur des éprouvettes élémentaires [±454]s. Il est clair
100
((( , # - 4 5!
que le jeu des paramètres ne sera pas unique et l’utilisation de méthodes d’optimisation
pour son identification pourrait améliorer ce processus.
Il nous reste à décrire l’évolution du glissement. La démarche suivie n’est pas identique
à celle présentée pour les variables d’endommagement. Dans ce cas, on utilisera un
formalisme non associé pour obtenir les évolutions des deux variables γ9 et γ:.
Contrairement aux évolutions complexes choisies pour les variables d’endommagement, les
constats expérimentaux [99] sur des composites bobinés nous orientent vers une évolution
plutôt linéaire de la variable de glissement en fonction de la déformation. Pour construire la
fonction seuil, on va décomposer la force thermodynamique correspondante en deux
parties (toujours dans les axes du matériau) :
Équation 78
Équation 79
l’expression est distincte de celle de la fonction seuil à partir duquel nous pourrons
Équation 80
Équation 81
101
((( , # - 4 5!
Équation 82
On notera que l’utilisation d’un formalisme non associé et d’une partition de la force
thermodynamique dans l’expression de la fonction seuil n’est pas une singularité de ce
modèle et peut être trouvée dans d’autres travaux [103].
d’une loi introduite explicitement. On peut choisir d’exprimer cette loi sous la forme d’une
équation différentielle comme on le montre ci@dessous :
Équation 83
terme de l’expression (qui a aussi la forme d’une contrainte). Le terme a une forme
quelconque (a priori), c’est@à@dire qu’il pourra être linéaire ou non linéaire en fonction des
besoins. Il est possible aussi de faire des modèles de viscosité encore plus complexes en
ajoutant plusieurs temps de relaxation (séries de Prony). Cependant, ce type de loi génère
un comportement 2 alors que notre modèle doit être élasto – plastique – visco –
endommageable et coupler la viscosité avec tous les autres phénomènes d’endommagement
102
((( , # - 4 5!
modélisés préalablement. Pour adapter cette loi à notre modélisation, on fera les choix
suivants :
existera un point d’inflexion qui marquera la transition d’un régime pour lequel
(qui est un régime stable) à un régime caractérisé par
Équation 84
103
((( , # - 4 5!
Modèle 1D
35
30 A B
Contrainte (Mpa)
25
20
15
10
5
C
0
0 0,005 0,01 0,015 0,02 0,025
Déformation
Sur cette figure, on peut voir clairement que tous les phénomènes qu’on essaie de
représenter sont présents. D’abord, on peut distinguer deux parties (A et B) qui
correspondent aux zones de comportement viscoélastique (pas d‘endommagement ni de
glissement) et viscoplastique endommageable respectivement. On voit que la viscosité joue
un rôle très important dans le début de la courbe où la déformation résiduelle qu’on
observe est due à la viscosité. Dans la partie B, on peut observer la dégradation du module
de décharge ainsi que l’augmentation de la déformation résiduelle. Dans le dernier cycle, on
a fait un maintien à charge nulle pour voir l’action de la viscosité (point C). On voit bien
que la déformation résiduelle évolue vers la valeur correspondant à la déformation
permanente réelle (glissement résiduel). Les évolutions de la variable , de sa dérivée et de
35 0,06
z 0,04
Z' (Mpa/s)
30
0,02
25
0
20
00,02 1 501 1001 1501 2001 2501
15 00,04
10 00,06
5 00,08
0 00,1
1 501 1001 1501 2001 2501 3001 t (s)
t (s)
On peut remarquer, comme on l’a dit précédemment, que la viscosité évolue dès le
début du chargement et que la contrainte visqueuse n’arrive jamais à atteindre la valeur de
la contrainte imposée (soit la valeur maximale soit la valeur zéro). Un autre détail est le
décalage entre le chargement imposé et la contrainte visqueuse. En effet, si nous regardons
104
((( , # - 4 5!
Jusqu’ici, tout le modèle présenté reste dans le cadre classique des modèles
d’endommagement déterministes. Néanmoins, il est bien connu que la rupture des fibres
n’est pas un phénomène déterministe mais stochastique [123] et [70]. Dans cette étude, on
veut modéliser le comportement du composite de la façon la plus réaliste possible et pour
ce faire, il est nécessaire d’intégrer de quelque façon que ce soit ce caractère aléatoire de la
rupture des fibres sans complexifier à outrance le modèle. La méthode envisagée est décrite
dans la suite.
ii. Placement des fibres dans le composite. On utilise le taux volumique des fibres
pour mesurer la répartition des zones plus riches en matrice ou des pores dans une
zone déterminée comme on peut voir sur la Figure 50:
105
((( , # - 4 5!
Figure 50. Micrographie d'une éprouvette en composite bobiné. Détail des zones
riches en matrice (cercles rouges).
Dans cette figure issue d’une éprouvette fabriquée par enroulement filamentaire (le
procédé sera détaillé dans le chapitre suivant), on peut observer que la densité des
fibres n’est pas la même selon l’endroit considéré. Il est évident que dans les zones
plus riches en matrice ou à forte porosité, les fibres seront plus sollicitées que leurs
voisines. Il est alors plus probable que la rupture apparaisse à cet endroit. Comme
on verra plus tard, la mesure du taux volumique local dépend de la taille de la zone
sur laquelle on le mesure (taille du volume élémentaire représentatif).
106
((( , # - 4 5!
Dans cette loi de mélange, on injecte le comportement stochastique des ruptures des
fibres d’une part et la variation de taux volumique des fibres via d’autre part. De cette
manière, on obtient une contrainte à rupture du pli qui nous servira pour modifier la valeur
On peut observer la simplicité de cette approche, qui est plutôt de nature numérique
puisqu’elle sera basée sur l’affectation de propriétés différentes en chaque point du
matériau. Dans la suite, on parlera de deux versions du modèle :
i. Déterministe : tous les points ont le même seuil d’endommagement des fibres.
ii. Probabiliste : chaque point est affecté d’un seuil différent selon la loi probabiliste ici
développée.
Le but de ce modèle est d’être utilisé dans un calcul de structures complexes comme
les réservoirs bobinés. Pour cette raison, il faut intégrer le modèle dans un code de calcul
par éléments finis. Dans le cadre de cette étude, on utilisera le code commercial ABAQUS,
très utilisé dans le domaine industriel. Ce code permet une haute flexibilité pour introduire
de nouveaux comportements (mécanique, thermique…) de matériaux ainsi que des
couplages multiphysiques. Le but est de montrer que ce modèle peut être un outil de
conception et de développement utilisé par les bureaux d’études des entreprises. Le
modèle est intégré dans le code via l’utilisation de la subroutine UMAT (User MATerial).
Le fonctionnement de cette subroutine est assez simple : le tenseur des contraintes est
calculé en chaque point Gauss à partir de plusieurs données d’entrée comme le tenseur des
incréments de déformation, l’incrément de temps, de température, etc… On peut déjà
remarquer l’avantage de la forme dans laquelle le modèle est construit. En fait, la variable
d’entrée du modèle est la déformation totale, ce qui cohérent avec la structure de la
subroutine. Cela nous permet de calculer directement le tenseur de contraintes sans avoir
besoin de résoudre un système d’équations non linéaires en chaque point Gauss. Une
conséquence de ce calcul direct est le gain en temps de calcul.
107
((( , # - 4 5!
La dépendance au maillage est un phénomène qui a été très étudié depuis l’apparition
des premiers modèles avec des lois de comportement adoucissantes. Cette dépendance se
manifeste de deux façons différentes :
108
((( , # - 4 5!
a) [121] b) [125]
Figure 52.a) Dépendance de la réponse en fonction de la taille de maille. b) Dépendance de la
zone en fonction de la forme et taille de maille.
L’équation d’équilibre (σij,j=0) impose que tous les éléments subissent la même
contrainte. Mais à cause des fluctuations numériques, un de ces éléments dépasse le seuil
d’endommagement et commence à se déformer beaucoup plus que les autres, qui subissent
une décharge élastique (puisque tous voient la même contrainte) comme l’illustre la Figure
54.
109
((( , # - 4 5!
Les méthodes les plus couramment utilisées pour traiter ce problème sont les modèles
à effet retard [126] et les modèles non – locaux [127–130], même s’il existe d’autres
approches comme le « crack band zone » mentionné préalablement [122]. La première
méthode consiste à introduire une « viscosité numérique » de telle sorte que l’effet de
l’endommagement prend du temps à se manifester. Cet intervalle de retard permet au
voisinage d’atteindre aussi le seuil et de s’endommager. On peut trouver plusieurs
méthodes de régularisation dans [131]. Un des problèmes de la régularisation visqueuse est
l’apparition d’une dépendance par rapport au temps complètement artificielle (Figure 54),
même si les ordres de grandeur sont inférieurs à ceux de la viscosité mécanique.
[131]
Figure 55. Influence de la vitesse sur la réponse.
110
((( , # - 4 5!
La deuxième famille des méthodes de régularisation est celle des modèles non locaux.
Dans cette formulation, les équations d’évolution ou les lois d’état liées à une ou plusieurs
variables d’état ne sont pas écrites en un point matériel mais font intervenir une moyenne
dans un volume déterminé [132–136]. Cette approche n’ajoute pas de dépendance par
rapport au temps mais elle introduit une ou plusieurs constantes du matériau : les longueurs
caractéristiques. Ces longueurs indiquent la taille du volume sur lequel on calcule la variable
non locale et qui aura une influence sur la forme et l’extension de la zone endommagée.
Équation 85
3.4.2 Convergence
Une autre voie d’amélioration réside dans le calcul de la matrice jacobienne. Les codes
de calcul utilisent par défaut la matrice jacobienne élastique, qui n’est plus adaptée
lorsqu’apparaissent des endommagements très rapides avec une diminution très brusque de
111
((( , # - 4 5!
la rigidité. Dans ce cas – là, le calcul de la matrice jacobienne exacte (ou approximée si cette
dernière n’est pas calculable) améliore la convergence de façon notable [138,139].
L’expression de la matrice jacobienne associée au modèle développé est montrée ci –
dessous :
Équation 86
Le calcul des contraintes à rupture est réalisé avant le premier incrément de temps. A
chaque point de Gauss, on génère un couple de distributions uniformes N(0,1) grâce à la
commande FORTRAN < =0 Chaque distribution contient 5000 individus (nombre
considéré suffisamment élevé pour que les tirages se distribuent sur tout l’intervalle (0,1)).
Une fois créées, ces deux familles de nombres uniformément distribués entre 0 et 1 feront
partie des données d’entrée pour les algorithmes de génération des distributions normale et
de Weibull décrits dans la suite. Ces algorithmes sont obtenus grâce à la Méthode de
Transformation qui nous permet de créer la distribution souhaitée à partir des données
fournies par une distribution uniforme entre 0 et 1. On sait que la probabilité de générer
un nombre entre et , notée , d’une distribution uniforme est donnée par :
Équation 87
112
((( , # - 4 5!
Équation 88
Imaginons que cette variable fait partie des distributions uniformes générées
probabilités :
Équation 89
Donc
Équation 90
+∞ vaut 1, on obtient :
Équation 91
Équation 92
113
((( , # - 4 5!
Équation 93
Équation 94
créée préalablement pour générer deux éléments qui vont suivre une distribution normale
standard. L’algorithme est le suivant. On calcule le terme :
Équation 95
qui est utilisé pour calculer les deux éléments appartenant à la distribution normale :
Équation 96
Équation 97
Équation 98
114
((( , # - 4 5!
Les moyennes et écart types de chaque tirage sont présentés dans le Tableau 3
Tirage 1 2 3 4
Moyenne 99.98 100.0283 99.9993 99.9882
Ecart type 0.9790 1.0597 1.0759 0.8983
Tableau 3. Moyenne et écart type pour chaque tirage.
On peut voir que, pour les quatre tirages, on obtient des valeurs de la moyenne et de
l’écart type proches de ce qu’on a introduit dans la loi en gardant toujours le caractère
stochastique.
Équation 99
115
((( , # - 4 5!
Équation 100
où est le facteur de forme et celui d’échelle. Un exemple des nombres générés avec cet
Figure 57. Données générées avec un facteur de forme de 4049 et d'échelle de 6.02.
116
((( , # - 4 5!
De cette façon, tout le processus est complètement aléatoire et on peut être sûr de ne
pas introduire de biais dans les valeurs concernées.
L’émission acoustique, dont les principes ont été rappelés dans le premier chapitre,
nous donne une estimation qualitative du moment d’apparition et de la position de
l’endommagement. La classification des échos selon le mode d’endommagement qui l’a
créé est beaucoup plus complexe et on ne peut pas être sûr des résultats obtenus. On
profite du fait que le modèle de comportement développé dans ce chapitre intègre plusieurs
modes d’endommagement pour essayer de créer des échos virtuels qu’on enregistrera en
utilisant les sorties de données spécifiques dans ABAQUS. Parmi toutes les variables
possibles, c’est l’énergie absolue qui semble la plus facile à différencier et identifier. On
montrera dans le Chapitre 5 des comparaisons entre la simulation de la libération d’énergie
acoustique et les enregistrements expérimentaux.
respectivement. La raison d’avoir choisi trois familles différentes s’appuie sur l’analyse des
échos lors des essais de traction. On montre ci@dessous une courbe -
correspondant à un essai de traction d’une éprouvette entaillée [±454]s.
117
((( , # - 4 5!
1.5
4
Energie Absolue
Energie Absolue
3
1
2
0.5
1
0 0
200 400 600 800 1000 1200 200 400 600 800 1000 1200
Temps (s) Temps (s)
a) b)
Figure 59. a) Energie absolue vs temps. b) zoom dans la zone de basse et moyenne énergie.
On observe dans cet essai des échos de moyenne et basse énergie (Figure 59b), ce qui
a été attribué dans la littérature à l’endommagement matriciel et diffus [50]. Au vu de cette
courbe, il semble raisonnable d’utiliser deux niveaux de signaux énergétiques qui seront
associés à et . La rupture de fibre aura un signal beaucoup plus énergétique, qu’on
seul coup. Cette option n’est pas légitime lorsqu’on travaille avec des endommagements de
type progressif puisqu’ils libèrent plusieurs échos dus aux décohésions, aux microfissures…
tout au long du processus de déformation. La démarche choisie dans le cadre de cette étude
est un compromis entre ces deux méthodes : on discrétise par exemple en trois intervalles
sur la Figure 60 la libération d’énergie lors de l’évolution de la variable d’endommagement.
118
((( , # - 4 5!
a) b)
Figure 60. Méthode choisie pour représenter des signaux énergétiques.
équivalente à l’énergie de l’événement type choisi pour chaque famille d’échos associés aux
endommagements (Figure 60b). Ces sont calculés comme le produit de l’énergie totale de
l’évènement type par l’intervalle des racines des forces thermodynamique correspondantes.
Équation 101
Où les valeurs montrées dans la Figure 60 sont les valeur moyennes de chaque
intervalle des forces thermodynamique comme nous voyons ci – dessous.
Équation 102
119
((( , # - 4 5!
pour chaque
3.5 Identification
Comme il a été commenté dans la section 3.1, les drapages utilisés aux fins
d’identification sont les suivants :
ii. [908] et [02, 902]s pour identifier la rupture matricielle et le processus de fissuration
dans le cas des plis à 90° placés entre des plis orientés dans la direction du
chargement.
iii. [±454]s pour identifier les conséquences des phénomènes liés au cisaillement
(endommagement diffus, déformation résiduelle et viscosité).
120
((( , # - 4 5!
Des modèles EF ont été construits pour simuler les essais de traction sur ces groupes
d’éprouvettes. On a utilisé des éléments de type coque à quatre nœuds et intégration
complète S4. Dans les futures simulations 3D (voir Chapitre 5), les propriétés hors@plan
seront assimilées à celles identifiées dans la direction transverse sur les éléments 2D . Les
conditions aux limites choisies sont un blocage des déplacements et des rotations d’un côté
sur un rectangle de dimensions 25mm x 50mm qui représente le mors fixe de la machine.
De l’autre côté, on impose un déplacement donné à la même vitesse que les essais
(0.5mm/s). Une particularité de ces modélisations est le fait d’avoir essayé de représenter la
structure bobinée en alternant des séquences de stratifié. Ainsi, on a dessiné sur
l’éprouvette le motif créé par le processus d’enroulement filamentaire en fonction de la
largeur de la mèche de composite utilisé. Deux types différents de séquence sont créés : un
premier noté +θ sur la figure suivante et qui représente l’empilement [+45, @45, +45, @
45…] et un deuxième noté −θ dont l’ordre est inversé [@45, +45, @45, +45…]. L’éprouvette
résultante est montrée sur la Figure 62.
+θ +θ
-θ -θ
121
((( , # - 4 5!
sorte que pour une valeur unité de l’endommagement, le terme affecté de la matrice de
rigidité effective deviendra nul. La seule identification de la matrice élastique reste
suffisante pour déterminer ces constantes d’effet de l’endommagement.
Test
1500 30
1000 20
500 10
0 0
0 0.005 0.01 0.015 0 1 2 3 4 5 6
Déformation axiale Déformation axiale @3
x 10
Figure 63. Identification des comportements à 0º et 90º.
Pour compléter la définition des lois d’évolution, il faut donner une valeur aux paramètres
proche du seuil initial pour imposer le caractère brusque propre à cet endommagement. Il
serait envisageable de préciser plus rigoureusement cette valeur en connaissant l’énergie de
rupture . Dans le cas de l’endommagement matriciel, on dispose, grâce à l’essai de
traction sur la séquence [02, 902]s, de données suffisantes pour déterminer précisément ce
paramètre. On utilise la courbe de la Figure 36 (normalisée par rapport à la valeur
122
((( , # - 4 5!
Fissuration Matricielle
1,2
normalisée
0,8
Sim
0,6
0,4
0,2
0
0 20 40 60 80 100
%Charge à rupture
Figure 64. Identification fissuration matricielle.
Traction [±454]s
180
Test 1
160
Test 2
Contrainte axiale (Mpa)
140
120 Test 3
100
80
60
40
20
0
0,00 2,00 4,00 6,00 8,00 10,00 12,00 14,00
Déformation axiale (%)
Figure 65. Courbes de traction [±45 4 ] s utilisées pour l'identification.
On exploite la courbe CRP (qui contient des décharges), de laquelle on peut extraire
des informations sur la dégradation de la rigidité en cisaillement et sur la déformation
résiduelle. On transforme préalablement la courbe montrée dans la figure ci –dessus en
termes de et . En ce qui concerne la rigidité, on mesure le module de décharge
obtenu à partir des deux points extrêmes de chaque boucle (lignes bleues) et on trace son
123
((( , # - 4 5!
0,6
4
0,5 Simulation
3
0,4
0,3 2
0,2 Experimentale
1
0,1 Simulation
0 0
0 0,02 0,04 0,06 0,08 0,1 0,12 0 1 2 3 4 5 6 7 8
γ12 totale (%) γ12 totale (%)
Figure 66. Identification de l'évolution de l'endommagement diffus et le glissement.
On constate dans cette figure que les deux cinétiques (endommagement et glissement)
sont capables de reproduire correctement le comportement observé dans les essais de
caractérisation. Il ne faut pas oublier que ces deux variables sont couplées et c’est ce
couplage qui permet une corrélation aussi bonne entre les résultats expérimentaux et la
simulation, surtout dans le cas de l’endommagement diffus où il y a un changement de
tendance très brusque autour de 1.5 % de déformation. C’est à cet instant du chargement
que le glissement commence à augmenter fortement.
Pour caler les paramètres liés au modèle visqueux, on a utilisé la courbe pour
essayer de reproduire au mieux les phases de descente où on voit une relaxation provoquée
par la viscosité. On montre la comparaison des courbes expérimentales et simulées dans la
Figure 67. Il est clair qu’avec un seul temps de relaxation et une fonction linéaire de la
déformation on n’arrive pas à reproduire totalement le comportement, surtout dans la
phase finale de l’essai.
124
((( , # - 4 5!
Évolution de la déformation
0,02
Simulation
0,015
Test
γ12
0,01
0,005
0
0 1000 2000 3000 4000 5000
t (s)
Tous les paramètres qui ont été identifiés sont présentés dans le tableau suivant.
8@0 @ 0> @ @ 0? A@
Tableau 4. Jeu de paramètres.
Pour vérifier que le comportement du modèle est satisfaisant, on compare les courbes
contrainte axiale – déformation axiale numérique et expérimentale sur la séquence [±454]s
(Figure 68).
Simulation
150
100
50
0
0 2 4 6 8 10 12 14
Déformation axiale (%)
Figure 68. Comparaison test @ simulation de la traction sur [±45 4 ] s bobiné.
125
((( , # - 4 5!
On peut voir que des boucles apparaissent dans la simulation mais d’une taille plus
petite que les boucles réelles. On constate la perte de rigidité, l’ouverture des boucles et
l’évolution d’une importante déformation résiduelle.
La section 3.4.3 avait décrit la démarche développée pour intégrer dans le modèle la
rupture probabiliste dans la direction des fibres. Deux sources de comportement
stochastique avaient été envisagées : la fibre seule et le taux volumique de fibres. Il faut
donc identifier ces deux sources pour pouvoir décrire correctement le comportement
probabiliste. Un point important de cette démarche d’identification est que les essais de
traction sur la séquence [08] ne sont plus des essais d’identification mais constituent une
première validation de l’approche. Dans la démarche probabiliste, on descend à l’échelle de
la fibre pour y identifier directement son comportement. Les micrographies d’éprouvettes
élaborées par enroulement filamentaire, quant à elles, fournissent des informations sur la
répartition des fibres et le taux volumique. On commence donc en montrant
l’identification du comportement des fibres.
Pour identifier la distribution de la contrainte à rupture des fibres, on s’est servi d’un
essai de traction sur des fibres réalisé par notre partenaire du projet, Mines ParisTech. La
longueur de ces fibres est de 25mm. On montre ci@dessous la probabilité cumulée de
rupture des fibres par rapport à la contrainte appliquée.
0.8
Probabilité
0.6
0.4
0.2
0
2500 3000 3500 4000 4500 5000 5500 6000
Contrainte (Mpa)
126
((( , # - 4 5!
Équation 103
utilisant le logiciel MATLAB, la valeur obtenue pour ces deux paramètres est de =4465.6
Identification Weibull
1
Test
0.8 Weibull identifié
Probabilité cumulée
0.6
0.4
0.2
0
2000 2500 3000 3500 4000 4500 5000 5500 6000
Contrainte (Mpa)
Pour identifier ce paramètre, on s’est basé sur l’étude des micrographies sur lesquelles
on peut voir des fibres perpendiculaires au plan d’observation. Afin de bénéficier d’une
large surface d’observation, on a utilisé des éprouvettes quasi@isotropes [02, ±452, 902]s
photographiées sur la tranche et on a exploité les plis à 90°. La méthode consiste à traiter
numériquement l’image pour la convertir en signal binaire de telle sorte que les fibres
(claires) ont une valeur de 1 et la matrice/porosité (sombre) a une valeur de 0. Cette
technique, qu’on qualifie de « seuillage » (« thresholding » en anglais) a été utilisée dans
d’autres travaux dont le but était l’étude de la distribution des taux volumiques [4]. Le
fonctionnement de cette technique est assez simple :
127
((( , # - 4 5!
biais dans les mesures. Ensuite, ces images sont converties en « niveaux de
gris ». Chaque pixel aura une valeur comprise entre 0 et 255.
ii. Une fois les images prises et converties, elles doivent être retouchées dans
le but de faciliter le choix d’un seuil de niveau de gris correct. Cette étape
est réalisée à l’aide d’un logiciel de traitement d’image type Photoshop® ou
Gimp®. On choisit un seuil de niveau de gris à partir duquel on divisera
l’image en « fibres » et « no@fibres ».
iii. Enfin, on doit choisir une fenêtre d’étude (rectangle) sur laquelle on viendra
mesurer le nombre de pixels qui correspondent aux fibres. Comme on
connait la taille de la zone d’étude (en pixels), le taux volumique est
directement obtenu. Si nous répétons ce processus dans plusieurs endroits
de l’image (minimum 40 – 50 fois), on obtiendra une série de valeurs de
taux volumique, distribuées en suivant une loi statistique à identifier.
Ce sont les trois étapes générales suivies dans la méthode de calcul du taux volumique
par analyse d’image. Pour automatiser le processus, on crée un code MATLAB qui
automatise la dernière partie de la méthode (calcul des taux). Dans notre étude, on a pris
plusieurs photographies différentes qu’on a assemblées en une seule image (Figure 71).
Comme on peut le voir dans la Figure 74, on part des images qui ont une dimension
0.54736 x 0.4012mm. L’assemblage final est composé de 16 photographies. La dimension
128
((( , # - 4 5!
de la zone étudiée est de 2.18 x 1.6152 mm. De cette façon, on dispose d’une aire de
matériau suffisante pour pouvoir étudier l’influence de la taille de la fenêtre d’étude dans la
distribution du taux volumique. En étudiant cette image avec plusieurs tailles de fenêtre
différentes, on obtient la probabilité suivante de trouver une certaine quantité de fibres
dans la fenêtre d’étude.
0.8
Probabilité cumulée
0.6
Vf(0.47mm) data
0.4 Vf(0.235mm) data
Vf(0.117mm) data
Vf(0.0585mm) data
0.2
Vf(0.02925mm) data
Vf(0.585mm) data
0
0.4 0.45 0.5 0.55 0.6 0.65 0.7 0.75 0.8
Vf
a) b)
Figure 72. Taux volumique en fonction de la taille de fenêtre utilisée.
On voit sur la Figure 72b la probabilité cumulée pour une large gamme de tailles de
zones d’étude. Ce comportement est en accord avec les attentes, c’est@à@dire que plus la
taille de la fenêtre est petite, plus on trouvera de dispersion dans la probabilité cumulée. Par
exemple, avec une taille de 0.02925 mm, on obtient une forte influence de la disposition de
clusters de fibres et de la présence de zones riches en matrice ou en porosité. C’est
pourquoi on trouve des valeurs se distribuant entre 0.4 et 0.8. Au contraire, la distribution
correspondant à la taille de fenêtre la plus grande (0.585mm) montre une dispersion
beaucoup plus faible en raison de l’effet homogénéisant de la taille de la fenêtre. Vu la
forme des distributions, il semble raisonnable utiliser une loi normale pour les identifier.
Un exemple de trois de ces distributions identifiées par une loi normale est montré dans la
figure ci –dessous (ID dénote la courbe identifiée).
129
((( , # - 4 5!
Identification Vf
1
Vf(0.235mm) data
Vf(0.235mm) ID
0.8 Vf(0.02925mm) data
Probabilité cumulée
Vf(0.02925mm) ID
Vf(0.585mm) data
0.6 Vf(0.585mm) ID
0.4
0.2
0
0.3 0.4 0.5 0.6 0.7 0.8 0.9
Vf
Cette figure montre que la loi normale est très bien adaptée à ces jeux de données
quelle que soit la taille de fenêtre utilisée pour étudier le matériau. Il suffit de modifier les
paramètres de la loi normale pour retrouver les courbes expérimentales, en fonction de la
taille de la fenêtre d’étude. On montre dans la figure suivante les évolutions de la moyenne
et de l’écart type de la loi normale en fonction de la taille de fenêtre.
0.1
0.585
Ecart type
0.08
Moyenne
0.58
0.06
0.575
0.04
0.57 0.02
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0 0.05 0.1 0.15 0.2 0.25 0.3 0.35
Aire (mm2) Aire (mm2 )
a) b)
Figure 74. Identification des paramètres de la loi normale et dépendance par rapport à la taille.
On voit sur cette figure que la moyenne de la loi normale oscille autour d’une valeur
proche de 0.5825. Par contre, l’écart type manifeste un comportement clairement
asymptotique qui semble tendre vers la valeur de 0.02. Pour s’assurer que cette valeur est
une asymptote, il faudrait étudier des domaines encore plus vastes, ce qui n’est pas
nécessaire dans le cadre de cette étude où la taille des éléments finis qui seront choisis dans
le calcul de structures est incluse dans la gamme de taille des domaines étudiés. En effet, la
distribution du taux volumique de fibres intervenant dans le calcul de la contrainte sera
choisie en fonction de la taille des éléments finis qui seront utilisés pour réaliser les calculs
130
((( , # - 4 5!
de structure. On établit ainsi un lien direct entre la taille de la fenêtre d’étude (Figure 72) et
la discrétisation spatiale d’une structure lors d’un calcul EF.
Les derniers paramètres qui nous restent à identifier sont les valeurs de l’énergie
acoustique associée aux trois modes d’endommagement capables d’émettre des échos
virtuels : la rupture des fibres, la rupture matricielle et l’endommagement diffus. On
utilisera les données acoustiques provenant de deux essais pour identifier les différentes
énergies. On remarquera le caractère qualitatif de la démarche suivie pour l’identification de
cet « outil ». Ces deux essais sont :
ii. Traction sur une éprouvette dont les fibres sont orientées à 0º. On a choisi
les données provenant d’un essai sur une éprouvette entaillée [02, ±202]s (ces
éprouvettes seront décrites en détail dans le prochain chapitre) sur laquelle
on observe principalement des ruptures des fibres.
On montre dans la Figure 75 la courbe d’énergie absolue de l’éprouvette [02, ±202]s sur
laquelle on identifie l’énergie correspondant à la rupture des fibres. Les courbes utilisées
pour l’identification des échos des endommagements diffus et matriciel ont été présentées
dans la section 3.4.4 (Figure 59).
131
((( , # - 4 5!
Figure 75. Energie absolue détectée dans un essai de traction sur une éprouvette entaillée [0 2 ,
±20 2 ]s.
Le rectangle rouge indique la zone d’énergie associée aux ruptures des fibres plus ou
moins isolées. Vers la fin de l’essai, des clusters plus grands de fibres commencent à
rompre en émettant des signaux extrêmement énergétiques. Finalement, les valeurs choisies
sont présentées dans le tableau suivant.
Energie absolue
Endommagement Nº d’émissions Intervalles
associée
=0.7
5.75 x 105
3 =0.9
B 4.4 x 105
=1 1.2039 x 106
# 1 =0.1 40000
=0.2
407
=0.4
4 463
55
=0.5 417
611
=SAT
On validera cette identification dans le chapitre 5 avec une comparaison entre l’énergie
enregistrée au cours des essais et celle obtenues dans la simulation.
3.6 Conclusion
132
((( , # - 4 5!
est construit en utilisant les outils de la théorie de représentation des fonctions tensorielles.
Cette théorie nous permet d’établir des combinaisons des invariants (quadratiques) qui
auront l’effet souhaité sur la rigidité du matériau. En effet, elle permet de relier directement
l’influence d’un système d’endommagement donné aux termes correspondants dans la
matrice de rigidité. L’anisotropie initiale puis induite par l’endommagement est incorporée
d’une façon naturelle par ce type de construction. En suivant la démarche propre à la
Thermodynamique des Processus Irréversibles, on obtient la contrainte par dérivation de
ce potentiel par rapport à la déformation. Si l’on dérive ce potentiel par rapport à chaque
variable d’endommagement ou déformation résiduelle, on obtient les forces
thermodynamiques correspondantes qui pilotent le phénomène associé. Avec ces forces
133
((( , # - 4 5!
Le but de ce modèle est d’être utilisé dans les calculs EF des structures complexes. Il a
été implémenté via une subroutine User MATerial (UMAT) dans le logiciel de calcul
éléments finis ABAQUS. Deux versions différentes ont été implémentées, une version 3D
complète qui nous donnera des informations très détaillées sur le comportement et les
modes d’endommagement et une version coque qui n’active pas l’endommagement hors
plan mais qui nous permet de gagner énormément en temps de calcul.
En ce qui concerne l’identification du modèle, elle reste assez simple car on a besoin
de trois essais (traction sur 0°, ±45°, 90° ou 0°@90º) pour caractériser les propriétés
134
((( , # - 4 5!
Cette démarche appliquée dans cette étude peut être utilisée pour n’importe quel
matériau privilégiant des directions d’endommagement spécifiques. Ainsi, on pourrait
utiliser le modèle pour étudier des stratifiés, des tissus 2D, 2D ½ et 3D, des composites à
matrice céramique, … On validera le modèle dans les chapitres suivants en comparant les
résultats des simulations d’éprouvettes entaillées dans un premier temps puis du réservoir
bobiné après avec les données expérimentale.
135
CHAPITRE IV:
" *##'+ # $
$'(""
(C , ! - ,
4.1 Introduction
La géométrie utilisée a été décidée dans le but d’assurer un état de contrainte complexe
et la présence de forts gradients. De cette façon, on s’attend à observer des
endommagements importants et de natures différentes sur une même éprouvette. La
confrontation des simulations du modèle développé dans le chapitre précédent à la réponse
des éprouvettes en présence d’une telle complexité de comportement constitue une
première étape de validation de l’approche, en prélude à son utilisation sur les réservoirs
bobinés.
Toutes ces éprouvettes ont été fabriquées par la procédure d’enroulement filamentaire
par le CEA Le Ripault. On fabrique d’abord un tube en bobinant un mandrin avec le
drapage [+θ/−θ] souhaité. On coupe le tube le long de l’axe pour pouvoir mettre à plat le
cylindre de matériau composite bobiné et on assemble les différentes orientations. On
effectue alors la cuisson pour obtenir une plaque de matériau composite bobiné avec le
drapage souhaité. Toutes ces éprouvettes ont été prélevées de telle sorte que la ligne
d’entrecroisement des mèches soit placée au milieu de l’éprouvette. On obtient ainsi
l’élément caractéristique de l’architecture bobinée (l’entrecroisement) placé dans la zone
critique (l’entaille) pour pouvoir étudier son influence sur le comportement.
138
(C , ! - ,
a) b)
Figure 78. a) Contour des éprouvettes. b) Types de motifs obtenus selon la largeur de mèche.
Les drapages fabriqués dans le cadre de cette étude et testés par la suite sont :
i. Les familles contenant des fibres orientées à 0º. Ces drapages ont été créés pour
pouvoir étudier des phénomènes de délaminage et l’influence de l’orientation des
plis croisés (+θ/−θ) sur la cinétique de ce type d’endommagement. De plus, le
drapage quasi – isotrope [02,±452,902]s nous permettra de coupler tous les modes
d’endommagement dans un même essai.
139
(C , ! - ,
ii. Le drapage [±454]s sur lequel on verra l’influence d’un fort gradient sur le
développement de l’endommagement diffus ainsi que la déformation résiduelle et
les phénomènes visqueux.
Ces drapages ont été testés en traction monotone à trois températures différentes, @
40ºC , 20ºC et 85ºC. La présentation des résultats de ces essais en température dans la suite
de ce chapitre répond à un objectif double : (i) les essais à température ambiante
constitueront une base de confrontation et de validation du modèle, (ii) même si la suite de
ce mémoire ne présentera pas de comparaison simulation / expérience sur les éprouvettes
entaillées à @40°C et 85°C, ces données expérimentales permettront de mieux appréhender
les spécificités des mécanismes de déformation et d’endommagement caractéristiques des
composites bobinés. Ces essais participent ainsi de l’amélioration de la connaissance du
comportement de ce type de matériau. Tous ces essais ont été multi@instrumentés avec
émission acoustique, corrélation d’image numérique, suivi optique, observations
micrographiques et rayons X post mortem. On commencera le chapitre par la description
du dispositif expérimental utilisé dans les essais de traction. Puis le comportement sera
étudié.
Tous les essais de traction ont été réalisés sur une machine INSTRON 4505
électromécanique équipée de mors auto@serrants. La vitesse de déplacement de la traverse,
identique pour tous les essais, est de 0.5mm/min. Cette machine est aussi équipée d’un
four pour les essais à basse et haute températures. La régulation de cette température est
faite par un contrôle automatique et des thermocouples. Le déplacement macroscopique de
l’éprouvette est assimilé au déplacement mesuré entre deux points placés à 45 mm de part
et d’autre du centre de l’entaille. Pour mesurer ce déplacement, on a utilisé un dispositif de
suivi optique de la position des points. On a également incorporé un dispositif d’émission
140
(C , ! - ,
acoustique (Euro Physical Acoustics) dans le but d’avoir le plus d’informations possibles
sur les processus d’endommagement subis par les éprouvettes. Les capteurs sont du type
Micro 80 et le logiciel utilisé est AE WIN, qui nous permet d’enregistrer et visualiser en
temps réel de nombreuses variables comme l’énergie cumulée, le nombre d’évènements, les
amplitudes, la localisation etc. Pour finir, les images nécessaires à la réalisation de la
corrélation d’images ont été prises avec une caméra HD Redlake EC 11000. Ces images
sont traitées par un logiciel de traitement d’image (GIMP ou Photoshop). Une fois les
images obtenues, on utilise un autre logiciel développé au sein de P’ (GRANU) qui réalise
la corrélation d’image et le calcul du déplacement pour chaque point de la surface. On
impose les déplacements obtenus à un maillage (ABAQUS) de la zone étudiée pour calculer
le champ de déformation. Le montage est schématisé sur les Figure 79, Figure 80 et Figure
81.
141
(C , ! - ,
(90mm) entre les deux points de suivi optique et ∆ l’allongement mesuré entre ces deux
points.
[02 , ±202]s
40000
35000
20ºC a
85ºC
30000
@40ºC
25000
Force (N)
20000
15000 b
10000
5000
0
0,00 0,50 1,00 1,50
Déformation axiale (%)
142
(C , ! - ,
[02 , ±302]s
35000
20ºC
30000
@40ºC
25000
85ºC
Force (N)
20000
15000
10000
5000 a
0
0,00 0,50 1,00 1,50 2,00
Déformation axiale (%)
Figure 83. Courbe de traction à trois temperatures [0 2 , ±30 2 ] s .
10000
5000
a
b
0
0,00 0,50 1,00 1,50 2,00
Déformation axiale (%)
Figure 84. Courbe de traction à trois températures [0 2 , ±45 2 , 90 2 ] s
On voit clairement dans ces trois figures l’effet de la température sur le comportement
mécanique. Dans les trois cas, on a des comportements très similaires pour les essais à @
40ºC et 20ºC et une chute de rigidité non négligeable dans le cas de la haute température.
En ce qui concerne le comportement à rupture, il varie selon le drapage considéré et aussi
selon la température. Les éprouvettes [02, ±202]s montrent un comportement clairement
linéaire jusqu’au point a (Figure 82). A ce moment de l’essai, des délaminages autour de
l’entaille apparaissent et, à @40ºC et 20ºC, la rupture des fibres survient, menant quasi@
instantanément à la ruine totale de l’éprouvette. En revanche, à 85ºC, on observe également
du délaminage, mais l’éprouvette est capable de résister alors que ce délaminage se propage
(points b) jusqu’à atteindre finalement la rupture des fibres.
143
(C , ! - ,
On finit l’analyse par les éprouvettes [02, ±452, 902]s qu’on qualifiera de quasi –
isotropes (Figure 84). L’essai à @40ºC montre un comportement proche de celui de la
famille [02, ±202]s car on observe une courbe très linéaire jusqu’au moment de la rupture
des fibres. L’essai à 20ºC montre par contre une non linéarité due très probablement à un
problème de suivi optique (seule une éprouvette, dans le lot des trois testées
systématiquement pour confirmer la reproductibilité des essais, a pu être exploitée à cette
température). Dans ce cas, on observe un premier délaminage juste avant la rupture finale
de l’éprouvette qui cause une diminution assez brusque de la rigidité comme le montrent
les deux flèches vertes (foncée et claire) qui représentent les raideurs pré@ et post@
délaminage. Finalement, la courbe de l’essai à 85ºC montre le même comportement que les
144
(C , ! - ,
autres drapages avec une forte chute de rigidité. Une particularité de cette séquence est la
grande différence entre la déformation à rupture à chaud et celle à froid/ambiante. En
résumé, on observe que les séquences contenant des plis à 0º montrent un comportement
qualitativement similaire avec une partie initiale linéaire suivie de sauts des déplacements
et/ou une rupture finale provoquée dans tous les cas par la rupture des fibres à 0º.
C’est le moment d’étudier maintenant les deux drapages formés seulement de plis
entre–croisés, à savoir [±454]s et [±67.54]s. On montre ci@dessous les courbes de traction
correspondant au premier drapage :
[±454 ]s
6000
20ºC
5000
a 85ºC
4000 @40ºC
Force (N)
3000
2000
1000
0
0,00 5,00 10,00
Déformation axiale (%)
Figure 85. Courbes de traction [±45 4 ] s À trois températures.
145
(C , ! - ,
a) b)
Figure 86. a) Rupture de l'éprouvette entaillée bobinée. b) rupture d'une éprouvette stratifiée.
146
(C , ! - ,
[±67,54 ]s
1600
1400 20ºC
1200 85ºC
1000 @40ºC
Force (N)
800
600
400
200
0
0,00 0,50 1,00 1,50 2,00 2,50
Déformation axiale (%)
Figure 87. Courbes de traction [±67.5 4 ] s à trois températures.
a) b)
Figure 88. DMTA dans une éprouvette à 0º (a) et 90º (b).
147
(C , ! - ,
élastiques à 90º dans le cas de réponses complètement visqueuses. La tangente de cet angle
représente le rapport entre les modules élastique et visqueux du matériau.
Équation 104
i. Les séquences contenant des plis externes à 0º montrent une légère diminution de
la rigidité, notamment provoquée par la dégradation de la rigidité des plis croisés
(±20°, ±30° et ±45°/90°).
148
(C , ! - ,
149
(C , ! - ,
Groupe I
Splitting Délaminage
90% de la
charge à rupture
Figure 89. Séquence des endommagements dans une éprouvette quasi @ isotrope.
On constate que le splitting apparait d’abord au fond des entailles. Après son
amorçage, les fissures commencent à se propager jusqu’à l’apparition des macro@
délaminages des zones au@dessus et en@dessous des entailles. Finalement, la rupture des
fibres arrive d’une façon brutale. La Figure 90 indique les types de dommage présents
autour des entailles des éprouvettes à [02, ±202]s après essai (dans tout le reste de ce
chapitre, les micro endommagements seront colorisés sur les micrographies lorsqu’il sera
nécessaire de les mettre en évidence).
150
(C , ! - ,
placée dans le centre de l’entaille mais sur la gauche (encadré en pointillé rouge). On trouve
ce phénomène (dissymétrie de la fissuration matricielle) dans tous les plis croisés du
Groupe I. Les simulations de l’endommagement de ce type de séquence (Chapitre 5)
reviendront sur cette spécificité. Il faut noter aussi que le délaminage apparait entre les plis
croisés, mais surtout entre les plis à 0º et les plis croisés adjacents (flèche blanche). Ce
phénomène est plus visible dans les éprouvettes quasi – isotropes comme on le montre ci –
dessous.
On voit bien que les fissures (en vert) des plis croisés ne sont pas placées dans le
centre de l’entaille (ligne jaune) mais des deux côtés (entourés en pointillé bleu). Seuls les
plis à 90º (pli central) subissent une distribution plus homogène des fissures. Si l’on se
concentre sur le délaminage (rouge), on constate que les interfaces des plis à 0º et à 90º
avec les plis adjacents correspondants sont beaucoup plus endommagés que les interfaces
entre plis croisés. On voit aussi le splitting dans le pli à 0º indiqué par la flèche blanche.
Dans le cas des éprouvettes [02, ±302]s, l’état extrêmement endommagé de l’éprouvette
empêche l’interprétation des observations post mortem, comme on voit dans la figure
suivante :
151
(C , ! - ,
Un détail important sur cette figure est la forte rupture par cisaillement qu’on trouve
dans les plis croisés (flèche rouge). On attendait un comportement de cette forme dans les
éprouvettes quasi – isotropes mais ce sont celles contenant des plis à 30º qui montrent une
rupture plus typique du cisaillement (celle qu’on voit dans les éprouvettes sans entaille à
±45º). Les délaminages observés sont aussi beaucoup plus brutaux que dans les autres
drapages du Groupe I (flèche bleue). Même si l’entaille agit comme un concentrateur de
contraintes, on a observé des endommagements placés assez loin de l’entaille et qui se sont
propagés de l’entaille vers les extrémités. Par exemple, dans les éprouvettes [02, ±202]s (qui
montrent moins d’endommagement matriciel ou de délaminage), on trouve des fissures très
loin de l’entaille Il semble raisonnable de penser qu’il s’agit de fissures générées dans le
fond d’entaille et qui se sont propagées jusqu’à arriver au bord. On peut vérifier cette
hypothèse en regardant une image post – mortem obtenue par radiographie X.
152
(C , ! - ,
Dans ce drapage, on observe que la zone délaminée coïncide avec la zone fortement
fissurée, ce qui est en accord avec la littérature [143] car les fissures matricielles provoquent
des micro @ délaminages en pointe de fissures qui coalescent pour former des zones macro
– délaminées. En ce qui concerne les autres drapages, on a vu sur la Figure 92 que les
éprouvettes contenant les plis à ±30º montrent un endommagement étendu le long de
toute l’éprouvette avec des zones complètement rompues par cisaillement (déchirées). On
voit bien cet effet lorsqu’on fait des analyses par rayons X (Figure 95).
Le drapage quasi – isotrope ne montre pas de zones déchirées comme dans le cas [02,
±302]s. Par contre, la présence d’un pli central à 90° par rapport à la charge fait apparaitre
de nouveaux types de fissures. Parmi ces fissures, on remarque des fissures intra –
laminaires au sein du pli central indiqués par les flèches bleues dans les encadrés en
pointillé rouge (Figure 96).
153
(C , ! - ,
1mm
a)
1mm
b)
Figure 96. Fissures intra @ laminaires dans le pli à 90º. a) entaille b) hors entaille.
Ces fissures intra – laminaires (flèches bleues) apparaissent dans l’entaille (Figure 96a)
et aussi sur les bords (Figure 96b). L’orientation observée est d’environ 45º par rapport à la
normale du plan du pli. Ces fissures connectent toujours des zones fortement délaminées.
L’observation par rayons X de ces éprouvettes nous donne le résultat suivant (Figure 97):
7.5mm
Comme on l’avait déjà constaté par des observations micrographiques, on vérifie que
le délaminage (ligne rouge) est moins développé dans ce drapage que dans le [02, ±302]s.
Par contre, grâce à la radiographie, on est capable d’observer une fissuration matricielle des
plis à ±45º assez importante autour de l’entaille (zone délimitée par la ligne bleue). On
remarque aussi le caractère non symétrique du délaminage et de la fissuration, beaucoup
plus forte à gauche de l’entaille. En ce qui concerne le pli orienté à 90º, il y a une forte
densité des fissures dans l’entaille, ce qui est tout à fait normal, mais on trouve aussi des
154
(C , ! - ,
fissures très loin de l’entaille (flèches bleues). Ces fissures éloignées amorcent sur le bord et
se propagent vers l’intérieur sans arriver à traverser toute l’éprouvette.
Alors que ce chapitre a présenté jusqu’à présent les modes d’endommagement des
éprouvettes entaillées à température ambiante, les lignes qui suivent abordent la question de
l’effet de la température sur le développement de l’endommagement. En effet, les essais en
température (même s’ils n’entreront pas dans la base de données expérimentale visant à
valider le modèle de comportement élaboré au chapitre précédent) vont apporter des
informations sur la compétition entre certains types de dégradation (notamment le
délaminage et la fissuration matricielle) et la viscosité de la matrice. Grâce à la technique
des rayons X, on peut comparer d’une façon qualitative les quantités de délaminage et de
fissuration présentes dans les éprouvettes. On montrera des comparaisons des
radiographies post – mortem des éprouvettes testées à @40ºC, 20ºC et 85ºC pour les
drapages [02, ±202]s et [02, ±452, 902]s. L’absence de radiographies du drapage à 30º est due
au fort endommagement de cette éprouvette. On commence d’abord en analysant le
drapage [02, ±202]s.
040ºC
20ºC
7.5mm
20ºC
85ºC
Figure 98. Comparaison des radiographies à trois températures [0 2 , ±20 2 ]s
155
(C , ! - ,
040ºC
7.5mm
20ºC
85ºC
Figure 99. Comparaison de radiographies à trois températures des éprouvettes quasi @ isotropes.
156
(C , ! - ,
La différence la plus notable entre les trois températures est le degré de fissuration
matricielle dans le pli central (90º). On constate que plus la température est élevée, plus le
nombre des fissures matricielles dans le pli central est faible, c’est qui est en accord avec
l’hypothèse de compétition entre l’endommagement matriciel et la viscosité. L’éprouvette à
@20ºC montre de l’endommagement matriciel dans toute sa partie utile avec des fissures qui
s’amorcent aux deux bords libres. On trouve aussi des fissures qui traversent l’éprouvette
loin des entailles (flèches bleues). On perd la symétrie de l’endommagement quand la
température augmente. On peut observer que les fissures se concentrent à gauche de
l’entaille (à 20°C) et surtout sur le bord qui subit le macro délaminage le plus fort, même s’il
reste encore quelque fissures isolées de l’autre côté (flèche verte). En augmentant encore
plus la température, on observe que l’endommagement est complètement concentré autour
de l’entaille et que la symétrie de l’endommagement (délaminage et fissuration) est de
nouveau présente. En ce qui concerne le délaminage, on trouve une forte concentration en
forme de X autour de l’entaille produite par le délaminage entre les plis croisés et aussi
entre ces derniers et le pli central. Comme dans les autres éprouvettes, les couches à 0º sont
totalement délaminées dans toute la longueur de l’éprouvette. Ces couches étaient tellement
séparées qu’on n’a pas pu imprégner correctement avec le liquide opacifiant.
80?D - 4 0@D - 4
157
(C , ! - ,
>AD - 4 @A0 D - 4
Figure 100. Champ de déformation de cisaillement d'une éprouvette quasi – isotrope a différents
moments du chargement.
On peut voir sur cette figure que le cisaillement se concentre dans les endroits où le
splitting apparait (Figure 89). L’état de cisaillement reste symétrique dans les deux
premières images, correspondant à 9.6% et 28.5% de la charge maximale. La troisième
image (47%), correspondant à la moitié du chargement, montre une légère concentration
de déformation dans l’entaille gauche. A 57.3% on voit clairement une forte déformation
localisée dans les points où le splitting est en train d’apparaitre.
Groupe II
158
(C , ! - ,
Une vue latérale de l’entaille montre la présence des fissures matricielles (en vert) dans
tous les plis ainsi que du délaminage, très concentré dans le centre de l’entaille. On se
servira des radiographies pour mieux étudier ce délaminage et l’influence de la température
sur la cinétique de l’endommagement observé. Dans la Figure 102, on peut reconnaitre
certaines tendances vues dans les éprouvettes du groupe I, surtout l’augmentation de la
fissuration des plis croisés avec la température. A basse température, les fissures sont
concentrées dans l’entaille avec quelques@unes qui ont traversé toute l’éprouvette (flèches
bleues). Le délaminage n’est pas très présent à cette température, dans la mesure où la
fissuration matricielle brusque a provoqué la rupture de l’éprouvette avant l’apparition du
délaminage, en raison de la cinétique très rapide provoquée par la fragilisation de la matrice
à basse température.
159
(C , ! - ,
040ºC
20ºC
85ºC
Figure 102. Radiographies post @ mortem des éprouvettes [±45 4 ] s a trois températures.
160
(C , ! - ,
17.5% 35%
51% 69%
Figure 103. Corrélation d’images sur des éprouvettes [±45 4 ] s à différentes étapes du chargement
exprimées en fraction de la charge à rupture.
Groupe III
Les éprouvettes du Groupe III, avec la séquence [±67.54]s, ne montrent pas trop de
différence entre les trois températures testées. On trouve toujours une rupture matricielle
très nette avec une zone déchirée post – rupture très concentrée autour de la zone de
rupture. La raison de la prépondérance de ce mode est la forte sollicitation de l’interface
fibre/matrice dans ce drapage. On montre dans la Figure 104 les trois radiographies
correspondant aux trois températures.
161
(C , ! - ,
040ºC
20ºC
85ºC
Figure 104. Radiographies post @ mortem des éprouvettes [±67.5 4 ] s a trois températures.
Il y a très peu de différences entre ces trois clichés. On peut voir quelques fissures
isolées qui s’amorcent loin de l’entaille à froid (flèche bleue), peut@être en raison du
caractère plus fragile de la matrice à @40ºC. Une autre fissure isolée est repérée à
température ambiante dans l’entaille, proche de la zone de rupture finale. Enfin, à haute
température, on ne trouve pas de fissures isolées. La zone de rupture finale semble être un
plus étendue que dans les autres cas. Il semble que l’interface fibre/matrice (fortement
sollicité dans cette séquence) n’est pas très sensible à la température.
162
(C , ! - ,
Comme on l’a décrit dans le premier chapitre, l’émission acoustique est une méthode
par laquelle on est capable de dire à quel moment et à quel endroit l’endommagement
apparait et, en utilisant les autres données déjà présentées, être capable de distinguer
(qualitativement) les modes d’endommagement impliqués. Dans cette étude, on va
présenter systématiquement des cumuls d’énergie absolue, ce qui nous permettra de
distinguer plus facilement des échos appartenant à chaque mode d’endommagement. On se
servira aussi des histogrammes d’évènements dans le temps et l’espace pour montrer l’effet
de l’entaille comme concentrateur d’endommagement. On va suivre la même démarche que
celle présentée dans la section précédente : d’abord, les éprouvettes du Groupe I seront
analysées à l’ambiante, puis on étudiera l’influence de la température. Après ce groupe, on
présentera les résultats correspondant aux groupes II et III.
x 10
4 [02 , ±202]s x 10
8
4 2
3 1.5
Load (N)
2 1
1 0.5
0 0
0 0.002 0.004 0.006 0.008 0.01 0.012 0.014
Axial strain
a) b)
163
(C , ! - ,
On peut voir sur la Figure 105a la courbe de cumul d’énergie acoustique (en vert) en
superposition avec la courbe de traction (en bleu) d’une éprouvette [02, ±202]s. On constate
qu’on commence à enregistrer des signaux presque dès le début (flèche rouge) même si on
ne voit pas d’effets à l’échelle macroscopique sur la courbe de traction. Puis, à un niveau de
déformation de 0.004, la pente de la courbe acoustique change brusquement ce qui indique
l’apparition d’un nouvel endommagement, très probablement le début de la rupture
matricielle des plis à 20º proche de l’entaille. Cette courbe est interrompue par des sauts
d’énergie intermittents. Ces sauts correspondent à des évènements très énergétiques. C’est à
ce moment@là que le splitting apparaît. Il semble donc raisonnable de considérer ce
phénomène comme l’émetteur de ces signaux. Chaque petit saut correspond à une
propagation du splitting. On note également une augmentation de la pente de la courbe
acoustique, témoin d’une multiplication des sites d’évolution de l’endommagement. Ce
phénomène est confirmé par l’histogramme 3D (Figure 105b), qui montre la concentration
des évènements dans le centre (avec un petit décalage) de l’éprouvette. Si nous regardons
cet histogramme par@dessus, on pourra voir une forme de triangle dans le plan temps –
position (Figure 106).
a) b)
Figure 106. a) Vue de dessus de l'histogramme 3D. b) Augmentation de la zone endommagée au
cours du temps.
Cette figure montre qu’à partir de la moitié de l’essai, l’aire du matériau qui émet des
échos commence à augmenter et la pente du cumul acoustique augmente par voie de
conséquence.
164
(C , ! - ,
La séquence [02, ±302]s montre un état post – mortem plus endommagé avec des plis à
0º complètement délaminés et des plis croisés rompus par cisaillement (déchirés).
4
x 10
[02 , ±302]s x 10
8
2 2
Load (N)
2 4 6 8 10 12 14 16
Axial strain x 10
@3
a) b)
Figure 107. a) Cumul d'énergie (vert) et courbe de traction (bleu) et b) histogramme 3D d'une
éprouvette [0 2 , ±30 2 ] s .
165
(C , ! - ,
x 10
4 [02 , ±452 , 902]s x 10
7
2.2 4.5
2 4
1.8 3.5
1.6 3
Load (N)
1.4 2.5
1.2 2
1 1.5
0.8 1
0.6 0.5
0.4 0
2 4 6 8 10
Axial strain @3
x 10
a) b)
Figure 109. a) Cumul d'énergie acoustique et b) histogramme 3d d'une éprouvette [0 2 , ±45 2 , 90 2 ] s .
Dans cette séquence, le changement de cinétique est manifeste par rapport aux
courbes correspondant aux éprouvettes [02, ±202]s et [02, ±302]s. La seule différence entre
cette séquence et les précédentes est la présence du pli à 90º qui va provoquer l’apparition
de la fissuration matricielle dans ce pli. La première partie de la courbe (augmentation
régulière de l’activité acoustique et légère non linéarité de la courbe force – déformation)
est due à cette fissuration matricielle à laquelle s’ajoute une contribution d’endommagement
diffus provenant des plis croisés. Le macro délaminage apparaît peu avant la charge à
rupture, ce qui provoque des sauts dans la courbe. Ces sauts sont plus progressifs que dans
166
(C , ! - ,
les cas précédents, ce qui pourrait indiquer une propagation moins brusque de cet
endommagement. Par rapport à l’histogramme 3D, on peut observer que la zone de
l’entaille concentre des signaux jusqu’à l’apparition du premier évènement acoustique
énergétique puis on note une extension vers le voisinage.
[02 , ±202]s
7
x 10
14
@ 40ºC
12 85ºC
20ºC
10
Abs Energy
0
0 0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016
Axial strain
167
(C , ! - ,
reçue. Par exemple, si nous regardons la courbe acoustique correspondant au drapage quasi
– isotrope, on obtient :
8 [0 , ±45 , 90 ] 7
x 10 2 2 2 s x 10
9
9
8 @ 40ºC @ 40ºC
8 85ºC
85ºC
7 20ºC
20ºC 7
6 6
Abs energy
5 5
4
4
3
3
2
2 1
1 0
Figure 111. Comparaison cumul d'énergie du drapage quasi @ isotrope et zoom de la partie basse
énergie.
Le drapage quasi – isotrope montre une grande différence entre l’essai à froid et celui à
85ºC mais cette fois@ci l’essai à température ambiante est plus proche de l’essai à chaud. Il
faut regarder dans la zone de très basse énergie pour s’apercevoir qu’à 20ºC l’activité se
déclenche en même temps qu’à @40ºC mais les énergies libérées n’ont pas le même niveau.
La différence entre ces deux drapages ([02, ±202]s et quasi – isotrope) est le mode
d’endommagement mis en jeu avec un rôle plus marqué de la fissuration matricielle et de
l’endommagement diffus dans l’éprouvette quasi – isotrope. C’est l’analyse des éprouvettes
du Groupe II qui nous donnera quelques pistes sur ces différences. On montre ci@dessous
la comparaison du comportement mécanique et acoustique et l’histogramme 3D
correspondant :
[±454]s x 10
6
6000 15
4000 10
Load (N)
2000 5
0 0
0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09
Axial strain
a) b)
Figure 112. a) Comportement mécanique et acoustique et b) histogramme 3d d'une éprouvette
[±45 4 ]s.
168
(C , ! - ,
7 [±45 ]
4
x 10 4 s
@ 40 ºC
85ºC
3
20ºC
Abs Energy
0
2 4 6 8 10 12
Axial strain
Cette figure permet de corroborer l’hypothèse que les plis à ±45º et l’endommagement
associé ont une forte influence puisque la signature acoustique enregistrée en fonction de la
température est très similaire à celle qu’on a retrouvée dans le drapage quasi – isotrope
(dont la dépendance à la température est différente de celle des drapages à 20º et 30º). Il est
clair que soit l’endommagement diffus, soit le glissement, soit la rupture matricielle, ou
l’ensemble de ces mécanismes joue un rôle très important dans les signatures acoustiques.
Cette conclusion souligne l’intérêt des tests réalisés sur une séquence quasi@isotrope, dans
169
(C , ! - ,
4.4 Conclusions
170
(C , ! - ,
• Groupe I : les éprouvettes qui contiennent des plis à 0º où la rupture des fibres
pilote la rupture finale mais avec une forte présence de délaminage.
171
(C , ! - ,
On montrera dans le chapitre suivant que cette ondulation a une influence sur la
cinétique du délaminage et même sur la rupture des fibres.
172
CHAPITRE V:
VALIDATION DU MODELE
ET
APPLICATIONA LA TENUE DE
RESERVOIRS HYPERBARES
C,C 4 .
5.1 Introduction
Ce chapitre est composé de deux parties bien différenciées: la première partie porte
sur la confrontation du modèle d’endommagement avec les données issues de certains
essais de traction sur éprouvettes entaillées montrés dans le chapitre précédente. Cette
procédure de validation nous permettra de vérifier les capacités du modèle par différentes
voies comme la comparaison des courbes de traction (information macroscopique,
présentant la force en fonction du déplacement enregistré entre les marqueurs déposés sur
l’éprouvette) et des charges à rupture, la comparaison des cartographies d’endommagement
avec des micrographies post mortem pour étudier les types et la localisation des
endommagements. Des radiographies (aussi post – mortem) nous permettront de
comparer l’extension des endommagements dans le plan de l’éprouvette. On se servira
aussi de l’émission acoustique virtuelle pour confronter modèle et expérience.
Une fois ce modèle validé, il sera utilisé pour étudier le comportement, notamment
l’éclatement, d’une structure complexe fabriquée par enroulement filamentaire : le réservoir
type IV haute pression. Cette structure sera d’abord présentée en détaillant sa composition,
le processus de fabrication et ses particularités. Puis, la modélisation par EF de cette
structure sera décrite en détail : géométrie, variation d’orientation dans le dôme, conditions
aux limites etc… Le cœur de cette partie sera la confrontation des résultats expérimentaux
et numériques en termes de comportement mécanique, de prédiction de l’éclatement et
d’autres phénomènes comme la fuite du liquide de pressurisation. Une étude comparative
entre le modèle déterministe et probabiliste sera faite à la fin du chapitre.
Dans cette première partie du chapitre, on exploitera les essais de traction effectués
sur des éprouvettes entaillées, et dont les résultats ont été présentés dans le chapitre
précédent, en les comparant aux résultats de simulations réalisées avec le modèle
d’endommagement développé au Chapitre 3. Notre but est d’analyser ses performances, ses
points forts et faibles, pour pouvoir poursuivre son amélioration dans le cadre des études
futures. La simulation du comportement de ces éprouvettes (et de leur état de contraintes
complexe) constitue un exemple de calcul de structures suffisamment contraignant pour
valider, ou non, la capacité du modèle à capter cet état et les différents endommagements
associés aux différentes séquences d’empilement.
176
C,C 4 .
[Link] Géométrie.
On rappelle la géométrie des éprouvettes entaillées qui seront modélisées par Eléments
Finis dans le logiciel commercial ABAQUS®. Il ne faut pas oublier que notre but ultime
est de développer un outil qui soit utilisable par des Bureaux d’Etudes et donc, aisément
intégrable dans un logiciel de calcul commercial.
L’épaisseur par pli est de 0.22mm environ. On a construit deux types de modèles
différents pour pouvoir utiliser les deux versions du modèle développées :
On montre ci –dessous les géométries résultantes dans les deux cas mentionnés :
a)
177
C,C 4 .
b)
Figure 116. a) Modèle EF coque et b) 3D d’une éprouvette [±45 4 ]s
On peut voir dans les deux cas la présence de l’entrecroisement pour essayer de rendre
le plus réaliste possible nos simulations. La démarche suivie est identique à celle expliquée
dans le Chapitre 3 pour la simulation des éprouvettes droites lors du processus
d’identification. Les orientations ont été inversées entre deux sections adjacentes pour
représenter les conséquences du procédé de bobinage. De cette manière, on retrouvera une
alternance des orientations [+θ,@θ, +θ,@θ…] et [@θ, +θ, @θ, +θ…] entre les sections
adjacentes. Cette technique a été utilisée par exemple dans les travaux de [144] pour
comparer les mesures DIC avec des calculs EF. En ce qui concerne le modèle complet
(3D), on doit remarquer que les interfaces n’ont pas été modélisées, comme dans le cas des
travaux de Cachan [79] ou [145], et toute l’éprouvette est un solide continu. De cette façon,
on gagne en temps de calcul et en simplicité du modèle.
La présence des éléments triangulaires (S6R) ou prismatiques (C3D15) est imposée par
la complexité de la géométrie et notamment par la présence des entrecroisements. Tous les
éléments utilisés sont quadratiques, à intégration réduite pour les éléments coque, complète
pour les éléments volumiques. Des tests de convergence du résultat ont été faits, d’abord
en utilisant des calculs élastiques pour étudier la convergence en contraintes, puis pour
calibrer les paramètres de régularisation visqueuse. Ces paramètres ont été identifiés de telle
178
C,C 4 .
sorte que la zone endommagée n’évolue pas en fonction de la taille de maille. On n’attend
pas une convergence totale en comportement mais la dépendance au maillage est fortement
réduite. Des images des maillages résultants sont présentées sur la Figure 117 :
a) b)
Figure 117. Maillage a) coque et b) 3D.
Pour le modèle 3D, on impose un élément par pli dans l’épaisseur. Ce choix, sur lequel
il sera peut@être nécessaire de revenir, est un compromis entre un nombre d’éléments
raisonnable (et des temps de calcul acceptables) et une discrétisation suffisamment fine de
l’épaisseur. La figure suivante présente une vue de l’entaille du modèle EF 3D de
l’éprouvette quasi – isotrope.
179
C,C 4 .
de se concentrer sur les empilements les plus complexes et les plus représentatifs (par
l’activation de nombreux modes d’endommagement) du comportement du bobiné.
Simulation [±454]s
7000
6000
5000
Force (N)
4000
3000 Test 1
2000 Test 2
Simulation Shell
1000
Simulation 3D
0
0 2 4 6
Déplacement (mm)
Figure 120. Courbes de traction numérique et expérimentale de l'éprouvette entaillée [±45 4 ] s .
180
C,C 4 .
Un autre paramètre macroscopique peut être comparé directement avec les données
obtenues lors des essais : l’énergie cumulée acquise par le système de mesure de l’émission
acoustique. Grâce aux émissions virtuelles implémentées dans le modèle
d’endommagement (chapitre 3 section 3.4.4 « Emission Acoustique Virtuelle), on est
capable de reproduire numériquement le phénomène d’émission acoustique (EA) et
d’améliorer la compréhension des données expérimentales via la comparaison avec l’EA
numérique. Pour pouvoir comparer ces deux signaux (numérique et expérimental), on a
effacé les données correspondant à la fin d’essai (émissions acoustiques post – rupture) car
ils n’apportent pas d’information sur l’amorçage et le développement de l’endommagement
et introduisent des perturbations importantes. En outre, on normalisera le temps d’essai et
l’énergie cumulée. De cette façon, on peut aisément comparer les allures des courbes, ce
qui nous permettra d’identifier et de commenter les changements de pentes ou les rapports
entre les énergies des différents modes d’endommagement. La comparaison est montrée
dans la Figure 122 ci – dessous.
181
C,C 4 .
Figure 122. Comparaison entre les mesures d'émission acoustique virtuelle (sim) et réelle (exp).
On remarque plusieurs détails intéressants dans cette figure. En premier lieu, les
débuts de l’émission acoustique coïncident dans les deux cas (flèche verte). Cela signifie que
l’identification de l’endommagement diffus (celui qui apparaît en premier dans cette
séquence) à partir de la seule courbe de diminution de rigidité est bien adaptée à la réalité.
On constate aussi que les changements de pente de la courbe d’EA sont assez bien
reproduits par le modèle (flèches rouges et bleues). Il faut noter que, dans le modèle, c’est
l’endommagement diffus qui génère ces émissions. Les changements de pente successifs
(alors qu’aucun nouveau mode d’endommagement n’est activé) s’expliquent par le type
d’évolution de cette variable d’endommagement (on rappelle que sa valeur tend vers une
asymptote). Dans la partie a de la courbe, seule une petite partie de l’échantillon autour de
l’entaille est endommagée et donc émet des signaux énergétiques. Quand cette partie est
proche de la saturation, d’autres zones de l’éprouvette situées au@dessus et en@dessous des
entailles commencent à émettre, ce qui marque le début de la partie b. Ainsi, les
changements de pente des courbes d’énergie cumulée peuvent témoigner d’une
augmentation de la taille de la zone endommagée et pas seulement de l’apparition d’un
autre type d’endommagement. La fin de la partie b sur la courbe simulée est marquée par
la présence de deux sauts liés au début de la rupture matricielle. La simulation de ces sauts
sont validés par la présence d’un saut au même moment sur la courbe expérimentale, avant
la croissance brutale de l’énergie.
182
C,C 4 .
Mais les capacités du modèle développé ne se limitent pas à fournir des résultats
macroscopiques comme les courbes force vs déplacement ou l’émission acoustique.
L’intérêt de l’approche réside dans la possibilité de suivre localement certaines grandeurs.
Par exemple, la Figure 123 compare le champ de déformation axiale du modèle avec celui
obtenu par DIC :
2017N
3000N
Figure 123. Comparaison du champ de déformation axiale experimental (gauche) et numerique
(droite).
On voit dans la Figure 123 que les champs de déformation axiale simulés sont assez
proches de la mesure expérimentale en ce qui concerne les valeurs maximales mais plus
éloignés pour les valeurs basses. La concentration de déformation en forme de triangle est
très bien reproduite par le modèle, mais, même si la corrélation d’image indique la présence
d’une zone de faible déformation placée au centre de l’éprouvette, le niveau de déformation
simulé dans cette même zone est supérieur (4.5e@3 contre 7.2e@3).
Toujours à l’échelle locale, le modèle est capable d’indiquer les niveaux des différents
d’endommagement. Par exemple, on peut observer les cartographies d’endommagement
183
C,C 4 .
diffus et de rupture matricielle pour pouvoir comparer la zone de rupture simulée avec celle
apparaissant sur les éprouvettes réelles. La Figure 124 montre la cartographie
d’endommagement diffus à plusieurs moments différents.
Figure 124. Niveau d'endommagement diffus dans les trois phases (a, b, c) définies sur la courbe
d'EA (Figure 122).
On voit sur la figure précédente les trois états correspondant aux différentes parties de
la courbe d’énergie acoustique. Dans le premier stade, l’endommagement diffus est proche
de la saturation dans les entailles. Une fois cette saturation atteinte, l’endommagement
commence à s’étendre vers le reste de l’éprouvette au cours du stade b. La dernière image,
qui montre un état plus homogène que les précédentes, est prise juste avant l’apparition de
la première rupture matricielle qui provoquera la ruine totale de l’éprouvette. Cette rupture
commence dans les zones très chargées des entailles et se propage à ±45º. Les zones de
rupture expérimentale et numérique (modèles coque et 3D) sont montrées dans la figure ci@
dessous.
a) b) c)
Figure 125. a) Rupture expérimentale, b) Modèle 3D (gauche : endommagement hors plan,
droite : endommagement matriciel), c) coque (endommagement matriciel). La couleur noire
représente la rupture totale de la matrice.
184
C,C 4 .
Cette figure montre d’abord que les simulations 2D et 3D prédisent la même zone de
rupture. Dans la réalité, la rupture matricielle choisit un chemin possible dans l’ensemble de
la zone prédite par le modèle. Même si les simulations de la Figure 125b et c sont très
proches en ce qui concerne l’endommagement matriciel, le modèle 3D fournit une
estimation plus fine de l’amorçage et de la propagation de l’endommagement, notamment
en donnant des indications sur le mode hors@plan. La partie gauche de la Figure 125b
indique la configuration d’endommagement hors – plan dans une vue latérale de l’entaille
du modèle 3D. Le modèle 3D prédit une rupture hors – plan plus étendue entre les plis
externes que dans l’intérieur de l’entaille (flèche blanche). Cette zone étendue coïncide avec
la formation d’une fissure et d’un délaminage au cours des essais (flèche blanche sur la
Figure 125a).
Une autre donnée fournie par le modèle est visible sur une coupe de l’éprouvette
(Figure 126).
Figure 126. Endommagement hors @ plan a l'intérieur de l'éprouvette juste avant la rupture finale
(Noir = Totalement endommagé).
Le rectangle rouge marque une zone d’endommagement hors – plan située au milieu
de l’éprouvette à l’endroit où le niveau de déformation axiale est minimal (voir Figure 123).
C’est dans cette zone qu’on trouve les entrecroisements (changements d’orientations dans
le modèle). Ces changements brusques d’orientation provoquent des états de contraintes de
cisaillement hors plan assez forts à la jonction des plis ±45° et ce sont ces contraintes de
cisaillement qui favorisent l’apparition de l’endommagement hors – plan. Le modèle
semble donc prédire une initiation de la rupture hors@plan au cœur de l’éprouvette, au
niveau des entrecroisements. Cette hypothèse n’a pas pu être vérifiée car cette partie se
déchire trop rapidement pendant la rupture finale.
185
C,C 4 .
Dans cette section, on va étudier l’éprouvette qui met en jeu la plus grande variété de
modes d’endommagement : l’éprouvette [02, ±452, 902]s, qu’on qualifie de « quasi –
isotrope ». Simuler son comportement est un test fortement discriminant, sur lequel on
évalue les performances du modèle pour pouvoir reproduire les phénomènes observés
pendant les essais, avant de l’utiliser pour simuler des structures beaucoup plus complexes
comme le réservoir bobiné. On suit la même démarche que celle de la section précédente
en analysant en premier des données macroscopiques comme la courbe de traction et
l’émission acoustique, puis on descendra à une échelle inférieure pour étudier les
différentes formes d’endommagement. On montre dans la suite les courbes correspondant
aux simulations 2D et 3D déterministes. Il faut préciser que l’empilement quasi@isotrope
contient des plis à 0º et, par conséquent, subit des ruptures de fibres. Il est donc légitime
d’utiliser la démarche probabiliste développée dans le chapitre 3, ce qui sera fait dans un
deuxième temps. On étudie en premier lieu la qualité des prédictions déterministes (Figure
127).
20000
Force (N)
15000
10000 Test
Simulation 2D
5000
Simulation 3D
0
0,00 0,20 0,40 0,60 0,80 1,00
Déplacement (mm)
On voit dans cette éprouvette que les deux versions du modèle (2D vs. 3D)
reproduisent très bien le comportement de l’éprouvette quasi – isotrope. Les charges à
rupture sont assez proches (Figure 128), même si la prédiction 3D montre l’erreur la plus
petite et la plus conservative.
186
C,C 4 .
Simulation 2D sta_3
Simulation 2D sta_4
15000
10000
5000
0
0,00 0,20 0,40 0,60 0,80
Déplacement (mm)
187
C,C 4 .
18000
Sta – 4 :@3
17000 Moyenne Sta 3@
Test Sta 0 1 Sta 0 2 Sta 0 3 Sta 0 4
On voit bien dans cette comparaison l’effet de l’incorporation d’une variabilité dans la
contrainte à rupture. La présence ou non d’un élément faible provoque les variations de
charge à rupture observées. En calculant la charge à rupture moyenne, on réussit à obtenir
une erreur de 0.51%, ce qui est tout à fait satisfaisant.
0,8 EA Sim
normalisée
0,6 EA Exp
0,4
0,2
0
0 0,2 0,4 0,6 0,8 1
Temps normalisé
Figure 131. Comparaison énergie cumulée expérimentale et simulée (séquence quasi @ isotrope).
On trouve aussi un bon accord entre l’énergie acoustique simulée et celle enregistrée
dans l’essai de traction. A partir de 0.4, la pente de l’énergie enregistrée augmente en raison
de la contribution de l’endommagement matriciel dans les plis à 90º et de l’extension de
l’endommagement diffus dans les plis croisés. Même si la pente de cette partie est
légèrement surestimée par le modèle, on arrive à bien capter le début puis l’évolution des
188
C,C 4 .
endommagements énergétiques (rupture des fibres, cercle vert) qui conduisent à la ruine de
l’éprouvette.
Figure 132. Comparaison fissuration matricielle post @ mortem expérimentale et numérique (noir
= complètement endommagé).
189
C,C 4 .
entailles. Le modèle prédit aussi une densité de fissures non nulles dans les zones éloignées
de l’entaille. Même si les radiographies témoignent bien de la présence de cette fissuration à
45° en dehors de la zone d’influence de l’entaille, le modèle semble en surestimer l’étendue.
En ce qui concerne la fissuration purement transverse, le pli à 90º (simulation) montre une
zone très fissurée dans le centre de l’éprouvette et une zone qui s’étend jusqu’au mors.
Dans cette zone éloignée, les bords montrent une couleur un peu plus foncée ce qui
signifie une densité des fissures plus élevée. Les radiographies confirment la pertinence des
prédictions dans ce pli puisque le centre de l’éprouvette est très fissuré mais on trouve
également des fissures horizontales loin de l’entaille. Ces fissures s’amorcent sur le bord
libre et continuent vers l’intérieur ce qui est en accord avec le gradient de densité des
fissures donné par le modèle.
190
C,C 4 .
On voit dans ces images que la zone délaminée simulée se concentre autour de
l’entaille. Il est intéressant de remarquer que le délaminage apparaissant entre les plis à 0º et
croisés d’une part, et 90° et croisés d’autre part, est beaucoup plus marqué que celui à
l’interface des plis croisés. Ce fait sera mis en évidence dans la suite, sur les micrographies
de la tranche de l’entaille. L’endommagement hors plan à l’interface du pli à 0º et du pli à
45° est important dans les coins des entailles, comme on peut également le repérer sur les
éprouvettes réelles (cercles rouges), mais aussi dans une zone elliptique au centre de
l’éprouvette. Dans cette zone, on remarque une partie non endommagée au centre de
l’éprouvette, similaire à celle observée sur les éprouvettes à ±45º. C’est vraisemblablement
la différence des déformations entre les plis à 0° et 45° qui provoque cette forme de
délaminage au centre. La présence d’un fort délaminage entre les plis croisés et le pli à 90°
peut également s’expliquer en remarquant que le pli à 90° subit une déformation dans la
direction latérale très faible (en raison de la faible valeur de son coefficient de Poisson), à la
différence des plis à ±45°, ce qui favorise la rupture de l’interface entre ces deux plis
(Figure 134). Cette figure montre aussi la tendance des plis à 0° à se séparer du reste de
l’éprouvette.
Figure 134. Déformée (grandissement x5) du centre de l'eprouvette entaillee quasi @ isotrope.
On note que l’étendue des délaminages expérimentaux est beaucoup plus grande que
celle simulée. Une explication de cet écart peut être reliée au fait que certains plis ne sont
pas plans. En effet, grâce à des observations micrographiques d’une autre séquence testée
pour le projet OSIRHYS IV ([02, ±452, 02]s), on a pu mettre en évidence une ondulation
systématique des plis adjacents aux plis croisés (Figure 135b). Le calcul apporte une
réponse quant à l’influence de ces ondulations sur l’amorçage et la propagation du
délaminage et même sur l’apparition des ruptures des fibres prématurées : on a réalisé une
simulation de traction sur un échantillon de matériau composite dans lequel on a modélisé
191
C,C 4 .
a b
Figure 135. Composite bobiné ondulé (a) modèle EF et b) micrographie)
Lorsque les plis ne présentent pas d’ondulation, l’état de contrainte local est
homogène. Les contraintes dans le sens des fibres ainsi que l’endommagement hors – plan
prédit par le modèle dans le cas de cette éprouvette sont analysés. .La Figure 136 représente
les contraintes locales dans la direction des fibres de cette séquence. La valeur maximale
σfmax des plis ondulés est comparée à la valeur σfmoy des plis non ondulés.
Figure 136. Endommagement hors plan (en haut) pour différentes valeurs de contraintes dans le
sens des fibres (en bas).
192
C,C 4 .
Les lignes qui précèdent ont analysé les cartographies d’endommagement par
comparaison avec les radiographies post – mortem. On souhaite, dans la suite, exploiter la
possibilité d’observer directement le fond des entailles pour y distinguer les différents types
d’endommagement dans l’épaisseur. Pour une meilleure visualisation de ces
endommagements, les fissures intralaminaires sont colorisées en vert et les fissures
interlaminaires en rouge. La Figure 137 montre une comparaison entre la fissuration
matricielle simulée et une micrographie après rupture.
Les deux carrés (noir et rouge) signalent une zone de fissuration matricielle dans le pli
central à 90º qu’on retrouve tant dans la micrographie que dans la simulation. On voit que
193
C,C 4 .
la partie droite du pli à 90° est davantage fissurée que la partie gauche. Cette dissymétrie se
retrouve sur la simulation En ce qui concerne la fissuration matricielle des plis croisés, on
constate (agrandissement de la micrographie à droite, avec des flèches bleues) que le pli
croisé directement adjacent au pli à 90° n’est pas fissuré alors que le pli suivant subit de la
fissuration. Cette alternance est retrouvée sur la simulation à droite : la zone indiquée par
les flèches bleues (zone de fort endommagement dans un pli croisé) est encadrée par deux
zones sans endommagement. Cette alternance est anti – symétrique (par rapport à un axe
vertical passant par le centre de l’entaille), tant dans la simulation que sur les micrographies,
puisque de l’autre côté de l’entaille (flèches bleues à gauche), on trouve la fissuration
matricielle dans les plis qui n’étaient pas fissurés et inversement.
194
C,C 4 .
a) 30% b) 30%
Figure 140. Champ de glissement plan γ 12 a) expérimental et b) numérique, pour une valeur de
force axiale appliquée de 30% de la charge maximale dans les deux cas.
Le modèle retrouve bien le même état de glissement (avec des ordres de grandeur
comparables), à l’origine de l’apparition de l’endommagement diffus et de la rupture
matricielle qui provoque le splitting.
Avec ces deux analyses de séquences différentes, on a montré les capacités du modèle
à reproduire le comportement macroscopique de ces structures entaillées via la simulation
de courbes charge – déplacement ou en utilisant une méthode innovante pour comparer
l’émission acoustique expérimentale avec celle implémentée numériquement. Cette
comparaison nous permet d’approfondir l’interprétation des données issues d’essais et
fournit des pistes afin de mieux calibrer certains paramètres du modèle (seuils, cinétique…).
195
C,C 4 .
Les deux versions du modèle (2D – 3D) donnent des résultats très satisfaisants. Cela nous
permet d’envisager une méthode de calcul dans laquelle les modèles 2D et 3D seraient
complémentaires. Les calculs de structures de grande taille seraient réalisés avec la version
2D (beaucoup plus économe en temps de calcul) et le modèle 3D serait utilisé pour des
études beaucoup plus fines dans les zones critiques. Cette possibilité de passer d’un modèle
à l’autre et la richesse des simulations 3D nous permettent de descendre à une échelle fine
pour aller étudier des modes d’endommagement au niveau du pli. Le modèle a montré sa
capacité à discerner les différents modes d’endommagement ainsi que leur localisation dans
le plan et dans l’épaisseur avec un très bon accord expérimental. On remarque la présence
d’asymétrie dans la fissuration transverse du pli central de la séquence quasi – isotrope qui
est capté par le modèle avec précision ainsi que l’alternance des plis croisés fissurés. La
présence prédominante du délaminage entre les plis à 0º et ±45º ou 90º et ±45º est aussi
très bien captée par le modèle même loin de l’entaille dans le cas du couple 90º/±45º où la
présence de fissuration matricielle dans le bord provoque l’amorçage du délaminage à cet
endroit. On met ainsi en évidence le couplage « indirect » entre ces deux endommagements.
Les outils introduits dans le modèle comme la rupture probabiliste permettent d’obtenir un
comportement très proche des éprouvettes réelles, ce qui nous rapproche du concept de
« virtual testing ». On trouve également une très bonne cohérence entre les données
fournies par l’EA virtuelle et celles issues de l’expérience. Cela nous ouvre une nouvelle
voie d’étude pour identifier les cinétiques des différents endommagements à l’aide de
l’émission acoustique ainsi que mieux comprendre les données expérimentales grâce au
modèle.
Dans cette deuxième partie du chapitre, on étudiera une structure beaucoup plus
complexe que celles étudiées précédemment : le réservoir bobiné hyperbare de type IV. On
détaillera la construction et les spécificités du modèle EF du réservoir. On fera une étude
de l’éclatement du réservoir et des différents phénomènes qui se manifestent pendant les
196
C,C 4 .
La prédiction de l’éclatement de ces structures est un type d’étude assez récent dans le
domaine des matériaux composites. On trouve des travaux sur l’éclatement des tubes
bobinés [147–151] et quelque travaux sur l´éclatement des réservoirs comme [147] et [152]
même si l’analyse des réservoirs reste assez simple. Des travaux plus développés comme
ceux de Thionnet [153] utilisent des modèles de comportement plus riches et plus
complexes (FE2), mais dont le coût numérique est un frein à la réalisation de calculs sur des
structures de taille réelle. Dans notre cas, on utilisera la version 3D du modèle
d’endommagement pour simuler un réservoir de type IV de 2 litres de capacité et une
architecture composite complexe. Les réservoirs de type IV sont composés de trois parties
et matériaux différents (Figure 141) :
197
C,C 4 .
• Très haute pression d’éclatement due aux facteurs de sécurité imposés (entre 2 et 3,
suivant la norme)
• Le réservoir doit résister à au moins 15000 cycles (de 20 à 875 bar) qui représentent
les cycles de vidange @ remplissage au cours de sa vie utile.
• Le réservoir doit assurer ses performances mécaniques dans un intervalle de
température de @40ºC à 85ºC.
Ces réservoirs sont fabriqués par la procédure d’enroulement filamentaire dans les
installations d’ALHYANCE du CEA Le Ripault par une procédure automatisée. Le liner
polymère sur lequel le composite sera bobiné est élaboré par la procédure de rotomoulage.
Le temps de fabrication du réservoir est compris entre 1 et 3 heures en fonction du drapage
imposé. Le savoir@faire du CEA a permis d’atteindre une dispersion en termes de masse
très faible, d’environ 18 grammes sur une moyenne de 1750 grammes de masse de matériau
composite, ce qui démontre la reproductibilité de la procédure de fabrication des
structures.
[Link] Géométrie
198
C,C 4 .
On va décrire dans la suite ces trois étapes de modélisation jusqu’à obtenir le modèle
final du réservoir bobiné. Comme on a vu préalablement, la première étape de construction
est la génération des embases et du liner. Ces deux parties sont construites avec les
dimensions présentées dans la figure suivante.
Avec ces données géométriques (en particulier le contour extérieur), on peut modéliser
complètement la structure sur laquelle le composite sera bobiné. La deuxième partie de la
modélisation commence à partir de ce point. La première et principale donnée d’entrée est
la structuration dans la virole, où les angles de chaque pli restent constants. Cette
structuration est une combinaison de plis hélicoïdaux (±θ) et circonférentiels (90º). On
présente dans le tableau suivant la structuration introduite dans le logiciel pour générer la
coque composite.
199
C,C 4 .
Le Tableau 8 montre la structuration choisie. Les angles des plis hélicoïdaux sont
croissants (θ1< θ2< θ3 …). Leur valeur n’est pas donnée pour des raisons de confidentialité.
Chaque angle indiqué dans le tableau correspond en fait à un couple ±θ puisque chaque
couche est composée par des fibres orientées dans les deux directions (bobinage allé et
retour). C’est une des particularités de la procédure d’enroulement filamentaire (Figure
144a). L’autre caractéristique des réservoirs bobinés est l’évolution progressive des angles
de chaque pli vers une orientation circonférentielle (90°) sur le dôme en raison de la
procédure même d’enroulement filamentaire.
a) [3] b)
Figure 144. a) Alternance d'orientations dans l'enroulement, b) Paramètres de calcul de
l'évolution des angles.
200
C,C 4 .
En effet, dans le dôme, l’angle d’un pli donné évolue depuis sa valeur initiale dans la
virole θ jusqu’à une valeur de 90º au point de retour où le pli contourne le pôle puis revient
vers la partie virole. Le calcul de l’angle est géré par le plugin d’une façon automatique
selon la formule suivante :
n
R R @ R0
θ ( R ) = sin 0@1
± δ
R R tl @ R 0
Équation 105
où R0 est le rayon au point de retour du pli, Rtl est le rayon au début du dôme et R est le
rayon courant (Figure 144b). Le paramètre δ contrôle la géodésie de l’enroulement. Si on
choisit δ=0, on obtient la formulation classique de Clairaut [154]. Ce paramètre de
glissement nous permettra de caler le point d’arrêt de la couche simulée pour la faire
coïncider avec l’arrêt observé dans le réservoir réel. L’angle de la couche n’est pas la seule
variable qui change le long du dôme. L’épaisseur aussi change quand on se déplace vers la
fin du dôme. Wound Composite calcule la valeur t de l’épaisseur en accord avec l’équation
suivante :
R + 2× BW ×
Rtl @ R0
Équation 106
Deux nouveaux termes apparaissent dans cette formule : ttl et BW. Le premier
paramètre est l’épaisseur de la couche à la fin de la virole (début du dôme) et le deuxième
est la largeur de la bande de composite, qui contrôle l’augmentation de l’épaisseur à la
proximité de l’arrêt.
201
C,C 4 .
Une fois la géométrie réalisée, une validation préalable consiste à comparer les masses
avec le réservoir de référence. On compare les masses des embases métalliques, du liner
polymère et de la coque composite. En effet, c’est la coque composite, et plus précisément
le dôme (par sa complexité géométrique), qui peut induire des différences par rapport à la
masse réelle de composite mesurée dans le réservoir. La comparaison des masses entre le
modèle EF final et le réservoir de référence est présentée dans le graphique suivant.
202
C,C 4 .
Masse (kg)
2
1,5
1
0,5
0
Modèle EF Réservoir réel
Liner 0,513 0,393
Embases 1,293 1,248
Composite 1,81997 1,819
On constate que les résultats des masses du modèle EF sont proches des données
issues du réservoir de référence. On trouve une équivalence quasi – exacte entre la masse
du composite modélisé et la masse réelle. La plus grande différence provient du liner, qui a
été modélisé avec l’épaisseur la plus grande mesurée sur les radiographies. Or, dans le
réservoir réel, l’épaisseur du liner n’est pas constante et présente des fluctuations assez
importantes.
[Link] Maillage
Une fois que la géométrie est créée et validée, on doit générer le maillage du réservoir
pour pouvoir affecter aux éléments finis les propriétés des matériaux correspondants. Cette
étape, qui ne pose aucune difficulté pour les embases et le liner, est plus ardue pour la
coque composite, en particulier pour les dômes. Comme on a expliqué précédemment, les
angles d’orientation du composite varient entre la valeur en virole et 90° au point d’arrêt. Il
est donc nécessaire d’associer une matrice de rigidité à chaque élément caractérisé par son
orientation. Le logiciel WOUND COMPOSITE crée un jeu de propriétés matériau pour
chaque angle θ de l’intervalle [0°,90°], de sorte que la matrice de rigidité associée à cet angle
corresponde à l’assemblage d’une couche à +θ et d’une couche à −θ (sans couplage entre la
flexion et la traction). Il faut remarquer aussi que la matrice de rigidité de chaque élément
est exprimée dans le repère axisymétrique du réservoir pour faciliter le calcul. Ensuite,
chaque élément fini du dôme est repéré dans le repère global et on lui affecte une valeur
d’angle dépendant de l’angle initial (en virole) et de sa position dans le dôme, selon
l’Équation 105 (Figure 148).
203
C,C 4 .
Figure 148. Exemple de la variation des propriétés dans une structure bobinée, en fonction de la
variation d’angle dans le dôme.
Dans cette figure, chaque couleur représente un angle différent. On peut voir, au début
du dôme, que la variation d’angle est très lente (plusieurs éléments avec la même couleur) ;
puis cette variation accélère jusqu’à avoir une valeur d’angle par élément fini. Le résultat de
cette assignation des propriétés est montré ci@dessous :
On présente ici les propriétés des matériaux utilisés dans le réservoir, à savoir : le liner
en polymère, les embases métalliques et la coque composite. Les propriétés du liner en PA6
204
C,C 4 .
PA6
E (MPa) ρ (g/cm3) υ
> 0@ 0
Tableau 10. Proprietes elastiques et masse volumique du polymere pa6 utilise dans le liner.
Dans le cas des embases, le matériau utilisé est un acier inoxydable 316L traité pour
améliorer sa résistance à la fissuration par hydrogène. On a choisi un comportement
élastoplastique suivant une loi bilinéaire. Les caractéristiques mécaniques élastiques sont
présentées dans le Tableau 11.
A316L
E (MPa) ρ (g/cm3) υ
0
Tableau 11. Propriétés élastiques et masse volumique de l'acier 316L (embases).
La loi plastique a été identifiée à partir d’un essai de traction fait par un de nos
partenaires. La courbe de traction à rupture et la loi identifiée sont montrées dans la Figure
150. Comme les déformations de l’embase resteront faibles, on identifie le comportement
plastique sur le seul début de la courbe de traction (Figure 150b).
600 400
Stress (Mpa)
500 350
300
400 250
300 200 Test
Test 150
200 Loi identifiée
100
100 50
0 0
0 0,2 0,4 0,6 0 0,02 0,04
Strain Strain
a) b)
Figure 150. a) courbe de traction a rupture de l'acier 316l et b) identification de la loi plastique.
E1 (MPa) E2 (MPa) E3 (MPa) υ12 υ23 G12/13 (Mpa) G23 (Mpa) ρ (g/cm3)
134400 8180 8180 0.36 0.38 4700 2920 1.55
Tableau 12. Propriétés mécaniques du composite bobiné dans un repère local.
Il faut remarquer que ces propriétés sont celles correspondant au repère local du
matériau en prenant la fibre comme direction 1. Dans les réservoirs, tous les calculs sont
faits dans le repère axisymétrique. Il est donc nécessaire de procéder à des changements de
bases pour pouvoir intégrer le modèle d’endommagement dans le calcul du réservoir. Le
schéma utilisé est le suivant :
De cette façon, le modèle fonctionne toujours dans le repère local mais on profite des
orientations axisymétriques générées par WOUND COMPOSITE MODELER.
Pour reproduire les conditions expérimentales de l’essai d’éclatement ainsi que les
mesures enregistrées, on a bloqué la totalité des degrés de liberté sur la surface à l’extrémité
d’une des deux embases et laissé l’autre extrémité libre. Le déplacement axial de l’embase
libre sous l’effet du chargement sera comparé aux mesures, ainsi que le déplacement radial
du centre de la virole. Un schéma récapitulatif est montré sur la Figure 151.
Le chargement appliqué est une rampe de pression exercée sur toute la surface interne
entre 0 bar et 2000 bar (Figure 141). Avec les conditions aux limites et le chargement, on a
206
C,C 4 .
Eclatement du réservoir
4
Déplacement Axial
3,5
Déplacement (mm)
(expérimental)
3
Déplacement Axial
2,5
(simulation)
2
Déplacement Radial
1,5
(expérimental)
1
Déplacement Radial
0,5
(simulation)
0
0 500 1000 1500
Pression (bar)
On observe une forte non linéarité dans le début du déplacement axial expérimental.
L’origine et la cause de cette non linéarité seront étudiées plus tard. Pour mieux comparer
les rigidités axiale et radiale, on a corrigé ce début non linéaire en obtenant la courbe
suivante.
207
C,C 4 .
Eclatement du réservoir
4
Déplacement Axial
3,5
Déplacement (mm)
(expérimental)
3
Déplacement Axial
2,5
(simulation)
2
Déplacement Radial
1,5
(expérimental)
1
Déplacement Radial
0,5
(simulation)
0
0 500 1000 1500
Pression (bar)
Figure 153. Courbes de déplacement expérimentale et simulée du réservoir de type IV à 20ºC
avec le déplacement axial corrigé.
Le point d’éclatement simulé se situe à une pression de 1600 bar, alors que la pression
d’éclatement expérimentale est de 1750 bar, ce qui nous donne une erreur d’environ @8%.
On peut observer sur le déplacement radial expérimental une perturbation qui commence
autour de 1500 bar. Cette perturbation est due aux premières ruptures de fibres dans des
couches circonférentielles. Ces ruptures ne provoquent pas la rupture définitive, dans la
mesure où les couches encore intactes reprennent les contraintes. Ce n’est que lorsqu’un
certain nombre de couches circonférentielles sont rompues qu’on assiste à l’explosion du
réservoir. Une autre information peut être déduite en analysant la forme de la courbe du
déplacement simulé : le mode d’éclatement. Il y a deux types d’éclatement :
• « sain » si les embases ne sont pas éjectées. Ce mode se traduit par une chute
brusque du déplacement axial sur la courbe simulée.
• « non sain » si les embases sont éjectées de la coque composite. Dans ce cas, le
déplacement axial est toujours croissant.
Les éclatements du type « sain » sont caractérisés par une rupture le long de l’axe du
réservoir et une ouverture des plis dans la direction circonférentielle. On verra dans la suite
des calculs que ce type d’éclatement est provoqué par une rupture des couches
circonférentielles dans le centre de la virole. Par contre, l’éclatement « non sain » est piloté
par des ruptures de fibres et du délaminage dans le dôme qui provoquent l’éjection des
embases. On voit un exemple de chaque rupture sur la Figure 154.
208
C,C 4 .
a) b)
Figure 154. a) Réservoir éclaté en mode sain et b) non sain.
Avant l’éclatement, on peut constater que les rigidités axiale et radiale sont assez bien
reproduites par la simulation. On retrouve également, par la simulation, la légère non
linéarité de la courbe du déplacement axial à partir de la moitié du chargement. Cette non
linéarité est provoquée par deux phénomènes différents :
Bouchon
On voit sur cette figure que la plasticité a déjà évolué à 600 bar (avec un maximum de
l’ordre de 2.10@3) dans trois endroits distincts : à proximité du bouchon scellant l’embase,
dans le col de l’embase et au niveau de la jonction liner@embase@composite. Ces trois points
constituent des zones assez critiques, notamment pour l’embase. Cette plasticité est très
209
C,C 4 .
influencée par la forme du dôme composite, ce qui nous indique l’importance de modéliser
proprement cette zone du réservoir bobiné.
La deuxième non linéarité, située au début de la courbe de déplacement axial, n’est pas
expliquée a priori par la simulation, qui prévoit un comportement initial linéaire. Une piste
possible, observée sur certaines radiographies de haute résolution avant éclatement (Figure
156), est la présence initiale d’un jeu entre l’embase/liner et la coque composite qui apparaît
lors du refroidissement du réservoir après le cycle de cuisson. Afin de vérifier cette
hypothèse par la simulation, la longueur de la virole du liner a été diminuée de 0.25mm. On
a donc ménagé un espace entre composite et embase/liner de cette valeur et on a ajouté un
contact entre les surfaces en regard. Le résultat obtenu est montré sur la Figure 157 ci –
dessous.
Déplacement Axial
3,5
(expérimental)
3
Déplacement Axial
2,5
(simulation)
2
Déplacement Radial
1,5
(expérimental)
1
Déplacement Radial
0,5
(simulation)
0
0 500 1000 1500
Pression (bar)
Figure 157. Simulation avec un jeu entre le composite et l'embase dans les dômes (le
comportement du composite est considéré ici comme élastique sans endommagement)
On observe dans cette figure que l’introduction de ce jeu nous permet de simuler
correctement la première partie de la courbe axiale.
210
C,C 4 .
1580 bar
1600 bar
1620 bar
Figure 158. Endommagement des fibres (une valeur de 1 représente des fibres rompues) (1580
bar, 1600 bar et 1620 bar) .
On voit dans cette séquence que la première couche qui atteint la rupture des fibres est
la couche circonférentielle la plus interne, c’est@à@dire la couche la plus chargée. Très
rapidement, les autres couches circonférentielles commencent à rompre jusqu’à atteindre la
plus externe. Le réservoir n’est alors plus capable de supporter la pression et le
déplacement radial augmente d’une façon divergente.
211
C,C 4 .
L’endommagement matriciel n’a pas une conséquence aussi brutale que la rupture des
fibres, mais on peut aisément suivre sa trace, notamment par l’analyse de l’émission
acoustique. La figure suivante montre une graphique amplitude d’événements acoustiques –
temps enregistré par le partenaire CEA sur l’ensemble du réservoir pendant l’essai
d’éclatement :
Figure 159. Emission acoustique pendant l’essai d'éclatement, (temps en s en abscisse, 2000
seconds = 650 bar).
Les évènements ont été classifiés dans deux groupes différents : le groupe bleu
correspondant aux évènements très énergétiques (typiquement la rupture des fibres) et le
groupe rouge lié plutôt aux ruptures matricielles / délaminage. On observe que c’est à
partir de 1500 secondes ( 500 bar) que le groupe rouge compte un grand nombre
d’événements. On peut utiliser la simulation pour faire le lien entre cette activité acoustique
et l’état d’endommagement matriciel, diffus et le délaminage.
En effet, on constate dans la Figure 160 que pour une pression de 600 bar,
l’endommagement matriciel et hors – plan commence à se développer, ce qui pourrait
expliquer l’augmentation significative des échos acoustiques. Il est intéressant d’observer
212
C,C 4 .
Pression
Endommagement matriciel Endommagement hors 0 plan
(bar)
1000
1200
1400
1600
Rupture
Tableau 13. Evolution des endommagements matriciel et hors @ plan au cours du chargement.
213
C,C 4 .
a) b)
Figure 161. a) Détail de l’endommagement hors @ plan et b) fuite détectée lors des essais.
Il est intéressant de remarquer que l’endommagement hors – plan génère trois voies
possibles de fuite parmi lesquelles deux ont été observés expérimentalement (b et c) : en
cas de défaillance de la jonction embase liner (matérialisée par le trait noir, Figure 161a), le
fluide servant à pressuriser le réservoir (en l’occurrence, de l’eau) peut soit traverser le
composite au droit de la jonction (a), soit diffuser dans l’extrémité du dôme (b), soit suivre
un chemin le long de l’interface embase / composite. On montre sur la Figure 161b un
exemple de fuite détectée lors des essais, semblant correspondre au chemin (c) (fuite en
extrémité de dôme, à l’interface avec l’embase). Ce type de simulations montre l’intérêt
d’implémenter des modèles complexes dans les calculs car eux seuls sont capables de
fournir une description fine de l’état de l’endommagement du réservoir.
Les simulations qui précèdent ont montré leur intérêt pour interpréter les profils
d’émission acoustique ainsi que les probables chemins de fuite. Mais alors que
l’endommagement des fibres avait une influence prépondérante sur le comportement et la
rupture globale du réservoir, on peut se demander si les endommagements matriciel et hors
plan ont un effet mesurable sur le comportement macroscopique du réservoir (en
particulier les courbes de déplacement). En ce qui concerne le délaminage, on n’observe
pas d’effet notable sur les courbes de déplacement. Ce n’est par contre pas le cas de la
rupture matricielle. En effet, un état avancé d’endommagement matriciel peut provoquer
une surcharge des fibres des plis circonférentiels car la charge supportée par les plis
endommagés est redistribuée entre les plis « sains ». Cette redistribution provoque une
activation prématurée des ruptures des fibres et donc des inflexions sur les courbes de
déplacement. Afin de vérifier et étudier qualitativement cette influence de
l’endommagement matriciel, on a mis en place deux calculs du même réservoir : un premier
214
C,C 4 .
dans lequel l’endommagement matriciel est activé avec un effet sur la rigidité plus sévère
que dans le cas standard et un second avec le seul endommagement des fibres. La
comparaison est montrée sur la Figure 162 ci@dessous.
2,5
Déplacement radial (Exp)
2
Déplacement axiale (sans
endommagement)
1,5
Déplacement radial (sans
endommagement)
1
Déplacement axial (avec
0,5 endommagement)
Déplacement radial (avec
0 endommagement)
0 500 1000 1500
Pression (bar)
Figure 162. Influence de la rupture matricielle dans le comportement du réservoir.
On constate une augmentation de presque 200 bar entre les deux simulations (avec et
sans endommagement matriciel). Si nous regardons aussi les courbes de déplacement, on
peut noter une petite différence entre les deux courbes de déplacement radial (ellipse en
pointillé) : la souplesse est augmentée lorsqu’on prend en compte l’endommagement
matriciel. Avec ce test, on montre l’importance de l’étude de la fissuration matricielle et de
l’effet de cet endommagement sur la redistribution des contraintes pour pouvoir prédire
correctement la pression d’éclatement.
Quelques réservoirs bobinés ont été instrumentés avec plusieurs fibres optiques
placées dans différentes couches pour mesurer la déformation dans le sens de la fibre, c’est@
à@dire dans la direction circonférentielle, les fibres optiques ayant une orientation à 90°. De
cette façon, on dispose d’un autre jeu de données (très locales) avec lequel comparer le
résultat de la simulation.
215
C,C 4 .
Figure 163. Comparaison des déformations de certaines couches (en rouge, à droite, couches
circonférentielles)
On trouve aussi un très bon accord entre les courbes expérimentales et simulées,
notamment dans les couches C12 et C17. On parvient à capter la non linéarité des
déformations à partir de 800 bar. Cette légère perte de linéarité est provoquée par le
développement de l’endommagement matriciel dans les couches hélicoïdales intérieures.
4
Simulation probabiliste
Déplacement axial (exp)
Déplacement (mm)
3,5
3 Déplacement radial (exp)
2,5
Déplacement axial (déterministe)
2
Déplacement radial (déterministe)
1,5
1 Déplacement axial (Probabiliste)
0,5
Déplacement radial (Probabiliste)
0
0 500 1000 1500
Pression (bar)
Figure 164. Comparaison entre simulations probabiliste et deterministe de l'éclatement du
réservoir.
216
C,C 4 .
On constate dans ces courbes l’effet 3D dans la redistribution des contraintes car la
divergence de la courbe est plus progressive. Dans ce cas, la pression d’éclatement évolue
de 1720 à 1750 bar, à comparer aux 1756 bar du réservoir de référence. On peut regarder
l’évolution des ruptures des fibres dans la première couche circonférentielle (la plus proche
217
C,C 4 .
1400bar 1600bar
1800bar
Figure 166. Endommagement des fibres de la première couche circonférentielle @simulation
probabiliste 3D (noir : fibre rompue).
On note qu’à 1400 bar la première couche commence à présenter des éléments
rompus, localisés en différents endroits sans coalescence notable. Au cours du chargement,
les zones endommagées finissent par se rejoindre jusqu’à arriver à un état complètement
dégradé où la couche ne résiste plus, d’autant plus que les couches supérieures sont
également fortement endommagées.
Cette section a démontré l’intérêt d’une démarche probabiliste dans le calcul des
structures composites, et du réservoir spécifiquement : non seulement elle permet de
simuler une pression d’éclatement supérieure (et plus conforme à la réalité) à celle obtenue
par une approche purement déterministe, mais elle donne accès (dans sa version 3D) au
scénario de rupture progressive des fibres dans les différentes couches. Cette approche
pourrait désormais être améliorée en affectant ce caractère probabiliste à d’autres propriétés
du matériau, telles que la rigidité.
218
C,C 4 .
5.4 Conclusions
Une fois le modèle validé, on l’a appliqué au calcul de l’éclatement du réservoir bobiné
type IV. Le modèle a montré des résultats satisfaisants, même dans sa version déterministe
(erreur d’environ 7% sur la pression d’éclatement). On a démontré l’intérêt d’utiliser des
modèles d’endommagement riches pour prédire des phénomènes complexes comme la
formation des chemins de fuite. D’après l’étude réalisée, on peut constater que la
modélisation axisymétrique du réservoir est suffisante pour décrire son comportement
même si l’effet probabiliste est un peu moins puissant dans cette modélisation.
219
CHAPITRE VI:
PERSPECTIVES
C( , 0
6.1 Introduction
222
C( , 0
a) b)
Figure 167. a) Cellule élémentaire, b) Chargements élémentaires pour calculer l'influence sur la
rigidité (Huchette 2005).
De cette façon, on pourrait calculer d’une façon beaucoup plus fine la valeur des
coefficients , et et éventuellement enrichir le potentiel de termes supplémentaires
puisés dans la base d’intégrité. Cette étude pourrait être étendue pour analyser d’autres
phénomènes et influences comme le couplage existant entre le délaminage et la fissuration
matricielle (endommagement hors – plan et matriciel dans notre modèle). Ce couplage
entre ces deux endommagements a été introduit par Laurin et al [29] dans le modèle de
l’ONERA pour les composites aéronautiques. Ils ont constaté expérimentalement dans des
essais de flexion quatre points sur des cornières que la présence de plis fissurés provoquait
une accélération non négligeable dans l’apparition du délaminage (et donc de la ruine de la
structure).
223
C( , 0
étude, on a utilisé une éprouvette unidirectionnelle (fibres parallèles à la longueur) avec les
dimensions suivantes (la longueur de 83,3 mm représente la distance entre les appuis) :
La simulation de l’essai de flexion trois points représenté sur la Figure 169 fournit la
courbe force en fonction du déplacement du rouleau représentée sur la Figure 170 :
224
C( , 0
Malgré une base d’identification réduite (essais de traction dans le sens des fibres et
dans la direction perpendiculaire aux fibres), on constate un bon accord entre la courbe
expérimentale et la simulation ne prenant en compte que l’endommagement hors – plan.
En fait, la simulation montre les quatre étapes observées dans l’essai réel :
225
C( , 0
# -2
Dans les années 70, [158] ont présenté un modèle de dégradation des propriétés
résiduelles linéaire assez simple dans lequel la résistance du matériau décroît linéairement en
fonction du nombre de cycles appliqués. Puis, d’autres auteurs ont introduit le concept de
taux de changement de la résistance résiduelle [159–161]. Ils ont supposé qu’il existe une
seule relation entre la résistance statique et cyclique. Cette relation implique que si un
spécimen est associé à un certain rang dans sa distribution de résistance statique, son
équivalent cyclique doit correspondre au même rang de la distribution de durée de vie.
Cette idée a été appelée plus tard SLERA (strength – life equal rank assumption). Sur le
même principe, Yang et al [162–168] ont présenté des modèles pour prédire la résistance
résiduelle à différents niveaux de charge et d’amplitude. Des courbes S – N basées sur
différents modèles de dégradation ont été présentés par [169,170].
Le modèle « d’interaction » a, quant à lui, été introduit par Adam et al [171]. Ce modèle
est purement phénoménologique et combine, dans un critère sous forme de loi puissance,
une résistance résiduelle normalisée et la durée de vie correspondante.
Des modèles adaptés aux composites qui subissent des chargements complexes ont été
aussi développés, comme le « modèle généralisé de dégradation des propriétés résiduelles »
[172]. Une grand partie de ces modèles (Broutman and Sahu, Hahn and Kim, Yang et al,
modèle d’interaction et celui de Sendeckyj) ont été comparés dans le papier [173].
226
C( , 0
# - !
On fournit ici quelques éléments concernant les modèles construits dans le cadre de la
mécanique de l’endommagement continu. Ce type de modèle essaie de reproduire la
dégradation progressive du matériau au cours du chargement. Les premiers modèles ont été
présentés par Gao [174]. Puis de nombreuses contributions ont apparu comme celles
développées par Hochard et al [175–177]. Ce dernier modèle est une adaptation du méso –
modèle de Cachan (décrit dans la section 2.3.2) pour des chargements cycliques. La
première hypothèse introduite pour passer du modèle statique vers le modèle de fatigue
suppose que les variables d’endommagement ont deux composantes différentes, une
statique et une autre cyclique. Dans le modèle de Hochard, la seule variable
d’endommagement qui subit une dégradation dans le modèle de fatigue est celle associé à
l’endommagement diffus .
suivante :
Équation 107
d’un cycle. correspond au seuil de fatigue. La Figure 172 et la Figure 173 montrent
227
C( , 0
Équation 108
On aura donc un terme du potentiel construit avec des variables dites statiques
les effets de fatigue. On notera que l’endommagement de fatigue agit de la même manière
228
C( , 0
sur la rigidité que l’endommagement statique. Seule l’évolution diffère. Certaines conditions
sont à ajouter comme :
Équation 109
Équation 110
Avec
Équation 111
229
C( , 0
Équation 112
230
C( , 0
Figure 175. Influence du paramètre n dans l'évolution des variables d'endommagement (d).
Le dernier détail qui nous reste à régler pour implémenter cette extension est l’effet du
cyclage sur la tenue résiduelle. Le changement de tenue résiduelle est traduit par une
diminution du seuil d’endommagement statique qui devra être fonction de la variable
d’endommagement cyclique. Pour prendre en compte cette influence de la variable cyclique
sur le seuil statique, on propose de modifier les expressions des seuils d’initiation
d’endommagement et de rupture, par exemple de la manière suivante :
Équation 113
Équation 114
On trouve dans la littérature des travaux sur des éprouvettes unidirectionnelles comme
[178] qui testent des UD à 0º. Afin de conclure d’une manière très qualitative sur les
capacités de l’approche proposée à traduire la réduction de résistance suite à un chargement
cyclique, on montre une courbe d’évolution de la contrainte à rupture en fonction du
231
C( , 0
nombre de cycles issue de ces travaux et le même type de courbe générée à partir du
modèle dans la Figure 176. On peut voir que les deux courbes se ressemblent beaucoup
avec un long intervalle sans influence du cyclage et une chute très brusque à partir d’un
seuil déterminé.
On a programmé cette extension en fatigue dans une UMAT et testé sur un unique
élément pour vérifier le comportement des variables de fatigue pour plusieurs chargements
différents. A titre d’illustration, la variable affectée sera celle associée à la rupture des fibres.
Ces chargements (traction) sont montrés ci – dessous :
Chargement
Cahrge (Mpa)
400
Chargement 1
300
Chargement 2
200
Chargement 3
100
0
0 1 2 3 4 5 6 7 8 9 10
Temps (s)
On constate que ces chargements ont des charges maximales différentes (1 et 2) ou des
amplitudes différentes (1, 2, et 3). La Figure 178 montre que ces différences influencent
l’évolution de la variable d’endommagement. Il faut garder en tête que ses simulations sont
purement qualitatives. En particulier, aucun paramètre n’a été identifié et que leurs valeurs
232
C( , 0
ont été choisies de telle sorte que l’endommagement s’active pour un faible nombre de
cycles.
Evolution de
l'endommagement
1
Endommagement
Chargement 1
0,8 Chargement 2
0,6 Chargement 3
0,4
0,2
0
0 50 100 150 200
Temps (s)
Pour finir avec cette étude du modèle cyclique, on analyse l’effet de l’endommagement
cyclique sur la tenue résiduelle. On montre ci@dessous une simulation d’une traction après
cyclage.
233
C( , 0
Cyclage 0 traction
Contrainte (Mpa)
1000
Cyclage Traction
800
600
400
200
0
0 0,005 0,01 0,015 0,02
Déformation
Il faut noter que dans ce cas, avec les valeurs de paramètres utilisées, la rupture statique
se produit à 2104 MPa. Après cyclage, nous trouvons une rupture autour de 1000 MPa.
Pour comprendre ce phénomène, on rappelle que le modèle est formulé en 5 0
Dans le cas statique, on identifie la valeur de la force thermodynamique, donc de la
déformation, qui correspond à la charge de rupture. Le cyclage ayant dégradé le matériau et
donc sa rigidité, la contrainte correspondant à la déformation à rupture en statique est
désormais beaucoup plus faible. Il sera nécessaire d’étudier expérimentalement cette
influence entre la perte de rigidité et la perte de tenue résiduelle pour pouvoir valider cette
influence ou introduire des modifications supplémentaires. Il reste encore beaucoup de
travail sur cette proposition de modèle incrémental pour être capable de reproduire
correctement les phénomènes observés dans les essais d’identification et les structures
simples. Néanmoins, on a implémenté le modèle tel qu’il a été développé dans une UMAT
pour l’appliquer à une structure. On a choisi de simuler le comportement des éprouvettes
entaillées, de séquence quasi – isotrope. On montre ci@dessous l’endommagement matriciel
dans le pli transverse à différents cycles.
234
C( , 0
On peut voir ici que l’endommagement matriciel est très développé dans le centre de
l’éprouvette et qu’il continue à évoluer vers l’extérieur en raison de la croissance du cumul
de force thermodynamique. On peut remarquer aussi que l’endommagement cyclique
n’atteint pas la valeur 1 sur le bord libre des entailles. Ce fait est une conséquence de
l’interaction entre l’endommagement statique et cyclique. En fait, l’endommagement
statique évolue lors de la mise en charge de l’éprouvette (avant l’application des cycles), ce
qui provoque un arrêt de la variable cyclique à une valeur plus petite que l’unité. On
rappelle que la somme des variables cyclique et statique doit être toujours comprise entre 0
et 1. Cette configuration d’endommagement matriciel a été retrouvée sur des observations
(micrographies des entailles) d’éprouvettes entaillées soumises à des chargements cycliques.
235
C( , 0
a b
236
C( , 0
On peut voir que les sections créées, associées chacune d’elles à une orientation
donnée, sont construites de telle sorte que le maillage soit propre, en améliorant la qualité
du modèle EF.
237
CONCLUSION
L’étude présentée dans ce manuscrit a eu pour objectif le développement d’un outil de
simulation du comportement mécanique des composites bobinés pour prédire avec
précision la pression d’éclatement des réservoirs bobinés de type IV. Après avoir présenté
quelques généralités sur les matériaux composites et s’être focalisé sur les aspects les plus
spécifiques des composites bobinés comme par exemple l’architecture entrecroisée, un
rapide tour d’horizon des modèles d’endommagement existant dans la littérature a montré
que la plupart des modèles ont été conçus pour étudier les composites stratifiés (de type
aéronautique) et appliqués sur ce type de composite. Les spécificités des matériaux
composites bobinés et de leur application aux réservoirs hyperbares constituent une limite à
ces approches, du moins dans leurs versions les plus communément utilisées : en effet, les
essais réalisés ont pu mettre en évidence une très forte non linéarité du comportement en
cisaillement ainsi qu’une viscosité non négligeable particulièrement visible lors des
décharges. Les modes d’endommagement associés aux profils épais des réservoirs sont
également un frein à l’utilisation directe des modèles construits pour les stratifiés minces.
On notera cependant que ces particularités de comportement sont loin d’être des obstacles
infranchissables pour des modèles tels que le méso@modèle de Cachan ou le modèle de
l’ONERA : des modifications sans doute assez indolores (mais qui restent à effectuer…)
leur permettraient de simuler correctement le comportement des réservoirs. Ce travail de
thèse n’a pas envisagé de procéder à ces modifications, mais de développer une approche
de modélisation alternative, résultant de la mise en synergie d’outils élémentaires que sont la
thermodynamique des processus irréversibles, l’approche à directions fixes et la théorie de
représentation de la fonction tensorielle. Le but du travail de thèse présenté ici a été de
démontrer que l’association de ces outils était à même de fournir un modèle dédié au
comportement des composites bobinés. La stratégie a été de prouver le caractère
suffisamment général de ces ingrédients, capables de prévoir la dégradation du composite
sans préjuger de la forme de la perte de rigidité.
241
perspectives pour une utilisation dans tout type de matériaux composites comme les
stratifiés, les CMC, les tissus 2D et 3D… c’est@à@dire les matériaux architecturés. Le
potentiel thermodynamique (le cœur du modèle) a la forme d’une énergie libre en
déformation. Ce fait nous permet d’implémenter le modèle dans les codes de calcul par
éléments finis industriels aisément en contournant l’utilisation de méthodes de résolution
d’équations non linéaires du type Newton Raphson (nécessaires quand le modèle est
formulé en contraintes) et donc épargner du temps de calcul. Le modèle ne traite pas que
de l’endommagement mais la déformation résiduelle et la viscosité de la matrice sont aussi
incorporées au potentiel. La déformation résiduelle a la même nature que les variables
d’endommagement : l’évolution de ces variables est pilotée par des forces
thermodynamiques et suit la règle de normalité. En ce qui concerne la viscosité, un modèle
très simple a été introduit avec le but de montrer qu’on pouvait coupler viscosité,
déformation résiduelle et endommagement. La variable visqueuse, qui a la forme d’une
contrainte, évolue selon un écart par rapport à l’état d’équilibre. Contrairement à
l’endommagement, ce mécanisme ne présente pas de seuil. Il reste cependant encore de
travail à faire dans cette partie du modèle comme l’introduction de plusieurs temps de
relaxation, des fonctions non linéaires, etc… On notera que ce modèle visco –
élastoplastique – endommageable est capable de reproduire le comportement d’éprouvettes
bobinées à ±45º (qui diffère de celui de leurs homologues stratifiées) avec leur forte non
linéarité, les fortes déformations résiduelles et l’apparition de boucles d’hystérésis.
Jusqu’à ici, le modèle développé reste déterministe, comme la plupart des modèles
présentés dans la littérature. On a essayé d’introduire la nature stochastique du
comportement des fibres pour mieux représenter le vrai comportement du composite.
Avec une loi de mélange classique, on a introduit les deux sources principales de
variabilité : la rupture de la fibre seule et la variabilité du taux volumique des fibres. Pour
cette dernière source, on a développé une méthode d’identification du taux volumique par
traitement d’images, qui permet de lier le comportement probabiliste à la taille des éléments
finis d’une façon approximée. Dans le but d’avoir plusieurs données différentes pour
comparer les résultats fournis par le modèle avec les résultats expérimentaux, on a ajouté
une signature acoustique virtuelle aux modes d’endommagement pour pouvoir comparer le
cumul de ces signatures avec celles enregistrées. Cette approche originale ouvre une
nouvelle voie dans l’étude et la compréhension des émissions acoustiques et pourrait être
utilisé dans l’identification des cinétiques des différents modes d’endommagement.
242
Finalement, le processus d’identification est décrit en montrant que quatre essais de
caractérisation du matériau (dans le sens des fibres, perpendiculaire aux fibres et sur
séquence ±45°, ainsi qu’une combinaison 0°/90°) suffisent à identifier l’ensemble des
paramètres (sous couvert de certaines hypothèses).
Ce modèle est validé par comparaison avec des essais de traction sur des éprouvettes
entaillées bobinées. Dans le cadre cette thèse, une importante base de données a été créée
avec l’objectif de disposer des données suffisantes pour bien comprendre le comportement
du matériau bobiné, la cinétique et l’interaction des modes d’endommagements. Les
éprouvettes testées ont une architecture purement bobinée, à l’image du matériau présent
dans le réservoir. Une procédure originale de fabrication a été mise en place par un de nos
partenaires. On a ajouté des entailles pour concentrer l’endommagement et les faire
apparaître simultanément de façon que l’état des contraintes du matériau soit suffisamment
complexe et constitue un cas d’étude discriminant pour le modèle d’endommagement. On
a montré que la haute température met en évidence l’importance de la viscosité dans le
comportement du composite bobiné et aussi l’apparition prématurée et étendue du
délaminage. Quelques pistes sur la cause de cette apparition ont été présentées, comme
l’influence des ondulations créées par les plis croisés et qui favorisent l’amorçage et la
propagation des ruptures hors – plan. Toutes les éprouvettes ont été contrôlées par
émission acoustique, micrographies et radiographies, ce qui nous a donné accès à un
ensemble de données permettant de mieux comprendre quel type d’endommagement
apparaît, à quel endroit et à quel moment de l’essai.
243
des données très macroscopiques comme les courbes force @ déplacement. Les résultats
montrent un très bon accord, que ce soit la version coque ou la version complète 3D. Mais
la richesse du modèle complet nous permet d’aller vers une échelle beaucoup plus fine et
prédire les types d’endommagement présents ainsi que leur localisation dans l’éprouvette,
dans le plan ou dans l’épaisseur. Les comparaisons avec des vue latérales des éprouvettes,
notamment le fond des entailles, montrent les performances du modèle pour arriver à
prédire des phénomènes comme la fissuration non centrée du pli transverse dans le cas
quasi – isotrope, le délaminage prédominant entre les plis 0º/±θ et 90º/±θ et le début du
phénomène de splitting. Grâce aux radiographies, on a pu comparer l’extension des
endommagements dans le plan du pli. Les prédictions restent aussi très proches des
observations expérimentales dans ce cas. On a constaté que la fissuration matricielle est très
bien prédite même loin des entailles où les radiographies montrent des fissures, notamment
sur le bord libre. Cette fissuration dans le bord libre est captée par le modèle qui montre
aussi un délaminage entre la couche transverse et les couches croisés adjacents. Ce fait nous
indique la présence d’un couplage indirect entre l’endommagement hors – plan et
l’endommagement matriciel. La fissuration matricielle n’est pas le seul endommagement qui
est reproduit d’une façon satisfaisante, l’endommagement hors – plan est aussi bien capté
au voisinage de l’entaille même si son étendue est bien plus large dans les éprouvettes
réelles.
Les comparaisons menées sur les éprouvettes entaillées permettent de valider les
performances du modèle pour reproduire l’endommagement du composite bobiné ainsi
que la rupture finale de ces structures contenant des concentrateurs de contraintes et donc
un état de contraintes très complexe. Cette validation nous garantit un certain degré de
confiance par rapport au modèle pour pouvoir simuler des structures encore plus
complexes comme le réservoir bobiné de type IV.
Le dernier test pour montrer la capacité du modèle est la simulation d’une vraie
structure industrielle : le réservoir bobiné. Deux modèles éléments finis ont été créés : un
axisymétrique et un autre 3D. Ces structures montrent un degré de difficulté assez
important en raison de sa nature multi matériaux (embases métalliques, liner polymère et
coque composite) et de l’architecture bobinée du matériau composite qui génère des
géométries et des orientations complexes au niveau du dôme. On a utilisé le plugin Abaqus
Wound Composite Modeler pour pouvoir prendre en compte cette géométrie et ces
orientations du pli variables dans les deux dômes. L’analyse de l’éclatement du réservoir
244
modélisé montre que le modèle déterministe arrive à bien reproduire le comportement du
réservoir en restant un peu conservatif dans la prédiction de l’éclatement (7% d’écart). Cet
éclatement est piloté complètement par la rupture des fibres mais est influencé par d’autres
endommagements, principalement la fissuration matricielle. Le délaminage n’a pas une
forte influence sur l’éclatement. Par contre, il peut provoquer une voie de fuite qui part de
la jonction embase – liner – composite et arrive jusqu’à l’extérieur. Dans le cas d’une
fracture du liner à cet endroit – là, l’hydrogène trouverait une voie pour s’échapper. Ce
scénario a été confirmé expérimentalement, plusieurs réservoirs ont montré des fuites avec
une grande quantité de liquide d’essai qui s’échappait par les endroits prédits par le modèle.
245
ANNEXES
' !
a) b)
Figure 182. Modes d'éclatement a) sain et b) non sain.
Il est très important de vérifier si le modèle est capable de prendre en compte ces deux
types d’éclatement, notamment pour son utilisation dans des démarches d’optimisation, où
l’éclatement sain est une de contraintes à respecter. On a créé deux architectures de
réservoirs différentes de celles utilisées dans le chapitre 5 pour étudier les capacités du
modèle à simuler ces deux éclatements et pour comprendre quels mécanismes ont le plus
d’influence dans une mode ou un autre. Pour mesurer l’épaisseur, on considère
249
' !
Modèle 1 Modèle 2
90 90
α2 α8
α3 90
90 α9
α4 90
60 α10
α5 90
α6 α11
90 90
La forme du dôme correspondant à ces deux architectures est montrée sur la figure
suivante:
a) b)
On peut observer dans cette figure que le modèle 1 présente un dôme beaucoup plus
fin que dans le cas du modèle 2. En regardant les deux architectures, il est clair qu’une
présence importante de couches avec un angle inférieure à 30º génère un dôme plus épais
que dans le cas où les couches sont plus orientées vers la direction circonférentielle.
4
3,5
Déplacement radial
3 modèle 1
2,5
2 Déplacement axial
1,5 modèle 2
1
0,5 Déplacement radial
modèle 2
0
0 500 1000 1500
Pression (bar)
En regardant cette figure, on constate que le modèle est capable de différencier les
deux types d’éclatement : sain pour le modèle 2 et non sain pour le modèle 1. Le modèle 1
montre un déplacement axial qui diverge, indiquant une éjection des dômes. En plus, on
observe que le déplacement radial reste linéaire sans montrer aucune fluctuation (pas de
rupture dans la virole). Le modèle 2 montre un comportement identique à celui analysé
dans le chapitre 5.
En analysant la rupture des fibres des deux réservoirs, on obtient le résultat suivant:
a) Modèle 1 b) Modèle 2
Figure 185. Rupture des fibres (rouge=totalement endommagé) pour a) modèle 1 et b) modèle 2.
On observe que le modèle 1 (rupture non saine) montre une rupture de fibres localisée
à l’extrémité du dôme et non pas dans la virole, comme pour le modèle 2. C’est ce
changement de localisation de la rupture des fibres qui fait basculer le type d’éclatement de
251
' !
sain vers le non sain. Si nous observons l’endommagement matriciel juste avant la rupture
des fibres, nous constatons un niveau nettement supérieur dans le modèle 1 (Figure 186).
a) Modèle 1 b) Modèle 2
Il est clair que dans le modèle 1, l’endommagement matriciel est concentré à proximité
des arrêts des plis, précisément où les fibres vont commencer à rompre. Cette rupture
matricielle provoquera une surcharge qui va favoriser encore plus la rupture des fibres dans
cette zone et donc, l’éclatement non sain.
Avec cette étude, on montre l’intérêt d’utiliser un modèle d’endommagement fin pour
pouvoir discerner ces deux modes d’éclatement et donc, pouvoir améliorer des futures
démarches d’optimisation des réservoirs bobinés.
252
BIBLIOGRAPHIE
[1] Laurin F. Introduction générale sur les matériaux composites. Colloque Aussois
MECAMAT, Aussois: 2011.
[5] Binienda WK, Wang Y. Residual stresses reduction in filament wound composite
tubes. Journal of Reinforced Plastics and Composites 1999;18:684 – 701.
[8] Gentz M, Armentrout D, Rupnowski P, Kumosa L, Shin E, Sutter JK, et al. In-
plane shear testing of medium and high modulus woven graphite fiber
reinforced/polyimide composites. Composites Science and Technology
2004;64:203–20.
[9] Manders PW, Chou TW, Jones RR, Rock JW. Statistical analysis of multiple
fracture in 0º/90º/0º glass–fiber/epoxy–resin laminates. Journal of Material
Sciences 1983;8:352–62.
[10] Parvisi A, Bailey JE. On multiple transverse cracking in glass-fiber epoxy cross-
ply laminates. Journal of Material Sciences 1978;13:2131–6.
[16] Garrett KW, Bailey JE. Multiple transverse fracture in 90◦ cross-ply laminates of
a glass fibre-reinforced polyester. Journal of Composite Materials 1977;18:339–
56.
[17] Boniface L, Smith PA, Bader MG, Rezaifard AH. Transverse ply cracking in
cross-plyCFRP laminates. Initiation or propagation controlled. Journal of
Composite Materials 1997;31:1080–112.
[18] Bâzant ZP. Scaling of Structural Strength. London: Hermes Penton Science;
2002.
[22] Gamstedt EK, Sjören BA. An experimental investigation of the sequence effect
in block amplitude loading of cross-ply composite laminates. International
Journal of Fatigue 2002;24:437–46.
[24] Crossman FW, Wang ASD. The dependence of transverse cracking and
delamination on ply thickness in graphite/epoxy laminates. In: Reifsnider K,
editor. Damage in composite materials, 1982, p. 118 139.
254
/ -
[27] Jamison RD, Schulte K, Reifsnider KL, Stinchcomb WW. Characterization and
analysis of damage mechanisms in tension-tension fatigue of graphite/epoxyde
laminates. In: Reifsnider K, editor. Effects of defects in composites materials,
1984, p. 21 – 55.
[28] Liu S, Kutlu Z, Chang FK. Matrix cracking and delamination in laminated
composite beams subjected to a transverse concentrated line load. Journal of
Composite Materials 1993;27:436 – 470.
[33] Laurin F, Carrere N, Huchette C, Maire J-F. A multiscale hybrid damage and
failure approach for strength predictions of composite structures. Journal of
Composite Materials,sous presse, 2013.
[35] Maimí P, Camanho PP, Mayugo JA, Turon A. Mechanics of Materials Matrix
cracking and delamination in laminated composites . Part I : Ply constitutive law
, first ply failure and onset of delamination. Mechanics of Materials
2011;43:169–85.
[36] Maimí P, Camanho PP, Mayugo J a., Dávila CG. A continuum damage model for
composite laminates: Part I – Constitutive model. Mechanics of Materials
2007;39:897–908.
[37] Maimí P, Camanho PP, Mayugo J a., Dávila CG. A continuum damage model for
composite laminates: Part II – Computational implementation and validation.
Mechanics of Materials 2007;39:909–19.
[38] Eitzen D., Wadley HN. Acoustic emission: Establishing the fundamentals.
Journal of Research of the National Bureau Of Standarts 1984;89:75 – 100.
255
/ -
[43] Papargyris AD, Cooke R., Papargyri S., Botis A. No TitleThe acoustics
behaviour of bricks in relation to their mechanical behaviour. Construction and
Building Materials 2001;15:361 – 369.
[44] Scholey JJ, Wilcox PD, Wisnom MR, Friswell MI. Composites : Part A
Quantitative experimental measurements of matrix cracking and delamination
using acoustic emission. Composites Part A 2010;41:612–23.
[46] Uenoya T. Acoustic emission and the fiber/matrix interface failure progress and
strength in glass fiber/epoxy laminates. 10th ICCM, 1995, p. 581 – 588.
[51] Wisnom MR, Hallett SR. Composites : Part A The role of delamination in
strength , failure mechanism and hole size effect in open hole tensile tests on
quasi-isotropic laminates. Composites Part A 2009;40:335–42.
256
/ -
[52] Hallett SR, Green BG, Jiang WG, Wisnom MR. Composites : Part A An
experimental and numerical investigation into the damage mechanisms in
notched composites. Composites Part A 2009;40:613–24.
[53] Sutton MA, Wolters WJ, Peters WH, Ranson WF, McNeill SR. Determination of
displacements using an improved digital correlation method. Image Vision and
Computing 1983;1:133 – 139.
[56] Kim J, Pierron F, Wisnom MR, Avril S. Composites : Part A Local stiffness
reduction in impacted composite plates from full-field measurements.
Composites Part A 2009;40:1961–74.
[58] Hashin Z, Rotem A. Fatigue failure criterion for fiber reinforced materials.
Journal of Composite Materials 1973;3:448.
[60] Tsai SW, Wu EM. A general theory of strengh of anisotropic materials. Journal
of Composite Materials 1971;5:55 – 80.
[61] Puck A, Schu H. Failure analysis of FRP laminates by means of physically based
phenomenological models. Composite Science and Technology 2002;62:1633–
62.
[62] Soden PD, Hinton MJ, Kaddour AS. Lamina properties , lay-up configurations
and loading conditions for a range of fibre-reinforced composite laminates.
Composites Science and Technology 1998;58:1011–22.
[63] Soden PD, Hinton MJ, Kaddour AS. A comparison of the predictive capabilities
of current failure theories for composite laminates. Composite Science and
Technology 1998;58:1225 – 1254.
[64] Hinton MJ, Kaddour AS, Soden PD. A comparison of the predictive capabilities
of current failure theories for composite laminates , judged against experimental
evidence. Composites Science and Technology 2002;62:1725–97.
257
/ -
[65] Kachanov L. On the time to failure under creep conditions. Izv AN SSSR, Otd
Tekhn Nauk 1958;8:26 – 31.
[66] Chaboche J. Continuum Damage Mechanics, parts I and II. Journal of Applied
Mechanics 1989;55:55 – 59, 55 – 65.
[67] Pachajoa ME, Frances MK, Lee JD. Stress and failure analysis of composites
structures. Engineering Fracture Mechanics 1995;50:883–902.
[69] Gasser a., Ladeveze P, Peres P. Damage modelling for a laminated ceramic
composite. Materials Science and Engineering: A 1998;250:249–55.
[72] Wang ASD, Crossman FW. Initiation and growth of transverse cracks and edge
delamination in composite laminates. Part I: an energy method. Journal of
Composite Materials 1980;14:71 – 87.
[74] Wang ASD, Chou PC, Lei S. A stochastic model for the growth of matrix cracks
in composite laminates. Journal of Composite Materials 1984;18:239–54.
[75] Fukunaga H, Chou TW, Peters PWM, Schulte K. Probabilistic failure strength
analyses of graphite/epoxy cross-ply laminates. Journal of Composite Materials
1984;18:339–56.
[77] Berthelot JM, Le Corre JF. Statistical analysis of the progression of transverse
cracking and delamination in cross-ply laminates. Composite Science and
Technology 2000;60:2659 – 2669.
258
/ -
[82] Thollon Y, Hochard C. A general damage model for woven fabric composite
laminates up to first failure. Mechanics of Materials 2009;41:820–7.
[83] Maire JF, Chaboche JL. A new formulation of continuum damage mechanics
(CDM) for composite materials. Aerospace Science and Technology 1997;1:247–
57.
[84] Chaboche J-L, Maire J-F. A new micromechanics based CDM model and its
application to CMC’s. Aerospace Science and Technology 2002;6:131–45.
[85] Laurin F, Carrère N, Maire J-F. A multiscale progressive failure approach for
composite laminates based on thermodynamical viscoelastic and damage models.
Composites Part A: Applied Science and Manufacturing 2007;38:198–209.
[86] Petipas C-R. Analyse et prévision du comportement à long terme des composites
à fibres de carbone/matrie organique". UFR Sciences et techniques de
la’Université de Franche - Comté, 2000.
[89] Irisarri F-X, Laurin F, Carrere N, Maire J-F. Progressive damage and failure of
mechanically fastened joints in CFRP laminates – Part I: Refined Finite Element
modelling of single-fastener joints. Composite Structures 2012;94:2269–77.
[90] Irisarri F-X, Laurin F, Carrere N, Maire J-F. Progressive damage and failure of
mechanically fastened joints in CFRP laminates – Part II: Failure prediction of an
industrial junction. Composite Structures 2012;94:2278–84.
259
/ -
[93] Coleman D. On the strength of classical fibres and fibre bundles. Journal of the
Mechanics and Physics of Solids 1958;7:60 – 70.
[102] Tagarielli VL, Minisgallo G, Mcmillan AJ, Petrinic N. The response of a multi-
directional composite laminate to through-thickness loading. Composites Science
and Technology 2010;70:1950–7.
260
/ -
[105] Lemaitre J, Chaboche JL. Mécanique des matériaux solides. In: Dunod editor
Paris: 1992.
[107] Krajcinovic D, Fonseca GU. The continuous damage theory of brittle materials -
Parts I and II. Journal of Applied Mechanics 1981;48:809 – 824.
[108] Onat ET, Leckie FA. Representation of mechanical behavior in the presence of
changing internal structure. Journal of Applied Mechanics 1988;55:1 – 10.
[112] Spencer AJM. Theory of invariants. In: Eringen C, editor. Contiuum Mechanics,
New York: Academic Press; 1971, p. 239 – 353.
[114] Spencer AJM. The formulation of constitutive equations for anisotropic solids.
In: Boehler JP, editor. Euromech 115, 1982, p. 3 – 26.
[118] Halm D, Dragon A, Charles Y. A modular damage model for quasi-brittle solids
– interaction between initial and induced anisotropy. Archive of Applied
Mechanics 2002;72:498–510.
261
/ -
[119] Maire JF, Lesne PM. A damage model for ceramic matrix composites. Aerospace
Science and Technology 1997;1:259–66.
[120] Dávila CG, Rose CA, Camanho PP. A procedure for superposing linear cohesive
laws to represent multiple damage mechanisms in the fracture of composites
2009:211–23.
[122] Bazant ZP, Oh BH. Crack band theory for fracture of concrete. Materials and
Structures 1983;16:155 – 177.
262
/ -
[132] Bazant ZP, Lin FB. Nonlocal smeared cracking model for concrete fracture.
Journal of Engineering Mechanics 1988;114:2493 – 2510.
[134] Jirasek M, Bazant ZP. Localization analysis of nonlocal model based on crack
interactions. Journal of Engineering Mechanics 1994;120:1521 – 1542.
[135] Bazant ZP, Xiang Y, Adley M, Prat PC, Akers SA. Microplane model for
concrete: II. Data delocalization and verification. Journal of Engineering
Mechanics 1996;122:255 – 262.
[136] Jirasek M. Nonlocal models for damage and fracture: comparison of approaches.
International Journal of Solids and Structures 1998;35:4133 – 4145.
[143] Hallett SR, Jiang W, Khan B, Wisnom MR. Modelling the interaction between
matrix cracks and delamination damage in scaled quasi-isotropic specimens.
Composites Science and Technology 2008;68:80–9.
[144] Arellano MT, Crouzeix L, Douchin B, Collombet F, Moreno HH, Velázquez JG.
Strain field measurement of filament-wound composites at ± 55º using digital
image correlation : An approach for unit cells employing flat specimens.
Composite Structures 2010;92:2457–64.
[145] Hallett SR, Green BG, Kin WJ, Cheung H, Wisnom MR. The open hole tensile
test : a challenge for virtual testing of composites 2009:169–81.
263
/ -
[146] O’Higgins RM, McCarthy M a., McCarthy CT. Comparison of open hole tension
characteristics of high strength glass and carbon fibre-reinforced composite
materials. Composites Science and Technology 2008;68:2770–8.
[148] Mulyukov RR. Rapid change of stresses in thickness direction in long orthotropic
tube under internal pressure and axial load. Acta Mechanica 2010;211:323–36.
[152] Xu P, Zheng JY, Liu PF. Finite element analysis of burst pressure of composite
hydrogen storage vessels. Materials and Design 2009;30:2295–301.
[153] Thionnet A, Bunsell A, Camara S, Allen DH. A simplified FE2 method to model
the burst of composite pressure vessels. Proc. JNC 17, 2011.
[154] Wang R, Jiao W, Liu W, Yang F. A new method for predicting dome thickness
of composite pressure vessels. Journal of Reinforced Plastic Composites
2010;29:3345 – 3352.
264
/ -
[160] Chou P, Croman R. Residual strength in fatigue based on strength - life equal
rank assumption. Journal of Composite Materials 1978;12:177 – 194.
[162] Yang J, Liu M. Residual strength degradation model and theory of periodic proof
tests for graphite/epoxy laminates. Journal of Composite Materials 1977;11:176 –
203.
[163] Yang J, Liu M. Fatigue and residual strength degradation for graphite/epoxy
composites under tension - compression cyclic loading. Journal of Composite
Materials 1978;12:19–39.
[164] Yang J, Sun C. Proof test and fatigue of unnotched composite laminates. Journal
of Composite Materials 1980;14:168 – 74.
[165] Yang J, Jones D. Effect of load sequence on the statistical fatigue of composites.
AIAA 1980;18:1525 – 1531.
[166] Yang J, Jones D. Load sequence effects on the fatigue of unnotched composite
laminates. ASTM STP 1981;723:213 – 232.
[168] Yang J, Miller R, Sun C. Effect of high load on statistical fatigue of unnotched
graphite/epoxy laminates. Journal of Composite Materials 1980;14:82 – 94.
[169] Sendeckyj G. Life prediction for resin - matrix composite materials. Composite
Materials Series 1991;4:431–83.
[170] Sendeckyj G. Fitting models to composite materials fatigue data. ASTM STP
1981;734:245–60.
[171] Adam T, Dickson R, Jones C, Reiter H, Harris B. A power law fatigue damage
model for fiber -reinforced plastic laminates. Instn Mech Engrs, 1986, p. 155–66.
[172] Diao X, Lessard L, Shokrieh M. Statistical model for multiaxial fatigue behavior
of unidirectional plies. Composite Science and Technology 1999;59:2025–35.
265
/ -
[174] Gao Z. A cumulative damage model for fatigue life of composite laminates.
Journal of Reinforced Plastic Composites 1994;13:128 – 141.
[176] Hochard C, Thollon Y. A generalized damage model for woven ply laminates
under static and fatigue loading conditions. International Journal of Fatigue
2010;32:158–65.
266
Caractérisation et modélisation de l’endommagement des composites bobinés.
Application à la prédiction de l’éclatement des réservoirs bobinés hyperbares
Un modèle d’endommagement dédié aux composites bobinés est développé à partir des outils de la
mécanique de l’endommagement continu, de la thermodynamique des processus irréversibles et de
la théorie de représentation des fonctions tensorielles. La particularité de ce modèle est l’utilisation
d’une approche à directions fixes de l’endommagement qui associe à chaque mode de dégradation
des variables internes scalaires et des tenseurs directionnels. La rupture des fibres (considérée
comme probabiliste), les fissurations tant matricielles, que hors – plan ou provoquées par le
cisaillement sont ainsi prises en compte. Le modèle est capable de reproduire, dans un contexte
tridimensionnel imposé par les fortes épaisseurs de composite, la perte de rigidité, l’interaction entre
anisotropies initiale et induite, la non linéarité du comportement, les déformations résiduelles et la
viscosité en cisaillement. Afin de valider cette approche, le comportement thermomécanique
d’éprouvettes issues de structures bobinées a été caractérisé grâce à des essais de traction multi –
instrumentés (vidéo – traction, émission acoustique,...). On montre que le modèle est capable non
seulement de simuler la réponse mécanique macroscopique de ces échantillons, mais également de
reproduire l’émission acoustique enregistrée, de distinguer les différentes formes
d’endommagement et de prédire précisément l’éclatement des réservoirs hyperbares.