La Tomographie Optique Diffuse Résolue en Temps (TOD-RT) est une technique d'imagerie clinique et préclinique en pleine croissance. Elle fournit les cartes d'absorption et de diffusion optique des organes explorés, et les paramètres physiologiques associés. La Tomographie Optique Diffuse de Fluorescence Résolue en Temps (TODF-RT) est basée sur la détection de photons de fluorescence. Elle permet de déterminer les cartes de la concentration et du temps de vie de sondes fluorescentes et d'accéder à une imagerie métabolique et moléculaire très importante pour le diagnostic et le suivi thérapeutique, en particulier en cancérologie. L'objectif de cette thèse était de réaliser des images 3D de TOD/TODF-RT sur des rongeurs en utilisant une technologie optique résolue en temps, les acquisitions étant réalisées à l'aide de fibres optiques disposées autour de l'animal et sans contact avec sa surface. Le travail a été mené en quatre étapes : 1- mise en place d'un dispositif d'imagerie de la surface de l'animal et reconstruction de son contour 3D, 2- modélisation de l'approche sans contact pour la résolution du problème direct, 3- traitement des mesures prenant en compte la réponse impulsionnelle de l'appareil, 4- établissement d'une méthode de reconstruction des images basée sur une sélection de points judicieusement choisis sur les profils temporels. Ces travaux ont permis d'obtenir des images optiques 3D de bonne qualité en réduisant la diaphonie entre l'absorption et la diffusion. Ces améliorations ont été obtenues tout en diminuant le temps de calcul, par comparaison avec les méthodes utilisant la totalité des profils temporels.

Ce travail de thèse est consacré à l'étude des renforts textiles techniques 2D à l'échelle mésoscopique. La méthode des éléments finis est utilisée pour résoudre un problème aux limites, fortement non linéaire, dans le domaine du renfort fibreux. Deux nouveaux modes de déformations pour le comportement des mèches de renforts ont été développés et caractérisés. Le premier mode est la compression longitudinale, qui permet de traduire la faible rigidité des mèches lorsqu'elles subissent une dilatation négative dans le sens des fibres. La relation conflictuelle sur le plan numérique entre la rigidité en tension, très forte, et la rigidité en compression longitudinale, très faible, peut être résolue via trois méthodes : réduction du pas de temps critique, addition de la contribution en tension ou avec une nouvelle stratégie pour l'actualisation du champ de contrainte. Le second mode de déformation est la dilatation transversale des mèches considérée comme conséquence directe de la compression longitudinale. Ce phénomène d'expansion de matière dans les directions transverses peut être observé avec un essai de compression longitudinale in-situ sous tomographie X et est communément ignorer. Un protocole expérimental a été mise en place pour mesurer cette dilatation transversale des mèches et un coefficient de Poisson a pu être caractérisé par identification inverse. Une campagne expérimentale a permis d'identifier les paramètres matériau du modèle mésoscopique et les résultats de simulations sont comparés aux images issues d'essai mécanique in situ sous tomographe. / This thesis is devoted to the mesoscopic study on the performance of textile reinforcements. F.E. simulation is carried out on a mesoscale model for the fibrous material, based on which two kinds of new deformation modes are developed. The first one is a longitudinal compression mode, which is used to reflect the small stiffness when the yarn is compressed longitudinally. The incompatibility problem between the small longitudinal compression stiffness and the large tension stiffness are solved by three different strategies: constraining the critical step time, adding the nonlinear tension part, or using a new strategy to update the stress. The second one is transverse expansion mode that could reflect the influence from longitudinal deformation to transverse deformation. This deformation could be found in tomography view but was ignored by the former researches. An experiment is designed to measure the expansion magnitude, and the geometrical inverse fitting process is applied to measure the value of the longitudinal-transverse Poisson ratio. The parameters of the mesoscale model are measured by a series of mechanical experiments and the simulation results are verified by the tomography methodology.

