CAPÍTULO 4

ANÁLISIS DE SISTEMAS DE 1 GDL CON

MÉTODOS NUMÉRICOS
RESUMEN

El Método 𝛽 de Newmark, es un clásico para encontrar la respuesta en el tiempo de un

sistema de un grado de libertad sometido a cualquier tipo de fuerza dinámica, sea esta pulsos o
el caso más importante la acción sísmica expresada mediante un registro de aceleraciones.
Además de ello todavía se lo utiliza no solo en el análisis lineal de estructuras sino también en el
análisis no lineal, con alguna variante.

Por este motivo es que se inicia el capítulo deduciendo las ecuaciones del Método 𝛽 de
Newmark, para tres casos de variación de la respuesta de aceleraciones que son: lineal,
constante y escalonada. En base a estos resultados se generaliza la respuesta con la variable 𝛽

Se aplica primero el Método 𝛽 de Newmark al caso en que la fuerza arbitraria viene

definida por pulso; se compara la respuesta con la solución exacta estudiada en el capítulo
anterior. Posteriormente se realizan varios ejemplos con pulsos rectangulares y triangulares.

Después se halla la respuesta en el tiempo de una estructura que es modelada como un

sistema de un grado de libertad, ante la componente N-S, obtenido en Manta en un perfil de
suelo C, del terremoto del 16 de abril de 2016. Se encuentra la respuesta en el tiempo y también
los valores máximos en valor absoluto.

En el Anexo A, se presentan algunos datos sismológicos del terremoto de Ecuador de

2016 y los epicentros, con sus áreas de ruptura de sismos registrados en la costa norte de
Ecuador, desde inicios del siglo XX, se observa que en términos redondos cada 20 años se
produce un sismo de magnitud mayor a 7, en la costa norte. Esta realidad se debe contrarrestar
con el diseño eficiente de las estructuras.

Para ese diseño eficiente se necesita contar con registros sísmicos de movimientos
fuertes de tectónica y condiciones de suelo similares a lugar del Proyecto, por esta razón se
presenta los mecanismos focales y epicentros de sismos interfase registrados en Ecuador (2) y
en Perú (4). Se indica en una tabla los nombres de los archivos que contienen los registros de
aceleración del suelo.
110 ANÁLISIS DE SISTEMAS DE 1 GDL CON MÉTODOS NUMÉRICOS

4.1 MÉTODO DE ACELERACIÓN LINEAL EN SISTEMAS DE 1GDL

El Método 𝛽 de Newmark, es un clásico para encontrar la respuesta en el tiempo de un

sistema de un grado de libertad, que se estudiará en este apartado, y de 𝑛 grados de libertad,
que se verá en otro capítulo posterior. Si 𝛽 = 1/6, la variación de la aceleración de la respuesta
es de tipo lineal, como se ilustra en la figura 4.1. Antes de proceder a la deducción numérica del
método, conviene indicar la ecuación diferencial del movimiento, en el rango elástico, para el
caso general es:

𝑚 𝑞̈ + 𝑐 𝑞̇ + 𝑘 𝑞 = 𝑄 (4.1)

Donde: 𝑚, 𝑐, 𝑘 son la masa, amortiguamiento y rigidez del sistema de un grado de

libertad; 𝑄 es la carga que actúa en el sistema; 𝑞, 𝑞,̇ 𝑞̈ , son el desplazamiento, velocidad y
aceleración. Se destaca que tanto 𝑄 como 𝑞, 𝑞,̇ 𝑞̈ son funciones del tiempo 𝑡. Dependiendo de la
forma de Q se puede pensar en tener una solución analítica exacta, si es que es factible pero
para el caso sísmico lo más fácil es pensar en una solución aproximada que se cumpla
únicamente en ciertos puntos discretos, en este caso la ecuación (4.1) toma la forma siguiente:

𝑚 𝑞̈ 𝑖 + 𝑐 𝑞̇ 𝑖 + 𝑘 𝑞𝑖 = 𝑄𝑖

Para el tiempo discreto 𝑖 + 1, se tiene:

𝑚 𝑞̈ 𝑖+1 + 𝑐 𝑞̇ 𝑖+1 + 𝑘 𝑞𝑖+1 = 𝑄𝑖+1

Si se resta estas dos últimas ecuaciones se tiene:

𝑚 Δ𝑞̈ + 𝑐 Δ𝑞̇ + 𝑘 Δ𝑞 = Δ𝑄 (4.2)

Donde:

Δ𝑞̈ = 𝑞̈ 𝑖+1 − 𝑞̈ 𝑖 Δ𝑞̇ = 𝑞̇ 𝑖+1 − 𝑞̇ 𝑖 Δ𝑞 = 𝑞𝑖+1 − 𝑞𝑖 Δ𝑄 = 𝑄𝑖+1 − 𝑄𝑖 (4.3)

Figura 4.1 Deducción del Método de Aceleración Lineal

Para el caso de aceleración lineal, se debe determinar la ecuación de la recta de la figura

4.1, para cuando el tiempo 𝑡𝑖 ≤ 𝑡 ≤ 𝑡𝑖+1 . Sea 𝑚𝑝 la pendiente de esta recta. (Lamar, 1981).

Δ𝑞̈ 𝑞̈ 𝑖+1 − 𝑞̈ 𝑖
𝑚𝑝 = =
Δ𝑡 𝑡𝑖+1 − 𝑡𝑖

La ecuación de una recta que pasa por un punto dado es:

 ROBERTO AGUIAR FALCONI 111

Δ𝑞̈
𝑞̈ (𝑡) − 𝑞̈ 𝑖 = (𝑡 − 𝑡𝑖 )
Δ𝑡
Δ𝑞̈ (4.4)
𝑞̈ (𝑡) = 𝑞̈ 𝑖 + (𝑡 − 𝑡𝑖 )
Δ𝑡

Al integrar la ecuación (4.4) entre el tiempo 𝑡𝑖 y un tiempo cualquiera 𝑡 que se encuentre

dentro del intervalo indicado en la figura 4.1, se tiene:
𝑡 𝑡 𝑡
Δ𝑞̈
∫ 𝑞̈ (𝑡) 𝑑𝑡 = ∫ 𝑞̈ 𝑖 𝑑𝑡 + ∫ (𝑡 − 𝑡𝑖 ) 𝑑𝑡
Δ𝑡
𝑡𝑖 𝑡𝑖 𝑡𝑖

Δ𝑞̈ (𝑡 − 𝑡𝑖 )2 (4.5)
𝑞̇ (𝑡) − 𝑞̇ 𝑖 = 𝑞̈ 𝑖 (𝑡 − 𝑡𝑖 ) +
Δ𝑡 2

Para el tiempo 𝑡 = 𝑡𝑖+1 se tiene:

Δ𝑞̈ (𝑡𝑖+1 − 𝑡𝑖 )2
𝑞̇ 𝑖+1 − 𝑞̇ 𝑖 = 𝑞̈ 𝑖 (𝑡𝑖+1 − 𝑡𝑖 ) +
Δ𝑡 2
Δ𝑡 (4.6)
Δ𝑞̇ = 𝑞̈ 𝑖 Δ𝑡 + Δ𝑞̈
2

