MÉTODOS NUMÉRICOS DE INTEGRACIÓN - ESPECTROS DE RESPUESTA ELÁSTICOS

FACULTAD DE INGENIERÍA Y CIENCIAS

ESCUELA DE INGENIERÍA CIVIL EN OBRAS CIVILES

Alumno: Raúl Sánchez Ortiz.

Profesor Cátedra: Gaspar Auad Álvarez.

Profesor Auxiliar: Williams Amaya Baldeon.

11 de octubre de 2024.

Santiago, Chile
 COC-3100
INGENIERÍA SÍSMICA

Índice
1. Introducción. ............................................................................................................................................................... 1

2. Problemas planteados. ............................................................................................................................................ 2

Problema 1. (Implementación de métodos numéricos) ....................................................................................... 2

Problema 2. .......................................................................................................................................................................... 14

Problema 3. .......................................................................................................................................................................... 21

3. Conclusiones............................................................................................................................................................. 30

Tabla de Ilustraciones.

Ilustración 1:Código de MATLAB para el método de Newmark. Fuente: Elaboración propia. ........................... 4
Ilustración 2:Código de MATLAB para el método de ZOH. Fuente: Elaboración propia........................................ 7
Ilustración 3:Código de MATLAB para la comparación de los métodos Newmark y ZOH. Fuente: Elaboración
propia. ...................................................................................................................................................................................................... 9
Ilustración 4: Resultados de comparación de ambos métodos respecto al desplazamiento, velocidad y
aceleración. Fuente: Elaboración propia. ............................................................................................................................... 11
Ilustración 5: Parte 1 del código de respuesta para pseudo aceleración y pesudo desplazamiento.
Elaboración: Fuente propia. ......................................................................................................................................................... 15
Ilustración 6:Parte 2 del código de respuesta para pseudo aceleración y pesudo desplazamiento.
Elaboración: Fuente propia. ......................................................................................................................................................... 16
Ilustración 7: Gráfica del espectro de respuesta de pseudo desplazamiento y pseudo aceleración a partir de
la condiciones del enunciado. Fuente: Elaboración propia. ............................................................................................ 18
Ilustración 8: Código para desplazamiento u(t) y ciclo del disipador. Fuente: Elaboración propia. ............. 24
Ilustración 9: Desplazamiento u(t) y ciclo del disipador. Fuente elaboración propia......................................... 25
Ilustración 10: Código para variación de energías en el sistema durante el evento sísmico. Fuente:
Elaboración propia........................................................................................................................................................................... 27
Ilustración 11: Gráfico para variación de energías en el sistema durante el evento sísmico. Fuente:
Elaboración propia........................................................................................................................................................................... 28
 COC-3100
INGENIERÍA SÍSMICA

1. Introducción.

El análisis del comportamiento dinámico de sistemas estructurales bajo cargas sísmicas es esencial en la
ingeniería sísmica. Una forma común de estudiar este fenómeno es mediante la idealización de la
estructura como un sistema de un grado de libertad, lo que permite analizar su respuesta ante
excitaciones externas, como los terremotos, en términos de desplazamiento, velocidad y aceleración. Para
llevar a cabo este análisis, se utilizan las propiedades del sistema, como su masa, rigidez y
amortiguamiento, con el fin de describir su respuesta dinámica bajo la acción de una excitación sísmica.

El objetivo principal de este trabajo es implementar y comparar dos métodos numéricos frecuentemente
empleados para resolver ecuaciones diferenciales de segundo orden que describen el movimiento de un
sistema SDOF: el método de Newmark y el método de Zero-Order Hold (ZOH). Estos métodos se aplicarán
para analizar la respuesta de un sistema cuando es sometido al registro sísmico de Concepción,
correspondiente al terremoto ocurrido el 27 de febrero de 2010. Asimismo, se busca incrementar la
experiencia en el uso del software MATLAB, herramienta ampliamente utilizada en la ingeniería para la
resolución de problemas complejos mediante métodos numéricos.

El método de Newmark es un esquema de integración numérica que se utiliza para resolver ecuaciones
de movimiento lineales y no lineales, y proporciona una solución aproximada para el desplazamiento, la
velocidad y la aceleración del sistema en cada instante de tiempo. El método ZOH, por su parte, se basa
en el cálculo de la exponencial de la matriz del sistema para obtener una aproximación de la respuesta
dinámica del sistema a lo largo del tiempo. Se analizarán y compararán los resultados obtenidos para los
desplazamientos, velocidades y aceleraciones, observando las similitudes y diferencias entre ambos
métodos. Adicionalmente, la implementación en MATLAB permitirá afianzar habilidades en el uso de este
software, crucial para la resolución de problemas numéricos en ingeniería.

1
Métodos numéricos de integración - Espectros de respuesta elásticos
 COC-3100
INGENIERÍA SÍSMICA

2. Problemas planteados.

Problema 1. (Implementación de métodos numéricos)

Problema 1(a).

Resolución.

En este problema se implementa el método de Newmark para resolver la ecuación diferencial de

segundo orden que describe el comportamiento dinámico de un sistema de un grado de libertad
(SDOF) sometido a una excitación externa. El sistema está caracterizado por su masa 𝑚,
constante de amortiguamiento 𝑐 y rigidez 𝑘, y se somete a una fuerza externa que varía en el
tiempo, representada por el vector 𝑝(𝑡). Este método de Newmark permite obtener los
desplazamientos, velocidades y aceleraciones del sistema en función del tiempo, siendo
especialmente útil en la dinámica estructural para sistemas sometidos a cargas sísmicas.

2
Métodos numéricos de integración - Espectros de respuesta elásticos
 COC-3100
INGENIERÍA SÍSMICA

En esta implementación, se ha utilizado el valor de 𝛽 =1/6, que corresponde a una integración

