振动落砂机系统的拟周期碰撞设计*

伍　新，文桂林，何莉萍，徐慧东，魏克湘

(1.湖南大学 汽车车身先进设计制造国家重点实验室，湖南 长沙　410082；2.湖南工程学院 机械工程学院，湖南 湘潭　411104)



振动落砂机系统的拟周期碰撞设计*

伍新1,2，文桂林1†，何莉萍1，徐慧东1，魏克湘2

(1.湖南大学 汽车车身先进设计制造国家重点实验室，湖南 长沙410082；2.湖南工程学院 机械工程学院，湖南 湘潭411104)

摘要：首先建立了落砂机系统周期运动的Poincaré映射，考虑到在设计过程中经典的Neimark-Sacker分岔临界准则需要直接计算特征值带来的局限性，利用不直接依赖于特征值计算的显式临界准则，获得了系统发生Neimark-Sacker分岔的两参数区域图，所获得的参数区域图有助于主动设计系统的拟周期碰撞运动.然后应用中心流形-正则形方法进一步分析了拟周期碰撞运动的稳定性.最后数值仿真表明在选定的系统参数处能产生稳定的拟周期碰撞运动.

关键词：落砂机；冲击振动；Neimark-Sacker分岔；拟周期碰撞运动；稳定性

碰撞振动在实际工程领域中普遍存在，由于碰撞和冲击过程中固有的不连续性造成的强非线性，使得系统的动力学响应十分复杂多变，产生丰富的非线性现象，如分岔和混沌现象等[1-6].

振动落砂机是一种利用碰撞振动原理对砂箱进行落砂的机械设备.这种周期碰撞的机械设备工作频率单一，导致生产效率低、能耗高.而拟周期碰撞是典型的非线性振动，与简单的周期碰撞相比，拟周期碰撞具有多频性，以及振动加振荡之复式激振品质，可提高系统的功能和效率，同时又没有混沌运动的不可预测性和初始条件敏感性，该特性在实际工程领域具有应用潜力，能有效解决上述存在的一些缺陷.

近来，国内外一些学者对振动落砂机系统非线性特性开展了理论和数值模拟分析.罗冠炜等[7]通过理论分析和数值仿真揭示了振动落砂机周期运动经概周期分岔和倍周期分岔通向混沌的演化过程.丁旺才等[8]使用中心流形-范式方法研究了振动落砂机系统在强共振下的两参数开折的局部动力学行为并通过数值仿真进一步揭示了系统共振点附近的Hopf分岔不变环面和次谐分岔4-4周期运动.

随着分岔理论和非线性动力学设计方法的发展，人们开始关注如何主动利用分岔特性来提高系统的功能与效率，通过主动选择系统的参数来设计出具有所期望特性的分岔解.然而上述文献中对振动落砂机系统的非线性现象的研究大部分是基于特征值的特性来描述的传统的分岔准则.对于一个四维多参数的落砂机碰撞振动系统，如果按照传统的Hopf分岔临界准则，逐点试算特征值是否满足分岔的存在条件，这对于通过设计来主动实现分岔解具有一定的局限性.针对传统分岔准则的不足，文桂林等[9]提出了新的离散系统Hopf分岔准则，建立的Hopf分岔准则是由一些系统参数构成的代数等式和不等式组成的显式分岔临界准则，并不依赖于特征值的计算，这更适合分岔参数机理分析和设计.

本文以振动落砂机系统为研究对象，针对传统的映射Hopf分岔临界准则在主动设计方面存在的局限性，基于Neimark-Sacker分岔(二次Hopf分岔)理论，使用不直接依赖于特征值计算的显式临界准则来设计系统的参数，使其产生拟周期碰撞运动.然后，应用中心流形-正则形方法来分析拟周期碰撞运动的稳定性.最后，通过选取适当的系统参数，数值实现了系统稳定的拟周期碰撞运动.

1力学模型及其运动方程

对于振动落砂机系统的设计，仅考虑垂直方向振动，将振动落砂机系统简化为图1所示的两自由度的质量弹簧阻尼器系统.其中，质量块M，m分别表示质量为M的振动基座和质量为m的砂箱(包括型砂和铸件)，振动基座和基础之间用刚度为K的线性弹簧和阻尼系数为C的阻尼器连接，振动基座受到简谐力Fsin(ωt+δ)的作用.图中X，Y分别表示基座和砂箱的位移.砂箱和基座不发生碰撞时，砂箱只受重力作用，当基座和砂箱位移相同即X=Y，并且相对速度不为零时，它们会发生垂直方向的正碰.为了方便计算，振动落砂机系统采用文献[10]中的无量纲运动微分方程来描述：

图1　振动落砂机的力学模型

