Infographie technique de l'interaction sol-structure avec modélisation du coefficient de Winkler et tassement élastique sous fondation.

L'analyse de l'Interaction Sol-Structure (ISS) constitue une étape cruciale dans le dimensionnement moderne des structures de génie civil. Ignorer la déformabilité du sol sous une fondation peut conduire à une sous-estimation dramatique des moments fléchissants dans la structure ou à une mauvaise évaluation des tassements différentiels.

Le modèle classique et le plus largement implémenté dans les logiciels de calcul de structures (tels que Robot Structural Analysis, SCIA Engineer ou ETABS) reste le modèle de fondation élastique de Winkler. Ce guide technique détaille les fondements théoriques, les limites méthodologiques et le processus de calcul rigoureux pour définir ce coefficient élastique.

1. Les Hypothèses Fondamentales du Modèle de Winkler

Proposé par Emil Winkler en 1867, ce modèle assimile le massif de sol situé sous une fondation à un système de ressorts élastiques linéaires, indépendants et continus.

Hypothèses physiques et mathématiques :

  • Linéarité élastique : La pression de contact p(x, y) en un point de l'interface sol-fondation est directement proportionnelle au déplacement vertical (tassement) s(x, y) en ce même point.
  • Indépendance des ressorts : Le tassement en un point A n’induit aucun déplacement en un point B adjacent, tant que la charge n'est pas directement appliquée en B. Le sol n'a donc pas de résistance au cisaillement mutuel entre ses mailles.
  • Comportement unilatéral ou bilatéral : En pratique numérique, les ressorts sont souvent paramétrés pour ne travailler qu’en compression (non-linéarité géométrique pour simuler le décollement de la fondation).

La loi constitutive fondamentale s'exprime par la formule :

p = ks · s

Où :

  • p est la pression de contact du sol (kN/m2 ou kPa)
  • s est le tassement vertical (m)
  • ks est le coefficient de réaction du sol ou coefficient élastique de Winkler (kN/m3)

2. Méthodologies pour Déterminer le Coefficient de Winkler (ks)

Le coefficient ks n'est pas une propriété intrinsèque du sol. Il dépend à la fois des caractéristiques du sol (module d'élasticité, coefficient de Poisson) et de la géométrie de la fondation (largeur, rigidité de la structure). Plusieurs approches permettent de l'isoler.

Méthode A : L'essai de chargement à la plaque (Norme NF P 94-117-1)

L'essai fournit directement un coefficient de réaction initial pour une plaque normalisée (souvent circulaire de diamètre 30 cm, noté k30). Pour l'extrapoler à une fondation réelle de largeur B, on applique les formules empiriques de Terzaghi :

  • Pour les sols cohérents (argiles) :
    ks = k30 · (0.3 / B)
  • Pour les sols pulvérulents (sables) :
    ks = k30 · [(B + 0.3) / (2B)]2

(Où B est exprimé en mètres)

Méthode B : L'approche analytique élastique (Formule de Vesic)

Pour s'affranchir des limites des essais à la plaque, Vesic (1961) a proposé une formulation rigoureuse basée sur la mécanique des milieux continus, reliant ks au module de déformation du sol (Es) et à la rigidité de la fondation (Ef If) :

ks = (1.3 / B) · 12√[(Es · B4) / (Ef · If)] · [Es / (1 - νs2)]

Où :

  • Es = Module de déformation élastique du sol (kPa)
  • νs = Coefficient de Poisson du sol
  • Ef = Module d'élasticité du matériau de la fondation (kPa)
  • If = Moment d'inertie de la section transversale de la fondation par unité de longueur (m4/m)
  • B = Largeur de la fondation (m)

Dans la majorité des applications pratiques de calcul de radiers ou de semelles filantes, le radical d'ordre 12 est proche de 1, ce qui permet de simplifier la formule :

ks ≈ Es / [B · (1 - νs2)]

3. Schéma Conceptuel et Équation Différentielle de l'Interaction

Le comportement cinématique d'une fondation superficielle (poutre sur sol élastique) est régi par l'équation différentielle du quatrième ordre d'Euler-Bernoulli, modifiée par le terme de réaction de Winkler :

Ef If · (d4s(x) / dx4) + q(x) - p(x) = 0

En remplaçant p(x) par ks · B · s(x), on obtient :

Ef If · (d4s(x) / dx4) + ks · B · s(x) = q(x)

Où q(x) représente la charge externe linéaire appliquée à la structure.

4. Exemple Pratique d’Application et de Calcul Réel

Énoncé du problème

Soit une semelle filante en béton armé supportant un voile porteur :

  • Largeur de la semelle (B) = 1.50 m
  • Épaisseur de la semelle (h) = 0.40 m
  • Matériau : Béton C25/30 → Ef = 32 x 106 kPa
  • Caractéristiques du sol (Sable moyennement dense) :
    • Module de déformation (Es) = 25 000 kPa
    • Coefficient de Poisson (νs) = 0.3

Étape 1 : Calcul des propriétés géométriques de la fondation

On calcule le moment d'inertie If par mètre linéaire de la semelle filante :

If = (1 · h3) / 12 = (1 · 0.403) / 12 = 0.005333 m4/m

Rigidité à la flexion de la fondation :

Ef If = 32 x 106 x 0.005333 = 170 656 kN·m2/m

Étape 2 : Calcul du coefficient de réaction ks par la formule de Vesic

Appliquons la formule complète de Vesic pour vérifier la sensibilité du terme de rigidité relative :

Terme sous la racine :
(Es · B4) / (Ef · If) = (25 000 × 1.504) / 170 656 = 126 562.5 / 170 656 ≈ 0.7416

La racine 12ème de 0.7416 donne environ 0.9754.

Calculons maintenant la valeur globale de ks :

ks = (1.3 / 1.50) × 0.9754 × [25 000 / (1 - 0.32)]

ks = 0.8667 × 0.9754 × (25 000 / 0.91)

ks = 23 225.3 kN/m3

Étape 3 : Interprétation et intégration MEF (Modèle Éléments Finis)

Si vous modélisez cette semelle filante sous un logiciel d'éléments finis à l'aide d'éléments filaires (poutres), la raideur linéaire du ressort à introduire (Klinéaire) s'obtient en multipliant ks par la largeur de contact B :

Klinéaire = ks · B = 23 225.3 × 1.50 = 34 837.95 kN/m/m

Signification physique : Une charge répartie uniforme de 34.84 kN/m provoquera un enfoncement élastique uniforme de la semelle de exactement 1 mm dans le sol.

5. Synthèse des Limites du Modèle et Recommandations

Bien que robuste et simple à implémenter, le modèle de Winkler présente l'inconvénient de surestimer les pressions au centre des radiers rigides et de sous-estimer le tassement de la cuvette périphérique. Le sol réel transfère les contraintes au-delà de la zone chargée.

Pour les projets d'envergure, il est recommandé de coupler cette approche avec une méthode de pseudo-couplage continu, en faisant varier le coefficient ks par zones (coefficients plus élevés sur les bords et les angles de la structure) afin de mieux simuler le comportement du milieu élastique continu de Boussinesq



DÉCOUVREZ NOS FORMATIONS SUR LES LOGICIELS DE GÉNIE CIVIL ET BTP

formations sur les logiciels de génie civil et btp