Respuesta Sísmica en Sistemas No Lineales
Respuesta Sísmica en Sistemas No Lineales
RESPUESTA SÍSMICA
A SISTEMAS NO LINEALES
9.1 INTRODUCCION.
Se ha visto que para un sistema linealmente elástico la cortante basal pico inducida por el movimiento del suelo
es: Vb=(A/g)·w donde w es el peso del sistema, A, es la aceleración espectral correspondiente a un periodo de
vibración natural y un amortiguamiento determinado. Sin embargo la mayoría de los edificios son diseñados para
cortantes basales menores que la cortante basal asociada con un temblor fuerte que puede ocurrir en el sitio. A
partir de la Figura 9.1, donde el coeficiente de cortante basal, A/g, es graficada para el espectro de diseño
correspondiente a la aceleración pico del suelo de 0.4g, además es comparado con el coeficiente de cortante basal
especificado en el Código Uniforme de la Edificación de 1997 (UBC 97). Se observa una gran disparidad, que
implica que los edificios diseñados a partir de las fuerzas propuestas por el código se deformarán más allá del
límite elástico cuando estén sujetos a movimientos del suelo representados por el espectro de diseño para 0.4g.
1.2
1.0
Espectro de Diseño Elástico
Coeficiente de Cortante Basal
ügo = 0.4·g
0.8
0.6
0.4
Código Uniforme de la Edificación
R = 4 a 12
0.2
0.0 1 2 3
Figura 9.1 Comparación entre coeficientes de cortante basal a partir del espectro de diseño y UBC
De este modo no es de sorprenderse que los edificios sufran daños durante un movimiento intenso del suelo. El
reto para el ingeniero es de diseñar las estructura de tal forma que el daño sea controlado dentro un rango
Respuesta sísmica a sistemas no lineales 95
aceptable; obviamente el diseño fracasará si el sismo causa daños severos los cuales no pueden ser reparados, o si
se produce el colapso de la estructura.
De este modo la importancia central en ingeniería sísmica es comprender la respuesta de las estructuras
deformadas dentro el rango inelástico durante un movimiento intenso del suelo. Este capítulo trata sobre este
importante punto. Se hace una introducción al sistema elastoplástico y se describen los parámetros
correspondientes a dicho sistema, se presenta la ecuación de movimiento y se identifican varios parámetros que
describen el sistema y la excitación. Entonces se comparan la respuesta de sistemas elásticos e inelásticos con el
objeto de comprender como la fluencia influye en la respuesta estructural.
Los resultados experimentales indican que el comportamiento cíclico de la relación fuerza-deformación de una
estructura depende principalmente del material y del sistema estructural.
fs
Real
fy
Idealizado
u
uy um
Figura 9.2 Curva fuerza deformación durante la carga inicial: real e idealización elastoplástica.
Considerar la relación fuerza-deformación para una estructura durante su carga inicial mostrada en la Figura 9.2.
Resulta conveniente idealizar esta curva por una relación fuerza-deformación elastoplástica debido a que esta
aproximación permite desarrollar el espectro de respuesta en forma similar a un sistema linealmente elástico. La
aproximación elastoplástica a la curva real de fuerza-deformación esta representada en la Figura 9.2, de tal forma
que las áreas bajo las dos curvas son las mismas para un valor de desplazamiento máximo, um. En el proceso
inicial de carga este sistema idealizado es linealmente elástico con una rigidez k mientras la fuerza no exceda fy.
La fluencia comienza cuando la fuerza alcanza el valor de fy, esfuerzo de fluencia. La deformación en la cual la
fluencia comienza es uy, deformación de fluencia. En la fluencia la fuerza es constante (la rigidez es cero).
En la Figura 9.3 se muestra un típico ciclo de carga, descarga y recarga para un sistema elastoplástico, en el cual
se observa claramente que cuando el sistema alcanza el estado elastoplástico existen deformaciones permanentes
que se incrementan en cada ciclo.
Respuesta sísmica a sistemas no lineales 96
fs
fy
k
1
u
uy um
k k
1 1
-f y
fs
f0
Sistema Elastoplástico
fy
u
uy u0 um
Se desea evaluar la deformación pico de un sistema elastoplástico debido a un movimiento sísmico del suelo y
comparar esta deformación con la deformación pico causada por la misma excitación en el sistema lineal
correspondiente. Este sistema elástico esta definido de tal forma que tiene la misma rigidez del sistema
elastoplástico durante su fase de carga inicial; ver Figura 9.4. Ambos sistemas tienen la misma masa y
amortiguamiento. De este modo el periodo natural de vibración del sistema lineal correspondiente y del sistema
elastoplástico bajo oscilaciones pequeñas (uuy) es el mismo.
Respuesta sísmica a sistemas no lineales 97
fy uy
fy (9.1)
f0 u0
donde f0 y u0 son valores pico de fuerza resistente y deformación en el sistema lineal correspondiente, inducidos
por un sismo. Se puede interpretar f0 como la fuerza requerida para que la estructura permanezca con su límite
linealmente elástico durante un movimiento del suelo. Si el esfuerzo normalizado del sistema es menor que uno,
el sistema se deformará mas allá de su límite linealmente elástico, por ejemplo si fy = 0.5 implica que el esfuerzo
de fluencia del sistema es la mitad de la fuerza requerida para que el sistema permanezca elástico durante el
movimiento del suelo. El esfuerzo normalizado del sistema que no se deforma más allá de su límite linealmente
elástico es igual a la unidad, porque dicho sistema se puede interpretar como un sistema elastoplástico con fy=f0.
Alternativamente fy puede relacionarse con f0 a través de el factor de reducción de fluencia definido por:
f 0 u0
Ry (9.2)
f y uy
obviamente Ry es el reciproco de f y ; Ry es igual a la unidad para sistemas lineales y es mayor que uno para
sistemas que se deforman en el rango inelástico.
El pico o máximo absoluto de deformación del sistema elastoplástico debido al movimiento del suelo es denotado
por um. Es significativo normalizar um relacionado con la deformación de fluencia del sistema de la siguiente
manera:
u
m (9.3)
uy
esta relación adimensional es llamada factor de ductilidad. Para sistemas que se deforman en el rango inelástico,
por definición, um excede a uy y el factor de ductilidad es mayor que la unidad. Para el sistema lineal
correspondiente el factor de ductilidad es uno si este sistema es interpretado como un sistema elastoplástico con
fy=f0. y sus relaciones pueden expresarse como:
um
fy (9.4)
u0 Ry
fs fs
1 1
u
uy 1
(a) (b)
donde la fuerza resistente f s (u , u ) para un sistema elastoplástico es mostrado en la Figura 9.3. La ecuación 9.5 es
resuelta por un procedimiento numérico1. Para identificar los parámetros del sistema que tienen influencia en la
respuesta de deformación, la ecuación anterior es dividida por m para obtener:
donde:
k c f s (u, u )
n f s (u, u ) (9.7)
m 2 m n fy
La cantidad n es la frecuencia natural del sistema inelástico vibrando dentro el rango linealmente elástico
(uuy); ésta también es la frecuencia natural del sistema lineal correspondiente. Similarmente, es la razón de
amortiguamiento del sistema basado en el amortiguamiento crítico 2·m·n del sistema inelástico vibrando en su
rango linealmente elástico; es también la razón de amortiguamiento del sistema lineal correspondiente. La
función f s (u, u) describe la relación fuerza-deformación en forma parcialmente adimensional como se muestra
en la Figura 9.5(a). La ecuación 9.6 indica que para üg(t), u(t) depende de tres parámetros del sistema: n, y uy.
La ecuación 9.6 es reescrita en términos de (t) u(t)/uy para identificar la influencia que tiene el factor de
ductilidad ; de la ecuación 9.3 sustituyendo u(t)=uy·(t), u (t ) u y (t ) y u(t ) u y
(t ) en la ecuación 9.6 y
dividiendo por uy se tiene:
ug (t )
2 n n 2 f s ( , ) n 2 (9.8)
ay
donde ay=fy/m puede ser interpretada como la aceleración de la masa para producir la fuerza de fluencia y
f s ( , ) es mostrada en la Figura 9.5(b). La relación de aceleraciones üg(t)/ay es la razón entre la aceleración del
suelo y la medida del esfuerzo de fluencia de la estructura.
Para üg(t) dada, el factor de ductilidad depende de tres parámetros del sistema: n, y f y
1
Anil K. Chopra, pp 249 [ref. 12]
Respuesta sísmica a sistemas no lineales 99
Para entender como la respuesta de un sistema SDF es afectada por la acción inelástica o la fluencia a
continuación se compara la respuesta de un sistema elastoplástico con su sistema lineal correspondiente a través
de las Figuras 9.6 y 9.7.
Deformación u 4
-4 u0=3.34 in
2
f s / w = -ü' / g
-2 f 0 / w=1.37
0 10 20 30
Tiempo, [s]
Figura 9.6 Respuesta de un sistema lineal con Tn=0.5 [s] y =0 para el movimiento del suelo de El Centro
A diferencia de un sistema estático, el sistema inelástico después de que empieza afluir no oscila alrededor de su
posición inicial de equilibrio. La fluencia provoca que el sistema se desplace de su posición inicial de equilibrio y
hace que éste oscile alrededor de una nueva posición de equilibrio hasta que éste sea afectado por otro ciclo de
fluencia. Por lo tanto una vez que el movimiento del suelo se ha detenido, el sistema entra en reposo en una
posición diferente de la inicial (existe deformación permanente). De este modo una estructura que ha sufrido
fluencia significativa durante un sismo puede no permanecer recta después de éste.
Como se puede ver intuitivamente, se espera que sistemas de bajo esfuerzo de fluencia, fluyen más
frecuentemente y por intervalos de tiempo mayores. Con mayor fluencia, la deformación permanente, up, de la
estructura tiende a incrementarse. Para valores de Tn y dados, la deformación pico, um, de los tres sistemas
elastoplásticos es menor que la deformación pico, u0, del sistema lineal correspondiente.
El factor de ductilidad para un sistema elastoplástico puede ser calculado usando la ecuación 9.4. Por ejemplo la
deformación pico de un sistema elastoplástico con fy = 0.25 y del sistema lineal correspondiente son: um=1.75
[in.] y u0=2.25 [in.] respectivamente. Sustituyendo los anteriores valores en la ecuación 9.4 da el factor de
ductilidad =(1.75/2.25)·(1/0.25)=3.11. Ésta es la demanda de ductilidad impuesta en el sistema elastoplástico
por el movimiento del suelo. Esto representa un requisito importante en el diseño del sistema en el sentido de que
la capacidad de ductilidad (habilidad de deformarse más allá del límite elástico) debe exceder a la demanda de
ductilidad.
Respuesta sísmica a sistemas no lineales 100
A continuación se examina como el periodo natural de vibración Tn tiene influencia en: (1) La demanda de
ductilidad en un sistema elastoplástico; (2) Los valores relativos de la deformación pico um y u0 del sistema
elastoplástico y del sistema lineal correspondiente, respectivamente; (3) Los valores relativos del esfuerzo de
fluencia fy del sistema elastoplástico y la fuerza pico f0 impuesta en el sistema elástico.
2
c
d
u, [in] 0
b e
a
g
-2 f
um=1.71 in (a)
0.3
f s / w = -ü' / g
b c f y / w=0.17
0 a g
d
- f y / w=-0.17
e f
-0.3
(b)
+ Fluencia b c
g
ad
- Fluencia e f
0 5 10
Tiempo, [s]
(c)
0.3
b c
fs/w
0 g a
d
f e
-0.3
-2 -1 0 1 2
Deformación, u, [in]
(d)
Figura 9.7 Respuesta de un sistema elastoplástico con Tn=0.5 [s] y =0, y fy = 0.125 para el movimiento del suelo de El Centro (a)
Deformación; (b) Fuerza resistente y aceleración; (c) Intervalos de tiempo de fluencia; (d) Relación fuerza-deformación
La Figura 9.9 es una grafica de um como una función de Tn para cuatro valores de fy = 1,0.5,0.25 y 0.125; u0 es el
mismo que um para fy = 1. En la Figura 9.10 esta graficado el factor de ductilidad versus Tn para los mismos
cuatro valores de fy; =1 si fy = 1. El histograma de respuesta presentado en la Figura 9.8 provee valores para
u0=2.25 [in.] y um=1.62, 1.75 y 2.07 [in.] para fy = 0.5,0.25 y 0.125 respectivamente. Dos de estos cuatro datos
son identificados en la Figura 9.9. La demanda de ductilidad para los tres sistemas elastoplásticos es 1.44, 3.11
y 7.36 respectivamente. Estos tres datos se identifican en la Figura 9.10, en estas graficas también se identifican
los valores de los periodos Ta, Tb, Tc, Td, Te y Tf, que definen la regiones espectrales.
Respuesta sísmica a sistemas no lineales 101
Figura 9.8 Respuesta de deformación y fluencia de cuatro sistemas para el movimiento del suelo de El Centro; Tn=0.5 [s] y =5%, y
f y 1, 0.5, 0.25, 0.125
Las Figuras 9.9 y 9.10 demuestran que para una excitación dada la demanda de ductilidad y la relación entre um y
u0 dependen de Tn y de fy. Para sistemas de periodos grandes (Tn>Tf) en la región sensitiva de desplazamiento
del espectro, la deformación um de un sistema elastoplástico es independiente de fy y es igual a u0. Cuando el
sistema es muy flexible (la masa permanece fija) la masa es estacionaria mientras el suelo se mueve; la masa
experimenta una deformación pico igual al desplazamiento pico del suelo independiente defy. De este modo
umu0ug0 y la ecuación 9.4 da 1 / f y o =Ry. Esto implica que para un dado el esfuerzo de fluencia de
diseño para un sistema elastoplástico con Tn>Tf es 1/ veces el esfuerzo requerido para que el sistema
permanezca elástico.
Para sistemas en la región sensitiva de aceleración, lo cual implica sistemas de periodos muy pequeños, la
demanda de ductilidad puede ser muy grande. Entonces los sistemas con periodos extremadamente pequeños
(Tn>Ta) deben ser diseñados para un esfuerzo de fluencia fy igual a f0 requerida para que el sistema permanezca
en el rango elástico, de otra forma la deformación inelástica y la demanda de ductilidad pueden ser excesivas.
Respuesta sísmica a sistemas no lineales 102
Regiones Espectrales
2 2
1 1
0.5 0.5
2.25/8.4=0.27
_
f y=0.125
um /u g 0 1.62/8.4=0.19
0.1 0.1
u0 /u g 0 o um /u g 0
_
f y=0.25
0.05 0.05
_
f y=0.5 _
0.01 uo /u go f y=1
0.01
0.005 0.005
T a=0.035
T b=0.125
T c= 0.5
T e= 10
T f = 15
T d= 3
0.001 0.001
0.02 0.05 0.1 0.5 1 5 10 50
Periodo natural de vibración T n, s
Figura 9.9 Deformación pico de sistemas elastoplásticos y sistema lineal correspondiente debido al movimiento de El Centro; Tn esta
variando; =5% y f y 1, 0.5, 0.25, 0.125 y ug0=8.4 in. Para un ug0=8.4 in.
Regiones Espectrales
Aceleración Velocidad Desplazamiento
Sensitiva Sensitiva Sensitiva
100
50
c
e
d
f
T
T
T
T
Demanda de ductilidad,
20
_
10 f y = 0.125
7.36 _
5 f y = 0.25
_
3.11 f y = 0.5
a
2
T
_
T
1.44 fy=1
1
Figura 9.10 Demanda de ductilidad para sistema elastoplásticos debido al movimiento de El Centro; Tn esta variando; =5% y
f y 1, 0.5, 0.25, 0.125
Respuesta sísmica a sistemas no lineales 103
Para propósitos de diseño se desea determinar el esfuerzo de fluencia fy (o deformación de fluencia uy) del
sistema, necesario para limitar la demanda de ductilidad, impuesta por el movimiento del suelo, a un valor
especifico.
9.6.1 Definiciones
El espectro de respuesta es trazado par las cantidades:
Dy u y Vy n u y Ay n2 u y (9.9)
notar que Dy es la deformación de fluencia uy del sistema elastoplástico, no es igual a la deformación pico. La
grafica de Dy versus Tn para valores fijos de ductilidad es llamada espectro de respuesta de fluencia-
deformación. Análogamente al capítulo anterior similares graficas de Ay y Vy son llamadas espectro de respuesta
de seudoaceleración y espectro de respuesta de seudovelocidad respectivamente.
Las cantidades de Dy, Vy, y Ay pueden ser presentadas en una grafica tetralogarítmica en la misma forma que para
un sistema lineal. Esto es posible por que estas cantidades están relacionadas a través de:
Ay
Vy n Dy (9.10)
n
Ay
fy w (9.11)
g
donde w es el peso del sistema.
El esfuerzo de fluencia fy de un sistema elastoplástico para un factor de ductilidad específico se puede obtener
utilizando fy y la ecuación 9.1. Para tener mayor precisión en este calculo, fy es obtenido a partir de un
procedimiento iterativo. A partir de pares de datos disponibles (fy, ) se asume una relación lineal entre log(fy )
y log (), y a través de una interpolación se obtienefy correspondiente a un especifico. Se calcula el histograma
de respuesta del sistema con estefy para determinar el factor de ductilidad. Si este factor es suficientemente
cercano, con un error del 1 %, al especificado; el valor defy se considera satisfactorio, de otra forma es
modificado hasta que lo sea.
4. Determinar la respuesta u(t) del sistema lineal con los valores de Tn y seleccionados. A partir de u(t)
determinar la deformación pico u0 y la fuerza pico f0=k·u0. Estos resultados se muestran para Tn=0.5 en
la Figura 9.8(a)
5. Determinar la respuesta u(t) de un sistema elastoplástico con los mismos valores de Tn y y la fuerza de
fluencia fy =fy · f0, con un valor defy < 1 seleccionado. A partir de u(t) determinar la deformación pico
um y el factor de ductilidad asociado a partir de la ecuación 9.4. Repetir dicho análisis para valores defy
suficientes para desarrollar una serie de puntos (fy,) que cubran el rango de ductilidad de interés. Estos
resultados se muestran para fy = 0.5, 0.25 y 0.125
6. a. Para una seleccionada determinarfy a partir de los resultados del paso 5 usando el procedimiento
de interpolación descrito en la sección 9.6.2. Si más de un valor defy corresponde a un valor particular
de se elegirá el mayor.
b. Determinar la correspondiente ordenada espectral para el valor defy calculado en el paso 6(a). La
ecuación 9.1 da uy a partir del cual Dy, Vy y Ay pueden ser calculados utilizando la ecuación 9.9. Estos
datos proveen un punto en el espectro de respuesta graficado en las Figuras 9.11 y 9.12
7. Repetir los pasos del 3 al 6 para un rango de valores de Tn validos para el valor de seleccionado en el
paso 6(a)
En las Figuras 9.11 y 9.12 se presenta el espectro de respuesta construido por este procedimiento para un sistema
elastoplástico con =5% para el movimiento del suelo de El Centro para =1,1.5,2,4 y 8.
0.8
f y /w = A y /g
0.6
=1
0.4
0.2
0
0 0.5 1 1.5 2 2.5 3
Figura 9.11 Espectro de respuesta para un sistema elastoplástico para el movimiento de El Centro: =1, 1.5, 2, 4 y 8; =5%
Respuesta sísmica a sistemas no lineales 105
50
10
10
0
20
10
1
10
V , [in/s]
5
1
0.
1
2
01
0.
0.
1
1
in
0.5 y,
A
1
y
u 00
·g
0.
0.
01
0.2
0.02 0.05 0.1 0.2 0.5 1 2 5 10 20 50
Figura 9.12 Espectro de respuesta para un sistema elastoplástico para el movimiento de El Centro: =1, 1.5, 2, 4 y 8; =5%
Considerar un sistema SDF a ser diseñado para una ductilidad, , admisible, basado en una deformación
admisible y en una capacidad de ductilidad que se pueden alcanzar para los materiales y materiales de diseño
seleccionados.
Se desea determinar el esfuerzo de fluencia de diseño y la deformación de diseño para el sistema sujeto a una
excitación dada. El valor Ay/g se lee del espectro de respuesta para el valor admisible de y los valores conocidos
de Tn y . La ecuación 9.11 proporciona el esfuerzo de fluencia, fy, necesario para limitar la demanda de
ductilidad a la ductilidad admisible. La deformación pico es:
um u y (9.12)
donde: uy=fy/k=Ay/(n)2
El factor de ductilidad y la deformación pico um representan los requisitos de diseño asociados con la fuerza de
diseño fy. De este modo el ingeniero deberá diseñar y detallar la estructura de acuerdo a la capacidad de ductilidad
y la capacidad de deformación que ésta posee.
El esfuerzo de fluencia de diseño fy que permite a un sistema SDF tener una deformación en el rango inelástico
es menor que el esfuerzo requerido por la estructura para permanecer en el rango elástico. El esfuerzo de fluencia
de diseño se reduce con el incremento del factor de ductilidad, esta aseveración es mostrada con mayor claridad
en la Figura 9.13, que no es otra cosa que las Figuras 9.11 y 9.12 graficadas en forma diferente.
Respuesta sísmica a sistemas no lineales 106
La implicación practica de estas observaciones es que la estructura puede ser diseñada para ser sismorresistente
haciéndola fuerte o dúctil; o diseñándola económicamente combinando ambas propiedades. Considerar de nuevo
un sistema SDF con Tn=0.5 [s] y =5% a ser diseñado para el movimiento de El Centro. Si este sistema es
diseñado para una fuerza f0=0.919·w o mayor, permanecerá dentro el rango linealmente elástico durante esta
excitación; de este modo no necesita ser dúctil. Por otro lado si ésta puede desarrollar un factor de ductilidad de
8, solo necesita ser diseñada para 12% de la fuerza f0 requerida para un comportamiento elástico.
Alternativamente puede ser diseñada para una fuerza igual al 37% de f0 y una capacidad de ductilidad de 2. Para
algunos tipos de materiales y miembros estructurales la ductilidad es difícil de alcanzar; para otras el proveerles
ductilidad es mucho más fácil que proveerles resistencia lateral y el diseño práctico refleja esto.
1 1
=1
= 1.5
0.5 2
=2
_
fy =4 Ry
0.2 5
=8
0.1 10
a
f
T
T
T
T
0.05
0.02 0.05 0.1 0.2 0.5 1 2 5 10 20 50
La resistencia normalizada para un factor de ductilidad especifico depende de la relación del amortiguamiento ,
pero esta dependencia no es fuerte; es usualmente ignorada en aplicaciones de diseño.