Dans ces travaux, nous nous intéressons à l'estimation et l'étude des intervalles cardiaques. L'objectif principal de cette thèse est donc de développer de nouveaux outils d'estimation de ces intervalles. Plus précisément, nous souhaitons concevoir de nouvelles techniques de traitement du signal pour extraire des ECG les intervalles P-R, qui sont très peu étudiés dans la littérature de part la difficulté d'extraire les ondes P, et dans un second temps les intervalles Q-T. Inspirés de la méthode itérative de Woody (1967) basée sur une technique de corrélation de chaque observation avec une moyenne des signaux réalignés, nous proposons tout d'abord un perfectionnement de cette méthode. En effet, à partir du même modèle d'observations, nous développons le critère de maximum de vraisemblance et nous obtenons une amélioration de la méthode de Woody au sens de l'optimalité. Le modèle d'observations de cette méthode est ensuite enrichi par une onde potentiellement parasite additionnée au bruit, et une généralisation de la méthode de Woody est proposée. En prenant l'exemple de l'onde T de l'ECG qui se superpose à l'onde P du battement suivant lorsque la fréquence cardiaque est élevée, différents modèles de l'onde T sont proposés. Grâce à ces deux méthodes proposées, les objectifs d'estimer et d'analyser les intervalles cardiaques, tels que les intervalles R-R, Q-T et P-R, sont atteints. Plus précisément, grâce à la généralisation de Woody, l'estimation des intervalles P-R à l'exercice est réalisable. L'application de ces méthodes aux enregistrements réels de différents types d'ECG (repos, exercice, fibrillation auriculaire,...), nous a révélé de nouveaux résultats concernant l'évolution des intervalles cardiaques, comme par exemple : l'existence d'un phénomène d'hystérésis en sens horaire dans la relation P-R/R-R, et une caractérisation possible des sujets en fonction de leur niveau d'entraînement à l'aide de la pente des intervalles P-R au début de la phase de récupération. Une modélisation de la réponse du Q-T à un changement de la période cardiaque a également été proposée s'inspirant du comportement électrique cellulaire.

Les matériaux magnétiques doux (Fe, Ni, Co, ... sous forme cristalline, polycristalline ou amorphe) sont utilisés en électrotechnique pour convertir l'énergie, guider les lignes de flux et transmettre des signaux. Leur principale caractéristique est de s'aimanter aisément. Ces matériaux sont aussi conducteurs et des courants sont induits en régime transitoire et périodique. Ces derniers peuvent être classiquement<br />diffusés à une échelle macroscopique mais aussi et surtout localement induits dans une microstructure magnétique en mouvement et réarrangement incessant. Ces effets d'amortissement et de pertes dynamiques ont plusieurs conséquences parfois difficilement prévisibles:<br />des pertes d'énergie, des temps de retard, de la distorsion de signal et de la rémanence. Nous nous intéressons à la recherche de modèles<br />fidèles dans le but de comprendre les propriétés des matériaux ferromagnétiques doux en régime statique et dynamique et de simuler avec précision le comportement des dispositifs électrotechniques. Nous choisissons de nous concentrer particulièrement sur l'effet des courants induits microscopiques en plus de ceux macroscopiques, puisqu'ils sont à l'origine des pertes en excès et de l'hystérésis dynamique qui sont observées expérimentalement.<br />Après avoir longuement étudié la problématique, les enjeux et les modèles existants; nous avons entrepris la construction d'une représentation des matériaux à l'échelle mésoscopique, intermédiaire entre l'échelle microscopique et macroscopique, utilisable aussi bien en Ingénierie des Matériaux qu'en simulation numérique type Eléments Finis. Pour ce, les équations de champs dans la matière ont été redérivées en présence des processus physiques dynamiques mis en jeu.<br />Ensuite nous nous sommes attachés à résoudre analytiquement certains problèmes simples pour appréhender les apports et limites de<br />notre modélisation. Il s'est agit de comparer la théorie à l'expérience et de caractériser des échantillons sur bancs de mesure normalisés.<br />Enfin, nous avons mis en oeuvre, à l'aide d'une méthode numérique (la Méthode des Eléments Finis), des formulations électromagnétiques 3-D et 2-D dédiées. Nous les avons testées sur des cas tests simples en les comparant systématiquement aux résultats standards et antérieurs, aux calculs analytiques et aux observations. Des premières configurations simples d'applications de sécurité électrique ont aussi été simulées et étudiées : un transformateur de courant et un actionneur électromécanique.

