L'objet de cette thèse est l'extension aux structures complexes de la méthode de Résolution Inverse Fenêtrée Filtrée (RIFF). L'idée principale se base sur le modèle Eléments Finis local et libre d'une partie de la structure étudiée. Tout d'abord, la méthode a été développée dans le cas des poutres. Les mesures de vibrations sont alors injectées dans le modèle Eléments Finis de la partie de poutre analysée. Les rotations sont estimées à l'aide de mesures de déplacements supplémentaires et des fonctions de forme sur le support élémentaire. La méthode étant sensible vis-à-vis des incertitudes de mesures, une régularisation a dû être développée. Celle-ci repose sur une double inversion de l'opérateur où une régularisation de type Tikhonov est appliquée dans la seconde inversion. L'optimisation de cette régularisation est réalisée par le principe de la courbe en L. A cause des effets de lissage induits par la régularisation, les moments ne peuvent être reconstruits mais ils apparaissent comme des ''doublets'' de forces. Ceci nous a conduit à résoudre le problème en supposant que seules des forces agissent sur la poutre. Enfin, une étude des effets de la troncature du domaine a été menée dans le but de s'affranchir des efforts de couplage apparaissant aux limites de la zone étudiée. Le cas des plaques a été considéré ensuite afin d'augmenter progressivement la complexité des modèles utilisés. L'approche Eléments Finis a permis d'intégrer à la méthode des techniques de condensation dynamique et de réduction par la méthode de Craig-Bampton. Le nombre de degrés de liberté est trop élevé pour permettre une estimation des rotations par mesures de déplacements supplémentaires, la condensation dynamique est employée afin de les supprimer dans le modèle théorique. Par ailleurs, la régularisation induisant une perte de résolution spatiale à cause de son effet de lissage, une procédure de déconvolution spatiale basée sur l'algorithme de Richardson-Lucy a été ajoutée en post traitement. Enfin, une application de la méthode à la détection de défauts a été envisagée de même que l'application de la méthode à l'identification des efforts appliqués par une pompe à huile sur un banc d'essais industriel. Le travail s'est donc appuyé sur des développements numériques et la méthode a été validée expérimentalement en laboratoire et en contexte industriel. Les résultats de la thèse fournissent un outil prédictif des efforts injectés par des sources de vibrations raccordées à une structure en s'appuyant sur un modèle Eléments Finis local et des mesures de déplacements vibratoires, le tout en régime harmonique. / The object of this thesis is the extension to complex structures of the RIFF method (Résolution Inverse Fenêtrée Filtrée). Considering a subpart of a structure, the main idea is to build a local Finite Element model using free boundary conditions. First, the general method was developed on beams. Vibration measurements are injected in the Finite Element model of the analysed part of the beam. Rotations are estimated using extra-displacement measurements and elementary shape functions. The method is highly sensitive towards errors present in measurements, so a regularisation had to be used. This one consists in a double inversion of the operator where a Tikhonov regularisation is applied when performing the second inversion. The regularisation parameter is tuned by the L-curve principle. Because of the smoothing effect of the Tikhonov procedure, moments cannot be reconstructed anymore at this stage, but they do still appear as sets of opposite forces. This setback led us to solve the problem by restricting it to forces only equations. At last, the study of the truncature of the domain was conducted in the aim to suppress coupling forces appearing at the limits of the studied area. Then, the case of plates was considered in order to increase progressively the models’ complexities. The Finite Element approach permitted us to implement dynamical condensation as well as Craig-Bampton reduction techniques. This allowed us to reduce the total number of degrees of freedom to be taken into account both from a numerical and an experimental standpoint. For example, dynamical condensation allows to eliminate rotations in the model. Besides, regularisation induces a lack of spatial resolution because of its smoothing effect. A spatial deconvolution technique was therefore developed; it is based on the Richardson-Lucy algorithm which is applied at a post-processing stage. At last, it was successfully proposed to extend the method to the application of detecting defaults present in the structure. The method was also validated on an industrial test bench in order to identify the forces applied by an oil pump taken from a truck’s engine. This phD thesis relied on numerical developments and the method was validated experimentally both in laboratory and industrial context. Main results provide a predictive tool to evaluate injected forces by vibration sources linked to a structure. It necessitates to inject vibratory displacements measurements into a Finite Element model.

