Dans ce travail de recherche, on aborde une étude expérimentale et numérique du comportement mécanique des roches argileuses soumis à des chargements mécaniques et à des sollicitations thermiques. L’esprit de cette étude provient de la pratique de l’exploitation des huiles lourdes avec la technique d’injection de vapeur à haute température où les roches de réservoir sont soumises à des sollicitations thermiques et hydromécaniques couplées. L’enjeu est d’étudier le comportement hydromécanique de ces matériaux soumis à des variations importantes de température afin d’évaluer la stabilité mécanique des réservoirs. L’étude expérimentale contient des modifications d'une cellule triaxiale autonome et auto compensé à haute température (250 C°) ainsi que le système de pilotage, le système de mesure et d’étalonnage de déformations. Ces modifications sont importantes pour effectuer des tests hydrostatiques, uniaxiaux et triaxiaux servaient à obtenir une base des données expérimentales, cette base caractérise l’effet thermique sur le comportement mécanique des roches argileuse.Le cadre général de la modélisation est d’abord proposé pour décrire le comportement mécanique d’argilite dans le cas isotrope. Après une analyse détaillée des données expérimentales, un modèle spécifique élastoplastique couplé à l’endommagement est élaboré pour décrire le comportement mécanique. Ensuite l’effet de la température est pris en compte. Les comparaisons entre les simulations numériques et les données expérimentales ont montré la capacité du modèle proposé pour la description du couplage hydromécanique et thermique. Afin de décrire le comportement des roches anisotrope, nous avons proposé une extension du modèle en y introduisant une formulation de tenseur de fabrique Cette formulation est exprimée en termes d'invariants couplé aux tenseurs de contraintes et d’orientation de chargement. Des essais en laboratoire sous différents chemins de sollicitations ont été modélisés, le modèle proposé semble décrire correctement les principales réponses mécaniques des matériaux. / We proposed, in this work, an experimental and numerical study of mechanical behavior of shale rocks subjected to mechanical and thermal loads.In the petroleum industry, during the production of heavy oil with the technique of steam water injection at high temperature, the cap rocks are subjected to coupled thermal and hydro-mechanical solicitations. The challenge is to study the hydro-mechanical behavior of these materials subject to large variations in temperature in order to assess the mechanical stability of the reservoir. The experimental study includes the modifications in a triaxial cell in ordre to support a high temperature (250° C). These modifications are very important for hydrostatic, uniaxial and triaxial tests, all these tests are used to obtain an experimental data base characterizing the thermal effect on the mechanical behavior of shale rocks.The modeling framework is proposed at first to describe the mechanical behavior of shale rock in isotropic case. After a detailed analysis of experimental data obtained in the experimental section, a specific coupled elastoplastic-damage model has been developed to describe the mechanical behavior of these shale materials. The effect of temperature is taken into account and a comparison between numerical simulations and experimental data have shown the ability of the proposed model for the description of thermo mechanical coupling. To describe the behavior of anisotropic rocks, we have proposed an extension of the fabric tensor model to present the initial anisotropy of shale rock. This formulation is expressed in terms of invariant stress tensor coupled with loading orientation. Laboratory tests under different stress paths were modeled, the proposed model seems able to describe correctly the main mechanical responses of shale materials.

