L’objectif principal de ce travail est la compréhension des mécanismes qui gouvernent la fissuration assistée par l’hydrogène dans les métaux, en s’appuyant sur l’analyse expérimentale de la propagation des fissures en atmosphère hydrogénant et de l’interaction entre hydrogène et défauts cristallins, et sur le développement d’un modèle de zone cohésive influencé par l’hydrogène.Des essais de propagation de fissure de fatigue ont été réalisés sous haute pression d'hydrogène gazeux sur le fer de pureté commerciale Armco. Les résultats montrent une forte influence de la pression, de la fréquence et de la valeur de ΔK sur la modification des modes de rupture, et, par conséquent, sur les vitesses de propagation. Afin d’identifier les paramètres physiques pertinents qui gouvernent les modes de rupture, une étude sur l’interaction entre hydrogène et défauts cristallins développés lors d’une sollicitation cyclique a été réalisée. Nous avons observé une augmentation de l’absorption totale d’hydrogène avec la déformation plastique cumulée, qui peut être attribuée à l’augmentation du piégeage de l'hydrogène par les dislocations générées au cours de la déformation. Ces données seront ensuite introduites dans un modèle pour reproduire la modification de la diffusion de l’hydrogène en pointe de fissure, et son effet sur la plasticité.Par ailleurs, des mesures de la déformation plastique hors plan en pointe de fissure en présence d’hydrogène ont permis de proposer une amélioration d’un modèle de zone cohésive en introduisant un effet de l’hydrogène sur le comportement plastique des éléments de volume. De plus, l'étude des composantes de la loi de diffusion de Krom a montré l'importance du gradient de contrainte hydrostatique sur la diffusion et l'accumulation de l'hydrogène en pointe de fissure. Le modèle prédit une forte dépendance de la propagation de fissures vis-à-vis de la diffusion de l’hydrogène en pointe de fissure, et est capable de simuler la propagation de fissure sous chargement statique, validant ainsi la superposition d’une composante de fissuration cyclique et d’une contribution statique (due à la présence d’hydrogène), et expliquant la transition des vitesses de propagation observée expérimentalement. / The main purpose of this work is to understand the mechanisms that govern hydrogen assisted cracking in metals, based on the experimental analysis of crack propagation data under gaseous hydrogen and the interaction between hydrogen and lattice defects on the one hand, and on the development of a cohesive zone model influenced by hydrogen on the other hand.Fatigue crack propagation tests were performed under high pressure of gaseous hydrogen on the Armco iron. The results show a strong influence of the pressure, the frequency and the ΔK value, on the modification of the failure modes and on the fatigue crack growth rates. In order to identify the physical parameters that govern the changing of the failure modes, a study on the interaction between hydrogen and the crystallographic defects developed during a cyclic loading was performed. We observe an increase in the total absorption of hydrogen with the cumulated plastic deformation, which can be attributed to the increase in the hydrogen trapping by the dislocations generated during the cyclic deformation. These data have to be introduced into a numerical model to reproduce the modification of the hydrogen diffusion at the crack tip, and its effect on plasticity.Moreover, measurements of the out-of-plane plastic deformation at the crack tip in presence of hydrogen have conducted to an improvement of the cohesive zone model by introducing an effect of hydrogen on the plastic behavior of the volume elements. In addition, the study of Krom diffusion law components has shown the importance of the hydrostatic stress gradient on the diffusion and accumulation of hydrogen at the crack tip. The model predicts a strong dependence of the crack propagation with respect to the hydrogen diffusion at the crack tip, and it is able to simulate the propagation under static load, thus validating the cyclic cracking and static cracking superposition, and explaining the transient regime in fatigue crack growth rates experimentally observed.

