Nous voulons déterminer les solitons spatiaux dans les bres optiques présentant un effet non-linéaire de type Kerr optique. Pour cela, nous proposons une nouvelle approche numérique basée sur la Méthode des Éléments Finis qui est particulièrement bien adaptée à l'étude de milieux inhomogènes. Un modèle scalaire du champ électrique dans les bres optiques est utilisé pour mettre en ÷uvre et valider notre méthode et comprendre la signication physique des nouvelles solutions dans un cas simple. Le champ est supposé harmonique en temps et suivant la direction d'invariance de la bre mais inhomogène selon sa section transverse. Plusieurs exemples dans les bres à saut d'indice et dans les bres optiques microstructurées (FOMs) à c÷ur plein de dimension transverse nie sont décrits. Pour chaque géométrie, une étude complète est réalisée pour obtenir et prouver numériquement l'existence d'une unique solution non-linéaire auto-cohérente ou soliton spatial de plus haute énergie atteignable avant l'autofocalisation instable. Le soliton spatial dépend du prol transverse ni de la structure, correspond au soliton de Townes dans le cas du milieu homogène mais est différent de celui-ci dans les bres optiques étudiées. Notre nouvelle approche de l'étude numérique de l'effet Kerr dans les bres optiques permet ainsi de généraliser le soliton de Townes. Cependant, le prol de l'indice de réfraction induit par l'effet Kerr optique montre que l'approximation du guidage faible n'est plus respectée et donc que les caractéristiques de la solution non-linéaire sont mal évalu ées dans le modèle scalaire. L'implémentation du cas vectoriel sans approximation est donc réalisée. La même étude que dans le cas scalaire est alors menée en consid érant toutes les composantes du champ électrique. Enn, pour étendre le champ d'application de notre méthode numérique, les solutions non-linéaires sont étudiées dans des FOMs avec un défaut creux au centre. Par analogie avec la physique du solide, les congurations donneur/accepteur sont dénies et les solutions non-linéaires sont analysées selon la taille du défaut. L'étude des effets Kerr positif et négatif à travers ces congurations nous permet de découvrir de nouveaux comportements non-linéaires dans ces structures.

Un certain nombre d'applications électrostatiques, notamment les dépoussiéreurs électrostatiques, mettent en oeuvre le phénomène de décharge couronne dans les gaz à partir d'électrodes ayant une forte courbure. Ces électrodes injectent des charges dans le gaz et la charge d'espace qui en résulte modifie la distribution du champ électrique. Le but de notre travail a été de développer un code numérique pour déterminer les distributions de champ électrique et de densité de charge dans la configuration lame-plan, la création de charge se faisant par effet couronne sur le bord de la lame. <br />Nous utilisons la méthode des éléments finis MEF pour résoudre l'équation de Poisson et la méthode des caractéristiques MC pour résoudre l'équation de conservation de la charge. Les deux équations couplées sont résolues par approximations successives en redéfinissant le maillage structuré à chaque itération pour éliminer la diffusion numérique. La redéfinition du maillage structuré est la partie la plus importante du travail. L'algorithme converge bien pour différentes distributions de la charge définies sur la lame injectrice. Les solutions numériques obtenues pour une loi d'injection imposée sur la lame plate donnent des résultats en très bon accord avec les mesures de densité de courant sur l'électrode plane d'un dispositif lame-plan.

