CAPÍTULO 3

SISTEMAS DE UN GRADO DE LIBERTAD

RESUMEN

Se inicia el capítulo con el estudio de Vibraciones Libres en un sistema de un grado de

libertad, para el efecto se analizan cuatro casos en función del valor del factor de
amortiguamiento. Estos casos son: sin amortiguamiento; subamortiguada; críticamente
amortiguada y sobre amortiguada. Las definiciones que se presentan son muy útiles en la
Ingeniería Sismo Resistente. En la parte final del capítulo se aplica el caso de vibración libre sin
amortiguamiento para modelar el comportamiento del impacto de las olas en una estructura,
durante un Tsunami.

Posteriormente se analiza el caso de Vibración Forzada con excitación armónica, como

una introducción al caso de cimentación de motores. Se estudia con bastante detenimiento el
factor de amplificación dinámica, se deducen sus ecuaciones y como aplicación, se presenta un
ejemplo muy práctico que consiste en determinar la frecuencia de vibración de un suelo; la
frecuencia de vibración de una estructura y luego se determina el factor de amplificación dinámica
en la estructura por efecto del suelo.

Luego se estudia la respuesta en el tiempo de un sistema de un grado de libertad ante

un escalón unitario; un pulso rectangular; la respuesta impulsiva y se termina con la Integral de
Convolución. De esta manera se orienta el estudio a encontrar la respuesta en el tiempo de un
sistema de un grado de libertad ante acciones sísmicas empleando la Integral de Duhamel. Como
caso de aplicación se halla la respuesta de una estructura de un piso sometida a un sismo
impulsivo como fue el de Northridge de 1994, con el propósito de ir conociendo cómo se
comportan las estructuras ante este tipo de sismos de corta duración.

Se presentan los siguientes programas en este capítulo: v_libre para el estudio de

vibraciones libres en sistemas de un grado de libertad; fad con el que se obtiene los factores de
amplificación dinámica y duhamel que halla la respuesta en el tiempo de un sistema de un grado
de libertad ante una acción sísmica empleando la Integral de Duhamel.

Finalmente, como caso práctico se analiza la rotura de vidrios de la fachada de un

edificio, modelando como un sistema de vibración forzada con excitación armónica y con
condiciones iniciales debido a que la losa que trabaja en voladizo tiene una deformación vertical,
la misma que se incrementó con las vibraciones producidas, por las máquinas, en la
compactación de la ampliación de la vía.
2 SISTEMAS DE UN GRADO DE LIBERTAD

3.1 RESPUESTA ANTE UN PULSO RECTANGULAR

Existen varias formas de encontrar la respuesta en el tiempo de un sistema de un grado

de libertad ante un pulso rectangular, una de ellas es resolviendo la ecuación diferencial como
se lo hace en el Ejemplo 1. Otra es por medio de la función de Heaviside que se verá más
adelante y una tercera forma es con métodos numéricos.

El estudio de la respuesta ante un pulso rectangular es importante para estudiar la

función delta de dirac.

 EJEMPLO 1

En un sistema de un grado de libertad sin amortiguamiento, actúa el pulso rectangular

indicado en la figura 3.1, si las condiciones iniciales son nulas, se pide:

i. Encontrar la ecuación que define la posición de la masa.

𝑘𝑔 𝑘𝑔 𝑠 2
ii. Si 𝐹0 = 2500 𝑘𝑔; 𝑘 = 1000 ; 𝑚 = 10 ; 𝑡1 = 0.9425 𝑠. Siendo 𝑡1 la duración del
𝑐𝑚 𝑐𝑚
pulso. Presentar la respuesta de desplazamientos.
iii. La relación entre el tiempo para el período del sistema vs el factor de amplificación 𝐷

𝑓(𝑡) = 𝐹0 0 ≤ 𝑡 ≤ 𝑡1

𝑓(𝑡) = 0 𝑡 ≥ 𝑡1

Figura 3.1 Pulso rectangular que actúa en el sistema de un grado de libertad.

 SOLUCIÓN

Fase 1 para 𝟎 ≤ 𝒕 ≤ 𝒕𝟏

𝑚 𝑞̈ + 𝑘 𝑞 = 𝐹0

La solución homogénea es:

𝑞ℎ = 𝐴 cos 𝑊𝑛 𝑡 + 𝐵 𝑠𝑒𝑛 𝑊𝑛 𝑡

La solución particular es:

𝐹0
𝑞𝑝 =
𝑘
 ROBERTO AGUIAR FALCONI 3
CEINCI-ESPE

𝐹0
𝑞 = 𝐴 cos 𝑊𝑛 𝑡 + 𝐵 𝑠𝑒𝑛 𝑊𝑛 𝑡 +
𝑘

𝑞̇ = − 𝐴 𝑊𝑛 sen 𝑊𝑛 𝑡 + 𝐵 𝑊𝑛 𝑐𝑜𝑠 𝑊𝑛 𝑡

 Condiciones iniciales

𝐹0
𝑡 = 0 𝑞0 = 0 →𝐴=−
𝑘

𝑡 = 0 𝑞̇ 0 = 0 → 𝐵 =0

 Respuesta en el tiempo para primera fase

𝐹0
𝑞(𝑡) = (1 − cos 𝑊𝑛 𝑡)
𝑘

Se recuerda que el desplazamiento lineal 𝑞𝑠𝑡 = 𝐹0 /𝑘

𝑞(𝑡) 𝐹0
= 1 − cos 𝑊𝑛 𝑡 𝑞𝑠𝑡 = 𝐷 = 1 − cos 𝑊𝑛 𝑡
𝑞𝑠𝑡 𝑘

Siendo 𝐷 el factor de amplificación por fuerza.

𝑞(𝑡) = 𝑞𝑠𝑡 𝐷

Condiciones Iniciales para segunda fase en 𝒕 = 𝒕𝟏

𝐹0
𝑞(𝑡1 ) = (1 − cos 𝑊𝑛 𝑡1 )
𝑘

𝐹0 𝑊𝑛
𝑞̇ (𝑡1 ) = 𝑠𝑒𝑛 𝑊𝑛 𝑡1
𝑘

Fase 2 para 𝒕 ≥ 𝒕𝟏

𝑚 𝑞̈ + 𝑘 𝑞 = 0

 Sea 𝜏 = 𝑡 − 𝑡1

 Condiciones iniciales

𝑞(𝜏 = 0) = 𝑞(𝑡1 )

𝑞 (𝜏 = 0) ̇ = 𝑞̇ (𝑡1 )

𝑞 (𝜏) = 𝐴 𝑠𝑒𝑛 𝑊𝑛 𝜏 + 𝐵 cos 𝑊𝑛 𝜏

𝑞̇ (𝜏 = 0)
𝑞(𝜏) = 𝑠𝑒𝑛 𝑊𝑛 𝜏 + 𝑞 (𝜏 = 0) cos 𝑊𝑛 𝜏
𝑊𝑛

𝑞̇ (𝑡1 )
𝑞(𝑡) = 𝑠𝑒𝑛 𝑊𝑛 (𝑡 − 𝑡1 ) + 𝑞 (𝑡1 ) cos 𝑊𝑛 (𝑡 − 𝑡1 )
𝑊𝑛
4 SISTEMAS DE UN GRADO DE LIBERTAD

