Captulo 2

EL FLAMEO EN PUENTES DE GRAN VANO

2.5.1.1 Ensayo aerodinmico

En este caso el modelo reducido del tablero del puente est en la cmara de ensayo con
todos sus movimientos impedidos. El objetivo que se persigue es identificar los
coeficientes aerodinmicos que relacionan las fuerzas ejercidas por el flujo de aire sobre
el tablero con el producto de la presin dinmica por la longitud de referencia o su
cuadrado en el caso del momento. En la figura 2.5.2 se muestra el criterio de signos
considerado.

Figura 2.5.2. Esquema de una ensayo aerodinmico.

Los coeficientes aerodinmicos as definidos vienen dados por las expresiones

siguientes:

CD =

D
1
U 2 B
2

; CL =

L
1
U 2 B
2

; CM =

M
1
U 2 B 2
2

(2.5.1)

donde D es la fuerza de empuje por unidad de longitud que acta sobre el tablero, L es
la fuerza de elevacin por unidad de longitud y M es el momento segn el eje
longitudinal del puente por unidad de longitud de tablero. Estas fuerzas por unidad de
longitud, que dependen del ngulo de ataque como se muestra en la figura 2.5.3,
quedan adimensionalizadas por el producto de la presin dinmica por la anchura del
tablero B . En el ensayo, el flujo de aire debe generarse para distintos ngulos de ataque
sobre el tablero, .

50

Captulo 2

EL FLAMEO EN PUENTES DE GRAN VANO

16

14.7

12.9

11.8

10.8

9.63

8.37

7.22

6.09

4.93

3.79

2.68

1.67

D (gr)
L (gr)
M (gr*m*10)

1800
1600
1400
1200
1000
800
600
400
200
0
-200
-400

0.63

D (gr), L (gr), M (gr*m*10)

ngulo de ataque: 5 grados.

U (m/s)

ngulo de ataque: -5 grados.

500

-1000

15,7

14

12,2

11,3

10,3

8,15

7,05

6,02

4,32

3,15

1,67

-500

0,89

0
0

D (gr), L (gr), M (gr*m*10)

1000

D (gr)
L (gr)
M (gr*m*10)

-1500
-2000
U (m/s)

Figura 2.5.3. Fuerzas aerodinmicas medidas en un ensayo aerodinmico del tablero del puente del Great
Belt para un ngulo de ataque de 5 y -5.

Los valores de CL , CD y CM obtenidos para cada valor del ngulo de ataque , segn
(2.5.1) pueden ser ajustados a curvas segn se indica en la figura 2.5.4, que muestra los
resultados obtenidos para el tablero del puente del Great Belt segn Rey[R2], o los
presentados para el tablero del puente de Tacoma por Simiu[S1] y Scanlan en la figura
2.5.5.

51

Captulo 2

EL FLAMEO EN PUENTES DE GRAN VANO

Coeficiente aerodinmico de arrastre: Cd

0.40
0.35
0.30

Cd

0.25
0.20
0.15
0.10
0.05

10

15

20

10

15

20

-2

-4

-6

-1
0

-1
5

-1
9

0.00

ngulo de ataque ()

Coef. aerodinmico de elevacin: Cl

1.00
0.80

0.20

-0.60
-0.80
ngulo de ataque ()

52

-2

-4

-6

0
-1

-1

-0.20
-0.40

0.00

-1

Cl

0.60
0.40

Captulo 2

EL FLAMEO EN PUENTES DE GRAN VANO

Coef. aerodinmico del momento: Cm

0.08
0.06
0.04

20

15

10

-2

-4

-6

-1
0

-0.02

-1
5

0.00
-1
9

Cm

0.02

-0.04
-0.06
-0.08
-0.10
ngulo de ataque ()

Figura 2.5.4. Coeficientes aerodinmicos del tablero del puente del Great Belt segn Rey[R2].

Figura 2.5.5. Coeficientes aerodinmicos del nuevo puente sobre el estrecho de Tacoma.

Los ensayos aerodinmicos suelen llevarse a cabo en instalaciones de menores

dimensiones que los tneles de viento de capa lmite, como son los tneles
aerodinmicos. Estos tneles pueden considerarse como los herederos de los tneles que

53

Captulo 2

EL FLAMEO EN PUENTES DE GRAN VANO

tradicionalmente se han venido usando en el mbito de la ingeniera aeronutica. Se

caracterizan por que su cmara de ensayos es ms corta que en el caso de los tneles de
capa lmite. En ellos se puede generar flujo laminar o turbulento, y para este ltimo caso
se utilizan sistemas de rejilla. Con el fin de generar flujos con diferentes intensidades de
turbulencia se modifica la seccin de las barras que constituyen la rejilla, as como la
separacin entre las mismas. Otro factor que tambin tiene influencia es la distancia
existente entre la rejilla y el modelo a ensayar. En la figura 2.5.6 se muestra una imagen
del tnel de viento aerodinmico del que dispone el Grupo de Mecnica de Estructuras
de la Escola Tcnica Superior de Enxeeiros de Camios, Canais e Portos de la
Universidade da Corua. En la figura 2.5.7 se muestra el esquema de una rejilla que
permite la generacin de flujo turbulento.

Figura 2.5.6. Tnel de viento aerodinmico de la E.T.S.E.C.C.P. de la Universidade da Corua

Figura 2.5.7. Dispositivo generador de turbulencia.

54

Captulo 2

EL FLAMEO EN PUENTES DE GRAN VANO

2.5.1.2 Ensayo aeroelstico

En este tipo de ensayo, el modelo seccional debe simular lo ms exactamente posible las
caractersticas dinmicas del prototipo. En consecuencia, se debe realizar un anlisis
dimensional y de semejanza de parmetros que permita determinar la masa apropiada
para el modelo, as como la posicin y rigidez de los muelles elsticos que lo soportan.
En la figura 2.5.8 se muestra una fotografa de un modelo seccional del tablero del
puente del Great Belt ensayado en el tnel de viento del Grupo de Mecnica de
Estructuras de la Escola Tcnica Superior de Enxeeiros de Camios, Canais e Portos
de la Universidade da Corua.