La justification de la stabilité d’un barrage-poids existant sur une fondation rocheuse nécessite une justification de la stabilité vis-à-vis du glissement à l’interface béton-roche. Le comportement au cisaillement de cette interface dépend de l’effort normal appliqué, des déformabilités de la roche et du béton, de la morphologie de la surface rocheuse et éventuellement de la propreté avant le coulage. Dans la littérature, un effet d’échelle a été mis en évidence : les évaluations de la résistance au cisaillement réalisées expérimentalement sur de petites interfaces en laboratoire donnent des valeurs différentes de celles estimées par rétro-analyse sur barrages existants. Le travail de cette thèse porte sur l'évaluation de cet effet d'échelle en tenant compte des facteurs influents, en particulier la morphologie de la surface rocheuse, l'objectif étant d'améliorer l'évaluation de la résistance au cisaillement à la base d'un barrage-poids.Tout d’abord, des essais de cisaillement directs ont été effectués sur plusieurs éprouvettes d’une même combinaison de différentes dimensions : à petite échelle sur des carottes et à des échelles intermédiaire et métrique sur des éprouvettes parallélépipédiques. Ces essais ont été complétés par des essais de traction directe permettant de caractériser l’adhésion initiale du contact béton-roche. Les résultats ont effectivement permis d'observer une variabilité des mesures de résistance au cisaillement selon l'échelle d'observation. Les essais usuels sur des carottes ont donné les mesures les plus élevées. A l'échelle métrique, les mesures de la résistance au cisaillement de l’interface béton-roche ont indiqué une valeur de cohésion non négligeable. A l'échelle intermédiaire, les valeurs de la résistance se sont réparties en deux groupes correspondant à deux morphologies différentes de la surface rocheuse. Du point de vue de l'effet d'échelle sur la résistance au cisaillement, si on représente les résultats par un critère linéaire, il a été constaté un effet d'échelle positif sur l'angle de frottement et un effet d'échelle négatif sur la cohésion.Un outil de description de la morphologie de la surface rocheuse basé sur une décomposition en fonctions sinusoïdales a été développé. Il a permis de décrire objectivement les caractéristiques de la morphologie de la surface rocheuse propres à chacun des deux groupes identifiés à l'échelle intermédiaire. Le premier groupe de surfaces rocheuses, avec les valeurs de résistance les plus élevées, se caractérise par une longueur d'onde principale de l'ordre de la taille de l'échantillon. Pour le second groupe, où les valeurs de résistance sont plus faibles, la longueur d'onde est comprise entre la moitié et les deux tiers de la taille de l'échantillon. Suite à ces observations expérimentales, les hypothèses suivantes ont été posées pour effectuer l'analyse des résultats : les essais sur carotte à petite échelle permettent de caractériser la résistance d'une interface béton-roche sans effet de morphologie de surface et l'échelle intermédiaire correspond à la surface élémentaire qui permet de caractériser les principaux effets de la géométrie de cette surface rocheuse sur la résistance. Sur cette base, un modèle 2D aux éléments finis de l'essai de cisaillement direct a été développé sous le logiciel Abaqus. Ce modèle a permis de reproduire le comportement au cisaillement des interfaces béton-roche à l’échelle intermédiaire. Ainsi il a pu être identifié que selon la morphologie, différents mécanismes d'initiation et de propagation de la rupture interviennent localement le long de l'interface et contribuent à moduler sa valeur de résistance globale, observée à l'échelle de l'essai. A l'échelle métrique, une décomposition de la surface en surfaces élémentaires et l'analyse de la morphologie de ces surfaces ont permis de retrouver les valeurs de résistance évaluées expérimentalement. / The shear strength of concrete-rock interface is a key factor to justify the stability of an hydraulic structure foundation. The shear behavior of this interface depends on the applied normal force, the deformability of rock and concrete, the rock surface morphology and the cleanliness of the rock surface before concrete casting. A scale effect was highlighted in the literature: the assessments of shear strength in a laboratory by achieving tests on small interfaces give values different from those estimated by back-analysis on the existing dams. This thesis work aims to assess the scale effect on the shear strength of concrete-rock interface taking into account some related factors, especially the rock surface morphology. The objective is to improve the shear strength assessment at the base of a gravity dam.Firstly, direct shear tests were performed on several prepared samples having different sizes: at small scale on cored samples usually used in the engineering field to assess the shear strength and at medium and metric scales on parallelepiped samples. To complete these tests, direct tensile tests were achieved to characterize the bonding of the concrete-rock contact. The results allowed a variability of shear strength measurements to be observed according to the observation scale. The usual tests on cored sample gave the highest values with high standard deviation. At metric scale, in the range of normal stresses to which the gravity dam foundation is mostly subjected, the concrete-rock interface shear strength showed a significant cohesion value. At medium scale, the shear strength values are divided into two groups corresponding to two different morphologies of the rock surface. Regarding the scale effect on the shear strength, if the results are fitted in a linear criterion, it was found a positive effect on the friction angle and negative effect on the cohesion.A descriptive tool based on sinusoidal functions was developed to characterize the rock surface morphology. It allowed to objectively describe the appropriate characteristics of the rock surface morphology for both groups identified at medium scale. The rock surfaces of the first group, which has the highest values of shear strength, are characterized relative to the best fit plane by a geometry having a main wavelength in the order of the sample size. For the second group, which has the lowest values of shear strength, the obtained wavelength is between half and two-thirds of the sample size.After these experimental observations and in order to accomplish the results analysis, the following assumptions were set: the tests on cored samples at small scale allowed the concrete-rock shear strength to be characterized independently of the rock surface morphology effect and the medium scale corresponds to the elementary surface which allows characterizing the main effects of the rock surface geometry on the interface shear behavior. On this basis, a 2D finite elements model was developed for direct shear test on Abaqus. This model reproduced the shear behavior of the concrete-rock interfaces at medium scale and therefore the stress distribution along the interface during the shearing process was studied. Thus it was identified that, according to surface morphology, different initiation and local failure propagation occur along the interface and contribute to adjust the overall shear strength of the interface. At scale 1, a decomposition of the surface by elementary surfaces and the analysis of the morphology of these surfaces allowed to recover the shear strength values assessed experimentally.

