Théorie de la plasticité
📝 Situation du Projet
Vous intégrez le Bureau d'Études Structures d'un grand major du BTP en charge de la conception du complexe logistique "Titan". Ce projet se caractérise par des portées importantes et une exigence stricte d'optimisation des tonnages d'acier pour réduire l'empreinte carbone et les coûts.
La structure principale est composée de poutres continues en acier supportant la toiture. Traditionnellement, le dimensionnement se fait au domaine élastique (ELU Élastique). Cependant, pour ce projet spécifique, le Maître d'Ouvrage a validé l'utilisation de la réserve de capacité plastique des matériaux, conformément à l'Eurocode 3. Cela permet d'exploiter la ductilité de l'acier au-delà de sa limite d'élasticité, acceptant la formation de rotules plastiques avant la ruine complète.
En tant qu'Ingénieur Structure Senior, vous devez déterminer la charge ultime théorique (\(q_{\text{u}}\)) que peut supporter la poutre de toiture avant de se transformer en mécanisme de ruine. Vous utiliserez le théorème cinématique de l'analyse limite pour justifier le dimensionnement du profilé IPE choisi.
"Attention, l'analyse plastique n'est valide que si la section est de Classe 1 (compacte). Assure-toi que le profilé IPE 300 peut développer une rotule plastique sans voilement local prématuré des semelles ou de l'âme."
Pour mener à bien cette étude de capacité portante ultime, nous nous basons sur les propriétés géométriques normalisées des profilés européens et les caractéristiques mécaniques de l'acier de construction standard.
📚 Référentiel Normatif
Eurocode 3 (EN 1993-1-1) Théorie de la Plasticité| GÉOMÉTRIE & MATÉRIAU | |
| Nuance d'Acier | S275 (Limite élastique \(f_{\text{y}} = 275 \text{ MPa}\)) |
| Hauteur de section (\(h\)) | 300 mm |
| Largeur de section (\(b\)) | 150 mm |
| Épaisseur âme (\(t_{\text{w}}\)) | 7.1 mm |
| Épaisseur semelle (\(t_{\text{f}}\)) | 10.7 mm |
| CARACTÉRISTIQUES DE FLEXION | |
| Module d'inertie élastique (\(W_{\text{el,y}}\)) | 557 cm³ |
| Module d'inertie plastique (\(W_{\text{pl,y}}\)) | 628 cm³ |
| Coefficient partiel de sécurité (\(\gamma_{\text{M0}}\)) | 1.0 (Hypothèse Eurocode) |
📐 Configuration du Système
- Type de poutre : Poutre continue sur 3 appuis (Hyperstatique de degré 1).
- Appuis : A (Simple), B (Intermédiaire), C (Simple).
- Portée par travée (\(L\)) : 6.00 m.
- Chargement : Uniformément réparti (\(q\)) sur les deux travées.
⚖️ Hypothèses de Calcul
E. Protocole de Résolution
Pour déterminer la charge de ruine, nous allons suivre une approche méthodique basée sur les théorèmes fondamentaux de l'analyse limite.
Calcul de la Capacité Plastique
Détermination du Moment Plastique de la section IPE 300, représentant la résistance maximale en rotation.
Identification du Mécanisme
Analyse de la structure pour localiser les positions potentielles des rotules plastiques et définir la cinématique de ruine.
Principe des Travaux Virtuels (PTV)
Égalisation du travail des forces extérieures et de l'énergie de dissipation interne pour isoler la charge limite.
Vérification & Conclusion
Comparaison avec la limite élastique pour quantifier le gain apporté par l'analyse plastique.
Théorie de la plasticité
🎯 Objectif
L'objectif primordial de cette première étape est de quantifier la capacité maximale de flexion de la poutre IPE 300 avant qu'elle ne devienne une "charnière" mécanique. Contrairement au calcul élastique où l'on s'arrête dès que la fibre la plus sollicitée atteint la limite d'élasticité (\(f_{\text{y}}\)), le calcul plastique suppose que l'intégralité de la section transversale a atteint la plastification (diagramme de contrainte rectangulaire bi-bloc). Cette valeur, notée \(M_{\text{pl,Rd}}\), est la "clé de voûte" de toute l'analyse cinématique qui suivra.
📚 Référentiel
EN 1993-1-1 Art 6.2.5Avant de lancer le calcul, il est crucial de comprendre la physique du phénomène. En flexion élastique, les contraintes varient linéairement de zéro (axe neutre) à \(f_{\text{y}}\) (fibres extrêmes). Le moment résistant dépend de \(W_{\text{el}}\). En régime plastique, on accepte que les fibres intérieures atteignent aussi \(f_{\text{y}}\) progressivement. Lorsque toute la section est plastifiée, le moment ne peut plus augmenter : une rotation libre se produit à moment constant. C'est la rotule plastique. La capacité de résistance dépend alors de \(W_{\text{pl}}\), qui est toujours supérieur à \(W_{\text{el}}\) (pour un IPE, le gain est d'environ 12 à 15%).
Dans la théorie de la plasticité parfaite : 1. Le matériau suit la loi de Hooke jusqu'à \(\sigma = f_{\text{y}}\). 2. Au-delà, la déformation \(\epsilon\) augmente à contrainte constante \(\sigma = f_{\text{y}}\) (palier plastique). 3. En flexion pure, la section se divise en deux zones (comprimée et tendue) intégralement plastifiées. L'axe neutre plastique divise la section en deux aires égales (pour une section homogène). 4. Le moment plastique correspond au couple formé par la résultante de compression et la résultante de traction : \(M_{\text{pl}} = F \cdot z = (A/2 \cdot f_{\text{y}}) \cdot z\), ou plus simplement \(M_{\text{pl}} = W_{\text{pl}} \cdot f_{\text{y}}\).
Schéma : Comparaison des Distributions de Contraintes
📋 Données d'Entrée :
| Paramètre | Valeur |
|---|---|
| Module plastique (\(W_{\text{pl,y}}\)) | 628 cm³ (\(628 \times 10^{-6} \text{ m}^3\)) |
| Limite élastique (\(f_{\text{y}}\)) | 275 MPa (\(275 \times 10^6 \text{ Pa}\)) |
L'erreur classique est le mélange des unités (cm et m, MPa et Pa). Travaillez toujours en unités SI (mètres, Newtons) ou convertissez directement : \(1 \text{ cm}^3 = 10^{-6} \text{ m}^3\) et \(1 \text{ MPa} = 1 \text{ N/mm}^2 = 10^6 \text{ N/m}^2\). Le résultat sera en Newton-mètre (Nm), facile à convertir en kNm.
📝 Calcul Détaillé
Nous appliquons directement la formule. C'est ce moment \(M_{\text{pl}}\) qui sera injecté dans l'équation de dissipation d'énergie interne (\(W_{\text{int}}\)) plus tard.
Manipulation détaillée de la formule :
Pour obtenir le moment résultant, on intègre les contraintes sur la section. En plasticité parfaite, l'effort normal résultant est nul, donc la zone comprimée égale la zone tendue.
1. Détermination de la résistance brute en Nm :
On multiplie le module plastique (converti en m³) par la limite d'élasticité (en Pa).
La valeur brute est de 172 700 Newtons-mètres.
2. Conversion en kNm :
Pour plus de lisibilité dans les calculs de structure, nous convertissons en kNm (division par 1000).
C'est la valeur de référence pour toutes les rotules plastiques de la structure.
✅ Interprétation Globale : La section IPE 300, en acier S275, peut développer un moment de flexion maximal de 172.7 kNm. Au-delà de cette valeur, la courbure augmente indéfiniment sans augmentation de moment : c'est la formation d'une "rotule plastique".
Comparons avec le moment élastique limite : \(M_{\text{el}} = W_{\text{el}} \cdot f_{\text{y}} = 557 \cdot 275 = 153.2 \text{ kNm}\). Le rapport \(M_{\text{pl}}/M_{\text{el}} = 172.7/153.2 \approx 1.13\). Ce ratio de 1.13 est typique pour les profils en I fléchis selon l'axe fort (le facteur de forme théorique est souvent entre 1.10 et 1.15). Le résultat est donc cohérent.
Ce calcul n'est valide QUE SI la section est de Classe 1 (ou éventuellement Classe 2 selon l'Eurocode). Si la section était de Classe 3 ou 4 (parois fines), le voilement local des semelles ou de l'âme se produirait AVANT l'atteinte du moment plastique complet, invalidant ce calcul et obligeant à utiliser un calcul élastique ou réduit.
🎯 Objectif
Il s'agit maintenant de prédire géométriquement comment la poutre va s'effondrer. En analyse limite, la ruine survient lorsque suffisamment de rotules plastiques se sont formées pour transformer la structure (initialement stable et hyperstatique) en un mécanisme (système instable à mouvements libres).
📚 Référentiel
Théorème Cinématique (Majorant)Notre poutre est continue sur 3 appuis (2 travées). Elle est hyperstatique de degré 1. Il faut donc former 2 rotules plastiques dans une travée pour créer un mécanisme de ruine local (mécanisme de poutre). Où vont-elles se former ? 1. Sur l'appui intermédiaire B, là où le moment élastique est maximal (en valeur absolue). 2. Dans la travée, proche du point de moment maximal positif. Comme le chargement est symétrique sur les deux travées, la ruine peut survenir indifféremment dans la travée de gauche ou de droite. Nous allons étudier le mécanisme de la travée AB.
Un mécanisme est une configuration cinématique où la structure devient mobile. Le nombre de rotules nécessaires \(n\) est donné par \(n = h + 1\), où \(h\) est le degré d'hyperstaticité. Ici, pour une travée encastrée-appuyée (équivalent de notre travée continue isolée), \(h=1\), donc il faut \(1+1=2\) rotules plastiques. Les rotules se forment aux points d'extremum du diagramme de moment fléchissant.
Pour un chargement uniforme sur une travée propped-cantilever (appuyée-encastrée), la rotule en travée ne se forme pas exactement au milieu, mais à environ \(0.414L\) de l'appui simple. Cependant, pour cet exercice pédagogique, nous ferons l'hypothèse simplificatrice standard d'une rotule à mi-travée (\(L/2\)). Cela simplifie la géométrie sans induire une erreur massive (l'erreur est sécuritaire pour le dimensionnement, mais ici nous cherchons la charge limite théorique, l'approximation surestime légèrement la capacité).
📋 Données Géométriques
| Paramètre | Valeur |
|---|---|
| Longueur Travée (\(L\)) | 6.00 m |
| Position Rotule 1 | Appui B (Encastrement) |
| Position Rotule 2 | Mi-travée (\(L/2 = 3.00\) m) |
Lorsque la rotule est supposée au milieu, la cinématique est symétrique par rapport à cette rotule dans la travée. Cela simplifie grandement le calcul des angles : \(\theta_{\text{gauche}} = \theta_{\text{droite}}\).
Fig 2. Cinématique du mécanisme de ruine (Travée AB)
📝 Calcul Détaillé
Manipulation détaillée de la formule :
Pour les petits angles, \(\tan(\theta) \approx \theta\). Dans le triangle rectangle formé par la position initiale, la position déformée et la verticale sous la rotule :
Pour l'angle au centre, la géométrie du triangle isocèle impose que l'angle extérieur est la somme des angles intérieurs opposés :
1. Relation Déplacement-Rotation :
On exprime le déplacement vertical \(\delta\) en fonction de l'angle \(\theta\) (supposé petit).
C'est la compatibilité géométrique de base.
2. Rotation des Rotules :
On identifie l'angle de rotation total subi par chaque section critique.
Notez bien que la rotule centrale "plie" deux fois plus que les appuis car elle connecte deux segments inclinés de \(\theta\).
✅ Interprétation Globale : Le mécanisme identifié est un "mécanisme de poutre" classique. La poutre se brise en deux segments rigides articulés entre eux et à l'appui B. L'appui A tourne librement.
Le mécanisme respecte bien les degrés de liberté : 1 degré de liberté cinématique (\(\theta\)) pour le mécanisme complet. Les signes des rotations sont compatibles avec le sens de la charge (descente).
Ne pas oublier que l'appui A est une rotule simple (physique). Il ne consomme pas d'énergie plastique car le moment y est nul par définition. Seules les sections où un moment plastique se développe et tourne consomment de l'énergie.
🎯 Objectif
C'est l'étape de calcul final. Le Théorème Cinématique stipule que pour tout mécanisme cinématiquement admissible, la charge calculée en égalant la puissance des forces extérieures à la puissance de dissipation interne est supérieure ou égale à la charge réelle de ruine. En cherchant le mécanisme le plus critique, on trouve la charge ultime.
📚 Référentiel
Équation : \(W_{\text{ext}} = W_{\text{int}}\)Nous allons écrire l'équation d'équilibre énergétique. 1. Travail Extérieur (\(W_{\text{ext}}\)) : C'est le travail de la charge répartie \(q\) descendant le long du déplacement de la poutre. Pour une charge uniforme sur un diagramme de déplacement triangulaire, le travail est équivalent à la charge totale (\(qL\)) multipliée par le déplacement moyen (\(\delta/2\)). 2. Travail Interne (\(W_{\text{int}}\)) : C'est l'énergie dissipée par l'acier en se plastifiant. Elle ne se produit QUE dans les rotules plastiques. Elle vaut le produit du Moment Plastique par l'angle de rotation de la rotule (\(M_{\text{pl}} \cdot \theta_{\text{rotule}}\)). L'appui A ne dissipe rien (rotule parfaite). En égalant les deux, \(\theta\) va se simplifier (le résultat ne dépend pas de l'amplitude du mouvement virtuel), et nous isolerons \(q_{\text{u}}\).
Le PTV stipule que pour un système en équilibre, le travail des forces extérieures pour tout déplacement virtuel compatible avec les liaisons est égal au travail des forces intérieures (dissipation). \(W_{\text{ext}} = \int \vec{F} \cdot \vec{du}\) et \(W_{\text{int}} = \sum \vec{M} \cdot \vec{d\theta}\). Ici, la dissipation est localisée uniquement aux rotules.
📋 Données pour le calcul :
| Variable | Expression |
|---|---|
| Moment Plastique | \(M_{\text{pl}} = 172.7\) kNm |
| Longueur | \(L = 6.00\) m |
| Rotation virtuelle | \(\theta\) (arbitraire) |
Pour calculer le travail d'une charge répartie sur un déplacement triangulaire, pas besoin d'intégrale complexe ! Utilisez la formule : \(W = \text{Résultante totale} \times \text{Déplacement du centre de gravité}\). Pour un triangle, le déplacement moyen est \(\delta/2\).
Schéma : Calcul Graphique du Travail (Aire du déplacement)
📝 Calcul Détaillé
Manipulation détaillée de la formule :
Calcul de l'intégrale du travail extérieur :
On remplace \(\delta\) par \((L/2)\theta\) :
1. Calcul du Travail Interne (Dissipation) :
On somme les produits \(M_{\text{pl}} \times \text{Angle}\) pour toutes les rotules plastiques (B et Travée).
On obtient une dissipation totale de 3 fois le moment plastique.
2. Calcul du Travail Extérieur :
Intégration de la charge \(q\) sur le champ de déplacement triangulaire défini par \(\delta\).
Le travail dépend du carré de la longueur, ce qui est classique pour une charge répartie.
3. Résolution de l'Équation d'Équilibre :
Nous égalons \(W_{\text{ext}} = W_{\text{int}}\) pour trouver la charge ultime \(q_{\text{u}}\).
Nous avons obtenu la formule analytique de la charge de ruine.
4. Application Numérique :
Avec \(M_{\text{pl}} = 172.7 \text{ kNm}\) et \(L = 6.0 \text{ m}\).
Le calcul numérique final donne la valeur limite.
✅ Interprétation Globale : La structure s'effondrera théoriquement (formation d'un mécanisme complet) lorsque la charge répartie atteindra 57.57 kN par mètre linéaire. Toute charge inférieure est stable.
L'unité est bien une charge linéique (kN/m). La formule \(12 M/L^2\) ressemble aux formules classiques RDM (\(8 M/L^2\)), mais avec un coefficient plus élevé (12 vs 8), ce qui indique une capacité supérieure, logique pour une analyse plastique.
Attention, ce calcul suppose que la rotule se forme exactement au milieu. Le calcul exact (optimisation de la position x) donnerait un facteur \(11.66\) au lieu de \(12\). Notre calcul surestime donc la résistance de \((12-11.66)/11.66 \approx 2.9\%\). Pour un pré-dimensionnement ou une vérification EXE rapide, c'est acceptable, mais pour un calcul final strict, il faudrait utiliser la position exacte.
🎯 Objectif
Pour justifier le choix de l'analyse plastique auprès du bureau de contrôle ou du client, il est impératif de quantifier le "gain" réalisé par rapport à un calcul élastique traditionnel. Combien de capacité portante avons-nous "débloqué" en acceptant la plastification ?
📚 Référentiel
Calcul élastique RDMDans un calcul élastique classique d'une poutre continue sur 2 travées, le moment maximal (négatif) se situe sur l'appui central B et vaut \(M_{\text{el,max}} = qL^2/8\). La ruine élastique est atteinte dès que ce moment touche la limite \(M_{\text{el,Rd}} = W_{\text{el}} \cdot f_{\text{y}}\). En comparant le \(q_{\text{el}}\) issu de cette formule avec notre \(q_{\text{u}}\) (qui est en réalité \(q_{\text{pl}}\)), nous verrons l'efficacité de la redistribution des efforts.
En élastique, la répartition des moments dépend de la rigidité (\(EI\)). En plastique, elle dépend de la résistance. Quand l'appui B plastifie, il ne prend plus de moment supplémentaire, mais il continue de porter son moment limite. Le surplus de charge est alors reporté vers la travée, qui n'est pas encore saturée. C'est ce phénomène de redistribution qui offre un gain de capacité.
Schéma : Redistribution du Moment (Élastique vs Plastique)
📋 Données d'Entrée :
| Paramètre | Valeur |
|---|---|
| Charge Ultime Plastique (\(q_{\text{u}}\)) | 57.57 kN/m |
| Module Élastique (\(W_{\text{el}}\)) | 557 cm³ |
Pour une poutre continue 2 travées égales chargées, le moment sur appui est \(0.125 qL^2\) (soit \(qL^2/8\)) et le moment max en travée est environ \(0.07 qL^2\). L'appui est toujours dimensionnant en élastique.
📝 Calcul Détaillé
Manipulation détaillée de la formule :
Le moment sur appui d'une poutre continue se démontre par l'équation des trois moments (Clapeyron) :
1. Calcul de la résistance élastique limite (Mel) :
On utilise le module élastique \(W_{\text{el}}\) et la limite d'élasticité.
C'est le moment où la première fibre plastifie.
2. Calcul de la charge limite élastique (qel) :
On cherche \(q\) tel que le moment sur appui (le plus sollicité) atteigne \(M_{\text{elastic}}\).
La poutre "élastique" ruine à 34 kN/m.
3. Calcul du Gain :
Comparaison finale.
Le bénéfice est considérable.
✅ Interprétation Globale : Le passage au calcul plastique permet une augmentation spectaculaire de la charge admissible (+69%). Cela est dû à deux facteurs : le facteur de forme (\(W_{\text{pl}}/W_{\text{el}} \approx 1.12\)) qui apporte environ 12% de gain, et surtout la redistribution des moments (formation de la rotule sur appui puis transfert vers la travée) qui apporte le reste.
Le gain est élevé mais réaliste pour une poutre hyperstatique. Pour une poutre isostatique, le gain se limiterait au rapport \(W_{\text{pl}}/W_{\text{el}}\) (soit 12%). C'est l'hyperstaticité qui permet la redistribution majeure.
Ce calcul est théorique. Dans la réalité : 1. Il faut vérifier l'effort tranchant (interaction M-V) qui peut réduire le moment plastique. 2. Le déversement (instabilité latérale) doit être impérativement empêché par des pannes anti-déversement, car les zones plastifiées n'ont plus de rigidité torsionnelle.
📄 Livrable Final (Note de Calculs EXE)
STRUCTURE
- Profilé : IPE 300 (Classe 1)
- Acier : S275 (\(f_{\text{y}} = 275\) MPa)
- Module Plastique \(W_{\text{pl,y}}\) : \(628 \text{ cm}^3\)
| Moment Résistant Plastique | 172.7 kNm |
Hypothèse : Mécanisme de ruine de poutre (Beam Mechanism) sur travée de rive.
Gain vs Élastique : +69%
Jean D.
Dr. A. Einstein
Laisser un commentaire