See discussions, stats, and author profiles for this publication at: https://www.researchgate.

net/publication/41905091

Ecuaciones diferenciales ordinarias y sus aplicaciones a la Ingeniería Civil

Book · July 2012

Source: OAI

CITATIONS READS

0 44,771

1 author:

Manuel Calixto
University of Granada
124 PUBLICATIONS 944 CITATIONS

SEE PROFILE

All content following this page was uploaded by Manuel Calixto on 23 December 2013.

The user has requested enhancement of the downloaded file.

 1

Ecuaciones Diferenciales Ordinarias

de Orden Superior

y
Aplicaciones a la Ingenier
a Civil

Manuel Calixto Molina

Departamento de Matem
atica Aplicada y Estad
stica,

Universidad Polit
ecnica de Cartagena,

Paseo Alfonso XIII, 52

30203 Cartagena

24 de Abril de 2002
2

Indice General

1 Teora sobre EDOs lineales de Orden Superior 5

1.1 Generalidades y estructura del espacio de soluciones . . . . . . . . . . . . . 5
1.1.1 Problemas de: valores iniciales (V.I.) y valores en la frontera (V.F.) 6
1.1.2 Ecuaciones homogeneas. Sistema fundamental de soluciones . . . . 7
1.1.3 Ecuaciones no homogeneas. Variacion de las constantes . . . . . . . 11
1.2 EDOs lineales con coe cientes constantes . . . . . . . . . . . . . . . . . . . 13
1.2.1 EDOs lineales homogeneas con coe cientes constantes . . . . . . . . 13
1.2.2 EDOs lineales no homogeneas con coe cientes constantes . . . . . . 14
1.3 Sistemas de EDOs lineales con coe cientes constantes . . . . . . . . . . . . 16
2 Aplicaciones a la Ingeniera Civil 19
2.1 Problemas de V.I.: vibraciones mecanicas y electricas . . . . . . . . . . . . 19
2.1.1 Oscilador armonico simple . . . . . . . . . . . . . . . . . . . . . . . 22
2.1.2 Oscilador armonico amortiguado . . . . . . . . . . . . . . . . . . . . 23
2.1.3 Oscilador armonico forzado: pulsaciones y resonancia pura . . . . . 25
2.1.4 Oscilador armonico amortiguado y forzado: factor de ampli cacion 28
2.1.5 Analogas electricas. Leyes de Kirchho . . . . . . . . . . . . . . . 30
2.2 Problemas de V.F.: exion y pandeo en vigas . . . . . . . . . . . . . . . . 30
2.2.1 Flexion en vigas: ecuacion diferencial de la curva elastica y echa
de exion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31
2.2.2 Pandeo en vigas: carga de Euler y modos de desviacion . . . . . . . 38
2.2.3 Derrumbe del puente colgante de Tacoma Narrows . . . . . . . . . 40

3
4 INDICE GENERAL
Cap
tulo 1

Teor
a sobre EDOs lineales de Orden
Superior

1.1 Generalidades y estructura del espacio de solu-

ciones
Una EDO lineal de n-esimo orden puede escribirse simbolicamente como:
an (x)y (n) + an 1 (x)y (n 1) + : : : a1 (x)y 0 + a0 (x)y = b(x) (1.1)
Esta ecuacion puede escribirse tambien como un sistema de n EDOs de orden uno de la
forma
~z0 = A(x)~z + ~b(x); (1.2)
0
donde ~z = (y; y ; : : : ; y (n 1)
),
0 1
0 1 0


0 0
B ... C
B 0 0 1 0 0 C
B ... C
B 0 0 0 0 0 C
B C
A(x) = B .. ... ... ... ... .. C
B . . C
B C
B ... C
@ 0 0 0 0 1 A
a0 (x)
an (x)
a1 (x)
an (x)
a2 (x)
an (x)


an 2 (x)
an (x)
an 1 (x)
an (x)
y
~b(x) = (0; 0;


; 0; b(x)=an (x))T ;
donde T indica trasposicion.
Esta traduccion de una EDO de orden n a un sistema equivalente de n EDOs de
orden 1 nos resultara util aqu, mas que desde un punto de vista practico, desde una
perspectiva teorica. Haremos uso extensivo de esta correspondencia entre ambos espacios
de soluciones.
Consideraremos dos tipos de problemas, dependiendo del tipo de restricciones sobre
la solucion general de (1.1)

5
6 CAPITULO 1. TEORIA SOBRE EDOS LINEALES DE ORDEN SUPERIOR

1.1.1 Problemas de: valores iniciales (V.I.) y valores en la fron-

tera (V.F.)
Problema de valores iniciales
Un problema de valores iniciales para una EDO lineal de orden n consiste en
Resolver : an (x)y (n) + an 1 (x)y (n 1) + : : : a1 (x)y 0 + a0 (x)y = b(x)
Sujeta a : y (x0 ) = y0 ; y 0(x0 ) = y 0 ; : : : ; y (n 1) (x0 ) = y (n 1) :
0 0 (1.3)
o en forma de sistema
~z0 = A(x)~z + ~b(x); ~z(x0 ) = ~z0 :
Para el problema de valores iniciales (1.3) tenemos un teorema de existencia y unicidad
de soluciones que nos dice lo siguiente:
Teorema 1.1.1. Sean an (x); an 1 (x); : : : ; a1 (x); a0 (x) y b(x) continuas en un intervalo
I  R , y sea an (x) 6= 0 8x 2 I . Si x = x0 es cualquier punto en el intervalo, existe una
unica solucion y (x) en dicho intervalo del problema de valores iniciales (1.3).
Observacion 1.1.2. Notese que si an (x) = 0 para alguna x en el intervalo I , la solucion
del problema lineal de valores iniciales (1.3) quiza no sea unica o incluso no exista; Por
ejemplo, es inmediato comprobar que la funcion y (x) = cx2 + x + 3 es una solucion del
problema de valores iniciales
x2 y 00 2xy 0 + 2y = 6; y (0) = 3; y 0(0) = 1
para x 2 ( 1; 1) y cualquier valor del parametro c. En otras palabras, no hay solucion
unica para este problema. Esto se debe a que, aunque los coe cientes an (x) son funciones
continuas en (1; 1), el coe ciente a2 (x) = x2 es cero en x = 0, precisamente donde se
han impuesto las condiciones iniciales.

Problema de valores en la frontera

Otro tipo de problema es resolver una EDO lineal en el que la variable dependiente y ,
o sus derivadas, esten especi cadas en puntos distintos. Un problema como
Resolver : a (x)y 00 + a (x)y 0 + a (x)y = b(x)
2 1 0
Sujeta a : y (x0 ) = y0 ; y (x1 ) = y1 ; (1.4)
se denomina problema de valores en la frontera. Las restricciones y (x0 ) = y0 ; y (x1 ) = y1 se
denominan condiciones en la frontera. Una solucion del problema anterior es una funcion
que satisface la ecuacion diferencial en algun intervalo I que contiene a x0 y x1 , cuya
gra ca pasa por los puntos (x0 ; y0 ) y (x1 ; y1 ). Estas condiciones en la frontera son casos
particulares de las condiciones generales:
0
0 y (x0 ) + 0 y (x0 ) = 0 ;
0
1 y (x1 ) + 1 y (x1 ) = 1 :
1.1. GENERALIDADES Y ESTRUCTURA DEL ESPACIO DE SOLUCIONES 7

Los ejemplos que siguen demuestran que aun cuando se satisfagan las condiciones del
teorema 1.1.1, un problema de valor en la frontera puede tener i) varias soluciones; ii)
solucion unica, o iii) ninguna solucion. En efecto la solucion general de la EDO lineal de
orden 2: y 00 + 16y = 0 es y = c1 cos(4x) + c2 sen(4x). Ahora:
i) Si imponemos y (0) = 0; y (=2) = 0 ) c1 = 0, con c2 arbitraria. Es decir, tenemos
una familia uniparametrica de soluciones y = c2 sen(4x) que pasan por los puntos
(0; 0) y (=2; 0) y que satisfacen la ecuacion y 00 + 16y = 0.
ii) Si imponemos y (0) = 0; y (=8) = 0 ) c1 = c2 = 0. Es decir, y = 0 es la unica
solucion de y 00 + 16y = 0 que pasa por estos puntos.
ii) Si imponemos y (0) = 0; y (=2) = 1 ) c1 = 0; 1; c2 =arbitraria. Es decir, el
problema no tiene solucion.

1.1.2 Ecuaciones homogeneas. Sistema fundamental de solucio-

nes
Una EDO lineal de orden n como (1.3) se dice homogenea si el lado derecho de la
igualdad es identicamente nulo: b(x) = 0. Veremos que la resolucion de una ecuacion no
homogenea como (1.3) pasa por la resolucion de la ecuacion homogenea asociada.

Operadores diferenciales
dy = Dy , donde el smbolo D se llama operador diferencial
A veces es util denotar dx
porque transforma una funcion diferenciable en otra funcionn. Las derivadas de orden
superior se pueden expresar como \potencias"de D, como dx d y = D n y . Las expresiones
n
polinomicas en D, como L = D2 + 3xD + 5 tambien son operadores diferenciales. En
general:
De nicion 1.1.3. Se de ne un operador diferencial de orden n como:
L = an (x)Dn + an 1 (x)Dn 1 +


+ a1 (x)D + a0 (x); (1.5)
donde ai (x); i = 0; 1; : : : ; n son funciones reales (o complejas) de x.
Denotemos por C (m) el conjunto de las funciones reales (o complejas) f : R ! R de
clase m. Obviamente C (m) es un espacio vectorial sobre el cuerpo R , y L es una aplicacion
L : C (m) ! C (l) ; m  n que hace corresponder a cada f 2 C (m) una funcion L(f ) de clase
C (l) . La aplicacion \multiplicacion por a0 (x)"hace corresponder a cada funcion f (x) la
funcion a0 (x)f (x). Veamos que la aplicacion L : C (m) ! C (l) ; m  n es una aplicacion
lineal entre espacios vectoriales.
Teorema 1.1.4. El operador L tiene la propiedad de linealidad
L( f (x) + g (x)) = L(f (x)) + L(g (x)); (1.6)
donde y son constantes y f; g 2 C (m) ; m  n. A causa de esta propiedad, el operador
diferencial de orden n (1.5) se denomina \ operador lineal".
8 CAPITULO 1. TEORIA SOBRE EDOS LINEALES DE ORDEN SUPERIOR

Demostracion: la demostracion radica en dos propiedades basicas de la diferenciacion:

1) D( f (x)) = Df (x), donde es una constante, y 2) D(f (x)+ g (x)) = Df (x)+ Dg (x).
Estas propiedades son claramente extensibles a la derivada n-esima Dn, incluyendo D0 
1, y a sus combinaciones lineales con coe cientes ak (x). 
Observacion 1.1.5. Notese que, al contrario que D, el operador L no tiene porque cumplir
la regla de Leibnitz D(f (x)g (x)) = D(f (x))g (x)+f (x)D(g (x)). Esto se debe a la presencia
del termino a0 (x) en (1.5).

Notacion operatorial de EDOs lineales

La EDO lineal (1.3) puede expresarse en forma compacta en terminos del operador
diferencial lineal (1.5) como:
L(y ) = 0;
para el caso homogeneo b(x) = 0, o bien en forma de un sistema de n EDOs de orden uno
como:
D~z = A(x)~z ) (DIn A(x))~z = ~0;
donde In denota la matriz identidad n  n.
As, resolver la ecuacion L(y ) = 0 signi ca encontrar el nucleo del operador L entre
las funciones reales n veces derivables C (n) . O, equivalentemente, encontrar el nucleo del
operador lineal DIn A(x) entre las funciones vectoriales derivables ~z : I ! R n .

Principio de superposicion
La propiedad de linealidad de L permite, dadas dos o mas soluciones de L(y ) = 0,
encontrar otra. Esta idea viene expresada de manera precisa en el siguiente teorema
Teorema 1.1.6. Sean y1 ; y2 ; : : : ; yk soluciones de la ecuacion diferencial homogenea de
orden n L(y ) = 0, donde x esta en un intervalo I . La combinacion lineal
y (x) = y (x) + 2 y2 (x) +


+
1 1 k yk (x);

donde i ; i = 1; 2; : : : ; k son constantes arbitrarias, tambien es una solucion cuando x

esta en el intervalo I .
Demostracion: la clave esta en la linealidad de L:
L(y ) = L( 1 y1 + 2 y2 +


+ k yk ) = 1 L(y1 ) + 2 L(y2 ) +


+ k L(yk ) = 0;
donde, en la ultima igualdad hemos usado que L(yi ) = 0; i = 1; : : : ; k, es decir, que
y1 ; y2 ; : : : ; yk son soluciones de la EDO lineal homogenea de orden n: L(y ) = 0. 
Corolario 1.1.7. Del teorema anterior se deduce inmediatamente que las soluciones de
L(y ) = 0 forman un espacio vectorial sobre R .
Veamos como encontrar la dimension de este espacio vectorial real.
1.1. GENERALIDADES Y ESTRUCTURA DEL ESPACIO DE SOLUCIONES 9

Dependencia e independencia lineal

Citaremos un par de conceptos basicos para estudiar EDOs lineales.
De nicion 1.1.8. Se dice que un conjunto de funciones, f1 (x); f2 (x); : : : ; fn(x) es lineal-
mente dependiente en un intervalo I si existen constantes , 1 ; 2 ; : : : ; n 2 R no todas
nulas, tales que
1 f1 (x) + 2 f2 (x) +


+ n fn (x) = 0; 8x 2 I:

Si el conjunto de funciones no es linealmente dependiente, se dice que es linealmente

independiente.
Para un conjunto de dos funciones, la dependencia lineal 1 f1 (x)+ 2 f2 (x) = 0 signi -
ca que una funcion es multiplo de la otra: f1 (x) = 21 f2 (x) (suponiendo que 1 6= 0). Por
ejemplo, las funciones f1 (x) = sen(2x) y f2 (x) = sen(x) cos(x) son linealmente dependien-
tes en I = ( 1; 1) porque f1 (x) = sen(2x) = 2 sen(x) cos(x) = 2f2 (x). Por otro lado,
las funciones f1 (x) = x y f2 (x) = jxj son linealmente independientes en I = ( 1; 1)
(aunque no en los intervalos ( 1; 0) y (0; 1), donde son dependientes).

Generalidades sobre las soluciones de EDOs lineales homogeneas

Antes de entrar en formulas explcitas para las soluciones, las cuales solo pueden pro-
porcionarse salvo casos excepcionales como el caso de coe cientes constantes, estudiemos
con detenimiento la estructura del espacio de soluciones.
Teorema 1.1.9. Las soluciones de L(y ) = 0, o equivalentemente de (DIn A(x))~z = ~0,
forman un espacio vectorial de dimension n sobre R .
Demostracion: hemos visto que la suma de soluciones y el producto de soluciones de
L(y ) = 0 por escalares son tambien solucion. Para conocer la dimension utilizaremos la
equivalencia L(y ) = 0 , (DIn A(x))~z = ~0. Fijemos una base ~v1 ; ~v2 ; : : :~vn de R n y
x0 2 I , y consideremos las soluciones ~z1 ; ~z2 ; : : : ~zn de (DIn A(x))~z = ~0 con condiciones
iniciales ~zi (x0 ) = ~vi ; i = 1; 2; : : : ; n.
Basta probar que ~z1 ; ~z2 ; : : : ~zn forman una base del espacio de soluciones (linealmente
independiente y generador). Para ello procedamos por reduccion al absurdo. Supongamos
que existen constantes 1 ; 2 ; : : : ; n 2 R no todas nulas, tales que
~z (x) + 2~z2 (x) +


+
1 1 n~
zn (x) = ~0; 8x 2 I:
Sustituyendo x por x0 ,
1 1 ~v + 2~v2 +


+ n~vn = ~0;
pero esto contradice que ~v1 ; ~v2 ; : : :~vn sea una base de R n , luego ~z1 ; ~z2 ; : : : ~zn son linealmente
independientes. Para ver que son un sistema generador, supongamos que ~z es una solucion
de (DIn A(x))~z = ~0 y pongamos ~z(x0 ) = ~v. Sean 1 ; 2 ; : : : ; n 2 R las coordenadas de
~v en la base ~v1 ; ~v2 ; : : :~vn de R n . Entonces
~z (x0 ) + 2~z2 (x0 ) +


+
1 1 n~
zn (x0 ) = ~v ;
10 CAPITULO 1. TEORIA SOBRE EDOS LINEALES DE ORDEN SUPERIOR

y por unicidad de soluciones del problema de Cauchy (Teorema 1.1.1)

~z (x) + 2~z2 (x) +


+
1 1 n~
zn (x) = ~z(x):
Por tanto ~z1 ; ~z2 ; : : : ~zn forman un sistema generador. 
De nicion 1.1.10. Llamaremos a cualquiera de las bases S = f~z1 ; ~z2 ; : : : ~zn g del espacio
de soluciones de (DIn A(x))~z = ~0 un sistema fundamental de soluciones. De las matrices
0 1 0 1
z1;0 z2;0


zn;0 y1 y2


yn
B z
1;1 z2;1


zn;1 CC B y1
B 0 y20


yn0 C
Z=B .. C 2 Mnn (R )
C
B .. .. .. =
.. C B .. .. ..
@ . . . . A @ . . . . A
z1;n 1 z2;n 1


zn;n 1 y(n 1)
1 y (n 1)
2


yn (n 1)

cuyas columnas ~z1 ; ~z2 ; : : : ~zn son un sistema fundamental S diremos que son matrices
fundamentales M [S ] = Z de (DIn A(x))~z = ~0 $ L(y ) = 0.
Entre las familias F = fy1 (x); y2 (x); : : : ; yn (x)g de soluciones de L(y ) = 0 distinguire-
mos las que constituyen un sistema fundamental de las que no a traves del Wronskiano,
que se de ne como el determinante
y1 (x) y2 (x)


yn(x)
y10 (x) y20 (x)


yn0 (x)
W [F ](x) = jM [F ](x)j = .. .. .. .. :
. . . .
y1(n 1) (x) y2(n 1) (x)


yn(n 1) (x)
Teorema 1.1.11. Sean y1 (x); y2 (x); : : : ; yn (x) soluciones de L(y ) = 0. Entonces las
siguientes a rmaciones son equivalentes:
(i) La familia F = f~z1 ; ~z2 ; : : : ~zn g, formada por los vectores con componentes ~zk =
(yk ; yk0 ; : : : ; yk(n 1) ), es un sistema fundamental;
(ii) existe x0 2 I tal que W [F ](x0) 6= 0;
(iii) para todo x 2 I; W [F ](x) =
6 0
Demostracion: en virtud del isomor smo lineal entre el espacio de soluciones de L(y ) =
0 y (DIn A(x))~z = ~0, podemos demostrar que (i))(ii). En efecto, supongamos
que M [F ](x) es una matriz fundamental pero que, sin embargo, existe x1 2 I tal que
det M [F ](x1 ) = W [F ](x1 ) 6= 0. Existiran entonces 1 ; 2 ; : : : ; n 2 R no todas nulas,
tales que 1~z1 (x1 ) + 2~z2 (x1 ) +


+ n~zn (x1 ) = ~0. De namos ~z(x) = 1~z1 (x) + 2 ~z2 (x) +



+ n~zn(x). Obviamente ~z(x) es solucion de (DIn A(x))~z = ~0 con la condicion inicial
~z(x1 ) = ~0, al igual que ~z(x) = ~0, luego, por unicidad, ~z(x) = ~0. Se contradice as que la
familia F sea libre.
Trivialmente (iii))(ii).
Para probar que (ii))(i), basta notar que ~vi = ~zi (x0 ); i = 1; 2; : : : ; n es una base de
n
R y razonar como en el Teorema 1.1.9 para concluir que F es un sistema fundamental.

1.1. GENERALIDADES Y ESTRUCTURA DEL ESPACIO DE SOLUCIONES 11

1.1.3 Ecuaciones no homogeneas. Variacion de las constantes