La rupture macroscopique d’un matériau intervient généralement lorsque des micro-défauts coalescent, plutôt que par la propagation catastrophique d’une seule fissure. Il est donc souhaitable d’étudier des configurations de rupture où de multiples fissures interagissent. Les paires de fissures en-passant (EP), où deux fissures parallèles croissent l’une vers l’autre, sont particulièrement intéressantes d’un point de vue applicatif. Cette configuration de rupture se retrouve aussi bien dans des situations naturelles (os, dorsales océaniques,…) qu’industrielles (génie civil, pièces métalliques,…). Malgré la diversité de tailles et de matériaux dans lesquels ces fissures existent, leurs trajectoires ont une forme typique en crochet quasi-universelle dont l’origine, résultant de l’interaction fissure-fissure répulsive puis attractive, est mal comprise. En particulier, le comportement répulsif initial semble mettre à mal la mécanique élastique linéaire de la rupture (MELR). Dans cette thèse, nous avons d’abord étudié les fissures EP dans le cadre de la MELR. L’étude de l’angle initial de déviation et la simulation de trajectoires a montré contre toute attente que la MELR permet de reproduire qualitativement la forme en crochet. Prédire précisément certaines caractéristiques, comme l’intensité de la phase répulsive, nécessite plus de finesse au niveau de la représentation du comportement matériau. Nous avons ensuite utilisé un modèle par champ de phase pour enrichir le modèle matériau. Les nouvelles trajectoires simulées étant fortement influencées par la longueur caractéristique du champ de phase, il est possible d’obtenir un modèle plus juste quantitativement. Une perspective intéressante reste de relier cette longueur à la microstructure du matériau / Macroscopic failure of a material happens generally through the coalescence of micro-defects rather than the catastrophic propagation of a single crack. It is therefore advisable to study fracture problems in which many cracks interact. The case of en-passant crack pairs (EP-cracks), two parallel and offset cracks approaching each other by propagating through their inner tips, presents a marked interest as these cracks can be found in various natural (bones, oceanic rifts,..) or industrial (civil engineering,…) situations. Despite the large variety of scales and materials in which these cracks are observed, their trajectories present a remarkably self-similar hook-shape. This shape result from the crack-crack interaction, first repulsive before becoming attractive, and its origin is poorly understood. In particular, the initial repulsive behaviour seems to question the validity of linear elastic fracture mechanics (LEFM). In this thesis, we first studied EP-cracks in the LEFM framework. The study of the initial kink angle and the simulation of crack paths showed against all expectations that LEFM is able to reproduce qualitatively the hook-shaped paths. Precise predictions of specific characteristics, such as the magnitude of repulsion, requires a more refined model of the material behaviour. We then used a phase-field model to augment the material representation. As they are strongly influenced by the characteristic length scale of the phase-field, the new simulated trajectories indicate that it is possible to develop a more quantitatively correct model. An attractive prospect is to link this characteristic length to the material microstructure