𝐹0 𝐹0
𝑞(𝑡) = 𝑠𝑒𝑛 𝑊𝑛 𝑡1 𝑠𝑒𝑛 𝑊𝑛 (𝑡 − 𝑡1 ) + (1 − cos 𝑊𝑛 𝑡1 ) cos 𝑊𝑛 (𝑡 − 𝑡1 )
𝑘 𝑘

𝑞 (𝑡)
= 𝑠𝑒𝑛 𝑊𝑛 𝑡1 𝑠𝑒𝑛 𝑊𝑛 (𝑡 − 𝑡1 ) − cos 𝑊𝑛 𝑡1 𝑐𝑜𝑠𝑊𝑛 (𝑡 − 𝑡1 ) + cos 𝑊𝑛 (𝑡 − 𝑡1 )
𝑞𝑠𝑡

𝑞 (𝑡)
= cos 𝑊𝑛 (𝑡 − 𝑡1 ) − [cos 𝑊𝑛 𝑡1 𝑐𝑜𝑠𝑊𝑛 (𝑡 − 𝑡1 ) − 𝑠𝑒𝑛 𝑊𝑛 𝑡1 𝑠𝑒𝑛 𝑊𝑛 (𝑡 − 𝑡1 )]
𝑞𝑠𝑡

𝑞 (𝑡)
= cos 𝑊𝑛 (𝑡 − 𝑡1 ) − cos(𝑊𝑛 𝑡1 + 𝑊𝑛 𝑡 − 𝑊𝑛 𝑡1 )
𝑞𝑠𝑡

𝑞 (𝑡)
= cos 𝑊𝑛 (𝑡 − 𝑡1 ) − cos 𝑊𝑛 𝑡
𝑞𝑠𝑡

𝐷 = cos 𝑊𝑛 (𝑡 − 𝑡1 ) − cos 𝑊𝑛 𝑡

𝑞(𝑡) = 𝑞𝑠𝑡 𝐷

La respuesta en desplazamientos, se indica en la figura 3.1. Se recuerda que el sistema

no tiene amortiguamiento y que 𝑡1 = 0.942 𝑠.

EL período de vibración del sistema es 𝑇 = 0.6283 𝑠. En la figura 3.2 se indica la relación

𝑡
entre y el factor de amplificación 𝐷. Para el ejemplo desarrollado la relación entre el tiempo 𝑡1
𝑇
con el período 𝑇, es 1.5

Figura 3.1 Respuesta en desplazamientos.

 ROBERTO AGUIAR FALCONI 5
CEINCI-ESPE

𝑡
Figura 3.2 Relación entre vs el factor de amplificación 𝐷
𝑇

3.2 RESPUESTA ANTE FUNCIÓN DE HEAVISIDE

La función de Heaviside o función escalón unitario, tiene un valor de cero y luego a partir
del tiempo igual a cero el valor de la fuerza 𝑓(𝑡) = 1. En la figura 3.3, se presenta esta función
que es muy utilizada no solo en dinámica de estructuras sino también en la evaluación de la
peligrosidad sísmica para entender la probabilidad de excedencia.

Figura 3.3 Función escalón unitario.

Se desea encontrar la respuesta en el tiempo de un sistema de 1 gdl, ante una fuerza

f (t ) definida por la función de Heaviside; considerando condiciones iniciales nulas. En la figura
3.4 se presenta el modelo de análisis y la ecuación diferencial del movimiento es:
6 SISTEMAS DE UN GRADO DE LIBERTAD

.. .
m q  c q  k q  f (t )

Figura 3.4 Modelo de un grado de libertad sometido a una fuerza 𝑓(𝑡)

f (t )  1 t0
Se consideran nulas las condiciones iniciales. Luego: 𝑞0 = 0; 𝑞̇ 0 = 0

La ecuación diferencial a resolver es:

.. . .. .
 1
m q  c q  k q 1  m q  c q  k q    0
 k

Se realiza el siguiente cambio de variable:

1
q z
k
Luego la ecuación diferencial se transforma en:
.. .
m z c z k z 0

Por el cambio de variable, las condiciones iniciales, son:

1
𝑧0 = − 𝑧̇0 = 0
𝐾

Por lo tanto la solución se ha transformado en un problema de vibración libre con

condiciones iniciales que se estudió en el apartado 1.1. Se denomina g (t ) a la solución del
escalón unitario. Las soluciones son:

 Caso sin amortiguamiento

z (t )  A cosWn t   B senWn t 

z (t )   A Wn senWn t   B Wn cosWn t 
.

Al reemplazar las condiciones iniciales, se tiene:

1
 A 0B
k
Luego:
 ROBERTO AGUIAR FALCONI 7
CEINCI-ESPE

z (t )   cosWn t 
1
k
Con el cambio de variable se tiene:

  cosWn t  
1 1 1
q(t )  z (t ) 
k k k
q(t ) 
1
1  cosWn t 
k

A la solución se denomina g (t ) . Luego:

g (t ) 
1
1  cosWn t  (3.1)
k
 Caso sub amortiguado

Al proceder en forma similar al caso de vibración libre sin amortiguamiento se obtiene:

1 𝜉
𝑔(𝑡) = [1 − exp(−𝜉 𝑊𝑛 𝑡) (cos 𝑊𝑎 𝑡 + 𝑠𝑒𝑛 𝑊𝑎 𝑡)]
𝑘 √1 − 𝜉 2
(3.2)

Siendo 𝑊𝑎 = 𝑊𝑛 √1 − 𝜉 2 , frecuencia natural del sistema amortiguado.

 Caso sobre amortiguado

1     
g (t )  1  exp  Wn t  coshW a t  senhW a t  (3.3)
k   2 1 
  

W a  Wn  2  1 (3.4)

 Caso críticamente amortiguado

1 𝑊𝑛 1
𝑔(𝑡) = − 𝑡 exp(−𝑊𝑛 𝑡) − exp(−𝑊𝑛 𝑡) (3.5)
𝑘 𝑘 𝑘

Si la fuerza actuante no fuera unitaria sino que tiene una magnitud F0 la respuesta en
el tiempo, es:

𝑞(𝑡) = 𝐹𝑜 𝑔(𝑡)

 EJEMPLO 2

Deducir la ecuación 𝑔(𝑡) para el caso críticamente amortiguado, 𝜉 = 1.

 SOLUCIÓN

La ecuación diferencial en z y las condiciones iniciales son:

8 SISTEMAS DE UN GRADO DE LIBERTAD

𝑚𝑧̈ + 𝑐𝑧̇ + 𝑘𝑧 = 0
1
𝑧0 = − → 𝑧̇0 = 0
𝑘

Si la ecuación se divide para 𝑚 se tiene:

𝑐 𝑘
𝑧̈ + 𝑧̇ + 𝑧 = 0
𝑚 𝑚

Se conoce que la frecuencia natural 𝑊𝑛 , y el factor de amortiguamiento 𝜉 valen:

𝑘 𝑐
𝑊𝑛 = √ 𝜉=
𝑚 2√𝑚 𝑘

Luego la ecuación diferencial queda.

𝑧̈ + 2𝜉𝑊𝑛 𝑧̇ + 𝑊𝑛2 𝑧 = 0

Para 𝜉 = 1 , se tiene:

𝑧̈ + 2𝑊𝑛 𝑧̇ + 𝑊𝑛2 𝑧 = 0

Se plantea la solución de la ecuación diferencial de la siguiente forma:

𝑧(𝑡) = 𝑎 𝑒 𝜆𝑡

