Influence of Shear Walls on Optimizing Seismic Internal Force Distribution of RC Structures
-
摘要: 为研究地震作用下钢筋混凝土框架结构倒塌机理,验证基于内力分配得到的抗倒塌措施的合理性,设计2组缩尺比例为1∶4的框架结构模型并开展振动台对比试验,单次双向输入地震动峰值加速度为1.0 g,分析地震作用下填充墙-框架结构变形模式和破坏特点,对比分析有、无落地剪力墙框架结构柱间内力分配规律。研究结果表明,填充墙影响框架结构整体变形及柱破坏模式,横向满砌填充墙约束结构扭转变形,即使在双向地震作用下,偏心结构也未发生扭转;由于窗下半高连续填充墙的约束作用,窗间柱抗侧刚度变大,在地震作用下承担的地震剪力是不受半高连续填充墙约束柱的6~8倍;设置落地剪力墙可优化底层柱间地震剪力分配,使框架柱地震剪力分配趋于均匀,避免结构因“凝震聚力”而发生倒塌,实现“大震不倒”。Abstract: In order to investigate the collapse mechanism of the RC frames and verify the rationality of the anti-collapse measures based on the internal force distribution under the earthquake, two models with the scale ratio of 1:4 were designed and the shaking table tests were conducted. In this paper, the structure deformation and column failure of the RC frame structures with infilled walls suffered from bidirectional earthquake reappeared in detail, and the internal force distribution of each column in the structures with or without shear walls was compared. Results show that the infilled walls have significant effects on the columns. The transverse filled walls can restrain the torsion of the structure. Thus, even under the bidirectional earthquake, the eccentric structures do not twist. Due to the restraints of the half-height continuous infilled walls under the window, the seismic shear force of the window column is 6-8 times that of unconstrained ones. The comparison of internal force distribution between columns based on strain of the column shows that the shear walls can optimize the distribution of seismic shear forces among the columns, make the internal force distribution of the columns on each axis balanced, and avoid the structure collapse.
-
Key words:
- RC frame /
- Earthquake collapse /
- Shear walls /
- Internal force distribution /
- Shaking table test
-
引言
2017年8月8日21时19分,四川省北部阿坝州九寨沟县发生MS 7.0地震。根据中国地震台网中心发布的信息,本次地震的震中位于33.20°N、103.82°E,震源深度10km。震中在九寨沟县城以西约39km处,东偏北方向距离陇南市105km,南距成都市285km。本次地震的最大烈度为Ⅸ度,共造成25人死亡,525人受伤,6人失联,受灾人数超过17万(中国地震局,2017)。地震导致7万多间房屋不同程度受损,其中76间房屋倒塌。此次地震获取到强震观测记录的台站数量有限,其中距离震中最近的九寨百河台站(台站编号51JZB),南北向地震记录的峰值加速度(PGA)达到了0.19g,是本次地震中观测到的最大加速度(图 1)。
大量的震害资料表明,地震动是造成结构破坏、诱发地质灾害(例如场地液化、地表裂缝和滑坡)的主要驱动力。地震动记录对于研究地震危险性、地震引发的地质灾害和结构抗震设计具有重要的意义。在缺少强震观测记录的情况下,可以通过有限的实际观测记录,借鉴科学合理的方法,综合考虑震源、传播路径和场地条件人工合成地震动,其结果可作为结构动力反应分析和灾后抗震设计的依据。
本文采用Wang等(2015)改进的随机有限断层方法合成2017年九寨沟地震的近场地震动。首先选取地震的有限断层模型和随机方法的输入参数,并通过与实际记录对比估计地震的应力降;采用随机有限断层法合成震源附近台站的加速度时程;并将得到的模拟地震动与实际记录和烈度图进行对比,验证结果的有效性。
1. 随机有限断层法
根据有限断层法的思想,断层破裂可以离散为N个子断层,每个子断层可以近似认为是1个点源。破裂过程从震源开始,逐渐向断层边缘扩散,当传播到每个子断层时,子断层发生破裂并产生地震波。将子断层产生的地震波按照到达的先后次序进行叠加,即可得到整个断层在某个观测点产生的地震动。Beresnev等(1997)将有限断层模型和随机点源方法(Boore,1983)结合,提出了随机有限断层法,可以用来模拟断层附近的高频地震动。Motazedian等(2005)进一步发展了这种方法,提出了基于动态拐点频率的随机有限断层法(EXSIM方法),解决了总能量受到子断层尺寸影响等问题。Wang等(2015)改进了Motazedian等(2005)的方法,采用一种双拐点震源谱,能够反映破裂过程中累积地震矩对震源谱的影响。
本文采用改进的随机有限断层法(Wang等,2015)模拟九寨沟地震动分布,综合考虑震源、传播路径、场地的影响,有限断层模型中第j行、第i列的子断层在场点产生地震动的加速度傅里叶谱可以表示为:
$$ {Y_{ij}}(f){\rm{ = }}E({M_{0ij}}, {f_{0ij}}, f)G({R_{ij}}){H_{ij}}{{\rm{e}}^{ - \frac{{{\rm{ \mathsf{ π} }}f{R_{ij}}}}{{Q(f)\beta }}}}S(f){{\rm{e}}^{ - {\rm{ \mathsf{ π} }}f\kappa }} $$ (1) 其中,f为频率(Hz);Rij为第j行、第i列子断层到观测点的距离(km);$E({M_{0ij}}, {f_{0ij}}, f) $代表震源谱(cm/s);$G({R_{ij}}) $为几何衰减函数;$ {H_{ij}}$为低频比例因子;$ Q(f)$为地壳品质因子;$ \beta $为剪切波速(km/s);$ S(f)$为与频率相关的场地放大因子;$ \kappa $为高频衰减因子。第j行、第i列子断层的震源谱按照式(2)计算:
$$ E({M_{0ij}}, {f_{0ij}}, f){\rm{ = }}\frac{{C{M_{0ij}}{{(2{\rm{ \mathsf{ π} }}f)}^2}}}{{{{[1 + {{(f/{f_{0ij}})}^a}]}^b}}} $$ (2) 其中,C为常数,按照Motazedian等(2005)给出的公式计算;${M_{0ij}} $是第j行、第i列子断层的地震矩(dyne-cm);a和b是控制震源谱形状和幅值的动态参数,可由断层已破裂部分的地震矩确定:
$$ \begin{array}{l} a{\rm{ = 6}}{\rm{.592}} - {\rm{0}}{\rm{.22}}\lg {M_0}(t)\\ b{\rm{ = 2}}/a \end{array} $$ (3) 其中,$ {M_0}(t)$表示在时间t断层上已破裂区域的地震矩,可根据地震的总地震矩,以滑动分布为权重分配求得。式(2)中拐点频率f0ij按照式(4)可表示为:
$$ {f_{0ij}} = 4.9 \times {10^6}\beta {\left({\mathit{\Delta }\sigma \mathit{/}{\mathit{M}_{0ij}}} \right)^{1/3}} $$ (4) 其中,$ \mathit{\Delta }\sigma $代表地震的应力降(bar)。
除震源谱模型外,式(1)中其他几项函数代表传播路径和场地条件对地震动傅里叶谱的影响,如几何衰减、场地放大因子和高频衰减因子。这些函数在相关文献(Boore,1983;Motazedian等,2005)中有详细的介绍,本文不再赘述。另外,Wang等(2015)采用上述改进的随机有限断层法合成了美国1994年Northridge地震的近场地震动,结果和实际记录较为符合,验证了方法的有效性。
2. 九寨沟强震的地震动模拟
2.1 震源模型及输入参数
基于远场P波和SH波的波形数据,王卫民等(2017)反演了九寨沟MS 7.0地震的震源破裂过程,得到了震源机制以及发震断层滑动分布的初步结果。结果显示地震矩为6.7×1025dyne-cm,对应的矩震级为MW 6.5;断层的走向、倾角和滑动角分别为148.5°、68.9°和-3.1°,属于左旋走滑型地震。基于此结果,本文得出的断层滑动分布见图 2,从图中可以看出,断层上的破裂滑动大部分集中在震源附近区域,即深度6—12km范围内,最大滑动量为85cm。此外,在震源东南约16km的位置,也有幅值小于40cm的破裂滑动。
表 1列出了随机有限断层方法的输入参数,其中大部分是其他学者根据四川及周边地区的大量地震记录分析得到的,反映了该地区地震动衰减特性。其中,几何衰减采用喻烟(2012)提出的汶川地震区地震动估计的经验模型,高频衰减因子采用Sun等(2013)分析汶川余震记录得到的结果,与距离相关的持时采用Atkinson等(1995)提出的分段模型。基于对九寨沟地震记录S波的傅氏谱反演结果,王宏伟等(2017)确定了这次地震震中区域的地壳品质因子为$ Q(f) = 84.9{f^{0.71}}$,体现了较强的地震波衰减特性,本文的地震动模拟中采用了此结果。对于场地效应,地表到30m深度的平均剪切波速(VS30)通常被用作场地类别划分的依据,例如美国地震减灾项目(NEHRP)。喻烟(2012)根据我国四川与甘肃地区的钻孔数据,假定钻孔底部到30m深度的波速与孔底一致,计算了台站场地的VS30。本文根据这些结果,确定了台站对应的NEHRP场地类别,并采用Boore等(1997)针对不同类别场地提出的场地放大因子。对于所有台站,断层破裂的脉冲面积百分比采用50%。本文通过对比模拟和观测的地震动,按照误差最小原则估算本次地震的应力降参数。
表 1 模型输入参数Table 1. Input parameters of model参数 值或表达式 断层走向和倾角 148.5°,68.9° 断层尺寸 54km×24km 子断层尺寸 2km×2km 剪切波速(β) 3.5km/s 介质密度(ρ) 2.8g/cm3 破裂速度(VR) 0.8×β 应力降 4.0MPa 品质因子 $ Q(f) = 84.9{f^{0.71}}$ 高频衰减因子(κ) 0.015 路径持时T=dR,参数d 0.0(R≤10km)
-0.16(10km<R≤70km)
-0.03(70km<R≤130km)
0.04(R>130km)几何衰减G(R)=Rb,参数b -1.0(R≤150km)
-0.5(R>150km)场地放大因子 参见Boore等(1997) 脉冲面积百分比 50% 2.2 应力参数估计
首先,选取震源距小于150km、记录的PGA大于10cm/s2的台站,将这些台站取得的记录作为估计应力降的依据。表 2给出了选取的8个台站信息,包括经纬度、震源距、VS30和相应的NEHRP场地类别。估计应力参数主要采用“试错法”。考虑1.0—6.0MPa的应力降,增量为0.5MPa,对于每个应力降,采用2.1节中的震源模型和输入参数合成8个台站的地震动加速度时程。考虑随机方法Gauss白噪声幅值的波动性,每个台站模拟30次,计算平均的反应谱。定义总体平均误差如下:
表 2 用于估计应力参数的台站信息Table 2. The information of stations selected for the estimation of stress drop台站名称 台站代码 纬度/°N 经度/°E 震源距/km VS30/m·s-1 NEHRP场地类别 岷县 62MXT 34.4 104.0 135.1 301 D 迭部 62DIB 34.1 103.2 117.8 C 沙湾 62SHW 33.7 104.5 82.5 384 C 九寨百河 51JZB 33.3 104.1 27.2 328 D 九寨勿角 51JZW 33.0 104.2 39.7 428 C 九寨永丰 51JZY 33.2 104.3 42.3 321 D 茂县叠溪 51MXD 32.0 103.7 133.7 268 D 平武木座 51PWM 32.6 104.5 89.8 373 C 注:迭部台由于缺少钻孔资料和VS30,场地类型假定为NEHRP-C。 $$ r(f) = \frac{1}{N}\sum\limits_{i = 1}^N {\lg \frac{{S{A_{{\rm{obs}}, i}}(f)}}{{S{A_{{\rm{sim}}, i}}(f)}}} $$ (5) 其中,$ S{A_{{\rm{obs}}, i}}(f)$,$ S{A_{{\rm{sim}}, i}}(f)$分别表示第i个台站记录和模拟的加速度反应谱(cm/s2),N表示台站数量。对于每个应力降,根据合成结果计算出总体平均误差。需要说明的是,由于随机方法通常用于模拟高频地震动,不一定能够体现出低频地震动的特征,因此本文考虑的频率范围是1—30Hz。根据上述准则,得到误差最小的应力降为4.0MPa。王宏伟等(2017)通过反演本次地震记录S波的傅里叶谱,确定了地震的应力降为3.854MPa,与本文通过“试错法”得到的结果较为接近。原因之一是本文的模拟采用了同一篇文献(王宏伟等,2017)中基于反演地震记录确定的品质因子$Q(f) $。图 3(a)给出了4.0MPa应力降误差随频率的变化以及误差±1个标准差的范围。从图中可以看出,在1—30Hz频率范围内,误差基本较小,绝对值小于0.2。作为对比,我们也计算了EXSIM方法得到的误差(图 3(b)),其中所有输入参数(包括应力降)和本文采用的方法相同。通过对比可以明显看出,EXSIM方法的误差在频率大于3Hz时明显大于本文的结果;同时在1Hz附近,EXSIM方法的平均误差为负,表示模拟地震动在该频段高估了观测地震动的反应谱幅值。由此可见,相比EXSIM方法,采用Wang等(2015)的方法得到的结果与本次地震记录更相符。
2.3 模拟地震动的时程和反应谱
按照上述随机有限断层方法,采用表 1的模型输入参数,分别模拟了6个强震台站的加速度时程。在本文选取的台站中,虽然大部分都有剪切波速的钻孔数据,但多个场地的钻孔未达到基岩,如果进行土层反应分析需要进行一定的假设或者延拓剪切波速剖面,存在较大的不确定性。因此,在模拟中根据台站的NEHRP场地类型,采用Boore等(1997)提出的NEHRP-C、NEHRP-D场地的放大系数,近似考虑场地效应的影响。图 4给出了6个台站模拟和观测地震动的加速度时程对比。为了便于进一步定量比较,表 3列出了模拟和观测的地震动PGA值。需要说明的是,随机有限断层法得到的地震动对应于任意水平分量。
表 3 模拟和观测地震动的峰值加速度Table 3. The PGA from synthetic and recorded ground motions台站代码 PGA(记录)/cm·s-2 PGA(模拟结果)/cm·s-2 EW分量 NS分量 51JZB 129.50 185.02 174.78 51JZW 73.84 91.70 134.90 51JZY 45.77 66.71 85.10 51MXD 12.21 23.32 11.13 62SHW 18.57 20.48 19.27 51PWM 18.60 20.87 21.26 通过对比模拟和记录的地震动可以发现,在峰值加速度最大的台站51JZB,模拟地震动的PGA值介于东西向和南北向记录的PGA值之间,与实际记录较为符合;在震源距较大的2个台站62SHW和51PWM,拟地震动的PGA与实际记录也很吻合;另外,51JZW和51JZY模拟地震动的PGA稍大于实际记录,51MXD模拟地震动的PGA小于地震记录,原因可能有以下3点:①采用Boore等(1997)的场地放大系数不能充分考虑场地条件对地震动的影响,如非线性土层反应和不规则地形影响;②本文采用的震源滑动分布模型是通过反演远场波形数据得到的,没有考虑近断层地震动和同震位移场,不一定能够完全反映出震源特征;③本文模拟采用的模型参数大多是基于分析汶川地震数据得到的,这些参数也会存在不同地震、不同地区之间的差异。
从图 4可以看出,大部分台站模拟地震动的持时在不同程度上小于实际记录。需要指出,随机有限断层法模拟的实际上是地震动的S波。一方面,在近场地震记录中P波与S波2种成分相互叠加,会导致记录的持时大于模拟地震动;另一方面,部分台站(如51JZY、51MXD和51PWM)的地震记录中面波成分很丰富,无法通过随机方法模拟出来;除此之外,近地表土层反应和不规则地形也可能造成地震记录的持时较长。解决这些问题需要更复杂的地震动模拟方法。
图 5给出了6个台站模拟和观测地震动的加速度反应谱。通过对比可以看出,在本文考虑的周期(T)范围内(0.03—1.0s),大部分台站模拟结果的反应谱与实际记录较为符合。51JZW台站模拟地震动的反应谱在短周期(T<0.1s)大于实际记录的反应谱,51MXD台站模拟地震动低估了T>0.1s地震记录的反应谱。这些差异说明本文采用的随机有限断层方法和模型输入参数仍有可改进之处。综合考虑随机性方法与确定性方法的混合方法,能够弥补随机性方法在模拟低频地震动和面波等方面的不足,但是需要提供更多的信息,如地壳介质波速模型等。
2.4 模拟的PGA分布图
为了进一步研究模拟地震动加速度峰值的空间分布,我们在发震断层附近以5km为间隔选取了1548个网格点。对于每个点,采用2.1节介绍的断层模型和输入参数模拟地震动,其中场地类型假定为NEHRP-C,并采用Boore等(1997)对此类场地提出的放大系数。考虑随机模拟结果的不确定性,每个点模拟20次,最后取平均的PGA作为结果。根据所有网格点的模拟PGA,画出该区域的PGA分布图(图 6(a))。从图中可以看出,模拟PGA在近断层区域的分布明显受到断层位置以及滑动分布的影响;模拟PGA最大的区域位于震中附近,最大值大约为850cm/s2。为了验证结果的有效性,我们对比了模拟PGA的分布图与中国地震局(2017)发布的本次地震烈度图(图 6(b)),可以看出,模拟PGA最大的区域与本次地震最大烈度区(Ⅸ度)的位置基本一致。除此之外,PGA和烈度的等值线都呈椭圆形,沿断层走向方向较长,垂直断层走向方向较短,说明本文模拟的PGA空间分布基本符合观测到的地震烈度分布情况。
3. 结论
本文采用改进的随机有限断层法合成了2017年九寨沟MS 7.0地震的近场地震动。结果表明:
(1)选取合适的震源模型和输入参数,通过“试错法”,按照误差最小原则估计本次地震的应力降为4.0MPa,与王宏伟等(2017)反演地震记录S波傅氏谱的得到结果3.854MPa基本一致。与EXSIM方法相比,Wang等(2015)的方法得到的结果在频域上与实际地震记录更为符合。
(2)按照美国NEHRP场地划分标准,依据VS30确定了断层附近台站的场地类型。在模拟地震动过程中,通过场地放大因子近似考虑场地效应的影响。最后,得到了强震台站的模拟加速度时程和反应谱,以及研究区域的模拟PGA分布图。部分台站的模拟地震动和实际记录在峰值加速度和反应谱方面较为符合,模拟的PGA分布图反映出地震烈度图的基本特征,验证了本文结果的有效性。
(3)部分台站的模拟结果在幅值、持时以及反应谱方面与实际记录存在不同程度的差异,说明本文采用的随机有限断层法和模型参数仍有可改进之处。
(4)对九寨沟强震的模拟结果说明,在缺少近断层记录的情况下,综合考虑震源、传播路径和场地条件人工合成地震动,总体上能够有效地重现实际地震记录的峰值和反应谱特征,为该地区的地震危险性分析和抗震设计提供一定依据。
致谢: 感谢国家强震动台网中心为本研究提供九寨沟地震的强震记录。 -
表 1 试验模型主要相似关系
Table 1. The similarity ratio of the experimental model
项目 相似比 模型Ⅰ
(人工质量13.4 t)模型Ⅱ
(人工质量11.6 t)长度lr 0.25 0.25 弹性模量Er 0.55 0.55 材料密度ρr 2.20 2.20 应力σr 0.55 0.55 时间tr 0.50 0.50 速度vr 0.50 0.50 加速度ar 1.00 1.00 表 2 不同时刻模型Ⅰ底层柱端应变
Table 2. The strain of the columns of model Ⅰ at different moments
时间t/s 柱端应变/με A/①轴柱 A/②轴柱 A/③轴柱 B/①轴柱 B/②轴柱 B/③轴柱 C/①轴柱 C/②轴柱 C/③轴柱 上端 下端 上端 下端 上端 下端 上端 下端 上端 下端 上端 下端 上端 下端 上端 下端 上端 下端 6.14 25 24 49 18 26 16 13 12 18 16 10 17 8 10 9 8 9 11 7.20 21 20 41 15 22 14 11 10 15 13 9 14 7 8 8 7 8 9 14.28 83 74 144 62 73 46 41 35 53 49 30 53 23 30 28 25 28 35 17.13 69 62 124 50 63 39 37 31 47 43 27 47 20 26 25 21 25 31 20.12 79 70 138 54 73 44 39 32 50 45 31 50 23 27 26 22 28 32 24.32 51 45 103 30 48 29 32 25 39 35 23 38 17 21 19 16 21 25 表 3 不同时刻模型Ⅱ底层柱端应变
Table 3. The strain of the columns of model Ⅱ at different moments
时间t/s 柱端应变/με A/①轴柱 A/②轴柱 A/③轴柱 B/①轴柱 B/②轴柱 B/③轴柱 C/①轴柱(带剪力墙) C/②轴柱 C/③轴柱(带剪力墙) 上端 下端 上端 下端 上端 下端 上端 下端 上端 下端 上端 下端 上端 下端 上端 下端 上端 下端 12.33 23 31 45 29 14 46 29 36 28 39 35 38 3 56 9 31 2 68 14.73 35 50 67 45 23 67 43 52 40 57 52 55 2 79 14 47 1 100 19.34 33 46 62 42 22 63 42 51 39 55 50 53 1 78 14 45 1 96 21.29 51 81 99 74 32 105 65 79 60 84 77 80 1 117 25 71 2 146 25.02 44 61 82 59 27 85 51 62 48 67 61 65 3 95 17 55 1 115 26.43 53 81 99 72 34 104 62 75 57 81 76 79 4 113 22 67 1 139 表 4 底层框架柱计算参数
Table 4. The column parameters of two models
编号 h/mm y/mm Iz/mm4 等效计算方法示意 A/①轴柱,A/③轴柱 380 40.00 5.12×106 A/②轴柱 380 40.00 3.41×106 B/①轴柱,B/③轴柱 740 40.00 5.12×106 B/②轴柱 740 40.00 3.41×106 C/②轴柱 740 40.00 3.41×106 C/①轴柱,C/③轴柱 740 40.00 3.41×106 C/①轴柱,C/③轴柱(带剪力墙) 740 88.22 3.94×107 注:表中右图为横墙对柱纵向抗侧刚度影响的等效计算方法示意,用于估算模型Ⅰ和模型Ⅱ的A/①轴柱、A/③轴柱、B/①轴柱、B/③轴柱与横墙组合体的惯性矩。 表 5 不同时刻模型Ⅰ底层柱地震剪力
Table 5. The shear force of the columns of model Ⅰ at different moments
时间t/s 地震剪力/N A/①轴柱 A/②轴柱 A/③轴柱 B/①轴柱 B/②轴柱 B/③轴柱 C/①轴柱 C/②轴柱 C/③轴柱 6.14 421 384 361 111 100 119 53 50 59 7.20 352 321 309 93 82 102 44 44 50 14.28 1 349 1 179 1 022 336 299 366 156 156 185 17.13 1 126 995 876 300 265 327 135 135 164 20.12 1 280 1 098 1 005 314 279 357 147 141 176 24.32 824 761 661 252 217 269 112 103 135 表 6 不同时刻模型Ⅱ底层柱地震剪力
Table 6. The shear force of the columns of model Ⅱ at different moments
时间t/s 地震剪力/N A/①轴柱 A/②轴柱 A/③轴柱 B/①轴柱 B/②轴柱 B/③轴柱 C/①轴柱
(带剪力墙)C/②轴柱 C/③轴柱
(带剪力墙)12.33 464 424 516 287 197 322 908 118 1 077 14.73 730 640 774 419 285 472 1 247 179 1 554 19.34 679 594 730 410 276 454 1 216 173 1 493 21.29 1 134 990 1 177 635 423 693 1 816 282 2 278 25.02 902 806 962 499 338 556 1 508 212 1 785 26.43 1 151 978 1 185 604 406 684 1 801 261 2 155 -
顾祥林, 黄庆华, 汪小林等, 2012. 地震中钢筋混凝土框架结构倒塌反应的试验研究与数值仿真. 土木工程学报, 45(9): 36—45Gu X. L. , Huang Q. H. , Wang X. L. , et al. , 2012. Spatial collapse responses of reinforced concrete frame structures under earthquakes-test study and numerical simulation. China Civil Engineering Journal, 45(9): 36—45. (in Chinese) 郭迅, 2010. 汶川地震震害与抗倒塌新认识. 土木建筑与环境工程, 32(S2): 28—29. 李小军, 于爱勤, 甘朋霞等, 2008. 汶川8.0级地震北川县城区灾害调查与分析. 震灾防御技术, 3(4): 352—362 doi: 10.3969/j.issn.1673-5722.2008.04.004Li X. J. , Yu A. Q. , Gan P. X. , et al. , 2008. Survey and analysis of the disaster and engineering damage of Beichuan county seat in MS8.0 Wenchuan earthquake. Technology for Earthquake Disaster Prevention, 3(4): 352—362. (in Chinese) doi: 10.3969/j.issn.1673-5722.2008.04.004 李英民, 罗文文, 韩军, 2013. 钢筋混凝土框架结构强震破坏模式的控制. 土木工程学报, 46(5): 85—92 doi: 10.15951/j.tmgcxb.2013.05.021Li Y. M. , Luo W. W. , Han J. , 2013. Control of failure mechanism for RC frame structure under strong earthquakes. China Civil Engineering Journal, 46(5): 85—92. (in Chinese) doi: 10.15951/j.tmgcxb.2013.05.021 林超, 郭子雄, 黄群贤等, 2018. 足尺砌体填充墙RC框架抗震性能试验研究. 建筑结构学报, 39(9): 30—37Lin C. , Guo Z. X. , Huang Q. X. , et al. , 2018. Experimental study on seismic behavior of full-scale infilled RC frames. Journal of Building Structures, 39(9): 30—37. (in Chinese) 刘伯权, 苏佶智, 马煜东等, 2019. 多层多跨钢筋混凝土平面框架拟静力倒塌试验研究. 土木工程学报, 52(8): 24—39 doi: 10.15951/j.tmgcxb.2019.08.003Liu B. Q. , Su J. Z. , Ma Y. D. , et al. , 2019. Pseudo-static collapse experiment of multi-story multi-span reinforced concrete plane frames. China Civil Engineering Journal, 52(8): 24—39. (in Chinese) doi: 10.15951/j.tmgcxb.2019.08.003 骆欢, 杜轲, 孙景江等, 2017. 地震作用下钢筋混凝土框架结构倒塌全过程振动台试验研究. 建筑结构学报, 38(12): 49—56Luo H. , Du K. , Sun J. J. , et al. , 2017. Shaking table test on complete collapse process of RC frame structure subjected to earthquake. Journal of Building Structures, 38(12): 49—56. (in Chinese) 曲哲, 和田章, 叶列平, 2011. 摇摆墙在框架结构抗震加固中的应用. 建筑结构学报, 32(9): 11—19Qu Z. , He T. Z. , Ye L. P. , 2011. Seismic retrofit of frame structures using rocking wall system. Journal of Building Structures, 32(9): 11—19. (in Chinese) 苏幼坡, 张玉敏, 王绍杰等, 2009. 从汶川地震看提高建筑结构抗倒塌能力的必要性和可行性. 土木工程学报, 42(5): 25—32 doi: 10.15951/j.tmgcxb.2009.05.010Su Y. P. , Zhang Y. M. , Wang S. J. , et al. , 2009. The necessity and feasibility of enhancing seismic design of structures based on the Wenchuan earthquake. China Civil Engineering Journal, 42(5): 25—32. (in Chinese) doi: 10.15951/j.tmgcxb.2009.05.010 王波, 2017. 典型RC框架结构地震倒塌机理研究. 三河: 防灾科技学院.Wang B., 2017. Experimental study on collapse mechanism of typical RC frame structures. Sanhe: Institute of Disaster Prevention. (in Chinese) 王波, 郭迅, 宣越, 2020. 基于新视角的震害分析——以四川长宁6.0级地震为例. 震灾防御技术, 15(3): 496—509 doi: 10.11899/zzfy20200303Wang B. , Guo X. , Xuan Y. , 2020. Analysis of earthquake damage from a new perspective: a case study of Changning M6.0 earthquake in Sichuan Province. Technology for Earthquake Disaster Prevention, 15(3): 496—509. (in Chinese) doi: 10.11899/zzfy20200303 王中阳, 车佳玲, 张尚荣等, 2018. 基于能量方法设计的RC框架结构易损性分析. 震灾防御技术, 13(3): 524—533Wang Z. Y. , Che J. L. , Zhang S. R. , et al. , 2018. Seismic fragility analysis of RC frame structure based on energy balance. Technology for Earthquake Disaster Prevention, 13(3): 524—533. (in Chinese) 徐超, 刘健, 温增平等, 2016. 汶川地震中某RC框架结构的震害分析与模拟. 震灾防御技术, 11(3): 542—551 doi: 10.11899/zzfy20160310Xu C. , Liu J. , Wen Z. P. , et al. , 2016. Damage analysis and simulation of a RC frame structure in the Wenchuan earthquake. Technology for Earthquake Disaster Prevention, 11(3): 542—551. (in Chinese) doi: 10.11899/zzfy20160310 许卫晓, 孙景江, 杜轲等, 2015. 阶梯墙框架结构体系的抗震性能分析. 工程力学, 32(2): 139—146Xu W. X. , Sun J. J. , Du K. , et al. , 2015. Seismic performance analysis of stepped wall-frame structure. Engineering Mechanics, 32(2): 139—146. (in Chinese) 闫培雷, 孙柏涛, 张昊宇, 2014. 芦山7.0级强烈地震钢筋混凝土框架结构教学楼震害. 土木工程学报, 47(S1): 24—28Yan P. L. , Sun B. T. , Zhang H. Y. , 2014. Seismic damage to RC frame teaching buildings in Lushan MS7.0 earthquake. China Civil Engineering Journal, 47(S1): 24—28. (in Chinese) 杨伟松, 郭迅, 许卫晓等, 2015. 增设翼墙RC框架结构的抗震性能研究. 振动与冲击, 34(24): 144—152Yang W. S. , Guo X. , Xu W. X. , et al. , 2015. Seismic behaviors of RC frames retrofitted with wing walls. Journal of Vibration and Shock, 34(24): 144—152. (in Chinese) 袁新星, 2018. RC框架结构“强柱弱梁”问题的探讨. 三河: 防灾科技学院.Yuan X. X. , 2018. Discussion on “strong column weak beam” of RC frame structure. Sanhe: Institute of Disaster Prevention. (in Chinese) 张望喜, 庞博, 徐帅, 2019. 基于IDA和易损性分析的RC框架结构底层柱端弯矩增大系数. 重庆大学学报, 42(5): 37—45Zhang W. X. , Pang B. , Xu S. , 2019. Research on bottom column end moment amplification factor of RC frame based on incremental dynamic analysis and fragility analysis. Journal of Chongqing University, 42(5): 37—45. (in Chinese) 中华人民共和国住房和城乡建设部, 中华人民共和国国家质量监督检验检疫总局, 2010. GB 50011—2010 建筑抗震设计规范. 北京: 中国建筑工业出版社.Ministry of Housing and Urban-Rural Development of the People’s Republic of China, General Administration of Quality Supervision, Inspection and Quarantine of the People’s Republic of China, 2010. GB 50011—2010 Code for seismic design of buildings. Beijing: China Architecture & Building Press. (in Chinese) 周晓洁, 李忠献, 续丹丹等, 2015. 柔性连接填充墙框架结构抗震性能试验. 天津大学学报(自然科学与工程技术版), 48(2): 155—166Zhou X. J. , Li Z. X. , Xu D. D. , et al. , 2015. Experiment on seismic behavior of flexible connection masonry infilled frame structure. Journal of Tianjin University (Science and Technology), 48(2): 155—166. (in Chinese) Penava D. , Sarhosis V. , Kožar I. , et al. , 2018. Contribution of RC columns and masonry wall to the shear resistance of masonry infilled RC frames containing different in size window and door openings. Engineering Structures, 172: 105—130. doi: 10.1016/j.engstruct.2018.06.007 Peng Q. M. , Zhou X. J. , Yang C. H. , 2018. Influence of connection and constructional details on masonry-infilled RC frames under cyclic loading. Soil Dynamics and Earthquake Engineering, 108: 96—110. doi: 10.1016/j.soildyn.2018.02.009 -