Mémoire Thèse AdlF

UNIVERSITE DE NANTES
ECOLE DOCTORALE
« MECANIQUE, THERMIQUE ET GENIE CIVIL »
DE NANTES
Année 2001
Thèse de DOCTORAT
Discipline : Sciences pour l’Ingénieur

Spécialité : Génie Civil
Présentée et soutenue publiquement par
Alexandre de la Foye
Le 9 juillet 2001
à l’Ecole Centrale de Nantes
CALCUL DE LA REPONSE DYNAMIQUE

DES STRUCTURES ELANCEES A LA TURBULENCE DU VENT
Jury :
WIELGOSZ Christian Professeur, Université de Nantes Président

BORRI Claudio Professeur, Université de Florence Rapporteur
GIBERT René Professeur associé, IUP d’Evry Rapporteur
BIETRY Jacques Ingénieur Invité
CREMONA Christian Ingénieur Invité
VIRLOGEUX Michel Ingénieur consultant Invité
GRILLAUD Gérard Docteur Examinateur
JOUVE Pierre Professeur, Université de Nantes Examinateur
OLIVARI Dominico Professeur, Institut Von Karman Examinateur
Directeur de thèse : JOUVE Pierre

Laboratoire de Génie Civil N°ED 0367-009
Sommaire
Sommaire 2
Notations 6
Introduction 13
Chapitre 1 Modélisation du vent 17
1.1 La couche limite atmosphérique 18
1.2 Caractéristiques du modèle statistique 18
1.3 Le vent moyen 20
1.4 Le vent turbulent 21

1.4.1 Ecarts types et intensité de turbulence 21
1.4.2 Coefficients de corrélation et échelles de turbulence 22
1.4.3 Densités (inter)spectrales de puissance 23
a) Densités spectrales de puissance (cas où m = n et P = P’) 23
b) Densités interspectrales et fonctions de cohérences 25
1.5 Conclusion 27
2
Chapitre 2 Les forces aérodynamiques 28
2.1 Forces aérodynamiques s’exerçant sur un corps mal profilé 30

2.1.1 Forces s’exerçant sur une section fixe plongée dans un flux laminaire 30
2.1.2 Forces s'exerçant sur une section fixe en flux turbulent: approche quasi-statique 34
2.1.3 Section mobile, flux incident variable : approche quasi-statique 42
2.2 Les forces aéroélastiques 43

2.2.1 Formulation mixte temps-fréquence : Scanlan et Tomko (1971) 43
2.2.2 Généralisation à un mouvement quelconque de la section 45
a) Formulation fréquentielle 45
b) Formulation temporelle des forces aéroélastiques 47
2.3 Le détachement tourbillonnaire 52
2.4 Conclusion 53
Chapitre 3 Mesure des forces aérodynamiques 54
3.1 Mesure des coefficients stationnaires 55
3.2 Mesure de l'admittance aérodynamique 60
3.3 Mesure des coefficients aéroélastiques par la technique des oscillations libres 62
3.3.1 Bases théoriques 63
3.3.2 Application pratique 66
3.4 Méthode des oscillations forcées 75

3.4.1 Instrumentation de la maquette 75
3.4.2 Cas des oscillations pures 77
3.4.3 Mesure des rapports d’amplitudes, fréquences d’oscillations et déphasages 79
3.4.4 Cas du mouvement couplé 80
3.4.5 Résultats 81
3.5 Conclusion 87
3
Chapitre 4 Calcul de la réponse dynamiques de l'ouvrage 89
4.1 Discrétisation en éléments finis 90

4.1.1 Généralités 90
4.1.2 Mise en équations 91
a) Cas d’un élément 91
b) Assemblage 92
4.2 Calcul des modes propres de la structure 93
4.3 La méthode spectrale 95

4.3.1 Principe 95
4.3.2 Définition des repères 96
4.3.3 Calcul de la matrice de densité spectrale des forces turbulentes 97
4.3.4 Calcul de la matrice de densité spectrale des forces aéroélastiques 100
4.3.5 Mise en équations 102
4.3.6 Analyse de la stabilité 103
4.3.7 Résolution complète 105
4.3.8 Résolution simplifiée 106
4.4 Simulation numérique d’un signal de vent tridimensionnel spatialement corrélé 110
4.4.1 Méthode de Shinozuka 111
4.4.2 Méthode de Jin, Lutes et Sarkani (1997) 113
a) Hypothèses et notations 113
b) Evaluation de la moyenne conditionnelle dans le cas général 114
c) Evaluation de la variable aléatoire Y 116
d) Algorithme de résolution 116
4.4.3 Comparaison entre les deux méthodes 119
4.5 La méthode temporelle 121

4.5.1 Mise en équation 121
a) Formulation des forces turbulentes 121
b) Formulation des forces aéroélastiques 122
c) Formulation simplifiée des forces aéroélastiques 128
d) Mise en équation 128
4.5.2 Résolution du système dans la base modale 129
4.6 Introduction de simulations temporelles dans le calcul spectral 130

4.6.1 Méthode de « linéarisation équivalente » des coefficients aérodynamiques 130
Combinaison des charges modales extrêmes en cas de charges statiques équivalents 133
4.7 Conclusion 137
4
Chapitre 5 Validation expérimentale 138
5.1 Présentation des maquettes 139
5.2 Caractérisation mécanique des maquettes 142

5.2.1 Mesure des fréquences, amortissements et masses généralisées 142
5.2.2 Mesure des déformées modales 142
5.2.3 Calcul des déformées 144
5.2.4 Récapitulatif données expérimentales utilisés pour le calcul numérique 144
5.3 Caractérisation du vent de la soufflerie 148

5.3.1 Dispositif expérimental 148
5.3.2 Résultats 149
a) Profils de vitesse moyenne et de turbulence 149
b) Echelles de turbulence et coefficients de cohérence 152
c) Récapitulatif des résultats 154
5.4 Caractérisation aérodynamique des tabliers 156
5.5 Mesure des déplacements 158
5.6 Résultats 160

5.6.1 Considérations sur les incertitudes de calcul 160
5.6.2 Calcul spectral 161
a) Cas du tablier trapézoïdal 162
b) Cas du tablier rectangulaire 162
c) Conclusion 163
5.6.3 Calcul temporel 167
5.7 Conclusion 174
Conclusion 175
Perspectives 177
Bibliographie 179
5
Notations
B Dimension caractéristique de tablier : en général sa largeur
Br Variable aléatoire, module de Ŷrn

c Coefficient d’amortissement
C1Lh , C 2 Lh , C 3 Lh , C 4 Lh , d 3 Lh , d 4 Lh Coefficients des fonctions indicielles pour les forces de
portance et les déplacements verticaux
C1Mh , C 2 Mh , C 3 Mh , C 4 Mh , d 3 Mh , d 4 Mh Coefficients des fonctions indicielles pour les forces de
moment et les déplacements verticaux
C1M , C 2 M , C 3 M , C 4 M , d 3 M , d 4 M Coefficients des fonctions indicielles pour les forces de
moment et les déplacements angulaires
C1L , C 2 L , C 3L , C 4 L , d 3 L , d 4 L Coefficients des fonctions indicielles pour les forces de
portance et les déplacements angulaires
C xu , C xv , C xw , C yu , C yv , C yw , C zu , C zv , C zw Coefficients de cohérence
CD, CL, CM Coefficients de traînée, portance, moment dans un repère lié
au vent
CX, CY, CM Coefficients de traînée, portance, moment dans un repère
local
D Epaisseur de tablier
FD, FL, M Forces de traînée, portance, moment dans un repère lié au
vent
FX , FY , M Z Forces de traînée, portance, moment dans un repère local
FX , FY , M Z Forces de traînée, portance, moment moyennes dans un

repère local
FXt , FYt , FZt , M Xt , M tY , M tZ Forces et moments liés à la turbulence du vent dans un

repère local
FXa , FYa , FZa , M aX , M aY , M aZ Forces et moments aéroélastiques dans un repère local
FL , FLh , FM , FMh Fonctions indicielles dans le domaine spectral
f L , f Lh , f M , f Mh Fonctions indicielles dans le domaine temporel

f Fréquence
fc Fréquence de coupure
6
fN Fréquence de Nyquist
g L , g Lh Fonctions indicielles de Lin
g nr , ĝ nr Moyenne conditionnelle des termes u nr , û nr

h Déplacement vertical
H 1* , H *2 , H *3 , H *4 , A1* , A *2 , A *3 , A *4 Coefficients aéroélastiques adimensionnels de Scanlan
H 1 , H 2 , H 3 , H 4 , A1 , A 2 , A 3 , A 4 Coefficients aéroélastiques de Scanlan

~ ~
hm ,n Déplacements vertical et angulaire sinusoïdaux de
pulsations m et n
ĥ, ̂ Transformées de Fourier des déplacements vertical et

angulaire3
Iu , Iv , Iw Intensités de turbulence du vent
J Jacobien
K Fréquence réduite (en fait pulsation réduite)
K 1, K 2 Fréquences réduites des déplacements verticaux et
angulaires
l eae , m eae Forces de portance et de moment aéroélastiques agissant sur
un élément
L Longueur d’un élément poutre
L ae , M ae Forces de portance et moment aéroélastiques
~ ~
L ae mn , M ae mn Forces de portance et moment aéroélastiques générés par
des mouvements de flexion et de torsion sinusoïdaux
~ ~ ~ ~
L ae m , M ae m , L ae n , M ae n Forces de portance et moment aéroélastiques dans générés
par des mouvements de flexion ou de torsion sinusoïdaux
Lh , M h , L , M Forces de portance et moment aéroélastiques dans le
domaine temporel
L h 3 , M h 3 , L h 4 , M h 4 , L 3 , M 3 , L  4 , M  4 Composantes des forces de portance et moment
aéroélastiques dans le domaine temporel
Lxu , Lxv , Lxw , Lyu , Lyv , Lyw , Lzu , Lzv , Lzw Echelles de turbulence
mg Masse généralisée
ndlt Nombre de degrés de liberté total de la structure
ne Nombre d’éléments poutres de la structure
P, P' Points de l’espace
7
R mp
jk Fonction de covariance entre la composante j du signal m et
la composante k du signal p
Re Nombre de Reynolds
Su ,Sv ,Sw Spectres des composantes u, v, w de la turbulence du vent
S mp
jk Densité interspectrale de puissance entre la composante j du
signal m et la composante k du signal p

St Nombre de Strouhal
T Période
t Temps
r
U Vecteur vitesse instantané du vent
U Norme du vecteur vitesse moyenne du vent
U ref Vitesse moyenne de référence
U 10 m Vitesse moyenne de référence à 10m

U Vitesse de vent instantanée
~
U Vitesse de vent instantanée dans un plan perpendiculaire à
un élément poutre
u, v , w Composantes de la turbulence du vent dans un repère lié au
vent puis dans un repère lié à un élément poutre
u L , vL , wL Composantes de la turbulence du vent dans un repère local
u nr , û nr rième composante du signal tridimensionnel généré au point n

dans les domaines temporel et spectral
nm nm
Wrjnm , Ŵrjnm Composante des matrices W et Ŵ
r r r
x , y, z Base de projection liée à un élément poutre
x1, y1, z1, rx1, ry1, rz1, x2, y2, z2, rx2, ry2, rz2 12 degrés de liberté d’un élément poutre
yP,zP Positions suivant les directions y et z du point P
Yrn , Ŷrn Variable gaussienne centrée, composantes de u nr , û nr
z0 Rugosité d’un site
z ref Altitude de référence
 Déplacement angulaire
u ,v ,w Fonctions de cohérence de Davenport
 Xu ,  Yu ,  Mu ,  Xv ,  Yv ,  Mv ,  Xw ,  Yw ,  Mw « Racines carrées » des fonctions d’admittance

aérodynamique
8
 Incidence du vecteur vent dans un repère lié à un élément poutre
 ki Phase associée au point k et à la fréquence if
 Amortissement rapporté au critique

 Densité de l’air
u , v , w Ecarts types des trois composantes de la turbulence du vent
 q ,  dyn
q , q
qs
Ecarts types total, dynamique et quasi statique des déplacements
modaux
 Temps utilisé dans les produits de convolutions
 Coordonnée paramétrique des éléments poutres
Xu , Xv , Xw Fonctions d’admittance aérodynamique de traînée
Yu , Yv , Yw Fonctions d’admittance aérodynamique de portance
Mu , Mv , Mw Fonctions d’admittance aérodynamique de moment
1 , 2 Pulsations des mouvements de flexion et torsion

 Pulsation
dP Vecteur des déplacements en un point d’un élément poutre dans le
repère local
d Vecteur des déplacements de la structure dans la base nodale et le
repère général
de Vecteur des déplacements dans un repère local
dn Vecteur des déplacements d’un élément poutre dans la base nodale et

le repère local
d ng Vecteur des déplacements d’un élément poutre dans la base nodale et
le repère général
E Vecteur contenant les modules d’Young des éléments poutres
fV , fS Vecteurs forces extérieures volumique et surfacique
f te Vecteur des buffeting forces dans un repère local
f aee Vecteur des forces aéroélastiques dans un repère local
f tne Vecteur des buffeting forces dans la base nodale et le repère local
f aene Vecteur des forces aéroélastiques dans la base nodale et le repère
local
f tng Vecteur des buffeting forces dans la base nodale et le repère général
9
Ftn , Ftm Vecteurs des buffeting forces globaux dans les bases nodale et modale
Fae Vecteur des forces aéroélastiques
Ft Vecteur des buffeting forces
F Vecteur des forces extérieures
G Vecteur contenant les modules de cisaillement des éléments poutres
IX , IY , IZ Vecteurs contenants les inerties surfaciques des éléments poutres
m ddl Vecteur contenant les masses et inerties ponctuelles affectées à
chaque degré de liberté de la structure

MX ,MY ,MZ Vecteurs contenant les moments de flexion et torsion des éléments
poutres
M Xi , M Yi , M Zi Vecteurs contenant les moments de flexion et torsion modaux des
éléments poutres
N Vecteur contenant les efforts normaux des éléments poutres
Ni Vecteur contenant les efforts normaux modaux des éléments poutres
q Vecteur des déplacements de la structure dans la base modale
S Vecteur contenant les sections des éléments poutres
TX , TY Vecteurs contenant les efforts tranchants des éléments poutres
TXi , TYi Vecteurs contenant les efforts tranchants modaux des éléments
poutres
un Vecteur vent tridimensionnel synthétique
x Vecteur position d’un point d’un élément poutre
xn Vecteur des coordonnées d’un nœud
Zj jième vecteur propre
Zj' Dérivée du jième vecteur propre par rapport à l’abscisse curviligne
 Vecteur des déformations élastiques
 Vecteur polynomial de l’incidence du vent
 Vecteur des forces aéroélastiques dans le domaine temporel
 eNc ,  eTcX ,  eTcY ,  eMc X ,  eMc Y ,  eMc Z Vecteurs des cas de charges extrêmes
A Matrice aéroélastique globale
10
Ae Matrice aéroélastique élémentaire exprimée dans un repère local
A ng Matrice aéroélastique élémentaire exprimée dans un repère général
B Matrice aéroélastique globale
Be Matrice aéroélastique élémentaire exprimée dans un repère local
B ng Matrice aéroélastique élémentaire exprimée dans un repère général
C ad Matrice des coefficients aérodynamiques
C12345 Matrice polynomiale des coefficients aérodynamiques
C Matrice des coefficients aérodynamique avec prise en compte de
l’admittance aérodynamique
C 1 , C 2  Matrice des coefficients aérodynamique avec prise en compte de
l’admittance aérodynamique
C aeng , K ng
ae Matrices amortissement et raideur aérodynamiques élémentaires
exprimées dans la base nodale et le repère général

C ae , K ae Matrices amortissement et raideur aérodynamiques globales
exprimées dans la base nodale et le repère général
Ĉ ae , K̂ ae Matrices amortissement et raideur aérodynamiques globales
exprimées dans la base modale

C eae , K eae Matrices amortissement et raideur aérodynamiques élémentaires
exprimées dans la base nodale et le repère local

D Matrice aéroélastique globale
De Matrice aéroélastique élémentaire exprimée dans un repère local
D ng Matrice aéroélastique élémentaire exprimée dans un repère général
E Matrice de rigidité
H Matrice de transfert
m , c, k Matrices masse, amortissement et raideur élémentaires
m j,cj,k j Matrices masse, amortissement et raideur élémentaires associées à
l’élément j
m j ', c j ', k j ' Matrices masse, amortissement et raideur élémentaires associées à
11
l’élément j insérées dans les matrices carrées nulles de dimension ndlt
M, C, K Matrices masse, amortissement et raideur globales exprimées dans la
base nodale
M̂, Ĉ, K̂ Matrices masse, amortissement et raideur globales exprimées dans la
base modale
N Matrice d’interpolation
P12 Matrice de changement de repère pour un vecteur 12 composantes
PV Matrice de changement de repère vent
PL Matrice de changement de repère local-général
S eV, e ' Matrice de densité interspectrale élémentaire de la turbulence du vent
S ng n 'g Matrice de densité interspectrale élémentaire des buffeting forces

f t ,f t
SF Matrice de densité interspectrale des forces aérodynamiques
X Matrice reliant  à u n
nm nm m m
W , Ŵ Matrices de pondération de u et û
 Matrice des vecteurs propres
12
Introduction
L’amélioration constante des technologies dans le domaine du génie civil et des matériaux permet de
concevoir des ouvrages sans cesse innovants, de plus en plus légers et élancés. Les ponts suspendus ou à
haubans en sont les représentants extrêmes. Il résulte de cette évolution une sensibilité accrue aux
sollicitations fluctuantes du vent. Les ouvrages peuvent alors adopter un comportement vibratoire très
marqué susceptible de les mener à la ruine ou tout simplement d'en rendre l'utilisation périlleuse.
Les bureaux d’études chargés du dimensionnement des ouvrages ont besoin d'introduire dans leurs
calculs les sollicitations induites par les tempêtes les plus sévères. Dans le cas de bâtiments massifs de
formes suffisamment simples, la prise en compte des effets du vent est estimée à partir des règles en
vigueur (Neige et Vent ou Eurocode) basées sur la détermination de charges statiques forfaitaires.
Dans le cas de formes complexes à comportement dynamique marqué, des études plus poussées et
personnalisées sont nécessaires, d’une part pour estimer les charges aérodynamiques s’exerçant sur
l’enveloppe, et d’autre part pour déterminer la réponse dynamique de la structure sous les sollicitations
de la turbulence du vent. Différentes approches peuvent être envisagées aujourd’hui :
 En ce qui concerne la détermination des efforts s’exerçant sur l’enveloppe d’une structure, il existe
deux types de démarches : la première consiste à mesurer en soufflerie les efforts générés par le vent
sur des maquettes rigides de section constante appelées "maquettes sectionnelles" ; la seconde
consiste à calculer numériquement ces mêmes efforts à l’aide de codes C.F.D. (Computational Fluid
Dynamics) basés sur la résolution discrète des équations de Navier-Stokes. Si cette seconde approche
est en plein essor et s’annonce prometteuse, la première reste aujourd’hui, et pour, sans doute,
quelques années encore, la plus opérationnelle.
13
 Pour estimer l’amplitude des déplacements, étape obligatoire avant le calcul des contraintes, deux
types de méthodes sont actuellement utilisées : la plus pragmatique consiste à mesurer directement
ces déplacements en soufflerie sur des maquettes dites "aéroélastiques" reproduisant, à échelle
réduite, les comportements aérodynamique et mécanique de la structure étudiée ; la seconde, plus
théorique, consiste, dans un premier temps, à mesurer les caractéristiques aérodynamiques de la
structure à l’aide de maquettes sectionnelles, pour estimer ensuite, à l’aide d’un code de calcul des
structures, la réponse de l'ouvrage complet à des conditions de vent données. C’est cette dernière
approche que nous avons développée et exposons dans le présent mémoire.
Il s’agit tout d’abord de s’intéresser à la modélisation du vent naturel turbulent dans la partie basse de
l'atmosphère: la couche limite atmosphérique. Pour définir ce processus tridimensionnel fortement
aléatoire, nous utilisons un certain nombre d’outils statistiques permettant d'intégrer les corrélations
adéquates dans le temps et dans l’espace. Nous sommes alors amenés à substituer la dimension
fréquentielle ou spectrale à la dimension temporelle. C’est l’objet de notre premier chapitre.
Dans un second chapitre, nous présentons des modèles qui permettent, à partir de la bonne connaissance
du vent au voisinage de la structure et de la caractérisation aérodynamique des profils transversaux de
l’enveloppe de l’ouvrage étudié, d’obtenir une estimation des forces s’exerçant sur la structure. Nous
verrons que des forces de différentes natures sollicitent dynamiquement les ouvrages :
 les forces fluctuantes directement liées au caractère turbulent du vent souvent appelées buffeting
forces; on y associe la notion d'admittance aérodynamique représentant l'impact différencié des
structures tourbillonnaires suivant leur taille;
 les forces dites aéroélastiques engendrées par les mouvements de la structure pouvant être
responsables d'un comportement instable de la structure;
 les forces liées à un phénomène aérodynamique propre aux corps mal profilés et élancés : le
détachement tourbillonnaire.
Nous consacrons le troisième chapitre aux différentes techniques expérimentales permettant de

caractériser les propriétés aérodynamiques des profils transversaux des structures élancées. Les méthodes
de mesure des buffeting forces et de l'admittance aérodynamique sont brièvement présentées. On
insistera plus particulièrement sur deux techniques expérimentales permettant de caractériser les forces
aéroélastiques : la technique des oscillations libres et celle des oscillations forcées. La dernière a été
l’objet d’une publication à l’occasion d’un congrès du Wind Engineering (de la Foye, 2000). Nous
présentons les résultats obtenus pour deux profils de tabliers. Ces résultats seront utilisés par la suite lors
de notre procédure de validation.
17
Le quatrième chapitre constitue la pierre angulaire de nos travaux. Nous y montrons comment, à partir
de la connaissance des forces aérodynamiques s’exerçant sur la structure, nous résolvons numériquement
le problème mécanique posé. Deux méthodes, s’appuyant sur une discrétisation en éléments finis de la
structure, sont développées et présentées :
 la méthode spectrale, de nature statistique, est basée sur la résolution des équations de la mécanique
dans le domaine fréquentiel;
 la méthode temporelle, semi-déterministe, est basée sur la génération d’un signal de vent synthétique
spatialement corrélé, et la résolution des mêmes équations, cette fois dans le domaine temporel.
Enfin, dans le dernier chapitre, nous présentons la procédure mise en place pour valider les codes
développés. Pour cela, nous plaçons en soufflerie deux maquettes aéroélastiques, représentant des travées
centrales de tabliers de ponts, et en mesurons la réponse dynamique pour différentes configurations de
vent. Des calculs sont alors effectués, avec, en conditions d’entrée, les caractéristiques du vent mesurées
en soufflerie, les caractéristiques mécaniques de la structure, et les propriétés aérodynamiques de la
structure présentées au chapitre 3. Nous procédons alors à une analyse comparative des différentes
méthodes de calcul numérique avec les résultats obtenus expérimentalement, et tentons de mettre en
évidence, d'une part l’influence des différents paramètres d’entrée sur les résultats obtenus
numériquement, puis, dans un deuxième temps, d'autre part les forces et limites d’application des
différentes approches.
Précisons que notre travail s’inscrit dans la continuité de nombreux travaux réalisés sur le sujet depuis
une quarantaine d'années. En effet, des approches spectrales ont été proposées dès 1962 par A. Davenport
pour calculer la réponse des structures élancées à la turbulence du vent. L'application au calcul spectral
du concept de forces aéroélastiques, issu de l'aéronautique, fut initiée en 1971 par R.H. Scanlan.
M. Attou développa au CSTB, dans le cadre d'une thèse de doctorat soutenue en 1994, un premier
module de calcul spectral en langage Fortran. La géométrie des structures était limitée à deux
dimensions, la vitesse moyenne du vent horizontale et les forces aéroélastiques négligées. Un module
temporel fut également élaboré dont les limites étaient similaires à celles du module spectral. Le signal
de vent, bidimensionnel, était généré par un programme basé sur l’algorithme de Shinozuka.
De 1994 à 1997, P. Bourcier, successeur de M. Attou, s'est essentiellement attaqué au problème de la

modélisation des effets aéroélastiques et de leur introduction dans des codes de calcul spectral et
temporel. Il n'a cependant pas finalisé ses travaux et laissé derrière lui des programmes difficilement
18
exploitables. Parallèlement, A. Patron-Solares développa au LCPC un ensemble de modules de calculs
spectraux et temporels intégrant les modèles de forces aéroélastiques proposés par Scanlan.
Notre tâche a d'abord consisté à reprendre en amont le travail entamé par P. Bourcier. Nous avons ensuite
fait le choix de redévelopper entièrement les outils spectraux et temporels sous l'environnement Matlab.
En outre, les éléments suivants ont été introduits:
 Une définition multizonale des caractéristiques du vent turbulent;

 Un générateur de processus temporel tridimensionel plus performant que le précédent basé sur les
travaux de Jin et al. (1997);
 Une méthode d'optimisation de la linéarisation des forces aérodynamiques dans l'optique de calculs
spectraux;
 La prise en compte de l'admittance aérodynamique dans les modules spectraux et temporels;
 La réelle prise en compte de l’incidence du vent et de la déformée statique ;
 Une méthode de combinaison des charges modales extrêmes et calcul de charges statiques
équivalentes à l'intention des bureaux d'études.
Notre travail est donc essentiellement axé sur le calcul numérique des déplacements des structures
élancées soumises à la turbulence du vent. Nous l’avons cependant inscrit dans le contexte plus général
de la prise en compte des effets du vent dans le dimensionnement des grands ouvrages, depuis l’étude
climatologique du site d’implantation jusqu’à l’estimation des contraintes extrêmes s’exerçant au sein de
la structure. Interviendront donc tout au long du mémoire les disciplines suivantes : climatologie,
aérodynamique, mécanique des vibrations et résistance des matériaux.
19
Chapitre 1 Modélisation
du vent
Tous les ouvrages du génie civil sont entièrement situés dans la plus basse couche de l'atmosphère
appelée couche limite atmosphérique. A cette altitude, le vent, très fluctuant dans le temps et dans
l'espace, est couramment assimilé à un processus aléatoire tridimensionnel. On est donc naturellement
amené à utiliser des approches statistiques pour modéliser ses fluctuations.
Ainsi, l'objet de ce premier chapitre est de présenter une série d'outils statistiques permettant de
caractériser le vent turbulent au voisinage des structures de génie civil. Nous verrons qu'il est alors
commode de substituer à la dimension temporelle, la dimension fréquentielle ou "spectrale" . Cette
notion essentielle reviendra ensuite comme un leit motiv dans notre travail.
17
Chapitre 1 Modélisation du vent
1.1 La couche limite atmosphérique
Pour une altitude supérieure à quelques centaines de mètres, la vitesse du vent résulte de l’interaction
entre les forces de Coriolis engendrées par la rotation de la terre et les forces créées par le champ de
pression atmosphérique. La turbulence est presque nulle. Pour une altitude inférieure à 300 mètres
environ, là où se trouvent les ouvrages d'art, les forces de frottement dues à la rugosité du sol et les
phénomènes thermiques deviennent prépondérants vis-à-vis des forces de Coriolis ; ils engendrent des
fluctuations de la vitesse du vent, dans le temps et dans l’espace, susceptibles d’exciter les structures
souples. Cette zone est couramment nommée couche limite atmosphérique.
En réalité, seul le champ de vent à proximité directe de l’ouvrage étudié intéresse notre travail.
Cependant, même ainsi réduite, la zone d’étude reste difficile à caractériser. On comprend bien qu’une
approche expérimentale ne permet d’accéder simultanément qu’à un nombre très limité de points de
l’espace. Les codes de calcul, résolvant les équations de la mécanique des fluides dans un espace et un
temps discrétisés, permettent bien, en théorie, d’accéder aux caractéristiques du vent en tout point du
temps et de l’espace ; seulement, la complexité géométrique de la structure étudiée, additionnée à celle
d’un environnement qu’on ne peut pas ignorer, rend, encore aujourd’hui, les temps de calcul
rédhibitoires. Des modèles mathématiques, essentiellement statistiques, ont donc été développés,
permettant de caractériser le vent par un nombre limité de paramètres, accessibles par la mesure, et jugés
suffisamment pertinents pour permettre d’en déduire la réaction d’un ouvrage soumis à son influence.
1.2 Caractéristiques du modèle statistique
L’analyse spectrale de la vitesse du vent dans la couche limite turbulente, caractérisant la distribution en
fréquence de son énergie, permet de mettre en évidence plusieurs échelles de fluctuation.
S
4 jours
1 an 1 min.
12 h.
1 h. 10 min.
1/T
Figure 1.1 - Spectre de Van der Hoven
18
La figure 1.1 montre l’allure d’un spectre de densité de puissance représentatif de la vitesse horizontale
du vent à 100 mètres au dessus du sol d’après Van der Hoven, 1957.
On constate qu’une partie importante de l’énergie est distribuée autour d'une période de l’ordre de la
minute, correspondant à l’échelle des mouvements turbulents. On observe également un trou d’énergie
pour des périodes comprises entre 10 minutes et une heure. Par conséquent, sur une période incluse dans
ce trou d’énergie, la vitesse du vent est localement assimilable à un processus stationnaire. On peut donc
aussi considérer la vitesse du vent comme la somme d’une composante moyenne et d’une composante
fluctuante. La figure 1.2 représente un enregistrement de la vitesse du vent sur 5 minutes.
Figure 1.2 – Signal de la composante longitudinale du vent sur cinq minutes
Pour des raisons de commodité, le vecteur vitesse du vent est modélisé par un processus supposé
r r r r
gaussien, stationnaire1 et ergodique2. Ainsi, par rapport à une base x , y, z  où x est parallèle au flux
r
moyen et z est généralement vertical, la vitesse du vent en tout point du site peut s'écrire :
r r r r
U  (U  u) x  v y  w z (1.1)
r
w
r r
U v
r r
z u
r
r
y r U
x
Figure 1.3 - Décomposition du vecteur vitesse du vent
1
Un processus aléatoire est dit stationnaire si sa structure de probabilité est invariante vis-à-vis d’un changement de
l’origine des temps.
2
L’hypothèse d’ergodicité permet de remplacer des moyennes d’ensemble par des moyennes temporelles sur une
réalisation unique du processus.
19
 U est la vitesse moyenne au point considéré ; elle caractérise la vitesse à laquelle sont convectés les
tourbillons;
 u, v, w sont trois processus gaussiens stationnaires centrés représentant les fluctuations du signal
autour de la moyenne ; ils caractérisent les variations de la vitesse du vent dues aux passages des
tourbillons.
1.3 Le vent moyen
En aval d’un site relativement homogène (étendue de mer, bois, ville), on peut considérer que la direction
moyenne du vent est identique en tout point et que son module ne dépend que de la hauteur. Différentes
lois de variation peuvent alors être envisagées :
La loi logarithmique, fondée sur la théorie de la couche limite donne :
  z  z 
U ( z )   ln   ln  ref   U (1.2)
  ref
  z0   z0 
où z est la hauteur du point considéré au dessus du sol, z0 le paramètre de rugosité exprimé en mètres,
U ref une vitesse de référence à la hauteur de référence zref.
A titre d’exemple, le tableau ci-dessous indique les ordres de grandeur proposés par Bietry et al. (1978) :
Site Classe z0 (m)

Océans et grands lacs I 0,001 - 0,01
Aéroport, rase campagne II 0,01 - 0,1
Bocage, banlieue III 0,1 - 0,5
Zones urbaines IV 0,5 - 1,5
Très grande ville V 1,5 - 2,5
Tableau 1.1 - Coefficients z0 caractérisant la rugosité du site
Lorsque la couche limite atmosphérique n’est pas bien développée, une loi puissance empirique peut être
utilisée :

 z 
U(z)  U ref   (1.3)
 z ref 
20
où Uref est la vitesse à la hauteur z ref , et  un exposant sans dimension dépendant de la rugosité du site.
Le tableau ci-dessous donne une idée de l’équivalence -z0 :
z0 (m) 0,005 0,07 0,3 1,0 2,5

 0,10 0,15 0,20 0,25 0,35
Tableau 1.2 – Coefficients 
1.4 Le vent turbulent
On note (u, v, w) le vecteur représentant les composantes fluctuantes de la vitesse du vent ; u représente
la composante longitudinale parallèle au flux moyen, w la composante verticale et v est la composante
latérale.
Il s’agit alors d’en décrire les variations temporelles et spatiales ; ces dernières jouent un rôle essentiel
dans l’étude des structures élancées dont la taille peut être largement supérieure aux dimensions
caractéristiques des tourbillons. Chacune des composantes u, v, w, en un point de l’espace donné, sera
modélisée par un processus gaussien centré, stationnaire et ergodique. Nous pourrons ainsi appliquer un
certain nombre d’outils statistiques permettant de les caractériser dans l’espace et dans le temps .
1.4.1 Ecarts types et intensité de turbulence
On caractérise tout d'abord l'énergie cinétique fluctuante massique moyenne du flux par l'écart type de
chaque composante turbulente:
T
1
 u m P  
u
2
m ( P, t ) dt um = u, v, w (1.4)
T 0
où T est la durée d’observation du signal. Pour T de l’ordre de 10 minutes, on a vu que le vent pouvait
être assimilé à un signal stationnaire. On constate expérimentalement que ces écarts types sont
proportionnels à la vitesse moyenne U . On remarque également que, pour des sites suffisamment
homogènes, ils sont indépendants de la position spatiale et, en particulier, ne dépendent pas de l’altitude
z. Par contre, ils augmentent avec la rugosité du sol.
21
On utilise plus souvent les notions alternatives d’intensités de turbulence. Adimensionelles, elles sont
définies par:
u  
Iu  , Iv  v , Iw  w (1.5)
U U U
1.4.2 Coefficients de corrélation et échelles de turbulence
D’une manière générale, afin de traduire la concomitance entre deux composantes de la vitesse du vent
en deux points P et P’ et à deux instants t et t+, on introduit la notion statistique d’inter-corrélation
définie par :
T
1
 u m ,u n P, P' ,   
