Ce travail de thèse a pour objet d’analyser le comportement en statique, en fatigue et en vibration linéaire et non linéaire des matériaux sandwichs en présence d’une décohésion de longueur variable. Une étude détaillée est d’abord menée pour caractériser le comportement mécanique en statique et en fatigue de ces matériaux. Les essais ont été conduits en flexion 3-points sur des poutres de ces matériaux pour plusieurs distances entre appuis et pour plusieurs longueurs de fissure. En vibration, une étude expérimentale de la réponse en fréquence à une impulsion, menée à l’aide d’un vibromètre laser a permis de mesurer les fréquences propres et les amortissements de ces matériaux autour de chaque pic de résonance en fonction de la longueur de fissure. Les résultats déduits de l’analyse expérimentale sont comparés à ceux obtenus à partir d’une analyse par éléments finis. Enfin, une méthode de vibration non linéaire a été appliquée pour caractériser le comportement des matériaux sandwichs endommagés par fissuration. Les paramètres non linéaires relatifs au décalage fréquentiel et à l’amortissement sont mesurés en faisant varier l’amplitude d’excitation. Cette étude a permis de montrer que les paramètres non linéaires sont plus sensibles à l’endommagement que les paramètres linéaires. / The aim of this work is to investigate the effects of debonding lengths on the static, fatigue, linear and nonlinear vibration behaviour of sandwich materials. First, a study was conducted in static and cyclic fatigue loading with various debonding lengths. Shear and flexural modulus in static tests were determined using the sandwich plate theory. The effects of debonding lengths on the stiffness, hysteresis loops and damping were studied for various numbers of cycles during fatigue tests. Then, modelling of the damping of a composite sandwich with debonding was established considering finite element analysis which evaluated the different energies dissipated in the material directions of the core and the skins. The effects of debonding variable lengths on natural frequencies and damping were studied numerically and compared with experimental results. Finally, the nonlinear vibration method was used to characterize the behaviour of sandwich beams with debonding. The nonlinear parameters corresponding to the elastic modulus and damping were determined for each frequency mode and each debonding length. The results showed that nonlinear parameters were much more sensitive to damage than linear parameters.

La forte sensibilité des ouvrages bois aux transmissions latérales, vibrations et autres sons basses fréquences représente un obstacle à leur développement. En réponse à cette problématique, le projet Vibracoubois, dans lequel s'inscrivent ces travaux, vise à développer la caractérisation ainsi que la modélisation des transferts vibro-acoustiques au sein de ces structures. L'objectif principal de ces travaux est donc d'élargir le champ de l'étude expérimentale et numérique à une structure complète constituée de différentes parois et de cavités couplées. Les jonctions, véritables nœuds des systèmes constructifs, sont placées au centre de l'étude. Pour ce faire, un protocole expérimental est développé, permettant d'effectuer des mesures en laboratoire mais également in-situ sur des structures réelles. Les résultats de campagnes d'essais effectuées au sein de deux maquettes d'ouvrages bois (ossature et bois massif) grandeur nature sont présentés. Par ailleurs, en vue de se rapprocher des jonctions réelles du bâtiment, des méthodes numériques de couplage incluant des effets locaux de dissipation sont développées. Une méthode hybride est présentée. Un modèle numérique de la maquette en bois massif, basé sur la méthode des éléments-finis, est créé. L'utilisation de cette méthode est justifiée de par la prépondérance des comportements modaux aux basses fréquences. Ce modèle intègre la méthode hybride en vue d'obtenir un couplage souple entre les parois, par opposition au couplage classique rigide des éléments-finis. Les résultats d'un processus de recalage portant sur les propriétés des parois et des jonctions ainsi que sur le couplage vibro-acoustique sont présentés. / Wooden buildings are very sensitive to flanking transmissions, vibrations and low frequency sounds, which interferes with their development. The Vibracoubois project, in which this research work is part of, aims to respond to this problem by improving the vibro-acoustic transfers characterization and modelling of these structures. This research work mainly focuses on expanding the experimental and the numerical scope of the study to an entire structure made of coupled walls and air cavities. The junctions, considered critical for the vibration propagation in these buildings, are placed at the center of this work. First, an experimental protocol is presented. It allows to measure both the vibro-acoustic behaviour of laboratory and real structures. The results of test campaigns run on two full-scale experimental mock-ups (timber frame and CLT) are then presented. Secondly, numerical coupling methods that include local dissipation effects are developed, in order to get closer to real junctions behaviour. A hybrid coupling method is presented. A numerical model of the CLT mock-up, based on the finite-elements method, is then created. The use of this method is justified because of the strong modal behaviour of these structures at low frequencies. This model implements the hybrid coupling method in order to create flexible junctions, as opposed to the classic rigid coupling of the finite-elements. Based on the CLT mock-up experimental measurements, the results of an adjustment process of the wall and the junction properties used in the numerical model is finally presented.

