Vargas, Pujades, Barbat, Hurtado, 2013

Documento descargado de http://http://www.elsevierciencia.es el 20/06/2013.
Copia para uso personal, se prohbe la transmisin de este documento por cualquier medio o formato.
Rev. int. mtodos numr. clc. diseo ing. 2013;29(2):6378
Revista Internacional de Mtodos Numricos para

Clculo y Diseo en Ingeniera
www.elsevier.es/rimni
ssmico de edicios
Evaluacin probabilista de la capacidad, fragilidad y dano
de hormign armado
Y.F. Vargas a, , L.G. Pujades a , A.H. Barbat a y J.E. Hurtado b
a
b
Universidad Politcnica de Catalu

na, Campus Norte, Jordi Girona 1-3, Edicio D2, 08034 Barcelona, Espa
na
Universidad Nacional de Colombia, Sede Manizales, Apartado 127, Manizales, Colombia
informacin del artculo
r e s u m e n
Historia del artculo:

Recibido el 27 de marzo de 2011
Aceptado el 28 de octubre de 2011
On-line el 23 de mayo de 2013
Actualmente en las zonas ssmicas existen estructuras altamente vulnerables, puesto que han sido cons sismorresistente o siguiendo normas obsoletas.
truidas sin seguir las especicaciones de normas de diseno
Muchas veces los mtodos para evaluar la vulnerabilidad de las estructuras no tienen en cuenta que su
comportamiento ssmico es dinmico y fuertemente no lineal y que, adems, las caractersticas estructurales y de la accin tienen grandes incertidumbres. En este artculo se propone evaluar la vulnerabilidad
de las estructuras, considerando que las propiedades mecnicas de los materiales y la accin ssmica son
variables aleatorias, mediante tcnicas avanzadas con base en el mtodo de Monte Carlo y en la dinmica
estocstica no lineal. Los resultados obtenidos con estas tcnicas se comparan con los correspondientes a una evaluacin de vulnerabilidad estndar, que utiliza modelos deterministas, a n de resaltar las
diferencias entre los 2 enfoques. La principal conclusin de este trabajo es que existe una necesidad de
abordar el problema de la evaluacin de la vulnerabilidad de estructuras desde una perspectiva probabilista que, al ser combinada con tcnicas avanzadas de anlisis del comportamiento estructural no lineal
esttico y dinmico, proporciona una poderosa herramienta que permite obtener informacin imposible de conseguir mediante modelos deterministas. Se incluyen resultados detallados obtenidos sobre un
edicio con forjados reticulares, que es una tipologa estructural muy utilizada en Espana.
2011 CIMNE (Universitat Politcnica de Catalunya). Publicado por Elsevier Espaa, S.L. Todos los
derechos reservados.
Palabras clave:
Fragilidad
ssmico
Dano
Incertidumbres
Edicios de hormign armado
Probabilistic assessment of the capacity, fragility and seismic damage

of reinforced concrete buildings
a b s t r a c t
Keywords:
Fragility
Seismic damage
Uncertainties
Reinforced concrete buildings
Currently, many structures existing in seismic areas are highly vulnerable because they have been built
without the use of seismic design codes or by using outdated codes. Often, methods for assessing the
vulnerability of the structures do not take into account that their seismic behavior is dynamic and highly
nonlinear and, moreover, that the structural characteristics and action have large uncertainties. This article aims to assess the vulnerability of structures taking into account that the mechanical properties of
materials and the seismic action are random variables, by using advanced techniques based on the Monte
Carlo method and on the nonlinear stochastic dynamics. The results obtained with these techniques
are compared with those corresponding to a standard vulnerability assessment, based on deterministic
models, in order to highlight the differences between both approaches. The main conclusion of this work
is the need to address the vulnerability assessment problem from a probabilistic perspective which, combined with advanced nonlinear static and dynamic structural analysis techniques, provides a powerful
tool giving information impossible to be captured by means of deterministic models. Finally, detailed
results obtained for a building with wafe slabs, which is a structural typology widely used in Spain, are
included and discussed.
2011 CIMNE (Universitat Politcnica de Catalunya). Published by Elsevier Espaa, S.L.
All rights reserved.
Autor para correspondencia.

Correos electrnicos: yeudy.felipe.vargas@upc.edu (Y.F. Vargas), lluis.pujades@upc.edu (L.G. Pujades), alex.barbat@upc.edu (A.H. Barbat),
jehurtadog@unal.edu.co (J.E. Hurtado).
0213-1315/$ see front matter 2011 CIMNE (Universitat Politcnica de Catalunya). Publicado por Elsevier Espaa, S.L. Todos los derechos reservados.
http://dx.doi.org/10.1016/j.rimni.2013.04.003
Documento descargado de http://http://www.elsevierciencia.es el 20/06/2013. Copia para uso personal, se prohbe la transmisin de este documento por cualquier medio o formato.
64
Y.F. Vargas et al / Rev. int. mtodos numr. clc. diseo ing. 2013;29(2):6378
1. Introduccin
Muchas de las variables involucradas en la evaluacin de la
vulnerabilidad y del riesgo ssmico de estructuras son altamente
aleatorias. La variabilidad de las caractersticas estructurales y de
la accin ssmica tiene especial relevancia, pero tambin la aleatoriedad de la respuesta del sistema, especialmente cuando el
comportamiento de este ltimo es no lineal. En lo referente a la
estructura, son aleatorias las variables relacionadas con las caractersticas de los materiales, la geometra y la dimensin de las
secciones, entre otras. Para cubrir esta variabilidad, las normas
sismorresistente recomiendan que se efecte un clculo
de diseno
determinista, utilizando valores disminuidos para la resistencia de
los materiales y valores aumentados para las acciones, mediante
coecientes de seguridad. De esta manera, se reconoce implcitamente la aleatoriedad tanto de las caractersticas mecnicas de los
materiales como de las propiedades de las acciones. Una desventaja
de este mtodo es que, cuando la respuesta de la estructura es no
lineal, el nivel de conanza asociado a la respuesta puede y suele
ser diferente al de las variables de entrada.
Por otra parte, existen en la actualidad estructuras en las zonas
ssmicas que no se disenaron

para resistir acciones dinmicas horizontales y, por tanto, fueron construidas sin seguir los requisitos de
sismorresistente e incluso, muchas veces, sin
una norma de diseno
control de los materiales. En el proyecto de otras muchas estruc que
turas se utilizaron normas que no cumplen criterios de diseno
hoy son generalmente aceptados. Por todo ello, las variables asociadas a la estructura y a su respuesta pueden tener incertidumbres
an mayores. Por supuesto, no debe olvidarse que si se proyecta
una estructura incluso con las normas actuales ms avanzadas y se
efectan controles de proyecto, de materiales y de obra, tendrn
que considerarse igualmente las incertidumbres en el proceso de
evaluacin de su vulnerabilidad.
Adems, todas las acciones, estticas y dinmicas, que se deben
considerar en el proyecto de las estructuras generan tambin
incertidumbres. En las normas, las acciones se suelen incrementar mediante coecientes que tratan de considerar, de una manera
aproximada, el efecto de dichas incertidumbres y aumentar as la
seguridad del proyecto. Sin embargo, el caso ssmico es muy especial, ya que la variabilidad de la accin es mucho mayor que la
esperada en las acciones estticas. Puede observarse tambin que
la accin ssmica denida por las normas ms antiguas puede diferir considerablemente de la que se prescribe en las normas actuales
aunque estas sean, a su vez, fuente de incertidumbre. Esta fuente de
aleatoriedad es, generalmente, la ms importante de todas, puesto
ssmica, como la
que todos los parmetros que describen la senal
aceleracin mxima en el suelo peak ground acceleration (PGA),
el contenido frecuencial del acelerograma y la duracin de la accin,
entre otros, son variables aleatorias. Partiendo de los espectros
de respuesta elstica correspondientes a la aceleracin ssmica de
buscan salvar esta aleatoriedad
una zona, las normas de diseno
proponiendo curvas suavizadas que representen los espectros de
aceleracin de la zona. En algunos pases los espectros se obtienen
a partir de estudios realizados en otros lugares, relacionando las
condiciones en las que se hizo el estudio inicial con las condiciones
propias de la zona. Esta forma de proceder puede ser desacertada,
como expone Bermdez [1], ya que si se comparan los espectros de
respuesta de los sismos que ocurren en la zona con el espectro que
proporciona la norma, los resultados pueden mostrar diferencias
signicativas. Es mejor prctica realizar estudios de microzonicacin para la zona estudiada y proponer los modelos adecuados [2].
Para conocer y prevenir el riesgo ssmico es necesario evaluar la
vulnerabilidad de las estructuras que existen actualmente en zonas
ssmicas. Con este propsito se han desarrollado diversos mtodos.
Uno de ellos es el mtodo del ndice de vulnerabilidad, en el cual
se dene la accin a partir de la EMS-98 mediante intensidades
macrossmicas y el comportamiento estructural mediante un ndice

de vulnerabilidad [3,4]. Otro mtodo muy difundido es el basado
en el espectro de capacidad, en el que la accin ssmica se dene
mediante espectros elsticos de respuesta y la vulnerabilidad o fragilidad del edicio mediante el espectro de capacidad; este ltimo
se calcula a partir de un anlisis esttico no lineal incremental (pushover analysis [PA]). Esta metodologa ha sido aplicada en varios
estudios anteriores [58]. Sin embargo, estos trabajos no tienen en
cuenta las incertidumbres en los parmetros que describen el problema, que deben tratarse como variables aleatorias. Por otra parte,
para evaluar el comportamiento de una estructura ante acciones
ssmicas fuertes resulta conveniente efectuar un anlisis dinmico
no lineal que, en muchos casos, se usa como referencia para comparar los resultados obtenidos con mtodos ms simplicados [911].
No obstante, en estos trabajos tampoco se trata el problema dinmico no lineal desde un enfoque probabilista.
El propsito de este artculo es evaluar la vulnerabilidad ssmica
de una tipologa estructural frecuentemente utilizada en Espana

