高烈度区复杂高层钢框架中心支撑结构性能化设计研究
0 引言
强震一直威胁着生命财产安全, 给人类社会带来了巨大灾难。随着经济高速发展和城市化进程的加速, 结构工程朝着体量更大、高度更高、体系更加复杂的方向发展, 随之而来的是一旦强震发生, 将有可能带来更为重大的社会负面影响和经济损失
目前, 部分振动控制装置已形成不同类型的产品
1 抗震性能评估
1.1 工程概况
某复杂高层建筑设计使用年限为50年, 结构安全等级为二级 (结构重要性系数γ0=1.0) , 丙类抗震设防, 抗震设防烈度为8度 (0.30g) , 设计地震分组为第二组, 场地类别为Ⅲ类, 特征周期为0.55s。值得注意的是, 本工程场地为近断层, 近场影响系数取为1.25。该建筑采用钢框架中心支撑结构体系, 平面呈现椭圆形, 地下4层, 地上27层, 建筑物地上结构总高95m。全楼钢材采用Q345B。主要框架柱截面为□600×600×40×40, 主要框架梁截面为HW600×300×12×20。建筑结构模型与平面图如图1所示。
1.2 精细化有限元模型的建立与地震动的选取
1.2.1 精细化有限元模型的建立
能否充分认识结构在地震下的损伤特性是结构性能化设计的重要基础。大量研究表明数值模拟已经成为了解各类结构地震损伤机理的重要手段, 其中大型商业通用有限元程序PERFORM-3D已被广泛应用于高层和超高层结构的抗震性能评估和指导上述结构的性能化设计
1.2.2 动力特性分析
合理可靠的精细有限元模型是研究结构抗震性能的重要基础, 本文通过对比设计软件PKPM和弹塑性分析软件PERFORM-3D的模态分析和反应谱分析的典型结果验证该模型在弹性阶段的可靠性。该结构模态分析获得的前6阶周期对比如表1所示, 从中可以看出两种软件计算所得的各阶周期整体吻合良好, 相对误差均不超过5%。两种软件反应谱分析计算所得的Y向的层间位移角分布对比如图4所示, 从中可以看出两者层间位移角分布基本一致, 相对误差不超过5%。结果表明本文基于PERFORM-3D建立的精细化有限元模型在弹性阶段具有较高的可靠性, 可用于本结构的抗震性能评估。
1.2.3 地震动的选取
按照《建筑抗震设计规范》 (GB 50011—2010) (2016年版)
PKPM与PERFORM-3D结构基本动力特性对比表1
阶数 |
1 | 2 | 3 | 4 | 5 | 6 | |
周期 /s |
PKPM |
2.67 | 2.47 | 2.31 | 0.78 | 0.72 | 0.70 |
PERFORM-3D |
2.65 | 2.54 | 2.33 | 0.79 | 0.74 | 0.71 | |
相对误差/% |
3.23 | 3.20 | 1.17 | 1.54 | 2.23 | 0.57 | |
振型 |
Y向 平动 |
X向 平动 |
扭转 |
Y向 平动 |
X向 平动 |
扭转 |
8度 (0.3g) 多遇地震作用下基底剪力与反应谱法对比结果表2
项目 |
X向 |
Y向 | |||
基底剪 力/kN |
相对误 差/% |
基底剪 力/kN |
相对误 差/% |
||
PKPM |
反应谱法 | 19 375 | — | 20 521 | — |
PERFORM -3D |
反应谱法 |
20 007 | — | 21 382 | — |
RSN12波 |
17 802 | -11.0 | 21 525 | 0.7 | |
RSN86波 |
19 867 | -0.7 | 22 364 | 4.6 | |
RGB波 |
17 682 | -11.6 | 19 778 | -7.5 | |
平均值 |
18 450 | -7.8 | 21 222 | -0.7 |
1.3 结构性能化评估
1.3.1 结构在地震动作用下的响应
采用上述地震波, 通过调整加速度峰值分别对结构进行多遇地震与罕遇地震作用下的动力时程分析, 即将地震动的峰值加速度调幅至375cm/s2和637.5cm/s2。多遇地震作用下X向最大层间位移角为1/324, 满足抗震规范1/250的限值要求;Y向最大层间位移角为1/249, 基本可视为满足规范1/250的限值要求。结构在罕遇地震作用下的层间位移角如图6所示。从图中可以看出, 在罕遇地震下结构X向最大层间位移角为1/80, Y向为1/80, 均满足规范1/50的限值要求。
1.3.2 结构构件塑性铰发展情况
本工程的结构分析采用基于性能的抗震设计方法。参考ASCE41-13
从表3中可以看出, 在3条地震波作用下, 结构6 258个梁构件中平均有13.02%达到IO性能等级, 意味着梁构件开始屈服, 这对于结构而言是可接受的;但值得注意的是, 罕遇地震下存在梁构件达到LS性能等级, 无法满足业主对于该高档商用办公楼的建筑功能需求;同时结构1 322个柱构件中的5.96%也达到了IO性能等级。由于本工程建筑功能较为重要, 经业主要求和专家委员建议采用金属剪切型阻尼器提高结构抗震性能。
各地震波下构件塑性铰个数表3
地震波 |
总塑性铰 |
梁塑性铰 | 柱塑性铰 | |||||||
IO |
LS | CP | IO | LS | CP | IO | LS | CP | ||
RSN12波 |
X向 |
1 118 | 0 | 0 | 1 042 | 0 | 0 | 76 | 0 | 0 |
Y向 |
672 | 11 | 0 | 603 | 11 | 0 | 69 | 0 | 0 | |
RSN86波 |
X向 |
1 057 | 0 | 0 | 969 | 0 | 0 | 88 | 0 | 0 |
Y向 |
779 | 6 | 0 | 682 | 6 | 0 | 97 | 0 | 0 | |
RGB波 |
X向 |
1 021 | 0 | 0 | 959 | 0 | 0 | 62 | 0 | 0 |
Y向 |
716 | 9 | 0 | 635 | 0 | 0 | 81 | 0 | 0 |
2 结构减震方案及其性能化评估
2.1 阻尼器方案
本工程依据抗震体系性能化评估的结果, 以减小结构层间位移角与减少结构构件损伤为目标, 考虑在1~23层适当位置沿结构两个主轴方向分别设置金属剪切型阻尼器, 平面布置图如图7所示。阻尼器性能参数如表4所示, 其恢复力模型如图8所示。采用PERFORM-3D中的剪切铰模型模拟金属剪切型阻尼器。
金属剪切型阻尼器参数及布置数量表4
屈服力 | 屈服位移 | 屈服后刚度比 | 极限位移 | 布置位置 | 布置总数量 |
500kN |
1mm | 0.01 | 60mm | 1~23层 | 92个 |
2.2 性能化评估
对减震结构进行大震弹塑性分析, 评估抗震性能, 包括层间位移角、能量分配 (尤其是阻尼器滞回耗能) 和构件塑性铰展开情况。
2.2.1 层间位移角
分析结果中RSN12波作用下结构响应最大, 减震前后结构的层间位移角对比如图9所示, 可以看出X向与Y向最大层间位移角分别减小29.8%与10.6%。同样, RSN86波作用下, 结构减震前后X向与Y向最大层间位移角分别减小5.8%与3.2%, RGB波作用下, 结构减震前后X向与Y向最大层间位移角分别减小18.9%与3.7%。
2.2.2 能量耗散情况
从分析结果可以看出, 减震前后结构总塑性耗能 (包括阻尼器耗能和构件塑性耗能) 、瑞雷阻尼耗能、应变能和动能占比基本保持不变。减震后结构在各条波下各部分能量与地震总输入能量的比值如表5所示。RSN12波作用下, X向和Y向阻尼器耗能分别占总塑性耗能的36.2%和21.1%;RSN86波作用下, X向和Y向阻尼器耗能分别占总塑性耗能的32.4%和20.1%;RGB波作用下, X向和Y向阻尼器耗能分别占总塑性耗能的36.15%和18.72%。由此可见, 阻尼器的耗能作用较明显, 可有效减少结构主要构件的塑性耗能及其损伤程度。
各地震波下各部分能量各部分占总能量之比表5
地震波 |
总能量 /kJ |
各部分占总能量之比/% |
|||||
阻尼器 耗能 |
构件塑 性耗能 |
瑞雷阻 尼能量 |
应变能 | 动能 | |||
RSN12波 |
X向 |
75 614 | 11.25 | 18.87 | 66.69 | 2.64 | 0.54 |
Y向 |
70 911 | 4.70 | 17.55 | 75.00 | 2.42 | 0.33 | |
RSN86波 |
X向 |
77 999 | 9.78 | 20.38 | 66.96 | 1.95 | 0.94 |
Y向 |
85 401 | 5.23 | 20.83 | 70.89 | 1.68 | 1.37 | |
RGB波 |
X向 |
48 753 | 10.89 | 19.23 | 58.61 | 3.12 | 8.15 |
Y向 |
57 362 | 5.64 | 24.49 | 58.82 | 2.40 | 8.64 |
2.2.3 塑性铰开展情况
减震结构的塑性铰开展情况如表6所示。从表中可以得出, 各条地震波下梁构件均未达到LS性能等级, 满足性能目标;达到IO性能等级的柱构件数量大幅度减少, 平均降幅为29.8%, 最大降幅达43.6%。由此可看到, 金属剪切型阻尼器的设置大大减轻地震作用下梁、柱构件的损伤程度, 提升了结构的抗震性能。
在原有结构尺寸不变情况下, 其原结构自振周期与增设金属剪切型阻尼器后的结构自振周期如表7所示。由表中可知, 增设金属剪切型阻尼器后, 因阻尼器自身具有初始刚度, 结构自振周期有所减小。
减震后各地震波作用下构件塑性铰个数表6
地震波 |
总塑性铰 |
梁塑性铰 | 柱塑性铰 | |||||||
IO |
LS | CP | IO | LS | CP | IO | LS | CP | ||
RSN12波 |
X向 |
1 019 | 0 | 0 | 968 | 0 | 0 | 51 | 0 | 0 |
Y向 |
645 | 0 | 0 | 592 | 0 | 0 | 53 | 0 | 0 | |
RSN86波 |
X向 |
1 076 | 0 | 0 | 1014 | 0 | 0 | 62 | 0 | 0 |
Y向 |
739 | 0 | 0 | 669 | 0 | 0 | 70 | 0 | 0 | |
RGB波 |
X向 |
856 | 0 | 0 | 821 | 0 | 0 | 35 | 0 | 0 |
Y向 |
683 | 0 | 0 | 620 | 0 | 0 | 63 | 0 | 0 |
减震前后结构基本动力特性表7
阶数 |
1 | 2 | 3 | 4 | 5 | 6 | |
自振周 期/s |
减震前 |
2.65 | 2.54 | 2.33 | 0.79 | 0.74 | 0.71 |
减震后 |
2.59 | 2.45 | 2.28 | 0.79 | 0.72 | 0.69 | |
振型 |
Y向 平动 |
X向 平动 |
扭转 |
Y向 平动 |
X向 平动 |
扭转 |
3 优化布设
3.1 阻尼器屈服力参数
阻尼器的关键设计参数 (如屈服力等) 及其沿结构高度方向的布置方案对结构的抗震性能存在较大影响, 因此本文通过设置不同的阻尼器关键设计参数和布置方案对其展开研究。
通过对阻尼器参数的调整, 来对结构进行优化设计, 选取如表8所示的5种初始刚度相同、屈服力不同的阻尼器, 对其进行大震弹塑性分析, 以RSN12波为例, 各部分能量所占地震总输入能量比值的对比如表9所示, 构件性能水平的对比如表10所示。
阻尼器参数表8
名称 |
屈服位移/mm | 初始刚度/ (kN/m) | 屈服力/kN |
Y400-0.8 |
0.8 | 500 | 400 |
Y500-1.0 |
1.0 | 500 | 500 |
Y600-1.2 |
1.2 | 500 | 600 |
Y800-1.6 |
1.6 | 500 | 800 |
Y1000-2.0 |
2.0 | 500 | 1 000 |
RSN12波下不同阻尼器各部分能量所占总能量之比表9
方向 |
阻尼器 名称 |
总能量 /kJ |
各部分占总能量之比/% |
||||
阻尼器 耗能 |
构件塑 性耗能 |
瑞雷阻 尼能量 |
应变能 | 动能 | |||
X向 |
Y400-0.8 |
75 528 | 12.23 | 18.82 | 65.94 | 2.53 | 0.48 |
Y500-1.0 |
75 614 | 11.25 | 18.87 | 66.69 | 2.64 | 0.54 | |
Y600-1.2 |
75 231 | 9.95 | 19.00 | 67.73 | 2.74 | 0.59 | |
Y800-1.6 |
75 072 | 6.83 | 20.36 | 69.34 | 2.83 | 0.63 | |
Y1000-2.0 |
75 373 | 4.46 | 21.28 | 70.75 | 2.87 | 0.65 | |
Y向 |
Y400-0.8 |
66 407 | 6.03 | 16.49 | 74.63 | 2.51 | 0.33 |
Y500-1.0 |
70 911 | 4.70 | 17.55 | 75.00 | 2.42 | 0.33 | |
Y600-1.2 |
67 737 | 3.64 | 17.50 | 75.95 | 2.54 | 0.36 | |
Y800-1.6 |
68 536 | 1.95 | 18.22 | 76.91 | 2.55 | 0.37 | |
Y1000-2.0 |
69 326 | 1.10 | 18.81 | 77.18 | 2.54 | 0.37 |
RSN12波下不同阻尼器塑性铰出铰情况表10
方向 |
阻尼器 名称 |
总塑性铰 |
梁塑性铰 | 柱塑性铰 | ||||||
IO |
LS | CP | IO | LS | CP | IO | LS | CP | ||
X向 |
Y400-0.8 |
1 025 | 0 | 0 | 974 | 0 | 0 | 51 | 0 | 0 |
Y500-1.0 |
1 019 | 0 | 0 | 968 | 0 | 0 | 51 | 0 | 0 | |
Y600-1.2 |
1 024 | 0 | 0 | 971 | 0 | 0 | 53 | 0 | 0 | |
Y800-1.6 |
1 028 | 0 | 0 | 975 | 0 | 0 | 53 | 0 | 0 | |
Y1000-2.0 |
1 034 | 0 | 0 | 981 | 0 | 0 | 53 | 0 | 0 | |
Y向 |
Y400-0.8 |
639 | 0 | 0 | 587 | 0 | 0 | 52 | 0 | 0 |
Y500-1.0 |
645 | 0 | 0 | 592 | 0 | 0 | 53 | 0 | 0 | |
Y600-1.2 |
651 | 0 | 0 | 598 | 0 | 0 | 53 | 0 | 0 | |
Y800-1.6 |
662 | 0 | 0 | 610 | 0 | 0 | 52 | 0 | 0 | |
Y1000-2.0 |
658 | 0 | 0 | 605 | 0 | 0 | 53 | 0 | 0 |
从表9与表10中可知, 随着阻尼器屈服力的增加, 结构总塑性耗能 (包括阻尼器耗能和构件塑性耗能) 有小幅度下降, 其中阻尼器耗能所占的比例有较大程度降低, 导致构件塑性耗能所占比例有所增加, 达到IO性能等级的构件数量有所增加。同时阻尼器屈服力的增加, 层间位移角也会有所增大, 最大相差可达3.6%。综上说明, 屈服力较小的阻尼器可消耗更多地震能量, 有效抑制塑性铰展开, 使层间位移角有所减小。
3.2 阻尼器布置位置
虽然上文提出的阻尼器初步布置设计方案已取得较好的减震效果, 但阻尼器在每层均布置, 数量较多, 经济性投入相对较大, 因此有必要进一步优化该方案, 以取得抗震性能和经济性的较好平衡。为此, 本文提出了如图10所示的4种新阻尼器布置方案。在更改阻尼器布置方案时, 其基本原则为地震作用下层间位移角相对较小的楼层去除阻尼器。其中减震方案1将1层、2层和3层 (共计3层) 阻尼器去除, 共减少阻尼器12个 (减少13.0%) ;减震方案2将1~5层 (共计5层) 的阻尼器去除, 共减少阻尼器20个 (减少21.7%) ;方案3将1层、3层、5层、16层、18层、20层和22层 (共计6层) 阻尼器去除, 值得注意的是, 由于大量减少阻尼器会使得结构刚度突变, 故采用阻尼器隔层布置的方案, 共减少阻尼器24个 (减少26.1%) ;减震方案4与减震方案3类似, 也采用阻尼器隔层布置的方案, 与之不同的是将层间位移角最小的1~5层阻尼器全部去除, 即将1~5层、16层、18层、20层和22层 (共计9层) 去除, 共减少阻尼器36个 (减少39.1%) 。
不同优化方案在RSN12波下的层间位移角最大值与原减震方案之比如图11所示, 能量耗散情况如表11所示, 塑性铰开展情况如表12所示。
从图11和表11, 12可以得出, 虽然随着阻尼器的减少, 阻尼器塑性耗能有所减少, 但层间位移角变化幅度较小, 最大增幅仍小于5%。值得注意的是, 减震方案3采用了将层间位移角较小的楼层进行隔层布置阻尼器的方法, 所以阻尼器塑性耗能充分, 塑性铰数量较其他减震方案有所减少, 最大可减少1.32%。
综上说明阻尼器优化布设可以使减震效果基本维持不变, 并使经济性得到大幅度提高。同时在层间位移角变化平缓且数值较小处间断布置阻尼器的方法可以为工程提供参考。
RSN12波下不同减震布置时各部分能量所占总能量之比表11
地震波 |
减震 方案 |
总能量 /kJ |
各部分占总能量之比/% |
||||
阻尼器 耗能 |
构件塑 性耗能 |
瑞雷阻 尼能量 |
应变能 | 动能 | |||
X向 |
原减震方案 |
79 746 | 10.89 | 19.16 | 66.92 | 2.52 | 0.51 |
减震方案1 |
80 213 | 7.68 | 23.20 | 66.37 | 2.43 | 0.32 | |
减震方案2 |
80 534 | 6.77 | 24.23 | 66.33 | 2.34 | 0.32 | |
减震方案3 |
79 762 | 8.97 | 22.53 | 65.79 | 2.40 | 0.31 | |
减震方案4 |
80 602 | 7.42 | 24.33 | 65.63 | 2.30 | 0.33 | |
Y向 |
原减震方案 |
70 911 | 4.70 | 17.55 | 75.00 | 2.42 | 0.33 |
减震方案1 |
69 255 | 2.25 | 19.66 | 75.02 | 2.91 | 0.16 | |
减震方案2 |
68 376 | 1.86 | 19.44 | 75.61 | 2.90 | 0.19 | |
减震方案3 |
69 205 | 3.21 | 19.05 | 74.68 | 2.89 | 0.17 | |
减震方案4 |
68 884 | 2.66 | 19.60 | 74.67 | 2.87 | 0.21 |
RSN12波下不同减震布置时塑性铰出铰情况表12
地震 波 |
减震 方案 |
总塑性铰 |
梁塑性铰 | 柱塑性铰 | ||||||
IO |
LS | CP | IO | LS | CP | IO | LS | CP | ||
X 向 |
原减震方案 |
1 019 | 0 | 0 | 968 | 0 | 0 | 51 | 0 | 0 |
减震方案1 |
1 059 | 0 | 0 | 989 | 0 | 0 | 70 | 0 | 0 | |
减震方案2 |
1 063 | 0 | 0 | 992 | 0 | 0 | 71 | 0 | 0 | |
减震方案3 |
1 049 | 0 | 0 | 979 | 0 | 0 | 70 | 0 | 0 | |
减震方案4 |
1 053 | 0 | 0 | 982 | 0 | 0 | 71 | 0 | 0 | |
Y 向 |
原减震方案 |
645 | 0 | 0 | 592 | 0 | 0 | 53 | 0 | 0 |
减震方案1 |
664 | 0 | 0 | 597 | 0 | 0 | 67 | 0 | 0 | |
减震方案2 |
655 | 0 | 0 | 588 | 0 | 0 | 67 | 0 | 0 | |
减震方案3 |
658 | 0 | 0 | 588 | 0 | 0 | 70 | 0 | 0 | |
减震方案4 |
661 | 0 | 0 | 591 | 0 | 0 | 70 | 0 | 0 |
故建议采用屈服力为400kN的阻尼器, 布置形式可采用减震方案3。采用减震方案3的结构在RSN12波下的层间位移角与原结构层间位移角对比如图12所示, 可以看出X向与Y向最大层间位移角分别减小25.2%与8.3%, 减震效果较为明显。
4 结论
本文采用金属剪切型阻尼器, 对一位于高烈度区的复杂钢框架中心支撑高层结构进行了性能化设计研究, 主要得出以下结论:
(1) 对传统抗震设计方案的性能化评估结果表明, 尽管其虽可满足规范相关要求, 但结构出现大量达到IO水准等级的构件和一定数量达到LS水准等级的构件, 抗震性能相对较差, 无法满足其较高的使用性能需求。
(2) 对采用金属剪切型阻尼器的减震设计方案与抗震方案进行对比分析结果表明, 在罕遇地震作用下, 结构层间位移角最大减幅可达29.8%, 阻尼器耗能状况良好, 所有构件均未达到LS水准等级, 达到IO水准等级的构件也有一定程度的减少, 结构性能有较大提升。
(3) 在减震初步设计方案基础上, 对阻尼器参数与阻尼器布置位置与数量进行了优化布设和分析比较, 结果表明:1) 阻尼器初始刚度相同时, 屈服力越大, 减震效果将变差, 结构出现塑性铰数量越多, 层间位移角也有所增大, 最大增幅为3.6%;2) 阻尼器宜尽量布置在结构层间位移角较大的位置, 其他位置减少阻尼器数量, 可在基本满足罕遇地震减震效果的前提下, 大大提高减震方案的经济性。
[2] 黄镇, 李爱群. 建筑结构金属消能器减震设计[M]. 北京:中国建筑工业出版社, 2015.
[3] 邱敏, 张学文, 陆中玏, 等. 土木工程结构振动控制的研究现状与展望[J]. 安全与环境工程, 2013, 20 (3) :14-18.
[4] 建筑消能阻尼器:JG/T 209—2007[S]. 北京:中国标准出版社, 2007.
[5] 建筑隔震橡胶支座:JG 118—2000[S]. 北京:中国标准出版社, 2000.
[6] SYMANS M D, CHARNEY F A, WHITTAKER A S, et al. Energy dissipation systems for seismic applications: current practice and recent developments[J]. Journal of Structural Engineering, ASCE, 2008, 134 (1) :3-21.
[7] 徐自然, 苏骏, 崔家春, 等. 上海浦东国际机场三期扩建工程卫星厅消能减震项目抗震性能化分析[J]. 建筑结构, 2017, 47 (12) :23-28.
[8] 陈巍, 范重, 张娜, 等. 厦门新机场结构耗能减震设计研究[J]. 建筑结构, 2017, 47 (8) :7-14.
[9] 于森林, 张盟, 李博宇, 等. 北京鲜活农产品项目消能减震设计[J]. 建筑结构, 2017, 47 (8) :64-68, 86.
[10] 张东方, 白羽.某装设金属阻尼器框架结构的减震控制分析[J].建筑结构, 2017, 47 (S2) :289-293.
[11] 陈学伟, 韩小雷, 林生逸, 等. 中洲中心二期结构抗震性能分析[J]. 建筑结构学报, 2010, 31 (1) :101-109.
[12] 韩小雷, 陈学伟, 林生逸, 等. 基于纤维模型的超高层钢筋混凝土结构弹塑性时程分析[J]. 建筑结构, 2010, 40 (2) :13-16.
[13] 郝效强, 李楚舒, 李立. PERFORM-3D在结构抗震性能分析与评估中的应用[J]. 建筑结构, 2012, 42 (S2) :202-206.
[14] 齐麟, 孙荣来, 黄兆纬, 等.基于性能的天津中信城市广场7号楼超限高层结构抗震分析[J].建筑结构, 2017, 47 (24) :26-30.
[15] 张沛洲, 欧进萍. 基于纤维模型的钢混框架结构拟静力试验数值模拟[J]. 建筑结构, 2013, 43 (18) :64-69.
[16] 韩小雷, 陈学伟, 林生逸. 基于宏观单元的结构非线性分析方法、算例及工程应用[J].工程力学, 2010, 27 (S1) :59-67.
[17] 杜轲, 孙景江, 许卫晓. 纤维模型中单元、截面及纤维划分问题研究[J]. 地震工程与工程振动, 2012, 32 (5) :39-46.
[18] 建筑抗震设计规范:GB 50011—2010[S]. 2016年版. 北京:中国建筑工业出版社, 2016.
[19] Seismic evaluation and retroft of existing buildings:ASCE 41-13[S]. Reston:American Society of Civil Engineers, 2013.