一类含间隙振动系统的周期运动稳定性、分岔与混沌形成过程研究

建立了两自由度含间隙振动系统对称周期碰撞运动的Poincar啨映射方程 ,讨论了该映射不动点的稳定性与局部分岔 .通过数值仿真研究了含间隙振动系统对称周期碰撞运动经叉式分岔、倍化分岔、“擦边”奇异性向混沌转迁的全局分岔过程     

碰撞振动系统的一类余维二分岔及T2环面分岔

建立了三自由度碰撞振动系统的动力学模型及其周期运动的Poincaré映射,当Jacobi矩阵存在两对共轭复特征值同时在单位圆上时,通过中心流形-范式方法将六维映射转变为四维范式映射.理论分析了这种余维二分岔问题,给出了局部动力学行为的两参数开折.证明系统在一定的参数组合下,存在稳定的Hopf分岔和T2环面分岔.数值计算验证了理论结果.     

齿轮传动系统共存吸引子的不连续分岔

大量的多吸引子共存是引起齿轮传动系统具有丰富动力学行为的一个重要因素.多吸引子共存时,运动工况的变化以及不可避免的扰动都可能导致齿轮传动系统在不同运动行为之间跳跃变换,对整个机器产生不良的影响.目前,一些隐藏的吸引子没有被发现,共存吸引子的分岔演化规律没有被完全揭示.考虑单自由度直齿圆柱齿轮传动系统,构建由局部映射复合的Poincaré映射,给出Jacobi矩阵特征值计算的半解析法.应用数值仿真、延拓打靶法和Floquet特征乘子求解共存吸引子的稳定性与分岔,应用胞映射法计算共存吸引子的吸引域,讨论啮合频率、阻尼比和时变激励幅值对系统动力学的影响,揭示齿轮传动系统倍周期型擦边分岔、亚临界倍周期分岔诱导的鞍结分岔和边界激变等不连续分岔行为.倍周期分岔诱导的鞍结分岔引起相邻周期吸引子相互转迁的跳跃与迟滞,使倍周期分岔呈现亚临界特性.鞍结分岔是共存周期吸引子出现或消失的主要原因.边界激变引起混沌吸引子及其吸引域突然消失,对应周期吸引子的分岔终止.     

一类两自由度干摩擦系统倍化分岔与混沌的数值模拟

本文利用Runge-Kutta数值积分法对一类两自由度干摩擦系统进行了数值仿真.文中列出了该系统中质块与传送带间处于黏附与滑移两种不同状态时质块的运动方程,同时给出了两种状态间的转换条件,在此基础上对系统进行数值模拟.通过不断变化传送带的速度对该系统的动力行为进行了分析,证实此系统存在稳定的周期运动及倍化分岔,得到了该系统周期运动的倍化分岔通向混沌的路径.对此系统分岔与混沌行为的研究为生产实践中的干摩擦现象的优化配置提供了理论依据.     

外激励作用下亚音速二维壁板分岔及响应研究

研究了亚音速流中二维壁板在外激励作用下的分岔和响应问题. 采用Galerkin方法将非线性运动控制方程离散为常微分方程组. 采用Runge-Kutta数值方法进行了数值计算,研究了壁板系统非单周期区在参数空间的分布情况. 结果表明: 在参数空间中, 非单周期区和单周期区会交替出现; 在不同的单周期区内, 系统运动轨线也在有规律的变化; 系统由单周期运动进入混沌运动是经过一系列周期倍化分岔产生的.      

一类碰撞振动系统在内伊马克沙克-音叉分岔点附近的局部两参数动力学

考虑一类具有对称性的三自由度碰撞振动系统.系统的庞加莱映射在一定条件下存在对称不动点,对应于系统的对称周期运动.根据对称性导出庞加莱映射P是另外一个隐式虚拟映射Q的二次迭代.推导了庞加莱映射对称不动点的解析表达式.根据映射不动点的稳定性及分岔理论,映射P的对称不动点发生内伊马克沙克-音叉(Neimark--Saker-pitchfork)分岔对应于映射Q发生内伊马克沙克-倍化(Neimark--Sakerflip)分岔.利用隐式虚拟映射Q,通过对范式作两参数开折分析,研究了映射P的对称不动点在内伊马克沙克-音叉分岔点附近的局部动力学行为.碰撞振动系统在这个余维二分岔点附近的局部动力学行为可能表现为投影后的庞加莱截面上的单一对称不动点、一对共轭不动点、单一对称拟周期吸引子以及一对共轭拟周期吸引子.数值模拟得到了内伊马克沙克-音叉分岔点附近的各种可能情况.内伊马克沙克-分岔和音叉分岔互相作用可能产生新的结果:对称不动点虽然首先分岔为两个共轭不动点,但是这两个共轭不动点是不稳定的,最终收敛到同一个对称拟周期吸引子.     

