Machine Translated by Google

Composites : Partie B 56 (2014) 660–669

Liste des contenus disponible sur ScienceDirect

Composites : Partie B

page d'accueil du journal : www.elsevier.com/locate/compositesb

Panneaux composites à rigidité variable : modélisation telle que fabriquée et son

influence sur le comportement à la rupture
b
O. Falcó a , JA Mayugo une,
, CS Lopés , N. Gascons
un
, A. Turon
un
, J. Costa
un

un
AMADE, École Polytechnique, Université de Gérone, Campus Montilivi s/n, 17071 Gérone, Espagne
b
IMDEA Materials ­ Institut d'études avancées des matériaux de Madrid, c/Eric Kandel, 2, Parque Científico y Tecnológico­Tecnogetafe, 28906 Getafe, Madrid, Espagne

informations sur l'article abstrait

Historique des articles : L'introduction de stratifiés remorqués dans les applications structurelles nécessite des méthodologies de conception
Reçu le 23 avril 2013 fiables pour prédire et simuler leur réponse mécanique. Les détails complets sur les mécanismes qui conduisent aux
Reçu sous forme révisée le 19 juillet 2013 dommages et à la défaillance de ces nouveaux panneaux ne sont pas encore connus. Dans ce travail, un outil d'analyse
Accepté le 2 septembre 2013
numérique a été développé pour des simulations structurelles dans un domaine tridimensionnel. Cet article présente une
Disponible en ligne le 11 septembre 2013
comparaison entre différentes configurations et approches de modélisation de ces panneaux. Des analyses par éléments
finis sont effectuées pour simuler la rupture du premier pli de panneaux dirigés par remorquage sous charge de traction.
Mots clés :
Les simulations montrent comment la fissuration de la matrice est affectée par les discontinuités dans l'angle des fibres
C. Modélisation informatique E.
entre les couches adjacentes et les zones sans fibres résultant des effets de fabrication ; tandis que la rupture en traction
Panneaux automatisés
à rigidité variable des fibres est directement influencée par les angles d'orientation de la fibre, dont la répartition dépend de la méthode de
coupe choisie : sur un côté ou sur les deux côtés.
2013 Elsevier Ltd. Tous droits réservés.

1. Introduction Stratifiés remorqués : le chevauchement et la méthode du remorquage [6,7]. Lors de

l'utilisation de la méthode de chevauchement, des régions plus épaisses sont créées
L'application efficace de nouveaux concepts composites, principalement dans à l'intérieur du stratifié. D'autre part, la méthode du tow­drop évite les zones de
l'ingénierie aérospatiale, nécessite non seulement un niveau avancé de technologie chevauchement en coupant les câbles lorsqu'ils se croisent, afin d'éviter toute
de fabrication, mais également des méthodologies de conception fiables, capables accumulation d'épaisseur. Lors de cette opération, de petites régions sans fibres
de prédire la réponse mécanique depuis l'apparition des dommages jusqu'à sont créées (Fig. 1b). Au cours du processus de durcissement du panneau, ces
l'effondrement structurel. Les développements récents en matière de placement zones de chute de câble deviennent des points potentiels d'initiation de fissuration
automatisé de fibres (AFP) permettent l'introduction de concepts de stratifiés non de la matrice, voire de délaminages, en raison de l'amplification locale des
conventionnels (NCL). Ces stratifiés comprennent des superpositions avec des contraintes interlaminaires [3]. Dans cet article, seule la méthode du remorquage est
séquences d'empilement dispersées [1] et des superpositions avec des fibres analysée.
dirigées par câble le long de chemins courbes ; également appelés panneaux à Dans la plupart des analyses par éléments finis (EF) précédentes des structures
rigidité variable (VSP) [2,3]. Les performances améliorées obtenues avec VSP, par placées en câble avec une épaisseur constante [8,2,7,9], les zones de chute de
rapport à leurs homologues à fibres droites, peuvent être attribuées à une capacité câble ou même les limites entre les différentes couches n'étaient pas prises en
supplémentaire à redistribuer les charges des régions traditionnellement fortement compte. En fait, les câbles sont coupés perpendiculairement aux directions des
chargées (par exemple, les environs des trous) vers les bords des panneaux, qui fibres, ce qui donne lieu à de petites zones triangulaires sans fibres. L'utilisation de
sont généralement soutenus par des éléments porteurs rigides. . câbles plus larges entraîne des zones sans fibres plus grandes, ce qui diminue
Dans le passé, différentes études ont été réalisées pour permettre aux encore la résistance du stratifié. Blom et coll. [10] ont présenté une étude théorique
concepteurs de mieux comprendre les besoins en matière d'adaptation des stratifiés [2,4]. et numérique sur l'influence de ces zones de chute de câble sur la résistance et la
Cependant, ces premières approximations ne tenaient pas compte des effets de rigidité des stratifiés à rigidité variable lorsque l'assise est coupée d'un seul côté. De
fabrication des panneaux remorqués et seules les configurations idéalisées ont été plus, Blom et al. [11] présentent une stratégie pour calculer les zones sans fibres
analysées. En conséquence, dans les panneaux réellement fabriqués, la répartition lorsque les deux côtés de la couche sont coupés. Cependant, dans un souci de
des angles des fibres sera différente des répartitions idéalisées. La distribution réelle simplification et d'efficacité, leurs analyses ont été modélisées à l'aide d'une
de l'angle des fibres provoquée par le processus de placement des fibres a été formulation de contraintes dans le plan implémentée dans les éléments de coque
décrite dans Blom et al. [5]. Il existe deux méthodes générales de fabrication pour dans un maillage régulier.

Dans les modèles FE avec maillage régulier, la discrétisation des zones de

Auteur correspondant. Tél. : +34 972418853. chute de câble dépend fortement de la densité du maillage et seule une description
Adresse e­mail : ja.mayugo@udg.edu (MYUGO). approximative de ces zones est obtenue. Résoudre

1359­8368/$ ­ voir la présentation 2013 Elsevier Ltd. Tous droits réservés. http://
dx.doi.org/10.1016/j.compositesb.2013.09.003
Machine Translated by Google

O. Falcon et coll. / Composites : Partie B 56 (2014) 660–669 661