(1)

(2)

2振动落砂机系统周期运动的Poincaré映射

由无量纲方程(1)和(2)可知，在系统未发生碰撞时，基座和砂箱会遵循方程(1)做连续运动.当发生碰撞时，基座和砂箱的速度会遵循冲击方程(2)发生突变，得到下一次做连续运动的初始值，因此系统会进行碰撞、连续运动、再碰撞的循环运动.为了使系统产生Neimark-Sacker分岔，求得一个周期冲击运动.其无碰撞部分的解析表达式如下[8]：

xp(θ)=e-ζ(θ/z)(b10cosη(θ/z)+b20sinη(θ/z))+

sin(θ+τ0)，

(3)

(b20ζ+b10η)sinη(θ/z)]+cos(θ+τ0)，

(4)

yp(θ)=b30+b40θ-e1θ2/2，

(5)

(6)

表达式(3)～(6)中，确定的周期运动的相位角和积分常数如下

(7)

(8)

b30=b10+sin(τ0)，b40=e1π.

(9) 其中，s=e-ζ(2π/z)sinη(2π/z),c=e-ζ(2π/z)cosη(2π/z).由表达式(7), 可知周期运动需满足下面的一个条件：

(10)

基于获得的周期解(3)～(6)，选择系统的Poincaré截面

S=R(mod2π).

这里，R表示实数，可建立如下的Poincaré映射

Xk+1=F(μ,Xk)=

(11)

3振动落砂机系统的拟周期碰撞运动

3.1振动落砂机系统Neimark-Sacker分岔的显式临界条件

为了设计出振动落砂机系统拟周期碰撞运动，主要任务就是确定适当的系统参数，使系统(11)发生Neimark-Sacker分岔.如果采用传统的分岔临界准则，需要在参数空间内通过逐点取值来计算和验证系统的特征值是否满足Neimark-Sacker分岔的临界准则，这种数值搜寻的方法具有一定的盲目性和不确定性，非常耗时.另外，虽然可以采用极点配置方法找到满足特征值分布条件的系统参数点，但该方法也是先确定特征值后再确定参数，确定的参数对于系统仍存在机理不明确问题.特别是对于橫截条件，由于需要求特征值对分岔参数的导数，极点配置方法无法解决.因此为了克服传统分岔临界准则的局限性，本文采用不直接依赖于特征值计算的映射Neimark-Sacker分岔的显式临界条件来获得系统参数.

Pμ(λ)=λ4+a1λ3+a2λ2+a3λ+a4.

(12)

这里选取μ=(ρ,β)，ai=ai(ρ,β)是与分岔参数ρ和β有关的实数，i=1,…，4.针对建立的振动落砂机系统周期运动的Poincaré映射(11), 有如下的引理.

引理1[9]映射(11)在μ=μ0处发生Neimark-Sacker分岔，当且仅当满足下列条件(Ⅰ)～(Ⅴ)

(Ⅱ) 1+a1+a2+a3+a4>0,1+a4>0，

1-a1+a2-a3+a4>0,1-a4>0,

(Ⅴ)cos(2π/m)≠ψ,m=3,4,5，…,

ψ=1-

其中条件(Ⅰ)保证有一对复共轭特征值位于单位圆上；条件(Ⅱ)保证其它的特征值位于单位圆内；条件(Ⅲ)保证映射不动点是合理存在的；条件(Ⅳ)保证在参数扰动下，位于单位圆上的特征值穿越单位圆的速度不为零；条件(Ⅴ)保证Neimark-Sacker分岔是非共振的.

3.2振动落砂机系统Neimark-Sacker分岔的存在性

ρ

图2中白色区域Ⅰ和Ⅱ内的点都满足引理1的条件(Ⅱ)-(Ⅲ)中的不等式，但在灰色区域Ⅲ和Ⅳ中至少有一个条件(Ⅱ)和(Ⅲ)中的不等式不成立.由曲线EA，AB，BC，CD和DE围成的白色区域Ⅰ除了满足条件(Ⅱ)-(Ⅲ)还满足Δ>0，因此白色区域Ⅰ为系统周期运动的稳定参数区域.曲线DE是由临界条件(Ⅰ)中的Δ=0得到，曲线DE上由条件(Ⅴ)得到的点R3和R4分别为系统出现3阶和4阶强共振点.点划线l和m由横截条件(Ⅳ)不等式左边的表达式取等号得到的，这样两曲线l和m和曲线DE的交点T1和T2不满足Neimark-Sacker分岔的横截条件(Ⅳ).由此在选取系统参数临界点时应该避开这些强共振点和非横截点.在由曲线DE，EF，FG和GD围成的白色区域Ⅱ内，在曲线DE的附近是出现系统拟周期碰撞运动的潜在区域.为了分析分岔解的稳定性，在分岔图的白色区域内的Neimark-Sacker分岔临界曲线DE上任取一点μ0=(ρ0,β0)=(0.6,0.954 41)作为临界分岔值.

