Asignatura: Ingeniería Sismorresistente
Unidad 2: Respuesta elástica de sistemas de un grado de libertad (GLD).
Movimiento en la base. Métodos numéricos para el cálculo de la respuesta
SESIÓN 7
Eje temático 7. Métodos numéricos para el cálculo de la respuesta-Newmark
Eje temático 8. Construcción de espectros de respuesta
Mg. Ing. Orlando Maldonado Salvatierra
https://chat.openai.com/chat
Equilibrio dinámico para sistemas de 1GDL
F(t) – k.x(t) – c.x'(t) – m.x''(t) = 0
No amortiguado c = 0
F(t) =0 Vibración libre
Amortiguado cǂ0
No amortiguado c = 0
F(t) ǂ 0 Vibración Forzada
Amortiguado cǂ0
Vibración armónica de sistemas amortiguados
En la figura 3.4 pueden observarse diversas funciones de carga. De éstas, nos
interesan por ahora, las periódicas y más particularmente las excitaciones
armónicas ya que mediante series de Fourier cualquier excitación periódica puede
llevarse a una suma de armónicas simples.
SISTEMAS DE UN GRADO DE LIBERTAD (1GDL)
Solución a la respuesta y características de la
respuesta libre.
En particular, se recurre al método de la transformada
de Laplace para determinar la solución de un sistema
lineal y la integración numérica se aplica para hallar la
solución de un sistema no lineal.
Vibraciones transitorias
F(t) – k.x(t) – c.x'(t) – m.x''(t) = 0
La determinación de la respuesta de un sistema de 1GDL que se
ve afectado por una excitación que no es ni periódica ni armónica
presenta un grado de complejidad mayor. En muchos casos
prácticos donde se tiene excitaciones que no se prestan a una
descripción matemática hay necesidad de recurrir a métodos
numéricos para obtener la solución.
Métodos numéricos
Los métodos numéricos son una sucesión de operaciones
matemáticas utilizadas para encontrar una solución numérica
aproximada a un problema determinado. Es decir, se trata de una
serie de cálculos para acercarnos lo más posible a una solución
numérica con una precisión razonablemente buena.
Los métodos numéricos son utilizados en ingeniería para facilitar
la resolución de problemas que conllevan una enorme cantidad
de cálculos, lo que permite ahorrar tiempo.
Excitación arbitraria
Cuando un sistema se somete a una excitación arbitraria en términos de
fuerza, es posible dividirla en una serie de impulsos que se aplica en el tiempo
y que tiene una duración
Excitación arbitraria
Integrando, el efecto de cada uno de estos impulsos diferenciales variando
Se obtiene:
Para un sistema no amortiguado
Para un sistema no amortiguado
Para un sistema amortiguado Para un sistema amortiguado
Estas integrales se conocen como integrales de convolución o de Duhamel, y
corresponden a la solución particular del sistema. Si hay condiciones iniciales
hay necesidad de adicionar la solución homogénea (transitoria).
Excitación en la base
Es el caso en el cual la excitación del sistema proviene de un movimiento en
su base es muy importante en la dinámica estructural, pues la excitación
sísmica induce este tipo de respuestas del sistema.
Excitación en la base
Cambio de variable:
La ecuación indica que un sistema al que se le introduce movimiento en su base es equivalente a
un sistema con su base fija al cual se le aplica una fuerza igual a la masa del sistema multiplicado
por el negativo de la aceleración del terreno.
Obtención de la respuesta dinámica
En la medida que la excitación se hace más compleja su tratamiento
por medio de la integral de convolución se hace mas difícil, si no
imposible. Por esta razón hay necesidad de recurrir a los métodos
numéricos de solución, los cuales utilizan diferentes técnicas y son
apropiados para distintos problemas especiales fácilmente
solucionables por medio del computador digital.
Métodos numéricos para la obtención
de la respuesta dinámica
-Método paso a paso en el tiempo
-Método basado en la interpolación de
la excitación
-Método de la diferencia central
-Método de Newmark
Integración numérica de la ecuación de
movimiento
Se trata de resolver en forma discreta la ecuación de movimiento para 1 GLD
dada por:
utilizando ecuaciones en diferencias que permitan obtener:
en un instante t+∆t en función, únicamente, de los valores en el instante t. Existen
varios métodos para plantear las ecuaciones en diferencias antes mencionadas,
pero solo nos dedicaremos en profundidad al de Newmark, dado que es uno de
los más difundidos en software de análisis dinámico mediante FEM.
Respuesta 1GDL
Método de Newmark (1 GLD)
Desarrollo y forma operativa
La ecuación de movimiento valuada en t=ti+1 toma la forma:
Reemplazando en ésta los valores en diferencias dados por (1) y (2), agrupando
los valores que permanecen constantes en cantidades denominadas por a0,…,
a7; se obtiene el algoritmo de integración paso a paso.
Algoritmo de integración paso a paso.
Algoritmo de integración paso a paso.
Respuesta 1GDL
Respuesta 1GDL
Respuesta 1GDL
Ejemplo
Hallar la respuesta en el tiempo del oscilador indicado en la figura 2.1, que tiene una masa
m = 0.004898 Tms2/cm, una frecuencia natural Wn = .6 2832 1/s y un coeficiente de
amortiguamiento ξ = 0.05. Ante el sismo del 9 de noviembre de 1974, registrado en el Perú.
El acelerograma fue obtenido a 80 Km., del epicentro sobre un suelo limo arcilloso. El
evento tuvo una magnitud de 6. La aceleración máxima del sismo, en valor absoluto fue de
117 gals (cm/s2). El incremento de tiempo con el cual fue obtenido el registro es ∆t = 0.02 s.
Solución
Ejemplo
Espectros de respuesta
Se define el espectro de respuesta como la respuesta máxima de un conjunto de
osciladores de 1 GDL que tienen el mismo amortiguamiento, sometidas a una
historia de aceleraciones dadas.
Construcción de espectros de respuesta
Los espectro de respuesta (para un acelerograma dado) es muy útil para el diseño
de estructuras donde solo los máximos son necesarios.
Para un sistema de 1GDL
La solución a esta ecuación es:
……….. (4)
Construcción de espectros de respuesta
y el máximo valor absoluto será Sd
Reemplazando
Puede obtenerse, para 1 GDL, la expresión exacta del espectro de
desplazamientos. Derivando una y dos veces se obtienen las expresiones para
velocidad y aceleración:
Construcción
de espectros de
respuesta
Construcción de espectros de respuesta
Construcción de espectros de respuesta
No obstante pueden hacerse ciertas hipótesis que simplifican las expresiones
obtenidas. Estas nuevas ecuaciones definen los “Seudo espectros de
respuesta”.
Hipótesis simplificativas:
i- en estructuras civiles, 0,002 < ξ < 0,2, por lo que puede reemplazarse ωv por ω
ii- la función coseno que aparece en el espectro de velocidades (al derivar ecuación 1)
puede reemplazarse, sin que ello implique grandes variaciones en los valores máximos,
por una tipo seno entonces:
Código en Matlab–espectro
de respuesta
Ejemplo
Hallar los espectros de
respuesta elástica de
desplazamiento, velocidad y
aceleración del sismo del 9
de noviembre de 1974, para
valores de ξ igual a: 0.05;
0.10 y 0.20.
Construcción de espectros de respuesta
Las ecuaciones anteriores definen los seudo espectros para desplazamientos
(Sd), velocidades (Sv) y aceleraciones (Sa).
Es importante notar que:
Por lo que un espectro
completo puede representarse
en un gráfico tipo logarítmico
como el de la figura.
Figura 6.6.4 Espectro de respuesta combinado D-V-A para el movimiento del terreno
de El Centro; ζ = 0, 2, 5, 10 y 20%. (Chopra)
Importancia De las formas espectrales
Los espectros de respuesta proporcionan información muy valiosa para el
proyectista estructural, ya que se puede inferir los edificios que van a estar
sujetos a mayores fuerzas sísmicas.
Respuesta elásticas espectrales
Para encontrar un espectro de diseño se deben clasificar los registros sísmicos
de acuerdo al lugar en que fueron registrados ya que la forma espectral depende
del tipo de suelo.
Una vez que se tienen clasificados los eventos se procede a obtener los
espectros de respuesta de cada uno de ellos, posteriormente se aplican las
estadísticas con las que se determina el espectro de diseño. Realmente es muy
sencillo encontrar un espectro de diseño lo difícil es tener una muestra de datos
que se la pueda considerar confiable
Respuesta elásticas espectrales
Es deseable que los registros sísmicos con los cuales se vayan a obtener los
espectros de diseño tengan una aceleración máxima de suelo considerable, por
lo menos que sean mayores al 10% de la aceleración de la gravedad.
En la mayor parte de países de Latinoamérica no se cuenta con una cantidad
suficiente de eventos fuertes por lo que han trabajado con sismos de
aceleraciones pequeñas normalizados a aceleraciones grandes, este
procedimiento no es correcto pero ante la ausencia de registros fuertes no
queda otra opción.
Ejemplo
Obtener un espectro de diseño a partir de los registros sísmicos indicados en la
tabla 3.1, que fueron sentidos o registrados en el Perú. En la última columna se
muestra el tipo de suelo en el cual se obtuvo el registro.
Ejemplo
En la tabla 3.1 se tiene un total de 17 registros, cantidad que es pequeña como para pensar en
separarlos de acuerdo al tipo de suelo en que fueron registrados, razón por la que se trabaja con
todos ellos. Cada uno de estos registros fue normalizado a 392 gals (0.4 g) de tal manera que los
registros se multiplicaron por un factor tal que la aceleración máxima sea la indicada. Los
espectros se obtuvieron para ξ = 0.05
Al trabajar con el espectro medio se tiene que la probabilidad de excedencia de las ordenadas
espectrales es del 50%. En efecto, se aprecia que existe una cantidad significativa de aceleraciones
que están sobre la curva media. Si se desea disminuir esta probabilidad de excedencia a la curva de
valores medios se deberá sumar una desviación estándar o más dependiendo de la probabilidad de
excedencia con la cual se desea trabajar.
Chávez et al (2012) efectuó una comparación de los espectros de pseudo-aceleración en
diferentes países en latino América en los que se encuentra el Perú, haciendo una adaptación de
los diferentes espectros consiguió una parametrización.
El factor de zona y de suelo (Z y S) y los periodos TP y TL, se obtienen del RNE, Norma E.0.30
Tabla N°1, Tabla N°3 y Tabla N°4, respectivamente. Para la zona Z3, le corresponde Z = 0.35,
para un suelo intermedio el factor de suelo es S= 1.15 y los periodos TP(s) = 0.6 y TL (s) = 2.0,
con estos datos se realiza el espectro de pseudo aceleración y espectro de desplazamiento.
Síntesis
En la medida que la excitación se hace más compleja hay necesidad
de recurrir a los métodos numéricos de solución, los cuales utilizan
diferentes técnicas y son apropiados para distintos problemas.
Los espectros de respuesta se obtienen con la respuesta máxima de
un conjunto de osciladores de 1 GDL que tienen el mismo
amortiguamiento, sometidas a una historia de aceleraciones dadas. Es
muy útil para el diseño de estructuras donde solo los máximos son
necesarios.
