Análisis del colapso del puente Tacoma Narrows

utilizando ecuaciones diferenciales

Jose M. Bedoya Salazar, Daniel Garcı́a Botero,
Santiago Gómez Echeverry.
Marzo de 2020

1 El Problema 1.5 Limitaciones

1.1 Planteamiento del Problema 2 Estado del Arte

1.2 Formulación del Problema 2.1 Antecedentes de la investi-
1.3 Objetivos gación
• Comprender de manera detallada el uso y apli- El Puente de Tacoma Narrows fue un puente col-
cación de las ecuaciones diferenciales sobre colapso gante de 1600 metros de longitud, con una distan-
del puente Tacoma Narrows. cia entre soportes de 850 m (el tercero más grande
• Identificar las diferentes formas para solucionar del mundo en la época en que fue construido). El
una ecuación no lineal de segundo orden. puente fue construido en 1940 y abierto al público
• Explicar por medio de modelos matemáticos, en Julio de ese mismo año. Después de 4 meses en
fenómenos fı́sicos a través de ecuaciones diferen- los cuales se evidenciaron las grandes vibraciones
ciales ordinarias. verticales ocasionadas por el viento en la capa de
asfalto, este colapsó a principios de Noviembre a
causa de las grandes magnitudes de las oscilaciones
generadas por una fuerte tormenta [1].
1.4 Justificación de la investigación
Este proyecto de aula nos da a conocer los difer-
entes usos que se pueden obtener con las ecua-
ciones diferenciales, tanto como en los fenómenos
fı́sicos como matemáticos; conociendo los difer-
entes métodos usados en clase para ası́ comprender
el método empleado en la ecuación diferencial para
explicar el colapso del puente Tacoma Narrows.
Además de comprender por qué se usó un método
matemático en un fenómeno fı́sico, dado que los
ingenieros no lograban creer que la causa del co-
lapso del puente era únicamente por el fenómeno Figura 1. Fotografı́as del puente Tacoma Narrows
de resonancia, por lo que se simplificó el modelo en 1940 antes y después de su colapso.
que conduce a una ecuación diferencial ordinaria Tomadas de [2] y [3].
no lineal; y ası́ dar a conocer la solución por medio Von Kármán fue el ingeniero encargado de deter-
de una ecuación diferencial. minar las causas del colapso del puente.

1
 orden que convenga. Una ecuación diferencial or-
dinaria de orden n(n ∈ N) en forma implı́cita es
una expresión del tipo:
F (x, y(x), y 0 (x), ..., y (n) (x)) = 0
Mientras que expresada en forma explı́cita adopta
la forma:
y (n) (x) = f (x, y(x), y 0 (x), ..., y (n−1) (x)) = 0
Y donde las solución general de estas ecuaciones
Figura 2. Fotografı́as del puente Tacoma Narrows es una familia de curvas cuya expresión analı́tica
en 1950 y en la actualidad, respectivamente. depende de n constante arbitraria; escribiéndose:
Tomadas de [5] y [6].
y = (y(x; C1 ), ..., Cn )

2.2 Bases teóricas en forma explı́cita, o

Se entiende como ecuación diferencial toda f (x, y; C1 , ..., Cn ) = 0

ecuación que relaciona los valores de una función de forma implı́cita. Dichas constantes se deter-
con los valores de sus derivadas minan una vez dado un conjunto de condiciones
iniciales:
Definición: Se denomina ecuación diferencial
(n−1)
(ED) a la ecuación que contiene derivadas de x0 , y(x0 ) = y0 , y 0 (x0 ) = y00 , ..., y (n−1) (x0 ) = y0
una o más variables dependientes respecto a una
o más variables independientes[1]. Por lo tanto, toda ecuación diferencial de orden n
(expresada de forma normal)
Se denomina ecuación diferencial ordinaria (EDO)
a una ecuación diferencial en la que la función de- y (n) (x) = f (x, y(x), y 0 (x), ..., y (n−1) (x)) = 0
pende solo de una variable. En este último caso, es equivalente a un sistema de n ecuaciones difer-
se dice que la ecuación está expresada en forma enciales de primer orden
normal o explı́cita si la derivada de orden supe-
rior aparece despejada como función de todos los dy
= y1 (x)
demás ingredientes de la ecuación. En caso con- dx
trario se dice que la ecuación está expresada en dy1
= y2 (x)
forma implı́cita. dx
..
.
2.2.1 Ecuación diferencial de orden supe- dyn−2
= yn−1 (x)
rior dx
dyn−1
Las ecuaciones diferenciales, a su vez, se pueden = f (x, y1 (x), ..., yn−1 (x))
dx
clasificar según su orden, donde este, es el orden
de la mayor derivada en la ecuación. Por ejemplo: donde
y(x), y1 (x)..., yn−1 (x)
d2 y dy x
dx2
+ 3( dx ) − 20y = e son funciones desconocidas.
De este modo, resolver la ecuación diferencial de
En la cual se puede ver que es una ecuación de orden n equivale a resolver el sistema anterior, esto
segundo orden debido a que el exponente de la es, a buscar las n funciones:
derivada es 2.
Si y = y(x), indica una función derivable hasta el y(x), y1 (x)..., yn−1 (x)

