Le procédé de soudage par friction et malaxage est un procédé récent qui a été développé au sein de l'institut de soudure britannique "The Welding Institute" au début des années 90. Ce procédé, utilisé généralement en aéronautique, est sans apport de matière et permet de souder principalement des alliages d'aluminium difficilement soudables par les procédés classiques de soudage. Il consiste à malaxer le matériau de base à l'aide d'un outil constitué d'un pion et d'un épaulement frottant sur les faces supérieures des tôles à souder. La modélisation de ce procédé est très complexe puisque ce dernier implique des couplages entre des phénomènes mécaniques, thermiques et métallurgiques. Le malaxage dans le procédé de soudage FSW est difficile à simuler à l'aide de la méthode des éléments finis (en lagrangien) puisque la zone proche de l'outil de soudage est le siège de grandes déformations. Donc le remaillage est nécessaire. Cependant, le remaillage est cher et très difficile pour les problèmes tridimensionnels. Par ailleurs, après un remaillage, il est nécessaire d'interpoler les champs (vitesses, contraintes,...) correspondant à la solution courante, ce qui peut introduire des erreurs supplémentaires dans le calcul (on parle de diffusion numérique). Nous proposons dans ce travail des modèles basés sur la méthode sans maillage dite "Smoothed Particle Hydrodynamics SPH" et la méthode des moindres carrés mobiles (Moving Least Square MLS) pour la simulation de ce procédé. Ces modèles sont formulés dans le cadre lagrangien et utilisent la forme forte des équations aux dérivées partielles. Le premier modèle basé sur SPH considère la zone de soudure comme un fluide non newtonien faiblement compressible et dont la viscosité dépend de la température. Ce modèle est proposé pour la simulation numérique du comportement thermomécanique d'un matériau soudé par le procédé FSW. Dans le deuxième modèle, un algorithme itératif implicite de premier ordre a été proposé, pour simuler le malaxage de la matière dans le cas d'un matériau viscoplastique, en utilisant la méthode MLS et la technique de collocation. Le troisième modèle est un algorithme implicite d'ordre élevé basée sur le couplage de la méthode MLS et la Méthode Asymptotique Numérique MAN. Cet algorithme permet de réduire le temps de calcul par rapport à l'algorithme itératif implicite de premier ordre. La validation de ces trois modèles proposés a été faite par le code industriel Fluent / Friction stir welding is a recent process that has been developed by the British Welding Institute TWI "The Welding Institute" since 1990s. This process, generally used in aerospace, does not need additional material and allows mainly joining plates of aluminum alloys which are difficult to weld by the classical welding processes. It consists in mixing the base material using a tool comprising a pin and a shoulder which heats the plates to be welded by friction. The modeling of this process is very complex since it involves the coupling between mechanical, thermal and metallurgical phenomena. The mixing in welding process FSW is difficult to simulate using finite element method in lagrangian framework since the area near the welding tool is submitted to large deformations. So remeshing procedure is often required. However, remeshing can be very expensive and difficult to perform for three-dimensional problems. Moreover, after remeshing step, it is necessary to interpolate the fields (velocities, constraints ...) corresponding to the current solution, which may lead to additional errors in the calculation (called numerical diffusion). We propose in this work models based on meshless methods called "Smoothed Particle Hydrodynamics SPH" and Moving Least Square method for the simulation of this welding process. These models are formulated in lagrangian framework and use the strong form of partial differential equations. The first model based on SPH considers the welding zone as a weakly compressible non-newtonian fluid and whose viscosity depends on the temperature. This model is proposed for the numerical simulation of thermo-mechanical behavior of a welded material by FSW process. The second model is a first order implicit iterative algorithm proposed to simulate material mixing in the case a visco-plastic behavior using the MLS method and the collocation technique. The third model is a high order implicit algorithm based on the coupling of MLS method and Asymptotic Numerical Method (ANM). This algorithm allows reducing the computation time by comparison with the first order implicit iterative algorithm. The validation of these three proposed models was done by the industrial code Fluent

