INTERACTION SOL STRUCTURE

1-INTRODUCTION

La dynamique de la structure permet d'évaluer les efforts, provenant des forces d'inertie
développées dans la structure lorsqu'elle est soumise à un mouvement de son support, qui sont
exercés par la structure sur sa fondation. De plus la dynamique des sols permet d’étudier les
mouvements du sol support lorsque les ondes sismiques se propagent dans celui-ci avant
d'atteindre l'ouvrage dont on cherche à étudier la réponse. La question se pose de savoir
comment ces deux phénomènes interagissent et dans quelle mesure le mouvement du support est
affecté par la réponse de la structure, dont la réponse sera elle même modifiée par le mouvement
du support. Le terme générique regroupant l'étude de ces phénomènes est désigné dans la
littérature sous le nom d'interaction sol-structure.
L’interaction action sol-structure traduit littéralement la modification du mouvement du sol ou
de la structure lors d’un séisme du fait de la présence de l’autre composant.
Cette interaction est bien entendu plus importante suivant :
 La nature du sol et ses caractéristiques
 Les caractéristiques de l’ouvrage
 Le mode de fondation
De nombreux aspects nécessitent une étude approfondie pour résoudre complètement un
problème d’interaction sol-structure. Ce sont :
 La définition de l’aléa sismique et du mouvement résultant
 L’étude du comportement du sol sous chargement cyclique
 L’évaluation de la réponse du sol en champ libre
 L’évaluation de la réponse dynamique des structures.
La figure1 illustre l'aspect fondamental de l'interaction; cet aspect est présenté ici dans le cas
d'une fondation sur pieux, partiellement enterrée dans le sol, mais les conclusions restent
applicables à tout type de fondation. Loin de la fondation, dans une région dénommée le champ
libre, les couches de sol sont traversées par des ondes sismiques dont la nature peut être
complexe : on y rencontre des ondes de volume, compression (P) et cisaillement (S), des ondes
de surface (Rayleigh, Love, Stoneley).
La nature des ondes est dictée par les caractéristiques de la source sismique mais également par
la géométrie et les caractéristiques mécaniques des terrains traversés.
Si l'on s'intéresse au mouvement de la fondation, les déformations du sol sont transmises à celle-
ci et engendrent un mouvement de la superstructure; même en l'absence de superstructure le
mouvement de la fondation est différent du mouvement du champ libre du fait des différences de

1
rigidité entre la fondation et le sol encaissant: le champ d'ondes incident est réfléchi et diffracté
par la fondation et donc modifie le mouvement total du sol au voisinage de celle-ci. Ce
phénomène est connu sous le nom d'interaction cinématique.
Par ailleurs, le mouvement induit sur la fondation développe des oscillations de la superstructure
et donc donne naissance à des forces d'inertie qui sont retransmises à la fondation sous forme de
forces et de moments. Ce phénomène est connu sous le nom d'interaction inertielle.
De toute évidence, le dimensionnement de la fondation doit tenir compte de ces deux
composantes de l'interaction.
Généralement, à tort, le terme interaction sol-structure ne désigne dans l'esprit des ingénieurs que
la part inertielle; il convient de garder à l'esprit que l'interaction cinématique peut dans certaines
configurations être significative, même si parfois elle peut être négligée.

Figure 1: Effet de l'interaction sol-structure sur un ouvrage

2
2 ILLUSTRATION DE L'EFFET DE L'INTERACTION SOL-STRUCTURE

2.1 Modèle analogique simplifiée

L'influence de l'interaction sol-structure sur la réponse d'un ouvrage peut être illustrée à l'aide du
modèle analogique de la figure2.
La structure est assimilée à une masse et un ressort, placés à une hauteur h au-dessus de la
fondation.
La liaison entre la structure et la fondation est réalisée par une barre rigide.
La fondation repose sur le sol et son interaction avec celui-ci est modélisée par le biais des
fonctions d'impédance qui seront définies plus loin. On admettra pour l'instant que les fonctions
d'impédance, c'est à dire les réactions exercées par le sol sur la fondation, peuvent être
représentées par un ensemble de ressorts et d'amortisseurs indépendants de la fréquence.
L'amortisseur rend théoriquement compte à la fois de l'amortissement radiatif, c'est à dire de la
dissipation d'énergie par les ondes s'éloignant de la fondation, et de l'amortissement propre du
matériau "sol", appelé amortissement matériel ou interne. Dans un souci de simplification de la
présentation, on supposera que l'amortissement matériel est négligeable devant l'amortissement
radiatif (comportement élastique du sol), ce qui est valide pour un milieu homogène et des
sollicitations sismiques d'amplitudes faibles à moyennes.

Figure2 : Modèle simplifié d'interaction sol-structure

