温度变化对悬索全局动力学特性影响_吴先强.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 3 2-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 1 2 0 2 2-1 2-2 5收到第1稿,2 0 2 3-0 4-0 3收到修改稿.*国家自然科学基金面上项目(1 2 2 7 2 1 3 9),福建省自然科学基金面上项目(2 0 2 2 J 0 1 2 9 0),福厦泉国家自主创新示范区协同创新平台项目(3 5 0 2 Z C Q X T 2 0 2 2 0 0 2),N a t i o n a lN a t u r a l S c i e n c eF o u n d a t i o no fC h i n a(1 2 2 7 2 1 3 9),N a t u r a l S c i e n c eF o u n d a t i o no fF u j i a nP r o v i n c e(2 0 2 2 J 0 1 2 9 0),C o l-l a b o r a t i v e I n n o v a t i o nP l a t f o r mP r o j e c t o fF u z h o u-X i a m e n-Q u a n z h o uN a t i o n a l S e l f-I n n o v a t i o nZ o n e(3 5 0 2 Z C Q X T 2 0 2 2 0 0 2).通信作者E-m a i l:y b z h a o h q u.e d u.c n温度变化对悬索全局动力学特性影响*吴先强 赵珧冰 郭智锐 陈林聪(华侨大学 土木工程学院,厦门 3 6 1 0 2 1)摘要 索是一类工程中常用的张力结构,其柔度大、阻尼轻,在各类外部荷载作用或端部位移激励下极易发生大幅振动,影响结构安全运营.已有研究表明悬索的振动特性对于温度变化极为敏感,因此本文同时考虑支座运动引起的参数共振以及模态间12内共振,基于全局分岔理论,系统探究温度变化对悬索全局动力学行为的影响.首先引入张力改变系数,建立考虑整体温度变化影响与受参数激励悬索的面内非线性运动微分方程.采用G a l e r k i n法进行离散,利用多尺度法得到该非线性系统直角坐标形式的平均方程,并基于坐标变换,将平均方程简化为规范形,采用能量相位法研究温度变化时悬索多脉冲混沌动力学行为.通过能量差函数的零点条件以及扰动系统下中心点的吸引域范围,分析激励幅值、阻尼系数和调谐参数的取值范围,并计算该四维系统的L y a p u n o v指数.研究结果表明:温度变化会影响系统S h i l n i k o v型多脉冲同宿轨道的产生;随着温度变化,多脉冲同宿轨道可能消失,导致系统的混沌运动转变为周期运动;受温度变化影响,动力系统可能展现出截然不同的动力学行为.关键词 悬索,温度变化,全局分岔,能量相位法,混沌中图分类号:O 3 2 2文献标志码:AI n f l u e n c e so fT e m p e r a t u r eC h a n g e so nG l o b a lD y n a m i c a lC h a r a c t e r i s t i c so fS u s p e n d e dC a b l e s*WuX i a n q i a n g Z h a oY a o b i n g G u oZ h i r u i C h e nL i n c o n g(C o l l e g eo fC i v i lE n g i n e e r i n g,H u a q i a oU n i v e r s i t y,X i a m e n 3 6 1 0 2 1,C h i n a)A b s t r a c t T h e s u s p e n d e dc a b l e i s a t y p ep f c o mm o n l yu s e d t e n s i o ns t r u c t u r e i ne n g i n e e r i n gs t r u c t u r e s.I th a sh i g h f l e x i b i l i t ya n d l i g h t d a m p i n g,a n d i t i sp r o n e t o l a r g ev i b r a t i o n su n d e r v a r i o u s e x c i t a t i o n s a n d/o rs u p p o r tm o t i o n s,w h i c he n d a n g e r st h es a f e t yo f t h ec a b l es t r u c t u r e s.T h ep r e v i o u ss t u d i e sh a v es h o w nt h a t t h ev i b r a t i o nc h a r a c t e r i s t i c s o f t h ed y n a m i c s y s t e ma r ev e r ys e n s i t i v e t o t e m p e r a t u r e c h a n g e s.T h e r e-f o r e,t h i sp a p e rc o n s i d e r sb o t ht h ep a r a m e t r i c r e s o n a n c e s c a u s e db y t h e s u p p o r tm o t i o n s a n d t h e t w o-t o-o n e i n t e r n a l r e s o n a n c e s,a n d t h e n i t s y s t e m a t i c a l l ye x p l o r e s t h e i n f l u e n c e so f t e m p e r a t u r e c h a n g e so n t h es u s p e n d e dc a b l esd y n a m i cb e h a v i o r sf r o mt h ep e r s p e c t i v eo fg l o b a lb i f u r c a t i o n s.F i r s t l y,t h et e n s i o nv a r i a t i o nc o e f f i c i e n t i s i n t r o d u c e d,a n dt h en o n l i n e a rd i f f e r e n t i a l e q u a t i o n so f t h e i n-p l a n em o t i o no f t h es u s p e n d e dc a b l es u b j e c t e dt op a r a m e t r i ce x c i t a t i o n i nt h e r m a l e n v i r o n m e n t s i se s t a b l i s h e d.T h eG a l e r k i nm e t h o d i su s e df o rd i s c r e t i z a t i o n,a n dt h em u l t i p l es c a l e sm e t h o d i sa d o p t e dt oo b t a i nt h ea v e r a g ee q u a-t i o no f t h en o n l i n e a rs y s t e mi nr e c t a n g u l a rc o o r d i n a t e s.B a s e do nt h ec o o r d i n a t et r a n s f o r m a t i o n,t h ea v-e r a g ee q u a t i o n i ss i m p l i f i e d i n t oan o r m a l f o r m.T h ee n e r g yp h a s em e t h o di su s e dt oe x p l o r et h em u l t i-第4期吴先强等:温度变化对悬索全局动力学特性影响p u l s ec h a o t i cd y n a m i c a lb e h a v i o ro f t h es u s p e n d e dc a b l ew h e nt h et e m p e r a t u r e i sc h a n g e d.T h r o u g ht h ez e r o-p o i n t c o n d i t i o no f t h e e n e r g y-d i f f e r e n c e f u n c t i o na n d t h e r a n g eo f a t t r a c t i o no f t h e c e n t e rp o i n t u n d e rt h ed i s t u r b a n c es y s t e m,t h ee x c i t a t i o na m p l i t u d e,d a m p i n gc o e f f i c i e n ta n dd e t u n i n gp a r a m e t e r so ft h es y s t e ma r e e x p l o r e d,a n d t h eL y a p u n o ve x p o n e n t so f t h e f o u r-d i m e n s i o n a l s y s t e ma r e c a l c u l a t e d.N u m e r-i c a l e x a m p l e ss h o wt h a t:t e m p e r a t u r ec h a n g e sa f f e c tt h eg e n e r a t i o no fs y s t e msS h i l n i k o v-t y p em u l t i-p u l s eo r b i t s.T h em u l t i-p u l s eo r b i t sm a yd i s a p p e a r c o n s i d e r i n g t e m p e r a t u r e e f f e c t s,a n d i t c a u s e s t h e s y s-t e msc h a o t i cm o t i o n s t r a n s f o r mi n t op e r i o d i cm o t i o n s.T h ed y n a m i c a l s y s t e mm a ye x h i b i t d i s t i n c t v i b r a-t i o nb e h a v i o r s i nt h e r m a l e n v i r o n m e n t s.K e yw o r d s s u s p e n d e dc a b l e,t e m p e r a t u r ec h a n g e,g l o b a lb i f u r c a t i o n,e n e r g y-p h a s em e t h o d,c h a o s引言由于索结构线形美观、柔度大、质量轻以及强度高等特点,被广泛运用于各类工程中.众所周知,各类索结构极易受到外界环境因素影响1,例如在太阳辐射等影响下,温度场会不断变化,从而引起结构热胀冷缩.国内外学者对索结构在温度场中的静、动力特性开展了大量研究2-6.此外,受支座位移运动引发的参数共振,将导致拉索展现出复杂的非线性动力学行为7,8.倘若进一步考虑温度变化影响,系统的振动特性会产生明显的定性和定量变化,进而影响索结构的安全性和稳定性.近年来随着全局分岔理论的不断完善,诸如混沌之类的复杂非线性动力现象被人们进一步了解和认知.张伟等9全面系统总结了高维非线性系统中多脉冲混沌理论,指出了全局分岔理论在工程应用中的现状和展望.尤其对于悬索这类典型的、同时包含平方和立方非线性的系统,全局分岔理论可以很全面地描述其非线性动力学行为.Z h a n g和T a n g1 0基于全局摄动方法研究了悬索在参数激励和外界环境激励耦合作用下的全局分岔现象和混沌行为.C h e n和X u1 1基于能量相位法研究了斜拉索的S h i l n i k o v同 宿 轨 道 以 及 混 沌 动 力 学 行 为.Z h a n g等1 2借助高维M e l n i k o v方法证明了张紧弦在12内共振情况下存在S h i l n i k o v同宿轨道.除悬索、斜拉索和弦等张力结构外,有学者通过广义M e l n i k o v方法和能量相位法研究了板、壳和天线等基本结构的混沌动力学和多脉冲轨道.A n和C h e n1 3研 究 了 气 动 热 弹 性 功 能 梯 度 材 料(F GM)截锥壳在复杂载荷作用下12内共振和主参数共振情况下全局分岔和多脉冲轨道.高美娟等1 4依据三阶剪切变形层合板理论,研究压电复合材料层合板在124内共振情况下的全局动力学.孙莹等1 5通过能量相位法研究热应力下环形桁架天线混沌运动.针对悬 索 这 类 典 型 的 工 程 柔 性 结 构,Z h a o等5,6基于局部分岔分析发现:温度变化会引起动力系统参数的微小改变,进而导致其整体非线性振动特性产生明显定性和定量的改变.倘若进一步采用全局分岔分析,可以从能量轨道的角度更好解释以及探究动力系统复杂的运动行为9.由此可见:一方面温度变化对支座位移运动引发的参数激励下悬索动力学行为影响有待进一步研究;另一方面基于全局分岔理论描述温度变化对系统动力学行为的影响更鲜有报道.因此本文将局部分岔分析6拓展到全局分岔分析,采用G a l e r k i n法和多尺度法并基于规范型理论和能量相位法,系统探究温度变化对悬索全局动力行为尤其是多脉冲轨道的影响,揭示温度变化下动力系统可能出现的大幅振动机理.1 运动方程及摄动分析1.1 数学模型如图1所示,悬索水平悬挂于O和B之间,以O为原点,建立直角坐标系O-x y,L和b分别代表悬索的跨度和初始垂度.本文考虑支座位移引发的共振响应,图中表示支座位移最大值,水平和垂直方向最大位移分量分别为x和y,支座运动A表示为:A(t)=c o s t,其中表示支座运动的频率.33动 力 学 与 控 制 学 报2 0 2 3年第2 1卷图1 悬索参数激励模型F i g.1 P a r a m e t r i ce x c i t a t i o nm o d e l o f s u s p e n d e dc a b l e s当环境温度发生整体改变时,基于增量热场理论,悬索将产生热应力构形2,3,此时张力改变,垂跨比随之变化,悬索的轴向和竖直方向位移分别用u(x,t)和v(x,t)表示.当不考虑温度变化影响时,利用H a m i l t o n变分原理,可得悬索面内运动微分方程1 6:m v t+cvv-H v+E A y+v()u+y v+12v 2=0(1)其中:点表示对t求导,撇表示对x求导,E A表示抗拉刚度,H表示初始水平张力,cv表示阻尼系数,悬索的静态构形用抛物线表示:y=4b x(L-x)/L2。在支座运动引发的参数激励作用下,悬索的总加速度v t由两部分组成:静态位移产生的加速度v 0以及动态应变产生的加速度v.由于垂度较小,分析时可以假设v 0是线性分布,此时v 0可表示为:v 0(x,t)=-(1-x)2yc o s t.为描述简便,引入以下无量纲参数:(x-,u-,v-)=(x,u,v)/L,cv=2m vt-=(t/L)H/m,-=L m/H=E A/H,=v-(2)其中:表示固有频率,表示静态拉伸系数,v示悬索垂向阻尼比.为表示简便,忽略上画线,可得支座运动激励下悬索面内运动方程:v+2v-v-(y+v)-xc o s t+10y v+12v 2dx=1-x()2yc o s t(3)已有研究表明5:温度变化会改变拉索张力以及垂跨比,其影响可通过引入张力改变系数2来表示.因此在无量纲方程的基础上,可得考虑温度变化时悬索受参数激励的运动微分方程为5,6:v+2v-2v-xc o s t+10y v 2+12v 2dx y 2+v=1-x()2yc o s t(4)采用G a l e r k i n法对方程进行离散:v x,t()=Nn=1qnt()nx()(5)其中:考虑温度影响下的模态函数n(x)参考文献5.本文考虑仅发生内共振的两个模态,令N=2,将方程(5)带入方程(4)中,方程两边同时分别乘模态函数n(x),并从0到1积分,可得两自由度的运动微分方程组:q 1+21q1+21q1+1 1q21+1 3q22+1 1q31+1 2q1q22=P1c o s t+K1c o s t q1(6)q 2+22q2+22q2+2 2q1q2+2 1q21q2+2 2q32=P2c o s t+K2c o s t q2(7)其中:相关系数如附录所示.1.2 摄动分析采用多尺度法对方程(6)、(7)进行求解,引入摄动参数以及三个时间尺度:T0=1t,T1=2t,T2=3t,将参数激励项和阻尼项摄动到二阶:Kn=2Kn,Pn=2Pn,1,2=21,2.外激励和固有频率间关系通过引入调谐参数1和2来表示:=21+1,=2+2.方程的解可以假设为qi(t,)=3j=1jqi,j,将qi(t,)带入式(6)、(7)中,并令两边i的系数相等.令qi(t,)两边摄动参数系数相等,求解可得qn,1的表达式:qn,i=An(T1,T2)einT0+A-n(T1,T2)e-inT0(n=1,2),其中A和A-为一对共轭复数.将qn,1带入qi(t,)中进行求解,通过消除长期项,可以得到方程(6)、(7)的平均方程:-2D2A1i1-21i1A1-31 1A21A-1+K1A-1ei 1T2/2-21 2A1A2A-2+1 021 1A1A1A-1321+41 11 3A1A2A-221-21 32 2A1A2A-221+212-21 32 2A1A2A-221-212=0(8)43第4期吴先强等:温度变化对悬索全局动力学特性影响P2ei2T2/2-22i2A2-2D2A2i2-22 1A1A-1A2-32 2A22A-2-1 32 2A2A2A-2422-21+21 12 2A1A2A-121+21 32 2A2A2A-221-22 2A1A-1A221+212-22 2A1A-1A221-212=0(9)引入直角坐标:Am(t)=pm(t)-iqm(t)eim/2,带入方程(8)、(9)中,分离虚部和实部,得直角坐标形式下的平均方程:p1=-1p1+(k-1/2)q1-1q1(p21+q21)-2q1(p22+q22)(1 0)q1=-1q1+(k+1/2)p1+1p1(p21+q21)+2p1(p22+q22)(1 1)p2=-2p2-2q2-3q2(p21+q21)-4q2(p22+q22)(1 2)q2=-2q2+2p2+3p2(p21+q21)+4p2(p22+q22)-f(1 3)其中:f=-P2/22,k=K1/41,其余系数如附录所示.方程(1 0)(1 3)确定了参数激励下悬索面内振动的相位和振幅,因此通过该方程组研究悬索在温度场中多脉冲混沌动力学行为.为了得到归一化方程,需对平均方程(1 0)(1 3)进行化简,使用M a p l e程序进行计算1 7,引入i=i(i=1,2),f=f以及线性变换:p1=(1-k+1/2)u,p2=2Ic o s,q1=1u+v,q2=2Is i n.方程(1 0)(1 3)可表示为可积的H a m i l t o n系统:u=v=H0/v+H1/v(1 4)v=u(2I 2-1)+1u3-21v=-H0/u-(H1/u+21v)(1 5)I=-22I-f2Is i n=H0+H1-22I(1 6)=2+3u2+42I-fc o s2I=-H0I-H1I(1 7)其中:-1=21+(k+1/2)(k+1/2-1).因此系统(1 4)(1 7)的哈密顿函数表示为:H0=v2/2-I2+3()u2+u2-1/2-1u4/4-2I-I24(1 8)H1=f2Ic o s(1 9)2 无扰动系统动力学假设系统(1 4)(1 7)中的扰动参数=0,那么余下系统中(1 4)、(1 5)不含,将I视作常数,因此可将系统对于(u,v)平面和(I,)平面进行解耦,得:u=vv=u(2I 2-1)+1u3(2 0)此时系统(2 0)的哈密顿函数简化为:H-u,v()=v2/2+u2/2-1u41/4(2 1)其中:=-1-2I 2.(a)稳定轨道:0(a)s t a b l eo r b i t s:0(b)异宿分岔:0,10(b)h e t e r o c l i n i cb i f u r c a t i o n:0,10图2 三维相空间流形结构图F i g.2 T h r e e-d i m e n s i o n a l p h a s es p a c em a n i f o l ds t r u c t u r ed i a g r a m s53动 力 学 与 控 制 学 报2 0 2 3年第2 1卷基于方程(2 1),绘制系统(2 0)在三维相空间中的流形结构图.如图2所示,当0,系统(2 0)将出现异宿分岔.当0,系统的唯一解为零解(u,v)=(0,0),该零解为鞍点.在曲线-1-2I 2上,即I=-1/22的零解通过叉型分岔由一个变为三个,不妨将三个解分别记为q0(I)=(0,0),q(I)=(D,0),由方程(2 0)的零解条件可以得到D:D=-1-22I()/11/2(2 2)观察图2(b)可知q0(I)为轨道中心,q(I)为轨道鞍点.由于异宿轨道可能产生混沌运动,故本文主要研究图2(b)这类情况.对于所有的I0,I1,方程(2 0)有一个鞍点q(I)和它自身的异宿轨道相连.四维空间(u,v,I,)由下式定义:M0=(u,v,I,)u=qI(),v=0,0II1,02(2 3)系统(2 0)异宿轨道的表达式:u(T1)=1t a n h()v(T1)=21s e c h2()(2 4)其中:=(2/2)T1.将无扰动系统限制在M0上再进行分析,对于任意I0,I1,满足DIHq0(I),I=0的值I被称为共振值,记作Ir,当I=Ir时可以得到此时连接异宿轨道的相位漂移角:=-2 2(3/1)(2 5)3 多脉冲混沌轨道由于不变流形M0在微小的非零扰动影响下是不动的,所以M0在小扰动的持续作用下变为M:M=(u,v,I,)u=q0(I),v=0,0Is,将Nf减去2k后重新定义Nf,验证降落点平移2k后是否还在鞍点s附近,如果Nf的能量大于鞍点s的能量,则:HD(0,Nf)HD(0,s)(3 7)c o sNf-c o ss2Ir(Nf-s)/f(3 8)当降落点落在吸引域之内,系统可能存在混沌运动.4 数值算例和分析本文的数值算例中,悬索的参数选取1 9:L=8 5 0 m,=8 3 7 6 k g/m3,E=1.7 9 41 01 1P a,A=0.1 1 5 9 m2,1=2=0.0 0 5,=1.21 0-5/,结合实际工程情况,温度变化 T=4 0 C.图3给出了温度变化以及I r v i n e参数(l o g2)对悬索的前六阶固有频率的影响.图中放大处,二阶正对称模态和一阶反对称模态频率之间呈现两倍关系,系统可能发生2:1内共振,温度变化下系统的各参数和系数大小如表1所示.图3 考虑温度变化影响下悬索的前六阶模态频率F i g.3 F i r s t s i xn a t u r a l f r e q u e n c i e so f s u s p e n d e dc a b l ec o n s i d e r i n gt e m p e r a t u r ec h a n g e s表1 考虑温度变化影响的悬索参数T a b l e1 P a r a m e t e r so f s u s p e n d e dc a b l e sc o n s i d e r i n gt e m p e r a t u r ec h a n g e s T/(C)l o g21/2/123405.0 8 12.0 0 04.0 2 51 0 0 4 2 1-3 2 8 3 7 8 6-3 3 2 5 9 0 51 6 7 9 0 1 4-4 04.9 8 81.9 6 53.9 7 11 0 0 0 2 4-2 3 6 9 0 9 5-2 4 1 0 6 1 91 2 2 1 1 7 64 05.1 9 22.0 3 34.1 1 01 0 3 4 6 9-2 6 2 7 1 4 9-2 6 7 0 1 1 91 3 5 1 4 1 4表2 考虑温度变化影响的L y a p u n o v指数和吸引子类型T a b l e2 L y a p u n o ve x p o n e n t sa n da t t r a c t o r s t y p e sc o n s i d e r i n gt e m p e r a t u r ec h a n g e s T/(C)1234A t t r a c t o r s t y p e s00.0 0 0 4 6-0.0 0 4 6 5-0.0 0 5 3 7-0.0 1 0 4 4P e r i o d i c-4 00.0 7 5 2 00.0 0 0 2 3-0.0 0 9 8 6-0.0 8 5 5 7C h a o t i c4 00.0 0 0 1 6-0.0 0 4 5 0-0.0 0 5 3 5-0.0 1 0 3 2P e r i o d i c 采用四阶龙格库塔法计算方程组(1 0)(1 3),并计算该四维系统的L y a p u n o v指数.假设支座横向振动幅值x=0.0 0 0 1,当调谐参数2=-0.1 0时,选取脉冲轨道数N=3.此时三组1均大于0,为了使系统在该环境下产生异宿轨道,必须让方程(2 2)中的D2为正数,即=-1-2I 20.图4表示不同温度场中悬索产生异宿轨道的条件.当 T=-4 0或0 C时,0的解集为闭集,而当 T=+4 0时,0的解集为开集.选取-4 0这组参数,当1在0.3 0 7 7 5,1.2 5 0 7 2 范围73动 力 学 与 控 制 学 报2 0 2 3年第2 1卷内,系统(2 0)产生异宿轨道.当1=0.4 8 3时,带入方程(3 4)得到满足能量差函数横截零点的激励f=1.4 7 2 41 0-4.此时共振值I=Ir=3.4 6 7 11 0-8,相位漂移值=2.0 8 6 2.通过求解方程(3 6),可得:Nf=6.2 4 0 9,s=3.1 5 9 5,c o s(Nf)-c o s(s)=1.9 9 8 9,2Ir(Nf-s)/f=0.0 2 7 6.系统参数满足方程(3 7)、(3 8),即可能出现混沌响应.图5分别表示 T=0,-4 0,4 0时,该动力系统的时程曲线、相位图和庞加莱截面.如图所示,温度变化 T为0以及4 0,图5(a)和5(c)时程曲线在一段时间运动后将逐渐趋于稳定,最大无量纲振幅减小约为0.0 0 0 2.二者相位图呈现为闭合圆环,庞加莱截面均表现为杂乱但有限的点集.由表2可知,这两种温度变化条件下,其L y a p u n o v指数均为(零,负,负,负),即可判定,当0及4 0时,系统吸引子类型均为周期型.图4 考虑温度变化影响的系统异宿轨道产生条件F i g.4 C o n d i t i o n s f o r t h eg e n e r a t i o no f s y s t e msh e t e r o c l i n i co r b i t sc o n s i d e r i n gt e m p e r a t u r ec h a n g e s图5 考虑温度变化影响的时程曲线,相位图,庞加莱截面:(a)0;(b)-4 0;(c)4 0F i g.5 T i m eh i s t o r yc u r v e,p h a s ed i a g r a m,P o i n c a r s e c t i o nc o n s i d e r i n gt e m p e r a t u r ec h a n g e s:(a)0;(b)-4 0;(c)4 0 然而,当 T=-4 0时,系数满足方程(3 7)-(3 8)条 件,即 图5(b)所 示 的 运 动 存 在 多 脉 冲S h i l n i k o v型轨道.对比图5中三组时程曲线,当 T=-4 0时,很长一段时间后运动仍没有趋于稳定,最大振幅约为0.0 0 1.此时图5(b)中相位图和庞加莱截面图均表现出系统存在无法预知的运动83第4期吴先强等:温度变化对悬索全局动力学特性影响行为.同理计算得到的-4 0的L y a p u n o v指数为(正,零,负,负),如表2所示,其中最大L y a p u n o v指数大于零,即判定该运动存在混沌吸引子,并且该混沌为多脉冲S h i l n i k o v型轨道混沌.由此可见,当 T=-4 0时,系统存在混沌吸引子,且伴随无序的大幅振动.此时动力系统的振动幅值难以准确预测,由此导致的大幅振动可能会影响结构的安全.当不考虑温度变化(0)或者变化升高(4 0)时,系统参数不再满足能量差函数的零解条件,S h i l n i k o v型多脉冲轨道消失,混沌运动也随之消失,最大振幅明显下降.由此可见温度变化有可能明显改变动力系统的吸引子类型以及系统的振动幅值.5 结论温度变化会直接改变悬索静态构形,导致其固有频率及非线性系数发生改变,影响频率间的公倍关系,进而引发系统全局动力学行为产生变化.当悬索存在S h i l n i k o v型多脉冲轨道并产生混沌运动时,温度变化会影响悬索系统中S h i l n i k o v型多脉冲轨道的产生条件.在异宿轨道参数限制下,相同的激励条件,温度变化可能会引起产生混沌响应的多脉冲轨道消失,导致动力系统展现出截然不同的动力学行为,改变系统的共振响应幅值,影响结构安全.参考文献1刘永健,刘江,张宁.桥梁结构日照温度作用研究综述 J.土木工程学报,2 0 1 9,5 2(5):5 9-7 8.L I U YJ,L I UJ,Z HAN GN.R e v i e wo ns o l a r t h e r-m a l a c t i o n so fb r i d g es t r u c t u r e sJ.C h i n aC i v i lE n-g i n e e r i n gJ o u r n a l,2 0 1 9,5 2(5):5 9-7 8.(i nC h i n e s e)2L E P I D IM,G A T TU L L IV.S t a t i ca n dd y n a m i cr e-s p o n s eo fe l a s t i c s u s p e n d e d c a b l e s w i t h t h e r m a le f f e c t sJ.I n t e r n a t i o n a l J o u r n a l o f S o l i d s a n dS t r u c-t u r e s,2 0 1 2,4 9:1 1 0 3-1 1 1 6.3T R E Y S S E D EF.F i n i t e e l e m e n tm o d e l i n go f t e m p e r-a t u r e l o a de f f e c t so nt h ev i b r a t i o no f l o c a lm o d e si nm u l t i-c a b l es t r u c t u r e sJ.J o u r n a l o fS o u n da n dV i-b r a t i o n,2 0 1 8,4 1 3:1 9 1-2 0 4.4MAL,X U H,MUNKH B AA T A RT,e t a l.A na c-c u r a t ef r e q u e n c y-b a s e dm e t h o df o ri d e n t i f y i n gc a b l et e n s i o nw h i l e c o n s i d e r i n ge n v i r o n m e n t a l t e m p e r a t u r ev a r i a t i o nJ.J o u r n a l o fS o u n da n dV i b r a t i o n,2 0 2 1,4 9 0:1 1 5 6 9 3.5Z HAO Y,P E NGJ,Z HAO Y,e ta l.E f f e c t so ft e m p e r a t u r ev a r i a t i o n so nn o n l i n e a rp l a n a rf r e ea n df o r c e do s c i l l a t i o n s a t p r i m a r y r e s o n a n c eo f s u s p e n d e dc a b l e sJ.N o n l i n e a rD y n a m i c s,2 0 1 7,8 9:2 8 1 5-2 8 2 7.6林恒辉,赵珧冰.温度变化对悬索非线性内共振响应特性影响 J.振动与冲击,2 0 2 1,4 0(8):1 6 5-1 7 2.L I N H H,Z HAO Y B.I n f l u e n c e so ft e m p e r a t u r eo ns u s p e n d e dc a b l e s n o n l i n e a rv i b r a t i o nc h a r a c t e r i s-t i c sc o n s i d e r i n gi n t e r n a lr e s o n a n c e sJ.J o u r n a lo fV i b r a t i o na n dS h o c k,2 0 2 1,4 0(8):1 6 5-1 7 2.(i nC h i-n e s e)7李凤臣,杨鸥,田石柱,等.考虑前2阶模态组合的拉索非线性参数共振研究 J.防灾减灾工程学报,2 0 1 5,3 5(2):2 4 9-2 5 5.L IFC,YAN GO,T I ANSZ,e t a l.N o n l i n e a rp a r-a m e t r i cr e s o n a n c e