L'archipel du Cap Vert constitue de 10 iles volcaniques appartient a la zone sahelienne qui s'etend de l'atlantique jusqu'a la mer rouge. Depuis, plusieurs decennies le Cap Vert est affecte par la desertification causee en grande partie par la recession climatique et l'erosion des sols. Ces facteurs, associes a la forte pression anthropique sur les ressources, a l'orographie accidentee et a des pluies tropicales parfois diluviennes, provoquent de serieuses pertes du patrimoine foncier. Cependant, depuis son Independance en 1975, le Gouvernement a mene un vaste programme d'arborisation, de restauration des terres et d'amenagement des cours d'eau. Pourtant, tres peu de recherches ont ete menees pour evaluer les actions de protection et de conservation des sols et des eaux. Par consequent, il n'existe quasiment pas de donnees sur la problematique de la degradation des terres ni de bilans. Dans le cadre de ce travail, nous avons etudie les differents facteurs qui controlent l'erosion hydrique des sols. Nous avons plus particulierement cherche a differencier les effets des activites humaines, notamment agricoles, de ceux des facteurs climatiques comme les precipitations et la generation des ecoulements. Nous avons egalement etabli les premiers bilans d'exportations de matieres en suspension et en solution dans le contexte de l'archipel du Cap Vert. L'etude a ete menee a l'echelle de trois bassins versants de l'ile de Santiago Cap-Vert. Ces trois bassins versant (Longueira, Grande et Godim) sont localises dans la partie centrale de l'ile de Santiago et representatifs des divers modes d'occupation du sol et des differents climats de l'ile. Il existe un gradient climatique entre les trois bassins versants. En effet, Longueira qui presente une superficie de 4,18 km2, une pente moyenne de 47 %, se localise dans une zone humide couverte a 69 % par une foret et une surface agricole de 17 %. Grande avec une superficie de 1,87 km2, se localise en zone sub humide pour une pente moyenne de 50 %, il est essentiellement agricole. Godim, avec une superficie de 2,0 km2, se localise en zone semi aride, il est particulierement agricole et sa pente moyenne est de 32 %. Pour ces trois bassins versants, les ecoulements de crue a l'exutoire ont ete mesures et echantillonnes de 2004 a 2009. Le bassin versant de Longueira a fait l'objet d'un suivi plus pousse, notamment en termes de frequence d'echantillonnage et de suivi des ecoulements hors crue. Sur chaque echantillon nous avons procede a la determination de la concentration des matieres en suspension ainsi qu'a l'analyse des elements majeurs. Les resultats obtenus montrent que l'erosion mecanique dans les 3 bassins versants est caracterisee par une forte variabilite spatiale et temporelle. Sur la periode 2005-2009, le bilan moyen annuel pour les bassins versants de Longueira, Grande et Godim est de : 4266, 157 et 10,1 t.km2.an-1 respectivement. La saison humide 2006 a ete la plus erosive pour l'ensemble des trois bassins versants et particulierement dans Longueira avec 2 crues exceptionnelles qui ont genere une concentration moyenne de matieres en suspension superieure a 100 g/l. En revanche, les saisons 2005 et 2008 ont ete dans l'ensemble peu erosives car les concentrations moyennes ne depasserent pas 20 g/l. Par ailleurs, il n'y a pas eu de lames d'eau ecoulees pour les saisons 2005 et 2007 pour le bassin de Godim. Sur le bassin de Longueira, l'etude des phenomenes d'hysteresis permet de caracteriser chaque crue et de montrer que l'evolution temporelle des exportations de matieres en suspension au cours de la saison est fortement influencee par les activites agricoles. En effet, la premiere crue provoque l'exportation massive des sediments disponibles et localises dans le lit du cours d'eau. En consequence, la seconde est moins exportatrice de sediments. Un mois apres les premieres pluies, les activites de sarclage diminuent la densite du couvert vegetal et destructurent la partie superficielle des sols, ce qui provoque a nouveau une tres forte exportation de sediments lors de la troisieme crue. Les resultats de l'erosion chimique sur le bassin de Longueira indiquent que le taux d'erosion chimique moyen s'eleve a 45 t.km2.an-1 avec une forte variabilite temporelle. En effet, les saisons les plus humides de 2006 et 2007 sont les plus exportatrices de matieres en solution, alors que 2005 a eu une faible exportation. L'utilisation du modele de melanges EMMA (End-Members Mixing Analysis) montre que les ecoulements hypodermique et profond, qui alimentent le cours d'eau en elements dissous, sont les principaux facteurs de l'erosion chimique. On montre ainsi que les ecoulements hors crue sont a l'origine de plus de 90% des flux d'erosion chimique. L'ecoulement superficiel, qui contribue a environ 70 % du debit du cours d'eau en crue, constitue un facteur de premier plan de l'erosion mecanique des sols.