Les équations d'équilibre dynamique du système s'obtiennent aisément à partir des équations de
Lagrange en prenant comme variables généralisées qi :
Le système de la figure 2 possède 3 degrés de liberté :
 le déplacement horizontal u de la masse m,
 le déplacement horizontal u0 de la fondation,
 la rotation θ de la fondation autour d'un axe horizontal.

3
Il est soumis à un déplacement horizontal du sol support, harmonique de pulsation ω et
d'amplitude ug.
Les équations d'équilibre dynamique du système s'obtiennent aisément à partir des équations de
Lagrange en prenant comme variables généralisées qi :
 q1= u , le déplacement relatif de la masse m par rapport à A
 q2=uo , le déplacement de la fondation
 q3= θ , la rotation de la fondation
On a la relation évidente entre le déplacement absolu u t de la masse m et les variables
précédentes :

(1)

Désignant par T l'énergie cinétique totale :

(2)

par V l'énergie potentielle :

(3)

et par δW le travail des forces non conservatives (forces d'amortissement) :

(4)

les équations de Lagrange s'écrivent :

(5)

soit avec les notations suivantes en tenant compte des relations entre accélération, vitesse et
déplacement

(6)

et en introduisant les pourcentages d'amortissement critique :

4
(7)

(8)

En introduisant les notations suivantes :

(9)

et en éliminant uo et θ entre les trois équations précédentes, il vient :

(10)

Tenant compte du fait que :

l'équation précédente devient :

(11)

Considérons maintenant un oscillateur simple à 1 degré de liberté de même masse m, de

pulsation propre et d'amortissement soumis à un déplacement harmonique de

pulsation ω à sa base (cas de la structure encastrée à sa base). La réponse harmonique de cet
oscillateur est :

(12)

L'oscillateur équivalent aura même réponse que la structure de la figure 2 si les équations
suivantes sont vérifiées :

5
(13a)

(13b)

(13c)

Les équations précédentes sont obtenues en égalant partie réelle et partie imaginaire des
équations (11) et (12) et pour (13c) en se plaçant à résonance Il résulte des équations
(13) que l'interaction sol-structure a pour effet :
 de diminuer la pulsation propre ωs de la structure base encastrée

 d'augmenter l'amortissement du système par rapport à la structure base encastrée

 de diminuer la sollicitation incidente effective à la base de la structure

Les conclusions précédentes sont visualisées sur la figure3 qui présente pour une fondation
circulaire reposant sur un semi espace élastique homogène, les variations relatives

en fonction des paramètres adimensionnels :

(14)

où r est le rayon de la fondation et ρ, Vs la masse volumique et la célérité des ondes S dans le

sol.
La figure 3 met clairement en évidence que l'influence de l'interaction sol-structure est d'autant
plus prononcée que le sol de fondation est mou (s croissant) ou que la structure est massive ( m
croissant).

6
7
 Figure.3 : Influence de l'interaction sol - structure
2.2 Exemple
L'exemple présenté ci dessous est celui de la centrale d'Humbolt Bay en Californie qui a subi en
1977 un séisme important au cours duquel les accélérations en champ libre et en quelques points
de l'ouvrage ont été enregistrées. La figure 4 montre les spectres de réponse calculés à partir des
enregistrements en champ libre et dans la structure à la même cote.

L'interaction sol structure est particulièrement marquée pour cet ouvrage partiellement enterré
dans le sol comme le montre la forte atténuation des accélérations spectrales pour les fréquences
supérieures au Hertz.

8
 Figure 4 : Spectres de réponse de la centrale d'Humbolt Bay

3 FORMULATION DE L'INTERACTION SOL STRUCTURE

Une analyse complète du phénomène d’interaction sol-structure doit prendre en compte :

 La variation des caractéristiques du sol avec la profondeur

 Le comportement non linéaire du sol

 Le caractère tridimensionnel du problème

 Le schéma complexe de propagation des ondes qui engendrent le mouvcement

 L’interaction avec les structures avoisinantes

Dans l’état actuel des connaissances, des simplifications doivent être Effectuées pour permettre
le traitement du problème. Ces modifications concernent surtout la schématisation du sol de
fondation (comportement linéaire du sol), la nature des ondes incidentes( ondes planes, de
cisaillement, de compression, de Rayleigh)…..

Avant d'examiner les différentes méthodes de prise en compte de l'interaction sol-structure, il est
utile de formuler de façon générale le problème. Cette formulation est orientée vers un traitement
par éléments finis du phénomène d'interaction. En effet, la complexité du problème est telle que
le recours aux méthodes numériques est pratiquement inévitable. Les équations du mouvement
sont obtenues par référence à la figure 5 qui schématise un ensemble sol-structure.
Désignant par M, C, K, les matrices de masse, amortissement et raideur du système, l'équation
du mouvement s'écrit :

9
(15)

Comme la source du mouvement (foyer du séisme) n'est généralement pas incluse dans le
modèle, le vecteur de charge Qf n'a de valeurs non nulles que sur la frontière extérieure du
modèle.
En l'absence de structure, l'équation du mouvement du champ libre est analogue de par sa forme
à l'équation (15); les indices f désignant les matrices masses, amortissement et raideur relative au
seul champ libre, cette équation s'écrit :

(16)

Posant :

(17)

l'équation (17) définit le déplacement d'interaction i U qui satisfait l'équation :

(18)

avec :

(19)

Le vecteur de charge Qi est déterminé à partir des déplacements en champ libre. Pour les
systèmes linéaires, on a alors le théorème de superposition illustré sur la figure 5: le problème
d'interaction est décomposé en la somme d'un problème de réponse du sol en champ libre (éq 16)
et d'un problème source (éq 18) où les forces appliquées Qi n'ont de composantes non nulles
qu'aux noeuds communs à la structure et au sol. Ce dernier problème est, par essence, analogue à
un problème de vibration de machine. Le déplacement total pour le problème d'interaction est
alors donné par l'équation (17).

10
 Figure 5 : Décomposition du problème d'interaction sol-structure

L'équation (19) met clairement en évidence le fait qu'il y a interaction dès qu'il y a différence de
masse ou de raideur entre le sol et la structure.
Supprimons pour simplifier le terme d'amortissement dans cette équation et restreignons le
problème à celui d'une structure posée à la surface du sol et soumise à la propagation verticale
d'ondes de volume (cisaillement ou compression). Dans ces conditions, en champ libre, tous les
points de la surface du sol sont animés d'un même mouvement. Si la fondation de l'ouvrage est
infiniment rigide, le dernier terme de l'équation (19) s'annule; le vecteur de charge Qi se réduit à :

(20)

Les forces Qi appliquées à la base de la structure engendrent un mouvement du support,

équivalent à un champ de forces d'inertie dans la structure. Par suite, l'interaction ne résulte que
des forces d'inertie développées dans cette structure. On lui donne le nom d'interaction inertielle,
dont l'effet a été illustré sur l'exemple du paragraphe 2.1.
A l'opposé, considérons une structure enterrée dont la masse est nulle hors du sol et égale (en
valeur et répartition) à celle du sol pour la partie en terre. Les forces Qi ont alors pour expression