acelerada con mayor amortiguamiento numérico, lo cual mejora la estabilidad en algunos casos.
El parámetro 𝛾 se mantiene en 𝛾 =0.5, lo que asegura una integración centrada, balanceando
estabilidad y precisión. El bucle de integración del código permite actualizar el desplazamiento
𝑢, la velocidad 𝑢𝑝 y la aceleración 𝑢𝑝𝑝 del sistema en cada instante de tiempo. La matriz de
rigidez efectiva 𝐾𝑒𝑓𝑓 se define como:

A partir de esta matriz, el desplazamiento en el siguiente paso de tiempo se calcula resolviendo:

donde 𝑝𝑒𝑓𝑓 es la fuerza efectiva, que considera las contribuciones del desplazamiento, la
velocidad y la aceleración en el instante anterior, junto con la fuerza externa aplicada.

3
Métodos numéricos de integración - Espectros de respuesta elásticos
 COC-3100
INGENIERÍA SÍSMICA

El código desarrollado para este método en MATLAB se presenta a continuación:

Ilustración 1:Código de MATLAB para el método de Newmark. Fuente: Elaboración propia.

Observaciones importantes: En la línea "n = length(p);", se define el tamaño del vector de fuerzas
aplicadas, lo que permite determinar la cantidad de pasos temporales. Luego, se inicializan los
vectores para desplazamiento, velocidad y aceleración ("u", "up" y "upp") con ceros,
correspondientes al número de pasos de la simulación. En la asignación de condiciones iniciales,
"u(1) = CI(1);" y "up(1) = CI(2);", se establecen los valores iniciales de desplazamiento y
velocidad, mientras que la aceleración inicial se calcula en "upp(1) = (p(1) - c up(1) - k u(1)) /
m;", lo que resuelve la ecuación de equilibrio para t=0.

En las líneas "a0 = 1 / (beta dt^2);" y "a1 = gamma / (beta dt);", se definen los coeficientes del
método de Newmark, los cuales determinan cómo se ponderan los desplazamientos, velocidades
y aceleraciones en cada paso de tiempo. Los coeficientes "a2" y "a3" también se establecen,
siendo "a2 = 1 / (beta dt);" y "a3 = 1 / (2 beta) - 1;", lo que asegura que la integración sea
precisa en cada iteración.

La matriz efectiva "K_eff", calculada en "K_eff = k + a0 m + a1 c;", combina la rigidez real "k" con
términos que incluyen la masa y el amortiguamiento del sistema, siendo clave para determinar
el desplazamiento en el siguiente paso.

El bucle "for i = 1
4
Métodos numéricos de integración - Espectros de respuesta elásticos
 COC-3100
INGENIERÍA SÍSMICA

" implementa la integración temporal del sistema. Dentro del bucle, "p_eff = p(i+1) + m (a0 u(i)
+ a2 up(i) + a3 upp(i)) + c (a1 u(i));" calcula la fuerza efectiva en el instante 𝑡𝑖 + 1,
considerando la inercia, el amortiguamiento y las fuerzas externas. El nuevo desplazamiento se
resuelve con "u(i+1) = p_eff / K_eff;", utilizando la rigidez efectiva para determinar cómo se
deforma el sistema en el siguiente instante.

La aceleración y velocidad se actualizan con las líneas "upp(i+1) = a0 (u(i+1) - u(i)) - a2 up(i)
- a3 upp(i);" y "up(i+1) = up(i) + dt ((1 - gamma) upp(i) + gamma upp(i+1));",
respectivamente. Aquí, la aceleración se calcula como una función del cambio en el
desplazamiento entre los pasos actuales y anteriores, ajustada por la velocidad, mientras que la
velocidad se actualiza considerando la aceleración media en el intervalo de tiempo.

Este enfoque asegura que el desplazamiento, la velocidad y la aceleración sean consistentes en

cada paso de tiempo, proporcionando una solución numérica estable y precisa de la ecuación de
movimiento del sistema bajo fuerzas externas.

5
Métodos numéricos de integración - Espectros de respuesta elásticos
 COC-3100
INGENIERÍA SÍSMICA

El método ZOH (Zero-Order Hold) es un enfoque utilizado para simular sistemas dinámicos,
especialmente en el análisis de estructuras bajo excitaciones externas. Este método es útil para
sistemas de un grado de libertad, ya que permite calcular la respuesta del sistema en función de
las fuerzas aplicadas a lo largo del tiempo.

6
Métodos numéricos de integración - Espectros de respuesta elásticos
 COC-3100
INGENIERÍA SÍSMICA

El sistema se puede describir mediante un conjunto de ecuaciones en el espacio de estado, donde

la dinámica del sistema está representada por las siguientes ecuaciones:

Donde;

Es así como a partir de este enfoque y uso del método se determina la siguiente configuración
para la función ZOH.

Ilustración 2:Código de MATLAB para el método de ZOH. Fuente: Elaboración propia.

7
Métodos numéricos de integración - Espectros de respuesta elásticos
 COC-3100
INGENIERÍA SÍSMICA

Observaciones importantes:
- En la línea "n = length(p);", se define el tamaño del vector de fuerzas aplicadas. Este paso
es crucial ya que determina cuántos pasos temporales se realizarán en la simulación. La
longitud del vector de fuerzas se utiliza para definir el número total de iteraciones en el
bucle que calcula la respuesta del sistema.
- Luego, se inicializa la matriz de estado "X" con ceros, donde la primera fila representa el
desplazamiento y la segunda fila representa la velocidad. Esta inicialización es esencial
para almacenar los resultados de cada paso temporal de la simulación.
- En la línea "X(:, 1) = X0;", se establece la condición inicial del sistema. Aquí, se asignan los
valores iniciales de desplazamiento y velocidad que son proporcionados en el vector
X0X0X0.
- La exponencial de la matriz AAA se calcula en "expA_dt = expm(A dt);". Esta operación
es fundamental ya que permite resolver el sistema de manera efectiva en el dominio del
tiempo. La exponencial de la matriz es utilizada para obtener el estado del sistema en el
siguiente instante de tiempo a partir del estado actual.
- El bucle "for i = 1" implementa la integración temporal del sistema. Dentro del bucle, la
expresión "X(:, i+1) = expA_dt X(:, i) + inv(A) (expA_dt - I) B p(i);" actualiza el estado
del sistema en el tiempo ti+1t_{i+1}ti+1. Esta fórmula combina el efecto de la matriz
exponencial (que captura la evolución del sistema) con el estado anterior y la entrada de
fuerza aplicada en el instante ti.
- Al final del bucle, se calcula la velocidad y la aceleración con la línea "Xp = A X + B p(:)';".
Esta línea aplica la matriz de estado para obtener la velocidad y la aceleración en el
tiempo actual, utilizando la relación entre el estado y las fuerzas aplicadas.

