L'ouvrage contient une modalité originale d'approche des objets fractals, tant du point de vue mathématique, que du point de vue physique. Par la rigueur de cette analyse, les concepts de fractal et de dimension fractale ont pu être correctement définis. Ainsi, de l'étude comparative des méthodes de calcul de la dimension fractale applicables aux surfaces de rupture, on relève comme méthodes les plus efficaces: la méthode Box Counting, la méthode des oscillations et la méthode de l'exposant de Hurst. En commençant par la méthode classique des îles (Slit Island), les applications de l'interférométrie sur la détermination de la dimension fractale des surfaces de rupture, dans le cas des éprouvettes Charpy pour l'acier XC 65 trempé et l'acier inox 316L, ont pu être mises en place. Cette nouvelle méthode peut mieux caractériser la topographie de la surface de rupture mesurée en interférométrie à l'aide de la dimension fractale. Pour réaliser le traitement des données obtenues à la suite de l'analyse de l'état des surfaces des éprouvettes basées sur l'interférométrie, il a été nécessaire d'élaborer un logiciel (Mesrug) permettant le traitement statistique d'une base de données complexes. Ce logiciel utilisé initialement en France et créé par le Professeur Maxence Bigerelle de l'Université de Technologies de Compiègne, a été développé et réactualisé pour le calcul des paramètres de rugosité de la section de rupture, pour les éprouvettes Charpy. Nos recherches ont, par ailleurs, montré que l'énergie à la rupture par un seul choc varie de façon inversement proportionnelle à la dimension fractale ; cette observation constituant une contribution majeure de l'ouvrage.

La propagation de fissures est le mécanisme fondamental responsable de la rupture catastrophique des matériaux fragiles. Celle-ci est décrite traditionnellement par la Mécanique Linéaire Elastique de la Rupture. Or, si ce cadre théorique apparait performant pour décrire des fissures lentes, il échoue largement à haute vitesse. En particulier, il ne permet pas de rendre compte des vitesses de rupture maximales observées expérimentalement, ni de la rugosité des faciès observée à haute vitesse. Pour explorer ces phénomènes, nous avons mis en place un dispositif expérimental qui permet d'étudier les mécanismes de rupture dans un matériau fragile modèle - nous avons opté pour le Plexiglas - sur une large gamme de vitesse, aux petites échelles d'espace et de temps. Ce dispositif nous a permis de mettre en évidence une nouvelle vitesse critique au delà de laquelle, la propagation de la fissure s'accompagne d'endommagements macroscopiques sous forme de nucléation et de croissance de microfissures en avant du front. Un scénario simple permet de prendre en compte cet endommagement pour quantifier la variation de l'énergie dite de fracture, i.e. l'énergie dissipée par le matériau lorsque la fissure se propage d'un incrément de surface, et expliquer la valeur anormalement basse de la vitesse limite de rupture observée dans les matériaux fragiles. Il explique aussi la nature grenue des faciès de rupture observés. Nous avons par ailleurs pu montrer qu'il était possible, à partir des faciès de rupture, de reconstruire de manière déterministe, à l'échelle du micromètre et de la microseconde la dynamique de propagation du front de fissure et le développement d'endommagement associé.