Las soluciones de la ecuacion lineal homogenea L(y ) = 0 $ (DIn A(x))~z = ~0 y de la
no homogenea L(y ) = b(x) $ (DIn A(x))~z = ~b(x) guardan una relacion muy estrecha,
como pone de mani esto el siguiente resultado:
Teorema 1.1.12. Sea S = fy1 (x); y2 (x); : : : ; yn(x)g un sistema fundamental de L(y ) = 0
e yp(x) una solucion cualquiera de L(y ) = b(x). Entonces la solucion general de L(y ) =
b(x) viene dada por
s:g:n:h:
z }| {
y (x) = yp(x) + c|1y1 (x) + c2 y2 (x{z) +


+ cn yn (x});
| {z }
s:p:n:h: s:g:h:

donde c1 ; c2 ; : : : ; cn son constantes arbitrarias.

Demostracion: sea y (x) la solucion general de la ecuacion no homogenea (s.g.n.h.)
L(y ) = b(x) e yp(x) una solucion particular de L(y ) = b(x) (s.p.n.h.). Si de nimos
u(x) = y (x) yp(x), por la linealidad de L se debe cumplir
L(u) = L(y ) L(yp ) = b(x) b(x) = 0:
Esto demuestra que u(x) es una solucion de la ecuacion homogenea L(y ) = 0; por
consiguiente, segun el Teorema 1.1.9, existen constantes c1 ; c2 ; : : : ; cn tales que u(x) =
c1 y1 (x) + c2 y2 (x) +


+ cn yn (x) = y (x) yp (x) . 

Principio de superposicion para ecuaciones no homogeneas

Respecto al calculo de soluciones particulares yp de L(y ) = b(x), a veces resulta util
el siguiente \principio de superposicion"para EDOs lineales no homogeneas (vease mas
adelante el metodo de los coe cientes indeterminados ).
Teorema 1.1.13. Sean yp1 ; yp2 ; : : : ; ypk soluciones particulares respectivas de
L(y ) = h1 (x); L(y ) = h2 (x); : : : ; L(y ) = hk (x)
en un cierto intervalo I , entonces yp = 1 yp1 + 2 yp2 +


+ k ypk , con i ; i = 1; : : : ; k
constantes reales, es solucion de

L(y ) = b(x) = 1 1 h (x) + 2 h2 (x) +


+ k hk (x):

Demostracion: en efecto, este resultado es consecuencia de la linealidad del operador

diferencial L

L( 1 yp1 + 2 yp2 +


+ k ypk ) = h (x) +

1 1 h (x) +


+
2 2 k hk (x) = b(x):
12 CAPITULO 1. TEORIA SOBRE EDOS LINEALES DE ORDEN SUPERIOR

Metodo de variacion de las constantes

Como en el caso de los sistemas lineales, la busqueda de la solucion general de la
ecuacion lineal no homogenea pasa por encontrar un sistema fundamental de soluciones
de la ecuacion homogenea y despues una solucion particular de la completa. Cuando ya
se conoce la s.g.h., el calculo de una solucion particular puede hacerse por el metodo de
variacion de las constantes. En este caso, suponemos que la solucion particular que busca-
mos podra escribirse, en terminos del sistema fundamental S = fy1 (x); y2 (x); : : : ; yn(x)g
del que ya disponemos, como

yp (x) = c1 (x)y1 (x) + c2 (x)y2 (x) +


+ cn (x)yn (x) (1.7)

para funciones adecuadas ci (x); i = 1; : : : ; n. Para comprenderlo mejor, hacemos uso del
isomor smo lineal entre el espacio de soluciones de L(y ) = 0 y (DIn A(x))~z = ~0. As,
la ecuacion (1.7) puede escribirse de manera equivalente como

~zp = Z (x)~c(x);
donde ~zp = (yp ; yp0 ; : : : ; yp(n 1) )T , Z es la matriz fundamental de nida en 1.1.10 y ~c(x) =
(c1 ; c2 ; : : : ; cn )T (x). Notese que, para una funcion ~zp de nida de esta manera ocurrira que

~zp0 (x) = Z 0 (x)~c(x) + Z (x)~c0 (x) = A(x)Z (x)~c(x) + Z (x)~c0 (x);

donde hemos utilizado que (DIn A(x))~zk = ~0; k = 1; : : : ; n para las columnas ~zk de
la matriz fundamental Z (x). Si queremos que ~zp cumpla la ecuacion no homogenea
~zp0 (x) = A(x)~zp (x) + ~b(x), se concluye que
Z Z !
W [S ](x)
Z (x)~c0 (x) = ~b(x) ) ~c(x) = Z (x)~b(x)dx =
1
dx;
W [S ](x)
donde, en la ultima igualdad, se ha utilizado el metodo de Cramer con
!
W [S ] = (W1 [S ]; W2 [S ]; : : : ; Wn[S ])
y Wj [S ] el determinante que resulta de sustituir la columna j en el Wronskiano W [S ] por
~b. As, la solucion general de (DIn A(x))~z = ~b(x) es
Z
~z(x) = Z (x)~ + Z (x) Z 1 (x)~b(x)dx;

o, si se imponen condiciones iniciales ~z(x0 ) = ~z0 :

Z x
~z(x) = Z (x)Z (x0 )~z0 + Z (x)
1
Z 1 (t)~b(t)dt: (1.8)
x0
1.2. EDOS LINEALES CON COEFICIENTES CONSTANTES 13

1.2 EDOs lineales con coe cientes constantes

Existen metodos de resolucion de ecuaciones lineales como (1.3) basados en desarrollos
en serie de potencias. No obstante, nosotros no trataremos estos metodos aqu y pasaremos
a abordar el caso mas sencillo en que los coe cientes aj (x) de (1.3) son constantes ak (x) =
ak . Sin perdida de generalidad podemos tomar an = 1, de manera que la ecuacion a
estudiar es
L(y ) = y (n) + an 1 y (n 1) +


+ a1 y 0 + a0 = b(x):
Comencemos por buscar un sistema fundamental de la EDO homogenea L(y ) = 0.

1.2.1 EDOs lineales homogeneas con coe cientes constantes

Para encontrar la solucion general de

L(y ) = y (n) + an 1 y (n 1)
+


+ a1 y 0 + a0 = 0; (1.9)

ensayaremos funciones del tipo y (x) = ex , que sustituidas en la ecuacion anterior L(y ) = 0
nos queda L(ex ) = ex L[] = 0, donde

L[] = n + an 1 n 1 +


+ a1  + a0 ; (1.10)

denota un polinomio de grado n en  al que denominaremos polinomio caracterstico de

(1.9). La ecuacion L[] = 0 se denomina ecuacion caracterstica.

Teorema 1.2.1. La solucion general de (1.9) es una combinacion lineal de las funciones
xk e x cos(!x); xk e x sen(!x);
donde = + i! recorre el conjunto de las races de (1.10) con !  0, y 0  k  m( ),
siendo m( ) la multiplicidad de .

Demostracion: basta demostrar que si es una raz de L[] de multiplicidad m y 0 

k  m, entonces xk e x es solucion de L[D](y ) = 0 (con L[D]  L).1 A este n, apreciese
que si L[] factoriza como L[] = L1 []L2 [] entonces si L2 [](y ) = 0 ) L[](y ) = 0.
Como nuestro objetivo es demostrar que si ( )m divide a L[] entonces L[D](xk e x ) = 0
para cada 0  k < m, bastara con probar que

(D )k+1(xk e x ) = 0 para cada k  0:

1 N
otese que si, en particular, = + i! tiene parte imaginaria ! = 0, entonces se tendra
6

0 = L[D](xk e x ) = L[D](xk e x cos(!x)+ ixk eax sen(!x)) = L[D](xk e x cos(!x))+ iL[D](xk eax sen(!x));

y por tanto L[D](xk e x cos(!x)) = 0 = L[D](xk e x sen(!x)).

14 CAPITULO 1. TEORIA SOBRE EDOS LINEALES DE ORDEN SUPERIOR

Esto es simple si razonamos por induccion. La a rmacion es obviamente cierta cuando

k = 0. Supuesta cierta para k 1 ! (D )k (xk 1 e x ) = 0, se tiene que:


(D )k+1 xk e x = (D )k (D )xk e x
= (D )k kxk 1 e x + xk e x xk e x



= k(D )k xk 1 e x = 0;
de modo que tambien sera cierta para k. Es inmediato comprobar tambien que dichas
soluciones son linealmente independientes. 

1.2.2 EDOs lineales no homogeneas con coe cientes constantes

Ya hemos visto un metodo general para el calculo de soluciones particulares de la
ecuacion no homogenea L(y ) = b(x) conocido un sistema fundamental de la ecuacion
homogenea L(y ) = 0, denominado metodo de variacion de las constantes [vease (1.8)].
Existe otro metodo mas directo para el caso en que el termino inhomogeneo b(x) sea
producto de exponenciales, polinomios en x, senos y cosenos.

Metodo de los coe cientes indeterminados

Teorema 1.2.2. Supongamos que
b(x) = e x (p(x) cos(!x) + q (x) sen(!x));
donde p(x) y q (x) son polinomios de grado a lo sumo k  0. Sea  = + i! . Entonces
se tienen las siguientes posibilidades:
(i) Si  NO es raz del polinomio caracterstico (1.10) entonces L(y ) = b(x) tiene una
solucion particular de la forma yp (x) = e x (r(x) cos(!x) + s(x) sen(!x)), con r(x)
y s(x) polinomios de grado a lo sumo k.
(ii) Si  es raz del polinomio caracterstico (1.10) con multiplicidad m, entonces L(y ) =
b(x) tiene una solucion particular de la forma yp (x) = e x xm (r(x) cos(!x)+s(x) sen(!x)),
con r(x) y s(x) polinomios de grado a lo sumo k.
Demostracion: utilizaremos la version compleja
L(y ) = y(n) + an 1 y(n 1)
+


+ a1 y0 + a0 y = ex v (x) = b(x); (1.11)
de la ecuacion real L(y ) = b(x), donde v (x) = p(x) iq (x) es un polinomio de grado menor
o igual que k con coe cientes complejos. Demostraremos que existe una solucion compleja
yp : R ! C de (1.11) de la forma yp(x) = xm ex w(x), donde m es la multiplicidad de 
como solucion de
n + an 1 n 1 +


+ a1  + a0 = 0
(m = 0 si  no es raz de la ecuacion) y w(x) = r(x) is(x) es un polinomio de gra-
do menor o igual que k con coe cientes complejos. La parte real yp (x)  <(yp(x)) =
1.2. EDOS LINEALES CON COEFICIENTES CONSTANTES 15

e x xm (r(x) cos(!x)+s(x) sen(!x)) sera pues solucion de L(y ) = b(x) con b(x) = <(b(x)) =
<(exv(x)) = e x(p(x) cos(!x) + q(x) sen(!x)).
Notese que la ecuacion (1.11) se puede reescribir en terminos de la descomposicion en
producto de monomios (D m
j ) j , con j una raz de (1.10) de multiplicidad mj , de la
siguiente forma:
(D 1 )m1 (D 2 )m2 : : : (D r)
mr (D )m y = ex v (x); (1.12)
donde admitimos la posibilidad de que  sea o no sea una raz, poniendo en este ultimo
caso: m = 0 ) (D )0 = 1. Procedamos por partes:
a) Comencemos por el caso m = 0, o sea, cuando  no es raz de (1.10). Bastara
probar que existe un polinomio w(x) de grado  k tal que yp(x) = ex w(x) es solucion de
(D y (x)
j ) = ex v~(x) (1.13)
para algun polinomio v~(x) de grado  k. En efecto, sustituyendo yp(x) en (1.13)
ex v~(x) = yp0 (x) j yp (x) = ex (w0 (x) + w(x) ) v~(x) = ( j )w(x) + w0(x):
j w (x))