Así de esta manera por medio del método ZOH se asegura que el desplazamiento y la velocidad
del sistema sean consistentes en cada paso de tiempo, proporcionando una solución precisa y
eficiente para la respuesta dinámica del sistema bajo la influencia de fuerzas externas. La
implementación del método permite una mejor comprensión de la dinámica estructural y facilita
el análisis de la respuesta a diferentes tipos de excitaciones.

8
Métodos numéricos de integración - Espectros de respuesta elásticos
 COC-3100
INGENIERÍA SÍSMICA

A partir de las funciones establecidas anteriormente se determina el siguiente código para

comparar ambos métodos en desplazamiento, velocidad y aceleración esta última solo para el
método Newmark. Se muestra a continuación el código obtenido:

Ilustración 3:Código de MATLAB para la comparación de los métodos Newmark y ZOH. Fuente: Elaboración propia.

9
Métodos numéricos de integración - Espectros de respuesta elásticos
 COC-3100
INGENIERÍA SÍSMICA

En la línea "ruta archivo = 'Conce_2010_Mw8.[Link]';", se especifica la ubicación del archivo, que

contiene los datos de aceleración registrados durante el terremoto. La función "load" se utiliza
para importar estos datos, lo que permite el acceso a la estructura que contiene las variables
necesarias para el análisis.

La variable "valores aceleracion" se inicializa con "valores aceleracion = [Link];", donde

se extraen los valores de aceleración del registro sísmico. Esta información es crucial ya que se
utilizará para calcular la fuerza sísmica que actúa sobre el sistema. Posteriormente, se aplanan
los datos con "valores aceleración = valores_aceleracion (:);", asegurando que el vector tenga la
forma correcta para los cálculos posteriores.

Los parámetros del sistema se definen a continuación. La masa del sistema se establece en "m =
1;", y el período natural se define como "Tn = 0.8;". Estos valores son fundamentales para
determinar las características dinámicas del sistema. La rigidez del sistema se calcula con "k =
(2 𝜋/ Tn)^2 m;", utilizando la relación entre el período natural y la rigidez, lo que permite
obtener el comportamiento elástico del sistema.

Además, se define la razón de amortiguamiento crítico como "xi = 0.02;", y se calcula el

coeficiente de amortiguamiento con "c = 2 xi sqrt(k m);". Este coeficiente es esencial para
modelar la respuesta del sistema ante excitaciones externas, ya que determina cómo se disipará
la energía en el sistema.

La fuerza sísmica se calcula en "p = -m valores aceleración;", donde se multiplica la masa del
sistema por los valores de aceleración obtenidos del registro sísmico. Este paso es crucial para
simular cómo las fuerzas inerciales afectan al sistema durante el evento sísmico.

Las condiciones iniciales del sistema se establecen a través de las variables "u0" y "up0", que
representan el desplazamiento y la velocidad iniciales, respectivamente. Ambas se inicializan en
cero, ya que el sistema comienza en reposo. El vector de condiciones iniciales se forma con "X0
= [u0; up0];", que se utilizará en las funciones de integración para calcular la respuesta del
sistema.

Los parámetros del método de Newmark se definen con "beta = 1/6;" y "gamma = 1/2;". El valor
de 𝛽vse establece en 1/6 para asegurar que el método de integración sea estable y preciso en la
simulación, mientras que 𝛾 se utiliza para ponderar la aceleración en el cálculo.

A continuación, se llama a las funciones que implementan los métodos de integración.

"Newmark_1gdl" se utiliza para calcular el desplazamiento, la velocidad y la aceleración del
sistema bajo el método de Newmark, utilizando todos los parámetros y condiciones iniciales
definidos previamente. Esta función realiza los cálculos necesarios para determinar la respuesta
dinámica del sistema ante la excitación sísmica.
10
Métodos numéricos de integración - Espectros de respuesta elásticos
 COC-3100
INGENIERÍA SÍSMICA

Simultáneamente, se llama a "ZOH_1gdl" para calcular el desplazamiento y la velocidad mediante

el método ZOH. Sin embargo, este método no calculará la aceleración, ya que su enfoque se centra
en la solución de la ecuación de movimiento sin considerar el efecto de la aceleración en cada
paso.

Finalmente, se genera un vector de tiempo con "t = (0: (p)-1) dt;", que abarca todos los pasos
temporales de la simulación. Este vector es fundamental para la visualización de los resultados,
permitiendo graficar la respuesta del sistema en función del tiempo.

Los gráficos se generan en tres subgráficas, donde se comparan los resultados de

desplazamiento, velocidad y aceleración. En la primera subgráfica, se muestra el desplazamiento
del sistema bajo ambos métodos. En la segunda subgráfica, se comparan las velocidades
calculadas por ambos métodos, y en la tercera subgráfica, se presenta únicamente la aceleración
calculada por el método de Newmark. Este enfoque de visualización permite analizar y comparar
fácilmente las respuestas del sistema bajo diferentes métodos de análisis, proporcionando
información valiosa sobre su comportamiento dinámico frente a la excitación sísmica.

Análisis de la gráfica obtenida:

Ilustración 4: Resultados de comparación de ambos métodos respecto al desplazamiento, velocidad y aceleración.

Fuente: Elaboración propia.

11
Métodos numéricos de integración - Espectros de respuesta elásticos
 COC-3100
INGENIERÍA SÍSMICA

Método de Newmark vs Método ZOH:

- Newmark (naranja): Este método es un esquema de integración numérica implícito que se
utiliza comúnmente en la dinámica estructural debido a su estabilidad incondicional para ciertos
valores de los parámetros. Para este análisis, se ha utilizado beta = 1/6, lo que corresponde al
método Newmark no amortiguado, conocido también como el método promedio constante de
aceleración. Este valor de beta garantiza estabilidad numérica para sistemas amortiguados y no
amortiguados.
- Zero-Order Hold (ZOH, azul): Este método asume que la entrada (en este caso la fuerza
sísmica) es constante durante cada intervalo de tiempo. A diferencia del método de Newmark,
ZOH es un método explícito que tiene algunas limitaciones de estabilidad dependiendo del paso
de tiempo y la rigidez del sistema. En sistemas fuertemente amortiguados o en caso de pasos de
tiempo grandes, el método puede ser más inestable en comparación con Newmark.

Análisis de la Discretización Temporal:

- Ambos métodos dependen de la discretización del tiempo para realizar la integración
numérica. El paso de tiempo (dt) es crítico para la precisión y estabilidad de ambos algoritmos.
- En la gráfica, ambos métodos se han implementado con el mismo (dt), lo que permite una
comparación justa. Sin embargo, el método Newmark generalmente es menos sensible a
variaciones en (dt), mientras que ZOH puede volverse inestable si el paso de tiempo es
demasiado grande, especialmente en sistemas con alta rigidez o bajo amortiguamiento.

Comparación del Comportamiento en el Tiempo:

- Desplazamiento y Velocidad:
- Las oscilaciones iniciales son muy pronunciadas. Esto es consistente con la baja razón de
amortiguamiento (xi = 0.02), lo que permite que las oscilaciones persistan por más tiempo.
- La buena coincidencia entre los métodos en desplazamiento y velocidad sugiere que el
sistema no es particularmente rígido, ya que, en sistemas con rigidez elevada o amortiguamiento
nulo, las diferencias entre los métodos podrían ser más notorias debido a que ZOH tiende a
subestimar o sobrestimar los máximos de respuesta.
- Aceleración:
- La aceleración calculada mediante Newmark es mucho más sensible a las pequeñas
diferencias numéricas entre los métodos. Aquí observamos una aceleración que alcanza
máximos altos, pero estabiliza rápidamente, lo que es coherente con un sistema que disipa
energía gradualmente gracias a la amortiguación.
- Sería interesante comparar la aceleración obtenida mediante ZOH, ya que las diferencias
podrían ser más evidentes en la aceleración que en el desplazamiento o velocidad, debido a que
ZOH suele simplificar la variación de las fuerzas aplicadas entre pasos de tiempo.

Estabilidad Numérica:
- El sistema tiene una baja razón de amortiguamiento crítico xi = 0.02, lo que significa que el
amortiguamiento es muy pequeño, permitiendo oscilaciones más prolongadas. En un sistema
con bajo amortiguamiento, el método de Newmark muestra una ventaja sobre ZOH, ya que tiende
a ser más estable en la simulación de sistemas con oscilaciones prolongadas, mientras que ZOH
puede sufrir de inestabilidad en estas condiciones si el paso de tiempo no es suficientemente
12
Métodos numéricos de integración - Espectros de respuesta elásticos
 COC-3100
INGENIERÍA SÍSMICA

pequeño.

Rigidez y Condiciones de Estabilidad:

- El método Newmark, siendo implícito, es incondicionalmente estable para ciertos valores de
beta y gamma. Para este caso, se utiliza beta = 1/6 y gamma = 0.5, lo que garantiza una buena
estabilidad en condiciones de baja amortiguación.
- ZOH, por su parte, depende del tamaño del paso de tiempo en relación con la frecuencia
natural del sistema. Si el paso de tiempo es demasiado grande en comparación con el período
natural, ZOH puede introducir errores de fase y amplitud, lo que en algunos casos lleva a una
respuesta artificialmente amplificada o suprimida.

Interpretación Física del Sistema:

- El sistema modelado, con Tn = 0.8 s y xi = 0.02, es representativo de estructuras con un
período relativamente corto y bajo amortiguamiento, lo que es tiempo máximo en estructuras
de baja altura o componentes no estructurales en edificaciones. Estas estructuras pueden ser
muy sensibles a excitaciones sísmicas.
- Las oscilaciones visibles durante los primeros segundos indican que el sistema está
respondiendo fuertemente a la entrada sísmica, pero la amortiguación, aunque baja, permite
disipar energía gradualmente, lo que se refleja en la disminución de las oscilaciones después de
los 50 segundos.

En resumen, la elección entre los dos métodos dependerá del nivel de precisión requerido y de
las características del sistema analizado, pero para este caso particular, ambos parecen ser
adecuados en términos de respuesta de desplazamiento y velocidad, con una ligera ventaja para
Newmark en la estabilidad general del análisis.

13
Métodos numéricos de integración - Espectros de respuesta elásticos
 COC-3100
INGENIERÍA SÍSMICA

Problema 2. (Espectro de respuesta)

Para calcular el Espectro de Respuesta, se utilizan las siguientes ecuaciones fundamentales en el análisis
dinámico de estructuras:

Es así como a partir del entendimiento de la formulas implicadas se obtiene el siguiente código:

14
Métodos numéricos de integración - Espectros de respuesta elásticos
 COC-3100
INGENIERÍA SÍSMICA

Ilustración 5: Parte 1 del código de respuesta para pseudo aceleración y pesudo desplazamiento. Elaboración: Fuente
propia.

15
Métodos numéricos de integración - Espectros de respuesta elásticos
 COC-3100
INGENIERÍA SÍSMICA

Ilustración 6:Parte 2 del código de respuesta para pseudo aceleración y pesudo desplazamiento. Elaboración: Fuente
propia.

Observaciones importantes: se comienza con la carga del archivo de datos que contiene el registro
sísmico mediante la línea "load('Conce_2010_Mw8.[Link]');". Esta acción permite acceder a la variable
que contiene la información de aceleración necesaria para el análisis. A continuación, se extraen los
valores de aceleración con "p = -[Link];", donde 𝑝 representa la fuerza inercial aplicada, y se asegura
que los datos sean un vector unidimensional.