𝑧̇ = 𝑎 𝜆 𝑒 𝜆𝑡

𝑧̈ = 𝑎 𝜆2 𝑒 𝜆𝑡

Al reemplazar en la ecuación diferencial, se tiene:

𝑎 𝜆2 𝑒 𝜆𝑡 + 2 𝑊𝑛 𝑎 𝜆 𝑒 𝜆𝑡 + 𝑊𝑛2 𝑎 𝑒 𝜆𝑡 = 0

𝑎 𝑒 𝜆𝑡 (𝜆2 + 2 𝑊𝑛 𝜆 + 𝑊𝑛2 ) = 0

𝜆2 + 2 𝑊𝑛 𝜆 + 𝑊𝑛2 = 0

−2 𝑊𝑛 ± √4𝑊𝑛2 − 4𝑊𝑛2
𝜆=
2

𝜆1 = 𝜆2 = −𝑊𝑛

Por lo tanto la respuesta en el tiempo es:

𝑧(𝑡) = (𝐴𝑡 + 𝐵)exp(−𝑊𝑛 𝑡)

𝑧(𝑡) = 𝐴𝑡 exp(−𝑊𝑛 𝑡) + 𝐵 exp(−𝑊𝑛 𝑡)

𝑧̇ (𝑡) = 𝐴[exp(−𝑊𝑛 𝑡) − 𝑡 𝑊𝑛 exp(−𝑊𝑛 𝑡)] − 𝑊𝑛 B exp(−𝑊𝑛 𝑡)

 ROBERTO AGUIAR FALCONI 9
CEINCI-ESPE

Las condiciones iniciales son:

1
𝑡=0 𝑧𝑜 = − → 𝑧̇𝑜 = 0
𝑘
1
− =𝐵
𝑘
𝑊𝑛 𝑊𝑛
0=𝐴+ → 𝐴=−
𝑘 𝑘

Luego,
𝑊𝑛 1
𝑧(𝑡) = (− 𝑡 − ) exp(−𝑊𝑛 𝑡)
𝑘 𝑘
1
𝑔(𝑡) = + 𝑧(𝑡)
𝑘
1 𝑊𝑛 1
𝑔(𝑡) = −( 𝑡 + ) exp(−𝑊𝑛 𝑡)
𝑘 𝑘 𝑘
1 𝑊𝑛 1
𝑔(𝑡) = − 𝑡 exp(−𝑊𝑛 𝑡) − exp(−𝑊𝑛 𝑡)
𝑘 𝑘 𝑘

 EJEMPLO 3

Encontrar la respuesta en el tiempo para la fuerza f (t ) que se indica en la figura 3.5 en

que la fuerza empieza en el tiempo T y tiene una magnitud F0 .

Figura 3.5 Fuerza escalón de magnitud F0 .

 SOLUCIÓN

Para un tiempo t  T se tiene que el tiempo de duración de la fuerza F0 es t  T .

Luego:

q(t )  F0 g t  T 
10 SISTEMAS DE UN GRADO DE LIBERTAD

 EJEMPLO 4

Resolver el sistema de un grado de libertad, sin amortiguamiento y con condiciones

iniciales nulas, del Ejemplo 1, sometido a un pulso rectangular, ver figura 3.6. Pero por intermedio
de la función de Heaviside.

𝑘𝑔 𝑘𝑔 𝑠 2
𝐹0 = 2500 𝑘𝑔; 𝑘 = 1000 ; 𝑚 = 10 ; 𝑡1 = 0.9425
𝑐𝑚 𝑐𝑚

k
m f(t)

Figura 3.6 Sistema de un grado de libertad sin amortiguamiento.

Sometido a un pulso rectangular.

 SOLUCIÓN

Escriba aquí la ecuación.

.. .
m q  c q  k q  F0 0  t  t1
.
q(0)  q(0)  0

Se resuelve la ecuación diferencial indicada, considerando condiciones iniciales nulas,

.
después se halla la respuesta en q(T ) y q(T ) que son las condiciones iniciales de la siguiente
ecuación diferencial que es válida para t  T .
.. .
mq  c q  k q  0 t T
 ROBERTO AGUIAR FALCONI 11
CEINCI-ESPE

Figura 1.20 Artificio para resolver un pulso rectangular.

La segunda forma de solución se presenta en forma gráfica en la figura 1.20 en que el

pulso rectangular es igual a una fuerza escalón de magnitud F0 más otra fuerza escalón pero
de magnitud negativa F0 y que empieza en el tiempo T .

La solución para el caso indicado en la figura 1.20 es la siguiente:

𝑞(𝑡) = 𝐹𝑜 𝑔(𝑡) 0<𝑡<𝑇

(1.33)
𝑞(𝑡) = 𝐹𝑜 𝑔(𝑡) − 𝐹𝑜 𝑔(𝑡 − 𝑇) 𝑡≥𝑇

1.5 RESPUESTA IMPULSIVA

Si en el pulso rectangular de la figura 1.18 se considera que la fuerza 𝐹0 = 1/𝑇. De tal

manera que el área del pulso rectangular valga la unidad. El Impulso de la fuerza valdrá.

1
𝐼= 𝑇=1
𝑇

Ahora si 𝑇 decrece, el tamaño de la fuerza para que el área sea la unidad, tiene que
aumentar. Para cuando 𝑇 → 0, la fuerza aplicada tiende al ∞. Este es el caso en que se da un
golpe a un sistema, a esto se llama Acción de Percusión y el modelo numérico se denomina
delta de dirac 𝛿(𝑡)

0 𝑡≠0
𝛿(𝑡) = { }
∞ 𝑡=0

En este modelo se cumple, lo siguiente:

12 SISTEMAS DE UN GRADO DE LIBERTAD

∞ 0+
∫ 𝛿(𝑡) 𝑑𝑡 = ∫ 𝛿(𝑡) 𝑑𝑡 = 1
−∞ 0−

Donde 0− es un instante antes de que se aplique el golpe (martillazo) y 0+ es un instante

después de aplicarse el golpe.

1.5.1 Respuesta en ausencia de condiciones iniciales

Se define al Impulso como el producto de la Velocidad por la Masa, también se conoce

con el nombre de Momentum. Ahora bien se desea encontrar la respuesta en el tiempo de un
sistema de un grado de libertad, como el indicado en la figura 1.21 si la fuerza 𝑓(𝑡) es un golpe.
León (1981).

Figura 1.21 Descripción del problema ante un Impulso Unitario

Como la fuerza (martillazo) se aplica en un instante de tiempo tan pequeño, el

desplazamiento inicial es cero. Pero la velocidad inicial no y esta es igual al Impulso dividido para
la masa. Luego las condiciones iniciales que se generan por el golpe dado al sistema son:

1
𝑞(0) = 0 𝑞̇ (0) =
𝑚

En la ecuación (1.33) se tiene la respuesta en el tiempo para un pulso rectangular, esta

vale para 𝑡 > 𝑇

𝑞(𝑡) = 𝐹0 𝑔(𝑡) − 𝐹0 𝑔 (𝑡 − 𝑇)
1
Pero 𝐹0 = (Impulso Unitario) y además 𝑇 → 0. Al reemplazar todo esto en la expresión
𝑇
anterior se encuentra luego de factorar 𝐹0 .

[𝑔(𝑡) − 𝑔(𝑡 − 𝑇)]