Poniendo v~(x) = v~0 + v~1 x +


+ v~k xk y w(x) = w0 + w1 x +


+ wk xk e igualando

coe cientes de igual grado:
v~0 = ( j )w0 + w1 ;
v~1 = ( j )w1 + 2w2 ;




v~k 1 = ( j )wk 1 + kwk ;
v~k = ( j )wk ;

que obviamente tiene solucion por ser  6= j ( no es raz). Procediendo de forma iterada
llegamos a que existe w(x) tal que yp (x) = ex w(x) es solucion de (1.12) con m = 0.
b) Analizamos a continuacion el caso en que  sea una raz con multiplicidad m = n )
mj = 0. Se trata de comprobar que yp (x) = xn ex w(x) es una solucion de (D )ny =
ex v (x) con w(x) un polinomio de grado  k. De hecho demostraremos bastante mas:
que para cada n  1 y l  0 la ecuacion
(D )n y = xl ex v (x) (1.14)
tiene una solucion de la forma yp(x) = xl+n ex w(x), con w(x) un polinomio de grado  k.
Razonemos por induccion sobre n y supongamos n = 1. La funcion yp(x) = xl+1 ex w(x)
sera solucion de (D )y = xl ex v (x) si y solo si
xl ex v (x) = xl+1 ex w(x) + [(l + 1)xl w(x) + xl+1 w0(x)]ex xl+1 ex w(x);
o sea,
v (x) = (l + 1)w(x) + xw0 (x);
16 CAPITULO 1. TEORIA SOBRE EDOS LINEALES DE ORDEN SUPERIOR

o igualando coe cientes de igual grado

v0 = (l + 1)w0 ;
v1 = (l + 2)w1 ;




vk 1 = (l + k)wk 1 ;
vk = (l + k + 1)wk ;
y despejando obtenemos wj = vj =(j + l + 1).
Para completar nuestro argumento por induccion supongamos que y~p(x) = xl+n ex w~ (x)
es solucion de ecuacion (D )n y = xl ex v (x) para un cierto n y demostremos entonces
que yp(x) = xl+n+1 ex w(x) es solucion de la ecuacion (D )n+1 y = xl ex v (x) utilizando
que, como ya hemos demostrado, (D )yp = xl+n ex w~(x) para algun polinomio w~(x) de
grado  k. En efecto:
(D )n+1 yp = (D )n[(D )yp]
= (D )n[x| l+n e{z ~(x)] = xl ex w(x)
x w
}
y~p

c) Acabamos de demostrar en el apartado b) que yp(x) = xm ex w(x) es solucion de la

ecuacion
(D )m y = ex w~(x)  y~p (x);
con w~(x) un polinomio de grado  k, y en el apartado a) que, a su vez, y~p(x) es solucion
de
(D m1 m2 : : : (D mr x
1 ) (D 2) r ) y = e v (x):
Por consiguiente, yp(x) = xm ex w(x) es solucion de (1.12).

1.3 Sistemas de EDOs lineales con coe cientes cons-

tantes
Consideremos el sistema
9
L11 [D]y1 (x) + L12 [D]y2 (x) +


+ L1n [D]yn (x) = b1 (x); >
>
>
L21 [D]y1 (x) + L22 [D]y2 (x) +


+ L2n [D]yn (x) = b2 (x); =
.. ! L(~y(x)) = ~b(x); (1.15)
. >
>
>
Ln1 [D]y1 (x) + Ln2 [D]y2 (x) +


+ Lnn [D]yn(x) = bn (x); ;

donde L denota ahora una matriz n  n cuyos coe cientes Lij [D] son operadores lineales
como los de nidos en (1.5). Para resolver este sistema diferencial lo inmediato es aplicar
el algoritmo de triangulacion de Gauss tratando el sistema formalmente como si las fun-
ciones incognita ~y = (y1; y2 ; : : : ; yn) fuesen numeros y sus coe cientes Lij [D] polinomios
1.3. SISTEMAS DE EDOS LINEALES CON COEFICIENTES CONSTANTES 17

en el \parametro"D. Sabemos que en este proceso de \triangulacion"se reemplaza una

la (o columna) por la que se obtiene al sumar a dicha la una \combinacion lineal"de
las restantes, en tendiendo aqu por P \combinacion lineal"la \multiplicacion"de las por
operadores diferenciales polinomiales j cj Dj arbitrarios (no necesariamente constantes).
En el siguiente captulo estudiaremos las ecuaciones de movimiento de sistemas de osci-
ladores armonicos acoplados, que pueden verse como casos particulares como los sistemas
de EDOs lineales de orden 2 de la forma
L(~y ) = MD2 ~y + CD~y + K~y = F~ (x);
donde M se denomina matriz de masas, K es la matriz de resortes o de constantes de
acoplamiento, C es la matriz de constantes de amortiguamiento, F~ es un vector de fuerzas
externas, e ~y(x) denota el vector de posicion ~y del sistema en funcion del tiempo x. En
vez de aplicar el algoritmo de triangulacion de Gauss, aqu es mas conveniente (desde el
punto de vista interpretativo) ensayar una solucion del tipo ~y(x) = e!x ~ para el sistema
homogeneo
(M2 + C + K )~ = ~0;
el cual tendra solucion si det(M2 + C + K ) = 0, que constituye una ecuacion carac-
terstica para el sistema. Un caso particularmente interesante es cuando C = 0 o C / K .
En tal caso, el sistema se convierte es una ecuacion de autovalores
M2 ~ = K~ ;
con autovalores j = i!j (\frecuencias naturales de vibracion") complejos puros (debido
a que M y K son en general matrices simetricas de nidas positivas) y autovectores ~j
(\modos normales de vibracion"). Por razones de tiempo, no trataremos aplicaciones
fsicas concretas de esta parte en el siguiente captulo.
18 CAPITULO 1. TEORIA SOBRE EDOS LINEALES DE ORDEN SUPERIOR
Cap
tulo 2

Aplicaciones a la Ingenier
a Civil

2.1 Problemas de V.I.: vibraciones mecanicas y electricas

Estudiaremos el movimiento unidimensional de una partcula de masa m sometida
a fuerzas de diferente naturaleza: elasticas o restauradoras Fe (t; x), viscosas o resistivas
Fv (t; x; x_ ) y otras fuerzas externas Fext (t) dependientes del tiempo. La ecuacion diferencial
que describe la evolucion de la posicion x(t) en funcion del tiempo viene dada por la
segunda ley de Newton:
mx = Fe (t; x) + Fv (t; x; x_ ) + Fext (t); (2.1)
que es una EDO de orden 2 cuyas condiciones iniciales
x(t0 ) = x0 ; v (t0 ) = x_ (t0 ) = v0 ;
consisten en especi car la posicion x(t) y la velocidad v (t) = dxdt(t)  x_ (t) en un instante
inicial t0 .
Pasemos a discutir distintas expresiones de dichas fuerzas y sus aproximaciones li-
neales, las cuales tomaremos como punto de partida para un tratamiento analtico de la
ecuacion (2.1)
1. Fuerzas recuperadoras o elasticas Fe . Hemos permitido, en general, fuerzas restau-
radoras Fe (t; x) dependientes de la posicion y del tiempo. Supongamos que x denota
un desplazamiento alrededor de una posicion de equilibrio y que Fe (t; x) admite un
desarrollo en serie alrededor de dicha posicion de equilibrio de la forma:
Fe (t; x) = k1 (t)x k3 (t)x3 + : : : ;
donde solo intervienen potencias impares del desplazamiento x debido que la fuerza
restauradora debe ser impar Fe (t; x) = Fe (t; x).
(a) Resorte lineal. Cuando los desplazamientos son peque~nos, podemos despreciar
ordenes superiores en el desarrollo de Fe y considerar
Fe (t; x) = k1 (t)x:

19
20 CAPITULO 2. APLICACIONES A LA INGENIERIA CIVIL

i. Ley de Hooke. El caso mas sencillo e ideal es aquel en el cual las carac-
tersticas fsicas del resorte no cambian con el tiempo, es decir, cuando
k1 (t) = k > 0 es independiente del tiempo, con lo cual la fuerza elastica
queda:
Fe (t; x) = kx:
ii. Resorte desgastable. Sin embargo, en el mundo real es logico esperar que el
resorte de debilite (o \pierda bro") conforme pasa el tiempo; por ejemplo,
de la forma k1 (t) = ke t ; k; > 0. En este caso, la EDO del sistema
\masa-resorte"
mx + ke t x = 0
se transforma en una
d2 x dx 2
Ec: de Bessel ( = 0) : s2 + s + s x = 0:
ds2 ds
mediante el cambio de variable
r
2 k t=2 :
s= e
m
La solucion de esta ecuacion se obtiene por desarrollo en serie de potencias
(nosotros no trataremos este metodo aqu) y puede escribirse en terminos
de las funciones de Bessel de primera y segunda clase de la forma:
r ! r !
2 k t=2 2 k t=2
x(t) = c1 J0 e + c2 Y0 e :
m m
iii. Resorte que se endurece. Cuando un resorte se somete a un ambiente en
que la temperatura decrece rapidamente, la constante elastica crece con el
tiempo, por ejemplo, de forma lineal k1 (t) = kt. En este caso, la EDO del
sistema \masa-resorte"viene descrito por la
Ec: de Airy : mx + ktx = 0:
que se resuelve por series de potencias. Dicha ecuacion tambien aparece
al estudiar la difraccion de la luz, la difraccion de las ondas de radio en
torno a la super cie de la Tierra, en aerodinamica y en el pandeo de una
columna vertical uniforme que se exiona bajo su propio peso.
(b) Resorte no lineal. Para desplazamientos no tan peque~nos, el muelle abandona
su comportamiento lineal y la fuerza elastica adopta (en el caso independiente
del tiempo) la forma Fe (x) = kx k3 x3 , donde hemos despreciado terminos
de orden superior. DuÆng estudio las oscilaciones forzadas de una masa pun-
tual sometida a una fuerza recuperadora no lineal modelada por la ecuacion
diferencial:
Ec: de DuÆng : mx + kx + k3 x3 = F0 cos !t:
Caben dos casos interesantes aqu:
 
2.1. PROBLEMAS DE V.I.: VIBRACIONES MECANICAS 
Y ELECTRICAS 21

i. Resorte duro. k3 > 0, para el cual la fuerza recuperadora es mayor (en

valor absoluto) que la del resorte lineal.
ii. Resorte suave. k3 < 0, para el cual la fuerza recuperadora es menor (en
valor absoluto) que la del resorte lineal.

F
x
Lineal

Duro
x
Blando

Figura 2.1: Fuerza Fe (x) = kx k3 x3 para un resorte lineal (k3 = 0), uno duro (k3 > 0) y
uno blando (k3 < 0)

2. Fuerzas resistivas o viscosas Fv . Supongamos que Fv (t; x; x_ ) admite un desarrollo

en serie de la forma:
Fv (t; x; x_ ) = c1 (t; x)x_ c2 (t; x)jx_ jx_ c3 (t; x)jx_ j2 x_ + : : : ;
donde los valores absolutos garantizan que la fuerza viscosa se oponga siempre al
movimiento.
(a) Fuerza viscosa lineal. Para velocidades peque~nas, se pueden despreciar terminos
de orden superior y tomar
Fv (t; x; x_ ) = c1 (t; x)x:
_
Caben varias posibilidades aqu:
i. Coe ciente de viscosidad constante. c1 (t; x) = c=constante> 0.
ii. Coe ciente de viscosidad variable. Dentro de las posibilidades destacamos
el caso c1 (t; x) = c(x2 a2 ) que da lugar a la
Ec: de van der Pol : mx + c(x2 a2 )x_ + kx = 0:
El hecho de que la fuerza viscosa sea Fv (jxj > jaj) < 0 y Fv (jxj < jaj) > 0
signi ca que la partcula se acelera para jxj < jaj y se retarda para jxj > jaj
22 CAPITULO 2. APLICACIONES A LA INGENIERIA CIVIL

(o, de otra forma, la partcula absorbe energa cuando jxj < jaj y disipa
cuando jxj > jaj) tendiendo, por tanto, permanecer en oscilacion estacio-
naria. Este es el resultado del Teorema de Lienard (vease [Simmons], pag.
518), el cual asegura para ecuaciones como la de van der Pol la existencia
de una unica trayectoria cerrada que rodea al origen en el plano de fases,
a la que tienden en forma de espirales todas las demas trayectorias cuando
t ! 1.
(b) Fuerza viscosa no lineal. Para velocidades no tan peque~nas, y dependiendo
del tipo de medio viscoso, a veces es necesario a~nadir nuevos terminos en el
desarrollo de Fv , por ejemplo:
Fv (x_ ) = c1 x_ c2 jx_ jx:
_
En la siguiente seccion solo consideraremos el caso lineal independiente del tiempo, con
lo cual, la ecuacion a estudiar es:

mx + cx_ + kx = Fext (t) ! x + 2 x_ + !02 x = Fext (t)=m;

q
donde hemos puesto 2 = c=m y !0 = mk por comodidad de calculo (vease mas ade-
lante). Como fuerzas externas Fext (t) utilizaremos, por simplicidad, funciones de tipo
periodico Fext (t) = Fext (t + T ).

2.1.1 Oscilador armonico simple

Comencemos considerando el caso en que no existe fuerza viscosa ( = 0) ni fuerzas
externas Fext (t) = 0. As, la ecuacion a estudiar es:
x + !02 x = 0:
Ensayando la solucion x(t) = et obtenemos la ecuacion caracterstica
2 + !02 = 0 )  = i!0 ;
luego, por el teorema 1.2.1, la solucion general es:
x(t) = c+ sen(!0 t) + c cos(!0 t) = A cos(!0 t Æ );
donde las constantes de integracion A; Æ representan la amplitud y el desfase, respectiva-
mente. Estas constantes se determinan a traves de las condiciones iniciales
s
  
) Æ = arctan( !vx0 ); A = x0 1 + xv!0 :
2
x(0) = x0
x_ (0) = v0 0 0 0 0