Le soudage par électrofusion est la technique majoritairement utilisée pour assembler les tubes et les accessoires en polyéthylène utilisés sur le réseau gazier. Ses très bons résultats initiaux ont été ternis ces dernières années par un certain nombre de dysfonctionnements relevés tant sur le terrain qu'en laboratoire. Ils se traduisent par des soudures de qualité médiocre du fait de la présence de zones de très faible cohésion. L'objectif de ce travail est de développer des outils experts tant expérimentaux que numériques permettant d'optimiser les conditions de soudage par électrofusion. L'originalité de notre étude a été de développer un modèle éléments finis tridimensionnel prenant en compte les différents mécanismes et phénomènes physiques sous-jacents, responsables de la formation de la soudure. Le modèle permet de calculer le taux de transformation de la matière, de prendre en compte les enthalpies de fusion et de cristallisation, de calculer le taux d'interdiffusion des macromolécules à l'interface entre les différents corps à souder pour finalement prédire la qualité de la soudure en fonction des conditions de chauffage imposées. Pour valider le modèle numérique, tout en facilitant l'accès aux grandeurs expérimentales, nous avons conçu et réalisé un dispositif expérimental se présentant sous la forme d'une géométrie relativement simple et plane mais respectant les caractéristiques d'un accessoire réel. La confrontation entre résultats numériques et expérimentaux a permis de démontrer les capacités du modèle numérique à reproduire fidèlement la réalité. Les différents outils ont alors été utilisés pour étudier l'influence des conditions de soudage sur la soudure et pour étudier le soudage de géométries plus complexes telles que celles rencontrées dans les pièces industrielles. / The Electrofusion welding process is widely used to join polyethylene components in gas distribution networks. Even trusty as a technique, the field feedback points out some divergences whose influence on the long term performance of the weld. One of the well-known consequences of these divergences is the “sticking” (aka “cold weld”) that is the result of an uncompleted or even inexistent interdiffusion of the macromolecules of the materials to join. Most numerical simulations are two-dimensional whereas the process is usually three-dimensional both in terms of heat transfer and mechanical aspects. The main objective of the work was to develop a 3D finite element model and to validate it by comparing to a real situation the temperature evolution and the thermally affected areas in a simple planar welding geometry with the same dimensional characteristics as a real fitting chosen to make easier the instrumentation. The numerical model takes into account the fitting parameters such as polyethylene thermal and mechanical properties (i.e. melting and crystallization kinetics, phase transition, thermal expansion) and the electrical and geometrical settings. It computes a criterion based on the macromolecular interdiffusion theory able to determine whether a good welding occurred or not at the end of the welding cycle. The computed results (temperature, melted and cold areas and fracture surfaces) were compared with experimental data and gave very good agreement in terms of temperature, liquid phase fraction distribution and fracture surfaces. Finally the numerical model and the experimental process were used to study the influence of welding conditions on the weld itself and to study the welding of complex geometries such as the industrial fittings.

Cette thèse traite de la modélisation du frittage à l’échelle du Volume Élémentaire Représentatif de la pastille de matériau. L’objectif est de développer des outils numériquesde compréhension des phénomènes physiques mis en jeu lors du frittage. Le domaine d’application ciblé est la fabrication du combustible nucléaire. Une approche multi-Échelle a été mise en oeuvre. Dans un premier temps une modélisation à l’échelle d’un empilement, basée sur la méthode des Éléments Discrets, a été adoptée. Différentes études utilisant cette approche ont été proposées dans la littérature ces dernières années. Tous ces travaux utilisent une méthode discrète explicite. Si certains résultats ont pu être validés expérimentalement,une des limites vient de l’utilisation des méthodes explicites dontle pas de temps critique est très petit. Afin d’augmenter le pas de temps, la masse des particules y est augmentée artificiellement de plusieurs ordres de grandeur. Or,il a été démontré que cette pratique conduit, dans certains cas, à une diminution du réarrangement des particules au sein de l’empilement. Dans cette thèse, une méthode Éléments Discrets implicite appelée Dynamique des Contacts a été adaptée au frittage. Elle permet l’utilisation d’un pas de temps très supérieur à celui des méthodes discrètes explicites et ne nécessite pas d’augmenter artificiellement la massedes particules. La comparaison entre la Dynamique des Contacts et la Méthode des Éléments Discrets explicite montre que notre approche conduit à une représentation plus fidèle du réarrangement. Une validation expérimentale par Microtomographie X ainsi qu’une étude paramétrique sur le frittage des poudres bidispersés sont également présentées pour montrer les possibilités de l’approche discrète appliquée au frittage.La seconde partie est consacrée à une modélisation à l’échelle de deux particules parla méthode des Éléments Finis. Ce modèle repose sur une approche mécanique et vise à représenter de façon plus précise le comportement de deux particules en contact. Les diffusions au joint de grains, en surface et en volume peuvent être représentées. Pour le moment, seules les diffusions en surface et au joint de grains ont été étudiées. Si certaines optimisations restent nécessaires pour que le code soit fonctionnel, plusieurs aspects apparaissent déjà déterminants, comme la courbure de la surface à proximité du joint de grains. A l’avenir, le modèle Dynamique des Contacts du frittage pourra être complété etamélioré grâce aux éléments apportés par le modèle mécanique à l’échelle du grain. / This thesis deals with the simulation of the sintering of nuclear fuel on a pellet scale. The goal is to develop numerical tools which can contribute to a better understandingof the physical phenomena involved in the sintering process. Hence, a multi scale approach is proposed. First of all, a Discrete Element model is introduced. It aims at modeling the motion of particles on a Representative Elementary Volume scale using an original Discrete Element Method. The latter is a Non Smooth Method called Contact Dynamics. Recently, there have been numerous papers about the simulation of sintering using Discrete Element Method. As far as we know, all these papers use smooth methods. Different studies show that the results match well experimental data. However, some limits come from the fact that smooth methods use an explicit scheme which needsvery small time steps. In order to obtain an acceptable time step, the mass of particles have to be dramatically increased. The Non Smooth Contact Dynamics uses an implicit scheme, thus time steps can be much larger without scaling up the mass of particles. The comparison between smooth and non smooth approaches shows thatour method leads to a more realistic representation of rearrangement. An experimental validation using synchrotron X-Ray microtomography is then presented, followedby a parametric study on the sintering of bimodal powders that aims at showing the capacity of this model.The second part presents a mechanical model on the sub-Granular scale, using a Finite Element method. This targets a better understanding of the behavior of twograins in contact. The model is currently being developped but the first results already show that some parameters like the shape of the surface of the neck are very sensitive.In the future, the Non smooth Contact Dynamics model of sintering may be improvedusing the results obtained by the sub-Granular scale mechanical model.