Les travaux de recherche présentés dans ce mémoire concourent à la modélisation et à la simulation de l'endommagement mécanique des structures, avec pour cibles les installations de production et de transport d'électricité d'EDF. L'objectif affiché est la prédiction de l'évolution de défauts, préexistants ou non, à l'échelle macroscopique de la structure, sous chargements quasi-statiques non cycliques. En particulier, les sollicitations répétées (fatigue) ou dynamiques (crash, séismes) sont exclues de ces travaux, ainsi que d'autres modes de ruine tels que le flambement ou les instabilités plastiques. Compte tenu de la difficulté à prédire l'évolution de défauts dans toute sa généralité, nous avons opté pour une stratégie de modélisation échelonnée, selon la nature et la richesse des informations recherchées. Trois niveaux d'attentes ont été répertoriés. · Prédiction de l'amorçage d'un défaut L'approche en taux de restitution d'énergie de la mécanique de la rupture linéaire a été étendue pour couvrir les situations élastoplastiques, l'apparition de défauts, les propagations instables (dynamiques). La prédiction de la rupture reste conservative car la ruine y est associée à l'amorçage d'un défaut. · Modélisation de la propagation d'un défaut selon un trajet prédéfini L'approche proposée vise à décrire l'initiation et la propagation d'un défaut, le long d'un trajet donné. Elle s'appuie sur des modèles de zones cohésives dans lesquels des forces d'interaction sont introduites entre les lèvres de la fissure. Les simulations numériques sont plus complexes que précédemment car tout le régime de propagation (non linéaire) est dorénavant modélisé. · Détermination du trajet de fissuration Lorsque le trajet de fissuration ne peut plus être anticipé, nous proposons une réponse fondée sur la mécanique continue de l'endommagement. Elle soulève des difficultés physiques et numériques de plusieurs ordres, auxquelles nous avons tenté d'apporter quelques réponses. Il nous semble malgré tout que ce n'est qu'à ce prix qu'on peut décrire l'initiation et la propagation de défauts dans toute sa généralité. Enfin, répondre à cet objectif, avec une finalité industrielle, requiert à la fois de formuler des modélisations appropriées et de développer les outils numériques nécessaires à leur mise en œuvre dans les codes de calcul des ingénieurs (Code_Aster dans notre cas). Leur emploi dans les pratiques industrielles impose en outre trois types d'exigence, - robustesse, au sens de la capacité à produire un résultat dans des situations variées ; - fiabilité, c'est-à-dire que le résultat n'est pas excessivement sensible aux choix de modélisations numériques (discrétisations spatiale et temporelle, entre autres) et donc qu'il est reproductible par des ingénieurs d'étude différents ; - performance, pour que l'utilisation des modèles soit compatible avec les délais des études industrielles. Ces préoccupations ont fait l'objet d'une attention particulière. Elles nécessitent souvent un compromis entre la finesse de description physique et le caractère opérationnel de la formulation.

Les travaux présentés dans ce mémoire s'inscrivent dans l'étude et l'amélioration des modèles d'endommagement continus régularisés (non locaux), l'objectif étant d'étudier la transition entre un champ d'endommagement continu défini sur l'ensemble d'une structure et un modèle discontinu de fissuration macroscopique.La première étape consiste en l'étude semi-analytique d'un problème unidimensionnel (barre en traction) visant à identifier une famille de lois d'interface permettant de basculer d'une solution non homogène obtenue avec un modèle continu à gradient d'endommagement vers un modèle discontinu de fissuration cohésive. Ce passage continu / discontinu est construit de telle sorte que l'équivalence énergétique entre les deux modèles soit assurée, et reste exacte quelque soit le niveau de dégradation atteint par le matériau au moment où cette transition est déclenchée.Cette stratégie est ensuite étendue au cadre 2D (et 3D) éléments finis dans le cas de la propagation de fissures rectilignes (et planes) en mode I. Une approche explicite basée sur un critère de dépassement d'une valeur " critique " de l'endommagement est proposée afin de coupler les modèles continus et discontinus au sein d'un même calcul quasi-statique par éléments finis. Enfin, plusieurs résultats de simulations menées avec cette approche couplée sont présentés.