como vivienda multifamiliar, concretamente la de los edicios de
hormign armado con forjados reticulares, que es conocida por su
ductilidad limitada [12,13]. En la evaluacin de su vulnerabilidad
se va a considerar la aleatoriedad de las variables que describen el
problema, lo que requiere el uso de tcnicas avanzadas basadas en
la simulacin por Monte Carlo y en el anlisis dinmico no lineal
de estructuras. Cabe destacar que cuando una estructura trabaja en
rgimen no lineal y se incorporan a su anlisis las incertidumbres
en sus caractersticas, la variabilidad de su respuesta es mucho
mayor que cuando el sistema es lineal. Tcnicas como las antes
mencionadas han sido usadas por Fragiadakis y Vamvatsikos [14]
para evaluar el comportamiento no lineal de estructuras mediante
el anlisis esttico no lineal, considerando las incertidumbres en las
propiedades mecnicas de los materiales y utilizando la simulacin
de datos por el mtodo de Monte Carlo. Por otra parte, Dolsek [15]
propone tener en cuenta la accin ssmica como variable aleatoria
a partir de acelerogramas reales en este tipo de estudios, pero no
las incertidumbres en las caractersticas estructurales. Siguiendo
el enfoque dado en estos trabajos, en este artculo se propone
evaluar la vulnerabilidad de la estructura considerando las incertidumbres tanto en las propiedades mecnicas de los materiales
como en la accin ssmica y tratndolas como variables aleatorias.
Como metodologa de evaluacin de la vulnerabilidad se usa la tcnica basada en el espectro de capacidad. Los resultados para un
edicio con forjados reticulares sern comparados con el enfoque
determinista clsico para analizar y cuanticar los niveles de incertidumbre cuando un fenmeno altamente aleatorio se reduce a uno
determinista.
2. Descripcin del edicio

La estructura analizada en este trabajo es de hormign armado,
para viviendas de tipo multifamiliar, tiene 8 niveles
tpica en Espana
y 6 vanos y est compuesta por forjados reticulares y pilares (g. 1).
En el modelado de este edicio los forjados reticulares se aproximan
por vigas de inercia equivalente. Adems, debido a su simetra en
planta, la estructura puede ser modelizada en 2D aproximndola
mediante un prtico virtual.
En el modelo de clculo, los elementos estructurales (vigas
equivalentes y pilares) siguen una ley de comportamiento elastoplstica en la que no se tienen en cuenta el endurecimiento ni
el ablandamiento. Las supercies de plasticacin estn denidas por el diagrama de interaccin momento-fuerza axial en los
pilares y por el de momento-deformacin angular en las vigas.
Las acciones aplicadas siguen las recomendaciones dadas por el
Eurocdigo 8 (EC8) [16] para estructuras de hormign armado.
El modelo estructural empleado es capaz de reproducir solo el
65
3.1. Anlisis esttico no lineal y curva de capacidad
Tabla 1
Caractersticas de las variables aleatorias de entrada
fc
fy
Valor medio (kPa)
Desviacin estndar (kPa)
Coeciente de variacin
30.000
420.000
1.500
21.000
0,05
0,05
fallo por momento ector y por cortante. Por ahora se ha usado

esta aproximacin para poder manifestar la gran incertidumbre
en la respuesta estructural. Las caractersticas mecnicas del hormign y del acero se denen a partir de los valores dados en
los planos. Como se ha mencionado en la introduccin, las nor proponen valores caractersticos de la resistencia
mas de diseno
de los materiales que se obtienen durante el control de calidad
tanto en su fabricacin como en la obra, a partir de ensayos de
compresin en muestras de hormign y de traccin en muestras
de acero. Mediante estos ensayos, la resistencia de los materiales se puede representar como una variable aleatoria. En el
caso de proyecto, el valor caracterstico de la resistencia de los
materiales garantiza que la probabilidad de ocurrencia, de una
muestra de resistencia menor a este, es del 5%. Sin embargo, las
condiciones en las que se hace el ensayo son muy diferentes a
las condiciones reales de la construccin, y por tanto las dispersiones en la resistencia a compresin del hormign (fc) y en la
resistencia de uencia del acero (fy) pueden ser mayores. La tabla
1 muestra el valor medio, la desviacin estndar y el coeciente de
variacin de las variables aleatorias consideradas, en el supuesto
de que sigan una distribucin normal.
El coeciente de variacin es la relacin que existe entre la desviacin estndar y el valor medio de la variable aleatoria. A partir de
los valores de la tabla 1, y suponiendo una funcin de distribucin
gaussiana de fc y fy, se calculan los valores caractersticos de resistencia de los materiales, que para el hormign ser fc = 27.525 kPa
y para el acero ser fy = 385.350 kPa.
3. Mtodo de evaluacin de la vulnerabilidad y del riesgo
ssmico basndose en el espectro de capacidad
Se describen a continuacin los aspectos metodolgicos de la
evaluacin de la vulnerabilidad de las estructuras de hormign
armado. Para mayor claridad se utiliza como ejemplo ilustrativo
la estructura descrita en el apartado anterior, y se supone que las
propiedades mecnicas de los materiales son las correspondientes a los valores caractersticos dados previamente. Es decir, en
este apartado se estudia la vulnerabilidad desde una perspectiva
determinista.
3.2. Espectro de capacidad

La resistencia de un edicio puede ser representada a partir de
la curva de capacidad que se muestra en la gura 2, que puede
ser transformada en un espectro de capacidad, el cual relaciona el
2000
Curva de capacidad determinista
1800
1600
Cortante en la base (kN)
Figura 1. Edicio estudiado, con H = 25,65 m, B = 24 m y perodo fundamental de

1,444 s.
Una herramienta frecuentemente usada para evaluar el comportamiento de las estructuras frente a acciones horizontales es
el anlisis esttico no lineal de carga incremental, comnmente
llamado pushover analysis (PA). Esta herramienta numrica consiste en aplicar una carga horizontal a la estructura de acuerdo
con un patrn de fuerzas determinado e ir incrementando su
valor hasta que se considere que la estructura ha colapsado. Utilizando este tipo de anlisis, se considera una variable de control
que, en este trabajo, es el desplazamiento en el techo de la
estructura, , que depender de la fuerzas aplicadas sobre esta
y cuya suma es igual al cortante en la base del edicio V. La
relacin entre las variables V y es la curva de capacidad. Como
se ha mostrado en otros trabajos, el patrn de carga es decir,
la forma que tenga la carga en funcin de la altura (triangular,
rectangular u otros patrones) tiene una inuencia relevante
sobre los resultados [9]. Un criterio muy usado consiste en
considerar que la forma de la carga es proporcional al primer
modo de vibracin. Sin embargo, cuando se realiza un clculo
no lineal incremental y el modelo estructural empieza a plasticar, es necesario recalcular el patrn de acciones debido
a que, por la prdida de rigidez de la estructura, la forma asociada
a la frecuencia del primer modo de vibracin cambia. En este
trabajo, debido al enfoque probabilista, el PA ser ejecutado repetidamente, con lo que conviene aplicar un procedimiento que dena
automticamente el patrn de carga y el lmite de carga horizontal.
Para resolver estos problemas, Satyarno [17] propone el anlisis
no lineal incremental adaptativo que establece el lmite de carga
horizontal como una funcin del decremento de la frecuencia
fundamental tangente, es decir, de la frecuencia asociada al primer
modo de vibracin, que se calcula en cada incremento de carga.
Por tanto, en cada paso se calcula el primer modo para determinar
la forma de la carga en la altura. Una descripcin detallada de
este procedimiento se encuentra en los manuales del programa
Ruaumoko [18], en el que se encuentra implementado el procedimiento de Satyarno. Realizando un anlisis esttico no lineal de la
estructura [19,20], se obtiene la curva de capacidad que se muestra
en la gura 2.
1400
1200
1000
800
600
400
200
0
0
0.05
0.1
0.15
0.2
0.25
0.3
0.35
0.4
Desplazamiento (m)
Figura 2. Curva de capacidad obtenida a partir del enfoque determinista.
66
Tabla 2
Distribucin de probabilidades para los dsi
Aceleracion Espectral (g)
0.3
P(ds1 ) = 0, 5
P(ds2 ) = 0, 5
P(ds3 ) = 0, 5
P(ds4 ) = 0, 5
0.25
0.2
(Dy,Ay)
P(ds3 )
P(ds4 )
0,119
0,5
0,866
0,988
0,012
0,135
0,5
0,881
0
0,008
0,104
0,5
espectrales correspondientes a los umbrales de los 4 estados de

dsi :
dano
0.1
ds1 = 0, 7Dy
0.05
ds2 = Dy
Espectro de capacidad
Representacin bilineal
0.05
0.1
0.15
0.2
0.25
0.3
(3)
ds3 = Dy + 0, 25(Du Dy)

0.35
Desplazamiento espectral (m)
ds4 = Du
De estas ecuaciones se obtiene que para la estructura modelada
son: ds1 = 8,8 cm, ds2 = 12,7 cm, ds3 = 14,8 cm y
los estados de dano
ds4 = 21,2 cm.
Figura 3. Espectro de capacidad y la representacin bilineal.
desplazamiento espectral con la aceleracin espectral mediante las

siguientes ecuaciones [21]:
sdi = i /PF1
(1)
Vi /W
1
(2)
sai =
P(ds2 )
0,5
0,896
0,992
1
(Du,Au)
0.15
P(ds1 )
El subndice i de las ecuaciones (1) y (2) hace referencia a los

incrementos de carga aplicados sobre la estructura durante el anlisis no lineal incremental. sdi es el desplazamiento espectral, i es
el desplazamiento en el techo de la estructura, PF1 es el factor de
participacin modal del primer modo de vibracin, sai es la aceleracin espectral, Vi es el cortante en la base, W es el peso de la
estructura y 1 es el coeciente de masa modal del primer modo
de vibracin.
Reemplazando en las ecuaciones (1) y (2) los valores de la curva
de capacidad mostrados en la gura 2, se obtiene el espectro de
capacidad de la gura 3. Dicho espectro puede ser representado en
un formato bilineal, lo que es til para denir los estados de dano.

Las hiptesis para construir el espectro de capacidad bilineal son:
1) El rea bajo la curva bilineal debe ser igual al rea de la curva
original. 2) Las coordenadas del punto de mximo desplazamiento
deben coincidir en las 2 curvas. 3) La pendiente del tramo inicial
debe ser igual en las 2 curvas. En la gura 3 se muestra tambin
la representacin bilineal del espectro de capacidad, que se puede
denir completamente mediante los puntos (Dy, Ay) y (Du, Au) mostrados en la gura 3. Estos puntos corresponden a la uencia y al
colapso de la estructura, respectivamente.
3.4. Curvas de fragilidad

