考虑摆动的层间隔震结构自振特性及随机响应分析张尚荣;毛明杰;刘芳兰;包超【摘要】采用铅芯叠层橡胶支座的层间隔震结构,由于结构自重和活载产生的竖向荷载较大,当叠层橡胶层的总厚度较大时,隔震支座在产生水平变形的同时,也会产生竖向压缩变形,导致上部子结构产生摆动效应,同时对下部子结构的响应也产生影响,进而影响层间隔震结构的隔震效果.根据层间隔震结构的自振特性将其简化为两质点模型和多质点模型,分析层间隔震结构考虑摆动效应时的固有自振特性,以及采用虚拟激励法,运用均匀调制非平稳随机响应分析方法,分析隔震层摆动对层间隔震结构的振动响应影响.结果表明,考虑层间隔震结构的摆动效应,将对各子结构的响应产生不同的影响;考虑摆动效应的不同质量比的层间隔震结构,自振特性变化规律存在一界限刚度比K;随着隔震层刚度比的增大,考虑摆动效应的层间隔震结构各子结构的响应逐渐趋近于未考虑摇摆特性结构的响应.【期刊名称】《兰州理工大学学报》【年(卷),期】2018(044)006【总页数】7页(P118-124)【关键词】层间隔震结构;摆动效应;自振特性;随机响应【作者】张尚荣;毛明杰;刘芳兰;包超【作者单位】宁夏大学土木与水利工程学院,宁夏银川 750021;宁夏大学土木与水利工程学院,宁夏银川 750021;宁夏大学土木与水利工程学院,宁夏银川750021;宁夏大学土木与水利工程学院,宁夏银川 750021【正文语种】中文【中图分类】TU352.12目前,关于隔震结构的研究大部分都是仅考虑隔震子结构的平动,而忽略了在平动过程中由于隔震层的不均匀压缩竖向变形引起的上部子结构的摆动.当上部子结构层间刚度较小,垂直荷载较大,且其所采用的多层橡胶垫的橡胶总厚度较大时,就有可能产生明显的竖向变形[1].当竖向变形不均匀时,就产生了明显的摆动现象. 国内外对隔震结构摆动方面的研究还较少.曾奔等[2]基于锥体模型,用子结构法对基础隔震结构体系进行平摆耦合的结构模型分析,推导了在地震作用下结构平摆耦合的传递函数;黄丽彬[3]研究了摆动对基础隔震结构的振动控制的影响;徐培淞[4]基于拟静力试验方法研究了铁路桩基础摇摆隔震桥墩的抗震性能.本次研究分别以层间隔震结构考虑隔震层摇摆特性的两质点简化模型和多质点模型为研究对象,分别分析隔震层的摆动效应对结构自振特性的影响;在时频域内从结构动力响应、时变方差和谱密度分布方面,分析隔震层摆动对结构响应的影响.研究结果可在隔震体系响应的控制、结构优化等方面提供指导作用.1 平摆动体系动力分析模型1.1 多质点平摆动体系动力分析模型图1所示为多质点平摆动体系层间隔震结构动力分析模型.图中,mi为隔震层的质量;ki为隔震层的水平总刚度;ci为隔震层阻尼系数;ms,i、ks,i、cs,i分别为下部子结构第i层的质量、刚度和阻尼系数;mu,i、ku,i、cu,i分别为上部子结构第i 层的质量、刚度和阻尼系数;Hu,i为上部子结构第i层至地面的垂直距离;Hi为隔震层至地面的垂直距离;θ为隔震层的摆动转角.根据图1所示多质点平摆动体系模型,上部子结构第i层的位移响应为图1 多质点平摆动体系动力分析模型Fig.1 Dynamic analysis model of multi-particle planar swing systemxU,i=xg+xu,i+xθ,i=xg+xu,i+θ(Hu,i-Hi)(1)式中:xg为地震引起的相对位移;xu,i、xθ,i分别为上部子结构第i层和隔震层摆动而引起的相对位移.隔震层位移响应为x=xg+xi(2)下部子结构位移响应为xS,i=xg+xs,i(3)由D’Alembert原理,各子结构的运动方程可表示为式中:为上部子结构第i层与地面的相对位移、速度和加速度反应;为地震输入加速度;mu、mu,i分别为上部子结构和上部子结构第i层的质量;ku、cu分别为上部子结构水平刚度和阻尼系数;为由于隔震层摆动而引起的上部子结构的相对位移、速度和加速度反应;为隔震层相对位移和速度响应;ki为隔震层水平刚度;为下部子结构与地面之间的水平相对位移、速度和加速度响应.由式(4~6),并结合下部子结构的运动方程可列出结构平摆动体系的运动方程为(7)式中:X={xs,1,xs,2,…,xs,n,xi,xu,1,xu,2,…,xu,n,θ},xs,i为下部子结构第i层位移;I为单位列向量;其中,转动刚度为橡胶支座数,kv,i为第i个支座竖向刚度.1.2 两质点模型考虑隔震层水平平动和摇摆转动,层间隔震结构两质点平摆动体系等效简化动力分析模型如图2所示.图2 两质点平摆动体系动力分析模型Fig.2 Dynamic analysis model of two-particle planar swing system图2中,mu为上部子结构质量;k0为隔震层水平刚度;ku、cu分别为上部子结构隔震层水平刚度和阻尼系数;ms为下部子结构质量;ks、cs分别为下部子结构水平刚度;hu、hs分别为上部子结构和隔震层至地面的垂直距离;Kθ为隔震层的转动刚度.有线弹性情况下,层间隔震结构两质点平摆动体系的运动方程可表达为:式中:转动刚度为橡胶支座的个数,kv,i为第i个支座的竖向刚度.2 平摆动体系自振特性根据两质点简化模型,定义上、下部子结构的质量比:μ=mu/ms(8)上、下部子结构的自振圆频率分别为假定层间隔震平摆动体系简化为两支座平面模型,其隔震层转动刚度与平动刚度之比为两质点平摆动无阻尼体系的自振频率特征方程为式中:由上式得到平摆动体系振型参与系数γn和模态质量参与系数ηn,分别为假定体系平面规则,体系水平宽度为L=1.0 m,隔震层刚心与结构平面形心重合,平摆动体系摆动角度较小.图3给出了不同质量比下的最优参数模型层间隔震平摆动体系的固有周期与未考虑平摆动体系固有周期比的关系.其中实线代表体系第一阶固有周期的比值,虚线代表体系第二阶固有周期的比值.由图中可以看出当层间隔震平摆动体系的刚度比kv /ks =250时,Tθ /T均趋近于1.0,即当层间隔震平摆动体系竖向刚度大于平动刚度250倍,两者前两阶固有周期接近.图3 平摇摆体系的固有周期特性Fig.3 Natural periodicity of swing system图4给出了不同质量比下层间隔震平摆动体系最优参数模型振型参与特性随隔震层刚度比变化规律.图5给出了不同质量比下层间隔震平摆动体系的最优参数模型模态质量参与特性随隔震层刚度比的变化规律.μ = 2.0时,体系一阶振型主要表现为上部子结构自由平动γ1Φp1=1.0,下部子结构基本不动γ1Φb1=0、γ1Φθ1=0;当刚度比小于50时,二阶振型表现为结构摇摆振动γ2Φθ2=1.0;当刚度比大于50时,二阶振型表现为下部子结构自由平动γ2Φp2=1.0;μ = 8.0时,体系一阶振型主要表现为上部子结构自由平动下部子结构基本不动当刚度比小于250时,二阶振型表现为结构的摇摆振动γ2Φθ2=1.0;当刚度比大于250时,二阶振型表现为下部子结构的自由平动图4 层间隔震平摆动体系的振型参与特性Fig.4 Modal participation features of inter-story vibro-isolation planar swing system图5 层间隔震平摆动体系的模态质量参与特性Fig.5 Modal mass participation features of inter-story vibro-isolation planar swing system图5表现的规律与图4所表现的振型参与特性有相同的规律特性,此处不再赘述.3 平摆动体系随机响应分析3.1 虚拟激励法假设地震作用随机过程为非平稳过程,采用均匀调制函数来模拟地震动的非平稳过程.假设均匀调制演变随机激励具有以下形式[5]:f(t)=g(t)x(t)(13)式中:g(t)为均匀调制函数;x(t)为零均值的平稳随机过程.构造虚拟激励:(14)则在t时刻所产生的结构响应为(15)其中,I(w,t)可由杜哈梅积分得到:I(w,t)=h(t-τ)g(τ)eiwτdτ(16)因此,结构响应的自功率谱密度为(17)则响应的时变方差为(18)3.2 随机地震动模型随机地震动模型采用Kanai-Tajimi谱,其表达式为[6]:其中:S0为基岩加速度自谱密度,谱强度S0=9.537 cm2/s3;ωg和ζg分别为场地土的卓越圆频率和阻尼比,Ⅱ类场地土的ωg=17.95 rad/s,ζg=0.72.均匀调制函数g(t)取三段式时间包络函数.其中,tb=1.0 s,tc=8.0 s,c1=0.60,t=20 s.4 算例分析根据文献[7-8]选取最优参数模型:上、下部子结构分别为12层和6层,层质量均为2.0×106 kg,层刚度均为4.2×109 N/m,阻尼比均为0.05,层高均为3.0 m;隔震橡胶垫采用直径为D=1 000 mm的铅芯橡胶支座,隔震层质量为2.0×106kg,层刚度为8.0×107 N/m(质量比μ =2.0,频率比fr =0.164 6,等效阻尼比ξ=0.25).隔震结构所在场地为8度设防,Ⅱ类场地.图6、图7对比分析了不同刚度比时,考虑摇摆与未考虑时隔震结构在非平稳随机激励下的响应时变方差,结果表明在考虑隔震结构刚度差异情况下产生的结构摇摆放大了结构响应,且随着摇摆刚度的增大,各响应逐渐趋近于未考虑摇摆时结构的响应.图6 摇摆隔震结构响应时变方差(kθ /ks=9.0×104)Fig.6 Response time-varied variance of swing vibro-isolation structure(kθ/ks =9.0×104)图7 摇摆隔震结构响应时变方差(kθ /ks =9.675×105)Fig.7 Response time-varied variance of swing vibro-isolation structure (kθ /ks =9.675×105)图8所示分别为层间隔震结构各子结构响应的功率谱密度.图中,K1和K2分别表示刚度比为9.0×104和9.675×105对应摇摆隔震结构,结构各响应功率谱密度曲线出现峰值的点对应结构自振频率点,在这些频率点上,地震动输入频率与结构自振频率相等产生共振,结构响应较大.当刚度比为K1时,即结构的刚度比为9.0×104时,摇摆隔震结构的响应大于未考虑摇摆时结构响应,随着刚度比增大,响应逐渐增大;当刚度比达到K2时,摇摆隔震结构响应趋近于未考虑摇摆隔震时结构的响应.图8 摇摆隔震结构响应功率谱密度Fig.8 Response power spectral density of swing vibro-isolation structure为验证随机分析结论的准确性,采用时程分析方法对其进行验证分析.选取EI-Centro波作为结构激励,地震动峰值调整为90gal,图9、图10所示分别为刚度比为9.0×104和9.675×105时所对应的摇摆隔震结构和未考虑摇摆隔震结构的上部子结构顶层位移及下部子结构顶层绝对加速度时程响应.当刚度比较小时,摇摆隔震结构分析结果大于未考虑摇摆时结构的响应;随着刚度比的增大,两者响应差距逐渐减小;当刚度比为kθ/ks = 9.675×105时,两种结构的响应近似相等.图9 摇摆隔震结构时程响应分析(kθ/ks=9.0×104)Fig.9 Time history response analysis of swing vibro-isolation structure (kθ/ks = 9.0×104)图10 摇摆隔震结构时程响应分析(kθ/ks = 9.675×105)Fig.10 Time history response analysis of swing vibro-isolation structure (kθ /ks =9.675×105)5 结论本次研究建立并分析了两自由度简化模型和多自由度平摆动层间隔震结构的动力特性、非平稳随机激励下结构的响应及时频域特性,得到以下结论.1) 不同质量比的层间隔震结构,在考虑摆动效应的情况下,自振特性变化规律存在一界限刚度比K.当刚度比小于K时,结构自振特性曲线存在明显变化趋势;当刚度比大于K时,结构自振特性曲线保持稳定.2) 考虑摆动效应的层间隔震结构,对各子结构响应都有不同程度的影响,对上部子结构顶层位移、隔震层位移和下部子结构顶层绝对加速度有不同程度的放大作用.3) 随着隔震层刚度比的增大,考虑摆动效应的层间隔震结构各子结构的响应逐渐趋近于未考虑摇摆特性结构的响应,且存在一临界刚度比K,K= 9.675×105,当隔震层刚度比大于此值,可以不考虑摆动效应的影响.致谢:本文得到宁夏大学引进人才科研项目(BQD2015004)的资助,在此表示感谢. 参考文献:【相关文献】[1] 周福霖.工程结构减震控制 [M].北京:地震出版社,1997.[2] 曾奔,周福霖,黄东阳.土-结构相互作用下考虑平摆耦合的基础隔震体系分析 [J].西安建筑科技大学学报,2008,40(4):538-543.[3] 黄丽彬,邹立华.基础隔震结构考虑摆动的振动控制研究 [J].世界地震工程,2010,26(1):212-218.[4] 徐培淞.铁路摇摆隔震桥墩抗震性能试验研究 [D].兰州:兰州交通大学,2015.[5] 林家浩,张亚辉.随机振动的虚拟激励法 [M].北京:科学出版社,2004.[6] 薛素铎,王雪生,曹资.基于新抗震规范的地震动随机模型参数研究 [J].土木工程学报,2003,36(5):5-10.[7] 张尚荣,谭平.基于NSGA-Ⅱ的层间隔震体系参数优化与能量响应分析 [J].合肥工业大学学报,2015,38(7):949-954.[8] 张尚荣,谭平.层间隔震体系可靠度的灵敏度分析 [J].振动、测试与诊断,2016,36(1):102-107.。