Calcul du Flambement Critique Global d'une Structure Élancée par Analyse Non-Linéaire (GMNIA) : Guide Pratique et Méthodologie

Le dimensionnement des structures élancées requiert une maîtrise approfondie des phénomènes d'instabilité élastoplastique. Parmi ces phénomènes, le flambement global représente un mode de ruine critique pour les poteaux, les mâts, les arcs et les treillis soumis à la compression.

Alors que l'analyse linéaire classique (première théorie) ou l'analyse de bifurcation élastique (analyse linéaire de flambement ou LBA) fournissent des estimations trop optimistes ou théoriques, l'analyse non-linéaire avancée (GMNIA : Geometrically and Materially Non-linear Analysis with Imperfections) s'impose aujourd'hui comme la méthode de référence pour capturer le comportement réel des structures élancées.

1. Fondements théoriques du flambement et limites des méthodes linéaires

L'approche classique d'Euler

Le flambement élastique d'une poutre droite idéale est régi par la charge critique d'Euler (N_cr), obtenue en résolvant l'équation différentielle de l'équilibre élastique :

N_cr = (π² × E × I) / (L_f)²

Où :

• E : module d'Young du matériau.
• I : moment d'inertie de la section par rapport à l'axe de flambement.
• L_f : longueur de flambement de la structure, fonction des conditions aux appuis.

Pourquoi l'analyse linéaire est insuffisante

Les structures réelles s'écartent des modèles théoriques par :

• La présence d'imperfections géométriques initiales (hors-plomb, courbures initiales).
• Des contraintes résiduelles issues des procédés de fabrication (laminage, soudage).
• Le comportement non-linéaire du matériau (plasticité) bien avant l'atteinte de la charge élastique critique.

2. Méthodologie de l'analyse non-linéaire (GMNIA)

La réalisation d'un calcul de stabilité par analyse non-linéaire suit une méthodologie rigoureuse en plusieurs étapes clés :

Étape A : Introduction des imperfections géométriques équivalentes

Les imperfections initiales doivent être modélisées explicitement dans le modèle éléments finis. Conformément aux Eurocodes (notamment l'Eurocode 3 pour l'acier), l'amplitude de la flèche initiale équivalente e₀ est choisie en fonction de la classe de la section et de la courbe de flambement.

Étape B : Prise en compte des contraintes résiduelles

Le modèle constitutif du matériau doit intégrer une loi de comportement élastoplastique avec écrouissage. Les contraintes résiduelles sont injectées soit par l'application de champs de contraintes initiaux, soit par l'utilisation de courbes contrainte-déformation adaptées.

Étape C : Résolution numérique

L'équilibre de la structure est résolu de manière incrémentale sous l'effet des charges. La méthode de Newton-Raphson est utilisée pour suivre le chemin d'équilibre.

3. Formulation mathématique des non-linéarités

Dans un code de calcul aux éléments finis, la matrice de raideur tangente globale [K_t] est décomposée ainsi :

[K_t] = [K_e] + [K_g] + [K_pl]

Où :

• [K_e] : raideur élastique.
• [K_g] : raideur géométrique (effets du second ordre).
• [K_pl] : raideur plastique.

4. Critères d'évaluation

L'analyse fournit une courbe d'équilibre décrivant la relation entre le déplacement caractéristique et le facteur de charge. Le sommet de la courbe (point limite) indique la capacité portante maximale.

Note technique : L'utilisation de logiciels de calcul par éléments finis (FEA) avancés est indispensable pour mener à bien ces analyses.