L'objectif de cette étude est de développer un modèle de comportement adapté pour reproduire le comportement mécanique de l'élastomère avec l'évolution de la température. Pour réaliser cette étude, un travail expérimental préliminaire est effectué pour mettre en évidence les différents phénomènes mécaniques présents pour le matériau avec l'évolution de la température. Le matériau étudié est un élastomère synthétique fluoro-carboné (FKM), le Viton. Cet élastomère possède d'excellentes performances aux hautes températures et aux agressions chimiques. Ses propriétés sont utilisées pour diverses applications comme l'étanchéité, l'amortissement, l'isolation, etc. La caractérisation mécanique du matériau est réalisée à l'aide des essais classiques de traction et compression ainsi que des essais de relaxation sur une plage de température allant de −40◦C à 100◦C. Les essais de chargement cyclique mettent en évidence une présence non négligeable de l'hyperélasticité pour les élastomères. Les essais de relaxation nous permettent de décrire la viscosité du matériau. De plus, des essais de chargement interrompus par des relaxations ont aussi été effectués, pour différents taux de déformation avec des temps de relaxations de 900s. Ils font apparaître un équilibre hystérétique démontrant la présence de l'hystérésis pour notre matériau. Dans cette étude, nous avons fait le choix de ne considérer que la plage de température de −8◦C à 100◦C pour reproduire correctement le comportement du matériau. Le modèle utilisé est une loi de comportement phénoménologique basée sur une décomposition additive des puissances. Le modèle est constitué de trois contributions de contrainte : l'hyperélasticité, la viscosité et l'hystérésis. Ce modèle sera nommé Hyperélasto-Visco-Hystérésis (HVH) où la contribution hyperélastique est prise en compte à l'aide de potentiel phénoménologique de Hart-Smith. La viscosité est décrite à l'aide du modèle de Maxwell généralisé et la contribution hystérétique est prise en compte à l'aide d'un modèle original. Au cours de cette étude nous considérons simplement le premier cycle de chargement de notre matériau. Une méthode d'identification a été mise en place pour permettre de mieux prendre en compte les différentes contributions du modèle. Les identifications des paramètres matériau sont réalisées pour des isothermes comprises dans la plage de température étudiée. La validation du modèle sera dans un premier temps testé à température ambiante sur un cas de chargement industriel d'un joint en relaxation. Pour valider le modèle en fonction de la température nous simulons différents cas de chargement en isotherme. Ensuite, la simulation d'un essai de relaxation en compression subissant des cycles de températures est réalisée pour confronter le modèle à un calcul dépendant de la température.