T u
0
m ( P, t ) u n ( P' , t  ) dt um, un = u, v, w (1.6)
vérifiant notamment la double inégalité   u m  u n   u m ,u n   u m  u n
Dans la plupart des cas, on se limite au cas où m = n ; on parle alors de corrélations spatiales ( = 0),
temporelles (P = P’), ou spatio-temporelles (  0, P  P’) pour une composante de vent donnée. On
définit, à partir de ces dernières, les notions d’échelles de turbulence. Celles-ci ont une double
signification:
 physique : les dimensions des tourbillons les plus énergétiques dans chaque direction ;
 mathématique : la distance, dans chaque direction, à partir de laquelle la corrélation entre deux
signaux devient négligeable.
On écrit ainsi:

U
Lxu m P  
um
2 0
u m ,u m (P, P, ) d um = u, v, w (1.7)

1
Lyu m P  
um
2 0
u m ,u m ( y P , y P  y, 0) dy um = u, v, w (1.8)

1
Lzu m P  
um
2 0
u m ,u m (z P , z P  z, 0) dz um = u, v, w (1.9)
22
On peut se demander pourquoi la définition de Lxu m diffère de celles des Lyu m , Lzu m . La raison est
double : d'abord il est délicat d'un point de vue expérimental de placer deux appareils de mesure l'un
derrière l'autre pour des problèmes évidents de sillage ; ensuite, les grandeurs Lyu m , Lzu m interviennent
effectivement par la suite pour caractériser des degrés de corrélation entre deux points de l'espace tandis
que la grandeur Lxu m intervient plutôt en tant qu'indicateur de corrélation temporelle en un point donné.
Inversement, la connaissance expérimentale des valeurs Lxu m , Lyu m , Lzu m , permet d’estimer la corrélation
spatiale par la formule empirique suivante:
 2 2 2 
  x   y   z  
 u m , u m (P, P' ,0)  exp         um = u, v, w (1.10)
  Lxu   Lyu   Lzu  
 m   m   m 
 
x  x P  x P ' y  y P  y P ' z  z P  z P '
1.4.3 Densités (inter)spectrales de puissance
Plutôt que de travailler directement avec les fonctions de corrélation définies en (1.6), on préfère
manipuler leurs transformées de Fourier, appelées densités (inter)spectrales de puissance. Ce sont elles
qui interviendront dans le calcul numérique de la réponse à la turbulence des structures élancées. Dans le
cas particulier où m = n et P = P’, on parle simplement de densité spectrale de puissance; celle-ci
représente la répartition fréquentielle de l’énergie turbulente en un point P donné et pour une composante
m donnée. Dans les cas plus généraux où m  n et/ou P  P’, les densités interspectrales de
puissance indiquent, pour chaque gamme de fréquence, avec quel degré de corrélation l’énergie cinétique
est fournie à deux signaux distincts.
a) Densités spectrales de puissance (cas où m = n et P = P’)
Dans ce cas précis, on définit la fonction densité spectrale de puissance par :

S u m (f )  2


 u m ,u m () e  2if d um = u, v, w (1.11)
Comme  u m , u m est une fonction paire, S u m est une fonction à valeurs réelles. Mais, comme toutes les
densités spectrales de puissance, elle vérifie surtout la propriété suivante :
23

S ( n ) df   u m
2
um um = u, v, w (1.12)
0
La forme typique du spectre de la composante u est représentée figure 1.4. On identifie trois régions
caractéristiques :
 une région de basses fréquences, associée aux grands tourbillons anisotropes (dont le comportement
est différent dans les trois directions de l'espace) qui dépendent fortement des conditions extérieures,
 une région de moyennes fréquences, concentrant l’essentiel de l’énergie, associée aux structures
tourbillonnaires issues de la désagrégation des gros tourbillons par les forces de cisaillement,
 une région de hautes fréquences, associée à une turbulence isotrope, dont l’énergie est convertie en
chaleur par dissipation visqueuse, appelée domaine d’équilibre universel.
Production Dissipation
d’énergie d’énergie
Transfert
d’énergie
f
Gros Tourbillons Domaine d’équilibre
tourbillons énergétiques universel
Figure 1.4 - Spectre de la composante longitudinale de la vitesse du vent
Différentes formulations sont proposées dans la littérature pour approcher cette distribution S u m .
Les plus couramment utilisées sont, pour la composante u, celles de :
4 Lux / U
Von Karman (1948) : S u (f )   2u (1.13)
1  70,7 fL u
x / U
2
 5 6
17 z / U
Kaimal (1972) : S u (f )   2u (1.14)
1  33 fz U
53
24
32z / U
Simiu (1974) : S u (f )   2u (1.15)
1  50 fz U
5 3
0,65.1200 / U 10 m
Davenport (1961) : S u (f )   2u (1.16)
1  1200f U 10 m 
2
 5 3
Pour v et w, les formulations proposées le plus généralement sont :
 2
  2

Lv 1  188,4 2f L x   Lw 1  188,4 2f L x  
v w
4 x 4 x
U   U   U   U   2
Von Karman: S v (f )    2v S w (f )   w
11 6 11 6
 2
  2

1  70,7 2f L x   1  70,7 2f x  
L
v w
  U     U  
 
(1.17)
4,2z / U 1,4z / U
Kaimal : S v (f )   2v S w (f )   2w (1.18)
1  9,5 fz U  53
1  5,3 fz U  53
2,5z / U 0,56z / U
Simiu : S v (f )   2v S w (f )   2w (1.19)
1  9,5 fz U  53
1  10 fz U  53
Si l'on s'amuse à intégrer ces fonctions sur 0   , on s'aperçoit que la relation (1.12) n'est pas vérifiée.
17 2
Par exemple, l'intégration analytique de (1.14) donne  u . Cela n'est pas bien grave dans la mesure où
22
on ne cherche jamais à modéliser l'intégralité du spectre mais uniquement la plage susceptible d'exciter
une structure souple.
b) Densités interspectrales et fonctions de cohérences
Entre deux points, et pour une composante m donnée, la densité interspectrale de puissance s’exprime :

S u m ( P, P ' , f )  2

 u m , u m (P, P' , )e  2 if d (1.20)
La partie imaginaire de cette fonction complexe, donnant des informations sur le déphasage entre les
deux signaux considérés, est généralement ignorée dans les calculs. On se contente donc de faire
25
intervenir son module, directement relié aux spectres des deux composantes par une fonction dite de
cohérence  u m :
2
S u m ( P, P ' , f )
 2
um ( P, P ' , f )  um = u, v, w (1.21)
S u m ( P, f ) S u m ( P ' , f )
Cette fonction est assez mal connue car les enregistrements de vent montrent trop de dispersion pour
bâtir un modèle précis. Néanmoins, la formulation généralement retenue est celle de Davenport (1968) :
fZ 2 x 2 y 2 z
 2u m (P, P' , f )  exp[ C ux m ( ) 2  C uy m ( ) 2  C zu m ( ) 2 ] um = u, v, w (1.22)
U Z Z Z
où Z est la hauteur moyenne entre les points P et P’. Les coefficients C ux m , C uy m , C zu m sont appelés
coefficients de cohérence. Adimensionnels, leurs valeurs habituelles sont de l’ordre de 10 à 30 pour la

composante u, de 5 à 15 pour v, et de 10 à 50 pour w.
En théorie, leur description nécessiterait une mesure par couple de points. Toutefois, pour des raisons
pratiques, ceci n'est jamais fait. On considère généralement un point fixe, l'autre point balayant l'espace.
Les valeurs retenues sont les moyennes des différentes valeurs des coefficients obtenues au cours de la
scrutation.
Pour un même point, et entre deux composantes m et n données, on peut également définir une fonction
de cohérence permettant de calculer la densité inter-spectrale de puissance. On remarque
expérimentalement que seules les composantes u et w présentent une fonction de cohérence non
négligeable. Celle-ci est définie par (Attou, 1994):
S u , w ( P, f )  0,3
 u ,w ( P, f )   (1.23)
S u ( P, f ) S w ( P, f )  0,247  0,344 log(z )  0,036 log(z )  2
1  0,4 (fLxu U ) 2
Entre deux points P et P’, on peut définir une fonction de cohérence par :
26
S u , w ( P, P ' , n )  Z x 2 y z 
 u ,w ( P, P' )   exp   C ux , w ( )  C uy ,w ( ) 2  C uz , w ( ) 2 
 (1.24)
S u ( P, n ) S w ( P ' , n )  2U Z Z Z 
avec C iu ,w  f (C iu , C iw , z, Lxu )
Ainsi, pour chaque couple de points de l’espace P et P’, on peut définir une matrice de densité inter-
spectrale de puissance s’écrivant :
 S u , u ' (S v , u ' ) S w , u ' 

 
S u (P, P' )  (S u , v ' ) S v , v ' (S w , v ' ) (1.25)
 S u , w ' (S v , w ' ) S w , w ' 
 
où tous les termes entre parenthèses sont supposés nuls faute de modèles appropriés.
1.5 Conclusion
Les modèles présentés dans ce chapitre apportent une réponse au problème de la modélisation de la
répartition fréquentielle et spatiale de l'énergie du vent. On verra dans les chapitres suivants que ce type
de modélisation est tout à fait adapté pour déterminer les effets du vent sur des structures élancées de
grandes dimensions.
Dans le cas des ouvrages exceptionnels pour lesquels l'effet du vent est dimensionnant (pont de Rion-
Antirion en Grèce ou Viaduc de Millau par exemple), les différents coefficients du modèle sont établis
par une étude spécifique pouvant inclure des mesures sur site, des essais sur maquette topographique en
soufflerie atmosphérique ou des simulations numériques.
Ces études peuvent être relativement coûteuses et donc, pour des projets de moindre envergure, ou moins
sensibles aux effets du vent, on peut se contenter d’estimer les différentes valeurs caractéristiques du
vent à partir d’études effectuées sur des sites similaires ou des données issues de la littérature.
27
Chapitre 2 Les forces
aérodynamiques
2 Les forces aérodynamiques
Après avoir défini un modèle permettant de qualifier le champ de vitesse du vent sur un site donné, nous
allons nous intéresser à la modélisation des efforts éoliens s’exerçant sur un ouvrage placé dans ce site.
Un tel ouvrage est soumis à un champ de pression dont l’intégration sur une surface donnée permet de
calculer des efforts résultants.
La première tâche affectée à l’outil numérique que nous allons développer sera de déduire des
caractéristiques spatio-temporelles du vent la répartition spatiale des efforts résultants et leur évolution
dans le temps.
Les structures élancées, dont font partie les ponts suspendus et à haubans, présentent l’avantage d’être
facilement modélisables via une discrétisation en éléments finis de type « poutre », de caractéristiques
mécaniques et aérodynamiques données ; les efforts résultants seront alors calculés par intégration
d’efforts linéiques sur ces poutres élémentaires.
Ces efforts linéiques sont préalablement mesurés en soufflerie sur de simples tronçons à échelle réduite,
de sections constantes, et de longueurs suffisamment importantes pour marginaliser les effets de bouts :
ce que l’on appellera par la suite des maquettes « sectionnelles ». Les sections des ouvrages vraie
grandeur étant généralement variables, autant de sections significatives que nécessaire devront être
testées, compte tenu de la précision de calcul attendue. Cette tâche pourra également être accomplie à
l’aide d’outils numériques de type CFD (Computational Fluid Dynamics) ; nous y reviendrons plus loin.
Dans le cas le plus trivial d’une structure fixe soumise à un flux incident laminaire et constant, on
modélise simplement les forces aérodynamiques, via l’introduction de coefficients adimensionnels,
caractéristiques de la forme de la section et de son incidence par rapport à la direction de l’écoulement.
Pour chaque incidence de vent, on établit ainsi une relation, généralement linéaire, entre la pression
dynamique amont et les forces aérodynamiques.
28
Chapitre 2 Les forces aérodynamiques
Dans le cas d’une structure réelle, l'écoulement incident est turbulent et génère donc des signaux de
forces fluctuants. On verra que l’utilisation des coefficients évoqués précédemment, bien que tentante,
est alors moins appropriée ; on introduit donc une nouvelle notion permettant de modéliser la
décroissance de corrélation spatiale du vent autour d'une tranche de structure à mesure que les
fluctuations deviennent plus rapides : l’admittance aérodynamique.
D'autre part, en entrant en vibration sous l'effet de la variation temporelle des forces extérieures,
l'ouvrage modifie la structure du flux qui l'environne et par conséquent les forces aérodynamiques. Ceci
sera modélisé par l’introduction de forces supplémentaires dites d'auto-excitation ou aéroélastiques selon
le modèle proposé par Scanlan.
Enfin, à l'instar de la plupart des corps mal profilés, les tabliers de ponts peuvent être soumis, lorsque le
flux incident n'est pas trop turbulent, à des sollicitations périodiques dues à un phénomène appelé
détachement tourbillonnaire. Nous présenterons quelques modèles permettant de modéliser ces
sollicitations. Cependant, ceux-ci ne seront pas intégrés dans les codes de calcul développés par la suite.
L'ensemble des modélisations présentées dans ce chapitre seront faites dans l'optique de leur utilisation
par les deux modules du code de calcul que nous développons : l'un, de nature statistique, basé sur une
résolution des équations de la mécanique dans le domaine fréquentiel, l'autre, de nature instationnaire,
basé sur une résolution des mêmes équations, mais dans le domaine temporel. Le premier présentera
comme avantage essentiel d'être économique en terme de temps de calcul et facile à manipuler. Le
second, plus lourd, permettra notamment de mieux prendre en compte le caractère non linéaire de
l'évolution des coefficients aérodynamiques en fonction de l'incidence du vent et d’appréhender les
phénomènes d'instabilités aérodynamiques.
29
2.1 Forces aérodynamiques s’exerçant sur un corps mal profilé
2.1.1 Forces s’exerçant sur une section fixe plongée dans un flux laminaire
Considérons une section de tablier fixe, de largeur B et de hauteur H, placée dans un écoulement
bidimensionnel laminaire, stationnaire et caractérisé par un vecteur vitesse amont U compris dans le plan
de la section (fig.2.1). Selon la forme du profil et l’incidence par rapport à l’écoulement, des zones
caractéristiques dites de recirculation sont générées par un décollement du flux au niveau des angles
d’attaque et se développent de façon plus ou moins stable et marquée.
Selon les cas, cette zone peut ressembler à un bulbe, caractérisé par un point de recollement sur le profil,
comme elle peut se confondre avec la zone dite de sillage se développant en aval de la section. Entre les
zones de recirculation ou de sillage et l’écoulement libre se développe une couche de cisaillement
caractérisée par un fort gradient vertical de vitesse et un taux de turbulence élevé. Le sillage, lui, est
composé de tourbillons plus ou moins bien organisés dont la dimension caractéristique est fixée par H.
Couche
de cisaillement
Zone de
recirculation
Sillage
Couche de cisaillement
B H
Figure 2.1 – Ecoulement laminaire autour d’un corps non profilé :
lignes de courant et pressions pariétales moyennes
30
Pour un nombre de Reynolds3 suffisamment grand (>104), ces tourbillons peuvent être convectés de
manière périodique et sont alors baptisés tourbillons de Von Karman : on parle alors d’échappement
tourbillonnaire. La périodicité de ce phénomène, dont le processus de création est encore mal connu, est
caractérisé par le nombre de Strouhal :
fH
St  (2.1)
U
où H est la hauteur du tablier et f la fréquence caractéristique du détachement tourbillonnaire. Pour un

tablier de pont ce nombre varie entre 0,1 et 0,2. Nous reviendrons à ce phénomène plus loin.
L’intégration des valeurs de pression sur une tranche de tablier, de section constante, fournit la valeur
des efforts résultants s’exerçant sur cette tranche : un effort de traînée FD (drag) dans la direction du flux
moyen, un effort de portance FL (lift) perpendiculaire au flux et un moment de torsion M autour de l'axe
perpendiculaire au plan de la section. On introduit alors les coefficients adimensionnels suivants :
FD FL M
CD  CL  CM  (2.2)
1 1 1
 B LU 2  B LU 2  B 2 LU 2
2 2 2
où  représente la masse volumique de l'air, B une dimension caractéristique de la section (en général la
largeur), L la longueur de la tranche, et U la vitesse du flux incident.
CL
CY
y
+
CX
CD

x CM
U
+
Figure 2.2 – Repères de projection des coefficients aérodynamiques
3
Le nombre de Reynolds Re = UH/, où U représente la vitesse moyenne du flux, H une dimension transversale
caractéristique de l’obstacle, et  la viscosité du fluide, traduit le rapport entre forces d’inertie et forces dissipatives.
Sa valeur conditionne largement le régime de l’écoulement, en particulier lorsque l’obstacle est dépourvu d’arêtes
vives.
31
Ces coefficients, appelés coefficients aérodynamiques stationnaires, peuvent être mesurés en soufflerie
sur des maquettes sectionnelles (cf. chapitre 3). Ils dépendent de la forme de la section, de son incidence
par rapport au flux et du nombre de Reynolds de l'écoulement au voisinage de la section.
On sait toutefois que pour des nombres de Reynolds élevés, et pour des sections peu profilées, cette
dernière dépendance est assez faible. En conséquence, les efforts générés sur la section sont
proportionnels au carré de la vitesse U du flux incident.
0.4
0.3
Cx
Cy
0.2 Cm
0.1
0 Re (105)
0 0.5 1 1.5
Figure 2.3 – Influence du nombre de Reynolds sur les coefficients stationnaires

du tablier du pont de Saint-Nazaire
On préfère le plus souvent exprimer les forces et coefficients associés dans un repère lié à la section
étudiée. Les CD et CL sont alors convertis en CX et CY par une simple combinaison linéaire :
C X ()  C D cos   C L sin  (2.3)
C Y ()  C D sin   C L cos  (2.4)
Les figures 2.4 à 2.6 illustrent l’évolution de ces coefficients en fonction de l’incidence du vent pour
trois types de profils de tabliers. On peut remarquer que pour un tablier profilé, les coefficients
aérodynamiques varient quasi linéairement avec l'incidence. Pour un viaduc ferroviaire, beaucoup plus
massif, ou un bipoutre, moins bien profilé, cette linéarité n'est plus du tout vérifiée, en particulier pour
CY.
Comme on le verra par la suite, le calcul dynamique par approche spectrale nécessite une expression
linéaire des coefficients aérodynamiques en fonction de l’incidence du vent. Cela nous amènera à
proposer une méthode permettant de linéariser « intelligemment » des courbes fortement non linéaires.
32
0.2
0 
-10 -8 -6 -4 -2 0 2 4 6 8 10
-0.2
-0.4
Cx
Cy
-0.6
Cm
-0.8
-1
-1.2
Figure 2.4 – Coefficients aérodynamiques d’un tablier profilé : viaduc de Millau
0.8
Cx
Cy 0.6
Cm
0.4
0.2
0 
-10 -8 -6 -4 -2 0 2 4 6 8 10
-0.2
Figure 2.5 – Coefficients aérodynamiques d’un tablier non profilé : viaduc de Verrières
0.8
0.6
0.4
0.2
0 
-10 -8 -6 -4 -2 0 2 4 6 8 10
-0.2
Cx
-0.4 Cy
-0.6 Cm
-0.8
Figure 2.6 – Coefficients aérodynamiques d’un tablier bipoutre : pont du Tage
33
2.1.2 Forces s'exerçant sur une section fixe en flux turbulent: approche quasi-statique
On ne peut évidemment pas se contenter d’une modélisation des forces aérodynamiques ne s’appliquant
qu’au cas trivial d’un vent de vitesse et d’incidence constantes. Toutefois, dans un premier temps, on se
limite à une approche dite « quasi-statique » selon laquelle, à un instant t donné, les efforts exercés par
un signal de vent fluctuant sont les mêmes que ceux qui seraient exercés par un signal constant de même
amplitude à cet instant t.
Cependant, appliquer cette hypothèse brutalement, revient à négliger l'impact du degré de turbulence sur
la façon dont l'écoulement moyen contournera l'obstacle ; les positions des points de décollement et de
recollement pourront en effet en être très affectées et par là-même les efforts moyens s'exerçant sur
l'obstacle. Des essais sur maquettes sectionnelles effectués en vent turbulent permettent de régler ce
problème en remplaçant les coefficients moyens "laminaires" par des coefficients moyens "turbulents".
Mais appliquer l'hypothèse de quasi-stationnarité, revient surtout à négliger la nature même de

l'écoulement turbulent. Celui-ci peut en effet être assimilé à un ensemble de structures tourbillonnaires
d'échelles variables dont le comportement face à un obstacle sera également variable. Les structures de
grande taille par rapport à l'obstacle seront peu perturbées par son franchissement ; à l'inverse, les petites
structures seront sérieusement désorganisées.
Ainsi, les coefficients moyens sont incapables de différencier l'impact des différentes échelles
tourbillonnaires sur l'obstacle. En d'autres termes, ils ne permettent pas de fournir des informations
précises sur la façon dont l'énergie éolienne sera convertie en énergie sollicitante dans le domaine
fréquentiel.
Heureusement, l’essentiel de l’énergie éolienne est généralement concentré dans des structures
tourbillonnaires dont les dimensions caractéristiques sont de l’ordre de quelques dizaines de mètres c’est-
à-dire plusieurs fois la taille d’un élément de pont classique. On peut donc raisonnablement assimiler le
flux turbulent à un flux « localement laminaire variable » et, par conséquent, accepter l’hypothèse de
quasi-stationnarité dans un premier temps.
Mais il faudra bien garder à l’esprit qu’à mesure que le taux de turbulence augmente, d’une part, et que le
spectre de puissance se décale vers les hautes fréquences, d’autre part, cette hypothèse perd de sa
pertinence. Nous seront alors amenés à introduire une nouvelle notion corrective : celle d’ admittance
aérodynamique…
34
Mais restons pour le moment dans le cadre de cette hypothèse et considérons le vecteur vitesse
r ~
instantanée U( t ) qu'aurait le flux au centre de gravité de la section, en l'absence de cette dernière, U la
r
norme de la composante de U( t ) comprise dans le plan de la section et (t) l'angle que fait cette
composante avec l’axe x du repère lié à la section. On écrit alors:
1 ~
FX ( t )  BU ( t ) 2 C X (( t )) (2.5)
2
1 ~
FY ( t )  BU ( t ) 2 C Y (( t )) (2.6)
2
1 ~
M Z ( t )  B 2 U( t ) 2 C M (( t )) (2.7)
2
r r r
On exprime à présent les composantes du vecteur vent dans une base ( x , y, z) liée à la structure, où
r r r
( x , y) définit une base du plan contenant la section considérée. Ainsi, le vecteur vitesse instantané U ( t )
s'écrit:
r r r r
U( t )  U L ( t ) x  VL ( t ) y  WL ( t )z (2.8)
on peut donc écrire:

 V (t ) 
( t )  arctan L  (2.9)
 U L (t) 
~
U( t )  U L ( t ) 2  VL ( t ) 2 (2.10)
~ r
En introduisant U à la place de U on néglige donc tous les effets de la composante WL(t) .
Par analogie avec la relation (1.1), le vecteur vitesse instantané peut s'écrire sous la forme:
r
U( t )  U L  u L ( t ), VL  v L ( t ), WL  w L ( t )  (2.11)
Les équations (2.5-7) deviennent donc:
Fx ( t ) 
1
2
 2 2

 U L  u L ( t )   VL  v L ( t )  BC x (( t )) (2.12)
2
 2 2

FY ( t )   U L  u L ( t )   VL  v L ( t )  BC Y (( t ))
1
(2.13)
2
 2 2

M Z ( t )   U L  u L ( t )   VL  v L ( t )  B 2 C M (( t ))
1
(2.14)
35
 V  v L (t ) 
( t )  arctan L  (2.15)
 U L  u L (t ) 
On suppose que la vitesse du vent "vue" par le tablier est la même que celle qu'on aurait en absence du
tablier ; or les tourbillons se déforment en entrant en contact avec ce dernier. Le champ de pression qui
en découle va alors être longitudinalement plus cohérent que le champ de vent incident. Ce phénomène
est négligé ici faute de modèles mathématiques appropriés.
La méthode spectrale, présentée plus loin, dont l’objet est la détermination des effets de la turbulence sur
les structures élancées nécessite que la relation liant les forces aérodynamiques au champ de vent
turbulent soit linéaire. Or, les intensités de turbulence, dans le cas de l'écoulement atmosphérique, sont
généralement inférieures à 20% donc uL et vL sont petits par rapport à U L , VL . On peut donc linéariser
les expressions ci-dessus au voisinage de u = 0 et v = 0 :
dC X
C X ( )  C X (  0 )  (  0 ) (2.16)
d
dC Y
C Y ( )  C Y (  0 )  ( 0 ) (2.17)
d
dC M
C M ( )  C M ( 0 )  (  0 ) (2.18)
d
avec   ( t )   0 .
on peut alors développer (t) comme suit :
VL  
( t )   0    arctan( ) (0,0) u L  (0,0) v L (2.19)
UL u L v L
 VL  UL
avec (0,0)   et (0,0)  (2.20)
u L U L  VL
2 2
v L U L  VL
2 2
En ne gardant que les termes du premier ordre, on obtient:
FX ( t )  FX  FtX ( t ) (2.21)
FY ( t )  FY  FtY ( t ) (2.22)
M Z ( t )  M Z  M tZ ( t ) (2.23)
36
1
FX  ( U L  VL ) BC x ( 0 )
2 2
avec (2.24)
2
1 C
FtX ( t )  ( U L u L ( t )  VL v L ( t ))BC X ( 0 )  B X ( 0 )( U L v L ( t )  VL u L ( t )) (2.25)
2 
1
FY  ( U L  VL )BC Y ( 0 )
2 2
(2.26)
2
1 C
FtY ( t )  ( U L u L ( t )  VL v L ( t ))BC Y ( 0 )  B Y ( 0 )( U L v L ( t )  VL u L ( t )) (2.27)
2 
1
MZ  ( U L  VL )B 2 C M ( 0 )
2 2
(2.28)
2
1 C M
M tZ ( t )  ( U L u L ( t )  VL v L ( t ))B 2 C M ( 0 )  B 2 ( 0 )( U L v L ( t )  VL u L ( t )) (2.29)
2 
Ainsi, dans la formulation linéarisée des forces aérodynamiques, chaque composante est la somme de
deux termes:
 un premier terme constant qui dépend des caractéristiques du vent moyen, de son incidence par
rapport à la section et de la valeur des coefficients aérodynamiques à l'incidence moyenne ;
 un second terme fluctuant proportionnel aux composantes turbulentes du vent couramment appelé
« buffeting force » ; c’est ce terme qui est à l’origine du comportement dynamique de la structure ;
c’est donc lui qui, dorénavant, retiendra l’essentiel de notre attention.
Par la suite, on exprimera les équations (2.25), (2.27), et (2.29) sous la forme matricielle :
  1 dC X   1 dC X  
 B C X U L  VL  B C X VL  UL  0
 FtX ( t )    2 d    2 d   
     
 FtY ( t )   1 dC 1 dC 
 B C Y U L  Y
VL  B C Y VL  Y
UL  0
 F (t )   2 d   2 d   u L (t) 
 
 tZ  0 0

0  v L (t)  (2.30)
 M tX ( t )    w (t ) 
   0 0 0  L 
 M tY ( t )   
 M (t )   0 0 0
 tZ   2  1 dC M   1 dC M  
B  C M U L  2 d VL  B  C M VL  2 d U L 
2
0
     
que l’on contractera en F t ( t )  C ad u L ( t ) .
Lorsque la section est correctement profilée, on a vu que les coefficients aérodynamiques variaient
quasiment linéairement avec l'incidence du vent ; dans ce cas, la formulation linéarisée des forces
37
aérodynamiques est pleinement justifiée. Pour les sections mal profilées, notre formulation est par contre
moins représentative de la réalité. Ainsi, à titre d’exemple, l’annulation de la dérivée du coefficient de
portance du tablier de Verrières (figure 2.5) autour de l’incidence 1° est peu significative des variations
du coefficient au voisinage de ce point, sauf à considérer d’infimes variations d’incidences autour de
cette position.
L’expérience montre cependant que l’incidence instantanée du vent varie de plusieurs degrés autour de
l’incidence moyenne. C’est pourquoi nous proposons une méthode permettant de remplacer les courbes
d’évolution des coefficients aérodynamiques en fonction de l’incidence du vent par des fonctions affines
équivalentes vis-à-vis de la réponse dynamique de nos structures. Nous y reviendrons dans le chapitre 4.
Les coefficients aérodynamiques stationnaires pourraient, si l’on poussait à fond la logique de

l’hypothèse de quasi-stationnarité, être mesurés dans des conditions d’écoulement laminaire puis
appliqués dans des cas d’écoulement turbulent. Cependant, comme on l'a déjà évoqué, on préfère mesurer
ces coefficients en régime turbulent. La figure 2.7 montre, pour le tablier du viaduc de Millau, les
coefficients CX, CY, CM, mesurés en conditions laminaires, et ceux mesurés en écoulement turbulent. On
peut remarquer que les différences sont faibles.
Il nous semble important d'appuyer sur le fait que, même dans cette configuration expérimentale, les
coefficients mesurés sont de simples rapports entre pressions dynamiques moyennes amont et efforts
moyens mesurés. Cela signifie que l’on rend compte, dans l’estimation des coefficients stationnaires
« turbulents », des caractéristiques moyennes de l’écoulement propres à son taux de turbulence (allures
moyennes des zones de recirculation, sillage,…) mais pas de l’impact différencié des différentes échelles
de turbulence.
0,2
0
-10 -8 -6 -4 -2 0 2 4 6 8 10
-0,2
-0,4 Cx uniforme
Cy uniforme
-0,6
Cm uniforme
-0,8 Cx turbulent
Cy turbulent
-1 Cm turbulent
-1,2
Figure 2.7 – Coefficients stationnaires du viaduc de Millau mesurés en vents uniforme et turbulent
38
De plus, même plongé dans un écoulement parfaitement laminaire, un corps mal profilé est soumis à une
turbulence générée par lui-même : son sillage ou « signature » aérodynamique.
Ces notions sont prises en compte via l’introduction d’une fonction corrective s’exprimant dans le
domaine fréquentiel : la fonction d’admittance aérodynamique. Cette fonction prend la valeur unité pour
les tourbillons de basses fréquences, et décroît à mesure qu’on se rapproche des échelles de hautes
fréquences. On modélise ainsi, pour une section donnée, l'évolution de la corrélation spatiale d'un champ
de pressions dynamiques en fonction de la taille des structures tourbillonnaires.
Ainsi, en réécrivant (2.30) dans le domaine spectral, et en introduisant les termes  ij () ( i = FX, FY, MZ,
j = u, v ), on obtient :
  1 dC X   1 dC X  
 B C X U L  VL   Xu () B C X VL  U L   Xv () 0
 FtX ()    2 d   2 d  
     
 FtY ()   1 dC 1 dC 
B C Y U L  Y
VL   Yu () B C Y VL  Y
U L   Yv () 0
 F ()    2 d   2 d   u L () 
 
 tZ  0 0

0  v L () 
 M tX ()    w () 
   0 0 0  L 
M
 tY ( )   
 M ()   0 0 0
 tZ   2 1 dC M  2 1 dC M  
B  C M U L  2 d VL   Mu () B  C M VL  2 d U L   Mv () 0
     
(2.31)
que l’on contractera en F t ()  C  () u L () .
Considérons deux points P et P’ d’une structure vibrante ; les fonctions d’admittance ij,kl’() sont les
produits conjugués des fonctions  ij () au point P et  kl () au point P’; elles apparaissent dans
l’expression reliant la matrice de densité interspectrale des forces à la matrice de densité interspectrale de
la turbulence pour cette paire de points :
S F ()  C  () S u , u ' () C  ' ()

t*
(2.32)
Dans le cas des ponts, seul le tablier, en raison de sa profondeur d'une part (désynchronisation des
tourbillons qui transitent) et de sa sensibilité d'autre part (forte contribution en terme de forces
généralisées), est généralement affecté de fonctions d’admittance aérodynamique. De plus, seules les six
fonctions d’admittance aérodynamiques « non couplées » du type ij() ( i = FX, FY, MZ, j = u, v ) sont
supposées non nulles.
39
Elles sont même souvent approchées par une seule et même fonction empirique. Plusieurs formulations
peuvent être utilisées (figure 2.8) :
 7 B 
 7 B 

2  1  e  2 U  
 2 U 
  
2
Davenport (1962) : 2
(2.33)
 7 B 
 
 2 U 
1
  
2
Vickery (Holmes,1975): (2.34)
2 B
1
 U
1
  
2
Liepmann (1952): (2.35)
B
1 
U
L'inconvénient majeur de ces formulations est qu'elles ne sont pas associées à un type de tablier
particulier. Il est donc difficile de faire un choix a priori. C'est pourquoi il est préférable de les mesurer
en soufflerie via une procédure que l’on explicitera dans le chapitre suivant.
1,0
Admittance aérodynamique
Vickery
0,8
Liepmann
0,6
Davenport
0,4
0,2
0,0
0 1 2
Fréquence réduite nB/U
Figure 2.8 – Modèles de fonctions d’admittance aérodynamique
Dans le domaine temporel, l’admittance aérodynamique intervient naturellement sous forme de produits
de convolution :
40
  1 dC X   1 dC X  
 B C X U L  VL   Xu () B C X VL  U L   Xv () 0
 FtX ( t )    2 d   2 d  
     
 FtY ( t )   1 dC Y 1 dC Y 
 F (t )   B C X U L  2 d VL   Yu () B C Y VL 
2 d
U L   Yv () 0
 u L ( t  ) 
t
       
 tZ 
 M tX ( t )  
0


0 0 0
0
 v L ( t  ) d
 w ( t  ) 
  0 0  L 
 M tY ( t )   
 M (t )   0 0 0
 tZ   2  1 dC   1 dC M  
B  C M U L  2 d VL   Mu () B  C M VL  2 d U L   Mv ()
M 2
0
     
t
que l’on contracte en

F t ( t )  C  () u L ( t  ) d
0
(2.36)
Les fonctions  ij ( t ) , transformées de Fourier inverses des fonctions  ij () , sont généralement
nommées fonctions indicielles. Des fonctions similaires interviendront dans la modélisation des forces
aéroélastiques. L’hypothèse selon laquelle toutes les fonctions d’admittance sont égales permet de
factoriser (2.36) :

F t ( t )   () C ad u L ( t  ) d
0
(2.37)
41
2.1.3 Section mobile, flux incident variable : approche quasi-statique
On considère maintenant une section en mouvement dans un flux turbulent. On suppose que le
mouvement modifie peu la structure du fluide environnent. On note x(t), y(t) et (t) les déplacements
horizontal, vertical, et en torsion de la section considérée.
Dans ce cas, l’équation 2.15 devient :
 V  v( t )  y& ( t ) 
( t )  arctan   ( t ) (2.38)
 U  u ( t )  x& ( t ) 
ce qui nous donne, par linéarisation, en supposant x& et y& petits devant U et V :
FX ( t )  FX  FXt ( t )  FXa ( t ) FY ( t )  FY  FYt ( t )  FYa ( t ) FM ( t )  FM  FMt ( t )  FMa ( t ) (2.39)
où les deux premiers termes de chacune des trois expressions sont identiques à ceux développés en 2.24-
29 et où les termes supplémentaires, correspondant aux forces issues des déplacements du tablier et dits
« aéroélastiques » s’expriment comme suit :
1 C X
Faex ( t )  B (  )( V L x& ( t )  U L y& ( t )   ( U L  V L ))
2 2
(2.40)
2 
1 C
FaeY ( t )  B Y (  )( VL x& ( t )  U L y& ( t )  ( U L  VL ))
2 2
(2.41)
2 
1 2 C M
FaeM ( t )  B (  )( VL x& ( t )  U L y& ( t )  ( U L  VL ))
2 2
(2.42)
2 
Ces forces, proportionnelles à des déplacements ou vitesses de déplacements, s’apparentent donc à des
raideurs et amortissements aérodynamiques. Cette première approche des forces aéroélastiques reste
correcte tant que l'hypothèse de quasi-stationnarité est vérifiée. Autrement dit, les formulations (2.40-42)
sont valables pour modéliser les déplacements de basses fréquences, mais perdent de leur pertinence à
mesure que l'on s'approche des hautes fréquences de vibration. C'est pourquoi d'autres modèles ont été
développés pour représenter les forces aéroélastiques sur l'ensemble du spectre des déplacements.
42
2.2 Les forces aéroélastiques
En première approximation, on a pu tenir compte du mouvement de la structure en remplaçant la vitesse

du vent par la vitesse relative par rapport à la section. Néanmoins ceci est insuffisant car les oscillations
de la section modifient la géométrie de l’écoulement et par conséquent les forces aérodynamiques.
Les premières formulations permettant de prendre en compte les phénomènes aéroélastiques sont issues
du domaine de l’aéronautique et font intervenir des fonctions dites indicielles (travaux de Wagner en
1925 et Theodorsen en 1935). Ces fonctions permettent d’exprimer l’incidence d’une perturbation de
type Dirac dans le signal des déplacements sur le signal des forces aérodynamiques. La prise en compte,
à un instant t, de tout l’historique des déplacements sur la valeur du signal des forces se fait via un simple
produit de convolution. Ces formulations ont été adaptées plus tard aux tabliers de ponts, notamment par
Bleich (1948), Bisplinghoff (1962), Scanlan (1971,1974) et Lin (1983).
Les formulations évoquées ci-dessus présentent l’inconvénient d’être difficilement appréhensibles

expérimentalement. Scanlan et Tomko (1971) ont proposé une formulation dite mixte, car s’exprimant à
la fois dans les domaines temporel et fréquentiel, ayant l’intérêt majeur d’être basée sur une famille de
modèles expérimentaux clairement définis : maquette sectionnelle rigide à un, deux ou trois degrés de
liberté en oscillations sinusoïdales. Cette formulation (ou plutôt ce type de formulation car il en existe
d’autres analogues), bien que non universelle, constitue néanmoins le nœud essentiel entre la
caractérisation expérimentale des propriétés aéroélastiques d’un tablier et le calcul numérique de son
comportement en vraie grandeur.
Ainsi, nous aborderons notre présentation de la modélisation des forces aéroélastiques par le modèle
proposé par Scanlan et Tomko. Nous montrerons ensuite comment généraliser ce modèle afin de
l’intégrer aux méthodes de calcul spectrale et temporelle.
2.2.1 Formulation mixte temps-fréquence : Scanlan et Tomko (1971)
Considérons un tablier de section constante à deux degrés de liberté : flexion + torsion de pulsations
respectives 1 et 2, placé dans un flux laminaire. Lorsque les deux modes sont excités, on observe
expérimentalement l’apparition de forces supplémentaires, proportionnelles aux déplacements et vitesses
de déplacements, et donc assimilables à des termes de raideur et d'amortissement. Ceci se traduit alors,
dans le cas d’oscillations libres, par une modification de la décroissance des signaux de déplacements et
un glissement en fréquence de ceux-ci.
43
Ces raideur et amortissement sont fonctions, outre de la forme de la section et de son incidence par
rapport à l’écoulement, de la vitesse du flux et de la pulsation des modes excités. Généralement, on
considère uniquement la dépendance vis-à-vis de la pulsation réduite Ki appelée couramment fréquence
réduite, définie comme suit :
i B
Ki  i = 1, 2 (2.43)
U
où i est la pulsation des oscillations, B une dimension caractéristique du tablier (en général sa largeur),
U la vitesse du vent. Une petite analyse dimensionnelle montre que cette variable représente le nombre
d’oscillations effectuées par le tablier pendant le trajet d’une particule d’air entre le bord d’attaque et le
bord de fuite. Abstraction faite de l’influence du nombre de Reynolds sur la géométrie de l’écoulement,
la dépendance au nombre K relève physiquement du bon sens.
Bien que le mouvement horizontal du tablier soit également soumis à de tels effets, on ne s'intéresse en
général qu'aux effets aéroélastiques de portance et de moment, qui sont les plus significatifs. En
supposant que ces forces supplémentaires dépendent linéairement des déplacements, on peut les écrire
sous la forme suivante:
1 2  
L ae   U B (H 1 h  H 2   H 3   H 4 h ) (2.44)
2
1 2  
M ae  U B 2 (A 1 h  A 2   A 3   A 4 h ) (2.45)
2
H 1 traduit l'amortissement aérodynamique en flexion, A 2 l'amortissement aérodynamique en torsion,
H 4 le glissement en fréquence en flexion, et A 3 le glissement en fréquence en torsion. Les autres termes
A 1 , A 4 , H 2 , H 3 traduisent les phénomènes de couplage entre modes.
On introduit généralement les coefficients adimensionnels H *i et A *i définis ci-dessous :
2
K K B 2 K
H 1  1 H 1* (K 1 ) H 2  2 H *2 (K 2 ) H3  K 2 H *3 (K 2 ) H 4  1 H *4 (K 1 ) (2.46)
U U B
2
K1 * K 2B * 2 K1
A1  A1 (K 1 ) A2  A 2 (K 2 ) A 3  K 2 A *3 (K 2 ) A4  A *4 (K 1 ) (2.47)
U U B
44
Ces nouveaux termes sont appelés coefficients aéroélastiques. Insistons sur le fait qu’ils dépendent de la
forme de la section, de son incidence, de la fréquence réduite, et … du nombre de Reynolds.
Remarquons que dans les publications de Scanlan sur le sujet (Scanlan, 1974), la variable K n’est
affectée d’aucun indice, ce qui pose régulièrement des problèmes dans les interprétations : rien n’indique
si la fréquence correspond au mouvement de flexion ou de torsion. Nous avons donc opté pour la
convention qui nous semblait la plus logique.
2.2.2 Généralisation à un mouvement quelconque de la section
a) Formulation fréquentielle
Les formules 2.44 et 2.45 ne fournissent une modélisation des effets aéroélastiques que pour un
mouvement purement sinusoïdal. On aimerait cependant généraliser cette formulation à un mouvement
quelconque.
Notons encore h(t) et (t) les déplacements verticaux et en torsion du même tablier à deux degrés de
liberté et supposons qu'ils admettent chacun une transformée de Fourier ĥ () et ˆ () . Considérons à
~ ~ ( t ) , d'amplitudes ĥ(  ) , ˆ (  ) , et de
présent deux mouvements purement sinusoïdaux h m ( t ) ,  n m n
pulsations  m ,  n.
~ ~ ( t )  ˆ ( ) e in t ces déplacements génèrent les

D’après 2.44 et 2.45, en posant h m ( t )  ĥ ( m ) e im t ,  n n
efforts aéroélastiques suivants:
~ ~
2 
~ 1 2 h (t) 2 ~ ( t )  K 2 H * ( K ) h m ( t ) 
~ ( t )  K 2 H * ( K )
L ae mn ( t )   U B i K m H 1* ( K m ) m  i K n H *2 ( K n ) 
B 
n n 3 n n m 4 m
2  B
(2.48)
~ ~
~ 1 2 2 h i (t ) h (t ) 
 2 *
M ae mn ( t )  U B i K i A 1 (K i )  i K j A 2 (K j )  j ( t )  K j A 3 (K j )  j ( t )  K i A 4 (K i ) i 
2 * ~ 2 * ~ 2 *
2  B B 

(2.49)
où Kn = nB/U et Km = mB/U.
45
Décomposant chacune des deux formulations ci-dessus en deux parties, suivant les indices i et j, on peut
considérer les quatre forces distinctes suivantes:
~ ~
~ 1 2  h (t) h (t) 
L ae m ( t )   U B  i K m H 1* ( K m ) m  K m H *4 ( K m ) m 
2 2
(2.50)
2  B B 
~ 1
2
2
 2 ~ ( t )  K 2 H * ( K )
L ae n ( t )   U B i K n H *2 ( K n )  n n 3 n
~ (t)
n  (2.51)
~ 1
2
2
2 ~ ( t)  K 2 A* ( K ) 
M ae n ( t )   U B 2 i K n A *2 ( K n )  n n 3 n
~ (t)
n  (2.52)
~ ~
~ 1 2  h (t) h (t) 
M ae m ( t )   U B 2  i K m A 1* ( K m ) m  K m A *4 ( K m ) m 
2 2
(2.53)
2  B B 
On remplace aisément les  par des  en divisant tous les termes de (2.50) et (2.51) par e im t et ceux
de (2.52) et (2.53) par e in t arrivant ainsi à une formulation fréquentielle des efforts aéroélastiques.
Considérant que h et  résultent de la sommation des signaux sinusoïdaux d’amplitudes ĥ , ̂ sur

l’ensemble du domaine fréquentiel, on fait la supposition que L ae , M ae résultent également de
l’intégration des signaux d’amplitudes L̂ ae , M̂ ae sur la même plage de fréquences: "principe de

superposition. On peut légitimement mettre en doute la validité de cette hypothèse ; en effet, les
équations de Navier-Stokes régissant l'écoulement du flux autour de la section n’ayant pas un caractère
linéaire, il n'y a aucune raison que les effets s'additionnent aussi simplement. Si on l'applique tout de
même, L̂ ae , M̂ ae sont alors les transformées de Fourier de L ae , M ae .
En introduisant les fonctions suivantes:
1 3 2 1
FL ()  B  (H3*  iH 2 * ) FLh ()  B 2  2 (H 4 *  iH1* ) (2.54)
2 2
1 1 3 2
FM ()  B 4  2 (A3*  iA 2 * ) FMh ()  B  (A 4 *  iA1* ) (2.55)
2 2
on pose alors:
L̂ ae ()  FLh ()ĥ ()  FL ()ˆ () (2.56)
M̂ ae ()  FMh ()ĥ ()  FM ()ˆ () (2.57)
46
Ainsi, la transformée de Fourier des déplacements fournit la transformée de Fourier des forces
aéroélastiques. Ceci se généralise facilement au cas où h et  sont deux processus aléatoires admettant
une fonction de densité spectrale de puissance.
b) Formulation temporelle des forces aéroélastiques
Les expressions (2.56) et (2.57), modélisant les forces aéroélastiques dans le domaine fréquentiel sous
forme de simples combinaisons linéaires, se transposent aisément dans le domaine temporel via
l’introduction de produits de convolution.
Soient f Lh , f L , f Mh , f M les transformées de Fourier inverses des fonctions FLh , FL , FMh , FM .
On écrit alors:
t t
L ae ( t ) 


f Lh ( t  )h ( )d 
f

L ( t  )( )d  L h ( t )  L  ( t ) (2.58)
t t
M ae ( t ) 

f Mh ( t   ) h (  ) d 
f

M ( t  )( )d  M h ( t )  M  ( t ) (2.59)
où Lh et L sont les forces de portance générées par les mouvements vertical et angulaire ; Mh et M sont
les moments générés par les mouvements vertical et angulaire.
Signalons qu’on retrouve là des expressions analogues à celles développées par Lin (1983), en cela
qu’elles font intervenir des fonctions de type indiciel. Toutefois, les formulations de Lin ne font pas
intervenir les fonctions déplacement mais leurs dérivées par rapport au temps. Ainsi, par exemple, la
force de portance s’exprime sous la forme suivante :
t t
L ae ( t ) 


g Lh ( t  )h& ()d 
g

L ( t  )& ()d (2.60)
Les fonctions indicielles gLh et gL ont le mérite d’avoir un sens physique clair.
En effet, gLh(t) h = g Lh (t) h& (0) t ( resp. gL (t)  = g L (t) & (0) t ) représente la force élémentaire
générée à l’instant t par une impulsion élémentaire h (resp. ) à l’instant t=0. Les efforts totaux
générés à un instant t s’obtiennent naturellement par intégration des efforts élémentaires, l’instant de
l’impulsion variant de - à t. Notre propos est illustré par les figures 2.9 et 2.10.
47
Si le sens physique des équations (2.58) et (2.59) paraît a priori moins évident, on remarque toutefois
qu’une simple intégration par parties des équations de Lin permet de retomber sur des équations de la
forme de celles issues de la transformée de Fourier des équations de Scanlan. Ce qui est plutôt rassurant.
0 t
Figure 2.9 – Allure d’une fonction indicielle
"forces" indicielles
déplacements force résultante
Figure 2.10 – Intégration des fonctions indicielles
Les coefficients H i et A i mesurés en soufflerie suffisent donc, en théorie, pour estimer les fonctions
indicielles permettant la prise en compte des efforts aéroélastiques dans le domaine temporel.
Cependant, afin de réduire les temps de calcul, on préfère ne pas calculer numériquement les
transformées de Fourier inverses des fonctions FLh , FL , FMh , FM puis le produit de convolution. On
approche donc ces quatre fonctions par des applications dont on connaît les transformées de Fourier
analytiques, suivant la forme proposée par Bucher (1988).
48
Chaque fonction F est paramétrée par six coefficients C1, C2, C3, C4, d3, d4, et a la forme suivante:
 
 4 
i
C
1 B
FLh ()  U 2  C 1Lh  iC 2 Lh   (2.61)
2  U
kLh
U 
 k 3
d kLh  i 
 B 
 
 4 
i
C
1 2 B 
FL ()  BU C 1L  iC 2 L  (2.62)
2  U
kL
U 
 k 3
d kL  i 
 B 
 
 4 
i
C
1 B
FMh ()  BU 2  C1Mh  iC 2 Mh   (2.63)
 
kMh
2 U U
 k 3 d kMh  i 
 B 
 
 4 
i
C
1 B
FM ()  B 2 U 2  C 1M  iC 2 M   (2.64)
2  U
kM
U 
 k 3
d kM  i 
 B 
Alors, en appliquant les transformées de Fourier inverses 4, on obtient :
 

2
 L h 3 ( t )  L h 4 ( t ) 
1 h(t ) B h (t )
L h ( t )  BU  C 1Lh  C 2 Lh (2.65)
2  B U B 
 
1  B 

L  (t )  BU 2  C 1L ( t )  C 2 L ( t )  L 3 ( t )  L  4 ( t )  (2.66)
2  U 
 

2
 M h 3 ( t )  M h 4 ( t ) 
1 h(t) B h(t)
M h ( t )  B U  C 1Mh
2
 C 2 Mh (2.67)
2  B U B 
 
1 2 2 B 

M  (t )  B U  C1M ( t )  C 2 M ( t )  M  3 ( t )  M  4 ( t )  (2.68)
2  U 
où on pose, avec k = 3, 4 :
4
Dans le cas présent, afin d’être en cohérence avec les formulations (2.54-55), la transformée de Fourier est définie
au facteur 1/2 près.
49
t  t
d Lhk U d U
h t      Lk  t   
 

L h k ( t )  C Lhk () e B
d L  k ( t )  C Lk ()e B
d (2.69)
B
 
t  t
d Mhk U d U
h t     Mk  t   
 

M h k ( t )  C Mhk () e B
d M  k ( t )  C Mk  () e B
d (2.70)
B
 
On montre facilement que les fonctions L h k ( t ), L  k ( t ), M h k ( t ), M  k ( t ) vérifient les relations suivantes:

 d U h  d Lk U 
L h k ( t )   Lhk L h k ( t )  C Lhk ( t ) L k (t)   L  k ( t )  C Lk ( t ) (2.71)
B B B

 d U h  d Mk U 
M h k ( t )   Mhk M h k ( t )  C Mhk ( t ) M k (t )   M  k ( t )  C Mk ( t ) (2.72)
B B B
Pour calculer les coefficients Ck et dk, on remanie les expressions (2.61-64) (resp.(2.54-55)) en
remplaçant les variables B, U et  par K (resp. les variables B et  par K et U) et on identifie terme à
terme ces expressions pour différentes valeurs de K. Cela donne :
 C Lhk d Lhk 1 K    
4 4
 1 d
1 1 C Lhk
H 1* K    KC Lh 2   H *4 K    C Lh1   (2.73)
K2  1  d Lhk 1 K 2    1 K  
2
 k 3
K2  k 3 Lhk
2 2
 C Lk d Lk 1 K    
4 4
 1 d
1 1 C Lk
H K   2
*
 KC L 2   H K   2
*
 C L1   (2.74)
 1  d Lk 1 K 2    1 K  
2 2 3 2 2
K  k 3
K  k 3 Lk
 C Mhk d Mhk 1 K    
4 4
 1 d
1 1 C Mhk
A 1* K    KC Mh 2   A *4 K    C Mh1   (2.75)
K2  1  d Mhk 1 K 2    1 K 2  
2
 k 3
K2  k 3 Mhk
2
 C Mk d Mk 1 K    
4 4
 1 d
1 1 C Mk
A *2 K    KC M 2   A *3 K    C M1   (2.76)
K2  1  d Mk 1 K 2    1 K  
2
 k 3
K2  k 3 Mk
2 2
On résout alors ce système d’équations en utilisant la méthode d’optimisation non linéaire par les
moindres carrés de Levenberg-Marquardt.
50
L'approximation des fonctions F par des formulations analytiques permet certes de manipuler des
expressions des forces aéroélastiques où n'apparaissent pas de produits de convolution, mais, en
contrepartie, nécessitent d’ajouter des équations différentielles supplémentaires à l'équation du
mouvement.
La forme temporelle des forces aéroélastiques donne a priori les efforts pour un mouvement quelconque.
Il ne faut pas pour autant oublier les nombreuses hypothèses qui ont été nécessaires pour exprimer ces
forces aéroélastiques:
 petite amplitude du mouvement permettant de supposer l'incidence constante,

 application du principe de superposition,
 flux incident constant.
Or, nous nous intéressons, dans ce qui suit, aux effets du vent turbulent sur les structures élancées ;
l'amplitude du mouvement n'est pas forcément faible et surtout, l'incidence n'est pas constante. On
supposera donc que les effets aéroélastiques, dans le cas d'un écoulement incident turbulent, sont
identiques à ceux créés par l'écoulement moyen autour de l'incidence moyenne.
En outre, on pourrait essayer d'estimer l'influence de la turbulence sur les effets aéroélastiques en
appliquant la même technique d'identification des coefficients H *i , A *i sur des signaux obtenus en
écoulement turbulent et en écoulement laminaire.
51
2.3 Le détachement tourbillonnaire
On a vu, au début de ce chapitre, que, dans certaines conditions, pouvait se produire un phénomène,
appelé détachement tourbillonnaire, caractérisé par un lâché de tourbillons de part et d’autre du profil,
alternés (i.e. en opposition de phase), à une fréquence n bien précise reliée à la vitesse U de
fSH
l’écoulement et à la hauteur du tablier par le nombre de Strouhal : St  .
U
Ce phénomène se traduit par la génération d’efforts fluctuants de fréquence 2n dans la direction du vent
et n dans la direction orthogonale. Les premiers sont généralement négligés tandis que les seconds font
l’objet de plusieurs formulations empiriques plus ou moins complexes.
Amplitude des Amplitude des

déplacements déplacements
t t
Spectre de puissance des Spectre de puissance des

déplacements déplacements
f propre
f dét. tourb. = f Strouhal f dét.tourb.= f propre
n n
Figure 2.11 - Signaux de déplacements verticaux lorsque la fréquence de détachement tourbillonnaire

est proche d’une fréquence propre structurelle : phénomène de synchronisation forcée
52
Dans le cas où la fréquence naturelle de détachement tourbillonnaire est suffisamment proche de la

fréquence propre d’un mode transversal de la structure considérée (fig.2.11), on assiste à un phénomène
de synchronisation forcée de la première sur la deuxième. Si l’écart est un peu plus important mais reste
suffisamment faible on observe sur le signal des déplacements le phénomène bien connu de
« battement ».
Il nous semble important de préciser que le signal des déplacements liés au détachement tourbillonnaire
n’a pas, en général, un caractère divergent ; l’expérience a montré que son amplitude dépassait rarement
la moitié de la hauteur du profil considéré. Il est clair que les oscillations maximales seront observées
pendant le phénomène de synchronisation forcée ; c’est donc dans ces conditions sécuritaires que seront
effectués les calculs de déplacements.
Différentes formulations empiriques, plus ou moins élaborées, ont été proposées afin de modéliser les
efforts dus au détachement tourbillonnaire. Cependant, celles-ci concernent essentiellement des formes
simples telles que les cylindres. Il n’existe pas véritablement de modèles permettant de représenter les
efforts exercés par le détachement tourbillonnaire sur des formes complexes telles qu’une section de
tablier de pont. Les effets du détachement tourbillonnaire ne seront donc pas pris en compte dans nos
calculs.
2.4 Conclusion
Nous avons, tout au long de ce chapitre, introduit un certain nombre de modèles mathématiques
permettant, à partir de la connaissance du vent, de représenter les forces aérodynamiques s’exerçant sur
une structure souple et élancée. Les buffeting forces, liées à l’action turbulente du vent ont été
modélisées à partir des coefficients stationnaires et de l’admittance aérodynamique. Les forces
aéroélastiques, liées aux oscillations de l’ouvrage, ont été modélisées par les coefficients instationnaires
de Scanlan dans le domaine spectral et par des fonctions indicielles dans le domaine temporel.
L’intérêt majeur des modèles présentés est de réduire la notion de forces aérodynamiques à un nombre
assez limité de coefficients et, surtout, des relations simples avec le champ de vitesse du vent. Cela
permet de s’affranchir de la résolution des équations de la mécanique des fluides lors de l’estimation de
la réponse dynamique de l’ouvrage ; on réduit ainsi considérablement les temps de calcul.
53
Chapitre 3 Mesure
des forces
aérodynamiques
3 Mesure des forces aérodynamiques
Dans le chapitre précédent, nous avons défini et modélisé les forces aérodynamiques susceptibles de
solliciter une structure élancée. Pour cela, nous avons introduit une série de coefficients adimensionnels
permettant de caractériser les propriétés aérodynamiques d’une section soumise à un écoulement
aéraulique bidimensionnel. Ceux-ci permettent donc, à partir de la connaissance du vent, des dimensions
de la structure et de l’historique de son comportement dynamique, de calculer les forces aérodynamiques
s’exerçant sur chaque section élémentaire de la structure.
Ces coefficients peuvent être estimés expérimentalement et/ou numériquement. Les méthodes
expérimentales sont basées sur la fabrication de maquettes rigides, de sections continues, appelées
maquettes sectionnelles, placées en soufflerie. Les méthodes numériques sont basées sur la résolution 2D
ou 3D des équations de Navier-Stokes dans une veine virtuelle contenant le profil à étudier. Ces
méthodes, en plein essor, ne sont pas encore réellement opérationnelles pour des applications au
domaine des ponts. Aussi nous consacrons ce chapitre exclusivement aux différentes techniques
expérimentales permettant de caractériser aérodynamiquement les tabliers de pont.
Ces techniques sont appliquées à deux types de profils : un profil "théorique" rectangulaire de ratio
largeur/épaisseur égal à huit et un profil "classique" trapézoïdal proche de celui du viaduc de Millau. Les
résultats présentés seront utilisés par la suite lors de la procédure de validation de notre code de calcul.
54
Chapitre 3 Mesure des forces aérodynamiques
3.1 Mesure des coefficients stationnaires
Les coefficients stationnaires sont généralement mesurés sur des maquettes sectionnelles fixées à deux
balances dynamométriques. Chaque balance est équipée de 16 jauges d’extensométrie, permettant de
mesurer 8 forces horizontales et 8 forces verticales. On obtient les coefficients de trainée, portance et
moment globaux par combinaison linéaire des trente-deux voies et adimensionalisation.
Figure 3.1 – Balance dynamométrique
Une deuxième technique consiste à équiper la maquette sectionnelle de prises de pression pariétales
reliées à des capteurs de pression et de calculer les efforts globaux par intégration de ces pressions sur la
surface de la maquette. En plus d’obtenir des informations sur l’allure des champs de pression (voir
figure 3.2), on peut accéder à des valeurs de coefficients « instantanés », la bande passante des capteurs
de pression étant plus large que celle des balances. Cette méthode présente toutefois les inconvénients
d’être plus lourde à mettre en place et moins précise en ce qui concerne la mesure de coefficients
moyens.
Figure 3.2 – Champs de pressions moyens et fluctuants autour du tablier trapézoïdal

pour les incidences –3°, 0°, +3°
55
Figure 3.3 – Capteurs de pression à l’intérieur d’une maquette sectionnelle
Dans tous les cas les maquettes sont fixées sur le "banc dynamique" de la soufflerie atmosphérique du
CSTB (figure 3.12). Ce banc est composé de deux parois verticales épaisses de 15 cm appelées
"flasques". Le profilage en amont permet de minimiser l'épaisseur de la couche limite sur les parois
internes et garantit ainsi un écoulement bidimensionnel. C'est à l'intérieur de ces flasques que sont placés
les balances et autres dispositifs métrologiques.
Nous avons mis en place un montage permettant de comparer les résultats obtenus par les deux
méthodes. Nos deux tabliers, équipés de 64 prises de pressions, ont été fixés à des balances
dynamométriques et soumis à un flux laminaire. Les mesures ont donc été effectuées simultanément avec
les balances et les capteurs de pression.
Les 2*32 signaux analogiques de pression sont multiplexés à une cadence de 20 kHz et envoyés en
rafales à une carte d’acquisition 16 bits à une cadence de 200 Hz . Les signaux analogiques issus du tube
de Pitot et des capteurs de déplacements, filtrés électroniquement à 100 Hz, sont également multiplexés
et placés à la suite des mesures de pressions avant d’être transmis à la carte d’acquisition, puis à une
station de travail. Les données numériques sont « synchronisées » pendant la phase de traitement par des
interpolations de type « Stirling ».
56
Figure 3.4 – Profils des maquettes sectionnelles des tabliers trapézoïdal et rectangulaire de ratio 8
0 ,1 2 0 ,4 0 ,2
0 ,0 8 0 ,0 0 ,1
-6 -3 0 3 6
C x int . C y int . C m int .

0 ,0 4 - 0 ,4 0 ,0
C x b a l. C y b a l. C m b a l.
-6 -3 0 3 6
0 ,0 0 - 0 ,8 - 0 ,1
-6 -3 0 3 6
- 0 ,0 4 - 1 ,2 - 0 ,2
Figure 3.5 – Coefficients aérodynamiques du tablier trapézoïdal mesurés par pressions et balances
0 ,1 6 0 ,8 0 ,0 3
0 ,0 2
0 ,1 5
0 ,4
0 ,0 1
C x int .
0 ,1 4 C y int . C m int .
C x b a l. 0 ,0 0 ,0 0
C y b a l. C m b a l.
-6 -3 0 3 6 -6 -3 0 3 6
0 ,1 3 - 0 ,0 1
- 0 ,4
- 0 ,0 2
0 ,1 2
-6 -3 0 3 6 - 0 ,8 - 0 ,0 3
Figure 3.6 – Coefficients aérodynamiques du tablier rectangulaire mesurés par pressions et balances
La comparaison des résultats obtenus en terme de coefficients stationnaires s'avère globalement

satisfaisante. Seul le coefficient de moment du tablier rectangulaire de ratio 8 posait un petit problème :
ses valeurs, très faibles en valeur absolue, sont de signes opposés selon que la mesure est globale ou
médiane.
57
Les deux mesures semblant « justes », nous avons décidé d’approfondir la question en étudiant
l’évolution des courbes de moment de tabliers rectangulaires autour du ratio 8. Nous espérions que les
tendances observées en faisant évoluer le ratio dans un sens ou dans l’autre nous permettraient de mieux
qualifier les propriétés aérodynamiques naturelles du ratio 8 et identifier les effets 3D.
Pour cela, nous avons placé sur les balances du banc dynamique une âme rectangulaire en contreplaqué
longue de 1,6 m sur laquelle nous avons pu fixer d’autres planches d’épaisseurs variables. Pour chaque
ratio, nous avons alors mesuré les coefficients de moment pour 3 vitesses de vent et 7 incidences.
Figure 3.7 – Etude paramétrique de tabliers rectangulaires
L’analyse des résultats, présentés figure 3.8, nous a amené à formuler l’explication suivante :
Deux différentes structures d’écoulement se développent le long d’une maquette sectionnelle

rectangulaire de ratio 8 : une première au centre, générant les coefficients de moment vus par les capteurs
de pression, une seconde, à chaque extrémité, générant des coefficients de moment proches de ceux
mesurés par les balances sur les ratios supérieurs. La première structure serait d’autant plus favorisée que
le profil est épais tandis que la seconde se développerait davantage quand le profil est élancé.
Dans le cas de la maquette sectionnelle de 1 m de long, les structures de bords sont prépondérantes et
constituent l’essentiel des forces mesurées par les balances. Dans le cas de la maquette sectionnelle de
1,6 m, la structure centrale se développe plus largement et prend le pas sur les structures latérales, ce qui
se traduit par une courbe intermédiaire entre les deux courbes du tablier de 1 m de long.
58
Tablier rectangulaire ratio 12 Tablier rectangulaire ratio 8
Cm
Cm
0,10
0,10
0,08
0,08
0,06
0,06
0,04
0,04
0,02 0,02
0,00 0,00 i
i
-6 -4 -2 0 2 4 6 -6 -4 -2 0 2 4 6
-0,02 -0,02
-0,04 -0,04
U=6m/s U=6m/s
-0,06 U=10m/s -0,06
U=10m/s
-0,08 U=15m/s -0,08 U=15m/s
-0,10 -0,10
Tablier rectangulaire ratio 9,5 Tablier rectangulaire ratio 7
Cm
Cm
0,10
0,10
0,08
0,08 U=6m/s
0,06 U=10m/s
0,06
U=15m/s
0,04 0,04
0,02 0,02
0,00 i 0,00 i
-6 -4 -2 0 2 4 6 -6 -4 -2 0 2 4 6
-0,02 -0,02
-0,04 -0,04
U=6m/s
-0,06 U=10m/s -0,06
-0,08 U=15m/s
-0,08
-0,10 -0,10
Tablier rectangulaire ratio 8,5 Tablier rectangulaire ratio 5,5
Cm
Cm
0,10
0,10
0,08 U=6m/s
0,08
0,06 U=10m/s
0,06
U=15m/s
0,04 0,04
0,02 0,02
0,00 i 0,00 i
-6 -4 -2 0 2 4 6 -6 -4 -2 0 2 4 6
-0,02 -0,02
-0,04 -0,04
U=6m/s
-0,06 U=10m/s -0,06
-0,08 U=15m/s
-0,08
-0,10 -0,10
Figure 3.8 – Evolution des coefficients de moment en fonction de l’incidence du vent pour des tabliers
rectangulaires de ratios variables
59
Les hypothèses formulées nous amènent à conclure que les résultats obtenus à partir des prises de
pressions sur le tablier de 1 m de long seraient bien représentatives de phénomènes « centraux » ; il ne
s’agirait donc effectivement pas d’erreurs de mesures. Par contre, la question qui se pose encore est la
suivante : les structures latérales proviennent-elles d’effets de bords ou se développent elles
naturellement autour d’une structure centrale ? En d’autres termes, peut on considérer la structure
centrale comme représentative de ce qui se passerait sur un tablier infiniment long ? Ou bien peut-il se
développer, par exemple, un phénomène « périodique » le long du tablier avec alternance de structures
« centrales » et « latérales »?
Pour répondre à cette question, on pourrait réaliser une étude paramétrique dont la variable serait la
longueur du tablier ; des couronnes de prises de pressions seraient alors disposées à intervalles réguliers.
Cependant, la faible valeur des coefficients de moment nous autorise à ne pas approfondir d’avantage un
problème dont les répercussions sur le calcul dynamique seront vraisemblablement faibles.
Par contre, on peut se demander si les phénomènes aéroélastiques seront ou non sujets à des problèmes
similaires. Faute de moyens pour répondre à cette interrogation, nous avons préféré, dans le cas du tablier
rectangulaire de ratio 8, ne pas utiliser les coefficients aéroélastiques issus des oscillations forcées dans
nos calculs de réponse dynamique mais plutôt ceux provenant de la méthode des oscillations libres,
garantissant un meilleur « moyennage » des phénomènes tout au long du tablier.
3.2 Mesure de l'admittance aérodynamique
Comme nous l'avons dit dans le chapitre 2, les coefficients aérodynamiques stationnaires perdent de leur
pertinence à mesure que l'on s'approche des sollicitations de hautes fréquences. La fonction d'admittance
aérodynamique (f.a.a.) permet de palier à ce problème. Dans l'absolu, il existe autant de fonctions
d'admittance aérodynamique que de couples effort / composante turbulente. Expérimentalement, il est
difficile de dissocier les effets des différentes composantes du vent. On se contente donc de mesurer trois
f.a.a. associées aux forces de portance, traînée et moment.
Ces trois f.a.a. peuvent être appliquées directement au calcul spectral. Leur introduction dans la méthode
temporelle nécessite, en principe, l’emploi de produits de convolution. Ces derniers entraînant des
surcoûts importants en temps de calcul, on préfèrera raisonner de manière fréquentielle dans la base
modale.
L'estimation de la f.a.a. consiste à diviser un spectre de forces "réaliste" c’est-à-dire directement mesurés
sur maquettes sectionnelles à l'aide de balances ou de prises de pressions par un spectre de forces "quasi-
60
statique" i.e. calculé d'après les spectres des deux composantes turbulentes du vent. Ces derniers sont
calculés à partir d'acquisitions effectuées à l'aide d'une sonde à fils chauds deux composantes.
Les résultats obtenus pour nos deux tabliers sont présentés figures 3.9 et 3.10. Deux f.a.a. ont été
estimées, correspondant aux forces de portance et de moment. Elles sont comparées à deux formulations
issues de la littérature.
Dans le cas du tablier rectangulaire, la f.a.a. associée au moment est particulièrement atypique puisque
qu’elle semble environ trois fois « supérieure à la normale ». Cela n’est en fait pas si surprenant compte
tenu des problèmes rencontrés dans la détermination des coefficients de moment : ce que l’on appelle ici
f.a.a. est vraisemblablement le produit entre la véritable f.a.a. et un facteur correctif intégrant l’erreur
d’estimation des coefficients de moment.
1,0
0 1 2 Admittance portance 3
Admittance moment
Vickery
Liepmann
0,1
0,0
Fréquence réduite fB/U
Figure 3.9 – Admittances aérodynamiques du tablier trapézoïdal
10,0 Admittance portance

Admittance moment
Vickery
Liepmann
1,0
0 1 2 3
0,1
0,0
Fréquence réduite fB/U
Figure 3.10 – Admittances aérodynamiques du tablier rectangulaire
61
3.3 Mesure des coefficients aéroélastiques par la technique des

oscillations libres
3.3.1 Bases théoriques
Les phénomènes d'instabilité induits par le vent sur un pont à haubans concernent généralement les
modes verticaux et les modes de torsion du tablier de l'ouvrage. On a vu que ces phénomènes pouvaient
être caractérisés par les coefficients aéroélastiques. La méthode expérimentale classique pour mesurer
ces coefficients, est la méthode des oscillations libres développée par R.H. Scanlan (1974). Elle consiste
à réaliser un tronçon représentatif d'une partie du tablier, que l’on suspend en soufflerie par un ensemble
de ressorts afin de représenter un mode de torsion et un mode de flexion verticale.
Figure 3.11 – Schéma de principe de la méthode des oscillations libres
Dans la mesure où ces essais sur tronçons ont pour but d'exprimer uniquement les caractéristiques
aérodynamiques du tablier sous la forme de coefficients sans dimensions, le modèle testé n'a plus besoin
d'être dynamiquement semblable au pont réel.
62
Toutefois, il est préférable que les caractéristiques dynamiques de la maquette restent proches de celles
imposées par la similitude. Cela permet la mesure en soufflerie d'une vitesse critique de flottement
directement transposable en vitesse réelle via les relations de similitude. En outre, les éventuels
problèmes de flottement couplé flexion/torsion sont bien mis en évidence.
Dans le cas où l’on ne considère que deux degrés de liberté (flexion h(t) et torsion (t)), les équations du
mouvement du tronçon d'un tablier de pont peuvent s'écrire avec les définitions de Scanlan :
1 2  h& B& 2 K2 * 
m (&h&  2 h  h h&  2h h ) = U B KH1*  KH *2  K H *3   H 4   L h (t ) (3.1)
2  U U B 
1 2 2  h& B& 2 K2 * 
j0 ( 
&&  2   &  2  ) = U B  KA1*  KA *2  K A *3   A 4   M  (t ) (3.2)
2  U U B 
m, j0 : masse et moment d'inertie par unité de longueur du modèle

 h ,  : amortissements structurels rapportés au critique pour chaque degré de liberté h et 
h,  : pulsations naturelles du mouvement à la vitesse U = 0
B
K= : fréquence réduite
U
Les forces L h ( t ), M  ( t ) sont induites par un écoulement turbulent. Celles-ci sont donc supposées nulles
lors d'essais en écoulement uniforme. Les équations prennent alors la forme :
&h&  2  h&   2 h  H h&  H h  H &  H  (3.3)

h h h 1 4 2 3
&&  2   &  2   A1 h&  A 4 h  A 2 &  A 3 

 (3.4)
Signalons que, dans un souci de simplification des écritures, la définition des termes Hi et A i (i=1,2,3)
est différente de celle employée dans le chapitre 2. Ici, on écrit :
B 2 H1* B 2  2 H *4 B3 H *2 B3  2 H *3

H1  ; H4  ; H2  ; H3  (3.5)
2m 2m 2m 2m
B3 A1* B3 2 A *4 B 4 A *2 B 4 2 A *3

A1  ; A4  ; A2  ; A3  (3.6)
2j 2j 2j 2j
63
Les coefficients aéroélastiques sont supposés indépendants de la nature du mouvement ; on peut donc
imposer successivement un mouvement de flexion verticale et un mouvement de torsion et ainsi obtenir
deux équations découplées :
&h&  2  h&   2 h  H h&  H h (3.7)

h h h 1 4

&&  2   &  2   A 2 &  A 3  (3.8)
Lorsque le tablier est en oscillations libres, l’équation du mouvement peut alors s'écrire sous la forme
suivante:
h  h 0 e 11t sin  1 1  1  t  1 

2
(3.9)
  
   0 e  22 t sin   2 1   2  t   2 

2
(3.10)
  
1,2 représente l’amortissement total du système "rapporté au critique" (amortissement structurel +

amortissement aérodynamique) et 1,2 la pulsation du mouvement. Les pulsations et amortissements du
système sont mesurés directement à partir des signaux temporels, période par période (fig.3.12) : les
pulsations par identification des passages à zéro, les amortissements par la formule classique :
   1 T  ln  X p max X p 1max  (3.11)
où T, Xpmax, Xp+1 max, représentent la pseudo-période et deux extremums successifs du signal.
D’autres méthodes plus sophistiquées comme la R.D.D. (Random Decrement Technique), dans le
domaine temporel, ou le Maximum de Vraisemblance, dans le domaine spectral, peuvent également être
utilisées mais présentent l’inconvénient de ne pas fournir d’informations sur l’évolution des deux
variables au cours du temps. En outre, elles sont plus appropriées en vent turbulent.
Dans le cas où i > 0, i.e. l’amortissement total est positif, les signaux de déplacements du tablier ont
l’allure présentée figure 3.12. Expérimentalement, on constate que l’amortissement i varie en fonction
de l’amplitude des oscillations. Ceci est essentiellement dû au fait que l’amortissement structurel, de
nature visqueuse pour les mouvements de grande amplitude, voit sa composante solide augmenter à
mesure que la taille des oscillations diminue.
64
Figure 3.12 – Deux excitations et lâchers : signal de déplacement ( > 0)
Supposons le mouvement de torsion bloqué ; la lecture de l'enregistrement du mouvement de flexion

permet, connaissant l’amortissement structurel h et la pulsation naturelle du système h, la détermination
de 1 et de 1, d'où celle de H1 et H4 par l'intermédiaire des équations :
H1  2(h h  11 ) (3.12)
H 4   2h  12  1 1 
2
(3.13)
De même,
A 2  2(   2  2 ) (3.14)
A 3  2  22  2  2 

2
(3.15)
Les coefficients H1 et A2 traduisent l'apport de l'amortissement aérodynamique à l'amortissement total,

respectivement en flexion et en torsion, tandis que les coefficients H4 et A3 traduisent le glissement en
fréquence induit par l'écoulement.
65
Lorsque H1, H4, A2 et A3 sont connus, le tablier est excité simultanément en flexion et en torsion à des
fréquences d’oscillations aussi proches que possible ; le système d'équations issu de (3.4), peut alors être
utilisé pour trouver les coefficients de couplage H 2 , H 3 , A 1 et A 4 .
Lorsque l'amortissement aérodynamique est négatif (coefficients positifs, H1 pour la flexion, A 2 pour la
torsion), celui-ci est soustrait à l'amortissement structurel. S'il devient supérieur en valeur absolue à
l'amortissement structurel, il y a alors flottement à un degré de liberté. Dans la très grande majorité des
cas, c'est ce flottement à un degré de liberté, notamment en torsion, qui a tendance à apparaître sur les
ponts. Le flottement couplé flexion/torsion ne peut apparaître que si les amortissements totaux en flexion
et en torsion restent positifs et les fréquences de vibrations en flexion et en torsion sont proches ; en
général, seuls les tabliers de pont très profilés, sont susceptibles de flottement couplé.
3.3.2 Application pratique
Nous présentons les résultats obtenus pour les deux maquettes sectionnelles rectangulaire et trapézoïdale.
Ces maquettes sont positionnées entre les flasques du banc dynamique.
Dans le cas des essais en flexion verticale, un câble de retenue permet de maintenir la maquette
horizontalement ; pour les essais en torsion, le mouvement est guidé par des pointaux. L’amortissement
structurel de base est ainsi minimal. Des dispositifs amortisseurs supplémentaires permettent de régler
l'amortissement structurel rapporté au critique des mouvements de flexion et torsion entre 0,1% et 1% .
Un jeu minimal de l'ordre de 3 mm, sépare l'extrémité de la maquette de la paroi du banc dynamique. Une
excitation initiale suivie d'un lâcher sont réalisés par un excitateur électromagnétique (figure 3.13) piloté
depuis le caisson de mesures. Cette excitation est réalisée à la fréquence de vibration en flexion verticale
avec une traction au centre du tablier, puis à la fréquence de torsion avec une traction latérale. Les
fréquences d’oscillations peuvent être modulées par intervention sur la raideur des ressorts, leur
écartement et/ou l’ajout de masses sur les bras situés à l’intérieur des flasques.
Les déplacements de la maquette sont mesurés par des capteurs optiques à triangulation placés dans les
flasques du banc dynamique. Les signaux analogiques sont discrétisés à une cadence de 500 Hz et filtrés
électroniquement à 200 Hz. La vitesse du vent est mesurée par un anémomètre à hélices et un tube de
Pitot situé en amont du banc dynamique. Un recalage a préalablement été effectué entre la pression
dynamique amont mesurée et la pression dynamique au niveau du bord d’attaque de tablier. Chaque
acquisition comporte 8096 points, ce qui correspond à une durée de 16 secondes.
66
Figure 3.13 –Banc dynamique et excitateur électromagnétique
Figure 3.14 – Vue de l’intérieur d’un flasque du banc dynamique
67
Les courbes d’évolution des coefficients aéroélastiques en fonction de la vitesse réduite (U/fB) peuvent
théoriquement être obtenues en fixant une fréquence de vibration et en ne faisant varier que la vitesse du
vent. Nous avons préféré faire varier ces deux paramètres afin d’estimer l’impact éventuel du nombre de
Reynolds sur la valeur des coefficients aéroélastiques. Nous avons donc tracé, pour chaque fréquence de
vibration, les courbes d’évolution des coefficients de Scanlan H *i , A *i en fonction de la vitesse réduite.
Précisons que cette vitesse réduite est définie à partir de la fréquence d’oscillation du système avec effets
aéroélastiques ; cela est en effet souvent source d’ambiguïtés.
L’analyse des résultats obtenus ne permet pas réellement de conclure sur l’importance de l’influence du
nombre de Reynolds. En effet, on n’observe pas d’évolution claire des courbes en fonction de la
fréquence propre du système. On observe par contre une forte dispersion des points. Celle-ci est
essentiellement imputable aux incertitudes de mesures, aux imperfections du système mécanique
(variations des raideurs et amortissements mécaniques en fonction de l’amplitude) et, dans le cas des
montages en flexion, à un comportement dynamique assez chaotique des maquettes pour les grandes
vitesses de vent.
Ce dernier point s’explique par le fait qu’à partir d’une certaine vitesse réduite (environ cinq pour le
profil trapézoïdal) l’écoulement de sillage autour du tablier (sa signature aérodynamique) devient
fortement turbulent. Cette turbulence excite non seulement le mode de flexion qui nous intéresse mais
également un mode de torsion autour de l’axe x, de fréquence légèrement supérieure. Les figures 3.15-18
montrent les signaux temporels (et transformées par ondelettes) obtenus pour le tablier rectangulaire en
flexion et en torsion pour deux vitesses de vent .
Une autre source d’erreurs, intrinsèque à la méthode des oscillations libres, est l’évolution naturelle de
l’incidence moyenne du tablier en fonction de la vitesse moyenne du vent. Ainsi, chaque série de points
présentée figure 3.19-26 correspond non seulement à une vitesse moyenne différente mais également à
une incidence de vent différente. En configuration de flexion, ce phénomène peut être minimisé en
écartant au maximum les lignes de ressort et donc en augmentant la raideur en torsion. En outre, les
tabliers symétriques sont a priori peu sujets à ce genre de problème. Par contre, dans le cas de tabliers
asymétriques en torsion libre, le problème est réel.
68
Figure 3.15 – Oscillations verticales du tablier rectangulaire de ratio 8 (U=5m/s)
Figure 3.16– Oscillations verticales du tablier rectangulaire de ratio 8 (U=15m/s)
69
Figure 3.17 – Oscillations en torsion du tablier rectangulaire de ratio 8 (U=5m/s)
Figure 3.18 – Oscillations en torsion du tablier rectangulaire de ratio 8 (U=15m/s)
70
0
0 2 4 6 8 10 12
-1
-2
4Hz
-3 5Hz
H1*
7Hz
-4 9Hz
10Hz
-5
-6
-7
Vitesse réduite
Figure 3.19 – Coefficient H1 du tablier trapézoïdal en fonction de la vitesse réduite