Fig. 1. Définition du parcours pour différentes configurations de variation d'angle de fibre. Le chemin de référence est représenté à l’aide d’une ligne tiret­point.

résoudre ce problème et obtenir une approche de modélisation plus réaliste, cela un chemin de référence basé sur des arcs de cercle [15]. Dans cet article, les deux
L'article propose une méthodologie de modélisation de VSP avec un des chemins linéairement variables et circulaires sont analysés.
engrener. Contrairement aux maillages réguliers, avec une technique de maillage
structuré, une description géométrique précise des zones sans fibres est
2.1. Trajectoires de variation d'angle linéaire
atteint, même avec des mailles plus grossières, limitées par la largeur d'un
remorquer. Récemment, une méthodologie de récupération des contraintes interlaminaires utilisant
On suppose que l'angle de la fibre varie linéairement à partir de la valeur
les éléments finis tridimensionnels ont été développés par Fagiano
T0 au centre du panneau à T1 à la dimension caractéristique d,
et coll. [12] et appliqué par Díaz et al. [13] aux panneaux à rigidité variable
le long de la direction x0 . Le système de référence local x0 y0 peut être pivoté par
sans tenir compte des effets de fabrication tels que la variation de l'angle des fibres et
rapport à un système de référence global xy d'un angle /,
les zones de résine de remorquage. Les deux ouvrages affirment qu'un
comme le montre la figure 1a.
L'évaluation des contraintes interlaminaires dans les stratifiés composites multicouches
Selon cette formulation, la famille des configurations curvilignes peut être représentée
est cruciale pour la prédiction correcte du début du délaminage. De plus, Fagiano et al.
par [/ ± hT0jT1i]. Semblable à la nomenclature conventionnelle des stratifiés, un signe ±
proposer la création d'un maillage
signifie qu'il y a deux
générateur pour modéliser des stratifiés dirigés par remorquage avec un
couches adjacentes avec des variations égales et opposées de l'angle d'orientation des
maille, permettant ainsi d'inclure des zones sans fibres.
fibres. Le chemin de la courbe de référence est la trajectoire de l'AFP
Cet article étend davantage l'idée du maillage structuré et développe un cadre de
centre de la tête de la machine, et est supposé continu et antisymétrique par rapport à
pré­traitement afin d'obtenir un modèle paramétrique.
l'origine. L'équation générale pour
VSP modélisé dans un domaine tridimensionnel. Cet outil contient
la trajectoire de variation d'angle linéaire est :
possibilités de générateur de maillage structuré, permettant ainsi des
Modèles FE pour des analyses fiables des dommages progressifs et des délaminages. jx0 j
La méthodologie proposée prend en compte l'effet uðx0 Þ ¼ / þ ðT1 T0Þ þ T0 ð1Þ
j
des discontinuités de l'orientation des fibres entre les assises,
Dans les premières implémentations de VSP, seules les définitions théoriques
ainsi que l'effet des zones sans fibres. Le générateur de maillage a
unidirectionnelles ont été utilisées pour introduire l'effet de la variation de rigidité sur
été implémenté au moyen de l'interface de script Abaqus
la réponse d'un stratifié. Dans ces idéalisés
(Dassault Systèmes, France) avec le langage de programmation Python [14]. Dans cet
plis à rigidité variable, l'angle des fibres varie continuellement en fonction de la
article, nous présentons une comparaison entre différentes configurations et approches
fonction de l'une des coordonnées planaires et toutes les fibres suivent le
de modélisation de VSP sous traction.
chemin de référence (voir Fig. 2a). Dans un stratifié à plis idéaux, non seulement
chargement. Ce travail comprend également une analyse post­traitement qui nous
les zones de résine sont omises, mais aussi les discontinuités entre les assises
permet de simuler la réponse mécanique des panneaux à direction
adjacentes. Le degré de discontinuité dépend du taux à
fibres, y compris la prédiction de la rupture du premier pli.
dont l'angle de la fibre varie et sur la largeur du parcours comme
montré par Blom et al. [11,5]. Une machine de placement de fibres typique
pose généralement des parcours de 4 à 32 traits chacun, qui peuvent être coupés
2. Définitions des panneaux dirigés par remorquage individuellement et réinitialisés pendant le processus de placement [16]. Le
la largeur du remorquage est normalement de 3,175 mm ou 6,35 mm, et la largeur du parcours
Deux configurations analytiques principales ont été proposées pour définir les varie selon le nombre de traits qu'il contient. Un plus réaliste
stratifiés d'écrans plats sur la base de chemins de fibres curvilignes : le chemin d'angle Le modèle devrait prendre en compte tous ces aspects du système AFP.
des fibres à variation linéaire (Fig. 1a) et la courbure constante. Du fait de ces conditions de fabrication, l'angle de la fibre
chemin (Fig. 1b). Les travaux initiaux de Gürdal et Olmedo [6] ont abouti à sera différent dans les deux directions planes orthogonales, car
la formulation d'une définition de pli dirigé par remorquage avec un minimum dans un parcours, les traits sont parallèles les uns aux autres. En revanche,
nombre de paramètres, et principalement pour la configuration avec variation linéaire­ la construction à plis idéaux ne considère que le chemin central, alors que
oreille. Plus tard, les mêmes auteurs ont adapté ces concepts à dans un pli manufacturé avec variation d'angle linéaire, les bords de l'assise
Machine Translated by Google

662 O. Falcon et coll. / Composites : Partie B 56 (2014) 660–669

Fig. 2. Direction de la fibre sur VSP avec variation d'angle linéaire (T0 = 10, T1 = 80). L'approche « idéale » en (a) considère uniquement le chemin central, tandis que l'approche « fabricable »
en (b) prend en compte la largeur du parcours.

sont obtenus en décalant le chemin de référence dans la direction y0 y ¼ y0 qðcosu cos T0Þ ð5Þ
(voir Fig. 1a). La couverture complète d'un pli sera obtenue par plusieurs X
péché ¼ péché T0 þ ðsin T1 péché T0Þ ð6Þ
passes de la tête de placement de fibres. En conséquence, des zones d'espacement ou d
de chevauchement sont créées.

2.3. Coupe de remorquage discrète

