L’objectif de ce travail est de modéliser la propagation des ondes cohérentes antiplanes dans une distribution aléatoire et uniforme de fissures fermées ou ouvertes, parallèles ou aléatoirement orientées. Elles contiennent un fluide visqueux où leurs lèvres sont libres de contraintes. Les mécanismes de diffusion multiple entre fissures sont pris en compte. Une première étude porte sur la réponse d’une couche endommagée à faces parallèles, sollicitée à incidence normale. Les champs de déplacement cohérent siégeant dans les trois régions de l’espace sont déterminés. Par suite, la masse volumique et la rigidité effectives de la couche, vue homogène par l’onde cohérente, sont définies. Les cas d’une distribution de fissures à concentration variable et d’une onde de surface (de type Love) en présence d’une couche endommagée sur un substrat sain sont alors traités. Une dernière application concerne les populations de fissures visqueuses à tailles variables présentes dans la croûte terrestre. / Abstract

Des panneaux sandwichs en nid d'abeilles sont largement utilisés, notamment dans l’industrie aérospatiale et aéronautique, à cause du très bon rapport entre rigidité en flexion et poids. Concernant leur modélisation, ils sont considérés classiquement comme de milieux homogénéisés équivalents afin d'éviter des modèles numériques prohibitifs en coûts de calculs. Cependant, des travaux précédents ont montré que, si le comportement dynamique en membrane des sandwichs peut être correctement représenté par des modèles homogénéisés classiques dans une large gamme de fréquences, ces mêmes modèles ne permettent malheureusement pas de bien décrire le comportement en flexion dans le domaine de hautes fréquences (HF). En effet, la couche centrale en nid d'abeilles joue un rôle important dans le comportement en flexion du sandwich, il est donc indispensable de la modéliser de manière appropriée. Or, lorsque les longueurs d’onde impliquées deviennent aussi petites que les longueurs caractéristiques des cellules du nid d’abeilles, cette microstructure cellulaire interagit fortement avec les ondes et génère des effets d’interaction non négligeables, qui ne sont malheureusement pas pris en compte par des modèles homogénéisés classiques. Dans le cadre de cette thèse, on s’intéresse donc à l'amélioration de l’analyse théorique et numérique de la propagation d’ondes élastiques HF dans ces panneaux composites. On exploite les caractéristiques périodiques du nid d'abeilles en utilisant sur une approche numérique basée sur la théorie des ondes de Bloch. En effet, en décomposant des solutions non périodiques sur une base composée de modes périodiques de Bloch, il est possible de développer des modèles numériques, qui considèrent des phénomènes de propagation des ondes à l’intérieur d’une seule cellule de base et captent toutes les interactions. Ces modèles numériques sont donc de taille raisonnable, par rapport aux dimensions souvent très importantes des structures industrielles. Des analyses théoriques et des outils de modélisation ont été développés pour des milieux périodiques composés de structures minces : poutres ou plaques. Notre approche a été développée et validée pour des structures périodiques uni- puis bi-dimensionnelles composées de poutres. Pour les cas 2D, la forme de la cellule est hexagonale ou rectangulaire. Nous avons aussi considéré des plaques sandwichs en nid d’abeilles. Pour toutes ces structures, en identifiant les valeurs propres et les modes propres de Bloch sur une cellule primitive pour tous les vecteurs d’onde de Bloch situés dans la première zone de Brillouin dans l’espace de phase, la relation de dispersion entre le vecteur d'onde de Bloch et la valeur propre est calculée. En analysant cette relation de dispersion, les résultats importants sont obtenus, tels que les bandes de fréquences passantes et bloquantes et les caractéristiques d'anisotropie et dispersives des structures périodiques, la comparaison quantitative entre les premiers modes de Bloch et ceux des modèles homogénéisés classiques en vue d’une définition précise du domaine validation en fréquence de ceux-derniers et la mise en évidence des modes de Bloch « rétro-propagatifs » munis d’une vitesse de groupe négative. / Honeycomb core sandwich panels are widely used in the aeronautic industry due to their excellent flexural stiffness to weight ratio. Generally, classical homogenized model is used to model honeycomb core sandwiches in order to have an efficient but not expensive numerical modeling. However, previous works have shown that, while the homogenized models could correctly represent the membrane waves’ behavior of sandwiches in a large frequency range, they could not give satisfying simulation results for the flexural waves’ behavior in the high frequency range (HF). In fact, the honeycomb core layer plays an important role in the propagation of the flexural waves, so that when the involved wavelengths become close to the characteristic lengths of honeycomb cells, the cellular microstructure starts interacting strongly with the waves and its effect should no longer be neglected, which is unfortunately not the case of the homogenized models. In the present work, we are interested in improving the theoretical and numerical analysis of HF elastic waves’ propagation in honeycomb core sandwich panels by a numerical approach based on the Bloch wave theorem, which allows taking into account the periodic characteristics of the honeycomb core. In fact, by decomposing non-periodic wave solutions into their periodic Bloch wave basis modes, numerical models are defined on a basic cell and solved in a efficient way, and provide a better description and so a better understanding of the interaction between HF wave propagation phenomena and the periodic structures. Our numerical approach is developed and validated in the cases of one-dimensional periodic beam structures, of two-dimensional periodic hexagonal and rectangular beam structures and of honeycomb core sandwich plates. By solving the eigenvalue problem of the Bloch wave modes in one primitive cell of the periodic structure for all the wave vectors located in the corresponding first Brillouin zone in the phase space, the dispersion relation between the wave vector and the eigenvalue is calculated. The analysis of the dispersion relation provides important results such as: the frequency bandgaps and the anisotropic and dispersive characteristics of periodic structures, the comparison between the first Bloch wave modes to those of the classical equivalent homogenized models and the existence of the retro-propagating Bloch wave modes with a negative group velocity.

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.