pour plusieurs fréquences d’oscillations
0
0 2 4 6 8 10 12
-2
-4
4Hz
-6
5Hz
H1*
-8 6Hz
8Hz
-10 9Hz
-12 10Hz
-14
-16
Vitesse réduite
Figure 3.20 – Coefficient H1 du tablier rectangulaire en fonction de la vitesse réduite

71
0
0 2 4 6 8 10 12
-1
-2
4Hz
H4*
5Hz
-3
7Hz
9Hz
-4
10Hz
-5
-6
Vitesse réduite
Figure 3.21 – Coefficient H4 du tablier trapézoïdal en fonction de la vitesse réduite

0
0 2 4 6 8 10 12
-2
4Hz
-4 5Hz
H4*
6Hz
8Hz
-6
9Hz
10Hz
-8
-10
Vitesse réduite
Figure 3.22 – Coefficient H4 du tablier rectangulaire en fonction de la vitesse réduite

72
0.1
0.0
0 2 4 6 8 10
4Hz
-0.1 5Hz
A2*
6Hz
7Hz
-0.2
8Hz
9Hz
-0.3 10Hz
-0.4
Vitesse réduite
Figure 3.23 – Coefficient A2 du tablier trapézoïdal en fonction de la vitesse réduite

0.1
0.0
0 1 2 3 4 5 6 7 8
-0.1 6Hz
7Hz
A2*
8Hz
-0.2 10Hz
11Hz
-0.3 12Hz
-0.4
Vitesse réduite
Figure 3.24 – Coefficient A2 du tablier rectangulaire en fonction de la vitesse réduite

