La simulation numérique du comportement et l'analyse de stabilité des ouvrages de géotechnique renforcés par inclusions linéaires "rigides" ou "souples", c'est-à-dire capables ou non de reprendre des efforts de flexion et cisaillement, 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 employées, ainsi qu'à leurs dimensions caractéristiques très petites vis à vis de celles de l'ouvrage dans son ensemble. Une approche alternative, fondée sur une modélisation dite "multiphasique", a été récemment développée, qui permet de s'affranchir des difficultés précédentes et d'aboutir ainsi à des méthodes de dimensionnement accessibles à l'ingénieur. On s'intéresse dans le cadre de ce travail aux développements théorique et numérique de cette approche en prenant en compte non seulement les effets de flexion et de cisaillement des inclusions, mais également les différents types d'interaction entre ces dernières et le sol environnant. La première partie du mémoire est consacrée à la construction des approches du calcul à la rupture et à la mise en œuvre numérique élasto-plastique du modèle multiphasique des ouvrages renforcés par inclusions "souples", qui est alors appliqué à l'analyse de la stabilité d'ouvrage en "terre armée". L'extension de l'approche du calcul à la rupture et du code éléments finis en élastoplasticité, développés dans la première partie, à la prise en compte des effets de flexion et de cisaillement fait l'objet de la seconde partie du mémoire. La version la plus complète du modèle multiphasique élaborée est appliquée à la simulation du comportement d'une fondation renforcée par inclusions "rigides", soumise à un séisme induisant un chargement latéral. Partant de l'idée d'optimiser le schéma de renforcement de l'ouvrage traité à la deuxième partie, en inclinant les inclusions, un modèle "triphasique" est formulé et développé. Ce modèle vise à décrire le comportement macroscopique d'un sol renforcé par deux familles d'inclusions ayant des orientations différentes.

L'objet de ce travail est l'analyse expérimentale et théorique du comportement mécanique anisotrope d'un grès des Vosges, dans les domaines des déformations élasto-plastiques, de la localisation des déformations et de la rupture. Un intérêt particulier porte sur l'influence de l'orthotropie de révolution initiale du matériau sur l'orientation des bandes de cisaillement apparaissant pour des états de compressions sous pressions de confinement. Un programme expérimental a été développé pour la caractérisation mécanique du comportement de la roche sous compressions isotropes, compressions orientées sous pressions de confinement et tractions simples orientées. Des méthodes expérimentales spécifiques pour les mesures de toutes les composantes du tenseur des déformations élasto-plastiques d'une roche poreuse sous pressions de confinement ont été conçues et développées dans le souci permanent d'obtenir des données expérimentales fiables. Les essais réalisés ont permis, d'une part, d'obtenir un ensemble de résultats réalistes et cohérents sur les caractéristiques élasto-plastiques orthotropes de révolution de la roche, et d'autre part, de mettre en évidence l'orientation privilégiée des bandes de cisaillement par rapport aux structures internes orientées du matériau. Une méthode de détection du seuil de localisation des déformations a également été développée, Parallèlement, l'orientation privilégiée des bandes de cisaillement apparaissant dans des roches orthotropes de révolution soumises à des compressions triaxiales de révolution a été modélisée. Cette modélisation a été réalisée par l'introduction des structures matérielles dans l'analyse de Rice de la localisation des déformations élasto-plastiques, grâce aux formes générales invariantes des lois de comportement anisotrope.

Ces travaux de thèse portent sur l'étude expérimentale des propriétés physiques et de la biomécanique des agrégats cellulaires embryonnaires. Le but de cette thèse était d'une part de mieux comprendre l'origine biologique de la viscosité et de la tension de surface tissulaire, d'autre part d'étudier quantitativement en détail l'élasticité cellulaire par des nouvelles mesures de rhéologie en cisaillement. Un premier chapitre concerne les mesures de tension de surface tissulaire par la méthode de compression et de viscosité tissulaire par analyse de la cinétique de fusion de deux agrégats en faisant varier comme paramètre principal la contractilité cellulaire que certains suspectent comme étant la principale origine biologique de ces paramètres. Nous utilisons le formalisme du DITH (Haris, 1976: Differential Interfacial Tension Hypothesis) pour interpréter les données. Le deuxième chapitre concerne les mesures rhéologiques en cisaillement à l'aide d'un rhéomètre commercial plan-plan sur plusieurs centaines ou milliers d'agrégats cisaillés ensembles. Nous montrons que les cellules deviennent moins rigides pour une déformation minimale d'environ 4%, mais sur l'échelle de l'heure les cellules sont capables de se rigidifier à nouveau. Ces expériences sont analysées à l'aide d'un modèle de ressorts qui cassent sous contrainte puis se ressoudent à contrainte nulle.