(21)

Elles ne résultent que de la différence de raideur pour la partie en terre, entre le sol et la
structure. Même sans différence de masse, il y a interaction; on lui donne le nom d'interaction
cinématique. Elle résulte de la raideur de la fondation qui l'empêche de suivre les mouvements
imposés par le sol. On a vu précédemment qu'elle était rigoureusement nulle pour certains cas;
elle peut être faible dans d'autres cas (fondations sur pieux souples) ou très importante (structure
raide fortement contreventée enterrée dans le sol).

11
Dans le cas le plus général, l'interaction résulte d'une interaction inertielle et d'une interaction
cinématique.
La figure 5 et les raisonnements précédents illustrent les deux grandes méthodes d'approche de
l'interaction sol-structure. La figure 5a correspond aux méthodes globales dont la solution est
obtenue par résolution directe de l'équation (15). Elles ne font intervenir aucune notion de
superposition et sont donc théoriquement adaptées aux problèmes non linéaires. Alternativement,
les méthodes de sous-structures s'appuient sur la décomposition des figures 5b-5c, ou sur des
décompositions analogues, pour résoudre le problème par étapes. Ces méthodes ne sont bien
entendu applicables qu'aux problèmes linéaires, justifiables de superposition.
3.1 Exemple d’interaction cinématique
On a signalé au paragraphe précédent qu'il existe des situations pour lesquelles l'interaction
cinématique est rigoureusement nulle et l'interaction sol-structure se réduit à sa part inertielle.
Cette situation est exceptionnelle; il existe par contre d'autres situations plus courantes pour
lesquelles l'interaction cinématique, sans être rigoureusement nulle, est négligeable. Prenons
l'exemple d'un pieu de fondation dans une couche de sol; si la raideur du pieu vis à vis du sol
encaissant est faible (cas d'un "spaghetti" dans un sol raide), on comprend que la raideur du pieu
ne modifie en rien le mouvement du sol par rapport à celui du champ libre. Pour illustrer ces
propos, considérons l'exemple de la figure6.

Figure 6 : Exemple d'interaction cinématique

La couche de sol est soumise à un mouvement harmonique à sa base, d'amplitude U0 et de

pulsation ω. On s'intéresse au mouvement résultant de la tête du pieu, Up, et plus exactement au
rapport de ce mouvement au mouvement du sol en champ libre  Uf : soit Iu ce rapport.
La figure 7 présente la variation de Iu en fonction d'un paramètre adimensionnel, proportionnel
au rapport des raideurs du pieu et du sol et inversement proportionnel à l'élancement du pieu.
Plus ce rapport est grand plus le pieu est raide vis à vis du sol encaissant.
L'examen de la figure 7 montre que pour des raideurs faibles du pieu, l'interaction cinématique
est effectivement nulle: le coefficient IU est égal à 1.0 et le pieu suit le mouvement du sol sans le
modifier. Lorsque la raideur augmente, la présence du pieu se ressent et le coefficient Iu devient
dif férent de 1.0: le pieu ne suit plus le mouvement du sol.

12
 Figure 7 : Interaction cinématique pour un pieu isolé

On notera qu'il existe des valeurs du paramètre FC pour lesquelles le coefficient Iu est supérieur
à 1.0 : le pieu amplifie le mouvement du sol en champ libre par suite de sa mise en résonance.

4 METHODES DE PRISE EN COMPTE DE L'INTERACTION SOL STRUCTURE

Parmi les méthodes de prise en compte de l'interaction sol structure on peut distinguer les
méthodes globales, qui résolvent, comme leur nom l'indique, le problème global de la figure5a,
et celles qui s'appuient sur une décomposition du système en sous systèmes, dans l'esprit de la
figure 5b-5c; ces méthodes sont désignées sous le nom générique de méthodes de sous structure.
4.1 Méthode globale
La méthode consiste à résoudre en une seule étape l'équation dynamique :

(22)

Où
U représente le vecteur des déplacements relatifs du système par rapport à l'assise
..
I un vecteur unité, donnant la direction de la sollicitation u g

M , K , C les matrices de masse, de raideur et d'amortissement du système.

Généralement la technique de résolution est basée sur la méthode des éléments finis qui présente
la plus grande flexibilité et permet en outre la prise en compte de phénomènes non linéaires tel
que le comportement anélastique des matériaux, le décollement ou le glissement des fondations.
La problématique de la résolution par une méthode globale est illustrée sur la figure 8 :
 Le mouvement sismique de dimensionnement est connu (spécifié) à la surface du sol, en
champ libre;
 Le mouvement est calculé à la base inférieure du modèle, choisie à une profondeur
suffisante pour que la présence d'une structure en surface n'affecte pas ce mouvement;
cette étape est connue sous le nom de déconvolution du mouvement sismique; 

13
  Le mouvement déconvolué est imposé uniformément à la base du système sol structure et
la réponse est calculée par résolution de l'équation (22).

Figure 8 : Schématisation d'un problème d'interaction sol structure en éléments finis

La difficulté de la résolution par éléments finis des problèmes dynamiques d'interaction sol
structure réside dans le traitement des conditions aux limites. On a vu au chapitre 9 qu'une onde
heurtant une surface libre était totalement réfléchie. Dans le cas présent, les limites du modèle
sont introduites de manière artificielle dans celui-ci et la réflexion des ondes à ces interfaces
contribue à maintenir l'énergie qu'elles transportent à l'intérieur du modèle, alors que dans la
réalité elles la transportent à l'infini (au moins partiellement si d'autres réflexions interviennent) :
ce transport d'énergie est désigné sous le vocable amortissement géométrique.
Il est possible d'annihiler ce phénomène en introduisant des frontières absorbantes dont le rôle
est de rétablir les conditions de contraintes et de déplacements à la frontière; dans le cas des
poutres ces frontières sont représentées par des amortisseurs dont les caractéristiques sont
fonction de celles du milieu extérieur au modèle. Pour le milieu tridimensionnel des frontières
analogues ont été développées; ces éléments ne constituent une solution exacte au problème que
dans le domaine fréquentiel; dans le domaine temporel elles ne représentent qu'une solution
approchée.
Une autre manière d'annihiler ce phénomène est de coupler le modèle éléments finis avec des
éléments infinis.