Ce travail, à l'interface entre les mathématiques appliquées et la physique, s'articule autour de la modélisation numérique des vésicules biologiques, un modéle pour les globules rouges du sang. Pour cela, le modéle de Canham et Helfrich est adopté pour décrire le comportement des vésicules. La modélisation numérique utilise la méthode des fonctions de niveau dans un cadre éléments finis. Un nouvel algorithme de résolution numérique combinant une technique de multiplicateurs de Lagrange avec une adaptation automatique de maillages garantit la conservation exacte des volumes et des surfaces. Cet algorithme permet donc de dépasser une limitation cruciale actuelle de la méthode des fonctions de niveau, à savoir les pertes de masse couramment observées dans ce type de problémes. De plus, les propriétés de convergence de la méthode des fonctions de niveau se trouvent ainsi grandement améliorées, comme l'indiquent de nombreux tests numériques. Ces tests comprennent notamment des problémes d'advection élémentaires, des mouvements par courbure moyenne ainsi que des mouvements par diffusion de surface. Concernant l'équilibre statique des vésicules, une condition générale d'équilibre d'Euler-Lagrange est obtenue à l'aide d'outils de dérivation de forme. En dynamique, le mouvement d'une vésicule sous l'action d'un écoulement de cisaillement est étudié dans le cadre des nombres de Reynolds élevés. L'effet du confinement est considéré, et les régimes classiques de chenille de char et de basculement sont retrouvés. Finalement, pour la premiére fois, l'effet des termes inertiels est étudié et on montre qu'au delà d'une valeur critique du nombre de Reynolds, la vésicule passe d'un mouvement de basculement à un mouvement de chenille de char. / This work, at the interface between the Applied Mathematics and Physics is connected about the numerical modelisation of biological vesicles, a pattern for the red blood cells. For this reason, the pattern of Canham and Helfrich is adopted to describe the behaviour of the vesicles. The numerical modelisation uses the Level Set method in finite element framework. A new algorithm of numerical resolution combining one technique of Lagrange multipliers with an automatic mesh adaptation ensures the accurate conservation of volumes and surfaces. Thus this algorithm enables to exceed an existing crucial restriction of the Level Set method, that's to say, the wastes of mass usually noticed in this kind of problems. Moreover, the proprieties of convergence of the Level Set method are thus much more improved, as shown in many numerical tests. Those tests chiefly include elementary problems of advection, motions by mean curvature just as motions by spread of surface. Concerning the static equilibrum of the vesicles, a mechanical equilibrum equation (Euler-Lagrange equation) of a vesicle membrane under a generalized elastic bending energy is obtained and the approach is based on shape optimization tools. In dynamics, the motion of a vesicle under the effect of a shear flow is elaborated in the frames of reference of high Reynolds numbers. The effect of confinement is respected, and the standard regimes of tank-treading and of tumbling motion are found again. Finally, for the first time, the effect of the inertia terms is elaborated and we show that beyond a critical value of the number of Reynolds the vesicle passes from a tumbling motion to a tank-treading motion.