Dans cette thèse, nous proposons une modélisation unifiée du couplage hydromécanique des massifs rocheux fracturés en se basant sur la caractérisation de réseaux de fractures discrètes (DFN). En particulier, nous proposons une extension de la méthode dite de ressorts de corps rigides (RBSM). Les principales contributions de ce travail portent sur les points suivants:(1) Un modèle de comportement mésoscopique, basé sur la méthode RBSM, est développé pour décrire la fissuration de la roche. Ce modèle permet de prendre en compte la distribution de fractures préexistantes, l'initiation, la propagation et la coalescence de nouvelles fissures. La roche intacte est représentée par un certain nombre de blocs rigides, répartis aléatoirement et reliés par des ressorts. La déformation macroscopique du massif est due à la déformation locale des interfaces entre blocs.(2) Nous avons ensuite proposé un modèle à double porosité pour modéliser le transfert hydraulique dans le massif à partir de la distribution de fractures discrètes. En combinaison avec le modèle RBSM, la variation de la perméabilité lors de la fissuration du massif est étudiée.(3) La formulation d’inégalité variationnelle pour le problème d’écoulement avec surface libre est étendue aux réseaux de fractures en 2D et 3D. Cette formulation permet d’aborder les problèmes de surfaces libres dans les massifs fracturés même des topologies et des frontières complexes.(4) Le comportement hydromécanique des massifs avec une grande densité de fractures est étudié. Un modèle de comportement est proposé pour décrire les relations non linéaires entre la contrainte normale et la déformation normale, entre le glissement tangentiel et la dilatation normale. / In this thesis, we have proposed a unified approach, based on the improved rigid body spring method (RBSM), for modeling the hydromechanical coupling in both intact and fractured rocks. The main contributions of the thesis are summarized as follows.(1) A mesoscopic model, based on the improved RBSM, is developed to account for the initial distribution of preexisting fractures, the whole process of initiation, propagation and coalescence of new cracks. The macroscopic deformation is due to the local deformation of interfaces between blocks.(2) A dual porosity model for seepage flow in rock mass is proposed, based on the discrete fracture network model and the permeability variation during rock fracturing is investigated.(3) The Variational Inequality (VI) formulation for seepage problem with free surface is extended to fracture networks both in 2D and 3D conditions. This formulation is in particular suitable to solve free surface problems in fracture networks even with very complex topology and boundaries.(4) By introducing a non-linear constitutive model for fractures, a hydromechanical coupling model is established for fracture networks. The governing mechanisms of hydromechanical coupling under different stress conditions are detected.

Dans l'industrie du pétrole et du gaz, on a estimé que plus de 500 millions de dollars US sont<br />perdus tous les ans, en raison des instabilités de puits (Chen, Tan et Haberfield, 2000, [7]).<br />Afin de limiter les problèmes d'instabilités relatifs à des excavations souterraines la prévision<br />de l'évolution des contraintes et des déformations à proximité de l'ouverture est essentielle.<br />Les essais sur des cylindres creux à parois épaisses sont un moyen facile, économique et<br />réaliste pour mieux comprendre les mécanismes de rupture de puits pétroliers. La première<br />partie de la présente recherche consistait à créer et développer un dispositif expérimental pour<br />réaliser des essais sur cylindres creux à parois épaisses sous deux conditions distinctes :<br />1) Essais sans fluide interstitiel (pour modéliser l'ouverture souterraine pendant et après le<br />creusement).<br />2) Essais avec écoulement radial convergent de fluide interstitiel (pour modéliser le puit<br />de pétrole pendant sa phase de production).<br />C'est un des avantages de ce dispositif expérimental puisque dans la majorité des cas,<br />cette possibilité d'écoulement radial n'existe pas dans les études rapportées dans la<br />littérature .La seconde partie de cette recherche est une campagne d'essais expérimentaux réalisée sur le<br />matériau artificiel appelé CPIR09. Son comportement est similaire à certain grès très poreux<br />qui constituent des réservoirs pétroliers.<br />Le mode de rupture observé dans les essais est une paire de plans de rupture radiaux<br />diamétralement opposés, parallèles à l'axe du trou et perpendiculaires à la paroi interne) qui<br />correspond à un mécanisme compactant. Les mesures montrent que les valeurs de la<br />perméabilité globale diminuent avec l'augmentation de pression externe. L'écoulement radial<br />convergent de fluide interstitiel n'a pas effet sur la forme de la rupture mais cet écoulement<br />enlève des particules de matière fracturée, ce qui crée une érosion à l'intérieure des bandes de<br />localisation. L'écoulement radial convergent de fluide interstitiel (en raison du différentiel de<br />pression interstitielle) peut légèrement diminuer la résistance des échantillons.

