This dissertation develops and implements an effective numerical scheme, the inclusion-based boundary element method (iBEM), to investigate the mechanical and multi-physical properties of the composites containing arbitrarily shaped particles. Besides the linear elasticity and transient heat conduction problems shown in the dissertation, it can be extended to other problems, such as potential flows and Stokes flows. Through the combination of conventional boundary element method (BEM) and the Eshelby's equivalent inclusion method (EIM), the local field is obtained through superposition of the domain integral of eigen-fields and boundary integral equations. Firstly, the boundary value problems of a composite containing various fully bonding phases of subdomains is introduced. Due to the continuity of displacement (potential) and traction (flux) at the interfaces between different material phases, the interfacial continuity equations are established, which can be solved with the multi-region BEM conventionally. Thanks to Eshelby's celebrated contribution, the material difference in inhomogeneity problems is simulated by an eigenstrain on the inclusion domain but with the same material properties as the matrix. Therefore, the boundary value problems with inhomogeneities can be transformed as domain integral of Green's function with the eigenstrain over the inclusion, where can be determined by the equivalent stress conditions in EIM. Hence, the algorithm of iBEM is formulated and established on the basis of boundary conditions and equivalent stress equations instead of various continuity constraint equations, which saves efforts in computational resources and pre/post-process. The domain integral of Green's function is the key to the algorithm of iBEM, as it bridges the inhomogeneities and the boundary. The closed-form expression of domain integrals for ellipsoidal / elliptical inclusions with polynomial eigenstrain, polygonal and polyhedral inclusions with constant eigenstrain have already existed in the literature. However, it is not applicable to arbitrary particles with varying eigenstrain. This dissertation derives the closed-form domain integrals for polygon and polyhedral inclusions with polynomial eigenstrain source terms, which creates feasibility to solve the local field and effective material properties for composites with arbitrary particles. Although the EIM with polynomial-form eigenstrain has been applied to simulate the material mismatch for ellipsoidal / elliptical inhomogeneities by using the Taylor's of eigenstrain field at the particle center, when it is extended to angular particles, the inaccuracy is significantly reduced due to the rapid and complicated eigenstrain variation in the neighborhood of vertices with the strong singular effects. Therefore, the domain discretization of an angular particle is proposed to tackle the complicated distribution of elastic fields, which keeps the features of exactness (no approximation of interior field) and 𝐂⁰ continuity of eigenstrain. Hereby, the iBEM is proposed to serve as an effective and powerful tool, which takes the advantages of both BEM and EIM. The interaction of inhomogeneities is considered in the process of constructing EIM equations, and boundary effects are taken into account as the contribution to displacement of the eigen-field over inhomogeneities, hence, a complete linear equation system can be established. For the inclusion problems with a prescribed eigenstrain, no domain discretization is required because the exact elastic solution is obtained given the specific dimension of the geometry. Regarding to inhomogeneity problems, 1) the ellipsoidal / elliptical shape is versatile, which could be switched to various of shapes by adjusting the aspect ratio and orientations; 2) though the angular subdomain requires discretization, this method is rapidly convergent and no mesh is needed for the matrix. Therefore, this method enables the simulation of thousands 3𝐷 and 2𝐷 arbitrary shaped particles in a desk-top computer and the effective moduli can be obtained through virtual experiments (i.e, uni-axial loading) or periodic boundary conditions. This method can be easily extended to multi-physical problems, such as transient hear transfer, steady state heat, through changing the fundamental solutions accordingly. Three major packages have been added to the iBEM software, as transient heat transfer, closed-form 2D/3D domain integrals, and domain discretization method. Some case studies demonstrate the capability and applications of this method and software. This main contributions of the PhD studies are as follows: 1) The closed-form domain integrals for polygonal and polyhedral inhomogeneities have been derived based on the gravitational potential theory and transformed coordinates. The solutions are verified with the classic solution of circular and spherical potentials with polynomial source terms (i.e, linear and quadratic) by using many triangular and tetrahedral elements. It enables to solve the inhomogeneity problems with arbitrary particles. 2) Due to the discontinuity on the surfaces and edges of the subdomains and strong singular effects on the vertices, the variation of eigenstrain field is complicated in the neighborhood of edges and vertices. The domain discretization approach is proposed to provide a rapid convergent and effective solution in the infinite space. Different from the Taylor's expansion, the eigenstrain is assigned exactly at the nodes with shape functions instead of at the centroid of the elements, therefore, a 𝐂⁰ continuity is enforced. Here 3-node, 6-node triangular elements and 4-node, 10-node tetrahedral elements are implemented in the code of iBEM, which agree well with FEM but with much fewer of elements. Other types of element are also implementable in the same fashion. 3) The discretization method is applied to investigate the stress singularities of a vertex on an isosceles triangle embedded in an unbounded matrix. Two types of stress singularities are investigated: when the load is applied to the triangular inclusion with the same stiffness as the matrix, the singularity is caused by the irregular load distribution, namely load singularity, and can be exactly evaluated by integral of the potentials on the source with Eshelby's tensor. The second singularity, namely material singularity, is caused by the stiffness mismatch between the triangular inhomogeneity and the matrix under a uniform far field stress, in which the material mismatch is simulated by an eigenstrain. The relationship between the load singularity and material singularity is investigated, and the linkages of these singularities with line distributed force, cracking, and point force are discussed. 4) A parametric study of accuracy on stress field for uniform, linear and quadratic eigenstrain fields was performed and case studies have been presented to demonstrate the capability of iBEM for virtual experiments of ellipsoidal / elliptical inhomogeneities. Subsequently, combining the domain discretization method, iBEM is also applied to study the local elastic fields of the angular inhomogeneities. The effective material behavior is obtained with either large number of particles or periodic boundary condition (PBC) and some interesting discoveries of microstructure-dependent material behavior are reported with the aid of virtual experiments. 5) The iBEM is extended to multiphysical problems. The temperature and hear flux fields of composite materials containing phase change materials (PCM) for energy efficient buildings is demonstrated. Different from the static EIM, the thermal property mismatch between PCM particle and matrix phase is simulated with a uniformly distributed eigen-temperature gradient field and a fictitious heat source on the particle. With the equivalent heat flux conditions and the specific heat-temperature relationship, the eigen-temperature gradient and fictitious heat source can be solved and temperature field of the bounded domain can be calculated. Verified with FEM and laboratory measurements of the transient heat transfer within a building block containing a PCM capsule. Parametric studies have also been conducted to study the influences of the PCM location and volume fraction on the temperature fields of composites with multiple particles. The virtual experiments demonstrate the energy saving and phase delay by using the PCM-concrete wall panel. In summary, the proposed iBEM algorithm bridges the gap between conventional EIM and BEM for virtual experiments of composites samples. The combination of shape functions and domain integrals of polygonal / polyhedral subdomain enables its application to arbitrary shaped particles. It serves as a powerful tool to conduct virtual experiments for composite materials with various geometry and investigate the effective moduli under uni-axial load of samples with large number of particles or under the periodic boundary condition. In the future, the iBEM will be implemented for time independent and dependent nonlinear behavior of composites, such as elastoplastic, viscoelastic, and dynamic elastic problems. In addition to the current parallel computing scheme, GPU can be employed to speed up particle - particle interactions.

