Dans l’article précédent sur le principe des contraintes effectives de Terzaghi, nous avons établi la distinction fondamentale entre contraintes totales et contraintes effectives dans un massif de sol. Cette compréhension nous permet maintenant d’aborder une question cruciale pour tout dimensionnement géotechnique : comment se propagent et se distribuent les contraintes induites par les charges de structure dans la profondeur du sol ?
La distribution spatiale des contraintes conditionne directement la prédiction des tassements, le dimensionnement des fondations et l’évaluation des interactions entre ouvrages voisins. Une mauvaise estimation peut conduire à des sous-dimensionnements dangereux ou des sur-dimensionnements économiquement pénalisants.
Les solutions analytiques classiques développées par Boussinesq (1885) et Westergaard (1938) permettent de calculer cette distribution pour différentes configurations de sol et de chargement. Ces méthodes, bien que centenaires, restent aujourd’hui les outils de référence pour les calculs géotechniques courants, recommandées par l’Eurocode 7 (NF EN 1997-1).
Dans cet article, vous découvrirez :
Cette maîtrise technique est indispensable pour tout ingénieur géotechnicien impliqué dans le dimensionnement des fondations superficielles ou l’analyse de stabilité.
Comment se Propagent les Contraintes dans un Massif de Sol ?
Réponse directe : Les contraintes induites par une charge en surface se propagent dans le sol selon une distribution tridimensionnelle décroissante avec la profondeur, formant un « bulbe de contraintes ». L’intensité de la contrainte verticale diminue proportionnellement à l’inverse du carré de la profondeur, selon les solutions de la théorie de l’élasticité.
La propagation des contraintes dans un massif de sol peut être modélisée mathématiquement en considérant le sol comme un milieu élastique, homogène et isotrope. Cette approche, bien qu’idéalisée, fournit des résultats suffisamment précis pour la plupart des applications pratiques en géotechnique.
Hypothèses Fondamentales du Milieu Élastique
Les solutions classiques de Boussinesq et Westergaard reposent sur les hypothèses suivantes :
Ces hypothèses sont raisonnables pour des sols relativement homogènes sous des niveaux de contraintes modérés, typiquement inférieurs à 50% de la résistance ultime du sol. Pour les contraintes géostatiques naturelles, ces conditions sont généralement satisfaites.
Notion de Bulbe de Contraintes
Le concept de bulbe de contraintes (ou bulbe isobare) est fondamental en géotechnique. Il représente le lieu géométrique des points où la contrainte verticale induite σz prend une valeur constante, exprimée comme un pourcentage de la charge appliquée.
Caractéristiques du bulbe :
Pour une semelle carrée de 3m × 3m avec q = 200 kPa, la zone d’influence significative s’étend jusqu’à environ 5 mètres de profondeur, ce qui correspond à 1.67B.
Importance pour le Dimensionnement Géotechnique
La connaissance précise de la distribution des contraintes permet de :
- Prédire les tassements : En intégrant les déformations sur la profondeur influencée (voir article suivant sur les tassements et consolidation)
- Optimiser le dimensionnement : Éviter les sur-dimensionnements coûteux
- Évaluer les interactions : Analyser l’influence d’ouvrages voisins
- Vérifier la stabilité : S’assurer que les contraintes induites restent admissibles
L’Eurocode 7 (NF EN 1997-1) recommande explicitement l’utilisation de ces méthodes pour les calculs de tassement aux états limites de service (ELS), soulignant leur importance géotechnique dans la pratique professionnelle.
Quelle est la Solution de Boussinesq pour une Charge Ponctuelle ?
Réponse directe : La solution de Boussinesq exprime la contrainte verticale σz sous une charge ponctuelle Q comme : σz = (3Q/2πz²) × [1/(1+(r/z)²)^(5/2)], où z est la profondeur et r la distance radiale horizontale. Cette formule suppose un sol élastique, homogène et isotrope sur un demi-espace infini.
En 1885, le physicien français Joseph Boussinesq a développé une solution analytique fondamentale pour calculer la distribution des contraintes dans un massif élastique semi-infini soumis à une charge ponctuelle verticale en surface.
Démonstration et Équation Générale
La solution de Boussinesq découle de la résolution des équations d’équilibre de la mécanique des milieux continus en coordonnées cylindriques, avec les conditions aux limites appropriées.
Solution de Boussinesq – Charge Ponctuelle
Équation de Boussinesq (1885)
• σz = Contrainte verticale induite à la profondeur z (kPa)
• Q = Charge ponctuelle appliquée en surface (kN)
• z = Profondeur sous la surface (m)
• r = Distance radiale horizontale depuis l’axe de la charge (m)
• KB = Coefficient d’influence de Boussinesq
• π = 3.14159…
🔍 Points clés de la formule
→ La contrainte décroît avec le carré de la profondeur (1/z²)
→ Dépend fortement du rapport r/z (distance/profondeur)
→ Sous l’axe de charge (r = 0) : σz = 3Q/(2πz²)
→ Symétrie axiale : indépendant de l’angle θ
Coefficients d’Influence Kb
Pour faciliter les calculs pratiques, on introduit le coefficient d’influence de Boussinesq Kb :
σz = Kb × (Q/z²)
Avec :
Kb = (3/2π) × [1 / (1 + (r/z)²)^(5/2)]
Pour r = 0 (sous la charge) :
Le tableau complet des valeurs de Kb en fonction de r/z est présenté en section 6.
📊 Calculateur Boussinesq – Charge Ponctuelle
Calculez la distribution des contraintes dans le sol
🔢 Données du Problème
📐 Calculs Détaillés
À 3 mètres de profondeur et 2 mètres de décalage horizontal, la charge de 500 kN induit une contrainte verticale de 13.7 kPa. Cette valeur représente l’incrément de contrainte à ajouter aux contraintes géostatiques naturelles pour calculer les contraintes effectives totales.
📍 Représentation Graphique
Abaques d’Application Pratique
Des abaques ont été développés pour éviter les calculs répétitifs. Le plus utilisé est le diagramme de Newmark, qui permet de déterminer graphiquement l’influence de charges réparties quelconques par découpage en éléments et sommation.
Pour les applications informatiques modernes, les logiciels de calcul géotechnique (Plaxis, Geostudio) intègrent directement ces solutions analytiques avec des maillages d’éléments finis pour des géométries complexes.
Comment Calculer les Contraintes sous une Charge Répartie ?
Réponse directe : Pour une charge répartie, on applique le principe de superposition : le massif est découpé en éléments infinitésimaux, la solution de Boussinesq est appliquée à chaque élément, puis les contributions sont intégrées. Pour une semelle rectangulaire, des formules analytiques intégrées existent, et une méthode simplifiée « 2:1 » offre une approximation rapide avec diffusion à 2 verticale pour 1 horizontale.
La majorité des charges en géotechnique sont réparties (semelles, radiers, remblais) plutôt que ponctuelles. La solution de Boussinesq peut être adaptée à ces cas par intégration.
Principe de Superposition pour le Calcul des Contraintes
Le principe de superposition, valide pour les systèmes linéaires élastiques, stipule que :
« L’effet total de plusieurs charges est égal à la somme des effets de chaque charge prise individuellement. »
Pour une charge répartie uniforme q (kPa) sur une surface A :
- Découper la surface en éléments infinitésimaux dA
- Chaque élément porte une charge dQ = q × dA
- Appliquer Boussinesq à chaque dQ
- Intégrer sur toute la surface chargée
Cette approche permet de traiter n’importe quelle géométrie de fondation, une capacité essentielle pour le dimensionnement des fondations superficielles.
📐 Contraintes sous Semelle Filante (Solution 2D)
Calcul de la distribution des contraintes pour fondations continues
📖 Définition : Semelle Filante
Une semelle filante (ou semelle continue) est une fondation dont la longueur L est très supérieure à sa largeur B (typiquement L > 5B). Cette configuration permet une analyse bidimensionnelle (2D) en déformation plane.
Exemples d’application : Murs de soutènement, fondations de bâtiments en bande, voies ferrées, pipelines enterrés.
📐 Formule Générale (Solution Intégrée)
Où :
• α et β = angles géométriques fonction de la position (x, z)
• α = angle sous-tendu par la semelle depuis le point considéré
• q = charge uniforme appliquée (kPa)
• Cette formulation suppose une déformation plane (état 2D)
🎯 Cas Simplifié : Sous le Centre de la Semelle (x = 0)
📊 Valeurs Caractéristiques
🔢 Paramètres d’Entrée
📊 Résultats du Calcul
À 2 m de profondeur (z = B), la contrainte représente environ 55% de la charge appliquée, conforme aux valeurs théoriques.
📉 Graphique de Décroissance : σz en fonction de z/B
✏️ Exercice d’Application : Testez Vos Connaissances
📝 Énoncé du Problème
Une semelle filante de largeur B = 3 m supporte une charge uniforme q = 150 kPa.
Question : Calculez la contrainte verticale σz sous le centre de la semelle à une profondeur z = 4.5 m (z = 1.5B).
σz = q × [1 – z/√(z² + (B/2)²)]
📐 Contraintes sous Semelle Rectangulaire
Distribution 3D des contraintes – Coefficient d’influence I
📖 Semelle Rectangulaire : Cas Général 3D
Pour une semelle rectangulaire de dimensions B × L (B = largeur, L = longueur) soumise à une charge uniforme q, la solution de Boussinesq est intégrée sur toute la surface de la fondation.
Contrairement à la semelle filante (2D), c’est un problème tridimensionnel (3D) nécessitant l’utilisation de coefficients d’influence I tabulés.
📐 Formule Générale – Coin de la Semelle
Où :
• σz = Contrainte verticale à la profondeur z (kPa)
• q = Charge uniforme appliquée (kPa)
• I = Coefficient d’influence fonction de m et n
• m = L/B (rapport longueur/largeur)
• n = z/B (rapport profondeur/largeur)
🎯 Cas du Centre de la Semelle
Pour calculer la contrainte sous le centre de la semelle, on utilise le principe de symétrie : le centre est le coin commun de 4 rectangles identiques B/2 × L/2.
📊 Table de Coefficients I (Das 2016)
Extrait pour m = 1.5 (semelle 3m × 2m)
| n = z/B | I (coin) | σz/q (%) |
|---|---|---|
| 0.0 | 0.2500 | 100% (surface) |
| 0.5 | 0.2285 | 91% |
| 1.0 | 0.1774 | 71% |
| 1.5 | 0.1349 | 54% |
| 2.0 | 0.1047 | 42% |
| 2.5 | 0.0843 | 34% |
📍 Vue 3D – Semelle Rectangulaire B × L
📊 EXEMPLE NUMÉRIQUE : Semelle 3m × 2m, q = 200 kPa
Données du problème :
- Dimensions : B = 2 m (largeur), L = 3 m (longueur)
- Charge uniforme : q = 200 kPa
- Point de calcul : Centre de la semelle, z = 2 m de profondeur
💡 Ce coefficient est obtenu par interpolation dans les tables de Boussinesq intégrées
✅ Résultat Final
Au centre de la semelle, à 2 mètres de profondeur (z = B), la contrainte induite représente 54% de la charge appliquée (107.9/200 = 0.54).
Cette valeur est cohérente avec les observations pratiques et permet de calculer le tassement par la méthode œdométrique.
⚠️ Vérification avec méthode simplifiée 2:1
L’écart de 44% montre que la méthode 2:1 sous-estime significativement les contraintes pour z = B, mais reste acceptable pour des estimations préliminaires à z ≥ 2B.
🔢 Paramètres d’Entrée
📊 Résultats du Calcul
Au centre, à z = B, la contrainte représente 54% de la charge appliquée.
✏️ Exercice d’Application : Testez Vos Connaissances
📝 Énoncé du Problème
Une semelle rectangulaire de dimensions B = 2.5 m × L = 4 m supporte une charge uniforme q = 180 kPa.
Question : Calculez la contrainte verticale σz sous le centre de la semelle à une profondeur z = 2.5 m (z = B).
Utilisez I ≈ 0.145 pour m = 1.6 et n’ = 2.0
σz(centre) = 4 × q × I
⚡ Méthode Simplifiée 2:1 (Trapézoïdale)
Diffusion des contraintes – Approche conservative rapide
📖 Principe de la Méthode 2:1
La méthode 2:1 (ou méthode trapézoïdale) est une approche simplifiée et conservative pour estimer rapidement la distribution des contraintes sous une fondation.
Principe : La charge se diffuse avec une pente de 2 vertical : 1 horizontal, soit un angle de 26.6° avec la verticale.
📐 Formule Générale
Où :
• σz = Contrainte verticale à la profondeur z (kPa)
• q = Charge uniforme appliquée (kPa)
• B = Largeur de la semelle (m)
• L = Longueur de la semelle (m)
• z = Profondeur sous la base de la fondation (m)
📍 Principe de Diffusion 2:1
✅ Avantages de la Méthode 2:1
- Calcul mental rapide : Formule simple, pas besoin de tables
- Conservative : Sous-estime généralement les contraintes (sécurité)
- Acceptable pour z ≥ 2B : Erreur réduite en profondeur
- Pédagogique : Visualisation intuitive de la diffusion
❌ Inconvénients et Limites
- Imprécise pour z < 2B : Erreur jusqu’à 50% près de la fondation
- Ne tient pas compte de la position : Même résultat centre/coin
- Non conforme EC7 : Pas pour calculs de justification finale
- Sous-estime : Peut conduire à des sur-dimensionnements
⚠️ Recommandations Eurocode 7
L’Eurocode 7 recommande d’utiliser cette méthode uniquement pour des estimations préliminaires, pas pour les calculs de justification finale. Pour les projets définitifs, utiliser Boussinesq ou des méthodes numériques.
🔢 Paramètres d’Entrée
📊 Résultats du Calcul
📊 Comparaison avec Boussinesq
À cette profondeur (z/B = 2.0), la méthode 2:1 est relativement précise.
La charge s’est diffusée sur une surface 7 fois plus grande, réduisant la contrainte à 14.3% de q.
⚖️ Comparaison Boussinesq vs Méthode 2:1
| Critère | Boussinesq (Exacte) | Méthode 2:1 (Simplifiée) |
|---|---|---|
| Formule | σz = 4 × q × I(m, 2n) Nécessite tables |
σz = q × BL/[(B+z)(L+z)] Calcul direct |
| Complexité | 🔴 Complexe | 🟢 Simple |
| Précision z < B | 🟢 Excellente (±5%) | 🔴 Mauvaise (-50%) |
| Précision z = B | 🟢 Très bonne (±10%) | 🟡 Moyenne (-40%) |
| Précision z ≥ 2B | 🟢 Bonne (±15%) | 🟢 Acceptable (-20%) |
| Usage EC7 | 🟢 Recommandé pour calculs définitifs | 🟡 Estimations préliminaires uniquement |
| Temps de calcul | ⏱️ 2-5 min (avec tables) | ⚡ 30 secondes |
📊 Exemple Comparatif : Semelle 3m × 2m, q = 200 kPa, z = 2m
Méthode Boussinesq
σz = 4 × 200 × 0.1349
σz = 107.9 kPa
54% de q
Méthode 2:1
σz = 200 × 6 / (4 × 5)
σz = 60 kPa
30% de q
La méthode 2:1 sous-estime significativement pour z = B. Acceptable uniquement pour z ≥ 2B ou estimations préliminaires.
✏️ Exercice d’Application : Testez Vos Connaissances
📝 Énoncé du Problème
Une semelle rectangulaire de dimensions B = 3 m × L = 5 m supporte une charge uniforme q = 150 kPa.
Question : En utilisant la méthode 2:1, calculez la contrainte verticale σz à une profondeur z = 6 m (z = 2B).
σz = q × (B × L) / [(B + z)(L + z)]
Quand Utiliser la Méthode de Westergaard pour le Calcul des Contraintes ?
Réponse directe : La solution de Westergaard doit être privilégiée pour les sols stratifiés avec couches horizontales de rigidités très différentes, typiquement argiles varvées ou sols avec intercalations raides. Elle suppose une rigidité infinie dans les plans horizontaux, réduisant la diffusion latérale et augmentant les contraintes verticales de 20 à 40% par rapport à Boussinesq selon Das (2016).
En 1938, le danois Harald Malcolm Westergaard a développé une solution alternative pour les sols présentant une anisotropie marquée, fréquente dans les dépôts sédimentaires stratifiés.
Distribution des Contraintes dans les Sols Stratifiés et Anisotropes
Contexte géologique : Certains sols présentent une structure en couches horizontales avec des propriétés mécaniques très différentes :
Dans ces configurations, la rigidité horizontale (dans les couches) est nettement supérieure à la rigidité verticale (perpendiculaire aux couches), créant un comportement anisotrope transverse.
Cette situation est fréquente dans : – Les argiles varvées du Canada et de Scandinavie – Les bassins sédimentaires français (Bassin Parisien) – Les dépôts alluviaux stratifiés
Formule de Westergaard pour Sols Anisotropes
Formule de Westergaard pour charge ponctuelle Q :
Ou avec coefficient d’influence Kw :
Avec :
Sous l’axe de la charge (r = 0) :
À comparer avec Kb = 0.4775 pour Boussinesq, soit une réduction de 33% sous la charge.
👉 Westergaard prédit environ 33 % de contrainte verticale en moins sous l’axe de la charge.
C’est exactement l’interprétation classique :
Comparaison Boussinesq vs Westergaard : Différences de Distribution des Contraintes
| Critère | Boussinesq | Westergaard |
|---|---|---|
| Hypothèse sol | Isotrope (E identique toutes directions) | Anisotrope transverse (rigidité horizontale >> verticale) |
| Type de sol | Sables, argiles homogènes | Argiles varvées, sols stratifiés |
| Diffusion latérale | Large (bulbe élargi) | Réduite (bulbe resserré) |
| σz sous charge (r=0) | Plus élevée (Kb=0.478) | Plus faible (Kw=0.318) |
| σz en profondeur | Décroissance modérée | Décroissance plus rapide |
| Usage EC7 | Standard, recommandé | Cas particuliers justifiés |
Retour d’expérience et littérature technique
De nombreuses études expérimentales et observations de terrain rapportées dans la littérature géotechnique montrent que, pour des sols fortement stratifiés ou anisotropes (argiles varvées, sols laminés), la solution de Westergaard fournit souvent une estimation plus réaliste de la diffusion des contraintes que la solution de Boussinesq, cette dernière ayant tendance à surestimer les contraintes verticales sous l’axe de la charge.
📊 EXEMPLE NUMÉRIQUE 3 : Sol Argileux Stratifié
📋 Données du problème :
- Site : Bassin Parisien, argiles vertes de Romainville (stratifiées)
- Charge ponctuelle : Q = 1000 kN (poteau)
- Point de calcul : z = 4 m, r = 1.5 m
- Reconnaissance : Essais pressiométriques montrant anisotropie Eh/Ev ≈ 3
⚖️ Comparaison Boussinesq vs Westergaard :
Méthode 1 Boussinesq (isotrope)
Kb = (3/2π) × [1/(1 + 0.375²)^2.5]
Kb = 0.4775 × [1/1.384]
Kb = 0.345
σz(Boussinesq) = 0.345 × (1000/16)
σz(Boussinesq) = 21.6 kPa
Méthode 2 Westergaard (anisotrope)
Kw = 1 / (π[1 + 2×0.375²]^1.5)
Kw = 1 / (π[1.281]^1.5)
Kw = 1 / (π × 1.449)
Kw = 0.220
σz(Westergaard) = 0.220 × (1000/16)
σz(Westergaard) = 13.7 kPa
📐 Écart calculé
💡 Interprétation :
Boussinesq surestime la contrainte de 58% pour ce sol stratifié. Cette différence majeure impacte directement le calcul des tassements et peut conduire à un sur-dimensionnement coûteux si Boussinesq est utilisé sans correction pour l’anisotropie du sol.
Recommandation projet : Pour ce site, Westergaard doit être utilisé, conformément à la norme NF P 94-261 qui prescrit son usage pour les sols avec anisotropie structurale marquée. La classification des sols issue des sondages doit justifier ce choix.
📊 TABLEAU 1 : Comparatif Boussinesq vs Westergaard
| Aspect | Boussinesq (1885) | Westergaard (1938) |
|---|---|---|
| Sol type | Homogène, isotrope | Stratifié, anisotrope |
| Diffusion | Large (bulbe élargi) | Étroite (bulbe resserré) |
| Coefficient sous charge | Kb = 0.478 | Kw = 0.318 (-33%) |
| Contrainte à z=B | ≈ 55% de q | ≈ 40% de q (-27%) |
| Précision sols sableux | Excellente (±10%) | Médiocre (sous-estime) |
| Précision argiles varvées | Médiocre (surestime 40%) | Excellente (±15%) |
| Profondeur influence | 1.5 à 2B | 1.2 à 1.5B (réduite) |
| Complexité calcul | Modérée | Similaire |
| Usage EC7 | Standard | Justifié par reconnaissance |
| Logiciels | Tous | Plaxis, Geostudio |
Quelles sont les Limites et Précautions d’Utilisation ?
Réponse directe : Les principales limites sont l’hypothèse d’élasticité linéaire (invalide près de la rupture), l’ignorance de la rigidité de la fondation (qui modifie la distribution), et l’hypothèse d’homogénéité (rare en pratique). L’Eurocode 7 recommande de les utiliser pour σz < 0.3σ’rupture et de considérer la stratigraphie réelle par calculs multicouches ou éléments finis pour projets critiques.
Bien que les solutions de Boussinesq et Westergaard soient largement utilisées, leurs hypothèses simplificatrices imposent des limites d’application qu’il est essentiel de connaître.
Hypothèse d’Élasticité Linéaire
Limite fondamentale : La loi de comportement élastique linéaire (σ = E × ε) n’est valable que pour de faibles niveaux de contraintes.
Critères de validité (CFMS 2024) :
Conséquence pratique : Près des charges concentrées ou sous des fondations fortement chargées, le sol développe des zones plastiques où Boussinesq sous-estime les contraintes et surestime la diffusion latérale.
Correction possible : Pour les calculs avancés, utiliser des modèles élasto-plastiques (Mohr-Coulomb, critère de rupture) dans des logiciels éléments finis.
Influence de la Rigidité de la Fondation
Les solutions analytiques supposent une charge appliquée directement à la surface du sol. En réalité, la rigidité de la fondation (semelle béton, radier) modifie significativement la distribution :
Fondation parfaitement flexible :
Fondation parfaitement rigide :
Cas réel (fondation semi-rigide) : La réalité se situe entre ces deux extrêmes. Le rapport de rigidité fondation/sol (facteur K) gouverne le comportement :
K = (Ef × If) / (Es × Bs³)
Où :
- Ef = module d’Young fondation (béton : 30 GPa)
- If = moment d’inertie fondation
- Es = module d’Young sol (sable : 20-50 MPa, argile : 5-20 MPa)
- Bs = dimension caractéristique
Correction (méthode de Steinbrenner) : Pour tenir compte de la rigidité, des facteurs correctifs sont appliqués aux solutions élastiques, tabulés dans la norme NF P 94-261.
Cas des Sols Multicouches
L’hypothèse d’homogénéité est rarement vérifiée en pratique. La stratigraphie réelle présente généralement :
Problématique : Boussinesq/Westergaard ne modélisent pas la réfraction des contraintes aux interfaces entre couches de rigidités différentes, analogue à la réfraction de la lumière en optique.
Solutions pratiques :
- Méthode stratifiée multicouche :
- Découper le massif en n couches homogènes
- Appliquer Boussinesq dans chaque couche avec module équivalent
- Itérer pour continuité des contraintes aux interfaces
- Programmable sous Excel (CFMS 2024)
- Modélisation éléments finis :
- Logiciels Plaxis, Geostudio, Cesar-LCPC
- Stratigraphie réelle issue des sondages
- Modèles de comportement avancés
- Obligatoire pour ouvrages critiques (EC7, catégorie géotechnique 3)
- Méthode de Steinbrenner :
- Facteurs correctifs pour bicouche (couche compressible sur substratum rigide)
- Tables dans Das (2016) et NF P 94-261
Recommandations Normatives (Eurocode 7)
L’Eurocode 7 (NF EN 1997-1, article 6.6) encadre l’usage de ces méthodes :
Pour calculs de tassement (ELS) : ✓ Boussinesq acceptable si sol relativement homogène sur profondeur d’influence ✓ Westergaard si argiles varvées ou anisotropie justifiée ✓ Méthodes numériques recommandées si hétérogénéité marquée
Pour vérification capacité portante (ELU) : ✗ Boussinesq/Westergaard NON ADAPTÉS (domaine plastique) ✓ Utiliser formules de capacité portante (Terzaghi, Meyerhof, EC7)
Catégories géotechniques :
- GC1 (projets simples) : Boussinesq acceptable avec méthode 2:1
- GC2 (standard) : Boussinesq/Westergaard avec reconnaissance adaptée
- GC3 (complexe) : Modélisation numérique obligatoire
Mission géotechnique (NF P 94-500) : Une mission G2 AVP minimum est requise pour définir les paramètres E, ν nécessaires aux calculs de distribution des contraintes.
📊 Tableau des Coefficients d’Influence Kb (Boussinesq)
Pour faciliter les calculs manuels ou sous Excel, voici le tableau complet des coefficients Kb en fonction du rapport r/z :
| r/z | Kb (Boussinesq) | σz/Q (pour z=1m) | r/z | Kb (Boussinesq) | σz/Q (pour z=1m) |
|---|---|---|---|---|---|
| 0.0 | 0.4775 | 0.478 | 1.2 | 0.1386 | 0.139 |
| 0.1 | 0.4657 | 0.466 | 1.4 | 0.1054 | 0.105 |
| 0.2 | 0.4329 | 0.433 | 1.6 | 0.0815 | 0.082 |
| 0.3 | 0.3849 | 0.385 | 1.8 | 0.0643 | 0.064 |
| 0.4 | 0.3295 | 0.330 | 2.0 | 0.0517 | 0.052 |
| 0.5 | 0.2733 | 0.273 | 2.5 | 0.0317 | 0.032 |
| 0.6 | 0.2214 | 0.221 | 3.0 | 0.0213 | 0.021 |
| 0.7 | 0.1762 | 0.176 | 3.5 | 0.0151 | 0.015 |
| 0.8 | 0.1386 | 0.139 | 4.0 | 0.0111 | 0.011 |
| 0.9 | 0.1083 | 0.108 | 4.5 | 0.0084 | 0.008 |
| 1.0 | 0.0844 | 0.084 | 5.0 | 0.0066 | 0.007 |
Mode d’emploi :
- Calculer r/z pour votre point d’intérêt
- Lire Kb dans le tableau (interpolation linéaire si nécessaire)
- Calculer σz = Kb × (Q/z²)
Exemple d’utilisation :
- Charge Q = 800 kN
- Point : z = 2.5 m, r = 2 m
- r/z = 2/2.5 = 0.8
- Kb = 0.1386 (lecture directe)
- σz = 0.1386 × (800/6.25) = 17.7 kPa
Observations :
- À r/z = 0 (sous charge) : Kb maximum = 0.4775
- À r/z = 2.0 : Kb = 0.0517 (10.8% du max) → limite zone d’influence
- À r/z = 3.0 : Kb = 0.0213 (4.5% du max) → influence négligeable
Ce tableau est conforme aux valeurs publiées par Bowles (1996) et Das (2016), références internationales en géotechnique.
❓ FAQ – 5 Questions Essentielles sur la Distribution des Contraintes
Cliquez sur chaque question pour afficher la réponse détaillée
Conclusion et Synthèse
Les méthodes de Boussinesq (1885) et Westergaard (1938) constituent des outils analytiques fondamentaux pour calculer la distribution des contraintes dans les massifs de sol soumis à des charges de structure.
À retenir (4 points clés) :
✓ Boussinesq : Solution standard pour sols homogènes isotropes, formule σz = Kb(Q/z²), applicable à la majorité des projets courants
✓ Westergaard : Alternative pour sols stratifiés anisotropes (argiles varvées), réduit les contraintes verticales de 20-40% par rapport à Boussinesq
✓ Superposition : Le principe de superposition permet de traiter charges réparties et multiples par intégration ou sommation, avec méthode 2:1 pour estimations rapides
✓ Limites : Hypothèses d’élasticité linéaire, homogénéité et semi-espace infini imposent des précautions, modélisation numérique nécessaire pour géométries complexes
Continuité de la formation :
Maintenant que vous maîtrisez le calcul des contraintes induites dans le sol, la question suivante se pose naturellement : comment ces contraintes se traduisent-elles en déformations et tassements ?
Dans le prochain article « Compressibilité et Tassements des Sols » (Module 2.2.1), vous découvrirez :
- Le mécanisme de compression et consolidation des sols
- Le calcul des tassements par la méthode œdométrique
- L’intégration des contraintes de Boussinesq pour prédire les déformations
- Les tassements différentiels et leurs conséquences structurelles
Les contraintes que vous savez maintenant calculer deviennent l’entrée du calcul de tassement, bouclant ainsi le cycle : charge → contrainte → déformation → tassement.
👉 Découvrez comment calculer les tassements à partir de la distribution des contraintes →
📚 Pour Aller Plus Loin
Approfondissez vos connaissances avec ces articles connexes :
- Principe des Contraintes Effectives de Terzaghi – Fondement théorique de la mécanique des sols
- Contraintes Géostatiques : Calcul et Distribution – État de contrainte naturel du sol
- Capacité Portante des Fondations Superficielles – Application au dimensionnement
- Essai Œdométrique et Consolidation – Détermination des paramètres de compressibilité
- Missions Géotechniques G1 à G5 – Cadre normatif de la reconnaissance
Cet article vous a été utile ?
Partagez-le avec vos collègues ingénieurs et n’hésitez pas à nous faire part de vos retours d’expérience sur l’application de ces méthodes à vos projets !
