Calcul de la Contrainte induite dans le sol
📝 Situation du Projet
Vous avez intégré le bureau d'études techniques "GeoTech Engineering" en tant qu'ingénieur spécialiste des interactions sol-structure. Nous travaillons actuellement sur l'extension majeure du site industriel "Logistique Sud", situé dans une zone sédimentaire complexe de la vallée du Rhône. Le projet prévoit la construction d'un nouveau hall de production pour l'industrie automobile.
L'élément critique de ce projet est l'implantation d'une presse hydraulique de 2000 tonnes (capacité nominale). Cette machine, véritable cœur de l'usine, possède des tolérances de déformation extrêmement faibles : le moindre tassement différentiel pourrait désaligner les vérins et stopper la production. La structure porteuse repose sur des semelles superficielles isolées en béton armé.
Les sondages géotechniques (G2 AVP) ont révélé une stratigraphie délicate : sous une couche superficielle de sables compacts et raides (bon sol) se cache, à partir de 5 mètres de profondeur, une lentille d'argile molle normalement consolidée. Cette couche est très sensible au tassement à long terme. Votre mission est de vérifier si la contrainte apportée par la fondation se diffuse suffisamment avant d'atteindre cette couche critique.
Vous devez modéliser la transmission des charges dans le sol pour quantifier l'incrément de contrainte verticale (\(\Delta \sigma_z\)) au toit de la couche d'argile. C'est cette valeur précise qui servira ensuite au calcul des tassements de consolidation. Vous utiliserez la méthode analytique de diffusion des contraintes (théorie de Boussinesq adaptée par Steinbrenner pour les charges rectangulaires).
"Attention, ne confondez pas la charge Q (Force) avec la contrainte q (Pression). Le calcul doit se faire à la verticale exacte du centre géométrique de la semelle. Vérifiez bien la profondeur z par rapport à la sous-face de la semelle, pas par rapport au terrain naturel."
Pour mener à bien cette vérification, nous nous appuyons sur un ensemble de données validées par la maîtrise d'œuvre. Ces paramètres sont figés et servent de base contractuelle à votre note de calcul.
📚 Cadre Normatif et Réglementaire
Les calculs doivent être menés conformément aux standards européens et français en vigueur pour garantir l'assurabilité de l'ouvrage.
Eurocode 7 (EN 1997-1) : Calcul Géotechnique DTU 13.1 : Fondations SuperficiellesLa fondation étudiée est une semelle isolée rigide sous un poteau central. Les dimensions ont été pré-dimensionnées pour satisfaire les critères de portance (ELU), mais le critère de tassement (ELS) reste à vérifier.
Concernant le chargement, nous considérons ici uniquement les charges à l'État Limite de Service (ELS), car le tassement est un phénomène lié aux déformations de service et non à la rupture. Cette charge inclut le poids propre de la machine, de la structure et les charges d'exploitation.
| GÉOMÉTRIE DE LA FONDATION (S1) | |
| Largeur de la semelle (\(B\)) | 2.40 m |
| Longueur de la semelle (\(L\)) | 3.60 m |
| Profondeur de la couche d'argile (\(z\)) | 5.00 m (sous la base de la semelle) |
| CHARGEMENT APPLIQUÉ (ELS) | |
| Charge Verticale de Service (\(Q_{\text{ser}}\)) | 2160 kN (216 Tonnes) |
| Type de sollicitation | Charge centrée verticale |
| Répartition supposée | Uniforme (Semelle Rigide) |
| Donnée | Symbole | Valeur | Unité |
|---|---|---|---|
| Charge Verticale | \(Q\) | 2160 | kN |
| Largeur Semelle | \(B\) | 2.40 | m |
| Longueur Semelle | \(L\) | 3.60 | m |
| Profondeur de calcul | \(z\) | 5.00 | m |
E. Protocole de Résolution
Pour déterminer rigoureusement l'impact de la fondation sur les couches profondes, nous appliquerons la méthodologie séquentielle suivante :
Calcul de la Contrainte de Contact (\(q\))
Transformation de la force concentrée (kN) en une pression surfacique uniformément répartie (kPa) sous la semelle.
Détermination des Paramètres Géométriques (\(m, n\))
Calcul des ratios adimensionnels normalisés liant les dimensions de la semelle à la profondeur d'investigation.
Calcul du Coefficient d'Influence (\(I_z\))
Application de la solution de Steinbrenner (intégrale de Boussinesq) pour évaluer la dissipation de la contrainte.
Calcul Final de la Contrainte Induite (\(\Delta \sigma_z\))
Détermination de la valeur de la surcharge verticale reçue par l'argile à 5 mètres de profondeur.
Calcul de la Contrainte induite dans le sol
🎯 Objectif
L'objectif de cette première étape est de convertir la charge ponctuelle (ou résultante) \(Q_{\text{ser}}\) transmise par le poteau de la structure en une action surfacique répartie sur le sol d'assise. En géotechnique, le sol ne réagit pas directement à une force globale en Newtons, mais à une pression locale exprimée en Pascals (Pa) ou Kilopascals (kPa). C'est cette pression de contact à l'interface béton/sol qui constitue la "source" de contrainte qui va ensuite se propager et s'atténuer en profondeur.
📚 Référentiel
- Mécanique des Milieux Continus : Définition du tenseur des contraintes.
- DTU 13.1 : Hypothèses de répartition des contraintes sous les semelles superficielles.
Avant de lancer le calcul, nous devons poser une hypothèse sur la rigidité de la semelle. Une semelle en béton armé de fondation est généralement considérée comme infiniment rigide par rapport au sol (qui est souple). Cela a une conséquence majeure sur la distribution des contraintes : pour un chargement centré, la semelle s'enfonce uniformément dans le sol (tassement constant), ce qui implique théoriquement une réaction du sol uniforme.
Bien que la réalité physique soit plus complexe (concentration des contraintes sur les bords pour une argile, ou au centre pour un sable), l'usage courant en ingénierie pour les calculs de tassement (ELS) est d'admettre une répartition rectangulaire uniforme de la pression sous la semelle (\(q = \text{Cte}\)). C'est cette hypothèse simplificatrice mais robuste que nous adoptons ici.
La contrainte normale moyenne \(\sigma\) (notée ici \(q\)) est définie physiquement comme le ratio d'une force \(F\) agissant perpendiculairement sur une surface \(S\).
Dans le système international, une force de 1 Newton sur 1 mètre carré équivaut à 1 Pascal. En géotechnique, les ordres de grandeur sont souvent en kilonewtons (kN) et mètres carrés (m²), donnant des résultats en kilopascals (kPa). \(1 \text{ kPa} = 1 \text{ kN/m}^2 = 1000 \text{ Pa}\).
📋 Données d'Entrée
| Paramètre | Symbole | Valeur |
|---|---|---|
| Charge Verticale | \(Q_{\text{ser}}\) | 2160 kN |
| Largeur Semelle | \(B\) | 2.40 m |
| Longueur Semelle | \(L\) | 3.60 m |
Vérifiez toujours vos unités avant de diviser. Il est fréquent de recevoir des charges en Tonnes ou en Méganewtons (MN). Convertissez tout systématiquement en kN et en mètres pour obtenir directement des kPa, l'unité standard des rapports de sol.
Calculs Détaillés
1. Calcul de la Surface de la Semelle (\(A\)) :
Nous calculons d'abord l'aire géométrique de l'interface béton-sol en multipliant la largeur par la longueur.
La semelle offre une surface d'appui de 8.64 m².
2. Calcul de la Pression de Contact (\(q\)) :
Nous divisons maintenant la charge totale par cette surface pour obtenir la densité de force.
La contrainte moyenne appliquée au sol superficiel est exactement de 250 kPa.
✅ Interprétation Globale
Le sol directement sous la semelle subit une pression de 250 kN/m². C'est une valeur conséquente (équivalente à environ 25 tonnes par m²), typique pour des fondations industrielles lourdes. Cette valeur \(q\) est notre condition à la limite supérieure (\(z=0\)) pour le calcul de diffusion.
Une contrainte admissible classique pour un sable compact varie entre 200 et 400 kPa. Avec 250 kPa, nous sommes dans une plage réaliste et acceptable pour un sol de bonne qualité en surface. Si nous avions trouvé 2500 kPa, la semelle aurait été sous-dimensionnée (poinçonnement immédiat).
Assurez-vous d'utiliser les charges à l'ELS (État Limite de Service) pour les calculs de tassement et de déformation. Les charges à l'ELU (État Limite Ultime, pondérées par 1.35G + 1.5Q) ne servent qu'à vérifier la sécurité contre la rupture du sol (portance) et conduiraient à surestimer les tassements réels.
🎯 Objectif
La théorie de Boussinesq (1885) et son intégration par Steinbrenner pour les surfaces rectangulaires reposent sur des abaques ou des formules complexes. Pour utiliser ces outils, il est impératif de normaliser la géométrie du problème. L'objectif est de transformer nos dimensions physiques (\(B, L, z\)) en ratios adimensionnels (\(m, n\)) qui serviront de clés d'entrée universelles pour le calcul du coefficient d'influence.
📚 Référentiel
- Théorie de l'Élasticité Linéaire : Solution de Boussinesq pour un milieu semi-infini.
- Abaques de Steinbrenner : Coefficient d'influence sous le coin d'une charge rectangulaire.
C'est ici que réside la subtilité technique de l'exercice. La formule de Steinbrenner (et les abaques associés) donne la contrainte verticale sous le COIN d'un rectangle chargé. Or, nous cherchons la contrainte sous le CENTRE de la semelle (là où elle est maximale).
Nous ne pouvons pas utiliser directement \(B\) et \(L\). Nous devons utiliser le principe de superposition : nous découpons fictivement notre semelle réelle en 4 petits rectangles identiques de dimensions \(B/2\) et \(L/2\). Le centre de la grande semelle devient alors le "coin commun" de ces 4 petits rectangles. Nous calculerons l'influence d'un petit rectangle et multiplierons le résultat par 4.
Dans un milieu élastique linéaire, les contraintes sont additives. La contrainte induite en un point \(M\) par plusieurs charges est égale à la somme des contraintes induites par chaque charge prise individuellement.
C'est ce qui nous permet de "découper" des surfaces complexes en rectangles élémentaires.
📋 Données d'Entrée
| Paramètre | Valeur |
|---|---|
| Largeur Totale (\(B\)) | 2.40 m |
| Longueur Totale (\(L\)) | 3.60 m |
| Profondeur (\(z\)) | 5.00 m |
Les paramètres \(m\) et \(n\) sont interchangeables dans la formule de Steinbrenner (le résultat est symétrique). Peu importe lequel est \(L'/z\) ou \(B'/z\). L'important est d'utiliser les dimensions du rectangle élémentaire (quart de semelle), pas de la semelle entière !
Calculs Détaillés
1. Dimensions des Rectangles Fictifs (\(B', L'\)) :
On divise les dimensions réelles par 2 pour obtenir le quart de semelle.
2. Calcul des Facteurs Adimensionnels (\(m, n\)) :
On rapporte ces dimensions réduites à la profondeur \(z = 5.00\) m.
Nous obtenons deux valeurs sans unité, inférieures à 1, ce qui indique que la profondeur (\(z\)) est grande par rapport aux dimensions transversales de la fondation.
✅ Interprétation Globale
Nous avons correctement traduit la géométrie physique du problème en paramètres mathématiques normalisés (\(m=0.36, n=0.24\)). Ces valeurs sont prêtes à être injectées dans la fonction d'influence. Le fait qu'elles soient faibles suggère déjà intuitivement que l'influence sera réduite à cette profondeur.
Si vous trouvez \(m\) ou \(n > 10\), cela signifie que vous êtes très proche de la surface (ou que la semelle est immense). Si \(m\) ou \(n < 0.1\), vous êtes très profond et l'influence sera quasi nulle. Ici, avec 0.24 et 0.36, nous sommes dans la zone de "diffusion active" du bulbe de contraintes.
L'erreur la plus fréquente est d'oublier de diviser \(B\) et \(L\) par 2. Si vous utilisez \(B/z\) et \(L/z\) directement, vous calculez la contrainte sous le coin de la semelle entière, et non sous son centre ! Le résultat serait alors 4 fois trop faible.
🎯 Objectif
Nous devons maintenant déterminer le "facteur d'atténuation" géométrique. Ce coefficient, noté \(I_z\) (ou parfois \(I_\sigma\)), représente la proportion de la contrainte de surface (\(q\)) qui parvient à traverser le sol jusqu'à la profondeur \(z\). C'est une valeur purement géométrique, comprise entre 0 et 1 (ou 0 et 100%), indépendante de la charge appliquée.
📚 Référentiel
- Formule de Steinbrenner (1934) : Solution analytique pour la contrainte verticale sous le coin d'une aire rectangulaire uniformément chargée.
Imaginez que la charge se diffuse comme la lumière d'un projecteur. Plus on s'éloigne, plus le faisceau s'élargit et plus l'intensité lumineuse par m² diminue. Le coefficient \(I_z\) mesure cette "dilution". Pour nos valeurs \(m=0.36\) et \(n=0.24\), nous sommes assez profonds (\(z=5m\)) par rapport à la taille de la semelle (\(B=2.4m\)). On s'attend donc à ce que \(I_z\) soit faible, probablement inférieur à 20%.
Le sol diffuse les charges latéralement. Les lignes d'égale contrainte (isobares) forment ce qu'on appelle un "bulbe". Au droit du centre de la fondation, la contrainte verticale \(\sigma_z\) est maximale mais décroît rapidement avec la profondeur \(z\), suivant approximativement une loi en \(1/z^2\) à grande distance (similaire à une charge ponctuelle).
📋 Données d'Entrée
| Paramètre | Valeur |
|---|---|
| Ratio \(m\) | 0.36 |
| Ratio \(n\) | 0.24 |
Cette formule est horrible à calculer à la main ! Dans la pratique, ou lors d'un examen, on utilise souvent des abaques de Fadum ou des tables numériques. Ici, nous effectuons le calcul numérique exact pour la précision, mais une lecture d'abaque donnerait environ 0.035 - 0.040.
Calculs Détaillés
1. Calcul du Coefficient Unitaire (\(I_{\text{coin}}\)) :
En injectant \(m=0.36\) et \(n=0.24\) dans la formule massive ci-dessus (substitution numérique) :
Ceci est l'influence d'un seul quart de la semelle sur le point central.
2. Calcul du Coefficient Global (\(I_z\)) :
Puisque notre point de calcul est au centre de 4 rectangles identiques, nous sommons les 4 contributions (superposition) en multipliant par 4.
Le coefficient d'influence total est de 0.152 (ou 15.2%).
✅ Interprétation Globale
Le résultat \(I_z = 0.152\) signifie que seulement 15.2% de la pression appliquée en surface est ressentie à 5 mètres de profondeur. La grande majorité de la contrainte a été dissipée latéralement dans les couches de sable supérieures. C'est une bonne nouvelle pour la couche d'argile profonde.
Une règle empirique (méthode "2:1") suggère que la contrainte diminue avec le carré de la profondeur. À \(z=5m\), pour une semelle de largeur \(B \approx 2.5m\), on est à une profondeur \(2B\). On s'attend classiquement à un coefficient entre 0.10 et 0.20. Notre valeur de 0.152 est parfaitement cohérente.
Attention aux radians vs degrés dans la fonction `arctan` de vos calculatrices ! La formule de Boussinesq utilise généralement des radians. Une erreur ici fausserait complètement le résultat.
🎯 Objectif
C'est l'étape de synthèse finale. Nous allons combiner la pression de surface initiale (\(q\)), qui représente la "source", avec le facteur d'atténuation (\(I_z\)), qui représente la "propagation", pour obtenir la valeur réelle de la surcharge (\(\Delta \sigma_z\)) au niveau de l'argile. C'est cette valeur précise qui, ajoutée au poids des terres existantes (\(\sigma'_{\text{v0}}\)), provoquera le tassement de consolidation.
📚 Référentiel
- Calcul des Tassements : Application directe de la contrainte induite.
Nous avons quantifié la source (250 kPa) et le filtre (15.2%). Il ne reste plus qu'à appliquer le filtre à la source. Cette valeur \(\Delta \sigma_z\) est cruciale : si elle est trop élevée par rapport à la contrainte de préconsolidation de l'argile, les tassements seront importants (domaine vierge). Si elle est faible, nous resterons dans le domaine pseudo-élastique (tassements faibles).
En géotechnique, on distingue l'état initial (poids des terres seules) de l'état final (poids des terres + surcharge ouvrage). La valeur \(\Delta \sigma_z\) représente uniquement le "delta", c'est-à-dire le supplément de pression apporté par la nouvelle construction.
📋 Données d'Entrée
| Paramètre | Valeur |
|---|---|
| Contrainte de contact (\(q\)) | 250 kPa |
| Facteur d'influence (\(I_z\)) | 0.152 (sans unité) |
Pour vérifier mentalement : si \(I_z\) est d'environ 15%, le résultat doit être environ 15% de 250. 10% de 250 font 25, 5% font 12.5. La somme fait 37.5. On s'attend à un résultat proche de 38.
Calcul Final
1. Application Numérique :
Multiplication de la pression de surface par le coefficient d'atténuation.
L'argile subit une surcharge effective de 38.0 kPa.
✅ Conclusion de l'Étude
La couche d'argile située à 5 mètres de profondeur ne perçoit qu'une fraction de la charge de la presse. Au lieu des 250 kPa appliqués en surface, elle ne supporte que 38 kPa supplémentaires. C'est une réduction drastique (division par ~6.5) due à la bonne capacité de diffusion des sables sus-jacents. Cependant, 38 kPa reste une charge significative pour une argile molle, susceptible de générer plusieurs centimètres de tassement. La vérification oedométrique est donc indispensable.
Le résultat est physiquement logique : la contrainte diminue avec la profondeur. Si nous avions trouvé une valeur supérieure à 250 kPa, il y aurait eu une erreur grave (création d'énergie). Si nous avions trouvé 0.1 kPa, cela aurait signifié une diffusion improbable.
Ce calcul suppose le sol homogène pour la diffusion (module d'Young constant). Si la couche de sable superficielle est beaucoup plus rigide que l'argile (ce qui est le cas ici), la diffusion réelle sera en fait plus large (effet de "dalle"), et la contrainte sur l'argile sera encore plus faible que 38 kPa (théorie multicouche de Westergaard ou Burmister). Notre calcul Boussinesq est donc sécuritaire (il surestime légèrement la contrainte).
📄 Livrable Final (Note de Calculs EXE)
SOLUTIONS
| Ind. | Date | Objet de la modification | Rédacteur |
|---|---|---|---|
| A | 20/01/2024 | Création du document / Première diffusion | Ing. T. Martin |
- Eurocode 7 : Calcul géotechnique (NF EN 1997-1)
- NF P 94-261 : Fondations superficielles
| Charge ELS (\(Q\)) | 2160 kN |
| Dimensions (\(B \times L\)) | \(2.40 \times 3.60\) m |
| Profondeur cible (\(z\)) | 5.00 m |
Calcul selon la théorie de l'élasticité (Boussinesq/Steinbrenner) pour un milieu semi-infini, homogène et isotrope.
Ing. T. Martin
Dr. P. Dubois
Laisser un commentaire