Au cours de ce travail de thèse, nous avons effectué les premières mesures de l'aimantation d'une seule particule ferromagnétique à basse température (0.1 - 6 K) à l'aide de SQUID hystérétiques continus extrêmement performants : la sensibilité de nos mesures est de 10^4 µB. Nous avons étudié des systèmes multiples et variés, tout d'abord des particules faites par lithographie électronique de forme elliptiques et de dimensions supérieures à 50 nm, mais aussi des fils magnétiques et des agrégats. <br /> En ce qui concerne les particules élaborées par lithographie électronique, les plus petites ont montré que la variation du champ de retournement en fonction de l'angle est en accord qualitatif avec le modèle de rotation uniforme de Stoner Wohlfarth. Cependant le déclenchement du retournement de l'aimantation s'effectue par un processus de nucléation. Les études de la dynamique de retournement de l'aimantation faites sur ces particules et les fils magnétiques ont montré que le retournement de l'aimantation est activé thermiquement et peut être approximativement décrit par le modèle de Kurkijärvi pour les températures comprises entre 1 et 6 K. En dessous de 1 K, ce modèle est mis en défaut ce qui pourrait être expliqué par un effet quantique. En ce qui concerne la probabilité de retournement de l'aimantation nous avons trouvé pour les plus petites particules que cette probabilité est proche d'une exponentielle. Cependant il existe de nouveau une déviation à plus basse température, probablement de même origine que la déviation observée sur le modèle de Kurkijärvi. Enfin, la largeur de la distribution des champs de retournement augment à basse température. Ceci est sans doute dû à l'existence de défauts (rugosité de surface, défauts cristallins, impuretés, etc.) qui créent de multiples vallées et cols sur la surface d'énergie potentielle décrivant le système, ceux-ci étant en fait révélés lorsque l'énergie thermique est du même ordre de grandeur ou inférieure à ces fluctuations de la surface d'énergie.

No description available.

Les actionneurs piézoélectriques font partie des outils les plus utilisés dans les applications à l'échelle micro/nano-métrique (micromanipulation, microassemblage, micropositionnement, etc…). Du point de vue fonctionnel, on distingue les actionneurs mono-axe (permettant d'obtenir la déflection suivant une direction) et les actionneurs multi-axes (pouvant fléchir suivant plusieurs directions). La notoriété des actionneurs piézoélectriques est due à un certain nombre de performances telles qu'une large bande passante (plus du kHz possible), une très bonne résolution (de l'ordre du nanomètre), une faible consommation en énergie électrique, une grande densité de force, une facilité d'alimentation et d'intégration, etc. Cependant, ces actionneurs sont caractérisés par des non-linéarités fortes (hystérésis et la dérive lente), des oscillations mal-amorties, et sont sensibles à la variation des conditions ambiantes (en particulier à la variation de la température). Pour les actionneurs multi-axes, il s'ajoute un problème des couplages entre les différents axes de l'actionneur. Cette thèse propose des stratégies innovantes de commande des actionneurs piézoélectriques multi-axes pour contrer les problèmes sus-mentionnés. Ces stratégies sont groupées en deux catégories. La première catégorie concerne les techniques de commande en boucle fermée. Ces techniques sont les plus adaptées pour garantir la robustesse et un niveau de précision élevé pour les actionneurs piézoélectriques. Cependant, à l'échelle micro/nano-métrique, ces techniques sont limitées par un manque d'espace suffisant pour installer des capteurs de position. La deuxième catégorie concerne la commande en boucle ouverte dont l'intérêt majeur est lié au fait qu'il n'y a pas besoin de capteurs pour la commande, ce qui constitue un avantage en terme de coût et facilité d'intégration. Dans cette thèse, nous proposons d'abord les techniques de modélisation et de commande en boucle ouverte multivariables. Ensuite, nous faisons une analyse des effets de la température sur les actionneurs piézoélectriques et nous proposons des techniques de commande en boucle ouverte et en boucle fermée de ces effets. Enfin, une stratégie de commande en boucle fermée par découplage, visant à obtenir des correcteurs d'ordre réduit pour les actionneurs multi-axes est proposée. Toutes ces techniques sont vérifiées et appliquées expérimentalement à un actionneur piézoélectrique de type tube. / Piezoelectric actuators are among the most used tools in many applications at micro/nano-scale (micromanipulation,microassembly, micropositioning, etc). From a functional perspective, there exist mono-axis actuators(which are made to bend in one direction) and multi-axis actuators (which provide deflections in different directions).The popularity of piezoelectric actuators is especially due to their high resolution (nanometric resolution),the large bandwidth (greater than 1kHz possible), the low electrical power consumption, the high force density,the ease of integration in positioning systems, etc. However, piezoelectric actuators are characterized by hysteresisand creep nonlinearities, badly damped vibrations and they are sensitive to the variation of ambient conditions(especially to the temperature variation). In addition, multi-axis actuators exhibit cross-couplings betweentheir axis. This thesis proposes novel strategies for modeling and control of multi-axis piezoelectric actuators,with the aim to counteract the aforementionned problems. These strategies are grouped into two categories.The first category concerns feedback control techniques. These techniques are the most suitable to ensurethe robustness and a high level of precision for piezoelectric actuators. However, at the micro/nanoscale, thesetechniques are limited by the lack of enough physical space to install feedback sensors. The second categoryconcerns the feedforward control techniques. The main advantage of these techniques is related to the factthat, in feedforward control schemes, feedback sensors are not needed for tracking. This allows to achieve ahigh degree of packageability and the cost reduction. In this thesis, we first propose multivariable modelingand feedforward control techniques. Then, we analyse the effects of temperature variation on piezoelectricactuators and we propose feedforward and feedback control techniques for these effects. Finally, a feedbackstrategy based on decoupling techniques with an aim to reduce the order of feedback controllers for multi-axispiezoelectric actuators, is proposed. All these modeling and control strategies are experimentally applied on apiezoelectric tube actuator.