Les structures de sphères creuses appartiennent à la famille des matériaux cellulaires qui ont récemment été étudiés pour leurs multiples propriétés. Dans le cas de cette thèse, le but des sphères creuses est de dissiper l’énergie d’impact d’un oiseau sur un cockpit d’avion. Elles sont développées dans le cadre du projet SAMBA (Shock Absorber Material for Birdshield Application) afin d’optimiser leur énergie spécifique absorbée (J/kg).Dans un premier temps, des essais quasi-statiques (v = 5 mm/min) et dynamiques (v = 2 m/s) de compression uni-axiale sont conduits à température ambiante sur une seule sphère creuse de diamètre 30 mm. Une propagation rapide de fissures macroscopiques est observée. Le formalisme de la Mécanique Élastique Linéaire de la Rupture (MELR) est utilisé pour estimer le taux de restitution d’énergie critique dynamique GIdc du matériau constitutif. La position du sommet de fissure est mesurée pendant la propagation de fissure à l’aide d’une caméra rapide. La Méthode des Éléments Discrets (DEM) permet de simuler la rupture dynamique en implémentant une technique de relâchement des nœuds. Le taux de restitution d’énergie GIdc peut être estimé à partir de l’histoire (position et temps) du sommet de fissure. Le modèle numérique montre que les structures sphériques dissipent une proportion importante de l’énergie par des effets dynamiques. A une même vitesse de propagation, plus l’épaisseur de coque est fine, plus les effets inertiels générés par la rupture sont importants et ce pour une même vitesse de propagation.Le modèle numérique DEM est ensuite employé pour reproduire la rupture dynamique sur une sphère creuse à l’aide d’un critère en contrainte seule ou un critère mixte en contrainte – énergie. Les bons résultats obtenus démontrent la capacité de la DEM à représenter la propagation de fissures en régime dynamique.Finalement, des essais numériques et expérimentaux multi-sphères sont réalisés afin évaluer le comportement des sphères creuses au sein d’un assemblage. / Hollow sphere structure (HSS) belongs to cellular solids that have been studied recently for its multiples properties. In our case, HSS aims to absorb soft impacts energy on an airliner cockpit. HSS is investigated through the SAMBA (Shock Absorber Material for Bird-shield Application) project because of its promises in term of specific energy dissipated (J/kg) during impact.First of all, quasi-static and dynamic (v = 5 mm/min to v = 2 m/s) uniaxial compression tests are conducted at room temperature on a single sphere (D = 30 mm). Rapid crack propagation (RCP) is observed to be predominant at macroscopic scale. The formalism of Linear Elastic Fracture Mechanics (L.E.F.M.) is therefore used to estimate the dynamic energy release rate GIdc . The crack tip location is measured during the crack propagation using a high speed camera. The Discrete Element Method (DEM) is used to simulate the dynamic fracture by implementing the node release technique. The dynamic energy release rate can be determined using an experimentally measured crack history. In spherical structures the numerical results reveal a high proportion of energy dissipated through inertial effects as well as a dependence of the thickness of the hollow sphere over the range of 0.04 mm to 1.2 mm.The DEM model Is then employed to reproduce the RCP according to two failure criterions: a stress criterion and a coupled stress-energy criterion. It reveals to be an interesting way to model the mechanical behavior of brittle materials.Eventually, experimental and numerical multi-spheres tests are performed to evaluate the behavior of brittle hollow spheres within an assembly.

Cette thèse s’intéresse à la propagation dynamique de fissure dans les membranes élastomères du point de vue expérimental. Elle a pour but d’identifier les paramètres qui gouvernent la cinématique de ces fissures se propageant à grande vitesse, afin d’en prédire la trajectoire. Fondé sur l’utilisation conjointe d’une caméra à haute résolution et d’une caméra rapide, le dispositif expérimental permet de mesurer des champs à partir de la corrélation d’images au cours de la propagation de la fissure. Mis en pratique sur un polyuréthane, ce dispositif permet de retrouver la configuration de référence de l’éprouvette pendant la propagation de fissure, préalable indispensable à l’étude mécanique du problème. En plus des champs cinématiques, la densité d’énergie élastique et les contraintes sont évaluées grâce à une loi de comportement hypérélastique Les résultats de ces essais constituent une large base de données sur la rupture dynamique de membranes élastomères. La méthode permet de réaliser une analyse cinématique et énergétique de la propagation stationnaire et instationnaire, toujours dans la configuration de référence. La propagation supersonique est observée pour les hauts niveaux de déformation de la membrane. Enfin,ces observations permettent une discussion sur l’utilisation de l’approche énergétique de la rupture dynamique et de la pertinence des mesures de champs actuelles pour caractériser ce type de propagation de fissure. / This PhD thesis tackles the issue of dynamic fracture of elastomer membranes from an experimental point of view. It aims at providing some insight to predict the trajectories of high speed cracks under large strain. An experimental procedure involving high resolution and high speed cameras is developed in order to perform full-field measurements based on digital image correlation during crack propagation. Tested with a highly stretchable elastomer (polyurethane), this set-up permits to retrieve the material configurations of the sample all along crack growth, which is a crucial step toward a complete mechanical analysis of the problem. In addition to the kinematic fields,both strain energy density and stress fields are estimated thanks to a hyperelastic model, which is issued from mechanical characterisation of the material. Results of these experiments provide a comprehensive database on dynamic fracture of membranes. The method is designed to perform kinematic and energetic analyses of both steady and unsteady crack propagation in the reference configuration. Supersonic crack growth is observed for large prescribed deformation of the membranes. Finally, these observations lead to a discussion on the energetic approach in dynamic crack growth and the current applicability of full-field measurements to characterise dynamic crack growth in elastomers.