<p>The capability of manipulation of the flow of mechanical energy in the form of mechanical waves (including acoustic and elastic waves) has always been a challenge and a critical part in various areas of engineering. The recent advances in topological acoustic/elastic metamaterials certainly open a new pathway to the manipulation of mechanical waves, especially for the novel scattering-immune wave-guiding capability, even in the presence of defects, disorders or sharp bends along the waveguide. In this Dissertation, the theoretical background and experimental evidence of various types of elastic-wave topological metamaterials including analogues to 2D quantum valley Hall effect (QVHE) materials, 2D quantum spin Hall effect (QSHE) topological insulators are presented. First, the formulation the elastic-wave analogue to QVHE materials in a general continuous elastic phononic structure (not limited to local resonant lattices, filling the gap in the literature) is proposed, and a strategy using pressurized cells to actively control the phononic lattice is presented. By finite prestrain and geometric nonlinear effect, the space inversion symmetry of the original hexagonal lattice is broken, resulting in distinct QVHE phases (characterized by valley Chern numbers) in lattice domains with opposite pressurization. With such mechanism, the edge-state path, i.e., the domain wall connecting lattices with distinct QVHE phases, can be real-time configured. Further more, edge states with tunable frequency-wavenumber dispersion can be created at the external boundaries of the lattice by appropriate pressurization of the outermost cells. An aluminum reticular sheet built with water-jet cutting is machined in the pre-deformed pattern with a Z-shape domain wall at the center, which spatially divides the sheet into two domains with opposite QVHE phases. Using piezoelectric transducers and laser Doppler vibrometry, the measured harmonic and transient responses confirm the back-scattering-immunity of the topological edge states, and the frequency-wavenumber dispersion matches the numerical prediction. A strategy is proposed for unidirectionally generating edge states along the domain wall using two off-phase transducers, which is also experimentally demonstrated. For elastic-wave analogue to QSHE topological insulators, we focus on the ``zone-folding'' method and propose a honeycomb 2D elastic beam network with periodically altered thickness with a generalized Kekule distortion pattern. Such framework provides a parametric space with exhaustive control in the topological phase diagram of waves in the lattice compared to earlier works in the literature. The effective Hamiltonian as well as the characterized topological phase are gauge dependent, particularly they change with different reference frames. This lead to ambiguity in the topological phase of such phononic crystal. Based on this argument, it is predicted that edge states could exist at a dislocation interface connecting two piece of phononic structures of the same pattern with relative displacement. Following the same idea, but considering the available fabrication options, a phononic plate with honeycomb groove pattern engraved on both sides is built, which the depth varied according to the Kekule pattern. With proper tuning of the parameters, it realizes an analogue to the QSHE topological insulator. With <em>ab initio</em> calculation of the Berry curvature (without involving any approximations such as the perturbative approach), a new topological invariant <em>local topological charge</em> is defined and evaluated as the counterpart of the Z₂ invariant in the classical-wave-zone-folding analogue. The local topological charge has intrinsic ambiguity and its value depends on the selected reference frame. However, its <em>change </em>according to changes in the parameters, under a consistent reference frame, is well-defined. Given the fact that shifting the reference frame by certain fractions of a lattice constant was equivalent to changing one of the parameters by a certain amount, it also lead to a well-defined change in the local topological charge, which indicates topological phase transition, and one can predict the existence of edge states at the displacement-dislocation interface between two neighboring lattices having the same pattern up to a rigid-body shifting. The phononic plate is machined by a CNC mill, and the experiment is carried out using piezoelectric transducers and laser Doppler vibrometry, which confirms the existence and robustness of the topological edge states at such dislocation interface connecting identical pattern, which was unprecedented in both quantum and classical systems. The final part of this Dissertation focuses on creating classical mechanical analogues to the 1D Kitaev superconducting model and Majorana-like bound states aimed at future acoustic-wave based computation.</p>