Si, pour les matériaux composites élaborés, de nombreuses études expérimentales ainsi que des modèles locaux de comportement ont été développés, la maîtrise du comportement des interfaces entre couches de surface ou d’assise de chaussées est actuellement un réel verrou scientifique. La méthode de dimensionnement française actuelle ne prend en compte, aux interfaces, que des conditions conventionnelles de collage ou de glissement parfait. Afin d’appréhender le comportement local de l’interphase\interface, les outils de la photomécanique apparaissent incontournables. La présente thèse propose une modélisation de l’interface rugueuse et endommageable par un modèle de zone cohésive en mode mixte. Ce manuscrit de thèse comporte trois chapitres. D’abord, le chapitre I présente un état de l’art sur les interfaces dans les matériaux et les structures et en particulier dans les couches de chaussées. Ensuite dans le chapitre II, une identification expérimentale des paramètres mécaniques et géométriques du modèle d’interface est proposée à travers des essais de traction et de cisaillement et des mesures de la texture (PMT et projections de franges). Les résultats obtenus (adhésion, rugosité,…) seront les paramètres d’entrée d’un modèle d’endommagement d’interface. Enfin, le chapitre III aborde la modélisation des interfaces entre couches de chaussées sous l’angle des modèles de zones cohésives avec la prise en compte de la rugosité géométrique. A l’issu de cette étude, une loi est proposée permettant de prendre en compte l’effet de la rugosité à une échelle locale dans une interface lisse équivalente à l’échelle globale. / Interface between bituminous layers is an important parameter for the pavement computational design.New pathologies in pavement structure require today rational methods taking into account theinterfaces behavior. Due to these concerns, the current study is based on a damage cohesive zonemodel (CZM) in mixed mode of the rough interfaces. The model was initially proposed by Allix-Ladevèze. This thesis presents a comprehensive interface modeling including delay effect, based ondamage energy release rate. The process of the present study is presented in three chapters. The firstchapter present the stat of art of interfaces. The second devoted to identify the parameters of theinterface model and material properties through advanced optical method as Digital Image Correlation(DIC) and (H-DIC). In the third chapter, a study of the influence of the elastic normal and tangentialstiffness and coupling parameters in the mixed mode on the debonding interfacial energy is presented.An analytical model provides relations between the interfaces stiffness, the coupling parameter of theCZM and the interfacial roughness. Then, a parametric numerical analysis is conducted to study theroughness effect on the interface constitutive law. Results show clearly the roughness influence in thiskind of structures. The damage behaviours predicted by the proposed model for pure mode I, puremode II and for mixed mode with taking into account of roughness are found in good agreement withexperimental results.

Les composants face au plasma constituent la première paroi en contact avec le plasma dans les machines de fusion telles que Tore Supra et ITER. Ces composants doivent supporter des flux importants et donc des températures élevées. Ils sont composés d'un matériau de protection, le composite carbone/carbone, d'un matériau de structure, le cuivre chrome zirconium, et d'un matériau, le cuivre doux, qui sert de couche d'adaptation entre le composite carbone/carbone et le cuivre chrome zirconium. L'utilisation de matériaux différents entraîne l'apparition de fortes contraintes résiduelles de fabrication, dues à la dilatation différentielle entre les matériaux, et compromet le fonctionnement durable des machines de fusion puisque l'endommagement apparu à la fabrication peut se propager. L'objectif de la thèse est donc d'étudier l'endommagement du composite carbone/carbone et de la liaison composite/cuivre sous les sollicitations que peuvent subir les composants face au plasma au cours de leur vie. Les comportements mécaniques du composite carbone/carbone et de la liaison composite/cuivre ont donc été étudiés afin de mettre en place les modèles les mieux adaptés à la description de leur comportement. Avec ces modèles, des calculs thermomécaniques ont été réalisés sur les composants face au plasma avec le code éléments finis Cast3M. La fabrication de ces composants engendre des contraintes élevées qui endommagent le composite carbone/carbone et la liaison composite/cuivre. L'endommagement se propage lors des retours à température ambiante et non lors de l'exposition à un flux thermique. Des géométries alternatives pour les composants face au plasma ont été étudiées pour réduire l'endommagement. La relation entre l'endommagement du composite carbone/carbone et sa conductivité thermique a également été démontrée. / Plasma facing components consitute the first wall in contact with plasma in fusion machines such as Tore Supra and ITER. These components have to sustain high heat flux and consequently elevated temperatures. They are made up of an armour material, the carbon/carbon composite, a heat sink structure material, the copper chromium zirconium, and a material, the OFHC copper, which is used as a compliant layer between the carbon/carbon composite and the copper chromium zirconium. Using different materials leads to the apparition of strong residual stresses during manufacturing, because of the thermal expansion mismatch between the materials, and compromises the lasting operation of fusion machines as damage which appeared during manufacturing may propagate. The objective of this study is to understand the damage mechanisms of the carbon/carbon composite and the composite/copper bond under solicitations that plasma facing components may suffer during their life. The mechanical behaviours of carbon/carbon composite and composite/copper bond were studied in order to define the most suitable models to describe these behaviours. With these models, thermomechanical calculations were performed on plasma facing components with the finite element code Cast3M. The manufacturing of the components induces high stresses which damage the carbon/carbon composite and the composite/copper bond. The damage propagates during the cooling down to room temperature and not under heat flux. Alternative geometries for the plasma facing components were studied to reduce damage. The relation between the damage of the carbon/carbon composite and its thermal conductivity was also demonstrated.

