# Mécanique de la rupture

La mécanique de la rupture tend à définir une propriété du matériau qui peut se traduire par sa résistance à la rupture fragile (fracture) ou ductile. Car si les structures sont calculées pour que les contraintes nominales ne dépassent pas, en règle générale, la limite d'élasticité du matériau et soient donc par voie de conséquence à l'abri de la ruine par rupture de type ductile ; elles ne sont pas systématiquement à l'abri d'une ruine causée par la présence d'une fissure préexistante à la mise en service ou créée en service par fatigue (comme lors de la catastrophe ferroviaire de Meudon) ou par corrosion sous contrainte.

## Modes de propagation d'une fissure

La propagation d'une fissure dans un volume peut se décomposer en trois composantes appelées mode :
* I - Ouverture selon la normale au plan de la fissure ;
* II - Ouverture parallèlement au plan de la fissure et perpendiculaire au front de fissure ;
* III - Ouverture parallèlement au plan de la fissure et parallèlement au front de fissure.
Chaque mode est généralement caractérisé individuellement bien que la propagation réelle d'une fissure soit la combinaison de ces 3 modes.

## Mode de rupture fragile

Le mode de rupture de type fragile peut survenir quand les sollicitations ont lieu dans des circonstances de basses températures, de grandes vitesses de chargement, ou de défauts préexistants ou créés pendant le service.
Les ruptures brutales dont il s'agit peuvent être classées en deux catégories. D'une part, les ruptures fragiles liées à l'absence de ductilité du matériau sollicité sous une certaine température (température de paroi minimum admissible) comme les aciers doux ; d'autre part, les ruptures ductiles sans prévenir, c'est-à-dire à très faible déformation plastique. Ces dernières peuvent affecter des matériaux à haute limite d'élasticité où il n'existe pas de dépendance très nette entre la ténacité et la température, c'est-à-dire où la rupture en charge est liée à la propagation quasi instantanée d'une fissure à partir d'un défaut préexistant.
Les essais classiques de ténacité (essais de résilience par exemple) ne permettent pas de définir une grandeur susceptible de prendre en compte le phénomène de rupture de type fragile dans les calculs.
Ce qui suit explique le concept de résistance ou d'endurance à la fracture, ainsi que les paramètres auxquels il est assujetti comme la température, le taux de déformation, la concentration de contraintes ou encore l'intensité de contrainte.
Des méthodes d'examen et d'essai mises à disposition des concepteurs pour apprécier l'aptitude à la rupture fragile sont présentées, ainsi qu'une approche de la mécanique linéaire élastique de rupture.
Dans la conception de structures en matériaux ductiles, l'aptitude à résister à la charge de manière sécurisée est basée sur l'analyse des contraintes de façon à assurer que la contrainte nominale reste confinée dans le domaine élastique du matériau. Des ruptures qui surviennent à l'intérieur de ce domaine élastique (en deçà de la limite d'élasticité du matériau) sont classées comme ruptures fragiles. Ces ruptures peuvent s'initier à partir de petits défauts de continuité au sein du matériau (ligne d'inclusions intermétalliques par exemple) ou de défauts de type fissure qui n'altèrent pas grandement la distribution de la contrainte nominale et qui ne sont habituellement pas pris en considération dans l'analyse de contrainte.
Pour un taux de chargement très élevé, une discontinuité peut fortement diminuer la bonne ductilité du matériau prédite par un essai de traction satisfaisant réalisé sur éprouvettes usinées (sans défaut de surface et aux contours arrondis) et provoquer la ruine de la structure par rupture fragile. C'est ce mode de rupture qui est à l'origine du tragique accident survenu à la plateforme pentagone Alexander Kielland.
Il devient par conséquent évident qu'une conception saine doit bannir toutes les discontinuités (singularités), ce qui reste difficile en construction soudée et surtout en construction moulée.
Pour beaucoup de structures traditionnelles comme les bateaux, les ponts et les équipements sous pression, l'expérience de la conception, des matériaux et des modes de fabrication ont permis d'établir une corrélation satisfaisante, pour le matériau, entre réussite à l'essai de flexion par choc standard sur éprouvette entaillée et bonne tenue en service.
Dans le domaine aéronautique, l'analyse de fatigue permet souvent de prévoir la durée de vie à l'amorçage des fissures. On utilise donc la phase de propagation de fissure pour réaliser des inspections en prévention de la rupture dont les intervalles sont définis par la mécanique de la rupture.
Une des motivations en faveur de l'application des concepts de la mécanique de la rupture et des essais à réaliser sur joints soudés ou pièces moulées est d'être en mesure de se prémunir, dès la phase de conception, des effets produits par les discontinuités présentes. Il est largement reconnu que les joints soudés et les pièces moulées comportent toujours un certain nombre de discontinuités, ce qui place le concepteur devant un dilemme. Le concepteur souhaitera toujours utiliser des joints soudés ou des pièces moulées exempts de discontinuités, ce qui n'est absolument pas réaliste.
L'approche pratique consiste à reconnaître la présence de ces discontinuités au sein des assemblages et de déterminer une taille critique à partir de laquelle une discontinuité devient préjudiciable, autrement dit à partir de laquelle il faut les rechercher et les éliminer. Mais comment décider à partir de quelle dimension critique une discontinuité a toutes les probabilités de devenir préjudiciable ?
Étant donné que la méthode conventionnelle d'essai de flexion par choc n'est plus appropriée, les méthodes d'essais issues de la mécanique de la rupture (essais de fracture), lorsqu'ils sont applicables, peuvent contribuer à établir une corrélation entre taille critique d'une discontinuité et contrainte de rupture pour un matériau ou une soudure donné et donc permettre une estimation directe de la taille des défauts acceptables pour différentes configurations et conditions de fonctionnement. Cependant, l'évolution permanente des matériaux (matériaux à haute limite d'élasticité), les conceptions de plus en plus complexes et les nouvelles technologies de soudage ou de moulage font que les ingénieries ne peuvent pas avoir le recul nécessaire sur le comportement de ces nouvelles conceptions et les contraintes attendues pas forcément attestées. En conséquence il y aura toujours une grande nécessité, pour le concepteur, à traiter de manière analytique le problème des discontinuités.

### L'idée de Griffith

La mécanique de la rupture a été inventée pendant la Première Guerre mondiale par l'ingénieur aéronautique anglais, A. A. Griffith, pour expliquer la rupture des matériaux fragiles. Le travail de Griffith a été motivé par deux faits contradictoires :
* la contrainte nécessaire pour rompre un verre courant est d'environ 100 MPa ;
* la contrainte théorique nécessaire à la rupture de liaisons atomiques est d'environ 10 000 MPa.
Une théorie était nécessaire pour concilier ces observations contradictoires. En outre, les expérimentations sur les fibres de verre que Griffith lui-même a menées suggèrent que la contrainte de rupture augmente d'autant plus que le diamètre des fibres est petit. Par conséquent il en déduit que le paramètre de résistance uniaxiale à la rupture $R_{r}$, utilisé jusqu'alors pour prédire les modes de défaillance dans le calcul des structures, pourrait ne pas être une valeur indépendante des propriétés du matériau.
Griffith suggère que la faiblesse de la résistance à la rupture observée dans ses expériences, ainsi que la dépendance de l'intensité de cette résistance, étaient due à la présence de défauts microscopiques préexistant dans le matériau courant.
Pour vérifier l'hypothèse de défauts préexistants, Griffith a introduit une discontinuité artificielle dans ses échantillons expérimentaux. La discontinuité artificielle était une forme de fissure débouchante plus importante que les autres discontinuités supposées préexistantes dans l'échantillon.
Les expériences ont montré que le produit de la racine carrée de la longueur de défauts (a) et la contrainte à la rupture (σf) était à peu près constante, ce qui est exprimé par l'équation :
    $(1)\qquad \sigma _{f}{\sqrt {a}}\approx C$

### Critère de Griffith : calcul simplifié

Si l'on suppose la présence d'une fissure de taille $a$ dans un matériau sous traction (contrainte $\sigma$), le calcul de la valeur $\sigma _{f}$ de la contrainte à partir de laquelle cette fissure grandit peut s'estimer simplement, cf. l'article Critère de Griffith.
L'explication de cette relation en termes de théorie de l'élasticité linéaire pose un problème de point singulier. En élasticité linéaire, la théorie prédit que la contrainte (et par conséquent l'effort) à l'extrémité d'une fissure dans un matériau idéalement élastique est infinie. Pour éviter ce problème, Griffith a développé une approche thermodynamique pour expliquer la relation qu'il a observée.
Le développement d'une fissure nécessite la création de deux nouvelles surfaces et donc une augmentation de l'énergie de surface. Griffith a trouvé une expression pour la constante C sur le plan de l'énergie de surface de la fissure en résolvant le problème de l'élasticité d'une fissure finie dans une plaque élastique. En bref, l'approche était la suivante :
* calculer l'énergie potentielle stockée dans un échantillon parfait sous une charge de traction uniaxiale,
* ajuster la contrainte de sorte que la charge appliquée ne déforme pas (plastiquement ?) l'échantillon, puis introduction d'une fissure dans l'éprouvette. La fissure relaxe la contrainte et, par voie de conséquence, relaxe l'énergie élastique aux abords des faces de la fissure. D'autre part, de par son existence, la fissure augmente l'énergie de surface totale de l'échantillon.
* calculer la différence d'énergie libre (énergie de surface - énergie élastique) en fonction de la longueur de la fissure. La rupture se produit lorsque l'énergie libre atteint une valeur maximale pour une longueur critique de fissure, au-delà de laquelle l'énergie libre diminue du fait de la propagation de la fissure (augmentation de l'énergie de surface) c'est-à-dire de l'augmentation de la longueur de la fissure jusqu'à provoquer la rupture. Grâce à cette méthode, Griffith a constaté que
    $(2)\qquad C={\sqrt {\cfrac {2E\gamma }{\pi }}}$
où :
    E est le module de Young du matériau, et
    γ est la densité d'énergie de surface du matériau
Les valeurs E = 62 GPa et γ = 1 J/m$^2$ donnent un excellent modèle pour déterminer la contrainte de rupture prédite par Griffith pour un matériau fragile.

### Apports de G. R. Irwin

L'œuvre de Griffith a été largement ignorée par la communauté des ingénieurs jusqu'au début des années 1950. Les raisons semblent être que, pour les matériaux employés dans la réalisation des structures, le niveau réel d'énergie nécessaire pour causer la rupture est de plusieurs ordres de grandeur supérieur à l'énergie de surface correspondante et que, dans les matériaux de construction il y a toujours des déformations plastiques en fond de fissure, ce qui rend l'hypothèse du milieu élastique linéaire avec contraintes infinies en pointe de la fissure tout à fait irréaliste F. Erdogan (2000).
La théorie de Griffith concorde parfaitement avec les données expérimentales sur des matériaux très fragiles tels que le verre. Pour des matériaux ductiles tels que l'acier, bien que la relation $\sigma _{y}{\sqrt {a}}=C$ soit toujours valable, l'énergie de surface (γ) prédite par la théorie de Griffith est souvent irréaliste. Un groupe de travail dirigé par G. R. Irwin à l'US Naval Research Laboratory (NRL), constitué durant la Seconde Guerre mondiale, a réalisé que la plasticité doit jouer un rôle important dans la rupture des matériaux ductiles.
Dans les matériaux ductiles (et même dans des matériaux qui semblent être fragiles), une zone plastique se développe en front de fissure. L'augmentation de la dimension de la zone plastique est fonction de l'augmentation de la charge jusqu'à ce que la fissure se propage libérant les contraintes en arrière du fond de fissure. Le cycle de chargement/libération de chargement plastique aux abords du front de fissure conduit à la dissipation d'énergie comme le ferait un traitement thermique de relaxation de contrainte. Par conséquent, un terme dissipatif doit être ajouté à la relation de l'équilibre énergétique tel qu'élaboré par Griffith pour les matériaux cassants. En termes physiques, de l'énergie supplémentaire est nécessaire pour que la propagation des fissures se produise dans les matériaux ductiles si on les compare aux matériaux fragiles.
La stratégie d'Irwin était de partitionner l'énergie :
* l'énergie stockée en déformation élastique (effet ressort) qui se libère lors de la propagation d'une fissure. Telle est la force motrice thermodynamique de rupture ;
* l'énergie dissipée qui comprend la dissipation plastique et l'énergie de création de surface (et toutes les autres forces dissipatives qui peuvent être au travail). L'énergie thermodynamique dissipée fournit une résistance à la rupture. L'énergie totale dissipée est donnée par :
    $(3)\qquad G=2\gamma +G_{p}$
où $\gamma$ est l'énergie de surface et $G_{p}$ est la dissipation plastique (ainsi que la dissipation provenant d'autres sources) par unité de surface de la fissure.
La version modifiée du critère énergétique de Griffith peut alors être écrite comme :
    $(4)\qquad \sigma _{f}={\sqrt {\cfrac {E~G}{\pi a}}}.$