Ahora, al integrar la ecuación (4.5) se obtiene: en primer lugar la ecuación que

proporciona el desplazamiento para el intervalo 𝑡𝑖 ≤ 𝑡 ≤ 𝑡𝑖+1 ; luego se halla el desplazamiento
en 𝑖 + 1 y finalmente el incremento de desplazamientos Δ𝑞.
𝑡 𝑡 𝑡 𝑡
Δ𝑞̈ (𝑡 − 𝑡𝑖 )2
∫ 𝑞̇ (𝑡)𝑑𝑡 = ∫ 𝑞̇ 𝑖 𝑑𝑡 + ∫ 𝑞̈ 𝑖 (𝑡 − 𝑡𝑖 ) 𝑑𝑡 + ∫ 𝑑𝑡
Δ𝑡 2
𝑡𝑖 𝑡𝑖 𝑡𝑖 𝑡𝑖

(𝑡 − 𝑡𝑖 )2 Δq̈ (𝑡 − 𝑡𝑖 )3
𝑞(𝑡) − 𝑞𝑖 = 𝑞̇ 𝑖 (𝑡 − 𝑡𝑖 ) + 𝑞̈ 𝑖 +
2 Δ𝑡 6

El desplazamiento en 𝑡 = 𝑖 + 1 vale:

Δ𝑡 2 Δ𝑞̈ Δ𝑡 2
𝑞𝑖+1 − 𝑞𝑖 = 𝑞̇ 𝑖 Δ𝑡 + 𝑞̈ 𝑖 +
2 6
De donde:

Δ𝑡 2 Δ𝑞̈ Δ𝑡 2 (4.7)
Δ𝑞 = 𝑞̇ 𝑖 Δ𝑡 + 𝑞̈ 𝑖 +
2 6

Ahora al reemplazar los incrementos de velocidad y desplazamiento en la ecuación (4.2)

se tiene:

Δ𝑞̈ Δ𝑡 2 Δ𝑞̈ Δ𝑡 2
𝑚 Δ𝑞̈ + 𝑐 (𝑞̈ 𝑖 Δ𝑡 + Δ𝑡) + 𝑘 (𝑞̇ 𝑖 Δ𝑡 + 𝑞̈ 𝑖 + ) = Δ𝑄
2 2 6

𝑐 𝑘 Δ𝑡 2 𝑞̈ 𝑖
(𝑚 + Δ𝑡 + ) Δ𝑞̈ = Δ𝑄 − 𝑐 Δ𝑡 𝑞̈ 𝑖 − 𝑘 𝑞̇ 𝑖 Δ𝑡 − 𝑘 Δ𝑡 2
2 6 2
 112 ANÁLISIS DE SISTEMAS DE 1 GDL CON MÉTODOS NUMÉRICOS

Sea 𝑀∗ la masa equivalente del sistema, definida de la siguiente manera:

𝑐 𝑘 Δ𝑡 2 (4.8)
𝑀∗ = 𝑚 + Δ𝑡 +
2 6

Por otro lado, sea Δ𝑄 ∗ el vector que contiene al incremento de cargas.

𝑘 Δ𝑡 2 (4.9)
Δ𝑄 ∗ = Δ𝑄 − 𝑘 Δ𝑡 𝑞̇ 𝑖 − 𝑞̈ 𝑖 ( + 𝑐 Δ𝑡)
2

De donde el incremento de aceleración se halla dividiendo el incremento de carga para

la masa equivalente del sistema.

Δ𝑄 ∗ (4.10)
Δ𝑞̈ =
𝑀∗
 Procedimiento de cálculo

i. Se determina la masa equivalente del sistema 𝑀∗

𝑐 𝑘 Δ𝑡 2
𝑀∗ = 𝑚 + Δ𝑡 +
2 6

ii. Se halla el incremento de carga Δ𝑄 ∗

𝑘 Δ𝑡 2
Δ𝑄 ∗ = Δ𝑄 − 𝑘 Δ𝑡 𝑞̇ 𝑖 − 𝑞̈ 𝑖 ( + 𝑐 Δ𝑡)
2

Δ𝑄 = 𝑄𝑖+1 − 𝑄𝑖
.
iii. Se calcula el incremento de aceleraciones Δ𝑞̈

Δ𝑄 ∗
Δ𝑞̈ =
𝑀∗

iv. Se encuentra el incremento de velocidad Δ𝑞̇

Δ𝑡
Δ𝑞̇ = 𝑞̈ 𝑖 Δ𝑡 + Δ𝑞̈
2

v. Se determina el incremento de desplazamiento

Δ𝑡 2 Δ𝑞̈ Δ𝑡 2
Δ𝑞 = 𝑞̇ 𝑖 Δ𝑡 + 𝑞̈ 𝑖 +
2 6

vi. Se obtiene el nuevo desplazamiento, velocidad y aceleración en 𝑡𝑖+1

𝑞𝑖+1 = 𝑞𝑖 + Δ𝑞 𝑞̇ 𝑖+1 = 𝑞̇ 𝑖 + Δ𝑞̇ 𝑞̈ 𝑖+1 = 𝑞̈ 𝑖 + Δ𝑞̈

vii. Los valores obtenidos en el tiempo 𝑡𝑖+1 se asignan a 𝑡𝑖 para el siguiente ciclo

𝑞𝑖 = 𝑞𝑖+1 𝑞̇ 𝑖 = 𝑞̇ 𝑖+1 𝑞̈ 𝑖 = 𝑞̈ 𝑖+1

 ROBERTO AGUIAR FALCONI 113

Dentro de un lazo, este último paso es automático. Para un nuevo incremento de tiempo
se repite desde el paso dos. Es importante destacar que en el Análisis Lineal, la masa equivalente
𝑀∗ se determina una sola vez.

 EJEMPLO 1

La estructura de un piso, indicada en la figura 4.2, ha sido modelada como un sistema

de un grado de libertad; tiene una masa 𝑚 = 10 kg s2/cm; una rigidez 𝑘 = 1000 kg/cm; un factor
de amortiguamiento 𝜉 = 0.1. Está sometida al pulso rectangular de 2500 kg, que tiene una
duración de 0.5 seg. Con los datos indicados el amortiguamiento 𝑐 = 20 kg*s/cm. Si las
condiciones iniciales de desplazamiento y velocidad son nulas se pide:

1. Aplicando el Método de aceleración lineal con un Δ𝑡 = 0.05 𝑠., presentar en una tabla,
los cálculos hasta un tiempo de 0.15 seg.
2. Calcular en forma analítica empleando la respuesta 𝑔(𝑡) estudiada en el capítulo
anterior y comparar los resultados con el método de aceleración lineal para Δ𝑡 = 0.05s
3. Comparar las 2 respuestas analítica y con el método de aceleración lineal para un Δ𝑡 =
0.01 𝑠.

Figura 4.2 Sistema de un grado de libertad sometida a un pulso rectangular.

En este ejercicio 𝑘𝑔 = 𝑘𝑔𝑓 (Es kilogramo fuerza)

 SOLUCIÓN

Para 𝑡 = 0, se tiene que el desplazamiento y velocidad son cero pero no la aceleración