L’industrie pétrolière s’intéresse désormais à des régions de la terre de plus en plus difficiles d’accès et il est essentiel de proposer des techniques permettant de garantir l’efficaité d’un forage. Parmi ces techniques, la Reverse Time Migration (RTM) est connue pour sa précision. Elle utilise les ondes réfléchies pour reconstruire une carte du sous-sol représentant les interfaces géophysiques. Elle peut être décrite en trois étapes : (i) propager le champs émis par les sources durant la campagne d’acquisition; (ii) pour chaque source, propager le champ enregistré par les récepteurs; (iii) obtenir une image du sous-sol en appliquant une condition d’imagerie à chaque pas de temps et pour chaque source. Cette technique requiert de très grosses capacités de calcul et il est encore difficile d’imager des milieux réalistes 3D, même avec l’aide du calcul haute performance. Nous avons choisi la méthode de Galerkine discontinue pour modéliser la partie propagation car elle permet d’obtenir des solutions précises et est adaptable au calcul parallèle. La quantité d’information à sauvegarder pour faire une corrélation étant importante, on se doit de trouver un algorithme de calcul d’images du sous-sol réduisant ce coût. Nous avons utilisé l’algorithme de Griewank, appelé “Optimal Checkpointing”. Ce problème de coût étant réglé, on se doit de considérer l’efficacité des ondes élastiques incluant des champs multiples pour améliorer la précision de l’image. La condition traditionnelle de J. Claerbout ne prend pas en compte les conversions d’ondes, et n’est alors surtout utile que dans le cas acoustique. De plus, les ondes P et S interagissant entre elles, il est intéressant de trouver une condition d’imagerie utilisant ce fait. Cela a été abordé dans le cadre de la méthode de l’état adjoint dans les travaux de A. Tarantola et J. Tromp et ce travail en propose utilisation dans le cadre de la RTM. Nous proposons une nouvelle condition d’imagerie prenant en compte les paramètres élastiques du milieu considéré et permettant de supprimer les artefacts numériques. Nous illustrons les images sur des cas industriels / Since a large number of sedimentary basins have been explored, oil exploration is now interested in investigating regions of the Earth which are hostile. Among existing methods for seismic imaging, Reverse Time Migration (RTM) is a technique known by industry to be efficient. The RTM uses reflected waves and is able to construct a map of the subsurface which is depicted by the interfaces limiting the geophysical layers. The algorithm of RTM can be described as a three-step procedure: (i) compute the wavefields emitted by the sources used during the seismic acquisition campaign; (ii) for each source, compute the so-called “backpropagated wavefield”, which is the wavefield obtained by using as sources the signals recorded at the receivers during the acquisition campaign and by reversing the time; (iii) get an image of the subsurface by applying an imaging condition combining the propagated and the backpropagated wavefields at each time step of the numerical scheme and for each source. This technique is computationnaly intensive and it is still difficult to image realistic 3D elastic media, even with the help of HPC. We have thus chosen to consider high-order Discontinuous Galerkin Methods which are known to be well-adapted to provide accurate solutions based upon parallel computing. As we need to correlate a lot of wavefields, we need to find an algorithm reducing the CPU time and the storage : this is the Griewank’s algorithm, so-called “Optimal Checkpointing”. The traditional imaging condition, proposed by J. Claerbout, does not take wave conversions into account and since P-wave and S-wave interact with each other, it might be relevant to use an imaging condition including these interactions. In fact, this has been done successfully by A. Tarantola and J. Tromp for seismology applications based upon the inversion of the global Earth. In this work, we propose a new imaging condition using the elastic parameters which attenuates numerical artifacts. We illustrate the properties of the new imaging condition on industrial benchmarks like the Marmousi model. In particular, we compare the new imaging condition with other imaging conditions by using as criteria the quality of the image.

