L’analyse du comportement des pieux est un problème complexe du fait de la diversité des phénomènes qui gouvernent le comportement du sol et en particulier celui se trouvant au voisinage du pieu. Ce dernier dépend particulièrement du procédé d’installation du pieu qui peut parfois engendrer des déformations de très grande amplitude dans le sol entre autre phénomènes. L’objectif de ce travail est de mettre en place un outil de modélisation pour évaluer le comportement des pieux sous chargements axiaux en développant des modèles représentant les phénomènes physiques pertinents et de les intégrer numériquement dans un logiciel aux éléments finis utilisant des algorithmes robustes et efficaces. Pour pouvoir modéliser la phase d’installation statique ou dynamique d’un pieu quelques problématiques doivent être considérées. Premièrement, étant donné que pendant l’installation d’un pieu des déformations finies (non-infinitésimal) ont lieu au niveau de l’interface, une formulation eulérienne avec un taux logarithmique des déformations a été adoptée pour prendre en compte le fait que l’hypothèse des déformations infinitésimales n’est plus valable. En plus, le modèle constitutive doit prendre en compte le comportement physique du sol lorsqu’il est soumis à des déplacements d’une magnitude élevée. Le dernier est constitué, entre autres, par le phénomène de l’écrasement des grains ce qui influence beaucoup le comportement volumique du sol et finalement sa résistance au cisaillement ou en d’autres termes le frottement mobilisé. Ce phénomène a été modélisé en introduisant une variable d’écrouissage supplémentaire au modèle de comportement élastoplastique de l’ECP. Les critères d’admissibilité thermodynamique ont été vérifiés pour le modèle constitutive original ainsi que pour le modèle révisé. Des simulations numériques ont été faites pour les deux types d’installation, monotone et pseudo-dynamique (cyclique) et les résultats ont été analysés en détail. Finalement, la dégradation de la résistance au cisaillement au niveau du fût des pieux est un phénomène typique qui se produit pendant le chargement cyclique des fondations composées par des pieux. Le modèle constitutive pour pouvoir bien reproduire ce phénomène n’est pas simple à définir car le chemin de chargement suivi par l’interface sol-pieu est très dépendant du comportement volumique du sol qui à son tours dépend de l’histoire de chargement et des conditions aux limites du problème. Une étude détaillée de toutes les composantes du comportement du pieu pendant ce type de chargement a été effectuée afin de mettre en évidence l’influence de l’histoire de chargement sur la résistance au cisaillement et l’apparition des phénomènes tels que la fatigue du frottement. / The analysis of pile behavior is a complex problem due to the diversity of the phenomena governing the soil behavior and particularly that of the neighboring soil. The objective of this work was to develop a valid modeling tool to evaluate piles’ behavior under axial loads by developing a pertinent mechanical model supported in a robust finite element program which would successfully reproduce the soil behaviour under extreme monotonic and cyclic shear strain. This is done to allow for the numerical modelling of the installation procedure of pile foundations and continued loading of high amplitude cyclic paths. In order to model the installation phase of a monotonic, jacked or dynamic pile foundation some issues must be addressed. Finite deformations take place whilst the pile is put in place requiring an adjustment in the mechanical formulation of the model at the interface level to take into account that the small deformations (rotations and strain) hypothesis is no longer valid. Moreover, the constitutive model must take into account the physical behaviour of the soil when subjected to high order of magnitude displacements. This includes the phenomenon of grain breakage, also referred to as particle crushing, which greatly influences the volumetric behaviour of soil as thus reflecting of shear stress mobilization. The elastoplastic ECP model has therefore been enhanced by introducing an internal variable taking into account the breakage mechanism. The thermodynamic admissibility criteria are verified for the original and revised constitutive models. Both monotonic and pseudo-dynamic installation procedures were numerically simulated and the results thoroughly analysed. Finally, the cyclic shear resistance degradation at the pile shaft is a commonly occurring phenomenon during continued cyclic loading of pile foundation (friction fatigue). The constitutive modelling of this phenomenon, however, is not a straightforward matter. The stress path followed by the thin layer at the soil-pile interface level is known to be directly related to the volumetric behaviour due to the boundary conditions of the problem. A comprehensive analysis of all the components of the behaviour of soil during this stage was object of study in this work.

