L'objectif de cette thèse est de développer des modèles d'homogénéisation analytiques de panneaux sandwichs en nid d'abeilles. A la différence des méthodes classiques, l'effet des peaux est pris en compte, conduisant à des propriétés mécaniques très différentes. Dans les cas des tractions, flexions, cisaillement dans le plan, cisaillements transversaux et torsion, différentes séries de fonctions analytiques sont proposées pour prendre en compte la redistribution des contraintes entre les parois du nid d'abeilles. Nous avons étudié l'influence de la hauteur du nid d'abeilles sur les propriétés élastiques. Les courbes des modules obtenues avec le modèle proposé sont bien bornées par les valeurs obtenues avec la théorie des poutres. Les contraintes d'interface sont également étudiées afin de comparer avec les modèles existant pour le problème de traction. De nombreux calculs numériques ont été réalisés avec nos H-modèles pour les problèmes de tractions, de flexions, de traction-flexion couplés, de cisaillement dans le plan, de cisaillement transversal et de torsion. De très bon accords ont été obtenus entre les résultats issus des H-modèles et ceux issus des calculs en éléments finis de coques en maillant complètement les panneaux sandwichs. Nos H-modèles ont été appliquées aux calculs de grandes plaques sandwichs industrielles en nid d'abeilles. La comparaison desrésultats entre les H-modèles et les calculs en éléments finis de coques du logiciel Abaqus sont en très bon accord. / The aim of this thesis is to develop an analytical homogenization model for the honeycomb core sandwich panels. Unlike conventional methods, the skin effects are taken into account, leading to a very different mechanical properties. In the cases of extensions, bendings, in-plane shear, transverse shears andtorsion, different analytical function series are proposed to consider the stress redistribution between the honeycomb walls. We have studied the influence of the height of the core on its homogenized properties. The moduli curves obtained by the present H-models are well bounded by the moduli values obtained by the beam theory. The interface stresses are also studied to compare with existing models for stretching problem. Many numerical computations with our H-models have been done for the problems of stretching, bending, in-plan and transverse shearing, as well as torsion. Very good agreement has been achieved between the results of the H-models and the results obtained by finite element simulations by completely meshing thesandwich panel with shell elements. Our H-models have been applied to the computations of industrial large sandwich panels with honeycomb core. The comparison of the results between the H-models and the simulations with Abaqus shell elements are in very good agreement.