ya que hay una fuerza de 2500 kg. A continuación se presenta el cálculo completo para 𝑡 = 0
𝑐𝑚
𝑚 𝑞̈ + 𝑐 𝑞̇ + 𝑘 𝑞 = 𝑄 → 10 𝑞̈ = 2500 → 𝑞̈ = 250
𝑠2

𝑐 𝑘 Δ𝑡 2 20 1000 ∗ 0.052 𝑘𝑔 𝑠 2
𝑀∗ = 𝑚 + Δ𝑡 + = 10 + 0.05 + = 10.9167
2 6 2 6 𝑐𝑚

𝑘 Δ𝑡 2 1000 ∗ 0.052
Δ𝑄 ∗ = Δ𝑄 − 𝑘 Δ𝑡 𝑞̇ 𝑖 − 𝑞̈ 𝑖 ( + 𝑐 Δ𝑡) = 0 − 0 − 250 ( + 20 ∗ 0.05) = −562.50 𝑘𝑔
2 2

Δ𝑄∗ 562.50 𝑐𝑚
Δ𝑞̈ = =− = −51.5267 2
𝑀∗ 10.9167 𝑠
Δ𝑡 0.05 𝑐𝑚
Δ𝑞̇ = 𝑞̈ 𝑖 Δ𝑡 + Δ𝑞̈ = 250 ∗ 0.05 − 51.5267 ∗ = 11.2118
2 2 𝑠

Δ𝑡 2 Δ𝑞̈ Δ𝑡 2 0.052 (−51.5267 ∗ 0.052 )

Δ𝑞 = 𝑞̇ 𝑖 Δ𝑡 + 𝑞̈ 𝑖 + = 0 + 250 + = 0.2910 𝑐𝑚
2 6 2 6
114 ANÁLISIS DE SISTEMAS DE 1 GDL CON MÉTODOS NUMÉRICOS

1. Tabla de cálculo

En la tabla 4.1 se presenta el cálculo hasta 0.15 s., y tiene por objetivo que el lector
compruebe los valores presentados.

Tabla 4.1 Resultados hasta 0.15 s.

𝑡 Δ𝑄∗ Δ𝑞̈ 𝑞̈ Δ𝑞̇ 𝑞̇ Δ𝑞 𝑞
(s.) (kg.) (cm/s2) (cm/s2) (cm/s) (cm/s) (cm) (cm)
0.00 -562.50 -51.5267 250.00 11.2118 0 0.2910 0
0.05 -1007.15 -92.2586 198.4733 7.6172 11.2118 0.7702 0.2910
0.10 -1180.43 -108.1314 106.2147 2.6074 18.8290 1.0292 1.0613
0.15 -1067.51 -97.7873 -1.9167 -2.5405 21.4365 1.0287 2.0904

2. Respuesta analítica

Para un pulso rectangular la respuesta en el tiempo, es la siguiente:

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

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

Para el presente ejemplo 𝑇 = 0.5 𝑠.; 𝐹𝑜 = 2500 𝑘𝑔. Para el caso sub amortiguado, la
respuesta 𝑔(𝑡) vale:
1 𝜉
𝑔(𝑡) = [1 − exp(−𝜉 𝑊𝑛 𝑡) (cos 𝑊𝑎 𝑡 + 𝑠𝑒𝑛 𝑊𝑎 𝑡)]
𝑘 √1 − 𝜉 2

Reemplazando valores se tiene:

𝑘 1000 1
𝑊𝑛 = √ = √ = 10
𝑚 10 𝑠

𝑊𝑎 = 𝑊𝑛 √1 − 𝜉 2 = 9,9499 1/s.

𝜉
= 0.1005
√1 − 𝜉 2

Para 𝑡 < 0.5 𝑠

2500
𝑞(𝑡) = [1 − exp(−𝑡) (cos 9.9499𝑡 + 0.1005 ∗ 𝑠𝑒𝑛 9.9499𝑡)]
1000

Para 𝑡 ≥ 0.5 𝑠

𝑞(𝑡) = 2.5 [1 − exp(−𝑡) (cos 9.9499𝑡 + 0.1005 ∗ 𝑠𝑒𝑛 9.9499𝑡)]

−2.5 [1 − exp[−(𝑡 − 0.5)] [cos 9.9499(𝑡 − 0.5) + 0.1005 ∗ 𝑠𝑒𝑛 9.9499(𝑡 − 0.5)]]

En la figura 4.3 se presenta la respuesta en forma analítica y la hallada con el método

de la aceleración lineal para un Δ𝑡 = 0.05 𝑠. Se aprecia que la aproximación que se halla con el
método de aceleración lineal es bastante buena.
 ROBERTO AGUIAR FALCONI 115

Figura 4.3 Comparación de respuestas: Solución analítica y aceleración lineal con Δ𝑡 = 0.05 𝑠.

3. Comparación de Respuestas: Analítica y con Aceleración Lineal para Δ𝑡 = 0.01 𝑠.

Cuando el Δ𝑡 tiende a cero la respuesta hallada con el método de aceleración lineal

(métodos numéricos) se aproxima muy bien a la respuesta analítica exacta y esto se debe a que
se tiene una mayor cantidad de puntos discretos en los cuales se cumple la ecuación diferencial.

En La figura 4.4 se presenta la respuesta en el tiempo hallada con la solución analítica

exacta y la que se encuentra con el método de la aceleración lineal para Δ𝑡 = 0.01 𝑠.
Prácticamente coinciden las dos curvas.

Figura 4.4 Comparación de respuestas: Solución analítica y aceleración lineal con Δ𝑡 = 0.01 𝑠.

El archivo Ej1_Newmark se encuentra el ejemplo desarrollado.

116 ANÁLISIS DE SISTEMAS DE 1 GDL CON MÉTODOS NUMÉRICOS

4.2 MÉTODO DE ACELERACIÓN CONSTANTE

El método de aceleración constante, considera que la respuesta de aceleraciones, en el

intervalo de tiempo 𝑡𝑖 y 𝑡𝑖+1 es constante, como se presenta en la figura 4.5; en este caso la
ecuación de la aceleración para un tiempo 𝑡𝑖 ≤ 𝑡 ≤ 𝑡𝑖+1 es:

Δ𝑞̈ (4.11)
𝑞̈ (𝑡) = 𝑞̈ 𝑖 +
2

Figura 4.5 Modelo de aceleración constante para la respuesta en el tiempo.

Para encontrar la expresión que define la velocidad en el intervalo, se debe derivar la

ecuación de aceleraciones.
𝑡 𝑡
Δ𝑞̈
∫ 𝑞̈ (𝑡) 𝑑𝑡 = ∫ (𝑞̈ 𝑖 + ) 𝑑𝑡
2 (4.12)
𝑡𝑖 𝑡𝑖

Δ𝑞̈
𝑞̇ (𝑡) − 𝑞̇ 𝑖 = 𝑞̈ 𝑖 (𝑡 − 𝑡𝑖 ) + (𝑡 − 𝑡𝑖 )
2