Durant la phase de construction des puits pétrolier, le ciment est coulé dans l'espace annulaire entre le cuvelage et la roche environnante. La gaine de ciment a pour but de tenir le cuvelage, garantir l'étanchéité des puits pétroliers, ou des réservoirs de stockage de CO2, et de protéger le cuvelage de la corrosion. Au cours de la vie du puits, cette gaine de ciment est soumise tout le long du puits à des sollicitations thermiques et mécaniques qui varient au cours du temps et qui peuvent modifier ses propriétés et altérer son étanchéité. L'objet de cette thèse est d'étudier l'effet de la température et des contraintes sur les propriétés mécaniques de la pâte de ciment en cours de prise et de la pâte de ciment durcie. L'approche est basée sur une étude expérimentale qui combine des essais calorimétriques, des mesures de vitesse des ondes et des essais oedométriques avec le système des cellules STCA (Slurry To Cement Analyzer) sur le ciment en cours de prise ainsi que des essais de compression uniaxiale et triaxiale sur la pâte de ciment durcie. Les résultats expérimentaux ont montré que la température et la pression accélèrent la cinétique d'hydratation et que la température affecte significativement les propriétés élastiques du matériau tandis que la pression ne les influence pas pour une gamme de pression limitée à 20MPa. Une modélisation de la cinétique d'hydratation associée à une technique d'homogénéisation est proposée afin d'interpréter les essais. On a mis aussi en évidence que lorsque la prise se fait sous contraintes mécaniques, des déformations irréversibles peuvent se développer dans la pâte de ciment, ce qui peut conduire à la formation d'un micro-annulaire entre la gaine ciment, le cuvelage et la formation géologique. Le comportement différé de la pâte de ciment durcie a été étudié à partir d'essais de fluage uniaxiaux et d'essais de compression isotrope. Les essais ont mis en évidence que le fluage sous chargement uniaxial est plus important pour un ciment hydraté à une température plus élevée, ce qui est attribué à une microstructure plus hétérogène. Un modèle visco-endommageable permettant de reproduire les phases de fluage primaire et tertiaire a été développé. Les essais de compression isotrope drainés et non-drainé isothermes sous forte contrainte ont montré un comportement différé avec hystérésis lors de cycles décharge-recharge. Ces essais ont été analysés à partir d'un modèle de comportement poro-visco-plastique. Le comportement élastoplastique à court terme a été abordé à l'aide des essais triaxiaux sous chargement déviatorique drainé. Ces essais ont mis en évidence que la température d'essai affecte fortement la surface de charge initiale et les déformations tandis qu'elle ne modifie pas significativement la surface de rupture. Un modèle de plasticité avec une surface de charge fermée et un écrouissage dépendant des déformations plastiques accumulées a été développé pour décrire ces essais. Enfin, une étude préliminaire sur les effets des cycles mécaniques et thermiques a été menée. Des cycles thermiques ne dépassant pas la température d'hydratation ne semblent pas affecter les propriétés mécaniques du matériau. Cependant, une dégradation très rapide avec le nombre de cycles mécaniques a été mise en évidence lorsque la contrainte dépasse 60% de la résistance en compression simple / During the construction phase of oil wells, a cement slurry is pumped into the annular space between the casing and the rock. The cement sheath aims to support the casing, provide zonal isolation in the well and reservoirs of CO2 storage and protect the casing against corrosion. During the life of the well, the cement is submitted to various thermal and mechanical solicitations along the well that can modify its mechanical properties and damage its sealing performance. The aim of this thesis is to study the effect of temperature and stresses on the mechanical properties of cement paste during hydration and in the hardened state. The used approach is based on an experimental study that combine the calorimetric tests, waves velocity measurement and oedometric tests in STCA system (Slurry To Cement Analyzer) on cement paste during hydration as well as the uniaxial and triaxial compression tests on hardened cement paste. The experimental results showed that temperature and pressure accelerate the kinetics of hydration. Temperature affects significantly the elastic properties of the material whereas the pressure does not modify them for a range of pressure limited to 20MPa. A hydration kinetics modelling associated to a homogenization method is used to interpret the tests. It is shown that for hydration under stress, the irreversible strains can evolve in the cement paste and conduct to the formation of a micro-annular between the cement sheath, the casing and geological formation. The time-dependent behaviour of hardened cement paste is studied using creep tests under uniaxial loading and also from the results of isotropic compression tests. The results show that the uniaxial creep is more important for cements hydrated at higher temperatures, which is attributed to a more heterogeneous microstructure. A visco-damaged model allowing to reproduce the primary creep and tertiary creep is developed and calibrated. Isothermal isotropic drained and undrained compression tests show a time-dependent behaviour with hysteresis during unloading-reloading cycles. These tests are analyzed on the basis of a poro-visco-plastic model. The elastoplastic behaviour in short terms is studied from triaxial tests under drained deviatoric loading. These tests bring to light that the test temperature affects highly the initial yield surface and the strains but it does not significantly modify the failure surface. A model of plasticity with a closed yield surface and hardening, depending on the accumulated plastic strains is developed to describe these tests. Finally, a preliminary study on the effect of mechanical and thermal loading cycles is performed. The thermal loading cycles with temperatures lower than the hydration temperature seem to do not affect the mechanical properties of the material. The mechanical loading cycles show a rapid degradation with the number of loading cycles when the axial stress exceeds 60% of the uniaxial strength