2.2. Chemins à courbure constante

Pour qu'un pli soit entièrement recouvert de câbles de fibres, un parcours donné de
Cette définition VSP est basée sur des arcs de cercle et présente les
largeur d'origine Cw, est décalée autant de fois que nécessaire d'une distance
avantage de réaliser des assises à courbure constante, dont
We = Cw/cos(T0) mesuré le long de l'axe de décalage de référence. Ce
la fabricabilité avec une machine de placement de remorquage est plus facile à
la distance doit être constante tout au long du parcours afin de garantir
vérifier et exécuter. Le chemin de référence ici est formulé en utilisant
une couverture complète du panneau sans chevauchements ni espaces de remorquage. Dans
les mêmes paramètres /, T0 et T1 que dans la configuration avec variation d'angle linéaire­
Dans le cas général des fibres dirigées, cela n'est possible que si le module discret
oreille.
les câbles de fibre sont facilement coupés d'un ou des deux côtés du parcours.
Sur la figure 1b, q représente le rayon de braquage qui doit obéir Puis, un parcours avec un nombre variable de traits, nt, chacun de largeur tw,
la condition q P qmin, où qmin est défini comme le minimum
a sa largeur variable totale définie par Cw = nt tw. Le minimum
rayon de braquage du parcours. Ce paramètre est la contrainte de fabrication la plus
la largeur du parcours est atteinte lorsque le dernier trait restant dans un parcours
importante en termes de conception et doit être
doit être coupé pour éviter les chevauchements. La coupe discrète du câble signifie que
pris en compte dans la machine AFP pour éviter le flambement local/
l'angle maximum de la fibre n'est pas de 90 mais en réalité
sertissage des câbles intérieurs, afin d'obtenir un stratifié avec un
umax = arccos(tw/Nous).
niveau de qualité acceptable. Lorsque le rayon de braquage est trop petit,
Dans les conceptions à courbure constante proposées par Blom et al. [dix]
les câbles intérieurs du parcours ont tendance à se froisser hors du plan, créant
(voir Fig. 1b), seuls les traits d'un des côtés du parcours sont
imperfections et une réduction de la résistance du stratifié durci. Les systèmes AFP
couper. Dans ce cas, le bord supérieur du parcours est parallèle au chemin central, tandis
disposent généralement d'un logiciel intégré pour évaluer
que le bord inférieur est obtenu en décalant le bord supérieur.
le rayon de braquage et signaler les parcours qui violent cette contrainte,
bord à une distance Nous. Il convient de noter que la quantité maximale de fibres
augmentant automatiquement le rayon lorsqu'il est plus petit que le
L'angle est obtenu au trait le plus élevé du parcours. Obtenir
rayon de braquage minimum [15]. Cependant, cela entraînera un écart indésirable s'il n'est
l'angle de fibre maximum de umax, le parcours de référence maximum
pas pris en compte lors de la définition des lignes centrales de la trajectoire d'origine.
L'angle au bord du panneau est :

d
Le rayon de braquage peut éventuellement être réduit en utilisant des
Tmax ¼ arc sinus 1 þ sinþT0Þ ð7Þ
les largeurs de remorquage, qui devraient être moins sujettes au froissement que les plus larges q1
ceux. Le rayon q, est défini par :
où q1 est le rayon associé à Tmax :
1
d
q¼ ð2Þ 2j þ Cw
ðsin T1 sin T0Þ q1 ¼ ð8Þ
2ð1 sinðT0ÞÞ
Le centre de courbure Pc de la trajectoire de la tête de référence (voir Fig. 1b) est
situé aux coordonnées :
3. Méthodologie pour modéliser les stratifiés dirigés par remorquage
xc ¼ q sin T0 yc ¼ ð3Þ
La méthodologie pour construire un modèle FE tridimensionnel pour un
q cos T0 ð4Þ
plaque plate dirigée par remorquage, comporte différentes étapes. Les principaux enjeux
Le long du chemin de référence, les expressions suivantes s'appliquent pour le il faut considérer comment calculer les angles d'orientation locaux des fibres,
coordonnée y et l'angle d'orientation u : la localisation des zones sans fibres et le type de modélisation
Machine Translated by Google

O. Falcon et coll. / Composites : Partie B 56 (2014) 660–669 663

technique employée. Les deux premiers points sont abordés dans les précédents P0 ¼ ½xP;ðyP n WeÞ ð9Þ

fonctionne [5,17,10,18,11]. Dans ce travail, nous les avons adaptés à notre

où:
objectif de construire un maillage structuré. Notre champ d'action s'est limité à
une configuration pour la méthode de remorquage avec paramètre de couverture yP yQ
n ¼ cellule ð10Þ
de 0% [15]. La méthodologie permet une variation linéaire de l'angle et un Nous
échelonnement. L'utilisation d'un paramètre de couverture différent de 0%
nécessitent une stratégie de modélisation différente.
La valeur de n représente le numéro de cours auquel appartient le point matériel
P, et elle est obtenue avec l'équation. 10. La fonction « plafond » arrondira au nombre
3.1. Orientation des matériaux entier le plus proche vers le positif.
infini. Dans cette expression, yQ correspond à la coordonnée y de

La procédure pour déterminer l'angle d'orientation local des fibres point Q, sur l'arête du parcours 0. L'angle d'orientation de la fibre u

prend en compte le chemin de référence défini comme un chemin à valeur unique a été calculé pour chaque élément fini dans un domaine de maillage régulier, comme

fonctionner dans un domaine à deux dimensions. Ici, l'orientation le montre la Fig. 3.

angle pour un point arbitraire P de coordonnées (x, y) situées à L'angle d'orientation des fibres peut être calculé analytiquement avec

le centre de gravité de l'élément est calculé. L'orientation des fibres en P est Éq. 5 évaluant le point Pf ou, géométriquement, avec la fonction 'angle (Pc, P0)'. Cette
fonction renvoie le sens anti­horaire
la même que celle du point P0 situé au cap de référence, ou
cours 0 (voir Fig. 3). Les coordonnées de ce point sont calculées angle entre l'axe des x et un segment défini par deux points quelconques,