* pour un matériau fragile comme le verre par exemple, le terme d'énergie de surface domine et $G\approx 2\gamma =2\,\,J/m^{2}$ ;
* pour un matériau ductile comme l'acier par exemple, le terme de dissipation plastique domine et $G\approx G_{p}=1000\,\,J/m^{2}$ ;
* pour les plastiques polymères proches de la température de transition de phase de vitrification, nous avons une valeur intermédiaire de $G\approx 2-1000\,\,J/m^{2}$.

### Facteur d'intensité de contrainte

Une autre réalisation importante du groupe de travail a été de trouver une méthode de calcul de la quantité d'énergie disponible pour une fracture au niveau de la contrainte asymptotique et les champs de déplacement autour d'un front de fissure dans un solide idéalement élastique, dans le cas d'un chargement en mode I.
    $(5)\qquad \sigma _{xy}\approx \left({\cfrac {K_{1}}{\sqrt {2\pi r}}}\right)~f_{xy}(\theta )$
avec « r » représentant la distance au front de fissure.
    $(6)\qquad \sigma _{x}\approx \left({\cfrac {K_{1}}{\sqrt {2\pi r}}}\right)~\cos {\cfrac {\theta }{2}}\left[1-\sin {\cfrac {\theta }{2}}\sin {\cfrac {3\theta }{2}}\right]$
    $(7)\qquad \sigma _{y}\approx \left({\cfrac {K_{1}}{\sqrt {2\pi r}}}\right)~\cos {\cfrac {\theta }{2}}\left[1+\sin {\cfrac {\theta }{2}}\sin {\cfrac {3\theta }{2}}\right]$
    $(8)\qquad \tau _{xy}\approx \left({\cfrac {K_{1}}{\sqrt {2\pi r}}}\right)~\sin {\cfrac {\theta }{2}}\cos {\cfrac {\theta }{2}}\cos {\cfrac {3\theta }{2}}$