When the point under consideration is not near to the source of vibration, the strains developed in the rock mass due to the passage of waves are usually of small magnitude, and within the elastic range. However, the rock mass may be subjected to a wide range of strain levels depending on the source of vibration and the wave frequency, even within the elastic limit. The present study is based on the two general conditions existing at field, long wave length propagation of waves and intermediate wavelength propagation of waves. When the wavelength of propagating wave is much longer than the joint spacing, it is referred to as long wavelength condition and is associated with propagation of low frequency waves across closely spaced joints. When wavelength of propagating wave is nearly equal to joint spacing, it is known as intermediate wavelength condition and is associated with propagation of high frequency waves. Long wave length propagation of waves has been studied by conducting laboratory experiments using Resonant Column Apparatus on developed plaster gypsum samples. The influence of joint types, joint spacing and joint orientation on wave propagation has been analyzed at three confining stresses under various strain levels. The wave velocities and damping ratios at various strain levels have been obtained and presented. Shear wave velocities are more dependent on confining stress than compression wave velocities across frictional joints whereas, compression wave velocities are more dependent on confining stress than shear wave velocities across filled joints. Wave velocities are at minimum and wave damping is at maximum across horizontal joints whereas wave velocities are at maximum and wave damping is at minimum across vertical joints. Shear wave velocity and shear wave damping are more dependent on joint orientations than compression wave velocity and compression wave damping. As Resonant Column Apparatus has some limitations in testing stiff samples, a validated numerical model has been developed using Discrete Element Method (DEM) that can provide resonant frequencies under torsional and flexural vibrations. It has been found from numerical simulations, that reduction of normal and shear stiffness of joint with increasing strain levels leads to wave velocity reduction in jointed rock mass. Intermediate wave length propagation of waves has been studied by conducting tests using Bender/ extender elements and the numerical simulations developed using 3DEC (Three Dimensional Distinct Element Code).Parametric study on energy transmission, wave velocities and wave amplitudes of shear and compression waves, has been carried out using the validated numerical model. The propagation of waves across multiple parallel joints was simulated and the phenomenon of multiple reflections of waves between joints could be observed. The transformations of obliquely incident waves on the joint have been successfully modeled by separating the transmitted transformed P and S waves. The frequency dependent behavior of jointed rocks has been studied by developing a numerical model and by applying a wide range of wave frequencies. It has been found that low frequency shear waves may involve slips of rock blocks depending on the strength of rock joint, leading to less transmission of energy; while low frequency compression waves are well transmitted across the joints. High frequency shear and compression waves experience multiple reflections and absorptions at joints.