La magnéto-tellurique (MT) (Cagniard 1953, Tikhonov 1950) est une technique d'exploration de la Terre basée sur des mesures de champs électromagnétiques (EM). Une source naturelle (non artificielle) harmonique en temps et située dans l'ionosphère (Weaver 1994) produit un champ EM régi par les équations de Maxwell. Les champs électromagnétiques sont enregistrés par plusieurs récepteurs placés sur la surface de la Terre. Ces mesures sont utilisées pour produire une image du sous-sol à partir d'un procédé d'inversion utilisant des méthodes numériques. Nous utilisons la méthode hp-FEM résultant d'une extension du travail de Demkowicz 2005. Nous avons développé un logiciel qui résout, pour la première fois, le problème MT avec des éléments finis auto-adaptatifs. La méthode hp-FEM permet des raffinements locaux, à la fois en taille h et en ordre p sur les éléments, ce qui est un avantage notoire puisque la combinaison de ces deux types de critères permet de mieux capter la présence de singularités, fournissant ainsi des erreurs de discrétisation faible. C'est donc une méthode très précise dont la convergence est exponentielle (Gui and Babuska 1986, Babuska and Guo 1996). En raison des défis d'implémentation encore non résolus (Demkowicz et al. 2002) et de la complexité technique des calculs hp-FEM en 3D, nous nous limitons, dans ce travail, à des calculs en 1D et 2D.Le domaine de calcul est tronqué par un matériau absorbant (Perfectly Matched Layer PML, Berenger 1994), qui est conçu pour s'adapter automatiquement aux propriétés physiques des matériaux. En particulier, il s'ajuste efficacement à l'interface air-sol, où le contraste entre la conductivité des matériaux atteint jusqu'à seize ordres de grandeur. Dans cette thèse, nous présentons également des résultats préliminaires pour la mise en place d'une technique dimensionnelle adaptative plus connue sous le nom de DAM (Dimensionally Adaptive Method (DAM)). Lorsque la distribution de la résistivité du sous-sol dépend de multiples variables spatiales, une analyse correcte de la dimensionnalité (Ledo 2005, Martí et al. 2009, Weaver and Agarwal 2000) rend parfois possible de considérer les différentes régions avec des dimensions spatiales différentes. Par exemple, il est parfois possible d’interpréter la distribution comme une formation unidimensionnelle plus quelques hétérogénéités en 2D (ou 3D). Basée sur cette interprétation, la DAM tire profit d’une telle situation. Ainsi, l'idée principale de cette méthode est d'effectuer l'adaptativité sur la dimension spatiale en commençant par un problème de faible dimension et en utilisant les résultats obtenus pour minimiser le coût des problèmes de dimension supérieure. Nous commençons l'inversion avec un modèle 1D. Les résultats de ce problème d'inversion 1D sont utilisés comme information a priori sur les modèles de dimension supérieure. Un avantage fondamental de cette approche est que nous pouvons utiliser les solutions des problèmes de dimension inférieure précédemment calculées comme composantes du terme de régularisation associé à un problème de dimension supérieure afin d'augmenter la robustesse de l'inversion. Cette thèse propose également une analyse numérique rigoureuse de divers aspects des problèmes MT. En particulier, nous avons: (a) étudié l'effet de la source, (b) effectué une analyse fréquentielle de sensibilité, (c) illustré l'augmentation du taux de convergence lorsque l'adaptativité hp est employée, (d) séparé les effets 1D et 2D dans la solution numérique et (e) exploré l'intérêt de considérer différentes variables pour effectuer l'inversion. / The magnetotelluric (MT) method is a passive exploration technique that aims at estimating the resistivity distribution of the Earth's subsurface, and therefore at providing an image of it. This process is divided into two different steps. The first one consists in recording the data. In a second step, recorded measurements are analyzed by employing numerical methods. This dissertation focuses in this second task. We provide a rigorous mathematical setting in the context of the Finite Element Method (FEM) that helps to understand the MT problem and its inversion process. In order to recover a map of the subsurface based on 2D MT measurements, we employ for the first time in Mts a multi-goal oriented self adaptive hp-Finite Element Method (FEM). We accurately solve both the full formulation as well as a secondary field formulation where the primary field is given by the solution of a 1D layered media. To truncate the computational domain, we design a Perfectly Matched Layer (PML) that automatically adapts to high-contrast material properties that appear within the subsurface and on the air-ground interface. For the inversion process, we develop a first step of a Dimensionally Adaptive Method (DAM) by considering the dimension of the problem as a variable in the inversion. Additionally, this dissertation supplies a rigorous numerical analysis for the forward and inverse problems. Regarding the forward modelization, we perform a frequency sensitivity analysis, we study the effect of the source, the convergence of the hp-adaptivity, or the effect of the PML in the computation of the electromagnetic fields and impedance. As far as the inversion is concerned, we study the impact of the selected variable for the inversion process, the different information that each mode provides,and the gains of the DAM approach.

