Ce travail est consacré à la formulation et à la résolution des problèmes de mécanique des milieux poreux dans la situation, courante dans le domaine pétrolier, et susceptible de se produire en génie civil et en géotechnique, dans laquelle le squelette subit une transformation géométrique finie. Comme pour les milieux continus monophasiques (ou "secs"), la prise en compte des transformations finies nécessite d'utiliser un formalisme plus complexe que celui employé dans le cadre familier des petites perturbations. Sur le plan théorique, la principale difficulté réside dans la formulation du comportement, abordée ici en s'appuyant sur l'étude de la thermodynamique du milieu. On présente d'abord (chapitre 1) la modélisation mécanique des milieux poreux due à Biot, et le cadre thermodynamique mis en place par Coussy pour l'étude des milieux continus ouverts. On étudie dans le chapitre 2 la formulation du comportement poroélastique, en portant une attention particulière aux milieux dont le constituant solide est incompressible. Le chapitre 3 compare, pour l'étude de la consolidation et de la compaction unidimensionnelles en poroélasticité, les résultats obtenus en transformation finie avec les résultats d'une modélisation en transformation infinitésimale. Le chapitre 4 est consacré à la formulation du comportement poroélastoplastique en transformation finie. On propose notamment un modèle qui généralise le modèle Cam-Clay aux milieux poreux en transformation finie. Le chapitre 5 présente les principes des méthodes de résolution numérique en poroélastoplasticité finie, et la résolution d'un problème académique simple. En pratique, l'étude des problèmes abordés dans ce travail montre que la prise en compte complète des transformations finies permet d'éviter de commettre des erreurs significatives sans augmenter sensiblement la difficulté de la résolution des problèmes.

Après un bref rappel concernant l'écrouissage isotrope, on étudie l'écrouissage anisotrope : on développe un modèle basé sur l'hypothèse que la variable interne tensorielle d'écrouissage est un tenseur des déformations plastiques et on l'applique pour modéliser l'anisotropie élastique et plastique induite. A l'aide du modèle développé on étudie l'écrouissage cinématique. On développe aussi un autre modèle qui suppose que la variable interne d'écrouissage n'est pas le tenseur des déformations plastiques; la détermination de cette variable tensorielle se faisant à partir de sa loi d'évolution