avec l'expression suivante : et constitue une autre façon d'obtenir l'angle u sans avoir à
calculer le point Pf. Cette méthode ne peut être utilisée qu'avec la constante
chemin de courbure car le centre de courbure du chemin de référence Pc
est toujours connu. En parcourant le maillage complet, élément par
élément, l'angle d'orientation des fibres du modèle à direction remorquée peut
être obtenu pour chaque centre d’élément.
Dans le cas de la configuration à courbure constante (CC), l'angle d'orientation
pour les points P, P0 et Pf est le même (voir Fig. 3) et
est calculé avec l'expression suivante :

uðPÞ ¼ uðP0Þ ¼ uðPfÞ ¼ angleðPc; P0Þ90 ð11Þ

Dans notre étude, le point Pf représente un pied­point de P0 sur la courbe de

référence (voir Fig. 3). Dans le chemin CC, ce point est facilement calculé,
tandis que dans le chemin LV, il est calculé à l'aide d'une méthode itérative non linéaire
solveur. Dans une configuration avec variation d'angle linéaire (LV), le
le centre de courbure Pc change continûment, et le point Pf, associé au point P0, doit
d'abord être localisé afin de calculer le
angle d'orientation (équation 1). Le problème se réduit au calcul d'un point pied situé
à la distance la plus courte de P0 sur un
courbe de référence f(x), qui peut être calculée en intégrant dy/
dx = bronzage (u(x)).
Simplifier l'équation. 1, pour le cas particulier de / = 0 et T0 – T1,
la fonction f(x) se lit comme suit :

ð12Þ
fðxÞ ¼ Z tanduðxÞÞdx ()
d ðT1 T0Þ
() fðxÞ ¼ lnðj cosðT0ÞjÞ ln cos x þ T0 ð13Þ
ðT1 T0Þ d

Une méthode de discrétisation a été utilisée pour approximer la

courbe f(x) avec une polyligne. La polyligne est représentée par un ensemble
S = {S1, S2, ..., Sn} de n segments de droite, tels que l'approximation
l'erreur E (f(x); S) est minimisée, selon la valeur définie dans le paramètre de dis­
crétisation. Le problème est simplifié pour trouver la distance orthogonale minimale,
du point P0 à un segment Si. Le
la pente kSi de chaque segment est connue, et il est très simple de trouver
la droite orthogonale de pente ð1=kSi Þ, passant par le point P0,
afin de déterminer le point équivalent à Pf (P0 f 2 Si : ,Pf). Ce
La méthode est simple, mais suffisamment rapide, précise et fiable pour résoudre
ce problème.

3.2. Étapes pour créer le modèle VSP

La méthodologie de prétraitement proposée est combinée avec

le progiciel d'analyse par éléments finis disponible dans le commerce, Simulia Abaqus.
Cette méthodologie est utile pour les techniques de modélisation de maillage structuré
et régulier. Pour le spécifique
cas avec un maillage régulier, les étapes 3 et 4 de la méthodologie sont
Fig. 3. Spectre de l'angle d'orientation u dans une configuration à constante simplifié et le maillage et les modèles d'assemblage sont obtenus
courbure avec drapage ±h0j45i et coupe d'un côté. indépendamment du processus de remorquage. L'approche de modélisation
Machine Translated by Google

664 O. Falcon et coll. / Composites : Partie B 56 (2014) 660–669

taille similaire, distinguant clairement les zones renforcées par des fibres des
zones en résine.
– Les assises sont modélisées à l'aide d'un maillage tridimensionnel avec des
éléments en brique, tandis que les zones triangulaires sans fibres peuvent
être modélisées correctement à l'aide d'éléments en brique et d'éléments
triangulaires­prismatiques juste à l'intersection des câbles courbes avec
l'assise de bord. En utilisant cette combinaison d’éléments, le maillage est
obtenu avec des éléments de bonne qualité de forme.
La localisation précise de ces zones au moyen d'une grille structurée
d'éléments de maillage contrairement à une technique de maillage régulière
est illustrée à la Fig. 6. Si une répartition précise des contraintes interlaminaires
est requise, il est également possible d'affiner le maillage dans le direction à
travers l'épaisseur, figure 5d.

Étape 4 : Calcul de l’angle d’orientation des fibres.

– Une fois les maillages créés, l'orientation locale du matériau est attribuée à
chaque élément, selon la procédure décrite dans la section 3.1.

Étape 5 : attribution des propriétés du matériau.