Le sujet de thèse que nous abordons concerne l'identification des lois de comportement élastoplastiques anisotropes en vue de leur utilisation pour la simulation numérique des procédés de mise en forme par déformation plastique de tôles minces métalliques d'emboutissage. Nous avons essentiellement contribué à la définition et à la mise en oeuvre de stratégies et techniques d'identification des lois de comportement à partir d'essais expérimentaux. Les essais classiquement utilisés pour l'identification des modèles : la traction simple dans les axes et hors axes, la traction plane et le gonfelement hydraulique sont présentés et analysés du point de vue homogénéité des déformations pour aboutir à la relation contrainte-déformation à partir des mesures globales force-déplacement. L'identification des modèles de comportement à partir des essais inhomogènes nécessite une simulation numérique couplée avec une méthode d'optimisation (méthode du Simplexe) pour minimiser l'écart entre les résultats expérimentaux et la réponse du modèle calculée par une méthode d'élements finis. Nous avons identifié des lois de comportement élastoplastiques anisotropes avec écrouissage isotrope. Ces lois sont basées en particulier sur le choix d'une ou deux fonctions "contrainte équivalente" définissant le critère de plasticité et le potentiel plastique (cadre de la normalité non associée) ayant la même structure que la fonction seuil. Plusieurs critères quadratiques et non quadratiques sont alors utilisés. Nous avons développé une technique spécifique d'analyse de sensibilité des essais par rapport aux differents paramètres à identifier de la loi de comportement. Cette technique est mise en oeuvre sur des exemples pratiques.

Le sujet de thèse que nous abordons concerne l'identification des lois de comportement élastoplastiques anisotropes en vue de leur utilisation pour la simulation numérique des procédés de mise en forme par déformation plastique de tôles minces métalliques d'emboutissage. Nous avons essentiellement contribué à la définition et à la mise en oeuvre de stratégies et techniques d'identification des lois de comportement à partir d'essais expérimentaux. Les essais classiquement utilisés pour l'identification des modèles : la traction simple dans les axes et hors axes, la traction plane et le gonfelement hydraulique sont présentés et analysés du point de vue homogénéité des déformations pour aboutir à la relation contrainte-déformation à partir des mesures globales force-déplacement. L'identification des modèles de comportement à partir des essais inhomogènes nécessite une simulation numérique couplée avec une méthode d'optimisation (méthode du Simplexe) pour minimiser l'écart entre les résultats expérimentaux et la réponse du modèle calculée par une méthode d'élements finis. Nous avons identifié des lois de comportement élastoplastiques anisotropes avec écrouissage isotrope. Ces lois sont basées en particulier sur le choix d'une ou deux fonctions "contrainte équivalente" définissant le critère de plasticité et le potentiel plastique (cadre de la normalité non associée) ayant la même structure que la fonction seuil. Plusieurs critères quadratiques et non quadratiques sont alors utilisés. Nous avons développé une technique spécifique d'analyse de sensibilité des essais par rapport aux differents paramètres à identifier de la loi de comportement. Cette technique est mise en oeuvre sur des exemples pratiques.

