Ce travail de thèse porte sur l'étude du comportement de l'interface sol-structure sous sollicitations cycliques. Pour cela, un important programme expérimental à la sonde-pieu, mise en place dans des massifs de sable siliceux en chambre d'étalonnage, a été réalisé. On s'est intéressé, particulièrement, à l'évolution du frottement latéral et de la résistance enpointe à grand nombre de cycles (100 000 cycles), en faisant varier l'état initial du massif (état de densité et contrainte de consolidation) et les paramètres de chargement (amplitude du déplacement cyclique). Lors des essais à déplacement contrôlé, un renforcement important du frottement latéral à grand nombre de cycles a été observé pour les faibles valeurs de résistance en pointe initiale. Ce renforcement n'a quasiment pas été observé avant par les auteurs, car l'accent a été mis sur la phase de dégradation (autour de 1000cycles). Il est attribué à une forte dilatance partiellement empêchée. En revanche, pour les fortes valeurs de résistance en pointe, le renforcement à grand nombre de cycles diminue considérablement. Cette diminution est liée à l'effet des particules fines créées lors du fonçage. Celles-ci jouent le rôle de cimentation/scellement de la surface latérale de la sonde pieu en diminuant sa rugosité. Par ailleurs, des essais à force contrôlée ont été réalisés afin d'étudier la stabilité des pieux. Enfin, des lois empiriques d'évolution du frottement latéral et de la résistance en pointe ont été proposées afin de reproduire les évolutions observées expérimentalement. Ces lois d'évolution ont été intégrées dans un modèle de calcul de pieu sous chargement cyclique de type t-z. Les premières simulations effectuées montrent un bon accord entre le modèle et les résultats expérimentaux à petit nombre de cycles

La caractérisation de la morphologie de surface des joints rocheux est d'un intérêt particulier pour une évaluation correcte de la résistance au cisaillement des discontinuités rocheuses. D'un point de vue géométrique, cette dernière est particulièrement caractérisée par une forme très irrégulière des aspérités, définissant ainsi une géométrie polygonale des profils des joints. Dans ce cadre, une méthode spécifique d'analyse spectrale et d'identification morphologique adaptée au problème de la caractérisation de la géométrie irrégulière des discontinuités rocheuses a été mise au point. Cette méthode repose sur une décomposition en série à partir de fonctions élémentaires triangulaires (méthode dénommée ATM). Sur le plan expérimental, les relevés des profils des épontes des joints rocheux ont été effectués au moyen d'un rugosimètre 3D en laboratoire, et d'un profilographe à règle pour les mesures in situ, mis au point dans le cadre de cette recherche. Nos travaux ont porté en premier lieu sur la caractérisation et la description quantitative de la morphologie des joints rocheux à différentes échelles. Diverses simulations sont présentées qui démontrent l'intérêt d'une telle approche notamment, en comparaison avec la méthode d'analyse spectrale par les transformées de Fourier. En second lieu, ont tire profit des informations spectrales pour transposer la géométrie réelle tridimensionnelle en image bidimensionnelle qualitativement équivalente du point de vue du comportement au cisaillement (dilatance). Ces travaux ont conduit à la mise au point d'un modèle numérique d'analyse et de synthèse spectrale, dénommé MORPHEE. Les caractéristiques et les performances de ce nouvel outil sont exposées et illustrées par l'application de la méthode d'analyse ATM à l'étude de la morphologie des profils de discontinuités rocheuses à différentes échelles. La méthode d'analyse a fait l'objet d'une validation expérimentale par comparaison avec des échantillons soumis à des essais de cisaillement en laboratoire.