L'objectif de cette thèse est d'évaluer les modifications du métabolisme secondaire des plantes contaminées aux éléments trace métalliques (ETM) et leurs conséquences sur les communautés bactériennes rhizosphériques associées incluant des bactéries présentant des phénotypes de MultiDrug Résistance (MDR). Nous nous sommes focalisés sur deux contextes de sols exposés aux métaux : la phytoremédiation de sites miniers au Vietnam et la reconversion de sols agricoles contaminés par la re-déposition atmosphérique d'activités métallurgiques en France. Nos résultats ont mis en évidence que la contamination par différents types de métaux (dont Cu et Pb principalement) a conduit à une altération de la production des métabolites secondaires des racines, tiges et feuilles de la plante hyperaccumulatrice Pteris vittata et que concernant les racines des tendances similaires dans les changements métaboliques ont pu être observés dans un autre type de contexte de pollution (Zn et Pb plus particulièrement). De même, les profils métaboliques des parties souterraines (racines et rhizomes) de Miscanthus x giganteus ont été modifiés par les concentrations en Pb, Cd et Zn des sols agricoles. Pour les deux plantes examinées des dérivés de l'acide chlorogénique ont été retrouvés en proportions augmentées dans les racines malgré des contextes de nature des sols et de pollutions métalliques très contrastés. Cependant, les dérivés de tanin catéchiques sont spécifiquement trouvés en proportions plus élevées dans les racines de P. vittata sous pression métallique. Ces polyphénols sont connus pour leur capacité à piéger les radicaux libres et leur pouvoir antioxydant et pourraient donc être impliqués dans l'adaptation de ces plantes au stress métallique en contribuant à limiter le stress oxydatif généré par les ETM. Au niveau des parties aériennes, nous n'avons étudié que le changement pour P. vittata et avons mis en évidence une proportion plus élevée de dérivés flavonoïdiques pour les plantes contaminées. Nos résultats de métagénomique nous permettent de conclure également sur un effet des ETM sur la diversité et la richesse spécifique des communautés bactériennes des sols étudiés : une forte contamination en Cu (10 fois la limite autorisée) a diminué la diversité et la richesse bactérienne, alors que pour des niveaux en ETM plus modérés incluant Cu, Pb et Zn, la diversité des communautés bactériennes rhizosphériques semble plus influencée par la plante ou la saison plutôt que par l'effet des ETM. Cet effet sur la composition bactérienne de la rhizosphère de P. vittata se traduit par un enrichissement de certains genres connus comme pathogènes opportunistes de l'homme, notamment Ralstonia, Acinetobacter, Burkholderia et Mycobacterium. En outre, le genre Cupriavidus, connu comme très résistant aux ETM est le seul genre spécifiquement associé à P. vittata qui ait été augmenté au sein de la communauté rhizosphérique pour les deux sites miniers étudiés par rapport aux sols rhizosphériques non pollués. Ce genre pourrait donc être impliqué dans le processus d'adaptation de cette plante au stress métallique. Quant aux communautés rhizosphériques de Miscanthus x giganteus, la sélection de Stenotrophomonas et Pseudomonas dans les sols agricoles contaminés a été observée. Dans le cadre de cette thèse nous avons également mis au point une méthode rapide pour tester l'impact de métabolites végétaux sur des souches pathogènes d'origine clinique et environnementale et également évaluer leur activité inhibitrice de pompes à efflux (IPE) de la famille des RND. Nos données ont ainsi permis de mettre en évidence des activités intéressantes et comparables à celle de l'inhibiteur de pompe à efflux PAßN pour des composés testés qui étaient extraits des racines de Fallopia x bohemica ou des dérivés de ces derniers. / The objective of this thesis is to evaluate the modification of plant secondary metabolism production contaminated with metallic trace elements (MTE) and its consequences on the associated rhizospheric bacterial communities including bacteria presenting MultiDrug Resistant (MDR) phenotypes. We have focused on two contexts of metals exposure: the phytoremediation of mining sites in Vietnam and the reconversion of agricultural soils contaminated by the atmospheric re-deposition of metallurgical activities in France. Our results highlighted that contamination by different types of metals (mainly Cu and Pb) has led to an alteration in the production of secondary metabolites in the roots, stems and leaves of the hyper-accumulating Pteris vittata and for roots, a similar trend in the metabolic changes could be observed in another type of pollution context (Zn and Pb more particularly). Similarly, the metabolic profiles of the underground parts (roots and rhizomes) of Miscanthus x giganteus were modified by the concentrations of Pb, Cd and Zn in agricultural soils. For the two plants examined chlorogenic acid derivatives have been found in increased proportions in the roots despite soil type and pollution context were highly contrasted. However, catechic tannin derivatives are specifically found in higher proportions in the roots of P. vittata under metal pressure. These polyphenols are known for their ability to scavenge free radicals and their antioxidant properties and thus could be involved in the adaptation of these plants to metallic stress by helping to limit the oxidative stress generated by MTE. At the level of the aerial parts, we studied only the change for P. vittata and evidenced higher proportions of flavonoid derivatives for contaminated plants. Our metagenomic results allow us to conclude also on the effect of MTE on the diversity and the specific richness of the bacterial communities of the studied soils: a high contamination of Cu (10 times the allowed limit) decreased dramatically bacterial richness and diversity, while for more moderate MTE levels including Cu Pb and Zn, the diversity of rhizosphere bacterial communities was more explained by plant or season effect rather than an effect of MTE. This effect on P.vittata rhizosphere bacterial composition is reflected by an enrichment in genera known as opportunistic human pathogens, including Ralstonia, Acinetobacter, Burkholderia and Mycobacterium. In addition, Cupriavidus, known as a highly resistant genus, is the only P. vittata specifically associated genus found in increased proportions at both mining sites compared to non-contaminated rhizosphere soils. This genus could then be involved in the adaptation process of this plant with metal stress. As for the rhizospheric communities of Miscanthus x giganteus, the selection of Stenotrophomonas and Pseudomonas in agricultural soils contaminated with MTE was observed. As a part of this thesis, we have also developed a rapid method for testing the impact of plant metabolites on pathogenic strains of clinical and environmental origin and their efflux pump inhibition (EPI) activity of RND family. Our data thus showed interesting and notable EPI activities comparable to that of the efflux pump inhibitor PAßN for tested compounds issued from Fallopia x bohemica roots or for their derivatives.