Submitted by Renata Lopes (renatasil82@gmail.com) on 2017-02-24T17:38:27Z No. of bitstreams: 1 jonathanestebanarroyosilva.pdf: 3851364 bytes, checksum: 7341b01ce42c37de611bb2df24f9012c (MD5) / Approved for entry into archive by Adriana Oliveira (adriana.oliveira@ufjf.edu.br) on 2017-03-06T19:29:14Z (GMT) No. of bitstreams: 1 jonathanestebanarroyosilva.pdf: 3851364 bytes, checksum: 7341b01ce42c37de611bb2df24f9012c (MD5) / Made available in DSpace on 2017-03-06T19:29:14Z (GMT). No. of bitstreams: 1 jonathanestebanarroyosilva.pdf: 3851364 bytes, checksum: 7341b01ce42c37de611bb2df24f9012c (MD5) Previous issue date: 2014-02-24 / FAPEMIG - Fundação de Amparo à Pesquisa do Estado de Minas Gerais / Este trabalho apresenta um novo esquema de marcha no tempo capaz de reduzir oscilações espúrias através de amortecimento numérico para problemas de propagação de ondas elásticas no âmbito da Aproximação Explícita de Green (\Explicit Green's Approach" (ExGA)) [1]. A expressão integral referente ao ExGA é escrita em termos das funções de Green e Degrau. Seus cálculos são realizados de forma independente por meio da formulação semi-discreta do MEF e o método Diferença Central. Devido ao princípio da causalidade, as funções de Green e Degrau possuem um suporte compacto ao redor dos pontos fonte para um intervalo de tempo suficientemente pequeno que é usualmente Empregado nos métodos explícitos clássicos de integração temporal aplicados à modelagem de propagação de ondas. Neste sentido, as funções de Green e Degrau em t = Δt podem ser eficientemente calculadas localmente através de subdomínios pequenos. Cada subdomínio local com sua respectiva submalha cobre somente pontos nodais onde os valores das funções de Green e Degrau são não nulos. A precisão e eficiência da metodologia proposta é demostrada ao analisar três exemplos numéricos. / This work presents a new time-marching scheme able to reduce spurious oscillations by means of numerical damping for elastic wave propagation problems in the framework of the Explicit Green's Approach (ExGA) [1]. The integral expression concerned with the ExGA is written in terms of the Green's and the Step response functions. Their computations are carried out independently by means of the semidiscrete FEM and the Central difference method. Due to the principle of causality, the Green's and Step response functions admit a compact support surround the source points for a small enough time step that is usually employed in common explicit time integration methods applied to wave propagation modeling. In this sense, the Green's and Step response functions at t = Δt can be e ciently computed locally through small subdomains. Each local subdomain with its respective submesh covers only nodes whose Green's and Step response function values do not vanish. The accuracy and e ciency of the proposed methodology are demonstrated by analyzing three numerical examples.

Le contexte scientifique de cette thèse est l'imagerie sismique dont le but est de reconstituer la structure du sous-sol de la Terre. Comme le forage a un coût assez élevé, l'industrie pétrolière s'intéresse à des méthodes capables de reconstituer les images de la structure terrestre interne avant de le faire. La technique d'imagerie sismique la plus utilisée est la technique de sismique-réflexion qui est basée sur le modèle de l'équation d'ondes. L'imagerie sismique est un problème inverse qui requiert de résoudre un grand nombre de problèmes directs. Dans ce contexte, nous nous intéressons dans cette thèse à la résolution du problème direct en régime harmonique, soit à la résolution des équations d'Helmholtz. L'objectif principal est de proposer et de développer un nouveau type de solveur élément fini (EF) caractérisé par un opérateur discret de taille réduite (comparée à la taille des solveurs déjà existants) sans pour autant altérer la précision de la solution numérique. Nous considérons les méthodes de Galerkine discontinues (DG). Comme les méthodes DG classiques sont plus coûteuses que les méthodes EF continues si l'on considère un même problème à cause d'un grand nombre de degrés de liberté couplés, résultat des approximations discontinues, nous développons une nouvelle classe de méthode DG réduisant ce problème : la méthode DG hybride (HDG). Pour valider l'efficacité de la méthode HDG proposée, nous comparons les résultats obtenus avec ceux obtenus avec une méthode DG basée sur des flux décentrés en 2D. Comme l'industrie pétrolière s'intéresse au traitement de données réelles, nous développons ensuite la méthode HDG pour les équations élastiques d'Helmholtz 3D. / The scientific context of this thesis is seismic imaging which aims at recovering the structure of the earth. As the drilling is expensive, the petroleum industry is interested by methods able to reconstruct images of the internal structures of the earth before the drilling. The most used seismic imaging method in petroleum industry is the seismic-reflection technique which uses a wave equation model. Seismic imaging is an inverse problem which requires to solve a large number of forward problems. In this context, we are interested in this thesis in the modeling part, i.e. the resolution of the forward problem, assuming a time-harmonic regime, leading to the so-called Helmholtz equations. The main objective is to propose and develop a new finite element (FE) type solver characterized by a reduced-size discrete operator (as compared to existing such solvers) without hampering the accuracy of the numerical solution. We consider the family of discontinuous Galerkin (DG) methods. However, as classical DG methods are much more expensive than continuous FE methods when considering steady-like problems, because of an increased number of coupled degrees of freedom as a result of the discontinuity of the approximation, we develop a new form of DG method that specifically address this issue: the hybridizable DG (HDG) method. To validate the efficiency of the proposed HDG method, we compare the results that we obtain with those of a classical upwind flux-based DG method in a 2D framework. Then, as petroleum industry is interested in the treatment of real data, we develop the HDG method for the 3D elastic Helmholtz equations.

Elastic waves describe particles vibrating in materials holding the property of elasticity. Particularly, several kinds of resistance in elasticity lead to the models of elastic waves with different damping mechanisms. In the thesis, the influence from friction, structural damping, Kelvin-Voigt damping on the linear and semilinear elastic waves in two or three dimensions are studied. Concerning the Cauchy problem for linear elastic waves, some qualitative properties of solutions including well-posedness, smoothing effect, propagation of singularities, energy estimates and diffusion phenomena, are derived by using WKB analysis associated with diagonalization procedures or the spectral theory. By constructing suitable time-weighted Sobolev spaces and using Banach's fixed point theorem, global (in time) existence of small data solutions to the weakly coupled systems for semilinear elastic waves with different damping terms have been proved. The main tools to treat the nonlinear terms in Sobolev spaces are some fractional tools in Harmonic Analysis. Finally, well-posedness and Lp-Lq estimates for elastic waves without any damping terms in three dimensions are analyzed by employing Riesz transform theory and stationary phase methods.

Un système SHM par ondes guidées a pour but d'évaluer l'intégrité d'une grande variété de structures fines, telles que les fuselages d'avions, les tuyaux, les réservoirs, etc. Un tel système est basé sur un réseau de capteurs piézoélectriques pour l'excitation et la mesure des ondes guidées. Cette méthode de SHM par ondes guidées est prometteuse pour l'inspection de structures de grande taille, ces ondes se propageant sur de grandes distances avec peu d'atténuation, tout en étant sensibles aux défauts surfaciques et subsurfaciques.Cette thèse présente les travaux menés dans le but de développer un système de SHM par ondes guidées capable de détecter, localiser et dimensionner efficacement les défauts dans des structures aéronautiques assimilables à des plaques, en matériaux composites ou en aluminium. Des simulations et des méthodes d'apprentissage sont utilisées pour déterminer les caractéristiques principales des ondes guidées propagées, notamment les vitesses de phase et de groupe ainsi que la fonction de Green 3D. Celles-ci sont ensuite utilisées pour traiter les signaux des ondes guidées afin de produire des images représentant l'intégrité des structures étudiées. Ce travail comprend également une étude approfondie des algorithmes d'imagerie DAS, MV et Excitelet, les plus prometteurs parmi ceux de la littérature, une évaluation de leurs performances par analyse statistique sur une grande base de données de résultats de simulations d'imagerie par ondes guidées et propose une méthode d'imagerie parcimonieuse.Alors que la détection et la localisation des défauts à partir de l'analyse des images est aisée, le dimensionnement du défaut est un problème plus complexe en raison de sa forte dimensionnalité et de sa non-linéarité. Il est démontré que ce problème peut être résolu par des méthodes d'apprentissage automatique sur une grande base de données de résultats de simulations d'imagerie par ondes guidées. Ces méthodes d'imagerie nécessitent cependant une référence, mesurée sur la structure dans un état supposé sain. Elles sont efficaces dans des conditions opérationnelles stationnaires mais sont sensibles aux variations de l'environnement et notamment aux fluctuations de température.Ce travail présente donc l'étude de la robustesse face aux effets thermiques des méthodes d'imagerie par ondes guidées et propose un modèle de détection de défauts capable d'analyser des résultats d'imagerie détériorés. Plusieurs techniques de compensation des effets thermiques sont étudiées et des améliorations sont proposées. Leur efficacité est validée pour les plaques d'aluminium mais des améliorations supplémentaires sont nécessaires pour les étendre aux plaques de composites. / A guided wave-based structural health monitoring (SHM) system aims at determining the integrity of a wide variety of plate-like structures, including aircraft fuselages, pipes, tanks etc. It relies on a sparse array of piezoelectric transducers for guided waves (GWs) excitation and sensing. With a number of benefits, these waves are standing out among other methods as a promising method for the inspection of large structures. They can propagate on significant distances with small attenuation while being sensitive to surface and subsurface defects.This thesis presents studies conducted with the purpose of developing such a GWs-based SHM system that is capable of efficient defect detection, localization and sizing aeronautical plate-like structures made of aluminum and composite materials. Simulation and data-driven approaches are presented for determining principal characteristics of propagating GWs, namely modal group and phase velocities, 3D Green's functions etc. in structures of interest. They are then used for GWs signals processing in order to compute images representing the integrity of studied structures. This work also provides a comprehensive overview of DAS, MV and Excitelet defect imaging algorithms, determines their performance using statistical analysis of an extensive dataset of simulated guided waves imaging (GWI) results and proposes a method for sparse defect imaging.While defect detection and localization are straightforward from the image analysis, the defect sizing is a more complex problem due to its high dimensionality and non-linearity. It is demonstrated that this problem can be solved by means of machine learning methods, relying on an extensive database of simulated GWI results. Aforementioned defect imaging methods are baseline demanding. They are efficient under stationary operational conditions but vulnerable to environmental variations, especially to the temperature fluctuation.Finally, this work presents studies on the robustness of GWI methods against thermal effects, and a defect detection model capable of analyzing deteriorated GWI results is proposed. Different techniques for thermal effects compensation are reviewed, and improvements are proposed. Their effectiveness is validated for aluminum plates but further improvements are required to translate these techniques to composite plates.

Le thème central de cette thèse est l'étude des résonances pour le problème intérieur en élastodynamique (géométries elliptique et sphéroïdale), et pour le problème de diffusion en acoustique (géométrie elliptique). On s'intéresse en particulier à la levée de dégénérescence des résonances liée à la brisure de symétrie de l'objet lors de la transition du disque circulaire vers le disque elliptique (2D), et de la sphère vers le sphéroïde (3D). Ce phénomène est étudié et interprété d'un point de vue théorique en prenant en compte les symétries de l'objet à l'aide de la théorie des groupes. Cette approche est complétée par une modélisation numérique et une partie expérimentale. En 2D, nous étudions le problème intérieur pour un disque elliptique élastique (étude des modes résonants) et le problème de la diffusion acoustique par des cylindres elliptiques élastiques. Ils sont traités à partir du formalisme modal combiné à la théorie des groupes dans le contexte vectoriel de l'élastodynamique. La levée de dégénérescence est observée théoriquement mais aussi expérimentalement en diffusion. La méthode simplifie considérablement le traitement numérique des problèmes étudiés, fournit une classification des résonances selon les 4 représentations irréductibles du groupe de symétrie C2v (associé à la géométrie elliptique) et donne une interprétation physique de la levée de dégénérescence en termes de brisure de symétrie. Une partie expérimentale en spectroscopie ultrasonore vient compléter l'étude théorique du problème de diffusion. Une série d'expériences en cuve est menée dans le cas de cylindres elliptiques de différentes excentricités en aluminium immergés dans l'eau, dans la bande de fréquence 0 ≤ kr ≤ 50, où kr est le nombre d'onde réduit dans le fluide. Les résultats expérimentaux présentent un très bon accord avec les résultats théoriques, la levée de dégénérescence est observée expérimentalement sur des fonctions de forme et mise en évidence sur des diagrammes angulaires. Le problème intérieur en 3D est traité expérimentalement à partir de la génération et la détection optiques d'ondes élastiques. Une série d'expérimentations sur des objets tridimensionnels (sphère, sphéroïdes oblates et prolates de différentes excentricités) en aluminium est réalisée. Ils sont mis en vibration par impacts laser et les mesures de vitesse et de fréquence s'effectuent par vibrométrie laser. On réalise ainsi une comparaison qualitative entre la théorie 2D et l'expérience 3D. Les mesures sont menées à la fois dans les domaines temporel et fréquentiel pour mettre en évidence la levée de dégénérescence d'une part, et l'onde de Rayleigh qui se propage sur la surface des objets d'autre part. Nous identifions deux trajets pour cette onde en géométrie sphéroïdale, l'un circulaire et l'autre elliptique.Enfin, dans le cadre des problèmes intérieurs 2D et 3D, on donne une interprétation en termes de rayons à travers la dualité entre le spectre des résonances et le spectre des longueurs des orbites périodiques (OPs), avec la mise en évidence du phénomène de conversion de mode et l'identification de l'onde de Rayleigh. Un phénomène, nouveau à notre connaissance, vient s'ajouter au phénomène de bifurcation de certaines orbites. Au cours de la déformation vers le disque elliptique, les orbites avec conversion de mode du disque circulaire se séparent en deux orbites dont les longueurs sont associées aux trajets minimal et maximal qu'elles parcourent. Cette observation s'interprète comme une conséquence du théorème de Fermat. Dans le cas du sphéroïde, on retrouve les orbites du disque circulaire dans le plan équatorial et celles du disque elliptique dans le plan méridien. Nous mettons également en évidence les pics associés aux deux trajets parcourus par l'onde de Rayleigh sur le spectre des OPs. / Resonances for the interior problem in elastodynamics and the acoustic scattering problem are considered in elliptical and spheroidal geometries. Interest is focused on the splitting up of resonances which occurs when the symmetry is broken in the transition from the circular disc to the elliptical one (2D), and from the sphere to the spheroid (3D). From the theoretical point of view, this physical phenomenon is studied and interpreted taking into account the symmetries of the object with the help of group theory. This approach is completed by a numerical modeling and an experimental part. As far as the two dimensional problems are concerned, the interior problem for an elliptical elastic disc (study of resonant modes) and the acoustic scattering problem for infinite elliptical elastic cylinders are studied combining modal formalism and group theory in the vectorial context of elastodynamics. The splitting up of resonances due to the circular symmetry breaking is not only theoretically observed but also experimentally for the scattering problem. The method significantly simplifies the numerical treatment of the problems studied, provides a full classification of resonances over the 4 irreducible representations of the symmetry group C2v (associated with the elliptical geometry) and gives a physical interpretation of the splitting up in terms of symmetry breaking of the symmetry group O(2) (invariance under rotation). An experimental part based on ultrasonic spectroscopy complements the theoretical study of the scattering problem. A series of tank experiments is carried out in the case of aluminum elliptical cylinders immersed in water, in the frequency range 0 ≤ kr ≤ 50, where kr is the reduced wave number in the fluid. The experimental results provide a very good agreement with the theoretical ones, the splitting up is observed on experimental form functions and the split resonant modes are identified on angular diagrams. The interior problem in 3D is studied by means of an experimental approach based on the optical generation and detection of elastic waves. A series of experiments is performed on three-dimensional objects in aluminium. These objects (sphere, prolate and oblate spheroids of various eccentricity) are excited by laser impacts, and the velocity and frequency measurements are carried out by laser vibrometry. Theory and experiments are qualitatively compared. The measurements are performed in both the frequency and time domains to highlight the splitting up of resonances on one hand, and the Rayleigh wave propagating on the surface of the objects on the other hand. We emphasize two paths for this surface wave in spheroidal geometry: a circular one in the equatorial plane and an elliptical one in the meridian plane. Finally, in the context of the interior problems in 2D and 3D, a physical interpretation of resonances in terms of geometrical paths is provided. Mode conversion is highlighted and the Rayleigh wave is identified on the periodic orbits lengths spectrum.In addition to the bifurcations of some periodic orbits, a phenomenon, new to our knowledge, appears. The orbits with mode conversion of the circular disc split in two orbits when the transition to the elliptic disc occurs. The lengths of these orbits are associated with the minimal and maximal travel paths. This observation is interpreted from Fermat's theorem.For the spheroid, orbits of the circular disc and those of the elliptical disc are recovered in the equatorial and meridian planes respectively. We also emphasize the peaks associated with the travel paths of Rayleigh wave in spheroidal geometry appearing on the periodic orbits spectrum.

Les instabilités paramétriques non linéaires (NL) ont été observées sur les ondes magnéto-élastiques dans le cas d’un couplage de trois quasi-phonons sous pompage électromagnétique. La théorie en prédit une dynamique supercritique explosive, mais limitée expérimentalement par le décalage de fréquence dû aux fortes nonlinéarités. La dynamique supercritique des instabilités paramétriques NL est étudiée dans deux matériaux antiferromagnétiques "plan facile" (AFEP): l’hématite α-Fe2O3 et le borate de fer FeBO3. Ces matériaux possèdent une très grande NL acoustique effective en raison du couplage magnéto-élastique élevé. Les mécanismes de limitation de la dynamique explosive ont été analysés à l'aide de l'approximation anharmonique. La compensation du décalage fréquentiel NL par une modulation de phase singulière du pompage a été proposée et théoriquement vérifiée, puis utilisée pour l’observation expérimentale de la dynamique supercritique explosive des excitations de trois quasi-phonons dans les résonateurs magnéto-élastiques. Les études sur FeBO3 ont été réalisées dans la gamme de température 77 K - 293 K où les paramètres magnéto-élastiques du cristal varient de façon significative. Un modèle fortement non linéaire des excitations de trois quasi-phonons dans les AFEPs a été développé. Les simulations numériques sont en accord avec les résultats expérimentaux. Les études théoriques de couplage de trois ondes magnéto-élastiques progressives ont été effectuées sur la base de modèles théoriques prenant en compte la non-linéarité cubique des cristaux AFEP réels. Les simulations numériques prévoient un comportement explosif et une localisation spatiale des triades générées / Recently discovered nonlinear parametric instabilities occur when nonlinear parameter of a system is modulated. These instabilities were reported on magnetoelastic waves as three quasi-phonon coupling under electromagnetic pumping. Theoretical studies predicted supercritical explosive dynamics of these instabilities. Experimentally such singular behavior is limited by strong nonlinear frequency shift.Presented work studies supercritical dynamics of nonlinear parametric instabilities in two easy plane antiferromagnets (AFEP): hematite α-Fe2O3 and iron borate FeBO3. These materials possess unprecedented effective acoustic nonlinearity due to high magneto-elastic coupling. Limiting mechanisms of explosive dynamics were analyzed with the help of anharmonic approximation. Nonlinear frequency shift compensation via singular pumping field phase modulation was suggested and theoretically approbated. This technique was successfully used for experimental observation and investigation of supercritical explosive dynamics of three quasi-phonon excitations in magnetoelastic resonators. Iron borate studies were performed in the temperature range 77 K – 293 K where magnetoelastic parameters of the crystal vary essentially. Strongly nonlinear model of three quasi-phonon excitations in AFEPs was developed. Numerical simulations of the model showed good agreement with experimental results.Theoretical studies of three travelling magnetoelastic waves coupling are performed on the basis of suggested theoretical models that take into account cubic nonlinearity of real AFEP crystals. Numerical simulations of the models suggest explosive behavior and spatial localization of generated triads