Cette étude fait partie d'un programme de recherche concernant la réparation de structures composites par collage de patchs externes. Les objectifs principaux de ce programme sont d'une part l'identification de l'ensemble des facteurs susceptibles d'influencer les performances à long terme de ce type de réparation, et d'autre part de déterminer dans quelle mesure l'utilisation de tels assemblages peuvent s'avérer une solution optimale. La conception d'un tel système passe obligatoirement par le développement d'un outil de simulation et de prédiction robuste du fait des divers mécanismes d'endommagement pouvant intervenir de fa?on très complexe et de la rupture finale du système résultant d'une propagation des zones endommagées. Cette étude compose d'une et d'autre l'aspect de la modélisation numérique, et l'aspect expérimental. Le dialogue entre les résultats numériques et expérimentaux permet, d'une part de comprendre les mécanismes d'endommagement et l'évolution de ce dernier dans le système réparé, d'autre part de valider le modèle numérique. En particulier, nos efforts ont été concentrés, en utilisant le logiciel LS-dyna, sur l'application des modèles de zone cohésive (MZC). Le comportement au délaminage d'un composite carbone/époxyde et de l'adhésif sont étudié avec les essais en mode I, mode II et mode mixte. Une étude paramétrique de MZC est effectuée. Le modèle de zone cohésive validé est utilisé pour modéliser le comportement en traction des composites réparés par collage de patches externes

Les composites stratifiés de plis unidirectionnels en carbone/époxy sont fortement utilisés pour alléger les structures aéronautiques tout en conservant de bonnes propriétés structurales. Toutefois, les avantages de ce type de matériau ne sont pas encore pleinement exploités de par le manque de confiance accordée aux modèles de prévision de l’endommagement, dont notamment ceux concernant le délaminage. C’est pourquoi l’objectif de cette thèse était de développer une stratégie de modélisation du délaminage adaptée aux structures composites stratifiées. Cette stratégie s’est appuyée sur le développement d’un modèle de zone cohésive prenant en compte les ingrédients nécessaires à la bonne description de l’amorçage et de la propagation de la fissure : (i) un critère d’amorçage avec un renforcement en compression/cisaillement hors-plan, (ii) une loi de propagation décrivant l’évolution de la ténacité en fonction de la mixité de mode et (iii) la prise en compte du couplage inter/intralaminaire. Pour identifier ce nouveau modèle, une procédure d’identification efficace, s’appuyant sur un essai de traction sur plaque rainurée, a été mise en place. Cette procédure d’identification a permis de démontrer que la ténacité semble indépendante (i) de l’orientation des plis adjacents à l’interface et (ii) de l’empilement étudié. De même, pour décrire l’évolution de la ténacité, une nouvelle loi de propagation adaptée au matériau carbone/époxy a été proposée. Pour finir, la stratégie de modélisation, complétée par une stratégie de calcul, a été appliquée sur différents cas structuraux pour mettre en avant ses apports et ses premières limites. / The carbon/epoxy laminated composites of unidirectional plies are strongly used in order to reduce the weight of aeronautical structures while at the same time proposing good structural properties. However, the advantages of this kind of material are not fully exploited due to the lack of confidence in damage models, like ones concerning delamination. Then, the purpose of this work was the development of a strategy to model delamination in laminated composite structures. This strategy was based on the development of a cohesive zone model taking into account the ingredients necessary to the well description of the onset of delamination and the crack growth: (i) an onset criterion with an out-of-plan compression/shearing reinforcement, (ii) a propagation law describing the evolution of the fracture toughness as a function of mixed mode ratio and (iii) the inter/intralaminar coupling. To identify this new model, an efficient identification procedure, basing on a tensile test on notched specimen, has been proposed. This identification procedure has demonstrated that the fracture toughness seems to be independent of (i) the orientation of plies closed to the interface and (ii) the stacking sequence. Furthermore, to describe the evolution of the fracture toughness, a new propagation law adapted to carbon/epoxy material has been proposed. Finally, the strategy to model delamination, completed by a calculation strategy, has been applied on several structural cases to prove its contributions and its first limitations.

