Seismic Vulnerability Analysis of Masonry Structures with Uncertain Parameters Based on DOE
-
摘要: 砌体结构参数具有离散性高且其受到地震作用后响应结果非线性强等特点,所以在砌体地震易损性研究中结构参数不确定性对结果产生的影响不容忽视。针对汶川地震中具有代表性的砌体结构,提出基于DOE(试验设计)方法考虑在砌体结构地震易损性研究中结构参数不确定性的影响。首先,使用Plackett-Burman方法开展结构参数的灵敏度分析,筛选出对砌体结构地震响应影响较大的3个结构参数;然后,根据筛选结果及地震参数PGA进行试验设计,进而建立结构参数与最大层间位移角的响应面回归模型;最后,通过进行蒙特卡罗模拟获得地震易损性曲线,并进一步评估结构参数不确定性对砌体地震易损性分析的影响程度。研究结果表明,弹性模量、密度和抗拉强度是对砌体结构易损性影响较大的3个参数;当地震动PGA为0.1 g时,相较于其他结构参数,弹性模量对结构地震响应的影响最显著;当PGA为0.2 g时,3个结构参数的影响由大到小依次为弹性模量、密度和抗拉强度,当PGA为0.4 g时,3个结构参数的影响由大到小依次为弹性模量、抗拉强度与密度。本文所提出的方法计算量少、精度高,可为有效解决砌体结构参数种类众多和材料非线性强等难题提供新思路。Abstract: Given the high variability of masonry structural parameters and the strong nonlinearity of the seismic response after earthquake activity, the uncertainty of structural parameters plays a significant role in the seismic vulnerability assessment of masonry structures. This study proposes the use of a Design of Experiments (DOE) method to account for the influence of structural parameter uncertainty, with a focus on masonry structures affected by the Wenchuan earthquake. First, the Plackett-Burman method was employed to perform a sensitivity analysis of structural parameters, identifying the three parameters that most significantly affect the seismic response of masonry structures. Based on these results and the peak ground acceleration (PGA) as the seismic parameter, the DOE was conducted, leading to the development of a response surface regression model between the structural parameters and the maximum inter-story drift angle. Finally, a Monte Carlo simulation was performed to derive the seismic vulnerability curve, further assessing the impact of structural parameter uncertainty on the seismic vulnerability analysis of masonry structures. The results indicate that: (1) the elastic modulus, masonry density, and tensile strength are the three parameters with the greatest impact on the seismic vulnerability of masonry structures; (2) when the local PGA is 0.1g, the elastic modulus has the most significant influence on the seismic response compared to other parameters. At a PGA of 0.2 g, the influence order is elastic modulus, masonry density, and tensile strength, while at a PGA of 0.4 g, the order changes to elastic modulus, masonry tensile strength, and masonry density. The proposed method in this study has low computational complexity and high accuracy, offering an effective approach for addressing the numerous structural parameters and strong material nonlinearity in masonry structures. This approach provides new ideas in improving seismic vulnerability assessments for such structures.
-
引言
应急演练是指各级人民政府及其部门企业事业单位、社会团体等组织相关单位及人员,依据有关应急预案,模拟应对灾害等突发事件的活动,对以角色扮演、实景模拟或计算机仿真模拟发生的突发事件进行响应和处置,以达到检验应急预案和协调机制、训练专业队伍、开展宣传教育等目的。
2003年应对SARS灾害以来,我国应急管理逐步向现代化发展。南方雨雪冰冻灾害、汶川8.0级地震、“12.31”上海踩踏事件、“8.12”天津港爆炸事故等重特大自然灾害和事故灾难的应对实践,推进了包括预案制定和演练等各项应急相关工作的开展,各级政府、企业事业单位、学校、社区等都开展了应对地震等自然灾害和安全生产事故的多类型应急演练,在一定程度上起到检验预案、完善准备、锻炼队伍、增强意识、提升能力的作用,应急演练已成为应急准备工作的重要环节和例行工作(张媛等,2014)。
社区是社会公共安全治理的基本单元。地震灾害应对实践表明,社区等基层组织在抢救生命、安置群众生活等方面发挥了重要作用。国家综合减灾救灾体制机制改革、国务院应急管理机构改革均要求推进基层自救互救能力和应急管理能力建设,加强对基层民众的自救互救知识培训和科普宣传教育。社区应急演练是提升灾害应对处置能力、开展科普宣传教育并有效提升民众意识的重要途径。近年来,很多社区都开展了针对地震等自然灾害的应急演练,主要内容包括人员疏散、自救互救、人员安置、次生灾害处置等。演练活动通常以实战形式开展,达到了锻炼基层队伍,提升基层组织、社区公众应急处置和防灾减灾意识的目的(张勤等,2014)。但社区地震应急演练仍存在很多问题,在发挥实效性方面存在很大差距,具体表现为:演练方案规划及设计专业程度不足,缺少地震应急领域相关专业的指导;由于涉及人员较多,组织耗时长、难度大,使开展演练的社区数量有限;主要以表演、宣讲的方式开展,针对社区应对灾害组织协调能力等的演练较少,不利于社区应急管理能力的提高;缺乏后续总结、对比、评价等环节,降低了演练实效性;由于社区组织的多样性和复杂性及所处环境的具体性和特殊性,社区应急演练缺乏实战性检验等。
针对上述问题,本文结合传统应急桌面演练及实战演练优势,提出基于在线平台的互联网+社区地震应急桌面演练新模式,并在青海省海西蒙古族藏族自治州(海西州)开展了试点实践,在提升演练专业性、扩大演练社会参与度、提升社区应对地震灾害的组织协调能力等方面取得了初步成效。目前,该模式已成功应用于青海省、江苏省、广州市、徐州市及海西州等20多个省市地区。
1. 互联网+社区地震应急桌面演练平台
互联网+社区地震应急桌面演练模式可实现多地同步联动演练,组织更省时省力,方案设计更灵活多样,由于具有一定的竞技性质,较传统桌面演练更贴近实战,同时该模式可进行全过程记录、评估,大大增强了演练实效性。应急桌面演练平台是在线云平台(图 1),可为应急桌面演练提供零部署、即开即用的实时应急演练服务。该平台基于先进的云服务理念,提供SaaS模式应急演练软件服务,具有较强的伸缩性和灵活性,支持计算机、Pad、手机等多终端应用,并可离线扩展,实现互联网+演练模式。该平台支持包括桌面演练、功能演练、综合演练在内的各种类型演练策划、设计、实施和评估,提供丰富的情景资源库、专业的演练脚本设计工具、全方位的推演过程控制、综合的演练评估工具和评估报告、三维可视化的互动体验,使演练组织者、参与者、观摩人员、评估专家等的深度参与成为可能,有助于提升演练的科学性和有效性,实现演练活动目标(仝鹏,2014)。
2. 海西州演练设计与实施
2018年5月15日至16日,海西州及其下辖市县73个社区管理人员、志愿者队伍负责人及社区民警等180余人采用互联网+社区地震应急桌面演练模式,同步开展不同地域桌面演练活动。演练分为4个时段,连续进行4场。由于演练具有分布式开展、可自行组织参加的便利,海西州各市、县(行委)参演组可根据各自时间安排参加演练。
本次演练定位为使社区相关管理人员进一步熟悉地震应急预案,掌握地震应急工作流程,加强突发事件应变处置及协调指挥能力,增强各种信息分析判断和决策能力,从而使地震等突发事件发生后,社区居民在管理人员的带领和引导下有序开展相关处置工作,有效应对灾害事件。桌面演练准备阶段主要包括演练方案制定、脚本编写、媒体信息制作、参演人员分组、培训等,具体流程如图 2所示。
2.1 演练方案编制
2.1.1 设计演练场景
演练场景设计是对假想事件按其发生过程进行的叙述性说明,即针对假想事件发展过程,设计出一系列突发和次生、衍生事件,让参演人员在演练过程中身临其境,对情景事件变化做出真实的应急反应。
本次演练对象主要为社区管理层及志愿者队伍负责人,为加强对地震灾害的体验和了解,场景设计采用了具有一定破坏性及人员伤亡的灾害背景。设定地震震级为6.5级,震源深度为10km,海西州普遍震感强烈,临近震中区域的城乡结合部发生少量建筑物或构件破坏,并出现人员伤亡情况。部分市政道路、线塔、管道等基础设施受到不同程度影响,局部区域出现供电、供水、供气中断及通信不畅的现象,震中区域附近的化工厂存在危险品泄漏风险。同时,此次地震还会引发大量社会舆论,备受关注,社会稳定面临严峻考验。
2.1.2 演练脚本编写
根据预设地震灾害背景,对地震发生后社区可能出现的真实情况进行设定,包括社区建(构)筑物破坏、通信不畅、物资短缺、人员恐慌、交通堵塞、毒气泄漏、人员伤亡、救援物资及装备短缺等,同时匹配相应的视频、图片、画外音等媒体资源辅助演练(李亦纲等,2007)。
演练脚本主要包括演练信息注入对象、演练事件类别、信息注入时间、事件信息来源及标题、事件信息发生地点、事件信息描述及预期应对方式,具体脚本格式如图 3所示。
脚本信息包含的“预期应对方式”是演练方案制定单位以演练脚本为依据,同时结合灾害处置案例及实际工作制定的用于参考的应对答案,“预期应对方式”并不是标准答案,仅供演练后期评估参考使用。
2.2 组织演练培训
培训是桌面演练必不可少的环节,除需对桌面演练整体日程及程序进行介绍外,还需对桌面演练进行详细说明,包括桌面演练的目的意义、演练形式及演练操作系统等,同时需对演练系统进行操作练习,并详细介绍演练所用设施设备用途及使用方法,使所有参演人员提前融入演练场景,正式演练开始后可快速进入角色,全情投入到演练中,令演练效果达到最佳。另外,在培训过程中还需重点强调演练纪律及注意事项,确保演练顺利进行。具体培训内容如图 4所示。
海西州社区地震应急桌面演练由于参演地区较分散,参演人员较多,参演时间不同,无法开展集中式培训,因此采用分散培训的方式。首先,制定桌面演练详细流程图及参演说明,下发至所有参演单位及相关人员,各单位分别开展演练前培训工作,所有技术支持由海西州地震局负责;桌面演练正式开始后,后参演的单位可在演练系统开放时进行操作练习;最后,对于个别有培训需求的参演单位通过在线方式由海西州地震局进行有针对性的培训。
2.3 具体演练流程
桌面演练地点确定并布置完成后,即可按照演练场次安排陆续开展演练。正式演练开始前需先进行身份注册,注册内容包括参演者所在区域、社区名称、社区任职或主要工作及电子邮箱等联系方式。具体演练流程如图 5所示。桌面演练正式开始后,演练系统会定时自动推出事件信息,每条信息推出后所有参演者进行充分讨论,讨论过程中须考虑社区实际情况及参演者自身工作职责。每条事件信息涉及某方面的灾害类型,处置过程中除能了解该灾害类型对社区的影响程度外,还可发现本社区存在的不足及薄弱环节,在制定灾害事件响应措施的同时也可对本社区各项工作方案的制定进行完善。部分参演组由多个社区组成,因此可在讨论中相互学习借鉴。
桌面演练结束后,可在最后环节提交本次演练参与感受及意见、建议,演练评估报告以邮件的方式发给各参演组,作为各社区后续评估总结资料,并永久保留。
3. 演练评估与数据分析
桌面演练数据可在演练全部结束后导出。海西州社区地震应急桌面演练共导出社区地震应急桌面演练应对及评估记录有效数据21套,数据内容主要包括演练场次、时间、演练分组、标题、演练事件内容及应对措施等。数据导出后,针对数据开展分析评估,结合本次演练整体情况完成评估报告,评估报告主要包括本次演练背景、开展情况及评估意见,另附评估记录表,评估记录表(部分)格式如表 1所示。
表 1 海西州社区地震应急桌面演练评估记录(部分)Table 1. Evaluation list of the tabletop exercise in haixi prefecture (Partial)序号 事件标题 事件内容 应对措施 评估结果 1 地震发生 地震发生 及时联系受困人员,进行评估,并上报请求支援,同时进行自救 合格 2 地震造成破坏 某年某月某日,某市发生强烈地震,某社区震感强烈,居民纷纷离开居民楼避震,在社区广场、公园等处聚集,有居民报告人员避震过程中被砸伤。社区部分居民楼等建筑出现不同程度的破坏,有建筑物构件、震碎的玻璃掉落,并发生砸坏车辆、砸伤人员的现象。地震发生在夏天,天气炎热 组织自救互救,保持正常秩序,了解伤亡情况,第一时间上报 优秀 3 楼梯间门扇卡死 有居民报告某栋楼楼梯间门框变形,门扇卡死无法打开,居民无法从楼梯间撤离,亟需帮助 对情况进行评估,同时上报进行救援,了解情况先组织自救 优秀 4 社区居民恐慌 在社区避难场所聚集的居民数量逐渐增加,由于天气炎热,部分老人、孩子出现身体不适,民众情绪出现波动,现场秩序面临失控 及时疏导,做好心理安慰,及时供给饮用水 优秀 本次桌面演练共21个参演组,事件信息应对数量共272条,其中应对优秀200条,占总量的74%,应对合格50条,占总量的18%,无效应对22条,占总量的8%,各组应对情况如图 6所示。
分析表明,此次海西州社区地震应急桌面演练所有参演社区完成了全部演练事件信息的应对任务,处置措施总体情况较好,体现出一定的地震灾害应对能力,但演练过程中也暴露了参演单位及参演者存在的不足,具体问题和改进建议如下:
(1)注意角色定位,社区发挥承上启下的过渡作用,应及时沟通上级和外部支援,专业问题寻求专业支持;
(2)进一步熟悉社区志愿者队伍及相关应急装备情况,未建立队伍的应及早建立,以便自救互救;
(3)明确社区是否制定针对地震灾害的应急预案或相关方案,未制定的尽快制定,并通过开展演练进行检验;
(4)进一步了解社区周边支持保障情况,如物资、视频、队伍支持等;
(5)注意震后社区次生灾害隐患的监控及排查等相关工作;
(6)提高科普宣传教育工作的重视程度,利用挂图等工具向居民普及地震现象、地震灾害、应急避险等知识;
(7)预期应对方式仅作为参考,具体处置措施应视社区实际情况而定。
最终的评估报告包括桌面演练背景、演练现场总体情况、演练总体评价、参演组评估统计表及演练应对与评估记录表。评估报告除包括专业分析评价外,还包括桌面演练原始数据及参演人员应对答案,该报告不仅可作为演练总结评估的依据,还可作为参演单位进行反馈对比的材料。根据评估报告,参演单位可对演练期间未能充分讨论或存有疑义的事件信息重新进行深入讨论,也可根据事件信息应对情况对本单位现有应急预案、工作流程等进行改进和完善,同时还可与其他社区进行横向对比,取长补短(郑通彦等,2014)。
4. 结论
通过开展海西州社区地震应急桌面演练活动,参演人员一致认为互联网+社区地震应急桌面演练模式形式新颖、流程清晰、情景设置针对性强、内容贴近实际,不仅充分体现了各市、县(行委)政府对突发事件应急管理和防震减灾工作的高度重视,更达到了锻炼队伍、检验预案、提升能力的目的,同时展示了社区工作人员良好的精神风貌和地震应急处置能力,对保障海西州公共安全起到了积极作用。
随着经济、社会的不断发展,社区功能不断完善,而社区的作用及需要承担的责任越来越大。为了能更好地应对类似地震的突发事件,加强应急管理意识、提高应急管理能力已成为所有社区必不可少的重要工作之一,而基于互联网的地震应急桌面演练也将成为提高社区地震应急处置能力的重要环节,并将在实践过程中不断发展、提升。
-
表 1 数值模型信息
Table 1. Numerical model information
建筑构件 单元类型 构件之间的相互作用 材料属性 悬挑梁 C3D8R Tie 设计规范混凝土本构模型 现浇楼板 C3D8R Tie 设计规范混凝土本构模型 钢筋 T3D2 Embedded region 双直线理想弹塑性模型 墙体 C3D8R Tie 整体式砌体结构本构模型 表 2 模型前2阶自振周期
Table 2. The first 2 vibration cycles of the model
振型 自振周期/s 振型特征 一阶 0.141 横向平动 二阶 0.115 纵向平动 表 3 结构参数及其分布
Table 3. Structural model parameters and their distribution
参数 平均值 变异系数 分布类型 砌体材料参数 阻尼比 0.05 0.3 正态 弹性模量/MPa 1 350 0.15 正态 抗压强度/kPa 2 448.29 0.17 对数正态 抗拉强度/kPa 235.938 0.2 对数正态 密度/(kg·m−3) 1 600 0.1 正态 混凝土材料参数 密度/(kg·m−3) 2 400 0.07 正态 结构几何参数 层高/m 3.3 0.05 — 墙厚/mm 240 0.05 — 表 4 输入参数及标准化
Table 4. Input parameters and standardization
输入参数 下限值 平均值 上限值 砌体阻尼比 输入值 0.035 0.05 0.065 标准化 −1 0 1 砌体弹性模量/MPa 输入值 1 147.5 1 350.0 1 552.5 标准化 −1 0 1 砌体抗压强度/kPa 输入值 2 032.08 2 448.29 2 864.50 标准化 −1 0 1 砌体抗拉强度/kPa 输入值 188.750 235.938 283.126 标准化 −1 0 1 砌体材料密度/(kg·m−3) 输入值 1 440 1 600 1 760 标准化 −1 0 1 混凝土密度/(kg·m−3) 输入值 2 328 2 400 2 568 标准化 −1 0 1 层高/m 输入值 3.135 3.300 3.465 标准化 −1 0 1 墙厚/mm 输入值 228 240 252 标准化 −1 0 1 表 5 Plackett-Burman试验设计
Table 5. A case of Plackett-Burman test design
项目 序号 1 2 3 4 5 6 7 8 9 10 11 12 砌体阻尼比 −1 1 −1 −1 1 1 −1 1 −1 −1 1 1 砌体密度 −1 −1 1 −1 −1 1 −1 1 1 1 1 −1 混凝土密度 −1 −1 1 1 1 1 −1 −1 1 −1 −1 1 砌体抗拉强度 −1 1 1 1 −1 1 1 −1 −1 −1 1 −1 砌体弹性模量 −1 1 −1 −1 1 1 1 −1 1 1 −1 −1 墙厚 −1 1 1 −1 −1 −1 1 1 1 −1 −1 1 层高 −1 −1 1 −1 1 −1 1 −1 −1 1 1 1 砌体抗压强度 −1 −1 −1 1 −1 1 1 1 −1 1 −1 1 表 6 中心复合试验设计
Table 6. A case of central composite test design
运行序 砌体弹性模量 砌体密度 砌体抗拉强度 1 1 1 −1 2 −1 −1 −1 3 1 0 0 4 0 1 0 5 1 −1 −1 6 −1 1 −1 7 0 0 0 8 1 1 1 9 −1 −1 1 10 −1 1 1 11 −1 0 0 12 1 −1 1 13 0 −1 0 14 0 0 1 15 0 0 −1 表 7 中心复合设计输出结果
Table 7. Output results from the central composite design
运行序 最大层间位移角(×10−7) 0.05 g 0.1 g 0.15 g ... 0.35 g 0.4 g $ {\hat y_{u1}} $ $ {\hat y_{\sigma 1}} $ $ {\hat y_{u2}} $ $ {\hat y_{u2}} $ $ {\hat y_{u3}} $ $ {\hat y_{\sigma 3}} $ ... ... $ {\hat y_{u7}} $ $ {\hat y_{\sigma 7}} $ $ {\hat y_{u8}} $ $ {\hat y_{\sigma 8}} $ 1 3 949.9 2 209.4 8 740.6 3 417.6 35 780 10 027 ... ... 29 832 7 952.2 35 779 1 002 2 2 998.8 1 969.57 7 020.4 2 665.2 28 183 6 577 ... ... 24 316 6 219.3 28 182 6 577.2 3 4 499.3 1 736.1 9 200.1 3733.5 39 739 12 563 ... ... 33 441 9 797.4 39 739 12 563 ... ... ... ... ... ... ... ... ... ... ... ... ... 15 4 235.6 1 558.1 8 759.3 2 412.7 39 003 12 186 ... ... 32 472 9 523.8 39 003 12 186 表 8 响应面模型检验结果
Table 8. Verification of response surface model
响应面标号 $ {R^2} $/% $ R_{\mathrm{A}}^2 $/% 1 99.59 98.31 2 98.61 97.52 3 98.36 97.62 4 98.76 97.32 5 97.79 96.80 6 98.26 97.82 7 97.12 96.32 8 95.30 93.25 表 9 极限状态限值
Table 9. Limitation value of limited state
极限状态 轻微破坏(LS1) 中等破坏(LS2) 严重破坏(LS3) 毁坏(LS4) $ {\theta _{\max }} $ 1/2 000 1/1 600 1/700 1/350 -
陈力波,黄才贵,谷音,2018. 基于改进响应面法的公路简支梁桥地震易损性分析. 工程力学,35(4):208−218. doi: 10.6052/j.issn.1000-4750.2017.06.0512Chen L. B., Huang C. G., Gu Y., 2018. Seismic vulnerability analysis of simply supported highway bridges based on an improved response surface method. Engineering Mechanics, 35(4): 208−218. (in Chinese) doi: 10.6052/j.issn.1000-4750.2017.06.0512 黄琼,2014. 基于响应面法的钢筋混凝土框架结构有限元模型修正研究. 兰州:兰州理工大学.Huang Q., 2014. Finite element model updating of reinforced concrete frame structure based on response surface method. Lanzhou:Lanzhou University of Technology. (in Chinese) 李龙师,2020. 松原查干花震区农村房屋抗震性能研究及震害预测. 哈尔滨:中国地震局工程力学研究所.Li L. S., 2020. Seismic Capacity Analysis and Earthquake damage prediction of rural houses in Chaganhua earthquake area of Songyuan City. Harbin:Institute of Engineering Mechanics,China Earthquake Administration. (in Chinese) 吕大刚,刘洋,于晓辉,2019. 第二代基于性能地震工程中的地震易损性模型及正逆概率风险分析. 工程力学,36(9):1−11,24.Lü D. G., Liu Y., Yu X. H., 2019. Seismic fragility models and forward-backward probabilistic risk analysis in second-generation performance-based earthquake engineering. Engineering Mechanics, 36(9): 1−11,24. (in Chinese) 宋帅,钱永久,钱聪,2018. 桥梁地震需求中随机参数的重要性分析方法研究. 工程力学,35(3):106−114. doi: 10.6052/j.issn.1000-4750.2016.11.0856Song S., Qian Y. J., Qian C., 2018. Research on methods for importance analysis of random parameter in bridge seismic demand. Engineering Mechanics, 35(3): 106−114. (in Chinese) doi: 10.6052/j.issn.1000-4750.2016.11.0856 孙柏涛,张桂欣,2012. 汶川8.0级地震中各类建筑结构地震易损性统计分析. 土木工程学报,45(5):26−30.Sun B. T., Zhang G. X., 2012. Statistical analysis of the seismic vulnerability of various types of building structures in Wenchuan M 8.0 earthquake. China Civil Engineering Journal, 45(5): 26−30. (in Chinese) 熊立红,吴文博,孙悦,2012. 汶川地震作用下约束砌体房屋的抗震能力分析. 土木工程学报,45(S2):103−108.Xiong L. H., Wu W. B., Sun Y., 2012. Seismic performance of confined masonry buildings during the Wenchuan earthquake. China Civil Engineering Journal, 45(S2): 103−108. (in Chinese) 杨卫忠,2008. 砌体受压本构关系模型. 建筑结构,38(10):80−82.Yang W. Z., 2008. Constitutive relationship model for masonry materials in compression. Building Structure, 38(10): 80−82. (in Chinese) 杨玉成,杨柳,高云学等,1982. 现有多层砖房震害预测的方法及其可靠度. 地震工程与工程振动,2(3):75−86.Yang Y. C., Yang L., Gao Y. X., et al., 1982. Method of damage prediction for existing multi story brick buildings and its reliability. Earthquake Engineering and Engineering Vibration, 2(3): 75−86. (in Chinese) 叶继红,江力强,2021. 考虑多重不确定性的我国多层冷成型钢结构地震风险评估. 土木工程学报,54(2):74−83,126.Ye J. H., Jiang L. Q., 2021. Seismic risk assessment of mid-rise cold-formed steel structures in China considering various uncertainties. China Civil Engineering Journal, 54(2): 74−83,126. (in Chinese) 于晓辉,2012. 钢筋混凝土框架结构的概率地震易损性与风险分析. 哈尔滨:哈尔滨工业大学.Yu X. H., 2012. Probabilistic seismic fragility and risk analysis of reinforced concrete frame structures. Harbin:Harbin Institute of Technology. (in Chinese) 于晓辉,吕大刚,2012. 考虑结构不确定性的地震倒塌易损性分析. 建筑结构学报,33(10):8−14.Yu X. H., Lü D. G., 2012. Seismic collapse fragility analysis considering structural uncertainties. Journal of Building Structures, 33(10): 8−14. (in Chinese) 张令心,李孟达,刘洁平等,2018. 考虑结构不确定性的框剪结构地震易损性分析. 自然灾害学报,27(4):112−118.Zhang L. X., Li M. D., Liu J. P., et al., 2018. Seismic fragility analysis of frame-shear wall structures considering structural parameter uncertainty. Journal of Natural Disasters, 27(4): 112−118. (in Chinese) 张尚荣,谭平,杜永峰等,2014. 基于响应面法的层间隔震结构地震易损性分析. 振动与冲击,33(15):42−48.Zhang S. R., Tan P., Du Y. F., et al., 2014. Seismic fragility analysis of inter-story isolation structures based on response surface method. Journal of Vibration and Shock, 33(15): 42−48. (in Chinese) 周强,2012. 砌体结构抗震试验及弹塑性地震反应分析. 哈尔滨:哈尔滨工程大学.Zhou Q., 2012. Seismic test and elastic-plastic seismic response analysis of masonry structure. Harbin:Harbin Engineering University. (in Chinese) 周强,闵全环,熊拥军等,2018. 江西典型村镇既有房屋抗震性能分析. 自然灾害学报,27(1):96−105.Zhou Q., Min Q. H., Xiong Y. J., et al., 2018. Seismic performance analysis of existing typical rural buildings in Jiangxi. Journal of Natural Disasters, 27(1): 96−105. (in Chinese) 朱震宇,2015. 砌体墙抗震性能分析模型及其工程应用. 上海:上海交通大学.Zhu Z. Y., 2015. Analysis model of seismic behavior of masonry walls and its engineering application. Shanghai:Shanghai Jiao Tong University. (in Chinese) Choudhury T., Kaushik H. B., 2018. Seismic response sensitivity to uncertain variables in RC frames with infill walls. Journal of Structural Engineering, 144(10): 04018184. doi: 10.1061/(ASCE)ST.1943-541X.0002190 Demirel İ. O., 2010. A nonlinear equivalent frame model for displacement based analysis of unreinforced brick masonry buildings. Ankara: Middle East Technical University. FEMA, 2009. Quantification of building seismic performance factors. Washington: Federal Emergency Management Agency. Freeny A., 1988. Empirical model building and response surfaces. Technometrics, 30(2): 229−231. doi: 10.1080/00401706.1988.10488371 Gokkaya B. U., Baker J. W., Deierlein G. G., 2016. Quantifying the impacts of modeling uncertainties on the seismic drift demands and collapse risk of buildings with implications on seismic design checks. Earthquake Engineering & Structural Dynamics, 45(10): 1661−1683. Loukas Y. L., 2001. A Plackett-Burman screening design directs the efficient formulation of multicomponent DRV liposomes. Journal of Pharmaceutical and Biomedical Analysis, 26(2): 255−263. doi: 10.1016/S0731-7085(01)00419-8 Montgomery D. C., 2001. Design and analysis of experiments. 5th ed. New York: John Wiley & Sons, Inc. Parisi F., Augenti N., 2012. Uncertainty in seismic capacity of masonry buildings. Buildings, 2(3): 218−230. doi: 10.3390/buildings2030218 Towashiraporn P., 2004. Building seismic fragilities using response surface metamodels. Atlanta: Georgia Institute of Technology. -