En déformation plane :
    $(9)\qquad \sigma _{z}=\nu (\sigma _{x}+\sigma _{y})$
En contrainte plane
    $(10)\qquad \sigma _{z}=0$
$K_{1}$ est le seul paramètre fournissant une caractéristique du champ de contrainte existant aux abords du fond de fissure. C'est le facteur d'intensité de contrainte (en déformation plane). La valeur de $K_{1}$ peut être calculée en analyse de contrainte en fond de fissure (à la pointe de la fissure). Des expressions de $K_{1}$ ont été déterminées pour un grand nombre de cas de chargement et de configuration de pièces. Toutes les expressions sont de la forme :
    $(11)\qquad K_{1}=\beta \sigma {\sqrt {\pi a}}$ en MPa·√m
avec
* $\beta$ facteur tenant compte de la géométrie de la fissure et de la répartition des contraintes ;
* $\sigma$ contraintes dans le matériau normales au plan de fissure et en l'absence de celle-ci.
Au niveau de chargement produisant la propagation brutale d'une fissure, c'est-à-dire entrainant la ruine de la pièce dans laquelle se trouve cette fissure, correspond une valeur particulière de $K_{1}$. Cette valeur particulière de $K_{1}$ est désignée par le symbole $K_{1c}$ et c'est une propriété du matériau au même titre que la limite d'élasticité. Ces deux propriétés varient avec la température, la vitesse de chargement et la structure métallurgique.
$K_{1c}$ caractérise la résistance du matériau à la propagation brutale de fissures.