14
Un des avantages principaux des méthodes globales est leur capacité à prendre en compte en
considération :
 Le caractère tridimensionnel du problème
 Des comportements non linéaires du sol ou de la structure
 Une configuration complexe des couches de sol
 Une variation des caractéristiques du sol avec la profondeur
 La promixité de structure avoisinante

4.2 Méthode de sous-structure

Cette catégorie de méthodes fait appel au principe de superposition. Ainsi il existent plusieurs
méthode de sous-structure suivant la manière dont s’effectuent la sous structuration : On
distinguent les méthodes dites de freontière et les méthodes dites de volume.
4.2.1 Méthode de frontière
L'idée consiste à analyser le problème d'interaction sol structure en plusieurs étapes successives,
chacune des étapes étant réputée plus facile à résoudre que le problème global. Pour des raisons
évidentes, les sous structures envisagées sont constituées d'une part par le sol et d'autre part par
la structure, comme indiqué sur la figure 9. On écrit les équations d'équilibre de chaque sous
système, puis les conditions de compatibilité à l'interface : continuité du déplacement et du
vecteur contrainte.

Figure 9 : Schématisation d'une méthode de sous structure

Pour des raisons qui apparaîtront naturelles dans la suite, le problème est traité ci-dessous dans le
domaine fréquentiel; par ailleurs pour ne pas alourdir la présentation, on omettra les termes
d'amortissement dans les équations d'équilibre. Il s'ensuit que les grandeurs, par exemple le
déplacement U , s'expriment sous la forme de leur transformée de Fourier :