Les collecteurs d'échappement des véhicules automobiles sont des pièces sollicitées en fatigue thermique de manière prépondérante. Pour la fabrication de ces pièces, les aciers inoxydables viennent aujourd'hui concurrencer la fonte traditionnellement utilisée.<br />Afin de qualifier les différentes nuances de sa gamme vis-à-vis de cette application, la société Ugine & ALZ, fabricant d'aciers inoxydables, a développé un essai original de fatigue thermique. Si cet essai apparaît représentatif des sollicitations thermomécaniques subies par les collecteurs en zones critiques, différents mécanismes de ruine sont observés entre les nuances ferritiques et austénitiques.<br />Le sujet de la présente thèse concerne l'endommagement en fatigue thermique d'aciers inoxydables pour l'application automobile collecteur d'échappement. Deux aciers inoxydables sont étudiés: une nuance ferritique nommée F17TNb (17%Cr stabilisé Ti et Nb) et une nuance austénitique nommée R20-12 (20%Cr et 12%Ni).<br />Le premier objectif est de comprendre la différence d'endommagement apparaissant en fatigue thermique entre les nuances ferritiques et austénitiques. Le second objectif est de développer un outil de dimensionnement des structures sollicitées en fatigue thermique.

L'objectif de cette thèse est d'étudier expérimentalement les phénomènes de vieillissement dus à la diffusion des atomes en solution dans les alliages d'aluminium et les instabilités qui leur sont associées comme le phénomène Portevin-Le Châtelier ou les bandes de Lüders et de proposer une modélisation de ces phénomènes dans le cadre de la thermodynamique des processus irréversibles. Une étude expérimentale détaillée est alors entreprise sur les alliages d'aluminium AA5083-H116 et AA5182-O. Le comportement du premier présente l'effet PLC de façon prononcée et les deux types d'instabilités sont observées simultanément pour le second. La corrélation d'images numériques et la thermographie infrarouge sont essentiellement employées pour détecter et caractériser les aspects spatiotemporels des instabilités observées. La déformation non homogène due à l'apparition et la propagation de bandes de localisation est mise en évidence. Ces bandes de déformation sont visualisées, permettant à leurs diverses caractéristiques (vitesse, orientation, largeur, vitesse de déformation à l'intérieur des bandes, augmentation de la température à l'intérieur des bandes) d'être mesurées sur des éprouvettes lisses (plates, cylindriques ou prismatiques). Dans le cas des éprouvettes plates, l'effet des épaisseurs de l'éprouvette a aussi été examiné. Les caractéristiques des bandes sont aussi analysées sur d'autres géométries d'éprouvette (entaillées avec divers types d'entailles) à des vitesses de déformation différentes pour exhiber leur morphologie en présence de chargements multiaxiaux. Des simulations non linéaires et tridimensionnelles ont été effectuées en utilisant le modèle de McCormick pour montrer comment la prise en compte des phénomènes de vieillissement, même partielle, permet de décrire les hétérogénéités et le mouvement des bandes de déformation ainsi que de prévoir leurs différentes caractéristiques. Enfin, en se basant sur les mécanismes physiques à la base des phénomènes de vieillissement et en soulignant les limites du modèle de McCormick, un modèle élasto-viscoplastique prenant en considération les phénomènes de vieillissement est proposé dans le cadre de la thermodynamique des processus irréversibles.

