Tratamiento de Acelerogramas
Tratamiento de Acelerogramas
Tratamiento de Acelerogramas
DE LA TIERRA
SISMOLOGA
E
INGENIERA SSMICA
(PRCTICAS)
Anexo VI
Prcticas de Sismologa e Ingeniera Ssmica
An. VI - P5 - Pgina 1
Figura 1 Estructura sometida a una fuerza
lateral.
PRACTICA 5.
TRATAMIENTO DE ACELEROGRAMAS.
1. OBJETIVO
Aprender a llevar a cabo una determinacin lo ms exacta posible de los valores de
aceleracin que se han producido en el lugar de registro durante un terremoto y obtener
aquellos parmetros de inters en la ingeniera sismolgica que se obtienen a partir de los
registros de aceleracin. Aprender a utilizar el programa DEGTRA (2000) (Ordaz y
Montoya, 2000) para llevar a cabo dichos objetivos.
2. TEORIA
Tanto el registro de acelerogramas como su tratamiento han ido cambiando a
medida que han evolucionado los sistemas de adquisicin y el almacenamiento y procesado
de los datos.
El conocimiento de las fuerzas que actan sobre un emplazamiento, debida a los
terremotos, el anlisis de sus efectos sobre las estructuras, el desarrollo de normas de
construccin y los mtodos para determinar el riesgo ssmico sobre las construcciones son
algunas de las principales motivaciones para el estudio de acelerogramas.
En ellos, el parmetro de mayor inters es la aceleracin del suelo, aunque tambin
son importantes, la velocidad y el desplazamiento. As si suponemos una estructura simple
con una masa concentrada en su tejado (Figura 1) y consideramos la fuerza lateral
concentrada en los elementos verticales, que seran las paredes o columnas. Aplicando las
leyes de Newton del movimiento tenemos que :
mx f x x t mz t && ( , & , ) &&( ) + =
ecuacin que nos da la respuesta al terremoto como
si la estructura estuviera fija por la base y sometida
a una fuerza dinmica , es la mz t &&( ) f x x t ( , & , )
fuerza restauradora del edificio que es funcin de la
desviacin relativa, x(t) y de su derivada. As, la
fuerza efectiva de carga debida al terremoto en un
punto de la masa es igual a tantas veces la masa por
la aceleracin del suelo. De esta ecuacin vemos
que los edificios construidos, por ejemplo, en
Prctica 5.
Tratamiento de acelerogramas
An. VI - P5 - Pgina 2
piedra, y caracterizados por grandes masas y baja resistencia, son particularmente
peligrosos durante los grandes terremotos, al menos en ausencia de otros factores.
Estos registros se denominan de campo cercano, ya que el estudio se realiza en la
regin epicentral, en la que los movimientos en superficie al paso de las ondas ssmicas son
ms intensos (movimientos fuertes del suelo), producindose desplazamientos, velocidades
y aceleraciones elevadas, aunque dependiendo fundamentalmente de la magnitud del sismo
y de su distancia, al epicentro.
Para obtener una determinacin lo ms exacta posible de los valores de aceleracin
que se han producido durante el terremoto en el lugar del registro, es necesario hacer un
cuidadoso tratamiento de los registros. Puesto que los datos que se recogen de los
acelergrafos estn en un formato propio del registrador (binario), ser necesario hacer una
conversin de los mismos a formato ASCII con el que, generalmente trabajaremos.
El acelerograma obtenido, presenta una serie de errores que debern ser corregidos
para que dicho registro sea apto para ser utilizado. En la Figura 2 presentamos el registro
en sus tres componentes, sin corregir, de un sismo de magnitud 4.3.
Figura 2 Acelerograma del terremoto de Ventas de Huelma (GR) (24/02/97) con mb
= 4.3 en el acelergrafo de Agrn (GR). Arriba componente N-S, centro: Componente
Vertical y Abajo: Componente E-W.
2.1 Correcciones sobre los acelerogramas
Correccin Instrumental: Consiste en la convolucin del registro con la funcin de
Anexo VI
Prcticas de Sismologa e Ingeniera Ssmica
An. VI - P5 - Pgina 3
transferencia del instrumento. Para ello, debemos conocer las constantes de calibracin del
instrumento, como son la sensibilidad del acelermetro, la frecuencia natural T
n
, y el
amortiguamiento >, obtenindose as un acelerograma corregido.
Correccin de la lnea base o correccin de ceros: A veces, ocurre que los valores
del acelerograma se encuentran desplazados respecto a la lnea cero de aceleracin. Este
error puede ocurrir porque el acelergrafo no est perfectamente nivelado en su
emplazamiento, o bien porque el sistema de registro provoque una deriva de los datos
respecto a la lnea de base. Aunque este error puede ser inapreciable en aceleracin, puede
ser muy importante cuando se obtienen la velocidad y el desplazamiento por integracin,
ya que estamos contabilizando el rea que hay entre la curva de aceleracin y la lnea de
base.
2.2 Filtrado del ruido
Una vez que el acelerograma ha sido corregido tanto de lnea base como de
correccin instrumental, habr que tener en cuenta la repercusin que el ruido tiene sobre
el acelerograma. Este ruido puede estar causado tanto por fenmenos naturales como por
el tratamiento de los datos. As podremos considerar que el acelerograma es la suma de la
seal ssmica ms el ruido,
a(t) = s(t) + r(t)
siendo a(t) el acelerograma sin corregir, s(t) la seal ssmica y r(t) el ruido. El ruido que
afecta al acelerograma (Figura 3), puede ser de dos tipos:
Figura 3 Izqda: Efecto del ruido de alta frecuencia; Dcha: Efecto del ruido de baja
frecuencia.
Prctica 5.
Tratamiento de acelerogramas
An. VI - P5 - Pgina 4
a) Ruido de alta frecuencia: Cuando aparece por encima de la banda de frecuencias en la
que trabajamos. Este ruido afecta fundamentalmente a los picos de aceleracin. Sin
embargo, no afecta mucho en el proceso de integracin, ya que el rea bajo la curva es
prcticamente la misma que para la seal corregida, dado que las desviaciones por encima
y por debajo se compensan.
b) Ruido de baja frecuencia: Cuando aparece por debajo de la banda de frecuencias en la
que trabajamos. Afecta al proceso de integracin bsicamente, ya que al integrar la
aceleracin para obtener velocidad y desplazamiento, estamos obteniendo el rea bajo la
curva.
Puesto que no queremos que estos errores nos afecten, limitaremos la informacin
a la banda de frecuencias deseada, en lugar de mantener informacin adicional a mayores
y menores frecuencias, que pueden contener gran nmero de errores. Para hacerlo, se
utilizan los filtros pasa baja para eliminar altas frecuencias, y los pasa alta para eliminar
las bajas. Los filtros pasa banda, eliminan el ruido de alta y baja frecuencia
simultneamente.
Lo ms importante a la hora de aplicar los filtros, es elegir la frecuencia de corte
adecuada, ya que una mala eleccin puede suponer eliminar parte de la seal que nos
interesa, o por el contrario, no eliminar el ruido de la seal.
El filtrado de la seal se hace multiplicando el espectro del acelerograma por una
funcin que reduzca el valor de la amplitud del espectro de Fourier fuera de la banda donde
estn contenidas las frecuencias importantes de nuestra seal. En esa banda, el filtro deja
las amplitudes inalterables.
En la prctica los filtros ideales no se pueden utilizar porque su funcin de
transferencia en el tiempo es de longitud infinita. Sin embargo, se utilizan aproximaciones
a ellos.
El programa Detra2000 permite usar los siguiente filtros:
a) Pasa baja hasta la frecuencia mxima Fmax:
H f
f
F
N
( )
max
=
+
|
\
|
.
|
1
1
2
Anexo VI
Prcticas de Sismologa e Ingeniera Ssmica
An. VI - P5 - Pgina 5
b) Pasa alta desde la frecuencia mnima Fmin:
H f
F
f
N
( )
min
=
+
|
\
|
.
|
1
1
2
c) Pasa banda entre las frecuencias mnima y mxima:
H f
f F F
f F F
N
( )
( )
min max
max min
=
+
|
\
|
.
|
1
1
2
2
d) Band Stop entre las frecuencias mnima y mxima
H f
f
F
F
f
N N
( )
min
max
=
+
|
\
|
.
|
+
+
|
\
|
.
|
1
1
1
1
2 2
e) Gaussiano entre las frecuencias mnima y mxima
H f
f f
S
con f
F F
S f F
C
C C
( ) exp ; ;
min max
min
=
|
\
|
.
|
(
=
+
=
1
2 2
2
f) Futterman con parmetro t*
H f ft it
f
f
c
( ) exp( *) exp( *log ) = 2
donde f
C
es la frecuencia de Nyquist, dada por: f
C
= 1 / 2)t
Este filtro es un filtro complejo por lo que tambin afecta a la fase de la seal.
2.3 Integracin de acelerogramas
Una vez que el registro temporal de aceleracin ha sido convenientemente corregido
para su posterior utilizacin, es preciso integrarlo para obtener los registros de velocidad
y desplazamiento (Figura 4). Si consideramos el dominio del tiempo, la integracin se hace
por mtodos numricos, teniendo en cuenta que:
Prctica 5.
Tratamiento de acelerogramas
An. VI - P5 - Pgina 6
r r r
v a dt v
t
o
= +
0
r r r
r v dt r
t
t
o
o
= +
Cada integracin suaviza los componentes ms altos del movimiento de las ondas,
por lo que los registros de velocidad y desplazamiento son ms simples y permiten una
interpolacin directa de los patrones de las ondas. As, los impulsos de los primeros
movimientos mayores de la onda S aparecen de forma inequvoca en los registros de
velocidad, mientras que en los registros de desplazamiento se puede apreciar la presencia
de un impulso de largo periodo asociado con el paso de la dislocacin cerca del
emplazamiento.
El pico de mxima velocidad nos da una clara medida del contenido de frecuencias
intermedias del terremoto. Los picos de desplazamiento, muestran las caractersticas de
largo perodo del movimiento. Por tanto, con los tres valores pico (de aceleracin,
velocidad y desplazamiento), tenemos una importante descripcin del movimiento del
suelo, mostrando cada uno una regin diferente del espectro de frecuencias.
Figura 4 Aceleracin, velocidad y desplazamiento, obtenidos a partir del acelerograma
registrado en Agrn (24/02/97), mb = 4.3.
En el dominio de la frecuencia, la integracin se reduce a dividir la aceleracin por
-iT para hallar la velocidad y por T
2
para hallar el desplazamiento:
Anexo VI
Prcticas de Sismologa e Ingeniera Ssmica
An. VI - P5 - Pgina 7
L a t A e
L v t
i
A e
L d t A e
i
i
i
[ ( )] ( )
[ ( )] ( )
[ ( )] ( )
( )
( )
( )
=
=
=
1
1
2
donde L [ ] indica transformada de Fourier.
2.5 Anlisis frecuencial
Una vez que se han obtenido los registros de velocidad y desplazamiento, es
necesario hacer un anlisis del contenido en frecuencias del acelerograma. Para ello, se
utilizan el espectro de Fourier y el espectro de respuesta.
Espectro de Fourier
El espectro de Fourier es un parmetro que nos proporciona un amplio
conocimiento acerca de cual es el contenido en frecuencias del acelerograma, de modo que
para cada periodo de oscilacin del acelereograma se representa la amplitud mxima que
le corresponde (Figura 5).
Matemticamente, el espectro de Fourier en amplitudes se define como:
| | | |
FS a d a sen d
t t
= +
( ) cos( ) ( ) ( )
0
2
0
2
De la interpretacin del espectro de Fourier, podemos deducir cuales son los
perodos ms importantes del terremoto ( que son los que lo van a caracterizar), sin ms que
buscar en qu lugares se alcanzan los picos de amplitud. Adems, podemos ver cual es la
frecuencia o periodo predominante de un registro, que es aquel para el cual el espectro de
Fourier de la aceleracin alcanza su valor mximo, y que adems mantiene una correlacin
directa con la distancia epicentral. Se suele representar con ejes logartmicos, de modo que
en el eje de abscisas se colocan los distintos periodos de oscilacin, y en el eje de ordenadas
se representan las amplitudes en cm/s.
Espectro de Respuesta
El espectro de respuesta de un acelerograma, representa la respuesta mxima de un
oscilador de un slo grado de libertad y de masa unitaria, con frecuencia T
o
y
amortiguamiento >, cuando se le somete a un movimiento del suelo dado por dicho
Prctica 5.
Tratamiento de acelerogramas
An. VI - P5 - Pgina 8
Figura 5 Espectro de Fourier del acelerograma registrado
en Agrn (GR) el 24/02/97, mb=4.3.
Figura 6 Modelo de un grado de libertad
sometido a una excitacin ssmica.
acelerograma (Figura 6).
Se suelen considerar distintos
periodos propios del oscilador y
distintos valores de amortiguamiento
de modo que para un periodo propio T
y un amortiguamiento > del oscilador,
el acelerograma provoca un
movi mi ent o cuyos pi cos de
a c e l e r a c i n , v e l o c i d a d y
desplazamiento se recogen en el
espectro de respuesta.
La solucin de las vibraciones
estacionarias de un sistema de un grado
de libertad, viene dada por la integral
de Duhamel, donde x(t) es funcin de
>, T y a(t):
x t a e sen t d
t
t
( ) ( ) [ ( )]
( )
=
1
0
donde:
= 1
2
que derivando nos da la historia de la respuesta en
velocidades:
& ( ) ( ) cos[ ( )] ( )
( )
x t a e t d x t
t
t
=
0
y derivando de nuevo, se obtiene la respuesta en aceleraciones:
| |
&&( ) ( ) ( ) [ ( )] & ( ) ( ) ( )
( )
x t a t a e sen t d x t x t
t
t
+ =
0
2
2
El espectro de respuesta en desplazamiento, velocidad o aceleracin, ser el
conjunto de mximos valores de las respuestas del sistema:
Anexo VI
Prcticas de Sismologa e Ingeniera Ssmica
An. VI - P5 - Pgina 9
S x t
S x t
S x t a t
d
r
v
r
a
r
( , ) ( )
( , ) & ( )
( , ) &&( ) ( )
max
max
max
=
=
= +
donde S
d
r
y S
v
r
son los mximos de desplazamiento y velocidad relativos y S
a
r
es el
mximo de aceleracin total del oscilador de un grado de libertad.
Puesto que estos valores requieren un largo proceso de clculo, se han hecho
simplificaciones para aligerarlo, obtenindose los pseudoespectros de respuesta, que vienen
dados por las siguientes expresiones:
S a e sen t d
S a e sen t d
S a e sen t d
d
t
t
v
t
t
a
t
t
( , ) ( ) [ ( )
( , ) ( ) [ ( )
( , ) ( ) [ ( )
( )
max
( )
max
( )
max
=
=
=
1
0
0
0
Como se puede ver, sin ms que calcular S
d
, podemos obtener S
v
y S
a
, ya que:
S S
S S
v d
a d
=
=
2
Tanto para calcular los espectros ssmicos de respuesta como los pseudoespectros
ssmicos de respuesta es necesario realizar el clculo numrico de la integral de Duhamel,
hacindose determinadas simplificaciones que dan una buena aproximacin y reducen
sensiblemente el coste numrico.
Una forma simple y elegante de representar el espectro de respuesta, es mediante
un diagrama trilogartmico de cuatro entradas en el que aparece el pseudoespectro de
velocidad frente al periodo T o frecuencia T, y en las lneas a 45 se representan los valores
constantes de pseudoaceleracin absoluta para pendiente 1, y desplazamientos constantes
relativos para pendiente -1.
Adems el espectro de respuesta se suele representar para distintos valores del
amortiguamiento (Figura 7).
No obstante, a veces es conveniente representar con ejes reales y por separado, los
Prctica 5.
Tratamiento de acelerogramas
An. VI - P5 - Pgina 10
Figura 7 Espectro de respuesta del acelerograma
registrado en Agrn (GR) el 24/02/97, mb = 4.3
distintos espectros de aceleracin,
velocidad o desplazamiento, y hacerlo
p a r a a q u e l l o s v a l o r e s d e
amortiguamiento que se consideran
necesarios (Figura 8).
Tambin se utilizan mucho los
grficos normalizados, en los que se
divide el valor del pseudoespectro de
aceleracin por el valor pico de
aceleracin en el origen, ya que estos
grficos normalizados son los que se
c o mp a r a n c o n l a n o r ma
sismorresistente de construccin
(Figura 9).
Figura 8 Espectro de aceleracin para un
amortiguamiento del 5% correspondiente al
acelerograma registrado en Agrn (GR) el
24/02/97, mb=4.3
Figura 9 Espectro de respuesta en aceleracin para
un amortiguamiento del 5% correspondiente al
acelerograma registrado en Agrn (GR) el
24/02/97, mb=4.3
2.6 Otros parmetros de inters en ingeniera
Duracin significante. Grficos de Husid
El concepto de duracin significante nos relaciona la duracin con la fase de
mxima energa del movimiento, y queda representado por el tiempo transcurrido entre los
instantes en los que se alcanza el 5% y el 95% del valor de la intensidad de Arias.
La intensidad de Arias es una medida de la energa del terremoto y se define segn
la integral siguiente:
Anexo VI
Prcticas de Sismologa e Ingeniera Ssmica
An. VI - P5 - Pgina 11
I
g
a t dt
A
t
=
2
2
0
( )
siendo a(t) la historia temporal de aceleracin, y t
f
el tiempo total del registro.
Por tanto, si calculamos la intensidad de Arias, y hallamos para que instantes se
alcanzan el 5% y el 95% de su valor, la diferencia entre los mismos nos dar la duracin
significante.
Para ver grficamente su significado (Figura 8), se utilizan los grficos de Husid,
que nos muestran la variacin de la intensidad de Arias frente al tiempo. Es un grfico de
valores acumulados que define la evolucin del movimiento del suelo en el emplazamiento
del acelergrafo.
El grfico se representa normalizado, es decir, colocando en el eje de abscisas los
valores del tiempo transcurrido en segundos, y en el eje de ordenadas los valores
acumulativos que representa el cociente:
I
I
g
a t dt
g
a t dt
t
A
t
t
f
=
2
2
2
0
2
0
( )
( )
Intensidad Espectral o Intensidad de Housner
Puesto que el espectro de respuesta define los valores mximos de la respuesta
ssmica de estructuras de un grado de libertad para un determinado rango de periodos, se
puede deducir que el rea bajo el espectro entre dos periodos cualesquiera constituye una
medida de la respuesta media de las estructuras en ese rango de periodos.
Por eso, Housner defini una medida de la energa local del movimiento que
provoca un terremoto, a la que llam intensidad espectral.
Los lmites de integracin y el valor de > se eligen de forma que incluyan el rango
de periodos ms habitual en el diseo sismorresistente de las estructuras de edificacin, por
lo que son normales los valores de 0.1 y 2.5 s, como lmites de integracin, y un valor de
> entre 0.05 y 0.2.
I S T dT
H v
=
( , )
.
.
0 1
2 5
Prctica 5.
Tratamiento de acelerogramas
An. VI - P5 - Pgina 12
siendo T el periodo del oscilador, S
v
el espectro de pseudovelocidad, y > el coeficiente de
amortiguamiento.
Figura 10 Grfico de Husid y definicin del concepto de duracin significante (Ds),
correspondiente al acelerograma registrado en Agrn (GR) el 24/02/97, mb = 4.3
Para la interpretacin de la intensidad espectral o intensidad de Housner, hay que
sealar que en el espectro de respuesta de cada periodo de vibracin se puede asociar a una
altura determinada de un edificio, por lo que podemos determinar que vibraciones pueden
hacer entrar en resonancia a una estructura determinada. Por lo tanto, la intensidad de
Housner se utiliza como un parmetro que mide el dao medio que sufren las estructuras.
Odograma
Grfica que describe la trayectoria de la partcula considerando que la seal
numerador define la coordenada en x, y la seal denominador define la coordenada en y.
3. METODOLOGA
1. Acceder a la base de datos de aceleracin del Instituto Geogrfico Nacional
(CDROM) y extraer el fichero ASCII con el registro del terremoto con epicentro en
Ventas de Huelma (GR) de m
B
= 4.3 obtenido en el acelergrafo de Agrn (GR),
el 24/02/97.
Anexo VI
Prcticas de Sismologa e Ingeniera Ssmica
An. VI - P5 - Pgina 13
2. Archivar la componente N-S en un fichero Ascii.
3. Crear un M-script que lea dicho fichero y lo escriba como una matriz en la que la
columna x correspondan a los tiempos y la columna y corresponda a la aceleracin
registrada.
4. Representar dicha matriz en Matlab y comprobar si el acelerograma necesita
correccin de lnea base. Para ello obtn el valor medio y la desviacin estndar de
la columna correspondiente a la aceleracin y comprueba que dicho valor sea
prcticamente cero. En caso contrario necesitaremos efectuar dicha correccin bien
en Matlab ( Ejecutar DETREND (y) donde y es la columna correspondiente a la
aceleracin o bien en DEGTRA (Pulsando el botn )).
5. Utilizar la herramienta filtrar de DEGTRA para llevar a cabo un filtrado paso-
banda con frecuencia mnima 0.150 Hz y frecuencia mxima 20 Hz. Archivar el
resultado obtenido y compararlo con los que se obtendran al aplicar los filtros
restantes.
6. Obtn a partir del acelerograma una vez realizado el filtrado paso-banda y por
integracin los registros correspondientes a velocidad y a desplazamiento.
Archiva los datos correspondientes a cada registro y obtn en Matlab o en
DEGTRA, un grfico con los tres registros.
7. Obtn el espectro de Fourier de la aceleracin, la velocidad y el desplazamiento
para diferentes valores de factor de suavizado y compara los resultados obtenidos.
Para el valor por defecto de suavizado, 6, representa en un diagrama trilogaritmico
dichos espectros de Fourier.
8. Obtn el pseudoespectro de respuesta de aceleracin, velocidad y
desplazamiento y representalos en un diagrama trilogartmico.
9. Obtener la duracin significante y el grfico de Husid .
10. Obtener la intensidad espectral o intensidad de Housner a partir de la ecuacin
propuesta en la teora.
11. Obtener el Odograma correspondiente al terremoto para las componentes N-S
y Vertical.
11. Repetir los pasos anteriores para las componentes E-W y vertical del acelerograma.
Prctica 5.
Tratamiento de acelerogramas
An. VI - P5 - Pgina 14
12. Escoger el registro obtenido en una estacin prxima cuyas caractersticas de suelo
sean muy diferentes a la del emplazamiento de Agrn. Obtener los espectros de
Fourier y de Respuesta y compararlos con los de la estacin de Agrn.
13. Realiza un informe con todos los resultados obtenidos y comenta las conclusiones
que se pueden extraer de cada uno de ellos.