含外激励van der Pol-Mathieu 方程的非线性动力学特性分析1)

黄建亮 王 腾 陈树辉

(中山大学应用力学与工程系,广州 510275)

引言

在工程中存在着很多可以用自激振动和参数激振联合作用的van der Pol-Mathieu 方程来描述的振动,例如,含有万向接头的转子系统的横向振动[1],卡盘作业过程中的参数激振[2],含有自振和参数激振的齿轮装置系统的振动[3],尘埃等离子体中的颗粒电荷的动力学行为[4],高层建筑结构在风荷载下的振动[5-6]等,都是可用van der Pol-Mathieu 方程来描述振动的典型例子.van der Pol-Mathieu 方程同时含有自激振动和参数激振,蕴含着丰富的动力学行为,多年来一直是众多学者关注点之一.

Tondl[7]首先分析了van der Pol-Mathieu 方程中自激振动和参数激振的相互作用,并在共振区域发现了周期响应.Kotera 和Yano[8]用两个频率的和分析了van der Pol-Mathieu 方程在参数共振区域的近似一阶和二阶的周期解.陈予恕和徐鉴[9]研究了van der Pol-Duffing-Mathieu 型系统主参数共振分岔解,得到该非线性参数激励系统依赖于物理参数变化的振动模式.Szabelski 和Warminski[10]分析了自激振动和参数激励对van der Pol-Mathieu 方程的影响,并且研究了附加外激励在同步区域内动力学行为的影响.Warminski 等[3]对含有自激振动和参数激励的两自由度系统进行分析,并得到了不同类型的响应,包含有周期响应,准周期运动响应和混沌.彭献和陈自力[11]引入参数变换,将强非线性系统转化为弱非线性系统,利用摄动思想分析得到了van der Pol-Mathieu方程的1/2 亚谐共振周期解.Belhaq 和Fahsi[12]和Pandey 等[13]分析了van der Pol-Mathieu-Duffing 系统的响应,得到该类系统可含有1:1 锁频,2:1 次谐波锁频和准周期运动响应.张琪昌等[14]利用改进的类Pad´e 方法计算了van der Pol-Duffing 方程的同异宿轨道.Warminski[15]研究了含有van der Pol 和Rayleigh函数在两个不同自激振动模型下具有时滞状态的自激振动,参数激振和强迫振动作用下的相互作用.许多学者也对各类含有参数激振的非线性系统进行研究[16-19],得到了不同的非线性振动特性和运动分岔.

早期对于van der Pol-Mathieu 方程的众多研究主要集中在含有一个基频的周期响应及其稳定性分析.可以利用不同的摄动方法求得这类方程的近似解析解[20].然而,对于van der Pol-Mathieu 方程来说,因有自激振动与参数激励振动的相互耦合作用,使得系统不仅有周期响应,而且还有准周期运动响应,甚至产生混沌,近年来受到众多学者的关注.Belhaq 和Houssni[21]为了构造准周期运动的近似解,提出了双摄动的思想,其方法包含了两个步骤:第一步利用广义的平均法将准周期系统变为周期减化系统;第二步是用多尺度法对周期减化系统构造出近似的准周期运动解.他们利用双摄动方法得到了同时含有二次和三次非线性项的参数激励和外激励的单自由度系统的准周期解.该双摄动方法的本质就是对周期减化系统在平衡点附近的周期解进行非线性近似,该方法可进一步推广到各类非线性系统的准周期运动分析中[5-6,12,22].Fan 等[23]也利用了双摄动方法分析了van der Pol-Mathieu 方程在有外激励和无外激励两种情况下的周期解和准周期运动近似解的包络线.然而,双摄动方法只能得到准周期运动包络线的最大和最小幅值,无法得到系统准周期运动的具体响应情况,更无法得到系统准周期运动的各个响应频率.Warminski 等[3]和Warminski[15]利用多尺度法分别分析了含有自激励和参数激励的两个自由度时滞系统,并利用数值法得到了两个时滞系统的准周期运动响应和混沌.Veerman 和Verhulst[24]利用平均法分析了van der Pol-Mathieu 方程,并得到了由1 阶和1/ε 阶基础周期构造而成的准周期运动响应,其中ε 是小量.上述的各种摄动法只能得到van der Pol-Mathieu 方程准周期运动近似解,有些方法得到的结果只能描述准周期运动的最大和最小振幅,特别是在邻近分岔点处用摄动法得到的准周期运动解与数值解相差甚大,据作者所知,迄今为止尚未有有效的摄动法能精确地计算并得到此类van der Pol-Mathieu 方程的精确准周期运动解.

