1. 绪论
蜂窝梁具有自重轻、节省钢材,方便管线穿越的同时还可以增加结构净高等优点 [1] ,近些年广泛应用。随着2001年911恐怖袭击,美国纽约世界贸易大厦发生倒塌后,钢结构的抗连续性倒塌性能引起工程界的广泛关注。
近些年多位学者对平面钢框架的抗连续性倒塌性能进行了研究,谢甫哲 [2] 通过对中柱失效下平面钢框架进行动力性能试验,研究结果表明框架梁在很短的时间内由弯曲为主的受力状态转变为受拉为主的悬链线受力状态。刘昕煜 [3] 提出一种新型钢–纤维复合筋材料,并以此设计一栋典型的三层RC框架结构,研究在不同拆柱工况下框架梁荷载系数、梁端轴力和弯矩随失效节点竖向位移的变化规律,并评估此新型结构的鲁棒性。孙玉婷 [4] 建立了普通钢框架结构和带有偏心支撑的钢框架结构的有限元模型,采用AP法对其进行静力和动力分析,研究结果表明边柱或中柱失效后,加入偏心支撑,底层其它柱的轴力基本没有明显变化,说明偏心支撑对结构的竖向承载力并没有显著提高。陆新征 [5] 等通过研究汶川地震中钢筋混凝土框架结构的震害机理,对钢筋混凝土框架结构的关键梁柱节点和框架柱进行了拟静力倒塌试验,结果表明节点区的破坏比相邻的框架梁和框架柱都要严重。众多研究结论表明:良好的延性能力可以有效提高结构的抗连续性倒塌能力,将蜂窝构件应用于钢结构中可能成为提高结构抗连续倒塌性能的有效方法,过往研究结果表明蜂窝构件具有良好的延性能力,但对其抗连续性倒塌动力性能研究较少。因此有必要对蜂窝梁平面钢框架倒塌能力进行研究。
本文基于拆除构件法,以端板连接蜂窝梁平面钢框架作为研究对象,通过有限元软件ABAQUS建立8个有限元模型对其进行有限元分析,研究中柱失效下不同端板厚度,开孔率,首孔距离及孔形等参数对其抗连续性倒塌动力性能的影响,并提出适用于端板连接蜂窝梁平面钢框架的动力放大系数。
2. 有限元模型建立与验证
2.1. 有限元模型建立
为研究端板连接蜂窝梁平面钢框架的抗连续性倒塌性能研究,如图1所示,以中柱失效下端板连接蜂窝梁平面钢框架为研究对象,研究结构的抗连续性倒塌动力性能。
Figure 1. Schematic diagram of failed column
图1. 失效柱示意图
2.1.1. 模型尺寸
试件名称:DZ60-14-250。试件尺寸:每层层高2500 mm,跨度均为4200 mm,梁、柱采用H型钢,梁截面尺寸为250 × 125 × 6 × 9 (mm),柱截面尺寸为300 × 200 × 8 × 12 (mm),螺栓为10.9级M20高强螺栓,梁与柱的连接采用高强螺栓连接,端板尺寸为420 × 200 × 14 (mm)。
2.1.2. 动力效应模拟方法
拆除构件法也叫做AP法 [6] ,是指在人为的选择在结构中一个或者多个构件损坏后,分析剩余结构在自重和荷载作用下的内力重分布,以评估结构的抗连续性倒塌能力。一般做法为:对柱失效前的原模型进行静力线性分析,得到柱失效前的内力。去除失效的柱并将该柱的内力反作用于失效点上,使剩余结构与原结构内力等效。本文通过准静态分析得到与原结构内力等效的初始状态,将反向作用的内力值在一定的时间内减小到0,计算得到柱失效后剩余结构的动力效应。
2.1.3. 网格划分及本构关系
蜂窝梁平面钢框架结构的构件均采用C3D8R实体单元,其中包括梁、柱、端板和螺栓构件。实体单元具有精度高,运算准确等优点 [7] ,为保证具有良好的计算精度与合适的时间长度,在网格划分中柱的网格尺寸为50 mm,梁的网格尺寸为50 mm,端板的网格尺寸为40 mm,螺栓的网格尺寸为8 mm,网格划分如下图2所示。
Figure 2. Mesh division of structural components
图2. 结构构件网格划分
模型中钢材本构关系采用双折线模型 [8] ,弹性模量为2.2 × 105 MPa,泊松比为0.3,屈服强度为370 MPa。
2.1.4. 接触及相互作用
在试验过程中采用与实际情况相符的建模过程,有限元模拟中设置三种接触,柱与端板接触,端板与蜂窝梁接触,螺栓与端板及柱接触,三种接触滑移公式选择小滑移,相互作用中切向作用采用“罚”,摩擦系数为0.3,法向作用定义为“硬接触”。
2.1.5. 荷载及边界效应
在端板连接蜂窝梁平面钢框架分析中,为保证非失效柱柱底稳定,将端板连接蜂窝梁平面钢框架结构中的钢柱柱底按X、Y、Z三个方向的位置与转角全部加以约束来模拟真实试验中的柱脚刚接。
进行试验时,如图3所示,在ABAQUS中施加命令,固定除中柱外其它四根柱子,向螺栓施加165,000 N螺栓荷载。根据DoD规范 [9] ,倒塌荷载按照(1.2 × DL + 0.5 × LL)的方法取值,其中的DL和LL分别代表恒载和活载,计算其在蜂窝梁上施加5.8 kN/m2。
2.2. 有限元模型验证
采用文献 [10] 中的平面钢框架结构中柱失效作用下抗倒塌动力试验验证有限元模拟方法的准确性。试件为一榀三层四跨平面钢框架结构,跨度为500 mm,总跨度2000 mm,首层层高510 mm,二三层层高400 mm,梁柱均为方钢管,梁的截面尺寸为20 * 1.2,柱的截面尺寸为30 * 2。破坏模式如图4。
如图5所示,失效柱极限的下降位移为97 mm,模拟结果得到的失效柱极限下降位移为95 mm,可见有限元模型计算结果与试验结果的位移时程曲线吻合良好,误差不到5%。失效柱旁框架梁形成塑性铰后破坏,与模拟中的形成塑性铰破坏模式吻合,表明该有限元建模具有良好的精度。
Figure 5. Comparison of vertical displacement time history curves of failure point
图5. 失效点竖向位移时程曲线对比
3. 端板连接蜂窝梁平面钢框架连续倒塌动力性能分析
破坏模式
破坏模式如图6所示,0.7 s时突然撤掉底层中柱,失效点迅速下落,失效柱两侧梁发生较大变形,随着失效点的竖向位移增加,失效柱两旁蜂窝梁首孔处逐渐形成塑性铰。在1.65 s时失效柱达到极限竖向位移,此时失效柱底的反作用力为0。从2 s开始一直到结束结构振动回弹,幅度逐渐减缓,并很快趋于稳定,此时结构内部传力路径发生改变,蜂窝梁由受弯状态转为受拉状态,使结构重新恢复稳定。
失效柱的位移时程曲线如图7所示,在试件底层中柱失效后,失效柱发生向下位移,在0.95 s后达到试件的极限竖向位移ΔVmax = 419.1 mm,并回弹到333 mm,此后振动衰减,结构的振动随时间最终稳定在370 mm附近。DOD2005 [11] 对于梁端竖向转角限值的规定超过6˚即可认定水平构件发生损坏,结构发生连续倒塌。试件DZ60-14-250中蜂窝梁的竖向极限转角为5.72˚,并未超过6˚,即试件DZ60-14-250没有发生连续倒塌。
Figure 7. Vertical displacement time history curve of DZ60-14-250 failed column
图7. DZ60-14-250失效柱竖向位移时程曲线
图8~10分别为试件DZ60-14-250中首层、第二层、顶层框架柱的水平位移。试件DZ60-14-250首层框架柱的极限水平位移Vmax = 8.4 mm,第二层框架柱的极限水平位移Vmax = 17.3 mm,顶层框架柱的极限水平位移Vmax = 25.7 mm,在振动衰减过后,最终稳定在5.8 mm,13.0 mm,20.0 mm。顶层框架柱的水平极限转角为0.38˚。按照DOD2005,对初级和中高级安全等级建筑中柱的极限允许塑形转角为0.0698 rad (4˚),很显然框架柱的极限转角远远小于4˚,因此端板和蜂窝梁协同工作可以有效的防止框架柱发生过大的水平位移。
Figure 8. Horizontal displacement of the first floor of DZ60-14-250
图8. DZ60-14-250首层水平位移
Figure 9. Horizontal displacement of the second layer of DZ60-14-250
图9. DZ60-14-250第二层水平位移
Figure 10. DZ60-14-250 top horizontal displacement
图10. DZ60-14-250顶层水平位移
4. 参数化分析
本章主要研究在不同的开孔率,首孔距离,端板厚度和开孔形式等参数下对于端板连接蜂窝梁平面钢框架的影响,分析蜂窝梁平面钢框架的抗连续性倒塌动力性能,各试件参数见表1。
4.1. 端板厚度
从图11中可以得到试件的极限竖向位移随端板厚度的增大而减小。试件DZ60-6-250的极限竖向位移Vmax = 538.1 mm,蜂窝梁竖向极限竖向转角为7.34˚,超过6˚,因此判定结构发生连续性倒塌。试件DZ60-10-250的极限竖向位移Vmax = 437.2 mm,DZ60-14-250的极限竖向位移Vmax = 419.7 mm,蜂窝梁竖向极限转角分别为5.96˚、5.72˚,并未超过6˚,结构并未发生连续倒塌。
Figure 11. Vertical displacement time history curves of failed columns with different end plate thicknesses
图11. 不同端板厚度的失效柱竖向位移时程曲线
4.2. 开孔率
从图12可以得到试件的极限竖向位移随开孔率的增大而增大,试件DZ70-14-250极限竖向位移Vmax = 425.7 mm,试件DZ60-14-250极限竖向位移Vmax = 419.1 mm,试件DZ50-14-250极限竖向位移Vmax = 401.6 mm,蜂窝梁在竖直方向的极限转角分别为5.47˚、5.72˚、5.86˚。试件DZ70-14-250的蜂窝梁极限竖向转角临近6˚,结构并未发生连续性倒塌。由此可以得到随着开孔率增大,结构的抗倒塌性能变差,但总体相差不大,试件DZ60-14-250与试件DZ70-14-250在振动结束趋于稳定时基本持平,在开孔率50%~70%时结构的整体抗倒塌情况基本相同,对于蜂窝梁平面钢框架结构而言,并不是开孔率越大越好。因此建议开孔率在50%~70%。
Figure 12. Vertical displacement time history curves of failed columns with different porosity ratios
图12. 不同开孔率的失效柱竖向位移时程曲线
4.3. 首孔距离
从图13可以得到试件的极限竖向位移随首孔距离的增大而减小。试件DZ60-14-375极限竖向位移Vmax = 425.0 mm,蜂窝梁竖向极限竖向转角为5.8˚。试件DZ60-14-125极限竖向位移Vmax = 394.1 mm,蜂窝梁竖向极限竖向转角为5.37˚,试件DZ60-14-125、DZ60-14-250和试件DZ60-14-375未发生连续倒塌。试件DZ60-14-125与试件DZ60-14-250极限位移相差不大,当首孔距离大于250 mm (1.0 h)时,失效点处竖向位移增大明显。因此建议首孔距离在0.5~1.5 h。
Figure 13. Vertical displacement time history curves of failed columns with different first hole distances
图13. 不同首孔距离的失效柱竖向位移时程曲线
4.4. 孔型
从图14可以得到试件DZY60-14-250的极限竖向位移大于试件DZ60-14-250。试件DZY60-14-375
Figure 14. Vertical displacement time history curves of failed columns with different hole types
图14. 不同孔型的失效柱竖向位移时程曲线
极限竖向位移Vmax = 441.8 mm,蜂窝梁竖向极限竖向转角为6.03˚,即判定结构发生连续性倒塌。所以正六边形孔端板连接蜂窝梁平面钢框架对于结构抗连续倒塌的性能更好一点。
5. 动力效应放大系数
结构倒塌是一个复杂的动力过程,在进行分析结构倒塌事故时通常有两种方法:一是进行动力分析法。此方法运算准确但是耗费时间较多。二是静力分析法。此方法需要将动力放大系数(λ)将荷载进行放大然后再用静力方法分析来模拟结构的动力反应。GSA2003 [12] 提出了在静力分析中将动力放大系数定在2.0来近似考虑倒塌过程中的动力效应。Ruth [13] 等对某一框架结构进行动力分析和静力分析,得出2.0的动力放大系数偏保守,因此有必要对动力效应放大系数进行研究。
动力放大系数可分为基于位移和基于荷载两种方法进行计算得到。其中位移法的做法为在相同荷载(P)的作用下对应动力分析法的倒塌结构位移和静力分析法得到的位移的比值。荷载法为相同的位移(U)对应下静力分析法得到的荷载与动力分析法得到的荷载的比值。本文采用位移法计算动力效应放大系数。
以DZ60-14-250为例探究在中柱失效不同情况下动力放大系数。动力放大系数(λ)取不同的值时,静力分析法产生的位移在某一值时与动力分析法产生的位移相等,则此时的动力放大系数(λ)即为最合理的动力放大系数。分别提取失效柱底端的位移,其中静力分析法失效柱底端的位移为,其中静力分析法失效柱底端的位移为
,即通过比较
两者位移比与λ的关系,从而确定动力放大系数。当
时,说明静力荷载在放大λ倍后的位移与动力荷载的位移相等。
Figure 15. Vertical displacement ratio curve of the failed column top when the middle column fails
图15. 中柱失效时失效柱柱顶竖向位移比值曲线
图15表明同一结构在中柱失效下得到的动力放大系数为1.17。所以2.0的动力放大系数偏保守,建议端板连接蜂窝梁平面钢框架结构的动力放大系数取1.2。
6. 结论
(1) 端板连接蜂窝梁平面钢框架在中柱失效作用下,失效柱两侧蜂窝梁随位移的增加形成塑性铰,试件的极限竖向位移ΔVmax = 419.1 mm,蜂窝梁的竖向极限转角为5.72˚,可有效防止结构发生连续倒塌。
(2) 本文对端板连接蜂窝梁平面钢框架进行了参数分析,通过分析结果给出了一些建议:端板厚度建议设置在10~14 mm,开孔率建议设置在50%~70%,首孔距离建议设置在0.5~1.5 h,正六边形孔端板连接蜂窝梁平面钢框架拥有更好的抗连续性倒塌动力性能。
(3) 基于位移法计算出试件DZ60-14-250的动力放大系数为1.17,规范GSA2003中建议的2.0相比于本文计算结果更为保守,因此建议采用位移法计算端板连接蜂窝梁平面钢框架的动力放大系数。