De l'automobile à l'aéronautique en passant par le biomédical, tout domaine d'applications utilise des capteurs dont certains sont à sortie fréquentielle. L'arrivée de nouveaux matériaux piézoélectriques a ouvert une porte vers des applications inaccessibles jusqu'à lors et essentiellement à hautes températures. L'orthophosphate de gallium est l'un de ces nouveaux matériaux piézoélectriques ayant de plus un fort facteur de couplage. Ce mémoire porte sur l'étude et la réalisation de capteurs en orthophosphate de gallium fonctionnant en ondes de volumes. La première partie consiste à étudier la sensibilité et stabilité thermique en élaborant un modèle analytique de poutre vibrant en élongation, flexion et torsion. Au cours de cette partie est démontrée l'existence de coupes compensées en température pour ces différents modes. Cette étude théorique est ensuite complétée par une analyse utilisant la méthode des éléments finis (MEF). La réalisation de résonateurs dans diverses orientations cristallographique est effectuée, dont celles compensées en températures. Les mesures associées à ces structures sont confrontées aux résultats théoriques. Vient ensuite l'application d'un biocapteur fonctionnant en cisaillement d'épaisseur. Un modèle éléments finis est élaboré pour étudier l'influence de certains paramètres sur le comportement du capteur. Une comparaison entre le modèle établi et les mesures expérimentales est également présentée. Ces résultats positifs ouvrent des perspectives d'évolution décrites en conclusion.

Ce travail s’inscrit dans le projet Fusible Ecologiquement et Economiquement Efficace destiné à la protection des onduleurs à commutation en tension (FE2E) porté par la société Mersen. Les fusibles étudiés sont constitués de lames d’argent noyées dans du sable de silice aggloméré. Lors d’un court-circuit, un arc électrique apparaît au niveau de rétrécissements dans les lames d’argent (« sections réduites »). Les coupures ultrarapides (moins de 100 microsecondes) ont été très peu étudiées. La structure formée après coupure (fulgurite) est constituée d’un canal d’arc entouré de silice fondue. C’est sur ces structures que cette étude a été réalisée.Après une étude bibliographique, un travail expérimental a été mené pour comprendre comment l’énergie de l’arc est absorbée par la matière. Des observations de coupures par radiographie in situ réalisées à l’ESRF ont permis d’établir des liens entre les phénomènes de création de fulgurite et les mesures électriques. Des observations et des mesures ont été réalisées. Elles ont conduit à établir des corrélations entre le volume, la nature et la structure de la masse isolante générée lors de la coupure et les caractéristiques électriques de l’arc. Les mesures ont en particulier permis d’obtenir des dimensions et des formes caractéristiques des zones affectées thermiquement lors des coupures aux grands di/dt.Grâce à ces informations, des modèles numériques ont été développés, avec pour objectifs de décrire la formation du canal d’arc et d’estimer l’effet du confinement de l’arc par le sable de silice. Des géométries 1D axisymétriques puis 2D ont été utilisées. L’arc électrique était considéré comme une densité de puissance et le sable comme un milieu continu équivalent. Deux cas extrêmes ont été testés, l’un dans lequel la matière vaporisée est confinée dans un canal d’arc fermé, l’autre dans lequel elle s’échappe instantanément. Cela a permis d’identifier les conséquences du confinement sur les phénomènes lors de la coupure.La corrélation entre modélisation et expérience a conduit à un modèle simple dont les résultats sont conformes aux mesures expérimentales et qui permettent de donner des éléments pour le dimensionnement des nouveaux fusibles. Les phénomènes à prendre en compte dans le futur ont été identifiés, comme la propagation de la silice liquide et l’écoulement du gaz dans le milieu granulaire et la formation de l’arc. / This work is part of a project called Environmentally and Economically Effective Fuses (FE2E) driven by the Mersen company. The studied fuses are silver blades embedded in agglomerated silica sand. In case of short circuit, an electric arc appears in shrinked parts (“notches”) of silver blades. The structure formed after the cut (fulgurite) consists in an arc channel surrounded by amorphous silica. Extremely fast cuts (less than 100 microseconds) have been poorly studied.After a bibliographical review, an experimental work was performed to understand how the arc power is absorbed by the surrounding material. Power cuts observations using in situ radiography performed at the ESRF allowed to establish links between fulgurite formation phenomena and electric measurements. Observations and measurements led to correlations between volume, nature and structure of the insulating mass generated during the power cut and electrical characteristics of the arc. In particular, these measurements provided characteristic sizes and shapes for heat affected zones during power cuts at high di/dtNumerical models were developed on this basis for describing the arc channel formation and determining the role of silica sand containment. 1D then 2D axisymmetric geometries were used. The electric arc was considered as a power source and the sand as an equivalent continuous medium. Two extreme cases were tested: when the vaporized material was confined to a closed arc channel and when it escapeed immediately from the arc channel. The results helped identify the consequences of gas containment in the arc channel during power cut.Modeling vs experience correlations led to a simple model which gives information for optimizing new fuses. Some phenomena to be taken into account in future works have been identified, such as the liquid silica penetration and gas flow in the granular packing and the early stages of arc formation.