Les modèles continus de plasticité cristalline appliqués aux calculs de polycristaux métalliques sont efficaces pour voir le comportement mécanique global du polycristal, à partir des lois de comportement du monocristal. On obtient ainsi également les champs de contraintes et déformations locaux dans les grains. Ceci est rendu possible à l'aide de simulations par éléments finis sur un volume élémentaire représentatif d'agrégat et les modèles d'homogénéisation numérique. Ces approches échouent néanmoins pour décrire les effets d'échelle, classiquement observés en métallurgie physique et dont l'archétype est l'effet de taille de grain. Un modèle de plasticité cristalline de Cosserat appliqué dans le cas du comportement élastoplastique d'aciers IF ferritiques est proposé dans ce travail. Il introduit dans sa formulation une loi de durcissement supplémentaire associé à la courbure de réseau. La simulation d'agrégats polycristallins permet de reproduire numériquement un effet analogue à la loi de Hall-Petch. Une autre limite de la plasticité cristalline est liée à l'étape clé de validation expérimentale locale. L'information expérimentale est, en général, disponible à la surface de l'éprouvette.<br />La géométrie des grains sous la surface est cependant inconnue. Cette information est le plus souvent introduite dans les calculs par extension des joints de grains perpendiculairement à la surface. L'écart, fréquemment observé entre les résultats du calcul et les résultats expérimentaux, peut être expliqué par l'erreur qu'introduit ce choix. On donne ici un minorant de cette erreur en considérant plusieurs agrégats ayant la même morphologie granulaire à la surface libre mais des morphologies tridimensionnelles distinctes. En élasticité la dispersion des contraintes, en un point donné de la surface, avec différentes morphologies de grains sous-jacents est de l'ordre de 30%. En élastoplasticité la dispersion peut aisément atteindre 50% de la valeur de la contrainte, ce qui amène à considérer avec prudence l'identification d'une loi de comportement à partir des seules mesures de surface.

Durant la phase de construction des puits pétrolier, le ciment est coulé dans l'espace annulaire entre le cuvelage et la roche environnante. La gaine de ciment a pour but de tenir le cuvelage, garantir l'étanchéité des puits pétroliers, ou des réservoirs de stockage de CO2, et de protéger le cuvelage de la corrosion. Au cours de la vie du puits, cette gaine de ciment est soumise tout le long du puits à des sollicitations thermiques et mécaniques qui varient au cours du temps et qui peuvent modifier ses propriétés et altérer son étanchéité. L'objet de cette thèse est d'étudier l'effet de la température et des contraintes sur les propriétés mécaniques de la pâte de ciment en cours de prise et de la pâte de ciment durcie. L'approche est basée sur une étude expérimentale qui combine des essais calorimétriques, des mesures de vitesse des ondes et des essais oedométriques avec le système des cellules STCA (Slurry To Cement Analyzer) sur le ciment en cours de prise ainsi que des essais de compression uniaxiale et triaxiale sur la pâte de ciment durcie. Les résultats expérimentaux ont montré que la température et la pression accélèrent la cinétique d'hydratation et que la température affecte significativement les propriétés élastiques du matériau tandis que la pression ne les influence pas pour une gamme de pression limitée à 20MPa. Une modélisation de la cinétique d'hydratation associée à une technique d'homogénéisation est proposée afin d'interpréter les essais. On a mis aussi en évidence que lorsque la prise se fait sous contraintes mécaniques, des déformations irréversibles peuvent se développer dans la pâte de ciment, ce qui peut conduire à la formation d'un micro-annulaire entre la gaine ciment, le cuvelage et la formation géologique. Le comportement différé de la pâte de ciment durcie a été étudié à partir d'essais de fluage uniaxiaux et d'essais de compression isotrope. Les essais ont mis en évidence que le fluage sous chargement uniaxial est plus important pour un ciment hydraté à une température plus élevée, ce qui est attribué à une microstructure plus hétérogène. Un modèle visco-endommageable permettant de reproduire les phases de fluage primaire et tertiaire a été développé. Les essais de compression isotrope drainés et non-drainé isothermes sous forte contrainte ont montré un comportement différé avec hystérésis lors de cycles décharge-recharge. Ces essais ont été analysés à partir d'un modèle de comportement poro-visco-plastique. Le comportement élastoplastique à court terme a été abordé à l'aide des essais triaxiaux sous chargement déviatorique drainé. Ces essais ont mis en évidence que la température d'essai affecte fortement la surface de charge initiale et les déformations tandis qu'elle ne modifie pas significativement la surface de rupture. Un modèle de plasticité avec une surface de charge fermée et un écrouissage dépendant des déformations plastiques accumulées a été développé pour décrire ces essais. Enfin, une étude préliminaire sur les effets des cycles mécaniques et thermiques a été menée. Des cycles thermiques ne dépassant pas la température d'hydratation ne semblent pas affecter les propriétés mécaniques du matériau. Cependant, une dégradation très rapide avec le nombre de cycles mécaniques a été mise en évidence lorsque la contrainte dépasse 60% de la résistance en compression simple