L'objectif de cette thèse est de caractériser par des approches expérimentales et numériques les effets hydrique, mécanique et hydromécanique dus à la désaturation et la resaturation de l'argilite de l'Est, roche hôte du futur laboratoire souterrain de stockage des déchets radioactifs à Bure (France). Nous présentons des approches expérimentales et numériques pour caractériser le phénomène de transfert hydrique dans l'argilite. Une méthode simple d'identification est proposée, en déterminant la diffusivité hydrique linéarisée à partir de mesures pondérales sur des échantillons (tube mince et plaque mince) soumis à des paliers d'humidité suivant un cycle de désaturation-resaturation. Le transfert hydrique est non-linéaire. Pour interpréter ce phénomène, on établit un modèle numérique non linéaire tenant compte des phénomènes physiques (la conduction hydraulique, la diffusion de vapeur, le changement de phase etc.). On étudie ensuite l'évolution du comportement physique et mécanique de la roche argileuse en fonction de l'humidité imposée, suivant un cycle de désaturation-resaturation par paliers successifs. La déformation hydrique, la vitesse de propagation des ondes ultrasoniques, les propriétés élastiques, les caractéristiques de rupture et les phénomènes différés dépendent de l'état hydrique de matériau. L'influence de la désaturation et de la resaturation sur "un modèle réduit de tunnel» est étudiée. Pour mettre en valeur l'anisotropie de la roche, on utilise des tubes épais dont l'axe est parallèle ou perpendiculaire à la stratification. Ces tubes sont soumis à des chargements hydriques par ventilation d'air à hygrométrie variable dans leur trou central. Nous utilisons un modèle poroélastique non-linéaire, pour interpréter les phénomènes hydromécaniques anisotropes observés. On montre que les chargements hydriques peuvent conduire à des ruptures de l'éprouvette qui respectent les directions d'anisotropie de la roche, et qui peuvent être interprétées par le modèle hydromécanique comme une violation d'un critère de rupture en contrainte totale de traction. Finalement, des calculs numériques pour les phénomènes engendrés par la désaturation et la resaturation à l'échelle des ouvrages souterrains sont réalisés, par des modèles numériques calibrés à partir de résultats expérimentaux. Des solutions analytiques et numériques par la méthode des éléments finis sont réalisées pour prévoir l'influence de la désaturation et de la resaturation sur le massif d'argilite dans les galeries souterraines.

Ce travail de thèse porte sur l'analyse des glissements de terrain, qui constituent un risque naturel majeur responsable de pertes humaines élevées ainsi que de grands dommages aux structures, des infrastructures et de l'environnement naturel. Cette question revêt une importance particulière, en raison du changement climatique, qui augmente le risque de fortes pluies ainsi que la sécheresse, et par conséquent le risque d'instabilité des pentes due à l'environnement. Généralement, l'analyse de la stabilité des talus est réalisée en utilisant l’approche de l'équilibre limite. Comme cette théorie ne tient pas compte du processus de mobilisation du frottement, elle pourrait conduire à une surestimation du facteur de sécurité. Une analyse fiable de la stabilité des talus, en particulier dans les sols hétérogènes soumis à l'action de l'eau, nécessite l'utilisation de méthodes numériques avancées. Deux méthodes ont été utilisées dans cette recherche: la méthode de couplage hydromécanique et la méthode dynamique non linéaire. / This research concerns analysis of landslides, which constitute a major natural risk responsible for high human losses as well as large damages to structures, infrastructure and natural environment. This issue becomes particularly important, because of the climate change, which increases the risk of heavy rains as well as severe drought and consequently the risk of slope instability due to the environment change. Generally, analysis of slope stability is conducted using the limit equilibrium theory. As this theory does not take into consideration the process of mobilization of the friction, it could lead to an overestimation of the safety factor. A reliable analysis of the slope stability, in particular in heterogeneous soils submitted to the water action, requires the use of advanced numerical methods. Two methods were used in this research: the coupled hydro-mechanical method and the nonlinear dynamic method.

