摘要
准确且切实可行的损伤评估方法可对震损结构的损伤程度进行量化,并为结构修复加固提供有效指导。针对实际工程条件以及已有损伤评估或识别方法在精度和可行性方面的不足,在传统弹塑性耗能差损伤评估方法上提出通过剪力比值消除质量影响的广义损伤评估方法,仅以地震信号和结构动力响应即可实现对结构损伤的动态评估。在此基础上,提出用前震信号和Kalman滤波识别技术推演结构理想弹性响应,并利用主震信号实现结构整体和楼层损伤快速评估的方法。以钢混框架结构为例进行验证,结果表明:所提出的震损结构整体及层间损伤评估方法能够准确反映结构损伤的动态演变过程以及薄弱层的分布。与传统弹塑性耗能差方法以及基于频率变化率的损伤评估方法对比,所提方法具有良好的工程实用性,可为实际震损结构损伤快速评定提供有效支持。
在中震和大震下,建筑结构的正常使用功能和整体性能将明显下降,严重时将造成大量人员伤亡和重大经济财产损
目前震损结构的损伤分析方法的研究和应用主要从损伤评
近年来,从结构健康监测理念发展而来的结构损伤识别为解决上述问题提供了新的思路和手段。在损伤识别方法中,广义的结构损伤识别被定义为对与结构力学性能相关的静、动力特征发生改变进行定性和定量的判断。结构地震损伤识别则是根据实时或准实时监测数据追踪地震过程中结构动力性能和动态特征参数的变化,发现结构的损伤位置和程度,以评估地震过程中结构的损伤状态。目前损伤识别方法主要包括:基于模态特征的方法,如模态频率、模态振型、模态阻尼、模态应变能以及模态柔度等与模态信息相关的损伤识别方法;基于频响函数的方法;基于现代信号分析的方法,如HHT、小波包变换等;基于物理参数识别的方法,如Kalman滤波和粒子滤波等;基于模型修正的识别方法;基于智能算法的识别方法,如模糊理论、遗传算法和神经网络
综上,损伤评估和损伤识别既有联系又有区别,前者注重损伤指标的实用性和结构整体评定效果、工程价值及社会影响,后者强调采用多种监测手段和识别算法对结构性能进行精准实时评定,二者的优势互补和协同发展是保障震损结构损伤分析研究不断发展的基石。
近年来,研究者提出了基于新型指标的震损评估方法,例如文[12]提出了基于弹塑性耗能差的损伤评估方法,认为结构的损伤程度可以由实际弹塑性状态(损伤状态)下的累积变形能与理想弹性状态(完好状态)下的累积变形能之差率表示。理论分析和试验验证表明:基于弹塑性耗能差的损伤评估方法机理明确、阈值严格,可以为静动力试验、动力时程分析及易损性分析等提供有效支持,同时能够反映损伤的动态演变过程。然而,与其他基于能量的损伤评估方法(如Park‑Ang模型等)类似,基于弹塑性耗能差的损伤评估方法仍需以地震作用下结构的剪力时程和位移时程作为基本数据。虽然地震中实际结构的位移响应可以由传感器测得,但囿于当前测试技术水平,仍难以直接测得剪力时程。此外,基于弹塑性耗能差的损伤评估方法需预测结构在同一激励且始终处于理想弹性状态下的响应,这也比较难实现。因此,上述方法的局限性要求基于弹塑性耗能差的损伤评估方法的实用性进一步提升。
损伤识别研究中基于物理参数识别的方法为解决上述问题提供了技术支持。以Kalman滤波为代表的物理参数识别方法能够根据部分实测响应信号、状态方程及观测方程对与结构损伤状态相关的未知物理量参数(如刚度、阻尼和动力响应等)进行识
针对目前实际震损结构损伤分析研究存在的问题和不足,本文在弹塑性耗能差方法的基础上提出广义弹塑性耗能差损伤评估模型,规避了实际结构基底剪力和楼层剪力无法直接测量的问题。针对实际结构在理想弹性状态下的地震响应难以获知的局限,提出基于Kalman滤波的理想弹性响应预测方法。在此基础上,提出基于广义弹塑性耗能差的震损结构整体及层间损伤评估方法,可根据结构的加速度及位移响应直接量化损伤程度和演变过程,具有良好的实用性和准确性。
对于如