Las curvas de fragilidad representan la probabilidad de que
en funcin del parmetro
se alcance o exceda un estado de dano
que representa la accin ssmica, que en este trabajo ser el
desplazamiento espectral. Para construir las curvas de fragilidad
se consideran las siguientes hiptesis: 1) En los desplazamientos
dsi , la probabiespectrales de los umbrales de cada estado de dano
es del 50%. 2) Las
lidad de que se iguale o exceda el estado de dano
curvas de fragilidad siguen una funcin de probabilidad lognormal
cumulativa descrita por la siguiente ecuacin:
P dsi /sd =
1
Ln
dsi
sd
(4)
dsi
donde sd es el desplazamiento espectral y dsi es la desviacin

ssmico
estndar del logaritmo natural de la variable dsi . 3) El dano
esperado en los edicios sigue una distribucin de probabilidad
binomial. Por tanto, a partir de las hiptesis 1 y 3 se obtiene la
distribucin de probabilidades para los dsi que se muestra en
la tabla 2.
Una vez calculadas estas probabilidades, la hiptesis 2 permite
calcular las curvas de fragilidad por medio de una tcnica de mnimos cuadrados. En la gura 4 se muestran las curvas de fragilidad
1
P(Nulo)
0.9
P(Leve)
0.8
3.3. Estados de da
no
0.7
P[ds/sd]
P(Moderado)
de la
En diferentes trabajos se ha propuesto calcular el dano
ds,
estructura partiendo de la denicin de unos estados de dano
en la estructura para un desque son una descripcin del dano
plazamiento espectral dado. Por ejemplo, en HAZUS 99 [22] y en
RISK UE [23] se denen 4 ds: leve, moderado, extensivo y completo.
Su descripcin depende de la tipologa estructural. Por ejemplo,
en el caso de estructuras de hormign armado, el ds leve se describe como el inicio de suras por momento o por cortante en las
vigas y columnas; en el ds completo se considera que la estructura
colapsa o que hay un inminente riesgo de colapso (HAZUS 99) [22].
RISK UE [23] propone denir los ds, de forma simplicada, a partir
del espectro de capacidad en su representacin bilineal. Partiendo
de los valores (Dy,Ay) y (Du,Au), se denen los desplazamientos
0.6
0.5
0.4
0.3
CF1
CF2
CF3
CF4
P(Extensivo)
0.2
0.1
P(Completo)
0
0
0.05
0.1
0.15
0.2
0.25
0.3
0.35

Figura 4. Curvas de fragilidad para el edicio estudiado.
0.4
1
67
3
Espectro de respuesta 1A del EC8
0.9
2.5
0.8
0.7
sa (g)
ID
0.6
0.5
1.5
0.4
1
0.3
0.2
0.5
0.1
0
ndice de dao
0.05
0.1
0.15
0.2
0.25
0.3
0.35
0.4
0
0
1.5
2.5
3.5
1.5
Espectro tipo 1 suelo A formato sasd
para los dsi correspondientes al edicio calculado. En esta misma

gura se muestra tambin cmo a partir de las curvas de fragilidad
se pueden calcular las probabilidades de ocurrencia de cada estado
de dano.
sa (g)
3.5. ndice de da
no esperado
A partir de las probabilidades de ocurrencia de cada estado de
P(dsi ) que, a su vez, se obtienen fcilmente a partir de las curdano
esperado (ID),
vas de fragilidad, se puede calcular el ndice de dano
medio normalizado y puede interpretarse
que es el estado de dano
global esperado en la estructura.
como una medida del dano
1
iP(dsi )
n
Periodo (s)
esperado calculado con el enfoque determinista.

Figura 5. ndice de dano
0.5
ID =
0.5
0
0
(5)
i=0
donde n es el nmero de estados de danos

considerados, en este
caso 4, y P(dsi ) es la probabilidad de que ocurra un dsi . En la gura 5
se puede apreciar el ID calculado con el enfoque determinista, a
partir de las curvas de fragilidad de la gura 4 y los desplazamientos
espectrales comprendidos entre 0 y 0,4 m. La curva de la gura 5
puede interpretarse como una curva de vulnerabilidad.
3.6. Demanda ssmica
En la metodologa basada en el espectro de capacidad es necesario calcular la demanda ssmica en trminos del desplazamiento
espectral sd. Existen diversas metodologas para evaluar la amenaza ssmica de una zona, desde mtodos basados en la intensidad
macrossmica, pasando por espectros de respuesta propuestos en
hasta la caracterizacin por medio de acenormativas de diseno,
lerogramas reales. En este trabajo, la accin ssmica se dene en
trminos del espectro de respuesta elstica en aceleracin con 5%
de amortiguamiento. Los resultados obtenidos utilizando esta denicin sern comparados con los del anlisis dinmico no lineal, en
el cual la accin ssmica es evaluada por medio de acelerogramas
reales tomados de bases de datos espanolas

y europeas. Para ello se
seleccionan acelerogramas, de forma que los espectros elsticos de
respuesta de los registros sean compatibles con el espectro usado
en la metodologa basada en el espectro de capacidad. Caracterizar
la accin ssmica desde el punto de vista de los espectros de respuesta resulta ser una buena alternativa para evaluar el riesgo de
la zona ya que, como se ver ms adelante, es posible calcular el
desplazamiento espectral mediante procedimientos simplicados.
0.05
0.1
0.15
0.2
0.25
sd (m)
Figura 6. a) Espectro normalizado tipo 1, suelo A (tomado del EC8). b) Espectro
tipo 1, suelo A, en formato sa-sd escalado a 0,5 g y espectro de capacidad determinista
del edicio objeto de estudio.
Actualmente, la mayora de los pases proponen, en sus normas

sismorresistente, espectros de diseno
en el formato de
de diseno
espectros de respuesta con 5% de amortiguamiento. Por ejemplo,
tiene la NCSE-02 y Europa el EC8, que es aplicable a los paEspana
ses de la Unin Europea. En este artculo se trabaja con los espectros
del EC8, en el que se proponen 2 tipos de espectros que se diferencian por la magnitud de la onda supercial Ms . Concretamente, el
EC8 recomienda usar los espectros tipo 1 si la Ms es mayor a 5,5; de
lo contrario, se usan los espectros tipo 2. Este hecho es interesante,
ya que permite calicar de una mejor manera las zonas con peligrosidad ssmica entre moderada y baja. En la gura 6a se muestra
el espectro tipo 1 para suelo A del EC8. Por otra parte, los espectros
de respuesta relacionan la aceleracin espectral sa con el perodo T
para un sistema con un grado de libertad. Mahaney et al. [24] proponen transformar el espectro de respuesta desde el formato sa-T a
un formato sa-sd por medio de la ecuacin (6). Este formato permite
comparar directamente la demanda y la capacidad de la estructura.
En la gura 6b se muestra el espectro tipo 1 suelo A del EC8 despus de aplicar esta transformacin y tambin la curva de capacidad
determinista del edicio obtenida con los valores caractersticos de
la resistencia de los materiales
sdi =
sai Ti2
42
(6)
68
a 1.5
Espectro de demanda
sa (g)
0.5
Punto de desempeo
0.05
0.1
0.15
0.2
0.25
sd (m)
b 1.5
sa (g)
En este trabajo se usa el programa Ruaumoko [18], que permite

calcular tanto la respuesta esttica no lineal como la dinmica no
lineal para estructuras de hormign armado y de acero, entre otros
materiales. El modelo de amortiguamiento usado en este artculo,
cuando se lleva a cabo un ADNL, es del tipo Rayleigh y se considera un valor del 5% del amortiguamiento crtico para los modos de
vibracin 1 y 8 de la estructura.
El anlisis dinmico incremental (ADI) propuesto por Vamvatsikos y Cornell [25] es una herramienta usada para evaluar el
comportamiento de las estructuras frente a acciones ssmicas por
medio de acelerogramas. El objetivo del ADI es obtener una medida
de la estructura mediante el incremento de la intensidad
del dano
de la accin escalndola para diferentes niveles de PGA. Vamvatsikos y Cornell [25] hacen una interesante analoga entre el PA y
el ADI, ya que en ambos procedimientos se incrementa la carga
aplicada a la estructura y se mide la respuesta del sistema desde el
punto de vista de una variable de control que puede ser el desplazamiento de techo o la deriva mxima de piso, entre otras. Este procedimiento permite obtener la respuesta de una estructura a medida
que aumenta la accin ssmica. En este artculo se utiliza un procedimiento similar al ADI cuando se calcula la respuesta por medio
de los mtodos simplicados propuestos en el ATC-40 (vanse los
apartados 4.5 y 4.6); de esta manera se puede hacer una comparacin de resultados en un intervalo de PGA y no en un solo punto.
4. Clculo de la vulnerabilidad de la estructura con el
enfoque probabilista
0.5
Punto de capacidad por demanda
0
0
0.05
0.1
0.15
0.2
0.25
En el apartado anterior se ha hecho una descripcin de mto fsico esperado de

dos para evaluar la vulnerabilidad y el dano
las estructuras debido a la accin ssmica de una zona. Ahora se
propone realizar esta evaluacin teniendo en cuenta que las variables de entrada del problema son aleatorias, particularmente las
propiedades mecnicas de los materiales y la accin ssmica.
sd (m)
Figura 7. a) Mtodo ALE para calcular el sd. b) Mtodo PA-8 para calcular el sd.
3.7. Punto de capacidad por demanda

Una vez calculada la capacidad de la estructura y la demanda a
la que estar sometida, es posible calcular el punto de capacidad
por demanda que dene el desplazamiento espectral sd que experimenta el edicio cuando se somete a una accin ssmica. Este sd
puede estimarse mediante mtodos simplicados propuestos en
el ATC-40. En este trabajo se utilizan 2 de estos mtodos. El primero es la aproximacin lineal equivalente (ALE), que consiste en
prolongar el tramo elstico del espectro de capacidad y cruzarlo
con el espectro elstico de demanda en formato sa-sd. La proyeccin de este punto de cruce sobre la curva de capacidad ser el
punto de capacidad por demanda. El segundo es un mtodo iterativo en el que el espectro de demanda es reducido debido a
la ductilidad de la estructura que se calcula a partir de la forma
bilineal del espectro de capacidad. Este mtodo se denomina procedimiento A en el ATC-40 y se encuentra en el captulo 8, y por
tanto en este artculo ser llamado PA-8. La gura 7a,b muestra el
procedimiento para calcular el sd con la ALE y el PA-8. Si se escala
el espectro de demanda para un PGA de 0,5 g, con la ALE se obtiene
un sd = 0,1765 y con el PA-8 un sd = 0,1685 m. Los ID asociados a
estos valores son 0,7144 y 0,6789, respectivamente.
3.8. Anlisis dinmico no lineal incremental
El anlisis dinmico no lineal (ADNL) permite evaluar la respuesta de una estructura sometida a una historia de aceleraciones.
4.1. Anlisis esttico no lineal incremental mediante el mtodo