Lam3/Tec3 est un modèle stationnaire éléments finis 3D de laminage dans lequel le maillage de la tôle est composé d'éléments hexaédriques. Des études passées ont montré que par lui-même, il est mal adapté pour simuler les défauts de planéité des tôles laminées. A l'heure actuelle, les éléments coques sont les mieux adaptés pour prédire le flambage des tôles sous contraintes résiduelles, dont font partie les défauts de planéité. Ainsi, l'objectif du présent travail est de perfectionner Lam3/Tec3 en le couplant itérativement avec un modèle éléments finis de flambage de coques sous contraintes résiduelles. Ce modèle, nommé « MAN », est basé sur la Méthode Asymptotique Numérique pour la résolution des problèmes non linéaires. Comparé aux modèles existants de la littérature, il possède une précision satisfaisante pour la détection de la charge critique de flambage et les modes correspondants; il peut aussi traiter des cas de champs de contraintes résiduelles et modes de flambage plus complexes; il possède enfin une capacité de calcul du post-flambage, ce qui le distingue aussi de ses concurrents. Nous avons comparé les résultats obtenus en mode découplé et en mode couplé itérativement, ainsi qu'avec Lam3/Tec3 additionné d'un modèle simplifié de flambage dû à Counhaye (plus simple et moins coûteux, mais sans capacité de post-flambement). Nous nous sommes particulièrement intéressés à la comparaison avec les mesures en ligne des profils transverses de contrainte longitudinale par "rouleaux de planéité" (la "planéité latente"), ainsi qu'à la question des interactions entre les champs mécaniques sous emprise et post-emprise. Les profils obtenus à la position des rouleaux de planéité sont en bon accord avec les mesures expérimentales, en particulier pour les cas de laminage présentant des défauts. Les résultats montrent l'importance de réaliser un calcul couplé, soit fortement soit itérativement, pour saisir les détails de ces profils de contrainte, capitaux à quelques MPa près pour le calcul des formes post-flambées, c'est-à-dire les défauts de planéité.

Avec l’essor de l’électronique embarquée, les câbles électriques constituentune part importante des pièces automobiles tandis que l’espace à bord n’a cessé de diminuer. Leur flexibilité requiert la prédiction de leur déformation durant leur montage afin d’éviter le contact avec d’autres pièces du véhicule et leur endommagement. Les outils actuels ne permettent pas une prédiction assez réaliste et précise de leur comportement, nécessaire dans un volume de travail très restreint. Les étapes de montage sont donc validées via la réalisation de maquettes réelles coûteuses. Cette thèsea pour but d’améliorer la simulation numérique de ces pièces souples. Nous proposonsici un code de simulation 3D basé sur un modèle de poutre géométriquement exact résolu par la méthode des éléments finis. Son originalité tient dans le couplage des quaternions pour modéliser les rotations 3D et de la méthode asymptotique numérique pour la continuation du système non linéaire qui lui confère une grande robustesse. Un banc d’essai permettant l’identification des paramètres homogénéisés nécessaires au modèle numérique et sa validation par comparaison de la géométrie finale et du chemin d’équilibre est présenté. Combinés à des développements analytiques sur les modèles de poutres avec cisaillement, les essais mènent à une évaluation critique du modèle deTimoshenko 3D pour la représentation des torons de câbles. / With on-board electronics expansion, electrical cables are an essential partof automotive pieces and the space on board has plummeted. Their flexibility requires to predict their deformation during vehicle assembly in order to avoid the contact with other pieces and damaging. Current numerical tools do not allow a realistic and accurate prediction, which is necessary in the obstructed car space. Assembly steps thus are validated on costly physical mock-ups. This thesis aims at improving numerical simulation of these flexible pieces. We herein propose a 3D algorithm based on a geometrically exact beam model solved by the finite element method. This work’s originality stands in coupling quaternions as rotational parameters and the asymptotic numerical method as nonlinear solver which results in a very robust algorithm. A test bench designed to identify the homogenized beam parameters of the numerical model and to validate it by offering a comparison on the final geometry and the equilibrium path is presented. Analytical developments on shear beams and the results of these experimental tests lead to a critical evaluation of the 3D Timoshenko model for representing stranded cables.

