Thse BELARBI
Thse BELARBI
Thse BELARBI
net/publication/290607344
CITATIONS READS
2 1,197
1 author:
Mohamed-Ouejdi Belarbi
Université de Biskra
34 PUBLICATIONS 124 CITATIONS
SEE PROFILE
Some of the authors of this publication are also working on these related projects:
Layerwise finite element model for composite sandwich plates. View project
Hygrothermal based analysis of sandwich FGM plates, beams and shells View project
All content following this page was uploaded by Mohamed-Ouejdi Belarbi on 17 January 2016.
Présentée par :
Mohamed Ouejdi BELARBI
i
i
DEDICACE
Je dédie ce travail à la
mémoire de mon père, mon exemple, bien aimé, qu’Allah ait son âme pour tous
ii
REMERCIEMENTS
Tout d’abord, j'exprime mes vifs remerciements à mon directeur de thèse, Monsieur
Abdelouahab TATI, Professeur à l'Université de Biskra, qui m'a encadré durant la période de
la thèse avec grande patience et bonne humeur. Son assistance, sa disponibilité, les
encouragements qu'il n'a cessé de me prodiguer et enfin la confiance qu'il m’a toujours
témoigné, m'ont été d'un grand secours pour l'accomplissement de ce travail. Je lui en suis très
reconnaissant.
Monsieur Fodil HAMADI, Professeur à l'Université de Bechar, qui a accepté lui aussi
d'examiner cette thèse. Je tiens à le remercier pour son extrême gentillesse.
Je tiens aussi à remercier en particulier mon collègue de travail: H. Ounis pour son
soutien moral et scientifique durant notre parcours doctoral.
iii
Publications
En plus de cette thèse de doctorat, ces travaux ont conduit aux publications et
communications suivantes:
Publications internationales
1. M. O. Belarbi, A. Tati, (In press), ''A new C0 finite element model for the analysis of
sandwich plates using combined theories'' International Journal of Structural Engineering.
4. A. Tati, M. O. Belarbi, '' A C0 rectangular element model for sandwich plates analysis
using combined theories " Submitted to Measurement journal, 2014.
Communications internationales
7. M. O. Belarbi, A. Tati, ''Étude du flambement des poutres sandwich par éléments finis''
Congrès Algérien de Mécanique (CAM), 25-28 novembre 2013, Mascara-Algérie.
iv
8. M. O. Belarbi, A. Tati, A. Khechai, ''Efficient layerwise finite element model for
multilayer sandwich plates analysis'' Sixth International Congress on Design and
Modeling of Mechanical Systems (CMSM’2015), Tunisia, March 23-25, 2015, Springer.
pp 305-314, DOI: 10.1007/978-3-319-17527-0_30
Communications nationals
10. M. O. Belarbi, A. Tati, H. Ounis, ''Finite Element for the Analysis of Sandwich Plates in
Bending'' JEMGM II, 22 et 23 Avril 2014, Biskra-Algeria.
11. H. Ounis, A. Tati, M. O. Belarbi, '' Thermal buckling behavior of laminated plates with
cut-outs'' JEMGM II, 22 et 23 Avril 2014, Biskra-Algeria.
12. M. O. Belarbi, A. Tati, H. Ounis, ''Étude numérique du flambage des poutres isotropes
et des poutres anisotropes'' JEGM 29 avril 2013, Biskra-Algérie.
13. M. O. Belarbi, A. Tati, ''Étude du Flambement des Poutres Sandwichs par MEF'
JEMGM’ 2013, Biskra-Algérie.
v
Résumé
Dans l’analyse des structures composites et notamment les structures sandwichs, la
détermination des déplacements et des contraintes interlaminaires aux interfaces requiert des
théories et des modèles d’éléments finis raffinés, capables de donner une meilleure
description aux réponses locales. Pour des raisons de coûts élevés en terme de mémoire et de
temps de calcul dans le cas des approches tridimensionnelles, les approches bidimensionnelles
de type layerwise ont été développées et ont été bien appliquées pour l’analyse des structures
composites sandwichs. Cependant, ces approches requièrent un temps de calcul relativement
long lorsqu’il s’agit d’un grand nombre de couches.
Dans ce travail de thèse, un nouveau modèle de type layerwise a été développé pour
l’analyse du comportement des plaques composites multicouches et notamment les plaques
sandwichs. Contrairement à tous les autres modèles layerwise présentés dans la littérature, le
nombre de degrés de liberté dans le présent modèle est indépendant du nombre de couches. Le
modèle proposé est issu d'une combinaison de deux théories, où la théorie de déformation en
cisaillement du premier ordre est adoptée pour modéliser les deux peaux, tandis que l’âme est
modélisée en utilisant une théorie de déformation en cisaillement du troisième ordre. Ces
deux théories combinées, satisfont les exigences de continuité des déplacements
interlaminaires ainsi que les conditions de compatibilité. Le nouveau modèle permet de
calculer les contraintes de cisaillement transversal d’une manière précise sans recours à des
facteurs de correction.
Sur la base du modèle proposé, deux éléments finis ont été développés pour l’analyse
statique et dynamique des plaques composites stratifiées et sandwichs. Ces éléments sont
géométriquement simples et ne possèdent que quatre nœuds et treize degrés de liberté par
nœud. La comparaison des résultats obtenus avec des solutions de références, déterminés
analytiquement, et ceux obtenus par d’autres modèles d’éléments finis disponibles dans la
littérature, a montré la performance et la précision du modèle proposé. Par ailleurs, une étude
paramétrique a également été menée pour mettre en évidence l'effet de certains paramètres,
sur les comportements de flexion et de vibration des plaques sandwichs multicouches.
Mots clés:
Composites multicouches, Sandwichs, Élément fini, Layerwise, Analyse statique, Dynamique
i
Abstract
In the analysis of composite structures, especially sandwich structures, the
determination of displacements and interlaminar stresses at the interfaces, requires theories
and finite elements models, able to give a better description of the local response. For reasons
of high cost in terms of memory and computation time in the case of three-dimensional
approaches, the “layerwise” two-dimensional approaches have been developed and have been
properly applied for the analysis of multilayer composites structures. However, these
approaches require a high computational time when dealing with a large number of layers.
In this thesis, a new layerwise model has been developed for the analysis of the
behavior of laminated composite plates and especially sandwich plates. Unlike all other
layerwise models, presented in the literature, the number of degrees of freedom in the present
model is independent of the number of layers. The proposed model is derived from a
combination of two theories; where the first-order shear deformation theory is adopted to
model the two face sheets, while the core is modeled using the third-order shear deformation
theory. This combined theory satisfies the continuity requirements of interlaminar
displacements as well as the compatibility conditions. The new model allows calculating the
transverse shear stresses in a precise manner without the use of correction factors.
Based on the proposed model, two finite elements have been developed for static and
dynamic analysis of laminated composite and sandwich plates. These elements are
geometrically simple and have only four nodes and thirteen degrees of freedom per node. The
comparison of the results obtained with those given by the three-dimensional elasticity
solutions and those obtained by other finite element models available in the literature; showed
the performance and accuracy of the proposed model. Moreover, a parametric study was also
conducted to investigate the effect of certain characteristic parameters on the bending and free
vibration behavior of laminated sandwich plates.
Keywords:
ii
الملخص
في تحليل الهياكل المركبة وخاصة هياكل السندويش ،تحديد االنتقاالت و االجهادات الموجودة ما
بين الطبقات ،يتطلب نظريات ونماذج من العناصر المحدودة المطورة ،القادرة على إعطاء وصف أفضل
لالستجابات المحلية .نظرا للتكلفة العالية من حيث الذاكرة وحساب الوقت في حالة نهوج ثالثية األبعاد،
فان نهوج ثنائية األبعاد من نوع اليرويز تم تطويرها و تطبيقها بشكل صحيح لتحليل الهياكل المركبة
متعددة الطبقات .لكن يجدر بنا ان نذكر ان هذه النهوج تتطلب حساب طويل نسبيا في حالة وجود عدد كبير
من الطبقات.
في هذه األطروحة ،تم تطوير نموذج جديد من نوع اليرويز من اجل تحليل سلوك الصفائح المركبة
و خاصة صفائح السندويش .خالفا لجميع نماذج اليرويز األخرى ،المقدمة في المراجع ،فان عدد درجات
الحرية في هذا النموذج هو مستقل عن عدد الطبقات .النموذج المقترح هو ناتج من مزيج نظريتين؛ نظرية
تشوه القص من الدرجة األولى التي اعتمدت لنمذجة اثنين من الطبقات القشرية ،في حين ان نظرية تشوه
القص من الدرجة الثالثة اعتمدت لنمذجة الطبقة الوسطى .هاتين النظريتين تلبيان متطلبات استمرارية
االنتقاالت بين الصفائح وكذلك شروط التوافق .النموذج الجديد يسمح بحساب إجهادات القص بطريقة
دقيقة من دون استخدام عوامل التصحيح.
بناءا على النموذج المقترح ،تم تطوير عنصرين محدودين من اجل تحليل السلوك الستاتيكي
والديناميكي للصفائح المركبة و صفائح السندويش .هذان العنصران يتميزان ببساطة الشكل الهندسي كما
انهما ال يحتويان اال على اربعة عقد و ثالثة عشر درجة حرية لكل عقدة .مقارنة النتائج المتحصل عليها
مع الحلول التحليلية وتلك التي تم الحصول عليها من قبل نماذج العناصر المحدودة ،المتاحة في المراجع،
أظهرت االداء الجيد و اثبتت دقة النموذج المقترح .من ناحية أخرى ،أجريت دراسة بارومتريه إلظهار
تأثير بعض المعلمات على سلوك االنحناء واالهتزاز لصفحات الساندويتش متعدد الطبقات.
الكلمات المفتاحية:
الهياكل المركبة من عدة طبقات ،سندويش ،العناصر المحدودة ،اليرويز ،السلوك الستاتيكي ،ديناميكي
iii
Table des matières
Résumé ……………………………………………………………………………... i
Abstract ………………...…………………………….……………………..………. ii
…………………………………………………………………………… ملخص..…. iii
Table des matières ……………………………………………………………….…. iv
Liste des figures ………………………………………………...……………...…… ix
Liste des tableaux ……………………………………………………….………….. xiv
Introduction générale ………………………………………………..……..……….. 1
1. Introduction ……………………………………………………………………… 1
2. Objective de la thèse …………………………………………….....……………. 3
3. Plan de la thèse …………………………………………………..…..………..… 3
Chapitre 1: Aperçu sur l’évolution des théories aux éléments finis pour la
modélisation des composites sandwichs
iv
1.2.2.2.1 Théories Zig-Zag …………………….……………… 26
1.2.2.2.2 Théories Global-Local d'ordre supérieur (GLHSDT) ... 30
1.2.2.3 Modèles à cinématique enrichie ……………………….………. 33
1.3 Modèles éléments finis des structures sandwichs multicouches …………..….. 35
1.3.1 Éléments finis basés sur les théories du premier ordre ………………..... 35
1.3.2 Éléments finis basés sur les théories d'ordre supérieur …………...…….. 37
1.3.3 Éléments finis basés sur l’approche Layerwise ………………………… 40
1.3.4 Éléments finis basés sur la théorie de Zig- Zag ……………………..….. 46
1.4 Conclusion ………………...…..…………………………….………………… 49
v
Chapitre 3 Développement d’un nouvel élément fini rectangulaire (RSFT52)
pour l’analyse des plaques sandwichs et multicouches
vi
Chapitre 4 Développement d’un nouvel élément fini isoparamétrique
(QSFT52) en introduisant la technique de déformation naturelle
vii
PARTIE III APPLICATION D'ELEMENT QSFT52 POUR L’ETUDE
DYNAMIQUE DES PLAQUES SANDWICHS ET DES STRATIFIEES
viii
Liste des figures
CHAPITRE 1
Figure 1.4 : Géométries déformées et non déformées d'un bord d'une plaque sous
l'hypothèse de Kirchhoff ……………………………………..…….. 10
Figure 1.5 : Géométries déformées et non déformées d'un bord d'une plaque
sous l'hypothèse de Reissner-Mindlin ………………………….…... 12
Figure 1.7 : Géométries de non déformées et déformées d'une plaque selon les
trois théories CLPT, FSDT et HSDT ……………………………..… 15
Figure 1.10: Géométrie et notations employées dans les théories en Zig-Zag …... 27
Figure 1.11 : Champ de déplacement des modèles Zig-Zag du premier ordre ….... 28
Figure 1.12 : Champ de déplacement des modèles Zig-Zag d’ordre supérieur ...… 28
ix
Figure 1.13 : Système des coordonnées, ordre de couche d’une plaque
multicouche …..................................................................................... 31
Figure 1.14 : Élément triangulaire à 16 nœuds et leurs degrés de liberté [215] …... 37
CHAPITRE 2
x
CHAPITRE 3
Figure 3.3 : Plaque stratifiée carrée à trois couches, simplement appuyée sous
chargement doublement sinusoïdal …………………...…………….. 79
xi
Figure 3.13 : Effet du rapport d’épaisseur a h sur le déplacement transversal
( w) d’une plaque sandwich (0 90 0 C 0 90 0) encastré et
simplement appuyée ………………………………………………... 103
CHAPITRE 4
Figure 4.1 : Transformation géométrique de l’élément réel vers l’élément parent. 108
Figure 4.4 : Plaque biaise simplement appuyée avec un maillage de m×n …..….. 120
xii
CHAPITRE 5
Figure 5.2 : Plaque sandwich carrée simplement appuyée avec des peaux
stratifiées 0 / 90 / C / 0 / 90 pour différents rapport d’épaisseur -
Comparaison des fréquences naturelles …………………………….. 162
Figure 5.3 : Plaque sandwich carrée simplement appuyée avec des peaux
stratifiées 0 / 90 / C / 0 / 90 pour différents rapport (hc h f ) -
Comparaison des fréquences naturelles ……….……………………. 162
Figure 5.4 : Plaque sandwich simplement appuyée ayant des peaux stratifiées
0 / 90 / C / 0 / 90 avec différents rapport d’aspect (a b) -
Comparaison des fréquences naturelles ………………..…………… 163
xiii
Liste des tableaux
CHAPITRE 1
CHAPITRE 3
Tableau 3.1 : Détails de conditions aux limites utilisées pour les exemples
numériques ………………………………………………………... 76
Tableau 3.4 : Propriétés des matériaux (normalisés) utilisés pour les plaques
stratifiées …………….………………..………..…………………. 79
xiv
couches (0 / 90 / 90 / 90 / 0) , simplement appuyée, sous chargement
uniforme …………..……………….………………………………
Tableau 3.9 : Propriétés des matériaux (normalisés) utilisés pour les plaques
sandwichs ……………………………………………………………… 84
Tableau 3.18: Déflexion ( w) d’une plaque sandwich rectangulaire ayant des 100
peaux stratifiées, sous chargement réparti uniforme …………...….
xv
CHAPITRE 4
Tableau 4.2 :
Déplacement transversal ( w) , contraintes planes xx , yy , xy et
contraintes de cisaillement transversal xz , yz , d’une plaque
sandwich symétrique ( f / c / f ) simplement appuyée sous
chargement doublement sinusoïdal ……………………………….. 124
Tableau 4.3 :
Déplacement transversal ( w) , contraintes planes xx , yy , xy et
contraintes de cisaillement transversal xz , yz , d’une plaque
sandwich symétrique ( f / c / f ) encastrée sous chargement
uniforme …………………………………………………………... 130
Tableau 4.4 : Propriétés mécaniques utilisés dans l’exemple 4.4.4 ……….…….. 131
Tableau 4.5 :
Déplacement transversal ( w) , contraintes planes xx , yy , xy et
contraintes de cisaillement transversal xz , yz , d’une plaque
sandwich symétrique ( f / c / f ) sous chargement uniforme …...…. 133
xvi
CHAPITRE 5
Tableau 5.4 : Propriétés mécaniques des matériaux utilisés dans l’exemple 5.4.4 149
Tableau 5.8 : Propriétés mécaniques des matériaux utilisés dans l’exemple 5.4.5 156
xvii
Tableau 5.10: Fréquence propre non-dimensionnelle ( ) , pour différents modes,
d’une plaque sandwich carrée simplement appuyée ayant des
peaux composites stratifiées non-symétriques 0 / 90 / C / 0 / 90
avec différents rapports d’épaisseurs (a h) .…................................ 160
Tableau 5.13: Propriétés mécaniques des peaux et de l’âme (HEREX - C70.130) 163
xviii
Introduction générale
Introduction générale
1. Introduction
Au cours des dernières décennies, le développement rapide des technologies modernes
requiert que l'on utilise des matériaux ayant des caractéristiques spécifiques élevées (bons
rapports rigidité-poids et résistance-poids, faible coût d'entretien, excellente durabilité et bien
d’amples qualités). Les matériaux composites sont des matériaux qui répondent à ces
exigences. Ils deviennent de plus en plus incontournables dans de nombreuses composantes
structurales comme la marine, le transport, la navigation, les sports, les applications
aérospatiales et ils ont été utilisés récemment dans les structures de génie civil comme des
compléments aux matériaux classiques (acier, bois et béton). Ces derniers temps, les
ingénieurs civils ont exploité les avantages d'utilisation des matériaux composites et
spécialement les plastiques renforcés par des fibres de verre dans la réparation du béton.
Parmi les composites les plus convoités, on peut distinguer deux grandes catégories de
matériaux composites structuraux: les stratifiés et les sandwichs. Les stratifiés sont des
matériaux couramment rencontrés dans les éléments structuraux. Ces derniers sont constitués
par l’empilement de plusieurs couches où chaque couche peut avoir une orientation différente,
i.e. les caractéristiques mécaniques varient selon les directions considérées. Les sandwichs
sont des structures composites dont l'utilisation connaît un essor important. Ces structures se
composent de trois couches principales dont les caractéristiques mécaniques sont différentes:
deux peaux fines, généralement stratifiées, possédant de très bonnes caractéristiques
mécaniques (module d’élasticité élevé ainsi qu’une grande résistance à la compression et à la
traction), collées sur une âme épaisse, légère et de faibles caractéristiques mécaniques
(mousse, bois) dont le rôle principal est de résister aux efforts de cisaillement. Ces différentes
couches sont liées entre elles par collage ou par soudure.
Malgré les nombreux avantages que présentent ces matériaux, leur comportement est
assez complexe à cause des problèmes spécifiques liés principalement à leur hétérogénéité.
Parmi ces complexités liées à l’hétérogénéité, l’effet de cisaillement transverse est
1
Introduction générale
prédominant et joue un rôle essentiel, car il est considéré comme étant la cause de plusieurs
modes de ruptures comme le délaminage au niveau des interfaces, fissure transverse, etc. Cet
effet devient beaucoup plus complexe dans le cas des structures sandwichs car la variation des
propriétés mécaniques entre le cœur et les peaux est très importante.
Toutes les théories et les approches développées à ce jour pour l’analyse des éléments
de structures composites et notamment de sandwichs, présentent des limites de performances
ou des difficultés dans les moyens de calcul. On cite à titre d’exemple, les modèles
monocouche équivalente (approximation globale), dont les variables sont introduites pour
toute la plaque/coque. Par conséquent, ces modèles sont incapables de prédire avec précision
le comportement local des structures multicouches y compris les endommagements.
L’approche monocouche équivalente englobe essentiellement trois théories bien connues: la
théorie classique des stratifiés, basée sur les hypothèses de Kirchhoff, qui néglige l'effet de
cisaillement transverse, la théorie de déformation en cisaillement du premier ordre, basée sur
les hypothèses de Reisssner-Mindlin et les théories de déformation en cisaillement d’ordre
supérieur basées sur l’approximation non linéaire des déplacements.
Pour ces raisons, des modèles plus précis, basés sur une description par-couche,
(layerwise, i.e. approximation locale), ont été proposés. Cette description est plus répondue
car chaque couche est traitée individuellement, en imposant les conditions de continuité en
déplacements ou en contraintes aux interfaces. Cependant, ces modèles exigent énormément
de temps de calcul, notamment lorsque la plaque stratifiée est constituée d’un grand nombre
de couche, puisque le nombre d'inconnues augmente d’une manière proportionnelle avec
l'augmentation du nombre de couches.
La méthode des éléments finis s'est imposée dans les dernières années comme étant une
méthode puissante et efficace, ce qui lui permet d'être largement utilisée dans l’analyse de
comportement complexe des structures sandwichs composites. La plupart des modèles
éléments finis, basés sur les théories sus-citées, sont pratiquement difficile à mettre en œuvre.
2
Introduction générale
Par conséquent, les ingénieurs doivent choisir entre les modèles éléments finis
classiques moins coûteux mais moins précis au niveau local et ceux raffinés et sophistiqués
mais trop coûteux en terme de temps de calcul. Ces difficultés ont amené les chercheurs à
développer d’autres modèles d’éléments finis assurant un meilleur compromis entre la
précision et le faible coût. Toutefois, ce domaine requiert des recherches importantes et plus
approfondies sur les modèles d’éléments finis qui sont capables d’analyser d’une manière plus
précise le comportement de ces structures.
2. Objectif de la thèse
L'objectif principal de cette thèse est le développement d'un nouveau modèle sandwich,
basé sur l’approche layerwise mais en gardant tous les avantages de l'approche monocouche,
e.g. le nombre de degrés de liberté est indépendant du nombre de couches. Sur la base de ce
nouveau modèle, on procède à la création de deux éléments finis en assurant un meilleur
compromis entre la précision et le faible coût. Ces éléments développés seront destinés à
l’analyse des différents types de comportements (statique et dynamique) des structures en
matériaux isotropes, composites stratifiées et plus spécialement les structures sandwichs. Par
ailleurs, un autre objectif de ce travail est d'étudier l'influence de différents paramètres sur la
variation du déplacement transversal ainsi que les fréquences naturelles d’une plaque
sandwich, tels que le rapport d’épaisseur, les combinaisons de conditions aux limites, le
rapport d’aspect, le degré d'anisotropie des sandwichs, etc.
3. Plan de la thèse:
Le travail de cette thèse est divisé en trois grandes parties:
Partie II : Développement d'un nouveau modèle d’éléments finis sandwich sur la base de
l’approche layerwise, articulée autour de trois chapitres :
3
Introduction générale
Partie III : Application de l’élément QSFT52 pour l’étude dynamique des plaques
sandwichs et des stratifiées
Finalement, le travail se termine par une conclusion générale résumant les résultats
obtenus et par des perspectives.
4
PARTIE I
SYNTHESE BIBLIOGRAPHIQUE
5
Chapitre 1 Analyse bibliographique
Chapitre 1
1.1 Introduction
6
Chapitre 1 Analyse bibliographique
Les premières analyses des structures sandwichs qui sont rapportées dans la littérature,
ont débuté juste après la deuxième guerre mondiale. En 1948, Libove et Batdorf [10] ont
développé une théorie de petit déplacement ‘small deflection theory’ pour les plaques
sandwichs anisotropes. Ce travail a été effectué sous le comité consultatif national pour
l'aéronautique ‘the National Advisory Committee for Aeronautics (NACA)’. Dans la même
année, Reissner [11] a fourni des solutions théoriques aux problèmes à grands déplacements.
Il suppose que les deux peaux ayant une épaisseur très petite, reprennent uniquement les
efforts membranaires, tandis que l’âme ne subit que les contraintes de cisaillement
transverse. En 1950, Reissner [12] a utilisé le même principe de Libove [10] pour étudier les
coques sandwichs. En revanche, Hoff [13] et Eringen [14] dérivent l’équation différentielle
pour étudier la flexion et le flambement des plaques sandwichs, en utilisant le principe des
déplacements virtuels. En 1951, Yen et al. [15] présentent une étude pour les plaques
sandwichs simplement supportées soumises à deux types de chargements (uniformément
réparti et concentré). Ils ont intégré les équations différentielles pour obtenir le déplacement
transversal. Une autre étude a été faite par Gerard [16] pour les plaques sandwichs épaisses.
D’autre part, Ericksen [17] a calculé la flèche et les contraintes des panneaux sandwichs sous
chargement normal. En 1958, l’auteur [18] a publié un rapport pour analyser les structures
sandwichs non-symétriques avec des peaux d’épaisseur inégale. Dans les années soixante,
Cheng [19] a présenté les mêmes théories de Reissner, mais dans cette étude, le matériau
utilisé pour l’âme était orthotrope. Lin et Yokota [20] ont calculé le moment de flexion et la
flèche d’un panneau sandwich sous deux chargements combinés (compression biaxiale et
pression). En 1967, Liaw et Little [21] ont considéré une plaque sandwich composée de (n)
membranes isotropes ayant différentes épaisseurs et (n-1) cœurs. Cette étude a été reprise par
Azar [22] pour les matériaux orthotropes.
7
Chapitre 1 Analyse bibliographique
Toutes les théories et les approches développées à ce jour pour l’analyse des éléments
de structures conçus à base de ces matériaux, telle que les plaques et les coques minces ou
épaisses, présentent des limites de performances ou des difficultés dans les moyens de calcul.
Dans ce contexte, de nombreux travaux d’investigations ont été menés à travers les années,
dans l’objectif d’améliorer les différentes approches afin de prévoir les réponses statiques et
dynamiques ainsi que la stabilité de ces structures. Selon lesquelles on peut distinguer deux
modèles principaux: le modèle tridimensionnel (3D) et le modèle bidimensionnel (2D).
L’approche tridimensionnelle vise à donner des résultats plus précis, utilisés notamment
comme référence, mais peu d’auteurs ont adopté cette approche à cause de son coût en terme
de temps de calcul et en capacité de mémoire pour le stockage de la base des données requise
qui est plus importante que celle dans l'analyse bi-dimensionnelle, (Pagano [23, 24], Jones
[25], Srinivas et Rao [26, 27], Noor et al.[28, 29] et Kant et Swaminathan [30]). De plus, les
solutions 3D sont difficiles à obtenir dans la majorité des cas de géométrie et de conditions
aux limites. Parmi les solutions d’élasticité 3D, on trouve l’article de Pagano [24] qui a
proposé une solution d'élasticité tridimensionnelle, sous forme de série sinusoïdale, pour le
problème de flexion des plaques stratifiées rectangulaires. Une autre solution a été présentée
par Noor et al. [28] pour le flambement et la vibration libre de panneaux sandwichs.
À ce jour, l'utilisation des modèles bidimensionnels (2D) est préférable dans la plupart
des applications liées à l’analyse des structures multicouches. En fait, ces modèles sont plus
pratiques que les modèles 3D en termes de simplicité dans la modélisation et d'efforts de
calcul nécessaires. Parmi ces modèles, deux familles peuvent être distinguées (Fig.1.2):
Approche monocouche équivalente (ESL)1;
Approche par-couche (LW)2.
1
Equivalent Single Layer Approach
2
Layer-Wise Approach
8
Chapitre 1 Analyse bibliographique
Eeq
Real Laminate
E1
E2
E3
E4
E5
Layer-Wise (LW)
E1
E2
E3
E4
E5
Dans l’approche monocouche équivalente, qui est la plus adoptée par les chercheurs, la
plaque multicouche hétérogène est traitée comme une seule couche homogène équivalente via
la technique d'homogénéisation, voir [31-34]. Les variables de cette approche sont introduites
pour toute la plaque. Cela va conduire à un nombre de degrés de libertés indépendant par
rapport au nombre de couches (Fig.1.3). Selon Reddy [35], les théories du modèle ESL
fournissent une description suffisamment précise de la réponse globale de stratifié (i.e. la
flèche transversale, la fréquence fondamentale de vibration, la charge critique du flambement
et le moment résultant). Les principaux avantages des modèles ESL sont leur simplicité
inhérente et leur faible coût de calcul, dû au petit nombre des variables dépendantes.
9
Chapitre 1 Analyse bibliographique
Dans cette approche, on peut distinguer trois théories bien connues, à savoir:
La théorie classique de plaque est la plus simple et la plus ancienne parmi les autres
théories. Elle est initiée par Kirchhoff [37] en 1850, puis elle a été poursuivie par Love [38]
au début du XXe siècle, voir [7, 39-47]. Cette théorie est applicable uniquement pour les
plaques/coques, minces. Cependant, elle a donné des résultats erronés quand elle a été
utilisée dans l’analyse des plaques épaisses. L’hypothèse principale de la théorie CLPT est
que les lignes normales au plan médian avant déformation restent droites et perpendiculaires à
la surface moyenne après déformation (Fig.1.4). Cette hypothèse implique la négligence des
déformations dues au cisaillement transverse (CT)6 ( xz yz 0) et les
déformations/contraintes normales ( z z 0) .
Fig. 1.4 Géométries déformées et non-déformées d'un bord d'une plaque sous
l'hypothèse de Kirchhoff [2]
3
Classical Lamination Plate Theory
4
First order Shear Deformation Theory
5
Higher order Shear Deformation Theory
6
Cisaillement Transverse
10
Chapitre 1 Analyse bibliographique
La cinématique de cette théorie peut s'écrire d’une manière linéaire sous la forme suivante:
w0
u x , y , z u0 x , y z
x
w0
v x, y, z v0 x, y z (1.1)
y
w x, y, z w0 x, y
où
w,x et w, y sont les rotations dues à la flexion (sans cisaillement) dans les deux directions.
Cependant, cette théorie n’est pas adéquate pour l'analyse des composites multicouches
et notamment les structures sandwichs (inadaptée pour la modélisation du cœur), à cause de
l’effet du cisaillement transverse. Cet effet joue un rôle très important sur le comportement de
ces structures, puisque les rapports entre les modules de Young et les modules de cisaillement
sont compris entre 20 et 50 dans les matériaux composites modernes et entre 2,5 et 3,0 dans
les matériaux isotropes [48]. En plus, l’effet du cisaillement est la cause de différents modes
de rupture comme le délaminage au niveau des interfaces, fissure transverse, etc., alors que,
une prédiction précise des contraintes de cisaillement interlaminaires est nécessaire dans
l'analyse des défaillances des composites sandwichs.
Étant donné l'importance de la question, Libove et Batdorf [10], Yang et al. [49] et
Whitney [50], ont étendu la théorie de Reissner [51], Mindlin [52], dans le cas de plaques
isotropes épaisses, pour introduire l’effet de cisaillement dans les plaques sandwichs et les
plaques stratifiées anisotropes. Cela est devenu ce qu'on appelle ''théorie de déformation en
cisaillement du premier ordre'' (FSDT) ou la théorie de la plaque de déformation en
cisaillement uniforme (USDPT)7 [4, 53-58].
La théorie du premier ordre (FSDT) peut être considérée comme une amélioration par
rapport à la théorie classique (CLPT). Cette théorie est basée essentiellement sur l’hypothèse
suivante de Reissner-Mindlin: les lignes normales au plan médian avant déformation restent
7
Uniform Shear Deformation Plate Theory
11
Chapitre 1 Analyse bibliographique
droites mais pas forcément perpendiculaires à la surface moyenne après déformation (à cause
de l’effet du cisaillement transverse), (Fig.1.5). Ils ont également supposé que la contrainte
normale z est négligeable par rapport aux autres composantes du tenseur de contraintes
u x, y, z u0 x, y zx x, y
v x, y, z v0 x, y z y x, y (1.2)
w x, y, z w0 x, y
x et y sont les rotations de la normale autour des axes (x, y), respectivement
Fig. 1.5 Géométries déformées et non déformées d'un bord d'une plaque sous l'hypothèse de
Reissner-Mindlin [2]
12
Chapitre 1 Analyse bibliographique
Fig.1.6 Variation du déplacement et des contraintes transversales à travers l’épaisseur pour les
deux théories: théorie classique (à droite) et la théorie du premier ordre (à gauche) [66]
Les facteurs de correction du CT (FC), notés k, sont ajoutés pour corriger le décalage
entre la répartition des contraintes de cisaillement transversal réelles (3D) et celles supposées
dans la théorie du premier ordre. Ces facteurs sont introduits comme des paramètres dans les
relations constitutives entre les forces de cisaillement transversal et les déformations de
cisaillement transversal. Pour les matériaux isotropes homogène, le premier concept de FC a
été présenté par Reissner [51, 53, 67]. Il a proposé une valeur de 5/6 qui est largement
utilisée, en employant une méthode de calcul basée sur des considérations d’équilibre statique
et d’équivalence énergétique. Puis, Mindlin [52] a proposé une valeur de k 2 12 , en
utilisant la méthode d'analyse dynamique pour la première fois (équivalence de fréquence
propre associée au mode de vibration en CT ou de vitesse de propagation d’ondes de flexion
[68]). Une autre méthode a été présentée par Ayad [69] pour calculer les facteurs de
correction. Il les évalue par une comparaison entre l’énergie de cisaillement obtenue via les
équations d’équilibre et celle obtenue par la théorie du premier ordre.
Par ailleurs, les facteurs de correction pour les composites stratifiées et les sandwichs
dépendent essentiellement des propriétés matérielles du matériau (module de Young,
coefficient de poisson, etc.), du rapport d’aspect, de la forme des fibres ainsi que leurs
orientations, de la géométrie du stratifiée, du nombre de stratification et les conditions aux
limites, ce qui rend ces facteurs difficiles à déterminer. Dans ce contexte, plusieurs méthodes
ont été proposées pour estimer les facteurs de correction dans ces matériaux. La méthodologie
la plus courante est basée sur la comparaison de certaines réponses (l’énergie de déformation
de CT, la vitesse de propagation d'une onde de flexion, la fréquence naturelle du mode de
13
Chapitre 1 Analyse bibliographique
vibration, etc.) prédite par la théorie du premier ordre et celles obtenues par la théorie
d’élasticité tridimensionnelle [49, 56, 61, 70-72]. Yu [73] a utilisé une méthode basée sur la
comparaison des fréquences naturelles (les fréquences obtenues par la théorie d’'élasticité 3D
et celles obtenues par la théorie du premier ordre) pour les plaques sandwichs.
Whitney [58] a publié une analyse précise pour calculer les facteurs de correction des
plaques stratifiées symétriques/antisymétriques et les plaques sandwichs. Il les identifie par
une méthode basée sur des considérations des flexions cylindriques. L’auteur a montré que,
les facteurs de correction de CT ne s’approchent pas de la valeur classique 5/6, quand le
nombre de couches du stratifié augmente. La même remarque a été trouvée par Madabhusi-
Raman et Davalos [64]. Dans l’article de Vlachoutsis [60], les facteurs de correction de CT
dans le cas sandwich, sont estimés à partir d’une expression analytique, il montre que le
rapport d’anisotropie (Ef/Ec) influe sur les valeurs des FC. Certains auteurs, dont Pai [48], ont
présenté une nouvelle formulation de FC, en comparant les résultantes des contraintes de
cisaillement exactes et l’énergie de déformation de cisaillement avec ceux déterminées par la
théorie du premier ordre. L’article de Birman et Bert [62] résume les différentes méthodes
concernant la détermination de ces facteurs pour les poutres/plaques isotropes et sandwichs.
Pour pallier les limitations rencontrées dans les deux théories précédentes (i.e.,
détermination des contraintes de CT et leurs distributions au sens de l’épaisseur de la plaque),
des théories raffinées dites ''théories de déformation en cisaillement d’ordre supérieur''
(HSDT), ont été proposées par de nombreux chercheurs [74-83]. Ces théories sont basées sur
une variation non-linéaire du champ de déplacement (Fig.1.7), afin d’obtenir une meilleure
représentation des déformations et des contraintes de cisaillement transverse sans recours à
l'utilisation de facteurs de correction de CT. En addition, la théorie de déformation en
cisaillement d’ordre supérieur permet de prendre en compte un gauchissement éventuel de la
section droite de la plaque lors de la déformation [55].
Selon Touratier [84], un grand nombre de modèles d’ordre supérieur a été développé en
utilisant une méthode basée sur le développement en série de puissances (série de Taylor qui a
été initialement utilisée par Hildebrand et al. [85] en 1949) du champ de déplacement. Basset
[86] en 1890, semble avoir été le premier à suggérer cette idée. L’équation (1.3) présente
l'expression générale de ce champ qui s’écrit comme suit:
14
Chapitre 1 Analyse bibliographique
Fig. 1.7 Géométries non déformées et déformées d'une plaque selon les trois théories
CLPT, FSDT et HSDT [87]
où z est la coordonnée normale au plan médian x, y,0 , les coordonnées d'un point M étant
x, y, z , M 0 x, y,0 ; i x, y, z
Certains auteurs, dont Whitney et Sun [79], ont tenté d’améliorer la cinématique de la
théorie FSDT en proposant une théorie de second ordre (SSDT)8, déduite par l’équation (1.3)
3 4 2 3
avec i i = 0 et z z 0 , pour l'analyse statique et dynamique des plaques
stratifiées orthotropes. Le champ de déplacement considéré est:
z2
u x, y, z u0 x, y zx x, y x x, y
2
z2
v x, y, z v0 x, y z y x, y y x, y (1.4)
2
w x, y, z w0 x, y zz x, y
8
Second-Order Shear Deformation Theories
15
Chapitre 1 Analyse bibliographique
Une année plus tard, les mêmes auteurs [88] ont développé une autre théorie de second
ordre, dans laquelle les composantes de déplacement dans le plan (u, v) varient linéairement,
et la composante de déplacement transversal (w) varie d’une manière quadratique à travers
2 3 4 3
l’épaisseur i i i = 0 et z 0 pour étudier les coques stratifiées cylindriques.
La même cinématique a été utilisée par Kwon et Akin [89] et Fares [90]. Nelson et Lorch [77]
ont proposé également une théorie de déformation en cisaillement d'ordre deux en utilisant un
polynôme quadratique pour les trois composantes de déplacement i i = 0 . Tessler
3 4
[91] a développé une nouvelle théorie de second ordre pour l'analyse élastostatique linéaire de
plaques composites stratifiées épaisses, où les déplacements plans sont linéaires et le
déplacement transversal est quadratique à travers l'épaisseur. Sacco et Reddy [92] ont
appliqué la même cinématique de Whitney et Sun [79] pour analyser les plaques multicouches
en flexion. Ces derniers ont constaté que l'inclusion des termes de second ordre ne s'est pas
améliorée de manière significative par rapport à la solution linéaire obtenue par la théorie de
FSDT, à cause de l’absence de la variation quadratique des contraintes de CT.
Plusieurs auteurs ont alors développé d’autres théories d'ordre supérieur de plaque en se
basant sur l’équation (1.3). En 1975, Reissner [93] a remarqué que les polynômes d'ordre trois
pourraient être capables de rapprocher la distribution parabolique des contraintes de
cisaillement transversal à travers l'épaisseur de la plaque, en modifiant le champ de
déplacement choisi initialement en 1945. En 1977, Lo et al. [75], reprennent cette fois de
façon complète le champ de déplacement de Reissner [93].
Levinson [94] et Murthy [81] ont repris l’idée de Reissner [93] et Lo et al. [75], pour
présenter une théorie modifiée de déformation en cisaillement du troisième ordre (TSDT)9
similaire, dans le but d’analyser les plaques isotropes et les composites stratifiées,
respectivement. Mais dans ces théories, les auteurs utilisent les équations d'équilibre qui sont
variationnellement incompatibles avec les déplacements. En 1984, Reddy [78] fut le premier à
obtenir les équations d'équilibre d'une manière compatible en utilisant le principe des
déplacements virtuels. Il développe un modèle plus simple et efficace pour les plaques
9
Third-Order Shear Deformation Theories
16
Chapitre 1 Analyse bibliographique
composites multicouches, en se basant sur la théorie TSDT i 0 et z = z z = 0 ,
4 1 2 3
où les composants de champ de déplacement membranaires sont des fonctions cubiques et le
déplacement transversal est supposé constant selon l'épaisseur de la plaque (Eq.1.5).
u x, y, z u0 x, y zx x, y z 2 x x, y z 3 x x, y
v x, y, z v0 x, y z y x, y z 2 y x, y z 3 y x, y (1.5)
w x, y, z w0 x, y
mentionnés dans la théorie du premier ordre, ils dénotent les déplacements et les rotations de
Reddy a utilisé les conditions de nullité des contraintes de cisaillement xz , yz sur les
faces supérieures et inférieures de la plaque qui ont été proposées pour la première fois par
Levinson [94] dans l’analyse des plaques isotropes et par Murthy [81] et Bert [96] dans
l’analyse des plaques stratifiées. Ces conditions permettent de réduire le nombre de variables
additionnelles de déplacement (élimination de et ). Donc, le champ de déplacement de
l’équation (1.5) contiendra des paramètres ayant une signification physique comme suit:
4 w
u x, y, z u0 x, y zx x, y z 3 2 x 0
3h x
4 w
v x, y, z v0 x, y z y x, y z 3 2 y 0 (1.6)
3h y
w x, y, z w0 x, y
La théorie de Reddy est devenue très populaire parce qu’elle donne une meilleure
précision pour les contraintes de cisaillement transversal, et se traduisent par une distribution
parabolique de ces contraintes, sans recours à des facteurs de correction. En plus, le nombre
de variables est le même que dans la théorie FSDT (Cinq variables indépendantes). Cette
17
Chapitre 1 Analyse bibliographique
théorie est aussi connue sous le nom de ''la théorie de la plaque de déformation en
cisaillement parabolique (PSDPT)10''.
Certains auteurs, dont Senthilnathan et al. [97], ont utilisé cette théorie de Reddy pour
l'analyse du flambage des plaques isotropes et composites. Ils ont enrichi la composante de
déplacement normal en introduisant des variables supplémentaires comme suit:
w x, y, z w f x, y wc x, y (1.7)
Lee et al. [98] ont utilisé autrement la cinématique de Reddy pour l’analyse non-
linéaire des plaques isotropes épaisses.
En outre, Rohwer [99] a présenté une étude comparative de sept modèles basés sur des
théories différentes (CLPT, SSDT et HSDT) pour l’analyse des plaques composites
multicouches. Il a aussi montré les avantages et les limites de ces modèles et a conclu d'après
les résultats obtenus que la théorie de la plaque de Reddy est une bonne alternative parmi
toutes les théories en cisaillement de déformation d'ordre supérieur dans la littérature. Selon
Ferreira et al. [100], la théorie de Reddy est le meilleur choix pour les théories HSDT et reste
toujours la théorie la plus simple et la plus populaire dans l'analyse des plaques composites
stratifiées.
Kant et Swaminathan [30] ont proposé des solutions analytiques pour l’analyse statique
des plaques composites et des sandwichs simplement appuyées. Ils ont présenté deux modèles
d’ordre supérieur où le champ de déplacement est défini de telle façon que les trois
composantes (u,v,w) soient cubiques pour le premier modèle et pour le deuxième modèle les
deux composantes membranaires (u, v) soient cubiques et la composante transversale (w) soit
constante sur l’épaisseur de la plaque. La comparaison a été faite avec des solutions
d'élasticité (3D) données par Pagano [24] et d’autres modèles trouvés dans la littérature [54,
78, 97]. Cette comparaison a montré la précision et la fiabilité de ces deux modèles dans le
calcul de déplacement transversal et les contraintes. Un résumé de quelques modèles d'ordre
supérieur discuté précédemment est donné dans le tableau 1.
10
Parabolic Shear Deformation Plate Theory
18
Chapitre 1 Analyse bibliographique
Dans la littérature, il existe un certain nombre de modèles basés sur des théories non-
polynomiales [84, 102-110]. Ces modèles permettent de décrire le gauchissement à travers
l'épaisseur de la plaque, et donnent ainsi une distribution adéquate pour les déformations et les
contraintes de cisaillement transversal. En 1877, Levy [103] a développé une théorie raffinée
pour les plaques isotropes épaisses pour la première fois, en utilisant des fonctions
sinusoïdales (i.e. fonctions trigonométriques) dans le champ de déplacement. Mais, l'efficacité
de cette théorie de plaque particulière n'a pas été validée pour plus d'un siècle [111]. Stein
[108] a évalué cette théorie pour l’appliquer dans les plaques isotropes de forme modifiée.
Plus tard, Stein et Jegley [112] ont étudié l’effet de cisaillement transverse dans la flexion
cylindrique de plaques composites stratifiées, en utilisant la théorie de la plaque de
déformation en cisaillement trigonométrique (TSDPT)11. Ils ont prouvé que cette théorie
prédit les contraintes avec plus de précision que les autres théories.
Touratier [84] a proposé un nouveau modèle, qui est largement utilisé dans les plaques
stratifiées et sandwichs, en se basant sur la théorie TSDPT. L’auteur a utilisé une fonction
sinusoïdale pour décrire la cinématique de son modèle où le champ de déplacement est donné
par les expressions suivantes :
11
Trigonometric Shear Deformation Plate Theory
19
Chapitre 1 Analyse bibliographique
w0 z
x, y sin x
h
u x, y , z u0 x , y z
x h
w0 z
x, y sin y
h
v x, y, z v0 x, y z (1. 8)
y h
w x, y, z w0 x, y
En 1997, Idlbi et al. [114] ont présenté une étude comparative entre la théorie de Reddy
(PSDPT), celle de Touratier [84] (TSDPT) et les théories classiques (CLPT et FSDT), pour la
flexion et la vibration des plaques multicouches (Cross-Play et Sandwich). Ils ont montré que
le modèle raffiné de type sinus est préférable, en particulier lorsque les exigences de
continuité interlaminaire sont incluses.
Un nouveau modèle basé sur une autre théorie nommée la théorie de la plaque de
déformation en cisaillement hyperbolique (HSDPT)12 a été développé par Soldatos [115] pour
l’analyse des plaques monocliniques homogènes. L’auteur a utilisé la fonction sinus
hyperbolique pour les composantes membranaires tandis que le déplacement transversal est
considéré comme constant à travers l'épaisseur. La cinématique de ce modèle est de la forme
suivante:
w0
x, y h sinh z cosh x
z 1
u x , y , z u0 x , y z
x h 2
w0
x, y h sinh z cosh y
z 1
v x, y, z v0 x, y z (1.9)
y h 2
w x, y, z w0 x, y
Soldatos et Timarci [116] ont unifié la théorie HSDPT pour les plaques composites
multicouches. Le principal avantage de cette théorie unifiée est la capacité de modifier la
répartition des contraintes transversales [117]. Il existe d’autres travaux notamment ceux de
Karama et al. [110] qui ont proposé une nouvelle fonction dite "exponentielle" pour les
12
Hyperbolic Shear Deformation Plate Theory
20
Chapitre 1 Analyse bibliographique
poutres composites stratifiées. Elle est plus tard étendue aux plaques stratifiées [118]. Les
fonctions exponentielles s'avèrent beaucoup plus riches que les fonctions trigonométriques
(sinus et cosinus) dans leur série de développement de Fourier [110]. Elle est communément
connue par la théorie de plaque de la déformation en cisaillement exponentielle (ESDPT)13
dont le champ de déplacement s’écrit comme suit:
w0
u x , y , z u0 x , y z x, y ze2 z h x
2
x
w
v x, y, z v0 x, y z 0 x, y ze y
2
2 z h
(1.10)
y
w x, y, z w0 x, y
En 2006, Aydogdu [117] a présenté une étude dans laquelle il a comparé les différentes
théories d’ordre supérieur (PSDPT, TSDPT, HSDPT et ESDPT) avec la théorie d’élasticité
tridimensionnelle. L’auteur a montré que le déplacement transversal et les contraintes sont
mieux prédits par la théorie ESDPT, tandis que les deux théories PSDPT et HSDPT donnent
des prédictions plus précises des fréquences naturelles et des charges critiques de flambement
par rapport à d'autres théories.
Généralement, afin d'augmenter la précision des résultats, une croissance dans le degré
du polynôme de l’équation (1.3) peut être une alternative. Cette technique est peu adoptée par
les chercheurs à cause de son coût élevé en termes de calcul. Dans ce contexte, plusieurs
simplifications ont été proposées afin de réduire le nombre de paramètres de déplacement.
Une de ces simplifications consiste à raccourcir les derniers termes de la série de Taylor en
introduisant une ”fonction de cisaillement” [119]. La forme du déplacement proposée à
travers l’épaisseur est alors la suivante :
w0
u x, y , z u0 x , y z f z x x, y
x
w
v x, y, z v0 x, y z 0 f z y x, y (1.11)
y
w x, y, z w0 x, y
13
Exponential Shear Deformation Plate Theory
21
Chapitre 1 Analyse bibliographique
f(z) peut être considérée comme la fonction de cisaillement qui permet de déterminer la
distribution des déformations et des contraintes de cisaillement transversal à travers
l'épaisseur. Selon cette fonction f(z), nous pouvons distinguer quelques modèles d’ordre
supérieurs importants dans la littérature qui sont chronologiquement décrits comme suit:
z h2 z 2
Ambarsumyan [44] f ( z)
2 4 3
5z 4 z 2
Kaczkowski [120], Panc [121] et Reissner [122] f ( z ) 1 2
4 3h
4z2
Levinson [94], Murthy [81] et Reddy [78] f ( z ) z 1 2
3h
h z
Levy [103], Stein [108] et Touratier [84] f ( z) sin
h
z 1
Soldatos [115] f ( z ) h sinh z cosh
h 2
2 z h
2
Karama et al.[110, 118] f ( z ) ze
z
Ferreira et al. [123] f ( z ) sin
h
2 z h ln
2
Aydogdu [102] f ( z) z , 0
2 z h
2
Mantari et al. [124] f ( z ) z
z m cos z h
Mantari et al. [106] f ( z ) sin e
h
z m cosh z h
Mantari et Guedes Soares [126] f ( z ) sinh e
h
rz 2r
Grover et al. [127] f ( z ) sinh 1 -z , r 3
h h r2 4
rh 4r
f ( z ) cot 1 - , r 0.46
z h 4r 2 1
Sahoo et Singh [128]
Les théories issues d’une approche monocouche équivalente sont incapables de décrire
et prédire avec précision le comportement local des structures sandwichs y compris les
22
Chapitre 1 Analyse bibliographique
Dans cette approche, chaque couche est traitée individuellement (Fig.1.8), en imposant
les conditions de continuité en déplacements ou en contraintes aux interfaces. On applique à
chaque sous structure la théorie du premier ordre ou une théorie d’ordre supérieur, ce qui
conduit à un nombre de variables dépendant du nombre de couches. Les équations de
mouvement sont écrites pour chaque couche constitutive, alors que les conditions d'interface
sur les déplacements et les variables de contraintes transversales, sont présentées comme étant
supplémentaires [129].
23
Chapitre 1 Analyse bibliographique
Le développement des théories basées sur l’approche par-couche (layerwise) a été lancé
par Kao et Ross [131], Srinivas [132], Sun et Whitney [133], Swift et Heller [134], Seide
[135] et d’autres auteurs [136-142]. Kao et Ross [131] ont fait une première tentative pour les
poutres sandwichs multi-cœurs. Par la suite, Srinivas [132] a présenté une théorie de stratifiée
discrète pour l’analyse des plaques multicouches. Dans son modèle, les composantes de
déplacements plans sont considérées linéaires pour chaque couche et le déplacement
transversal est supposé constant sur l’épaisseur de la plaque. Sun et Whitney [133] ont utilisé
la même idée dans une autre procédure pour l’obtention de la réponse dynamique des plaques
stratifiées. Ils ont considéré la cinématique de la théorie FSDT pour chaque couche
séparément. Le champ de déplacement pour kth couche est exprimé comme suit:
u k x, y, z u0k x, y zxk x, y
v k x, y, z v0k x, y z yk x, y (1.12)
wk x, y, z w0k x, y
En 1987, Reddy [136] a proposé une approximation par-couche, plus fine, des champs
de déplacement suivant l’épaisseur d’une multicouche (Fig.1.9). Ce modèle est basé sur l'idée
que les déplacements peuvent être représentés comme une combinaison linéaire des produits
de fonctions des coordonnées dans le plan et les fonctions de l'épaisseur. Le champ de
déplacement correspondant dans la kème couche s'écrit comme suit:
24
Chapitre 1 Analyse bibliographique
n
u k
x, y, z U kj x, y jk z
j 1
n
v k x, y, z V jk x, y jk z (1.13)
j 1
n
wk x, y, z W jk x, y jk z
j 1
k k k
où u , v et w représentent les composantes des déplacements de la couche k dans les
Cho et al. [140, 141] ont présenté une théorie dite "Cho et al. Higher order, Individual-
Layer Theory" pour l’analyse des contraintes thermiques dans les plaques composites
stratifiées. Ils ont utilisé la cinématique suggérée par Lo et al. [75] pour l’appliquer à chaque
couche séparément. Cette théorie satisfait la compatibilité des déplacements et les conditions
de l’équation d’équilibre de contraintes de cisaillement au niveau de l’interface. Le champ de
déplacement s’écrit comme suit:
Plus tard, un nouveau modèle basé sur la théorie layerwise d’ordre supérieur a été
proposée par Lee et al. [148] pour la flexion des plaques composites stratifiées. Les auteurs
ont supposé que les composantes des déplacements dans le plan possèdent une variation
cubique pour chaque couche, ainsi, une variation quadratique pour les contraintes de
cisaillement avec des valeurs nulles aux surfaces libres (Eq.1.15). Dans ce modèle, la
continuité des déplacements plans et les contraintes de cisaillement transversal au niveau des
interfaces est assurée. On note également que le nombre de variables étant le même que celui
de la théorie de Mindlin (FSDT).
w0 4z3
u k x, y, z u0k z z 2 xk
x 3h
w0 4z3
v k x, y, z v0k z z 2 yk (1.15)
y 3h
wk x, y, z w0k x, y , k 1,...., n
Pour remédier au problème de temps de calcul élevé dans les théories de plaques à
couche discrète, Di Sciuva [149, 150], Murakami [151], Cho et Parmerter [152, 153], Lu et
Liu [154] et d’autres auteurs, ont présenté des modèles raffinés de type Zig-Zag (ZZT)14 dont
les déplacements plans ont une variation linéaire [149, 150, 155-157] (Figs.1.10, 1.11) ou
d’ordre supérieur [151-153, 158-166] (Fig.1.12), par morceaux, dans la direction de
l’épaisseur. L’avantage principal du champ de déplacements des modèles Zig-Zag réside dans
l’introduction de la perturbation de l’interface sans augmenter le degré de l’ordre des
équations fondamentales de la théorie de couche équivalente [119].
14
Zig-Zig Theories
26
Chapitre 1 Analyse bibliographique
Dans ces modèles, les conditions de continuité des déplacements et des contraintes de
cisaillement transversal aux niveaux d’interfaces (entre les peaux et l’âme) sont respectées
dans quelques travaux. Et le plus important, est que le nombre de variables cinématiques est
indépendant du nombre de couches (comme dans l’approche ESL). L’idée fondamentale dans
le développement de théories en Zig-Zag consiste à prendre un certain modèle de déplacement
(ou de contraintes) dans chaque couche, puis d'utiliser les conditions de compatibilité et de
l’équilibre d'interface pour réduire le nombre de variables inconnues [129].
Les théories LMT et AMT décrivent l’effet Zig-Zag par l’introduction de conditions de
continuité des contraintes de cisaillement transverses à travers les équations constitutives de
chaque couche et les relations entre la déformation et le déplacement. Par contre, la RMT
utilise les équations indépendantes dans la théorie [168].
15
Lekhnitskii Multilayered Theory
16
Ambartsumian Multilayered Theory
17
Reissner Multilayered Theory
27
Chapitre 1 Analyse bibliographique
En 1935, Lekhnitskii [169] a été le premier à proposer la théorie de Zig-Zag pour les
structures multicouches. Elle est obtenue par la résolution d’un problème d’élasticité des
poutres cantilevers. Dans cette théorie, l’auteur décrit la forme en Zig-Zag d'un champ de
déplacement dans le sens de l'épaisseur et l'équilibre des contraintes transversales aux
interfaces. En outre, des contributions exceptionnelles dans l’analyse des sous structures
(plaques/coques) composites ont été établies par Ambartsumian [170, 171]. En 1969, Whitney
[50] a par la suite étendu la théorie multicouche d’Ambartsumian aux problèmes de plaques
stratifiées anisotropes, symétriques et non-symétriques. Sur la base de la théorie AMT, Di
Sciuva [149] a développé une théorie de type Zig-Zag pour les plaques multicouches
anisotropes. Dans ce modèle, les composantes des déplacements plans possèdent une
distribution linéaire par morceaux dans le sens de l'épaisseur de la plaque, ainsi la continuité
28
Chapitre 1 Analyse bibliographique
des déplacements et des contraintes de cisaillement transversal aux interfaces est satisfaite (la
déformation transversale normale est négligée). Deux ans après, l’auteur a appliqué son
modèle (LZZ)18 aux problèmes de flexion, vibration libre et flambage, des plaques composites
multicouches [150].
u x, y, z u0 x, y zx x, y M z u zz x, y
v x, y, z v0 x, y z y x, y M z v zz x, y (1.16)
w x, y, z w0 x, y
2 zk
M z 1 k avec k
k
hk
k est la coordonnée adimensionnelle d’un point situé sur la couche k, u zz et v zz sont les
amplitudes effectives de l’effet du zig-zag, ayant le sens de déplacements.
Une année plus tard, Toledano et Murakami [173, 174] ont développé cette théorie en
ajoutant des termes du deuxième et troisième ordre au champ de déplacement de Murakami,
afin d’enrichir leur cinématique (Eq.1.17).
Dans le même contexte, Cho et Parmerter [152, 153] ont fourni une théorie efficace de
type Zig-Zag en combinant la théorie d'ordre supérieur de Reddy avec celle de Di Sciuva
18
Linear Zig-Zag model
19
Murakami’s Zig-Zag Function
29
Chapitre 1 Analyse bibliographique
[150]. Par ailleurs, Di Sciuva et Gherlone [176] ont proposé une théorie de zig-zag d’ordre
supérieur (HZZT)20, capable de calculer les contraintes de cisaillement transversal à partir des
équations constitutives. Cette théorie conduit à une variation parabolique des contraintes de
cisaillement par morceaux qui s'annulent sur les surfaces supérieure et inférieure.
Pour avoir plus d’informations et de détails sur les différentes approches, théories et
modèles d’éléments finis, des revues bibliographiques exhaustives, servant à analyser le
comportement des structures composites multicouches et notamment les structures sandwichs,
sont accessibles et sont données par [9, 35, 36, 59, 95, 114, 117, 129, 177-186].
Étant donné que les théories mentionnées ci-dessus présentent des limitations dans la
prédiction des contraintes de cisaillement transversal directement à partir des équations
constitutives, des théories nommées ‘’théorie Global-Local d'ordre supérieur’’ (GLHSDT),
basées sur des hypothèses de superposition double (global-local), ont été proposées pour
surmonter ces limitations. Dans cette théorie, le nombre de variables est indépendant du
nombre de couches (pour plus de détails, on peut consulter les deux références [185, 186]). Li
et Liu [187] furent les premiers à initier des travaux concernant cette approche. Les auteurs
ont proposé un modèle, où les conditions de continuités de déplacements et ceux des
contraintes de cisaillement interlaminaires sont satisfaites (Eq.1.19). Ainsi, les conditions de
nullité des contraintes de cisaillement sur les deux faces libres de la plaque sont également
imposées (Eq.1.20). Ce modèle est capable de prédire avec précision les contraintes de
cisaillement transversal directement à partir des équations constitutives. Dans les travaux
[185, 188-194], les auteurs ont utilisé la technique de superposition double (global-local),
pour prédire les contraintes interlaminaires dans les structures multicouches.
xz z zn 0 et xz z z0 0 (1.20)
20
Higher order Zig-Zag Theory
21
Global–Local Higher order Shear Deformation Theory
30
Chapitre 1 Analyse bibliographique
u k x, y, z uG x, y, z u Lk x, y, z
v k x, y, z vG x, y, z vLk x, y, z (1.21)
wk x, y, z wG x, y, z
avec
uG x, y, z u0 x, y zu1 x, y z 2 u2 x, y z 3u3 x, y
vG x, y, z v0 x, y zv1 x, y z 2 v2 x, y z 3 v3 x, y (1.22)
wG x, y, z w0 x, y zw1 x, y z 2 w2 x, y
et
k ak z bk (1.24)
où
Fig.1.13 Système des coordonnées, ordre de couche d’une plaque multicouche [195]
31
Chapitre 1 Analyse bibliographique
Deux ans après, les auteurs Li et Liu [196] ont procédé à une amélioration sur le modèle
précédent, en proposant un autre modèle raffiné (GLTSDT)22, basé sur la technique de
superposition globale-locale (The 1,2-3 Double-Superposition Theory), dont les termes par-
couche "local" du troisième ordre sont superposés avec des termes "global" du troisième ordre
des déplacements plans.
u k x, y, z uG x, y, z uLk x, y, z uˆLk x, y, z
v k x, y, z vG x, y, z vLk x, y, z vˆLk x, y , z (1.25)
wk x, y, z wG x, y, z
D’après l’équation (1.25), le terme local se compose de deux parties, la première partie
u k
L , vLk contient les termes de premier et de deuxième ordre, tandis que la seconde partie
uˆ k
L , vˆLk contient le terme de troisième ordre, comme suit:
uLk x, y, z u1k x, y k u2k x, y k2
vLk x, y, z v0k v1k k v2k x, y k2
(1.26)
uˆLk x, y, z u3k x, y k3
vˆLk x, y, z v3k x, y k3
uG x, y, z u0 x, y zu1 x, y z 2 u2 x, y z 3u3 x, y
vG x, y, z v0 x, y zv1 x, y z 2 v2 x, y z 3 v3 x, y (1.27)
wG x, y, z w0 x, y
Dans ce modèle, les auteurs ont imposé l’hypothèse de la continuité des déplacements et
des contraintes de cisaillement aux interfaces, qui permet de réduire considérablement le
nombre de variables inconnus. Toutefois, quand la théorie de double-superposition 1,2-3 est
comparée avec les théories par-couche, les coefficients d’ordre supérieur (en termes de z) sont
liés aux propriétés propres des couches et à la coordonnée z suivant l’épaisseur au lieu de
variables inconnues [195].
Certains auteurs, dont Zhen et Wanji [193] ont montré que la théorie global-local du
troisième ordre ou la théorie 1, 2-3 en double superposition de Li et Liu [196], donne des
22
Global–Local Third order Shear Deformation Theory
32
Chapitre 1 Analyse bibliographique
solutions très satisfaisantes aux problèmes de flexion, quand le nombre de couches est
inférieur à six. Cependant, la théorie GLTSDT présente des difficultés dans la prédiction des
contraintes de cisaillement lorsque le nombre de couches est supérieur à six. En addition,
l'exactitude de la solution de cette théorie diminue lorsqu'elle est appliquée à l'analyse des
plaques multicouches non-symétriques [195].
Dans le cadre de développement des éléments finis, certains chercheurs ont adopté la
théorie global-local d’ordre supérieur (GLHSDT), pour obtenir une bonne description des
structures sandwichs [189-192, 197, 198]. Citons par exemple, Wu et al. [189] qui ont étendu
la théorie de 1,2-3 double superposition proposée par Li et Liu [196] aux problèmes de
plaques stratifiées alternées (angle-ply). Ces derniers ont ajouté les termes Q45k de l’angle plis,
puisqu’ils ne sont pas nuls (Eq.1.32). En partant de cette théorie, global-local, les auteurs ont
formulé deux éléments finis de plaque, un élément quadrilatère à quatre nœuds (PQ) et un
autre élément triangulaire à trois nœuds (PT). Les éléments ont treize degrés de liberté par
nœud et le nombre d'inconnues est indépendant du nombre de couches. Les résultats
numériques trouvés par ces éléments montrent que les contraintes planes et les contraintes de
cisaillement transversal peuvent être calculées avec précision par les équations constitutives
directes, tandis que la contrainte normale zz est calculée à partir des équations d'équilibre.
La théorie a été appliquée dans l'analyse des comportements dynamiques et flambage des
plaques composites stratifiées et sandwichs [199, 200].
D’autres chercheurs ont tenté d’enrichir et d’améliorer les cinématiques des modèles
précités (CLPT, FSDT, etc.), en ajoutant des termes capables de donner une bonne description
des structures multicouches. Tessler et al. [201] ont présenté un nouveau modèle basé sur la
théorie du premier ordre enrichie pour les plaques sandwichs et les composites stratifiées. Ils
ont ajouté des fonctions zig-zag linéaires par morceau (Piecewise-Linear Zig-Zag functions),
qui permettent de fournir une représentation plus réaliste des états des déformations et des
contraintes de cisaillement. Le nombre des inconnus est relativement faible par rapport aux
autres théories. Cependant, la formulation de ce modèle n’impose pas les conditions de
33
Chapitre 1 Analyse bibliographique
u x, y, z u0 x, y z x x, y x k z x x, y
k
v x, y, z v0 x, y z y x, y y z y x, y
k k
(1.29)
w x, y, z w0 x, y
k k
où les fonctions zig-zag ( x , y ) sont définies par:
x k x, y, z 1 / 2 1 k u0 k 1 1 / 2 1 k u0 k
(1.30)
x, y, z 1 / 2 1 v
y
k k
1 / 2 1 v
0 k 1
k
0 k
avec
k z z k 1
h 1 1,1
k
k 1,...,n
Oh et al. [202] ont élaboré une nouvelle formulation d’éléments finis basée sur la
théorie du premier ordre améliorée (EFSDT)23 pour étudier avec précision le comportement
des composites stratifiées et des structures sandwichs. Pour cela, ils ont choisi le champ de
déplacement du modèle de Cho [153] (HZZT), pour enrichir leurs cinématiques. Par ailleurs,
Xiaohui et al. [203] ont développé un modèle efficace basé sur la théorie d’ordre supérieur
global-local, pour analyser les structures stratifiées et les sandwichs. La partie locale de ce
modèle est représentée par une cinématique de type zig-zag d’ordre supérieur. Les conditions
de continuité des contraintes interlaminaires et la nullité des contraintes de cisaillement sur les
surfaces libres des peaux supérieure et inférieure, sont satisfaites. On note également que ces
contraintes sont évaluées à partir de l'intégration des équations d'équilibre tridimensionnel. En
partant de ce modèle, les auteurs ont formulé un élément fini triangulaire à continuité (C0)
défini par six nœuds et 7ddl/nœud.
d’interface des couches. Dans la même période, Sedira et al. [205] ont proposé un élément
fini de forme quadrilatère possédant 4 nœuds, formulé sur la base du modèle variationnel
discrète de Mindlin (DDM)24, pour les plaques multicouches/sandwichs. Ils ont introduit la
fonction Zig-Zag de Murakami dans le modèle DDM afin d’améliorer la précision de
contraintes planes et de cisaillement transversal. Le modèle est caractérisé par une variation
linéaire par couche du déplacement, ce qui permet de satisfaire aux exigences de continuité
des contraintes.
u k x, y, z u0 x, y z x x, y z 2 x x, y z 3x x, y S k x x, y
v k x, y, z v0 x, y z y x, y z 2 y x, y z 3 y x, y S k y x, y (1.31)
wk x, y, z w0 x, y zw1 x, y z 2 x, y
Les solutions variationnelles ou analytiques exactes pour les problèmes qui ont des
géométries complexes, des conditions aux limites arbitraires ou bien dans le cas d’une non-
linéarité, sont terriblement difficiles à obtenir. Par conséquent, il faut recourir à des méthodes
numériques d'analyse capables de résoudre de tels problèmes. La méthode des éléments finis
s'est imposée dans les dernières décennies comme étant une méthode puissante et efficace, ce
qui lui permet d'être largement utilisée dans l’analyse du comportement complexe de
structures composites sandwichs.
24
Displacement Discrete Mindlin
25
Functionally Graded Material (FGM): est la version continue de la plaque stratifiée. Plus précisément, grâce à
de nouveaux procédés, il est possible de faire varier presque continument dans l’épaisseur les propriétés du
matériau constitutif. Ce nouveau concept est une façon de répondre aux difficultés soulevées par les
concentrations de contraintes présentées par les plaques stratifiées [123].
35
Chapitre 1 Analyse bibliographique
Depuis de nombreuses années, le développement des éléments finis basés sur la théorie
FSDT a attiré l'attention de nombreux chercheurs, comme Abel et Popov [207], Monforton et
Schmit [208], Ahmed [209], Khatua et Cheung [210, 211], Cook [212, 213] et Fazio et Ha
[214], sont considérés comme étant les premiers à utiliser la méthode des éléments finis dans
l’analyse du comportement des structures composites sandwichs. Monforton et Schmit [208]
ont formulé un élément fini rectangulaire à base du modèle déplacement, en adoptant la
théorie du premier ordre pour analyser les plaques et les coques sandwichs. Dans cet élément,
les variables de déplacement (membranaires et transversales) des deux peaux sont interpolées
par un polynôme d’Hermite.
En 1990, Lardeur [70] a développé deux éléments finis baptisés DST et DSQ (Discrete
Shear Triangular et Quadrilateral), pour l’analyse statique, dynamique et du flambage linéaire
des structures composites multicouches. Ces éléments sont construits selon la formulation
variationnelle mixte modifiée de Reissner, qui permet d’éviter le problème de verrouillage.
L’auteur a utilisé les équations d’équilibre pour calculer les déformations et les contraintes de
cisaillement transverse.
Dans le travail de Sheikh et al. [215], un élément fini triangulaire, possédant 16 nœuds
et 55 ddl, a été développé pour l’analyse des plaques composites multicouches de différentes
géométries. La distribution des degrés de liberté est présentée sur la figure 1.14. Dans cet
complet d’ordre quatre, tandis que les déplacements membranaires u et v et les deux
36
Chapitre 1 Analyse bibliographique
Dans le cadre de sa thèse doctorale, Tafla [216] a élaboré, à l’aide d’une formulation
mixte-hybride, un modèle élément fini quadrilatérale (MiSP/ml) à 4 nœuds et 3 ddl/nœud,
basé sur la théorie FSDT pour l’analyse des structures stratifiées et sandwichs. Le principal
avantage dans ce travail est la non nécessité de recourir aux facteurs de correction pour avoir
une distribution quadratique des contraintes de CT (approximation quadratique des
contraintes de cisaillement à travers l’épaisseur, qui exige la continuité C1). L’élément est une
extension de l’élément MiSP isotrope, proposé en 1993 par Ayad [69] sur la base du modèle
mixte projeté en cisaillement (MiSP: Mixed with Shear Projection). Le verrouillage en CT est
évité via une méthode nommée, ‘Assumed Natural Strain’ (ANS).
D’une manière similaire, Sakami [217] a développé dans son travail de doctorat, un
modèle d’éléments finis pour les plaques et coques multicouches, basé sur une nouvelle
approche variationnelle appelée ‘’Discret Displacement Mindlin’’ (DDM). Le modèle
d’éléments finis élaboré (DMQP/ml, DMQS/iso, DMQS/ml) est l’extension du modèle de
Katili [218], proposé initialement pour les plaques/coques homogènes isotropes. L’auteur a
introduit d’une manière discrète deux hypothèses (cinématique et mécanique) modifiées de
Mindlin dans leur formulation, ainsi, elle a utilisé la méthode ANS, pour éliminer le problème
de verrouillage.
Au début des années quatre-vingt, un grand nombre de chercheurs ont fait d'importantes
contributions à l'élaboration de nouveaux modèles d’éléments finis efficaces et capables de
donner une bonne description des structures multicouches, en se basant sur la théorie HSDT.
On commence par Kant et al. [219] qui sont les premiers à avoir publié des travaux
37
Chapitre 1 Analyse bibliographique
concernant la formulation d’éléments finis à base de cette théorie. Par la suite, Engblom et
Ochoa [220] ont adopté la théorie du deuxième ordre pour développer un élément fini destiné
à l’analyse statique des plaques stratifiées. Il s’agit de l’élément QHD40 qui est un
quadrilatère à 8 nœuds et 7 ddl/nœud. Les équations d’équilibre sont utilisées pour calculer
les contraintes de cisaillement transversal, tandis que les contraintes planes sont évaluées par
les équations constitutives.
Vers la fin des années 80, Pandya et Kant [221-224] ont étendu le model de Kant et al.
[219], pour présenter une série de travaux, destinés à l’étude du comportement des plaques
symétriques de type sandwich et starifiées. Les auteurs ont formulé un élément fini C0
isoparamétrique de type Lagrange possédant neuf nœuds, confiné pour l'analyse statique, en
se basant sur divers champs de déplacements (Eq.1.32-33).
u x, y, z z x x, y z 3 x x, y
v x, y, z z y x, y z 3 y x, y (1.32)
w x, y, z w0 x, y
u x, y, z z x x, y z 3 x x, y
v x, y, z z y x, y z 3 y x, y (1.33)
w x, y, z w0 x, y z 2 y x, y
Kant et Kommineni [225] ont également élaboré un élément fini, basé sur la théorie du
troisième ordre (Eq.1.5), pour l’analyse linéaire et non-linéaire géométrique, des plaques
composites FRC (Fibre Reinforced Composite) et des sandwichs. L’élément est un
quadrilatéral de type Lagrange à continuité C0, défini par 9 nœuds et 9 ddl par nœud. Par la
suite, Manjunatha et Kant [226, 227] ont présenté deux éléments finis isoparamétriques (9 et
16 nœuds), basés sur les deux modèles d’ordre supérieur précédents (HOST5 et HOST6),
pour une évaluation précise des contraintes interlaminaires dans les plaques composites
stratifiées et sandwichs. Les auteurs ont utilisé les équations constitutives pour évaluer les
contraintes planes, tandis que les contraintes de cisaillement transversal sont estimées en
utilisant les équations d’équilibre. Suite à cela, ils ont remarqué d’après les tests effectués,
38
Chapitre 1 Analyse bibliographique
Nayak et al. [231] ont décrit deux éléments finis quadrilatères isoparamétriques C0,
contenants 4 et 9 nœuds avec 7 ddl/nœud, afin d’étudier le comportement en flexion des
plaques: isotropes, composites stratifiées et sandwichs. Ils ont utilisé la méthode ANS qui a
été proposée en 1984 par Dvorkin et Bathe [234] et Huang et Hinton [235], pour éviter le
problème de verrouillage en cisaillement ainsi que les modes parasites. Les deux éléments
développés donnent des résultats précis et fiables (pour les déplacements et les contraintes)
notamment l’élément à 9 nœuds. Au cours de la même période, Sheikh et Chakrabarti [233]
ont publié un travail concernant le développement d’un élément fini triangulaire non
conforme pour l'analyse statique de plaques composites stratifiées. Celui-ci est défini par 6
nœuds contenant u, v, w, x , y , x et y comme des degrés de liberté à chacun . Un concept
semblable a été appliqué par Rezaiee-Pajand et al. [232] pour concevoir un nouvel élément
triangulaire à 13 nœuds. Les auteurs ont employé un polynôme d’interpolation approprié pour
le déplacement transversal, en plus, ils ont pris la dérivée normale du déplacement transversal
39
Chapitre 1 Analyse bibliographique
élément, où les nœuds (1, 2, 3, 13) possèdent 5 ddl u, v, w, x , y , les nœuds (4, 6, 7, 9, 10,
12) contiennent 6 ddl (u, v, w, x , y , w / n) et les nœuds (5, 8, 11) ont seulement un degré de
liberté w n .
Récemment, Tu et al. [236] ont formulé un élément fini rectangulaire ayant 9 ddl dans
chacun de ses neuf nœuds pour l’analyse du comportement flexionnel et vibratoire des
plaques composites starifiées et sandwichs, en utilisant la cinématique de Pandya and Kant
[237]. Dans cette étude aussi, les auteurs ont effectué une étude paramétrique, dans le but de
montrer l’effet de certains paramètres (rapport d’aspect, angle d’orientation des fibres et le
degré d’anisotropie) sur les contraintes planes, les contraintes de cisaillement transversal, les
déplacements et les fréquences fondamentales.
En 2013, Goswami et Becker [238] ont élaboré un élément fini de forme rectangulaire à
4 nœuds avec 11 ddl/nœud, pour analyser les plaques isotropes et les plaques
composites/sandwichs. La formulation de cet élément a été basée sur la théorie de troisième
ordre (TSDT). Les auteurs ont pris les deux déformations angulaires ( xz et yz ) comme des
degrés de liberté. Dans la même année, Lee and Kim [239] ont développé un élément fini de
forme quadrilatéral à 4 nœuds et 7 ddl/nœud (u, v, w, x , y , w x, w y) , pour l’analyse des
plaques composites multicouches. Les auteurs ont introduit la méthode ANS pour pallier le
problème de verrouillage en CT.
De nombreux travaux de recherche ont été effectués dans le but de développer des
modèles d’éléments finis avec des cinématiques plus adéquates que celles qui ont précédé.
40
Chapitre 1 Analyse bibliographique
Ces éléments sont basés sur l’approche layerwise, et sont capables de donner une bonne
description du comportement des structures sandwichs (dans le cas des réponses locales aux
niveaux des interfaces). En 1982, Spilker [240] a développé un élément fini isoparamétrique à
8 nœuds pour l’analyse des plaques multicouches minces et épaisses, en utilisant une
formulation hybride en contraintes. Les composantes des déplacements plans sont continus et
linéaires pour chaque couche, et toutes les composantes de contrainte sont incluses et
interpolées de façon indépendante dans chaque couche. La même formulation a été adoptée
par Fazio et al. [241]. Dans la référence [242], un élément fini isoparamétrique C0 à 9 nœuds,
formulé via une approche mixte, est élaboré pour prédire les contraintes interlaminaires et les
déplacements dans les plaques sandwichs épaisses. Les auteurs ont proposé, pour chaque
couche, un champ cubique et quadratique pour les déplacements dans le plan et hors-plan,
respectivement. Dans leur modèle, uniquement la continuité des déplacements est assurée.
En outre, Lee et Fan [243] ont décrit un nouveau modèle dans lequel la théorie du
premier ordre (FSDT) est utilisée pour modéliser les peaux, tandis que les déplacements au
niveau de l’âme sont exprimés en fonction des déplacements des deux peaux (Eq.1.34). Dans
ce modèle, les contraintes de cisaillement transversal varient linéairement, et la déformation
selon l’épaisseur du cœur est constante. À partir de ce modèle, les auteurs ont formulé un
élément fini isoparamétrique à 9 nœuds et 5 ddl/nœud, afin d’étudier la flexion et la vibration
des plaques sandwichs.
uc
u0b u0t
t bxb t txt z u b
0 u0t t b b
x
t txt
h 2h 2 4
vc
v0b v0t
t b yb t t yt z vb
0 v0t t
b b
y
t t yt
(1.34)
h 2h 2 4
wc
w b
0 w0t z w
b
x wxt
h 2
Par ailleurs, certains auteurs proposent des modèles éléments finis tridimensionnels,
citons par exemple Oskooei et Hansen [244] qui ont développé un nouveau modèle basé sur
l’approche layerwise, pour l’analyse statique des plaques sandwichs ayant des peaux
composites stratifiées. Dans ce modèle, les peaux sont modélisées comme des plaques de
Riessner-Mindlin, tandis que l’âme est modélisée en utilisant la théorie du troisième ordre, où
le champ de déplacement est cubique pour les composantes u et v et quadratique pour la
41
Chapitre 1 Analyse bibliographique
composante w, dans le sens de l’épaisseur. En partant de ce modèle, les auteurs ont construit
un élément fini tridimensionnel (3D) contenant 240 degrés de liberté (Fig.1.17).
Évoquons également l’élément de Ramtekkar et al. [245, 246] qui est proposé
notamment pour l’évaluation précise des contraintes interlaminaires dans les plaques
stratifiées et sandwichs. Il s’agit d’un élément brique à 18 nœuds avec 6 ddl/nœud (trois
déplacements ( u, v, w ) et trois contraintes ( xz , yz , z )), comme le montre la figure 1.18.
Les auteurs utilisent une cinématique basée sur l’approche layerwise, où les composantes du
champ de déplacement u, v et w sont des fonctions cubiques pour chaque couche. Les
exigences de continuité des déplacements ainsi que les contraintes de cisaillement transversal
dans ce modèle sont satisfaites.
Par la suite, Demiray et al. [247, 248] ont proposé un nouveau modèle pour analyser la
non-linéarité statique des plaques sandwichs. Dans ce modèle, les peaux sont représentées par
la théorie classique de Kirchhoff, et l’âme est représentée par la théorie du deuxième ordre
(les composantes u et v sont quadratiques et la composante w est linéaire au sens de
l’épaisseur de la plaque). Sur la base de ce modèle, les auteurs ont développé un élément fini
prismatique (3D) défini par 6 nœuds: les nœuds (1, 2, 3) sont situés dans la face supérieure et
les nœuds (4, 5, 6) dans la face inférieure, où chaque nœud contient 6 ddl (Fig. 1.19). La non-
linéarité géométrique est prise en compte dans ce modèle. Pour éviter le problème de
verrouillage en cisaillement, les auteurs ont utilisé la technique d'intégration réduite.
tc t f t
v1t u1t x3 1
2
tc t f t
v2t u2t x3 2 (1.35)
2
v3t u3t
43
Chapitre 1 Analyse bibliographique
tc t f b
v u x3 1
b b
1 1
2
tc t f b
v2b u2b x3 2 (1.36)
2
v3 u3
b b
4 x 2
tc t tf 3 c
v u1 u1 c u1 u1 x3 1 1 c 1 1 x3
1 t 1 t
c b b b t b
1 1
tc 2
1
2 t 4 2t
4 x 2
tc t tf 3 c
v u2 u2 c u2 u2 x3 2 2 c 2 2 x3
1 t 1 t
c b b b t b
1 2 (1.37)
2
2 t 4 2t tc 2
v3c u3t u3b c u3t u3b x3
1 1
2 t
En outre, Linke et al. [249] ont développé un élément fini (3D) contenant 44 degrés de
liberté, comme le montre la figure 1.20. Chaque peau contient 5 ddl par nœud, tandis que
l’âme contient seulement un degré de liberté. Cet élément est proposé pour l’analyse statique
et l’analyse de la stabilité des plaques sandwichs. La formulation de cet élément est basée sur
le modèle de trois couches, où les peaux sont décrites par des éléments de plaques classiques
basés sur l’hypothèse de Love-Kirchhoff et l’âme est décrite par la théorie du troisième ordre
(polynôme cubique pour les composantes u et v, et quadratique pour la composante w).
Citons aussi, Ramesh et ses collègues [250] qui ont proposé un élément fini triangulaire
(HE-45-LT) de 45 nœuds avec 7 ddl/nœud (Fig.1.21), pour prévoir des contraintes
interlaminaires dans les plaques composites stratifiées. La construction de cet élément est
basée sur deux théories, la théorie du troisième ordre (TSDT) et la théorie layerwise de Reddy
[136]. L'élément est capable de prédire avec précision les forces de cisaillement, les moments
de torsion et spécialement les contraintes de cisaillement transversal à travers l'épaisseur.
Cependant, l'analyse des plaque de plusieurs couches nécessite énormément de temps de
calcul, lorsque le nombre de couches est important car la variation de ce dernier est
proportionnelle au nombre d’inconnues.
Durant la même période, Ćetković et Vuksanović [251] ont utilisé le modèle layerwise
de Reddy [136] pour étudier le comportement statique, dynamique et de flambage des plaques
composites stratifiées et sandwichs. Dans ce modèle, les composantes des déplacements plans
ont une variation linéaire par couche (piecewise), et la composante transversale est constante,
à travers l’épaisseur de la plaque. Les auteurs ont achevé leur travail par une étude
paramétrique, afin d’étudier l’effet de l’anisotropie, l’orientation des fibres, le nombre de
couches et le rapport d’aspect sur les déplacements, les contraintes planes, les contraintes de
cisaillement transversal, les fréquences fondamentales et la charge critique de flambement. Ils
ont trouvé que ces paramètres jouent un rôle important dans le comportement de plaques
multicouches.
45
Chapitre 1 Analyse bibliographique
Darılmaz [252] a publié un travail dans lequel un élément fini solide hybride
tridimensionnel (Fig. 1. 22) à 8 nœuds et 48 ddl a été employée dans l'analyse de plaques
sandwichs. La formulation de cet élément est basée sur le principe variationnel de Hellinger-
Reissner (développée déjà en 2005 par le même auteur [253]). L’objectif de l’investigation a
été de prédire avec précision les contraintes hors-plan ( xz , yz , z ).
Dans un autre ordre d’idées, Mantari et al. [104] ont développé un nouveau modèle
layerwise, dans lequel les composantes du déplacement plan (u et v) sont exprimées par des
fonctions tangentielles trigonométriques, tandis que la composante normale (w) est constante
selon l’épaisseur de la plaque. En partant de ce modèle, un élément fini quadrilatère
isoparamétrique de continuité C 0, défini par 4 nœuds et 7 ddl/nœud
u , v , w , w x, w y,
0 0 0
c
x , yc , est utilisé pour étudier la flexion des panneaux multicouches
épais. Les résultats obtenus par ce modèle sont satisfaisants pour la flèche et les contraintes
qu’elles soient planes ou cisaillement.
Pour améliorer les modèles précédents et corriger certains problèmes rencontrés dans
les travaux cités ci-dessus, plusieurs auteurs ont axé leurs travaux vers la recherche des
éléments finis de plus en plus élaborés en utilisant la théorie dite zig-zag. Leurs objectifs est
d’obtenir un meilleur compromis entre la précision et le faible coût de calcul. Di Sciuva [254,
255] a élaboré deux éléments finis conformes de continuité C1 pour les plaques multicouches.
Ces derniers possèdent dix degrés de liberté par nœud. Il s'agit de l'élément quadrilatère à
quatre nœuds (RHQ40) et de l'élément triangulaire à trois nœuds (RHT30). L’auteur a utilisé
46
Chapitre 1 Analyse bibliographique
En 2003, Topdar et al. [163] ont développé un élément fini pour l’analyse statique des
plaques sandwichs, en utilisant le modèle zig-zag d’ordre supérieur (HZZT)26 de Cho et
Parmerter [153]. L’élément est rectangulaire à quatre nœuds avec huit degré de liberté
26
Higher order Zig-Zag Theory
47
Chapitre 1 Analyse bibliographique
nu 1 nl 1
u u0 S k
x z zk H z zk Txk z k H z k x z 2 x z 3
k 0 k 0
nu 1 nl 1
v v0 S yk z zk H z zk Tyk z k H z k y z 2 y z 3 (1.38)
k 0 k 0
w w0 x, y
D’une manière similaire au modèle de Cho, Chakrabarti et Sheikh [158, 260] ont
développé un nouvel élément fini triangulaire à 6 nœuds et 7 ddl/nœud pour l’analyse statique
des plaques sandwichs. Les déplacements dans le plan varient selon une loi cubique alors que
le déplacement hors plan est supposé constant à travers l’épaisseur de la plaque. La variation
des contraintes de cisaillement transversal à travers l’épaisseur est représentée par des
fonctions paraboliques pour chaque morceau, alors qu’elles sont nulles sur les deux faces de
la plaque (supérieure et inférieure). L’élément formulé est non conforme dû à la condition de
continuité de pente normale qui n’est pas satisfaite. Le seul problème signalé dans cet élément
est que ce dernier nécessite une continuité C1 de déplacement transversal aux interfaces.
48
Chapitre 1 Analyse bibliographique
En 2008, Pandit et al. [162] ont amélioré le modèle zig-zag de Chakrabarti et Sheikh
[260] afin d’étudier les plaques sandwichs avec une âme molle (soft core). Le déplacement
transversal est susceptible de varier comme une fonction quadratique dans la couche centrale,
tandis qu’il reste constant dans les deux faces. En partant de ce modèle, les auteurs ont
développé un élément fini quadrilatère à continuité C0 défini par 9 nœuds et 11 ddl/nœud. En
plus, les auteurs ont pallié le problème rencontré ci-dessus pour la continuité C1 en utilisant
des fonctions de formes séparées pour définir les dérivés de déplacement transversal.
En 2013, Khandelwal et al. [160] ont développé un modèle d’élément fini efficace à
continuité (C0) destiné à l'analyse statique des plaques sandwichs à âme molle. Ce modèle est
élaboré en utilisant deux théories combinées: la théorie de zig-zag d’ordre supérieur (HZZT)
et la méthode de moindres carrés (LES)27. L’élément est un quadrilatère isoparamétrique à 9
nœuds et 11ddl/nœud, formulé sur la base de la théorie HZZT. Cette théorie a fourni des
résultats très satisfaisants pour les contraintes planes et le déplacement transversal par rapport
aux autres théories 2D existantes dans la littérature, alors que les contraintes de cisaillement
transversal sont évaluées d’une manière précise par la méthode LES. D’autres scientifiques,
dont Singh et al. [261] et Chalak et al. [262] ont pris l’idée de Cho et Parmerter [153] et
Pandit et al. [162] pour formuler leurs modèles d’éléments finis des plaques multicouches.
1.4 Conclusion
La revue bibliographique a montré que les modèles basés sur l’approche monocouche
équivalente deviennent moins efficaces lorsqu’il s’agit d’étudier les problèmes liés à
l’endommagement (i.e. délaminage) et aux effets locaux (i.e. estimation des contraintes
interlaminaires). Par conséquent, l'utilisation d'une description par-couche devient
indispensable.
27
The Least Square Error method
49
Chapitre 1 Analyse bibliographique
Actuellement, les éléments finis créés sur la base de l’approche par-couche (LW, ZZT
and GLHSDT) gagnent du terrain, parce qu’ils fournissent une représentation
cinématiquement correcte des contraintes interlaminaires ainsi que les déplacements au
niveau d’interface.
La partie suivante portera sur la présentation d’un nouveau modèle de type layerwise
pour la modélisation des plaques sandwichs et stratifiées. Ce modèle, qui a pour objectif
l’enrichissement de la bibliothèque des éléments finis existants, est destiné non seulement
pour l’analyse des plaques sandwichs, mais également pour l’étude des plaques isotropes et
des composites multicouches.
50
PARTIE II
51
Chapitre 2 Proposition d’un nouveau modèle layerwise pour l’analyse des structures multicouches
Chapitre 2
Proposition d’un nouveau modèle layerwise
pour l’analyse des structures sandwichs et des
multicouches
2.1 Introduction
Dans la partie bibliographie précédente, nous avons vu que les modèles issus d’une
approche monocouche équivalente, présentent des difficultés lorsqu'il s'agit de la prédiction et
de l’estimation des déplacements et des contraintes interlaminaires aux interfaces. Donc,
l'utilisation d'une description par-couche devient indispensable pour avoir une bonne
prévision de l’état des contraintes. Cette description est plus répandue car chaque couche est
traitée individuellement, en imposant les conditions de continuité en déplacements et en
contraintes au niveau des interfaces.
Par conséquent, ce chapitre a pour objectif de présenter un nouveau modèle, basé sur
l’approche layerwise, pour analyser le comportement des composites multicouches et
notamment les plaques sandwichs qu’elles soient symétriques ou non-symétriques.
Contrairement à tous les autres modèles layerwise présentés dans la littérature, le nombre de
degrés de liberté est indépendant du nombre de couches. Le modèle proposé est issu d'une
combinaison de deux théories; la théorie de déformation en cisaillement du premier ordre
(FSDT) est adoptée pour modéliser les deux peaux, tandis que l’âme est modélisée en utilisant
la théorie de déformation en cisaillement du troisième ordre (TSDT). Ces deux théories
combinées, satisfont les exigences de continuités des déplacements interlaminaires ainsi que
les conditions de compatibilité. Le nouveau modèle permet de calculer les contraintes de
cisaillement transversal d’une manière précise en utilisant directement les équations
constitutives plutôt que les équations d'équilibre et sans recours à des facteurs de correction.
52
Chapitre 2 Proposition d’un nouveau modèle layerwise pour l’analyse des structures multicouches
Les plaques sandwichs sont des structures qui se composent de trois couches
principales, (figure 2.1) dont deux peaux minces (supérieure-inférieure) d’épaisseurs ht ,
L'épaisseur totale de la plaque h est la somme des épaisseurs de toutes les couches. Le plan
x, y du système de coordonnées coïncide avec le plan moyen de la plaque.
On suppose dans tous les cas qu'il n'y a pas de glissement aux interfaces, i.e. les couches sont
parfaitement collées.
53
Chapitre 2 Proposition d’un nouveau modèle layerwise pour l’analyse des structures multicouches
uc u0 z xc z 2 xc z 3 xc
vc v0 z yc z 2 yc z 3 yc (2.1)
wc w0
où u0 , v0 et w0 sont respectivement, les composantes de déplacement dans le plan et hors
plan, du plan moyen de la plaque sandwich. xc et yc représentent les rotations du premier
ordre de la sections autour des axes y et x, respectivement. Les rotations xc , yc , xc et yc
sont des termes d’ordre supérieur, où x et y symbolisent les termes du second ordre
c c
xc yc xc yc
1 x1 , y1 , xy1 x
,
y
,
y
x
xc yc xc yc
2
2 2
x , y , xy 2
,
x y y
,
x
c c c c
3 x3 , y3 , xy3 x , y , x y (2.3)
x y y x
w
s 0 yz 0 , xz 0 yc 0 , xc 0
y x
w
, 2 , 2
s
1 1
yz
1
xz
c
y
c
x
, 3 ,3
s
2 2
yz
2
xz
c
y
c
x
54
Chapitre 2 Proposition d’un nouveau modèle layerwise pour l’analyse des structures multicouches
xx x 0 x1 x 2 x 3
0 1 2 2 3 3
yy y z y z y z y
0 1 2 3
xy xy xy xy xy
(2.4)
0 1 2
yz yz yz 2 yz
0 z 1 z 2
xz xz xz xz
La théorie utilisée pour modéliser les deux peaux est la théorie du premier ordre (définie
dans la section 1.1.2.1.2). Les composantes de déplacement (u, v et w) sont exprimées par les
expressions suivantes:
u x, y, z u0 x, y z xf x, y
f
v x, y, z v0 x, y z yf x, y , f t, b
f
(2.5)
w x, y, z w0 x, y
f
h h
ut uc c z c xt
2 2
h h
vt vc c z c yt (2.6)
2 2
wt w0
respectivement.
h hc c hc c hc c
2 3
uc c u0 2 x 4 x 8 x
2
avec, (2.7)
h hc c hc c hc c
2 3
vc c v0 2 y 4 y 8 y
2
55
Chapitre 2 Proposition d’un nouveau modèle layerwise pour l’analyse des structures multicouches
hc c hc c hc c h
2 3
ut u0 x x x z c xt
2 4 8 2
h h2 h3 h
vt v0 c yc c yc c yc z c yt (2.8)
2 4 8 2
wt w0
h h
ub uc c z c xb
2 2
h h
vb vc c z c yb (2.9)
2 2
wb w0
respectivement.
h hc c hc c hc c
2 3
uc c u0 2 x 4 x 8 x
2
avec, (2.10)
h hc c hc c hc c
2 3
vc c v0 2 y 4 y 8 y
2
h h2 h3 h
ub u0 c xc c xc c xc z c xb
2 4 8 2
hc c hc c hc c h
2 3
vb v0 y y y z c yb (2.11)
2 4 8 2
wb w0
ut hc x hc x hc x
u0 hc xt
c 2 c 3 c
t
(z )
x x 2 x 4 x 8 x 2 x
xx
v v h y hc2 y hc3 y h y
c c c t
tyy t 0 c (z c )
y y 2 y 4 y 8 y 2 y
ut vt u0 v0 hc xc y hc2 xc y
c c
xy
t
y x y x 2 y x 4 y x
(2.12)
h hc
3 c c t t
(z )
c x y x y
8 y x 2 y x
w0
tyz yt
y
w0
xzt xt
x
D’une autre manière, l’équation (2.12) peut être écrite sous forme matricielle comme suit:
yy (2.13)
t m y 2 y 4 y 8 y 2 y
xy
u0 v0 hc xc y hc2 xc y hc3 xc y hc xt y
c c c t
y x 2 y
x 4 y x 8 y x 2 y x
xt
t f x w0
xx xt
t
t
x
tf yyt f z y
, c t
t xz
(2.14)
y
yz w
t f y 0 t
xy t t
y
x y
y x
57
Chapitre 2 Proposition d’un nouveau modèle layerwise pour l’analyse des structures multicouches
ub hc x hc x hc x
u0 hc xb
c 2 c 3 c
b
(z )
x 2 x 4 x 8 x
xx
x 2 x
v v h y hc2 y hc3 y h y
c c c b
byy b 0 c (z c )
y y 2 y 4 y 8 y 2 y
xyb
y x y x 2 y x 4 y x
hc3 xc yc
hc xb y
b
(z )
8 y x 2 y x (2.15)
w0
byz yb
y
w0
xzb xb
x
De même, l’équation (2.15) peut être sous-divisée en trois termes; membranaire, flexionnelle
et cisaillement transversal. Soit:
yy
y
(2.16)
x 2 y 4 y 8 y 2 y
xy
b m
u0 v0 hc x y hc x y hc x y hc x y
c c 2 c c 3 c c b b
y x 2 y x 4 y x 8 y x 2 y x
xb
x w0
xxb f xb
b
b x
bf yyb f z y
, cb b
xz
(2.17)
b f y yz w0 yb
xy b b y
x y
y x
58
Chapitre 2 Proposition d’un nouveau modèle layerwise pour l’analyse des structures multicouches
La théorie élémentaire des plaques fait l'hypothèse que les contraintes normales zz sont
hypothèse est étendue à la théorie des stratifiés. Les relations contraintes-déformations dans le
système de coordonnées locales des fibres sont données par :
E1 12 E2 E2
Q11 , Q12 , Q22 , Q66 G12 , Q44 G23 , Q55 G13 (2.19)
1 12 21 1 12 21 1 12 21
avec E1 , E2 sont des modules d'Young du matériau dans les directions 1 et 2, G12 , G13 et G23
Fig. 2.3 Les axes de référence X , Y , Z ; et les axes de fibre du couche (1,2,3)
La loi de Hooke généralisée dans le système de coordonnées globales (Fig.2.3) peut être
59
Chapitre 2 Proposition d’un nouveau modèle layerwise pour l’analyse des structures multicouches
c 2 s 2 2sc 0 0
2 2
s c 2sc 0 0
T sc sc c 2 s 2 0 0 (2.20)
0 0 0 c s
0 s c
0 0
D’où,
c cos et s sin , dans laquelle est l'angle de la fibre comme représenté sur la figure
2.3.
Les relations contraintes-déformations, pour une couche k, exprimées dans le repère global,
s’écrivent donc:
où Q ij sont les coefficients de la matrice de rigidité d'une couche k donnée. Ils sont exprimés
Les deux peaux (supérieure et inférieure) sont considérées comme des plaques en
matériaux composites stratifiées. Par conséquent, les relations contrainte-déformation d’une
couche kth dans le système de coordonnées global peuvent être exprimées par:
60
Chapitre 2 Proposition d’un nouveau modèle layerwise pour l’analyse des structures multicouches
k k k
xxf Q11 Q12 0 0 Q16 xxf
f f
yy Q21 Q22 yy
0 0 Q26
i.e. Q
yz 0 yz
f f k k k
0 Q44 Q45 0
(2.23)
f 0 f
xzf 0 Q54 Q55 0 xzf
xy Q61 Q62 0 0 Q66 xy
f top, bottom
L’âme est considérée comme un matériau composite orthotrope et les relations ''contrainte-
déformation'' sont données sous la forme suivante:
Les efforts résultants dans l’âme sont obtenus par l’intégration des contraintes à travers son
épaisseur.
Le vecteur des résultantes en membrane noté N x, y , d’une couche k, est défini par:
N x hc x
2
N ( x, y ) N y y dz (2.25)
N xy h2c xy
M x hc x
2
M ( x, y ) M y y z dz (2.26)
M xy h2c xy
N x hc
2 x 2
N ( x, y ) N y y z dz (2.27)
- hc
N xy 2 xy
M x hc
2 x 3
M ( x, y ) M y y z dz (2.28)
- hc
M xy 2 xy
hc
Vx 2
xz
V x, y dz (2.29)
Vy hc yz
2
Sx Rx
yz z, z dz
2
xz
S 2
(2.30)
y Ry hc
2
En substituant les équations (2.4) et (2.23) dans les expressions des efforts résultants
(Eq.2.25-2.30), on obtient:
62
Chapitre 2 Proposition d’un nouveau modèle layerwise pour l’analyse des structures multicouches
hc 0 1 2 3
M xy
2 xy xy xy xy
hc
hc
S y 2 yz0 1yz 2 yz
2
z Q 0 Q z 1 Q z 2 dz (2.36)
Sx
hc xz xz xz
2
hc
Ry 2
yz0 1yz 2 yz
2
z Q 0 Q z 1 Q z 2 dz
2
(2.37)
Rx hc xz xz xz
2
hc
Qij 1, z, z , z , z , z , z dz
2
Q 1, z, z , z , z dz
2
N x A11 A12 A16 xx0 B11 B12 B16 1xx D11 D12 D16 xx2 E11 E12 E16 xx3
0 3
N y A21 A22 A26 yy B21 B22 B26 1yy D21 D22 D26 yy2 E21 E22 E26 yy
N A 0
A66 xy B61 B66 1
D66 2 3
xy 61 A62 B62 xy D61 D62 xy E61 E62 E66 xy
(2.39)
M x B11 B12 B16 xx0 D11 D12 D16 1xx E11 E12 E16 xx2 F11 F12 F16 xx3
3
M y B21 B22 B26 yy0 D21 D22 D26 1yy E21 E22 E26 yy2 F21 F22 F26 yy
M B 0 D
B66 1 E
D66 2 F
E66 F66 3
xy 61 B62 xy 61 D62 xy 61 E62 xy 61 F62 xy
(2.40)
63
Chapitre 2 Proposition d’un nouveau modèle layerwise pour l’analyse des structures multicouches
N x D D12 D16 xx0 E11 E12 E16 1xx F11 F12 F16 xx2 G11 G12 G16 xx3
11
3
N
y 21 D D22 D26 yy0 E21 E22 E26 1yy F21 F22 F26 yy2 G21 G22 G26 yy
D 0 E
D66 1 F
E66 2 G
F66 G66 3
61 61 61 61 G62
D62 xy E62 xy F62 xy xy
N xy
(2.41)
M x E E12 E16 xx0 F11 F12 F16 1xx G11 G12 G16 xx2 H11 H12 H16 xx3
11 3
M y E21 E22 E26 yy0 F21 F22 F26 1yy G21 G22 G26 yy2 H 21 H 22 H 26 yy
E E66 xy0 F61 F66 1xy G61 G62 G66 xy2 H 61 3
H 66 xy
M xy 61
E62 F62 H 62
(2.42)
Il est intéressant de noter que ( B, D, E, F , G, H ) sont des matrices de couplage. Les termes
de couplage de même ordre s’annulent pour une âme symétrique.
Vy A44 A45 yz0 B44 B45 1yz D44 D45 yz2
(2.43)
Vx A54 A55 xz0 B54 B55 1xz D54 D55 xz2
Les équations constitutives de l’âme peuvent être écrites sous la forme contractée suivante:
N A B D E 0
M
B D E F 1
(2.46)
N D E F G 2
M E F G
H 3
V A B D s 0
s s s
s 1
S B D E s
s s
(2.47)
R s 2
D E s F s s
Selon la théorie de déformation en cisaillement du premier ordre (FSDT), les expressions des
résultantes de la peau supérieure s'écrivent:
64
Chapitre 2 Proposition d’un nouveau modèle layerwise pour l’analyse des structures multicouches
N t At Bt 0 mt
t t
M B Dt 0 tf (2.48)
Tt 0 t
Act
0 c
par:
hc
ht k 1
n layer h
1, z, z dz Q 1, z, z dz
2
(k ) (k )
Aijt , Bijt , Dijt Q ij 2
ij
2
(i, j 1, 2, 6)
hc k 1 h
k
2
(2.49)
hc
ht
k 1
2 n layer h
Q
(k ) (k )
t
Aij Q ij dz ij dz (i, j 4, 5)
hc k 1 h
k
N b Ab Bb 0 mb
b b b
M B Db 0 f (2.50)
Tb 0 Acb b
0 c
hc
k 1
n layer h
1, z, z dz Q 1, z, z dz
2
(k ) (k )
Aijb , Bijb , Dijb Q ij 2
ij
2
(i, j 1, 2, 6)
h k 1 h
k
c hb
2
avec,
hc
(2.51)
k 1
2 n layer h
Q
(k ) (k )
b
Aij Q ij dz ij dz (i, j 4, 5)
h k 1 h
k
c hb
2
2.3 Conclusion
Dans ce chapitre, nous avons présenté la formulation théorique d’un nouveau modèle de
plaque sandwich multicouches, basée sur l’approche layerwise, i.e. combinaison de deux
théories: (i) la théorie de déformation en cisaillement du premier ordre (FSDT) pour la
modélisation des deux peaux; et (ii) la théorie du troisième ordre (TSDT) pour modéliser
l’âme. Contrairement à tous les autres modèles layerwise présentés dans la littérature, le
nombre de degrés de liberté est indépendant du nombre de couches. Ceci a pour avantage la
réduction du temps de calcul lorsque le nombre de couches augmente. Sur la base du modèle
65
Chapitre 2 Proposition d’un nouveau modèle layerwise pour l’analyse des structures multicouches
proposé, nous allons aborder, dans le chapitre III, la formulation et la validation d’un nouvel
élément fini de forme rectangulaire pour l’analyse statique des plaques sandwichs, qu’elles
soient symétriques ou non-symétriques, ainsi que les plaques isotropes et les composites
multicouches.
66
Chapitre 3 Développement d’un nouvel élément fini rectangulaire (RSFT52)
Chapitre 3
3.1 Introduction
Dans ce chapitre, nous allons alors procéder au développement d’un nouvel élément
fini, bi-dimensionnelle (2D), construit sur la base du modèle layerwise, décrit dans le chapitre
précédent [Réf]. Cet élément est destiné à l’analyse statique des plaques sandwichs, qu’elles
soient symétriques ou non-symétriques, ainsi que les composites multicouches. L’élément
développé est de forme rectangulaire, nommé RSFT52 (Rectangular Sandwich First Third
with 52-degrees of freedom), défini par quatre nœuds et treize degrés de liberté par nœud.
Chaque nœud contient: deux degrés de liberté de rotation pour chaque peau, six degrés de
liberté de rotation aussi pour l’âme, alors que les trois degrés de liberté de type translation
sont communs aux trois couches. Après la formulation, une série de tests ont été entrepris
dans l’objectif de valider l’élément développé (RSFT52).
67
Chapitre 3 Développement d’un nouvel élément fini rectangulaire (RSFT52)
n
x , y N i x, y i (3.1)
i 1
où i ui vi wi xic yic xic yic xic yic xit yit xib yib est le vecteur des déplacements
T
Les fonctions d’interpolation bi-linéaire classique N i de type Lagrange ont les expressions
suivantes:
68
Chapitre 3 Développement d’un nouvel élément fini rectangulaire (RSFT52)
1
N1 a x b y
4ab
1
N2 a x b y
4ab
(3.2)
1
N3 a x b y
4ab
1
N4 a x b y
4ab
L’approximation des variables cinématiques des trois couches peut être exprimée comme suit:
b. L’âme:
4 4
xc x, y N i x, y xic , yc x, y N i x, y yic
i 1 i 1
4 4
xc x, y N i x, y xic , yc x, y N i x, y yic
i 1 i 1
c. Peau supérieure:
4 4
d. Peau inférieure:
4 4
B
e e
B i e
1 e 1 2 e 2
; i e
B
3 e
i e
e e e
B
3 0 0
; s s i e
(3.7)
B
e
B i e
1 e 1 2 e 2 e
s s
; s s i e
69
Chapitre 3 Développement d’un nouvel élément fini rectangulaire (RSFT52)
Où les matrices B1 , B 2 , B3 , B 0 , B1 et B 2 peuvent être obtenues à
s s s
N i
0 0 0 0 0 0 0 0 0 0 0 0
x
N i
B1 0 0 0 0 0 0 0 0 0 0 0 0
y
N i N i
0 0 0 0 0 0 0 0 0 0 0
y x
N i
0 0 0 0 0 0 0 0 0 0 0 0
x
2 N i
B 0 0 0 0 0 0 0 0 0 0 0 0
y
N i N i
0 0 0 0 0 0 0 0 0 0 0
y x
N i
0 0 0 0 0 0 0 0 0 0 0 0
x
N i
B 3 0 0 0 0 0 0 0 0 0 0 0 0 (3.8)
y
N i N i
0 0 0 0 0 0 0 0 0 0 0
y x
N i
0 0 x Ni 0 0 0 0 0 0 0 0 0
B 0
0 0 N i 0 0 0 0 0 0 0 0
s
0 Ni
y
0 0 0 0 0 2 Ni 0 0 0 0 0 0 0
B1
s
0 0 0 0 0 0 2 Ni 0 0 0 0 0 0
0 0 0 0 0 0 0 3Ni 0 0 0 0 0
B 2
s
0 0 0 0 0 0 0 0 3Ni 0 0 0 0
70
Chapitre 3 Développement d’un nouvel élément fini rectangulaire (RSFT52)
où
N i N i N i N i
0 0 0 c 0 d 0 e 0 f 0 0 0
x x x x
N i N i N i N i
Bmt 0 0 0 0 c 0 d 0 e 0 f 0 0
y y y y
N i N i N i N i N i N i N i N i
0 0 0 c y c
x
d
y
d
x
e
y
e
x
f
y
f
x
0 0
avec : (3.10)
2 3
hc h h hc
c , d c
, e c
, f
2 4 8 2
N i
0 0 0 0 0 0 0 0 0 0 0 0
x
N i
B f 0 0 0 0 0 0 0 0 0
t
0 0 0
y
N i N i
0 0 0 0 0 0 0 0 0 y x
0 0
N i
0 0 x 0 0 0 0 0 0 Ni 0 0 0
Bct
0 0 N i 0 0 0 0 0 0 0 Ni 0 0
y
où,
N i N i N i N i
0 0 0 f 0 d 0 h 0 0 0 c 0
x x x x
N i N i N i N i
Bm 0 0 0
b
0 f 0 d 0 h 0 0 0 c
y y y y
N i N i N i N i N i N i N i N i
0 0 0 f
y
f
x
d
y
d
x
h
y
h
x
0 0 c
y
c
x
avec :
hc hc2 hc3 hc
f , d , h , c
2 4 8 2
71
Chapitre 3 Développement d’un nouvel élément fini rectangulaire (RSFT52)
N i
0 0 0 0 0 0 0 0 0 0 0 0
x
N i
B f 0 0 0 0 0 0 0 0 0 0 0
b
0 (3.12)
y
N i N i
0 0 0 0 0 0 0 0 0 0 0 y x
N i
0 0 x
0 0 0 0 0 0 0 0 Ni 0
Bc
b
0 0 N i
0 0 0 0 0 0 0 0 0 Ni
y
Pour établir la relation entre les forces et les déplacements, on a utilisé le principe du
travail virtuel, soit:
U W 0 (3.13)
U U c + U t U b (3.14)
où
vc vt vb
(3.15)
L'expression de travail virtuel des efforts interne (Eq. 3.15) peut s’écrire comme suit:
72
Chapitre 3 Développement d’un nouvel élément fini rectangulaire (RSFT52)
hc
xxc yyc yyc xyc xyc xzc xzc yzc yzc dz dAc
2
U c
xx
Ac
hc
2
hc
ht
xxt yyt yyt xyt xyt xzt xzt yzt yzt dz dAt
2
t
xx
(3.16)
At
hc
2
hc
xxb yyb yyb xyb xyb xzb xzb yzb yzb dz dAb
2
b
xx
Ab
hc
hb
2
Le travail virtuel effectué par les forces externes est donné par:
b a
W f x, y w dxdy
b a
(3.17)
Dans laquelle f x, y est la charge statique transversale répartie sur la surface de l’élément.
La substitution des expressions des efforts résultants (Eqs 2.48–2.49) dans l'expression
de travail virtuel dans l’âme (Eq.3.15), nous permet d’écrire :
b a
b a
T
2
s F s s 2 dxdy f w dxdy 0
b a
b a
K e c 0 T
A B0 B0 B B1 B0 D B 2
T T
B
b a
s s s s s s
1 T
2 2 T
0 2 T
1
B E B B D B B E B
s
s
s s
s
s s
s
s
Les mêmes étapes ont été suivies pour élaborer la matrice de rigidité des deux peaux. On aura
donc:
a. Peau supérieure
b a
K e Bmt A t Bmt Bmt B t B tf B tf B t Bmt
t T T T
b a
membranaire couplage membrane-flexion couplage flexion-membrane
(3.22)
+ B tf D B tf Bct Ac Bct dxdy
T t T t
flexionnelle cisaillement
b. Peau inférieure
b a
Ke b B b T A b B b B b T B b B b B b T B b B b
m m m f f m
b a
membranaire couplage membrane-flexion couplage flexion-membrane
(3.23)
+ B bf D B bf Bcb Ac Bcb dxdy
T b T b
flexionnelle cisaillement
74
Chapitre 3 Développement d’un nouvel élément fini rectangulaire (RSFT52)
La matrice de rigidité globale s’obtient par l’assemblage des matrices de rigidité élémentaires.
Dans cette section, nous allons évaluer les performances de l'élément développé
(RSFT52), en termes de précision, de convergence et de stabilité, à travers une série de tests.
Des plaques isotropes, composites stratifiées et sandwichs symétriques et non symétriques
ainsi que différents cas de chargement, géométrie et conditions aux limites ont été considérés.
Les résultats obtenues seront comparés avec les solutions obtenus analytiquement et ceux
obtenus par d’autres modèles d’éléments finis disponibles dans la littérature.
Les conditions aux limites considérées dans cette étude sont présentées dans le tableau
3.1. Les notations de ces conditions sont:
CSCS: Deux bords opposés encastrés et deux autres bords simplement appuyés
L’arrangement du maillage utilisé dans tous les exemples est illustré dans la figure 3.2.
Le premier exemple concerne une plaque isotrope carrée simplement appuyée avec
différents rapports d’épaisseur a h . Quant au deuxième exemple, il concerne une plaque
75
Chapitre 3 Développement d’un nouvel élément fini rectangulaire (RSFT52)
Tableau 3.1 Détails de conditions aux limites utilisées pour les exemples numériques.
b
w0 xc xc xc xt xb 0 at y
2
Simplement appuyé SSSS
b
w0 yc yc yc yt yb 0 at x
2
w0 xc yc xc yc 0
Encastré CCCC
xc yc xt yt xb yb 0
a
Encastré à x
Encastré - Simplement 2
CSCS
appuyé b
Simplement appuyé à y
2
a
Encastré à x
2
Encastré-Libre CFCF
b
Libre à y
2
76
Chapitre 3 Développement d’un nouvel élément fini rectangulaire (RSFT52)
Les résultats obtenus par l’élément développé (RSFT52) sont en bon accord avec la
solution analytique de Kant [82]; et ceux trouvés par des modèles d’éléments finis, basés sur
des théories du premier ordre [263, 264] et des théories d’ordre supérieur [219].
Tableau 3.2 Déflexion normalisé ( w) d’une plaque carrée isotrope simplement appuyée.
wG
w ; xz xz (3.26)
hp0 p0
transversal ( xz ) sont représentés dans le tableau 3.3. La comparaison a été faite avec la
solution analytique de Kant et al. [265]. Les résultats obtenus par l’élément RSFT52 montrent
les performances de la formulation proposée.
77
Chapitre 3 Développement d’un nouvel élément fini rectangulaire (RSFT52)
a b b
a a w , ,0 xz 0, , 0
2 2 2
b h
Présent élément Élasticité [265] Présent élément Élasticité [265]
Dans cette section, trois exemples ont été considérés pour mettre en évidence la
performance du modèle proposé dans l’analyse du comportement des plaques stratifiées.
Dans le premier exemple, une plaque starifiée carrée simplement appuyée (SSSS),
considérée (figure 3.3). La plaque est constituée de trois couches identiques (0 / 90 / 0) dont
les caractéristiques mécaniques sont présentées dans le tableau 3.4. Différents types de
maillages (4×4, 6×6, 8×8, 10×10, 12×12, 14×14 et 16×16), ainsi que différents rapports
w 100 E2 h3
w (3.27)
a 4 q0
78
Chapitre 3 Développement d’un nouvel élément fini rectangulaire (RSFT52)
Fig.3.3 Plaque stratifiée carrée à trois couches, simplement appuyée sous chargement
doublement sinusoïdal
Tableau 3.4 Propriétés des matériaux (normalisés) utilisés pour les plaques stratifiées.
Propriétés
Location E1 E2 G12 G13 G23 12 21
mécaniques
Plaques stratifiées All layer 25E E 0.500E 0.50E 0.20E 0.25
FEM-T6-HSDT: Élément fini triangulaire à six nœuds, basé sur la théorie d’ordre
supérieur, présenté par Sheikh et Chakrabarti [233].
FEM-T45-LW: Élément fini triangulaire à quarante-cinq nœuds, basé sur la théorie
layerwise, présenté par Ramesh et al. [250].
FEM-T45-TSDT: Élément fini triangulaire à quarante-cinq nœuds, basé sur la théorie
du troisième ordre, présenté par Ramesh et al. [250].
FEM-T6-HZZT: Élément fini triangulaire à six nœuds, basé sur la théorie de zig-zag
d’ordre supérieur, présenté par Chakrabarti et Sheikh [260].
FEM-Q4-RTOST: Élément fini quadrilatère à quatre nœuds, basé sur la théorie du
troisième ordre raffinée , présenté par Kulkarni et Kapuria [266].
FEM-Q9-HZZT: Élément fini quadrilatère à neuf nœuds, basé sur la théorie de zig-zag
d’ordre supérieur , présenté par Khandelwal et al.[160].
79
Chapitre 3 Développement d’un nouvel élément fini rectangulaire (RSFT52)
FEM-SFR8O: Éléments fini hexaèdre à huit nœuds, basé sur le concept " Space Fiber
Rotation", présenté par Meftah [267]
D’après les résultats qui figurent sur le tableau 3.5, on remarque que le modèle
d’élément fini développé RSFT52 donne des résultats très satisfaisants, en terme de stabilité
et rapidité de convergence, notamment dans le cas des plaques épaisses (h a 0.25) , où les
résultats sont en excellent accord avec les solutions de références, à titre d’exemple pour un
maillage de 4×4, on trouve une erreur de 0.0273 %.
On considère une plaque starifiée carrée simplement appuyée et soumise à une charge
uniformément répartie. Les propriétés mécaniques sont celles de l’exemple 3.3.3.1. Deux
types de stratifications croisées, dont lesquelles, trois couches (0 / 90 / 0) et cinq couches
(0 / 90 / 90 / 90 / 0) , ont été étudiées. Dans les deux cas, les couches de même orientation ont
la même épaisseur. Différents maillage (4×4, 6×6, 8×8, 10×10, 12×12, 14×14 et 16×16) ont
été considérés. Dans les tableaux 3.6 et 3.7, on présente les résultats du déplacement
transversal au centre de la plaque pour les deux stratifications. Les résultats obtenus sont
normalisés par la formule précédente (3.27).
Les résultats obtenus sont comparés avec ceux issus de la solution analytique donnée
par la théorie d’ordre supérieur de Reddy [78], la méthode EFG28, basée sur la théorie FSDT,
présentée par Belinha et al. [268], et à ceux obtenus par d’autres modèles numériques de la
littérature [163, 233, 267, 269]. Par exemple: Xiao et al. [270] ont utilisé le logiciel
commercial ABAQUS pour modéliser ce type de plaque, en choissisant un élément fini brique
(3D) à 20 nœuds, avec un maillage de 4800 éléments, pour le quart de la plaque. Les auteurs
ont également utilisé deux éléments MQ-MLPG et TPS-MLPG, qu’ils ont développé sur la
base de '' la théorie d’ordre supérieur et la déformation normale de la plaque '' (HOSNDPT)29;
et la formulation de point fini (MLPG)30 basée sur la méthode de la fonction de base radiale
multiquadrique (RBFs)31. La comparaison montre la précision et l’efficacité de l’élément
RSFT52, dans l’analyse des plaques stratifiées multicouches.
28
Element-Free Galerkin
29
Higher Order Shear and Normal Deformable Plate Theory
30
Meshless Local Petrov–Galerkin
31
Radial Basis Functions
80
Chapitre 3 Développement d’un nouvel élément fini rectangulaire (RSFT52)
Tableau 3.5 Déplacement transversal non-dimensionnel ( w) d’une plaque stratifiée carrée à trois
couches (0 / 90 / 0) simplement appuyée sous chargement doublement sinusoïdal.
h
Rapport d’épaisseur
Références Modèles EF a
0.25 0.1 0.05
Fig.3.4 Convergence du déplacement transversal d'une plaque carrée starifiée à trois couches,
Tableau 3.6 Déplacement transversal non-dimensionnel ( w) d’une plaque stratifiée carrée à trois
couches (0 / 90 / 0) , simplement appuyée, sous chargement uniforme.
a
Rapport d’épaisseur
Références Modèles EF h
4 10 20
Tableau 3.7 Déplacement transversal ( w) d’une plaque stratifiée carrée à cinq couches
(0 / 90 / 90 / 90 / 0) , simplement appuyée, sous chargement uniforme.
a
Rapport d’épaisseur
Références Modèles EF h
5 10 20
82
Chapitre 3 Développement d’un nouvel élément fini rectangulaire (RSFT52)
4Q E1 E2 1 2 23 a
ww avec Q 4G12 et S (3.28)
4
12 S h q0 1 12 21 h
Les résultats obtenus sont en accord avec la solution analytique, fournie par Pagano
[271] et ceux obtenus par des modèles d’éléments finis de Lardeur [70], Engblom et Ochoa
[220] Sakami [217] et Sedira [195]. On peut remarquer aussi que le présent modèle est
applicable pour les deux types de plaques, minces ou épaisses.
Tableau 3.8 Déplacement transversal non-dimensionnel ( w) d’une plaque stratifiée carrée à neuf
couches (0 / 90 / 0 / 90 / 0 / 90 / 0 / 90 / 0) simplement appuyée, sous chargement doublement
sinusoïdal.
83
Chapitre 3 Développement d’un nouvel élément fini rectangulaire (RSFT52)
appuyée, soumise à une charge doublement sinusoïdale (Fig. 3. 6) a été considérée. Les
propriétés matérielles de l’âme et celles des peaux sont représentées dans le tableau 3.9.
L’épaisseur de chaque peau est égale à 0.1 h et celle de l’âme est 0.8 h , où h étant
l’épaisseur totale de la plaque. La plaque est analysée pour différents rapports d’épaisseurs
16×16).
Tableau 3.9 Propriétés des matériaux (normalisés) utilisés pour les plaques sandwichs.
Propriétés
E1 E2 G12 G13 G23 12 21
mécaniques
Âme 0.04E 0.04E 0.016E 0.06E 0.06E 0.25
Plaques sandwich
Peaux 25E E 0.500E 0.5E 0.20E 0.25
100 E2 h3 a b
w w , (3.29)
a 4 q0 2 2
xz
, yz
h
q0 a
xz , yz (3.30)
85
Chapitre 3 Développement d’un nouvel élément fini rectangulaire (RSFT52)
yz a b
xz 0,
a b a
Référence Modèles EF ,0 , 0, 0 w , , 0
h 2 2 2 2
86
Chapitre 3 Développement d’un nouvel élément fini rectangulaire (RSFT52)
yz a b
xz 0,
a b a
Référence Modèles EF ,0 , 0, 0 w , , 0
h 2 2 2 2
87
Chapitre 3 Développement d’un nouvel élément fini rectangulaire (RSFT52)
88
Chapitre 3 Développement d’un nouvel élément fini rectangulaire (RSFT52)
Cet exemple a été proposé et étudié pour la première fois par Azar [22], puis repris par
d’autres auteurs, dont Khatua et Cheung [210], Topdar et al. [163] et Chakrabarti et Sheikh
de chaque peau est de 0.7112 mm. L'étude a été faite en considérant trois types de matériaux:
Matériau 1:
Peaux: E11 E22 68.8 GPa , G13 G23 27.52 GPa et 13 0.3
L’âme: E11 E22 6.88 1012 GPa , G13 G23 0.2064 GPa et 13 0.3
Matériau 2:
Peaux: E11 68.8 GPa, E22 27.52 GPa, G12 12.9 GPa, G13 G23 27.52 GPa et 13 0.3
89
Chapitre 3 Développement d’un nouvel élément fini rectangulaire (RSFT52)
Matériau 3:
Peaux: E11 E22 68.8 GPa , G13 G23 27.52 GPa et 13 0.3
L’âme: E11 E22 6.88 1012 GPa , G13 G23 0.3131 GPa et 13 0.3
Dans le dernier cas (matériau 3), la plaque sandwich est considérée encastrée sur les
quatres côtés (CCCC). Les résultats du déplacement transversal au centre de la plaque, avec
différents maillages (4×4, 8×8, 12×12 et 16×16), sont rapportés dans le tableau 3.11. Les
résultats obtenus par l’élément RSFT52 présentent une très bonne concordance avec la
solution analytique donnée par Azar [22] et ceux obtenus par les modèles d’éléments finis de
Khatua et Cheung [210], Topdar et al. [163] et Chakrabarti et Sheikh [260].
Cet exemple a été choisi dans le but d’évaluer la polyvalence du modèle d’élément fini
développé dans l’analyse des plaques sandwichs ayant des peaux composites stratifiées. Il
s’agit d’une plaque sandwich stratifiée non-symétrique 0 / 90 / C / 0 / 90 soumise à une
charge sinusoïdale (figure 3.10). L’épaisseur de l’âme est de 0.8 h , tandis que celle de chaque
couche de stratifiée elle est de 0.05 h , h étant l'épaisseur totale de la plaque. Les propriétés
mécaniques des matériaux utilisés pour l’âme et les deux peaux sont résumées dans le tableau
3.9. Dans ce problème, trois types de conditions aux limites (CCCC, SCSC et SSSS) avec six
rapports d'épaisseur h a 0.01, 0.05, 0.1, 0.2, 0.25 et 0.5 , ont été considérés.
a b
Références Modèles EF w 100 mm , ,0
2 2
91
Chapitre 3 Développement d’un nouvel élément fini rectangulaire (RSFT52)
élevé que ceux de plus faibles rapports d'épaisseurs, i.e. 0.01 à 0.05, respectivement, pour les
deux conditions aux limites (CCCC et SCSC). La cause principale derrière tout cela est due à
l’effet de la flexibilité transversale de l’âme et aux effets de la déformation en cisaillement,
qui sont plus prononcés dans le cas des rapports d'épaisseurs élevés. Le même genre
d'observations peut également être noté pour la contrainte de cisaillement transversal non-
3.3.4.4 Plaque sandwich carrée avec des peaux stratifiées alternés (angle-ply)
Ce problème a été traité par Chakrabarti et Sheikh [158] et Khandelwal et al. [160]. Il
s’agit d’une plaque sandwich carrée ayant des peaux stratifiées alternées (angle-ply)
f ( 90) C f ( 90) , simplement appuyée sur les quatres cotés et soumise à une
charge uniformément répartie. Les mêmes propriétés du problème précédent (3.3.4.3) qu’elles
soient géométriques ou mécaniques, ont été adoptées.
L'analyse est effectuée pour trois types de maillage (8×8, 12×12 et 16×16), trois
rapports d'épaisseur h a 0.05, 0.1 et 0.2 ainsi que pour trois angles d’orientation
92
Chapitre 3 Développement d’un nouvel élément fini rectangulaire (RSFT52)
yz a b
xz 0,
h b a
Références Modèles EF ,0 , 0, 0 w , , 0
a 2 2 2 2
93
Chapitre 3 Développement d’un nouvel élément fini rectangulaire (RSFT52)
yz a b
xz 0,
h b a
Références Modèles EF ,0 , 0, 0 w , , 0
a 2 2 2 2
94
Chapitre 3 Développement d’un nouvel élément fini rectangulaire (RSFT52)
yz a b
xz 0,
h b a
Références Modèles EF ,0 , 0, 0 w , , 0
a 2 2 2 2
0.25 Présent élément (8×8) RSFT52 0.1515 0.1436 5.6655
Présent élément (12×12) RSFT52 0.1563 0.1384 5.6329
Présent élément (16×16) RSFT52 0.1579 0.1326 5.6225
Khandelwal et al. [160] FEM-Q9-HZZT 0.1747 0.1327 5.6638
Chalak et al. [262] FEM-Q9-HZZT 0.1667 0.1907 5.6453
95
Chapitre 3 Développement d’un nouvel élément fini rectangulaire (RSFT52)
h
Références Modèles EF xz yz w
a
0
0.2 Présent élément (8×8) RSFT52 0.3033 0.3033 6.3228
Présent élément (12×12) RSFT52 0.3242 0.3242 6.2940
Présent élément (16×16) RSFT52 0.3343 0.3343 6.2909
Pagano [24] Elasticity solution 0.3436 0.3413 6.2981
Khandelwal et al. [160] FEM-Q9-HZZT (Const)* 0.3982 0.3993 6.3001
Khandelwal et al. [160] FEM-Q9-HZZT (Equil)** 0.3455 0.3332 6.3001
Chakrabarti et Sheikh [158] FEM-T6-HZZT 0.3482 - 6.3016
30
0.2 Présent élément (8×8) RSFT52 0.3185 0.3185 5.9480
Présent élément (12×12) RSFT52 0.3382 0.3382 5.9326
Présent élément (16×16) RSFT52 0.3472 0.3472 5.9308
Khandelwal et al. [160] FEM-Q9-HZZT (Const) 0.4376 0.4321 5.9579
Khandelwal et al. [160] FEM-Q9-HZZT (Equil) 0.3064 0.2951 5.9579
Chakrabarti et Sheikh [158] FEM-T6-HZZT 0.3762 - 5.9463
96
Chapitre 3 Développement d’un nouvel élément fini rectangulaire (RSFT52)
h
Références Modèles EF xz yz w
a
45
0.2 Présent élément (8×8) RSFT52 0.2955 0.2955 5.6998
Présent élément (12×12) RSFT52 0.3144 0.3144 5.6615
Présent élément (16×16) RSFT52 0.3231 0.3231 5.6516
Khandelwal et al. [160] FEM-Q9-HZZT (Const) 0.4612 0.4484 5.6329
Khandelwal et al. [160] FEM-Q9-HZZT (Equil) 0.3399 0.3187 5.6329
Chakrabarti et Sheikh [158] FEM-T6-HZZT 0.2197 - 5.6079
3.3.4.5 Comparaison avec des résultats expérimentaux sur des plaques sandwichs
rectangulaires ayant des peaux stratifiées de type CFRP
97
Chapitre 3 Développement d’un nouvel élément fini rectangulaire (RSFT52)
constitués de plis croisés en fibre de carbone-époxy (CFRP)32, tandis que l’âme est un
matériau en nids-d’abeilles en aluminium (Aluminum Honeycomb Core). Les spécifications
et les propriétés mécaniques des matériaux utilisés dans l'analyse sont présentées,
respectivement dans les tableaux 3.14 à 3.17. Quatre types de séquence d'empilement dont les
désignations sont SPl, SP2, SP3 et SP4 ont été adoptés pour les faces de la plaque sandwich
(figure3.11). L'épaisseur de chaque couche est de 0,125 mm, tandis que l’épaisseur de l’âme
est de 10 mm pour les cas SP1 et SP2, et de 7 mm pour les cas SP3 et SP4. La plaque est
soumise à un chargement réparti uniforme d'intensité P = 1,01 KPa. Dans ces tests, les auteurs
ont contrôlé la pression appliquée par un manomètre, et la déflexion au centre de la plaque a
été mesurée en utilisant des transformateurs différentiels à variation linéaire (LVDT)33. Une
technique holographique a été utilisée afin de visualiser les modes de flexion.
Z Z Z Z
30 0
30 0 30 0
30 0 30 90
30 0
X X X X
30 0 30 0
30 0 30 90
30 0 30 0
SP1 SP 2 SP3 SP 4
Par ailleurs, les auteurs ont fourni des solutions analytiques, basées sur la méthode de
Rayleigh–Ritz, pour le même problème de plaque, en utilisant deux types de conditions aux
limites, simplement appuyée (SSSS) et encastrée (CCCC). Les résultats du déplacement
32
Carbon Fiber–Reinforced Plastic
33
Les capteurs de déplacement LVDT (Linear Variable Differential Transformer) ou transformateur différentiel
à variation linéaire sont des capteurs qui servent à déterminer avec une bonne précision la position d’un point.
Les capteurs de déplacement LVDT sont constitués par deux enroulements électromagnétiques indépendants
dont le primaire est alimenté en courant électrique. Lorsque la tige subit un mouvement, celle-ci provoque un
changement de flux magnétique. Ce changement engendre alors un courant électrique dans l’enroulement
secondaire. Ce courant secondaire permet de quantifier la position des points [lien]
98
Chapitre 3 Développement d’un nouvel élément fini rectangulaire (RSFT52)
transversal ( w) , obtenus par le présent élément sont représentés dans le tableau 3.18, en
utilisant un maillage de 16×16. La comparaison a été faite avec des solutions analytiques
données par Kanematsu et al. [273], avec des solutions basées sur des travaux expérimentaux
présentés par les mêmes auteurs [273] et avec les modèles d’éléments finis de Lee et Fan
[243] et Nayak et al. [231]. La comparaison a montré l’efficacité et la fiabilité du modèle
d’élément fini développé (RSFT52) dans les analyses de ce type de structures.
Matériaux Spécification
Peaux (CFRP) TORAYKA Pre-preg S305-125
Âme (A1 Honeycomb) Al 3/8-5056-0007
adhésif REDUX 312 UL/REDUX 540
99
Chapitre 3 Développement d’un nouvel élément fini rectangulaire (RSFT52)
Tableau 3. 18 Déflexion ( w) d’une plaque sandwich rectangulaire ayant des peaux stratifiées, sous
chargement réparti uniforme
Le rapport d’épaisseur (a h)
Le rapport d’aspect (a b)
Les conditions aux limites.
Le degré d'orthotropie de la peau ( E1 E2 ) .
L'analyse est effectuée pour une plaque sandwich symétrique soumise à une charge
sinusoïdale. Le matériau utilisé est un sandwich constitué d’une âme en mousse HEREX C70-
100
Chapitre 3 Développement d’un nouvel élément fini rectangulaire (RSFT52)
130 et deux peaux stratifiées de plis croisés 0 90 0 en fibre de verre / polyester.
L'épaisseur de chaque pli est de 0,125 mm, tandis que l’épaisseur de l’âme est de 10 mm. Les
propriétés mécaniques de ces matériaux sont présentées dans le tableau 3.19. Les valeurs du
Peaux en fibre de
24.51 7.77 3.34 3.34 1.34 0.078 1800
verre / polyester
Âme en mousse de
0.1036 0.1036 0.05 0.05 0.05 0.32 130
HEREX - C70.130
Les courbes de la figure 3.12 représentent l’effet de l'épaisseur de l’âme (hc ) sur le
Les deux courbes de la figure 3.13, montrent l’effet du rapport d’épaisseur a h sur le
d’aspect a b .
101
Chapitre 3 Développement d’un nouvel élément fini rectangulaire (RSFT52)
montrée sur la figure 3.15 pour des valeurs différentes de hc h . À partir de cette figure, on
remarque, pour différentes valeurs de l’épaisseur de l’âme, que la variation du déplacement
Fig. 3. 12 Effet de l'épaisseur de l’âme (hc ) sur le déplacement transversal ( w) d’une plaque
sandwich (0 90 0 C 0 90 0) avec différentes conditions aux limites.
102
Chapitre 3 Développement d’un nouvel élément fini rectangulaire (RSFT52)
103
Chapitre 3 Développement d’un nouvel élément fini rectangulaire (RSFT52)
3.4 Conclusion
Dans ce chapitre, un nouvel élément fini de forme rectangulaire à quatre nœuds avec
treize degrés de liberté par nœud, a été développé. Il est formulé sur la base du modèle
layerwise que nous avons proposé et présenté dans le chapitre précédent. Cet élément est
destiné à l’analyse statique des plaques sandwichs, qu’elles soient symétriques ou non-
symétriques, ainsi que les plaques isotropes et les composites stratifiées multicouches.
Par ailleurs, on a aussi montré à travers une étude paramétrique, l'effet de certains
paramètres, tels que les caractéristiques géométriques et mécaniques, sur la variation du
déplacement transversal d’une plaque sandwich ayant des peaux stratifiées.
105
Chapitre 4 Développement d’un nouvel élément fini isoparamétrique (QSFT52)
Chapitre 4
4.1 Introduction
Par conséquent, pour généraliser notre élément fini RSFT52 dans l’analyse des
structures de formes quelconques, ce chapitre a été consacré au développement d’un nouvel
élément quadrilatéral, nommé QSFT52 (Quadrilateral Sandwich First Third with 52- degrees
of freedom), à quatre nouds et 13 ddl par nœud. Il est construit sur la base du modèle
layerwise, présenté dans le chapitre II. La transformation a été accomplie en introduisant la
notion d’isoparamétrie. En outre, la méthode des déformations naturelles de substitution
106
Chapitre 4 Développement d’un nouvel élément fini isoparamétrique (QSFT52)
(ANS)34 est introduite dans la présente formulation afin d'atténuer d'éventuels phénomènes de
verrouillage en cisaillement.
4.2.1 Définition
Le contenu théorique de cette partie est issu d’une petite synthèse des références
consultées par: Belarbi [68] et Tati [275].
L’utilisation des éléments finis isoparamétriques dont les côtés peuvent être courbes
(appelés “éléments distordus”), ont été introduits pour la première fois par Ergatoudis et al.
[276]. Le concept d’isoparamétrie permet à la fois de systématiser l’utilisation des variables
réduites, et de générer aussi des éléments finis à bords courbes; il est donc très avantageux
pour certaines discrétisations. Ce concept est basé sur la transformation entre le domaine
physique de l’élément et un domaine simple sans dimensions, en utilisant les coordonnées
intrinsèques. Cette transformation permet de créer facilement les fonctions d’interpolation des
déplacements. Donc, la construction d’un élément fini fait intervenir deux types
d’interpolations: l’interpolation géométrique et l’interpolation des déplacements. Un élément
est dit iso-paramétrique si les interpolations de la géométrie et des déplacements sont
identiques. Les propriétés de l’élément isoparamétrique distordu pourront être établies
moyennant la connaissance :
La figure 4.1 représente un élément parent carré ainsi qu’un élément distordu qui en
dérive.
34
Assumed Natural Strains method
107
Chapitre 4 Développement d’un nouvel élément fini isoparamétrique (QSFT52)
ii) La transformation doit être définie uniquement par les coordonnées xi , yi des 4
nœuds de l’élément distordu.
iii) Si le point P’ coïncide avec un nœud de l’élément parent, le point P correspondant est
le nœud de même numéro de l’élément distordu.
x ξ
f
y η (4.1)
108
Chapitre 4 Développement d’un nouvel élément fini isoparamétrique (QSFT52)
Représentation de la géométrie
4
x , Ni , xi
i 1
(Interpolation géométrique) (4.2)
4
y , N i , yi
i 1
4
u , Ni , ui
i 1
(Interpolation des déplacements) (4.3)
4
v , N i , vi
i 1
où:
1
N1 , 1 1
4
1
N 2 , 1 1
4 (4.4)
1
N 3 , 1 1
4
1
N 4 , 1 1
4
109
Chapitre 4 Développement d’un nouvel élément fini isoparamétrique (QSFT52)
Les dérivées des fonctions d'interpolation géométriques seront calculées par la formule
suivante:
Ni N i N i
x x x
(4.5)
Ni N i N i
y y y
Ni Ni
x x x
N
Ni
(4.6)
i
y y y
Les dérivées , , et se déterminent à partir de la matrice Jacobienne inversée
x y x y
x y
J x y
(4.7)
Ni Ni
4
xi yi
J N Ni
(4.8)
i 1 i
xi yi
1 1
dxdy det J d d
1 1
(4.9)
110
Chapitre 4 Développement d’un nouvel élément fini isoparamétrique (QSFT52)
4.3.1 Description
111
Chapitre 4 Développement d’un nouvel élément fini isoparamétrique (QSFT52)
4
N i i ,i i (4.10)
i 1
où i ui vi wi xic yic xic yic xic yic xit yit xib yib est le vecteur des déplacements
T
L’approximation des variables cinématiques de trois couches peut être exprimée comme suit:
b. L’âme:
4 4
xc , N i , xic , yc , N i , yic
i 1 i 1
4 4
c
x , N i , c
xi
, c
y , N i , yic
i 1 i 1
c. Peau supérieure:
4 4
d. Peau inférieure:
4 4
, Ni ,
b
x
b
xi
, , Ni , yib
b
y
(4.14)
i 1 i 1
B B
e e e
B i e ;
1 e 1 2 e 2 3 e 3
i e
; i e
(4.15)
0 e
B B
i e ;
e e
B
0 e 1 e 1 2 e 2
s s s s i e
; s s i e
112
Chapitre 4 Développement d’un nouvel élément fini isoparamétrique (QSFT52)
Où les matrices de déformation B1 , B 2 , B3 , B 0 , B1 et B 2 peuvent être
s s s
N i
0 0 0 0 0 0 0 0 0 0 0 0
x
1 N i
B 0 0 0 0 0 0 0 0 0 0 0 0
y
N i N i
0 0 0 x y
0 0 0 0 0 0 0 0
N i
0 0 0 0 0 0 0 0 0 0 0 0
x
2 N i
B 0 0 0 0 0 0 0 0 0 0 0 0
y
N i N i
0 0 0 0 0 x y
0 0 0 0 0 0
N i
0 0 0 0 0 0 0 0 0 0 0 0
x
N i
B 3 0 0 0 0 0 0 0 0 0 0 0 0 (4.16)
y
N i N i
0 0 0 0 0 0 0 y x
0 0 0 0
N i
0 0 x
Ni 0 0 0 0 0 0 0 0 0
0
B
0 0 N i
0 0 0 0 0 0 0 0
s
0 Ni
y
0 0 0 0 0 2 Ni 0 0 0 0 0 0 0
B1
s
0 0 0 0 0 0 2 Ni 0 0 0 0 0 0
0 0 0 0 0 0 0 3Ni 0 0 0 0 0
B 2
s
0 0 0 0 0 0 0 0 3 Ni 0 0 0 0
113
Chapitre 4 Développement d’un nouvel élément fini isoparamétrique (QSFT52)
d’où
N i N i N i N i
0 0 0 c 0 d 0 e 0 f 0 0 0
x x x x
N i N i N i N i
Bmt 0 0 0 0 c 0 d 0 e 0 f 0 0
y y y y
N i N i N i N i N i N i N i N i
0 0 0 c y c
x
d
y
d
x
e
y
e
x
f
y
f
x
0 0
avec : (4.18)
2 3
hc h h hc
c , d c
, e c
, f
2 4 8 2
N i
0 0 0 0 0 0 0 0 0 0 0 0
x
N i
B tf 0 0 0 0 0 0 0 0 0 0 0 0
y
N i N i
0 0 0 0 0 0 0 0 0 y x
0 0
N i
0 0 x
0 0 0 0 0 0 Ni 0 0 0
Bc
t
0 0 N i
0 0 0 0 0 0 0 Ni 0 0
y
Même façon pour la peau inférieure, les matrices déformation-déplacement généralisées sont
données par :
d’où,
114
Chapitre 4 Développement d’un nouvel élément fini isoparamétrique (QSFT52)
N i N i N i N i
0 0 0 f 0 d 0 h 0 0 0 c 0
x x x x
N i N i N i N i
Bm 0 0 0
b
0 f 0 d 0 h 0 0 0 c
y y y y
N i N i N i N i N i N i N i N i
0 0 0 f
y
f
x
d
y
d
x
h
y
h
x
0 0 c
y
c
x
avec :
hc hc2 hc3 hc
f , d , h , c
2 4 8 2
N i
0 0 0 0 0 0 0 0 0 0 0 0
x
N i
B f 0 0 0 0 0 0 0 0 0 0 0
b
0
y
N i N i
0 0 0 0 0 0 0 0 0 0 0 y x
(4.20)
N i
0 0 x
0 0 0 0 0 0 0 0 Ni 0
Bc
b
0 0 N i
0 0 0 0 0 0 0 0 0 Ni
y
Dans l’analyse des plaques minces, en utilisant des théories tenant compte de l’effet de
cisaillement transverse, un phénomène appelé blocage ou verrouillage, en cisaillement
("Locking") pourrait survenir. La raison de ce blocage est que l'élément, lorsqu'il est chargé
en flexion pure, présentera une surestimation de la rigidité en cisaillement, i.e. le
comportement en cisaillement devient prédominant devant les autres comportements. Afin de
remédier à ce problème, Dvorkin and Bathe [234] et Huang and Hinton [235] ont proposé la
méthode dite: ''méthode des déformations naturelles de substitution (ANS)'' comme une
solution pour l’éviter ou du moins pour en limiter ses effets (pour plus de détails, le lecteur
pourra consulter les articles [231, 239, 277, 278]).
Dans ce chapitre, nous allons introduire cette méthode, dans les deux peaux, pour
atténuer l’éventuel phénomène de verrouillage en cisaillement. Les déformations naturelles de
substitution sont dérivées en utilisant des fonctions d'interpolation, basées sur le polynôme de
115
Chapitre 4 Développement d’un nouvel élément fini isoparamétrique (QSFT52)
Lagrange et les valeurs de déformation angulaire aux points (sampling points) où il n'existe
pas de verrouillage [239].
f A
Fig. 4.3 Position du point de prélèvement pour les composantes: xz
f A
(à gauche) et yz (à
droite)
Où indique la position du point de prélèvement comme illustré dans la figure 4.3. Les
fonctions P et Q sont interpolés d’une manière linéaire, deux à deux, comme suit:
1 1
P1 1 , P2 1
2 2 (4.23)
1 1
Q1 1 , Q2 1
2 2
t A
4
s B s i
t
i 1
(4.24)
B b A
4
b
s
i 1
s
i
116
Chapitre 4 Développement d’un nouvel élément fini isoparamétrique (QSFT52)
t A b A
où B s et B s sont les matrices de déformations naturelles de substitution de peau
supérieure et inférieure, respectivement.
U U c +U t U b (4.25)
où
U c c dv t t dv b b dv
T T T
v v v
(4.26)
Le travail virtuel effectué par les forces externes est donnés par:
1 1
W f ( x, y) w det J d d
1 1
(4.27)
117
Chapitre 4 Développement d’un nouvel élément fini isoparamétrique (QSFT52)
1 1
det J d d f w det J d d 0
1 1
T
2 2
s F s
s
1 1
Les composantes de la matrice de rigidité et des vecteurs forces, de l’âme, sont évaluées
par intégration numérique en utilisant (2x2) points de Gauss, soit:
1 1
B 0 T A B 0 B 0 T B B 1 B 0 T D B 2
K c
1 1
s s s s s s
s s s s s s
Les mêmes étapes ont été suivies pour élaborer la matrice de rigidité globale de deux peaux,
on aura donc:
a. Peau supérieure
1 1
K t B t T A t B t B t T B t B t B t T B t B t
m m m f f m
1 1
membranaire couplage membrane-flexion couplage flexion-membrane
(4.32)
+ B tf D B tf Bct Ac Bct det J d d
T t T t
flexionnelle cisaillement
b. Peau inférieure
1 1
1 1
membranaire couplage membrane-flexion couplage flexion-membrane
(4.33)
+ B f D B f Bc Ac Bc det J d d
b T b b b T b b
flexionnelle cisaillement
La matrice de rigidité globale s’obtient par l’assemblage des matrices de rigidité élémentaires.
Le but de cette partie est de montrer, à travers une série d’applications, la performance
de l’élément QSFT52 dans la résolution des problèmes de flexion des plaques composites
stratifiées et sandwichs, de forme rectangulaires et biaises (skew plates). Cet élément permet
d’évaluer le déplacement transversal, les contraintes de cisaillement transversal ainsi que les
contraintes normales. Dans ces tests, différents cas de chargements, de géométries et de
conditions aux limites ont été considérés. Les résultats numériques obtenus sont comparés
119
Chapitre 4 Développement d’un nouvel élément fini isoparamétrique (QSFT52)
avec les solutions analytiques et ceux obtenus par d’autres modèles d’éléments finis
disponibles dans la littérature.
dont les caractéristiques mécaniques sont les mêmes que celles adoptées dans l'exemple
3.3.3.1 (chapitre III). Dans la présente analyse, trois angles de biaises 30 , 45 et 60 avec
deux rapports d'épaisseur h / a 0.1 et 0.2 sont considérés.
120
Chapitre 4 Développement d’un nouvel élément fini isoparamétrique (QSFT52)
121
Chapitre 4 Développement d’un nouvel élément fini isoparamétrique (QSFT52)
doublement sinusoïdal
x , y , xy h2
q0 a 2
x , y , xy (4.35)
xz
, yz
h
q0 a
xz , yz (4.36)
100 E 2 h3 w
w 4 (4.37)
a q0
Les résultats obtenus sont rapportés dans le tableau 4.2 pour différents types de
maillage. Ces derniers sont comparés avec la solution d’élasticité-3D donnée par Pagano [24],
d’une part, et d’autre part avec des modèles d’éléments finis basés sur des théories différentes,
connues dans la littérature scientifique (MFEM-3D-LW, FEM-Q9-TSDT, FEM-Q9-IHZZT,
FEM-Q8-HZZT, FEM-Q9-HZZT, FEM-Q9-HSDT):
MFEM-3D-LW: Élément fini mixte 3D, basé sur l’approche layerwise, présenté par
Ramtekkar et al. [246]
FEM-Q9-TSDT: Élément fini quadrilatère à neuf nœuds, basé sur la théorie du
troisième ordre, présenté par Tu et al. [236]
FEM-Q9-IHZZT: Élément fini quadrilatère à neuf nœuds, basé sur la théorie de zig-
zag d’ordre supérieur amélioré, présenté par Pandit et al. [162]
FEM-Q8-HZZT: Élément fini quadrilatère à huit nœuds, basé sur la théorie de zig-zag
d’ordre supérieur, présenté par Singh et al. [261]
122
Chapitre 4 Développement d’un nouvel élément fini isoparamétrique (QSFT52)
FEM-Q9-HZZT: Élément fini quadrilatère à neuf nœuds, basé sur la théorie de zig-zag
d’ordre supérieur , présenté par Khandelwal et al. [160]
FEM-Q9-HSDT: Élément fini quadrilatère à neuf nœuds, basé sur la théorie d’ordre
supérieur, présenté par Manjunatha et Kant [227]
Les figures 4.5 et 4.6, illustrent les distributions de contraintes planes, suivant
l'épaisseur de la plaque sandwich, pour les deux rapports d’épaisseur h a 0.1 et 0.25 ,
respectivement.
Le tableau 4.2 et les figures 4.5 et 4.6, montrent que l’élément développé QSFT52
donne des résultats très satisfaisants par rapport aux valeurs de référence, spécialement en ce
qui concerne les valeurs de contraintes de cisaillement transversal.
Cet exemple a notamment été traité par Pandit [162]. Il s’agit d’une plaque sandwich
carrée encastrée dans les quatres côtés et soumise à une charge uniformément répartie. Les
mêmes propriétés géométriques et mécaniques du problème précédent ont été adoptées.
L’étude a été faite pour différents rapports d’épaisseur a / h 100, 50, 20, 10 et 4 avec un
rapport d’épaisseur a h ont été présentées sur la figure 4.7. D'après cette figure, on observe
une excellente corrélation entre les résultats obtenus par l’élément développé (QSFT52) et à
ceux obtenus par l’élément de Pandit et al. [162]. On peut noter que l’élément développé est
applicable aux problèmes de flexion des plaques sandwichs qu’elles soient minces ou épaisses
(aucun phénomène de verrouillage en cisaillement).
123
Chapitre 4 Développement d’un nouvel élément fini isoparamétrique (QSFT52)
Tableau. 4. 2 Déplacement transversal ( w) , contraintes planes xx , yy , xy et contraintes de cisaillement transversal xz , yz , d’une plaque sandwich
symétrique ( f / c / f ) , simplement appuyée, sous chargement doublement sinusoïdal
a b *
xx , , * yy , , * xz 0, ,0 * yz , 0, 0 * xy 0,0, *
h a b h a b h b a h
Références Modèles EF w , ,0
a 2 2 2 2 2 2 2 2 2 2 2
0.5 Présent élément (8×8) QSFT52 23.4200 2.2463 0.3931 0.1766 0.1363 0.2283
Présent élément (10×10) QSFT52 23.3493 2.2886 0.4014 0.1803 0.1388 0.2330
Présent élément (12×12) QSFT52 23.3106 2.3119 0.4060 0.1823 0.1401 0.2356
Présent élément (14×14) QSFT52 23.2875 2.3260 0.4088 0.1835 0.1409 0.2372
Présent élément (16×16) QSFT52 23.2725 2.3351 0.4106 0.1843 0.1414 0.2382
Pagano [24] Elasticity solution* 21.6531 2.6530 0.3920 0.1850 0.1400 0.2340
Ramtekkar et al. [246] MFEM-3D-LW - 2.6840 0.3960 0.1860 0.1420 0.2360
Kant et Kommineni [225]** FEM-Q9-TSDT 21.3707 2.7985 - - - 0.2371
0.25 Présent élément (8×8) QSFT52 7.7651 1.3281 0.2399 0.2255 0.1038 0.1386
Présent élément (10×10) QSFT52 7.7526 1.3563 0.2452 0.2303 0.1052 0.1417
Présent élément (12×12) QSFT52 7.7457 1.3719 0.2481 0.2330 0.1059 0.1434
Présent élément (14×14) QSFT52 7.7416 1.3813 0.2499 0.2346 0.1064 0.1444
Présent élément (16×16) QSFT52 7.7388 1.3875 0.2511 0.2357 0.1067 0.1451
Pagano [24] Elasticity solution 7.5962 1.5160 0.2595 0.2390 0.1072 0.1440
Pandit et al. [162] FEM-Q9-IHZZT 7.6552 1.5218 0.2506 0.2520 0.1156 0.1468
Tu et al. [236] FEM-Q9-TSDT 7.5610 1.5518 0.2483 0.2447 0.1184 0.1459
Singh et al. [261] FEM-Q8-HZZT 7.8556 1.5480 - 0.2611 - 0.1671
Khandelwal et al. [160] FEM-Q9-HZZT 7.5873 1.5316 0.2674 0.2538 0.1192 -
Chalak et al. [262] FEM-Q9-HZZT 7.5822 1.5306 0.2581 0.2436 0.1147 0.1445
Ramtekkar et al. [246] MFEM-3D-LW - 1.5700 0.2600 0.2400 0.1080 0.1450
Wu et Lin [242] MFEM-3D-LW - 1.5480 0.2413 0.2497 - 0.1339
Pandya et Kant [223] FEM-Q9-HSDT 0.6947 1.5230 0.2414 0.2750 - 0.1419
Manjunatha et Kant [227] FEM-Q9-HSDT 7.1596 - - 0.2370 0.1040 -
Kant et Kommineni [225] FEM-Q9-TSDT 7.1502 1.4989 - - - 0.1428
Kant et Swaminathan [30] HSDT-Anal 7.0551 1.5137 0.2648 - - 0.1379
124
Chapitre 4 Développement d’un nouvel élément fini isoparamétrique (QSFT52)
a b
xx , , * yy , , * xz 0, ,0 * yz , 0, 0 * xy 0,0, *
h a b h a b h b a h
Références Modèles EF w , , 0 *
a 2 2 2 2 2 2 2 2 2 2 2
0.2 Présent élément (8×8) QSFT52 5.5650 1.2038 0.1972 0.2424 0.0893 0.1157
Présent élément (10×10) QSFT52 5.5600 1.2311 0.2017 0.2476 0.0903 0.1183
Présent élément (12×12) QSFT52 5.5571 1.2461 0.2042 0.2505 0.0908 0.1197
Présent élément (14×14) QSFT52 5.5554 1.2552 0.2057 0.2523 0.0911 0.1206
Présent élément (16×16) QSFT52 5.5554 1.2612 0.2067 0.2535 0.0913 0.1212
Pagano [24] Elasticity solution* 5.4746 1.3704 0.2094 0.2569 0.0918 -
Khandelwal et al. [160] FEM-Q9-HZZT 5.4464 1.3617 0.2216 0.2530 0.1025 -
0.1 Présent élément (8×8) QSFT52 2.1964 1.0484 0.1020 0.2815 0.0538 0.0662
Présent élément (10×10) QSFT52 2.2036 1.0777 0.1047 0.2880 0.0534 0.0680
Présent élément (12×12) QSFT52 2.2075 1.0939 0.1062 0.2916 0.0532 0.0690
Présent élément (14×14) QSFT52 2.2099 1.1038 0.1071 0.2938 0.0531 0.0696
Présent élément (16×16) QSFT52 2.2115 1.1103 0.1077 0.2952 0.0530 0.0699
Pagano [24] Elasticity solution 2.2004 1.1531 0.1104 0.3000 0.0530 0.0707
Pandit et al. [162] FEM-Q9-HZZT 2.2002 1.1483 0.1086 0.3158 0.0570 0.0709
Tu et al. [236] FEM-Q9-TSDT 2.2027 1.1466 0.1105 0.3181 0.0532 0.0715
Singh et al. [261] FEM-Q8-HZZT 2.2389 1.1594 - 0.3237 - 0.0759
Khandelwal et al. [160] FEM-Q9-HZZT 2.1786 1.1539 0.1184 0.3185 0.0598 -
Chalak et al. [262] FEM-Q9-HZZT 2.1775 1.1528 0.1143 0.3058 0.0575 0.0705
Ramtekkar et al. [246] MFEM-3D-LW - 1.1590 0.1110 0.3030 0.0550 0.0720
Wu et Lin [242] MFEM-3D-LW - 1.2100 0.1115 0.3240 - 0.0713
Pandya et Kant [223] FEM-Q9-HSDT 0.2023 1.1660 0.1052 0.3400 - 0.0692
Kant et Kommineni [225]** FEM-Q9-TSDT 2.0864 1.1657 - - - 0.0692
Kant et Swaminathan [30] HSDT-Anal 2.0798 1.1523 0.1100 - - 0.0685
Nayak et al. [231] FEM-Q4-HSDT - 1.1410 0.1034 0.3465 0.0574 0.0685
Nayak et al. [231] FEM-Q9-HSDT - 1.1510 0.1043 0.3506 0.0580 0.0689
125
Chapitre 4 Développement d’un nouvel élément fini isoparamétrique (QSFT52)
a b
xx , , yy , , * xz 0, ,0 * yz , 0, 0 * xy 0,0, *
h a b h a b h b a h
Références Modèles EF w , , 0
a 2 2 2 2 2 2 2 2 2 2 2
0.05 Présent élément (8×8) QSFT52 1.1766 0.9895 0.0627 0.2944 0.0430 0.0459
Présent élément (10×10) QSFT52 1.1953 1.0292 0.0652 0.3025 0.0407 0.0477
Présent élément (12×12) QSFT52 1.2058 1.0516 0.0666 0.3069 0.0393 0.0487
Présent élément (14×14) QSFT52 1.2122 1.0653 0.0674 0.3097 0.0385 0.0493
Présent élément (16×16) QSFT52 1.2164 1.0743 0.0680 0.3115 0.0379 0.0497
Pagano [24] Elasticity solution* 1.2264 1.1100 0.0700 0.3174 0.0361 0.0511
Pandit et al. [162] FEM-Q9-HZZT 1.2254 1.1055 0.0694 0.3342 0.0392 0.0509
Singh et al. [261] FEM-Q8-HZZT 1.2424 1.1161 - 0.3429 - 0.0536
Khandelwal et al. [160] FEM-Q9-HZZT 1.2128 1.1113 0.0769 0.3374 0.0415 -
Chalak et al. [262] FEM-Q9-HZZT 1.2121 1.1103 0.0742 0.3272 0.0399 0.0508
Ramtekkar et al. [246] MFEM-3D-LW - 1.1150 0.0700 0.3170 0.0360 0.0510
Wu et Lin [242] MFEM-3D-LW - 1.1730 0.0724 0.3530 - 0.0525
Kant et Kommineni [225]** FEM-Q9-TSDT 1.1947 1.1246 - - - 0.0506
Kant et Swaminathan [30] HSDT-Anal 1.1933 1.1110 0.0705 - - 0.0504
0.02 Présent élément (8×8) QSFT52 0.7347 0.8099 0.0419 0.2801 0.0646 0.0329
Présent élément (10×10) QSFT52 0.7962 0.8979 0.0465 0.2937 0.0547 0.0365
Présent élément (12×12) QSFT52 0.8341 0.9523 0.0493 0.3019 0.0483 0.0387
Présent élément (14×14) QSFT52 0.8588 0.9877 0.0511 0.3072 0.0441 0.0401
Présent élément (16×16) QSFT52 0.8756 1.0119 0.0524 0.3107 0.0412 0.0411
Pagano [24] Elasticity solution 0.9348 1.0990 0.0569 0.3230 0.0306 0.0446
Pandit et al. [162] FEM-Q9-HZZT 0.9341 1.0948 0.0566 0.3403 0.0333 0.0445
Singh et al. [261] FEM-Q8-HZZT 0.9458 1.1050 - 0.3617 - 0.0465
Chalak et al. [262] FEM-Q9-HZZT 0.9248 1.0997 0.0611 0.3300 0.0321 0.0443
126
Chapitre 4 Développement d’un nouvel élément fini isoparamétrique (QSFT52)
127
Chapitre 4 Développement d’un nouvel élément fini isoparamétrique (QSFT52)
128
Chapitre 4 Développement d’un nouvel élément fini isoparamétrique (QSFT52)
sinusoïdale (h a 0.25)
Tableau 4.3 Déplacement transversal ( w) , contraintes planes xx , yy , xy
et contraintes de cisaillement transversal xz , yz , d’une plaque sandwich
symétrique ( f / c / f ) , encastrée, sous chargement uniforme
a b yz
xx , , yy , , xz 0, ,0 xy 0,0,
a a b h a b h b a h
Références Modèles EF w , , z , 0, 0
h 2 2 2 2 2 2 2 2 2 2 2
130
Chapitre 4 Développement d’un nouvel élément fini isoparamétrique (QSFT52)
propriétés mécaniques sont représentées dans le tableau 4.4. L’épaisseur totale de la plaque
(h) est de 20.4724 mm, où l’épaisseur de chaque peau est de 0.7112 mm. Deux types de
matériaux avec deux types de conditions aux limites (SSSS et CCCC) ont été étudiés. Pour la
modélisation, trois types de maillage ont été utilisés (8×8, 12×12 et 16×16).
élément, sont comparés à ceux obtenus par la solution analytique de Azar [22] (pour le
matériau I) et ceux obtenus à partir des modèles d’éléments finis rencontrés dans la littérature
basés sur la théorie de zig-zag (tableau 4.5). Les résultats obtenus sont en très bon accord avec
les solutions de références de Khatua et Cheung [210], Topdar et al. [163] et Chakrabarti et
Sheikh [260]. Ces derniers confirment la fiabilité et la robustesse de l’élément développé
QSFT52 en termes de convergence et de stabilité.
4.4.5 Plaque sandwich carrée avec des peaux stratifiées alternés (angle-ply)
On considère pour cette étude, une plaque sandwich carrée ayant des peaux stratifiées
131
Chapitre 4 Développement d’un nouvel élément fini isoparamétrique (QSFT52)
de la plaque. Les propriétés mécaniques des matériaux utilisés pour l’âme et les deux peaux
sont les mêmes que celles adoptées dans l'exemple 3.3.4.4 (chapitre III). L'analyse est
effectuée pour trois rapports d'épaisseur h a 0.05, 0.1 et 0.2 et trois angles d’orientation<
(0 , 30 et 45 ) .
Dans le tableau 4.6, on présente les résultats du déplacement transversal, les contraintes
planes au centre de la plaque, et les contraintes de cisaillement transversal aux milieux des
côtés. Les résultats obtenus par le présent élément (QSFT52) sont comparés à ceux obtenus
avec des modèles d’éléments finis de Chakrabarti et Sheikh [158] et Khandelwal et al. [160],
en utilisant la théorie de zig-zag d’ordre supérieur (FEM-Q9-HZZT, FEM-T6-HZZT). D’après
les résultats qui figurent sur le tableau 4.5, on remarque que l’élément développé QSFT52
donne des résultats très satisfaisants et totalement comparables avec les résultats de
références.
4.5 Conclusion
Dans ce qui suit, l’élément QSFT52 sera utilisé dans l’analyse dynamique des plaques
qu’elles soient isotropes, composites stratifiées ou sandwichs. Et cela fait l’objet du chapitre
suivant.
132
Chapitre 4 Développement d’un nouvel élément fini isoparamétrique (QSFT52)
Tableau 4.5 Déplacement transversal ( w) , contraintes planes xx , yy , xy
et contraintes de cisaillement transversal xz , yz , d’une plaque sandwich
symétrique ( f / c / f ) , sous chargement uniforme
a b h a b h a b h b
yz , 0, 0
a h
Références Modèles EF w , , 100 xx , , yy , , xz 0, ,0 xy 0,0,
2 2 2 2 2 2 2 2 2 2 2
2
Cas I Présent élément (8×8) QSFT52 1.8679 1.5146 1.5146 0.0244 0.0244 -0.9700
Présent élément (12×12) QSFT52 1.8745 1.5588 1.5588 0.0262 0.0262 -1.0396
Présent élément (16×16) QSFT52 1.8767 1.5742 1.5742 0.0270 0.0270 -1.0704
Khatua et Cheung [210] FEM solution 1.8697 - - - - -
Azar [22] Analytical solution* 1.8780 - - - - -
Chakrabarti et Sheikh [260]** FEM-T6-HZZT 1.8750 1.5817 1.5817 0.0284 - -
Chakrabarti et Sheikh [260] FEM-T6-RFSDT 1.8682 1.5816 1.5816 0.0297 - -
Chakrabarti et Sheikh [260] FEM-T6-HSDT 1.7233 1.5766 1.5766 0.0045 - -
Topdar et al. [163] FEM-Q4-RHZZT 1.8793 - - 0.0232 0.0232 -
Cas II Présent élément (8×8) QSFT52 0.8813 0.7906 0.7906 0.0221 0.0221 -0.1802
Présent élément (12×12) QSFT52 0.8750 0.8149 0.8149 0.0246 0.0246 -0.1762
Présent élément (16×16) QSFT52 0.8717 0.8141 0.8141 0.0252 0.0252 -0.1664
Khatua et Cheung [210] FEM solution 0.8707 - - - - -
Chakrabarti et Sheikh [260] FEM-T6-HZZT 0.9535 0.8916 0.8916 0.0558 - -
Chakrabarti et Sheikh [260] FEM-T6-RFSDT 0.8880 0.8225 0.8225 0.0475 - -
Chakrabarti et Sheikh [260] FEM-T6-HSDT 0.7779 0.8076 0.8076 0.0106 - -
Folie [272] FEM solution 0.8814 - - - - -
133
Chapitre 4 Développement d’un nouvel élément fini isoparamétrique (QSFT52)
Tableau 4.6 Déplacement transversal ( w) , contraintes planes xx , yy , xy et contraintes de cisaillement transversal xz , yz , d’une plaque sandwich
carrée avec des peaux stratifiées f ( 90) c f ( 90) , simplement appuyée et soumise un chargement uniforme
a b
xx , , yy , , a
h a b h a b h b
Références Modèles EF w , ,0 xz 0, , 0.4h yz ,0, 0.4h
a 2 2 2 2 2 2 2 2 2 2
0
0.05 Présent élément QSFT52 1.6766 1.6197 1.6197 0.3408 0.3408
Pagano [24] Elasticity solution 1.6957 1.5851 1.5692 0.3560 0.3563
Khandelwal et al. [160] FEM-Q9-HZZT (Const)× 1.7107 1.6265 1.5653 0.4515 0.5230
Khandelwal et al. [160] FEM-Q9-HZZT (Equil)×× 1.7107 1.6265 1.5653 0.3663 0.3612
Chakrabarti et Sheikh [158] FEM-T6-HZZT 1.7126 1.5988 - 0.3732 -
0.1 Présent élément QSFT52 2.6149 1.7320 1.7320 0.3337 0.3337
Pagano [24] Elasticity solution 2.6168 1.6004 1.5794 0.3496 0.3503
Khandelwal et al. [160] FEM-Q9-HZZT (Const) 2.6295 1.6537 1.5829 0.4142 0.4142
Khandelwal et al. [160] FEM-Q9-HZZT (Equil) 2.6295 1.6537 1.5829 0.3489 0.3397
Chakrabarti et Sheikh [158] FEM-T6-HZZT 2.6296 1.6249 - 0.3612 -
0.2 Présent élément QSFT52 6.2940 2.1538 2.1538 0.3242 0.3242
Pagano [24] Elasticity solution 6.2981 1.7093 1.7523 0.3436 0.3413
Khandelwal et al. [160] FEM-Q9-HZZT (Const) 6.3001 1.8111 1.7328 0.3982 0.3993
Khandelwal et al. [160] FEM-Q9-HZZT (Equil) 6.3001 1.8111 1.7328 0.3455 0.3332
Chakrabarti et Sheikh [158] FEM-T6-HZZT 6.3016 1.7792 - 0.3482 -
30
0.05 Présent élément QSFT52 1.2199 0.7884 0.7884 0.3465 0.3465
Khandelwal et al. [160] FEM-Q9-HZZT (Const) 1.2450 0.7985 0.7666 0.4637 0.5079
Khandelwal et al. [160] FEM-Q9-HZZT (Equil) 1.2450 0.7985 0.7666 0.3592 0.3480
Chakrabarti et Sheikh [158] FEM-T6-HZZT 1.2381 0.7653 - 0.3603 -
134
Chapitre 4 Développement d’un nouvel élément fini isoparamétrique (QSFT52)
a b
xx , , yy , , a
h a b h a b h b
Références Modèles EF w , ,0 xz 0, , 0.4h yz ,0, 0.4h
a 2 2 2 2 2 2 2 2 2 2
45
0.05 Présent élément QSFT52 1.0671 0.4366 0.4366 0.3215 0.3215
Khandelwal et al. [160] FEM-Q9-HZZT (Const) 1.0773 0.4479 0.4228 0.5309 0.5309
Khandelwal et al. [160] FEM-Q9-HZZT (Equil) 1.0773 0.4479 0.4228 0.3335 0.3335
Chakrabarti et Sheikh [158] FEM-T6-HZZT 1.0615 0.4247 - - -
0.1 Présent élément QSFT52 2.0035 0.4991 0.4991 0.3223 0.3223
Khandelwal et al. [160] FEM-Q9-HZZT (Const) 1.9950 0.4698 0.4422 0.4445 0.3320
Khandelwal et al. [160] FEM-Q9-HZZT (Equil) 1.9950 0.4698 0.4422 0.3375 0.4592
Chakrabarti et Sheikh [158] FEM-T6-HZZT 1.9764 0.4444 - 0.2203 0.3232
PARTIE III
136
Chapitre 5 Vibration libre des plaques sandwichs et des composites stratifiées
Chapitre 5
5.1 Introduction
5.2 Revue des travaux antérieurs sur l'analyse dynamique des plaques multicouches et
des sandwichs
Les premières études sur le comportement dynamique de plaques stratifiées ont été
entreprises vers les années 1960 par Yu [282, 283], puis d’autres chercheurs, dont Srinivas et
137
Chapitre 5 Vibration libre des plaques sandwichs et des composites stratifiées
Rao [284], Chan et Cheung [285], Noor [286], Ibrahim et al. [287], Reddy et Kuppusamy
[288] et Wang et al. [289], qui ont effectué une analyse dynamique sur des plaques
multicouches/sandwichs, en proposant des solutions analytiques, basées sur les équations
d'élasticité tridimensionnelle.
Cependant, les solutions analytiques rencontrées dans les travaux cités ci-dessus, sont
disponibles uniquement pour des cas de géométries et de conditions aux limites simples. Par
conséquent, il est donc naturel de chercher d'autres techniques numériques telles que la
méthode des éléments finis pour résoudre les problèmes dynamiques des structures
composites multicouches et sandwichs ayant des géométries complexes et pour divers
conditions aux limites.
138
Chapitre 5 Vibration libre des plaques sandwichs et des composites stratifiées
De nombreux modèles d’éléments finis, basés sur des diverses théories, ont été proposés
par les chercheurs. Khatua et Cheung [210] sont les premiers à prédire les fréquences
naturelles de plaques sandwichs, en utilisant deux éléments finis de type triangulaire et
rectangulaire. D’autre part, basant sur la théorie d'ordre supérieur de Reddy [78], Nayak et al.
[277] ont élaboré deux éléments finis quadrilatères isoparamétriques C0 (quatre et neuf
nœuds), avec 7 ddl/nœud, afin d’étudier le comportement dynamique en vibration libre des
plaques: isotropes, orthotropes, composites stratifiées et sandwichs. La matrice consistante de
masse a été adoptée dans cette formulation. Ils ont utilisé la méthode ANS pour éviter le
problème de verrouillage en cisaillement ainsi que les modes parasites. Les auteurs ont achevé
leur travail par une étude paramétrique, afin d’étudier l’effet de l’anisotropie, le rapport
d’aspect, le rapport longueur/épaisseur et le nombre de couches, sur la variation des
fréquences naturelles. Par ailleurs, Khare et al. [292] ont utilisé un élément fini
isoparamétrique quadrilatère à continuité C0 pour étudier les vibrations libres des plaques
sandwichs. La formulation de cet élément a été basée sur deux théories d’ordre supérieur,
nommées HOST11 et HOST9. Une année plus tard, Garg et al. [293] ont présenté l’étude des
vibrations libres des plaques biaises (skew) en matériaux composites stratifiés et sandwichs,
en utilisant un élément fini isoparamétrique C0 , formulé sur la base de la théorie du troisième
ordre, où les composants de champ des déplacements membranaires sont des fonctions
cubiques et le déplacement transversal est supposé constant sur l'épaisseur de la plaque.
Basant sur l’approche layerwise, Lee et Fan [243] ont présenté un élément fini
isoparamétrique à 9 nœuds avec 5 ddl/nœud, afin d’étudier la flexion et la vibration des
plaques sandwichs. Cet élément a été formulé sur la base d’un nouveau modèle de type
layerwise, où la théorie du premier ordre est utilisée pour modéliser les deux peaux, tandis
que les déplacements au niveau de l’âme sont exprimés en fonction des déplacements des
deux peaux. En outre, Nabarrete et al. [294] ont publié un travail concernant le
développement d’un nouveau modèle de type layerwise pour l’analyse du comportement en
vibration des plaques sandwichs avec des peaux stratifiées. Dans ce modèle, les peaux sont
modélisées comme des plaques de Riessner-Mindlin, tandis que l’âme est modélisée en
utilisant la théorie du troisième ordre, où le champ de déplacement est cubique pour les
composantes u, v et quadratique pour la composante w, dans le sens de l’épaisseur. En partant
de ce modèle, les auteurs ont construit un élément fini tridimensionnel (3D) contenant 240
degrés de liberté. Citons aussi l’article de Roque et ses collègues [295] qui a été concentré sur
139
Chapitre 5 Vibration libre des plaques sandwichs et des composites stratifiées
Basant sur la théorie de zig-zag d’ordre supérieur, Chakrabarti et Sheikh [298] ont
exploité le potentiel de leur élément, développé précédemment par les auteurs présents [299],
dans des analyses de vibration libre des plaques multicouches et sandwichs avec des peaux
stratifiées. Cet élément est un triangle à six nœuds avec 7 ddl/nœud. Ce dernier a présenté des
résultats satisfaisants, cependant, le seul problème signalé est que cet élément nécessite la
continuité C1 de déplacement transversal aux interfaces. D’autre part, Kulkarni et Kapuria
[300] ont exploité le potentiel des éléments DKQ35 et IDKQ36, développés respectivement par
Batoz et Tahar [301] et Jeyachandrabose et al. [302], pour analyser le comportement
dynamique des plaques stratifiées composites et sandwichs. Ces deux éléments ont été
formulés en se basant sur la théorie de zig-zag du troisième ordre (ZIGT)37. Le principe de
Hamilton est adopté afin de calculer la matrice de masse consistante, la matrice de rigidité
élémentaire et le vecteur de charge. De plus, les auteurs ont tenté d’étudier l'effet des
fonctions d'interpolation des éléments DKQ et IDKQ sur la précision des éléments de type
quadrilatères, basés sur les théories d’ordre supérieur et zig-zag. En 2013, Chalak et al. [303]
35
Discret Kirchhoff Quadrilateral
36
Improved Discret Kirchhoff Quadrilateral
37
Third-order Zig-zag Theory
140
Chapitre 5 Vibration libre des plaques sandwichs et des composites stratifiées
ont élaboré un modèle d’élément fini efficace, à continuité (C0), destiné pour l'analyse
dynamique des plaques sandwichs à âme molle et peaux stratifiées. L’élément est un
quadrilatère isoparamétrique à 9 nœuds et 11ddl/nœud, formulé sur la base de la théorie de
zig-zag d’ordre supérieur (HZZT). Ce modèle a fourni des résultats très satisfaisants pour les
fréquences naturelles par rapport aux autres théories 2D citées en littérature. Une année plus
tard, d’autres scientifiques, dont Sahoo et Singh [304] ont présenté une nouvelle théorie de
type zig-zag trigonométrique (TZZT)38, dans laquelle, un élément fini isoparamétrique à huit
nœuds avec sept variables nodales, a été développé. Ce modèle est destiné pour l’étude des
réponses des structures composites stratifiées et sandwichs, aux vibrations et flambement.
T Tt Tc + Tb (5.1)
où
chaque couche du sandwich (i t , c, b) ; et ( ) est une dérivée seconde par rapport au temps.
..
D’après le champ de déplacement de l’âme qui est défini par l’équation (2.1), l'énergie
cinétique de l’âme peut être exprimée comme suit:
38
Trigonometric Zig-zag Theory
141
Chapitre 5 Vibration libre des plaques sandwichs et des composites stratifiées
Tc I u u
0 0 0 v0 v0 w0 w0 I1 u0 xc xc u0 v0 yc yc v0
0
I 3 u0 xc xc xc xc xc xc u0 v0 yc yc yc yc yc yc v0 (5.3)
I 4 xc xc xc xc xc xc yc yc yc yc yc yc
I 5 xc xc xc xc yc yc yc yc I 6 xc xc yc yc dxdy
Où les moments d’inerties Ii (0,1,2,3,4,5 et 6) de l’âme sont définis par:
hc
I 0 , I1 , I 2 , I 3 , I 4 , I 5 , I 6 c 1, z, z 2 , z 3 , z 4 , z 5 , z 6 dz
2
(5.4)
hc
2
hc c hc2 c hc3 c h
Tt I 0t u0 u0 u0 x u0 x u0 x u0 c xt
0 2 4 8 2
h h2 h3 h
+v0 v0 v0 c yc v0 c yc v0 c yc v0 c yt w0 w0
2 4 8 2
hc c hc2 c c hc3 c c hc4 c c hc2 c t
+I 0t x u0 x x x x x x x x
2 4 8 16 4
hc c hc2 c c hc3 c c hc4 c c hc2 c t
+ y v0 y y y y y y y y
2 4 8 16 4
h
2
h 3
h
4
h 5
h
3
+I 0t c xc u0 c xc xc c xc xc c xc xc c xc xt
4 8 16 32 8
hc2 c hc3 c c hc4 c c hc5 c c hc3 c t
+ y v0 y y y y y y y y
4 8 16 32 8
h3 h4 h5 h6 h4
+I 0t c xc u0 c xc xc c xc xc c xc xc c xc xt
8 16 32 64 16
h3 h4 h5 h6 h4 (5.5)
+ c yc v0 c yc yc c yc yc c yc yc c yc yt
8 16 32 64 16
142
Chapitre 5 Vibration libre des plaques sandwichs et des composites stratifiées
h h2 h3 h4 h2
+I 0t c xt u0 c xt xc c xt xc c xt xc c xt xt
2 4 8 16 4
h h2 h3 h4 h2
+ c yt v0 c yt yc c yt yc c yt yc c yt yt
2 4 8 16 4
h h2 h3 h
+I1t u0 xt c xc xt c xc xt c xc xt c xt xt
2 4 8 2
hc c t hc c t hc c t hc t t
2 3
+v0 y y y y y y y y
t
y
2 4 8 2
h h2 h3 h
+I1t xt u0 c xt xc c xt xc c xt xc c xt xt
2 4 8 2
h h2 h3 h
+ yt v0 c yt yc c yt yc c yt yc c yt yt
2 4 8 2
+ I 2t xt xt yt yt dxdy
Où les moments d’inerties Iit (0,1 et 2) de la peau supérieure, pour une stratification de n
couches, sont donnés par:
hc
ht k 1
n layer h
D’une manière similaire, l’énergie cinétique de la peau inférieure est exprimée à partir
de l’expression du champ de déplacement (2.11) telle que:
hc c hc2 c hc3 c h
Tb I 0b u0 u0 u0 x u0 x u0 x u0 c xb
0 2 4 8 2
hc c hc2 c hc3 c h
+v0 v0 v0 y v0 y v0 y v0 c yb w0 w0
2 4 8 2
hc c hc2 c c hc3 c c hc4 c c hc2 c b
+I 0b x u0 x x x x x x x x (5.7)
2 4 8 16 4
hc c hc2 c c hc3 c c hc4 c c hc2 c b
+ y v0 y y y y y y y y
2 4 8 16 4
h2 h3 h4 h5 h3
+I 0b c xc u0 c xc xc c xc xc c xc xc c xc xb
4 8 16 32 8
143
Chapitre 5 Vibration libre des plaques sandwichs et des composites stratifiées
h2 h3 h4 h5 h3
+ c yc v0 c yc yc c yc yc c yc yc c yc yb
4 8 16 32 8
h3 h4 h5 h6 h4
+I 0b c xc u0 c xc xc c xc xc c xc xc c xc xb
8 16 32 64 16
hc3 c hc4 c c hc5 c c hc6 c c hc4 c b
+ y v0 y y y y y y y y
8 16 32 64 16
h h2 h3 h4 h2
+I 0b c xb u0 c xb xc c xb xc c xb xc c xb xb
2 4 8 16 4
h h2 h3 h4 h2
+ c yb v0 c yb yc c yb yc c yb yc c yb yb
2 4 8 16 4
h h2 h3 h
+I1b u0 xb c xc xb c xc xb c xc xb c xb xb
2 4 8 2
h h2 h3 h
+v0 yb c yc yb c yc yb c yc yb c yb yb
2 4 8 2
h h2 h3 h
+I1b xb u0 c xb xc c xb xc c xb xc c xb xb
2 4 8 2
h h2 h3 h
+ yb v0 c yb yc c yb yc c yb yc c yb yb
2 4 8 2
+ I 2b xb xb yb yb dxdy
Où les moments d’inerties Iib (0,1 et 2) de la peau inférieure, pour une stratification de n
couches, sont donnés par:
hc
k 1
n layer h
M e N
1 1
m(t ) N N m( c ) N + N m(b ) N det J d d (5.9)
T T T
e 1 1
144
Chapitre 5 Vibration libre des plaques sandwichs et des composites stratifiées
Ni
Ni
Ni
Ni
Ni
Ni
N Ni , i 1 à 4 (5.10)
Ni
Ni
Ni
Ni
Ni
Ni
t2
(U - W - T ) dt 0 (5.11)
t1
Après avoir évalué la matrice de rigidité et la matrice de masse pour tous les éléments,
les équations du mouvement, pour l'analyse de la vibration libre, peuvent être établies sous la
forme de problème aux valeurs propres généralisées comme suit:
145
Chapitre 5 Vibration libre des plaques sandwichs et des composites stratifiées
5.4.1 Plaque isotrope carrée avec différents types de conditions aux limites
été utilisé. Les résultats de la fréquence naturelle ( ) , normalisés par la formule (5.14), sont
présentés dans le tableau 5.1 avec différents types de conditions aux limites.
a (5.14)
G
On peut remarquer que les résultats obtenus par le présent élément QSFT52 sont en bon
accord avec ceux obtenus par Ferreira [305] qui a utilisé un élément quadrilatère à quatre
nœuds formulé sur la base de la théorie du premier ordre.
5.4.2 Plaques composites stratifiées carrées de trois, cinq et neuf couches 0 / 90 / ... / 0
Dans cet exemple, l’effet du nombre de couches (n) ainsi que l’effet de l’anisotropie
( E11 E22 ) sur la fréquence fondamentale (1 ) d’une plaque stratifiée carrée simplement
appuyée sont étudiées. Trois types de stratifications croisées de trois (0 / 90 / 0) , cinq
(0 / 90 / 0 / 90 / 0) et neuf (0 / 90 / 0 / 90 / 0 / 90 / 0 / 90 / 0) couches, en gardant la même
G12 G13
0.6 ,
E22 E22
G23
0.5 ,
E22
12 0.25 , 1.0
146
Chapitre 5 Vibration libre des plaques sandwichs et des composites stratifiées
Tableau 5.1 Fréquences propres ( ) d’une plaque carrée isotrope avec différents types de
conditions aux limites
Conditions Présent élément Ferreira [305]
Modes m N
aux limites QSFT52 FEM-Q4-FSDT
Les résultats de la fréquence fondamentale qui sont obtenus à l'aide d'un maillage de
12×12 sont rapportés dans le tableau 5.2 pour les trois types de stratifications. D’après ce
tableau, on observe une excellente corrélation entre les résultats obtenus par l’élément
développé (QSFT52) et les résultats obtenus par la solution d’élasticité-3D, donnée par Noor
[286], ainsi que ceux obtenus par des modèles d’éléments finis basés sur différentes théories:
147
Chapitre 5 Vibration libre des plaques sandwichs et des composites stratifiées
fondamentale (1 ) d’une plaque stratifiée carrée simplement appuyée, nous avons illustré les
résultats réunis dans le tableau 5.2 sur la figure 5.1. D’après cette dernière, on remarque que
les trois courbes ont pratiquement la même allure et que les valeurs des fréquences
fondamentales augmentent avec l’augmentation du rapport d’anisotropie.
cœur sont:
G23
0.2668 , 12 0.3 et 1.0
E11
h2
100 (5.15)
E11
Les résultats de la première fréquence propre obtenue dans la présente analyse sont
présentés dans le tableau 5.3 avec les solutions d’élasticité 3D de Srinivas et Rao [284] et
celles de Di Sciuva [150]. Les résultats de l’élément QSFT52 sont en bon accord avec les
résultats d’élasticités 3D.
148
Chapitre 5 Vibration libre des plaques composites stratifiées et sandwichs
Tableau 5.2 Effet de l’anisotropie ( E11 E22 ) sur la fréquence fondamentale (1 ) de plaques stratifiée carrées 0 / 90 / ... / 0 simplement appuyée
E11
Références Modèles EF Nbr de couches E22
3 10 20 30
Fig.5.1 Effet du nombre de couches et le rapport d’anisotropie ( E11 E22 ) sur la fréquence
fondamentale (1 ) d’une plaque stratifiée carrée simplement appuyée
Tableau 5.3 Fréquence propre fondamentale (1 ) d’une plaque sandwich carrée f /c/ f
simplement appuyée
5.4.4 Plaque sandwich carrée simplement appuyée ayant des peaux composites
stratifiées symétriques
Dans ce test, deux plaques sandwichs carrées simplement appuyées avec des peaux
sont étudiées. L’épaisseur de l’âme est de 0.8 h , tandis que celle de chaque peau est de 0.1 h
, où h est l'épaisseur totale de la plaque. Les propriétés mécaniques des matériaux utilisés
pour l’âme et les deux peaux sont résumées dans le tableau 5.4. L’étude a été menée pour
150
Chapitre 5 Vibration libre des plaques sandwichs et des composites stratifiées
c
100 a (5.16)
E11 f
Tableau 5.4 Propriétés mécaniques des matériaux utilisés dans l’exemple 5.4.4
Les fréquences propres obtenues pour les six premiers modes sont rapportées dans les
tableaux 5.5 et 5.6. Les résultats obtenus par l’élément QSFT52 sont comparés avec la
solution d’élasticité 3D de Kulkarni et Kapuria [300], les solutions analytiques de Wang et al.
[289], basées sur la méthode de p-Ritz, et avec ceux obtenus par les modèles d’éléments finis
de Chakrabarti et Sheikh [298] et de Kulkarni et Kapuria [300], basés sur la théorie HZZT.
Il apparaît clairement qu’ à partir du tableau 5.5, les résultats de l'élément développé
sont en excellent accord avec les résultats trouvés dans la littérature, à titre d’exemple:
Pour un rapport d’épaisseur a h 6.67 , l’erreur est de 0.0038 % pour la fréquence
conditions aux limites différentes, i.e. CCCC et SCSC. Les fréquences propres non-
dimensionnelles obtenues pour les six premiers modes sont résumées dans le tableau 5.7 avec
différents rapports d’épaisseur (a / h 5, 10 and 20) . On peut observer, en comparaison avec
des modèles d’éléments finis basés sur la théorie HZZT [298, 300, 303, 309], que le présent
élément présente des performances très satisfaisantes.
151
Chapitre 5 Vibration libre des plaques composites stratifiées et sandwichs
Tableau 5.5 Fréquence propre non-dimensionnelle, pour différents modes, d’une plaque sandwich carrée simplement appuyée ayant des peaux composites
stratifiées symétriques 0 / 90 / C / 90 / 0
a Présent élément Kulkarni et Kapuria [300] Kulkarni et Kapuria [300] Wang et al. [289] Chakrabarti et Sheikh [298]
Modes
h QSFT52 3D-Elasticity FEM-Q4-TSDT Anal-Solution FEM-T6-HZZT
152
Chapitre 5 Vibration libre des plaques sandwichs et des composites stratifiées
Tableau 5.6 Fréquence propre non-dimensionnelle, pour différents modes, d’une plaque sandwich carrée simplement appuyée ayant des peaux composites
stratifiées symétriques 0 / 90 / 0 / 90 / 0 / 90 / 0 / 90 / C / Sym
a Présent élément Kulkarni et Kapuria [300] Kulkarni et Kapuria [300] Chakrabarti et Sheikh [298]
Modes
h QSFT52 3D-Elasticity FEM-Q4-TSDT FEM-T6-HZZT
153
Chapitre 5 Vibration libre des plaques sandwichs et des composites stratifiées
Tableau 5.7 Fréquences propres non-dimensionnelles d’une plaque sandwich carrée ayant des peaux composites stratifiées symétriques 0 / 90 / C / 90 / 0
avec différents types de conditions aux limites
a Conditions Présent Khandelwal et al. [309] Chalak et al. [303] Kulkarni et Kapuria [300] Chakrabarti et Sheikh [298]
Modes
h aux limites QSFT52 FEM-Q9-HZZT FEM-Q9-HZZT FEM-Q4-HZZT FEM-T6-HZZT
154
Chapitre 5 Vibration libre des plaques sandwichs et des composites stratifiées
a Conditions Présent Khandelwal et al. [309] Chalak et al. [303] Kulkarni et Kapuria [300] Chakrabarti et Sheikh [298]
Modes
h aux limites QSFT52 FEM-Q9-HZZT FEM-Q9-HZZT FEM-Q4-HZZT FEM-T6-HZZT
155
Chapitre 5 Vibration libre des plaques composites stratifiées et sandwichs
5.4.5 Plaque sandwich carrée simplement appuyée ayant des peaux composites
stratifiées non-symétriques 0 / 90 / C / 0 / 90
On considère une plaque sandwich carrée simplement appuyée ayant des peaux
mécaniques des matériaux utilisés pour l’âme et les deux peaux sont résumées dans le tableau
5.8. Afin d’évaluer la performance de l’élément développé QSFT52 dans l’analyse des
vibrations libres des plaques sandwichs minces ou épaisses, différents rapport d’épaisseur
(a h) variant de 2 à 100 ont été considérés. Le rapport de l’épaisseur de l’âme à l’épaisseur
Tableau 5.8 Propriétés mécaniques des matériaux utilisés dans l’exemple 5.4.5
Les résultats des fréquences propres sont normalisés par la formule 5.17.
b2 c
(5.17)
h E22 f
Les résultats des fréquences propres obtenus par le présent élément sont rapportés dans
le tableau 5.9 avec ceux de la littérature.
156
Chapitre 5 Vibration libre des plaques sandwichs et des composites stratifiées
FEM-Q9-HZZT: Élément fini quadrilatère à neuf nœuds, basé sur la théorie de zig-
zag d’ordre supérieur, présenté par Chalak et al. [303]
Pour les deux rapports d’épaisseur (a h 10 et 100) , les fréquences propres obtenues
pour les six premiers modes sont rapportées dans le tableau 5.10.
fréquence fondamentale sont présentés dans le tableau 5.11. On peut constater, d’après les
tableaux 5.9, 5.10 et 5.11 ainsi que les figures 5.2 et 5.3, que les résultats de l’élément
développé QSFT52 sont en excellent accord avec ceux obtenus par les solutions de
références. On peut également noter que le présent élément se comporte bien, que ce soit en
plaques minces ou épaisses (aucun phénomène de verrouillage en CT).
5.4.6 Vibration libre des plaques sandwichs stratifiées en utilisant une mousse de PVC
dans l’âme:
Dans cet exemple, deux plaques sandwichs stratifiées avec différentes séquences
considérées. Le matériau utilisé est un sandwich constitué d’une âme en mousse PVC
(polychlorure de vinyle) et deux peaux stratifiées en fibre de verre / polyester. Les propriétés
mécaniques de ces matériaux sont présentées dans le tableau 5.13 [274]. L'épaisseur de
chaque pli est de 0,125 mm tandis que l’épaisseur de l’âme est de 10 mm. Les caractéristiques
géométriques de la plaque sont: a h 10 , a b 1 et hc h 0.88 .
157
Chapitre 5 Vibration libre des plaques sandwichs et des composites stratifiées
Les valeurs des fréquences propres sont normalisées par la relation suivante:
a2 c
(5.17)
h Ec
Les résultats des fréquences propres obtenues pour les quatres premiers modes sont
présentés dans le tableau 5.14. La comparaison a été faite avec des solutions analytiques de
type layerwise, présentées par Jam et al. [291], ainsi qu’avec des éléments finis (Q4 et Q9) de
Nayak et al. [277]. Les résultats obtenus par l’élément montrent que les valeurs des
fréquences propres sont comparables avec ceux obtenus via la littérature.
5.4.7 Vibration libre des plaques composites starifiées biaises ( skew plate)
Dans cet exemple, le problème d’une plaque composite starifiée biaise symétrique a été
étudié afin de montrer l’influence de la distorsion géométrique sur les fréquences naturelles.
La plaque est constituée de cinq couches 90 / 0 / 90 / 0 / 90 ayant la même épaisseur. Dans la
présente étude, deux conditions aux limites i.e. SSSS et CCCC, avec différentes angles de
biaises 0 ,15 ,30 et 45 sont considérées. Les caractéristiques mécaniques de chaque
Les valeurs des fréquences propres sont normalisées par la relation suivante:
b2
(5.18)
2
E2
Les résultats des fréquences propres ( ) , obtenues pour les sept premiers modes, sont
rapportés dans le tableau 5.15 pour un rapport d'épaisseur a / h 10 . On remarque que les
résultats obtenus par le présent élément QSFT52 sont également bons en comparaison avec
ceux obtenus via la solution analytique de Wang [313], basée sur la méthode de Rayleigh–
Ritz, ainsi que les deux modèles d’éléments finis (FEM-Q9-FSDT, FEM-Q9-TSDT) de Garg
et al. [293]. Ces résultats sont légèrement dépendant de l'angle d'inclinaison (la valeur de
fréquence propre s’éloigne de celle de la solution 3D avec l'augmentation de l’angle ).
158
Chapitre 5 Vibration libre des plaques composites stratifiées et sandwichs
Tableau 5.9 Fréquence fondamentale non-dimensionnelle (1 ) avec différents rapports d’épaisseurs (a h) pour une plaque sandwich carrée simplement
appuyée ayant des peaux composites stratifiée non-symétriques 0 / 90 / C / 0 / 90
a Présent Rao et al. [310] Sahoo et Singh [304] Chalak et al. [303] Zhen et al. [296] Zhen et al. [296] Kant et Swaminathan [290]
h QSFT52 3D-Elasticity FEM-Q8-TZZT FEM-Q9-HZZT FEM-Q8-GLHSDT GLHSDT-Anal HSDT-Anal
159
Chapitre 5 Vibration libre des plaques sandwichs et des composites stratifiées
Tableau 5.10 Fréquence propre non-dimensionnelle ( ) , pour différents modes, d’une plaque sandwich carrée simplement appuyée ayant des peaux
composites stratifiées non-symétriques 0 / 90 / C / 0 / 90 avec différents rapports d’épaisseurs (a h)
a Présent Rao et Desai [312] Chalak et al. [303] Zhen et al. [296] Zhen et al. [296] Kant et Swaminathan [311]
Modes
h QSFT52 LW-Anal FEM-Q9-HZZT FEM-Q8-GLHSDT GLHSDT-Anal HSDT-Anal
160
Chapitre 5 Vibration libre des plaques sandwichs et des composites stratifiées
Tableau 5.11 Fréquence fondamentale non-dimensionnelle (1 ) avec différents rapports (hc h f ) pour une plaque sandwich carrée simplement appuyée
hc Présent Rao et al. [310] Sahoo et Singh [304] Chalak et al. [303] Zhen et al. [296] Zhen et al. [296] Kant et Swaminathan [290]
hf QSFT52 3D-Elasticity FEM-Q8-TZZT FEM-Q9-HZZT FEM-Q8-GLHSDT GLHSDT-Anal HSDT-Anal
Tableau 5.12 Fréquence fondamentale non-dimensionnelle (1 ) avec différents rapports d’aspects (a b) pour une plaque sandwich carrée simplement
appuyée ayant des peaux composites stratifiées non-symétriques 0 / 90 / C / 0 / 90 , en considérant le rapport a h 10 et hc h f 10
a Présent Rao et al. [310] Chalak et al. [303] Zhen et al. [296] Zhen et al. [296] Kant et Swaminathan [290]
b QSFT52 3D-Elasticity FEM-Q9-HZZT FEM-Q8-GLHSDT GLHSDT-Anal HSDT-Anal
161
Chapitre 5 Vibration libre des plaques composites stratifiées et sandwichs
Fig.5.2 Plaque sandwich carrée simplement appuyée avec des peaux stratifiées
0 / 90 / C / 0 / 90 pour différents rapports d’épaisseur - Comparaison des fréquences
naturelles
Fig.5.3 Plaque sandwich carrée simplement appuyée avec des peaux stratifiées
0 / 90 / C / 0 / 90 pour différents rapports (hc h f ) - Comparaison des fréquences naturelles
162
Chapitre 5 Vibration libre des plaques sandwichs et des composites stratifiées
Tableau 5.13 Propriétés mécaniques des peaux et de l’âme (HEREX - C70.130) [274]
Peaux en fibre de
24.51 7.77 3.34 3.34 1.34 0.078 1800
verre / polyester
Âme en mousse de
0.1036 0.1036 0.05 0.05 0.05 0.32 130
HEREX - C70.130
163
Chapitre 5 Vibration libre des plaques composites stratifiées et sandwichs
Tableau 5.14 Fréquences propres non-dimensionnelles ( ) d’une plaque sandwich carrée avec différentes séquences d'empilement dans les peaux
164
Chapitre 5 Vibration libre des plaques sandwichs et des composites stratifiées
Tableau 5.15 Fréquences propres non-dimensionnelles ( ) , pour différents modes, d’une plaque carrée composite stratifiée biaise 90 / 0 / 90 / 0 / 90
avec deux types de conditions aux limites (SSSS et CCCC)
Fréquences (Hz)
Skew angle Références Modèles EF
Mode 1 Mode 2 Mode 3 Mode 4 Mode 5 Mode 6 Mode 7
Simplement appuyée (SSSS)
Présent élément QSFT52 1.6020 3.0936 3.9511 4.7734 5.3915 6.5139 6.6478
Wang [313] Anal Solution 1.5699 3.0371 3.7324 4.5664 5.1469 6.0343 6.1820
0°
Garg et al. [293] FEM-Q9-TSDT 1.5703 2.8917 3.8041 4.5314 4.7273 5.8888 6.2125
Garg et al. [293] FEM-Q9-FSDT 1.5699 3.0372 3.7325 4.5664 5.1475 6.0347 6.1825
Présent élément QSFT52 1.6556 3.1025 4.0721 4.7873 5.4713 6.4469 6.8284
15° Garg et al. [293] FEM-Q9-TSDT 1.6877 3.0458 4.0264 4.4818 5.1533 5.8277 6.4910
Garg et al. [293] FEM-Q9-FSDT 1.6874 3.1413 3.9600 4.6073 5.4606 6.0880 6.3125
Présent élément QSFT52 1.8572 3.1609 4.4985 4.8378 5.7640 6.4343 7.4741
Wang [313] Anal Solution 2.0844 3.5127 4.6997 4.8855 6.2494 6.3118 7.2533
30°
Garg et al. [293] FEM-Q9-TSDT 2.0840 3.4023 4.7176 4.7674 5.9852 6.1177 7.1966
Garg et al. [293] FEM-Q9-FSDT 2.0884 3.5147 4.7033 4.8864 6.2514 6.3149 7.2570
Présent élément QSFT52 2.3630 3.3989 4.9769 5.4698 6.5139 6.6802 8.4955
Wang [313] Anal Solution 2.8825 4.2823 5.5868 6.1808 6.9022 7.8390 8.2114
45°
Garg et al. [293] FEM-Q9-TSDT 2.8925 4.1906 5.4149 6.2868 6.6302 7.7433 7.8398
Garg et al. [293] FEM-Q9-FSDT 2.8932 4.2852 5.5886 6.1874 6.9035 7.8366 8.1971
165
Chapitre 5 Vibration libre des plaques sandwichs et des composites stratifiées
Fréquences (Hz)
Skew angle Références Modèle EF
Mode 1 Mode 2 Mode 3 Mode 4 Mode 5 Mode 6 Mode 7
Encastrée (CCCC)
Présent élément QSFT52 2.4317 3.8296 4.3472 5.2745 5.9913 6.8568 7.0147
Wang [313] Anal Solution 2.3820 3.7383 4.1122 5.0366 5.6127 6.2393 6.5591
0°
Garg et al. [293] FEM-Q9-TSDT 2.3687 3.5399 4.2134 4.9852 5.1895 6.2801 6.4666
Garg et al. [293] FEM-Q9-FSDT 2.3820 3.7383 4.1122 5.0366 5.6133 6.2397 6.5595
Présent élément QSFT52 2.4610 3.8253 4.4573 5.2916 6.0549 6.9328 7.0249
15° Garg et al. [293] FEM-Q9-TSDT 2.4663 3.6255 4.4247 4.9395 5.5603 6.2086 7.2133
Garg et al. [293] FEM-Q9-FSDT 2.4750 3.7872 4.3418 5.0599 5.8762 6.4405 6.5178
Présent élément QSFT52 2.5812 3.8521 4.8330 5.3560 6.2860 6.9186 7.6454
Wang [313] Anal Solution 2.7921 4.0566 5.0220 5.2906 6.5698 6.6042 7.4485
30°
Garg et al. [293] FEM-Q9-TSDT 2.8001 3.9557 5.1079 5.1675 6.3624 6.4236 7.5067
Garg et al. [293] FEM-Q9-FSDT 2.7922 4.0568 5.0222 5.2909 6.5702 6.6047 7.4487
Présent élément QSFT52 2.9190 4.0216 5.5212 6.9025 7.1661 8.9485 8.9815
Wang [313] Anal Solution 3.4738 4.7393 5.9554 6.3750 7.1961 8.0031 8.4468
45°
Garg et al. [293] FEM-Q9-TSDT 3.5215 4.7129 5.8789 6.5290 7.0275 7.9963 8.1687
Garg et al. [293] FEM-Q9-FSDT 3.4739 4.7396 5.9553 6.3745 7.1920 7.9978 8.4213
166
Conclusion générale & perspectives
1. Conclusion générale:
Dans cette étude, un nouveau modèle a été développé pour l’analyse du comportement
des plaques composites multicouches et notamment les plaques sandwichs. Le modèle
proposé est basé sur l’approche layerwise, mais en gardant tous les avantages de l'approche
monocouche. Contrairement à tous les autres modèles layerwise présentés dans la littérature,
le nombre de degrés de liberté dans ce modèle est indépendant du nombre de couches. Les
deux peaux et l’âme sont modélisées individuellement en utilisant, respectivement, la théorie
de déformation en cisaillement du premier ordre et celle du troisième ordre. Ces deux théories
combinées, satisfont les exigences de continuités des déplacements interlaminaires ainsi que
les conditions de compatibilité. Le nouveau modèle permet de calculer les contraintes de
cisaillement transversal, d’une manière précise en utilisant les équations constitutives plutôt
que les équations d'équilibre et sans recours à des facteurs de correction.
Sur la base de ce nouveau modèle, deux éléments finis, baptisés RSFT52 et QSFT52
(Rectangular & Quadrilateral Sandwich First Third with 52-degrees of freedom), ont été
développés pour l’analyse statique et dynamique des plaques composites multicouches et des
sandwichs. Ces éléments sont géométriquement simples et ne possèdent que quatre nœuds et
treize degrés de liberté par nœud. Chaque nœud contient: deux degrés de liberté de rotation
pour chaque peau, six degrés de liberté de rotation aussi pour l’âme, alors que les trois degrés
de liberté de type translation (3ddl/nœud) sont communs aux trois couches. La formulation de
ces deux éléments est basée sur le principe du travail virtuel.
cas de chargement, de géométrie et des conditions aux limites. Les résultats obtenus ont été
comparés avec des solutions analytiques de références, et ceux obtenus par des modèles
d’éléments finis basés sur différentes théories (HSDT, HZZT, LW, GLHSDT). Pour l’analyse
statique, les exemples numériques montrent que les éléments finis développés sont capables
de donner d'excellents résultats en particulier pour les contraintes de cisaillement transverse.
Ces dernières ont été calculées directement en utilisant les équations constitutives plutôt que
les équations d’équilibre. Ainsi, dans le cas de l’analyse dynamique, les résultats de l’élément
développé QSFT52 sont en excellent accord avec ceux trouvés dans la littérature, où l’erreur
n’atteint pas 0.004 % par comparaison avec la solution d’élasticité-3D.
On peut conclure que l'utilisation des éléments finis développé ainsi que la combinaison
de la théorie de déformation en cisaillement du premier et celle du troisième ordre (FSDT,
TSDT), pour modéliser respectivement, les deux peaux et l’âme, ont montré une bonne
précision et une convergence rapide pour l’analyse les différents problèmes. Finalement, le
potentiel des éléments finis développés peut être utilisé pour générer un certain nombre de
nouveaux résultats, qui sont censés être utiles pour les futures recherches dans ce domaine.
2. Perspectives:
Les modèles d’éléments finis développés RSFT52 et QSFT52 peuvent être considérés
comme les plus simples parmi les éléments finis existants, quelle que soit l’approche ou le
modèle utilisé. Les résultats obtenus sont fort encourageants et permettent d’augurer
favorablement sur les possibilités réelles d’élargir ce modèle à l’étude d’autres phénomènes.
168
Conclusion générale & perspectives
Peut-on généraliser ce modèle aux problèmes des coques minces et épaisses? Cette
question mérite également une investigation.
169
Références bibliographiques
Références bibliographiques
[1] Reddy, J.N., Mechanics of laminated composite plates: theory and analysis. Vol. 1.
1997: CRC press Boca Raton.
[2] Reddy, J.N., Mechanics of laminated composite plates and shells: theory and analysis.
2003: CRC press.
[3] Kollár, L.P. and G.S. Springer, Mechanics of composite structures. 2003: Cambridge
university press.
[4] Vinson, J.R., The Behavior of Sandwich Structures of Isotropic and Composite
Materials: Jack R. Vinson, 1999: CRC Press.
[5] Zenkert, D., An introduction to sandwich construction1997: Engineering Materials
Advisory Services.
[6] Plantema, F.J., Sandwich construction. Vol. 966. 1966: Wiley New York.
[7] Jones, R.M., Mechanics of composite materials. Second Editions ed1998: CRC Press.
[8] Carlsson, L.A. and G.A. Kardomateas, Structural and failure mechanics of sandwich
composites. Vol. 121. 2011: Springer.
[9] Kreja, I.,(2011), A literature review on computational models for laminated composite
and sandwich panels, Central European Journal of Engineering, 1 (1) 59-80.
[10] Libove, C. and S. Batdorf, A general small-deflection theory for flat sandwich plates,
1948, DTIC Document.
[11] Reissner, E.,(1948), Finite Deflections of Sandwich Plates, Journal of the Aeronautical
Sciences (Institute of the Aeronautical Sciences), 15 (7) 435-440.
[12] Reissner, E., Small bending and stretching of sandwich-type shells, 1950, National
Advisory Committee for Aeronautics. Rpt 975.
[13] Hoff, N.J., Bending and buckling of rectangular sandwich plates.1950: National
Advisory Committee for Aeronautics.
[14] Eringen, A.C. Bending and buckling of rectangular sandwich plates. in Proceedings of
the First US National Congress of Applied Mechanics. 1951.
170
Références bibliographiques
[15] Yen, K.T., S. Gunturkun, and F.V. Pohle, Deflections of a simply supported
rectangular sandwich plate subjected to transverse loads1951: National Advisory
Committee for Aeronautics.
[16] Gerard, G.,(1951), Note on Bending of Thick Sandwich Plates, Journal of the
Aeronautical Sciences (Institute of the Aeronautical Sciences), 18 (6) 424-426.
[17] Ericksen, W.S.,( 1953, US Forest Products Lab, Madison WI, FPL-1828), The
bending of a circular sandwich plate under normal load.
[18] Ericksen, W.S. and H. March,(1958, rapport FPL n°1583-B. ), Effects of shear
deformation in the core of a flat rectangular sandwich panel: compressive buckling of
sandwich panels having dissimilar facings of unequal thickness.
[19] Cheng, S., On the theory of bending of sandwich plates, 1961, DTIC Document:
Proceedings of the 4 U.S. National Congress of A. M. p. 511–518.
[20] Lin, T.H. and L.T. Yokota,(1965), Deflections and bending moments of rectangular
sandwich panels with clamped edges under combined biaxial compressions and
pressure, AIAA Journal, 3 (6) 1162-1164.
[21] Liaw, B.D. and R.W. Little,(1967), Theory of bending multi-layer sandwich plates,
AIAA Journal, 5 (2) 301-304.
[22] Azar, J.J.,(1968), Bending theory for multilayer orthotropic sandwich plates, AIAA
Journal, 6 (11) 2166-2169.
[23] Pagano, N.,(1969), Exact solutions for composite laminates in cylindrical bending,
Journal of composite materials, 3 (3) 398-411.
[24] Pagano, N.,(1970), Exact solutions for rectangular bidirectional composites and
sandwich plates, Journal of composite materials, 4 (1) 20-34.
[25] Jones, A.T.,(1970), Exact natural frequencies for cross-ply laminates, Journal of
Composite materials, 4 (4) 476-491.
[26] Srinivas, S. and A. Rao,(1971), A three-dimensional solution for plates and laminates,
Journal of the Franklin Institute, 291 (6) 469-481.
[27] Srinivas, S., C. Joga Rao, and A. Rao,(1970), An exact analysis for vibration of
simply-supported homogeneous and laminated thick rectangular plates, Journal of
Sound and Vibration, 12 (2) 187-199.
[28] Noor, A., J. Peters, and W. Burton,(1994), Three‐Dimensional Solutions for Initially
Stressed Structural Sandwiches, Journal of Engineering Mechanics, 120 (2) 284-303.
[29] Noor, A.K. and W.S. Burton,(1990), Three-dimensional solutions for
antisymmetrically solutions for antisymmetrically laminated anisotropic plates,
Journal of applied mechanics, 57 (1) 182-188.
171
Références bibliographiques
[30] Kant, T. and K. Swaminathan,(2002), Analytical solutions for the static analysis of
laminated composite and sandwich plates based on a higher order refined theory,
Composite Structures, 56 (4) 329-344.
[31] Lascoup, B., et al.,(2010), Homogenization of the core layer in stitched sandwich
structures, Composites Science and Technology, 70 (2) 350-355.
[32] Lewiński, T.,(1993), On recent developments in the homogenization theory of elastic
plates and their application to optimal design: Part I, Structural optimization, 6 (1) 59-
64.
[33] Rabczuk, T., et al.,(2004), Homogenization of sandwich structures, International
Journal for Numerical Methods in Engineering, 61 (7) 1009-1027.
[34] Takano, N., M. Zako, and N. Kikuchi,(1995), Stress analysis of sandwich plate by the
homogenization method, Materials science research international, 1 (2) 82-88.
[35] Reddy, J.N.,(1993), An evaluation of equivalent-single-layer and layerwise theories of
composite laminates, Composite structures, 25 (1–4) 21-35.
[36] Carrera, E.,(2003), Theories and Finite Elements for Multilayered Plates and Shells:A
Unified compact formulation with numerical assessment and benchmarking, Archives
of Computational Methods in Engineering, 10 (3) 215-296.
[37] Kirchhoff, G.,(1850), Über das Gleichgewicht und die Bewegung einer elastischen
Scheibe, Journal für die reine und angewandte Mathematik, 40 51-88.
[38] Love, A.E.H., A treatise on the mathematical theory of elasticity1934: Cambridge
University Press.
[39] Reissner, E. and Y. Stavsky,(1961), Bending and Stretching of Certain Types of
Heterogeneous Aeolotropic Elastic Plates, Journal of Applied Mechanics, 28 (3) 402-
408.
[40] Librescu, L., Elastostatics and kinetics of anisotropic and heterogeneous shell-type
structures. Vol. 2. 1975, Noordhoff, Leyden, Netherlands
[41] Reddy, J.N.,(1989), On refined computational models of composite laminates,
International Journal for Numerical Methods in Engineering, 27 (2) 361-382.
[42] Stavsky, Y.,(1965), On the theory of symmetrically heterogeneous plates having the
same thickness variation of the elastic moduli, Topics in applied mechanics, 105.
[43] Rohwer, K., S. Friedrichs, and C. Wehmeyer,(2005), Analyzing laminated structures
from fibre-reinforced composite material–an assessment, Technische Mechanik, 25 (1)
59-79.
[44] Ambarsumyan, S.,(1958), K tieorii izgiba anisotropnvch plasinok, Izv. Akad Nauk
SSR, 5, 69-77.
[45] Dong S.B., P.K.S.e.T.R.L.,(1962), On the Theory of Laminated Anisotropic Shells and
Plates, Journal of the Aerospace Sciences, 29 (8) 969-975.
172
Références bibliographiques
173
Références bibliographiques
[63] Chatterjee, S.N. and S.V. Kulkarni,(1979), Shear Correction Factors for Laminated
Plates, AIAA Journal, 17 (5) 498-499.
[64] Madabhusi-Raman, P. and J.F. Davalos,(1996), Static shear correction factor for
laminated rectangular beams, Composites Part B: Engineering, 27 (3–4) 285-293.
[65] Laitinen, M., H. Lahtinen, and S.G. Sjölind,(1995), Transverse shear correction factors
for laminates in cylindrical bending, Communications in numerical methods in
engineering, 11 (1) 41-47.
[66] Carrera, E. and S. Brischetto,(2008), A Survey With Numerical Assessment of
Classical and Refined Theories for the Analysis of Sandwich Plates, Applied
Mechanics Reviews, 62 (1) 010803-010803.
[67] Reissner, E.,(1972), A consistent treatment of transverse shear deformations in
laminated anisotropic plates, AIAA Journal, 10 (5) 716-718.
[68] Belarbi, M.T., Développement de nouveaux éléments finis basés sur le modèle en
déformation. Application linéaire et non linéaire, 2000, Thèse de Doctorat d’état,
Université de Constantine.
[69] Ayad, R., Eléments finis de plaque et coque en formulation mixte avec projection en
cisaillement. Thèse de doctorat, Université de Technologie de Compiègne, France.
1993.
[70] Lardeur, P., Développement et évaluation de deux nouveaux éléments finis de plaques
et coques composites avec influence du cisaillement transversal, 1990, Thèse de
Doctorat, University of Technology of Compiègne,France.
[71] Noor, A.K. and W.S. Burton,(1989), Assessment of Shear Deformation Theories for
Multilayered Composite Plates, Applied Mechanics Reviews, 42 (1) 1-13.
[72] Bert, C.W.,(1973), Simplified analysis of static shear factors for beams of
nonhomogeneous cross section, Journal of Composite materials, 7 (4) 525-529.
[73] Yu, Y.-Y.,(1959), Simple Thickness-shear Modes of Vibration in Infinite Sandwich
Plates, (26) 679–681.
[74] Krishna Murty, A. and S. Vellaichamy,(1987), On higher order shear deformation
theory of laminated composite panels, Composite Structures, 8 (4) 247-270.
[75] Lo, K., R. Christensen, and E. Wu,(1977), A high-order theory of plate deformation—
part 2: laminated plates, Journal of Applied Mechanics, 44 (4) 669-676.
[76] Lo, K., R. Christensen, and E. Wu,(1977), A high-order theory of plate deformation—
Part 1: Homogeneous plates, Journal of Applied Mechanics, 44 (4) 663-668.
[77] Nelson, R.B. and D.R. Lorch,(1974), A Refined Theory for Laminated Orthotropic
Plates, Journal of Applied Mechanics, 41 (1) 177-183.
[78] Reddy, J.N.,(1984), A simple higher-order theory for laminated composite plates,
Journal of applied mechanics, 51 (4) 745-752.
174
Références bibliographiques
[79] Whitney, J. and C. Sun,(1973), A higher order theory for extensional motion of
laminated composites, Journal of Sound and Vibration, 30 (1) 85-97.
[80] Bhimaraddi, A. and L.K. Stevens,(1984), A Higher Order Theory for Free Vibration of
Orthotropic, Homogeneous, and Laminated Rectangular Plates, Journal of Applied
Mechanics, 51 (1) 195-198.
[81] Murthy, M.,(1981), An improved transverse shear deformation theory for laminated
antisotropic plates, NASA Technical, 1-37.
[82] Kant, T.,(1982), Numerical analysis of thick plates, Computer methods in applied
mechanics and engineering, 31 (1) 1-18.
[83] Ren-Huai, L. and H. Ling-Hui,(1991), A simple theory for non-linear bending of
laminated composite rectangular plates including higher-order effects, International
journal of non-linear mechanics, 26 (5) 537-545.
[84] Touratier, M.,(1991), An efficient standard plate theory, International Journal of
Engineering Science, 29 (8) 901-916.
[85] Hildebrand, F., E. Reissner, and G. Thomas, Notes on the foundations of the theory of
small displacements of orthotropic shells1949: NACA Technical Note No. 1833,
Washington, DC.
[86] Basset, A.B.,(1890), On the extension and flexure of cylindrical and spherical thin
elastic shells, Philosophical Transactions of the Royal Society of London. A, 181 (6)
433-480.
[87] Nguyen-Van, H., Development and application of assumed strain smoothing finite
element technique for composite plate/shell structures, Thèse de Doctorat, University
of Southern Queensland, 2009, 295 pp.
[88] Whitney, J. and C.-T. Sun,(1974), A refined theory for laminated anisotropic,
cylindrical shells, Journal of Applied Mechanics, 41 (2) 471-476.
[89] Kwon, Y. and J. Akin,(1987), Analysis of layered composite plates using a high-order
deformation theory, Computers & structures, 27 (5) 619-623.
[90] Fares, M.,(1993), A refined two-dimensional theory of anisotropic plates, Mechanics
Research Communications, 20 (4) 319-327.
[91] Tessler, A.,(1993), An improved plate theory of {1, 2}-order for thick composite
laminates, International Journal of Solids and Structures, 30 (7) 981-1000.
[92] Sacco, E. and J. Reddy,(1992), On first‐and second‐order moderate rotation theories of
laminated plates, International Journal for Numerical Methods in Engineering, 33 (1)
1-17.
[93] Reissner, E.,(1963), On the derivation of boundary conditions for plate theory,
Proceedings of the Royal Society of London. Series A. Mathematical and Physical
Sciences, 276 (1365) 178-186.
175
Références bibliographiques
[94] Levinson, M.,(1980), An accurate, simple theory of the statics and dynamics of elastic
plates, Mechanics Research Communications, 7 (6) 343-350.
[95] Reddy, J. and D. Robbins,(1994), Theories and computational models for composite
laminates, Applied mechanics reviews, 47, 147.
[96] Bert, C.W.,(1984), A critical evaluation of new plate theories applied to laminated
composites, Composite Structures, 2 (4) 329-347.
[97] Senthilnathan, N., et al.,(1987), Buckling of shear-deformable plates, AIAA Journal,
25 (9) 1268-1271.
[98] Lee, K., et al.,(1989), A simple higher-order non-linear shear deformation plate theory,
International journal of non-linear mechanics, 24 (2) 127-137.
[99] Rohwer, K.,(1992), Application of higher order theories to the bending analysis of
layered composite plates, International Journal of Solids and Structures, 29 (1) 105-
119.
[100] Ferreira, A., C. Roque, and R. Jorge,(2006), Static and free vibration analysis of
composite shells by radial basis functions, Engineering Analysis with Boundary
Elements, 30 (9) 719-733.
[101] Krishna Murty, A.,(1987), Flexure of composite plates, Composite Structures, 7 (3)
161-177.
[102] Aydogdu, M.,(2009), A new shear deformation theory for laminated composite plates,
Composite structures, 89 (1) 94-101.
[103] Levy, M.,(1877), Mémoire sur la théorie des plaques élastiques planes, Journal de
mathématiques pures et appliquées, 219-306.
[104] Mantari, J., A. Oktem, and C. Guedes Soares,(2012), A new trigonometric layerwise
shear deformation theory for the finite element analysis of laminated composite and
sandwich plates, Computers & Structures, 94 45-53.
[105] Mantari, J.L., A.S. Oktem, and C. Guedes Soares,(2011), Static and dynamic analysis
of laminated composite and sandwich plates and shells by using a new higher-order
shear deformation theory, Composite structures, 94 (1) 37-49.
[106] Mantari, J., A. Oktem, and C. Guedes Soares,(2012), A new higher order shear
deformation theory for sandwich and composite laminated plates, Composites Part B:
Engineering, 43 (3) 1489-1499.
[107] Soldatos, K.P.,(1992), A general laminated plate theory accounting for continuity of
displacements and transverse shear stresses at material interfaces, Composite
Structures, 20 (4) 195-211.
[108] Stein, M.,(1986), Nonlinear theory for plates and shells including the effects of
transverse shearing, AIAA Journal, 24 (9) 1537-1544.
176
Références bibliographiques
[109] Sahoo, R. and B.N. Singh,(2013), A new inverse hyperbolic zigzag theory for the
static analysis of laminated composite and sandwich plates, Composite structures, 105
(0) 385-397.
[110] Karama, M., K. Afaq, and S. Mistou,(2003), Mechanical behaviour of laminated
composite beam by the new multi-layered laminated composite structures model with
transverse shear stress continuity, International Journal of Solids and Structures, 40 (6)
1525-1546.
[111] Ghugal, Y.M. and A.S. Sayyad,(2010), A static flexure of thick isotropic plates using
trigonometric shear deformation theory, J Solid Mech, 2 (1) 79-90.
[112] Stein, M. and D.C. Jegley,(1987), Effects of transverse shearing on cylindrical
bending, vibration, and buckling of laminated plates, AIAA Journal, 25 (1) 123-129.
[113] Nguyen, V.T., Modélisation globale et locale des structures multicouches par éléments
finis de plaque. Thèse de Doctorat, Ecole des Ponts ParisTech. , 2004, 216 pp.
[114] Idlbi, A., M. Karama, and M. Touratier,(1997), Comparison of various laminated plate
theories, Composite Structures, 37 (2) 173-184.
[115] Soldatos, K.P.,(1992), A transverse shear deformation theory for homogeneous
monoclinic plates, Acta Mechanica, 94 (3-4) 195-220.
[116] Soldatos, K. and T. Timarci,(1993), A unified formulation of laminated composite,
shear deformable, five-degrees-of-freedom cylindrical shell theories, Composite
Structures, 25 (1) 165-171.
[117] Aydogdu, M.,(2006), Comparison of various shear deformation theories for bending,
buckling, and vibration of rectangular symmetric cross-ply plate with simply
supported edges, Journal of Composite materials, 40 (23) 2143-2155.
[118] Karama, M., K. Afaq, and S. Mistou,(2009), A new theory for laminated composite
plates, Proceedings of the Institution of Mechanical Engineers, Part L: Journal of
Materials Design and Applications, 223 (2) 53-62.
[119] Nguyen, D.T., Benchmark d'un modèle layer wise de multicouches et implémentation
du modèle dans Abaqus. Thèse de Doctorat ,Université Paris-Est., 2012, 160 pp.
[120] Kaczkowski, Z.,(1968), Plates, Statical calculations, Arkady, Warsaw.
[121] Panc, V., Theories of elastic plates. Prague: Academia, 1975: Springer.
[122] Reissner, E.,(1975), On transverse bending of plates, including the effect of transverse
shear deformation, International Journal of Solids and Structures, 11 (5) 569-573.
[123] Ferreira, A., C. Roque, and R. Jorge,(2005), Analysis of composite plates by
trigonometric shear deformation theory and multiquadrics, Computers & structures, 83
(27) 2225-2237.
177
Références bibliographiques
[124] Mantari, J., A. Oktem, and C. Guedes Soares,(2011), Static and dynamic analysis of
laminated composite and sandwich plates and shells by using a new higher-order shear
deformation theory, Composite Structures, 94 (1) 37-49.
[125] Mantari, J., A. Oktem, and C. Guedes Soares,(2012), A new trigonometric shear
deformation theory for isotropic, laminated composite and sandwich plates,
International Journal of Solids and Structures, 49 (1) 43-53.
[126] Mantari, J. and C. Guedes Soares,(2012), Analysis of isotropic and multilayered plates
and shells by using a generalized higher-order shear deformation theory, Composite
Structures, 94 (8) 2640-2656.
[127] Grover, N., D. Maiti, and B. Singh,(2013), A new inverse hyperbolic shear
deformation theory for static and buckling analysis of laminated composite and
sandwich plates, Composite Structures, 95 667-675.
[128] Sahoo, R. and B. Singh,(2013), A new shear deformation theory for the static analysis
of laminated composite and sandwich plates, International Journal of Mechanical
Sciences, 75 324-336.
[129] Carrera, E.,(2002), Theories and finite elements for multilayered, anisotropic,
composite plates and shells, Archives of Computational Methods in Engineering, 9 (2)
87-140.
[130] Saeedi, N., Une approche très efficace pour l'analyse du délaminage des plaques
stratifiées infiniment longues. Thèse de Doctorat ,Université Paris-Est., 2012, 223 pp.
[131] Kao, J.S. and R.J. Ross,(1968), Bending of multilayer sandwich beams, AIAA Journal,
6 (8) 1583-1585.
[132] Srinivas, S.,(1973), A refined analysis of composite laminates, Journal of Sound and
Vibration, 30 (4) 495-507.
[133] Sun, C.-T. and J. Whitney,(1973), Theories for the dynamic response of laminated
plates, AIAA Journal, 11 (2) 178-183.
[134] Swift, G.W. and R.A. Heller,(1974), Layered beam analysis, Journal of the
Engineering Mechanics Division, 100 (2) 267-282.
[135] Seide, P.,(1980), An improved approximate theory for the bending of laminated plates.
Chapter book.
[136] Reddy, J.N.,(1987), A generalization of two-dimensional theories of laminated
composite plates, Communications in Applied Numerical Methods, 3 (3) 173-180.
[137] Reddy, J.N., E.J. Barbero, and J.L. Teply,(1989), A plate bending element based on a
generalized laminate plate theory, International Journal for Numerical Methods in
Engineering, 28 (10) 2275-2292.
178
Références bibliographiques
[138] Barbero, E.J., J.N. Reddy, and J. Teply,(1990), An accurate determination of stresses
in thick laminates using a generalized plate theory, International Journal for Numerical
Methods in Engineering, 29 (1) 1-14.
[139] Robbins, D. and J. Reddy,(1993), Modelling of thick composites using a layerwise
laminate theory, International journal for numerical methods in engineering, 36 (4)
655-677.
[140] Cho, K., A. Striz, and C. Bert,(1989), Thermal stress analysis of laminate using
higher-order theory in each layer, Journal of thermal stresses, 12 (3) 321-332.
[141] Cho, K., C. Bert, and A. Striz,(1991), Free vibrations of laminated rectangular plates
analyzed by higher order individual-layer theory, Journal of Sound and Vibration, 145
(3) 429-442.
[142] Pagano, N.J. and S.R. Soni,(1983), Global-local laminate variational model,
International Journal of Solids and Structures, 19 (3) 207-228.
[143] Ren, J.G.,(1986), A new theory of laminated plate, Composites Science and
Technology, 26 (3) 225-239.
[144] Kassapoglou, C. and P.A. Lagace,(1986), An Efficient Method for the Calculation of
Interlaminar Stresses in Composite Materials, Journal of Applied Mechanics, 53 (4)
744-750.
[145] Kassapoglou, C. and P.A. Lagace,(1987), Closed form solutions for the interlaminar
stress field in angle-ply and cross-ply laminates, Journal of Composite materials, 21
(4) 292-308.
[146] Wan-Lee, Y.,(1994), Interlaminar stress analysis of composite laminates using a
sublaminate/layer model, International Journal of Solids and Structures, 31 (11) 1549-
1564.
[147] Yin, W.L.,(1994), Free-Edge Effects in Anisotropic Laminates Under Extension,
Bending and Twisting, Part I: A Stress-Function-Based Variational Approach, Journal
of Applied Mechanics, 61 (2) 410-415.
[148] Lee, K., et al.,(1990), An improved zig-zag model for the bending of laminated
composite plates, Composite Structures, 15 (2) 137-148.
[149] Di Sciuva, M.,(1984), A refined transverse shear deformation theory for multilayered
anisotropic plates, Atti. Accad. Sci. Torino, 118 279–295.
[150] Di Sciuva, M.,(1986), Bending, vibration and buckling of simply supported thick
multilayered orthotropic plates: an evaluation of a new displacement model, Journal of
Sound and Vibration, 105 (3) 425-442.
[151] Murakami, H.,(1986), Laminated composite plate theory with improved in-plane
responses, Journal of applied mechanics, 53 (3) 661-666.
179
Références bibliographiques
[152] Cho, M. and R.R. Parmerter,(1992), An efficient higher-order plate theory for
laminated composites, Composite Structures, 20 (2) 113-123.
[153] Cho, M. and R. Parmerter,(1993), Efficient higher order composite plate theory for
general lamination configurations, AIAA journal, 31 (7) 1299-1306.
[154] Lu, X. and D. Liu,(1992), An Interlaminar Shear Stress Continuity Theory for Both
Thin and Thick Composite Laminates, Journal of Applied Mechanics, 59 (3) 502-509.
[155] Di Sciuva, M.,(1987), An Improved Shear-Deformation Theory for Moderately Thick
Multilayered Anisotropic Shells and Plates, Journal of Applied Mechanics, 54 (3)
589-596.
[156] Averill, R.C.,(1994), Static and dynamic response of moderately thick laminated
beams with damage, Composites Engineering, 4 (4) 381-395.
[157] Cho, Y. and R. Averill,(2000), First-order zig-zag sublaminate plate theory and finite
element model for laminated composite and sandwich panels, Composite Structures,
50 (1) 1-15.
[158] Chakrabarti, A. and A.H. Sheikh,(2005), Analysis of laminated sandwich plates based
on interlaminar shear stress continuous plate theory, Journal of engineering mechanics,
131 (4) 377-384.
[159] He, L.-H.,(1994), A linear theory of laminated shells accounting for continuity of
displacements and transverse shear stresses at layer interfaces, International Journal of
Solids and Structures, 31 (5) 613-627.
[160] Khandelwal, R., A. Chakrabarti, and P. Bhargava,(2013), An efficient FE model based
on combined theory for the analysis of soft core sandwich plate, Computational
Mechanics, 51 (5) 673-697.
[161] Pandit, M., A.H. Sheikh, and B.N. Singh,(2010), Analysis of laminated sandwich
plates based on an improved higher order zigzag theory, Journal of sandwich
structures and materials, 12 (3) 307-326.
[162] Pandit, M.K., A.H. Sheikh, and B.N. Singh,(2008), An improved higher order zigzag
theory for the static analysis of laminated sandwich plate with soft core, Finite
elements in analysis and design, 44 (9) 602-610.
[163] Topdar, P., A.H. Sheikh, and N. Dhang,(2003), Finite element analysis of composite
and sandwich plates using a continuous inter-laminar shear stress model, Journal of
sandwich structures and materials, 5 (3) 207-231.
[164] Icardi, U.,(2001), Higher-order zig-zag model for analysis of thick composite beams
with inclusion of transverse normal stress and sublaminates approximations,
Composites Part B: Engineering, 32 (4) 343-354.
[165] Kapuria, S., P. Dumir, and N. Jain,(2004), Assessment of zigzag theory for static
loading, buckling, free and forced response of composite and sandwich beams,
Composite Structures, 64 (3) 317-327.
180
Références bibliographiques
181
Références bibliographiques
[181] Noor, A.K., W.S. Burton, and C.W. Bert,(1996), Computational models for sandwich
panels and shells, Applied mechanics reviews, 49 155.
[182] Zhang, Y. and C. Yang,(2009), Recent developments in finite element analysis for
laminated composite plates, Composite Structures, 88 (1) 147-157.
[183] Altenbach, H.,(1998), Theories for laminated and sandwich plates, Mechanics of
composite materials, 34 (3) 243-252.
[184] Ghugal, Y. and R. Shimpi,(2002), A review of refined shear deformation theories of
isotropic and anisotropic laminated plates, Journal of Reinforced Plastics and
Composites, 21 (9) 775-813.
[185] Wanji, C. and W. Zhen,(2008), A selective review on recent development of
displacement-based laminated plate theories, Recent patents on mechanical
engineering, 1 (1) 29-44.
[186] Chen, W. and P. Jia,(2012), Interlaminar stresses analysis and the limit state function
approximating methods for composite structure reliability assessment: A selected
review and some perspectives, Journal of Composite materials, .
[187] Li, X. and D. Liu,(1995), A laminate theory based on global–local superposition,
Communications in numerical methods in engineering, 11 (8) 633-641.
[188] Shariyat, M.,(2010), A generalized global–local high-order theory for bending and
vibration analyses of sandwich plates subjected to thermo-mechanical loads,
International Journal of Mechanical Sciences, 52 (3) 495-514.
[189] Wu, Z., R. Chen, and W. Chen,(2005), Refined laminated composite plate element
based on global–local higher-order shear deformation theory, Composite Structures,
70 (2) 135-152.
[190] Wu, Z., et al.,(2012), A higher order finite element including transverse normal strain
for linear elastic composite plates with general lamination configurations, Finite
Elements in Analysis and Design, 48 (1) 1346-1357.
[191] Zhen, W. and C. Wanji,(2010), A global-local higher order theory including
interlaminar stress continuity and C0 plate bending element for cross-ply laminated
composite plates, Computational Mechanics, 45 (5) 387-400.
[192] Zhen, W. and C. Wanji,(2010), A C0-type higher-order theory for bending analysis of
laminated composite and sandwich plates, Composite Structures, 92 (3) 653-661.
[193] Zhen, W. and C. Wanji,(2007), A study of global–local higher-order theories for
laminated composite plates, Composite Structures, 79 (1) 44-54.
[194] Chen, R.-g. and W.-j. Chen,(2000), Higher-order shear deformation theory and
triangular plate element based on global-local superposition, Acta Materiae
Compositae Sinica. (China), 17 96-102.
182
Références bibliographiques
183
Références bibliographiques
[208] Monforton, G. and L. Schmit. Finite element analysis of sandwich plates and
cylindrical shells with laminated faces. in Proceedings of the Second Conference of
Matrix Methods in Structural Mechanics, AFFSL-TR-68-150. 1968.
[209] Ahmed, K.,(1971), Static and dynamic analysis of sandwich structures by the method
of finite elements, Journal of Sound and Vibration, 18 (1) 75-91.
[210] Khatua, T. and Y. Cheung,(1973), Bending and vibration of multilayer sandwich
beams and plates, International journal for numerical methods in engineering, 6 (1)
11-24.
[211] Khatua, T.P. and Y.K. Cheung,(1972), Triangular element for multilayer sandwich
plates, Journal of the Engineering Mechanics Division, 98 (5) 1225-1238.
[212] Cook, R.D.,(1972), Two hybrid elements for analysis of thick, thin and sandwich
plates, International journal for numerical methods in engineering, 5 (2) 277-288.
[213] Cook, R.D.,(1975), Finite element buckling analysis of homogeneous and sandwich
plates, International Journal for Numerical Methods in Engineering, 9 (1) 39-50.
[214] Fazio, P.P. and K.H. Ha,(1974), Sandwich plate structure analysis by finite element,
Journal of the Structural Division, 100 (6) 1243-1262.
[215] Sheikh, A.H., S. Haldar, and D. Sengupta,(2002), A high precision shear deformable
element for the analysis of laminated composite plates of different shapes, Composite
Structures, 55 (3) 329-336.
[216] Tafla, A., Eléments finis mixtes-hybrides naturels sans facteurs correctifs du
cisaillement pour les plaques et les coques composites multicouches. Thèse de
Doctorat, University of Reims Champagne Ardenne, 2007, 132 pp.
[217] 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). Thèse de Doctorat, University of Reims Champagne Ardenne, 2008, 188 pp
[218] Katili, I., Formulation et évaluation de nouveaux éléments finis pour l'analyse linéaire
des plaques et coques de forme quelconque. Thèse de doctorat, Université de
Technologie de Compiègne, France, 1993.
[219] Kant, T., D. Owen, and O. Zienkiewicz,(1982), A refined higher-order C plate bending
element, Computers & structures, 15 (2) 177-183.
[220] Engblom, J.J. and O.O. Ochoa,(1985), Through‐the‐thickness stress predictions for
laminated plates of advanced composite materials, International journal for numerical
methods in engineering, 21 (10) 1759-1776.
[221] Pandya, B. and T. Kant,(1988), Flexural analysis of laminated composites using
refined higher-order C° plate bending elements, Computer Methods in Applied
Mechanics and Engineering, 66 (2) 173-198.
184
Références bibliographiques
185
Références bibliographiques
[236] Tu, T.M., L.N. Thach, and T.H. Quoc,(2010), Finite element modeling for bending
and vibration analysis of laminated and sandwich composite plates based on higher-
order theory, Computational Materials Science, 49 (4) S390-S394.
[237] Pandya, B. and T. Kant,(1988), Finite element analysis of laminated composite plates
using a higher-order displacement model, Composites Science and Technology, 32 (2)
137-155.
[238] Goswami, S. and W. Becker,(2013), A New Rectangular Finite Element Formulation
Based on Higher Order Displacement Theory for Thick and Thin Composite and
Sandwich Plates, World, 3 194-201.
[239] Lee, S.J. and H.R. Kim,(2013), FE analysis of laminated composite plates using a
higher order shear deformation theory with assumed strains, Latin American Journal
of Solids and Structures, 10 (3) 523-547.
[240] Spilker, R.,(1982), Hybrid‐stress eight‐node elements for thin and thick multilayer
laminated plates, International journal for numerical methods in engineering, 18 (6)
801-828.
[241] Fazio, P., K. Gowri, and K. Ha,(1987), Rectangular hybrid elements for the analysis of
sandwich plate structures, Canadian Journal of Civil Engineering, 14 (4) 455-460.
[242] Wu, C.-P. and C.-C. Lin,(1993), Analysis of sandwich plates using a mixed finite
element, Composite Structures, 25 (1) 397-405.
[243] Lee, L. and Y. Fan,(1996), Bending and vibration analysis of composite sandwich
plates, Computers & Structures, 60 (1) 103-112.
[244] Oskooei, S. and J. Hansen,(2000), Higher-order finite element for sandwich plates,
AIAA journal, 38 (3) 525-533.
[245] Ramtekkar, G., Y. Desai, and A. Shah,(2002), Mixed finite-element model for thick
composite laminated plates, Mechanics of Advanced Materials and Structures, 9 (2)
133-156.
[246] Ramtekkar, G., Y. Desai, and A. Shah,(2003), Application of a three-dimensional
mixed finite element model to the flexure of sandwich plate, Computers & Structures,
81 (22) 2183-2198.
[247] Demiray, S., W. Becker, and J. Hohe,(2004), A triangular v. Kármán type finite
element for sandwich plates with transversely compressible core, Computer methods
in applied mechanics and engineering, 193 (23) 2239-2260.
[248] Demiray, S., W. Becker, and J. Hohe, A Triangular Finite Element for Sandwich
Plates Accounting for Transverse Core Compressibility, Sandwich Structures 7:
Advancing with Sandwich Structures and Materials, 2005, Springer. p. 301-310.
[249] Linke, M., W. Wohlers, and H.-G. Reimerdes,(2007), Finite element for the static and
stability analysis of sandwich plates, Journal of sandwich structures and materials, 9
(2) 123-142.
186
Références bibliographiques
[250] Ramesh, S.S., Wang, CM, Reddy, JN, Ang, KK, (2009), A higher-order plate element
for accurate prediction of interlaminar stresses in laminated composite plates,
Composite Structures, 91 (3) 337-357.
[251] Ćetković, M. and D. Vuksanović,(2009), Bending, free vibrations and buckling of
laminated composite and sandwich plates using a layerwise displacement model,
Composite structures, 88 (2) 219-227.
[252] Darılmaz, K.,(2012), Analysis of sandwich plates: A three-dimensional assumed stress
hybrid finite element, Journal of sandwich structures and materials, 14 (4) 487-501.
[253] Darilmaz, K.,(2005), A hybrid 8-node hexahedral element for static and free vibration
analysis, Structural Engineering and Mechanics, 21 (5) 571-590.
[254] Di Sciuva, M.,(1993), A general quadrilateral multilayered plate element with
continuous interlaminar stresses, Computers & structures, 47 (1) 91-105.
[255] Di Sciuva, M.,(1995), A third‐order triangular multilayered plate finite element with
continuous interlaminar stresses, International Journal for Numerical Methods in
Engineering, 38 (1) 1-26.
[256] Carrera, E.,(1996), C0 Reissner-Mindlin multilayered plate elements including zig-zag
and interlaminar stress continuity, International Journal for Numerical Methods in
Engineering, 39 (11) 1797-1820.
[257] Carrera, E.,(1998), A refined multilayered finite-element model applied to linear and
non-linear analysis of sandwich plates, Composites Science and Technology, 58 (10)
1553-1569.
[258] Averill, R. and Y.C. Yip,(1996), Development of simple, robust finite elements based
on refined theories for thick laminated beams, Computers & structures, 59 (3) 529-
546.
[259] Cho, Y.-B. and R.C. Averill,(1997), An improved theory and finite-element model for
laminated composite and sandwich beams using first-order zig-zag sublaminate
approximations, Composite Structures, 37 (3) 281-298.
[260] Chakrabarti, A. and A.H. Sheikh,(2004), A new triangular element to model inter‐
laminar shear stress continuous plate theory, International journal for numerical
methods in engineering, 60 (7) 1237-1257.
[261] Singh, S.K., Chakrabarti, A., Bera, P. and Sony, J. S. D., (2011), An efficient C0 FE
model for the analysis of composites and sandwich laminates with general layup, Latin
American Journal of Solids and Structures, 8 (2) 197-212.
[262] Chalak, H.D., Chakrabarti, A., Iqbal, M. A. and Sheikh, H, A., (2012), An improved
C0 FE model for the analysis of laminated sandwich plate with soft core, Finite Elem.
Anal. Des., 56 20-31.
187
Références bibliographiques
[263] Rao, G.V., J. Venkataramana, and I. Raju,(1974), A high precision triangular plate
bending element for the analysis of thick plates, Nuclear Engineering and Design, 30
(3) 408-412.
[264] Yuan, F.G. and R.E. Miller,(1989), A cubic triangular finite element for flat plates
with shear, International journal for numerical methods in engineering, 28 (1) 109-
126.
[265] Kant, T., Gupta, Avani B., Pendhari, S. S., Desai, Yogesh M.,(2008), Elasticity
solution for cross-ply composite and sandwich laminates, Composite Structures, 83 (1)
13-24.
[266] Kulkarni, S. and S. Kapuria,(2007), A new discrete Kirchhoff quadrilateral element
based on the third-order theory for composite plates, Computational Mechanics, 39 (3)
237-246.
[267] Meftah, K., Modélisation numérique des solides par éléments finis volumiques basés
sur le concept SFR (Space Fiber Rotation). Thèse de Doctorat, Université Mohamed
Khider–Biskra, 2013.
[268] Belinha, J. and L.M.J.S. Dinis,(2006), Analysis of plates and laminates using the
element-free Galerkin method, Computers & structures, 84 (22–23) 1547-1559.
[269] Ghosh, A. and S. Dey,(1992), A simple finite element for the analysis of laminated
plates, Computers & Structures, 44 (3) 585-596.
[270] Xiao, J., et al.,(2008), Analysis of thick composite laminates using a higher-order
shear and normal deformable plate theory (HOSNDPT) and a meshless method,
Composites Part B: Engineering, 39 (2) 414-427.
[271] Pagano, N.J. and H.J. Hatfield,(1972), Elastic Behavior of Multilayered Bidirectional
Composites, AIAA Journal, 10 (7) 931-933.
[272] Folie, G.,(1970), Bending of clamped orthotropic sandwich plates, Journal of the
Engineering Mechanics Division, 96 (3) 243-265.
[273] Kanematsu, H.H., Y. Hirano, and H. Iyama,(1988), Bending and vibration of CFRP-
faced rectangular sandwich plates, Composite Structures, 10 (2) 145-163.
[274] Meunier, M. and R. Shenoi,(1999), Free vibration analysis of composite sandwich
plates, Proceedings of the Institution of Mechanical Engineers, Part C: Journal of
Mechanical Engineering Science, 213 (7) 715-727.
[275] Tati, A., Analyse du flambage des plaques stratifiees en matériaux composites par la
méthode des éléments finis. Thèse de Doctorat, Université Mohamed Khider de
Biskra, 2007, 144 pp.
[276] Ergatoudis, I., B.M. Irons, and O.C. Zienkiewicz,(1968), Curved, isoparametric,
“quadrilateral” elements for finite element analysis, International Journal of Solids and
Structures, 4 (1) 31-42.
188
Références bibliographiques
[277] Nayak, A., S. Moy, and R. Shenoi,(2002), Free vibration analysis of composite
sandwich plates based on Reddy's higher-order theory, Composites Part B:
Engineering, 33 (7) 505-519.
[278] Lee, S.,(2004), Free vibration analysis of plates by using a four-node finite element
formulated with assumed natural transverse shear strain, Journal of sound and
vibration, 278 (3) 657-684.
[279] Kabir, H.R.H.,(1995), A shear-locking free robust isoparametric three-node triangular
finite element for moderately-thick and thin arbitrarily laminated plates, Computers &
Structures, 57 (4) 589-597.
[280] Chalak, H.D., et al.,(2014), C0 FE model based on HOZT for the analysis of laminated
soft core skew sandwich plates: Bending and vibration, Applied Mathematical
Modelling, 38 (4) 1211-1223.
[281] Dey, P., A. Sheikh, and D. Sengupta,(2014), A new element for the analysis of
composite plates, Finite Elements in Analysis and Design, 82 62-71.
[282] Yu, Y.-Y.,(1960), Flexural vibrations of elastic sandwich plates, Journal of Aerospace
Sciences, (27) 272–282.
[283] Yu, Y.-Y.,(1960), Simplified Vibration Analysis of Elastic Sandwich Plates, Journal
of the Aerospace Sciences, 27 (12) 894-900.
[284] Srinivas, S. and A.K. Rao,(1970), Bending, vibration and buckling of simply
supported thick orthotropic rectangular plates and laminates, International Journal of
Solids and Structures, 6 (11) 1463-1481.
[285] Chan, H.C. and Y.K. Cheung,(1972), Static and dynamic analysis of multi-layered
sandwich plates, International Journal of Mechanical Sciences, 14 (6) 399-406.
[286] Noor, A.K.,(1973), Free vibrations of multilayered composite plates, AIAA Journal,
11 (7) 1038-1039.
[287] Ibrahim, I.M., M.N.F. Rizk, and A. Farah,(1981), Dynamic analysis of unbalanced
anisotropic sandwich plates, Journal of the Engineering Mechanics Division, 107 (2)
405-418.
[288] Reddy, J.N. and T. Kuppusamy,(1984), Natural vibrations of laminated anisotropic
plates, Journal of sound and vibration, 94 (1) 63-69.
[289] Wang, C., et al.,(2000), Free vibration of skew sandwich plates with laminated
facings, Journal of sound and vibration, 235 (2) 317-340.
[290] Kant, T. and K. Swaminathan,(2001), Analytical solutions for free vibration of
laminated composite and sandwich plates based on a higher-order refined theory,
Composite Structures, 53 (1) 73-85.
189
Références bibliographiques
[291] Jam, J.E., B. Eftari, and S.H. Taghavian,(2010), A new improved high-order theory for
analysis of free vibration of sandwich panels, Polymer Composites, 31 (12) 2042-
2048.
[292] Khare, R.K., A.K. Garg, and T. Kant,(2005), Free vibration of sandwich laminates
with two higher-order shear deformable facet shell element models, Journal of
Sandwich Structures and Materials, 7 (3) 221-244.
[293] Garg, A.K., R.K. Khare, and T. Kant,(2006), Free vibration of skew fiber-reinforced
composite and sandwich laminates using a shear deformable finite element model,
Journal of Sandwich Structures and Materials, 8 (1) 33-53.
[294] Nabarrete, A., S.F. M. De Almeida, and J.S. Hansen,(2003), Sandwich-plate vibration
analysis: three-layer quasi-three-dimensional finite element model, AIAA Journal, 41
(8) 1547-1555.
[295] Roque, C., A. Ferreira, and R. Jorge,(2006), Free vibration analysis of composite and
sandwich plates by a trigonometric layerwise deformation theory and radial basis
functions, Journal of Sandwich Structures and Materials, 8 (6) 497-515.
[296] Zhen, W., C. Wanji, and R. Xiaohui,(2010), An accurate higher-order theory and C0
finite element for free vibration analysis of laminated composite and sandwich plates,
Composite Structures, 92 (6) 1299-1307.
[297] Elmalich, D. and O. Rabinovitch,(2012), A high-order finite element for dynamic
analysis of soft-core sandwich plates, Journal of Sandwich Structures and Materials,
14 (5) 525-555.
[298] Chakrabarti, A. and A.H. Sheikh,(2004), Vibration of laminate-faced sandwich plate
by a new refined element, Journal of Aerospace Engineering, 17 (3) 123-134.
[299] Chakrabarti, A. and A. Sheikh. A new triangular element based on a layer wise zigzag
theory for the analysis of composite panels having sandwich construction. in
Proceedings of the International Conference on Ship and Ocean Technology. 2002.
[300] Kulkarni, S. and S. Kapuria,(2008), Free vibration analysis of composite and sandwich
plates using an improved discrete Kirchhoff quadrilateral element based on third-order
zigzag theory, Computational Mechanics, 42 (6) 803-824.
[301] Batoz, J.L. and M.B. Tahar,(1982), Evaluation of a new quadrilateral thin plate
bending element, International Journal for Numerical Methods in Engineering, 18 (11)
1655-1677.
[302] Jeyachandrabose, C., J. Kirkhope, and L. Meekisho,(1987), An improved discrete
Kirchhoff quadrilateral thin‐plate bending element, International Journal for
Numerical Methods in Engineering, 24 (3) 635-654.
[303] Chalak, H.D., et al.,(2013), Free vibration analysis of laminated soft core sandwich
plates, Journal of Vibration and Acoustics, 135 (1) 011013.
190
Références bibliographiques
[304] Sahoo, R. and B. Singh,(2014), A new trigonometric zigzag theory for buckling and
free vibration analysis of laminated composite and sandwich plates, Composite
Structures, 117 316-332.
[305] Ferreira, A.J., MATLAB codes for finite element analysis: solids and structures. Vol.
157. 2008: Springer.
[306] Marjanović, M. and D. Vuksanović,(2014), Layerwise solution of free vibrations and
buckling of laminated composite and sandwich plates with embedded delaminations,
Composite Structures, 108 9-20.
[307] Vuksanović, D.,(2000), Linear analysis of laminated composite plates using single
layer higher-order discrete models, Composite Structures, 48 (1–3) 205-211.
[308] Owen, D.R.J. and Z.H. Li,(1987), A refined analysis of laminated plates by finite
element displacement methods—II. Vibration and stability, Computers & Structures,
26 (6) 915-923.
[309] Khandelwal, R.P., A. Chakrabarti, and P. Bhargava,(2013), Vibration and buckling
analysis of laminated sandwich plate having soft core, International Journal of
Structural Stability and Dynamics, 13 (08).
[310] Rao, M., et al.,(2004), Natural Vibrations of Laminated and Sandwich Plates, Journal
of Engineering Mechanics, 130 (11) 1268-1278.
[311] Kant, T. and K. Swaminathan,(2001), free vibration of isotropic, orthotropic, and
multilayer plates based on higher order refined theories, Journal of sound and
vibration, 241 (2) 319-327.
[312] Rao, M.K. and Y.M. Desai,(2004), Analytical solutions for vibrations of laminated
and sandwich plates using mixed theory, Composite Structures, 63 (3–4) 361-373.
[313] Wang, S.,(1997), Free vibration analysis of skew fibre-reinforced composite laminates
based on first-order shear deformation plate theory, Computers & Structures, 63 (3)
525-538.
[314] Malekzadeh, K. and A. Sayyidmousavi,(2009), Free vibration analysis of sandwich
plates with a uniformly distributed attached mass, flexible core and different boundary
conditions, Journal of Sandwich Structures and Materials, 12(6) 709–732.
191