多自由度参激系统稳定性的数值分析

提出多自由度周期参激系统稳定性的数值直接法。通过将扰动方程表示成状态方程形式,再根据Flo-quet理论将扰动解表示成指数特征分量与周期分量之积,并将其周期分量与系统周期系数展成Fourier级数,导出一系列代数方程,建立矩阵特征值问题,从而由数值求解特征值可直接确定参激系统的稳定性。该方法可用于一般周期参激阻尼系统,特征值矩阵不含逆子阵。应用于斜拉索在支座周期运动激励下的参激振动不稳定性分析,数值结果表明该方法的有效性。     

Duffing系统在双参数平面上的动力学特性分析

将单参数最大Lyapunov指数的计算推广到双参数平面上,数值计算Duffing系统在双参数平面上的最大Lyapunov指数,得到系统在参数平面上周期运动、混沌运动、各种分岔曲线的参数区域;结合系统单参数分岔图、相图、庞加莱截面图讨论了系统在参数平面上的分岔混沌过程以及阻尼系数对系统双参数特性的影响。结果表明:在双参数平面上系统出现了周期跳跃、周期倍化分岔、叉式分岔等复杂的分岔曲线,而且这些分岔曲线随阻尼系数的增加不断发生着复杂变化;得到系统在以往单参数分岔过程中很少出现的分岔曲线相交、嵌套、演变等特殊现象;阻尼系数对系统双参数耦合动力学特性有重要的影响。本文对工程中其它多参数系统的参数耦合特性的研究具有一定的参考价值。     

雨刮器非线性摩擦振动的稳定性与分岔特性

针对雨刮器建立2自由度非线性摩擦振动动力学模型,基于复模态理论计算复特征值并进行稳定性及其对刮刷速度的依赖性分析;通过数值计算分析摩擦振动对刮刷速度的分岔特性,并利用相轨迹、庞加莱映射、频谱特性分析不同刮刷速度下的非线性振动现象.研究发现:摩擦-速度特性的负斜率是导致系统不稳定的根本原因,增大刮刷速度有利于提高系统的稳定性;在高、低刮速区,随着刮刷速度的下降,系统振动形态遵循周期→准周期→混沌的演化规律,并会伴随显著的粘滑振动;仅高速区的周期振动和非振动条件下,刮刷时无附加的粘滑振动.     

两自由度机械臂λ^2=1下的Hopf分岔与混沌研究

研究了一类周期系数力学系统因周期运动失稳而产生Hopf分岔及混沌问题.首先根据拉格朗日方程给出了该力学系统的运动微分方程,并确定其周期运动的具有周期系数的扰动运动微分方程,再根据Floquet理论建立了其给定周期运动的Poincaré映射,根据该系统的特征矩阵有一对复共轭特征值从-1处穿越单位圆情况,分析该Poincaré映射不动点失稳后将发生次谐分岔、Hopf分岔、倍周期分岔,而多次倍周期分岔将导致混沌.并用数值计算加以验证.结果表明,随着分岔参数的变化,系统的周期运动可通过次谐分岔形成周期2运动,进而发生Hopf分岔形成拟周期运动,并再次经次谐分岔、倍周期分岔形成混沌运动.     

RESEARCH IN STABILITY OF PERIODIC MOTION,BIFURCATION AND CHAOS IN A VIBRATORY SYSTEM WITH A CLEARANCE

A two-degrees-of-freedom vibratory system with a clearance or gap is under consideration based on the Poincard map. Stability and local bifurcation of the period-one doubleimpact symmetrical motion of the system are analyzed by using the equation of map. The routes from periodic impact motions to chaos, via pitchfork bifurcation, period-doubling bifurcation and grazing bifurcation, are studied by numerical simulation. Under suitable system parameter conditions, Neimark-Sacker bifurcations associated with periodic impact motion can occur in the two-degrees-of-freedom vibro-impact system.     