### Critère en énergie

Lorsque la variation d'énergie élastique est supérieure à la variation d'énergie superficielle, il y a propagation de la fissure ce qui correspond à une diminution de l'énergie libre du système. Dans une plaque de dimension infinie en matériau idéalement élastique contenant une fissure, celle-ci se propage lorsque :
    $(12)\qquad \sigma =\sigma _{r}={\sqrt {\frac {2E\gamma _{s}}{\pi a}}}$
Avec
    E = Module de Young, et
    $\sigma _{r}=$ contrainte de rupture.
$2\gamma _{s}$ correspond à l'énergie nécessaire pour créer une surface unité de fissuration. Il s'agit en fait d'une énergie critique que l'on note : $G_{1c}$ en $J.m^{-2}$.
On peut donc réécrire (12) comme suit :
    $(13)\qquad \sigma =\sigma _{r}={\sqrt {\frac {EG_{1c}}{\pi a}}}$
Les deux critères $K_{1c}$ et $G_{1c}$ caractérisent la propagation brutale d'une fissure. Ces deux critères sont liés par les relations :
Pour un état de déformations planes, cas des structures épaisses ;
    $(14)\qquad K_{1c}={\sqrt {\frac {EG_{1c}}{1-\nu ^{2}}}}$
avec
    $\nu =$ module de Poisson.