本文针对含有外激励的van der Pol-Mathieu 方程进行研究,主要是发现了单自由度的van der Pol-Mathieu 方程准周期运动的频谱含有均匀边频带的新特性。此新特性与之前研究分析的多自由度非线性系统中内共振引起的准周期运动的频谱特性[25-27]相类似,即都含有均匀的边频带,不同之处在于多自由度非线性系统中的准周期运动是由于不同振动模态之间在内共振条件下相互作用产生的,而单自由度的van der Pol-Mathieu 方程的准周期运动是由于自激振动与参数激振耦合产生的,仅是van der Pol 方程中的自激振动或仅是Mathieu 方程中的参数激振并不能产生准周期运动.根据此准周期运动频谱含有均匀边频带的特性,它包含了两个基频,一个是已知激励频率ω;另一个是事先未知的频率ωd,即为边频带中相邻两个频率之间的距离,那么准周期运动中所有的频率成份都可表示为这两个基频的线性组合.因此,本文利用传统的增量谐波平衡法(IHB 法)分析单自由度含有外激励的van der Pol-Mathieu 方程的周期响应,并推广两时间尺度的IHB 法,其中一个时间尺度是快时间尺度τ1=ωt;另外一个是慢时间尺度τ2=ωdt(ω ≫ωd),应用于分析此van der Pol-Mathieu方程的准周期运动响应.

1 含外激励van der Pol-Mathieu 方程的周期解及其稳定性

对于含有外激励van der Pol-Mathieu 方程,有3种激励共同作用,即自激励,参数激励和外激励,可描述为下列的微分方程

其中,x是因变量,t是时间,k1,k2和k3是常数,是线性固有频率,ω 是激励频率,f是外激励幅度.方程(1)包含了自激励,即van der Pol 项−(k1−k2x2)dx/dt;参数激励,即Mathieu 项k3cos(2ωt)x和简谐外激励fcos ωt.对于此方程,首先利用传统的IHB 法确定其周期解,然后利用Floquet 理论,并结合精细的Hsu 法判断其周期解的稳定性及分岔,最后得到了Saddlenode 和Hopf 这两种分岔类型.

1.1 传统的IHB 法

对于含外激励的van der Pol-Mathieu 方程,可以利用传统的IHB 法来确定其周期解.引入一个新的时间变量

方程(1)可变为

其中,′表示对时间τ 的导数.

传统的IHB 法包含两个过程,增量过程和谐波平衡过程.增量过程即Newton-Raphson 迭代过程,是对微分方程(3)进行线性化.设x0,ω0,f0是方程(3)的解,则其邻近点可表示为

式中,∆x,∆ω,∆f为增量,将式(4)代入方程(3),并略去高阶小量后可得到以∆x,∆ω,∆f为未知量的增量方程

为不平衡力.如果x0,ω0,f0是方程(3)的准确解,那么R=0.

传统 IHB 法的第二个过程是谐波平衡过程(Galerkin 过程).因方程(1)是自激励和参数激励的系统,包含了奇次非线性项,所以对于其周期解,可设

这里,ak和bk是傅里叶系数,nc和ns分别是cosine和sine 谐波项的项数

对式(7)和式(8)进行微分,得

将式(7),式(8),式(12)和式(13)代入方程(5),并利用Galerkin 过程平衡谐波项,得

积分上式并整理归并为以∆A,∆ω 和∆f为未知量的代数方程组