Le but de cette thèse est de développer les outils numériques et semi-analytiques qui nous permettent d’étudier des guides non-linéaires complexes et réalistes qui pourraient être fabriqués et caractérisés expérimentalement.Nous présentons une étude complète d’une version améliorée du guide d’onde plasmonique non-linéaire à une dimension, en ajoutant deux couches tampons de diélectrique linéaire entre le cœur non-linéaire isotrope et les deux gaines métalliques. Ses couches réduisent les pertes et permettent leur diminution avec la puissance contrairement aux guides simples. De plus, les principaux modes plasmoniques non-linéaires peuvent présenter une transition spatiale vers des modes de types différents, qui peut être contrôlée par la puissance.Par la suite, nous étudions un nouveau guide plasmonique non-linéaire à fente utilisant un métamatériau, soit dans le cœur non-linéaire soit dans les gaines linéaires. Nous avons mis au point une méthode semi-analytique et une méthode numérique afin d’étudier les solutions stationnaires non-linéaires dans ce nouveau guide non-linéaire anisotrope. Nous avons montré analytiquement et numériquement que le cœur non-linéaire anisotrope peut être conçu afin d’atteindre de forts effets non-linéaires à faible puissance. Pour certains métamatériaux dans les gaines linéaires, la figure de mérite de ce guide d’onde augmente de plus de 50 fois par rapport aux guides isotropes.Pour conclure, nous présentons une nouvelle méthode basée sur la méthode des éléments finis de type vectoriel couplée à l’algorithme à puissance fixée pouvant calculer rigoureusement les effets non-linéaires dans des guides d’onde plasmoniques 2D. / The main goal of this PhD is to develop the semi-analytical and the numerical tools that allow us to study complicated and realistic nonlinear plasmonic waveguides which can be fabricated and characterized experimentally.First, we present a full study of an improved version of the one-dimensional nonlinear plasmonic slot waveguide, by adding two linear dielectric buffer layers between the isotropic nonlinear core and the two metal claddings. These additional layers reduce the overall losses and allow the losses to decrease with the power for some configurations unlike the usual slot. Furthermore, the main plasmonic modes can exhibit nonlinear spatial modal transitions towards new families of modes that can be controlled with the power.Second, we propose and study new one-dimensional nonlinear plasmonic slot waveguides with metamaterial regions either in the nonlinear core or in the linear claddings. For the metamaterial nonlinear core, we developed semi-analytical and fully numerical methods in order to study the nonlinear stationary solutions propagating in this anisotropic nonlinear waveguide. We have demonstrated both analytically and numerically that the anisotropic nonlinear core can be designed in order to achieve strong nonlinear effects at low power.For the structures with metamaterial linear claddings, the figure of merit can be extremely enhanced by more than 50 times compared with the simple one. Finally, we present the full derivation of a new nonlinear full vectorial finite element method based on the fixed power algorithm in order to quantify rigorously the nonlinear characteristics of realistic two-dimensional nonlinear plasmonic structures.

Les structures maçonnées sont largement connues en génie civil comme constituant d’une partie des bâtiments, mais également en tant que garnissages réfractaires dans des structures utilisées à hautes températures, par exemple en sidérurgie. Malheureusement, les outils actuels ne sont pas suffisamment puissants pour prédire le comportement de ces structures avec l’apparition de fissures et pour tenir compte du comportement non linéaire d’un des deux constituants (le mortier par exemple). Ce travail de thèse contribue à la modélisation multi-échelles des maçonneries classiques et des garnissages réfractaires avec un faible coût numérique grâce à la technique d’homogénéisation périodique. Les techniques de modélisation et de simulation du comportement des maçonneries sont présentées et développées. L’influence des lois d’interface entre briques et mortier, des paramètres géométriques et matériels, ainsi que de la densité des fissures sur le comportement effectif des maçonneries est étudié. Trois approches (une extension analytique de Cecchi et Tralli, une approche numérique et un modèle micromécanique) sont proposées pour la détermination du comportement effectif d’une cellule périodique dans le cas de maçonneries avec mortier viscoélastique microfissuré et briques saines élastiques ou rigides. Les résultats des calculs sur deux exemples de maçonneries (1D et 2D) ont confirmé que l’approche multi-échelle est une solution appropriée avec une grande capacité à exprimer le comportement des maçonneries viscoélastiques microfissurées. Ce travail, limité au cas sans propagation de fissures, peut être étendu aux mortiers à comportement viscoplastique. / Masonry structures are widely used in civil engineering as part of buildings or in refractory linings of structures working at high temperatures, like in steel industry. Unfortunately, the present tools are not powerful enough to predict the behavior of these structures at their micro-cracked state and/or if one of their constituents behaves nonlinearly (e.g. the mortar). This research contributes to the multi-level modeling of classical masonries and refractory linings with a low numerical cost using basically the periodic homogenization technique. Modeling and simulation techniques of masonry behavior are presented and developped. The influence of interface law between bricks and mortar, of geometrical and material parameters, and of crack density on the effective masonry behavior is studied. Three approaches (analytical extension of Cecchi and Tralli, numerical approach and micromechanical modeling) were proposed to determine the effective behavior of a periodic masonry cell with micro-cracked viscoelastic mortar and safe elastic or rigid bricks. The results obtained on two examples of masonry (1D and 2D) confirmed that the multi-scale approach is a suitable solution with a great ability to model the effective behavior of microcracked viscoelastic masonry. This work, actually limited to the case without crack propagation, could be extended to mortars with viscoplastic behavior.