Les procédés de trempe sont largement répandus dans l'industrie en particulier dans le domaine de l'automobile, du nucléaire et de l'aérospatiale car ils ont un impact direct sur la microstructure, les propriétés mécaniques et les contraintes résiduelles de pièces critiques. La trempe est un processus fortement non-linéaire à cause des couplages forts entre la mécanique des fluides, les transferts thermiques aux différentes interfaces, les transformations de phase du solide et l'ébullition du milieu de trempe. Malgré les progrès effectués par la simulation numérique, ce procédé reste extrêmement difficile à modéliser.Dans ce travail, nous proposons le développement d'outils numériques permettant la simulation réaliste de ce procédé à l'échelle industrielle. La mécanique des fluides est simulée en utilisant une méthode d'Elements Finis stabilisés permettant de considérer des écoulements à haut nombre de Reynolds. Les transferts thermiques sont calculés directement sans l'utilisation de coefficients de transferts empiriques, en utilisant le couplage fort entre le fluide et le solide. Nous avons développé un modèle de changement de phase pour l'eau permettant de considérer les différents régimes d'ébullition. Une formulation unifiée des équations de Navier-Stokes, considérant une phase compressible et une phase incompressible a été développée afin de prendre en compte plus précisément la dynamique de la vapeur et de l'eau. Une procédure dynamique d'adaptation anisotrope de maillage, permettant une description plus fine des interfaces et une prise en compte plus précise des caractéristiques des écoulements est utilisée.Des exemples numériques exigeants ainsi qu'une validation expérimentale permettent d'évaluer la précision et la robustesse des outils proposés.Les outils développés permettent ainsi l'optimisation du mode opératoire du procédé, des ressources consommées et servent ainsi d'outils prospectifs pour la conception de produits. / Quenching processes of metals are widely adopted procedures in the industry, in particular automotive, nuclear and aerospace industries, since they have direct impacts on changing mechanical properties, controlling microstructure and releasing residual stresses of critical parts. Quenching is a highly nonlinear process because of the strong coupling between the fluid mechanics, heat transfer at the interface solid-fluid, phase transformation in the metal and boiling. In spite of the maturity and the popularity of numerical formulations, several involved mechanisms are still not well resolved.Therefore we propose a Direct Numerical Simulation of quenching processes at the industrial scale dealing with these phenomena. The fluid mechanics is simulated using a Finite Element Method adapted for high convective flows allowing the use of high stirring velocity in the quenching bath. Heat transfers are computed directly without the use of heat transfer coefficients but using the strong coupling between the fluid and the solid. We use a phase change model for the water that models all boiling regimes. A unified formulation of the Navier-Stokes equations, taking into account a compressible gas and an incompressible liquid is developed to model more accurately the vapor-water dynamics. A dynamic mesh adaptation procedure is used, increasing the resolution in the description of the interfaces and capturing more accurately the features of the flows.We assess the behavior and the accuracy of the proposed formulation in the simulation of time-dependent challenging numerical examples and experimental results.These recent developments enable the optimization of the process in terms of operating conditions, resources consumed and products conception.