Pour un état de contraintes planes, cas des structures minces ;
    $(15)\qquad K_{1c}={\sqrt {EG_{1c}}}$

## Mécanique de la rupture élasto-plastique
### Intégrale J
#### Définition

L'intégrale curviligne J représente un moyen de calculer le taux de restitution de l'énergie de déformation ou de travail (énergie) par unité de surface de zone rompue au sein d'un matériau. Le concept théorique de l'intégrale J a été développé, de façon indépendante, en 1967 par Cherepanov et en 1968 par James Rice. Ces travaux mettent en évidence que le contour délimitant la zone plastique aux abords du front de fissure (appelé J) est indépendant du profil (contour) de la fissure.
Par la suite, des méthodes expérimentales ont été élaborées pour permettre la mesure des propriétés de rupture critiques à partir d'échantillons à l'échelle du laboratoire pour des matériaux dans lesquels la dimension des prélèvements est insuffisante pour garantir la validité des hypothèses de la mécanique linéaire élastique de la rupture, et d'en déduire une valeur critique de l'énergie de rupture $J_{1c}$.
La quantité $J_{1c}$ définit le point à partir duquel se forme une zone plastique dans le matériau au moment de la propagation et pour un mode de chargement$^,$.
L'intégrale J est équivalente au taux de restitution de l'énergie de déformation d'une fissure dans un solide soumis à une charge constante. Cela est vrai, dans des conditions quasi statiques, tant pour les matériaux linéairement élastiques que pour les échantillons expérimentés à petite échelle en passe de céder en front de fissure.

#### Intégrale-J à deux dimensions

L'intégrale-J à deux dimensions a été définie à l'origine par (voir figure ci contre)
    $(16)\qquad J=\int _{\Gamma }\left(W~\mathrm {d} x_{2}-\mathbf {t} \cdot {\cfrac {\partial \mathbf {u} }{\partial x_{1}}}~\mathrm {d} s\right)$
Avec :
    $W(x_{1},x_{2})$ est la densité d'énergie de déformation,
    $x_{1},x_{2}$ sont les deux directions,
    $\mathbf {t} =\mathbf {n} \cdot {\boldsymbol {\sigma }}$ est le vecteur de traction,
    $\mathbf {n}$ est la normale à la courbe $\Gamma$,
    $\sigma$ est le tenseur de contrainte de Cauchy, et
    $\mathbf {u}$ est le vecteur de déplacement.