Cette thèse a pour objet l'étude des propriétés de transport électronique des nanotubes de carbone soumis à des pressions hydrostatiques de l'ordre du gigapascal. L'originalité de ce travail réside dans le fait d'étudier ces propriétés sur des nanotubes individuels. Ceci permet de simplifier la géométrie du système et de faire apparaître des comportements propres à chaque type de nanotubes. Le cas d'un nanotube multifeuillets composé d'un tube externe à faible bande interdite et d'un tube interne métallique a permis d'observer i) l'évolution sous pression de la barrière de Schottky aux contacts or-nanotube, ii) l'évolution de la résistance inter-feuillet, mettant en évidence une transition associée au changement de section du nanotube, iii) la diminution de l'hystérésis en tension de grille. D'autre part, une étude systématique sur des nanotubes métalliques permet de faire ressortir un comportement général pour le contact or-nanotube sous pression, indépendant de la chiralité du nanotube et du milieu transmetteur de pression. Nombre de ces effets peuvent être utilisés pour des applications électroniques ou électro-mécaniques, tels que des capteurs de pression miniatures et environ dix fois plus sensibles que certains standards actuels. Un modèle simple de calculs par la méthode des liaisons fortes est aussi mis en œuvre afin de prédire l'évolution des propriétés électroniques des nanotubes de carbone sous pression en fonction de leur chiralité. Ce modèle permet de prédire d'autres phénomènes qui pourraient être observés grâce à une étude approfondie et systématique utilisant la méthode expérimentale développée au cours de cette thèse / This thesis focuses on the electronic transport properties of carbon nanotubes under hydrostatic pressures as high as one gigapascal. The originality of this work is the study of these properties on individual nanotubes. This simplifies the geometry of the system and allows studying the behaviour of each type of nanotubes. The case of a multiwalled nanotube made of at least an external semiconducting tube and an internal metallic one led us to observe i) a pressure induced change of the Schottky barrier at the gold-nanotube contacts, ii) a evolution of the intershell resistance, featuring a transition associated to the change of the nanotube cross section, iii) the decrease of the gate voltage hysteresis. Additionally, a systematic study on metallic nanotubes allowed pointing out a general behaviour of the gold-nanotube contact under pressure, which is independent on the nanotube chirality and on the pressure transmitting medium. Many of these effects can be useful to design electronic or electro-mechanical devices, such as miniature pressure sensors that would be about ten times more sensitive than some of today's standards. A simple tight-binding model is also applied to predict the evolution of the carbon nanotube electronic properties under pressure with respect to their chirality. This model allows predicting other phenomena, which could be observed in the context of a deeper and systematic study using the experimental method that was developed in this thesis.