15
(23)

Dans la suite on notera les transformées de Fourier avec le symbole ˜ au dessus de la variable et
on omettra l'indice n. Chaque équation écrite ci après s'applique donc à toute harmonique de la
décomposition de Fourier.
Dans ces conditions les équations des différents sous systèmes s'écrivent :
 Structure

(24)

 Sol

(25)

Dans les équations ci dessus, comme indiqué sur la figure 9, on a désigné par l'indice B les
déplacements de la structure, l'indice F ceux de l'interface sol-structure et l'indice S ceux du sol.
De plus, les équations ont été partitionnées de façon à isoler dans chaque sous système les
équations faisant intervenir les degrés de liberté communs: ces équations sont repérées par les
indices FB lorsqu'elles appartiennent au sous système structure et par les indices FS lorsqu'elles
appartiennent au sous système sol. Enfin on a tenu compte des conditions de compatibilité
exprimées ci dessus.
Considérons maintenant le cas du sous système sol en l'absence de la structure; son équation
d'équilibre s'écrit de façon analogue à l'équation (25) :

(26)

Dans l'expression (26) u* représente le déplacement d'interaction cinématique. Appelant

ui le déplacement d'interaction défini par :

(27)

par soustraction des équations (25) et (26) il vient :

16
(28)

Le système d équations (28) peut être utilisé pour éliminer tous les degrés de liberté qui
n'appartiennent pas à l'interface sol-structure. Ce processus appelé condensation permet de
relier les déplacements de l'interface aux réactions en ces noeuds :

(29)

Dans l'équation (29) la matrice s'appelle la matrice d'impédance de la fondation. On notera

que les quantités intervenant dans (29), et en particulier , dépendent de la fréquence.

A ce stade, la réaction du sol est inconnue; on élimine cette grandeur en reportant (29)
dans (24) et en tenant compte de (27); l'équation d'équilibre de la structure devient :

(30)

4.2.1.1 Frontière rigide

Si la fondation est rigide le champ des déplacements de celle-ci peut s'exprimer en termes d'un
mouvement de corps rigide défini par rapport à un point quelconque, par exemple son centre :

(31)

où est la matrice de transformation et le vecteur des déplacements et rotations du centre

de la fondation.

De même, les forces nodales sont reliées aux forces et moments appliqués à la
fondation en son centre:

(32)

Avec ces définitions l'équation (30) prend la forme :

17
(33)

La matrice : représente la matrice d'impédance de la fondation rigide. Notant

que l'équation précédente représente le mouvement d'une structure

reliée à un support par la matrice d'impédance et soumise à un mouvement de ce support

défini par , qui rappelons le représente le mouvement d'interaction cinématique.

Figure 10 : Théorème de superposition

Ainsi, dans l'hypothèse d'une fondation rigide, il est pertinent de scinder le problème global
en trois sous-problèmes :
 détermination du mouvement d'une fondation rigide sans masse soumise à la sollicitation
sismique; cette étape représente la solution de l'équation (26);
 détermination de la matrice d'impédance de la fondation (éq 28);
 calcul de la réponse dynamique de la structure reliée à la matrice d'impédance et sujette à
son support au mouvement d'interaction cinématique (éq 33).
Dans la mesure où la fondation est parfaitement rigide, cette démarche est rigoureusement
identique à celle conduisant à la résolution du système global en une étape (éq 22). Cette

18
décomposition est connue sous le nom de théorème de superposition de Kausel et est illustrée sur
la figure 10. Son intérêt apparaît clairement s'il est possible de simplifier une des trois étapes du
calcul.
Le problème de diffraction (étape a) existe toujours sauf pour une structure fondée en surface et
soumise à la propagation verticale d'ondes de volume; dans ce cas, la résolution de la première
étape est identique à celle de la réponse d'un profil de sol en champ libre puisque l'interaction
cinématique est nulle.
La solution à la deuxième étape peut être évitée, pour certaines configurations, en utilisant les
résultats de fonctions d'impédances publiés dans la littérature.
La troisième étape est, en tout état de cause, indispensable; elle est cependant plus simple et plus
familière aux ingénieurs car elle procède de l'analyse dynamique classique des structures.
4.2.1.1 Frontière flexible
Dans le cas général, la fondation n’est pas infiniment rigide ; le problème à résoudre est celui
plus complexe définit par l’équation 30. La structure, reliée aux impédances du sol, est soumise à
un système de sollicitations multiples à chacune de ses liaisons avec le sol. Dans le principe, la
résolution des trois étapes que précédemment est nécessaire :
 Diffraction du champ libre
 Impédance
 Réponse de la structure
