端部激励相位差对悬索亚谐波共振响应的影响_李聪.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 6 7-0 0 9D O I:1 0.6 0 5 2/1 6 7 2-6 5 5 3-2 0 2 3-0 4 5 2 0 2 2-1 1-0 8收到第1稿,2 0 2 3-0 3-0 6收到修改稿.*国家自然科学基金资助项目(5 1 8 0 8 0 8 5),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 1 8 0 8 0 8 5).通信作者E-m a i l:s u n c s h n u.e d u.c n端部激励相位差对悬索亚谐波共振响应的影响*李聪 孙测世 赵碧航(重庆交通大学 土木工程学院,重庆 4 0 0 0 7 4)摘要 悬索受多个端部激励后可能产生复杂的非线性动力学现象.为研究端部激励相位差对悬索亚谐波响应的影响,建立两端受激励的悬索模型,得到无限维离散运动方程,采用多尺度法研究面内亚谐波共振,并通过数值积分加以验证.文中给出了多激励系统等效激励幅值和响应相位的一般表达式,定性定量分析了面内两两激励的不同组合对响应幅值和相位的影响.结果表明:端部激励相位差对悬索响应幅值的影响以2 为周期变化,关于对称.端部激励间的相位差可以使响应相位在一个名义单激励系统的响应相位基础上发生移动,当激励相位差逐渐由0转变为的同时,相频曲线也非同步地位移了.激励相位差不改变系统的软硬性质,但会影响响应幅值和共振区宽度,且其共振区宽度关于=0对称.关键词 悬索,端部激励,相位差,亚谐波共振,非线性振动中图分类号:O 3 2 2;TH 1 1 3文献标志码:AI n f l u e n c eo fE n dE x c i t a t i o nP h a s eD i f f e r e n c eo nS u b h a r m o n i cR e s o n a n t o fS u s p e n d e dC a b l e*L iC o n g S u nC e s h i Z h a oB i h a n g(C o l l e g eo fC i v i lE n g i n e e r i n g,C h o n g q i n gJ i a o t o n gU n i v e r s i t y,C h o n g q i n g 4 0 0 0 7 4,C h i n a)A b s t r a c t C o m p l e xn o n l i n e a rd y n a m i cp h e n o m e n aw i l lo c c u rw h e nt h es u s p e n d e dc a b l ei ss u b j e c t e dt om u l t i p l ep a r a m e t r i c e x c i t a t i o n s.I no r d e r t os t u d y t h e e f f e c t so f p a r a m e t r i c e x c i t a t i o nu n d e rd i f f e r e n t p h a-s e so nt h eh a r m o n i c r e s p o n s eo f s u s p e n d e dc a b l e,t h em o d e lo fs u s p e n d e dc a b l ew i t hp a r a m e t r i ce x c i t a-t i o na tb o t he n d s i se s t a b l i s h e d,a n dt h e i n f i n i t e-d i m e n s i o n a ld i s c r e t ee q u a t i o n so fm o t i o na r eo b t a i n e d.T h es u b h a r m o n i c r e s o n a n c e i np l a n e i s s t u d i e db yu s i n g t h em u l t i p l e s c a l e sm e t h o dv e r i f i e db yn u m e r i c a li n t e g r a t i o n.T h eg e n e r a l e x p r e s s i o n so f e q u i v a l e n t e x c i t a t i o na m p l i t u d ea n dr e s p o n s ep h a s eo fm u l t i-e x c i-t a t i o ns y s t e mi sg i v e n,a n dt h ee f f e c t so fd i f f e r e n tc o m b i n a t i o n so f i n-p l a n et w oe x c i t a t i o n so nt h er e-s p o n s ea m p l i t u d e a n dp h a s e a r eq u a l i t a t i v e l ya n a l y z e d.T h e r e s u l t s s h o wt h a t t h e e f f e c t o fm u l t i-p o i n t e x-c i t a t i o n sp h a s ed i f f e r e n c eo n t h e s u s p e n d e dc a b l e r e s p o n s e a m p l i t u d ev a r i e sw i t h2 a s ap e r i o d i c v a r i a t i o na n d i s s y mm e t r i c a b o u t.T h e s e e x c i t a t i o n sm a k e t h e r e s p o n s ep h a s eo f t h e s y s t e ms h i f t b a s e do n t h e r e-s p o n s ep h a s eo f t h e s u s p e n d e dc a b l e.Wh e n t h ep h a s ed i f f e r e n c e e x c i t a t i o nb e t w e e n t h e mg r a d u a l l yc h a n-g e s f r o m0t o,t h ep h a s e-f r e q u e n c yc u r v e sa l s os h i f tb y.T h ee x c i t a t i o np h a s ed i f f e r e n c ed o e sn o tc h a n g e t h es o f t a n dh a r dp r o p e r t i e so f t h es y s t e m,b u ta f f e c t st h er e s p o n s ea m p l i t u d ea n dt h ew i d t ho ft h er e s o n a n c er e g i o n,a n dt h ew i d t ho f t h er e s o n a n c er e g i o n i ss y mm e t r i c a l a b o u t=0.K e yw o r d s s u s p e n d e dc a b l e,p a r a m e t r i ce x c i t a t i o n,p h a s ed i f f e r e n c e,s u b h a r m o n i cr e s o n a n c e,n o n l i n e a rv i b r a t i o n动 力 学 与 控 制 学 报2 0 2 3年第2 1卷引言悬索受到的激励形式主要是风、雨激励等外激励和桥塔竖向的位移、车辆引起主梁的振动导致索的端部激励.实际上在悬索桥的使用中,悬索必然受到多点多激励,它的非线性响应与激励的幅值、频率、相位和激励方向有关.从索的非线性运动方程出发,悬索两端的端部激励虽不改变其二次和三次非线性项系数,但是会出现在参数激励和强迫振动项中,从而影响索的振动响应性质.悬索在多点多激励作用下可能产生大幅振动和非线性动力特性1-4,较为常见的有超谐波共振5、主参数共振6和主共振7.由于多个端部激励可能具有不同的初始状况,如激励幅值、相位等,所以有必要确定激励相位差对响应的影响.已有相关研究表明:激励相位差不仅能改变悬索的振动模态和幅频响应特性8,而且能加剧与激励频率相近的面内、外模态振动.在任意端部激励组合下,激励相位差不仅可使斜拉索非线性振动出现定量变化,还可改变内共振的表现形式9.E l-A t t a r等1 0分析了考虑初始相位差的悬索在两种垂直/横向支承激励下的非线性响应,发现相位差对幅频曲线有显著影响.Z a g h a r i1 1,1 2等发现最大/最小响应振幅取决于不同激励源之间的相位差.本文研究不同激励组合下的激励相位差对悬索稳态响应幅值和相位的影响.激励与响应的相位关系有助于进一步研究振动控制或能量收集器策略1 3,1 4.在悬索非线性动力学模型中引入上下端两向激励的幅值和相位,通过数值分析,绘制各工况下的幅频曲线和相频曲线,分析在不同端部激励两两组合下的悬索非线性振动特性.1 力学模型1.1 基本假定本文研究了小垂跨比(f=b/l1/8)的悬索结构并建立如图1所示的两端面内激励下悬索的非线性振动模型.悬索的两端锚固都受到水平、竖向端部激励,分别用U(t)、V(t)表示,其中用下标a、b表示左右锚点,激励幅值和相位分别为zs、s(s=1,2,3,4),数字1到4是各个方向激励的激发,假设端部激励的激励频率大小一致,悬索长取l和悬索垂度取b.考虑到悬索振动方向与端部激励方向的一致性,此处需建立局部振动坐标系o-x y,其中坐标原点o为左端锚固点,x为悬索的索向坐标,y为索面内垂直索向下的坐标,另外,各方向对应位移分别用u(x,t)、v(x,t)表示.并且做如下假设:悬索的抗弯刚度、扭转刚度和剪切刚度足够小以至于可以忽略不计;悬索只承受拉力;悬索在振动过程中的轴向应变足够小;只考虑几何非线性,不考虑其他非线性.图1 两端面内激励悬索模型F i g.1 S u s p e n d e dc a b l em o d e lw i t h i n-p l a n e e x c i t a t i o na tb o t he n d s1.2 动力平衡方程鉴于 悬 索 的 几 何 非 线 性 和 两 端 铰 接,根 据H a m i l t o n能量守恒变分原理,可得到悬索拟静态面内非线性动力学平衡方程:m u+cuu-E A u+y v+v 2+w 22=0(1)m v+cvv-H v-E Ay+v()u+y v+v 2+w 22=0(2)式中:上标一点表示对时间t求导,右上一撇表示对x坐标求导;w为面外方向对应位移;m为悬索单位长度质量;cu、cv分别为悬索沿u、v方向上的阻尼系数;H为悬索的水平张力;E为弹性模量;A为悬索横截面面积;y x()=4f x1-x()为悬索静态构型,其中f为垂跨比.假设悬索的轴向振动是准静态的,忽略高阶无穷小,整理简化悬索的无量纲控制方程得:m v+cvv-v 2-2(y+v)U+10y v+v 2+w 22dx=0(3)式中:U=Ub(t)-Ua(t).在无量纲过程中,为了简化方程的表达式,省略了“*”号,使用了以下变量.x*=xl,y*=yl,v*=vl,w*=wl,=E AH,86第4期李聪等:端部激励相位差对悬索亚谐波共振响应的影响t*=lHmt,U*=Ul,c*v=l1mHcv1.3 离散模型为了将系统从偏微分方程离散到常微分方程,将面内位移函数v看作是准静态运动和模态运动两部分的组合:v(x,t)=Nn=1n(t)n(x)+Vb(t)x+Va(t)(1-x)(4)其中,n(t)是第n个广义时间坐标;n(x)是第n个模态形状;N是平面内的模态数;Vb(t)和Va(t)分别表示左、右端的竖向激励(如图1);等式(4)右边第一项为模态位移,后两项为端部激励产生的准静态位移.悬索的面内模态形状可以通过求解线性化的运动方程方程得到,即面内模态函数为:n(x)=knc o s(n x)-1+t a nn2s i n(n x),n=1,3,n(x)=2 s i n(n x),n=2,4,(5)式中:kn=2nc o s2n2(2+c o sn)n-3 s i nnt a nn2=n2-11 6n23将式(4)代入式(3)后通过G a l e r k i n方法,采用模态离散得到端部激励下悬索面内的振动微分方程:k+2v kk+2kk=22Nn=1Nr=1QkPn r+2QrPn k()nr-22Nn=1Pn kn-22Nn=1Nr=1Nj=1Pn kPr jnrj+22Qk-R v-2v kRv(6)式中:2v k=10cv2kdx,k=1,2,;Qn=10ny dx,Pr j=10 r jdx=2Ub(t)-Ua(t)+Vb(t)-Va(t)2Rv=10Vb(t)x+Va(t)(1-x)kdx2 多尺度摄动基于多尺度法引入无量纲参数,假设方程(6)的广义时间坐标近似解展开前三阶为:k(t;)k1(T0,T1,T2)+2k2(T0,T1,T2)+3k3(T0,T1,T2)(7)式中:T0,T1,T2是3个不同尺度的时间变量,其变化的速度随时间t一次减慢一个数量级.由链式法则对时间t求导可转变为对Ti(i=1,2,3,)求偏导数,可得:ddt=dT0dtT0+dT1dtT1+dT2dtT2+=D0+D1+D2+d2dt2=D20+2 D0D1+2(2D0D2+D21)+(8)其中:Dn/Tn(n=1,2,3,)设zi2zi(i=1,2,.);v k2v k(k=1,2,);将式(7)-式(8)代入常微分方程组(6),然后根据的阶次大小顺序排列可得:1:D20k1+2kk1=0(9)2:D20k2+2kk2=22Nn=1Nr=1(QkPn r+2QrPn k)nr-2D0D1k1-R v+2QkUb(t)-Ua(t)(1 0)3:D20k3+2kk3=-2D0D1k2-2D0D2k1-2vD0k1-D21k1+22Nn=1Nr=1(QkPn r+2QrPn k)(n1r2+n2r1)-22Nn=1Nr=1Nj=1Pn kPr jn1r1j1-2Nn=1Pn kn1Ub(t)-Ua(t)(1 1)2.1 模态分析激励频率中引入调谐参数,使得:=2m+2(1 2)根据式(9)可设:m1=Am(T1,T2)eimT0+c ck1=0,km(1 3)式中:c c代表各项的共轭复数.96动 力 学 与 控 制 学 报2 0 2 3年第2 1卷把式(1 3)代入式(1 0)得:D20k2+2kk2=-2D0D1(AmeimT0+c c)+22(QkPmm+2QmPm k)(AmeimT0+c c)2+22Qk4s=1Hsei(T0+s)+c c+224s=1Vsei(T0+s)+c c(1 4)式中:H1=z1,H2=-z2,H3=H4=0,V1=V2=0;V3=z310 x ndx,V4=z410(1-x)ndx将式(1 2)代入式(1 4)得:D20k2+2kk2=-2D0D1(AmeimT0+c c)+22(QkPmm+2QmPm k)(AmeimT0+c c)2+22Qk4s=1Hsei(2mt+T2+s)+c c+(2m+2)224s=1Vsei(2mt+T2+s)+c c(1 5)对式(1 5)化简长期项:D1Am=0(1 6)随后,对式(1 4)求解得:k2=E1k(A2me2 imT0+c c)+E2kAmA-m+4s=1fsei(To+s)2k-2+c c,k(1 7)式中:E1k=22QkPmm+2QmPm k2k-42m,kmE2k=2QkPmm+2QmPm k2kfs=22QkHs+22Vs,s=1,2,3,4这里值得注意的是,当k=2m,导致E1k分母为0,非线性系统将出现无解.把式(1 2)、式(1 3),式(1 6)、式(1 7)代 入 式(1 1),并令k=m,化简长期项得:2 imD2Am+vAm()+e mA2mA-m-4s=1F-sA-mei T2+s()=0,(1 8)式中:e m=322P2mm-224Nr=132r-82m(2r-42m)2r(QrPmm+2QmPm r)2F-s=2Nr=1(QrPm m+2QmPm r)2fs2r-2-12Pm mHs以极坐标的形式表示复变函数Am:Am=12ameim(1 9)将式(1 9)代入式(1 8),并设=T2-2m,a m=0和=0,化简并分离虚实部得:2mvam=am4s=1F-ss i n(+s)mam+e ma3m4=am4s=1F-sc o s(+s)(2 0)将式(2 0)中的虚部和实部求平方和得频率响应方程:(2mv)2+14e ma2m-m2=F2(2 1)式中:F2=4s=1F-ss i ns()2+4s=1F-sc o ss()2(2 2)2.2 响应分析(1)响应幅值可以从亚谐波共振面内频率响应方程式(2 1)得出频响方程解析解:a2m=4me m4e mF2-(2vm)2(2 3)若响应幅值要为稳态解,从式(2 3)可以得出F2-(2vm)20.当系统偏硬时,即e m0且满足F2vm,系统满足上诉的条件时才能产生持续不断的运动.当F2-(2vm)2/m时,系统会有两个稳态解;当|F2-(2vm)2/m时,系统会有一个稳态解.当系统偏软时,即e m2vm,当-F2-(2vm)2/m时,系统会有两个稳态解;当|F2-(2vm)2/m时,系统会有一个稳态解.从上可以看出,偏软和偏硬非线性系统的稳态解范围是关于=0成对称状态,并且响应幅值am与参数F存在一定的关系,本文希望以后进一步分析.(2)响应相位根据式(2 2)可设:Fs i n-=4s=1F-ss i n(+s)Fc o s-=4s=1F-sc o s(+s)(2 4)07第4期李聪等:端部激励相位差对悬索亚谐波共振响应的影响根据式(2 0)和式(2 2),虚部与实部相除得:t a n-=2mvm+14e ma2m(2 5)式(2 5)与初始相位s=0的单激励非线性振动的相位频率方程相同,这里称之为名义单激励相位差系统,其中F和-是名义单激励系统的幅值和相位.若Fc o s-0,则结合方程组(2 4)中两式相除得:=-a r c t a n4s=1F-ss i ns4s=1F-sc o ss=-k(2 6)式中的是悬索在多个不同初始相位的端部激励下的响应相位,本文将关系式(2 6)称为多激励相位差系统,它反映了所有端部激励对面内亚谐波共振下响应相位的影响.观察发现,响应相位因各个激励的初始相位s和参数F-s的存在而发生不同程度的偏移(相移值),其中参数F-s综合考虑了索的物理性质,激励幅值、频率以及模态形状.3 数值分析本节分析了多个端部激励相位差对悬索亚谐波响应幅值和相位的影响,但没有研究分岔问题,故不区分稳态解的稳定性,希望以后能够进一步研究.对式(2 1)用数值方法求解,为了验证多尺度法解的正确性,对式(6)采用龙格-库塔法验证.采用悬索的物理参数如下1 6:横截面积A=7.0 6 91 0-2m2,单位长度质量m=5 5 1.3 8 2 g/m,弹性模量:E=2.01 01 1P a,悬索长l=2 0 0 m,垂跨比f=0.0 0 5,无量纲阻尼比v k=0.0 0 2.对于面内索激励组合下,其相位差对索响应幅值的影响是以为对称,2 为周期变化9,则本文将讨论3种以两两激励组合相位差为0的第一阶面内对称模态(m=1)的频率响应曲线和相频曲线.3.1 工况1:两个水平端部激励固定左端激励相位a u(2)=/4的初始相位和两个支承的无量纲激发振幅分别为au=0.0 0 0 1和bu=0.0 0 0 3,改变右端激励相位b u(1),区间为/4,+/4.仅展示左端激励相位2=/4的初始相位,右端激励相位1=/4以及两个支承的无量纲激发振幅分别为au=0.0 0 0 1和bu=0.0 0 0 3的龙格库塔法值,如图2所示,龙格库塔法的值与多尺度计算的结果吻合良好,当右端激励相位1从/4增长到5/4,内侧的幅频曲线向外侧平移.在工况1激励作用下,激励相位差不会改变系统的软硬特性,但它会对响应幅值和共振区宽度产生影响.当激励相位1=/4时,其龙格库塔法最大响应幅值为0.1 5 4,共振区的范围为-0.1 0 7,0.1 0 7(共振区呈=0对称);当1从/4增长到5/4时,即相位差b u从0变化到,系统的响应幅值和共振区逐渐增大,成正比关系.若悬索保持它的物理参数和激励幅值不变,所有的频率响应曲线都会在曲线1=/4(下界)和1=5/4(上界)之间.图2 工况1:幅频响应曲线F i g.2 C o n d i t i o n1:a m p l i t u d e-f r e q u e n c yr e s p o n s ec u r v e s在工况1的幅频响应曲线基础上绘制相频曲线,如图3所示,名义单激励相频曲线和多激励相频曲线的响应相位几乎都不随激励频率的改变而变化,而为一个稳定值.从图3(a)可以看出,无论是较大的稳态解还是较小的稳态解,所对应曲线1=/4和1=5/4的响应相位-差值都很小且相同,其差值为=0.0 3 9.从图3(b)可以看出,响应相位分别被激发,正如式(2 6),激励相位差加剧了响应相位-的平移.系统的响应相位分别激励到2.4 3 3,1.3 3 6,0.5 1 3,-0.1 4 8和-0.7 4 7,响应相位的1=/4(下界)值2.4 3 3与1=5/4(上界)值-0.7 4 7之差为3.1 8,其值减去差值=0.0 3 9约为,即相位差的变化区间为,而响应相位也随着移动了.17动 力 学 与 控 制 学 报2 0 2 3年第2 1卷(a)名义单激励相频曲线(a)N o m i n a l s i n g l ee x c i t a t i o np h a s e-f r e q u e n c yc u r v e s(b)多激励相频曲线(b)M u l t i-e x c i t a t i o np h a s e-f r e q u e n c yc u r v e s图3 工况1:相频曲线F i g.3 C o n d i t i o n1:p h a s e-f r e q u e n c yc u r v e s3.2 工况2:两个竖向端部激励固定左端激励相位a v(4)=/4的初始相位和两个支承的无量纲激发振幅分别为av=0.0 0 0 1和bv=0.0 0 0 3,改变右端激励相位b v(3),区间为/4,+/4.仅展示左端激励相位4=/4的初始相位,右端激励相位3=/4以及两个支承的无量纲激发振幅分别为av=0.0 0 0 1和bv=0.0 0 0 3的龙格库塔法值,如图4所示,两者吻合良好,在两竖向激励作用下,龙格库塔法的最大响应幅值激励到0.0 0 2 9后就迅速跳跃到0值,响应幅值远小于两水平激励下的响应幅值0.1 5 4,说明悬索中的水平激励作用远大于竖向激励的作用.但与其他工况相比,工况2的数值验证效果较差,其原因可能是竖向激励下的响应集中于=0附近,在数值计算中取离散的初始值进行扫频,而工况2下的非线性动力系统对初始值很敏感.同时,响应幅值在两竖向激励作用下能够被激励起的范围很小,共振区很窄,大约在-0.0 0 5,0.0 0 5(共振区呈=0对称),相对于两水平激励作用下共振区范围的1/2 0.图4 工况2:幅频响应曲线F i g.4 C o n d i t i o n2:a m p l i t u d e-f r e q u e n c yr e s p o n s ec u r v e s(a)名义单激励相频曲线(a)N o m i n a l s i n g l ee x c i t a t i o np h a s e-f r e q u e n c yc u r v e s(b)多激励相频曲线(b)M u l t i-e x c i t a t i o np h a s e-f r e q u e n c yc u r v e s图5 工况2:相频曲线F i g.5 C o n d i t i o n2:p h a s e-f r e q u e n c yc u r v e s从图4看出,当右端激励相位3从/4增长到5/4,外侧的幅频曲线向内侧平移,当激励相位差增大时,响应幅值和共振区共同减小,成反比关系.但与工况1相比,两者的幅频曲线平移方向刚27第4期李聪等:端部激励相位差对悬索亚谐波共振响应的影响好相反,当两竖向激励达到最小值时,而两水平激励达到最大值,从中可以看出两者互为反向的振动.工况2与工况1相比,其相位差没有改变频率响应曲线的弯曲方向,但影响了使响应幅值更小,共振宽度更窄.当3从/4增长到5/4时,即相位差从0变化到,系统的响应幅值和共振区逐渐减小;与工况1相反的是,当相位差为0或2 时,响应幅值取得最大值,而当相位差为时,响应幅值取得最小值.绘制相频曲线,如图5所示,响应相位-影响较大,相对于工况1的相频曲线敏感很多,随激励频率的变化而变化,并且存在上下两条曲线交会于一点,然后响应相位-达到一个稳定值/2,呈9 0,其中=0.从图5(b)可以看出,响应相位分别激励到0.7 8 5,0.1 8 8,-0.4 6 3,-1.2 7 1和3.9 2 7,3=5/4响应相位值3.9 2 7与3=/4响应相位值0.7 8 5之差为3.1 4 2,其值约为,即名义单激励响应相位经过两竖向激励相位差平移了.与两水平激励相位差类似的是,其响应相位存在一个向上移动的范围,但与两水平激励相位差不同的是,3=5/4的响应相位并没有随着3=的响应相位继续往负数响应相位移动,而是发生了一个相位的跳跃,直接从负数响应相位转换成一个正数响应相位,并且从相位差的步长大小可以看出这个转换是迅速的,为非线性变化.与工况1不同,工况2的相位差下的一阶响应相位会随激励频率改变而改变,相位值对激励频率的改变较敏感.3.3 工况3:一个水平和一个竖向端部激励固定左端激励相位a u(2)=/4的初始相位和两个支承的无量纲激发振幅分别为au=0.0 0 0 1和bv=0.0 0 0 3,改变右端激励相位b v(3)/4,+/4.仅展示左端激励相位2=/4的初始相位,右端激励相位3=/4以及两个支承的无量纲激发振幅分别为au=0.0 0 0 1和bv=0.0 0 0 3的龙格库塔法值,如图6所示,两者吻合良好,随激励频率的增长,响应幅值和共振区范围也增大,这个范围为-0.0 3,0.0 3(共振区呈=0对称).当右端激励相位3从/4增长到5/4,外侧的幅频线向内侧平移,与工况1的情况相反.竖向激励的相位与水平激励的相位相比,前者对幅频曲线的平移规律更有一定的影响.并且工况3的共振区范围大小位于工况1和工况2之间,对于工况1的共振区而言,其共振区范围约为其的1/3;而对于工况2的共振区而言,其共振区范围为其的6倍.可见,与竖向激励之间的相位差相比,水平激励之间的相位差会加大共振区范围.图6 工况3:幅频响应曲线F i g.6 C o n d i t i o n3:a m p l i t u d e-f r e q u e n c yr e s p o n s ec u r v e s(a)名义单激励相频曲线(a)N o m i n a l s i n g l ee x c i t a t i o np h a s e-f r e q u e n c yc u r v e s(b)多激励相频曲线(b)M u l t i-e x c i t a t i o np h a s e-f r e q u e n c yc u r v e s图7 工况3:相频曲线F i g.7 C o n d i t i o n3:p h a s e-f r e q u e n c yc u r v e s37动 力 学 与 控 制 学 报2 0 2 3年第2 1卷当相位差b v=b v-a u=0和=0时,工况3最大响应幅值0.0 0 5 3位于同种情况下的工况1和工况2响应幅值之间.此时工况1的响应幅值0.0 0 6 7(最小幅值),约为工况3的响应幅值的1.2 6倍;而工况2响应幅值0.0 0 2 7(最大幅值),约为其响应幅值的0.5.当相位差值b v=和=0时,工况3的响应幅值为0.0 0 4 6(最小幅值)且位于两水平激励和两竖向激励的响应幅值之间.工况1的响应幅值0.0 0 9 4(最大幅值),约为其响应幅值的2倍;工况2的响应幅值0.0 0 2(最小幅值),约为其响应幅值的0.5倍.由此可见,当激励间相位差b v=的响应幅值与无相位差b v=0相比,在激励组合中的水平激励相位差与竖向激励相位差诱导的响应幅值相比,前者占绝对主导位置.从上分析可知,竖向激励相位差会抑制水平激励相位差所引起系统的响应幅值,但是水平激励对响应幅值的影响远大于竖向激励.同理,绘制出相频曲线,如图7所示,可以观察得出此组合下的激励相位差变化对响应相位-影响不大,由局部放大图观察得知,与两水平激励作用下类似,所有曲线最终交会于一点,但-值变化不明显,即-值以水平激励作用为主导.如图7(b)所示,响应相位分别激励到-0.6 4 5,-0.6 1 0,-0.7 4 1,-0.7 1 4和-0.7 0 2,3=5/4响应相位值-0.6 1 0与3=/4响应相位值-0.6 4 5之差为0.0 3 5,其值减去(0.0 3 5)为0,并且响应相位都为负值.这与以上2种工况激励相位差不同,它们的响应相位移动约为,并且它的响应相位与以上两种组合工况的初始增长方向相反.综上,该工况下的相位差与响应幅值图呈反比,与两水平激励相位差的相位差与响应幅值图呈正比相反,故竖向激励会抑制水平激励对响应幅值的影响.因竖向激励对名义单激励相位值更加敏感,若激励组合中存在一竖向激励,其相频曲线会汇聚于一点,但这个点的激励频率大小并不一致,且相位的大小还是以水平激励为主.4 结论(1)承受多端部激励的悬索,激励间的相位差可以使响应相位在名义单激励系统的响应相位基础上发生移动,相移值与不同激励组合的相位差、激发模态、激励幅值和初始相位有关.(2)当激励相位差从0逐渐变化到时,相频曲线也发生不同程度的移动,最终相位变化大小为.激励相位差不改变系统的软硬性质,但会影响响应幅值和共振区宽度,且共振区宽度关于=0对称.(3)水平激励对响应幅值大小和共振区宽度的影响远大于竖向激励.系统主要靠参数振动激发,但竖向激励会抑制水平激励作用效果.响应相位对水平激励初始相位的敏感程度远小于对竖向激励初始相位的敏感程度.参考文献1丛云跃,康厚军,郭铁丁,等.C F R P索斜拉桥面内自由振动的多索梁模型及模态分析 J.动力学与控制学报,2 0 1 7,1 5(6):4 9 4-5 0 4.C ONGYY,KAN G HJ,GUOTD,e t a l.Am u l-t i p l ec a b l e-b e a m m o d e la n d m o d a la n a l y s i so ni n-p l a n e f r e ev i b r a t i o no f c a b l e-s t y e db r i d g ew i t hC F R Pc a b l e sJ.J o u r n a l o fD y n a m i c sa n dC o n t r o l,2 0 1 7,1 5(6):4 9 4-5 0 4.(i nC h i n e s e)2康厚军,赵跃宇,蒋丽忠.参数振动和强迫振动激励下超长拉索的面内非线性振动 J.中南大学学报(自然科学版),2 0 1 1,4 2(8):2 4 3 9-2 4 4 5.KAN G HJ,Z HAOYY,J I AN GLZ.I n-p l a n en o n-l i n e a rv i b r a t i o no f s u p e r l o n gs t a yc a b l e su n d e rp a r a-m e t r i ca n df o r c e de x c i t a t i o n sJ.J o u r n a lo fC e n t r a lS o u t hU n i v e r s i t y(S c i e n c ea n dT e c h n o l o g y),2 0 1 1,4 2(8):2 4 3 9-2 4 4 5.(i nC h i n e s e)3陈杰夫,康厚军,苏潇阳,等.谐波激励下斜拉桥面内非线性振动试验研究 J.固体力学学报,2 0 1 9,4 0(3):2 4 8-2 5 9.CHE NJF,KANG HJ,S HU XY,e ta l.N o n l i n-e a rd y n a m i c so fap a r a m e t r i c a l l ya n de x t e r n a l l ye x-c i t e dt h i np l a t eJ.C h i n e s eJ o u r n a lo fS o l i d M e-c h a n i c s,2 0 1 9,4 0(3):2 4 8-2 5 9.(i nC h i n e s e)4吴娟,钱有华.一类弦-梁耦合非线性振动系统的动力学数值模拟研究 J.动力学与控制学报,2 0 1 8,1 6(5):4 0 3-4 1 0.WUJ,Q I AN Y H.N u m e r i c a ls i m u l a t i o nr e s e a r c ho nd y n a m i c so fas t r i n g-b e a mc o u p l e dn o n l i n e a rv i-b r a t i o ns y s t e mJ.J o u r n a lo fD y n a m i c sa n dC o n-t r o l,2 0 1 8,1 6(5):4 0 3-4 1 0.(i nC h i n e s e)5赵跃宇,李永鼎,王连华,等.悬索的超谐波共振与13内共振分析 J.动力学与控制学报,2 0 0 7,547第4期李聪等:端部激励相位差对悬索亚谐波共振响应的影响(2):1 1 2-1 1 7.Z HAOYY,L IYD,WANGL H,e ta l.O no n e-t o-t h r e e i n t e r n a l r e s o n a n c e s a n ds u p e r h a r m o n i c r e s o-n a n c e so f s u s p e n d e dc a b l e sJ.J o u r n a l o fD y n a