Dans le cadre de la lutte contre les nuisances sonores et vibratoires, cette thèse porte sur la modélisation numérique des structures amorties par dispositifs piézoélectriques shuntés. La première partie du travail concerne la modélisation par éléments finis de structures en vibrations avec des pastilles piézoélectriques shuntées. Dans un premier temps, une formulation éléments finis originale, qui utilise des variables électriques globales (différence de potentiel et charge dans chaque pastille piézoélectrique), est analysée et validée. Dans un second temps, différentes stratégies de réduction de modèle basées sur la méthode de projection modale sont proposées pour résoudre le problème électromécanique discrétisé par éléments finis à moindre coût. La convergence de ces modèles d’ordre réduits est ensuite analysée pour les cas de shunts résistif et résonant. La deuxième partie du travail est consacrée à l’optimisation du système électromécanique, dans le but de maximiser l’amortissement apporté par les dispositifs piézoélectriques shuntés. Pour cela, une procédure d’optimisation topologique, basée sur la méthode SIMP (Solid Isotropic Material with Penalization method), est développée pour déterminer les géométries et les emplacements optimaux des pastilles piézoélectriques. Cette procédure permet de maximiser le coefficient de couplage électromécanique modal entre les éléments piézoélectriques et la structure hôte, ceci de façon indépendante du choix des composants du circuit électrique. Les avantages de l’approche proposée sont mis en avant à travers un exemple de validation et un cas d'application industrielle. Enfin, la dernière partie du travail propose une approche numérique pour modéliser et optimiser la réduction du rayonnement acoustique de plaques minces dans le domaine des basses fréquences avec des éléments piézoélectriques shuntés. Cette approche est valable pour n’importe quelle plaque mince bafflée et non trouée, indépendamment des conditions aux limites. Un exemple d’application concernant l’atténuation du rayonnent acoustique d’une plaque avec renforts est présenté et analysé. / Passive structural vibration and noise reduction by means of shunted piezoelectric patches is addressed in this thesis. The first part of the work concerns the finite element modeling of shunted piezoelectric systems. Firstly, an original finite element formulation, with only a couple of electric variables per piezoelectric patch (the global charge/ voltage), is analyzed and validated. Secondly, several reduced order models based on a normal mode expansion are proposed to solve the electromechanical problem. The convergence of these reduced order models is then analyzed for a resistive and a resonant shunt circuits. In the second part of the work, the concept of topology optimization, based on the Solid Isotropic Material with Penalization method (SIMP), is employed to optimize, in terms of damping efficiency, the geometry of piezoelectric patches as well as their placement on the host elastic structure. The proposed optimization procedure consists of distributing the piezoelectric material in such a way as to maximize the modal electromechanical coupling factor of the mechanical vibration mode to which the shunt is tuned, independently of the choice of electric circuit components. Numerical examples validate and demonstrate the potential of the proposed approach for the design of piezoelectric shunt devices. Finally, the last part of the work concerns the numerical modeling of noise and vibration reduction of thin structures in the low frequency range by using shunted piezoelectric elements. An efficient approach that can be applied to any thin continuous plates in an infinite baffle, independently of the boundary conditions, is proposed. An application example of a thin plate with reinforcements is presented and analyzed.