L'énergie de déformation est donnée par :
    $(17)\qquad W=\int _{0}^{\epsilon }{\boldsymbol {\sigma }}:\mathrm {d} {\boldsymbol {\epsilon }}~;~~{\boldsymbol {\epsilon }}={\tfrac {1}{2}}\left[{\boldsymbol {\nabla }}\mathbf {u} +({\boldsymbol {\nabla }}\mathbf {u} )^{T}\right]~.$
L'intégrale-J aux abords d'un front de fissure est fréquemment exprimée, sous sa forme générale (et en notation d'Einstein), par :
    $(18)\qquad J_{i}=\lim _{\epsilon \rightarrow 0}\int _{\Gamma _{\epsilon }}\left(Wn_{i}-n_{j}\sigma _{jk}~{\cfrac {\partial u_{k}}{\partial x_{i}}}\right)\mathrm {d} \Gamma$
Avec :
    $J_{i}$ est la composante de l'intégrale-J pour une ouverture de fissure dans la direction $x_{i}$,
    $\epsilon$ est une petite région aux abords du front de fissure.
En utilisant le théorème de Green nous pouvons montrer que cette intégrale prend la valeur zéro quand la courbe $\Gamma$ est fermée, que la région ainsi délimitée ne contient pas de point singulier et forme une surface de genre 0. Si les faces de la fissure ne sont pas sous tension, alors l'intégrale-J est indépendante.
Rice a également démontré que la valeur de l'intégrale-J représente le taux de restitution d'énergie pour la propagation des fissures planes.
L'intégrale-J a été développée pour résoudre des difficultés rencontrées dans le calcul des contraintes aux abords d'une fissure dans un matériau linéairement élastique. Rice a montré qu'en mode de chargement constant et sans atteindre l'adaptation plastique, l'intégrale-J peut aussi être utilisée pour calculer le taux de restitution d'énergie dans un matériau plastique.
Pour montrer l'indépendance de l'intégrale-J, nous devons d'abord montrer qu'elle prend la valeur zéro sur une courbe fermée délimitant une région de genre 0. Prenons simplement l'expression $J_{1}$ donnée par :
    $J_{1}:=\int _{\Gamma }\left(Wn_{1}-n_{j}\sigma _{jk}~{\cfrac {\partial u_{k}}{\partial x_{1}}}\right)\mathrm {d} \Gamma$
Nous pouvons écrire que :
    $J_{1}=\int _{\Gamma }\left(W\delta _{1j}-\sigma _{jk}~{\cfrac {\partial u_{k}}{\partial x_{1}}}\right)\;n_{j}\;\mathrm {d} \Gamma$
D'après le théorème de Green nous avons :
    $\int _{\Gamma }f_{j}~n_{j}~\mathrm {d} \Gamma =\int _{A}{\cfrac {\partial f_{j}}{\partial x_{j}}}~\mathrm {d} A$
En utilisant ce résultat nous pouvons exprimer $J_{1}$ comme :
    ${\begin{aligned}J_{1}&=\int _{A}{\cfrac {\partial }{\partial x_{j}}}\left(W\delta _{1j}-\sigma _{jk}~{\cfrac {\partial u_{k}}{\partial x_{1}}}\right)\mathrm {d} A\\&=\int _{A}\left[{\cfrac {\partial W}{\partial x_{1}}}-{\cfrac {\partial \sigma _{jk}}{\partial x_{j}}}~{\cfrac {\partial u_{k}}{\partial x_{1}}}-\sigma _{jk}~{\cfrac {\partial ^{2}u_{k}}{\partial x_{1}\partial x_{j}}}\right]~\mathrm {d} A\end{aligned}}$
Avec :
    $A$ l'aire contenue dans le chemin $\Gamma$.
Maintenant, s'il n'y aucune force présente, l'état d'équilibre (conservation du moment linéaire) impose que :
    ${\boldsymbol {\nabla }}\cdot {\boldsymbol {\sigma }}=\mathbf {0} \qquad \implies \qquad {\cfrac {\partial \sigma _{jk}}{\partial x_{j}}}=0~.$
Aussi :
    ${\boldsymbol {\epsilon }}={\tfrac {1}{2}}\left[{\boldsymbol {\nabla }}\mathbf {u} +({\boldsymbol {\nabla }}\mathbf {u} )^{T}\right]\qquad \implies \qquad \epsilon _{jk}={\tfrac {1}{2}}\left({\cfrac {\partial u_{k}}{\partial x_{j}}}+{\cfrac {\partial u_{j}}{\partial x_{k}}}\right)~.$
Donc :
    $\sigma _{jk}{\cfrac {\partial \epsilon _{jk}}{\partial x_{1}}}={\tfrac {1}{2}}\left(\sigma _{jk}{\cfrac {\partial ^{2}u_{k}}{\partial x_{1}\partial x_{j}}}+\sigma _{jk}{\cfrac {\partial ^{2}u_{j}}{\partial x_{1}\partial x_{k}}}\right)$
À partir de l'équilibre du moment angulaire nous avons $\sigma _{jk}=\sigma _{kj}$. Par conséquent :
    $\sigma _{jk}{\cfrac {\partial \epsilon _{jk}}{\partial x_{1}}}=\sigma _{jk}{\cfrac {\partial ^{2}u_{j}}{\partial x_{1}\partial x_{k}}}$
L'intégrale-J peut s'écrire :
    $J_{1}=\int _{A}\left[{\cfrac {\partial W}{\partial x_{1}}}-\sigma _{jk}~{\cfrac {\partial \epsilon _{jk}}{\partial x_{1}}}\right]~\mathrm {d} A$
Maintenant, pour un matériau élastique la contrainte peut être dérivée de la fonction de l'énergie de stockée $W$ en utilisant :
    $\sigma _{jk}={\cfrac {\partial W}{\partial \epsilon _{jk}}}$
Donc, en différenciant :
    $\sigma _{jk}~{\cfrac {\partial \epsilon _{jk}}{\partial x_{1}}}={\cfrac {\partial W}{\partial \epsilon _{jk}}}~{\cfrac {\partial \epsilon _{jk}}{\partial x_{1}}}={\cfrac {\partial W}{\partial x_{1}}}$
C'est pourquoi nous obtenons $J_{1}=0$ pour une courbe fermée de genre 0 une région sans point singulier (point où la contrainte serait infinie).
Considérons le contour $\Gamma =\Gamma _{1}+\Gamma ^{+}+\Gamma _{2}+\Gamma ^{-}$. Dans la mesure où cette courbe est fermée et délimite une région de genre 0, l'intégrale-J sur cette courbe prend la valeur zéro (cf paragraphe ci-dessus).
    $J=J_{(1)}+J^{+}-J_{(2)}-J^{-}=0$
On admet que les intégrales antihoraires autour de la pointe de la fissure sont de signe positif. Or, puisque les faces de la fissure sont parallèles à l'axe $x_{2}$, la composante normale à ces faces $n_{1}=0$. Aussi, tant que ces surfaces ne sont pas sous tension ; $t_{k}=0$ et il vient :
    $J^{+}=J^{-}=\int _{\Gamma }\left(Wn_{1}-t_{k}~{\cfrac {\partial u_{k}}{\partial x_{1}}}\right)\mathrm {d} \Gamma =0$
Donc,
    $J_{(1)}=J_{(2)}$
et l'intégrale-J est une courbe indépendante du contour.

#### Integrale-J et rupture ductile

Pour les matériaux isotropes présentant une transition ductile / fragile marquée, l'intégrale-J peut être directement liée au mode de rupture ductile.
Dans le cas d'une déformation plane dans les conditions de chargement correspondant au mode I, la relation est donnée par :
    $(19)\qquad J_{Ic}=G_{Ic}=K_{Ic}^{2}\left({\frac {1-\nu ^{2}}{E}}\right)$
Avec :
    $G_{Ic}$ le facteur de relaxation de l'énergie de déformation critique,
    $K_{Ic}$ le facteur d'intensité de contrainte en rupture ductile sous chargement en mode I,
    $\nu$ le coefficient de Poisson, et
    E le module de Young du matériau.
Dans le cas d'un chargement en mode II, la relation entre l'intégrale-J et la rupture ductile en mode II, ($K_{IIc}$) est donnée par :
    $(20)\qquad J_{IIc}=G_{IIc}=K_{IIc}^{2}\left[{\frac {4(1-\nu ^{2})}{3E}}\right]$
Dans le cas d'un chargement en mode III, la relation est donnée par :
    $(21)\qquad J_{IIIc}=G_{IIIc}=K_{IIIc}^{2}\left({\frac {1+\nu }{E}}\right)$

### Essais COD (CTOD)

L'essai COD (Crack Tip Opening Displacement) sert à mesurer l'aptitude d'une structure sous contrainte à propager ou non une fissure à partir d'un défaut existant.
Les nombreuses transformations métallurgiques qui se produisent lors du soudage sont à l'origine de microstuctures complexes qui peuvent avoir des propriétés mécaniques significativement différentes. De plus, la valeur de la ténacité peut être fortement influencée par la position d'un défaut. Aussi, lorsqu'un défaut est détecté dans un assemblage soudé trois attitudes peuvent être envisagées :
* le conserver en l'état en continuant à exploiter le matériel ;
* le réparer, sachant qu'une réparation peut être à l'origine d'un autre défaut encore plus sévère ;
* remplacer la partie défectueuse.
Tout l'avantage de l'essai COD réside dans la possibilité de déterminer la taille d'une indication maximum admissible n'entraînant pas le ruine de l'assemblage en service et ce, dans un contexte précis.
D'autre part, les essais COD peuvent contribuer à l'obtention de garantie de tenue en service sans nécessité d'effectuer un traitement thermique après soudage (TTAS) souvent imposé par les codes de construction à partir d'une certaine épaisseur soudée pour diminuer le niveau de contraintes induits dans l'ouvrage pendant sa construction. Par exemple, les essais COD pourraient permettre de qualifier l'emploi d'acier thermo-mécanique dans la réalisation d'équipements sous pression brut de soudage quelle que soit l'épaisseur soudée (les aciers thermo-mécaniques voient leurs propriétés mécaniques dégradées après TTAS).
Il existe différentes méthodes d'analyse de nocivité, basées sur des méthodes développées pour des structures homogènes et adaptées au cours de ces dernières années aux assemblages soudés.
Parmi les plus connues on peut citer :
* la méthode R6 modifiée (méthode des deux critères),
* la méthode ETM-MM (engineering treatment model for mismatch welds),
* la méthode EMM (equivalent material method),
* la recommandation PD 6493 du British Standard, intitulée Guidance on methods for assessing the acceptability of flaws in fusion welded structures.
Une fois établie la taille maximum admissible de l'indication, celle-ci demeure une indication pouvant éventuellement faire l'objet d'un suivi et, d'autre part, toute indication de taille supérieure à cette limite sera cataloguée comme un défaut à corriger.

#### Le « mismatching »

Les différences de caractéristiques mécaniques entre le métal de base (MB), la zone fondue (ZF) et la zone affectée thermiquement (ZAT) peuvent invalider les résultats obtenus en COD. Ces différences peuvent être mise en évidence, en première approximation, par des filiations de dureté dans les différentes zones voire complétées par des essais mécaniques.
Le « mismatching » s'exprime par le rapport entre la limite d'élasticité de la ZF et celle du MB. Il se note « M ».
Quand M est supérieur à 1 on parle d'over-matching et dans le cas contraire d'under-matching.
La détermination de la ténacité de la zone fondue est influencée par l'effet de mismatching, la largeur de la ZF par rapport à celle de l'éprouvette, le mode de chargement (flexion ou traction), la position de la pointe de la fissure par rapport au cordon, l'évolution de la plasticité dans le ligament.
Par exemple, dans le cas de l'overmatching la plasticité se développe et reste confinée dans le métal le plus mou. Actuellement les codes ne tiennent pas compte de l'effet de mismatching. Ils considèrent le joint soudé comme un matériau homogène prenant en compte les caractéristiques les plus pénalisantes parmi celles constituant la soudure. Cela conduit en général à des résultats conservatifs.