– Les propriétés d'élasticité et de résistance du polymère renforcé de fibres de
carbone (CFRP) sont introduites et attribuées à chaque région particulière du
Fig. 4. Méthodologie de modélisation d'un VSP avec une approche de maillage structuré.
modèle. Ces régions sont associées aux plis extérieurs et encastrés.
Auparavant, la localisation des zones sans fibres dans le système de matériaux
se compose d'une séquence d'étapes, comme le montre la figure 4 et expliquée était définie selon la procédure développée par Blom et al. [dix]. De plus, les
comme suit : propriétés globales de la résine sont attribuées aux éléments situés dans ces
zones.
Étape 1 : Définition du panneau.
– La méthodologie commence par la définition du panneau à rigidité variable.
Les données d'entrée sont liées aux spécifications géométriques (T0, T1, d), Étape 6 : Processus d’assemblage pour construire le stratifié.
aux contraintes de fabrication et au drapage du stratifié, selon la nomenclature – Les plis sont assemblés dans un stratifié en tenant compte de sa séquence
[/ ± hT0jT1i . Les paramètres de conception incluent l'épaisseur nominale du de jalonnement. A noter que des configurations hybrides (stratifiés avec des
pli, les dimensions du panneau et la configuration VSP adoptée (CC ou LV). plis de fibres à la fois droits et curvilignes) peuvent également être envisagées
Les paramètres de fabrication (qmin, tw, nt) se réfèrent au rayon de braquage (par exemple, [±45/± < 0j80i4]s).
minimum, à la largeur des câbles et au nombre de câbles placés par la – La nature propre des maillages structurés provoque des maillages non
machine AFP, selon la méthode de fabrication adoptée (Tow­drop ou Overlap). conformes dus à la non­coïncidence des nœuds (voir Fig. 6a). Les interactions
de propriétés entre toutes les surfaces liées, qui incluent les interfaces des plis
ainsi que l'interface entre les couches adjacentes, sont définies. Une
formulation de contrainte cinématique est définie entre les interfaces de plis
Étape 2 : Création de la géométrie du parcours. ainsi qu'entre toutes les limites de couche. Pour les analyses de propagation
– La procédure de génération du modèle EF de chaque pli est similaire au des dommages, le comportement interlaminaire entre les différents plis peut
processus de fabrication, où le pli est formé par un ensemble successif de être modélisé au moyen d'un cadre de contact cohésif basé sur la surface.
couches décalées. Chaque cours représente une partie­instance du pli et
également un sous­domaine indépendant dans le modèle global. La géométrie
finale du pli est construite par l'assemblage des éléments générés.
Étape 7 : Conditions aux limites.
cours. – La dernière étape de la méthodologie proposée est la spécification du type
de conditions aux limites et de leur localisation sur le domaine de simulation.
Étape 3 : Génération de maillage. L'application des conditions aux limites définies dans ce travail sera expliquée
– Chaque « partie de cours » est partitionnée dans l'esquisse structurée, selon plus en détail dans la section 4.
les paramètres de conception du VSP. Dans cette opération, la forme, la taille
moyenne et la position de chaque élément fini seront contrôlées. De plus, le
processus d'ensemencement est entièrement contraint en fonction des faces Le tableau 1 résume les différentes configurations et stratégies de modélisation
partitionnées du modèle. qui peuvent être simulées avec l'ensemble d'outils actuel. Des trajets de référence à
Les étapes du processus de génération de maillage sont illustrées à la figure courbure constante ou des trajets avec variation linéaire de l'orientation des fibres sont
5. Premièrement, le câble est subdivisé longitudinalement en considérant une autorisés. Dans les deux cas, la modélisation peut prendre en compte que les
taille d'élément pratique (dans cet exemple, chaque câble a été subdivisé en panneaux sont fabriqués en coupant les câbles de fibres sur un ou deux côtés de la
deux éléments). Ensuite, les points d'intersection Pi entre les subdivisions et trajectoire de la tête de machine.
le bord du parcours sont calculés, comme le montre la figure 5a. Deuxièmement, Dans ce tableau, le signe astérisque indique que le calcul des zones sans fibres a été
les parcours sont divisés transversalement à partir des points d'intersection. mis en œuvre principalement lorsque les câbles de chute sont coupés d'un seul côté.
Chaque secteur résultant s'étend sur un angle hi, comme le montre la figure La stratégie de coupe du remorquage, qui consiste à larguer des remorquages des
5b. Troisièmement, en fonction de la taille moyenne des éléments, l'angle hi deux côtés d'un parcours, est en cours de développement. De même, différentes
est subdivisé pour compléter la discrétisation des parcours, Fig. 5c. Le résultat stratégies de modélisation peuvent être utilisées en fonction des exigences de
de ce processus est que toute la couche est divisée en éléments triangulaires simulation. Pour les analyses de rupture où les effets locaux sont importants, un
et quadrangulaires de domaine 3D détaillé peut être nécessaire, tandis que pour l'évaluation de la rigidité et
du flambement, des analyses 2D idéalisées sont recommandées.
Machine Translated by Google

O. Falcon et coll. / Composites : Partie B 56 (2014) 660–669 665

(un)

(b)

(c)

(d)
Fig. 5. Génération de maillage à partir d'une esquisse structurée en « remorquage ».

Fig. 6. Génération de maillage d'un VSP. (a) Maillage non conforme. (b) Modélisation précise des zones sans fibres avec une approche de maillage structuré. (c) Zones sans fibres dans un
domaine de maillage régulier.

recommandé. Bien que les outils de modélisation permettent de nombreuses options ; variation de l'angle de la fibre tive de T0 à T1 du trajet de la tête de référence pour
A titre d'exemple, dans cet article nous considérons pour le modèle à maillage un rayon de braquage de q = 212,13 mm. Pour cette analyse particulière
structuré une configuration à courbure constante du trajet de référence et les câbles T0 = 0 et T1 = 45 et la séquence d'empilement de stratifié symétrique et équilibrée
ont été coupés uniquement sur un côté du parcours. pour le panneau était de [± < 0j45i/± < 45j0i]s. La frontière
Les conditions de l'analyse sont représentées sur la figure 7. Ces conditions ont
4. Géométrie du modèle et discrétisation par éléments finis été choisi compte tenu d'un essai de traction réalisé avec une force contrôlée
déplacement du bord supérieur du panneau jusqu'à 3% du panneau
Des tests virtuels ont été réalisés sur un panneau de 150 à 150 mm sous charge hauteur, tandis que le bord inférieur était fixe. L'analyse FE a été
de traction. Les dimensions du panneau permettent une représentation réalisé avec intégration temporelle implicite (Abaqus/Standard).
Machine Translated by Google

666 O. Falcon et coll. / Composites : Partie B 56 (2014) 660–669