de Monte Carlo
En este trabajo se efecta el anlisis esttico no lineal incremental 10.000 veces, considerando como variables aleatorias la
resistencia a compresin del hormign, fc, y la resistencia a traccin del acero, fy. Los valores de las variables aleatorias pueden
verse en la tabla 2. Como se muestra en la gura 1, la estructura
est compuesta por 54 pilares, que se dividen en 7 grupos segn
las caractersticas de la seccin transversal. Estos pilares se unen a
forjados reticulares que son modelados como vigas que tienen una
inercia equivalente. Puesto que los forjados son iguales en todos los
niveles, solo habr un tipo de seccin transversal. Para cada seccin transversal se genera una muestra aleatoria de resistencia del
hormign y del acero segn el procedimiento expuesto en Kalos y
Whitlock [26]. Por tanto, en cada ejecucin del PA la resistencia de
los elementos estructurales es diferente. Estas variables generadas
siguen una distribucin de probabilidad gaussiana.
En la gura 8 se muestran las 10.000 curvas de capacidad obtenidas, la curva de capacidad media y la curva de capacidad obtenida
con los valores caractersticos en el anlisis determinista. El
valor de desplazamiento mximo determinista es excedido por
el 82,5% de los valores mximos de desplazamiento de todas las
curvas de capacidad, lo que quiere decir que la probabilidad de que
de colapso antes de este
la estructura alcance el estado de dano
valor es del 17,5%. Por otra parte, las curvas de capacidad calculadas
tienen caractersticas que se distribuyen de forma aleatoria, como
la rigidez inicial Ko o el desplazamiento ltimo du. En la gura 9a
estas caractersticas se relacionan con las variables de entrada fc
y fy, respectivamente. Se puede observar que el desplazamiento
69
Tabla 3
Matriz de correlacin entre las variables aleatorias del sistema
2000
1800
1600
fy
fc
q
Ko
dy
du
1400
1200
1000
600
Curva de capacidad media

Curva de capacidad determinista
Todas las curvas
200
0
0
0.05
0.1
0.15
0.2
0.25
0.3
0.35
0.4
Desplazamiento en el techo (m)

Figura 8. Las 10.000 curvas de capacidad obtenidas a partir del enfoque probabilista, la curva de capacidad media y la curva de capacidad obtenida con el enfoque
determinista.
ltimo de la estructura tiende a ser levemente mayor a medida que

la resistencia de plasticacin del acero, fy, aumenta. La gura 9b
muestra la dependencia entre la rigidez inicial de la estructura y la resistencia a compresin del hormign, fc. Para medir
Ko
dy
du
0,000
1,000
0,137
0,551
0,178
0,098
0,126
0,137
1,000
0,477
0,203
0,973
0,006
0,551
0,477
1,000
0,713
0,318
0,335
0,178
0,203
0,713
1,000
0,028
0,205
0,098
0,973
0,318
0,028
1,000
el grado de dependencia entre las variables de entrada y salida

se calcula la matriz de correlacin, que tiene por componentes el
coeciente de correlacin entre cada par de variables aleatorias.
Antes de calcular esta matriz se obtiene la representacin bilineal
para todas las curvas de capacidad, mediante el procedimiento
descrito en el apartado 3.2. De esta manera se obtienen el desplazamiento espectral de plasticacin, denominado dy, y el
desplazamiento ltimo de la estructura, denominado du. Ntese
que dy y du se calculan sobre la curva de capacidad y no sobre
el espectro de capacidad, y por tanto son diferentes a Dy y Du
mostrados en la gura 3. Otras variables aleatorias que se deben
calcular son el factor de ductilidad de la estructura q:
q=
du
dy
Correlacin entre fy y du
i,j =
0.45
du (m)
0.35
0.3
0.25
0.2
a)
0.1
3.9
4.1
4.2
4.3
4.4
Cov(xi , xj )
(8)
xi xj
donde ij dene la correlacin entre las variables aleatorias xi y xj ,

Cov es la covarianza y xj es la desviacin estndar de xj . La matriz
de correlacin entre todas las variables aleatorias consideradas se
muestra en la tabla 3.
En general, 2 variables aleatorias son independientes cuando el
correspondiente elemento ij de la matriz de correlacin tiende a
cero, aunque tambin existe la posibilidad de que no se relacionen linealmente. Destaca la correlacin entre la rigidez inicial y el
desplazamiento de plasticacin y entre la ductilidad y el desplazamiento de capacidad ltima. Al calcular la matriz de correlacin
se observa la riqueza de la informacin obtenida cuando se hace el
clculo del PA usando el mtodo de Monte Carlo.
0.4
0.15
(7)
y la rigidez del tramo elstico de la curva de capacidad, denominada

Ko. Una vez denidas las variables aleatorias de salida, se procede
a calcular la matriz de correlacin mediante la siguiente ecuacin:
0.5
4.5
x 105
fy (kPa)
8600
4.2. El espectro de capacidad bilineal como variable aleatoria

Las ecuaciones (1) y (2) permiten calcular los espectros de capacidad a partir de la curva de capacidad, transformando las variables
dy y du en Dy y Du. Aplicando las ecuaciones (3) se obtienen los
dsi como variables aleatorias. La
umbrales de los estados de dano
gura 10 muestra los resultados obtenidos y la tabla 4 muestra el
valor medio dsi y la desviacin estndar dsi de cada estado de
y los estados de dano
calculados con el enfoque determinista
dano
dsi,det .
Correlacin entre fc y Ko
8500
8400
8300
Ko (kPa)
fc
1,000
0,000
0,126
0,006
0,335
0,205
800
400
fy
8200
8100
Tabla 4
Valores medios y desviaciones estndar de los umbrales de los estados de dano

considerados como las variables aleatorias. La tercera la muestra los coecientes
de variacin y la cuarta los umbrales correspondientes al espectro de capacidad
determinista
8000
7900
7800
7700
2.75
2.8
2.85
2.9
2.95
3.05
3.1
3.15
fc (kPa)
Figura 9. a) Correlacin entre du y fy. b) Correlacin entre Ko y fc.
3.2
x 104
ds (cm)
ds (cm)
cv
dsdet (cm)
ds1
ds2
ds3
ds4
9,2
0,16
0,0164
8,8
13,2
0,22
0,0164
12,7
15,6
0,51
0,0315
14,8
22,7
1,6
0,0706
21,2
70
0.35
1
0.9
0.3
0.7
0.2
P[ds/sd ]
Aceleracin espectral (g)
0.8
0.25
ds1
ds2
ds3
0.15
ds4
0.6
0.5
0.4
dsi
0.1
ds1, det
0.3
Curvas de fragilidad del ds1
ds2, det
0.05
0.2
ds3, det
ds4, det
0.1
0
0.05
0.1
0.15
0.2
0.25
0.3
0
0
0.35
0.1
0.2
0.3
0.4
0.5
0.6
obtenidos con el enfoque probabilista y el determinista.

Figura 10. Estados de dano
4.3. Curvas de fragilidad probabilistas

Una vez calculados los dsi como variables aleatorias, se calculan las curvas de fragilidad siguiendo las hiptesis descritas en
el apartado 3.4. La gura 11 muestra las curvas de fragilidad calculadas, donde se puede apreciar que la variabilidad es mucho
debido
mayor a medida que se incrementa el estado de dano,
probablemente a la no linealidad creciente en la respuesta de la
estructura.
Tal como se ha visto, cada curva de fragilidad queda completamente denida si se conocen dsi y dsi . En la gura 11 se puede
leve y
ver que las curvas de fragilidad asociadas al estado de dano
moderado muestran poca dispersin tanto en dsi como en dsi ; por
Tabla 5
obtenidos con el mtodo
Probabilidad de que los umbrales de los estados de dano
calculado con el mtodo determinista
probabilista sean menores al estado de dano
P[ds1 < ds1,det ]
P[ds2 < ds2,det ]
P[ds3 < ds3,det ]
P[ds4 < ds4,det ]
0,4%
0,4%
2%
17,5%
0.5
ds1
ds2
ds3
0.4
ds4
En la gura 10 se aprecia que los umbrales ds1 y ds2 estn fuertemente correlacionados con la aceleracin espectral porque los
desplazamientos espectrales asociados a estos estn muy cerca del
campo de comportamiento lineal de la estructura. Por otra parte,
los umbrales ds3 y ds4 muestran un comportamiento ms aleatorio.
Destaca la alta aleatoriedad asociada al umbral del estado de dano

de colapso. Se recuerda que en las simulaciones cuyos resultados
se muestran en la gura 10, solo se han tenido en cuenta como
variables aleatorias las propiedades mecnicas de los materiales
que constituyen la estructura y no las fuerzas aplicadas.
Es importante observar en la gura 10 que, a medida que el
dsi aumenta, los dsi,det cubren un menor
umbral del estado de dano
intervalo de los valores obtenidos con el enfoque probabilista, es
decir, aumenta la probabilidad de que se produzca un dsi menor
que dsi,det . Esto muestra que no es suciente enfocar un problema
no lineal desde una perspectiva determinista cuando las variables
de entrada son aleatorias, ya que el nivel de conanza asignado
a dichas variables no ser el mismo que tendr la respuesta del
sistema modelado en trminos probabilistas. En la tabla 5 se muestran la probabilidades de que el dsi probabilista sea menor que
dsi,det .
Se observa que la probabilidad aumenta con el estado de dano.

