Review of Seismic Fragility Analysis of Building Structure
-
摘要: 地震易损性是指在不同强度水平的地震作用下建筑结构发生各种破坏的条件概率。本文对地震易损性分析的发展历程进行回顾和总结,介绍了地震易损性的相关概念、研究目的和意义以及5个表达工具。随后,针对国内外学者对建筑结构的易损性分析方法开展了系统的梳理和分析,特别是从基于震害调查的经验分析法和基于数值模拟的理论计算法方面展开了重点综述,并简单介绍了一些不常用的方法。最后指出了国内外对建筑结构易损性分析研究的不足和未来的研究方向。Abstract: Seismic fragility is the probability that structural demand equals or exceeds a limit state conditional on a seismic ground motion given a specified intensity.In this paper, the development of seismic fragility analysis is summarized.The related concepts, purpose and significance and five expressions of seismic fragility are introduced.Then we review the empirical analysis method based on the seismic damage investigation, as well as the theoretical calculation method based on numerical simulation.Some methods that are not commonly used are analyzed too.Finally, the deficiencies and future research directions of fragility analysis of building structures are discussed.
-
Key words:
- Seismic fragility /
- Building structure /
- Numerical simulation /
- Damage investigation
-
引言
随着社会经济的快速发展,城市化程度的不断提高,城市人口的急剧增加,城市地上空间的开发与利用已接近极限,由此促进了城市地下空间的开发与利用,促进了以地铁盾构隧道为代表的地下结构工程的大规模兴建,下层盾构隧道(简称下层隧道)近距离穿越上层盾构隧道(简称上层隧道)的工程也越来越多。如北京地铁10号线三元桥站至亮马桥站盾构隧道左右线并行,最小净距仅为1.7m,6号线南锣鼓巷站至东四站区间上下行隧道叠落,间距仅为2.1m。其抗震问题成为城市工程抗震和防灾减灾研究的重要组成部分(权登州等,2015)。
针对地下结构抗震问题的设计和计算方法,从力学特性上可分为动力时程法和拟静力计算法两大类。动力时程法能够模拟地震波输入过程中各个时刻结构的内力和变形状态,准确地反映地下结构的地震响应特性。当然,在进行动力分析时需要注意动力加载方式和人工边界、力学阻尼及波的传播,以便完整地模拟土体与结构在地震过程的惯性作用以及相互作用。拟静力法主要有地震系数法(施仲衡,1997)、自由场变形法(Hashash等,2001)、反应位移法(川岛一彦,1994)、Pushover分析方法(刘晶波等,2008)以及整体式反应位移法(刘晶波等,2010)等。其中反应位移法的基本原理是通过地层位移来计算地下结构的地震响应,在一定程度上反映地下结构地震反应的特点。该方法理论清晰、计算简便,在众多的设计规范和工程中得到了广泛应用。
目前反应位移法只适用于单个结构,但在强震作用下,小间距盾构隧道夹土层的地震动反应对结构的内力和变形影响很大,现采用的计算模型存在一定缺陷,会增大计算结果的误差率。本文通过改进反应位移法的计算模型,对与下层隧道呈不同角度、间距、埋深等多种工况下的上层隧道横断面进行抗震研究,并与准确的动力有限元的模拟数值进行对比分析,同时对两种算法生成的数据进行曲线拟合,得到的拟合理论模型对提高反应位移法计算结果的精度具有重要的实际意义,也为后续深入研究奠定了基础。
1. 反应位移法概述
在地震动作用下,盾构隧道变形主要受周围地层的约束。反应位移法的基本原理是将地震荷载作用下地层周围的剪力、结构自身惯性力等施加于结构,同时把地层在地震时产生的位移差通过地基弹簧以静载的形式作用于结构上,从而求得结构的内力和变形等。
反应位移法最重要的一步是在结构周围设置地基弹簧,以模拟隧道与地层之间的相互作用,其系数可按计算公式或通过有限元方法求解。但隧道间距较小时,夹土层在强震作用下对结构产生不可忽略的影响,若采用上述方法只分析单个隧道结构,未考虑夹土层区域弹簧参数变化的影响,将使结构内力的计算结果产生误差(唐金良,2014)。
2. 改进反应位移法
2.1 改进模型
为克服反应位移法的计算模型未考虑夹土层弹簧变化这一不足,在反应位移法内容合理的基础上,改进反应位移法的计算模型。
基于有限元结构分析软件SAP84,改进计算模型:两条隧道夹土层等效为连接结构上两点的直弹簧,如图 1所示。结构的受力分析与其周围其他地基弹簧按反应位移法设置(蒋英礼等,2013)。
2.2 夹土层弹簧刚度
结合地基弹簧常数的近似取值法(林皋,1990),拟定夹土层弹簧的计算公式如下:
$$k=\frac{3G\ \cdot A}{S}$$ (1) 式中,k为夹土层弹簧刚度;G为与地震震动最大应变幅度相应的地基土的剪切模量;A为夹土层弹簧反映的土体面积;S为两条隧道结构的间距。
2.3 拟合理论模型
采用上述改进模型对两个结构进行计算,不可避免会产生一些偏差,为简单有效地解决这个问题,通过与准确的动力时程法进行对比分析,提出了引入拟合理论模型来提高计算精度,使改进方法与动力时程法的结果相吻合。
在Matlab软件中对两种算法生成的数据进行拟合时,对函数对话框“Results”中的参数SSE(拟合误差方差的平方和,其值越小,拟合结果越好)与R-square(置信区间,越趋近1,拟合结果越好)进行不断对比和分析,同时参考Adjusted R-square(调整置信区间,越趋近1,拟合度越高),选出较合理的理论模型(唐家德等,2008;史立新等,2007)。
3. 算例
地铁隧道的相对位置主要存在并行、斜穿与下穿这3类情况,本章利用改进模型,在多种埋深、间距工况下,研究隧道的横向地震响应特性,并与动力时程法的数值结果进行对比分析。
基于Midas-NX有限元分析软件对算例进行分析。为提高计算精度(兰景岩等,2012),动力时程法所使用的地震波依据《北京地铁16号线工程场地地震安全性评价报告》,采用提供的50年超越概率为10%的加速度时程曲线,见图 2。土体遵从经典的Mohr-Coulomb破坏准则,盾构隧道结构采用弹性模型。进行静力计算时,土层两侧采用水平约束的人工天然边界;进行动力计算时,在模型底部竖直方向施加地震动的加速度边界,四周为粘弹性动力人工边界;动力分析是在静力计算获得应力场后进行的。
3.1 计算工况
以北京地铁17号线03标段区间盾构隧道为背景,衬砌外径为6m,内径为5.4m,厚度为0.3m。通过查阅并统计北京市地铁盾构隧道资料,设定上下层隧道所呈角度θ=0°、30°、60°、90°,间距S=2m、4m、6m、8m、10m、12m、14m,上层隧道埋深H=8m、10m、12m、14m、16m。隧道尺寸与计算简图详见图 3。盾构隧道周围土层简化为代表北京地区常见的均一粉质黏土(施春花等,2009),结构与土体的参数如表 1所示。
表 1 土体、隧道计算参数Table 1. Physical parameters of soils and tunnels名称 模量/Pa 弹性模量 泊松比 重度/kN·m-3 粘聚力/Pa 内摩擦角/° 压缩模量 动剪切模量 土层 9.1×107 8.92×107 0.33 19.4 31 18 衬砌 3×1010 0.20 25 3.2 数值结果及分析
盾构隧道设计主要依据结构的弯矩值、轴力值及直径变化率,其值越大,结构变形越明显。在进行变化规律分析时,选取最大值进行数据处理。
3.2.1 弯矩分析
(1)对比分析
呈不同角度的两个隧道的变形模态基本一致,以上层隧道为例进行论述。上层隧道弯矩的最大值随埋深和间距呈规律性变化,见图 4,其中以动力时程法的模拟结果作参照,验证改进反应位移法的准确性。由图 4易知两种算法下弯矩最大值的变化规律基本一致:
① 通过图 4(a)与(b)可知,埋深对弯矩值的影响最大,间距次之;与下层隧道位于同一角度、间距工况下的上层隧道,其弯矩最大值与单孔隧道的最大值均随埋深的增加而增加(杨宁波,2013)。
② 通过图 4(c)与(d)可知,在0°并行的情况下,两隧道的上层隧道在同一埋深。随着间距增大,弯矩最大值由先小于单孔隧道,然后增大到与单孔隧道近乎相同。在地震动作用下,当两条隧道间距较小时,由于下层隧道-夹土层-上层隧道惯性作用增强,提高整体结构稳定性,从而减小了上层隧道的差异位移,对上层隧道起到抗震的效果。随间距逐渐增大,抗震效果也随之减弱。随埋深的增加,上层隧道弯矩最大值趋于单孔隧道的速率越来越快。
③ 通过图 4(f)—(j)可知,在30°、60°并行与90°叠落的情况下,当上层隧道在同一埋深下,随着间距增大,弯矩最大值由先大于单孔隧道,然后减小到与单孔隧道近乎相同。下层隧道的开挖造成临近上层隧道的一定范围内土层应力的释放,在强震中,松弛的地层明显放大了上层隧道的地震反应特性,使其弯矩值出现大于单孔隧道的可能。弯矩最大值在隧道间距处于4—6m时出现峰值,该范围是设计双线隧道的薄弱区域。
在横断面抗震分析中应考虑内力值的分布区域(以30°并行的上层隧道为例,如图 5所示),弯矩最大值主要分布在上部拱腰两侧。
采用改进反应位移法计算时,在同一埋深与角度下,当隧道间距S≥8m时,弯矩值的变化不大,并且趋近于单孔隧道,这说明当两条隧道间距S≥8m时,可忽略夹土层地基弹簧的变化对上层隧道弯矩值的影响,按现有的反应位移法进行计算。
(2)拟合理论模型选取
基于Matlab软件进行计算,经残差曲线的对比分析,弯矩值拟合理论模型采用函数Cubic Polynomial(公式2)进行拟合,效果良好,如图 6所示。
$$f(x)={{p}_{1}}{{x}^{3}}+{{p}_{2}}{{x}^{2}}+{{p}_{3}}x+{{p}_{4}}$$ (2) 式中,f(x)为动力时程法弯矩值;x为改进反应位移法弯矩值;p1、p2、p3、p4为修正系数,其均值与取值范围分别为-1.00×10-5(-4.88×10-5,2.88×10-5)、1.16×10-3(-1.42×10-2,1.65×10-2)、1.20(-0.76,3.15)、-8.58(-88.69,71.53)。
拟合优度方面,拟合误差平方和(SSE)为0.04;置信区间(R-square)为0.95;调整置信区间(Adjusted R-square)为0.95。
3.2.2 轴力分析
(1)对比分析
轴力与弯矩的分析方法一致,即对每种工况下的最大值进行研究。两种算法下轴力最大值的变化规律基本一致,同时与改进位移法弯矩最大值的变化规律也一致,以30°并行的上层隧道为例,见图 7。
在地震动作用下,轴力最小值分布在隧道顶部区域,最大值基本上分布在结构下部拱腰附近,集中现象比较明显,以30°并行的上层隧道为例,见图 8。
(2)拟合理论模型选取
基于Matlab软件进行计算分析,轴力值拟合理论模型采用Sum of Sin Functions函数(公式3)进行拟合,效果良好,曲线见图 9。
$$f(x)={{a}_{1}}\ \cdot \ \sin ({{b}_{1}}x+{{c}_{1}})+{{a}_{2}}\ \cdot \ \sin ({{b}_{2}}x+{{c}_{2}})$$ (3) 式中,f(x)为动力时程法轴力值;x为改进反应位移法轴力值;a1、b1、c1、a2、b2、c2为修正系数,其均值与取值范围分别为1.51×103(5.69×102,2.45×103)、6.75×10-4(-5.366×10-5,1.40×10-3)、6.25(5.98,6.51)、28.46(14.42,42.5)、1.41×10-2(9.02×10-3,1.92×10-2)、-4.90(-10.42,0.63)。
在拟合优度方面,拟合误差平方和(SSE)为0.07;置信区间(R-square)为0.96;调整置信区间(Adjusted R-square)为0.96。
3.2.3 直径变化率分析
(1)对比分析
隧道属于柔性结构,选用隧道直径变化率作为隧道横向抗震的性能指标。直径变化率为${{k}_{0}}=\frac{\Delta {{D}_{\max }}}{D}$,其中$\Delta {{D}_{\text{max}}}$为直径的最大变化值(姚毅超,2012)。根据改进反应位移法生成的数据,分析隧道横截面变形量,得到直径的最大变化值,从而求出直径最大变化率,计算结果如表 2所示。
表 2 直径变化率(‰)Table 2. Diameter variation rate(‰)工况 埋深/m 8 10 12 14 16 单孔隧道 间距/m 1.93 1.72 1.54 1.31 1.19 2 1.28 1.19 1.06 0.97 0.91 4 1.30 1.20 1.07 1.00 0.93 6 1.28 1.21 1.04 0.98 0.88 0°并行 8 1.26 1.20 1.03 0.99 0.96 10 1.29 1.17 1.07 0.96 0.94 12 1.27 1.17 1.08 0.95 0.92 14 1.26 1.18 1.06 0.97 0.90 2 1.46 1.36 1.27 1.16 1.04 4 1.42 1.37 1.28 1.17 1.07 6 1.45 1.39 1.25 1.15 1.05 30°并行 8 1.47 1.37 1.25 1.14 1.07 10 1.46 1.34 1.29 1.17 1.02 12 1.44 1.38 1.3 1.16 1.05 14 1.45 1.38 1.18 1.15 1.06 2 1.63 1.57 1.39 1.22 1.07 4 1.62 1.55 1.37 1.24 1.08 6 1.65 1.53 1.40 1.24 1.05 60°并行 8 1.67 1.58 1.35 1.27 1.09 10 1.60 1.48 1.38 1.21 1.08 12 1.64 1.50 1.38 1.20 1.10 14 1.63 1.55 1.37 1.22 1.06 2 1.73 1.64 1.44 1.20 1.12 4 1.77 1.66 1.45 1.18 1.11 6 1.87 1.67 1.47 1.25 1.13 90°叠落 8 1.63 1.52 1.42 1.22 1.12 10 1.58 1.52 1.41 1.23 1.17 12 1.69 1.55 1.46 1.17 1.13 14 1.83 1.57 1.42 1.22 1.14 易知直径变化率受间距影响不大,随埋深增加而减小,且变化幅度很小;单孔隧道直径变化率最大,而在与下层隧道呈90°叠落、60°、30°、0°并行的情况下,上层隧道的直径变化率依次减小。该变化规律与动力时程法模拟的上层隧道的地震动反应特性类似,见图 10。
(2)拟合理论模型选取
基于Matlab软件与上述理论分析,用拟合理论模型的Fourier函数来拟合直径变化率,效果良好。曲线见图 11。
$$f\left( x \right)={{a}_{0}}+{{a}_{1}}\ \cdot \ \cos \left( wx \right)+{{b}_{1}}\ \cdot \ \sin (wx)+{{a}_{2}}\ \cdot \ \cos (2wx)+{{b}_{2}}\ \cdot \ \sin (2wx)$$ (4) 式中,f(x)为动力时程法直径变化率;x为改进反应位移法直径变化率;a0、a1、b1、a2、b2、w为修正系数,其均值与取值范围分别为6.91×10-2(-5.83×10-2,1.97×10-1)、5.75×10-2(-2.46×10-2,1.40×10-1)、-2.30×10-2(-0.61,0.56)、1.76×10-2(-0.26,0.29)、1.51×10-2(-0.21,0.24)、32.39(-14.13,78.91)。
在拟合优度方面,拟合误差平方和(SSE)为2.78×10-3;置信区间(R-square)为0.97;调整置信区间(Adjusted R-square)为0.96。
4. 结语
(1)本文改进了反应位移法计算模型,通过在双线隧道之间设置夹土弹簧来反映土体-结构间的相互作用,提出了该弹簧刚度的计算公式。在此基础上,将上层隧道横断面的计算结果与准确的动力时程法模拟数据进行对比分析。
(2)在与下层隧道呈不同角度、间距、埋深的多种工况下,上层隧道横断面的抗震研究结果表明,两种方法所得的弯矩与轴力的最大值的变化规律基本一致:埋深对内力值的影响最大,间距次之;隧道水平并行时,上层隧道在同一埋深下,随着间距增大,内力值由先小于单孔隧道,然后增大到与单孔隧道近乎相同;隧道由斜穿变化到竖直下穿时,上层隧道在同一埋深下,随着间距增大,内力值由先大于单孔隧道,然后减小到与单孔隧道近乎相同,并且在4—6m时出现峰值,该范围是设计双线隧道的薄弱区域。同时,上层隧道直径最大变化率受间距影响不大,随埋深增加而减小,且隧道所呈角度越大变化率越大。
(3)弯矩最大值主要分布在上部拱腰两侧;轴力最小值分布在隧道顶部区域,最大值基本上分布在结构下部拱腰附近,集中现象比较明显。
(4)当隧道间距S≥8m时,内力最大值变化不大,并趋近于单孔隧道,这说明隧道间距S≥8m时,可忽略夹土层地基弹簧的变化对上层隧道内力值的影响,按现有的反应位移法进行计算。
(5)基于Matlab软件与准确的动力时程法,对改进模型的弯矩、轴力最大值和直径最大变化率分别采用Cubic Polynomial、Sum of Sin Functions、Fourier拟合理论模型进行修正来减少误差,拟合效果良好。
综上所述,本文提出的改进反应位移法较为合理,适用于小间距盾构隧道的抗震分析与设计。
-
表 1 破坏概率矩阵(单位:%)
Table 1. Damage probability matric(unit:%)
破坏等级 地震烈度 Ⅵ Ⅶ Ⅷ Ⅸ ≥Ⅹ 基本完好 88 74 55 35 15 轻微破坏 12 23 33 30.5 20.5 中等破坏 0 2 10.3 1.5 0.2 严重破坏 0 0 1.5 7.5 1.5 完全破坏 0 0 0.2 1.5 7.5 表 2 不同易损性分析方法优缺点的对比
Table 2. Comparison of different methods of seismic fragility analysis
方法 优点 缺点 震害调查法 简单易行;可呈现实际震害情况;可展示结构实际易损性 依赖专家经验,主观性强;数据难以获得;具有区域性和结构专有性;观察结果主观性强 解析法 可信度高;可考虑多数不确定性;客观性强 计算量大;需求模型中分散度较大 混合法 包含震害资料;计算量小 耗时且计算费用高;需要选择适当的解析方法;定义破坏状态有主观性 实验法 可提供真实破坏条件;可信度高 实验费用高且耗时 -
党王桢, 梁兴文, 2015.RC框架结构地震倒塌易损性分析的实用方法.震灾防御技术, 10(4):986-995. http://zzfy.eq-j.cn/zzfyjs/ch/reader/view_abstract.aspx?flag=1&file_no=20150417&journal_id=zzfyjs 冯耀华, 2017.基于变形和能量的木结构古建筑易损性分析.西安:西安建筑科技大学. 高志刚, 许金余, 白二雷等, 2013.灰色统计法在建筑物(群)地震易损性分析中的应用.四川建筑科学研究, 39(4):205-208, 242. doi: 10.3969/j.issn.1008-1933.2013.04.047 韩建平, 魏世龙, 张鑫, 2018.基于谱匹配的近远场地震动作用下RC框架结构易损性分析.工程科学与技术, 50(3):149-157. http://d.old.wanfangdata.com.cn/Periodical/scdxxb-gckx201803018 何志明, 陈清军, 2018.考虑竖向地震作用效应的地下典型空间结构易损性性能指标探讨.力学季刊, 39(1):117-125. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=lxjk201801011 李静, 陈健云, 温瑞智, 2012.框架结构群体震害易损性快速评估研究.振动与冲击, 31(7):99-103. doi: 10.3969/j.issn.1000-3835.2012.07.022 林庆利, 林均岐, 刘金龙, 2017.汶川地震公路桥梁易损性研究.振动与冲击, 36(4):110-118. http://d.old.wanfangdata.com.cn/Periodical/zdycj201704018 刘如山, 胡少卿, 邬玉斌等, 2009.基于地震动参数的结构易损性表达方法研究.地震工程与工程振动, 29(6):102-107. http://d.old.wanfangdata.com.cn/Periodical/dzgcygczd200906012 乔亚玲, 闫维明, 郭小东, 2005.建筑物易损性分析计算系统.工程抗震与加固改造, 27(4):75-79. doi: 10.3969/j.issn.1002-8412.2005.04.019 邱舒睿, 高惠瑛, 2015.青海省农居地震灾害易损性研究.震灾防御技术, 10(4):969-978. http://zzfy.eq-j.cn/zzfyjs/ch/reader/view_abstract.aspx?flag=1&file_no=20150415&journal_id=zzfyjs 苏亮, 向畅, 蒋亦庞, 2017.基于模糊相似理论的地震易损性分析方法.华中科技大学学报(自然科学版), 45(8):74-80. http://d.old.wanfangdata.com.cn/Periodical/hzlgdxxb201708015 孙柏涛, 张博, 闫培雷, 2018.基于IDA方法的高层剪力墙结构地震易损性分析.山西建筑, 44(18):27-29. doi: 10.3969/j.issn.1009-6825.2018.18.015 王玉梅, 王仁杰, 王涛, 2017.设防框架的梁铰模型和能力易损性的试验研究.工程力学, 34(10):168-177. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=gclx201710019 徐超, 温增平, 2017.不同设防标准RC框架结构基于易损性分析的抗震性能评估.震灾防御技术, 12(4):845-857. http://zzfy.eq-j.cn/zzfyjs/ch/reader/view_abstract.aspx?flag=1&file_no=20170413&journal_id=zzfyjs 杨硕, 2016.非住宅类RC框架结构地震易损性分析混合方法研究.哈尔滨:中国地震局工程力学研究所. 姚新强, 2016.天津农居易损性与抗震能力分布研究.哈尔滨:中国地震局工程力学研究所. 尹之潜, 1996.结构易损性分类和未来地震灾害估计.中国地震, 12(1):49-55. http://www.cnki.com.cn/Article/CJFDTOTAL-ZGZD601.006.htm 于晓辉, 吕大刚, 王光远, 2008.土木工程结构地震易损性分析的研究进展.见: 第二届结构工程新进展国际论坛论文集.大连: 同济大学, 中国建筑工业出版社, 61-62. 于晓辉, 吕大刚, 2016.基于云图-条带法的概率地震需求分析与地震易损性分析.工程力学, 33(6):68-76. 张桂欣, 孙柏涛, 陈相兆等, 2018.北京市建筑抗震能力分类及地震灾害风险分析.地震工程与工程振动, 38(3):223-229. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dzgcygczd201803026 Braga F., Dolce M., Liberatore D., 1982.A statistical study on damaged buildings and an ensuing review of the MSK-76 scale.In: Proceedings of the Seventh European Conference on Earthquake Engineering.Athens, Greece: ECEE, 431-450 Casciati F., Faravelli L., 1989.Seismic vulnerability via knowledge based expert systems.In: Brebbia CA, ed., Structural Repair and Maintenance of Historical Buildings.Southampton: Computational Mechanics Pub., 299-307. Giannini R., Pagnoni T., Pinto P.E., et al., 1996.Risk analysis of a medieval tower before and after strengthening.Structural Safety, 18(2-3):81-100. doi: 10.1016/0167-4730(96)00007-0 Giovinazzi S., Lagomarsino S., 2004.A macroseismic method for the vulnerability assessment of buildings.In: Proceedings of the 13th World Conference on Earthquake Engineering.Vancouver: WCEE. Giovinazzi S., 2005.The vulnerability assessment and the damage scenario in seismic risk analysis.Braunschweig: Technical University Carolo-Wilhelmina at Braunschweig. Kiremidjian A.S., 1985.Subjective probabilities for earthquake damage and loss.Structural Safety, 2(4):309-317. doi: 10.1016/0167-4730(85)90017-7 Leon F., Atanasiu G.M., Gâlea D., 2006.Using data mining techniques for the management of seismic vulnerability.Key Engineering Materials, 326-328:501-504. doi: 10.4028/www.scientific.net/KEM.326-328 Lind N C., 1995.A measure of vulnerability and damage tolerance.Reliability Engineering & System Safety, 48(1):1-6. http://cn.bing.com/academic/profile?id=764f2dd68db29210e25ba3e53de04c97&encoded=0&v=paper_preview&mkt=zh-cn Lu Z., Yu Y., Woodman N.J., et al., 1999.A theory of structural vulnerability.The Structural Engineer, 18(77):17-24. http://d.old.wanfangdata.com.cn/OAPaper/oai_arXiv.org_1008.1770 MöllerO., Foschi R.O., Rubinstein M., et al., 2010.Estimating structural seismic vulnerability:an approach using response neural networks.Structure and Infrastructure Engineering, 6(1-2):63-75. doi: 10.1080/15732470802663797 Rossetto T., Elnashai A., 2003.Derivation of vulnerability functions for European-type RC structures based on observational data.Engineering Structures, 25(10):1241-1263. doi: 10.1016/S0141-0296(03)00060-9 SpenceR.J.S., Coburn A.W., Pomonis A., 1992.Correlation of ground motion with building damage: The definition of a new damage-based seismic intensity scale.In: Proceedings of the 10th World Conference on Earthquake Engineering.Rotterdam: WCEE, 551-556. Whimtan R.V., Reed J.W., Hong S.T., 1973.Earthquake damage probability matrices.In: Proceedings of the 5th World Conference on Earthquake Engineering.Rome: WCEE, 2531-2540. Yamaguchi N., Yamazaki F., 2002.Fragility curves for buildings in Japan based on damage surveys after the 1995 Kobe earthquake.In: Proceedings of the 12th European Conference on Earthquake Engineering.Auckland: ECEE, 2451-2458. -