Ce mémoire de thèse présente l'étude de la propagation et de l'arrêt d'une fissure de clivage dans un acier de cuve REP. Une bonne compréhension des phénomènes en jeu nécessite de solides données expérimentales, ainsi qu'un outil de modélisation performant. La méthode X-FEM est implantée dans Cast3m, permettant de simuler efficacement la propagation de fissure par éléments finis. Deux techniques concernant l'actualisation des fonctions de niveau, ainsi qu'une intégration non conforme sont proposées. La campagne d'essais de propagation de fissure concerne trois géométries : des éprouvettes CT, et des anneaux comprimés en mode I et mixte. On relève les vitesses de propagation. Des fractographies désignent le clivage comme responsable de la ruine. Un modèle de propagation basé sur la contrainte principale en pointe de fissure est identifié. La contrainte critique de clivage dépend de la vitesse de sollicitation. Ce modèle permet de prédire précisément par la simulation numérique, le comportement de la fissure observé expérimentalement.

L’objectif de cette thèse est d’approfondir la connaissance sur la génération et la propagation des ondes à travers l’interface de contact afin de contrôler leurs effets sur le frottement, contrôler les instabilités et réduire les phénomènes d'usure. Le travail est organisé en deux parties principales. La première partie est axée sur une analyse non-linéaire par éléments finis en grandes transformations; la rupture dynamique à l’interface de contact avec frottement, qui sépare deux corps (isotropes et élastiques) en condition de pré-charge statique, est simulée. On définit comme étant une rupture une zone, initialement en adhérence, qui change son état en devenant glissante. Les propriétés des ruptures sont analysées pour une surface plane entre deux matériaux différents en fonction de l’énergie de nucléation; l’effet de la rugosité de surface est ensuite analysée. En outre, la rupture "cumulatives" dans les aspérités et les conditions pour le couplage ou le non-couplage entre les ondes qui se propagent dans les deux corps sont étudiées. Dans la deuxième partie, l'amorçage du glissement entre deux corps en contact constitués de matériaux différents et séparés par une interface avec frottement, est simulée. L’évolution, en fonction du temps, des forces globales (normale et tangentielle) a été mise en relation avec les phénomènes locaux qui se déroulent à l’interface. L'analyse montre comment les micro-slips à l'interface, en agissant comme ruptures localisées, déclenchent les macro-slips entre les deux corps. L'interaction entre les dynamiques locale du contact et globale du système a été aussi étudiée. Enfin, une analyse paramétrique est menée en fonction de plusieurs paramètres (loi de contact, coefficient de frottement, amortissement matériau, charge normale, etc.). Les résultats mettent en avant le rôle-clé des micro-slips et des précurseurs (propagations d'ondes détectables qui ont lieu pour des valeurs de la force tangentielle globale inférieures à la valeur prévue par la loi de frottement) dans l'amorçage du macro-slip entre les deux corps. En fonction de leur distribution et de leur intensité, l'évolution des forces de contact change en passant d'un comportement de type stick-slip à un glissement continu. La dynamique locale au contact (propagation des ondes et ruptures) a été reliée au comportement global du système (stick-slip, glissement continu, vibrations induites); l'effet des paramètres du contact et du système sur le transfert d'énergie vibrationnelle entre le contact glissant et le système a également été examinée: en fonction de leurs valeurs, on peut avoir différentes modalités d'excitation du système (par une distribution de micro-slips ou par des macro-slips) et différentes processus de propagation et dissipation d'énergie. Les résultats numériques obtenus dans les deux parties de la thèse sont cohérents avec les résultats expérimentaux de la littérature. / This thesis is addressed to the understanding of the mechanisms at the origin of the contact wave fields at frictional interfaces and its relationship with the local characteristics of the surfaces in contact, as well as with the global dynamics and macroscopic frictional behaviour of the system. The aim of this work is to provide insights on the generation and propagation of the waves through the contact both to avoid instabilities and to control their effect on friction. The work is organized in two main parts. The first part presents the development of a non-linear finite element analysis in large transformations of the dynamic rupture at the interface with contact friction separating two bodies (isotropic and elastic) without relative motion. A rupture is considered when an initially sticking zone shifts in sliding state. The properties of the obtained ruptures are analyzed for a flat interface between dissimilar materials in function of the nucleation energy; then the effect of the interface roughness is analyzed. The differentiated rupture inside the asperities and the conditions for coupling and uncoupling between the waves radiating in the two bodies have been also investigated. In the second part, the analysis deals with the sliding onset between two bodies in contact. The sliding between two bodies made of different isotropic elastic materials and separated by a frictional interface is simulated. The evolution along the time of the global normal and tangential forces is analyzed, relating it to the local phenomena occurring at the interface. This part tries to investigate how micro-slips at the interface, acting as distributed ruptures, trigger the macro-slips between the two bodies. The interaction between local and global dynamics is also studied. Finally a numerical parameter space study is carried out, as a function of several system parameters (contact law, friction coefficient, material damping, normal load, translational velocity and regularization time). The results show the key role of the micro-slips and precursors (detectable wave propagations that occur at tangential global force well below the critical value expected by the friction law) in triggering the macro-slip between the two bodies. Depending on their distribution and magnitude the evolution of the contact forces passes from stick-slip-like behaviour to continuous sliding. The local dynamics at the contact (wave and rupture propagation) is linked to the global behaviour of the system (stick-slip, continuous sliding, induced vibrations); the effect of the contact and system parameters on the transfer of vibrational energy between the sliding contact and the system is investigated. The numerical results obtained by the two parts of the work show a good agreement with experimental results in literature.

