Calcul des Contraintes Principales (État Plan)
📝 Situation du Projet
Au sein du bureau d'études international "Structure & Conception", vous êtes mandaté sur le chantier de rénovation du Viaduc de la Vallée (Ouvrage V12). Ce pont historique en treillis métallique, datant des années 1960, doit subir une requalification pour autoriser le passage de convois exceptionnels de classe E. La structure, exposée depuis des décennies aux cycles de gel/dégel et aux vibrations du trafic autoroutier, fait l'objet d'une surveillance accrue.
Les inspections par ultrasons et magnétoscopie ont révélé des micro-fissures superficielles sur les goussets d'assemblage des nœuds inférieurs (zone N-24), à la jonction entre la membrure tendue et les diagonales comprimées. Un modèle numérique global a permis d'extraire le tenseur des contraintes local en un point critique \(P\), situé en bordure de soudure. Ce point est soumis à un état de contrainte plan complexe, combinant la traction axiale de la membrure, la flexion secondaire induite par l'excentrement des assemblages, et le cisaillement transmis par les rivets.
Votre expertise est requise pour transformer cet état de contrainte brut (exprimé dans le repère géométrique de la pièce) en contraintes principales. Cette étape est indispensable pour appliquer le critère de Von Mises et vérifier si l'acier S355 J2 en place conserve une marge de sécurité suffisante (Coefficient > 1.5) sous les nouvelles charges d'exploitation pondérées (ELU). Une défaillance de ce nœud entraînerait une redistribution d'efforts catastrophique pour la stabilité globale de la travée centrale.
En tant qu'Ingénieur Structure Senior, vous devez déterminer analytiquement les contraintes principales (\(\sigma_1, \sigma_2\)) et leur orientation (\(\alpha_{\text{p}}\)) à partir des relevés tensiométriques. Vous conclurez ensuite sur la viabilité de l'assemblage en calculant le critère de Von Mises et en le comparant à la limite élastique de l'acier.
"Attention : L'état de contrainte initial est fourni dans le repère local (x, y) lié à la géométrie de la membrure. Une erreur fréquente consiste à confondre cet état avec les contraintes principales. Vérifiez scrupuleusement les signes des contraintes normales (Traction > 0, Compression < 0) avant tout calcul."
Les paramètres ci-dessous définissent le cadre normatif et matériel de l'étude. Ils sont issus des notes de calculs préliminaires (Réf. NC-V12-001) et des certificats matières.
📚 Référentiel Normatif
Eurocode 3 (EN 1993-1-1)Théorie de l'Élasticité Linéaire| PROPRIÉTÉS MÉCANIQUES | |
| Limite d'élasticité (\(f_y\)) | 355 \(\text{ MPa}\) |
| Résistance à la rupture (\(f_u\)) | 510 \(\text{ MPa}\) |
| CONSTANTES ÉLASTIQUES | |
| Module de Young (\(E\)) | 210 000 \(\text{ MPa}\) |
| Coefficient de Poisson (\(\nu\)) | 0.3 |
📐 Géométrie Locale
- Épaisseur du gousset (\(t\)): 20 mm
- Distance à la fibre neutre: 145 mm
- Type de soudure: Angle à pénétration partielle
⚖️ Tenseur des Contraintes (Point P)
Valeurs extraites au nœud N-24 sous combinaison ELU fondamentale :
| Donnée | Symbole | Valeur | Unité |
|---|---|---|---|
| Contrainte Normale (Axe X) | \(\sigma_x\) | + 80 | \(\text{MPa}\) |
| Contrainte Normale (Axe Y) | \(\sigma_y\) | - 40 | \(\text{MPa}\) |
| Contrainte de Cisaillement | \(\tau_{xy}\) | + 50 | \(\text{MPa}\) |
| Limite Élastique Acier | \(f_y\) | 355 | \(\text{MPa}\) |
E. Protocole de Résolution
Pour garantir la fiabilité du diagnostic structurel, nous suivrons une méthodologie rigoureuse basée sur la diagonalisation du tenseur des contraintes.
Calcul des Invariants Géométriques
Détermination du centre \(C\) (contrainte moyenne) et du rayon \(R\) (cisaillement maximal) du Cercle de Mohr associé à l'état de contrainte.
Détermination des Contraintes Principales
Calcul analytique des contraintes normales extrêmes (\(\sigma_1\) et \(\sigma_2\)) correspondant aux valeurs propres du tenseur.
Analyse de l'Orientation Principale
Calcul de l'angle \(\alpha_{\text{p}}\) définissant la direction des facettes principales par rapport au repère géométrique du gousset.
Validation du Critère de Ruine
Calcul de la contrainte équivalente de Von Mises et vérification du coefficient de sécurité vis-à-vis de la limite élastique.
Calcul des Contraintes Principales (État Plan)
🎯 Objectif Scientifique
L'objectif primordial de cette première étape est de s'affranchir du système de coordonnées local pour identifier les propriétés intrinsèques de l'état de contrainte. En Résistance des Matériaux, les contraintes \(\sigma_x\), \(\sigma_y\) et \(\tau_{xy}\) ne sont que des "ombres" projetées sur des plans spécifiques. Nous cherchons ici à construire le Cercle de Mohr, un outil graphique et analytique puissant qui synthétise l'infinité des états de contrainte possibles en un seul objet géométrique défini par son centre \(C\) et son rayon \(R\).
📚 Référentiel Théorique
Théorie du Cercle de Mohr Invariants du TenseurAvant de se lancer dans les calculs, il faut comprendre la signification physique de \(C\) et \(R\). Le centre \(C\) représente la contrainte moyenne, c'est-à-dire la pression hydrostatique (ou tension isotrope) qui tend à changer le volume de la matière sans la déformer (changement de forme). Le rayon \(R\), quant à lui, est directement lié à l'énergie de distorsion : plus le rayon est grand, plus le cisaillement maximal possible est important, et donc plus le risque de glissement des plans atomiques (plasticité) est élevé. Calculer ces deux valeurs revient à "prendre l'empreinte digitale" de l'état de contrainte, indépendamment de l'orientation de notre jauge de mesure.
Dans un repère orthonormé \((\sigma, \tau)\), les facettes verticales et horizontales sont représentées par deux points diamétralement opposés sur le cercle : \(X(\sigma_x, -\tau_{xy})\) et \(Y(\sigma_y, \tau_{xy})\).
- Centre : Le centre \(C\) est le milieu du segment \([XY]\). Son abscisse est donc la moyenne des abscisses :
- Rayon : Le rayon est la distance \(CX\). Dans le triangle rectangle formé par \(C\), \(X\) et la projection de \(X\) sur l'axe horizontal, le côté adjacent vaut \((\sigma_x - \sigma_{\text{moy}}) = \frac{\sigma_x - \sigma_y}{2}\) et le côté opposé vaut \(\tau_{xy}\). Le théorème de Pythagore donne alors :
Moyenne arithmétique des contraintes normales.
Hypoténuse du triangle de construction.
📋 Données d'Entrée pour le Calcul
| Paramètre | Valeur |
|---|---|
| Contrainte Normale X (\(\sigma_x\)) | + 80 \(\text{ MPa}\) |
| Contrainte Normale Y (\(\sigma_y\)) | - 40 \(\text{ MPa}\) |
| Cisaillement (\(\tau_{xy}\)) | + 50 \(\text{ MPa}\) |
Attention aux signes ! Une contrainte de compression est négative. Ici, \(\sigma_y = -40\). Lors du calcul de la différence \(\sigma_x - \sigma_y\), on obtient \(80 - (-40) = 120\). C'est une source d'erreur très fréquente : sous-estimer le rayon en oubliant que la différence s'agrandit quand les signes sont opposés.
📝 Application Numérique Détaillée
Schéma : Construction géométrique du cercle à partir des points X et Y.
1. Détermination de la position du Centre \(C\) :
Nous calculons la contrainte moyenne, qui correspond à l'abscisse du centre du cercle sur l'axe \(\sigma\).
Interprétation : Le centre du cercle est situé à +20 MPa. L'état de contrainte global est donc centré sur une légère traction.
2. Détermination du Rayon \(R\) :
Le rayon représente l'amplitude maximale de variation des contraintes normales et la valeur maximale du cisaillement.
Interprétation : Le rayon est de 78.1 MPa. Cela implique que le cisaillement maximal sur une facette inclinée sera de 78.1 MPa.
✅ Interprétation Globale de l'Étape 1
Nous avons défini géométriquement l'enveloppe de résistance de notre matériau au point P. Le cercle est centré en (20, 0) et a un rayon de 78.1. Cela signifie que les contraintes normales oscilleront entre \(20 - 78.1\) et \(20 + 78.1\) selon l'angle d'observation.
Le rayon obtenu (78.1) est supérieur à la contrainte de cisaillement initiale (50), ce qui est logique : le plan de coupe initial n'est pas le plan de cisaillement maximal. La moyenne (20) est bien comprise strictement entre -40 et +80.
Vérifiez toujours que le terme sous la racine est positif (somme de carrés). Si vous trouvez un rayon négatif ou imaginaire, vous avez fait une erreur de calcul.
🎯 Objectif Scientifique
Il s'agit maintenant de trouver les valeurs extrêmes des contraintes normales. Dans la conception mécanique, les matériaux fragiles cassent souvent sous l'effet de la contrainte principale majeure (traction max), tandis que les instabilités (flambement local) peuvent être initiées par la contrainte principale mineure (compression max). Ce calcul permet de connaître les "pires" cas de traction et de compression que subit réellement la matière, au-delà des axes géométriques.
📚 Référentiel Théorique
DiagonalisationValeurs PropresGéométriquement, les contraintes principales correspondent aux intersections du cercle de Mohr avec l'axe horizontal (\(\tau = 0\)). En effet, par définition, un plan principal est un plan où le cisaillement est nul (propriété fondamentale des vecteurs propres). Puisque nous connaissons le centre \(C\) et le rayon \(R\), les points extrêmes sont simplement situés à une distance \(R\) de part et d'autre de \(C\). C'est beaucoup plus rapide que de diagonaliser la matrice 2x2.
Les contraintes principales \(\sigma_1\) et \(\sigma_2\) sont les solutions de l'équation caractéristique du tenseur des contraintes :
La résolution de cette équation du second degré :
mène exactement aux mêmes formules que l'approche géométrique de Mohr : \(\lambda = \text{Moyenne} \pm \text{Rayon}\).
Les extrema sont situés à une distance \(R\) de part et d'autre du centre \(C\).
📋 Rappel des Valeurs Intermédiaires
| Paramètre | Valeur |
|---|---|
| Centre (\(\sigma_{\text{moy}}\)) | 20 \(\text{ MPa}\) |
| Rayon (\(R\)) | 78.10 \(\text{ MPa}\) |
Pour vérifier rapidement vos calculs : la somme des contraintes principales \(\sigma_1 + \sigma_2\) doit toujours être égale à la somme des contraintes initiales \(\sigma_x + \sigma_y\). C'est le premier invariant du tenseur (la Trace).
📝 Application Numérique Détaillée
Schéma : Visualisation des contraintes principales (intersections avec l'axe horizontal).
1. Calcul de la Contrainte Principale Majeure (\(\sigma_1\)) :
C'est la contrainte la plus positive (traction maximale) subie par le gousset. Elle détermine le risque de rupture par traction.
Interprétation : La traction maximale réelle dans le matériau est de 98.1 MPa, ce qui est supérieur à la contrainte de traction initiale (\(\sigma_x = 80\) MPa).
2. Calcul de la Contrainte Principale Mineure (\(\sigma_2\)) :
C'est la contrainte la plus négative (compression maximale) subie par le gousset.
Interprétation : La compression maximale est de -58.1 MPa, valeur plus sévère en magnitude que la compression initiale (\(\sigma_y = -40\) MPa).
✅ Interprétation Globale de l'Étape 2
L'analyse révèle que le matériau subit en réalité des contraintes plus sévères que celles mesurées initialement selon les axes x et y. La traction monte à près de 100 MPa et la compression à près de 60 MPa.
On vérifie l'invariant de la trace (la somme des éléments diagonaux de la matrice est constante).
Initial : \(\sigma_x + \sigma_y = 80 - 40 = 40\).
Principal : \(\sigma_1 + \sigma_2 = 98.1 - 58.1 = 40\).
L'invariant est respecté, le calcul est correct.
Ne jamais confondre \(\sigma_1\) avec la contrainte équivalente. \(\sigma_1\) est une composante vectorielle réelle, alors que Von Mises sera un scalaire énergétique.
🎯 Objectif Scientifique
Savoir que la contrainte maximale vaut 98.1 MPa est nécessaire mais pas suffisant. Pour un ingénieur structure, il est vital de savoir dans quelle direction cette force maximale s'exerce. Cela permet de positionner correctement les soudures, d'orienter les fibres de renfort ou de prédire la direction de propagation d'une éventuelle fissure (qui s'ouvrira perpendiculairement à \(\sigma_1\)). L'objectif est de trouver l'angle \(\alpha_{\text{p}}\) de rotation du repère.
📚 Référentiel Théorique
Cercle de MohrTrigonométrieSur le cercle de Mohr, les angles sont doublés. L'angle entre le point représentatif de l'état initial \(X\) et le point principal \(P_1\) (sur l'axe horizontal) est \(2\alpha_{\text{p}}\). Nous allons utiliser la géométrie du triangle rectangle formé par le point \(X(\sigma_x, -\tau_{xy})\), le centre \(C\) et la projection de \(X\) sur l'axe horizontal pour trouver cet angle double, puis en déduire l'angle physique réel.
Considérons le triangle rectangle formé par le centre \(C\), le point \(X\) et sa projection \(H\) sur l'axe horizontal. L'angle en \(C\) est \(2\alpha_{\text{p}}\).
- Le côté opposé est la hauteur de \(X\), soit \(\tau_{xy}\).
- Le côté adjacent est la distance \(CH = \sigma_x - \sigma_{\text{moy}} = \frac{\sigma_x - \sigma_y}{2}\).
- La tangente est donc :
Cette formule découle directement de la simplification de la fraction obtenue par le raisonnement géométrique.
📋 Données d'Entrée
| Paramètre | Valeur |
|---|---|
| Différence Normales (\(\sigma_x - \sigma_y\)) | 120 \(\text{ MPa}\) |
| Cisaillement (\(\tau_{xy}\)) | 50 \(\text{ MPa}\) |
Si le dénominateur est nul (\(\sigma_x = \sigma_y\)), l'angle double vaut 90°, donc l'angle principal vaut 45°. C'est le cas du cisaillement pur.
📝 Application Numérique Détaillée
Schéma : L'angle 2α représente la rotation sur le cercle pour aller de l'état X à l'état principal σ1.
1. Calcul de la tangente de l'angle double :
Injectons les valeurs initiales dans la formule.
2. Extraction de l'angle principal \(\alpha_{\text{p}}\) :
Nous utilisons la fonction arctangente pour trouver \(2\alpha_{\text{p}}\), puis nous divisons par 2.
Interprétation : Il faut tourner l'élément de matière de +19.9° dans le sens anti-horaire (sens trigonométrique) pour que la facette initialement verticale subisse la contrainte principale majeure \(\sigma_1\).
✅ Interprétation Globale de l'Étape 3
L'analyse directionnelle montre que la contrainte maximale ne s'exerce pas horizontalement (selon la membrure) mais est déviée d'environ 20° vers le haut à cause de l'action des diagonales (cisaillement positif).
L'angle est positif, ce qui est cohérent avec le signe positif du cisaillement \(\tau_{xy}\) qui "tire" le point X vers le haut du cercle, donc vers une rotation anti-horaire pour rejoindre l'axe principal.
Attention au quadrant ! Si le dénominateur \((\sigma_x - \sigma_y)\) était négatif, il faudrait ajouter 180° à l'angle double (fonction `atan2` en programmation). Ici, le dénominateur est positif, donc l'angle est direct.
🎯 Objectif Scientifique
C'est l'étape de décision finale. Les aciers ductiles comme le S355 ne rompent pas simplement quand une contrainte atteint un seuil, mais quand la densité d'énergie de distorsion atteint une valeur critique. Le critère de Von Mises est une formule énergétique qui combine les deux contraintes principales en une seule valeur scalaire, la contrainte équivalente \(\sigma_{\text{vm}}\). C'est cette valeur unique que l'on compare à la limite d'élasticité du matériau pour valider la tenue de la structure.
📚 Référentiel Théorique
Critère de Von MisesEurocode 3 (ELU)Nous avons un état de contrainte mixte (traction dans une direction, compression dans l'autre). Cette situation est souvent plus critique que la traction simple car les contraintes de cisaillement (responsables du glissement plastique) s'additionnent. La formule de Von Mises en contraintes principales pour un état plan (\(\sigma_3 = 0\)) est l'outil parfait pour quantifier cette sévérité. Si \(\sigma_{\text{vm}} < f_y\), la déformation reste élastique et réversible. Sinon, la plastification commence.
Le critère complet de Von Mises s'écrit :
Pour un état plan de contraintes, la troisième contrainte principale \(\sigma_3\) est nulle (surface libre). En remplaçant \(\sigma_3\) par 0 et en développant les carrés, on obtient la formule simplifiée ci-dessous.
Formulation simplifiée utilisant les contraintes principales planes.
📋 Rappel des Données Calculées
| Paramètre | Valeur |
|---|---|
| Principale 1 (\(\sigma_1\)) | 98.1 \(\text{ MPa}\) |
| Principale 2 (\(\sigma_2\)) | -58.1 \(\text{ MPa}\) |
| Limite Élastique (\(f_y\)) | 355 \(\text{ MPa}\) |
Dans un état de cisaillement pur, Von Mises vaut \(\sqrt{3} \times \tau\). Ici, avec un état mixte, le calcul complet est indispensable.
📝 Application Numérique Détaillée
1. Calcul de la Contrainte Équivalente \(\sigma_{\text{vm}}\) :
Nous remplaçons les termes dans la formule quadratique.
Interprétation : L'état de contrainte complexe équivaut, en termes de dangerosité plastique, à une traction simple de 136.7 MPa.
2. Calcul du Taux de Travail (Ratio) :
Nous comparons la contrainte équivalente à la capacité limite du matériau.
Interprétation : Le matériau n'est sollicité qu'à 38.5% de sa capacité élastique maximale.
✅ Conclusion Finale
La contrainte équivalente de Von Mises (136.7 MPa) est très largement inférieure à la limite d'élasticité de l'acier S355 (355 MPa). Le coefficient de sécurité réel est de \(1/0.385 \approx 2.6\), ce qui est bien supérieur aux coefficients de sécurité partiels habituels de l'Eurocode (1.0 ou 1.1). La pièce est donc validée.
La valeur de 136.7 est supérieure à la plus grande contrainte principale (98.1), ce qui est normal car la compression transverse "aide" au cisaillement et aggrave l'état de contrainte au sens de Von Mises.
Regardez attentivement le terme croisé \(- \sigma_1 \sigma_2\). Comme \(\sigma_2\) est négative (compression), le produit \(\sigma_1 \sigma_2\) est négatif. Avec le signe moins devant, ce terme devient positif (\(+5699.6\)). C'est ce mécanisme qui rend les états de contrainte "Traction + Compression" beaucoup plus sévères au sens de Von Mises que la simple somme des carrés. Ne jamais négliger ce signe !
📄 Livrable Final (Note de Calculs EXE)
| Ind. | Date | Objet de la modification | Rédacteur |
|---|---|---|---|
| A | 10/10/2024 | Hypothèses initiales | J. Dupont |
| B | 24/10/2024 | Ajout vérification Von Mises | A. Martin |
- Eurocode 3 : Calcul des structures en acier (EN 1993-1-1)
- Matériau : Acier de construction S355 J2
| Contrainte Normale X | + 80 \(\text{ MPa}\) |
| Contrainte Normale Y | - 40 \(\text{ MPa}\) |
| Cisaillement XY | + 50 \(\text{ MPa}\) |
Résultats issus de la diagonalisation analytique du tenseur des contraintes (Méthode de Mohr).
L'Ingénieur Stagiaire
Chef de Projet Structure
Laisser un commentaire