这里,上标T表示为矩阵的转置,KA是n×n的矩阵,和Rω是n×1 的矩阵,n=nc+ns.本文只考虑在某一固定外激励幅值下的频率响应曲线,那么f取固定值,∆f=0,于是方程(15)成为

用传统的IHB 法求解时,可事先给定一个初值,然后利用增量过程和谐波平衡过程追踪出所有的解.在求解过程中,可采用振幅增量,频率增量或弧长增量,具体可参见Cheung 等[28],陈树辉[29]的研究结果.

1.2 周期解的稳定性与分岔

利用传统的IHB 法确定含外激励van der Pol-Mathieu 方程的周期解后,需要分析其稳定性.设x0为已求得的解,给x0一个小的扰动∆x

将式(21)代入方程(1),略去高阶小量,并注意x0满足方程(1),可得

方程(22)称为扰动方程,即从已知的平衡位置扰动而得的方程.方程(22)的稳定性特征可以用Floquet理论来分析.将方程(22)重新写为状态方程

由于x0是τ 的周期为T=2π 的函数,所以Q21和Q22也是周期为T=2π 的函数.

对于方程(23),存在着一组基础解系

这一基础解系可用矩阵来表示

可以看到,Y满足矩阵方程

由于Q(τ+T)=Q(τ),Y(τ+T)也是基础矩阵解,因此,它可表示为

式中,P为非奇异常数矩阵,称为转移矩阵.根据Floquet 理论,方程(23)的稳定性准则与矩阵P的特征值λi有关.如果转移矩阵P的所有特征值的模都小于1,则方程(23)的运动是有界的,因而此周期解是稳定的;如果转移矩阵P中至少有一个特征值的模大于1,则方程(23)的运动是无界的,此周期解也就不稳定.

采用数值方法求解方程(30)的基础解系,其中有效地计算转移矩阵P是稳定性分析的关键.Hsu[30-31]和Hsu 和Cheng[32]提出一个近似求解转移矩阵P的有效方法,其主要思想是把一个周期等分为若干时间间段,在每一时间段上进行积分.Friedmann 等[33]给出了该方法的具体推导表达式.Huang等[34]在此基础上,结合了精细积分算法[35],提出了精细的Hsu 法,有效地减小了计算转移矩阵P的舍入误差.

1.3 含外激励van der Pol-Mathieu 方程周期解的响应特性

在用传统的IHB 求解过程中,取谐波项项数nc=ns=3,式(7)可写为

图1 所示为当k1=0.01,k2=0.01,k3=0.05 和f=0.01 时van der Pol-Mathieu 方程周期响应的频率响应曲线ω-A1,高阶的谐波项较小而没有在此显示,图1(b)是图1(a)中标示区域的放大图,其中实线表示稳定的周期解,虚线表示不稳定的周期解,实圆点表示Saddle-node 分岔,实三角形点表示Hopf 分岔.图1 中4 个临近分岔点的响应点的Floquet 乘子如表1 所示.从表1 可以看出,在复平面上其中的一个Floquet 乘子从+1 方向穿出单位圆,导致了Saddlenode 分岔(S1和S2),从而出现跳跃现象;而一对共轭的Floquet 乘子穿出单位圆,导致了Hopf 分岔(H1和H2),从而出现了准周期运动.

图1 含外激励van der Pol-Mathieu 方程周期响应的频率响应曲线,ω-A1其中k1=0.01,k2=0.01,k3=0.05 和f=0.01,(b)为(a)中标示区域的放大图Fig.1 Frequency response curve ω-A1of periodic response of the van der Pol-Mathieu equation with external excitation with k1=0.01,k2=0.01,k3=0.05,and f=0.01,(b)is an enlarged view of a zone highlighted in(a)