Une méthode d’analyse des topographies résiduelles, obtenues par essais de nanoindentation Berkovich sur trois échantillons de nickel polycristallins cubiques à faces entrées (CFC), a été développée dans cette thèse. L’objectif de la méthode est d’évaluer la richesse de l’information contenue dans les empreintes pour l’identification de tout ou partie des paramètres d’une loi de plasticité monocristalline. La création d’une base de données de topographies résiduelles, mesurées par microscopie à force atomique (AFM), constitue la partie expérimentale du travail. Les distributions et dimensions de bourrelets montrent une grande sensibilité à l’orientation relative indenteur/grain et au taux d’écrouissage. Les topographies obtenues s’avèrent être de véritables « empreintes digitales » du mécanisme de plasticité à l’échelle du grain. L’élaboration sous le code ZéBuLon d’un modèle éléments finis (EF) 3D de l’essai de nanoindentation Berkovich, intégrant la loi de Méric-Cailletaud, permet de retrouver les observations expérimentales. Une étude numérique confirme la sensibilité de la topographie à l’orientation relative indenteur/grain et aux paramètres plastiques, notamment aux coefficients de la matrice d’interaction des dislocations présentes sur les systèmes de glissement. Afin d’évaluer la richesse du contenu informatif des empreintes, un indice d’identifiabilité est proposé. Son calcul est basé sur la multi-colinéarité des vecteurs de sensibilité des topographies résiduelles aux paramètres de la loi. Il permet de quantifier, en fonction des données topographiques prises en compte, le caractère mal posé du problème d’identification paramétrique. Les résultats obtenus montrent que l’identification de quatre à cinq paramètres de la loi de Méric-Cailletaud est envisageable en exploitant seulement deux empreintes. Ces travaux ouvrent la voie à l’identification du comportement à l’échelle du cristal, guidée par l’identifiabilité paramétrique. / In this thesis, a method is developed to analyse the residual topographies obtained by Berkovich nanoindentation tests on three face-centered cubic (FCC) polycrystalline nickel samples. The purpose is to measure the information richness of imprints for identifying all or part of parameters of a single crystal plasticity law. The experimental part consists in creating a residual topography database by atomic force microscopy (AFM) measurements on the samples. Pile-up distributions and sizes are largely sensitive to the indenter/grain relative orientation and the hardening rate. The topographies are true “fingerprints” of the plasticity mechanism at the grain scale. A 3D finite element (FE) modelling of nanoindentation test is developed using the code ZeBuLon and making use of the Méric-Cailletaud law. Numerical results show a good agreement with experimental observations and are largely sensitive to the indenter/grain relative orientation and the plastic parameters, including the interaction matrix coefficient specifying the interaction between dislocations on different slip systems. To measure the imprint information content, an identifiability index is proposed. Its calculation is based on the multicollinearity among the sensitivity vectors of topographies to the law parameters. According to the considered topographies, it measures if the numerical model updating problem is ill-posed. The results show that four to five parameters of the Méric-Cailletaud law can be identified by considering two topographies. This work paves the way for identifying the material behaviour at the grain scale using parametric identifiability

