Modeling crack propagation under extreme loading in Mindlin-Reissner shells using X-FEM

Actuellement, les méthodes classiques (éléments finis, endommagement, critère de rupture) pour analyser la tenue des structures sous des chargements extrêmes sont très dépendantes de la taille de maille du modèle et nécessitent à la fois un savoir-faire spé- cifique dans le domaine et des études de sensibilité au maillage. De nouvelles approches basées sur la méthode des éléments finis étendus permettent de traiter des propagations de fissure sur des structures de petites tailles et volumiques. Cependant, la propagation sur de grandes longueurs avec des modèles volumiques demande une puissance de calcul importante, souvent inaccessible dans le cadre industriel. Cette thèse a pour but de cou- pler des éléments finis de coque avec la méthode des éléments finis étendue (X-FEM). On peut ainsi diminuer la taille des modèles et gagner en temps de calcul. La fissure peut éga- lement évoluer librement dans le maillage. Après avoir fait le choix d’un élément fini de coque simple et de bonne qualité, l’objectif est de modifier cet élément afin de permettre la description d’une fissure au sein même de celui-ci. Ensuite, l’enjeu est d’adapter les critères de propagation qui existent déjà pour des modèles plans ou volumiques pour les matériaux dits "ductiles" afin de les utiliser dans le cadre d’une modélisation coque. Ces critères sont basés sur l’analyse des champs de contrainte et déformation sur un demi- disque aval à la pointe de fissure. Le calcul de la contrainte équivalente extraite de ces champs servant de seuil pour déclencher ou non la propagation est un point clef de ce travail. Cette étude se place dans le cadre de la plasticité généralisée et fait l’hypothèse d’une fissure initialement traversante dans l’épaisseur de la coque. La phase d’amorçage de la fissure n’est pas prise en compte et le défaut initial est supposé préexistant au sein de la structure. En vue de valider le couplage coque/X-FEM et le critère de propagation, des essais de fissuration sur des structures minces sont réalisés et présentés dans ce document.In shipbuilding industry, classical methods to analyze the behavior of structures under extreme loadings are very dependent on the size of the mesh. Moreover, propagation over long lengths with volumetric models requires huge processing power, often inaccessible within this framework. In order to manage these issues and due to the geometry to be considered, a coupling between shell finite element and the extended finite element method (X-FEM) using an adapted propagation criterion is proposed. The developments are made in the fast explicit dynamic finite element code EUROPLEXUS, CEA Saclay. For shell structures involving significant thickness such as submarines, Mindlin-Reissner theory is needed to enable shear strain. Therefore, locking-free element are used to avoid the numerical issue of shear-locking that appears when the shell becomes too thin. The fracture of Mindlin-Reissner plates based on the X-FEM discrete approximation framework is studied by Dolbow and Belytschko with the MITC4. A four node shell element using the same formulation is here only enriched with a step function along the crack line to take into consideration the discontinuity of the displacement field across the crack. The calculation remains accurate without the asymptotic enrichment functions near the crack-tip, as long as the mesh is refined near the crack tip. The numerical integration issue for elements cut by the crack is solved by a partitioning strategy developed by Elguedj. Since the crack is contained in the shell for which the mid plane's position is entirely known, only one information left is needed to locate it. Therefore, a crack is represented by several line segments on the three-dimensional mesh. Only through thickness cracks are considered so far. As regards to the crack propagation, a local criteria proposed by Haboussa is used based on the calculation of mechanical equivalent quantities in the vicinity of the crack tip. The maximum of the equivalent stress tensor near the crack tip is used to decide if the crack propagates as well as its propagation direction, and the Kaninen equation gives the crack velocity