L'étanchéité composite géomembrane(GM)-géosynthétique bentonitique (GSB) mise en place dans les barrières de fond d'installations de stockage de déchets non dangereux (ISDnD) peut être sujette à des transferts advectifs liés à l'existence de défauts dans la GM. Les lixiviats peuvent percoler dans le GSB, pénétrer dans le sol et les nappes phréatiques sous-jacentes ce qui peut nuire à l'environnement. Il est donc important de comprendre les mécanismes de transferts dans les étanchéités composites GM-GSB et de les quantifier afin de connaitre, maitriser et minimiser l'impact des transferts advectifs à travers la barrière vers l'environnement. Cependant, l'inaccessibilité de la GM rend difficile l'estimation réelle des fuites à travers l'étanchéité composite. La présente thèse évalue via une démarche expérimentale et numérique les transferts advectifs à travers les étanchéités composites et contribue à améliorer la compréhension des mécanismes de transfert en fonction des sollicitations extérieures. Le but est de bien cerner la problématique des transferts advectifs à travers les étanchéités composites GM-GSB, combler les vides des précédentes études et mettre en évidence les principaux paramètres à prendre en compte (contrainte de confinement, hétérogénéité du GSB, qualité de contact à l'interface GM-GSB, altération chimique et physique des GSB durant leur durée de service sur site). Leur influence sur l'étanchéité de l'ensemble et sur l'évolution des caractéristiques des matériaux utilisés est étudiée. L'approche expérimentale a consisté à acquérir des données sur la rétention d'eau et à introduire la notion des surfaces d'états des GSB sous le confinement généré par les déchets. Il a été démontré que le confinement réduit le gonflement du GSB lors de son hydratation ce qui permet de diminuer sa conductivité hydraulique à saturation. Ce résultat renforce les recommandations appelant à confiner rapidement les GSB après leur mise en place afin qu'ils acquièrent rapidement toutes leurs capacités d'étanchéité. Cette étude expérimentale a été complétée par une quantification numérique des transferts à travers les étanchéités composites en prenant en compte le caractère hétérogène du GSB. Cette étude numérique de prédiction des fuites a fourni de précieux renseignements sur la réduction des débits de fuite qui s'est avérée dépendante des courbes de rétention respectives du géotextile et de la bentonite. En effet, en plus de contenir la bentonite et de promouvoir la résistance aux efforts de traction, le géotextile contribue à la réduction des débits de fuite grâce sa faible perméabilité une fois désaturé. Des expérimentations de mesure du débits de fuites et de quantification de la transmissivité d'interface pour différentes combinaisons d'étanchéités composites GM-GSB ont permis de conclure que ni les sollicitations extérieures conduisant à l'altération chimique, physique et hydraulique des GSB, conduisant à augmenter la conductivité hydraulique du GSB de 4 à 5 ordres de grandeurs, ni la qualité de la GM (en lien avec la rugosité, rigidité, épaisseur) n'affectent significativement les transferts dans les étanchéités composites GM-GSB en régime permanent. Un effet est décelé sur le régime transitoire. Une approche synthétique est enfin donnée et permet une analyse globale des transferts dans les étanchéités composites GM-GSB afin de prévoir au mieux l'impact des transferts dans les barrières vers l'environnement. Il a été confirmé que le gonflement de la bentonite intervient dans la réduction du débit de fuite à travers les étanchéités composites indépendamment de la forme granulométrique de la bentonite. Par ailleurs, il a été conclu que la phase transitoire durant laquelle une réduction des débits de fuites à travers la barrière d'étanchéité composite est observée, doit être prise en compte dans les calculs de fuites pour une prédiction réaliste et rigoureuse des fuites à travers la barrière. / Geomembrane (GM)-geosynthetic clay liner placed in bottom landfill barriers could be faced to advective transfers caused by the appearance of GM defects. Leachate could percolate naturally through the GCL; penetrate the soil and the ground water which could result in environmental damage. It's therefore important to understand leakage transfer mechanism though GM-GCL composite liners and quantify them in order to be aware, to master and to minimize advective transfer impact through the barrier to the environment. However, the GM is not accessible in bottom land fill liner which makes it harder to correctly estimate leakage through the composite liner. The present work used to evaluate advective transfer through GM-GCL composite liner via an experimental and numerical approach. This work used also to ameliorate the comprehension of leakage transfer mechanisms as a function of external solicitations compared to the literature. The purpose of this study was thus to properly identify advective transfer problem through composite liner by highlighting the main parameters affecting advective transfers through composite liners (confining stress, heterogeneity of the GCL composition, contact quality at the GM-GCL interface, chemical and physical alteration during its service life). Those parameters influence the whole barrier performances and material characteristics evolution were studied. Experimental program was dealing with acquiring water retention data of GCL by introducing state surface concept under the confining stress generated by the waste. It has been found that confining stress used to reduce GCL swelling facilities while hydrated and consequently lead to the decrease of the saturated hydraulic conductivity. This result emphasizes on landfill conception recommendations based on rapidly covering the GCL in in order to acquiring its watertightness capabilities. This experimental study was reinforced by a numerical computation study dealing with water transfer through composite liner due to a GM defect and a hydraulic head with considering the heterogeneity of the GCL. This numerical study highlighted new phenomena regulating flow rate through composite liners consisting of geotextile deaturation due to high suction performed by the bentonite as part of the GCL. Indeed, in addition to containing the bentonite and providing tensile shear stresses, the geotextile contribute to reduce the flow rate through composite liner thanks to its high hydraulic conductivity while desaturated. Interface transmissivity tests were also carried for different composite liners combinations. Measurements concluded that neither external solicitations resulting from chemical and physical alteration conducting to the increase of the hydraulic conductivity of the GCL nor the quality of the GM (in relation with its roughness, rigidity and thickness) significantly affect advective transfer through composite liners at the steady state. However an effect was highlighted during transient state. A synthetic approach was thus given to summarize composite liners transfers mechanism and anticipate environmental impact of its leakage. It was later confirmed that the bentonite swelling contributes also to flow rate reduction through composite liner regardless of the bentonite nature and granulometry. In addition, it has been concluded that the transient state within which flow rate reduction though has been observed must be taken into consideration to realistically predict flow rate leakage through composite liners.