Figura 2.5.8. Modelo seccional del Gran Belt preparado para realizar un ensayo aeroelstico

El objetivo final de este tipo de ensayos es determinar los coeficientes de flameo, flutter

derivatives en ingls, que relacionan las fuerzas aeroelsticas y los movimientos del
tablero. Como se ver ms adelante estas expresiones son de la forma siguiente:

D=

B
1

w
w

U 2 B KP1* + KP2* x + K 2 P3* x + K 2 P4* + KP5* + K 2 P6*

2
B
U
U
B
U

55

(2.5.2)

Captulo 2

L=

EL FLAMEO EN PUENTES DE GRAN VANO

B
1

w
w

U 2 B KH1* + KH 2* x + K 2 H 3* x + K 2 H 4* + KH 5* + K 2 H 6*
2
B
U
U
B
U

(2.5.3)

B
1

w
w

U 2 B 2 KA1* + KA2* x + K 2 A3* x + K 2 A4* + KA5* + K 2 A6*

2
B
U
U
B
U

(2.5.4)

M=

donde D, L y M son las tres componentes de la fuerza aeroelstica mientras que los
valores Hi*, Ai*, Pi* i=1,...,6 representan a los 18 coeficientes de flameo. Adems, K es
la frecuencia reducida, K = B / U ; siendo la frecuencia de la respuesta. En la figura
2.5.9 se representa el criterio de signos utilizado.

U
v

x
L

B
Figura 2.5.9 Fuerzas y movimientos en el ensayo aeroelstico.

En la figura 2.5.10 se adjuntan los grficos de algunos de los coeficientes de flameo

obtenidos para diversos puentes, as como los de un perfil alar, en los que U es la
velocidad del viento incidente, B el ancho del tablero y n la frecuencia de oscilacin.
Adems, en la figura 2.5.11 se muestra un posible esquema de sustentacin del modelo
seccional para realizar ensayos aeroelsticos.
Este tipo de ensayos fue propuesto por primera vez por Scanlan[S2] y Tomko en 1971.
Sin embargo, al emplear solamente dos grados de libertad, el giro segn el eje
longitudinal del tablero y el movimiento vertical, no era posible obtener todos los

56

Captulo 2

EL FLAMEO EN PUENTES DE GRAN VANO

coeficientes de flameo. En el ao 1977 Ibrahim[I1] y Mikulcik desarrollan un mtodo

numrico con el que se analizan los resultados de sistemas con respuestas que se
amortiguan exponencialmente. Basndose en este mtodo, Sarkar[S5] y Jones pusieron a
punto en 1994 una metodologa para obtener simultneamente los 18 coeficientes de
flameo. Asimismo, Jurado[J1], Len y Hernndez tambin han trabajado en la obtencin
del conjunto completo de coeficientes de flameo.

Figura 2.5.10. Ejemplos de coeficientes aeroelsticos

57

Captulo 2

EL FLAMEO EN PUENTES DE GRAN VANO

Figura 2.5.11. Esquema de sustentacin en un ensayo aeroelstico

Las ventajas que el Committee on Wind Effects de la ASCE[C1] destaca en relacin con
los ensayos sobre modelos seccionales son:

Bajo coste de los modelos seccionales.

La escala a la que se construye el modelo suele ser lo suficientemente grande como

para reproducir pequeos detalles de la geometra de la seccin con precisin.

Tanto las propiedades geomtricas y dinmicas del modelo seccional como el

sistema de soporte pueden modificarse fcilmente.

En lo que se refiere a los principales inconvenientes se citan los siguientes:

El modelo no permite reproducir el comportamiento real del puente en las

proximidades de estribos y torres. Sin embargo, en estos puntos lo habitual es que el

tablero est coaccionado o que las posibles distorsiones sean de pequea relevancia.
La metodologa hbrida tambin permite el estudio del fenmeno del bataneo; para ello
se definen las fuerzas por unidad de longitud del tablero debidas a la turbulencia del
viento como sigue:

58

Captulo 2

EL FLAMEO EN PUENTES DE GRAN VANO

Lb =

u dC

1
u
U 2 B 2CL y + L + CD z
U d
2
U

Db =
Mb =

u dC

1
U 2 B 2CD y + D
U d
2

uz

U

u dC

1
U 2 B 2 2CM y + M
U d
2

(2.5.5)

uz

U

donde el subndice b hace referencia al fenmeno del bataneo, representa el ngulo

de ataque del viento incidente, uy(x,t) y uz(x,t) representan la variacin respecto a la
velocidad media del viento incidente segn la componente transversal y vertical al
tablero respectivamente, mientras que los coeficientes CL , CD y CM se obtienen
experimentalmente, igual que en el caso del flameo. Por ltimo, se debe realizar el
anlisis estructural del puente bajo la accin conjunta de las fuerzas de flameo y de
bataneo.

2.5.2 Fase computacional. Clculo de la velocidad de flameo

Desde la dcada de los aos 60, el profesor Robert H. Scanlan desarroll junto a
numerosos investigadores como Tomko[S2], Sabzevari[S6], Gade[S7], Lin[S8], Huston[S9],
Kumarasena[K1] o Singh[S10], entre otros, toda una serie de investigaciones sobre la
modelizacin matemtica de las inestabilidades aeroelsticas en puentes de gran vano.
Scanlan puede ser considerado por tanto, el padre de la aeroelasticidad en el mbito de
la ingeniera civil. En la actualidad su testigo ha sido tomado por el profesor Nicholas
P. Jones[J2], actualmente investigador principal en este campo en la Johns Hopkins
University en Baltimore, Maryland, USA.
En el ao 1995 Jain[J3] desarroll junto a Scanlan y Jones una ampliacin de esta teora,
aadiendo, a los dos grados de libertad clsicos del flameo (w, x), el grado de libertad
correspondiente al movimiento horizontal en el sentido de la direccin del viento v, y su
primera derivada temporal v . La inclusin de este grado de libertad en el anlisis del
flameo tiene una importancia capital ya que los tableros de puentes de gran vano
presentan una gran flexibilidad en esa direccin debido a su gran longitud. Tanto es as,