73
4Hz
3
5Hz
6Hz
7Hz
A3*
2
8Hz
9Hz
10Hz
1
0
0 2 4 6 8 10
Vitesse réduite
Figure 3.25 – Coefficient A3 du tablier trapézoïdal en fonction de la vitesse réduite

3
6Hz
7Hz
A3*
2 8Hz
10Hz
11Hz
1 12Hz
0
0 1 2 3 4 5 6 7 8
Vitesse réduite
Figure 3.26 – Coefficient A3 du tablier rectangulaire en fonction de la vitesse réduite

74
3.4 Méthode des oscillations forcées
Contrairement à la méthode des oscillations libres, les forces aéroélastiques sont directement mesurées
par des capteurs de pression et non déduites des équations de la dynamique. Cela permet de s’affranchir
des caractéristiques mécaniques du système, qui occasionnaient des incertitudes supplémentaires. Un
traitement conjoint des signaux de forces et de déplacements de la maquette sectionnelle permet de
calculer les coefficients aéroélastiques des équations de Scanlan.
3.4.1 Instrumentation de la maquette
Les maquettes utilisées sont celles des oscillations libres. Elle sont cependant équipées d’une couronne
de 64 prises de pression reliées par des tubes de vinyle à deux boîtiers composés chacun de 32
membranes de silicone piézoélectriques. Chaque tube est équipé d’un « restricteur » chargé de filtrer les
phénomènes de résonance pneumatique. Pour minimiser les effets d’accélération, les boitiers sont fixés
« couchés » à l’intérieur des maquettes sectionnelles. On verra que, malgré cette précaution, les capteurs
sont assez sensibles à l’accélération.
La maquette sectionnelle est excitée sinusoïdalement, via un système articulé, par un moteur électrique à
courant continu placé dans un des flasques du banc dynamique (voir figure 3.28). Trois types de
mouvements peuvent être générés: flexion verticale, torsion, et les deux couplés à la même fréquence.
Différentes amplitudes d’oscillations peuvent être imposées dépendant de la position du bras sur le
disque en sortie du moteur ; on peut ainsi en étudier l’influence sur l’allure du champ de pression. De
plus, on peut faire varier la fréquence d’oscillation par simple réglage de l’alimentation du moteur, et
ainsi balayer une large gamme de fréquences réduites pour une vitesse de vent donnée. Les déplacements
sont mesurés par deux capteurs à triangulation lasers placés à l’intérieur des flasques du banc dynamique.
Les 2*32 signaux analogiques de pression sont multiplexés à une cadence de 20 kHz et envoyés en
rafales à une carte d’acquisition 16 bits à une cadence de 200 Hz. Les signaux analogiques issus du tube
de Pitot et des capteurs de déplacements, filtrés électroniquement à 100 Hz, sont également multiplexés
et placés à la suite des mesures de pressions avant d’être transmis à la carte d’acquisition, puis à une
station de travail. Les données numériques sont « synchronisées » pendant la phase de traitement par des
interpolations de type « Stirling ».
75
Figure 3.27– Partie haute du système articulé : bras de transmission, capteurs de déplacements…
Figure 3.28 – Partie basse du système articulé : moteur électrique, excentrique…
76
3.4.2 Cas des oscillations pures
Contrairement à la méthode des oscillations libres, la méthode des oscillations forcées permet de mesurer
tous les coefficients aéroélastiques, y compris les termes couplés, à partir d’oscillations pures. Si l’on
prend le cas d’un mouvement de flexion verticale, l’équations (3.1) peut être réécrite sous la forme:
F sin(t  )  C H14 H 1* H cos t   C H14 H *4 H sin t  (3.16)
où F et H représentent respectivement les amplitudes des forces de portance et des déplacements
verticaux et C H14  1 B 2  2 . On en déduit:

2
H 1*  F C H14 H  sin  (3.17)
H *4  F C H14 H  cos  (3.18)
De la même façon, considérant l’équation (3.2), et prenant C A14  1 B 3  2 , on peut écrire:

2
A 1*  F C A14 H  sin  (3.19)
A *4  F C A14 H  cos  (3.20)
Les quatre autres coefficients sont calculés de la même manière à partir des oscillations en torsion. H est
remplacé par A, amplitude des oscillations en torsion. On introduit alors les paramètres
C H 23  1 2 B 3  2 et C A 23  1 2 B 4  2 , respectivement pour les coefficients couplés et non couplés.
Les équations (3.17), (3.18), (3.19) et (3.20) montrent que les calculs des coefficients aéroélastiques
peuvent se résumer à des calculs de rapports d’amplitudes, de fréquences d’oscillations, et de
déphasages.
On doit cependant penser à retirer des signaux de forces les composantes dues à l’accélération. Pour cela,
on effectue des mesures sans vent pour les différentes amplitudes et fréquences d’oscillations. On déduit
de l’analyse comparative des signaux de déplacements et de "forces", une fonction de transfert. En
principe, cette fonction de transfert est un simple scalaire de la forme C. puisque les signaux de
déplacements et d’accélérations sont en phase. Par la suite, on retranchera simplement au signal de forces
le signal de déplacement multiplié par le scalaire en question. La figure 3.28 montre des exemples de
signaux de déplacements, accélération et forces de portance + accélérations pour une excitation à 5 Hz et
un vent à 12 m/s.
77
Figure 3.29 – Signaux de déplacements, accélération et portance pour une excitation à 5Hz
78
3.4.3 Mesure des rapports d’amplitudes, fréquences d’oscillations et déphasages
Avant de calculer les rapports d’amplitudes, les signaux de forces et de déplacements sont
numériquement filtrés passe-bas. On a vérifié que ce filtre n’occasionnait pas de déphasage parasite.
Puis, comme l’amplitude des oscillations varie légèrement avec le temps, à cause des imperfections du
système mécanique, les rapports sont calculés période par période puis moyennés. Les fréquences et
déphasages entre signaux sont estimés à partir du calcul de la fonction de densité interspectrale entre le
signal des forces et celui des déplacements. Afin d’améliorer la précision des résultats, un spectre
représentant la réponse d’un oscillateur linéaire non amorti à un degré de liberté est ajusté aux valeurs
brutes par une méthode d’optimisation basée sur les moindres carrés.
1000
Spectre de puissance
100
Données brutes
Modèle 1d.d.l.
10
1
3 4 5 6
0,1
Fréquence
Figure 3.30 – Détermination de la fréquence d’un signal par ajustement

d’un modèle d’oscillateur linéaire à 1 d.d.l. par une méthode moindres carrés
On peut quand même signaler que plus la fréquence d’oscillations augmente, moins les déphasages sont
estimés avec précision ; tout d’abord parce que le nombre de points par période diminue, ensuite parce
que les déphasages électroniques et mécaniques introduits par les différents maillons du système
d’acquisition sont difficiles à caractériser.
3.4.4 Cas du mouvement couplé
Au moment du traitement des essais en flexion forcée, nous avons noté l’existence de légères oscillations
en torsion, à la même fréquence que les oscillations en flexion. Les déplacements induits par ce
mouvement de torsion en rive étaient de l’ordre de 10% des mouvements générés par le flexion verticale.
On a donc dû traiter ce cas en tant que mouvement couplé en introduisant les coefficients mesurés
pendant les oscillations en torsion pure. Notons au passage qu’il serait tout à fait envisageable et
intéressant, après avoir amélioré le système mécanique du montage en flexion, de comparer les résultats
issus des oscillations simples et ceux issus des oscillations couplées.
79
Donc, dans le cas où les oscillations en flexion sont couplées à des oscillations en torsion, à la même
fréquence, l’équation (3.1) est réécrite sous la forme :
   
F sin t    C1 A H *2 cos t  H *3 sin t  C 2 H H 1* cost     H *4 sin t   (3.21)
où C1  1 2 B 3  2 et C 2  1 2 B 2  2 .
que l’on peut condenser en:
F sin t    F1 sin t  1   F2 sin t   2    (3.22)
Supposant que nous connaissons les valeurs de H *2 et H *3 et que nous cherchons à calculer H 1* et H *4 , la
procédure de calcul est la suivante :
 F
 Calcul des rapports d’amplitudes des signaux : et
F H
F1 F1 A F 2 2
   avec 1  C1 H *2  H *3
F A F A
 1  a tan H *2 H *3 
2
F2 F F
  1  12  2 1 cos  1 
F F F
 F 
 sin   1 sin 1 
  2  a tan  F 
 F
cos   1 cos 1 
 F 
F2 F 1
 H 1*    sin  2
F H C2
F2 F 1
 H *4    cos  2
F H C2
Les autres paires de coefficients sont calculées de la même manière. Deux paires de coefficients doivent
être préalablement estimées à partir d’un montage simple afin d’en calculer une autre issue d’un montage
couplé.
80
3.4.5 Résultats
Trois maquettes sectionnelles ont été testées dans le cadre d’un partenariat avec le LCPC, l’INRIA et
l’ENPC : deux profils rectangulaires de ratios 4 et 8 et le profil trapézoïdal (figures 3.31 et 3.32). Seuls
les résultats du tablier trapézoïdal et du tablier rectangulaire de ratio 8 sont présentés ici, le tablier de
ratio 4 n’étant pas étudié par la suite.
Les mesures ont été effectuées en vent uniforme pour deux vitesses de vent moyennes : 13 m/s et 18 m/s.
Trois amplitudes d’oscillations ont été testées en torsion : 8°, 4°, 2°, deux en flexion : 13 mm et 30 mm
de crête à crête. La gamme de fréquence d’oscillations était de 2 à 7 Hz. Les résultats sont présentés
figures 3.33 à 3.36.
Figure 3.31 – Implantation des couronnes de pressions du tablier trapézoïdal
Figure 3.32 –Coupe sur une des couronnes de pressions du tablier trapézoïdal
81
0
0 5 10 15 20 25
-5
-10
A mp l. 13mm
H1 *
-15
A mp l. 30mm
-20
-25
-30
Vites s e réduite
40
35
30
25
A mpl. 1°
20
H2*
A mpl. 2°
15
10
0
0 5 10 15 20 25 30
-5
Vites s e réduite.
0
0 5 10 15 20 25 30
-20
-40
A mpl. 1°
H3*
-60
A mpl. 2°
-80
-100
-120
Vites s e réduite
0
0 5 10 15 20 25
-2
H4 *
-4
A mpl. 13mm
-6 A mpl. 30mm
-8
-10
Vites s e réduite
Figure 3.33 – Coefficients aéroélastiques de portance du tablier trapézoïdal obtenus par la méthode des
oscillations forcées (en couleur les résultats issus des oscillations libres)
82
A mpl. 13mm
3
A mpl. 30mm
A1*
2
0
0 5 10 15 20 25
Vites s e réduite
5
A mp l. 1°
A2 *
4 A mp l. 2°
0
0 5 10 15 20 25 30
Vites s e réduite
30
25
20
A mp l. 1°
A3 *
15 A mp l. 2°
10
0
0 5 10 15 20 25 30
Vites s e réduite
0,9
0,8
0,7
0,6
0,5 A mp l. 13mm
A4*
0,4 A mp l. 30mm
0,3
0,2
0,1
0
0 5 10 15 20 25
Vites s e réduite
Figure 3.34 – Coefficients aéroélastiques de moment du tablier trapézoïdal obtenus par la méthode des
oscillations forcées (en couleur les résultats obtenus par la méthode des oscillations libres)
83
0
0 2 4 6 8 10 12 14 16
-5
H1*
-1 0 Am p l.1 3 m m
Am p l.3 0 m m
-1 5
-2 0
Vi te s s e ré du i te
0
0 2 4 6 8 10 12 14 16
-5
-1 0
Am p l.1 °
H 2*
Am p l.2 °
-1 5
-2 0
-2 5
40
30
H 3*
20 Am p l.1 °
Am p l.2 °
10
0
0 2 4 6 8 10 12 14 16
-1 0 2 4 6 8 10 12 14 16
-2
-3
-4
H 4*
-5
-6
Am p l.1 3 m m
-7 Am p l.3 0 m m
-8
-9
-1 0
Figure 3.35 – Coefficients aéroélastiques de portance du tablier rectangulaire obtenus par la méthode
des oscillations forcées (en couleur les résultats obtenus par la méthode des oscillations libres)
84
0
0 2 4 6 8 10 12 14 16
Am p l.1 3 m m
A1 *
-1 Am p l.3 0 m m
-2
3
A2 *
Am p l.1 °
2 Am p l.2 °
0
0 2 4 6 8 10 12 14 16
3
A3 *
Am p l.1 °
2 Am p l.2 °
0
0 2 4 6 8 10 12 14 16
2
A4 *
Am p l.1 3 m m
1 Am p l.3 0 m m
0
0 2 4 6 8 10 12 14 16
Figure 3.36 – Coefficients aéroélastiques de moment du tablier rectangulaire obtenus par la méthode
des oscillations forcées (en couleur les résultats obtenus par la méthode des oscillations libres)
85
Si les coefficients aéroélastiques du tablier trapézoïdal semblent peu dépendants de l’amplitude des
oscillations, cela n’est pas le cas pour le tablier rectangulaire. On peut donc en conclure que, pour ce
dernier, les formulations de Scanlan perdent de leur pertinence pour des mouvements de grande
amplitude. Au contraire, le tablier trapézoïdal semble plus influencé par le nombre de Reynolds que le
tablier rectangulaire; cela est particulièrement net pour les coefficients H1, H4, A1, A4 et se traduit par
deux séries de points légèrement décalées, correspondant aux deux vitesses de vent imposées.
Dans le cas des oscillations verticales, on observe une dispersion assez importante des résultats, en
particulier du côté des coefficients couplés. Cela est partiellement imputable au mouvement de torsion
parasite, dont l’évaluation s’est avérée peu précise.
La comparaison des résultats obtenus par les méthodes libre et forcée ne nous permet pas de tirer de
conclusion bien franche. En effet, si les coefficients H1 et A3 sont peu sensibles à la méthode employée,
l’estimation des coefficients H4 du tablier rectangulaire et A2 du tablier trapézoïdal est franchement
différente selon la technique utilisée. La surestimation du H4 en oscillations libres peut être partiellement
attribuée à l’influence de l’amplitude des oscillations, plus faibles qu’en forcé. Le coefficient A2, lui,
illustre bien quelques points faibles de la méthode des oscillations forcées : précision moyenne dans les
hautes fréquences et, surtout, mauvaise appréhension des phénomènes physiques . En effet, une erreur
absolue modérée peut induire un changement de signe sans qu’on puisse se rendre compte visuellement
du caractère stable ou instable du profil.
N’oublions pas que les phénomènes observés par les deux méthodes ne sont pas totalement identiques :
dans le cas des oscillations libres, la structure « s’adapte » sans arrêt à l’écoulement ; dans le cas des
oscillations forcées c’est plutôt l’écoulement qui s’adapte aux mouvements de la structure. Notamment,
le phénomène de « glissement en fréquence » n’est absolument pas reproduit.
Par ailleurs, les incidences du tablier ne sont pas les mêmes dans les deux méthodes. La méthode des
oscillations forcées assure une incidence moyenne constante bien connue, en l’occurrence, ici, 0°. Avec
la méthode des oscillations libres, l’écoulement moyen impose une incidence moyenne différente, a
priori plus instable. Cela peut expliquer que le tablier trapézoïdal s’avère plus instable d’après cette
méthode (incursion du coefficient A2 vers les ordonnées positives).
Enfin, il convient de ne pas opposer en bloc deux méthodes qui sont aussi complémentaires. La méthode
des oscillations libres est en effet plus adaptée aux basses vitesses réduites pour lesquelles les
amortissements aérodynamiques sont négatifs (pas de divergence) et relativement faibles en valeur
absolue (bonne excitabilité et nombre significatif de périodes pendant la décroissance du signal). La
86
méthode des oscillations forcées est moins appropriée à faible vitesse (sensibilité des capteurs) mais tout
à fait adaptée lorsque l’amortissement aérodynamique est élevé (l’amplitude des signaux a peu
d’importance) ou devient positif (mouvement est contrôlé).
3.5 Conclusion
Les techniques expérimentales présentées dans ce chapitre peuvent aujourd’hui encore être considérées
comme les méthodes de référence de caractérisation aérodynamique des tabliers de pont. Les méthodes
numériques, généralement basées sur la résolution des équations de Navier Stokes, sont en plein essor
mais encore fortement limitées par la puissance des ordinateurs actuels. On peut imaginer qu’un jour
viendra où celles-ci pourront se substituer, totalement ou partiellement, à l’expérimentation. Les
avantages sont évidents : coûts réduits, souplesse d’utilisation, non-intrusivité, accès à de nombreuses
données non accessibles expérimentalement,…
Nous avons consacré l’essentiel du chapitre à deux méthodes de mesure des coefficients aéroélastiques :
la méthode des oscillations libres et la méthode des oscillations forcées. Elles ont été appliquées à deux
profils de tablier de pont : un trapézoïdal et un rectangulaire de ratio 8. Seuls les résultats obtenus en
oscillations libres seront utilisés dans le cadre de la procédure de validation présentée au chapitre 5. La
raison est double :
 sens physique plus clair,

 gamme de vitesses réduites balayée plus appropriée.
87
Chapitre 4 Calcul de la
réponse dynamique
de l’ouvrage
4 Calcul de la réponse dynamique de l’ouvrage
Après avoir caractérisé les forces aérodynamiques s’exerçant sur une structure élancée, on se propose
d’en calculer la réponse. Deux approches sont développées, basées sur une discrétisation en éléments
finis de la structure.
La première consiste à raisonner de manière statistique en se plaçant dans l'espace des fréquences : on
parle alors d'approche spectrale. Cette méthode présente comme intérêt majeur d'être économique en
temps de calculs. Elle implique cependant de linéariser l'évolution des forces aérodynamiques en
fonction de l'incidence du vent, ce qui peut parfois se révéler très hasardeux. Nous proposons néanmoins
une méthode de "linéarisation équivalente" des coefficients aérodynamiques permettant de réduire ce
problème. Enfin, la méthode spectrale, basée sur des hypothèses de stationnarité, ne permet pas d'étudier
finement les problèmes d'instabilités aérodynamiques.
La seconde approche consiste à résoudre les équations de la mécanique dans le domaine temporel : on
parle alors d'approche temporelle. Cette approche, plus gourmande en temps de calcul, permet
néanmoins une meilleure modélisation des forces aérodynamiques. Contrairement à la méthode spectrale,
elle permet de prendre en compte leur non-linéarité. En outre, elle permet mieux d'appréhender les
phénomènes d'instabilités aérodynamiques.
89
Chapitre 4 Calcul de la réponse dynamique de l’ouvrage
4.1 Discrétisation en éléments finis
4.1.1 Généralités
La structure est modélisée selon la méthode des éléments finis. Son élancement lui permet d’être
facilement discrétisée en éléments de type « poutre 3D ». La figure 4.1 montre, à gauche, une maquette
du viaduc de Millau en phase de construction, à droite, sa modélisation en éléments poutres par la
méthode des éléments finis.
Figure 4.1 - Modélisation par les éléments finis du viaduc de Millau en phase de construction
Chaque élément poutre est défini par deux nœuds, soit 12 degrés de libertés (fig.4.2). Elle est en outre
caractérisée par un certain nombre de propriétés mécaniques (masse, inertie, raideur, amortissement) et
aérodynamiques (coefficients stationnaires et instationnaires).
Figure 4.2 - Modélisation d’un élément poutre 3D
90
Il est alors possible d'appliquer à ce type de modèle les outils classiques de la mécanique des vibrations.
Dans les quelques pages qui suivent nous rappelons brièvement la nature de ces outils. L'objectif n'est en
aucun cas de proposer quelque chose de nouveau mais de préciser les notions et notations utilisées par la
suite.
4.1.2 Mise en équations
a) Cas d’un élément
Considérons un élément quelconque d’une structure. Notons x le vecteur position d’un point de cet
élément et x n le vecteur contenant les coordonnées de ses nœuds. On introduit alors une coordonnée
paramétrique  liée à notre élément, et une fonction d’interpolation matricielle N reliant les deux
vecteurs :
x (  )  N ( ) x n    0 ;1  (4.1)
Si d P représente le vecteur déplacement en chaque point de l’élément et d n le vecteur contenant les

coordonnées des déplacements des nœuds, une approximation dite « isoparamétrique » permet d’écrire,:
d P ( )  N (  ) d n (4.2)
A la fonction matricielle d’interpolation N ( ) , on associe la matrice jacobienne J   permettant
d’effectuer le changement de variables, des coordonnées cartésiennes aux coordonnées paramétriques, au

cours d’une intégration . Dans le cas d’un élément poutre, la matrice jacobienne est un scalaire égal à L,
longueur de la poutre. Ainsi, pour une fonction quelconque f(x), on écrit :
 f (x) dx   J() f () d

élément élément
(4.3)
On introduit également la matrice X reliant le vecteur des déformations  au vecteur d n :
  X dn (4.4)
91
Ainsi, on peut écrire le principe des travaux virtuels pour un élément poutre :
 d &d& P dv 
 cd d& P dv 
d d
P *T P *T P *T P *T
H  dv  f v dv  f s ds  0
*T
(4.5)
élément élément élément élément élément
où H est la matrice de rigidité,  la masse volumique, c le coefficient d’amortissement, f s et f v les
vecteurs forces extérieures surfaciques et volumiques.
En reportant (4.2) et (4.4) dans (4.5), on obtient l’équation matricielle suivante :
m &d& n  c d& n  k d n  f n  0 (4.6)
m
 N
T
avec N dv la matrice masse de l’élément, (4.7)
élément
c
cN
T
N dv la matrice amortissement, (4.8)
élément
k
X
T
H X dv la matrice raideur, (4.9)
élément
fn 
 N f dv   N f ds
T T
v s le vecteur des forces équivalentes. (4.10)
élément élément
b) Assemblage
Il s’agit maintenant d’écrire l’équation matricielle régissant l’ensemble de la structure : c’est l’étape de
l’assemblage. Ainsi, si on note d le vecteur contenant l’ensemble des degrés de liberté de la structure i.e.
6 fois le nombre de nœuds utilisés pour sa discrétisation, on peut écrire :
M &d&  C d&  K d  F  0 (4.11)
où M, C, K sont les matrices de masse, d'amortissement et de raideur de la structure discrétisée obtenues

par assemblage des matrices m, c, k propres à chacun des élément.
92
La procédure d’assemblage des M, C, K est la suivante :
Notons ne le nombre d’éléments, ndlt le nombre de degrés de liberté total, m j , c j , k j les matrices masse,
amortissement et raideur associées à un élément j (j = 1, ne). Considérons maintenant les matrices carrées
m j ' , c j ' , k j ' de dimensions ndlt, contenant les valeurs des matrices m j , c j , k j , aux intersections entre
lignes et colonnes correspondant à des degrés de liberté appartenant à l’élément i, et des 0 partout
ailleurs. Alors, les matrices M , C , K peuvent être définies comme les sommes des matrices m j ' , c j ' ,
k j' :
M m '
j
j C c 'j
j K k '
j
j (4.12)
4.2 Calcul des modes propres de la structure
Reconsidérons l’équation (4.11) appliquée à un mouvement libre sans amortissement. On écrit alors :
M &d&  K d  0 (4.13)
Par analogie avec un système à un degré de liberté, on cherche des solutions de la forme :
d ( t )  Z e it (4.14)
où Z est un vecteur colonne constant.
On obtient ainsi l’équation suivante :
 M 2
 K  Z e it  0 (4.15)
L’existence d’une solution non triviale implique :
dét  M2  K   0 (4.16)
93
Il existe en général ndlt valeurs distinctes j2 solutions de cette équation algébrique en 2. Il reste à
résoudre l’équation matricielle qui suit :
K Z j  j M Z j
2
(4.17)
On obtient ainsi ndlt solutions à l’équation (4.15) caractérisées par les « vecteurs propres » (ou
« déformées ») Z j et les « valeurs propres » associées  j . Ces ndlt « modes propres » constituent une
2
base de l’espace des solutions de (4.15).
Les Z j étant définis à un facteur près, on les normalise généralement en fixant à 1 leur plus grande
composante. On introduit alors les notions de matrices masse et raideur généralisées M̂ et K̂ :
 M   M̂
T
(4.18)
 K   K̂
T
(4.19)
où  est une matrice dont les colonnes sont constituées des vecteurs propres Z j .
Dans la pratique, un nombre limité des premiers modes propres suffit à représenter les vibrations d’une
structure. On a alors intérêt à passer de la base nodale à une base modale réduite pour diminuer le rang
des systèmes à résoudre et par conséquent les durées de calcul.
94
4.3 La méthode spectrale
Les premières applications de la théorie des vibrations aléatoires à l’étude des structures soumises à
l’action de la turbulence du vent sont dues à Liepmann (1952). Elles concernaient particulièrement les
ailes d’avion minces. C’est Davenport (1961) qui, le premier, étendit les travaux de Liepmann à l’étude
des structures élastiques dans la couche limite atmosphérique. Un an plus tard, il utilisa ce nouvel outil
pour l’analyse des effets du vent sur les ponts suspendus. Il ne s'intéressait alors qu'aux problèmes de
turbulence. En 1974, Scanlan introduit les forces aéroélastiques dans le calcul spectral des ponts.
4.3.1 Principe
La méthode spectrale consiste à assimiler le mouvement de la structure étudiée à la superposition de

modes propres significatifs, et à évaluer pour chacun d'eux la densité spectrale de puissance des
déplacements. Chaque signal modal, supposé gaussien, peut alors être caractérisé par son écart-type,
calculé par simple intégration de son spectre.
Les étapes principales de cette approche sont :

 Discrétisation de la structure en éléments finis (poutres à 12 ddl) et calcul des modes propres;
 Caractérisation aérodynamique de l'enveloppe de la structure ;
 Caractérisation probabiliste du vent par l’évaluation de la matrice de densité interspectrale de
puissance de la composante turbulente ;
 Evaluation de la matrice densité spectrale de puissance des forces turbulentes par l’intermédiaire des
coefficients aérodynamiques dans la base nodale ;
 Constitution de la matrice densité spectrale de puissance des forces turbulentes dans la base modale ;
 Résolution de l’équation fondamentale de la dynamique dans la base modale et déduction de la
densité spectrale de puissance des oscillations pour chaque mode considéré ;
 Détermination des écarts-types des déplacements modaux par intégration des spectres ;
 Déduction des sollicitations extrêmes.
