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

现有参激系统的动力稳定性问题研究主要集中在主不稳定区域上。为获得组合不稳定区域,基于Floquet方法,采用Bolotin方法在不同周期数下设解形式,结合特征值分析法得到确定多自由度参激系统动力不稳定区域的数值解法。对一个两自由度受周期轴向力的旋转轴系算例的稳定性分析,发现通过增加设解近似项数可获得高阶不稳定区域,且各阶不稳定区域边界随近似次数的增加逐渐趋于稳定,此外,增大阻尼可使各不稳定区域边界变得更加平滑。本文方法可用于一般多自由度周期参激阻尼系统,是一种简明易操作的直接数值解法。     

基于可观测状态的轴承-转子系统周期解计算及稳定性分析

分析了轴承-转子系统的稳定性和分岔,基于系统可观测状态信息给出1种求解系统周期解及识别周期解稳定性的方法,同时将该方法与Floquet理论相结合分析系统周期解的稳定性及失稳分岔形式,将转速作为分岔参数分析系统响应的周期、拟周期、多解共存和跳跃现象.结果表明,采用该方法计算系统周期解及稳定性时,利用系统可观测稳态和瞬态信息,即可求解出系统Jacobian矩阵而无需实时求解轴承非线性油膜力的Jacobian矩阵.与传统PNF方法相比,该方法不仅具有很高的精度而且可以节约计算量,同时可以预测追踪随控制参数变化的系统周期解及其稳定性,可用于指导轴承-转子系统的非线性动力学设计.     

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

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

一类双自由度碰振系统运动分析

基于Poincare映射方法对一类两自由度碰撞系统进行了分析。经过详细的理论演算得到单碰周期n的次谐运动的存在性判据和稳定性条件,给出计算Jacobi矩阵特征值的公式。数值模拟表明,该方法具有令人满意的结果。此外,还讨论了当不满足所提出的单碰周期n次谐运动的存在性条件时,可能会出现的运动形式。     

一类三维随机动力系统的最大Lyapunov指数

基于一维扩散过程的奇异边界理论,使用摄动方法研究了白噪声参激的一类余维二分岔系统的最大Lyapunov指数渐近表达式和数值解,主要讨论了一维相扩散过程同时存在两类奇异边界以及FPK方程存在平稳解的一般性条件. 通过对参激噪声作用项系数矩阵的分析,给出了不变测度的解析解及其相应的Monte Carlo数值仿真结果,并导出了一维相扩散过程P分岔点的确定方法. 对于一类特殊情形,给出了最大Lyapunov指数的渐近表达式;对于参激噪声作用项系数矩阵的一般情形,则给出了系统最大Lyapunov指数的数值结果.      

基于Fourier级数的时变周期系数Riccati微分方程精细积分

结合Fourier级数展开方法,本文提出了基于精细积分的时变周期系数Riccati微分方程求解高效算法.首先,利用Fourier级数展开方法将周期系统表示成三角级数形式,在一个积分步内使用精细积分方法得到对应Hamilton系统状态转移矩阵的表达式.然后,通过Riccati变换的方法,得到含有状态转移矩阵的时变周期系数Riccati微分方程解的递推格式.本文方法充分利用了方程本身的周期性特点,文中的数值算例表明算法具有计算效率高、结果可靠等优势.     

粘性可压混合层时间稳定性对称紧致差分求解

基于可压扰动方程组的一阶改型 ,将高精度对称紧致格式引入边值法数值线性稳定性分析。对所获非线性离散特征值问题给出了一个通用形式二阶迭代局部算法 ,实现了时间模式和空间模式的统一求解 ,并将扰动特征值及其特征函数同时得到。据此分析了可压平面自由混合层时间稳定性 ,涉及二维 /三维扰动波、粘性 /无粘扰动波、第一 /第二模态、特征函数、伪特征值谱等。研究表明 ,压缩性效应和粘性效应对最不稳定扰动波数和增长率呈相似的减抑作用 ;在 Mc=1附近 ,从高波数段开始 ,粘性效应可强化二维不稳定扰动波由第一模态向第二模态的过渡     