Un grand nombre de systèmes thermomécaniques industriels se trouve confronté à des régimes transitoires plus ou moins rapides suivant la fréquence de fonctionnement. L'amélioration de leurs performances nécessite l'utilisation de nouvelles structures du type multimatériaux ou barrière thermique. En effet, ces matériaux peuvent être de type multicouche en associant plusieurs couches rangées de façon à améliorer le comportement mécanique et thermique d'un système ou alors constitués d'un substrat revêtu d'une succession de couches minces obtenues par projection thermique par exemple.Dans un système donné, ces matériaux subissent généralement des sollicitations cycliques qui peuvent être d'origine thermique et/ou mécanique. Il est donc nécessaire de mieux connaître leur comportement thermomécanique en régimes élastique et plastique. Ainsi, l'étude présentée dans ce travail, limitée ici à des conditions périodiques uniquement d'origine thermique, traite de l'évolution de l'endommagement d'un matériau sous une ou plusieurs formes de fatigue thermique.L'origine de la sollicitation imposée provient d'une condition de flux périodique (sous forme d'échelon, de triangle ou de sinus) prenant en compte les pertes par convection. Sur le plan mécanique, le matériau est supposé fixe sur l'une de ses deux extrémités et libre de se déformer sur l'autre. Les contraintes et les déformations mécaniques dans le matériau proviennent essentiellement des différences des coefficients de dilatation thermique et des gradients de température dans le matériau. La nature variable et transitoire du comportement thermique du matériau permet de suivre l'évolution de la distribution des contraintes et des déformations au sein du matériau.L'étude de son endommagement est menée selon les cas, soit sur des modèles établis directement à partir du comportement thermo élastique soit sur des modèles nécessitant l'étude thermo-élastoplastique. Dans les deux cas, comme la plupart des modèles d'endommagement (Lemaître et Chaboche) rencontrés dans la littérature ne sont valides que sur des matériaux uniformes et homogènes, une recherche de matériau équivalent du multi-matériau étudié était nécessaire. L'équivalence entre le matériau réel et le matériau équivalent repose sur un critère d'équivalence thermique. Les modèles étudiés fournissent dans les deux cas, l'évolution de l'endommagement du matériau, en fonction des paramètres géométriques et aussi de la forme des sollicitations thermiques imposées telles que le coefficient d'échange par convection, l'amplitude et la période du flux imposé.Une application de ces modèles sur un exemple de moteur à combustion interne est proposée à la fin de ce mémoire. Elle montre une prédiction du nombre de cycles (durée de vie) du cylindre moteur en fonction des conditions de fonctionnement utilisées.

Nous présentons deux approches différentes pour décrire le couplage hydromécanique des géomatériaux. Dans une approche de type phénoménologique nous traitons le milieu poreux comme un milieu continu équivalent dont les interactions entre la phase fluide et le squelette solide constituent le couplage du mélange à l'échelle macroscopique. En caractérisant le comportement de chaque phase nous arrivons à décrire le comportement couplé du milieu couplé saturé.Nous utilisons cette approche pour modéliser des essais expérimentaux faits sur un cylindre creux pour une roche argileuse (argile de Boom). Les résultats expérimentaux montrent de façon claire que le comportement de cette roche est fortement anisotrope. Nous avons choisi de modéliser ces essais en utilisant une lois de comportement élasto-plastique pour laquelle la partie élastique est transversalement isotrope.Le problème aux conditions aux limites étudié met en évidence des déformations localisées autour du forage intérieur. Afin de décrire de façon objective le développement de ces bandes de cisaillement nous avons considéré un milieu continu local de type second gradient qui permet d'introduire une longueur interne. De ce fait nous avons pu étudier le problème d'unicité en montrant qu'un changement de la discrétisation temporelle du problème aux limites peut conduire à des solutions différentes.Dans la deuxième approche étudiée nous caractérisons la microstructure du matériau avec des grains et un réseau de canaux pour la phase fluide. À l'aide d'un processus numérique d'homogénéisation nous arrivons à calculer numériquement la contrainte du mélange et le flux massique. Cette méthode d'homogénéisation numérique a été implémentée dans un code aux éléments finis afin d'obtenir des résultats macro. Une validation de l'implentation est proposée pour des calculs en mecanique pure et en hydromécanique.