La compréhension et l'amélioration de l'environnement vibratoire des charges utiles demande la mise au point de démarches prédictives maîtrisées qui permettent de comprendre les phénomènes de transmission des ondes de chocs d'origine pyrotechnique dans le lanceur Ariane5. Plus particulièrement, la maîtrise du comportement transitoire des coques sandwichs en nid d'abeilles, principaux constituants de l'Adaptateur de Charges Utiles – structure porteuse des satellites, est nécessaire pour prédire les vibrations au pied des équipements électroniques des satellites et des lanceurs. Cette problématique présente un caractère multi-échelle tant d'un point de vue temporel (charge mobile supersonique, temps d'analyse) que spatial (dimensions des structures du lanceur, taille des cellules en nid d'abeilles, longueurs d'ondes liées aux hautes fréquences). Celui-ci a été traité dans cette thèse en s'appuyant d'une part, sur une qualification à la fois analytique et numérique des modèles classiques homogénéisés des plaques sandwichs en nid d'abeilles pour la gamme de fréquence mise en jeu et d'autre part, sur une application des stratégies de remaillage adaptatif pour la propagation des ondes développées dans le cadre de la méthode de Galerkin espace-temps discontinue en temps. Deux catégories de modèles de plaques épaisses ont été ainsi construites dans le but d'enrichir la cinématique classique de plaques épaisses de Mindlin-Reissner qui s'est avérée être insuffisante pour correctement représenter le comportement dynamique hors-plan des plaques sandwich en nid d'abeilles. Ainsi ont été analysés les modèles dits monocouches basés sur un enrichissement de la cinématique par ajout de degrés de liberté dans l'épaisseur, et les modèles multicouches composés d'une superposition de trois plaques avec une homogénéisation séparée des matériaux. Il a été montré que ces deux sortes de modèles améliorent la description des phénomènes de hautes fréquences, notamment ceux de flexion et de cisaillement transverse qui sont plus délicats à retranscrire. Toutes les études numériques ont été effectuées avec un code éléments finis qui emploie des solveurs adaptatifs dynamiques basés sur la méthode de Galerkin espace-temps discontinue en temps. Cette méthode d'intégration en temps introduit un amortissement numérique dépendant du pas de temps et qui peut interférer avec un amortissement physique susceptible d'être introduit dans un modèle numérique et conduire au final à un amortissement total différent de celui qui est attendu. Cette interaction a été analysée et mise en évidence dans ce travail à travers l'introduction de l'amortissement de Rayleigh dans les modèles de propagation de chocs. Les outils et les modèles de propagation ainsi développés ont été validés sur plusieurs structures académiques et industrielles. Des comparaisons avec des données expérimentales sur des structures industrielles de grande taille, plus particulièrement sur un Adaptateur de Charges Utiles d'Ariane5, sont effectuées et soulignent la cohérence de notre approche ainsi que la fiabilité et l'efficacité des modèles de propagation proposés.

Dans le cadre linéaire, les modèles de Mindlin-Reissner pour les plaques épaisses et de Naghdi pour les coques épaisses sont les plus utilisés. Il est connu que la discrétisation par éléments finis de ces modèles conduit à un phénomène de verrouillage numérique quand l’épaisseur tend vers zéro. Il s’agit du verrouillage en cisaillement dans le cas des plaques et du verrouillage en cisaillement et en membrane dans le cas des coques. Il existe quelques éléments finis qui permettent d’éviter ces difficultés ou du moins de les réduire. L’élément DKMQ pour les plaques et sa version DKMQ24 pour les coques, sont des éléments de bas ordre, basés sur une formulation mixte, qui ont été proposés il y a quelques années afin d’éviter ces phénomènes de verrouillage. Dans cette thèse, on s’est attaché à évaluer numériquement les performances de ces éléments. Outre les cas tests classiques, on s’est focalisé sur l’analyse de la condition inf-sup discrète pour l’élément DKMQ. Nous avons étudié également le test de la s-norme proposé par Bathe, pour l’élément DKMQ24. Enfin, nous avons effectué une analyse d’erreur a posteriori pour les éléments DKMQ et DKMQ24, en utilisant l’estimateur d’erreur Z2 (dû à Zienkiewicz et Zhu), associé aux techniques de recouvrement de la moyenne, de projection ou encore SPR. Les résultats obtenus ont permis de quantifier les performances de ces deux éléments finis pour les problèmes de verrouillage, et d’en dégager les limites. Deux applications importantes de ces éléments DKMQ et DKMQ24 ont été ensuite présentées, la première concerne la simulation des poutres à parois minces à section ouverte et la seconde le calcul des plaques composites. / In the linear case, the Mindlin-Reissner model for thick plates and the Naghdi model for thick shells are commonly used. The finite element discretization of these models leads to numerical locking phenomenon when the thickness approaches zero : shear locking for plates and both shear and membrane locking for shells. There are some finite elements that could reduce or even eliminate this phenomenon. DKMQ element for plates or DKMQ24 element for shells, are low-order elements, based on a mixed formulation, introduced a few years ago to prevent the numerical locking phenomenon. In this thesis, we concentrated on numerical evaluation of the performance of these elements. Besides the classical benchmark tests, we also focused on the analysis of discrete inf-sup condition for DKMQ element. We studied the s-norm test proposed by Bathe for DKMQ24 element. Finally, we performed a posteriori error estimation for DKMQ and DKMQ24 elements, using the error estimator Z2 (proposed by Zienkiewicz and Zhu), associated with the averaging, projection or SPR recovery methods. The results obtained have enabled us to quantify the performance of these two finite elements for locking problems, and to identify their limits. Two important applications of these elements DKMQ and DKMQ24 were then presented ; the first one concerns thin-walled beams with open cross-section and the second one composite plates.

La compréhension et l’amélioration de l’environnement vibratoire des charges utiles demande la mise au point de démarches prédictives maîtrisées qui permettent de comprendre les phénomènes de transmission des ondes de chocs d’origine pyrotechnique dans le lanceur Ariane5. Plus particulièrement, la maîtrise du comportement transitoire des coques sandwichs en nid d’abeilles, principaux constituants de l’Adaptateur de Charges Utiles – structure porteuse des satellites, est nécessaire pour prédire les vibrations au pied des équipements électroniques des satellites et des lanceurs. Cette problématique présente un caractère multi-échelle tant d’un point de vue temporel (charge mobile supersonique, temps d’analyse) que spatial (dimensions des structures du lanceur, taille des cellules en nid d’abeilles, longueurs d’ondes liées aux hautes fréquences). Celui-ci a été traité dans cette thèse en s’appuyant d’une part, sur une qualification à la fois analytique et numérique des modèles classiques homogénéisés des plaques sandwichs en nid d’abeilles pour la gamme de fréquence mise en jeu et d’autre part, sur une application des stratégies de remaillage adaptatif pour la propagation des ondes développées dans le cadre de la méthode de Galerkin espace-temps discontinue en temps. Deux catégories de modèles de plaques épaisses ont été ainsi construites dans le but d’enrichir la cinématique classique de plaques épaisses de Mindlin-Reissner qui s’est avérée être insuffisante pour correctement représenter le comportement dynamique hors-plan des plaques sandwich en nid d’abeilles. Ainsi ont été analysés les modèles dits monocouches basés sur un enrichissement de la cinématique par ajout de degrés de liberté dans l’épaisseur, et les modèles multicouches composés d’une superposition de trois plaques avec une homogénéisation séparée des matériaux. Il a été montré que ces deux sortes de modèles améliorent la description des phénomènes de hautes fréquences, notamment ceux de flexion et de cisaillement transverse qui sont plus délicats à retranscrire. Toutes les études numériques ont été effectuées avec un code éléments finis qui emploie des solveurs adaptatifs dynamiques basés sur la méthode de Galerkin espace-temps discontinue en temps. Cette méthode d’intégration en temps introduit un amortissement numérique dépendant du pas de temps et qui peut interférer avec un amortissement physique susceptible d’être introduit dans un modèle numérique et conduire au final à un amortissement total différent de celui qui est attendu. Cette interaction a été analysée et mise en évidence dans ce travail à travers l’introduction de l’amortissement de Rayleigh dans les modèles de propagation de chocs. Les outils et les modèles de propagation ainsi développés ont été validés sur plusieurs structures académiques et industrielles. Des comparaisons avec des données expérimentales sur des structures industrielles de grande taille, plus particulièrement sur un Adaptateur de Charges Utiles d’Ariane5, sont effectuées et soulignent la cohérence de notre approche ainsi que la fiabilité et l’efficacité des modèles de propagation proposés. / Reliable and efficient numerical models for the pyrotechnic shock wave propagation in structures of the Ariane5 launcher are necessary for a good understanding and a predictive analysis of the payload vibration environment. More precisely, the correct modeling of the dynamic behaviour of the honeycomb sandwich shells, the main material composing the payload adaptor, is essential to control the vibration environment of the payload and the embarked electronic equipments and so to prevent them from damages caused by the shock wave propagation. The topic is obviously a multi-scale problem from both temporal and spatial points of view : short time intervals imposed by supersonic moving loads vs. large total time interval that the slowest waves need to travel throughout the adaptor ; very short wavelengths of high frequency waves, and very small size of the honeycomb cells vs. large structure dimensions. To take into account all involved space-time scales in a reliable and efficient way, the herein study is based both on the analytical and numerical qualification of the classical homogenized models of honeycomb sandwich shells for the frequency range introduced by the pyrotechnic shock wave, and on a dynamic solver based on the well-known space-time discontinuous Galerkin method, allowing the use of adaptive remeshes for the wave propagation. The classical Mindlin-Reissner’s kinematics of thick plates being inefficient to correctly represent the dynamic out-of-plane behaviour of the honeycomb sandwich plates, two kinds of its enrichment are considered : One-layered models based on an enrichment of the kinematics by adding degrees of freedom in the thickness, and multi-layered models composed of a superposition of three plates with separated material homogenisations. It has been shown theoretically and numerically that, both types of enrichment allow more precise descriptions of flexure and transverse shear modes in the high frequency range. However, the multi-layered models give much more promising results, as the important role played by the honeycomb core for the transverse shear behaviour of the whole sandwich is not “smeared” in a one-layered homogenized model. All the numerical studies were conducted with a finite element code which uses a dynamic solverbased on the time discontinuous space-time Galerkin method. The built-in numerical damping of this solver can interfere with a physical damping potentially introduced by the numerical model and results in a global damping totally unexpected. This interaction has been analysed and underlined in this work thanks to the introduction of the Rayleigh damping in the shock wave propagation models. Theoretical and numerical tools and propagating models thus developed have been validated on several academic and industrial structures. Comparison with experimental data on large size industrial structures, especially a real size payload adaptor, is performed and emphasizes the coherence of our approach and the reliability and the efficiency of the proposed propagating models.

La compréhension et l'amélioration de l'environnement vibratoire des charges utiles demande la mise au point de démarches prédictives maîtrisées qui permettent de comprendre les phénomènes de transmission des ondes de chocs d'origine pyrotechnique dans le lanceur Ariane5. Plus particulièrement, la maîtrise du comportement transitoire des coques sandwichs en nid d'abeilles, principaux constituants de l'Adaptateur de Charges Utiles - structure porteuse des satellites, est nécessaire pour prédire les vibrations au pied des équipements électroniques des satellites et des lanceurs. Cette problématique présente un caractère multi-échelle tant d'un point de vue temporel (charge mobile supersonique, temps d'analyse) que spatial (dimensions des structures du lanceur, taille des cellules en nid d'abeilles, longueurs d'ondes liées aux hautes fréquences). Celui-ci a été traité dans cette thèse en s'appuyant d'une part, sur une qualification à la fois analytique et numérique des modèles classiques homogénéisés des plaques sandwichs en nid d'abeilles pour la gamme de fréquence mise en jeu et d'autre part, sur une application des stratégies de remaillage adaptatif pour la propagation des ondes développées dans le cadre de la méthode de Galerkin espace-temps discontinue en temps. Deux catégories de modèles de plaques épaisses ont été ainsi construites dans le but d'enrichir la cinématique classique de plaques épaisses de Mindlin-Reissner qui s'est avérée être insuffisante pour correctement représenter le comportement dynamique hors-plan des plaques sandwich en nid d'abeilles. Ainsi ont été analysés les modèles dits monocouches basés sur un enrichissement de la cinématique par ajout de degrés de liberté dans l'épaisseur, et les modèles multicouches composés d'une superposition de trois plaques avec une homogénéisation séparée des matériaux. Il a été montré que ces deux sortes de modèles améliorent la description des phénomènes de hautes fréquences, notamment ceux de flexion et de cisaillement transverse qui sont plus délicats à retranscrire. Toutes les études numériques ont été effectuées avec un code éléments finis qui emploie des solveurs adaptatifs dynamiques basés sur la méthode de Galerkin espace-temps discontinue en temps. Cette méthode d'intégration en temps introduit un amortissement numérique dépendant du pas de temps et qui peut interférer avec un amortissement physique susceptible d'être introduit dans un modèle numérique et conduire au final à un amortissement total différent de celui qui est attendu. Cette interaction a été analysée et mise en évidence dans ce travail à travers l'introduction de l'amortissement de Rayleigh dans les modèles de propagation de chocs. Les outils et les modèles de propagation ainsi développés ont été validés sur plusieurs structures académiques et industrielles. Des comparaisons avec des données expérimentales sur des structures industrielles de grande taille, plus particulièrement sur un Adaptateur de Charges Utiles d'Ariane5, sont effectuées et soulignent la cohérence de notre approche ainsi que la fiabilité et l'efficacité des modèles de propagation proposés.