L’objectif de cette thèse a été d’utiliser une approche numérique pour étudier l’impact des dommages d'irradiation sur la microstructure du combustible nucléaire composé d'un mélange d’oxyde de plutonium et d’uranium (MOX). Cette approche, réalisée à l'échelle atomique, repose sur l'utilisation de potentiels empiriques développés dans la littérature pour l'oxyde mixte $(U,Pu)O_2$ dans l’approximation des interactions d’ions rigides.Une première étape a été une analyse critique des propriétés structurales, thermodynamiques et mécaniques prédites par 5 potentiels de la littérature. Les calculs de dynamique moléculaire ont été réalisés à l'aide du code LAMMPS, sur l'ensemble de la gamme de composition de $UO_2$ à $PuO_2$ et à des températures comprises entre 300 K et le point de fusion. Nous montrons que les potentiels les plus satisfaisants sont ceux développés par Cooper et Potashnikov. Ces deux potentiels reproduisent correctement la stabilité thermodynamique des phases ainsi que l'évolution en température des paramètres de maille et de la chaleur spécifique. Cependant le comportement mécanique des oxydes est différent selon le potentiel choisi. Tout d'abord, on note que les constantes élastiques et le facteur d'anisotropie obtenus par le potentiel Cooper sont plus fidèles aux recommandations de la littérature. Ensuite, nous montons que la propagation d'une fracture induit une transformation de phase en pointe de fissure avec le potentiel de Cooper alors qu'une fissuration fragile est observée avec le potentiel Potashnikov.Une seconde étape a été l'étude des dommages d'irradiation dans les oxydes mixtes en se limitant à l'utilisation des deux potentiels de Cooper et de Potashnikov. Des calculs de dynamique moléculaire de cascades de déplacement ont été réalisés à différentes énergies, température et compositions en plutonium. Ensuite, afin d'évaluer les dommages d'irradiation sur des temps plus longs que ceux accessibles en dynamique moléculaire, la méthode d’accumulation de défauts a été utilisée. Pour les deux potentiels, l’évolution des dommages primaires avec l’augmentation de la dose suit globalement les mêmes étapes que celles trouvées précédemment dans $UO_2$. Tout d’abord, les défauts ponctuels sont créés. Par la suite, ils se regroupent et forment de petites boucles de Frank, qui après se transforment en dislocations parfaites. Cependant, la cinétique de la recombinaison de défauts ponctuels est significativement plus lente avec le potentiel de Cooper ce qui conduit à la création de petites régions désordonnées pour les cascades d'énergie élevée. L’effet de la teneur en plutonium est analysé. Nous montrons en particulier que la densité de dislocations crée diminue lorsque la teneur en plutonium augmente.Bien que la dynamique moléculaire a été décrite comme un microscope moléculaire en raison de sa capacité à décrire avec précision des systèmes atomiques, elle présente un inconvénient majeur, celui lié aux temps de l’ordre de la femto-seconde nécessaires pour résoudre les vibrations atomiques. Cela limite le temps total de simulation approximativement quelques nanosecondes. Afin de stimuler les processus, tels que la diffusion de cations, un autre outil de calcul est nécessaire. Les techniques de Monte Carlo (KMC) atomiques peuvent simuler de plus longtemps que la dynamique moléculaire . Cependant, pour que KMC fonctionne avec précision, il est nécessaire de connaître a priori tous les mécanismes de transition entre les états possibles. Pour cette raison, la method de Monte Carlo cinétique adaptative (aKMC) est choisie pour surmonter ces limitations. Cette méthode détermine les états de transition disponibles pendant la simulation. De cette manière, elle entraîne le système dans des états imprévus via des mécanismes complexes. La puissance de cet outil s'est révélée efficace pour découvrir la recombinaison de cations sur de plus longues périodes de temps que la DM. / The objective of this doctoral research is to use a numerical approach to study the impact of irradiation damage on the microstructure of the mixed uranium-plutonium oxide fuel (MOX). This numerical approach comprises mainly the use of Molecular Dynamics (MD) using empirical potential. Several empirical potentials for $(U,Pu)O_2$ can be found in the literature. The results of these potentials can exhibit significant differences. For this reason an extensive assessment of the main empirical potential found in the literature had to be performed.Five empirical interatomic potentials were assessed in the approximation of rigid ions and pair interactions for the $(U_{1-y},Pu_y)O_2$ solid solution. Simulations were carried out on the structural, thermodynamics and mechanical properties over the full range of plutonium composition, meaning from pure $UO_2$ to pure $PuO_2$ and for temperatures ranging from 300 K up to the melting point. The best results are obtained by potentials referred as Cooper and Potashnikov. The first one reproduces more accurately recommendations for the thermodynamics and mechanical properties exhibiting ductile-like behaviour during crack propagation, whereas the second one gives brittle behaviour at low temperature.From our results from the empirical potentials assessment, we can move to the radiation damage using only two potentials (Cooper and Potashnikov). In order to know the main source of defect during irradiation, MD displacement cascades were simulated. This revealed the damage created due to varying projectile energies. In addition, the Frenkel pair accumulation method was chosen to investigate the dose effect. This method circumvents the highly computing time demanding accumulation of displacement cascade by directly creating their final states, i.e. mainly point defects. Overall, results obtained with both potentials show the same trend. However, kinetics of point defect recombination are significantly slower with Cooper potential implying creation of small disordered region with high energy displacement cascades. The evolution of the primary damage with increasing dose follows the same steps as those found previously in pure $UO_2$. First, point defects are created. Subsequently, they cluster and form small Frank loops, which in turn transform and grow into unfaulted loops. We demonstrate also that increasing temperatures accelerate the production of dislocations shifting their creation to lower doses. The effect of the plutonium content is also evidenced, especially with Cooper potential. It shows that the dislocation density decreases when the plutonium content increases.Although, MD has been described as a molecular microscope due to its ability to discribe accuratily systems of atoms, it has a large drawback that is the short time steps of the order of femto-seconds needed to resolve the atomic vibrations. This limits the time typically few microsecond. In order to invetigate processess, such as, cation diffusion and rare-event annihilation of defects after cascaces, another computational tool is required. Atomistic or object kinetic Monte Carlo (KMC) techniques can run for longer timescales than MD. However, for KMC to work accurately, all of the possible inter-state transition mechanisms and their associated rates need to be known a priori. For this reason, the adaptive kinetic Monte Carlo (AKMC) is chosen to overcome these limitations. This method determines the available transition states during simulation. In this way, it takes the system into unforeseen states via complex mechanisms. The power and range of this tool proved to be efficient to discover cation Frenkel pair recombination over a longer periods of time than MD.