La modélisation de la déformation, de la rupture et de la fracturation des milieux rocheux représente un enjeu majeur dans de nombreux domaines scientifiques et pratiques, notamment pour l'exploitation des réservoirs fracturés. Cependant, cette modélisation se heurte à la mauvaise connaissance des lois constitutives, ce qui rend particulièrement importantes les approches expérimentales et en particulier la modélisation physique, qui est au coeur de ce travail.<br /> Nous avons développé une technique originale d'une telle modélisation, basée sur un nouveau matériau Crack1 à l'échelle physique par rapport à une roche réservoir typique: le calcaire. Un programme étendu d'expérimentation avec les modèles de Crack1 a été réalisé avec un dispositif de chargement poly-tridimensionnel. Ce dispositif combine à la fois la simplicité et l'efficacité des solutions réalisées pour la réduction et l'élimination complète de la friction aux bords des modèles. Les résultats principaux peuvent être résumés de la façon suivante:<br />1.Des réseaux de diaclases ont été reproduits pour la première fois dans des conditions de chargement homogènes.<br />2.Les diaclases se forment, dans nos expériences, sous une compression triaxiale. Elles ne sont donc pas des fractures en Mode I. <br />3.Les figures fractographiques visibles sur les surfaces des diaclases naturelles et celles obtenues dans les modèles sont très similaires, ce qui suggère que la similarité physique est assurée pas seulement à macro-échelle, mais aussi à micro-échelle.<br />4.L'espacement S entre les diaclases ne dépend pas de l'épaisseur du modèle , contrairement au concept de « saturation » largement adopté dans la littérature et les modèles de réservoirs. Il a été démontré que S est contrôlé par l'état des contraintes imposé et la déformation accumulée du modèle.<br />5.Nous avons pu également reproduire dans des modèles les couloirs fracturés dont la formation est contrôlée par la rigidité des contacts aux limites de l'unité mécanique affectée par la fracturation.<br />6.L'augmentation de la pression moyenne appliquée au modèle résulte en un changement du style de la fracturation, qui évolue de manière continue des diaclases vers des fractures «obliques » conjuguées.<br />7.L'analyse des conditions de cette transition à partir de la théorie de la bifurcation de la déformation montre qu'elle est associée à la réduction du facteur de dilatance β de Crack1 avec l'augmentation de la pression p (ce type de dépendance β (p) est également connu pour les roches réelles).<br />Les résultats obtenus confirment donc que les lois constitutives contrôlent directement la fracturation. Des études expérimentales plus larges doivent désormais être menées pour contraindre ces lois dans toute leur complexité et avec une précision nécessaire.

Nous étudions expérimentalement l'effondrement d'une colonne granulaire dans un liquide visqueux. Contrairement au cas sec, la morphologie des dépôts n'est principalement plus contrôlée par le rapport d'aspect initial du tas mais par la fraction volumique initiale de la masse granulaire. Deux régimes différents sont identifiés selon l'empilement initial. L'empilement lâche donne lieu à des dépôts minces et longs et la dynamique est rapide. Une surpression du liquide est mesurée sous de la colonne. Pour l'empilement dense, l'étalement final est deux fois moindre, le mouvement est lent et une dépression interstitielle est mesurée. Ces observations suggèrent que la dynamique de l'effondrement granulaire dans un fluide est fortement affectée par le comportement de la dilatance du milieu granulaire.Nous développons ensuite un modèle théorique basé sur des équations diphasiques moyennées dans l’épaisseur prenant en compte les mécanismes de dilatance. L’étude dimensionnelle de notre modèle permet de montrer que l’effondrement d’une colonne est contrôlé par trois paramètres sans dimension : le rapport d’aspect de la colonne, la fraction volumique initiale, et le nombre de grains dans l’épaisseur. On montre également que le temps caractéristique met en compétition le frottement visqueux et la gravité.De ce modèle, nous développons un algorithme de résolution lagrangien. Cette approche, grossière mais robuste, permet d’implanter s sans trop de difficulté. Le code est validé sur des configurations simples sur plan incliné avant de simuler l’effondrement de colonnes granulaires immergées dans la même gamme de paramètres que nos expériences. / Nous étudions expérimentalement l'effondrement d'une colonne granulaire dans un liquide visqueux. Contrairement au cas sec, la morphologie des dépôts n'est principalement plus contrôlée par le rapport d'aspect initial du tas mais par la fraction volumique initiale de la masse granulaire. Deux régimes différents sont identifiés selon l'empilement initial. L'empilement lâche donne lieu à des dépôts minces et longs et la dynamique est rapide. Une surpression du liquide est mesurée sous de la colonne. Pour l'empilement dense, l'étalement final est deux fois moindre, le mouvement est lent et une dépression interstitielle est mesurée. Ces observations suggèrent que la dynamique de l'effondrement granulaire dans un fluide est fortement affectée par le comportement de la dilatance du milieu granulaire.Nous développons ensuite un modèle théorique basé sur des équations diphasiques moyennées dans l’épaisseur prenant en compte les mécanismes de dilatance. L’étude dimensionnelle de notre modèle permet de montrer que l’effondrement d’une colonne est contrôlé par trois paramètres sans dimension : le rapport d’aspect de la colonne, la fraction volumique initiale, et le nombre de grains dans l’épaisseur. On montre également que le temps caractéristique met en compétition le frottement visqueux et la gravité.De ce modèle, nous développons un algorithme de résolution lagrangien. Cette approche, grossière mais robuste, permet d’implanter s sans trop de difficulté. Le code est validé sur des configurations simples sur plan incliné avant de simuler l’effondrement de colonnes granulaires immergées dans la même gamme de paramètres que nos expériences.