Tableau 1 Tableau 2
Stratégies de conception et de modélisation pour VSP. * indique que le calcul des zones sans fibres a été mis Propriétés du système CFRP Hexply AS4/9773 et de la résine époxy [3].
en œuvre principalement lorsque les câbles de chute sont coupés d'un seul côté.
Propriétés élastiques du CFRP
Conception/fabrication E11 = 129,8 GPa ; E22 = 9,2 GPa ; E33 = 9,2 GPa ;
Configuration – Variation linéaire (BT) m12 = 0,36 ; m13 = 0,36 ; m23 = 0,487 ;
– Courbure constante (CC) G12 = 5,1 GPa ; G13 = 5,1 GPa ; G23 = 3,1 GPa
Coupure de cours – Remorquage tombant d’un côté
Points forts des plis CFRP
– Remorquage tombant des deux côtés
XT = 2070MPa ; XC = 1 160 MPa ; YT = 29,0 MPa ;
YC = 157,9 MPa ; SL = 91 MPa
Modélisation/simulation
Une résine époxy
Niveau de détail – Répartition idéale des fibres
E22 = 9,2 GPa ; m12 = 0,36
– Distribution discrète des cours
– Zones « sans fibre » incluses⁄
Discrétisation du maillage – Maille régulière
– Maille structurée
Domaine de simulation ­2D
le gradient de contrainte à travers l'épaisseur ne peut pas être décrit avec précision
– 3D
avec ce modèle, mais pour l'application d'un mésomodèle pour
simulation de dommages ou de défaillances des plis, c'est­à­dire un modèle qui prédit le
apparition de fissures traversantes, état de contrainte constant
on suppose qu'il traverse l'épaisseur des plis. Cependant, ce maillage
est suffisant pour obtenir une discrétisation adéquate des cours et
zones sans fibres dans les mailles structurées et régulières, et
pour permettre des comparaisons entre eux. La taille de l'élément a été conservée
constante pour tous les modèles afin d’éviter les effets de taille de maillage sur le
résultats. Aussi, pour exclure l’influence des conditions aux limites
ou effets de bord libre, la surface des panneaux sélectionnée pour la contrainte
analyses a une distance de décalage de 16 mm par rapport au bord libre.
L'un des avantages d'un modèle tridimensionnel est qu'il
inclut l'effet des contraintes hors plan, permettant ainsi leur
inclusion dans les critères d’échec. La prédiction de défaillance entre le
modèles ont été calculés au moyen du critère de défaillance tridimensionnel pour
les stratifiés renforcés de fibres proposé par Catalanotti
et coll. [20]. Ce critère prend en compte toutes les composantes du
tenseur de contraintes et « l'effet in situ », qui se caractérise par
résistances in situ à la traction transversale plus élevées YT et résistances au cisaillement
est

(SxTest; SL ) dans le pli lorsqu'il est contraint par des plis de

est

orientations des fibres. La défaillance des zones sans fibres n'a pas été
considéré.

Fig. 7. Géométrie des panneaux à rigidité variable et conditions aux limites. 5. Résultats et discussion