Les transformations qui conservent l'ensemble des solutions des équations de Navier-Stokes (NS) sont appelées les symétries de NS. Elles forment un groupe de Lie dénommé groupe de symétrie de NS. Ce groupe jouent un rôle important dans la description de la physique des équations (loi de conservation, loi de paroi, ...). Ainsi, les modèles de turbulence devraient être invariant sous l'action de ce groupe. Dans la première partie de la thèse, on effectue alors une analyse de quelques modèles de sous-maille courants sous l'angle des symétries, puis, on construit une classe de modèles de sous-maille qui, d'une part, respectent le groupe de symétrie de NS et, d'autre part, sont conformes au second principe de la thermodynamique. Un modèle très simple de la classe est alors testé et validé numériquement. L'analyse et la construction de modèles sont également étendues au cas de la convection thermique. Dans la seconde partie de la thèse, on explore la possibilité d'intégrer la LES (simulation des grandes échelles) dans un algorithme de la famille MAN (méthode asymptotique numérique). La MAN est une technique numérique de perturbation, qui consiste à calculer la solution sous forme d'une série entière. Dans un premier temps, on construit et on teste un algorithme associant la MAN et la LES, avec l'aide d'une technique d'homotopie. Face aux limites de ce premier algorithme, on étudie dans un second temps l'utilisation d'un autre algorithme où on effectue un développement en série temporelle. Pour augmenter le domaine de validité de la série obtenue, ou bien pour calculer une solution analytique à partir de la série lorsque celle-ci diverge, on propose d'effectuer la méthode de resommation de Borel-Laplace. Dans les exemples numériques, on applique cette méthode à des modèles réduits issus des équations de Navier-Stokes.

Cette thèse traite de la manière dont le musicien contrôle les régimes d'oscillation d'un instrument de la famille des cuivres. L'’association d’'un instrument et d'un musicien forme un système dynamique complexe. Le couplage non linéaire entre les lèvres et la colonne d’'air peut mener à un régime stationnaire, des auto-oscillations périodiques ou quasi-périodiques. De nombreux paramètres de contrôle permettent à un musicien de choisir le régime d’oscillation.Deux approches sont utilisées pour comprendre l’'émergence des régimes de l’'instrument et leur contrôle par le musicien. D'’une part, on étudie un modèle physique de cuivre, consistant en une valve à 1 degré de liberté couplée non linéairement à une expression modale de l’'impédance d’'entrée de l’'instrument. Les caractéristiques des régimes de ce modèle sont étudiés par analyse de stabilité linéaire, simulation numérique et équilibrage harmonique. L’'objectif est d'’estimer la capacité de ce modèle à reproduire le comportement de l’'instrument. Un outil expérimental est également développé pour mesurer simultanément des paramètres de contrôle employés par les musiciens et des variables de jeu.Les résultats améliorent la compréhension des régimes d’oscillation des cuivres. Les mesures expérimentales donnent les plages de variation des paramètres de jeu. Les résultats des méthodes d'’analyse se complètent, donnant une vision globale de l’'apparition des régimes oscillants du modèle et des informations détaillées sur certains points de fonctionnement. Certaines situations de jeu particulières sont étudiées: note pédale du trombone, ’influence de l’'utilisation d’'une sourdine et nature des sons multiphoniques. / This thesis focuses on how a musician controls and chooses the oscillation regimes of brass instruments.The association of an instrument and a musician is a complex dynamical system. The nonlinear coupling between the lips and the resonator can lead to a stationary regime, periodic or quasi-periodic self-sustained oscillations. Through several control parameters, a musician can select the desired regime.Two approaches are developed to understand the emergence of the instrument oscillation regimes and how a musician controls it. A physical model of brass instrument is studied. It consists in a one-degree-of-freedom lip valve nonlinearly coupled to a modal fit of the input impedance of the instrument. The features of several oscillation regimes of this model are studied through linear stability analysis, numerical simulation and harmonic balance. Goals are the assessment of the ability of this model to reproduce the behavior of the modelled instrument. An experimental setup is also developed to measure simultaneously playing parameters used by a musician and instrument variables.Results provide a better understanding of the behavior of brass instrument oscillation regimes. Measurements provide ranges of the control parameters. Results of the analysis methods of the model are complementary. A global vision of the occurrence of certain oscillation regimes is provided, along with more detailed information for some operating points. Some specific playing situations are examined in detail, such as the trombone pedal note, the influence of a passive or active mute on the latter, and the nature of multiphonic sounds.

