La propagation des ondes sismiques dans les milieux poreux multiphasiques présente des enjeux nombreux, tant sur le plan environnemental (risques naturels, géotechnique, pollutions de nappes...) que pour les réservoirs (aquifères, hydrocarbures, stockages de CO2...). L'utilisation des ondes sismiques pour étudier ces milieux se justifie par le fait qu'en se propageant, les ondes sont déformées par le milieu qu'elles traversent et contiennent ainsi des informations aux capteurs sur les phases fluides et solides et sur le squelette poreux. Ce travail de thèse s'intéresse aux caractéristiques des ondes sismiques dans les milieux multiphasiques (plusieurs phases fluides et solides), depuis la description physique jusqu'à la caractérisation des paramètres constitutifs par inversion, en passant par la modélisation numérique 2D de la propagation. La première partie du travail a consisté à décrire la physique des milieux multiphasiques (phase par phase et leurs intéractions dynamiques) en utilisant des méthodes d'homogénéisation pour se ramener à un milieu équivalent défini par sept paramètres. Ainsi, dans des milieux simple porosité saturés et dans des milieux plus complexes (double porosité, partiellement saturés ou visco-poroélastiques), je peux calculer la propagation des ondes sismiques sans approximation. En effet, j'utilise une méthode numérique dans le domaine fréquence-espace qui permet de prendre en compte tous les termes qui dépendent de la fréquence sans approximation. La discrétisation spatiale utilise une méthode d'éléments finis discontinus (Galerkin discontinu) qui permet de considérer des milieux hétérogènes.Je montre notamment que les attributs sismiques (vitesses et atténuations) des milieux poreux complexes sont fortement dispersifs et les formes d'ondes complètes, calculées sans approximation, sont fortement dépendantes de la description physique du milieu. La caractérisation des paramètres poroélastiques s'effectue par inversion. Une méthode en deux étapes a été proposée : la première consiste en une inversion ''classique'' (tomographie, inversion des formes d'ondes complètes) des données (sismogrammes) pour obtenir des paramètres macro-échelles (attributs sismiques). La seconde étape permet de reconstruire, à partir des paramètres macro-échelles, les paramètres poroélastiques micro-échelles. Cette étape d'inversion utilise une méthode d'optimisation semi-globale (algorithme de voisinage). Une analyse de sensibilité montre qu'en connaissant a-priori certains paramètres, on peut inverser avec précision les paramètres du squelette poroélastique ou retrouver la nature du fluide saturant, à partir des vitesses de propagation. En revanche, pour retrouver la saturation en fluide, il est préférable de connaître les atténuations. Deux applications réalistes (monitoring de réservoir et hydrogéophysique) mettent en oeuvre ce type d'inversion en deux étapes et démontrent qu'à partir de données estimées par des méthodes classiques d'imagerie, on peut remonter à certains paramètres poroélastiques constitutifs.

La présente étude a permis de mettre au point une méthode itérative de calcul non-linéaire des structures super-cavitantes. Le coefficient de pression au bord d'attaque de la structure est obtenu avec une grande précision, ce qui est essentiel pour la prévision du domaine de cavitation. En écoulement 2D le cas de l'aile super-cavitante à extrados mouillé est traité et la méthode permet de déterminer le point de décollement de la cavité à l'extrados. Le calcul de l'aile à extrados complètement dénoyé est rendu ainsi possible. Le cas des cavités très courtes est obtenu avec une grande précision. Ces trois problèmes n'avaient jamais été résolus auparavant. En écoulement 3D, pour des ailes à envergure finie, une méthode de résolution matricielle spécifique a été mise au point. C'est une méthode de Gauss Seidel par blocs, chaque bloc étant inversé par la méthode de Householder (orthogonalisation de Schmidt). La méthode nécessite des temps de calcul importants mais la richesse des indications qu'elle fournit pourra la faire préférer à d'autres méthodes moins onéreuses sur le plan informatique mais ne permettant pas d'avoir une description de l'écoulement aussi précise. Un champ d'application particulièrement intéressant a été celui des structures sub-cavitantes fonctionnant avec une cavité partiellement développée à partir du bord d'attaque. Les résultats sont très précis et comparés à certaines valeurs expérimentales.

Une méthode de calcul non-linéaire d'écoulement autour de structures minces par la méthode des singularités est présentée. Le concept de minimisation est développé et appliqué à des configurations sub-cavitantes, partiellement cavitantes et super-cavitantes. Ce concept est relatif aux erreurs numériques commises sur les débits de fuite. Le calcul non-linéaire de la géométrie de la poche de cavitation est obtenu par une procédure itérative. La méthode est validée par différentes comparaisons des résultats avec les résultats de référence issus de la littérature. Le cas de la grille d'aubes a été traité, les résultats comparés à ceux obtenus avec une méthode d'éléments finis. La comparaison est excellente. Le code de calcul a été étendu au calcul 3D de l'écoulement autour d'une aile partiellement cavitante et super-cavitante en présence d'une surface libre ainsi qu'au calcul 3D d'une hélice sub-cavitante, suivie de son sillage, en écoulement tournant. Les résultats sont confrontés aux résultats expérimentaux.

L'objectif de cette recherche est de construire un modèle réduit de petite dimension pour prévoir les réponses dynamiques dans une bande BF sur les parties rigides d'un véhicule automobile complet. Un tel modèle réduit "léger" est une aide à la phase de conception en "Avant Projet" de ces véhicules qui ont la particularité de présenter de nombreux modes élastiques locaux en BF dues à la présence de nombreuses parties flexibles et d'équipements. Pour la construction du modèle réduit, nous avons introduit une base non usuelle de l'espace admissible des déplacements globaux. La construction de cette base requiert la décomposition en sous-domaines du domaine de la structure qui peut présenter une très grande complexité géométrique et dont les modèles EF font intervenir de très nombreux types d'éléments finis. Cette décomposition en sous-domaines a été réalisée par la Fast Marching Method que nous avons due étendre pour pouvoir traiter la complexité des modèles EF des véhicules automobiles. Puis les équations matricielles du modèle EF sont projetées sur cette base. Afin de prendre en compte les incertitudes sur les paramètres du modèle, les incertitudes de modèle induites par les erreurs de modélisation et enfin les incertitudes liées à la non prise en compte des contributions locales dans le modèle réduit des déplacements globaux, un unique modèle probabiliste non paramétrique de ces trois sources d'incertitude a été implémenté sur le modèle réduit construit avec les vecteurs propres globaux. Les paramètres de dispersion de ce modèle probabiliste ont été identifiés en utilisant le principe du maximum de vraisemblance et des réponses obtenues à l'aide d'un modèle stochastique de référence qui inclut des informations expérimentales résultant de travaux précédents. Le modèle réduit stochastique, pour la prévision des déplacements globaux sur les parties rigides dans la bande BF qui a été développé, a été validé sur un modèle de structure automobile "nue" puis a été appliqué avec succès sur un modèle complet de véhicule automobile

Ce travail porte sur l'étude expérimentale et théorique du comportement des gouttelettes d'eau (instabilité, déformation, coalescence), déposées à la surface d'un isolateur composite haute tension, en présence d'un champ électrique. La déformation des gouttelettes d'eau est étudiée en fonction de plusieurs paramètres tels que le volume, le nombre et la conductivité des gouttelettes ainsi que leurs positions par rapport aux électrodes. Un intérêt particulier est porté à la charge accumulée sur la surface d'une gouttelette. L'influence des gouttelettes sur la répartition du champ électrique ainsi que le facteur d'amplification du champ électrique par une méthode numérique sont étudiés. Une formulation mathématique de la déformation d'une gouttelette d'eau posée sur une surface d'isolateur est développée. Les mécanismes ainsi que les paramètres impliqués dans la coalescence d'une paire de gouttelettes d'eau sous l'effet du champ électrique sont également traités. Un modèle mathématique permettant d'analyser la coalescence d'une paire de gouttelettes d'eau posées sur la surface d'un isolateur en présence du champ électrique est proposé. Ce modèle est basé sur une approche à interface diffuse qui consiste essentiellement en un couplage entre les équations de Navier-Stokes, le calcul de la force électrique et une équation de type Cahn-Hilliard décrivant l'évolution de l'interface en prenant en compte les phénomènes de tensions superficielles. La démarche suivie et les résultats présentés dans ce travail ouvrent de multiples perspectives tant d'un point de vue expérimental que d'un point de vue de la modélisation et de la simulation numérique des phénomènes physiques intervenant sur des isolateurs de lignes de transport et de distribution de l'énergie électrique dans des conditions d'humidité (pluie, rosée ...).

Afin de surmonter les difficultés de la méthode des éléments finis (MEF) liées à la forte dépendance de la solution au maillage, des méthodes sans maillage ont été développées durant les 20 dernières années. Ces techniques ont l’avantage de procurer des approximations très régulières, capables de répondre de manière plus satisfaisante aux exigences croissantes des applications. Cependant, certaines caractéristiques intrinsèques à la plupart de ces approches rendent leur mise en œuvre difficile : souvent des techniques supplémentaires spécifiques doivent être mises en place pour imposer les conditions aux limites et traiter les discontinuités physiques. Récemment, la méthode des éléments naturels (MEN) est apparue, se basant sur les concepts de diagramme de Voronoï et de voisins naturels. C’est une approche capable d’associer les avantages d’une approximation très régulière propre aux méthodes sans maillage et une mise en œuvre quasiment identique à la MEF. Cette thèse porte sur l’étude et le développement de la MEN dédiée aux applications du génie électrique. Le but principal de ce travail exploratoire est l’étude des limitations ainsi que des avantages et du potentiel de la MEN et ses concepts sous-jacents. Les analyses de performances de la méthode sont présentées. Sur les points ouverts tels que l’intégration numérique, la montée en ordre et l’interpolation vectorielle, des développements originaux sont proposés. / In order to overcome the limitations related to the finite element method’s (FEM) narrow dependency of the solution on the mesh, meshless or meshfree methods were developed over the last 20 years. These techniques present the advantage of yielding very smooth approximations, being able to respond more adequately to the increasing demands of applications. However, some intrinsic features of most of these approaches make the implementation difficult, often requiring additional specific techniques for the imposition of the boundary conditions and the treatment of physical discontinuities. Recently, the natural element method (NEM) was developed. This approach, based on the Voronoi diagram and the “natural neighbors” concepts, combines the advantages of very smooth approximations and a FEM-like implementation. This thesis focuses on the study and development of the NEM, dedicated to electrical engineering applications. The main purpose of this exploratory work is the study of the limitations, benefits and the potential of the NEM and its underlying concepts. Several analyses of NEM’s performance are presented. As far as the numerical integration, higher order approximations and the vector interpolation are concerned, original developments are proposed.

En l’espace d’une quarantaine d’année, les connaissances sur l’angiogenèse tumorale ont littéralement explosé. Dans les années 1970, Judah Folkman démontre que les tumeurs ont besoin d’être vascularisées pour continuer à proliférer. Peu de temps après, les protagonistes principaux de l’angiogenèse tumorale ont été découverts, de même que les mécanismes dans lesquels ils sont impliqués. La décennie suivante marque le début des recherches sur les molécules à visée anti-angiogénique et c’est en 2004 que le bevacizumab (Avastin, Roche), premier médicament anti-angiogénique utilisé en oncologie, voit le jour. Parallèlement à cela, l’essor de la modélisation à effets-mixtes couplée aux progrès des outils informatiques ont permis de développer des méthodes d’analyses de données de plus en plus performantes. Dès 2009 L’agence de régulation FDA (Food and Drug Administration) aux Etats-Unis a identifié le rôle central de la modélisation numérique pour mieux analyser les données d’efficacité et de toxicité, préclinique et clinique en cancérologie. Le but de ce projet est d’étudier les effets de différents inhibiteurs de l’angiogenèse sur la dynamique tumorale, en se basant sur une approche populationnelle. Les modèles développés seront des modèles à base d’équations différentielles ordinaires intègrant des données et des informations issues de la littérature. L’objectif de ces modèles est de caractériser la dynamique des tailles tumorales chez les animaux et chez les patients, afin de comprendre les effets des traitements anti-angiogéniques et apporter un soutien pour le développement de ces molécules ou pour la prise de décisions thérapeutiques par les cliniciens. / Within the last 40 years, knowledge of tumor angiogenesis has literally exploded. In the seventies, Judah Folkman demonstrated that tumors need to be vascularized to continue to proliferate. Shortly after, the main protagonists of tumor angiogenesis have been discovered, as well as the mechanisms in which they are involved. The next decade is the beginning of the research on molecules with anti-angiogenic effects and in 2004 bevacizumab (Avastin, Roche), the first antiangiogenic drug used in oncology, was available for treating solid cancer patients. Along with this, the increasing interest of mixed-effects modeling coupled with advances in computer tools allowed developing more efficient methods of data analysis. In 2009, the regulatory agency FDA (Food and Drug Administration) in the United States has identified the central role of numerical modeling to better analyze the efficacy and toxicity preclinical and clinical oncology data. The aim of this project is to study the effects of different angiogenesis inhibitors on tumor dynamics, based on a population approach. The developed models are models based on ordinary differential equations and that integrate data and information from the literature. The objective of these models is to characterize the dynamics of tumor sizes in animals and patients in order to understand the effects of anti-angiogenic treatments and provide support for the development of these molecules, or to help clinicians for therapeutic decision making.

Les futures missions d'exploration spatiale nécessitent le stockage d'ergols cryogéniques sur de longues durées.Sous l'effet d'entrées thermiques résiduelles, l'ergol se vaporise et le réservoir s'auto-pressurise, pouvant entraîner la rupture du réservoir pour des missionssuffisamment longues. Cette thèse s'intéresse à un système de contrôle, appelé Thermodynamic Venting System (TVS), reposant sur l'injection d'un jet sous-refroidi dans le réservoir.L'injection entraîne la condensation de la vapeur, la déstratification du bain liquide et donc une baisse de pression dans le réservoir.L'étude expérimentale a permis de mettre en place une technique originale d'isolation active générant une condition de paroi à flux de chaleurnet nul. Une base de données expérimentales, d'auto-pressurisation et de contrôle TVS, a été constituée avec ce nouveau dispositif.Elle a notamment permis de valider un modèle thermodynamique homogène permettant de prédire l'évolution de température et depression dans le réservoir. Ce modèle a été étendu pour de façon à décrire le comportement de tous les éléments constitutifs du système TVS.Un outil de dimensionnement du système complet a ainsi été mis en place.En le couplant à une plateforme d'optimisation, un système TVS optimal a pu être établi pour une mission de démonstration.Enfin une étude numérique a permis de mettre en évidence, pour notre cas d'étude, les faiblesses des modèles de changement de phase disponibles dans les solveursCFD commerciaux. Une modélisation prédictive du changement de phase a été mise en place dans un code de calcul recherche puis validée sur uncas académique 1D. / Future operations in space exploration require the ability to store cryogens for long duration. Residual heat loads induce cryogenic propellant vaporization andtank self-pressurization (SP), eventually leading to storage failure for long enough mission duration.This thesis focuses on a control strategy, called Thermodynamic Venting System (TVS), based on a recirculating liquid subcooled injection. The injection results inan ullage condensation, a liquid bath destratification and thus a tank pressure reduction.Experimentally, an original active insulation technique has been set up, yielding a net zero heat flux wall boundary condition. A data base ofself-pressurisation and TVS control experiments has been gathered with this new aparatus.It was used to validate an homogeneous thermodynamic model providing a fast prediction of tank temperature and pressure during control.This model has been extended to discribe the TVS system behaviour including all its components. This full system design tool has been coupled with an optimisationplatform and an optimal TVS design has been established for a demonstration mission.Furthermore, a numerical study has evidenced the weakness of commercial CFD software to simulate phase change, for TVS configuration.A predictive phase change formulation has been set up in a home-made software and validated on a 1D academic case.

L'activité des volcans andésitiques, tels que le Mont St Helens (États-Unis), Montserrat (Antilles) ou encore le Merapi (Indonésie), alterne entre des périodes relativement calmes, avec coulées de lave et formation d'un dôme, et des événements explosifs parfois très violents. Prévoir les transitions entre ces deux régimes est essentiel pour assurer la sécurité des populations voisines, mais demeure actuellement un vrai défi. Or les données expérimentales et les observations de terrain montrent que l'explosivité du magma est étroitement liée à son contenu en gaz. L'objectif de cette thèse est d'améliorer notre compréhension de l'évolution de ce contenu en gaz et de son influence sur l'activité volcanique, en nous appuyant sur des simulations numériques, l'analyse de données expérimentales ainsi que sur l'interprétation de données de déformation enregistrées au Merapi.Une part importante de ce travail réside dans le développement et l'amélioration de modèles d'écoulement en 2D pour prendre en compte le dégazage dans la partie supérieure du conduit, en régime transitoire. Nous présentons un modèle d'écoulement du gaz en temps qui tient compte des pertes en gaz aux bords du conduit et à sa sortie, selon les conditions présentes dans la roche encaissante et le dôme. Nous proposons également une adaptation des modèles de conduit permettant de coupler complètement l'écoulement du gaz avec celui du magma pour étudier l'évolution des conditions dans le conduit en régime transitoire. À partir de simulations de l'évolution du dégazage lors de l'emplacement d'un dôme, nous identifions les para-mètres contrôlant les pertes en gaz. Nos résultats montrent que ces pertes sont extrêmement sensibles à l'évolution de la perméabilité du magma et des gradients de pression autour du conduit en réponse au poids du dôme. La perméabilité du dôme a quant à elle peu d'influence. Au cours de la croissance du dôme, les pertes en gaz diminuent en profondeur. En haut du conduit, la pression du gaz augmente de quelques dizaines de MPa. Ces effets sont associés à une augmentation de l'explosivité du magma et de l'aléa volcanique en cas d'effondrement du dôme.Bien que la perméabilité du magma exerce un fort contrôle sur la perte de gaz, comme l'ont montré nos résultats, son évolution dans le conduit est peu contrainte. Les lois de perméabilité utilisées actuellement ne sont pas en accord avec l'ensemble des mesures réalisées sur des échantillons de magmas riches en silice. Dans le but d'améliorer notre compréhension du développement de la perméabilité dans le conduit, nous avons cherché à éclaircir le lien entre perméabilité, conditions d'écoulements, et caractéristiques géométriques du réseau de bulles connectées. Nous proposons une formulation du seuil de percolation, moment exact où le magma devient perméable compatible avec un grand nombre d'échantillons naturels et expérimentaux. Nous présentons aussi une nouvelle loi de perméabilité en accord avec la plupart des observations existantes, que nous avons intégrée à notre modèle 2D de dégazage. Nos résultats montrent qu'en fonction du nombre de bulles dans le magma et de la distribution de leurs tailles, l'importance des pertes en gaz et par conséquent les conditions d'écoulement dans le conduit varient d'effusives à explosives.Enfin, afin d'évaluer l'utilité des données de déformation pour suivre l'évolution des conditions d'écoulement, nous utilisons des modèles d'écoulement simples couplés à de la déformation élastique en 3D pour retrouver la déformation observée au sommet du Merapi peu avant l'éruption de 2006. Bien que ces modèles permettent de mieux comprendre les déplacements observés, le peu de données, associé à la complexité géologique et rhéologique du sommet, ainsi qu'à celle des processus physiques intervenant dans le conduit font qu'il est difficile de contraindre les conditions d'écoulement grâce à la déformation dans ce cas précis. / At silicic volcanoes, such as Mount St Helens (United States), Montserrat (British West Indies), or Merapi (Indonesia), periods of relative quiescence, with lava flows and dome emplacement, alternate with explosive, sometimes very violent events. Forecasting the effusive/explosive transitions, which is essential for the safety of nearby populations, remains currently a real challenge. However, experimental as well as field observations provide evidence that magma gas content is a major clue for understanding explosivity. This thesis, based on numerical simulations, experimental samples analysis, as well as on the interpretation of ground deformation data recorded at Merapi volcano, aims at improving our understanding of gas loss evolution, and its impact on the eruptive regime.A major part of this work consisted in developing and improving 2D axisymmetric conduit flow models for integrating gas loss in transient conditions. We provide a time-dependent model for gas flow in the upper conduit, that accounts for gas loss both at the conduit walls and at its top, depending on conditions in the surrounding rock and dome. We also propose an adaptation of conduit flow models allowing for full coupling between magma and gas flow in 2D that should be used to further investigate flow conditions evolution during transient regimes. From time-dependent gas flow simulations in the case of an effusive dome emplacement, we identify controlling parameters for gas loss. Our results provide evidence that gas loss is extremely sensitive to the evolution of magma permeability and of pressure gradients around the conduit due to dome loading, whereas, contrary to the common idea, dome permeability has almost no influence. Along with dome growth, gas loss decreases at depth, thus causing an increase in the magma gas content. At the top of the conduit, this results in an increase in gas pressure by a few tens of MPa, thus increasing the likelihood of magma explosivity and hazard in the case of a rapid decompression due to dome collapse.Although magma permeability plays a major role for gas extraction, as revealed by our results, its evolution within the conduit is poorly constrained. Currently used permeability laws fail in reassembling the whole dataset of permeability measurements from natural and experimental silicic samples. In order to improve our understanding of permeability development in the conduit, we worked on linking permeability and flow conditions with geometrical parameters that characterise the connected bubble network, based on experimental samples analysis. We propose an expression for the percolation threshold, i.e. the very moment when magma becomes permeable, that succeeds in classifying a wide dataset of natural and experimental samples. We also develop a new permeability law that reassembles most of the existing observations, and implement it within our gas flow 2D model. Results show that depending on the number of bubbles within the magma and on their size distribution, gas loss and then magma flow conditions evolve from effusive to explosive conditions.Eventually, we evaluate the applicability of monitoring flow conditions from observed ground deformation by using simplified conduit flow models, coupled with elastic deformation in 3D, to interpret ground deformation recorded in the near field at Merapi a few days before the 2006 eruption. Although conduit flow models provide important clues for interpreting observed displacements, the sparsity of field observations together with the complexity of the volcano summit geology, rheology and processes happening in the conduit make it very complex to constrain flow conditions from observed deformation.

Une des méthodes les plus populaires et largement utilisée pour l’amélioration des sols compressibles est la technique d'inclusions rigides verticales qui constitue une alternative intéressante aux techniques plus traditionnelles telles que le pré-chargement, la mise en place de drains verticaux, le remplacement du sol en place, etc.Le Projet National ASIRI (Amélioration du Sol par Inclusions Rigides) a permis d’apporter des réponses sur le comportement de cette technique de renforcement et d’établir des recommandations de dimensionnement et de réalisation. Celui-ci comprend un volet expérimental (expérimentations en vraie grandeur, modèles physiques réduits à 1g, études en chambre d’étalonnage) et un volet numérique. Cependant, les recommandations publiées par ASIRI sont limitées aux cas de chargement monotone. Pourtant, différents cas de structures sous chargement cyclique sont couramment rencontrés dans la pratique: chargement des vagues sur les structures offshore, force du vent, charge sismique, la charge de trafic, remplissage et vidange de réservoirs, chargement et déchargement des zones de stockage… Cela nécessite alors de comprendre le comportement sous chargement cyclique et/ou dynamique de ces ouvrages.Ce travail de thèse constitue une contribution à la compréhension du comportement de ce type d’ouvrage, et s’intéresse plus particulièrement à la modélisation des mécanismes qui se développent dans le matelas de transfert de charge lors d’un chargement cyclique vertical quasi-statique. Il se fait suivant deux approches complémentaires : Expérimentale : par des essais sur un nouveau modèle physique réduit tridimensionnel à 1g et à l’échelle 1/10ème, Numérique : par la modélisation numérique tridimensionnelle en milieu continu utilisant le logiciel FLAC3D.Dans un premier temps, le travail expérimental réalisé nous a permis d’étudier le comportement de cette technique de renforcement sous chargement monotone et cyclique mettant en œuvre la méthode de corrélations d’images. La modélisation physique réalisée nous a permis d’étudier l’influence de l’épaisseur de la plateforme de transfert de charge granulaire et des conditions aux limites (cas remblai/cas dallage).Dans un deuxième temps, les résultats expérimentaux nous ont servi comme une base de données pour la validation d’un modèle numérique en milieu continu. La première étape de la modélisation numérique consiste à valider un modèle numérique à la même échelle réduite du modèle physique. La validation de l’approche de modélisation numérique à échelle réduite nous a permis, dans une deuxième étape, de réaliser une étude paramétrique afin de déterminer l’influence des différents paramètres sur le comportement du système. / One of the most popular and widely used methods for soft soil improvement is the reinforcement using vertical rigid inclusions. It constitutes an interesting alternative to the other traditional techniques such as preloading, vertical drains, replacing soil etc...The national French project ASIRI (Amélioration du Sol par Inclusions Rigides) allowed to provide answers about the behavior of this technique and to establish recommendations for design and construction. It included an experimental part (full-scale, 1g physical model experiments and calibration chamber studies) and a numerical part. However, the recommendations established by ASIRI are limited to the case of monotonic loading. Yet, various cases of structures under cyclic loading are commonly encountered in practice: waves loading on offshore structures, wind strength, seismic load, traffic load, filling and emptying of tanks, loading and unloading of storage areas... This then requires the understanding of the behavior of this technique under cyclic and/or dynamic loading.The work done in this thesis constitutes a contribution to understanding the behavior of this type of structures, and is particularly interested in modeling the mechanisms that develop in the load transfer mattress under a quasi-static vertical cyclic loading . Two complementary approaches have been followed: Experimental: by performing experimental tests using a new 1g tridimensional physical model with a scale factor of 1/10, Numerical: by performing numerical tridimensional simulations in a continuous media using the software FLAC3D.At first, the experimental work has allowed us to study the behavior of this reinforcement technique under monotonic and cyclic loading implementing the digital image correlation method. A parametric study performed with the physical modeling allowed us to study the influence of the granular load transfer platform (LTP) and the boundary conditions (presence of a rigid slab at the surface of the LTP).Secondly, the experimental results have served as a database for the validation of a numerical model in continuous medium. The first step of the numerical modeling is to validate a numerical model at the same reduced scale of the physical model. The validation of the numerical modeling approach at the reduced scale allowed us in a second step to perform a parametric study in order to determine the influence of different parameters on the behavior of the system.