在含有外激励van der Pol-Mathieu 方程中,由于非线性的影响,周期解的频率响应曲线出现了多值性,从而出现了一些重要而且有趣的动力学现象,如图1 所示的跳跃现象.假设外激励幅值f不变而频率ω 慢慢地变化,先考察向右扫频的过程,周期解从分岔点H1出发,当频率ω 逐渐增大时,振幅A1开始变大,在ω=1.0 附近达到最大值,过了最大值后振幅A1变小到达分岔点S1,如果ω 继续增大,则振幅A1突然从点S1跳跃到点P1(如图1(a)所示),然后沿曲线逐渐变小,最后到达点H2.反之,考察向左扫频的过程,周期解从分岔点H2出发,当频率ω 逐渐减小时,振幅A1开始变大,在ω=1.0 附近也达到最大值,过了最大值后振幅A1变小到达分岔点S2,如果ω 继续减小,则振幅A1突然从点S2跳跃到点P2(如图1(b)所示),然后沿曲线逐渐变小,最后到达点H1.

表1 频率ω 在4 个分岔点S1,S2,H1和H2临近点响应的Floquet 乘子λ1和λ2,其中=√Table 1 Floquet multipliers λ1and λ2for the frequency ω near the four bifurcation points S1,S2,H1,and H2,where=√

表1 频率ω 在4 个分岔点S1,S2,H1和H2临近点响应的Floquet 乘子λ1和λ2,其中=√Table 1 Floquet multipliers λ1and λ2for the frequency ω near the four bifurcation points S1,S2,H1,and H2,where=√

图2 ω=1.0 时含外激励van der Pol-Mathieu 方程含有3 个周期解,其中k1=0.01,k2=0.01,k3=0.05 和f=0.01:(a)第一个稳定周期解的时间历程图;(b)第二个稳定周期解的时间历程图;(c)3 个不同周期解的相图;(d)两个稳定周期解的频谱图Fig.2 Three different periodic solutions of the van der Pol-Mathieu equation with external excitation for ω=1.0 with k1=0.01,k2=0.01,k3=0.05,and f=0.01:(a)Time history of the first stable periodic response;(b)time history of the second stable periodic response;(c)phase plane portraits of the three solutions;(d)Fourier spectra of the two stable periodic response

在图1 中从点S1到S2的区域,存在着两个稳定的周期解和一个不稳定的周期解.图2 所示为频率ω=1.0 时van der Pol-Mathieu 方程含有3 个不同解的周期响应,其中图2(a)和图2(b)是两个稳定周期解的时间历程图,其中实线表示用Runge-Kutta(RK)数值法求得的解,圆圈表示用传统的IHB 法得到的结果,可以看出,两个方法得到的结果非常吻合;图2(c)所示为3 个解的相图,可以看出3 条都是封闭的曲线,进一步表明系统的响应为周期响应;图2(d)所示为两个稳定周期解的频谱图,可以看出|A1| ≫|A2| ≫|A3|,系统的响应主要以第一阶响应为主,其频率是参数激励频率2ω 的一半,与外激励频率一样,因此系统的响应为自激振动,1:2 参数激振和外激励的联合作用.

2 含外激励van der Pol-Mathieu 方程的准周期运动

对于含外激励van der Pol-Mathieu 方程的周期响应,其频率成份之间是可约的,例如在本文中,系统的频率成份为ω,3ω 和5ω.在上节分析得知,系统的周期解经Hopf 分岔后会产生准周期运动.对准周期运动来说,其频率成份至少含有两个或以上的不可约频率.事实上,本文首次发现了van der Pol-Mathieu 方程的准周期运动的新特性,即其频率成份是由落在ω,3ω 和5ω 附近的边频带组成,并且这些边频带里的频率是等相距的,即为均匀的边频带.此准周期运动的新特性,正是由于van der Pol-Mathieu 方程中的自激振动和参数激振相互作用产生的.如若只考察van der Pol 方程中的自激振动或只考察Mathieu 方程的参数激振,并不会产生准周期运动.根据此准周期运动的新特性,其频谱中含有两个基频,一个是已知的频率ω,另一个可把它看作是均匀边频带里相邻频率的距离ωd,且ωd事先是未知的,此时准周期运动中所有的频率都可由这两个基频线性组合而成.根据此特性,发展了传统的IHB 法,引入两个时间尺度,使之适用于精确求解含外激励van der Pol-Mathieu 方程的准周期运动,并且可以得到此准周期运动所有的频率成份.