Le plissement est généralement observé dans les structures minces ayant un comportement de type membrane. Ces structures minces ne supportent pas d'effort de flexion et sont donc sollicitées en traction. Dans cette thèse, nous avons développé une technique de réduction de modèle pour la modélisation du plissement des membranes minces. Cette technique, basée sur les séries de Fourier à double échelle, permet de déduire d'un modèle complet de membrane, un modèle réduit capable de prendre en compte les instabilités globales et locales. Les valeurs critiques de charge et longueur d'onde sont déterminées analytiquement puis numériquement. Des exemples numériques nous ont permis de valider le modèle numérique par rapport au modèle analytique. Les modèles numériques étudiés prennent en compte le modèle complet et le modèle réduit de la membrane. Le modèle complet est simulé dans Abaqus et résolu numériquement à l'aide de la méthode de la longueur d'arc et le modèle réduit est implémenté dans Matlab et résolu numériquement à l'aide de la méthode asymptotique numérique. Nous avons étudié le comportement de la membrane sous sollicitation mécanique, thermique et thermo-mécanique. Les résultats obtenus montrent que le modèle réduit est capable de se substituer au modèle complet dans la détermination des contraintes critiques et longueurs d'onde correspondantes. Le gain en temps de calcul obtenu est important, ceci grâce à la très faible densité de maillage requis par le modèle réduit. Le modèle réduit est très sensible aux conditions aux bords de la membrane et requiert d'avoir une longueur d'onde des plis quasiment constante dans la largeur de la membrane / Wrinkling is an instability phenomenon generally observed in thin structures with membrane's behavior. Those thin structures have no rigidity to flexion and are therefore used in traction. In this thesis, we developed a reduction model's technique for the modeling of wrinkling phenomenon in thin membranes. This technique, based on the double scale Fourier series, allow us to deduce from a full membrane model, a reduced membrane model that is able to take into account the global and local instability of the structure. The critical load and critical wavelength are determined analytically on one side, then numerically on the other side. Numerical exemples are conducted to validate the numerical model towards the analytical one. Numerical models studied take into account both full and reduce membrane models. The full model is simulated in Abaqus and solved numerically using the arc length method and the reduced model is implemented in Matlab and solved numerically using the asymptotic numerical method. We studied the membrane behavior under mechanical, thermal and thermo-mechanical loading. The results obtained show that the full membrane model can be replaced by the reduced one in determining critical loads and corresponding wavelengths. The gain in computation time obtained is important, due to the coarse mesh required by the reduced model. The reduced model is very sensitive to membrane's boundaries conditions and requires to have a quasi constant wavelength along the membrane width

Ces travaux s'articulent autour du calcul des solutions périodiques dans les systèmes dynamiques non linéaires, au moyen de méthodes numériques de continuation. La recherche de solutions périodiques se traduit par un problème avec conditions aux limites périodiques, pour lequel nous avons implémenté deux méthodes d'approximation : - Une méthode spectrale dans le domaine fréquentiel : l'équilibrage harmonique d'ordre élevé, qui repose sur une formulation quadratique des équations. Nous proposons en outre une formulation originale permettant d'étendre cette méthode aux cas de non-linéarités non rationnelles. - Une méthode pseudo-spectrale par éléments dans le domaine temporel : la collocation à l'aide fonctions polynômiales par morceaux. Ces méthodes transforment le problème continu en un système d'équations algébriques non linéaires, dont les solutions sont calculées par continuation à l'aide de la méthode asymptotique numérique. L'ensemble de ces outils, intégrés au code de calcul MANLAB et complétés d'une analyse linéaire de stabilité, sont alors utilisés pour l'étude des régimes périodiques d'une classe particulière de systèmes dynamiques non linéaires : les instruments de musique auto-oscillants. Un modèle physique non-régulier de clarinette est étudié en détail : à partir de la branche de solutions statiques et ses bifurcations, on calcule les différentes branches de solutions périodiques, ainsi que leur stabilité et leurs bifurcations. Ce modèle est ensuite adapté au cas du saxophone, pour lequel on intègre une caractérisation acoustique expérimentale, afin de mieux tenir compte de la géométrie complexe de l'instrument. Enfin, nous étudions un modèle physique simplifié de violon, avec une non-régularité liée frottement de Coulomb. Cette dernière application illustre ainsi la polyvalence des outils développés face aux différents types de non-régularité.