La liquéfaction des sols saturés lors des séismes est l’un des problèmes les plus importants auxquels sont confrontés les ingénieurs. Il n’y a guère eu de séisme majeur sans au moins quelques cas de liquéfaction. Des tassements, des basculements de bâtiments, des écoulements latéraux, des cônes de liquéfaction et des instabilités de pentes, sont certaines de ses manifestations. La conception sismique des centrales nucléaires et autres installations critiques comprend systématiquement une évaluation du risque de liquéfaction.Dans ce cadre, des méthodes de nature entièrement empirique sont couramment utilisées en ingénierie. Ces approches procurent des marges à la conception et des limitations d’utilisation. Pour exploiter ces marges en situation de réévaluation, on a recours à des calculs transitoires non linéaires avancés dans lesquels on doit modéliser finement la loi de comportement du sol pour mettre en évidence les montées de pression interstitielle.Ces derniers calculs sont coûteux en termes de temps et de compétences numériques. L’objectif de ce travail de recherche est notamment de réduire les conservatismes en vigueur lors de l’utilisation de la méthode simplifiée sans pour autant mettre en œuvre d’emblée les méthodes les plus sophistiqués. On propose pour cela une méthode de complexité intermédiaire qui élargit l’applicabilité des modèles semi-empiriques pour une analyse plus fine du risque sismique.Dans un premier temps, en partant d’un calcul linéaire équivalent conventionnel, une nouvelle approche pour la prise en compte de la montée de pression interstitielle est proposée sous le nom « X-ELM ». Le modèle de comportement employé est basé sur la relation entre la pression interstitielle et les déformations volumiques plastiques. La nouvelle approche «X-ELM » est utilisée pour modéliser la réponse des sols pour le séisme de Tōhoku (Mw=9.0) à la ville d’Urayasu au Japon. Le modèle a été appliqué sur douze profils de sols différents. L’étude de ces cas rend possible la validation du modèle par comparaison des résultats des calculs aux observations in situ. Le modèle peut donc être considéré comme un outil fiable pour la prédiction de déclenchement de liquéfaction des sols saturés.Ensuite, un outil de prédiction rapide a été conçu en se basant sur des approximations de processus aléatoire, sur les propriétés mécaniques de base du sol et sur les caractéristiques du chargement sismique. Outre son ampleur, une caractéristique importante du signal sismique d’entrée est sa durée qui peut conduire à de fortes non linéarités et à un état de liquéfaction étendu. En considérant donc la durée de phase forte, le spectre de réponse, la fréquence propre du modèle et les caractéristiques de densification du sol, l’outil de prédiction proposé procure des estimations rapides du taux de montée de pression interstitielle et du tassement pré-liquéfaction sans devoir exécuter des calculs transitoires.Enfin, un modèle 2D de barrage est étudié, en examinant l’influence de la montée de pression interstitielle et celle des déformations de cisaillement sur la réponse sismique de l’ouvrage. Un calcul linéaire équivalent adapté aux situations bidimensionnelles est élaboré et le prédicteur est employé pour évaluer la montée de pression interstitielle. On compare les résultats de la simulation aux observations in situ, piézométriques et accélérométriques.En conclusion, ce travail de recherche fournit des méthodes et outils de calculs numériques performants et accessibles aux ingénieurs pour l’évaluation sismique des profils de sols et des ouvrages en terre tels que digues ou barrages / The liquefaction of saturated soils during earthquakes is one of the most important problems facing engineers. There has hardly been a major earthquake without at least some cases of liquefaction. Settlements, tilting of buildings, lateral flows, sand boilings and slope instabilities have been some of its manifestations. The seismic design of nuclear power plants and other critical facilities systematically includes a liquefaction risk assessment.In this context, fully empirical methods are commonly used in engineering. These approaches provide design margins and limitations of use. To exploit these margins in a re-evaluation situation, we use advanced nonlinear transient calculations in which the soil behavior must be finely modeled to highlight the pore-water pressure build-up. These last calculations are expensive in terms of time and numerical skills. The objective of this research work is to reduce the conservatisms in force when using the simplified method without necessarily implementing the most sophisticated methods from the outset. To this end, we propose a method of intermediate complexity that broadens the applicability of semi-empirical models for a more detailed analysis of seismic risk.First, starting from a conventional equivalent linear calculation, a new approach for taking into account excess pore pressure is proposed under the name "X-ELM". The behavioral model employed is based on the relationship between pore pressure and plastic volumetric deformations. The new "X-ELM" approach is used to model soil response in the city of Urayasu,Japan during the Tohoku earthquake (Mw = 9.0). The model has been applied to twelve different soil profiles. The study of these cases makes possible the validation of themodel by comparing the results of the calculations with the observations in situ. The model can therefore be considered as a reliable tool for the prediction of liquefaction triggering of saturated soils.Then, a prediction tool was designed based on random process approximations, the basic mechanical properties of the soil and the characteristics of the seismic loading. In addition to its magnitude, an important feature of the input seismic signal is its duration which can lead to strong nonlinearities and an extended liquefaction state. Considering the strong phase duration, the response spectrum, the natural frequency of the model and the characteristics of soil compaction, this tool provides fast estimations of the rate of pore pressure build-up and pre-liquefaction settlement without having to perform transient calculations.Finally, a 2D dam model is studied, by examining the influence of excess pore pressure and that of the shear strains on the seismic response of the structure. An equivalent linear computation adapted to two-dimensional situations is elaborated and the predictor is used to evaluate pore pressure increase. Simulation results are compared with in situ, piezometric and accelerometric observations.In conclusion, this research work provides methods and tools of numerical computation that are efficient and accessible to engineers for the seismic evaluation of soil profiles and earth structures such as dikes or dams