Mais chacune est plus complexe que dans le cas d’une fondation rigide. En particulier, on ne
peut plus avoir recours à des solutions publiées pour la détermination des impédances
4.2.2 Méthode de volume
La méthode de volume permet d’éviter la résolution du problème de diffraction. L’interaction
entre le sol et la structure est prise en compte, non seulement à l’interface entre le sol et la
structure , mais à tous les éléments (nœuds) de la structure sous la surface du sol, comme cela
apparaît sur la figure 11
Les équations du problème peuvent alors être formulée comme pour les méthodes de sous-
structure envisagées précédemment en procédant exactement de la même façon.
La solution du problème d’ interaction est réduite alors aux trois étapes suivantes :
 Détermination de la réponse de champ libre. Aucun problème de diffraction n’existe à ce
stade, si les couches de sol sont horizontales. Le problème peut alors être résolu, pour
n’importe quel type d’ondes, à l’aide des méthodes classiques.
 Détermination de la matrice d’impédance. Cette dernière contient plus de terme que celle
des méthodes de frontière 
 Détermination classique de la réponse dynamique de la structure. Il s’agit là aussi d’un
problème classique

19
 = +

Structure moins
Système global Fondation sol excavé

Figure 11 : Méthode de volume

Cette méthode à été incorporée dans un programme d’interaction de calcul SASSI(1981) qui
permet de prendre le caractère tridimensionnel du problème pour une stratigraphie régulière
autour de la structure.
4.3 Méthode hybride
Exception faite du cas où il existe une solution analytique, ou une solution déjà publiée pour une
configuration proche de celle étudiée, la résolution du problème d’impédance constitue une étape
longue et coûteuse dans la méthode de sous structure. La solution au problème d’impédance
nécessite la résolution, pour un grand nombre de fréquence, de l’équation (29)
Pour s’affranchir de cette étape, l’idée consiste à séparer le sol en un champ libre et un champ
lointain (figure 12)
Le champ lointain est modélisé à l’aide d’une matrice d’impédance. En d’autre termes, le
concept de sous-structures est étendu de façon à inclure le champ proche dans le modèle de
structure.
Le problème consiste alors à définir les coefficients de la matrice d’impédance du champ
lointain.
Un deuxième difficulté réside dans l’obtention, le long de la frontière séparant le champ proche
du champ lointain, de la solution du problème de diffraction du champ incident.

20
 Figure 12 Méthode Hybride
4.4 Comparaison des différentes méthodes
Un des arguments majeurs avancés pour l’utilisation des méthodes de sous-strucures, de
préférences aux méthodes globales, est leur facilité d’emploi alliée à un coût moindre. Cela est
certainement vrai pour des ouvrages fondés superficiellement, soumis à la propagation d’ondes
de volume. Dans ce cas, on a vu que l’interaction cinématique est nulle, supprimant l’étape du
problème de diffraction, et que des solutions publiées existent pour le problème d’impédance.. Si
par ailleurs le mouvement est défini à la surface du sol la seule étape reste l’analyse dynamique
de la structure.
Par ailleurs les méthodes de sous-structures présentent l’avantage de permettre certaines
modifications sans qu’il soit nécessaire de reprendre toute l’analyse. Par exemple une
modification au niveau de la structure ne requiert qu’une nouvelle analyse dynamique de la
structure. Ou encore un modification des caractéristiques du séisme de projet permet de
conserver la solution du problème d‘impédance.
Mais dés que l’ouvrage est partiellement ou totalement enterré, la résolution du problème de
diffraction devient complexe. Des méthodes simplifiées pour la prise en compte de l’interaction
cinématique, ont été proposée dans le cas de fondation rigide. Mais dans le cas le plus général les
méthodes de sous structures deviennent moins compétitives que les méthodes globales. Ces
dernières permettent en outre de traiter les problèmes non linéaires.
La figure 13 résume les différentes étapes nécessaires à la résolution d’un problème d’interaction
sol structure par les méthodes décrites ci-dessus. On constate que si les étapes de réponse de site
en champ libre et d’analyse dynamique de la structure sont commune à toutes les méthodes, la
réponse globale et la méthode de volume des sous-structures sont les plus attrayantes, exception
faite du cas des fondations superficielles sur un sol relativement homogène.

21
 Figure 13 Comparaison des méthodes d’interaction sol-structure

5 IMPEDANCE D'UNE FONDATION SUPERFICIELLE