图1 结构或构件的变形能示意图
Fig. 1 Diagram of deformation energy of structure or member
在单调荷载作用下,基于弹塑性耗能差的损伤计算公式如下:
(1) |
在动力作用下,基于弹塑性耗能差的结构损伤计算公式如下:
(2) |
式中 DE(t)表示t时刻结构的损伤指数;EE(t)和EF(t)分别表示t时刻结构在理想弹性状态下和实际弹塑性状态下的累计变形能;FFi和ΔuFi分别表示i时刻结构在实际弹塑性状态下的基底剪力及顶部位移增量;FEi和ΔuEi分别表示i时刻结构在理想弹性状态下的基底剪力及顶部位移增量;FRi和ΔuRi分别表示i时刻结构在真实弹性状态下的基底剪力和顶部位移增量;β表示将真实弹性状态下的响应放大至与实际弹塑性状态相同的系数。
由于无法直接测得地震下实际结构的剪力值及质量,既有基于弹塑性耗能差的损伤评估方法及目前众多与力或者能量相关的损伤评估方法在实际工程应用中受到限制。此外,该方法还要解决如何推演或预测实际结构在理想弹性状态下响应的问题。因此,为了实现该方法的实用化,亟需提出有效的改进方法。
根据结构动力学理论,单自由度结构在地震下的基底剪力时程为其质量与结构绝对加速度响应的乘积。通过分析
(3) |
式中 m表示结构或构件的质量,表示i时刻结构或构件在弹塑性状态下的加速度响应,表示i时刻结构或构件在理想弹性状态下的加速度响应。
由
损伤程度 | 损伤现象描述 | 损伤指数 |
---|---|---|
轻微损伤 | 少数承重构件有微小残余变形,个别非结构构件破坏明显,非结构构件发生不同程度破坏 | 0~0.25 |
中等破坏 | 多数承重构件有轻微残余变形,部分承重构件有明显残余变形,个别非结构构件严重破坏 | 0.25~0.55 |
严重破坏 | 多数承重构件严重破坏或部分倒塌 | 0.55~0.80 |
倒塌 | 多数承重构件倒塌 | 0.80~1.00 |
对于实际震损结构,仅能测得地震下的弹塑性响应,不能直接获取结构在同一地震作用下且始终处于理想弹性状态时的响应,为此还需提出切实可行的预测方法。由于主震发生前通常会发生振幅较小的前震,而此时结构一般处于弹性状态,因此可以利用相关前震信号和结构响应监测信号进行等效多自由度结构的物理参数识别,进而建立相关动力学模型。在获取主震下的激励和结构响应信号后,用该弹性模型快速推演出针对主震的结构理想弹塑性响应,从而解决上述问题。具体的物理参数识别方法可以基于Kalman滤波方法开展,下文将进行深入研究。
为了克服传统物理参数识别方法在不完备信息条件下多自由度结构识别能力不足的局限,基于Kalman滤波原理,本文认为可根据结构在前震下的响应对结构等效刚度及等效质量进行识别,建立结构等效模型,进而实现主震下的结构响应推演。与传统思路不同,本方法不严格要求等效模型与实际结构的质量、刚度相一致,只要求确保在响应层面能够与实际结构在前震下的响应相吻合,从而能够利用等效模型计算实际结构在理想弹性状态下的响应。对于实际的多自由度结构,可用串联层间模型来表征,其动力方程为:
(4) |
式中 u,和分别表示结构的位移、速度和加速度;M,C和K分别表示结构的质量矩阵、阻尼矩阵和刚度矩阵;F表示外荷载向量;I表示单位对角矩阵;表示动力荷载加速度。
阻尼矩阵C可用Rayleigh阻尼表示,其形式为:
(5) |
式中 a和b分别表示比例系数,可由下式确定:
(6) |
式中 ωi和ωj分别表示结构第i和j阶频率,可通过结构的频率方程求得;ξ为结构阻尼比,其值可根据中国《建筑结构抗震设计规范
根据结构动力方程以及实测的结构位移响应u和速度响应,结合Kalman滤波方法,对结构的质量和刚度进行识别。此时状态空间方程可以描述为:
(7) |
以三自由度串联结构体系为例,将状态方程中的质量和刚度看作未知参数,此时Kalman识别中的全部状态量包括:各自由度相对基础的水平位移u1,0,u2,0和u3,0;各自由度相对基础的水平速度u4,0,u5,0和u6,0;各层刚度u7,u8和u9;各层质量u10,u11和u12。此时三自由度串联结构体系的状态方程可以表示为:
(8) |
式中 ci,j表示将卡尔曼滤波识别的质量和刚度代入Rayleigh阻尼中计算得到的阻尼矩阵C中第i行、第j列的数值。
利用Kalman滤波算法和体系状态方程,并根据地震前震信号和结构各层响应信号可以识别出近似的结构各层等效质量和等效刚度,建立实际框架结构的等效弹性模型,进而实现推演主震时结构理想弹性状态下的响应。整体计算流程如