Les procédés de traitement de surface tels que le grenaillage sont utilisés pour introduire des contraintes résiduelles de compression dans les matériaux métalliques. Cela permet de prévenir l'initiation et la propagation des fissures afin d'augmenter la durée de vie en fatigue des pièces mécaniques. Les limites et inconvénients des procédés actuels généralement utilisés pour cette fin sont connus et ont été mis en avant dans plusieurs études. Il s'agit d'une modification importante de l'état de surface (rugosité), une profondeur traitée limitée, des difficultés de mise en œuvre, etc. Aussi, de nouveaux procédés permettant d'obtenir des résultats équivalents voire meilleurs tout en évitant ces problèmes suscitent un intérêt grandissant. Le traitement de surface par cavitation (cavitation peening) et la mise en compression par impulsions électromagnétiques font partie de ces procédés innovants dont la modélisation a fait l'objet de cette thèse. Le traitement de surface par cavitation consiste à générer un jet d'eau submergé contenant des bulles de cavitation et dirigé vers la surface à traiter. La principale difficulté dans la modélisation du cavitation peening est la détermination du chargement mécanique sur la surface provenant de l'implosion des bulles de cavitation. L'effondrement des bulles peut se faire de manière sphérique ou non-sphérique suivant leur proximité de la paroi solide. Des modèles analytiques et numériques ont donc été élaborés pour étudier la dynamique des bulles et déterminer l'onde de pression due à une bulle sphérique d'une part, et la pression d'impact du micro-jet (bulle asphérique) d'autre part. Ces deux types de sollicitations ont été comparés et un modèle macroscopique pour le traitement par cavitation a été dérivé, en s'appuyant sur une simulation numérique du jet de cavitation. Une comparaison avec des résultats expérimentaux a été effectuée et a montré la pertinence du modèle proposé. En ce qui concerne le traitement par impulsions électromagnétiques, il consiste à générer un puissant champ magnétique transitoire par induction dans la pièce de matériau conducteur. Les contraintes résiduelles proviennent de la déformation plastique créée par les forces de Laplace à l'intérieur du matériau. Un modèle numérique couplant l'électromagnétisme et la mécanique des solides a été mis en place pour la simulation du procédé. Grâce à une étude d'optimisation, ce modèle a permis d'obtenir les paramètres électriques nécessaires pour arriver à mettre en compression des matériaux de type alliage de base nickel ou d'aluminium par cette méthode. Les profondeurs de la mise en compression calculées sont théoriquement par calcul plus élevées que celles obtenues avec les méthodes classiques de traitement de surface. Une étude de l'influence des différents paramètres a été faite et a montré que le procédé peut être contrôlé de manière relativement aisée en variant l'intensité et la pseudo-fréquence du courant traversant l'inducteur. / Surface treatments methods like shot peening are used to introduce compressive residual stresses in metallic materials. The compressive stresses prevent the initiation and growth of cracks and hence improve the fatigue life of mechanical parts. The drawbacks and limitations of the existing processes generally used for this purpose are known and have been highlighted in many studies. These are, among others, an important surface modification (roughness), a limited compressive depth, difficulties in execution, debris and contamination problems, etc. Therefore, the interest in new surface treatment methods, which permit to obtain equivalent or even better compressive results while avoiding the previous problems, are growing. Cavitation peening and electromagnetic pulse peening are part of these innovative processes which modeling is the aim of this PhD thesis. Cavitation peening is a process of surface treatment which acts by the generation of cavitation bubbles near the workpiece surface. The modeling of this process is challenging because of the complexity of cavitation phenomenon and the main problematic point is the determination of the mechanical loading on the material due the bubbles collapse. In this thesis, an approach of modeling for cavitation peening based on the study of the dynamics of cavitation bubbles is proposed. Spherical and aspherical collapse of bubbles near a solid surface are studied by some numerical and analytical models. These two sources of loading pressures have been compared and a macroscopic model for cavitation peening have been derived by associating the numerical simulation of the cavitation jet and the localization of the cavitation zone. The comparison between the final residual stress profile calculated with the proposed model and the experimental results were satisfactory. Electromagnetic pulse peening (EMP) is a contactless process of surface treatment which could be used to introduce compressive residual stresses in conductive materials, by the generation of a high transient electromagnetic field. Laplace forces induced in the material by magnetic induction are the source of the material plastic deformation and compressive residual stresses introduction. To predict the EMP results, a numerical model have been built for the process simulation. The model, based on finite element method, coupled successfully electromagnetic and mechanical phenomena by using a sequential-coupled approach. It was proven theoretically, by the study of a reference case, that compressive residual stresses could be induced in metallic materials like nickel-based (Inconel) or aluminum alloys by the means of the present process. It was also shown by the calculations that much higher compressive depths than those of conventional peening processes could be achieved. The parametric study exhibited the influence of the maximum current intensity and frequency which affect both the compressive depth and the maximum residual stress.