2
Donde cada una de ellas tendrá una constante ar- 2.2.2 El fenómeno de la resonancia
bitraria y el total de estas se determinará a partir
Se conoce como resonancia, el estado en el que en-
de las condiciones iniciales.
tra un sistema elástico cuando este tiende a vibrar
• Linealidad de las ecuaciones diferenciales:
u oscilar con una máxima amplitud [7]. La man-
Se dice que una ecuación diferencial de n-ésimo
era más sencilla de explicar este fenómeno es me-
orden tal como:
diante un sistema sencillo masa-resorte, donde m
0 (n)
F (x, y(x), y (x), ..., y (x)) = 0 es la masa sujeta al resorte con constante elástica
K; por ende, para describir la dinámica de este
es lineal si F es lineal en sistema se tiene como partida la segunda ley de
Newton, donde
y, y 0 , ..., y (n) d2 y
dt2
Significando ası́ que una EDO de n-ésimo orden es
lineal cuando: es la aceleración.
d2 y
F (x, y(x), y 0 (x), ..., y (n) (x)) = 0 m = −Ky (1)
dt2
es
Y solucionando ası́ esta ecuación diferencial an-
an (x)y n + an−1 (x)y (n−1) + ...
terior, se propone como solución en función del
0
... + a1 (x)y + a0 (x)y − g(x) = 0 tiempo:
Por lo tanto, una EDO no lineal es simplemente y(t) = Acosw0 t (2)
una ED donde las funciones no lineales de la vari- Es decir, un movimiento armónico simple.
able dependiente o de sus derivadas como sen(y) • Frecuencia natural del sistema: Si sustitu-
o ey y no pueden aparecer en una ecuación lineal imos esta solución en (1) y despejando w0 , se tiene
Por ejemplo: que la frecuencia angular con que el sistema se
y 0 + sen(y) = 0 moverá es: r
Donde hay una función no lineal de y K
w0 = (3)
m
d4 y Viéndose ası́ que esta frecuencia angular no de-
+ y2 = 0
dx4 pende de la amplitud si no de la constante elástica
Donde el exponente de y es diferente de 1. y de la masa, por tanto, este sistema tiene una sola
Las ecuaciones de n-ésimo orden lineales, tienen frecuencia que si acoge inmediatamente el sistema
la propiedad de que la combinación lineal de sus este deja de oscilar libremente, denominándose ası́
soluciones también es una solución. Mientras “frecuencia natural del sistema”.
que las ecuaciones no lineales no cuentan con • Oscilaciones forzadas: Sobre el sistema an-
esta propiedad de superposición debido a que es- teriormente montado, se aplica sobre la masa otra
tas ecuaciones no lineales desafı́an virtualmente fuerza F , la cual será periódica y con una amplitud
la solución con métodos analı́ticos [1]; sin em- F , frecuencia angular w y en dirección paralela a
bargo, siempre es posible analizar de modo cualita- la del eje del resorte. Por sumatoria de fuerzas,
tivo y numérico una ecuación diferencial no lineal. tenemos:
Cuando hablamos de ecuaciones diferenciales no
d2 y
lineales de orden superior, hablamos de, quizás, m 2 = −Ky + F cos(wt) (4)
las más importante si se quiere modelar un sis- dt
tema fı́sico, debido a que medida que se ajusta y al igual que en el caso anterior, a esta ecuación
a un modelo matemático al sistema, también au- diferencial se propone la solución
menta la probabilidad de que este sea un sistema
no lineal. y = Acos(wt)

3
y sustituyendo, tenemos: Figura 3. Vórtices de Von Kármán.
Tomada de [9].
−mAw2 cos(wt) = −KAcos(wt) + F cos(wt) (5)

Y si despejamos la amplitud A de la oscilación:

F
A= (6)
K − mw2
Pero de acuerdo a (3), K = mw02 ası́ que susti-
tuyendo este valor en la anterior relación se obtiene
finalmente que
F0
A= Figura 4. Otra vista de los vórtices de Von
m(w02 − w2 )
Kármán.
Por lo anterior, si la frecuencia de la fuerza externa Tomada [10].
w, igual a la frecuencia natural del sistema w0 , la
amplitud tiende a infinito; o en términos fı́sicos,
el sistema va a vibrar a una amplitud extremada-
mente grande. • Oscilación forzada debido a vórtices:
• Vórtices de Von Kármán: Un vórtice de Von Este fenómeno se puede ver en numerosas estruc-
Karman es un remolino causado por una sepa- turas que son expuestas a corrientes de viento,
ración no estacionaria del flujo de un fluido sobre como por ejemplo a los vórtices de Von Kármán;
un cuerpo y donde su frecuencia depende de la ve- donde, como se mencionó anteriormente, su fre-
locidad de flujo [8]. Este fenómeno recibe su nom- cuencia depende de la magnitud de su velocidad
bre en honor al ingeniero y experto en dinámica del viento. Por consiguiente, si esta frecuencia se
de fluidos Theodore Von Kármán; pionero en su encuentra cerca a la frecuencia de flexión o nat-
investigación. ural de la estructura, sobre la cual está siendo
En la naturaleza es posible presenciar este sometida la corriente de viento, puede llegar a exis-
fenómeno en algunas ocasiones, como por ejem- tir el riesgo de que la estructura entre en un estado
plo, cuando el viento sopla sobre el mar este se de resonancia, donde el desplazamiento o amplitud
perturba al topar con las islas, montañas, o hasta solamente estarı́a limitado por el amortiguamiento
estructuras grandes; dejando ası́ una serie de este- de la estructura.
las de igual forma que los vórtices de Von Kármán,
las cuales resultan visibles gracias a las nubes. Lo
anterior mencionado se puede ver en las siguientes 2.2.3 Definición de términos básicos
figuras:
En términos fı́sicos, una “fuerza externa” hace ref-
erencia a las fuerzas que realiza un cuerpo sobre
otro diferente, es decir, representan la interacción
del sistema con el exterior del mismo; o en otros
términos, es la capacidad de un cuerpo de modi-
ficar la velocidad de desplazamiento, movimiento
y/o estructura de otro cuerpo diferente.
El término “flexión” en términos estructurales e
ingenieriles se refiere a al tipo de deformación que
presenta un elemento estructural alargado en una
dirección perpendicular a su eje longitudinal.

4
2.2.4 Sistema de Hipótesis [5] Seattle Daily Journal Of Commerce,
”This Century’s Top Ten Construction
“El acercamiento de la frecuencia de una fuerza
Projects”, 1999. [En lı́nea]. Disponible:
aplicada a una estructura a la frecuencia natural
https://www.djc.com/special/century/10060860.
de la misma estructura, producirá daños, defor-
htm. [Acceso: 29- Mar- 2020].
mación y/o rompimiento de la misma”
[6] M. Cruz, ”Puente de Tacoma Narrows”,
“A medida que se ajusta más un modelo
Makinando, 2015. [En lı́nea]. Disponible:
matemático a uno fı́sico, este presentará un com-
https://makinandovelez.wordpress.com/2015/
portamiento no lineal”.
05/14/puente-de-tacoma-narrows/. [Acceso: 29-
Mar- 2020].
3 Bibliografı́a [7] J. Peralta, P. Reyes López and A. Godı́nez
Muñoz, ”El fenómeno de la resonancia”, Latin-
[1] D. Zill y W. Wright, Ecuaciones diferenciales American Journal of Physics Education, vol. 3,
con problemas con valores en la frontera, 8th ed. no. 3, pp. 612, 618, 2009. [Acceso: 29 March
Distrito Federal: CENGAGE Learning, 2013, pp. 2020].
13-16. [8] G. Sibera Quilez, ”Estudio de la calle de von
[2] J. Zasky, ”Tacoma Narrows Bridge Collapse Kármán en un cilindro mediante PIV”, Ingeniero
- Failure magazine”, Failuremag.com, 2000. [En Superior Aeronáutico, Universitat Politècnica de
lı́nea]. Disponible: http://failuremag.com/article/ Catalunya, 2013.
tacoma-narrows-bridge-collapse. [Acceso: 29- [9] ”Vórtices de Von Karman: las nubes-
Mar- 2020]. remolinos más bonitas están en Canarias”, El
[3] J. Velasco, ”El puente de Tacoma, Tiempo Hoy, 2020. [En lı́nea]. Disponible en:
el viento y la resonancia”, Hipertex- https://www.eltiempohoy.es/elcielo/meteorologia/
tual, 2012. [En lı́nea]. Disponible: Vortices-Von-Karman-nubes-remolinos-Canarias-
https://hipertextual.com/2012/02/puente-de- fenomeno-meteorologia02231775126.html. [Ac-
tacoma-viento- resonancia. [Acceso: 29- Mar- ceso: 30- Mar- 2020].
2020]. [10] ”Dos vistas de vórtices de von Kármán”,
[4] A. Baldomero Valiño, ”¿SON EVITABLES Tiempo.com - Meteored, 2020. [Online].
LOS DESASTRES EN LAS OBRAS? La im- Disponible en: https://www.tiempo.com/ram/
portancia de las Ecuaciones Diferenciales en la 357642/dos-vistas-vortices-von-karman/. [Acceso:
Mecánica de la Construcción.”, ResearchGate, 30- Mar- 2020].
2011. [En lı́nea]. [Accessed: 29- Mar- 2020].