Fatigue des composites tissés : étude et modélisation
Fatigue des composites tissés : étude et modélisation
&OWVFEFMPCUFOUJPOEV
%0$503"5%&-6/*7&34*5²%&506-064&
%ÏMJWSÏQBS
Institut Supérieur de l’Aéronautique et de l’Espace (ISAE)
1SÏTFOUÏFFUTPVUFOVFQBS
Thomas ROUAULT
le mardi 18 juin 2013
5JUSF
6OJUÏEFSFDIFSDIF
Institut Clément Ader
%JSFDUFVS T EFʾÒTF
Jury :
Je pense d'abord à Christophe Bouvet mon directeur de thèse, qui a assuré seul
mon encadrement scientifique. Ses nombreuses idées et sa détermination à
l'aboutissement de ces travaux, ainsi que ses conseils m'ont beaucoup aidé pendant ces
3 années.
Je remercie mes rapporteurs, Zoheir Aboura et Eric Martin, ainsi que Frédéric
Laurin, qui ont accepté et honoré la lourde tâche de relire consciencieusement ce
mémoire, et pour leur participation à mon jury. J'adresse également mes remerciements
à Jean-Jacques Barrau pour avoir présidé ce jury.
La fin de ces travaux s’est ponctuée par une éprouvante période de rédaction que
j’ai achevée sur le site de Marignane, et j'adresse un grand merci à l'équipe du service
"pale" du bureau d'étude, à ses permanents, ses anciens, et à ses stagiaires, pour leur
bonne humeur et leur sympathie.
Je n'oublie pas ceux qui ont suivi cette période de thèse de plus ou moins loin, ou
qui sont venu assister à ma soutenance. Merci à Thibaud, Hindu, Olivier, Pascal, Jpeg,
Edouard, Damien, Jean pour les bons moments passés ensemble.
4
Table des matières
5
[Link]. Validation du montage .................................................................................................... 63
[Link]. Comportement du matériau en cisaillement ......................................................... 65
[Link]. Endommagement et pseudo-plasticité..................................................................... 67
[Link]. Comportement en fatigue .............................................................................................. 69
[Link]. Endommagement microscopique ............................................................................... 70
[Link]. Etude de l’endommagement par thermographie infrarouge .......................... 73
2.3.4. Discussion ..................................................................................................................................... 74
2.4. Essais de couplage d’endommagement ...................................................................... 75
2.4.1. Principe .......................................................................................................................................... 75
[Link]. Aperçu des études sous chargement cyclique variable ..................................... 75
[Link]. Description des essais successifs ................................................................................ 77
2.4.2. Couplage traction-traction ..................................................................................................... 77
[Link]. Essais de fatigue en traction sens 1 ........................................................................... 77
[Link]. Couplage (ε1 – ε2) quasi-statique................................................................................. 81
[Link]. Couplage (ε1 – ε2) fatigue ................................................................................................ 82
2.4.3. Couplage traction-cisaillement ............................................................................................ 85
2.4.4. Couplage cisaillement-traction ............................................................................................ 87
2.4.5. Discussion ..................................................................................................................................... 89
6
Chapitre IV – Modélisation numérique de la propagation de coupure........139
4.1. Bibliographie .....................................................................................................................140
4.1.1. Mécanique linéaire de la rupture et méthodes analytiques .................................. 140
[Link]. Eléments de mécanique de la rupture ................................................................... 140
[Link]. Direction de propagation ............................................................................................ 145
4.1.2. Méthodes empiriques ........................................................................................................... 146
4.1.3. Modélisations multi-échelles ............................................................................................. 148
4.1.4. Endommagement continu ................................................................................................... 152
4.1.5. Endommagement discret .................................................................................................... 154
4.1.6. Eléments finis avancés.......................................................................................................... 156
[Link]. X-FEM .................................................................................................................................. 156
[Link]. Formulation péridynamique ..................................................................................... 159
4.2. Principe du modèle à l’échelle de la mèche ............................................................ 161
4.2.1. Aspects pratiques du maillage........................................................................................... 162
4.2.2. Rôle des éléments ................................................................................................................... 165
4.2.3. Cas de maillages à plusieurs plis ...................................................................................... 165
4.3. Comportement du stratifié sain ..................................................................................166
4.4. Prise en compte de l’évolution du comportement mécanique ........................168
4.4.1. Endommagement en fatigue .............................................................................................. 169
[Link]. Données expérimentales ............................................................................................. 169
[Link]. Matrice de raideur de l’élément d’interface ........................................................ 171
[Link]. Relation raideur – endommagement ..................................................................... 173
[Link]. Loi d’endommagement en fatigue ........................................................................... 175
[Link]. Déformation prise en compte dans le calcul de la raideur ............................ 178
[Link]. Cumul d’endommagement ......................................................................................... 180
4.4.2. Pseudo-plasticité..................................................................................................................... 183
[Link]. Loi d’endommagement en cisaillement ................................................................ 184
[Link]. Modèle de pseudo-plasticité endommageable ................................................... 185
4.5. Rupture des éléments .....................................................................................................186
4.5.1. Calcul de la durée de vie d’une mèche ........................................................................... 187
4.5.2. Pilotage des éléments de rupture .................................................................................... 188
4.5.3. Propagation de l’endommagement ................................................................................. 190
4.6. Eléments d’interface 0-45 ............................................................................................. 191
4.6.1. Principe ....................................................................................................................................... 191
4.6.2. Identification de la matrice de raideur .......................................................................... 193
[Link]. Matrice d’interpolation I ............................................................................................. 194
[Link]. Matrice de rigidité KM ................................................................................................... 195
[Link]. Matrice d’extrapolation des efforts J ...................................................................... 196
[Link]. Assemblage de la matrice ........................................................................................... 197
[Link]. Dégradation – rupture .................................................................................................. 197
4.6.3. Discussion .................................................................................................................................. 198
4.7. Intégration sur un modèle générique .......................................................................198
4.7.1. Aspects pratiques de l’algorithme ................................................................................... 198
4.7.2. Modélisation d’une structure............................................................................................. 202
7
Chapitre V – Comparaison essais - calculs .............................................................205
5.1. Introduction – Méthodologie .......................................................................................206
5.1.1. Précisions concernant les modèles ................................................................................. 206
5.1.2. Dépouillement .......................................................................................................................... 207
5.2. Essais de propagation en traction ..............................................................................207
5.2.1. Cas du drapage [0/90]n ........................................................................................................ 207
[Link]. Eprouvette de 50 mm de large.................................................................................. 208
[Link]. Eprouvettes de 30 mm de large ............................................................................... 214
[Link]. Discussion.......................................................................................................................... 218
5.2.2. Cas du drapage [±45]2........................................................................................................... 218
5.2.3. Cas du drapage [0;45]s .......................................................................................................... 226
5.3. Essais de propagation en cisaillement......................................................................231
5.3.1. Cas du drapage [0/90]4 ........................................................................................................ 231
5.3.2. Cas du drapage [0;45]s .......................................................................................................... 236
5.4. Synthèse et discussion ....................................................................................................240
Références bibliographiques......................................................................................249
8
Chapitre I – Introduction générale
Les pales d’hélicoptères sont des pièces dont la rupture, même partielle, est
catastrophique pour l’appareil, ce qui les place dans la catégorie des composants les plus
critiques de l’hélicoptère. Une pale peut schématiquement se décrire par une structure poutre,
orientée selon l’axe de l’envergure, et constituée d’un longeron et d’un arêtier de rigidité
importante, de remplissages non-travaillants, d’une ou plusieurs nervures (éventuellement), et
recouverte d’une peau appelée revêtement (Figure I.1).
Les travaux présentés ici se sont concentrés sur un matériau donné, un tissu pré-
imprégné verre/époxy employé comme revêtement de pale mais également sur d’autres pièces
de l’hélicoptère (poutre de queue, dérive) ou pour des applications telles que les pales
d’éoliennes.
Outre son chargement de service, le revêtement d’une pale est soumis aux aléas de
l’environnement extérieur, tel que les risques d’impacts [Navarro12], ou de foudre. Il peut
également être le siège de défauts de matériau, de production, de maintenance ou de
concentrations de contraintes inattendues. Ces phénomènes peuvent être à l’origine de rupture
des plis dans l’épaisseur. On désignera par le terme « coupure », la rupture du pli à travers son
9
épaisseur, c'est-à-dire avec rupture de mèches, pour lever l’ambigüité avec les « fissures » qui,
pour les matériaux composites, concernent dans la majorité des cas le délaminage (rupture
interlaminaire), les décollements fibre/résine ou la rupture matricielle dans les mèches. On
parle également de rupture translaminaire ou « through-the-thickness crack » en anglais pour
désigner les coupures.
Ces travaux s’inscrivent dans cette problématique générale, peu étudiée dans la
littérature. Ils font suite et s’appuient sur les travaux de thèse de M. Bizeul [Bizeul09] portant
sur la même problématique de l’étude expérimentale et numérique du phénomène de
propagation de coupure sur revêtement mince et sous chargement de vol.
Le chargement de vol d’une pale est décomposé en un effort centrifuge constant et des
charges aérodynamiques. Ces dernières sont cycliques pour des raisons intrinsèques à la
mécanique du vol de l’hélicoptère et de par les vibrations engendrées par les parties tournantes.
Sur la structure on a coutume de négliger les efforts tranchants et on décompose les charges
aérodynamiques en trois moments selon les axes du repère lié à la pale : battement, trainée et
torsion (Figure I.2).
Figure I.2 – Définition des axes du repère pale et des actions mécaniques s’exerçant sur une
pale en vol.
10
Chapitre I – Introduction
Sur le revêtement, cela se traduit par des sollicitations planes de traction pour l’effort
centrifuge, les moments de battement et de trainée, et de cisaillement pour le moment de
torsion. La démarche de l’étude expérimentale de la propagation de coupure sur un revêtement
de pale a été de considérer séparément les sollicitations de traction et de cisaillement
(Figure I.3), de les analyser expérimentalement, et de proposer un modèle capable de reproduire
le comportement expérimental. Différents cas de drapages et de sollicitations ont été étudiés,
permettant ainsi d’évaluer les résultats d’un modèle sur un panel varié de situations.
On peut préciser ici que la structure de la pale, en particulier la présence d’un longeron
concentrant la majorité de la rigidité de l’éprouvette, a pour conséquence la conservation de la
déformation dans le revêtement en cas de propagation de coupure. En dehors du voisinage de la
pointe de coupure, siège de concentrations de contraintes, la déformation du revêtement n’est
11
pas affectée par la propagation de la coupure. Ainsi, il est apparu judicieux de réaliser l’étude
expérimentale à déformation imposée, et d’utiliser la grandeur « déformation » plutôt que
« contrainte » dans les différentes analyses.
Pour les étudier, des essais de caractérisation sur éprouvettes non entaillées
accompagnent l’étude de la propagation de coupure afin d’obtenir des données matériau. Des
essais sous sollicitations quasi-statiques et cycliques en traction et en cisaillement ont été menés
ainsi que des observations microscopiques des éprouvettes endommagées. D’autres essais
originaux ont été réalisés pour étudier le couplage entre deux types d’endommagement
différents. Plus précisément, des éprouvettes ont subi successivement différents types de
sollicitation afin d’évaluer l’effet de l’endommagement obtenu sous une première sollicitation
donnée sur le comportement sous une seconde sollicitation.
Les études expérimentales ont orienté une stratégie de modélisation qui constitue un
second volet de ces travaux de thèse.
Elle repose sur une représentation originale du matériau à l’échelle des mèches en
s’appuyant sur des constats expérimentaux. La modélisation du comportement du matériau à
travers une représentation semi-discrète de l’endommagement s’inspire des observations
expérimentales. Ce modèle tire son originalité du critère de rupture qui permet la propagation :
il s’appuie sur la tenue en fatigue des fibres, et est identifié uniquement sur des essais de fatigue
dans l’axe des fibres. En particulier il ne nécessite pas l’introduction de paramètres obtenus par
identification inverse ou de valeur de taux de restitution d’énergie critique nécessitant des essais
de propagation coûteux à mettre en œuvre, et à exploiter, en comparaison des essais à réaliser
pour obtenir une courbe de fatigue S-N ou ε-N.
12
Chapitre I – Introduction
Ces travaux de thèse ont ainsi comporté une partie expérimentale, et une partie
modélisation numérique, et le plan de ce mémoire suit ce découpage. Il est organisé en 4
chapitres :
Les chapitre II et III traitent de l’étude expérimentale menée sur le matériau. Le Chapitre
II, présente l’étude de l’endommagement en traction, cisaillement, et le couplage possible entre
les endommagements. L’étude du comportement du matériau entaillé sous sollicitation cyclique,
est ensuite présentée dans le Chapitre III.
13
Chapitre II – Comportement mécanique
et endommagement
chapitre2
15
2.1. Introduction
L’ensemble des travaux de thèse détaillés dans ce mémoire concerne un unique
matériau, composé de fibres de verre et d’une matrice époxy ordonnées suivant une armature
tissée. Utilisé comme revêtement d’une pale d’hélicoptère il est soumis à des sollicitations
cycliques combinées de traction et de cisaillement. La présence d’un défaut géométrique
détourne les chemins d’efforts et provoque ainsi une concentration de contrainte en pointe
d’entaille. Les champs de contraintes et de déformations perturbés et amplifiés peuvent
provoquer un endommagement important sous diverses formes et affectent la tenue mécanique
du matériau.
16
Chapitre II – Comportement mécanique et endommagement
Paramètres
Les points de tissages sont alignés selon des directions particulières qui ne sont pas
symétriques par rapport à l’un des axes d’orthotropie (voir Figure II.2), ce qui a conduit à des
voilages de plaques après polymérisation. [Bishop89] constate cet effet, et l’insuffisance de la
notation couramment utilisée pour décrire la séquence d’empilement. Elle propose une notation
plus précise présentée Figure II.2 pour expliquer le voilage de plaques obtenu pour certains
empilements équilibrés présentant a priori une symétrie miroir. Le premier et second nombre
indiquent les orientations de renforts majoritaires sur les faces inférieures et supérieures du pli
respectivement. Les crochets sous les nombres indiquent le tissage entre les deux directions. Les
lettres sur les chiffres indiquent la nature (« chaîne » ou « trame ») des mèches correspondant à
l’orientation indiquée. Dans le motif du satin de 8, les points de tissages sont alignés selon un axe
orienté à environ 18° de l’axe des mèches. La Figure II.2 permet de constater que ces directions
changent de signe en pivotant le pli de 90° ou en le retournant pour permuter les faces visibles
et cachées mais sans changer l’orientation.
17
Figure II.2 - Différentes configurations de drapages possibles avec un pli 0/90, et
introduction des notations de [Bishop89]. La couleur utilisée pour représenter les mèches
caractérise sa direction « chaîne » ou « trame ».
Pour les drapages de notre étude, la symétrie miroir n’est pas appliquée, les plis sont
drapés sans permuter les faces visibles et cachées dans le but de limiter le voilage des plaques.
Nous n’utiliserons pas les notations de [Bishop89] et les drapages utilisés seront notés :
- [0/90]n signifiant n plis au total, comportant chacun des renforts dans les directions
0° et 90°. Les éprouvettes notées [0]n (resp. [90]n), indiquent que l’effort s’applique
dans la direction chaîne (resp. trame).
- [±45]n signifiant n plis au total comportant chacun des mèches dans les directions
+45° et -45°.
18
Chapitre II – Comportement mécanique et endommagement
Figure II.3 - Mise en évidence du phénomène d'imbrication des plis par comparaison
essai/simulation de l’épaisseur moyenne d’un pli en fonction du nombre total de plis [Lomov03].
19
Figure II.4 – cycle de polymérisation du pré-imprégné verre-époxy étudié.
On distingue les deux directions de tissages par les termes « chaîne » et « trame », et le
tissu est qualifié d’équilibré lorsque les mèches de chaîne et de trame contiennent le même
nombre de fibres moyen.
Le comportement des tissus composites dans le cas de traction dans un des deux axes de
renforts a déjà été étudié par de nombreux auteurs, dans le cas de sollicitations quasi-statiques.
Dans le cas où les mèches sont constituées de fibres de verre, et de matrice organique la courbe
σ-ε typique est donnée Figure II.5.
20
Chapitre II – Comportement mécanique et endommagement
Pour les tissus à base de fibres de carbone, ce phénomène est moins marqué, voir
inexistant où même inverse selon les constituants. A titre d’exemple, [Daggumati13] observe une
inflexion dans l’autre sens (rigidification). En effet, les fibres des tissés ont la particularité de
pouvoir se déplier avec l’accroissement de la sollicitation et le développement de
l’endommagement matriciel. Ce dépliement peut parfois s’observer pour les unidirectionnels
[VanDreumel77] dont les fibres présentent un défaut d’alignement par rapport à la direction de
chargement. La concurrence entre les phénomènes de dégradation de la résine et ceux
d’alignement des fibres sur la réponse contrainte/déformation semble tourner en faveur d’une
baisse de rigidité pour les renforts en verre ([Alif97], [Kumagai04], [Osada03], [Steeves04]
[John01] [Lomov08]), mais pas forcément pour les renforts en carbone. [Gao99], [Daggumati13],
[Alif97] n’observent pas d’affaiblissement de rigidité en traction, contrairement à [Karahan11]
qui observe une diminution du module d’Young tangent sur sergé carbone/époxy. Les fibres de
carbone étant plus rigides que les fibres de verre, la matrice et a fortiori sa rupture a moins
d’effet sur la raideur totale du composite. De plus [Alif97] rappelle que les fibres de verre
travaillent sur une plage de déformations plus élevées que le carbone.
La rupture finale du tissu correspond à la rupture fragile et brutale des fibres pour des
déformations typiquement comprises entre 1,5 et 3% pour les tissus de verre.
Par suite, de nombreux travaux se sont concentrés sur l’endommagement des tissus à
l’échelle microscopique. Dans le cas de traction quasi-statique, on peut citer [Manger97],
[Gao99], [Osada03], [Daggumati10], [Karahan11] et [DeGreef11a] qui décrivent la nature des
endommagements constatés, ainsi que leurs chronologies. Les types de dommage relevés sont :
21
en résine créent des zones de sur-contraintes propices à la fissuration [DeCarvalho11]. Ainsi la
fissuration matricielle apparait dans des zones bien particulières, généralement dans les zones
d’ondulation des mèches transverses. [Daggumati10] classe ces phénomènes d’endommagement
en fonction de leur échelle: l’échelle microscopique concerne les ruptures d’interface entre fibres
et résine ou la rupture de fibres isolées. L’échelle mésoscopique se rapporte aux fissures de la
taille des mèches, et aux méta-délaminages. Enfin, on peut associer les délaminages étendus,
ainsi que les ruptures longitudinales de plusieurs mèches à l’échelle macroscopique.
Figure II.6 – Illustration des modes d’endommagement observés sur les composites tissés
[Alif97].
Les différents auteurs qui se concentrent sur l’endommagement des tissus, cherchent à
mesurer les seuils des différents types d’endommagement, caractériser leurs formes, et
l’influence de certains paramètres d’essais. Ces endommagements apparaissent avec des
chronologies différentes. L’initiation se fait généralement par fissuration dans les mèches
transverses au niveau des zones d’ondulation. A l’aide de corrélation d’image, [Daggumati11]
parvient à corréler le champ de déformation au motif de tissage et montre que la déformation
maximale se situe au centre des points de tissages, lieu où il détecte les premières fissurations
matricielles (Figure II.7). Par cette méthode, il parvient également à détecter une discontinuité
temporelle nette dans la déformation locale signe d’apparition brutale de micro-fissuration.
[Karahan11] utilise la même méthode et obtient des résultats similaires sur traction dans l’axe et
hors axes des mèches d’un sergé.
22
Chapitre II – Comportement mécanique et endommagement
Figure II.7 - Champ de déformation longitudinale (selon x) obtenu par corrélation d’images
pour une déformation moyenne de 1,19% sur un satin de 5 carbone/PPS. Les zones de
déformations maximales se superposent aux points de tissages [Daggumati11].
[Gao99] a étudié un tissu satin de carbone en traction dans l’axe des fibres, et l’influence
du nombre de plis sur l’évolution des coefficients élastiques et de l’endommagement. Cet
endommagement est détecté par émission acoustique, rayon X, ou microscopie. Les
microfissures détectées traversent entièrement les mèches dans leur épaisseur. Elles se
propagent à l’interface entre les deux directions de renforts créant un méta-délaminage dans les
zones d’ondulations, selon le scénario également décrit par [Pandita01] en fatigue. La baisse de
rigidité constatée, l’évolution du coefficient de Poisson et les déformations résiduelles sont plus
importantes pour des éprouvettes de 2 plis que pour 4 ou 6 plis. Les mèches transverses du pli
externe se fissurent plus tôt en déformation et ces plis extérieurs ont un rôle plus important
dans la tenue mécanique qu’en présence de 4 ou 6 plis. De plus, la séparation des directions de
renfort décharge totalement les mèches du pli extérieur, ce qui, d’après l’auteur, peut expliquer
cette différence.
23
Figure II.8 – (a) : schéma du mécanisme de fissuration transverse dans un tissu, (b) :
différence de géométrie entre taffetas (plain) et satin, rendant plus difficile la mise en flexion des
mèches transverses du satin ce qui retarde la fissuration [Osada03].
24
Chapitre II – Comportement mécanique et endommagement
[Bizeul09] a rassemblé des données de fatigue en traction sur tissu [0/90]n à base de
fibres de verre disponibles dans la littérature (Figure II.9).
Figure II.9 – Synthèse des durées de vie en fatigue sur des tissus à base de fibres de verre
pour des drapages [0/90]n [Bizeul09].
Ces essais ont pour la plupart été réalisés à contrainte imposée pour des rapports de
charge R = σmin/σmax généralement de 0,1, pour des stratifiés verre/époxy. Seuls [Fujii93],
[Echtermeyer95] et [Dyer98] ont utilisé une résine polyester. Les essais de [Echtermeyer95] ont
par ailleurs été réalisés avec un rapport de charge R=-1 (traction – compression alternée).
A une contrainte donnée, la durée de vie du stratifié ou le nombre de cycle à rupture, Nf,
peut présenter une dispersion d’un demi ordre de grandeur (série [Pandita01] sur la Figure II.9)
et parfois supérieure. Néanmoins, la corrélation entre Nf et la sollicitation mécanique
25
(représentée à titre d’exemple par σmax) s’approxime généralement bien par une loi
logarithmique de type Wöhler (II.1) ou puissance de type Basquin (II.2) (Figure II.9):
(II.1)
(II.2)
pour lesquels k d’une part, et σ6 et γ d’autre part sont des constantes à identifier. Pour Nf
en mégacycles (Mcy), σ6 représente la contrainte correspondant à une rupture en 106 cycles.
Les essais à déformations imposées sont relativement rares dans la littérature. On peut
néanmoins citer les travaux de [VanPaepegem01a], qui a développé un banc de flexion cyclique
(Figure II.10), permettant l’étude à déplacement imposé. Il génère des sollicitations non-
uniformes dans l’épaisseur (flexion) et dans la longueur de l’éprouvette, pour étudier le
comportement en fatigue d’un tissu verre/époxy. Il exploite cet essai structural pour valider un
modèle d’endommagement en fatigue.
26
Chapitre II – Comportement mécanique et endommagement
Figure II.11 – Evolution caractéristique du module d’Young en fatigue pour des stratifiés
[0/90]n selon [Demers98a].
Ce dernier fait parti des auteurs qui évoquent la notion d’endommagement critique
(dans le sens dégradation du module), ou d’état caractéristique d’endommagement, et qui
correspond à la fin de la seconde phase, au-delà de laquelle la rupture a lieu brutalement. Il
étudie l’évolution de l’endommagement en fonction du nombre de cycles de fatigue avec R=-1
sur deux tissus de verre à résine phénolique ou polyester. L’auteur montre que si l’amplitude de
chargement est élevée, le module d’Young s’affaiblit régulièrement dans un premier temps,
jusqu’à atteindre une valeur critique mesurée à 0,77.E0, à partir de laquelle, il décroit
brutalement (Figure II.12a au-delà de 90MPa). Cette valeur d’endommagement critique apparait
indépendante du chargement, de sorte que si le chargement est plus faible, et que le module
27
d’Young n’atteint pas cette valeur, la rupture n’est pas observée (Figure II.12b en deçà de
41MPa). On relève par ailleurs pour ces chargements, un temps d’initiation de
l’endommagement en début d’essai pour lequel le module est stable.
28
Chapitre II – Comportement mécanique et endommagement
La résistance résiduelle est l’autre grandeur estimée notamment pour être utilisée en
modélisation. Les stratifiés à base de fibres de carbone présentent peu, voire pas de sensibilité à
la fatigue. Par exemple [Daggumati13] ne constate aucun abattement dans la résistance quasi-
statique après fatigue à des rapports de contrainte S = σmax / σUTS > 0,8 sur 106 cycles. La
résistance résiduelle des stratifiés composés de fibres de verre est quant à elle plus sensible à la
fatigue. Cependant, pour des chargements faibles n’excédant pas le knee-point (Figure II.5),
[Osada03] note que la résistance résiduelle n’est pas impactée par la fatigue. Ce point particulier,
seuil de l’endommagement quasi-statique semble correspondre à la limite d’endurance
([Kumagai05]).
29
pendant la première étape d’affaiblissement du module. En plus de la fissuration transverse, du
méta-délaminage, et du délaminage, [Daggumati13] observe sur un satin de 5 carbone/PPS, de
nombreuses ruptures de fibres longitudinales (Figure II.14). Les différents auteurs s’accordent
sur le fait que l’endommagement des tissus est spécifiquement lié au motif de tissage. La
position des fissures dans les mèches transverses (au centre de celles-ci, sur les bords, ou entre
les mèches) diffère selon les matériaux, sans doute du fait de la nature des constituants, de leur
teneur, et de la forme des mèches qui induit une répartition des contraintes propre à leurs
géométries.
Figure II.14 – Endommagement en fatigue d’un satin de 5. (a), fissuration des mèches
transverses, (b), rupture de fibres axiales et méta-délaminage entre les deux directions de renforts
[Daggumati13].
Il ressort des différentes sources que la rupture en fatigue apparait dominée par les
fibres. [Gagel06] précise que la diminution de rigidité et la rupture en fatigue sont des
mécanismes différents et découplés. [Curtis89] évoque également ce découplage entre
propriétés fibres et résine.
Cependant la rupture des fibres dépend de la contrainte qu’elles doivent soutenir, qui
résulte de l’état d’endommagement du matériau. [Kawai96] étudie la fatigue de stratifiés de 9
plis de tissu taffetas carbone/époxy. Certaines de ses éprouvettes ont présenté du délaminage et
il observe qu’elles sont nettement plus endurantes que celles qui n’ont pas généré ce type
d’endommagement. Il conclue que le délaminage relaxe les concentrations de contraintes et
retarde la rupture des fibres en fatigue.
30
Chapitre II – Comportement mécanique et endommagement
Des moyens de visualisation autres que la microscopie optique ou le MEB ont été
employés pour observer l’endommagement local. [Lambert11] utilise la tomographie 3D pour
contrôler l’endommagement et les porosités d’un non-tissé verre/époxy. Cette technique lui
permet de reconstruire un volume d’éprouvette de 25x25x6,5 mm3, d’en détecter la fissuration
matricielle, le délaminage, et la rupture de fibres, et d’avoir un aperçu tridimensionnel de ces
défauts et des porosités (Figure II.15a).
Figure II.15 – (a) : Image d’un stratifié verre/époxy obtenue par tomographie [Lambert11].
(b) : Champ de déphasage entre signal mécanique et signal thermique (en °) d’un sergé 2x2
verre/époxy en fatigue de traction, à différents instants [Fruehmann10]. Ces cartes révèlent le
motif du sergé, ainsi que le développement d’une fissure dans une des mèches transverses à la
direction de la charge.
31
effet sur la façon dont le matériau se dégrade et la vitesse à laquelle les directions de renforts se
dissocient. [Kawai96] cite l’importance des propriétés de la résine et de l’interface fibre/matrice
ainsi que l’état de contraintes résiduelles obtenu durant la fabrication. Mais l’effet de ces
paramètres n’est pas toujours intuitif et est difficile à cerner. D’une part, sur résine époxy, il
constate que les éprouvettes présentant un délaminage démontrent une meilleure tenue en
fatigue, d’autre part, la résine nylon (thermoplastique), qui s’endommage moins, augmente
également la tenue en fatigue de ses éprouvettes. D’après [Curtis89], sur unidirectionnel de
carbone, l’utilisation d’une résine époxy plus tenace améliore la tenue en fatigue à haut niveau
de chargement, mais au contraire la diminue à bas niveau de chargement. Il explique ce
phénomène par l’apparition de mécanismes d’endommagement différents avec la résine plus
tenace. [Song98] compare le comportement en fatigue de deux satins à fibres de carbone, l’un
avec une résine époxy l’autre avec une résine PEEK. Il met en évidence un mode de ruine
différent dans les deux cas : la fissuration matricielle dans les mèches transverses, pour la résine
époxy se trouve être 5 à 10 fois plus importante en densité que pour la résine PEEK. Quant au
composite PEEK, il se dégrade en fatigue par rupture de fibres dont la densité atteint des valeurs
supérieures d’un ordre de grandeur au composite époxy. [Dyer98] a également travaillé sur
l’influence de la résine, sur la fatigue utilisant des fibres de verre et une résine polyester (PE) ou
polyuréthane-vinyl-ester (PU-VER). Cette influence est sensible uniquement à faible niveau de
chargement : la résine PU-VER plus tenace donne une durée de vie supérieure.
La modification de résine peut se faire par adjonction d’autre constituants tels que les
nanotubes de carbones ou de particules d’élastomères. [Higashino95] montre que l’ajout de
« caoutchouc nitrile » (NBR) dans la matrice d’un stratifié carbone/époxy améliore sa tenue en
fatigue. L’auteur remarque d’abord que les trajets des fissures transverses sont plus sinueux
lorsque des particules de caoutchouc jalonnent leurs parcours. Puis il argumente que la tenue en
fatigue du tissu est gouvernée par les méta-délaminages. Plus ils progressent, et plus les
directions de renforts se découplent, et donc plus les mèches longitudinales supportent de
charge et auront tendance à rompre prématurément. La résine renforcée est plus tenace en
mode I ce qui freine le développement du méta-délaminage en fatigue et explique la meilleure
tenue en fatigue avec cette matrice. [Manjunatha10] est arrivé au même résultat sur stratifié
verre/époxy. Des essais ont été conduits sur la résine seule afin d’évaluer son comportement en
quasi-statique et fatigue. Les particules de caoutchouc abaissent la résistance en quasi-statique,
mais améliore significativement la durée de vie en fatigue. Ce phénomène est attribué à la
cavitation des particules caoutchouteuses qui augmente la ténacité de la résine. Le même effet se
répercute sur le stratifié, avec une meilleure tenue en fatigue malgré une résistance quasi-
statique plus faible.
Enfin l’influence du motif élémentaire a été étudié par [Nishikawa06] qui compare la
réponse en fatigue de plusieurs stratifiés taffetas carbone/époxy avec des largeurs de mèches
s’échelonnant entre 0,05 et 20 mm (mèches étalées). Ces dernières se rapprochent de
l’agencement unidirectionnel, tant la longueur d’ondulation est faible (< 1 mm) par rapport à la
longueur où les fibres sont rectilignes (~20 mm). L’auteur étudie également l’influence de
32
Chapitre II – Comportement mécanique et endommagement
l’épaisseur des mèches sur l’apparition de fissuration matricielle. Les mèches transverses fines
retardent la fissuration matricielle.
33
d’une baisse de rigidité obtenue par différents mécanismes de fissuration matricielle tendant à
découpler les directions de renforts. Cet endommagement est dans sa nature, sa localisation et
sa cinétique, fonction du matériau et du niveau de chargement. Il finit par atteindre un état de
saturation avant la rupture finale qui apparait brutalement au bout d’un nombre de cycles
déterminé par la tenue des fibres. D’après les différentes sources, cet endommagement a un
effet, tantôt bénéfique et tantôt néfaste, sur la durée de vie selon le matériau et particulièrement
sa matrice.
Figure II.17 – Comportement du tissu en traction dans l’axe des renforts. (a) :
comportement sous chargement monotone, (b) : évolution de l’endommagement, (c) : évolution de
la déformation anélastique [Bizeul09].
34
Chapitre II – Comportement mécanique et endommagement
Ces essais ont été réalisés selon la norme EN-ISO14129 à déplacement imposé à la
vitesse de 0,5mm/min. Une différence est observée selon que la traction s’effectue dans l’axe des
mèches de chaîne ou de trame. Le module d’Young est légèrement plus faible dans la direction
trame (Tableau II.2), et sur la courbe contrainte/déformation l’inflexion typique des tissus de
verre (knee-point) a lieu pour une déformation inférieure (Figure II.17a).
Tableau II.2 – Résultats de caractérisation du matériau en traction [Bizeul09].
D’autres auteurs observent cette différence de l’ordre de quelques pourcents sur tissu
équilibré, dont [VanPaepegem03] et [Kumagai04], qui l’expliquent par l’ondulation plus
importante des mèches de la direction trame par rapport à la direction chaîne. En traction dans
la direction trame un état de contrainte local différent du cas de traction chaîne entrainerait un
endommagement précoce dans les mèches transverses, et une rupture prématurée des fibres. En
comparant le module d’Young E(ε) à une déformation ε donnée par rapport au module d’Young
initial E0, on définit l’endommagement longitudinal d1 comme la perte de rigidité par rapport à la
rigidité initiale de l’éprouvette :
(II.3)
35
Figure II.18 – Géométrie des éprouvettes de traction utilisées
36
Chapitre II – Comportement mécanique et endommagement
Figure II.19 – Clichés MEB d’échantillons d’éprouvettes de traction ayant subi différentes
déformations maximales ε : A gauche du schéma de l’éprouvette, les coupes « planes », à droite les
coupes de « profil » de normale x, en bas, la courbe de comportement en traction de l’éprouvette.
Les points de couleurs correspondent aux déformations maximales atteintes par les éprouvettes
sur lesquelles les clichés ont été réalisés.
37
Il est difficile d’affirmer que les points de tissages favorisent l’apparition des fissures,
étant donné que les premières fissures n’ont pas forcément été observées aux points de tissage
sur les coupes de profil, et que les coupes planes ne révèlent pas de fissures qui seraient limitées
à une zone d’ondulation. Au contraire, sur les coupes planes, on n’observe pas de fissure
transverse au niveau des points de tissage. Rappelons cependant que l’ondulation de la mèche
ne permet pas de suivre son plan médian, et que les micro-fissures observées se poursuivent
peut-être au niveau des points de tissage situés sous le plan de coupe. Cependant, les fissures
observées sur les échantillons « plan » s’étendent systématiquement sur une longueur
importante, pratiquement entre deux points de tissages. On peut donc penser que la propagation
d’une micro-fissure est instable dans une certaine mesure, et qu’une fois initiée dans une mèche
transverse, elles se propagent instantanément sur son épaisseur, mais également dans sa
longueur.
Des essais de traction ont été réalisés sur des éprouvettes de dimensions données
Figure II.18. Ils ont été filmés à une fréquence de 50Hz par une caméra FLIR SC7000 MW d’une
résolution maximale de 320 x 256 pixels. La fenêtre de visualisation permet d’observer
l’éprouvette sur une zone d’environ 30 x 30 mm². Au champ de température mesuré par la
caméra, on soustrait le champ de température initial pour obtenir le champ de température
relative, ie la variation ΔT de température au cours de l’essai. Lors d’un essai de traction
homogène sans phénomènes dissipatifs, on peut montrer que la température diminue (resp.
augmente) avec l’augmentation (resp. la diminution) de la déformation [Désoyer10].
L’apparition de phénomènes dissipatifs de type écrouissage ou pseudo-écrouissage,
viscoélasticité ou fissuration, engendre une dissipation d’énergie (plastique, par frottement, ou
de création de surfaces) par création de chaleur qui se traduit par une élévation de température.
Ces sources de chaleur ont été déterminées quantitativement à partir de l’équation de la chaleur
et sous certaines hypothèses [Lisle12] :
38
Chapitre II – Comportement mécanique et endommagement
Figure II.20 – Résultats d’un essai de traction filmé par caméra infrarouge sur la zone
encadrée de l’éprouvette en bas à gauche. (a)-(d), vue de l’évolution de la température par rapport
à la température initiale, ΔT, à différents instants. (e), évolution temporelle de ΔT moyennée sur 4
zones différentes définies par des rectangles de couleur en (a).
Le traitement analytique pour obtenir les sources de chaleur est détaillé par [Lisle12], et
nous nous contenterons d’en montrer le résultat en tant que cumul des sources de chaleur, c'est-
à-dire la quantification des phénomènes dissipatifs sommés dans le temps, en Figure II.21. On
retrouve ici, comme sur les images en température, l’apparition de phénomènes dissipatifs
localisés, qui apparaissent sur l’ensemble de la surface observée. Des agrandissements sur
lesquels on a superposé un motif de tissage, permettent de constater que les dimensions de ces
zones d’émissions de chaleur s’accordent avec la dimension des mèches. Ces sources sont
locales, étendues selon la longueur des mèches, ce qui corrobore les conclusions des
39
observations MEB. Elles se multiplient au cours du temps et tendent à recouvrir l’ensemble de la
zone observée (fissures dans toutes les mèches transverses), jusqu’à la rupture finale.
Figure II.21 – Evolution des sources de chaleur cumulées dans le temps (a-d) et
agrandissement de ces champs (e-h) sur lesquels le motif de l’armure satin de 8 a été superposé.
[Link]. Bilan
Le tissu verre/époxy étudié apparait se comporter conformément aux matériaux
similaires, avec également ses propres spécificités. L’étude des propriétés macroscopiques, les
différentes coupes micrographiques ainsi que les images de thermographie infrarouge
permettent de confirmer le scénario d’endommagement suivant : le comportement est linéaire
élastique jusqu’au « knee-point » à partir duquel des fissures matricielles apparaissent peu à peu
dans les mèches transverses, s’accompagnant d’une diminution progressive de la rigidité du
stratifié. Ces micro-fissures ne s’initient pas forcément dans les zones d’ondulation, mais sont
liées au motif de tissage. Elles traversent la mèche dans son épaisseur, et se propagent
vraisemblablement instantanément dans le sens longitudinal entre deux points de tissage. Les
mèches se fissurent les unes après les autres, sur l’ensemble de l’éprouvette de sorte qu’en fin
d’essai chaque mèche transverse est fissurée en un ou deux endroits. Très peu de ces fissures se
propagent à l’interface entre les deux directions de renforts (méta-délaminages) contrairement
à d’autres matériaux. On peut penser que la fissuration matricielle engendre un léger dépliement
des fibres qui, combiné à la présence de débris dans les fissures pourrait expliquer l’apparition
de faibles déformations inélastiques.
40
Chapitre II – Comportement mécanique et endommagement
cyclique sur éprouvettes saines ont donc été réalisés dans le but d’étudier les mécanismes
d’endommagement en fatigue, la manière dont ils se traduisent macroscopiquement, et
d’alimenter le modèle numérique en données quantitatives d’endommagement du matériau et
de tenue en fatigue.
[Link]. Principe
Ces essais ont été réalisés en accord avec la norme [ASTM D3479] pour des déformations
imposées (procédure B). La géométrie des éprouvettes est donnée Figure II.18, les drapages
[0/90]2 et [90/0]2 permettent de distinguer les directions « chaîne » et « trame ». Un chargement
sinusoïdal a été appliqué par une machine hydraulique INSTRON 8500 pour des déformations
imposées εmax s’échelonnant entre 3.10-3 et 12.10-3 mm/mm. Un rapport de charge (R=εmin/εmax)
unique a été adopté pour ces essais dans le but d’en limiter le nombre. Sa valeur a été fixée à 1/3,
en cherchant à se rapprocher du chargement typique d’une pale en le reproduisant par un signal
d’amplitude constante. La fréquence d’essai est de 10 ou 20Hz selon le chargement et les limites
techniques de la machine et de son asservissement. Les performances en fatigue des composites
à base de tissus de verres sont connus pour être dépendantes du taux de déformation et par
conséquent de la fréquence d’essai ([Curtis89] qui cite d’autres auteurs). Ce dernier
recommande d’utiliser un taux de déformation constant. Sans aller à une telle précision, une
fréquence décroissante en fonction de la déformation a été appliquée dans les essais de fatigue.
Deux ou trois éprouvettes ont été testées pour chaque niveau de déformation et pour
chaque configuration, « chaîne » ou « trame ». La rupture a systématiquement eu lieu dans la
partie courante mais au voisinage immédiat du talon. La température a été mesurée par un
thermocouple fixé sur la surface de l’éprouvette pour quelques essais préliminaires. Celle-ci
n’augmente que de quelques degrés, et nous supposerons que cela n’affecte pas les propriétés
mécaniques de la résine.
(II.4)
41
élevés cette troisième phase n’apparait pas. Rappelons également que les études en déformation
imposée sont minoritaires dans la littérature. En contraintes imposées, la déformation appliquée
augmente en cours d’essai, alors qu’en déformation imposée la contrainte apparente diminue
progressivement. En fin d’essai, notamment lorsque l’endommagement se localise, il peut
s’accélérer contrairement au cas de déformation imposée.
42
Chapitre II – Comportement mécanique et endommagement
43
Figure II.22 – Evolution du module d’Young en fatigue sur éprouvettes sens chaîne (C) et
sens trame (T). Résultats d’essais pour différentes déformations maximales εmax imposées.
Seuls les essais à εmax = 3.10-3 n’ont pas été menés à rupture, l’un d’entre eux atteignant
14.106 cycles. Pour les autres niveaux de déformation, les résultats de ces évolutions
d’endommagement sont synthétisés en fonction du nombre de cycles normé par le nombre de
cycles à rupture en Figure II.23 pour les éprouvettes sens chaîne. Les éprouvettes les plus
faiblement chargées (6.10-3) présentent une courbure légèrement plus douce que les autres,
mais dans l’ensemble, les cinétiques d’endommagement relativement au nombre de cycles à
rupture sont relativement similaires d’un essai à l’autre. Mis à part l’éprouvette C1, le module
d’Young du matériau tend vers une valeur d’environ 80% du module d’Young initial avant sa
rupture. Cette observation s’accorde avec les conclusions de [Echtermeyer95] (voir [Link]) qui
met en évidence des endommagements critiques de 23% et 27% pour des tissus de verre à
matrice phénolique et polyester respectivement.
44
Chapitre II – Comportement mécanique et endommagement
Figure II.23 – Evolution du module d’Young en fonction du nombre de cycles normé par le
nombre de cycles à rupture N/Nf pour des essais de fatigue dans le sens chaîne à différents niveaux
de déformation.
On peut ici réitérer les analyses précédentes. A bas niveau de déformation (3.10-3),
l’endommagement sature à une valeur de 5 à 10% (selon la direction chaîne ou trame qui n’est
pas discriminée sur le graphe). Au-delà d’une certaine déformation située entre 3.10-3 et 6.10-3,
l’endommagement de fatigue maximal est supérieur à 20%. Il existe un seuil d’endommagement
en fatigue à partir duquel il est possible de conduire le matériau à son état d’endommagement
caractéristique, autour de 20 à 23%.
45
Figure II.24 – Endommagement d1f mesuré en fonction de la déformation imposée en
fatigue pour différentes valeurs du nombre de cycles. L’endommagement à rupture de chaque
essai est également reporté, le pointillé noir en représente la moyenne.
46
Chapitre II – Comportement mécanique et endommagement
Figure II.25 – Durée de vie en fatigue à déformation imposée pour un rapport de charge R
de 1/3.
Ces résultats ont d’abord été approchés par une loi de type Basquin, relation entre une
contrainte de fatigue à un nombre de cycles par une loi puissance (Eq. (II.2)). Adapté à la fatigue
à déformation imposée, elle relie la déformation maximale εmax au nombre de cycles N (en
Mégacycles):
ε
ε (II.5)
(II.6)
- ε0 = 24,4.10-3
- k = 0,128
47
présentées Figure II.25 s’avèrent conservatives. [Daggumati13] a été confronté à la même
difficulté avec des géométries d’éprouvettes rectangulaires. L’auteur a alors dupliqué certains
essais sur des éprouvettes de type haltère, et démontre que ce changement de géométrie n’a pas
d’incidence sur la cinétique d’endommagement mesurée ([Daggumati13], [DeBaere11]). Par
contre il ne donne pas d’information concernant l’effet de ce changement de géométrie sur la
durée de vie.
Figure II.26 – Clichés MEB d’une éprouvette de traction chargée en fatigue à une
déformation maximale de 8.10-3 jusqu’à rupture.
Cette figure révèle les mèches transverses à la direction de la charge. Chacune d’elle est
traversée par une ou deux fissures qui s’étendent entre deux points de tissages que l’on peut
deviner. Elles peuvent être localisées au centre de la mèche comme en périphérie. Leur position
48
Chapitre II – Comportement mécanique et endommagement
dans la largeur de la mèche n’est donc pas gouvernée par les concentrations de contrainte dues à
la forme (elliptique) de la mèche car elle induirait la même position de la fissure dans sa largeur.
La position est sans doute pilotée par des effets plus locaux, les zones riches en fibres ou en
résine réparties aléatoirement dans la mèche, qui affectent le champ de déformation local et
déclenchent l’initiation de fissures.
Sur la Figure II.26c on aperçoit les fibres sectionnées par le plan de coupe au niveau de
leur ondulation. Elles se prolongent hors du plan de coupe, et on voit la fissure se prolonger
également le long de ces fibres, au moins sur quelques centaines de µm. On ne peut pas, à ce
stade, conclure si l’extrémité de fissure que l’on aperçoit Figure II.26d constitue la pointe réelle
de fissure, ou si elle se poursuit hors du plan de coupe. Au vu de la fissure visible un peu plus
loin dans son prolongement Figure II.26a on penche plus facilement pour la deuxième
hypothèse. On peut constater que dans cette zone riche en résine (Figure II.26d), l’aspect de
cette fissure est relativement accidenté. Il est cependant difficile de comparer les faciès des
fissures en quasi-statique et en fatigue à ce niveau, mais globalement, les images MEB obtenues
sont très similaires (voir Figure II.19).
[Link]. Bilan
Ces essais de fatigue, réalisés sur une large plage de déformation imposée ont mis en
évidence la différence chaîne/trame à bas niveau de déformation qui s’estompe pour des
sollicitations plus importantes. Un seuil d’endommagement en fatigue qui approche le knee-
point (seuil d’endommagement en quasi-statique) a été révélé. En dessous de ce seuil,
l’endommagement est contenu à un niveau faible, et au-delà, il converge vers une valeur
comprise entre 20 et 23%, d’autant plus rapidement que la déformation imposée est élevée. A
déformation très élevée, la rupture de fibres peut prendre le pas sur l’endommagement, dans le
sens où les fibres rompent en fatigue avant que la résine n’ait eu le temps d’atteindre son niveau
maximal d’endommagement. L’état microscopique du matériau après fatigue est proche de celui
que l’on observe en quasi-statique : chaque mèche est traversée par une ou deux fissures. Enfin,
la direction de renfort (chaîne ou trame), n’affecte pas le nombre de cycles à rupture, et
l’ensemble des résultats de durée de vie est fidèlement représentée par une loi de type Wöhler.
49
En effet, sous hypothèse d’homogénéité des contraintes, un effort imposé selon la
direction y (Figure II.27a) correspond à une contrainte axiale homogène dans l’éprouvette σy.
Les tenseurs de contrainte σ et de déformation ε s’écrivent alors :
(II.7)
(II.8)
Figure II.27 – (a) illustration des repères structure (x,y) et matériau (1,2) sur un essai de
traction sur drapage [±45]. (b), illustration de la rigidification par rotation des fibres à haute
déformation lors d’essais de traction sur stratifié [±45]2s pour différentes éprouvettes
Carbone/PEEK [Kellas92] (L’auteur explique la dispersion à déformation élevée par le procédé de
fabrication).
50
Chapitre II – Comportement mécanique et endommagement
Les propriétés de cisaillement peuvent être mesurées via cet essai qui est largement
exploité dans cette optique (Par exemple [Kellas92] [Payan04] [Marin09]). Il a été normalisé
pour l’étude des composites à matrice polymère [ASTMD3518]. La contrainte de cisaillement τ12
et la distorsion angulaire γ12 dans le repère d’orthotropie s’écrivent :
(II.9)
(II.10)
où F est l’effort appliqué et S la section de l’éprouvette.
Un nombre important d’autres essais de cisaillement plan ont été employés pour
caractériser les composites. [Lessard97] en a dénombré plus de 30, et ils constituent le sujet de
plusieurs revues bibliographiques ([Tarnopol’skii00], [Lee86]). La Figure II.28 en représente
quelques exemples. Certains sont utilisés pour déterminer uniquement la rigidité sans la
résistance (Plate Twist Test), d’autres uniquement la résistance (Four-Point Shear test, Double
Notched Shear test). Certains essais exigent des moyens d’essais spécifiques (torsion de tube),
ou des éprouvettes complexes et onéreuses. La diversité de ces essais laisse à penser qu’aucun
d’entre eux n’est applicable à toutes les études de cisaillement plan. Parmi les difficultés
expérimentales inhérentes à l’étude du cisaillement plan, on peut citer :
- L’obtention d’un état de déformation homogène qui n’apparait que sur une partie
limitée de l’éprouvette (par exemple Double Notched Shear test).
- La nécessité d’explorer des niveaux de déformation plus élevés que pour l’étude de la
traction dans l’axe des renforts, ce qui peut occasionner des non-linéarités
géométriques.
- La détermination des valeurs de contrainte et de déformation n’est pas directe et est
tributaire d’hypothèses sur la géométrie de ces champs et sur les coefficients
élastiques [Odegard11].
- Le flambement est généralement rencontré en cisaillement, incitant à l’utilisation
d’éprouvettes épaisses [DeBaere08a].
- La rupture peut apparaitre en dehors des parties courantes (Par exemple, essai
ARCAN [ElHajjar04]).
51
Figure II.28 – Exemples d’essais permettant l’étude du cisaillement plan des stratifiés
[Tarnopol’skii00].
52
Chapitre II – Comportement mécanique et endommagement
Figure II.29 – Comportement pour un essai de traction cyclé à 45° d’un tissu 2,5D
carbone/époxy. (a), courbe contrainte/déformation, (b), détail de l’évolution des déformations
résiduelles, (c), recouvrement des déformations résiduelles au cours du temps [Couégnat08].
La résistance en contrainte est bien entendu plus faible que lorsque la charge est alignée
avec un axe de renfort, mais la déformation à rupture est beaucoup plus élevée. Sur UD on peut
obtenir une rupture finale sans rupture de fibre, mais par délaminage, ce que l’agencement tissé
empêche.
53
40% en appliquant la compression transversale. En fatigue, en adimensionnant la contrainte par
cette contrainte à rupture en quasi-statique, les résultats « uniaxial load » et « biaxial load » sont
similaires.
Figure II.30 –Synthèse des durées de vie en fatigue sur des tissus à base de fibres de verre
pour des drapages [±45]n [Bizeul09]. Une loi de type Basquin a été identifiée sur les essais de
[Pandita01].
La dispersion expérimentale apparait plus importante que dans le cas de traction axiale.
[Degallaix02] indique qu’elle peut s’étaler sur deux décades de nombre de cycles à un niveau de
chargement donné. [Dyer98] qui compare deux résines pour deux drapages, [0/90]2s et [±45]4,
montre que l’influence de la résine sur les performances en fatigue est plus importante pour le
drapage désorienté.
54
Chapitre II – Comportement mécanique et endommagement
Figure II.31 - Perte de module d’un stratifié de tissu verre/époxy [±45] 5 sollicité en traction
ondulée à un niveau maximum de 50% de la contrainte de rupture et à différentes fréquences
[Pandita01].
[Karahan11] adopte la même approche, sur sergé carbone/époxy. Il détecte les premiers
évènements acoustiques à partir de 0,88% en traction hors-axe, contre 0,21% en traction dans
l’axe. Mais à partir de ce seuil, des fissures sont observées à la fois dans les mèches de chaîne et
de trame à partir de 1,2%, puis rapidement entre les deux directions de renforts (méta-
délaminage) et entre les plis. Les méta-délaminages ou délaminages s’expliquent par les deux
directions de renforts qui pivotent ce qui induit du cisaillement hors plan dans l’interface entre
les deux directions de mèches.
55
En plus d’étudier le matériau en cisaillement, l’objectif est d’obtenir un montage capable
de générer la propagation d’une coupure en mode II, en sollicitant un stratifié entaillé en
cisaillement plan. Ainsi il doit satisfaire les spécifications suivantes :
Un essai de cisaillement au rail modifié a été retenu pour cette étude. La Figure II.32
présente un schéma de principe de la version initiale telle que décrit dans la norme [ASTM
D4255]. Une éprouvette rectangulaire plane est encastrée à 3 rails par serrage. Les deux rails
externes constituent un cadre fixe tandis que le rail central est en liaison glissière avec le bâti. Le
chargement est dans la majorité des cas appliqué par compression sur le rail central. Cet essai ne
nécessite pas d’équipement particulièrement complexe, les éprouvettes sont relativement
simples à préparer (forme simple, pas de talons) et la zone est suffisamment étendue pour
autoriser une mesure de déformation, et insérer une fissure.
56
Chapitre II – Comportement mécanique et endommagement
L’essai a par la suite été normalisé en plusieurs temps jusqu’à la version [ASTM D4255]
de 2001, qui vérifie les spécifications 1 et 3 citées plus haut. Des études analytiques
([Whitney71] [Garcia80]) et numériques ([Lee86] [Lakshminarayana84] [DeBaere08b]
[MohseniShakib09]) démontrent un état de cisaillement pur, homogène, sur une zone étendue et
des effets de bords limités dans la plupart des cas. [Lakshminarayana84] a reporté des effets de
bords sur les 4 bords des rectangles sollicités qui s’étendent sur environ 10% de la longueur des
cotés. Ces effets de bords se manifestent par une diminution du cisaillement plan au voisinage
des bords libres, accompagnée de contraintes de traction et de compression parasites localisées
au niveau des coins (Figure II.33b), et qui peuvent initier une rupture prématurée de
l’échantillon le long des mors.
Figure II.33 – Simulation EF des champs de cisaillement τ12 (a) et de contrainte σ11 (b) dans
un des deux rectangles sollicités en cisaillement [DeBaere08a].
57
[Link]. Essais préliminaires
La géométrie des éprouvettes de cisaillement est reportée Figure II.34. Les parties
hachurées symbolisent les zones encastrées par serrages aux trois rails. Des essais préliminaires
ont d’abord été réalisés sur des éprouvettes [0/90]2 sur un montage de cisaillement au rail
standard tel que représenté Figure II.32. L’effort est transmis par une liaison ponctuelle entre un
plan et une bille dont le sommet est à l’aplomb du plan moyen de l’éprouvette. Un chargement
monotone quasi-statique a été appliqué, et la déformation a été mesurée à la fois par des
rosettes de trois jauges et par stéréo-corrélation d’images. Les résultats ont montré l’apparition
de flambement à des déformations de cisaillement faibles (de l’ordre de 1%), ainsi qu’un
déplacement hors plan d’environ 1 mm causé par un jeu dans la liaison glissière du montage.
Le flambement sur cisaillement au rail a été mis en évidence par [MohseniShakib09]. Son
étude numérique indique que les modes de flambement ainsi que les contraintes critiques
associées dépendent fortement des conditions aux limites appliquées. Ainsi, si le déplacement
hors plan du rail central est libre, la contrainte critique de flambement est divisée par 10. En
conséquence pour satisfaire le point 4, il semble indispensable de guider précisément le rail
central. Pour ce faire, un appui fixe a été ajouté au cadre pour guider le rail central en sa partie
inférieure. Par ailleurs, le nombre de plis des éprouvettes a été doublé, toujours pour retarder le
flambement.
58
Chapitre II – Comportement mécanique et endommagement
Figure II.35 – Données expérimentales des essais préliminaires sur montage non modifié :
(a) charge en fonction de la distorsion angulaire pour le premier et le dernier cycle de chargement,
(b) signaux des premiers cycles, et signal écrêté en fin d’essai.
Dans le montage [MohseniShakib09], le rail central coulisse à travers un cadre selon une
liaison pivot glissante double (Figure II.36b). Cette solution présente l’avantage d’être plus
facilement adaptable sur une machine de traction/compression standard.
59
Figure II.36 – Solution technologique proposée par [DeBaere08b] (a), et
[MohseniShakib09] (b) de montage de cisaillement au rail autorisant la fatigue à rapport de charge
négatif.
Afin de tirer parti du montage existant, la solution que nous avons retenue est
légèrement différente. Elle utilise une liaison ponctuelle bilatérale qui autorise les 3 rotations de
l’espace, les 2 translations du plan xz (Figure II.37), et permet de transmettre un effort dans les
deux sens de la direction y. Le schéma CAO du montage entier et du détail de cette liaison sont
montrés Figure II.37. Celle-ci est réalisée par une bille maintenue entre deux plans parallèles.
Les deux plans sont fixes l’un par rapport à l’autre après serrage de la pièce 2 sur la pièce 1 par
deux vis de fixation non représentées. La pièce faisant office de bille a été usinée dans une masse
d’acier sous forme de deux portions de sphère sur un parallélépipède percé de deux trous de
fixations. Ces trous, oblongs, permettent un réglage du sommet de la bille à l’aplomb du plan
médian de l’éprouvette dont la position dépend de l’épaisseur de celle-ci. Le positionnement est
entaché d’erreurs de quelques dixièmes de millimètre, mais le transfert d’effort dans le plan de
l’éprouvette est réalisé par le mors central dont le déplacement est guidé précisément en deux
endroits ce qui contrarie au maximum les déplacements hors plan de l’éprouvette. Ces zones de
guidages sont recouvertes par un antiadhésif de manière à limiter au maximum les efforts de
liaison.
60
Chapitre II – Comportement mécanique et endommagement
Figure II.37 – Schéma CAO du montage de cisaillement au rail modifié. (a) : montage entier,
le dispositif d’encastrement de l’éprouvette n’est pas représenté, (b) : vue en coupe du détail de la
liaison ponctuelle bidirectionnelle.
Le cadre est bridé sur un vérin hydraulique tandis que la partie supérieure de la liaison
ponctuelle est encastrée sur le plateau supérieur de la machine d’essai (Figure II.37a). La
Figure II.38 montre trois vues du montage réel au cours d’un essai. La géométrie des éprouvettes
est imposée par les dimensions du montage et n’a pas évoluée par rapport au montage standard
(Figure II.34).
61
Figure II.38 – Clichés du montage installé sur une machine de traction/compression. (a) :
vue de face, (b) : vue de la face arrière montrant l’extensomètre et une bride, (c) : détail de la
liaison ponctuelle bidirectionnelle.
[Link]. Instrumentation
Pour satisfaire l’exigence 5, les deux zones de guidage ont été recouvertes d’antiadhésif
pour empêcher les liaisons de transmettre un chargement axial, de manière à ce que la cellule
d’effort de la machine de traction/compression mesure l’effort qui transite réellement dans
l’éprouvette. La contrainte de cisaillement τ se déduit alors de l’effort mesurée F :
(II.11)
(II.12)
62
Chapitre II – Comportement mécanique et endommagement
2.3.3. Résultats
Le montage a été fixé sur une machine hydraulique de traction/compression de marque
SCHENCK équipée d’une cellule d’effort de 450kN. Elle peut être pilotée sur un des trois canaux
suivants : effort, déplacement, déformation. Les tests ont été réalisés sur éprouvettes [0/90]4,
d’abord pour valider le montage puis pour caractériser le comportement en cisaillement du
matériau.
63
verticalement (Figure II.41). Ces résultats mettent en évidence l’importance de mesurer la
déformation dans la partie centrale des zones sollicitées.
Figure II.40 – Résultats de stéréo-corrélation pour une éprouvette sous chargement quasi-
statique à une déformation homogène de cisaillement γxy = 0,0125. (a) : déplacement hors plan de
l’éprouvette, (b) (c) (d) : champs de déformations γxy, εyy et εxx respectivement.
Le second effet de bord, les concentrations de contraintes aux coins des zones sollicitées,
induit une initiation de la rupture de l’éprouvette (voir Figure II.41b) comme indiqué par
d’autres auteurs (par ex. [Lessard95]). Par conséquent, comme la rupture est localisée le long
des mors pour des contraintes de traction/compression, cet essai ne fournit pas de valeur de
résistance en cisaillement plan.
64
Chapitre II – Comportement mécanique et endommagement
La courbe de traction sur drapage [±45]2 est légèrement différente, ce qui peut
s’expliquer par la présence de traction dans les axes des renforts qui s’ajoute à l’état de
cisaillement. La rigidification par réorientation des fibres est bien observée à partir de γ=0,08, ce
qui ne permet pas la mesure d’une contrainte ou d’une déformation à rupture en cisaillement
par cet essai non plus.
65
Figure II.42 – (a), résultats d’essai de cisaillement au rail et de traction à ±45° quasi-
statique monotone. (b), pour l’essai de cisaillement au rail, la rupture est initiée dans les coins des
zones sollicitées, sièges de contraintes de traction/compression parasites.
66
Chapitre II – Comportement mécanique et endommagement
Une seconde séquence, constituée d’une charge jusqu’à la déformation γi suivie d’une
décharge, et d’une charge jusqu’à une valeur de déformation opposée -γi a été appliquée sur une
éprouvette (Figure II.43b). Deux points particuliers ont été mis en évidence sur le graphe
contrainte/déformation : les branches de charge vers les déformations croissantes sont
concourantes au point (0,0075 mm/mm ; 22±1,5 MPa) et les branches de charge vers les
déformations décroissantes se coupent au point (-0,008 mm/mm ;-24±1 MPa). Le
comportement du polyamide PA66 étudié par [Blès02] présente la même propriété sans qu’une
explication n’ait été avancée. Ces essais avaient surtout pour vocation d’identifier une loi de
comportement du matériau en cisaillement en vue d’une implémentation dans un modèle
numérique.
(II.13)
67
où G peut prendre la valeur tangentielle ou sécante du module de cisaillement. En
modélisation afin de reproduire le comportement de cisaillement de manière simplifiée,
l’hystérésis sera négligé et les boucles seront modélisées par des droites dirigées par le module
sécant [Payan04]. De plus, la mesure du module sécant est plus directe et plus robuste que celle
du module tangent, elle est donc couramment adoptée pour quantifier l’endommagement
([Totry09]).
Tableau II.3 – Module de cisaillement G0 obtenu expérimentalement.
E1 3550
E2 4040
E3 3860
Moyenne 3820
Ecart type/moyenne 6%
La déformation résiduelle, ou déformation anélastique, est notée γp. Elle été estimée
comme l’intersection de la tangente à l’origine des différentes branches avec l’axe horizontal,
c'est-à-dire après la relaxation de l’éprouvette à contrainte nulle. Les évolutions des
déformations résiduelles γP et des modules GT et GS normalisés par rapport au module initial G0,
sont tracées Figure II.44 pour deux éprouvettes différentes. En accord avec la représentation
Figure II.43a, la décroissance est plus importante avec le critère de module sécant (jusqu’à 60%)
qu’en module tangent (20%). Les résultats en module sécant sont également moins dispersés. La
déformation anélastique γp s’approxime très bien par une évolution quadratique en fonction de
la déformation totale.
68
Chapitre II – Comportement mécanique et endommagement
Figure II.44 – Evolution des grandeurs GT, GS, et γp en fonction de la déformation maximale
de cisaillement appliquée pour deux éprouvettes « 1 » et « 2 ».
69
Figure II.45 – Evolution du module de cisaillement tangent en fonction du nombre de cycles
pour une distorsion maximale γmax=0,012. La ligne en tiret épais matérialise l’endommagement
obtenu en sollicitation quasi-statique c'est-à-dire après 1 cycle.
Par ailleurs, les allures des courbes sont différentes d’un essai à l’autre. La tendance est
soit une décroissance progressive (courbes A3 et B2) similaire aux courbes de fatigue en
traction, soit au contraire une chute brutale, suivie d’une augmentation progressive de la rigidité
sur environ 105 cycles (A4 et B1). Ces essais n’ont pas été poursuivis jusqu’à rupture, car comme
sous sollicitation monotone, les concentrations de contraintes au bord des mors ne permettent
pas d’obtenir une rupture en fatigue de cisaillement.
70
Chapitre II – Comportement mécanique et endommagement
Figure II.46 – Visualisation de quelques cycles de fatigue à déformation imposée sur une
éprouvette avant sa visualisation au MEB.
Figure II.47 – Clichés MEB d’une éprouvette de cisaillement après 2.103 cycles à une
distorsion maximale de 0,024 puis 2.103 cycles à une distorsion maximale de 0,036.
Après cet essai, l’éprouvette a été démontée et découpée pour observation au MEB (voir
procédé en [Link]). Une coupe dans le plan des mèches et une coupe « profil », perpendiculaire à
la direction chaîne ont été observées (Figure II.47).
71
Malgré la déformation imposée importante, les déformations anélastiques atteintes et la
diminution importante de rigidité, le matériau ne présente pas ou très peu d’endommagement
apparent à cette échelle. La fissuration matricielle, généralement bien révélée par MEB
n’apparait pas sur les vues de la coupe dans le plan des mèches (Figure II.47c-d). Sur la coupe
« profil », normale au plan des mèches, on aperçoit de la fissuration matricielle (Figure II.47a),
mais confinée, et qui ne traverse généralement pas l’épaisseur de la mèche. Des décollements
fibre/matrice isolés sont également observés (Figure II.47b). Compte tenu de la déformation
importante, ces quelques fissures isolées semblent insuffisantes pour justifier l’altération du
comportement macroscopique, la diminution de rigidité, et les déformations anélastiques de
cisaillement. Ces observations contrastent avec le comportement en traction où la
microfissuration apparait pour des déformations bien inférieures (cf Figure II.19).
(II.14)
Figure II.48 – Schématisation du cisaillement d’un échantillon de tissu dans le repère (x,y)
se traduisant en traction/compression dans un repère tourné de 45°.
Sous l’hypothèse (très forte) de déformation homogène, la résine résiste donc à des
déformations répétées de traction/compression de 18.10-3 sans fissurer, alors qu’en traction
homogène elle se fissure à partir d’une déformation homogène de 6 à 8.10-3.
72
Chapitre II – Comportement mécanique et endommagement
de contraintes modérées. Un calcul par éléments finis à l’échelle microscopique estimant les
champs de contraintes dans la résine pourrait confirmer ou infirmer cette hypothèse.
73
pseudo-écrouissage. Ce dernier dissipe de l’énergie sous forme de chaleur. L’augmentation de
température est homogène malgré des effets de bords à proximité des rails. On ne détecte pas
d’éruptions localisées de température, perceptibles en traction, et ce même en appliquant des
échelles de température plus sensibles à ces images. La rupture de l’éprouvette s’initie aux
alentours de 70 s, en dehors de la zone de visualisation. La température n’est pas affectée par
cette rupture. Elle commence à décroitre une fois que l’essai est interrompu.
2.3.4. Discussion
Dans le cadre d’une étude de propagation de coupure pour différentes sollicitations
planes, l’étude du comportement du stratifié en cisaillement s’est révélée une problématique
plus délicate que la traction. Elle a conduit au développement d’un essai de cisaillement au rail
innovant, capable d’appliquer des chargements réversibles, et pilotable en déformation grâce à
l’adaptation d’un extensomètre à la mesure du cisaillement. L’utilisation du montage pour des
essais de propagation de coupure est exposée dans la partie suivante.
74
Chapitre II – Comportement mécanique et endommagement
75
- Le type de chargement est variable au cours d’un cycle, il s’agit alors de fatigue
multiaxiale non proportionnelle (la fatigue multiaxiale proportionnelle
([Quaresimin10]), ne rentrant pas dans la catégorie des chargements variables).
La première catégorie d’essais a pour but de formuler des lois calculant la durée de vie
d’un stratifié, ou d’évaluer l’influence de l’ordre dans lequel le chargement est appliqué. Le
critère le plus couramment utilisé, en raison de sa simplicité, est la loi de Miner. Elle introduit
une variable d’endommagement sur la durée de vie appelée fraction de vie (de Miner), notée Dm.
(II.15)
où ni est le nombre de cycles effectués au chargement Si pour lequel la durée de vie est Ni.
Ce critère prévoit la rupture lorsque Dm atteint 1.
Pour les composites cette loi est reconnue inappropriée ([Kawada05]). [Post08] indique
même qu’elle peut être non conservative. Une de ses limitations concerne la non prise en compte
de l’effet de séquence de chargement, c'est-à-dire qu’un bloc B1 de N1 cycles à une sollicitation S1,
suivi d’un bloc B2 de N2 cycles à une sollicitation S2>S1 est équivalent à la séquence « B2 suivie de
B1 ». Pour les composites cette équivalence n’est en général pas vérifiée.
La seconde catégorie d’essais est bien représentée pour les matériaux métalliques
(notamment [Calloch97]) mais a très peu été étudiée sur composites. Dans une large revue
bibliographique d’essais et de modèles de fatigue multiaxiale sur composites, [Quaresimin10]
cite le déphasage entre deux directions de chargement comme un paramètre important de la
tenue en fatigue, mais ne donne que très peu de résultats. On peut également citer l’étude de
[DeMonte10] sur tubes en fibres courtes verre/polyamide qui semble montrer une influence
faible du déphasage entre les sollicitations de traction et de torsion. Il s’avère que le cas le plus
pénalisant, c'est-à-dire le plus conservatif en fatigue multiaxiale est celui d’un déphasage nul, ie
la fatigue multiaxiale proportionnelle.
Les essais présentés ici appartiennent à une troisième catégorie, très rarement étudiée,
celle du couplage entre 2 endommagements consécutifs obtenus par deux types de sollicitations
différentes. La seule étude recensée qui traite d’essais successifs est celle de [Yaniv89], bien que
l’approche soit différente. Il a développé une méthode de mesure de l’endommagement en
fatigue de traction par un essai de cisaillement. En remarquant que le module de cisaillement est
sensible à la fissuration transverse et longitudinale, il propose d’estimer l’endommagement
longitudinal de ses éprouvettes de traction, d1, à partir d’une mesure du module de cisaillement.
Pour cela il soumet des éprouvettes à de la fatigue de traction, puis les adapte sur un montage de
76
Chapitre II – Comportement mécanique et endommagement
L’étude de l’effet de la fatigue de traction dans une direction (ε1) sur le comportement en
traction dans l’autre direction (ε2), puis sur le cisaillement (γ12) est d’abord présentée. Enfin le
couplage (γ12 – ε1) est évalué. Le couplage (γ12 – ε2) est équivalent à (γ12 – ε1) si l’on ne tient pas
compte de la distinction chaîne/trame. Le Tableau II.4 résume cette démarche.
Tableau II.4 – Synthèse des types de couplage étudiés. La sollicitation (i) est
systématiquement cyclique (fatigue). Le type de sollicitation indiqué dans les cellules est celui de
la sollicitation (ii).
(i)
fatigue ε1 fatigue γ12
(ii) .
ε1 statique
ε2 statique et fatigue
γ12 statique
77
Figure II.51 – Géométrie de l’éprouvette de traction utilisée pour l’étude des couplages
traction-traction, et traction-cisaillement.
Les deux premières éprouvettes, E1 et E2, ont été pilotées à εmax = 6.10-3, jusqu’à une
rupture prématurée par rapport aux résultats à ce niveau de déformation sur éprouvettes
standards (Figure II.22b). L’évolution du module d’Young pour les 4 essais est reportée
Figure II.52. Pour les deux premiers essais (E1 et E2), le module diminue rapidement, sans
atteindre la deuxième phase où la progression est plus lente. La fissuration matricielle n’a donc
pas eu le temps de se développer suffisamment, avant que les fibres ne rompent en fatigue. Cette
rupture finale de l’éprouvette est apparue prématurément à un endommagement limité à 13%
seulement contre environ 20% pour les éprouvettes de taille standard au même chargement.
Pour le troisième essai, E3, un chargement à 4.10-3 a été appliqué sur 5.105 cycles puis
augmenté à 5.10-3. L’éprouvette a finalement rompu sans atteindre un endommagement
supérieur à 7%. Le chargement faible s’est révélé peu endommageant pour la résine mais très
pénalisant pour les fibres puisque le nombre de cycles effectués à 5.10-3 est très faible. On note la
rupture de pente importante qui apparait dans l’évolution du module d’Young, pour une
augmentation de la déformation de seulement 1.10-3. Il semble qu’un seuil d’endommagement en
fatigue soit ici mis en évidence.
78
Chapitre II – Comportement mécanique et endommagement
Figure II.52 – Evolution du module d’Young en fatigue pour les éprouvettes de grande taille
E1 à E4, et leur comparaison avec des éprouvettes de drapage équivalent de taille standard.
En termes de durée de vie, ces éprouvettes de grandes tailles se sont révélées peu
endurantes. Nous ne disposons pas de point de comparaison pour les éprouvettes E3 et E4 mais
les éprouvettes E1 et E2 ont été réalisées dans les mêmes conditions que C8 et C9. La différence
peut s’expliquer par l’utilisation d’une machine différente dont les mors exercent une pression
sur une partie seulement des talons. Cela se traduit par une surcontrainte néfaste pour le
matériau. D’autre part les dimensions des talons sont différentes ce qui a pu accentuer les
concentrations de contraintes.
Cela met également en lumière un effet d’échelle, que l’on peut appréhender par la
théorie statistique de Weibull sur la distribution des défauts et la probabilité de défaut d’autant
plus importante que la pièce est étendue.
79
Enfin, la flexion parasite selon l’axe 3 (normal à l’éprouvette), qui peut apparaitre au
moment du serrage des mors génère des déformations axiales supplémentaires proportionnelles
à la largeur de l’éprouvette. En éprouvette large, elles pourraient être suffisantes pour que la
déformation au bord de l’éprouvette soit supérieure au seuil d’endommagement d’un coté, qui
s’endommage alors rapidement, tandis que l’autre soit en deçà et s’endommage peu. De plus,
l’endommagement est calculé à partir de la contrainte moyennée sur la section, et de la
déformation prélevée par l’extensomètre au centre de l’éprouvette qui ne doit pas être impactée
par la flexion. L’endommagement mesuré doit donc correspondre à la moyenne des
endommagements sur la section de l’éprouvette entre les couteaux de l’extensomètre. Cela
pourrait expliquer, dans le cas d’une composante de flexion significative que les pentes
d’endommagement des éprouvettes E1 et E2 soient plus faibles que C8 et C9. Cette flexion
expliquerait également que la rupture finale, pilotée par la déformation maximale dans
l’éprouvette apparaisse plus tôt.
Toutefois, ces résultats sont à relativiser, car les essais réalisés n’ont pas été reproduits.
Néanmoins l’objectif de cette campagne, à savoir endommager notablement le matériau afin
d’étudier son influence sur une autre sollicitation a été atteint pour l’éprouvette E4.
Suite à ces essais, les éprouvettes E2 et E4, ainsi qu’une éprouvette saine E5, ont été
découpées à la scie à disque diamantée (voir Figure II.53) en éprouvettes secondaires de
traction dans le sens 2 et de cisaillement. Des talons ont été collés sur les éprouvettes de traction
sens 2 (couplage (ε1 – ε2)) (voir Figure II.54a), et les éprouvettes de couplage (ε1 - γ12) ont été
percées (Figure II.54b).
80
Chapitre II – Comportement mécanique et endommagement
Figure II.54 – Géométrie des éprouvettes de couplage (ε1 – ε2) (a) et (ε1 – γ12) (b).
(II.16)
81
Tableau II.5 – résultats d’essais de couplage (ε1 - ε2) quasi-statiques.
82
Chapitre II – Comportement mécanique et endommagement
l’endommagement préalable d1 sur les résultats de fatigue dans la direction 2. De manière plus
détaillée, à déformation modérée, les éprouvettes E4-4, E4-5, et E4-6 ont tout de même tendance
à s’endommager un peu plus rapidement que les éprouvettes saines. Cette observation est à
relativiser par rapport à l’incertitude, mais serait cohérente avec l’idée que les fissures déjà
existantes dans une direction favorisent la concentration de contraintes, et se comportent
comme des amorces de fissures pour l’autre direction.
83
Figure II.56 - Evolution du module d’Young tangent en fatigue, comparaison entre
éprouvettes saines (préfixes E5) et préalablement endommagées dans le sens 1 à hauteur de 19%
(préfixe E4) pour différents niveaux de déformation maximale.
84
Chapitre II – Comportement mécanique et endommagement
Il apparait sur la Figure II.57 que la réponse de l’éprouvette E5-1 concorde avec celle de
la géométrie standard jusqu’à une distorsion d’environ 2%. Au-delà, le changement de géométrie
de l’éprouvette semble avoir un effet. Des phénomènes d’effets de bords qui s’intensifient avec la
déformation pourraient expliquer ces écarts, puisque relativement à la taille de l’éprouvette,
ceux-ci sont plus étendus sur l’éprouvette de petites dimensions. L’endommagement d1 semble
avoir une influence sur le comportement en cisaillement. Plus d1 est important, plus la rigidité de
l’éprouvette diminue. Cet effet est moins marqué sur la raideur initiale de cisaillement que sur le
pseudo-écrouissage apparent.
85
Certaines modélisations sur l’endommagement des stratifiés UD [Ladeveze90] où à
renforts tissés lient l’endommagement transverse à l’endommagement de cisaillement par une
relation de proportionnalité ([Miot09] [Hochard10]):
(II.17)
[Miot09] décompose le tissu en la superposition de 2 plis d’unidirectionnels et a
identifié le coefficient c à 1 pour les plis unidirectionnels équivalents par un essai de traction sur
un tissu déséquilibré [±22°]. Cette même valeur est ensuite utilisée sur deux autres tissus
équilibrés. Elle suppose alors que la fissuration matricielle a quantitativement le même effet sur
la rigidité en traction transverse que sur celle en cisaillement.
Sur le présent matériau le couplage est moins important. Avec cette décomposition
simplifiée d’un tissu en deux nappes d’unidirectionnels, seule une nappe est endommagée, le
cisaillement dans l’autre nappe ne devrait donc pas être affecté. On devrait alors se trouver avec
86
Chapitre II – Comportement mécanique et endommagement
un coefficient de couplage de l’ordre de c = 0,5. Avec le peu de données à disposition Tableau II.7,
le tissu de verre étudié semble présenter un coefficient de couplage c < 0,3.
Figure II.58 – Illustration du découpage des éprouvettes de couplage (γ 12 - ε1) après essai
de fatigue en cisaillement.
Les modules d’Young initiaux E0 sont reportés Tableau II.8. La contrainte à rupture n’a
pas été relevée, car l’utilisation de talons en aluminium d’épaisseur trop importante a engendré
des ruptures prématurées de toutes les éprouvettes (saine ou non) à proximité immédiate, ou en
dessous des talons.
Les valeurs obtenues pour le module d’Young sont toutes en accord avec la valeur
nominale du matériau sans que l’endommagement préalable n’ait d’effet. Le couplage (γ12 –ε1)
semble inexistant pour notre matériau, ce qui amènerait à la valeur de c infiniment grande dans
ce cas (voir Eq. (II.17)).
87
Figure II.59 – Comparaison des réponses contrainte/déformation en traction entre des
éprouvettes saines et préalablement sollicitées en fatigue de cisaillement.
Tableau II.8 – Comparaison des modules d’Young en traction entre éprouvettes saines et
préalablement sollicitées en fatigue de cisaillement.
88
Chapitre II – Comportement mécanique et endommagement
2.4.5. Discussion
Les essais présentés dans cette partie ont eu pour vocation de couvrir un large panel de
couplages possibles entre les différentes sollicitations planes. Toutes les configurations n’ont pas
été étudiées pour deux raisons :
Le temps que nécessite une campagne d’essai avec les différentes étapes (essai de fatigue
- découpage - second essai) a limité les couplages de type fatigue-fatigue qui ont été assez peu
abordés.
Mais surtout, pour les essais déjà réalisés, des tendances marquées ont été dégagées et
elles sont en cohérence avec l’idée que l’on peut se faire du scénario d’endommagement suite à
l’étude sous comportement de traction et de cisaillement précédemment présentée.
Un réseau de fissures orientées selon une direction n’a intuitivement pas de raison de
provoquer de baisse de rigidité dans une direction orthogonale et cette idée à été confirmée par
les essais (ε1 - ε2). Les fibres transverses (direction 2) ne sont pas affectées par la sollicitation ε1
et leur résistance en traction quasi-statique ou fatigue est inchangée. Par contre la fissuration
matricielle a une incidence sur la rigidité de cisaillement comme classiquement supposé dans les
modèles d’endommagement. Elle dégrade les transferts des efforts dans la résine qui supporte
en grande partie la sollicitation. Quantitativement cette chute de rigidité est limitée car la
fissuration ne concerne qu’une seule direction de renforts. Dans l’autre direction les transferts
d’efforts par cisaillement de la résine ne sont pas dégradés.
Par contre, pour le matériau étudié, un endommagement de cisaillement ne crée que très
peu de fissuration matricielle. L’affaiblissement du module de cisaillement résulte sans doute de
phénomènes de fissuration à l’échelle des chaînes de polymères. Tant que la résine n’est pas
fissurée, le comportement dans l’axe des renforts n’est pas altéré.
Le couplage (ε1 – γ12) est perceptible alors que le couplage (γ12 –ε1) est nul. Ce constat
nous empêche de définir un couplage biunivoque entre les endommagement d2 (ou d1) et d12 et a
fortiori une loi linéaire de type (II.17). La raison physique de cette limite réside dans la
différence de nature entre les endommagements de traction et de cisaillement.
89
Chapitre III – Etude expérimentale de la
propagation de coupure
91
3.1. Introduction
L’étude de la propagation de coupure en fatigue constitue le cœur de ce travail de thèse.
On rappelle que l’on utilisera le terme « coupure » pour désigner les ruptures du stratifié à
travers l’épaisseur, avec nécessairement rupture de fibres dans le cas des tissus. Cette
appellation permet de lever l’ambigüité avec les « fissures », qui s’accordent plus couramment
avec la rupture de matrice ou d’interface de différentes natures : microfissuration matricielle,
décohésion fibre/matrice, méta-délaminage ou délaminage. Dans ce chapitre, le stratifié entaillé
comportant entre 2 et 4 plis a été étudié pour différents drapages standards, et sous
sollicitations de traction ou de cisaillement. On présentera d'abord une revue des études
comparables recensées dans la littérature en se focalisant sur les techniques d'essais et
d'acquisition, les types de phénomènes détectés, les résultats caractéristiques observés, et de
préférence en s’orientant sur des matériaux à renforts tissés. La méthodologie des essais réalisés
et de leurs dépouillements est ensuite présentée, suivie des résultats obtenus en traction puis en
cisaillement.
Les trous et les perçages sont également des accidents géométriques couramment
rencontrés sur les composites, pour des applications évidentes d’assemblages boulonnés ou
rivetés. La fatigue sur éprouvettes percées retient l’attention de certains auteurs, et il se trouve
que le stratifié troué présente des similitudes de comportement avec le stratifié entaillé. Ainsi au
92
Chapitre III – Etude expérimentale de la propagation de coupure
vu du peu de résultats sur stratifiés entaillés, cette bibliographie a été élargie aux stratifiés qui
présentent des accidents de géométrie en général, siège de concentrations de contraintes.
93
Figure III.2 – Effet d’un trou sur la courbe S-N de stratifiés tissés et non-tissés à base de
fibres de carbone [Bishop89].
94
Chapitre III – Etude expérimentale de la propagation de coupure
Figure III.3 – (a), Courbes S-N d'éprouvettes saines et trouées. (b), C-scan d’une éprouvette
trouée après différents nombres de cycles de fatigue à 0,85σr [Wang02].
95
Figure III.4 – Résistance résiduelle en fonction du nombre de cycles de fatigue subis par
l’éprouvette [Bizeul09].
Encore une fois, l’effet bénéfique de la fatigue sur éprouvette percée est plus
spectaculaire sur unidirectionnels. [Ramani77] mesure plus de 30% d’augmentation de la
contrainte résiduelle sur éprouvette trouée de drapage [0;±45;90]2s d’unidirectionnels après
fatigue. [Wang02] a étudié l’influence du niveau de contrainte et du nombre de cycles
(Figure III.5) sur l’augmentation de la résistance résiduelle. Il cite la concurrence entre des
phénomènes de dégradation qui résultent en la redistribution des contraintes de manière plus
uniforme (appelés « wear-in »), et qui tendent donc à augmenter la durée de vie et la résistance,
et d’autre part la dégradation du matériau qui l’affaiblit, dénommé « wear-out », le conduisant à
la rupture finale. Il montre que plus le niveau de contrainte est élevé, plus la résistance
résiduelle est augmentée (voir Figure III.5), ce qui suggère qu’à niveau de contrainte élevé les
phénomènes de wear-in prennent le pas sur le wear-out. L’endommagement apparait plus tôt, et
se développe plus rapidement, ce qui relaxe rapidement les contraintes, et laisse les fibres à des
niveaux de chargement plus faibles.
96
Chapitre III – Etude expérimentale de la propagation de coupure
97
Figure III.6 – (a) Evolution de la longueur de coupure en fonction du nombre de cycles pour
6 éprouvettes identiques et pour la même sollicitation cyclique. (b) Vitesse de propagation en
fonction de la longueur de coupure pour les mêmes essais. Les éprouvettes qui initient vite
propagent vite et réciproquement [Bizeul09].
Enfin, on retrouve la même idée dans les résultats de [Kawai96] dont l’étude a été
mentionnée au Chapitre II. Il compare l’influence des résines époxy et nylon dans les
performances d’éprouvettes de tissu de carbone saines et entaillées en fatigue. La résine époxy
se fissure plus facilement que la résine nylon thermoplastique et les éprouvettes
correspondantes ont tendances à délaminer plus facilement. Ainsi, bien que l’auteur montre une
meilleure tenue en fatigue sur éprouvette saine des éprouvettes carbone/nylon, en présence
d’une entaille, un endommagement intervient rapidement en pointe de coupure pour le stratifié
carbone/époxy et induit du wear-in qui soulage les fibres. Finalement en configuration entaillée,
c’est le stratifié carbone/époxy qui présente une meilleure tenue à la propagation.
98
Chapitre III – Etude expérimentale de la propagation de coupure
nouvelle surface dA, et donc faire propager une coupure de la longueur da = dA/b ou a est la
longueur de coupure et b l’épaisseur de l’éprouvette.
(III.1)
On distingue les modes I, II et III de rupture, par GI, GII et GIII qui partitionnent le taux de
restitution d’énergie. Le mode I est considéré comme le plus critique sur les métaux et
correspond au mode le plus largement étudié. Sous chargement quasi-statique, la propagation a
lieu lorsque le taux de restitution d’énergie atteint une valeur critique Gc, qui est idéalement une
propriété matériau mais dans la réalité dépend de paramètres géométriques de l’éprouvette
(notamment drapage et longueur de coupure) et du mode de rupture. [Laffan12] précise qu’il
s’agit d’une valeur structurale, en comparant le taux de restitution d’énergie critique des fibres
seules (7 J/m²) à celles de stratifiés qui atteignent généralement des valeurs 1000 fois plus
élevées. En effet, dans le composite, la rupture de fibre s’accompagne d’une zone étendue
d’endommagement (décollement fibre/résine, fissuration matricielle).
(III.2)
(III.3)
où CI dépend des coefficients élastiques du matériau.
En fatigue, la vitesse de propagation da/dN peut être reliée à une puissance du taux de
restitution d’énergie ou d’un facteur d’intensité de contrainte selon une loi de Paris :
(III.4)
Elle peut s’écrire en facteur d’intensité de contrainte KI, ou en variation d’une des deux
grandeurs au cours d’un cycle de fatigue (ΔKI ou ΔG).
Cette loi empirique est largement utilisée en dimensionnement, s’applique très bien dans
de nombreux cas sur métallique et en propagation de délaminage, et trouve un écho important
chez les différents auteurs cherchant à identifier une loi sur leurs résultats expérimentaux.
[Shindo05] a réalisé des essais de propagation sur tissu verre/époxy à partir d’un essai
CT à température ambiante et cryogénique. Il mesure la longueur de coupure par la méthode de
99
la complaisance couplée à une méthode par éléments-finis, et le taux de restitution d’énergie est
obtenu numériquement par calcul de l’intégrale J de Rice. Les résultats obtenus sont cohérents
avec une loi de Paris (Figure III.7).
Figure III.7 – Identification d’une loi de Paris sur des essais de propagation à différentes
températures [Shindo05].
[Bizeul09] constate que cette loi est valable pour un drapage quasi-isotrope. Par contre
sur drapage [0/90] dont les courbes de propagation sont présentées Figure III.6, la dispersion
expérimentale se retrouve dans les droites de Paris tracées Figure III.8. Les taux de restitution
d’énergie indiqués par cette figure ont été calculés par une méthode numérique (VCCT), sur un
modèle qui tient compte du temps d’initiation en tant que diminution de la rigidité de
l’éprouvette. Malgré cette prise en compte, les éprouvettes suivent des lois de Paris de pentes
semblables mais décalées verticalement d’un facteur 10.
Figure III.8 – Identification d’une loi de Paris pour 6 essais sur éprouvettes [0/90] 2 réalisés
selon les mêmes conditions [Bizeul09].
100
Chapitre III – Etude expérimentale de la propagation de coupure
En marge des études sur composites, on peut citer les travaux de [Trappe12] qui
concernent une résine époxy pure. L'auteur a réalisé des éprouvettes de Single Edge Notch
modifiées (Figure III.10a) pour observer la vitesse de propagation d'une coupure pour des
échantillons d’une résine ayant subit une post-cuisson de 15h à 50° dans un cas et 60° dans
l’autre. Il montre que malgré des propriétés identiques de rigidité et ténacité en quasi-statique
pour les deux résines, la vitesse de propagation en fatigue mesurée est augmentée d'un facteur
101
10 (Figure III.10b). Sans que l’auteur n’avance d’explication, il illustre ainsi la différence de
mécanismes de rupture en statique et en fatigue sur la résine époxy.
Dans le même ordre d’idées, [Hansen99] a réalisé des essais de fatigue sur éprouvettes
7781/G913 quasi-isotropes pré-impactées à 9 ou 20 J. La rupture en quasi-statique a lieu près
du talon et non dans la zone impactée quelle que soit l’énergie d’impact. En fatigue, les
éprouvettes hautes énergies, perforées par l’impacteur à 20J, montrent au niveau de l’impact
une propagation progressive de l’endommagement jusqu’à rupture. Pour l’impact de faible
énergie (9J), présentant simplement un BVID, en fatigue la rupture a lieu près des talons à
niveau de contrainte élevée et au niveau de l’impact à niveau de contrainte faible. Finalement,
alors qu’il n’avait pas affecté le comportement en quasi-statique, l’impact devient donc critique
en chargement cyclique, et plus spécialement à chargement modéré.
Figure III.10 - (a), éprouvette de propagation de résine époxy. (b), résultat de propagation
de type loi de Paris [Trappe12].
102
Chapitre III – Etude expérimentale de la propagation de coupure
(III.5)
103
surtout de la largeur de l’éprouvette, qui est souvent assez faible afin de retarder le flambement.
Pour forcer la propagation en mode II dans des plis à 0°, [Laffan13] a intercalé de nombreux plis
à 90°, c'est-à-dire dans l’axe de la pré-entaille. Il utilise pour cela des éprouvettes de drapage
[904/0/904]2s qui atteignent 4,68 mm d’épaisseur. Il montre des faciès de rupture différents en
fonction du taux de mode II avec une zone de délaminage d’autant plus étendue que le mode II
est important.
Figure III.12 - Essai Compact Tension Shear utilisé pour la propagation translaminaire en
mode II et mixte I+II par [Laffan13].
En mode II pur, [Tan88] a utilisé des éprouvettes de cisaillement au rail sur lesquelles il a
pratiqué 2 entailles centrées. Son étude concerne les chargements quasi-statiques uniquement.
Il observe la propagation de coupures à 45° par rapport à la direction initiale de l’entaille, et
fréquemment le délaminage des plis qui conduit à du flambage local.
104
Chapitre III – Etude expérimentale de la propagation de coupure
Il conclut que l’éprouvette de type « écharde » est critique par rapport aux deux autres,
en cohérence avec les valeurs de taux de restitution d’énergie calculées par éléments finis pour
chaque configuration. La vitesse de propagation a été tracée en fonction du taux de restitution
d’énergie en mode I pour le cas « écharde », elle se révèle en bon accord avec une loi de Paris.
Des difficultés ont été éprouvées pour suivre la propagation des « fissures de bord » et « fissures
centrales », le front de coupure a été le plus difficile à suivre de par la présence de pontages de
fibres et de dissymétrie entre les deux faces de l’éprouvette. Le matériau est un unidirectionnel
verre/époxy et les géométries sont telles que la propagation ne traverse pas les fibres, mais
l’auteur précise qu’une approche similaire peut être adoptée notamment pour des tissus.
3.2.6. Bilan
Les composites présentent une bonne résistance aux accidents géométriques de par leur
capacité à s’endommager et à redistribuer les contraintes sur une zone plus étendue. Cet
endommagement peut être néfaste par exemple en compression, mais s’avère bénéfique pour la
durée de vie en traction cyclique. Celle-ci est fortement liée à la tenue en fatigue des fibres. En
vue d’une étude expérimentale on retiendra que le temps d’initiation est généralement long, et
représente la plus grande partie de l’essai. Puis la propagation se fait à une vitesse qui dépend
du taux de restitution d’énergie. L’utilisation d’une loi de Paris est fréquente mais peut se révéler
hasardeuse dans certains cas. L’état d’endommagement du composite dépend de sa nature et de
sa géométrie, et des phénomènes spécifiques à la fatigue peuvent se révéler. Il apparait que les
zones d’ondulation des mèches favorisent la rupture et guident la trajectoire de la coupure.
Enfin, la propagation en mode II apparait plus délicate, elle soulève des difficultés
expérimentales liées au flambage.
105
3.3. Méthodes expérimentales
3.3.1. Procédure commune aux essais de propagation en fatigue
Les essais de propagation en traction et en cisaillement ont été réalisés selon une
procédure commune. L’éprouvette est préalablement entaillée par une scie à fil diamanté de
diamètre 0,17mm. Un chargement de fatigue ondulé à déformation imposée est appliqué. Il est
caractérisé par la déformation maximale imposée εmax et le rapport de charge R = εmin/εmax fixé à
1/3. L’essai est divisé en séquences. Une séquence est constituée d’une rampe de chargement
quasi-statique jusqu’à la déformation εmax, d’un palier de quelques secondes, puis d’un nombre
ΔN de cycles à une fréquence f de 10 ou 20Hz selon les cas. Enfin une décharge quasi-statique
jusqu’ à effort nul termine la séquence. L’essai est arrêté lorsque la coupure s’est propagée sur
l’ensemble de l’éprouvette ou qu’elle s’est stabilisée.
Le suivi de coupure est effectué par l’acquisition d’images à l’aide d’une caméra CCD.
Celle-ci enregistre des images numériques de résolution 1280x960 px², en nuances de gris codés
sur 16 bits. Une source d’éclairage halogène placée derrière l’éprouvette révèle
l’endommagement sur le tissu de verre. On exploite ici, à l’instar de [Marissen06] ou [Shindo05],
la propriété de relative transparence des composites à base de fibres de verre ainsi que le
phénomène de blanchiment de résine qui permet l’utilisation de cette technique dite de rétro-
éclairage.
La procédure d’acquisition d’images est automatisée par une routine LabView®. Celle-ci
surveille une sortie tout-ou-rien active uniquement lors du palier à déformation maximale de
chaque séquence (Figure III.14). La sortie déclenche d’abord l’activation d’un relai alimentant
l’halogène, puis la prise de vue et son enregistrement.
106
Chapitre III – Etude expérimentale de la propagation de coupure
107
ont été déterminées suite au dépouillement des clichés de l’éprouvette au cours de l’essai. Une
routine Matlab® a été développée dans ce but. Le principe de cette routine est schématisé de
manière simplifiée Figure III.15. Elle permet d’effectuer la soustraction pixel par pixel entre
deux images « a » et « b », de façon à en révéler la différence, ou l’évolution du faciès de rupture
de l’éprouvette entre deux instants. Ensuite, une étape de seuillage et d’amplification du
contraste permet d’obtenir une image exploitable D1. Un filtrage de type Wiener est ensuite
appliqué à l’image obtenue. Il s’agit d’un filtre passe-bas utilisé sur des signaux à la fois bruités et
flous, qui présente l’avantage de préserver les gradients. Pour étudier nos résultats nous
chercherons à comparer l’image à un nombre de cycles N par rapport à l’image initiale « 0 », ou
l’évolution entre l’image N et l’image N+ΔN. La résolution de la caméra et la taille du champ
observé induisent une taille de pixel comprise entre 0,03 et 0,05 mm. Au cours de l’essai, le
champ de vision initial peut se décaler progressivement du fait de l’environnement d’essai, et
des perturbations vibratoires sur un durée de plusieurs jours. Pour améliorer la précision des
mesures, une étape de recalage des images est préalablement effectuée en cherchant le couple
(Δx, Δy) qui maximise le produit de convolution entre les deux images d’entrée. Le principe de la
routine de soustraction est finalement résumé Figure III.15 et un exemple de l’interface
utilisateur est présenté Figure III.16.
108
Chapitre III – Etude expérimentale de la propagation de coupure
109
[Bizeul09] a également démontré l’influence d’une autre bande d’unidirectionnel
carbone située en amont de la coupure et analogue à un arêtier. Sur les drapage [±45]2, on
n’observe pas ou peu de propagation avec le longeron seul. L’ajout d’un arêtier sur les drapages
[±45]2 uniquement, permet d’observer la propagation à chaque fois que ces deux bandes
d’unidirectionnel carbone supplémentaires sont présentes. Les talons sont collés par-dessus
l’ensemble stratifié verre/époxy entaillé + renforts carbone.
110
Chapitre III – Etude expérimentale de la propagation de coupure
Sur la Figure III.19, sont comparées les évolutions des longueurs de coupure au cours du
temps. La longueur considérée n’inclut pas l’entaille préalable de 7 mm. Un scénario commun à
toutes les éprouvettes peut être dégagé. Dans une première phase, la propagation est lente. Des
longueurs de coupure non nulles sont détectées avant 105 cycles, et la coupure atteint
généralement 5 mm entre 105 et 3.105 cycles. A partir de cette longueur, elle se propage
rapidement, jusqu’à atteindre une stabilisation autour de 30 mm, où elle entre dans une dernière
phase de stabilisation. A ce moment là, la coupure rejoint le longeron en carbone, dont la raideur
importante s’oppose à l’ouverture. On montre par un modèle par éléments finis linéaire que le
taux de restitution d’énergie diminue à l’approche du longeron.
111
Figure III.19 – Evolution de la longueur de coupure en propagation de traction.
Comparaison entre les éprouvettes de 2 ou 4 plis, sollicitées dans la direction chaîne ou trame.
112
Chapitre III – Etude expérimentale de la propagation de coupure
vitesses de propagation en fonction de la longueur de coupure pour les essais 4 plis et pour 2
essais 2 plis (Figure III.20).
Par ailleurs, une différence de comportement entre les éprouvettes chaîne et les
éprouvettes trame a été constatée par [Bizeul09]. Elle est perceptible pour les éprouvettes 2 plis
Figure III.19, mais n’apparait pas sur les éprouvettes 4 plis. [Bizeul09] note cette différence
particulièrement marquée dans le cas d’éprouvettes plus étroites de largeur 30 mm. Sur la
Figure III.21 sont reportés ses résultats pour 6 éprouvettes trames (ie [90]2). Il compare à ces
résultats les évolutions des éprouvettes [0]2 ayant propagé le plus rapidement, et le plus
lentement. Les éprouvettes de trame s’initient plus tard et propagent plus lentement que les
éprouvettes sens chaîne.
113
Figure III.21 – Evolution de la longueur de coupure (notée « a ») en fonction du nombre de
cycles pour des éprouvettes à mèches de trame travaillantes de largeur 30 mm [Bizeul09].
La longueur de coupure des graphes présentés notamment Figure III.19 et Figure III.20
est déterminée par technique de soustraction d’images (voir 3.3.2). Les images de la Figure III.22
comparent la progression de l’endommagement évalué par soustraction d’images entre les
cycles N1 et N2 indiqués sur la flèche chronologique. A titre d’exemple, la première image est
issue de la différence entre les images à 0 et 2500 cycles. Pour permettre une comparaison des
différents résultats, les paramètres de recalage, de seuillage, de contraste, et de filtre de Wiener
sont conservés identiques (voir Figure III.14).
114
Chapitre III – Etude expérimentale de la propagation de coupure
Une zone d’endommagement très nette apparaissant dès les premiers cycles est révélée
par la différence entre les deux premières images de l’essai. Par la suite, aucune évolution n’est
constatée jusqu’à N = 10000 cycles. Puis un évènement plus localisé est détecté en pointe de
l’entaille initiale. Si la mèche en pointe d’entaille avait rompu avant 10000 cycles, cette zone
aurait été totalement déchargée et n’aurait pas subi d’endommagement ultérieur. Ce contraste
détecté entre 10000 et 15000 cycles puis entre 15000 et 20000 cycles correspond
115
vraisemblablement à des ruptures de mèche, c'est-à-dire à l’initiation de la coupure. Un scénario
identique a été constaté sur les autres essais réalisés dans les mêmes conditions. Il rappelle les
images de [Marissen06] sur drapages quasi-isotropes de plis unidirectionnels, entaillés, et
chargés de manière monotone. Le rétro-éclairage permet donc de distinguer une zone
d’endommagement en amont de l’entaille, indépendamment des premières ruptures de fibres
qui apparaissent plus tard.
116
Chapitre III – Etude expérimentale de la propagation de coupure
117
3.4.3. Drapage [±45]2
On rappellera ici les résultats obtenus par [Bizeul09] pour des drapages [±45]2 (aucun
essai supplémentaire n’a été réalisé lors de cette thèse). Trois essais ont été réalisés pour ce
drapage sur la géométrie de la Figure III.17 avec la présence d’arêtier ; ils se sont tous les 3
propagés dans l’axe de la pré-entaille. L’initiation de la coupure est très rapide, avec un fort
développement en début d’essai, qui s’atténue progressivement (Figure III.24). La progression
semble s’interrompre lorsque la coupure atteint 25 à 28 mm, c'est-à-dire à 5 ou 8 mm de la
bande d’UD carbone faisant office de longeron. Ces courbes présentent des paliers, ou de légères
oscillations attestant du caractère discret de la propagation comme dans le cas des drapages
[0/90]4.
Par contre, sur l’ensemble de l’essai, le comportement est radicalement différent par
rapport au cas [0/90]4. Les essais sont plus longs d’un rapport 10, et les vitesses de propagation
des éprouvettes [±45]2 (reportées Figure III.25) présentent une décroissance très marquée. En
outre la dispersion apparait plus faible pour ce drapage, qui résulte sans doute d’un temps
d’initiation nul pour chaque éprouvette.
118
Chapitre III – Etude expérimentale de la propagation de coupure
Plus précisément, on peut voir sur cette figure, que les cas à 0° et 90° sont équivalents en
termes d’orientations privilégiées des points de tissage (Figure III.27a-b). Tous deux présentent
deux axes caractéristiques à 18,4° et 45° par rapport à l’horizontale. Par contre pour les cas à
+45° et -45°, les angles caractéristiques sont 0° et ±63,4° dans un cas (Figure III.27c), et ±26,6° et
90° dans l’autre (Figure III.27d). N’ayant pas conscience de cette distinction au moment du
drapage des 4 éprouvettes, et le tissu étant équilibré, les directions +45° et -45° du rouleau de
tissu n’ont pas été distinguées. Néanmoins, nous avons constaté après polissage post mortem
des plis à ±45° de chaque éprouvette, que le motif est dans la configuration 26,6° (Figure III.27d)
pour les éprouvettes [45;0]s-1 et [0;45]s-1 et 0° (Figure III.27c) pour les éprouvettes [45;0]s-2 et
[0;45]s-2, en accord avec les directions de propagation observées.
119
Figure III.26 – Cliché de 4 éprouvettes de drapage quasi-isotrope à la fin de l’essai de
propagation en traction, obtenu par rétro-éclairage. Les croix rouges représentent les directions
privilégiées des points de tissages pour les plis à +45° (c) et (d) ou -45° (a) et (b).
Figure III.27 - Motif de tissage de l'armure satin de 8 orienté dans différentes directions. En
termes de directions privilégiées par rapport à l’horizontale, les directions 0° et 90° sont
équivalentes à la symétrie près contrairement aux directions +45° et -45°.
120
Chapitre III – Etude expérimentale de la propagation de coupure
L’angle de propagation d’environ 26° sur 2 des essais semble indiquer que ce sont les plis
à ±45° qui pilotent la direction de propagation pour au moins deux essais. Pour les deux
éprouvettes qui se propagent à 0°, la même conclusion serait hasardeuse. Une direction
privilégiée des points de tissages est certes dans la direction 0° (Figure III.27c) mais cette
direction de propagation est également observée pour toutes les éprouvettes [0/90]4.
121
Figure III.29 – Evolution de la vitesse de propagation en fonction de la longueur de coupure
pour des éprouvettes quasi-isotropes et une éprouvette [90]4.
Pour conclure, en dépit du faible nombre d’essais, il apparait que la disposition du motif
de tissage des plis à ±45° par le choix de la direction +45° ou -45° dans le drapage, a une
influence sur la direction de propagation et sur sa vitesse, et que le cas Figure III.27c est le plus
pénalisant sur la durée de vie. Dans l’autre cas, les plis à ±45° semblent guider la propagation,
bien que sa vitesse correspondent clairement à celles que l’on mesure sur drapages [0/90]n. Le
lien entre la direction de drapage et la direction de propagation soulève une question de
causalité. Il est difficile de dire si la vitesse de propagation plus faible dans la configuration
Figure III.27d découle de la direction de l’agencement des points de tissages, ou de la direction
prise par la coupure qui elle-même est causée par l’agencement des points de tissages.
122
Chapitre III – Etude expérimentale de la propagation de coupure
une seule pointe de coupure. La solution que nous avons retenue, utilise donc une pré-entaille
débouchante de longueur a0, dont la position est paramétrée par le rapport b/w (Figure III.30).
Le paramètre b/w prend les valeurs 0,2 ou 0,5 selon les cas, tandis que a0 a été fixé à 25 mm, de
manière à s’exempter des effets de bords, et que la dissymétrie de l’éprouvette ne conduise par à
une déviation trop importante du rail central. A l’aide d’un calcul par éléments finis de l’essai de
cisaillement au rail, en imposant le rail central rigide mais libre dans le plan, on peut estimer,
pour une déformation de cisaillement appliquée γ = 12.10-3, que le déplacement selon x du rail
central est limité à 3 µm. La Figure III.31 représente le maillage et les conditions limites
imposées pour ce calcul. Ce modèle a permis de vérifier que l’entaille ne perturbe pas
l’homogénéité du champ de déformation dans la zone non-entaillée. Cette zone non entaillée a
donc été utilisée pour mesurer la déformation via l’extensomètre adapté au cisaillement, collé
sur cette partie de l’éprouvette (Figure III.30).
Les pré-entailles des éprouvettes ont été réalisées à l’aide d’une mini scie circulaire, et
par une scie à fil diamanté de 0,17 mm de diamètre pour les derniers millimètres.
123
Figure III.31 – Champ de distorsion angulaire γxy et conditions limites du modèle par
éléments finis de l’essai de cisaillement au rail permettant d’estimer la rotation du rail central.
124
Chapitre III – Etude expérimentale de la propagation de coupure
Figure III.32 – Résultats d’essai de propagation d'un stratifié [0/90] 4 entaillé sous
chargement de cisaillement quasi-statique. Différents états sont représentés correspondant à des
niveaux de déformation croissants. Le champ de déformation dans le repère (1,2) obtenu par
stéréo-corrélation est représenté pour l’état 3.
(III.6)
125
Le critère suppose la propagation de la coupure dans la direction θc pour lequel R(θc) est
maximum. [Beuth89] a comparé les résultats de ce critère à des résultats d’essais de
sollicitations planes et pour diverses directions de stratifiés unidirectionnels. Il définit la
résistance du pli dans la direction θ par une loi elliptique en fonction des résistances en traction
à 0°, Xt et à 90°, Yt:
(III.7)
Dans le cas d’un tissu équilibré, cette hypothèse est inappropriée puisqu’elle suppose
une résistance identique dans toutes les directions (car XT = YT). En réalité cette résistance est
minimale à 45°, la direction où les fibres sont le moins sollicitées. Sans connaitre le champ de
contrainte en pointe de coupure, le critère de ratio de contrainte normale pourrait être plus
facilement atteint pour θ = 45°.
126
Chapitre III – Etude expérimentale de la propagation de coupure
127
Figure III.35 - Evolution de l'endommagement en propagation sous chargement de
cisaillement de l’éprouvette S2 obtenue par soustraction d’images. (a), entre image initiale et 10 6
cycles, (b) entre 106 et 5.106 cycles, (c) entre 5.106 et 12.106 cycles.
Dans les deux cas, S1 et S2, la direction de propagation est guidée par le motif de tissage
(Figure III.36). La zone de propagation de l’éprouvette S2 a été prélevée et polie pour
observation MEB (Figure III.36d). Il apparait clairement que le motif de tissage impose le chemin
à la coupure. La coupure elle-même reproduit un motif. Selon une période spatiale de 3 mèches
horizontales, un point de tissage est rompu, puis la progression se fait plus ou moins le long de la
mèche verticale suivante, avant de se réorienter vers le point de tissage suivant. Si l’on considère
le sens d’application du cisaillement, en descendant avec un léger angle par rapport à l’axe
vertical, la coupure se propage quasiment perpendiculairement à la direction de traction
maximale (45°). On portera surtout l’attention sur le fait que 4 plis sont superposés mais que
leurs points de tissage n’ont aucune raison d’être alignés, et également que les positions de leurs
motifs de tissage relativement les uns aux autres varient selon l’éprouvette. Comme pour le cas
de traction, le fait de superposer 4 plis au lieu de 2 devrait avoir pour effet de moyenner cette
incertitude. Néanmoins, la forme tridimensionnelle de la coupure doit être beaucoup plus
complexe que l’image obtenue par une seule coupe peut nous montrer. On peut imaginer que
lorsque la coupure se dirige vers le bas de l’image MEB, elle est peut-être guidée par le point de
tissage d’un pli inférieur. Cette hypothèse pourrait être confirmée par des coupes plus profondes
de l’échantillon, voire l’utilisation de tomographie qui pourrait apporter la visualisation de la
coupure dans l’espace, et améliorer notre compréhension à ce sujet. Ces déviations de la coupure
que nous avons nommées « coupures secondaires » ou « subcracks » (terme employé par
[Mandell75]) traversent en fait réellement des mèches horizontales. Elles s’accompagnent de
fissuration matricielle se développant dans les mèches verticales, générant l’endommagement
révélé par les images optiques (Figure III.34). Il apparait en bas de l’image MEB qu’elles
traversent l’intervalle entre deux mèches avec un angle de 45° donc perpendiculairement à la
direction de traction maximale.
128
Chapitre III – Etude expérimentale de la propagation de coupure
Figure III.36 – Mise en évidence du trajet de coupure privilégié par les points de tissage
dans le cas de la propagation de coupure en cisaillement. (a) et (b), clichés de l’éprouvette après
propagation. (c), schéma de l’armure satin de 8. (d), vue post-mortem au MEB de l’éprouvette S2
qui met en évidence le trajet de la coupure à travers les points de tissage et les fissures
secondaires à intervalles réguliers.
129
Aucune fissure n’a pu être observée au-delà de 10 mm à la verticale de la coupure, alors
que la soustraction d’images révèle un changement de teinte étendu sur environ 20 mm. La zone
d’altération de couleur identifiée optiquement a été reportée sur la Figure III.37 en pointillés
rouges. Le blanchiment de résine est attribué au phénomène de craquelures (crazes en anglais
qui donne le crazing) ou à la cavitation des nodules de chaînes polymères [Hiver02]. Ces
phénomènes sont difficilement visibles aux grossissements employés ici, mais ces images
permettent de conclure que la fissuration matricielle principale ne suffit pas à expliquer le
blanchiment de résine observé.
Finalement un dernier groupement d’images MEB (Figure III.38) nous montre un faciès
de rupture complexe, très endommagé, par de la micro-fissuration. La coupure traverse les
mèches, par l’intermédiaire de subcracks qui se propagent simultanément. Elles évoluent
quasiment longitudinalement le long des mèches et les traversent par endroit. Elles rejoignent la
mèche suivante en se propageant à 45° (Figure III.38). Dans ces zones inter-mèches riches en
résine, que l’on peut considérer comme isotropes, la direction de propagation selon le critère de
contrainte normale, σθθ, maximale semble s’appliquer.
130
Chapitre III – Etude expérimentale de la propagation de coupure
Un autre essai a été réalisé dans la configuration suivante : 12.106 cycles à γmax = 6.10-3
puis la déformation a été augmentée à 12.10-3, jusqu’à rupture. Sur la Figure III.39, on peut
constater que l’endommagement s’est très peu développé pendant la phase à déformation faible,
puis il s’étend de manière importante dès lors que la déformation a été augmentée à 12.10-3. Le
phénomène de « subcrack » apparait alors, ce qui semble montrer que le chemin de la coupure
est influencé par le niveau de sollicitation. Pour γmax = 6.10-3, la zone d’endommagement visible
est très faible (Figure III.39a) et la coupure se propage en ligne droite. A γmax = 12.10-3, la
coupure emprunte un chemin accidenté analogue à celui de la Figure III.36, et la zone
d’endommagement est beaucoup plus importante (Figure III.39b).
131
Figure III.39 – Evolution de l’endommagement au cours de l’essai S3 : déformation imposée
de γmax = 6.10-3 pendant 12 Mcy puis γmax = 12.10-3 jusqu’à rupture.
La coupure S3 présente une croissance régulière sur 12.106 cycles (γmax = 6.10-3), à un
niveau inférieur à 10-6 mm/cy, évidemment bien plus faible que les éprouvettes S1 et S2. En
augmentant la sollicitation à γmax = 12.10-3, la vitesse de propagation retrouve les valeurs des
éprouvettes S1 et S2, et conserve la même tendance à la baisse à l’approche du mors.
132
Chapitre III – Etude expérimentale de la propagation de coupure
133
montre par calcul numérique que celui-ci augmente, et pour rester dans le domaine des
propagations à grand nombre de cycles, la déformation imposée a été ramenée à γmax = 10.10-3.
Figure III.42 –Résultat de l’essai S4. (a) cliché de la coupure en fin d’essai. (b)-(e), résultats
de soustraction d’images entre les différents instants indiqués. La trace rouge marque la position
de la pré-entaille.
134
Chapitre III – Etude expérimentale de la propagation de coupure
135
entravées par les plis à ±45° et s’étendent beaucoup moins. La concentration de contrainte en
bord de trou se relaxe peu en fatigue et les éprouvettes sont alors moins endurantes.
On peut citer également l’observation de [Marissen06] qui réalise des essais de traction
sur éprouvettes quasi-isotropes entaillées. Il note l’effet « bridging » qu’ont les plis à 45° sur les
plis à 0° (dans l’axe du chargement) et qui atténuent les concentrations de contrainte.
Les vitesses de propagation des deux essais quasi-isotropes et de l’essai S4 réalisés dans
les mêmes conditions ont été reportées Figure III.45. L’allure générale de décroissance de la
vitesse de propagation avec l’avancée de la coupure est retrouvée. Les éprouvettes quasi-
isotropes propagent plus lentement que l’éprouvette orientée.
136
Chapitre III – Etude expérimentale de la propagation de coupure
137
pour le cas [0;45]s ce qui d’ordinaire, se traduit par une propagation plus rapide. Cependant, une
déformation imposée ne se traduit pas par un champ de déplacement identique pour tous les
drapages [Whitney71] et cela peut se traduire par un état de déformation différent dans les
mèches à 0° et 90° selon qu’elles se situent dans un drapage [0/90]4 ou [0;45]s.
Bien que le nombre d’essais soit réduit, l’orientation du pli extérieur au drapage semble,
comme dans le cas de la traction, n’avoir aucune incidence sur la vitesse de propagation.
3.6. Synthèse
La propagation de coupure a été étudiée sur des éprouvettes de 2 ou 4 plis de tissu
verre/époxy pour différents drapages, géométries d’éprouvettes et en sollicitation de traction et
de cisaillement. Le phénomène de propagation présente beaucoup de dispersion dans certains
cas (traction largeur 30 mm) et très peu dans d’autres (notamment traction sur [±45]). Les cas
de propagation dispersés semblent être ceux présentant un nombre de cycles d’initiation élevé.
Le phénomène de rupture en fatigue est lié à la microstructure du matériau de par
l’endommagement, la fissuration matricielle inhérente à la taille de la mèche, et la manière dont
la coupure se propage. La propagation est fortement influencée voire guidée par le motif de
tissage. En traction, on a constaté que la coupure s’arrête dans la zone entre deux mèches et la
traverse brutalement. En cisaillement le chemin de la coupure est plus complexe. On note la
cohabitation entre de l’endommagement sous forme de fissuration matricielle et de
l’endommagement sous forme de pseudo-plastification et, sans doute, de micro-crazing. Cet
endommagement est marqué par l’influence entre les plis des différentes orientations,
notamment pour le drapage quasi-isotrope, qui réduit considérablement l’endommagement
visible par blanchiment de résine.
138
Chapitre IV – Modélisation numérique de
la propagation de coupure
139
Le coût important d’une campagne d’essais de propagation, ainsi que le temps nécessaire
à leur réalisation incite à une réduction du nombre d’essais. En parallèle, dans le domaine
aéronautique, les autorités de certification exigent de justifier la tolérance aux dommages dans
des situations plus nombreuses qu’auparavant et dont certaines sont trop complexes pour être
représentées fidèlement par un essai. Combiné à cela, l’expansion relativement récente des
outils numériques ainsi que le progrès des performances de calcul, le développement de
modélisations numériques suscite une demande croissante dans le milieu industriel. Concernant
les matériaux composites, de nombreuses stratégies de modélisation voient le jour, pour
appréhender les différents problèmes de dimensionnement. Elles sont fonction de la
problématique étudiée (dégradation, tenue en fatigue, rupture multiaxiale, endommagement,
fatigue, impact) et du volume à représenter s’étalant de l’échantillon élémentaire à la structure.
Dans un premier temps (4.1), nous nous intéresserons aux approches employées pour
représenter la rupture dans les composites présentant un défaut géométrique, en portant une
attention particulière pour les stratifiés à renforts tissés. La dépendance de la rupture à la zone
d’endommagement nous conduit à nous intéresser également à la représentation continue ou
discrète de l’endommagement dans le stratifié.
Dans les parties suivantes, nous proposerons un modèle développé pour répondre à la
problématique de propagation de coupure en fatigue, adapté au problème, au matériau
composite tissé, et aux observations expérimentales. Son principe est décrit dans la partie 4.2. Il
met en jeu plusieurs types d’éléments qui sont détaillés dans les parties 4.3 à 4.6. Enfin, la
partie 4.7 donne un aperçu synthétique des aspects pratiques de l’algorithme utilisé. La méthode
a ensuite été implémentée dans le logiciel commercial SAMCEF© et l’ensemble des résultats
obtenus font l’objet du Chapitre V.
4.1. Bibliographie
4.1.1. Mécanique linéaire de la rupture et méthodes analytiques
[Link]. Eléments de mécanique de la rupture
Les premières approches de représentation de la rupture abordent des cas simples de
matériau isotrope à comportement élastique linéaire. Le taux de restitution d’énergie est
classiquement défini, (équation (III.1)) comme la variation d’énergie potentielle dEp lors d’une
avancée infinitésimale de fissure da. L’idée de cette notion est attribuée à Griffith en 1920 qui
relie cette grandeur à l’énergie de création de nouvelle surface dA = bda, (avec b la largeur du
front de fissure) c'est-à-dire l’énergie nécessaire à la rupture des liaisons atomiques sur dA.
140
Chapitre IV – Modélisation numérique de la propagation de coupure
Pour les décrire, on utilise classiquement un repère cylindrique centré sur la pointe de
fissure et dirigé par son axe (Figure IV.2). Les travaux d’Irwin ont en particulier mis en évidence,
que quel que soit le chargement, les champs de contrainte et de déformation présentent une
singularité en r-1/2, et a ainsi démontré l’insuffisance d’un critère de résistance traditionnel pour
la détermination de la rupture.
141
(IV.1)
Par calcul, on recense de nombreuses méthodes pour déterminer les facteurs d’intensité
de contrainte et/ou de taux de restitution d’énergie que l’on peut classer en 4 grandes familles :
Méthodes analytiques :
Les méthodes analytiques découlent des développements asymptotiques effectués par
différents auteurs pour de nombreux cas simplifiés de géométrie de pièce fissurée, et de
chargement. Ils font l’objet de recueils (par exemple ceux de [Tada73] et de [Murakami87]), et
sont fréquemment considérés comme des cas de référence pour évaluer les méthodes
numériques, bien qu’ils traitent de cas idéalisés en faisant appel à des hypothèses de dimensions
infinies, ou de chargements ponctuels. Pour des chargements plus complexes, on a recours au
théorème de superposition. La somme de plusieurs sollicitations conduit à la somme des
facteurs d’intensité de contrainte via l’hypothèse de linéarité de comportement. Les facteurs
d’intensité de contrainte peuvent ensuite être reliés au taux de restitution d’énergie par :
(IV.2)
Cette méthode s’étend aux matériaux orthotropes qui introduisent des termes
supplémentaires dans les expressions des facteurs d’intensité de contrainte. [Sih65] donne les
expressions des taux de restitution d’énergie pour un matériau orthotrope homogène :
142
Chapitre IV – Modélisation numérique de la propagation de coupure
(IV.3)
Approche locale :
L’approche locale consiste en la résolution numérique du problème par éléments finis,
qui donne un champ de contrainte en pointe d’entaille sur lequel le facteur d’intensité de
contrainte peut être identifié par l’équation (IV.1). La singularité est bien entendu mal
représentée avec des éléments classiques ce qui a mené à l’introduction d’éléments destinés à
représenter spécifiquement les singularités par des points d’intégration judicieusement placés
(par exemple éléments de Barsoum). Dans ce cas, ou si le maillage utilisé est suffisamment fin,
on peut approcher le champ de contrainte par une fonction en r-1/2 et ainsi déterminer un facteur
d’intensité de contrainte par l’équation (IV.1).
Différence d’énergie :
D’autres méthodes que l’on trouve sous le nom de VCE pour « Virtual Crack Extension »
ou VCCT « Virtual Crack Closure Technique » (voir [Krueger02]) reposent sur une « différence
d’énergie » entre deux longueurs de fissure et sont implémentées dans les logiciels par éléments
finis standards. Le principe est de calculer l’énergie potentielle Ep stockée dans une structure
chargée présentant une entaille de longueur a, puis dans la même structure présentant une
longueur a+Δa (correspondant à une nouvelle surface de fissure ΔS). L’avancée Δa est obtenue
soit en dupliquant un nœud et en relâchant un effort entre les deux nœuds (VCCT) ou par
introduction d’une perturbation des nœuds du front de fissure (VCE). On obtient ainsi une
estimation du taux de restitution d’énergie :
(IV.4)
En outre, cette méthode fournit une partition des modes de rupture en déplaçant les
nœuds de la propagation virtuelle selon l’une des trois directions. Cette méthode requiert un
maillage fin en pointe d’entaille, et dans la pratique on recommande de faire plusieurs calculs
pour différentes extensions. Pour appliquer cette méthode en orthotrope, la direction de
143
l’entaille et sa direction de propagation doivent être confondues avec un axe de symétrie
élastique du matériau [Schmitt02].
Figure IV.3 – Illustration de la méthode d’avancée des fissures par déplacement des nœuds
au voisinage immédiat de la fissure [François93].
Intégrales de contour :
Enfin, d’autres méthodes utilisent une intégrale de contour sur un chemin s’appuyant sur
les lèvres de la fissure. La plus répandue, l’intégrale J (égale au taux de restitution d’énergie G),
est attribuée généralement à [Rice68] et parmi les variantes d’écriture on la notera :
(IV.5)
où 1 est la direction de la fissure, We, la densité d’énergie élastique (σij) le tenseur des
contraintes, (ui) le vecteur déplacement, et (ni) la normale au contour s’appuyant sur les lèvres
et dirigée vers l’extérieur (Figure IV.4).
144
Chapitre IV – Modélisation numérique de la propagation de coupure
D’autres critères ont été ensuite été avancés. [Gope10] cite parmi eux :
(IV.6)
Ce critère a été évalué avec succès par [Beuth89] en traction sur des éprouvettes
d’unidirectionnels orientés à 0°, 15°, 75° et 90° par rapport à la direction de la charge, et sur des
145
essais Iosipescu à 0°, 15°, 30° et 45°. Dans tous les cas la propagation s’est produite le long des
fibres, donc perpendiculairement à la direction la plus faible du point de vue de la résistance.
Ces deux critères empiriques sont relativement limités et s’adaptent mal au tissu. Dans
ce cas, on a Ex1 = Ex2 et Xt = Yt donc pas d’orthotropie visible pour ces critères qui considèrent les
propriétés selon les axes de symétrie élastiques. En particulier la traction dans les directions
±45° sollicite la matrice, de résistance moindre, ce qui n’est pas pris en compte par ces deux
critères.
146
Chapitre IV – Modélisation numérique de la propagation de coupure
Figure IV.5 – Comparaison entre les prévisions de rupture et les données expérimentales
en traction et en cisaillement pour différentes longueurs d’entaille et pour un drapage [0/90/±45] s
(a) et [0/90]2s (b) d’après [Tan88].
[Miot10] a également utilisé le même critère non local. Après une identification de la
longueur caractéristique à partir de deux configurations d’éprouvettes trouées, les résultats
obtenus sont cohérents avec les essais pour différentes géométries d’entailles (éprouvettes
doublement entaillées) et de drapages.
147
Figure IV.6 - Illustration de la zone identifiée par [Hochard06] pour le critère non-local,
utilisé pour déterminer la rupture en bord de défaut.
Ces critères sont simples à mettre en œuvre mais leur inconvénient majeur est leur
dépendance aux paramètres d’essai. Ainsi, [Chang84] a réalisé des essais de traction et de
compression sur éprouvettes trouées et a identifié des distances caractéristiques différentes
pour le PSC, en traction et compression (d’un rapport 4), différentes également selon les
drapages et selon les diamètres de trous pour un même matériau.
148
Chapitre IV – Modélisation numérique de la propagation de coupure
Dans un premier temps, différents auteurs ont cherché à prendre en compte la géométrie
des tissus en représentant l’ondulation des mèches ainsi que leur forme (souvent elliptique)
dans le calcul de propriétés élastiques (par exemple [Kriz85] ou [Naik92]). [Scida99] a pris en
compte l’angle des mèches avec la direction du pli afin de déduire les coefficients élastiques par
un modèle analytique se basant sur une loi des mélanges adaptée aux tissus. Ce modèle a été
évalué sur de nombreux tissus présentant plusieurs armures de tissage. [Takeda12] a également
pris en compte la section des mèches et leurs parcours dans un maillage à l’échelle
microscopique pour déterminer les coefficients élastiques à différentes températures.
[Searles01] a utilisé une échelle plus fine en modélisant un volume élémentaire d’un tronçon de
fibre entourée de résine. Les coefficients élastiques du tissu ont été déduits à partir des
géométries de mèches, des formes de section et des données matériau des constituants.
Au-delà des données élastiques, il est possible de rendre compte des sur-contraintes
induites dans les zones d’ondulation ou entre les mèches. [Melro12] a réalisé un maillage
microscopique d'un motif élémentaire de satin de 5 composé de mèches entourées de résine
(Figure IV.8). A partir de coefficients élastiques des constituants élémentaires (fibres de carbone
et résine époxy) l'auteur a dans un premier temps calculé les coefficients élastiques et les
résistances selon les différentes sollicitations pour une mèche en prenant en compte la
dispersion dans la répartition des fibres. Puis, à l’échelle supérieure, les mèches et la résine sont
maillées séparément. Les éléments des mèches sont dotés d’un comportement élastique
endommageable dans le sens transverse et un comportement élasto-plastique endommageable
en cisaillement est attribué à la résine. Ce modèle fournit des champs de contrainte et de
déformation locaux permettant d'obtenir l'évolution de la raideur et de l'endommagement sous
sollicitation de traction, cisaillement, ou traction-cisaillement.
149
Figure IV.8 – (a) VER d’une mèche tenant compte de la dispersion des fibres, portion de
maillage à l’échelle du motif de tissage, (c) et (d) champ des variables d’endommagement
transverse du composite en traction à deux niveaux de charge [Melro12].
Ensuite, l’architecture de deux tissus 2,5D est reproduite par un maillage représentatif
d’une mèche (Figure IV.9b). L’auteur introduit une modélisation discrète de l’endommagement
(décohésion entre les constituants ou fissuration) dans les cellules élémentaires par séparation
des éléments en doublant les nœuds traversés par les fissures. L’effet des différents types
d’endommagement à l’échelle microscopique est ensuite estimé sur les propriétés élastiques du
fil, puis du pli. En présupposant la forme des fissures, la cinétique d’endommagement est
obtenue par identification inverse. Cette cinétique d’endommagement inclut également le critère
de rupture des fibres.
Le modèle a été testé sur plusieurs drapages et certaines structures, dont un essai de
traction sur une éprouvette présentant 2 entailles décalées (Figure IV.9c-d).
150
Chapitre IV – Modélisation numérique de la propagation de coupure
151
Le modèle a été comparé uniquement à un résultat de traction simple et représente bien
la baisse de rigidité du matériau.
Figure IV.10 – (a), maillage fin d’une cellule élémentaire de mèche, (b), modèle binaire
d’éléments poutres superposés à un élément résine [Römelt12].
Cette modélisation a été adaptée aux tissus (par ex. [Payan02], [Hochard06]).
[Hochard06] modélise le stratifié à l’échelle du pli, et fait l’hypothèse que la rupture du stratifié
apparait pour la rupture du premier pli. Cette hypothèse est à la fois conservative et adaptée au
tissu. Trois variables d’endommagement d1, d2 et d12 sont définies pour chaque pli pour
caractériser les raideurs E1, E2 et E12 par rapport aux raideurs initiales E10, E20 et E120.
(IV.7)
La fraction de fibres dans la direction 0° est notée δ et chaque pli de tissu est considéré
comme la superposition de deux plis d’unidirectionnels fictifs à 0° et à 90° pondérés par δ et
affectés d’une variable d’endommagement transverse (d2*0° et d2*90°) et de cisaillement (d12*0° et
d12*90°). Généralement, les plis UD ne présentent pas d’endommagement dans la direction
longitudinale autre que la rupture totale. Les variables d’endommagement s’écrivent alors en
fonction du module transverse de l’UD équivalent E20*:
152
Chapitre IV – Modélisation numérique de la propagation de coupure
(IV.8)
En contraintes planes l’énergie élastique dans chaque pli UD fictif est alors écrite
classiquement :
(IV.9)
(IV.10)
(IV.11)
avec as, bs, l, m, Y0, c, des constantes matériau à identifier sur des essais homogènes. Ce
modèle a ensuite été généralisé aux sollicitations de fatigue ([Hochard10]) en ajoutant à d2UD*
définie précédemment la variable df déduite de l’expression :
(IV.12)
La rupture est ensuite pilotée par des critères non-locaux sur les forces
thermodynamiques des plis comme précisé en 4.1.2.
153
composites à matrice céramique et des tissus 3D ([Maire98], [Marcin10]). Ces modèles visent à
reproduire le comportement du matériau et son endommagement en prenant en compte les
déformations thermiques ainsi que la viscoélasticité. [Rakotoarisoa11] a élargi le modèle
d’endommagement aux sollicitations de fatigue. Ses variables d’endommagement comprennent
une part statique et un cumul d’endommagement en fatigue selon une loi analogue à l’équation
(IV.12).
Figure IV.11 – Faciès de rupture d’une éprouvette de traction [45; 90; -45; 0]s [Wisnom10].
154
Chapitre IV – Modélisation numérique de la propagation de coupure
Les éléments cohésifs ont également (mais moins fréquemment) été employés pour
modéliser la fissuration matricielle et l’endommagement dans le pli. Pour modéliser un essai de
traction trouée, la démarche de [Hallett09] a été de relever la position des fissures en essai puis
de mailler l’éprouvette en plaçant des éléments cohésifs sur les chemins potentiels de
propagation de la fissure (Figure IV.13). Ces éléments s’ajoutent aux éléments cohésifs de
délaminage. Cette technique lui permet d’obtenir un scénario d’endommagement en très bon
accord avec les essais. L’inconvénient réside dans la dépendance au maillage, en particulier du
placement des éléments cohésifs.
155
Figure IV.14 – Illustration des travaux de [McLaughlin00] selon 2 situations de délaminage
et fissuration transverse combinés.
Par ailleurs ce type de modélisation a également été appliqué à l’étude de l’impact par
[Bouvet09] qui utilise des éléments d'interface interlaminaires représentant le délaminage et
intralaminaires représentant la fissuration matricielle.
156
Chapitre IV – Modélisation numérique de la propagation de coupure
Figure IV.15 – Illustration des enrichissements des degrés de liberté d’après [Moës01] sur
un maillage régulier (a), et quelconque (b).
L’ensemble des nœuds est noté I. La fissure est représentée numériquement via un
enrichissement par des degrés de liberté associés à des fonctions de formes discontinues (sur
l’ensemble des nœuds repassés en bleu noté J) pour les éléments traversés par la fissure, et
asymptotiques (sur l’ensemble des nœuds repassés en rouge noté K) pour décrire les champs de
contrainte et de déformation singuliers en pointe de fissure.
(IV.13)
Par suite, le déplacement d’un point x s’exprime :
(IV.14)
Les fonctions Fl des nœuds de K sont définies comme combinaisons linéaires des
fonctions asymptotiques en r1/2 de la mécanique linéaire élastique de la rupture :
(IV.15)
157
bj et ckl constituent les degrés de libertés supplémentaires correspondant aux 2 types
d’enrichissement (discontinu et asymptotique).
Cette méthode s’est rapidement popularisée et a progressivement été intégrée dans les
principaux logiciels de calcul commerciaux. Les travaux concernant cette méthode s’appliquent
majoritairement aux matériaux métalliques. Coté composite, [Hettich08] propose une méthode à
une échelle microscopique pour mailler une résine renforcée par des inclusions. La Figure IV.16
montre le maillage utilisé. Il représente une résine et des inclusions. La propriété
d’indépendance au maillage des X-FEM conduit à utiliser des mailles de taille constante.
Certaines mailles représentent alors les 2 constituants à la fois. Pour ce faire L’auteur utilise une
formulation d’ensembles de niveaux, c'est-à-dire une fonction indicatrice qui permet à l’intérieur
d’un même élément d’attribuer à chaque point x les matériaux de l’un ou l’autre des constituants
(résine ou inclusion). Un volume élémentaire d’une inclusion entourée de résine est maillé. Les
interfaces inclusions/résine sont considérées et un processus multi-échelle est ensuite envisagé
pour considérer des structures.
Figure IV.16 – Maillage d’un matériau à plusieurs phases (matrice renforcée par des
inclusions) en X-FEM et utilisation de la formulation d’ensembles de niveaux pour définir les
géométries des constituants [Hettich08].
158
Chapitre IV – Modélisation numérique de la propagation de coupure
Figure IV.17 - Illustration de la nécessité d’un nombre important de points de Gauss par
éléments appartenant au support des nœuds j [Asadpoure06].
Une méthode moins répandue, notée E-FEM [Oliver96] repose sur un enrichissement des
éléments et non des nœuds, comparée à la méthode X-FEM [Oliver06], sur des matériaux
isotropes. Ces méthodes X-FEM et E-FEM ont très peu été appliquées aux matériaux composites
ou alors avec une représentation rudimentaire du matériau. [Hettich08] et [Wang12] ont
considéré les deux matériaux homogènes séparément. [Asadpoure06] a réalisé des simulations
sur matériaux orthotropes, mais en se restreignant à un comportement linéaire non
endommageable. Ces méthodes numériques sont confrontées aux difficultés spécifiques au
composite, notamment dans la définition d’un critère de rupture et d’un critère de direction de
propagation.
(IV.16)
159
Figure IV.18 – Notations utilisées et grille de particules [Askari06].
(IV.17)
Lorsque les deux particules sont éloignées d’une distance supérieure à l’horizon, elles
n’interagissent pas (f = 0). Sinon l’effort d’interaction entre les deux particules s’apparente à une
élasticité c jusqu’à un allongement maximal s0. La fonction f s’écrit alors :
(IV.18)
160
Chapitre IV – Modélisation numérique de la propagation de coupure
pour la modélisation de la fissure car les ruptures de liens sont traitées comme les autres
déformations, et sans nécessiter de loi de propagation ou de critère de direction.
[Askari06] a appliqué la méthode sur composite. Des liens de rigidités c différentes ont
été définis selon les positions relatives des particules, qu’elles soient à 0°, 90° ou 45°, ou entre
deux plis. La Figure IV.20 montre des résultats de simulation pour une éprouvette de traction
fissurée pour différents drapages. Le modèle est capable de rendre des faciès différents selon
que la rupture soit dominée par la rupture matricielle, le délaminage ou la rupture de fibre.
Figure IV.20 – Faciès de rupture obtenus par modélisation péridynamique d’un essai de
traction présentant une entaille centrée. De gauche à droite, le nombre de plis à 0° par rapport aux
plis à 90° décroit [Askari06].
La méthode péridynamique est prometteuse en raison des lois de rupture que l’on n’a
pas à identifier car elles sont traitées comme les autres déformations, de la possibilité de multi-
fissuration, et de l’indépendance théorique au maillage. Cependant un nombre important de
paramètres est à identifier surtout pour le cas du composite, et a fortiori du tissu pour lequel le
nombre de configurations d’interaction est encore plus important.
161
propagation, de par sa rupture brutale, ainsi que son rôle dans l’endommagement sur la
répartition des microfissures et leur densité.
Ainsi, dans le modèle proposé dans ces travaux, on a fait le choix de mailler le pli de tissu
mèche par mèche. S’appuyant sur des constats expérimentaux, on fait l’hypothèse qu’une mèche
ne peut rompre qu’entièrement. La propagation s’effectue par pas spatial correspondant à la
largeur de mèche et on choisit par conséquent de n’utiliser qu’un seul élément dans la largeur de
la mèche. Il est donc apparu judicieux de mailler chaque mèche par une rangée de quadrangles.
162
Chapitre IV – Modélisation numérique de la propagation de coupure
Figure IV.21 – Maillage d’une mèche par des éléments surfaciques séparés par des
éléments de rupture.
Deux portions de mèches sont représentées en Figure IV.22. Les éléments de rupture ont
été étirés, et les deux mèches séparées l’une de l’autre pour la compréhension du schéma. Dans
la pratique les nœuds 2 et 3 sont initialement confondus, de même que les nœuds 6,7,10 et 11, et
les nœuds 14 et 15.
Figure IV.22 – Illustration d’une portion de maillage dans une direction, distendu pour une
meilleure clarté. En configuration saine, les nœuds 2 et 3 sont géométriquement confondus, de
même que les nœuds 6,7,10 et 11, et les nœuds 14 et 15.
Chaque direction de renfort est maillée séparément et les deux nappes notées « chaîne »
et « trame » sont superposées. Les éléments des deux nappes sont liés 2 à 2 comme schématisé
sur la Figure IV.23. L’interface est réalisée par 4 éléments à 2 nœuds, représentés par des
ressorts.
163
Figure IV.23 – Liaison chaîne/trame réalisée par 4 éléments d’interface représentés par
des ressorts (maillage distendu pour la compréhension).
La superposition des deux maillages dans les directions chaîne et trame conduit au
maillage représenté Figure IV.24, toujours de façon distendue pour la clarté du schéma. Au total,
les points géométriques des nœuds du maillage coïncident géométriquement avec 8 nœuds
superposés (4 par direction de maillage).
Figure IV.24 – Schéma du maillage d’un pli, distendu pour la compréhension. En pratique,
le maillage est plan, et les éléments chaîne et trame sont initialement confondus. Les éléments
d’interface relient 2 nœuds initialement confondus, et les éléments de rupture sont de surface
nulle.
164
Chapitre IV – Modélisation numérique de la propagation de coupure
On a observé dans l’étude expérimentale qu’en traction dans la direction des renforts, le
matériau s’endommage par fissuration matricielle dans les mèches transverses à la direction de
la charge. Cette fissuration dégrade les transferts d’efforts entre deux mèches transverses de
même direction, et entre les deux directions de renforts. Pratiquement la quantification de ces
deux types de dégradations semble expérimentalement hors de portée, et suggérerait de faire
appel à une modélisation microscopique d’un volume représentatif d’une zone de tissu fissurée.
Le but ici n’est pas de dissocier l’effet d’une fissure sur la dégradation des liaisons chaîne/chaîne,
trame/trame ou chaîne/trame et de tenter d’obtenir un champ de déformation à l’échelle
microscopique, mais de représenter l’effet de cet endommagement sur la perte de raideur du
stratifié amenant une redistribution des efforts, pour finalement quantifier une déformation à
l’échelle de la mèche.
Les éléments d’interface sont dotés d’une raideur évolutive. Initialement celle-ci est
infinie, et c’est sa diminution qui découple les directions de renforts et permet ainsi de simuler
l’endommagement.
165
En revanche, le drapage [0;45]n nécessitera des éléments d’interface spécifiques (détaillé
en 4.6) pour assurer la liaison entre un maillage orienté à 0/90 et un maillage à ±45. On pourrait
imaginer que cette méthode se transpose aux stratifiés constitués de plis de mêmes directions,
comme discuté en 4.6.3.
- e = 0,31 mm
- w = 0,5 mm
- Ell = Ett = E0 = 21 GPa (tissu équilibré)
- νlt = νtl = ν0 = 0,15
- Glt = G0 = 3,5 GPa
En configuration saine, les éléments d’interface et de rupture sont idéalement de raideur
infinie. Pour des considérations numériques, on adopte la valeur k0 = 106 N/m suffisamment
grande par rapport à la raideur des éléments surfaciques (calcul détaillé en section 4.4).
Les éléments surfaciques sont élastiques orthotropes linéaires. Ils présentent une
épaisseur que l’on fixe à e/2, pour que l’épaisseur d’un élément chaîne et d’un élément trame
corresponde à celle d’un pli. Le comportement orthotrope 2D des éléments surfaciques nécessite
l’identification de 4 paramètres. On choisit d’utiliser les paramètres E1, module dans la direction
de la mèche, E2, module dans la direction transverse (voir Figure IV.25a), G12, module de
cisaillement et ν12, coefficient de Poisson.
On doit avoir l’équivalence entre la matrice de rigidité du pli et celle qui résulte de la
superposition d’un élément orienté à 0° et d’un élément à 90°, ce qui donne la relation :
166
Chapitre IV – Modélisation numérique de la propagation de coupure
(IV.19)
avec
(IV.20)
(IV.21)
(IV.22)
Pour traduire cette hypothèse, [Bizeul09] a calculé le module d’Young E1, en prenant en
compte la fraction volumique de fibres longitudinales sur une section de pli Vfibres_longi, et le
module d’Young des fibres Efibre :
(IV.23)
Le module d’Young des fibres de verre Efibre est une donnée bibliographique, tandis que la
fraction volumique des fibres longitudinales Vfibres_longi a été déterminée à partir d’observations
micrographiques et de seuillages de sections de mèches :
- Efibres = 74 GPa
- Vfibres_longi = 0,19
- e = 0,31 mm
On déduit des équations (IV.22) et (IV.23) les valeurs des coefficients élastiques des
éléments surfaciques :
- E1 = 28 GPa
- E2 = 14 GPa
- G12 = 3,5 GPa
- ν12 = 0,23
167
Le module E1 représente le comportement des fibres seules. Le module E2 représente
donc le complément, à savoir les fibres transverses et la résine qu’elle soit dans les mèches
chaîne, trame ou entre les mèches. Lorsque les deux directions sont totalement découplées,
l’endommagement atteint est 1-E2/(E1+E2) = 1/3. Par conséquent les fibres longitudinales seules
représentent 2/3 de la rigidité longitudinale du pli, et l’ensemble fibres transverses et résine y
contribue pour 1/3. Le module de cisaillement G12 est simplement pris égal à Glt. Chaque élément
contribue de la même manière à la rigidité de cisaillement.
168
Chapitre IV – Modélisation numérique de la propagation de coupure
Pour résumer, par rapport au comportement sain décrit en 4.3, les aspects suivants ont
été pris en compte :
(IV.24)
L’évolution de l’endommagement a été mesurée sur les essais de fatigue présentés dans
la section 2.2.4. Le module d’Young considéré est mesuré sur la pente de la courbe
contrainte/déformation sur la plage de déformation [10-3 ; 3.10-3] en accord avec la norme
[ASTM3039]. Il s’agit donc dans ce cas d’un module d’Young tangent. Cette méthode permet de
comparer un module après N cycles au module nominal E0 ; ce dernier étant classiquement
évalué par la pente à l’origine (sur la plage [10-3; 3.10-3] d’après la norme [ASTM3039]).
169
L’introduction d’une variable d’endommagement doit permettre d’associer une
contrainte à une déformation en cohérence avec la réponse du matériau. Mais le comportement
légèrement non-linéaire que l’on peut constater Figure IV.26 montre que le module tangent est
inapproprié pour remplir ce rôle. Sur cette figure, on voit en trait plein le signal de
contrainte/déformation d’un essai de fatigue, lors de la première charge et pour une charge
après N cycles de fatigue. Les courbes en tirets montrent les pentes des modules tangents, Etan
des deux charges mesurées sur [10-3; 3.10-3]. En utilisant ces pentes, on surestime la contrainte
endurée par l’éprouvette à la déformation maximale de fatigue. Par contre, le module sécant
permet d’associer une contrainte à une déformation donnée en accord avec le comportement du
matériau. Comme les déformations anélastiques ne sont pas prises en compte en traction, le
module sécant Esécant est calculé tel que défini Figure IV.26, c'est-à-dire par la pente entre
l’origine du repère et le point de la courbe de charge à déformation maximale.
Dans le cas présenté Figure IV.26 la baisse de rigidité calculée via le module sécant sera
environ 50% plus importante que celle obtenu en utilisant le module tangent.
Figure IV.26 – Illustration de la différence entre module tangent Etan et module sécant Esécant
utilisé dans la définition de l’endommagement.
Les essais de fatigue présentés en 2.2.4 ont été de nouveau dépouillés en mesurant
l’endommagement pour des nombres de cycles compris entre 102 et 106 cycles à partir du
module sécant. Les résultats sont reportés Figure IV.27, affectés de courbes de tendance. Les
résultats des essais de fatigue réalisés pour l’étude de couplage (ε1 – ε2) ont été ajoutés afin
d’augmenter le nombre de points, notamment pour les hautes déformations (15.10-3). Les points
« fatigue 1 cy » sont calculés à partir du module sécant de la première charge comparé au
module initial E0 (qui lui est tangent) pour chaque essai de fatigue. En comparaison du graphe
Figure II.24 dépouillé en module tangent, les tendances des résultats restent identiques. La prise
en compte du module tangent modifie uniquement les valeurs des endommagements obtenus
170
Chapitre IV – Modélisation numérique de la propagation de coupure
qui sont uniformément supérieures de 40 à 50%. On note également que la dispersion relative
est importante quel que soit le niveau de chargement et le nombre de cycles.
Une asymptote d’endommagement maximum est constatée entre 33 et 35%. Cette valeur
est à rapprocher des parts de rigidité portées par les différents constituants. Le calcul présenté
en 4.3 montre que sous des hypothèses simplificatrices les fibres longitudinales reprennent
environ 2/3 de la rigidité, donc si la fissuration matricielle est telle que les mèches transverses
et la résine ne reprennent plus aucun effort, un endommagement de 33% devrait être atteint. Ce
constat est à modérer car le raisonnement pour obtenir la valeur de E1 correspond à une
représentation idéalisée des mèches et du pli qui néglige certains aspects (présence de points de
tissage, dépliements des fibres sans résine notamment).
171
Figure IV.28 – Elément d’interface utilisé pour lier les deux directions de maillage entre les
nœuds N1 et N2. (a) : cas où le maillage est orienté suivant le repère structure (Oxy). (b) : cas
général.
(IV.25)
(IV.26)
(IV.27)
Dans la pratique parmi les drapages désorientés par rapport au repère structural, seul le
cas θ=45° est utilisé pour notre étude:
(IV.28)
172
Chapitre IV – Modélisation numérique de la propagation de coupure
Soit une mèche de chaîne sous un effort F de traction (voir Figure IV.29), et (c) et (t)
deux éléments superposés, représentant respectivement la direction chaîne et la direction
trame. On raisonne sur un modèle simplifié en une dimension des deux éléments surfaciques liés
par 2 éléments d’interface.
Figure IV.29 – Représentation d’une mèche en traction pour établir la relation raideur –
endommagement, et réduction en un problème 1D sur un ensemble de deux éléments surfaciques
superposés liés par deux éléments ressorts.
(IV.29)
L’effort F se décompose en une force dans l’élément (c), Fc, et une force dans l’élément
(t), Ft qui s’expriment en fonction des déformations εc et εt:
(IV.30)
(IV.31)
173
(IV.32)
(IV.33)
La loi de comportement du ressort pour une raideur k non nulle donne :
(IV.34)
(IV.35)
(IV.36)
(IV.37)
(IV.38)
La relation précédente, permet via les équations (IV.29) et (IV.22) d’établir la relation
entre d et k:
(IV.39)
(IV.40)
La relation (IV.39) est tracée Figure IV.30 en échelle logarithmique. On note deux
asymptotes sur cette courbe. Lorsque k tend vers l’infini, plus précisément k/e >> E2,
l’endommagement est nul ce qui correspond au matériau sain. Au contraire, une raideur très
faible (k/e << E2) donne l’endommagement maximal E2/2E0 = 0,33 (lorsque seules les fibres
longitudinales reprennent des efforts). Entre les deux, la plage de variation de k ayant une
influence significative sur l’endommagement est comprise entre 102 et 105 N/mm.
174
Chapitre IV – Modélisation numérique de la propagation de coupure
(IV.41)
(IV.42)
En d’autres termes, cette loi relie le logarithme de kf à la déformation de fatigue ε par une
loi linéaire, dont la pente p(N) varie linéairement en fonction du nombre de cycles (en
logarithme). Par ailleurs la raideur ne peut dépasser sa valeur initiale k0. La loi (IV.41)-(IV.42)
nécessite l’identification de 4 constantes que sont ε0f, k0f, ap et bp, ce qui est raisonnable pour une
fonction à deux variables. De manière à retrouver les résultats expérimentaux de la Figure IV.27,
ceux-ci ont été identifiés à :
- ε0f = 0,7.10-3
- k0f = 105 N/mm
- ap = -60
- bp = 70
En utilisant l’équation (IV.39), on peut tracer la loi obtenue en termes de courbes
d’endommagement dans un graphe (d – ε), (Figure IV.31). Elle représente globalement
l’évolution de l’endommagement en fonction de la déformation et du nombre de cycles, mais
175
avec une précision moyenne. Des lois d(ε,N) plus complexes permettant d’obtenir une fonction
plus conforme aux résultats expérimentaux ont été envisagées mais afin de limiter le nombre de
paramètres matériau, la loi (IV.41)-(IV.42) a été adoptée.
(IV.43)
Les deux constantes εs et ηs ont été identifiées différemment en sens chaîne et trame.
Dans un soucis de simplification, cette distinction n’est pas prise en compte ici, et les coefficients
identifiés sont :
176
Chapitre IV – Modélisation numérique de la propagation de coupure
calculée par (IV.43) apparait inférieure à la raideur calculée par la loi de raideur en fatigue à
nombre de cycles faible (102).
(IV.44)
L’endommagement d résultant des relations (IV.44) et (IV.39) est alors tracé sur un
graphe (ε – N) (Figure IV.33). On note sur la partie gauche de cette figure, la zone de déformation
pour laquelle l’endommagement est très faible, quel que soit le nombre de cycles. A mesure que
la déformation augmente, le nombre de cycles nécessaire à un endommagement conséquent
diminue. Les changements de pentes aux alentours de 103 cycles correspondent au changement
de régime (ks ou kf) dans la loi (IV.44).
177
Figure IV.33 – Endommagement d en fonction de la déformation de fatigue ε et du nombre
de cycles N.
178
Chapitre IV – Modélisation numérique de la propagation de coupure
Cependant, l’étude expérimentale a montré que les fissures se créent dans les mèches
transverses et semblent se propager instantanément sur une distance de quelques largeurs de
mèche. Le calcul de l’endommagement doit donc prendre en compte des déformations dans des
éléments éloignés de l’élément d’interface. Au vu de l’étude expérimentale sur le matériau, le
caractère instantané de cette propagation est avéré mais son étendue est difficile à quantifier. On
peut estimer d’après les clichés MEB (Figure II.26) qu’elle doit être inférieure à la distance entre
deux points de tissage, soit 7 largeurs de mèches. Devant la difficulté à quantifier cette étendue,
la valeur de 4 largeurs de mèche a été fixée comme compromis (voir Figure IV.35).
Sur cette figure les éléments surfaciques sont numérotés de 1x à 16x (resp. 1y à 16y)
pour les éléments de la direction x (resp. y). On note εix (resp. εiy), la valeur maximale de
déformation longitudinale sur l’élément ix (resp. iy). Le ressort situé au centre de la figure et
représenté en rouge admet un couple de raideur (kx,ky). Le calcul de ces raideurs doit prendre en
compte les déformations dans les 8 éléments surfaciques repassés en rouge, ce que l’on peut
appréhender par un exemple : si l’élément 3x subit une déformation longitudinale ε3x
importante, la fissuration matricielle apparait et se propage instantanément au niveau des
éléments 7 et 11. Ainsi, l’élément d’interface rouge doit tenir compte de la déformation
longitudinale de l’élément 3x, dans le calcul de sa raideur. Par contre une déformation
longitudinale importante dans l’élément 3y génère une fissure qui se propage dans la mèche
transverse qui dans ce cas, est portée par les éléments 1, 2, 3, et 4. Elle n’affecte pas l’élément
d’interface rouge. Finalement, l’expression des déformations εkx et εky prises en compte dans le
calcul des raideurs kx et ky de l’élément d’interface s’écrit :
ε ε ε ε (IV.45)
ε ε ε ε (IV.46)
179
Figure IV.35 – Zone de prise en compte de la déformation dans le calcul des raideurs kx et ky
de l’élément d’interface repassé en rouge.
D’un point de vue pratique, on utilise les déformations extrapolées aux nœuds pour
obtenir les déformations maximales. Chaque élément d’interface a besoin de connaitre les 8
éléments surfaciques « rouges » qui pilotent les valeurs de ses raideurs kx et ky.
Ce cumul se distingue d’une loi de cumul en fatigue de type Miner qui porte sur la durée
de vie en fatigue, et la rupture. On cherche ici à représenter l’évolution d’une variable
180
Chapitre IV – Modélisation numérique de la propagation de coupure
d’endommagement, grandeur mesurable qui modélise une perte de rigidité, et non d’une
variable empirique permettant d’estimer une durée de vie.
La méthode utilisée pour ces travaux s’inspire d’un modèle cumulatif d’endommagement
en fatigue présenté par [VanPaepegem02b] et [VanPaepegem01b]. Des essais de fatigue lui ont
permis d’identifier une loi de variation d’endommagement D fonction du chargement et de
l’endommagement lui-même :
σ (IV.47)
σ (IV.48)
(IV.49)
Dans le cas général, la déformation ε(N) n’est pas constante mais varie avec le nombre
de cycles N et une intégration numérique sera nécessaire :
(IV.50)
Cette déformation sera supposée constante entre N0 et N0+ΔN ce qui est acceptable tant
que ΔN reste faible :
(IV.51)
(IV.52)
où Ñ représente le nombre de cycles nécessaire pour obtenir un endommagement d(N0)
sous une déformation ε(N0), soit :
181
(IV.53)
Ñ se déduit des relations (IV.40), (IV.41) et (IV.42) :
(IV.54)
(IV.55)
(IV.56)
Cette démarche suppose évidemment que la déformation ε(N) soit constante entre N0 et
N0+ΔN. Dans la pratique cette déformation évolue à cause de l’avancée de la coupure et de
l’évolution de l’endommagement lui-même.
L’avancée de la coupure se fait de manière discrète par rupture brutale d’une mèche à
l’issue d’un certain nombre de cycles et on adopte comme pas de temps ΔN le nombre de cycles
nécessaire à la rupture de la mèche suivante. Pendant ce laps de temps ΔN, l’endommagement
lui-même devrait évoluer, et éventuellement atténuer les concentrations de contraintes, et donc
la déformation ε qui intervient dans le calcul du même endommagement. Un pas de temps plus
faible que le nombre de cycles nécessaire à la rupture de la mèche pourrait être adopté pour
affiner le calcul de l’endommagement. Mais au vu des incertitudes sur la mesure de
l’endommagement, et afin de limiter le temps de calcul, cette précision du calcul de
l’endommagement est jugée suffisante.
182
Chapitre IV – Modélisation numérique de la propagation de coupure
On note sur ce graphe que c’est la phase à déformation maximale (ici 0,008) qui est la
plus endommageante. Lorsque la déformation augmente fortement (entre les phases 1 et 2, et
entre 2 et 3) la prise en compte du cumul d’endommagement n’est pas indispensable, car (10 4,5 -
103) ≃ 104,5 et (105 - 103,5) ≃ 105. Lorsque la déformation diminue fortement, l’endommagement
n’évolue quasiment plus (phase 4).
4.4.2. Pseudo-plasticité
Le comportement en cisaillement est fortement non-linéaire avec notamment une
diminution de la raideur et une pseudo-plasticité importante. Ces deux aspects sont pris en
compte dans le comportement des éléments surfaciques par une loi de pseudo-plasticité
endommageable de type puissance identifiée sur les essais réalisés.
183
[Link]. Loi d’endommagement en cisaillement
L’évolution du module de cisaillement sécant GS a été dépouillée en fonction de la
déformation γ (cf Figure II.44). Cette grandeur rapportée au module de cisaillement initial
permet de définir l’endommagement de cisaillement :
(IV.57)
La valeur choisie pour le module de cisaillement initial, G12 = 3500 MPa, est légèrement
inférieure au module mesuré expérimentalement. Ce choix permet de retarder le seuil
d’endommagement, et donc de représenter un comportement linéaire sur une plus grande plage
de déformation. L’évolution de d12(γ) est présentée Figure IV.37 approximée par une loi
exponentielle :
(IV.58)
- d12max = 0,62
- η12 = 40
- γ0 = 0,008
(IV.59)
184
Chapitre IV – Modélisation numérique de la propagation de coupure
(IV.60)
Bien que ces mécanismes de pseudo-plasticité soient fondamentalement différents de la
plasticité sur matériaux métalliques, une loi classique de plasticité est utilisée (par exemple
[Lemaitre09]). Un écrouissage cinématique parait représenter au mieux le comportement du
stratifié (Figure II.43). On définit alors une variable d’écrouissage cinématique X(γp) pour écrire
la fonction de charge :
(IV.61)
avec :
(IV.62)
et :
- A = 105 MPa
- τy = 26 MPa, limite d’élasticité
- α = 0,28
La Figure IV.38 compare la courbe contrainte/déformation expérimentale à la loi
identifiée. La valeur de G12 initiale est inférieure à la rigidité expérimentale, mais permet de
considérer le comportement du matériau comme élastique linéaire sur une plage de
déformation plus étendue.
185
Concernant les aspects pratiques, la déformation pseudo-plastique γp est déterminée en
utilisant un algorithme de Newton-Raphson :
(IV.63)
avec :
(IV.64)
la prédiction élastique, et
(IV.65)
la contrainte actualisée.
Les éléments de rupture comportent 4 nœuds et relient deux éléments surfaciques d’une
même mèche (Figure IV.39). Ils assurent une raideur k théoriquement infinie entre les nœuds 1
et 2, et entre 3 et 4, selon un comportement binaire.
186
Chapitre IV – Modélisation numérique de la propagation de coupure
On choisit donc de définir une variable de cumul de Miner Dm, appelée fraction de vie :
(IV.66)
(IV.67)
Comme expliqué en [Link], la durée de vie des essais de fatigue réalisés est sous-estimée
d’un facteur 2 à 3, et d’autres essais de fatigue en nombre plus réduit, ont permis de dégager les
coefficients suivants :
- ε6 = 0,0076
- γ = 0,129
Ces essais reportés en jaune sur la Figure IV.40 ont été réalisés pour 3 géométries
d’éprouvettes différentes (éprouvettes de couplage (ε1 – ε2), éprouvettes 4 plis de taille standard,
et éprouvettes 2 plis de taille standard à talons biseautés) présentant des concentrations de
contraintes faibles au niveau des talons. Ne s’agissant pas d’une campagne d’essais homogènes
sur un même type d’éprouvette, il est apparu plus approprié de présenter les coefficients
conservatifs en [Link] bien que la modélisation utilise les paramètres plus représentatifs du
comportement du matériau.
187
Figure IV.40 – Loi de Basquin prise en compte dans la modélisation et identifiée sur des
essais présentant des concentrations de contraintes faibles, comparée à la loi de Basquin
conservative présentée en [Link].
Dans la pratique les fractions de vie sont des variables qui dépendent des déformations,
et elles sont donc évaluées dans les éléments surfaciques. Dès qu’une fraction de vie de Miner
atteint 1 dans un élément surfacique, l’élément de rupture le plus proche est rompu.
D’un point de vue plus pratique, pour déterminer la rupture de la Tième mèche, on cherche
l’élément jT pour lequel la rupture en fatigue nécessite le nombre de cycles le plus faible. Pour
cela les fractions de vie de Miner Dm(j) et les déformations εj des éléments j sont considérées, et
l’élément jT est celui qui vérifie:
(IV.68)
Le nombre de cycles NT est celui nécessaire à la Tième rupture d’élément qui a lieu au
niveau de l’élément jT.
188
Chapitre IV – Modélisation numérique de la propagation de coupure
Figure IV.41 – Elément de rupture, et variables de Miner affectées aux points de Gauss
associés.
(IV.69)
avec :
(IV.70)
(IV.71)
189
(IV.72)
telle que :
(IV.73)
(IV.74)
190
Chapitre IV – Modélisation numérique de la propagation de coupure
191
Figure IV.43 – Intersection entre deux mèches matérialisée par le losange noir.
Une technique courante pour lier deux maillages consiste à projeter chaque nœud d’un
maillage esclave sur un maillage maitre et à appliquer des conditions cinématiques à ces nœuds
esclaves. Le déplacement d’un nœud esclave est alors imposé par ceux des nœuds de l’élément
maitre sur lequel il a été projeté et par ses fonctions de forme. Le choix des rôles des maillages se
fait généralement naturellement en attribuant le rôle de maitre au maillage le plus grossier, qui
va imposer les déplacements aux nœuds esclaves.
L’approche employée ici consiste à réaliser l’interface entre les deux éléments contenant
M1 et M2 en liant ces deux points par un élément d’interface à 2 nœuds. La raideur de l’élément
est prise théoriquement infinie pour réaliser la liaison rigide entre les deux maillages.
Dans la pratique, il est possible de connaitre les déplacements des deux points fictifs M1
et M2 à partir des déplacements des nœuds 1 à 8 des deux éléments via les fonctions de formes
de (a) et (b) ainsi que deux paramètres géométriques (que l’on note y90 et y45 sur la Figure IV.44).
Les efforts résultants au niveaux de ces 2 nœuds peuvent alors être extrapolés aux 8 nœuds des
éléments (a) et (b). On simule alors une rigidité entre M1 et M2 sans créer de nœud
supplémentaire.
192
Chapitre IV – Modélisation numérique de la propagation de coupure
Figure IV.44 – Représentation des 8 nœuds supports de l’élément d’interface 0-45 entre les
éléments (a) et (b), et définition des paramètres géométriques associés.
Pour cela on définit donc un élément d’interface à 8 nœuds notés 1 à 8 dont l’aspect est
remarquablement atypique et la géométrie initiale dépend de deux paramètres y90 et y45, et on
lui attribut une matrice de rigidité K dans le plan équivalente à un ressort de translation de
raideur k et de rotation de raideur c entre les nœuds fictifs M1 et M2. Cette rigidité est fonction
des deux paramètres géométriques y90 et y45.
(IV.75)
(IV.76)
de manière à obtenir deux coefficients variant entre -1 et 1.
193
(IV.77)
Avec les notations de la Figure IV.45, les déplacements (uM1,vM1) du point M1 pour un
élément d’ordre 1 s’écrivent :
(IV.78)
(IV.79)
La rotation pourrait être évaluée de façon similaire. Néanmoins, l’élément fait partie
d’une mèche avec un axe privilégié, qui peut être défini par la fibre neutre de la mèche. On écrit
alors la rotation θM1 de M1 comme la rotation de la fibre neutre. Celle-ci s’écrit au premier ordre :
(IV.80)
194
Chapitre IV – Modélisation numérique de la propagation de coupure
(IV.81)
On a :
(IV.82)
Le changement de base n’a pas d’effet sur la matrice A, car celle-ci est diagonale par bloc.
Concernant la rotation de l’élément (b) :
(IV.83)
Finalement :
(IV.84)
(IV.85)
De cette façon,
(IV.86)
195
[Link]. Matrice d’extrapolation des efforts J
La matrice J répartit les efforts et les moments calculés en M1 et M2 sur les 8 nœuds de
l’élément d’interface :
(IV.87)
Les efforts, Fx et Fy, et le moment M, doivent être distribués sur les 4 nœuds. On
extrapole les efforts sur les 4 nœuds en accord avec les fonctions de forme en utilisant la
transposée de la matrice d’interpolation :
(IV.88)
(IV.89)
On répartit les moments suivant l’axe transverse de la mèche, en accord avec la rotation
prise en compte. Le problème consiste à distribuer le moment M en 4 efforts selon x F1x à F4x
équivalents tels que représentés sur la Figure IV.46 :
196
Chapitre IV – Modélisation numérique de la propagation de coupure
(IV.90)
(IV.91)
(IV.92)
(IV.93)
avec :
(IV.94)
(IV.95)
(IV.96)
(IV.97)
avec :
(IV.98)
197
(IV.99)
(IV.100)
La dégradation de ces valeurs matérialise un délaminage entre un pli 0/90 et un pli ±45.
Le délaminage n’est pas observé en essai, que le drapage soit orienté ou quasi-isotrope. Les
raideurs des éléments d’interface 0-45 ne sont donc pas dégradées. Le seul cas où elles le sont a
lieu au moment de la rupture d’une mèche. On choisit d’annihiler localement la liaison entre les
deux plis au niveau de la rupture pour éviter des pontages aberrants de la coupure.
Dans la pratique, la rupture d’une mèche a lieu entre deux éléments surfaciques, et on
imposera des raideurs k et c nulles pour les éléments d’interface 0-45 situés sur l’un des deux
éléments surfaciques en question, c'est-à-dire pour ceux qui présentent un point M qui
appartient géométriquement à l’un des deux éléments surfaciques.
4.6.3. Discussion
L’approche adoptée pour réaliser l’interface entre les maillages 0° et 45° est adaptée à la
stratégie de modélisation choisie. Elle est tributaire de la taille de la mèche, et a conduit au
développement d’un élément d’interface original afin de ne pas discriminer une direction. Cet
élément à 8 nœuds, utilise uniquement des nœuds déjà existants des éléments surfaciques (a) et
(b), il ne rajoute donc pas de degré de liberté supplémentaire au maillage.
Lier deux mèches en un seul point n’est bien sûr pas équivalent à un collage cinématique.
On peut notamment imaginer des cas de deux éléments seuls, pour lesquels un point de transfert
d’effort est insuffisant pour que les déformations de l’un soient transmises à l’autre. Mais ce cas
est uniquement théorique. Dans la pratique, chaque maillage comporte plusieurs mèches, et les
différentes liaisons rigides des interfaces 0-45 imposent aux champs de déformations des deux
drapages à s’accorder l’un à l’autre.
Cette méthode pourrait également se généraliser à l’interface entre deux plis d’une
même direction. On pourrait l’adapter pour réaliser l’interface entre deux plis de direction 0/90
décalés c'est-à-dire pour lesquels les mèches ne se superposent pas.
198
Chapitre IV – Modélisation numérique de la propagation de coupure
d’éléments d’interface adjacents (4.5.3) et d’éléments d’interface 0-45 ([Link]). Les éléments
doivent pouvoir se situer dans l’espace et connaitre les éléments proches. Cela est rendu
possible par une « mise en données », et des routines d’éléments utilisateurs adaptées.
199
Tableau IV.1 – Nomenclature associée à la modélisation.
Variable Définition
200
Chapitre IV – Modélisation numérique de la propagation de coupure
Figure IV.48 - Algorithme simplifié d'une simulation d’un essai de propagation. Les
notations utilisées sont en Tableau IV.1, lorsque le numéro d’itération i est omis, la variable
correspondante est prise au pas de temps convergé.
201
Cet algorithme admet en entrée :
2) Le pas de temps Δt peut éventuellement être variable. Il prend des valeurs comprises
entre 0,01 et 0,1.
3) Le calcul de k(j,t,i) puis de F(j,t,i) se fait pour chaque élément d’interface en fonction
du champ de déformation et du nombre de cycles calculé au pas de temps convergé
précédent. La fonction f donnant k est détaillée en 4.4.1.
5) jT est l’élément qui vérifie l’équation (IV.68). Le temps d’initiation est ici ajouté à ΔNT
pour les pas de temps précédant la rupture de la première mèche (t<1).
7) Pour la rupture, on prend en compte l’état convergé de T-0,1, c'est-à-dire pour lequel
le chargement est maximal et l’équilibre déformation/endommagement est atteint
(voir Figure IV.47). L’élément de rupture est rompu en lui affectant une raideur nulle
(cf 4.5.2).
Pour représenter un problème de structure présentant une crique, seule une petite zone
(jusqu’à quelques dizaines de cm²) peut raisonnablement être maillée à l’aide de ce maillage.
202
Chapitre IV – Modélisation numérique de la propagation de coupure
L’interface entre les deux maillages est assurée par collage cinématique des nœuds
esclaves du maillage fin « MEM » sur les arêtes du maillage grossier « contour ».
La Figure IV.49 présente un exemple de maillage qui illustre cette démarche. Un tronçon
de pale est supposé présenter une crique (une coupure) à l’intrados et à l’extrados, à proximité
du bord de fuite. Le maillage de son revêtement est représenté en éléments membranes
classiques, sur l’ensemble du tronçon sauf au niveau de la coupure, que l’on maille sur les deux
faces par deux maillages « MEM ». Le drapage du revêtement est [±45]n, le maillage à l’échelle de
la mèche est donc orienté en conséquence. Les pré-entailles ont été repassées en rouges.
Ce type de maillage a été réalisé afin de vérifier la faisabilité d’intégration d’un modèle
MEM sur un maillage de pale. La propagation a ainsi été simulée numériquement sur une
structure mais les résultats obtenus n’ont pas pu être comparés à des données expérimentales.
203
Chapitre V – Comparaison essais - calculs
205
5.1. Introduction – Méthodologie
Le modèle présenté au chapitre IV a été éprouvé dans différentes configurations qui
correspondent aux essais de propagation décrits dans le chapitre III. La présentation de ces
simulations fait l’objet du présent chapitre. Il est divisé selon les configurations, classées d’abord
par type d’essai, puis par drapage. Ainsi, la propagation en traction est évaluée en 5.2, tandis que
les simulations en cisaillement sont présentées en 5.3. Enfin une discussion est proposée en fin
de chapitre (5.4). Avant ces développements, quelques précisions communes à tous les cas
concernant les maillages, les propriétés matériau implémentées et les procédures de
dépouillement sont présentées.
- Evl = 21 GPa
- Evt = 21 GPa
- Gvlt =3,5 GPa
- νvlt = 0,15
Les coefficients élastiques suivants sont utilisés pour le comportement de
l’unidirectionnel carbone/époxy (pour représenter les longerons et arêtiers) :
Les coordonnées des nœuds, et la table de connectivités sont générés par une routine
fortran. Le reste du fichier d’instruction contient les données matériau et de chargement. Le
206
Chapitre V – Comparaison essais-calculs
calcul est mené pour un certain nombre de pas de temps, correspondant au nombre de ruptures
de mèches souhaité Tmax. Les calculs ont été réalisés avec le logiciel SAMCEF©, sur des modèles
comportant au total entre 30000 et 200000 DDL et des temps de calcul de l’ordre de quelques
heures CPU.
5.1.2. Dépouillement
A l’issue du calcul, on évalue les résultats en observant la direction empruntée par la
coupure, et en calculant sa vitesse de propagation. Les aspects qualitatifs et quantitatifs de
l’évolution de la vitesse de propagation sont évalués. L’idéal étant de pouvoir obtenir des valeurs
de vitesse de propagation de l’ordre de celles mesurées expérimentalement mais également de
pouvoir reproduire les accélérations et décélérations. Les déformées obtenues ainsi que l’ordre
de rupture des mèches sont également observés.
207
Figure V.1 – Géométrie des maillages utilisés pour la simulation des essais de propagation
en traction sur des drapages [0/90]n. (a), éprouvette de largeur 30 mm, (b), éprouvette de largeur
50 mm.
208
Chapitre V – Comparaison essais-calculs
Figure V.2 – Faciès de rupture en fin de simulation de la propagation en traction sur une
éprouvette [0/90]4 de largeur 50 mm. En bas, représentation des éléments de rupture seuls.
On constate que la coupure s’est globalement propagée dans l’axe de la pré-entaille, bien
qu’elle ne présente pas un faciès rectiligne. Par endroit elle se décale de son axe initial, mais tend
systématiquement à revenir à l’horizontal ce qui traduit une certaine stabilité dans la
propagation.
Sur la Figure V.3 sont présentées les deux directions de maillage (0° et 90°) autour des
lèvres de la fissure après une propagation d’environ 15 mm ainsi que le champ de déformation
longitudinale dans les mèches. Naturellement, les mèches à 90° sont les plus chargées et
présentent les déformations les plus élevées. La déformation maximale que présente la mèche
en pointe de coupure est de l’ordre de 18.10-3. Cette pointe de coupure ne constitue pas une
singularité car le maillage utilisé est semi-discret et l’endommagement des éléments d’interface
par dégradation de raideur donne une pointe de coupure qui n’est pas singulière, et donc non
susceptible de générer des champs asymptotiques infinis.
209
longitudinale sauf au niveau de la pointe de coupure, qui a tendance à leur imposer une
déformée de flexion.
Figure V.3 – Déformation longitudinale dans les mèches orientées à 0° et 90° pour la
simulation d’un essai de propagation en traction sur éprouvette [0/90] 4 de largeur 50 mm, après
une propagation d’environ 15 mm (déformée x3).
La baisse de rigidité des éléments ressorts a été tracée pour la même longueur de
coupure (Figure V.4). On rappelle que chaque élément d’interface présente une rigidité
différente dans les deux directions, k1 et k2. A partir de l’équation (IV.39) on peut estimer les
champs des variables d1 et d2 sur l’ensemble du maillage MEM (voir Figure V.4). Ceux-ci
fournissent un ensemble de valeurs aux nœuds qui ont été interpolées pour la visualisation.
210
Chapitre V – Comparaison essais-calculs
211
Figure V.4 – Champ des variables d’endommagement dans les directions du repère
d’orthotropie calculées à partir des raideurs des éléments d’interface pour la simulation d’un essai
de propagation en traction sur éprouvette [0/90] 4 de largeur 50 mm, après une propagation
d’environ 15 mm. Les axes représentent les coordonnées spatiales en mm.
212
Chapitre V – Comparaison essais-calculs
Figure V.5 – Images obtenues par soustraction d’images de la propagation en traction sur
une éprouvette [0/90]4 après environ 15 mm de propagation.
Enfin, le nombre de cycles correspondant à chaque rupture de mèche a été relevé afin de
déduire la vitesse de propagation qui a été comparée aux résultats de deux éprouvettes
(Figure V.6). Qualitativement on retrouve bien l’allure caractéristique correspondant à cette
géométrie et ce drapage, à savoir une propagation qui s’accélère sur les 10 premiers millimètres,
se stabilise et décroit en fin d’essai de par la présence du longeron. Les valeurs de vitesse
obtenues pour cette configuration sont également très satisfaisantes sauf pour le début de
l’essai. L’initiation est plus difficile à reproduire car elle est dépendante des conditions initiales.
En particulier la position de la pré-entaille par rapport aux mèches des différents plis en pointe
de coupure est un paramètre important qui n’est pas être pris en compte par la modélisation
actuelle.
213
Figure V.6 – Evolution de la vitesse de propagation de coupure en traction obtenue pour
deux éprouvettes [0/90]4 de largeur 50 mm et pour la simulation correspondante.
Les résultats en termes de vitesse de propagation sont reportés en Figure V.7. Il apparait
comme dans le cas précédent (éprouvettes de 50 mm de large) une bonne représentation
qualitative de l’évolution de la vitesse de propagation en accord avec l’essai pour le temps
d’initiation nul. Par contre, aucune différence significative n’est constatée entre les deux
simulations, celle avec un temps d’initiation important ne reproduit pas la diminution de vitesse
de propagation d’un facteur 10 obtenue expérimentalement.
214
Chapitre V – Comparaison essais-calculs
La comparaison entre les deux simulations est poursuivie Figure V.8 qui présente le
faciès de rupture obtenu en fin de simulation, par la visualisation des éléments rompus, ainsi que
sur la Figure V.9 qui compare les déformations longitudinales dans les mèches.
Figure V.8 – Visualisation des éléments de rupture à l’issue de la simulation d’un essai de
propagation en traction sur éprouvette [0/90]2 de largeur 30 mm pour deux nombres de cycles
d’initiation.
On constate que la coupure s’est globalement propagée dans l’axe de la pré-entaille, bien
qu’elle ne présente pas un faciès rectiligne. La pointe de coupure présente des champs de
déformation quasiment identiques. Les zones de sur-contrainte (ou sur-déformation) s’étendent
215
sur la même distance et culminent à des valeurs proches, pour la mèche en pointe de coupure, ce
qui explique la similitude entre les vitesses de propagations.
Figure V.9 – Déformation longitudinale dans les mèches orientées à 90° pour la simulation
d’un essai de propagation en traction sur éprouvette [0/90] 2 de largeur 30 mm après une
propagation d’environ 10 mm, pour deux nombres de cycles d’initiation.
Pourtant, la Figure V.10 montre que le temps d’initiation a bien été pris en compte dans
le calcul de l’endommagement. Le champ de la variable d1 est identique dans les deux cas, mais
pour un nombre de cycles d’initiation Ninit = 4.106, l’endommagement d2 est plus important que
pour Ninit = 0. Plus particulièrement, dans les quelques millimètres en amont de la pointe de
coupure, d2 prend des valeurs similaires, mais au-delà, le temps d’initiation se traduit par des
valeurs supérieures d’environ 0,05.
216
Chapitre V – Comparaison essais-calculs
conséquent très similaire au cas Ninit = 0, et la vitesse de propagation (qui découle directement
du champ de déformation) n’apparait donc pas impactée par le temps d’initiation.
Figure V.10 – Champ des variables d’endommagement dans les directions du repère
d’orthotropie calculées à partir des raideurs des éléments d’interface sur des éprouvettes [0/90] 2
pour 2 nombres de cycles d’initiation. Les axes représentent les coordonnées spatiales en mm, et
les flèches indiquent la position de la pointe de coupure.
217
0,1) n’a aucune influence sur D(N+ΔN). Cela signifie que quel que soit le temps d’initiation, tant
que l’endommagement engendré est faible, il n’aura que peu d’influence sur l’état
d’endommagement au voisinage de la pointe de coupure, donc sur le champ de déformation en
pointe de coupure, et donc sur la vitesse de propagation.
Sur le champ d2 de la Figure V.10 sont tracées les limites de la zone d’endommagement
supérieur à 0,2. L’effet du temps d’initiation est visible pour les 5 premiers millimètres de
propagation environ. Au-delà l’endommagement est aussi étendu avec ou sans temps
d’initiation.
[Link]. Discussion
La modélisation présentée donne de bons résultats sur deux géométries différentes
d’éprouvette [0/90]n. La direction de propagation horizontale est retrouvée et les zones
d’endommagement obtenues sont comparables à celles observées expérimentalement. On
retrouve à la fois l’allure des vitesses de propagation obtenues expérimentalement, mais
également l’ordre de grandeur.
218
Chapitre V – Comparaison essais-calculs
Figure V.11 – Maillages utilisés pour la simulation des essais de propagation en traction sur
des drapages [±45]2. (a), éprouvette de largeur 30 mm, (b), éprouvette de largeur 50 mm.
La Figure V.12 présente le faciès de rupture à la fin des deux simulations, par l’affichage
des éléments de rupture. Ceux qui ont rompu présentent une surface non nulle et marquent le
chemin de la coupure qui se propage dans l’axe de la pré-entaille pour ce drapage. Sur [±45]2 on
ne bénéficie plus de l’effet « fermeture éclair » décrit en [Link] des mèches horizontales qui
permettent de guider la propagation, cette direction est donc plus délicate à obtenir
numériquement.
219
Figure V.12 - Visualisation des éléments de rupture à l’issue de la simulation d’un essai de
propagation en traction sur éprouvette [±45] 2 de largeur 30 mm (a) et 50 mm (b), et comparaison
avec le faciès de rupture d’une éprouvette de largeur 50 mm (c).
220
Chapitre V – Comparaison essais-calculs
Figure V.13 – Champ des variables d’endommagement dans les directions du repère
d’orthotropie calculées à partir des raideurs des éléments d’interface après une propagation
d’environ 10 mm sur une éprouvette [±45]2 de largeur 30 mm. Les axes représentent les
coordonnées spatiales en mm, et les flèches indiquent la pointe de coupure.
221
Figure V.14 – Clichés obtenus par soustraction d’images de la propagation en traction sur
une éprouvette [±45]2 de largeur 50 mm après environ 10 mm de propagation.
Le champ de déformation longitudinale dans les deux directions de renforts est tracé
Figure V.15. La concentration de contraintes (ou de déformations) dans les deux maillages est
orientée à 45°. Elle est localisée sur une zone plus restreinte que pour le cas [0/90]n, ce qui est
confirmé sur les agrandissements de la Figure V.16.
222
Chapitre V – Comparaison essais-calculs
Figure V.15 – Déformation longitudinale dans les mèches orientées à +45° et -45° pour la
simulation d’un essai de propagation en traction sur éprouvette [±45]2 de largeur 50 mm, après
une propagation d’environ 15 mm (déformée x1).
Les images Figure V.16a-b illustrent le glissement des mèches, les unes par rapport aux
autres, rendu possible par la modélisation semi-discrète de l’endommagement par éléments
d’interface dégradables. Les mèches de pointe de coupure présentent un champ de flexion
important dans la largeur. Une flexion importante s’ajoute à la traction. La Figure V.16c met
l’emphase sur l’une des mèches de pointe de coupure, qui suit une double courbure.
Manifestement, si une telle déformée était imposée dans la réalité, la mèche en pointe de
coupure s’endommagerait (décohésion fibre/matrice, fissuration matricielle) de sorte à atténuer
cette flexion.
Le comportement d’une mèche en flexion semble donc mal représenté par une rangée de
quadrangles de degré 1, en particulier parce que les mèches sont elliptiques, ondulent et
pivotent autour de leur axe (Figure II.19). Il faudrait alors modifier les coefficients élastiques des
éléments surfaciques pour prendre en compte la géométrie des mèches, ou recourir à des
éléments de type poutre de manière à mieux représenter la flexion de la mèche.
223
Comme il est délicat d’estimer l’endommagement d’une mèche en flexion, nous avons
considéré qu’une flexion trop importante endommageait la mèche et que la déformation
maximale n’était pas représentative de la réalité. Dans la pratique, la déformation de traction
semble plus représentative de la réalité et a donc été prise en compte pour le calcul de la vitesse
de propagation (voir 4.5.2).
Figure V.16 – Champ de déformation longitudinale pour des éprouvettes [±45] 2. (a) et (b) :
pour largeur 30 mm après environ 10 mm de propagation. (c), double courbure des mèches en
front de fissure (éprouvette largeur 50 mm) après 30 mm de propagation environ (déformée x3).
Ces résultats de vitesse de propagation sont comparés aux résultats expérimentaux sur
les Figure V.17 et Figure V.18.
La vitesse est bien reproduite dans le cas 50 mm. La décroissance est certes moins
marquée dans la simulation mais la vitesse moyenne de propagation est bien retrouvée. L’allure
présente quelques oscillations. Par contre, pour l’éprouvette de largeur 30 mm, la vitesse de
propagation est surestimée par la modélisation d’un facteur 5 environ.
224
Chapitre V – Comparaison essais-calculs
L’étude plus précise des éprouvettes fournit un élément d’explication de cette différence
quantitative importante entre les configurations. En effet, les plis utilisés pour les éprouvettes de
30 mm et 50 mm, n’ont pas été découpés selon la même orientation. Par polissage post-mortem,
on aperçoit que les 3 éprouvettes de largeur 50 mm ont été drapées avec des plis à +45°
(Figure V.19b) et présentent donc des points de tissages orientés selon l’axe de la pré-entaille.
225
Par contre les 2 éprouvettes de largeur 30 mm ont été drapées avec des plis à -45° et présentent
donc un alignement des points de tissage selon un axe incliné de 26° par rapport à l’horizontale
(Figure V.19a). Dans la section 3.4.4, une différence significative dans la vitesse de propagation
pour des drapages quasi-isotropes a été observée suivant l’orientation des plis à ±45°. En
particulier, les drapages contenant des plis à -45° présentant un axe des points de tissages à 26°
par rapport à la pré-entaille ont donné des vitesses de propagation environ 5 fois plus lentes que
lorsque ces plis sont orientés à +45°. Ces conclusions ne sont pas définitives de par le nombre
d’essais limité, et ne sont pas forcément transposable en l’état sur drapage [±45]2, mais
pourraient expliquer la différence essai/calcul sur largeur 30 mm (Figure V.18).
Figure V.19 – Faciès de rupture des éprouvettes de drapages [±45] 2. (a) : largeur 30 mm,
(b), largeur 50 mm.
L’angle de 26° dans le cas 30 mm ne peut pas être obtenu par modélisation puisqu’il
découle du motif de tissage qui n’est pas reproduit. On peut donc considérer que la direction de
propagation horizontale obtenue numériquement est cohérente avec les essais même pour la
largeur 30 mm.
226
Chapitre V – Comparaison essais-calculs
Figure V.20 – Faciès de rupture en fin de simulation de la propagation en traction sur une
éprouvette [0;45]s mis en évidence par le champ de déplacement selon l’axe y pour les mèches à
90° et à 45°.
Les éléments d’interface à 8 nœuds liant les deux maillages sont représentés en bleu
Figure V.21 superposés au maillage 0/90 (a), et au maillage ±45 (b). Cette figure permet de
constater la géométrie singulière des éléments d’interface et leur distribution spatiale
irrégulière de par la manière dont ils sont adaptés aux deux maillages. Ces orientations peuvent
avoir créé les conditions d’un décalage de la direction de propagation par rapport à l’horizontale.
227
Figure V.21 – (a), maillage non déformé à 0/90 et éléments d’interface 0-45, (b), maillage
non déformé à ±45 avec les mêmes éléments d’interface 0-45.
Les déformations des mèches en pointe d’entaille sont plus importantes pour ce drapage
avec des valeurs dépassant 0,02. Enfin la zone de surcontrainte est plus tourmentée que pour les
simulations à une seule direction de maillage.
228
Chapitre V – Comparaison essais-calculs
Figure V.22 – Déformation longitudinale dans les mèches orientées à 0°, 90°, 45° et -45°
pour la simulation d’un essai de propagation en traction sur éprouvette [0;45]s de largeur 50 mm,
et après une propagation d’environ 15 mm (déformée x1).
Concernant les interfaces 0/90 et +45/-45, ce maillage met en jeu des raideurs dans 4
directions. Les cartes d’endommagement sont représentées Figure V.23. Les variables d1 et d2
ont été calculées à partir des raideurs des éléments d’interface entre les directions 0 et 90,
tandis que d1_45 et d2_45 correspondent à l’endommagement entre les directions +45 et -45.
Sur les différentes cartes on retrouve l’aspect général observé pour les essais de
propagation sur les autres drapages avec une zone d’endommagement qui s’étend sur quelques
millimètres autour de la coupure. La direction 2, celle de la traction, présente un gradient
d’endommagement plus important ; l’endommagement selon « 45_2 » semble plus important
que selon « 45_1 » malgré une direction de propagation qui s’oriente plutôt perpendiculairement
à l’axe 45_1.
Autour de cette zone proche de la coupure, la zone faiblement endommagée présente des
irrégularités, sur des zones allongées dont la largeur est celle des mèches et ce pour toutes les
directions. La répartition inhomogène des éléments d’interface 0-45 est sans doute à l’origine de
ce bruitage de l’endommagement. Tant que celui-ci reste faible et éloigné de la zone de
propagation, on peut penser que son influence sur le résultat reste limitée.
229
Figure V.23 – Champ des variables d’endommagement dans les directions des repères
d’orthotropie des plis 0/90 et ±45 calculées à partir des raideurs des éléments d’interface après
une propagation d’environ 15 mm sur une éprouvette [0;45]s. Les axes représentent les
coordonnées spatiales en mm, et les flèches indiquent la position de la pointe de coupure.
230
Chapitre V – Comparaison essais-calculs
L’inclinaison obtenue par modèle n’est pas naturelle car le motif de tissage n’est pas pris
en compte. Elle ne peut pas être rapprochée de celle que l’on obtient expérimentalement. La
prise en compte de ce motif est discutée en fin de chapitre.
231
Figure V.25 – Maillage utilisé pour la simulation des essais de propagation en cisaillement
sur des drapages [0/90]4. La pointe de pré-entaille sur le maillage MEM a été repassée en rouge.
Cette simulation a été menée jusqu’au pas de temps Tmax = 30. En fin d’essai, le faciès de
rupture obtenu numériquement est présenté Figure V.26b sous forme de champ de déplacement
selon l’axe y. L’évolution de la coupure se fait selon une direction cohérente avec celle obtenue
expérimentalement (a).
On remarque alors que sur les 15 premières ruptures, une seule (n°9) concerne une
mèche horizontale. La coupure se propage ainsi d’abord dans le maillage vertical avant de se
propager dans les mèches horizontales.
232
Chapitre V – Comparaison essais-calculs
Figure V.26 – (a), faciès de rupture en fin d’essai (b), faciès de rupture de la simulation
correspondante, (c) rappel du maillage d’un pli 0/90 mettant en évidence les éléments de rupture,
(d), évolution du faciès de rupture en cours de simulation par numérotation des éléments de
rupture rompus. Les éléments de rupture bleus (rouges) correspondent à des rupture de mèches
verticales (horizontales).
La déformée après 15 ruptures de mèches est montrée Figure V.27. Une seule mèche à 0°
a rompu, les autres présentent une double inflexion, dans une zone vraisemblablement très
endommagée puisque les glissements constatés entre les mèches sont très importants.
Cette déformée n’est pas physiquement acceptable. Lorsque la coupure se propage dans
le tissu, elle coupe des mèches dans les deux directions. On retrouve ici la déformée des mèches
constatée en traction sur drapage quasi-isotrope. Elle semble mettre en évidence une lacune
dans le critère de rupture utilisé puisque certaines mèches sont fortement chargées en flexion
sans rompre. En effet, le critère de rupture porte sur la déformation aux points de Gauss et non
sur la déformation extrapolée aux nœuds. Une partie seulement de la flexion est prise en compte
(voir 4.5.2). La Figure V.26d montre que ces mèches horizontales finissent par rompre dans un
second temps (n° 20, 23-25, 27-28 notamment).
233
Figure V.27 – Visualisation du champ de déformation longitudinale au pas de temps T
correspondant à 15 ruptures de mèches (déformée x3).
234
Chapitre V – Comparaison essais-calculs
Pour la variable d’endommagement d2, une zone sous la coupure se dessine, avec une
intensité moindre. Sa forme est similaire à la zone d’endommagement située sous la coupure de
la vue expérimentale mais semble moins prononcée en simulation (d2 ~ 0,05).
235
Figure V.29 – Comparaison des variables d’endommagement avec l’endommagement
obtenu par soustraction d’images expérimentales. Les axes représentent les coordonnées spatiales
en mm, et les flèches indiquent la position de la pointe de coupure.
236
Chapitre V – Comparaison essais-calculs
l’influence du motif de tissage. Compte tenu de cet élément qui n’est pas représenté, la direction
globale obtenue par modèle est satisfaisante.
Figure V.30 – Simulation des essais de propagation en cisaillement sur des drapages
[0;45]s. (a), géométrie du maillage (b), déformée, (c), image expérimentale.
Globalement des déformations élevées (>0,02) sont obtenues tant pour ce pas de temps
que pour d’autres. Cela se traduit par une vitesse de propagation rapide tracée en Figure V.32,
qui apparait éloignée des résultats expérimentaux.
237
Figure V.31 - Déformation longitudinale dans les mèches orientées à 0°, 90°, 45° et -45°
pour la simulation d’un essai de propagation en cisaillement sur éprouvette [0;45]s après 45
ruptures de mèche (déformée x2).
238
Chapitre V – Comparaison essais-calculs
Enfin, la Figure V.33 montre les cartes des variables d’endommagement dans les 4
directions. Le même phénomène de bruitage qu’en propagation de traction sur drapage quasi-
isotrope apparait. Celui-ci est plus important qu’en traction (voir Figure V.23), avec des
variations d’endommagement supérieures à 0,1 entre deux mèches adjacentes pour d1 et d45_2.
Ces 2 variables sont pilotées par les déformations longitudinales dans les éléments surfaciques
des mèches à 0° et -45° respectivement, et on retrouve des perturbations similaires sur les
champs de déformation observés notamment dans ces mèches à 0° et -45° (Figure V.31a et d).
Les éléments d’interface 0-45 sont manifestement à l’origine de ces perturbations. Leur
répartition hétérogène affecte les champs de déformation sur une amplitude acceptable pour la
traction, mais trop élevée pour le cisaillement. On rappelle que leur comportement est binaire
(état sain ou état rompu) et leur non-endommagement pourrait donc être à l’origine de ces sur-
contraintes.
239
Figure V.33 – Champ des variables d’endommagement dans les directions des repères
d’orthotropie des plis 0/90 et ±45 calculées à partir des raideurs des éléments d’interface après
une propagation en cisaillement d’environ 15 mm sur une éprouvette [0;45]s. Les axes
représentent les coordonnées spatiales en mm, et les flèches indiquent la position de la pointe de
coupure.
240
Chapitre V – Comparaison essais-calculs
isotropes ont été étudiés, la prise en compte des différents endommagements a été
rationnalisée, et le caractère saccadé des vitesses de propagation a été lissé. Le champ
d’application de la modélisation a ainsi été étendu à différentes configurations de sollicitations
planes, correspondant à différents essais réalisés en laboratoire.
On note tout d’abord que sans prendre en compte le motif de tissage dans la
modélisation, les directions de propagation sont généralement bien obtenues. Seul le cas de la
traction [0;45]s donne une propagation qui dévie de l’axe de l’entaille initiale alors que
numériquement elle ne devrait pas en détourner. Dans les autres configurations, les directions
de propagation correspondent globalement aux directions expérimentales, au motif de tissage
près.
En traction, l’unique configuration qui donne des résultats numériques éloignés des
vitesses mesurées en essai concerne l’éprouvette [±45]2 de largeur de 30 mm pour laquelle la
vitesse simulée est surestimée sur l’ensemble de l’essai, d’un rapport 5 environ. Cependant, la
distinction qui existe entre l’orientation +45° et -45°, dans les directions des points de tissage,
pourrait expliquer cette différence. Cette hypothèse mériterait une comparaison expérimentale
sur même géométrie d’éprouvette entre des drapages [+45/-45]2 et [-45/+45]2 voire sur des
éprouvettes à un pli ([+45/-45] et [-45/+45]).
Sur les deux configurations à [0;45]s, les éléments d’interface entre les deux directions
remplissent bien leur rôle de transmission des efforts d’interface, n’alourdissent pas la
modélisation par des degrés de liberté supplémentaires, mais leur répartition irrégulière
perturbe les champs de déformation et donc d’endommagement, de manière non négligeable
dans le cas du cisaillement. Ils sont peut-être à l’origine de la direction de propagation déviée
pour le cas de traction [0;45]s. L’interaction entre les plis à 0° et 45° est difficile à appréhender,
en particulier l’influence de l’un sur l’endommagement de l’autre, dont la manifestation est
perceptible en 3.5.3. Sans ces informations, le bruitage des résultats constaté amène à revoir la
définition de l’élément d’interface 0-45 de manière à ce que sa répartition sur le maillage MEM
241
soit plus uniforme. L’introduction d’un endommagement dans la liaison 0-45 est également
envisageable afin de réduire ce bruit.
Enfin, l’effet du temps d’initiation sur la vitesse de propagation n’a pas été retrouvé par
cette modélisation. La loi d’endommagement définie, et en particulier les hypothèses
considérées pour le calcul de son évolution sous chargement variable, en est peut-être la cause.
Néanmoins, l’ensemble de ces résultats est encourageant pour poursuivre dans cette voie
et envisager des perspectives de développement adaptées au comportement du matériau. De
plus, le point de vue pragmatique actuellement adopté dans ce travail doit être maintenu afin de
cibler le modèle pour une utilisation industrielle.
L’étude expérimentale a montré que le motif de tissage et la taille de la mèche sont deux
éléments influençant les phénomènes de propagation de coupure, étant donné que la
propagation s’effectue par rupture mèche par mèche, et que cette propagation suit un chemin
guidé par le motif de tissage. Le modèle proposé ne tient compte que du premier des deux
aspects, la taille de la mèche. Les différents résultats de direction de propagation, et les
différences observées entre les plis orientés à +45° et -45°, amènent à se pencher sur la question
du motif de tissage et de son éventuelle prise en compte. On peut alors considérer plusieurs
approches.
Par ailleurs, bien qu’il présente des limites, le principe de l’élément 0-45 a montré des
caractéristiques intéressantes de par son caractère adaptable à tous les maillages. En effet, ce
242
Chapitre V – Comparaison essais-calculs
type d’interface pourrait se généraliser à l’interface entre deux directions de renforts d’un même
pli (liaison réalisée par les éléments d’interface à deux nœuds). Cela permettrait de simplifier le
modèle en limitant le nombre de types d’éléments. En outre, cette définition pourrait être
généralisée pour deux plis d’une même direction. Au lieu de considérer un problème à deux plis
de même direction comme deux fois le même problème, les plis pourraient être maillés
indépendamment et liés par ce type d’interface. Cela permettrait par exemple d’évaluer l’effet
d’un décalage des mèches, (et éventuellement des motifs de tissage) entre les plis d’une même
direction.
Enfin, outre l’effet du temps d’initiation, la dispersion obtenue sur certains essais amène
à considérer l’introduction d’incertitudes statistiques dans les paramètres du modèle. Une
variable aléatoire pourrait être ajoutée aux variables d’endommagement et de Miner. Ainsi, le
lancement successif de plusieurs simulations permettrait d’évaluer la robustesse des résultats,
en termes de direction et vitesse de propagation, en fonction des sollicitations, géométrie, et du
nombre de plis.
243
Conclusions et perspectives
245
théoriquement applicable à d’autres composites tissés, en utilisant leur taille de mèche pour le
maillage et leurs paramètres identifiés sur éprouvettes saines.
Les résultats sur l'effet du temps d'initiation ne sont pas en accord avec ceux observés
expérimentalement. Par conséquent, ils remettent en question certaines hypothèses du modèle,
telles que la définition de la loi de cumul d’endommagement en fatigue et son indépendance au
chemin suivi.
L’idée de la prise en compte du motif de tissage a également été évoquée. Elle pourrait
permettre d’obtenir des faciès de rupture et des directions de propagation plus conformes aux
essais, et de rendre compte des éventuels effets du motif de tissage sur la vitesse de propagation.
246
Conclusions et perspectives
dispersion dans la vitesse de propagation sur ce même cas pourrait être étudiée en comparant
des éprouvettes préalablement endommagées en fatigue à des éprouvettes saines avant de
pratiquer les entailles. On évaluerait ainsi l’effet de l’endommagement sur la vitesse de
propagation.
Aussi, l’effet de l’orientation des plis à ±45° sur les drapages quasi-isotropes est à
confirmer sur drapage [±45]n pour différentes géométries. Si cette direction a une influence sur
la vitesse de propagation, le choix de la direction la moins pénalisante pour le drapage des pales
serait un progrès à moindre coût pour leur durée de vie en présence d’une coupure.
Enfin, ce modèle a pu être intégré sur une structure de pale, bien que les résultats
obtenus n’aient pas été confrontés à des données expérimentales. Il serait indispensable pour
l’amélioration et la validation du modèle de poursuivre dans cette voie, et d’évaluer sa
robustesse sur un essai de structure réelle.
247
Références bibliographiques
Alif N, Carlsson LA (1997) Failure mechanisms of woven carbon and glass composites. Astm
Spec Tech Publ 471–493.
Allix O, Ladeveze P (1992) Interlaminar interface modelling for the prediction of delamination.
Compos Struct 22:235–242.
Amijima S, Fujii T, Hamaguchi M (1991) Static and fatigue tests of a woven glass fabric
composite under biaxial tension-torsion loading. Composites 22:281–289.
Askari E, Xu J, Silling S (2006) Peridynamic analysis of damage and failure in composites. 44th
Aiaa Aerosp. Sci. Meet. Exhib. Reno Nev. Rest. Va Aiaa
ASTM (1996) ASTM D3479/D3479M-96: Test Method for Tension-Tension Fatigue of Polymer
Matrix Composite Materials.
ASTM (2000) ASTM D3039/D3039M-00: Standard test method for tensile properties of
polymer matrix composite materials".
ASTM (2001) ASTM D3518/D3518M-01: Standard Test Method for In-Plane Shear Response of
Polymer Matrix Composite Materials by Tensile Test of 45° Laminate.
ASTM (1994) ASTM D4255/D4255M-83: Standard Guide for Testing In-Plane Shear Properties
of Composite Laminates (Two- and Three-Rail Shear Test).
Bathias C (1991) Fatigue des matériaux composites à matrice polymère et fibres longues. Tech
Ingénieur Plast Compos AM6:A7760.1–A7760.11.
Belmonte HM., Manger CI., Ogin S., et al. (2001) Characterisation and modelling of the notched
tensile fracture of woven quasi-isotropic GFRP laminates. Compos Sci Technol 61:585–
597.
Belytschko T, Moës N, Usui S, Parimi C (2001) Arbitrary discontinuities in finite elements. Int J
Numer Methods Eng 50:993–1013.
Beuth Jr. JL, Herakovich CT (1989) Analysis of crack extension in anisotropic materials based on
local normal stress. Theor Appl Fract Mech 11:27–46.
249
Bishop S (1989) Strength and failure of woven carbon-fibre reinforced plastics for high
performance applications. Compos Mater Ser 3:173–207.
Bizeul M, Bouvet C, Barrau JJ, Cuenca R (2010) Influence of woven ply degradation on fatigue
crack growth in thin notched composites under tensile loading. Int J Fatigue 32:60–65.
Bizeul M, Bouvet C, Barrau JJ, Cuenca R (2011) Fatigue crack growth in thin notched woven
glass composites under tensile loading. Part I: Experimental. Compos Sci Technol
71:289–296.
Blès G, Gadaj SP, Nowacki WK, Tourabi A (2002) Experimental study of a PA66 solid polymer in
the case of cyclic shear loading. Warszawa, pp 155–174
Bouvet C, Castanié B, Bizeul M, Barrau J-J (2009) Low velocity impact modelling in laminate
composite panels with discrete interface elements. Int J Solids Struct 46:2809–2821.
Buczek MB, Herakovich CT (1985) A Normal Stress Criterion for Crack Extension Direction in
Orthotropic Composite Materials. J Compos Mater 19:544–553.
Carloni C, Nobile L (2004) Fracture response of cracked orthotropic plates. Convegno Igf Xvii
Chang F-K, Scott RA, Springer GS (1984) The Effect of Laminate Configuration on Characteristic
Lengths and Rail Shear Strength. J Compos Mater 18:290–296.
Chrysochoos A, Maisonneuve O, Martin G, et al. (1989) Plastic and dissipated work and stored
energy. Nucl Eng Des 114:323–333.
Coker EG (1912) An Optical Determination of the Variation of Stress in a Thin Rectangular Plate
Subjected to Shear. Proc R Soc Lond Ser 86:291–319.
Curtis PT (1989) The fatigue behaviour of fibrous composite materials. J Strain Anal Eng Des
24:235–244.
Daggumati S, De Baere I, Van Paepegem W, et al. (2013) Fatigue and post-fatigue stress–strain
analysis of a 5-harness satin weave carbon fibre reinforced composite. Compos Sci
Technol 74:20–27.
250
Daggumati S, De Baere I, Van Paepegem W, et al. (2010) Local damage in a 5-harness satin
weave composite under static tension: Part I – Experimental analysis. Compos Sci
Technol 70:1926–1933.
Daggumati S, Voet E, Van Paepegem W, et al. (2011) Local strain in a 5-harness satin weave
composite under static tension: Part I – Experimental analysis. Compos Sci Technol
71:1171–1179.
De Baere I, Van Paepegem W, Degrieck J (2008a) Comparison of the modified three‐rail shear
test and the [(+45°,−45°)]ns tensile test for pure shear fatigue loading of carbon fabric
thermoplastics. Fatigue Fract Eng Mater Struct 31:414–427.
De Baere I, Van Paepegem W, Degrieck J (2008b) Design of a modified three-rail shear test for
shear fatigue of composites. Polym Test 27:346–359.
De Carvalho NV, Pinho ST, Robinson P (2011) An experimental study of failure initiation and
propagation in 2D woven composites under compression. Compos Sci Technol 71:1316–
1325.
De Greef N, Gorbatikh L, Godara A, et al. (2011a) The effect of carbon nanotubes on the damage
development in carbon fiber/epoxy composites. Carbon 49:4650–4664.
De Greef N, Gorbatikh L, Lomov SV, Verpoest I (2011b) Damage development in woven carbon
fiber/epoxy composites modified with carbon nanotubes under tension in the bias
direction. Compos Part Appl Sci Manuf 42:1635–1644.
251
Dyer KP, Isaac DH (1998) Fatigue behaviour of continuous glass fibre reinforced composites.
Compos Part B Eng 29:725–733.
El-Hajjar R, Haj-Ali R (2004) In-plane shear testing of thick-section pultruded FRP composites
using a modified Arcan fixture. Compos Part B Eng 35:421–428.
Erdogan F, Sih GC (1963) On the crack extension in plates under plane loading and transverse
shear. J Basic Eng 85:19–27.
Fiedler B, Hojo M, Ochiai S, et al. (2001) Failure behavior of an epoxy matrix under different
kinds of static loading. Compos Sci Technol 61:1615–1624.
Fruehmann RK, Dulieu-Barton JM, Quinn S (2010) Assessment of fatigue damage evolution in
woven composite materials using infra-red techniques. Compos Sci Technol 70:937–946.
Gagel A, Lange D, Schulte K (2006) On the relation between crack densities, stiffness
degradation, and surface temperature distribution of tensile fatigue loaded glass-fibre
non-crimp-fabric reinforced epoxy. Compos Part Appl Sci Manuf 37:222–228.
Gao F, Boniface L, Ogin SL, et al. (1999) Damage accumulation in woven-fabric CFRP laminates
under tensile loading: Part 1. Observations of damage accumulation. Compos Sci Technol
59:123–136.
Goel A, Chawla KK, Vaidya UK, et al. (2009) Characterization of fatigue behavior of long fiber
reinforced thermoplastic (LFT) composites. Mater Charact 60:537–544.
Gope PC, Sharma SP, Srivastava AK (2010) Prediction of Crack Initiation Direction for Inclined
Crack Under Biaxial Loading by Finite Element Method. J Solid Mech 2:257–266.
Gutkin R, Laffan ML, Pinho ST, et al. (2011) Modelling the R-curve effect and its specimen-
dependence. Int J Solids Struct 48:1767–1777.
Hallett SR, Green BG, Jiang WG, Wisnom MR (2009) An experimental and numerical
investigation into the damage mechanisms in notched composites. Compos Part Appl Sci
Manuf 40:613–624.
252
Hansen U (1999) Damage Development in Woven Fabric Composites during Tension-Tension
Fatigue. J Compos Mater 33:614–639.
Harris CE, Morris DH (1986) Effect of laminate thickness and specimen configuration on the
fracture of laminated composites. Compos Mater Test Des Astm Stp 893:177–195.
Hettich T, Hund A, Ramm E (2008) Modeling of failure in composites by X-FEM and level sets
within a multiscale framework. Comput Methods Appl Mech Eng 197:414–424.
Higashino M, Takemura K, Fujii TJ (1995) Strength and damage accumulation of carbon fabric
composites with a cross-linked NBR modified epoxy under static and cyclic loadings.
Compos Struct 32:357–366.
Hiver JM, Dahoun A, G’Sell C (2002) Endommagement plastique dans les mélanges de
polymères: mesure in-situ en temps réel dans la striction.
Hochard C, Lahellec N, Bordreuil C (2007) A ply scale non-local fibre rupture criterion for CFRP
woven ply laminated structures. Compos Struct 80:321–326.
Hochard C, Payan J, Bordreuil C (2006) A progressive first ply failure model for woven ply CFRP
laminates under static and fatigue loads. Int J Fatigue 28:1270–1276.
Hochard C, Thollon Y (2010) A generalized damage model for woven ply laminates under static
and fatigue loading conditions. Int J Fatigue 32:158–165.
John S, Herszberg I, Coman F (2001) Longitudinal and transverse damage taxonomy in woven
composite components. Compos Part B Eng 32:659–668.
Jones CJ, Dickson RF, Adam T, et al. (1984) The Environmental Fatigue Behaviour of Reinforced
Plastics. Proc R Soc Lond Math Phys Sci 396:315–338.
Kaddour AS, Al-Salehi FAR, Al-Hassani STS, Hinton MJ (1995) Simultaneous determination of
in-plane shear and transverse moduli of unidirectional composite laminae at different
strain rates and temperatures. Compos Sci Technol 53:431–444.
Kawai M, Matsuda Y (2012) Anisomorphic constant fatigue life diagrams for a woven fabric
carbon/epoxy laminate at different temperatures. Compos Part Appl Sci Manuf 43:647–
657.
Kawai M, Morishita M, Fuzi K, et al. (1996) Effects of matrix ductility and progressive damage
on fatigue strengths of unnotched and notched carbon fibre plain woven roving fabric
laminates. Compos Part Appl Sci Manuf 27:493–502.
253
Kensche CW (2006) Fatigue of composites for wind turbines. Int J Fatigue 28:1363–1374.
Kriz RD, Muster WJ (1985) Mechanical-damage effects in woven laminates at low temperatures.
Mater Stud Magn Fusion Energy Appl Low Temp-Viii 49–86.
Krueger R (2002) The Virtual Crack Closure Technique: History, Approach and Applications.
Ladeveze P, LeDantec E (1992) Damage modelling of the elementary ply for laminated
composites. Compos Sci Technol 43:257–267.
Laffan MJ, Pinho ST, Robinson P (2013) Mixed-mode translaminar fracture of CFRP: Failure
analysis and Fractography. Compos Struct 95:135–141.
Laffan MJ, Pinho ST, Robinson P, McMillan AJ (2012) Translaminar fracture toughness testing of
composites: A review. Polym Test.
Lakshminarayana HV (1984) A Symmetric Rail Shear Test for Mode II Fracture Toughness
(GIIC) of Composite Materials—Finite Element Analysis. J Compos Mater 18:227–238.
Lambert J, Chambers AR, Sinclair I, Spearing SM (2011) 3D damage characterisation and the
role of voids in the fatigue of wind turbine blade materials. Compos Sci Technol.
Lee S, Munro M (1986) Evaluation of in-plane shear test methods for advanced composite
materials by the decision analysis technique. Composites 17:13–22.
Lemaitre J, Chaboche J-L, Benallal A, Desmorat R (2009) Mécanique des matériaux solides -
3ème édition. Dunod
Lessard LB, Eilers OP, Shokrieh MM (1997) Modification of the three-rail shear test for
composite materials under static and fatigue loading. Compos Mater Test Des Ed Am Soc
Test Mater Thirteeth:217–233.
Lessard LB, Eilers OP, Shokrieh MM (1995) Testing of In-Plane Shear Properties under Fatigue
Loading. J Reinf Plast Compos 14:965–987.
Lisle T, Bouvet C, Pastor ML, et al. (2012) Damage assessment of thin woven composite
subjected to quasi-static tensile loading using infrared thermography.
254
Lomov SV, Ivanov DS, Truong TC, et al. (2008) Experimental methodology of study of damage
initiation and development in textile composites in uniaxial tensile test. Compos Sci
Technol 68:2340–2349.
Lomov SV, Verpoest I, Peeters T, et al. (2003) Nesting in textile laminates: geometrical
modelling of the laminate. Compos Sci Technol 63:993–1007.
Maillet I, Michel L, Rico G, et al. (2012) A new test methodology based on structural resonance
for mode I fatigue delamination growth in an unidirectional composite. Compos. Struct.
Maire JF, Chaboche JL (1997) A new formulation of continuum damage mechanics (CDM) for
composite materials. Aerosp Sci Technol 1:247–257.
Mandell JF (1975) Fatigue crack propagation rates in woven and nonwoven fiber glass
laminates. Astm Spec Tech Publ 515–527.
Mandell JF, Reed RM, Samborsky DD (1992) Fatigue of fiberglass wind turbine blade materials.
Sandia National Laboratories
Mandell JF, Samborsky DD, Sutherland HJ (1999) Effects of materials parameters and design
details on the fatigue of composite materials for wind turbine blades. Ewec-Conf.-. pp
628–633
Manger CIC, Ogin SL, Smith PA, Greaves RP (1997) Damage Development in plain weave GFRP.
Murray L. Scott, Gold Coast, Queensland, Australia, pp 58–66
Manjunatha CM, Taylor AC, Kinloch AJ, Sprenger S (2010) The Tensile Fatigue Behavior of a
GFRP Composite with Rubber Particle Modified Epoxy Matrix. J Reinf Plast Compos
29:2170–2183.
Mao H, Mahadevan S (2002) Fatigue damage modelling of composite materials. Compos Struct
58:405–410.
Marín JC, Barroso A, París F, Cañas J (2009) Study of fatigue damage in wind turbine blades. Eng
Fail Anal 16:656–668.
Melro AR, Camanho PP, Andrade Pires FM, Pinho ST (2012) Numerical simulation of the non-
linear deformation of 5-harness satin weaves. Comput Mater Sci 61:116–126.
Metcalfe AG, Schmitz GK (1972) Mechanism of stress corrosion in E-glass filaments. Glass
Technol 13:5–16.
255
Miot S (2009) Rupture de structures composites stratifiées sous chargements statique et de
fatigue. Thèse de doctorat.
Miot S, Hochard C, Lahellec N (2010) A non-local criterion for modelling unbalanced woven ply
laminates with stress concentrations. Compos Struct 92:1574–1580.
Moës N, Dolbow J, Belytschko T (1999) A FINITE ELEMENT METHOD FOR CRACK GROWTH
WITHOUT REMESHING. Int J Numer Methods Eng 131–150.
Mohseni Shakib SM, Li S (2009) Modified three rail shear fixture (ASTM D 4255/D 4255M) and
an experimental study of nonlinear in-plane shear behaviour of FRC. Compos Sci Technol
69:1854–1866.
Mourot G (2007) Effet de taille lié à l’endommagement des mortiers. Thèse de doctorat.
Bordeaux I
Murakami Y, Fukushima Y, Toyama K, Matsuoka S (2008) Fatigue crack path and threshold in
Mode II and Mode III loadings. Eng Fract Mech 75:306–318.
Muskhelishvili NI, Radok JRM (1953) Some basic problems of the mathematical theory of
elasticity. Cambridge Univ Press
Naderi M, Kahirdeh A, Khonsari MM (2012) Dissipated thermal energy and damage evolution of
Glass/Epoxy using infrared thermography and acoustic emission. Compos Part B Eng
43:1613–1620.
Naik NK, Shembekar PS (1992) Elastic behavior of woven fabric composites: I—Lamina
analysis. J Compos Mater 26:2196–2225.
Naik RA, Patel SR, Case SW (2001) Fatigue Damage Mechanism Characterization and Modeling
of a Woven Graphite/Epoxy Composite. J Thermoplast Compos Mater 14:404–420.
Navarro P, Aubry J, Marguet S, et al. (2012) Experimental and numerical study of oblique
impact on woven composite sandwich structure: Influence of the firing axis orientation.
Compos Struct 94:1967–1972.
Nishikawa Y, Okubo K, Fujii T, Kawabe K (2006) Fatigue crack constraint in plain-woven CFRP
using newly-developed spread tows. Int J Fatigue 28:1248–1253.
Nuismer RJ, Whitney JM (1975) Uniaxial failure of composite laminates containing stress
concentrations. Fract Mech Compos Astm Stp 593:117–142.
Odegard GM, Bandyopadhyay A (2011) Physical aging of epoxy polymers and their composites.
J Polym Sci Part B Polym Phys 49:1695–1716.
Oliver J (1996) Modelling Strong Discontinuities in Solid Mechanics Via Strain Softening
Constitutive Equations. Part 2: Numerical Simulation. Int J Numer Methods Eng 39:3601–
3623.
256
Oliver J, Huespe AE, Sánchez PJ (2006) A comparative study on finite elements for capturing
strong discontinuities: E-FEM vs X-FEM. Comput Methods Appl Mech Eng 195:4732–
4752.
Osada T, Nakai A, Hamada H (2003) Initial fracture behavior of satin woven fabric composites.
Compos Struct 61:333–339.
Oterkus E, Barut A, Madenci E (2010) Damage growth prediction from loaded composite
fastener holes by using peridynamic theory. Proc.
Pandita SD, Huysmans G, Wevers M, Verpoest I (2001) Tensile fatigue behaviour of glass plain-
weave fabric composites in on- and off-axis directions. Compos Part Appl Sci Manuf
32:1533–1539.
Parks DM (1974) A stiffness derivative finite element technique for determination of crack tip
stress intensity factors. Int J Fract 10:487–502.
Pettersson KB, Neumeister JB, Gamstedt EK, Öberg H (2006) Stiffness reduction, creep, and
irreversible strains in fiber composites tested in repeated interlaminar shear. Compos
Struct 76:151–161.
Philippidis TP, Vassilopoulos AP (2000) Fatigue design allowables for GRP laminates based on
stiffness degradation measurements. Compos Sci Technol 60:2819–2828.
Plank R, Kuhn G (1999) Fatigue crack propagation under non-proportional mixed mode
loading. Eng Fract Mech 62:203–229.
Post NL, Case SW, Lesko JJ (2008) Modeling the variable amplitude fatigue of composite
materials: A review and evaluation of the state of the art for spectrum loading. Int J
Fatigue 30:2064–2086.
Quaresimin M, Susmel L, Talreja R (2010) Fatigue behaviour and life assessment of composite
laminates under multiaxial loadings. Int J Fatigue 32:2–16.
Wisnom MR (1995) The effect of fibre rotation in ±45° tension tests on measured shear
properties. Composites 26:25–32.
Ramani SV, Williams DP (1977) Notched and unnotched fatigue behavior of angle-ply
graphite/epoxy composites. Fatigue Filamentary Compos Mater Astm Stp 636:27–46.
Rakotoarisoa C, Laurin F, Hirsekorn M, et al. (2011) Prévision de la durée de vie en fatigue des
composites à matrice organique tissés interlock. Comptes-Rendus 17èmes Journées Natl.
Sur Compos. Jnc17.
257
Reifsnider KL, Talug A (1980) Analysis of fatigue damage in composite laminates. Int J Fatigue
2:3–11.
Rice JR (1968) A path independent integral and the approximate analysis of strain
concentration by notches and cracks. DTIC Document
Römelt P, Cunningham PR (2012) A multi-scale finite element approach for modelling damage
progression in woven composite structures. Compos Struct 94:977–986.
Schellekens JCJ, De Borst R (1991) Numerical simulation of free edge delamination in graphite-
epoxy laminates under uniaxial tension. Compos Struct 647–657.
Searles K, Odegard G, Kumosa M (2001) Micro- and mesomechanics of 8-harness satin woven
fabric composites: I — evaluation of elastic behavior. Compos Part Appl Sci Manuf
32:1627–1655.
Shin C., Wang C. (2002) A comparison of as-fatigue and re-consolidation residual properties for
notched quasi-isotropic [0/45/90/−45]2S and cross-ply [0/90]4S AS4/PEEK composite
laminates. Compos Part Appl Sci Manuf 33:1519–1528.
Shindo Y, Inamoto A, Narita F (2005) Characterization of Mode I fatigue crack growth in GFRP
woven laminates at low temperatures. Acta Mater 53:1389–1396.
Shindo Y, Takano S, Horiguchi K, Sato T (2006) Cryogenic fatigue behavior of plain weave
glass/epoxy composite laminates under tension–tension cycling. Cryogenics 46:794–
798.
Sih GC (1974) Strain-energy-density factor applied to mixed mode crack problems. Int J Fract
10:305–321.
Sih GC, Paris PC, Irwin GR (1965) On cracks in rectilinearly anisotropic bodies. Int J Fract Mech
1:189–203.
Silling SA (2000) Reformulation of elasticity theory for discontinuities and long-range forces. J
Mech Phys Solids 48:175–209.
Silling SA, Askari E (2005) A meshfree method based on the peridynamic model of solid
mechanics. Comput Struct 83:1526–1535.
258
Smith EW, Pascoe KJ (1989) Biaxial Fatigue of a Glass-Fibre Reinforced Composite. I. Fatigue
and Fracture Behaviour. Mech Eng Publ Biaxial Multiaxial Fatigue 367–396.
Sneddon In (1946) The distribution of stress in the neighbourhood of a crack in an elastic solid.
Proc R Soc Lond Ser Math Phys Sci 187:229–260.
Song D-Y, Otani N (1998) Approximate estimation of fatigue strength of polymer matrix
composites by material properties. Mater Sci Eng 254:200–206.
Steeves CA, Fleck NA (2004) Collapse mechanisms of sandwich beams with composite faces and
a foam core, loaded in three-point bending. Part II: experimental investigation and
numerical modelling. Int J Mech Sci 46:585–608.
Tada H, Paris PC, Irwin GR (1973) The stress analysis of cracks handbook. Del Research
Corporation
Talreja R (1985) Transverse cracking and stiffness reduction in composite laminates. J Compos
Mater 19:355–375.
Tan S, Kim R (1988) Fracture of composite laminates containing cracks due to shear loading.
Exp Mech 28:364–372.
Tenchev RT, Falzon BG (2006) A pseudo-transient solution strategy for the analysis of
delamination by means of interface elements. Finite Elem Anal Des 42:698–708.
Thiruppukuzhi SV, Sun CT (2001) Models for the strain-rate-dependent behavior of polymer
composites. Compos Sci Technol 61:1–12.
Toubal L, Karama M, Lorrain B (2006) Damage evolution and infrared thermography in woven
composite laminates under fatigue loading. Int J Fatigue 28:1867–1872.
Trappe V, Günzel S, Jaunich M (2012) Correlation between crack propagation rate and cure
process of epoxy resins. Polym Test 31:654–659.
259
Turon A, Davila CG, Camanho PP, Costa J (2007) An engineering solution for mesh size effects in
the simulation of delamination using cohesive zone models. Eng Fract Mech 74:1665–
1682.
Van Dreumel WHM, Kamp JLM, TU Delft: Aerospace Engineering (1977) Non Hookean
behaviour in the fibre direction of carbonfibre composites and the influence of fibre
waviness on the tensile properties.
Van Paepegem W, Degrieck J (2001) Experimental set-up for and numerical modelling of
bending fatigue experiments on plain woven glass/epoxy composites. Compos Struct
51:1–8.
Van Paepegem W, Degrieck J (2003) Modelling damage and permanent strain in fibre-
reinforced composites under in-plane fatigue loading. Compos Sci Technol 63:677–694.
Van Paepegem W, Degrieck J (2002) Effects of Load Sequence and Block Loading on the Fatigue
Response of Fiber-Reinforced Composites. Mech Adv Mater Struct 9:19–35.
Van Paepegem W, Degrieck J (2002) A new coupled approach of residual stiffness and strength
for fatigue of fibre-reinforced composites. Int J Fatigue 24:747–762.
Van Paepegem W, Degrieck J, De Baets P (2001) Finite element approach for modelling fatigue
damage in fibre-reinforced composite materials. Compos Part B Eng 32:575–588.
Wang CM, Shin CS (2002) Residual properties of notched [0/90]4s AS4/PEEK composite
laminates after fatigue and re-consolidation. Compos Part B Eng 33:67–76.
Whitney JM, Stansbarger DL, Howell HB (1971) Analysis of the rail shear test applications and
limitations. J Compos Mater 5:24–34.
Wisnom MR (2010) Modelling discrete failures in composites with interface elements. Compos
Part Appl Sci Manuf 41:795–805.
Xiao J, Bathias C (1994) Fatigue behaviour of unnotched and notched woven glass/epoxy
laminates. Compos Sci Technol 50:141–148.
Yaniv G, Daniel IM, Lee JW (1989) Method for Monitoring In-Plane Shear Modulus in Fatigue
Testing of Composites,. Northwestern Univ Evanston IL Dept Of Civil Engineering
Zrida M, Laurent H, Rio G, et al. (2009) Experimental and numerical study of polypropylene
behavior using an hyper-visco-hysteresis constitutive law. Comput Mater Sci 45:516–
527.
260
Liste des figures
261
Figure II.16 – Comparaison de la tenue en fatigue d’un stratifié quasi-isotrope de taffetas carbone/époxy à différents
rapports de charge [Kawai12]...................................................................................................................................................................... 33
Figure II.17 – Comportement du tissu en traction dans l’axe des renforts. (a) : comportement sous chargement
monotone, (b) : évolution de l’endommagement, (c) : évolution de la déformation anélastique [Bizeul09]. ....... 34
Figure II.18 – Géométrie des éprouvettes de traction utilisées ................................................................................................................... 36
Figure II.19 – Clichés MEB d’échantillons d’éprouvettes de traction ayant subi différentes déformations maximales ε :
A gauche du schéma de l’éprouvette, les coupes « planes », à droite les coupes de « profil » de normale x, en
bas, la courbe de comportement en traction de l’éprouvette. Les points de couleurs correspondent aux
déformations maximales atteintes par les éprouvettes sur lesquelles les clichés ont été réalisés. ........................... 37
Figure II.20 – Résultats d’un essai de traction filmé par caméra infrarouge sur la zone encadrée de l’éprouvette en bas
à gauche. (a)-(d), vue de l’évolution de la température par rapport à la température initiale, ΔT, à différents
instants. (e), évolution temporelle de ΔT moyennée sur 4 zones différentes définies par des rectangles de
couleur en (a)........................................................................................................................................................................................................ 39
Figure II.21 – Evolution des sources de chaleur cumulées dans le temps (a-d) et agrandissement de ces champs (e-h)
sur lesquels le motif de l’armure satin de 8 a été superposé. ....................................................................................................... 40
Figure II.22 – Evolution du module d’Young en fatigue sur éprouvettes sens chaîne (C) et sens trame (T). Résultats
d’essais pour différentes déformations maximales εmax imposées. ............................................................................................ 44
Figure II.23 – Evolution du module d’Young en fonction du nombre de cycles normé par le nombre de cycles à rupture
N/Nf pour des essais de fatigue dans le sens chaîne à différents niveaux de déformation. ........................................... 45
Figure II.24 – Endommagement d1f mesuré en fonction de la déformation imposée en fatigue pour différentes valeurs
du nombre de cycles. L’endommagement à rupture de chaque essai est également reporté, le pointillé noir en
représente la moyenne. ................................................................................................................................................................................... 46
Figure II.25 – Durée de vie en fatigue à déformation imposée pour un rapport de charge R de 1/3........................................ 47
Figure II.26 – Clichés MEB d’une éprouvette de traction chargée en fatigue à une déformation maximale de 8.10 -3
jusqu’à rupture. ................................................................................................................................................................................................... 48
Figure II.27 – (a) illustration des repères structure (x,y) et matériau (1,2) sur un essai de traction sur drapage [±45].
(b), illustration de la rigidification par rotation des fibres à haute déformation lors d’essais de traction sur
stratifié [±45]2s pour différentes éprouvettes Carbone/PEEK [Kellas92] (L’auteur explique la dispersion à
déformation élevée par le procédé de fabrication). ........................................................................................................................... 50
Figure II.28 – Exemples d’essais permettant l’étude du cisaillement plan des stratifiés [Tarnopol’skii00]. ........................ 52
Figure II.29 – Comportement pour un essai de traction cyclé à 45° d’un tissu 2,5D carbone/époxy. (a), courbe
contrainte/déformation, (b), détail de l’évolution des déformations résiduelles, (c), recouvrement des
déformations résiduelles au cours du temps [Couégnat08]. ......................................................................................................... 53
Figure II.30 –Synthèse des durées de vie en fatigue sur des tissus à base de fibres de verre pour des drapages [±45] n
[Bizeul09]. Une loi de type Basquin a été identifiée sur les essais de [Pandita01]. ........................................................... 54
Figure II.31 - Perte de module d’un stratifié de tissu verre/époxy [±45] 5 sollicité en traction ondulée à un niveau
maximum de 50% de la contrainte de rupture et à différentes fréquences [Pandita01]. ............................................... 55
Figure II.32 – Schéma de principe de l’essai de cisaillement standard. L’éprouvette est maintenue encastrée sur les 3
rails. Les deux rails externes sont fixes tandis que le troisième coulisse selon l’axe y et génère un état de
cisaillement pur (aux effets de bord près) sur les deux zones visibles de l’éprouvette. .................................................. 56
Figure II.33 – Simulation EF des champs de cisaillement τ12 (a) et de contrainte σ11 (b) dans un des deux rectangles
sollicités en cisaillement [DeBaere08a]................................................................................................................................................... 57
Figure II.34 – Géométrie de l’éprouvette de cisaillement. .............................................................................................................................. 58
Figure II.35 – Données expérimentales des essais préliminaires sur montage non modifié : (a) charge en fonction de la
distorsion angulaire pour le premier et le dernier cycle de chargement, (b) signaux des premiers cycles, et
signal écrêté en fin d’essai. ............................................................................................................................................................................. 59
262
Figure II.36 – Solution technologique proposée par [DeBaere08b] (a), et [MohseniShakib09] (b) de montage de
cisaillement au rail autorisant la fatigue à rapport de charge négatif. ..................................................................................... 60
Figure II.37 – Schéma CAO du montage de cisaillement au rail modifié. (a) : montage entier, le dispositif
d’encastrement de l’éprouvette n’est pas représenté, (b) : vue en coupe du détail de la liaison ponctuelle
bidirectionnelle.................................................................................................................................................................................................... 61
Figure II.38 – Clichés du montage installé sur une machine de traction/compression. (a) : vue de face, (b) : vue de la
face arrière montrant l’extensomètre et une bride, (c) : détail de la liaison ponctuelle bidirectionnelle. .............. 62
Figure II.39 – Vue CAO de l’extensomètre adapté à la mesure du cisaillement. .................................................................................. 63
Figure II.40 – Résultats de stéréo-corrélation pour une éprouvette sous chargement quasi-statique à une déformation
homogène de cisaillement γxy = 0,0125. (a) : déplacement hors plan de l’éprouvette, (b) (c) (d) : champs de
déformations γxy, εyy et εxx respectivement. ............................................................................................................................................ 64
Figure II.41 – Déformations longitudinales, transversales et de cisaillement le long des lignes bleue et rouge pour une
distorsion angulaire γ = 0,0125. .................................................................................................................................................................. 65
Figure II.42 – (a), résultats d’essai de cisaillement au rail et de traction à ±45° quasi-statique monotone. (b), pour
l’essai de cisaillement au rail, la rupture est initiée dans les coins des zones sollicitées, sièges de contraintes de
traction/compression parasites. ................................................................................................................................................................. 66
Figure II.43 – Réponse contrainte/déformation d’éprouvettes [0/90]4 en cisaillement pour deux séquences de
chargement différentes. Les droites en pointillés matérialisent les pentes correspondant aux modules tangent
(en rouge) et sécant (en bleu) mesurés pour la dernière charge. ............................................................................................... 67
Figure II.44 – Evolution des grandeurs GT, GS, et γp en fonction de la déformation maximale de cisaillement appliquée
pour deux éprouvettes « 1 » et « 2 ». ......................................................................................................................................................... 69
Figure II.45 – Evolution du module de cisaillement tangent en fonction du nombre de cycles pour une distorsion
maximale γmax=0,012. La ligne en tiret épais matérialise l’endommagement obtenu en sollicitation quasi-
statique c'est-à-dire après 1 cycle. ............................................................................................................................................................. 70
Figure II.46 – Visualisation de quelques cycles de fatigue à déformation imposée sur une éprouvette avant sa
visualisation au MEB. ........................................................................................................................................................................................ 71
Figure II.47 – Clichés MEB d’une éprouvette de cisaillement après 2.10 3 cycles à une distorsion maximale de 0,024
puis 2.103 cycles à une distorsion maximale de 0,036. .................................................................................................................... 71
Figure II.48 – Schématisation du cisaillement d’un échantillon de tissu dans le repère (x,y) se traduisant en
traction/compression dans un repère tourné de 45°. ...................................................................................................................... 72
Figure II.49 – Courbe contrainte/déformation d’un échantillon de résine époxy en traction (a) et en cisaillement (b)
[Fiedler01]. ............................................................................................................................................................................................................ 73
Figure II.50 - (a)-(c), vue de l’évolution de la température par rapport à la température initiale, ΔT à différents
instants. (d), évolution temporelle de ΔT moyennée sur la zone définie par un rectangle en (a) comparée à
l’évolution de la contrainte en fonction du temps. ............................................................................................................................. 74
Figure II.51 – Géométrie de l’éprouvette de traction utilisée pour l’étude des couplages traction-traction, et traction-
cisaillement............................................................................................................................................................................................................ 78
Figure II.52 – Evolution du module d’Young en fatigue pour les éprouvettes de grande taille E1 à E4, et leur
comparaison avec des éprouvettes de drapage équivalent de taille standard. .................................................................... 79
Figure II.53 – Faciès de rupture et découpage de 2 éprouvettes de traction de grande taille. ..................................................... 80
Figure II.54 – Géométrie des éprouvettes de couplage (ε1 – ε2) (a) et (ε1 – γ12) (b)........................................................................... 81
Figure II.55 – (a), comportement contrainte/déformation pour 2 éprouvettes saines et 2 éprouvettes préalablement
endommagées dans la direction 1 à hauteur de 19%. (b), évolution du module d’Young tangent rapporté au
module d’Young initial, en fonction de la déformation maximale atteinte pour les 6 essais réalisés. ...................... 82
Figure II.56 - Evolution du module d’Young tangent en fatigue, comparaison entre éprouvettes saines (préfixes E5) et
préalablement endommagées dans le sens 1 à hauteur de 19% (préfixe E4) pour différents niveaux de
déformation maximale. .................................................................................................................................................................................... 84
263
Figure II.57 – Comparaison entre les réponses contrainte/déformation des éprouvettes de couplage en cisaillement, et
d’une éprouvette de taille standard sous chargement monotone. ............................................................................................. 86
Figure II.58 – Illustration du découpage des éprouvettes de couplage (γ12 - ε1) après essai de fatigue en cisaillement. 87
Figure II.59 – Comparaison des réponses contrainte/déformation en traction entre des éprouvettes saines et
préalablement sollicitées en fatigue de cisaillement. ........................................................................................................................ 88
Figure III.1 - Eprouvettes utilisées pour la mesure de la ténacité translaminaire. (a), compact tension (CT), (b), flexion
4 points (4PB), (c), double edge notched tension, (d), extended compact tension (ECT), (e), center notched
tension, (f), single edge notched tension [Laffan12]. ........................................................................................................................ 93
Figure III.2 – Effet d’un trou sur la courbe S-N de stratifiés tissés et non-tissés à base de fibres de carbone [Bishop89].
..................................................................................................................................................................................................................................... 94
Figure III.3 – (a), Courbes S-N d'éprouvettes saines et trouées. (b), C-scan d’une éprouvette trouée après différents
nombres de cycles de fatigue à 0,85σr [Wang02]. .............................................................................................................................. 95
Figure III.4 – Résistance résiduelle en fonction du nombre de cycles de fatigue subis par l’éprouvette [Bizeul09]. ........ 96
Figure III.5 - Augmentation de la résistance résiduelle post-fatigue d’éprouvettes trouées en fonction du niveau de
chargement appliqué et du nombre de cycles de fatigue effectués. Matériau AS4/PEEK drapage [0 ;90]4s
[Wang02]. ............................................................................................................................................................................................................... 97
Figure III.6 – (a) Evolution de la longueur de coupure en fonction du nombre de cycles pour 6 éprouvettes identiques
et pour la même sollicitation cyclique. (b) Vitesse de propagation en fonction de la longueur de coupure pour
les mêmes essais. Les éprouvettes qui initient vite propagent vite et réciproquement [Bizeul09]. .......................... 98
Figure III.7 – Identification d’une loi de Paris sur des essais de propagation à différentes températures [Shindo05]. 100
Figure III.8 – Identification d’une loi de Paris pour 6 essais sur éprouvettes [0/90] 2 réalisés selon les mêmes
conditions [Bizeul09]..................................................................................................................................................................................... 100
Figure III.9 – Eprouvette CT de propagation de coupure en fatigue d’un stratifié verre/époxy mettant en évidence la
zone de blanchiment de résine en pointe d’entaille (a), ainsi que le trajet de la coupure influencé par le motif
de tissage (b)...................................................................................................................................................................................................... 101
Figure III.10 - (a), éprouvette de propagation de résine époxy. (b), résultat de propagation de type loi de Paris
[Trappe12]. ......................................................................................................................................................................................................... 102
Figure III.11 - Illustration d'une propagation de coupure sous chargement de traction, l’apparition d’une longueur
caractéristique de « ligament » et la présence de « subcracks » parallèles à l’axe de chargement à intervalles
réguliers [Mandell75]. ................................................................................................................................................................................... 103
Figure III.12 - Essai Compact Tension Shear utilisé pour la propagation translaminaire en mode II et mixte I+II par
[Laffan13]. ........................................................................................................................................................................................................... 104
Figure III.13 – Géométries d'éprouvettes utilisées par [Schmitt02] (a), éprouvette « écharde », (b), éprouvette à
« fissure de bord », (c), éprouvette à « fissure centrale ». ............................................................................................................ 105
Figure III.14 – Principe de pilotage et d’automatisation de la prise de vue des essais de propagation en fatigue. ......... 107
Figure III.15 – Procédure de la routine de dépouillement par différence d’images. ...................................................................... 108
Figure III.16 – Interface utilisateur de la routine de soustraction d’images, exemple d’utilisation avec deux images
brutes A et B en entrée et le rendu des opérations de soustraction, seuillage et contraste (image D1), puis de
filtre de Wiener (image D2). ....................................................................................................................................................................... 109
Figure III.17 – Géométrie de l’éprouvette de propagation en traction. En haut, la géométrie de base utilisée pour les
drapages [0/90]n et quasi-isotrope [0;45]s. En bas, l’adjonction d’UD carbone faisant office d’arêtier est utilisée
pour les drapages [±45]n.............................................................................................................................................................................. 110
Figure III.18 – Eprouvette à la fin de l’essai de propagation en traction. (a), zone de propagation de coupure. (b), cliché
obtenu par rétro-éclairage permettant de visualiser un endommagement le long de la coupure par contraste
avec les zones saines. ..................................................................................................................................................................................... 111
264
Figure III.19 – Evolution de la longueur de coupure en propagation de traction. Comparaison entre les éprouvettes de
2 ou 4 plis, sollicitées dans la direction chaîne ou trame. ............................................................................................................ 112
Figure III.20 - Evolution de la vitesse de propagation en fonction de la longueur de coupure pour des éprouvettes de 2
ou 4 plis sollicitées dans la direction chaîne ou trame. ................................................................................................................. 113
Figure III.21 – Evolution de la longueur de coupure (notée « a ») en fonction du nombre de cycles pour des
éprouvettes à mèches de trame travaillantes de largeur 30 mm [Bizeul09]. .................................................................... 114
Figure III.22 – Estimation du développement de l’endommagement au début d’un essai de propagation de coupure
obtenu par soustraction des images correspondant au nombre de cycles indiqué. Elles ont été obtenues pour
les mêmes paramètres de post-traitement. Le segment rouge marque la pré-entaille de 7 mm. ............................ 115
Figure III.23 – Développement de l’endommagement obtenu par soustraction des images correspondant au nombre
de cycles indiqué. A titre d’exemple la première image est issue de la différence entre les images 0,8.10 5 et
1.105 cycles. Ces images ont été obtenues pour les mêmes paramètres de post-traitement. .................................... 117
Figure III.24 – Evolution mesurée de la longueur de coupure en fonction du nombre de cycles pour 3 éprouvettes
[±45]2 [Bizeul09]. ............................................................................................................................................................................................ 118
Figure III.25 –Evolution de la vitesse de propagation en fonction de la longueur de coupure pour 3 éprouvettes [±45]2
[Bizeul09]. ........................................................................................................................................................................................................... 119
Figure III.26 – Cliché de 4 éprouvettes de drapage quasi-isotrope à la fin de l’essai de propagation en traction, obtenu
par rétro-éclairage. Les croix rouges représentent les directions privilégiées des points de tissages pour les plis
à +45° (c) et (d) ou -45° (a) et (b). ........................................................................................................................................................... 120
Figure III.27 - Motif de tissage de l'armure satin de 8 orienté dans différentes directions. En termes de directions
privilégiées par rapport à l’horizontale, les directions 0° et 90° sont équivalentes à la symétrie près
contrairement aux directions +45° et -45°.......................................................................................................................................... 120
Figure III.28 – Evolution mesurée de la longueur de coupure en fonction du nombre de cycles pour des éprouvettes à
revêtement quasi-isotrope, et faciès de rupture correspondants. .......................................................................................... 121
Figure III.29 – Evolution de la vitesse de propagation en fonction de la longueur de coupure pour des éprouvettes
quasi-isotropes et une éprouvette [90]4. ............................................................................................................................................. 122
Figure III.30 – Géométrie de l’éprouvette de propagation de coupure en cisaillement adoptée pour cette étude. La
largeur w est constante, imposée par le montage, le paramètre a0 a été fixé à 25 mm, et b/w vaut 0,2 ou 0,5
selon les cas. ....................................................................................................................................................................................................... 123
Figure III.31 – Champ de distorsion angulaire γxy et conditions limites du modèle par éléments finis de l’essai de
cisaillement au rail permettant d’estimer la rotation du rail central. .................................................................................... 124
Figure III.32 – Résultats d’essai de propagation d'un stratifié [0/90]4 entaillé sous chargement de cisaillement quasi-
statique. Différents états sont représentés correspondant à des niveaux de déformation croissants. Le champ
de déformation dans le repère (1,2) obtenu par stéréo-corrélation est représenté pour l’état 3. .......................... 125
Figure III.33 - Illustration des notations de contrainte et résistance normales en pointe d’entaille. .................................... 125
Figure III.34 – Illustration de la propagation de coupure en cisaillement pour les éprouvettes S1 et S2. A gauche, les
faciès de la pointe d’entaille avant l’essai. Au centre, l’éprouvette après 10 6 cycles, la zone d’endommagement
apparente sous la coupure a été repassée en pointillés rouges. A droite, la coupure en fin d’essai. ...................... 127
Figure III.35 - Evolution de l'endommagement en propagation sous chargement de cisaillement de l’éprouvette S2
obtenue par soustraction d’images. (a), entre image initiale et 10 6 cycles, (b) entre 106 et 5.106 cycles, (c) entre
5.106 et 12.106 cycles. .................................................................................................................................................................................... 128
Figure III.36 – Mise en évidence du trajet de coupure privilégié par les points de tissage dans le cas de la propagation
de coupure en cisaillement. (a) et (b), clichés de l’éprouvette après propagation. (c), schéma de l’armure satin
de 8. (d), vue post-mortem au MEB de l’éprouvette S2 qui met en évidence le trajet de la coupure à travers les
points de tissage et les fissures secondaires à intervalles réguliers. ...................................................................................... 129
Figure III.37 – Cliché MEB de la zone d’endommagement. On observe la présence de fissuration matricielle dans les
mèches verticales jusqu’à 10 mm à l’aplomb de la coupure. Ces fissures se propagent avec un léger angle par
265
rapport à la direction des mèches. Le grossissement montre que ces fissures traversent des fibres en se
propageant. ......................................................................................................................................................................................................... 130
Figure III.38 – Cliché MEB de la coupure entière. On aperçoit de la fissuration matricielle à 45° dans les zones riches en
résine (agrandissement bleu). La pointe de coupure est visible dans l’agrandissement rouge. .............................. 131
Figure III.39 – Evolution de l’endommagement au cours de l’essai S3 : déformation imposée de γmax = 6.10-3 pendant
12 Mcy puis γmax = 12.10-3 jusqu’à rupture. ........................................................................................................................................ 132
Figure III.40 - Résultats expérimentaux de propagation de coupure : longueur de coupure en fonction du nombre de
cycles...................................................................................................................................................................................................................... 133
Figure III.41 - Résultats expérimentaux de propagation de coupure: vitesse de propagation en fonction de la longueur
de coupure. ......................................................................................................................................................................................................... 133
Figure III.42 –Résultat de l’essai S4. (a) cliché de la coupure en fin d’essai. (b)-(e), résultats de soustraction d’images
entre les différents instants indiqués. La trace rouge marque la position de la pré-entaille. .................................... 134
Figure III.43 – Propagation de coupure en cisaillement avec drapage quasi-isotrope. ................................................................ 135
Figure III.44 - Evolution de l'endommagement de l'éprouvette QI2 par comparaison entre 2 images à des nombres de
cycles donnés. .................................................................................................................................................................................................... 136
Figure III.45 – Comparaison des vitesses de propagation des éprouvettes quasi-isotropes et de l’éprouvette S4 de
drapage [0/90]4 et éprouvées dans les mêmes conditions. ........................................................................................................ 137
Figure III.46 - Agrandissement de 3 images de post-traitement successives issues de la soustraction numérique entre
les instants N et N+ΔN pour N = (a) : 1,4.105, (b) : 1,5.105 (c) : 1,6.105 cycles et ΔN=104 cycles. La flèche rouge
représente un repère qui pointe vers le même point géométrique correspondant à la pointe de coupure
mesurée sur l’image (a). ............................................................................................................................................................................... 137
Figure IV.1 – Illustration des 3 modes de rupture, en fonction de la direction de déplacement des lèvres de la fissure.
.................................................................................................................................................................................................................................. 141
Figure IV.2 – Repère cylindrique associé à un problème de mécanique de la rupture. ................................................................ 141
Figure IV.3 – Illustration de la méthode d’avancée des fissures par déplacement des nœuds au voisinage immédiat de
la fissure [François93]................................................................................................................................................................................... 144
Figure IV.4 - Illustration de la méthode de l'intégrale J. ............................................................................................................................... 145
Figure IV.5 – Comparaison entre les prévisions de rupture et les données expérimentales en traction et en cisaillement
pour différentes longueurs d’entaille et pour un drapage [0/90/±45] s (a) et [0/90]2s (b) d’après [Tan88]. .... 147
Figure IV.6 - Illustration de la zone identifiée par [Hochard06] pour le critère non-local, utilisé pour déterminer la
rupture en bord de défaut. .......................................................................................................................................................................... 148
Figure IV.7 – Illustration de différentes échelles de modélisation [Couegnat08]. ........................................................................... 149
Figure IV.8 – (a) VER d’une mèche tenant compte de la dispersion des fibres, portion de maillage à l’échelle du motif
de tissage, (c) et (d) champ des variables d’endommagement transverse du composite en traction à deux
niveaux de charge [Melro12]. .................................................................................................................................................................... 150
Figure IV.9 – Illustration des travaux de [Couegnat08]. (a), cellule élémentaire représentative d’une mèche, (b),
maillage du tissu, (c), illustration du maillage d’une éprouvette présentant 2 entailles décalées, (d), faciès de
rupture de l’éprouvette correspondante. ............................................................................................................................................ 151
Figure IV.10 – (a), maillage fin d’une cellule élémentaire de mèche, (b), modèle binaire d’éléments poutres superposés
à un élément résine [Römelt12]. .............................................................................................................................................................. 152
Figure IV.11 – Faciès de rupture d’une éprouvette de traction [45; 90; -45; 0]s [Wisnom10]. ................................................. 154
Figure IV.12 – Illustration d’un modèle de zone cohésive [Vandellos11]. ........................................................................................... 155
Figure IV.13 – Principe de modélisation d’un essai de traction troué [Hallett09]. ......................................................................... 155
Figure IV.14 – Illustration des travaux de [McLaughlin00] selon 2 situations de délaminage et fissuration transverse
combinés. ............................................................................................................................................................................................................. 156
266
Figure IV.15 – Illustration des enrichissements des degrés de liberté d’après [Moës01] sur un maillage régulier (a), et
quelconque (b). ................................................................................................................................................................................................. 157
Figure IV.16 – Maillage d’un matériau à plusieurs phases (matrice renforcée par des inclusions) en X-FEM et
utilisation de la formulation d’ensembles de niveaux pour définir les géométries des constituants [Hettich08].
.................................................................................................................................................................................................................................. 158
Figure IV.17 - Illustration de la nécessité d’un nombre important de points de Gauss par éléments appartenant au
support des nœuds j [Asadpoure06]...................................................................................................................................................... 159
Figure IV.18 – Notations utilisées et grille de particules [Askari06]...................................................................................................... 160
Figure IV.19 – Fonction de comportement des liens [Silling05]. ............................................................................................................. 160
Figure IV.20 – Faciès de rupture obtenus par modélisation péridynamique d’un essai de traction présentant une
entaille centrée. De gauche à droite, le nombre de plis à 0° par rapport aux plis à 90° décroit [Askari06]. ....... 161
Figure IV.21 – Maillage d’une mèche par des éléments surfaciques séparés par des éléments de rupture. ...................... 163
Figure IV.22 – Illustration d’une portion de maillage dans une direction, distendu pour une meilleure clarté. En
configuration saine, les nœuds 2 et 3 sont géométriquement confondus, de même que les nœuds 6,7,10 et 11,
et les nœuds 14 et 15. .................................................................................................................................................................................... 163
Figure IV.23 – Liaison chaîne/trame réalisée par 4 éléments d’interface représentés par des ressorts (maillage
distendu pour la compréhension). .......................................................................................................................................................... 164
Figure IV.24 – Schéma du maillage d’un pli, distendu pour la compréhension. En pratique, le maillage est plan, et les
éléments chaîne et trame sont initialement confondus. Les éléments d’interface relient 2 nœuds initialement
confondus, et les éléments de rupture sont de surface nulle. .................................................................................................... 164
Figure IV.25 – Schéma pour le calcul des constantes élastiques E1 et E2. ............................................................................................ 166
Figure IV.26 – Illustration de la différence entre module tangent Etan et module sécant Esécant utilisé dans la définition
de l’endommagement. ................................................................................................................................................................................... 170
Figure IV.27 – Résultats expérimentaux de l’endommagement en fonction de la déformation de fatigue et du nombre
de cycles obtenus en considérant le module tangent..................................................................................................................... 171
Figure IV.28 – Elément d’interface utilisé pour lier les deux directions de maillage entre les nœuds N1 et N2. (a) : cas où
le maillage est orienté suivant le repère structure (Oxy). (b) : cas général. ......................................................................... 172
Figure IV.29 – Représentation d’une mèche en traction pour établir la relation raideur – endommagement, et
réduction en un problème 1D sur un ensemble de deux éléments surfaciques superposés liés par deux
éléments ressorts............................................................................................................................................................................................. 173
Figure IV.30 – Evolution de l’endommagement d en fonction de la raideur des ressorts k. ....................................................... 175
Figure IV.31 – Comparaison des résultats expérimentaux d’endommagement avec le modèle identifié pour des
nombres de cycles entre 103 et 106. ....................................................................................................................................................... 176
Figure IV.32 – Comparaison de la loi de raideur des ressorts en statique et en fatigue en fonction de la déformation
longitudinale de fatigue. ............................................................................................................................................................................... 177
Figure IV.33 – Endommagement d en fonction de la déformation de fatigue ε et du nombre de cycles N. ......................... 178
Figure IV.34 – Visualisation de la zone d’interface représentée par l’élément d’interface marqué par un ressort rouge.
.................................................................................................................................................................................................................................. 179
Figure IV.35 – Zone de prise en compte de la déformation dans le calcul des raideurs kx et ky de l’élément d’interface
repassé en rouge. ............................................................................................................................................................................................. 180
Figure IV.36 – Exemple de l’évolution de l’endommagement pour un élément d’interface sous la séquence de
chargement décrite plus haut. Les numéros en noir indiquent l’une des 4 phases de la séquence. ....................... 183
Figure IV.37 – Comparaison entre la loi d’endommagement identifiée et les mesures expérimentales
d’endommagement calculées par rapport à la valeur G12 = 3500 MPa.................................................................................. 184
Figure IV.38 - Comparaison entre la loi de pseudo-plasticité identifiée et un essai de cisaillement. ..................................... 185
Figure IV.39 – Schéma d’un élément de rupture. ............................................................................................................................................. 186
267
Figure IV.40 – Loi de Basquin prise en compte dans la modélisation et identifiée sur des essais présentant des
concentrations de contraintes faibles, comparée à la loi de Basquin conservative présentée en [Link]. ........... 188
Figure IV.41 – Elément de rupture, et variables de Miner affectées aux points de Gauss associés. ........................................ 189
Figure IV.42 – Illustration de la propagation de l’endommagement suite à la rupture d’une mèche. Les raideurs des
éléments d’interface de la zone repassée en rouge sont annulées suite à la rupture de la mèche via l’élément
rompu. ................................................................................................................................................................................................................... 191
Figure IV.43 – Intersection entre deux mèches matérialisée par le losange noir. ........................................................................... 192
Figure IV.44 – Représentation des 8 nœuds supports de l’élément d’interface 0-45 entre les éléments (a) et (b), et
définition des paramètres géométriques associés. ......................................................................................................................... 193
Figure IV.45 –Schéma de support pour le calcul de la matrice d’interpolation I. ............................................................................ 194
Figure IV.46 – Distribution du moment sur les 4 nœuds. ............................................................................................................................ 196
Figure IV.47 – Fonction de chargement utilisée dans la modélisation. ................................................................................................. 199
Figure IV.48 - Algorithme simplifié d'une simulation d’un essai de propagation. Les notations utilisées sont en
Tableau IV.1, lorsque le numéro d’itération i est omis, la variable correspondante est prise au pas de temps
convergé. .............................................................................................................................................................................................................. 201
Figure IV.49 – Exemple de maillage utilisé pour l’étude de la propagation de coupure. Le revêtement d’un tronçon de
pale de drapage [±45]n est représenté par des éléments membranes classiques. Seule une petite zone d’étude
est maillée avec le maillage spécifique à l’étude de la propagation. ....................................................................................... 203
.1 – Géométrie des maillages utilisés pour la simulation des essais de propagation en traction sur des
Figure V
drapages [0/90]n. (a), éprouvette de largeur 30 mm, (b), éprouvette de largeur 50 mm. .......................................... 208
Figure V.2 – Faciès de rupture en fin de simulation de la propagation en traction sur une éprouvette [0/90] 4 de largeur
50 mm. En bas, représentation des éléments de rupture seuls. ............................................................................................... 209
.3 – Déformation longitudinale dans les mèches orientées à 0° et 90° pour la simulation d’un essai de
Figure V
propagation en traction sur éprouvette [0/90]4 de largeur 50 mm, après une propagation d’environ 15 mm
(déformée x3). ................................................................................................................................................................................................... 210
Figure V.4 – Champ des variables d’endommagement dans les directions du repère d’orthotropie calculées à partir des
raideurs des éléments d’interface pour la simulation d’un essai de propagation en traction sur éprouvette
[0/90]4 de largeur 50 mm, après une propagation d’environ 15 mm. Les axes représentent les coordonnées
spatiales en mm. ............................................................................................................................................................................................... 212
.5 – Images obtenues par soustraction d’images de la propagation en traction sur une éprouvette [0/90] 4
Figure V
après environ 15 mm de propagation. .................................................................................................................................................. 213
Figure V.6 – Evolution de la vitesse de propagation de coupure en traction obtenue pour deux éprouvettes [0/90]4 de
largeur 50 mm et pour la simulation correspondante. ................................................................................................................. 214
Figure V.7 – Evolution de la vitesse de propagation de coupure en traction obtenue pour deux éprouvettes [0/90]2 de
largeur 30 mm présentant des nombres de cycles d’initiation différents et pour les simulations
correspondantes. ............................................................................................................................................................................................. 215
Figure V.8 – Visualisation des éléments de rupture à l’issue de la simulation d’un essai de propagation en traction sur
éprouvette [0/90]2 de largeur 30 mm pour deux nombres de cycles d’initiation. .......................................................... 215
Figure V.9 – Déformation longitudinale dans les mèches orientées à 90° pour la simulation d’un essai de propagation
en traction sur éprouvette [0/90]2 de largeur 30 mm après une propagation d’environ 10 mm, pour deux
nombres de cycles d’initiation................................................................................................................................................................... 216
.10 – Champ des variables d’endommagement dans les directions du repère d’orthotropie calculées à partir
Figure V
des raideurs des éléments d’interface sur des éprouvettes [0/90]2 pour 2 nombres de cycles d’initiation. Les
axes représentent les coordonnées spatiales en mm, et les flèches indiquent la position de la pointe de
coupure. ................................................................................................................................................................................................................ 217
268
Figure V.11 – Maillages utilisés pour la simulation des essais de propagation en traction sur des drapages [±45] 2. (a),
éprouvette de largeur 30 mm, (b), éprouvette de largeur 50 mm. ......................................................................................... 219
Figure V.12 - Visualisation des éléments de rupture à l’issue de la simulation d’un essai de propagation en traction sur
éprouvette [±45]2 de largeur 30 mm (a) et 50 mm (b), et comparaison avec le faciès de rupture d’une
éprouvette de largeur 50 mm (c). ............................................................................................................................................................ 220
.13 – Champ des variables d’endommagement dans les directions du repère d’orthotropie calculées à partir
Figure V
des raideurs des éléments d’interface après une propagation d’environ 10 mm sur une éprouvette [±45] 2 de
largeur 30 mm. Les axes représentent les coordonnées spatiales en mm, et les flèches indiquent la pointe de
coupure. ................................................................................................................................................................................................................ 221
.14 – Clichés obtenus par soustraction d’images de la propagation en traction sur une éprouvette [±45] 2 de
Figure V
largeur 50 mm après environ 10 mm de propagation. ................................................................................................................. 222
.15 – Déformation longitudinale dans les mèches orientées à +45° et -45° pour la simulation d’un essai de
Figure V
propagation en traction sur éprouvette [±45]2 de largeur 50 mm, après une propagation d’environ 15 mm
(déformée x1). ................................................................................................................................................................................................... 223
Figure V.16 – Champ de déformation longitudinale pour des éprouvettes [±45]2. (a) et (b) : pour largeur 30 mm après
environ 10 mm de propagation. (c), double courbure des mèches en front de fissure (éprouvette largeur 50
mm) après 30 mm de propagation environ (déformée x3). ....................................................................................................... 224
Figure V.17 –Evolution de la vitesse de propagation de coupure en traction obtenue pour trois éprouvettes [±45] 2 de
largeur 50 mm et pour la simulation correspondante. ................................................................................................................. 225
Figure V.18 – Evolution de la vitesse de propagation de coupure en traction obtenue pour trois éprouvettes [±45]2 de
largeur 30 mm et pour la simulation correspondante. ................................................................................................................. 225
Figure V.19 – Faciès de rupture des éprouvettes de drapages [±45]2. (a) : largeur 30 mm, (b), largeur 50 mm. ............ 226
Figure V.20 – Faciès de rupture en fin de simulation de la propagation en traction sur une éprouvette [0;45] s mis en
évidence par le champ de déplacement selon l’axe y pour les mèches à 90° et à 45°. ................................................... 227
Figure V.21 – (a), maillage non déformé à 0/90 et éléments d’interface 0-45, (b), maillage non déformé à ±45 avec les
mêmes éléments d’interface 0-45. .......................................................................................................................................................... 228
Figure V.22 – Déformation longitudinale dans les mèches orientées à 0°, 90°, 45° et -45° pour la simulation d’un essai
de propagation en traction sur éprouvette [0;45]s de largeur 50 mm, et après une propagation d’environ 15
mm (déformée x1). ......................................................................................................................................................................................... 229
Figure V.23 – Champ des variables d’endommagement dans les directions des repères d’orthotropie des plis 0/90 et
±45 calculées à partir des raideurs des éléments d’interface après une propagation d’environ 15 mm sur une
éprouvette [0;45]s. Les axes représentent les coordonnées spatiales en mm, et les flèches indiquent la position
de la pointe de coupure. ............................................................................................................................................................................... 230
Figure V.24 - Evolution de la vitesse de propagation de coupure en traction obtenue pour quatre éprouvettes [0;45] s et
pour la simulation correspondante. ....................................................................................................................................................... 231
Figure V.25 – Maillage utilisé pour la simulation des essais de propagation en cisaillement sur des drapages [0/90] 4.
La pointe de pré-entaille sur le maillage MEM a été repassée en rouge. .............................................................................. 232
Figure V.26 – (a), faciès de rupture en fin d’essai (b), faciès de rupture de la simulation correspondante, (c) rappel du
maillage d’un pli 0/90 mettant en évidence les éléments de rupture, (d), évolution du faciès de rupture en
cours de simulation par numérotation des éléments de rupture rompus. Les éléments de rupture bleus
(rouges) correspondent à des rupture de mèches verticales (horizontales). .................................................................... 233
Figure V.27 – Visualisation du champ de déformation longitudinale au pas de temps T correspondant à 15 ruptures de
mèches (déformée x3)................................................................................................................................................................................... 234
.28 - Evolution de la vitesse de propagation de coupure en cisaillement obtenue pour deux éprouvettes
Figure V
[0/90]4 et pour la simulation correspondante. ................................................................................................................................. 235
269
.29 – Comparaison des variables d’endommagement avec l’endommagement obtenu par soustraction
Figure V
d’images expérimentales. Les axes représentent les coordonnées spatiales en mm, et les flèches indiquent la
position de la pointe de coupure. ............................................................................................................................................................. 236
.30 – Simulation des essais de propagation en cisaillement sur des drapages [0;45]s. (a), géométrie du
Figure V
maillage (b), déformée, (c), image expérimentale. .......................................................................................................................... 237
Figure V.31 - Déformation longitudinale dans les mèches orientées à 0°, 90°, 45° et -45° pour la simulation d’un essai
de propagation en cisaillement sur éprouvette [0;45]s après 45 ruptures de mèche (déformée x2). ................... 238
.32 – Evolution de la vitesse de propagation de coupure en cisaillement obtenue pour deux éprouvettes
Figure V
[0;45]s et pour la simulation correspondante. .................................................................................................................................. 239
Figure V.33 – Champ des variables d’endommagement dans les directions des repères d’orthotropie des plis 0/90 et
±45 calculées à partir des raideurs des éléments d’interface après une propagation en cisaillement d’environ
15 mm sur une éprouvette [0;45]s. Les axes représentent les coordonnées spatiales en mm, et les flèches
indiquent la position de la pointe de coupure. .................................................................................................................................. 240
270
Propagation de coupure en fatigue sur composites tissés – Etude expérimentale et modélisation
Résumé
Les pales d’hélicoptère sont des structures composites soumises à un chargement cyclique
multiaxial, et leur criticité impose de porter une attention particulière à la tolérance aux
dommages. Leur revêtement peut potentiellement présenter des criques suite à certains
évènements (impact, défaut, foudre). Ces travaux se focalisent sur un matériau de revêtement
donné (tissu de verre) et concernent l’étude de la propagation de coupure (crique) sous
chargement cyclique. Les sollicitations de service ont amené à considérer la traction et le
cisaillement plan.
Une étude expérimentale a été menée afin d’étudier les modes d’endommagement du matériau
et sa résistance à la propagation de coupure pour différentes sollicitations (en traction et en
cisaillement) et pour les drapages les plus courants. Elle a permis de dégager les mécanismes
d’endommagement mis en jeu, et a fourni un ensemble important de propriétés matériau et de
données quantitatives de vitesse de propagation.
Elle a par ailleurs guidé vers une modélisation par éléments finis adaptée à l’architecture du
matériau, et la manière dont il se dégrade en fatigue. Ce modèle repose sur un maillage à
l’échelle de la mèche, et la prédiction de la propagation est obtenue par l’utilisation d’une courbe
de fatigue S-N. La simulation a été évaluée par comparaison des faciès de rupture, des vitesses
de propagation et de l’étendue des zones d’endommagement avec les essais réalisés sur
éprouvettes.
Mots-clés
matériau composite tissé, fissure translaminaire, propagation en fatigue, modélisation
Fatigue Crack Growth in woven composites – Experimental study and numerical modeling
Abstract
Helicopter blades consist of composite structures which have to sustain multi-axial cyclic
loading. Because of their criticality, damage tolerance has to be considered carefully. Their skin
is subjected to environmental events like impact, flaw, lightning which can cause through-the-
thickness cracks. The present work focuses on one given skin material (woven glass fabric) and
concerns the study of the through-the-thickness crack growth under cyclic loading. In-flight
loading lead to consider tension and shear.
An experimental study has been carried out to study damage in the material and its crack
growth resistance under different loadings (tension and shear) and for usual stacking sequences.
It highlighted damage mechanisms and provided an important set of material data and crack
growth speeds.
Besides, this led to a finite element approach adapted to the woven fabric architecture, and
damage feature under fatigue loading. This modeling is based on a bundle scale mesh, a semi-
discrete damage modeling and an S-N curve to predict fiber failure. Numerical simulations of
crack growth tests were carried out, and results were compared with experiments in terms of
crack direction, crack growth speed, and size of damaged area.
Keywords
woven composite materials, translaminar crack, fatigue crack growth, finite element analysis