A partir de ce schéma de base, de nombreuses variations peuvent être envisagées, plus ou moins précises,
plus ou moins économiques en temps de calcul. Elles interviennent essentiellement aux niveaux suivants:
 Problème du couplage lors de la résolution du système dans la base modale ;
 Prise en compte de la déformée statique ;
 Linéarisation des coefficients aérodynamiques ;
 Intégration du spectre des déplacements.
95
4.3.2 Définition des repères
Une base « locale » est associée à chaque élément poutre ; c’est dans cette base que sont exprimés les
différents coefficients aérodynamiques dans le fichier d’entrée. Les axes x et y sont respectivement
associés aux efforts de traînée et de portance tandis que l’axe z est dirigé dans l’axe de la barre. Cela
correspond aux conventions généralement employées en soufflerie pour les maquettes sectionnelles. La
première composante du vecteur x et la troisième composante du vecteur y sont choisies positives.
z

Figure 4.3 – Repère local associé à un élément poutre
On définit également une base associée au vent ; l’axe x est dans le sens et la direction du vecteur vitesse
de vent moyen, tandis que les axes y et z sont respectivement associés aux fluctuations turbulentes v et
w . Le repère vent, ainsi que les repères locaux, sont définis au sein d’un repère général par des matrices
de transformation PV et PL :
 U V  UW 
 
 U V W U 2  V2 U 2  V2 U 2  V2  W2 
2 2 2
 V U  VW 
Pv   (4.20)
2 
 U V W U 2  V2 U V U V W 
2 2 2 2 2 2 2
 W U2  V2 
0
 U2  V2  W2 U V U V W 
2 2 2 2 2 

r
où U, V, W représentent les trois composantes du vecteur vent moyen U dans la base générale.
 signe ( Y )  Y  X Z 
  signe ( Y )  X 
 X 2  Y 2 X  Y
2 2

 signe ( Y )  X  Y Z
PL    signe ( Y )  Y  (4.21)
 X 2  Y 2 X 2  Y 2 
 
 0 X 2  Y 2  signe ( Y )  Z 
 
où X, Y, Z, représentent les coordonnées du vecteur directeur de l’élément dans le repère général.
96
4.3.3 Calcul de la matrice de densité spectrale des forces turbulentes
La vitesse du vent agissant sur un élément est exprimée dans une base quelconque définie par sa matrice
de passage PV . Soit PL la matrice de passage liée à la base locale. Alors, les composantes fluctuantes de
la vitesse du vent s'expriment, dans la base locale:
 uL  u
   
 v L   PL PV  v 
T
(4.22)
w  w
 L  
u, v, w étant définis dans la base liée au vent.
Exprimant (2.31) sous la forme suivante :
C    2 C ad (4.23)
où C ad est définie en (2.30) et  2 est une matrice diagonale contenant les termes  FX ,  FY ,  M Z , « racines
carrées » des fonctions d’admittance aérodynamique liées aux forces de traînée, portance et moment,
on déduit que le vecteur des forces par unité de longueur de l'élément s’écrit, dans la base locale :
u
 
f  C  PL
e T
t PV  v  (4.24)
w
 
Les coefficients aérodynamiques et leurs dérivées par rapport à l'incidence du vent sont recalculés pour
chaque élément en fonction de la direction du vent et de la déformée moyenne de la structure. Pour cela,
les coefficients et leurs dérivées sont entrés sous forme de polynômes du 4ème degré, fonctions de
l'incidence relative du vent.
Ces polynômes sont obtenus par optimisation non linéaire par les moindres carrés autour des valeurs
discrètes calculés par le programme de linéarisation équivalente des coefficients aérodynamiques. Nous
reviendrons sur les principes de cette méthode à la fin du chapitre. Les méthodes de calcul basées sur une
discrétisation en éléments finis nécessitent une formulation des efforts aux nœuds de la structure et non
97
sur ses éléments. On introduit donc une nouvelle base de projection, de dimension 12, associée aux
degrés de liberté des deux noeuds de l’élément poutre. Ainsi, le vecteur représentant les déplacements de
l’élément s’écrit sous la forme suivante:
d n  x 1 rz 2 
T
y1 z1 rx 1 ry1 rz1 x2 y2 z2 rx 2 ry 2 (4.25)
Le vecteur des déplacements en un point quelconque de l'élément, défini par sa coordonnée paramétrique
 [0;1], est alors égal à N () d n .
Le vecteur des forces élémentaires exprimé dans la base nodale de l’élément, s’exprime donc :
u
 

 N ( ) PV  v  J ( ) d
ne T T
f t C  PL (4.26)
e
w
 
ce qui donne, dans la base générale :
u
 

P N ( ) C  PL PV  v  J ( ) d
ng T T T
f t 12 (4.27)
e
w
 
où P12 est la matrice de changement de base général  local pour un vecteur à 12 composantes:
 PL T  0 33
 0 33
 03 3

 
 0 P  T
 0  0
P12   33 L 33 3 3 
(4.28)
 0  0 P  T
 0 
 33 33 L 3 3

 0 33  0  0 P  T
33 33 L 
En optant pour un schéma d’intégration linéaire, on écrit :
1   0 0 0 0 0  0 0 0 0 0
 0 1  0 0 0 0 0  0 0 0 0
 
 0 0 1  0 0 0 0 0  0 0 0
N ()    (4.29)
 0 0 0 1  0 0 0 0 0  0 0
 0 0 0 0 1  0 0 0 0 0  0
 
 0 0 0 0 0 1   0 0 0 0 0 
98
ce qui permet de réécrire (4.27) sous la forme :
u
 
 L  P12 N1 2 C  PL
ng T T T
f t PV  v  (4.30)
w
 
L étant la longueur de l’élément poutre.
Remarquons bien que le vecteur des forces élémentaires agissant sur l'élément est une transformation
linéaire des composantes turbulentes agissant sur l'élément. On peut maintenant construire le vecteur des
forces global Ftn par assemblage des vecteurs f tng . La jème coordonnée du vecteur Ftn s'exprime de la
façon suivante :
Ftn ( j )  f
e / ie
t
ng
 n  j, e  (4.31)
Pour tous les éléments e auquel appartient le degré de liberté j, on ajoute à la composante n(j,e) (numéro
local dans l'élément e du degré de liberté de numéro global j du vecteur) . C'est la procédure classique
d'assemblage.
Ce vecteur peut à présent être exprimé dans la base modale. Si l’on note Z j le vecteur propre du jème
mode, la jème coordonnée du vecteur des forces généralisées s’écrit :
Ftm ( j)  Z j Ftn (4.32)
u
 
On a donc exprimé le vecteur des forces généralisées comme une combinaison linéaire des vecteurs  v  .
w
 
u
 
Les vecteurs  v  étant aléatoires centrés stationnaires, ils admettent une matrice de densité spectrale
w
 
S V  . Le vecteur F admet par conséquent une matrice de densité spectrale, fonction quadratique des
matrices S V  .
99
On peut réécrire plus simplement (4.27) sous la forme:
u
 

C e (  )  v  J (  ) d
ng
f t (4.33)
w
élément
 
Soit S eV,e ' () la densité inter-spectrale des composantes turbulentes (u,v,w) exprimée dans le repère vent
et agissant sur e, et (u', v', w') exprimée dans le repère vent et agissant sur e' (celle ci est directement
donnée par le modèle de vent). On écrit alors :
 C () S

S ng n 'g ()  e
e,e '
V () C e ' ( ' ) J ( ) J (' ) d d ' (4.34)
f t ,f t
e e'
De manière analogue au cas vectoriel, la matrice de densité spectrale globale S Ft () s’obtient par
assemblage des S ng n 'g () :

f t ,f t
S Ft (  )  j , j'  S
e / je e ' / j'e '
ng n 'g
f t ,f t
 n ( j , e) , n ( j' , e' )  (4.35)
La matrice de densité spectrale des forces généralisée s’obtient alors en multipliant à gauche et à droite
par la matrice des déformées modales  :
S Ft ()   S Ft () 
T
(4.36)
4.3.4 Calcul de la matrice de densité spectrale des forces aéroélastiques
Considérons encore un élément poutre défini en chacun de ses points par sa coordonnée paramétrique .
Le vecteur des forces aéroélastiques exprimé dans un repère local vaut, à l’instant t :
 0 
 e 
 l ae (, t ) 
 0 
f aee (, t )    (4.37)
 0 
 0 
 
 m e (, t ) 
 ae 
100
avec:
1 2 K KB * K2 * 
l eae (, t )   U B H 1* h& , t   H 2 & , t   K 2 H *3 , t   H 4 h , t  (4.38)
2 U U B 
1 K KB * K2 * 
m eae (, t )   U B 2  A1* h& , t   A 2 & , t   K 2 A *3 , t   A 4 h , t 
2
(4.39)
2 U U B 
Les équations (4.38) et (4.39) peuvent être réécrites sous la forme matricielle suivante :
f aee (, t )  C eae ( ) d& e (, t )  K eae ( ) d e (, t ) (4.40)
0 0 0 0 0 0 
 K * KB * 
0 U
H1 0 0 0
U
H2 
0 0 0 0 0 0 
1
avec C eae  U 2 B  (4.41)
2 0 0 0 0 0 0 
0 0 0 0 0 0 
 KB * KB2 * 
0 A1 0 0 0 A2 
 U U 
0 0 0 0 0 0 
0 K 2 BH * 0 0 0 K H 
2 *
 2 3
1  0 0 0 0 0 0 
K eae  U 2 B  (4.42)
2 0 0 0 0 0 0 
0 0 0 0 0 0 
 2 * 
0 K A 4 0 0 0 K BA *3 
2
On écrit alors, dans la base nodale :
   
 
f aeen ( t )   N   C eae N  J   d  d& n   N   K eae N  J   d  d n
T T
(4.43)
 e   e 
On en déduit l’expression de deux nouvelles matrices raideur et amortissement dans le repère général :
 

C aeng  P12  N   C eae N   J   d  P12
T T
(4.44)
 e 
101
 

K aeng  P12  N   K eae N   J  d  P12
T T
(4.45)
 e 
que l’on assemble suivant la procédure classique :
C ae (  )  C
e
ng
ae   (4.46)
K ae (  )  K
e
ng
ae   (4.47)
4.3.5 Mise en équations
Les déplacements d’une structure soumise à l’action turbulente du vent sont donc régis dans le temps par
l’équation différentielle matricielle suivante :
M &d& ( t )  C d ( t )  K d ( t )  Fae ( U , d& , d, t )  Ft ( U, V, t ) (4.48)
avec M : matrice des masses
C : matrice d'amortissement
K : matrice de raideur
d : vecteur des déplacements
Fae : vecteur des forces aéroélastiques
Ft : vecteur des forces turbulentes
Dans le domaine fréquentiel, l’équation peut être développée sous la forme suivante (i² = -1) :
( 2 M  i C  K ) d ()   i  C ae U,   K ae U,   d ()  Ft ( U, V, ) (4.49)

 
Ce système de rang n, où n est le nombre de degrés de liberté de la structure i.e. 6 fois le nombre de
nœuds peut être ramené à un système de rang inférieur m en effectuant une projection sur une base
constituée des m premiers vecteurs propres.
102
Introduisant  la matrice de dimension n x m composée des m premiers vecteurs propres disposés en
colonnes, on peut alors écrire :

d ()   q () (4.50)
où q est le vecteur des coordonnées généralisées.
En insérant l’équation (4.50) dans l’équation (4.49) et en multipliant à gauche tous les termes par  , on
T
obtient l’équation suivante:
( 2 M̂  i Ĉ  K̂ ) q ()   i  Ĉ ae U,   K̂ ae U,   q ()   Ft 

T
(4.51)
 
M̂   M  matrice des masses généralisée (diagonale)

T
avec
Ĉ   C 
T
matrice d’amortissement généralisée (non diagonale)
K̂   K 
T
matrice des raideurs généralisée (diagonale)
Ĉ ae U,    C ae U,  
T
matrice d’amortissement aérodynamique (non diagonale)
K̂ ae U,    K ae U,  
T
matrice des raideurs aérodynamique (non diagonale)
4.3.6 Analyse de la stabilité
Avant de chercher à déterminer la réponse de notre structure à l’action turbulente du vent, il convient
d’en vérifier la stabilité. On constate dans la pratique que les phénomènes d’instabilité se développent
d’autant plus facilement que le taux de turbulence est faible. C’est pourquoi, dans le cadre d'une
approche sécuritaire, l’étude de la stabilité d’un tablier de pont, qu’elle soit expérimentale ou numérique,
doit être effectuée en vent uniforme.
Le caractère instable d’un tablier est essentiellement caractérisé par sa vitesse critique de flottement,
vitesse au-delà de laquelle les oscillations de la structure prennent un caractère divergent. Il convient
cependant de bien distinguer deux aspects du phénomène d’instabilité : son développement proprement
dit, mais également son initiation. En effet, toute instabilité doit être initiée par une impulsion
suffisamment puissante pour pouvoir ensuite se développer. Dans la nature, un minimum de turbulence
103
(ou éventuellement un détachement tourbillonnaire) est donc nécessaire pour observer un comportement
divergent.
Ce qui nous intéresse, c’est de savoir à quel moment la structure sera potentiellement susceptible
d’adopter un comportement divergent. La question de l’initiation de cette divergence ne nous intéresse
donc pas ici.
Reprenant l’équation (4.40), et imposant Ft ( )  0 , l’instabilité se définit mathématiquement comme le
moment, où, pour une vitesse donnée U, il existe une valeur de  annulant la matrice G() définie par:
G ()   2 M̂  i  Ĉ  Ĉ ae ( U , )    K̂  K̂ ae ( U , )  (4.52)
   
La vitesse la plus faible, permettant de satisfaire le critère d’instabilité, sera supposée correspondre à la
vitesse critique de flottement. La méthode la plus appropriée pour déterminer cette vitesse est une
méthode du type « point fixe » nommée pK-F. Elle a été développée par Namini (1992). Supposant que
les m pulsations j des m premiers modes ont été calculées, on peut suivre l'algorithme proposé par
Patron Solares (1998):
On recherche des solutions de la forme s = (+i) à l’équation :
G (s)   s 2 M̂  s  Ĉ  Ĉ ae Im(s)    K̂  K̂ ae Im(s)   0 (4.53)

   
La procédure qui suit est réalisée pour des vitesses croissantes Uk, (k=0,…, f) :
1. Calcul de deux valeurs initiales de s :
s j1  U k   0,01  j  i  j (4.54)
s j2  U k   i  j (4.55)
2. Calcul des coefficients instationnaires de Scanlan correspondant aux fréquences réduites :
B Im( s j1 )
K j1  (4.56)
Uk
104
B Im( s j2 )
K j2  (4.57)
Uk
3. Calcul du déterminant de G(s) correspondant à chaque valeur de s j1  U k  et s j2  U k 
4. Définition d’une nouvelle valeur de s en utilisant un schéma linéaire d’interpolation :
s j2 ( U k ) G (s j1 )  s j1 ( U k ) G (s j2 )
s j3 ( U k )  (4.58)
G (s j1 )  G (s j2 )
5. Répétition des étapes 3-4 jusqu’à satisfaction du critère de convergence :
G (s j2 )  tolérance1
Pour chaque nouvelle itération, s j1  U k   s j2  U k  et s j2  U k   s j3  U k 
6. Lorsque le test de l’étape 5 est vérifié, la réponse oscillatoire aéroélastiquement modifiée a été
trouvée. En outre, si la condition :
Re(s j2 )  tolérance 2
est vérifiée, le flottement intervient pour la pulsation Im(s j2 ) et la vitesse Uk est une vitesse critique
de flottement.
« tolérance 1 » et « tolérance 2 » sont des seuils de convergence à définir.
4.3.7 Résolution complète
On cherche maintenant à calculer la réponse d’une structure élancée soumise à l’action turbulente du
vent. L’équation (4.51) permet, pour une pulsation  donnée, de déterminer la réponse de la structure
étudiée dans la base modale :
q ()  H ()  Ft 

T
(4.59)
105
où H () est la matrice de transfert du système définie par :
1
 
H ()    2 M̂  i Ĉ  Ĉ ae ( U, )    K̂  K̂ ae ( U, )   (4.60)
     
Dans le cas général, la matrice H() n’est pas diagonale ; son inversion pour différentes valeurs de 
représente donc des temps de calcul relativement importants. C’est néanmoins la seule façon de faire
intervenir l’intégralité des phénomènes aéroélastiques dans le domaine spectral.
La densité spectrale de puissance de chaque coordonnée généralisée s’obtient en multipliant à droite

l’équation (4.59) par son expression adjointe :
S q ()  H()  S Ft ()  H () T*

T
(4.61)
Une fois connue la matrice de densité spectrale des déplacements modaux, on peut retourner dans le
domaine temporel et estimer leurs variances :
  q1q1 2 K  q1qm 
2
  1 
 M
 2
O M 
K  qmqm 
2 2  S () d
0
q (4.62)
 qmq1
La matrice variance des déplacements nodaux peut alors être calculée de la manière suivante:
  d1d1 2 K  d1d n 2    q1q1 2 K  q1q m 2 

    T
 M O M    M O M  (4.63)
 d d 2 K  d d 2   q q 2 K  q q 2 
 n1 n n   m1 m m 
Les termes diagonaux représentent les carrés des écarts-types de chaque degré de liberté de la structure ;
les autre termes représentent les termes de covariance des différentes paires de degrés de liberté.
4.3.8 Résolution simplifiée
Le principe consiste à reprendre l’équation matricielle (4.51) et à en éliminer les termes non diagonaux.
Cela revient à éliminer tout phénomène de couplage structurel et aérodynamique entre les différents
106
modes de la structure. En ce qui concerne l’amortissement structurel, cette simplification ne change en

réalité pas grand chose dans la mesure où les termes couplés sont très difficiles à calculer et donc très
mal connus. Par contre, les termes aéroélastiques n’apparaissent plus que sous forme d’amortissement
modal ajouté à l’amortissement structurel et un glissement en fréquence.
Cette simplification est tout à fait justifiée dans le cas où les fréquences des différents modes sont bien
détachées. Dans le cas de ponts suspendus ou à haubans "vraie grandeur", lorsque les fréquences des
modes significatifs sont proche, les termes de couplage peuvent prendre un rôle plus important, et
l'hypothèse simplificatrice perdre de sa pertinence.
Le système matriciel (4.51) se résume donc à m systèmes indépendants à un degré de liberté :
( 2 m g  i c g  k g ) q ()  f tg  (4.64)
où m g , c g , k g , f tg sont respectivement les masse, amortissement, raideur et force généralisées.
On réécrit sous forme normalisée :
f tg 
( 2  2 i   f  f   f ) q   h  q  
2
(4.65)
mg
avec  f   0   et  f   s   a
Deux solutions s’offrent maintenant à nous. La première consiste à calculer Sq() pour  variant de « 0 à
 » puis à estimer l’écart-type des déplacements modaux par l’intégration classique :
1 
q 
S q () d
2
(4.66)
2 0
La seconde consiste à considérer la solution q comme la somme d’une contribution quasi-statique

(excitations de basses fréquences) et d’une contribution dynamique (excitation à la fréquence de
résonance).
La contribution dynamique est calculée en remplaçant l’excitation réelle par un bruit blanc d’intensité
égale à la densité de puissance spectrale de l’excitation à la fréquence naturelle du système. Cela
donne (Preumont, 1990):
107
2  S f  f 
 dyn
q  (4.67)
2  m g f
2 3
Les conditions sous lesquelles l’approximation se justifie sont :
 Les valeurs de S f  sont telles que la contribution dominante de l’intégrale provient du proche
voisinage de f ;
 S f  varie lentement autour de la pulsation naturelle f.
La contribution quasi-statique est calculée en assimilant l'ensemble du spectre des sollicitations à une
sollicitation statique, c'est-à-dire associé à une admittance mécanique unitaire.
f
2
2
 qsq  (4.68)
m g f
2
Cette approximation est a priori beaucoup plus grossière que la précédente ; en outre une partie de la
contribution quasi statique est déjà incluse dans la contribution dynamique. Mais cela n’a pas de
répercussions fâcheuses sur la précision des résultats. En effet, l'importance de la contribution quasi-
statique par rapport à la contribution dynamique est généralement faible.
Un schéma de principe de la méthode spectrale simplifiée est présenté figure 4.4.
108
SU

0
AA
Calcul de la Calcul de la
composante composante
quasi-statique dynamique

0
SF SF SF
SF0 SF0 SF0

  
0 0
AM AM AM
  
0 0
SD SD SD
  
0 0 0
Figure 4.4 – Schéma de principe de la méthode spectrale : calcul des écarts-types

des composantes dynamique et quasi-statique des déplacements modaux
SU : densité spectrale de puissance du vent
SF : densité spectrale de puissance des forces
SD : densité spectrale de puissance des déplacements
AA : admittance aérodynamique AM : admittance mécanique
0 : pulsation de résonance du mode considéré
109
4.4 Simulation numérique d’un signal de vent tridimensionnel

spatialement corrélé
Le calcul temporel de la réponse d’un structure élancée sous l’action du vent nécessite une formulation
temporelle des forces aérodynamiques. En particulier, on doit pouvoir imposer, sur chaque élément de la
structure, un signal de forces turbulentes respectant des propriétés statistiques bien précises. Cela nous
amène à générer un signal temporel tridimensionnel de vent dont les propriétés sont exprimées dans le
domaine fréquentiel sous forme de matrices de densité spectrale.
Trois raisons nous amènent à insérer ce sujet entre la présentation de la méthode spectrale et la
présentation de la méthode temporelle :
 Tout d’abord, il ne s’agit pas de calcul temporel au sens où nous l’entendons mais plutôt de
génération de conditions d’entrée à destination du calcul temporel ;
 ensuite, il nous semble que ce sujet constitue une véritable transition entre les approches spectrales et
temporelles ;
 enfin, la portée des méthodes présentées par la suite dépasse largement le cadre des structures
souples. Les applications sont multiples : génération de vents traversiers pour le calcul de la tenue au
vent des trains, étude du comportement dynamique des grues à tour, génération de conditions limites
pour des calculs C.F.D. type L.E.S ou D.N.S., …
De nombreuses méthodes numériques existent permettant de générer des signaux multidimensionnels

spatialement corrélés. La plus utilisée dans le domaine du « wind engineering » est celle développée par
Shinozuka et Lin en 1972. Cette méthode présente surtout l’avantage d’être simple à mettre en œuvre et
peu coûteuse en temps de calculs. Néanmoins ses performances sont un peu décevantes dans le domaine
des basses fréquences.
C’est pourquoi nous proposons également une autre méthode de simulation, plus récente, aux fondements
théoriques plus solides, mais également plus lourde en terme de temps de calculs. Cette méthode, dite de
« simulation conditionnelle », est basée sur l’algorithme proposé par Jin, Lutes et Sarkani (1997).
110
4.4.1 Méthode de Shinozuka
Cette méthode permet de générer en n points de l’espace n processus corrélés uk avec, en entrée, un pas
de temps t, une durée de simulation T, et les variables permettant de définir les spectres et interspectres
entre chacun des signaux. Dans notre cas, il s'agit de la vitesse moyenne du vent, des écart-types et
échelles de la turbulence, et des coefficients de la cohérence.
Dans un premier temps, on détermine la fréquence de coupure fc telle que la quasi-totalité de l’énergie du
processus u(t) soit localisée dans la bande fréquentielle [-fc +fc] i.e. telle que
 fc 
fc 
S u (f ) df   S u (f ) df

(4.69)
avec  suffisamment grand. Pour un vent moyen vraie grandeur, une valeur de 0,95 pour  donne une
fréquence de coupure égale à environ 3Hz. Dans la pratique, un autre critère permettant de définir la
fréquence de coupure est la fréquence du dernier mode de la structure considéré. Dans le cas d’un pont
comme celui de Rion-Antirion, la gamme de fréquence des 20 premiers modes considérés est incluse
dans l’intervalle [0,1Hz 0,5Hz]. Autant dire que bien avant 3 Hz il ne se passera plus rien
dynamiquement parlant.
La fréquence de Nyquist définie par
1
fN  (4.70)
2 t
doit, d’après la règle d’échantillonnage de Shannon, vérifier l’inégalité suivante :
f N  fc (4.71)
Bien entendu, plus fN est élevé plus t est petit et donc, pour une durée de simulation donnée, plus le
nombre de points à simuler sera élevé. Dans la pratique, le compromis généralement fait consiste à
transformer (4.71) en égalité, ce qui permet d’écrire :
1
t  (4.72)
2 fc
111
L’espace des fréquences [-fc +fc] est alors discrétisé suivant le pas :
2f c
f  (4.73)
n
Pour chaque pas en fréquence on calcule la matrice de densité spectrale du vent turbulent Sif  d’après
les équations (1.13). Cette matrice étant réelle et symétrique, on peut lui appliquer la décomposition de
Cholesky :
Sif   H if  H if 

T
(4.74)
où H if  est une matrice triangulaire inférieure.
Les éléments de la matrice H if  peuvent être considérés comme les coefficients d’une combinaison
linéaire dans une base de n processus unitaires et uniformément distribués (bruit blanc par exemple). Ce
processus unitaire peut être écrit, pour chaque fréquence if sous la forme d’un vecteur colonne
X (if ) de dimension n :
 ei1i 
 
 M 
X (if )  e iki  (4.75)
 
 M 
e ini 
 
 ki est la phase associée au point k pour la ième fréquence if . C’est une variable aléatoire uniformément
distribuée sur [0 2]. Ceci nous assure la convergence du processus simulé u vers un processus gaussien
quand n tend vers l’infini.
Soit Y (if ) le vecteur colonne défini par :
Y (if )  H (if ) X (if ) (4.76)
Le processus temporel uk est alors défini par la transformée de Fourier inverse de Y k (if ) .
Afin d’améliorer l’aspect des processus ainsi générés, on peut multiplier préalablement chaque matrice
Sif  par une variable aléatoire ‫א‬i admettant une répartition gaussienne, unitaire et centrée.
112
4.4.2 Méthode de Jin, Lutes et Sarkani (1997)
Cette méthode de simulation appartient à la famille des simulations « conditionnelles » ; elle consiste à
générer un nouveau signal temporel dont on connaît les propriétés statistiques intrinsèques ainsi que le
degré de corrélation avec une série de signaux temporels existants ; ces signaux existants pouvant eux-
mêmes être synthétiques ou issus de données expérimentales.
a) Hypothèses et notations
Supposons connus n-1 processus temporels tridimensionnels en n-1 points de l’espace. Il s’agit
 
maintenant d’en générer un nième; notons u n ( t )   u 1n ( t ), u n2 ( t ), u 3n ( t )
T
ce nouveau processus vectoriel
et u n ( t )   u 1n ( t ), u n2 ( t ), u 3n ( t ) 
T
la valeur prise à l’instant t par ce processus; on suppose également
connues ses fonctions densité spectrale et interspectrales de puissance que l’on note respectivement
S n () et S nm (), m variant de 1 à n-1.
La composante u n
1 ( t ) est simulée en premier, à partir de la connaissance des signaux existants
u m
( t ) ; on génère ensuite  u ( t ) n
2 à partir de  u 1n ( t ) et des u m
( t ), et enfin u n
3 ( t )à partir de
u n
1 ( t ),  u ( t ), et des  u ( t ).
n
2
m
Les processus simulés étant supposés gaussiens, la distribution conditionnelle5 de chaque variable
aléatoire u nr ( t ) est connue si sa moyenne et sa variance conditionnelles sont connues. L’ensemble du
processus  u nr ( t ), lui, est entièrement défini par sa moyenne conditionnelle et sa fonction de covariance
conditionnelle (ou sa fonction densité spectrale de puissance conditionnelle si l’on raisonne dans le
domaine fréquentiel).
Rappelons également que, dans le cas de processus gaussiens corrélés, la moyenne d’ensemble
conditionnelle d’une variable aléatoire, représentant la valeur prise par un processus à un instant t, peut
être considérée comme une combinaison linéaire des valeurs prises par les autres processus, mais que la
covariance conditionnelle de cette même variable ne doit obéir qu’à la fonction de covariance
conditionnelle propre au processus considéré et ne dépendant pas des valeurs prises par les autres
processus.
5
On parle de distribution conditionnelle d’un processus (resp.d’une variable aléatoire) lorsque la valeur prise par ce
processus à un instant t (resp. cette variable lors d’un tirage t) est dépendante des valeurs prises par d’autres
processus resp. variables) à cet instant t (resp. lors de ce tirage t).
113
On peut ainsi écrire la composante u nr ( t ) sous la forme suivante :
u nr ( t )  g nr ( t )  Yrn ( t ) (4.77)
où g nr ( t ) représente la moyenne conditionnelle et Yrn ( t ) une variable gaussienne centrée indépendante

des composantes précédentes.
La relation (4.77) peut être réécrite dans le domaine spectral par simple transformée de Fourier :
û nr ()  ĝ nr ()  Ŷrn () (4.78)
b) Evaluation de la moyenne conditionnelle dans le cas général
La moyenne conditionnelle du vecteur u n ( t ) peut être formulée ainsi :



E u 1n ( t )  u 1 ( t ),  u 2 ( t ),...,  u n 1 ( t )  

g n ( t )   g 1n ( t ), g n2 ( t ), g 3n ( t )  = 
T

E u n2 ( t )  u 1 ( t ),  u 2 ( t ),...,  u n 1 ( t ),  u 1n ( t )   (4.79)

 
 E u 3n ( t )  u 1 ( t ),  u 2 ( t ),...,  u n 1 ( t ),  u 1n ( t ),  u n2 ( t )




Sa dépendance linéaire par rapport aux signaux déjà simulés peut s’exprimer sous la forme
générale suivante:
n 

1
g n (t )  W (  ) u m ( t   ) d
nm
(4.80)
2 m 1

Les W ( ) sont les matrices de pondération 3*3, pleines pour m < n ; les termes de W () , eux, sont
nm nn
nuls, excepté sous la diagonale.
Il s’agit maintenant de déterminer les valeurs des 9  ( n  1)  3 fonctions de pondération W jknl ( ) . Pour
cela, on développe les fonctions de corrélation entre g nr ( t ) et les différents signaux générés
précédemment :
114
n 1 3  r 1 
R np n

rk d ( )  E g r ( t ) u k ( t  ) 
p
 1
2   Wrjnm () R mp

jk (    ) d 
1
2  W 
nn
rj () R npjk (  ) d
m 1 j1 j1
(4.81)
avec pour p = 1,…, n-1 k = 1, 2, 3
pour p = n k = 1,…, r-1
On applique la transformée de Fourier à l’équation (4.81), ce qui conduit à la formulation algébrique

suivante :
n 1 3 r 1
S () 
np
rk  Ŵ
m 1 j1
nm
rj () S () 
mp
jk  Ŵj1
nn
rj () Snpjk () (4.82)
avec pour p = 1,…, n-1 k = 1, 2, 3

pour p = n k = 1,…, r-1
On se retrouve ainsi, pour une valeur de  donnée avec un système linéaire à 3  ( n  1)  r  1
équations et autant d’inconnues Ŵrjnm () et Ŵrjnn () . Sous forme matricielle, cela donne :
 S11 11 S12 11 S13 11 L L S1n 11 1

S n2 11 1
S 3n 11 1
S1n 11 L S nr 11 1
  Ŵ nr 11   S nr 11 
 11   n1   n1 
 S1 2 S12 12 S13 12 L L S1n 21 1
S n2 21 1
S 3n 21 1
S1n 21 L S nr 11 2   Ŵ r 2   S r 2 
 S11 13 S12 13 S13 13 L L S1n 31 1
S n2 31 1
S 3n 31 1
S1n 31 L S nr 11 3   Ŵ n 1   S n 1 
   r3   r3 
 M M M O M M M M M M   M   M 
 M M M O M M M M M M   M   M 
 1 n 1     
1 n 1 1 n 1 n 1 n 1 n 1 n 1 n 1 n 1 n n 1
S1 1 S 21 S 31 L L S1 1 S 21 S 31 S 11 L S nr 1n 11   Ŵ nr 1n 1   S rn 1n  1 
S1 n 1  n n 1   n n 1 
 12 S 1 n 1
2 2 S 1 n 1
3 2 L L S n 1 n 1
1 2 S n 1 n 1
2 2 S n 1 n 1
3 2 S n n 1
1 2 L S nr 1n 21   Ŵ r 2  S r 1 
S11 3n 1 S 1 n 1
2 3 S 1 n 1
3 3 L L S n 1 n 1
13 S n 1 n 1
2 3 S n 1 n 1
3 3 S n n 1
13 L S nr 1n 31   Ŵ nr 3n 1  S rn 1 n 1 
 1n 1 n 1 n n 1 n n 1 n n 1 n n n   n n   nn 
 S1 1 S 21 S 31 L L S 11 S 21 S 31 S 11 L S nr 1n 1   Ŵ r 1   S r 1 
 M M M L L M M M M O M   M   M 
 1n 1 n 1 n n 1 n n 1 n n 1 n n n n n
  n n   nn 
S1 r 1 S 2 r 1 S 3 r 1 L L S 1 r 1 S 2 r 1 S 3 r 1 S 1 r 1 L S r 1 r 1   Ŵ r r 1   S r r 1 
(4.83)
Une fois ce système résolu, on calcule aisément la moyenne conditionnelle dans le domaine fréquentiel :
 Ŵ
1
ĝ () 
n nm
() û m () (4.84)
2 m 1
115
Une simple transformée de Fourier inverse permet de remonter au signal temporel g n ( t ) .
c) Evaluation de la variable aléatoire Y
Les termes Ŷrn () peuvent alors être calculés en insérant (4.84) dans (4.78) :
n 3
Ŷ ()  û () 
r
n n
r  Ŵ
m 1 k 1
nm
rk () û mn
kr () (4.85)
En multipliant chaque membre de l’équation par son conjugué, puis en appliquant une moyenne
d’ensemble, on obtient l’expression de la densité spectrale de puissance de Yrn ( t ) :
 T  nn
 
n 3
E Ŷ ()Ŷ () 
r
n
r
n*
S rr () 
2  
m 1 k 1
Ŵrknm () S mn
kr () 

(4.86)
Il nous reste à générer un signal aléatoire respectant la fonction densité spectrale de puissance déduite de
(4.86) ; pour cela, on choisi d’utiliser la forme suivante (Lutes et al. 1996):
Ŷrn ()  B r e j r r = 1, 2, 3 (4.87)
où r est une variable aléatoire uniformément distribuée sur [0, 2] et Br une variable aléatoire,
indépendante de r suivant une distribution de Rayleigh avec :
  
E B r  E Ŷrn  Ŷrn 
2 *
 (4.88)
d) Algorithme de résolution
On comprend bien que la résolution du système (4.83) va occuper l’essentiel du temps de calcul. C’est
donc à ce niveau qu’il nous faut chercher à optimiser la procédure algorithmique.
Le système (4.83) peut s’écrire sous la forme synthétique suivante :
A J  () X J  ()  a J  ()

T
(4.89)
116
où J représente le rang du système matriciel ; si N représente le nombre total de points, J varie de 1 à

3N-1. Si l’on considère que la complexité de résolution du système (4.89) est de l’ordre de J3, la
complexité de la résolution de l’ensemble des systèmes, pour J variant de 1 à 3N-1, est a priori de l’ordre
de (3N)4. Cependant, les auteurs proposent un algorithme tirant partie de la spécificité du problème et
permettant de ramener l’ordre de la complexité à (3N)2.
L’astuce consiste à profiter de la relation existant entre les matrices A J  () et les vecteurs a J  () pour
T
différentes valeurs de J. En particulier, A J 1 () se déduit de A J  () par l’adjonction d’une ligne et
T T
d’une colonne dont les J premiers éléments sont issus de a J  () (conjugués pour la colonne). En fait, seul
un nouveau scalaire est rajouté :
 A J  () a J  () 
A TJ 1 ()   *T  (4.90)
a () a J 1 ()
 J  
Pour résoudre (4.90), on commence par appliquer une décomposition de Cholesky à la
matrice A J  () (qui est hermitienne) :