Les cavités salines représentent une technique prometteuse de stockage massif d’énergie, notamment pour les énergies renouvelables dont la production est par nature intermittente et imprévisible. Historiquement utilisées pour le stockage saisonnier d’hydrocarbures (méthane, pétrole...), les cavités salines sont aujourd’hui sollicitées pour le stockage de nouveaux fluides (hydrogène, dioxyde de carbone...) avec des scenarii plus exigeants. Les méthodes de dimensionnement des cavités doivent être mises à jour pour répondre aux nouveaux défis de la transition énergétique.Cette thèse propose une nouvelle méthodologie de dimensionnement des cavités salines, basée sur le développement d’un nouveau modèle constitutif pour le sel gemme incluant des critères de dilatance et de traction. Ce nouveau modèle permet d’ajuster avec un unique jeu de paramètres de nombreux essais de laboratoire différents, en particulier courts et longs.Des simulations couplées thermo-mécaniques de cavités, remplies de méthane ou d’hydrogène, et du sel gemme environnant sont réalisées pour différents scenarii d’exploitation, classiques ou se rapprochant des nouveaux besoins liés à la transition énergétique. On étudie en particulier les effets de la durée et de l’amplitude des cycles, du débit d’injection ou de soutirage. Les résultats obtenus avec la nouvelle méthodologie sont comparés avec ceux de la méthodologie classique. / Salt caverns are a promising technique for massive energy storage, especially in the case of the intermittent and unpredictable renewable energy. Historically used for seasonal storage of hydrocarbons (methane, oil...), they are potentially operated with increasingly demanding scenarios for the storage of other fluids (hydrogen, carbon dioxide...). Design methods need to be updated to rise to the new challenges of the energy transition.This thesis proposes a new methodology for salt cavern design, based on the development of a new rheological model including a dilatancy and a tensile criteria. This new model allows to fit numerous different laboratory tests with a single parameter set, in particular short- and long-term tests.Thermo-mechanical numerical simulations of caverns, filled with either methane or hydrogen, and the surrounding rock salt are performed under various cycling scenarios which are classical or closer to the needs associated with the energy transition. Effects of cycle duration, amplitude and mass flow are especially investigated. Results obtained with the new and the classical methodologies are compared.

Le présent travail à caractère théorique et numérique étudie le comportement mécanique des interfaces sol granulaire-structure ainsi que son application aux pieux isolés sous chargement statique. A partir du modèle de comportement élastoplastique d'interface bidimensionnel "MEPI 2D", un modèle tridimensionnel dit "MEPI 3D" est développé afin de mettre en évidence l'effet des contraintes de confinement tridimensionnelles régnant à l'interface au cours du cisaillement. Ces modèles considèrent un écrouissage déviatorique, prennent en compte l'état caractéristique (notion de contractance-dilatance) et intègrent l'état critique aux grands déplacements tangentiels. Les prévisions du comportement de l'interface à l'aide de MEPI 3D sont satisfaisantes. Parallèlement à cette contribution rhéologique, une méthodologie générale de modélisation numérique par éléments finis d'un essai de chargement de pieu est proposée. Les modèles MEPI 2D et MEPI 3D sont implantés dans le code de calcul aux éléments finis CESAR-LCPC. L'analyse numérique par CESAR-LCPC d'un pieu modèle en chambre d'étalonnage (Cermes, type 1), de pieux modèles sur site réel (Imperial College Pile, type 2) et d'un pieu en vraie grandeur (caisson du LCPC, type 3) est menée considérant les modélisations axisymétriques (MEPI 2D, MEPI 3D) et tridimensionnelles (MEPI 3D). En comparant les prévisions numériques aux données expérimentales, il est montré que le comportement des pieux est correctement reproduit. Par ailleurs, la comparaison de l'influence des modèles MEPI 2D et MEPI 3D est abordée.

