La méthode isogéométrique (IGA) récemment proposée en tant que méthode numérique générique offre de réelles perspectives dans l’unification des modèles géométriques et computationnel. La méthode isogéométrique est intiment liée à la méthode des éléments finis (FEM) étant donné que la méthode est basée sur le même cadre variationnel. Cette méthode a montré dans de nombreuses circonstances de très bonne qualités numériques notamment avec des maillages grossiers (précision numérique, capacité à supporter de grandes déformations…). Notre objectif final dans ce travail est de fournir un environnement de base, numérique et logiciel, pour la simulation de problèmes à champs et physiques multiples pour des pièces élastomériques de type industriel. Dans ce contexte, les points numériques à développer pour l’IGA sont le traitement de l’incompressibilité et le caractère multi-champs du problème thermique dans la formulation de Galerkin. Ainsi dans ce travail nous proposons en premier, un paradigme objet de l’IGA intégré au sein d’une architecture orientée objet en Java, initialement con?ue pour résoudre des problèmes multi-champs couplés en transformations finies. L’approche proposée s’appuie pleinement sur le contexte variationnel existant dans le code dans le cadre des éléments finis pour réduire les développements pour MEF et IGA (une formulation développée en IGA tourne en MEF et vice versa). Dans un second temps, nous avons étudié le problème de l’incompressibilité pour notamment réduire le verrouillage numérique existant toujours sur l’IGA standard. Par un souci de simplicité, nous adoptons des formulations mixtes à 2 champs (déplacement/pression). Afin d’essayer de satisfaire la condition inf-sup en relachant la contrainte sur le déplacement, nous avons développé deux idées de la littérature (naturelle en NURBS) qui consiste à soit dupliquer une fois les n?uds intérieurs du patch des déplacements ou subdiviser les éléments du patch des déplacements. Nous avons étendu ce type d’éléments aux transformations finies. Enfin, et de manière originale, nous avons adopté la même stratégie pour les problèmes à 2-champs pour la thermomécanique. Différentes simulations à petites et grandes déformations confirment le potentiel de l’approche. Enfin, nous évaluons l’ensemble sur un modèle quasi-incompressible thermo-visco-élastique de type Zener sur des éprouvettes classiques dans un contexte physique complexe. / Recently proposed as a general purpose numerical method, the Isogeometric Analysis (IGA) offers great perspective to bridge the gap between CAD and CAE. The IGA is closely related to the finite element method (FEM) as the method is based on the same variational framework. Moreover, this method has shown in many circumstances to be have a better accuracy than the FEM (large mesh distortions…). Our final aim in this work is to simulate complex multiphysics problems for elastomers industrial parts. As matter of fact, the two main numerical issues in this context is the incompressibility/quasi-incompressibility of the material and the thermochemical coupling in Galerkin formulations. First, we propose, a programming paradigm of the IGA in an existing Java object-oriented hierarchy initially designed for solving multi-fields coupled problems at finite strains. We develop an approach that fully take benefit of the original architecture to reduce developments for both FEM and IGA (one problem developed in FEM can be run in IGA and vice versa). Second, we investigate volumetric locking issues persisting for low order NURBS element observed with standard displacement formulation as finite elements. To cure the problem, we adopt two-fields mixed formulation (displacement/pressure) for the sake of simplicity and target at assessing different discretizations in stability (inf-sup condition). The basic idea is to first to increase the internal knot’s multiplicity or to subdivide the patch for displacements. These ideas that are directly inspired from patches properties, have been found in the literature for the Stokes problem and extended to large strain in solid mechanics. The comparison between the two-fields mixed formulation and a strain projection method is lead at small and large strains. At last, we originally adopt a similar strategy for thermomechanical problem at small and large strains. In the context two-fields formulation, displacement/temperature, the LBB stability condition must be fulfilled to guaranty stability. Thus, we investigate the choices of patches for two-fields formulation displacement/temperature fields for IGA applied to thermoelasticity. Several numerical results for thermomechanical problems at small and finite strains, linear and nonlinear have been presented. At last, an incompressible viscous thermo-hyperelastic model is evaluated in the IGA framework with the proposed approach.

Cette étude a pour objectif principal d’établir une stratégie de modélisation robuste, fiable et performante pour décrire des propagations de fissures d’échelle centimétrique en régime ductile dans des composants industriels. Le modèle d’endommagement de GTN écrit en grandes déformations est utilisé pour modéliser l’endommagement ductile. Ce modèle conduit généralement à une localisation de la déformation, conformément à l’expérience. L’échelle caractéristique de ce phénomène est introduite dans les équations de comportement via l’adoption d’une formulation non locale.Sur le plan numérique, ce modèle non local rend bien compte de la localisation dans une bande d’épaisseur donnée lorsqu’on raffine suffisamment le maillage. Par ailleurs, le problème de verrouillage numérique associé au caractère initialement isochore de la déformation plastique est limité en utilisant une formulation à base d’éléments finis mixtes. Enfin, la distorsion des éléments totalement cassés (i.e. sans rigidité apparente), qui pourrait nuire à la bonne convergence des simulations numériques, est traitée par une régularisation viscoélastique.L’ensemble de ces ingrédients sont appliqués pour simuler la propagation de fissure dans un milieu infini plasticité confinée), de sorte à établir un lien avec les approches globales en J-Δa. L’émoussement, l’amorçage et la (grande) propagation de fissure sont bien prédits. Le modèle est également appliqué à une tuyauterie métallique testée en grandeur réelle dans le cadre du projet européen Atlas+. Après une phase d’identification des paramètres sur éprouvette, les réponses globales et locales d’autres éprouvettes et du tube sont confrontés aux résultats expérimentaux. Ces résultats illustrent le degré de robustesse, de fiabilité et de performance qu’on peut attendre du modèle. / The major goal of this work is to establish a robust, reliable and efficient modeling technique so as to describe ductile tearing over a distance of several centimeters in industrial cases. The GTN damage model expressed in the context of finite strains is chosen to model ductile damage. Generally, the model leads to strain localization in agreement with experimental observations. The characteristic length scale of this phenomenon is introduced into the constitutive equations through the use of a nonlocal formulation.On a numerical ground, the nonlocal model controls the width of the localization band as soon as the mesh is sufficiently refined. Besides, the issue of volumetric-locking associated with plastic incompressibility is handled using a mixed finite element formulation. Finally, the distortion of broken elements (i.e. without any stiffness), which may affect the computational convergence of numerical simulations, is treated using a viscoelastic regularization.The improved GTN model is applied to simulate crack propagation under small-scale yielding conditions, so as to establish a relation with the global (J-Δa) approach. Crack tip blunting, crack initiation and (large) crack propagation are well captured. The model is also applied to a full-scale metallic pipe in the framework of the UE project Atlas+. After a phase of parameter calibration based on the experimental results on some small specimens, the global and local responses of other small specimens and of the full-scale pre-cracked pipe are compared with the experimental results. The results illustrates the robustness, the reliability and the efficiency of the current model.