水平地震作用下浮放柱体动态响应研究.pdf

【目的】浮放柱体是传统木、石建筑的典型特征,是其核心受力构件之一,揭示浮放柱体在水平地震作用下的响应特性至关重要。【方法】基于传统木、石建筑结构特征与柱顶构件偏心实际情况建立了柱体二维模型,推证了不同动态响应触发临界条件,依据触发临界条件阐释了柱体动态响应触发敏感度。【结果】结果表明:水平地震作用下决定柱体动态响应类型的因素包含基础平面与柱底面间静摩擦系数、柱体长细比、柱顶构件偏心率;在可被触发的抬升、先抬升后滑移及抬升耦合滑移三种柱体动态响应中,抬升最易被触发,抬升后的摇摆运动是柱体的主要运动形式;柱体长细比及柱顶构件偏心率的改变均能影响柱体抬升响应触发条件,但后者的影响更显著;柱体几何特征、柱顶构件质量与其偏心程度均对柱体摇摆运动状态有影响。研究为传统建筑状态评估与保护提供必要理论依据。关键词:浮放柱体;偏心;水平地震作用;动态响应;临界条件中图分类号:TU 3 6 6.2 文献标识码:AD O I:1 0.1 6 3 5 5/j.c n k i.i s s n 1 0 0 7-9 4 3 2 t y u t.2 0 2 3.0 4.0 1 6 文章编号:1 0 0 7-9 4 3 2(2 0 2 3)0 4-0 7 1 6-0 9D y n a m i c R e s p o n s e o f F r e e S t a n d i n g C o l u m n s u n d e r H o r i z o n t a l E a r t h q u a k eG A O C h a o,L I Y u x i n(I n s t i t u t e o f S c i e n c e a n d T e c h n o l o g y R e s e a r c h,C h i n a T h r e e G o r g e s C o r p o r a t i o n,B e i j i n g 1 0 0 0 3 8,C h i n a)A b s t r a c t:【P u r p o s e s】F r e e s t a n d i n g c o l u m n s a r e o n e o f t h e t y p i c a l f e a t u r e s o f t r a d i t i o n a l t i m-b e r a n d s t o n e b u i l d i n g s.C o l u m n s a r e t h e k e y c o m p o n e n t s f o r t r a n s m i t t i n g v e r t i c a l l o a d s i n t r a d i-t i o n a l s t r u c t u r e s.D e a d w e i g h t o f t h e s u p e r s t r u c t u r e t r a n s m i t s t o c o l u m n s t h r o u g h t h e c o m p o n e n t a t t h e t o p o f t h e c o l u m n s.T h e d y n a m i c r e s p o n s e c h a r a c t e r i s t i c s o f f r e e s t a n d i n g c o l u m n s s u b j e c-t e d t o h o r i z o n t a l e a r t h q u a k e s a r e a n a l y z e d i n t h i s r e s e a r c h.【M e t h o d s】A t w o-d i m e n s i o n a l r i g i d b o d y m o d e l o f a c o l u m n w a s e s t a b l i s h e d o n t h e b a s i s o f t h e s t r u c t u r a l c h a r a c t e r i s t i c s a n d e c c e n-t r i c i t y o f t h e c o m p o n e n t a t t h e t o p o f t h e c o l u m n.T h e t r i g g e r i n g c r i t i c a l c o n d i t i o n s f o r d i f f e r e n t d y n a m i c r e s p o n s e s w e r e d e d u c e d.T h e t r i g g e r s e n s i t i v i t y o f t h e d y n a m i c r e s p o n s e s w a s e x p l a i n e d a c c o r d i n g t o t h e t r i g g e r i n g c r i t i c a l c o n d i t i o n s.【F i n d i n g s】R e s u l t s s h o w t h a t f a c t o r s d e t e r m i n i n g t h e t y p e o f d y n a m i c r e s p o n s e a r e s t a t i c f r i c t i o n c o e f f i c i e n t b e t w e e n t h e f o u n d a t i o n a n d t h e b o t t o m s u r f a c e o f t h e c o l u m n,t h e s l e n d e r n e s s r a t i o o f t h e c o l u m n,a n d t h e e c c e n t r i c i t y o f t h e c o m p o n e n t a t t h e t o p o f t h e c o l u m n.T h e r e a r e t h r e e t y p e s o f d y n a m i c r e s p o n s e t h a t c a n b e t r i g g e r e d o n t h e c o l u m n,i n c l u d i n g u p l i f t,u p l i f t a n d t h e n s l i p,a n d s i m u l t a n e o u s u p l i f t-s l i p.Am o n g t h e m,u p l i f t i s t h e m o s t l i k e l y t o b e t r i g g e r e d,a n d r o c k i n g a f t e r u p l i f t i s t h e m a i n f o r m o f m o t i o n o f t i m b e r c o l u m n s.T h e t r i g g e r i n g c o n d i t i o n s o f t h e c o l u m n u p l i f t c a n b e a f f e c t e d b y t h e s l e n d e r n e s s r a t i o o f t h e t i m b e r c o l u m n a n d t h e e c c e n t r i c i t y o f t h e c o m p o n e n t a t t h e t o p o f t h e c o l u m n,b u t t h e l a t-t e r h a s a m o r e s i g n i f i c a n t e f f e c t.T h e g e o m e t r i c c h a r a c t e r i s t i c s o f t h e c o l u m n,t h e m a s s o f t h e c o l u m n t o p m e m b e r a n d i t s e c c e n t r i c i t y h a v e e f f e c t s o n t h e r o c k i n g m o t i o n o f t h e c o l u m n.T h i s s t u d y p r o v i d e s t h e n e c e s s a r y t h e o r e t i c a l b a s i s f o r t h e s t a t e a s s e s s m e n t a n d p r o t e c t i o n o f a n c i e n t t i m b e r s t r u c t u r e s.K e y w o r d s:f r e e s t a n d i n g c o l u m n;e c c e n t r i c i t y;h o r i z o n t a l e a r t h q u a k e;d y n a m i c r e s p o n s e;c r i t i c a l c o n d i t i o n 以浮放方式置于基础平台之上的物体称为浮放物体1,物体与基础平台间无典型、有效的机械连接。地震是影响传统建筑存续的关键因素之一2。在水平地震作用下,浮放物体或保持静止,或被触发摇摆、滑移、摇摆-滑移耦合等运动3。柱体是建筑结构核心受力构件之一,柱体在地震作用下的稳定性直接决定了传统建筑的整体稳定性,准确地预测其在外荷载作用下的响应是科学保护传统建筑的前提4,因而开展水平地震作用下浮放柱体动态响应研究是传统建筑状态评估与修缮保护的必要环节。近年来与浮放柱体力学特性有关的研究已广泛开展。WANG e t a l5考虑柱顶荷载的影响对柱脚节点的受力性能进行了研究,通过理论分析阐释了木柱摇摆机理和柱脚受力状态,建立了柱脚节点的解析模型。姚侃等6提出了柱与础石摩擦滑移临界条件,建立了柱脚隔振体系模型。MA E D A e t a l7基于对日本木结构古建筑的木柱摇摆特性理论分析指出,木柱的刚度和强度增加并不会减小地震作用下其位移的响应。针对浮放物体动力学特性的体系化研究已持续逾半个世纪。HOU S N E R8将浮放物体简化为刚性二维块体模型,应用牛顿第二定律推导了分段微分方程以描述浮放物体摇摆动态响应,被学界普遍认为是现代推动浮放物体动力特性分析研究进程的重要节点。S HE NTON e t a l9通过解析形式研究了水平激励下浮放物体的多种运动形式。T AN-I GU CH I1 0经推证得出,低频长持时地震会增加浮放物体发生大滑移与倾覆的风险。KOUNA D I S1 1对浮放物体的稳定性进行了研究,在假设块体与基础间无滑移的基础上分析了正弦激励地震作用下的不同倾覆模式。杨维国等1 2将浮放物体简化为木制块体,利用振动台试验与数值模拟相结合的方法,揭示了宽高比越大试件越稳定,同时阐明宽高比的变化对浮放物体的滑移量无影响,但摇摆角随宽高比的增大而减小。上述研究揭示了柱体摇摆、滑移等动态响应特征,为相关研究提供了必要支撑。但分析所涉动态响应的选取、分析目标的确立大多基于概念分析,缺乏定量的科学判断。而相关的试验研究也仅对试验现象进行了定性归纳,未能系统阐释不同类型动态响应的特征和临界条件。此外,尽管上述研究中大都考虑了柱顶荷载,但均忽略了荷载偏心的可能性。为此,本文系统分析了存在柱顶构件偏心情况的浮放柱体在水平地震作用下的动态响应,阐明了不同动态响应的触发临界条件,揭示了柱体的动态响应特征。1 基本假定与简化模型部分传统建筑的柱体以浮放方式竖立于基础之上,不设置典型、有效的机械连接,以国外石质传统建筑与国内木构古建筑最具代表性,如图1所示。此类建筑的柱体既可在基础平台上滑动,也可自由抬升、摇摆,形成滑动铰支座节点。梁柱节点间缝隙1 3在传统建筑中常见,这使得节点在承受荷载之前处于松散状态,导致柱体在承受荷载之初不受节点约束。石梁、木梁与斗栱等柱顶构件作为传力单元,将上部结构自重传递至柱体,致使柱顶构件与柱体顶部接触面间产生极大的摩擦力阻碍柱顶构件滑动。由于结构布局、建造误差、尺寸不规等因素导致柱顶构件常偏心布置,即柱顶构件的质心与柱体质心连线不与地面垂直是普遍存在的。（a）佛光寺东大殿（中国）（b）锡德古城（土耳其）图1 典型浮放传统建筑F i g.1 T y p i c a l f r e e s t a n d i n g t r a d i t i o n a l b u i l d i n g s石质构件由于刚度较大,难于在震后发生塑性破坏。而木构传统建筑在经历强震后尽管存在一定程度的位形改变,但木构件大都无明显塑性变形,更鲜有构件坏损出现1 4。因而针对其稳定问题分析717 第4期 高 潮,等:水平地震作用下浮放柱体动态响应研究时可忽略构件变形问题。水平地震作用下柱体滑移属于平动,摇摆属于转动轴周期性变换的定轴转动。由屋盖、梁构件等组成的上部结构的刚度极大且质量远大于柱体自重,这使得柱体摇摆时上部结构始终与地面平行,且柱顶构件与柱头之间的相对滑移极小1 5。基于所述浮放传统建筑结构特点及动力响应特征,当不考虑结构平面外位形变化时,对柱体动态响应分析做出如下基本假定:1)分析模型由柱体、基础、柱顶构件三部分组成,建筑上部结构的全部自重集中于柱顶构件并由其传递至柱体;2)全部构件均为二维匀质刚体;3)忽略柱顶构件底面与柱头间相对滑移;4)柱脚抬升时,柱顶构件始终与地面保持平行。取最具代表性的偏心情况作为分析对象,即柱顶构件一段偏移歪闪出柱顶范围之外,另一端位于柱顶面之内,简化后的分析模型如图2(a)所示。当柱脚分别向负向、正向抬升时,受柱顶构件偏心情况影响,柱顶构件的支点分别为图2(b)、图2(c)中的点N、M.图中其他符号的含义如下:(0)为上部结构质量与柱体质量比;L、R分别为点N、M与柱脚点L、R的连线同竖直方向的夹角;dL、dR分别为点N、M与柱脚点L、R连线的距离;L、R分别为柱顶构件重心投影与柱脚点L、R的连线同竖直方向的夹角。此外,lL、lR分别代表柱顶构件质心与柱脚点L、R连线的距离。定义柱顶构件质心相对于柱顶中点的偏移率为e,则e=b/rs i n.(1)式中:r为柱重心与柱脚点(点L或R)连线的距离;为上述连线与竖直方向的夹角;b(|b|rs i n)为柱顶构件偏移量,即柱顶构件重心点B与柱顶部中心点A在水平方向的距离,当点B在点A左侧(向X轴负方向偏移动)时b为负,此时-1e0;当点B在点A右侧(向X轴正方向偏移动)时b为正,CmgCmgLRLRzg zg LR基础平面XY柱体柱顶构件er sinmgNNAACRLdLdRrBMBBMmgA（a）物理量标注（b）负向抬升（c）正向抬升图2 柱体分析模型F i g.2 A n a l y s i s m o d e l o f t h e c o l u m n此时0e1.2 动态响应判定与控制方程在水平地震作用下,柱体可以被触发的动态响应包括原位静止、抬升、滑移和滑移-抬升4种。下面分别对其受力状态展开分析,并得到不同动态响应的判定条件。2.1 静止状态令作用于地面的水平地震作用加速度为z g,且假设其初始方向与X轴正方向相同。柱体受过质心的重力m g、惯性力mz g,柱顶构件自重m g、惯性力mz g,柱底面中心的支反力FN及与柱体位形变化趋势相反的静摩擦力fs,受力状态如图3(a)所示。mzg er sinB（a）静止状态（b）滑移状态Ber sinmzg AACCmzg mzg mgmgfsfsLRFNLRFNu图3 柱体静止及滑移状态受力分析图F i g.3 F o r c e a n a l y s i s o f s t a t i o n a r y a n d s l i p令柱体柱脚抬升产生虚位移,则由虚功原理可得mz grc o s-m g r ss i n+mz glLc o s-m g lLs i n=0.(2)化简得z g=rs i n+lLs i nLrc o s+lLc o sLg(3)式中lLs i nL=(1+e)rs i n.(4)lLc o sL=2rc o s.(5)将式(4)、(5)代入式(3)可得柱体抬升的临界水平加速度z g=1+(1+e)1+2gt a n(6)由水平和竖直方向柱体力的平衡条件FX=0、FY=0可得fs-(1+)mz g=0,FN-(1+)m g=0.(7)式中:fs为静摩擦力,FN为支持力,且fs满足fsfs,m a x=s(1+)m g.(8)817太 原 理 工 大 学 学 报 第5 4卷 式中:fs,m a x为摩擦力,s为静摩擦系数。综合式(6)、(8)可得柱体静止的临界条件为z g1+(1+e)1+2gt a n.(9)z gsg.(1 0)式(9)中,右侧因式(1+e)/(1+2)取值与偏心率e(-1e1)有关而与荷载比无关,故该因式取值范围为01+(1+e)1+21.(1 1)传统建筑中,木柱、石柱的高径比大都居于71 0之间1 6-1 7,即满足关系0.1 0 t a n0.1 5.(1 2)则式(9)右侧满足范围01+(1+e)1+2gt a n0.1 5g.(1 3)综合石、木两类传统建筑石基与木柱、石基与石柱间的静摩擦系数,s应满足6,1 80.4s0.8.(1 4)则式(1 0)右侧满足范围0.4gsg0.8g.(1 5)对比式(1 3)、(1 5)可得1+(1+e)1+2gt a nsg.(1 6)恒成立。故若柱体不抬升则一定静止,即柱体静止临界条件唯一,如式(9)所示。同理,若水平地震加速度的初始方向与X轴负方向相同,则柱体静止的临界条件式(9)变形为z g1+(1+e)1+2gt a n.(1 7)2.2 滑移状态当水平地震作用产生的惯性力增加到超过柱脚接触面最大静摩擦力但不导致柱体抬升时柱体发生滑移。滑移柱体动态响应如图3(b)所示,图中u为滑移距离。由式(1 0)可知,柱体发生滑移临界条件为z gsg.(1 8)结合柱体不被触发抬升的临界条件式(9)可得柱体滑移临界条件为sgz g1+(1+e)1+2gt a n.(1 9)由式(1 6)可知,式(1 9)恒不成立,z g的初始方向与X轴负方向相同时亦然。故只滑移不抬升状态恒不可能发生。由达朗贝尔原理,当柱体受荷向左侧滑动时相当于受到向右的惯性力,根据其在水平方向力的平衡关系即可得到柱体此时的动态响应方程kg+u-z g=0.(2 0)式中:k为动摩擦系数。2.3 抬升状态令水平地震加速度的初始方向与X轴正方向相同,柱体被触发抬升但其所受惯性力不超过柱脚接触面的最大静摩擦力,抬升柱体动态响应如图4(a)所示,图中为柱脚抬升角度。在柱体抬升过程中,柱顶构件与柱头的接触方式由线接触转换为点接触,接触点为柱顶构件底角N点。mzg B（a）抬升ACmgLRFN+L+dLmzg mgNmzg BACmgLRFN+L+dLmzg mgNfsrmzg BACmgLRFN+L+dLmzg mgNfsru（b）抬升后滑移（c）抬升同时滑移图4 与抬升相关的状态受力分析图F i g.4 A n a l y s i s m o d e l o f t h e c o l u m n a s s o c i a t e d w i t h u p l i f t由式(6)可得柱体抬升临界条件为z g1+(1+e)1+2gt a n.(2 1)结合柱体不发生滑移的临界条件式(1 0)可得柱体仅抬升不滑移的临界条件为1+(1+e)1+2gt a nz gsg.(2 2)柱脚抬升后柱体开始摇摆,有唯一自由度,系统部分有势力,其动态响应控制方程可由拉格朗日方程求得ddt La g -La g=Q.(2 3)式中:La g为拉格朗日量,Q为广义力,La g=T-V,其中T为动能,V为势能。其中动能T=12J2+12m(x2B+y2B).(2 4)式中:J为柱体绕柱脚点L或R在平面内转动的转动惯量。J=4m r2/3.(2 5)xB、yB分别为柱顶构件在水平、竖直方向的位移:xB=dLs i nL-s i n(L+).(2 6)yB=dLc o s(L+)-c o sL.(2 7)则其在水平、竖直方向的位移xB、yB分别为xB=-dLc o s(L+).(2 8)917 第4期 高 潮,等:水平地震作用下浮放柱体动态响应研究yB=-dLs i n(L+).(2 9)将式(2 5)、(2 8)和(2 9)代入式(2 4)可得柱体及柱顶构件动能为T=23m r22+12m d2L2.(3 0)势能为V=m g rc o s(+)-c o s+m g dLc o s(L+)-c o sL.(3 1)广义力为Q=W=-m r z gc o s(+)-m dLz gc o s(L+).(3 2)式中:W是水平地震作用下在位移增量为 时柱体和柱顶质点的虚功。将式(3 0)-(3 2)代入式(2 3)可得柱体摇摆控制方程(4r2/3+d2L)-g rs i n(+)-g dLs i n(L+)=-r z gc o s(+)-dLz gc o s(L+).(3 3)考虑式(1 2)柱体几何尺寸因素可将式(3 3)化简为(4r2/3+d2L)-g r(+)-g dL(L+)+(r+dL)z g=0.(3 4)相似地,若水平地震加速度的初始方向与X轴负方向相同,则柱体被触发抬升的临界条件变为1+(1-e)1+2gt a nz gsg.(3 5)化简后的摇摆响应控制方程为(4r2/3+d2R)-g r(-)+g dR(R-)+(r+dR)z g=0.(3 6)2.4 滑移-抬升状态当滑移、抬升2种动态响应耦合时,根据其发生的先后顺序可分为3种不同情况:先滑移后抬升、先抬升后滑移、滑移与抬升同步发生。但由式(1 6)可知,滑移先于抬升的动态响应恒不能发生,故仅对后2种情况予以分析。2.4.1 先抬升后滑移若水平地震加速度的初始方向与X轴正方向相同,且当水平地震作用产生的惯性力增加到导致柱体抬升时柱体开始抬升,随水平地震作用持续增加至惯性力增加到超过柱脚接触面最大静摩擦力,柱体开始滑移,并在此后保持抬升与滑移2种响应叠加。柱体的动态响应如图4(c)所示。柱体抬升临界条件同式(2 1).对柱体即将滑移的瞬间进行受力分析,由平衡条件FX0、FY=0可得fs,m a x-(1+)mz g+m r c o s(+)+2s i n(+)+m dL c o s(L+)+2s i n(L+)0,FN-(1+)m g-m r s i n(+)-2c o s(+)-m dL s i n(L+)-2c o s(L+)s1+(1+)g+r s i n(+)-2c o s(+)+dL s i n(L+)-2c o s(L+)+11+r c o s(+)+2s i n(+)+dL c o s(L+)+2s i n(L+).(3 8)类似地,若水平地震加速度的初始方向与X轴负方向相同,则柱体发生先抬升后滑移的临界条件为z g1+(1-e)1+2gt a n.(3 9)z gs1+(1+)g+r s i n(-)-2c o s(-)+dR s i n(R-)-2c o s(R-)+11+r c o s(-)+2s i n(-)+dR c o s(R-)+2s i n(R-).(4 0)2.4.2 滑移与抬升同步发生由式(1 6)可知,若柱体滑移发生则一定耦合抬升发生,则滑移与抬升同步发生的临界条件即为柱体滑移的临界条件,如式(1 8)所示。柱体的动态响应如图4(c)所示。滑移-抬升响应耦合时柱体有平动、转动2个自由度。令当柱体抬升角度时滑移位移为u,系统仍为部分有势力,其动态响应方程可由拉格朗日方程求得ddt La gqi -La g qi=Qi.(4 1)式中:qi(i=1,2)为广义坐标和u.柱体与柱顶构件的动能为027太 原 理 工 大 学 学 报 第5 4卷 T=(2r2/3+d2L/2)m2+1+2m u2+m r uc o s(+)+dLm uc o s(L+).(4 2)势能同式(3 1)一致。虚功为 W=-m r z gc o s(+)-m dLz gc o s(L+)+-(1+)mz g-fk u.(4 3)式中的滑动摩擦力fk满足fk=k(1+)m g+m r s i n(+)-2c o s(+)+m dL s i n(L+)-2c o s(L+).(4 4)式(4 4)中,滑动摩擦力fk是柱体抬升角度的因变量。为简化,在计算滑动摩擦力时可不考虑由于柱体抬升导致的竖直方向的荷载变化,而以等效摩擦系数k E替代k(k Ek)1 9,则简化后的滑动摩擦力为fk=k E(1+)m g.(4 5)将式(3 1)、(4 2)、(4 3)、(4 5)代入式(4 1)可得到柱体摇摆、滑移时的控制方程分别为(4r2/3+d2L)+r u c o s(+)-r us i n(+)+dLu c o s(L+)-dLus i n(L+)+r us i n(+)+dLus i n(L+)-g rs i n(+)-g dLs i n(L+)=-r z gc o s(+)-dLz gc o s(L+).(4 6)(1+)u+r c o s(+)-r2s i n(+)+dL c o s(L+)-dL2s i n(L+)=-(1+)z g-k E(1+)m g.(4 7)同前所述,考虑柱体几何特征可将上述控制方程简化为(4r2/3+d2L)+r u-r u(+)+dLu-dLu(L+)+r u(+)+dLu(L+)-g r(+)-g dL(L+)+(r+dL)z g=0.(4 8)(1+)u+r -r2(+)+dL-dL2(L+)=-(1+)z g-k E(1+)m g.(4 9)同理,若地震作用的初始方向与X轴负方向相同,则柱体摇摆、滑移的控制方程分别为(4r2/3+d2R)+r u-r u(-)+dRu-dRu(R-)+r u(-)+dRu(R-)+g r(-)+g dR(R-)+(r+dR)z g=0.(5 0)(1+)u+r -r2(-)+dR-dR2(R-)=-(1+)z g-k E(1+)m g.(5 1)3 初始动态响应影响因素分析不失一般性地,取水平地震加速度初始方向与X轴正方向相同的情况予以分析,则柱体可以被触发的响应临界条件分别为:静止(见式(9);抬升(见式(2 2);滑移及抬升滑移同步发生(见式(1 8)。比较各临界条件可得,柱体初始动态响应判定涉及基础平台顶面与柱体底面间静摩擦系数s、柱体长细比倒数t a n、上部结构与柱体自重比值及其偏心率e共4个参数。在不同动态响应判定条件中,的增减不显著影响动态响应判定结果是显而易见的,仅需对s、t a n和e 3个参数进行分析。由于中国木构传统建筑存世量大,上部结构的自重更易量化,故本节中的取值以相关研究2 0,2 1为基础,依据上部结构自重的概率分布情况取中值=7 61 9.式(1 8)表明,滑移临界加速度与静摩擦系数s成线性关系。依s的取值范围,在其取值范围内可绘制如图5(a)所示的s-z g曲线。当加速度值取值于临界值曲线上方区域时,柱体将发生滑移伴随摇摆响应。抬升响应临界条件中,水平加速度z g是柱体长细比t a n及柱顶构件偏心率e的二元函数。由于、e的取值范围较小,在0.1 00.1 5、-1.0e1.0范围内分别生成及e的随机数5 0 0个即可覆盖及e绝大部分取值的可能。将、e、代入式(6)可得抬升临界水平加速度,利用插值法绘制3维曲面如图5(b)所示。曲面之下的水平加速度不能触发柱体响应,柱体保持静止;曲面之上的水平加速度引起柱体摇摆,直至水平加速度突破图5(a)滑移临界曲线后,滑移将同步被触发。对比图5(a)、5(b)可知,浮放柱体最易被触发的响应是抬升,随即进入摇摆运动状态,即摇摆是柱体最主要的运动形式。由图5(b)可见,随e从-1至1逐渐增大,柱顶构件质心逐渐远离柱脚转动轴L点,激发柱体抬升所需力矩增大;当从0.1 0逐渐增大至0.1 5时,柱体体型由细长发展至短粗,自身稳定性增强,柱脚抬升所需的水平加速度同样增大。由图5(b)对比长细比-水平加速度、偏心率-水127 第4期 高 潮,等:水平地震作用下浮放柱体动态响应研究平加速度可知,后者的斜率大于前者,即抬升临界水平加速度对柱顶构件偏心率更敏感。0.450.550.400.550.600.650.700.750.80s0.800.750.700.650.600.550.500.450.40zg?/?g（a）滑移临界0.110.120.100.130.140.15-1-0.500.51e/?rad0.160.140.120.100.080.060.040.020zg?/?g（b）抬升临界图5 滑移及抬升临界水平加速度F i g.5 C r i t i c a l h o r i z o n t a l a c c e l e r a t i o n s o f s l i p a n d u p l i f t4 摇摆响应求解大量的地震波都包含长持时的脉冲作用,因而研究结构在正弦激励下的响应有助于掌握结构的动力学特性2 2。本节选取柱体在水平激励下最容易被触发的摇摆响应进行分析,求解其在正弦激励下的摇摆运动方程。作用于柱体的正弦激励为一个完整周期,形如z g(t)=aps i n(pt+),-/pt(2-)/p.(5 2)z g(t)=0,其他 (5 3)其中,t(t0)为自变量时间,ap、p、分别是激励的幅值、频率和初相。若激励开始作用于大地的瞬间(t=0)恰好可以引起柱脚抬升,则由式(2 1)、(3 4)可得aps i n=1+(1+e)1+2gt a n,(t)0.(5 5)将式(5 2)代入式(3 3)、(3 6),可得柱体摇摆角度为(t)=C1s i nh(LP t)+C2c o sh(LP t)+aps i n(pt+)g(2p/2LP2+1)-LdL/r+dL/r+1,(t)0.(5 7)其中,Ci(i=1,2,3,4)为积分常数,P是一个表征柱体动力性能的参数,单位为r a d/s;L、R为柱体几何和物理特征参数。上述参数表达式分别为P=3g/4r.(5 8)L=dL/r+13d2L/4r2+1.(5 9)R=dR/r+13d2R/4r2+1.(6 0)对式(5 6)、(5 7)求导可得柱体摇摆角速度(t)=C1LPc o sh(LP t)+C2LPs i nh(LP t)+appc o s(pt+)g(2p/2LP2+1),(t)0.(6 2)对式(6 1)、(6 2)求导可得柱体摇摆角加速度(t)=C12LP2s i nh(LP t)+C22LP2c o sh(LP t)-ap2ps i n(pt+)g(2p/2LP2+1),(t)0.(6 4)柱体初始为静止状态,在激励开始作用的瞬间抬升角度、角速度均为0,即t=0时满足(0)=0与(0)=0成立,分别代入式(5 6)、(5 7)、(6 1)、(6 2)可求得积分常数C1=-appc o sLP g(2p/2LP2+1).(6 5)C2=-aps i ng(2p/2LP2+1)+LdL/r+dL/r+1.(6 6)C3=-appc o sRP g(2p/2RP2+1).(6 7)C4=-aps i ng(2p/2RP2+1)-RdR/r+dR/r+1.(6 8)综上,柱体的摇摆角度、角速度、角加速度等运227太 原 理 工 大 学 学 报 第5 4卷 动特征受柱体几何特征、柱顶构件质量与其偏心程度共同影响。5 结论基于浮放传统建筑结构特征将柱体与柱顶构件简化为2维匀质刚体,考虑柱顶构件的偏心情况,研究了柱体在水平地震作用下的动态响应,得到如下结论:1)水平地震作用下,柱体可以被触发的动态响应为静止、抬升、先抬升后滑移及抬升滑移叠加,不可能发生的响应为无抬升仅滑移或先滑移后抬升。2)决定柱体动态响应类型的因素包含基础平面与柱体底面间静摩擦系数、柱体长细比、柱顶构件偏心率。3)水平地震作用可触发的各类柱体动态响应中,抬升最容易被触发,柱体继而进入以摇摆为主的运动。4)相较于柱体长细比,柱顶构件偏心率的改变能更显著影响柱体抬升响应的触发条件。5)柱体摇摆运动状态受柱体几何特征、柱顶构件质量与其偏心程度共同影响。6)考虑柱顶构件及其偏心情况的柱体动态响应研究是传统建筑预防性保护的基础,可为传统建筑状态评估提供理论支撑。参考文献:1 张润森,王琪.浮放物体平面多刚体动力学建模与算法研究J.力学学报,2 0 1 7,4 9(6):1 3 7 0-1 3 7 9.Z HAN G R S,WAN G Q.R e s e a r c h o n m o d e l i n g a n d n u m e r i c a l m e t h o d o f f r e e s t a n d i n g b o d y o n p l a n a r r i g i d m u l t i b o d y d y n a m-i c sJ.C h i n e s e J o u r n a l o f T h e o r e t i c a l a n d A p p l i e d M e c h a n i c s,2 0 1 7,4 9(6):1 3 7 0-1 3 7 9.2 GAO C,HE J X.D y n a m i c r e s p o n s e o f w o o d e n c o l u m n s i n t r a d i t i o n a l t i m b e r s t r u c t u r e s u n d e r h o r i z o n t a l e a r t h q u a k e