Pour illustrer la notion d'impédance d'une fondation, grandeur essentielle pour le calcul sismique
d'une structure par une méthode de sous structure, considérons le cas simple d'une fondation
circulaire rigide reposant à la surface d'un semi espace élastique, homogène et isotrope (figure
14)

Figure 14 : Impédance d'une fondation superficielle circulaire

22
La fondation est caractérisée par son rayon r0, et le semi espace par sa masse volumique ρ et
deux paramètres de comportement, par exemple le module de cisaillement G et le coefficient de
Poisson ν; on utilisera également la célérité des ondes de VS pour définir le semi espace.

Par définition (éq. 29), l'impédance de la fondation est égale à la réaction exercée sur la
fondation sans masse lorsqu'elle est soumise à des déplacements harmoniques unitaires dirigés
suivant l'un quelconque de ses degrés de liberté. La fondation étant sans masse l'impédance
représente également le quotient d'une force directement appliquée à la fondation (qui est égale à
la réaction du sol) par le déplacement résultant. Une fondation rigide possédant six degrés de
liberté, la matrice d'impédance a pour dimension 6x6; si la fondation est de forme
quelconque, les différents degrés de liberté sont couplés et la matrice d'impédance est pleine. Si
la fondation possède des symétries, certains des termes de couplage (termes hors diagonale)
disparaissent. Dans le cas de la fondation circulaire, il existe 4 degrés de liberté: les translations
horizontale et verticale, la rotation autour d'un axe horizontal et la rotation autour d'un axe
vertical. La translation verticale et la rotation autour d'un axe vertical sont totalement découplés
des autres degrés de liberté; par contre, en toute théorie, le déplacement horizontal et la rotation
autour d'un axe horizontal sont couplés entre eux; cependant, pour la fondation superficielle le
terme de couplage est négligeable et il est licite de considérer que la matrice d'impédance est une
matrice diagonale de dimension 4x4.
Chaque terme de la matrice représente donc le quotient de la force appliquée par le déplacement
résultant suivant le même degré de liberté. Dans la suite du paragraphe on raisonnera donc sur un
des termes de la matrice que l'on notera K, par exemple celui correspondant au déplacement
vertical, sachant que les considérations qui sont développées sont également applicables aux
autres termes; on dénommera ce terme par le vocable générique impédance.
Avant d'étudier l'impédance de la fondation, il est utile de faire le parallèle avec l'oscillateur à un
degré de liberté pour mieux comprendre la structure de l'impédance. Considérons un tel
oscillateur soumis à une force harmonique Peiωt. Le déplacement résultant vaut :

(34)

Par définition l'impédance de l'oscillateur simple s'écrit :

(35)

dans laquelle on a utilisé les notations usuelles de l'oscillateur simple : ω n et ξ désigne la

pulsation propre et le pourcentage d'amortissement critique.

23
L'examen de la relation (35) montre que l'impédance est le produit :
 d'un terme correspondant à la raideur statique k;
 d'un terme qui représente la partie dynamique; ce terme comporte une partie réelle et une
partie imaginaire qui provient du fait que le déplacement est déphasé par rapport à la
force appliquée.
Ce déphasage est lié à la dissipation d'énergie du système.
On notera également que la partie réelle de l'impédance peut devenir négative à haute fréquence.
Par analogie avec l'équation (35) on montre que l'impédance de la fondation peut s'écrire sous la
forme générale :

(36)

Avec

(37)

Elle se compose d'un terme multiplicatif qui est la raideur statique (raideur à fréquence nulle) de
la fondation KS et d'une partie représentant la contribution dynamique. Cette contribution
dynamique comporte une partie réelle et une partie imaginaire. Les coefficients a0, k1 et c1 sont
sans dimension et dépendent de la pulsation ω.
Il est instructif de donner une interprétation physique à la notion d'impédance. Par définition
l'équation (36) donne la réaction R du sol exercée sur la fondation, soit :

(38)

Tenant compte du fait que pour une sollicitation harmonique :

(39)

l'équation (38) peut s'écrire :

(40)

La réaction du sol se compose donc de deux termes, proportionnels au déplacement et à la

vitesse de la fondation. On peut donner à ces termes la signification physique d'un ressort et d'un
amortisseur de caractéristiques :

24
(41)

qui dépendent de la pulsation w. Par ailleurs on notera qu'il existe une dissipation d'énergie dans
le système, représentée par l'amortisseur, bien que le milieu sur lequel repose la fondation soit
élastique, donc non dissipatif. Cette dissipation provient du transport de l'énergie dans le milieu
par les ondes issues du mouvement de la fondation: il s'agit de l'amortissement géométrique,
La figure 15 présente les fonctions d'impédance de la fondation superficielle circulaire sur le
semi espace élastique, homogène, isotrope. Seules sont présentées les parties dynamiques de
l'impédance pour chacun des degrés de liberté de la fondation. On vérifie comme indiqué par les
relations (41) que les termes k1 et c1 dépendent plus ou moins fortement de ω et comme on l'a
signalé pour l'oscillateur simple, qu'ils peuvent devenir négatifs.