图2 等效模型建立流程
Fig.2 Establishing process of equivalent model
为实现对震损结构损伤程度的精确量化,结合广义弹塑性耗能差以及等效弹性模型的建立方法,本文提出基于Kalman滤波及广义弹塑性耗能差的震损结构整体及层间损伤的评价方法及流程。其核心思想是根据前震下真实结构各层的地震响应,借助Kalman滤波建立一个始终处于理想弹性状态的结构等效模型,由该模型计算震损结构在理想弹性状态下的响应,通过比较震损结构在理想弹性状态和实际弹塑性状态下累计变形能的差异性评估震损结构整体和层间的损伤程度。
结构的顶层位移是各层相对位移量累加之和,可通过地震中结构顶层位移或位移角的变化反映结构整体的抗震性能。本文建议将建筑结构的顶层相对位移作为整体性能指标,通过对比结构在实际弹塑性状态以及理想弹性状态下的累计变形能差值评估结构的整体损伤程度,相应的计算公式如下:
(9) |
式中 DT,E(t)表示震损结构在第t时刻的整体损伤程度;uT,Fi和∆uT,Fi分别表示实测所得震损结构在第i时刻的顶层位移和顶层位移增量;uT,Ei和∆uT,Ei分别表示根据等效模型得到的震损结构在理想弹性状态下第i时刻的顶层位移和顶层位移增量。
按类似思路,结构层间损伤可以通过结构各层层间响应的变化来表征,相应的层间损伤计算式为:
(10) |
式中 Dn,E(t)表示震损结构的第n层在第t时刻的层间损伤程度;un,Fi和un,Ei分别表示结构在实测和理想弹性状态下第n层第i时刻的位移;Δun,Fi和Δun,Ei分别表示第i时刻结构第n层在实测和理想弹性状态下的相对位移增量;Δ(un,Fi-un-1,Fi)表示i时刻结构第n层和第n-1层在实测状态下的相对位移增量的差值;Δ(un,Ei-un-1,Ei)表示i时刻结构第n层和第n-1层在理想弹性状态下的相对位移增量的差值。
以上方法的计算流程如

图3 震损结构整体及层间损伤的评价流程
Fig.3 Evaluation process of overall and interlayer damage of seismic-damaged structure
考虑实际双向水平地震作用下,结构在两方向上均可产生损伤。若需要进行结构多维损伤评估,可借助动力传感器获取结构在两个主轴方向上的地震动和结构响应,并分别采用上述方法计算震损结构两个主轴方向上的损伤程度,再根据下式综合量化评价结构的整体损伤程度:
(11) |
式中 DT,Ex和DT,Ey分别表示震损结构在x和y两个主轴方向上的整体损伤程度。
整体损伤和层间损伤作为结构损伤在不同层面的量化方式具有密切的内在关联,整体损伤是层间损伤的宏观综合表征,层间损伤是整体损伤的细节表现。考虑各层间损伤对综合损伤的贡献,本文提出以层间累计应变能作为楼层重要性系数的结构整体损伤程度加权计算方法,公式如下:
(12) |
式中 m表示建筑结构的总楼层数;En表示震损结构第n层的累计弹塑性应变能;Dn,E表示震损结构第n层的层间损伤指数。
为验证本文所提基于Kalman滤波及广义弹塑性耗能差的震损结构损伤评估方法的准确性及实用性,以7度区某五层和某十层钢筋混凝土框架结构为例进行分析。两个建筑结构的平面布置及立体构造分别如图

图4 结构平面布置图(单位:mm)
Fig. 4 Layout of structural plan(Unit: mm)