Ce travail vise à accélérer les calculs lors de la simulation numérique du laminage à pas de pèlerin. Pour ce faire, nous nous sommes penchés sur la méthode Multimaillages Multiphysiques Parallèle (MMP) implémentée au sein du code Forge®, et destinée à accélérer les calculs pour des procédés de mise en forme à faible où la déformation est très localisée sur une petite zone du domaine. Dans cette méthode, un maillage localement raffiné dans la zone de déformation et plus grossier sur le reste du domaine est utilisé pour résoudre les équations mécaniques, alors qu'un maillage uniformément raffiné est retenu pour le calcul thermique. Le calcul mécanique étant le plus coûteux, la réduction du nombre de nœuds de son maillage permet d'obtenir des accélérations très importantes. Le maillage de calcul thermique étant uniformément fin, il sert aussi de maillage de stockage pour les champs calculés pour les deux physiques. Pour appliquer la efficacement méthode MMP au laminage à pas de pèlerin, plusieurs aspects importants ont été pris en compte. Tout d'abord la géométrie complexe du tube nécessite le développement d'une technique de déraffinement spéciale afin d'assurer un déraffinement maximal tout en garantissant un maillage convenable pour des calculs. Une technique de déraffinement de maillage utilisant une métrique anisotrope cylindrique a été alors développée. Ensuite, avec la loi de comportement élastoplastique utilisée, des perturbations importantes sont observées sur les contraintes dues aux diffusions numériques engendrées par les différents types de transports des champs P0 (constants du maillage thermique vers le maillage mécanique. Pour y remédier, une approche combinant deux techniques a été développée. La première consiste à effectuer la réactualisation des variables d'état directement sur le maillage mécanique plutôt que sur le maillage thermique et de les transporter ensuite. La deuxième technique est l'utilisation d'un opérateur de transport P0 basé sur un recouvrement super convergent (SPR) et la construction de champs d'ordre supérieur recouvrés. De bonnes accélérations sont alors obtenues sur les cas de laminage à pas de pèlerin étudiés, allant jusqu'à un facteur 6,5 pour la résolution du problème thermomécanique. Les accélérations globales de simulation vont jusqu'à un facteur 3,3 sur un maillage contenant environ 70 000 nœuds en séquentiel. En parallèle les performances chutent légèrement, mais elles restent semblables (2,7). / This work aims at speeding-up the calculations of the numerical simulation of the cold pilgering process. To this end, it is focused on a Parallel Multiphysics Multimesh (MMP) method that bas been implemented in the Forge® code; this method is dedicated to speeding-up the calculations for processes in which the deformation is localized within a small area of the computational domain. A locally refined mesh is used to solve the mechanical equations while a uniformly refined mesh serves as basic mesh to store state variables and is preferred for thermal calculations. The mechanical computations being the most expensive, reducing the number of nodes of its mesh provides high speed-ups. To effectively apply MMP method to cold pilgering process, many important aspects have been taken into account. Firstly the complex geometry of the tube requires the development of a special mesh coarsening technique, in order to ensure a maximum coarsening while guaranteeing a suitable mesh for calculations. A technique using a cylindrical anisotropic metric is then introduced. Afterwards, with the elastoplastic behaviour law used for the considered process, inaccuracies were observed on the stress field. They are mainly due to the numerical diffusion generated by the different transfers operations of P0 variables (constant per element) from the thermal mesh to the mechanical one. To remedy this issue, an approach combining two techniques has been developed. Firstly state variables are directly updated on the mechanical mesh, instead of doing it on the thermal mesh before transferring them on the the mechanical mesh. The second technique consists in using a P0 transfer operator based on super convergent recovery (SPR) technique, which improves the accuracy of the transported field through introduction of higher order recovered fields. High speed-ups are obtained on the studied cold pilgering cases, up to a factor 6,5 for the resolution of the thermomechanical problem, and the global simulation speed-up is up to a factor of 3,2, on a mesh with about 70 000 nodes in sequential calculations. For parallel calculations performances slightly drop but remain quite good.

Des modèles élastoplastiques écrouissables incorporant des lois de dilatance réalistes ont récemment été mis au point pour décrire le comportement des sables. La capacité portante en frottement des fondations sous charge axiale dépend largement de la valeur des contraintes normales au fût après mise en place et chargement ; ces contraintes peuvent, a priori, varier fortement avec la dilatance. Les études théoriques dans ce domaine se sont limitées jusqu'à présent aux cas de milieux élastiques et de milieux dilatants modélisés d'une manière relativement simple : milieu linéaire dilatant ou milieu élastique-parfaitement plastique. Ces comportements ne correspondent qu'imparfaitement au comportement réel des sols dilatants. Dans le travail présenté ici, on a étudié les modèles élastoplastiques écrouissables de Nova et al. à une seule surface de charge et de Vermeer à deux surfaces de charge. Ces modèles sont capables de représenter les principales caractéristiques des matériaux granulaires dilatants, au moins quand ils sont soumis aux trajets de chargement monotones des essais classiques de laboratoire. Ces deux modèles ont été programmés dans le code d'éléments finis ROSALIE et leurs réponses dans le cas d'un essai triaxial de compression ont été examinées. Ils ont ensuite été appliqués aux problèmes des pieux sous charge axiale. Les résultats des calculs effectués ont mis en évidence l'importance de la prise en compte du glissement avec frottement à l'interface sol-pieu pour tous les modèles de sols utilisés. Le comportement du sol, qui est identique pour les modèles élastique et élastique-parfaitement plastique, est différent pour les modèles élastoplastiques écrouissables, les déplacements axiaux étant plus élevés dans ces derniers cas. Une certaine augmentation de volume élastique a, par ailleurs, été relevée au voisinage du pieu pour tous les calculs. La dilatance plastique existe dans tous les calculs élastoplastiques mais elle se manifeste généralement assez tard, à peu près en même temps que le glissement relatif sol-pieu. L'effet conjugué de la dilatance plastique et des augmentations de volumes élastiques induit en moyenne une certaine diminution de la contrainte radiale sur le fût par rapport au calcul élastique. Au vu de ces résultats, et pour les modèles utilisés, il semble donc que la dilatance des sols granulaires n'accroisse pas la capacité portante en frottement, tout du moins par augmentation de la contrainte radiale.

