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.
-
引言
2016年5月,习近平总书记在“科技三会”上明确提出科技创新、科学普及“一体两翼”的重要论断,将科普工作提升到前所未有的高度(任福君,2019)。同时,我国作为地震多发、灾害严重的国家(薄景山等,2019),开展地震科普工作对提升公众防震减灾意识和科学素养,促进经济社会发展,实现重大地震灾害事件的主动防御、科学避灾、有效减灾具有重要意义(蓝姝等,2017)。
开展公众问卷调查与分析,既是明确地震科普工作方向的重要依据,又是提升地震科普工作实效的重要途径。张吕等(2019)对云南地震科普工作开展了问卷调查,并提出相应的对策建议;魏本勇等(2013)针对北京市公众地震灾害认知开展了调查,并对其特点进行初步分析;刘子一等(2015)以初中生为例,针对上海市中学生地震科普工作现状开展了调查研究;张文佳等(2014)以江西省萍乡地区为例,针对少震弱震区公众防震减灾意识现状开展了调查分析;王林等(2015)基于媒体开展地震科普公众认知度专题调查研究。
天津市是全国地震重点监视防御区,也是我国唯一一个遭受过Ⅷ度地震破坏的特大城市。天津市历年来高度重视地震科普工作,特别是自2019年在全国率先推行全域地震科普的新模式以来,创新开展了一系列卓有成效的工作,促进社会公众防震减灾意识与地震科学素养的普遍提高,在地震科普工作方面反映出一定典型性和代表性,同时具备较好的社会基础。基于此,本研究以天津市为例,对该市公众地震科普认知与需求开展问卷调查和相应研究,旨在为提高地震科普工作科学性、针对性和有效性提供参考。
1. 研究方法
1.1 调查问卷的设计
本研究以公众地震科普认知与需求为主题,邀请天津市地震科普教育、地震灾害防治等领域专家,自行设计调查问卷。参考邹文卫等(2011)对北京市社会公众地震科普认知与需求的研究成果,并结合与专家的深度访谈,将公众地震科普认知与需求划分为2个要素:①公众对地震科学知识的认知与需求,主要体现为公众对地震监测预报、震害防御、应急救援、政策法规等领域知识的掌握程度与关注程度;②公众对地震科普工作的认知与需求,主要体现为公众对地震科普工作内容、形式、渠道等方面的满意度与行为偏好。问卷内容主要包括以下部分:
(1)公众对地震科学知识的认知现状。开展地震科普工作首先应准确把握公众认知现状,并针对公众认知的薄弱环节开展有针对性的工作。为此,在整体调查公众地震科学知识普及程度的基础上,从公众对地震灾害的恐惧心理、对地震自救互救技能的认知、对避震知识的认知、对地震预报知识的认知、对地震谣传误传的认知、对地震逃生知识教育和演练的认知6个重点领域综合调查公众认知水平。
(2)公众对地震科学知识的需求情况。公众需求是决定地震科普内容的重要依据,为此,调查了公众对地震科学知识的总体关注程度及对地震灾害报道、地震预报知识、地震科技知识、地震政策法规等方面的关注偏好。
(3)公众对地震科普工作的认知与评价。主要包括公众对地震科普工作的满意度、对地震科普通俗性和易获取性的评价、对地震科普作用的评价、对地震科普工作的改进建议等,其目的在于调查地震科普工作的有效性与实用性。
(4)公众对地震科普工作的需求与偏好。主要包括公众关注地震科普工作的目的、对地震科普内容的倾向性、偏好的地震科普形式与接收渠道等,其目的在于从受众的角度调查地震科普服务的精准化与个性化。
(5)调查者的基本情况。主要包括受访者性别、年龄、受教育程度、职业与地域等,其目的在于分析不同人群的认知与需求差异,从而便于开展不同人群的细分研究。
1.2 调查方法
1.2.1 调查对象
为保证调查的多样性与全面性,本研究根据年龄、学历、职业等不同,选取天津市机关单位工作人员、企业员工、事业单位工作人员、大学生、中小学生、社区居民、农民、驻津部队、残障人士等特殊人群作为调查对象,可较充分地反映天津市不同群体地震科普认知与需求现状。
1.2.2 抽样方法
本研究遵循随机抽样的原则,选取天津市防震减灾科普示范学校、地震安全示范社区、有关政府部门和企事业单位为调查点,采取现场随机调查与入户访问调查的方式,最终现场回收1517份问卷。其中有效问卷1500份,有效问卷率达98.88%,有效问卷数量满足99%的可信度和5%的抽样误差条件。
1.2.3 数据分析
所有问卷回收、审查、编码后,使用SPSS21.0统计软件对数据进行处理与分析,并采用描述性统计方法,根据样本性别、年龄、文化程度、职业等不同特征,对不同人群的数据进行差异性统计与分析。
2. 调查结果
2.1 调查样本信息
有效调查问卷分析结果表明(表 1),受访者性别比例基本相同,其中,男性占49.26%,女性占50.74%,女性占比略高于男性;受访者主要为中青年人,其中,25岁以下占25.80%,25—35岁占28.73%,36—45岁占25.20%,45岁以上占20.27%;受访者文化程度主要集中于初中、高中、大学专科和大学本科,共占89.30%;受访者以事业单位人员、社区居民和企业员工居多,分别占20.00%、20.00%和16.67%。
表 1 调查样本基本情况Table 1. Basic data sheet of interviewees in the study性别 年龄 文化程度 职业 构成 占比/% 构成 占比/% 构成 占比/% 构成 占比/% 男 49.26 <18 10.53 小学 6.73 社区居民 20.00 女 50.74 18—25 15.27 初中 18.57 企业员工 16.67 26—35 28.73 高中 19.78 中小学生 10.00 36—45 25.20 大学专科 10.77 大学生 10.00 46—55 11.53 大学本科 40.17 事业单位 20.00 56—65 6.67 硕士及以上 3.97 特殊人群 6.67 ≥66 2.07 农民 6.67 机关 6.67 驻津部队 3.33 2.2 公众对地震科学知识的认知现状
(1)公众地震科学知识的普及程度(图 1)。整体来看,天津市公众地震科学知识普及程度较高,总体普及率达72.75%,但群体性差异较大,其中中小学生占比最高,达84.00%;驻津部队、企事业单位和政府机关工作人员次之,均超过60%;社区居民和农民普及率为50%—60%;残障人士等特殊人群普及率最低,仅为35.00%,与其他群体存在较大差距。
(2)公众对地震灾害的恐惧心理。调查结果显示,42.99%的公众对地震灾害存在恐惧心理。大多数人认为所居住的房屋建筑处于中等设防水平,对居住环境安全缺乏足够的信心,且普遍认为办公场所抗震设防等级略高于居住的房屋建筑。
(3)公众对地震自救互救技能的认知。据统计,汶川地震20min内被救出的伤员存活率为98%,30min内为90%,1h内为63%,48h内为36%,这说明如果公众能掌握一定自救互救技能,可大大降低地震灾害导致的人员伤亡。本次调查中仅有54.93%的公众表示掌握自救互救知识,仅有4.4%的公众了解应急包内应放有逃生设备和应急照明设备。
(4)公众对避震知识的认知。整体来看,公众对自身震后避震能力的评价为中上等(3.59分,满分5分),83.33%的公众能正确选择避震方式。具体而言,对于震后家庭避震空间的选择,选择桌底为避震空间的占71.27%,选择小开间房屋为避震空间的占43.13%,选择桌子旁边为避震空间的仅占12.13%,14.40%的公众错误地选择将阳台作为安全的避震空间。对于室外避震空间的选择,13.73%的公众错误地选择迅速返回室内,3.27%的公众错误地选择抱住电线杆。
(5)公众对地震预报的认知。13.52%的公众认为应了解地震预报的复杂性;64.32%的公众认为地震部门应加强地震预报工作,提高预警预报工作的及时性;22.16%的公众认为应大力普及地震预报知识。另外,调查结果显示,公众对地震预报发布权限的认知不清,误认为地震部门有发布地震预报权限的比例高达65.38%。
(6)公众对地震谣传误传的认知。28.07%的公众认为地震谣传是无稽之谈,不足为信;33.00%的公众根据传言认为存在地震发生的可能;27.07%的公众将传言记在心里,并为可能发生的地震做准备;23.67%的公众表示会恐慌,周围的人都信了,应去安全的地方躲避。
(7)公众对地震逃生知识教育与演练的认知。调查结果显示,大部分公众能正确认识开展地震逃生教育与演练的必要性和重要性,48.63%的公众认为开展防震避震演练非常有必要,47.56%的公众认为有必要,但仍有3.82%的公众认为无所谓。
2.3 公众对地震科学知识的需求情况
(1)公众对地震科学知识的关注度。调查结果显示,10.37%的公众经常关注地震科学知识,54.11%的公众偶尔关注,29.77%的公众很少关注,5.75%的公众从未关注。调查数据表明,由于地震灾害相对于气象、洪水等灾害而言,发生概率较小,导致公众出现了麻痹心理,对地震科学知识的关注度相对较低。
(2)公众对地震科学知识的关注点(图 2)。调查结果显示,公众对国内外重大地震灾害的新闻报道关注度最高,分别为55.53%和48.00%;公众对地震预报知识、地震科技知识的关注度次之,分别为39.67%和35.33%;公众对震灾防御和自救互救、地震领域政策法规的关注度最低,分别为28.00%和21.13%。
2.4 公众对地震科普工作的认知与评价
(1)公众对当前地震科普工作的满意度(图 3)。总体而言,公众对当前地震科普工作持满意态度,16.61%的公众表示非常满意,31.61%的公众表示满意,32.55%的公众表示一般,19.22%的公众表示不太满意或非常不满意。具体来看,公众对地震科普作品数量的满意度为3.44分(满分5分,下同),对地震科普作品质量的满意度为3.61分,对地震科普作品内容生动性和丰富性的满意度为3.65分。
(2)公众对地震科普工作通俗性与易获取性的评价。在通俗性方面,30.57%的公众表示均能看懂,49.26%的公众表示大部分能看懂,18.69%的公众表示仅能看懂一部分,另有1.48%的公众表示基本看不懂。考虑受访者在文化程度、职业、年龄等方面的差异,目前地震科普工作在通俗性方面总体处于中等偏上水平。在易获取性方面,24.58%的公众表示基本均能获取,52.59%的公众表示仅能获取一部分,18.48%的公众表示仅能获取较少的内容,另有4.38%的公众表示无法获取。
(3)公众对地震科普工作作用的评价。调查结果显示,28.50%的公众认为地震科普工作非常有帮助,59.73%的公众认为有帮助。具体而言,公众认为最有必要宣传的是防震避震自救知识,其次是地震领域政策法规、政府有关地震灾害的应急预案与措施、震情灾情信息等。
(4)公众对改进地震科普工作的建议。公众认为地震科普工作应提高参与性、互动性与趣味性,占37.27%;其次为应及时更新地震科学知识,占33.80%;另外,还有33.27%的公众认为应使地震科普内容更贴近社会公众,增强地震科学知识对公众生活、出行的具体指导。
2.5 公众对地震科普工作的需求与偏好
(1)公众关注地震科普工作的目的。调查结果显示,公众关注地震科普工作主要考虑保护自身生命财产安全,62.07%的公众为保护自己与家人关注地震科普工作。此外,45.13%的公众为扩展自己的知识面关注地震科普工作。
(2)公众对地震科普内容的倾向性。最受公众喜爱的地震科普内容为防灾减灾自救知识,占63.93%;其次为地震预报知识方面的内容,占44.27%。另外,对于地震灾害预警的具体含义、地震监测仪器常识、地震应急预案与措施、地震名词术语、地震成因分析等方面的知识,均有20%以上的公众表示喜爱。
(3)公众对地震科普形式与接受渠道的诉求。调查结果显示,公众最喜爱的科普形式为科普视频短片和电视科普节目,占37.13%;其次为网站科普文章、微信和微博,占比分别为36.93%、29.40%。在地震科普接受渠道方面,电视仍是公众的首选渠道,占53.87%,其次为网站和微信,分别占47.47%、37.73%。
3. 结论与建议
通过对天津市公众的调查与分析,得出以下结论:
(1)公众地震科学知识普及率总体较高,但不同人群之间存在较大的差异性;
(2)公众对地震科学知识的掌握程度不够全面,在自救互救技能、地震预报知识、地震谣传辨识等方面存在明显短板;
(3)公众对地震科普工作存在较高期望,公众主动学习地震科学知识的目的主要为保护自身与家人的安全;
(4)公众对目前地震科普工作的满意度整体处于中等偏上水平,但在参与性、互动性与趣味性等方面仍存在较大的提升空间,在科普形式与渠道方面应进一步优化。
基于上述结论,提出以下工作建议:
(1)加强重点宣传,提高地震科普工作的针对性。以地震应急、地震预报、地震预警等为重点,策划开展专题科普宣传活动,促进公众地震科学素养的整体提升。尤其是目前公众对地震谣传的辨识能力普遍偏低,易造成社会恐慌,给经济社会发展带来负面影响,这应作为地震科普工作的突出重点。
(2)顺应公众需求,提高地震科普工作的实效性。进一步突出地震安全的主题,大力弘扬“生命至上、安全第一”的先进理念,特别是须加大对公众地震安全知识的普及,使地震科普工作更贴近群众生活。围绕公众最关注的知识领域,加大科普作品创作力度,制作公众喜闻乐见、通俗易懂的科普作品。通过举办社会性的地震科普作品创作大赛、地震科学知识竞赛、新媒体线上线下互动等,广泛动员和邀请社会公众参与到地震科普工作中,提高地震科普的互动性和参与性。
(3)针对弱势群体,提升地震科普服务的个性化水平。针对目前地震科学知识认知水平较低的人群,制定专项科普方案与措施,创新方式、补齐短板。本次调查中接受地震科普占比最低的为残障人士等特殊群体,建议有关部门与街道、社区对接,建立残障人士等特殊人群信息库,采取入户发放宣传资料或以社区为单位组织小型科普宣传讲座的方式,加强对残障人士等特殊人群的地震科学知识普及。
(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. -