Ce travail de recherche porte sur le comportement des milieux poreux (triphasiques), plus particulièrement les sols non saturés sous sollicitations hydro-mécaniques. Un modèle constitutif élastoplastique couplé est développé. Ce modèle original est formulé selon les principes suivants: une loi constitutive est développée pour décrire le comportement de chaque phase (squelette solide, liquide, et gaz). Ensuite, des relations de couplage sont ajoutées entre chacune des phases. Pour le comportement du squelette solide, une loi élastoplastique non associée est adoptée, avec deux surfaces de charges, en cisaillement et en compression. La partie hydrique est décrite par une formulation qui permet de prendre en compte l'effet d'hystérésis. Ce modèle a été enrichi par une relation de couplage hydromécanique qui permet d'exprimer la pression d'entrée d'air en fonction de la porosité. Ensuite, le couplage complet se fait avec la contrainte effective de Bishop en utilisant une nouvelle définition du paramètre de succion χ grâce à laquelle, les différents phénomènes présents dans la réponse des milieux poreux sous différentes sollicitations peuvent être reproduits. Ce modèle est validé par une confrontation à des données expérimentales issues de la littérature sur différents types de sol (sable, limon,...). Le modèle est implanté dans le code aux éléments finis Cast3M. L'analyse de problèmes particuliers, tels que la mise en œuvre d'un cas test d'un sol d'assise soumis à un cycle pluvial, ainsi que l'étude de la stabilité d'une pente, permette de montrer la capacité du modèle à reproduire le comportement des milieux poreux non saturés.

Cette thèse rapporte les résultats d'études expérimentales effectuées à l'?domètre avec imposition de succion par la méthode osmotique sur des matériaux gonflants compactés lâches et naturels denses. Plusieurs cycles de séchage/humidification ont été appliqués sur ces matériaux sous trois faibles charges mécaniques constantes. Pendant ces cycles, les éprouvettes manifestent un retrait cumulé pour le sol lâche et un gonflement cumulé pour le sol dense. Les résultats montrent que dans les deux cas, les déformations volumiques convergent vers un état d'équilibre où le sol présente un comportement réversible. A la fin des cycles de succion, un cycle de chargement/déchargement a été effectué sous les succions constantes. Les valeurs de la pression de préconsolidation p0(s), de l'indice de compression vierge lambda(s) et de l'indice de compression élastique lambda dépendent directement des chemins de contrainte suivis.<br>L'ensemble des résultats expérimentaux permet de déterminer les surfaces de charge : la limite de séparation de micro/macro (Lm/M) ; la surface de chargement-effondrement (LC) et la surface de comportement saturé (SCS). La succion limite entre la micro- et la macrostructure (Lm/M) dépend parfaitement de la structure interne et du diamètre qui délimite les deux familles de pores. L'évolution de la pression de préconsolidation en fonction de la succion imposée est présentée par la surface LC. Les courbes de compressibilité sous différentes succions convergent vers la courbe correspondant à l'état saturé sous de fortes contraintes appliquées. La pression à partir de laquelle, le sol continue son chemin sur la courbe du comportement normalement consolidé est appelée la pression de saturation (Psat). Plus la succion imposée est élevée, plus la charge nécessaire pour atteindre cette pression de saturation est importante. La surface SCS présente la variation de la pression de saturation en fonction de la succion imposée. Nous pouvons considérer que les surfaces de charge SCS et LC sont uniques pour les sols denses cependant elles se superposent à la fin des cycles de succion pour les sols lâches. Les cycles hydriques augmentent aussi la limite (Lm/M) entre la micro- et la macrostructure pour les deux sols.