59

Captulo 2

EL FLAMEO EN PUENTES DE GRAN VANO

que puentes como el del Great Belt, o el puente sobre el estrecho de Akashi, presentan
modos laterales entre sus primeros modos de vibracin, y con frecuencias menores que
los modos de torsin.
El anlisis multimodal consiste en expresar los movimientos de una estructura como
una combinacin lineal de un conjunto de los modos naturales de vibracin de la
misma. ste mtodo fue aplicado en el problema del flameo de puentes por Jain[J3],
quien expres los movimientos del tablero como una combinacin lineal de los m
primeros modos naturales de vibracin.
Seguidamente se presenta la formulacin del flameo multimodal siguiendo la notacin
establecida por Jurado[J4] y Hernndez. Esta notacin se adecua a la empleada
habitualmente en ingeniera de estructuras, y difiere de la empleada tradicionalmente
por otros autores como Scanlan[S7], Jones[J2] o Jain[J3], que siguen la notacin acuada en
el mbito de la ingeniera aeronutica. Autores como Dyrbye[D3] y Hansen emplean
notaciones diferentes de las dos antedichas que no aaden claridad. Igualmente, se pasa
revista a la metodologa matricial que permite formular el problema del flameo
multimodal ya empleada por investigadores como Astiz[A4], Katsuchi[K2], Jurado[J4] y
Hernndez, Mendes[M5] o Chen[C2] y Kareem.
El anlisis aeroelstico de un tablero de un puente es un anlisis dinmico en el que las
fuerzas actuantes son dependientes de los movimientos del propio tablero dentro del
flujo de viento.
Los modelos de clculo habitualmente utilizados en este tipo de estructuras, tanto
puentes atirantados como puentes colgantes, son modelos estructurales tridimensionales
de elementos barra, como el que se muestra en la figura 2.5.12 Cabe destacar que en
stos el tablero se modela como elementos barra tridimensionales a flexin, por lo que
los nudos poseen 6 grados de libertad, tres desplazamientos y tres giros, y 6 esfuerzos,
esfuerzo axil, dos esfuerzos cortantes, dos momentos flectores y un momento torsor. La
compatibilizacin de movimientos entre los movimientos del tablero y las pndolas y
cables principales, se realiza mediante la definicin de una serie de barras ficticias
transversales a dicho tablero, que conectan a ste con los extremos de dichas pndolas o
tirantes en el caso de puentes atirantados. Este tipo de modelo estructural puede ser
60

Captulo 2

EL FLAMEO EN PUENTES DE GRAN VANO

hallado en multitud de trabajos. Algunos de ellos son los realizados por Larsen[L7],
Cobo[C3], Namini[N2] y Aldbrecht, Astiz[A4] , Jurado[J5] et al., Semio[S11] y Mendes o
Miyata[M6] et al.

Figura 2.5.12. Modelo estructural tridimensional de un puente colgante

Las fuerzas aeroelsticas que ejerce el viento sobre el tablero se descomponen en tres
acciones por unidad de longitud:

Una fuerza de arrastre, drag, que produce una flexin lateral del tablero en el

plano XY. Esta fuerza en la notacin de Scanlan se denominaba como D, la notacin

aqu utilizada ser py.

Una fuerza de elevacin, lift, que produce una flexin vertical del tablero en el

plano XZ. Esta fuerza en la notacin de Scanlan se denominaba como L, siguiendo la

notacin aqu utilizada ser pz.

Un momento, moment, alrededor del eje X del tablero que produce movimientos

de tipo torsional a lo largo del mismo. Esta fuerza en la notacin de Scanlan se

denominaba como M, la notacin aqu utilizada ser mx.
Los movimientos de un nudo cualquiera de la estructura, siguiendo la notacin
tradicional de clculo de estructuras estarn definidos mediante el vector
61

Captulo 2

EL FLAMEO EN PUENTES DE GRAN VANO

u = ( u, v, w, x , y , z )

(2.5.6)

cuyas componentes representan los tres desplazamientos y tres giros de cada uno de los
nudos de la barra. Con esta notacin las expresiones de las fuerzas aeroelsticas, vase
la figura 2.5.13, resultan

py = D =

B
1

w
w

U 2 B KP1* + KP2* x + K 2 P3* x + K 2 P4* + KP5* + K 2 P6*

2
B
U
U
B
U

pz = L =

B
1

w
w

U 2 B KH1* + KH 2* x + K 2 H 3* x + K 2 H 4* + KH 5* + K 2 H 6* (2.5.8)
2
B
U
U
B
U

mx = M =

(2.5.7)

B
1

w
w

U 2 B 2 KA1* + KA2* x + K 2 A3* x + K 2 A4* + KA5* + K 2 A6* (2.5.9)

2
B
U
U
B
U

U
v

x
L

B
Figura 2.5.13. Movimientos y fuerzas aeroelsticas en el tablero

Dado que estas fuerzas aeroelsticas son una combinacin lineal de los movimientos y
velocidades de oscilacin, trabajando con un modelo estructural de puente con el tablero
discretizado en elementos barra, se puede definir de forma matricial un vector de cargas
actuando en cada elemento barra e de la forma

62

Captulo 2

EL FLAMEO EN PUENTES DE GRAN VANO

f ae = K aeu e + Caeu e

(2.5.10)

donde el subndice a hace referencia al origen aeroelstico de los trminos. Las fuerzas
en los dos nudos extremos de cada barra estn relacionadas con sus movimientos y
velocidades por las ecuaciones (2.5.7), (2.5.8), (2.5.9), llevndose cada uno de ellos la
mitad de la carga distribuida sobre el elemento. De este modo se puede establecer la
expresin matricial siguiente para cada uno de los nudos i=1,2 en los extremos de barra
f ai = K ai ui + Cai ui

(2.5.11)

que desarrollada queda

Fx
0
0

Fy
py
0
Fz le pz le 0

= =
M x 2 mx 2 0
My
0
0

0
0
Mz
0

