A Projection Method for an Elasto-plasticity Model with Linear Kinematic Hardening

1. Problème Étudié

L'article traite d'un système dynamique d'élasto-plasticité pour des matériaux présentant un écrouissage cinématique linéaire de type Melan–Prager, couplé à une viscosité de type Kelvin–Voigt.

Le modèle physique décrit le comportement d'un corps occupant un domaine borné $\Omega \subset \mathbb{R}^d$. Les inconnues principales sont :

• La vitesse $v = \partial u / \partial t$ (où $u$ est le déplacement).
• Le tenseur des contraintes $\sigma$.
• La variable interne de backstress $\alpha$.

Le système est régi par les équations suivantes :

1. Équation du mouvement (avec viscosité) : $\rho \partial_t v = \eta \text{div} E(v) + \text{div} \sigma + F$.
2. Loi de comportement élasto-plastique : Une inclusion variationnelle liant la dérivée temporelle de la contrainte à la déformation et à la règle d'écoulement plastique.
3. Évolution du backstress : $\partial_t \alpha = a \xi$, où $\xi$ est le taux de déformation plastique et $a > 0$ est un paramètre d'écrouissage.

Caractéristique clé du modèle :
La contrainte admissible n'est pas fixe. L'ensemble de contraintes $K(t; \alpha)$ est défini comme un ensemble convexe fermé $\tilde{K}(t)$ (défini par un critère de von Mises avec une borne d'écoulement $g(t, x)$) translaté par la variable interne $\alpha(t)$.
$$K(t; \alpha) = \tilde{K}(t) + \alpha(t)$$
Cette dépendance de l'ensemble de contraintes par rapport à une inconnue ($\alpha$) confère au système une structure quasi-variationnelle. De plus, la borne d'écoulement $g$ peut dépendre du temps et de l'espace sans hypothèse de régularité spatiale stricte ni de borne inférieure strictement positive.

2. Méthodologie

Pour construire des solutions, les auteurs utilisent la méthode de Rothe (discrétisation temporelle) combinée à une méthode de projection.

Schéma de discrétisation proposé (P)

Le temps est discrétisé en pas $\Delta t = T/N$. À chaque pas de temps $n$, le calcul se décompose en trois étapes :

1. Étape de calcul élasto-visqueux (Essai) :
   On résout un sous-problème linéaire pour obtenir une vitesse $v_n$ et une contrainte d'essai $\sigma^*_n$. Cette étape ignore temporairement la contrainte de plasticité, agissant comme une prédiction élastique-visqueuse.
   $$S \frac{\sigma^*_n - \sigma_{n-1}}{\Delta t} = E(v_n) + h_n$$
   (où $h_n$ regroupe les termes de force et de conditions aux limites).

2. Mise à jour du backstress :
   Le backstress $\alpha_n$ est mis à jour en utilisant une approximation du taux de déformation plastique $\xi_n$.

3. Étape de Projection :
   La contrainte d'essai, une fois corrigée par le backstress, est projetée sur l'ensemble de contraintes translaté. L'équation clé est :
   $$\sigma_n - \alpha_n = P_{g_n}(\sigma^*_n - \alpha_n)$$
   où $P_{g_n}$ est l'opérateur de projection ponctuelle sur l'ensemble défini par la borne $g_n$. Cette projection garantit que la contrainte finale respecte le critère de plasticité.

Résolution du couplage :
Bien que le système semble implicite et non linéaire à chaque pas (car $\sigma_n$ et $\alpha_n$ sont couplés), les auteurs montrent qu'il peut être résolu explicitement point par point (presque partout dans $\Omega$) en exploitant la géométrie de la projection et la relation linéaire entre les parties déviatoriques des tenseurs.

3. Résultats Principaux

Les auteurs établissent les résultats théoriques suivants pour le problème continu (Problème 1.1) :

• Stabilité du schéma discret :
  Ils démontrent des estimations d'énergie uniformes (indépendantes de $\Delta t$) pour les solutions discrètes dans des espaces de Sobolev appropriés ($L^\infty(0,T; H)$, $L^2(0,T; V)$, $H^1(0,T; H)$, etc.). Ces estimations reposent sur l'inégalité de Korn, la positivité définie du tenseur de compliance $S$, et une inégalité de Gronwall discrète.

• Existence de solution faible :
  En utilisant la compacité (théorème d'Aubin-Lions et théorème d'Ascoli-Arzelà) et le passage à la limite lorsque $\Delta t \to 0$, ils prouvent l'existence d'une solution faible au sens variationnel.

  • La limite des solutions discrètes satisfait l'équation du mouvement et l'inclusion variationnelle.
  • La contrainte $\sigma(t) \in K(t; \alpha(t))$ est préservée à la limite grâce à la fermeture de l'ensemble convexe et la convergence forte de la borne $g$.

• Unicité de la solution :
  L'unicité est prouvée par un argument d'énergie. En considérant la différence entre deux solutions hypothétiques et en testant les équations avec cette différence, on obtient une inégalité différentielle qui force la différence à être nulle, grâce à la monotonie des opérateurs et la positivité des paramètres physiques ($\eta, a > 0$).

4. Contributions Clés

1. Faiblissement des hypothèses sur la borne d'écoulement $g$ :
   Contrairement à des travaux antérieurs (comme [2]), ce papier ne nécessite pas que $g$ soit spatialement continu ni qu'il soit strictement borné inférieurement par une constante positive ($0 < c \le g$). Cela rend le modèle applicable à des situations où la limite d'écoulement peut s'annuler localement ou présenter des discontinuités.

2. Approche par projection pour l'écrouissage cinématique :
   L'article adapte une méthode de projection (habituellement utilisée pour des ensembles de contraintes fixes ou prescrits) à un cadre quasi-variationnel où l'ensemble dépend de la solution elle-même. Cela permet de traiter le couplage complexe entre le mouvement et l'écrouissage de manière constructive.

3. Base constructive pour l'approximation numérique :
   La structure du schéma (P) est directement compatible avec une discrétisation spatiale par éléments finis. Le schéma proposé fournit une base algorithmique directe pour la simulation numérique, évitant la nécessité de résoudre des problèmes d'optimisation non linéaires complexes à chaque itération temporelle grâce à la forme explicite de la projection.

5. Signification et Perspectives

Ce travail est significatif car il fournit une preuve rigoureuse d'existence et d'unicité pour un modèle dynamique d'élasto-plasticité avec écrouissage cinématique, dans un cadre très général concernant la borne d'écoulement.

Perspectives futures mentionnées :

• Discrétisation complète : Étendre l'analyse à une méthode entièrement discrète (Espace-Temps) en utilisant des éléments finis et prouver la convergence de ces approximations.
• Généralisation : Étendre la méthode de projection à d'autres critères d'écoulement (ex: Tresca), à des ensembles convexes plus généraux, ou à des paramètres d'écrouissage matriciels.
• Cas non régularisé ($\eta = 0$) : Aborder le cas purement élasto-plastique sans viscosité, ce qui est beaucoup plus difficile car la viscosité fournit actuellement le contrôle nécessaire sur la régularité spatiale de la vitesse ($H^1$).
• Adoucissement (Softening) : Étudier le cas où le paramètre d'écrouissage $a$ est négatif, ce qui brise la structure monotone/coercive utilisée pour l'unicité.

En résumé, cet article combine analyse fonctionnelle avancée et conception de schémas numériques pour résoudre un problème complexe de mécanique des milieux continus, offrant à la fois des garanties théoriques solides et une voie pratique pour la simulation.