Dans les milieux fracturés, l'analyse des causes menant un versant fracturé stable à évoluer vers de grands<br />mouvements de terrain met en évidence le rôle de l'eau comme facteur aggravant et déclencheur des instabilités<br />de versants. Les interactions à long terme eau / roche sont interprétées à partir de la composition chimique de<br />l'eau des sources, qui renseigne sur l'état d'altération du versant (facteur aggravant) et sur le passage des ondes<br />de crue dans le massif (facteur déclencheur).<br />Ainsi, l'hydrochimie permet de reconstituer les écoulements dans les versants instables. Une approche à<br />plusieurs échelles des relations entre les écoulements et les structures gravitaires, souligne le rôle de drain de ces<br />structures. Les comparaisons entre les sites permettent de discuter l'évolution à long terme des effets de l'eau sur<br />la déformation des milieux fracturés. Les mouvements gravitaires lents (basculement, effondrement)<br />décompriment le versant, favorisent l'infiltration des eaux de surface et provoquent l'apparition d'aquifères<br />perchés dans les secteurs décomprimés. Avec la déstabilisation, l'effet saisonnier des variations de pressions<br />d'eau diminue, mais l'influence de la charge, liée au volume d'eau infiltrée, augmente fortement. Cette charge<br />s'applique en partie haute du versant et peut déclencher une rupture.<br />De plus les équilibres hydrogéochimiques observés dans les zones stables sont rompus lorsque le massif est<br />soumis à une déformation. La dynamique de dissolution des sulfates évolue avec les déformations du massif.<br />Ainsi, l'hydrochimie enregistre les déformations actuelles, mêmes faibles, des versants et constitue un outil<br />performant pour la surveillance des zones instables.

Les sols gonflants compactés sont largement utilisés dans de nombreuses applications, notamment en géotechnique et en géotechnique de l'environnement. Après leur mise en place, ces matériaux risquent d'être soumis à des sollicitations hydriques et/ou mécaniques complexes qui peuvent entraîner de fortes modifications de leurs propriétés hydromécaniques susceptibles de mettre en cause leur fonctionnement initial. Il est donc important de prévoir le comportement hydromécanique de ces matériaux compte tenu des applications sensibles auxquelles ils sont destinés.<br />L'équipe de Barcelone a développé un code de calcul aux éléments finis (Code_Bright) pour l'analyse des problèmes couplés (THM) dans des géomatériaux, avec l'intégration du modèle élastoplastique (BBM) pour les sols non saturés non gonflants. Ce modèle est capable de décrire correctement le comportement hydromécanique des sols non saturés mais ne tient pas compte de quelques aspects particuliers observés sur des sols gonflants. Un deuxième modèle (BExM) a été alors proposé pour prendre en considération ces aspects. <br />L'objectif principal du travail présenté ici est : (i) d'implanter le modèle élastoplastique BExM pour les sols gonflants non saturés dans le code de calcul aux éléments finis Code_Bright ; (ii) de valider le modèle numérique à partir d'essais réalisés au laboratoire sur des sols gonflants ; et (iii) d'appliquer ce modèle aux simulations des problèmes pratiques. Pour atteindre cet objectif, une nouvelle famille de procédures numériques adaptées au modèle BExM a été introduite dans le code. L'équation de la fonction de charge de ce modèle pour un état donné de contrainte déviatorique a été déterminée de manière à faciliter les calculs de ses dérivées. La validation numérique du modèle a été vérifiée par des simulations numériques d'essais œdométriques à succion contrôlée, réalisées sur trois sols gonflants différents. Les résultats des simulations numériques de ces essais ont montré les aptitudes du modèle numérique implanté à reproduire correctement les résultats expérimentaux. Enfin, le modèle a été appliqué à deux problèmes pratiques : le stockage de déchets radioactifs dans les couches géologiques profondes et l'effet du retrait-gonflement du sol sur une fondation superficielle. Les résultats obtenus ont montré la capacité du modèle implanté à modéliser des problèmes couplés hydromécaniques.