Les méthodes d'homogénéisation numériques permettent de lever les limitations associées aux approches analytiques ou semi-analytiques pour prédire le comportement des matériaux hétérogènes dont les phases sont décrites par un comportement non linéaire, en particulier pour des morphologies de géométries complexes, ou des comportements dépendant de l'histoire, ou en grandes déformations. Dans ce travail, plusieurs contributions aux méthodes d'homogénéisation numériques pour traiter les problèmes non linéaires sont proposées. Dans une première partie, nous introduisons une méthode permettant de réaliser l'homogénéisation des matériaux dont les phases sont élastoplastiques parfaites. La technique est une extension de la méthode Non Uniform Transformation Field Analysis (NTFA) dans laquelle un algorithme de type prédiction-correction est introduit pour actualiser les variables internes. Dans une seconde partie, une méthode pour réaliser l'homogénéisation des composites non linéaires élastiques en petites et grandes déformations est proposée. Celle-ci se base sur la construction de potentiels associés au comportement effectif sous forme d'une base de données composée de valeurs discrètes interpolées à partir de valeurs discrètes obtenus par calculs éléments finis réalisés sur un Volume Élémentaire Représentatif (VER). Un schéma original d'interpolation est introduit, utilisant la notion de réseaux de neurones artificiels pour limiter le nombre de valeurs pré-calculées et choisir les points aléatoirement dans l'espace du chargement et des paramètres. Il est possible d'introduire des paramètres associés à la microstructure dans le comportement tels que la fraction volumique ou des paramètres de comportement des phases, en vue de possible optimisation de la réponse associée des structures non linéaires / Computational homogenization methods allow circumventing issues associated to analytical or semi-analytical approaches, for predicting the effective behavior of heterogeneous materials whose phase are described by a nonlinear behavior, more precisely when microstructures have complex morphologies or a history-depedendent behavior, or when considering finite strains. In this work, several contributions to computational homogenization methods for modeling nonlinear heterogeneous materials are proposed. In a first part, we introduce a method for homogenizing meterials whose phases are perfectly elastoplastic. The technique extends the Non Uniform Transformation Field Analysis (NTFA) method by développing a return-mapping algorithm to update macroscopic internal variables. In a second part, a method is introduced for the homogenization of nonlinear composites whose phases are nonlinearly elastic, at both small and finite strains. The approach is based on the construction of potentials associated to the effective behavior in the form of a database whose discrete values are computed through Finite Element Computations realized on a Representative Volume Element (RVE). An original interpolation scheme is introduced, which is based on artificial Neural Networks to reduce the number of computed values which can be distributed randomly in the parameter space. It is then possible to introduce parameters associated to the microstructure in the constitutive law, such as volume fraction or constitutive parameters of local phases, for optimization of the related non linear structure response

On étudie une approche de type multimécanimes intermédiaire entre deux classes de modèles: les modèles micromécaniques et les modèles macroscopiques. Le but visé par l'application de ces modèles est d'essayer de reproduire simultanément les phénomènes les plus classiques au sens du comportement des matériaux, mais aussi de faire en sorte que les éléments du comportement généralement moins bien décrits soient également pris en compte. Le couplage entre les variables cinématiques des différents mécanismes permet de décrire: (1) l'effet normal à la vitesse de chargement et l'effet inverse à cette vitesse, (2) l'interaction entre plasticité et fluage, (3) du rochet à haute contrainte moyenne et de l'accommodation ou de l'adaptation quand cette contrainte est moins élevée. Par ailleurs, on s'est intéressé à une méthode d'intégration numérique pour minimiser les temps de calcul dans l'analyse de structures viscoplastique sous chargements cycliques par la méthode des EF. Cette méthode dite de sauts de cycles, est basée sur des pas de calcul qui se chiffrent en nombre de cycles. Cet algorithme a été appliqué à des exemples industriels complexes. La méthode de sauts de cycles permet une diminution du temps de calcul d'un facteur 2 à 10. Enfin on s'est intéressé à l'étude de problèmes d'optimisation. Des exemples ont été traités : l'optimisation de forme et l'optimisation de matériau. La méthodologie qu'on a adoptée consiste à utiliser deux logiciels différents dans leur version standard. Le premier logiciel est le code de calcul par éléments finis 2ebulon, et la deuxième est le code d'identification de paramètres Sidolo. Ce dernier est destiné initialement à identifier des paramètres de modèles physique sur des expériences réelles. Les expériences, dans la procédure d'identification, sont remplacées, dans le cas de l'optimisation, par des objectifs à atteindre selon des critères donnés.