leve y moderado se pasa al 2%
Del 0,4% para los estados de dano
severo y al 17,5% para el de colapso.
para el estado de dano
Figura 11. Curvas de fragilidad obtenidas con el enfoque probabilista.
0.3
0.2
0.1
0.1
0.2
0.3
0.4
dsi (m)
Figura 12. Correlacin entre dsi y dsi y anlisis de correlacin.
Tabla 6
Coeciente de correlacin entre dsi y dsi
i,j
ds1 y ds1
ds2 y ds2
ds3 y ds3
ds4 y ds4
0,1967
0,2007
0,8171
0,9720
extensivo y de colapso
el contrario, respecto a los estados de dano
se observa una importante dispersin en ambos parmetros dsi y
dsi . La dispersin en el parmetro dsi , particularmente para las
curvas de fragilidad del estado de colapso, se reeja en la pendiente de las curvas que decae a medida que aumenta el umbral
En la gura 12 se muestra la correlacin entre dsi
del estado de dano.
y dsi . La tabla 6 muestra los coecientes de correlacin, que ponen
de maniesto el incremento del grado de correlacin a medida
que dsi aumenta. Es importante observar que el intervalo en
el que ucta la variable dsi tambin aumenta cuanto mayor
sea dsi , lo que no es ms que otra conrmacin de la
mayor aleatoriedad e incertidumbre asociadas a la respuesta
cuando la estructura empieza a comportarse en rgimen no lineal.
Es interesante buscar una forma compacta de representar las
curvas de fragilidad obtenidas mediante el enfoque probabilista sin
perder la riqueza de la informacin. Para ello ser necesario evaluar
si las variables aleatorias dsi y dsi responden a alguna distribucin de probabilidad, lo que permitira almacenar los resultados
en trminos de las variables que denen dichas distribuciones.
Para hacer esta comprobacin, las variables de salida se agrupan
20
0.1
0.2
0.3
0.5
0
0.1
ds1 (m)
20
0.1
0.2
0.3
P[sd<ds3]
0
0.1
0.2
0.3
0.1
ds2 (cm)
ds2
ds3 (cm)
ds3
ds4 (cm)
ds4
0,267
0,004
12,8
0,2
0,218
0,015
16,1
0,5
0,231
0,03
22,7
1,6
0,28
0,035
0.3
0.9
0.8
0.5
0.7
0
0.1
0.2
0.3
ds3 (m)
P[sd<ds4]
0.2
20
0.1
ds1
9,2
0,2
a
0
ds3 (m)
40
ds1 (cm)
ds2 (m)
20

0.5
ds2 (m)
40
0.3
P[sd<ds2]
40
0.2
ds1 (m)
P[sd > dsi ]
Tabla 7
Valores medios y desviaciones estndar de las variables aleatorias que denen las
curvas de fragilidad
P[sd<ds1]
40
0.2
0.3
71
0.5
0
0.1
ds4 (m)
0.2
0.6
0.5
0.4
0.3
0.3
ds4 (m)
0.2
Variable no paramtrica
Curvas de fragilidad medias

Curvas de fragilidad deterministas
Curvas de fragilidad confianza del 95%
Aproximacin a variable aleatoria normal
0.1
Figura 13. Histogramas para los dsi y comparacin con una variable aleatoria gaussiana.
0
0
0.1
0.2
0.3
0.4
0.5
0.6

0.2
Desviacin estndar CF1
0.18
0.16
0.14
0.12
0.1
0.08
0.06
0.04
0.02
0
0
0.1
0.2
0.3
0.4
0.5
0.6
1
1
P[<ds ]
30
Desviacin estndar (m)
en histogramas y los valores obtenidos se dividen por el nmero

total de muestras. Con los histogramas normalizados se calcula
la funcin de probabilidad cumulativa de la variable, y esta se
compara con la funcin de probabilidad cumulativa dada por un
modelo clsico que, en este caso, es la funcin de probabilidad
cumulativa gaussiana. Las guras 13 y 14 muestran los resultados
obtenidos a partir de estas comparaciones para las variables dsi
y dsi , respectivamente. De esta comparacin se deduce que estas
variables responden adecuadamente a una funcin de probabilidad
gaussiana. Los valores medios y las desviaciones estndar de las
variables aleatorias comparadas en las guras 13 y 14 se muestran
en la tabla 7. Con esto se conocen las caractersticas de las variables aleatorias de salida y se pueden obtener curvas de fragilidad
para cualquier nivel de conanza. En la gura 15a se muestran las
curvas de fragilidad medias, las curvas de fragilidad para un nivel
de conanza del 95% a partir de las variables de salida y las curvas
20
10
0
0
0.1
0.2
ds
0.3
0
0
0.4
0.1
0.2
ds
0.3
0.4
Figura 15. a) Curvas de fragilidad obtenidas con el enfoque probabilista y determinista; adems se muestran las curvas de fragilidad para un nivel de conanza del
95% b) Desviacin estndar de las curvas de fragilidad probabilistas.
5
P[<ds ]
40
0.5
20
0
0
0.1
0.2
ds
0.3
4
0.1
0.2
ds
0.3
0.4
3.5
20
0
0
0.1
0.2
ds
0.3
0
0
0.4
0.1
0.2
ds
0.3
0.4
20
0
0
0.1
0.2
0.3
0.4
2.5
P[<ds ]
40
3
0.5
sa (g)
P[<ds ]
40
0.5
0
0
0.4
Espectro tipo 1, suelo A

Espectro tipo 1, suelo D
4.5
1.5
0.5
1
0
0
0.1
ds
0.2
0.3
0.4
ds
0.5
Variable no paramtrica
Aproximacin a variable aleatoria normal
Figura 14. Histogramas para los dsi y comparacin con una variable aleatoria gaussiana.
0.5
1.5
2.5
3.5
Periodo (s)
Figura 16. Espectros considerados en el anlisis tomado del EC8 para los tipos 1 y
2, suelos tipos A y D.
72
7
Espectro Eurocdigo 8, tipo 1 suelo A
Media espectros reales
Media espectros reales+Desviacin estndar
Periodo fundamental del edificio, enfoque determinista
1
Terremoto representativo espectro 1A
0.8
0.6
0.4
Aceleracin (g)
sa (g)
0.2
0
0.2
0.4
0.6
1
0.8
0
0.2
0.4
0.6
0.8
1.2
1.4
1.6
1.8
10
20
30
Periodo (s)
7
Espectro Eurocdigo 8, tipo 1 suelo D
40
50
60
70
80
Tiempo (s)
1
Terremoto representativo espectro 1D
0.8
0.6
0.4
Aceleracin (g)
sa (g)
0.2
0
0.2
0.4
0.6
1
0.8
0
0.2
0.4
0.6
0.8
1.2
1.4
1.6
1.8
10
20
30
Periodo (s)
Espectro Eurocdigo, tipo 2 suelo A

40
50
60
70
80
Tiempo (s)
1
Terremoto representativo espectro 2A
0.8
0.6
0.4
Aceleracin (g)
sa (g)
0.2
0
0.2
0.4
0.6
1
0.8
0
0.2
0.4
0.6
0.8
1.2
1.4
1.6
1.8
10
20
30
Periodo (s)
h
Espectro Eurocdigo 8, tipo 2 suelo D
50
60
70
80
1
Terremoto representativo espectro 2D
0.8
0.6
0.4
Aceleracin (g)
sa (g)
40
Tiempo (s)
0.2
0
0.2
0.4
0.6
1
0.8
0
0.2
0.4
0.6
0.8
1.2
1.4
1.6
1.8
Periodo (s)
10
20
30
40
50
60
70
80
Tiempo (s)
Figura 17. Espectros de respuesta elstica con 5% de amortiguamiento crtico tomados del EC8; adems se muestra el valor medio de los espectros correspondientes a
acelerogramas reales que son compatibles y la respectiva desviacin estndar. A la derecha de cada espectro se muestra un ejemplo de los acelerogramas seleccionados.
de fragilidad calculadas con el enfoque determinista con el n de

comparar los resultados obtenidos.
En la comparacin de la gura 15a puede verse que las curvas de
fragilidad con un nivel de conanza del 95%, obtenidas a partir del
enfoque probabilista, tienen siempre ordenadas mayores que las
de las curvas de fragilidad que resultan del mtodo determinista.

La diferencia entre estas curvas se acenta a medida que el estado
y el desplazamiento espectral aumenta, lo que demuestra
de dano
la importancia de trabajar este problema desde un punto de
vista probabilista. Por otra parte, para obtener una medida de los
Tal como se ha mencionado anteriormente, es importante caracterizar la accin ssmica teniendo en cuenta su aleatoriedad. En
el apartado 3.7 se ha visto que, para aplicar el mtodo basado
en el espectro de capacidad, es necesario denir la accin ssmica en trminos de espectros elsticos de respuesta con el 5% de
amortiguamiento crtico. La gura 16 muestra los 4 espectros seleccionados (2 por cada tipo del EC8). Puesto que uno de los propsitos
de este trabajo es comparar los resultados obtenidos mediante
variables deterministas y probabilistas, la accin ssmica tambin
debe denirse en trminos probabilistas. Para ello, se utilizan regis
tros de acelerogramas reales de 2 bases de datos, una espanola
y la
otra europea, y se calculan los espectros de respuesta elstica para
cada registro. Los valores obtenidos se comparan con los espectros
suavizados dados por los cdigos y, despus, se seleccionan los que
se aproximen a los propuestos en el EC8. Para cada uno de ellos se
han obtenido 20 registros de aceleracin. En la gura 17 se muestran los espectros elsticos seleccionados del EC8, el valor medio de
los espectros de los acelerogramas reales seleccionados y su valor
medio ms una desviacin estndar; adems, para cada espectro se
muestra uno de los 20 registros seleccionados para representar la
accin ssmica.
A partir de los espectros seleccionados y con el perodo de
vibracin fundamental del edicio, se puede hacer una valoracin
preliminar de los registros que van a producir ms desplazamiento
En la gura 17, que muestran los espectros
y, por tanto, ms dano.
de respuesta, se representa el perodo fundamental del edicio y
se obtiene la ordenada del espectro; cuanto mayor es el valor de
en el edicio. Sin embargo,
dicha ordenada, mayor ser el dano
como el perodo es tambin una variable aleatoria, en la gura 17 se
muestra el perodo medio obtenido con el enfoque probabilista. En
este caso, se dibuja el perodo del edicio sobre el espectro de respuesta normalizado. De aqu se obtiene que los acelerogramas que
la estructura son los compatibles con el espectro
ms van a danar
causarn son los comtipo 1 suelo D y que los que menos dano
patibles con el espectro tipo 2 suelo A. Comparando los registros
con espectros de respuesta compatibles con los espectros tipo 1
suelo A y tipo 2 suelo D del Eurocdigo, se observa que las ordenadas asociadas al perodo fundamental son muy parecidas, por
lo que habra que evaluar los perodos asociados a los modos de
vibracin superiores al fundamental. Sin embargo, hay que tener en
cuenta 2 aspectos en el espectro tipo 2 suelo D: 1) A medida que el
perodo decrece, la aceleracin espectral crece ms rpido que en
el caso del espectro tipo 1 suelo A. 2) La meseta de aceleracin constante del espectro tipo 2 suelo D tiene ordenadas mayores que en
el caso del espectro tipo 1 suelo A. Esto sugiere que los registros
compatibles con el espectro tipo 2 suelo D sern ms perjudiciales
para el edicio que los espectros tipo 1 suelo D. Esta consideracin preliminar ser constatada y conrmada usando los 2 mtodos
simplicados del ATC-40, ALE y PA-8, para estimar el punto de
capacidad por demanda y tambin usando el ADNL.
0.4
0.35
0.3
0.25
0.2
0.15
0.1