2.1 两时间尺度的IHB 法

相较于传统的IHB 法,两时间尺度的IHB 法在求解含外激励van der Pol-Mathieu 方程的准周期运动过程中有3 点改进之处,具体推导如下.

(1)第一点改进是引入两个时间变量代替式(2)中针对于周期解的一个时间变量τ=ωt,如下那么x(t)是τ1和τ2的函数,即x(t)=x(ωt,ωdt)=x(τ1,τ2).记

同样,利用增量过程,对于第一步的两时间尺度IHB法,设x0,ω0,ωd0,f0是方程(36)的解,则其邻近点可表示为

将式(37)代入方程(36),并忽略高阶小量,得到以∆x,∆ω,∆ωd,∆f为未知量的增量方程

(2)第二点改进是将x0展开为含有两个时间变量τ1和τ2的二重傅里叶级数形式

式中,Nm是傅里叶级数中保留到最高阶谐波项的项数,和是谐波项系数.在含外激励van der Pol-Mathieu 方程的准周期运动中,发现了其频谱的边频带是集中在频率ω 及其奇数倍频上,即ω,3ω,5ω,···,而且边频带内相邻频率之间的距离是相等的,也就是说,所有含外激励van der Pol-Mathieu方程的准周期运动的频率成份可表示为ω,ω±ωd,ω±2ωd,···,ω ±m1ωd,3ω,3ω ± ωd,3ω ± 2ωd,···,3ω ±m2ωd,···,(2nc−1)ω,(2nc−1)ω ± ωd,(2nc−1)ω ±2ωd,···,(2nc−1)ω±mncωd,其中2m1+1,2m2+1,···,2mnc+1 是对应于ω,3ω,···,(2nc−1)ω 边频带中谐波项的项数.那么式(41)可表示为

2.2 准周期运动的频率响应曲线

图3 含外激励van der Pol-Mathieu 方程周期响应的频率响应曲线ω-A1和准周期运动(QP)的频率响应曲线ω-1,−1,ω-1,0和ω-1,1,其中k1=0.01,k2=0.01,k3=0.05 和f=0.01,(b)为(a)中标示区域的放大图Fig.3 Frequency response curves ω-A1of periodic response and ω-1,−1,ω-1,0,and ω-1,1of quasi-periodic(QP)motion of the van der Pol-Mathieu equation with external excitation with k1=0.01,k2=0.01,k3=0.05,and f=0.01,(b)is an enlarged view of a zone highlighted in(a)

2.3 不同频率ω 点的准周期运动

为了考察含外激励van der Pol-Mathieu 方程的准周期运动特性,在图3 中选取3 个不同频率ω 点来分析,第一个点是ω=0.987 35,即在左侧曲线上临近于分岔点H1;第二个点是ω=0.985,即在左侧比第一个点较远离分岔点H1;第三个点是ω=1.02,即在右侧曲线上的一点.图4 至图6 所示分别为在图3 中上述3 个点的时间历程图,频谱图和庞加莱截面图.表2 所示为利用两时间尺度IHB 法求出在图3 中3 个点ω=0.987 35,0.985,1.02 处的ωd和从表2 可以看出,在系统的准周期运动中,当频率ω越远离Hopf 分岔点时,均匀边频带的频率间的间距ωd越大.