0
l 0
+ e
2 0
0

*
4
2
*
6
2 *
6

pK P
pK H
pBK A
0
0

0 ui

0 vi
pBK P
0 wi
pBK H
+
0 xi
pB K A
0
0 yi

0
0 zi
0 0 ui

U 0 0 vi
U 0 0 wi

U 0 0 xi
0 0 yi

0 0 zi
0

*
6
2
*
4
2 *
4

pK P
pK H
pBK A
0
0

pBKP1* U
pBKH 5* U

pBKP5* U
pBKH1* U

pB 2 KP2*
pB 2 KH 2*

pB 2 KA5* U
0
0

pB 2 KA1* U
0
0

pB 3 KA2*
0
0

*
3
2
*
3
2
2 *
3

0
0
0
0
0
0

(2.5.12)

donde le es la longitud del elemento y p = U 2 2 es la presin dinmica del viento a

velocidad U. Con esta notacin matricial, el vector de fuerzas aeroelsticas para un
elemento barra del tablero segn la ecuacin (2.5.10) resulta
0 u1 Ca1 0 u1
f
K
f ae = a1 = K aeu e + Caeu e = a1
+

fa 2
0 K a 2 u 2 0 Ca 2 u 2

63

(2.5.13)

Captulo 2

EL FLAMEO EN PUENTES DE GRAN VANO

Ensamblando para todos los elementos que constituyen el tablero, se obtiene el vector
global de fuerzas aeroelsticas fa, as como las matrices globales aeroelsticas Ka y Ca,
que cumplen
f a = K a u + Ca u

(2.5.14)

siendo la dimensin de los vectores fa y u, el nmero de grados de libertad en el tablero

n=6N, siendo N el nmero de nudos en que se ha discretizado el tablero, y la dimensin
de las matrices aeroelsticas Ka y Ca de 6Nx6N, aunque observando sus desarrollos en
(2.5.12) slo poseen 3Nx3N trminos no nulos en general.
Con esta notacin matricial, el sistema de ecuaciones que gobierna el comportamiento
dinmico de un tablero sometido a unas fuerzas de origen aeroelstico es
Mu + Cu + Ku = f a = K a u + Ca u

(2.5.15)

siendo M, C y K las matrices de masa, amortiguamiento y rigidez del tablero. Pasando

los trminos de las fuerzas aeroelsticas al lado izquierdo de la igualdad resulta
Mu + ( C - Ca ) u + ( K - K a ) u = 0

(2.5.16)

que representa el movimiento de un sistema amortiguado en vibracin libre, en el que

las matrices de rigidez y amortiguamiento han sido modificadas por las fuerzas
aeroelsticas inducidas por el viento.
Aplicando el anlisis modal a los movimientos del tablero, la solucin del problema
definido en la ecuacin anterior puede ser aproximada como una combinacin lineal de
los m modos ms significativos
m

u = r qr = q
r =1

64

(2.5.17)

Captulo 2

EL FLAMEO EN PUENTES DE GRAN VANO

donde los modos de vibracin r pueden ser calculados para el puente completo, en
teora lineal o teniendo en cuenta la no linealidad geomtrica propia de los puentes de
gran vano, y slo incluyen los trminos correspondientes a los grados de libertad del
tablero. es la matriz modal de orden nxm, siendo n el nmero de grados de libertad
del tablero y m el nmero de modos de vibracin considerado, definida del siguiente
modo:
1,1

.
= .

n ,1

1,m

.
.

.
. n ,m

(2.5.18)

y cada elemento del vector q representa la participacin de cada uno de los modos de
vibracin en el vector de movimientos u.
Introduciendo la ecuacin (2.5.17) en (2.5.16) y premultiplicando por T se obtiene:
T Mq + T ( C - Ca ) q + T ( K - K a ) q = 0

(2.5.19)

Definiendo las siguientes matrices reducidas de dimensin mxm

M R = T M

C R = T ( C - Ca )

K R = T ( K - K a )

(2.5.20)

y teniendo en cuenta que los modos de vibracin se normalizan con respecto a la masa,
de modo que la matriz de masas reducida coincide con la matriz unitaria ( M R = I ), la
expresin (2.5.19) resulta

Iq + C R q + K R q = 0

(2.5.21)

donde las componentes de las matrices CR y KR, suponiendo un amortiguamiento tipo

Rayleigh[H7], verifican para i, j=1,...,m

65

Captulo 2

EL FLAMEO EN PUENTES DE GRAN VANO

n
n

CR ,ij = ri Ca ,rssj con i j

r =1
s =1

n
n

CR ,ii = 2ii ri Ca ,rssi

r =1
s =1

K R ,ij

= ri K a ,rssj con i j
r =1
s =1

(2.5.22)

n
n

K R ,ii = i2 ri K a , rssi
r =1
s =1

siendo i y i la frecuencia natural de vibracin y el coeficiente de amortiguamiento

clsico del modo i , respectivamente.
Para resolver el problema formulado en (2.5.21) para la respuesta aeroelstica del
tablero, se supone una oscilacin amortiguada en el tiempo de modo que

q(t) = we t

(2.5.23)

por lo que la ecuacin (2.5.21) puede ser reformulada como

( Iw + C w + K w ) e
2

=0

(2.5.24)

Utilizando la identidad siguiente

Iw + Iw = 0

(2.5.25)

y combinando ambas expresiones se obtiene

I 0 2 w C R

0 I w I
Definiendo

66

K R w t

e = 0
0 w

(2.5.26)

Captulo 2

EL FLAMEO EN PUENTES DE GRAN VANO

w
w =

-C
A= R
I

-K R

(2.5.27)

la ecuacin (2.5.26) puede escribirse del modo siguiente

( A I ) w e t = 0

(2.5.28)

que constituye un problema no lineal de valores y vectores propios.

La solucin de este problema de valores propios al que conduce el anlisis del
comportamiento aeroelstico de un tablero de puente afectado por un viento de
velocidad genrica U, proporciona en general 2m valores propios complejos de la forma

j = j + i j

j = 1, , m

(2.5.29)