El resultado es un movimiento sinusoidal de periodo T = 2=!0. Las orbitas en el espacio

de las fases x_ (x) son elipses o vortices ya que
x(t)2 x_ (t)2
+ 2 2 = 1:
C2 C !0
 
2.1. PROBLEMAS DE V.I.: VIBRACIONES MECANICAS 
Y ELECTRICAS 23

2.1.2 Oscilador armonico amortiguado

Introduzcamos ahora una fuerza viscosa en el oscilador armonico simple, de forma que
la ecuacion a estudiar es:
x + 2 x_ + !02 x = 0:
Ensayando la solucion x(t) = et obtenemos la ecuacion caracterstica
q
 + 2  + ! = 0 )  =
2 2
0  2 !02 :
Caben diferentes posibilidades aqu dependiendo de los valores relativos de y !0 . Estu-
diemos cada uno por separado.

a) Oscilador armonico subamortiguado ( < !0 )

p
En este caso tenemos dos races complejas  =  i !02 2 =  i! y la
solucion general puede escribirse como:
q
x(t) = e t (c
1 sen(!t) + c2 cos(!t)) = e| {zt A} cos(!t Æ ); !  !02 2 :
A~(t)

Las constantes de integracion A; Æ (amplitud y desfase) se determinan a traves de las

condiciones iniciales
s
  2  
x(0) = x0 x0 + v0 x0 + v0
x_ (0) = v0 ) A = x0 1+
!x0
; Æ = arctan
!x0
:
p
El resultado es un movimiento oscilatorio de frecuencia !  !02 2 (menor que la

frecuencia natural !0 ) y amplitud decreciente A(t) = e A ! 0 cuando t ! 1 (vease la

~ t
gura 2.2). Las orbitas en el espacio de las fases x_ (x) son espirales (vease la gura 2.3).

b) Oscilador armonico crticamente amortiguado ( = !0 )

En este caso tenemos una raz real doble  = y la solucion general puede escribirse
como:
x(t) = c1 e t + c2 te t :
Las constantes de integracion c1 ; c2 se determinan a traves de las condiciones iniciales

x(0) = x0 ) c1 = x0 ; c2 = v0 + x0:
x_ (0) = v0
El resultado es un movimiento no oscilatorio en el que x(t) decrece con el tiempo de forma
exponencial (vease la gura 2.4). Notese que la partcula puede pasar por la posicion de
equilibrio a lo sumo una vez. Las orbitas en el espacio de las fases x_ (x) son nodos lmite.
24 CAPITULO 2. APLICACIONES A LA INGENIERIA CIVIL

x
t

Figura 2.2: La amplitud de las oscilaciones decrece con el tiempo para un oscilador subamor-
tiguado
v
x


Figura 2.3: Orbita espiral en el plano de fases \velocidad-posicion"para un oscilador suba-
mortiguado

c) Oscilador armonico sobreamortiguado ( > !0 )

p
En este caso tenemos dos races reales  =  2 !02 y la solucion general
puede escribirse como:
 p2 !02 t
p2 !02 t

x(t) = c e t + c+ e+ t = e t c+ e +c e ;
Las constantes de integracion c se determinan a traves de las condiciones iniciales

x(0) = x0
x_ (0) = v0 ) c+ = v0  x0 ; c =
v0 + + x0
+ 
:
+

El resultado es un movimiento no oscilatorio en el que x(t) decrece con el tiempo de forma

exponencial, al igual que en el caso crtico anterior. No obstante, la mayor viscosidad para
el oscilador sobreamortiguado hace que la amplitud vaya a cero mas lentamente (vease la
gura 2.5).
 
2.1. PROBLEMAS DE V.I.: VIBRACIONES MECANICAS 
Y ELECTRICAS 25
x
t

Figura 2.4: La amplitud de las oscilaciones decrece con el tiempo para un oscilador
crticamente amortiguado

x
t

ΒΩ0

t
ΒΩ0

Figura 2.5: Comparacion entre el oscilador crticamente amortiguado y sobreamortiguado

Las orbitas en el espacio de las fases x_ (x) son nodos.

2.1.3 Oscilador armonico forzado: pulsaciones y resonancia pu-

ra
Supongamos primeramente que no existe amortiguamiento e introduzcamos una fuerza
externa de tipo sinusoidal de la forma Fext (t) = F0 cos(!e t), donde !e denota la frecuencia
de la fuerza externa. La ecuacion a estudiar es:
F0
x + !02x = cos(!et):
m
Aplicando el metodo de los coe cientes indeterminados (teorema 1.2.2), ensayaremos una
la solucion particular del tipo xp (t) = A1 cos(!et) + A2 sen(!et). Introduciendo esta
solucion particular en la ecuacion diferencial obtenemos que A1 = !F020 =m
!e2 ; A2 = 0, con lo
26 CAPITULO 2. APLICACIONES A LA INGENIERIA CIVIL

cual, la solucion general de la ecuacion no homogenea es la suma de la solucion general

de la homogenea (s.g.h.) mas una particular de la no homogenea (s.p.n.h.):
F0 =m
x(t) = c|1 cos !0 t {z
+ c2 sen !0}t + cos !et :
!02 !e2
s:g:hh | {z }
s:p:n:h:
Caben estudiar aqu dos fenomenos interesantes dependiendo de los valores relativos de
la frecuencia natural !0 y la frecuencia externa !e .

a) Pulsaciones !e  !0

x
t

Figura 2.6: Fenomeno de las pulsaciones

Tomemos, por ejemplo, como condiciones iniciales

x(0) = 0 F0 =m
x_ (0) = 0 ) x(t) =
!02 !e2
(cos !e t cos !0 t)
   
F0 =m !0 !e !0 + !e
= 2 2 sen t sen t ;
!0 !e2 2 2
| {z }
A(t)
donde se ha utilizado la identidad trigonometrica: cos a cos b = 2 sen 12 (a b) sen 21 (a +
b). La solucion anterior se interpreta como una oscilacion rapida de frecuencia !0 +2 !e
(la media aritmetica de la frecuencia natural y la externa) modulada o envuelta por una
oscilacion lenta de frecuencia peque~na !0 2 !e cuando las frecuencias externa y natural son
muy parecidas !e  !0 (<pero distintas!). Vease gura 2.6.
Este fenomeno recibe el nombre de latidos o pulsaciones y es la base de la frecuencia
modulada (FM) en las emisoras de radio, donde el sonido audible (de menor frecuencia)
modula o envuelve a las ondas de radio (de alta frecuencia). Tambien es facilmente
detectable durante la \a nacion de una guitarra".
 
2.1. PROBLEMAS DE V.I.: VIBRACIONES MECANICAS 
Y ELECTRICAS 27