图4(a)所示为在分岔点H1附近频率ω=0.987 35 时分别用四阶的RK 数值法和用含有350个谐波项的两时间尺度IHB 法计算得到的含外激励van der Pol-Mathieu 方程准周期(QP)运动的时间历程图,其中k1=0.01,k2=0.01,k3=0.05 和f=0.01.可以看出,准周期运动具有‘拍’ 现象,且有明显的包络线.图4(b)是图4(a)中标示区域的放大图,可以看出不管是在振幅较大的区域,还是在如图4(a)中振幅较小的区域,两时间尺度IHB 法得到的结果与RK 法得到的结果完全一致.图4(c)和图4(d)是频率ω=0.987 35 点上的准周期运动的频谱图,其中图4(d)是在频率ω 附近的频谱放大图.从图4(c)可以看出,频率成份主要集中在ω,3ω 和5ω 附近,且集中在ω 的幅值远大于其他幅值;在图4(d)上标示出在靠近ω 的5 个频率及其幅值,可以得到相邻频率间的间距为ωd=8.846 28×10−4,最大的幅值在频率为ω+ωd=0.988 24 上.

图4 在分岔点H1附近频率ω=0.987 35 时含外激励van der Pol-Mathieu 方程的准周期运动,其中k1=0.01,k2=0.01,k3=0.05 和f=0.01Fig.4 Quasi-periodic motion of the van der Pol-Mathieu equation with external excitation with k1=0.01,k2=0.01,k3=0.05,and f=0.01 at the freqeuncy ω=0.987 35 near the bifurcation point H1

图4 在分岔点H1附近频率ω=0.987 35 时含外激励van der Pol-Mathieu 方程的准周期运动,其中k1=0.01,k2=0.01,k3=0.05 和f=0.01(续)Fig.4 Quasi-periodic motion of the van der Pol-Mathieu equation with external excitation with k1=0.01,k2=0.01,k3=0.05,and f=0.01 at the freqeuncy ω=0.987 35 near the bifurcation point H1(continued)

图5 频率ω=0.985 时含外激励van der Pol-Mathieu 方程的准周期运动,其中k1=0.01,k2=0.01,k3=0.05 和f=0.01Fig.5 Quasi-periodic motion of the van der Pol-Mathieu equation with external excitation with k1=0.01,k2=0.01,k3=0.05,and f=0.01 at the parametric excitation freqeuncy ω=0.985

图6 频率ω=1.020 时含外激励van der Pol-Mathieu 方程的准周期运动,其中k1=0.01,k2=0.01,k3=0.05 和f=0.01Fig.6 Quasi-periodic motion of the van der Pol-Mathieu equation with external excitation with k1=0.01,k2=0.01,k3=0.05,and f=0.01 at the freqeuncy ω=1.020

表2 利用两时间尺度IHB 法求出在图3 中三个点ω=0.987 35,0.985,1.02 处的ωd和1,−1,1,0,1,1Table 2 A prior unknown ωdand three amplitudes 1,−1,1,0,and1,1in frequency response curves of the van der Pol-Mathieu equation with external excitation that are calculated by the IHB method with two time-scales at the three points ω=0.987 35,0.985,1.02 in Fig.3

表2 利用两时间尺度IHB 法求出在图3 中三个点ω=0.987 35,0.985,1.02 处的ωd和1,−1,1,0,1,1Table 2 A prior unknown ωdand three amplitudes 1,−1,1,0,and1,1in frequency response curves of the van der Pol-Mathieu equation with external excitation that are calculated by the IHB method with two time-scales at the three points ω=0.987 35,0.985,1.02 in Fig.3

图4(e)是取不同的谐波项的两时间尺度IHB 法和采用RK 法得到的频率ω=0.987 35 时含外激励van der Pol-Mathieu 方程准周期运动的庞加莱截面图.图中所示为RK 法得到的和取不同谐波项的两时间IHB 法得到的庞加莱截面图,可以看出随着m1越大,其结果是收敛的,且越接近于RK 法的结果;当取到m1=70,m2=15,m3=1 时(即谐波总项数˜n=350),两时间尺度IHB 法得到的结果与RK 法得到的结果相吻合,这也说明了含外激励van der Pol-Mathieu 方程在频率ω=0.987 35 点上的准周期运动至少含有175 个频率.表3 所示为用RK 法和取不同谐波项项数的两时间尺度IHB 法得到在ω=0.987 35 处的Amax和Amin的对比,其中Amax和Amin分别是准周期运动包络线幅值的最大值和最小值.从表3 可以看出,随着谐波项项数的增加,两时间尺度IHB 的结果越接近于RK 法的结果,进一步说明了在某些频率点上,两时间尺度IHB 法需要取足够的谐波项项数,才能得到精确的结果.