Se establecen los parámetros del sistema: la masa se define como "m = 1000;", el período natural como
"T = periodos(i);", y la rigidez se calcula con "k = m * w^2;". La razón de amortiguamiento se define en
"xi = xi_values(j);", y el coeficiente de amortiguamiento se calcula con "c = 2 * xi * sqrt(m * k);", lo cual
es crucial para determinar cómo se comportará el sistema bajo la excitación sísmica.

El vector de tiempo se construye con "tiempo = (0)' * dt;", que abarca todos los pasos temporales de la
simulación y es fundamental para la visualización de los resultados. Para calcular el espectro de
respuesta, se inicializan las matrices para almacenar los resultados de pseudo-desplazamiento y pseudo-
aceleración: "pseudo_desplazamiento = zeros(length(xi_values), n_periodos);" y "pseudo_aceleracion =
zeros(length(xi_values), n_periodos);".

16
Métodos numéricos de integración - Espectros de respuesta elásticos
 COC-3100
INGENIERÍA SÍSMICA

El bucle principal recorre los diferentes valores de amortiguamiento, donde se llaman a las funciones que
implementan el método ZOH para obtener las respuestas del sistema. Se calcula el pseudo-
desplazamiento como "pseudo_desplazamiento (j, i) = max (X(1, :));" y la pseudo-aceleración como
"pseudo_aceleracion(j, i) = max(abs(Xp(2, :)));". Estos cálculos permiten extraer la máxima respuesta en
desplazamiento y aceleración para cada periodo y cada razón de amortiguamiento.

Finalmente, se grafican los resultados obtenidos, mostrando en la primera subgráfica el pseudo-

desplazamiento y en la segunda la pseudo-aceleración, facilitando la comparación entre los diferentes
niveles de amortiguamiento en función de los períodos naturales. Las líneas "xlabel ('Período (s)');" y
"ylabel('Pseudo-desplazamiento (m)');" definen las etiquetas de los ejes, y "title('Espectro de Respuesta
- Pseudo-desplazamiento');" proporciona el título correspondiente a cada gráfico.

17
Métodos numéricos de integración - Espectros de respuesta elásticos
 COC-3100
INGENIERÍA SÍSMICA

Análisis de la gráfica obtenida a partir del código:

Ilustración 7: Gráfica del espectro de respuesta de pseudo desplazamiento y pseudo aceleración a partir de las
condiciones del enunciado. Fuente: Elaboración propia.

Espectro de Pseudo-desplazamiento:

- Comportamiento general: En la gráfica del pseudo-desplazamiento, se observa que el valor

máximo se alcanza para períodos alrededor de 1.5 a 2 segundos, lo que indica la frecuencia natural
del sistema con mayor respuesta dinámica. La respuesta disminuye a medida que el período
aumenta, estabilizándose después de 4 segundos.
- Influencia del amortiguamiento: A medida que aumenta el amortiguamiento (ξ), la amplitud de la
respuesta disminuye. Esto se puede ver en las curvas donde ξ=0.02 (azul) tiene una mayor
amplitud en comparación con ξ=0.05 (rojo) y ξ=0.07 (verde). Este comportamiento es esperado,
ya que el amortiguamiento tiende a disipar energía, reduciendo la amplitud de las oscilaciones en
la estructura.
- Periodo crítico: La respuesta máxima ocurre en un período específico cercano a 1.5 s para todos
los valores de ξ, aunque la magnitud de la respuesta varía. Esto indica que el sistema es más

18
Métodos numéricos de integración - Espectros de respuesta elásticos
 COC-3100
INGENIERÍA SÍSMICA

vulnerable a excitaciones de este rango de períodos, lo cual puede ser un período resonante del
sistema. La diferencia entre las curvas indica la importancia de la amortiguación en mitigar los
efectos de resonancia.
-
Espectro de Pseudo-aceleración:

- Comportamiento general: La gráfica de pseudo-aceleración muestra un patrón similar al pseudo-

desplazamiento, pero en este caso, la aceleración tiende a ser más pronunciada para los períodos
cortos, con una reducción gradual a medida que el período aumenta.

- Influencia del amortiguamiento: Nuevamente, se observa que para valores de ξ más altos, las
amplitudes disminuyen considerablemente. La curva correspondiente a ξ=0.02 presenta los
valores más altos de pseudo-aceleración en comparación con las otras dos curvas (ξ=0.05 y
ξ=0.07). Esto es consistente con la observación de que un menor amortiguamiento permite
mayores amplitudes de respuesta dinámica.

- Frecuencia alta y baja: Para períodos menores a 1 segundo, las aceleraciones son significativas,
indicando que, en estructuras con períodos cortos, las fuerzas inerciales pueden ser críticas. Este
comportamiento es común en edificaciones rígidas y elementos no estructurales. Para períodos
mayores a 4 segundos, la respuesta en pseudo-aceleración tiende a ser pequeña, reflejando una
menor respuesta dinámica en estructuras flexibles.

Relación entre pseudo-desplazamiento y pseudo-aceleración:

- Máximo en el pseudo-desplazamiento: A medida que aumenta el período, la pseudo-aceleración

decrece de manera más rápida que el pseudo-desplazamiento, lo cual es consistente con las
ecuaciones de movimiento, donde la aceleración está inversamente relacionada con el período
cuadrado.
- Resonancia: El máximo en pseudo-desplazamiento alrededor de los 2 segundos está
correlacionado con un aumento en pseudo-aceleración, pero los máximos no ocurren
exactamente en el mismo período. Este desfase se debe a que la resonancia máxima en aceleración
tiende a ocurrir en períodos menores debido a las relaciones entre la frecuencia y las fuerzas
inerciales.

Conclusiones respecto al comportamiento del período:

- Sistemas con bajos períodos (T<1.5s): experimentan aceleraciones más significativas, lo que
sugiere que para estructuras rígidas (como edificios de baja altura), las fuerzas inerciales son más
pronunciadas, particularmente si el amortiguamiento es bajo.
- Sistemas con períodos intermedios (1.5s<T<3s): muestran los máximos más altos en pseudo-
desplazamiento, lo que sugiere que estas estructuras (edificios de altura media) son más
vulnerables a grandes desplazamientos cuando están sometidas a excitaciones con períodos en
este rango. El aumento del amortiguamiento reduce significativamente el desplazamiento.
- Sistemas con períodos largos (T> 4s): tienen una respuesta moderada en desplazamiento y

19
Métodos numéricos de integración - Espectros de respuesta elásticos
 COC-3100
INGENIERÍA SÍSMICA

aceleración. Esto indica que, para estructuras flexibles o de gran altura, los efectos de las
aceleraciones son menos críticos, pero el desplazamiento podría ser aún considerable si el
amortiguamiento no es lo suficientemente alto. Este análisis demuestra que, dependiendo del
período natural de la estructura y del nivel de amortiguamiento, se pueden tomar decisiones de
diseño informadas para mitigar los efectos de movimientos sísmicos mediante la adición de
amortiguamiento o cambios en la rigidez estructural.

20
Métodos numéricos de integración - Espectros de respuesta elásticos
 COC-3100
INGENIERÍA SÍSMICA

Problema 3. (Estructura con disipador lineal)

21
Métodos numéricos de integración - Espectros de respuesta elásticos
 COC-3100
INGENIERÍA SÍSMICA

3a.

La ecuación de movimiento está dada por medio de la siguiente expresión:

𝑚𝑢̈ (𝑡) + 𝑐𝑢̇ (𝑡) + 𝑘𝑢(𝑡) + 𝑓𝑑 (𝑡) = −𝑚𝑢̈ 𝑔 (𝑡) (𝑎)

Se determina que 𝑓𝑑 (𝑡) es la fuerza del disipador y y se caracteriza por medio de la siguiente expresión:

𝑓𝑑 (𝑡) = 𝑐𝑑 𝑧̇ (𝑡) = 𝑘𝑑 (𝑢(𝑡) − 𝑧(𝑡)) (𝑏)

Para determinar la ecuación de primer orden se considera el siguiente vector de estado 𝑋(𝑡) como:

𝑢(𝑡)
𝑋(𝑡) = [𝑢̇ (𝑡)]
𝑧(𝑡)
Donde:

𝑋1 (𝑡) = 𝑢(𝑡) (desplazamiento de la masa)

𝑋2 (𝑡) = 𝑢̇ (𝑡) (velocidad de la masa)
𝑋3 (𝑡) = 𝑧(𝑡) (deformación del disipador)

Der 𝑋(𝑡), se obtiene:

𝑢̇ (𝑡)
̇
𝑋(𝑡) = [𝑢̈ (𝑡)]
𝑧̇ (𝑡)
A partir de (𝑏) se despeja 𝑧̇ (𝑡) y se obtiene:

𝑘𝑑 (𝑢(𝑡) − 𝑧(𝑡))
𝑧̇ (𝑡) = (𝑐)
𝑐𝑑
Por medio de la expresión en (𝑎), se tiene:

𝑚𝑢̈ (𝑡) = −𝑐𝑢̇ (𝑡) − 𝑘𝑢(𝑡) − 𝑓𝑑 (𝑡) − 𝑚𝑢̈ 𝑔 (𝑡)

Se reemplaza 𝑓𝑑 (𝑡) de la expresión (𝑏), de la cual se obtiene:

𝑚𝑢̈ (𝑡) = −𝑐𝑢̇ (𝑡) − 𝑘𝑢(𝑡) − 𝑐𝑑 𝑧̇ (𝑡) − 𝑚𝑢̈ 𝑔 (𝑡) (𝑑)

Ya que se tiene que:

𝑢̇ (𝑡) = 𝑋2 (𝑡)

Por lo tanto, se obtiene la siguiente igualdad,

22
Métodos numéricos de integración - Espectros de respuesta elásticos
 COC-3100
INGENIERÍA SÍSMICA

𝑋̇1 (𝑡) = 𝑋2 (𝑡)

Utilizando la expresión generada en (𝑑), se tiene:

𝑚𝑢̈ (𝑡) = −𝑐𝑋2 (𝑡) − 𝑘𝑋1 (𝑡) − 𝑐𝑑 𝑋̇3 (𝑡) − 𝑚𝑢̈ 𝑔 (𝑡) (𝑒)

Usando la ecuación del disipador se debe expresar 𝑋̇3 (𝑡) de esta manera según lo obtenido en (c) queda
como:

𝑘𝑑 (𝑋1 (𝑡) − 𝑋3 (𝑡))

𝑧̇ (𝑡) = 𝑋̇3 (𝑡) =
𝑐𝑑

Posteriormente se reemplaza en 𝑋̇3 (𝑡) con lo obtenido en (𝑒), se obtiene que:

𝑚𝑢̈ (𝑡) = −𝑐𝑋2 (𝑡) − 𝑘𝑋1 (𝑡) − 𝑘𝑑 (𝑋1 (𝑡) − 𝑋3 (𝑡)) − 𝑚𝑢̈ 𝑔 (𝑡)

Luego se normaliza la ecuación por medio del factor 𝑚 y se tiene que 𝑋̇2 (𝑡) de esta manera finalmente se
obtiene que:

𝑘 𝑐 𝑘𝑑 (𝑋1 (𝑡) − 𝑋3 (𝑡))

𝑋̇2 (𝑡) = − 𝑋1 (𝑡) − 𝑋2 (𝑡) − + 𝑢̈ 𝑔 (𝑡)
𝑚 𝑚 𝑚

Se definen la matriz 𝐴 y el vector 𝐵 de la siguiente forma:

0 1 0
𝑘 𝑘𝑑 𝑐 𝑘𝑑
− −
𝐴= 𝑚 𝑚 𝑚 𝑚
𝑘𝑑 𝑘𝑑
0 −
[ 𝑐𝑑 𝑐𝑑 ]
0
1
𝐵=[ ]
𝑚
0
Finamente la forma matricial de las ecuaciones de movimiento se define como una ecuación diferencial
de primer orden que se aprecia en la siguiente expresión:

0 1 0
𝑘 𝑘𝑑 𝑐 𝑘𝑑 0
− − 1
𝑋̇(𝑡) = 𝑚 𝑚 𝑚 𝑚 𝑋(𝑡) + [ ] 𝑢̈ 𝑔 (𝑡)
𝑘𝑑 𝑘𝑑 𝑚
0 − 0
[ 𝑐𝑑 𝑐𝑑 ]

23
Métodos numéricos de integración - Espectros de respuesta elásticos
 COC-3100
INGENIERÍA SÍSMICA

3b. Se definen los parámetros del sistema, incluyendo la masa, la rigidez, el amortiguamiento y los
parámetros del disipador. La rigidez del disipador es una fracción de la rigidez del sistema principal,
mientras que el coeficiente de amortiguamiento del disipador es el doble del del sistema. Se carga el
registro sísmico desde un archivo, y la fuerza sísmica se calcula en función de la aceleración del suelo
durante el sismo. El sistema se modela mediante matrices en el espacio de estado, donde las variables de
estado incluyen el desplazamiento, la velocidad y el desplazamiento del disipador. La integración
temporal se realiza utilizando el método de ZOH, que permite calcular el estado del sistema en cada paso
de tiempo, y los resultados se almacenan en una matriz para registrar el desplazamiento, la velocidad y
el desplazamiento del disipador en cada instante. Finalmente, la fuerza del disipador se calcula como la
diferencia entre el desplazamiento del sistema y el del disipador, multiplicada por la rigidez del disipador.

Es así como de lo reseñado con anterioridad a partir del enunciado se determina el siguiente código:

Ilustración 8: Código para desplazamiento u(t) y ciclo del disipador. Fuente: Elaboración propia.

24
Métodos numéricos de integración - Espectros de respuesta elásticos
 COC-3100
INGENIERÍA SÍSMICA

Observaciones importantes del código: Comienza con la definición de los parámetros del disipador. En
las líneas "kd = k / 5;" y "cd = 2 * c;", se establecen la rigidez y el amortiguamiento del disipador. Estos
valores son cruciales para determinar cómo se comportará el sistema frente a la excitación sísmica, ya
que el disipador lineal está diseñado para absorber parte de la energía del movimiento.

Las ecuaciones de movimiento del sistema se expresan en términos de espacio de estado mediante las
matrices "A" y "B". En la línea "A = [0 1 0; -k/m -c/m kd/m; kd/cd 0 -kd/cd];", se define la matriz del
sistema que relaciona el desplazamiento, la velocidad y el desplazamiento del disipador. Esta matriz
captura la dinámica del sistema y permite calcular la respuesta temporal del mismo. Por otro lado, "B =
[0; -1/m; 0];" representa la influencia de las fuerzas externas aplicadas en el sistema.

Durante la integración temporal, la línea "X = expm(A * dt) * X + inv(A) * (expm(A * dt) - eye(3)) * B *
p(i);" utiliza la matriz exponencial para avanzar el estado del sistema en cada paso de tiempo. Esta
expresión es fundamental, ya que permite calcular cómo evolucionan el desplazamiento del sistema y la
fuerza del disipador en respuesta a la fuerza sísmica.

Finalmente, la fuerza del disipador se calcula en la línea "fd = kd * (u - z);", donde se relaciona la diferencia
entre el desplazamiento del sistema y el desplazamiento del disipador. Esta relación es crucial para
entender cómo el disipador reduce el movimiento del sistema principal.

Análisis de la gráfica obtenida por medio del código ejecutado en el software Matlab:

Ilustración 9: Desplazamiento u(t) y ciclo del disipador. Fuente elaboración propia.

Se observa el desplazamiento 𝑢(𝑡) del sistema a lo largo del tiempo. Esta gráfica muestra el
comportamiento oscilatorio del desplazamiento de una estructura sometida a una excitación sísmica,
reflejando la respuesta dinámica de un sistema con un disipador de energía. Durante los primeros

25
Métodos numéricos de integración - Espectros de respuesta elásticos
 COC-3100
INGENIERÍA SÍSMICA

segundos, la amplitud del desplazamiento incrementa, lo cual sugiere una excitación inicial fuerte,
seguida de una disminución gradual conforme el sistema se estabiliza y la energía se disipa a través del
tiempo. La oscilación tiende a reducirse alrededor del valor cero, lo que indica que el disipador está
funcionando eficazmente al reducir las vibraciones hasta casi eliminarlas después de cierto tiempo.

En el segundo gráfico de esta figura, se presenta la relación entre la fuerza del disipador 𝑓𝑑(𝑡) y el
desplazamiento 𝑧(𝑡). La forma cerrada de las curvas elípticas en este plano revela el comportamiento
cíclico típico de un disipador de energía que sigue una ley constitutiva no lineal, donde la fuerza disipada
es proporcional al desplazamiento. Las elipses más amplias sugieren mayores niveles de disipación en las
primeras oscilaciones, mientras que las más estrechas hacia el centro muestran cómo el sistema se
estabiliza al disipar la energía almacenada en el desplazamiento.

El gráfico muestra la evolución del desplazamiento del sistema y la relación entre la fuerza del disipador
y el desplazamiento. En el ámbito de la ingeniería sísmica, esto es fundamental, ya que:

- Disminución de la amplitud del desplazamiento: La gráfica de desplazamiento u(t) refleja cómo

la estructura es capaz de reducir sus oscilaciones con el tiempo, gracias a la incorporación de un
disipador de energía. Esto es crucial para reducir el daño estructural durante y después de un
evento sísmico.
- Ciclos de disipación de energía: La gráfica de fuerza vs desplazamiento (curvas elípticas) indica
un ciclo de histéresis, lo cual es característico de los dispositivos de disipación. Este
comportamiento asegura que la energía sísmica se disipa de manera controlada, evitando que se
almacene en la estructura y cause deformaciones permanentes o fallos.

