吴应雄, 王兆樑, 祁 皑, 颜学渊
(福州大学 土木工程学院,福州 350108)
隔震结构是通过设置水平刚度较小的隔震层来减少向上部结构传递的地震能量[1-3]。目前我国建成的隔震建筑已经超过一千栋,主要以多层建筑为主[4],其隔震装置绝大部分采用叠层橡胶隔震支座。叠层橡胶隔震支座分为不提供水平屈服力的橡胶隔震支座(LNR)和提供水平屈服力的铅芯橡胶支座(LRB)。对处于基本风压较大地区的隔震建筑,为满足抗风设计的要求,通常的设计方法是通过增加铅芯橡胶支座的数量,来提高隔震层的总屈服剪力抵抗风荷载[5-6]。因此产生的问题是,铅芯橡胶支座数量的增加,隔震层的水平刚度相应增大,结构隔震效果降低,这就削弱了隔震结构的优势。
目前,通过设置钢板抗风支座来满足隔振层抗风设计要求在工程中已有所应用,工程应用表明,钢板抗风支座具有安装简单、维修和更换方便、造价低廉等优点。我国现行GB50011-2010建筑抗震设计规范[7](简称《抗规》(10版))中又未涉及抗风支座与隔震支座共同工作的设计规定。 对于隔震层抗风承载力由铅芯隔震支座和钢板抗风支座两者共同承担的组合隔震形式,两者的设置数量、位置及设计参数对上部结构的减震效果的影响,目前的研究分析较少,钢板抗风支座工作机理也是一个需要研究的问题。因此本文以某实际隔震工程为背景,对上述问题进行分析与讨论,为工程应用提供合理建议。
背景工程为福建厦门市某5层的幼儿园,图1、图2给出底层平面图和剖面图。钢筋混凝土框架结构,基础隔震设计,填充墙为加气混凝土砖。建筑总长度(X向)为54.30 m,最小宽度(Y向)为14.40 m,结构高度为18.90 m,高宽比为1.32。建筑抗震设防烈度为7度,基本地震加速度为0.15 g,地震分组为第二组,场地类别为Ⅱ类,特征周期为0.40 s,基本风压为0.80 kN/m2,地面粗糙程度为B类。表1给出结构主要设计信息。
图1 底层平面图(单位:mm)
表1 结构设计信息
注:层6为屋顶梯间。
本文采用国际通用有限元分析软件ETABS建立结构模型,有限元分析模型如图3所示。分析模型包括了隔震层和上部结构的主要结构信息。框架梁、柱和楼板分别采用空间杆系单元和膜单元模拟,隔震支座采用软件自带的Isolator1连接单元模拟,橡胶隔震支座采用线性模型模拟,铅芯橡胶支座采用空间双向耦合的非线性恢复力模型模拟。表2给出隔震支座的型号及主要参数。
选取适合于Ⅱ类场地土的三组实测地震波El Centro波、Taft波、唐山波和一组根据场地地质条件人工合成的地震波。按照《抗规》(10版)规定验算四组地震波,验算结果表明,上述四组地震波可用于本工程设计,设计取地震响应的包络值为时程分析的代表值。地震作用沿X向和Y向分别输入上述四组波,双向输入的地震波记录按X∶Y的比例为1∶0.85。
图3 有限元分析模型
表2 隔震支座型号及主要参数
注:r表示隔震支座剪应变,橡胶剪切模量为G=0.46 N·mm-2。
《抗规》(10版)规定,结构水平向减震系数β为多遇地震作用下,隔震结构与抗震结构楼层最大水平剪力的比值,本工程考虑两种β值,分别是β值≤0.53,可满足设防烈度降半度要求;β值<0.40,可满足设防烈度降半度且可降低上部结构的抗震措施。以此建立在同比条件下的两种隔震方案。
优先考虑满足抗风设计要求进行隔震层设计的方案,同时要求水平向减震系数β值≤0.53,以此建立隔震方案一模型(简称隔震一),经过多轮优化设计,最终采用11个LNR500和33个LRB500隔震支座,支座平面布置如图4所示。
优先考虑上部结构减震系数较小,减震效果较显著的方案,同时要求水平向减震系数β值<0.40,以此建立隔震方案二模型(简称隔震二)。经过多轮优化设计,最终采用26个LNR400、18个LRB500隔震支座及2个抗风支座,隔震装置平面布置如图5所示。
隔震一和隔震二结构模型的计算都是经过多轮的时程分析,设计指标满足《抗规》(10版)的相关规定,分析结果详见本文第四节。经过综合比较,工程选择隔震方案二实施。
图4 隔震方案一支座平面布置 (单位:mm)
图5 隔震方案二支座平面布置 (单位:mm)
表3给出了两种方案的隔震层最不利方向(Y向)抗风承载力设计值验算结果。
表3 隔震层Y向抗风验算
注:表中VRW为抗风装置的水平承载力设计值;rW为风荷载分项系数1.4;VWR为风荷载作用下隔震层的水平剪力标准值。
由表3可知,两种隔震方案都满足抗风设计要求。
抗风支座由抗风钢板和连接钢板组成,其一端连接于上部结构(即隔震层框架梁上),另一端与下支墩拉梁相连接。抗风支座采用Q235软钢。将抗风支座设计成中间截面小,上下截面大的缺口型,保证结构在遭遇中震时抗风支座能够在设计屈服面屈服。经多次有限元模型优化设计,确定抗风钢板厚度为10 mm,有效截面长度100 mm,高度为260 mm,每组4块钢板,每个抗风支座的屈服力设计值为542 kN。
(1)模型基本情况:模型由4块10 mm厚的受力钢板及上下连接板组成,图6给出了抗风支座的安装现场实体模型。分析模型采用功能强大的有限元软件Abaqus建立,钢材采用Solid单元模拟,钢材本构模型采用多折线强化模型如图7所示。图7中强化阶段应变 为0.05,极限应变 为0.1,屈服强度设计值 为210 MPa,强化阶段强度设计值 为325 MPa,极限强度设计值 为380 MPa,泊松比为0.3。
图6 抗风支座实体模型
(2)边界条件及加载情况:钢板的3个平动自由度和3个转动自由度全部约束;加载方式采用静力加载。水平荷载直接施加在上连接钢板的单个网格单元会出现局部应力集中,且运算不收敛,因此在水平力施加之前,对上连接板所有单元耦合水平加载方向的平动自由度,施加一个平行于受力钢板的水平力542 kN。加载过程仅对耦合后的主节点加载,这种加载形式能很好的避免应力集中,且操作简便。
(3)计算结果:对Abaqus有限元模型进行静力加载分析,图8和图9分别给出了抗风支座在加载最终阶段的截面应力分布和位移情况。
图8 抗风支座应力
图9 抗风支座位移
由图8可得,加载后,抗风支座最大的Mises应力达到384 MPa,超过受力钢板的极限强度设计值,钢板已发生破坏。本文计算为了防止钢板破坏后计算不收敛,将受力钢板的极限强度设计值设定为400 MPa。因此抗风支座能够满足隔震结构的抗风要求,又能保证在中震时屈服。由图9可得,加载后,抗风支座的位移与加载方向一致,极限变形值为5.26 mm。
根据ETABS计算结果,隔震方案二的隔震层质心坐标为(29.40 m,11.20 m),刚心坐标为(29.47 m,11.60 m),偏心率为(Ex=0.13%,Ey=2.01%),偏心率小于3%的要求[8]。经多次试算,结合工程实际,同时考虑到抗风支座沿结构外围布置有利于结构的抗扭,最终确定两个抗风支座的位置坐标分别为(53.90 m,18.50 m)和(5.20 m,5.35 m),使得两个抗风支座承受的水平剪力对质心的矩能够基本自平衡。
隔震结构的运动微分方程为:
(1)
质量矩阵M、刚度矩阵K、阻尼矩阵C可以写成如下形式:
在ETABS建立的有限元模型中,选用plastic1单元模拟抗风支座,采用的Y向等效刚度为Abaqus设计计算结果(即抗风支座的设计剪力抗力与其极限变形之比)反算出来的等效刚度,X向与Z向的刚度对相应方向的刚度影响较小,忽略不计,假设为0值。抗风支座的阻尼亦假设阻尼为0值。
表4给出了隔震二结构在小震作用下隔震层的剪力设计峰值。表5给出了隔震二结构在小震及中震作用下隔震层Y向的位移峰值。
表4 小震作用下隔震层剪力 (单位:kN )
由表4得,隔震二结构在小震作用下隔震层的最大剪力值为1 972 kN小于隔震二结构的抗风承载力设计值2 310 kN;由表5得,隔震二结构在小震作用下隔震层Y向的位移峰值为3.58 mm小于抗风支座的极限变形值5.26 mm,说明抗风支座未破坏,参与工作,满足变形协调条件。在中震作用下隔震层Y向的最大位移值为36.39 mm大于抗风支座的极限变形值5.26 mm,说明抗风支座破坏,退出工作,内力重新分配到结构就相当于去除抗风支座的影响,在ETABS中通过删除plastic1单元模拟的抗风支座实现。结构在中震和大震作用下内力重新分配的结果即为结构在不设抗风支座的隔震结构的计算结果,如4.4节中所述。
表5 隔震层Y向位移(单位:mm)
结构在多遇(罕遇)地震作用下,分别取隔震支座100%(250%)剪应变所对应的等效水平刚度和等效阻尼比。对两种隔震方案结构模型与抗震结构模型进行模态分析,表6给出了模态分析结果。
由表6得,在设防烈度中震及大震作用下,隔震一结构较抗震结构的基本周期延长分别为3倍以上,隔震二较隔震一基本周期更长,按照加速度反应谱理论,周期越长,结构地震作用力越小,可见隔震二较隔震一隔震效果更明显。
表6 结构基本周期(单位:s)
分析抗震和两种隔震结构模型在设防烈度下各条地震波激励下的响应得到:①在中震作用下,隔震一和隔震二结构楼层剪力水平向减震系数分别为0.52和0.36,符合上述隔震方案的建立条件,可满足设防烈度降半度要求;②在大震作用下,隔震一和隔震二结构隔震层的位移分别为91 mm及143 mm都小于隔震支座的水平位移限值208 mm;抗震、隔震一和隔震二结构的上部结构最大层间位移角分别为1/104、1/305和1/343,相比较于抗震和隔震一结构,隔震二结构有更高的安全储备。③在大震作用下,隔震一和隔震二结构楼层加速度减震率最小值分别为57%和74%,隔震二结构的减震效果较隔震一结构好。
图10给出了中震下结构最大的楼层剪力,图11给出了大震下结构最大的层间位移响应,图12给出了大震下结构各楼层最大的加速度响应,图13给出了大震下结构顶层(五层)绝对加速度时程曲线。图10-图12中楼层0层表示隔震层。
图10 中震作用下楼层剪力
图13 大震作用下顶层绝对加速度时程
由图10可得,抗震结构楼层剪力最大,比较抗震结构,隔震结构的楼层剪力都有大幅度的降低,隔震二结构的楼层剪力最小。由图11可得,对于隔震层上部结构,隔震一与隔震二的层间位移相近,且远小于抗震结构层间位移;对于隔震层下部结构,隔震结构的层间位移也都在允许值之内,且都有较大的安全度。由图13可得,相比较于抗震结构,隔震结构顶层的加速度反应有较大幅度的减少。从各个指标综合分析,隔震二结构减震效果较隔震一结构略好。
(1) 本文对抗风支座设计进行了计算分析尝试,验证了通过叠层橡胶隔震支座和抗风支座的组合隔震形式能够实现背景工程的上部结构楼层剪力减震效果适中、隔震层位移在合理范围和隔震层抗风承载力满足要求三者的平衡。
(2) 基于Abaqus对抗风支座进行有限元模拟,验证了抗风支座能够满足设计要求,即结构在正常使用条件和小震作用下,抗风支座参与工作,提供水平抗力,隔震层不屈服;在结构遭遇中震及大震时,抗风支座能屈服并破坏,退出工作,不影响上部结构的减震效果。钢板抗风支座受力机理清晰,构造明确,可最大限度发挥隔震技术的优势。
对于基本风压较大地区的隔震结构,建议采用由叠层橡胶支座和抗风支座组合而成的隔震形式,并合理的选择铅芯隔震支座和钢板抗风支座的数量和布置来降低结构的地震响应,提高隔震结构隔震效果,从而进一步提高结构的安全性。
参 考 文 献
[1]Petti L, Giannattasio G, De Iuliis M, et al. Small scale experimental testing to verify the effectiveness of the base isolation and tuned mass dampers combined control strategy[J]. Smart Structures and Systems, 2010, 6(1):57-72.
[2]Jung H J, Eem S H, Jang D D, et al. Performance analysis of a smart base-isolation system considering dynamics of mrelastomers[J]. Journal of Intelligent Material Systems and Structures, 2011, 22(13): 1439-1450.
[3]商昊江,祁 皑.高层隔震结构减震机理探讨[J].振动与冲击,2012,31(4):8-12.
SHANG Hao-jiang, QI Ai. Isolation mechanism of high-rise isolation structrue [J]. Journal of Vibration and Shock, 2012, 31(4): 8-12.
[4]王曙光,陆伟东,刘伟庆,等.昆明新国际机场航站楼基础隔震设计及抗震性能分析[J].振动与冲击,2011,30(11):260-265.
WANG Shu-guang, LU Wei-dong, LIU Wei-qing, et al. Seismic isolation design and dynamic analysis of Kunming new international airport terminal [J]. Journal of Vibration and Shock, 2011, 30(11):260-265.
[5]杜东升,王曙光,刘伟庆,等.宿迁某高层办公楼基础隔震设计方法[J].建筑结构,2010,40(11):49-53.
DU Dong-sheng, WANG Shu-guang, LIU Wei-qing, et al. Design of base isolation for a high-rise office building in Suqian [J]. Building Structure, 2010, 40(11):49-53.
[6]Karabork T . Performance of multi-storey structures with high damping rubber bearing base isolation systems[J]. Structural Engineering and Mechanics, 2011, 39(3):399-410.
[7]GB50011-2010.建筑抗震设计规范[S].北京:中国建筑工业出版社,2010.
[8]日本建筑学会.隔震结构设计[M].刘文光,译.北京:地震出版社,2006.