Méthodes de solution itérative pour les études non linéaires

Etudes statiques non linéaires

Dans une analyse statique non linéaire, l'ensemble d'équations de base à résoudre à n'importe quel pas de temps, t+Δt, est :

t+Δt{R} - t+Δt{F} = 0,

où :

t+Δt{R} = Vecteur de chargements nodaux appliqués en externe

t+Δt{F} = Vecteur de forces nodales générées en interne.

Dans la mesure où les forces nodales internes t+Δt{F} dépendent des déplacements nodaux au temps t+Δt, t+Δt{U}, une méthode itérative doit être utilisée. Les équations suivantes représentent les grandes lignes de base d'un schéma itératif permettant de résoudre les équations d'équilibre à un certain pas de temps, t+Δt,

{ΔR}(i-1) = t+Δt{R} - t+Δt{F}(i-1)

t+Δt[K](i) {ΔU}(i) = {ΔR}(i-1)

t+Δt{U}(i) = t+Δt{U}(i-1) + {ΔU}(i)

t+Δt{U}(0) = t{U}; t+Δt{F}(0) = t{F}

où :

t+Δt{R} = Vecteur de chargements nodaux appliqués en externe

t+Δt{F}(i-1) = Vecteur de forces nodales générées en interne à l'itération (i)

{ΔR}(i-1) = Vecteur de chargement déséquilibré à l'itération (i)

{ΔU}(i) = Vecteur de déplacements nodaux incrémentiels à l'itération (i)

t+Δt{U}(i) = Vecteur de déplacements totaux à l'itération (i)

t+Δt[K](i) = Matrice Jacobienne (raideur tangente) à l'itération (i).

Il existe différents schémas pour l'exécution des itérations ci-dessus. Une description succincte de deux méthodes de type Newton est fournie ci-dessous :

Méthodes de solution itérative - Schéma de Newton-Raphson (NR)

Dans ce schéma, la matrice de raideur tangentielle est formée et décomposée à chaque itération dans un pas particulier (voir la figure ci-dessous). Le taux de convergence de la méthode NR est particulièrement élevé et est quadratique. Toutefois, dans la mesure où la raideur tangentielle est formée et décomposée à chaque itération, qui peuvent être particulièrement chères pour les grands modèles, il peut être plus avantageux d'utiliser une autre méthode itérative.

Méthodes de solution itérative - Schéma de Newton-Raphson modifié (MNR)

Dans ce schéma, la matrice de raideur tangentielle est formée et décomposée au début de chaque étape (ou comme spécifié lors de la définition des propriétés de l'étude) et est utilisée dans toutes les itérations (voir la figure ci-dessous).