Al encontrar la velocidad en 𝑡𝑖+1 y recordando que Δ𝑞̇ = 𝑞̇ 𝑖+1 − 𝑞̇ 𝑖 , y la ecuación de Δ𝑡

se tiene:

Δ𝑞̈ (4.13)
Δ𝑞̇ = 𝑞̈ 𝑖 Δ𝑡 + Δ𝑡
2

Ahora al integrar la ecuación (4.12) se encuentra la ecuación que define el

desplazamiento pero antes escribamos de nuevo esa ecuación pasando el término 𝑞̇ 𝑖 al lado
derecho.
𝑡 𝑡
Δ𝑞̈
∫ 𝑞̇ 𝑑𝑡 = ∫ (𝑞̇ 𝑖 + 𝑞̈ 𝑖 (𝑡 − 𝑡𝑖 ) + (𝑡 − 𝑡𝑖 ) ) 𝑑𝑡
2
𝑡𝑖 𝑡𝑖

(𝑡 − 𝑡𝑖 )2 Δ𝑞̈ (𝑡 − 𝑡𝑖 )2
𝑞(𝑡) − 𝑞𝑖 = 𝑞̇ 𝑖 (𝑡 − 𝑡𝑖 ) + 𝑞̈ 𝑖 +
2 2 2

De donde:

Δ𝑡 2 Δ𝑡 2 (4.14)
Δ𝑞 = 𝑞̇ 𝑖 Δ𝑡 + 𝑞̈ 𝑖 + Δ𝑞̈
2 4
 ROBERTO AGUIAR FALCONI 117

Al reemplazar (4.13) y (4.14) en la ecuación diferencial escrita en forma incremental,

definida en (4.2) se obtiene:

𝑐 𝑘 Δ𝑡 2 𝑘 Δ𝑡 2 (4.15)
(𝑚 + Δ𝑡 + ) Δ𝑞̈ = Δ𝑄 − 𝑘 𝑞̇ 𝑖 Δ𝑡 − (𝑐 Δ𝑡 + ) 𝑞̈ 𝑖
2 4 2

De tal manera, que en el Método de Aceleración Constante, cambian las siguientes

ecuaciones

𝑐 𝑘
𝑀∗ = 𝑚 + Δ𝑡 + Δ𝑡 2 (4.16)
2 4
Δ𝑡 2 Δ𝑡 2
Δ𝑞 = 𝑞̇ 𝑖 Δ𝑡 + 𝑞̈ 𝑖 + Δ𝑞̈
2 4

4.3 MÉTODO DE ACELERACIÓN ESCALONADA

En la figura 4.6 se presenta la variación de aceleración para el método de aceleración
escalonada. Nótese que el punto medio entre el intervalo 𝑡𝑖 y 𝑡𝑖+1 se ha denominado 𝑡̅, de tal
manera que Δ𝑡/2 = 𝑡̅ − 𝑡𝑖 = 𝑡𝑖+1 − 𝑡̅. Básicamente se tienen dos ecuaciones para la variación de
la aceleración y son:

 Para 𝑡𝑖 ≤ 𝑡 ≤ 𝑡̅
𝑞̈ (𝑡) = 𝑞̈ 𝑖

 Para 𝑡̅ ≤ 𝑡 ≤ 𝑡𝑖+1
𝑞̈ (𝑡) = 𝑞̈ 𝑖+1

Figura 4.6 Variación escalonada de la aceleración.

Por lo tanto se debe trabajar en los dos intervalos de aceleración, y en el primer intervalo
se debe determinar la velocidad y desplazamientos en 𝑡̅ que es igual a la mitad del intervalo. A
continuación se presenta el desarrollo numérico para el primer intervalo, como se lo ha realizado
en los otros apartados.

 Para 𝑡𝑖 ≤ 𝑡 ≤ 𝑡̅
𝑞̈ (𝑡) = 𝑞̈ 𝑖
118 ANÁLISIS DE SISTEMAS DE 1 GDL CON MÉTODOS NUMÉRICOS

𝑡 𝑡

∫ 𝑞̈ (𝑡) 𝑑𝑡 = ∫ 𝑞̈ 𝑖 𝑑𝑡
𝑡𝑖 𝑡𝑖

𝑞̇ (𝑡) − 𝑞̇ 𝑖 = 𝑞̈ 𝑖 (𝑡 − 𝑡𝑖 ) → 𝑞̇ (𝑡) = 𝑞̇ 𝑖 + 𝑞̈ 𝑖 (𝑡 − 𝑡𝑖 ) (4.17)

Para 𝑡 = 𝑡̅

Δ𝑡
𝑞̇ 𝑡̅ = 𝑞̇ 𝑖 + 𝑞̈ 𝑖
2

Ahora se debe integrar la ecuación (4.17) para encontrar el desplazamiento.

𝑡 𝑡

∫ 𝑞̇ (𝑡) 𝑑𝑡 = ∫(𝑞̇ 𝑖 + 𝑞̈ 𝑖 (𝑡 − 𝑡𝑖 ) )𝑑𝑡

𝑡𝑖 𝑡𝑖

(𝑡 − 𝑡𝑖 )2
𝑞(𝑡) − 𝑞𝑖 = 𝑞̇ 𝑖 (𝑡 − 𝑡𝑖 ) + 𝑞̈ 𝑖
2

Para 𝑡 = 𝑡̅

Δ𝑡 Δ𝑡 2
𝑞𝑡̅ = 𝑞𝑖 + 𝑞̇ 𝑖 + 𝑞̈ 𝑖
2 8

 Para 𝑡̅ ≤ 𝑡 ≤ 𝑡𝑖+1
𝑞̈ (𝑡) = 𝑞̈ 𝑖+1
𝑡 𝑡

∫ 𝑞̈ (𝑡) 𝑑𝑡 = ∫ 𝑞̈ 𝑖+1 𝑑𝑡
𝑡̅ 𝑡̅

El límite inferior de la integral es 𝑡̅

𝑞̇ (𝑡) − 𝑞̇ 𝑡̅ = 𝑞̈ 𝑖+1 (𝑡 − 𝑡̅) → 𝑞̇ (𝑡) = 𝑞̇ 𝑡̅ + 𝑞̈ 𝑖+1 (𝑡 − 𝑡̅)

Al reemplazar la velocidad en el tiempo 𝑡̅ se encuentra

Δ𝑡
𝑞̇ (𝑡) = 𝑞̇ 𝑖 + 𝑞̈ 𝑖 + 𝑞̈ 𝑖+1 (𝑡 − 𝑡̅)
2

Ahora al evaluar la velocidad en 𝑡 = 𝑡𝑖+1 , se tiene:

Δ𝑡 Δ𝑡
𝑞̇ 𝑖+1 = 𝑞̇ 𝑖 + 𝑞̈ 𝑖 + 𝑞̈ 𝑖+1
2 2

De donde:
Δ𝑡
Δ𝑞̇ = (𝑞̈ 𝑖 + 𝑞̈ 𝑖+1 )
2

Al reemplazar 𝑞̈ 𝑖+1 = 𝑞̈ 𝑖 + Δ𝑞̈ , en esta última ecuación, se halla:

Δ𝑡 Δ𝑡
Δ𝑞̇ = (Δ𝑞̈ + 2 𝑞̈ 𝑖 ) = 𝑞̈ 𝑖 Δ𝑡 + Δ𝑞̈ (4.18)
2 2
 ROBERTO AGUIAR FALCONI 119

Ahora al integrar la ecuación de la velocidad, se halla:

𝑞̇ (𝑡) = 𝑞̇ 𝑡̅ + 𝑞̈ 𝑖+1 (𝑡 − 𝑡̅)

𝑡 𝑡

∫ 𝑞̇ (𝑡) 𝑑𝑡 = ∫(𝑞̇ 𝑡̅ + 𝑞̈ 𝑖+1 (𝑡 − 𝑡̅)) 𝑑𝑡

𝑡̅ 𝑡̅

(𝑡 − 𝑡̅)2
𝑞(𝑡) − 𝑞𝑡̅ = 𝑞̇ 𝑡̅ (𝑡 − 𝑡̅) + 𝑞̈ 𝑖+1
2

Para 𝑡 = 𝑡𝑖+1

Δ𝑡 Δ𝑡 2
q(t) = 𝑞𝑡̅ + 𝑞̇ 𝑡̅ + 𝑞̈ 𝑖+1
2 8

Al reemplazar el desplazamiento 𝑞𝑡̅ y la velocidad 𝑞̇ 𝑡̅ y finalmente al evaluar el

desplazamiento en el tiempo 𝑖 + 1, se tiene:

Δt Δt 2 Δt Δt Δt 2
q𝒊+𝟏 = qi + q̇ i + q̈ i + [q̇ i + q̈ i ] + q̈ i+1
2 8 2 2 8

De donde, el incremento de desplazamientos es:

3 2 Δ𝑡 2
Δ𝑞 = 𝑞̇ 𝑖 Δ𝑡 + 𝑞̈ 𝑖 Δ𝑡 + 𝑞̈ 𝑖+1
8 8

De igual manera se deja que el lector deduzca que:

Δ𝑡 2 Δ𝑡 2 Δ𝑡 2 (4.19)
Δ𝑞 = 𝑞̇ 𝑖 Δ𝑡 + (Δ𝑞̈ + 4 𝑞̈ 𝑖 ) = 𝑞̇ 𝑖 Δ𝑡 + 𝑞̈ 𝑖 + Δ𝑞̈
8 2 8

Al reemplazar el incremento de velocidad y el incremento de desplazamiento en la

ecuación (4.2) se tiene:

Δ𝑡 Δ𝑡 2
𝑚 Δ𝑞̈ + 𝑐 (Δ𝑞̈ + 2 𝑞̈ 𝑖 ) + 𝑘 [𝑞̇ 𝑖 Δ𝑡 + (Δ𝑞̈ + 4 𝑞̈ 𝑖 )] = Δ𝑄
2 8

De donde se encuentra:

𝑀∗ Δ𝑞̈ = Δ𝑄 ∗
Siendo:
𝑐 Δ𝑡 𝑘 Δ𝑡 2 (4.20)
𝑀∗ = 𝑚 + +
2 8

𝑘 Δ𝑡 2
Δ𝑄 ∗ = Δ𝑄 − 𝑘 𝑞̇ 𝑖 Δ𝑡 − (𝑐 Δ𝑡 + ) 𝑞̈ 𝑖
2

4.4 MÉTODO 𝜷 DE NEWMARK

Se han deducido por separado, los tres casos que contempla el Método 𝛽 de Newmark.
Ahora se presentan estos tres en la forma que lo hizo el Profesor Nathan Newmark. Pero solo se
indican las ecuaciones que varían en los tres casos y son:
120 ANÁLISIS DE SISTEMAS DE 1 GDL CON MÉTODOS NUMÉRICOS

Δ𝑡 2 (4.21)
Δ𝑞 = 𝑞̇ 𝑖 Δ𝑡 + 𝑞̈ 𝑖 + 𝛽 Δ𝑞̈ Δ𝑡 2
2

𝑐 Δ𝑡
𝑀∗ = 𝑚 + + 𝛽 𝑘 Δ𝑡 2 (4.22)
2

Para el método de aceleración lineal 𝛽 =1/6 (0.1667); para el método de aceleración

constante 𝛽 = 1/4 (0.25); y para el método de aceleración escalonada 𝛽 = 1/8 (0.125). De tal
manera que se puede indicar que 𝛽 varía entre:

0.125 < 𝛽 < 0.25

Esto para el caso en que no se introduce amortiguamiento al método, tema que estudiará
en el siguiente apartado. En principio se puede indicar que los valores de 𝛽 se hallan entre 0.125
y 0.25; de tal manera que si un ejemplo es inestable con un determinado valor de 𝛽 se puede
trabajar con otro valor en el intervalo dado y si continúa siendo inestable considerar lo que se
indica en el siguiente apartado.

 EJEMPLO 2

Encontrar la respuesta en desplazamientos de la estructura del Ejemplo 1, ante el pulso

rectangular de 2500 kg., que tiene una duración de 0.5 s. Trabajar con Δ𝑡 = 0.05 𝑠.

Presentar en un solo gráfico las respuestas considerando el valor de 𝛽 para: aceleración

lineal; aceleración constante y aceleración escalonada.

 SOLUCIÓN

𝑚 = 10 kg s2/cm; 𝑘 = 1000 kg/cm; 𝑐 = 20 kg/cm.

En la figura 4.7 se presentan las tres respuestas que se obtienen con los diferentes
valores de 𝛽. La respuesta para el método de aceleración lineal y escalonada, prácticamente es
la misma; para aceleración constante existe una pequeña diferencia. El archivo Ej2_Newmark
contiene el Ejemplo desarrollado.

Figura 4.7 Respuesta ante un Pulso Rectangular

 ROBERTO AGUIAR FALCONI 121

 EJEMPLO 3

La estructura del Ejemplo 1, que tiene: 𝑚 = 10 kg s2/cm; 𝑘 = 1000 kg/cm; 𝑐 = 20 kg/cm;

𝜉 = 0.10, está sometida a la carga triangular indicada en la figura 4.8, se pide:

1. Encontrar la respuesta en el tiempo, empleando el Método 𝛽 de Newmark, utilizando Δ𝑡 =

0.01 𝑠. Presentar las tres primeras filas de cálculo para el método de aceleración lineal.
2. En un gráfico presentar la respuesta en el tiempo, considerando el valor de 𝛽 para:
aceleración lineal; aceleración constante y aceleración escalonada. Considerar Δ𝑡 =
0.05 𝑠.
3. Utilizando la Integral de Duhamel presentar la respuesta en el tiempo y comparar los
resultados con el Método de aceleración lineal, en una gráfica. Considerar Δ𝑡 = 0.01 𝑠.

Figura 4.8 Sistema de un grado de libertad sometido a una fuerza triangular

 SOLUCIÓN

1) Método de aceleración lineal

En la tabla 4.2 se presenta la respuesta en el tiempo para las tres primeras filas de
cálculo. El objetivo es que el lector compruebe estos valores.