Cette thèse concerne la modélisation du comportement mécanique de poudres ductiles au cours de leur mise en forme par le procédé de compression à froid. L’approche utilisée consiste à modéliser la poudre par un échantillon numérique de 50 particules sphériques distribuées aléatoirement. Chaque particule est maillée en éléments finis avec un comportement élasto-plastique et interagit par contact avec ses voisines. Sous l’effet du chargement, les particules se déforment et subissent de grandes transformations.La méthode développée consiste à imposer aux frontières de l'échantillon numérique des conditions aux limites en déplacement ou en force afin de simuler une réponse continue équivalente en termes de contrainte-déformation pour plusieurs chemins de chargement depuis les faibles valeurs de densité relative jusqu'aux fortes valeurs. Outre le fait que la méthode permet d’appliquer par la simulation des chemins de chargement inaccessibles aux techniques expérimentales, elle fournit de nombreuses indications sur la phénoménologie du comportement aux frontières de l'échantillon en relation avec la morphologie des particules déformées plastiquement et l’évolution des surfaces de contact inter-particules.Cette réponse obtenue sera étudiée au moyen des notions de l'élasto-plasticité classique, à savoir la surface de charge et la loi d'écoulement. L'accent est mis sur l’analyse de l’écoulement plastique et l'existence du potentiel plastique. Les résultats montrent que la direction d’écoulement est approximativement unique pour les états de contrainte éloignés du point de chargement. En revanche, la direction d’écoulement cesse d’être unique à proximité de point de chargement, ce qui révèle le caractère anisotrope de l’écrouissage. Ce comportement, propre aux poudres ductiles, est intimement lié à la déformation des grains. La dégradation anisotrope des propriétés mécaniques résultant du mécanisme de dilatance est également étudiée. / This thesis is concerned with the modelling of the mechanical behaviour of ductile powders during cold compaction process. A numerical method was implemented, in which a powder sample was assimilated to an assembly of 50 meshed particles in a cubic box which was submitted to compressive loadings simulated by means of the finite element method. Each particle was meshed and assigned an elastic-plastic constitutive behaviour. Particles interacted through mechanical, frictional contact and underwent large transformation during loading.The method developped in the thesis consisted in imposing boundary conditions as displacements or forces on the bounding walls of the numerical sample. The mechanical response of the sample was averaged to obtain the equivalent response of a continuum in terms of stress and strain. Various loading paths, including loading paths that are technically not attainable by experimental means, were applied up to different values of relative densities, from low (60%) to high values (98%). Much information on the phenomenological behaviour of the numerical sample could be obtained, in relation with the morphology of deformed particles and the evolution of contact surfaces.The results were interpreted based on the concepts of classical elasto-plasticity, i.e., yield surface and flow rule. The main focus of the study was the analysis of incremental plastic flow and the existence of a plastic potential. Results showed that the flow rule postulate, based in the plastic potential, could be considered as valid for stress states relatively far from the loading point. In the vicinity of the loading point, the direction of the plastic strain increment vector ceases to be unique. This behaviour, which is an original feature of ductile powders mechanical behaviour, was attributed to the anisotropic strain-hardening processes at stake, related to the anisotropic formation of contact surfaces between particles. The drop in mechanical properties related to dilatancy was also studied.