Pour prémunir les structures spatiales d'impacts hyper-véloces, le secteur aérospatial est continuellement à la recherche de matériaux toujours plus performants. Dans cette optique, les verres métalliques massifs se présentent comme de potentiels éléments de blindages spatiaux. De récentes études ont mis en exergue une meilleure résistance à la pénétration de ces matériaux comparativement aux blindages actuels. Les impacts par lanceurs permettent d'étudier et caractériser le comportement sous chocs des matériaux. Cependant, les vitesses des projectiles se retrouvent actuellement limitées à 10 km/s, correspondant aux niveaux d'impacts hyper-véloces les plus modérés. Pour s'affranchir de cette limitation, il est possible de se tourner vers les lasers de puissance. Ces dispositifs permettent de générer des ondes de choc dont les niveaux de pression et de vitesse de déformation sont supérieurs aux lanceurs. Les travaux menés ont permis d'étudier et de caractériser le comportement et l'endommagement de plusieurs verres métalliques ternaires ZrCuAl sous choc laser. Plusieurs campagnes expérimentales ont été réalisées sur les installations du Laboratoire pour l'Utilisation des Lasers Intenses (LULI2000 et ELFIE). Une partie de l'équation d'état des nuances étudiées a été obtenue à la fois par choc laser et compression isentropique. Les processus d'endommagement, l'influence des vitesses de déformation et de composition sur la rupture ont été étudiés. Pour des régimes de vitesse de déformation supérieurs de trois ordres de grandeur à ceux disponibles dans la littérature, il a été mis en évidence que les verres métalliques étudiés présentaient une limite à la rupture cinq à dix fois supérieure. / Space industry is always searching for efficient materials to protect space structures from high-velocity impacts. In this context, bulk metallic glasses appear as suitable elements of space debris shielding assemblies. Recent studies revealed a higher tolerance to impact of metallic glasses compared to materials currently used in shield assemblies. Gas-gun and powder launchers are usually used to study and characterize the dynamic and shock behaviour of materials. However, projectiles velocities are currently limited to 10 km/s, corresponding to the lowest high-velocity impacts levels. To overcome this limitation, one may consider shock waves induced by high-power laser facilities, whose pressure and strain rate levels can exceed those induced by canons. Hence, this work enabled to study and to characterize the dynamic and damage behaviours of several compositions of ternary ZrCuAl bulk metallic glasses subjected to shock waves induced by laser irradiation. Several experimental campaigns have been conducted on various laser facilities of the Laboratoire pour l'Utilisation des Lasers Intenses (LULI2000 and ELFIE). A part of the equation of state of the studied compositions was established using both shock waves and quasi-isentropic compressions. Damage processes and the composition and strain rate effects on fracture were also studied. For strain rate levels higher of three or more orders of magnitude than those available in the literature, it was shown that studied bulk metallic glasses displayed a five to ten times higher dynamic tensile limit.