y adems, si un valor propio es complejo, tambin es valor propio su complejo

conjugado

j = j i j

j = 1, , m

(2.5.30)

A cada par de valores propios complejos conjugados le corresponde un par de vectores

propios complejos conjugados entre s, de dimensin 2m

w = w R + iw I

w = w R iw I

(2.5.31)

Cada pareja de vectores propios conjugados junto con sus correspondientes valores
propios conjugados, define una forma de oscilacin amortiguada del tablero, que se
expresa como
u = we

jt

= ( w R + iw I ) e

( j + i j ) t

(2.5.32)

donde w es un vector de dimensin m que es parte de la definicin de w , ya indicada

anteriormente
67

Captulo 2

EL FLAMEO EN PUENTES DE GRAN VANO

w
w =

(2.5.33)

Si se separan las partes real e imaginaria de (2.5.32) se obtiene:

t
t
u = w R cos ( j t ) w I sen ( j t ) e j + i w R sen ( j t ) + w I cos ( j t ) e j

(2.5.34)

expresin que puede ser comparada con la que define las vibraciones libres
amortiguadas de un oscilador lineal con un grado de libertad:
u ( t ) = A cos (a t ) B sen (a t ) e at

(2.5.35)

en la que A y B son constantes que dependen de las condiciones iniciales de la

oscilacin, es la frecuencia natural, a = 1 a2 es la frecuencia de la respuesta
amortiguada y a es el amortiguamiento. Comparando estas dos ltimas expresiones,
(2.5.34) y (2.5.35), se puede establecer el siguiente paralelismo: La parte imaginaria de
los valores propios j juega el papel de la frecuencia amortiguada a en el oscilador
lineal, mientras que la parte real j se asocia al producto ( a ) , quedando as
definidos, en un tablero excitado con cargas aeroelsticas, los parmetros caractersticos
de la respuesta del sistema: la frecuencia de respuesta j-sima, aj = j y el coeficiente
de amortiguamiento j-simo, aj =

2j + 2j

Paso 1: Resolucin del problema de autovalores para una velocidad de viento U

La resolucin del problema de valores propios

( A I ) w = 0

en busca de las m

parejas de valores propios conjugados j = j i j ha de realizarse siguiendo un

proceso iterativo, puesto que en la construccin de la matriz A intervienen los
coeficientes de flameo, que son funcin de la frecuencia reducida del sistema
K = B a U . Pero la frecuencia asociada a cada modo aj j no se conoce hasta
que se ha resuelto el problema de valores propios. Un ejemplo de tcnica iterativa que

68

Captulo 2

EL FLAMEO EN PUENTES DE GRAN VANO

resuelve este problema se puede encontrar en Jurado[J6]. En la figura 2.5.14 se muestra

un diagrama de flujo en el que queda recogido el procedimiento iterativo propuesto. Los
significados de los valores de utilizados son:

jp : Valor de partida para la frecuencia de vibracin del modo j.

l : Parte imaginaria de uno cualquiera de los 2m valores propios obtenidos con
A ( jp ) I w = 0

md : Parte imaginaria del valor propio que hace mnima la diferencia jp l .

j : Frecuencia de respuesta del modo j obtenida despus de la convergencia. Este valor
cumple: A ( jp ) ( j j ) I w = 0 .

69

Captulo 2

EL FLAMEO EN PUENTES DE GRAN VANO

Figura 2.5.14. Diagrama de flujo para la obtencin de la respuesta aeroelstica.

70

Captulo 2

EL FLAMEO EN PUENTES DE GRAN VANO

Paso 2: Obtencin de la velocidad crtica de flameo

Atendiendo a la ecuacin (2.5.34), si para una determinada velocidad de viento U, son

negativas todas las partes reales j de cada uno de los valores propios complejos
obtenidos tras la resolucin del problema de valores propios definido en (2.5.28), el
amortiguamiento efectivo del tablero ser siempre positivo, y se obtiene una respuesta
aeroelstica que se atena exponencialmente con el tiempo, tendiendo a la estabilidad.
Por lo tanto la condicin crtica de flameo incipiente se produce para la velocidad ms
baja de viento con la que se obtiene un valor propio complejo con parte real nula, es
decir, la transicin entre amortiguamientos positivos (oscilaciones atenuadas) y
amortiguamientos negativos (oscilaciones de amplitud creciente). Por tanto,

j= 0

min

j
= min( aj ) = min
2
j + 2j

=0

(2.5.36)

identifica la condicin crtica de flameo, donde adems la parte imaginaria de dicho

valor propio j es la frecuencia de vibracin con la que se producen las oscilaciones en
la situacin de flameo.
De nuevo, el proceso de bsqueda de la velocidad de flameo se convierte en un proceso
iterativo, en el que se parte de una velocidad de viento U suficientemente baja, de modo
que j < 0 j=1,...,m. El proceso prosigue incrementando progresivamente la velocidad
U hasta que se verifica la condicin de flameo definida en (2.5.36). Dicha velocidad de
viento Uf producir la inestabilidad por flameo en el tablero. En la figura 2.5.15 se
adjunta el diagrama de flujo en el que se muestra el proceso a seguir para determinar la
velocidad de flameo elaborado por Jurado[J6].

71

Captulo 2

EL FLAMEO EN PUENTES DE GRAN VANO

Figura 2.5.15. Diagrama de flujo para la obtencin de la velocidad de flameo.

72

Captulo 2

EL FLAMEO EN PUENTES DE GRAN VANO

La metodologa en la que se combina el empleo de ensayos aeroelsticos en tneles de