表3 用RK 法和取不同谐波项项数的两时间尺度IHB 法得到在ω=0.987 35 处的Amax和Amin的对比,其中Amax和Amin分别是准周期运动包络线幅值的最大值和最小值Table 3 Comparison of Amaxand Aminfrom the RK method and the IHB method with two time-scales with different harmonic terms at the point ω=0.987 35,where Amaxand Aminare the maximal and minimal value of envelops of the quasi-periodic motions,respectively

图5 和图6 分别是在另外两个点ω=0.985和ω=1.02 上的准周期运动响应情况,得到与点ω=0.987 35 上相似的准周期运动特征,其中利用两时间尺度IHB 法得到的结果所用到的谐波项项数为=70(m1=10,m2=5,m3=1).图5(a)和图6(a)的左部分是在[0,2000]紧缩时间尺度上的时间历程图,具有明显的包络线,右部分是在[2960,3000]扩张时间尺度上的时间历程图,从图中可以看出两时间尺度IHB 法得到的结果与RK 法得到的结果完全一致.

观测在不同频率ω 点的准周期运动,可以得到一些有趣的动力学现象.比较3 个点准周期运动的时间历程图4(a),图5(a)和图6(a)可以得知,准周期运动的包络线似乎都有‘周期’,且其‘周期’长度恰似为2πω/ωd,即分别为7012.78,758.17 和414.27,说明越接近分岔点的准周期运动的包络线的‘周期’长度越长;比较3 个点准周期运动的频谱图4(c),图5(b)和图6(b)可以得知,越远离分岔点的准周期运动的边频带的带宽越宽且其频率分布越稀疏,此原因是由于边频带相邻频率间的间距变大;比较3 个点准周期运动的庞加莱截面图4(e),图5(d)和图6(d)可以得知,利用两时间尺度IHB 法需要不同的谐波项项数以满足计算的精度,如在频率ω=0.987 35 点上需要350 个谐波项,而在其他两个点上只需70 个谐波项,表明了在不同点上的准周期运动含有不同的频率数.

3 结论

本文基于含外激励van der Pol-Mathieu 方程的非线性系统,针对系统在自激励,参数激励和外激励3种激励共同作用下,利用传统的IHB 法和两时间尺度IHB 法分别分析了系统的周期响应和准周期运动响应,得到的结果与数值RK 法得到的结果完全吻合,并得到以下结论.

(1)在含外激励van der Pol-Mathieu 方程的周期响应中含有两种类型的分岔,Saddle-node 分岔和Hopf 分岔.Saddle-node 分岔会导致跳跃现象,即系统响应从一个稳态的周期响应跳跃到另一个稳态的周期响应;Hopf 分岔会导致准周期运动,即系统响应从周期响应变化到准周期运动响应.另外,在两个周期解的Saddle-node 分岔点区域间含有两个稳定的周期解.

(2)含外激励van der Pol-Mathieu 方程中,由于自激振动和参数激振联合作用下,会产生准周期运动响应,并发现其新特性.在准周期运动的频谱中，频率ω,3ω,5ω 附近含有边频带,且边频带是均匀的，即边频带中相邻频率的间距ωd是一样的；在准周期运动的时间历程图中，准周期的包络线也有‘周期’,其‘周期’为2πω/ωd.

(3)在不同点上的准周期运动响应有不同的动力学特征.越靠近Hopf 分岔点的准周期运动响应,其边频带的带宽越小,且其频率分布越密集,同时其包络线的‘周期’越长.另外,在不同点上的准周期运动响应所含的频率数目也不同.

此外,针对准周期运动响应含有大量的频率成份,如何有效提高两时间尺度IHB 法的计算效率将是进一步的研究工作.