Avec la réduction de taille des composants microélectroniques, le monosiliciure de nickel (NiSi) a été largement utilisé dans les transistors CMOS (Complementary-Metal-Oxide-Semiconductor) en tant que contacts pour les sources, drains et grilles. Cependant, NiSi se dégrade lors du recuit à haute température. Il apparait essentiel d'étudier la séquence de formation de phases et la stabilité du monosilicide en présence d’éléments d'alliage. Les réactions à l'état solide entre des films de Ni allié en W et/ou Pt et un substrat de Si ont été étudiées principalement par diffraction des rayons X (DRX) in-situ et sonde atomique tomographique (SAT). L'analyse combinatoire de la réaction entre des films Ni avec différents composition gradients et le Si a permis de comprendre la séquence de formation. Les concentrations des éléments d'alliage (W et Pt) ont un impact significatif sur la séquence de formation des siliciures de Ni et sur la cinétique de formation des phases. Le mécanisme d'agglomération des films minces de NiSi a également été étudié au cours de cette thèse. Un film de Ni pur de 15 nm a été déposé sur un substrat de Si enrichi de multicouches de Si isotopique. Des analyses SAT ont été effectuées sur l'échantillon de Ni/Si (isotope) après un recuit à 600°C. En observant la distribution des isotopes de Si, un mécanisme d'agglomération de NiSi a été proposé selon que les isotopes de Si restent sous forme de multicouches ou qu'ils se mélangent. Cette étude rendue possible grâce à la capacité unique de la SAT d'observer les isotopes en 3D et à l'échelle atomique apporte une meilleure compréhension de l'agglomération de films minces poly-cristallin d'intermétallique / With the downscaling of microelectronic devices, NiSi has been widely used in complementary-metal-oxide-semiconductor (CMOS) transistors as contact on source, drain and gate. However, NiSi suffers from degradation upon annealing at high temperatures. Adding alloying elements is an effective method to increase the stability of nickel monosilicide but the formation sequence of Ni silicides is substantially modified. Therefore, the studies of the phase formation sequence and the stability of monosilicide are of great importance.The solid-state reactions between Ni films alloyed with W and/or Pt and Si substrates were studied mainly by in-situ X-ray diffraction (XRD) and atom probe tomography (APT) using combinatorial analysis of co-deposited gradient films. The phase sequence was monitored by in-situ XRD and APT was used to examine the silicides and reveal the redistribution of alloying elements. The content of alloying elements (W and Pt) has a large impact on the phase sequence of Ni silicides and the kinetics of phase formation. The basic agglomeration mechanism of NiSi thin films was studied. A 15nm pure Ni film was deposited on a Si substrate enriched with isotope multilayers. APT analyses were performed on the sample of Ni/Si (isotope) after an annealing at 600°C. By observing the distribution of Si isotopes (30Si, 29Si and 28Si), whether they maintain a multilayer structure or are mixed together, a mechanism of the agglomeration of NiSi was proposed. This was possible because of the unique capability of APT to observe isotopes in 3D at the atomic scale and it allows a better understanding and to control of the agglomeration of poly-crystalline compound thin films