T
A J  ()  L J  ()L J  ()

T *T
(4.91)
où L J  () est une matrice triangulaire inférieure.
Le système (4.91) peut ainsi être décomposé en deux sous-systèmes triangulaires résolus
successivement :
L J  () Z J  ()  a J  () (4.92)
L J  ()X J  ()  Z J  ()

*T
(4.93)
On vérifie facilement les propriétés suivantes :
 La matrice L J 1 () se déduit de L J  () par adjonction d’une nouvelle ligne,
 Les J premiers éléments de cette nouvelle ligne sont ceux de Z J  () .

*T
117
Ainsi, connaissant L J  () , on peut déduire L J 1 () par la formulation qui suit :
 L J  () 0 
L J 1 ()   *T
Z () l J 1 ()
(4.94)
 J  
Reste à calculer le terme l[J+1](). Pour cela on écrit la matrice A J 1 () en combinant (4.91) et (4.94).
T
Cela donne :
 L J  () L*tJ  () L J  () Z J  () 

A J 1 ()   *t
t
2  (4.95)
 Z J  ()L J  () Z J  () Z J  ()  l J 1 () 
*t *t
En comparant (4.95) et (4.90), on déduit l[J+1]() :
l J 1 ()  a J 1 ()  Z *TJ  () Z J  () 2

(4.96)
On peut ainsi décomposer la procédure numérique, pour une valeur de  donnée, comme suit :
 On applique la décomposition de Cholesky à A 2  () puis on calcule Z 2  () et X 2  () ;

T
 Pour J variant de 3 à 3N-1 :
on utilise (4.94) et (4.96) pour calculer L J  () d’après L J 1 () , Z J 1 () , et l[J]() ;
T T
on en déduit Z J  () et X J  () i.e. les fonctions de pondération W() ;
 Pour r variant de 2 à 3, on utilise (4.84) pour générer ĝ 1r () ;
On calcule ensuite Ŷr1 () d’après (4.87) ;
On en déduit û 1r () à partir de (4.78) ;
 Pour n variant de 2 à N, et r variant de 1 à 3, on utilise (4.84) pour générer ĝ nr () ;
118
On calcule ensuite Ŷrn () d’après (4.87) pour r = 1 à 3;
On en déduit û nr () à partir de (4.78).
On réitère l’ensemble de cette procédure sur une plage de pulsation jugée nécessaire et suffisante pour le
cas étudié. La problématique est évidemment la même que pour la méthode de Shinozuka. On obtient
ainsi 3N processus dans le domaine spectral. Par simples transformées de Fourier inverses, on en déduit
les 3N processus dans le domaine temporel.
4.4.3 Comparaison entre les deux méthodes
Afin de justifier pleinement l’utilisation de l’algorithme de Sarkani au détriment de celui de Shinozuka, il

convient de réaliser un test comparatif entre les deux méthodes. La configuration utilisée pour cette
comparaison est la suivante:
 Deux points situés à la même altitude avec un écartement variable dans la direction transversale;
 Vitesse moyenne du vent U  10 m / s ;
 Echelle longitudinales de la turbulence Lxu  0,5 m , Lxw  0,1 m ;
 Coefficients de la cohérence C yu  10 , C yw  10 ;
 Ecarts types de la turbulence  u  1,5 m / s ,  w  1 m / s
On compare d’une part les courbes de densité spectrale obtenues à l’aide des méthodes de Shinozuka et
Sarkani avec les courbes théoriques issues des formulations de Karman, d’autre part la fonction de
cohérence entre les deux signaux avec la fonction de cohérence théorique de Davenport. Les résultats
présentés figure (4.5) et (4.6) correspondent à deux signaux de composante longitudinale de vent espacés
de 0,4 m. Ils sont tout à fait représentatifs de l'ensemble des résultats obtenus.
En ce qui concerne la fidélité aux conditions spectrales imposées, on constate que l'algorithme de Sarkani
est légèrement plus performant que celui de Shinozuka ; cependant la différence n'est pas frappante. Par
contre, la fonction de cohérence est très bien simulée par la méthode de Sarkani alors que celle de
Shinozuka la sous-estime très largement.
Notons que la performance de la méthode de Sarkani n'est pas très surprenante, étant donné que la base
"déterministe" de chaque signal est construite d'après les contraintes en terme d'interspectres avec les
autres signaux.
119
0.4
Jin, Lutes & Sarkani
Shinozuka
0.3 Karman
n.S(n)/sigma²
0.2
0.1
0
0 5 10 15 20
f (Hz)
Figure 4.5 – Comparaison des spectres obtenus à partir des méthodes de Shinozuka et Sarkani
1
Jin, Lutes & Sarkani
Shinozuka
0.8
Davenport
0.6
Coherence
0.4
0.2
0
0 5 10 15 20
f (Hz)
Figure 4.6 – Comparaison des cohérences obtenues à partir des méthodes de Shinozuka et Sarkani
120
4.5 La méthode temporelle
4.5.1 Mise en équation
a) Formulation des forces turbulentes
La formulation des forces turbulentes s’exprime directement à partir de la connaissance du signal de vent
tridimensionnel précédemment calculé. On a vu que la linéarisation des forces aérodynamiques était peu
pertinente pour certains profils de tabliers. Mais si cette approximation était nécessaire pour la méthode
spectrale, il n’en est plus de même pour la méthode temporelle. On peut en effet, à chaque instant, et
pour chaque élément, calculer l’incidence du vent dans le repère local et assigner de nouvelles valeurs
aux coefficients aérodynamiques.
Pour cela, on approche les courbes d’évolutions des coefficients stationnaires par des fonctions
polynomiales du 4ème degré. Nos trois coefficients stationnaires s’écrivent donc sous la forme :
C X (  )  C 1X  C 2 X   C 3 X  2  C 4 X  3  C 5 X  4 (4.97)
C Y (  )  C 1Y  C 2 Y   C 3 Y  2  C 4 Y  3  C 5 Y  4 (4.98)
C M (  )  C 1M  C 2 M   C 3 M  2  C 4 M  3  C 5 M  4 (4.99)
Le vecteur buffeting forces instantané s’écrit donc :
1 ~
f te ( t )  BU 2 C 12345 ( t ) (4.100)
2
( t )   1   2 4 
T
avec 3 (4.101)
 C 1X C 2X C 3X C 4X C 5X 
C C 2Y C 3Y C 4Y C 5 Y 
 1Y
 0 0 0 0 0 
et C12345   (4.102)
 0 0 0 0 0 
 0 0 0 0 0 
 
C1M C 2M C 3M C 4M C 5 M 
121
On exprime alors ce vecteur dans le repère général et la base nodale de l’élément :
T
1
f tng ( t )  L  P12 N   f te ( t )
T
(4.103)
 2
Dans l’absolu, la fonction d’admittance aérodynamique devrait intervenir ici sous forme de produits de
convolution. On préfère cependant, dans un soucis d’économie de temps de calcul, la faire intervenir
mode par mode, sous forme fréquentielle. Ainsi, pour chaque mode k, on réécrit l’équation (4.103) sous
la forme :
T 1
1
f tkng ( t )  L  P12 N   2 ( U, B,  k ) f te ( t )
T
(4.104)
2
 
1
où  2 U, B,  k est une matrice diagonale (3*3) contenant les racines carrées des admittances
aérodynamiques mesurées pour les forces de trainée, portance, et moment. On obtient alors chacune des
forces généralisées en multipliant à gauche par la transposée du vecteur propre correspondant :
f tk ( t )  Z k f tkng ( t )
T
(4.105)
b) Formulation des forces aéroélastiques
Considérons dans un premier temps le cas de l’élément poutre. En un point de cet élément d’abscisse ,
les déplacements verticaux et angulaires peuvent s’écrire sous la forme :
h ( , t )  y (  )  h ( t ) (4.106)
(, t )  r z ( )  ( t ) (4.107)
Le vecteur des forces aéroélastiques s’exerçant sur l’élément s’écrit alors, d’après (2.65-68) :
f aee ( t )  0 L ae ( t ) 0 0 0 M ae ( t ) 
T
(4.108)
avec :
122
L ae ( t ) 
 Lh(t)  L(t)J() d 
e
1 C C   1 C C  
2  B U   e   2
  B B   e 
BU 2  Lh1 h ( t )  Lh 2 h& ( t )   y( )J ( )d  + BU 2  L1 ( t )  L 2 & ( t )   rz ( )J( )d 


+
 Lh (t)  Lh (t)J()d   L (t)  L (t)J()d
e
3 4
e
3 4 (4.109)
M ae ( t ) 
 Mh(t)  M(t)J() d 
e
1 2 2 BC M1   1 C C  
2
B U  C M1 ( t ) 
 U   e   2
  B U   e 
& ( t )   rz () J( )d  + B 2 U 2  Mh1 h ( t )  Mh 2 h& ( t )   y( )J ( )d 


+
 Mh (t)  Mh (t) J()d   M (t)  M (t) J()d
e
3 4
e
3 4 (4.110)
où les
 Lh (t) J()d ,  L (t) J()d ,  Mh (t) J()d ,  M (t) J()d
e
k
e
k
e
k
e
k vérifient les
relations :
 
 Lh ( t ) J ( )d  =  d U Lh ( t ) J ()d  C Lhk

e
k


Lhk
B
e
k
B  y()J()d h& (t)
e
(4.111)
 
 L ( t ) J ( )d  =  d U L ( t ) J ( )d  C

e
k


Lk
B
e
k
 Lk
e

rz ( )J ( )d & ( t ) (4.112)
 
 Mh ( t ) J ( )d  =  d U Mh ( t ) J ( )d  C Mhk

e
k


Mhk
B
e
k
 B  y()J()d h& (t)
e
(4.113)
  U

 M ( t ) J ( )d  =  d

e
k


Mk
B
e
 e

M k ( t ) J ( )d  C Mk rz ()J ( )d & ( t ) (4.114)
123
On obtient ainsi la formulation vectorielle suivante, en introduisant le vecteur des déplacements d dans
un repère lié à l’élément:
f aee ( t )  C eae d&  K eae d  A e  e (4.115)
& e  D e  e  B e d& (4.116)
avec :
d(t )   d 1 (t ) h ( t ) d 3 ( t ) d 4 ( t ) d 5 (t ) (t )
T
(4.117)



 Lh (t) J() d 
e
3

  Lh (t) J() d 
e
4

  L (t) J() d 
e
3
 

e  
e
L ( t ) J () d 
4
 (4.118)

 e
Mh ( t ) J ( ) d 
3


 e
Mh ( t ) J ( )d 
4

 


 M (t) J() d 
e
3

  M (t) J() d 
e
4
0 0 0 0 0 0 
 C Lh 2 C L 2 B 
0

U e 
y (  ) J (  ) d 0 0 0
U e 
rz ()J ( )d 

1 2 0 0 0 0 0 0
C ae  BU 
e  (4.119)
2 0 0 0 0 0 0 
0 0 0 0 0 0 
 C M 2 B 2 
C B
0

Mh 2
U e 
y (  ) J (  ) d 0 0 0
U erz ( )J ( )d 

0 0 0 0 0 0 
 C Lh1 
0

B e 
y (  ) J (  ) d 0 0 0

C L1 rz ( )J ( )d 
e

1 2 0 0 0 0 0 0
K eae  BU   (4.120)
2 0 0 0 0 0 0 
0 0 0 0 0 0 
 
 e 
0 C Mh1 y( )J ( )d 0 0 0

C M1 B rz ()J ( )d 
e 
124
0 0 0 0 0 0 0 0
1 1 1 1 0 0 0 0

0 0 0 0 0 0 0 0
Ae    (4.121)
0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0
 
0 0 0 0 1 1 1 1
d Lh 3 0 0 0 0 0 0 0 
 0 d Lh 4 0 0 0 0 0 0 

 0 0 d L 3 0 0 0 0 0 
 
U 0 0 0 d L 4 0 0 0 0 
D 
e
(4.122)
B 0 0 0 0 d Mh 3 0 0 0 
 
 0 0 0 0 0 d Mh 4 0 0 
 0 0 0 0 0 0 d M 3 0 
 
 0 0 0 0 0 0 0 d M 4 
 C Lh 3 
0

B  y (  ) J (  ) d
e
0 0 0 0 

C Lh 4
0
 B  y() J() d
e
0 0 0 0 

0

0 0 0 0 C L 3

e
rz ( ) J ( ) d 

 
0
Be  
0 0 0 0 C L 4

e
rz ( ) J () d 
 (4.123)
C Mh 3
0
 B e
y (  ) J (  ) d 0 0 0 0 

0 C Mh 4 


B  y() J() d
e
0 0 0 0 

0

0 0 0 0 C M 3

e
rz ( ) J ( ) d 

0

0 0 0 0 C M 4

e
rz ( ) J ( ) d 

On peut maintenant exprimer nos équations dans la base nodale de l’élément et le repère général en
introduisant les matrices bien connues N, P12 et en intégrant par la méthode du point milieu. d ng devient
le vecteur colonne (12*1) contenant les déplacements aux deux nœuds de l’élément dans le repère
général,  e est inchangé :
f aeng ( t )  C ng & ng  K ng d ng  A ng  e
ae d ae (4.124)
& e  D  e  B d& ng
ng ng
(4.125)
125
avec :
  0 0 0 0 0 0   
   
  C Lh 2 C L 2 B  
 0 0 0 0  
  
U U
  
 1 1
T
2 0 0 0 0 0 0   1 
 P12  L  N   BU  
 N   P12
T
C aeng (4.126)
2 2 0 0 0 0 0 0 
     2 
  0 0 0 0 0 0   
  0 C Mh 2 B 2
C M 2 B   
  0 0 0  
   U U   
  0 0 0 0 0 0   
   C Lh1   
  0 0 0 0 C L1   
 T B  
1 1 2 0 0 
T
 P12  L  N  BU 
0 0 0 0
  N 1   P
  2   12
K aeng (4.127)
2 2 0 0 0 0 0 0 
     
   
 0 0 0 0 0 0 

 0 C Mh1 
  0 0 0 C M1 B  
T
1
A ng  P12 N  A e (4.128)
2
 C Lh 3 
0 B
0 0 0 0 
 C Lh 4 
0 0 0 0 0 
 B 
0 0 0 0 0 C L 3 
0 0 0 0 0 C L 4   1  T
B ng L N  P12 (4.129)
0 C Mh 3
0 0 0 0  2
 B 
 C Mh 4 
0 0 0 0 0 
 B 
0 0 0 0 0 C M 3 
0 0 0 0 0 C M 4 

D ng  D e (4.130)
126
Il nous reste maintenant à assembler les matrices élémentaires. Les matrices C ae et K ae sont assemblées
suivant la procédure habituelle. Les matrices A g et B g , de dimensions respectives (ndlt, 8*ne) et (8* ne,
ndlt), sont assemblées suivant la procédure décrite figures (4.7) et (4.8). La matrice, de dimension (8* ne,
8* ne), est composée des matrices élémentaires D ng positionnées sur sa diagonale principale.
Matrice Ang (12, 8) de l’élément k : Matrice globale A (ndlt, 8*nel) :
élt 1 élt k élt nel

nœud i noeud 1
nœud i
nœud j nœud j
nœud nnd
Figure 4.7 – Procédure d’assemblage de la matrice A
Matrice Bng (8,12) de l’élément k : Matrice globale B (8*nel, ndlt) :
nœud i
nœud j
élt 1
nœud i elt k
elt nel
nœud j
Figure 4.8 – Procédure d’assemblage de la matrice B
On obtient alors la formulation matricielle suivante :
Fae ( t )  C ae d&  K ae d  A  (4.131)
&  D   B d& (4.132)
127
c) Formulation simplifiée des forces aéroélastiques
La formulation précédente n’a en fait d’intérêt que lorsque les fréquences de vibration des différents
modes propres de l’ouvrage sont proches ou lorsque la résolution des équations de la dynamique sont
effectuées dans la base nodale. En outre, si les coefficients des fonctions indicielles sont mal estimés, on
peut assister à des problèmes d’instabilités numériques. Dans le cas où les fréquences des modes sont
bien détachées et la résolution effectuées dans la base modale, on peut se contenter de la détermination
d’une matrice raideur et d’un matrice amortissement aérodynamiques diagonales similaires à celles
définies en (4.51), à ceci près que chaque terme est associé à sa pulsation modale. C’est cette méthode
que nous utilisons dans le cadre de notre procédure de validation.
d) Mise en équation
Les équations du mouvement s’écrivent donc, dans la base nodale, en tenant compte des forces
turbulentes et aéroélastiques :
M &d&( t )   C  C ae  d& ( t )   K  K ae  d ( t )  A   Ft ( t ) (4.133)

   
& ( t )  D  ( t )  B d& ( t ) (4.134)
et, dans la base modale :
M̂ &q&( t )   Ĉ  Ĉ ae  q& ( t )   K̂  K̂ ae  q ( t )   A    Ft ( t )
T T
(4.135)
   
& ( t )  D  ( t )  B  q& ( t ) (4.136)
Dans le cas d’une formulation simplifiée des forces aéroélastiques, on écrit simplement :
M̂ &q&( t )   Ĉ  Ĉ ae  q& ( t )   K̂  K̂ ae  q ( t )   Ft ( t )
T
(4.137)
   
où Ĉ ae et K̂ ae sont les matrices amortissement et raideur aérodynamiques définies en c). Elles n’ont donc
rien à voir avec celles intervenant en (4.135).
128
4.5.2 Résolution du système dans la base modale
L’intégration du système d’équations différentielles (4.135-136) est effectué à l’aide d’un schéma
temporel implicite de Newmark pour les termes de déplacements et un schéma explicite pour le terme .
Ces schémas permettent d’écrire les équations reliant les valeurs des variables à l’instant t aux valeurs
des variables à l’instant t+t :
t 2
q ( t  t )  q( t )  t q& ( t )   &q&(t )  &q&( t  t ) (4.138)
4
t
q& ( t  t )  q& ( t )   &q&( t )  &q&(t  t ) (4.139)
2
 ( t  t )   ( t )  t & ( t ) (4.140)
En injectant ces équations dans le système d’équations (4.135-136), on trouve :
 4 2   4 2 
 2 M̂   Ĉ  Ĉ ae    K̂  K̂ ae   q ( t  t )   2 M̂   Ĉ  Ĉ ae   q ( t ) 
 t t       t t   
 4
 
 T T

 M̂   Ĉ  Ĉ ae   q& ( t )  M̂ &q&( t )   Ft ( t  t )   A  ( t )  t D  ( t )  t B  q& ( t )  (4.141)
 t
Une fois calculé  ( t  t ) , on calcule q ( t  t ) , puis &q&( t  t ) à partir de (4.138), et enfin q& ( t  t ) à
partir de (4.139) . On réitère cette procédure le nombre de pas de temps désiré.
Dans le cas où la formulation simplifiée des forces aéroélastiques est employée, la procédure
d’intégration temporelle est strictement la même, à ceci près que les termes A, B, D,  n’interviennent
plus.
129
4.6 Introduction de simulations temporelles dans le calcul spectral
L’objectif est d’introduire des procédures de calcul temporel dans le module de calcul spectral afin
d’améliorer les performances de ce dernier. Bien entendu, ces apports ne restent intéressants que dans la
mesure où ils n’occasionnent pas de surcoûts significatifs en terme de temps de calcul. Deux procédures
sont ainsi introduites. La première permet d’optimiser la linéarisation des courbes d’évolution des
coefficients stationnaires en fonction de l’incidence du vent. La seconde permet de combiner les
déplacements modaux extrêmes afin d’en déduire les déplacements nodaux extrêmes.
4.6.1 Méthode de « linéarisation équivalente » des coefficients aérodynamiques
On a vu que l’une des faiblesses majeures de l’approche spectrale était l’obligation de représenter
l’évolution des coefficients aérodynamiques en fonction de l’incidence du vent par des droites affines.
Nous proposons ici une méthode permettant d’optimiser cette linéarisation Le but est de trouver une
expression linéaire de l’évolution des coefficients aérodynamiques équivalente à l’évolution réelle des
coefficients en terme de densité spectrale des déplacements résultante.
Supposons qu’on s’intéresse au calcul spectral des déplacements d’un tablier de pont projetés sur un
mode k caractérisé par sa masse généralisée Mk, sa pulsation k, et son amortissement k. On considère
alors le cas d’une maquette sectionnelle virtuelle à un degré de liberté et de caractéristiques mécaniques
équivalentes.
 u(t ) 
On simule alors une trajectoire de vent unidimensionnelle   respectant une distribution spectrale
 w(t) 
appropriée. A partir de ce signal, on peut construire un signal de forces (portance, traînée, ou moment).
 u(t ) 
Pour cela, il suffit, pour chaque série   , de calculer l’angle d’incidence instantané du vent, puis de
 w(t) 
se reporter à la courbe d’évolution de la force en question en fonction de l’angle d’incidence. On obtient
ainsi un signal de forces réaliste car non linéarisé. Ce signal n’est donc a priori plus gaussien !
On peut alors calculer la densité spectrale de puissance du signal de forces, puis, à l’aide de la fonction
de transfert simplifiée (4.142), la variance des déplacements résultants :
S Ftemp. ( k )
k 
2
(4.142)
2 k  k
3
130
Il nous reste maintenant à trouver une formulation linéaire de l’évolution des forces aérodynamiques en
fonction de l’incidence du vent c’est-à-dire un couple ,   tel que la relation suivante :
1
F( t )  B( Uu ( t )  Vv( t ))  B( Uv( t )  Vu ( t )) (4.143)
2
fournisse un signal de forces générant le même écart type de déplacements. En fait, on voit bien que
seule la valeur SF(k) a besoin d’être reproduite.
Si on écrit l’expression de la densité spectrales des forces à partir de (4.143) :
S F (  k )  x S u ( k )  y S v (  k ) (4.144)
avec
2
  1 
x  B  U   V   (4.145)
  2 
2
  1 
y  B V   U  (4.146)
  2 
Il reste à estimer par une méthode d’optimisation linéaire les valeurs de x et y, dont on déduit  et  :
U x V y
 (4.147)
BU 2  V 2 ) 
U y V x
 (4.148)
1 2 BU 2  V 2 ) 
Les figures (4.9) et (4.10) montrent les résultats obtenus pour les tabliers trapézoïdal et rectangulaire.
L'incidence de vent moyenne est de 3°. A gauche, figurent les signaux temporels aléatoires simulés et
leurs spectres respectifs. A droite, sont représentés les courbes d'évolution des coefficients
aérodynamiques en fonction de l'incidence. Les courbes noires correspondent aux coefficients réels,
tandis que les courbes en couleur représentent les coefficients linéarisés.
Les résultats sont logiquement très dépendants du type de vent et de la gamme de fréquences utilisée
pour l'optimisation. Aussi il convient de se placer dans les mêmes conditions que celles qui seront
imposées dans le calcul spectral. Dans le cas présent, les caractéristiques du vent sont de 8%, 5%, 0,5 m,
0,1 m, respectivement pour les écarts types et les échelles de la turbulence. La gamme de fréquence est
de 5 à 10 Hz
131
Précisons que pour chaque couple C k , dC k d le nombre de solutions est infini, y compris les plus
saugrenues d'un point de vue physique. On s'arrange pour ne conserver que les solutions physiquement
acceptables.
−3
x 10
12
6 0.1
U(m/s)
10 4
0.08
2
8 0.06
5 10 15 20 25 30 35 40 2 5 10 20
Cx
t(sec.)
−4
f(Hz) 0.04
x 10
1.5 8
0.02
1
W(m/s)
6
0.5 0
0 4
−0.5
2 −6 −4 −2 0 2 4 6
5 10 15 20 25 30 35 40 2 5 10 20
Incidence (deg.)
t(sec.) f(Hz)
10 100
Inc.(deg.)
5
50
0 0
0
5 10 15 20 25 30 35 40 2 5 10 20
Cy
t(sec.) f(Hz)
−6
x 10
−0.5
0.15 3
0.1 2
Cx
0.05 1
−1
0
5 10 15 20 25 30 35 40 2 5 10 20 −6 −4 −2 0 2 4 6
t(sec.) f(Hz) Incidence (deg.)
−4
x 10
0.5
3
0 2 0.15
Cy
−0.5 1
0.1
0
5 10 15 20 25 30 35 40 2 5 10 20
t(sec.) f(Hz) 0.05
Cm
−6
x 10
8 0
0.15
6
0.1
−0.05
Cm
0.05 4
0 2
−0.05
0
−0.1
5 10 15 20 25 30 35 40 2 5 10 20
−6 −4 −2 0 2 4 6
t(sec.) f(Hz)
Incidence (deg.)
Figure 4.9 – Linéarisation des coefficients aérodynamiques du tablier trapézoïdal

pour une incidence de vent moyenne de 3°
−3
x 10
8
12 0.15
6
U(m/s)
10 4 0.14
2
8
0.13
Cx
5 10 15 20 25 30 35 40 1 2 5 10 20
t(sec.) f(Hz)
x 10
−4 0.12
1.5
8 0.11
1
W(m/s)
0.5 6
0 4
0.1
−0.5
2 −6 −4 −2 0 2 4 6
5 10 15 20 25 30 35 40 1 2 5 10 20
Incidence (deg.)
t(sec.) f(Hz)
10
100
Inc.(deg.)
5 0.6
50
0
0.4
−5 0
5 10 15 20 25 30 35 40 1 2 5 10 20 0.2
Cy
t(sec.) f(Hz)
x 10
−6
0
6
0.2 −0.2
4
Cx
0.15
2
−0.4
0.1
0
5 10 15 20 25 30 35 40 1 2 5 10 20 −6 −4 −2 0 2 4 6
t(sec.) f(Hz) Incidence (deg.)
−4
x 10
3
0.5
2
0.02
Cy
0
1
−0.5 0
5 10 15 20 25 30 35 40 1 2 5 10 20 0.01
t(sec.) f(Hz)
Cm
−7
x 10
0.02 0
0.01 2
Cm
0
1
−0.01
−0.01
0
5 10 15 20 25 30 35 40 1 2 5 10 20
−6 −4 −2 0 2 4 6
t(sec.) f(Hz)
Incidence (deg.)
Figure 4.10 – Linéarisation des coefficients aérodynamiques du tablier rectangulaire

pour une incidence de vent moyenne de 3°
132
4.6.2 Combinaison des charges modales extrêmes en cas de charges statiques équivalents
La stricte connaissance des déplacements modaux extrêmes induits par le vent ne présente un intérêt que
dans l’optique d’une évaluation du confort des usagers de l’ouvrage étudié. Ce qui, dans l’absolu,
intéresse les bureaux d’études chargés de dimensionner les structures, c’est de connaître les contraintes
extrêmes sollicitant les matériaux occasionnées par les vibrations de l’ouvrage. Dans la pratique, il est
plus commode pour ces bureaux d’études de travailler avec des cas de charges statiques équivalents qui
seront combinés avec les autres cas de charges : poids propre, exploitation, neige, … dans le cadre de
dimensionnements ELU ou ELS.
C’est dans cette optique que nous proposons une méthode permettant, à partir de la connaissance des
densités spectrales de puissance des déplacements modaux issus de calculs spectraux ou temporels, de
déterminer une série de charges statiques fictives qui seraient un reflet fidèle des charges dynamiques
réelles.
La validité d’un certain nombre d’outils de R.D.M. utilisés par la suite est théoriquement limitée à des
études statiques. On suppose néanmoins que les fréquences de vibration propres aux ouvrages d’arts tels
que les ponts suspendus ou à haubans (0,1-0,5 Hz) sont suffisamment basses pour que les approximations
effectuées ne soient pas trop grossières.
On introduit les vecteurs E , G , S , I X , I Y , I Z contenant respectivement les modules d’Young,
modules de cisaillement, sections, inertie suivant l’axe x, inertie suivant l’axe y et inertie suivant l’axe z
en chaque nœud de la structure discrétisée. On introduit également les coordonnées curvilignes sk
associées à chacune des parties clés composant l’ouvrage (tablier, pile,…) ; on définit alors, pour chaque
mode propre considéré, le vecteur Z i ' (i=1,n) qui est la dérivée du vecteur propre Z i par rapport aux
coordonnées curvilignes définissant la structure (fig. 4.11).
s1
s2
s3
Figure 4.11 – Définition des coordonnées curvilignes
133
Pour chaque élément e de la structure on définit le vecteur élémentaire Z ie ' par la relation suivante :
Z ie1 PL  Z ie 2 PL
T T
Z '
e
i (4.149)
L
Z ie1 et Z ie 2 sont des vecteurs colonnes de dimension 6 extraits du vecteur Z i et contenant les déformées
des deux nœuds de l’élément e ; PL est la matrice de changement de repère local général et L la longueur
de l’élément e. Le vecteur Z i ' est obtenu par assemblage des Z ie ' . Contrairement au vecteur Z i , il n’est
pas défini au niveau des nœuds de la structure mais au niveau du milieu de ses éléments (croix fig. 4.11).
La raison est double : d’abord, les caractéristiques mécaniques de la structure discrétisée sont définies
élément par élément ; ensuite, la définition d’un repère local est immédiate.
A partir de Z i ' , E , G , S , I X , I Y , I Z , on peut calculer les vecteurs N i , M Xi , M Yi , M Zi , TXi , TYi ,
de dimension ne, contenant respectivement les efforts normaux, moments de flexion suivant x, moments
de flexion suivant y, moments de torsion suivant z et efforts tranchants en chaque élément pour un
déplacement généralisé i unitaire. Les vecteurs N i , M Xi , M Yi , M Zi , TXi , TYi , contiennent des termes
de la forme :
N i (e)  E (e)  Z ei ' 6  (e  1)  3 (4.150)
TXi (e)  G (e) S(e)  Z ie ' 6  (e  1)  1 (4.151)
TYi (e)  G (e) S(e)  Z ei ' 6  (e  1)  2 (4.152)
M Xi (e)  E (e) I X (e)  Z ie ' 6  (e  1)  4  (4.153)
M Yi (e)  E(e) I Y (e)  Z ie ' 6  (e  1)  5 (4.154)
M Zi (e)  G (e) I Z (e)  Z ie ' 6  (e  1)  6  (4.155)
Remarquons que les efforts tranchants sont définis à un facteur près, fonction de la géométrie de la
section. Cela n’a aucune importance dans la mesure où notre but n’est pas de calculer les efforts en
chaque point de la structure mais uniquement de les « maximiser ».
134
On cherche alors à construire des signaux temporels réalistes N i (t), M Yi (t), M Zi (t), M Xi (t), TXi , TYi ,
représentant les fluctuations des différents efforts sous les sollicitations du vent. Pour cela, on
uniquement besoin de signaux temporels réalistes des déplacements modaux q i t  .
Dans le cas où un calcul temporel a été effectué, il suffit de réutiliser les signaux enregistrés. Dans le cas
où seul un calcul spectral a été effectué, il s’agit de reproduire artificiellement le comportement
dynamique de l’ouvrage à partir de la connaissance des écarts types, pulsation et amortissement de
chaque mode de vibration. On génère donc les signaux aléatoire q i t  à partir de la fonction densité
spectrale de puissance d’un oscillateur linéaire à 1 degré de liberté (figure 4.12) :
2 i  i  i 
3 2
S i ()  (4.156)
 i
2
 2 
2
 4 i  2  i
2 2
Fréquence: 5Hz , amortissement: 0.5%

3
−1
−2
−3
0 2 4 6 8 10 12 14 16 18 20
t(sec.)
1.5
0.5
0
0 5 10 15 20 25
f(Hz)
Figure 4.12 - Signal aléatoire respectant la fonction densité spectrale de puissance

d’un oscillateur linéaire à un degré de liberté amorti
On combine alors les N i , M Yi , M Zi , M Xi , TXi , TYi , avec les q i ( t ) pour obtenir des signaux
temporels N (t), M Y (t), M Z (t), M X (t), TX ( t ) , TY ( t ) :
135
N(t )   q (t) N
i 1
i i (4.157)
n
M X (t)   q (t ) M
i 1
i Xi (4.158)
n
M Y (t)   q (t ) M
i 1
i Yi (4.159)
n
M Z (t)   q (t ) M
i 1
i Zi (4.160)
n
TX ( t )   q (t ) T
i 1
i Xi (4.161)
n
TY ( t )   q (t ) T
i 1
i Yi (4.162)
Généralement, le dimensionnement des ouvrages d’art est effectué sur la base d’une vitesse de vent
cinquantenale, correspondant à la valeur, moyennée sur 10 minutes, susceptible d’être atteinte sur une
période de 1 an avec une probabilité égale à 2 %. On effectue donc un nombre nt de simulations d’une
durée de 10 minutes.
Pour chaque simulation s, on identifie les valeurs maximales prises par chacun des efforts F (F = N, MX,
MY, MZ, TX, TY), en un nombre nec d’éléments clés ec et les instants t eFsc correspondants. On moyenne
alors les valeurs maximales obtenues sur les nt simulations. On construit alors des vecteurs « cas de
charges extrêmes »  eNc ,  eTcX ,  eTcY ,  eMc X ,  eMc Y ,  eMc Z , par combinaison linéaire des forces d’inerties
modales à des instants t eFsc correspondant à des valeurs extrêmes proches des moyennes calculées:
 eNc  m ddl   q (t
i 1
i
ec
Ns ) Z i  i2 (4.163)
n
 eTcX  m ddl   q (t
i 1
i
ec
TX s ) Z i  i2 (4.164)
n
 eTcY  m ddl   q (t
i 1
i
ec
TY s ) Z i  i2 (4.165)
n
 ec
MX  m ddl   q (t
i 1
i
ec
M Xs ) Z i  i2 (4.166)
n
 eMc Y  m ddl   q (t
i 1
i
ec
M Ys ) Z i  i2 (4.167)
n
 eMc Z  m ddl   q (t
i 1
i
ec
M Zs ) Z i  i2 (4.168)
136
où m ddl représente le vecteur colonne contenant les masses et inerties massiques ponctuelles affectées à
chaque degré de liberté et  symbolise la multiplication vectorielle terme à terme.
On obtient donc, en théorie, 6*nec cas de charges équivalents. Dans la pratique, le nombre de cas de
charges associés à chaque élément clé sera inférieur à 6 ; par exemple, l’effort tranchant ne présente
d’intérêt qu’au niveau des liaisons ; l’effort normal pourra lui être supposé marginal à côté des efforts de
flexion … Ces cas de charges peuvent alors être intégrés, un par un, dans le cadre de dimensionnements
ELU ou ELS avec les coefficients de sécurité appropriés.
4.7 Conclusion
Ce chapitre constitue le cœur de notre travail. Nous y avons exposé les principes fondamentaux
permettant de développer des codes de calcul basés sur les approches spectrale et temporelle. Il convient
de rappeler que le calcul numérique ne constitue que le maillon central dans le dimensionnement au vent
d’une structure élancée.
Le premier maillon, consistant à caractériser expérimentalement les propriétés aérodynamiques de