L'objectif de ce travail est d'obtenir un outil numérique permettant de simuler et d'optimiser le procédé de mise en forme par compression et frittage de poudres métalliques. La démarche proposée repose sur une étude expérimentale du comportement des matériaux, une modélisation phénoménologique et l'utilisation du modèle dans un code éléments finis. Une étape de validation des calculs à partir de pièces fabriquées en milieu industriel complète la démarche. La première partie de ce manuscrit est consacrée à l'étude de poudres métalliques (mélanges fer-graphite et carbure de tungstène-cobalt). Une vaste étude expérimentale a plus particulièrement été entreprise sur le mélange fer-graphite. Les conditions expérimentales à adopter pour obtenir le comportement intrinsèque du matériau ont été étudiées. Une loi élastoplastique isotrope à une variable d'écrouissage a été ensuite retenue pour simuler la compression d'une pièce de référence complexe. La confrontation entre résultats numériques et expérimentaux montre quelques écarts qui sont discutés. Les déformations au cours du frittage sont étudiées sur les mélanges de poudres de carbure de tungstène et de cobalt à partir d'essais de dilatométrie. Les résultats expérimentaux ont permis d'identifier les paramètres d'une loi phénoménologique simple. La seconde partie de ce travail est consacrée à l'étude de pâtes d'hexaferrite de strontium destinées à la fabrication d'aimants permanents. Le comportement du matériau a été étudié sous diverses sollicitations mécaniques et thermiques. Nous avons montré que la simulation de la compression d'un segment ferrite à partir d'une approche monophasique était mal adaptée parce qu'elle n'était pas suffisamment représentative des phénomènes observés. Enfin les retraits anisotropes au cours du frittage ont été calculés à partir d'un modèle thermoélastique. L'agrément obtenu entre les résultats numériques et expérimentaux valide les calculs proposés.