viento aerodinmicos, con la formulacin propia de la aeroelasticidad computacional
que se ha desarrollado, es una metodologa hbrida, puesto que se debe recurrir a datos
experimentales para la identificacin de los coeficientes de flameo que permiten el
clculo numrico de los parmetros aeroelsticos de cada puente. Esta metodologa
presenta ventajas e inconvenientes con respecto a la completamente experimental
descrita con anterioridad. Entre las ventajas se encuentra el hecho de que las
instalaciones y el ensayo en s son ms baratos en trminos econmicos, pero se pierde
la clara visualizacin de las deformaciones del modelo que se observan en los ensayos
de los modelos de puentes completos en tneles de viento de capa lmite. Este problema
ha sido solucionado por Hernndez[H3][H4][H5][H6], Jurado y Mosquera al transferir los
resultados de la deformacin del puente obtenidos de los estudios aeroelsticos usando
una metodologa hbrida, a un modelo digital del puente completo. Vase la figura
2.5.16 en la que utilizando los resultados aeroelsticos que arroja la metodologa
hbrida, y por medio de un modelo digital de gran realismo del puente completo, se
visualiza el comportamiento dinmico y aeroelstico del mismo.

Figura 2.5.16. Imagen digital de la deformacin aeroelstica del puente sobre el estrecho de Messina.

73

Captulo 2

EL FLAMEO EN PUENTES DE GRAN VANO

La metodologa hbrida descrita en este apartado presenta una serie de ventajas que
merece la pena resaltar:

Requiere instalaciones de bajo coste en su fase experimental, con los beneficiosos

efectos que se comentaron en apartados anteriores.

La fase computacional se ve favorecida por el actual desarrollo de las capacidades

de los ordenadores digitales, lo que permite hacer cada vez modelos estructurales ms
refinados con anlisis multimodales en los que se considera un elevado nmero de
modos.

Permite un enfoque multidisciplinar incorporando otros mbitos cientficos como la

visualizacin avanzada, ya indicada anteriormente, o la aplicacin de tcnicas de

programacin distribuida u optimizacin numrica.
Por todos estos aspectos, esta metodologa ha sido la utilizada en esta investigacin y
algunos de los nuevos enfoques indicados han sido implementados satisfactoriamente.

2.6

REFERENCIAS

[A1] Astiz M. A. [1996] Wind Related Behaviour of Alternative Suspension Systems.

15th IABSE Congress Report, Copenhagen, pp. 1079-1090.
[A2] Anagnostopoulos P. [1997] Computer-aided Flow Visualization and Vorticity
Balance in the laminar Wake of a Circular Cilinder. Journal of Fluids and Structures,
Vol 11, pp 33-72.
[A3] Alan G. Davenport Wind Engineering Group, The Boundary Layer Wind Tunnel
Laboratory and The University of Western Ontario, London, Ontario, Canada. [2004]

http://blow.blwtl.uwo.ca/blwtlhom.htm
[A4] Astiz M. A. [1999] Flutter Stability of Very Long Suspension Bridges. Journal of
Bridge Engineering, Vol. 3, N 3, August, pp.132-139.

74

Captulo 2

EL FLAMEO EN PUENTES DE GRAN VANO

[B1] Bisplinghoff R. L., Ashley H., Halfman R. L. [1996] Aeroelasticity. Dover

Publications, Inc. New York.
[B2] Braza M., Chassaing P., Minh H. H. [1986] Numerical Study and Physical
Analysis of the Pressure and Velocity fields in the near Wake of a Circular Cylinder.
Journal of Fluid Mechanics, Vol 165, pp. 79-130.
[C1] Committee on Wind Effects of the Committee on Dynamic Effects of the
Structural Division of the American Society of Civil Engineers [1987] Wind Loading
and Wind-Induced Structural Response. American Society of Civil Engineers (ASCE).
[C2] Chen X., Kareem A. [2002] Advanced Analysis of Coupled Buffeting Response of
Bridges: a Complex Modal Decomposition Approach. Probabilistic Engineering
Mechanics. Vol 17, N 2, pp. 201-213.
[C3] Cobo del Arco D. [1998] An analysis of wind stability. Improvements to the
response of suspension bridges. Tesis doctoral de la Universitat Politcnica de
Catalunya. Escola Tcnica Superior DEnginyers de Camins, Canals i Ports. (Director:
Prof. ngel Aparicio)
[D1] Dawes W. N. [1993] Simulating Unsteady Turbomachinery Flows on Unstructed
Meshes wich adap both in Time and Space. International Gas Turbine and Aeroengine
Congress and Exposition, Cincinnati, Ohio. 1993.
[D2] Davenport A. G., Isyumov N. [1968] The application of the boundary-layer wind
tunnel to the prediction of wind loading. Proceedings of the International Research
Seminar on Wind Effects on Buildings and Structures. Univ. of Toronto Press, Toronto,
pp.201-230.
[D3] Dyrbye C., Hansen S. O. [1996] Wind Loads on Structures. John Wiley & Sons.
[F1] Fujino Y., Yoshida Y. [2002] Wind-Induced Vibration and Control of TransTokyo Bay Crossing Bridge. Journal of Structural Engineering. Vol. 128, N 8, August,
pp. 1012-1025.
75

Captulo 2

EL FLAMEO EN PUENTES DE GRAN VANO

[F2] Frandsen J. B., McRobie F. A. [2004] Bridge Deck Aeroelasticity using FE.

http://www.fsi.lsu.edu/~frandsen/approach.html
[F3] Farquharson F. B., Smith F. C. Vincent G. S. [1949-1954] Aerodynamic Stability
of Suspension Bridges with special reference to the Tacoma Narrows Bridge. University
of Washington Engineering Experiment Station. Bulletin n 116, Parts I to V.
[F4] Fung Y. C. [1955] An Introduction to the Theory of Eeroelasticity. Dover
Publication, Inc. New York.
[F5] Frandsen J. B., McRobie F. A. [1999] Computational Aeroelastic Modelling to
guide long-span Bridge cross-section Design. 10th International Conference on Wind
Engineering. Copenhagen. 1999.
[F6] Force Technology home page. [2004] http://www.force.dk/gb/default.htm
[H1] Holmes J. D., [2001] Wind Loading of Structures. Spon Press. London.
[H2]

Honshu-Shikoku

Bridge

Authority

home

page.

[2004]