Dans les configurations considérées dans ce travail, un minimum de Déterminer l’influence des différentes configurations et maillages
quatre plis sont nécessaires pour l’équilibrage global du stratifié. Ce approches de modélisation sur la réponse mécanique et la rupture
c'est parce que l'équilibrage nécessite une variation symétrique de la rigidité comportement des panneaux remorqués, diverses simulations ont été réalisées.
sur les axes x et y. Dans ce travail, tous les stratifiés analysés contiennent 8 plis, Le premier cas correspond à un dimensionnement à courbure constante
avec une épaisseur nominale de pli t = 0,2 mm, modélisé avec une approche IDEAL où l'orientation des fibres à
et sont globalement équilibrés et symétriques. Quant à la fabrication chaque point suit la courbe de référence. Le deuxième cas correspond à un plan
spécifications, la largeur du remorquage, le paramètre tw prend une valeur de de variation de fibre linéaire (LV), tel que décrit dans
3,175 mm et chaque parcours est formé de 16 traits. Ces paramètres de fabrication Article 2.1. Dans le cas particulier du LV, les deux côtés du parcours
ont été supposés simplifier les tests virtuels et sont coupés et la zone sans fibres n'est pas prise en compte. Ces
correspondent à un environnement de fabrication possible. Les propriétés du deux premières simulations ont été modélisées avec un maillage régulier. L'autre
système de matériaux CFRP et les propriétés linéaires du les modèles analysés utilisent les configurations à courbure constante (CC)
Les zones de résine sont présentées dans le tableau 2. Ces propriétés matérielles ont été avec distribution de cours discrète, décrite dans la section 2.2 et avoir
attribué uniquement aux régions composites CFRP, alors que les zones sans été modélisé dans un domaine de maillage régulier (CC­Regular) et structuré (CC­
fibres étaient considérées comme constituées d'une résine isotrope Struc­tured) avec et sans zones sans fibres incluses. Dans
matériel afin de simplifier le modèle. Des études expérimentales récentes ont Dans cet article, les modèles avec configuration (CC) présentent des discontinuités
démontré que ces zones comprennent également certaines fibres angulaires plus importantes entre les assises adjacentes que les modèles avec configuration (CC).
des couches adjacentes [19]. (LV), où les deux côtés ont été coupés.
Les modèles éléments finis tridimensionnels à maillage régulier Les parcours de référence sont définis par les angles T0 = 0 et
se composent de 70 688 éléments de brique 3D continus à intégration réduite T1 = 45. Dans un pli idéal, les angles obtenus sont évidemment limités
(C3D8R). Un nombre similaire (71 040 éléments) a été entre T0 et T1. Cependant, les angles locaux des fibres sont également affectés.
utilisé dans le maillage structuré. Dans le modèle de maillage structuré, le par le mode de fabrication et peut atteindre des valeurs supérieures à
les parcours ont été modélisés avec des éléments hexaédriques C3D8R avec les idéaux (voir Fig. 8a­b). La valeur maximale de l'angle d'orientation est davantage
intégration réduite alors que, près des bords libres, des éléments en coin influencée par la coupe du parcours que par
de type C3D6 ont été sélectionnés pour compléter le maillage. La taille des la configuration adoptée (BT ou CC). Les valeurs sont plus élevées dans
éléments hexaédriques était de 1,6 à 1,6 mm en moyenne, et chaque pli stratifié a le modèle n'a été coupé que d'un côté car la distance moyenne jusqu'à
été modélisé avec un seul élément solide. Donc, la courbe de référence est plus grande. Le parcours fabricable
Machine Translated by Google

O. Falcon et coll. / Composites : Partie B 56 (2014) 660–669 667

(un) (b)

(c) (d)
Fig. 8. (a, b) Angles d'orientation des fibres le long de AA0 et BB0 respectivement, (c) indice de rupture des fibres dans un modèle structuré CC à un niveau de charge de traction de 81,2 kN, (d) contrainte (r11) tracée
le long du chemin BB0 .

Les définitions produisent des discontinuités notables dans les angles d’orientation dans 81,2 kN. La figure 8c montre la distribution de l'indice de rupture en traction des fibres.
les limites des parcours. Les variations d'angle sont en hausse du modèle CC­Structured.
à 15% et 19% le long du trajet vertical (AA0 ) respectivement dans les modèles LV et CC, Les différences dans l'indice de rupture en traction des fibres entre les
alors que dans le modèle IDEAL l'angle est constant Les modèles sont principalement liés à la répartition de l'angle d'orientation des fibres sur
(voir Fig. 8a). le parcours. Les principales différences dans r11 tracées le long
La rigidité globale calculée présente de légères différences entre les différents modèles. BB0 , montrent une dépendance aux angles d'orientation locaux dans chaque
La valeur obtenue pour le modèle IDEAL est de 21 kN/mm. Cette valeur augmente d'environ modèle. Cette influence est évidente lorsque la figure 8b est comparée à
4% pour le LV Figure 8d. Les valeurs d'angle le long du trajet BB0 ont une évolution qui
configuration et entre 7% et 8% pour les modèles CC sans est similaire aux valeurs de contrainte des fibres (r11) (tracées sur la figure 8d). Plus haut
zones sans fibres. Lorsque les zones sans fibre sont incluses, la les valeurs des angles d'orientation des fibres impliquent des valeurs plus élevées de la fibre
la rigidité du panneau diminue de 3 % par rapport à un modèle similaire sans force de tension. Par conséquent, les résultats montrent que les modèles qui prennent
zones sans fibres. Il n’y a pas de différences significatives dans le monde prendre en compte les effets du processus de fabrication, prédire la fibre
rigidité du panneau en comparant le maillage régulier avec le maillage structuré. échec à des valeurs inférieures aux modèles idéaux. Les auteurs supposent que
L'indice de rupture en traction des fibres et l'indice de rupture en traction de la matrice le modèle avec l'approche de la couche idéale (IDEAL) fournit une surprédiction de la
pour différents modèles ont été analysés. Rupture de traction de la matrice valeur de charge critique pour la rupture de la fibre. L'indice de rigidité et de rupture du
est le premier mode de défaillance observé et il devrait se produire à premier pli se rapporte à la valeur de la charge critique
niveaux de charge similaires dans tous les modèles simulés. L'échec de la matrice de 81,2 kN sont repris dans le tableau 3.
L'indice atteint des valeurs critiques à des niveaux de charge proches de 25 kN. Sur le La fissuration matricielle a été analysée pour les modèles à constante
d'autre part, la rupture en traction des fibres apparaît à des charges plus élevées en fonction configuration de courbure (CC) avec des zones riches en résine, utilisant à la fois la
sur le modèle analysé. Par exemple, la valeur de la traction des fibres les mailles régulières et structurées. Ils ont été comparés à
Il a été fait référence à l'indice de défaillance par rapport à l'indice de défaillance du premier l'approche IDEAL et une configuration sans zones sans fibre.
pli du modèle avec CC­Structured à une valeur de charge de Comme prévu, les couches les plus critiques sont les couches externes. La matrice
Machine Translated by Google

668 O. Falcon et coll. / Composites : Partie B 56 (2014) 660–669

Tableau 3
L'indice de rigidité et de rupture du premier pli se réfère à la valeur de charge critique de 81,2 kN.

Approche de configuration et de modélisation Rigidité globale (kN/mm) Indice de traction des fibres normalisé (%)

IDÉAL 20,87 74,8

BT 21,73 75,9
CC­Régulier 21,72 95,7
CC­Structuré 21,85 100,0

LaRC_, FMT
(Moyenne : 75 %)

1.000
0,917
0,833
0,750
0,667
0,583
0,500
0,417
0,333
0,250
0,167
0,083
0,000

(un) (b) (c) (d)

Fig. 9. Indice de rupture en traction de la matrice au niveau de la couche externe pour un VSP idéal et réalisable à une charge de 25,2 kN. (a) Approche IDÉALE, (b) configuration avec parcours de remorquage discrétisés
dans un domaine à maillage régulier, (c) configuration avec parcours discrétisés et zones sans fibres incluses et (d) configuration avec maillage structuré. Uniquement le composite CFRP
Les sections ont été analysées.

indice de rupture en traction pour les couches externes dans toutes les approches sous le prend en compte les problématiques de fabrication des VSP, une analyse plus précise
la même charge de 25,2 kN est représentée en 9. Dans l'approche IDEAL (Fig. 9a), la représentation numérique a été obtenue par rapport à
l'indice de défaillance de la matrice est constant dans la direction verticale. Dans l’approche de la couche idéale. Les étapes de prétraitement des deux
le modèle sans zones sans fibres (Fig. 9b), l'indice de rupture présente une localisation et le maillage structuré sont générés automatiquement. Si la taille moyenne des
évidente dans les zones avec des valeurs d'angle élevées ou éléments est similaire, le nombre d'éléments dans le maillage est
avec des discontinuités d'orientation des matériaux (situées entre les de la même ampleur et le coût de calcul est pratiquement
cours adjacents). Dans les modèles CC­régulier et CC­structuré équivalent. Bien que les modèles simulés ne montrent aucune différence significative
(Fig. 9c et 9d), l'indice de rupture en traction de la matrice présente une distribution en ce qui concerne la rigidité globale, ils se comportent différemment lorsque les
avec un motif très similaire, qui est plus complexe en raison de la critères de rupture du premier pli sont pris en compte. La fibre
influence des zones de résine. l'indice de rupture en traction pour la rupture en traction longitudinale, montre une
Comme prévu, dans le VSP, la rigidité non uniforme dans le plan dépendance significative sur les angles d'orientation obtenus selon
La distribution entraîne de grands gradients de contraintes qui agissent comme des au processus de conception du VSP. Le mode de défaillance initial est lié à
discontinuités et des zones de concentration de contraintes. De tels gradients fissuration de traction de la matrice et est fortement influencée en raison de la
pourraient contribuer à l'amplification des contraintes interlaminaires et à la discontinuité entre les couches et également amplifiée par les chutes de câble.
possibilité de délaminage dans ces zones. En ce qui concerne l'échec Il convient de souligner qu’aucune différence significative n’a été constatée
critères, aucune différence significative n’est appréciée entre les se trouve entre les maillages réguliers et structurés. Cependant,
CC­regular et le modèle CC­structuré. Tout le potentiel du le maillage structuré devrait être plus adapté à l'étude
une approche de maillage structuré sera utile dans la simulation du la progression des dommages interlaminaires dans les zones de contact entre
comportement de dommage progressif non linéaire dans VSP. couches adjacentes. La méthodologie proposée pour générer des modèles FE pour
VSP avec un maillage structuré et prendre en compte
compte des contraintes de fabrication, est un outil puissant pour les études futures
6. Conclusions lorsqu'il s'agit d'aborder des simulations plus avancées telles que les dommages
intralaminaires progressifs et/ou les dommages par délaminage
Dans le présent travail, une méthodologie pour modéliser les véhicules à direction remorquée prédictions.
des panneaux dans un domaine tridimensionnel ont été proposés. UN
générateur de maillage structuré pour la modélisation paramétrique de Remerciements
panneaux de fibres dirigés, et qui peuvent être utilisés pour la conception paramétrique
et l'analyse FE, a été développée. Différentes configurations et Les auteurs reconnaissent le soutien financier du ministère de l'Économie et de
les approches de modélisation ont été analysées, et leur influence dans la Compétitivité de l'Espagne dans le cadre du projet
la prédiction de certains cas de rupture de panneaux sous chargement de traction MAT2009­07918. Le premier auteur remercie la Generalitat
ont été évalués. En développant un outil informatique qui de Catalogne pour la subvention FI (AGAUR, IUE/2365/2009).
Machine Translated by Google