时变张力作用下轴向运动梁的分岔与混沌

轴向运动结构的横向参激振动一直是非线性动力学领域的研究热点之一. 目前研究较多的是轴向速度摄动的动力学模型,参数激励由速度的简谐波动产生. 但在工程应用中,存在轴向张力波动的运动结构较为广泛,而针对轴向张力摄动的模型研究较少. 本文研究了时变张力作用下轴向变速运动黏弹性梁的分岔与混沌. 考虑随着时间周期性变化的轴向张力,计入线性黏性阻尼,采用Kelvin模型的黏弹性本构关系,给出了梁横向非线性 振动的积分--偏微分控制方程. 首先应用四阶Galerkin截断方法将控制方程离散化,然后采用四阶Runge-Kutta方法计算系统的数值解,进而确定其动力学行为. 基于梁中点的横向位移和速度的数值结果,仿真了梁沿平均轴速、张力摄动幅值、张力摄动频率以及黏弹性系数变化的倍周期分岔与混 沌运动,并且通过计算系统的最大李雅普诺夫指数来识别其混沌行为. 结果表明:较小的平均轴速有助于梁的周期运动,梁在临界速度附近容易发生倍周期分岔与混沌行为. 随着张力摄动幅值的增大,梁的振动幅值的混沌区间不断增大. 较小的黏弹性系数和张力摄动频率更容易使梁发生混沌运动. 最后,给出时程图、频谱图、相图以及Poincaré 映射图来确定梁的混沌运动.      

Symmetry, cusp bifurcation and chaos of an impact oscillator between two rigid sides

Both the symmetric period n-2 motion and asymmetric one of a one-degree- of-freedom impact oscillator are considered.The theory of bifurcations of the fixed point is applied to such model,and it is proved that the symmetric periodic motion has only pitchfork bifurcation by the analysis of the symmetry of the Poincarémap.The numerical simulation shows that one symmetric periodic orbit could bifurcate into two antisymmet- ric ones via pitchfork bifurcation.While the control parameter changes continuously, the two antisymmetric periodic orbits will give birth to two synchronous antisymmetric period-doubling sequences,and bring about two antisymmetric chaotic attractors subse- quently.If the symmetric system is transformed into asymmetric one,bifurcations of the asymmetric period n-2 motion can be described by a two-parameter unfolding of cusp, and the pitchfork changes into one unbifurcated branch and one fold branch.     

Bifurcation and chaos of a simple walking model driven by a rhythmic signal

The walk of animals is achieved by the interaction between the dynamics of their mechanical system and the central pattern generator (CPG). In this paper, we analyze dynamic properties of a simple walking model of a biped robot driven by a rhythmic signal from an oscillator. In particular, we examine the long-term global behavior and the bifurcation of the motion that leads to chaotic motion, depending on the model parameter values. The simple model consists of a hip and two legs connected at the hip through a rotational joint. The joint is driven by a rhythmic signal from an oscillator, which is an open loop. In order to analyze the bifurcation, we first obtained approximate solutions of the walking motion and then constructed discrete dynamics using the Poincaré map. As a result, we found that consecutive period-doubling bifurcations occur as the model parameter values change, and that the walking motion leads to chaotic motion over the critical value of the model parameters. Moreover, we approximately obtained the period-doubling solutions and the critical value by employing a Newton-Raphson method. Our analytical results were verified by the numerical simulations.     

计及短周期误差的直齿轮副近周期运动及其辨识