Le béton présente une rupture quasi-fragile du fait de la présence d’une zone d’élaboration de la rupture principalement microfissurée et de taille conséquence, se développant en avant du front de fissure. L’objectif de ces travaux consiste à mettre en évidence le processus d’évolution de la zone d’élaboration et l’étendue de cette dernière, et ce, sur la base d’une campagne d’essai Wedge Splitting suivie de simulations numériques aux éléments-finis utilisant le modèle d’endommagement isotrope de Fichant – La Borderie à l’échelle mésoscopique. Expérimentalement, les courbes de résistance obtenues dans le cadre de la Mécanique Linéaire Élastique de la Rupture équivalente attestent d’un phénomène de confinement précoce de la zone d’élaboration manifestement lié à la géométrie testée et aux propriétés du béton étudié.Les simulations numériques obtenues sur la base du modèle d’endommagement sont en accord avec les résultats expérimentaux et conduisent également à soupçonner ce phénomène de confinement. Par la suite, la simulation numérique est à nouveau utilisée sur cette même géométrie de spécimen mais en réduisant les propriétés de rupture de la matrice cimentaire afin de diminuer la taille de la zone d’élaboration. Les résultats de cette nouvelle simulation montrent un développement libre de la zone d’élaboration suivie d’une propagation auto-similaire de la fissure principale attestée par la présence d’un régime plateau de la courbe de résistance correspondante.Une analyse détaillée du champ d’endommagement puis du champ d’énergie restituée est réalisée et permet de mettre en évidence un critère de développement de la zone d’élaboration fondée sur une valeur maximale du taux local de restitution d’énergie d’endommagement. / Concrete exhibits a quasibrittle fracture due to the existence of a large fracture process zone (FPZ), mainly microcracked, which develops ahead of the crack front. The aim of the current work consists in highlighting the FPZ development and its extent. For that purpose,an experimental campaign using Wedge Splitting Test was carried out and followed by finite element simulation using Fichant – La Borderie isotropic damage model acting at the mesoscale of concrete. Experimental analysis exhibits that by the use from Resistance curves estimated in the framework of equivalent Linear Elastic Fracture Mechanics, the used geometry combined to the studied concrete properties are subjected to a confinement of FPZ. Numerical simulations achieved by the damage model are in agreement with experimental results and also seem to show FPZ confinement. There after, numerical simulations are used again on the same specimen geometry but by decreasing fracture properties of cementitious matrix, in order to minimize the FPZ size. Numerical results exhibit that the FPZ was now freely developed and followed by a self-similar propagation of the main crack as shown by the existence of a “plateau” value on the corresponding Resistance curve. A numerical analysis is performed and leads to the propositionof a FPZ development criteria based on a maximal value of the local damage energy release rate.