Tabla 4.2 Resultados hasta 0.15 s., empleando Método de Aceleración Lineal con Δ𝑡 = 0.01 𝑠.
𝑡 Δ𝑄 ∗ Δ𝑞̈ 𝑞̈ Δ𝑞̇ 𝑞̇ Δ𝑞 𝑞
(s.) (kg.) (cm/s2) (cm/s2) (cm/s) (cm/s) (cm) (cm)
0.00 -112.50 -11.1248 250.00 2.444 0 0.0124 0
0.01 -134.1625 -13.2670 238.8752 2.3224 2.4444 0.0362 0.0124
0.02 -154.07 -15.2356 225.6082 2.1799 4.7668 0.0588 0.0486

2) Respuesta en el tiempo con Δ𝑡 = 0.05 𝑠.

Al encontrar la respuesta con Δ𝑡 = 0.01 𝑠., los tres métodos reportan los mismos
resultados por este motivo se cálculo con Δ𝑡 = 0.05 𝑠., para ver alguna diferencia; los resultados
se presentan en la figura 4.9.
122 ANÁLISIS DE SISTEMAS DE 1 GDL CON MÉTODOS NUMÉRICOS

Figura 4.9 Respuesta en el tiempo ante fuerza triangular. Δ𝑡 = 0.05 𝑠.

3) Integral de Duhamel

Aplicando las ecuaciones descritas en el Capítulo 3, al caso de que la excitación sea una
fuerza en lugar de una aceleración, se tiene:

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

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

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

Para la fuerza triangular su ecuación es:

𝜏
𝐹(𝜏) = 𝐹0 (1 − )
𝑡𝑑

Donde 𝐹0 = 2500 𝑘𝑔; 𝑡𝑑 = 0.5 𝑠. 𝑊𝑛 = √𝑚𝑘 = √1000

10
1
= 10 ; 𝑊𝑎 = 𝑊𝑛 √1 − 𝜉 2 = 9,9499 1/s.
𝑠

En la programación primero se halla las variables 𝐴(𝑡)y 𝐵(𝑡) utilizando la función trapz
de Matlab. La integración se realiza tantas veces como incrementos de tiempo se hayan
considerado. Luego se utiliza la ecuación (4.23) para encontrar la respuesta en el tiempo, la
misma que se presenta en la figura 4.10

El archivo Ej3_Newmark resuelve el Ejemplo 3.

 ROBERTO AGUIAR FALCONI 123

Figura 4.10 Comparación de la respuesta con Integral de Duhamel y Aceleración Lineal.

 EJEMPLO 4

Resolver analíticamente el problema anterior, mediante la Integral de Duhamel para el

caso de que 𝜉 = 0, graficar la respuesta en el tiempo hallada con la Integral de Duhamel y con
el Método de Aceleración Lineal, con Δ𝑡 = 0.01 𝑠.

Figura 4.11 Sistema de un grado de libertad sin amortiguamiento ante un pulso triangular.

 SOLUCIÓN

La respuesta en el tiempo, de un sistema de un grado de libertad, sin amortiguamiento,

ante un pulso triangular, figura 4.11, aplicando la Integral de Duhamel, es la siguiente. Paz
(1992).

Para 𝑡 ≤ 𝑡𝑑
𝐹0 𝐹0 𝑠𝑒𝑛 𝑊𝑛 𝑡
𝑞(𝑡) = (1 − cos 𝑊𝑛 𝑡) + ( − 𝑡)
𝑘 𝑘 𝑡𝑑 𝑊𝑛
(4.26)
La velocidad es:

𝐹0 𝑊𝑛 𝐹0 𝐹0
𝑞̇ (𝑡) = 𝑠𝑒𝑛 𝑊𝑛 𝑡 + cos 𝑊𝑛 𝑡 −
𝑘 𝑘 𝑡𝑑 𝑘 𝑡𝑑
124 ANÁLISIS DE SISTEMAS DE 1 GDL CON MÉTODOS NUMÉRICOS

Se necesita conocer el desplazamiento y velocidad en 𝑡 = 𝑡𝑑 ya que estos valores serán

condiciones iniciales para hallar la respuesta en el tiempo en 𝑡 > 𝑡𝑑 .

𝐹𝑜 𝑠𝑒𝑛 𝑊𝑛 𝑡𝑑
𝑞𝑡𝑑 = ( − cos 𝑊𝑛 𝑡𝑑 )
𝑘 𝑊𝑛 𝑡𝑑

𝐹𝑜 𝐶𝑜𝑠 𝑊𝑛 𝑡𝑑 1
𝑞̇ 𝑡𝑑 = (𝑊𝑛 𝑠𝑒𝑛 𝑊𝑛 𝑡𝑑 + − )
𝑘 𝑡𝑑 𝑡𝑑

Para 𝑡 ≥ 𝑡𝑑 se debe resolver la siguiente ecuación diferencial.

𝑚 𝑞̈ + 𝑘 𝑞 = 0

Esto es debido a que para 𝑡 ≥ 𝑡𝑑 , la fuerza es cero. La solución de esta ecuación

diferencial homogénea con condiciones iniciales fue estudiada en el Capítulo 1 y es la siguiente:

𝐹𝑜 𝐹𝑜 (4.27)
𝑞(𝑡) = {𝑠𝑒𝑛 𝑊𝑛 𝑡 − 𝑠𝑒𝑛 𝑊𝑛 (𝑡 − 𝑡𝑑 )} − cos 𝑊𝑛 𝑡
𝑘 𝑊𝑛 𝑡𝑑 𝑘

En la figura 4.12 se presenta la respuesta de desplazamientos hallada con la Integral de

Duhamel y la que se encuentra con el Método de Aceleración Lineal, para un Δ𝑡 = 0.01 𝑠. Se nota
que éste último método se aproxima muy bien a la respuesta analítica hallada con la Integral de
Duhamel.

Figura 4.12 Respuestas encontradas con Integral de Duhamel y Aceleración Lineal.

4.5 RESPUESTA ANTE UN SISMO

La ecuación diferencial que gobierna el problema de un sistema de un grado de libertad

ante una acción sísmica, es la siguiente:

𝑚 𝑞̈ + 𝑐 𝑞̇ + 𝑘 𝑞 = −𝑚 𝑢̈ 𝑔 (4.28)
 ROBERTO AGUIAR FALCONI 125

Donde 𝑚, 𝑐, 𝑘 son la masa, amortiguamiento y rigidez; 𝑞, 𝑞̇ , 𝑞̈ son el desplazamiento,

velocidad y aceleración; 𝑢̈ 𝑔 es la aceleración del suelo. El procedimiento de cálculo al aplicar el
Método 𝛽 de Newmark, es el siguiente.

i. Se determina la masa equivalente del sistema 𝑀∗

𝑐 Δ𝑡
𝑀∗ = 𝑚 + + 𝛽 𝑘 Δ𝑡 2
2

ii. Se halla el incremento de carga Δ𝑄 ∗

𝑘 Δ𝑡 2
Δ𝑄∗ = Δ𝑄 − 𝑘 Δ𝑡 𝑞̇ 𝑖 − 𝑞̈ 𝑖 ( + 𝑐 Δ𝑡)
2