http://www.hsba.go.jp/bridge/akdesign.htm
[H3] Hernndez S. [2001] Virtual Visualization of Aeroelastic Behaviour of Suspension
Bridges. IABSE Conference in Cable Supported Bridges. Challenging Technical Limits.
Seoul, 12-14 Junio 2001.
[H4] Hernndez S., Jurado J. A., Mosquera A. [2002] Multidisciplinary and
multiobjective aeroelastic design of long span bridges. In-Vento 2002. International
Workshop on Long Span Bridge Aerodynamics. Milan, September 15-18, 2002.
[H5] Hernndez S. Jurado J. ., Mosquera A. [2003] Virtual Wind Tunnels for
Aeroelastic Design of Cable Supported Bridges. 11th International Conference on Wind
Engineering. Lubbock, Texas, USA. 2-5 Junio 2003.

76

Captulo 2

EL FLAMEO EN PUENTES DE GRAN VANO

[H6] Hernndez S., Jurado J. . [2004] Computer Animations of the Aeroelastic

Behaviour of the Proponed Messina Bridge. In-Vento 2004. 5 Convegno Nazionale di
Ingegneria del Vento. Reggio Calabria. 21-23 Junio 2004.
[H7] Hart G. C., Wong K. [2000] Structural Dynamics for Structural Engineers. John
Wiley & Sons, Inc.
[I1] Ibrahim S. R., Mikulcik E. C. [1977] A method for the direct identification of
vibration parameters from the free response. The Shock and Vibration Bull. 47, (Part
4).
[J1] Jurado J.., Len A., Hernndez S. [2005] Wind Tunnel Control Software for
Identification of Flutter Derivatives on Bridge Section Tests. The fourth European &
African Conference on Wind Engineering. Prague, 11-15 July. Paper #108.
[J2] Jones N. P., Jain A., Katsuchi H. [1998] Advances (and challenges) in the
prediction of long-span bridges response to wind. Bridge Aerodynamics. Edit. By
Larsen A. & Soren E. A. A. Balkema/Rotterdam/Brookfield.
[J3] Jain A. [1996] Multi-Mode Aeroelastic and Aerodynamic Analysis of Long-Span
Bridges. A dissertation submitted to The Johns Hopkins University in conformity with
the requirements for the degree of Doctor of Philosophy.
[J4] Jurado J. A., Hernndez S. [1999] Design of Ultra-Long Span Bridges With
Aeroelastic Constraints. Optimization in Industry II, Junio, Banff, (Alberta) Canada.
[J5] Jurado J. ., Mosquera A., Nieto F., Hernndez S. [2004] Configuraciones
ptimas de Cables en Puentes Atirantados segn Anlisis de Sensibilidad de su
Comportamiento Aeroelstico. Congresso de Mtodos Computacionais em Engenharia.
Lisboa, Portugal. 31 Mayo-2 Junio 2004.
[J6] Jurado J. A. [2001] Anlisis Aeroelstico y de Sensibilidad del Fenmeno de
Flameo en Puentes Soportados por Cables. Tesis doctoral de la Universidade da
Corua. Escola Tcnica Superior de Enxeeiros de Camios, Canais e Portos.
77

Captulo 2

EL FLAMEO EN PUENTES DE GRAN VANO

[K1] Kumarasena T. [1989] Wind response prediction of long-span bridges. PhD thesis,
The Johns Hopkins University, Baltimore, MD, USA.
[K2] Katsuchi H., Saeki S., Miyata T., Sato H. [1998] Analytical assessment in wind
resistant design of long-span bridges in Japan. Proceedings of the international
symposium on advances in bridge aerodynamics. Copenhagen (Denmark), 10-13 May.
Pp87-98.
[L1] Larose G. L., Franck N. [1997] Early Wind Engineering experiments in Denmark.
Journal of Wind Engineering an Industrial Aerodynamics, 72: 493-9
[L2] Lecointe, Y., Piquet J. [1989] Flow Structure in the Wake of an Oscillating
Cylinder. Journal of Fluid Engineering. Vol 111, pp. 139-148.
[L3] Larsen A. [1998] Advances in Aeroelastic Analyses of Suspension and Cablestayed Bridges. Journal of Wind Engineering and Industrial Aerodynamics. Vol 74-76,
pp. 73-90.
[L4] Larsen A., Walther J. H. [1998] Discrete Vortex Simulation of Flow around Five
generic Bridge Deck Sections. Journal of Wind Engineering and Industrial
Aerodynamics, Vol 77-78, pp. 591-602.
[L5] Larsen A., Walther J. H. [1996] Aeroelastic Analysis of Bridge Girder Sections
based on Discrete Vortex Simulations. 2nd International Conference on Computational
Wind Engineering. Colorado, USA. 1996
[L6] Lopes A. V., Cunha A., Simoes L. M. C. [2004] Modelo Computacional de
Anlise Aeroelstica das Condioes de Utilizaao de Estructuras Esbeltas. Congresso
de Metodos Computacionais em Engenharia. 31 Mayo-2 Junio. Lisboa, Portugal.
[L7] Larsen A. [1993] Aerodynamic aspects of the final design of the 1624 m
suspension bridge across the Great Belt. Journal of Wind Engineering and Industrial
Aerodynamics, 48, pp.261-285.

78

Captulo 2

EL FLAMEO EN PUENTES DE GRAN VANO