Il est connu depuis des décennies que la propagation des fissures de fatigue dans les matériaux élastoplastiques est très sensible à l’histoire du chargement car le comportement non-linéaire du matériau peut avoir une grande influence sur les vitesses de propagation. Cependant, le calcul brut de millions de cycles de fatigue avec des comportements matériaux non-linéaires sur des structures tridimensionnelles réalistes conduirait à des temps de calcul prohibitifs. Ainsi, nous proposons de coupler deux approches de réduction de modèle a priori et a posteriori, afin de diminuer considérablement le coût de calcul de ce type de problèmes. Tout d’abord, considérant l’hypothèse de plasticité confinée, une stratégie de réduction de modèle a posteriori du comportement plastique de la structure fissurée est proposée. Le modèle réduit ainsi obtenu fournit incrémentalement l’état plastique autour du front de fissure, duquel est déduite la vitesse instantanée de la fissure. De plus, une seconde approche de réduction de modèle, a priori cette fois, est aussi mise en place afin d’accélérer encore plus les temps de résolution du problème global. Cette approche a priori consiste à construire incrémentalement —et sans calculs préalables— une base réduite spécifique à chaque cas-test, en extrayant de l’information des champs de déplacement de la structure au cours du temps et pendant la propagation éventuelle de la fissure. Ainsi, les champs de déplacement solutions de la géométrie fissurée réactualisée sont vus comme une combinaison linéaire de cette base réduite de vecteurs. La méthode numérique considérée ici est la méthode des éléments finis. De fait, pendant la propagation de la fissure, la discrétisation spatiale du modèle doit être réactualisée afin d’être conforme avec le front de la fissure. Dans ce but, une technique spécifique de déformation de maillage est utilisée, et permet de discrétiser la géométrie variable du modèle avec des maillages de même topologie. Cette technique de déformation de maillage apparaît comme une étape clé de la stratégie de réduction de modèle. Finalement, une approche adaptative est construite autour de cette stratégie. Elle permet de garantir la qualité des résultats obtenus par rapport à un critère de précision donné. La précision et l’efficacité de cette stratégie globale sont démontrées à travers de nombreux exemples bidimensionnels et tridimensionnels dans le cadre de propagation de fissure en model, de même que pour un exemple industriel d’une pièce fissurée d’hélicoptère. / It has been known for decades that fatigue crack propagation in elastic-plastic media is very sensitive to load history since the nonlinear behavior of the material can have a great influence on propagation rates. However, the raw computation of millions of fatigue cycles with nonlinear material behavior on tridimensional structures would lead to prohibitive calculation times. In this respect, we propose a global model reduction strategy, mixing both the a posteriori and a priori approaches in order to drastically decrease the computational cost of these types of problems. First, the small scale yielding hypothesis is assumed, and an a posteriori model reduction of the plastic behavior of the cracked structure is performed. This reduced model provides incrementally the plastic state in the vicinity of the crack front, from which the instantaneous crack growth rate is inferred. Then an additional a priori model reduction technique is used to accelerate even more the time to solution of the whole problem. This a priori approach consists in building incrementally and without any previous calculations a reduced basis specific to the considered test-case, by extracting information from the evolving displacement field of the structure. Then the displacement solutions of the updated crack geometries are sought as linear combinations of those few basis vectors. The numerical method chosen for this work is the finite element method. Hence, during the propagation the spatial discretization of the model has to be updated to be consistent with the evolving crack front. For this purpose, a specific mesh morphing technique is used, that enables to discretize the evolving model geometry with meshes of the same topology. This morphing method appears to be a key component of the model reduction strategy. Finally, the whole strategy introduced above is embedded inside an adaptive approach, in order to ensure the quality of the results with respect to a given accuracy. The accuracy and the efficiency of this global strategy have been shown through several examples; either in bidimensional and tridimensional cases for model crack propagation, including the industrial example of a helicopter structure.