L’objectif de cette thèse est d’utiliser la géochimie des spéléothèmes à travers une étude multi-proxy afin de reconstituer des changements environnementaux et climatiques dans le Sud-Ouest de l’Europe au cours de l’Holocène. Ce travail basé sur trois sites répartis entre l’Espagne et la France couvre les 14 000 dernières années. Cette région du sud-ouest de l’Europe, située à la confluence des hautes et des basses latitudes en fait une région particulièrement sensible à la fois à l’influence climatique de la région Nord-Atlantique mais également des régions Méditerranéenne et subtropicales Africaines et Asiatiques. Ce contexte particulier permet d’étudier idéalement l’impact des différents mécanismes climatiques agissant sur les modifications de son environnement. La répartition des sites d’étude a pour but de comparer les sites sous influence majeure Atlantique (Candamo, Nord Espagne), Méditerranéenne (Salamandre, SE France) ou sous influence mixte (Nerjà, sud Espagne). La compréhension des comportements géochimiques d’éléments clés, a constitué la première partie de cette étude. Après avoir établi un cadre temporel bien contraint (14C, U/Th), la géochimie élémentaire (éléments traces et terres rares) a été couplée à l’étude isotopique (87Sr/86Sr, d18O, d13C). Le rapport 87Sr/86Sr traditionnellement utilisé comme traceur de sources du Sr a été couplé à l’étude des variations élémentaires des éléments alcalins (Li, Rb, Cs) afin de déterminer des changements de sources des éléments (roche-mère ou de sol) dans les spéléothèmes. Le comportement de ces alcalins s’est également révélé lié à des modes de transports différents (particulaire, dissous ou liés à des complexes organiques) et l’utilisation du d13C et des terres rares (fractionnement terres rares légères et terres rares lourdes et anomalie en cerium) a permis de confirmer ou non ces hypothèses. La deuxième partie de ce travail a été d’appliquer ces mécanismes géochimiques à la détermination des processus hydrologiques et environnementaux et de les replacer dans le cadre climatique propre à chaque site d’étude. Ainsi, les éléments alcalins (Li, Rb, Cs) se sont révélés de bons indicateurs de phénomènes d'altération pouvant être liés à des changements climatiques (comme la transition entre le Bölling-Alleröd et le Younger-Dryas dans le site de la Salamandre), et/ou environnementaux (période de déforestation ou de mise en place de végétation à l’Holocène moyen dans le cas du site de Nerjà ou de Candamo). Ces événements pouvant donc se dérouler en contexte très différents (climat sec, déforestation, apports détritiques en climat humide par lessivage), le couplage à d’autres proxis est alors nécessaire pour préciser les conditions d'enrichissement de ces éléments alcalins. / The purpose of this thesis is to use speleothem geochemistry through a multi-proxy study to reconstruct environmental and climatic variations in southwestern Europe during the Holocene. This work, based on three sites distributed between Spain and France, covers the last 14,000 years. Due to its location at the confluence of the high and low latitudes, this region of south-western Europe is especially sensitive to both the influence of the North Atlantic climate but also the Mediterranean and subtropical African and Asian climate. This particular context allows to study ideally the impact of the different climatic mechanisms driving the modifications of its environment. The choice of the study sites’ distribution aims to compare sites with main Atlantic influence (Candamo, N Spain), Mediterranean (Salamandre cave, SE France) and mixt one (Nerjà, S Spain). The understanding of geochemical behavior was the first part of this study. After the establishment of a chronological frame (14C, U/Th), elemental geochemistry (trace elements and rare earth elements) and isotopic (87Sr/86Sr, d18O, d13C) was coupled. The 87Sr/86Sr isotope ratio, traditionally used as a source tracer for Sr, has been coupled to the study of elementary variations of alkali elements (Li, Rb, Cs) in order to determine changes in the contribution of bedrock and soil in speleothems. The behavior of these alkalis has also been shown to be related to different modes of transport (particulate dissolved or bound to organic complexes) and the use of d13C and rare earths (light rare earth and heavy rare earth fractionation and cerium anomaly) made it possible to confirm or not these hypotheses. The second part of this work was to apply these geochemical mechanisms to the determination of the hydrological and environmental processes and to place them in the specific climatic context of each study site. Thus, the alkali elements (Li, Rb, Cs) have proved to be good indicators of weathering phenomena that may be related to climate change (such as the transition between Bölling-Alleröd and Younger-Dryas in the Salamandre site), and/or environmental changes (period of deforestation or establishment of vegetation in the Middle Holocene in the case of the sites of Nerjà or Candamo). These events can thus take place in very different contexts (dry climate, deforestation, detrital inputs during wet climate by leaching) the coupling to other proxis is then necessary to specify the conditions of enrichment of these alkali elements.

