SEDIRA 2013 Archivage
SEDIRA 2013 Archivage
SEDIRA 2013 Archivage
THÈSE EN CO-TUTELLE
Présentée à l'université de Biskra
Pour Obtenir le diplôme de
Doctorat en Sciences
Université Mohamed Khider - Biskra
par :
Lakhdar SEDIRA
i
Remerciements
ii
permanents, techniciens et secrétaires, pour leurs collaboration et l'ambiance
amicale dans les groupes.
L'intégration des éléments développés dans le code ABAQUS n'aurait pas été
possible sans l'implication de mes collègues: Messieurs Wajdi ZOUARI et
Mustapha ASSARAR que je souhaite également remercier pour tout le travail
réalisé et l'accueil convivial à IUT de Troyes.
Je voudrais remercier toute l'équipe de l'école doctorale: directeurs,
responsables et secrétaires, d'avoir veillé à la bonne formation dans notre
parcours doctoral, pour leur préséance et leur soutien administratif. Je suis
très reconnaissant à Madame Estelle ODINOT de l'effort considérable dédié à
la gestion de notre parcours de la thèse.
Merci également aux personnels du CROUS de Reims, qu'ils étaient très
accueillants pendant mon séjour de dix-huit mois. Un remerciement particulier
à la responsable du service des relations internationales, Madame Sylvie
PERIQUET, qui a veillé à nous assurer les meilleures conditions de vie, sa
préséance a été remarquable dans notre intégration culturelle et sociale.
iii
Table des matières
iv
2.2.5 Matrice de rigidité élémentaire ....................................................................58
2.3 CALCUL DES CONTRAINTES PLANES ET DE CISAILLEMENT TRANSVERSAL ............59
2.4 FORMULATION DE LA MATRICE MASSE ELEMENTAIRE DU MODELE PLAQUE
COMPOSITE ..................................................................................................................61
2.5 CONCLUSION ........................................................................................................62
CHAPITRE 3 VALIDATION NUMERIQUE DU MODELE DE PLAQUES
COMPOSITES DMQPZ ................................................................ 63
3.1 INTRODUCTION ....................................................................................................63
3.2 RESULTATS DES CAS TESTS STANDARDS DE PLAQUES MULTICOUCHES
STRATIFIEES/SANDWICH ...............................................................................................64
3.2.1 Plaques stratifiées composites à 3 et 9 couches ............................................64
3.2.2 Plaque carrée sandwich sous chargement uniforme ....................................70
3.2.3 Plaque carrée et rectangulaire avec deux types de charges ..........................73
3.2.4 Plaque carrée sandwich (f/c/f), charge uniforme........................................80
3.2.5 Plaque carrée sandwich (f/c/f) sous chargement sinusoïdal .....................84
3.2.6 Plaque carrée sandwich (f(0,90),c, f(90,0)) sous charge uniforme ...............86
3.3 TESTS DE VIBRATIONS LIBRES DES STRUCTURES COMPOSITES .............................88
3.3.1 Plaque carrée isotrope simplement supportée ..............................................88
3.3.2 Plaque carrée sandwich simplement supportée ...........................................90
3.3.3 Carton ondulé ...............................................................................................91
3.3.4 Plaque carrée multicouche croisée (0/90/90/0) ..........................................96
v
5.2 APPLICATIONS AUX STRUCTURES COQUES ISOTROPES .......................................126
5.2.1 Cylindre pincé avec diaphragme ................................................................126
5.2.2 Poutre vrillé sous chargement concentré ....................................................130
5.2.3 Hémisphère pincé sous chargement diamétralement opposées ..................132
5.2.4 Paraboloïde hyperbolique ...........................................................................133
5.3 APPLICATIONS AUX STRUCTURES COQUES COMPOSITES (STRATIFIES
MULTICOUCHES, SANDWICHS) ....................................................................................134
5.3.1 Cylindre orthotrope sous pression ..............................................................134
5.3.2 Plaque sandwich non-symétrique ..............................................................136
5.3.3 Cylindre multicouche croisé .......................................................................138
5.4 ETUDE DE VIBRATIONS LIBRES DE STRUCTURES COMPOSITES............................139
5.4.1 Vibrations libres d’une aube de compresseur isotrope ...............................139
5.4.2 Vibrations libres de panneaux composites .................................................140
CHAPITRE 6 MODELE DISCRET DE COQUE NON-LINEAIRE
GEOMETRIQUE. FORMULATION ET EVALUATION ........ 144
6.1 INTRODUCTION ..................................................................................................144
6.2 FORMULATION VARIATIONNELLE LAGRANGIENNE ACTUALISEE A CHAQUE
ITERATION (FLAI) ......................................................................................................144
6.2.1 La matrice de rigidité linéaire ...................................................................146
6.2.2 La matrice de rigidité géométrique de l'élément DMQS ............................146
6.2.3 Vecteur de forces résiduelles .......................................................................147
6.3 PRESENTATION ET DISCUSSION DES RESULTATS ................................................147
6.3.1 Poutre console sous chargement de flexion.................................................148
6.3.2 Panneau cylindrique articulé .....................................................................151
6.3.3 Hémisphère pincé .......................................................................................156
6.3.4 Panneau encastré .......................................................................................160
CONCLUSION GENERALE ET PERSPECTIVES .............................................. 162
B.1 COMPOSANTES DES VITESSES D'UN POINT QUELCONQUE D'UNE PLAQUE ........181
vi
Notations
vii
HT Effort tranchant dû au zigzag
σ Tenseur des contraintes de Cauchy
S Tenseur des contraintes de Piola-Kirchhoff de seconde espèce
E Tenseur des déformations de Green-Lagrange
F Tenseur gradient de déformations
L Tenseur gradient de déplacements
Wint Travail virtuel interne
Wext Travail virtuel externe
R Vecteur de forces résiduelles (hors-équilibres)
«*» Grandeur virtuelle
ou « ,x » Dérivation partielle par rapport à x
x
dx Différentielle ou incrément
x Variation
x Accroissement
Ni Fonctions d’interpolation
J Jacobien
j Jacobien inverse
k Matrice de rigidité élémentaire
k t
e
Matrice tangente (élémentaire)
k l
e
Matrice de rigidité linéaire
k e
σ Matrice de rigidité géométrique (des contraintes initiales).
M Matrice masse
Valeur propre
Fréquence propre
Masse volumique
E Module d’Young
Coefficient de Poisson
E1 , E2 , E3 Modules d’Young dans les directions d’orthotropie 1, 2 et 3
G12 , G13 , G23 Modules de cisaillement dans les plans1-2, 1-3 et 2-3
12 , 13 , 23 Coefficients de Poisson
hk Epaisseur de la couche k
nC Nombre des couches
viii
Liste des figures
ix
Figure 3.10. Distribution de la contrainte plane suivant l’épaisseur de la plaque
multicouche, ....................................................................................77
Figure 3.11 Distribution de la contrainte de cisaillement suivant l’épaisseur de la
plaque, ...........................................................................................78
Figure 3.12 Distribution de la contrainte de cisaillement suivant l’épaisseur de la
plaque, ...........................................................................................78
Figure 3.13 Plaque sandwich à un seul cœur f/c/f .......................................................81
Figure 3.14 Contrainte du cisaillement transverse au point B, cas 1. .....................83
Figure 3.15 Contrainte du cisaillement transverse au point B, cas 2. .....................83
Figure 3.16 Contrainte du cisaillement transverse au point , maillage
. .............................................................................................................87
Figure 3.17 Plaque carrée isotrope simplement supporté, influence de sur la
fréquence propre ......................................................................................89
Figure 3.18 Montage expérimentale de vibration libre d’une plaque rectangulaire en
carton ondulé ................................................................................................91
Figure 3.19. Accélération du point P. Plaque de longueur L = 200 mm ........................92
Figure 3.20. FFT de l’accélération du point P. Plaque de longueur L = 200 mm ..........93
Figure 3.21 Vibrations libres du carton ondulé ( ), les six premiers modes
propres. .........................................................................................................94
Figure 3.22 Variation de la première pulsation propre pour différents longueurs de la
plaque. ..........................................................................................................95
Figure 3.23 Convergence du modèle DDM pour la plaque multicouche encastrée,
fréquence , S=5 .........................................................................................99
Figure 3.24 Convergence du modèle DDM pour la plaque multicouche encastrée,
fréquence , S=5 .........................................................................................99
Figure 3.25 Convergence du modèle DDM pour la plaque multicouche encastrée,
Fréquence , S=10 ....................................................................................100
Figure 3.26 Convergence du modèle DDM pour la plaque multicouche encastrée,
Fréquence , S=10 ....................................................................................100
Figure 4.1 Discrétisation géométrie d'une coque avec un .... élément quadrilatère courbe
....................................................................................................................102
Figure 4.2 vecteurs de bases .........................................................................................103
Figure 4.3 Présentation des rotations sur les côtés .....................................................107
Figure 5.1 Cylindre pincé avec diaphragmes, données géométriques et mécaniques ..127
Figure 5.2 Cylindre pincé avec diaphragmes: déformée globale d'un huitième du
cylindre (12×12). .........................................................................................127
Figure 5.3 Cylindre pincé. l'allure du déplacement normal le long des côtés DC et
BC. ..............................................................................................................128
Figure 5.4 Cylindre pincé avec diaphragmes, déformée globale, maillage avec l'élément
DMQS (32×100) sur Abaqus. .....................................................................129
Figure 5.5 Poutre vrillée. Données................................................................................130
Figure 5.6 Hémisphère pincé. Données ........................................................................132
x
Figure 5.7 Paraboloïde hyperbolique, données du problème .......................................133
Figure 5.8 Paraboloïde hyperbolique, convergence du déplacement normal .........134
Figure 5.9 Cylindre orthotrope. Géométrie et données. ...............................................135
Figure 5.10 Plaque carrée sandwich non symétrique simplement supportée sous charge
uniformément répartie, données du problème. ...........................................136
Figure 5.11 Distribution de la contrainte suivant l'épaisseur de la plaque sandwich
non symétrique . ........................................................................137
Figure 5.12 Cylindre simplement supporté sous chargement doublement sinusoïdal,
données géométriques. ................................................................................138
Figure 5.13 Modèle simplifié d' aube de compresseur, données géométriques et
mécaniques .................................................................................................140
Figure 5.14 Panneau composite, données géométriques et matérielles. .......................141
Figure 5.15 Fréquences propres de panneaux cylindriques composites. ......................142
Figure 6.1 Poutre console sous chargement de flexion, géométrie et chargement. .......148
Figure 6.2 Panneau cylindrique articulé, géométrie et données. .................................151
Figure 6.3 Convergence du déplacement aux points A et B (h=12.7 mm). ...................152
Figure 6.4 Convergence du déplacement au point A (h=6.35 mm) ..............................154
Figure 6.5 Hémisphère pincé. Géométrie et données ....................................................156
Figure 6.6 Déplacement UA, maillage ...............................................................157
Figure 6.7 Déplacement (-UB), maillage ...........................................................158
Figure 6.8 Déformées de l'hémisphère pincé ................................................................159
Figure 6.9 Panneau cylindrique articulé, géométrie et données. .................................160
Figure 6.10 Panneau encastré, influence du nombre d'éléments. ................................161
Figure A.1 Couche orthotrope .......................................................................................177
Figure A.2 Organigramme de calcul des facteurs de correction de CT par la méthode
totale. ..........................................................................................................179
Figure A.3 Organigramme de calcul des facteurs de correction de CT par la méthode
flexion .........................................................................................................180
xi
Liste des tableaux
Tableau 2.1 Les fonctions de forme utilisées dans l'élément quadrilatère ....................45
Tableau 3.1 Plaque carrée simplement supportée sous chargement doublement
sinusoïdal, 3 couches (0/90/0), contraintes planes maximales ................66
Tableau 3.2. Plaque carrée simplement supportée sous chargement doublement
sinusoïdal, 3 couches (0/90/0), déplacement et contraintes de CT
maximales ..................................................................................................67
Tableau 3.3 Plaque carrée simplement supportée sous chargement doublement
sinusoïdal, 9 couches, contraintes planes maximales, maillage (6x6) ......67
Tableau 3.4 Plaque carrée simplement supportée sous chargement doublement
sinusoïdal, 9 couches, déplacement et contraintes de CT maximales. ......68
Tableau 3.5 Plaque carrée sandwich sous chargement uniforme, déplacement .......71
Tableau 3.6. Plaque carrée sandwich simplement supportée sous chargement
uniforme, contrainte . ...........................................................................72
Tableau 3.7 Flèche et contraintes dans la plaque multicouche (0/90/0) sous
chargement. sinusoïdal, , ...............................................................74
Tableau 3.8 Flèche et les contraintes dans la plaque rectangulaire multicouche
(0/90/0) sous chargement sinusoïdal. ......................................................76
Tableau 3.9 Plaque carrée composite à 3 couches (0/90/0), de même épaisseur
simplement supportée sous chargement uniforme. Déplacement
transversal et des contraintes maximums. ................................................80
Tableau 3.10 plaque carrée sandwich (f/c/f). Résultats de la flèche au centre de la
plaque. ........................................................................................................82
Tableau 3.11 Plaque carrée sandwich (f/c/f). Résultats des contraintes de CT
maximales ..................................................................................................82
Tableau 3.12. Déplacement transverse, contraintes du cisaillement dans la plaque
carrée sandwich (f/c/f) soumise à une charge sinusoïdale .......................85
Tableau 3.13. Contraintes planes dans la plaque carré sandwich (f/c/f), simplement
supportée soumise à une charge sinusoïdale. ............................................85
Tableau 3.14 Plaque carrée sandwich ( f(0,90), c, f(90,0) ) sous chargement uniforme.
déplacement transversal au point C, contrainte plane et contraintes de CT
maximales. .................................................................................................87
Tableau 3.15. Plaque carrée isotrope simplement supportée, ν=0.3, Influence de a/h
sur la fréquence fondamentale (maillage 6 6) ..........................................88
Tableau 3.16. Plaque carrée isotrope simplement supportée : ν=0.3, a/h=10, les quatre
premières fréquences propres .....................................................................89
Tableau 3.17 Fréquence de la vibration libre de la plaque sandwich à 3 couches. .......90
xii
Tableau 3.18 Caractéristiques géométriques et mécanique du carton ondulé utilisé
pour l’essai vibration libre .........................................................................91
Tableau 3.19 Résultats numériques et expérimentaux des pulsations propres. ............95
Tableau 3.20 Vibrations libres du carton ondulé ( ), les six premières
fréquences propres. .....................................................................................95
Tableau 3.21 Fréquences naturelles , de la plaque carrée multicouche (0/90/90/0)
encastrée. ....................................................................................................98
Tableau 5.1 Cylindre pincé avec diaphragmes, convergence de la flèche normale .
..................................................................................................................128
Tableau 5.2 Cylindre pincé avec diaphragmes, ...........................................................129
convergence du déplacement tangentiel ...................................................................129
Tableau 5.3 Poutre vrillée, résultats des déplacements en pour ..................131
Tableau 5.4 Poutre vrillée, résultats des déplacements en pour ..............131
Tableau 5.5 Hémisphère pincé. Convergence de ......................................................133
Tableau 5.6 Cylindre orthotrope sous pression interne. Déplacement maximal et
moments. ..................................................................................................135
Tableau 5.7 Plaque carré sandwich non symétrique simplement supportée sous charge
uniforme- Comparaison des déplacements. .............................................137
Tableau 5.8 Cylindre simplement supporté, flèche au point C. ...................................138
Tableau 5.9 Fréquences propres ( ) de vibration libre d’une aube de
compresseur. influence de maillage .........................................................140
Tableau 5.10 Fréquences propres ( ) de vibrations libres de panneaux
composites. ...............................................................................................143
Tableau 6.1 Influence du maillage sur la convergence du déplacement transversal w
..................................................................................................................149
Tableau 6.2 Influence du maillage sur la convergence du déplacement axial U.........149
Tableau 6.3 Erreurs sur le déplacement axial U pour les différents pas (10 éléments)
..................................................................................................................150
Tableau 6.4 Erreurs sur le déplacement transversal W pour les différents pas (10
éléments) ..................................................................................................150
Tableau 6.5 Panneau cylindrique articulé ( , maillage .................153
Tableau 6.6 Panneau cylindrique articulé ( , maillage .............154
Tableau 6.7 Panneau cylindrique articulé ( , maillage ............155
Tableau 6.8 Hémisphère pincé, déplacement ( ), maillage ..........................157
Tableau 6.9 Hémisphère pincé, déplacement ( ), maillage ........................158
xiii
RESUME
Mots clés:
Eléments finis, Mindlin Discret, Plaques, Coques, Zigzag, Non linéarité géométrique
xiv
ABSTRACT
This doctoral thesis deals with the finite element formulation and evaluation of
a modified Mindlin's discrete variational model for static, dynamic, linear and non-
linear composite plates and shells analysis. Including additional terms of zigzag type,
in order to improve the accuracy of stress, the model has been reformulated to take
into account the linear picewise of displacement variation. Consequently, two finite
plate and shell elements with four nodes, called DMQP and DMQS (Discrete
Quadrilateral Mindlin Plates and Shells respectively), enhanced by quadratic field
rotations, have been developed and validated under REFLEX and ABAQUS codes.
Both elements including the zigzag effect have been also developed in a second
version, and validated through several static and dynamic test problems known from
the literature, highlighting the independence towards the transverse shear correction
and in particular the stress accuracy with respect to the initial model without the
zigzag effect.
The satisfactory results of this model found through cases of linear isotropic
shell tests, motivated us to extend this approach to the non-linear geometric
applications. An isoparametric curve element of shell has been developed for this
purpose, where small elastic deformation assumptions of and large displacements and
moderate rotations are adopted. It is geometrically simple and has only four nodes at
corners and 6 DOF/node. The elementary calculation of the tangent stiffness matrix
consists in combining the linear part of the curved shell element (DMQS) with that of
the membrane Q4 non-linear part. An Updated Lagrangian Formulation at each
Iteration (ULFI) is used with Newton-Raphson resolution Method. Some standard
tests of nonlinear geometrical shell structures are presented; they show a very good
convergence and global behavior better than such elements.
Keywords:
Finite elements, Discrete Mindlin, Plates, Shells, Zigzag, geometric nonlinearity
xv
Introduction générale
1 Introduction
Le présent travail doctoral s’inscrit dans le cadre d’un programme d’échange et
de collaboration scientifique entre le groupe de recherche GRESPI/LMN (LISM par la
suite) de l'université de Reims et le laboratoire de Génie Mécanique de l’Université de
Biskra (Algérie). Le thème principal de la collaboration est le développement de
nouveaux éléments finis dits « spéciaux » à valeur ajoutée pour l’analyse du
comportement mécanique 2D et 3D des structures composites. Il s’agit pour notre part
de développer des éléments finis spéciaux d’ordre supérieur, plus simples et aussi
précis que ceux existants dans la littérature et dans les codes industriels de calcul des
structures.
Dans la thématique de modélisation, et jusqu’à l’avant dernière décennie, tous
les chercheurs ont surmonté, d’une manière ou d’une autre, les difficultés numériques
fréquemment gênantes rencontrées dans les modèles éléments finis issus des théories
correspondantes. L’examen des modèles 2D existants permit de conclure que la
tendance actuelle s'attachera à limiter le nombre de degrés de liberté nécessaire à une
modélisation fiable et réaliste avec des formes cinématiques ou mécaniques réputées
simples et robustes. Cela se traduit par le maintien des formes standards, à savoir le
quadrilatère et le triangle à quatre et trois nœuds respectivement, comme support de
tout développement futur d'un nouveau modèle élément fini. Par contre, cette plate
forme ne peut à leur seule représenter des variations brusques des contraintes
suivant l’épaisseur, du fait d'une discrétisation non fiable liée au manque des entités
géométriques (nœuds) à travers l'épaisseur et de la platitude des variables
cinématiques suivant l’épaisseur. Toutefois, l'exigence de décrire la forme parabolique
du profil des contraintes ramène la question sur les modèles mixtes ou les modèles
d'ordre supérieur qui nécessitent en eux-mêmes une certaine affinité entre les couches
qu'elle s'avère difficile à se concrétiser avec un seul point suivant la totalité de
l'épaisseur de la structure.
Ainsi, nous présentons dans cette thèse doctorale, une contribution à la
modélisation du comportement des plaques et coques multicouches/sandwichs avec un
développement des éléments finis basés sur l'aspect global des quelques approches
pour l'analyse linéaire et non-linéaire, statique et dynamiques des structures
isotropes et composites. Les premiers travaux ont consisté à élaborer une synthèse
1
Introduction générale
2 Objectifs de la thèse
Après une évaluation globale des études antérieures dans la littérature pour
l'analyse des coques, nous constatons que la plupart des formulations sont basées sur
des fonctionnels mixtes ou en déplacement avec des éléments finis d'ordre inférieurs.
Ces formulations sont limitées à l'analyse des matériaux isotropes avec peu
d'applications sur les composites stratifiés, leur performance est fiable et demeure
limitée aux réponses globales. Cela nous a motivé à adopter la notion des termes de
zigzag comme valeur ajoutée aux modèles développées récemment au niveau du
groupe de recherche GRESPI. Ce travail porte sur l'amélioration du modèle discret en
déplacement de Mindlin, en proposant deux éléments finis améliorées, pour le calcul
et l'analyse linéaire et non linéaire des structures plaques et coques composites. Plus
précisément, notre objectif dans ce travail se résume ainsi:
Valoriser le modèle DDM par la reformulation et la mise en œuvre de deux
éléments finis de plaques et coques multicouches DMQP et DMQS (riches en
flexion et projetés en cisaillement).
Pallier à l'insuffisance liée à la limite en précision des contraintes
(essentiellement du cisaillement transverse) inhérentes généralement aux
modèles du premier ordre, observée sur certains problèmes 2D/3D des
plaques/coques composites. Cela est traduit par l'insertion des termes de
zigzag, issus de l'approche de Murakami, dans le modèle initiale. Deux
éléments finis raffinés avec zigzag sont développés à cet effet.
Elargir le domaine d'application du présent modèle aux applications non
linéaires géométrique, en adoptant une formulation Lagrangienne Actualisée à
chaque Itération (FLAI).
2
Introduction générale
3 Plan de la thèse
La thématique de recherche consiste au développement des éléments finis
fondés sur le concept discret de Mindlin et qui tient en compte l'effet de zigzag comme
termes d'ordre élevé. Le mémoire de thèse débute par une introduction générale,
s'articule ensuite autour de six chapitres et s'achève par une conclusion et des
perspectives. Dans le premier chapitre, nous décrivons selon une synthèse historique
récente, les fameuses tendances concernant la formulation des théories de plaques et
coques composites et les éléments finis correspondantes. L'étude repose sur les
travaux réalisés durant la dernière décennie, faisant chaque fois le recours à leurs
descriptions authentiques fondées sur des théories antérieures. L’évaluation de la
performance et l’efficacité numériques de ces contributions demeure dans un cadre
général de comparaison et qui ne dépassent pas ce qu’ont été signalé par les auteurs
eux mêmes.
Le deuxième chapitre est consacré à la formulation théorique et l'aspect
variationnel du nouveau modèle de plaques DDM concrétisé par un élément fini
quadrilatère à quatre nœuds. Le champ de déplacement de cet élément est enrichi par
des termes quadratiques introduites dans les variables cinématiques de rotations, ce
qui y nécessite des entités géométriques sur les mi- cotés, généralement gênantes
pour l'ingénieur. Par conséquent, nous avons adopté deux hypothèses: l'une
cinématique et l'autre mécanique afin de rendre la géométrie de l'élément plus simple
et ne comporte que des nœuds aux sommets. D'autre part, en mesure d'améliorer la
réponse locale du modèle proposé, et qui possède toutefois un aspect du premier ordre.
Une attention particulière est portée à l'incorporation de la fonction de zigzag de
Murakami dans le but d'attribuer des termes supérieurs au modèle primitif.
Rappelons que la formulation de tous les modèles se traduit par une géométrie
identique. Dans cette partie, seulement la cinématique et les lois d'équilibre seront
attribués à ce nouveau modèle, en étalant à la fin, la formulation des contraintes et le
calcul de la matrice masse correspondantes.
Le troisième chapitre est consacré à la validation du nouveau modèle, en
mettant l'accent sur les structures multicouches stratifiées ou sandwich, avec
différentes dispositions des couches et des matériaux. La validation du présent
modèle sur le comportement vibratoire des structures réelles est mis en évidence par
des tests sur le carton ondulé (sandwich) qui a été abondamment commenté dans des
travaux précédentes [4, 6-8]. La réponse harmonique de quelques structures en
vibration libre est testée également pour les deux versions du modèle DDM.
Les résultats obtenus dans le chapitre 3 et le succès de la mise en œuvre du
nouveau concept du modèle DDM avec zigzag sur les plaques composites nous ont
donnés les ambitions d'adopter la fonction de zigzag dans le modèle de coque
correspondant. Cela est illustré dans le quatrième chapitre avec une description
détaillée de la géométrie et l'écriture tensorielle du gradient de déformation. La
formulation variationnelle décrite dans cette section nous a conduits à définir les
3
Introduction générale
4
Chapitre 1 Théories aux éléments finis des plaques et
coques composites. Synthèse bibliographique
1.1 Introduction
L’abondance des théories développées ces dernières décennies pour la
modélisation et l'analyse des plaques et coques composites est importante, des
théories ayant but de réduire le problème tridimensionnel à un problème
bidimensionnel et de respecter les conditions issues de la mécanique des milieux
continus. Selon ces exigences qui ne sont jamais réalisées simultanément se révèle
une pléthore de théories et de formulations qui s’adaptent bien à la modélisation par
éléments finis.
Les théories des plaques/coques des structures multicouches peuvent être
développées en faisant d’une part, le choix approprié de trois points suivants : la
méthode d’élimination de la variable d’épaisseur z (axiomatique, asymptotique), le
choix des variables inconnues (formulation contrainte, déplacement ou mixte) et
finalement selon la description (approximation) des variables (modèle monocouche
équivalente dit par fois modèle global, modèle par- couche, modèle global- local…etc.).
, et d'autre part les différentes conditions cinématiques et mécaniques à remplir au
niveau des interfaces, à l’intérieur des couches et inter éléments.
Ainsi, leur classification devient de plus en plus compliquée et par fois
chevauchante. Dans un bref aperçu, on a constaté deux aspects fondamentaux faisant
la différence entre les différentes approches ainsi que leur classification : la
description des variables (déplacements ou contraintes) suivant la direction de
l’épaisseur, et leur aptitude à satisfaire les exigences de continuité. Selon ces critères,
plusieurs combinaisons ont été faites, citant, en premier lieu, les théories globales en
déplacement, dont les variables sont introduites pour toute la plaque/coque, d’où ils
font généralement référence à leur surface moyenne. Elles se scindent en: les théories
classiques des multicouches, les théories de déformation en cisaillement du premier
ordre, les théories de Valsov- Reddy et les théories de déformation en cisaillement
d’ordre supérieur. D’autre coté, des théories globales améliorées provenant de la prise
en compte de la condition de continuité des déplacements dans les théories globales.
Une autre théorie appelée «théorie de zig zag» dérive également du modèle global, en
imposant la continuité des contraintes sur les interfaces, tout en vérifiant que les
faces supérieures et inférieures sont libres en contraintes de cisaillement, des autres
théories sont développées par la suite, nommées les théories par- couche, dans
lesquelles, chaque couche est vue comme une plaque/coque indépendante. Par
5
Chapitre 1. Théories aux éléments finis des plaques et coques composites.
Synthèse bibliographique
6
Chapitre 1. Théories aux éléments finis des plaques et coques composites.
Synthèse bibliographique
7
Chapitre 1. Théories aux éléments finis des plaques et coques composites.
Synthèse bibliographique
classique des plaques (CPT)2 qui se base sur les fameuses hypothèses de kirchhoff. En
conséquence, des approches de modélisation qui exigent la prise en compte des
déformations du cisaillement transverse ont constitué le sujet d'une sérieuse
recherche dans les quatre dernières décennies. Il convient de mentionner ici les
travaux de Bert [15], Noor et al. [16, 17].
Quelques approches de modélisation sont des extensions d'approches similaires
à celles utilisées dans les plaques isotropes. Celles appartenant à cette catégorie
nommées « Smeared Laminate Models », sont des modèles basés sur des hypothèses
globales (pour la totalité de la stratification). Ils sont caractérisés par une distribution
linéaire ou non-linaire des déplacements plans ou des contraintes de cisaillement
transverse suivant la direction de l'épaisseur. Ces approches considèrent toutes les
couches comme une monocouche anisotrope équivalente [18].
Citant par exemple les deux théories largement discutées dans la littérature
spécialisée, et qui sont fondées sur une distribution linéaire globale: la théorie
classique de plaques (CPT) basée sur la négligence des déformations de cisaillement
transverse et des déformations et des contraintes normales [19] ainsi que la théorie de
déformation en cisaillement du premier ordre (FSDT3) dont les déformations du
cisaillement transverse sont supposées constantes a travers l'épaisseur [20] et la
composante de déformation normale et celle de la contrainte transverse sont
négligées. En revanche les théories de déformation en cisaillement d'ordre supérieur
(HSDT4) constituent de leur part des exemples de modèles globales basés sur une
distribution non- linéaire [21-28]. Du point de vue conditions de contact, des
contraintes de cisaillement incompatibles entre toutes les couches adjacentes
résulteraient comme conséquence.
z
z
ESLM LW
8
Chapitre 1. Théories aux éléments finis des plaques et coques composites.
Synthèse bibliographique
(1.1)
9
Chapitre 1. Théories aux éléments finis des plaques et coques composites.
Synthèse bibliographique
(1.2)
Nous remarquons que les variables additionnelles, par rapport au cas du FSDT, sont
introduites dans le développement du déplacement transverse.
A fin de s'affranchir de limitations liées aux théories de déformation en
cisaillement du premier ordre qui considèrent une variation constante de déformation
transverse à travers l’épaisseur de la plaque, et qui nécessitent des facteurs
arbitraires de correction, les théories d'ordre supérieur ont été développées. Elles
impliquent des termes d'ordre élevé correspondent au développement de Taylor dans
les champs des déplacements, chaque terme additionnel correspond à une variable
dépendante. Contrairement aux théories du premier ordre, ces théories sont fondées
sur une variation cubique (ou plus) des déplacements dans le plan afin d’obtenir une
variation parabolique (ou plus) des déformations de cisaillement suivant l’épaisseur
de la plaque. Dans le cas général, selon Carrera [29], les théories d’ordre supérieur
sont basées sur des champs de déplacements de la forme suivante :
(1.3)
Avec Ni sont les ordres de développements utilisées dans les variables de
déplacements. Les champs de déplacements relatifs à la théorie HSDT se sont
reportés sur la figure 1.3.
Hildebrand et al. [35] ont eu la primeur d'introduire cette approche dans les
théories des plaques et coques. Nelson et Lorch [36], Librescu [37], ont présenté des
déplacements d'ordre élevé basés sur les mêmes théories pour l'analyse des plaques
multicouches. Lo et al. [28, 38] ont proposé une solution approchée pour les plaques
multicouches selon le modèle en déplacement d'ordre supérieur et qui prend en
considération l'effet de la déformation normale transverse. Levinson [23] et Murthy
[24] ont proposé également des théories de déformation en cisaillement d'ordre trois,
dont les équations d'équilibre ne sont pas variationnellement compatibles avec les
déplacements. D'après Bickford [39] et Reddy [25, 40], recourant dans leurs analyses à
la théorie celle du premier ordre déjà utilisé par Whitney et Pagano [20], les
déplacements proposés ont été de point de vue cinématique, les plus fréquents dans
l'ensemble des déplacements d'ordre supérieur [41].
Ayant basé sur la cinématique développée par Levinson, Bickford [39] a
proposé une théorie des poutres d'ordre supérieur qui contourne le désavantage de
l'approche originale. Une autre forme cohérente avec les modèles cinématiques
originalement proposés par Levinson et Murthy a été formulée ultérieurement par
Reddy [25]. Elle est complètement simple, respecte les conditions de nullité des
contraintes de cisaillement sur les faces supérieures et inférieures de la plaque, et
n’inclut plus des variables additionnelles. Cette théorie n’est en quelque sorte que
l’extension du raffinement de la théorie de Reissner-Mindlin établi par Valsov [34].
Elle est communément connue par la théorie de Valsov-Reddy du troisième ordre dont
les déplacements sont ainsi :
10
Chapitre 1. Théories aux éléments finis des plaques et coques composites.
Synthèse bibliographique
(1.4)
3,z
3,z
u3, x
13
u3, x
1, x, 1, x,
u1
90°
u1
90°
FSDT
CLT
u1 ( x, z ) u1 z 1
u1 ( x, z ) u z u3,x
13 u 3, x
1
3,z
u1
HSDT
u1 ( x, z) u1 z i u1i , i 1, n
Le lecteur peut se référer au travail de Noor et Burton [16], dans lequel ils ont
présenté une liste complète des références concernant les théories de déformation en
cisaillement du premier et d'ordre supérieur. Elles comportent l'analyse statique,
dynamique et de flambement des structures composites multicouches.
11
Chapitre 1. Théories aux éléments finis des plaques et coques composites.
Synthèse bibliographique
Rohwer [46] a regroupé dans une étude comparative diverses théories d'ordre
supérieur pour l'analyse de flexion des plaques composites multicouches, où il a met
en valeur les avantages et les inconvénients de différentes théories. Les investigations
numériques complètes faites par l'auteur ont montré que la théorie de Murthy et celle
de Reddy sont plus avantageuses parmi les théories d'ordre supérieur. Celle de Reddy,
demeure la plus fameuse et la plus simple dans l'analyse des plaques composites,
selon Ferreira [47].
Reprenant la démarche de Kant [48], où il a proposé autrement une
formulation analytique et des solutions basées sur des théories raffinées servant à
l'analyse des problèmes statiques des plaques multicouches et sandwich simplement
supportées. Cette théorie s'inspire de la solution d'élasticité, qui reflète une variation
parabolique suivant l'épaisseur de la plaque de contraintes de cisaillement, et qui
implique des fonctions cubiques des déplacements plans, dont une forme non linéaire
de la déformation normale transverse est exigée. Une comparaison a été faite avec des
solutions d'élasticité tridimensionnel montre que cette théorie prévoit le déplacement
transversal et les contraintes planes d'une façon plus fiable que les autres théories
mentionnées dans son écrit.
Sur la base de la même théorie, Matsunaga [49] s'est investi dans l'analyse de
distributions des contraintes et de déplacement dans les plaques stratifiées croisées
(cross-ply) et sandwich, soumises à une pression latérale. Les contraintes normales et
de cisaillement sont obtenus par intégration des équations d'équilibre tridimensionnel
suivant la direction de l'épaisseur. L'ordre de l'approximation de la théorie peut être
estimé en corrélation avec le rapport du coté de la plaque et son épaisseur d'une part,
et le nombre de couches et l'ordre d'empilement d'autre part. Le travail traite l'impact
de différentes conditions de continuité du déplacement sur la théorie
tridimensionnelle par couche de Pagano [50]. La continuité de tous les déplacements
(plans et normaux) sur les interfaces est de type C0, et celle de monocouche utilisée
(ordre supérieur est généralement supérieure à C1.
Selon la même démarche et dans un travail subséquent de comparaison avec
les théories 3D par-couche, l'effet des charges thermiques a été cette fois-ci étudié.
Matsunaga [51] a évalué les contraintes inter laminaires et les déplacements dans les
mêmes structures. L’auteur a montré que la théorie adoptée (2D monocouche) est
entièrement acceptable. S’en rapportant à Wu et al. [52], ils ont constaté d’autant que
la théorie adoptée comporte un nombre des variables indépendant de celui de couches,
elle ne satisfait pas les conditions de continuité de contraintes du cisaillement
transverse entre les couches. Même remarque a été signalée sur la base des tests
numériques effectués dans le travail de Zenkour [53]. Les solutions analytiques
présentées sont basées sur la théorie globale d'ordre supérieur, elles correspondent
aux plaques stratifiées croisées.
Swaminathan et Patil [54] ont présenté de leur part des formulations et des
solutions analytiques pour l'analyse statique des stratifiés alternés antisymétriques
(antisymitric angle-ply) et des plaques sandwichs simplement supportées. Les
12
Chapitre 1. Théories aux éléments finis des plaques et coques composites.
Synthèse bibliographique
5 Cisaillement transverse
13
Chapitre 1. Théories aux éléments finis des plaques et coques composites.
Synthèse bibliographique
non exhaustive d’éléments finis du premier ordre, en particulier pour les plaques
épaisses en flexion/CT, vaut ainsi d’être citée :
Les éléments DSQ et DST de Lardeur [63] sont basés sur un modèle
variationnel mixte modifié. L’auteur calcule la déformation de CT en partant des
équations d’équilibre ( constant, cubique, quadratique). Ensuite, il ajoute à
chaque côté de l’élément deux variables cinématiques et deux rotations au
milieu . Il obtient ainsi huit rotations et quatre déplacements transverse pour
un quadrilatère qui seront éliminés par les équations d’équilibre. Les deux éléments
sont à labri de verrouillage et ne présentent pas de modes parasites. Cependant, la
formulation théorique reste complexe et nécessite le calcul des facteurs de correction
pour chaque type de problème [6] ;
Selon une formulation mixte en contrainte de type Hellinger-Reissner,
Bouabdalallah [64] a formulé un élément quadrilatéral à 4 nœuds. Toutes les
contraintes sont définies dans un seul vecteur et s'expriment en termes
indéterminés . Les modes parasites ont été contrôlés par au moins 14 paramètres et
qui sont éliminés par la suite par condensation statique au niveau élémentaire.
Quoique l'élément a donné des bons résultats, cette technique reste limitée à
l’application aux structures courbes essentiellement cylindriques et fait recours
d'autre part aux facteurs de correction du CT pour bien simuler le comportement
mécanique de ce type des structures.
L’élément fini proposé par Katili [1] pour les plaques est basé sur la théorie des
plaques épaisses de Reissner/Mindlin (théorie du premier ordre). La fonctionnelle
utilisée est de type mixte modifié de Hu-Washizu. Les variables cinématiques (w,
βx,et βy) de continuité C0 sont interpolées de manière quadratique en introduisant des
variables , au milieu du côté k. Ces variables sont éliminés par une technique
utilisant les équations d’équilibre. L’élément fini, baptisé DST-BK surmonte le
verrouillage en CT et ne présente aucun mode parasite. Néanmoins, la correction du
cisaillement transverse est un calcul incontournable.
Auricchio [65] a formulé un élément à 4 nœuds et 5ddl/nœud, baptisé ELM4.
Celui-ci est obtenu via une approche mixte améliorée. La formulation variationnelle
adoptée introduit des modes incompatibles pour améliorer la précision sur les
déformations planes.
L’élément « CTMQ20 » de Cen [66], est un élément fini simple de plaque en
flexion/CT à 4 nœuds et 5ddl/nœud, basé sur la théorie du premier ordre. Il est
proposé pour l’analyse des plaques composites et stratifiées à disposition arbitraire.
Cet élément est construit selon la procédure suivante : Les fonctions de variation de la
rotation et de la déformation du CT le long de chaque côté de l’élément sont
déterminées à l’aide de la théorie des poutres de Timochenko. Les champs
élémentaires des déplacements plans et des rotations ainsi que le champ des
déformations du CT sont déterminés par l'utilisation d'une technique d’interpolation
14
Chapitre 1. Théories aux éléments finis des plaques et coques composites.
Synthèse bibliographique
améliorée. En outre, une procédure hybride simple est également proposée pour
améliorer les solutions de contrainte.
Belinha et al. [67] ont étendu puis utilisé la méthode EFGM (Element Free
Galerkin Method) pour l’analyse des plaques anisotropes et stratifiées, en considérant
la théorie de Reissner-Mindlin (FSDT). Les fonctions d’approximation sont calculées
sur la base d’une approche au sens des moindres carrés (MLS: Moving Least Square).
Le verrouillage en CT est évité en utilisant de polynômes d’interpolation appropriés.
La continuité étant bien assurée selon les auteurs.
Moleiro et al. [68] ont récemment présenté un modèle d’élément fini mixte à 4
nœuds et 5 ddl/nœud pour l’analyse statique des plaques composites stratifiées. La
formulation est basée sur un principe variationnel au sens des moindres carrés que
les auteurs considèrent comme une approche alternative aux modèles éléments finis
issus de formulations variationnelles faibles classiques.
Suite à la série des éléments finis proposés et développés depuis 1993 par Ayad
[57-60, 69] pour les structures plaques et coques isotropes, une tenue importante a été
consacrée aux structures composites. En se basant sur le modèle mixte projetée en
cisaillement (MiSP : Mixed with Shear Projection), le développement se poursuivait.
Dans un premier lieu, un élément fini mixte-hybride baptisé MiSP4/ml a été
développé par l’auteur et son groupe [6], [70]. C’était une adaptation du modèle MiSP4
isotrope au cas des plaques stratifiées et sandwich. Cela a permis à ce modèle de
gagner du terrain dans la modélisation des plaques composites, quelque soit le
nombre de couches et sans recourir aux facteurs de correction. L’élément fini proposé
est de forme quadrilatéral à 4 nœuds et 3 ddl/nœud. Il est libre de tout verrouillage en
CT et présente des performances appréciables en matière de précision,
comparativement aux éléments issus des théories d’ordre supérieur, souvent
complexes et coûteux en temps de calcul. Les développements qui viennent en second
lieu étaient deux modèles d’éléments finis isoparamétriques courbes à 4 nœuds pour
l’analyse des coques composites multicouches. Le premier modèle est mixte hybride
naturel, baptisé NHMiSP4/ml (Natural Hybrid Mixed with Shear Projection 4-
node/multilayer). La partie membrane est complètement hybride naturelle, au sens de
Pian modifiée, et la partie flexion/CT est mixte-hybride, offrant toutes les deux un
avantage de taille à l’élément, en terme de précision sur les contraintes, avec une
sensibilité moindre aux distorsions de maillages. Le second modèle, baptisé
HMiSP4/Q4/ml, développe une formulation en déplacement pour la membrane et une
formulation mixte-hybride en flexion/CT. Ce modèle a été formulé pour tester son
aptitude à reproduire de bons résultats en dépit d’une représentation moins riche de
la membrane. Un des avantages de ces deux modèles de coques est la prise en compte
dans la formulation théorique du gauchissement dans les coques. Leur adaptation aux
coques multicouches suit la même démarche que pour la plaque multicouche. Les
facteurs de correction du CT sont également évités.
Basant sur l’hypothèse du Mindlin [32] proposée pour les plaques isotropes, un
nouveau modèle variationnelle en déplacement dit de Mindlin discret (DDM:
15
Chapitre 1. Théories aux éléments finis des plaques et coques composites.
Synthèse bibliographique
Displacement Discrete Mindlin) est développé également par Ayad [3]. Sa projection a
été faite sur deux éléments finis de plaque et de coque courbe géométriquement
simple (4 et 3 nœuds) [3]. Ils prennent en compte de l’effet du cisaillement transverse
à travers l’épaisseur. Le modèle DDM introduit de manière discrète deux hypothèses
de Mindlin. La première hypothèse est cinématique, elle consiste à introduire sous la
forme d’une intégrale de contour une équation cinématique de la déformation de CT.
Elle permet l’élimination du verrouillage en CT sans introduire des fonctions bulles ou
sans recourir à l’intégration réduite ou sélective. Il s’agit de l’approche des
déformations de CT de substitution, connue sous le nom de « ANS method : Assumed
Natural Strains ». La seconde hypothèse fait appel à deux lois de comportement, l’une
en flexion et l’autre en CT, et deux équations d’équilibre d’une plaque en flexion/CT.
Elle a pour principal avantage une élimination locale des degrés de liberté de rotation,
introduits initialement au milieu des bords élémentaires par le biais d’une
approximation quadratique des rotations de la normale à la surface moyenne. Ces
hypothèses ont été modifiées localement afin d'être utilisées pour les plaques/coques
stratifiées composites. Le modèle de plaque dérivé, baptisé DDM, a été appliqué au
début sur un triangle à 3 nœuds (DMTP: Discrete Mindlin Triangle for plate) [71],
ensuite dupliqué sous un autre élément à quatre nœuds avec 5 ddl/nœud pour les
coques composites. Les tests ont été effectués sur des problèmes standards, statique et
dynamiques. Ils montraient des réponses globales précises avec un cout réduit, qui
rend le modèle DDM comme alternative aux éléments finis de Kirchhoff discrets,
connu pour leur simplicité et leur robustesse.
Une reformulation en contraintes hybrides est adoptée par Daghia et al. [72]
comme généralisation de celle présentée par De Miranda et Ubertini [73] pour les
plaques isotropes. L’approximation du déplacement transverse est améliorée en y
reliant avec les rotations nodales, une approximation convenable est choisie
cependant pour les contraintes résultantes. Les contraintes transverses à travers
l'épaisseur du multicouche sont reconstruites a posteriori en utilisant l'équilibre
tridimensionnel. Pour améliorer les résultats, les contraintes résultantes entrant dans
le processus de reconstruction sont d'abord récupérées à l'aide d'une procédure basée
sur le rapiéçage super-convergent (superconvergent patch-based procedure) nommée
recouvrement par compatibilité dans les patchs (Recovery by Compatibility in
Patches), et qui est bien étendue dans cette contribution pour les plaques
multicouches. Ce traitement préliminaire est très efficace du point de vue calcul et
généralement utile soit pour évaluer avec précision les contraintes résultantes soit
pour estimer l'erreur de discrétisation. En effet, dans le même contexte, il joue
également un rôle clé dans la prédiction, de manière efficace, les profils de contraintes
du cisaillement. Il garantit la convergence globale de l’ensemble de la stratégie de
reconstruction, qui ne nécessite pas de correction pour s’accommoder aux défauts
d'équilibre.
Ainsi, comme extension de l'élément 9βQ4, proposée récemment par De
Miranda et Ubertini [73], un nouveau élément a été reformulé afin d'inclure le
16
Chapitre 1. Théories aux éléments finis des plaques et coques composites.
Synthèse bibliographique
17
Chapitre 1. Théories aux éléments finis des plaques et coques composites.
Synthèse bibliographique
Le lecteur peut se référer aux articles [17, 81], dans lesquels les éléments ne
satisfont pas toutes les exigences dans la liste précitée; Ils ne sont pas formulé
principalement sur le critère de C0 continuité. Pour Carrera, seulement les tentatives
suivantes peuvent être mentionnées à l'égard des éléments finis HSDT qui incluent la
condition C0 : Le premier est celui réalisé par Pandya et Kant [82-84], dont les auteurs
ont étendu la théorie élaborée par Kant [45] pour présenter des formulations éléments
finis de type C0, respectivement pour les plaques stratifiées symétriques et non
symétriques et les panneaux sandwichs. Dans ces travaux, les résultats ont été
améliorés par rapport au modèle Reissner-Mindlin standard, ayant utilisé un
développement d'ordre supérieur pour le déplacement à travers l'épaisseur. L'élément
résultant représente une amélioration vis à vis aux résultats de FSDT mais n'inclut
pas l'effet du zigzag et n'accomplit pas la continuité inter- laminaire, à savoir, les
conditions d'équilibre pour les contraintes sur chaque interface.
Dans ce cadre, Phan et Reddy [85] ont fait une tentative d’améliorer un
élément conforme des plaque minces. L'élément de base est formulé par Bogner et al.
[86], étant rectangulaire à quatre nœuds d’où les degrés de liberté par nœud sont :
, et (le champ de est une interpolation bi- cubique).
Concernant l'élément développé par Phan et al. [85], les variables des champs
supplémentaires utilisées sont les déplacements plans , et les rotations de la
normale , . Elles sont rapprochées par des fonctions bi-
linéaires. Comme les déformations en cisaillement transverse varient à travers
l’épaisseur de la plaque, et correspondent à celles du plan moyen de la plaque
[87]. Ainsi, Ils ont fait une tentative dans laquelle les composantes de déformation en
cisaillement transverse sont prises au lieu des rotations de la normale .
L’élément, quoiqu' il est assez simple, sa forme ne le permet pas d’être utilisé que
dans des plaques de formes rectangulaires.
Basant sur la théorie de Reddy d'ordre supérieur, Nayak et al. [88] ont formulé
deux éléments finis (quatre et neuf nœuds) avec sept d.d.l/noeud. Les auteurs ont
adopté une interpolation propre à la déformation du cisaillement transverse. Le
phénomène de blocage a été maîtrisé par les techniques de Dvorkin et Bathe [89], et
Huang et Hinton [90] respectivement pour les éléments à quatre et neuf nœuds. Les
éléments sont validés avec des tests de vibration libre des plaques multicouches et
sandwich. Ceci a montré l'applicabilité de tous les éléments de type Reddy à l'analyse
dynamique. Les mêmes auteurs [91] ont exploité par la suite l'élément à neuf nœuds
dans des analyses du flambage et de vibration libre des plaques multicouches et
sandwich soumises à des contraintes initiales.
Le succès de la dernière théorie a permis à Sheikh et Chakrabarti [87] de
formuler un élément triangulaire d'ordre supérieur, à six nœuds et six ddl/nœud. La
formulation est basée sur des champs de déplacement indépendants. De point de vue
numérique, l'élément s'avère entièrement performant.
Khare et al. [92] ont formulé un élément coque d'ordre supérieur à facette
plane, avec neuf degré de liberté par nœud, en se basant sur un élément quadrilatère
18
Chapitre 1. Théories aux éléments finis des plaques et coques composites.
Synthèse bibliographique
19
Chapitre 1. Théories aux éléments finis des plaques et coques composites.
Synthèse bibliographique
20
Chapitre 1. Théories aux éléments finis des plaques et coques composites.
Synthèse bibliographique
couche k+1
couche k
face supérieure
(Top)
kxxb1
kxxb1 kzxb1
k kzyb1 kzzb1
xxt
kzyb face inférieure
kzzt (bottom)
k
zxt
kxxt
Figure 1.5 Notation détaillée des contraintes sur les interfaces [29].
21
Chapitre 1. Théories aux éléments finis des plaques et coques composites.
Synthèse bibliographique
originales de ces trois scientifiques. Une note historique, récemment signalée par
Carrera [105], dont il a proposé d'authentifier les trois principales théories selon leurs
fondateurs, ainsi : Lekhnitskii-Ren, Ambartsumian-Whitney-Rath-Das et finalement
Reissner-Murakami-Carrera.
(1.5)
Quatre fonctions indépendantes de x, y ont été introduites pour décrire les contraintes
du cisaillement transverse. Les constantes de couche sont fonction
parabolique de la coordonnée de l'épaisseur z. Leurs formules pourraient être
mentionnées dans le travail original. La formulation (relation 1.5) représente un
champ continu de contrainte du cisaillement transverse, qui est parabolique dans
chaque couche. Les champs du déplacement sont obtenus en intégrant la relation du
déformation/déplacement et en imposant des conditions de compatibilité pour le
déplacement à l'interface. Le champ du déplacement final prend la forme de zigzag
suivante:
22
Chapitre 1. Théories aux éléments finis des plaques et coques composites.
Synthèse bibliographique
(1.6)
(1.7)
23
Chapitre 1. Théories aux éléments finis des plaques et coques composites.
Synthèse bibliographique
(1.8)
(1.9)
(1.10)
Il s'avère clairement que les déplacements et qui sont linéaires par couche
engendrent des coefficients et indépendantes de la coordonnée z, (relation 1.10).
En effet, les conditions homogènes des contraintes de cisaillement transverse sur les
faces supérieures et inférieures ne sont pas satisfaites.
(1.11)
suivant la même démarche et dans un travail ultérieur [116], Di Sciuva a proposé une
formulation d'ordre élevé qui permet une variation non- linéaire (cubique) de
24
Chapitre 1. Théories aux éléments finis des plaques et coques composites.
Synthèse bibliographique
(1.12)
Avec
(1.13)
Dans le même contexte d’une théorie semblable, baptisée EHOPT6, avec un nombre
des variables préliminaires identique à celui du premier ordre, Cho et al. [117] ont
combiné la théorie de plaque d’ordre supérieur de Reddy [25] et celle du zigzag par
couche améliorée par Di Sciuva [118]. Ainsi, ils ont proposé une formulation (équation
1.14) qui reflète une variation cubique des déplacements plans et qu’elle prend en
compte outre la distribution parabolique de déformation transverse, la forme zigzag
linéaire. Cette formulation vérifie la discontinuité des déformations entre les couches
(selon l'exigence de continuité des contraintes), où les déformations normales
transverses sont négligées.
(1.14)
(1.15)
Les conditions aux limites des contraintes de cisaillement nulles sur les faces
supérieures et inférieures seraient satisfaites en écrivant:
(1.16)
25
Chapitre 1. Théories aux éléments finis des plaques et coques composites.
Synthèse bibliographique
(1.17)
7 Eléments finis
26
Chapitre 1. Théories aux éléments finis des plaques et coques composites.
Synthèse bibliographique
27
Chapitre 1. Théories aux éléments finis des plaques et coques composites.
Synthèse bibliographique
z
-1 1
z4=0 h4
z3=0 h3 z=0 Ω
z2=0 h2
z1=0 h1
(1.19)
(1.20)
28
Chapitre 1. Théories aux éléments finis des plaques et coques composites.
Synthèse bibliographique
4z 3 4z 3 8z 3
u u0 w,x z φx ζk ( 1)k u M
3h2 3 h 2 3 h h 2
k
4z 3
4z
3
8z 3
v v0 2 w, y z 2 φy ζ k ( 1)k v M
2
(1.22)
3h 3h 3h k h
w w0
29
Chapitre 1. Théories aux éléments finis des plaques et coques composites.
Synthèse bibliographique
30
Chapitre 1. Théories aux éléments finis des plaques et coques composites.
Synthèse bibliographique
En 2003, Topdar, Sheikh et Dhang [134] ont constaté, pour formuler leur
élément fini des plaques multicouches et sandwich, que la théorie proposée par Cho et
Parmerter [117] (relation 1.14) est la plus adéquate à adopter parmi les théories
raffinées de déformation en cisaillement d’ordre supérieur (RHSDT). Telle théorie
satisfait complètement l’exigence de la continuité inter-élément (continuité C1 du
déplacement transverse le long des cotés de l’élément, comme c'est exigée dans la
théorie de Reddy du troisième ordre). Le concept de l’élément conforme des plaques
minces proposé par Bogner et al. [86] a été utilisé également pour développer cet
élément rectangulaire à quatre nœuds, d’où il reflète des bonnes performances en
termes d’exactitude des résultats.
Dans la même année, et dans le but de valoriser la théorie de Reddy d’ordre
supérieur avec une formulation en déplacement, Cheikh et Chakrabarti [87] ont
montré que la technique adoptée, est ainsi applicable pour une forme triangulaire
avec un élément conforme des plaques minces, ayant trois nœuds aux sommets et
trois aux mi- cotés avec des degrés de liberté indépendants : u, v, w, x , y , x et y. De
point de vue numérique, l’élément est entièrement performant, libre du problème de
blocage en cisaillement et ne possède aucun mode parasite. L’auteur a adopté cette
théorie pour la modélisation des plaques multicouches, avec une application très
ample. En revanche, il a retourné dans quelques contributions ultérieures - dans les
travaux réalisés en collaboration avec Pandit et al. [135, 136] - à l’approche de Cho et
Parmerter concernant les formulations qui se rapportent aux plaques sandwich à
cœurs légèrement compressible.
31
Chapitre 1. Théories aux éléments finis des plaques et coques composites.
Synthèse bibliographique
Les modèles déplacement définis par l’équation (1.23) requièrent l’adoption des
critères stricts afin de satisfaire les conditions de compatibilité sur chaque interface.
Une généralisation de ces types de théories a été donnée par Nosier et al. [143] et par
Reddy [43], qui ont exprimé les variables de déplacement dans la direction de
l’épaisseur en termes des polynômes de Lagrange,
32
Chapitre 1. Théories aux éléments finis des plaques et coques composites.
Synthèse bibliographique
Les valeurs sur l'interface sont utilisées comme variables inconnues (les deux
premiers termes du développement, c.-à-d. et ), ce qui permet facilement
de faire une liaison entre les conditions de compatibilité sur chaque interface. En
effet, et coïncident avec les polynômes linéaires de Lagrange, cependant
devraient être une base indépendante de polynômes qui commence par le facteur
(parabolique).
33
Chapitre 1. Théories aux éléments finis des plaques et coques composites.
Synthèse bibliographique
34
Chapitre 1. Théories aux éléments finis des plaques et coques composites.
Synthèse bibliographique
Donc, les théories par couche sont des techniques prometteuses qu'ont été développées
afin d'améliorer l'exactitude des contraintes transverses dans les structures
multicouches. Elles sont basées sur les hypothèses en déplacements pour chaque
couche individuelle, d'où les distributions globales et locales de déplacements et de
contraintes sont bien représentées. Prenant le cas de la reformulation faite par Lu
and Liu [170] (la relation 1.27), les champs de déplacements s'écrivent explicitement
pour la couche :
u αk ( x , y, z ) u0kα ( x , y) u1kα ( x , y)z u2kα ( x , y)z 2 u3kα ( x , y)z 3 , α 1,2
(1.28)
w( x , y, z ) w 0 ( x , y) w1 ( x , y)z w2 ( x , y)z 2
Bien que cette théorie donne d’excellents déplacements et contraintes dans et hors
plan, elle souffre contrairement aux théories d'ordre supérieur, du nombre total de
degrés de liberté qui dépend du nombre de couches. Cet inconvénient peut ne pas être
gênant dans quelques formulations analytiques appliquées dans l'analyse des plaques
simples, cependant, il constitue un critère important dans la sélection d'une telle
théorie dans l'analyse des structures complexes par des méthodes de discrétisation.
Comme alternative, en 1995, Li et Liu [169] ont fait un compromis dans leur
formulation entre les théories par couche et d'ordre supérieur, dans laquelle un
nombre limité des termes dépendants de couches est considéré. Cela s'exprime par les
relations suivantes:
35
Chapitre 1. Théories aux éléments finis des plaques et coques composites.
Synthèse bibliographique
Li et al. ont procédé par la suite à assurer les conditions inter laminaires et
d'homogénéité (relations 1.30 et 1.31) afin de minimiser le nombre des variables
dépendantes du nombre de couches. La démarche est affectée aux stratifications
croisées (cross-ply) seulement en utilisant les équations constitutives
contraintes/déformations dans chaque couche. Le nombre total de degrés de liberté se
réduit finalement à trois.
uk 1
z zk
uk
z zk
τ xz
k 1
τ xz k
z zk
avec k 2,3,4,..., n (1.30)
z zk
h
τ xz 0 (1.31)
2
Les déplacements sont représentés par des polynômes impliquant des termes des
couches indépendants (relation 1.32), dont les coefficients sont fonctions des
propriétés de couches et des coordonnées.
uk u0 ( R1k A1k z )u1 ( R2k A2k z z 2 )u2 ( R3k A3k z z 3 )u3
( R4k A4k z )w0,x ( R5k A5k z )w1,x ( R6k A6k z )w2,x (1.32a)
wk ( x , y, z ) w0 ( x , y) w1 ( x , y)z w2 ( x , y)z 2
Avec
Q55
1
2 k 3 k
A1k
Q55k
, A2
k
Q55k j 2
α j z j , A k
3 α j z 2j , A4k A1k 1,
Q55k j 2
(1.32b)
1 1 k
A5k A2k , A6k A3k , Rik ( Ail 1 Ail )zl , αk Q55k 1 Q55k
2 3 i 2
36
Chapitre 1. Théories aux éléments finis des plaques et coques composites.
Synthèse bibliographique
et
u Lk ( x , y, z ) u 0k u1k ξ k
(1.35)
v Lk ( x , y, z ) v 0k v1k ξ k
37
Chapitre 1. Théories aux éléments finis des plaques et coques composites.
Synthèse bibliographique
z
nème couche kème couche
ζ
h/2
zk+1 zk k
x
h/2
1ière couche
v vG v v
k k ~ k
(1.36)
L L
w wG
La composante locale possède deux groupes, l’un ( u L ou v L ) ayant deux termes locaux,
k k
l’autre ( u~L ou ~v L ) avec un seul terme local. Prenant le cas dans lequel le premier
k k
groupe consiste en deux termes du premier et deuxième ordre, alors que le deuxième
groupe comporte le terme du troisième ordre. La présente combinaison représente
donc la théorie de double-superposition 1,2-3, dont les composantes locales s’écrivent
ainsi :
u Lk ( x , y, ξ k ) u1k ( x , y )ξ k u 2k ( x , y)ξ k2
v Lk ( x , y, ξ k ) v1k ( x , y)ξ k v 2k ( x , y)ξ k2
(1.37)
~ k ( x , y, ξ k ) u k ( x , y)ξ 3
u L 3 k
38
Chapitre 1. Théories aux éléments finis des plaques et coques composites.
Synthèse bibliographique
Dans la théorie de double-superposition 1,2-3 [172], Il est constaté que le nombre total
de variables est indépendant du nombre de couches. En comparant avec la théorie de
déformation en cisaillement du même ordre, il y a six variables additionnelles,
indépendantes de couches. Toutefois, quand elle est comparée avec les théories par-
couche, les coefficients d’ordre supérieur (en termes de z) sont liées aux propriétés
propres des couches et à la coordonnée z suivant l’épaisseur au lieu de variables
inconnues. En plus, basant sur les résultats antérieurs [169-171], les distributions des
déplacements plans à travers l’épaisseur prendraient une forme zigzag cependant les
contraintes de cisaillement transverse pourraient être dépendantes de couches.
En utilisant la technique de la superposition des composantes du déplacement
globales et locales, plusieurs travaux d’analyse et de développement sur la présente
théorie pour les composites stratifiées ont été effectués [56, 168, 173-176].
Néanmoins, Wu et al. [174] ont noté que la théorie globale-locale du troisième
ordre (à savoir la théorie 1, 2-3 en double superposition, appelée par Li et Liu [172],
«les théories généralisées des structures multicouches basées sur l'hypothèse de
double superposition») éprouverait des difficultés à prévoir les contraintes de
cisaillement transverses avec précision si le nombre de couches est supérieur à cinq.
En particulier, l'exactitude de cette théorie diminue lorsqu'elle est appliquée à
l'analyse des plaques multicouches non-symétriques. Cependant, ils ont constaté que
d'autres théories globales-locales d'ordre supérieur peuvent présenter des contraintes
planes adéquates et permettre d'avoir des contraintes de cisaillement transverse
satisfaisantes par l'emploi directe des équations constitutives. L'utilisation des
champs de déplacement selon les théories «global-local» a acquis toutefois un large
domaine d'application dans d'analyse des composites multicouches [56, 168, 173-176].
39
Chapitre 1. Théories aux éléments finis des plaques et coques composites.
Synthèse bibliographique
continuité. Les deux éléments se dérivent du fameux élément non conforme à neuf
paramètres baptisé BCIZ [171] et qui ne vérifie que la continuité C0. Pour assurer une
faible continuité C1 du déplacement w, les auteurs ont adopté le concept de Kirchhoff
discret pour formuler le premier élément nommé PDKT, quant à l'autre élément
triangulaire raffiné PRT9, la méthode proposée par Chen et al. ([178, 179]) est utilisée.
En 2006, Chen et al. [176] ont proposé une théorie généralisée (nommée la
théorie d’ordre supérieur: mème ordre global- 1,2-3 ordre local) pour les plaques
composites multicouches. le développement est basé sur une théorie d’ordre supérieur
et autre locale d’ordre 1,2-3, proposée par Li et Liu. Afin d’estimer son applicabilité
sur des structures composites avec différent nombre de couches, ils ont utilisée un
champ du déplacement composé d’un polynôme d’ordre m ( ) de coordonnée
globale z suivant l'épaisseur et d’une série de puissance d’ordre 1, 2-3 de la
coordonnée locale k au niveau de chaque couche. Ils ont montré que la présente
théorie reste encore efficace quand le nombre de couches est compris entre 5 et 14. Il
ont montré aussi que l’élément triangulaire amélioré de Chen et Cheung ([179, 181])
possède des performances élevées ayant comparé avec les éléments reconnus.
Pour éviter l’utilisation des fonctions du déplacement de continuité C1, le
travail qu’a été fait récemment par Xiaohui, Wanji et Zhen [182], avait pour but le
développement d’une nouvelle théorie de zigzag qui exige au niveau d’implémentation
éléments finis que des fonctions de forme C°. La fonction zigzag proposée satisfait a
priori les conditions de continuité inter-laminaire de contraintes du cisaillement ainsi
que les conditions d'homogénéité. En outre, le nombre total d'inconnues impliquées
dans le modèle est indépendant du nombre de couches. Une comparaison vis à vis
d'autres théories de zigzag a montré que les dérivées premières du déplacement
transverse ont été ressorties à partir du champ de déplacement plan du modèle
proposé, de sorte que des fonctions d'interpolation de type C° uniquement sont
nécessaires à la mise en œuvre « élément fini».
1.6 Conclusion
Toute continuité des contraintes inter laminaire nécessite forcément une
discontinuité des déformations transverses. Cette condition pourrait se traduire par
un champ de déplacement ou moins linéaire par couche, échelonné ou sous forme de
zigzag. Trois différentes approches sont disponibles et qui sont largement employées
dans la formulation des théories ainsi que les modèles numériques correspondants.
Elles se résument aux travaux de Lekhnitskii, Ambartsumian et Reissner, dont
l'objectif est de reproduire une cinématique et un état de contraintes assez proches de
la réalité dans le comportement des structures plaques/coques composites.
Afin de montrer l’effet de la fonction C0 de la forme zigzag de Murakami sur les
théories classiques et d’ordre supérieur, plusieurs auteurs [5, 123, 124, 137], ont conclu
qu'un modèle en zigzag constitue un bon outil pour augmenter les performances des
deux théories. En particulier et à travers des évaluations numériques, ils ont montré
40
Chapitre 1. Théories aux éléments finis des plaques et coques composites.
Synthèse bibliographique
qu’il est plus commode d'améliorer une théorie des plaques en introduisant une
fonction zigzag que la raffiner en ajoutant deux ou trois termes d'ordre supérieur.
Cependant, pour bien représenter les effets locaux ou d'analyser les plaques épaisses,
un modèle avancé (par couche ou 3D) est requis en telles conditions.
L’avantage de la formulation en zigzag en utilisant la fonction de Murakami
est qu’on peut l’exploiter pour les deux composantes du déplacement autant que pour
le troisième. Cependant, il est difficile d'appliquer les théories de Lekhnitskii et de
Ambartsumian si l'effet du zigzag est exigé pour le déplacement transversal où
l’ajout des efforts additionnels est nécessaire [124].
Du point de vue modélisation numérique, des éléments finis généralement
conformes et de haute précision ont été formulés et élaborés par Di Sciuva. Ces
éléments ne requièrent que cinq déplacements généralisés pour décrire la cinématique
de déformation de la plaque stratifiée. Ils satisfont à priori les conditions de
continuité géométrique et des contraintes sur les interfaces. Néanmoins, ils soufrent
d'un défaut fâcheux, qui est l'exigence de la continuité du degré de liberté de la
flèche transverse afin que son interpolation du type Hermiten soit impliquée dans
les modèles éléments finis. Généralement la préséance de dix degrés de liberté par
nœud rend difficile l'implémentation de ces éléments dans des codes de calcul
commerciaux qui n'admettent que six degrés de liberté.
Plusieurs scientifiques ont trouvé que la théorie proposé en 1993 par Cho et
Parmerter pour formuler leurs éléments finis des plaques multicouches/sandwichs est
la plus attractive à adopter parmi les théories raffinées de déformation en
cisaillement d’ordre supérieur (RHSDT). Elle satisfait la majorité des exigences de
continuité réclamées dans une formulation envisagée. Afin d’améliorer les
performances de leurs résultats, le concept des éléments conformes des plaques
minces a été utilisé dans certains cas.
Un examen sur les théories de zigzag existantes, nous a permis de conclure que
selon l'intégration des équations d'équilibre tridimensionnel, des dérivées du
déplacement d'ordre élevé devraient être calculées (généralement la continuité C1 des
fonctions du déplacement transverse est nécessaire) pour prédire les contraintes de
cisaillement transverses avec précision.
Afin de surmonter ces exigences, nous évoquons particulièrement ici, un
modèle élément fini tridimensionnel basé sur la théorie dite de sous-couche
« sublaminate theory», proposée par Cho et ses collaborateurs. Ce modèle combine à
la fois les avantages de la théorie par-couche et celle de zigzag d’ordre supérieur. Le
fait qu’il n’inclut que les degrés de liberté standards, le rend compatible avec les
exigences recommandées dans les codes éléments finis commerciaux.
Vu les performances de quelques théories et éléments précités, qui se sont
inspirées par fois d’elles mêmes, et les conditions de continuités exigées dans les
structures composites, les scientifiques ont toujours le souci de préserver les
approches classiques ou de premier ordre. Toutefois, il est évident que le modèle en
41
Chapitre 1. Théories aux éléments finis des plaques et coques composites.
Synthèse bibliographique
déformation du premier ordre seul est incapable de présenter d’une façon fiable les
effets locaux dans les plaques composites. Donc, La nécessité d’un modèle qui inclut la
forme zigzag autant que des termes d’ordre élevé et qui respecte les conditions de
continuité, rend les théories dites « order global-local theories » plus compétitives
et délimitent toute formulation future envisagée.
42
Chapitre 2 Modèle Discret en Déplacements avec une
fonction Zigzag pour les plaques composites.
Formulation théorique
2.1 Introduction
Dans le début des années 1980, de nombreux travaux ont été consacrés à la
formulation dite « Discret Kirchhoff », dont plusieurs éléments finis ont été développés
avec succès par la suite, notamment ceux basés sur le concept de Batoz et BenTahar
[183]. Malgré leur puissance, ces éléments étaient, malheureusement incapables de
représenter un comportement adéquat de plaques épaisses. La théorie de Reissner-
Mindlin qui tient en compte l'effet de cisaillement transverse s'est avéré une
alternative à cet effet, dont plusieurs auteurs préservaient ses hypothèses dans leurs
développements, en citant les éléments simples et robustes pour les plaques isotropes,
issus de la même école, formulés par Lardeur [63], Katili [1, 2] et Ayad [57]. Bien que
cette formulation conduise souvent à un champ de contrainte de cisaillement constant
et fait recours aux paramètres de correction, elle reste convaincante comme support
pour toute future formulation améliorée pour les plaques composites.
Ces dernières décennies, plusieurs travaux de recherche se sont focalisés sur la
formulation des modèles éléments finis qui fournissent des prévisions précises sur
l'état de contraintes. L'adoption de l'approche de Mindlin constitue l'un de nos axes
d'intérêt, où un modèle variationnel Discret en Déplacement de Mindlin, labélisé
DDM « Displacement Discrete Mindlin » est développé [3]. Celui ci prend en compte
l'effet de cisaillement transverse intégrant le caractère multicouche des structures
composites. Outre que les éléments développés passaient les patch-tests, ils sont à
l'abri de tous problèmes numériques (verrouillage, distorsion ... etc.) [4]. Afin
d'appliquer ce concept pour les structures composites, en utilisant des éléments finis
simples ayant que des nœuds aux sommets, deux hypothèses de Mindlin modifiées
localement ont été introduites au modèle isotrope. La première hypothèse est
cinématique, elle consiste à introduire, sous forme d’une intégrale de contour, une
équation cinématique de déformation du cisaillement transverse. Elle permet
l'élimination de verrouillage en cisaillement, sans introduire des fonctions bulles, ni
faire le recours à l'intégration réduite. Telle méthode est communément connue
« ANS » (Assumed Naturel Strain), développée par Bathe et al. [62]. La seconde
hypothèse est mécanique. Elle intègre l'aspect multicouche de la plaque à travers
deux lois de comportement, en flexion et en cisaillement, et deux équations
43
Chapitre 2. Modèle Discret en Déplacements avec une fonction Zigzag pour les
plaques composites. Formulation théorique
44
Chapitre 2. Modèle Discret en Déplacements avec une fonction Zigzag pour les
plaques composites. Formulation théorique
N
i P
k
u1M et u2M sont les variables généralisées associées à la fonction zigzag M(z) (relation
1.18a,b).
45
Chapitre 2. Modèle Discret en Déplacements avec une fonction Zigzag pour les
plaques composites. Formulation théorique
βx 4 βxi _ 8 Ck
Ni Pk Δβsk (2.2c)
β y i 1 β yi k 5 Sk
, sont les rotations nodales aux sommets du quadrilatère, Ni et Pk sont des
fonctions d’interpolations prescrites en termes des fonctions indépendantes linéaires
et quadratiques respectivement (tableau 2.1). Les coefficients et sont les cosinus
directeurs du vecteur reliant le nœud i par j du côté k, faisant un angle (figure 2.2).
Les rotations tangentielles sont liées aux rotations cartésiennes par la relation
suivante:
(2.3)
(2.4)
(2.5)
46
Chapitre 2. Modèle Discret en Déplacements avec une fonction Zigzag pour les
plaques composites. Formulation théorique
Elimination des Δβ sk 4
γ ξz 8 s 4 γs7
γ s8 γ ξz 7 3
1 7 1
8 s 3 Ae
s 6
5
2
2 s Z, w
γ ξz 5 γ s 5 γ s 6 γ ξz 6 βx β
M M y
wi , βxi , βyi , u , u1 2 u2M Y
sk X M
u 1
(2.6)
Il se résume à:
; (2.7)
et sont respectivement les déformations de flexion et du cisaillement
transverse couplés avec des déformations dues au zigzag.
47
Chapitre 2. Modèle Discret en Déplacements avec une fonction Zigzag pour les
plaques composites. Formulation théorique
Soit:
(2.8)
Où, est une courbure de flexion, est une déformation due au zigzag. Elles
sont définies par:
(2.9)
(2.10)
(2.11)
avec
0 N i , x 0 0 0
B1bβ 0 0 Ni , y 0 0 i 1,4 (2.12)
0 Ni , y Ni , x 0 0
et
0 0 0 N i , x 0 ...
B2bβ 0 0 0 0 N i , y ... i 1,4
0 0 0 N i , y N i , x ...
(2.13)
d'où
Ni , x j11Ni , ξ j12 Ni , η
(2.14)
Ni , y j21Ni , ξ j22 Ni , η
et
C k Pk ,x
BbΔβ S k Pk , y k 5,6,7,8 (2.16)
C k Pk , y S k Pk ,x
48
Chapitre 2. Modèle Discret en Déplacements avec une fonction Zigzag pour les
plaques composites. Formulation théorique
avec
Pk , x j11Pk , ξ j12Pk , η
Pk , y j21Pk , ξ j22Pk , η (2.17)
et
Δβ n Δβ s5 Δβ s6 Δβ s7 Δβ s8 (2.19)
γ xz
γ
s γ 0
2
( 1) l u M (2.20)
γ
yz h k
avec
u u M
1
M
M
(2.21)
u 2
49
Chapitre 2. Modèle Discret en Déplacements avec une fonction Zigzag pour les
plaques composites. Formulation théorique
où
γξz
γ
ξ
γηz (2.24)
est l’inverse de la matrice Jacobéenne (rel. 2.15).
γξz 12 (1 η)γξς
l
5 2 (1 η)γξς7
1 l
1 (2.25)
γηz 2 (1 ξ )γης 8 2 (1 ξ )γης 6
l 1 l
L5 L7
γξς5 γs5 , γξς7 γ s7
2 2 (2.26)
L6 L8
γης 6 γs6 , γης 8 γs 8
2 2
Posant
L5 J11(1 η) L6 J12 (1 ξ ) L7 J11(1 η) L8 J12 (1 ξ )
N γ 1 (2.27)
4 L5 J 21(1 η) L6 J 22 (1 ξ ) L7 J 21(1 η) L8 J 22 (1 ξ )
50
Chapitre 2. Modèle Discret en Déplacements avec une fonction Zigzag pour les
plaques composites. Formulation théorique
avec
(2.29)
M x ,x M xy , y Tx 0
M xy ,x M y , y Ty 0
(2.30)
H x ,x H xy , y H xz 0
H xy ,x H y , y H yz 0
Moments de flexion
Mx h2 σx
M M y z σ y dz H f χ1 H fz χ 2 (2.31)
M xy h 2 τ xy
51
Chapitre 2. Modèle Discret en Déplacements avec une fonction Zigzag pour les
plaques composites. Formulation théorique
Hx h 2 σx
H P H y M ( z ) σ y dz H fz χ1 H zz χ 2
(2.32)
H xy h 2 τ xy
H xz h 2 2( 1)l τ xz
HT dz 2H c 01 γ0 H c11 u
M
(2.34)
H yz h 2 hl τ
yz
avec
h2
H z H ( z )dz 1 ( z
nc
f
2 3
l 1 z l3 )H
h 2 3 l 1
h2
H H z M ( z ) H ( z )dz 2( 1)
l
nc
1 3 1 2
( z l 1 z l ) z 0l ( z l 1 z l ) H (2.36a-c)
3 2
fz zf
h 2 h l 1 l 3 2
h2 l nc
4 1 3
H M
zz
2
( z ) H ( z )dz ( z l 1 z l ) z 0l ( z l 1 z l ) z 0l .hl H
3 2 2 2
h 2 l 1 hl 3
h2 nc
H H
c τ ( z )dz hl H τ
h 2 l 1
h2
( 1) l nc
H
c 01 H τ ( z )dz 2.( 1) l H τ (2.37a-c)
h 2 hl l 1
h2 nc
4
H M
c11
2
( z )H τ ( z )dz H τ
h 2 l 1 hl
52
Chapitre 2. Modèle Discret en Déplacements avec une fonction Zigzag pour les
plaques composites. Formulation théorique
Avec
~ 2( 1) M l
γsk (w,s βs us )k
hl (2.39)
L'approximation de est déjà définie dans l'expression (2.4), cependant, on propose une
fonction linéaire suivant pour le terme de zigzag , cette variation est semblable à celle
utilisée pour définir , (l'équation 2.5).
w j wi 1 2 ( 1)l M
γsk ( βsi βsj ) Δβsk (usi usjM )
Lk 2 3 hl (2.40)
soit:
~ 2( 1)l M 2( 1)l M
γnk sk w,x βx u1 ck w,y β y u2 (2.41)
hl hl
2( 1) l
γ nk s k .w, x c k .w, y s k . β x c k . β y
~ s k .u1M c k .u2M (2.42)
hl
2( 1) M l
γ nk s k .w, x c k .w, y β nk
~ un (2.43)
hl
Posant
53
Chapitre 2. Modèle Discret en Déplacements avec une fonction Zigzag pour les
plaques composites. Formulation théorique
Avec :
H P H
τk k τ ( z )Pk
1
(2.48)
cosθk sin θk
P P
k k
1
(2.49)
sin θk cosθk
est l'angle comprise entre l'axe du repère cartésien et l'axe s reliant le nœud i et
j du bord k (figure 2.2).
en substituant les déformations de cisaillement (rel. 2.46) dans (2.47), les efforts
tranchants définis dans le repère tangentiel au coté k seront:
h
54
Chapitre 2. Modèle Discret en Déplacements avec une fonction Zigzag pour les
plaques composites. Formulation théorique
d'où:
γ s0
1 Ts
nc usiM usjM
H ck ( 1) H τk M
l
M
γ
k
n 0 T
k
n l 1 u
ni unj k
(2.52)
La substitution de la déformation s 0 , équation (2.52), dans l'expression (2.46) conduit
à l'expression suivante:
Nous procédons à écrire dans le repère tangentiel les relations d’équilibre (§ 2.2.3)
sur la surface moyenne reliant les efforts tranchants en fonction des moments,
(figure 2.6).
soit:
M Ms Mn M sn (2.54)
avec
(2.56)
, (2.57)
(2.58)
avec
(2.59)
55
Chapitre 2. Modèle Discret en Déplacements avec une fonction Zigzag pour les
plaques composites. Formulation théorique
nc usiM usjM
γ sk φkml Δβsk H ck
inv
11
inv
H ck12
( 1) l
H τk M M
l 1 uni unj (2.60)
( 1) M
l
usi usjM
hl
avec φkml est le facteur d’influence du CT, défini par :
H ckinv11 H ckinv12
φkml 8( H f 11 H f 13 ) (2.61)
L2k L2k
la combinaison des équations (2.40) et (2.60) permet de définir les incréments des
rotations, il en résulte :
wi w j 1
Lk 2 ( βsi βsj )
3
Δβsk (2.62)
2 3φkml inv nc u M
u M
H ck11 H ck12 ( 1) H τk M
inv l si sj
M
l 1 uni unj
Remarques
1. Les rotations des mi- cotés sont définies sur le plan moyen de la plaque.
L'expression (2.62) montre leur dépendance du caractère globale de plaque.
2. Les termes de déplacement transverse et associés à l'écriture des variables
, équation (2.62), sont explicitement introduits lors de l'évaluation des courbures
.
56
Chapitre 2. Modèle Discret en Déplacements avec une fonction Zigzag pour les
plaques composites. Formulation théorique
avec
3
2 3φ5ml 0 0 0
3
0 0 0
2 3φ6ml
An 3 (2.64)
0 0 0
2 3φ7ml
3
0 0 0 ml
2 3φ8
et
2 L5 C5 S5 0 0 2 L5 C5 S5 0 0
0 0 0 0 0 2 L6 C6 S6 0 0
0 0 0 0 0 0 0 0 0 0
2 L8 C8 S8 0 0 0 0 0 0 0
AG 0 1 (2.65)
2 0 0 0 0 0 0 0 0 0 0
2 L6 C6
C6 0 0 0 0 0 0 0
2 L7 C7 S7 0 0 2 L7 C7 S7 0 0
0 0 0 0 0 2 L8 C8 S8 0 0
(2.66)
Avec
F1k Hck11 Hτk11 .Ck Hτk12 .Sk Hck12 Hτk 21 .Ck Hτk 22 .Sk ,
inv inv
k 5,6,7,8 (2.67)
F2k Hck11 Hτk11 .Sk Hτk12 .Ck Hck12 Hτk 21 .Sk Hτk 22 .Ck ,
inv inv
k 5,6,7,8 (2.68)
et
nc
H τk ( 1)l H τk , k 5,6,7,8 (2.69)
l 1
57
Chapitre 2. Modèle Discret en Déplacements avec une fonction Zigzag pour les
plaques composites. Formulation théorique
s γ φ n G0 G1 G1 G2 n (2.70)
hl
avec
φ5ml 0 0 0
φ ml
Aφ 0 6 0 0
(2.71)
0 0 φ7ml 0
0 0 0 φ8ml
et
0 0 0 C5 S5 0 0 0 C5 S5 0 0 0 0 0 0 0 0 0 0
1 0 0 0 0 0 0 0 0 C6 S6 0 0 0 C6 S6 0 0 0 0 0
AG 2
2 0 0 0 0 0 0 0 0 0 0 0 0 0 C7 S7 0 0 0 S7 S7
0 0 0 C8 S8 0 0 0 0 0 0 0 0 0 0 0 0 0 C8 S8
(2.72)
L'appariation du terme en fonction de l et les variables additionnels de zigzag
(équation 2.70), montre que la déformation en cisaillement est liée au couches, Telle
forme est différente à celle du modèle habituel du premier ordre par une allure
constante par couche et qui reflète une discontinuité entre les couches.
e
Wint est le travail virtuel interne élémentaire, Wexte est le travail élémentaire des
forces extérieures.
soit:
Wint
e
( εs* σ s γs* τ s ) dV (2.74)
Ve
k k k
e
fz cz (2.76)
58
Chapitre 2. Modèle Discret en Déplacements avec une fonction Zigzag pour les
plaques composites. Formulation théorique
D’où
B f T H f B f B f T H fz Bz
k fz dA (2.77)
e Bz H fz B f Bz H zz Bz
T T
A
Avec
B f B1bβ B1bΔβ An AG 0 AG1 (2.78)
Bz B2bβ
H , H , H (rel. (2.36a-c))
f fz zz
H , H
c c 01 , H (rel. (2.37a-c))
c11
Il en résulte:
59
Chapitre 2. Modèle Discret en Déplacements avec une fonction Zigzag pour les
plaques composites. Formulation théorique
χ A H A M
1 1 P 2
(2.82)
χ B H B M
2 1 P 2
avec
A A H
1
1
fz
1
B H H A
1 fz
1
f 1 (2.83)
B H I H A
2 fz
1
f 2
A I H H H H fz
1
zz fz
1
f
1
Afin d'éviter toute éventuelle singularité de ces matrices, le schéma de calcul tiendrait
en compte l'absence de l'effet de zigzag dans la version simplifiée du modèle DMQP.
En se basant sur les équations d'équilibre en contraintes, les contraintes de
cisaillement s'écrivent:
alors
σ zA (z)H A (z)M M (z)B (z)H B (z)M
1 P 2 1 P 2 (2.86)
avec
1
2
dz , D2
2 A131 A123
dz
t 2 A231 A223 A222 A233 t A122 A133
1
2
dz , E 2 2 B131 B123
dz
t 2 B2 B2
31 23
B222 B233 t B122 B133
60
Chapitre 2. Modèle Discret en Déplacements avec une fonction Zigzag pour les
plaques composites. Formulation théorique
et
T T T x y M x ,x M xy , y M xy ,x M y , y
H H
T Tx H Ty H x ,x H xy , y H xy ,x H y , y
(2.89a-d)
λ M M x ,x M xy , y M xy ,x M y , y M y ,x M x ,y
λH H x ,x H xy , y H xy ,x H y , y H y ,x H x ,y
(2.90)
on aura donc
(2.91)
(2.92)
61
Chapitre 2. Modèle Discret en Déplacements avec une fonction Zigzag pour les
plaques composites. Formulation théorique
(2.93)
Rappelons qu’une procédure itérative est nécessaire pour un calcul des valeurs et les
modes propres, en utilisant la méthode inverse pour le calcul de la petite valeur et
l’orthogonalisation de Gram Shmidt pour les modes supérieurs ( ).
(le code Reflex 1.0, [188, 190]).
Sachant que
(2.95)
2.5 Conclusion
En se basant sur des hypothèses modifiées de Mindlin discrète, un élément fini de
plaque multicouche à quatre nœuds, travaillant en flexion/cisaillement est formulé.
L'élément possède en outre de trois degrés de liberté traditionnelles, à savoir un
déplacement transverse et deux rotations, deux variables additionnelles. La démarche
utilisée est celle de Reissner-Murakami-Carrera, basée sur l'ajout d'une fonction de
zigzag. Les termes additionnels ayant un sens de déplacement permettent d'impliquer
un déplacement linéaire par morceaux dans une seconde version du modèle.
Ce choix de telle fonction du zigzag se justifie, selon Carrera [124], du point de vue
implémentation par : l’incorporation de la fonction de zigzag de Murakami dans les
modèles existants de plaque/coques/poutres requiert le même effort que celui
nécessaire à inclure un terme additionnel d’ordre élevé. D’autre coté, du point de vue
numérique, l’incorporation de cette fonction conduit à des améliorations significatives
qui ne peuvent être obtenues par d’autres fonctions.
62
Chapitre 3 Validation numérique du modèle de plaques
composites DMQPz
3.1 Introduction
La thématique du développement des éléments finis d'ordre supérieur se base
principalement sur des termes secondaires attribués à la partie globale du champ de
déplacement. La fiabilité du modèle standard qu'ils en dérivent est généralement
garantie dans le domaine isotrope. Selon l'approche utilisée et d'autant que l'ordre de
l'interpolation est élevé, la précision ne diffère que légèrement d'un élément à un
autre. Cependant les termes secondaires qui sont des valeurs ajoutées servent à
améliorer la réponse locale à travers l'épaisseur de la structure (contraintes et
déformations dans le sens transversal). Par conséquent, afin d'analyser le modèle
proposé, nous mettons en relief le modèle en zigzag avec une validation qui se focalise
surtout sur les plaques stratifiées/sandwichs. On évoquera brièvement l'énoncé de
quelques tests standards et patch test concernant les structures isotropes, élaborés
antérieurement [4]:
Le premier test, est un patch-test mécanique, proposé par Irons [191] et
considéré comme une condition suffisante de convergence. Ce test a permet de vérifier
l’aptitude du modèle DMQP avec un maillage grossier à reproduire un état de
contrainte constant, en imposant un nombre minimum de ddl nécessaire à
l'élimination des mouvements de corps rigide.
D’autre part, afin d’étudier l’influence de la distorsion géométrique sur la
précision du présent modèle, le problème de la plaque biaise (angle aigu= 30°) a été
traité. La plaque est simplement supportée sur tout le contour et soumise à une
charge uniforme. La convergence de la flèche et des moments maximum et minimum
au centre des plaques sont confrontés avec des solutions de plaques minces calculées
par Morley [192]. Une solution tridimensionnelle obtenue par Babuska et al. [193]
ainsi que des résultats de l’élément fini MITC4 [62] ont été également utilisés. Dans le
cas très mince (L/h = 1000), la convergence vers la solution de plaque mince [192] est
obtenue par l’ensemble des éléments avec des vitesses différentes. Cependant, dans le
cas limite (L/h=100) nous avons constaté une convergence vers une solution autre que
celle de Morley pour les plaques minces. Leur solution 3D est très proche de celle
obtenue par l’élément DMQP.
Le test de la plaque circulaire isotrope épaisse et mince (R/h=2.5 et 50) est un
test encore bien vu car il met en évidence la distorsion géométrique des éléments.
Deux conditions aux limites ont été examinées, l’évaluation de la flèche et le moment
63
Chapitre 3. Validation numérique du modèle de plaques composites DMQPz
64
Chapitre 3. Validation numérique du modèle de plaques composites DMQPz
; sur CD sur AD
Rappelons que certains aspects des solutions exactes donnés par Pagano et al. [194]
sont basés sur la formulation développée par le premier auteur dans la référence [50],
ainsi que ceux analogues de la théorie classique des plaques. Tous les résultats sont
normalisés avec les valeurs suivantes :
, ,
Avec ,
65
Chapitre 3. Validation numérique du modèle de plaques composites DMQPz
DMQPz
5 DMQP
E.O
DSQ
MiSP4/ml
Elasticité
4 CPT
w
0
0 5 10 15 20 25 30 35 40 45 50
a/h
66
Chapitre 3. Validation numérique du modèle de plaques composites DMQPz
67
Chapitre 3. Validation numérique du modèle de plaques composites DMQPz
0.5
DMQP
0.4 DMQPz
0.3
0.2
0.1
0.05 0.1 0.15 0.2 0.25 0.3
0
-0.1 yz (a/2, 0,z)
-0.2
-0.3
-0.4
-0.5
68
Chapitre 3. Validation numérique du modèle de plaques composites DMQPz
-0.2
-0.3
-0.4
-0.5
0.5
DMQP
0.4 DMQPz
0.3
0.2
0.1
0.05 0.1 0.15 0.2 0.25 0.3
0
-0.1 yz (a/2,0,z)
-0.2
-0.3
-0.4
-0.5
69
Chapitre 3. Validation numérique du modèle de plaques composites DMQPz
DMQP
0.5
DMQPz
0.4
0.3
0.2
0.1
0.05 0.1 0.15 0.2 0.25
0
-0.1 xz (0,a/2,z)
-0.2
-0.3
-0.4
-0.5
70
Chapitre 3. Validation numérique du modèle de plaques composites DMQPz
C Modèles Erreur. %
1 DMQPz 180.48 0.3148
DMQPml 181.23 0.6352
DSQ 181.42 0.2044
Heterosis 183.99 1.6239
MiSP4/ml 182 0.5247
Elasticité 181.05
10 DMQPz 41.95 0.0954
DMQPml 42.00 0.5727
DSQ 42.04 0.3102
Heterosis 41.92 0.0239
MiSP4/ml 43.845 4.6170
Elasticité 41.91
50 DMQPz 16.82 0.4179
DMQPml 16.87 1.0746
DSQ 16.89 0.8358
Heterosis 16.85 0.5970
MiSP4/ml 18.8 12.239
Elasticité 16.75
Tableau 3.5 Plaque carrée sandwich sous chargement uniforme, déplacement
71
Chapitre 3. Validation numérique du modèle de plaques composites DMQPz
. L’analyse est faite pour un seul rapport , et un maillage de 6×6 a été utilisé
pour le quart de la plaque. Les résultats sont donnés sous la forme adimentionnelle
suivante :
D’après les résultats de la flèche au centre de plaque (tableau 3.5), les deux
versions du modèle DMQP sont en bon accord avec la solution d’élasticité. Un examen
minutieux des erreurs relatives, montre que le modèle avec zigzag apporte une valeur
ajoutée au modèle DMQP. D’ailleurs cette remarque est également confirmée au
niveau des contraintes planes, calculées au centre de la plaque. Le modèle avec zigzag
dans ce test montre une excellente précision qui se voie clairement pour le cas C=50.
En revanche le modèle DMQP ressemble aux autres éléments du premier ordre, qui
soufrent généralement d’inexactitude au niveau des contraintes planes dans le cas des
structures fortement anisotropes.
72
Chapitre 3. Validation numérique du modèle de plaques composites DMQPz
73
Chapitre 3. Validation numérique du modèle de plaques composites DMQPz
a/h Modèles
DMQP/ml .392 .613 .0439 .305 .218 2.238
DMQPz .687 .524 .0480 .247 .206 1.972
E. O. .391 .572 .0448 .308 .251 -
4 HSDT [116] .734 .503 .050 .286 .209 1.922
RFSDT [116] .620 .661 .045 .262 .254 1.891
Elasticité[50] .755 .556 .0505 .256 .217 -
.801 .534 .0511
74
Chapitre 3. Validation numérique du modèle de plaques composites DMQPz
,
Les valeurs du déplacement transversal, les contraintes planes et de
cisaillement transverse, pour les deux formes de la plaque sont présentées dans les
tableaux 3.7 à 3.9. Pour une plaque modérément épaisse a/h = 4, la précision de la
contrainte plane obtenue par différents éléments du premier ordre (DMQP/ml,
DSQ) apparaît affectée pour ce rapport d’élancement. Les valeurs obtenues par
l’élément DMQPz sont plus proches aux celles de référence. En revanche, les deux
modèles DDM montrent une bonne performance concernant la flèche au centre et les
autres contraintes planes ( et ), particulièrement pour les rapports a/h = 50 et
100 (tableau 3.8). Généralement, pour les petits rapports, l'élément DMQPz se
comporte mieux que le modèle DMQP, dont la performance reste comparable à celle
donnée par le modèle raffiné d’ordre supérieur. En comparaison avec l'approche par
morceaux, les résultats obtenus par le présent modèle en zigzag, pour la structure
rectangulaire (b/a=3) (tableau 3.8), sont encore similaires à ceux de la théorie de
zigzag raffinée; basée sur la théorie de déformation en cisaillement du premier ordre
(FSDTZZ), développé par (Fares et al. 2008). Les valeurs de contraintes de
cisaillement transverses demeurent, dans l’ensemble satisfaisantes, même mieux que
celles obtenues par DSQ, EO et HSDT.
Pour la plaque carrée multicouche modérément épaisse ( ), l'allure du
déplacement au point B est tracée sur la figure 3.9. Les courbes illustrent les
valeurs rapportées aux deux versions du modèle DMQP confrontées avec les valeurs
75
Chapitre 3. Validation numérique du modèle de plaques composites DMQPz
exactes [50] et celles de la référence [20]. Un déplacement linéaire par morceaux sous
forme de zigzag est représenté avec un changement du signe de la pente entre les
couches.
a/h Modèles
76
Chapitre 3. Validation numérique du modèle de plaques composites DMQPz
0.5
0.3
0.1 Exact
DMP Qz
z/h
DMQP
-0.1 FSDT
-0.3
-0.5
-0.012 -0.008 -0.004 0 0.004 0.008 0.012
u(a/2,0,z)
0.5
0.3
Exact
0.1 DMQP z
z/h
DMQP
FSDT
-0.1
-0.3
-0.5
-1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1
σxx(a/2, a/2,z)
77
Chapitre 3. Validation numérique du modèle de plaques composites DMQPz
0.5
0.3
Exact
DMQP z
0.1 DMQP
FSDT
z/h
-0.1
-0.3
-0.5
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4
σxz(0, a/2,z)
0.5
Exact
DMQP z
0.3 DMQP
FSDT
0.1
z/h
-0.1
-0.3
-0.5
0 0.05 0.1 0.15 0.2 0.25
σyz(a/2,0,z)
Les valeurs relatives au modèle DMQPz ne sont obtenues que par simple
addition du terme linéaire standard dû à la rotation et le terme de zigzag qui
78
Chapitre 3. Validation numérique du modèle de plaques composites DMQPz
79
Chapitre 3. Validation numérique du modèle de plaques composites DMQPz
a/h=10 a/h=4
80
Chapitre 3. Validation numérique du modèle de plaques composites DMQPz
Cas 2
Peaux :
Cœur :
hf
h/2
hc
h/2
hf
La flèche au centre de la plaque (a/2,a/2) avec la solution analytique ([200], [50], [201]
,[202]) et les résultats de l'élément fini de Khatua et al. [203] sont présentés sur le
tableau (3.10). Les valeurs obtenues sont également comparées avec les modèles du
premier et d’ordre supérieur développés par Topdar et al. [134]: PFSDT, PHSDT et
PRHSDT (§ 1.4.2.1). Les résultats afférents aux modèles DMQP et DMQPz sont
presque identiques et possèdent un excellent agreement avec la solution analytique,
voire meilleur, comparés à ceux obtenus par Khatua et Topdar dans le deuxième cas.
Par contre, les valeurs obtenues par les deux modèles PFSDT et PHSDT sont
inférieures aux celles de référence et présente une erreur flagrante.
81
Chapitre 3. Validation numérique du modèle de plaques composites DMQPz
Matériau Modèles
82
Chapitre 3. Validation numérique du modèle de plaques composites DMQPz
0.5
0.4 DMQPz
DMQP
0.3
0.2
0.1
-5 5 10 15 20 25 yz (a/2,0,z)
z/h
-0.1
-0.2
-0.3
-0.4
-0.5
0.5
0.4
DMQP
0.3 DMQPz
0.2
0.1
20 40 60 80 100 120 140 160 180 200
z/h
yz (a/2,0,z)
-0.1
-0.2
-0.3
-0.4
-0.5
83
Chapitre 3. Validation numérique du modèle de plaques composites DMQPz
84
Chapitre 3. Validation numérique du modèle de plaques composites DMQPz
Tableau 3.13. Contraintes planes dans la plaque carré sandwich (f/c/f), simplement
supportée soumise à une charge sinusoïdale.
85
Chapitre 3. Validation numérique du modèle de plaques composites DMQPz
Peau
25.× 1.× .25
(GPa)
Cœur
.25
(GPa)
Les valeurs obtenues à certains points sont calculés sous forme adimensionnelle
suivante:
,
,
86
Chapitre 3. Validation numérique du modèle de plaques composites DMQPz
Modèles Maillage
DMQPz 6×6 1.4318 1.6759 0.356 0.3494
16×16 1.4304 1.6722 0.3606 0.3522
DMQPml 8×8 1.4544 1.68 0.362 0.386
50 19×19 1.4512 1.664 0.368 0.358
PRHSDT[204] 8×8 1.4505 1.6723 0.3625 0.3460
12×12 1.4504 1.6704 0.373 0.3575
Exacte[50] 1.4507 1.6702 0.3735 0.3578
DMQPz 7×7 1.5804 1.6539 0.3541 0.3444
16×16 1.5789 1.6503 0.3581 0.3484
DMQPml 8×8 1.7162 1.655 0.359 0.3505
20 19×19 1.71 1.6425 0.3635 0.3545
PRHSDT[204] 8×8 1.7063 1.6693 0.3613 0.3456
12×12 1.7063 1.6678 0.3718 0.3571
Exacte[50] 1.6998 1.6682 0.3701 0.3551
DMQPz 7×7 2.1106 1.6056 0.3432 0.3462
16×16 2.1085 1.6038 0.3473 0.3507
DMQPml 8×8 2.652 1.61 0.35 0.346
10 19×19 2.636 1.6 0.355 0.353
PRHSDT[204] 8×8 2.6228 1.6855 0.3609 0.3469
12×12 2.6228 1.6854 0.3713 0.3582
Exacte[50] 2.5964 1.6914 0.3629 0.3490
Tableau 3.14 Plaque carrée sandwich ( f(0,90), c, f(90,0) ) sous chargement uniforme.
déplacement transversal au point C, contrainte plane et contraintes de CT maximales.
0.5 DMQP
DMQPz
0.4
0.3 S=50
0.2
0.1
0.05 0.1 0.15 0.2 0.25 0.3 0.35
xz (0,a/2,z)
-0.1
-0.2
-0.3
-0.4
-0.5
87
Chapitre 3. Validation numérique du modèle de plaques composites DMQPz
88
Chapitre 3. Validation numérique du modèle de plaques composites DMQPz
Nous avons pris en suite le cas pour analyser les quatre premières
fréquences propres. Les résultats sont reportés dans le tableau (3.16) ( et
représentent le nombre des bosses et les creux d’ondes suivant x et y respectivement).
D’après les valeurs obtenues, le modèle DMQP montre un bon accord avec les
solutions de références. Cette performance est constatée pour tous les cas. Quant au
modèle avec zigzag, quoiqu’il s’avère un peu rigide, il converge lentement mais mieux
que le modèle mixte hybride MiSP4/ml. Cela est remarqué notamment pour le
troisième et le quatrième mode.
6
fréquence 1
Fréquence propre
4 DMQP
DMQPz
DSQ
CLT
Elasticité
0
1 5 9 13 17
a/h
89
Chapitre 3. Validation numérique du modèle de plaques composites DMQPz
Nous limiterons notre étude au cas de la plaque carrée, avec un rapport d’élancement
.
La fréquence propre fondamentale est donnée par la formule suivante :
90
Chapitre 3. Validation numérique du modèle de plaques composites DMQPz
Nous avons constaté, d’après les résultats numériques que le modèle DDM
donne dans les deux cas de matériaux des valeurs de fréquence de bonne précision.
Même pour des valeurs faibles du module élastique de noyau, la version avec zigzag
maintient un bon accord avec la théorie exacte et demeure meilleure que le modèle
mixte (MiSP4/ml).
Géométrie (mm) :
Propriétés
homogénéisées :
91
Chapitre 3. Validation numérique du modèle de plaques composites DMQPz
92
Chapitre 3. Validation numérique du modèle de plaques composites DMQPz
mauvaise évaluation des fréquences propres, et d’autre part, nous rappelons que le
comportement du carton ondulé est sensiblement lié aux variations de la température
et de l’humidité ambiantes.
93
Chapitre 3. Validation numérique du modèle de plaques composites DMQPz
Mode 2
0.515+
0.412 to 0.515
0.309 to 0.412
0.206 to 0.309
0.103 to 0.206
0 to 0.103
-0.103 to 0
-0.206 to -0.103
-0.309 to -0.206
Mode 1 -0.412 to -0.309
0.4015+ -0.515 to -0.412
0.365 to 0.4015 -0.618 to -0.515
0.3285 to 0.365
0.292 to 0.3285
0.2555 to 0.292
0.219 to 0.2555
0.1825 to 0.219
0.146 to 0.1825
0.1095 to 0.146
0.073 to 0.1095
0.0365 to 0.073
0 to 0.0365
Mode 4
0.568333+
Mode 3 0.454667 to 0.568333
0.341 to 0.454667
0.29075+
0.227333 to 0.341
0.2225 to 0.29075
0.113667 to 0.227333
0.15425 to 0.2225
0 to 0.113667
0.086 to 0.15425
-0.113667 to 0
0.01775 to 0.086
-0.227333 to -0.113667
-0.0505 to 0.01775
-0.341 to -0.227333
-0.11875 to -0.0505 -0.454667 to -0.341
-0.187 to -0.11875 -0.568333 to -0.454667
-0.25525 to -0.187 -0.682 to -0.568333
-0.3235 to -0.25525
-0.39175 to -0.3235
-0.46 to -0.39175
Mode 6
1.093E-3+
9.455E-4 to 1.093E-3
Mode 5 7.982E-4 to 9.455E-4
0.344083+
0.252167 to 0.344083
6.510E-4 to 7.982E-4
0.16025 to 0.252167 5.038E-4 to 6.510E-4
0.0683333 to 0.16025 3.565E-4 to 5.038E-4
-0.0235833 to 0.0683333 2.092E-4 to 3.565E-4
-0.1155 to -0.0235833 6.200E-5 to 2.092E-4
-0.207417 to -0.1155 -8.525E-5 to 6.200E-5
-0.299333 to -0.207417 -2.325E-4 to -8.525E-5
-0.39125 to -0.299333 -3.798E-4 to -2.325E-4
-0.483167 to -0.39125 -5.270E-4 to -3.798E-4
-0.575083 to -0.483167
-0.667 to -0.575083
Figure 3.21 Vibrations libres du carton ondulé ( ), les six premiers modes
propres.
94
Chapitre 3. Validation numérique du modèle de plaques composites DMQPz
Fréquences
DMQPz QSR8-ABAQUS
Propres (rad/s)
444.83 444.83
1280.82 1265.8
2772.16 2746.10
4452.42 4370.70
5145.97 5083.40
5564.80009 5552.2
500
450
DMQP
DMQPz
400 MiSP/ml
DMTS
Expérience
350
Pulsation propre (rd/s)
300
250
200
150
100
50
0
90 110 130 150 170 190 210 230 250 270
Longueur L(mm)
95
Chapitre 3. Validation numérique du modèle de plaques composites DMQPz
avec
96
Chapitre 3. Validation numérique du modèle de plaques composites DMQPz
valeurs plus proches au modèle 3D. Le deuxième modèle avec zigzag DMQPz
ressemble à celui de ZIGT FE dans le cas où la structure est modérément épaisse, or
une sur-rigidité se révélait pour un rapport d'élancement au delà de 20.
L’examen graphique de convergence des fréquences naturelles ( et ) pour
deux rapports d'aspect, et montre que les valeurs de références sont un
peu distinctes. Généralement les résultats de la théorie du troisième ordre dite TOT
de Reddy constitue un compromis entre les différentes valeurs de référence. Cette
remarque est bien illustrée sur les figues 3.23. à 3.26. Pour le mode dominant, tous les
résultats montrent une convergence rapide est stable autour de la valeur obtenue par
le premier maillage ( ). Quelque soit le rapport , un maillage raffiné est
nécessaire pour le calcul des modes supérieur. Comme nous l'avons déjà mentionné,
les résultats du modèle DMQP sont comparables à ceux obtenus par l’élément 3D FE
du code ABAQUS, tandis que les graphes du modèle avec zigzag se rapprochent de
ceux du modèle ZIGT FE. D'autre part, le modèle avec DMQP avec zigzag s’avère plus
précis pour les structures épaisses. En revanche, si on admet les valeurs d'éléments
finis 3D comme référence, le modèle DMQP reflète généralement la même allure et
maintient un bon accord de précision.
97
DMQP DMQPz ZIGT FE TOT FE 3D FE
ABAQUS
5 11.1933 11.1138 11.0860 11.0731 11.0662 13.0298 12.9428 12.9123 12.8981 12.8905 12.4229 12.2336 12.1585 12.1197 12.0965 11.9516 11.4860
18.0946 17.8383 17.7486 17.7073 17.6848 20.6718 20.2412 20.0921 20.0234 19.9860 20.4106 19.8051 19.5636 19.4335 19.3528 19.0208 18.2956
19.0108 18.6000 18.4574 18.3915 18.3559 22.4415 22.1036 21.9854 21.9313 21.9016 22.0082 21.4534 21.2649 21.1762 21.1266 20.3969 19.3319
23.9499 23.4652 23.2965 23.2182 23.1761 27.8006 27.2482 27.0556 26.9670 26.9187 28.0214 26.8915 26.5140 26.3278 26.2179 25.4967 24.2519
27.0131 26.2057 25.9228 25.7924 25.7218 31.8491 30.5296 30.0726 29.8617 29.7480 31.6794 30.0722 29.5378 29.2653 29.0997 28.5190 27.3411
29.9232 28.6748 28.2414 28.0421 27.9340 34.5520 33.4675 33.0877 32.9124 32.8169 35.5384 33.9248 33.4451 33.2334 33.1200 30.9182 28.9117
10 17.5909 17.4449 17.3939 17.3703 17.3574 20.0191 19.8606 19.8050 19.7793 19.7652 18.4292 18.2471 18.1725 18.1342 18.1118 18.2744 17.7502
29.9004 28.9889 28.6764 28.5329 28.4555 31.9473 31.0161 30.6973 30.5513 30.4719 29.8103 29.4010 29.1333 29.1333 29.0729 28.9047 28.2032
31.7090 31.1787 30.9914 30.9052 30.8586 38.6041 37.9149 37.6717 37.5585 37.4973 34.7486 34.0295 33.6352 33.6352 33.5629 33.8184 32.4505
40.3478 39.3123 38.9537 38.7890 38.6988 46.4249 45.2748 44.8744 44.6897 44.5894 42.4472 41.5875 41.1102 41.1102 41.0151 41.0769 39.6697
49.7556 46.7393 45.7327 45.2721 45.0237 52.1211 49.0289 47.9751 47.4937 47.2343 47.7942 46.5290 45.8210 45.8210 45.6649 44.9767 43.7521
49.9554 48.3473 47.7562 47.4827 47.3341 62.4141 59.4705 58.4651 58.0052 57.7569 55.6810 53.7990 52.9209 52.9209 52.7698 52.7390 50.4550
20 23.8427 23.6385 23.5662 23.5325 23.5142 25.5770 25.3554 25.2766 25.2399 25.2199 24.0778 23.9687 23.9165 23.8871 23.8689 24.1130 23.7212
38.7178 37.3545 36.8909 36.6790 36.5639 40.1636 38.7978 38.3326 38.1209 38.0063 36.9491 36.8143 36.7516 36.7140 36.6889 36.7473 36.3153
50.2437 49.2302 48.8684 48.6995 48.6078 57.4447 56.1813 55.7275 55.5160 55.4001 51.6234 50.8759 50.5601 50.3919 50.2908 51.6151 49.6061
59.7250 58.0365 57.4481 57.1761 57.0284 66.3100 64.3267 62.5038 61.6795 61.2376 59.3967 58.8536 58.6198 58.4891 58.4074 59.1384 57.6707
68.1189 62.7402 60.9519 60.1430 59.7086 69.8169 64.4417 63.7848 63.4805 63.3150 60.3323 59.7326 59.5236 59.4102 59.3371 59.0253 58.2579
82.7835 77.7755 76.1016 75.3421 74.9343 88.5631 83.4531 81.7408 80.9628 80.5445 76.9830 76.1557 75.8823 75.7354 75.6414 75.9905 74.4577
98
Chapitre 3. Validation numérique du modèle de plaques composites DMQPz
15
Fréquence naturelle
DMQP
DMQPz
14
ZIGT FE
3D FE ABAQUS
TOT FE
13
1
12
576 éléments
11
10
9
8 10 12 14 16 18 20 22 24
38
Fréquence naturelle
37 DMQP
DMQPz
36 ZIGT FE
3D FE ABAQUS
35 TOT FE
34
33
6
32
31
30
29
576 éléments
28
27
26
8 10 12 14 16 18 20 22 24
99
Chapitre 3. Validation numérique du modèle de plaques composites DMQPz
22
Fréquence naturelle
DMQP
21 DMQPz
ZIGT FE
3D FE ABAQUS
TOT FE
20
19
1
18 576 éléments
17
16
15
8 10 12 14 16 18 20 22 24
66
Fréquence naturelle
DMQP
63
DMQPz
ZIGT FE
TOT FE
60 3D FE ABAQUS
57
6
54
51 576 éléments
48
45
42
8 10 12 14 16 18 20 22 24
100
Chapitre 4 Formulation théorique du modèle de coques
composites à effet zigzag
4.1 Introduction
Comme c'est déjà discuté dans le premier chapitre, et compte tenu le caractère
global ou local utilisé dans la description du champ de déplacement (ou de contrainte),
on peut distinguer différentes approches qui se basent généralement sur la théorie de
Kirchoff pour les structures minces ou sur la théorie de Reissner-Mindlin pour les
structures épaisses. Nous pouvons les classer soit en tant que approches globales,
dont la structure est considérée homogène équivalente, soit comme étant locales (ou
combinée: globale-locale) dites «par couche», dans lesquelles le champ de déplacement
ou de contrainte est dépendant de couche et varie linéairement ou plus suivant
l'épaisseur. Compte tenu de la performance numérique liée proportionnellement au
cout du calcul élevé, ces dernières approches sont généralement réservées aux
applications spécifiques dont la réponse locale est réclamée et/ou le comportement
local est dominant dans certains endroits.
Néanmoins, le modèle classique reste plus répandu dans la modélisation des
structures composites géométriquement complexes, d'où l'introduction d'une
correction sur les composantes de la rigidité en cisaillement est parfois nécessaire
pour améliorer le comportement à l'effort tranchant. Les résultats obtenus dépendent
essentiellement du choix des coefficients de correction, l'estimation des contraintes
dans quelques structures épaisses est incertaine et nécessite un modèle minutieux.
Afin de surmonter cet anomalie tout en gardant le caractère opératoire de ce
modèle homogène équivalent, on se propose comme but, d'étendre la notion de zigzag,
qui reflète un champ de déplacement par couche (déjà discutée en deuxième chapitre),
sur le modèle élément fini DDM pour les structures coques multicouches. Donc, la
formulation élément fini de deux versions pour la modélisation des coques composites
sera présentée dans la suite de travail. Les éléments sont des quadrilatères courbes à
quatre nœuds: basés sur le concept de Mindlin discret, ils tiennent compte de l'aspect
multicouche de la structure d'une part et l'effet de zigzag selon la démarche de
Murakami d'autre part.
101
Chapitre 4. Formulation théorique du modèle de coques composites à effet zigzag
Figure 4.1 Discrétisation géométrie d'une coque avec un élément quadrilatère courbe [6]
soit ( une base covariante (naturelle) au point p définie par deux vecteurs
tangents aux directions ( constant) et ( constant) et un troisième vecteur
perpendiculaire à leur plan, équation (4.2).
, , (4.2)
soit
(4.3)
un élément différentiel au point s'écrit:
(4.4)
sa longueur est donnée par le terme , avec:
(4.5)
102
Chapitre 4. Formulation théorique du modèle de coques composites à effet zigzag
; (4.6)
On définit également une autre base contravariante (duale) en même point , notée
( , figure (4.2), ou sous forme matricielle .
avec
(4.8)
; (4.9)
;
(4.10)
(4.11)
(4.12)
avec
(4.13)
103
Chapitre 4. Formulation théorique du modèle de coques composites à effet zigzag
(4.14)
avec: ; ;
cette matrice est obtenue par une rotation rigide du repère autour de l'axe
de telle sorte que le vecteur coïncide avec la normale après la rotation. étant
supposé non parallèle à . Si le dénominateur de A est nul, c-à-dire , alors la
matrice prendra la forme suivante:
, (4.15)
avec
, (4.16)
; (4.17)
; (4.18a,b)
avec
; (4.19a,b)
; ; (4.20a)
; (4.20b)
La matrice est une base covariante au point q de la coque, elle est linéaire en z et
ses vecteurs sont tangents à la surface située à une distance z de la surface
moyenne.
Un élément d'aire est défini par:
(4.21)
104
Chapitre 4. Formulation théorique du modèle de coques composites à effet zigzag
(4.24)
(4.26)
, (4.27)
H et K sont respectivement les courbures moyenne et totale (Gaussienne). La
substitution des relations (4.7) et (4.25) dans l'équation (4.22), permet ainsi d'écrire
l'élément de volume :
(4.28)
En utilisant la relation (4.23), l'inverse de s'écrit:
(4.29)
avec
; (4.30)
avec
(4.32)
105
Chapitre 4. Formulation théorique du modèle de coques composites à effet zigzag
(4.33)
(4.35)
avec
(4.36)
et
; (4.37a,b)
(4.38)
(4.40)
(4.42)
106
Chapitre 4. Formulation théorique du modèle de coques composites à effet zigzag
est définie par une interpolation quadratique incomplète, faisant apparaitre aux
milieux de chaque côté un accroissement de rotation . L'approximation nodale est
basée sur les fonctions de forme (tableau 2.1), de telle sorte que:
(4.43)
(4.46b)
En d'autres termes
(4.47)
avec
(4.48)
(4.49)
(4.50)
107
Chapitre 4. Formulation théorique du modèle de coques composites à effet zigzag
La confrontation des deux équations (4.51) et (4.52), nous permet d'aboutir finalement
à
(4.53)
avec (4.54)
ainsi les composantes de déformations virtuelles sont définies dans la base curviligne:
(4.55)
avec
(4.56)
(4.58)
avec
(4.59)
Où
(4.60)
(4.61)
(4.62)
108
Chapitre 4. Formulation théorique du modèle de coques composites à effet zigzag
(4.63)
(4.64)
(4.65)
En regroupant les termes ayant les mêmes coefficients, on obtient:
(4.66)
avec
(4.67)
(4.68)
(4.69)
(4.70)
(4.71)
d'où
; (4.72)
; (4.73)
; (4.74)
(4.75a)
109
Chapitre 4. Formulation théorique du modèle de coques composites à effet zigzag
(4.75b)
(4.76)
avec
(4.78)
(4.81)
avec
(4.82)
110
Chapitre 4. Formulation théorique du modèle de coques composites à effet zigzag
(4.83)
4.4.1.1. Déformations
La partie linéaire du vecteur de déformation est ainsi:
(4.84)
(4.86)
avec
(4.89)
et
et ;( rel. 3.37b)
; (4.90)
111
Chapitre 4. Formulation théorique du modèle de coques composites à effet zigzag
(4.91)
avec
(4.92)
avec
(4.94)
soit
; (4.95)
4.4.1.2. Déformations
Afin de limiter la complexité de l'expression de la déformation du cisaillement
due aux termes d'ordre deux, figurant dans le produit (4.83), nous adoptons dans la
suite du développement une forme simplifiée de , c'est ainsi que seuls les termes
constants et constants par couche seront retenus.
soit
on aura donc:
(4.96)
avec
(4.97)
soit
; (4.98a,b)
avec
(4.99)
112
Chapitre 4. Formulation théorique du modèle de coques composites à effet zigzag
on aura
(4.100)
avec
(4.101)
la condition du cisaillement constant par côté (dans le plan) se traduit par une
intégrale de contour, en utilisant la composante tangentielle du vecteur de
déformation, c'est ainsi:
(4.102)
on obtient
(4.103)
(rel. 2.60)
La combinaison des relations (4.103) et (2.60) nous permet de déduire les rotations
aux milieux des cotés.
d'où
(4.104)
avec
; , (4.105)
, (4.106)
(4.107)
; (4.108)
(4.109)
113
Chapitre 4. Formulation théorique du modèle de coques composites à effet zigzag
avec
(4.110)
(4.111)
; ; (4.112)
(4.113)
(4.114)
114
Chapitre 4. Formulation théorique du modèle de coques composites à effet zigzag
(4.115)
(4.116)
avec
(4.117)
et
(4.118)
(4.119)
avec sur
le travail virtuel interne est donné par:
(4.121)
(4.124)
115
Chapitre 4. Formulation théorique du modèle de coques composites à effet zigzag
avec
; (4.125)
; (4.126)
; (4.127)
; (4.128)
(4.129)
; (4.130)
et
(4.131)
(4.132)
(rel. 2.36a)
(4.133)
(rel. 2.36b)
(rel. 2.36c)
et
(4.134)
; (4.135)
; (4.136)
; (4.137)
avec
(rel. 2.37a)
(rel. 2.37b)
(rel. 2.37c)
116
Chapitre 4. Formulation théorique du modèle de coques composites à effet zigzag
Rappelons que la nullité des matrices constitutives liées aux termes de zigzag conduit
directement à une formulation simplifiée de l'élément DQMQS à 6 ddl/nœud. Les
conditions aux limites imposées devraient tenir en compte la suppression des lignes et
des colonnes liées aux variables nodales de zigzag. Une version du modèle DMQS sans
zigzag a été également élaborée afin d'optimiser le temps de calcul dans le cas
échéant.
(4.138)
avec
; ; (4.139)
Le même développement est retenu pour les termes virtuels, nous obtenons
(4.140)
(4.141)
est une valeur caractéristique associé à une rigidité de flexion (par exemple
. Pour le cas d'un multicouche, et peuvent être
remplacées par et .
est un coefficient petit ( dépendant de la précision de l'ordinateur
comme il peut être défini en fonction de la géométrie de l'élément [3].
L'ajout de la matrice à la matrice de rigidité élémentaire nous permet
d'éviter la singularité de sa forme globale . Concernant le modèle avec zigzag, la
même démarche a été adoptée lors de la projection des déplacements et dans le
repère cartésien global. Ceci conduit à l'apparition d'un degré de liberté additionnel,
également fictif suivant l'axe perpendiculaire au plan moyen. Contrairement au
sixième variable des modèles classiques, le degré de liberté du modèle en zigzag -
comme d'autres modèle d'ordre supérieur - ayant un sens physique un peu ambigu, en
revanche, son importance demeure indispensable sur le plan numérique.
117
Chapitre 4. Formulation théorique du modèle de coques composites à effet zigzag
d'où ; (4.143)
; : matrice unitaire
concernant le modèle avec zigzag, la matrice s'étend à la forme suivante:
(4.144)
(4.145)
avec
(4.147)
118
Chapitre 4. Formulation théorique du modèle de coques composites à effet zigzag
(4.149)
Or, les dérivées des contraintes planes (équation 4.146) seront définies à partir des
efforts généralisés , et .
soit
donc
(4.150)
(4.151)
Voyant les équations d'équilibre (2.30) données par Varadan et al. [189], qui révèle
l'entité définie par:
(4.152a)
avec
(4.152b)
et vice-versa, les allongements et les courbures peuvent être définis en fonction des
efforts généralisés, ils s'expriment ainsi:
(4.153a-b-c)
Par conséquent, ces trois expressions (4.153a-b-c) servent à définir en fonction des
efforts résultants. La substitution de leurs dérivées dans (4.148) conduit alors à:
119
Chapitre 4. Formulation théorique du modèle de coques composites à effet zigzag
(4.154)
avec
; (4.155a)
; (4.155b)
; (4.155c)
; (4.155d)
; (4.155e)
; (4.155f)
; (4.155g)
; (4.155h)
; (4.155i)
; (4.155j)
; (4.155k)
; (4.155l)
; (4.155m)
; (4.155n)
; (4.155o)
; (4.155p)
; (4.155q)
120
Chapitre 4. Formulation théorique du modèle de coques composites à effet zigzag
(4.155r)
avec
(4.156a)
(4.156b)
(Hypothèse)(4.156c)
et
(4.157a)
(4.157b)
(4.157c)
avec
, (4.159a)
, (4.159b)
(4.159c)
où
(4.160a)
(4.160b)
(4.160c)
121
Chapitre 4. Formulation théorique du modèle de coques composites à effet zigzag
(4.160d)
avec
; (4.161a,b)
; (4.162a,b)
et
(4.163a,b)
avec
; (4.164a,b)
Les dérivées secondes d'une fonction de forme sont généralement données par:
(4.165a)
(4.165b)
rappelons que cette forme est valable pour les fonctions utilisées dans
l'interpolation des rotations (Tableau 2.1). Cependant ces dérivées se réduisent à une
forme plus simple (4. 167c) quand il s'agit des fonctions de forme bilinéaire
(4.165c)
122
Chapitre 4. Formulation théorique du modèle de coques composites à effet zigzag
(4.166)
avec
(4.168a-b)
(4.168c)
on aura
(4.169a)
(4.169b)
(4.169c)
123
Chapitre 4. Formulation théorique du modèle de coques composites à effet zigzag
(4.170)
pour discrétiser cette énergie, nous pouvons réécrire les variables cinématiques sous
cette forme:
(4.171a-b)
(4.172a-b)
avec
– –
– – (4.172c)
– –
(4.172d)
et
(4.172c)
(4.172c)
(4.174)
124
Chapitre 4. Formulation théorique du modèle de coques composites à effet zigzag
avec
(4.175a)
(4.175b)
(4.175c)
(4.175d)
(4.175e)
125
Chapitre 5 Validation numérique des modèles de coques
composites avec et sans effet Zigzag
5.1 Introduction
Dans ce chapitre nous allons évaluer les performances des éléments de coques
DMQS (sans effet zigzag) et DMQSz (avec effet zigzag) à travers des tests standards
des coques isotropes, stratifiées et sandwichs. Selon les formulations détaillées dans
le chapitre précédent, l'étude comporte deux aspects: la partie statique dans laquelle
on traite la réponse des éléments formulés au comportement des structures
composites plate non-symétriques, coques cylindriques, vrillées, hémisphériques...etc.
d'une part, et la réponse aux vibrations libres de quelques panneaux de structures
cylindriques. Ces tests sont considérés indispensable à la validation des éléments finis
de coque selon NAFEMS (National Agency for Finite Element Methods and
Standards) [219]. Ils ont pour but de tester la convergence en augmentant le nombre
d'éléments et de vérifier la capacité de ces éléments à simuler des comportements
jugés complexes par l'ingénieur.
126
Chapitre 5. Validation numérique des modèles de coques composites
avec et sans effet Zigzag
Dans les tableaux 5.1 et 5.2 nous présentons respectivement les résultats pour
le déplacement normal et déplacement tangentiel en fonction de nombre
d'éléments par côté (N). Les solutions analytiques de Flügge [220] et Lindberg et al
[221].sont données par:
- Déplacement : =164.24
- Déplacement : =4.11
Figure 5.2 Cylindre pincé avec diaphragmes: déformée globale d'un huitième du
cylindre (12×12).
127
Chapitre 5. Validation numérique des modèles de coques composites
avec et sans effet Zigzag
50
-50
Wref
-100
-150 EhW/P
DMQS -165.63
-200
0 0.5 1 1.5 2 2.5 3
Longueur du côté DC (m)
100
50
0
Wref
-50
-100
EhW/P
-150
DMQS
-165.63
-200
B C
Points du contour
Figure 5.3 Cylindre pincé. L'allure du déplacement normal le long des côtés DC et
BC.
128
Chapitre 5. Validation numérique des modèles de coques composites
avec et sans effet Zigzag
Figure 5.4 Cylindre pincé avec diaphragmes, déformée globale, maillage avec
l'élément DMQS (32×100) sur Abaqus.
129
Chapitre 5. Validation numérique des modèles de coques composites
avec et sans effet Zigzag
Outre les solutions exactes citées dans [223], celles de Batoz et al.[224] fondées
sur la théorie des poutres sont proposées comme valeurs de référence. Elles sont
obtenues par un élément fini de poutre vrillée, en considérant un modèle variationnel
mixte où les équations d'équilibre interne sont satisfaites. Les résultats de deux
modèles DMQS et DMQSz sont également confrontés à ceux données par les éléments
DKT-CST [218], [69] et à l'élément mixte hybride naturel NHMiSP4/ml (la
partie membrane est basée sur une formulation hybride au sens de Pian) [6].
A part l'incapacité de l'élément DKQ24 à reproduire la réponse appropriée
dans le premier rapport d'élancement ( ) tableau 5.3, et en dépit du petit
dépassement constaté au niveau de l'élément quand il s'agit du chargement
hors plan (plan ), l’ensemble des éléments converge de façon monotone vers la
solution de référence pour les deux cas de charge.
Dans le cas du chargement hors plan, d'où l'effet de membrane au voisinage de
l'encastrement est dominant. L'élément DMQS (riche en flexion) donne des résultats
130
Chapitre 5. Validation numérique des modèles de coques composites
avec et sans effet Zigzag
131
Chapitre 5. Validation numérique des modèles de coques composites
avec et sans effet Zigzag
Données géométriques:
Propriétés mécaniques:
132
Chapitre 5. Validation numérique des modèles de coques composites
avec et sans effet Zigzag
Maillage
Eléments
N=2 N=4 N=6 N=8 N=10 N=12
DMQS 0.0960 0.0979 0.0960 0.0947 0.0940 0.0937
NHMiSP/ml - 0.1075 0.1027 0.9997 - 0.09733
0.084 0.093 0.093 0.093 0.093 0.093
DKQ24 0.096 0.098 0.096 0.095 0.094 0.094
Analytique [223] 0.094
Simo et al. [225] 0.093
Tableau 5.5 Hémisphère pincé. Convergence de
Données géométriques:
et
Propriétés mécaniques:
133
Chapitre 5. Validation numérique des modèles de coques composites
avec et sans effet Zigzag
134
Chapitre 5. Validation numérique des modèles de coques composites
avec et sans effet Zigzag
Nous constatons que les déplacements obtenus par les éléments DMQS et
DMQSz restent les plus proches à la solution de référence par apport aux autres
éléments, avec une erreur relative de l'ordre de 1,35%. L'effort normal au centre du
cylindre est également obtenu avec une bonne précision pour les deux versions du
modèle DDM. Les moments restent en bon accord avec les valeurs de référence, avec
une légère sensibilité due au cisaillement transverse. Pour approuver l'absence de la
fonction zigzag (une droite qui rejoint les deux interfaces de la section plane) dans ce
type des structures (une seule couche orthotrope), nous avons simulé la section comme
ayant trois couches fictives. Les résultats obtenus par le modèle avec zigzag,
développé sans correction et ceux du modèle primitif «DMQS» sont quasi-identiques,
d'où cette similitude réfère à la rationalité de cet hypothèse.
Y R = 0.508 m - L = 0.508 m,
Z 3 h = 0.0254 m
2 Matériau : 1 couche à 0°/x
E1 = 1.379 1010 Pa
1 E2 = 5.1711 1010 Pa
G12 = 8.6184 109 Pa
G13 = G23 = 4.3092 109 Pa
12= 0.0667
Conditions aux limites
sur AB : U=V=W=X=Y=Z=0
sur BC : W=X=Y=0
X sur AD : U=Y=Z=0
sur DC : V=X=Z=0
Chargement :
Pression interne p = 1 Pa
Théo. des
DMQS DMQSz DSQ24 DST18 DMTS MiSP4-Q4 Coques
minces
- Flèches max. aux
D: 0.681 D : 0.616
points D et C 0.6720 0.6719 0.677 0.684 0.663
C: 0.684 C : 0.685
( )
- Moment au point
D 2.323 2.318 2.14 2.01 - - 2.517
)
- Moment au point
A 8.229 8.050 6.72 4.24 - - 7.287
)
- Effort
circonférentiel
0.468 0.468 0.471 0.447 - - 0.461
au point D,
)
135
Chapitre 5. Validation numérique des modèles de coques composites
avec et sans effet Zigzag
Géométrie:
Matériau: (orthotrope), 3 couches (0°/0°/0°)
peau:
Cœur:
Conditions aux limites: sur AB, sur BC, sur CD,
sur DA
Chargement uniforme unitaire.
Figure 5.10 Plaque carrée sandwich non symétrique simplement supportée sous
charge uniformément répartie, données du problème.
136
Chapitre 5. Validation numérique des modèles de coques composites
avec et sans effet Zigzag
Srinivas [11] et celle de l'élément hétérosis conçu par Owen et al. [197]. Cependant,
tous les éléments à l'exception de l'élément mixte hybride NHMiSP/Q4/ml et celui
DMQS avec zigzag, révèlent des valeurs assez proches à la valeur de référence. Cette
performance se justifie par la correction exclusive de leur partie du cisaillement
transverse faite par des coefficients spécifiques et qui ne
pourrait pas être concrétisée par la forme zigzag ou le modèle mixte. Néanmoins, la
précision du modèle DMQSz concernant les contraintes planes se montre au niveau de
la section inférieure au plan moyen de la plaque, dont les résultats sont assez proches
aux valeurs de référence que les autres éléments.
D'autre part, il est bien montré sur la figure (5.11) la dissymétrie de l'allure de
la contrainte de cisaillement par rapport à l’axe moyen. Les courbes définies par les
deux modèles basés sur le concept «DDM» montrent des valeurs légèrement
distinguées par rapport à ceux des autres éléments finis.
Modèles
Solution 3D
34.55 63.76 50.24 5.04 -0.34 -3.40 -43.77
[11]
[11, 230]
Tableau 5.7 Plaque carré sandwich non symétrique simplement supportée sous charge
uniforme- Comparaison des déplacements.
0.5
0.4
DMQS
0.3 DMQSz
Héterosis
DSQ24
0.2 HMiSP4/Q4/ml
0.1
2 4 6 8 10 12 14 16
xz
-0.1
-0.2
-0.3
-0.4
-0.5
137
Chapitre 5. Validation numérique des modèles de coques composites
avec et sans effet Zigzag
Sol.
EHOST Elasticité
DMQSz DMQS NHMiSP4/ml analytique
[232] 3D [231]
[233]
16X20 20X16 16X20 20X16 20X20 20X6
10 1.1516 1.1917 1.3472 1.3888 1.4805 1.373 1.369 1.38
50 0.7396 0.7482 0.7466 0.7553 0.8045 0.7613 0.762 0.7622
100 0.6137 0.6188 0.6150 0.6201 0.7460 0.6256 0.6255 0.6261
138
Chapitre 5. Validation numérique des modèles de coques composites
avec et sans effet Zigzag
Le tableau (5.8) montre les résultats du modèle ainsi développé et qui est
défini par ses deux variantes DMQS et DMQSz. Des résultats sont également fournis
dans la référence [232], issus d'un élément quadrilatère à double courbure à 9 nœuds
et 5 ddl/noeud, proposé pour les coques composites stratifiées symétriques. Sa
formulation est basée sur une théorie d'ordre supérieur qui tient compte de l'effet de
zigzag (EHOST: Efficient Higher Order Shell Theory). Tous les résultats sont
comparés aux solutions analytiques de Ren et al. [233] et celles de l'élasticité 3D [231].
Les résultats sont présentés pour la flèche au centre du cylindre, en
considérant les rapports de selon la valeur adimensionnelle suivante:
Nous observons que le modèle DMQS donne de bons résultats par rapport aux
valeurs de référence. En particulier, l'élément primitif DMQS semble plus précis que
celui avec zigzag, Cette qualité est également confirmée vis à vis de l'élément
NHMiSP4/ml dans ce type de problème. Comme on s'y attendait, les résultats
convergent mieux vers la solution analytique pour les différents maillages quand le
rapport devient grand. Du fait que la structure subit un fort gradient de la charge et
qui nécessite un maillage plus raffiné, une meilleure précision est remarquée lorsque
l'accroissement de la charge est bien représenté sur la circonférence. Cela est exprimé
par un maillage de 20×16 éléments, tableau (5.8).
139
Chapitre 5. Validation numérique des modèles de coques composites
avec et sans effet Zigzag
C
F L
D B
E
S
R A
Z
Géométrie : R=609.6 mm , L=S=304.8 mm, h=3.048 mm
Y Matériau: E=206.843 Gpa; ν=0.3; ρ=7.833x103Kg/m3
Conditions aux limites : u=v=w=θx=θy= θz 0 sur DC
X
Maillage Elément/Mode 1 2 3 4 5 6
140
Chapitre 5. Validation numérique des modèles de coques composites
avec et sans effet Zigzag
D
L
A C
S
R B
Z
Y
X
Le tableau 5.10 regroupe les cinq premières valeurs de fréquences propres, des
représentations graphiques fréquences-modes sont également illustrées sur la figure
5.15. Nous pouvons constater une bonne cohérence entre nos résultats et ceux
obtenus numériquement par Bouabdalahh [64] et Crawley [234]. Ces résultats se
montrent plus fiables pour le calcul de la première fréquence dont le premier mode est
dominant. D'autre part, nous pouvons conclure, que les modèles numériques sont en
bon accord avec les valeurs empiriques en dépit de la légère différence qui se montre
pour les autres fréquences.
141
Chapitre 5. Validation numérique des modèles de coques composites
avec et sans effet Zigzag
142
Chapitre 5. Validation numérique des modèles de coques composites
avec et sans effet Zigzag
143
Chapitre 6 Modèle discret de coque non-linéaire
géométrique. Formulation et évaluation
6.1 Introduction
Dans tout les problèmes discutés précédemment, on a implicitement supposé
que tous deux, déplacements et déformations développés dans les structures sont
petits, en termes pratiques cela veut dire que la géométrie des éléments demeure
fondamentalement inchangée durant le processus de chargement et, au premier ordre
infinitésimal. Une approximation linéaire des déformations peut être utilisée.
En pratique de telles suppositions manquent fréquemment de régularité, les
structures coques soumises à des charges externes, présentent des comportements
très différents lorsque la non-linéarité géométrique devient importante. L'analyse
non-linéaire et la localisation de la limite et les points de bifurcation sont des étapes
importantes d'une étude complète et rigoureuse. Une étude bibliographique très
détaillée sur la non-linéarité est présentée par Frey [235]. Il a présenté les progrès,
qui ont permis aux chercheurs de résoudre les problèmes non-linéaires. Argyris et al
[236] ont été parmi les premiers qui se sont penchés sur la non-linéarité géométrique
pour laquelle ils ont proposé et adopté le principe. Dans la littérature, on peut trouver
de nombreux articles consacrés aux éléments de coque à facettes planes et qui sont
fondés généralement sur le concept du Kirchhoff. Vu le caractère épais de quelques
structures coques et qui exige la prise en compte de l'effet du cisaillement transverse,
les hypothèses de Mindlin demeurent une nécessité dans les développements futurs.
Selon cette thématique, un modèle coque basé sur un concept discret de Mindlin
modifié a été développé récemment [4]. La bonne performance de ce modèle à travers
les cas-tests linéaires de coque et les résultats satisfaisants, obtenus, dans le domaine
non-linéaire, sur la version "triangle à trois nœuds" du modèle DDM [208], nous ont
motivé pour étendre notre approche aux applications non-linéaires. Nous proposons
dans les sections suivantes de développer le modèle non-linéaire géométrique de coque
courbe à 4 nœuds DMQS pour un matériau homogène isotrope, avec l'hypothèse des
petites déformations élastiques, en formulation lagrangienne actualisée à chaque
itération.
144
Chapitre 6. Modèle discret de coque non-linéaire géométrique.
Formulation et évaluation.
Avec E * est le tenseur des déformations virtuelles de Green -Lagrange (avec CR
comme configuration de référence). Le tenseur S est celui des contraintes de Piola-
Kirchhoff de seconde espèce.
Wint s'écrit en fonction des tenseurs gradients de déformation réelle et virtuelle F et
F * respectivement ainsi:
Avec les hypothèses des petites déformations, d'où on peut écrire S σt et
E ε L L .
* *
t
1
2
*
t t
ΔWint est représenté comme la somme de deux termes: une partie linéaire et une
deuxième partie dite géométrique.
ΔWint ΔW l ΔW g (6.4)
Avec
k k k
e
t
e
l
e
σ (6.6)
145
Chapitre 6. Modèle discret de coque non-linéaire géométrique.
Formulation et évaluation.
ΔW g Tr Lit σti ΔLit dV i un* kσi Δun
T
(6.7)
Vei
En utilisant les propriétés des traces [218], la définition de la matrice géométrique est
donnée, selon [237], par:
3
kσi M 0i k EF1i M 0i k dA i , k 1,3
T
(6.8)
k 1 Aei
t1 N j ,x 0 0 t2 N j ,x 0 0
M i
0 1
i
t1 N j , y 0 0 j 1, 4 , M 0 2 t2 N j , y 0 0 j 1, 4
0 00 0 00
n N j ,x 0 0
M i
0 3 n N j , y 0 0 j 1, 4
0 00
(6.9)
avec
0 0
N j ,x N j ,ξ C11 N j ,η C21
(6.10)
0 0
N j , y N j ,ξ C12 N j ,η C22
146
Chapitre 6. Modèle discret de coque non-linéaire géométrique.
Formulation et évaluation.
h2 N xi i
N xy Txi
EF1i σti dz N xyi i
N yy
Tyi (6.11)
h 2
Txi Tyi 0
Rappelons que les efforts normaux uniquement sont pris en considération dans le
calcul de la matrice géométrique (équation 6.8).
W i un* r i U *n R i (6.12)
Ve
147
Chapitre 6. Modèle discret de coque non-linéaire géométrique.
Formulation et évaluation.
z
y
L
P
x h
A
148
Chapitre 6. Modèle discret de coque non-linéaire géométrique.
Formulation et évaluation.
149
5 pas 10 pas 40 pas 100 pas
Tableau 6.3 Erreurs relatives sur le déplacement axial U pour les différents pas (10 éléments)
150
Chapitre 6. Modèle discret de coque non-linéaire géométrique.
Formulation et évaluation.
D'autre part, pour étudier l'influence du nombre de pas sur la précision de l'élément
DMQS, Les mêmes paramètres du problème précédent ont été adoptés en utilisant un
maillage de 10 éléments. Quatre incréments de charge identiques de 5, 10, 40 et 100
ont été utilisés à cet effet. Les résultats numériques des déplacements U et W au
point A et ceux relatifs aux éléments déjà mentionnés, et qui sont figurés par leurs
erreurs relatives par rapport aux valeurs de Timoshenko et Gere [238] sont reportés
dans les tableaux (6.3) et (6.4). Les résultats obtenus par le modèle DMQS présentent
des valeurs de même ordre que les autres éléments triangulaires présentés par
Hammadi [237], rappelons que ces éléments ont montré une bonne précision par
rapport aux deux meilleurs éléments de l'étude de Fafard, reportée dans [239].
Considérons les remarques précitées, L'élément DMQS montre une excellente
aptitude à représenter les grands déplacements et les grandes rotations avec une
formulation lagrangienne actualisée à chaque itération.
151
Chapitre 6. Modèle discret de coque non-linéaire géométrique.
Formulation et évaluation.
Les éléments que nous avons retenus dans cette étude sont les éléments triangulaires
à trois nœuds: DKT18 et le T3 15 développés par Hammadi et Batoz [237] ainsi que
l'élément quadrilatère DKQ24 [241] à 4 nœuds et 6ddl/nœud. Prenant comme
référence les résultats obtenus par Surana [242] par son élément quadrilatère à 16
nœuds. Du fait que la solution non linéaire et les points critiques sont bien connues,
nous utilisons dans le premier cas un pilotage en déplacement avec les mêmes
paramètres que ceux de Surana [242] qui a utilisé 12 incréments de déplacements
(WA=2,5 mm) et un pilotage en longueur d'arc pour le deuxième cas.
Pour le premier cas, les résultats juxtaposés de l'évolution de la charge en
fonction du déplacement de deux points A et B obtenus par le présent modèle sont
illustrés sur la figure 6.3. Nous constatons que l'élément DMQS modifié donne des
résultats concordants avec tous les éléments triangulaires et qui sont très proches à
ceux donnés par Surana [242]. Nous notons la présence d'une manière évidente du
phénomène de retour en charge défini par deux points limites.
La figure 6.4 montre les courbes charge-déplacement vertical du point central A pour
les différents éléments, ces courbes confirment une forte instabilité présentée par
deux points limites en charge et deux points limites en déplacement «snap through/
snap-back». Excepté le «snap-back» mal présenté dans les résultats de Surana [242],
Les résultats obtenus par le présent modèle sont proches de référence [242] et
concordent avec ceux des éléments triangulaires de Hammadi [237]. il est à noter que
l'allure de convergence de l'élément quadrilatère de Batoz [241] diffère légèrement de
celle d'autres éléments au voisinage du point d'inflexion (quasi-identique) de toutes
les courbes
4
DMQS (Point A) DMQS (Point B)
DKT15 (Point A) DKT3 (Point B)
T3G15 (Point A) T3G15 (Point B)
Chargement central P (KN)
0
0 5 10 15 20 25 30
Déplacement (mm)
152
DMQS DKT15 [237] T3 [237] Surana [242]
153
Chapitre 6. Modèle discret de coque non-linéaire géométrique.
Formulation et évaluation.
1.5
0.5
0.25
-0.25
-0.5
0 5 10 15 20 25 30 35
Déplacement W (mm)
154
Chapitre 6. Modèle discret de coque non-linéaire géométrique.
Formulation et évaluation.
0.0000 .0000 0.0 0.0000 0.0000 0.0000 0.0000 0.0000 0.0 0.0000
-2.4397 .1953 -1.7 -0.150 -3.4395 -0.2590 -3.4502 -0.2663 -2.5 -0.2101
-4.7239 .3268 -3.5 -0.265 -6.5276 -0.4103 -6.5631 -0.4198 -5.0 -0.3479
-6.8243 .4202 -4.9 -0.345 -9.2189 -0.5086 -9.2800 -0.5178 -7.5 -0.4528
-8.7334 .4899 -6.8 -0.410 -11.5160 -0.5694 -11.5940 -0.5761 -10.0 -0.5319
-10.4514 .5414 -8.4 -0.475 -13.3970 -0.5897 -13.4743 -0.5910 -12.5 -0.5817
-11.9754 .5752 -9.8 -0.520 -14.7935 -0.5611 -14.8463 -0.5541 -15.0 -0.5555
-13.2892 .5888 -11.1 -0.555 -15.6995 -0.4858 -15.7238 -0.4701 -16.0 -0.4694
-14.3627 .5789 -12.2 -0.580 -16.2825 -0.3809 -16.2948 -0.3582 -17.0 -0.1264
-15.1722 .5438 -13.1 -0.595 -16.6943 -0.2627 -16.6971 -0.2344 -18.0 0.3701
-15.7456 .4875 -14.0 -0.600 -16.9561 -0.1395 -16.9264 -0.1067 -19.0 0.3486
-16.1586 .4168 -14.9 -0.585 -16.9831 -0.0166 -16.8616 0.0186 -20.0 0.3168
-16.4762 .3380 -15.5 -0.565 -16.5853 0.0969 -16.2688 0.1275 -21.0 0.2636
-16.7277 .2548 -16.1 -0.525 -15.6016 0.1817 -15.1549 0.2025 -22.0 0.2205
-16.9116 .1696 -16.4 -0.480 -14.6166 0.2506 -14.4795 0.2817 -23.0 0.1542
-17.0029 .0838 -16.7 -0.415 -14.3570 0.3371 -14.8100 0.3648 -24.0 0.0744
-16.9544 -.0006 -16.9 -0.350 -15.6401 0.3847 -16.7449 0.3839 -25.0 -0.0198
-16.6933 -.0800 -17.0 -0.290 -18.0355 0.3637 -19.4317 0.3327 -26.0 -0.1298
-16.1456 -.1479 -17.1 -0.225 -20.9427 0.2698 -22.5282 0.1871 -27.0 -0.2565
-15.3684 -.1992 -17.1 -0.150 -24.1609 0.0620 -25.8250 -0.1048 -28.0 -0.4011
-14.6811 -.2477 -17.0 -0.090 -27.5001 -0.3102 -29.1605 -0.5909 -29.0 -0.5647
-14.3442 -.3073 -16.8 -0.020 -30.8339 0.8845 -32.4612 -1.3050 -30.0 -0.7485
-14.5221 -.3628 -16.4 0.035
-15.7111 -.3850 -16.0 0.085
-17.3654 -.3741 -15.3 0.130
-19.2891 -.3311 -14.8 0.155
-21.4100 -.2457 -14.2 0.195
-23.6644 -.1013 -14.0 0.240
-25.9920 .1194 -14.0 0.295
-28.3467 .4302 -14.3 0.345
-30.7021 .8418 -15.0 0.370
-16.1 0.380
-17.3 0.375
-18.7 0.350
-20.3 0.305
-21.8 0.230
-23.5 0.120
-25.2 -0.025
-26.8 -0.210
-28.5 -0.445
155
Chapitre 6. Modèle discret de coque non-linéaire géométrique.
Formulation et évaluation.
156
Chapitre 6. Modèle discret de coque non-linéaire géométrique.
Formulation et évaluation.
6
DMQS DKT15
DKT18 T3G15
5 Q4G24
Déolacement (mm)
0
0 50 100 150 200 250 300 350 400
Charge F(N)
157
Chapitre 6. Modèle discret de coque non-linéaire géométrique.
Formulation et évaluation.
12
DMQS DKT15
DKT18 T3G15
Q4G24
Déplacement (mm) 10
0
0 50 100 150 200 250 300 350 400
Charge F (N)
158
Chapitre 6. Modèle discret de coque non-linéaire géométrique.
Formulation et évaluation.
Configuration initiale
Configuration initiale (azimut=0°, elévation=60°)
(azimut=0°, élévation=-60°)
159
Chapitre 6. Modèle discret de coque non-linéaire géométrique.
Formulation et évaluation.
Les charges correspondantes à chaque pas obtenus par le modèle DMQS sont
confrontées à celles calculées par les éléments triangulaires DKT15 et T315 [237] et
l'élément quadrilatère [245] (avec 13 incréments de valeur égale à 0,8575 m). Les
courbes reportées sur la figure (6.10) montrent une excellente convergence du modèle
DMQS qui se déduit facilement de la concordance des résultats exprimés pour les
deux types de maillage. En revanche, la convergence de l'élément T315, dans ce cas,
160
Chapitre 6. Modèle discret de coque non-linéaire géométrique.
Formulation et évaluation.
s'avère un peu lente que celle des autres éléments. Toutefois, tous les éléments
convergent vers les résultats de référence [245] pour un maillage raffiné.
10
9
DMQS 4 x 6 DKT15 10 x 12
DKT15 4 x 6 T3G15 10 x 12
8 T3G15 4 x 6 Brendel & Ramm
DMQS 10 x 12
7
6
Charge (KN)
0
0 2 4 6 8 10 12
W/h au centre A
161
Conclusion générale et perspectives
162
Conclusion générale
déplacement linéaire par morceaux est formulé. Les deux éléments finis
correspondants baptisés DMQPz et DMQSz comportent outre les variables
cinématiques du modèle standard, des termes additionnels de type déplacement.
Cette dernière formulation est élaborée sans faire recours à la technique de
correction. Les tests de validation ont montrés une bonne cohérence de la deuxième
version du modèle DDM avec les modèles d'ordre supérieur et ceux du premier ordre
améliorés.
En dépit du caractère global (monocouche équivalente) de la majorité des modèles du
premier ordre, le modèle DMQP s’avère avantageux dans la représentation des
contraintes planes et une performance légèrement meilleure est montrée sur la
version en zigzag. Contrairement au modèle des plaques, la contribution de l'effet de
zigzag sur le comportement des coques composites est plus complexe à analyser. En
effet, le modèle DMQS seul possède une bonne capacité à reproduire le comportement
global de ce type des structures. Cependant lorsque la fonction de zigzag est ajoutée,
le comportement des coques multicouches ne permet pas de distinguer explicitement
l'effet des termes secondaires sur le modèle primitif. Or, comme les autres modèles
d'ordre supérieur, dont la performance est spécifique, le modèle DMQSz dépend des
conditions géométriques et mécaniques de la structure à étudier (rapport
d'élancement longueur/épaisseur, nombre de couches, rapport d'hétérogénéité, type de
charges extérieures, condition aux limites...).
Finalement, et comme c’est prévu au début de la thèse, les routines
développées dans le code Reflex sont implémentées dans le code de calcul industriel
ABAQUS. On s'est contenté essentiellement de l'aspect linéaire du modèle DMQS
multicouches/sandwichs. Une interface graphique a été élaborée également avec le
post- traitement du logiciel ABAQUS, ce travail ainsi que l'adaptation du modèle à la
géométrie non-linéaire est le couronnement de cette thèse.
En guise de perspectives à ce travail doctoral de recherche, nous proposons
d'étendre l'analyse non-linéaire aux cas des structures composites multicouches, en
intégrant non seulement la fonction "Zigzag" mais aussi les aspects liés au
comportement endommageable des matériaux sous sollicitations statiques et
dynamiques (fatigue et choc).
163
Références bibliographiques
[1] Katili I. Formulation et évaluation de nouveaux éléments finis pour l'analyse linéaire des
plaques et coques de forme quelconque [Ph.D Thesis]: University of Technology of
Compiègne,France; 1993.
[2] Katili I. A new discrete Kirchhoff-Mindlin element based on Mindlin–Reissner plate theory
and assumed shear strain fields—part II: An extended DKQ element for thick-plate bending
analysis. International Journal for Numerical Methods in Engineering. 1993;36:1885-908.
[3] Ayad R. Contribution à la Modélisation numérique pour l’analyse des solides et des
structures, et pour la mise en forme des fluides non newtoniens. Application à des matériaux
d'emballage [H.D.R]: University of Reims Champagne Ardenne; 2002.
[4] Sakami S. Modélisation numérique des structures composites multicouches à l’aide d’une
approche discrète au sens de Mindlin. Le modèle DDM (Displacement Discrete Mindlin)
[Doctorat Thesis]: University of Reims Champagne Ardenne; 2008.
[5] Murakami H. Laminated composite plate theory with improved in-plane responses. Journal
of Applied Mechanics. 1986;53:661-6.
[6] Tafla A. Eléments finis mixtes-hybrides naturels sans facteurs correctifs du cisaillement
pour les plaques et les coques composites multicouches [Ph.D Thesis]. France: University of
Reims Champagne Ardenne; 2007.
[7] Aboura Z, Talbi N, Ayad R. Une méthode d’homogénéisation adaptée au carton ondulé.
research report University of Technologie of Compiègne; 2002.
[8] Ayad R, Erre D, Hendou H, Odof S, Prudhomme JC. Adaptation du modèle DDM aux
structures orthotropes. Application aux vibrations libres d’une plaque en carton ondulé.
16ème Congrès Français de Mécanique. Nice, France.
[9] Robbins DH, Reddy JN. Modelling of thick composites using a layerwise laminate theory.
International Journal for Numerical Methods in Engineering. 1993;36:655-77.
[10] Robbins DH, Reddy JN. Analysis of piezoelectrically actuated beams using a layer-wise
displacement theory. Computers & Structures. 1991;41:265-79.
[11] Srinivas S. A refined analysis of composite laminates. Journal of Sound and Vibration.
1973;30:495-507.
[12] Seide P. An improved approximate theory for the bending of laminated plates. Mechanics
today. 1980;5:451-66.
[13] Barbero EJ. On a generalized laminate theory with application to bending, vibration, and
delamination buckling in composite laminates. Virginia Polytechnic Inst. and State Univ.,
Blacksburg, VA (USA); 1989.
[14] Barbero E, Reddy J, Teply J. An accurate determination of stresses in thick laminates
using a generalized plate theory. International Journal for Numerical Methods in Engineering.
1990;29:1-14.
[15] Bert CW. A critical evaluation of new plate theories applied to laminated composites.
Composite Structures. 1984;2:329-47.
[16] Noor AK, Burton WS. Assessment of Shear Deformation Theories for Multilayered
Composite Plates. Applied Mechanics Reviews. 1989;42:1-13.
[17] Noor AK, Burton WS. Assessment of computational models for multilayered anisotropic
plates. Composite Structures. 1990;14:233-65.
[18] Di Sciuva M. A general quadrilateral multilayered plate element with continuous
interlaminar stresses. Computers & Structures. 1993;47:91-105.
164
Références bibliographiques
[19] Jones RM. Mechanics of Composite Materials. New York: Mc Graw-Hill; 1975.
[20] Whitney J, Pagano N. Shear deformation in heterogeneous anisotropic plates. Journal of
Applied Mechanics. 1970;37:1031.
[21] Giavotto V. Sulla meccanica di pannelli anisotropi cd eterogenei. Istituto Lombardo
(Memorie Sc. Mat.) XXV; 1969. p. 437-80.
[22] Whitney JM, Sun C. A higher order theory for extensional motion of laminated
composites. Journal of Sound and Vibration. 1973;30:85-97.
[23] Levinson M. An accurate, simple theory of the statics and dynamics of elastic plates.
Mechanics Research Communications. 1980;7:343-50.
[24] Murthy M. An improved transverse shear deformation theory for laminated anisotropic
plates. 1981;TP-1903.
[25] Reddy JN. Simple higher-order theory for laminated composite plates. Journal of Applied
Mechanics, Transactions ASME. 1984;51:745-52.
[26] Krishna Murty AV, Vellaichamy S. On higher order shear deformation theory of
laminated composite panels. Composite Structures. 1987;8:247-70.
[27] Librescu L, Khdeir A. Analysis of symmetric cross-ply laminated elastic plates using a
higher-order theory: Part I--Stress and displacement. Composite Structures. 1988;9:189-213.
[28] Lo KH, Christensen RM, Wu EM. A High-Order Theory of Plate Deformation---Part 2:
Laminated Plates. Journal of Applied Mechanics. 1977;44:669-76.
[29] Carrera E. Theories and finite elements for multilayered, anisotropic, composite plates
and shells. Archives of Computational Methods in Engineering. 2002;9:87-140.
[30] Medwadowski SJ. A refined theory of elastic orthotropic plates. Journal of Applied
Mechanics. 1958;25:437-43.
[31] Yang PC, Norris CH, Stavsky Y. Elastic wave propagation in heterogeneous plates.
International Journal of Solids and Structures. 1966;2:665-84.
[32] Mindlin RD. Influence of rotatory inertia and shear on flexural motions of isotropic, elastic
plates. Journal of Applied Mechanics. 1951;18:31-8.
[33] Koiter WT. The Theory of Thin Elastic Shells. Proceedings of First IUTAM Symposium,
North-Holland1960.
[34] Vlasov BF. On the equations of Bending of plates. Dokla Ak Nauk Azerbeijanskoi-SSR.
1957;3:955-79.
[35] Hildebrand F, Reissner E, Thomas G. Notes on the foundations of the theory of small
displacements of orthotropic shells. D.C Washington: NACA; 1949.
[36] Nelson RB, Lorch DR. A Refined Theory for Laminated Orthotropic Plates. Journal of
Applied Mechanics. 1974;41:177-83.
[37] Librescu L. Elastostatics and kinetics of anisotropic and heterogeneous shell-type
structures: Springer; 1975.
[38] Lo KH, Christensen RM, Wu EM. A High-Order Theory of Plate Deformation---Part 1:
Homogeneous Plates. Journal of Applied Mechanics. 1977;44:663-8.
[39] Bickford WB. A consistent higher order beam theory. Development in Theoretical and
Applied Mechanics. 1982;11:137-50.
[40] Reddy JN. A refined nonlinear theory of plates with transverse shear deformation.
International Journal of Solids and Structures. 1984;20:881-96.
[41] Shi G. A new simple third-order shear deformation theory of plates. International Journal
of Solids and Structures. 2007;44:4399-417.
[42] Reddy JN. Energy and variational methods in applied mechanics: With an introduction to
the finite element method: Wiley (New York); 1984.
[43] Reddy JN. Mechanics of laminated composite plates- Theory and analysis. Boca Raton,
FL: CRC Press. 1997.
165
Références bibliographiques
[44] Reddy JN. Analysis of functionally graded plates. International Journal for Numerical
Methods in Engineering. 2000;47:663-84.
[45] Kant T. Numerical analysis of thick plates. Computer Methods in Applied Mechanics and
Engineering. 1982;31:1-18.
[46] Rohwer K. Application of higher order theories to the bending analysis of layered
composite plates. International Journal of Solids and Structures. 1992;29:105-19.
[47] Ferreira AJM, Roque CMC, Jorge RMN. Static and free vibration analysis of composite
shells by radial basis functions. Engineering Analysis with Boundary Elements. 2006;30:719-
33.
[48] Kant T, Swaminathan K. Analytical solutions for the static analysis of laminated
composite and sandwich plates based on a higher order refined theory. Composite Structures.
2002;56:329-44.
[49] Matsunaga H. Assessment of a global higher-order deformation theory for laminated
composite and sandwich plates. Composite Structures. 2002;56:279-91.
[50] Pagano NJ. Exact Solutions for Rectangular Bidirectional Composites and Sandwich
Plates. Journal of Composite Materials. 1970;4:20-34.
[51] Matsunaga H. A comparison between 2-D single-layer and 3-D layerwise theories for
computing interlaminar stresses of laminated composite and sandwich plates subjected to
thermal loadings. Composite Structures. 2004;64:161-77.
[52] Zhen W, Wanji C. An efficient higher-order theory and finite element for laminated plates
subjected to thermal loading. Composite Structures. 2006;73:99-109.
[53] Zenkour AM. Analytical solution for bending of cross-ply laminated plates under thermo-
mechanical loading. Composite Structures. 2004;65:367-79.
[54] Swaminathan K, Patil SS. Higher order refined computational model with 12 degrees of
freedom for the stress analysis of antisymmetric angle-ply plates – analytical solutions.
Composite Structures. 2007;80:595-608.
[55] Cho M, Kim JS. Improved Mindlin plate stress analysis for laminated composites in finite
element method. AIAA journal. 1997;35:587-90.
[56] Wanji C, Zhen W. A selective review on recent development of displacement-based
laminated plate theories. Recent Patents on Mechanical Engineering. 2008;1:29-44.
[57] Ayad R, Dhatt G, Batoz JL. A new hybrid‐mixed variational approach for Reissner–
Mindlin plates. The MiSP model. International Journal for Numerical Methods in
Engineering. 1998;42:1149-79.
[58] Ayad R, Rigolot A, Talbi N. An improved three-node hybrid-mixed element for
Mindlin/Reissner plates. International Journal for Numerical Methods in Engineering.
2001;51:919-42.
[59] Ayad R, Rigolot A. An improved four-node hybrid-mixed element based upon Mindlin's
plate theory. International Journal for Numerical Methods in Engineering. 2002;55:705-31.
[60] Ayad R, Rigolot A. Un nouvel élément mixte-hybride naturel pour les coques courbes. Le
modèle NHMiSP4. 16ème Congrès Français de Mécanique. Nice- France.
[61] Briossilis D. The four node Co Mindlin plate bending elements reformulated. Part I:
formulation. Computer Methods in Applied Mechanics and Engineering. 1993;107:23-43.
[62] Bathe K-J, Dvorkin EN. A four-node plate bending element based on Mindlin/Reissner
plate theory and a mixed interpolation. International Journal for Numerical Methods in
Engineering. 1985;21:367-83.
[63] Lardeur P. Développement et évaluation de deux nouveaux éléments finis de plaques et
coques composites avec influence du cisaillement transversal [Ph.D Thesis]: University of
Technology of Compiègne,France; 1990.
[64] Bouabdallah M. Modélisation de coques cylindriques raidies, isotropes et composites [Ph.D
Thesis]: University of Technology of Compiègne,France; 1992.
166
Références bibliographiques
167
Références bibliographiques
[87] Sheikh AH, Chakrabarti A. A new plate bending element based on higher-order shear
deformation theory for the analysis of composite plates. Finite Elements in Analysis and
Design. 2003;39:883-903.
[88] Nayak AK, Moy SSJ, Shenoi RA. Free vibration analysis of composite sandwich plates
based on Reddy's higher-order theory. Composites Part B: Engineering. 2002;33:505-19.
[89] Dvorkin EN, Bathe KJ. A continuum mechanics based four-node shell element for general
nonlinear analysis. Engineering Computations. 1984;1:77-88.
[90] Huang H, Hinton E. A nine node Lagrangian Mindlin plate element with enhanced shear
interpolation. Engineering Computations. 1984;1:369-79.
[91] Nayak AK, Moy SSJ, Shenoi RA. A higher order finite element theory for buckling and
vibration analysis of initially stressed composite sandwich plates. Journal of Sound and
Vibration. 2005;286:763-80.
[92] Khare RK, Kant T, Garg AK. Free vibration of composite and sandwich laminates with a
higher-order facet shell element. Composite Structures. 2004;65:405-18.
[93] Pervez T, Seibi AC, Al-Jahwari FKS. Analysis of thick orthotropic laminated composite
plates based on higher order shear deformation theory. Composite Structures. 2005;71:414-22.
[94] Goswami S. A C° plate bending element with refined shear deformation theory for
composite structures. Composite Structures. 2006;72:375-82.
[95] Carrera E, Antona E. A Class of Two-dimensional Theories for Anisotropic Multilayered
Plates Analysis [Memorie Scienze Fisiche]: Accademia delle Scienze di Torino; 1995.
[96] Carrera E. CZ° requirements—models for the two dimensional analysis of multilayered
structures. Composite Structures. 1997;37:373-83.
[97] LIBRESCU L. Refined geometrically nonlinear theories of anisotropic laminated shells.
Quarterly of Applied Mathematics 1987;55:1-22.
[98] Reply by the author to C.T. Sun. AIAA journal. 1994;32:2135-6.
[99] Carrera E. Developments, ideas, and evaluations based upon Reissner's Mixed Variational
Theorem in the modeling of multilayered plates and shells. Applied Mechanics Reviews.
2001;54:301-28.
[100] Lekhnitskii S. Strength calculation of composite beams1935.
[101] Ambartsumian SA. On a theory of bending of anisotropic plates. Investiia Akad Nauk
SSSR, Ot Tekh Nauk1958.
[102] Ambartsumian SA. Theory of anisotropic plates: Tech. Pub. Co; 1969.
[103] Reissner E. On a certain mixed variational theory and a proposed application.
International Journal for Numerical Methods in Engineering. 1984;20:1366-8.
[104] Reissner E. On a certain mixed variational theorem and on laminated elastic shell
theory. Proceedings of the Euromech-Colloquium 1986. p. 17–27.
[105] Carrera E. Historical review of Zig-Zag theories for multilayered plates and shells.
Applied Mechanics Reviews. 2003;56:287-308.
[106] Lekhnitskii SG. Anisotropic Plates1968.
[107] Ren JG. A new theory of laminated plate. Composites Science and Technology.
1986;26:225-39.
[108] Ren JG. Bending theory of laminated plate. Composites Science and Technology.
1986;27:225-48.
[109] Stürzenbecher R, Hofstetter K. Bending of cross-ply laminated composites: An accurate
and efficient plate theory based upon models of Lekhnitskii and Ren. Composite Structures.
2011;93:1078-88.
[110] Whitney JM. The Effect of Transverse Shear Deformation on the Bending of Laminated
Plates. Journal of Composite Materials. 1969;3:534-47.
[111] Rath BK, Das YC. Vibration of layered shells. Journal of Sound and Vibration.
1973;28:737-57.
168
Références bibliographiques
[112] Di Sciuva M. Bending, vibration and buckling of simply supported thick multilayered
orthotropic plates: An evaluation of a new displacement model. Journal of Sound and
Vibration. 1986;105:425-42.
[113] Di Sciuva M. A refinement of the transverse shear deformation theory for multilayered
orthotropic plates. l’Aerotecnica Missili Spazio. 1984;63:84-92.
[114] Di Sciuva M. A refinement of the transverse shear deformation theory for multilayered
anisotropic plates. Atti del Dipartimento di Ingegneria Aeronaurica e Spaziale del Politecnico
di Torino 1983;5.
[115] Di Sciuva M. A refined transverse shear deformation theory for multilayered anisotropic
plates. Atti Academia delle Scienze di Torino. 1984;118:279-95.
[116] Di Sciuva M. Multilayered anisotropic plate models with continuous interlaminar
stresses. Composite Structures. 1992;22:149-67.
[117] Cho M, Parmerter RR. Efficient higher order composite plate theory for general
lamination configurations. AIAA journal. 1993;31:1299-306.
[118] Di Sciuva M. Further refinement in the transverse shear deformation theory for
multilayered composite plates. Atti Accademia delle Scienze di Torino. 1990;124:248-68.
[119] Oh J, Cho M, Kim JS. Enhanced lower-order shear deformation theory for fully coupled
electro-thermo-mechanical smart laminated plates. Smart Materials and Structures.
2007;16:2229-41.
[120] Cho M, Oh J. Higher order zig-zag plate theory under thermo-electric-mechanical loads
combined. Composites Part B: Engineering. 2003;34:67-82.
[121] Cho M, Oh J. Higher order zig-zag theory for fully coupled thermo-electric–mechanical
smart composite plates. International Journal of Solids and Structures. 2004;41:1331-56.
[122] Oh J, Cho M. A finite element based on cubic zig-zag plate theory for the prediction of
thermo-electric-mechanical behaviors. International Journal of Solids and Structures.
2004;41:1357-75.
[123] Demasi L. Refined multilayered plate elements based on Murakami zig–zag functions.
Composite Structures. 2005;70:308-16.
[124] Carrera E. On the use of the Murakami's zig-zag function in the modeling of layered
plates and shells. Computers & Structures. 2004;82:541-54.
[125] Di Sciuva M. An improved shear-deformation theory for moderately thick multilayered
anisotropic shells and plates. ASME, Transactions, Journal of Applied Mechanics.
1987;54:589-96.
[126] Di Sciuva M. Development of an anisotropic, multilayered, shear-deformable rectangular
plate element. Computers & Structures. 1985;21:789-96.
[127] Di Sciuva M. Evaluation of some multilayered, shear-deformable plate elements.
Computers & Structures. 1986;24:845-54.
[128] Di Sciuva M. A third-order triangular multilayered plate finite element with continuous
interlaminar stresses. International Journal for Numerical Methods in Engineering.
1995;38:1-26.
[129] Reddy JN. A generalization of two-dimensional theories of laminated composite plates.
Communications in Applied Numerical Methods. 1987;3:173-80.
[130] Reddy JN. On refined computational models of composite laminates. International
Journal for Numerical Methods in Engineering. 1989;27:361-82.
[131] Cho M, Parmerter R. Finite element for composite plate bending based on efficient
higher order theory. AIAA journal. 1994;32:2241-8.
[132] Cho M, Parmerter RR. An efficient higher-order plate theory for laminated composites.
Composite Structures. 1992;20:113-23.
[133] Cho YB, Averill RC. First-order zig-zag sublaminate plate theory and finite element
model for laminated composite and sandwich panels. Composite Structures. 2000;50:1-15.
169
Références bibliographiques
[134] Topdar P, Sheikh AH, Dhang N. Finite Element Analysis of Composite and Sandwich
Plates Using a Continuous Inter-laminar Shear Stress Model. Journal of Sandwich Structures
and Materials. 2003;5:207-31.
[135] Pandit MK, Singh BN, Sheikh AH. Buckling of laminated sandwich plates with soft core
based on an improved higher order zigzag theory. Thin-Walled Structures. 2008;46:1183-91.
[136] Pandit MK, Sheikh AH, Singh BN. An improved higher order zigzag theory for the static
analysis of laminated sandwich plate with soft core. Finite Elements in Analysis and Design.
2008;44:602-10.
[137] Toledano A, Murakami H. A high-order laminated plate theory with improved in-plane
responses. International Journal of Solids and Structures. 1987;23:111-31.
[138] Rao KM, Meyer-Piening HR. Analysis of thick laminated anisotropic composite plates by
the finite element method. Composite Structures. 1990;15:185-213.
[139] Carrera E. Sullo sviluppo di un elemento finito computazionalmente efficiente e capace
di effetti HSDT. Aerotecnica Missili Spazio. 1994;73:14-23.
[140] Savithri S, Varadan T. Accurate bending analysis of laminated orthotropic plates. AIAA
journal. 1990;28:1842-4.
[141] Hsu T, Wang JT. A theory of laminated cylindrical shells consisting of layers of
orthotropic laminate. AIAA journal. 1970; 8:2141-6.
[142] Cho KN, Bert CW, Striz AG. Free vibrations of laminated rectangular plates analyzed by
higher order individual-layer theory. Journal of Sound and Vibration. 1991;145:429-42.
[143] NOSIER A, KAPANIA R, REDDY J. Free vibration analysis of laminated plates using a
layerwise theory. AIAA journal. 1993;31:2335-46.
[144] Reddy JN. On the Generalization of Displacement-Based Laminate Theories. Applied
Mechanics Reviews. 1989;42:S213-S22.
[145] Sun CT, Whitney J. Theories for the dynamic response of laminated plates. AIAA
journal. 1973;11:178-83.
[146] Chaudhuri RA. An equilibrium method for prediction of transverse shear stresses in a
thick laminated plate. Computers & Structures. 1986;23:139-46.
[147] Chaudhuri RA, Seide P. An approximate semi-analytical method for prediction of
interlaminar shear stresses in an arbitrarily laminated thick plate. Computers & Structures.
1987;25:627-36.
[148] Owen D, Li Z. A refined analysis of laminated plates by finite element displacement
methods--I. Fundamentals and static analysis. Computers & Structures. 1987;26:907-14.
[149] Reddy J, Barbero E, Teply J. A plate bending element based on a generalized laminate
plate theory. International Journal for Numerical Methods in Engineering. 1989;28:2275-92.
[150] Hinrichsen R, Palazotto A. Nonlinear finite element analysis of thick composite plates
using cubic spline functions. AIAA journal. 1986;24:1836-42.
[151] Epstein M, Glockner PG. Nonlinear analysis of multilayered shells. International
Journal of Solids and Structures. 1977;13:1081-9.
[152] Epstein M, Huttelmaier H. A finite element formulation for multilayered and thick
plates. Computers & Structures. 1983;16:645-50.
[153] Barker R, Pryor Jr C. A finite-element analysis including transverse shear effects for
applications to laminated plates. AIAA journal. 1971;9:912-7.
[154] Di Sciuva M, Gherlone M. A global/local third-order Hermitian displacement field with
damaged interfaces and transverse extensibility: analytical formulation. Composite
Structures. 2003;59:419-31.
[155] Cho M, Kim JS. Four-noded finite element post-process method using a displacement
field of higher-order laminated composite plate theory. Computers & Structures. 1996;61:283-
90.
170
Références bibliographiques
[156] Chakrabarti A, Sheikh AH. A new triangular element to model inter‐laminar shear
stress continuous plate theory. International Journal for Numerical Methods in Engineering.
2004;60:1237-57.
[157] Chakrabarti A, Sheikh AH. Analysis of laminated sandwich plates based on interlaminar
shear stress continuous plate theory. Journal of Engineering Mechanics. 2005;131:377-84.
[158] Kapuria S. An efficient coupled theory for multilayered beams with embedded
piezoelectric sensory and active layers. International Journal of Solids and Structures.
2001;38:9179-99.
[159] Kapuria S, Dumir P, Ahmed A. An efficient coupled layerwise theory for static analysis
of piezoelectric sandwich beams. Archive of Applied Mechanics. 2003;73:147-59.
[160] Kapuria S, Bhattacharyya M, Kumar A. Assessment of coupled 1D models for hybrid
piezoelectric layered functionally graded beams. Composite Structures. 2006;72:455-68.
[161] Kapuria S, Achary GGS. A coupled zigzag theory for the dynamics of piezoelectric hybrid
cross-ply plates. Archive of Applied Mechanics. 2005;75:42-57.
[162] Kapuria S, Achary GGS. Nonlinear coupled zigzag theory for buckling of hybrid
piezoelectric plates. Composite Structures. 2006;74:253-64.
[163] Kapuria S, Achary G. Nonlinear zigzag theory for electrothermomechanical buckling of
piezoelectric composite and sandwich plates. Acta Mechanica. 2006;184:61-76.
[164] Kapuria S, Dumir PC, Jain NK. Assessment of zigzag theory for static loading, buckling,
free and forced response of composite and sandwich beams. Composite Structures.
2004;64:317-27.
[165] Kapuria S, Dumir PC, Ahmed A. An efficient higher order zigzag theory for composite
and sandwich beams subjected to thermal loading. International Journal of Solids and
Structures. 2003;40:6613-31.
[166] Kapuria S, Achary GGS. An efficient higher order zigzag theory for laminated plates
subjected to thermal loading. International Journal of Solids and Structures. 2004;41:4661-84.
[167] Icardi U. Eight-noded zig-zag element for deflection and stress analysis of plates with
general lay-up. Composites Part B: Engineering. 1998;29:425-41.
[168] Icardi U. C° plate element for global/local analysis of multilayered composites, based on a
3D zig-zag model and strain energy updating. International Journal of Mechanical Sciences.
2005;47:1561-94.
[169] Li X, Liu D. Zigzag theory for composite laminates. AIAA journal. 1995;33:1163-4.
[170] Lu X, Liu D. An interlaminar shear stress continuity theory for both thin and thick
composite laminates. Journal of Applied Mechanics. 1992;59:502-9.
[171] Li X, Liu D. A laminate theory based on global–local superposition. Communications in
Numerical Methods in Engineering. 1995;11:633-41.
[172] Li X, Liu D. Generalized laminate theories based on double superposition hypothesis.
International Journal for Numerical Methods in Engineering. 1997;40:1197-212.
[173] Zhen W, Wanji C. Free vibration of laminated composite and sandwich plates using
global–local higher-order theory. Journal of Sound and Vibration. 2006;298:333-49.
[174] Zhen W, Wanji C. A study of global–local higher-order theories for laminated composite
plates. Composite Structures. 2007;79:44-54.
[175] Zhen W, Wanji C. A higher-order displacement model for stress concentration problems
in general lamination configurations. Materials & Design. 2009;30:1458-67.
[176] Chen W, Cheung YK, Wu Z. Augmented higher order global–local theory and refined
triangular element for laminated composite plates. Composite Structures. 2007;81:341-52.
[177] Wu Z, Chen R, Chen W. Refined laminated composite plate element based on global–
local higher-order shear deformation theory. Composite Structures. 2005;70:135-52.
[178] Chen W. Nine-parameter triangular thin plate bending element by using refined directed
stiffness method. Journal of Dalian University of Technology. 1993;33:289-95.
171
Références bibliographiques
[179] Cheung YK, Wanji C. Refined nine-parameter triangular thin plate bending element by
using refined direct stiffness method. International Journal for Numerical Methods in
Engineering. 1995;38:283-98.
[180] Wanji C, Cheung YK. Refined non-conforming quadrilateral thin plate bending element.
International Journal for Numerical Methods in Engineering. 1997;40:3919-35.
[181] Chen W, Cheung YK. The non-conforming nine-parameter triangular element for thin
plate bending and vibration analysis. Proceedings of the International Conference on CMSGE.
Hong Kong 1994. p. 1580–3.
[182] Xiaohui R, Wanji C, Zhen W. A new zig-zag theory and C° plate bending element for
composite and sandwich plates. Archive of Applied Mechanics. 2011;81:185-97.
[183] Batoz J-L, Tahar MB. Evaluation of a new quadrilateral thin plate bending element.
International Journal for Numerical Methods in Engineering. 1982;18:1655-77.
[184] Carrera E, Demasi L. Classical and advanced multilayered plate elements based upon
PVD and RMVT. Part 1: Derivation of finite element matrices. International Journal for
Numerical Methods in Engineering. 2002;55:191-231.
[185] Carrera E, Demasi L. Classical and advanced multilayered plate elements based upon
PVD and RMVT. Part 2: Numerical implementations. International Journal for Numerical
Methods in Engineering. 2002;55:253-91.
[186] Sedira L, Ayad R, Hecini M. Amélioration du modèle linéaire DDM par introduction
d'une fonction Zig-Zag. Application aux plaques stratifiées. Séminaire Internationale sur les
Technologies Mécaniques, Tlemcen (Algérie)05-06 Décembre 2009.
[187] Sedira L, Ayad R, Sabhi H, Hecini M, Sakami S. An enhanced discrete Mindlin finite
element model using a zigzag function. European Journal of Computational Mechanics/Revue
Européenne de Mécanique Numérique. 2012;21:122-40.
[188] Dhatt G, Batoz JL. Modélisation des structures par éléments finis: poutres et
plaques1990.
[189] Ali JSM, Bhaskar K, Varadan TK. A new theory for accurate thermal/mechanical
flexural analysis of symmetric laminated plates. Composite Structures. 1999;45:227-32.
[190] Batoz JL, Dhatt G. Modélisation des structures par éléments finis: solides élastiques:
Presses Université Laval; 1990.
[191] Irons B, Loikkanen M. An engineers' defence of the patch test. International Journal for
Numerical Methods in Engineering. 1983;19:1391-401.
[192] Morley LSD. Skew plates and structures: Pergamon Press;[distributed in the Western
Hemisphere by Macmillan, New York]; 1963.
[193] Babuška I, Scapolla T. Benchmark computation and performance evaluation for a
rhombic plate bending problem. International Journal for Numerical Methods in Engineering.
1989;28:155-79.
[194] Pagano N, Hatfield HJ. Elastic behavior of multilayered bidirectional composites. AIAA
journal. 1972;10:931-3.
[195] Engblom JJ, Ochoa OO. Through-the-thickness stress predictions for laminated plates of
advanced composite materials. International Journal for Numerical Methods in Engineering.
1985;21:1759-76.
[196] Srinivas S, Joga Rao CV, Rao AK. An exact analysis for vibration of simply-supported
homogeneous and laminated thick rectangular plates. Journal of Sound and Vibration.
1970;12:187-99.
[197] Owen DRJ, Figueiras JA. Anisotropic elasto-plastic finite element analysis of thick and
thin plates and shells. International Journal for Numerical Methods in Engineering.
1983;19:541-66.
[198] Fares ME, Elmarghany MK. A refined zigzag nonlinear first-order shear deformation
theory of composite laminated plates. Composite Structures. 2008;82:71-83.
172
Références bibliographiques
[199] Xiao JR, Gilhooley DF, Batra RC, Gillespie Jr JW, McCarthy MA. Analysis of thick
composite laminates using a higher-order shear and normal deformable plate theory
(HOSNDPT) and a meshless method. Composites Part B: Engineering. 2008;39:414-27.
[200] Plantema FJ. Sandwich construction: Wiley New York,, USA; 1966.
[201] J. AJ. Bending theory for multilayered orthotopic sandwich plates. AIAA journal.
1968;6:2166-9.
[202] Reddy JN, Chao WC. A comparison of closed-form and finite-element solutions of thick
laminated anisotropic rectangular plates. Nuclear Engineering and Design. 1981;64:153-67.
[203] Khatua TP, Cheung YK. Bending and vibration of multilayer sandwich beams and
plates. International Journal for Numerical Methods in Engineering. 1973;6:11-24.
[204] Topdar P, Chakraborti A, Sheikh AH. An efficient hybrid plate model for analysis and
control of smart sandwich laminates. Computer Methods in Applied Mechanics and
Engineering. 2004;193:4591-610.
[205] Reddy JN. Free vibration of antisymmetric, angle-ply laminated plates including
transverse shear deformation by the finite element method. Journal of Sound and Vibration.
1979;66:565-76.
[206] Levinson M. Free vibrations of a simply supported, rectangular plate: An exact elasticity
solution. Journal of Sound and Vibration. 1985;98:289-98.
[207] Aboura Z. Etude du processus de délaminage mode I, mode II et mode mixte (I+II) de
matériaux composites à renforts tissé à différentes vitesses de sollicitation [Doctorat Thesis]:
University of Technology of Compiègne; 1993.
[208] Talbi N, Ayad R. An efficient shell finite element based upon discrete Mindlin
hypothesises for non linear structural analysis. Proceedings of the Fifth World Congress on
Computational Mechanics. Vienne- Autriche.
[209] Hendou H. Analyse statique et dynamique des coques composites à l’aide du modèle
MiSP. Application au carton ondulé. UTC Compiègne; 2002.
[210] Uflyand YS. The propagation of waves in the transverse vibrations of bars and plates.
Prikl Mat Mekh. 1948;12:287-300.
[211] Kulkarni S, Kapuria S. Free vibration analysis of composite and sandwich plates using
an improved discrete Kirchhoff quadrilateral element based on third-order zigzag theory.
Computational Mechanics. 2008;42:803-24.
[212] Jeyachandrabose C, Kirkhope J, Meekisho L. An improved discrete Kirchhoff
quadrilateral thin-plate bending element. International Journal for Numerical Methods in
Engineering. 1987;24:635-54.
[213] Kapuria S, Kulkarni S. An improved discrete Kirchhoff quadrilateral element based on
third-order zigzag theory for static analysis of composite and sandwich plates. International
Journal for Numerical Methods in Engineering. 2007;69:1948-81.
[214] Kärger L, Wetzel A, Rolfes R, Rohwer K. A three-layered sandwich element with
improved transverse shear stiffness and stresses based on FSDT. Computers & Structures.
2006;84:843-54.
[215] Kim J-S, Cho M. Buckling analysis for delaminated composites using plate bending
elements based on higher-order zig-zag theory. International Journal for Numerical Methods
in Engineering. 2002;55:1323-43.
[216] Kulkarni S, Kapuria S. A new discrete Kirchhoff quadrilateral element based on the
third-order theory for composite plates. Computational Mechanics. 2007;39:237-46.
[217] ABAQUS/STANDARD User’s Manual2002.
[218] BATOZ JL, Dhatt G. Modélisation des structures par éléments finis: Coques. Paris:
Hermès; 1992.
[219] Methods NAfFE, Standards. The standard NAFEMS benchmarks: NAFEMS; 1990.
[220] Flügge W. Stresses in shells: Springer Verla; 1960
173
Références bibliographiques
[221] Lindberg G, Olson M, Cowper G. New developments in the finite element analysis of
shells. Q Bull Div Mech Eng and Nat Aeronautical Establishment, National Research council
of Canada. 1969;4.
[222] CHO M, ROH HEEY. Development of geometrically exact new shell elements based on
general curvilinear co-ordinates. International Journal for Numerical Methods in Engineering.
2003;56:81-115.
[223] Macneal RH, Harder RL. A proposed standard set of problems to test finite element
accuracy. Finite Elements in Analysis and Design. 1985;1:3-20.
[224] Batoz JL, Triki S. Développement d'un élément fini de poutre vrillée, rapport interne.
Octobre 1991. p. UTC/LG2MS/MNM.
[225] Simo J, Fox D, Rifai M. On a stress resultant geometrically exact shell model. Part II:
the linear theory; computational aspects. Computer Methods in Applied Mechanics and
Engineering. 1989;73:53-92.
[226] Brebbia C, Connor J. Fundamentals of finite element techniques for structural
engineers. Wiley (New York); 1974.
[227] Brebbia C, Connor J. Stiffness matrix for shallow rectangular shell element(Stiffness
matrix for shallow rectangular shell element and application of finite element displacement
method). AMERICAN SOCIETY OF CIVIL ENGINEERS, ENGINEERING MECHANICS
DIVISION, JOURNAL. 1967;93:43-65.
[228] Kraus H. Thin elastic shells: John Wiley & Sons; 1967.
[229] Chatelain J. Analyse non linéaire des coques minces isotropes et composites par
éléments finis quadrilatéraux. Université de Technologie de Compiègne; 1987.
[230] Paccola RR. Análise não linear física de placas e cascas anisotrópicas laminadas
acopladas ou não com meio contínuo tridimensional viscoelástico através da combinação entre
o MEC eo MEF (2004). Physical non-linear analysis of anisotropic laminated plates and shells
coupled or not with three-dimensional viscoelastic medium by BEM/FEM coupling.
[231] Varadan T, Bhaskar K. Bending of laminated orthotropic cylindrical shells—an elasticity
approach. Composite Structures. 1991;17:141-56.
[232] Cho M, Kim JS. A postprocess method for laminated shells with a doubly curved nine-
noded finite element. Composites Part B: Engineering. 2000;31:65-74.
[233] Ren J. Analysis of simply-supported laminated circular cylindrical shell roofs. Composite
Structures. 1989;11:277-92.
[234] Crawley EF. The natural modes of graphite/epoxy cantilever plates and shells. Journal of
Composite Materials. 1979;13:195-205.
[235] Frey F. L'analyse statique non linéaire des structures par la méthode des éléments finis
et son application à la construction métallique: Université de Liège, Faculté des sciences
appliquèes, Laboratoire de mećanique des matériaux et statique des constructions; 1978.
[236] Argyris JH, Dunne PC, Malejannakis G, Schelkle E. A simple triangular facet shell
element with applications to linear and non-linear equilibrium and elastic stability problems.
Computer Methods in Applied Mechanics and Engineering. 1977;10:371-403.
[237] HAMMADI F. Formulation et évaluation d'éléments finis à continuité C° de la géométrie
pour l'analyse linéaire et non linéaire des coques [Thèse de Doctorat]: Université de
Compiègne; 1998.
[238] Timoshenko S, Gere JM. Mechanics of materials. New York: Van Nostrand Reinhold Co.;
1972.
[239] Fafard M. Eléments finis en mécanique des solides. Cours GCI-63244. Université de
Laval, Québec, Canada1993.
[240] Sze K, Liu X, Lo S. Popular benchmark problems for geometric nonlinear analysis of
shells. Finite Elements in Analysis and Design. 2004;40:1551-69.
[241] Batoz JL, Cantin G, CA. NPSM. Geometrically Nonlinear Analysis of Shell Structures
Using Flat DKT Shell Elements: Defense Technical Information Center; 1985.
174
Références bibliographiques
[242] Surana KS. Geometrically nonlinear formulation for the curved shell elements.
International Journal for Numerical Methods in Engineering. 1983;19:581-615.
[243] Ammar S, Dhatt G, Fafard M. Exact stability model of space frames. Computers &
Structures. 1996;60:59-71.
[244] Pol P. Modélisation du comportement élastoplastique de coques minces par éléments
finis. Université de Technologie de Compiègne; 1992.
[245] Brendel B, Ramm E. Linear and nonlinear stability analysis of cylindrical shells.
Computers & Structures. 1980;12:549-58.
[246] Batoz J-L, Dhatt G. Plaques et coques par éléments finis: analyse linéaire et non
linéaire:[conférences] 2-4 juin 1986: Institut pour la Promotion des Sciences de l'Ingénieur;
1986.
[247] BATOZ J, JAAMEI S. Etude de différentes formulations lagrangiennes pour l'analyse
des poutres et des coques minces élastiques en grandes rotations.
[248] JAAMEI S. Etude de différentes formulations lagrangiennes pour l'analyse non linéaire
des plaques et des coques minces élasto-plastiques en grandes rotations [PhD]: Universitie de
Technologie de Compiègne; 1986.
[249] Fafard M. Calcul automatique des configurations pré et post-flambement en calcul non
linéaire des structures [Thèse de doctorat]. Ottawa: Université Laval, Québec, Canada; 1987.
[250] Quesnel T. Facteurs de correction de cisaillement transverse dans les plaques composites
[Mémoire d'Etudes approfondies]: Université de P. et M. Curie, Université de Technologie de
compiègne; 1990.
175
Annexe A
(A.1)
avec
et sont les contraintes planes et de cisaillement, définies par
et (A.2)
(A.3)
; (A.4a,b)
; (A.5a,b)
; (A.6 a,b)
, (A.7)
avec
176
Annexe A
, , , , (A.8a-e)
avec
, (A.11a,b)
Des calculs primitifs ont montré qu’il n y’a pas une méthode prédominante sur
le résultat obtenu. En effet, la méthode totale gagne beaucoup du terrain du fait
qu’elle s’adapte d’autant plus pour différents types des stratifiés ou des sandwichs,
que ce soit, symétriques ou non symétriques. Ces coefficients corrigent
177
Annexe A
178
Annexe A
Début
‐ ‐ ‐
2 NC
‐
‐
‐ ‐
Ecrire et ‐ ‐ ‐
-
- - - -
Fin
179
Annexe A
Début
Calcul de la matrice
‐ ‐ ‐
NC
Calcul de la matrice de la
couche
= Calcul de la matrice
=
‐ ‐
‐
Ecrire et ‐ ‐ ‐
‐
-
- - - -
Fin
180
Annexe B
(B.1)
(B.2)
(B.3)
(B.4)
(B.5)
(B.6)
(B.7)
(B.8)
(B.9)
(B.10)
avec
(B.11)
(B.12)
, (B.13)
(B.14)
181