惯性式冲击振动落砂机周期倍化分岔的反控制

在不改变惯性式冲击振动落砂机系统平衡解结构的前提下,考虑碰撞振动系统的Poincaré映射的隐式特点以及经典的映射周期倍化分岔临界准则给反控制带来的困难,基于不直接依赖于特征值计算的周期倍化分岔显式临界准则,研究了落砂机系统周期倍化分岔的反控制。论文首先对落砂机系统施加线性反馈控制,得到受控闭环系统的Poincare映射,并应用不直接依赖于特征值计算的周期倍化分岔显式临界准则,获得了系统发生周期倍化分岔的控制参数区域。然后应用中心流形-正则形方法分析了周期倍化分岔的稳定性。最终采用数值仿真验证了在任意指定的系统参数点通过控制能产生稳定的周期倍化分岔解。     

深水钻井隔水管在平台升沉运动下的参激振动响应分析

建立了深水钻井隔水管在平台升沉运动作用下的参激振动控制方程;利用振型叠加法将控制方程转化成马蒂厄方程形式,并用摄动法对马蒂厄方程进行求解;考虑海水的阻尼作用,得到了隔水管参激振动的马蒂厄不稳定区。最后对参激振动方程进行了实例求解,得到了参激振动下的隔水管位移,并在此基础上对参激振动的稳定性及隔水管强度进行了分析。结果表明:在给定的平台升沉周期与升沉振幅条件下,隔水管外径、壁厚、顶张力的变化均有可能引发隔水管产生参激共振;在不产生参激共振的情况下,隔水管Mises应力随隔水管外径与壁厚的增大而减小,随顶张力与升沉补偿装置轴向刚度系数的增大而增大;最大Mises应力均出现在隔水管顶部,沿水深近似成线性分布。     

广义特征值问题的并行块Jacobi-Davidson方法及应用

给出了对称矩阵广义特征值问题AX=λBX的并行块Jacobi-Davidson方法.该方法使用投影技术将大型矩阵特征值问题转变成低维子空间中矩阵特征值问题,并利用Neumann级数展开对校正方程进行预处理.该方法可同时并行计算广义特征值问题的几个极端特征对,具有良好的并行性.将这一方法应用于某型号机翼及挂架的结构动力分析并行计算,在IBM-P650并行计算机上的数值试验结果表明,在相同迭代精确度的条件下,Jacobi-Davidson方法比子空间迭代法使用较少的迭代次数和运算时间,并具有更高的加速比和并行效率.     

A new numerical technique for the analysis of parametrically excited nonlinear systems

A new computational scheme using Chebyshev polynomials is proposed for the numerical solution of parametrically excited nonlinear systems. The state vector and the periodic coefficients are expanded in Chebyshev polynomials and an integral equation suitable for a Picard-type iteration is formulated. A Chebyshev collocation is applied to the integral with the nonlinearities reducing the problem to the solution of a set of linear algebraic equations in each iteration. The method is equally applicable for nonlinear systems which are represented in state-space form or by a set of second-order differential equations. The proposed technique is found to duplicate the periodic, multi-periodic and chaotic solutions of a parametrically excited system obtained previously using the conventional numerical integration schemes with comparable CPU times. The technique does not require the inversion of the mass matrix in the case of multi degree-of-freedom systems. The present method is also shown to offer significant computational conveniences over the conventional numerical integration routines when used in a scheme for the direct determination of periodic solutions. Of course, the technique is also applicable to non-parametrically excited nonlinear systems as well.     

Moment Lyapunov exponent for three-dimensional system under real noise excitation