Ce travail consiste en la mise au point d'une méthodologie simplifiée, utilisable industriellement, de réalisation virtuelle de procédés en formage de pièces massives avec la prise en compte de l'endommagement. Ceci permettra de prévoir l'apparition de l'endommagement ductile en cours de formage des pièces. On pourra agir sur les paramètres technologiques pertinents du procédé afin de retarder cet endommagement pour obtenir des pièces saines, ou au contraire favoriser celui-ci afin de simuler des procédés de coupe. Une formulation théorique générale du couplage comportement-endommagement est introduite, prenant en compte une loi d'évolution de l'endommagement ductile. Sur le plan numérique un soin particulier a été apporté à l'intégration locale des équations d'évolution couplées. Deux modèles sont proposés et discutés en détails: le couplage « fort » et le couplage « faible ». L'algorithme d'intégration utilisé est implicite et utilise une méthode classique de prédiction élastique-correction plastique pour le calcul des incréments des variables internes par un schéma de Newton-Raphson. Le calcul de la matrice tangente consistante prend en compte l'influence de l'endommagement ductile. Pour valider cette méthodologie, une large gamme de procédés dont certains sont issus de l'industrie, tels que le découpage de tôle, la compression de cylindres, le forgeage et l'extrusion d'un croisillon, la mise en forme d'un écrou six-pans, le filage,… est présentée, et comparée avec des essais expérimentaux quand cela est possible. La capacité de cette méthodologie à prédire correctement l'initiation et la propagation de l'endommagement en mise en forme est clairement démontrée

Ce travail a pour objet l'analyse limite des structures par la méthode des éléments finis. Lorsqu'une structure atteint sa charge limite, certaines de ses composantes sont dans la phase inélastique de leur comportement, alors que dans les parties les plus critiques, du fait de la localisation des déformations inélastiques, se produit la rupture du matériau. Les effets de localisation sont, dans les matériaux fragiles liés à l'apparition et au développement de macro fissures alors qu'ils sont, dans les matériaux ductiles, gouvernés par les bandes de cisaillement localisées. L'étude de la charge limite est ainsi reliée à la modélisation du comportement inélastique standard du matériau mais également à la modélisation des effets localisés correspondant au comportement adoucissant des matériaux. Le comportement inélastique standard du matériau est, dans ce travail, décrit par des modèles élastoplastiques, élastoviscoplastiques ou élastiques non linéaires. Tous les modèles de comportement sont définis en termes d'efforts généralisés. Un certain nombre d'approches mathématiques et d'algorithmes numériques sont disponibles mais sont bien souvent inefficaces et manquent de précision. Ainsi, nous utilisons une approche développée plus récemment s'appuyant sur une méthode d'éléments finis enrichis de discontinuités. Nous avons développé de nouvelles formulations d'éléments standards prenant en compte des cinématiques et des descriptions des champs de déplacements discontinus complexes. Plusieurs formulations d'éléments finis ont été développées pour l'analyse de différents composants structurels. Nous présentons, dans un premier temps, un élément fini dédié à l'analyse limite des plaques en béton armé. La formulation d'un élément de plaque élastoplastique et élastoviscoplastique écrite en efforts généralisés associée à une procédure commune d'intégration sont présentées ensuite. Un élément de coque non linéaire, faisant intervenir une fonction seuil à deux surfaces incluant à la fois un écrouissage isotrope et un écrouissage cinématique est ensuite présenté. Les deux derniers éléments finis développés dans ce travail sont dédiés à la modélisation de la rupture localisée dans les poutres planes et les solides bidimensionnels. L'élément de poutre d'Euler-Bernouilli est enrichi par une discontinuité en rotation. Une stratégie s'appuyant sur l'analyse préalable, par un modèle raffiné, d'une partie de la structure est proposée afin d'obtenir les paramètres du modèle constitutif de la poutre. Enfin, nous présentons la formulation d'un élément quadrangulaire à discontinuité forte dont la cinématique permet de prendre en compte des sauts de déplacements linéaires dans les deux directions normale et tangentielle le long de la surface de discontinuité. Des résultats numériques montrent que les éléments développés ainsi que les algorithmes associés constituent un outil efficace et robuste d'analyse de la charge limite et de la rupture des structures. Parmi les exemples, nous présentons la simulation de la propagation d'une fissure dans un matériau fragile ainsi que le développement d'une bande de cisaillement dans un matériau ductile. Les codes numériques associés aux formulations présentées dans ce travail ont été générés par l'outil de programmation symbolique et d'optimisation de code AceGen. Les performances des éléments sont présentés à travers un grand nombre d'exemples numériques réalisés à partir du code AceFem.