b) Resonancia pura !e = !0

x
t

Figura 2.7: Crecimiento lineal en el tiempo de la amplitud para una resonancia pura

Cuando la frecuencia externa !e coincide exactamente con la frecuencia natural del sis-
tema !0 (en ausencia de rozamiento), el metodo de los coe cientes indeterminados (teore-
ma 1.2.2) nos sugiere una solucion particular del tipo xp (t) = t(A1 cos(!e t) + A2 sen(!et)).
Introduciendo esta solucion particular en la ecuacion diferencial obtenemos que A1 =
0; A2 = F20!=m
e , con lo cual, la soluci
on general de la ecuacion no homogenea es la suma
de la solucion general de la homogenea (s.g.h.) mas una particular de la no homogenea
(s.p.n.h.):
F0 =m
x(t) = c|1 cos !0 t {z
+ c2 sen !0}t + t cos !0 t :
|
2 ! 0
{z }
s:g:hh
s:p:n:h:
Tomemos, por ejemplo, como condiciones iniciales
   F0
x(0) = x0 F0 =m v0 2m!
x_ (0) = v0 ) x(t) = x0 + 2! t cos !0t + ! 0 sen !0t:
0 0
| {z }
A(t)

En este caso, el coe ciente A(t) es lineal en el tiempo y ello supone que la amplitud de
las oscilaciones crezca inde nidamente (vease la gura 2.7).
Este fenomeno se denomina resonancia pura, y puede conllevar la rotura del resorte
cuando la amplitud de las oscilaciones alcanza el lmite elastico del material del que esta
fabricado dicho resorte. Esta es la vertiente \negativa"de la resonancia, la cual puede
presentarse como un efecto \pernicioso"en las vibraciones de ciertas estructuras en cons-
truccion. No obstante, este fenomeno tambien subyace a la fabricacion de sintonizadores,
como los de una radio, que \ ltranuna determinada frecuencia (vease seccion 2.1.5 sobre
la aplicacion a la teora de circuitos).
28 CAPITULO 2. APLICACIONES A LA INGENIERIA CIVIL

2.1.4 Oscilador armonico amortiguado y forzado: factor de am-

pli cacion

Ρ
r
4
Β1
3 1
Β 
2
2 1
Β 
4
1
1 Β 
8
r
1 2 3 4 5

Figura 2.8: Dependencia del factor de ampli cacion (r) con r  !!0e para distintos valores
de viscosidad

Supongamos ahora que existe amortiguamiento e introduzcamos una fuerza externa

de tipo sinusoidal de la forma Fext (t) = F0 cos(!e t), donde !e denota la frecuencia de la
fuerza externa. La ecuacion a estudiar es:
F0
x + 2 x_ + !02 x = cos(!e t):
m
Aplicando el metodo de los coe cientes indeterminados (teorema 1.2.2), ensayaremos una
la solucion particular del tipo xp (t) = A1 sen(!e t) + A2 cos(!et). Introduciendo esta
solucion particular en la ecuacion diferencial anterior obtenemos que
2 !eF0 (!02 !e2 )F0 =m
xp (t) = sen(! e t ) + cos(!e t)
(!02 !e2 )2 + 4 2!e2 (!02 !e2 )2 + 4 2 !e2
F0 =m F0
= p 2 cos( ! e t Æ ) =  (r ) cos(!e t Æ );
(!0
|
!e2 )2 + 4 2 !e2
{z } | {z k}
A(!e ) A(r)

donde se ha de nido:
1 !e
Factor de ampli cacion : (r) = p ; r  !0
;
(1 r2 )2 + 4 2 r2
 
2 !e
Desfase : Æ = arctan 2 :
!0 2 !e
 
2.1. PROBLEMAS DE V.I.: VIBRACIONES MECANICAS 
Y ELECTRICAS 29

En la gura 2.8 tenemos un gra co de la dependencia del factor de ampli cacion 

en funcion del cociente r  !!0e . Observamos que (r) presenta un maximo para r =
1 ) !e = !0 , o sea, cuando la frecuencia de la fuerza externa coincide con la frecuencia
natural de oscilador del resorte. Dicho maximo es mas pronunciado conforme el coe ciente
de rozamiento es mas peque~no, tendiendo a la resonancia pura (estudiada en el apartado
anterior) en el lmite ! 0. En este caso, la amplitud A(!e ) de las oscilaciones es maxima.
Es decir: 
Si
!e  !0 ) A(!e) >> 1 ! Resonancia:
<< 1
La solucion general de la ecuacion no homogenea sera la suma de la solucion general de

x
t

Ωe 4Ω0
t
Ωe Ω0

Figura 2.9: Comparacion de la amplitud para un movimiento amortiguado ( = 0:3) y forzado

con !e = 4!0 y en resonancia !e = !0

la homogenea (s.g.h.) mas la solucion particular de la no homogenea (s.p.n.h.) calculada

anteriormente. Es decir:

xg:n:h:(t) = xg:h: (t) + xp (t) = |e t B cos(!t + ') + A(r ) cos(! t Æ});

e
{z } | {z
Parte transitoria Parte estacionaria

donde B y ' son constantes arbitrarias a jar por las condiciones iniciales x(0) = x0 y
x_ (0) = v0 ; recuerde
p que la frecuencia para el movimiento amortiguado no forzado se de na
como: !  !02 on general de la homogenea xg:h:(t) se ve fuertemente
2 . La soluci
t
amortiguada por el factor e , de manera que es insigni cante pasado un cierto tiempo
t >> 1= y supone, por lo tanto un termino transitorio o solucion transitoria. Pasado ese
tiempo, la solucion general de la no homogenea tendera a la solucion particular, es decir,
xg:n:h:(t)  xp (t); t >> 1= , que se denomina parte estacionaria o solucion de estado
estacionario.
30 CAPITULO 2. APLICACIONES A LA INGENIERIA CIVIL

2.1.5 Analogas electricas. Leyes de Kirchho

La ecuacion diferencial que modela las vibraciones electricas (variacion de la carga
q (t) en el tiempo) en un circuito RCL, con una resistencia de valor R, un condensador
de capacidad C , una bobina de inductancia L y una fuente de alimentacion de fuerza
electromotriz E (t) en serie, viene dada por la segunda ley de Kirchho (\la suma de las
cadas de tension o voltajes en los elementos de un circuito debe ser igual a la fuerza
electromotriz aplicada"):
1
L q
 + R q
_ + q = E (t);
|{z} |{z} C
|{z}
VL VR
VC
donde VL = Lq denota la cada de tension en la bobina, VR = Rq_ la cada de tension en
la resistencia y VC = C1 q la cada de tension en el condensador. Todo lo dicho para las
vibraciones mecanicas puede traducirse directamente a las electricas sin mas que hacer la
siguiente identi cacion:
CONCEPTO MECANICO  ! CONCEPTO ELECTRICO 
desplazamiento x(t) ! carga q (t)
masa m ! inductancia L
viscosidad c ! resistencia R
cte. resorte k ! capacitancia 1=C
fuerza externa Fext (t) ! fuerza electromotriz E (t)

i
L
+

E (t) C

Figura 2.10: Circuito RCL

2.2 Problemas de V.F.: exion y pandeo en vigas

En los problemas de valor en la frontera que seguidamente exponemos se combinan
conceptos de interes en Geometra, Mecanica de Medios Continuos y Teora de Distribu-
  Y PANDEO EN VIGAS
2.2. PROBLEMAS DE V.F.: FLEXION 31

ciones como: curvatura, tensor de esfuerzos y deformaciones, Ley de Hooke, distribucion

de cargas, etc.

2.2.1 Flexion en vigas: ecuacion diferencial de la curva elastica

y echa de exion

 x -
Y6
Fibra neutra

X-


| {z }
Elemento diferencial de viga

Figura 2.11: Viga encorvada y bra neutra

En la gura 2.11 mostramos el aspecto general que presenta una viga horizontal al
encorvarse bajo las cargas que soporta, las cuales se suponen verticales. Imaginemos
la viga descompuesta en delgadas laminas horizontales. Debido al encorvamiento, las
laminas situadas en la region superior de la viga se encuentran comprimidas, en tanto
que las situadas en la region inferior estan estiradas. Ambas regiones estan separadas por
una capa cuyas bras no estan ni estiradas ni comprimidas; esa capa recibe el nombre
de capa o super cie neutra. Segun la gura 2.11, la interseccion de la super cie neutra
con el plano XY nos de ne una curva llamada bra neutra. Consideremos ahora un
elemento diferencial longitudinal de la viga, tal como se ilustra en la gura 2.12, sometido
a momentos exores de signo opuesto M (x + dx) = M (x dx) (para que el elemento no
gire) en los extremos. El elemento de viga experimenta una exion tal que su bra neutra
toma la forma de un arco de circunferencia de radio R. Otra bra situada a una distancia
y por encima de la neutra tendra un radio R y ; entonces, la deformacion longitudinal
unitaria (o porcentaje en cambio de longitud), de nida como el cociente: ((longitud de
la bra contrada) - (longitud de la bra neutra))/(longitud de la bra neutra), sera:
(R y )d Rd y
xx = = ; (2.2)
Rd R
donde d es el angulo substendido por el elemento de viga de longitud dx en el centro de
curvatura C . El signo negativo en la expresion (2.2) nos indica que las bras situadas por
32 CAPITULO 2. APLICACIONES A LA INGENIERIA CIVIL

C


d

M (x dx) Y
j R 6
9M x ( + dx)
Fibra
jj 
)
jj
q
 1
11
Fibra neutra
YY i 
y
-
xx (x; y )
X

Figura 2.12: Elemento longitudinal diferencial de viga y esfuerzo normal y momento exor
que soporta debido a la accion del material situado a izquierda y derecha

encima de la bra neutra (y > 0) estan comprimidas (xx < 0) en tanto que las situadas
por debajo (y < 0) estan estiradas (xx > 0). La ley de Hooke generalizada dice que la
deformacion longitudinal unitaria xx (x; y ) es proporcional al esfuerzo normal xx (x; y )
(compresor o tensor) que actua en la direccion del eje X sobre una bra situada a una
altura y por encima de la bra neutra:
E
xx (x; y ) = Exx (x; y ) = y;
R(x)
donde la constante de proporcionalidad E es el modulo de Young del material.
La fuerza normal total que actua sobre una seccion transversal S (en el plano Y Z )
de la viga como resultado de la accion del material situado a la izquierda sobre el situado
a la derecha de la seccion debe ser nula:
Z Z
E
F = xx dydz = ydydz = 0;
S R(x) S
ya que estamos suponiendo que todas las cargas que soporta la viga son verticales, de
modo que noR puede haber fuerza neta horizontal en ninguna seccion transversal de la viga.
Puesto que S ydydz = 0, sera necesario tomar el origen del eje vertical Y coincidiendo
con el centro de area o centroide de la seccion transversal S para cada valor de x, lo
que equivale a decir que la bra neutra pasa por los centroides de todas las secciones
transversales.
Aun cuando la distribucion de esfuerzos normales xx sobre una seccion transversal
de la viga representa una fuerza normal neta nula, no pasa lo mismo con el momento de
  Y PANDEO EN VIGAS
2.2. PROBLEMAS DE V.F.: FLEXION 33

dicha distribucion:
Z Z
E EIz
M (x) = yxxdydz = y 2dydz = ;
S R(x) | S {z } R(x)
Iz
R
donde Iz  S y 2dydz es el analogo del \momento de inercia "de la seccion transversal,
con relacion al eje Z que pasa por la bra neutra, supuesta dicha seccion una \lamina
delgada de densidad super cial unidad".
Se trata de determinar la forma de la bra neutra o curva elastica y (x). Recordemos
que, de acuerdo con la Geometra Diferencial, la curvatura  (el inverso del radio de
curvatura R) de una curva plana y (x) en un punto generico x de la misma se calcula
como:
1 y 00 (x)
(x) = =
R (1 + y 0(x)2 )3=2
 y00(x);
donde se ha supuesto que, como sucede an la mayor parte de los casos de interes practico, la
exion que experimenta la viga es muy peque~na, de manera que y 0(x)2 puede considerarse
despreciable frente a la unidad. As, la ecuacion diferencial de la curva elastica queda:
1
y 00 (x) = M (x):
EIz
Teniendo en cuentaR que la derivada del momento exor M (x) es igual a menos la fuerza
cortante F (x) = S xy dydz , y que la derivada de la fuerza cortante es igual a la densidad
de carga q (x) (peso por unidad de longitud),
M 0 (x) = F (x); F 0 (x) = q (x);
podemos representar tambien la ecuacion diferencial de la curva elastica como:
1
y IV (x) = q (x); (2.3)
EIz
donde y IV denota la derivada cuarta.
Hay que distinguir entre:
1. Cargas distribuidas, dF (x) = q (x)dx = w(x)dx, tales como el propio peso de la viga,
que se aplican de forma uniforme a lo largo de la viga (en general, funciones w(x)
continuas a trozos).
P
2. Cargas concentradas, dF (x) = q (x)dx = i Wi Æ (x xi )dx, tales como las reacciones
verticales Vi en los apoyos, que se aplican como cargas Wi localizadas en los puntos
xi .
Hemos introducido la \distribucion"(\<que no funcion!") delta de Dirac Æ (x xi )dx para
dar un tratamiento matematico uni cado de cargas distribuidas y concentradas como
distribuciones dF (x) de peso sobre la viga. Fsicamente, las cargas concentradas Wi estan
34 CAPITULO 2. APLICACIONES A LA INGENIERIA CIVIL

localizadas en un entorno peque~no de anchura a alrededor de los puntos de aplicacion xi , de

manera que (\a efectos practicos") las densidades de carga concentrada pueden escribirse
como

X
a;
1
si x 2 [xi a; x + a]
i 2
w(x) = Wi Æa (x xi ); donde Æa (x xi ) = 0; si x 62 [xi
2
a; x + a] ;
i 2 i 2

de manera que el area que substiende la \funcion meseta"Æa (x xi ) es siempre igual a

uno (vease gura 2.13). La abstraccion matematica consiste en hacer tender a ! 0, es
Æ0 (x)
6 Æa3 (x)

Æa2 (x)

Æa1 (x)

a1 a2 a3 a3 a2 a1
- x

Figura 2.13: Delta de Dirac como \lmite de funciones meseta"

decir, localizar mas y mas las cargas Wi hasta que estas actuen en un solo punto xi . As,
la delta de Dirac tiene sentido como una distribucion, al hacer promedios (integrales),
cumpliendo que, para una funcion f (x) continua en un punto x0 se tiene:
Z
f (x)Æ (x x0 )dx = f (x0 ):

Todas estas de niciones se pueden formalizar aun mas y esto constituye parte de la doc-
trina denominada Teora de Distribuciones. Para nuestros nes sera su ciente con saber
que la fuerza cortante en un punto x de la viga puede calcularse como:
Z x Z x X Z x X
F (x) = F (0) + q (x)dx = (w(l) + Wi Æ (l xi ))dx = w(l)dl + Wi ;
0 0 i 0 x>xi

donde hemos puesto F (0)  W0 y x0 = 0. Conseguimos as rebajar en uno el orden de la

ecuacion diferencial (2.3), quedando:
1
y 000 (x) = F (x);
EIz
  Y PANDEO EN VIGAS
2.2. PROBLEMAS DE V.F.: FLEXION 35

Recordando igualmente que M 0 (x) = F (x), el momento exor M puede calcularse como
la integral:
Z x X Z x
M (x) = M (0) F (l)dl = M (0) + xi Wi + lw(l)dl xF (x);
0 x>xi 0

que introducida en la ecuacion diferencial anterior rebaja el orden a dos:

1
y 00 (x) = M (x);
EIx
Recuerde que F (x) tiene discontinuidades de salto en las posiciones xi donde estan
localizadas las cargas concentradas Wi . Por lo tanto la curva elastica y (x) sera, en gene-
ral, una funcion continua con derivada segunda y 00 continua (es decir, de clase C 2 ), con
puntos angulosos donde estan localizadas las cargas concentradas. Si no existen cargas
concentradas, la funcion y (x) sera de clase C 3 .
Y
6 Condiciones de contorno:
` X
...................................... -
Soporte simple y(0) = y(`) = y 00
(0) = y (`) = 0
00

.....................................
Viga empotrada y(0) = y(`) = y (0) = y (`) = 0
0 0

.....................................
Viga en voladizo
y(0) = y (0) = y (`) = y (`) = 0
0 00 000

Figura 2.14: Condiciones de contorno para vigas con soporte simple, empotradas y en voladizo

La solucion general de la ecuacion (2.3) depende de 4 constantes arbitrarias que pueden

jarse con condiciones de contorno (vease gura 2.14). Tenemos tres tipos de condiciones
de contorno, dependiendo del tipo de sujecion que tenga la viga en sus extremos (x = 0; `):
soporte simple, soporte interconstruido (empotrada) y en voladizo (vease gura 2.14).
Veamos el signi cado fsico de las mismas:
1. La condicion y 0(0) = 0 signi ca que la viga esta empotrada en el extremo x = 0.
2. La condicion y 00 (`) = 0 signi ca que la viga no esta sometida a momentos exores
en el extremo x = `, es decir, esta apoyada o libre.
3. La condicion y 000(`) = 0 signi ca que no existen fuerzas cortantes en el extremo
x = `, es decir, la viga esta libre en ese extremo.
36 CAPITULO 2. APLICACIONES A LA INGENIERIA CIVIL

A dichas condiciones de contorno deben a~nadirse las

condiciones de empalme : y (xi ) = y (x+i ); y 0(xi ) = y 0(x+i );
a izquierda y derecha de los puntos xi donde se aplican las cargas concentradas. Veamos
algunos ejemplos:

Viga con soporte simple

Y
6
x
z }| {
x=2
z }| { W

6 6 - X

?
V1 = !` + W=2 ? ( x;y) V2 = !` + W=2
!x

0 <x<`

Figura 2.15: Viga en soporte simple sometida a su propio peso y a una carga W concentrada
en x = `

Supongamos que tenemos una viga horizontal de longitud 2` que esta apoyada sobre
sus extremos (vease gura 2.15). Se trata de encontrar la ecuacion de su curva elastica
y (x) y su de exion maxima o echa de exion ymax cuando:
1. La viga esta sometida a su propio peso por unidad de longitud (carga) w =cte.
5w`4
Solucion: EIy (x) = w`x3 =6 wx4 =24 w`3 x=3; ymax = 24 EI
2. La viga esta sometida a su propio peso y a una carga W localizada en el centro
(vease gura 2.15). Solucion: EIy (x) = w(`x3 =6 x4 =24 `3 x=3) + W
4 3 12 (3`x
2
j`
xj 6` x + ` ); ; ymax = 24EI 6EI .
3 2 3 5w` W `
  Y PANDEO EN VIGAS
2.2. PROBLEMAS DE V.F.: FLEXION 37

Viga con soporte interconstruido

y6

 x2
x22
-
 x1 -  -
qM x21
 -  x2 2l - )
0O R M0 -x
6 6
wx?1 Q1 Q2
wx?2 ?
w`+2 W w`+2 W
W
Figura 2.16: Viga con soporte interconstruido sometida a su propio peso y a una carga W
concentrada en x = `=2

Una viga horizontal de longitud ` esta empotrada en sus extremos (vease gura 2.16).
Encontrar la ecuacion de su curva elastica y su de exion maxima o echa de exion
cuando:

1. La viga esta sometida 2 a su propio peso por unidad de longitud (carga) w =cte.
Solucion: EIy (x) = wx (2 `x `2
x 2
) ; y = y ( `=2) = w`4
24 max 384EI

2. La viga esta sometida a su propio peso y a una carga W localizada en el centro

(vease gura 2.16). Solucion:
 wx2 
242 (2`x `2 x2 ) + W
48 (4x
3
3`x2 ); 0  x  `=2
EIy (x) = wx
24(2`x `2 x2 ) + W
48 (`
3
6`2 x + 9`x2 4x3 ); `=2  x  `

w`4 + 2W `3
ymax = y (`=2) =
384EI
38 CAPITULO 2. APLICACIONES A LA INGENIERIA CIVIL

y 6

 `
-
M j
0
x
- ` x
2 -
-
x

6 O

Q(x; y )

V = w` + W ? R
w(` x)

?
W
Figura 2.17: Viga en voladizo sometida a su propio peso y a una carga W concentrada en
x=`

Viga en voladizo
Una viga horizontal de longitud ` esta apoyada en su extremo izquierdo y libre en
el derecho (vease gura 2.17). Encontrar la ecuacion de su curva elastica y su de exion
maxima o echa de exion cuando:
1. La viga esta sometida a su propio peso por unidad de longitud (carga) w =cte.
w (4`x3 6`2 x2 x4 ); y
Solucion: EIy (x) = 24 w`4
max = 8EI

2. La viga esta sometida a su propio peso y a una carga W localizada en el extremo

w (4`x3
libre (vease gura 2.17). Solucion: EIy (x) = 24 6`2 x2 x4 ) + W (x3 =6
`x2 =2); ; ymax = ( 8w`EI4 + W3EI` )
3

2.2.2 Pandeo en vigas: carga de Euler y modos de desviacion

En el siglo XVIII Leonhard Euler fue uno de los primeros matematicos en estudiar un
problema de valores propios al analizar como se curva una columna elastica sometida a
una fuerza axial de compresion (vease gura 2.18) Hemos visto que la relacion entre el
momento exor M (x) y el radio de curvatura R(x) en un punto x de la viga viene dada
por:
EIz
M (x) = :
R(x)
  Y PANDEO EN VIGAS
2.2. PROBLEMAS DE V.F.: FLEXION 39

6Y
y + dy ............................
y ..........................
..
..
.. ..
 - . .
x x + dx
 -X
F F

Figura 2.18: Pandeo de una viga sometida a una fuerza axial F de compresion

Supongamos primeramente que la curvatura es peque~na, de manera que R(1x)  y 00(x).

El momento exor M (x) a una distancia x del extremo izquierdo de la viga es igual al
producto de la fuerza F por el brazo de momento (es decir, la ordenada correspondiente):
M (x) = F y (x);
de manera que, para pandeos peque~nos la ecuacion a estudiar es:
F
y 00 (x) + y (x) = 0; y (0) = 0 = y (`)
EIz
que tiene soluciones:
r
n Fn
yn (x) = A sen(n x); n = = ; n = 1; 2; 3; : : : :
` EIz
Es decir, se trata de un de un problema de autovalores n (cargas crticas) y autofunciones
yn (modos de desviacion) cuyo signi cado fsico es el siguiente: la columna se desva solo
cuando la fuerza de compresion tiene uno de los valores:
 2 EIz 2
Fn = n ; n = 1; 2; 3; : : : :
`2
Estas fuerzas se llaman cargas crticas. La curva de de exion y (x) que corresponde a
2 z
la mnima carga crtica, F1 =  `EI2 (carga de Euler ), es y1 (x) = A sen(x=`), que se
conoce como primer modo de desviacion. Las curvas de de exion correspondientes a
cargas superiores Fn ; n = 2; 3; : : : se corresponden con modos de desviacion donde la
columna tiene algun tipo de restriccion fsica o gua en xn = Ln (vease la gura 2.19). La
columna no se desviara para otras cargas intermedias que no sean estas cargas crticas
(o autovalores) Fn . Este es un problema parecido al de la cuerda giratoria (vease la
referencia [Zill], pag. 228), donde la cuerda adopta ciertas formas (modos de desviacion)
para ciertas frecuencias crticas de rotacion.
40 CAPITULO 2. APLICACIONES A LA INGENIERIA CIVIL

6y
-x n=1

6
- n=2

6
- n=3

Figura 2.19: Modos de desviacion para cargas crticas F1 ; F2 y F3

Consideremos ahora curvaturas  no necesariamente peque~nas. Recordemos que la

curvatura de de na como el cociente:
1 d
= = ;
R(x) ds
con ds un diferencial de arco de curva y  = arctan(y 0(x)) (vease gura 2.18). Si derivamos
nuevamente con respecto a s la ecuacion anterior y tenemos en cuenta que dy=ds = sen ,
obtenemos nalmente que la ecuacion diferencial de la curva elastica es:
d2  F
= sen :
ds 2 EIz
Esta ecuacion es formalmente identica a la de un <pendulo simple!. Se trata de una
ecuacion no lineal de orden dos cuyo estudio se sale de los objetivos de esta memoria.
Trate el lector de imaginar la forma cualitativa de las curvas elasticas solucion de esta
ecuacion (utilice la intuicion fsica. . . ).

2.2.3 Derrumbe del puente colgante de Tacoma Narrows

El 1 de julio de 1940, el puente colgante (de unos 12 metros de ancho y casi un kilometro
de largo) \Tacoma Narrows", estado de Washington (EEUU), se termino y abrio al tra co.
Casi de inmediato se observo que cuando el viento soplaba en direccion transversal a la
carretera, originaba grandes oscilaciones verticales en la plataforma o \tablero". Estas
ondulaciones dieron nalmente al traste con el puente el 7 de noviembre de ese mismo a~no,
cuando una racha intensa de viento hizo que las oscilaciones aumentaran hasta niveles
nunca vistos. Pronto las oscilaciones provocaron un movimiento torsional del puente y su
rotura.
Se pidio a Theodor von Karman, conocido ingeniero, que determinara la causa del
derrumbe. El  y sus colaboradores [AKW] dictaminaron que el viento, al soplar perpendi-
cularmente a la carretera, se separaba formando vortices alternos arriba y abajo del tablero
  Y PANDEO EN VIGAS
2.2. PROBLEMAS DE V.F.: FLEXION 41

y con ello estableca una fuerza vertical que actuaba sobre el puente y que causo las osci-
laciones. Otras personas propusieron que la frecuencia de esa fuerza periodica coincida
con la frecuencia natural del puente, llegando a la resonancia, a las grandes oscilaciones y
a la destruccion, como se describe en la gura 2.7. Durante casi 50 a~nos se propuso que la
resonancia fue la causa del derrumbe del puente, aunque el grupo de von Karman lo re-
chazo diciendo que \es muy improbable que la resonancia con vortices alternos desempe~ne
una funcion importante en las oscilaciones de los puentes colgantes"[AKW].
Como se discutio en las secciones 2.1.3 y 2.1.4, la resonancia es un fenomeno lineal.
Ademas, para que se presente debe haber una coincidencia exacta entre la frecuencia de
la fuerza externa y la natural del puente (en ausencia de amortiguamiento), lo cual es
siempre improbable.
Si la resonancia no origino el derrumbe, >cual fue la causa?. Las explicaciones recientes
[LM, CM, DH] sostienen que fueron los efectos no lineales, y no la resonancia lineal, los
factores principales que provocaron las grandes oscilaciones. En su teora interviene la
siguiente ecuacion en derivadas parciales
EIz yxxxx = q (x) Mytt cyt ky + + Fext (x; t); (2.4)
con condiciones de contorno
y (0; t) = y (`; t) = yxx(0; t) = yxx(`; t) = 0;
para la de exion vertical y (x; t) (direccion positiva hacia abajo) de una viga de masa M
por unidad de longitud, modulo de Young E y momento de inercia de la seccion transversal
Iz , sometida a las siguientes fuerzas:
1. La distribucion de peso q (x) o carga por unidad de longitud.
2. La fuerza inercial Mytt .
3. La fuerza amortiguadora cyt causada por la viscosidad del aire.

4. La fuerza recuperadora o elastica ky = +
k maxfy; 0g = ky si y > 0 que
0 si y < 0
ejercen los cables sobre el tablero cuando dichos cables estan tensos (y > 0).
5. La fuerza externa Fext (x; t) debida a los vortices de von Karman.
Lazer y McKenna introducen la \fuerza recuperadora a un lado solo" ky + como un me-
canismo no lineal que puede dar cuenta de los efectos destructivos del viento en el colapso
del puente de Tacoma Narrows. De hecho, demuestran que pueden ocurrir respuestas
subarmonicas de gran amplitud para excitaciones de peque~na amplitud en un amplio es-
pectro de valores de los parametros (viscosidad, masa, etc). Este mecanismo de \banda
ancha"parece mas plausible que el basado en la resonancia.
Notese que la ecuacion (2.4) es una ecuacion diferencial en derivadas parciales (EDP)
y no una EDO. No obstante, pueden hacerse aproximaciones y separaciones de variables
42 CAPITULO 2. APLICACIONES A LA INGENIERIA CIVIL

que llevan esta EDP a una EDO. En el puente de Tacoma se observo que fue el primer
modo de oscilacion (sin nodos) el que aparentemente perdio estabilidad y dio lugar a las
oscilaciones torsionales, de manera que es plausible el considerar:
x x
y (x; t) = u(t) sen( ); Fext (x; t) = F0 sen(!t) sen( ); q (x) = W 0 sen(`= ):
` `
Despues de reescalamientos y ciertas simpli caciones se llega a la siguiente ecuacion dife-
rencial ordinaria de nida a trozos:
u + 2Æ u_ + ( + 1)u = W + B sen(!t); u > 0;
u + 2Æ u_ + u = W + B sen(!t); u < 0:
El estudio de esta ecuacion no lineal cae fuera del alcance de esta memoria. La busqueda
de soluciones y el estudio cualitativo y numerico de dicha ecuacion sigue siendo objeto de
investigaciones actualmente.
 Bibliograf
a

[NOR] S. Novo, R. Obaya y J. Rojo, Ecuaciones y sistemas diferenciales, McGraw-Hill,

Madrid (1995).
[VJ] V. Jimenez, Ecuaciones diferenciales, Universidad de Murcia (1999).
[Ayres] Frank Ayres, Ecuaciones diferenciales, McGraw Hill (1991).
[Zill] Denis G. Zill, Ecuaciones diferenciales con aplicaciones de modelado, Int. Thomson
Publishing (1997).
[Simmons] G. Simmons/J. Robertson, Ecuaciones diferenciales con aplicaciones y notas
historicas, McGraw-Hill (1993)
[Edwards-Penney] C.H. Edwards/D.E. Penney, Ecuaciones diferenciales elementales y
problemas con condiciones en la frontera, Prentice-Hall Hispanoamericana (1994).
[Makarenko] A. Kiseliov, M. Krasnov, G. Makarenko, Problemas de ecuaciones diferen-
ciales ordinarias, Editorial Mir Moscu (1984).
[Ortega] M.R. Ortega, Lecciones de Fsica: mecanica III, Tipografa Catolica SCA (1988).
[Ortiz] I. Ortiz Berrocal, \Elasticidad", Litoprint-Pricam S.A. (1985).
[AKW] O.H. Amann, T. von Karman and G.B. Wooddru , The Failure of the Tacoma
Narrows Bridge, Washington, DC: Federal Works Agency, 1941.
[LM] A.C. lazer and P.J. McKenna, Large amplitude periodic oscilations in suspension
bridges: some new connections with nonlinear analysis, SIAM Reviews 32, 537-578
(1990).
[CM] A.R. Champneys and P.J. McKenna, On solitary waves of a piecewise-linear sus-
pension bridge model, Preprint (1996).
[DH] S.H. Doole and S.J. Hogan, A piecewise linear suspension bridge model: nonlinear
dynamics and orbit continuation, Dynamics and Stability of Systems 11(1), 1996.

43

View publication stats