La problématique de la stabilisation des flammes non-prémélangées reste primordiale. Il faut pour la résoudre déterminer l’importance relative des phénomènes aérodynamiques, thermiques et chimiques intervenant dans les mécanismes de stabilisation. La démarche expérimentale utilisée pour cela au cours de cette thèse porte une attention particulière sur l’influence des conditions aux limites thermiques, et comporte deux volets à travers lesquels le rôle des transferts thermiques est mis en exergue. D’abord, la zone d’attachement d’une flamme stabilisée derrière la lèvre du brûleur est examinée, en considérant les couplages entre le positionnement du bout de flamme à proximité du brûleur et la température de la lèvre. Différentes régions ont ainsi pu être identifiées selon le comportement du bout de flamme, qui évolue depuis une nature diffusive vers une nature propagative à l’approche des limites aérodynamiques de stabilité. Par ailleurs, une modification des propriétés thermiques du brûleur a permis de mettre en évidence une évolution, avec la température de la lèvre, du rôle relatif des modes de coincement thermique et chimique de la flamme par la paroi du brûleur. Ensuite, l’étude concerne non plus un état stabilisé de la flamme, mais les transitions entre les différents régimes de combustion, et plus particulièrement le décrochage d’une flamme attachée. L’examen des conditions conduisant à la déstabilisation de la flamme est un moyen d’apporter des éléments-clefs de compréhension quant aux couplages et aux équilibres aérothermochimiques prévalant préalablement au décrochage. Une évolution du processus de décrochage a ainsi été mise en avant avec l’augmentation de la température initiale des réactants, en lien avec l’évolution de phénomènes transitoires d’extinction locale de la zone de réaction. / Non-premixed flame stabilization is still an important issue in combustion. Addressing this issue requires to evaluate the relative importance of aerodynamic, thermal and chemical phenomena involved in the stabilization mechanisms. This thesis develops to this end an experimental approach, with a particular focus on the influence of thermal boundary conditions, examining the role of heat transfer in a twofold analysis. At first, the attachment zone of a rim-stabilized jet-flame is investigated, by careful consideration of the coupling existing between the burner lip temperature and the flame attachment location relative to the burner. Several regions have been identified according to the flame leading edge behavior, which evolves from diffusive to propagative closer to the aerodynamic stability limits. Besides, by modifying the burner thermal properties, a change has been shown in the relative roles of thermal and chemical quenching of the flame by the burner wall, depending on the burner lip temperature. Secondly, the attention is directed to transitions between different combustion regimes, namely attached and lifted flames. In particular, beyond the stable state of an attached flame, its lifting process is investigated. Examining in which conditions destabilization of the flame occurs indeed appears to be an ideal opportunity to gain insight into the aerothermochemical coupling and equilibriums existing prior to lift-off. Thus, the lifting process has been shown to be modified by the reactant initial temperature, in close relationship with the change in occurrence of localized transitory extinction events of the reaction zone.

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. / This work focuses on the behaviour of porous triphasic media, particularly on unsaturated soils subjected to hydromechanical loading. A coupled elastoplastic constitutive model has been developed. This original model is formulated according to the following principles: (1) a constitutive law describing the behaviour of different phases (solid skeleton, liquid and gas). (2) coupling relationships between each phase. For the behaviour of the solid skeleton, a non associated elastoplastic constitutive law is adopted, with two loading surfaces: shear surface and compression cap surface. The hydric part is discribed using a formulation which allows to take into account the hysteresis effect. This model has been extended using a hydromechanical coupling relation between the air entry value and the porosity. Then the coupling is completed with the Bishop effective stress, using a new definition for the suction parameter χ. Using this formulation, the various phenomena present in the porous media behaviour under different loading can be reproduced. The developed model has been validated through a comparison with experimental data on different types of soil (sand, silt,…). This model is implemented in the french finite element code Cast3M. The analysis of specific problems, such as (1) the study of shallow foundation subjected to cyclic rain event, as well as (2) the study of slope stability, show the model capacity to reproduce the behaviour of unsaturated porous media.