0.05
0
0.2
0.4
0.6
0.8
1.2
1.4
1.6
1.8
PGA (g)
0.1
0.09
0.08
4.4. La accin ssmica como variable aleatoria
0.07
0.06
0.05
0.04
0.03
Espectro tipo 1 suelo A
Espectro tipo 1 suelo D
Espectro tipo 2 suelo A
Espectro tipo 2 suelo D
0.02
0.01
0
0.05
0.1
0.15
0.2
0.25
0.3
0.35
0.4

Figura 18. a) Desplazamiento espectral obtenido con la ALE. b) Desviacin estndar
del desplazamiento espectral obtenido con la ALE.
5
Espectro de respuesta en formato sasd
4.5
4
Aceleracin espectral (g)
tramos en que las curvas de fragilidad tienen mayor dispersin,

en la gura 15b se muestra la desviacin estndar en cada curva
de fragilidad como una funcin del desplazamiento espectral, y
adems puede observarse que la mayor incertidumbre ocurre
alrededor del dsi de las curvas de fragilidad medias.
73
3.5
3
2.5
2
1.5
1
4.5. Aproximacin lineal equivalente, enfoque probabilista

Los resultados obtenidos con el enfoque probabilista por medio
de la metodologa ALE se muestran en la gura 18. Como se
mencion en el apartado 3.8, en este trabajo se quiere evaluar la
respuesta de la estructura como una funcin del PGA siguiendo
el procedimiento ADI. En la ALE es suciente escalar los espectros
0.5
0
0.5
1.5
2.5

Figura 19. Muestra de un espectro de respuesta en formato sa-sd para un sismo
compatible con el espectro tipo 1 suelo D dado por el EC8.
74
aplicar el PA-8 tambin se utilizan espectros de acelerogramas reales, al igual que en el apartado anterior cuando se aplic la ALE. Sin
embargo, cuando dichos espectros son transformados al formato
sa-sd, se obtienen funciones que no tienen una forma sencilla, tal
como puede verse en la gura 19. Este hecho complica la automatizacin del proceso iterativo que es de ms fcil aplicacin para
denidos
espectros suavizados, como son los espectros de diseno,
mediante funciones analticas sencillas.
0.4
0.35
0.3
0.25
0.2
ALE 1A
PA8 1A
ALE 1D
PA8 1D
ALE 2A
PA8 2A
ALE 2D
PA8 2D
0.15
0.1
0.05
0
0.2
0.4
0.6
0.8
1.2
1.4
1.6
1.8
1600
PGA (g)
0.1
0.09
ALE 1A
PA8 1A
ALE 1D
PA8 1D
ALE 2A
PA8 2A
ALE 2D
PA8 2D
0.08
0.07
0.06
2000
1800
1400
1200
1000
800
0.06 g
0.12 g
0.18 g
0.24 g
0.30 g
0.36 g
600
400
200
0
0
0.05
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.04
2000
0.03
1900
0.05
0.1
0.15
0.2
0.25
0.3
0.35
0.4

Figura 20. a) Desplazamiento espectral obtenido con el PA-8 y comparado con la
ALE. b) Desviacin estndar del desplazamiento espectral calculado con el PA-8 y
comparado con la ALE.
de respuesta para un solo PGA, ya que el resultado es lineal, por

lo que solo es necesario extender una recta desde el origen que
pase por este punto, con lo que se consigue que el desplazamiento espectral sea una funcin del PGA. Como los espectros de
respuesta y los de capacidad son variables aleatorias, las rectas
que relacionan el PGA y el desplazamiento espectral tambin
sern aleatorias. La gura 18a muestra solo las rectas medias que
permiten comprobar el anlisis hecho a partir de los espectros
de respuesta y el perodo fundamental del edicio. La gura 18b
muestra la desviacin estndar de estas rectas, que tambin vara
de forma lineal. La desviacin estndar es mayor para el espectro
tipo 2 suelo A y menor para el espectro tipo 1 suelo D, por lo que
habra que evaluar la variacin de la relacin entre el PGA y el
desplazamiento espectral cuando se quiere trabajar con curvas
que garanticen niveles de conanza superiores a las asociadas al
valor medio.
1800
1700
1600
1500
0.36 g
1400
0.2
0.25
0.3
0.35
0.4
0.45
0.5
0.55
0.6
0.65
1100
1000
0.01
0.02
900
800
700
600
4.6. Procedimiento PA-8

500
El otro mtodo simplicado propuesto en el ATC-40 que ser

usado en este trabajo es el PA-8, en el que el espectro de demanda
se reduce por la ductilidad de la estructura. Dicho mtodo ha sido
estudiado rigurosamente por Fajfar [27], que concluye que es una
buena aproximacin para obtener el desplazamiento espectral de
una estructura teniendo en cuenta su comportamiento no lineal. Al
0.06 g
400
0.05
0.06
0.07
0.08
0.09
0.1
0.11

Figura 21. Resultados del ADNL probabilista. Zoom a los resultados para valores del
PGA de 0,06 y 0,36 g.
En la gura 20 se comparan los resultados obtenidos con los

mtodos PA-8 y ALE. Esta comparacin permite establecer que con
el ALE se obtienen resultados que estn del lado de la seguridad,
puesto que nunca se subestima el desplazamiento espectral. En el
siguiente apartado se avanza en estas comparaciones y discusin
usando el anlisis dinmico no lineal
Despus de aplicar los mtodos simplicados de evaluacin del
ssmico, se establece que el dano
esperado de la estructura
dano
modelada es muy alto en zonas ssmicas compatibles con el espectro de respuesta tipo 1, suelo D. Por ejemplo, para un PGA de 0,3 g, el
desplazamiento espectral asociado es aproximadamente de 0,25 m,
valor que corresponde a una alta probabilidad de colapso. Por otra
parte, si la zona ssmica es compatible con el espectro de res esperado es menor. Esto se reeja
puesta tipo 2, suelo A, el dano
en el hecho de que para que el colapso sea probable debera ocurrir
un terremoto con un PGA superior a 1,5 g, valor que es altamente
improbable, y ms teniendo en cuenta que este tipo de espectros
corresponden a sismos de baja magnitud.
5. Anlisis dinmico no lineal. Enfoque probabilista
En este apartado se muestra el anlisis dinmico no lineal desde
una perspectiva probabilista. En la gura 20 se aprecia que las acciones ssmicas que producen mayores desplazamientos espectrales
en la estructura para un PGA menor son las compatibles con el
espectro tipo 1, suelo tipo D. Para los otros espectros, tendra que
escalarse la aceleracin hasta valores que son poco probables en

son
donde las PGA para un periodo de retorno de 500 anos
Espana,
inferiores a 0,25 g (NCSE-02, 2002). Por tanto, en este artculo solo
se hace el anlisis dinmico incremental (ADI) para los sismos compatibles con el espectro mencionado, suciente para el objetivo de
este trabajo. De acuerdo con los resultados obtenidos mediante los
procedimientos propuestos en el ATC-40, para los sismos compatibles con el espectro tipo 1 y suelo tipo D, se escala el PGA hasta
0,36 g para llegar al desplazamiento de colapso. El PGA se incrementar en intervalos de 0,06 g hasta llegar a 0,36 g. Para mantener
el enfoque probabilista, las propiedades mecnicas de los materiales se han tratado como variables aleatorias, tal y como se hizo en
el anlisis esttico pero generando solo 100 muestras aleatorias de
las propiedades mecnicas para el conjunto de estas variables ya
que, como se ver en los resultados, la mayor fuente de aleatoriedad proviene de la accin ssmica. Dicha accin se representa por
20 acelerogramas seleccionados entre datos reales. Al nal se ejecutarn 2.000 anlisis dinmicos no lineales por cada PGA. En la
gura 21 se muestran los resultados obtenidos despus de ejecutar
todos los ADNL. Se puede apreciar cmo evoluciona el desplazamiento mximo en el techo y el cortante en la base a medida que
aumenta el PGA. Si se hace una ampliacin de los resultados obtenidos con el primer PGA, que es de 0,06 g, se puede ver claramente que
la principal fuente de aleatoriedad es la accin ssmica, ya que se
7
Espectro Eurocdigo tipo 1 suelo D
Espectro medio
Espectro caracteristico
Espectros reales
0.4
6
ALE
PA8
ADNL
Dy, Dy fiabilidad 95%
Du, Du fiabilidad 95%
0.35
0.3
0.25
sa (g)
75
0.2
0.15
2
0.1
1
0.05
0
0
0
0.05
0.1
0.15
0.2
0.25
0.3
0.35
0.4
0.5
1.5
0.07
ALE
PA8
ADNL
0.06
0.05
0.04
0.03
0.02
0.01
2.5
3.5
0.25
0.3
0.35
0.4
0.4
Ale, enfoque determinista
PA8, enfoque determinista
ADNL, enfoque determinista
Dy, Dy fiabilidad 95%
Du, Du fiabilidad 95%
0.35
Periodo (s)
PGA (g)
0.3
0.25
0.2
0.15
0.1
0.05
0.05
0.1
0.15
0.2
0.25
0.3
0.35
0.4
PGA (g)
Figura 22. a) Relacin entre el PGA y el desplazamiento espectral obtenido con todos
los mtodos simplicados y el ADNL; enfoque determinista. b) Comparacin de la
desviacin estndar.
0.05
0.1
0.15
0.2
PGA (g)
Figura 23. a) Espectro caracterstico obtenido a partir de los acelerogramas reales. b) Resultados del enfoque determinista utilizando el espectro caracterstico y el
acelerograma articial generado a partir de aquel.
76
pueden distinguir con claridad 20 estructuras de nubes asociadas a

