Analyse de l’Interaction Sol-Structure
📝 Contexte Géologique et Enjeux de l'Étude
📄 OBJET : RAPPORT D'ALERTE GÉOTECHNIQUE - MISSION G3
- Projet : Silo Céréalier "Alpha"
- Localisation : Extension du Port de Commerce
- Phase : Étude d'Exécution (EXE) - Statut : URGENT
Suite aux récents sondages pressiométriques réalisés sur la parcelle d'extension du port, la stratigraphie du site a été entièrement réévaluée. Le terrain naturel, particulièrement complexe, se décompose de la manière suivante :
- 🟤 En surface : Une couche de remblais compactés servant de plateforme de travail.
- 🟠 En profondeur : Une épaisse couche d'argile molle saturée, reconnue comme hautement compressible.
- 🪨 À la base : Un substratum rocheux de type calcaire, réputé totalement indéformable.
Le maître d'ouvrage prévoit l'édification d'un imposant silo de stockage, reposant sur un radier circulaire rigide en béton armé. La nature très déformable de la couche d'argile soulève de sérieuses inquiétudes concernant l'interaction sol-structure.
Un tassement excessif ou asymétrique sous le poids cyclique des céréales pourrait provoquer :
- Une ruine par poinçonnement du sol de fondation (Défaut de capacité portante).
- Une fissuration critique de la coque structurelle (Désordre d'usage grave).
- Ouvrage mitoyen : Le site jouxte, à seulement \( 15 \text{ m} \), un quai centenaire en maçonnerie. Toute déformation du sol risque d'induire des poussées parasites critiques.
- Hétérogénéité des données : L'ingénieur devra faire preuve de discernement pour extraire la valeur caractéristique de cohésion à partir des profils hétérogènes du rapport de laboratoire.
Votre objectif est de garantir la stabilité de l'ouvrage selon les normes Eurocode 7. Vous devrez vérifier la capacité portante à l'ELU et évaluer précisément les tassements de consolidation à l'ELS.
Pour faciliter les calculs manuels, le fluage secondaire de l'argile et la souplesse exacte de la superstructure seront négligés. Le radier est modélisé comme infiniment rigide.
Voici la synthèse des paramètres de calcul, extraits de la campagne d'investigation et de la descente de charges.
📚 Référentiels Normatifs Appliqués
[NF EN 1997-1] Eurocode 7 [NF EN 1992-1-1] Eurocode 2🟤 CARACTÉRISTIQUES DU SOL / MILIEU
-
Couche 1 : Remblais Compactés (de 0,0 à -2,0 m)
- Épaisseur de la couche : \( h_{\text{remblai}} = \mathbf{2{,}0} \text{ m} \)
- Poids volumique apparent : \( \gamma_{\text{remblai}} = \mathbf{19{,}0} \text{ kN/m}^3 \)
- Angle de frottement interne : \( \varphi'_{\text{remblai}} = \mathbf{30}^{\circ} \)
-
Couche 2 : Argile Molle Compressible (de -2,0 à -10,0 m)Extrait Brut du Rapport de Laboratoire (Sondage SC-03) : Le toit de la couche est détecté à -2,0 m, base sur substratum à -10,0 m.
- Poids volumique saturé : \( \gamma_{\text{sat}} = \mathbf{18{,}0} \text{ kN/m}^3 \)
- Cohésion non drainée retenue : \( c_{\text{u}} = \mathbf{40{,}0} \text{ kPa} \) (Plage scissomètre 35-48 kPa)
- Indice des vides initial : \( e_{\text{0}} = \mathbf{0{,}9} \) (Sol normalement consolidé)
- Indice de compression : \( C_{\text{c}} = \mathbf{0{,}25} \)
-
Conditions Hydrauliques (Nappe Phréatique)
- Niveau initial stabilisé : Au toit de l'argile (\( z = \mathbf{-2{,}0} \text{ m} \))
- Poids volumique de l'eau : \( \gamma_{\text{w}} = \mathbf{10{,}0} \text{ kN/m}^3 \)
⬇️ CHARGEMENTS ET GÉOMÉTRIE
-
Ouvrage (Le Silo et son Radier)
- Diamètre de la fondation : \( D = \mathbf{15{,}0} \text{ m} \)
- Profondeur d'encastrement : \( D_{\text{f}} = \mathbf{2{,}0} \text{ m} \) (Radier fondé directement sur le toit de l'argile)
- Action permanente (Poids propre) : \( G = \mathbf{5~000} \text{ kN} \)
- Action d'exploitation (Blé) : \( Q = \mathbf{10~000} \text{ kN} \)
3. Séquence de Résolution Géotechnique
Afin de justifier l'ouvrage selon les règles de l'art, l'étude est découpée en quatre phases analytiques.
Question 1 : Évaluation de l'État Contraint Initial
Avant tout apport de charge externe, il est indispensable de modéliser le sol au repos. 👉 Déterminez la contrainte verticale totale, la pression interstitielle et la contrainte verticale effective au centre exact de la couche d'argile molle (à \( z = -6{,}0 \text{ m} \)).
Question 2 : Vérification de la Capacité Portante (ELU)
L'ouvrage doit présenter une sécurité absolue contre la rupture du sol. 👉 Calculez la capacité portante limite du sol à court terme (comportement non drainé). Vérifiez ensuite le critère de portance à l'ELU sous la combinaison fondamentale d'actions (\( 1{,}35 \cdot G + 1{,}5 \cdot Q \)).
La contrainte de rupture \( q_{\text{rupture}} \) pour une fondation sur sol fin saturé (\( \varphi_{\text{u}} = 0 \)) est donnée par : \( q_{\text{rupture}} = c_{\text{u}} \cdot N_{\text{c}} \cdot s_{\text{c}} + q \cdot N_{\text{q}} \)
| Paramètre Géotechnique | Symbole | Valeur / Formule |
|---|---|---|
| Facteur de portance (Cohésion) | \( N_{\text{c}} \) | \( \pi + 2 \approx \mathbf{5{,}14} \) |
| Facteur de portance (Surcharge) | \( N_{\text{q}} \) | \( \mathbf{1{,}0} \) |
| Facteur de portance (Poids) | \( N_{\gamma} \) | \( \mathbf{0{,}0} \) |
| Facteur de forme (Cohésion) | \( s_{\text{c}} \) | \( \mathbf{1{,}2} \) |
| Facteur de forme (Surcharge) | \( s_{\text{q}} \) | \( \mathbf{1{,}0} \) |
Question 3 : Estimation du Tassement de Consolidation (ELS)
L'argile subira une compression lente et différée dans le temps. 👉 Évaluez l'accroissement de contrainte \( \Delta \sigma_{\text{v}} \) au centre de la couche d'argile via la méthode de diffusion géométrique (2:1). Calculez ensuite le tassement primaire de consolidation \( s_{\text{c}} \) sous la combinaison ELS quasi-permanente (\( G + Q \)).
Question 4 (Finale) : Instabilité Hydraulique et Flottabilité
🚨 ALÉA DE CHANTIER DÉCLARÉ : Une violente tempête a provoqué une submersion marine exceptionnelle sur l'extension du port. L'eau a non seulement saturé les remblais, mais la nappe s'est établie temporairement au-dessus du niveau du sol, avec une charge d'eau mesurée à \( z = \mathbf{+1{,}0} \text{ m} \).
👉 Travail d'expertise demandé :
1. Évaluez l'impact physique direct de cet aléa sur la contrainte effective du sol de couverture.
2. Vérifiez la stabilité globale de l'ouvrage (silo totalement vide) vis-à-vis du risque de soulèvement par flottabilité.
3. L'ouvrage étant en péril, proposez une parade technique pérenne.
La stabilité est assurée si l'action permanente déstabilisante de l'eau (\( G_{\text{dst}} \)) reste inférieure à l'action permanente stabilisante de l'ouvrage (\( G_{\text{stb}} \)), pondérées par leurs coefficients partiels.
| Nature de l'action permanente | Symbole | Valeur (UPL) |
|---|---|---|
| Action déstabilisante (Eau) | \( \gamma_{\text{G,dst}} \) | \( \mathbf{1{,}0} \) |
| Action stabilisante (Poids) | \( \gamma_{\text{G,stb}} \) | \( \mathbf{0{,}9} \) |
Analyse de l’Interaction Sol-Structure
🎯 Objectif Technique
L'objectif central de cette première étape est d'établir avec précision le profil des contraintes initiales in-situ (notre point zéro géotechnique) à la profondeur névralgique de \( z = -6{,}0 \text{ m} \). Sans la connaissance absolue de cet état naturel de référence, toute estimation du tassement futur sous les charges du silo serait scientifiquement caduque.
En mécanique des sols, l'eau interstitielle ne porte pas les charges structurelles de manière pérenne ; seul le squelette granulaire s'en charge. Mon raisonnement pas-à-pas se décompose ainsi :
- 🪨 Évaluer la charge globale : Calculer le poids géostatique total descendant (terres superposées).
- 💧 Mesurer l'opposition : Calculer la force de flottaison ascendante générée par la nappe phréatique.
- ⚖️ Trouver l'équilibre : Soustraire l'eau du poids total pour isoler l'effort réellement encaissé par la matrice argileuse.
C'est la loi fondatrice de la géotechnique moderne. Elle stipule que toutes les déformations mesurables d'un sol saturé sont exclusivement gouvernées par les modifications de la contrainte effective.
Afin d'obtenir la contrainte effective, nous devons d'abord exprimer la contrainte verticale totale, qui correspond à la somme des poids des différentes couches de sol superposées :
- \( \sigma_{\text{v}} \) : Contrainte verticale totale au point d'étude (en kPa).
- \( \gamma_{\text{i}} \) : Poids volumique de la couche \( \text{i} \) (en kN/m³).
- \( h_{\text{i}} \) : Épaisseur traversée de la couche \( \text{i} \) (en m).
Ensuite, nous modélisons la pression interstitielle de l'eau présente dans les pores de l'argile (que nous identifierons visuellement en bleu) :
- \( \color{#0284c7}{u} \) : Pression interstitielle de l'eau (en kPa).
- \( \color{#0284c7}{\gamma_{\text{w}}} \) : Poids volumique de l'eau (en kN/m³).
- \( \color{#0284c7}{h_{\text{w}}} \) : Hauteur de la colonne d'eau au-dessus du point d'étude (en m).
Enfin, nous manipulons ces deux grandeurs par la relation de Terzaghi en soustrayant l'eau de la contrainte totale pour isoler notre inconnue (identifiée en brun terre) :
- \( \color{#92400e}{\sigma'_{\text{v0}}} \) : Contrainte verticale effective initiale (en kPa).
Origine rigoureuse des paramètres pour le calcul à \( z = -6{,}0 \text{ m} \) :
Couche 1 (Remblais secs) :- \( \gamma_{\text{remblai}} = \mathbf{19{,}0} \text{ kN/m}^3 \) ➔ Source : Énoncé, Section 2 (Mesure in-situ).
- \( h_{\text{remblai}} = \mathbf{2{,}0} \text{ m} \) ➔ Source : Énoncé, Section 2 (Épaisseur de 0,0 à -2,0 m).
- \( \gamma_{\text{sat}} = \mathbf{18{,}0} \text{ kN/m}^3 \) ➔ Source : Énoncé, Section 2 (Rapport SC-03).
- \( h_{\text{argile}} = \mathbf{4{,}0} \text{ m} \) ➔ Source : Déduction géométrique (Toit à -2,0 m, point d'étude à -6,0 m. L'épaisseur traversée est \( 6{,}0 - 2{,}0 = 4{,}0 \text{ m} \)).
- \( \color{#0284c7}{\gamma_{\text{w}}} = \mathbf{10{,}0} \text{ kN/m}^3 \) ➔ Source : Énoncé, Section 2 (Valeur normative).
- \( \color{#0284c7}{h_{\text{w}}} = \mathbf{4{,}0} \text{ m} \) ➔ Source : Déduction géométrique (Nappe à -2,0 m, point d'étude à -6,0 m).
Dessinez systématiquement l'axe des profondeurs (Z) orienté vers le bas ! Une erreur classique d'étudiant consiste à utiliser la profondeur absolue (\( 6{,}0 \text{ m} \)) pour calculer la pression de l'eau, en omettant que la nappe phréatique ne démarre qu'à la cote \( -2{,}0 \text{ m} \).
Exécution de la Note de Calculs
💡 Démonstration des Unités : En multipliant un poids volumique (\( \text{kN/m}^3 \)) par une épaisseur (\( \text{m} \)), les mètres s'annulent partiellement pour donner une pression surfacique (\( \text{kN/m}^2 \)). Or, par définition physique du Système International : \( 1 \text{ kN/m}^2 = 1 \text{ kPa} \).
Étape A Contrainte verticale totale (\( \sigma_{\text{v}} \))Nous remplaçons les termes littéraux par leurs valeurs et sommons l'impact du remblai et de l'argile pour obtenir la charge globale :
Nous multiplions la hauteur de la colonne d'eau par son poids volumique pour connaître la force de flottaison :
Nous retranchons l'eau de la contrainte totale pour isoler la charge supportée par les grains solides :
INTERPRÉTATION : Avant même le premier coup de pelleteuse, le squelette granulaire situé au cœur de l'argile encaisse une pression stabilisée de \( 70 \text{ kPa} \). C'est notre pression de référence de base ; l'apparition de tout tassement futur dépendra du dépassement de ce seuil spécifique.
L'ordre de grandeur est rigoureusement validé. Une règle de chantier empirique rapide stipule qu'un sol déjaugé (immergé sous la nappe) pèse environ \( 10 \text{ kPa/m} \). Ici, nous cumulons \( 2 \text{ m} \) de sol sec (\( \approx 40 \text{ kPa} \)) et \( 4 \text{ m} \) de sol immergé (\( \approx 40 \text{ kPa} \)). Le total intuitif valide pleinement notre résultat théorique rigoureux.
Attention à ne jamais appliquer le poids volumique saturé comme s'il s'agissait d'un sol sec ! Omettre de soustraire la sous-pression de l'eau mènerait à surévaluer la résistance initiale du sol de plus de 50%, ce qui ruinerait la justesse de l'ensemble du modèle géotechnique.
❓ Pourquoi ne pas utiliser directement la formule du sol déjaugé (\( \gamma' = \gamma_{\text{sat}} - \gamma_{\text{w}} \)) pour gagner du temps ?
Bien que cette méthode soit mathématiquement juste et plus rapide, dissocier systématiquement les contraintes totales de l'hydraulique est une exigence absolue en bureau d'études. Cela permet de sécuriser la note de calcul face aux redoutables fluctuations futures de la nappe phréatique (comme nous le traiterons dans l'aléa de la Question 4).
🎯 Objectif Technique
L'enjeu primordial est de garantir que le radier ne cèdera pas par poinçonnement dans l'argile molle lors d'une sollicitation extrême. Nous devons certifier que la pression transmise par le silo reste inférieure à la capacité ultime du terrain selon les directives de l'Eurocode 7.
L'argile molle est un sol fin extrêmement peu perméable. Lors du remplissage rapide du silo, l'eau prisonnière dans les pores n'aura pas le temps de s'évacuer. Mon plan d'action stratégique consiste à :
- 📉 Agréger les charges : Estimer la charge totale maximale descendante en appliquant les pondérations sécuritaires de l'Eurocode.
- 🎯 Répartir l'effort : Transformer cet effort global (en kN) en contrainte surfacique (en kPa) via l'aire géométrique du radier.
- 🛡️ Vérifier le sol : Calculer la résistance maximale du sol en condition non drainée (à court terme), où seule la cohésion interne retient la structure.
L'État Limite Ultime (ELU) de portance exige formellement que l'action de calcul appliquée à la base de la fondation (\( q_{\text{ELU}} \)) soit strictement inférieure à la contrainte de rupture géotechnique du sol (\( q_{\text{rupture}} \)).
La pondération sécuritaire des actions structurelles s'exprime en sommant les charges permanentes et variables multipliées par leurs coefficients respectifs :
- \( V_{\text{ELU}} \) : Effort vertical pondéré (en kN).
- \( \gamma_{\text{G}} \), \( \gamma_{\text{Q}} \) : Coefficients partiels de sécurité.
- \( G \), \( Q \) : Actions permanentes et variables.
Pour obtenir la contrainte, il convient de manipuler l'effort en le divisant par l'aire du cercle (\( A = \frac{\pi \cdot D^2}{4} \)) :
- \( \color{#dc2626}{q_{\text{ELU}}} \) : Contrainte appliquée sous la base (Pression déstabilisante en rouge).
- \( D \) : Diamètre de la fondation (en m).
Enfin, la résistance de rupture en condition non drainée (Eurocode 7, Annexe D) est dictée par la théorie de la plasticité (combinaison de la cohésion et de la butée latérale) :
- \( \color{#16a34a}{q_{\text{rupture}}} \) : Capacité portante limite du sol (Résistance en vert).
- \( c_{\text{u}} \) : Cohésion non drainée de l'argile (en kPa).
- \( N_{\text{c}}, N_{\text{q}} \) : Facteurs de portance adimensionnels.
- \( s_{\text{c}} \) : Facteur de forme géométrique.
- \( q \) : Surcharge géostatique latérale liée à l'encastrement (\( q = \gamma_{\text{remblai}} \cdot D_{\text{f}} \)).
D'où proviennent nos paramètres de calcul à l'ELU ?
- \( G = \mathbf{5~000} \text{ kN} \) et \( Q = \mathbf{10~000} \text{ kN} \) ➔ Source : Énoncé, Section 2 (Descente de charges).
- \( \gamma_{\text{G}} = \mathbf{1{,}35} \) et \( \gamma_{\text{Q}} = \mathbf{1{,}50} \) ➔ Source : Norme Eurocode 7 (Combinaison fondamentale STR/GEO).
- \( D = \mathbf{15{,}0} \text{ m} \) ➔ Source : Énoncé, Section 2.
- \( c_{\text{u}} = \mathbf{40{,}0} \text{ kPa} \) ➔ Source : Énoncé, Section 2 (Valeur prudente du scissomètre).
- \( q = \mathbf{38{,}0} \text{ kPa} \) ➔ Source : Calcul de la surcharge d'encastrement (\( 19{,}0 \text{ kN/m}^3 \times 2{,}0 \text{ m} \)).
- \( N_{\text{c}} = \mathbf{5{,}14} \), \( s_{\text{c}} = \mathbf{1{,}2} \), \( N_{\text{q}} = \mathbf{1{,}0} \) ➔ Source : Énoncé, Section 3 (Tableau normatif D.1).
Utilisez bien le poids total du remblai (\( 19{,}0 \text{ kN/m}^3 \)) pour calculer la butée latérale \( q \), car la nappe phréatique est située sous ce remblai. Si l'eau affleurait en surface, il aurait été impératif d'utiliser le poids déjaugé !
Exécution de la Note de Calculs
💡 Démonstration des Unités : Lors du calcul de la contrainte \( q_{\text{ELU}} \), la division d'une force surfacique en kilonewtons (\( \text{kN} \)) par une surface en mètres carrés (\( \text{m}^2 \)) forme nativement le rapport \( \text{kN/m}^2 \), soit très exactement la définition du kilopascal (\( \text{kPa} \)).
Étape A Calcul de la charge totale pondérée ELU (\( V_{\text{ELU}} \))Nous injectons les charges avec leurs coefficients de sécurité pour modéliser le cas d'usage le plus défavorable :
Nous divisons cet effort extrême par l'aire du radier pour trouver la pression unitaire sous la fondation :
Nous appliquons la théorie de l'Eurocode pour trouver la pression maximale absolue que le sol peut encaisser :
INTERPRÉTATION : La pression extrême et destructrice injectée par le silo dans le sol est de \( 123 \text{ kPa} \). L'argile, confinée par les remblais et sa propre cohésion interne, peut encaisser jusqu'à \( 284 \text{ kPa} \). Le radier joue parfaitement son rôle de "raquette à neige" empêchant l'enfoncement par poinçonnement.
Le critère de rupture ELU est brillamment validé avec un facteur de sécurité d'environ 2,3. L'argile molle ne cédera pas à court terme. Cependant, une portance saine ne protège en rien contre les redoutables déformations lentes à long terme, ce qui motive l'analyse urgente de la Question 3.
Dans les sols argileux étudiés à court terme (\( \varphi_{\text{u}} = 0 \)), le terme de portance lié à la largeur de la fondation (\( N_{\gamma} \)) est totalement nul. Élargir un radier n'augmente absolument pas la résistance unitaire intrinsèque du sol, cela augmente uniquement son aire géométrique de répartition.
❓ Pourquoi utilise-t-on la contrainte totale pour la surcharge d'encastrement (\( q \)) ?
À court terme (en conditions non drainées), il est impératif de raisonner en contraintes totales car la surpression interstitielle engendrée par le poids de la construction n'a pas encore eu le temps physique de se dissiper. C'est l'approche en "bloc" imposée par l'Eurocode.
🎯 Objectif Technique
L'objectif est d'anticiper la compression oedométrique irréversible de l'argile sur plusieurs décennies. Nous devons calculer le tassement induit par l'expulsion extrêmement lente de l'eau interstitielle, sous la combinaison d'utilisation quotidienne du silo (État Limite de Service - ELS).
Mon raisonnement logique et séquentiel est le suivant :
- Poids réel : Calculer la charge descendante de service (sans pondérations artificielles de sécurité).
- Effet de profondeur : Estimer l'"écho" de cette charge au centre de la couche d'argile via la méthode géométrique 2:1. Plus l'onde de charge descend, plus elle s'étale et perd en intensité.
- Loi de déformation : Appliquer la loi logarithmique de Terzaghi pour transformer cette surcontrainte en tassement physique (en mètres).
Le sol étant attesté "normalement consolidé", sa déformation verticale est proportionnelle au logarithme du rapport entre l'état de contrainte final (\( \sigma'_{\text{v0}} + \Delta \sigma_{\text{v}} \)) et l'état initial au repos (\( \sigma'_{\text{v0}} \)), pondéré par sa signature de compressibilité (\( C_{\text{c}} \)).
L'évaluation ELS nécessite de sommer les charges structurelles brutes (sans majoration) :
- \( V_{\text{ELS}} \) : Effort vertical total en service (en kN).
La méthode de diffusion 2:1 modélise l'élargissement de la zone d'influence d'une largeur initiale \( D \), augmentée de la profondeur relative \( z' \) :
- \( \color{#dc2626}{\Delta \sigma_{\text{v}}} \) : Accroissement de contrainte transmis au point d'étude (Surcharge en rouge).
- \( z' \) : Profondeur relative mesurée sous la base du radier (en m).
La loi de tassement se formule par l'équation intégrée de Terzaghi :
- \( s_{\text{c}} \) : Tassement absolu de consolidation (en m).
- \( e_{\text{0}} \) : Indice des vides initial (adimensionnel).
- \( C_{\text{c}} \) : Indice de compression (adimensionnel).
Les données d'entrée du calcul de tassement :
- \( V_{\text{ELS}} = \mathbf{15~000} \text{ kN} \) ➔ Source : Énoncé, Section 2 (Somme non pondérée \( G + Q \), soit \( 5000 + 10000 \)).
- \( D = \mathbf{15{,}0} \text{ m} \) ➔ Source : Énoncé, Section 2.
- \( z' = \mathbf{4{,}0} \text{ m} \) ➔ Source : Déduction (Base du radier à -2,0 m, centre argile à -6,0 m. Distance \( = 4{,}0 \text{ m} \)).
- \( h_{\text{argile}} = \mathbf{8{,}0} \text{ m} \) ➔ Source : Énoncé, Section 2.
- \( e_{\text{0}} = \mathbf{0{,}9} \) et \( C_{\text{c}} = \mathbf{0{,}25} \) ➔ Source : Énoncé, Section 2 (Rapport Oedométrique).
- \( \color{#92400e}{\sigma'_{\text{v0}}} = \mathbf{70{,}0} \text{ kPa} \) ➔ Source : Résultat rigoureux de la Question 1.
Soyez extrêmement vigilant sur la variable \( z' \). Dans le théorème de diffusion géométrique, \( z' \) désigne la distance d'épanouissement sous la base de la fondation, et non la cote altimétrique absolue lue depuis la surface topographique !
Exécution de la Note de Calculs
💡 Démonstration Transparente des Unités de Tassement : Dans l'équation \( s_{\text{c}} \), les indices mécaniques (\( C_{\text{c}} \) et \( e_{\text{0}} \)) et le terme logarithmique (\( \log_{\text{10}} \)) sont de purs scalaires sans unité. Le résultat mathématique emprunte donc intégralement l'unité de la variable géométrique d'entrée \( h_{\text{argile}} \), soit des mètres (\( \text{m} \)). Une simple conversion en fin de calcul (\( \times 100 \)) donnera une lecture plus parlante en centimètres.
Étape A Accroissement de contrainte (\( \color{#dc2626}{\Delta \sigma_{\text{v}}} \))Nous injectons les valeurs dans le modèle géométrique : le diamètre d'influence de l'onde de charge passe de \( 15 \text{ m} \) (sous le radier) à \( 15 + 4 = 19 \text{ m} \) (au cœur de l'argile).
Nous appliquons la contrainte trouvée dans la loi de déformation logarithmique pour obtenir l'affaissement du sol :
Pour une lecture de chantier lisible, nous convertissons les mètres en centimètres (\( 0{,}257 \text{ m} \times 100 \text{ cm/m} \)) :
INTERPRÉTATION : Sous le tonnage du silo en service continu, l'argile molle va refermer ses porosités par consolidation hydrique. L'ouvrage s'enfoncera verticalement de près de 26 centimètres sur le long terme (généralement mesuré sur une décennie).
Un tassement de \( 25{,}7 \text{ cm} \) est qualifié de structurellement pathologique. Il exigera obligatoirement la mise en œuvre de raccords souples et flexibles pour la tuyauterie industrielle du silo, et conforte pleinement les craintes exprimées dans l'énoncé vis-à-vis du quai centenaire mitoyen.
L'écueil classique des étudiants est d'activer la touche du logarithme népérien (\( \ln \)) au lieu du logarithme décimal (\( \log_{\text{10}} \)). L'indice de compression \( C_{\text{c}} \) géotechnique est strictement établi sur une échelle logarithmique en base 10.
❓ Pourquoi utiliser une méthode géométrique 2:1 simplifiée plutôt que les équations rigoureuses de Boussinesq ?
La théorie d'élasticité linéaire de Boussinesq implique des calculs d'intégrales très lourds à gérer manuellement. La méthode 2:1 offre une approximation redoutable (souvent précise à 5% près), devenue le standard de l'industrie pour évaluer rapidement la descente de charge sous un radier homogène.
🎯 Objectif Technique face à l'imprévu
L'objectif est d'expertiser la stabilité globale de l'ouvrage face à une submersion marine exceptionnelle. Nous devons quantifier la perte d'efficacité du sol de remblai (déjaugeage) et prouver analytiquement si le silo vide risque d'être détruit sous l'effet de la poussée d'Archimède.
Mon raisonnement en temps de crise se fait en deux étapes :
- Étape Sol : L'inondation à \( z = +1{,}0 \text{ m} \) noie totalement les remblais. Ces derniers perdent leur poids apparent à cause de la poussée de l'eau, réduisant dangereusement leur pouvoir confinant.
- Étape Structure : Le radier ancré à \( -2{,}0 \text{ m} \) subit une pression ascendante féroce sur sa face inférieure (3 mètres d'eau). Si le silo est vide, seule la tôle s'oppose à ce navire en devenir. Il faut confronter ces deux forces vectorielles avec les pondérations réglementaires.
La norme (NF EN 1997-1) exige de démontrer de façon absolue que la résultante des actions déstabilisantes (l'eau) demeure strictement inférieure à la résultante des actions stabilisantes (le poids de la structure).
Nous calculons d'abord la chute de la contrainte effective du remblai noyé au moyen du poids volumique déjaugé :
- \( q' \) : Nouvelle contrainte effective d'encastrement (en kPa).
Nous établissons par une simple soustraction la charge hydraulique réelle (\( \color{#0284c7}{h_{\text{w}}} \)) frappant la sous-face du radier :
La force de flottaison est quantifiée en multipliant la pression d'eau par l'aire du radier circulaire (\( A = \frac{\pi \cdot D^2}{4} \)) :
- \( \color{#dc2626}{G_{\text{dst}}} \) : Action permanente déstabilisante de l'eau (Poussée ascendante en rouge).
Le critère normatif de stabilité UPL confronte l'eau contre la gravité en y appliquant les pondérations :
- \( \color{#16a34a}{G_{\text{stb}}} \) : Action permanente stabilisante (Poids mort en vert).
- \( \gamma_{\text{G,dst}}, \gamma_{\text{G,stb}} \) : Coefficients partiels (UPL).
Les entrées de crise UPL :
- \( \color{#16a34a}{G_{\text{stb}}} = \mathbf{5~000} \text{ kN} \) ➔ Source : Énoncé, Section 2 (Poids mort du silo métallique).
- \( D = \mathbf{15{,}0} \text{ m} \) ➔ Source : Énoncé, Section 2 (Géométrie).
- \( \color{#0284c7}{h_{\text{w}}} = \mathbf{3{,}0} \text{ m} \) ➔ Source : Déduction. Eau à \( +1{,}0 \text{ m} \), radier à \( -2{,}0 \text{ m} \). Différentiel de \( 1{,}0 - (-2{,}0) = 3{,}0 \text{ m} \).
- \( \gamma_{\text{G,dst}} = \mathbf{1{,}0} \) et \( \gamma_{\text{G,stb}} = \mathbf{0{,}9} \) ➔ Source : Énoncé, Section 3 (Tableau A.15 - Coefficients Eurocode).
Le coefficient minorateur de \( 0{,}9 \) sur le poids stabilisant permet d'absorber l'incertitude liée aux épaisseurs réelles de béton ou d'acier de l'ouvrage. En géotechnique de sécurité, on ne parie jamais sur 100% de la gravité pour contrer l'océan !
Exécution de la Note de Calculs et Solutions
💡 Démonstration des Unités : Pour \( G_{\text{dst}} \), nous multiplions la pression de l'eau (\( \text{kN/m}^3 \times \text{m} = \text{kN/m}^2 \)) par l'aire du radier (\( \text{m}^2 \)). Les \( \text{m}^2 \) s'annulent naturellement. Le produit est parfaitement homogène à une force pure en kilonewtons (\( \text{kN} \)).
Étape A Évaluation de la nouvelle contrainte effective du remblai (\( q' \))La submersion annule une grande partie du poids apparent du remblai, prouvant l'impact catastrophique de l'eau :
Nous insérons la pression de l'eau et l'aire de la fondation pour obtenir la force de soulèvement totale :
Nous appliquons la confrontation de l'Eurocode en pondérant les valeurs :
PROPOSITION TECHNIQUE DE PARADE : L'inéquation est violée, le silo est virtuellement arraché du sol. La parade technique d'urgence consiste à créer un ancrage pour clouer le radier au sol. Solution recommandée : Perforer le radier et sceller des micropieux ou des tirants d'ancrage dans le substratum rocheux sain (à \( -10{,}0 \text{ m} \)), travaillant en traction pure pour bloquer le soulèvement.
Cet événement climatique démontre qu'un dimensionnement exclusivement gravitaire (fondé sur le seul poids mort de l'ouvrage) était notoirement insuffisant pour ce site exposé. L'intégration de fondations profondes travaillant en traction va inévitablement bouleverser l'économie initiale du projet, mais demeure l'unique alternative de pérennité.
Le clouage du radier par micropieux va encastrer rigidement l'ouvrage dans le substratum rocheux. Conséquence fâcheuse : le tassement de \( 25{,}7 \text{ cm} \) validé à l'étape 3 sera totalement bloqué ! Le modèle mécanique passe brusquement d'un radier superficiel "souple" à un radier sur pieux "ultra-rigide". Les efforts de cisaillement internes dans le béton vont s'envoler, obligeant le pôle Génie Civil à reprendre l'intégralité du ferraillage.
❓ Pourquoi le sol de remblai pesant sur les bords du radier n'aide-t-il pas à maintenir l'ouvrage plaqué ?
Lors de la submersion, le remblai pesant sur les débords perd lui-même la moitié de son poids sous l'eau (comme démontré à l'Étape A). De surcroît, le modeste frottement latéral sol/béton généré sur seulement \( 2 \text{ m} \) d'encastrement est insignifiant face aux milliers de kilonewtons de la terrible poussée d'Archimède.