L'optimisation des systèmes de forage nécessite une meilleure compréhension des vibrations indésirables comme le stick-slip. Ce phénomène vibratoire, qui affecte principalement les outils PDC (Polycristalline Diamond Compact), met en péril l'intégrité des équipements de forage et réduit considérablement la vitesse de pénétration de l'outil. Plusieurs travaux ont été menés ces dernières années pour déterminer ses origines. Les observations réalisées en fond de puits montrent que ces oscillations s'accompagnent systématiquement d'une baisse du couple à l'outil en fonction de sa vitesse de rotation. De nombreux groupes de recherche attribuent cette baisse de performance à l'occurrence du stick-slip.L'objectif de ce travail est de développer un modèle élémentaire de coupe qui permet d'analyser l'effet de la vitesse de coupe sur la forabilité des roches dans des conditions opératoires réalistes. Dans le cadre de cette thèse, nous avons réalisé une série d'essais de coupe en utilisant des taillants et des outils à échelle réelle dans trois roches de propriétés hydromécaniques différentes, et ceci à pression atmosphérique et sous pression de fluide. Les essais réalisés à pression atmosphérique montrent que les efforts élémentaires de forage augmentent avec la vitesse de coupe. Sous pression de boue, cet effet dépend largement de la perméabilité de la roche. En effet, nous avons observé que l'effet de la vitesse est relativement faible dans les formations de faible et de moyenne perméabilité sous pression de boue de 20 MPa. En revanche, cet effet augmente d'un ordre de grandeur dans les roches très perméables.Afin de comprendre ces observations, nous avons développé un modèle hydromécanique d'interaction taillant-roche construit à partir de la théorie de la poroélastoplasticité. D'abord, le problème est résolu analytiquement en s'inspirant des travaux existants. Par la suite, nous avons apporté une résolution numérique aux éléments finis des équations de la promécanique appliquées à la coupe des roches sous pression de boue. Les deux modèles montrent que le phénomène de dilatance génère une baisse de la pression de pore qui augmente la résistance de la roche au forage. Cette chute de pression dépend de la vitesse de coupe ainsi que des caractéristiques hydrodynamiques de la roche. Les résultats théoriques ont été comparés aux nombreux résultats expérimentaux obtenus dans le cadre de ce travail. / The optimization of the drilling practice requires a better understanding of drillstring harmful vibrations such as stick-slip. This form of torsional vibrations is a typical problem of PDC (Polycristalline Diamond Compact) drillbits. It can reduce the rate of penetration drastically and can raise fatigue of the drilling devices. Many attempts were carried out in the last years in order to determine the causes of stick-slip phenomenon. Field observations show that torque on bit decreases as a function of bit velocity during stick-slip oscillations. Hence, it is widely believed that this decreasing relationship is the root cause of stick-slip.The purpose of this work is to examine cutting speed influence on rock drillability as a function of operating conditions and hydromechanical properties of the drilled formation. For this, a set of drilling tests was performed in three sedimentary rocks of different permeability using a full scale PDC drillbit and a single PDC cutter. In the first step, dry tests were carried out at atmospheric pressure. As previously observed in literature, single-cutter tests showed that drilling forces increase with cutting velocity. In a second step, we performed the same experiments at 20 MPa bottom-hole pressure. It appears that rate effect on cutting forces in the medium and low-permeability rocks is relatively low. By contrast, rate effect in the highly permeable rock increases by one order of magnitude in comparison with dry experiments.In order to understand this phenomenon, a steady state solution of the cutting model is derived in the framework of the theory of poroelastoplasticity. The problem is firstly solved analytically using some assumptions derived from previous works. Then, a numerical resolution based on finite element method is presented to solve the fully coupled problem ensuring the satisfaction of poro-material physics basic equations. Using these two different approaches, we show that pore pressure in shear-dilatant rocks decreases as a function of cutting velocity depending on rock permeability and interstitial fluid properties. This change has a hardening effect resulting in an increase of rock drilling resistance. Comparison between theory and experience shows good agreements.