O. Falcon et coll. / Composites : Partie B 56 (2014) 660–669 669

Les références [10] Blom A, Lopes CS, Kromwijk PJ, Gürdal Z, Camanho PP. Un modèle théorique pour étudier l'influence
des zones de chute de câble sur la rigidité et la résistance des stratifiés à rigidité variable. Compos
Mater 2009 ; 43 : 403­25.
[1] Sebaey TA, González EV, Lopes CS, Blanco N, Costa J. Résistance aux dommages et tolérance aux
[11] Blom AW, Stickler BF, Gürdal Z. Conception et fabrication d'une coque cylindrique à rigidité variable.
dommages des stratifiés CFRP dispersés : conception et optimisation.
Dans : Actes de la 30e conférence internationale SAMPE Europe ; 2009.
Compos Struct 2013 ; 95 : 69­76.
[2] Gürdal Z, Tatting BF, Wu KC. Panneaux à rigidité variable : effets de la variation de rigidité sur les
[12] Fagiano C, Abdalla M, Kassapoglou C, Gürdal Z. Récupération de contrainte interlaminaire pour les
réponses dans le plan et au flambement. Compos Partie A : Appl Sci Manuf 2008 ; 39 : 911–22.
éléments finis tridimensionnels. Compos Sci Technol 2010 ; 70 : 530–8.
[13] Díaz J, Fagiano C, Abdalla M, Gürdal Z, Hernández S. Une étude des contraintes interlaminaires dans
[3] Lopes CS, Gürdal Z, Camanho PP. Adaptation pour la résistance des panneaux composites en fibres
les plaques à rigidité variable. Compos Struct 2012 ; 94 : 1192–9.
directionnelles avec découpes. Compos Partie A : Appl Sci Manuf 2010 ; 41 : 1760–7.
[14] Documentation en ligne ABAQUS version 6.11. Manuel d'utilisation de l'analyse 35.1.10.
[4] Duvaut G, Terrel G, L.F, Verijenko VE. Optimisation des fibres renforcées
SIMULIA Inc. ; 2010.
matériaux composites. Compos Struct 2000 ; 48 : 83–9.
[15] Gürdal Z, Tatting BF, Wu KC. Technologie de placement de remorquage et problèmes de fabrication
[5] Blom A, List J, Stickler P, Gürdal Z. Optimisation d'un cylindre composite en flexion en adaptant les
pour les structures composites stratifiées. Dans : 46e conférence AIAA/ASME/ASCE/AHS/ASC sur
propriétés de rigidité dans la direction circonférentielle. Compos Partie B : Eng 2010 ; 41 : 157­65.
les structures, la dynamique structurelle et les matériaux ; 2005.
[16] Campbell F. Technologie de fabrication des matériaux de structure aérospatiale. Elsevier SA; 2006.
[6] Gürdal Z, Olmedo R. Stratifiés composites avec des orientations de fibres spatialement variables :
concept de panneau à rigidité variable. Dans : Actes de la 33e conférence AIAA/ASME/ASCE/AHS/
[17] Waldhart CJ, Ribbens C, Gürdal Z. Analyse de stratifiés à rigidité variable, à fibres parallèles et à câble
ASC sur les structures, la dynamique des structures et les matériaux ; 1992. p. 798­808.
placé. Dans : Actes de la 37e conférence AIAA/ASME/ASCE/AHS/ASC sur les structures, la
dynamique structurelle et les matériaux (SDM). Salt Lake City, UT ; 1996.
[7] Lopes CS, Camanho P, Gürdal Z, Tatting BF. Analyse de rupture progressive de panneaux composites
remorqués à rigidité variable. Structure des solides 2007 ; 44 : 8493–516.
[18] Blom A, Tatting BF, Hol JM, Gürdal Z. Définitions du chemin des fibres pour les coques coniques
[8] Wu KC, Gürdal Z, Starnes JH. Réponse structurelle des panneaux à rigidité variable, chargés en
adaptées élastiquement. Compos Partie B : Eng 2009 ; 40 : 77­84.
compression et placés en remorquage. Dans : Actes de la 43e conférence AIAA/ASME/ASCE/AHS/
[19] Croft K, Lessard L, Pasini D, Hojjati M, Chen J, Yousefpour A. Étude expérimentale de l'effet du
ASC sur les structures, la dynamique des structures et les matériaux. Denver, Colorado ; 2002. AIAA
placement automatisé des fibres induit des défauts sur les performances des stratifiés composites.
2002­1512.
Compos Partie A : Appl Sci Manuf 2011 ; 42 : 484–91.
[9] Lopes CS, Gürdal Z, Camanho P. Panneaux composites à rigidité variable : améliorations du flambage
[20] Catalanotti G, Camanho P, Marques AT. Critères de rupture tridimensionnelle pour les stratifiés renforcés
et de la rupture du premier pli par rapport aux stratifiés à fibres droites. Structure informatique 2008 ;
de fibres. Compos Struct 2013 ; 95 : 63­79.
86 : 897­907.