-
引言
-
建筑结构在地震、风等动力荷载作用下,会产生较大的振动响应,导致结构的变形和破坏[1].随着现代制造技术的快速发展,建筑结构更加轻柔化和复杂化,其振动现象日益显著.例如,摩天大厦中的悬挑、不规则等设计会造成结构刚度突变,诱发大幅的持续振动[2].为了有效抑制结构的振动水平,人们提出了多种控制方法,如非线性吸振器等.
-
非线性吸振器能够有效拓展减振带宽和提升鲁棒性.目前,非线性吸振器的种类和形式不断更新,如分段线性式[3-4]、碰撞式[5-6]、轨道式[7-9]等.轨道式吸振器通过设计轨道形状提供不同形式的非线性作用力,并在土木工程中得到了应用.例如,轨道吸振器可以将输入建筑结构的能量从低阶模态向高阶模态转化,从而快速有效地衰减结构响应[7].Lu等人提出了一种轨道非线性能量汇,给出了与质量比和轨道形状相关的优化设计过程,并在五层钢框架结构上进行了一系列振动试验,有效验证了轨道吸振器的良好减振性能[8].此外,单边碰撞的轨道非线性吸振器不仅能够吸收和耗散建筑结构的能量,还可以在结构刚度和输入能级变化时表现出很强的鲁棒性[9].
-
近年来,主动控制备受关注,能够提升控制效果和适应复杂工况.在控制过程中,传感器、滤波器、控制器、作动器等环节中存在时滞效应.Olgac首次提出了时滞动力吸振器的概念[10].徐鉴教授团队长期从事时滞吸振器的研究,取得了许多重要进展,并得到了实际应用[11-13].Zheng等人研究了具有时滞加速度反馈的主动质量阻尼器对建筑结构的振动控制问题,给出了控制增益和时滞的参数空间中的稳定区域[14].Bellizzi等人考察了时滞主动阻尼非线性吸振器,基于复变量平均法和几何奇异摄动法预测了周期响应和准周期响应.结果表明,较小的时滞会导致系统状态发生重大变化,如抑制准周期响应区域和降低振动水平[15].Mohanty等人提出了具有时滞反馈的压电片堆式主动非线性吸振器,给出了基础简谐激励和外扰力共同作用下的振动控制效果,可有效提升安全性和经济性[16].
-
本文基于轨道型非线性吸振器,引入时滞位移反馈控制,提出一种时滞轨道非线性吸振器,并应用于建筑结构的振动抑制.
-
1 动力学模型
-
如图1所示,建筑结构质量、刚度、阻尼系数分别为m1,k1,c1,其顶部置放吸振器.吸振器由带有轨道的小车和小球构成,其中m2为小车质量,m3为小球质量,c3为小球与小车间的阻尼系数.轨道小球视为刚体,可沿轨道自由滑动,其轨道形状满足y3(t)=a[x3(t)-x2(t)]4,其中y3为小球的竖向位移,x2和x3分别为小车和小球的水平位移,a为轨道形状系数[8].建筑结构和吸振器间的弹簧刚度系数为k2,阻尼系数为c2.采用时滞位移反馈控制,控制力取为Ax2(t-τ),其中A为反馈增益系数,τ为时滞量.
-
图1 含时滞轨道非线性吸振器的建筑结构
-
Fig.1 A building structure with a nonlinear delayed vibration absorber
-
系统的总动能T表达式为:
-
其中x1为框架结构的水平位移,v3为轨道小球的速度,满足,其中和分别为小球的水平和竖直方向速度.系统的总势能V表达式为:
-
根据拉格朗日方程,则有:
-
其中,Qi为外激励和控制力;D为耗散函数[17],即
-
将式(1)和式(2)代入式(3)中,整理后可得:
-
其中f(t)=f0cos(ωt),f0为激励幅值,ω为激励频率.小车和小球之间的非线性力为FT=FT1+FT2+FT3,其中,g=9.8N/kg为重力加速度.
-
系统参数取为:m1=30kg,k1=2×104N/m,c1=0.05Ns/m,m2=3kg,k2=3470N/m,c2=0.6Ns/m,m3=3kg,a=200m-3,c3=0.8Ns/m[7].激励幅值f0分别取为1N和10N,当激励频率ω=23.5rad/s时计算在t=70T内系统响应,其中T为外激励周期.图2为FT,FT1,FT2以及FT3的变化曲线.如图2(a)所示,FT的主要部分为FT1,而FT2和FT3在FT中占比很小.如图2(b)所示,FT2和FT3的影响较大,在FT1,FT2以及FT3的同时作用下,FT在出现了峰值现象,即当相对位移增大到某个值后,FT会随其的增大而减小,这种现象有助于在发生较大响应时减轻由于作用力过大而造成减振设备的损害[7].
-
图2 不同外激励幅值作用下非线性力的变化曲线
-
Fig.2 The curve of nonlinear force for different excitation amplitudes
-
2 幅频特性
-
本文采用谐波平衡法(HBM)求解系统动力学响应[16].为简化起见,引入轨道小球相对轨道小车的位移xn(t)=x3(t)-x2(t).方程可改写为:
-
其中.式(5)的谐波近似解设为如下形式:
-
通过将式(6)代入式(5)中计算残差力,并借助伽辽金条件可得到非线性微分方程组.由于a1(t),b1(t),a2(t),b2(t),a3(t),b3(t)均是关于时间t的慢变函数,忽略高阶项,从而得到关于和的微分方程组.其中,S=a1(t),b1(t),a2(t),b2(t),a3(t),b3(t)T,矩阵N=(npq),Z={z1,z2,z3,z4,z5,z6}T,npq和zq表达式见附录,p,q=1,2,···,6.根据定常解条件=0,可得非线性代数方程组N-1Z=0.借助弧长延拓法求解该代数方程组,可给出位移幅值.通过计算雅可比矩阵的特征值可判断定常解的稳定性.若所有特征值实部为负则稳态解稳定,为正则不稳定.
-
图3 主系统的位移幅频曲线
-
Fig.3 The amplitude-frequency curve of the primary system
-
图3为当k2=3470N/m,a=200m-3,f0=10N,A=-80N/m,τ=0.05s时,建筑结构的位移的幅频曲线.为了验证谐波平衡法的准确性,采用龙格库塔法(RKM)进行验证.由图3可知,龙格库塔法与谐波平衡法的结果在稳定周期解的区域基本吻合,说明两种方法的结果相一致.此外,主系统位移幅频曲线的共振峰向左偏转,表现出软弹簧特性.同时,在共振峰附近出现了不稳定的解.
-
3 参数讨论
-
采用纯被动控制时,考察吸振器参数对于幅频曲线的影响情况.如图4(a)所示,随着弹簧刚度系数由3700N/m增大到5200N/m,幅频响应曲线向右移动,并且向左偏转的程度有所减弱.此外,在稳定区域内,周期运动的位移振幅的峰值从0.0535m增大到0.0763m,而不稳定区域内的共振峰幅值从0.2039m增大到0.5665m.图4(b)给出了不同轨道形状系数对应的主系统的幅频曲线.当时,吸振器退化为线性情形.当时,系统的幅频曲线向左偏转,共振峰附近也出现了不稳定解.随着轨道形状系数的增大,幅频曲线向左偏转的程度也增大.
-
图4 不同弹簧刚度系数k2与轨道形状系数a时,主系统的位移幅频曲线
-
Fig.4 Effects of the spring stiffness k2 and the coefficient of track shape a on the amplitude-frequency curve of the primary system
-
下面考察位移反馈控制,首先考虑无时滞情形.系统的弹簧刚度系数k2与轨道形状系数a分别取为3470N/m和200m-3,其他参数保持不变.如图5所示,反馈增益的值从零变为-500N/m时,稳定周期运动振幅的峰值分别为:0.0477m,0.0431m,0.0379m,0.0387m,0.0347m.不稳定的共振峰幅值从0.1658m降低到0.1000m.与无反馈控制相比,受控系统的稳定周期运动的幅值的峰值与不稳定区域的共振峰幅值都有所降低.
-
图5 不同反馈增益A 的主系统位移幅频曲线
-
Fig.5 Effects of the feedback gain A on the amplitude-frequency curve of the primary system
-
图6 不同时滞下的主系统位移幅频曲线
-
Fig.6 Effects of the time delay on the amplitude-frequency curve of the primary system
-
考虑控制中的时滞效应,而其余系统参数保持不变.图6给出了不同时滞量情形下主系统的位移幅频曲线,其中.如图6所示,无时滞时,幅频曲线的共振峰表现出向左偏转特性,并且在共振峰附近出现不稳定区域.当τ=0.01s时,共振峰基本不偏转,不稳定区域消失,类似于线性情形.此时,共振峰的幅值降至0.0230m.随着时滞量的继续增大,幅频曲线的峰值进一步降低,并与龙格库塔法的计算结果基本吻合.这说明,引入时滞可以消除幅频响应曲线中的不稳定区域,减弱系统的软弹簧特性,并可提高振动抑制效果.
-
图7 主系统位移随时滞量变化的分叉图
-
Fig.7 The bifurcation diagram of the displacement of the primary structure with the change of time delay
-
图8 在主结构位移随外激励幅值变化的分叉图
-
Fig.8 The bifurcation diagram of the displacement of the primary system as the excitation amplitude varies
-
图7给出了不同反馈增益时,主结构的位移在随时滞量变化的分叉图,其中外激励频率确取为23.1rad/s.Poincaré截面选取为P={(τ,x1(t)):(y1(t)=0,(t)>0)},其中y1(t)= (t).如图7所示,时滞量位于0.002s和0.08s之间,主系统的响应为周期运动,而继续增加时滞量会导致混沌运动.随着反馈增益绝对值的增加,当0.002s<τ<0.012s时,系统周期运动的振幅不断增加,而当0.033s<τ<0.08s时,系统周期运动的振幅减小,进一步抑制了振动.这表明,反馈增益和时滞量均会影响周期运动的幅值.此外,当时滞量很小时,即约为0~0.002s,随着反馈增益绝对值的增大,周期振荡可变为混沌运动.
-
图8给出了主结构的位移在Poincaré截面上随激励幅值变化的分叉图,其中Poincaré截面为P={(f0,x1(t)):(y1(t)=0,(t)>0)},外激励频率为23.1rad/s.如图8所示,当A=-500N/m,τ=0时,即无时滞的位移反馈控制,在f0>6.36N情形下主系统会进入混沌状态,而时滞取0.03s可使系统由混沌变为周期运动.保持时滞量不变,反馈增益分别取-400N/m,-600N/m,-700N/m时,主系统的位移幅值会随着激励幅值的增大而增加.这说明,在激励幅值较大时,时滞能够抑制混沌响应,改善减振效果.并具有一定的鲁棒性.
-
4 结论
-
本文提出了一种时滞轨道非线性吸振器,并将其安装在建筑结构顶部.通过谐波平衡法结合弧长延拓法给出了幅频响应曲线,并与龙格库塔法的结果进行了对比验证.结果表明,系统的幅频响应曲线共振峰会出现向左偏转现象,并且共振峰附近存在不稳定区域.在被动控制时,轨道非线性造成了幅频曲线共振峰附近的软弹簧特性与不稳定性.考虑时滞反馈控制后,可消除幅频曲线中的不稳定区域,降低甚至消除共振峰的软弹簧特性,并能够降低共振峰的幅值.此外,时滞反馈控制还可有效调控复杂运动,改善振动抑制效果.
-
附录:
-
参考文献
-
[1] THENOZHI S,YU W.Advances in modeling and vibration control of building structures [J].Annual Reviews in Control,2013,37(2):346-364.
-
[2] 李子懿,肖从真,李建辉,等.复杂高层建筑结构抗震设计方法研究现状与展望 [J].建筑科学,2019,35(9):120-125.LI Z Y,XIAO C Z,LI J H,et al.Research status and prospects of seismic design method for complicated high-rise building structure [J].Building Science,2019,35(9):120-125.(in Chinese)
-
[3] GOURDON E,LAMARQUE C H.Energy pumping with various nonlinear structures:Numerical evidences [J].Nonlinear Dynamics,2005,40(3):281-307.
-
[4] 姚红良,张钦,杨沛然,等.分段线性刚度非线性能量阱的参数优化方法 [J].东北大学学报(自然科学版),2019,40(12):1732-1738.YAO H L,ZHANG Q,YANG P R,et al.Optimization method of nonlinear energy sinks with piecewise linear stiffness [J].Journal of Northeastern University(Natural Science),2019,40(12):1732-1738.(in Chinese)
-
[5] GENDELMAN O.Analytic treatment of a system with a vibro-impact nonlinear energy sink [J].Journal of Sound and Vibration,2012,331(21):4599-4608.
-
[6] LI T,SEGUY S,BERLIOZ A.Optimization mechanism of targeted energy transfer with vibro-impact energy sink under periodic and transient excitation [J].Nonlinear Dynamics,2017,87(4):2415-2433.
-
[7] WANG J,WIERSCHEM N,SPENCER J B F,et al.Experimental study of track nonlinear energy sinks for dynamic response reduction [J].Engineering Structures,2015,94:9-15.
-
[8] LU X,LIU Z,LU Z.Optimization design and experimental verification of track nonlinear energy sink for vibration control under seismic excitation [J].Structural Control and Health Monitoring,2017,24(12):e2033.
-
[9] LI W,WIERSCHEM N E,LI X,et al.Numerical study of a single-sided vibro-impact track nonlinear energy sink considering horizontal and vertical dynamics [J].Journal of Vibration and Acoustics,2019,141(6):061013.
-
[10] OLGAC N,HOLM-HANSEN B T.A Novel active vibration absorption technique:delayed resonator [J].Journal of Sound and Vibration,1994,176(1):93-104.
-
[11] 赵艳影,徐鉴.时滞非线性动力吸振器的减振机理 [J].力学学报,2008,40(1):98-106.ZHAO Y Y,XU J.Mechanism analysis of delayed nonlinear vibration absorber [J].Chinese Journal of Theoretical and Applied Mechanics,2008,40(1):98-106.(in Chinese)
-
[12] SUN Y,XU J.Experiments and analysis for a controlled mechanical absorber considering delay effect [J].Journal of Sound and Vibration,2015,339:25-37.
-
[13] MENG H,SUN X,XU J,et al.The generalization of equal-peak method for delay-coupled nonlinear system [J].Physica D:Nonlinear Phenomena,2020,403:132340.
-
[14] ZHENG Y,HUANG J,SUN Y,et al.Building vibration control by active mass damper with delayed acceleration feedback:multi-objective optimal design and experimental validation [J].Journal of Vibration and Acoustics,2018,140(4):1-7.
-
[15] BELLIZZI S,CHUNG K-W,SAMPAIO R.Response regimes of a linear oscillator with a nonlinear energy sink involving an active damper with delay [J].Nonlinear Dynamics,2019,97(2):1667-1684.
-
[16] MOHANTY S,DWIVEDY S K.Nonlinear dynamics of piezoelectric-based active nonlinear vibration absorber using time delay acceleration feedback [J].Nonlinear Dynamics,2019,98(2):1465-1490.
-
[17] 陈峻峰,刘昆,梁文杰,等.磁悬浮飞轮储能系统机电耦合非线性动力学研究 [J].动力学与控制学报,2013,11(3):225-234.CHEN J,LIU K,LIANG W,et al.Study on nonlinear dynamics of electromechanical coupling in flywheel energy storage system based on active magnetic bearings [J].Journal of Dynamics and Control,2013,11(3):225-234.(in Chinese)
-
摘要
本文研究了具有时滞轨道非线性吸振器的建筑结构的动力响应及其控制效果.采用谐波平衡法和弧长延拓法给出了简谐激励作用下主结构的幅频响应曲线,并与龙格库塔法的结果进行了对比.考察了不同参数情况下幅频曲线的变化情况,揭示了系统的复杂运动现象.结果表明,轨道非线性会导致幅频曲线向左偏转.时滞反馈控制能够降低主结构的位移幅值,并可抑制混沌响应.
Abstract
In this paper, the dynamical responses and control performance of a building structure with a delayed track nonlinear vibration absorber are studied. The amplitude-frequency response curves of the system under external excitation are given by harmonic balance method and arc-length continuation method. The fourth-order Runge-kutta algorithm is performed to validate the obtained results. The properties of the amplitude-frequency response for different system parameters are analyzed and complex dynamical motions are shown. It is shown that the track nonlinearity can lead to the left-leaning phenomenon. In addition, the delayed feedback control can reduce the displacement amplitude of the primary structure and suppress the complicated responses.