Les equations de Maxwell se resolvent aisement lorsque le domaine d'etude est regulier, mais lorsqu'il existe des singularites geometriques (coins rentrants en 2D, coins et aretes rentrants en 3D), le champ electromagnetique est localement non borne au voisinage de ces singularites. Nous nous interessons a la resolution des equations de Maxwell dans des domaines bornes, singuliers, a l'aide de methodes d'elements finis continus. En pratique, cela permet de modeliser des instruments de telecommunication comme les guides d'onde, les filtres a stubs. Nous analysons tout d'abord le probleme quasi-electrostatique 2D, afin de maitriser la discretisation en espace. Nous presentons trois methodes de calcul (formulations augmentees mixtes) qui donnent des resultats numeriques tres convaincants : - Une version epuree de la methode du complement singulier (conditions aux limites essentielles). - La methode de regularisation a poids : on introduit un poids qui depend des distances aux singularites geometriques (conditions aux limites essentielles). - La methode avec conditions aux limites naturelles. Nous etudions ensuite la generalisation de ces methodes aux domaines 3D. Nous detaillons la resolution des equations de Maxwell instationnaires en domaines singuliers 3D par la methode de regularisation a poids, et nous donnons des resultats numeriques inedits.

Nous modelisons les processus de glissement instable, en tenant compte de la geometrie souvent complexe des systemes de failles et des lois de frottement observees en laboratoire. Nous assimilons les failles a des interfaces de discontinuite dans un milieu elastique parfait, la croute, et utilisons une loi de frottement dependant du glissement deduite d'experiences qui montrent que la resistance du materiau diminue du seuil statique au seuil dynamique, proportionnellement au glissement. <br /><br />Nous proposons deux methodes numeriques. La premiere a pour but de simuler l'evolution temporelle spontanee d'un reseau de faille, soumis initialement a un champ de contraintes donne et auquel nous appliquons une perturbation initiale en vitesse. Elle utilise un schema de type Newmark en temps, et une discretisation spatiale en elements finis avec decomposition de domaine. Elle se revele capable de capturer efficacement les instabilites du glissement, et en particulier la phase d'initiation, qui precede la propagation de la rupture dynamique et qui se caracterise par une forme auto-similaire et une croissance exponentielle du glissement au cours du temps. Des experiences numeriques montrent que l'interaction, sur des segments de faille presentant un recouvrement significatif, se manifeste par l'existence de « zones d'ombre » dans lesquelles les contraintes sont dechargees et le glissement inhibe. En cas de recouvrement important de deux segments de faille, on observe une dissymetrie des profils de glissement, correspondant a la disparition de la singularite de contraintes a l'une des pointes de faille.<br /><br />Le deuxieme schema numerique realise l'analyse spectrale non lineaire du probleme de l'initiation « pseudo-linearise » autour de la position d'équilibre dans laquelle le système de failles est initialement au seuil de resistance statique. La non-linearite du probleme provient de la prise en compte des zones d'ombre dont on ne connait pas la geometrie a priori. Cette analyse permet de trouver le mode qui porte la signature de l'initiation, c'est-a-dire la forme auto-similaire remarquee plus haut. La version statique de cette analyse modale fournit un critere de stabilite des reseaux de failles, c'est-a-dire la valeur limite du taux d'affaiblissement au-dela de laquelle un episode de glissement donnera lieu a un evenement sismique.<br /><br />Nous faisons l'hypothese que le mode non lineaire statique, qui caracterise un comportement en affaiblissement a la limite de la stabilite, peut etre utilise pour decrire le glissement cumule a l'echelle tectonique sur un reseau de failles normales particulier en Afar, dont nous connaissons les glissements mesures en surface. Nous montrons qu'un choix judicieux du profil d'affaiblissement « equivalent » a l'echelle tectonique permet un bon accord entre le glissement observe et le mode statique. Nous en tirons des conclusions en termes d'interaction, mais aussi de propagation et/ou branchement des segments de failles. Enfin, nous decrivons deux autres applications : <br /><br />1) l'influence de la fracturation secondaire (endommagement) sur la forme du glissement<br />2) les parametres geometriques favorables a l'apparition d'une zone de relai entre deux segments de faille se propageant l'un vers l'autre.