26
Métodos numéricos de integración - Espectros de respuesta elásticos
 COC-3100
INGENIERÍA SÍSMICA

3c. Se examina la variación de la energía en el sistema a lo largo del tiempo. Para ello, se calculan varias
formas de energía, incluyendo la energía cinética, que se determina a partir de la velocidad del sistema
1 1
mediante la fórmula 2 𝑚𝑢̇ 2 . Además, se evalúa la energía elástica, que se calcula como 2 𝑘u2 , donde 𝑢 es
el desplazamiento del sistema. También se incluye la energía viscosa, representando la energía disipada
por el amortiguamiento, que se obtiene integrando la fuerza de amortiguamiento multiplicada por la
velocidad a lo largo del tiempo. Por otro lado, se determina la energía del disipador, que se calcula a partir
del trabajo realizado por la fuerza del disipador en función del desplazamiento. Por último, se evalúa la
energía del input sísmico, que se obtiene del trabajo realizado por las fuerzas sísmicas sobre el sistema
durante el movimiento. Cada uno de estos componentes energéticos se calcula en cada paso de tiempo y
se almacena para su análisis posterior. A partir de lo anterior se obtiene el siguiente código que se
muestra a continuación:

Ilustración 10: Código para variación de energías en el sistema durante el evento sísmico. Fuente: Elaboración propia.

27
Métodos numéricos de integración - Espectros de respuesta elásticos
 COC-3100
INGENIERÍA SÍSMICA

Finalmente, los resultados se grafican, permitiendo observar cómo varían la energía cinética, la energía
elástica, la energía viscosa, la energía del disipador y la energía del input sísmico a lo largo del tiempo, lo
que proporciona una comprensión más profunda del comportamiento dinámico del sistema estructural
ante la excitación sísmica.

Ilustración 11: Gráfico para variación de energías en el sistema durante el evento sísmico. Fuente: Elaboración propia.

En el gráfico anterior se presenta la variación de las energías con el tiempo, y en ella se muestran diversas
formas de energía: la energía cinética, energía viscosa, energía elástica, energía del disipador y la energía
de input sísmico. El gráfico ilustra cómo la energía sísmica aplicada al sistema se transforma y se
distribuye en las diferentes formas de energía a medida que el sistema responde al evento sísmico. La
energía de input sísmico (línea negra) se incrementa rápidamente al inicio, lo cual representa el impulso
inicial de la excitación sísmica sobre el sistema. Posteriormente, esta energía se transfiere a otras formas
de energía. La energía elástica (línea verde) crece a un ritmo estable, indicando que una parte importante
de la energía sísmica se almacena en la estructura mediante deformaciones elásticas. La energía cinética
(línea roja) muestra picos en las primeras etapas, lo que refleja el movimiento acelerado del sistema. Es
así como, la energía viscosa y la energía disipada por el disipador (líneas morada y azul, respectivamente)
crecen de manera continua, evidenciando la capacidad del disipador para absorber y disipar la energía
mecánica, lo que contribuye a reducir las oscilaciones.

Finalmente, muestra la variación de diferentes formas de energía en el tiempo, aporta información clave
sobre el comportamiento energético del sistema ante un evento sísmico:

- Energía disipada y controlada: La energía cinética y la energía de input sísmico, que inicialmente
se presentan de forma elevada, son absorbidas en su mayoría por el sistema de disipación. Esto

28
Métodos numéricos de integración - Espectros de respuesta elásticos
 COC-3100
INGENIERÍA SÍSMICA

es fundamental, ya que la disipación efectiva reduce el impacto de las aceleraciones y, por lo tanto,
disminuye las fuerzas internas en la estructura.
- Distribución energética: El hecho de que la energía viscosa y la energía del disipador crezcan
continuamente indica que gran parte de la energía sísmica es absorbida por estos dispositivos.
Esto evidencia su capacidad para evitar que la energía se acumule en la estructura en forma de
deformaciones permanentes o fallos.
- Control de la deformación: La energía elástica también se acumula en el sistema, lo cual es
indicativo de la deformación elástica que puede ser recuperada. Esto demuestra que la estructura
trabaja dentro de un rango elástico, es decir, que las deformaciones no son irreversibles, un
objetivo crucial en el diseño sismorresistente.

29
Métodos numéricos de integración - Espectros de respuesta elásticos
 COC-3100
INGENIERÍA SÍSMICA

3. Conclusiones.

Se destaca que se implementaron y compararon los métodos de integración numérica de Newmark y

Zero-Order Hold (ZOH) para analizar la respuesta dinámica de un sistema de un grado de libertad (SDOF)
sometido a la excitación sísmica del terremoto de Concepción de 2010. Ambos métodos fueron efectivos
en la obtención de desplazamientos y velocidades, pero el método de Newmark mostró una estabilidad
superior, particularmente en sistemas con bajo amortiguamiento, debido a su naturaleza implícita e
incondicionalmente estable para los parámetros seleccionados (β = 1/6 y γ = 1/2). El ZOH, siendo un
método explícito, demostró ser más sensible a la discretización temporal, pudiendo volverse inestable
con pasos de tiempo grandes. Los resultados indicaron que las oscilaciones iniciales del sistema, con una
razón de amortiguamiento crítico ξ = 0.02, fueron más pronunciadas, lo que refleja la baja disipación de
energía en el sistema. Además, se observó que ambos métodos arrojaron respuestas similares en
términos de desplazamiento y velocidad, pero las diferencias se hicieron más evidentes al comparar las
aceleraciones, donde el método de Newmark presentó resultados más coherentes con la física del
sistema. Por último, el análisis permitió identificar la sensibilidad del sistema a ciertas frecuencias de
resonancia y la importancia del amortiguamiento en la mitigación de las respuestas dinámicas más
críticas.

30
Métodos numéricos de integración - Espectros de respuesta elásticos