Cette contribution, en s’appuyant sur expérimentation et modélisation numérique vise à une meilleure compréhension du comportement des dalles en béton armé sous sollicitations de cisaillement. Une campagne expérimentale a été réalisée sur des dalles épaisses à pleine échelle de centrales nucléaires. Ces dalles sans armatures d’effort tranchant sont soumises à une sollicitation de cisaillement en chargement quasi-statique. Les essais sont réalisés en faisant varier différents paramètres qui peuvent influencer le comportement au cisaillement. Sont ainsi étudiés : résistance en compression du béton, épaisseur, taux d’armatures longitudinales et transversales, taille des granulats, longueur de la plaque de chargement. L’influence des efforts de membrane, de compression ou de traction, sur le comportement au cisaillement a également été analysée. Les résultats des essais sont ensuite comparés aux prédictions des codes de calcul. Ces résultats ont d’abord permis d’apporter une réponse aux divergences qui existent entre l’Eurocode 2 et l’Annexe Nationale Française quant à la prédiction du cisaillement. Ont également été évalués le niveau de précision donné par d’autres normes de dimensionnement au cisaillement: la norme américaine ACI 318-14, le code nucléaire AFCEN ETC-C 2010, le fib-Model Code 2010 et l’approche par la théorie de la fissure critique de cisaillement CSCT. Ensuite est évalué la possibilité d’analyses non-linéaire par élément finis (EF) pour reproduire le phénomène du cisaillement dans les dalles. Un modèle de béton élastoplastique avec endommagement est combiné à une analyse quasi-statique à schéma de résolution explicite. Des lois de comportement non linéaires appropriées du béton avec des comportements post-pic associés à un critère énergétique ont été considérées. La bonne concordance entre le modèle proposé et les résultats expérimentaux en termes de résistance au cisaillement et de modes de rupture permet de valider la modélisation proposée. Une étude paramétrique a été réalisée sur la base du modèle proposé avec les mêmes propriétés mécaniques de béton. Des lois simplifiées permettant d’estimer les capacités en cisaillement en fonction des différents paramètres étudiés sont finalement proposées. / This study, based on experiments and numerical modeling, aims at a better understanding of the shear behavior of reinforced concrete slabs. An experimental campaign was carried out on full-scale thick slabs typical of nuclear power plant slabs. These slabs without shear reinforcement are subjected to a quasi-static shear loading. The tests are carried out by varying different parameters that can influence the shear behavior: the concrete compressive strength, the slab depth, the bottom longitudinal and transverse reinforcement ratio, the concrete aggregate size, the loading plate length. The influence on shear behavior of compression or tension membrane forces has also been analyzed. The results of tests are then compared with the predictions of the calculation codes. These results first of all helped to answer the differences between the Eurocode 2 and the French National Annex concerning the prediction of the shear capacity of reinforced concrete slabs. The level of accuracy given by other shear dimensioning standards was also assessed: The American standard ACI 318-14, the AFCEN ETC-C 2010 code used for nuclear buildings, the fib-Model 2010 and the Critical Shear Crack Theory. Next, we evaluate the possibilities of a non-linear finite element analysis (EF) to reproduce the phenomenon of shear in slabs. An elastoplastic concrete model with damage was used and combined with a quasi-static analysis using an explicit resolution scheme. Appropriate nonlinear behavior laws of concrete with post-peak behaviors associated with an energy criterion were considered. The good agreement between the proposed model and the experimental results in terms of shear strength and failure modes allowed validating the proposed modeling. A parametric study was conducted based on the numerical proposed model with the same mechanical properties of concrete. Simplified laws allowing estimating the shear capacities according to the different parameters studied are proposed.

