Critère de Rupture de Mohr-Coulomb
📝 Situation du Projet et Enjeux Structurels
Dans le cadre de l'aménagement critique d'une nouvelle infrastructure industrielle lourde, le bureau d'études géotechniques de pointe est expressément mandaté pour valider la portance absolue du sol d'assise. En effet, la future superstructure, lourdement composée de charpentes métalliques massives et de machines vibrantes, transmettra des charges colossales et permanentes au terrain naturel de fondation. Si ces sollicitations mécaniques viennent à dépasser la capacité intrinsèque de résistance au cisaillement du sol, nous risquons d'assister à un poinçonnement brutal, entraînant la ruine complète de l'ouvrage. La sécurité globale du projet repose donc intégralement sur notre modélisation mathématique rigoureuse du comportement mécanique du terrain.
Par conséquent, une vaste campagne de reconnaissance a été déclenchée. Des prélèvements intacts par carottage ont été réalisés in-situ à une profondeur exacte de \(4 \text{ mètres}\), correspondant à la base prévue des futures semelles de fondation. Ces échantillons cylindriques de sable argileux ont été précieusement rapatriés en laboratoire pour y subir une série d'essais triaxiaux de révolution. Il est impératif de souligner que ces tests destructifs ont été conduits sous des conditions expérimentales rigoureuses dites "Consolidées Drainées" (CD). C'est pourquoi, notre premier objectif scientifique et technique consiste à exploiter minutieusement les données de pressions issues de cette campagne d'essais afin d'en extraire les véritables paramètres intrinsèques de résistance de la matrice granulaire.
Dans un second temps, une fois ces limites physiques clairement identifiées, modélisées et quantifiées, il nous incombera de projeter ce modèle théorique sur le cas pratique de notre chantier. Nous devrons ainsi évaluer la stabilité réelle de l'ouvrage in-situ face aux colossales descentes de charges calculées par les ingénieurs structure, garantissant de cette façon la pérennité décennale de l'installation.
En tant qu'Ingénieur Géotechnicien Principal, vous avez l'impérieuse responsabilité de déterminer les caractéristiques mécaniques fondamentales du massif de sol (la cohésion effective et l'angle de frottement interne) en exploitant le célèbre critère de rupture de Mohr-Coulomb. Par la suite, vous êtes tenu d'appliquer ces paramètres fraîchement calculés pour vérifier méthodiquement si l'état de contrainte titanesque généré par la future fondation de l'usine franchit la limite de rupture par cisaillement.
"Prenez garde aux conditions expérimentales ! Les essais sont du type Consolidé Drainé (CD). Cela certifie que l'eau a eu le temps de s'échapper des pores du sol durant le chargement. En conséquence, la pression interstitielle est maintenue strictement nulle (\(u = 0\)) tout au long du processus jusqu'à la rupture. Les contraintes totales mécaniques appliquées par les pistons sont donc mathématiquement et physiquement égales aux contraintes effectives agissant sur le squelette solide (\(\sigma = \sigma'\)). Traitez rigoureusement vos unités en kilopascals (\(\text{kPa}\))."
L'ensemble des paramètres scientifiques exposés ci-dessous définit le socle expérimental de notre étude. Il est crucial de noter que les contraintes mécaniques relevées par les capteurs de la presse correspondent exclusivement à l'état de rupture critique. Il s'agit de la fraction de seconde exacte où le squelette granulaire de l'échantillon capitule et subit une déformation de cisaillement brutale et irréversible, marquant la naissance d'un plan de glissement interne.
📚 Référentiel Mécanique et Normatif
Norme NF EN ISO 17892-9 (Essai Triaxial)Théorie de Rupture de Mohr-CoulombAfin de cartographier rigoureusement l'enveloppe de résistance de notre terrain d'assise, le laboratoire géotechnique a méthodiquement détruit deux carottes de sol géologiquement identiques. Il est fondamental de comprendre que ces tests visent à reproduire artificiellement différentes conditions d'enfouissement pour mesurer la résistance maximale de la terre avant sa dislocation absolue.
Dans un premier temps, l'Essai Triaxial N°1 a été initié pour simuler un enfouissement relativement superficiel. L'éprouvette a été placée dans la chambre pressurisée sous une pression d'eau radiale constante de \(50 \text{ kPa}\). Ce confinement, modélisant mathématiquement la contrainte principale mineure (\(\sigma'_{3,1}\)), maintient les grains serrés et reproduit fidèlement l'étreinte naturelle du sol. Par la suite, le piston d'acier vertical a écrasé l'échantillon à une vitesse millimétrique. L'effondrement structurel complet, caractérisé par un plan de cisaillement net, a été enregistré par les capteurs au moment précis où la contrainte verticale effective majeure (\(\sigma'_{1,1}\)) atteignait l'intensité critique de \(300 \text{ kPa}\).
Dans un second temps, pour comprendre l'évolution du frottement interne face à une contrainte de profondeur plus importante, l'Essai Triaxial N°2 a été mené sur une éprouvette jumelle. Le fluide de la cellule vitrée a cette fois appliqué une étreinte latérale féroce de \(150 \text{ kPa}\) en guise de contrainte mineure (\(\sigma'_{3,2}\)). En effet, grâce à ce confinement supérieur qui empêche littéralement les grains de glisser les uns sur les autres, le piston vertical a dû développer un effort mécanique bien plus monumental avant de vaincre la cohésion de la matrice. La rupture de ce second échantillon s'est finalement déclenchée sous une contrainte axiale écrasante, soit une contrainte majeure (\(\sigma'_{1,2}\)) culminant à \(700 \text{ kPa}\).
L'implantation de la future superstructure industrielle sur le terrain naturel va irrémédiablement bouleverser l'équilibre géostatique originel. Notre démarche d'ingénieur impose de confronter les capacités intrinsèques du sol, qui seront déduites des essais en laboratoire, à l'agression réelle qu'il subira inévitablement sur le chantier.
D'une part, focalisons-nous sur l'interface précise de contact entre la base de la semelle en béton et le sol d'assise, située à \(4 \text{ mètres}\) de profondeur. À ce niveau stratigraphique, le poids mort des couches de terres sus-jacentes exerce une pression latérale continue. L'analyse détaillée du tenseur géostatique naturel, fournie par la campagne de sondages préalables, révèle que la contrainte effective de confinement in-situ projetée (\(\sigma'_{3,\text{in-situ}}\)) s'établit très exactement à \(100 \text{ kPa}\). Cette pression latérale agit inéluctablement comme un bouclier et un renfort naturel indispensable pour la matrice granulaire face aux charges futures.
D'autre part, intéressons-nous à l'action destructrice potentielle de l'ouvrage. Le bâtiment lui-même va injecter une énergie descendante colossale dans le sous-sol. En agrégeant méticuleusement le poids propre des voiles de béton armé, la lourde charpente de toiture et les charges d'exploitation cycliques des machines industrielles, la note de calculs structurelle nous a formellement transmis le pic de sollicitation verticale. L'effort extrême induit par la fondation frappera le sol sous une contrainte principale effective majeure (\(\sigma'_{1,\text{bâtiment}}\)) estimée à \(450 \text{ kPa}\).
| Description Physique | Variable Mathématique | Mesure Relevée | Unité Internationale |
|---|---|---|---|
| Étreinte hydrostatique radiale (Essai 1) | \(\sigma'_{3,1}\) | \(50\) | \(\text{kPa}\) |
| Charge d'écrasement vertical à la ruine (Essai 1) | \(\sigma'_{1,1}\) | \(300\) | \(\text{kPa}\) |
| Étreinte hydrostatique radiale (Essai 2) | \(\sigma'_{3,2}\) | \(150\) | \(\text{kPa}\) |
| Charge d'écrasement vertical à la ruine (Essai 2) | \(\sigma'_{1,2}\) | \(700\) | \(\text{kPa}\) |
| Renfort latéral géostatique (Chantier) | \(\sigma'_{3,\text{in-situ}}\) | \(100\) | \(\text{kPa}\) |
| Agression verticale transmise par l'usine | \(\sigma'_{1,\text{bâtiment}}\) | \(450\) | \(\text{kPa}\) |
E. Protocole de Résolution
Voici la méthodologie séquentielle recommandée pour mener à bien cette étude géotechnique complète, allant de l'interprétation des données brutes jusqu'à la vérification structurelle.
Analyse Géométrique des Cercles de Mohr
Calculer les rayons et les coordonnées des centres pour chaque état de contrainte à la rupture issu des essais triaxiaux.
Détermination du Critère de Mohr-Coulomb
Exploiter l'enveloppe de rupture tangente aux deux cercles pour en déduire l'angle de frottement interne (\(\varphi'\)) et la cohésion effective (\(c'\)).
Calcul du Plan de Rupture
Déterminer mathématiquement l'inclinaison théorique (\(\theta\)) du plan de glissement qui se forme à l'intérieur de l'échantillon lors de la rupture.
Vérification de la Stabilité In-Situ
Utiliser les paramètres intrinsèques obtenus pour évaluer la sécurité de la fondation sous les charges réelles du futur bâtiment.
Critère de Rupture de Mohr-Coulomb
🎯 Objectif
Dans cette première étape analytique, notre but absolu est de transformer les données de pressions brutes issues des pistons de la presse triaxiale en grandeurs géométriques exploitables mathématiquement. En effet, le plan de Mohr (constitué des axes de contrainte normale et de cisaillement) est l'espace de référence universel dans lequel le critère de rupture s'exprime graphiquement par une droite tangente. C'est pourquoi nous devons d'abord déterminer avec une extrême précision la position sur l'axe des abscisses (le centre) et la taille (le rayon) des deux cercles de rupture représentant la fin de vie de nos deux échantillons.
📚 Référentiel
Représentation Graphique de Mohr des Contraintes 2DAvant de foncer tête baissée dans les équations de résistance, prenons le temps d'observer notre jeu de données expérimental et de comprendre sa traduction géométrique. Chaque essai triaxial nous fournit les contraintes principales majeure (\(\sigma'_1\)) et mineure (\(\sigma'_3\)) à l'instant exact et fatal de la rupture. Sur un axe unidimensionnel, ces deux valeurs constituent les deux extrémités (les deux pôles) du cercle de Mohr. Par conséquent, la distance totale séparant ces deux points représente rigoureusement le diamètre du cercle. En divisant ce diamètre par deux, nous obtiendrons le rayon. Ensuite, en ajoutant ce rayon à la valeur de départ (\(\sigma'_3\)), nous tomberons inexorablement sur la coordonnée du centre du cercle.
La démonstration est purement arithmétique. Le diamètre (\(D\)) est la différence entre la contrainte maximale et la contrainte minimale : \(D = \sigma'_1 - \sigma'_3\). Le rayon (\(R\)) est par définition la moitié de ce diamètre, d'où \(R = (\sigma'_1 - \sigma'_3)/2\). Pour localiser le centre (\(C\)) situé à mi-chemin, on peut calculer la moyenne arithmétique des deux extrêmes. La formule se dérive ainsi : \(C = \sigma'_3 + R = \sigma'_3 + (\sigma'_1 - \sigma'_3)/2\). En mettant sous un dénominateur commun, on obtient \(C = (2\sigma'_3 + \sigma'_1 - \sigma'_3)/2\), ce qui se simplifie magnifiquement en \(C = (\sigma'_1 + \sigma'_3)/2\).
📋 Données d'Entrée
| Paramètre Analytique | Valeur Lues (Laboratoire) |
|---|---|
| Contrainte Mineure (Essai 1) - \(\sigma'_{3,1}\) | \(50 \text{ kPa}\) |
| Contrainte Majeure (Essai 1) - \(\sigma'_{1,1}\) | \(300 \text{ kPa}\) |
| Contrainte Mineure (Essai 2) - \(\sigma'_{3,2}\) | \(150 \text{ kPa}\) |
| Contrainte Majeure (Essai 2) - \(\sigma'_{1,2}\) | \(700 \text{ kPa}\) |
Une erreur grossière et récurrente chez les jeunes concepteurs consiste à confondre la contrainte déviatorique brute (notée \(q = \sigma'_1 - \sigma'_3\)), qui est lue directement sur la presse, avec le rayon du cercle. N'oubliez jamais que le rayon de Mohr est obligatoirement le demi-déviateur (\(q/2\)) ! Diviser par deux est une étape non négociable.
📝 Calcul Détaillé
Nous appliquons à présent de manière très systématique nos équations démontrées sur les données brutes de chaque essai pour en extraire l'essence géométrique.
A. Opérations sur l'Essai N°1 (Faible Confinement)
Nous isolons la somme et la différence des pressions de \(300 \text{ kPa}\) et \(50 \text{ kPa}\) :
Exécution du Calcul du Centre \(C_1\) :La pression isostatique moyenne subie lors de la première rupture culmine donc à \(175 \text{ kPa}\).
Exécution du Calcul du Rayon \(R_1\) :Le cisaillement destructeur subi par la première éprouvette a atteint un pic mortel de \(125 \text{ kPa}\).
B. Opérations sur l'Essai N°2 (Fort Confinement)
Nous reproduisons la même démarche logique avec les contraintes majorées de \(700 \text{ kPa}\) et \(150 \text{ kPa}\) :
Exécution du Calcul du Centre \(C_2\) :✅ Interprétation Globale
Nous constatons de manière indiscutable que lorsque la pression de confinement (\(\sigma'_3\)) augmente, le centre du cercle se translate puissamment vers la droite, et surtout, le rayon gonfle de manière très significative. En conséquence, cette observation empirique prouve que notre sable argileux voit sa résistance au cisaillement croître prodigieusement sous l'effet de la pression. C'est l'illustration parfaite du comportement frottant d'un milieu granulaire : plus on l'écrase, plus il devient dur à cisailler.
Un contrôle de rétro-ingénierie rapide est toujours bienvenu. Nous savons que le sommet droit du cercle (\(\sigma'_1\)) doit être égal au centre plus le rayon (\(C + R\)). Vérifions pour le second essai : \(425 + 275 = 700 \text{ kPa}\). Ce résultat retombe exactement sur notre contrainte de rupture axiale lue sur la machine. L'architecture de nos cercles est mathématiquement irréprochable.
Assurez-vous absolument que toutes les contraintes intégrées dans ces équations sont des contraintes effectives (notées avec un repère prime). Si les essais avaient été menés en condition non-drainée (CU ou UU) sans mesure de la pression interstitielle \(u\), traiter ces valeurs totales conduirait à un cercle de Mohr erroné, décalé vers la droite, générant un faux sentiment de sécurité.
🎯 Objectif
Maintenant que nos deux cercles de rupture sont mathématiquement localisés et dimensionnés, notre mission critique est d'en extraire l'équation paramétrique de l'enveloppe de rupture. Cette droite fuyante, parfaitement tangente aux parties supérieures des cercles, matérialise la limite physique absolue, le plafond de verre infranchissable pour le matériau. En effet, l'identification de son équation affine nous livrera les deux paramètres intrinsèques ultimes du sol : l'angle de frottement interne effectif (\(\varphi'\)) et la cohésion géotechnique (\(c'\)).
📚 Référentiel
Critère de Ruine de Mohr-Coulomb (Mécanique des Milieux Continus)Comment déduire l'angle de cette droite sans la dessiner ? Il faut visualiser l'agencement géométrique. Relions les centres \(C_1\) et \(C_2\) de nos deux cercles aux points de tangence sur l'enveloppe. Par définition, un rayon touchant un point de tangence forme un angle droit (\(90^\circ\)) avec la droite. Nous sommes donc en présence de triangles rectangles. Par conséquent, si nous traçons une ligne horizontale fictive partant de \(C_1\) et s'arrêtant sous le centre \(C_2\), nous découpons un immense triangle rectangle. L'hypoténuse de ce triangle est la distance entre les centres (\(C_2 - C_1\)), et son côté opposé à l'angle de l'enveloppe correspond à la différence d'élévation des points de tangence, soit la différence des rayons (\(R_2 - R_1\)). Le sinus de l'angle est ainsi mis en évidence.
Une fois l'angle \(\varphi'\) défini, il faut remonter à la cohésion (\(c'\)), qui est l'ordonnée à l'origine (\(\tau\) lorsque \(\sigma' = 0\)). Si l'on prolonge l'enveloppe de rupture vers la gauche jusqu'à croiser l'axe horizontal, elle frappe l'abscisse à une distance négative égale à \(-c' \times \cot(\varphi')\). Si nous formons un triangle rectangle partant de cette "origine virtuelle" jusqu'au centre \(C_1\), l'hypoténuse globale posée sur l'axe des abscisses mesure \((c' \cot(\varphi') + C_1)\). Le côté opposé est le rayon \(R_1\). Dès lors, la loi de trigonométrie impose que : \(\sin(\varphi') = \frac{R_1}{C_1 + c' \cot(\varphi')}\). En développant algébriquement cette égalité (\(C_1 \sin(\varphi') + c' \cos(\varphi') = R_1\)), on isole très facilement la variable \(c'\) pour obtenir la formule analytique de projection.
📋 Données d'Entrée
| Coordonnée Isostatique (Centre) | Contrainte Déviatorique Maximale (Rayon) |
|---|---|
| Cercle 1 : \(C_1 = 175 \text{ kPa}\) | Cercle 1 : \(R_1 = 125 \text{ kPa}\) |
| Cercle 2 : \(C_2 = 425 \text{ kPa}\) | Cercle 2 : \(R_2 = 275 \text{ kPa}\) |
Les fonctions trigonométriques inverses (\(\arcsin\)) sont extrêmement sensibles aux arrondis. Je vous conseille vivement de conserver au minimum quatre décimales, voire la fraction exacte (\(3/5\)), lors des calculs intermédiaires de \(\sin\) et \(\cos\), pour éviter de propager une grave erreur d'arrondi sur le calcul final de la cohésion qui suit.
📝 Calcul Détaillé
Nous allons dérouler pas-à-pas la relation des pentes différentielles pour verrouiller l'angle, puis isoler algébriquement la cohésion structurelle.
A. Détermination de l'Angle de Frottement Interne (\(\varphi'\))
Nous substituons les grandeurs de nos deux cercles dans le ratio différentiel de Thalès pour isoler le rapport trigonométrique :
Connaissant cette proportion exacte, nous déclenchons la fonction arc sinus inverse pour révéler l'angle physique :
Extraction de l'angle réel :B. Détermination de la Cohésion Effective (\(c'\))
Avant de substituer dans la formule de projection, complétons notre arsenal trigonométrique. Nous savons que \(\sin(\varphi') = 0.6\). Par la loi fondamentale \(\cos^2 + \sin^2 = 1\), on déduit que \(\cos(\varphi') = \sqrt{1 - 0.6^2} = \sqrt{0.64} = 0.8\).
Pour prouver l'homogénéité de notre modèle et confirmer que la cohésion est bien une caractéristique intrinsèque (invariable selon le confinement), nous allons projeter l'équation sur nos deux états de rupture consécutifs :
Projection sur l'Essai N°1 (Cercle 1) :Bilan de la projection : Nous obtenons une égalité parfaite (\(c'_1 = c'_2\)). La démonstration est totale, la cohésion de la matrice granulaire est scellée à \(25 \text{ kPa}\) de manière universelle pour cette couche de sol.
✅ Interprétation Globale
Les paramètres intrinsèques qui conditionneront tout le reste de notre ingénierie sont désormais scellés. L'angle de frottement interne élevé de \(36.87^\circ\) dénote un matériau très granulaire et anguleux, bénéficiant d'un excellent enchevêtrement mécanique sous charge. Par ailleurs, la cohésion effective non négligeable de \(25 \text{ kPa}\) vient confirmer le constat visuel : il s'agit d'un sable "argileux". La matrice d'argile agit comme un liant, une légère colle interne, qui offre au sol une petite résistance structurelle même lorsqu'il est complètement à l'air libre, sans aucun confinement.
Faisons appel à notre bon sens géotechnique. Les valeurs trouvées (\(\varphi'=37^\circ\) et \(c'=25 \text{ kPa}\)) sont-elles réalistes pour le sol décrit ? Absolument. Un sable propre, lavé de toute particule fine, aurait affiché une cohésion strictement nulle (\(c'=0\)). À l'inverse, une argile molle pure aurait présenté un angle de frottement rachitique (souvent sous les \(20^\circ\)). Nos résultats hybrides cadrent magnifiquement avec la classification de sable argileux dense.
Si la soustraction mathématique de la cohésion vous avait renvoyé une valeur négative (ex: \(-10 \text{ kPa}\)), les alarmes auraient dû retentir. Dans la nature, un sol ne peut physiologiquement pas posséder une cohésion négative (ce qui signifierait que les grains se repoussent magnétiquement). Un tel résultat mettrait en lumière une aberration expérimentale majeure (bulle d'air coincée, fuite membranaire lors de l'essai) ou une violation de l'hypothèse de linéarité de l'enveloppe de Mohr-Coulomb aux très fortes pressions.
🎯 Objectif
Lorsqu'un massif de sol atteint son état de stress critique sous le poids d'un ouvrage, la rupture géologique ne se produit pas sous la forme d'une explosion chaotique. Au contraire, elle s'organise mathématiquement et se propage le long d'une ligne de faiblesse bien spécifique, que l'on nomme pompeusement le plan de glissement. Notre objectif ciblé ici est de prédire par le calcul algébrique l'inclinaison spatiale théorique de ce plan cisaillant, notée \(\theta\), par rapport au plan principal majeur (qui est l'horizontale, sous la semelle, dans notre cas de fondation).
📚 Référentiel
Théorie du Pôle de Mohr et Cinématique de RankinePour lier le cercle bidimensionnel à la géométrie réelle en 3D, nous utilisons la règle de transformation des plans. Le point critique de rupture sur le cercle correspond au point de tangence exact avec l'enveloppe. Géométriquement, l'angle formé au centre du cercle entre l'axe horizontal et le rayon pointant vers ce point de tangence vaut toujours \(90^\circ + \varphi'\) (puisque la tangente est perpendiculaire au rayon). Cependant, la mécanique des milieux continus démontre de manière implacable que tout angle balayé au centre du cercle de Mohr correspond au double de l'angle physique réel d'inclinaison du plan de contrainte dans le massif de sol. L'équation de base est donc \(2\theta = 90^\circ + \varphi'\). Il suffit de diviser par deux pour clore le problème.
Contrairement aux matériaux hautement ductiles ou parfaitement cohésifs (comme l'acier doux de construction ou l'argile non drainée) qui vont cisailler préférentiellement et systématiquement à un angle neutre de \(45^\circ\) sous compression axiale, les sols sont massivement dotés de frottement interne. Cette rugosité particulaire "dévie" mathématiquement la trajectoire du plan critique. Le frottement force le plan à s'incliner de manière nettement plus abrupte d'une valeur strictement équivalente à la moitié de l'angle de frottement intrinsèque \(\varphi'\).
📋 Données d'Entrée
| Paramètre Analytique Nécessaire | Valeur Intégrée |
|---|---|
| Frottement interne effectif verrouillé (\(\varphi'\)) | \(36.87^\circ\) |
Faites très attention à votre référentiel topographique ! L'angle \(\theta\) calculé par cette formule est toujours l'angle formé par rapport au plan principal majeur (le plan où agit \(\sigma'_1\)). Dans le cas d'une fondation verticale, \(\sigma'_1\) écrase le sol verticalement. Le plan sur lequel s'applique cette charge est donc logiquement le sol horizontal. L'angle \(\theta\) est par conséquent mesuré directement depuis la ligne d'horizon.
📝 Calcul Détaillé
L'opération s'effectue par une division élémentaire et une addition séquentielle dans l'équation cinématique. Nous prenons soin de traiter la décimale angulaire avec respect.
A. Substitution et Résolution Angulaire
Nous incorporons la valeur caractérisée de \(36.87^\circ\) dans la demi-somme de Rankine :
✅ Interprétation Globale
La conclusion visuelle est sans appel : la rupture ne cisaillera absolument pas le sous-sol à un angle modéré de \(45^\circ\). En réalité, la surface de glissement destructrice traversera le bulbe de pression de la fondation sous un angle abrupt et plongeant de \(63.4^\circ\). C'est une donnée spatiale d'une valeur inestimable, car elle conditionnera directement le dimensionnement de la largeur de sécurité des radiers ou la profondeur d'ancrage nécessaire pour des tirants de soutènement éventuels visant à enjamber ce coin de glissement fatal.
Vérification expresse : Pour un sol pulvérulent classique ou frottant, l'angle \(\theta\) est obligatoirement compris entre \(45^\circ\) (cas d'un fluide parfait ou argile pure saturée, \(\varphi'=0\)) et \(90^\circ\) (matériau infiniment résistant). Notre score de \(63.4^\circ\) tombe confortablement dans le spectre réaliste et typique des sables fortement compactés.
Le piège absolu guette les distraits : de nombreux techniciens omettent de diviser l'angle \(\varphi'\) par deux dans la précipitation algébrique, ou confondent l'angle du pôle sur le cercle de Mohr (qui vaudrait ici \(90^\circ+37^\circ=127^\circ\)) avec l'angle de déclinaison physique de la faille dans le terrain. Ce facteur de division par \(2\) entre le cercle et la réalité est la traduction mathématique primordiale de la mécanique de Mohr.
🎯 Objectif
Le triomphe théorique réalisé dans les étapes précédentes ne suffit en aucun cas à rassurer le maître d'ouvrage. Il est grand temps d'appliquer notre modèle mathématique à l'agressivité du monde réel du chantier. Notre colossale infrastructure industrielle génère un état de contrainte massif sur le terrain d'assise originel. Notre ultime mission d'expert est de vérifier, avec une certitude mathématique absolue et indéniable, si la semelle de fondation provoquera le cisaillement mortel et l'affaissement du terrain environnant. Pour ce faire, nous allons calculer la "contrainte limite admissible" (le plafond de rupture absolu sous le bâtiment), pour ensuite la confronter de front à l'impact destructeur réel via un Facteur de Sécurité global.
📚 Référentiel
Théorie des États d'Équilibre Limite et Facteurs de Sécurité (Eurocode 7)Plutôt que de tracer un énième cercle de Mohr in-situ pour y chercher visuellement une collision avec l'enveloppe, optons pour l'élégance et la robustesse de l'algèbre. Nous pouvons dériver l'équation universelle d'état limite à partir de la relation géométrique déjà démontrée : \(\sin(\varphi') = \frac{\sigma'_1 - \sigma'_3}{\sigma'_1 + \sigma'_3 + 2c' \cot(\varphi')}\). En effectuant un produit en croix pour supprimer le dénominateur, on obtient : \(\sigma'_1 \sin(\varphi') + \sigma'_3 \sin(\varphi') + 2c' \cos(\varphi') = \sigma'_1 - \sigma'_3\). Ensuite, en factorisant les termes en \(\sigma'_1\) d'un côté de l'égalité, nous arrivons à \(\sigma'_1(1 - \sin(\varphi')) = \sigma'_3(1 + \sin(\varphi')) + 2c' \cos(\varphi')\). En divisant le tout par \((1 - \sin(\varphi'))\), on isole \(\sigma'_1\). Les puissantes identités trigonométriques nous permettent alors de remplacer le bloc \(\frac{1 + \sin(\varphi')}{1 - \sin(\varphi')}\) par la fonction simplifiée \(N_\varphi = \tan^2(45^\circ + \varphi'/2)\). L'équation devient une simple ligne mathématique prête à l'emploi pour calculer la limite de rupture !
L'équation majestueuse issue de notre démonstration repose sur le fameux facteur d'influence géométrique \(N_\varphi\). Dans la mécanique des sols, ce terme agit littéralement comme un multiplicateur exponentiel de résistance. Il quantifie la façon dont l'angle de frottement interne des grains parvient à décupler et majorer l'influence positive du confinement naturel (\(\sigma'_3\)) pour s'opposer de plus en plus fermement à l'écrasement vertical.
📋 Données d'Entrée Restituées
| Données Environnementales du Chantier | Valeurs Imposées |
|---|---|
| Contrainte horizontale naturelle de confinement (\(\sigma'_{3,\text{in-situ}}\)) | \(100 \text{ kPa}\) |
| Pression verticale destructrice projetée par l'usine (\(\sigma'_{1,\text{bâtiment}}\)) | \(450 \text{ kPa}\) |
| Couples Paramétriques Validés Précédemment | \(\varphi' = 36.87^\circ\) | \(c' = 25 \text{ kPa}\) |
L'erreur fatale à ce stade décisionnel serait de réutiliser aveuglément les valeurs de confinement de \(50 \text{ kPa}\) ou \(150 \text{ kPa}\) issues des tests triaxiaux en laboratoire ! Ces valeurs de piston n'étaient là que pour calibrer la machine. Pour la vérification finale de conception in-situ, vous devez impérativement utiliser le confinement géostatique réel existant sous la future semelle, soit les \(100 \text{ kPa}\) fournis par les sondages, qui représentent la véritable ceinture de protection du sol naturel.
📝 Calcul Détaillé
Nous enclenchons la procédure de vérification en évaluant la composante trigonométrique \(N_\varphi\), pour faciliter grandement le déploiement de l'équation globale du sol étape par étape.
A. Résolution Numérique du Coefficient \(N_\varphi\)
Rappelons le résultat du calcul d'inclinaison : l'angle dans la tangente est \(45^\circ + 36.87^\circ/2 = 63.435^\circ\). La tangente de cet angle vaut rigoureusement \(2.0\). Élevons immédiatement ce terme au carré :
Le multiplicateur de butée du terrain est fixé à la valeur nette et sans unité de \(4.0\).
B. Détermination Analytique de la Capacité Limite Admissible
Nous injectons méthodiquement ce coefficient d'amplification dans la loi magistrale d'état limite, associé au confinement réel de chantier (\(100 \text{ kPa}\)) et à l'effet de colle argileux (\(25 \text{ kPa}\)) :
Ce plafond vertigineux de \(500 \text{ kPa}\) représente la frontière de ruine. Le mur de briques absolu de la matrice terrestre in-situ.
C. Calcul Sanction du Facteur de Sécurité Global
Il est temps de confronter la capacité de survie ultime du terrain (\(500 \text{ kPa}\)) face à la descente de charge implacable prévue par le plan d'architecture de l'usine (\(450 \text{ kPa}\)) :
✅ Interprétation Globale du Coefficient
Puisque le ratio mathématique aboutit à un facteur de sécurité strictement supérieur à la barre fatidique d'équilibre de \(1.0\) (\(1.11 > 1.0\)), la mécanique garantit que la fondation ne pulvérisera pas le terrain instantanément le jour de l'inauguration. L'état de contrainte titanesque imposé de force par l'édifice s'insère de justesse, sur le fil du rasoir, à l'intérieur du domaine de non-rupture délimité par la droite intransigeante de Mohr-Coulomb.
Les ordres de grandeur des pressions calculées (\(500 \text{ kPa}\) de résistance vs \(450 \text{ kPa}\) de charge) valident le réalisme de la substitution algébrique. Par ailleurs, le facteur \(F_s\), représentant la division de deux contraintes de même unité (\(\text{kPa} / \text{kPa}\)), se révèle être une grandeur parfaitement adimensionnelle, ce qui confirme l'homogénéité physique implacable de notre processus de vérification d'ingénierie.
Rangez vos sourires, la catastrophe est proche ! Bien que l'équation d'équilibre tienne encore mathématiquement debout, un ingénieur assermenté du Génie Civil ou un Bureau de Contrôle Technique ne signera absolument JAMAIS une validation structurelle avec un facteur de sécurité famélique de \(1.11\) ! Face aux inévitables hétérogénéités cachées du sous-sol, la pratique normative stricte (dictée par l'Eurocode 7 et les DTU fondations) impose que le coefficient de sécurité sur la portance superficielle soit au strict minimum de \(3.0\). Par conséquent, le projet d'infrastructure, dans sa conception actuelle, est techniquement irrecevable et virtuellement mort-né. Une refonte intégrale du système de transmission des charges (par élargissement majeur des radiers ou par ancrage profond sur pieux) est formellement requise avant tout coulage de béton !
📄 Livrable Final (Note de Calculs EXE)
| Ind. | Date | Objet de la modification | Rédacteur |
|---|---|---|---|
| A | 20/10/2023 | Intégration des données triaxiales de laboratoire | Bureau Sols |
| B | 24/10/2023 | Calcul analytique et alerte de sécurité structurelle | Ingénieur Principal |
L'interprétation de la tangente de rupture des cercles de Mohr (Méthode CD) a abouti aux paramètres du massif :
| Paramètre de Résistance | Valeur Approuvée |
|---|---|
| Angle de frottement interne effectif (\(\varphi'\)) | \(36.87^\circ\) |
| Cohésion apparente effective (\(c'\)) | \(25.00 \text{ kPa}\) |
| Coefficient de portance / butée (\(N_\varphi\)) | \(4.00\) |
Vérification rigoureuse de l'intégrité du sol sous l'application des charges d'infrastructure.
Ingénieur Géotechnicien (G2)
Directeur Bureau d'Études Sols
Laisser un commentaire