The pth moment Lyapunov exponent of a two-codimension bifurcation system excited parametrically by a real noise is investigated. By a linear stochastic transformation, the differential operator of the system is obtained. In order to evaluate the asymptotic expansion of the moment Lyapunov exponent, via a perturbation method, a ralevant eigenvalue problem is obtained. The eigenvalue problem is then solved by a Fourier cosine series expansion, and an infinite matrix is thus obtained, whose leading eigenvalue is the second-order of the asymptotic expansion of the moment Lyapunov exponent. Finally, the convergence of procedure is numerically illustrated, and the effects of the system and the noise parameters on the moment Lyapunov exponent are discussed.     

Parametric Identification of a Base-Excited Single Pendulum

A harmonic balance based identification algorithm was applied to the simulated single pendulum with horizontal base-excitation. The purpose of this simulation was to examine the applicability of the algorithm on parametrically excited, whirling chaotic systems. Modifications were adopted to adapt to the whirling systems. The system was supposed to be unknown except only the excitation frequency. Linear interpolation functions and the Fourier series functions were tested to approximate unknown nonlinear functions in the governing differential equation. After extracting unstable periodic orbits, all of the parameters were simultaneously identified. By direct comparison, Poincaré section plots and reconstructed phase portrait techniques, it was shown that the identified system had similar dynamical characteristics to the original simulated pendulum, which implies the effectiveness of the examined algorithm.     

Perturbation analysis in parametrically excited two-degree-of-freedom systems with quadratic and cubic nonlinearities

Based on temporal rescaling and harmonic balance, an extended asymptotic perturbation method for parametrically excited two-degree-of-freedom systems with square and cubic nonlinearities is proposed to study the nonlinear dynamics under 1:2 internal resonance. This asymptotic perturbation method is employed to transform the two-degree-of-freedom nonlinear systems into a four-dimensional nonlinear averaged equation governing the amplitudes and phases of the approximation solutions. Linear stable analysis at equilibrium solutions of the averaged equation is done to show bifurcations of periodic motion and homoclinic motions. Furthermore, analytical expressions of homoclinic orbits and heteroclinic cycles for the averaged equation without dampings are obtained. Considering the action of the damping, the bifurcations of limit cycles are also investigated. A concrete example is further provided to discuss the correctness and accuracy of the extended asymptotic perturbation method in the case of small-amplitude motion for the two-degree-of-freedom nonlinear system.     

Non-smooth stability analysis of the parametrically excited impact oscillator

The aim of this paper is to give a Lyapunov stability analysis of a parametrically excited impact oscillator, i.e. a vertically driven pendulum which can collide with a support. The impact oscillator with parametric excitation is described by Hill's equation with a unilateral constraint. The unilaterally constrained Hill's equation is an archetype of a parametrically excited non-smooth dynamical system with state jumps. The exact stability criteria of the unilaterally constrained Hill's equation are rigorously derived using Lyapunov techniques and are expressed in the properties of the fundamental solutions of the unconstrained Hill's equation. Furthermore, an asymptotic approximation method for the critical restitution coefficient is presented based on Hill's infinite determinant and this approximation can be made arbitrarily accurate. A comparison of numerical and theoretical results is presented for the unilaterally constrained Mathieu equation.     

Bifurcation Analyses in the Parametrically Excited Vibrations of Cylindrical Panels

We consider parametrically excited vibrations of shallow cylindrical panels. The governing system of two coupled nonlinear partial differential equations is discretized by using the Bubnov–Galerkin method. The computations are simplified significantly by the application of computer algebra, and as a result low dimensional models of shell vibrations are readily obtained. After applying numerical continuation techniques and ideas from dynamical systems theory, complete bifurcation diagrams are constructed. Our principal aim is to investigate the interaction between different modes of shell vibrations under parametric excitation. Results for system models with four of the lowest modes are reported. We essentially investigate periodic solutions, their stability and bifurcations within the range of excitation frequency that corresponds to the parametric resonances at the lowest mode of vibration.     

Bifurcation Analysis of Parametrically Excited Rayleigh–Liénard Oscillators