Ce mémoire de thèse à été réalisée dans le cadre d'une collaboration scientifique avec "La Manufacture Française des Pneumatiques Michelin". Il porte sur l'analyse mathématique et numérique de la convergence et de la stabilité de formulations mixtes ou hybrides de problèmes d'optimisation sous contrainte avec la méthode des multiplicateurs de Lagrange et dans le cadre de la méthode éléments finis étendus (XFEM). Tout d'abord, nous essayons de démontrer la stabilité de la discrétisation X-FEM pour le problème d'élasticité linéaire incompressible en statique. Le deuxième axe, qui représente le contenu principal de la thèse est dédié à l'étude de certaines méthodes de multiplicateur de Lagrange stabilisées. La particularité de ces méthodes est que la stabilité du multiplicateur est assurée par l'ajout de termes supplémentaires dans la formulation faible. Dans ce contexte, nous commençons par l'étude de la méthode de stabilisation de Barbosa-Hughes appliquée au problème de contact unilatéral sans frottement avec XFEM cut-off. Ensuite, nous construisons une nouvelle méthode basée sur des techniques de projections locales pour stabiliser un problème de Dirichlet dans le cadre de X-FEM et une approche de type domaine fictif. Nous faisons aussi une étude comparative entre la stabilisation avec la technique de projection locale et la stabilisation de Barbosa-Hughes. Enfin, nous appliquons cette nouvelle méthode de stabilisation aux problèmes de contact unilatéral en élastostatique avec frottement de Tresca dans le cadre de X-FEM. / This Ph.D. thesis was done in collaboration with "La Manufacture Française des Pneumatiques Michelin". It concerns the mathematical and numerical analysis of convergence and stability of mixed or hybrid formulation of constrained optimization problem with Lagrange multiplier method in the framework of the eXtended Finite Element Method (XFEM). First we try to prove the stability of the X-FEM discretization for incompressible elastostatic problem by ensured a LBB condition. The second axis, which present the main content of the thesis, is dedicated to the use of some stabilized Lagrange multiplier methods. The particularity of these stabilized methods is that the stability of the multiplier is provided by adding supplementary terms in the weak formulation. In this context, we study the Barbosa-Hughes stabilization technique applied to the frictionless unilateral contact problem with XFEM-cut-off. Then we present a new consistent method based on local projections for the stabilization of a Dirichlet condition in the framework of extended finite element method with a fictitious domain approach. Moreover we make comparative study between the local projection stabilization and the Barbosa-Hughes stabilization. Finally we use the local projection stabilization to approximate the two-dimensional linear elastostatics unilateral contact problem with Tresca frictional in the framework of the eXtended Finite Element Method X-FEM.