图5 结构三维示意图
Fig.5 Three-dimensional diagram of structure
首先,采用九寨沟白河地震台站监测的汶川地震波并将其加速度峰值调至0.15g后输入至五层建筑结构模型中进行时程分析。如

图6 汶川地震波记录
Fig.6 Seismic wave record of Wenchuan
随后,提取五层建筑结构在前20 s内的位移响应、加速度响应以及地震动信号,结合前述方法对结构各层的等效质量和等效刚度进行识别,部分楼层等效质量和等效刚度的识别结果如

图7 卡尔曼滤波识别结果
Fig.7 The recognition result of Kalman filter
根据得到的等效刚度和等效质量建立五层结构的等效层模型,并向等效模型输入汶川地震波进行时程分析,原结构和等效结构典型楼层的位移、速度、加速度响应对比如

图8 汶川地震波作用下结构第三、五层响应对比
Fig.8 Responses comparison of the structure’s third and fifth floors under Wenchuan wave

图9 Loma波作用下结构三、五层响应对比
Fig.9 Responses comparison of the structure’s third and fifth floors under Loma wave
采用相同的方法,建立十层结构的等效模型。之后,在五层结构和十层结构的有限元模型中分别输入加速度峰值为0.15g和0.30g的汶川波、Loma波、Taft波、El‑Centro波等7条不同的地震波,计算结构在不同地震作用下的响应,并将相同的地震波输入至五层及十层结构的等效模型中获取结构在理想弹性状态下的响应。根据建筑结构在理想弹性状态下的响应(利用等效模型获取)和建筑结构在实际弹塑性状态下的响应(由有限元模型直接获取),并通过本文所提基于Kalman滤波及广义弹塑性耗能差的震损结构损伤计算方法计算结构在不同加速度峰值以及不同地震波下的整体及层间损伤程度,同时与传统弹塑性耗能差方法的结果进行对比。典型地震动下结构整体损伤的对比如图

图10 五层结构整体损伤指数 (PGA=0.30g)
Fig.10 Overall damage index of 5-story structure (PGA=0.30g)

图11 十层结构整体损伤指数(PGA=0.30g)
Fig.11 Overall damage index of 10-story structure (PGA=0.30g)

图12 五层结构层间损伤指数 (PGA=0.30g)
Fig.12 Interlayer damage index of 5-story structure (PGA=0.30g)

图13 十层结构层间损伤指数 (PGA=0.30g)
Fig.13 Interlayer damage index of 10-story structure (PGA=0.30g)
从图
由基于瞬时频率变化率以及基于小波包分解的频率变化率的结构整体损伤指数可以看出,这两种方法得到的结构最终损伤程度值与弹塑性耗能差方法以及广义弹塑性耗能差方法的结果大致相同。但由于损伤评价机理不同,使得损伤的演变规律存在一定的差异,其中基于瞬时频率变化率的损伤演变过程与弹塑性耗能差方法更为接近。但囿于现有信号分析技术,实际结构的瞬时频率尚难以准确提取,且环境噪声对信号分析精度也存在影响,因此,与频率相关的结构损伤评估方法实用性相对较差。
图
在加速度峰值分别为0.15g及0.30g的7条不同地震波下的结构整体及层间损伤统计平均值结果,如图

图14 结构整体损伤指数统计值
Fig.14 Statistical value of structural overall damage index