La réaction alcali-silice est une pathologie des bétons qui se caractérise par une réaction chimique entre la silice réactive contenue dans certains granulats et les ions contenus dans la solution interstitielle. Elle entraîne la formation de gels silico-calco-alcalins, le gonflement voire la chute des caractéristiques des bétons affectés. Ce travail a pour but de fournir des outils d'aide à la compréhension et à l'expertise des ouvrages atteints. Il est composé de trois parties complémentaires. La première vise à prédire le gonflement du béton à partir de la connaissance des mécanismes réactionnels locaux. La seconde modélise l'influence de l'eau sur la réaction. La dernière vise à reproduire les effets structuraux induits par la réaction (endommagement, déplacements, chute des caractéristiques mécaniques...).

Ce travail concerne le développement d'une plateforme de prototypage virtuel dédiée à la simulation numérique des problèmes en grandes transformations thermomécaniques rapides. Dans ces situations, le matériau est soumis à de grandes déformations, de forts taux de déformations et une forte élévation de température résultant principalement de la déformation plastique. Cette sorte de sollicitations, pouvant être très rapides et fortement non linéaire, est généralement rencontrée lors de simulations numériques de problèmes d'usinage ou d'impact et de crash. Deux voies de recherche sont explorées dans ce travail: la première concerne l'aspect expérimental alors que la seconde se focalise sur le développement du code de calcul éléments finis dynamique DynELA. Concernant l'approche expérimentale, l'objectif principal est de développer un ensemble d'outils (dispositif expérimental et logiciels utilitaires dédiés) dans le but d'identifier les paramètres constitutifs pour divers matériaux soumis à des grandes déformations et de forts taux de déformation. Concernant l'aspect numérique, l'objectif est de développer un code de calcul par éléments finis explicite basé sur une approche de Programmation Orientée-Objets en C++. Ce travail a été réalisé au sein du Laboratoire Génie de Production (LGP) de l'Ecole Nationale d'Ingénieurs de Tarbes (ENIT).

L'objet de ce travail de these est l'etude des effets thermiques dus a l'interaction laser-matiere dans le cas de materiaux metalliques irradies avec des impulsions femtosecondes. Deux approches sont developpees. La premiere utilise un traitement numerique a deux dimensions du Modele a Deux Temperatures (MDT). Apres definition d'une Zone Affectee Thermiquement (ZAT) sur la base d'arguments de physique des materiaux, l'etendue de celle-ci est estimee par simulation numerique dans la direction perpendiculaire au faisceau laser. La deuxieme approche est de type experimentale. Les techniques de microscopie electronique en transmission et a balayage sont utilisées pour mesurer l'evolution des microstructures en bordure de la zone irradiee. Les regimes femtoseconde et nanoseconde sont compares dans l'ensemble des resultats de ce travail.