Le travail de cette thèse a été réalisé dans le cadre d'une collaboration entre EDF R&D et le LMA de Marseille (CNRS). Le but était de développer des outils théoriques et numériques pour le calcul de modes non-linéaires de structures industrielles possédant des non-linéarités localisées à jeu. La méthode de calcul utilisée est une combinaison de la méthode d'équilibrage harmonique (EH) et de la méthode asymptotique numérique (MAN), appelée EHMAN. Elle est réputée pour sa robustesse sur les problèmes réguliers. L'enjeu de ce travail de thèse est de l'appliquer sur des problèmes non-réguliers régularisés de type butée à jeu pour lequel un grand nombre d'harmonique est nécessaire. Des améliorations ont été apportées à la méthode de base pour rendre effectif le traitement de modèles à "grand" nombre de degrés de liberté (DDL). Les développements réalisés pendant la thèse ont été capitalisés par la création de nouveaux opérateurs dans Code_Aster.Une étude approfondie d'un système à 2 degrés de liberté a permis de faire émerger quelques caractéristiques des systèmes non-linéaires à jeu. Celles-ci ont servi entre autre à établir une méthodologie pour l'étude de systèmes à grand nombre de DDL. Pour finir, la potentialité des modes non-linéaires comme outil de diagnostic vibratoire est démontrée avec l'étude d'un tube cintré de générateur de vapeur. Le calcul des modes non-linéaires a monté l'existence d'une interaction entre un mode hors-plan (basse fréquence) et un mode plan (haute fréquence) expliquant des régimes vibratoires non-standards. Ce résultat, impossible à obtenir avec les outils de l'analyse modale linéaire, est confirmé expérimentalement. / This work is a collaboration between EDF R&D and the Laboratory of Mechanics and Acoustics. The objective is to develop theoretical and numerical tools to compute nonlinear normal modes (NNMs) of structures with localized nonlinearities.We use an approach combining the harmonic balance and the asymptotic numerical methods, known for its robustness principally for smooth systems. Regularization techniques are used to apply this approach for the study of nonsmooth problems. Moreover, several aspects of the method are improved to allow the computation of NNMs for systems with a high number of degrees of freedom (DOF). Finally, the method is implemented in Code_Aster, an open-source finite element solver developed by EDF R&D.The nonlinear normal modes of a two degrees-of-freedom system are studied and some original characteristics are observed. These observations are then used to develop a methodology for the study of systems with a high number of DOFs. The developed method is finally used to compute the NNMs for a model U-tube of a nuclear plant steam generator. The analysis of the NNMs reveals the presence of an interaction between an out-of-plane (low frequency) and an in-plane (high frequency) modes, a result also confirmed by the experiment. This modal interaction is not possible using linear modal analysis and confirms the interest of NNMs as a diagnostic tool in structural dynamics.