l’ouvrage, a fait l’objet du troisième chapitre. Il pourrait néanmoins, à terme, être intégré à un code de
calcul temporel sous forme de module C.F.D. résolvant les équations de Navier Stokes pas de temps par
pas de temps. Cela n’est cependant pas encore envisageable avant plusieurs années.
Le dernier maillon concerne la détermination de cas de charges statiques équivalents utilisables par les
bureaux d’études, à partir de la connaissance des réponses modales. Nous avons, dans ce chapitre,
proposé une méthode basée sur la simulation temporelle des efforts dynamiques s’exerçant au sein de la
structure. Les cas de charges ainsi déterminés peuvent alors être affectés de coefficients de sécurité, et
combinés séparément aux autres cas de charges classiques (exploitation, neige, …) dans le cadre de
dimensionnements ELS ou ELU.
137
Chapitre 5 Validation
expérimentale
5 Validation expérimentale
Après avoir présenté les principes théoriques utilisés pour l’élaboration de notre code de calcul, il s’agit
d’en vérifier la pertinence. Nous avons donc mis en place une procédure de validation expérimentale
basée sur le passage en soufflerie atmosphérique de deux maquettes aéroélastiques représentant des
travées centrales de tabliers de ponts. Les deux profils de tabliers étudiés sont ceux déjà présentés
chapitre 3 : un profil trapézoïdal et un profil rectangulaire de ratio largeur/épaisseur égal à 8.
Le premier est représentatif de la famille des tabliers « profilés » : vitesse critique de flottement élevée et
évolutions des coefficients aérodynamiques en fonction de l’incidence relativement linéaire. Le second,
bien qu’atypique du point de vue géométrique, représente bien le comportement aérodynamique d’un
tablier mal profilé : vitesse critique faible et coefficients aérodynamiques fortement non linéaires.
Deux configurations de vent sont étudiées : turbulente et uniforme. La première est a priori favorable à la
prévision de la réponse aux buffeting forces mais moins fiable en ce qui concerne l’appréhension des
forces aéroélastiques ; au contraire, la seconde doit théoriquement bien se prêter à la prédiction du
comportement instable de la structure mais être plus dure à modéliser en matière de « buffeting forces »
en raison du rôle non négligeable du sillage.
Le principe de la procédure consiste, on l’aura deviné, à mesurer la réponse de nos deux maquettes au
vent de la soufflerie, effectuer les calculs correspondants avec les approches spectrale et temporelle, et
comparer les résultats obtenus.
138
Chapitre 5 Validation expérimentale
5.1 Présentation des maquettes
L’objectif était de fabriquer des maquettes représentant le comportement dynamique de travées centrales
« fictives » de tabliers de ponts. L’utilisation de la technologie « taut-tube » était envisageable mais ne
fut pas retenue. En effet, l’intérêt majeur de ce genre de maquettes est de pouvoir régler avec précision
les fréquences de vibrations des deux premiers modes de flexion et de torsion. Cela est fort appréciable,
dans le cadre d’une étude, lorsqu’il s’agit de reproduire le comportement dynamique d’une structure
réelle : les fréquences de vibrations à échelle réduite sont alors imposées.
Cela n’étant pas notre cas, nous avons opté pour une technologie moins sophistiquée, mais beaucoup plus
facile à réaliser et économique, basée sur l’utilisation d’un profilé standard du commerce. Un avantage
non négligeable lié à ce procédé est de pouvoir reproduire des modes de torsion excentrés, plus
représentatifs du comportement mécanique réel des tabliers de pont que les modes habituellement
représentés par les maquettes de type taut-tube.
Notre cahier des charges était le suivant :
 Echelle de la maquette comprise entre le 100ème et le 200ème pour s’accorder aux échelles des
structures tourbillonnaires reproductibles par la soufflerie en configuration turbulente,
 Fréquences de vibration des premiers modes comprises entre 2 et 10 Hz pour que la gamme de
vitesses soufflerie (0 - 20 m/s) corresponde à une gamme de vitesses réduites raisonnable, compte-
tenu de l’échelle de la maquette,
 Premiers modes de flexion verticale et de torsion aussi proches que possible en fréquence, pour
éviter que l’un des deux soit marginalisé par l’autre,
 Masses généralisées et amortissements modaux réalistes par rapport à la vraie grandeur,
 Flèche statique minimale,
 Longueur maximale.
Un profilé ouvert était nécessaire pour garantir une fréquence en torsion raisonnable. Le U s’est imposé
face au T pour des raisons de facilités de montage des éléments. Un U « plat » aurait été idéal pour
rapprocher au maximum les modes de torsion et de flexion. Cependant, pour minimiser la flèche statique,
nous avons opté pour un U « carré ».
Les deux maquettes sont donc basées sur le même principe structurel : un profilé U 30*30*3 mm en
aluminium de longueur totale 295 cm encastré de chaque coté sur une paire de poteaux de sections
139
carrées en aluminium de 30*30 mm hauts de 50 cm. La longueur libre entre les deux encastrements est
égale à 249 mm. Chaque paire de poteaux est fixée au sol en bois de la soufflerie par 6 vis via une platine
en aluminium de 3 cm d’épaisseur. Des nervures et des tirants assurent les rigidités longitudinales et
transversales des poteaux.
Quinze éléments de tabliers à l’échelle du centième sont fixés sur le profilé U. Seuls onze d’entre eux
participent réellement au comportement dynamique de la maquette ; les quatre autres sont fixes et ont
pour unique fonction d’assurer la continuité aérodynamique et minimiser les effets de bouts. Ces
éléments, larges de 19,5 cm et espacés de 1 mm, sont en résine polyuréthane de densité 0,5. Outre sa
densité, cette résine présente l’avantage de garantir un bon aspect de surface et une bonne usinabilité. La
liaison avec le profilé est assurée par une paire de vis et des cales en aluminium de 2 cm de large posées
sur les trois faces extérieures du U. On peut donc considérer cette liaison comme quasi linéaire et
interférant peu sur la raideur de l’ensemble.
Les éléments de résine, fabriqués « ouverts » sont refermés sur leur face inférieure par des plaques
d’aluminium. Pour rapprocher d’avantage les modes de torsion et de flexion, des barres d’acier ont été
fixées sur la face supérieure des plaques d’aluminium, à l’intérieur du U, et des trous circulaires ont été
percés aux extrémités du profil rectangulaire.
Figure 5.1 – Vues de dessus et face de l’ossature de la maquette aéroélastique

avec deux éléments de tablier
140
Figure 5.2 – Vues de dessus, face et coté d’un pied de maquette et coupes sur les éléments de tabliers
Figure 5.3 – Eléments de maquettes aéroélastiques pendant la phase de fabrication
141
5.2 Caractérisation mécanique des maquettes
Pour caractériser les modes propres de la maquette, on a opté pour une double approche numérique/
expérimentale. L’intérêt était surtout de comparer les résultats obtenus par les deux méthodes, les
mesures faisant office de référence.
On s’est arrangé pour que les masses et inerties linéiques des deux profils soient aussi proches que
possible. Ainsi, le calcul des modes propres, mesure des déformées et mesure des masses généralisées
sont uniquement effectués pour le profil rectangulaire, plus facile à manipuler. Seuls les fréquences et
amortissements sont mesurés dans les deux cas. Supposant les déformées des deux maquettes identiques,
on estime les masses généralisées du tablier trapézoïdal à partir de celles du tablier rectangulaire et des
fréquences mesurées.
5.2.1 Mesure des fréquences, amortissements et masses généralisées
Les signaux de déplacements sont capturés à l’aide de capteurs de déplacements à triangulation laser. Les
mesures des fréquences et amortissements modaux sont issus de traitements des signaux dans le domaine
temporel. On a ainsi pu caractériser l’évolution de ceux-ci en fonction de l’amplitude de vibration. Celle-
ci s’avère d’ailleurs fortement linéaire. Les résultats sont représentés figure 5.15.
L’estimation des masses généralisées est effectuée par ajouts successifs de masses aux points de
déformées maximales et mesure des fréquences d’oscillations correspondantes (fig. 5.14). La justesse des
résultats résulte d’un bon compromis dans le dosage de ces masses : suffisamment élevées pour que les
mesures soient précises mais assez discrètes pour ne pas modifier l’allure des déformées modales.
5.2.2 Mesure des déformées modales
Le tablier est excité verticalement à une de ses extrémités par un moteur linéaire alimenté par un
générateur de fréquence amplifié. Chaque mode est passé en revue par simple ajustement de la fréquence
d’excitation. Les mesures sont effectuées par des capteurs de déplacements à triangulation laser sur la
partie opposée du tablier. Un capteur fixe, mesurant les déplacements verticaux en quart de travée,
fournit les signaux de référence. On s’affranchit ainsi des problèmes de variation d’intensité de
l’excitation au cours du temps. Deux autres capteurs, montés sur un support mobile, mesurent
respectivement les déplacements horizontaux et verticaux (figures 5.4-5.6). Les amplitudes des
déformées sont alors déduites des écarts types des déplacements mesurés rapportés à la mesure de
référence.
142
Figures 5.4, 5.5, 5.6 – Dispositif de mesure des déformées modales
143
5.2.3 Calcul des déformées
Dans un premier temps, le calcul des modes propres a été effectué à partir d’un logiciel commercial
éléments finis, en modélisant le profilé U par une série d’éléments poutres. Les trois premiers modes
obtenus étaient alors :
 un mode de pure torsion aux alentours de 5,8 Hz ;

 un mode de pure flexion verticale vers 7,9 Hz ;
 un mode de pure flexion horizontale à 10,0 Hz.
Des mesures effectuées sur la vraie maquette ont pourtant vite montré que les déformées des modes réels
de la structure étaient plus complexes que ceux calculés. Ainsi, le premier et le troisième mode
ressemblaient plutôt à deux modes de torsion dont les centres étaient situés respectivement au-dessus et
au-dessous du tablier. Le problème venait probablement du fait que le calcul ne tient pas compte de
l’excentrement du centre de cisaillement du profil U par rapport à son centre d’inertie.
On a donc modélisé le profil U par 15*3 plaques planes solidaires. Les résultats (fig. 5.7-5.12 et 5.13) se
sont avérés nettement meilleurs bien que non exactement conformes aux mesures, en particulier en ce qui
concerne les composantes horizontales des déformées des modes de torsion. On a quand même pu obtenir
des informations non accessibles par la mesure comme, par exemple, les composantes en torsion autour
de l’axe x du mode de flexion verticale.
5.2.4 Récapitulatif des fréquences, amortissements et masses généralisées utilisés pour le

calcul numérique
On rappelle que les masses généralisées du tablier trapézoïdal n’ont pas été directement mesurées mais
déduites de celles du tablier rectangulaire. Les valeurs des fréquences propres et amortissements
(rapportés au critique) correspondent à des oscillations de faibles amplitudes (0,5-2mm).
Tablier rectangulaire Tablier trapézoïdal

F(Hz) (%) Mg (kg.m²) F(Hz) (%) Mg (kg.m²)
Mode de torsion 1 5,79 0,25 1,76 5,5 0,25 1,95
Mode de flexion verticale 7,94 0,20 7,69 7,41 0,20 8,82
Mode de torsion 1bis 10,08 0,35 2,18 9,57 0,35 2,42
Mode de torsion 2 12,25 0,20 1,51 11,63 0,20 1,67
Tableau 5.1 – Fréquences propres, amortissements et masses généralisées utilisées pour le calcul
144
Figures 5.7 et 5.8 – Déformées des modes de vibration 1 et 2 obtenues par le calcul
145
Déformée mode T1 Déformée mode T2
1,0 1,0
0,8 0,8
Verti. exp.
0,6 Verti. exp.
0,6 Hori. exp.
0,4 Hori. exp.
0,4 Verti. num.
0,2 Verti. num .
0,2 Hori. num.
0,0 Hori. num .
0,0 -0,2 0 1 2 3 4 5 6
0 2 4 6 -0,4
Déformée mode V1 Déformée mode T3
1,0 1,0
0,8
0,5 Verti. exp.
0,6 Verti. exp.
Hori. exp.
0,4 Verti. num. 0,0
Verti. num .
0,2 0 1 2 3 4 5 6
-0,5 Hori. num .
0,0
0 2 4 6 -1,0
Déformée mode T1b Déformée mode V2
1,0 1,0
0,8 Verti. exp. 0,8
0,6 Hori. exp. 0,6 V erti. exp.
0,4 Verti. num. 0,4 V erti. num.
0,2 Hori. num. 0,2
0,0 0,0
0 2 4 6 0 2 4 6
Figure 5.13 – Comparaison des déformées modales expérimentales et numériques

du tablier rectangulaire(moitié de travée)
Mode de torsion 1 Mode de torsion 2
2
2
M(kg.m²)
M(kg.m²)
1 1
0 0
0,0 0,2 0,4 0,6 0,0 0,2 0,4 0,6
dm (kg) dm (kg)
Mode de flexion 1 Mode de torsion 3
10 2
M(kg.m²)
M(kg.m²)
5 1
0 0
0,0 0,2 0,4 0,6 0,8 0,0 0,2 0,4 0,6
dm (kg) dm (kg)
Mode de torsion 1bis Mode de flexion 2
3 20
M(kg.m²)
M(kg.m²)
15
2
10
1 5
0 0
0,0 0,2 0,4 0,6 0,8 0,0 0,2 0,4 0,6
dm (kg) dm (kg)
Figure 5.14 - Détermination des masses généralisées du tablier rectangulaire
146
Amortissement torsion 1 Fréquence torsion 1
0,7 5,84
Amortis s ement/crit.
0,6 5,82
0,5 5,80
5,78
0,4
f(Hz)
5,76
0,3
5,74
0,2 5,72
0,1 5,70
0 5,68
0 1 2 3 4 5 6 7 8 0 1 2 3 4 5 6 7 8
Amplitude (mm) Amplitude (mm)
Amortissement flexion 1 Fréquence flexion 1
0,8 8,00
0,7
0,6 7,95
0,5
f(Hz)
0,4 7,90
0,3
0,2 7,85
0,1
0 7,80
0 1 2 3 4 5 6 0 1 2 3 4 5 6
Am ortissement torsion 1bis Fréquence torsion 1bis
1,2 10,20
10,15
1,0
10,10
0,8 10,05
f(Hz)
0,6 10,00
0,4 9,95
9,90
0,2 9,85
0,0 9,80
0 0,5 1 1,5 2 2,5 0 0,5 1 1,5 2 2,5
0,6 12,35
0,5
12,30
0,4
12,25
f(Hz)
0,3
12,20
0,2
0,1 12,15
0 12,10
0 0,5 1 1,5 2 0 0,5 1 1,5 2
1 19,30
0,8 19,20
19,10
0,6
f(Hz)
19,00
0,4
18,90
0,2 18,80
0 18,70
0 0,2 0,4 0,6 0,8 1 1,2 0 0,2 0,4 0,6 0,8 1 1,2 1,4
Figure 5.15 – Evolution de l’amortissement et de la fréquence des 5 premiers modes propres en fonction
de l’amplitude des oscillations pour le tablier rectangulaire
147
5.3 Caractérisation du vent de la soufflerie
5.3.1 Dispositif expérimental
Les calculs dynamiques requièrent une connaissance fine de l’écoulement aéraulique au voisinage des
maquettes étudiées. Outre le profil de vitesse moyenne le long des tabliers, il nous faut mesurer les
intensités de turbulence, les échelles caractéristiques, et les coefficients de cohérence pour les
composantes horizontale et verticale du vent.
On caractérise donc le vent de la soufflerie atmosphérique à l’aide de deux sondes à fils chauds 2
composantes. Ces sondes sont fixées à 50 cm du sol sur des supports orientés dans l’axe de la veine sur
une ligne correspondant au bord d’attaque des maquettes aéroélastiques. Une sonde est placée au milieu
de la veine, en un point « 0 », l’autre est déplacée successivement en 12 points espacés de 20 cm,
numérotés de 1 à 6 et 1’à 6’ de part et d’autre de l’axe de la veine. Par la suite, on parlera respectivement
de points « à gauche » et « à droite » de l’axe central. Chacun des points correspond au milieu du bord
d’attaque d’un élément de tablier.
6 5 4 3 2 1 0 1’ 2’ 3’ 4’ 5’ 6’
2,4 m
4,0 m
Figure 5.16 – Vue en plan du montage fil chaud dans la veine atmosphérique
Un tube de Pitot, positionné légèrement en amont à 1,5 m du sol permet de mesurer la pression
dynamique de référence et le coefficient de recalage associé, utilisés au cours des manipulations
ultérieures.
148
Deux configurations sont étudiées :
 configuration turbulente : les rugosités placées en amont génèrent un taux de turbulence de l’ordre de
10 % longitudinalement et 6 % verticalement ; 4 vitesses moyennes de vent sont imposées : 6 m/s,
8 m/s, 10 m/s, 12 m/s ;
 configuration « laminaire » : la veine est vide et le taux de turbulence est inférieur à 2 % ; les vitesses
de vent générées sont de 6 m/s et 10 m/s.
Figures 5.17 et 5.18 - Supports et sondes à fils chauds dans la soufflerie en configuration turbulente
5.3.2 Résultats
a) Profils de vitesse moyenne et de turbulence
Les résultats en terme de gradients latéraux de vitesses moyennes et de taux de turbulence, dans les deux
configurations, sont présentés figures 5.19 à 5.24.
En configuration turbulente, les gradient de vitesse moyenne présentent des écarts maximaux de l’ordre
de ±7%. Les gradients de turbulence s’inscrivent eux dans des fourchettes de l’ordre de (8±3)% et
(6±1,5)%, respectivement pour les composantes u et w , i.e. ±37% et ±25% autour des valeurs moyennes.
L’importance des écarts nous amène à identifier cinq zones distinctes délimitées en rouge sur les
figures 5.19 à 5.24. L’allure particulière des courbes ne permet pas réellement de réduire ce nombre.
149
U(y)/U(0)
1,1 6m/s
8m/s
10m/s
1,0 12m/s
0,9 y(m)
-1,2 -1 -0,8 -0,6 -0,4 -0,2 0 0,2 0,4 0,6 0,8 1 1,2
Figure 5.19 – Gradient latéral de vitesse moyenne en vent turbulent
Iu(m/s)
12% 6m/s
10% 8m/s
8% 10m/s
6% 12m/s
4%
2%
0% y(m)
-1,2 -1 -0,8 -0,6 -0,4 -0,2 0 0,2 0,4 0,6 0,8 1 1,2
Figure 5.20 – Gradient latéral d’intensité de turbulence Iu en vent turbulent
Iw(m/s)
8% 6m/s
8m/s
6%
10m/s
4% 12m/s
2%
0% y(m)
-1,2 -1 -0,8 -0,6 -0,4 -0,2 0 0,2 0,4 0,6 0,8 1 1,2
Figure 5.21 – Gradient latéral d’intensité de turbulence Iw en vent turbulent
150
U(y)/y(0)
1,02
1,00 6m/s
10m/s
0,98
0,96 y(m)
-1,2 -1 -0,8 -0,6 -0,4 -0,2 0 0,2 0,4 0,6 0,8 1 1,2
Figure 5.22 – Gradient latéral de vitesse moyenne en vent laminaire
Iu(m/s)
2%
6m/s
10m/s
1%
0% y(m)
-1,2 -1 -0,8 -0,6 -0,4 -0,2 0 0,2 0,4 0,6 0,8 1 1,2
Figure 5.23 – Gradient latéral d’intensité de turbulence Iu en vent laminaire
Iw(m/s)
2%
6m/s
10m/s
1%
0% y(m)
-1,2 -1 -0,8 -0,6 -0,4 -0,2 0 0,2 0,4 0,6 0,8 1 1,2
Figure 5.24 – Gradient latéral d’intensité de turbulence Iw en vent laminaire
151
En configuration uniforme, les gradients de vitesse moyenne et d’intensité de turbulence sont beaucoup
plus faibles qu’en configuration turbulente. On peut donc raisonnablement admettre qu’une unique zone
suffit à représenter les caractéristiques statistiques du vent sur toute la largeur de la veine. On a toutefois
indiqué en rouge une possibilité de partition en cinq zones.
b) Echelles de turbulence et coefficients de cohérence
Les mesures effectuées permettent en outre d’estimer les fonctions densités spectrales et interspectrales
des composantes u et w pour les paires de points suivantes : (0,1), (0,2), (0,3), (0,4), (0,5), (0,6), (0,1’),
(0,2’), (0,3’), (0,4’), (0,5’), (0,6’). Il s’agit alors d’approcher ces fonctions par les modèles
analytiques inclus dans notre code de calcul : Karman et Davenport, respectivement pour les densités
spectrales et interspectrales (cf. chapitre 1). Ces fonctions sont paramétrées par un certain nombre de
variables qu’il faut donc déterminer : échelles de turbulences Lyu , Lyw pour les spectres et coefficients de
cohérence C uy , C wy pour les interspectres (ou, plus précisément, les fonctions de cohérence).
La détermination des échelles de turbulence ne pose pas de problème majeur. On s’arrange toutefois pour
que les courbes théoriques soient particulièrement bien ajustées pour la gamme de fréquences 2 - 15 Hz
contenant l’essentiel de l’énergie susceptible d’exciter nos maquettes.
En configuration turbulente, les valeurs de Lyu oscillent entre 0,4 et 0,6 m au point 0. On constate donc
que le vent de la soufflerie n’est pas, selon nos mesures, parfaitement stationnaire. On adoptera
néanmoins, sans trop de regrets, la valeur de 0,5 m, soit 500 m vraie grandeur. Ces valeurs diminuent
symétriquement et linéairement à mesure que l’on s’éloigne de l’axe central de la veine. On peut ainsi
estimer à 0,4 m les échelles de turbulence à 0,6 m de l’axe et 0,3 m à 1,2 m. On assigne donc les valeurs
0,5, 0,4, 0,3 m respectivement aux zones centrale, intermédiaires et extrêmes définies en 5.1.2a). Pour la
même configuration, les valeurs de Lyw varient de 0,075 m à 0,1 m. On ne constate pas de variations
significatives en fonction de la position par rapport à l’axe central.
En configuration uniforme, les valeurs de Lyu sont toujours proches de 0,1 m. On observe toutefois de
légères variations de chaque coté de l’axe : entre 0,1 m et 0,15 m à gauche et entre 0.075 m et 0.1 m à
droite. On attribue aux cinq zones, en partant de la gauche, les valeurs 0,125, 0,125, 0,1, 0,075, 0,1. Les
valeurs de Lyw oscillent entre 0,075 m et 0,1 m. En réalité, l’allure des signaux ne se prête pas assez à
l’ajustement de fonctions de Karman pour que se dégage une vraie tendance. On assigne donc la valeur
moyenne de 0,09 m à toutes les zones.
152
La détermination des coefficients de cohérence se révèle plus délicate. Autant les formulations de
Karman semblent particulièrement appropriées pour modéliser le spectre d’un vent en soufflerie, autant
les formulations de Davenport semblent inefficaces. La détermination du coefficient optimum pour un
écart y donné est déjà assez hasardeuse : on est en effet obligé de surestimer les degrés de corrélations
dans les basses fréquences et de les sous-estimer dans les hautes fréquences. Ces coefficients se révèlent
surtout totalement obsolètes lorsqu’on doit les appliquer à un écart y différent.
C’est la raison pour laquelle nous avons introduit un nouveau modèle, directement inspiré de celui de
Davenport, mais plus adapté à notre cas. Ce modèle n’a, en aucun cas, la prétention d’être universel mais
simplement de transférer efficacement des données expérimentales bien précises vers notre code de
calcul. Deux modifications ont été apportées à la fonction de cohérence de Davenport, dont on rappelle la
formulation appliquée à notre étude :
  n y C iy 
 i (y, n )  exp 
2
i=u, v (5.1)
 U 
 
La première modification consiste à « diminuer » l’influence du terme y en le remplaçant par y .
La seconde consiste à introduire deux fonctions « limites » de la forme (5.1), représentant, l’une les
basses fréquences, l’autre les hautes fréquences. La nouvelle fonction de cohérence résulte d’une
moyenne pondérée entre ces deux fonctions limites, dont le coefficient varie en fonction de la fréquence.
On introduit un terme supplémentaire nc, fréquence caractéristique pour laquelle les deux fonctions
limites ont le même poids. On écrit ainsi :
  n y C y   y C 2 i 
 exp  n
y
exp   1  0,5 n n c
n 1i
 i (y, n )  0,5
2 nc
 U   
  U  (5.2)
   
Remarquons que cette nouvelle fonction tend bien vers 1 quand ny tend vers 0 et vers 0 quand ny tend
vers l’infini. Le problème est que les termes situés dans les exponentielles ne sont pas adimensionnels.
Cela peut être artificiellement résolu en ajoutant au dénominateur un terme homogène à la racine carrée
d’une distance : dimension caractéristique de la soufflerie, échelle de turbulence… En fait, le sens
physique de ce terme n’a d’importance que dans l’optique d’une généralisation de notre formulation… ce
qui n’est pas le cas… On laisse donc l’expression (5.2) telle quelle.
En configuration turbulente, on estime les couples de coefficients à (5, 25) et (4, 20), respectivement
pour les composantes u et w. Notons que les interspectres liés à w sont plus faibles donc plus bruités et
plus difficiles à modéliser. En configuration laminaire, on trouve les couples de coefficients (6, 25) et
(7, 28) respectivement pour les composantes u et w.
153
c) Récapitulatif des résultats
En configuration turbulente :
y y y y
Zones U / U ref I u (%) I w (%) Lyu (m) Lyw (m) C1 u C2u C1 w C2w
1 1,04 6,5 5,5 0,3 0,09 5 25 4 20

2 0,97 10 6,5 0,4 0,09 5 25 4 20
3 1 8 5 0,5 0,09 5 25 4 20
4 1,02 6,5 4,5 0,4 0,09 5 25 4 20
5 0,97 8 5,5 0,3 0,09 5 25 4 20
Tableau 5.2 – Modélisation du vent en configuration turbulente
5 4 3 2 1
Figure 5.25 – Localisation des zones de vent en configuration turbulente
En configuration uniforme :
y y y y
Zones U / U ref I u (%) I w (%) Lyu (m) Lyw (m) C1 u C2u C1 w C2w
1 1 1,3 1,5 0,125 0,09 6 25 7 28

2 1,005 0,8 0,75 0,125 0,09 6 25 7 28
3 1 1,1 1,25 0,1 0,09 6 25 7 28
4 0,985 1,1 1,25 0,075 0,09 6 25 7 28
5 0,97 1,5 1,5 0,1 0,09 6 25 7 28
Tableau 5.3 – Modélisation du vent en configuration laminaire
5 4 3 2 1
Figure 5.26 – Localisation des zones de vent en configuration uniforme
154
0.2
Spectre signal 1 0.18 Interspectre

0.25 Spectre signal 2 Davenport C=15
0.16
Karman L=0.6m Davenport C=30
Karman L=0.4m 0.14 Dav. mod. C=(5,25)
0.2
S(m²/s²) 0.12
S(m²/s²)
0.15 0.1
0.08
0.1
0.06
0.04
0.05
0.02
0 0
2 5 10 20 50 100 2 5 10 20 50 100
f(Hz) f(Hz)
Figure 5.27 – Traitement de la composante u aux points 0 et 1 pour U=6m/s en vent turbulent
0.2

0.16
0.2
0.12
S(m²/s²)
0.15 S(m²/s²) 0.1
0.08
0.1
0.06
0.04
0.05
0.02
0 0
2 5 10 20 50 100 2 5 10 20 50 100
f(Hz) f(Hz)
0.1

0.08
0.2
0.06
S(m²/s²)
S(m²/s²)
0.15 0.05
0.04
0.1
0.03
0.02
0.05
0.01
0 0
2 5 10 20 50 100 2 5 10 20 50 100
f(Hz) f(Hz)
0.4 0.1

0.35
Spectre signal 2 Davenport C=15
0.08
0.3 Karman L=0.4m Davenport C=30
0.25
0.06
S(m²/s²)
S(m²/s²)
0.2 0.05
0.04
0.15
0.03
0.1
0.02
0.05
0.01
0 0
2 5 10 20 50 100 2 5 10 20 50 100
f(Hz) f(Hz)
155
5.4 Caractérisation aérodynamique des tabliers
Les valeurs des coefficients aérodynamiques stationnaires et instationnaires ont déjà été présentés au
chapitre 3. Il s’agit à présent de les modéliser en vue du calcul numérique.
Les coefficients aérodynamiques stationnaires sont linéarisés autour des incidences –5°, -4°, -3°, -2°, -1°,
0°, 1°, 2°, 3°, 4°, 5°, puis approchés par des polynômes du quatrième degré. Les résultats, obtenus en
vent uniforme et turbulent, sont présentés figures 5.31 et 5.32.
0.9
0.8
0.08 0.8 0.08
0.7
0.7
dCx/di
dCx/di
0.06 0.06
Cx
Cx
0.6 0.6
0.04 0.04
0.5
0.5
0.02 0.4 0.02
−5 0 5 −5 0 5 −5 0 5 −5 0 5
Incidence (deg.) Incidence (deg.) Incidence (deg.) Incidence (deg.)
6.5
0 7 0
−0.2 −0.2 6
6
dCy/di
dCy/di
Cy
Cy
−0.4 −0.4
5 5.5
−0.6 −0.6
−0.8 4 −0.8
5
−5 0 5 −5 0 5 −5 0 5 −5 0 5
0.1 1.5 0.1 1.5
0.05 0.05
dCm/di
dCm/di
1
Cm
Cm
1
0 0
0.5
−0.05 −0.05
0.5
−5 0 5 −5 0 5 −5 0 5 −5 0 5
Figure 5.31 – Linéarisation des coefficients aérodynamiques du tablier trapézoïdal

A gauche : vent uniforme ; à droite : vent turbulent
0.145 0.145
0.3
0.3
0.14 0.14 0.25
dCx/di
dCx/di
Cx
Cx
0.2
0.2
0.135 0.1 0.135
0.15
−5 0 5 −5 0 5 −5 0 5 −5 0 5
0.5 0.5
8 7
6
dCy/di
dCy/di
6
Cy
Cy
0 0
4 5
−0.5 −0.5 4
2
−5 0 5 −5 0 5 −5 0 5 −5 0 5
0.01 0.3 0.01

0.2
dCm/di
dCm/di
Cm
Cm
0 0.2 0
0.15
0.1
−0.01 −0.01
0.1
−5 0 5 −5 0 5 −5 0 5 −5 0 5
Figure 5.32 – Linéarisation des coefficients aérodynamiques du tablier rectangulaire

A gauche : vent uniforme ; à droite : vent turbulent
156
Les coefficients instationnaires de Scanlan sont eux approchés par les fonctions indicielles présentées
chapitre 2. Ces fonctions sont ajustées autour des données issues des essais en oscillations libres (figures
5.33 et 5.34).
0 −1
−2
−0.1
−3
H1*
−0.2
A2*
−4
−0.3 −5
−6
−0.4
−7
−0.5 −8
0 1 2 3 4 5 6 7 8 9 0 2 4 6 8 10
Vitesse réduite Vitesse réduite
3
−1
2.5
2 −2
H4*
A3*
1.5 −3
1
−4
0.5
0 −5
0 1 2 3 4 5 6 7 8 9 0 2 4 6 8 10
Figure 5.33 – Fonctions indicielles caractérisant le tablier trapézoïdal
0.1 0
0
−5
−0.1
H1*
A2*
−0.2 −10
−0.3
−15
0 1 2 3 4 5 0 1 2 3 4 5 6 7 8 9 10
2
−2
1.5
−4
H4*
A3*
1
−6
0.5
−8
0
0 1 2 3 4 5 0 1 2 3 4 5 6 7 8 9 10
Figure 5.34 – Fonctions indicielles caractérisant le tablier rectangulaire
157
5.5 Mesure des déplacements
Les déplacements sont mesurés en quatre points : deux à mi-travée et deux à quart-de-travée (figures 5.35
et 5.36). Les écarts types des déplacements modaux sont obtenus par combinaison linéaire des quatre
voies, analyse spectrale et ajustement de spectres d’oscillateurs linéaires à 1 d.d.l. par les moindres
carrés. On en déduit également fréquences et amortissements modaux.
Figures 5.35 et 5.36 – Capteurs de déplacements et étiquettes de visée
Torsion 1 Flexion 1
0.1 0.018
Points expérimentaux Points expérimentaux
0.09 5.07Hz, 9.53e−001mm, 4.7% 0.016 7.48Hz, 3.54e−001mm, 2.4%
0.08 8.95Hz, 4.31e−001mm, 2.1%
0.014
0.07
0.012
S(mm²/s²)
S(mm²/s²)
0.06
0.01
0.05
0.008
0.04
0.006
0.03
0.02 0.004
0.01 0.002
0 0
2 5 10 20 50 2 5 10 20 50
f(Hz) f(Hz)
−3
Torsion 2 x 10 Flexion 2
6
0.045
Points expérimentaux Points expérimentaux
0.04 5.08Hz, 6.39e−001mm, 4.7% 5 7.45Hz, 2.01e−001mm, 2.5%
0.035 8.97Hz, 3.32e−001mm, 2.3%

11.39Hz, 3.07e−001mm, 1.9%
4
0.03
S(mm²/s²)
S(mm²/s²)
0.025
3
0.02
0.015 2
0.01
1
0.005
0 0
2 5 10 20 50 2 5 10 20 50
f(Hz) f(Hz)
Figure 5.37 – Analyse spectrale des déplacements modaux : tablier trapézoïdal en vent turbulent à 10m/s
158
La figure 5.38 montre le développement d’une instabilité en torsion pure du tablier trapézoïdal en vent
uniforme à une vitesse de 20 m/s. Cela correspond, d’après les échelles de similitude (1/100 en distance
et 25/1 en fréquence), à une vitesse vraie grandeur de l’ordre de 80 m/s.
Figure 5.38 – Instabilité en torsion pure du tablier de Millau à 20m/s (24 images/seconde)
159
5.6 Résultats
5.6.1 Considérations sur les incertitudes de calcul
Avant de porter un jugement sur la comparaison entre résultats numériques et expérimentaux, il convient
d’avoir une idée des incertitudes entachant les uns et les autres. Les premiers sont liées aux erreurs
d’estimation des différents paramètres utilisés pour le calcul. Les seconds concernent simplement la
précision de la mesure des déplacements par les capteurs lasers.
En se référant à la formulation (4.67), on peut proposer une estimation de l’incertitude relative

concernant l’estimation numérique de l’amplitude des déplacements . On écrit ainsi :
 q 1 S F  1  m g 3 
    (5.3)
q 2 S F  2  mg 2 
l’incertitude relative de la densité spectrales des buffeting forces pouvant être décomposée comme suit :
S F  C ad S u U
2  2 (5.4)
S F  C ad Su U
où Cad représente « l’ensemble » des coefficients aérodynamiques.
Le tableau 5.4 rassemble les incertitudes relatives estimées « à la louche » des différents paramètres.
L’incertitude concernant les coefficients aérodynamiques est estimée essentiellement à partir des écarts
obtenus entre mesures par balances et mesures par prises de pression (figures 3.5-3.6). Elle ne prend en
compte ni l’incertitude liée à la détermination de l’admittance aérodynamique, ni les problèmes
spécifiques liés à la détermination des coefficients de moment du tablier rectangulaire, ni les effets de
Reynolds.
L’incertitude dans la détermination du spectre du vent est estimée d’après les figures 5.27-30.
L’incertitude liée à l’amortissement ne prend en compte que l’amortissement aérodynamique qui est
prépondérant sur l’essentiel de la gamme de vitesses ; elle est estimée d’une part d’après les
« impressions » ressenties pendant les mesures et traitements des coefficients aéroélastiques, d’autre part
d’après l’évaluation de la pertinence des fonctions indicielles (figures 5.33 et 5.34).
160
Les incertitudes concernant masses généralisées et pulsations proviennent de l’analyse des nuages de
points expérimentaux (figures 5.14 et 5.15).
Paramètre Cad Su U  mg 
P/P 5% 5% 1% 20% 5% 0,5%
Tableau 5.4 – Incertitudes relatives des paramètres utilisés par les calculs spectraux et temporels
Des expressions (5.3) et (5.4), on estime l’incertitude relative dans l’estimation des déplacements, liée
aux erreurs expérimentales, à 25% environ.
Si l’on estime l’incertitude entachant la mesure des déplacements (détermination des gains, position du
faisceau laser) à 5%, on arrive à une incertitude relative totale de 30%, uniquement liée à l’ensemble des
mesures effectuées pendant l’ensemble de la procédure de validation.
5.6.2 Calcul spectral
Les calculs par la méthode spectrale ont été effectués pour les deux tabliers et les deux types de vent, soit
quatre configurations. Pour chacune d'elles, les réponses des quatre premiers modes de vibration sont
présentées figures 5.39 à 5.42. On y trouve :
 les résultats issus des essais expérimentaux;

 les résultats obtenus par le calcul spectral "complet", c'est à dire avec prise en compte des forces
aéroélastiques et de l'admittance aérodynamique;
 les résultats obtenus sans prise en compte de l'admittance aérodynamique;
 les résultats obtenus par le calcul spectral "quasi statique", c'est à dire sans tenir compte ni de
