Les géomatériaux sont présents dans la nature sous des formes très diverses : sols et rochesin situ solides, argiles ductiles, boues quasiment liquides, etc... La géomécanique s’intéresseà la compréhension du comportement solide des géomatériaux. Cependant, il arrive que sousl’effet de conditions extérieures particulières, des terrains initialement solides se transformenten fluides : c’est ce qui se produit par exemple lors de coulées de boues ou de débris. Dans untel contexte, il existe peu d’outils numériques capables de modéliser les différentes phases ducomportement. Il semble de plus, qu’il n’existe pas à ce jour de modèle constitutif satisfaisantpour décrire une telle transition.Nos travaux s’intéressent de manière générale à la transition solide-fluide dans le comportementdes géomatériaux, et à l’élaboration d’un modèle constitutif décrivant la phase solide, laphase fluide, ainsi que la transition entre les deux. Nous avons choisi dans ce cadre de menerles calculs en nous basant sur la méthode MEFPIL (Méthode aux Éléments Finis avec desPoints d’Intégration Lagrangiens) qui a déjà montré de fortes potentialités pour décrire descomportements très variés (dont des comportements à variables d’histoire), dans un mêmemodèle.Après avoir implanté et validé la première loi élasto-plastique dans Ellipsis (code basé sur laMEFPIL), nous avons pu introduire dans ce code le modèle de transition. Celui-ci se base surl’évolution du comportement solide élasto-plastique vers un comportement fluide, visqueux àseuil, et ce, au moment de la rupture matérielle détectée par le critère du travail du secondordre.Après quelques applications du modèle de transition solide-fluide sur des cas simples et homogènes(en considérant la loi élasto-plastique Plasol et loi visqueuse de Bingham), nous avonsappliqué ce modèle à la modélisation des coulées de boue de Sarno et Quindici (Italie, 1998).Les premiers modèles montrent la possibilité de décrire les trois phases de ce mouvement deterrain (l’initiation, la propagation et l’immobilisation), et nous avons pu étudier l’effet dedifférents paramètres sur l’arrêt contre un ouvrage de protection. / Geomaterials are present in nature in many forms : solid soil or rock, soft clay, almost liquidmud, etc... Geomechanics deals with the understanding the solid behavior of geomaterials.However, solid ground can happen, under specific external conditions, to turn into fluid : asfor example during mudflows or debris flows. In such a context, there are few numerical toolsable of modeling the different phases of the behavior. Furthermore, it seems that there is,nowadays, no satisfactory constitutive model to describe such a transition.Our work concerns, in a general way, solid-fluid transition in geomaterials behavior and thedevelopment of a constitutive model describing both the solid phase, fluid phase, and thetransition between the two. In this framework, we chose to carry out calculations with theFEMLIP numerical method (Finite Element Method with Lagrangians Integration Points)which has shown a strong potential to describe a wide variety of behaviors (including historydependant behavior), in a unique model.Having implemented and validated the first elasto-plastic law in Ellipsis (FEMLIP basedcode), we have introduced in this code the solid-fluid transition model. This last is based onthe evolution, at the failure state detected by the second order work criterion, of the solidelasto-plastic behavior towards a viscous fluid behavior, exhibiting a yield stress.After validation of the solid-fluid transition model in homogeneous cases (considering Plasolelasto-plastic law and Bingham viscous one), we applied this model to the modeling of Sarnoand Quindici mudflows (Italy, 1998). The first models shows the possibility to describe thethree phases of the flow (initiation, propagation and immobilization), and we could study theeffect of various parameters on the stop against a protection work.

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.