Dans cette thèse, des techniques d’homogénéisation multi-échelles sont proposées pour l’analyse des vibrations des matériaux composites viscoélastiques. Dans la première partie, la Méthode Asymptotique à Deux Echelles (MADE) est proposée pour la modélisation des vibrations des longues structures sandwichs viscoélastiques répétitives. Pour ce type de structures les pulsations amorties correspondant aux modes propres de vibration sont regroupées en paquets bien distincts. La MADE décompose le problème initial de grande taille en deux problèmes de petites tailles. Le premier est défini sur quelques cellules de base et le second est une équation différentielle d’amplitude à coefficients complexes. La résolution de ces problèmes permet de déterminer les propriétés amortissantes correspondant aux modes de début et de fin de paquet de la structure tout en évitant la discrétisation de toute la structure. Pour les structures dont les coeurs ont un module d’Young dépendant de la fréquence, le problème non linéaire formulé sur les cellules de bases est résolu par l’approche diamant. Les modèles ADF et à dérivées fractionnaires ont été considérés dans les tests numériques. En utilisant la MADE, on évite la discrétisation de toute la structure, ce qui permet donc de réduire considérablement le temps de calcul ainsi que l’espace mémoire CPU nécessaires. L’approche proposée a été validée en comparant les résultats à ceux de la simulation éléments finis basée sur la discrétisation de toute la structure, et utilisant l’approche diamant. Dans la seconde partie de cette thèse, la méthode des éléments finis multi-échelles (EF2) a été développée pour le calcul des propriétés modales des structures à matériaux hétérogènes viscoélastiques en terme de fréquences amorties et amortissements modaux. Dans le principe de l’approche EF2, le problème de vibration est formulé à deux échelles : l’échelle de la structure globale (échelle macroscopique) et l’échelle d’un VER minutieusement choisi (échelle microscopique). Le problème à résoudre à l’échelle microscopique est un problème non linéaire alors que le problème à résoudre à l’échelle macroscopique est un problème linéaire. La non linéarité à l’échelle microscopique est introduite par la dépendance en fréquence du module d’Young des matériaux des phases viscoélastiques. Le problème non linéaire ainsi généré à l’échelle microscopique est résolu grâce à la MAN et ses outils de différentiation automatique réalisés sous Matlab, Fortran et C++. Un outil numérique, générique, robuste, peu coûteux en temps de calcul et espace mémoire CPU, de résolution des problèmes de vibrations non amorties des structures composites viscoélastique est ainsi mis en place. Le modèle viscoélastique à module constant ainsi que des modèles à modules dépendant de la fréquence notamment le modèle ADF et le modèle à dérivées fractionnaires ont été considérés pour les tests numériques de validation. Les comparaisons avec les résultats ABAQUS ont confirmé l’efficacité du code propos é. Le modèle est ensuite utilisé pour le calcul des propriétés amortissantes des structures sandwichs viscoélastiques à coeur composite. Les capacités de la nouvelle approche à concevoir des structures sandwichs viscoélastiques à coeur composite et à haut pouvoir amortissant ont été testées avec succès à travers l’étude de l’influence des différents paramètres des inclusions sur les propriétés amortissantes d’une structure sandwich viscoélastique à coeur composite / In this thesis, multiscale homogenization techniques are proposed for vibration analysis of structures with viscoelastic composite materials. In the first part, the Double Scale Asymptotic Method is proposed for vibration modeling of large repetitive viscoelastic sandwich structures. For this kind of structures, la eigenfrequencies are closely located in well separated packets. The DSAM splits the initial problem of large size into two problems of relatively small sizes. The first problem is posed on few basic cells, and the second one is an amplitude equation with complex coefficients. The resolution of these equations permits to compute the damping properties that correspond to the beginning and the end of every packets of eigenmodes. In case of structure with frequency dependent Young modulus in the core, the diamant approach is used to solve the nonlinear problem posed on basic cells. The ADF and fractional derivative models are considered in numerical tests. By using the DSAM, one avoid the discretization of the whole structure, and the computation time and needed CPU memory are thus reduced. The proposed method is validated by comparing its results with those of the direct finite element method using the diamant approach. In the second part of this thesis, the multiscale finite element method (FE2) is proposed for computation of modal properties (resonant frequency and modal loss factors) of structures with composite materials. In the principle of the (FE2) method, the vibration problem is formulated at two scales: the scale of the whole structure (macroscopic scale) and the scale of a Representative Volume Element (RVE) considered as the microscopic scale. The microscopic problem is a nonlinear one and the macroscopic problem is linear. The nonlinearity at the microscopic scale is introduced by the frequency dependence of the Young modulus of the viscoelastic phases. This nonlinear problem is solved by the Asymptotic Numerical Method and its automatic differentiation tools realizable in Matlab, Fortran or C++. From this approach, numerical tool that is generic, flexible, robust and inexpensive in term of CPU time and memory is proposed for vibration analysis of viscoelastic structures. The constant Young modulus and frequency dependent Young modulus are considered in validation tests. The results of numerical simulation with ABAQUS are used are reference. The model is then used to compute the modal properties of sandwich structure with viscoelastic composite core. To test the capacities of the proposed approach to design sandwich viscoelastic structure with high damping properties, the influence of parameters of the inclusions are studied

