Calcul du gradient hydraulique critique
📝 Situation du Projet
Dans le cadre de l'aménagement urbain du quartier "Rives de Seine", le projet prévoit la construction d'un parking souterrain de deux niveaux (R-2). Ce site est particulièrement complexe d'un point de vue géotechnique : il est situé en bordure immédiate du fleuve, dans une zone constituée d'alluvions sableuses récentes. Le niveau de la nappe phréatique est quasi-affleurant, fluctuant directement avec les marées et les crues du fleuve voisin.
Pour permettre le terrassement à sec de la fouille, l'entreprise générale a mis en œuvre une enceinte étanche constituée d'un rideau de palplanches métalliques de type AZ-18. Cependant, lors des phases d'excavation, une différence de charge hydraulique importante se crée entre l'extérieur (nappe haute) et l'intérieur de la fouille (fond de fouille asséché). Cette différence de potentiel génère un écoulement vertical ascendant sous la fiche du rideau.
Ce flux d'eau ascendant exerce une force de traînée sur les grains de sable du fond de fouille. Si cette force dépasse le poids effectif des terres, nous assistons à un phénomène de rupture catastrophique appelé "Renard Hydraulique" ou "Boulance". Le sol perd toute sa résistance au cisaillement et se comporte comme un fluide bouillant, entraînant l'inondation brutale du chantier et la ruine potentielle des soutènements. Votre rôle est de quantifier ce risque avec précision.
En tant qu'Ingénieur Géotechnicien au sein du bureau d'études, vous êtes mandaté pour valider la stabilité hydraulique du fond de fouille à la cote finale d'excavation. Vous devez mener une analyse déterministe en calculant d'abord le gradient critique théorique du matériau en place, puis en estimant le gradient réel induit par la géométrie de l'écran. La comparaison de ces deux valeurs vous permettra de définir le coefficient de sécurité global et de statuer sur la nécessité de mesures compensatoires (approfondissement de fiche, bouchon injecté, ou pompage).
"Attention, le phénomène de renard est insidieux. Si le gradient hydraulique effectif dépasse le gradient critique, la contrainte effective s'annule : le sable se comporte comme un liquide lourd. La marge de sécurité doit être stricte (FS > 1.5 minimum selon l'Eurocode)."
Les paramètres ci-dessous proviennent de la synthèse des reconnaissances géotechniques (Sondages CPT, essais en laboratoire triaxiaux) et des plans d'exécution des parois moulées. Ces données constituent le socle technique intangible pour votre note de calcul.
📚 Référentiel Normatif & Physique
L'étude s'inscrit strictement dans le cadre réglementaire européen et français pour les ouvrages géotechniques.
Eurocode 7 (EN 1997-1) : Calcul Géotechnique Loi de Darcy (Hydraulique souterraine)Le sol en place est un sable alluvial propre, homogène et non plastique. Les essais d'identification en laboratoire ont permis de déterminer les propriétés intrinsèques des grains et de l'arrangement granulaire :
- Densité des grains (\(G_s\)) : Mesurée au pycnomètre à hélium, elle représente le rapport entre la masse volumique du squelette solide et celle de l'eau. Une valeur de 2.65 est standard pour des sables siliceux (quartz).
- Indice des vides (\(e\)) : Ce paramètre sans dimension traduit la compacité du sol. Ici, \(e=0.65\) indique un sable moyennement dense, laissant un volume de vide significatif pour la circulation de l'eau.
| IDENTIFICATION DU SABLE | |
| Densité des grains solides (Specific Gravity) | \(G_s = 2.65\) |
| Indice des vides (Void Ratio) | \(e = 0.65\) |
| FLUIDE INTERSTITIEL | |
| Poids volumique de l'eau | \(\gamma_w = 9.81\) kN/m³ |
📐 Géométrie Hydraulique
Les dimensions critiques pour l'analyse de stabilité proviennent des plans de phasage. Elles définissent le chemin que l'eau doit parcourir et la "force motrice" (différence de hauteur d'eau) qui la pousse.
- Charge Hydraulique (\(H\)) : C'est le moteur de l'instabilité. Elle correspond à la différence d'altitude entre le niveau piézométrique extérieur et le fond de la fouille. Pour le cas le plus défavorable (crue), elle est fixée à 4.50 m.
- Fiche Mécanique (\(D\)) : C'est la longueur de palplanche enfoncée sous le fond de fouille. Elle joue le rôle de "frein" hydraulique en allongeant le parcours de l'eau. Elle est actuellement de 3.50 m.
- Longueur de Percolation (\(L_{\text{flow}}\)) : C'est la distance effective parcourue par une particule d'eau le long de la ligne de courant la plus courte (longeant le rideau). Pour cet exercice, nous utiliserons l'approximation conventionnelle \(L_{\text{flow}} = D + \frac{D}{2}\).
Cette vue microscopique illustre le combat qui se joue à l'échelle du grain de sable. Deux forces s'opposent : la gravité (qui maintient le grain en place) et la force de courant (qui tente de le soulever). L'état critique est atteint lorsque ces deux forces s'équilibrent parfaitement.
| Donnée | Symbole | Valeur | Unité |
|---|---|---|---|
| Densité relative | \(G_s\) | 2.65 | [-] |
| Indice des vides | \(e\) | 0.65 | [-] |
| Charge Hydraulique Différentielle | \(H\) | 4.50 | [m] |
| Profondeur Fiche | \(D\) | 3.50 | [m] |
E. Protocole de Résolution
Pour garantir la stabilité de l'ouvrage, nous allons procéder par étapes successives, en partant des propriétés intrinsèques du sol pour arriver au coefficient de sécurité global.
Calcul des Poids Volumiques
Déterminer le poids volumique saturé du sol (\(\gamma_{\text{sat}}\)) à partir des paramètres \(G_s\) et \(e\).
Gradient Critique Théorique (\(i_c\))
Calculer la valeur limite du gradient hydraulique pour laquelle la contrainte effective s'annule.
Gradient Réel de l'Ouvrage (\(i\))
Estimer le gradient hydraulique effectif généré par la géométrie actuelle de la fouille et du rideau.
Coefficient de Sécurité (\(F_s\))
Comparer le gradient critique au gradient réel pour statuer sur la stabilité (Norme > 1.5).
Calcul du gradient hydraulique critique
🎯 Objectif Scientifique
L'objectif de cette première étape est de caractériser l'état physique du sol dans son environnement naturel inondé. Avant de parler de forces ou de gradients, nous devons déterminer la "densité" réelle du matériau lorsqu'il est saturé d'eau. La grandeur recherchée est le Poids Volumique Saturé (\(\gamma_{\text{sat}}\)). Cette valeur est fondamentale car elle conditionne directement le poids du squelette solide qui s'opposera à la pression de l'eau. Sans cette valeur précise, aucun calcul de stabilité ultérieur n'est possible.
📚 Référentiel Technique
Mécanique des Milieux Poreux Relations Volumétriques (Diagramme des phases)Nous sommes face à un problème d'hydraulique souterraine : le sol est considéré comme biphasique (grains solides + eau), sans air. Les données d'entrée fournies par le laboratoire sont des paramètres intrinsèques "séparés" : la densité relative des grains (\(G_s\)) et l'indice des vides (\(e\)). Notre tâche intellectuelle consiste à "reconstruire" mentalement un volume unitaire de sol (1 m³) pour en déduire son poids total. Nous savons que les vides sont remplis d'eau. Le poids total sera donc la somme du poids des grains et du poids de l'eau contenue dans les pores. C'est une application directe du changement de variables d'état.
En mécanique des sols, un volume de sol \(V\) est composé d'un volume de solides \(V_s\) et d'un volume de vides \(V_v\). L'indice des vides est défini par \(e = V_v / V_s\). Si l'on considère un volume de solides unitaire \(V_s = 1\), alors le volume des vides est \(e\) et le volume total est \(1+e\).
Dans un état saturé, tous les vides sont remplis d'eau. Le poids de l'eau est donc \(e \cdot \gamma_w\) et le poids des solides est \(G_s \cdot \gamma_w\). Le poids volumique total est le rapport Poids Total / Volume Total.
Étape 1 : Données d'Entrée
| Paramètre | Symbole | Valeur |
|---|---|---|
| Densité des grains | \(G_s\) | 2.65 |
| Indice des vides | \(e\) | 0.65 |
| Poids vol. eau | \(\gamma_w\) | 9.81 kN/m³ |
Avant même de calculer, ayez en tête les ordres de grandeur. Pour des sables et graviers, \(\gamma_{\text{sat}}\) se situe presque toujours entre 19 et 21 kN/m³. Une valeur de 15 kN/m³ serait celle d'une vase très lâche, et 24 kN/m³ celle d'un béton. Si vous sortez de cette plage, recalculez !
Étape 2 : Application Numérique Détaillée
1. Substitution des valeurs :
Nous remplaçons les symboles littéraux par les valeurs numériques fournies dans l'énoncé.
Le terme (3.30 / 1.65) représente la densité relative globale du mélange sol+eau.
2. Calcul final :
On effectue la division (qui donne exactement 2) puis la multiplication finale.
Interprétation : Chaque mètre cube de ce sol, une fois gorgé d'eau, pèse 19.62 kiloNewtons (soit environ 2 tonnes-force). C'est ce poids qui va devoir résister à la pression ascendante.
✅ Conclusion de l'étape
Nous avons déterminé que le sol saturé pèse 19.62 kN/m³. Cette valeur est la base de tout l'équilibre vertical. C'est la "masse" (en termes de force) que l'eau devra soulever pour provoquer la rupture.
Le résultat obtenu (19.62 kN/m³) est parfaitement cohérent avec un sable alluvial compact (\(G_s=2.65\)). Il est supérieur au poids de l'eau (9.81), ce qui est logique puisque le sol coule dans l'eau. Il est inférieur à la densité de la roche (26 kN/m³) à cause des vides.
Ne confondez pas \(\gamma_{\text{sat}}\) (poids volumique saturé) avec \(\gamma_d\) (poids volumique sec) ou \(\gamma_{\text{hum}}\) (poids volumique humide naturel). Dans un calcul de boulance sous la nappe, c'est impérativement le poids saturé qu'il faut utiliser initialement.
🎯 Objectif Scientifique
Cette étape vise à déterminer la limite de résistance absolue du sol. Nous cherchons à calculer le Gradient Critique (\(i_c\)). C'est une valeur théorique, propre au matériau, qui correspond à la pente hydraulique exacte pour laquelle les forces de frottement de l'eau compensent intégralement le poids apparent des grains. À cet instant précis, la contrainte effective s'annule (\(\sigma' = 0\)). Le sol est en état d'apesanteur apparente : c'est le seuil mathématique de la liquéfaction statique (boulance).
📚 Référentiel
Principe de Terzaghi (Contrainte Effective) Équilibre limite statiqueIl est crucial de comprendre que le gradient critique ne dépend pas de la hauteur d'eau dans le fleuve, ni de la profondeur de la fouille. C'est une propriété intrinsèque du sol, au même titre que sa densité. Un sol "lourd" (forte densité de grains, faible porosité) aura un gradient critique plus élevé et sera donc plus difficile à soulever qu'un sol "léger". Notre calcul va consister à comparer le poids du sol sous l'eau (poids déjaugé) au poids de l'eau elle-même.
Un grain de sol immergé est soumis à son poids (vers le bas) et à la poussée d'Archimède (vers le haut). La résultante est le "Poids Déjaugé" (\(\gamma'\)). L'eau en mouvement exerce une troisième force : la force d'écoulement \(j = i \cdot \gamma_w\). L'équilibre critique est atteint quand la force d'écoulement (vers le haut) égale le poids déjaugé (vers le bas).
1. Démonstration du Poids Déjaugé
Le principe d'Archimède s'applique au volume total :
2. Démonstration du Gradient Critique
On écrit l'équilibre des forces verticales (Poids = Force d'écoulement). À la limite critique, la contrainte effective est nulle :
Étape 1 : Données nécessaires
| Paramètre | Valeur |
|---|---|
| Poids vol. saturé \(\gamma_{\text{sat}}\) (Calculé en Q1) | 19.62 kN/m³ |
| Poids vol. eau \(\gamma_w\) | 9.81 kN/m³ |
Pour la grande majorité des sols minéraux courants, le gradient critique est très proche de 1.0 (souvent entre 0.9 et 1.1). Si vous trouvez 0.5 ou 2.0, vérifiez vos calculs, c'est très improbable pour un sol naturel standard.
Étape 2 : Calculs Détaillés
1. Calcul du Poids Volumique Déjaugé (\(\gamma'\)) :
Nous soustrayons la poussée d'Archimède (représentée par \(\gamma_w\)) au poids saturé total.
Le sol "pèse" donc 9.81 kN/m³ une fois immergé.
2. Calcul du Gradient Critique (\(i_c\)) :
Nous divisons ce poids déjaugé par le poids volumique de l'eau pour obtenir le gradient limite.
Interprétation : Le résultat est exactement 1.00. Cela signifie que pour soulever ce sol, l'eau doit perdre une quantité d'énergie (charge) égale à la hauteur de sol traversée.
✅ Conclusion de l'étape
Le gradient critique de ce sable est de 1.00. C'est la valeur maximale théorique que le sol peut supporter sans se fluidifier. Cette valeur servira de référence pour le calcul de sécurité final.
La valeur de 1.00 est canonique. Elle s'explique par le fait que \(G_s \approx 2.65\) et \(e \approx 0.65\), ce qui mène mathématiquement à ce résultat rond pour \(\gamma_w \approx 10\). C'est un ordre de grandeur parfaitement validé.
Ce gradient est une limite théorique ultime. Il ne faut JAMAIS dimensionner un ouvrage pour fonctionner à \(i = i_c\). Dès que l'on approche cette valeur, des déformations importantes et des remaniements du sol se produisent avant même la rupture totale.
🎯 Objectif Scientifique
Après avoir déterminé la résistance du sol (Question 2), nous devons maintenant évaluer l'intensité de l'agression hydraulique. Le gradient hydraulique réel \(i\) mesure la perte d'énergie de l'eau par mètre parcouru le long du rideau de palplanches. C'est cette valeur qui va déterminer si la force d'écoulement est suffisante pour déstabiliser le fond de fouille.
📚 Référentiel
Loi de Darcy (1856) Hydraulique des SolsL'eau cherche toujours le chemin le plus court pour dissiper son énergie potentielle. Elle part du niveau haut (+H) et contourne l'obstacle (le rideau) pour remonter au niveau bas (0). Le gradient n'est pas constant dans tout le massif de sol : il est maximal là où les lignes de courant se resserrent, c'est-à-dire juste à la pointe de la fiche et à la remontée le long du rideau côté fouille. Pour cet exercice, nous utilisons une approche simplifiée conservatrice en calculant un gradient moyen sur la ligne de courant la plus dangereuse (celle qui longe l'acier).
Selon Bernoulli appliqué aux milieux poreux, la charge totale \(h\) diminue le long du trajet de l'eau à cause des frottements visqueux. Le gradient \(i\) est la dérivée de la charge par rapport à la distance : \(i = -dh/dl\). Dans une approche discrète simplifiée, c'est le rapport entre la différence de hauteur d'eau (potentiel) et la longueur du parcours.
Le gradient est le rapport de la perte de charge totale sur la longueur du chemin de drainage :
Avec \(\Delta H\) la différence de niveau d'eau et \(L_{\text{flow}}\) la longueur du chemin de percolation.
Étape 1 : Hypothèses & Données
| Paramètre | Valeur |
|---|---|
| Charge hydraulique \(H\) | 4.50 m |
| Fiche du rideau \(D\) | 3.50 m |
| Formule Longueur \(L_{\text{flow}}\) (Approximation) | \(D + D/2\) |
L'approximation \(L = 2D\) ou \(L = D + D/2\) est souvent utilisée en pré-dimensionnement. En réalité, le gradient de sortie exact dépend de la largeur de la fouille et de l'épaisseur de la couche perméable. Une simulation numérique par éléments finis donnerait une valeur plus précise, mais cette méthode manuelle est excellente pour une première approche sécuritaire.
Étape 2 : Calculs Détaillés
1. Estimation de la longueur de percolation (\(L_{\text{flow}}\)) :
On estime la longueur efficace que l'eau doit parcourir pour dissiper sa charge lors de la remontée critique. L'approximation considère la descente et une partie de la remontée sous la fiche.
2. Calcul du Gradient Réel (\(i\)) :
On divise la charge motrice (4.50m) par la longueur résistante (5.25m).
Interprétation : Le gradient obtenu est de 0.857. Cela signifie que pour chaque mètre parcouru, l'eau perd 85.7 cm de charge. C'est une force d'entraînement très intense.
✅ Conclusion de l'étape
Le gradient réel est de 0.86. Cela signifie que l'eau exerce une pression dynamique très proche de la limite de résistance du sol (1.00). L'ouvrage est très sollicité.
Le gradient est élevé (> 0.5), ce qui est typique des fouilles profondes dans des nappes hautes. La valeur est inférieure à 1, donc théoriquement stable, mais elle est dangereusement proche de l'unité.
Ce calcul suppose un sol parfaitement homogène. La présence d'une lentille plus perméable localement pourrait concentrer le flux et augmenter drastiquement ce gradient réel par endroits.
🎯 Objectif Scientifique
C'est l'étape de décision finale et critique. Nous devons comparer la capacité de résistance du sol (calculée en Q2) à l'action déstabilisatrice de l'eau (calculée en Q3). Le rapport de ces deux valeurs nous donnera le Coefficient de Sécurité (\(F_s\)). Ce chiffre unique synthétise le risque et dictera si le chantier peut ouvrir ou si des travaux de confortement sont indispensables.
📚 Référentiel
Eurocode 7 - Vérification HYD (Hydraulic Heave) Recommandations professionnelles (Clouterre / Écrans)En ingénierie structurelle, des coefficients de sécurité de 1.35 ou 1.50 sont courants. En géotechnique hydraulique, l'incertitude est bien plus grande : on ne voit pas ce qui se passe sous terre, la perméabilité peut varier d'un facteur 10 à 100 sur quelques mètres. C'est pourquoi les normes exigent des marges de sécurité importantes. Un coefficient calculé juste au-dessus de 1.0 signifie en réalité un risque de ruine quasi-certain à cause des imperfections locales.
Le facteur de sécurité \(F_s\) est défini comme le rapport entre les forces résistantes (gravité) et les forces motrices (écoulement). Dans notre approche en gradients, cela revient simplement à : \(F_s = i_c / i\). Selon l'Eurocode 7 (approche DA1/combinaisons), ou les méthodes traditionnelles, la valeur minimale requise pour le Renard est généralement fixée entre 1.5 et 2.0.
Étape 1 : Récapitulatif des Gradients
| Type | Variable | Valeur |
|---|---|---|
| Gradient Critique (Résistance) | \(i_c\) | 1.00 |
| Gradient Réel (Action) | \(i\) | 0.86 |
Si \(F_s < 1.5\) : Ouvrage considéré comme instable ou à risque inacceptable.
Étape 2 : Calcul du Coefficient de Sécurité
1. Ratio Résistance / Action :
Nous effectuons le rapport des deux gradients.
2. Résultat numérique :
Interprétation : Nous avons 17% de "marge" théorique avant la catastrophe. C'est très peu.
✅ Interprétation & Conclusion
Le coefficient de sécurité obtenu est de 1.17. Bien que mathématiquement supérieur à 1 (ce qui implique que le sol ne devrait pas se soulever instantanément selon la théorie parfaite), cette valeur est très inférieure au seuil réglementaire de 1.5 requis par l'Eurocode 7 pour les états limites ultimes hydrauliques.
Un \(F_s\) de 1.17 signifie que si le niveau d'eau extérieur monte de seulement quelques dizaines de centimètres (crue, marée) ou si la densité du sol est localement un peu plus faible, la rupture est immédiate. L'ouvrage est dans une zone de fragilité critique.
Le dimensionnement actuel n'est pas acceptable. Il est impératif de modifier le projet pour augmenter la sécurité. Les solutions techniques possibles sont :
- Augmenter la fiche \(D\) des palplanches (pour augmenter \(L_{\text{flow}}\)).
- Réaliser un bouchon de sol injecté ou de béton immergé en fond de fouille.
- Rabattre la nappe extérieure (pompage) pour diminuer \(H\).
📄 Livrable Final (Note de Calculs EXE)
| Ind. | Date | Objet de la modification | Rédacteur |
|---|---|---|---|
| A | 24/10/2024 | Vérification renard hydraulique - Version Initiale | Ing. T. Martin |
- Eurocode 7 - Calcul Géotechnique
- DTU 14.1 - Travaux de cuvelage
| Charge Hydraulique (H) | 4.50 m |
| Fiche Palplanche (D) | 3.50 m |
| Poids vol. saturé sol | 19.62 kN/m³ |
Analyse du risque de boulance en fond de fouille.
L'Ingénieur Stagiaire
L'Expert Senior
Laisser un commentaire