Afin d'optimiser la mise aux milles au laminage à froid, la thèse consiste à développer un modèle prédictif de variation de largeur à la fois précis et rapide pour des utilisations en temps réel. Des efforts ont commencé en 1960s en développant des formules empiriques. Par la suite, la Méthode des Bornes Supérieures (MBS) est devenue la plus connue. [Oh 1975] utilisant le champ de vitesse 3D "simple" prédit bien la variation de largeur au laminage en conditions d'un tandem finisseur. [Komori 2002] a proposé une combinaison des champs fondamentaux et obtenu une structure informatique peu dépendante aux champs de vitesse. Néanmoins, seuls deux champs fondamentaux ont été introduits qui forment un sous-ensemble de la famille 3D "simple". [Serek 2008] a étudié des champs de vitesse quadratique qui inclue la famille "simple" et donne des meilleurs résultats avec un temps de calcul plus long. Le premier résultat de la thèse est un modèle 2D (MBS) avec des champs de vitesse oscillante. Ce modèle aboutit à une vitesse optimale qui oscille spatialement le long de l'emprise. Les résultats (puissance, vitesse...) sont plus proches des ceux de Lam3-Tec3 que la MBS 2D "simple". Pour une modélisation 3D, nous avons choisi la MBS avec la vitesse 3D "simple" et obtenu un très bon accord avec les expériences réalisées sur des produits étroits à Arcelor Mittal [64]. En outre, un nouveau modèle MBS est développé pour une bande bombée et des cylindres droits. Les résultats montrent que la variation de largeur diminue avec la bombée de la bande et correspondent bien à ceux de Lam3-Tec3. Cependant, la MBS admet un comportement rigide-plastique tandis qu'au laminage des bandes larges les déformations élastique et thermique ont des impacts importants sur la déformation plastique. Les modèles existant prenant en compte ces phénomènes [23,64] sont couteux en temps de calcul. Ainsi, l'idée est de décomposer la variation de la largeur de plastique en trois termes : les variations de largeur totales, élastique et thermique à travers la zone de plastique déterminés par trois nouveaux modèles simplifiés. Les deux premiers permettent d'estimer les variations de largeur élastique et plastique avant et après l'emprise. Ils donnent aussi les conditions aux limites au modèle d'emprise qui est en effet la MBS avec le champ de vitesse 3D "simple" permettant d'estimer la variation de la largeur totale. En outre, avec les puissances de déformation et de dissipation plastique de frottement données par le même modèle, la variation de largeur thermique est également obtenue. Le modèle de variation de largeur est donc appelée UBM-Slab combiné, très rapide (0,05 s) et prédit avec précision la largeur de variation par rapport à Lam3-Tec3 (<6%) / In order to save material yields in cold rolling process, the thesis aims at developing a predictive width variation model accurate and fast enough to be used online. Many efforts began in the 1960s in developing empirical formula. Afterward, the Upper Bound Method (UBM ) became more common. [Oh 1975]'s model with 3D "simple" velocity field estimates well the width variation for finishing mill rolling conditions. [Komori 2002] proposed a combination of fundamental ones to obtain a computer program depending minimally on the assumed velocity fields. However, only two fundamental fields were introduced and formed a subset of the "simple" family. [Serek 2008] studied a quadratic velocity family that includes the "simple" one and leads to better results with a higher computing time. Focusing on UBM , the first result of the thesis is a 2D model with an oscillating velocity field family. The model results to an optimum velocity that oscillates spatially throughout the roll-bite. The optimum power and the velocity field are closer to Lam3-Tec3 results than the "simple" one. For 3D modelling, we chose the 3D "simple" UBM and carried a comparison to the experiments performed at Arcelor Mittal using narrow strips [64]. A very good agreement is obtained. Further, a new UBM model is developed for a crowned strip with cylindrical work-rolls. It shows that the width variation decreases as a function of the strip crown and the results match well those of Lam3-Tec3 . However, the UBM considers only a rigid-plastic behaviour while in large strip rolling, the elastic and thermal deformations have important impacts on the plastic one. There exist some models considering these phenomena [23,64] but they are all time-consuming. Thus, the idea is to decompose the plastic width variation into three terms : total, elastic and thermal width variations through the plastic zone that are determined by three new models. The simplified roll-bite entry & exit models allow estimating the elastic and plastic width variations before and after the roll-bite. They give equally the longitudinal stresses defining the boundary conditions for the roll-bite model which is indeed the 3D "simple" UBM approximating the total width variation term. Moreover, with the plastic deformation and friction dissipation powers given by the same model, the thermal width variation term is also obtained. The width variation model, called UBM-Slab combined is very fast (0.05s) and predicts accurately the width variation in comparison with Lam3-Tec3 (<6%)

