Effet de la température sur les fréquences propres des plaques

stratifiées en matériaux composites par la méthode des éléments finis

MASMOUDI Amina¹ , TATI Abdelouahab²
¹ Laboratoire de Génie Civil LRGC, Université Mohamed KHEIDER Biskra, BP 145, 07000 Biskra, Algérie
² Laboratoire de Génie Energétique et Matériaux, LGEM. Université Mohamed KHEIDER Biskra, BP 145, 07000 Biskra, Algérie

Abstract
4. Résultats numériques Effet de nombre des couches :
Le travail présenté a pour but d’étudier l’effet de la température sur les fréquences propre des plaques composites stratifiées en utilisant un élément fini rectangulaire à quatre nœuds avec
cinq degrés de liberté basé sur la théorie du premier ordre de cisaillement transverse . Les résultats obtenus, à travers une série d’exemples, sont comparés avec les solutions analytiques et
ceux obtenus via des modèles d’éléments finis disponibles dans la littérature. Cette comparaison a montré la bonne performance et la précision de l’élément utilisé.
A/Validation numérique d’élément fini utilisé : Cette section nous allons présenter l’effet du nombre des couches et la température sur
les fréquences des plaques stratifiées carrées (0/90/…./90/0). Les caractéristiques de la
Ce travail a été achevé par une étude paramétrique qui a pour but de mettre en évidence les effets de nombre des couches sur la variation de fréquences propres des plaques stratifiées plaque sont données Tableau 03
soumises à un chargement de température. Tableau 02 présente la convergence de la température critique de flambement d'une
plaque stratifiée simplement appuyée de 10 couches avec stratification [+𝜃°/−𝜃°/+𝜃°/ … ./
− 𝜃°] 𝜃=15°. Par comparaison avec les résultats obtenus par l’analyse 3D de NOOR on Interprétation des résultats
1. Introduction remarque que les valeurs convergent avec les premiers nombres d’éléments et restent
constantes ce qu’indique la stabilité de la plaque. Les fréquences propres diminuent
Récemment les plaques en matériaux composites sont largement utilisées en différents domaines de construction et l'industrie grâce aux avantages qu'elles présentent. Elles sont souvent avec l’augmentation de la
appliquées dans un environnement avec température élevée qui causée des changements dans les comportements mécaniques des plaques. Les résultats de cette étude montrent l’influence Tableau 1 : Les caractéristiques mécaniques de la plaque. température. Lorsque la
de nombre des couches, la précontrainte et les contraintes thermique sur le comportement vibratoire et flambement des plaques avec l’indépendance et la dépendance de propriété du
matériau à la température. E11/E22 G12/E22 G13/E22 G23/E22 𝝑𝟏𝟐 𝜶𝑳 /𝜶𝟎 𝜶𝑻 /𝜶𝟎 𝝆 a/h température donnée est égale à la
température critique, les
15 0.2 0.2 0.3356 0.3 0.015 1 1 100
2. Etude théorique : fréquences propres tendent vers
3. Formulation de l’élément fini
Tableau 02 : Température critique de plaque stratifiée simplement appuyée
zéro.
Les efforts d'une plaque et les moments sont
La cinématique adoptée est celle de la liés aux déformations et aux courbures par les L’élément choisi, dans ce travail, est un élément Figure 02 : Variation du carré de fréquence de plaque
rectangulaire à 4 nœuds avec cinq degré de encastrée en fonction de température T0/Tcr avec
théorie de premier ordre de cisaillement expressions suivantes: Présent 4×4 8×8 12 × 12 14 × 14 16 × 16 Noor[2] variation du nombre des couches
des stratifiés qui est l'extension de la liberté 𝒖, 𝒗, 𝒘, 𝝋𝒙 , 𝝋𝒚 par nœud (figure1). Les
élément Interprétation des résultats
théorie de Reissner-Mindlin. Cette longueurs 2a et 2b sont les dimensions de
l’élément dans les directions x et y dans cet Températu 0.0012 0.0012 0.0011 0.0011 0.0011 0.00115
théorie tient en compte des La variation de fréquence de plaque
déformations dues au cisaillement ordre. re critique
encastrée en fonction du nombre des
transverse [1]. La résultante de force thermique est définies Figure 01: géométrie et variables nodales de l’élément. couches avec température imposée.
La cinématique de cette théorie peut par:
On remarque que les valeurs de
s'écrire sous la forme suivante: B/Etude paramétrique fréquence augmentent en
Le champ de déplacements d’un point M(x,y) de la plaque peut être exprimé par
augmentant du nombre des couches
Après la vérification de performance de l’élément, nous allons passer à une étude
les résultantes en cisaillement est donnée par jusqu’au nombre de couches est
Où δαi vecteur de déplacement paramétrique pour voir l’influence de différents paramètres et la température
égale à 9.
L’état de déformation à n’importe quel critique sur la vibration des plaques composites stratifiées.
point dans la plaque est donné par les Nα (x, y) les fonctions de forme bilinéaire de Lagrange associé au nœud , Le maillage adopté pour cette étude est (16×16). Figure 03 : variation de fréquence de plaque
deux vecteurs : respectivement. encastrée en fonction nombre de couche.
Les conditions aux limites 𝑢 = 𝑣 = 𝑤 = 𝜑𝑥 = 𝜑𝑦 = 0 au X=0,a
L’équation constitutive peut également être
La relation contrainte-déformation dans écrite sous forme contractée suivant 𝑢 = 𝑣 = 𝑤 = 𝜑𝑥 = 𝜑𝑦 = 0 au Y=0,a 5. Conclusion
une couche k est donnée par :
cette étude nous a permis de montrer la bonne performance de l’élément utilisé.
Les valeurs propres de fréquence sont normalisées par la formule : Après la validation de l’élément on a passé pour étudier l’effet du nombre des
𝜌 𝜔
ഥ𝑎
2
𝑇0 couches, sur les fréquences naturelles des plaques stratifiées. Cette analyse nous
𝜔
ഥ = 𝜔𝑎 2
. =1− a permis de tirer
𝐸22 ℎ2 𝜔
ഥ0 𝑇𝑐𝑟
-Les fréquences propres de vibration diminuent avec l’augmentation de la
: Coefficients de dilatation thermique L’énergie potentielle totale de l’élément est donné par 𝜋 =𝑈−𝑊 température et s’annulent lorsque la température atteint température critique.
𝑈 est l’énergie potentielle de déformation et W est le travail des forces extérieures. L’énergie -En augmentant le nombre des couches les valeurs de la fréquence croient
potentielle totale peut être écrite en termes des efforts résultants et les déformations Tableau 03: Les caractéristiques mécaniques de la plaque
jusqu’à 9 couches pour rester pratiquement inchangées quelle que soit la
correspondantes comme suit : température
E11/E22 G12/E22 G13/E22 G23/E22 𝝑𝟏𝟐 𝜶𝑳 /𝜶𝟎 𝜶𝑻 /𝜶𝟎 𝝆 a/h
Vibrations libres d’une plaque soumise aux contraintes thermiques :
Où A est l’air de l’élément.
Le principe de Hamilton principale.
40 0.6 0.6 0.5 0.25 1,14 × 10−6 11,4 × 10−6 1 10 Référence
En utilisant les équations précédentes on aura :
En utilisant l’équation de Lagrange, l’équation of mouvement peut être exprimé par
[1] Reddy J.N., 2004, Mechanics of laminated composite plates and shells: theory and
𝑀 :matrice masse , 𝐾𝑇 : matrice de rigidité totale incluant l’effet thermique : L’annulation de la seconde variation de l’énergie potentielle totale par rapport à au vecteur des analysis, by CRC press.
déplacements 𝑞 , on obtient le problème de valeur propre suivant : [2] Noor, A.K. and W.S. Burton, Three-dimensional solutions for the free vibrations and
En substituent 𝑞ሷ par −𝜔2 𝑞 dans l’équation précédant on aura le problème de valeur propre
buckling of thermally stressed multilayered angle-ply composite plates. Journal of Applied
suivant : En utilisant le facteur de charge, 𝑁𝑡 = 𝑁𝑡0 Mechanics, 1992. 59(4): p. 868-877.
Et la température critique Δ𝑇cr = λcr . ΔT