cada acelerograma, cada nube contiene 100 puntos correspondientes a los resultados asociados a las propiedades mecnicas como
variable aleatoria. Vale la pena observar que, a medida que aumenta
el PGA, estas nubes pierden su organizacin, su dispersin
aumenta de forma signicativa y, en consecuencia, tambin se
incrementa la desviacin estndar; es decir, cuando el PGA
aumenta, la inuencia de las propiedades mecnicas de los materiales en la respuesta tambin lo hace. Este es un resultado importante
de este estudio.
Para comparar estos resultados con los obtenidos previamente por las metodologas del ATC-40, se calcula el valor
medio y la desviacin estndar para cada PGA. En los cambios de
pendiente de la desviacin estndar se aprecian variaciones signicativas del comportamiento de la estructura, como exponen Vargas
et al. [28]. La gura 22a muestra con todas las metodologas usadas
el promedio del desplazamiento esperado y la gura 22b muestra
la desviacin estndar. En la gura 22a se muestran tambin los
desplazamientos de plasticacin Dy y de capacidad ltima Du,
junto con sus respectivos intervalos de conanza del 95%.
Esto quiere decir que las aproximaciones propuestas en el
ATC-40 no subestiman el desplazamiento espectral esperado si se
comparan con los resultados obtenidos con el anlisis dinmico
no lineal. Sin embargo, el PA-8 da una mejor estima del desplazamiento espectral que la ALE si se comparan los resultados con el
ADNL.
Para calcular con el ADNL en trminos deterministas la relacin

entre el PGA y el desplazamiento espectral esperado, es necesario
disponer de una accin ssmica que sea representativa para los acelerogramas seleccionados en el anlisis aleatorio. Por tanto, a partir
del espectro medio y la desviacin estndar de los espectros reales
seleccionados y usados en el anlisis dinmico no lineal se procede
a calcular un espectro que ser llamado espectro caracterstico, cuya
probabilidad de ser excedido es del 5%, tal y como se muestra en
la gura 23a. Con el espectro caracterstico as denido se genera
un acelerograma sinttico a partir del procedimiento de Gasparini
y Vanmarcke [29], usando la envolvente trapezoidal propuesta
por Hou [30]. Como en el modelo de clculo no se considera la
degradacin de la rigidez de los elementos, la duracin de la senal

generada no inuye en la respuesta de la estructura [31], por lo
que la duracin del acelerograma sinttico se ha denido como
la duracin media de los 20 acelerogramas usados en el ADNL
probabilista. Para calcular el desplazamiento espectral esperado
se usan, adems del ADNL, las 2 aproximaciones del ATC-40.
Utilizando el enfoque determinista con el espectro caracterstico y
el correspondiente acelerograma sinttico, se calculan resultados
para los mismos PGA considerados en el estudio probabilista. La
gura 23b muestra los resultados obtenidos despus de aplicar
los 3 mtodos mencionados y el enfoque determinista, usando la
0.9
0.9
0.8
0.8
0.7
0.7
0.6
0.6
ID
ID
6. Comparacin de los resultados obtenidos con el enfoque

probabilista y determinista
0.5
0.5
0.4
0.4
0.3
0.3
0.2
0.2
ID Curvas de fragilidad medias
ID Curvas de fragilidad deterministas
ID Curvas de fragilidad 95%
0.1
ID curvas de fragilidad medias

ID curvas de fragilidad deterministas
ID curvas de fragilidad 95%
0.1
0
0
0
0.05
0.1
0.15
0.2
0.25
0.3
0.35
0.05
0.1
0.2
0.9
0.8
0.8
0.7
0.7
0.6
0.6
0.5
0.3
0.35
0.5
0.4
0.4
0.3
0.3
0.2
0.2
ID curvas de fragilidad medias
ID Curvas de fragilidad deterministas
ID Curvas de fragilidad 95%
0.1
0.25
0.9
ID
ID
0.15
PGA (g)
PGA (g)
ID ALE, fiabilidad del 95%

ID PA8, fiabilidad del 95%
ID ADNL, fiabilidad del 95%
0.1
0
0
0.05
0.1
0.15
0.2
PGA (g)
0.25
0.3
0.35
0.05
0.1
0.15
0.2
0.25
0.3
0.35
PGA (g)
esperado a partir de la ALE. b) ndice de dano

esperado a partir del PA-8. c) ndice de dano
esperado a partir del ADNL. d) ndice de dano
esperado
Figura 24. a) ndice de dano
por todas las metodologas para un nivel de conanza del 95%.
accin compatible con el espectro caracterstico, tal como se ha

