Hazards and Their Effects of Earthquakes with MS 6.0 and Above Since 2010 in Qinghai Province
-
摘要: 近十多年来中国大陆强震主要围绕青藏高原内部块体边界发育，如2008年汶川8.1级地震、2010玉树7.1级地震、2021年玛多7.4级地震、2022年门源6.9级地震等。其中青海省位于青藏高原东北部，块体构造复杂，变形强烈，自2010年以来6.0级及以上地震共计发生8次，占中国大陆6.0级及以上地震约22%，为地震频发区。本文基于地震现场资料、应急管理部、中国地震台网中心、美国地质勘探局（USGS）等不同信息源，对2010年以来青海6.0级及以上地震灾害情况进行汇总，并对主要震例震害特征及其抗震设防特点进行分析讨论，为后续青海地区震害防御提供相关建议。Abstract: In the past decade, strong earthquakes in mainland China have mainly developed around the boundary of the Tibetan Plateau, including Wenchuan 8.1 earthquake in 2008, Yushu 7.1 earthquake in 2010, Maduo 7.4 earthquake in 2021, Menyuan 6.9 earthquake in 2022, etc. Qinghai province is located in the northeast of the Qinghai-Tibet plateau, with complex block structure and strong deformation. Since 2010, a total of 8 earthquakes with magnitude 6.0 or above occurred, accounting for about 22% of the earthquakes of magnitude 6 or above in mainland China. In this paper,we collected the earthquake data from field investigation,as well as from the Ministry of Emergency Management, the China Earthquake Networks Center, the United States Geological Survey and other different information sources, summarized the earthquake hazards of earthquakes with magnitude 6.0 or above since 2010 in Qinghai region, and analyzed the characteristics of major earthquake damage and seismic fortification to provide relevant suggestions for seismic disaster prevention in Qinghai in the future.1)
2 2 青海省地震局，2016. 2016年1月21日青海门源6.4级地震灾害直接经济损失评估结报告（青海灾区）.2)4 3 青海省地震局，2016. 2016年10月17日青海杂多6.2级地震灾害直接经济损失评估结报告.3)1 4 应急管理部救灾和物资保障司，2021. 应急管理部发布2021年上半年全国自然灾害情况.4)3 5 应急管理部救灾和物资保障司，2022. 应急管理部发布2022年1月全国自然灾害情况. -
引言
高强钢有着其它钢材所不具备的高强度、高延性以及高耗能特征，且其自重小，回收利用率高，在大跨度、超高层以及不规则结构中得到了广泛的应用（Fukumoto等，1985）。目前，国内外学者对高强钢的研究也越来越深入，甚至已经出版了与之对应的设计规范（Rasmussen等，1992；国际桥梁与结构工程协会，2010）。自20世纪60年代至今，国外学者进行了很多关于高强钢压弯构件稳定性和残余应力影响的研究工作，通过试验验证发现：对于残余应力的分布情况，高强钢和普通钢大致相同，材料种类对其的影响极小。Usami等（1984）对25根不同长度、宽厚比和截面形式的高强钢构件进行了试验研究，并结合理论分析研究了梁、柱局部稳定性的问题。Rasmussen等（1995）为澳大利亚的钢结构设计找到1条合理的柱子骨架曲线，对一批屈服强度为690MPa的“工”字型及“箱”型焊接柱进行了轴压及压弯试验研究。而国内对于Q460的研究主要集中于构件的力学性能和稳定性承载力方面，也有一些学者考虑了残余应力的影响（班慧勇等，2011）。施刚等（2012）对一批Q460“工”字型钢柱进行了低周反复荷载试验，结果表明，当翼缘和腹板宽厚比超过限值时，试件的承载力及其退化程度并没有明显的变化，且构件的极限位移角均明显大于规定的1/50（0.02）。这说明Q460钢构件极限承载力较好，具有良好的抗震性能。建议在实际设计中，当轴压比较大时，适当的减小板件宽厚比。李国强等（2013）对一批轴压比为0.3的Q460焊接“工”字型和“箱”型柱进行了低周反复荷载试验，分析表明宽厚比对构件的滞回性能有较大的影响。
虽然对高强钢构件的研究已经取得了较大的进展，但仍有诸多方面需要完善，特别是国内的研究起步相对较晚。目前《建筑抗震设计规范》（GB 50011-2010）（中华人民共和国住房和城乡建设部等，2010）中对于翼缘宽厚比、腹板高厚比及箱型截面壁板宽厚比的限值仅适用于普通钢材Q235，对于Q460及以上高强钢在实际工况中的应用仍没有明确的理论与试验支撑。
基于目前研究现状，以Q460钢为材料，对框架中的主要承重构件进行了精细化的有限元模型分析。根据现行的《建筑抗震设计规范》，以轴压比、翼缘宽厚比、腹板高厚比和壁板宽厚比为变量，设计了一系列“工”字型框架柱和“箱”型框架柱，并对其进行了循环往复荷载作用下的模拟试验，分析了框架柱的抗震性能，研究了其承载能力、耗能能力、变形能力及刚度退化程度，以期对Q460钢在实际工程中的应用提供一些参考。
1. 有限元计算模型
采用ABAQUS有限元软件建立Q460柱的分析模型，并用已有的实验数据验证模型的准确性，对构件的抗震性能进行分析。
1.1 有限元模型建立
在水平地震荷载下，柱的弯矩随荷载方向的变化而发生周期性的变化，反弯点一般出现在底层柱高约2/3位置处。因此，截取反弯点以下一部分为子结构，作为柱试验的研究构件，并简化为底端固接、顶部可以水平移动的悬臂端。计算简图如图 1（a）所示。
本文采取常用的“工”字型截面和“箱”型截面柱进行有限元建模分析，“工”字型和“箱”型截面以及截面中的尺寸b、h0、bf、t、tf、tw、B和D的含义如图 1（b）、（c）所示。
有限元模型所使用的Chaboche钢材塑性本构模型采用了Von Mises流动法则，如图 2所示，该模型包含了非线性随动强化部分和等向强化部分，是1种混合强化模型（Chaboche，1986）。
等向强化部分定义了屈服面${\sigma ^0}$的大小，表示如下：
$$ {\sigma ^0} = \sigma \left| {_0} \right. + {Q_\infty }\left({1 - {{\rm{e}}^{ - b{\varepsilon ^{ - pl}}}}} \right) $$ (1) 其中，$\sigma \left| {_0} \right.$为等效塑性应变；-pl为零时的材料应力；${Q_\infty }$为屈服面的最大变化值；b为屈服面大小随塑性应变增加的变化率；基于滞回环数据，通过曲线拟合可得材料参数${Q_\infty }$和b。
模型的随动强化部分则定义了背应力${\alpha _k}$，公式如下：
$$ {\alpha _k} = \frac{{{C_k}}}{{{\gamma _k}}}\left({1 - {{\rm{e}}^{ - {\gamma _k}{\varepsilon ^{pl}}}}} \right) + {\alpha _{k, 1}}{{\rm{e}}^{ - {\gamma _k}{\varepsilon ^{pl}}}} $$ (2) $$ \alpha = \sum\limits_{k = 1}^n {{\alpha _k}} $$ (3) 其中，n是背应力的个数，Ck和${\gamma _k}$是通过循环试验数据校准的材料参数，Ck是初始随动强化模量，${\gamma _k}$是随动强化模量随着塑性变形增加而减小的变化率，${\alpha _{k, 1}}$表示在第1个数据点的第k个背应力（第K个背应力的初始值）。为满足精度要求，采用式（3）所示的多背应力叠加方式可得到更为准确的拟合曲线。
试验模拟所需的Q460参数取自石永久课题组在钢材试验基础上所得到的本构参数（王萌等，2013），Q235和Q345相关参数取自文献（王萌，2013）。Q460、Q345和Q235钢材本构参数如表 1所示。
表 1 钢材本构模型参数Table 1. Steel constitutive model parameters钢材材性 $\sigma \left| {_0} \right.$/ N·mm-2 ${C_{{\rm{k}},1}}$/ N·mm-2 ${\gamma _1}$ ${C_{{\rm{k}},2}}$/ N·mm-2 ${\gamma _2}$ ${C_{{\rm{k}},3}}$/ N·mm-2 ${\gamma _3}$ ${C_{{\rm{k}},4}}$/ N·mm-2 ${\gamma _4}$ Q460 474 4797 156 3794 145 1498 107 Q345 429 7993 175 6773 116 2854 34 1450 29 Q235 407 6013 173 5024 120 3026 32 990 35 应用表 1中的参数在ABAQUS中的Combined模块和Cycle hardening材料属性进行循环荷载下的模拟试验。
模型均采用八节点六面体单元，并设置为当单元材料达到损伤破坏后将会被自动删除。网格划分采用结构化网格划分技术，用网格试验的方法来确定合理的网格密度，有限元模型网格划分图如图 3所示。研究的柱为从门式刚架反弯点中截取的部分子结构，边界条件假定为：柱脚固接，固定其6个自由度；柱顶为沿X方向为可以自由移动的滑动支座，在柱顶截面中心的正上方20mm处设置耦合点施加轴向力和横向水平位移，并由该耦合点控制柱顶表面所有点的自由度和刚性。在施加轴力的同时，采用侧向位移控制的方法施加循环荷载，观察柱的受力性能。
图 3 “工”字型柱、“箱”型柱网格划分Figure 3. Division diagram of I-type column and box-type column grid根据图 4所示的加载示意图，轴向力N在第一分析步时施加竖向荷载，构件产生初始轴向变形；施加沿横向的水平位移，最大位移设置为0.05—0.1倍的柱高。
1.2 有限元模型验证
对于Q460钢材，按《建筑抗震设计规范》对钢材的翼缘[b/tf]和腹板宽厚比限值[h0/tw]进行计算，表 2分别为根据规范换算得到在一级抗震设防情况下的限值。表 3为位移荷载增量，屈服前第一级位移荷载取为0.5δy，循环1圈；从第二级开始每1级循环3圈，且位移大小分别为δy，2δy，3δy，4δy，……，直至构件发生破坏。
表 2 试件参数Table 2. Test piece parameters试件编号 H/mm B/mm tf/mm tw/mm b/tf h0/tw [b/tf] [h0/tw] n N/kN H-0-1 200 150 12 10 5.8 17.6 7.1 30.7 0.2 493.1 H-0-2 300 180 12 10 7.1 27.6 7.1 30.7 0.2 651.4 H-0-3 300 220 12 10 8.8 27.6 7.1 30.7 0.2 739.7 表 3 试件的位移荷载Table 3. Displacement load of specimen试件编号 ${\delta _1}$/mm ${\delta _2}$/mm ${\delta _3}$/mm ${\delta _4}$/mm ${\delta _5}$/mm ${\delta _6}$/mm ${\delta _7}$/mm H-0-1 9 18 36 54 72 90 108 H-0-2 6 12 24 36 48 60 72 H-0-3 6 12 24 36 48 60 72 本文对施刚等（2012）的试验对象进行数值模拟，并将试验结果与计算结果对比，验证有限元模型的正确性。模拟结果与试验的对比如图 5—7所示。
从图 5—7和表 4的误差分析可以看出，模拟结果所得曲线和试验曲线基本一致，且模拟结果的饱满度更高。通过结果对比分析可知，采用本文所述的数值模型以及材料参数，可以较好地模拟Q460构件在循环荷下的受力性能。
表 4 ABAQUS计算与试验结果的最大承载力对比Table 4. Comparison of maximum bearing capacity between calculation and test构件编号 ABAQUS结果/kN 试验结果/kN 相差/% H-0-1 132 139 -5 H-0-2 232 222 5.1 H-0-3 305 294 3.7 2. 试件设计及加载工况
为研究Q460柱的抗震性能，本文以轴压比、翼缘宽厚比以及腹板高厚比为变量，对工程实际中常用到的“工”字型与“箱”型截面柱进行数值模拟分析。轴压比变化范围0.1—0.6。对柱截面壁板的限值，以《建筑抗震设计规范》中一级抗震等级为例，“工”字型截面翼缘宽厚比不应超过10，腹板高厚比不应超过43，“箱”型截面壁板宽厚比的限值为33，换算为对应于Q460钢材时的限值分别为7.1、30.7、23.6，考虑到Q460钢材刚度大，所以设计中翼缘宽厚比的限值取到9，对于长细比而言不应大于42。“工”字型、“箱”型截面试件设计参数如表 5、表 6所示。
表 5 “工”字型截面试件设计参数Table 5. Design parameters of I-type specimen试件编号 b/mm h0/mm tw/mm tf/mm b/tf h0/tw λ L/mm n N/N H-1-1 178 250 10 12 7 25 30 1213 0.1 311.5 H-1-2 178 250 10 12 7 25 30 1213 0.2 623.0 H-1-3 178 250 10 12 7 25 30 1213 0.3 934.5 H-1-4 178 250 10 12 7 25 30 1213 0.4 1246.0 H-1-5 178 250 10 12 7 25 30 1213 0.5 1557.6 H-1-6 178 250 10 12 7 25 30 1213 0.6 1869.1 H-2-1 130 250 10 12 5 25 30 1183 0.4 1034.1 H-2-2 154 250 10 12 6 25 30 1199 0.4 1140.1 H-2-3 202 250 10 12 8 25 30 1224 0.4 1352.0 H-2-4 226 250 10 12 9 25 30 1234 0.4 1458.0 H-3-1 178 230 10 12 7 23 30 1088 0.4 1209.2 H-3-2 178 280 10 12 7 28 30 1700 0.4 1301.2 H-3-3 178 300 10 12 7 30 30 1524 0.4 1338.0 H-3-4 178 320 10 12 7 32 30 1648 0.4 1374.8 H-4-1 130 250 10 12 5 25 30 1183 0.2 517.0 H-4-2 154 250 10 12 6 25 30 1199 0.2 570.0 H-4-3 202 250 10 12 8 25 30 1524 0.2 676.0 H-4-4 226 250 10 12 9 25 30 1534 0.2 729.0 H-5-1 178 230 10 12 7 23 30 1088 0.2 604.6 H-5-2 178 280 10 12 7 28 30 1400 0.2 650.6 H-5-3 178 300 10 12 7 30 30 1524 0.2 669.0 H-5-4 178 320 10 12 7 32 30 1648 0.2 687.4 H-6-1 130 250 10 12 5 25 30 1183 0.6 1551.1 H-6-2 154 250 10 12 6 25 30 1199 0.6 1710.1 H-6-3 202 250 10 12 8 25 30 1224 0.6 2028.0 H-6-4 226 250 10 12 9 25 30 1234 0.6 2187.0 H-7-1 178 230 10 12 7 23 30 1088 0.6 1813.9 H-7-2 178 280 10 12 7 28 30 1400 0.6 1951.9 H-7-3 178 300 10 12 7 30 30 1524 0.6 2007.1 H-7-4 178 320 10 12 7 32 30 1648 0.6 2062.3 表 6 “箱”型截面试件设计参数Table 6. Designing parameters of box-type specimen试件编号 B/mm D/mm t/mm B/t [B/t] λ L/mm n N/N X-1-1 300 300 20 15 23.6 30 1573 0.1 533.6 X-1-2 300 300 20 15 23.6 30 1573 0.2 1067.2 X-1-3 300 300 20 15 23.6 30 1573 0.3 1600.8 X-1-4 300 300 20 15 23.6 30 1573 0.4 2134.4 X-1-5 300 300 20 15 23.6 30 1573 0.5 2668 X-1-6 300 300 20 15 23.6 30 1573 0.6 3201.6 X-2-1 260 260 20 13 23.6 30 1288 0.4 1840 X-2-2 280 280 20 14 23.6 30 1430 0.4 1987.2 X-2-3 320 320 20 16 23.6 30 1716 0.4 2281.6 X-2-4 340 340 20 17 23.6 30 1859 0.4 2428.8 X-3-1 260 260 20 13 23.6 30 1288 0.2 920 X-3-2 280 280 20 14 23.6 30 1430 0.2 993.6 X-3-3 320 320 20 16 23.6 30 1716 0.2 1140.8 X-3-4 340 340 20 17 23.6 30 1859 0.2 1214.4 X-4-1 260 260 20 13 23.6 30 1288 0.6 2760 X-4-2 280 280 20 14 23.6 30 1430 0.6 2980.8 X-4-3 320 320 20 16 23.6 30 1716 0.6 3422.4 X-4-4 340 340 20 17 23.6 30 1859 0.6 3643.2 在对单一参数变化分析过程中，为避免不同加载方式对试验结果的影响，对所有构件均采用相同的位移加载制度，加载制度参考文献（施刚等，2012），如图 8所示。
3. 模拟结果与分析
以Q460钢柱的抗震性能为研究目标，依据现行《建筑抗震设计规范》，设计并完成了30个“工”字型柱及18个“箱”型柱的模拟试验，通过提取设置在柱顶截面正上方20mm处的耦合点处的模型分析数据，对不同参数下Q460钢柱的承载力、耗能能力、延性特征等抗震性能进行了分析。
3.1 滞回分析
各试件的滞回曲线如图 9和图 10所示。由图可知，所有的试件曲线形状均符合梭形的特点。滞回环形状饱满，表明柱具有良好的变形性能和耗能能力。在加载初期，构件处于弹性阶段，变形可以完全恢复，滞回曲线并不饱满，此时滞回环包围的面积很小；在加载中期，即塑性发展阶段，滞回环十分饱满，耗能增加，卸载后出现残余变形；加载后期（退化阶段）滞回曲线又变得不太饱满，表明钢框架柱的耗能主要发生在加载中期的塑性发展阶段，且存在强度衰减和刚度退化现象。试件的耗能能力通常以滞回曲线滞回环包围的面积来衡量，图 11和图 12为构件的能量耗散系数（E）-循环次数（n）曲线对比。
从图中11和图 12可以看出，随着位移及加载次数的增加，耗散能量系数逐渐增大，且在加载后期，能量耗散系数达到3左右。从图 11可以看出“工”字型柱随着轴压比的增大，同一滞回环下的能量耗散系数越小，但影响不明显。不同翼缘宽厚比的试件同一滞回环下的能量耗散系数基本一致。不同腹板高厚比对能量耗散系数有很明显的影响，当腹板高厚比越大时，能量耗散系数越小。从图 12可以得到轴压比对“箱”型柱与“工”字型柱的影响规律是一致的，但壁板宽厚比越大，能量耗散系数越小，且影响较为明显，分析原因是由于箱型柱主要靠壁板来耗散能量。
3.2 承载力分析
骨架曲线基本上能直观地反映出结构或构件的屈服荷载和屈服位移、极限荷载和极限位移等特征点，也能清晰地反映构件的强度、刚度和变形性能等。Q460柱的骨架曲线如图 13和图 14所示，构件在往复荷载作用下均经历了弹性、屈服和极限3个阶段，正反向与骨架曲线基本一致。
根据骨架曲线可以得到构件各特征点对应的荷载值，从图 13（a）和图 14（a）可以明显看出，极限承载力随轴压比增大而变小，且之后下降趋势一致。从图 13（b）可以看出，随着翼缘宽厚比的增大，承载力会随之增大，但构件的刚度退化较为严重，最终不同翼缘宽厚比的骨架曲线会趋于一点。当翼缘宽厚比到达某一值时，翼缘宽厚比的影响变小，而其它影响因素的敏感度会增加。从图 13（c）可以看出，当轴压比n为0.4，翼缘宽厚比为7时，达到最大承载力后，承载力随着腹板高厚比的增大，下降趋势也越明显。当水平位移均为20mm左右时，构件均达到极限荷载。腹板高厚比大小越接近规范限值，承载力的下降趋势更加明显，因此构件延性性能也越差，已不能满足弹塑性变形的要求，因此在设计的过程中，建议腹板高厚比不宜过大。对于壁板宽厚比对“箱”型截面柱承载力的影响，从图 14（b）可以看出，随着壁板宽厚比的增加，箱型柱的极限承载力随之增大，且之后下降趋势一致。
3.3 延性特征分析
延性是反映结构、构件或材料非弹性变形能力的1个度量指标。延性的大小可以通过位移延性系数、屈服位移转角、最大位移转角以及极限位移转角来衡量（郭子雄等，2004），本文对于所有构件采用位移角和延性系数作为性能指标。构件破坏的时刻应该为最大承载力下降至0.85Pmax时，此时所对应的位移角为极限位移角${\theta _u}$，即位移与柱高的比值$\left({{\mathit{\Delta }_u}/L} \right)$。因各构件骨架曲线无明显屈服点，故采用如下方法进行确定：从原点中引出OB切线与极限荷载Pmax相交于点B，做AB直线与位移-荷载曲线交于A点，连接OA与极限荷载Pmax相交于C点，C点所对应的位移就是屈服位移。E点就是屈服点。计算简图如图 15所示。
表 7为所有模拟试件的延性比取值情况，图 16和17为不同轴压比、翼缘宽厚比、腹板高厚比和壁板宽厚比下的极限位移角与屈服位移角取值情况。对比分析可以看出，所有构件的延性系数均大于1，但相对于Q235柱的延性系数可以达到4，表明Q460柱虽然有较高的承载力，但是不能承受较大的塑性变形。屈服位移角与极限位移角有较大的差异，表明Q460钢柱在屈服之后仍然有良好的变形能力。现行《建筑抗震设计规范》中规定钢结构弹塑性位移角达到1/50时，结构或构件发生破坏。但是，该限值是通过大量混凝土柱试验得到的，对于钢结构，甚至高强钢并不一定适用。下面就“工”字型和“箱”型2种截面柱的延性和极限位移角分别进行分析，并给出极限位移角的建议值。
表 7 试件延性比计算结果Table 7. Ductility ratio result of specimen from calculation截面形状 试件编号 延性比 试件编号 延性比 试件编号 延性比 “工”字型 H-1-1 2.5256 H-3-1 2.5495 H-5-3 1.8175 H-1-2 2.5596 H-3-2 1.9412 H-5-4 1.8245 H-1-3 2.6382 H-3-3 2.4116 H-6-1 1.9092 H-1-4 2.5034 H-3-4 2.2301 H-6-2 2.0477 H-1-5 2.4784 H-4-1 1.9758 H-6-3 2.6084 H-1-6 1.3040 H-4-2 1.9735 H-6-4 1.9906 H-2-1 2.8007 H-4-3 1.5396 H-7-1 2.4463 H-2-2 2.8804 H-4-4 1.5725 H-7-2 1.8307 H-2-3 2.6508 H-5-1 2.0780 H-7-3 1.6720 H-2-4 2.3611 H-5-2 1.8176 H-7-4 1.7206 “箱”型 X-1-1 2.8794 X-2-1 2.4012 X-3-3 1.9370 X-1-2 2.7504 X-2-2 1.6506 X-3-4 1.8464 X-1-3 2.4262 X-2-3 1.8997 X-4-1 2.0211 X-1-4 2.1488 X-2-4 1.6812 X-4-2 1.8577 X-1-5 2.0218 X-3-1 2.0817 X-4-3 1.6825 X-1-6 1.6233 X-3-2 2.1513 X-4-4 1.9984 
图 16 “工”字型柱位移角的统计结果Figure 16. Statistical results of displacement angle of I-type column对表 7和图 16综合对比分析，可以看出“工”字型截面柱的极限位移角和延性比随着轴压比的增加而减小，分析原因可知当轴压比增大时，限制了柱脚的转动能力。除H-1-3构件外，所有构件的极限位移角都远远超过规定的0.02（1/50），甚至有的可以达到0.04以上。极限位移角随着翼缘宽厚比的增加呈下降的趋势，分析原因为，在往复荷载作用下，当翼缘宽厚比增大，更易发生屈曲变形，导致柱子承受竖向荷载能力减弱而发生破坏。通过对表 7中H-2-1—H-2-4、H-4-1—H-4-4以及H-6-1—H-6-4的柱延性比的分析，表明翼缘宽厚比对于延性比的影响不明显。
而翼缘宽厚比限值的规定为7.1，考虑到Q460钢材的强度高、刚度大，本文通过试验表明当翼缘宽厚比为9时，依然有良好的变形能力，因此建议在抗震设计中，Q460钢柱翼缘宽厚比的限值取为9。极限位移角和延性比随着腹板宽厚比的增大呈下降趋势，分析原因为，腹板高厚比越大，水平荷载作用下柱更易发生整体失稳。因此，在抗震设计中，腹板高厚比选择应合理且不宜过大。
通过对极限位移角和延性比的分析，考虑到轴压比对于极限位移角的影响较大，将轴压比以0.2、0.4和0.6为界分为4个阶段，并且考虑翼缘以及腹板的综合影响，得出如表 8所示的极限位移角的取值范围。分析可得，极限位移角取值范围为[0.0193，0.0414]，平均值为0.0305，标准差为0.00565。建议取0.03作为Q460“工”字型柱构件破坏时的极限位移角。
表 8 “工”字型柱位移角取值Table 8. The recommended value of displacement angle of I-type columnn ${b_1}/{t_{\rm{f}}} \le 7\;且\;{h_0}/{t_{\rm{w}}} \le 28$ ${b_1}/{t_{\rm{f}}} > 7\;且\;{h_0}/{t_{\rm{w}}} > 28$ n≤0.2 θ＞0.035 0.035≥θ≥0.03 0.2＜n≤0.4 0.0375≥θ≥0.0325 0.035≥θ≥0.03 0.4＜n≤0.6 0.0325≥θ≥0.0275 0.03≥θ≥0.025 n＞0.6 0.03≥θ≥0.0275 0.0275≥θ≥0.025 对表 7和图 17对比分析可看出，“箱”型柱的极限位移角和延性比整体上随着轴压比的增加而减小，即轴压比增大，柱脚的转动能力变差。所有柱子的极限位移角均远远超过规范所规定的取值0.02（1/50），甚至有的可以达到0.04以上。并且随着翼缘宽厚比的增加，极限位移角有下降的趋势。
图 17 “箱”型柱位移角的统计结果Figure 17. Statistical results of displacement angle of box-type column通过对不同轴压和壁板宽厚比的极限位移角的分析，可以看出轴压比对于极限位移角的影响较大，将轴压比以0.2、0.4和0.6为界分为4个阶段，并且考虑壁板宽厚比的影响，得出如表 9所示的极限位移角的范围。分析可得，极限位移角取值范围为[0.0226，0.0426]，平均值为0.0321，标准差为0.0055，建议取0.032作为Q460“箱”型柱构件破坏时的极限位移角。
表 9 “箱”型柱层间位移角的取值Table 9. The recommended value of displacement angle of box-type columnn $B/t \le 15$ $B/t > 15$ n ≤0.2 θ＞0.035 0.035≥θ≥0.04 0.2＜n ≤0.4 0.04≥θ≥0.035 0.035≥θ≥0.03 0.4＜n ≤0.6 0.035≥θ≥0.03 0.03≥θ≥0.025 n＞0.6 0.03≥θ≥0.025 0.025≥θ≥0.02 3.4 刚度退化分析
构件刚度退化曲线可以反映构件累积损伤的影响。构件在往复荷载作用下，由于材料的屈服和累积损伤等影响，会出现刚度退化的现象。
“工”字型截面柱刚度退化情况如图 18—20所示，“工”字型截面柱随着荷载位移的加大刚度呈下降趋势，且下降趋势较一致。由图 18可以看出，在初始加载时刚度均为9左右，不同轴压比在弹性阶段对刚度退化影响较小，但是，当构件进入弹塑性直到破坏的整个过程中，轴压比越大的构件刚度退化越为严重。从图 19可以看出，在高、中、低3种轴压比下，在初始循环加载时，初始刚度相差较大，随着构件进入弹塑性阶段直至破坏，构件的刚度呈现不同程度的降低，当位移达到0.55时刚度均小于1，且相差很小，趋于同一刚度值。因此在实际工程中，使用高强钢时不宜采用翼缘宽厚比较小的构件。对图 20分析表明，对于不同腹板高厚比的柱，初始刚度有较小的差异且随着荷载的增加，刚度退化呈现不同程度的降低。对于腹板高厚比小的柱，虽初始刚度较高，但刚度退化较快；相反，腹板高厚比较大的构件，刚度退化比较缓慢。
图 18 “工”字型柱不同轴压比对应的刚度退化Figure 18. Stiffness degradation of I-type column under different axial compression ratios
图 19 “工”字型柱不同翼缘宽厚比对应的刚度退化Figure 19. Stiffness degradation of I-type column under different flange width-thickness ratio
图 20 “工”字型柱不同腹板高厚比对应的刚度退化Figure 20. Stiffness degradation of I-type column under different web thickness ratio“箱”型截面柱刚度退化情况如图 21和图 22所示，与“工”字型柱相比较，“箱”型截面柱有着较高的初始刚度，刚度退化缓慢，在塑性阶段仍有较高的刚度。
图 21 “箱”型柱不同轴压比对应刚度退化Figure 21. Stiffness degradation of box-type column under different axial compression ratios
图 22 “箱”型柱不同壁板宽厚比对应刚度退化Figure 22. Stiffness degradation of box-type column under different panel wide heavy ratios图 21表明，当其它条件保持一致时，“箱”型柱的刚度随着轴压比的增大退化更为严重，与“工”字型柱变化规律一致。图 22表明，随着壁板宽厚比的增大，“箱”型柱的刚度退化越缓慢，因此在抗震设计中，可以适当采用较大的壁板宽厚比，以避免更为严重的刚度退化现象。
4. 结论
本文利用ABAQUS有限元软件对48根不同轴压比、翼缘宽厚比、腹板高厚比、壁板宽厚比的“工”字型和“箱”型截面压弯构件进行了低周反复荷载作用下的抗震性能分析，得出如下结论：
（1）轴压比较大的构件，其能量耗散系数相对较小。对于“工”字型框架柱，不同翼缘宽厚比对能量耗散系数的影响较小，而不同腹板高厚比和壁板宽厚比的试件，腹板高厚比和壁板宽厚比越大，能量耗散系数越小，且影响较为明显。
（2）“工”字型和“箱”型柱的极限承载力随着轴压比的减小及壁板宽厚比的增大而增大，并且下降趋势一致。承载力随着翼缘宽厚比的增大而增大，最终骨架曲线趋于一点，当翼缘宽厚比到达某一值时，翼缘宽厚比的影响将变小，而其它影响因素的敏感度会增加。构件达到最大承载力后，随着腹板高厚比接近限值，承载力下降趋势变明显。
（3）相比Q235钢柱的延性系数可以达到4，Q460的延性系数仅在2左右，表明其在强震作用下虽然有较高的承载力，但是却不能够承受较大的塑性变形。当轴压比一定时，构件的延性系数会随壁板宽厚比的增大而减小。极限位移角随着轴压比、翼缘宽厚比和腹板高厚比的增大而减小，并且轴压比的影响较为明显。根据数值模拟结果分析给出极限位移角的建议值：“工”字型柱为0.03，“箱”型柱为0.032；为避免翼缘屈曲造成柱子失稳，规范规定翼缘宽厚比不应超过7.1。考虑到高强钢的碳含量高、刚度较大，本文通过模拟试验给出的建议值为9。
（4）构件刚度随轴压比的增大，退化程度加剧。随翼缘宽厚比的增大，构件初始刚度差值增大，刚度退化现象严重但最终趋于同一刚度，因此在实际应用中，建议Q460强钢柱不宜用于翼缘宽厚比过小的构件。对于腹板高厚比越小的柱，虽然初始刚度较高，但刚度退化相对较快，因此在实际应用中，要综合考虑初始刚度以及刚度退化的因素，选择腹板高厚比适宜的构件。不同壁板宽厚比的“箱”型柱初始刚度差别不大，且随着宽厚比的增大刚度退化越缓慢，抗震设计中，可以适当采用较大的壁板宽厚比，以避免更为严重的刚度退化现象。
-
图 2 2010年以来青海6.0级及以上地震M-T图
Figure 2. M-T map of Qinghai earthquakes of magnitude 6.0 or above since 2010
图 8 主要震例中国地震动区划参数对比
Figure 8. Comparison of main earthquake examples-china earthquake zoning parameters
表 1 2010年以来青海6.0级及以上地震
Table 1. List of the Qinghai earthquakes of magnitude 6.0 or above since 2010
日期 北京
时间经度
/（°）纬度/（°） 震级/级 震源深度/km 震中位置 人员伤亡 2010-04-14 07:49 96.7 33.1 7.1 14 玉树藏族自治州玉树市 2698人遇难，270人失踪，
约10000人不同程度受伤2010-04-14 09:25 96.6 33.2 6.3 30 玉树藏族自治州玉树市 玉树地震余震 2010-10-17 15:14 96.6 33.1 6.2 9 玉树藏族自治州玉树市 玉树地震余震 2016-01-21 01:13 101.62 37.65 6.4 10 海北藏族自治州门源县 9人不同程度受伤 2016-10-17 15:14 94.93 32.81 6.2 9 玉树藏族自治州杂多县 1人遇难（地震惊吓） 2021-05-22 02:04 98.34 34.59 7.4 17 果洛藏族自治州玛多县 19人不同程度受伤 2022-01-08 01:45 101.26 37.77 6.9 10 海北藏族自治州门源县 10人不同程度受伤 2022-03-26 00:21 97.33 38.50 6.0 10 海西藏族自治州德令哈市 无人员伤亡，震中周边50 km未有村庄分布 
下载: 导出CSV
-
曹生奎, 刘峰贵, 张海峰等, 2005. 青海高原地震重灾区的灾害特点及成因探析. 灾害学, 20（1）: 76—79Cao S. K. , Liu F. G. , Zhang H. F. , et al. , 2005. An analysis on characteristics and causes of earthquake disaster in heavy disastered area in the Qinghai Plateau. Journal of Catastrophology, 20(1): 76—79. (in Chinese) 陈华静, 李一行, 宫玥等, 2020. 地震灾害风险防治协同机制研究. 震灾防御技术, 15（4）: 731—738Chen H. J. , Li Y. H. , Gong Y. , et al. , 2020. Research on the cooperative mechanism of earthquake disaster risk prevention and control. Technology for earthquake Disaster Prevention, 15(4): 731—738. (in Chinese) 盖海龙, 李智敏, 姚生海等, 2022.2022年青海门源MS6.9地震地表破裂特征的初步调查研究. 地震地质, 44（1）: 238—255Gai H. L. , Li Z. M. , Yao S. H. , et al. , 2022. Preliminary investigation and research on surface rupture characteristics of the 2022 Qinghai Menyuan MS6.9 earthquake. Seismology and Geology, 44(1): 238—255. (in Chinese) 高孟潭, 2021. 中国地震区划技术的发展与展望. 城市与减灾, （4）: 7—12Gao M. T. , 2021. Development and prospect on technology of seismic zoning in China. City and Disaster Reduction, (4): 7—12. (in Chinese) 郭迅, 2009. 汶川大地震震害特点与成因分析. 地震工程与工程振动, 29（6）: 74—87Guo X. , 2009. Characteristics and mechanism analysis of the great Wenchuan earthquake. Journal of Earthquake Engineering and Engineering Vibration, 29(6): 74—87. (in Chinese) 胡朝忠, 杨攀新, 李智敏等, 2016.2016年1月21日青海门源6.4级地震的发震机制探讨. 地球物理学报, 59（5）: 1637—1646Hu C. Z. , Yang P. X. , Li Z. M. , et al. , 2016. Seismogenic mechanism of the 21 January 2016 Menyuan, Qinghai MS6.4 earthquake. Chinese Journal of Geophysics, 59(5): 1637—1646. (in Chinese) 李红, 2011. 玉树地震房屋震害的几点启示. 青海师范大学学报（自然科学版）, 27（3）: 87—89Li H. , 2011. Several views from Yushu earthquake damages on the houses. Journal of Qinghai Normal University (Natural Science Edition), 27(3): 87—89. (in Chinese) 李华玥, 文鑫涛, 陈雅慧等, 2021.2020年国外地震灾害及其影响综述. 震灾防御技术, 16（3）: 583—588Li H. Y. , Wen X. T. , Chen Y. H. , et al. , 2021. Review of worldwide earthquake disasters and the impacts in 2020. Technology for Earthquake Disaster Prevention, 16(3): 583—588. (in Chinese) 李巧萍, 程新宇, 2020. 青海防震减灾综合能力快速提升. 防灾博览, （2）: 12—15. 李鑫, 殷翔, 姚生海等, 2021 a. 青海玛多7.4级地震重灾区房屋震灾调查及分析. 地震工程学报, 43（4）: 896—902Li X. , Yin X. , Yao S. H. , et al. , 2021 a. Investigation and analysis of earthquake disasters of houses in the stricken areas of Maduo M7.4 earthquake in Qinghai Province. China Earthquake Engineering Journal, 43(4): 896—902. (in Chinese) 李鑫, 姚生海, 殷翔等, 2021 b. 青海玛多7.4级地震极震区震灾调查及分析. 震灾防御技术, 16（3）: 429—436Li X. , Yao S. H. , Yin X. , et al. , 2021 b. Investigation and analysis of earthquake disaster in the extreme earthquake area of MS7.4 earthquake in Qinghai Province. Technology for Earthquake Disaster Prevention, 16(3): 429—436. (in Chinese) 李鑫, 李智敏, 盖海龙等, 2022. 青海门源MS6.9地震极震区震害调查与防灾建议. 震灾防御技术, 17（1）: 84—94Li X. , Li Z. M. , Gai H. L. , et al. , 2022. Investigation and prevention suggestion of earthquake disaster in the extreme earthquake area of MS6.9 earthquake in Menyuan County, Qinghai Province. Technology for Earthquake Disaster Prevention, 17(1): 84—94. (in Chinese) 李智敏, 盖海龙, 李鑫等, 2022.2022年青海门源MS6.9级地震发震构造和地表破裂初步调查. 地质学报, 96（1）: 330—335Li Z. M. , Gai H. L. , Li X. , et al. , 2022. Seismogenic fault and coseismic surface deformation of the Menyuan MS 6.9 earthquake in Qinghai, China. Acta Geologica Sinica, 96(1): 330—335. (in Chinese) 梁姗姗, 刘敬光, 邹立晔等, 2017.2016年10月17日青海杂多MS6.2地震震源机制解测定. 国际地震动态, （9）: 12—17Liang S. S. , Liu J. G. , Zou L. Y. , et al. , 2017. Determination of the focal mechanism solution of the October 17 2016 MS6.2 earthquake in Zadoi Qinghai.Recent Developments in World Seismology,(9):12—17.(in Chinese) 罗开海, 保海娥, 左琼, 2018. 我国建筑抗震设防水准的历史沿革、现状及展望. 地震工程与工程振动, 38（4）: 41—47Luo K. H. , Bao H. E. , Zuo Q. , 2018. The historical evolution, current situation and prospect of seismic precautionary level of buildings in China. Earthquake Engineering and Engineering Dynamics, 38(4): 41—47. (in Chinese) 潘家伟, 李海兵, Chevalier M. L. 等, 2022.2022年青海门源MS6.9地震地表破裂带及发震构造研究. 地质学报, 96（1）: 215—231Pan J. W. , Li H. B. , Chevalier M. L. , et al. , 2022. Coseismic surface rupture and seismogenic structure of the 2022 MS6.9 Menyuan earthquake, Qinghai Province, China. Acta Geologica Sinica, 96(1): 215—231. (in Chinese) 秦松涛, 李智敏, 谭明等, 2010. 青海玉树7.1级地震震害特点分析及启示. 灾害学, 25（3）: 65—70Qin S. T. , Li Z. M. , Tan M. , et al. , 2010. Analysis of damage characteristics of the M 7.1 Yushu earthquake of Qinghai and the enlightenments. Journal of Catastrophology, 25(3): 65—70. (in Chinese) 邵楠, 2011. 青海玉树地震建筑震害浅析. 工程抗震与加固改造, 33（4）: 136—140Shao N. , 2011. Analysis of building damages in Yushu earthquake of Qinghai. Earthquake Resistant Engineering and Retrofitting, 33(4): 136—140. (in Chinese) 孙鑫喆, 徐锡伟, 陈立春等, 2012.2010年玉树地震地表破裂带典型破裂样式及其构造意义. 地球物理学报, 55（1）: 155—170Sun X. Z. , Xu X. W. , Chen. L. C. , et al. , 2012. Surface rupture features of the 2010 Yushu earthquake and its tectonic implication. Chinese Journal of Geophysics, 55(1): 155—170. (in Chinese) 文鑫涛, 郑通彦, 2018.2016年中国大陆地震灾害损失述评. 灾害学, 33（3）: 141—144Wen X. T. , Zheng T. Y. , 2018. Review of earthquake damage losses in mainland China in 2016. Journal of Catastrophology, 33(3): 141—144. (in Chinese) 文鑫涛, 李华玥, 段乙好等, 2021.2020年中国大陆地震灾害损失述评. 震灾防御技术, 16（4）: 651—656Wen X. T. , Li H. Y. , Duan Y. H. , et al. , 2021. Earthquake disasters loss on Chinese mainland in 2020. Technology for Earthquake Disaster Prevention, 16(4): 651—656. (in Chinese) 姚生海, 盖海龙, 殷翔等, 2021. 青海玛多MS7.4地震地表破裂带的基本特征和典型现象. 地震地质, 43（5）: 1060—1072Yao S. H. , Gai H. L. , Yin X. , et al. , 2021. The basic characteristics and typical phenomena of the surface rupture zone of the Maduo MS7.4 earthquake in Qinghai. Seismology and Geology, 43(5): 1060—1072. (in Chinese) 殷翔, 李鑫, 马震等, 2021. 青海玛多MS7.4地震震害特点分析. 地震工程学报, 43（4）: 868—875Yin X. , Li X. , Ma Z. , et al. , 2021. Characteristics of seismic disasters caused by the Maduo MS 7.4 earthquake in Qinghai Province. China Earthquake Engineering Journal, 43(4): 868—875. (in Chinese) 詹艳, 梁明剑, 孙翔宇等, 2021.2021年5月22日青海玛多MS7.4地震深部环境及发震构造模式. 地球物理学报, 64（7）: 2232—2252Zhan Y. , Liang M. J. , Sun X. Y. , et al. , 2021. Deep structure and seismogenic pattern of the 2021.5. 22 Madoi (Qinghai) MS7.4 earthquake.Chinese Journal of Geophysics,64(7):2232—2252.(in Chinese) 张雪亭, 2006. 青海省大地构造格架研究. 北京: 中国地质大学（北京）.Zhang X. T. , 2006. Study on the tectonic framework of Qinghai. Beijing: China University of Geosciences Beijing). (in Chinese) -
下载:
点击查看大图
计量
- 文章访问数: 240
- HTML全文浏览量: 44
- PDF下载量: 32
- 被引次数: 0