[M1] Meseguer J., Sanz A., Perales J. M., Pindado S. [2001] Aerodinmica Civil.
Cargas de Viento en las Edificaciones. McGraw-Hill/Interamericana de Espaa, S. A.
U. Madrid.
[M2] Morgenthal G. [2000] Fluid-Structure Interaction in Bluff-Body Aerodynamics
and Long-Span Bridge Design: Phenomena and Methods. Technical Report. University
of Cambrige, Department of Engineering.
[M3] Morgenthal G. [2000] Comparison of Numerical Methods for Bridge-Deck
Aerodynamics. Thesis submitted to the University of Cambridge in partial fulfilment of
the requirements for the degree of Master of Philosophy.
[M4] Mosquera A., [2003] Sensibilidad Analtica de la Vibracin Libre y del Flameo en
Puentes Colgantes de Gran Vano en Servicio y en Fase Constructiva. Tesis doctoral de
la Universidade da Corua. Escola Tcnica Superior de Enxeeiros de Camios, Canais
e Portos.
[M5] Mendes P. [1997] Estabilidade Aerodinmica de Pontes. Clculo da Velocidade
Crtica de Flutter. Relatorio ICIST, DT N19/97.
[M6] Miyata T., Yamada H., Katsuchi H. [2003] Comparative Analysis of Messina
Bridge International Benchmark study. 11th International Conference on Wind
Engineering. Lubbock, Texas, USA. 2-5 Junio 2003.
[N1] Nieto F., Jurado J. ., Hernndez S. [2002] Ensayos Aeroelsticos de Modelos
Seccionales de Tableros de Puentes. Universidade da Corua. Escola Tcnica Superior
de Enxeeiros de Camios, Canais e Portos. Ed.: Hernndez S.
[N2] Namini A., Aldbrecht P. [1992] Finite Element-Based Flutter Analysis of Cable
Suspended Bridges. Journal of Structural Engineering. Vol. 118, N 6, June, pp.15091526.

79

Captulo 2

EL FLAMEO EN PUENTES DE GRAN VANO

[P1] Preidikman S., Mook D. T. [1998] On the Development of a Passive-Damping

System for Wind-Excited Oscillations fo Long-Span Bridges. Journal of Wind
Engineering and Industrial Aerodynamics. Vol 77, pp. 443-456.
[R1] Rodi W. [1997] Comparison of LES and RANS calculations of the flow around
bluff bodies. Journal of Wind Engineering and Industrial Aerodynamics. Vol 69-71, pp
55-75.
[R2] Rey J. [2000] Puesta a Punto de un Tnel de Viento y Realizacin de Ensayos
Aerodinmicos de Tableros de Puentes. Proyecto Tcnico. Universidade da Corua.
Escola Tcnica Superior de Enxeeiros de Camios, Canais e Portos.
[S1] Simiu E., Scanlan R. H. [1996] Wind Effects on Structures: Fundamentals and
applications to Design. (3rd Ed.). John Wiley & Sons, Inc.. New York.
[S2] Scanlan R. H., Tomko J. J. [1971] Airfoil and Bridge Deck Flutter Derivatives. J.
Eng. Mech. Div., ASCE, 97, N EM6, Diciembre, pp1717-1737.
[S3] Son J. S., Hanratty T. J. [1969] Numerical Solution for the Flow around a Cylinder
at Reynolds Number of 40, 200, 500. Journal of Fluid Mechanics. Vol 35, pp. 369-386.
[S4] Steggel N., Rockliff N. [1997] Simulations of the Effects of Body Shape on Lock-in
Characteristics in Pulsating Flow by the Discrete Vortex Method. Journal of Wind
Engineering and Industrial Aerodynamics. Vol 69-71, pp. 317-329.
[S5] Sarkar P. P., Jones N. P., Scanlan R. H. [1994] Identification of aeroelastic
parameter of flexible bridges. J. Eng. Mech. Vol. 120, N 8, August, pp.1718-1742.
[S6] Scanlan R. H., Sabzevari A. [1969] Experimental aerodynamic coefficients in the
analytical study of suspension bridge flutter. J. Eng. Mech. Vol. 11, N 3, pp.234-242.
[S7] Scanlan R. H., Gade R. H. [1977] Motion of suspended bridge spans under gusty
wind. J. Str. Div., ASCE. Vol.103, N 9, pp.1867-1883.

80

Captulo 2

EL FLAMEO EN PUENTES DE GRAN VANO

[S8] Scanlan R. H., Lin W. H. [1978] Effects of turbulence on bridge flutter derivatives.
J. Eng. Mech., Vol. 104, N 4, pp.719-733.
[S9] Scanlan R. H., Huston D. [1985] Sensitivity of Bridge Decks to Turbulent Wind.
Proc., Asia Pacific Sym. On Wind Eng., University of Roorkee, India, pp.167-172.
[S10] Singh L. [1997] Experimental Determination of Aeroelastic and Aerodynamic
Parameters of Long-Span Bridges. A dissertation submitted to The Johns Hopkins
University in conformity with the requirements for the degree of Doctor of Philosopgy.
Baltimore, Maryland, USA.
[S11] Semio C., Mendes P. [1998] Avaliazao das condizoes de estabilidade
aerodinmica em pontes. Jornadas portuguesas de engenharia de estructuras. Lisboa,
Portugal.
[T1] Theodorsen T. [1935] General Theory of Aerodynamic Instability and the
Mechanism of Flutter. NACA Rept. 496.
[T2] Theodorsen T., Garrick I. E. [1940] Mechanism of Flutter, a Theoretical and
Experimental Investigation of the Flutter Problem. NACA. Rept. 685.
[T3] Tukiyyah G. M., Reed D. A. [1998] 3D Discrete Vortex Methods and Their
Parallel Implementation. Bridge Aerodynamics, Larsen A. & Sren E. (Ed.). Balkena,
Rotterdam.
[U1] University of Washington. University of Washington Libraries Digital Collections.

www.lib.washington.edu/specialcoll/tnb/default.html
[V1] Vincent G. S. [1958] Golden Gate Bridge Vibrational Studies. Journal of the
Structural Division, ASCE, paper 1817.
[W1] Wardlaw R. L. [1992] The Improvement of Aerodynamic Performance.
Aerodynamics of Large Bridges, Larsen A. (Ed.). Balkema, Rotterdam.

81

Captulo 2

EL FLAMEO EN PUENTES DE GRAN VANO

[W2] Washington State Departament of Transportation home page. www.wsdot.wa.gob

[W3] Walther J. H., [1998] Discrete Vortex Methods in Bridge Aerodynamics and
prospects for parallel computing techniques. Bridge Aerodynamics, Larsen A. & Sren
E. (Ed.). Balkena, Rotterdam.
[Z1] Zhou Z., Chen A., Xiang H. [2003] Identification of Aeroelastic Parameter of
Flexible Bridge Decks by Random Discrete Vortex Method. 11th International
Conference on Wind Engineering. 2-5 June. Lubbock, Texas. USA.

82