Cette thèse est une contribution à l'amélioration de la connaissance des comportements mécaniques des tissus biologiques, en particulier dans le contexte accidentologique. L'objectif de cette étude biomécanique est de faire le lien entre les notions de lésions utilisées en clinique et celles d'endommagement et de rupture utilisées en mécanique. Elle porte sur la modélisation du comportement des ligaments dans l'articulation du genou humain, et s'intéresse plus particulièrement à l'insertion du ligament dans l'os. Les lésions qui touchent cette structure se produisent soit dans le ligament, soit dans une région proche de la zone de transition entre l'os et le ligament. Ce problème est original sur les plans anatomique et mécanique. Sur le plan anatomique, il n'existe pas de description fine (c'est-à-dire à l'échelle microscopique) de l'architecture de cette transition, qui paraît brutale à l'oeil nu et qui est un lieu possible de lésion ligamentaire. Sur le plan mécanique, il s'agit de l'étude d'un tissu biologique, qui comprend de plus la transition entre un tissu dur minéralisé (l'os) et un tissu mou non minéralisé à fibres longues (le ligament). La méthodologie a donc été de comprendre le fonctionnement de l'articulation du genou et les lésions et traumatismes dont elle est le lieu, de connaître la composition et l'organisation des tissus qui constituent cette transition, au moyen d'un étude histologique, de décrire son comportement mécanique dans le fonctionnement du ligament, par le développement d'un protocole expérimental, et enfin de développer un modèle de comportement de l'insertion ligamentaire. L'étude porte sur les ligaments croisé postérieur et latéral externe, ligaments dont les mécanismes lésionnels dans le cas des accidents de la route ont été plus particulièrement étudiés. L'architecture de la zone de transition en microscopies optique et électronique a été décrite et a montré qu'elle est la superposition d'un front de minéralisation et d'un changement de la structure du tissu (présence de cartilage fibreux). Cette transition se fait sur une longueur d'environ 300 micromètres. Les essais expérimentaux de traction ont été réalisés sur la structure insertion ligamentaire - ligament - insertion ligamentaire prélevée sur sujets d'anatomie. La sollicitation a été appliquée soit dans la direction des fibres ligamentaires, soit une direction réaliste du point de vue physiologique. Les protocoles développés permettent de faire des essais cycliques sans endommager les tissus et des essais à rupture sous sollicitations quasi-statiques (1 mm/s et 20 mm/s) et dynamiques (0.5 m/s et 1 m/s). Pour les essais dans la configuration physiologique, les genoux ont été testés en extension complète (station érigée) et à 120 degrés de flexion (position de conduite). Les résultats expérimentaux montrent que, pour les vitesses et amplitudes testées, les ligaments dissipent de l'énergie par frottement interne et que leur comportement est très dépendant de l'angle de flexion du genou. La rupture se produit toujours par arrachement osseux au niveau de l'os cortical en quasi-statique. En dynamique, pour la vitesse testée, il y a toujours décohésion entre fibres dans le ligament, mais la rupture se produit le plus souvent par arrachement de l'os spongieux en profondeur. L'étude du comportement du ligament (influence de l'orientation de l'insertion ligamentaire) et de la micro-structure de l'insertion nous ont menés à choisir un modèle d'interface pour décrire son comportement. Deux modèles de zones cohésives (couplant frottement, adhésion et éventuellement endommagement) ont été développés pour prédire les lésions par arrachement osseux et par décohésion entre les fibres ligamentaires.

Les travaux présentés dans ce mémoire s'inscrivent dans l'étude et l'amélioration des modèles d'endommagement continus régularisés (non locaux), l'objectif étant d'étudier la transition entre un champ d'endommagement continu défini sur l'ensemble d'une structure et un modèle discontinu de fissuration macroscopique.La première étape consiste en l'étude semi-analytique d'un problème unidimensionnel (barre en traction) visant à identifier une famille de lois d'interface permettant de basculer d'une solution non homogène obtenue avec un modèle continu à gradient d'endommagement vers un modèle discontinu de fissuration cohésive. Ce passage continu / discontinu est construit de telle sorte que l'équivalence énergétique entre les deux modèles soit assurée, et reste exacte quelque soit le niveau de dégradation atteint par le matériau au moment où cette transition est déclenchée.Cette stratégie est ensuite étendue au cadre 2D (et 3D) éléments finis dans le cas de la propagation de fissures rectilignes (et planes) en mode I. Une approche explicite basée sur un critère de dépassement d'une valeur « critique » de l'endommagement est proposée afin de coupler les modèles continus et discontinus au sein d'un même calcul quasi-statique par éléments finis. Enfin, plusieurs résultats de simulations menées avec cette approche couplée sont présentés. / The present work deals with the study and the improvement of regularized (non local) damage models. It aims to study the transition from a continuous damage field distributed on a structure to a discontinuous macroscopic failure model.First, an analytical one-dimensional study is carried out (on a bar submitted to tensile loading) in order to identify a set of interface laws that enable to switch from an inhomogeneous solution obtained with a continuous gradient damage model to a cohesive zone model. This continuous / discontinuous transition is constructed so that the energetic equivalence between both models remains ensured whatever the damage level reached when switching.This strategy is then extended to the bi-dimensional (and tri-dimensional) case of rectilinear (and plane) crack propagation under mode I loading conditions, in a finite element framework. An explicit approach based on a critical damage criterion that allows coupling both continuous and discontinuous approaches is then proposed. Finally, results of several simulations led with this coupled approach are presented.

L'objectif de ce travail de thèse est de mettre en place une simulation thermo-mécanique couplée d'une aube revêtue permettant de modéliser l'écaillage de la barrière-thermique qui survient dans les conditions de service de l'aube. La barrière thermique est un revêtement isolant déposé à la surface du substrat monocristallin base Nickel AM1 constitutif de l'aube préalablement recouverte d'une sous-couche. Le mode de dégradation dominant dans ces systèmes est la création de fissures qui résultent de l'accroissement des ondulations hors-plan d'une couche intermédiaire d'oxyde formée en service entre la céramique et la sous-couche. En vue de modéliser ce phénomène d'écaillage, un ensemble d'outils numériques permettant de réaliser un calcul 3D par éléments finis thermo-mécanique couplé de l'aube revêtue est développé au sein du code de calcul par éléments finis Z-set. L'insertion d'éléments de zone cohésive mécanique et thermique au niveau de l'interface barrière-thermique/substrat permet de tenir compte simultanément des changements dans le processus de transert de charge et des variations du flux de chaleur causés par l'amorçage et la propagation d'une fissure interfaciale. L'élément fini d'interface mixte de Lorentz qui repose sur un Lagrangien augmenté, est mis en oeuvre. Afin de tenir compte des propriétés structurelles du revêtement, la modélisation de la barrière thermique est réalisée au moyen d'éléments de coque thermo-mécaniques reposant sur l'approche dite “Continuum Based”. Ces éléments sont développés puis validés dans le cadre de la thèse. La méthode utilisée pour réalier le couplage thermo-mécanique est l'algorithme partitioné CSS (Conventional Serial Staggered) sous-cyclé à pas de couplage fixe dont on montre les limitations dans le cas d'une simulation impliquant la propagation d'une fissure. L'introduction de pas de couplage adaptatifs contrôlés au moyen d'une variable interne du problème mécanique a permis de contourner ces limitations. L'ensemble des briques numériques est validé sur des cas tests de complexité croissante. Des cas d'applications effectués sur des géométries tubulaires à gradient thermique de paroi sont réalisés afin de tester le modèle couplé sur des structures et des chargements proches des conditions de service de l'aube. Enfin, des calculs thermo-mécaniques couplés sur aube revêtue sont présentés. / The purpose of this thesis is to perform a coupled thermomechanical simulation of the failure of thermal barrier coatings for turbine blades under service conditions. The thermal barrier coating is an insulating component applied to the single crystal Nickel-based superalloy AM1 substrate which is covered with a bond coat beforehand. The main degradation mode of those systems is due to the initiation and propagation of cracks caused by the out-of-plane undulation growth of an oxide layer formed in service. A set of numerical tools is implemented into the Finite Element code Z-set in order to perform a 3D thermomechanically coupled simulation of the failure of thermal barrier coatings for turbine blades. Inserting thermomechanical cohesive zone elements at the interface between the coating and the substrate makes it possible to account for the changes in the load transfer and the variations in the heat flux as a consequence of interface degradations. The mixed finite interface element of Lorentz based on an Augmented Lagrangian is used. The thermal barrier coating is modelled by means of thermomechanical shell elements implemented using the Continuum-Based approach to take advantage of the structural properties of the coating layer. Moreover, the partitionned CSS (Conventional Serial Staggered) algorithm used to couple thermal and mechanical problems is assessed. The limitations of sub-cycling with constant coupling time-step are shown through a simulation with crack propagation. The introduction of adaptative time-stepping allows to circumvent that issue. The numerical tools are assessed on test cases with increasing complexity. Numerical simulations on cylindrical tube with a thermal through-thickness gradient are performed with realistic loading sequences. Finally, thermomechanical simulations on turbine blades covered with thermal barrier coating are shown.