Figure 15 : Impédance de la fondation circulaire superficielle

25
Si les termes de l'impédance étaient constants la modélisation de celle ci dans le calcul
dynamique de la structure, dernière étape du théorème de superposition de la figure 10, serait
immédiate: on la représenterait pour chaque degré de liberté par un ressort et un amortisseur (éq
41). Dans le cas général cette représentation n'est pas possible et oblige pour la résolution des
problèmes d'interaction sol structure à se placer dans le domaine fréquentiel, comme cela a été
fait pour l'établissement du théorème de superposition.
Il existe cependant des configurations pour lesquelles l'approximation par un ressort et un
amortisseur (un pour chaque degré de liberté) indépendants de la fréquence est possible: c'est le
cas en particulier pour la fondation circulaire étudiée dans ce paragraphe. En comparant la
réponse "exacte" de cette fondation à celle d'une fondation sur ressort et amortisseur, il est
possible de reproduire de façon satisfaisante le réponse de la fondation par un choix approprié
des caractéristiques du ressort et de l'amortisseur. A titre d'exemple le tableau 1 donne un choix
possible des paramètres et la figure 16 compare la solution exacte à la réponse obtenue avec le
modèle simplifié dans le cas de la sollicitation verticale.

Tableau 1 : Coefficients d'impédance pour la fondation superficielle circulaire

26
 Figure 16 : Réponse d'une fondation circulaire à une force verticale
Dans les expressions du tableau 1 les quantités I désignent les moments d'inertie massiques en
rotation de la fondation autour des axes vertical et horizontal.
La simplification ci dessus ne doit pas faire perdre de vue que pour des configurations plus
complexes, soit en termes de stratification du sol, soit de géométrie de la fondation,
l'approximation par des ressorts et amortisseurs indépendants de la fréquence n'est plus possible.
La figure 17 illustre ces propos en montrant la variation des parties réelle et imaginaire de
l'impédance verticale d'une fondation circulaire superficielle fondée en surface d'une couche
d'épaisseur limitée surmontant un substratum infiniment rigide.

Figure 17 : Impédance d'une fondation circulaire sur un mono-couche

27
On note la variation très erratique des deux termes avec la fréquence et surtout le fait que pour
des fréquences faibles, inférieures à une valeur f 0, c1 est nulle: il n'existe aucune dissipation
d'énergie dans le système. Ce phénomène s'explique en notant que les ondes diffractées par la
fondation sont totalement réfléchies lorsqu'elles heurtent le substratum rigide; l'énergie n'est pas
transportée à l'infini dans le semi espace mais est renvoyée vers la fondation. Ce n'est que
lorsque f devient supérieure à f0 que la dissipation d'énergie prend place; on peut montrer que f0
correspond à la fréquence propre de vibration de la couche (VP/4H pour le mode vertical, VS/4H
pour le mode horizontal) et qu'il y a alors création d'ondes de surface se propageant
horizontalement dans la couche de sol supérieure. Il est bien évident que des variations telles que
celles de la figure 17 ne peuvent être prises en compte par des ressorts et amortisseurs
indépendants de la fréquence. Il est cependant possible, au prix d'une légère complication, de
rendre compte de ces variations avec des modèles rhéologiques simples. En reprenant l'exemple
de l'oscillateur à un degré de liberté on note que la variation avec la fréquence de la partie réelle
de l'impédance est induite par la présence de la masse de l'oscillateur. L'idée est donc d'adjoindre
au modèle de ressort et amortisseur une masse additionnelle dont le rôle est de traduire cette
variation. Certains auteurs ont justifié l'introduction d'une masse additionnelle en arguant d'une
masse de sol "attachée" à la fondation et vibrant en phase avec elle; il ne s'agit en fait que d'un
artifice mathématique permettant de mieux traduire le comportement global de la fondation. La
figure 18 présente deux exemples de modèles rhéologiques simples avec masses additionnelles
donnant naissance à des termes d'impédance variables avec la fréquence. Les paramètres de
divers éléments entrant dans le modèle sont calés de façon à reproduire la variation exacte de
l'impédance. La figure 18 montre un tel exemple dans le cas d'une fondation réelle d'un ouvrage:
la courbe en trait continu représente la solution exacte et les symboles les prédictions du modèle.
On notera en particulier que le modèle permet de représenter la valeur négative de la partie réelle
de l'impédance au delà de la fréquence de 1.5Hz.

Figure 18 : Modèles rhéologiques simplifiés pour l'impédance

28
Figure 18 : Exemple d'application du modèle de la figure 17
Partie réelle de l'impédance verticale

29