denido ms arriba. De nuevo, se observa que las metodologas
simplicadas ALE y PA-8 son conservadoras respecto al anlisis
dinmico no lineal; adems, el PA-8 sigue siendo una mejor estima
para el clculo del desplazamiento espectral que la ALE, si se
comparan los resultados con el ADNL.
En la gura 24 se pueden ver los ID calculados a partir de los
procedimientos del ATC-40 y el ADNL. En cada gura se compara el
ID obtenido despus de utilizar las curvas de fragilidad medias y las
curvas de fragilidad para un nivel de conanza del 95% (g. 15a).
Tambin se muestran los resultados obtenidos mediante el enfoque
determinista. En estas guras se puede apreciar la gran dispersin
que existe en los resultados, que se ve reejada en las diferencias signicativas entre el ID medio correspondiente a un nivel
de conanza del 50% en las curvas de fragilidad y el ID calculado
para un nivel de conanza del 95%. Los resultados conrman que
suponer niveles de conanza del 95% en las variables de entrada
y enfocar el problema desde un punto de vista determinista no
garantiza que las variables de salida tengan este mismo nivel de
especonanza. En la gura 24 puede verse que el ndice de dano
rado que resulta del enfoque determinista es menor que el obtenido
en trminos probabilistas para un nivel de conanza del 95% en
la variable de salida. La diferencia entre estos enfoques es ms
acentuada cuando se usan procedimientos que tienen en cuenta
el comportamiento no lineal de la estructura, es decir, el PA-8 y
el ADNL. Por otra parte, las aproximaciones propuestas en el ATC40 no subestiman el desplazamiento espectral si se comparan los
resultados con el ADNL. En la gura 24d se comparan los ID obtenidos con los 3 mtodos para un nivel de conanza del 95%. Se
han logrado resultados similares para edicios de la misma tipologa estructural pero con diferentes alturas, como puede verse
en [32].
77
Es importante resaltar que si el sistema tiene un comportamiento no lineal, no se garantiza que el nivel de conanza de la
respuesta obtenida con el enfoque probabilista sea el mismo que el
asignado a las variables de entrada cuando se hace un clculo determinista. Este hecho es visible en las diferencias signicativas que
esperado, ID, obtenidos a partir del
existen entre los ndices de dano
desplazamiento espectral calculado mediante el anlisis dinmico
esperado obtenido con
no lineal. Concretamente, el ndice de dano
un enfoque determinista llega a ser en algunos casos un 20% menor
que el que se obtiene con un enfoque probabilista.
Para el clculo del desplazamiento esperado se ha hecho un estudio comparativo sobre los enfoques determinista y probabilista, y
los mtodos simplicados del ATC-40, como la aproximacin lineal
equivalente y el que tiene en cuenta la reduccin por la ductilidad
de la estructura. De esta comparacin se concluye que, si se utilizan como referencia los resultados del anlisis dinmico no lineal,
el procedimiento simplicado con base en la reduccin del espectro de demanda por la ductilidad de la estructura proporciona una
mejor estima que la aproximacin lineal equivalente. Sin embargo,
cabe anotar que las 2 metodologas simplicadas estn del lado de
la seguridad.
Otro resultado relevante y muy interesante del estudio probabilista realizado en este artculo es la matriz de correlacin entre
las variables de entrada y de salida, que proporciona una valiosa
de
informacin que puede ser muy til tanto en la fase de diseno
las estructuras nuevas como en el caso de la evaluacin del riesgo
ssmico de las estructuras existentes.
Finalmente, es importante recalcar que en la evaluacin de la
vulnerabilidad y del riesgo ssmico, sea cual sea la metodologa
utilizada, es importante seguir un enfoque que tenga en cuenta el
comportamiento no lineal de la estructura, la aleatoriedad de las
propiedades mecnicas de los materiales y la enorme incertidumbre asociada a la accin ssmica.
7. Conclusiones
En este trabajo se ha propuesto una metodologa de evaluacin
probabilista de la vulnerabilidad y del riesgo ssmico de estructuras de hormign armado utilizando el mtodo de Monte Carlo.
Se ha tenido en cuenta el carcter aleatorio no solo de las propiedades mecnicas de los materiales sino tambin de la accin
ssmica. Esta ltima ha sido considerada mediante acelerogramas reales compatibles con los espectros de respuesta elstica
dados en el EC8. Para ello, se han considerado los 2 tipos de
espectros del EC8 para 2 tipos de suelo, para abarcar dife ssmica. Asimismo,
rentes contenidos frecuenciales de la senal
se han evaluado la vulnerabilidad y el riesgo mediante un
enfoque determinista a partir de los valores caractersticos de
las propiedades mecnicas de los materiales y deniendo una
accin ssmica caracterstica a partir de acelerogramas reales.
ssmico esperado se evala mediante un ndice calEl dano
culado a partir de las curvas de fragilidad que, a su vez, son
funciones que dependen del desplazamiento espectral. Dicho
desplazamiento ha sido calculado por medio de mtodos simplicados propuestos en el ATC-40. Los resultados obtenidos con
estos mtodos han sido comparados con el anlisis dinmico
no lineal.
Una importante conclusin es que, a pesar de trabajar con herramientas avanzadas de anlisis estructural, como son el anlisis
esttico incremental y el anlisis dinmico no lineal, estos procedimientos muestran incertidumbres signicativas cuando se tiene
en cuenta la aleatoriedad de las variables asociadas al problema. Es
necesario subrayar que en este trabajo se han considerado coe para las variables de
cientes de variacin relativamente pequenos
entrada en relacin con las incertidumbres que pueden existir en
estructuras que no han sido proyectadas de acuerdo con criterios
sismorresistentes.
Agradecimientos
Este trabajo ha sido parcialmente nanciado por el Instituto
(IGC), por el Ministerio de Ciencia e
Geogrco de Cataluna
por la Comisin Europea y por una beca
Innovacin de Espana,
para los siguientes proyectos de investigacin: CGL-2005-04541C03-02/BTE, CGL2008-00869/BTE, CoPASRE, CGL2011-29063,
INTERREG: POCTEFA 2007-2013/73/08 y MOVEFT7-ENV2007-1-211590. El autor Jorge Eduardo Hurtado agradece a la
Universidad Nacional de Colombia, Sede Manizales, la nanciacin
para la
de una pasanta en la Universidad Politcnica de Cataluna
realizacin de esta investigacin.
Bibliografa
[1] C. Bermdez, Vulnerabilidad ssmica de edicios de acero [tesis doctoral],
Departamento de Ingeniera del Terreno, Cartogrca y Geofsica, Universidad
2010.
Politcnica de Cataluna,
[2] J.E. Hurtado, Modelizacin estocstica de la accin ssmica, Centro Internacional
1999.
de Mtodos Numricos en Ingeniera (CIMNE), IS-33, Barcelona, Espana,
[3] A.H. Barbat, F. Ypez Moya, J.A. Canas, Damage scenarios simulation for risk
assessment in urban zones, Earthq. spectra 2 (3) (1996) 371394.
[4] A.H. Barbat, U. Mena, F. Ypez, Evaluacin probabilista del riesgo ssmico en
ing. 14 (2) (1998) 247268.
zonas urbanas, Rev. int. mtodos numr. clc. diseno
[5] B. Borzi, R. Phino, H. Crowley, Simplied pushover analysis for large-scale
assessment of RC buildings, Eng. Struct. 30 (2008) 804820.
[6] A.H. Barbat, L.G. Pujades, N. Lantada, R. Moreno, Seismic damage evaluation in
urban areas using the capacity spectrum method: application to Barcelona, Soil
Dyn. Earthquake Eng. 28 (2008) 851865.
[7] N. Lantada, L.G. Pujades, A.H. Barbat, Vulnerability index and capacity spectrum
based methods for urban seismic risk evaluation. A comparison, Nat. Hazard.
Rev. 51 (2009) 501524.
[8] L.G. Pujades, A.H. Barbat, R. Gonzlez-Drigo, J. Avila, S. Lagomarsino, Seismic
performance of a block of buildings representative of the typical construction
in the Eixample district in Barcelona (Spain), Bull. Earthquake Eng. 10 (1) (2012)
331349.
78
[9] A.M. Mwafy, A.S. Elnashai, Static pushover versus dynamic collapse analysis of
RC buildings, Eng. Struct. 23 (2001) 407424.
[10] M. Poursha, F. Khoshnoudian, A.S. Moghadam, A consecutive modal pushover
procedure for estimating the seismic demands of tall buildings, Eng. Struct. 31
(2009) 591599.
[11] S.P. Kim, Y.C. Kuruma, An alternative pushover analysis procedure to estimate
seismic displacement demands, Eng. Struct. 30 (2008) 37933807.
[12] J.C. Vielma, A.H. Barbat, S. Oller, Seismic performance of wafed-slabs oor
buildings, Proc. Inst. Civ. Eng. Struct. Build. 162 (SB3) (2009) 169182.
[13] J.C. Vielma, A.H. Barbat, S. Oller, Seismic safety of limited ductility buildings,
Bull. Earthquake Eng. 8 (1) (2010) 135155.
[14] M. Fragiadakis, D. Vamvatsikos, Estimation of uncertain parameters using static
pushover methods, Safety Reliability and Risk of Structures, Infrastructures and
Engineering Systems (2010).
[15] M. Dolsek, Effects of uncertainties on seismic response parameters of reinforced
concrete frames, Safety Reliability and Risk of Structures, Infrastructures and
Engineering Systems (2010).
[16] Eurocdigo 8, Design of structures for earthquake resistance. Part 1: General
rules, seismic actions and rules for building, 2004.
[17] I. Satyarno I, Pushover analysis for the seismic assessment of reinforced concrete buildings [doctoral thesis], Department of Civil Engineering, University
of Canterbury, 1999.
[18] A.J. Carr, Ruaumoko-Inelastic Dynamic Analisys Program, Department of
Civil Engineering, University of Canterbury, Christchurch, New Zealand,
2000.
[19] P.A. Mata, S. Oller, A.H. Barbat, Static analysis of beam structures under nonlinear geometric and constitutive behaviour, Comput. Meth. Appl. Mech. Eng.
196 (2007) 44584478.
[20] J. Faleiro, S. Oller, A.H. Barbat, Plastic-damage seismic model for reinforced
concrete frames, Comput. Struct. 86 (7-8) (2008) 581597.
[21] ATC-40, Seismic evaluation and retrot of concrete buildings, Applied Technology Council, Redwood City, California, 1996.
[22] HAZUS 99, Earthquake Loss Estimation Methodology Earthquake. Technical
Manual, Vol. 1, Federal Emergency Management Agency (FEMA). Washington
D.C., 1999.
[23] RISK UE Project of the European Commission (2004) An advanced approach
to earthquake risk scenarios with applications to different European towns.
Contract number: EVK4-CT-2000-00014.
[24] J.A. Mahaney, T.F. Paret, B.E. Kehoe, S.A. Freeman, The capacity spectrum method for evaluating structural response during the Loma Prieta earthquake,
National Earthquakes Conference, Memphis (1993).
[25] D. Vamvatsikos, C.A. Cornell, The incremental dynamic analysis, Earthquake
Eng. Struct. Dyn. 31 (3) (2001) 491514.
[26] M. Kalos, P.A. Whitlock, Monte Carlo Methods, Vol. I, Basics, 1986.
[27] P. Fajfar, Capacity spectrum method based on inelastic demand spectra, Earthquake Eng. Struct. Dyn. 28 (1999) 979993.
[28] Y.F. Vargas, L.G. Pujades, A.H. Barbat, J.E. Hurtado, Probabilistic Assessment of
the Global Damage in Reinforced Concrete Structures, en: Proceeding of the
14th European Conference on Earthquake Engineering, Ohrid, 2010.
[29] D. Gasparini, E.H. Vanmarcke, Simulated earthquake motions compatible with
prescribed response spectra. M.I.T. Department of Civil Engineering. Research
report R76-4, Order No. 527, 1976.
[30] S. Hou, Earthquake simulation models and their applications, M.I.T. Department
of Civil Engineering. Research Report. (1968), R68-17.
[31] J. Hancock, J. Bommer, A state of knowledge Review of the inuence of strong
motion duration on structural damage, Earthq. spectra 22 (3) (2006) 827845.
[32] Y.F. Vargas, L.G. Pujades, A.H. Barbat, J.E. Hurtado, Incremental dynamic analysis
and pushover analysis. A probabilistic comparison, en: III ECCOMAS Thematic
Conference on Computational Methods in Structural Dynamics and Earthquake
Engineering, Corfu, Greece, 2011.

Vargas, Pujades, Barbat, Hurtado, 2013

Cargado por

Copyright:

Formatos disponibles

Vargas, Pujades, Barbat, Hurtado, 2013

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Vargas, Pujades, Barbat, Hurtado, 2013

Cargado por

Copyright:

Formatos disponibles

Documento descargado de http://http://www.elsevierciencia.es el 20/06/2013.

Rev. int. mtodos numr. clc. diseo ing. 2013;29(2):6378

Revista Internacional de Mtodos Numricos para

Universidad Politcnica de Catalu

informacin del artculo

Historia del artculo:

Probabilistic assessment of the capacity, fragility and seismic damage

Autor para correspondencia.

ssmicas que no se disenaron

macrossmicas y el comportamiento estructural mediante un ndice

de una tipologa estructural frecuentemente utilizada en Espana

2. Descripcin del edicio

3.1. Anlisis esttico no lineal y curva de capacidad

Valor medio (kPa)

Desviacin estndar (kPa)

fallo por momento ector y por cortante. Por ahora se ha usado

3.2. Espectro de capacidad

Cortante en la base (kN)

Figura 1. Edicio estudiado, con H = 25,65 m, B = 24 m y perodo fundamental de

Aceleracion Espectral (g)

espectrales correspondientes a los umbrales de los 4 estados de

ds3 = Dy + 0, 25(Du Dy)

Desplazamiento espectral (m)

Figura 3. Espectro de capacidad y la representacin bilineal.

desplazamiento espectral con la aceleracin espectral mediante las

El subndice i de las ecuaciones (1) y (2) hace referencia a los

un formato bilineal, lo que es til para denir los estados de dano.

3.4. Curvas de fragilidad

donde sd es el desplazamiento espectral y dsi es la desviacin

Desplazamiento espectral (m)

Desplazamiento espectral (m)

para los dsi correspondientes al edicio calculado. En esta misma

esperado calculado con el enfoque determinista.

donde n es el nmero de estados de danos

reales tomados de bases de datos espanolas

Actualmente, la mayora de los pases proponen, en sus normas

En este trabajo se usa el programa Ruaumoko [18], que permite

Punto de capacidad por demanda

En el apartado anterior se ha hecho una descripcin de mto fsico esperado de

3.7. Punto de capacidad por demanda

4.1. Anlisis esttico no lineal incremental mediante el mtodo

Cortante en la base (kN)

Curva de capacidad media

Desplazamiento en el techo (m)

ltimo de la estructura tiende a ser levemente mayor a medida que

el grado de dependencia entre las variables de entrada y salida

donde ij dene la correlacin entre las variables aleatorias xi y xj ,

y la rigidez del tramo elstico de la curva de capacidad, denominada

4.2. El espectro de capacidad bilineal como variable aleatoria

Valores medios y desviaciones estndar de los umbrales de los estados de dano

Aceleracin espectral (g)

Curvas de fragilidad del ds2

Curvas de fragilidad del ds3

Curvas de fragilidad del ds4

Desplazamiento espectral (m)

Desplazamiento espectral (m)

obtenidos con el enfoque probabilista y el determinista.

4.3. Curvas de fragilidad probabilistas

P[ds2 < ds2,det ]