3.3振动落砂机系统拟周期碰撞运动的稳定性

振动落砂机系统出现的拟周期碰撞运动的稳定性，也即Neimark-Sacker分岔解(不变圈)的稳定性取决于映射(11)的非线性项.采用中心流形-范式方法[10]或者频域方法[11]都可以分析Neimark-Sacker分岔的稳定性.本文使用投影法[12]来分析分岔的稳定性.

取坐标变换

Yk=Xk-X*,ν=μ-μ0.

(13)

(13)式中的X*为映射(11)的不动点，μ0为临界分岔参数值.

映射(11)经过坐标变换(13)变换成

(14)

则通过变量变换后的新映射(14)的不动点和分岔点都转化为了零点.这样映射(11)发生Neimark-Sacker分岔后产生的拟周期碰撞运动的稳定性可由引理2来确定.

(15)

(15)式中的g21，g20，g11和g02可通过中心流形-正则化方法得到，它们具有下面的形式：

(16)

根据(15)和(16)，在μ=μ0计算得到

α(0)=-0.089<0.

因此，根据引理2可判断系统会产生一个稳定的Neimark-Sacker分岔解，即系统稳定的拟周期碰撞运动.

3.4数值实验

为了验证上述理论分析和研究振动落砂机系统在分岔点附近的动力学行为，其它三个参数不变，临界参数ρ0不变，变化分岔参数β，在落砂机的分岔临界点附近设置了6组参数扰动值，并做了相应的数值仿真，文中的数值仿真都采用4 000次碰撞.在Poincaré映射分岔图的白色区域Ⅰ内的分岔临界曲线DE附近取一组分岔参数μ=μ0+Δμ=(ρ0,β0-0.02)，其中Δμ是临界参数扰动量，设置映射初始值X=X*+ΔX，其中ΔX=(0,0,0,0.001)T为不动点扰动量，在该参数点处系统处于稳定的周期运动，即Poincaré映射上一个不动点，如图3(a)所示.在Poincaré映射分岔图的白色区域Ⅱ内并且充分接近曲线DE的参数区域内取分岔参数μ= (ρ0,β0+0.004)，在该参数点处系统处于稳定的拟周期碰撞运动，即Poincaré映射上一个不变圈，如图3(b)所示.继续变化分岔参数值，当取分岔参数μ= (ρ0,β0+0.045 6)时，拟周期运动失稳，产生锁相运动，如图3(c)所示.当取参数μ= (ρ0,β0+0.046)时系统退出锁相运动，又产生拟周期吸引不变圈，如图3(d)所示.当取控制参数μ= (ρ0,β0+0.1)时系统产生拟周期吸引不变圈，不变圈幅值更大，如图3(e)所示.当继续扰动参数值至μ= (ρ0,β0+0.22)时系统经锁相转迁为混沌运动，如图3(f)所示.通过以上仿真分析，振动落砂机系统展示出了丰富的动力学行为，在Hopf分岔临界点附近作参数扰动，当参数扰动量足够小的时候由Hopf分岔所产生的不变圈在形状上类似于一个椭圆.随着参数扰动量的增大，不变圈不断增大，其形状也变得越来越不规则.仿真显示只有当系统的扰动参数很大时才发生混沌，说明系统的拟周期碰撞运动具有强的鲁棒性和较大的稳定域.因此，在设计系统参数时，在白色区域Ⅱ邻域内靠近分岔临界曲线DE处取系统参数可以产生稳定的拟周期碰撞运动.

xk

xk

xk

xk

xk

xk

4结论

1)基于主动利用Neimark-Sacker分岔解特性的思想，通过选定合适的系统参数，设计出了稳定的拟周期碰撞的振动落砂机系统.

2)利用显式的Neimark-Sacker分岔临界准则获得了落砂机系统产生拟周期碰撞运动的两参数区域图，此参数区域具有较大的分岔可行域范围，可保障所产生的拟周期碰撞运动具有较大的稳定域和较强的鲁棒性.

3)数值分析实现了落砂机系统产生的稳定的拟周期碰撞运动并调查了附近锁相和混沌等动力学行为.

参考文献

[1]WAGG D J. Periodic sticking motion in a two-degree-of-freedom impact oscillator[J]. International Journal of Non-Linear Mechanics, 2005, 40(8): 1076-1087.