Les liquides confinés présentent beaucoup de comportements fascinants, très différents de ceux qui sont observés dans leur volume. Le confinement peut induire un déplacement de l'équilibre des phases (par exemple de la transition liquide-vapeur, aussi appelé condensation capillaire), il peut modifier la température de transition vitreuses des polymères, ou bien imposer un ordre dans l'arrangement moléculaire du fluide. Les modifications des propriétés mécaniques des liquides aux interfaces sont particulièrement importantes au niveau des applications. Cependant au niveau de la compréhension, le simple cas des liquides newtoniens est toujours sujet à controverse, avec d'une part des simulations numériques montrant que la viscosité ne devrait pas être modifiée pour des confinements supérieurs à quelques tailles moléculaires, et d'autre part des expériences non unanimes, montrant parfois des modifications qualitatives des propriétés rhéologiques sous confinement. Récemment nous avons montré que les méthodes d'impédance hydrodynamique en géométrie sphère-plan constituent une méthode privilégiée, non-intrusive et non-ambigüe, pour aborder la nano-mécanique des liquides aux interfaces (1,2). S'agissant d'interphases, cad de couches fluides dont les propriétés sont modifiées par la proximité d'un solide, il est possible d'accéder à leur module sans contact, donc sans la perturbation apportée par une seconde surface. S'agissant de l'effet du confinement sur la rhéologie, nous avons montré que la déformation élastique à l'échelle du pico-mètre des surfaces confinantes, donne une forte modification de la rhéologie apparente du fluide, même en l'absence de tout effet intrinsèque. Le sujet de thèse vise à mettre en oeuvre les méthodes d'impédance hydrodynamique pour étudier la rhéologie de solutions de polymères confinés. On étudiera plus précisément deux systèmes modèles d'importance fondamentale aussi bien que pratique : les brosses de polymères greffés, dont les propriétés mécaniques sont un enjeu dans les applications de lubrification aussi bien que pour les écoulements biologiques, et les solutions de polymères hydro-solubles d'intérêt pour la récupération assistée du pétrole, en vue de comprendre les effets de fluidification sous confinement et de faire la part entre modification de la viscosité et couche de déplétion induite par l'écoulement. Au niveau instrumental, un des enjeux de la thèse sera de mettre en oeuvre les mesures d'impédance hydrodynamique sur deux types d'instruments complémentaires au niveau de l'échelle de la sonde: l'appareil de mesure de forces dynamique (SFA) du Liphy, et l'AFM à détection interférométrique développé à l'Institut Néel. Ces différentes échelles d'investigation devront permettre de préciser les propriétés moyennes mécaniques moyennes des liquides confinés et leurs gradients au voisinage de la paroi. Une perspective du travail sera de mettre en regard les propriétés mécaniques et rhéologiques intrinsèques des brosses polymères déterminées directement sur SFA ou AFM, avec leur propriétés fonctionnelles: propriétés de lubrification des contacts frottants, ou de modification des écoulements des dans micro-canaux. Ceci sera poursuivi sur la plateforme expérimentale mise en place par Lionel Bureau au Liphy : SFA de friction, systèmes micro-fluidiques à visée biomimétique (parois fonctionnalisées par des brosses polymères). L'enjeu sera alors de comprendre comment les propriétés mécaniques et rhéologiques des brosses déterminent celles des systèmes dans lesquels elles interviennent. / Liquids confined present many fascinating behaviors very different from those observed in their volume. Confinement can induce a shift in the balance of phases (eg the liquid-vapor transition, also called capillary condensation), it can change the glass transition temperature of the polymer, or impose order on the molecular arrangement of fluid. The changes in the mechanical properties of liquid interfaces are particularly important in applications. However the level of understanding, the simple case of Newtonian liquids is still controversial, with one hand, numerical simulations show that the viscosity should not be changed for some higher molecular sizes containment, and secondly non-unanimous experiences, sometimes showing qualitative changes in rheological properties under confinement. Recently we have shown that the methods of hydrodynamic impedance sphere-plane geometry is a privileged, non-intrusive method and unambiguous, to discuss the mechanics of nano-liquid interfaces (1,2). As interphase, ie fluid layers whose properties are modified by the proximity of a solid, it is possible to accede their contactless module, so without the disturbance caused by a second surface.S As regards the effect of confinement on the rheology, we have shown that the elastic deformation across the pico meter of confining surfaces, gives a strong modification beyond apparent rheology of the fluid, even in the absence of any intrinsic effect. The thesis aims to implement the hydrodynamic impedance methods to study the rheology of polymer solutions confined. We specifically consider two models of fundamental importance as well as practical systems: brushes grafted polymer whose mechanical properties are an issue in lubrication applications as well as for biological flows and solutions of water-soluble polymers interest in enhanced oil recovery, in order to understand the effects of thinning containment and to distinguish between changes in viscosity and depletion layer induced by the flow. At the instrumental level, one of the challenges of the thesis is to implement the hydrodynamic impedance measurements on two complementary instruments at the level of the probe: the measuring dynamic power (SFA) of Liphy, and AFM interferometric detection developed at the Institut Néel. These different scales of investigation will help to clarify the medium average mechanical properties of liquids confined and their gradients near the wall. A view of work will be to look mechanical and rheological properties of polymer brushes intrinsic determined directly on SFA or AFM with their functional properties: lubricating properties of sliding contacts, or modification of the flow in microchannels. This will continue on the implementation by Lionel Bureau Liphy experimental platform: SFA friction advised biomimetic micro-fluidic systems (walls functionalized polymer brushes). The challenge will be to understand how the mechanical and rheological properties of brushes determine those systems in which they operate.

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.