De nombreux phénomènes physiques peuvent être étudiés au moyen de modélisations et de simulations numériques, courantes dans les applications scientifiques. Pour être calculable sur un ordinateur, des techniques de discrétisation appropriées doivent être considérées, conduisant souvent à un ensemble d’équations linéaires dont les caractéristiques dépendent des techniques de discrétisation. D’un côté, la méthode des éléments finis conduit généralement à des systèmes linéaires creux, tandis que les méthodes des éléments finis de frontière conduisent à des systèmes linéaires denses. La taille des systèmes linéaires en découlant dépend du domaine où le phénomène physique étudié se produit et tend à devenir de plus en plus grand à mesure que les performances des infrastructures informatiques augmentent. Pour des raisons de robustesse numérique, les techniques de solution basées sur la factorisation de la matrice associée au système linéaire sont la méthode de choix utilisée lorsqu’elle est abordable. A cet égard, les méthodes hiérarchiques basées sur de la compression de rang faible ont permis une importante réduction des ressources de calcul nécessaires pour la résolution de systèmes linéaires denses au cours des deux dernières décennies. Pour les systèmes linéaires creux, leur utilisation reste un défi qui a été étudié à la fois par la communauté des matrices hiérarchiques et la communauté des matrices creuses. D’une part, la communauté des matrices hiérarchiques a d’abord exploité la structure creuse du problème via l’utilisation de la dissection emboitée. Bien que cette approche bénéficie de la structure hiérarchique qui en résulte, elle n’est pas aussi efficace que les solveurs creux en ce qui concerne l’exploitation des zéros et la séparation structurelle des zéros et des non-zéros. D’autre part, la factorisation creuse est accomplie de telle sorte qu’elle aboutit à une séquence d’opérations plus petites et denses, ce qui incite les solveurs à utiliser cette propriété et à exploiter les techniques de compression des méthodes hiérarchiques afin de réduire le coût de calcul de ces opérations élémentaires. Néanmoins, la structure hiérarchique globale peut être perdue si la compression des méthodes hiérarchiques n’est utilisée que localement sur des sous-matrices denses. Nous passons en revue ici les principales techniques employées par ces deux communautés, en essayant de mettre en évidence leurs propriétés communes et leurs limites respectives, en mettant l’accent sur les études qui visent à combler l’écart qui les séparent. Partant de ces observations, nous proposons une classe d’algorithmes hiérarchiques basés sur l’analyse symbolique de la structure des facteurs d’une matrice creuse. Ces algorithmes s’appuient sur une information symbolique pour grouper les inconnues entre elles et construire une structure hiérarchique cohérente avec la disposition des non-zéros de la matrice. Nos méthodes s’appuient également sur la compression de rang faible pour réduire la consommation mémoire des sous-matrices les plus grandes ainsi que le temps que met le solveur à trouver une solution. Nous comparons également des techniques de renumérotation se fondant sur des propriétés géométriques ou topologiques. Enfin, nous ouvrons la discussion à un couplage entre la méthode des éléments finis et la méthode des éléments finis de frontière dans un cadre logiciel unique. / Many physical phenomena may be studied through modeling and numerical simulations, commonplace in scientific applications. To be tractable on a computer, appropriated discretization techniques must be considered, which often lead to a set of linear equations whose features depend on the discretization techniques. Among them, the Finite Element Method usually leads to sparse linear systems whereas the Boundary Element Method leads to dense linear systems. The size of the resulting linear systems depends on the domain where the studied physical phenomenon develops and tends to become larger and larger as the performance of the computer facilities increases. For the sake of numerical robustness, the solution techniques based on the factorization of the matrix associated with the linear system are the methods of choice when affordable. In that respect, hierarchical methods based on low-rank compression have allowed a drastic reduction of the computational requirements for the solution of dense linear systems over the last two decades. For sparse linear systems, their application remains a challenge which has been studied by both the community of hierarchical matrices and the community of sparse matrices. On the one hand, the first step taken by the community of hierarchical matrices most often takes advantage of the sparsity of the problem through the use of nested dissection. While this approach benefits from the hierarchical structure, it is not, however, as efficient as sparse solvers regarding the exploitation of zeros and the structural separation of zeros from non-zeros. On the other hand, sparse factorization is organized so as to lead to a sequence of smaller dense operations, enticing sparse solvers to use this property and exploit compression techniques from hierarchical methods in order to reduce the computational cost of these elementary operations. Nonetheless, the globally hierarchical structure may be lost if the compression of hierarchical methods is used only locally on dense submatrices. We here review the main techniques that have been employed by both those communities, trying to highlight their common properties and their respective limits with a special emphasis on studies that have aimed to bridge the gap between them. With these observations in mind, we propose a class of hierarchical algorithms based on the symbolic analysis of the structure of the factors of a sparse matrix. These algorithms rely on a symbolic information to cluster and construct a hierarchical structure coherent with the non-zero pattern of the matrix. Moreover, the resulting hierarchical matrix relies on low-rank compression for the reduction of the memory consumption of large submatrices as well as the time to solution of the solver. We also compare multiple ordering techniques based on geometrical or topological properties. Finally, we open the discussion to a coupling between the Finite Element Method and the Boundary Element Method in a unified computational framework.