La thèse traite, à l'aide de modèles simples, quelques phénomènes présents en mécanique de la rupture en dynamique et en milieu hétérogène. Dans un premier temps, nous nous sommes intéressés aux phénomènes de propagation brutale sous chargement quasi-statique (on entend par chargement quasi-statique, lorsque la vitesse du chargement tend vers 0). Pour cela nous avons développé un modèle de décollement de film (soumis à un déplacement imposé) à travers une interface ayant une discontinuité de ténacité. Le film commence par se décoller de façon stable, puis lorsqu'il se décolle à travers la discontinuité de ténacité, le décollement devient brutal sur une certaine longueur. La solution dynamique du modèle a été développée et nous avons obtenu que lorsque le chargement du film devient quasi-statique alors la solution dynamique converge en tous points vers une solution quasi-statique particulière que l'on nomme "métastable". Dans un deuxième temps, nous avons cherché à construire le même résultat pour une structure 2D présentant une discontinuité de ténacité à travers le trajet de fissuration. Pour cela nous avons utilisé un algorithme de type déboutonnage de noeuds pour construire la solution dynamique. Nous avons obtenu que la longueur de la propagation brutale donnée par l'approche numérique est plus faible que celle obtenue avec la solution "métastable". Ces différentes considérations nous ont permis d'introduire quelques grandeurs couramment utilisées dans l'industrie comme la ténacité à l'arrêt.<br> Dans un troisième temps nous avons cherché à analyser l'influence d'un défaut sur la propagation d'une fissure. Pour cela nous avons utilisé le modèle de film précédemment défini et la discontinuité de ténacité a été remplacé par un défaut. On montre que le décollement du film est brutal à travers le défaut, redevient stable en sortie de défaut et converge (lorsque la vitesse du chargement tend vers zéro) vers une courbe limite que l'on nomme "couche limite". Ensuite après un temps suffisamment long, la couche limite regagne la solution quasi-statique. Nous avons montré que cette couche limite peut s'interpréter d'un point de vue énergétique par un mécanisme de relaxation de l'énergie cinétique. Différentes interprétations de ce phénomène en statique et en dynamique sont présentées et en paticulier des ténacités apparentes sont introduites. Dans un quatrième temps nous nous sommes intéressés à l'influence d'un milieu contenant une infinité de défauts. Le même modèle de film est utilisé avec une infinité de défauts. Différentes solutions quasi-statiques sont présentées, en particulier la solution métastable et celle obtenue avec une approche de type minimum global de l'énergie. En dynamique deux approches différentes sont étudiées: une première qui permet de construire la solution limite (lorsque la vitesse du chargement tend vers zéro) et une seconde plus directe et qui consiste à prendre en compte toutes les ondes. Pour la seconde méthode, un passage à la limite (vitesse de chargement) sera effectué de façon numérique. Au final on obtient que lorsque la taille des défauts tend vers zéro et que le chargement devient quasi-statique alors une ténacité équivalente du milieu apparait.