The acoustics of small cavities raises interest of the scientific community since it involves particular damping mechanisms. In fluid dynamics, when a small perturbation is propagating within a Newtonian and heat-conducting fluid bounded by a rigid and isothermal surface, viscous and thermal dissipative mechanisms are generated near the walls. Such effects can have significant impact on the acoustic behaviour of the system.<p>Several types of practical applications can be cited, among which: hearing aids, micro-electro-mechanical systems (transducers, microphones and loud-speakers), absorbing materials made of thin capillary nets or small pores, dissipative silencers, thermo-acoustic heat exchangers, or any kind of device bringing into play small resonant cavities filled with a dissipative fluid (micro-acoustics).<p><p>This study focuses on appropriated reductions of the physical equations, in order to enhance the efficiency of the numerical resolution without adversely affecting the accuracy. Moreover, the proposed strategies lead to numerically stable systems as they involve only one scalar partial order differential equation (or equivalent fluid equation). The emphasis is put on the physical aspect of those reductions, their range of applicability, benefits and drawbacks.<p>Two new reduced models are proposed to estimate the visco-thermal acoustic wave propagation. A first extension deals with waveguide geometries and relax the hypothesis of the fluid at rest. The second original formulation addresses visco-thermal acoustics in 3D arbitrary geometries. This model is based on different considerations coming from existing techniques as well as the estimation of a wall-distance field.<p><p>A second part aims at studying the acoustic behaviour of biphasic materials and more specifically poro-elastic materials. This type of acoustic component is widely used in industry because of their good absorbing properties in the medium- and high-frequency <p>ranges.<p>A preliminary bibliographic research deals with the derivation of the set of partial order differential equations that account for both fluid/structure interactions and the anisotropy of a given poro-elastic material. It has been shown that transversely orientated capillary materials (for instance catalyst substrates) can be simulated using the proposed reduction technique.<p>At last, the modelling of the acoustic transmission between two domains separated by perforated or micro-perforated plates or thin plates of poro-elastic materials is discussed. The analogy between the rigid perforated plate models with an equivalent fluid formulation has been presented. As a result, this model has been extended in order to account for flexural effects of the solid part.<p><p><p>Ce travail porte sur l'étude de certains phénomènes d'amortissements intervenant dans l'acoustique des petites cavités. En méchanique des fluides, lorsqu'une petite perturbation se propage au sein d'un fluide newtonien et caloporteur borné par un mur rigide et isotherme, ces mécanismes dissipatifs particuliers se localisent aux abords des parois et jouent un rôle significatif dans certaines situations.<p>Parmi les exemples d'applications pratiques, il est possible de citer les appareils d'aide auditive, les systèmes microélectromécaniques (transducteurs, microphones et haut-parleurs), les matériaux absorbants constitués de fins réseaux capillaires ou de pores aux dimensions réduites, les systèmes de silencieux, d'échangeurs de chaleur thermo-acoustiques ou tout autre appareil mettant en jeu des cavités résonantes aux dimensions réduites (micro-acoustique).<p><p>L'étude proposée ici se focalise sur des stratégies de réduction appropriées des équations physiques, ceci afin d'améliorer l'efficacité du modèle tout en conservant une précision acceptable. Les techniques présentées aboutissent à des systèmes numériquement stables mettant en jeu une seule équation scalaire (ou équation fluide équivalent). Ainsi, l'accent est porté sur l'aspect physique des réductions, leurs domaines d'application, avantages et inconvénients.<p>Deux modèles originaux sont proposés afin de prédire la propagation acoustique visco-thermique. Une première extension permet d'évaluer la pression acoustique au sein de géométries particulières de type guides d'onde en présence d'un écoulement hydrodynamique. La seconde formulation présentée s'intéresse à l'acoustique dans des domaines 3D arbitraires. Cette méthode se base sur des considérations conjointes de modèles réduits existants ainsi que sur l'estimation d'un champ de distance à la plus proche paroi.<p><p>Dans une seconde partie, nous nous proposons d'étudier le comportement acoustique de matériaux biphasique et plus précisément les matériaux poro-élastiques (très utilisés dans l'industrie en raison de leurs caractéristiques absorbantes dans les domaines des moyennes et hautes fréquences).<p>Une étude bibliographique préliminaire nous a permis d'exprimer l'ensemble des équations aux dérivées <p>partielles modélisant à la fois les interactions fluide/structure et l'anisotropie générale des matériaux. <p>Cette réflexion nous a permis d'aboutir à un modèle de matériau isotrope transverse intéressant, combinant le modèle fluide proposé et la formulation acousto-élastique équivalente. Ainsi la modélisation de structures capillaires orientées (comme les matériaux utilisés dans les catalyseurs automobiles) s'en trouve grandement simplifiée.<p>Enfin, la transmission acoustique intervenant entre deux domaines fluides séparés par une plaque perforée ou micro-perforée ou bien une couche de matériau poreux a été étudiée. L'analogie entre les modèlisations existantes et un modèle générique de fluide équivalent a été mise en évidence. Pour finir, cette formulation a été étendue afin de prendre en compte les effets de flexion de la partie solide. / Doctorat en Sciences de l'ingénieur / info:eu-repo/semantics/nonPublished

Cette thèse porte sur le thème général de la fiabilité des circuits microélectroniques. Le but de notre travail fut de développer un outil de simulation multi-physiques pour la conception des circuits intégrés fiables qui possède les caractéristiques innovatrices suivantes : • (i) L’intégration dans un environnement de conception microélectronique standard, tel que l’environnement Cadence® ; • (ii) La possibilité de simulation, sur de longues durées, du comportement des circuits CMOS analogiques en tenant compte du phénomène de vieillissement ; • (iii) La simulation de plusieurs physiques (électrique-thermique-mécanique) couplées dans ce même environnement de CAO en utilisant la méthode de simulation directe. Ce travail de thèse a été réalisé en passant par trois grandes étapes traduites par les trois parties de ce manuscrit. / This thesis was carried out under the theme of the microelectronics Integrated Circuits Reliability. The aim of our work was to develop a multi-physics simulation tool for the design of reliable integrated circuits. This tool has the following innovative features : • (i) The integration in a standard microelectronics design environment, such as the Cadence® environment ;• (ii) The possibility of efficient simulation, over long periods, of analog CMOS circuits taking into account the aging henomenon ; • (iii) The simulation of multiple physical behaviours of ICs (electrical-thermalmechanical) coupled in the same environment using the direct simulation method. This work was carried out through three main stages detailed in the three parts of this Manuscript.

Les cordes textiles torsadées à fibres continues sont couramment utilisées en tant que renforts dans une grande variété d'applications industrielles et plus spécifiquement comme renforts dans le caoutchouc des pneumatiques. Les cordes concernées sont fabriquées à partir de fibres de nylon, de polyester ou d’aramide en les torsadant ensemble (de 200 à 1000 fibres de même matériau) pour obtenir ce qu’on appelle un surtors. L’assemblage de deux ou trois surtors en les torsadant ensemble permet de construire un retors qui sera dit hybride si les surtors initiaux sont composés de fibres différentes. L’ensemble des configurations donne une large gamme de propriétés mécaniques possibles impliquant à la fois les effets propres au comportement du matériau, à l’orientation locale des fibres ainsi qu’aux pressions transverses induisant des déséquilibres dans le cas de la présence de différentes natures de fibres. Cette thèse présente une étude structurelle précise des orientations de fibres par microtomographie aux rayons X dans les surtors et les retors ainsi qu’un outil original permettant d’analyser les fréquences spatiales caractéristiques des trajectoires de fibre. Un modèle mécanique de surtors en extension de la littérature est enrichi à partir de ces travaux appuyés d’une étude expérimentale supplémentaire de leur contraction latérale sous traction. L’influence du désordre du paquet de fibres initial sur le comportement mécanique du surtors final est analysée numériquement par la méthode des éléments finis en simulant le processus physique ayant lieu au cours de la prise de torsion. Enfin, un modèle mécanique général est proposé pour prédire le comportement mécanique en extension des retors hybrides constitué de trois surtors pouvant chacun afficher des torsions différentes entre elles. / Twisted continuous filament yarns are commonly used as reinforcements in a wide variety of industrial applications and more specifically as reinforcements in tyre rubber. The yarns concerned are made of nylon, polyester or aramid fibres by twisting them together (from 200 to 1000 fibres of the same material) to obtain what is called a single yarn. The assembly of two or three single yarns by twisting them together makes it possible to build a multiply yarn that will be called hybrid if the initial single yarns are made of different fibres. All configurations give a wide range of possible mechanical properties involving the effects of the material's behaviour, the local orientation of the fibres and the transverse pressures inducing an imbalance in the case of the presence of different types of fibres. This thesis presents a precise structural study of fiber orientations by X-ray microtomography of single and multiply yarns and an original tool to analyze the spatial frequencies characteristic of fiber trajectories. A mechanical model of single yarn in extension from the literature is enriched from this work supported by an additional experimental study of their lateral contraction under extension. The influence of the disorder of the initial fibre bundle on the mechanical behaviour of the final single yarns is analysed numerically by the finite element method by simulating the physical process taking place during twisting. Finally, a general mechanical model is proposed to predict the mechanical behaviour in extension of the hybrid multiply yarns, consisting of three single yarns showing different torsions.

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.