Lors du forage des puits d'hydrocarbure, une pâte de ciment est coulée dans l'espace annulaire entre le cuvelage en acier et les formations géologiques traversées. Pompée à l'état liquide, cette pâte de ciment fait sa prise le long du puits sous différentes conditions de température et de pression. La gaine de ciment ainsi mise en place a pour principales fonctions de promouvoir l'étanchéité pour protéger le casing contre la corrosion, de fournir le support mécanique pour assurer la stabilité du puits et d'isoler les différents fluides dans les couches traversées. Au cours de sa vie dans le puits, depuis le forage à la complétion et de la production à l'abandon, la gaine de ciment est soumise à différentes sollicitations mécaniques et thermiques qui peuvent l'endommager et altérer ses principales fonctions. La réponse de la pâte de ciment soumis à ces sollicitations dépend non seulement des conditions d'hydratation mais aussi de l'histoire des chargements précédemment appliqués. La prédiction du comportement de la gaine de ciment doit donc se faire à l'aide d'une modélisation numérique qui nécessite une loi de comportement pour la pâte de ciment. Le but de cette thèse est d'établir une loi de comportement de la pâte de ciment en cours d'hydratation pendant le jeune-âge (les 144 premières heures). Pour ce faire, des essais calorimétriques, de mesures de vitesse des ondes et des essais œdométriques ont été réalisés sur une pâte de ciment pétrolier classe G (w/c = 0,44) en cours de prise. Les conditions d'hydratation explorées vont de 7 à 30°C pour les températures et de 0,3 à 45MPa pour les pressions. Les résultats expérimentaux ont montré que les déformations volumiques de la pâte de ciment dues à son hydratation (retrait macroscopique) sont considérablement influencées par la contrainte sous laquelle la pâte de ciment s'hydrate. Plus la contrainte d'hydratation est élevée, plus élevé est le retrait macroscopique à 144 heures. Inversement, les déformations irréversibles dues à un cycle de chargement mécanique à cet âge sont moins importantes pour les contraintes plus élevées. Les résultats ont également montrés qu'au cours de la prise du ciment, il existe un temps critique à partir duquel l'application des cycles de chargement mécanique crée des déformations résiduelles dans la pâte de ciment. Ce temps critique arrive à un degré d'hydratation relativement constant, compris entre 0,18 et 0,20. Le modèle « Boundary Nucleation and Growth » a été utilisé pour étudier la dépendance de ce temps critique à la pression et à la température. Pour la modélisation du retrait macroscopique et de la réponse contrainte – déformation de la pâte de ciment, un modèle élasto-plastique chemo-poro-mécanique couplé, prenant en compte la désaturation du milieu, a été développé. Ce modèle utilise une surface de charge fermée de type Cam-Clay et une loi plastique associée. La loi d'écrouissage dépend des déformations volumiques plastiques et du degré d'hydratation. Les paramètres du modèle ont été évalués pour simuler le retrait macroscopique de la pâte de ciment hydratée sous différentes contraintes et températures. A un degré d'hydratation donnée, le modèle permet également de simuler la réponse contrainte-déformation due à un chargement mécanique / When drilling oil & gas well, cement slurry is pumped between the casing and the rock formation. This cement slurry sets at different conditions of temperature and pressure. The role of this cement sheath is to provide zonal isolation of different fluid along the well, to protect the casing against corrosion and to provide mechanical support. During the life of the well, from drilling to completion, production and P&A (plug and abandonment), the cement sheath is submitted to various mechanical and thermal loading that can potentially damage its properties and alter its performance. The behavior of cement paste submitted to theses solicitations depends both on the hydration condition and the loadings previously applied on the cement paste. The prediction of cement sheath behavior should be done by numerical modeling, which needs a constitutive law for cement paste. The purpose of the present work is to establish a constitutive law of cement paste during its hydration at early age (first 144 hours). The approach is based on combined calorimetric, wave velocities and oedometric tests on an oil-well class G cement paste with water-to-cement ratio equals 0.44. The hydration conditions explored are 7 to 30°C for temperature and 0.3 to 45MPa for pressure. The experimental results showed that the volumetric strain due to cement hydration (macroscopic shrinkage) depends considerably on the hydration pressure. At 144 hours of hydration, the macroscopic shrinkage increases with the hydration pressure increase. But, the residual strain due to application of mechanical cycle at this age is less for cement hydrated under higher pressure. The experimental results revealed that during the hydration there is a critical time after which, the application of mechanical loading can potentially induce residual strain in cement paste. This time is reached at constant hydration degree between 0.18 and 0.20. The Boundary Nucleation and Growth model was used to model the pressure and temperature dependence of this critical time. A coupled elasto-plastic chemo-poro-mechanical model is developed to simulate the macroscopic shrinkage of cement paste hydrated at different conditions of temperature and pressure. A modified Cam-Clay type yield surface with associate flow rule is used. The hardening law depends both on the degree of hydration and on the plastic volumetric strain. At constant degree of hydration, the developed model permits to simulate the stress – strain behavior of cement paste due to the mechanical loading