Pendant sa mission, un lanceur est soumis à des excitations large bande, sévères, aérodynamiques, de provenances diverses, qui peuvent mettre en danger la survivabilité de la charge utile et de l’équipement électronique du véhicule, et par conséquent le succès de la mission. Les structures aérospatiales sont généralement caractérisées par l’utilisation de matériaux composites exotiques des configurations et des épaisseurs variantes, ainsi que par leurs géométries largement complexes. Il est donc d’une importance cruciale pour l’industrie aérospatiale moderne, le développement d’outils analytiques et numériques qui peuvent prédire avec précision la réponse vibroacoustique des structures larges, composites de différentes géométries et soumis à une combinaison des excitations aéroacoustiques. Récemment, un grand nombre de recherches ont été menées sur la modélisation des caractéristiques de propagation des ondes au sein des structures composites. Dans cette étude, la méthode des éléments finis ondulatoires (WFEM) est utilisée afin de prédire les caractéristiques de dispersion des ondes dans des structures composites orthotropes de géométries variables, nommément des plaques plates, des panneaux simplement courbés, des panneaux doublement courbés et des coques cylindriques. Ces caractéristiques sont initialement utilisées pour prédire la densité modale et le facteur de perte par couplage des structures connectées au milieu acoustique. Par la suite, la perte de transmission (TL) à large bande des structures modélisées dans le cadre d’une analyse statistique énergétique (SEA) dans un contexte ondulatoire est calculée. Principalement en raison de la complexité géométrique importante de structures, l’utilisation des éléments finis (FE) au sein de l’industrie aérospatiale est souvent inévitable. L’utilisation de ces modèles est limitée principalement à cause du temps de calcul exigé, même pour les calculs dans la bande basses fréquences. Au cours des dernières années, beaucoup de chercheurs travaillent sur la réduction de modèles FE, afin de rendre leur application possible pour des systèmes larges. Dans cette étude, l’approche de SOAR est adoptée, afin de minimiser le temps de calcul pour un système couplé de type structurel-acoustique, tout en conservant une précision satisfaisante de la prédiction dans un sens large bande. Le système est modélisé sous diverses excitations aéroacoustiques, nommément un champ acoustique diffus et une couche limite turbulente (TBL).La validation expérimentale des outils développés est réalisée sur un ensemble de structures sandwich composites orthotropes. Ces derniers sont utilisés afin de formuler une approche couche équivalente unique (ESL) pour la modélisation de la réponse spatiale du panneau dans le contexte d’une approche de matrice de raideur dynamique. L’effet de la température de la structure ainsi que du milieu acoustique sur la réponse du système vibroacoustique est examiné et analysé. Par la suite, un modèle de la structure SYLDA, également fait d’un matériau sandwich orthotrope, est testé principalement dans le but d’enquêter sur la nature de couplage entre ses divers sous-systèmes. La modélisation ESL précédemment développée est utilisé pour un calcul efficace de la réponse de la structure dans la gamme des basses et moyennes fréquences, tandis que pour des fréquences plus élevées, une hybridisation WFEM / FEM pour la modélisation des structures discontinues est utilisé. / During its mission, a launch vehicle is subject to broadband, severe, aeroacoustic and structure-borne excitations of various provenances, which can endanger the survivability of the payload and the vehicles electronic equipment, and consequently the success of the mission. Aerospace structures are generally characterized by the use of exotic composite materials of various configurations and thicknesses, as well as by their extensively complex geometries and connections between different subsystems. It is therefore of crucial importance for the modern aerospace industry, the development of analytical and numerical tools that can accurately predict the vibroacoustic response of large, composite structures of various geometries and subject to a combination of aeroacoustic excitations. Recently, a lot of research has been conducted on the modelling of wave propagation characteristics within composite structures. In this study, the Wave Finite Element Method (WFEM) is used in order to predict the wave dispersion characteristics within orthotropic composite structures of various geometries, namely flat panels, singly curved panels, doubly curved panels and cylindrical shells. These characteristics are initially used for predicting the modal density and the coupling loss factor of the structures connected to the acoustic medium. Subsequently the broad-band Transmission Loss (TL) of the modelled structures within a Statistical Energy Analysis (SEA) wave-context approach is calculated. Mainly due to the extensive geometric complexity of structures, the use of Finite Element(FE) modelling within the aerospace industry is frequently inevitable. The use of such models is limited mainly because of the large computation time demanded even for calculations in the low frequency range. During the last years, a lot of researchers focus on the model reduction of large FE models, in order to make their application feasible. In this study, the Second Order ARnoldi (SOAR) reduction approach is adopted, in order to minimize the computation time for a fully coupled composite structural-acoustic system, while at the same time retaining a satisfactory accuracy of the prediction in a broadband sense. The system is modelled under various aeroacoustic excitations, namely a diffused acoustic field and a Turbulent Boundary Layer (TBL) excitation. Experimental validation of the developed tools is conducted on a set of orthotropic sandwich composite structures. Initially, the wave propagation characteristics of a flat panel are measured and the experimental results are compared to the WFEM predictions. The later are used in order to formulate an Equivalent Single Layer (ESL) approach for the modelling of the spatial response of the panel within a dynamic stiffness matrix approach. The effect of the temperature of the structure as well as of the acoustic medium on the vibroacoustic response of the system is examined and analyzed. Subsequently, a model of the SYLDA structure, also made of an orthotropic sandwich material, is tested mainly in order to investigate the coupling nature between its various subsystems. The developed ESL modelling is used for an efficient calculation of the response of the structure in the lower frequency range, while for higher frequencies a hybrid WFEM/FEM formulation for modelling discontinuous structures is used.