Cette thèse s’intéresse au développement d’éléments finis solide–coques dédiés à la modélisation de structures multicouches sollicitées en vibrations. En effet, la plupart des modèles multicouches dans la littérature présentent des limitations dans certaines configurations géométriques et matérielles. Face à ce constat et dans un souci de proposer un outil moins coûteux en temps de calcul, nous avons proposé une approche basée sur le concept solide–coques. Il s’agit d’éléments finis tridimensionnels dont le comportement a été amélioré par l’Assumed Strain Method. Dans un premier temps, nous avons formulé le problème de vibrations de structures sandwichs à cœur viscoélastique. La dépendance en fréquence a ainsi été prise en compte en utilisant une loi constitutive complexe. Pour résoudre le problème discrétisé, la Méthode Asymptotique Numérique, couplée à l’homotopie, et utilisant l’approche DIAMANT, a été adoptée pour les excellents résultats qu’elle offre par rapport aux autres méthodes. Des tests ont permis de valider les modèles proposés et de montrer l’avantage par rapport aux éléments ayant la même cinématique. Poursuivant nos travaux, et dans un souci d’augmenter l’amortissement, nous nous sommes orientés vers un contrôle actif des vibrations. Pour ce faire, deux éléments finis piézoélectriques ont été formulés. Il s’agit des éléments SHB8PSE et SHB20E qui sont des extensions des éléments finis SHB8PS et SHB20, respectivement. Le couplage électromécanique a consisté en l’ajout d’un degré de liberté à chacun des nœuds des dits éléments. Quelques exemples en statique et en vibrations menés sur des structures multicouches allant de simples poutres aux structures présentant des non-linéarités géométriques ont permis de valider les éléments solide–coques proposés. Pour finir, une synthèse des acquis des chapitres 2 et 3 a permis de proposer une modélisation de structures multicouches comprenant des couches élastiques, viscoélastiques et piézoélectriques. À l’amortissement passif provenant du pouvoir amortissant des matériaux viscoélastiques, on ajoute un contrôle actif qui découle du courant électrique généré au cours de la déformation des couches piézoélectriques. Ainsi, un filtre a été installé entre les capteurs et actionneurs. Ce filtre permet d’amplifier ou d’atténuer le potentiel électrique généré dans le but de réduire les amplitudes de vibrations. Pour résoudre le problème résultant nous avons étendu le solveur utilisé au chapitre 2. Pour valider les modèles proposés, des tests de contrôle actif–passif ont été menés sur des structures plaques multicouches. Enfin, quelques lois de contrôle découlant de filtres ont permis de montrer comment cette procédure permet de réduire ou même d’éviter l’amplification des vibrations / This thesis deals with the development of solid–shell finite elements for vibration modeling of multilayer structures. Indeed, most of multilayer models in the literature show some limitations in certain geometric and material configurations. Considering these restrictions and in order to develop a more efficient calculation tool, we proposed an approach based on the solid–shell concept. This consists of three-dimensional finite elements enhanced through the Assumed Strain Method. First of all, we have formulated the problem of vibrations of sandwich structures with viscoelastic core. The frequency dependence has been taken into account by using a complex constitutive law. To solve the discretized problem, the Asymptotic Numerical Method, coupled with the homotopy technique and the DIAMANT toolbox approach, was adopted due to the excellent results it provides compared to other methods. Benchmark tests have validated the models and highlighted their advantages over existing elements having the same kinematics. In order to increase damping properties, we directed our attention towards an active vibration control. For this purpose, two piezoelectric finite elements have been developed. These finite elements SHB8PSE and SHB20E are extensions, of the elements SHB8PS and SHB20, respectively. The electromechanical coupling consisted in adding an electrical degree of freedom to each node of these elements. A variety of test problems both in static and vibration analysis conducted on multilayer structures ranging from simple beams to structures involving geometric nonlinearities allowed validating the proposed solid–shell elements. Finally, combining the achievements made in chapters 2 and 3, we proposed a modeling approach for multilayer structures composed of elastic, viscoelastic and piezoelectric layers. Active control is introduced using the piezoelectric properties in order to improve the reduction in vibration amplitudes. Thus, a filter has been mounted between the sensors and actuators. This filter allows amplifying or attenuating the generated electric potential in order to reduce the vibration amplitudes. To solve the resulting problem, we extended the resolution method used in chapter 2. To validate the proposed models, active–passive control tests have been conducted on multilayer plate structures. Finally, some control laws, associated with filters, have shown how this procedure can allow reducing or even avoiding amplification of vibrations