Bilan du Groupe de Travail AFGC sur les éléments finis

La non-linéarité d’un problème de mécanique vient de ce que les coefficients apparaissant dans l’équation d’équilibre dépendent de la déformation du solide à l’équilibre. Autrement dit, l’équation d’équilibre est généralement une équation de type implicite.

Il existe plusieurs catégories de non-linéarités pour les problèmes de mécanique statique :

• Les non-linéarités matérielles : c’est le cas dès lors que la loi de comportement n’est pas linéaire ou que la réponse dépend de l’histoire du chargement. Autrement dit, la contrainte n’est pas une fonction linéaire de la déformation. L’exemple le plus fréquent en génie civil est celui des matériaux sollicités au-delà de leur limite élastique, et qui développent un comportement élasto-plastique. Ce comportement se caractérise par une dépendance de la raideur du matériau vis-à-vis de son état de contrainte.
• Les non-linéarités géométriques : c’est le cas dès lors qu’on travaille en grands déplacements ou en grandes déformations. Dans le premier cas, on ne peut plus écrire le problème en négligeant les changements de configurations. Dans le second cas, on ne peut plus approcher la déformation par un simple gradient du déplacement.
• Les non-linéarités des conditions aux limites : cela peut être le cas, lorsque chargement est progressif, lorsqu’il y a contact potentiel entre 2 corps, dans le cas de force suiveuse, mais aussi lorsqu’on simule le phasage de construction, ou le lançage d’un tablier, le creusement d’une galerie, la mise en remblai, etc.

Ces types de non-linéarités peuvent être couplés, si le code de calcul le permet, mais cela complexifie d’autant la résolution du problème.

Lorsqu’on résout le problème élément fini, on cherche le champ de déplacement u, tel que les efforts intérieurs L_int soient égaux aux efforts extérieurs L_ext :

, qui est un problème non-linéaire en fonction de u.

Pour résoudre un problème de statique non-linéaire, on utilise généralement un algorithme incrémental. Pour cela le problème est paramétré en fonction d’un paramètre t (qui est un pseudo-temps, contrairement aux calculs dynamiques). Ce paramètre sert à indexer les chargements successifs imposés à la structure. Plus précisément, on va chercher l’état d’équilibre correspondant aux chargements successifs F₁, F₂, …

Ce découpage conduit à résoudre une succession de problèmes quasi-linéaires comme illustrée sur la figure ci-dessous, et à déterminer l’état de la structure à l’instant t (déplacements, déformations, contraintes) en connaissant la solution à l’état t-1. Plus le nombre de découpage est important, meilleure est la précision.

Principe du paramétrage de problème en fonction de t

A chaque incrément t_i, le problème discrétisé à résoudre est : K_i x q_i =F_i où q_i est le vecteur déplacement inconnu sous le chargement imposé F_i. Alors que dans le cas linéaire vu au chapitre 1, la matrice K était explicite, lorsque le problème est non-linéaire, K_i est une matrice dont les termes dépendent implicitement de la valeur de q_i. Déterminer q_i ne peut donc pas se faire directement en inversant la matrice K.

La méthode la plus courante pour résoudre cette équation non-linéaire est d’utiliser un algorithme de type Newton. Le principe est de construire une bonne approximation de la solution de l’équation en considérant son développement de Taylor au premier ordre . On part donc d’un point initial (suffisamment proche de la solution) et on construit par récurrence la suite .

A chaque itération, on évalue le vecteur résidu F(qk) et on poursuit tant que celui-ci n’est pas inférieur (en norme) à une valeur arbitrairement proche de 0. Ce critère de convergence doit être choisi avec soin en fonction de la norme utilisée par le code de calcul (voir §3.3 pour plus de détails).

Remarque : Avec la méthode de Newton, à chaque itération, on doit calculer la matrice tangente au point considéré . Le coût de calcul de cette matrice peut être coûteux en temps. Si l’utilisation de cette matrice permet d’avoir une convergence quadratique (donc en peu d’itérations), il n’est pas indispensable d’utiliser cette matrice. D’autres stratégies peuvent être appliquées pour estimer cette matrice. On parle alors de méthode de quasi-Newton. Ainsi il est ainsi envisageable d’utiliser la matrice tangente sans la réactualiser à chaque itération mais aussi d’utiliser la matrice élastique (figure b) ou la matrice sécante dans les cas de modèle endommagement. Une illustration des itérations successives selon la matrice utilisée est proposée sur la figure ci-dessous.

Illustration de la méthode de Newton ou quasi-Newton (matrice élastique)

En règle générale, l’utilisation de la matrice tangente permet une convergence plus rapide (en moins d’itérations), mais les variantes peuvent être plus efficaces ou plus robustes selon les situations.

La méthode de résolution étant itérative, le processus doit être interrompu lorsque le critère d’arrêt est atteint, c’est-à-dire que l’on a vérifié qu’une grandeur (plusieurs définitions sont possibles) devient négligeable. L’algorithme global est le suivant :

en désignant les incréments, i indexant les itérations de Newton et ε étant une valeur positive, proche de 0.

Remarque : l’algorithme de Newton est utilisé pour résoudre l’équilibre à chaque pas de temps. Il peut être également utilisé pour trouver la contrainte en chaque point de Gauss (à toutes les itérations du problème de Newton à l’échelle globale) lorsque la loi de comportement le nécessite.

A suivre : C.2 Pourquoi exécuter des calculs non-linéaires ?

Lien vers la bibliographie

Laissez nous vos commentaires pour améliorer le site !

Retour vers l'accueil