Δ𝑄 = −𝑚 (𝑢̈ 𝑖+1 − 𝑢̈ 𝑖 )

iii. Se halla el incremento de aceleraciones Δ𝑞̈

Δ𝑄 ∗
Δ𝑞̈ =
𝑀∗

iv. Se encuentra el incremento de velocidad Δ𝑞̇

Δ𝑡
Δ𝑞̇ = 𝑞̈ 𝑖 Δ𝑡 + Δ𝑞̈
2

v. Se determina el incremento de desplazamiento

Δ𝑡 2
Δ𝑞 = 𝑞̇ 𝑖 Δ𝑡 + 𝑞̈ 𝑖 + 𝛽Δ𝑞̈ Δ𝑡 2
2

vi. Se obtiene el nuevo desplazamiento, velocidad y aceleración en 𝑡𝑖+1

𝑞𝑖+1 = 𝑞𝑖 + Δ𝑞 𝑞̇ 𝑖+1 = 𝑞̇ 𝑖 + Δ𝑞̇ 𝑞̈ 𝑖+1 = 𝑞̈ 𝑖 + Δ𝑞̈

vii. Los valores obtenidos en el tiempo 𝑡𝑖+1 se asignan a 𝑡𝑖 para el siguiente ciclo

4.5.1 Programa lineal

El programa lineal, halla en forma gráfica, la respuesta en el tiempo de un sistema de

un grado de libertad, ante una acción sísmica, definida por su acelerograma, aplicando el Método
𝛽 de Newmark.

Antes de utilizarlo se debe grabar el archivo que contiene únicamente las aceleraciones
del sismo en formato ASCII. El nombre tiene extensión .dat. Después de ello cuando se
encuentra en la modalidad consola se carga el acelerograma y después se ejecuta lineal, de la
siguiente manera:

[d,v,a] = lineal (p,m,c,k,dt,beta)

 p es el nombre del archivo que contiene el acelerograma.

 m es la masa del sistema de 1 gdl.
 c es el amortiguamiento del sistema de 1 gdl.
 k es la rigidez del sistema de 1 gdl.
126 ANÁLISIS DE SISTEMAS DE 1 GDL CON MÉTODOS NUMÉRICOS

 dt es el incremento de tiempo con el cual se desea hallar la respuesta. El mismo

que tiene que ser igual al incremento de tiempo con el cual se obtuvo el
acelerograma.
 beta Valor de 𝛽 del Método, varía entre 0.125 y 0.25; el valor de 0.125 es para
aceleración escalonada; 𝛽 =0.1667 para aceleración lineal y 0.25 para
aceleración constante.

Una vez que se ejecuta lineal aparecen cuatro gráficas, la primera de ellas es el
acelerograma, que es dato. La segunda la respuesta en el tiempo de los desplazamientos, la
tercera de las velocidades y la última de las aceleraciones.

4.5.2 Sismo de Ecuador de 2016

En el Anexo A, se presentan datos sismológicos del terremoto del 16 de abril de 2016,

de magnitud de momento 7.8; en este apartado se presenta en la figura 4.13 los acelerogramas
registrados en la ciudad de Manta, en un perfil sísmico de suelo C.

(a) (b)

(c)
Figura 4.13 Registros de aceleraciones en Manta del terremoto de 2016; a) Componente E-W;
b) Componente N-S; c) Componente Vertical

El archivo rec_manta contiene los tres registros con los siguientes nombres: ew, ns, v,
para las componentes este-oeste (figura 4.13 a); norte-sur (figura 4.13 b); y vertical (figura 4.13
c). EL incremento de tiempo tiene el nombre Dt, vale 0.01 s. Los registros de aceleraciones están
en gals.

Los valores máximos en valor absoluto (PGA) de los registros se indican en la tabla 4.3.
 ROBERTO AGUIAR FALCONI 127

Tabla 4.3 Aceleraciones máximas de los registros de Manta en gals y en fracción de la

gravedad
Ciudad E-W N-S V
Manta 397.12 (0.41 g) 514.26 (0.52 g) 162.06 (0.17 g)

 EJEMPLO 5

Se desea encontrar la respuesta en el tiempo del pórtico plano indicado en la figura 4.14,
ante la componente N-S del registro de Manta (Terremoto del 16 de Abril de 2016). Los datos
que definen el sistema de un grado de libertad para un factor de amortiguamiento del 5% son:

𝑇 𝑇 𝑠2 𝑇𝑠
𝑘 = 1807.7 𝑚 = 3.06 𝑐 = 7.437
𝑚 𝑚 𝑚

𝑊𝑛 = 24.306 𝐻𝑧 𝑇 = 0.258 𝑠

Figura 4.14 Marco Plano sometido a componente N-S de registro de aceleraciones de Manta.

Presentar también las respuestas máximas de desplazamiento, velocidad y aceleración,

en valor absoluto.

 SOLUCIÓN

Se deja al lector el cálculo de los datos indicados se trabajó con un módulo de elasticidad
𝐸 = 2323800 𝑇/𝑚2 . La respuesta en el tiempo considerando 𝛽 = 1/6, se indican en la figura
4.15.

Figura 4.15 Registro de aceleraciones del movimiento del suelo y respuestas de

desplazamientos, velocidad y aceleración del sistema.
128 ANÁLISIS DE SISTEMAS DE 1 GDL CON MÉTODOS NUMÉRICOS

La gráfica superior es el registro de aceleraciones de la componente N-S en 𝑚/𝑠 2 ; luego

se tiene las respuestas de: desplazamientos (m), velocidad (m/s) y aceleración (m/s2)

Se denominan 𝑆𝑑 , 𝑆𝑣 , 𝑆𝑎 , las respuestas máximas en valor absoluto de desplazamiento,

velocidad y aceleración. Estos valores son:
𝑐𝑚 𝑐𝑚
𝑆𝑑 = 0.62 𝑐𝑚 𝑆𝑣 = 21.70 𝑆𝑎 = 822.40
𝑠 𝑠2

Las respuestas máximas interesan para encontrar los espectros de respuesta, tema que
se trata en el siguiente capítulo. Para ello se debe tener presente que la estructura analizada
tiene un período 𝑇 = 0.258 𝑠

La solución se encuentra en el programa: Ej1_sismo_Newmark

4.5.2 Sismos importantes registrados en Perú y Ecuador con registros

Perú y Ecuador se han visto afectados por sismos interplaca tipo thrust (subducción de
bajo ángulo de buzamiento) desde hace mucho tiempo, pero en la mayor parte de los casos no
se dispone de registros (acelerogramas). En la figura 4.16 se presenta 4 sismos fuertes
registrados en el Perú y son los registrados en 1966, 1970, 1974 y 2007. También se presenta
dos sismos registrados en el Ecuador, que son el del 4 de agosto de 1998, conocido como el
sismo de Bahía de Caráquez y el terremoto de Pedernales del 16 de abril de 2016.

Los sismos de 1966, 1970 y 1974, fueron utilizados en la Normativa Sísmica de 1997;
los registros que se disponen se indican en la tabla 4.4 y fueron registrados a distancias
considerables del epicentro. Los cuatro sismos, que se tienen registros, incluido el de 2007; en
Perú causaron gran daño en estructuras de adobe o mampostería no reforzada, dejando una
gran cantidad de muertos.