Durant un tremblement de terre, la rupture grandit et se propage sur la faille. Lorsqu'elle s'arrête, le séisme atteint sa taille finale. Comprendre ce qui, dans la nature, détermine la capacité à se propager ou à s'arrêter de la rupture est fondamental en sismologie, puisqu'il existe des séismes de toutes les tailles, mais que seuls les plus grands sont dévastateurs. Dans cette thèse, nous étudions l'impact de différentes façons d'arrêter ou de perturber la propagation de la rupture, à travers des études numériques dynamiques. <br /><br />Tout d'abord, nous avons inclus une limite à l'élasticité du milieu entourant la faille, afin de simuler la propagation de la rupture dans un milieu fracturé, qui dissipe une partie de l'énergie libérée. Nous avons, pour la première fois, inclus et étudié l'impact de cette dissipation dans un modèle de rupture 3D. La rupture, dans ces conditions, est beaucoup plus sensible aux barrières, et s'arrête plus facilement. La cinématique de la rupture est remarquablement modifiée (vitesse de propagation plus lente, vitesse de dislocation maximale limitée). Les mouvements engendrés en surface sont atténués. La plasticité est, de fait, un phénomène crucial à prendre en compte dans la modélisation de la rupture sismique.<br /><br />Dans une deuxième partie, nous avons étudié l'impact d'une hétérogénéité spatiale de résistance à la rupture sur la faille. L'introduction de l'hétérogénéité permet d'obtenir des profils de glissement dont la forme se rapproche des formes observées dans les cas naturels. La propagation et l'arrêt de la rupture perdent leur caractère prédictible lorsque ce sont de petites barrières qui arrêtent la rupture progressivement. On peut obtenir une grande variété de tailles d'événements pour une même statistique de taille des barrières. L'obtention d'une loi puissance de type Gutemberg-Richter sur toute la gamme des tailles d'événements est conditionnée par l'augmentation progressive de l'énergie de fracturation avec la taille de la rupture, d'une façon similaire à ce qui est obtenu en considérant un comportement plastique du milieu. <br /><br />Enfin, l'étude des relations glissement final - taille de l'aspérité rompue, dans des modèles lisses de type aspérité/barrière, a montré que la dynamique contrôlait une partie de la loi d'échelle du glissement maximum, et que la segmentation des failles modifie sensiblement la loi d'échelle.