[2]WEN G L, XIE J H. Period-doubling bifurcation and non-typical routes to chaos of a two-degree-of-freedom vibro-impact system[J]. Journal of Applied Mechanics, 2001, 68(4): 670-674.

[3]LUO A C J, CHEN L D. Periodic motions and grazing in a harmonically forced piecewise linear oscillator with impacts[J]. Chaos Solitons & Fractals, 2005, 24(2):567-578.

[4]AGUIAR R R, WEBER H I. Mathematical modeling and experimental investigation of an embedded vibro-impact system[J]. Nonlinear Dynamics, 2011, 65(3): 317-334.

[5]LUO G W, XIE J H. Bifurcations and chaos in a system with impacts[J]. Physica D, 2001, 148(3/4):183-200.

[6]张思进，尹磊磊，文桂林. 一类拟Hamilton 碰振系统的全局分岔及多解共存现象分析[J]. 湖南大学学报：自然科学版，2014，41(10): 55-61.

ZHANG Si-jin, YIN Lei-lei, WEN Gui-lin.Analysis of the global bifurcations and coexistence of multiple solutions for a general quasi-Hamiltonian vibro-impact system [J]. Journal of Hunan University: Natural Sciences,2014, 41(10): 55-61. (In Chinese)

[7]罗冠炜，谢建华. 冲击振动落砂机的周期运动与分岔[J]. 机械工程学报，2003，39(1): 74-78.

LUO Guan-wei, XIE Jian-hua. Stability and bifurcation of periodic motions of the inertial shakers[J]. Journal of Mechanical Engineering, 2003, 39(1): 74-78. (In Chinese)

[8]丁旺才，谢建华，李万祥. 碰撞振动系统强共振下的两参数动力学分析[J]. 计算力学学报，2004，21(6): 658-664.

DING Wang-cai, XIE Jian-hua, LI Wan-xiang. Dynamics analysis of two-parameter family in a case of strong resonance of vibro-impact system[J]. Chinese Journal of Computational Mechanics, 2004, 21(6): 658-664. (In Chinese)

[9]WEN G L. Criterion to identify Hopf bifurcations in maps of arbitrary dimension[J].Physical Review E, 2005, 72(2): 026201-026204.

[10]罗冠炜, 谢建华. 碰撞振动系统的周期运动和分岔[M]. 北京: 科学出版社, 2004 : 187-188.

LUO Guan-wei, XIE Jian-hua. Periodic motion and bifurcation of vibro-impact system[M]. Beijing: Science Press, 2004: 187-188. (In Chinese)

[11]D'AMICO M B, MOIOLA J L, PAOLINI E E. Controlling bifurcations in maps via a frequency domain approach[J]. Dynamics of Continuous Discrete and Impulsive Systems Series B-applications & Algorithms, 2003, 10(6): 781-798.

[12]KUZNETSOV Y A. Elements of applied bifurcation theory [M]. 2nd ed. New York: Springer-Verlag, 1998: 185-187.

[13]IOOSS G. Bifurcation of maps and applications [M]. Amsterdam: North-Holland,1979:28-31.

Design of Quasi-periodic Impact Motion of an Impact Shaker System

WU Xin1，2, WEN Gui-lin1†, HE Li-ping1, XU Hui-dong1, WEI Ke-xiang2

(1.State Key Laboratory of Advanced Design and Manufacturing for Vehicle Body, Hunan Univ, Changsha, Hunan410082,China；2.School of Mechanical Engineering, Hunan Institute of Engineering , Xiangtan, Hunan411104,China)

Abstract:The Poincaré map of periodic motion was established, and the bifurcation diagram of two parameters was obtained by using an explicit critical criterion without using eigen values when the difficulties caused by the classical Neimark-Sacker bifurcation critical criterion described by the properties of eigenvalues were considered. The bifurcation diagram helps to proactively design the quasi-periodic impact motion of the system. Then, the stability of the quasi-periodic impact motion was further analyzed by utilizing the center manifold and normal formal theory. Finally, numerical experiments verify that the stable quasi-periodic impact motion can be generated at chosen parameters points.

Key words:shaker system; impact vibration; Neimark-Sacker bifurcation; quasi-periodic impact motion; stability

中图分类号：O322；TB123

文献标识码：A

作者简介：伍新(1976-)，男，湖南新化人，湖南大学博士研究生，湖南工程学院讲师†通讯联系人，E-mail:glwen@hnu.edu.cn

基金项目：国家杰出青年科学基金资助项目(11225212)；国家自然科学基金资助项目(11172100,11002052),National Natural Science Foundation of China(11172100,11002052) ；湖南省自然科学基金资助项目(2016JJ4027)

收稿日期：2015-01-30

文章编号：1674-2974(2016)04-0038-06