Figura 4.16 Sismos interplaca tipo thrust que se disponen de registros sísmicos.
Fuente: Aguiar et al. (2018)

En la figura 4.16 se indica la ubicación del epicentro y su mecanismo focal; es importante

conocer que el epicentro no es nada más que el punto donde se inicia la rotura de las placas
tectónicas. Para entender el daño que estos sismos han ocasionado se debe mirar la forma de
propagación de la rotura.

Es importante conocer que se cuenta con registros de aceleraciones, para que el diseño
sísmico de estructuras localizadas en el Ecuador, se lo haga con registros de Perú y Ecuador,
 ROBERTO AGUIAR FALCONI 129

en lugar de utilizar registros de otras regiones del mundo, que a lo mejor su tectónica no tiene
nada que ver con la de Ecuador.

Tabla 4.4 Acelerogramas cuyos registros se tienen en sistema CEINCI-LAB

Fecha Latitud Longitud Profundidad 𝑴𝒘 Δ𝒕 Archivo Unidad
Focal (km) (s)
17/10/66 -10,72 -78,70 37 8.1 0.02 rec_7035 gals
rec_7036
31/05/70 -9,36 -78,87 64 7.9 0.02 rec_7038 gals
rec_7039
3/10/74 -12,50 -77,90 21 8.0 0.02 rec_7050 gals
rec_7051
15/8/07 -13,67 -76,76 40 7.9 0.01 rec_7081 gals
ux,uy,uv
16/4/16 0o22’16’’ 79o56’24’’ 20 7.8 0.01 rec_manta gals
rec_portoviejo
rec_chone
rec_pedernales

Se debe verificar el comportamiento de una estructura no solo ante sismos interplaca de

bajo ángulo de buzamiento sino también ante sismos asociados a fallas geológicas que se
denominan corticales. Todo depende del lugar de ubicación de la estructura.

En los estudios de peligrosidad sísmica por métodos probabilísticos se determina la

desagregación donde se puede ver de donde proviene la mayor amenaza sísmica. La
desagregación determina tres factores: Magnitud, distancia al plano de ruptura y la incertidumbre.
Dependiendo del período para el cual se obtenga la desagregación se verá la fuente que es más
importante considerar. Ordaz et al. (2013)

Un aspecto fundamental, es conocer el perfil sísmico del suelo, en que se obtuvieron los
registros sísmicos. Para el caso de los registros del terremoto de 2016 de Ecuador, se debe
indicar que los registros de Manta fueron hallados en perfil sísmico de suelo C y los de Portoviejo
en perfil de suelo D. Parece ser que el registro de Chone fue hallado en perfil de suelo D o E.

REFERENCIAS.

1. Aguiar R., Rivas A., Iza E., (2018), “Determinación de espectros y sismos para el diseño
de Hospital de Quinindé con aisladores sísmicos”, Revista Internacional de Ingeniería de
Estructuras, 23 (1), 73-100.
2. Hayes, G. P.,Wald, D. J. & Johnson, R. L., (2012), “Slab1.0: a three-dimensional model
of global subduction zone geometries. J. Geophys”. Res. Solid Earth, 117, B01302.
3. Lamar S., (1981), Clases magistrales de Dinámica de Estructuras, Post Grado en
Ingeniería Sismo Resistente, Universidad Central de Venezuela, Caracas.
4. Nocquet, J., Jarrin P., Vallée M., Mothes P., Grandin R., Rolandone F., Delouis B., Yepes
H., Font Y., Fuentes D., Régnier M., Laurendeau A., Cisneros D., Hernandez S., Sladen
A., Singaucho J., Mora H., Gómez J., Montes L., Charvis P., (2016), Supercycle at the
Ecuadorian subduction zone revealed after the 2016 Pedernales earthquake. Nature
Geoscience Journal. 10, 145-149.
5. Ordaz M., Martinelli F., D’Amico V., and Maletti C., (2013), “CRISIS2008: A flexible tool
to perform probabilistic seismic hazard assessment”, Seismological Research Letters, 84
(3), 495-504.
130 ANÁLISIS DE SISTEMAS DE 1 GDL CON MÉTODOS NUMÉRICOS

ANEXO A

EL TERREMOTO DEL 16 DE ABRIL DE 2016

Durante el siglo XX la costa norte de Ecuador soportó una serie de sismos de alta
magnitud, el gran terremoto en 1906 (Mw 8.8) ocasionó la ruptura de aproximadamente 500 km
de longitud de la placa principal; en las siguientes décadas en la misma zona se rompieron
diversos sub-segmentos a causa de sismos suscitados de manera secuencial de sur a norte en:
1942, Jama (Mw 7.8); 1958, Esmeraldas (Mw 7.7); 1979 Sur de Colombia (Mw 8.2); 1998 Bahía
de Caráquez (Mw 7.2). Se añade a esta lista el terremoto del 16 de abril de 2016 (Mw 7.8). En
110 años se han registrado, en la costa norte del Ecuador, 6 sismos de magnitud mayor a 7, lo
que significa que cada 20 años (redondeando) se tiene un sismo muy fuerte en el Ecuador y para
eso se debe tener estructuras seguras.

Figura A1 Sismos interplaca tipo thrust registrados en la costa de Ecuador, en el período

1906-2016. Fuente: Hayes et al. (2016)

La línea negra de figura A1, con triángulos señala la fosa de subducción; las estrellas
negras muestran los lugares de los epicentros de los grandes terremotos del siglo XX; los círculos
blancos gruesos delimitan el área de liberación de energía de los terremotos; la línea blanca
paralela a la costa ecuatoriana muestra la extensión de ruptura generada por el gran terremoto
de 1906 (Mw 8.8); la estrella amarilla marca el epicentro del terremoto del 16 de abril de 2016 en
Pedernales.

El terremoto de Pedernales en 2016 (Mw 7.8) supone una nueva etapa de eventos
sísmicos, debido a que se encuentra en el mismo lugar del segmento ya roto en el terremoto de
1942 (Mw 7.8) estos acontecimientos fueron el resultado del empuje de fallas poco profundas en
el límite de las placas de Nazca y Sudamérica; en el lugar del terremoto, la placa de Nazca
subduce hacia el este por debajo de la placa Sudamericana a una velocidad de 47 mm / año.
Nocquet et al. (2016)
El sismo de Pedernales rompió un segmento de cerca de 100 km de longitud por 40 km
de ancho de la placa principal en profundidades entre 15 y 30 km, esto dio lugar a un
deslizamiento sísmico de norte-sur de uno a seis metros lo que ocasionó como resultado la
ruptura de dos asperezas adyacentes distintas. La segunda aspereza de 40×40 km2 se rompió
en un lapso de tiempo de 10 s con un deslizamiento superior a 5 m y generó más de la mitad de
la liberación total de la energía sísmica, esta parte de la ruptura situada a 20 km de profundidad
 ROBERTO AGUIAR FALCONI 131

de la zona poblada de la costa oeste hace relación al extenso daño sufrido de Pedernales a
Canoa. (Nocquet et al. 2016).