A parametrically excited Rayleigh–Liénard oscillator is investigatedby an asymptotic perturbation method based on Fourier expansion and timerescaling. Two coupled equations for the amplitude and the phase ofsolutions are derived and the stability of steady-state periodic solutionsas well as parametric excitation-response and frequency-response curvesare determined. Comparison with the parametrically excited Liénardoscillator is performed and analytic approximate solutions are checkedusing numerical integration. Dulac's criterion, thePoincaré–Bendixson theorem, and energy considerations are used in order to study the existence and characteristics of limit cycles of the twocoupled equations. A limit cycle corresponds to a modulated motion forthe Rayleigh–Liénard oscillator. Modulated motion can be also obtainedfor very low values of the parametric excitation, and in this case, anapproximate analytic solution is easily constructed. If the parametricexcitation is increased, an infinite-period bifurcation is observed because the modulation period lengthens and becomes infinite, while themodulation amplitude remains finite and suddenly the attractor settlesdown into a periodic motion. Floquet's theory is used to evaluatethe stability of the periodic solutions, and in certain cases,symmetry-breaking bifurcations are predicted. Numerical simulationsconfirm this scenario and detect chaos and unbounded motions in theinstability regions of the periodic solutions.     

Nonlinear dynamics of idler gear systems

This work examines the nonlinear, parametrically excited dynamics of idler gearsets. The two gear tooth meshes provide two interacting parametric excitation sources and two possible tooth separations. The periodic steady state solutions are obtained using analytical and numerical approaches. Asymptotic perturbation analysis gives the solution branches and their stabilities near primary, secondary, and subharmonic resonances. The ratio of mesh stiffness variation to its mean value is the small parameter. The time of tooth separation is assumed to be a small fraction of the mesh period. With these stipulations, the nonsmooth separation function that determines contact loss and the variable mesh stiffness are reformulated into a form suitable for perturbation. Perturbation yields closed-form expressions that expose the impact of key parameters on the nonlinear response. The asymptotic analysis for this strongly nonlinear system compares well to separate harmonic balance/arclength continuation and numerical integration solutions. The expressions in terms of fundamental design quantities have natural practical application.     

Arbitrary amplitude periodic solutions for parametrically excited systems with time delay

Periodic solutions for parametrically excited system under state feedback control with a time delay are investigated. Using the asymptotic perturbation method, two slow-flow equations for the amplitude and phase of the parametric resonance response are derived. Their fixed points correspond to limit cycles (phase-locked periodic solutions) for the starting system. In the system without control, periodic solutions (if any) exist only for fixed values of amplitude and phase and depend on the system parameters and excitation amplitude. In many cases, the amplitudes of periodic solutions do not correspond to the technical requirements. On the contrary, it is demonstrated that, if the vibration control terms are added, stable periodic solutions with arbitrarily chosen amplitude and phase can be accomplished. Therefore, an effective vibration control is possible if appropriate time delay and feedback gains are chosen.     

Stability and torsional vibration analysis of a micro-shaft subjected to an electrostatic parametric excitation using variational iteration method

In this article stability and parametrically excited oscillations of a two stage micro-shaft located in a Newtonian fluid with arrayed electrostatic actuation system is investigated. The static stability of the system is studied and the fixed points of the micro-shaft are determined and the global stability of the fixed points is studied by plotting the micro-shaft phase diagrams for different initial conditions. Subsequently the governing equation of motion is linearized about static equilibrium situation using calculus of variation theory and discretized using the Galerkin’s method. Then the system is modeled as a single-degree-of-freedom model and a Mathieu type equation is obtained. The Variational Iteration Method (VIM) is used as an asymptotic analytical method to obtain approximate solutions for parametric equation and the stable and unstable regions are evaluated. The results show that using a parametric excitation with an appropriate frequency and amplitude the system can be stabilized in the vicinity of the pitch fork bifurcation point. The time history and phase diagrams of the system are plotted for certain values of initial conditions and parameter values. Asymptotic analytically obtained results are verified by using direct numerical integration method.