La thèse porte sur la réponse d'une plaque couplée à un liquide et soumise au champ de pression créé par une détonation aérienne. Celle-ci correspond au mode d'explosion le plus violent et se caractérise par la propagation d'une onde de choc qui génère sur la plaque un champ de pression mobile dont l'intensité et la vitesse de propagation sont élevées. L'objectif est d'étudier la réponse dynamique du système couplé pendant la durée de propagation de l'onde sur la plaque. Après avoir présenté le chargement de détonation, la mise en équation du problème est exposée. Celle-ci est adaptée au contexte de dynamique rapide imposé par la sollicitation extérieure. La flexion de la plaque est étudiée selon la théorie de Mindlin Reissner et prend en compte les non linéarités géométriques. Les non linéarités matérielles sont décrites par la loi de comportement élastoplastique de Prandtl-Reuss avec écrouissage isotrope. La dynamique du fluide est décrite par les équations d'Euler linéarisées. L'étude analytique permet d'obtenir des solutions au problème d'une bande infinie, reposant sur un domaine liquide non borné, soumise à un chargement mobile stationnaire. Les solutions stationnaires décrivent la réponse du système couplé au voisinage du front de chargement. La résolution numérique du problème réel d'une plaque couplée est obtenue à partir d'un schéma aux différences finies d'ordre 2 en temps et en espace. L'intégration temporelle des équations est obtenue par un schéma explicite. L'étude expérimentale présente le banc d'essai et des exemples de réponse de plaques en contact avec de l'eau, soumises à des détonations. Ces réponses sont comparées aux solutions numériques.

La technique de renforcement des sols par inclusions linéaires rigides, c'est-à-dire capables de reprendre des efforts de flexion et cisaillement, connaît depuis quelques décennies un essor important bien que la simulation et le calcul de tels ouvrages reste un problème difficile, en raison de la très forte hétérogénéité du sol renforcé combinée au grand nombre d'inclusions mis en jeu dans de tels ouvrages. Une modélisation multi phasique a récemment été développée permettant d'aborder le dimensionnement de tels ouvrages de géotechnique avec des temps de calcul considérablement réduits par rapport à ceux d'une modélisation numérique classique, où les inclusions de renforcement doivent être discrétisées individuellement. On s'intéresse dans le cadre de ce travail à la mise en oeuvre de ce modèle dans le cadre du comportement élastoplastique des différents constituants, lorsque les effets de flexion et de cisaillement des inclusions sont pris en compte. Après avoir développé un certain nombre de solutions analytiques de référence, un outil numérique basé sur la méthode des éléments finis est mis au point, puis appliqué à la simulation du comportement de différentes structures renforcées, telles que les radiers de fondation sur groupes de pieux, ou la stabilisation de pentes instables par l'incorporation de telles inclusions.

