These Karmi Yacine
These Karmi Yacine
These Karmi Yacine
Thèse
Présentée pour obtenir le diplôme de
Doctorat
Spécialité : Mécanique Avancée
Par :
KARMI Yacine
Thème :
Etude vibro-acoustique des structures avec matériaux
viscoélastiques sous chargement dynamique
Je tiens d’abord à rendre grâce à ALLAH, le tout puissant qu’il m’a accordé, la patience
et la volonté pour terminer ce travail.
Je voudrais aussi remercier Pr. Doauadji Ali de l’INSA de Lyon, de m’avoir accueilli
dans son laboratoire, pour son aide durant mon séjour de stage et pour ses conseils inoubliables.
Mes remerciements également vont aux membres du jury le Professeur Lagred Ahmed,
le Professeur Djamaa Mohamed Cherif, le Professeur Guedri Abdelmoumene, Dr. Dib Amar
et Dr. Nehal Abdelaziz pour l’intérêt qu’ils portent pour ce travail ainsi que pour le temps qu’ils
ont dépensé pour examiner ce présent manuscrit et aussi les remarques et les suggestions
enrichissantes qui m’ont aidé à améliorer le contenu de cette thèse.
Bien sûr, j’adresse mes reconnaissances particulières à mes très chers parents pour leur
soutien et aide sans faille qu’ils m’ont apporté durant tout le cycle de ma scolarité.
Enfin, je remercie tous ceux qui ont participé de près ou de loin, à la réalisation de ce travail.
i
Dédicaces
KARMI Yacine
ii
ملخص
في هذا العمل ،نحن مهتمون بدراسة السلوك االهتزازي الصوتي لبعض الهياكل املركبة على شكل عوارض و ألواح التي
ً
جديدا يمكن أن يحسن يتكون قلبها من مادة مرنة لزجة .يعتبر استخدام املواد اللزجة املطاطية املدمجة في الهياكل اقتر ً
احا
االستقرار الوظيفي للهياكل واملنشآت في مختلف القطاعات (الصناعة ،الهندسة املدنية ،الطيران ،الفضاء ، )...باإلضافة
إلى تقديم حل تنافس ي وموثوق للمشاكل التي تسببها االهتزازات الهيكلية والضوضاء الصوتية الناتجة عنها .ينتج عن هذا
اضطرابات أو خلل في مختلف املنشأت .يتم تنفيذ نمذجة الظواهر االهتزازية والصوتية للهياكل ذات التخميد السلبي بواسطة
طرق عددية متقدمة ،مع مراعاة :قوى القصور الذاتي الطولي والدوراني ؛ عدم التناسق الهندس ي خصائص املواد اإلجهادات
واللخطية لقوانين سلوك اللزوجة املرنة.
الحرارية ا
يعتمد النهج العددي املقدم على مبدأ القوى االفتراضية املرتبطة بمبدأ هاملتونا ) (Hamiltonمن أجل إنشاء املعادلة
العامة التي تحكم الحركة االهتزازية و وصف دينميكيتها .مما يوفر نظام معادلة خطية عالي الترتيب ،يتكيف مع أنواع قوى
التحميل املختلفة ،الظروف الحدودية والظروف البيئية ،أي التأثير الحراري على سلوك املنشئة .يتم إكمال هذا النهج من
ً ً
تنافسيا لنظام املعادالت الذي تم الحصول عليه من حيث وقت خلل التحليل الشكلي ( ، )Analyse modaleمما يتيح حل
الحساب دون فقدان الدقة ،مع تجنب معدل تنفيذ مرتفع للخوارزمية املطبقة في ماتلب ( .)Matlabمن أجل بناء هذا
القاعدة الشكلية ،يجب إجراء حل ملشكلة القيمة الذاتية املعقدة وغير الخطية .علوة على ذلك ،فإن نماذج اللزوجة املرنة
غير خطية ،ويفسر ذلك عدم خطية مصفوفة الصلبة .لهذا الغرض ،يتم تطبيق خوارزمية طريقة التقارب العددي
جنبا إلى جنب مع طريقة لتفاضل التلقائي ( .)Différentiation Automatiqueلذلك(ً )Méthode Asymptotique Numérique
،تم تكييف مخطط تكامل نيومارك ) (Newmarkلحل معادلة الحركة املخفضة بعد استخدام القاعدة الشكليةً .
أخيرا ،تم
حل مشكلة االهتزازات الصوتية من خلل تطبيق نظرية اقتران بنية السوائل التي تم صياغتها من خلل الجمع بين طريقة
العناصر املحدودة( )Méthode des Eléments Finisوطريقة العناصر الحدودية ( .)Méthode des Eléments de frontièresتم
فحص العديد من املعلمات الصوتية مثل عامل فقدان نقل الصوت ( )Perte de Transmission Acoustiqueعلى نطاق واسع
من أجل تقييم فعالية التحكم الخامد.
كلمات املفتاحية :التحكم الخامد ،التخميد ،املواد اللزجة املرنة ،الساندويتش ،االهتزازية .
iii
Abstract
This research focuses on the study of the vibratory and vibro-acoustic behavior of passive
damped beam and plate sandwich structures, in which the core is made of a viscoelastic material. To
use viscoelastic materials integrated into structures is considered a new suggestion that can improve
the functional stability of structures and systems in different sectors (industry, civil engineering,
aeronautic, aerospace ...), as well as, offering a competitive and reliable solution for the problems
caused by structural vibrations and the acoustic noises generated by them. This results in
disturbances or dysfunctions of the various installations. The modeling of the vibratory and vibro-
acoustic phenomena of passively damped structures is carried out by advanced numerical methods,
taking into account: longitudinal and rotational inertia forces; geometrical asymmetry; material
properties; thermal stresses and the non-linearity of the viscoelastic behavior laws.
The presented numerical approach is mainly based on the principle of virtual works associated
with Hamilton's principle in order to establish a general equation governing the vibratory motion
and describing its dynamics. The discretization of the motion equation is accomplished using the
finite element method (FEM), providing a high-order linear equation system, adapted to the various
loading types; boundary conditions and environmental conditions, namely, the thermal effect on the
structure behavior. This approach is completed by the modal analysis (MA), permitting a
competitive resolution of the obtained system of equations in terms of calculation time without
losing accuracy, while avoiding a high execution rate of the algorithm implemented in Matlab. In
order to build this modal basis, a resolution of the complex and nonlinear eigenvalue problem must
be performed. Moreover, the viscoelastic models are strongly nonlinear, this is explained by the
nonlinearity of the rigidity matrix. For this purpose, the algorithm of the numerical asymptotic
method (MAN) combined with the automatic differentiation method (MDA), is applied. Therefore,
the Newmark integration scheme has been adapted to solve the reduced equation of motion after
using the modal basis from the eigenvalue problem. Finally, the solution of the vibro-acoustic
problem has been approached by applying a fluid-structure coupling theory which formulation is
performed by combining the finite element method with the boundary element method (BEM).
Several acoustic parameters such as the acoustic transmission loss factor (STL) have been
extensively examined in order to evaluate the effectiveness of the passive damping treatment.
iv
Keywords: passive control, damping, viscoelastic materials, sandwich, vibro-acoustics.
Résumé
Dans le présent manuscrit, nous nous intéressons à l’étude du comportement vibroacoustique des
structures poutres et plaques sandwichs à amortissement passif dont le cœur est constitué d’un
matériau viscoélastique. L'utilisation des matériaux viscoélastiques intégrés aux structures est
considérée comme une nouvelle suggestion qui peut améliorer la stabilité fonctionnelle des
structures et installations dans différents secteurs (industrie, génie civil, aéronautique,
aérospatial…), offrant également une solution compétitive et fiable face aux problèmes provoquant
les vibrations structurelles et les bruits acoustiques engendrés. Ce qui en résulte des perturbations ou
de dysfonctionnements des diverses installations. La modélisation des phénomènes vibratoires et de
vibro-acoustiques des structures à amortissement passif est effectuée par des méthodes numériques
avancées, parfaitement adaptées, tenant en compte : les forces d’inerties longitudinales et
rotationnelles, l’asymétrie géométrique, les propriétés matérielles, les contraintes thermiques et la
non linéarité des lois du comportement viscoélastique.
L’approche numérique présentée est basée principalement sur le principe de travaux virtuels
associé au principe de Hamilton dans le but d’établir l’équation générale gouvernant le mouvement
vibratoire et décrivant sa dynamique. La discrétisation de l’équation du mouvement est accomplie
par la méthode des éléments finis (MEF), permettant d’avoir un système d’équations linéaires
d’ordre supérieur, adapté aux différents types de chargement ; des conditions aux limites ; ainsi
qu’aux conditions environnementales, à savoir, l’effet thermique influant sur le comportement des
structures étudiées. Cette approche est complétée par l’analyse modale (AM), permettant une
résolution compétitive du système d’équations obtenues en matière de temps de calcul sans pour
autant perdre de précision, tout en évitant un taux d’exécution élevée de l’algorithme implémenté
sous Matlab. Afin de construire cette base modale, une résolution du problème aux valeurs propres
complexe et non linéaire doit être effectuée. Par ailleurs, les modèles viscoélastiques sont fortement
non linéaire, ceci s’explique par la non linéarité de la matrice de rigidité. À cet effet, l’algorithme de
la méthode asymptotique numérique (MAN) combinée avec la méthode de différenciation
automatique (MDA), est appliqué. Ainsi, le schéma d’intégration de Newmark a été adapté pour
résoudre l’équation réduite du mouvement après avoir utilisé la base modale issue du problème aux
valeurs propres. Pour terminer, la résolution du problème vibro-acoustique a été abordée en
v
appliquant une théorie de couplage fluide-structure dont la formulation est effectuée en combinant
la méthode des éléments finis à la méthode des éléments de frontières (BEM). Plusieurs paramètres
acoustiques tels que le facteur de perte de transmission acoustique (STL) ont été largement examinés
afin d’évaluer l’efficacité du traitement passive d’amortissement.
vi
Table des matières
viii
Table des matières
ix
Table des matières
x
Table des figures
xi
Liste des tableaux
xii
Table des figures
Figure III.A.5. Variations des fréquences naturelles et des facteurs de perte de la poutre sandwich en
appuis simples pour différentes valeurs de 𝜃 avec 𝜂𝑣 = 0.1 et ℎ𝐻 = 0.1 ((a) fréquence naturelle ; (b)
facteur de perte). .................................................................................................................................. 97
Figure III.A.6. Comparaison des réponses fréquentielles au milieu de la poutre sandwich en appuis
simples de cœur viscoélastique et faces composites obtenues pour différentes valeurs de 𝜃 avec 𝜂𝑣 =
0.1 et ℎ𝐻 = 0.1.................................................................................................................................... 97
Figure III.A.7. Comparaison des réponses fréquentielles à l’extrémité de la poutre sandwich en porte-
à-faux de cœur viscoélastique et faces composites obtenues pour différentes valeurs de 𝜃 avec 𝜂𝑣 = 0.1
et ℎ𝐻 = 0.1. ......................................................................................................................................... 98
Figure III.A.8. Variations des fréquences naturelles et des facteurs de perte de la poutre sandwich en
appuis simples pour différentes valeurs de ℎ𝐻 avec 𝜂𝑣 = 0.1 et 𝜃 = 0° ((a) fréquence naturelle ; (b)
facteur de perte). .................................................................................................................................. 99
Figure III.A.9. Comparaison des réponses fréquentielles au milieu de la poutre sandwich en appuis
simples de cœur viscoélastique et faces composites obtenues pour différentes valeurs de ℎ𝐻 avec 𝜂𝑣 =
0.1 et 𝜃 = 0°. ....................................................................................................................................... 99
Figure III.A.10. Comparaison des réponses fréquentielles à l’extrémité de la poutre sandwich en porte-
à-faux de cœur viscoélastique et faces composites obtenues pour différentes valeurs de ℎ𝐻avec 𝜂𝑣 =
0.1 et 𝜃 = 0°. ..................................................................................................................................... 100
Figure III.A.11. Variations des fréquences naturelles et des facteurs de perte pour les trois premières
modes de la poutre sandwich en porte-à-faux obtenues pour différentes valeurs du rapport d'épaisseur
h3/h1 ((a) fréquence naturelle ; (b) facteur de perte).......................................................................... 101
Figure III.A.12. Comparaison des réponses en fréquence de la poutre sandwich en appuis simples
obtenues pour différentes valeurs du rapport d'épaisseur h3/h1. ........................................................ 102
Figure III.A.13. Comparaison des réponses en fréquence de la poutre sandwich en porte-à-faux
obtenues pour différentes valeurs du rapport d'épaisseur h3/h1. ........................................................ 102
Figure III.A.14. Réponses fréquentielles de déplacement de la poutre sandwich PCLD en appuis
simples sans et avec les cœurs viscoélastiques (PVB, ISD112 à 27° et ISD112 à 20°). ..................... 105
Figure III.A.15. Réponses fréquentielles de vitesse de la poutre sandwich PCLD en appuis simples sans
et avec les cœurs viscoélastiques (PVB, ISD112 à 27° et ISD112 à 20°) .......................................... 105
Figure III.A.16. Variations des fréquences naturelles et des facteurs de perte de la poutre sandwich en
appuis simples ((a) Cœur PVB ; (b) ISD 112 à 27°) .......................................................................... 106
xiii
Liste des tableaux
Figure III.A.17. Réponses fréquentielles de déplacement de la poutre sandwich PCLD à cœur PVB en
appuis simples pour différentes valeurs de rapport d’épaisseur ......................................................... 107
Figure III.A.18. Réponses fréquentielles de déplacement de la poutre sandwich PCLD à cœur PVB en
appuis simples pour différentes valeurs de rapport d’épaisseur ((a) PVB ; (b) ISD112 27°) ............. 108
Figure III.B.1. Variation des fréquences propres et des facteurs de pertes en fonction de facteur de perte
viscoélastique pour une plaque sandwich en appuis simples. ............................................................ 111
Figure III.B.2. Réponses en déplacement de la plaque sandwich en appuis simples pour différentes
valeurs de facteur de perte viscoélastique. ......................................................................................... 111
Figure III.B.3. Réponses en vitesse de la plaque sandwich en appuis simples pour différentes valeurs
de facteur de perte viscoélastique. ..................................................................................................... 112
Figure III.B.4. Variation des fréquences propres et de facteurs de pertes en fonction de rapport
d’épaisseur ℎ2ℎ3 pour une plaque sandwich en appuis simples ........................................................ 113
Figure III.B.5. Réponses en déplacement de la plaque sandwich en appuis simples pour différents
rapports d’épaisseur h2/h3 viscoélastique .......................................................................................... 114
Figure III.B.6. Réponses en vitesse de la plaque sandwich en appuis simples pour différents rapports
d’épaisseur h2/h3 ............................................................................................................................... 114
Figure III.B.7. Comparaison entre les réponses fréquentielles de déplacement de la plaque sandwich
sans et avec le cœur viscoélastique (PVB, ISD112 20° et ISD112 27°) ............................................. 116
Figure III.B.8. Comparaison entre les réponses fréquentielles de vitesse de la plaque sandwich sans et
avec le cœur viscoélastique (PVB, ISD112 20° et ISD112 27°) ........................................................ 116
Figure III.B.9. Variation des fréquences propres et des facteurs de perte pour différentes températures
du cœur viscoélastique DYAD606 .................................................................................................... 118
Figure III.B.10. Réponses fréquentielles de déplacement de la poutre sandwich à cœur viscoélastique
DYAD606 considéré pour différentes températures. ......................................................................... 118
Figure III.B.11. Réponses fréquentielles de vitesse de la poutre sandwich à cœur viscoélastique
DYAD606 considéré pour différentes températures .......................................................................... 119
Figure IV.A.1. Charge mobile définie dans le modèle d'éléments finis ............................................ 128
Figure IV.A.2. Algorithme du schéma d’intégration de Newmark ................................................... 129
Figure IV.A.3. Réponses dynamiques de la poutre sandwich à cœur viscoélastique en PVB sous
chargement mobile pour différentes valeurs de ℎ/𝐻.......................................................................... 131
xiv
Table des figures
xv
Liste des tableaux
Figure V.2. Réponses fréquentielles de la pression acoustique de la plaque sandwich avec viscoélastique
indépendant de la fréquence............................................................................................................... 154
Figure V.3. Réponses fréquentielles de la perte de transmission acoustique « STL » de la plaque
sandwich avec viscoélastique indépendant de la fréquence ............................................................... 154
Figure V.4. Réponses fréquentielles de la pression acoustique de la plaque sandwich à cœur
viscoélastique indépendant en appuis simples pour différentes valeurs de facteur de perte ((a)- ΔT=0°C,
ΔT=10°C, ΔT=20°C, ΔT=30°C)........................................................................................................ 155
Figure V.5. Réponses fréquentielles de STL de la plaque sandwich à cœur viscoélastique indépendant
en appuis simples pour différentes valeurs de facteur de perte ((a)- ΔT=0°C, ΔT=10°C, ΔT=20°C,
ΔT=30°C) .......................................................................................................................................... 156
Figure V.B.6. Réponses fréquentielles de la pression acoustique de la plaque sandwich avec
viscoélastique à cœur viscoélastique en DYAD606 ........................................................................... 157
Figure V.B.7. Réponses fréquentielles de STL de la plaque sandwich avec viscoélastique à cœur
viscoélastique en DYAD606.............................................................................................................. 157
Figure B.1. Exemple d’une transformation isoparamétrique plane ................................................... 170
Figure C.1. Exemple d’un morceau de solution ................................................................................ 173
Figure D.1. Réponses de la plaque sandwich à cœur viscoélastique ISD112 à 27° sous chargement
mobile de mouvement rectiligne et vitesse constante pour 0 < 𝑋/𝜏 < 0.4 (vibrations forcées) ........ 176
Figure D.2. Réponses de la plaque sandwich à cœur viscoélastique ISD112 à 27° sous chargement
mobile de mouvement rectiligne et vitesse constante pour 0.4 < 𝑋/𝜏 < 0.85 (vibrations forcées) .. 177
Figure D.3. Réponses de la plaque sandwich à cœur viscoélastique ISD112 à 27° sous chargement
mobile de mouvement rectiligne et vitesse constante pour 0.85 < 𝑋/𝜏 < 1 (vibrations forcées)...... 178
Figure D.4. Réponses de la plaque sandwich à cœur viscoélastique ISD112 à 27° sous chargement
mobile de mouvement rectiligne et vitesse constante pour 1 < 𝑋/𝜏 < 2.4 (vibrations libres) .......... 179
Figure D.5. Réponses de la plaque sandwich à cœur viscoélastique ISD112 à 27° sous chargement
mobile de mouvement rectiligne et vitesse constante pour 2.5 < 𝑋/𝜏 < 3.6 (vibrations libres) ....... 180
Figure D.6. Réponses de la plaque sandwich à cœur viscoélastique ISD112 à 27° sous chargement
mobile de mouvement rectiligne par morceaux (mouvement zigzag) et vitesse constante pour 0<
𝑋/𝜏 < 0.27 (vibrations forcées) ........................................................................................................ 181
xvi
Table des figures
Figure D.7. Réponses de la plaque sandwich à cœur viscoélastique ISD112 à 27° sous chargement
mobile de mouvement rectiligne par morceaux (mouvement zigzag) et vitesse constante pour 0.29 <
𝑋/𝜏 < 0.58 (vibrations forcées) ........................................................................................................ 182
xvii
Liste des tableaux
xviii
Liste des symboles, abréviations et notations
Symboles
𝜎𝑘𝑙 : tenseur de contrainte
𝜂 : facteur de perte
𝜙 : déphasage
𝐸 : module élastique
𝐸 : module de cisaillement
𝐷𝛼 : opérateur de dérivation fractionnaire
T : Température
ΔT : variation de la température
Tg : température de transition
𝛼𝑇 : facteur de décalage
𝛼𝑇 𝜔 : fréquence réduite
𝑢𝑖 : champ de déplacement
𝑤 : déplacement transversal
xix
Liste des symboles, abréviations et notations
Ω𝑠 : domaine solide
Ω𝑓 : domaine fluide
∂Ω : surface limite
(𝜉, 𝜂, 𝜁) : coordonnées locales
xx
Liste des symboles, abréviations et notations
ℎ : épaisseur
𝐿 : longueur
𝑙, 𝑏 : largeur
𝑁𝑖 : force normale
𝑀𝑖 : moment fléchissant
𝑇 : force de cisaillement
𝑁𝑇 : contrainte thermique
𝜏𝑐 : contrainte de cisaillement
𝛼𝑥 , 𝛼𝑦 , 𝛼𝑥𝑦 : coefficient d’expansion thermique
𝛼 : paramètre de chemin
𝜆 = 𝜔2 = Ω2 1 iη : pulsation
𝜀 : paramètre de précision
𝑅(𝑈, 𝜆) : erreur résiduelle
𝐵𝑚 , 𝐵𝑓 , 𝐵𝑓 , 𝐵𝑐 : matrices d’interpolations
𝐷𝑚 , 𝐷𝑓 , 𝐷𝑐 : matrices de comportment
𝑃(𝑥, 𝑡) : charge
𝛿(𝑥) : fonction de Dirac
Abbreviation
MEF : méthode des éléments finis
AM : Analyse modale
MAN : méthode asymptotique numérique
MDA : méthode de différenciation automatique
BEM : méthode des éléments de frontières
xxi
Liste des symboles, abréviations et notations
Notations
𝜕∎
∎,𝑥 = : dérivée spatiale
𝜕𝑥
𝜕∎
∎̀,𝑥 = : dérivée temporaire
𝜕𝑡
∎∗ : module complexe
∎′ : partie réelle du module
[∎] : matrice
[∎]𝑒 : matrice élémentaire
{∎} : vecteur
𝛿∎ : quantité virtuelle
̅
∎ : conjugué d’un nombre complexe
∎𝑇 : transposé
xxii
Introduction générale
Introduction générale
Le problème de vibrations est devenu de nos jours une préoccupation majeure en raison de
son effet sur le fonctionnement des systèmes et des installations. Les vibrations peuvent
engendrer des perturbations ou des instabilités dynamiques qui peuvent mener aux dégâts
structurels et fonctionnels majeurs. Par conséquent, la maitrise des phénomènes vibratoires est
une phase très importante lors de la conception des systèmes mécaniques. Le problème de
vibrations peut être résolu parfois par l’augmentation du niveau d’amortissement de la structure
par l’intégration d’un traitement amortissant ou par l’élimination des sources de vibrations.
La première est la technique d’isolation qui consiste à éviter la transmission des vibrations
d’un système à un autre sans faire éliminer les sources de vibrations. La deuxième est une
technique de suppression de vibration qui consiste à éliminer les vibrations par les vibrations
elle-même, ce principe actif est réalisé en ajoutant à la structure des capteurs et des actionneurs
reliés à un système de contrôle permettant de rétablir le comportement de la structure au
comportement régulier. Les éléments de cette technique sont généralement fabriqués à partir des
céramiques piézo-électriques transformant l’énergie de déformation en énergie électrique.
Alors que la troisième technique est une technique d’amortissement passif qui consiste à
intégrer depuis la conception des matériaux possédant des propriétés dissipatives des énergies
vibratoires. Ces matériaux possèdent un comportement viscoélastique capable de transformer
l’énergie de vibration en énergie thermique ce qui permet de dissiper l’énergie de déformation
sous forme de chaleur. Souvent, une couche viscoélastique est associée à deux parements
(métalliques, composites...) constituant une tôle sandwich à fort pouvoir amortissant tout en
gardant les caractéristiques mécaniques du sandwich. L’amortissement de la structure est apporté
par le polymère viscoélastique présent dans la tôle sandwich grâce à sa propriété de passer d’un
état peu rigide (caoutchouteux) à un état rigide (vitreux) en fonction des conditions de
chargement températures, fréquences…). Dans cette région de transition entre les deux états, les
1
Introduction générale
propriétés élastiques varient très rapidement, dans laquelle ces matériaux sont capables de
dissiper l’énergie vibratoire.
Les résultats des études effectuées ont montré que la technique d’amortissement passif est
généralement très efficace quand il s’agit de moyennes et hautes fréquences. Alors qu’en
revanche, la technique d’amortissement actif est plus efficace aux faibles fréquences. Pour cette
raison, la technique dite hybride « Passive + active) est mise en place afin de prendre en
considération toutes les plages de fréquences.
Le développement d’un système fiable et robuste de contrôle des vibrations, est devenu une
nécessité face aux exigences fonctionnelles des structures, le contrôle passif reste efficace en
matière d’atténuation des vibrations structurelles et des bruits acoustiques. La simulation
numérique des structures amorties par des matériaux viscoélastiques est une étape importante
dans la conception pour prédire l'efficacité de ce traitement amortissant. Cependant, cette étape
requiert une bonne connaissance du comportement du matériau dans la zone de transition
vitreuse. Cela a besoin de mettre en place les outils nécessaires à la modélisation de leur
comportement vibratoire.
Dans le contexte de l’étude des vibrations structurelles, notre travail vise à étudier le
comportement vibratoire et acoustique des structures sandwich à cœur viscoélastique sous
l’action des chargements dynamiques. Ce problème est devenu d’une importance fondamentale
alors que les structures de transport notamment les ponts sur lesquels circulent des véhicules ou
des trains soumissent aux charges mobiles différentes lors de passages des véhicules et trains.
Cependant, durant ces dernières décennies, les matériaux viscoélastiques ont connu une
grande évolution dans plusieurs domaines tels que le génie civil, l'aéronautique et l'industrie
automobile, en raison de leurs propriétés mécaniques spécifiques. Les matériaux viscoélastiques
contribuent dans l’atténuation des vibrations structurelles et des bruits acoustiques. Nombreuses
études ont été effectuées pour étudier leur comportement, la difficulté se réside dans leurs
propriétés qui sont dépendantes de la fréquence et de la température ainsi que dans la présence
de la non-linéarité géométrique.
2
Introduction générale
Nombreux chercheurs ont étudié l’impact de différents paramètres dynamiques tels que le
type de mouvement, la vitesse de la charge et l’amortissement de la structure mais, peu de travaux
se sont intéressés aux structures munies de matériaux viscoélastiques sous l’action de la charge
mobile. En outre, plusieurs études ont été munies pour caractériser le comportement vibratoire
des structures par matériaux viscoélastique dont nombreux chercheurs se sont intéressés à la
caractérisation des propriétés dynamiques des structures munies de matériaux viscoélastiques en
utilisant des approches analytiques et numériques. La plupart des chercheurs n’ont considéré
généralement que des modèles viscoélastiques complexes indépendants de la fréquence.
3
Introduction générale
Organisation de la thèse
Les travaux réalisés au cours de cette thèse sont présentés en cinq chapitres, organisés comme
suit :
Au chapitre I, tout d’abord, nous présentons un rappel des principales techniques
d’atténuation des vibrations avant de nous focaliser sur la technique du contrôle passif par
l’intégration des matériaux viscoélastiques dont les principales configurations sont présentées
dans ce chapitre. Ainsi, un rappel de la théorie de la viscoélasticité linéaire des matériaux non
vieillissants est effectué auquel nous présentons :
4
Introduction générale
numérique présentée y est considérée comme étant une approche typique bien destinée à résoudre
le problème de vibrations des structures sandwich à cœur viscoélastique. Cette approche est
validée à la fin de ce chapitre par une étude comparative avec les résultats des autres travaux liés.
Dans le chapitre III, une analyse des vibrations libres et forcées des structures munies par
matériaux viscoélastique est abordée en utilisant l’approche numérique développée au chapitre
II dont plusieurs configurations des sandwichs y sont examinées. Les réponses forcées obtenues
dans le domaine fréquentiel, mettent en examen plusieurs modèles de comportement
viscoélastique entre ceux indépendants de la fréquence et autres dépendants de la fréquence et de
la température.
Le comportement dynamique des différentes configurations des structures sandwich à cœur
viscoélastique sous l’action d’une charge mobile à vitesse constante est abordé dans le chapitre
IV. L’effet de l’asymétrie géométrique, des propriétés mécaniques, de la vitesse de la charge
mobile de la structure poutre sandwich ont été examinés en considérant différents modèles de loi
de comportement viscoélastique. D’autre part, pour la plaque, deux types de mouvements sont
considérés tandis que le premier est un mouvement linéaire suivant une ligne droite et le
deuxième est un mouvement en zigzag avec deux modèles de comportement viscoélastique l’un
est dépendant de la fréquence et l’autre est dépendant du temps.
Le dernier chapitre porte sur l’analyse du comportement vibro-acoustique des structures
sandwich étudiées dans les chapitres précédents auquel nous examinons les différentes réponses
acoustiques. L’objectif est d’examiner l’efficacité de traitement passif d’amortissement tout en
examinant des lois de comportement dépendant non seulement de la fréquence mais, aussi de la
température. Nous utilisons une approche numérique basée sur la méthode des éléments finis et
la méthode des éléments de frontières afin de formuler le problème du couplage entre une
structure vibrante et le domaine fluide, ainsi que sur la méthode asymptotique numérique pour
les deux types de structures poutre et plaque sandwichs.
5
Chapitre I
Etat de l’art
Sommaire
I.1. Introduction
L’étude du comportement des structures par matériaux viscoélastiques remonte au début des
années 1950 durant lesquelles des recherches ont été réalisées sur les structures poutres et plaques
7
Chapitre I
composites par matériaux viscoélastiques. Les principales recherches reviennent aux travaux de
Liénard [1], Kerwin [2], Ross [3], Mead [4], Ungar [5], Yu [6], DiTaranto [7], Mead et
Markus [8]-[9], Yan et Dowell [10], Rao [11], Rao et al [12], Kristensen et al [13], Kumara et
Singh [14], et Sher et Moreira [15].
Les vibrations des structures sandwich à cœur viscoélastique à module constant ont été
largement examinées dont nombreuses approches analytiques y ont été utilisées. En se basant sur
les travaux de Kerwin [2]et Ross [3]qui ont utilisé une expression analytique du facteur de perte
en fonction des caractéristiques d’une poutre sandwich en partant d’une équation différentielle
ordinaire de quatrième ordre, d’autres modèles analytiques ont été proposés par Ungar [5], Yu [6]
pour caractériser les propriétés amortissantes des sandwichs viscoélastiques. Ensuite,
DiTaranto [7] a défini une équation décrivant les propriétés amortissantes (pulsation amortie,
facteur de perte) avec la prise en compte de différentes conditions aux limites. Rao [11] a utilisé
le principe énergétique de Hamilton pour formuler l'équation gouvernant le mouvement de
vibration, une résolution analytique de l’équation de mouvement a été proposée par l’auteur en
caractérisant les propriétés amortissantes de poutres sandwichs. Ainsi, Rao [16] a examiné
l'application des matériaux viscoélastiques pour contrôler les vibrations et de bruit acoustique
dans les automobiles et les avions commerciaux.
Toutefois, plusieurs approches numériques basées sur la méthode des éléments finis ont été
considérées pour les structures sandwich afin de mieux prendre en compte des géométries plus
complexes. Les auteurs ont considéré plusieurs modèles cinématiques pour leurs études décrivant
le champ de déplacement de trois couches simultanément entre-autres Love [17], Mindlin [18],
Touratier [19] pour les structures plaques correspondantes aux théories d’Euler-Bernoulli et
Timoshenko adoptées aux poutres sandwichs respectivement en ce qui concerne les deux
premiers modèles. Récemment, des modèles cinématiques décrivant le champ de déplacement
« couche par couche » afin d’améliorer la précision, ce qui conduit à des modèles de type zigzag
qui combine à la fois le modèle cinématique de Rao [11]et celui de Reddy [20]. Wang et
al [21]ont utilisé une approche lagrangienne pour étudier les propriétés mécaniques sous une
charge d'impact à grande vitesse du pare-brise viscoélastique d'un avion en tenant compte des
effets de la vitesse et de la température. Le comportement dynamique a été obtenu par des
8
Chapitre I
Les matériaux viscoélastiques sont intégrés aux structures, selon différentes modalités, nous
citons dans ce cadre :
Ce type de structure, présenté par la figure I.3 est la plus simple et qui consiste à associer une
couche ou un patch du matériau viscoélastique à la surface libre de la structure. Lorsque la couche
de base soumet à une vibration en flexion, la couche d'amortissement vibre produisant une
contrainte traction-compression dans la couche viscoélastique. Selon le mécanisme de dissipation
d'énergie du matériau viscoélastique, lorsque le matériau d'amortissement soumet aux contraintes
alternées, l'énergie mécanique est convertie en énergie thermique. D’après les résultats de
recherches faites sur la configuration FLD [27]-[30], les performances d'amortissement de la
structure deviennent meilleures quand la couche viscoélastique est plus épaisse dont le facteur
de perte d'amortissement est plus important.
9
Chapitre I
Une autre approche d’amortissement utilisant des couches de contrainte aux faces de
sandwich permettant d’améliorer l’amortissement de la structure. Lorsque la structure vibre, la
couche d'amortissement viscoélastique est sous contrainte par les déplacements longitudinaux
des couches de face, cependant, la contrainte et la déformation sont déphasées, ce qui permet de
dissiper l'énergie vibratoire et réduire par conséquent les vibrations et les bruits engendrés. Cette
méthode du contrôle des vibrations, est simple et significative et elle est largement appliquée
pour améliorer les performances amortissantes des structures.
La configuration PCLD (figure I.4) est une application efficace, qui est surtout utilisée dans
l'industrie (automobile, aéronautique, génie civil…) sur laquelle plusieurs recherches ont été
menées étudiant l’efficacité de ce traitement entre autres Cai et al [31] qui ont utilisé une
approche analytique pour examiner la réponse vibratoire de la poutre en utilisant la méthode
énergétique de Lagrange. Le modèle de Mead et Markus [9] a été utilisé pour décrire les relations
cinématiques entre les trois couches. Aussi, d’autres ont utilisé ce traitement d’amortissement
avec une couche partiellement contraignante sous forme d’un patch [32]-[36].
10
Chapitre I
Le principe du contrôle actif des vibrations, donné par la figure I.5 consiste généralement à
recueillir le signal vibratoire de la structure par des capteurs, puis transmettre le signal à
l'ordinateur tout en utilisant un algorithme de contrôle avant d’envoyer un signal de commande
au système d'actionnement pour contrôler les vibrations de la structure. Le grand avantage du
contrôle actif réside à son efficacité aux basses fréquences et sa facilité à s'adapter à
l'environnement extérieur et de suivre automatiquement le changement des caractéristiques
vibratoires [37]-[44].
Ces dernières années, certains chercheurs ont proposé une nouvelle méthode de contrôle, soit
le contrôle hybride « Actif - Passif ». Le mécanisme de cette méthode (figure I.6) consiste à
ajuster activement la caractéristique d'amortissement de la structure grâce à l'amortissement de
la couche de contrainte. Lorsque la structure vibre, la couche d'amortissement viscoélastique
subit à une déformation par cisaillement et dissipe l'énergie de vibration. Pendant ce temps, la
couche piézoélastique augmente la déformation en cisaillement de la couche viscoélastique pour
11
Chapitre I
Un matériau viscoélastique est un polymère dont son module est fortement lié à la fréquence
ou la température, cette propriété offre à la structure une grande capacité à atténuer les vibrations.
Toutefois, nombreuses études ont été menées pour étudier les variations des paramètres
d’amortissement tels, que le module et le facteur de perte des matériaux viscoélastiques, les
études ont montré que ces types de matériaux peuvent prendre quatre états différents, selon la
température et la fréquence. Les zones distinctes du comportement des matériaux viscoélastiques
sont identifiées par la figure I.7. La première est le domaine vitreux qui correspond aux plus
faibles températures et aux plus hautes fréquences et qui se caractérise par un module de stockage
atteignant sa valeur maximale et un facteur de perte atteignant sa valeur minimale. Dans cette
12
Chapitre I
zone, l'énergie thermique est faible, empêchant ainsi tout mouvement d'extension ou de rotation
des chaînes macromoléculaires de polymère. Par conséquent, le polymère est rigide caractérisé
par un module de stockage élastique élevé et un facteur de perte faible. La deuxième zone
correspond au domaine de transition viscoélastique qui se caractérise par une diminution rapide
du module stockage et une augmentation importante du facteur de perte qui atteint sa valeur
maximale au milieu de cette zone correspondant à la température de transition vitreuse du
polymère traduisant sa grande capacité d’amortissement. L’analyse moléculaire montre que
l’amortissement viscoélastique, dû principalement à l’énergie dissipée sous forme de chaleur lors
du glissement des segments de chaînes, ce qui se traduit par un retard à la déformation.
Cependant, entre le domaine de transition viscoélastique et le domaine d’écoulement, le
polymère présente un état caoutchoutique, où le module de stockage et le facteur de perte sont
relativement faibles et leurs variations avec la température et la fréquence sont très faibles. Le
dernier domaine correspond au domaine d’écoulement dans lequel le matériau se comporte
comme un fluide visqueux dont la déformation devient permanente.
13
Chapitre I
Essai de fluage
σ ε
σ0 ε0
t1 t t1 t
Essai de relaxation
14
Chapitre I
déformation. On observe une diminution de la contrainte dans le temps qui tend vers une limite
asymptotique. Ce type d’évolution de la contrainte s’appelle la relaxation de contrainte.
ε σ
ε0 σ0
t1 t t1 t
Figure I. 9. Essai de relaxation
Pour un essai de recouvrance (figure I.10), un créneau de contrainte entre deux instants t1 et
t2 est appliqué tout en suivant l’évolution de la déformation. Cette expérience est similaire à
l’essai de fluage pour t1≤t<t2 auquel on observe un saut instantané de la déformation à t=t1 avant
qu’elle de commencer de croitre pour t<t2, alors pour t>t2 la déformation décroît jusqu’à qu’elle
prend son état initial.
15
Chapitre I
n
16
Chapitre I
Cette intégrale est considérée comme la somme d’une réponse de l’élasticité instantanée et
une autre de l’effet de mémoire qui tient compte l’histoire de la déformation ε𝑘𝑙 .
R ijkl t ,
ij R ijkl t , t ε kl t ε kl d I.5
t
Le même principe peut être appliqué pour la déformation :
I t ,
t , t kl t kl d
t
εij I
ijkl
(I.6)
ijkl
Pour un matériau non vieillissement, cela veut dire que ses paramètres n’évoluent pas dans
le temps signifiant que sa réponse à l’instant t pour une sollicitation à l’instant t est la même
réponse, à l’instant t+ pour la même sollicitation appliquée à l’instant t+ les fonctions de fluage
et de relaxation ne dépendent donc que de l’intervalle de temps séparant l’effet et la cause [51] :
R ijkl t ,
ij R ijkl 0 ε kl t ε kl d
t
(I.7)
I ijkl t ,
0 kl t kl d
t
εij I (I.8)
ijkl
ij s R * s ijkl ε kl s (I.9.b)
Où s est le paramètre de l’espace de Laplace et de Carson. 𝜀(𝑠), 𝜎(𝑠), 𝜀 ∗ (𝑠) et 𝜎 ∗ (𝑠) sont la
transformée de Laplace et de Carson de 𝜀(𝑡) et de 𝜎(𝑡), respectivement. ℜ∗ (𝑠) et 𝜀𝑖𝑗 (𝑠) sont
respectivement les tenseurs d’ordre 4 de module et de complaisance opérationnels. Grâce aux
hypothèses de linéarité et de non-vieillissement du comportement viscoélastique, les équations
Eq. (I.9.a) et Eq. (I.9.b) sont analogues aux relations contrainte-déformation des matériaux
élastiques décrits par la loi de Hooke généralisée.
17
Chapitre I
L’une des méthodes utilisées pour étudier et caractériser les propriétés mécaniques
viscoélastiques, consiste à appliquer une charge harmonique dynamique sinusoïdale et mesurer
sa réponse.
ij s ij 0cos t ij 0 Re eit (I.10)
ij ij 0 R ijkl t , Re iei d
t
(I.11)
La partie réelle du module ℜ′𝑖𝑗𝑘𝑙 , est le module de stockage qui traduit la capacité stockage
de l’énergie de déformation et la partie imaginaire ℜ′′ est le module de dissipation qui traduit la
18
Chapitre I
On peut remarquer que la contrainte qui résulte du chargement harmonique, est également
harmonique Eq. (I.19), de même fréquence, mais déphasée par rapport à la déformation
harmonique Eq. (19) [51]. Le facteur de perte est défini par l’angle de ce déphasage 𝜙/𝜔 (figure
I.12).
19
Chapitre I
Modèles rhéologiques
un modèle rhéologique élémentaire consiste à faire une combinaison d’éléments mécaniques
permettant de représenter le comportement mécanique du matériau (figure I.13).
(a) (b)
Pour le ressort qui représente un élément élastique idéal, la contrainte est proportionnelle à la
déformation et sa relation est définie par la loi de Hook généralisée, le module E correspond au
module élastique du matériau. Cependant, l’amortisseur représente un élément visqueux
dissipatif dont la contrainte est proportionnelle à la vitesse de déformation et sa relation est
définie par la loi de frottement visqueux de Newton. Les modèles rhéologiques standards
fréquemment utilisés sont présentés ci-dessous avec leurs lois de comportement et leurs modules
opérationnels.
Modèle de Maxwell
𝜂 𝑬𝟎 η 𝑖𝐸0 𝜔𝜂
𝜎 (𝑡) + 𝐸 𝜎̇ (𝑡) = 𝜂𝜀̇ (𝑡) 𝑅∗ (𝑖𝜔) = 𝐸
0 0 +𝑖𝜔𝜂
η
Figure I. 15. Modèle de Kelvin-Voigt
𝑬𝟎
Modèle de Zener
𝑬𝟏
𝐸1 (𝐸0 +𝑖𝜔𝜂)
(𝐸0 + 𝐸1 )𝜎(𝑡) + 𝜂𝜎̇ (𝑡) = 𝐸1 (𝐸0 𝜀(𝑡) + 𝜂𝜀̇(𝑡)) 𝑅 ∗ (𝑖𝜔) =
𝐸0 +𝐸1 +𝑖𝜔𝜂
η
Figure I. 16. Modèle de Zener
20
Chapitre I
𝑬𝟏 𝑬𝟐 𝑬𝟑 𝑬𝒏
𝑬𝟎
η0
η1 η2 η3 ηn
𝛈𝟐 𝛈𝐧
𝛈𝟏
Modèles fractionnaires
La construction des modèles fractionnaires consiste à intégrer aux modèles standards
linéaires un élément mécanique « spring-pot » (figure I.19) ayant un comportement rhéologique
entre un ressort élastique et un amortisseur visqueux.
𝑬; 𝝉; 𝜶
21
Chapitre I
ED 0 1 (I.20)
d 1 d t
D
d (I.21)
dt Γ 1 dt 0 t
+∞
où Γ(𝑧) = ∫0 𝑒 −𝑥 𝑥 𝑧−1 𝑑𝑥 est la fonction gamma définie pour chaque réel strictement positive
z.
N
d k N
d k t
t ak E k k
t b (I.22)
k 1 dt k k 1 dt
le module complexe correspondant est donné par :
E k 1bk i
N k
R
*
(I.23)
1 k 1ak i
N k
22
Chapitre I
𝑬𝟏
𝑬𝟏 ; 𝝉; 𝜶
Figure I. 21. Modèle de Kelvin-Voigt fractionnaire
Modèle de Zener 𝑬𝟎
𝑬𝟏
𝐸𝑞 𝐸𝑞 (1+𝑖𝜏𝜔)𝛼
𝜎 + 𝜏 𝛼 𝐷 𝛼 (𝜎) = 𝐸𝑞 𝜀 + 𝐸𝑞 𝜏 𝛼 𝐷 𝛼 (𝜀) 𝑅 ∗ (𝑖𝜔) = 𝐸
𝑞
𝐸0 1+𝐸 (𝑖𝜏𝜔)𝛼
0
𝑬𝟎 ; 𝝉; 𝜶
Figure I. 22. Modèle de Zener
𝐸1 𝐸2
avec 𝐸𝑞 =
𝐸1 +𝐸2
Les essais expérimentaux [52] montrent que pour des températures différentes, et une
fréquence fixe, la variation du module de stockage et du facteur de perte en fonction de la
température présente la même allure que ceux obtenu à température fixe en faisant varier la
fréquence (figure I.23), ce qui permet de superposer les courbes à différentes températures. La
superposition s’effectue par une translation horizontale selon un facteur dit facteur de décalage
𝛼 𝑇 correspondant à un changement d’échelle de fréquence dite fréquence réduite 𝛼 𝑇 𝜔 sur la
nouvelle courbe superposée. Cela permet de relever intuitivement les fréquences réduites pour
chaque température.
Figure I. 23. Superposition des courbes du module de stockage et du facteur de perte à différentes températures T = 0 – 50°C
du matériau viscoélastique BI2F [52]
23
Chapitre I
L’obtention d’une nouvelle valeur du module complexe ou du facteur de perte pour une
fréquence ωj et une température Tk quelconques se fait en trois étapes, dont la première consiste
à repérer un point d’intersection P de la droite horizontale ωj et de la droite oblique isotherme
Tk, Ensuite, l’abscisse du point P fournit la valeur de αTωj correspondant à Tk, Finalement, les
valeurs de 𝐸′ et 𝜂 sont obtenus sur les courbes maitresses en fonction de (Tk).
Par conséquent, l’évolution du module et du facteur peut être exprimée mathématiquement
par l’interpolation des courbes de maîtresse sous une forme générale donnée par :
24
Chapitre I
C10 T T0
log T (I.25)
C20 T T0
𝐶10 et 𝐶20 sont des constantes dépendantes du matériau sont données en fonction de la température
de référence 𝑇0. Ces deux constantes peuvent prendre des valeurs dépendamment du type de
polymère, généralement, les études ont montré que ces valeurs peuvent être restreintes dans une
plage définie en fonction de la température de transition Tg [52] - [54].
Ea 1 1
log αT (I.26)
R T T0
où, 𝑹 = 𝟖. 𝟑𝟏𝟒 × 𝟏𝟎−𝟑 𝒌𝑱𝒎𝒐𝒍−𝟏 𝑲−𝟏 est la constante des gaz parfaits et 𝑬𝒂 est l’énergie
d’activation de la relaxation.
25
Chapitre I
ces modèles précités négligent les contraintes de cisaillement dans les trois couches en raison de
leurs faibles effets par rapport aux autres types de contraintes.
ui x, z, t u0i x, t z zi w, x I.27
Le cisaillement de la couche viscoélastique, dû à sa faible rigidité et aux mouvements
longitudinaux des couches de faces assure la fonction d’amortissement ou de dissipation de
l’énergie de déformation à la structure. Pour cette raison la prise en considération de la
déformation de cisaillement de la couche viscoélastique est très nécessaire pour mieux estimer
les propriétés ou la réponse de la structure sous différents chargements. Plusieurs modèles
cinématiques, y ont été proposés permettent de prendre en considération la déformation de
cisaillement, ces modèles qui reposent sur la formulation continue permettent de traiter les
couches du sandwich comme étant une seule couche ce qui permet de simplifier l’étude de ces
structures (figure I.26), ces modèles y ont montré une grande efficacité notamment pour les
structures stratifiées multicouches. Parmi ces modèles, on peut citer les modèles linéaires de
Timoshenko [57] et de Mindlin [58] correspondants respectivement aux structures poutres et
plaques, donnés par l’Eq.(I.28), les modèles cubiques de Reddy [59], Reissner [64] Eqs.(I.29-
I.31) et le modèle sinusoïdal de Touratier [60] Eq.(I.32). Ces approches sont mieux adaptées pour
décrire la continuité de déplacement aux interfaces de contacts entre les différentes couches de
la structure.
ui x, z, t u0i x, t z zi w, x z I.28
26
Chapitre I
4z3
ui x, z, t u0i x, t z zi w, x z 2 I.29
3h
1 zh2 z 3
ui x, z, t u0i x, t z zi w, x I.30
2 4 3
5 4z3
ui x, z, t u0i x, t z zi w, x z 2 I.31
4 3h
h z
ui x, z, t u0i x, t z zi w, x sin I.32
h
Récemment, des modèles cinématiques décrivant le champ de déplacement « couche par
couche » ont été développés afin d’améliorer la précision, ce qui conduit à des modèles de type
zigzag qui combine à la fois plusieurs modèles cinématiques. Ces modèles cinématiques zigzag
sont relativement coûteux vu le nombre de variables indépendants dans chaque couche, pourtant,
cela permet l'obtention de résultats plus précis. Nombreux auteurs ont utilisé ce genre des
modèles, tels que Bilasse [26], Z. Huang et al [61], Trindade [62], Daya et al [25]. Pour un
sandwich viscoélastique, plusieurs modèles zigzags ont été proposés en utilisant généralement
un modèle linéaire aux faces de sandwich et un modèle d’ordre supérieur pour le cœur, parmi ces
modèles, celui de Rao [11] (figure I.27) basé sur le théorème de déformation de premier ordre
dont la théorie d’Euler-Bernouilli est appliquée aux faces composites de sandwich et la théorie
de Timoshenko au cœur viscoélastique dont le champ de déplacement est donné par l’Eq.(33).
Euler-Bernoulli model
Timoshenko model
Euler-Bernoulli model
Figure I. 27. Champ de déplacement donné par le modèle linéaire en zigzag de Rao [11]
ui x, z, t u0i x, t z zi w, x i 1 et 3
u2 x, z, t u02 x, t z w, x z i2 (I.33)
Cependant, Bickford [63] a présenté un nouveau modèle (figure I.28) dont le champ de
déplacement du cœur est d’ordre supérieur à celui considéré par Rao [11], le modèle utilisé par
27
Chapitre I
Bickford [63] applique la théorie de modèle cubique de Reddy [59] toute en gardant le modèle
linéaire d’Euler-Bernouilli aux faces du sandwich.
Euler-Bernoulli model
Reddy model
Euler-Bernoulli model
Figure I. 28. Champ de déplacement donné par le modèle linéaire en zigzag de Bickford [63]
ui x, z, t u0i x, t z zi w, x i 1, 2
4z3
u2 x, z, t u02 x, t z zi w, x z 2 i 3 (I.34)
3h
Plusieurs théories ont été aussi développées afin d’améliorer l’efficacité du modèle zigzag
[64-69]. Par conséquent, la plupart des résultats ont montré que ces modèles de type zigzag
présentent plus de précision avec un coût de calcul acceptable. Ainsi, la simplicité
d’implémentation et d’adaptation permet l’amélioration des théories actuelles pour l’étude des
structures sandwich [[70]-[73]]. Cela permet d’établir des modèles plus complexes en combinant
différents modèles d’ordre supérieur. Les travaux de Ossadzow [75], Ossadzow [76]et
Karama [76]ont utilisé le modèle sinusoïdal de Touratier [60] dans le modèle zigzag afin
d’améliorer les effets de cisaillement. Ensuite, une fonction exponentielle y a été intégrée par
Afaq et al. [77]-[78] donnant un modèle cinématique non linéaire d’ordre supérieur (figure 1.29).
Les structures de transport notamment les ponts sur lesquels circulent des véhicules ou des
trains sont subis à des charges mobiles différentes lors de passages des véhicules. Ce type de
charge entraine généralement des vibrations qui peuvent réduire les performances de la structure.
Ce problème est d'une importance fondamentale ce qui permet d’attirer l’attention de chercheurs
28
Chapitre I
29
Chapitre I
résultats de ces études montrent que la vitesse de la charge mobile et l’orientation de la fibre
influent significativement sur le comportement dynamique de la poutre.
Cependant, face aux effets indésirables de différents types de chargement et de contraintes
structurelles et acoustiques liées aux exigences de fonctionnement des structures, l'absence d'une
approche d'évaluation numérique pour analyser le problème des vibrations des structures
sandwich à amortissement passif est la principale motivation de la présente recherche.
L'importance et l'efficacité des matériaux viscoélastiques "Passive Vibration Control" nécessite
d'explorer plus profondément le comportement des matériaux viscoélastiques dans la zone de
transition viscoélastique, en particulier sous l'effet de la charge mobile.
Une grande partie des sources sonores rayonnent de l'énergie par l'action des surfaces solides
vibrantes sur le fluide environnant. Certaines de ces sources, telles que les pianos, les haut-
parleurs de radio, sont généralement souhaitables à l’inverse des blocs de moteurs à combustion
interne et des roues de train qui sont évidemment non souhaitables. Ainsi, le confort acoustique
des passagers est une considération majeure dans la conception des véhicules routiers,
ferroviaires, maritimes et aériens, dans lesquels une grande partie du bruit est due à la vibration
de l'enveloppe structurelle. Par ailleurs, le bruit rayonné vers l'extérieur par les véhicules au sol
est une préoccupation environnementale majeure. De même le concepteur de navires militaires
et de bateaux de pêche doit réduire le rayonnement sonore des structures de la coque afin de
minimiser, respectivement, les chances de détection ou la perturbation des poissons.
Par conséquent, le contrôle des bruits et des vibrations, est devenu l’une des préoccupations
majeures. Les ingénieurs doivent comprendre les mécanismes physiques du rayonnement sonore
par les structures vibrantes afin de spécifier et de mettre en œuvre de manière appropriée les
mesures rentables pour le contrôler.
Le problème vibro-acoustique implique une interaction entre la structure et les fluctuations
de pression des fluides environnants. La réponse de la structure est régie par l'excitation, les
conditions limites et le couplage avec le milieu environnant. Ce couplage dépend à la fois de la
géométrie du problème et les propriétés de la structure et des fluides.
Bien que le mécanisme de génération du son par vibration de surface, l'efficacité du
rayonnement par rapport à l'amplitude des vibrations varie largement d’une source à autre. Pour
30
Chapitre I
qu'une surface vibrante puisse émettre un son, il doit non seulement être capable de modifier la
densité du fluide avec lequel est en contact, mais doit le faire de manière à produire des
changements de densité dans la région du fluide éloignée de la surface. Les problèmes couplés
fluide-structure peuvent être classés en trois catégories principales, selon le type d’interaction
fluide-structure. Le premier est le problème intérieur qui concerne le comportement d’interaction
d’un fluide délimité par une surface (figure I.30.a), ce type de problème est similaire à un
problème de cavité. Le deuxième appelé problème extérieur concerne l’interaction d’une
structure avec fluide infini ou semi-fini (figure I.30.b), le problème est similaire au processus de
rayonnement sonore par vibration des surfaces solides, ainsi que la charge fluidique associée.
Le troisième problème vibro-acoustique traite l’interaction intérieure et extérieure fluide-
structure, on peut citer l’exemple du rayonnement sonore d'une structure couplée à une cavité
interne.
Le rayonnement sonore des plaques rectangulaires est important dans de nombreux domaines
de l'ingénierie où le contrôle du bruit est nécessaire. En particulier, la puissance sonore qui est
considérée l’un des paramètres appropriés permettant de caractériser la force du son généré par
une structure vibrante. Les rayonnements acoustiques des plaques rectangulaires ont fait l’objet
de nombreuses recherches. Le but de ces études, est généralement la prédiction et la réduction
des bruits, émettent par la structure, ce qui nécessite une meilleure compréhension des
phénomènes du couplage entre fluide et structure, par conséquent, nombreuses méthodes et
logiciels de prédiction sont développées. La première recherche dans ce domaine remonte à 1945,
31
Chapitre I
où Rayleigh [90] a présenté une approche d’analyse du champ sonore rayonné par une plaque
plane placée dans un baffle infini. Ensuite, Wallace [91] a présenté une solution asymptotique
pour les rendements de rayonnement de différents modes de plaque à basse fréquence. Lomas et
Hayek [92] ont développé la solution de la fonction de Green pour les vibrations en régime
permanent d'une plaque rectangulaire appuis élastiquement et couplée à un milieu acoustique
semi-infini. Cremer et Heckl [93]ont présenté le rayonnement sonore d'une source planaire en
utilisant l'approche de la transformée de Fourier dans l'espace K (espace du nombre d'onde).
Williams [94] a proposé une expansion en série de la puissance ascendante de l'onde numéro 𝐾
pour la puissance acoustique rayonnée à partir d'une source plane. Le rayonnement sonore des
disques circulaires et annulaires minces a été examiné par Lee et Singh [95] en se concentrant
sur les modes de vibration en flexion. Liu et Li [96] ont étudié la perte de transmission acoustique
(STL) des panneaux sandwich constitués de matériaux orthotropes dans un environnement
thermique. Les résultats de la simulation numérique sont validés par des essais expérimentaux.
Ainsi, Tao et al [97] ont présenté une expression mathématique pour le rayonnement sonore
d’une fine plaque isotrope infinie couplée à des fluides et soumis à une excitation multipoint.
Lin [98] a étudié le champ acoustique de rayonnement des plaques minces sous la sollicitation
de flexion. La distribution de la pression de la source dans un champ acoustique lointain a été
obtenue de manière analytique et son champ acoustique proche a également été calculé
numériquement. Zhou et al [99] ont calculé la pression acoustique à champ lointain à partir d'une
fine plaque inhomogène stratifié en contact avec un fluide soumis à une excitation ponctuelle.
Musha [100] a présenté une nouvelle méthode de calcul numérique basée sur l'intégrale de
Rayleigh, qui a été utilisée pour obtenir le diagramme de rayonnement instantané en utilisant la
transformée en ondelettes harmoniques pour une plaque vibrante. Lee et Singh [101]ont présenté
deux solutions analytiques pour le rayonnement modal de la pression acoustique en champ
lointain à partir des modes structurels radiaux d'un disque annulaire épais à limites libres. Ainsi,
les rayonnements dans des conditions limites différentes, sont pris en compte par nombreuses
recherches [102]-[103].
Récemment, de nouvelles recherches ont été établies portant des approches numériques qui
permettent de résoudre des problèmes plus complexes. Arenas et al [104] ont analysé le champ
sonore d'un transducteur de forme elliptique sur un baffle infini. Legault et Atalla [105] ont
32
Chapitre I
développé un modèle théorique pour analyser la transmission du son à travers un double panneau
périodiquement couplé à des isolateurs de vibrations. Abdelhamid et al [106] ont combiné des
approches physiques et subjectives pour évaluer la qualité sonore des plaques impactées. Jeyaraj
et al [107] utilisent des méthodes de simulation combinant la méthode des éléments finis avec la
méthode des éléments de frontière (FEM-BEM) pour étudier les réponses vibratoires et
acoustiques des plaques isotropes dans un environnement thermique. Ils ont constaté que les
fréquences naturelles diminuent tandis que la réponse de déplacement de la structure augmente
avec l'augmentation de la température de la plaque. Jeyaraj et al [108] ont étudié numériquement
les caractéristiques vibratoires et la réponse acoustique d'une plaque composite renforcée de
fibres dans un environnement thermique avec la considération des propriétés d'amortissement
intrinsèques du matériau composite. Ils ont constaté que la réponse vibratoire de la structure
diminue avec l'augmentation uniforme de la température pour les deux matériaux Verre-Epoxy
et PEEK IM7, mais le rayonnement acoustique de la plaque ne diminue que marginalement en
raison de l'interaction entre la rigidité réduite et l’amortissement amélioré. Une approche
analytique a été présentée par Geng et Li [109] pour étudier les vibrations et les caractéristiques
acoustiques d'une plaque mince rectangulaire isotrope dans des conditions thermiques uniformes,
les mêmes conclusions précédentes y ont été également tirées d’où la première fréquence
naturelle est plus sensible aux changements de température. Les résultats montrent que les
fréquences naturelles du panneau diminuent et que les pics du STL ont tendance vers les basses
fréquences avec l'augmentation de la température. De plus, Zhao et al [110] ont présenté une
étude numérique des plaques composites à fibres renforcées dans un environnement thermique
en utilisant la théorie des plaques laminées classiques. Cette étude a montré que les amplitudes
de résonance diminuent avec l'augmentation de la température. Par ailleurs, Geng et Li [111] ont
étudié les caractéristiques vibratoires acoustiques d'une plaque mince isotrope rectangulaire dans
des conditions thermiques. L’étude montre que les courbes de réponse des vibrations des plaques
et l'efficacité du rayonnement acoustique se dirigent vers les basses fréquences avec
l'augmentation de la température de la structure, ce qui est vérifié avec des simulations établies
par l'approche combinée de la méthode des éléments finis (FEM) et la méthode des éléments finis
de frontière BEM « Boundary element method».
33
Chapitre I
Cependant, nous remarquons que peu de travaux de recherche traitant le comportement vibro-
acoustique avec matériaux viscoélastiques, nous citons dans ce cadre les travaux de Wang et
al. [112], Foin et al. [113], Larbi et al. [114], Assaf [115], Assaf et al. [116] et Abdennadher et
al. [117]. Ainsi, Das et al [118] ont contrôlé la réponse acoustique d'une plaque composite
laminée à l'intérieur de la cabine d'un véhicule afin d'améliorer le confort des opérateurs dans les
véhicules de transport. Narayanan et ses collaborateurs [119] ont établi une expression de forme
fermée du facteur de perte de transmission en matière de paramètres du cœur viscoélastique pour
un panneau sandwich amorti. Dans ce contexte, Wang et al [112]] ont examiné les effets pondérés
du panneau sandwich afin d'optimiser la perte de transmission acoustique pour une gamme de
fréquences élevées.
La plupart de ces recherches ont montré l’efficacité des matériaux viscoélastiques en matière
d’amortissement des radiations acoustiques. Les avantages d’amortissement passif sont
nombreux en raison de la facilité de mise en œuvre avec un coût très réduit, ainsi qu’il n’exige
aucune surveillance. Ce traitement d’amortissement peut être considéré comme une solution la
plus adéquate qui peut faire face aux effets indésirables de différents types de chargement et de
contraintes structurelles et acoustiques, liées aux exigences de fonctionnement des structures.
L’étude vibro-acoustique implique une analyse du problème structurel couplé à un fluide fini
ou infini. La réponse des équations fondamentales régissant la propagation linéaire des ondes
sonores dans les domaines des fluides (acoustique) est un couplage de l'ensemble du système
soumis aux excitations et aux conditions limites données [120].
34
Chapitre I
Dans ce cas, la pression acoustique 𝑝(𝑥, 𝑡) dans le volume extérieur de densité 𝜌0 , due à une
distribution de la source acoustique doit satisfaire l'équation d’Helmholtz:
p k 2 p 0 I.35
où, 𝑘 = 𝜔/𝑐0 est le nombre d’onde avec 𝑐0 est la vitesse du son.
Cette équation doit être associée aux conditions initiales et limites. Tout d’abord, la continuité
des déplacements normaux à la surface d’interaction fluide-structure ∂Ωs donnée par l'équation
d'Euler :
p
I.36
s
0 2u.n sur
n
En outre, pour les problèmes de rayonnement acoustique dans un milieu infini. La condition
de rayonnement de Sommerfeld doit être satisfaite pour assurer la disparition de radiations à
l’infini :
p
lim r ikp 0 (I.37)
r
r
Pour modéliser du couplage fluide-structure, il est nécessaire de tenir compte de la charge du
fluide. En d'autres termes, nous devons résoudre simultanément les équations couplées régissant
la propagation des ondes dans le fluide et dans la structure. La formulation de l’équation du
couplage fluide-structure se diffère en fonction le type d’interaction entre les deux domaines.
35
Chapitre I
Pour les problèmes intérieurs dans lesquels nous nous intéressons au comportement du domaine
fluide délimité par la structure. En fait, la formulation appliquée pour résoudre ce problème, est
généralement basée sur la méthode des éléments finis pour les deux domaines du fluide et la
structure. Pour les problèmes extérieurs, le domaine fluide est fini (délimité) ou infini (non
délimité). La résolution du problème couplé, est réalisée par une méthode mixte des éléments
finis pour la structure et une méthode des éléments finis de frontière pour le fluide. Ainsi, pour
les problèmes intérieurs/extérieurs, la structure couple deux domaines de fluide dont le domaine
intérieur est limité par la surface de la structure et le domaine extérieur est non limité. Une
combinaison des méthodes des éléments finis pour la structure et le domaine intérieur du fluide
et des éléments finis de frontière pour le fluide extérieur, y est appliquée.
Dans ce travail, nous nous intéressons aux problèmes extérieurs auxquels nous étudions les
radiations acoustiques engendrées par une structure vibrante dans un domaine fluide infini ou
semi-fini. Nous présentons dans la section suivante, la formulation de l’équation gouvernant le
comportement vibro-acoustique.
Nous considérons une structure Ω𝑠 couplée à un milieu extérieur, illustrée par la figure I.32,
la surface limite de la structure, est dénotée 𝜕𝛺 et 𝑛̅𝑠 est la normale à la surface 𝜕𝛺 et 𝑛̅𝑓 est la
normale extérieure à Ω𝑓 . Le problème du rayonnement acoustique est défini par le système
d'équations suivant :
36
Chapitre I
Ʈ∞
ns
Ωf R
nf Ωs
∂Ωs,f
∂Ωs,u
2 2u 0
s dans s I.38.a
I.38.b
s
.n F sur s ,u
I.38.c
s f
.n p n dans s, f
2
p x k p x 0 f I.38.d
2
sur
p
I.38.e
f
0 2un sur s , f
n
où 𝜎, 𝑢 et 𝜌𝑠 est le tenseur de contrainte, le champ de déplacement de la structure et la densité de
la structure, respectivement. 𝑝 est la pression acoustique dans le fluide, 𝜌0 est la densité du fluide,
et 𝑘 = 𝜔⁄𝑐0 est le nombre d'ondes.
37
Chapitre I
: dV u udV
f
s
2
p n udS F udS (I.39.a)
s s s s ,u
La discrétisation par les éléments finis des équations I.39.a et I.39.b, implique la discrétisation
du domaine en sous-domaines 𝛺 ≈ 𝛺1 ∪ 𝛺2 . . .∪ 𝛺𝑛𝑒 en utilisant un système des coordonnées
locales sans dimension, associé à chaque élément (Annexe B).
T : , , r x , , , y , , , z , , p (I.40)
La construction de cette transformation repose sur la représentation des nœuds par des
fonctions de forme pour une transformation isoparamétrique (Annexe B) :
38
Chapitre I
x , , N i , , x e (I.41.a)
y , , N , , y
i
e
(I.41.b)
z , , N , , z
i
e
(I.41.c)
Pour évaluer l'intégrale élémentaire par rapport à l'élément de référence, nous devons
exprimer les gradients et les intégrales de l'espace physique en termes des coordonnées locales
(Annexe B). La géométrie différentielle classique fournit les outils nécessaires pour effectuer
cette transformation. Enfin, le volume élémentaire 𝑑𝑉 de l'élément physique 𝛺 𝑝 est lié au volume
élémentaire de l'élément de référence 𝛺 𝑟 par :
J d d d (I.42)
Par conséquent, les dérivées de ces matrices dans le repère physique sont données par:
39
Chapitre I
N1 N1
x
N1 1 N1
J (I.45)
y
N1 N1
z
La discrétisation par les éléments finis des équations I.39.a et I.39.b, implique les matrices et
les vecteurs suivants:
: ( u)dV : ( u) J d d d u B C B U n J d d d
T
p p p
s s s
u K U n
e
(I.46)
N1 N 2 N n
...
N N 2 N n
B 1
... (I.47)
N1 N 2 N n
...
u udV u u J d d d u Nu Nu U n J d d d
T
s s s
pe pe pe
s s s
u M U n (I.48)
F udS F u (II.49)
s ,u
u N p N u Pn
f f T
p n udS p n u J d d
pe pe pe
s s s
u Cup Pn
e
(I.50)
1 1
p pdV p N p , , N p , , Pn J d d d
T
0 pe 0 pe
f f
1
p S S Pn J d d d p H e Pn
T
(I.51)
0 pe
f
1 1
p N p N p Pn J d d d p Q Pn
T
p pdV (I.52)
0c02 s
0c02
pe
f
40
Chapitre I
K Cup M 0 U 0
- 2 n (I.56)
0 H Cup Q Pn 0
T T
À partir de la deuxième ligne de l'équation (I.56), il est possible d'écrire Irons [121] :
Pn 2 H
1
C
up
U n QPn
T
(I.57)
41
Chapitre I
s = s K s
T
(I.59.a)
I = s K s
T
(I.59.b)
U n s U nm (I.59.c)
Fn s F (I.59.d)
K f
T
f = f (I.60.a)
I = f H f
T
(I.60.b)
Cup ,n = s Cup f
T
(I.60.c)
Pn s Pnm (I.60.e)
s 2 - 2 I Cup ,n
U n Fn
1 2 P 0 (I.61)
- I
T
Cup ,n 2 f n
42
Chapitre I
La résolution du système réduit est maintenant moins couteuse, mais la base modale du
couplage fluide-structure est généralement tronquée à un nombre fini de modes. Cela introduit
des erreurs dans la solution dans certains cas spécifiques.
Dans cette section, une formulation variationnelle est associée au problème du rayonnement
acoustique extérieur provoqué par une structure vibrante (figure I.33). La formulation utilisée est
basée sur le couplage de la méthode des éléments finis (MEF) pour la structure et la méthode des
éléments de frontière (BEM) pour le fluide extérieur. Soit Ω le domaine externe d'une structure
vibrante dont la limite est désignée par ∂Ω (figure II.4). Le problème du rayonnement acoustique
est régi par le système d'équations (I.62).
2 p k 2 p 0 sur f I.62.a
p
I.62.b
f
0 un sur s , f
2
n
p
lim r ikp 0 (I.62.d)
r r
La représentation intégrale directe de l'équation (I.62) peut être écrite comme suit :
f
p 2G ( x, y )dV k 2 G( x, y) pdV
f f
p ( x, y ) dV 0 (I.63.b)
43
Chapitre I
𝐺𝑥𝑦 est la solution élémentaire appelée fonction de Green représente le champ sonore en un point
y dû à une source ponctuelle unitaire située en un point x et satisfaisant toutes ou certaines
conditions limites homogènes (Dirichlet, Newmann, Robin) du problème (Annexe A). Il est
important de noter que la fonction de Green peut être choisie en fonction du problème qui nous
intéresse. En pratique, ce choix implique une fonction de Green qui peut être exprimée
analytiquement et qui permet de simplifier les équations.
La fonction de Green satisfait l'équation d’Helmholtz pour une source ponctuelle unitaire est
représentée par la fonction delta de Dirac:
2G k 2G 0 sur f I.64.a
G
I.64.b
f
0 2un sur s , f
n
G
lim r ikG 0 (I.64.e)
r
r
La solution de l’équation (I.64) est une fonction de Green donnée par Wright [122]:
eikr
G( x, y ) I.65
4 r
avec 𝑟 = √(𝑥1 − 𝑦1 )2 + (𝑥3 − 𝑦3 )2 + (𝑥3 − 𝑦3 )2 est la distance entre le potentiel créé au point
y et la source acoustique située au point x. Nous établissons des équations d'intégrale aux
frontières pour le problème extérieur en se basant sur le champ de pression et son gradient normal
en faisant la soustraction de deux équations Eq.(I.63.a- I.63.b), cela conduit à l'identité :
G( x, y ) p( y )
p ( x)
f
p( y )
n
dS G( x, y )
n
dS (I.67)
f
Cette équation est une représentation intégrale du champ 𝑝(𝑥) qui est également connue sous
le nom d'équation de Kirchhoff Helmholtz. Les équations I.66 et I.67 permettent de calculer le
champ de pression total, résultant de la diffusion du champ de pression incidente par le corps
44
Chapitre I
occupant le volume 𝛺𝑠 . Cette équation est valable pour tous les points, qui n'appartiennent pas à
∂Ωs . Pour 𝑥 ∈ ∂Ωs, l’équation (I.67) devient :
G( x, y ) G0 ( x, y ) p( y )
p ( x)
f
p( y ) n
n
p( x) dS
f
G( x, y )
n
dS (I.68)
avec:
1
G0 ( x, y) I.69
4 r
Les équations I.67 et I.68 mettent en relation les valeurs de la pression 𝑝(𝑥) en un point 𝑥
avec la pression 𝑝(𝑦) sur la surface et sa dérivée normale 𝜕𝑝(𝑦)⁄𝜕𝑛 sur la surface. Cependant,
ces valeurs de surface doivent être connues, les conditions aux limites de la surface 𝜕𝛺
fournissent soit la pression, soit sa dérivée normale, soit une combinaison des deux (Annexe A).
Une configuration importante qui se produit dans les problèmes vibro-acoustiques, est le cas
d'une structure encastrée d'un baffle infini rigide, l’importance de ce type de cette configuration
se réside dans l’aptitude d’examiner nombreuses propriétés telles que la perte de transmission du
son à travers la structure et l’interaction d'un champ sonore avec un objet au-dessus d'un plan
rigide. Considérons le cas de la formulation directe pour le problème extérieur en tenant en
compte le baffle rigide.
45
Chapitre I
avec:
eikr eikr '
Gb ( x, y) I.71
4 r 4 r '
Les indicateurs vibro-acoustiques, peuvent être calculés, une fois que les vecteurs 𝑈𝑛𝑚 et 𝑃𝑛𝑚
sont retrouvées. Nous présentons ci-dessous les principaux indicateurs qui peuvent caractériser
les radiations acoustiques induisent par une structure vibrante.
La puissance injectée par un champ de pression acoustique 𝑝 appliqué sur la structure, est donnée
par :
1
1
in pin n .vdS i Pin,n Cup U n dS
2
f
2
T *
I.72
Ou 𝑣 ∗ est la conjugué du champ de vitesse.
La puissance sonore rayonnée par une structure vibrante qui génère une pression sonore
rayonnée 𝑝 est donnée par:
1
1
1
rad pn .vdS i Pn Cup U n dS 2 U n Z rad U n dS
2
f
2
T *
2
*
I.73
où [𝑍𝑟𝑎𝑑 ] est la matrice d'impédance de rayonnement de la structure.
46
Chapitre I
L'efficacité de rayonnement est une mesure de la capacité de la structure à émettre du son. Elle
est définie comme :
rad
rad I.74
0C0 SVn2
avec 𝜌0 𝑐0 𝑆𝑉𝑛 correspond à la puissance acoustique rayonnée par la structure de vitesse moyenne
quadratique normale 𝑉𝑛2 .
La perte de transmission acoustique d'une structure désignée par STL est un indicateur de sa
capacité d'isolation acoustique lorsqu'elle est soumise à une excitation acoustique. Elle est définie
comme :
1
STL 10 log I.75
𝜋
avec 𝜏 = 𝜋 𝑡 est le facteur de transmission et 𝛱𝑡 est la puissance sonore transmise équivalente à
𝑖𝑛
la puissance de radiation 𝛱𝑟𝑎𝑑 donnée par l’Eq.(I.73). 𝛱𝑖𝑛 est la puissance sonore incidente. On
suppose qu'une onde sonore plane uniforme variant harmonieusement dans le temps frappe
obliquement le fond du panneau sandwich avec un angle d'élévation incident 𝜙𝑖 et un angle
d'azimut 𝜃𝑖 .
Ai S cos(i )
2
in (i , i ) I.76
2 0C0
1.8. Positionnement du sujet
Le développement d'un système de contrôle des vibrations, fiable et robuste est devenu une
nécessité face aux exigences fonctionnelles des structures. Cependant, le contrôle passif atténue
très efficacement les vibrations de la structure et les bruits acoustiques. La simulation numérique
de structures amorties par des matériaux viscoélastiques est une étape importante dans la
conception permettant de prédire l'efficacité de ce traitement d'amortissement. Cependant, cette
étape nécessite une bonne connaissance du comportement des matériaux viscoélastiques et la
prise en compte de tous les facteurs significatifs influant sur leurs lois de comportement.
L’étude bibliographique a montré que très peu de travaux qui se sont intéressés à
l’application et l’évaluation de comportement des structures sandwich à cœur viscoélastique dans
l’atténuation des vibrations provoquées par chargements mobiles ou par bruits acoustiques.
47
Chapitre I
Cependant, les systèmes de contrôle à amortissement passif par matériaux viscoélastiques ont
montré leur efficacité d’amortissement et d’atténuation des amplitudes de vibrations dans
nombreuses applications aux différents secteurs [97]-[106]. Cela nous a motivé à penser
d’examiner et évaluer le comportement de structures sandwich à cœur viscoélastique en premier
lieu les réponses vibratoires sous l’effet de l’action des chargements mobiles et en deuxième lieu
les réponses vibroacoustic.
Ce présent travail est divisé en trois grandes parties dans le but d’évaluer les propriétés
amortissantes et d’examiner les réponses vibratoires structurelles et acoustiques des structures à
contrôle passif d’amortissement.
La première partie est consacrée à l’étude des propriétés amortissantes des structures par
matériaux viscoélastiques, ainsi, nous y examinons leurs réponses fréquentielles dans le but de
construire une base modale servant à la détermination des réponses de vibrations structurelles et
acoustiques.
La deuxième partie, nous réalisons une évaluation des vibrations sous chargement mobile
des structures par une configuration passive à laquelle un modèle numérique décrivant
l’évolution de la charge mobile, appliquée à la structure, est intégré à la méthode des éléments
finis dont les détails sont reportés au chapitre 3. Le schéma d’intégration directe de Newmark
couplé à une méthode itérative sera utilisé pour résoudre l’équation caractéristique de
mouvement vibratoire des structures sandwich.
Ensuite, les réponses acoustiques sont examinées afin de tester la validité de traitement passif
d’amortissement face aux chargements acoustiques appliqués à la structure en étudiant le
comportement vibro-acoustique des structures avec traitement PCLD « Passive Constrained
Layer Damping » dans un environnement thermique avec la prise en considération de la
dépendance en fréquence et en température des différentes lois de comportement viscoélastiques.
Ainsi, en se basant sur l’approche numérique basée sur la méthode des éléments finis dans la
discrétisation et la formulation de l’équation de mouvement et la méthode asymptotique
numérique dans le but d’approximer la solution du problème complexe aux valeurs propres et
construire la base modale.
48
Chapitre II
Sommaire
Introduction ......................................................................................................................................... 51
Partie A. Poutres sandwichs ................................................................................................................ 52
II.A.1. Formulation mathématique ....................................................................................................... 53
II.A.2. Formulation par éléments finis ................................................................................................. 57
II.A.3. Méthode asymptotique numérique et problème aux valeurs propres ........................................ 60
II.A.4. Validation de l’approche numérique ........................................................................................ 64
II.A.4.1. Vibrations libres ................................................................................................................. 64
II.A.4.2. Réponses fréquentielles ...................................................................................................... 66
Partie B. Plaques Sandwichs ................................................................................................................ 68
II.B.1. Formulation mathématique ....................................................................................................... 69
II.B.2 Formulation par les éléments finis ............................................................................................. 74
II.B.3. Validation de l’approche numérique ......................................................................................... 79
II.B.3.1. Vibrations libres ................................................................................................................. 79
II.B.3.2. Réponses fréquentielles ...................................................................................................... 86
Conclusions ......................................................................................................................................... 88
49
Chapitre II
Introduction
Dans ce chapitre, nous présentons une approche numérique d’analyse des vibrations
linéaires, libres et forcées dans le domaine fréquentiel des structures sandwich munies de
matériaux viscoélastiques sous l’effet de la température avec la prise en compte de l’asymétrie
géométrique. De ce fait, une nouvelle théorie d’ordre élevé a été adoptée. Ainsi, l’équation
de mouvement gouvernant le comportement vibratoire libre et forcé est obtenue en combinant
la formulation variationnelle du principe de Hamilton avec le principe de déplacement virtuel
et en utilisant la configuration du modèle zigzag proposé par Rao [11] dont la théorie d’Euler-
Bernoulli est appliquée aux faces du sandwich et la théorie de Timoshenko au cœur
viscoélastique. Cependant, nous appliquons la méthode des éléments finis pour la
discrétisation de l’équation de mouvement et obtenir le système matriciel décrivant le
comportement vibratoire libre et forcé dans le domaine fréquentiel.
En outre, la difficulté de résoudre le problème aux valeurs propres en raison de la
dépendance en fréquence de la matrice de rigidité nous ramène à utiliser la méthode
asymptotique numérique associée à la procédure de différenciation automatique.
L’approche numérique présentée dans notre travail est considérée comme une approche
typique bien destinée à résoudre le problème des vibrations, quelle que soit la nature de la
structure étudie poutre, plaque ou même coque. Nous nous limitons dans ce chapitre comme
déjà mentionné à présenter les vibrations fréquentielles de deux types de structures « poutre
et plaque » sandwich à cœur viscoélastiques. L’approche numérique est validée dans ce
chapitre en comparant nos résultats aux résultats fournis par d’autres références.
51
Chapitre II. Partie A
Les vibrations linéaires des poutres viscoélastiques sont menées dans le cadre de petites
déformations. Les hypothèses considérées par Bilasse [26] sont modifiées afin de prendre en
considération l’effet des inerties longitudinales et rotationnelles ainsi que l’asymétrique du
sandwich.
Les sections planes transverses initialement normales au plan moyen demeurent planes
et normales au plan moyen durant la flexion pour les couches de face.
La contrainte normale transverse dans la direction (Oz) est très faible par rapport à la
contrainte normale axiale dans la direction (Ox).
Les effets d’inertie axiale, rotationnelle et transversale sont considérés pour la poutre
sandwich.
L’effet du cisaillement est considéré seulement pour la couche centrale viscoélastique.
La continuité du champ de déplacement entre les couches de sandwich est parfaitement
considérée en appliquant le modèle zigzag [11].
La configuration et la déformation géométrique de la poutre sandwich à cœur viscoélastique
asymétrique considérée dans ce travail, sont présentées par la figure II.A.1 et la figure II.A.2,
respectivement, où ℎ1 , ℎ2 et ℎ3 sont l’épaisseur de la couche supérieure, la couche centrale
et la couche inférieure, respectivement.
53
Chapitre II. Partie A
Les champs de déplacements et de déformations pour les couches de face sont donnés
par :
ui x, z, t u0i x, t z zi w, x
wi x, z, t w x, t i 1,3 II.1
ni x, z, t u0i , x z zi w, xx
𝑢, 𝑤 et 𝜀𝑛 sont respectivement le déplacement longitudinal, transversal et la déformation
normale, avec 𝑧1 = ℎ1 ⁄2, 𝑧2 = 0 et 𝑧3 = −ℎ3 ⁄2 correspondent aux épaisseurs des couches,
supérieure, centrale et inférieure, respectivement. 𝑢0𝑖 , représente le déplacement longitudinal
du plan moyen du n-ième couche du sandwich.
u 2 x , z , t u0 x , t z
w2 x, z , t w x, t
(II.2)
n 2 x, z, t u0, x z
c 2 x, z, t x, t w, x
54
Chapitre II
h1 w h h w
u01 usup u0 2 1
x 2 x
2 2
II.3
h w h h w
u03 uinf 3 u0 2 3
2 x 2 2 x
où 𝑢𝑠𝑢𝑝 et 𝑢𝑖𝑛𝑓 sont les déplacements associés aux faces, supérieure et inférieure de la
couche viscoélastique. 𝑢0 = 𝑢02 correspond au déplacement du plan moyen de la couche
viscoélastique.
U K W II.4
où 𝛱 est l'énergie potentielle totale, 𝑈,𝐾 et 𝑊 sont l'énergie de déformation, l’énergie
cinétique et l'énergie des forces extérieures, respectivement, données par :
3
1
U i i dV II.5.a
i 1
V 2
3
1
K iuiui i wi wi dV II.5.b
i 1 2 V
W fi v wi dV fi s wi dS
V s
II.5.c
avec 𝜎𝑖 , 𝑓𝑖 𝑣 , 𝑓𝑖 𝑠 and 𝜌𝑖 sont respectivement le tenseur de contrainte, la force volumique, la
force de surface et la densité de la n-ième couche du sandwich. Selon la formulation
variationnelle, par l'application du principe de déplacement virtuel, Eq.(II.4) peut être écrit
sous la forme :
T T
T L
N1 u1, x N 2 u2, x N3 u3, x M 1 w, xx M 2, x M 3 w, xx T ( w, x )
dxdt 0 (II.6)
0 0 1 S1u1 u1 3 S3u3 u3 2 S 2 u2 u2 ( 1 S1 3 S3 2 S 2 ) w w P ( x , t ) w
55
Chapitre II. Partie A
Ni E j nj ds Ei Si ui , x
si
avec i 1, 3 II.7.a
M i z z j E j nj ds Ei I i w, xx avec i 1,3 II.7.b
si
N2 x, t E2 n 2 ds II.8.a
s2
T x, t c 2 ds II.8.c
s2
L f t * g t L f t L g t f s g s (II.10.a)
L f t sf s II.10.b
on obtient :
N 2 x, t S 2 E2* * u, x II.11.a
M 2 x, t I 2 E2* * w, xx II.11.b
T x, t
S2
2 1 2
E2* * w, x II.11.c
56
Chapitre II
En substituant la relation de continuité entre la couche centrale et les couches de face Eq.(II.3)
avec les termes de forces Eqs. (II.7 ; II.11; II.12) dans l’équation de mouvement Eq.(II.6), on
obtient :
h2 h1 h2 h1
E1S1 u0 w u0, x , x w, xx E2 ( ) S 2u0, x u0, x
*
2 2 2 2
h h h h
E3 S3 u0 2 1 w u0, x 2 , x 1 w, xx E1 I1w, xx w, xx
2 2 2 2
T L
S
0 0 E2 ( ) I 2 , x, x E2 I 2 w, xx w, xx 2(1 22 ) E2 ( )( w, x )( w, x ) dxdt 0
* *
( 1S1 3 S3 2 S 2 ) w w 1S1 u0 2 1 w
h h h2 h1
0, x
u ,x w, xx
2 2 2 2
3 S3 u0 h2 h1 w
0, x
u
h2
h1
w
, xx S u u P ( x , t ) w
,x 2 2 0 0
2 2 2 2
(II.13)
𝒘𝟏 𝜷𝟏 , 𝜽𝟏 𝜷𝟐 , 𝜽𝟐 𝒘𝟐
𝒖𝟏 𝒖𝟐
57
Chapitre II. Partie A
U i
w N wU e
i
U N uU e Avec U e i 1, 2 (II.14)
N U e wi
i
N w 0 0 1 2 x 3 x
3 2
x 2x x L
3 2
e 0 0 2 x 3x x x L
3 2 3 2
e
Nu 1- x 0 0 0 x 0 0 0
(II.15)
N 0 1- x 0 0 0 x 0 0
𝑥
avec x̅ = 𝐿 , le système matriciel élémentaire décrivant le comportement vibratoire de la
𝑒
poutre sandwich peut être obtenu en remplaçant les équations Eq.(II.14) dans Eq.(13), on
obtient :
M U K U e F
e e e
(II.16)
L
1S1 3 S3 0 N u N u
h2 2
N N dx h2 1S1 3 S3 0 N N u dx
e e
T T L T
4
h2
3 S3 h3 1 S1h1 0 N w, x N u dx 1S1h1 2 S 2 h2 0 N w, x N dx
e e
L T L T
M
e
2
1
4 1 S1h1 3 S3 h3 0 N w, x N w, x dx 2 S 2 0 N u N u dx
e e
L T L T
2 2
E S E S E S L N T N dx
e
2 2 0 w, xx w, xx
1 1 3 3
(II.17)
58
Chapitre II
L
e
h2 2
3 3 0 u,x u,x
T T
1 1
E S E S N N N , x N , x dx
4
e e
h2 E1S1 E3 S3 0 N , x N u , x dx E3 S3 h3 E1S1h1 0 N w, xx N u , x dx
L T L T
E * ( ) S L N T N dx h2 E S h E S h L N T N dx
e e
2 0 u,x u ,x 3 3 3 0 w , xx ,x
2 1 1 1
K 2
e
E1 I1 E3 I 3 N w, xx N w, xx dx 1 E1S1h12 E3 S3 h32 N w, xx N w, xx dx
e e
L T L T
0
4 0
E2 ( ) I 2 0 N , x N , x dx
e
L T
*
E2* ( ) S 2 L
e
0 N N 2 N N w, x N w, x N w, x dx
T T T
2 1 c
(II.18)
Le
F P x, t N w x dx
e T
(II.19)
0
Cependant, afin d’étudier les vibrations libres et établir la base modale du sandwich à
cœur viscoélastique, nous devons résoudre le problème de valeurs propres suivants :
K M U 0
2
II.22
59
Chapitre II. Partie A
Afin de résoudre le problème complexe de Vibrations libres Eq. (II.22), les techniques de
la méthode asymptotique numérique y sont appliquées, on commence par la décomposition
du module 𝐸2∗ (𝜔) et de la matrice de rigidité en deux parties sous la forme :
L
e
h2 2
1 1 3 3 0 u , x u , x
T T
E S E S N N N , x N , x dx
4
e e
h2 E1S1 E3 S3 0 N , x N u , x dx E3 S3 h3 E1S1h1 0 N w, xx N u , x dx
L T L T
E S L N T N dx h2 E S h E S h L N T N dx
e e
0 2 0 u,x u ,x 3 3 3 0 w, xx , x
1 1 1
K0
2
e
E1 I1 E3 I 3 N w, xx N w, xx dx 1 E1S1h12 E3 S3 h32 N w, xx N w, xx dx
e e
L T L T
0
4 0
e
E0 I 2 0 N , x N , x dx
L T
E0 S 2 L
e
0 N N 2 N N w, x N w, x N w, x dx
T T T
2 1 c
(II.24)
S L N T N dx I L N T N dx
e e
2 0 u,x u,x 2 0 ,x ,x
Kc
e
(II.25)
S2 L N T N 2 N T N N T N dx
e
2 1 c 0 w, x w, x w, x
60
Chapitre II
K K0 E Kc II.26
L’introduction de la matrice de rigidité globale dans le problème aux valeurs propres
Eq.(II.23) permet d’exprimer le problème sous la forme :
n 1
N
2 0 n II.28.b
n
n 1
N
E E0 En II.28.c
n
n 1
où, 𝑈0 , 𝜆0 et 𝐸0 sont calculés à partir du problème linéaire aux valeurs propres suivant :
K M U 0
0
2
(II.29)
A0 U 0 0 (II.31.a)
A0 K0 max E Kc M (II.31.b)
le premier problème est identiquement satisfait vu que (𝑈0 , 𝜆0 ) est la solution de l’Eq.(II.31).
Problème d’ordre n +1
61
Chapitre II. Partie A
A0 U n Fn (II.32.a)
n 1 n 1 n 1
Fn max Ek K c {U n k } Ek K c {U n k 1} k M {U j k }
k 1 k 0 k 1 (II.32.b)
max En K c {U 0 } n M {U 0 }
Le problème Eq.(II.32) est mal posé puisque le système comporte n équations et n+1
inconnus. Une condition supplémentaire est nécessaire à chaque ordre pour le résoudre.
Cependant, la propriété d’orthogonalité des modes propres, offre une autre solution
supplémentaire permet de résoudre le problème. Cette propriété donnée par :
T
Un U0 0 (II.33)
La matrice 𝐴0 est irréversible, ce qui implique que l’Eq. (II.32) à chaque ordre 𝑛 ≥ 1 ne
possède pas de solution si, et seulement si la condition de solvabilité suivante est satisfaite :
n 1 n 1 n 1
Fn U E U Kc {U nk } Ek U 0 Kc {U nk 1} k U 0 M {U j k }
T T T T
0 max k 0
k 1 k 0 k 1
nous remarquons dans ce cas que nous avons obtenu une équation implicite due à l’apparition
des termes 𝜆𝑛 et 𝐸𝑛 simultanément. Pour résoudre ce type de problème, la formule de
récurrence de Faà di Bruno pour la différenciation de haut degré des fonctions composées.
Cette formule permet de décrire analytiquement le n-ième coefficient de Taylor de
(𝑓 ∘ 𝑔)𝑛 [125]:
dn n
1 m
dx n
f x f x
n
n ! x m
f x0 Βn,m ( x1 , x 2 , x 3 , x n m1 ) (II.35)
m 1
La formule Eq.(II.36) montre que le coefficient de Taylor d’ordre n peut être décomposé
à une équation qui dépend explicitement aux paramètres initiaux :
f n f1 x1 1 xn f1 xn 0 (II.36)
{𝑓1|𝑥1 =1 } et {𝑓1|𝑥𝑛 =0} sont respectivement la dérivée d’ordre n=1 de 𝑓 évaluée pour 𝑥1 = 1 et
la dérivée d’ordre n de 𝑓 évalué pour 𝑥𝑛 = 0, respectivement. Nous appliquons une
décomposition du module complexe, on obtient :
62
Chapitre II
En E1 1 1 n E1 n 0 (II.37)
L’insertion de l’Eq. (II.37) dans Eq. (II.34) permet d’exprimer après développement
le coefficient de Taylor 𝜆𝑛 en fonction des termes aux ordres inférieurs :
n-1 n-1
α j Ek U 0 [K c ] U n k Ek U 0 [K c ] U n k 1
T T
λn k 1 k 0
U 0
T
[M] U 0 α j E1 λ1 1 U 0 [K c ] U 0
T
λk U 0 [M] U j k α j E1 λn 0 U 0 [K c ] U n k 1
n-1
T T
k 1
(II.38)
U 0 [M] U 0 α j E1 λ1 1 U 0 [K c ] U 0
T T
A0 U 0 U 0 Fn
T (II.39)
U 0 0 0
j 1 j max (II.42)
63
Chapitre II. Partie A
Le processus se termine lorsque la valeur de 𝑎𝑗+1 > 1, à cette étape les valeurs propres et
les vecteurs propres sont calculés par :
N
U s (1 j ) nU n (II.43.a)
n 0
N
s (1 j ) n n (II.43.b)
n 0
R
Dans le but de valider notre approche numérique, les résultats obtenus d’une poutre
sandwich à cœur viscoélastique intercalée entre deux couches d’aluminium pour différentes
conditions aux limites, sont comparés avec ceux obtenus par Bilasse [26]. Nous rappelons
que cette dernière étude néglige l’effet des inerties longitudinales et rotationnelles alors que
notre modèle proposé prend en considération l’ensemble des inerties agissant sur le
comportement vibratoire des poutres sandwichs. Le modèle du comportement viscoélastique
est considéré avec un module 𝐸 indépendant de la fréquence donnée par :
E E0 1 iηv (II.45)
64
Chapitre II
Les fréquences naturelles et les facteurs de perte correspondant aux trois premiers modes
de la poutre sandwich en appuis simples sont présentés au tableau II.A.2. Les résultats y sont
obtenus et comparés pour différentes valeurs de facteur de perte viscoélastique ηv=0.1, 0.6,
1.0 et 1.5. Les résultats montrent que lorsque le facteur d’amortissement du cœur
viscoélastique augmente, la pulsation amortie augmente impliquant l’amélioration de
l’amortissement structurel. Cependant, nous pouvons constater l’efficacité de notre nouvelle
formulation qui tient en compte les effets des inerties longitudinales et rotationnelles, par
conséquent, les résultats restent très acceptables vu de leurs coïncidences avec ceux obtenus
par la référence [26].
Tableau II.A. 2. Fréquences naturelles et facteur de perte de la poutre sandwich en appuis simples pour
différentes valeurs de 𝜂𝑐
Présente formulation Billasse [26]
𝜂𝑣 N° mode ω (Hz) 𝜂 ⁄𝜂𝑣 ω (Hz) 𝜂 ⁄𝜂𝑣
1 148.41 0.3501 148.51 0.3502
0.1 2 488.27 0.1953 488.48 0.1958
3 1034.50 0.1067 1034.75 0.1071
1 150.60 0.33283 150,71 0,3328
0.6 2 489.54 0.19388 489,76 0,1943
3 1035.20 0.10654 1035,44 0,1069
1 154.30 0.3053 154.42 0.3052
1.0 2 491.83 0.1914 492.07 0.1918
3 1036.40 0.1062 1036.69 0.1065
1 160.58 0.2629 160,72 0,2626
1.5 2 496.23 0.1866 496,50 0,1870
3 1038.80 0.1055 1039,13 0,1059
Le tableau A.3 rapporte les propriétés amortissantes correspondantes aux trois premiers
modes de vibration d’une poutre Encastrée-Libre pour différentes valeurs du facteur de perte,
où ηv=0.1, 0.6, 1.0 et 1.5. La précision des résultats, est illustrée par l’erreur résiduelle
𝑅(𝑈, 𝜆) = ‖([𝐾 ] − 𝜔2 [𝑀])𝑈‖ qui était inférieur à 0.5 × 10−5 .
65
Chapitre II. Partie A
Tableau II.A.3. Fréquences naturelles et facteur de perte de la poutre sandwich en porte-à-faux pour différentes
valeurs de 𝜂𝑣
Présente formulation Billasse [26]
𝜂𝑣 N° mode ω (Hz) 𝜂 ⁄𝜂𝑣 𝑅 (𝑈, 𝜆) ω (Hz) 𝜂 ⁄𝜂𝑣 𝑅(𝑈, 𝜆)
1 64.37 0.277 1.20e-06 64,10 0,281 3,28e-4
0.1 2 300.59 0.239 2.30e-06 296,70 0,242 1,21e-3
3 754.19 0.151 8.62e-07 744,50 0,154 1,58e-3
1 65.77 0.242 1.61e-06 65.50 0.246 3.28e-4
0.6 2 303.14 0.229 2.84e-06 299.20 0.232 1.21e-3
3 755.94 0.150 8.70e-07 746.30 0.153 1.58e-3
1 67.68 0.199 1.29e-06 67,50 0,202 3,28e-4
1.0 2 307.12 0.214 2.14e-06 303,10 0,218 1,21e-3
3 759.03 0.148 8.55e-07 749,40 0,15 1,58e-3
1 70.12 0.150 1.38e-06 69.80 0.153 3.28e-4
1.5 2 313.31 0.194 2.19e-06 308.80 0.197 1.21e-3
3 764.80 0.144 9.43e-07 754.00 0.146 1.58e-3
II.A.4.2. Réponses fréquentielles
Nous discutons dans cette section les réponses dynamiques des poutres sandwichs du type
PCLD (Passive Constrained Layer Damping) sous une charge ponctuelle harmonique. Pour
cette raison, l’utilisation de la propriété d’orthogonalité des modes propres obtenus par la
résolution de l’Eq.(II.22) permet de réduire l’Eq.(II.21) en système d’équations couplées.
m k () A f
2
(II.46)
Le traitement d'amortissement passif avec une couche contrainte (PCLD) est largement
utilisé dans nombreuses applications d’ingénierie pour objectif d’améliorer les capacités
d'amortissement de la structure en augmentant la contrainte de cisaillement dans la couche
viscoélastique. Cette technique est généralement réalisable en associant une fine couche de
matériaux viscoélastiques intercalée entre deux couches de parement dont la couche
supérieure est la couche contraignante constituant par conséquent un sandwich à fort pouvoir
amortissant. Les propriétés du sandwich sont rapportées par le tableau II.A.4.
Tableau II.A. 4. Propriétés mécaniques et géométriques du sandwich PCLD
Face supérieure Cœur viscoélastique Face inférieure
Module de Young(𝐺𝑃𝑎) 𝐸𝑓 = 49 𝐸𝑐 = 2𝐺(1 + 𝜐𝑐 ) 𝐸𝑓 = 70
Soft G = 0.895 + 1.3067𝑖
Module de cisaillement (𝑀𝑃𝑎) / /
Hard 𝐺 = 9.89 + 14.4394𝑖
Coefficient de Poisson 𝜐𝑓 = 0.3 𝜐𝑐 = 0.3 𝜐𝑓 = 0.3
Masse volumique(𝐾𝑔⁄𝑚3 ) 𝜌𝑓 = 750 𝜌𝑐 = 968.1 𝜌𝑓 = 2110
Epaisseur 𝑚𝑚 ℎ𝑓 = 2 ℎ𝑐 = 1 ℎ𝑓 = 4
Longueur (𝑚) 𝐿 = 0.4
Largeur(𝑚) 𝑙 = 0.03
66
Chapitre II
« soft » et « hard » et pour une poutre en porte-à-faux. Les figures de A.4.a et A.4.b, illustrent
les deux réponses obtenues par l’approche analytique [45] et notre modèle d’éléments finis à
l’extrémité de la poutre.
(a)
(b)
Figure II.A. 4. Comparaison des réponses fréquentielles à l’extrémité de la poutre sandwich en porte-à-faux
entre la présente approche et l’approche analytique de la référence [45] a- hard cœur ; b- soft cœur.
Nous pouvons constater que les résultats des deux approches sont très proches en mépris
des valeurs de fréquences propres obtenues par Cai et al [45] qui sont un peu moins élevées
par rapport à nos résultats. En effet, ceci, prouve clairement l’efficacité de la présente
approche proposée.
67
Chapitre II. Partie B
Dans le but d’étudier le comportement linéaire libre et forcé des plaques sandwich à cœur
viscoélastiques, nous considérons une structure sandwich rectangulaire présentée dans la
figure B.1, constitue des couches de face non homogène et une couche d’un matériau
viscoélastique non vieillissant. Ainsi, le comportement linéaire des structures plaques, est
mené dans le cadre de petites déformations dont la théorie de Kirchoff-love [56] est appliqué
aux faces extérieures de la plaque et la théorie de Mindlin [58] au cœur viscoélastique. Dans
notre approche numérique, la modélisation du comportement des plaques sandwichs
asymétriques, est basée sur les hypothèses suivantes :
La continuité de champ de déplacement entre les trois couches de la plaque, est
parfaitement considérée en admettant un glissement nul des couches de faces extérieures.
La contrainte normale transverse dans la direction (Oz) est très faible par rapport à la
contrainte normale axiale dans la direction (Ox).
Les effets d’inertie axiale, rotationnelle et transversale sont considérés pour la plaque
sandwich.
La prise en compte du cisaillement, est adoptée seulement pour la couche centrale
viscoélastique.
La configuration géométrique de la plaque sandwich à cœur viscoélastique asymétrique
considérée dans ce travail, sont présentées par la figure B.1, respectivement, où ℎ1 ,
ℎ2 , et ℎ3 sont les épaisseurs de la couche supérieure, la couche centrale et la couche
inférieure, respectivement.
69
Chapitre II. Partie B
w
ui x, y, z , t u0i x, y, t z zi
i 1,3 (II.47.a)
x
w
vi x, y, z, t v0i x, y, t z zi i 1,3 (II.47.b)
x
𝑢 et 𝑣 sont respectivement les déplacements longitudinaux dans la direction (Ox) et (Oy),
respectivement, avec 𝑧1 = ℎ1 ⁄2, et 𝑧3 = −ℎ3 ⁄2 correspondent aux épaisseurs des couches,
supérieure et inférieure, respectivement. 𝑢0𝑖 , représente le déplacement longitudinal du plan
moyen de la couche i du sandwich.
Le déplacement transversal 𝑤 est donné par :
wi x, y, z, t w x, y, t (II.48)
Le champ de déformation linéaire, selon la théorie des plaques adoptée pour les couches
de face, est donné par :
u0i 2w
2
x x xi
v0i 2w
mi ; fi z 2 ; ti T yi with i 1,3 (II.49)
x y
u0i v0i 2w xyi
y 2
x
xy
avec 𝜀𝑚 , 𝜀𝑓 et 𝜀𝑡 représentent les champs de déformation en membrane, en flexion et la
déformation due à la variation de la température, respectivement.
En utilisant la loi de Hook, les contraintes dues aux efforts membranes 𝜎𝑚 , aux moments
fléchissant 𝜎𝑓 sont exprimés par :
Exi hi E yi hi vxyi
0 w
1 vxyi v yx ,i 1 vxyi v yxi 2 x
h2 E hv E yi hi w avec i=1,3 (II.50)
fi D fi fi i xi i yxi 0 2
12 1 vxyi v yxi 1 vxyi v yxi y
0 0 Gxyi hi 2 w
xy
Exi hi E yi hi vxyi
0 ui
1 vxyi v yxi 1 vxyi v yxi x
E hv E yi hi vi
mi Dmi mi xi i yx ,i
0 avec i=1,3 (II.51)
1 vxyi v yxi 1 vxy ,i v yxi x
0 0 Gxyi hi ui vi
y x
70
Chapitre II
Les contraintes thermiques appliqué aux couches de face extérieure sont données par :
Exi hi E yi hi vxyi
0
1 vxyi v yxi 1 vxy ,i v yx ,i
E hv E yi hi xi
Nti Dmi ti T 0 yi avec i=1,3 (II.52)
xi i yxi
où 𝛼𝑥𝑖 , 𝛼𝑦𝑖 et 𝛼𝑥𝑦𝑖 sont les coefficients d’expansion thermique pour la n-ième couche de face
de la plaque sandwich. 𝛥𝑇 est la variation de la température 𝛥𝑇 = 𝑇 − 𝑇0, avec 𝑇0 est la
température ambiante.
u0 x
x x xi
v0 y
x
m2 ; f 2 z ; ti T yi ; (II.54)
y y y
u0 v0 y xyi
x
y x y
x
Les contraintes due aux efforts membranes 𝜎𝑚 , aux moments fléchissant 𝜎𝑓 ainsi qu’aux
efforts tranchants 𝜎𝑐 dû au cisaillement dans la couche centrale considéré homogène sont
exprimés par :
71
Chapitre II. Partie B
u2
1 v2 0 x
E2* h2 v2
m2 Dm 2 t 2 0
1 v2 2
v2 1 (II.55)
x
1 v2
0 0 u v
2
2
2
y x
x
1 x
v2 0
E2* hi 3
f2 D f 2 fi
y
v2 1 0 (II.56)
12 1 v2 2 y
1 v2
0 0 y x
2
y
x
w
x
E hi 1 0
*
x (II.57)
c Dc 2
2 1 v2 0 1 w
y
y
La contrainte thermique pour la couche centrale est donnée par :
1 v2 0 x2
E2* h2
Nt 2 Dm 2 t 2 T v2 1 0 y2 (II.58)
1 v2
2
1 v2 xy 2
0 0
2
où 𝛼𝑥𝑖 , 𝛼𝑦𝑖 et 𝛼𝑥𝑦𝑖 sont les coefficients d’expansion thermique pour la n-ième couche de face.
ΔT est la variation de la température 𝛥𝑇 = 𝑇 − 𝑇0, avec 𝑇0 est la température ambiante.
c. Equation de continuité
Les équations décrivant la continuité des champs de déplacement entre la couche centrale
et les couches de faces sont données par :
h2 h w
u01 u0 1 II.59.a
2 2 x
h2 h1 w
v01 v0 II.59.b
2 2 x
h2 h3 w
u03 u0 II.59.c
2 2 x
h h w
v03 v0 2 3 II.59.d
2 2 x
72
Chapitre II
u0 x 2w
2
x x x
v0 h2 y
h1
2w
m1 2 (II.60.a)
x 2 y 2 y
u0 v0 y 2w
y x x 2
y
x
xy
u0 x w
x 2
x x
v0 h2 y h3 w
m3 2 (II.60.b)
x 2 y 2 y
u0 v0 y w
y x x 2
y x xy
d. Equation de mouvement
W fi v wi dV fi s wi dS fi t wi dS (II.62.c)
V s s
fi t Tf 2 Dm 2 t 2 (II.63)
73
Chapitre II. Partie B
Πdt U K W dt
T T
0 0
(II.64)
Notre modèle numérique utilisé dans ce présent travail repose sur la méthode des éléments
finis dans l’ordre de formuler en premier lieu le problème aux valeurs propres et de formuler
en deuxième lieu le système d’équations décrivant le comportement vibratoire forcé sous
chargement harmonique ponctuel. L’élément considéré est un élément quadrangulaire
linéaire (figure II.B.2) à quatre nœuds dont le nombre de degrés de liberté par nœud est quatre
(7) qui sont les déplacements longitudinaux 𝑢 et 𝑣, le déplacement transversal 𝑤, les rotations
de la normale de la couche centrale 𝛽𝑥 et 𝛽𝑦 et les rotations 𝜃𝑥 = 𝑑𝑤⁄𝑑𝑦 et 𝜃𝑦 = 𝑑𝑤⁄𝑑𝑥 .
Node 4 Node 3
y
O x
Node 1 Node 2
Figure II.B. 2. Elément quadrangulaire linéaire à quatre nœuds
74
Chapitre II
x 20 21 x 22 y 23 xy (II.65.d)
y 24 25 x 26 y 27 xy (II.65.e)
T
w w
U w u0 v0 x y P (II.66)
x y
avec {𝛼 } est le vecteur de coefficients d’interpolations et [𝑃] est la matrice d’interpolations
définie comme suit :
1 𝑥 𝑦 𝑥 2 𝑥𝑦 𝑦 2 𝑥 3 𝑥 2 𝑦 𝑥𝑦 2 𝑦 3 𝑥 3 𝑦 𝑥𝑦 3 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 1 0 2𝑥 𝑦 0 3𝑥 2 2𝑥𝑦 𝑦 2 0 3𝑥 2 𝑦 𝑦3 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
[𝑃] = 0 0 0 0 0 0 0 0 0 0 0 0 1 𝑥 𝑦 𝑥𝑦 0 0 0 0 0 0 0 0 0 0 0 0
000 0 0 0 0 0 0 0 0 0 0 0 0 0 1 𝑥 𝑦 𝑥𝑦 0 0 0 0 0 0 0 0
000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 𝑥 𝑦 𝑥𝑦 0 0 0 0
[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 𝑥 𝑦 𝑥𝑦]
(II.67)
le nombre de degrés de liberté par élément est considéré 28 et le vecteur des déplacements
nodaux élémentaire est donné par :
T
w wi
U wi i
e
u0i v0i xi yi avec i 1, 2,3, 4 (II.68)
x y
par conséquent, le vecteur de déplacements élémentaires devient :
U e C (II.69)
75
Chapitre II. Partie B
C
1
Ue (II.70)
En remplaçant cette dernière équation Eq. (II.70) dans Eq.(69), le vecteur de champ des
déplacements devient :
U P C U e N U e
1
(II.71)
avec [𝑁] = [𝑃][𝐶 ]−1 est la matrice des fonctions d’interpolation 7 × 28 correspond aux
𝜕𝑤 𝜕𝑤
différents champs des déplacements {𝑤 𝑢0 𝑣0 𝛽𝑥 𝛽𝑦 }, donnée par :
𝜕𝑥 𝜕𝑦
T
N Nw N w N w Nu Nv N x Ny (II.72)
x y
L’expression de champ de déplacements en fonction des déplacements nodaux devient
donc :
par conséquent, les champs de déformations peuvent être écrits en fonction des vecteurs
d’interpolation tout en remplaçant le nouveau champ de déplacements Eq.(II.71) dans les
Eqs.(II.49 ; II.54), on obtient :
h2 h
m1 Bm Bf 2 1 Bf (II.74.a)
2 2
h2 h
m3 Bm Bf 2 1 Bf (II.74.b)
2 2
Nu 2 Nw N x
x x
2
x
N
2
N w
N y
Bm v
; Bf 2 ; Bf 2 (II.75.a)
x y x
Nu N v N 2 N x N y
y x 2 w
xy
y x
N w N w
x N x x N
Bc1 ; Bc 2 ; BN u (II.75.a)
N w N N w Nv
y y
y
76
Chapitre II
M U K U e F
e e e
(II.76)
4 4
1 T
B f Dm1h1 Dm 3 h3 B f B f Dm1 Dm 3 B f Bm Dm 2 Bm B f 2 D f 2 B f 2
2 2 T T T T
4
BT D B BT D B BT D B BT D B
c2 c c2 c2 c c1 c1 c c2 c1 c c1
(II.77.a)
h h
11 31 N
h h
B
T
B N
h2 2
4
T
4
Bc 2 Bc 2 1 1 3 3 Bc 2 BN BN Bc 2
T T
h h 1S1h1 3 S3 h32 h2
2 2 2
2 4
S h2 S h2
1 1 1 3 3 3
4
Bc1 Bc1 2 h2 BN BN 1h1h2 3h3h2 Bc1 Bc1
T T T
(II.77.b)
F P x, t N
a b
e T
w dxdy (II.77.c)
a b
M U K U F vt (II.78)
77
Chapitre II. Partie B
où, [𝑀] et [𝐾 ] sont respectivement les matrices de masse de force nodale de rigidité globales,
{𝐹 } est le vecteur nodal des forces extérieures, U est le vecteur de déplacement nodal.
La réponse en amplitudes fréquentielles de la plaque sandwich à cœur viscoélastique, peut
être obtenue en admettant l’expression suivante du vecteur des déplacements nodaux :
U U eit (II.79)
K M U 0
2
(II.80)
E2 E0 E K K0 E Kc
e e e
(II.81)
4 4
1 T
B f Dm1h1 Dm 3 h3 B f B f Dm1 Dm 3 B f Bm Dm 2 Bm B f 2 D f 2 B f 2
2 2 T T T 0 T 0
4
BT D 0 B BT D 0 B BT D 0 B BT D 0 B
c2 c c2 c2 c c1 c1 c c2 c1 c c1
78
Chapitre II
1 v2 1 v2
0 0 0 0
2 2
w
1 v 0
E0 hi 1 0 x 0 E0 h2
x 2
Dc
0
; Nt 2 2 2
v 1 0 (II.83.b)
2 1 v2 0 1 w 1 v2
y 1 v2
y 0 0
2
1 v2 0 1 v2 0
h hi 3
Dmc 2 2 2 v2 1 0 ; D cf 2 v2 1 0
(II.83.c)
1 v2 12 1 v2 2
1 v2 1 v2
0 0 0 0
2 2
w 0
x x
1 v2
h 1 0 h
D cf 2 i
; Ntc2 2 2 v2 1 0 (II.83.c)
2 1 v2 0 1 w 1 v2
y 1 v2
y 0 0
2
La matrice de rigidité globale du comportement, est obtenue après avoir réalisé
l’assemblage des matrices élémentaires :
K K0 E Kc (II.84)
Nous y rappelons que les différentes étapes de la méthode asymptotique numérique sont
appliquées avec le même principe utilisé pour la poutre sandwich avec cœur viscoélastique,
les fréquences et les vecteurs propres sont obtenus avec les expressions Eq.(II.38) et
Eq.(II.39), respectivement.
Les résultats des modes propres pour une plaque sandwich à cœur viscoélastique, obtenus
par la présente approche sont présentés et comparés avec les résultats obtenus par Xiangyang
et Kaiping [126]. Le modèle de comportement du cœur viscoélastique est considéré dans
79
Chapitre II. Partie B
cette section avec un module 𝐸 indépendant de la fréquence donné par l’Eq.(II.45). Les
propriétés mécaniques et géométriques de deux plaques sandwichs sont rapportées dans le
tableau II.B.1.
Nous pouvons constater que les résultats de trois approches se coïncident quoique le
troisième mode de Xiangyang et Kaiping [126] soit un peu différent de celui obtenu par notre
approche et du logiciel commercial Abaqus.
80
Chapitre II
81
Chapitre II. Partie B
Figure II.B. 3. Comparaison des modes propres entre les résultats de la présente approche, de la référence [126] et ceux du logiciel commercial « Abaqus ».
82
Chapitre II
Les résultats de fréquences naturelles des deux approches pour différentes valeurs de
températures Δ𝑇 = 0°, 25°, 50° et 75° sont rapportés dans le tableau II.B.3. Nous observons une
légère différence entre les deux approches. Cela est dû à l’approche analytique utilisée par
Xiangyang et Kaiping [126] qui est basée sur les fonctions trigonométriques pour exprimer les
vecteurs propres, cela peut provoquer des erreurs au niveau de l’estimation des fréquences
naturelles.
Tableau II.B. 2. Comparaison des fréquences naturelles obtenues avec les résultats de la référence [126] et ceux
du logiciel commercial Abaqus.
Δ𝑇 (°) 0 25 50 75 100
𝜂𝑣 N° mode ω (Hz) 𝜔 (𝐻𝑧) ω (Hz) 𝜔 (𝐻𝑧) ω (Hz)
1 228.10 196.26 158.12 107.39 45.14
2 390.01 353.83 313.50 267.39 176.24
3 439.92 391.04 335.11 268.12 210.80
Présente formulation 0.1
4 543.27 491.16 432.81 365.92 281.92
5 568.79 523.73 474.39 419.57 345.71
6 680.29 614.00 539.64 453.70 355.90
1 235.09 204.06 167.65 120.72 32.22
2 397.25 361.73 322.32 272.78 183.70
3 443.16 394.62 339.22 277.37 223.55
Référence [126] 0.1
4 548.44 496.68 438.85 372.14 290.50
5 573.55 528.77 479.84 425.31 347.81
6 681.60 615.37 541.09 454.84 362.67
Ensuite, nous validons les résultats obtenus de la plaque présentée dans le tableau II.B.4 avec
les résultats obtenus par Bilassse et al [26]. Les figures II.B.4 et II.B.5 présentent les modes
propres réels et imaginaires obtenus par les deux approches auxquelles nous y remarquons une
coïncidence entre les résultats. Les fréquences naturelles, sont aussi comparés et donnés dans le
tableau II.B.4. Les résultats sont presque identiques à ceux obtenus par la référence [26]. Cela
valide l’efficacité de notre approche proposée.
83
Chapitre II. Partie B
Figure II.B. 4. Comparaison des modes propres réels entre les résultats de la présente approche et l’approche de
Bilasse et al [26].
84
Chapitre II
Figure II.B. 5. Comparaison des modes propres imaginaires entre la présente approche et l’approche de Bilasse
et al [26]
85
Chapitre II. Partie B
Tableau II.B. 4. Validations des fréquences naturelles avec les résultats de Bilasse et al [26].
Dans le but de valider les réponses forcées, les réponses de déplacement et de vitesse obtenues
par notre approche, sont comparées et présentées dans les figures II.B.6 et II.B.7 avec les résultats
de l’approche analytique développée par Xiangyang et Kaiping [126]. Nous y pouvons remarquer
que nos réponses obtenues correspondent bien les réponses de Xiangyang et Kaiping [126].
86
Chapitre II
Figure II.B. 6. Comparaison des réponses fréquentielles de déplacement entre la présente approche et l’approche
de Xiangyang [126]
3
Figure II.B. 7. Comparaison des réponses fréquentielles de vitesse entre la présente approche et l’approche de
Xiangyang [126]
87
Chapitre II. Conclusions
Conclusions
Dans ce chapitre, nous nous sommes intéressé à l’étude du comportement des poutres et
plaques sandwichs à cœur viscoélastique. Une approche numérique basée sur une théorie d’ordre
élevé en considérant les inerties longitudinales et rotationnelles ainsi que l’asymétrie des poutres
sandwichs à cœur viscoélastique a été présentée afin de caractériser leurs réponses vibratoires.
L’objectif principal était donc de valider notre approche numérique en comparant les résultats
de vibrations libres et forcées avec d’autres recherches établies.
Les résultats de l’étude comparative ont montré l’efficacité de nos approches numériques
unidimensionnelle et bidimensionnelle adoptées pour caractériser le comportement des structures
poutres et plaques sandwichs viscoélastiques.
88
Chapitre III
Sommaire
Introduction ......................................................................................................................................... 90
Partie A. Poutres sandwichs ................................................................................................................. 91
III.A.1. Résultats numériques et analyses ............................................................................................ 92
III.A.1.1. Modèle du comportement viscoélastique à module constant ............................................. 93
III.A.1.1.1. Etude paramétrique ................................................................................................... 93
a. Effet du facteur de perte viscoélastique................................................................................ 93
b. Effet de l’orientation des fibres des couches composites ..................................................... 96
c. Effet du rapport d’épaisseur ................................................................................................. 98
d. Effet de l'asymétrie de la structure sandwich ..................................................................... 100
III.A.1.2. Modèle de comportement à module variable en fréquence.............................................. 102
III.A.1.2.1. Etude comparative ........................................................................................................ 102
III.A.1.2.2. Etude de l’effet d’épaisseur viscoélastique ................................................................... 105
Partie B. Plaques Sandwichs .............................................................................................................. 109
III.B.1. Résultats numériques et analyses .......................................................................................... 110
III.B.1.1. Vibrations des plaques sandwichs PCLD à cœur viscoélastique ..................................... 110
III.B.1.2. Modèle de comportement à module variable ................................................................... 114
III.B.1.2.1. Vibrations des plaques sandwichs PCLD à cœur viscoélastique.............................. 115
III.B.1.2.2. Vibrations des plaques sandwichs à cœur viscoélastique DYAD606 ...................... 117
Conclusions ....................................................................................................................................... 120
89
Chapitre III
Introduction
Ce chapitre porte sur l’étude des vibrations libres et forcées des structures sandwich à
cœur viscoélastique en utilisant l’approche numérique basée sur la méthode des éléments finis
et la méthode asymptotique numérique. Cette approche a été construite en utilisant une
approche d’ordre élevé en tenant en compte des inerties longitudinales et rotationnelles ainsi
que l’asymétrie des structures sandwich.
Dans la première partie, nous examinons tout d’abord les fréquences naturelles et les
facteurs de perte caractérisant les propriétés amortissantes des poutres sandwichs munies de
matériaux viscoélastiques issus du problème aux valeurs propres. Ensuite, nous abordons les
vibrations forcées dans le domaine fréquentiel de plusieurs configurations adoptées aux
poutres sandwichs sous l’effet d’une charge ponctuelle.
91
Chapitre III. Partie A
Dans ce chapitre, nous étudions le comportement dynamique libre et forcé des poutres
sandwichs à cœur viscoélastique en utilisant plusieurs modèles de lois viscoélastiques
dépendante et indépendante de la fréquence et de la température avec différentes
configurations de la poutre sandwich présentée dans la figure III.A.1.
Les réponses dynamiques des poutres sandwichs ont été examinées sous une charge
ponctuelle harmonique de la forme :
respectivement. De ce fait, nous proposons un modèle d’éléments finis avec 100 éléments, le
degré de liberté nodale est quatre (4) qui sont le déplacement longitudinal 𝑢, le déplacement
transversal 𝑤, la rotation 𝑑𝑤⁄𝑑𝑥 et la rotation de la normale de la couche centrale 𝛽. Le
paramètre de précision de l’Eq. (II.41) est 𝜖 = 10−8 .
Différentes configurations des poutres telles que l’orientation de fibres des couches de
face, le rapport d’épaisseur et le facteur de perte viscoélastique seront évaluées afin d’évaluer
leurs effets sur le comportement dynamique.
92
Chapitre III
E2 E0 1 iηv (III.2)
avec E0 est le module de stockage et ηv est le facteur de perte viscoélastique. Ce modèle est
largement utilisé pour étudier le comportement viscoélastique.
L’effet de différents paramètres de configurations des poutres tels que le facteur de perte
viscoélastique 𝜂𝑣 , l’orientation de fibres 𝜃, le rapport d’épaisseur h/H et l’asymétrie de la
poutre sandwich sur les propriétés amortissantes et les réponses vibratoires de la structure,
sont considérés dans cette étude paramétrique. Les propriétés mécaniques et géométriques de
la poutre considérée sont données par le tableau III.A.1.
Les fréquences naturelles et les facteurs de perte de la structure poutre de trois premiers
modes correspondants aux différentes valeurs de facteur de perte viscoélastique pour une
poutre sandwich Appuie-Appuie (A-A) et autre Encastrée-Libre (E-L) avec 𝜃 = 0° et
h/H=0.1 sont rapportés par le tableau III.A.2. On observe l’augmentation du facteur de perte
de la structure lorsque le facteur interne viscoélastique augmente. Les résultats sont tout à fait
logiques en raison de l’importance de la partie imaginaire du module de Young qui traduit la
93
Chapitre III. Partie A
capacité des matériaux viscoélastiques à amortir les vibrations notamment lorsque cette partie
imaginaire est considérable. La figure III.A.2 valide les résultats pour la même poutre
sandwich en appuis simples tout en remarquant la proportionnalité entre le facteur de perte
interne viscoélastique et le facteur de perte de la structure néanmoins les fréquences propres
restent constantes.
Tableau III.A. 2. Fréquences naturelles et facteur de perte de la poutre sandwich pour différentes valeurs de
𝜂𝑣 obtenues pour 𝜃 = 0°etℎ⁄𝐻 = 0.1.
A-A E-L
𝜂𝑐 N° mode ω (Hz) 𝜂 ω (Hz) 𝜂
1 182.17 2.87e-03 67.652 9.13e-03
0.1 2 620.63 9.94e-02 407.09 1.87e-03
3 720.48 9.32e-04 1132.3 5.89e-04
1 182.17 8.62e-03 67.663 2.74e-02
0.3 2 620.85 2.98e-01 407.10 5.61e-03
3 720.39 2.59e-03 1132.30 1.77e-03
1 182.18 1.72e-02 67.70 5.47e-02
0.6 2 621.31 5.96e-01 407.12 1.12e-02
3 720.31 4.78e-03 1132.30 3.53e-03
1 182.19 2.87e-02 67.79 9.07e-02
1.0 2 720.27 7.67e-03 407.17 1.87e-02
3 622.22 9.91e-01 1132.3 5.89e-03
1 182.23 4.31e-02 67.963 1.35e-01
1.5 2 720.27 1.13e-02 407.26 2.80e-02
3 623.91 1.48e+00 1132.4 8.83e-03
Figure III.A. 2. Variations des fréquences naturelles et des facteurs de perte de la poutre sandwich en appuis
simples pour différentes valeurs de 𝜂𝑣 obtenues avec 𝜃 = 0° et ℎ⁄𝐻 = 0.1 ((a) fréquence naturelle ; (b)
facteur de perte)
94
Chapitre III
{0.0,0.1,0.6,1.0,1.5} sont présentés dans les deux figures III.A.3 et III.A.4 correspondants
aux conditions aux limites Appuie-Appuie et porte-à-faux, respectivement. On remarque que
les pics d’amplitudes obtenues pour la valeur de 𝜂𝑣 = 0 sont très grands à ceux obtenus pour
des valeurs plus élevées de 𝜂𝑣 > 0.
Figure III.A. 3. Comparaison des réponses fréquentielles au milieu de la poutre sandwich en appuis simples
de cœur viscoélastique et faces composites obtenues pour différentes valeurs de 𝜂𝑣
Figure III.A. 4. Comparaison des réponses fréquentielles à l’extrémité de la poutre sandwich en porte-à-faux
de cœur viscoélastique et faces composites obtenues pour différentes valeurs de 𝜂𝑣
Pour la poutre en appuis simples, en ce qui concerne le premier pic d’amplitude qui
correspond la première fréquence propre 𝜔 = 182, la zone zoomée des amplitudes, illustre
que les pics d’amplitudes diminuent progressivement avec l’augmentation de 𝜂𝑣 . Les résultats
valident tout ce qui est obtenu auparavant, d’après lesquels les matériaux viscoélastiques sont
capables d’atténuer les amplitudes de vibrations grâce à leurs propriétés amortissantes
95
Chapitre III. Partie A
Les quatre cas étudiés pour 𝜃 = {0. 30. 60 . 90} reportés par le tableau.III.A.3 pour les
deux poutres considérées montrent que les fréquences amorties varient en fonction de
l’orientation des fibres dont on remarque que les fréquences naturelles atteignent des valeurs
élevées pour 𝜃 = 0 et 𝜃 = 60 alors qu’elles atteignent de basses valeurs avec la configuration
𝜃 = 90 et 𝜃 = 30. Les variations des propriétés amortissantes, pulsation amortie et facteur
de perte structurelle correspondant aux trois premiers modes de vibration sont présentés dans
la figure III.A.5 pour une poutre sandwich simplement appuyée à laquelle on observe que les
fréquences propres atteignent des valeurs maximales pour l’orientation 0°. 50° et 60°, et elles
atteignent des valeurs faibles pour le reste d’orientation.
Tableau III.A. 3. Fréquences naturelles et facteur de perte de la poutre sandwich pour différentes valeurs
de 𝜃 obtenues pour 𝜂𝑣 = 0.1 et ℎ⁄𝐻 = 0.1.
A-A E-L
𝜃° N° mode (
ω Hz ) 𝜂 (
ω Hz ) 𝜂
1 182.18 1.72e-02 67.70 5.47e-02
0 2 621.31 5.96e-01 407.12 1.12e-02
3 720.31 4.78e-03 1132.30 3.53e-03
1 59.731 1.29e-01 25.679 1.68e-01
30 2 215.6 4.72e-02 127.33 8.00e-02
3 473.17 2.31e-02 336.13 3.66e-02
1 121.61 3.77e-02 46.933 1.01e-01
60 2 474.2 1.04e-02 270.01 2.42e-02
3 620.38 5.91e-01 744.78 8.06e-03
1 26.891 1.68e-01 62.968 1.20e-01
90 2 134.99 7.41e-02 29.54 4.19e-02
3 358.26 3.25e-02 505.12 2.05e-02
96
Chapitre III
Figure III.A. 5. Variations des fréquences naturelles et des facteurs de perte de la poutre sandwich en appuis
simples pour différentes valeurs de 𝜃 avec 𝜂𝑣 = 0.1 et ℎ⁄𝐻 = 0.1 ((a) fréquence naturelle ; (b) facteur de
perte)
Cependant, les pics d’amplitudes pour les configurations 30° et 90° sont plus élevés par
rapport à ceux obtenus pour 0° et 60° notamment le premier pic. Les mêmes conclusions
s’appliquent pour la poutre en porte-à-faux indiquant que les configurations avec 0° et 60°
des couches de faces sont les meilleures solutions possibles qui peuvent apporter à la structure
les meilleurs niveaux de renforcement. La perte de la rigidité due à l’orientation des fibres
des couches de face, est la principale cause de l’augmentation des amplitudes des réponses
vibratoires, des études montrent que la rigidité d’un composite monocouche réduit lorsque
l’orientation de fibres est différente de 0°ou 60° ce qui démontre l’augmentation des
amplitudes pour les angles 30°ou 90°.
Figure III.A. 6. Comparaison des réponses fréquentielles au milieu de la poutre sandwich en appuis simples
de cœur viscoélastique et faces composites obtenues pour différentes valeurs de 𝜃 avec 𝜂𝑣 = 0.1 et ℎ⁄𝐻 =
0.1.
97
Chapitre III. Partie A
Figure III.A. 7. Comparaison des réponses fréquentielles à l’extrémité de la poutre sandwich en porte-à-faux
de cœur viscoélastique et faces composites obtenues pour différentes valeurs de 𝜃 avec 𝜂𝑣 = 0.1 et ℎ⁄𝐻 = 0.1
98
Chapitre III
Figure III.A. 8. Variations des fréquences naturelles et des facteurs de perte de la poutre sandwich en appuis
simples pour différentes valeurs de ℎ⁄𝐻 avec 𝜂𝑣 = 0.1 et 𝜃 = 0° ((a) fréquence naturelle ; (b) facteur de
perte).
Figure III.A. 9. Comparaison des réponses fréquentielles au milieu de la poutre sandwich en appuis simples
de cœur viscoélastique et faces composites obtenues pour différentes valeurs de ℎ⁄𝐻 avec 𝜂𝑣 = 0.1 et 𝜃 = 0°
99
Chapitre III. Partie A
Figure III.A. 10. Comparaison des réponses fréquentielles à l’extrémité de la poutre sandwich en porte-à-faux
de cœur viscoélastique et faces composites obtenues pour différentes valeurs de ℎ⁄𝐻avec 𝜂𝑣 = 0.1 et 𝜃 = 0°.
100
Chapitre III
Figure III.A. 11. Variations des fréquences naturelles et des facteurs de perte pour les trois premiers modes de
la poutre sandwich en porte-à-faux obtenues pour différentes valeurs du rapport d'épaisseur h3/h1 ((a)
fréquence naturelle ; (b) facteur de perte)
101
Chapitre III. Partie A
Figure III.A. 12. Comparaison des réponses en fréquence de la poutre sandwich en appuis simples obtenues
pour différentes valeurs du rapport d'épaisseur h3/h1.
Figure III.A. 13. Comparaison des réponses en fréquence de la poutre sandwich en porte-à-faux obtenues
pour différentes valeurs du rapport d'épaisseur h3/h1.
Dans cette section, nous examinons les réponses fréquentielles des poutres sandwichs à
cœur viscoélastique obtenues pour différentes lois viscoélastiques dépendant de la fréquence
et pour différentes conditions thermiques.
Les deux premiers modèles viscoélastiques dans notre présent travail sont considérés
dépendants de la fréquence. Deux poutres sandwichs y sont étudiées avec deux cœurs
viscoélastiques différents, le module opérationnel du premier matériau viscoélastique
polyvinyle de Butyral (PVB) [127] considéré à la température 20° est décrit par un modèle
de dérivée fractionnaire Eq. (III.3) et le deuxième du ISD112 considéré à 20 °et 27° [128]
102
Chapitre III
par le modèle de Maxwell généralisée Eq. (III.4). Le tableau 4 présente les propriétés
mécaniques et géométriques de la structure sandwich viscoélastique.
Tableau III.A. 5. Propriétés mécaniques et géométriques du sandwich PCLD à cœur viscoélastique en PVB
𝜐2
Coefficient de Poisson 𝜐1 = 0.28 𝜐2 = 0.4 𝜐3 = 0.28
= 0.29
𝜌2
Masse volumique (𝐾𝑔⁄𝑚3 ) 𝜌1 = 1520 𝜌2 = 999 𝜌3 = 1520
= 1600
ℎ1 = ℎ𝑢 ℎ3 = ℎ𝑏
Epaisseur (m) ℎ2 = ℎ𝑐 = 0.001
= 0.002 = 0.004
Longueur (𝑚) 𝐿 = 0.6
Largeur (𝑚) 𝑏 = 0.02
3
∆𝑗
𝐸2 ∗ (𝜔) = 2𝐺0 (1 + υ2 ) (1 + ∑ ) (III. 4)
ω − iΩ𝑗
𝑗=1
20°𝐶 27°𝐶
j 𝐺0 ∆𝑗 Ω𝑗 𝐺0 ∆𝑗 Ω𝑗
1 2.8164 31.1176 0.746 468.7
2 5.11×104 13.1162 446.4542 5×105 3.265 4742.4
3 45.4655 5502.5318 43.248 71532.5
Les résultats obtenus des fréquences naturelles et des facteurs de perte des poutres
sandwichs à cœur viscoélastique pour les trois modèles de cœur viscoélastique PVB, ISD112
à 27° et ISD112 à 20° avec les conditions, Appuyée-Appuyée (A-A) et Encastrée-Libre (E-
L), sont présentés dans le tableau III.A.7. Les résultats de l’erreur résiduelle 𝑅(𝑈, 𝜆) =
‖([𝐾 ] − 𝜔2 [𝑀])𝑈‖ illustrent l’efficacité de l’approche numérique utilisée dont 𝑅(𝑈, 𝜆) était
généralement inférieur à 0.5 × 10−3 .
103
Chapitre III. Partie A
Tableau III.A. 7. Fréquences naturelles et facteur de perte des poutres sandwichs avec cœur viscoélastique
dépendant de la fréquence (PVB, ISD112 à 27° et ISD112 à 20°)
PVB Core ISD112 Core at 27° ISD112 Core at 20°
ω (Hz) 𝜂 𝑅(𝑈, 𝜆) ω (Hz) 𝜂 𝑅(𝑈, 𝜆) ω (Hz) 𝜂 𝑅(𝑈, 𝜆)
1 68.08 5.35e-3 0.1501e-4 39.66 0.1600 0.1528e-4 37.19 0.2687 0.1777e-4
2 264.32 1.29e-2 0.4555e-5 141.78 0.1354 0.18403e-4 141.13 0.1231 0.2162e-4
3 570.63 2.05e-2 0.8957e-5 308.16 0.1166 0.18170e-4 308.12 0.0824 0.2272e-4
A-A
104
Chapitre III
Figure III.A. 14. Réponses fréquentielles de déplacement de la poutre sandwich PCLD en appuis simples sans
et avec les cœurs viscoélastiques (PVB, ISD112 à 27° et ISD112 à 20°).
Figure III.A. 15. Réponses fréquentielles de vitesse de la poutre sandwich PCLD en appuis simples sans et
avec les cœurs viscoélastiques (PVB, ISD112 à 27° et ISD112 à 20°)
105
Chapitre III. Partie A
Ainsi, l’interprétation des résultats pour la poutre sandwich PCLD à cœur PVB dévoile
que la structure devient à fort pouvoir amortissant lorsque l’épaisseur de la couche centrale
viscoélastique est épais vu que la variation du facteur de perte est proportionnelle au rapport
d’épaisseur tout en remarquant l’importance de ce facteur de perte pour les grandes valeurs
d’épaisseur viscoélastique. En revanche, les résultats obtenus de la poutre sandwich PCLD à
cœur ISD112 montrent que la valeur élevée du facteur de perte est atteinte avec la plus faible
valeur de rapport d’épaisseur h2/h3=1/4. Le facteur décroit progressivement avec
l’augmentation de l’épaisseur de la couche viscoélastique jusqu’au rapport h 2/h3=1 avant qu’il
devienne proportionnel au-delà de cette valeur. Cela permet de dire que l’amortissement de
la structure est considérable lorsque l’épaisseur de la couche viscoélastique est totalement
différente de celle de la couche inférieure pour une poutre en porte-à-faux.
Figure III.A. 16. Variations des fréquences naturelles et des facteurs de perte de la poutre sandwich en appuis
simples ((a) Cœur PVB ; (b) ISD 112 à 27°)
106
Chapitre III
Figure III.A. 17. Réponses fréquentielles de déplacement de la poutre sandwich PCLD à cœur PVB en appuis
simples pour différentes valeurs de rapport d’épaisseur
107
Chapitre III. Partie A
Figure III.A. 18. Réponses fréquentielles de déplacement de la poutre sandwich PCLD à cœur ISD112 27° en
appuis simples pour différentes valeurs de rapport d’épaisseur
108
Chapitre III. Partie B
Dans la présente section, nous examinons le comportement libre et forcé des plaques
sandwichs avec cœur viscoélastique considéré indépendant de la fréquence décrit auparavant par
l’Eq. (III.2). Le tableau III.B.1 présente les propriétés mécaniques et géométriques de la plaque
sandwich.
Tableau III.B. 1. Propriétés mécaniques et géométriques du sandwich à cœur viscoélastique
Face supérieure Cœur viscoélastique Face inférieure
𝐸11 = 130 𝐸11 = 130
Module de Young (GPa) 𝐸22 = 9.5 𝐸0 = 2.67 × 10−3 𝐸22 = 9.5
𝐺12 = 6.0 𝐺12 = 6.0
Coefficient de Poisson 𝜐1 = 0.3 𝜐2 = 0.49 𝜐3 = 0.3
Masse volumique (𝐾𝑔⁄𝑚3 ) 𝜌1 = 1600 𝜌2 = 999 𝜌3 = 1600
Epaisseur (m) ℎ1 = 0.002 ℎ2 = 0.001 ℎ3 = 0.004
−6 𝛼𝑥1 = −0.3 𝛼𝑥3 = −0.3
Coefficient d’expansion thermique (°/𝐶) × 10 𝛼𝑐 = 18
𝛼𝑦1 = 28.1 𝛼𝑦3 = 28.1
Longueur (𝑚) 𝐿 = 0.4
Largeur(𝑚) 𝑏 = 0.36
Les résultats de la variation des fréquences propres et des facteurs de perte de la plaque
sandwich simplement appyée qui correspondent aux trois premiers modes sont présentés pour
différentes valeurs de facteur de perte viscoélastique dans la figure III.B.1. Nous observons que
la fréquence propre varie très peu face à la variation du facteur de perte viscoélastique. Par contre,
le facteur de perte structurel du sandwich varie significativement en fonction du facteur de perte
viscoélastique, d’où l’augmentation de ce dernier provoque l’augmentation de l’amortisement
structurel. Cela signifie que le facteur de perte viscoélastique a un effet consdérable dans
l’atténuation des vibrations structurelles. Cela est très évident dans les figures III.B.2 et III.B.3,
dans lesquelles les différentes réponses de déplacement et de vitesse correspondent aux
110
Chapitre III
Figure III.B. 1. Variation des fréquences propres et des facteurs de pertes en fonction de facteur de perte
viscoélastique pour une plaque sandwich en appuis simples.
Figure III.B. 2. Réponses en déplacement de la plaque sandwich en appuis simples pour différentes valeurs de
facteur de perte viscoélastique.
111
Chapitre III. Partie B
Figure III.B. 3. Réponses en vitesse de la plaque sandwich en appuis simples pour différentes valeurs de facteur de
perte viscoélastique.
112
Chapitre III
Figure III.B. 4. Variation des fréquences propres et de facteurs de pertes en fonction de rapport d’épaisseur
ℎ2 ⁄ℎ3 pour une plaque sandwich en appuis simples
Ces résultats sont plus clairs dans les figures III.B.5 et III.B.6 auxquelles les réponses de
déplacement et de vitesse sont présentées, tout en remarquant l’augmentation des amplitudes
vibratoires quand les valeurs de rapport d’épaisseur augmentent. Cela peut être interprété par
l’augmentation de la contrainte de cisaillement dans la couche viscoélastique provoquant par
conséquent l’augmentation de la capacité dissipative de la couche viscoélastique et par
conséquent la capacité dissipative de toute la structure. Cependant, cette contrainte peut
provoquer une augmentation des amplitudes vibratoires à cause de la rigidité réduite de la
structure.
113
Chapitre III. Partie B
Figure III.B. 5. Réponses en déplacement de la plaque sandwich en appuis simples pour différents rapports
d’épaisseur h2/h3 viscoélastique
Figure III.B. 6. Réponses en vitesse de la plaque sandwich en appuis simples pour différents rapports d’épaisseur
h2/h3
Les deux premiers modèles viscoélastiques dans notre présent travail sont considérés
dépendants de la fréquence. Deux poutres sandwichs y sont étudiées avec deux cœurs
viscoélastiques différents, le module opérationnel du premier matériau viscoélastique polyvinyle
de Butyral (PVB) considéré à la température 20° est décrit par un modèle de dérivée
fractionnaire Eq.(III.3) et le deuxième du ISD112 considéré à 20 °et 27° par le modèle de
114
Chapitre III
Tableau III.B. 2. Propriétés mécaniques et géométriques du sandwich à cœur viscoélastique en PVB et ISD112
Face composite Cœur viscoélastique Face composite
𝐸11 = 130 × 109 𝐸11 = 130 × 109
PVB ISD112
Module de Young (Pa) 𝐸22 = 9.5 × 109 𝐸22 = 9.5 × 109
𝐺12 = 6.0 × 109 Eq.(37) Eq.(38) 𝐺12 = 6.0 × 109
Coefficient de poisson 𝜐1 = 0.3 𝜐2 = 0.49 𝜐3 = 0.3 𝜐1 = 0.3
𝜌3
Masse volumique (𝐾𝑔⁄𝑚3 ) 𝜌1 = 1570 𝜌2 = 999 𝜌1 = 1570
= 1600
Epaisseur (m) ℎ1 = 0.002 ℎ2 = 0.001 ℎ3 = 0.002
𝛼𝑥1 = −0.3 × 10−6
Longueur (𝑚)
𝛼𝑦1 = 28.1 × 10−6
Longueur (𝑚); Largeur (m) 𝐿 = 0.4 ; 𝑏 = 0.4
Nous étudions maintenant le comportement libre et forcé des plaques sandwichs avec
différents cœurs viscoélastiques dépendants non seulement de la fréquence, mais aussi de la
température. Par conséquent, les résultats obtenus des fréquences naturelles et des facteurs de
perte de la plaque sandwich à cœur viscoélastique pour les trois modèles de cœur viscoélastique
PVB, ISD112 à 20° et ISD112 à 27° sont présentés dans le tableau III.B.3.
Tableau III.B. 3. Fréquences naturelles et facteur de perte de plaque sandwichs simplement appuyée avec cœur
viscoélastique dépendant de la fréquence (PVB, ISD112 à 20° et ISD112 à 27°)
PVB ISD112 à 20° ISD112 à 27°
ω (Hz) 𝜂 ω (Hz) 𝜂 ω (Hz) 𝜂
146.0192 0.0089 89.7080 0.3099 53.0320 0.3121
244.4866 0.0078 160.9029 0.3379 91.6162 0.4228
437.7459 0.0079 280.8981 0.3041 151.9834 0.4323
482.4585 0.0205 275.3187 0.1240 150.5588 0.1778
546.0745 0.0198 323.4942 0.1694 176.1514 0.2510
687.2256 0.0178 420.1072 0.2041 223.9344 0.3194
Les différentes réponses de déplacement et de vitesse pour les trois cœurs viscoélastiques
PVB, ISD112 à 27° et ISD112 à 20° sont comparés avec ceux obtenus avec un cœur constitue
d’une mousse d’un matériau isotrope dont les paramètres Ec=106 Pa, υ=0.3, ρ=40 kg/m3 [126]. Les
figures III.B.7 et III.B.8 présentent les réponses de déplacement transversal et de vitesse obtenues
115
Chapitre III. Partie B
au point d’observation (0.2 (m) ,0.18 (m), 0 (m)). Par conséquent, nous pouvons déduire les
mêmes observations déduites auparavant auxquelles l’avantage de l’application des matériaux
viscoélastiques est très significatif dans l’atténuation des amplitudes de vibrations.
2
Figure III.B. 7. Comparaison entre les réponses fréquentielles de déplacement de la plaque sandwich sans et avec
le cœur viscoélastique (PVB, ISD112 20° et ISD112 27°)
Figure III.B. 8. Comparaison entre les réponses fréquentielles de vitesse de la plaque sandwich sans et avec le
cœur viscoélastique (PVB, ISD112 20° et ISD112 27°)
116
Chapitre III
Dans cette section, nous réalisons une étude du comportement vibratoire des plaques
sandwichs PCLD sous l’effet de la température dont le modèle viscoélastique est considéré
dépendant de la fréquence et de la température. Le module opérationnel du matériau
viscoélastique DYAD606 [26] décrit par le modèle de Maxwell généralisé Eq. (38) est considéré
pour différentes températures T=𝟏𝟎°, 𝟐𝟓°, 𝟑𝟎° et 𝟑𝟖°. Les paramètres du module de Young
viscoélastique dépendant fortement de la fréquence, sont donnés par le tableau III.B.4. Les
propriétés mécaniques et géométriques du sandwich sont présentées dans le tableau III.B.2.
5
∗ ∆𝑗
𝐸2 (𝜔) = 2𝐺0 (1 + υ2 ) (1 + ∑ ) (III. 5)
ω − iΩ𝑗
𝑗=1
10°𝐶 25°𝐶 30°𝐶 38°𝐶
𝐺0(Pa) 5.94×106 2.02×106 2.09×106 1.74×106
J ∆𝑗 Ω𝑗 ∆𝑗 Ω𝑗 ∆𝑗 Ω𝑗 ∆𝑗 Ω𝑗
1 5.88 5.85 9.89 58.18 5.40 73.06 1.15 27.02
2 13.66 2345.09 13.14 6.75 14.15 453.34 3.55 213.35
3 8.94 331.70 18.94 403.00 1.43 8.83 11.79 1257.50
4 6.47 50.65 35.06 3097.38 28.33 3406.80 24.41 7585.29
5 34.52 25033.79 165.97 57244.00 128.85 52781.28 113.12 92517.87
Dans le but d’étudier l’effet de la température sur le comportement des plaques sandwichs
PCLD avec cœur viscoélastique en DYAD606, les différentes réponses du sandwich en appuis
simples et en porte-à-faux sont obtenues pour les températures considérées auparavant pour le
cœur viscoélastique T=0°, 25°, 30° et 25°. Afin de faciliter l’interprétation des résultats, les
résultats sont tracés et présentés dans la figure III.B.9. On observe que les fréquences naturelles
diminuent avec l’augmentation de la température, cela est dû aux contraintes de compression
internes induit par l’augmentation de la température élevée, ce qui affaiblit la rigidité structurelle.
En revanche, l’amortissement structurel traduit par le facteur de perte augmente avec
l’augmentation de la température, ce qui permet d’augmenter la capacité dissipative de la
structure. Généralement, les matériaux viscoélastiques sont très sensibles aux variations de la
température auxquelles ces matériaux deviennent plus souples et plus amortissants pour les
hautes températures.
117
Chapitre III. Partie B
Figure III.B. 9. Variation des fréquences propres et des facteurs de perte pour différentes températures du cœur
viscoélastique DYAD606
Les résultats obtenus des réponses de déplacement et de vitesse des plaques sandwichs sont
présentés dans les figures III.B.10 et III.B.11, respectivement. Ces résultats mettent en évidence
la diminution des amplitudes vibratoires quand la température du cœur viscoélastique augmente.
Figure III.B. 10. Réponses fréquentielles de déplacement de la poutre sandwich à cœur viscoélastique DYAD606
considéré pour différentes températures.
118
Chapitre III
Figure III.B. 11. Réponses fréquentielles de vitesse de la poutre sandwich à cœur viscoélastique DYAD606
considéré pour différentes températures
119
Chapitre III. Conclusions
Conclusions
Dans ce chapitre, l’approche numérique présentée au deuxième chapitre, a été appliquée pour
étudier le comportement des structures sandwich munies de matériaux viscoélastiques. Deux
types de structure sont étudiés dont la première concerne les poutres sandwichs et la deuxième
concerne les plaques sandwichs. Nous présentons les principaux résultats et conclusions tirées
de ce chapitre.
Poutre sandwich
L’approche numérique basée sur la théorie d’ordre élevé en considérant les inerties
longitudinales et rotationnelles ainsi que l’asymétrie des poutres sandwichs à cœur
viscoélastique a été en conséquence présentée afin de caractériser les réponses vibratoires des
poutres sandwichs à cœur viscoélastique. L’objectif principal était donc d’examiner
l’efficacité de l’amortissement des matériaux viscoélastiques décrits par des lois de
comportement dépendant et indépendant de la fréquence et de la température. Ainsi,
différents modèles viscoélastiques, ont été examinés en tenant en compte leurs propriétés
mécaniques spécifiques ainsi que différentes configurations de la structure sandwich
particulièrement celles avec traitement PCLD.
La première partie était consacrée à l’étude des vibrations libres et forcées des poutres
sandwichs à cœur viscoélastique indépendant de la fréquence. Une étude paramétrique y a
été effectuée dans le but de chercher la configuration optimale de la structure qui assure un
niveau d’amortissement élevé, par conséquent, les conclusions suivantes peuvent être tirées:
- Les fréquences naturelles atteignent des valeurs élevées pour 𝜃 = 0° et 𝜃 = 60° alors
qu’elles atteignent de basses valeurs pour la fréquence avec la configuration 𝜃 = 90°
120
Chapitre III
et 𝜃 = 30°. Cependant, les pics d’amplitudes des réponses fréquentielles pour les
configurations 30° et 90° sont plus élevés par rapport à ceux obtenus pour 0° et 60°
notamment pour les premiers pics obtenus.
La deuxième partie abordait les vibrations libres et forcées des poutres sandwichs à cœur
viscoélastique dépendant non seulement du temps, mais aussi de la température. De ce
fait, la méthode des éléments finis combinée à la méthode asymptotique numérique a été
utilisée dans ce présent travail afin de déterminer en premier lieu les propriétés
amortissantes des sandwichs viscoélastiques en résolvant le problème complexe aux
valeurs propres et d’obtenir en deuxième lieu les réponses structurelles et acoustiques
fréquentielles correspondantes aux vibrations forcées sous une charge ponctuelle à
amplitude unitaire. D’après les résultats obtenus, on peut conclure :
- L’étude comparative des résultats des sandwichs sans et avec les trois cœurs
viscoélastiques PVB, ISD112 à 27° et ISD112 à 20° a montré que l’amortissement
apporté par les matériaux viscoélastiques étudiés est significatif en ce qui concerne
l’atténuation des pics d’amplitudes particulièrement pour les hautes fréquences à
l'inverse aux réponses obtenues sans couche viscoélastique. Ainsi, l’étude a illustré
que les cœurs ISD112 sont plus efficaces en matière d’amortissement par rapport au
cœur PVB après avoir remarqué la diminution des amplitudes des réponses
structurelles et acoustiques.
121
Chapitre III. Conclusions
Plaque sandwich
Le même principe d’analyse des poutres sandwichs a été suivi dont la première partie
était consacrée à l’étude de comportement des plaques sandwichs à cœur viscoélastique
indépendant de la fréquence. De ce fait, une étude paramétrique y a été effectuée dans le
but de faire sortir la configuration optimale de la structure plaque. D’après les résultats
obtenus, les conclusions suivantes peuvent être tirées :
122
Chapitre III
D’un autre côté, les vibrations libres et forcées des plaques sandwichs à cœur
viscoélastique dépendant du temps et de la température dont les principales conclusions
sont presque identiques à celles obtenues pour les poutres sandwichs :
- L’atténuation des d’amplitude pour les sandwichs avec les couches viscoélastiques
sont très significative en termes d’atténuations des amplitudes vibratoires par rapport
à ceux sans amortissement viscoélastique.
- Les cœurs viscoélastiques en ISD112 ont plus de capacités dissipatives par rapport
au cœur en PVB.
- Les amplitudes de réponses forées des plaques sandwichs avec cœur en DYAD606
qui fortement dépendant de la fréquence et de la température, diminuent lorsque la
température de ce dernier augmente. Cela signifie que le cœur viscoélastique a plus
de capacité dissipative lorsqu’il devient plus souple.
Ces conclusions mettent en évidence la perte de rigidité due à l’orientation des fibres des
couches de face des poutres sandwichs qui est la principale cause de l’augmentation des
amplitudes des réponses vibratoires. Ainsi, les résultats reflètent le fort pouvoir amortissant
de la structure lorsque l’épaisseur de la couche centrale viscoélastique devient plus mince
pour les deux types de structure. Toutefois, la structure est plus efficace et plus résistante à la
charge dynamique lorsque l'épaisseur de la couche inférieure, est différente de celle de la
couche supérieure.
123
Chapitre IV
Sommaire
Introduction ....................................................................................................................................... 125
Partie A. Poutres sandwichs ............................................................................................................... 126
IV.A.1. Equation du mouvement ....................................................................................................... 127
IV.A.2. Résultats et discussions ......................................................................................................... 130
IV.A.2.1. Poutre sandwich à cœur viscoélastique en PVB ..................................................... 130
IV.A.2.2. Poutre avec traitement PCLD ................................................................................. 132
Partie B. Plaques Sandwichs .............................................................................................................. 135
IV.B.1. Equation du mouvement ....................................................................................................... 136
IV.B.1.1. Modèle viscoélastique dépendant de la fréquence ........................................................... 137
IV.B.1.2. Modèle viscoélastique dépendant du temps..................................................................... 138
IV.B.2. Charge mobile ....................................................................................................................... 140
IV.B.3. Résultats et discussions ......................................................................................................... 141
IV.B.3.1. Plaques sandwichs à cœur viscoélastique dépendant de la fréquence .............................. 142
IV.B.3.2. Plaques sandwichs à cœur viscoélastique dépendant du temps ....................................... 146
Conclusions ....................................................................................................................................... 148
Chapitre IV
Introduction
Dans ce présent chapitre, nous adoptons la même formulation mathématique présentée aux
chapitres précédents pour étudier le comportement vibratoire des structures sandwich à cœur
viscoélastique. Le passage du modèle complexe fréquentiel au modèle réel défini dans le
domaine temporel est réalisé par la définition de la matrice d’amortissement en tant que partie
imaginaire de la matrice de rigidité en utilisant la propriété d’orthogonalité des modes propres.
Cela permet de déterminer le coefficient d’amortissement et le facteur de perte.
Nous étudions tout d’abord les vibrations sous l’action de la charge mobile des poutres
sandwichs. Ensuite, nous abordons le comportement dynamique des plaques sandwich.
Différentes configurations sont étudiées avec différents modèles viscoélastiques dépendants du
temps et de la fréquence.
125
Chapitre IV. Partie A
126
Chapitre IV
Nous y rappelons que nous avons utilisé la même formulation mathématique appliquée dans
le chapitre précédent.
M U K U F vt IV.1
Ou [𝑀] et [𝐾] sont respectivement les matrices de masse et de rigidité globales, {𝐹 } est le
vecteur de force nodale, U est le vecteur de déplacement nodal.
M U [ K ]R i K U F vt
I
IV.2
En utilisant la propriété 𝑈̇ = 𝑖𝜔𝑈, l’équation (IV.2) devient :
M U C U K U F vt
R
IV.3
[𝐾]𝐼
où [𝐶 ] = 𝜔
est la matrice d’amortissement équivalent [128]-[129]. Cependant, la poutre
sandwich simplement appuyée est soumise à une charge mobile à vitesse constante comme
montre la figure IV.A.1, la charge mobile dynamique apparaissant dans l'équation de mouvement
est définie par :
P x, t P0 x vt IV.4
où δ désigne la fonction Delta de Dirac et P0 désigne la force concentrée de magnitude constante,
par conséquent, le vecteur de force nodal {F} devient:
127
Chapitre IV. Partie A
F 0 . . . fi . . 0 0
T
i 1, 2,.....n IV.5.a
P0 N w vt pour t
avec fi P0 x vt N w vt dx IV.5.b
L
0
0 pour t
avec 𝜏 = 𝐿⁄𝑣 est le temps nécessaire pour que la charge traverses complètement la poutre, 𝑣 est
la vitesse constante de la charge mobile et 𝑛 est le nombre des nœuds.
128
Chapitre IV
Début
Accélération initiale
𝑢̈ 0 = 𝑀−1 ( 𝐹0 − 𝐶𝑢0̇ − 𝐾𝑢0 )
Incrémentation du temps
𝑡𝑛+1 = 𝑡𝑛 + Δt
Oui
If 𝑡𝑛+1 < 𝑡𝑚𝑎𝑥
Non
Estimation initiale
∗
Solution finale
𝑢̇ 𝑛+1 = 𝑢̇ 𝑛 + (1 − 𝛾)𝑢̈ 𝑛
∗
1 𝑈𝑠 = 𝑈𝑛
𝑢𝑛+1 = 𝑢𝑛 + Δt𝑢̇ 𝑛 + ( − 𝛽)Δt 2𝑢̈ 𝑛
2 𝑈̇𝑠 = 𝑈̇𝑛
𝑈̈𝑠 = 𝑈̈𝑛
Correction
129
Chapitre IV. Partie A
Tableau IV.A. 1. Propriétés mécaniques et géométriques de la poutre sandwich à cœur viscoélastique en PVB
130
Chapitre IV
validé par la figure IV.A.7 dont les vibrations libres pour ℎ⁄𝐻 = 0.1 sont disparues plus
rapidement par rapport aux autres valeurs d’épaisseur.
Figure IV.A. 3. Réponses dynamiques de la poutre sandwich à cœur viscoélastique en PVB sous chargement
mobile pour différentes valeurs de ℎ⁄𝐻.
Figure IV.A. 4. Réponses de facteur d’amplification dynamique de la poutre sandwich à cœur viscoélastique pour
différentes valeurs de ℎ⁄𝐻
131
Chapitre IV. Partie A
La configuration de la poutre sandwich avec un traitement PCLD est présentée dans la figure
IV.A.9 dont la face supérieure est considérée comme étant la couche contraignante alors que la
face inférieure est considérée comme étant la couche de base. Le cœur viscoélastique de la couche
centrale est constitué de matériau polymère ISD112 dont le module viscoélastique est décrit par
le modèle de Maxwell généralisé et considéré pour deux températures différentes 20° et 27°. Le
module de Young dépendant fortement de la fréquence, est donné par le tableau IV.A.3. Les
propriétés mécaniques et géométriques des sandwichs sont présentées dans le tableau IV.A.4.
3
∆𝑗
𝐸2 ∗ (𝜔) = 2𝐺0 (1 + υ2 ) (1 + ∑ ) (IV. 7)
ω − iΩ𝑗
𝑗=1
20°𝐶 27°𝐶
j 𝐺0 ∆𝑗 Ω𝑗 𝐺0 ∆𝑗 Ω𝑗
1 2.8164 31.1176 0.746 468.7
2 5.11×104 13.1162 446.4542 5×105 3.265 4742.4
3 45.4655 5502.5318 43.248 71532.5
Tableau IV.A. 3. Propriétés mécaniques et géométriques du sandwich viscoélastique avec un module fortement
dépendant de la fréquence
132
Chapitre IV
Dans le but d’étudier l’effet de différents rapports d’épaisseur sur le comportement vibratoire
de la poutre sandwich, les vibrations forcées des sandwichs avec cœur viscoélastique à 20° et 27°
sous chargement mobile, sont y évaluées. La figure IV.A.10 présente les réponses forcées pour
deux sandwichs considérés à deux températures du cœur viscoélastique obtenues pour différentes
valeurs de rapport d’épaisseur de la couche viscoélastique. On remarque que les amplitudes
maximales des réponses pour le sandwich avec cœur viscoélastique à 20° sont obtenues pour les
rapports d’épaisseur h2/h3= 1/2, 3/2 et 1 dont le déplacement transversal atteint 𝑤 = 13.17 ×
10−3 , 13.83 × 10−3 et 13.13 × 10−3 , respectivement. Cependant, l’amplitude minimale a été
enregistré pour h2/h1= 1/4 dont le déplacement atteint 𝑤 = 12.88 × 10−3 . Les résultats du cœur
à 27° sont similaires à ceux de 20° dont les amplitudes de déplacements pour h 2/h1= 1/4, 1/2, 3/2
et 1, sont atteints les valeurs 𝑤 = 12.46 × 10−3 , 12.74 × 10−3 , 12.27 × 10−3 et 12.86 × 10−3 ,
respectivement. Cela révèle que l’amortissement devient plus efficace pour les faibles épaisseurs
de la couche viscoélastique.
Figure IV.A. 6. Réponses forcées de la poutre sandwich en appuis simples avec le cœur viscoélastique ISD112
considéré aux températures 20° et 27° obtenues pour différentes valeurs de rapport d’épaisseur entre la couche
viscoélastique et la couche d’inférieure h2/h3 ((a)- ISD112 à 20°, (b)- ISD112 à 27°).
Les réponses en amplitude de deux poutres sandwichs sous chargement mobile, obtenues pour
différentes valeurs de rapport d’épaisseur entre la couche inférieure de base et la couche
supérieure contraignante sont présentées dans les figures IV.A.11.a et IV.A.11.b, respectivement
pour les deux cœurs viscoélastiques considérés à 20 et 27°. Les amplitudes de déplacement
maximales du cœur à 20° correspondants aux quatre rapports d’épaisseur h3/h1= 1/4, 1/2, 1, et 2,
sont atteintes respectivement les valeurs 𝑤 = 0.1217, 0.1063, 0.0582 et 0.01288 alors que
pour le cœur à 27° les déplacements sont atteints les valeurs 𝑤=
133
Chapitre IV. Partie A
0.1022, 0.0851, 0.0472 et 0.1217. On remarque que les amplitudes maximales de déplacement
sont obtenues avec un rapport de h3/h1=1/4 pour les deux températures du cœur viscoélastique
alors que les amplitudes minimales sont obtenues pour h 3/h1=2. Ces résultats montrent que la
structure devient à fort pouvoir amortissant lorsque l’épaisseur de la couche de base est
considérable et strictement supérieure de la couche supérieure.
Figure IV.A. 7. Réponses forcées de la poutre sandwich en appuis simples avec le cœur viscoélastique ISD112
considéré aux températures 20° et 27° obtenues pour différentes valeurs de rapport d’épaisseur entre la couche
inférieure de base et la couche supérieure ((a)- ISD112 à 20° ; (b)- ISD112 à 27°)
134
Chapitre IV
Nous abordons dans cette section le comportement vibratoire sous l’action d’une charge
mobile des plaques sandwich à cœur viscoélastique dont plusieurs modèles sont examinés. Le
premier modèle est celui indépendant de la fréquence
Nous y rappelons que la même formulation par la méthode des éléments finis définie au
chapitre précédent est utilisée pour obtenir l’équation suivante de mouvement :
M U K U F vt IV.8
où [𝑀] et [𝐾 ] sont respectivement les matrices de masse et de rigidité globales, U est le vecteur
de déplacement nodal. {𝐹 (𝑣𝑡)} est le vecteur nodal défini en fonction de la charge mobile
bidirectionnelle.
136
Chapitre IV
Dans le but de passer de la formule définie dans le domaine fréquentiel Eq.(IV.8) au modèle
défini dans le domaine temporel et de trouver la réponse sous l’effet de la charge mobile, nous
décomposons la matrice de rigidité dépendante de la fréquence sous la forme suivante :
K K R iK I IV.9
avec 𝐾 𝑅 et 𝐾 𝐼 sont les parties réelle et imaginaire de la matrice de rigidité. En remplaçant la
nouvelle matrice 𝐾 (𝜔) dans (IV.8), on obtient :
M U ( K ]R i K
I
U F vt IV.10
Cette écriture peut être exprimée sous la forme de l’Eq. (IV.11) si nous considérons la
définition de la matrice d’amortissement équivalente définie en fonction de la partie imaginaire
de la matrice de rigidité.
M U C U K U F vt IV.11
R
Tout d’abord, nous devons résoudre en premier lieu le problème aux valeurs propres afin de
construire la base modale en utilisant le même algorithme de la méthode asymptotique
numérique. Nous résumons dans l’organigramme de la figure IV.B.2 les principales étapes de
cette méthode qui consiste à exprimer les variables du problème en fonction des séries de Taylor
et les remplacer dans le problème aux valeurs propres décomposé. La résolution du problème se
fait à partir d'un problème initial linéaire 𝑹(𝑼, 𝝀) = [[𝑲𝑹 ] − 𝝎𝟐𝟎 [𝑴]] 𝑼𝟎 pour 𝒂𝒋 = 𝟎 jusqu'au
problème non linéaire 𝑹(𝑼, 𝝀) = [[𝑲𝑹 ] + 𝑬𝒏 (𝝎)[𝑲𝑰 ] − 𝝎𝟐𝒏 [𝑴]] 𝑼𝒏 pour 𝒂𝒋 = 𝟏. On applique
donc une procédure de continuation qui consiste à définir une nouvelle tranche de solution à
partir d'un point de départ (𝑼𝒋 , 𝝀𝒋 ). La solution est obtenue à chaque itération en résolvant un
système d'équations linéaires. Le processus se termine lorsque la valeur 𝒂𝒋+𝟏 > 𝟏. Le principe
de la méthode est bien détaillé dans le chapitre II, section II.A.3. Ensuite, nous utilisons le
schéma d’intégration de Newmark présenté dans la figure IV.A.2 combiné avec la méthode de
137
Chapitre IV. Partie B
superposition modale pour obtenir la réponse temporelle de la plaque sandwich sous l’action de
la charge mobile.
Pour un comportement linéaire non vieillissant des matériaux viscoélastiques dont le principe
de superposition de Boltzmann est applicable (voir section I.2.2), la superposition des
sollicitations implique la superposition analogue des réponses. Ainsi, la réponse d'un matériau
138
Chapitre IV
0
avec 𝐷(𝑡) est la fonction de relaxation. Cette intégrale est considérée comme la somme d’une
réponse de l’élasticité instantanée et une autre de l’effet de mémoire qui tient compte l’histoire
de la déformation ε. Cette écriture peut être exprimée aussi sous la forme:
t
D t ε kl t d IV.13
0
M U K U K D t U d F vt IV.14
f c0
avec 𝐾𝑓 et 𝐾 𝑐0 sont les matrices de rigidité des couches de face et du cœur viscoélastique
respectivement. L’application de la méthode de Newmark pour résoudre ce type d’équation
nécessite quelques arrangements. De ce fait, le schéma implicite d’intégral direct à accélération
moyenne également appelée méthode pas à pas, est y appliquée pour remplacer les dérivées
temporelles en utilisant les relations suivantes :
1
U n1 U n tU n t 2U n t 2U n1 IV.15.a
2
¨ ¨
un1 U n 1 U n tU n1 IV.15.b
Les paramètres du schéma d’intégrale de Newmark à accélération moyenne sont définis
comme :
1 1
; IV.16
2 4
De ce fait, avec ces valeurs de 𝛾 et 𝛽, le schéma d’intégrale devient inconditionnellement
stable. En remplaçant ces valeurs dans la relation IV.15, le déplacement et la vitesse peuvent être
exprimés :
1 ¨ 1 ¨
U n1 U n tU n t 2 U n t 2 U n1 IV.17.a
4 4
139
Chapitre IV. Partie B
1 ¨ 1 ¨
U n1 U n U n tU n 1 IV.17.b
2 2
L’expression de l’accélération et la vitesse en fonction du déplacement permet de les
réécrire comme suit :
4
U n1 U n tU n U n 1 U n tU n 1 IV.18.a
t 2
2
U n1 U n
U n1 U n IV.18.a
t
Ainsi, l’intégrale apparaît dans l’équation de mouvement peut être exprimé en fonction
du n-ième pas de l’intégrale sous la forme :
t n 1 t
D t U d D n 1 t U d IV.19
0 0
en remplaçant les relations Eqs. (IV.18) et Eq. (IV.20) dans l’équation de mouvement Eq.
(IV.14), on obtient :
n 1 t
M U n tU n 2 U n1 U n tU n1 K f U n K c0 D n 1 t p U F
4
n 1
t
p
p 0
(IV.21)
cela permet de déterminer le déplacement en résolvant le système d’équations suivant :
4 4 ¨
f c0 c0
t 2
M K K D 0 U n 1 Fn 1 t 2
M K D 0 U n M U n
4 D n 1 t
IV.22
c0
c0
t M K D 0 U n K t U n D nt U n D t U n
2
P N v t , v t pour t
avec fi P0 x vxt y v yt N w x, y dx 0 w x y IV.23.c
L
0 pour t
0
𝛿 est la fonction Delta de Dirac, 𝑃0 l’amplitude de la force, 𝑣𝑥 et 𝑣𝑦 sont les vitesses constantes
de la charge mobile suivant les axes 𝑥 et 𝑦, respectivement et 𝑛 est le nombre des nœuds.
140
Chapitre IV
𝜏 est le temps nécessaire pour que la charge parcoure tous les nœuds du chemin défini par :
l
nele
(IV.24)
v v
avec 𝑛𝑒 et 𝑙𝑒 est le nombre d’éléments et la longueur d’élément à parcourir par la charge mobile
et 𝜏𝑙 est la distance parcourue par la charge mobile.
Selon le principe utilisé auparavant pour la poutre sandwich, la force nodale 𝑓𝑖 (𝑖 =
1, 2, . . . , 𝑛) se déplace sur tous les nœuds du chemin, lorsque 0≤t≤τ comme montre la figure
IV.B.3 à laquelle la valeur de la charge est f(i=j)=P0 au j-ème nœud pour 𝑡𝑗 = (𝑗 − 1)𝑙𝑒 ⁄𝑣, tandis
que f(i≠j)=0 sur le reste des nœuds.
Dans cette section, nous examinons les réponses vibratoires des plaques sandwich à cœur
viscoélastique avec deux modèles viscoélastiques dont le premier est dépendant de la fréquence
et le deuxième dépendent du temps. Les deux modèles sont examinés avec la même configuration
de la plaque sandwich dont les propriétés mécaniques et géométriques sont présentées dans le
tableau IV.B.1.
141
Chapitre IV. Partie B
X
P0 P0 X
yp b
L
L
Y
Y
142
Chapitre IV
instant notée 𝑿 est normalisée par rapport à la distance totale parcourue par la charge mobile
𝜏𝑙 = 𝑛𝑒 𝑙𝑒 . Il est très évident que les amplitudes des réponses augmentent quand la vitesse de la
charge augmente aussi. Cela peut être expliqué par la vitesse critique de résonance à laquelle les
amplitudes peuvent atteindre de grandes valeurs.
Figure IV.B. 5. Réponses de la plaque sandwich à cœur viscoélastique ISD112 à 27° en fonction de la vitesse de la
charge mobile ((a)- Mouvement rectiligne ; (b)- Mouvement en zigzag)
Dans le but d’expliquer le phénomène de résonance, les réponses de RMS (Root Mean
Square) et de Rmax (réponse maximale) en fonction de la vitesse de la charge mobile sont
présentés dans la figure IV.B.6 et figure IV.B.7, respectivement. Ces réponses traduisent l’effet
de résonance sur les réponses de vibrations. Par conséquent, nous observons une augmentation
des amplitudes jusqu’à ce qu’elles atteignent des valeurs maximales correspondantes aux vitesses
critiques avant de décroitre toute en éloignant de ces valeurs. Cependant, nous y observons que
les réponses de RMS diminuent lentement au-delà de la vitesse critique à cause de l’effet de la
charge mobile à vitesse élevée lorsqu’elle traverse la plaque, cela provoque généralement des
oscillations à amplitudes élevées dans la région des vibrations libres.
143
Chapitre IV. Partie B
Figure IV.B. 6. Réponses RMS de la plaque sandwich à cœur viscoélastique ISD112 à 27° en fonction de la vitesse
de la charge mobile ((a)- Mouvement rectiligne ; (b)- Mouvement en zigzag)
Figure IV.B. 7. Réponses Rmax de la plaque sandwich à cœur viscoélastique ISD112 à 27° en fonction de la
vitesse de la charge mobile ((a)- Mouvement rectiligne ; (b)- Mouvement en zigzag)
De ce fait, afin de valider cette dernière interprétation, nous examinons les réponses de
déplacement avec différentes vitesses. Les plages de vitesse sont choisies pour comparer les
réponses avec des vitesses inférieures et supérieures à la vitesse critique sachant que la vitesse
critique peut être calculée pour le premier mode dominant comme suit :
vc l ω0 (IV.25)
avec ω0 est la fréquence propre.
La vitesse critique pour la trajectoire rectiligne est aux alentours de vc=60 m/s alors que pour
la trajectoire en zigzag considéré, la vitesse critique atteint vc=80 m/s. Cependant, à partir des
résultats obtenus présentés dans la figure IV.B.8, nous pouvons remarquer que les amplitudes
des déplacements sont clairement plus faibles dans les deux régions des vibrations forcées et
144
Chapitre IV
libres lorsque la vitesse est inférieure à la vitesse critique. Cependant, les amplitudes atteignent
des valeurs maximales avec la vitesse critique. En revanche, les amplitudes commencent à
diminuer au-delà de la vitesse critique, mais les oscillations libres s’atténuent lentement par
rapport aux oscillations induisent par des vitesses inférieures à la vitesse critique.
Figure IV.B. 8. Réponses de la plaque sandwich à cœur viscoélastique ISD112 sous la charge mobile de
magnitude 100 N obtenues pour différentes valeurs de vitesse ((a)- Mouvement rectiligne ; (b)- Mouvement en
zigzag)
Dans le but de comparer les réponses de la plaque sandwich avec différents types de cœurs
viscoélastiques, nous examinons trois cœurs viscoélastiques présentés auparavant qui sont
ISD112 à 27°, ISD112 à 20° et PVB. Les réponses de déplacement transverse pour un
mouvement rectiligne de la charge mobile sont présentées dans la figure IV.B.9. Il est très clair
que les amplitudes des réponses obtenues avec les cœurs ISD112 à 27° et ISD112 à 20° sont plus
élevées par rapport à celles obtenues avec le cœur en PVB. Par conséquent, les mêmes
conclusions sont obtenues pour les réponses de RMS et Rmax présentées dans la figure IV.B.10
dont les réponses des cœurs en ISD112 sont plus élevées par rapport aux réponses de la plaque
sandwich avec le cœur viscoélastique en PVB. Ces résultats illustrent que le matériau
viscoélastique PVB est plus rigide (module élastique élevé) par rapport aux autres cœurs
viscoélastiques en ISD112. Cela montre que l’atténuation des amplitudes de vibrations, ne
dépend pas seulement du module de dissipation, mais aux aussi de la partie élastique du cœur
viscoélastique qui renforce la rigidité globale de la structure.
145
Chapitre IV. Partie B
Figure IV.B. 9. Réponses sous chargement mobile de la plaque sandwich avec les trois types de cœurs
viscoélastiques ISD112 à 27°, ISD112 à 20° et PVB
Figure IV.B. 10. Comparaison des réponses de RMS (a) et de Rmax (b) de la plaque sandwich avec les trois types
de cœurs viscoélastiques ISD112 à 27°, ISD112 à 20° et PVB
La deuxième partie de cette section concerne le comportement des plaques sandwich à cœur
viscoélastique dépendant du temps. Nous y considérons la plaque sandwich présentée dans le
tableau IV.4. La fonction de relaxation du cœur dépendant du temps est donnée par [32]:
146
Chapitre IV
structure particulièrement dans la zone de vibrations libres. Cela confirme ce qui a été constaté
auparavant que l’augmentation de l’épaisseur de la couche viscoélastique n’implique pas
forcément le renforcement la capacité dissipative de la structure.
Figure IV.B. 11. Réponses sous chargement mobile du sandwich à cœur dépendant du temps
Figure IV.B. 12. Réponses de Rmax de la plaque sandwich à cœur dépendant du temps
147
Chapitre IV. Conclusions
Conclusions
Dans ce chapitre, une approche numérique a été présentée pour étudier les vibrations forcées
sous chargement mobile des structures sandwich avec âme viscoélastique. La méthode des
éléments finis a été utilisée pour discrétiser l’équation gouvernant le mouvement de vibration.
Les réponses dynamiques sont obtenues en utilisant la méthode de superposition modale et le
schéma d’intégration implicite de Newmark pour résoudre l’équation de mouvement. Ainsi,
différentes lois de comportement viscoélastique y ont été étudiées entre autres ceux dépendants
de la fréquence. Cependant, la difficulté de résoudre le problème de valeur propre en raison de
la dépendance en fréquence de la matrice de rigidité nous a ramené à utiliser la méthode
asymptotique numérique afin d’obtenir la base modale.
A partir des travaux réalisés dans ce chapitre, plusieurs conclusions peuvent être tirées.
Poutre sandwich
Pour la poutre sandwich, les différentes analyses du comportement sous l’effet de la charge
mobile nous ont permis de conclure :
148
Chapitre IV
Plaque sandwich
De ce fait, l’analyse des résultats prouve que l’amélioration des propriétés amortissantes,
nécessite non seulement de maitriser les paramètres de chargement tel que la vitesse de la charge
mobile, mais aussi les propriétés mécaniques et géométriques de la structure.
149
Chapitre V
Sommaire
Introduction
Dans ce chapitre, une approche numérique est présentée pour étudier les réponses vibro-
acoustiques des plaques sandwichs avec cœur viscoélastique dans le but d’analyser et d’examiner
l’amortissement structurel apporté par la couche viscoélastique en considérant différents modèles
viscoélastiques décrits par différentes lois de comportement. L’étude numérique est basée sur la
méthode des éléments de frontières et la méthode des éléments de finis dans la formulation et la
discrétisation de l’équation de mouvement établie par le principe de Hamilton ainsi que sur la
méthode asymptotique numérique afin de résoudre le problème aux valeurs propres et de
construire une base modale de la structure.
151
Chapitre V
Dans cette partie, nous présentons une approche pour étudier les réponses vibro-acoustiques
des plaques sandwich à cœur viscoélastique. La distribution de la pression acoustique dans le
domaine du fluide semi-plan faisant face à la couche supérieure de la plaque sandwich, peut être
obtenue à partir de la solution de l'équation de Helmholtz [120]:
2 2 2
2 k 2 p x, y , z , 0 V.1
x y z
2 2
La condition de rayonnement de Sommerfeld doit être satisfaite pour que l'amplitude de
l'onde disparaisse à l'infini :
] lim p 1 p 0 (V.2)
r
r c0 r
p x, z,
| 0 2v t , (V.3)
r
h
z
2
Figure V.1. Illustration schématique de la transmission du son à travers un panneau dans baffle plan et infini
rigide ((a) vue latérale ; (b) vue globale)
152
Chapitre V
Cette étude porte sur deux modèles viscoélastiques présentés auparavant au chapitre III dont
le premier présenté à la section III.A.1.1 est indépendant de la fréquence et le deuxième est un
cœur viscoélastique en DYAD606 présenté à la section III.A.1.2.2. La plaque sandwich est
soumise à une pression acoustique d’incidence d’amplitude ps=1 Pa. Ainsi, L’étude comporte
l’effet de la variation de température sur les différentes réponses acoustiques avec une
température ambiante T0=25 °C.
V.B.2.1. Modèle de comportement à module constant
Dans cette section, nous examinons les réponses acoustiques des plaques sandwichs avec
cœur viscoélastique considéré indépendant de la fréquence, décrit auparavant par l’Eq. (III.2).
Nous considérons les mêmes propriétés mécaniques et géométriques de la plaque sandwich
présentées au tableau III.B.1.
L’effet de la variation de la température sur la pression acoustique est présenté dans la figure
V. 2, la figure montre que les pics de pression acoustique ont tendance à s’atténuer et à se diriger
vers les fréquences plus basses avec l'augmentation de la charge thermique. La raison est
évidemment, que la température adoucit la rigidité de la structure et diminue les fréquences
153
Chapitre V
Figure V. 2. Réponses fréquentielles de la pression acoustique de la plaque sandwich avec cœur viscoélastique
indépendant de la fréquence
Figure V. 3. Réponses fréquentielles de la perte de transmission acoustique « STL » de la plaque sandwich avec
cœur viscoélastique indépendant de la fréquence
Afin d’étudier les réponses acoustiques sous l’effet de différents facteurs d’amortissement
viscoélastique, les différentes réponses de pression acoustique pour différentes valeurs de facteur
154
Chapitre V
d’amortissement ηv=0.1, 0.5, 1.0 et 1.5, sont présenté dans la figure V. 4 avec différentes valeurs
de température de chargement ΔT=0, 10, 20 et 30 °C. Il est apparu que les pics d’amplitudes
diminuent avec l’augmentation de facteur d’amortissement dont une atténuation complète a été
observée pour ηv=1.5. Ainsi, les réponses de STL (figure V.5) montrent que la structure devient
plus amortissante quand les valeurs du facteur de perte viscoélastique sont élevées, atténuant par
conséquent l’effet d’irradiation acoustique.
155
Chapitre V
Figure V. 5. Réponses fréquentielles de STL de la plaque sandwich à cœur viscoélastique indépendant en appuis
simples pour différentes valeurs de facteur de perte ((a)- ΔT=0°C, ΔT=10°C, ΔT=20°C, ΔT=30°C)
156
Chapitre V
Figure V. 6. Réponses fréquentielles de la pression acoustique de la plaque sandwich avec viscoélastique à cœur
viscoélastique en DYAD606
Figure V. 7. Réponses fréquentielles de STL de la plaque sandwich avec viscoélastique à cœur viscoélastique en
DYAD606
157
Chapitre V. Conclusion
Conclusions
Dans ce chapitre, une approche numérique basée sur une théorie d’ordre élevé des plaques
sandwichs à cœur viscoélastique a été présentée afin de caractériser les réponses vibro-
acoustiques des poutres sandwichs à cœur viscoélastique. L’objectif principal était donc
d’examiner l’efficacité de l’amortissement des matériaux viscoélastiques décrits par des lois de
comportement dépendant non seulement de la fréquence, mais aussi à la température. Par
conséquent, différents modèles viscoélastiques y ont été examinés en tenant en compte leurs
propriétés mécaniques spécifiques.
D’après les résultats obtenus, les conclusions suivantes peuvent être tirées :
158
Bilan, conclusions et perspectives
159
Bilan, conclusions et perspectives
Bilan
L’objective de cette étude était la maitrise des phénomènes vibro-acoustiques des structures
mécaniques munies de matériaux viscoélastiques en examinant les réponses vibro-acoustique,
ainsi, afin d’évaluer l’efficacité de traitement passif à amortir les vibrations et les radiations
acoustiques. Les matériaux viscoélastiques utilisés dans les structures composites contribuent
dans l’atténuation des vibrations qui peuvent provenir de divers chargements dynamiques. Ces
matériaux viscoélastiques sont capables de transformer l’énergie de vibration en énergie
thermique ce qui permet de dissiper l’énergie de déformation sous forme de chaleur.
Par conséquent, le comportement vibratoire libre et forcée ainsi que le comportement vibro-
acoustique des structures poutres et plaques sandwichs à cœur viscoélastique avec différentes
lois de comportement, différentes configurations, différentes conditions aux limites et sous divers
chargements dynamiques (mobile, concentré harmonique), ont été caractérisés et analysés en
tenant en compte aussi l’effet de la température.
L’approche numérique utilisée dans ces travaux de thèse est basée généralement sur la
méthode des éléments finis et la méthode des éléments de frontières pour la formulation de
l’équation de mouvement et la méthode asymptotique numérique combinée à la méthode
différenciation automatique pour obtenir la base modale servant à la résolution de l’équation de
mouvement.
Tout d’abord, la validation de l’approche numérique a été réalisée dans le chapitre II en
comparant les résultats avec d’autres références. Ensuite, une étude paramétrique dans le
domaine fréquentiel a été effectué au chapitre III pour évaluer l’amortissement structurel des
sandwichs munis de matériaux viscoélastiques et chercher la configuration optimale des poutres
et plaques sandwich. Au chapitre IV, les réponses dynamiques des poutres et plaques sandwichs
à cœur viscoélastique sous chargement dynamique mobile sont obtenues en considérant
différentes lois viscoélastiques dépendant non seulement de la fréquence mais aussi du temps et
plusieurs configurations des structures sandwich. Cependant, pour la plaque sandwich deux types
de mouvement sont examinés pour la structure plaque dont le premier est linéaire selon une ligne
160
Bilan
Conclusions
Les structures abordées durant les travaux réalisés sont deux types dont la première est une
poutre sandwich et la deuxième est une plaque sandwich. Les deux types de structures sont munis
de matériaux viscoélastiques. Par conséquent, les chapitres ont traité chacun de structure de
manière indépendante tout en adoptant des théories et des applications adéquates en tenant en
compte les propriétés géométriques pour chaque type de ces structures. Les conclusions de
chaque chapitre sont aussi divisées en deux parties dont la première concerne les études établies
des structures poutres sandwichs et la deuxième pour les structures plaques. A partir des travaux
réalisés, nombreuses conclusions peuvent y être tirées concernent les deux types de structures
étudiées.
161
Conclusions
L’efficacité de traitement passif d’amortissement était très évidente pour les différentes
réponses vibratoires et vibro-acoustique après avoir réalisé étude comparative de différents types
de cœurs viscoélastiques (PVB, ISD112 à 20° et ISD112 à 27°) avec un matériau purement
élastique. D’un autre côté, les propriétés d’amortissement se diffèrent d’un matériau
viscoélastique à un autre dont les cœurs viscoélastiques en ISD112 ont plus de capacités
dissipatives par rapport au cœur en PVB.
Ainsi, l’effet de la vitesse de la charge mobile est très important lorsque cette vitesse est près
de la vitesse critique tout en observant une augmentation significative des amplitudes vibrations
de différentes structures.
Perspectives
Les travaux de cette thèse de recherche sont réalisés pour surmonter les problèmes liés aux
vibrations structurelles et les bruits acoustiques qui peuvent altérer les fonctions des systèmes en
appliquant la technique d’amortissement passif aux structures sandwich. Tandis que ces travaux
ne sont qu’une goutte d’eau dans l’océan et ils ouvrent la voie de nombreux axes de recherche.
162
Perspectives
De point de vue numérique aussi, l’efficacité de traitement passif abordé est très efficace
quand il s’agit de moyennes et hautes fréquences à la réciproque de la technique d’amortissement
actif qui est plus efficace aux faibles fréquences. L’application de la technique hybride « Passif
+ Actif) sera plus efficace tenant en compte toutes les plages de fréquences. L’élaboration d’un
modèle interactif de la technique hybride s’apparaît comme un objectif très important à atteindre.
Ainsi, l’optimisation des paramètres de conception est un axe aussi important à explorer, une
procédure d'optimisation numérique est une approche décisive pour obtenir une véritable
configuration optimale globale du système non seulement à l’échelle macroscopique, mais aussi
à l’échelle microscopique en parlant sur les micros/nano structures. Cette étude repose sur
l’optimisation des paramètres matériels et géométriques ainsi que les conditions de
fonctionnement de ces structures.
Le troisième axe qui semble aussi d’importance est le comportement non-linéaire des
structures sandwich. En réalité, plusieurs travaux ont été portés sur cet axe de recherche, mais
toujours avec contraintes et restrictions empêchant à maitriser le véritable problème. L’utilisation
d’une approche basée sur une base modale linéaire et à un seul mode doit être évitée en utilisant
par exemple les outils du numérique avancé. Il serait envisageable d’introduire du concept de
l’intelligence artificielle « Artificiel Intelligence AI» tels que les machines d’apprentissage
« Machine Learning ML » dans la résolution du problème est une étape importante qui peut
affranchir tous types de difficultés rencontrées en raison de l’efficacité de ces méthodes ainsi
qu’ils représentent une tendance actuelle de la recherche.
Du point de vue expérimental, les essais expérimentaux dans le domaine d’ingénierie sont
toujours nécessaires pour valider n’importe quelle approche numérique basée sur la théorie
fondamentale. Les modèles du comportement viscoélastique sont fortement dépendants du
temps, de la fréquence et de la température ainsi que leurs données expérimentales sont
pertinentes et rares, il est nécessaire de mener des essais expérimentaux caractéristiques pour
avoir des modèles plus précis dans le but de mieux maitriser dans les approches numériques
163
Annexe
ANNEXE A
ANNEXE A
U K W 0 A.1
où 𝛱 est l'énergie potentielle totale, 𝑈,𝐾 et 𝑊 sont l'énergie de déformation, l’énergie cinétique
et l'énergie des forces extérieures, respectivement.
Le principe de Hamilton est une forme intégrale dans le temps du principe du travail virtuel.
Il est obtenu par une transformation intégrale de l'expression :
T T
165
ANNEXE A. Principe de formulation variationnelle
La condition aux limites de Dirichlet, est la solution 𝑢(𝑥) associée à une équation
différentielle ou à une équation aux dérivées partielles sur les frontières d’un domaine compact,
par exemple :
u( x 0) A.3
La condition aux limites de type Neumann, est la dérivée de la solution 𝑢′(𝑥) associée à une
équation différentielle ou à une équation aux dérivées partielles sur les frontières d’un domaine
compact, par exemple :
u
A.4
x x 0
La condition aux limites de type Mix, est donnée en fonction de la solution 𝑢(𝑥) et la dérivée
normale de la solution 𝑢′(𝑥) associée à une équation différentielle ou à une équation aux dérivées
partielles sur les frontières d’un domaine compact, par exemple :
u
u A.5
x
Nabla noté ∇, est un opérateur utilisé en mathématiques, notamment en calcul vectoriel, comme
opérateur différentiel vectoriel présenté par un vecteur de dérivés partiels.
166
ANNEXE A
, , A.6
x1 y1 z1
grad (u ) u A.6.a
div(u ) .u A.6.b
rot (u ) u A.6.c
Premier ordre
Pour 𝑢 ⃗ 𝜙, l’équation (A.7.f) devient :
⃗ = 𝑔∇
V S
Deuxième ordre
Pour 𝑢
⃗ = 𝑔∇ ⃗⃗⃗⃗⃗ 𝑔, l’équation (A.7.f) devient :
⃗ 𝜙 − 𝜙∇
V S
La distribution de Dirac ou la fonction Delta noté 𝛿(𝑥) est une fonction qui prend une valeur
infinie en 𝑥 = 0, et la valeur zéro partout ailleurs dont l'intégrale sur ℝ est égale à 1.
167
ANNEXE B. Méthode des élément finis
ANNEXE B
La méthode des éléments finis est l’un des outils de mathématiques appliquées permettant de
rechercher une solution approchée d’un problème décrit par une équation différentielle (ED) ou
une équation aux dérivées partielles (EDP) sur un domaine compact avec conditions aux bords
et/ou à l'intérieur du domaine du type Dirichlet (valeurs aux limites) ou Neumann (gradients aux
limites) ou de Robin (relation gradient/valeurs sur les frontières). La méthode des éléments finis
est basée généralement sur la formulation variationnelle.
B.1. Discrétisation
La discrétisation consiste à subdiviser le domaine 𝛺 en un ensemble d'éléments disjoints 𝛺𝑒 .
Chaque élément 𝛺𝑒 est défini par un ensemble de nœuds dont le nombre est 𝑛𝑛𝑒 .
ne
i1
e (B.1)
B.1. Approximation
En éléments finis, le champ de variable est remplacé par un champ décrit en fonction de
nombre finis des nœuds :
nnt
u r ( x, y, z ) u app ( x, y, z ) Ni ( x, y, z ) un (B.2)
i 1
168
ANNEXE B
du modèle d’éléments finis. Ces fonctions doivent être continuées et doivent satisfaire les valeurs
nodales (figure B.2) :
1 pour j=i
Ni ( x j , y j , z j ) (B.3)
0 pour j i
La construction des fonctions de forme (interpolations) dépend aussi de la nature du problème
et la précision des résultats demandés. Plusieurs méthodes sont existes pour construire ces
fonctions en se basant sur le triangle du Pascale afin de choisir le degré des polynômes
d’interpolation.
Ainsi, les polynômes de Lagrange sont l’une des méthodes d’approximation les plus utilisés
en éléments finis pour construire des fonctions de forme. Elles y offrent une qualité très
acceptable. Ces polynômes sont donnés pour une approximation unidimensionnelle par :
x xi
Ni inn1,e i j (B.4)
xi x j
Où 𝑛 est le degré du polynôme.
x , , Ti , , x e (B.5.a)
y , , T , , y
i
e
(B.5.b)
z , , T , , z
i
e
(B.5.c)
Le nombre des points géométriques, est déterminant dans le choix de type de transformation.
Pour cela, nous pouvons distinguer trois types de transformations, la première est la
169
ANNEXE B. Méthode des élément finis
transformation sub-paramétrique, dans ce cas le nombre des points géométriques est inférieur au
nombre de nœuds, la deuxième est la transformation isoparamétrique dont le nombre est égal au
nombre de nœuds. Le troisième type est la transformation super-paramétrique dont le nombre
des points géométriques, est supérieur au nombre de nœuds.
x , , N i , , x e (B.7.a)
y , , N , , y
i
e
(B.7.b)
z , , N , , z
i
e
(B.7.c)
Transformation iso-paramétrique
Point géométrique
Nœud
y Elément physique
𝜂
Elément local
x
Figure B. 1. Exemple d’une transformation isoparamétrique plane
170
ANNEXE B
N1 N1
x
N1 N1
J (B.8)
y
N1 N1
z
ng N1
ng
N1
ng
N1 N N ng x1 y1 z1
i 1
xi i 1
yi
i 1
zi
1 N 2
...
x2 y2 z 2
ng
N1
ng
N1
ng
N1 N1 N 2 N ng . .
J xi yi zi ... (B.9)
. .
i 1 i 1 i 1
ng N1
ng
N1
ng
N1 N1 N 2 N ng . .
i 1
xi
i 1
yi
i 1
zi
...
xng yng z3
avec (𝑥𝑖 , 𝑦𝑖 , 𝑧𝑖 ) sont les points géométriques de l’élément réel et 𝑛𝑔 est le nombre de nœuds
de l’élément.
Puisque les champs de déplacements sont définis en fonction des coordonnées de l’élément
local, les formes d’intégration doivent être effectuées en fonction de l’élément parent. Le volume
élémentaire 𝑑𝑉𝑝 de l'élément physique est lié au volume élémentaire 𝑑𝑉𝑝 de l'élément de
référence par :
Le Le
(...)dV (...) J d d d
0
p
0
(B.11)
B.4. Assemblage
La procédure d'assemblage consiste à établir la somme des matrices et des vecteurs du
comportement. La continuité de la variable entre les éléments doit être assurée puisque la
procédure de maillage est reposée sur les liens de connectivité entre les éléments de maillage.
171
ANNEXE C. Méthode Asymptotique Numérique (MAN)
ANNEXE C
Nombreux problèmes en mécanique peuvent être discrétisés et présentés sous la forme d’un
problème de résidus suivant:
R(U , ) 0 (C.1)
Où 𝑅 est un système d’équations d’ordre 𝑛 × 𝑛 dont (𝑈, 𝜆) sont les solutions. Les méthodes
asymptotiques numériques y sont largement appliquées [130]-[131]. Elles reposent sur les
techniques de développement en séries entières tronquées à ordre bien défini les courbes de
solutions (𝑈, 𝜆) par rapport à un paramètre de chemin noté 𝛼.
N
U U j U n C.2
n
n 1
N
n C.3
j n
n 1
L’idée est de générer une séquence de morceaux de solutions dans le cadre d’une technique
appelée « technique de perturbation » en commençant d’un problème initial 𝛼𝑗 = 0 vers un
problème non linéaire pour 𝛼𝑗 = 1.
172
ANNEXE C
Chaque morceau de solution n’est valide que pour 𝛼 ∈ [0, 𝛼𝑚𝑎𝑥 ] , où 𝛼𝑚𝑎𝑥 est le rayon de
convergence du morceau de solution, donnée par :
1
U1 1
N
max (C.4)
U n
Pour que le morceau soit la solution du système, ce dernier doit vérifier l’équation (C.1) qui
s’écrit sous la forme :
avec:
Rn 0 pour n 1 C.6
173
ANNEXE C. Méthode Asymptotique Numérique (MAN)
R R
Up p Fpnl pour n 1 C.7
U j j
où le second membre 𝐹𝑝𝑛𝑙 dépend des termes d’ordre inférieur du couple (𝑈, 𝜆).
U1 12 0 C.9
2
U p .U1 p 1 0 C.10
par conséquent, la résolution du problème Eq.(C.1) passe par le calcul des membres 𝐹𝑝𝑛𝑙 et la
résolution de 𝑁 équations du système Eq.(C.7).
174
ANNEXE D. Réponses 3D de la plaque sandwich en ISD112 à 27° sous chargement mobile
ANNEXE D
Dans cette section, les réponses de déplacement de tous les nœuds du maillage de la plaque
sandwich étudiée à la section IV.B.3.1, sont présentées pour chaque instant dont 0< 𝑡 < 5𝜏𝑙 .
175
ANNEXE D. Réponses 3D de la plaque sandwich en ISD112 à 27° sous chargement mobile
Figure D. 1. Réponses de la plaque sandwich à cœur viscoélastique ISD112 à 27° sous chargement mobile de mouvement rectiligne et vitesse constante pour 0 < 𝑋/𝜏 < 0.4 (vibrations forcées)
176
ANNEXE D
Figure D. 2. Réponses de la plaque sandwich à cœur viscoélastique ISD112 à 27° sous chargement mobile de mouvement rectiligne et vitesse constante pour 0.4 < 𝑋/𝜏 < 0.85 (vibrations forcées)
177
ANNEXE D. Réponses 3D de la plaque sandwich en ISD112 à 27° sous chargement mobile
Figure D. 3. Réponses de la plaque sandwich à cœur viscoélastique ISD112 à 27° sous chargement mobile de mouvement rectiligne et vitesse constante pour 0.85 < 𝑋/𝜏 < 1 (vibrations forcées)
178
ANNEXE D
Figure D. 4. Réponses de la plaque sandwich à cœur viscoélastique ISD112 à 27° sous chargement mobile de mouvement rectiligne et vitesse constante pour 1 < 𝑋/𝜏 < 2.4 (vibrations libres)
179
ANNEXE D. Réponses 3D de la plaque sandwich en ISD112 à 27° sous chargement mobile
Figure D. 5. Réponses de la plaque sandwich à cœur viscoélastique ISD112 à 27° sous chargement mobile de mouvement rectiligne et vitesse constante pour 2.5 < 𝑋/𝜏 < 3.6 (vibrations libres)
180
ANNEXE D
Figure D. 6. Réponses de la plaque sandwich à cœur viscoélastique ISD112 à 27° sous chargement mobile de mouvement en zigzag et vitesse constante pour 0 < 𝑋/𝜏 < 0.27 (vibrations forcées)
181
ANNEXE D. Réponses 3D de la plaque sandwich en ISD112 à 27° sous chargement mobile
Figure D. 7. Réponses de la plaque sandwich à cœur viscoélastique ISD112 à 27° sous chargement mobile de mouvement en zigzag et vitesse constante pour 0.29 < 𝑋/𝜏 < 0.58 (vibrations
forcées)
182
ANNEXE D
Figure D. 8. Réponses de la plaque sandwich à cœur viscoélastique ISD112 à 27° sous chargement mobile de mouvement en zigzag et vitesse constante pour 0.58 < 𝑋/𝜏 < 0.85 (vibrations
forcées)
183
ANNEXE D. Réponses 3D de la plaque sandwich en ISD112 à 27° sous chargement mobile
Figure D. 9. Réponses de la plaque sandwich à cœur viscoélastique ISD112 à 27° sous chargement mobile de mouvement rectiligne par morceaux (mouvement zigzag) et vitesse constante pour
0.9 < 𝑋/𝜏 < 1 (vibrations forcées)
184
ANNEXE D
Figure D. 10. Réponses de la plaque sandwich à cœur viscoélastique ISD112 à 27° sous chargement mobile de mouvement en zigzag et vitesse constante pour 1 < 𝑋/𝜏 < 2.2 (vibrations libres)
185
ANNEXE D. Réponses 3D de la plaque sandwich en ISD112 à 27° sous chargement mobile
Figure D. 11. Réponses de la plaque sandwich à cœur viscoélastique ISD112 à 27° sous chargement mobile de mouvement en zigzag et vitesse constante pour 2.2 < 𝑋/𝜏 < 3.35 (vibrations
libres)
186
Références
Références
Références
[1]. P. Liénard. Etude d’une méthode de mesure du frottement intérieur de revêtements plastiques
travaillant en flexion. Rech Aéronaut.; 20:11-22, 1951.
[2]. E. M. Kerwin. Damping of flexural waves by a constrained viscoelastic layer. The Journal of
the Acoustical Society of America, 31(7): 952–962, 1959.
[3]. D. Ross, E. M. Kerwin. Damping of flexural vibration in plates by free and constrained
viscoelastic layers. Technical report, BBN Report 623, 1959.
[4]. D.J. Mead, D.C.Ae. The effect of a damping compound on jet-efflux excited vibrations: an
article in two parts presenting theory and results of experimental investigation. Aircraft
engineering and aerospace technology. Aircraft Engineering and Aerospace Technology, 32:
106 – 13, 1960.
[5]. E. E. Ungar. Loss factors of viscoelastic ally damped beam structures. The Journal of the
Acoustical Society of America, 34: 1082–1089, 1962.
[6]. Y. Y. Yu. Damping of flexural vibrations of sandwich plates. Journal of Aerospace Science,
29:790–803, 1962.
[7]. R.A. DiTaranto. Theory of vibratory bending for elastic and viscoelastic layered finite length
beams. Transactions of the American Society of Mechanical Engineers, Journal of Applied
Mechanics, 87:881–886, 1965.
[8]. D.J. Mead, S. Markus. The forced vibration of a three-layer, damped sandwich beam with
arbitrary boundary conditions. Journal of Sound and Vibration, 10: 163-75, 1969.
[9]. D.J. Mead, S. Markus. Loss factors and resonant frequencies of encastre damped sandwich
beams. Journal of Sound and Vibration, 12: 99-112, 1970.
[10]. M.J. Yan, EH Dowell. Governing Equations for Vibrating Constrained-Layer Damping
Sandwich Plates and Beams. Journal of Applied Mechanics, 39: 1041-6, 1972.
[11]. D.K. Rao. Frequency and loss factors of sandwich beams under various boundary
conditions. Mechanica Engineering Science, 20:271–282, 1978.
[12]. M.D. Rao, R. Echempati, S. Nadella. Dynamic analysis and damping of composite
structures embedded with viscoelastic layers. Composites Part B: Engineering, 28:547-54,
1997.
[13]. R.F. Kristensen, K.L. Nielsen, L.P. Mikkelsen. Numerical studies of shear damped
composite beams using a constrained damping layer. Composite Structures, 83: 304-11,
2008.
[14]. N. Kumar, S.P. Singh. Experimental study on vibration and damping of curved panel
treated with constrained viscoelastic layer. Composite Structures, 92:233-43, 2010.
[15]. B.R. Sher, Moreira RAS. Dimensionless analysis of constrained damping treatments.
Composite Structures, 99: 241-54, 2013.
[16]. M.D. Rao. Recent applications of viscoelastic damping for noise control in automobiles
and commercial airplanes. Journal of Sound and Vibration, 262: 457-74, 2003.
188
Références
[17]. A.E.H. Love. The mathematical theory of elasticity, 4th edition. Cambridge University
Press, 1927.
[18]. R.D. Mindlin. Influence of rotatory inertia effect and shear deformation in flexural motions
of isotropic elastic plates. Journal of Applied Mechanics, 18: 1031-1036, 1951.
[19]. M. Touratier. An efficient standard plate theory. International Journal of Engineering
Science, 29: 901-916, 1991.
[20]. J. N. Reddy. A simple higher-order theory of laminated composite plate. Journal of Applied
Mechanics, 51: 745-752, 1984.
[21]. J. Wang, Y. Xu, W. Zhang. Finite element simulation of PMMA aircraft windshield against
bird strike by using a rate and temperature dependent nonlinear viscoelastic constitutive
model. Composite Structures, 108: 21-30, 2014.
[22]. A. Arikoglu, I. Ozkol. Vibration analysis of composite sandwich beams with viscoelastic
core by using differential transform method. Composite Structures, 92: 3031-3039, 2010.
[23]. L. Irazu, M.J. Elejabarrieta. The effect of the viscoelastic film and metallic skin on the
dynamic properties of thin sandwich structures. Composite Structures, 176: 407-41 2017.
[24]. E.M. Daya, M. Potier-Ferry. A numerical method for nonlinear eigenvalue problems
application to vibrations of viscoelastic structures. Computers and Structures, 79 533-541,
2001.
[25]. E.M. Daya, L. Azrar, M. Potier-Ferry. An amplitude equation for the non-linear vibration
of viscoelastically damped sandwich beams. Journal of Sound and Vibration, 271 : 789-
813, 2004.
[26]. M. Bilasse, E.M.Daya, L.Azrar. Linear and nonlinear vibrations analysis of viscoelastic
sandwich beams. Journal of Sound and Vibration, 329 : 4950-4969, 2010.
[27]. S. Yi, M.F. Ahmad, H.H. Hilton. Dynamic responses of plates with viscoelastic free layer
damping treatment. Journal of Vibration and Acoustics.; 118:362-367, 1996
[28]. B.C. Nakra. Vibration control in machines and structures using viscoelastic damping.
Journal of Sound and Vibration, 211: 449-66, 1998.
[29]. Y.C. He, W.L. Chen, S.W. Sun, Lina Hao. Finite element research on damping of
viscoelastic free layer damping sheet. Applied Mechanics and Materials, 472: 56-61, 2014.
[30]. F. Cortes, M.J. Elejabarrieta. Structural vibration of flexural beams with thick
unconstrained layer damping. International Journal of Solids and Structures, 45: 5805-13,
2008.
[31]. C. Cai, H. Zheng, G.R. Liu. Vibration analysis of a beam with PCLD Patch. Applied
Acoustics, 65:1057-1076, 2004.
[32]. H. Arvin, M. Sadighi, A.R. Ohadi. A numerical study of free and forced vibration of
composite sandwich beam with viscoelastic core. Composite Structures, 92:996-1008,
2010.
[33]. L.H. Chen, S.C. Huang. Vibration attenuation of a cylindrical shell with constrained layer
damping strips treatment. Computers Structure,79:1355-62, 2001.
189
Références
[34]. B. Khalfi, A. Ross. Influence of partial constrained layer damping on the bending wave
propagation in an impacted viscoelastic sandwich. International Journal of Solids and
Structures, 50: 4133-44, 2013.
[35]. B. Khalfi, A. Ross. Transient response of a plate with partial constrained viscoelastic layer
damping. International Journal of Mechanical Sciences, 68: 304-12, 2013.
[36]. H. Zheng, C. Cai, X.M. Tan. Optimization of partial constrained layer damping treatment
for vibrational energy minimization of vibrating beams. Computers Structures, 82:2493-
507, 2004.
[37]. A.M. Baz. Optimization of energy dissipation characteristics of active constrained layer
damping. Smart Materials and Structures, 6: 360, 1997.
[38]. S. Tliba, H. Abou-Kandil, C. Prieur. Active vibration damping of a smart flexible structure
using piezoelectric transducers: H_infinity design and experimental results. IFAC Proc.,
38: 760, 2005.
[39]. T. Owen, R. Massimo, B. Amr. Attenuation of wave propagation in fluid-loaded shells with
periodic shunted piezoelectric rings. Smart Materials and Structures, 14:594-604, 2005.
[40]. K. El Harti, M. Rahmoune, M. Sanbi, R. Saadani, M. Bentaleb, M. Rahmoune. Finite
element model of vibration control for an exponential functionally graded timoshenko
beam with distributed piezoelectric sensor/actuator. Actuators, 8: 19, 2019.
[41]. S. Herold, D. Mayer. Adaptive piezoelectric absorber for active vibration control.
Actuators, 5: 7, 2016.
[42]. A. Khan, H.S. Kim. Active vibration control of a piezo-bonded laminated composite in the
presence of sensor partial debonding and structural delaminations. Sensors, 19 : 540, 2019.
[43]. S.C. Her, C.S. Lin. Vibration analysis of composite laminate plate excited by piezoelectric
actuators. Sensors, 13: 2997:3013, 2013.
[44]. B. Zhong, B. Deng, H. Zhao. Simulation Model and Method for Active Torsional Vibration
Control of an HEV. Applied Science, 9 : 34, 2019.
[45]. C. Cai, H. Zheng, K. Chung, Z. J Zhang. Vibration analysis of a beam with an active
constraining layer damping patch. Smart Materials and Structures, 15 (2006) 147–156.
[46]. K. Byungjun, R. Emiliano, W. Timothy, M. Brian. Vibration control of beams and plates
with hybrid active-passive constrained layer damping treatments. In Proceedings of the
International Congress of Sound & Vibration, Athens, Greece, 10-14 July, 2016.
[47]. N. Kumar, S.P. Singh, Vibration and damping characteristics of beams with active
constrained layer treatments under parametric variations. Materials & Design, 30 : 4162-
4174, 2009.
[48]. M.C. Ray, J.N. Reddy. Active damping of laminated cylindrical shells conveying fluid
using 1-3 piezoelectric composites. Composite Structures, 98: 26-271, 2013.
[49]. S.C. Kattimani, M.C. Ray. Vibration control of multiferroic fibrous composite plates using
active constrained layer damping. Mechanical Systems and Signal Processing, 106: 334-
354, 2018.
190
Références
[50]. R.S. Kumar, M.C. Ray. Active control of geometrically nonlinear vibrations of doubly
curved smart sandwich shells using 1–3 piezoelectric composites. Composite Structures,
105 : 173–187, 2013.
[51]. J. Landier. Modélisation et étude expérimentale des propretés amortissantes des tôles
sandwich. Thèse de doctorat, Université de Metz, 1993.
[52]. G. Kergourlay. Mesure et prédiction vibro-acoustique de structures viscoélastiques
Application à une enceinte acoustique. Thèse de doctorat, Ecole Centrale Paris, 2004.
[53]. J.L. Halary. Traitement pratique des résultats des essais de viscoélasticité dans
l’introduction à la mécanique des polymères. C. G'sell, J. M. Haudin, Institut National
Polytechnique de Lorrain, 1998.
[54]. X.Q. Zhou, D.Y. Yu, X.Y. Shao, S.Q. Zhang, S. Wang, Research and applications of
viscoelastic vibration damping materials: A review, Composite Structures, 136: 460-480,
(2015)
[55]. M.L. Williams, R.F. Landel, J.D. Ferry. The temperature dependence of relaxation
mechanisms in amorphous polymers and other glass forming liquids. Journal of American
Chemical Society, 77: 3701–3707, 1955.
[56]. E. H. Love. The mathematical theory of elasticity, 4th edition. Cambridge University Press,
1927.
[57]. S.P. Timoshenko. On the transverse vibration of bars of uniform cross-section,
Philosophical Magazine Series 6, 43: 125-131, 1922.
[58]. R. D. Mindlin. Influence of rotatory inertia effect and shear deformation in flexural motions
of isotropic elastic plates. Journal of Applied Mechanics, 18: 1031-1036, 1951.
[59]. J. N. Reddy. A simple higher-order theory of laminated composite plate. Journal of Applied
Mechanics, 51: 745-752, 1984.
[60]. M. Touratier. An efficient standard plate theory. International Journal of Engineering
Science, 29: 901-916, 1991.
[61]. Z. Huang, Z. Qin, F.Chu. Vibration and damping characteristics of sandwich plates with
viscoelastic core. Journal of Vibration and Control, 22(7): 1876-1888, 2016.
[62]. M. Trindade. Contrôle hybride actif–passif des vibrations de structures par des matériaux
piézoélectriques et viscoélastiques : poutres sandwich/multicouches intelligentes. Thèse de
doctorat, 2000.
[63]. J. N. Reddy. A simple higher-order theory of laminated composite plate. Journal of Applied
Mechanics, 51: 745–752, 1984.
[64]. E. Reissner. On certain mixed variational theorem and a proposed application. International
Journal for Numerical Methods in Engineering, 20: 366-1368, 1984.
[65]. J.G. Ren. A new theory of laminated plate. Composite Science and Technology, 26: 225-
239, 1986.
191
Références
[66]. J.G. Ren, E. Hinton. The finite element analysis of homogeneous and laminated composite
plates using a simple higher-order theory. Communications in Applied Numerical Methods,
2: 217-228, 1986.
[67]. D.H. Robbins, J.N. Reddy. Modeling of thick composite using a layer wise laminate theory.
International Journal for Numerical Methods in Engineering, 36: 655-677, 1993.
[68]. W.L. Yin. Free-edge effects in anisotropic laminated under extension, bending and
twisting, Part I: A stress function Based variational approach. Journal of Applied
Mechanics, 61: 410, 1994.
[69]. W.L. Yin. Interlaminar stress analysis of composite laminates using a sublaminate/ layer
model. International Journal of Solids and Structures, vol. 31, no. 11, pages 1549-1564,
1994.
[70]. H. Hu, S. Belouettar, M. Potier-Ferry, E.M. Daya. Review and assessment of various
theories for modeling sandwich composites. Composite Structures, 84(3) :282–292, 2008.
[71]. H. Murakami. Laminated composite plate theory with improved in-plane responses. Journal
of Applied Mechanics, 53 :661–666, 1986.
[72]. E. Carrera. On the use of murakami’s zig-zag function in the modeling of layered plates
and shells. Computers and Structures, 82: 541–554, 2004.
[73]. S. Brischetto, E. Carrera, L. Demasi. Improved bending analysis of sandwich plates using
a zig-zag function. Composites Structures, 89 :408–415, 2009.
[74]. C. Ossadzow, M. Touratier. An improved shear-membrane theory for multilayered shells.
Composite Structures, 52: 85-95, 2001.
[75]. C. Ossadzow. Modélisation de coques composites multicouches. PhD thesis, Ecole
Normale Superieure de Cachan, 1996.
[76]. M. Karama, B. Abou Harb, S. Mistou, S. Caperaa. Bending, buckling and free vibration of
laminated composite with a transverse shear stress continuity model. Composite Part B, 29
(3) : 223-234, 1998.
[77]. K.S. Afaq. Développement d'un nouveau modèle pour les structures composites
"multicouches et sandwichs" avec prise en compte du cisaillement transverse et des effets
de bord. Thèse de doctorat, Université Toulouse III - Paul Sabatier, 2003.
[78]. K.S. Afaq, M. Karama, S. Mistou. Un nouveau modèle raffiné pour les structures
multicouches. In Comptes-rendus des 13 èmes Journées Nationales sur les Composites,
pages 289-292. Strasbourg, March 2003.
[79]. Y. Khadri, S. Tekili, E. M. Daya, A. Daouadji, B.Merzoug. Effects of rail joints and train’s
critical speed on the dynamic behavior of bridges, Mechanika, 19(1): 46-52, 2013.
[80]. Y. Khadri, S. Tekili, E. M. Daya, A. Daouadji, B. Merzoug. Dynamic analysis of train–
bridge system and riding comfort of trains with rail irregularities, Journal of Mechanical
Science and Technology, 27(4): 951-962, (2013).
[81]. S. Tekili, Y. Khadri, B. Merzoug, E. M. Daya, A. Daouadji. Free and forced vibration of
beams strengthened by composite coats subjected to moving loads. Mechanics of
Composite Materials, 52: 789-798, 2017.
192
Références
[82]. M. Abu-Hilal, H. S. Zibdeh, Vibration analysis of beams with general boundary conditions
traversed by a moving force. Journal of Sound and Vibration, 229(2): 377-388 (2000).
[83]. F. Giuntaa, G. Muscolinoa, A. Sofic and I. Elishakoffb. Dynamic analysis of Bernoulli-
Euler beams with interval uncertainties under moving loads. X International Conference on
Structural Dynamics, EURODYN 2017.
[84]. Y.Chen, Yiming Fu, Jun Zhong and Chang Tao. Nonlinear dynamic responses of fiber
metal beam subjected to moving harmonic loads resting on tensionless elastic foundation.
Composites Part B, 131: 253-259, 2017.
[85]. C. Tao, Y.M. Fu, H.L Dai, Nonlinear dynamic analysis of fiber metal laminated beams
subjected to moving loads in thermal environment, Composite Structures, 140:410-416,
2016,
[86]. V. Sarvestan, H.R. Mirdamadi, M. Ghayour, Vibration analysis of cracked Timoshenko
beam under moving load with constant velocity and acceleration by spectral finite element
method, International Journal of Mechanical Sciences, 122: 318-330, 2017.
[87]. Z. Kiral and B. G. Kiral. Dynamic analysis of a symmetric laminated composite beam
subjected to a moving load with constant velocity. Journal of Reinforced Plastics and
Composites, 27(1): 19-32, 2008.
[88]. F.G. Yuan, R.E. Miller. A new finite element for laminated composite beams. Computers
Structures, 31: 737-745, 1989.
[89]. V. Kahya. Dynamic analysis of laminated composite beams under moving loads using finite
element method. Nuclear Engineering and Design, 243: 41-48, 2011.
[90]. L. Rayleigh. The Theory of Sound, 2nd ed., Reprinted by Dover, New York, (1945).
[91]. C.E. Wallace. Radiation resistance of a rectangular panel, Journal of the Acoustical Society
of America, 51: 946-952, 1972.
[92]. N.S. Lomas, S.I. Hayek. Vibration and acoustic radiation of elastically supported
rectangular plates, Journal of Sound and Vibration, 2: 1-25, 1977.
[93]. L. Cremer, M. Heckl. Structure-Borne Sound, second ed., Springer- Verlag, New York,
1987.
[94]. E.G. Williams. A series expansion of the acoustic power radiated from planar sources.
Journal of the Acoustical Society of America, 73: 1520-1524, 1983.
[95]. M.R. Lee, R. Singh. Analytical formulations for annular disk sound radiation using
structural modes. Journal of the Acoustical Society of America, 95(6): 3311–3323, 1994.
[96]. D. Liu. X. Li. An Overall View of Laminate Theories Based on Displacement Hypothesis,
Journal of Composite Materials, 30(14): 1539-1561, 1996.
[97]. J.S. Tao, G.R. Liu, K.Y. Lam. Sound radiation of a thin infinite plate in light and heavy
fluids subject to multi point excitation. Applied Acoustics, 62: 573–587, 2001.
[98]. S. Lin. Study on the radiation acoustic field of rectangular radiators in flexural vibration.
Journal of Sound and Vibration, 254(3): 469-479, 2002.
[99]. L. Zhou, H. Zheng, K.C. Hung. Sound radiation from a thin infinite plate in contact with
a layered inhomogeneous fluid, Applied Acoustics, 63: 1177-1192, 2002.
193
Références
194
Références
[115]. S. Assaf. Modélisation par la méthode des éléments finis du comportement vibratoire des
poutres et plaques sandwiches : Métal matériaux viscoélastiques. Thèse de doctorat,
Université de Technologie de Compiègne, 1991.
[116]. S. Assaf, M. Guerich, P. Cuvelier. Vibration and acoustic response of damped sandwich
plates immersed in a light or heavy fluid. Computer and Structures, 88:870-878, 2010.
[117]. M. Abdennadher, S. Fenina, L. Hammami, M. Haddar. Vibro-acoustic analysis of a double
sandwich panels system. International Journal of Engineering Systems Modelling and
Simulation. 6: 3-9, 2005.
[118]. S. Das, S. Chandra, A.G. Niyogi. Vibro-acoustic analysis inside cross ply laminated
composite cavity. Procedia Engineering. Procedia Engineering, 144: 1389-96, 2016.
[119]. S. Narayanan, R.L Shanbhag. Sound transmission through a damped sandwich panel.
Journal of Sound and Vibration, 80(3):315-27, 1982.
[120]. N. Atalla, F. Sgard. Finite Element and Boundary Methods in Structural Acoustics and
Vibration. CRC Press, Taylor & Francis Group, 2015.
[121]. B.M. Irons. Role of part inversion in fluid-structure problems with mixed variables. AIAA
Journal, 7: 568, 1970.
[122]. M.C.M. Wright. In lecture notes on the mathematics of acoustics. 1ère Edition. London,
Imperial College Press. 2005.
[123]. M. Latifi, M. Kharazi, H.R. Ovesy, Effect of Integral Viscoelastic Core on the Nonlinear
Dynamic Behaviour of Composite Sandwich Beams with Rectangular Cross Sections,
International Journal of Mechanical Sciences, 123: 141-150, 2017.
[124]. E. Carrera, F. A. Fazzolari, M. Cinefra. Thermal Stress Analysis of Beams Plates and
Shells, Coupled and uncoupled variational formulations, Elsevier, 2017.
[125]. M. Bilasse, Modélisation numérique des vibrations linéaires et non linéaires des structures
sandwich à âme viscoélastique. Thèse de doctorat, Université Paul Verlaine Metz, 2010.
[126]. L. Xiangyang, K. Yu. Vibration and acoustic responses of composite and sandwich panels
under thermal environment, Composite Structures, 131: 1040-1049, 2015.
[127]. M. Haberman. Design of high loss viscoelastic composites through micromechanical
modeling and decision based material by design, PhD Thesis, Woodruff School of
Mechanical Engineering, 2007.
[128]. M. Trindade, A. Benjeddou, R. Ohayon. Modeling of frequency dependent viscoelastic
materials for active–passive vibration damping. Journal of Vibration and Acoustics,122:
169-174, 2000.
[129]. S. Adhikari, Damping models for structural vibration, PhD thesis, Cambridge University,
UK, (2000).
[130]. B. Cochelin, N. Damil, M. Potier-Ferry. Méthode Asymptotique Numérique. Hermès
Science publications, France, 2007
[131]. N. Damil. An iterative method based upon Padé approximants. Communications in
Numerical Methods in Engineering, 15: 701-708, 1999.
195
Références
[132]. S. Tekili, Y. Khadri, Y. Karmi. Dynamic analysis of sandwich beam with viscoelastic core
under moving loads. MECHANIKA. 26 (4): 325330, 2020.
[133]. Y. Karmi, SabihaTekili, Youcef Khadri. Free vibrations analysis of sandwich beams with
viscoelastic core, Rev. Sci. Technol., Synthèse 39, 2019.
[134]. Y. Karmi, Y. Khadri, S. Tekili, A. Daouadji, E.M. Daya, Effect of Mechanical and
geometrical properties on dynamic behavior of asymmetrical composite sandwich beam
with viscoelastic core. UPB Scientific Bulletin, Series D: Mechanical Engineering, 2020.
[135]. Y. Karmi, Y. Khadri, S. Tekili, A. Daouadji, E.M. Daya, Dynamic analysis of composite
sandwich beams with frequency dependent viscoelastic core under the action of moving
load. Mechanics of Composite Materials. 25(6), 2020.
196