l'admittance aérodynamique, ni des forces aéroélastiques.
Les résultats expérimentaux sont représentés sous forme de cercles dont le rayon, dans la représentation
en échelle logarithmique, correspond à une incertitude de 20%.
161
a) Cas du tablier trapézoïdal
Les réponses des maquettes sont globalement sous-estimées pour les faibles vitesses de vent et
surestimées pour les hautes vitesses. Cela est partiellement imputable à l'évolution de l'amortissement
structurel en fonction de l'amplitude des déplacements: l'amortissement structurel retenu pour les calcul
est surévalué pour les petites amplitudes et sous-évalué pour les grandes amplitudes. Cela est surtout vrai
à faibles vitesse et amplitude lorsque l'amortissement aérodynamique est faible par rapport à
l'amortissement structurel. Ajoutons qu'à faible vitesse les mesures peuvent être surestimées du fait de
l'importance relative croissante des oscillations générées par les vibrations du sol de la soufflerie.
Les résultats obtenus en vent uniforme ne sont pas plus mauvais que ceux obtenus en vent turbulent. On
aurait pu s'attendre à une sous-estimation des réponses en vent laminaire du fait de la non prise en
compte des effets de sillage. On peut en conclure que ces effets ne sont pas trop importants.
La vitesse critique de flottement est significativement sous-estimée en vent uniforme et très largement
sous estimée en vent turbulent. La première explication réside dans la mauvaise prédiction des
phénomènes aéroélastiques. Rappelons en effet que les coefficients de Scanlan ont été déterminés dans
des conditions où le taux de turbulence était de l'ordre de 0,5%. Ici, les taux de turbulence pour la
composante longitudinale du vent sont de 1,5% et 8% respectivement pour les configurations uniformes
et turbulentes. Or, les phénomènes d'instabilité aérodynamique ont d'autant plus de mal à se développer
que le taux de turbulence est fort, même s’ils sont également plus facilement initiés. Cela n'est
absolument pas pris en compte dans nos calculs, ce qui explique la sous-estimation des vitesses critiques
de flottement. Ajoutons que la superposition de différents modes de vibration ne peut que gêner le
développement de phénomènes aérodynamiques organisés. Cela non plus n'est pas pris en compte dans
nos calculs. La seconde explication réside dans la sous-estimation de l'amortissement structurel du
premier mode pour les grandes amplitudes. La vitesse permettant d'annuler l'amortissement total est donc
sous-estimée.
Les résultats obtenus par une approximation "quasi statique" montrent bien l'importance des forces
aéroélastiques et admittance aérodynamique. En particulier, l'amortissement aérodynamique joue ici un
grand rôle du fait de la faiblesse relative des amortissements structurels.
b) Cas du tablier rectangulaire
Le comportement des modes de torsion est globalement mal simulé par le calcul : les déplacements sont
généralement sous-estimés. L'effet du sillage est sans doute ici plus important que dans le cas du tablier
trapézoïdal ; on constate notamment le rôle essentiel joué par le détachement tourbillonnaires. Rappelons
162
surtout les difficultés rencontrées lors de la détermination expérimentale des coefficients de moment.
L'introduction d'une admittance aérodynamique "fourre-tout" permet certes de corriger partiellement la
mauvaise connaissance des moments s'exerçant tout au long du tablier. Néanmoins, ce terme correctif
n'est strictement valable que dans le cas d'un tablier d'un mètre de long. On a en effet pu voir l'importance
des effets tridimensionnels liés à ce type de profil.
La réponse du mode de flexion est, par contre, bien prédite par le calcul. Cela peut être mis au crédit de
la plus forte prédominance des forces aéroélastiques par rapport au tablier trapézoïdal (fig. 3.19 et 3.20)
qui réduit le problème lié à la mauvaise estimation de l'amortissement structurel.
La vitesse critique de flottement en vent uniforme est également correctement prédite. Une explication,
sans doute insuffisante, est que les forces de traînée, jouent ici un rôle majeur. Le premier mode, dont la
déformée horizontale est très importante, est donc favorisé et moins perturbé par les autres modes. En
vent turbulent, la vitesse critique de flottement est par contre logiquement sous-estimée.
c) Conclusion
Compte tenu des considérations énoncées en 5.5.1 et des difficultés rencontrées autour de la
détermination du moment du tablier rectangulaire, on peut affirmer que les résultats obtenus par la
méthode spectrale sont tout à fait satisfaisant dans une gamme de vitesse « intermédiaire ». A basse
vitesse, les calculs sont globalement nettement surestimés, spécialement en vent uniforme. On peut
avancer deux hypothèses :
 Les effets de Reynolds deviennent importants : on a déjà observé des effets non négligeables lors de
l’étude paramétrique des coefficients de moments (figure 3.8) ;
 La turbulence du vent est mal modélisée à faible vitesse : les mesures de turbulence, faites pour des
vitesses comprises entre 6 m/s et 12 m/s, ont été appliquées brutalement à l’ensemble de la gamme
de vitesses .
Enfin, les instabilités sont assez mal prédites. Les raisons principales sont :
 la non prise en compte de l’impact de la turbulence sur les coefficients aéroélastiques ;

 l’absence de modélisation de la cohérence latérale des forces aéroélastiques ;
 les limites du principe de superposition des modes.
163
Mode 1
E.t. déplacements (mm)

0
10
−1
10
Points expérimentaux
−2
10 Calcul spectral
C.S. sans admittance
C.S. quasi−statique
5 10 15 20 25
U(m/s)
Mode 2
0
10
−1
10
−2
10
Calcul spectral
5 10 15 20 25
U(m/s)
Mode 3
1
10
0
10
−1
10
−2
10 Points expérimentaux
Calcul spectral
5 10 15 20 25
U(m/s)
Mode 4
0
10
−1
10
−2
10
Calcul spectral
5 10 15 20 25
U(m/s)
Figure 5.39 – Résultats calculs spectraux : tablier trapézoïdal en vent uniforme
164
Mode 1

1
10
0
10
−1 Points expérimentaux
10 Calcul spectral
−2
10 C.S. quasi−statique
5 10 15 20 25
U(m/s)
Mode 2
1
10
0
10
−1
10
−2 Calcul spectral
10 C.S. sans admittance
5 10 15 20 25
U(m/s)
Mode 3
1
10
0
10
−1
10
−2
10 Calcul spectral
5 10 15 20 25
U(m/s)
Mode 4
1
10
0
10
−1
10
−2
10 Calcul spectral
5 10 15 20 25
U(m/s)
Figure 5.40 – Résultats calculs spectraux : tablier trapézoïdal en vent turbulent
165
Mode 1
1
10
0
10
−1
10
−2 Calcul spectral
10
2 4 6 8 10 12 14
U(m/s)
Mode 2
0
10
−1
10
−2
Calcul spectral
2 4 6 8 10 12 14
U(m/s)
Mode 3
−1
10
−2
10
Calcul spectral
2 4 6 8 10 12 14
U(m/s)
Mode 4
−1
10
−2
10
Calcul spectral
2 4 6 8 10 12 14
U(m/s)
Figure 5.41 – Résultats calculs spectraux : tablier rectangulaire en vent uniforme
166
Mode 1
1
10
0
10
−1
10
Calcul spectral
−2 C.S. sans admittance
10
2 4 6 8 10 12 14 16 18 20
U(m/s)
Mode 2
0
10
−1
10
Calcul spectral
−2
10 C.S. sans admittance
2 4 6 8 10 12 14 16 18 20
U(m/s)
Mode 3
0
10
−1
10
−2 Points expérimentaux
10 Calcul spectral
2 4 6 8 10 12 14 16 18 20
U(m/s)
Mode 4
0
10
−1
10
−2
Calcul spectral
2 4 6 8 10 12 14 16 18 20
U(m/s)
Figure 5.42 – Résultats calculs spectraux : tablier rectangulaire en vent turbulent
167
5.6.3 Calcul temporel
Pour le calcul temporel, on a généré une série de signaux de vent, en 24 points de l’espace, pour
différentes valeurs de vitesses moyennes de référence comprises entre 1 et 25 m/s. Chaque signal est
composé de 8192 points espacés de 1/100e de seconde. Chaque période du quatrième mode est ainsi
représentée par environ 8 points. Deux signaux de la composante longitudinale d’un vent turbulent sont
représenté figure 5.43. On a également représenté leurs spectres et interspectres, réels et imposés en
condition d’entrée.
Les résultats issus du calcul temporel sont présentés figures 5.47 à 5.50. On y a représenté :
 les résultats issus des essais expérimentaux;

 les résultats obtenus par le calcul temporel "complet", c'est à dire avec prise en compte des forces
aéroélastiques et de l'admittance aérodynamique;
 les résultats obtenus sans prise en compte des forces aéroélastiques;
 un rappel des résultats obtenus auparavant par le calcul spectral.
Signal 1: point 10 , composante 1

12
U(m/s)
10
8
0 10 20 30 40 50 60 70 80
t(sec.)
Signal 2: point 11 , composante 1

12
U(m/s)
10
8
0 10 20 30 40 50 60 70 80
t(sec.)
0.4
Spectre 1
0.35
Spectre 2
Interspectre
0.3
Sp. Karman 1
0.25 Sp. Karman 2
Intersp. Dav.
S(m²/s²)
0.2
0.15
0.1
0.05
0
1 2 5 10 20 50
f(Hz)
Figure 5.43 – Signaux de vent turbulent générés en deux points distants de 20 cm
168
Une bonne partie des commentaires formulées à propos du calcul spectral sont applicables au calcul
temporel. On peut cependant apporter deux nuances.
Tout d’abord les déplacements à faible vitesse sont moins sous-estimés. Cela relativise (sans toutefois les
contredire) les analyses proposées à ce sujet à propos du calcul spectral puisque les arguments invoqués
alors sont toujours valables ici. A vrai dire, ne pouvant expliquer ce fait, il nous est difficile de prétendre
que le calcul temporel est plus précis à faible vitesse que le calcul spectral. Les raisons, inconnues, sont
peut-être de mauvaises…
Par contre la meilleure prédiction du comportement instable des tabliers peut en partie être mise au crédit
d’une meilleure prise en compte de l’impact de la turbulence sur les phénomènes aéroélastiques. On voit
notamment sur la figure 5.44, représentant les signaux obtenus avec et sans prise en compte des effets
aéroélastiques que l’apport de ces derniers ne se résume pas à une simple changement d’échelle. On
pourrait presque parler de « mélange intime » entre les buffeting forces et les forces aéroélastiques,
représentatifs d’une « partie de la réalité ». Si l’on s’intéresse maintenant aux figures 5.45 et 5.46,
représentant les oscillations des deux premier modes à 2, 10 et 19 m/s en configurations laminaire et
turbulente, on s’aperçoit que dans le premier cas l’instabilité à 19 m/s est immédiate, alors que dans le
second cas les buffeting forces semblent retarder l’apparition de cette même instabilité. Le fait que les
vitesses critiques obtenues en vent laminaire et en vent turbulent soient similaires provient uniquement
du fait que les temps de simulations sont suffisamment longs pour que, quelque soit le degré de
turbulence, les instabilités potentielles se muent en instabilités réelles. Cependant, on peut très bien
imaginer qu’avec un taux de turbulence supérieur et une durée de simulation égale, la vitesse de 19 m/s
soit déclarée inférieure à la vitesse critique de flottement.
4
Déplacements (mm)
−2
10 20 30 40 50 60 70 80
Temps (sec.)
Déplacements (mm)
−2
−4
10 20 30 40 50 60 70 80
Temps (sec.)
1
Figure 5.44 – Signaux de déplacement des deux premiers mode du tablier trapézoïdal (mode 1 en bleu,
mode 2 en vert) en vent turbulent pour une vitesse de 10 m/s avec et sans effets aéroélastiques
169
0.06
Déplacements (mm)
0.04
0.02
0
−0.02
−0.04
2 4 6 8 10 12 14 16 18 20
Temps (sec.)
4
Déplacements (mm)
−2
2 4 6 8 10 12 14 16 18 20
Temps (sec.)
100
Déplacements (mm)
50
−50
−100
−150
2 4 6 8 10 12 14 16 18 20
Temps (sec.)
Figure 5.45 – Signaux de déplacement des deux premiers mode du tablier trapézoïdal
en vent turbulent pour 2 m/s, 10 m/s, 19 m/s
Déplacements (mm)
0.02
0.01
1 2 3 4 5 6 7
Temps (sec.)
1.5
Déplacements (mm)
0.5
1 2 3 4 5 6 7
Temps (sec.)
100
Déplacements (mm)
50
−50
−100
1 2 3 4 5 6 7
Temps (sec.)
Figure 5.46 – Signaux de déplacement des deux premiers mode du tablier trapézoïdal
en vent uniforme pour 2 m/s, 10 m/s, 19 m/s
170
Mode 1
0
10
−1
10
−2
10 Calcul temporel
C.T. sans forces aéroélastiques
Calcul spectral
5 10 15 20 25
U(m/s)
Mode 2
−1
10
−2
10
Calcul temporel
Calcul spectral
5 10 15 20 25
U(m/s)
Mode 3
1
10
0
10
−1
10
−2
Calcul temporel
Calcul spectral
5 10 15 20 25
U(m/s)
Mode 4
−1
10
−2
10
Calcul temporel
Calcul spectral
5 10 15 20 25
U(m/s)
Figure 5. 47 – Résultats calculs temporels : tablier trapézoïdal en vent uniforme
171
Mode 1
1
10
0
10
−1
10
Calcul temporel
−2 C.T. sans forces aéroélastiques
10 Calcul spectral
5 10 15 20 25
U(m/s)
Mode 2
0
10
−1
10
−2
10 Calcul temporel
Calcul spectral
5 10 15 20 25
U(m/s)
Mode 3
0
10
−1
10
−2
10 Calcul temporel
Calcul spectral
5 10 15 20 25
U(m/s)
Mode 4
0
10
−1
10
−2
Calcul temporel
Calcul spectral
5 10 15 20 25
U(m/s)
Figure 5.48 – Résultats calculs temporels : tablier trapézoïdal en vent turbulent
172
Mode 1
1
E.t. déplacements (mm) 10
0
10
−1
10
−2
10 Calcul temporel
Calcul spectral
2 4 6 8 10 12 14 16
U(m/s)
Mode 2
−1
10
−2
10
Calcul temporel
Calcul spectral
2 4 6 8 10 12 14 16
U(m/s)
Mode 3
0
10
−1
10
−2
10
Calcul temporel
Calcul spectral
2 4 6 8 10 12 14 16
U(m/s)
Mode 4
−1
10
−2
10
Calcul temporel
Calcul spectral
2 4 6 8 10 12 14 16
U(m/s)
Figure 5.49 – Résultats calculs temporels : tablier rectangulaire en vent uniforme
173
Mode 1
1
E.t. déplacements (mm) 10
0
10
−1
10
Calcul temporel
−2
10 C.T. sans forces aéroélastiques
Calcul spectral
2 4 6 8 10 12 14 16
U(m/s)
Mode 2
0
10
−1
10
−2 Calcul temporel
10 C.T. sans forces aéroélastiques
Calcul spectral
2 4 6 8 10 12 14 16 18 20
U(m/s)
Mode 3
0
10
−1
10
−2
10 Calcul temporel
Calcul spectral
2 4 6 8 10 12 14 16 18 20
U(m/s)
Mode 4
0
10
−1
10
−2
Calcul temporel
Calcul spectral
2 4 6 8 10 12 14 16 18 20
U(m/s)
Figure 5.50 – Résultats calculs temporels : tablier rectangulaire en vent turbulent
174
5.7 Conclusion
Ce chapitre, consacré à la procédure de validation expérimentale des codes développés, était basé sur le
passage en soufflerie atmosphérique de deux maquettes aéroélastiques représentant des travées centrales
de tabliers de ponts. Les déplacements expérimentaux ont été comparés aux déplacements calculés par
les approches spectrales et temporelles. Compte tenu des incertitudes liées à la détermination des
différents paramètres d’entrée, les résultats se sont avérés tout à fait satisfaisants.
Des études complémentaires permettraient cependant de mieux définir les limites des approches
spectrales et temporelles. Des maquettes de type taut-tube pourraient notamment être employées pour
étudier les phénomènes de couplage flexion-torsion. On pourrait ensuite envisager l’étude de structures
plus complexes comme un pont en phase de construction (figure 4.1).
175
Conclusion
L’objectif de ce travail était de développer deux outils de calcul permettant d’estimer la réponse des
structures élancées aux sollicitations turbulentes du vent : le premier basé sur une approche spectrale, le
second basé sur une approche temporelle.
Pour cela, nous avons d’abord réalisé une synthèse bibliographique des différents modèles disponibles
aujourd’hui pour représenter, d’une part le vent naturel dans la couche limite atmosphérique, d’autre part
les forces aérodynamiques s’exerçant sur les ouvrages d’art. C’était l’objet des deux premiers chapitres.
Puis, nous avons présenté différentes méthodes expérimentales permettant d’estimer les efforts
aérodynamiques s’exerçant sur des maquettes sectionnelles. Nous nous sommes particulièrement
intéressé aux deux principales techniques permettant de mesurer les forces aéroélastiques sur les tabliers
de pont : la technique des oscillations libres et la technique des oscillations forcées. L’ensemble des
méthodes expérimentales présentées ont été appliquées à deux types de sections utilisées ensuite lors de
la procédure de validation du code.
Nous avons ensuite exposé les bases théoriques des approches spectrale et temporelle. Les deux sont
basées sur une discrétisation de la structure en éléments de type poutre et l’utilisation d’outils issus de la
mécanique des vibrations.
Dans le cas de l’approche spectrale, le vent et les forces aérodynamiques sont utilisées sous forme de
densités spectrales de puissance. Les forces aéroélastiques, elles, interviennent sous forme de matrices
raideur et amortissement aérodynamiques. On aboutit à un système matriciel exprimé dans la base
modale, dans lequel des termes de couplage aérodynamique peuvent intervenir. Ajoutons que la
déformée statique de la structure ainsi que l’incidence moyenne du vent sont désormais parfaitement pris
en compte.
L’approche temporelle repose elle sur la génération d’un signal de vent aléatoire tridimensionnel
respectant certaines conditions de corrélations spatio-temporelles. On a pour cela développé un nouveau
module de calcul basé sur l’algorithme de Jin et al., remplaçant l’ancien module basé sur l’algorithme de
Shinozuka.
176
Conclusion
Les forces aéroélastiques interviennent sous forme de fonctions indicielles. Celles-ci sont à l’origine d’un
système de deux équations matricielles couplées : l’une étant l’équation classique de la mécanique,
l’autre reliant une partie des forces aéroélastiques aux déplacements de l’ouvrage. Ce système est résolu
pas de temps par pas de temps à l’aide d’un schéma implicite de Newmark pour la première équation,
d’un schéma explicite pour la seconde.
Une méthode d’optimisation de la linéarisation des coefficients aérodynamiques, à destination des

calculs spectraux et basée sur des simulations temporelles, a été introduite. On a également proposé une
méthode de combinaison des charges modales extrêmes et conversion en charges statiques équivalentes,
à destination des bureaux d’études chargés du dimensionnement des ouvrages. Enfin les caractéristiques
statistiques du vent peuvent désormais être définies de façon zonale.
Nous avons mis en place une procédure de validation des deux approches, basée sur le passage en
soufflerie de deux maquettes aéroélastiques représentant des travées centrales de tabliers de ponts. Les
résultats numériques se révèlent globalement en bon accord avec les données expérimentales, compte
tenu des incertitudes liées à la détermination expérimentales des caractéristiques mécaniques et
aérodynamiques de nos maquettes. On a toutefois noté une sous-estimation des déplacements pour les
faibles vitesses de vent, vraisemblablement due à des effets de Reynolds et à une mauvaise modélisation
de la turbulence. Les vitesses critiques de flottement s’avèrent également sous-évaluées : sont en cause le
principe de superposition des phénomènes aéroélastiques modaux, l’évolution de l’amortissement
structurel en fonction de l’amplitude, la mauvaise prise en compte de l’impact de la turbulence sur le
développement de phénomènes aéroélastiques et, enfin, l’absence de modélisation de la cohérence
latérale des forces aéroélastiques. On a cependant pu noter que l’impact de la turbulence sur le
développement des instabilités aérodynamiques était mieux appréhendé par la méthode temporelle.
177
Perspectives
L’objectif de ce travail était non seulement de développer un outil de calcul, mais également de situer ce
dernier dans son contexte expérimental. Ainsi, on a vu qu’une bonne partie des écarts entre l’estimation
numérique des déplacements et leur détermination expérimentale pouvaient être attribuée aux erreurs de
mesure des différents paramètres d’entrée. L’amélioration de l’outil numérique passe d’abord par
l’amélioration de la détermination de ces paramètres d’entrée.
En particulier, il serait souhaitable d’améliorer la mesure des amortissements structurels et

aérodynamiques, de bien caractériser leur évolution en fonction de l’amplitude des oscillations, et, pour
le second, en fonction du nombre de Reynolds et du taux de turbulence. On améliorerait ainsi d’une part
la prévision de la réponse à la turbulence, d’autre part la prévision des instabilités aérodynamiques.
Il serait également utile de développer la mesure des coefficients aéroélastiques en oscillations couplées :
soit sur des maquettes sectionnelles (un mode de flexion et un mode de torsion) soit sur des maquettes
souples (plusieurs modes de flexion et de torsion). Cela permettrait de réduire les erreurs liées au
principe de superposition des effets aéroélastiques modaux.
L’influence du Reynolds et de la turbulence sur les coefficients aérodynamiques stationnaires mériterait

aussi aujourd’hui d’être mieux caractérisée. La fonction d’admittance aérodynamique permet certes de
prendre en compte l’impact relatif des différentes échelles de la turbulence mais rien ne dit que celle-ci
n’évolue pas avec l’allure du spectre de la turbulence (i.e. de la vitesse moyenne et des intensités et
échelles de turbulence).
Un point important, qui n’a pas été abordé dans ce mémoire, concerne la caractérisation de la cohérence
latérale des forces aérodynamiques. L’expérience montre en effet que celle-ci est souvent plus importante
que celle du vent. La non prise en compte de l’accroissement de la cohérence latérale conduit à une sous-
estimation de la réponse des modes primaires.
L’ensemble des améliorations effectuées sur le plan expérimental doivent naturellement être
accompagnées de modèles mathématiques permettant leur prise en compte dans le cadre d’un calcul
numérique. Ceux-ci sont a priori plus faciles à définir dans l’optique de calculs spectraux,
178
l’expérimentation se ramenant souvent à des concepts statistiques. Il nous semble cependant primordial
de les appliquer également au calcul temporel, seul à pouvoir appréhender réellement le développement
des instabilités aérodynamiques.
En ce qui concerne l’évaluation des coefficients aérodynamiques stationnaires et instationnaires, de

nombreux outils numériques existent pouvant potentiellement remplacer tout ou partie des essais en
soufflerie. Les plus courants sont basés sur la résolution des équations de Navier-Stokes.
Les codes de calculs de type k- sont conçus pour appréhender les phénomènes aérodynamiques
moyens ; les fluctuations liées à la turbulence ne sont pas calculées mais modélisées par un terme unique.
Ces outils présentent l’intérêt d’être économiques en temps de calcul ; cependant, ils manquent
cruellement de précision pour prédire les décollements et recollement des écoulements, du fait d’une
surproduction d’énergie cinétique aux points d’impacts. Or la connaissance de ces phénomènes est
primordiale pour l’estimation des coefficients aérodynamiques stationnaires.
Les codes de type L.E.S.(Large Eddy Simulation), calculant l’histoire des structures turbulentes jusqu’à
de très petites échelles, permettent d’appréhender les phénomènes aérodynamiques instationnaires.
Cependant, on est souvent obligé de se ramener à des calculs bidimensionnels pour limiter les temps de
calcul. La turbulence, de nature tridimensionnelle, est alors forcément moins bien appréhendée. En outre
l’analyse des couplages fluide-structure nécessitent l’emploi de maillages mobiles, sources de difficultés
supplémentaires.
De nouvelles techniques de calcul des écoulements sont apparues il y a quelques années pour répondre
aux besoins spécifiques du Wind Engineering, regroupées sous l’appellation de Discret Vortex Method.
Celles-ci sont basées sur un certain nombre de simplifications (écoulement bidimensionnel, non
visqueux) permettant notamment de réduire le nombre de variables à un seul terme: la vorticité. En outre,
elles ne nécessitent aucun maillage, ce qui facilite l’étude des couplages fluide-structure. Ces méthodes
simplifiées peuvent constituer un compromis intéressant entre précision et rapidité des calculs.
L’utilisation de la C.F.D. dans des calculs complets de couplage fluide-structure semble, à terme,
inéluctable. La L.E.S. constituera probablement, un jour, la meilleure solution. La Discret Vortex
Method semble cependant être, à moyen terme, l’approche la plus adaptée. On peut en effet d’ores et déjà
envisager d’appliquer à chacun des éléments de la structure un calcul DVM bidimensionnel. Il s’agirait
alors d’introduire des « degrés de couplage » appropriés entre ces calculs parallèles.
179
Bibliographie
Attou M. (1994) – Réponse dynamique des structures élancées à la turbulence du vent. Thèse de
doctorat, Ecole Centrale de Nantes.
Attou M., Grillaud G. (1993) – Dynamic response of structures to wind turbulence. Proceedings of the
2nd European Conference on structural Dynamics, Trondheim, Norway, pp. 1035-1042.
Barré C. et Brebion H. (1990) – Coefficients aérodynamiques du tablier du pont de Saint-Nazaire.

Rapport CSTB, Réf. EN-AS90.6C.
Barré C. (1996) – Coefficients aérodynamiques des versions béton et mixte du tablier et des piles du
viaduc de Verrières. Rapport CSTB, Réf. EN-ECA96. 7C.
Barré C. (1998) – A75 – Grand viaduc de Millau – Coefficients aérodynamiques du tablier et des piles.
Rapport CSTB, Réf. EN-AEC98.37C.
Biétry J., Sacré C., Simiu E. (1978) – Mean wind profiles and change of terrain roughness. ASCE, Vol.
104, n° ST10
Biétry J. et Grillaud G. (1984) – Wind studies for the Normandy bridge. Actes de la conférence de
l'AFPC sur les ponts suspendus et à haubans, Deauville, France, Vol.2, pp. 29-44.
Bisplinghoff R.L. and Ashley H. (1962) - Principles of Aeroelasticity. John Wilwy and sons, Inc., New
York, USA.
Bleich F. (1948) – Dynamic instability of truss-stiffned suspension bridges under wind action. Proc.,
ASCE, Vol. 74, No.8, pp. 1269-1314.
180
Bibliographie
Bucher C.G., Lin Y.K. (1988) – Stochastic stability of bridges considering coupled modes. Journal of
Engineering Mechanics, ASCE, Vol. 114, n° 12, pp. 2055-2071.
Davenport A. G. (1961a) – The application of statistical concepts to wind loading of structures.

Proceedings of Institution of Civil Engineers, 19, pp. 449-472.
Davenport A. G. (1961b) – The spectrum of horizontal gustiness near the ground in high winds. Journal
of the Royal Meteorology Society, 87, pp. 194-211.
Davenport A. G. (1962) – The buffeting response of a suspension bridge by storm winds. Journal of the
Structural Division, ASCE, 88, pp. 233-268.
De La Foye A., Flamand O. and Grillaud G. (2000) – Measurement of aerodynamic derivatives with the
forced oscillation technique. Abstract for BBAA IV congress, Bochum.
De La Foye A. and Flamand O. (2001) – Time domain simulation of vibrating bridge deck under wind
action. Abstract for EAWE congress, Eindhoven.
Delaunay D. (1994) – Les caractéristiques du vent sur le site du grand viaduc de Millau. Rapport CSTB,
Réf. EN-CLI 94.16 C
Dubigeon S. (1993) – Mécanique des vibrations, systèmes à plusieurs degrés de liberté. Cours de l’Ecole
Centrale de Nantes.
Flamand O. et Grillaud G. (1998) – Viaduc de Millau : comportement au vent de la pili P2 en phase de

construction. Rapport CSTB, Réf. EN-AEC98.89C.
Flamand O. and De La Foye A., (2001) – The forced oscillation technique and L.ES. simulation. Abstract
for EAWE congress, Eindhoven.
Gandemer J. (1981) – Aérodynamique. Sciences du bâtiment, Vol. 2. Centre Scientifique et Technique du

Bâtiment, France.
181
Bibliographie
Gibert R.J. (1988) – Vibrations des structures. Interactions avec les fluides, sources d’excitation
aléatoires. Editions Eyrolles.
Grillaud G., Bourcier P., Barré C., Flamand O. (1997) – Wind action on the Vasco da Gama cable stayed
bridge. 2nd European and African Conference on Wind Engineering, Genova, Italy, Vol. 2, pp. 1449-
1456.
Grillaud G., Flamand O., Barré C. (1991) – Comportement au vent du pont de Normandie :étude en
soufflerie sur maquette aéroélastique à échelle 1/200e . Rapport CSTB, Réf. EN-AS 91.5C.
Holmes J. D. (1975) – Prediction of the response of a cable-stayed bridge to turbulence. 4th Conference
of wind effects on buildings and structures, Heathrow, England, pp. 187-197.
Jin S., Lutes L.D. and Sarkani S. (1997) – Efficient simulation of Multidimensional Random Fields.
Journal of Engineering Mechanics, ASCE, Vol. 123, No.10, pp. 1082-1089.
Kaimal J. C., Wyngaard J., Izumi Y., et Cote O. R. (1972) – Spectral characteristics of surface-layer
turbulence. Journal of the Royal Meteorology Society, 98, pp. 563-589.
Larsen A. et al. (1992) – Aerodynamics of large Bridges. Larsen Editor, Balkema.
Liepmann H. W. (1952) – On the application of statistical concepts to the buffeting problem. Journal of
Aeronautical Science, Vol. 19, No. 12, pp. 793-800.
Lin Y. K. and Yang J. N. (1983) – Multimode bridge response to wind excitations. Journal of
Engineering Mechanics Division, ASCE, Vol. 100, No.4, pp. 657-672.
Loiseau H. and Szechenyi E. (1973) – Etude du comportement aéroélastique du tablier d’un pont à
haubans. Document ONERA, Châtillon.
Morgenthal G. (2000) – Comparison of Numerical Methods for Bridge-Deck Aerodynamics. Master of

Philosophy, University of Cambridge.
Namini A., Albrecht P., Bosch H. (1992) - Finite element-based flutter anaysis of cable-suspended
bridges. Journal of Structural Engineering, ASCE, Vol. 118, No. 118, No. 6, pp. 1509-1526.
182
Bibliographie
Patron-Solares A. (1998) – Modélisation numérique de la tenue au vent des ponts souples. Thèse de
doctorat, Ecole Nationale des Ponts et Chaussées.
Preumont A. (1990) – Vibrations aléatoires et analyse spectrale. Presses Polytechniques et

Universitaires Romandes.
Scanlan R.H. (1984) – Role of indicial functions in buffeting analysis of bridges. Journal of Structural
Engineering, ASCE, Vol. 110, n° 7, pp. 1433-1446.
Scanlan R. H. and Tomko J. J. (1971) – Airfoil and bridge flutter derivatives. Journal of Engineering
Mechanics, ASCE, Vol. 109, No. 2, pp. 586-603.
Scanlan R. H. , Béliveau J. G., and Budlong K. S. (1974) – Indicial aerodynamic functions for bridge
decks. Journal of Engineering Mechanics Division, ASCE, Vol. 100, No. 4, pp. 657-672.
Simiu E. (1974) – Wind spectra and dynamic alongwind response. Journal of Structural Engineering,
ASCE, Vol. 119, No. 2, pp. 365-382.
Shinozuka M. (1971) – Simulation of multivariate and multidimensional random processes. Journal of

the Acoustical Society of America, Vol. 49, No. 1, Part.2 , pp. 357-367.
Shinozuka M. , Jan C.M. (1972) – Digital simulation of random process and its applications. Journal of
Sounds and Vibrations, Vol. 25, No. 1, pp. 111-128.
Vaussy P. (1993) – Résistance des matériaux. Cours de l’Ecole Centrale de Nantes, Tome 2.
Virlogeux M. (1992) – Wind design and analysis for the Normandy Bridge. Aerodynamics of Large
Bridges, A. Larsen (ed.), pp.183-215.
Von Karman T. (1948) – Progress in the statistical theory of turbulence. Proc. Nat. Acad. Sci.,
Washington, D.C.,88, pp. 530-539.
183

Mémoire Thèse AdlF

Transféré par

Droits d'auteur :

Formats disponibles

Mémoire Thèse AdlF

Transféré par

Informations du document

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Mémoire Thèse AdlF

Transféré par

Droits d'auteur :

Formats disponibles

UNIVERSITE DE NANTES

Discipline : Sciences pour l’Ingénieur

Présentée et soutenue publiquement par

CALCUL DE LA REPONSE DYNAMIQUE

WIELGOSZ Christian Professeur, Université de Nantes Président

Directeur de thèse : JOUVE Pierre

Chapitre 1 Modélisation du vent 17

1.1 La couche limite atmosphérique 18

1.2 Caractéristiques du modèle statistique 18

1.3 Le vent moyen 20

1.4 Le vent turbulent 21

2.1 Forces aérodynamiques s’exerçant sur un corps mal profilé 30

2.2 Les forces aéroélastiques 43

2.3 Le détachement tourbillonnaire 52

Chapitre 3 Mesure des forces aérodynamiques 54

3.1 Mesure des coefficients stationnaires 55

3.2 Mesure de l'admittance aérodynamique 60

3.4 Méthode des oscillations forcées 75

4.1 Discrétisation en éléments finis 90

4.2 Calcul des modes propres de la structure 93

4.3 La méthode spectrale 95

4.5 La méthode temporelle 121

4.6 Introduction de simulations temporelles dans le calcul spectral 130

4.7 Conclusion 137

5.1 Présentation des maquettes 139

5.2 Caractérisation mécanique des maquettes 142

5.3 Caractérisation du vent de la soufflerie 148

5.4 Caractérisation aérodynamique des tabliers 156

5.5 Mesure des déplacements 158

5.6 Résultats 160

5.7 Conclusion 174

B Dimension caractéristique de tablier : en général sa largeur

Br Variable aléatoire, module de Ŷrn

FX , FY , M Z Forces de traînée, portance, moment moyennes dans un

FXt , FYt , FZt , M Xt , M tY , M tZ Forces et moments liés à la turbulence du vent dans un

FXa , FYa , FZa , M aX , M aY , M aZ Forces et moments aéroélastiques dans un repère local

FL , FLh , FM , FMh Fonctions indicielles dans le domaine spectral

f L , f Lh , f M , f Mh Fonctions indicielles dans le domaine temporel

g nr , ĝ nr Moyenne conditionnelle des termes u nr , û nr

H 1 , H 2 , H 3 , H 4 , A1 , A 2 , A 3 , A 4 Coefficients aéroélastiques de Scanlan

ĥ, ̂ Transformées de Fourier des déplacements vertical et

signal m et la composante k du signal p

U 10 m Vitesse moyenne de référence à 10m

u nr , û nr rième composante du signal tridimensionnel généré au point n

Yrn , Ŷrn Variable gaussienne centrée, composantes de u nr , û nr

z0 Rugosité d’un site

z ref Altitude de référence

 Xu ,  Yu ,  Mu ,  Xv ,  Yv ,  Mv ,  Xw ,  Yw ,  Mw « Racines carrées » des fonctions d’admittance

 Amortissement rapporté au critique

Yu , Yv , Yw Fonctions d’admittance aérodynamique de portance

Mu , Mv , Mw Fonctions d’admittance aérodynamique de moment

1 , 2 Pulsations des mouvements de flexion et torsion

dn Vecteur des déplacements d’un élément poutre dans la base nodale et

fV , fS Vecteurs forces extérieures volumique et surfacique

f te Vecteur des buffeting forces dans un repère local

f aee Vecteur des forces aéroélastiques dans un repère local

f aene Vecteur des forces aéroélastiques dans la base nodale et le repère

Fae Vecteur des forces aéroélastiques

Ft Vecteur des buffeting forces