Lors du forage des puits d'hydrocarbure, une pâte de ciment est coulée dans l'espace annulaire entre le cuvelage en acier et les formations géologiques traversées. Pompée à l'état liquide, cette pâte de ciment fait sa prise le long du puits sous différentes conditions de température et de pression. La gaine de ciment ainsi mise en place a pour principales fonctions de promouvoir l'étanchéité pour protéger le casing contre la corrosion, de fournir le support mécanique pour assurer la stabilité du puits et d'isoler les différents fluides dans les couches traversées. Au cours de sa vie dans le puits, depuis le forage à la complétion et de la production à l'abandon, la gaine de ciment est soumise à différentes sollicitations mécaniques et thermiques qui peuvent l'endommager et altérer ses principales fonctions. La réponse de la pâte de ciment soumis à ces sollicitations dépend non seulement des conditions d'hydratation mais aussi de l'histoire des chargements précédemment appliqués. La prédiction du comportement de la gaine de ciment doit donc se faire à l'aide d'une modélisation numérique qui nécessite une loi de comportement pour la pâte de ciment. Le but de cette thèse est d'établir une loi de comportement de la pâte de ciment en cours d'hydratation pendant le jeune-âge (les 144 premières heures). Pour ce faire, des essais calorimétriques, de mesures de vitesse des ondes et des essais œdométriques ont été réalisés sur une pâte de ciment pétrolier classe G (w/c = 0,44) en cours de prise. Les conditions d'hydratation explorées vont de 7 à 30°C pour les températures et de 0,3 à 45MPa pour les pressions. Les résultats expérimentaux ont montré que les déformations volumiques de la pâte de ciment dues à son hydratation (retrait macroscopique) sont considérablement influencées par la contrainte sous laquelle la pâte de ciment s'hydrate. Plus la contrainte d'hydratation est élevée, plus élevé est le retrait macroscopique à 144 heures. Inversement, les déformations irréversibles dues à un cycle de chargement mécanique à cet âge sont moins importantes pour les contraintes plus élevées. Les résultats ont également montrés qu'au cours de la prise du ciment, il existe un temps critique à partir duquel l'application des cycles de chargement mécanique crée des déformations résiduelles dans la pâte de ciment. Ce temps critique arrive à un degré d'hydratation relativement constant, compris entre 0,18 et 0,20. Le modèle « Boundary Nucleation and Growth » a été utilisé pour étudier la dépendance de ce temps critique à la pression et à la température. Pour la modélisation du retrait macroscopique et de la réponse contrainte – déformation de la pâte de ciment, un modèle élasto-plastique chemo-poro-mécanique couplé, prenant en compte la désaturation du milieu, a été développé. Ce modèle utilise une surface de charge fermée de type Cam-Clay et une loi plastique associée. La loi d'écrouissage dépend des déformations volumiques plastiques et du degré d'hydratation. Les paramètres du modèle ont été évalués pour simuler le retrait macroscopique de la pâte de ciment hydratée sous différentes contraintes et températures. A un degré d'hydratation donnée, le modèle permet également de simuler la réponse contrainte-déformation due à un chargement mécanique / When drilling oil & gas well, cement slurry is pumped between the casing and the rock formation. This cement slurry sets at different conditions of temperature and pressure. The role of this cement sheath is to provide zonal isolation of different fluid along the well, to protect the casing against corrosion and to provide mechanical support. During the life of the well, from drilling to completion, production and P&A (plug and abandonment), the cement sheath is submitted to various mechanical and thermal loading that can potentially damage its properties and alter its performance. The behavior of cement paste submitted to theses solicitations depends both on the hydration condition and the loadings previously applied on the cement paste. The prediction of cement sheath behavior should be done by numerical modeling, which needs a constitutive law for cement paste. The purpose of the present work is to establish a constitutive law of cement paste during its hydration at early age (first 144 hours). The approach is based on combined calorimetric, wave velocities and oedometric tests on an oil-well class G cement paste with water-to-cement ratio equals 0.44. The hydration conditions explored are 7 to 30°C for temperature and 0.3 to 45MPa for pressure. The experimental results showed that the volumetric strain due to cement hydration (macroscopic shrinkage) depends considerably on the hydration pressure. At 144 hours of hydration, the macroscopic shrinkage increases with the hydration pressure increase. But, the residual strain due to application of mechanical cycle at this age is less for cement hydrated under higher pressure. The experimental results revealed that during the hydration there is a critical time after which, the application of mechanical loading can potentially induce residual strain in cement paste. This time is reached at constant hydration degree between 0.18 and 0.20. The Boundary Nucleation and Growth model was used to model the pressure and temperature dependence of this critical time. A coupled elasto-plastic chemo-poro-mechanical model is developed to simulate the macroscopic shrinkage of cement paste hydrated at different conditions of temperature and pressure. A modified Cam-Clay type yield surface with associate flow rule is used. The hardening law depends both on the degree of hydration and on the plastic volumetric strain. At constant degree of hydration, the developed model permits to simulate the stress – strain behavior of cement paste due to the mechanical loading

