拉索-阻尼器系统的改进实模态分析方法_陈志.pdf

第2 1卷第4期2 0 2 3年4月动 力 学 与 控 制 学 报J OUR NA LO FD YNAM I C SAN DC ON T R O LV o l.2 1N o.4A p r.2 0 2 3文章编号:1 6 7 2-6 5 5 3-2 0 2 3-2 1(4)-0 0 7-0 0 8D O I:1 0.6 0 5 2/1 6 7 2-6 5 5 3-2 0 2 3-0 3 8 2 0 2 2-1 1-0 9收到第1稿,2 0 2 3-0 3-2 2收到修改稿.*国家自然科学基金资助项目(5 2 2 7 8 3 0 7),湖南省自然科学基金资助项目(2 0 2 2 J J 3 0 2 6 1),N a t i o n a lN a t u r a lS c i e n c eF o u n d a t i o no fC h i n a(5 2 2 7 8 3 0 7)a n dt h eN a t i o n a lN a t u r a lS c i e n c eF o u n d a t i o no fH u n a nP r o v i n c e(2 0 2 2 J J 3 0 2 6 1).通信作者E-m a i l:c e_h a n y a n 1 6 3.c o m拉索-阻尼器系统的改进实模态分析方法*陈志1 彭文林2 禹见达1 韩艳2 田静莹2(1.湖南科技大学 土木工程学院,湘潭 4 1 1 2 0 1)(2.长沙理工大学 土木工程学院,长沙 4 1 0 1 1 4)摘要 阻尼减振是拉索振动控制的主要措施,目前拉索-阻尼器系统的主要分析方法为复模态法.虽然采用复模态法可以精确分析拉索-阻尼器系统的阻尼比和频率等主要特征参数,但采用复模态法求出的拉索复振型由于不包含减振所需要的参数而被研究较少,而系统的振型对研究阻尼器的减振效果又具有非常重要的意义.文中首先采用复模态法求解了拉索-阻尼器系统,基于复振型的物理意义,引入相位函数,提出了一种改进的实数分离变量法,分析了阻尼力在拉索中的传递特征;建立了拉索-阻尼器系统的改进实模态分析方法,采用该方法求解得到的非线性运动方程通解与复模态法一致.基于该新方法计算得到的实模态振型,分析了阻尼器对拉索实际振型以及相位的影响.关键词 拉索-阻尼器系统,复模态,分离变量法,相位差中图分类号:U 4 4 1+.3文献标志码:AA nI m p r o v e dM e t h o dt oR e a lM o d a lA n a l y s i so fC a b l e-D a m p e rS y s t e m s*C h e nZ h i1 P e n gW e n l i n2 Y uJ i a n d a1 H a nY a n2 T i a nJ i n g y i n g2(1.S c h o o l o fC i v i lE n g i n e e r i n g,H u n a nU n i v e r s i t yo fS c i e n c ea n dT e c h n o l o g y,X i a n g t a n 4 1 1 2 0 1,C h i n a)(2.S c h o o l o fC i v i lE n g i n e e r i n g,C h a n g s h aU n i v e r s i t yo fS c i e n c ea n dT e c h n o l o g y,C h a n g s h a 4 1 0 1 1 4,C h i n a)A b s t r a c t D a m p e ra r e t h em a i nm e a s u r ef o rc a b l ev i b r a t i o nc o n t r o l,a n dc o m p l e xm o d a la n a l y s i s i st h em a i na p p r o a c ho f c a b l e-d a m p e r s y s t e ma n a l y s i s.T h e c o m p l e x f r e q u e n c i e sd e r i v e d f r o mc o m p l e xm o d a l a-n a l y s i sa l l o wa na c c u r a t ea n a l y s i so f s y s t e msd a m p i n gr a t i oa n d f r e q u e n c y,w h i c h i so fm a j o r i n t e r e s t o fb r i d g er e s e a r c h e r s.T h e r e i s f e ws t u d i e s f o r i n v e s t i g a t i n gt h ec o m p l e xm o d es h a p eo fc a b l es i n c e i td o e sn o t c o n t a i nt h ek e yp a r a m e t e r sf o rv i b r a t i o ns u p p r e s s i o n.T h i sp a p e rf o c u s e so nt h ea n a l y s i so fc a b l esp h a s ec o n t a i n e d i nt h ec o m p l e xm o d a l s h a p e.F i r s t l y,t h ec o m p l e xm o d a l a n a l y s i s i su s e dt oa n a l y s i s t h ec a b l e-d a m p e rs y s t e m,a n dt h ep h y s i c a lm e a n i n go fc o m p l e x m o d es h a p ei se x p l a i n e d.S u b s e q u e n t l y,b a s e do n i t sp h y s i c a lm e a n i n g s,t h er e a ls e p a r a t i n gv a r i a b l em e t h o di s i m p r o v e db ya d d i n gp h a s ef u n c-t i o n,t h u s i t i sp o s s i b l e t oe x p l a i nh o wt h ed a m p i n g f o r c e i s t r a n s l a t e d i nt h e c a b l e.I nt h i sp r o c e s s,as e to fn o n l i n e a re q u a t i o n sa r eo b t a i n e da n dt h e i rg e n e r a ls o l u t i o n sa r es o l v e d,a n dt h er e s u l t so b t a i n e da r ec o n s i s t e n tw i t h t h o s eb y t h e c o m p l e xm o d a l a n a l y s i s.F i n a l l y,t h e e f f e c t o f t h ed a m p e r o n t h e c a b l es r e a lm o d a l s h a p ea n dp h a s e i sa n a l y z e db a s e do nt h er e s u l t so b t a i n e d.K e yw o r d s c a b l e-d a m p e rs y s t e m,c o m p l e xm o d a l,s e p a r a t ev a r i a b l em e t h o d,p h a s ed e v i a t i o n动 力 学 与 控 制 学 报2 0 2 3年第2 1卷引言作为斜拉桥的关键承重构件,大跨度斜拉索极易在环境荷载下诱发大幅振动1-3,如不加以抑制,则会导致拉索锚固端疲劳甚至破坏.目前拉索控制振动的主要方法有气动措施、结构措施以及阻尼耗能措施4-6,其中安装阻尼器的耗能减振方法对不同机理和模态的振动均有抑制效果,因此应用较广.阻尼减振技术在拉索上应用至今,目前已有许多阻尼器设计方法,最早的拉索阻尼器设计方式可追溯至K o v a c s7,但该方法只能考虑一阶模态的最优值.P a c h e c o等8采用G a l e r k i n方法分析了张紧弦模型,得到了安装有阻尼器的拉索模态阻尼比,并对阻尼参数无量纲化后得出了拉索各模态通用的阻尼器设计曲线.P a c h e c o的设计曲线方便易用,但由于G a l e r k i n法不能考虑阻尼器对拉索模态的影响,并且忽略了垂度、拉索抗弯刚度等因素,因此计算结果偏大,有试验证明实测拉索最大阻尼比只达到了理论上限的6 0%9.拉索-阻尼器系统具有较强的非比例阻尼特性,这是由阻尼器安装的不均匀性引起的.复模态理论1 0可以精确分析非比例阻尼系统,K r e n k1 1在拉索-阻尼器系统中引入该理论,假设拉索振型为复模态振型,并认为复模态振型在阻尼器安装位置处不连续,以此提出了拉索-阻尼器系统的复模态分析方法,并在后续考虑了拉索垂度影响1 2.M a i n和J o n e s1 3-1 5根据复模态方法,重点研究了阻尼器对拉索结构模态频率的影响,发现相比于线性阻尼,非线性阻尼器可以在更大的模态范围内实现最优的阻尼性能.K r e n k和H g s b e r g1 6采用复模态法分析了阻尼器并联刚度和集中质量的影响,发现集中质量对阻尼器的效果有增强作用,即“负刚度增强效应”.近年来,一些被动负刚度装置以及惯性装置陆续出现1 7-1 8,为了分析这些装置对拉索阻尼器的增强效应,复模态法被广泛采用.S u n等1 9应用复模态理论研究了两根拉索之间连接调谐惯质阻尼器后的振动特性.W a n g等2 0基于复模态方法从理论上研究了惯性质量对垂度拉索阻尼器效果的增强效应,并设置了模型试验进行验证.S h i等2 1用复模态法分析了调谐惯质阻尼器对拉索的减振效果.孙浚杰等2 2采用复模态方法理论研究了惯性质量对拉索线性黏滞阻尼器减振效果的增强效应.刘菁等2 3使用复模态分析,对某实际拉索进行了惯质阻尼器参数设计,并分析了减振效果.程志鹏等2 4基于能量等效线性化方法和复模态法,分析了负刚度非线性黏滞阻尼器对斜拉索减振的增强效应.X u等2 5采用磁流变阻尼器与可编程的智能控制装置实现了伪被动负刚度,并在斜拉索的理论分析中,使用复模态考察了阻尼系数和负刚度系数对斜拉索振动主动控制效果的影响.除此以外,若要精确考虑集中刚度、阻尼等因素对拉索模 态的影响,也 需要用到复 模态.N a mH o a n g等2 6应用复模态分析了弯曲刚度对拉索阻尼器的影响;J a v a i dA h m a d等6建立了双索之间连接连杆件的双索网络模型,采用复模态理论分析了拉索长度比、连杆的位置、拉索频率比和拉索质量-张力比对双索网络模型的振动特性影响;杨超等2 7现场观测了超长拉索的高阶模态振动现象,随后根据复模态法研究了适用于拉索多模态振动的阻尼器参数方案.孙利民等2 8采用复模态法分析了垂度拉索在安装双阻尼器后的效果,并在实桥现场试验验证了理论计算精确性.Y o u n e s p o u r等2 9采用复模态方法分析了垂度拉索在安装了横向弹性支承后的振型与频率变化,讨论了支撑刚度、支撑位置以及拉索的弯曲刚度等参数的影响.对于安装了阻尼器的拉索,桥梁研究者更关注拉索的阻尼比以及振动频率,多采用复模态法分析复频率,没有过多深入研究复模态振型,但其对研究阻尼器的减振效果又不容忽视.采用常规复模态法得到的复振型可以表示拉索的振幅与相位,并且其中相位与拉索各点的运动衰减息息相关,但由于是复数分析,其物理含义不直观,所以没有办法研究阻尼器对振型的影响.本文基于复振型与分离变量法提出了改进实模态分析方法,该方法可以清晰地描述拉索-阻尼器系统运动衰减与相位的相关性,也能自洽地解释阻尼器对拉索的振型影响.本文主要结构如下:第一部分为拉索-阻尼器系统的常规分析方法;第二部分为改进的实模态分析方法与求解;第三部分为阻尼器对拉索振型与相位影响的实例分析;最后在第四部分对全文进行了总结.1 拉索-阻尼器系统的常规分析方法1.1 实数分离变量法本文采用均匀张紧弦模型表示拉索,忽略拉索8第4期陈志等:拉索-阻尼器系统的改进实模态分析方法本身的阻尼以及重力影响,如图1所示,引入函数表示集中力,则拉索自由运动微分方程可写为:2u(x,t)t2+cu(l1,t)t(x-l1)=T2u(x,t)x2(1)式中为拉索的每延米质量,T为拉索索力,L为拉索长度,l1为阻尼器的安装位置,u(x,t)表示t时刻x位置处斜拉索的振动位移,c为阻尼器的阻尼黏性系数,为狄拉克函数.图1 张紧弦-阻尼器模型F i g.1 T a u t c a b l ew i t he x t e r n a l v i s c o u sd a m p e r除去阻尼器安装点和两端点以外,拉索不再受外力,因此可将拉索按两端点和阻尼安装点分为两段,这两段拉索必定满足无外加激励下弦振动方程:2u(x,t)t2-T2u(x,t)x2=0(2)对方程(2)的解作分离变量法假设u(x,t)=Y(x)F(t),并且考虑实际物理意义,认为Y(x),F(t)都为实数,代入方程(2)得到:d2F(t)dt2+KF(t)=0(3)式中K=-Td2Y(x)dx2/Y(x),式(3)可视为无阻尼振子的自由运动方程,其求得的时间函数必然不会随时间衰减,这与实际事实相违背,因此传统的实数分离变量法不再适用于拉索-阻尼器系统求解.问题也由此提出:拉索-阻尼器系统的位移应当如何假设与求解?并且引申出以下问题:拉索上安装的阻尼器怎样影响其他无外力作用的点自由衰减?即阻尼器产生的阻尼力如何在拉索中传递.1.2 复模态法有关拉索-阻尼器系统的精确求解,K r e n k早在文献1 1 给出了解答,后续许多拉索系统研究都是基于此方法.下面应用复模态法对拉索-阻尼器系统进行分析,由上一部分可知,拉索运动无法分离成实数空间函数与实数时间函数乘积的形式,既然如此,不妨假设它们均为复数,即:u(x,t)=Y(x)eit(4)式中=|(1-2+i)为复频率,Y(x)为复振型,为系统阻尼比.将式(4)代入式(2)得到:d2Y(x)dx2+2Y(x)=0(5)式中=/T为复波数,方程通解为:Y(x)=A1s i nx+A2c o sx(6)A1、A2为待定系数,需要根据边界条件确定,阻尼器安装点由于阻尼力的存在,拉索复振型导数并不连续,因此需要分为两段,如图2所示,右下标为1代表左边索段,为2代表右边索段.假索的振动幅度较小,则索轴线与水平轴的夹角为u/x,因此得到这两段索的边界条件与阻尼器安装点的平衡条件:u1(0,t)=0u2(0,t)=0u1(l1,t)=u2(l2,t)Tu1(l1,t)x+u2(l2,t)x+cu(l1,t)t=0(7)复振型求得:Y(x)=Y0s i nx1s i nl1(0 x1l1)Y0s i nx2s i nl2(0 x2l2)(8)图2 拉索分段与安装点受力平衡F i g.2 B a l a n c eo f f o r c e sa tm o u n t i n gp o i n t s式中Y0为阻尼器安装点的位移复幅值,需要根据拉索自由振动的初始状态求解,再将式(8)代入,安装点的平衡条件得到:c o tl1L+c o t(l2)L+ic Tm=0(9)超越方程(9)称为波数方程,根据方程(9)可以求解出复波数,从而实现复模态方法的闭合分析.桥梁研究者主要关心拉索的阻尼比以及振动频率,9动 力 学 与 控 制 学 报2 0 2 3年第2 1卷因此只需要求解方程(9)就可以分析出阻尼器对拉索的阻尼比和频率影响.但由于复数的物理意义并不明确,对于阻尼器产生的阻尼力如何在拉索中传递这一点,只分析波数方程无法给出清晰的解释.2 改进实模态求解方法2.1 基于改进实数分离变量法的系统建模由上一部分内容可知,复模态解u(x,t)能够满足方程(2)及其边界条件,那么其实部与虚部也会满足要求,我们对式(4)作稍许变化:u(x,t)=Y(x)eit=|Y(x)|eia r g(x)ieit=|Y(x)|e-|tei t+a r g(x)(1 0)式中|为取模运算,a r g(x)为复振型Y(x)的辐角,=|1-2为自由振动频率.应用欧拉公式,可以发现式(1 0)的实部或者虚部就是随着时间变化振幅指数衰减的谐波,除了各点的运动幅值不同外,各点的运动相位也不同.模为拉索实际振幅,辐角为拉索各点的运动相位3 0.由复振型的物理意义,我们对实数分离变量法进行改进,在时间函数内添加一个空间相关的相位函数(x):u(x,t)=Y(x)F(t+(x)(1 1)此时式(1 1)也一定能够适用于拉索-阻尼器系统求解,令=t+(x),注意:2u(x,t)t2=Y(x)d2F()d2(1 2 a)2ux2=d2Y(x)dx2F()+2dY(x)dxdF()dd(x)dx+Y(x)d2F()d2d(x)dx2+dF()dd2(x)dx2(1 2 b)代入方程(2)并整理得到:M-d2F()d2+C-dF()d+K-F()=0(1 3 a)M-=-Td(x)dx2C-=-T2dY(x)dxd(x)dx/Y(x)+d2(x)dx2K-=-Td2Y(x)dx2/Y(1 3 b)由方程(1 3)可以发现,对于系统上无外力段的任意一点x,都可视为有阻尼振子,为振子的广义位移,当然其前提条件是d(x)/dx0,即(x)常数,若令d(x)/dx=0,则方程(1 3)退化为式(3).图3 拉索微元及等效模型F i g.3 D i f f e r e n t i a l u n i to f c a b l ea n d i t se q u i v a l e n tm o d e l提取某x0处拉索微元如图3所示,对于拉索上任一不受外力的微元,其上索力横向不平衡量与惯性力相等,如果微元两端无相位差(d(x)/dx=0),索力横向不平衡量只提供刚度回复力,可等效为无 阻 尼 振 子(图3(b);如 果 存 在 相 位 差(d(x)/dx0),索力横向不平衡量不仅提供刚度恢复力,还提供阻尼力与惯性力,等效模型如图3(c).由此可见相位差d(x)/dx在阻尼器阻尼力传递中扮演着重要的角色,并且会对微元的惯性力和刚度回复力产生影响,间接导致微元振动幅值的变化,从而影响拉索振型.2.2 基于改进实数分离变量法的实模态求解方法根据复模态法,拉索-阻尼器系统的自由振动为:各点阻尼比、频率一致、相位差不同的周期衰减谐波,因此将方程(1 3)中的F()设为:F(g)=e-s i n()(1 4)为了后续书写方便,用对数减幅表示衰减,仍然为自由振动圆频率.将其代入式(1 3)并考虑三角函数的系数为0,得到如下方程组:Y+(2+2)2Y=22 Y=-+22-12Y(1 5)式中,分别表示对x求一阶导数和二阶导数,=T/.令Y(x)=y(x)e(x),可使方程变为:01第4期陈志等:拉索-阻尼器系统的改进实模态分析方法y+2-22y=22y(1 6 a)2 y+22y=0(1 6 b)式中y(x)即为实际的振幅,式(1 6)是一组非线性的常微分方程,若要求其通解,需要做一些变化,令=f(x)g(x),那么=f(x)g(x)-g(x)g(x),观察式(1 6 b),发现若:f(x)g(x)=-22(1 7 a)式(1 6 b)则变为:g(x)g(x)=2y y(1 7 b)因此可得到:g(x)=y2(x),f(x)=-22y2(x)dx,令(x)=y2(x)dx,可以推出(x)以及y(x)与(x)关系:(x)=-22(x)(x)(x)=y2(x)=2y y (x)=2y y+2y 2(1 8)再将式(1 6 a)变化为关于(x)的非线性常微分方程:2 -2+42-22 2-1 62242=0(1 9)方程(1 9)不难求出 特解,只需要 令(x)=es x,可得特征方程:s2-22s2+22=0(2 0)(x)的四个特解,即:1(x)=e2x,2(x)=e-2x,3(x)=ei2 x,4(x)=e-i2 x(2 1)一般来说,非线性方程通解无法表示为特解的组合叠加,但这四个特解存在一种特殊的组合使其能够表达方程(1 9)的通解,将其变化为三角函数和双曲函数并整理得到如下通解:(x)=A2c o s2 x+B2s i n2 x+D2s i n h2x(2 2)式中s i n h以及下文的c o s h、t a n h分别为双曲正弦、双曲余弦和双曲正切函数;A、B、D为待确定实常数,并且必须满足A2+B2=D2,具体的取值需要根据边界条件和某一点初始振幅确定.实振型函数y(x)与相位差d(x)dx由(x)表示:y(x)=(x),d(x)dx=-22(x)(x)(2 3)由d(x)dx由(x)的 原函数得到 相位函数(x):(x)=1a r c t a nt a n x+12a r c t a nABt a n hx+n+C(2 4)式中C为待定常数,需要根据某一点初始相位确定,n=0,1,2,保证函数连续.拉索各点的位移因此表示为:u(x,t)=y(x)e-ts i n(t+(x)(2 5)式(2 5)中有A、B、C、D四个待定参数,两段不受外力的拉索都可用此结果表示,因此有8个待定系数以及两个假设的系统参数、,一共1 0个未知量,在已知边界条件(2 6 a)、平衡条件(2 6 b)、阻尼器安装点初始状态(2 6 c)和系数关系(2 6 d)的情况下,可得到如下1 0个方程:y1(0)=0,y2(0)=0y1(l1)=y2(L-l1)1(l1)=2(L-l1)(2 6 a)dy1(l1)dx+dy2(L-l1)dx=cTy1(l1)d1(l1)dx+d2(L-l1)dx+cT=0(2 6 b)y1(l1)=y0,1(l1)=0(2 6 c)A21+B21=D21,A22+B22=D22(2 6 d)式中下标为1表示阻尼器左边段拉索,下标为2表示阻尼器右边段拉索,Y0为阻尼器安装点的初始幅值,0为阻尼器安装点的初始相位.通过式(2 6)中包含的1 0个方程,最终可以得出8个待定系数的解析值,以及两个超越方程:s i n2 x1Lc o s h2x2L+s i n2 x2Lc o s h2x1L=s i n2(2 7 a)11动 力 学 与 控 制 学 报2 0 2 3年第2 1卷s i n h2 x1Lc o s2x2L+s i n h2 x2Lc o s2x1L+s i n h2c o s h2-c o s2=c T(2 7 b)系统参数、的确定则需要根据该超越方程组求解,这两个超越方程其实就是复波数方程(9)的实部与虚部,与M a i n在文献1 3 中求得的结果一致,方程(2 7 a)表达了、之间的关系,M a i n在文中称之为“相位方程”,方程(2 7 b)则表现了阻尼器对系统的影响.3 实例分析由式(2 6 b)可知,拉索的振幅函数与相位函数在阻尼器安装点处不可导,整根拉索被分为了两段,并且两段的边界相似,均为一端固定一端阻尼约束.因此为了方便分析,作如图4所示的理想模型考虑,选取阻尼器安装点的左半边索段L,并在索段右边提供一个不变的索力使结构受力平衡,索力方向始终与索右端截面垂直.图3 拉索-阻尼器半边模型F i g.4 H a l f-s i d e dm o d e l o f t h ec a b l e-d a m p e r s y s t e m取拉索左端固定点位置的相位(0)=0,根据边界条件求出振幅与相位函数表达式,并对振幅函数最大值归一,得到正则振幅函数y-(x-)与正则相位函数-(x-):y-(x-)=N c o s h(2 x-)-c o s(2 x-)-(x-)=(x-)=a r c t a nt a n(x-)t a n h(x-)+2n+12(2 8)式中x-=x/l,N 为最大值归一算子,n=0,1,2,用于保证-(x-)连续,无量纲阻尼与无量纲频率 根据方程(2 7)取值如下:小阻尼情况(c2T):=k2(k=2,4,6,)=a r c c o s hc2+Tc2-T(2 9 b)其中 只能取/2的整数倍,并且奇数倍只能在小阻尼(c2T)的情况下选取,拉索为反对称模态.由式(2 8)可知y-(x-)、-(x-)只与、有关,这里我们只关心阻尼器为拉索附加阻尼对拉索振幅与相位的影响,假设可以通过调节T、c使任意变化,图5、图6分别作出了=0和=0.1时y-(x-)、-(x-)、d-(x-)/dx-的函数图像;=0表示阻尼器对拉索不产生影响,=0.1表示阻尼器赋予拉索阻尼并使拉索运动衰减.(a)=0,=/2 (b)=0,=(c)=0,=3/2 (d)=0,=2 图5 无附加阻尼时拉索振幅、相位曲线F i g.5 Am p l i t u d ea n dp h a s ec u r v eo f t h ec a b l ew i t h o u ta d d i t i o n a l d a m p i n g由图5、图6可以看出,阻尼器为拉索结构附加的阻尼会影响拉索运动的振型与相位.在阻尼器附加阻尼为0的情况下,拉索微元两端的相位差d(x)/dx都为0,但是在驻点处相位突变了,因21第4期陈志等:拉索-阻尼器系统的改进实模态分析方法此运动方向突变为反向,相位差在驻点处为无穷大,这与经典的拉索模态分析结果一致.但在阻尼器附加阻尼不为0的情况下,拉索不再存在不动的驻点,相位也不是发生突变,而是有一个连续的变化过程,拉索总体表现为行波与驻波的叠加.(a)=0.1,=/2 (b)=0.1,=(c)=0.1,=3/2 (d)=0.1,=2 图6 附加阻尼时拉索振幅、相位曲线F i g.6 Am p l i t u d ea n dp h a s ec u r v eo f t h ec a b l ew i t ha d d i t i o n a l d a m p i n g4 总结本文对拉索-阻尼器系统进行了研究,首先采用复模态法也即复数分离变量法求解了该系统;随后依据复振型的物理意义对实数分离变量法进行了改进,得出了阻尼器通过相位差传递阻尼力的结论,并得到了一组有关拉索-阻尼器系统实际振幅和相位函数的非线性方程组;最后求解了该方程组,补充了复模态法对应的实模态求解方法,并简要分析了阻尼器对拉索振型与相位的影响.通过以上分析,本文得出了如下结论:(1)复振型的模为拉索实际振幅,辐角为拉索各点的运动相位.(2)在拉索阻尼器系统中,拉索通过相位差传递阻尼器产生的阻尼力,并且相位差还会对拉索的振型造成影响.(3)对于无阻尼拉索,其2阶以上的振型存在驻点,相位在驻点处突变了,相位差在该点为无穷大.(4)而安装了集中阻尼器的拉索,阻尼器为拉索结构附加了非对称阻尼,导致拉索不存在静止的驻点,相位在原驻点处也不发生突变,而是有一个连续的变化过程,表现出了行波与驻波叠加的特性.本文从拉索相位的角度分析了阻尼器对拉索产生的影响,该方法能够解释阻尼力的传递并描述阻尼器对振型的影响,可以为拉索的阻尼器设计提供参考.参考文献1N IYQ,WAN GXY,CHE NZQ,e t a l.F i e l do b-s e r v a t i o n so f r a i n-w i n d-i n d u c e dc a b l ev i b r a t i o n i nc a-b l e-s t a y e d D o n g t i n g L a k e B r i d g eJ.J o u r n a lo fW i n d E n g i n e e r i n g a n d I n d u s t r i a l A e r o d y n a m i c s.2 0 0 7,9 5:3 0 3-3 2 8.2J I NGH Q,HUAN GFY,HEXH,e t a l.W a k e-i n-d u c e dv i b r a t i o n so f t a n d e mf l e x i b l e c a b l em o d e l s i naw i n dt u n n e lJ.O c e a nE n g i n e e r i n g,2 0 2 1,2 3 3.3L I UZW,L ISQ,WAN GLH,e ta l.E x p e r i m e n-t a l i n v e s t i g a t i o no nh i g h-m o d ev o r t e x-i n d u c e dv i b r a-t i o no f a f l e x i b l e s t a yc a b l e i ns m o o t h f l o wJ.J o u r-n a l o fB r i d g eE n g i n e e r i n g.2 0 2 2,2 7(8):0 4 0 2 2 0 6 8.4CHE N Z Q,WAN G X Y,KOJ M,e ta l.MRd a m p i n gs y s t e mf o rm i t i g a t i n gw i n d-r a i n i n d u c e dv i-b r a t i o no nD o n g t i n gL a k ec a b l e-s t a y e db r i d g eJ.W i n da n dS t r u c t u r e s,2 0 0 4,7:2 9 3-3 0 4.5K L E I S S LK.C a b l ea e r o d y n a m i cc o n t r o l:w i n dt u n-n e ls t u d i e sD.C o p e n h a g e n:T e c h n i c a lU n i v e r s i t yo fD e n m a r k,2 0 1 3.6AHMA DJ,CHE N GSH.E f f e c t o f c r o s s-l i n ks t i f f-n e s so nt h ei n-p l a n ef r e ev i b r a t i o nb e h a v i o u ro fat w o-c a b l en e t w o r kJ.E n g i n e e r i n g S t r u c t u r e s,2 0 1 3,5 2:5 7 0-5 8 0.7KOVA C SI.Z u rF r a g ed e rS e i l s c h w i n g u n g e nu n dd e rs e i l d m p f u n gJ.B a u t e c h n i k,1 9 8 2,5 9(1 0):3 2 5-3 3 2.8P A CHE C O B M,F U J I NO Y,S U L E KH A.E s t i-m a t i o nc u r v ef o rm o d a ld a m p i n gi ns t a yc a b l ew i t hv i s c o u sd a m p e rJ.J o u r n a lo fS t r u c t u r a lE n g i n e e r-i n g,1 9 9 3,1 1 9(6):1 9 6 1-1 9 7 9.9Z HOU H,S UNL,X I NGF.D a m p i n go f f u l l-s c a l es t a yc a b l ew i t hv i s c o u sd a m p e r:e x p e r i m e n ta n da-n a l y s i sJ.A d v a n c e si nS t r u c t u r a lE n g i n e e r i n g,2 0 1 4;1 7(2):2 6 5-2 7 4.1 0F O S SKA.C o o r d i n a t e sw h i c hu n c o u p l e t h e e q u a t i o n31动 力 学 与 控 制 学 报2 0 2 3年第2 1卷o fd a m p e dl i n e a rd y n a m i cs y s t e mJ.J o u r n a lo fA p p l i e dM e c h a n i c s-A S ME,1 9 5 8,2 5,(1):3 6 1-3 6 4.1 1K R E NKS.V i b r a t i o no f a t a u t c a b l ew i t ha ne x t e r n a ld a m p e rJ.J o u r n a l o fA p p l i e sM e c h a n i c s,2 0 0 0,6 7(4):7 7 2-7 7 61 2K R E N KS,N I E L S E NS.V i b r a t i o no fas h a l l o wc a b l ew i t ha v i s c o u s d a m p e rC.P r o c e e d i n go f t h eR o y a l S o-c i e t yo fL o n d o nS e r i e sA:M a t h e m a t i c a lP h y s i c a la n dE n g i n e e r i n gS c i e n c e s,2 0 0 2,4 5 8:3 3 9-3 5 7.1 3MA I NJA,J O N E SNP.F r e ev i b r a t i o n so f t a u tc a b l ew i t ha t t a c h e dd a m p e r:l i n e a rv i s c o u sd a m p e rJ.J o u r n a lo fE n g i n e e r i n g M e c h a n i c s,2 0 0 2,1 2 8(1 0):1 0 6 2-1 0 7 1.1 4MA I NJA,J ON E SNP.F r e ev i b r a t i o no f t a u t c a b l ew i t ha t t a c h e sd a m p e r:n o n l i n e a rd a m p e rJ.J o u r n a l o fE n g i n e e r i n gM e c h a n i c s,2 0 0 2,1 2 8(1 0):1 0 7 2-1 0 8 1.1 5MA I N J A,J ON E S N P.E v a l u a t i o no fv i s c o u sd a m p e r sf o rs t a yc a b l ev i b r a t i o n m i t i g a t i o nJ.