𝑞(𝑡) = lim
𝑇→0 𝑇
𝑑 𝑔(𝑡)
𝑞(𝑡) =
𝑑(𝑡)

A la respuesta a un Impulso Unitario se le denomina ℎ(𝑡). Luego:

(1.34)
d g(t)
h(t)=
d(t)
 ROBERTO AGUIAR FALCONI 13
CEINCI-ESPE

1.5.2 Casos Particulares

En ausencia de condiciones iniciales, los casos particulares para cuando a un sistema

de un grado de libertad se da un golpe, la respuesta del sistema ℎ(𝑡) se halla derivando la función
𝑔(𝑡) con respecto al tiempo. Con lo que se halla.

 Sistema Sin Amortiguamiento 𝝃 = 𝟎

1 (1.35)
ℎ(𝑡) = 𝑆𝑒𝑛 (𝑊𝑛 𝑡)
𝑚 𝑊𝑛

 Sistema Subamortiguado 𝝃 < 1

1 (1.36)
h(t)= e-ξ Wn t Sen (Wa t)
m Wa

𝑊𝑎 = 𝑊𝑛 √1 − 𝜉 2

 Sistema Sobreamortiguado 𝝃 > 1

(1.37)
1
ℎ(𝑡) = ̂𝑎 𝑡)
𝑒 −𝜉 𝑊𝑛 𝑡 𝑆𝑒𝑛ℎ (𝑊
𝑚𝑊̂𝑎

̂a = Wn √ξ2 - 1
W

 Sistema Críticamente Amortiguado 𝝃 = 𝟏

𝑡
ℎ(𝑡) = exp(−𝑊𝑛 𝑡)
𝑚 (1.38)

1.6 INTEGRAL DE CONVOLUCIÓN

Sea 𝐹(𝑡) una función de fuerzas arbitrarias en función del tiempo, como la que aparece
a la izquierda de la figura 1.22, la misma que actúa en un sistema de un grado de libertad. Se
desea conocer la respuesta en el tiempo.
14 SISTEMAS DE UN GRADO DE LIBERTAD

f(t)

t
Δτ

Figura 1.22 Excitación arbitraria en el tiempo.

Para encontrar la respuesta del sistema, en el tiempo, la excitación se discretiza en una

serie de impulsos, cada uno de ellos tiene un intervalo de tiempo ∆𝜏. Mientras más pequeño es
el valor de ∆𝜏 los impulsos se van a aproximar a la función 𝐹(𝑡).

La Fuerza Impulsiva valdrá ∆𝜏 𝑓(𝜏). Siendo 𝑓(𝜏) el valor de la fuerza del pulso
rectangular que debe ser aproximada al valor de 𝐹(𝑡) en ese instante de tiempo. A la derecha de
la figura 1.22 se presenta la discretización del Primer Impulso y del segundo Impulso, debiendo
discretizarse en igual forma todos los demás impulsos

La respuesta de desplazamiento para un pulso cualquiera será: ∆𝜏 𝑓(𝜏) ℎ(𝑡 − 𝜏). Para
todos los pulsos la respuesta se convierte en una sumatoria.

𝑞(𝑡) = lim ∑ 𝑓(𝜏) ℎ(𝑡 − 𝜏) 𝑑𝜏

∆𝜏→0
 ROBERTO AGUIAR FALCONI 15
CEINCI-ESPE

𝑡
𝑞(𝑡) = ∫0 𝑓(𝜏) ℎ(𝑡 − 𝜏) 𝑑𝜏 (1.39)

 EJEMPLO

Presentar con cierto detalle la respuesta en el tiempo, aplicando la Integral de

Convolución, para un sistema de un grado de libertad sin amortiguamiento, ante una fuerza
arbitraria como la indicada en la figura 1.22. Considerar para la explicación un pulso rectangular
que se halla a un tiempo 𝜏

 SOLUCIÓN

Para el instante de tiempo 𝜏, la fuerza vale 𝑓(𝜏)𝑑𝜏. Si se tratará de un pulso rectangular

la respuesta es:

𝑓(𝜏)𝑑𝜏
𝑞(𝑡) = 𝑓(𝜏)𝑑𝜏 𝑔(𝑡 − 𝜏) = [1 − 𝑐𝑜𝑠𝑊𝑛 (𝑡 − 𝜏)]
𝑘

Ahora como no se tiene un Pulso sino un Impulso la respuesta es:

𝑑 𝑓(𝜏)
𝑞(𝑡) = ℎ(𝑡) = [𝑓(𝜏)𝑑𝜏 𝑔(𝑡 − 𝜏)] = 𝑠𝑒𝑛𝑊𝑛 (𝑡 − 𝜏)𝑑𝜏
𝑑𝑡 𝑚 𝑊𝑛

Finalmente, al sumar las contribuciones de todos los Impulsos, se encuentra que la

solución particular vale.
𝑡
𝑓(𝜏)
𝑞(𝑡) = ∫ 𝑠𝑒𝑛𝑊𝑛 (𝑡 − 𝜏)𝑑𝜏
𝑜 𝑚 𝑊𝑛

La solución total es igual a la homogénea más la particular, vale:

𝑡
𝑞̇ 0 𝑓(𝜏)
𝑞(𝑡) = 𝑠𝑒𝑛 𝑊𝑛 𝑡 + 𝑞0 cos 𝑊𝑛 𝑡 + ∫ 𝑠𝑒𝑛𝑊𝑛 (𝑡 − 𝜏)𝑑𝜏
𝑊𝑛 𝑜 𝑚 𝑊𝑛

 EJEMPLO

Presentar la solución particular y la solución total para un sistema sub amortiguado,

sometido a una fuerza arbitraria pero ya no con detalle sino únicamente sus resultados.

 SOLUCIÓN PARTICULAR
𝑡
1
𝑞(𝑡) = ∫ 𝑓(𝜏) 𝑒 −𝜉𝑊𝑛 (𝑡−𝜏) 𝑠𝑒𝑛 𝑊𝑎 (𝑡 − 𝜏) 𝑑𝜏
𝑚 𝑊𝑎 0

 SOLUCIÓN TOTAL
𝑡
𝑞̇ 0 + 𝜉 𝑊𝑛 𝑞0 1
𝑞(𝑡) = [𝑞0 cos 𝑊𝑎 𝑡 + 𝑠𝑒𝑛 𝑊𝑎 𝑡] + ∫ 𝑓(𝜏) 𝑒 −𝜉𝑊𝑛 (𝑡−𝜏) 𝑠𝑒𝑛 𝑊𝑎 (𝑡 − 𝜏) 𝑑𝜏
𝑊𝑎 𝑚 𝑊𝑎 0

𝑊𝑎 = 𝑊𝑛 √1 − 𝜉 2
16 SISTEMAS DE UN GRADO DE LIBERTAD

A la Solución particular se conoce con el nombre de Integral de Duhamel.

1.7 INTEGRAL DE DUHAMEL

Para el caso de que un sistema de 1 gdl., está sometida a una acción sísmica definida
por su acelerograma 𝑎(𝜏). La fuerza 𝑓(𝜏) = 𝑚 𝑎(𝜏). Por lo tanto, la solución particular para el
caso sub amortiguado queda.

1 𝑡 (1.49)
𝑞(𝑡) = ∫ 𝑎(𝜏) 𝑒 −𝜉 𝑊𝑛 (𝑡−𝜏) 𝑠𝑒𝑛 𝑊𝑎 (𝑡 − 𝜏)𝑑𝜏
𝑊𝑎 0

Al resolver las funciones exponencial y trigonométrica, indicada en la ecuación (1.49) se

obtiene:
𝑡
1 (1.50)
𝑞(𝑡) = ∫ 𝑎 (𝜏)𝑒 −𝜉 𝑊𝑛 𝑡 𝑒 𝜉 𝑊𝑛𝜏 [𝑠𝑒𝑛 𝑊𝑎 𝑡 cos 𝑊𝑎 𝜏 − cos 𝑊𝑎 𝑡 𝑠𝑒𝑛 𝑊𝑎 𝜏] 𝑑𝜏
𝑊𝑎 0

La ecuación (1.50) puede escribirse de la siguiente manera:

(1.51)
𝑒 −𝜉 𝑊𝑛𝑡
𝑞(𝑡) = [𝐴(𝑡)𝑠𝑒𝑛 𝑊𝑎 𝑡 − 𝐵(𝑡) cos 𝑊𝑎 𝑡]
𝑊𝑎

𝑡
(1.52)
𝐴(𝑡) = ∫ 𝑎 (𝜏)𝑒 𝜉𝑊𝑛 𝜏 cos 𝑊𝑎 𝜏 𝑑𝜏
0

𝑡
(1.53)
𝜉𝑊𝑛 𝜏
𝐵(𝑡) = ∫ 𝑎(𝜏)𝑒 sen 𝑊𝑎 𝜏 𝑑𝜏
0

Las integrales indicadas en las ecuaciones (1.52) y (1.53) pueden resolverse con
Métodos Numéricos pero en este capítulo, no se lo va a realizar de esa manera sino que se
pretende encontrar una solución analítica lo más exacta posible ya que en el próximo capítulo se
resuelve con el Método 𝛽 de Newmark.

Es lo más exácta posible ya que la eceleración del sismo, se la discretiza y se considera

que la variación entre dos instantes de tiempo 𝑡𝑖−1 y 𝑡𝑖 es lineal, como se muestra en la figura
1.21. Para un instante de tiempo en el intervalo indicado la variación de la aceleración 𝑎(𝜏) está
definida por la ecuación de una recta, con pendiente 𝑚

𝑎(𝑡𝑖 ) − 𝑎(𝑡𝑖−1 )
𝑚=
𝑡𝑖 − 𝑡𝑖−1
𝑎(𝜏) = 𝑎(𝑡𝑖−1 ) + 𝑚 (𝜏 − 𝑡𝑖−1 ) (1.54)
 ROBERTO AGUIAR FALCONI 17
CEINCI-ESPE

Figura 1.21 Variación lineal de la excitación.

Al reemplazar la ecuación (1.54) en (1.52) y (1.53), se obtiene las variables 𝐴, 𝐵 en el

tiempo discreto 𝑡𝑖 . (Barbat y Caicedo, 1992),

𝑒 𝜉 𝑊𝑛 𝜏
(𝜉 𝑊𝑛 𝑎(𝜏) + 𝑚 (1 − 2𝜉 2 )) cos 𝑊𝑎 𝜏 +
𝐴(𝑡𝑖 ) = 𝐴(𝑡𝑖−1 ) + [ 𝑊𝑎 ] (1.55)
𝑊𝑛2 (𝑊𝑛 𝑎(𝜏) − 2 𝜉 𝑚) 𝑠𝑒𝑛 𝑊𝑎 𝜏
𝑊𝑛

𝑒 𝜉 𝑊𝑛 𝜏
(𝜉 𝑊𝑛 𝑎(𝜏) + 𝑚 (1 − 2𝜉 2 )) sen 𝑊𝑎 𝜏 −
𝐵(𝑡𝑖 ) = 𝐵(𝑡𝑖−1 ) + [ 𝑊𝑎 ] (1.56)
𝑊𝑛2 (𝑊𝑛 𝑎(𝜏) − 2 𝜉 𝑚) 𝑐𝑜𝑠 𝑊𝑎 𝜏
𝑊𝑛

Los límites de integración de las ecuaciones (1.55) y (1.56) son desde 𝑡𝑖−1 hasta 𝑡𝑖

El programa Duhamel encuentra la respuesta en el tiempo de: desplazamiento,

velocidad y aceleración de un sistema de un grado de libertad sometido a una acción sísmica.
Los datos del programa, son:

 p Archivo que contiene las aceleraciones del sismo.

 dt Incremento de tiempo del archivo de aceleraciones.
 m Masa del sistema.
 k Rigidez del sistema.
 zeda Factor de amortiguamiento.

function Duhamel (p,dt,m,k,zeda)

%
% Programa para encontrar la respuesta en el tiempo
% de un sistema de un grado de libertad aplicando la
% Integral de Duhamel
%
% Por: Roberto Aguiar Falconi
% CEINCI-ESPE
% 12 de octubre de 2011
%-------------------------------------------------------------
% duhamel(p,m,k,zeda)
%-------------------------------------------------------------
% p: Archivo que contiene al acelerograma
% dt: Incremento de tiempo del acelerograma entre dos puntos discretos
% m: Masa del sistema de un grado de libertad
% k: Rigidez del sistema de un grado de libertad
% zeda: Factor de amortiguamiento
%-------------------------------------------------------------
w=sqrt(k/m);wa=w*(1-zeda*zeda); % Frecuencias sin y con amortiguamiento
np=length(p); % Numero de puntos del acelerograma
18 SISTEMAS DE UN GRADO DE LIBERTAD

q(1)=0;qp(1)=0; %Condiciones iniciales nulas

Asum=0; Bsum=0;
for i=2:np
a1=p(i-1); %Aceleración en tiempo discreto i-1
t1=dt*(i-1); %Tiempo
a2=p(i); %Aceleración en tiempo discreto i
t2=dt*i;t(i)=t2; %Tiempo
s=(a2-a1)/(t2-t1); % Pendiente de variación de la aceleración (sismo)
valA1=Aintegral(s,a2,t2,zeda,w,wa);valA2=Aintegral(s,a1,t1,zeda,w,wa);
Asum=Asum+valA1-valA2;
valB1=Bintegral(s,a2,t2,zeda,w,wa);valB2=Bintegral(s,a1,t1,zeda,w,wa);
Bsum=Bsum+valB1-valB2;
q(i)=exp(-zeda*w*t2)/wa*(Asum*sin(wa*t2)-Bsum*cos(wa*t2));
qp(i)=exp(-zeda*w*t2)*(Asum*cos(wa*t2)+Bsum*sin(wa*t2))-w*zeda*q(i);
qpp(i)=-w*w*q(i)-2*zeda*w*qp(i);
end
figure(1)
plot (t,q,'LineWidth',1.5);xlabel ('Tiempo'); ylabel ('Desplazamiento')
figure(2)
plot (t,qp,'LineWidth',1.5);xlabel ('Tiempo'); ylabel ('Velocidad')
figure (3)
plot(t,qpp,'LineWidth',1.5);xlabel ('Tiempo'); ylabel ('Aceleracion')
% fin
function [valor]=Aintegral(m,a,t,zeda,w,wa)
fac1=(zeda*w*a+m*(1-2*zeda*zeda))*cos(wa*t);
fac2=(wa/w)*(w*a-2*zeda*m)*sin(wa*t);
valor=exp(zeda*w*t)/(w*w)*(fac1+fac2);
return
function [valor1]=Bintegral(m,a,t,zeda,w,wa)
fac11=(zeda*w*a+m*(1-2*zeda*zeda))*sin(wa*t);
fac22=(wa/w)*(w*a-2*zeda*m)*cos(wa*t);
valor1=exp(zeda*w*t)/(w*w)*(fac11+fac22);
return

 EJEMPLO 11

Encontrar la respuesta en el tiempo, de la estructura indicada en la figura 1, ante el

registro obtenido en la estación Sylmar, durante el sismo de Northridge del 17 de enero de 1994,
que tuvo una magnitud de 6.7. El acelerograma tuvo una aceleración máxima de 826.76 gals. Se
considera un factor de amortiguamiento 𝜉 = 0.05

 SOLUCIÓN

Como complemento a la información sísmica se debe indicar que el registro de Sylmar

fue obtenido en un suelo de bastante resistente (Clasificación tipo B según el Servicio
Sismológico de los Estados Unidos con velocidades de la onda de corte comprendidas entre 360
y 750 m/s), la estación se encontraba a 18 km., de la zona epicentral. Es un sismo impulsivo
cuya fase intensa tuvo una duración de 2 seg.

4 (5 + 4) 𝑇 𝑠2 𝑇 𝑠2
𝑚= = 3.67 = 0.0367
9.8 𝑚 𝑐𝑚

Se obtuvo la matriz de rigidez lateral utilizando el programa rlaxinfi. Aguiar (2008). El

valor obtenido es:
 ROBERTO AGUIAR FALCONI 19
CEINCI-ESPE

𝑇 𝑇
𝑘 = 4388.6 = 43.886
𝑚 𝑐𝑚

Figura 1.22 Estructura de Ejemplo 8.

La rigidez lateral obtenida significa que si se aplica una fuerza lateral de 43.886 T., a la
altura del primer piso la estructura se desplaza 1 cm. En el capítulo IV, de este libro se presenta
el cálculo de la Matriz de Rigidez Lateral.

Se hizo el cambio de unidades, debido a que el acelerograma está en gals (cm/s2). Las
respuestas en el tiempo se obtuvieron con el Programa Duhamel y se obtuvieron las respuestas
que se indican en la figura 1.23.
20 SISTEMAS DE UN GRADO DE LIBERTAD

Figura 1.23 Acelerograma Sylmar. Respuestas de: desplazamientos; velocidad; aceleración.

 ROBERTO AGUIAR FALCONI 21
CEINCI-ESPE

Se presentan cuatro figuras, la primera es el acelerograma de Sylmar, a partir de los 20

seg., la aceleración del suelo es prácticamente nula (sismo). La segunda figura es la respuesta
de desplazamiento relativo de la estructura con respecto al suelo, nótese que el desplazamiento
máximo no llega a 1 cm., esto es consecuencia del sismo impulsivo que prácticamente no le da
tiempo a la estructura a desplazarse pero eso si le introduce una gran cantidad de energía por
esto la velocidad relativa es apreciable y que decir con la aceleración relativa de la estructura
con respecto al suelo que se presentan en la tercera y cuarta figura.

1.8 TRANSFORMADA DE FOURIER

Sea 𝑓(𝑡) una función en el dominio del tiempo. Se define la 𝐹(𝜃) la Transformada de
Fourier en el dominio de las frecuencias de la siguiente manera.
∞ (1.40)
𝐹(𝜃) = ∫−∞ 𝑓(𝑡) 𝑒 −𝑖𝜃𝑡 𝑑𝑡

Por otra parte, la transformada inversa de Fourier, viene definida por:

1 ∞
𝑓(𝑡) =
2𝜋
∫−∞ 𝐹(𝜃)𝑒 𝑖𝜃𝑡 𝑑𝜃 (1.41)

Donde 𝑓(𝑡) es una función continua en el dominio del tiempo; 𝐹(𝜃) es una función
continua en el dominio de las frecuencias; 𝜃 frecuencia angular analógica medida en radianes o
grados/seg; 𝑡 es la variable tiempo; 𝑖 = √−1.

1.9 FUNCIÓN DE TRANSFERENCIA

La ecuación diferencial del movimiento, de un sistema de un grado de libertad, ante
acciones sísmicas definidas por un acelerograma 𝑎(𝑡), es la siguiente.

𝑚 𝑞̈ + 𝑐 𝑞̇ + 𝑘 𝑞 = −𝑚 𝑎(𝑡) = 𝑓(𝑡) (1.42)

La transformada de Fourier de la excitación 𝑓(𝑡) viene definida por:

∞ ∞ (1.43)
𝐹(𝜃) = ∫−∞ 𝑓(𝑡) 𝑒 −𝑖𝜃𝑡 𝑑𝑡 = −𝑚 ∫−∞ 𝑎(𝑡) 𝑒 −𝑖𝜃𝑡 𝑑𝑡 = −𝑚 𝐴(𝜃)

Donde 𝜃 es la frecuencia de excitación; 𝐴(𝜃) la transformada de Fourier de la aceleración

𝑎(𝑡). Ahora, la transformada de Fourier de la respuesta 𝑞(𝑡), es:
∞
𝑄(𝜃) = ∫−∞ 𝑞(𝑡) 𝑒 −𝑖𝜃𝑡 𝑑𝑡 (1.44)

La transformada inversa de la respuesta como de la excitación, son:

1 ∞
𝑞(𝑡) =
2𝜋
∫−∞ 𝑄(𝜃)𝑒 𝑖𝜃𝑡 𝑑𝜃

∞ (1.45)
1
𝑓(𝑡) = ∫ 𝐹(𝜃)𝑒 𝑖𝜃𝑡 𝑑𝜃
2𝜋 −∞

Se define la función de transferencia 𝐻(𝜃) en el campo complejo como la relación entre

la respuesta (salida) 𝑄(𝜃) y la excitación (entrada) 𝐹(𝜃). Barbat y Canet (1994)
22 SISTEMAS DE UN GRADO DE LIBERTAD

𝑄(𝜃) (1.46)
𝐻(𝜃) =
𝐹(𝜃)
Por lo tanto, la respuesta en el dominio de las frecuencias viene dada por:

𝑄(𝜃) = 𝐻(𝜃) 𝐹 (𝜃) (1.47)

Al remplazar (1.47) en la primera ecuación de (1.45) se halla.

1 ∞
𝑞(𝑡) =
2𝜋
∫−∞ 𝐻(𝜃) 𝐹(𝜃)𝑒 𝑖𝜃𝑡 𝑑𝜃 (1.48)

1.10 EJERCICIOS COMPLEMENTARIOS

Con el propósito de reforzar ciertos temas que se han expuesto en este capítulo se
resuelven algunos ejemplos que pueden ser de interés para el lector.

 EJEMPLO 12

Las olas que destruyeron Japón, (11 de marzo de 2011, magnitud de 8.8), tenían una
altura de 10 m., cuando en la costa se produjo un movimiento sísmico de 2.5 Hz. Hallar la altura
ℎ a la cual se desplaza la plataforma para producir estas olas; establecer el modelo matemático
de lo sucedido sabiendo que los mares tiene un peso específico de 2.3 T/m3 (corteza más agua).
Si se desplaza 1 m3, de agua, que fuerza se generan en las estructuras, si se modela como un
caso de vibración libre sin amortiguamiento. (Merchán, 2011).

 SOLUCIÓN

En la figura 1.24 se presenta un modelo matemático muy sencillo sobre la generación de

un Tsunami. En el sismo se rompe la corteza y empieza a vibrar, esta vibración es la que genera
las olas las mismas que se desplazan hacia la costa, como se observa en la figura 1.24.

Figura 1.24 Modelo elemental para la generación de un Tsunami. (Tomado de Wikipedia).

Si el movimiento de la plataforma marina, se modela como un caso de vibración libre sin

amortiguamiento (aproximación), la altura ℎ es de 20 m. Debido a que 10 m., se desplaza hacia
arriba y 10 m., se desplaza hacia abajo, la plataforma.

El modelo numérico de vibración libre sin amortiguamiento, en sistemas de un grado de

libertad se ha presentado en este capítulo. La respuesta de desplazamientos, es la siguiente:

𝑞(𝑡) = 𝐶 𝑠𝑒𝑛 (𝑊𝑛 𝑡 + 𝜙)

Donde 𝐶 es la amplitud; 𝑊𝑛 es la frecuencia natural de vibración; 𝜙 es el ángulo de fase;

𝑡 es la variable tiempo. Se deriva dos veces para sacar la aceleración.
 ROBERTO AGUIAR FALCONI 23
CEINCI-ESPE

𝑞̇ (𝑡) = 𝐶 𝑊𝑛 cos(𝑊𝑛 𝑡 + 𝜙) ̈ = −𝐶 𝑊𝑛2 𝑠𝑒𝑛 (𝑊𝑛 𝑡 + 𝜙)

𝑞 (𝑡)

La fuerza que se genera en las estructuras es igual a masa por aceleración. Pero
interesa la fuerza máxima 𝐹 en valor absoluto, esta vale:

𝐹 = 𝑚 𝐶 𝑊𝑛2

1 ∗ 2.3 𝑇 𝑠2
𝑚= = 0.2347
9.8 𝑚

𝐹 = 0.2347 ∗ 10 ∗ 2.52 = 14.66 𝑇.

La fuerza de impacto en la estructura, es bastante alta. Ahora pensemos si la cantidad

de agua que llevó la ola es mayor o la frecuencia es mayor, la fuerza de impacto será mucho
más alta.

 EJEMPLO 13

Si la estructura del Ejemplo 8 se encuentra en el perfil de suelo indicado en la figura 1.25.

Encontrar el Factor de Amplificación Dinámica considerando la frecuencia natural del suelo y la
frecuencia natural de la estructura.

Figura 1.25 Perfiles de suelo en que se halla la estructura del Ejemplo 8.

 SOLUCIÓN

Con los datos indicados se obtiene en primer lugar la velocidad de la onda de corte del
perfil de suelo 𝑉𝑆 ; luego se encuentra el período de vibración suelo 𝑇𝑆 ; y posteriormente la
frecuencia de vibración del suelo 𝜔

ℎ1 ℎ2 𝐻
+ =
𝑉𝑠1 𝑉𝑠2 𝑉𝑠

4 2 6 𝑚
+ = ⟹ 𝑉𝑠 = 286.4
280 300 𝑉𝑠 𝑠

4𝐻 4∗6
𝑇𝑠 = = = 0.084 𝑠𝑒𝑔.
𝑉𝑠 286.4
24 SISTEMAS DE UN GRADO DE LIBERTAD

2𝜋 2𝜋 1
𝜔= = = 74.97
𝑇𝑠 0.084 𝑠

Por otra parte con los datos del ejemplo 8,se halla la frecuencia natural de la estructura.

𝑘 4388.6 1
𝑊𝑛 = √ = √ = 34.566
𝑚 3.673 𝑠

Finalmente, se obtiene la relación entre las frecuencias 𝑟 y el factor de amplificación

dinámica 𝛼. (Ecuación 1.25).

𝜔 74.97
𝑟= = = 2.17
𝑊𝑛 34.566

1 1
   0.28
1  r   2  r 
2 2 2
1  2.17   2  0.05  2.17
2 2 2

La estructura analizada no tiene problemas de amplificación de las ondas sísmicas por

efecto del tipo de suelo.

1.11 APLICACIÓN PRÁCTICA

El Edificio Administrativo de la Universidad Nacional del Chimborazo, en enero de 2012,
presentaba varias rajaduras en los vidrios de las ventanas de la fachada, que se indica en la
figura 1.26, debido a que se hallan en el extremo del voladizo y las vibraciones que se generaron
durante la ampliación de la vía a Guano.

Figura 1.26 Edificio Administrativo de la Universidad Nacional del Chimborazo .

En la figura 1.27 se presenta solamente la parte de la estructura que trabaja en voladizo;

las columnas se encuentran cada 5.0 m., la losa es de 25 cm., de peralte, alivianada con los
bloques que están indicados en la gráfica superior y el voladizo es de longitud variable. El gráfico
inferior corresponde al modelo de cálculo, que es una malla espacial, en la cual se considera que
las columnas están completamente empotradas y los restantes nudos tienen tres grados de
libertad, dos rotaciones y un desplazamiento vertical. Aguiar (2006).

Debido al peso de los elementos (carga muerta) y de las personas (carga viva) que en
ella laboran, los máximos desplazamientos verticales se producen entre los ejes: F-E-D. En la
figura 1.28 se indican los corrimientos verticales en el borde, los mismos que están en función
 ROBERTO AGUIAR FALCONI 25
CEINCI-ESPE

de la carga vertical que sobre ellos gravita, donde se tiene mayor cantidad de hormigón se
presentan desplazamientos más grandes, de ahí la importancia de no dejar vigas banda de
considerable dimensión en los extremos de un voladizo debido a que no se requiere y lo único
que ocasionan es un mayor desplazamiento vertical.

Figura 1.27 Macizado de la losa que trabaja en voladizo y malla de cálculo.

Figura 1.28 Desplazamientos verticales entre los ejes F-E

Se va a ver el comportamiento del nervio que tiene 2.65 mm., de deformación vertical y
que está con una línea más gruesa en la figura 1.28. Este nervio tiene 10 cm., de ancho y 25
cm., de peralte. Para el estudio de cómo se incrementan los desplazamientos verticales por las
vibraciones producidas por una máquina, se modela este elemento como un sistema de un grado
de libertad como el de la figura 1.29 sometido a una excitación armónica.

Figura 1.29 Vibración forzada armónica en un sistema de un grado de libertad.

26 SISTEMAS DE UN GRADO DE LIBERTAD

La ecuación diferencial para el caso de vibración forzada armónica, sin amortiguamiento

(hipótesis de cálculo) es la siguiente:

𝑚 𝑞̈ + 𝑘 𝑞 = 𝐹0 𝑠𝑒𝑛 Ω𝑡

Donde 𝑚, 𝑘 son la masa y rigidez del modelo; 𝑞, 𝑞̈ son el desplazamiento y aceleración

que es función del tiempo 𝑡, como se aprecia a la derecha de la figura 1.29; 𝐹0 , Ω es la fuerza y
frecuencia de la excitación que generan vibraciones en la estructura.

El sistema tiene condiciones iniciales y son la deformación vertical en el extremo del

voladizo por efecto de la carga vertical 𝑞0 , se considera que la velocidad es nula en 𝑡 = 0. La
solución de la ecuación diferencial es:

𝑞 = 𝑞ℎ + 𝑞𝑝

Donde 𝑞ℎ es la solución homogénea y 𝑞𝑝 es la solución particular, que valen.

𝑞ℎ = 𝐴 𝑠𝑒𝑛 𝑤𝑛 𝑡 + 𝐵 cos 𝑤𝑛 𝑡
𝐹0
𝑞𝑝 = 2 𝑠𝑒𝑛 Ω𝑡
(𝑘−𝑚 Ω )

𝐹0
𝑞(𝑡) = 𝐴 𝑠𝑒𝑛 𝑤𝑛 𝑡 + 𝐵 cos 𝑤𝑛 𝑡 + 2 𝑠𝑒𝑛 Ω𝑡
(𝑘 − 𝑚 Ω )

Donde 𝐴, 𝐵 son constantes de integración que se obtienen en función de las condiciones

iniciales; 𝑤𝑛 es la frecuencia de vibración del sistema.

𝑘 3𝐸𝐼
𝑤𝑛 = √𝑚 𝑘= 𝐿3

Donde 𝐸 es el módulo de elasticidad del hormigón; 𝐼 es el momento de inercia de la

viga; 𝐿 es la longitud del voladizo. Al reemplazar las condiciones iniciales en 𝑡 = 0, se hallan las
constantes de integración.
𝐹0 Ω
𝐴= − 2
(𝑘 − 𝑚Ω ) 𝑤𝑛

𝐵 = 𝑞0

Finalmente la respuesta en el tiempo del sistema es:

𝐹0 Ω 𝐹0
𝑞(𝑡) = − 2 𝑠𝑒𝑛 𝑤𝑛 𝑡 + 𝑞0 cos 𝑤𝑛 𝑡 + 2 𝑠𝑒𝑛 Ω𝑡
(𝑘 − 𝑚Ω ) 𝑤𝑛 (𝑘 − 𝑚 Ω )

En la tabla 1.2 se presentan los datos con los que se resuelve el problema de vibración
forzada con condiciones iniciales y en la figura 1.30 se indica la respuesta en el tiempo.

Tabla 1.2 Datos del nervio analizado, para el cálculo de las deformaciones verticales.
Dato Valor
Ancho y peralte del nervio 𝑏 = 10 𝑐𝑚., ℎ = 25 𝑐𝑚.
Condiciones iniciales 𝑞0 = 2.65 𝑚𝑚., 𝑞̇ 0 = 0
Longitud de nervio y masa puntual 𝑇 𝑠2
𝐿 = 1.65 𝑚. , 𝑚 = 0.1624
𝑚
Fuerza y Frecuencia de excitación armónica 𝐹𝑜 = 0.5 𝑇. , Ω = 0.5 Hertz
 ROBERTO AGUIAR FALCONI 27
CEINCI-ESPE

En la figura 1.30 se aprecia que la deformación vertical del extremo del nervio, antes de
que operen los equipos viales es 2.65 mm., y cuando entran en funcionamiento se incrementan
hasta tener alrededor de 6 mm., tanto hacia arriba como hacia abajo.

Por otra parte, el desplazamiento vertical máximo, para edificios que tienen elementos
no estructurales que pueden sufrir daño es. ACI (318 S-05).

𝐿
𝑞𝑚𝑎𝑥 =
480

Siendo 𝐿, la longitud del voladizo. En este caso 𝐿 = 1.65 𝑚., con lo que el desplazamiento
máximo permitido es 3.44 mm., cantidad que si fue sobrepasada con el movimiento vibratorio
durante la construcción de la vía.

El desplazamiento de 2.65 mm., que se obtuvo en el cálculo por efecto de las cargas
verticales, se da cuando termina la construcción, luego con el paso del tiempo, en los voladizos
estos desplazamientos verticales se incrementan por lo que se denomina Flujo Plástico del
Hormigón. Un libro clásico que trata muy bien este tema es el de Park y Paulay (1979)

Figura 1.30 Respuesta en el tiempo del desplazamiento vertical de un nervio mientras está en
funcionamiento las máquinas viales.

Por lo tanto, cuando se tiene una ventana bajo un voladizo se debe dejar una separación
entre el marco de la ventana y la viga o losa, para permitir esta deformación vertical instantánea,
en el transcurso del tiempo y que puede incrementarse con las vibraciones de equipos; esta
separación que puede estar entre 5 y 10 mm., debe ser llenada con un material amortiguante,
como una goma.

REFERENCIAS
28 SISTEMAS DE UN GRADO DE LIBERTAD

1. ACI (2005), Requisitos de Reglamento para Concreto Estructural y Comentario (ACI

318S-05). American Concrete Institute, 490 p.
2. Aguiar R., (2008), Análisis Sísmico de Edificios, Centro de Investigaciones Científicas.
Escuela Politécnica del Ejército, 322 p., Quito, Ecuador.
3. Aguiar R., Almazán J. L., Dechent P., Suárez V., (2008), Aisladores de base
elastoméricos y FPS, Centro de Investigaciones Científicas. Escuela Politécnica del
Ejército, 292 p., Quito, Ecuador.
4. Aguiar R., (2006), Análisis Estático de Estructuras, Centro de Investigaciones Científicas.
Escuela Politécnica del Ejército. Colegio de Ingenieros Civiles de Pichincha, 161 p.,
Quito.
5. Aguiar R., (1981), Vibraciones Mecánicas, Materia de Curso de Post Grado en Ingeniería
Sismo Resistente. Universidad Central de Venezuela, Caracas.
6. Barbat A., Canet J., (1994), Estructuras sometidas a acciones sísmicas. Cálculo por
ordenador, Centro Internacional de Métodos Numéricos en Ingeniería. Universidad
Politécnica de Cataluña. Segunda Edición, 809 p., Barcelona, España.
7. Barbat A., Caicedo C., (1992), “Dos soluciones exactas para la ecuación del movimiento”,
Trabajo presentado en el Curso de Post grado Ingeniería Sísmica y Dinámica Estructural.
Universidad Politécnica de Cataluña, 8 p., Barcelona, España.
8. Caiza P., (2009), Dinámica de Estructuras, Materia de Curso de Post grado en Ingeniería
Estructural. Universidad de Illinois en Urbana-Champaign, USA.
9. Constantinou M., (2013), Structural Dynamic and Earthquakes Engineering, Notas de
Clase del Curso dictado en 2013 en la Universidad de Buffalo, USA.
10. Chopra A., (1996), Dynamics of structures: Theory and applications to earthquakes
engineering, Prentice Hall, N.J.
11. Chopra A., (2001), Dynamics of structures: Theory and applications to earthquakes
engineering, Prentice Hall, N.J.
12. Jerome P. (2004), Structural Dynamic, Notas de clase del Curso 2004-2005 dictado en
la Universidad de Michigan, USA.
13. León J., (1981), Introducción a las Vibraciones Mecánicas, Materia del Curso de Post
grado en Ingeniería Sismo Resistente. Universidad Central de Venezuela, Caracas,
Venezuela.
14. Merchán P. (2011), Vibraciones Libres. Clase de Física en la Carrera de Mecatrónica.
Escuela Politécnica del Ejército, Quito, Ecuador.
15. Newmark N., and Hall W., (1982), Earthquake Spectra and Design. California, United
States of America: Earthquake Engineering Research Institute.
16. Park R., Paulay T., (1979), Estructuras de concreto reforzado, Editorial Limusa, 796 p.,
México.