Les glissements de terrain sont les risques naturels fréquents dans le monde. D'une part, ils concernent un comportement alterné pour les géomatériaux : le comportement solide décrit classiquement par la géomécanique, et celui fluide après la rupture. D'autre part, ils peuvent être induits par plusieurs facteurs. Par exemple, la géologie, la topographie, la végétation, la variationclimatique, et les activities humaines. Dans les régions tropiques et tempérées, il est bien connu que les glissements de terrain surviennent fréquemment après les pluies intensives, à cause de la saturation des sols non saturés.Comme jusqu'à présent, il n'y a aucun modèle constitutif satisfaisant pour décrire une telle transition pour les géomatériaux non saturés, nos travaux se concentrent d'abord à un modèle unifié, qui associe un modèle hydro-élasto-plastique pour les géomatériaux non saturés, un modèle visqueux de Bingham, et un critère de transition entre les phases solide et fluide. Le modèle unifié permet dedécrire complètement le comportement avant et après la rupture, induite par des conditions hydromécaniques, pour les géomatériaux non saturés dans un cadre unifié.Ensuite, la Méthode des Eléments Finis avec les Points d'Intégration Lagrangiens (MEFPIL), qui bénéficie des avantages des approches Lagrangiennes (suivre les variables internes) et de celles Eulériennes (traiter les grandes transformations), est améliorée et choisie pour effectuer les calculs dans ce cadre. Nous avons introduit une nouvelle formulation hydro-élasto-plastique avec le couplage hydromécanique basée sur la MEFPIL, et l'implanté dans l'Ellipsis (le code basé sur la MEFPIL). En outre, pour inverser la matrice non symétrique de la formulation, un nouveau solver basé sur la méthode stabilisée du gradient bi-conjugué (BiCGSTAB) a été aussi implanté.Enfin, plusieurs benchmarks sont proposés pour valider le comportement mécanique des géomatériaux non saturés décrit par le modèle, et un glissement de terrain heuristique induit par la pluie est simulé, nous pouvons étudier les effets de paramètres divers sur l'infiltration d'eau et sur la rupture. En outre, une simulation du glissement de terrain à Shien (Chine, 2012), avec des paramètres réels est présentée. / Landslides are common natural disasters all over the world. On one hand, they implicate an alternate behaviour for geomaterials: the solid behaviour dealt with classically by geomechanics, and the fluid behaviour after failure. On the other hand, they can be induced by several factors. For example,geology, topography, vegetation, climatic variation, and human activities. In tropic and tempered regions, it is well known that the landslides occur frequently after intense rainfalls, because of saturation effect of unsaturated soils.As there is, up to now, no satisfactory constitutive model to describe such a transition for unsaturated geomaterials, our work is thus focused firstly on a unified model, that associates an hydro-elastoplastic model for unsaturated geomaterials, a Bingham's viscous law, and a transition criterion between solid and fluid states. The model allows to describe both solid and fluid behaviours for unsaturated geomaterials in a unified framework, and is made possible to simulate completely the behaviour of unsaturated geomaterials before and after failure, induced by hydromecanical conditions.Secondly, the Finite Element Method with Lagrangian Integration points (FEMLIP), that benefits from both the Lagrangian approaches (track internal variables) and Eulerian approaches (solve large transformations), is ameliorated and chosen to carry out calculations in this framework. We have introduced a new hydro-elasto-plastic FEMLIP formulation with hydromecanical coupling, and implemented it in Ellipsis (FEMLIP based code). Moreover, in order to inverse nonsymmetric matrixin the formulation, a new solver based on the biconjugate gradient stabilized method (BiCGSTAB) has been also implemented.Finally, several benchmarks are proposed to validate the model for the main features of unsaturated geomaterials, and a heuristic rainfall-induced landslide is simulated, we could study the effects of various parameters on the water infiltration and on the failure. Moreover, a simulation of the Shien landslide in China (2012) with real parameters is presented.