图15 结构层间损伤指数统计值
Fig.15 Statistical value of structural interlayer damage index
可以发现:随着地震动幅值增加,结构的整体损伤程度及层间损伤程度不断加剧;受地震波的频谱特性和随机性影响,同一结构在不同地震波下的损伤程度也具有一定的随机性。对于五层结构,由广义弹塑性耗能差方法得到的整体损伤的统计平均结果与传统方法结果十分接近,其最大偏差不超过±20%,精度良好。对于十层结构,在某些地震下两种方法的整体损伤值偏差较大。这主要是因为随着楼层数的增加,未知物理量增多,Kalman滤波的识别精度有所下降。由两种方法获得的楼层损伤统计平均值比较接近,且各楼层的损伤程度排序和薄弱层的判断依然是基本准确的。由于震后应急评估和抗震救援着重于从宏观层面量化结构整体损伤程度,以判断震损结构的可修性,而对层间损伤程度的关注度相对较小,因此广义弹塑性耗能差方法的精度在统计意义上是满足工程需求的。
五层结构在PGA=0.3g的三向汶川波作用下,由
针对目前损伤评估及损伤识别方法在实际工程运用中存在的诸多不足,本文在传统基于弹塑性耗能差的损伤评估方法基础上,用剪力比的形式将结构质量效应消除,提出广义弹塑性耗能差损伤评估方法。利用前震信号识别且利用主震信号进行推演的思路,提出了基于Kalman滤波及广义弹塑性耗能的震损结构整体及层间损伤的评价方法。其主要研究结论如下:
(1) 提出广义弹塑性耗能差损伤评估方法,根据结构或构件在正常运营及损伤状态下的位移、加速度监测信号能够快速量化其损伤程度。
(2) 采用Rayleigh阻尼,结合Kalman滤波识别结构特征物理参数的原理,提出根据前震下的结构响应识别结构等效质量和等效刚度并建立等效弹性模型的方法和流程,力求与实际结构在前震下的响应保持一致,由此计算实际结构在理想弹性状态下的响应。
(3) 针对震损结构提出基于Kalman滤波及广义弹塑性耗能差的震损结构整体及层间损伤评价方法,该方法仅以结构的真实弹塑性响应及由等效弹性模型获取的理想弹性响应作为评判依据,摆脱了目前诸多损伤评价或识别方法需已知结构质量或剪力时程变化的桎梏,具有较高的实用性。
(4) 通过对五层和十层钢筋混凝土框架结构进行分析证明本文所提方法准确性较高,能够反映结构整体损伤的动态演变过程以及结构薄弱层的分布情况,宜在实际工程中推广应用。
综上所述,利用本文所提基于Kalman滤波及广义弹塑性耗能差的损伤评估方法,能够根据结构在前震和主震下的位移及加速度响应及时、准确、合理地量化震损结构的损伤程度,规避了现有损伤评估或识别方法需已知结构质量或剪力等难以实测物理量等不足。该方法能够为实际震损结构损伤评估及修复加固提供有效指导或建议。
参考文献
李宏男, 肖诗云, 霍林生. 汶川地震震害调查与启示 [J]. 建筑结构学报, 2008, 29(4):10-19. [百度学术]
Li Hongnan, Xiao Shiyun, Huo Linsheng. Damage investigation and analysis of engineering structures in the Wenchuan earthquake[J]. Journal of Building Structures, 2008, 29(4):10-19. [百度学术]
张超, 翁大根, 彭林海. 震损钢筋混凝土结构减震加固研究进展 [J]. 工程抗震与加固改造, 2012, 34(3):112-120+107. [百度学术]
Zhang Chao, Weng Dagen, Peng Linhai. Research review on seismic retrofit of earthquake-damage RC structure by using energy dissipation technology [J]. Earthquake Resistant Engineering and Retrofitting, 2012, 34(3):112-120+107. [百度学术]
李宏男, 何浩祥. 利用能力谱法对结构地震损伤评估简化方法 [J]. 大连理工大学学报, 2004, 43(2):267-270. [百度学术]
Li Hongnan, He Haoxiang. A simplified approach for seismic damage assessment based on capacity spectrum method [J]. Journal of Dalian University of Technology, 2004, 43(2):267-270. [百度学术]
王丰, 李宏男. 基于简化IDA的结构地震损伤评估方法 [J]. 工程力学, 2018, 35(12):194-202. [百度学术]
Wang Feng, Li Hongnan. Seismic damage assessment procedure of structures by using simplified IDA [J]. Journal of Engineering Mechanics. 2018, 35(12):194-202. [百度学术]
宗周红, 牛杰, 王浩. 基于模型确认的结构概率损伤识别方法研究进展 [J]. 土木工程学报, 2012, 45(8):121-130. [百度学术]
Zong Zhouhong, Niu Jie, Wang Hao. A review of structural damage identification methods based on the finite element model validation [J]. Journal of Civil Engineering, 2012, 45(8):121-130. [百度学术]
Kim Jeong-Tae, Ryu Yeon-Sun, Cho Hyun-Man, et al. Damage identification in beam-type structures: frequency-based method vs mode-shape-based method [J]. Engineering Structures, 2003, 25(1) :57-67. [百度学术]
Park Y, Ang A H. Mechanistic seismic damage model for reinforced concrete [J]. Journal of Structural Engineering, 1985, 111(4):722-739. [百度学术]
Kratzig W B, Meyer I F, Meskouris K. Damage evolution in reinforced concrete members under cyclic loading [C]. Proceedings of the 5th International Conference on Structural Safety and Reliability, San Francisco, 1989:795-802. [百度学术]
杜修力, 欧进萍. 建筑结构地震破坏评估模型 [J]. 世界地震工程, 1991, 7(3):52-58. [百度学术]
Du Xiuli, Ou Jinping. Earthquake damage evaluation model of building structure [J]. World Earthquake Engineering, 1991, 7(3):52-58. [百度学术]
王东升, 冯启民, 王国新. 考虑低周疲劳寿命的改进Park-Ang地震损伤模型[J]. 土木工程学报, 2004, 37(11):41-49. [百度学术]
Wang Dongsheng, Feng Qimin, Wang Guoxin. A modified Park-Ang seismic damage model considering low-cycle fatigue life [J]. Journal of Civil Engineering, 2004, 37(11):41-49. [百度学术]
安永辉. 基于振动信息的结构损伤识别的几种新方法 [D]. 大连:大连理工大学, 2013. [百度学术]
An Yonghui. Several new methods for vibration-based damage identification of structures [D]. Dalian: Dalian University of Technology, 2013. [百度学术]
何浩祥, 陈奎, 范少勇. 基于弹塑性耗能差率的地震损伤评估模型及分析方法 [J]. 振动工程学报, 2018, 31(3):382-390. [百度学术]
He Haoxiang, Chen Kui, Fan Shaoyong. Seismic damage model based on differential ratio of elastic plastic dissipated energy and application [J]. Joural of Vibration Engineering, 2018, 31(3):382-390. [百度学术]
Wu Meiliang, Smyth Andrew W. Application of the unscented Kalman filter for real-time nonlinear structural system identification[J]. Journal of Structural Control and Health Monitoring, 2007, 14(7): 971-990. [百度学术]
何浩祥, 吕永伟, 韩恩圳. 基于静动力凝聚及扩展卡尔曼滤波的连续梁桥损伤识别 [J]. 工程力学, 2015, 32(7):156-163. [百度学术]
He Haoxiang, Lü Yongwei, Han Enzhen. Damage detection for continuous girder bridge based on static-dynamic condensation and EKF [J]. Journal of Engineering Mechanics, 2015, 32(7):156-163. [百度学术]
Yang J N, Lin S, Huang H, et al. An adaptive extended Kalman filter for structural damage identification[J]. Structural Control and Health Monitoring, 2006, 13(4): 849-867. [百度学术]
吴新亚, 周丽. 基于有限范围自适应卡尔曼滤波的结构损伤识别方法 [J]. 振动工程学报, 2007, 20(4): 401-406. [百度学术]
Wu Xinya, Zhou Li. A finite-horizon adaptive Kalman filter for structural damage identification [J]. Journal of Vibration Engineering, 2007, 20(4): 401-406. [百度学术]
Lei Y, Chen F, Zhou H. An algorithm based on two-step Kalman filter for intelligent structural damage detection [J]. Structural Control and Health Monitoring, 2015, 22 (4): 694-706. [百度学术]
李杰,陈隽. 未知输入条件下的结构物理参数识别研究 [J]. 计算力学学报, 1999, 16 (1):32-40. [百度学术]
Li Jie, Chen Jun. Identification of structural parameters with unknown input [J]. Journal of Computational Mechanics, 1999, 16(1): 32-40. [百度学术]
Lei Y, Liu C, Jiang Y Q, et al. Substructure based structural damage detection with limited input and output measurements [J]. Smart Structures & Systems, 2013, 12(6):619-640. [百度学术]
Carrella A, Ewwins D J. Identifying and quantifying structural nonlinearities in engineering applications from measured frequency response functions [J]. Mechanical Systems & Signal Processing,2011, 25(3):1011-1027. [百度学术]
何浩祥, 胡一凡, 吴山. 基于耗能差指数和SIR模型的结构易损性及可恢复能力评估 [J]. 土木工程学报, 2020, 53(4):11-22. [百度学术]
He Haoxiang, Hu Yifan, Wu Shan. Evaluation of structural vulnerability and resilience based on energy dissipation difference index and SIR model [J]. Journal of Civil Engineering, 2020, 53(4):11-22. [百度学术]