齿轮副中的齿距偏差等短周期误差使系统出现复杂的周期运动, 影响齿轮传动的平稳性. 将该类复杂周期运动定义为近周期运动, 采用多时间尺度Poincaré映射截面对其进行辨识. 为研究齿轮副的近周期运动, 引入含齿距偏差的直齿轮副非线性动力学模型, 并计入齿侧间隙与时变重合度等参数. 采用变步长4阶Runge-Kutta法数值求解动力学方程, 由所提出的辨识方法分析不同参数影响下系统的近周期运动. 根据改进胞映射法计算系统的吸引域, 结合多初值分岔图、吸引域图与分岔树状图等研究了系统随扭矩与啮合频率变化的多稳态近周期运动. 研究结果表明, 齿轮副中的短周期误差导致系统的周期运动变复杂, 在微观时间尺度内, 系统的Poincaré映射点数呈现为点簇形式, 系统的点簇数与实际运动周期数为宏观时间尺度的Poincaré映射点数. 短周期误差导致系统在微观时间尺度内的吸引子数量增多, 使系统运动转迁过程变复杂. 合理的参数范围及初值范围可提高齿轮传动的平稳性. 该辨识与分析方法可为非线性系统中的近周期运动研究奠定理论基础.      

Stability and bifurcation for a flexible beam under a large linear motion with a combination parametric resonance

Stability and bifurcation behaviors for a model of a flexible beam undergoing a large linear motion with a combination parametric resonance are studied by means of a combination of analytical and numerical methods. Three types of critical points for the bifurcation equations near the combination resonance in the presence of internal resonance are considered, which are characterized by a double zero and two negative eigenvalues, a double zero and a pair of purely imaginary eigenvalues, and two pairs of purely imaginary eigenvalues in nonresonant case, respectively. The stability regions of the initial equilibrium solution and the critical bifurcation curves are obtained in terms of the system parameters. Especially, for the third case, the explicit expressions of the critical bifurcation curves leading to incipient and secondary bifurcations are obtained with the aid of normal form theory. Bifurcations leading to Hopf bifurcations and 2-D tori and their stability conditions are also investigated. Some new dynamical behaviors are presented for this system. A time integration scheme is used to find the numerical solutions for these bifurcation cases, and numerical results agree with the analytic ones.     

Complex dynamics of a harmonically excited structure coupled with a nonlinear energy sink

Nonlinear behaviors are investigated for a structure coupled with a nonlinear energy sink. The structure is linear and subject to a harmonic excitation, modeled as a forced single-degree-of-freedom oscillator. The nonlinear energy sink is modeled as an oscillator consisting of a mass,a nonlinear spring, and a linear damper. Based on the numerical solutions, global bifurcation diagrams are presented to reveal the coexistence of periodic and chaotic motions for varying nonlinear energy sink mass and stiffness. Chaos is numerically identified via phase trajectories, power spectra,and Poincaré maps. Amplitude-frequency response curves are predicted by the method of harmonic balance for periodic steady-state responses. Their stabilities are analyzed.The Hopf bifurcation and the saddle-node bifurcation are determined. The investigation demonstrates that a nonlinear energy sink may create dynamic complexity.     

Delay equations with fluctuating delay related to the regenerative chatter

The mathematical models representing machine tool chatter dynamics have been cast as differential equations with delay. In this paper, non-linear delay differential equations with periodic delays which model the machine tool chatter with continuously modulated spindle speed are studied. The explicit time-dependent delay terms, due to spindle speed modulation, are replaced by state-dependent delay terms by augmenting the original equations. The augmented system of equations is autonomous and has two pairs of pure imaginary eigenvalues without resonance. The reduced bifurcation equation is obtained by making use of Lyapunov-Schmidt Reduction method. By using the reduced bifurcation equations, the periodic solutions are determined to analyze the tool motion. Analytical results show both modest increase of stability and existence of periodic solutions near the new stability boundary.     

Stability and bifurcation for limit cycle oscillations of an airfoil with external store

In this paper, stability and local bifurcation behaviors for the nonlinear aeroelastic model of an airfoil with external store are investigated using both analytical and numerical methods. Three kinds of degenerated equilibrium points of bifurcation response equations are considered. They are characterized as (1) one pair of purely imaginary eigenvalues and two pairs of conjugate complex roots with negative real parts; (2) two pairs of purely imaginary eigenvalues in nonresonant case and one pair of conjugate complex roots with negative real parts; (3) three pairs of purely imaginary eigenvalues in nonresonant case. With the aid of Maple software and normal form theory, the stability regions of the initial equilibrium point and the explicit expressions of the critical bifurcation curves are obtained, which can lead to static bifurcation and Hopf bifurcation. Under certain conditions, 2-D tori motion may occur. The complex dynamical motions are considered in this paper. Finally, the numerical solutions achieved by the fourth-order Runge–Kutta method agree with the analytic results.