Nous nous intéressons à l'écoulement de milieux granulaires cohésifs soumis à des vibrations horizontales en vue d'une application au remplissage des moules de presse. L'effet des vibrations est analysé pour un grain ou un milieu granulaire déposé sur un plan incliné et un milieu granulaire contenu dans une conduite verticale de section constante. Pour un grain polyédrique, trois régimes sont identifiés : stick, stick-slip et slip. Un modèle basé sur le bilan des forces appliquées à un grain a permis de préciser les limites de ces trois domaines. La vitesse des grains tend vers une asymptote lorsque l'inclinaison est inférieure à 0,7 fois l'angle de frottement entre les grains et le plan. Par ailleurs, l'inertie du grain conduit à une réduction de l'amplitude de son oscillation lorsque l'accélération augmente. La comparaison modèle-expérience permet d'estimer le coefficient de frottement. Pour un milieu granulaire de faible épaisseur, nous retrouvons les mêmes tendances. Les milieux granulaires épais sont le lieu de forts cisaillements provoquant leur étalement très rapide. En conduite verticale, une compétition s'établit entre les déplacements horizontal et vertical des particules. Selon l'inertie, les vibrations appliquées pourront conduire soit à une dilatance soit à une compaction. Cette dernière provient du déplacement de la conduite et du milieu granulaire en sens opposé pendant une partie d'une période de vibration. Si le phénomène de compaction atteint le centre de la conduite, il se forme des arches qui bloquent l'écoulement. Les observations faites à l'échelle mésoscopique permettent d'interpréter les vitesses d'écoulement mesurées à l'échelle macroscopique. / With the target to improve the feeding step during the moulding process, we study the flow of cohesive granular matter submitted to vibration. Vibration effects are first analyzed on the sliding motion of a single particle on inclined plane and also on granular matter sample deposited on an inclined plane or in a vertical funnel. For a single particle, three regimes are identified: stick, stick-slip and slip regimes. A simple model based on the movement equation of one grain allows to determinate the limits between the three regimes. The grain velocity reaches asymptote when the plane inclination is below 0.7 times the friction angle between the grain and the substrate. Otherwise, when the acceleration increases, the grain inertia leads to reduce the transverse amplitude of the grain oscillation motion. The comparison between experimental and numerical results allows to estimate a value for the friction coefficient parameter. For a granular sample, we find the same tendencies. Thick granular layers are submitted to high shear, causing a fast spreading. In vertical funnel, there is a competition between vertical and horizontal motions. Depending on granular inertia, submitted vibrations could lead to dilatancy due to the shear or lead to compaction. This compaction is due to an opposite motion direction between the sample and the funnel during one part of the period vibration. If the compaction reaches the center of the funnel, arches are formed and jammed the flow. Observations realized at mesoscopic scale allow to interpret flow velocities measured at macroscopic scale.

Ce travail de thèse porte sur l'étude du comportement de l'interface sol-structure sous sollicitations cycliques. Pour cela, un important programme expérimental à la sonde-pieu, mise en place dans des massifs de sable siliceux en chambre d'étalonnage, a été réalisé. On s'est intéressé, particulièrement, à l'évolution du frottement latéral et de la résistance enpointe à grand nombre de cycles (100 000 cycles), en faisant varier l'état initial du massif (état de densité et contrainte de consolidation) et les paramètres de chargement (amplitude du déplacement cyclique). Lors des essais à déplacement contrôlé, un renforcement important du frottement latéral à grand nombre de cycles a été observé pour les faibles valeurs de résistance en pointe initiale. Ce renforcement n'a quasiment pas été observé avant par les auteurs, car l'accent a été mis sur la phase de dégradation (autour de 1000cycles). Il est attribué à une forte dilatance partiellement empêchée. En revanche, pour les fortes valeurs de résistance en pointe, le renforcement à grand nombre de cycles diminue considérablement. Cette diminution est liée à l'effet des particules fines créées lors du fonçage. Celles-ci jouent le rôle de cimentation/scellement de la surface latérale de la sonde pieu en diminuant sa rugosité. Par ailleurs, des essais à force contrôlée ont été réalisés afin d'étudier la stabilité des pieux. Enfin, des lois empiriques d'évolution du frottement latéral et de la résistance en pointe ont été proposées afin de reproduire les évolutions observées expérimentalement. Ces lois d'évolution ont été intégrées dans un modèle de calcul de pieu sous chargement cyclique de type t-z. Les premières simulations effectuées montrent un bon accord entre le modèle et les résultats expérimentaux à petit nombre de cycles / We study in this present work the behavior of the soil-structure interface under large number of cycles (100 000 cycles). An important program with a probe-pile jacked into the beds of sand was carried out in calibration chamber. We interested particularly on the evolution of local skin friction and tip resistance at different initial state of beds (initial density and confining pressure) and the parameters of the load (cyclic displacement amplitude). We conducted an important part of the experimental work under axial cyclic displacement controlled test which allows studying the large number of cycles. After the degradation phase, already observed by several authors (about 1000 cycles), an important re-increase in the skin friction at large cyclic number was observed at low values of initial tip resistance.This re-increase is attributed to a high dilation of sand around the probe. However, for high values of initial tip resistance there is almost not re-increase in the skin friction. This is due tothe formation of a shear band around the probe made of crushed sand. These reduce theroughness of the probe by cementing/sealing. In addition to the displacement-controlled tests, we conducted the load-controlled tests in order to study the stability of the pile. The results showed a good agreement with the displacement-controlled tests. Finally, empirical model of evolution of skin friction and tip resistance have been proposed in order to reproduce the experimental results. This model was incorporated into a computational model of pile under cyclic axial loading. The first simulations showed a good agreement with the experimental results at low number of cycles