大规模动力系统改进的快速精细积分方法

提出一种针对大规模动力系统的改进的快速精细积分方法（FPIM）。以精细积分方法为基础,利用大规模动力系统矩阵的稀疏性和动力问题的物理特性,分析了矩阵指数的特殊结构,并基于此给出一种计算大规模动力系统矩阵指数及其动力响应的高效率方法。     

Krasnoselskii不动点与中立型多变时滞随机动力系统的指数p稳定性

探讨了一类中立型多变时滞随机动力系统,并通过Krasnoselskii不动点方法,给出了该系统零解指数p稳定的条件;同时对所得结果进行了严格的证明,形成了中立型多变时滞随机动力系统零解指数p稳定性定理。本文首次采用Krasnoselskii不动点方法研究了一类中立型多变时滞随机动力系统零解的指数p稳定性,所得结果改进和推广了已有文献中的结论;并且根据多变时滞随机动力系统时滞τ_i(t)和δ_i(t)的特点,分别引入对应的函数h_i(s),i=1,2,L,n来构造算子,相比已有文献中的方法更加灵活实用。     

横向磁场激励下铁磁梁式板的混沌运动分析

在磁体力分布的磁弹性理论模型和磁场准静态假定模式基础上，对于处在周期时变磁场 中的不可移简支铁磁梁式板非线性磁弹性动力特性进行定性与定量分析．首先利用磁场的摄 动技术和结构变形的模态法，导出了关于模态坐标的非线性动力方程；然后利用Ｍｅｌｎｉｋｏｖ方 法，从理论上给出这一磁弹性动力系统可能出现混沌运动的必要条件及参数范围；最后采用变 步长Ｒｕｎｇｅ－Ｋｕｔｔａ数值积分方法对其磁弹性相互作用的混沌现象进行了定量搜索与模拟，并 利用其轨迹的Ｐｏｉｎｃａｒｅ截面图与Ｌｉａｐｕｎｏｖ指数加以判断．结果表明磁弹性简支梁式板在横 向周期时变磁场中存在混沌吸引子，且在机械阻尼很小时其混沌吸引子表现出稠的特性．     

非线性动力系统线性模型数值计算的Taylor变换法

将非线性动力系统化为连续变化的线性系统,并导出任意自治或非自治非线性动力系统的瞬时线性化方程,该线性方程的连续变化描述了系统的全部复杂动力行为.进一步采用Taylor变换法求解系统的线性化方程,得到一种非线性动力系统数值计算的新方法,避免了指数矩阵展开的乘积运算.计算实例表明该方法在不增加计算机时的前提下,精度高于传统的Houbolt,Wilson-θ及Newmark-β等方法.计算了Duffing方程和van　Pol方程的混沌及周期特性.     

有界随机噪声激励下碰撞系统的最大Lyapunov指数

为了研究单自由度线性单边碰撞系统在有界随机噪声参数激励下的最大 Lyapunov 指数和稳定性问题，用 Zhuravlev 变换将碰撞系统转化为连续的非碰撞系统，然后用随机平均法得到了关于慢变量的随机微分方程。在没有随机扰动的情形下，给出了系统最大Lyapunov指数的值；在有随机扰动的情形下，通过求解FPK方程得到了系统的不变测度和最大Lyapunov指数的解析表达式。研究结果表明：随着系统阻尼项、有界随机噪声带宽、碰撞恢复系数的减少和有界随机噪声振幅的增大，最大Lyapunov指数增加；当随机激励的中心频率等于系统固有频率的两倍时，系统的Lyapunov指数达到最大，从而使系统变得更不稳定。根据系统的Lyapunov指数得到了系统稳定的充分必要条件，即当Lyapunov指数大于零时系统几乎必然不稳定，而当Lyapunov指数小于零时系统几乎必然稳定，Lyapunov指数等于零为系统的稳定性分叉点，并讨论了相应的稳定性分叉问题。     

Liapunov函数的计算机构造及其应用

本文提出了一种易于编程借助计算机实现的构造线性离散时间系统Ｌｉａｐｕｎｏｖ函数的方法。该方法迥避了直接求解离散Ｌｉａｐｕｎｏｖ矩阵方程，而是利用相应矩阵的负定性条件，通过逐次改变解矩阵中某些元素之值，构造出系统的Ｌｉａｐｕｎｏｖ函数。作为实例，分析讨论了某鱼雷在深度控制系统操纵下纵平面运动的稳定性     

两级悬浮EMS型磁悬浮控制系统的非线性动力学特性

在考虑二级悬浮弹簧的非线性特性的基础上，建立了两级悬浮EMS型磁悬浮控制系统的非线性动力学模型，给出了控制参数G1，G2的稳定性条件，进一步讨论了该磁悬浮系统在外界激励下的分叉行为及混沌动力特性，并利用Poincare映射，功率谱分析及最大Lyapunov指数等混沌运动的统计特征描述了该状态下控制系统的混沌运动特性。     

周期结构动力响应的高效数值方法

提出一种计算周期结构动力响应的高效率算法. 以精细积分方法为基础, 利用周期结构的对称性和动力问题的物理特性, 分析了周期结构对应矩阵指数的特殊结构, 并基于此给出一种计算周期结构对应矩阵指数的高效率方法. 在高效和精确计算周期结构对应矩阵指数的基础上, 得到了周期结构动力响应的高效率和高精度算法. 数值算例表明, 该方法效率高且节省存储要求.     

多维磁浮柔性转子控制系统分岔与控制器设计

讨论了多维悬浮柔性转子控制系统局部及全局分岔问题,首先建立了该复杂系统动力学模型,应用中心流形和求规范形综合方法,得到此系统非半简双零特征值问题的规范形及其普适开折,并进一步讨论了此控制系统的分岔 行为（余维二分岔）及稳定性;给出了为实现稳定控制,控制器参数、转子系统结构参数的相互关系及稳定控制域,即给出分岔 参数条件、分岔曲线、转迁集,最后,给出此柔性转子控制系统的数值仿真结果。     

基于位移型Gurtin变分原理计算动力响应的逐步积分法

本文利用位移型Ｇｕｒｔｉｎ变分原理,在时间域上采用三次Ｈｅｒｍｉｔｅ插值函数进行离散,给出了一种计算结构动力响应的逐步积分方法。通过稳定性分析研究了该方法的稳定区情况表明,当１．６４≤θ≤２．０８时,该方法的数值计算精度很高,但是条件稳定积分格式。当θ≥４．１时,该方法是无条件稳定的积分格式,精度较高。     

Stability analyses on dynamic control for a moving body levitated magnetically over elastic guideways

Based on the Floquet theory on ordinary differential equations with periodically variable coefficients and the bifurcation approach to nonlinear equations, a numerical approach to determining the stability region of control parameters is established for a dynamic control system composed of a moving body levitated magnetically over flexible guideways. The system is nonlinearly coupled among the elastic deformation of guideways, disturbance the levitation position of the body and electromagnetic control forces. The numerical simulation is given for the system in the case in which the control parameters are taken within and out of the stability region respectively. Results show that this approach is reliable. This project is supported by NSFC(No. 19725207), the Pre-research Project of the Committee of Science and Technology for Defence of China and the Science Foundation of Education Ministry of China for Ph. D. Programmes.     

Almost sure stochastic stability of a viscoelastic double-beam system

This paper investigates the dynamic stability of a viscoelastic double-beam system under parametric excitations. It is assumed that the two beams, made from Voigt–Kelvin material, are simply supported and continuously joined by a Winkler elastic layer. Each pair of axial forces consists of a constant part and a time-dependent stochastic function. In the case of “non-white” excitations, by using the direct Liapunov method, bounds of the almost sure stability of the double-beam system as a function of retardation time, bending stiffness, stiffness modulus of the Winkler layer, variances of the stochastic forces and the intensity of the deterministic components of axial loading are obtained. Numerical calculations are performed for the Gaussian process with a zero mean, as well as a harmonic process with a random phase. When the excitations are wideband noises, almost sure stability is obtained within the concept of the Liapunov exponent. White noise and Ornstein–Uhlenbeck processes are considered as models of wideband noises.     

Stability of a spinning liquid-filled spacecraft

Summary The stability of a spinning liquid-filled spacecraft has been investigated in the present paper. Using Galerkin's method, the attitude dynamic equations have been given. The Liapunov direct method was employed to obtain a sufficient condition for stability. Three kinds of characteristic modals were investigated: free motion of inviscid fluid, slosh motion and non-slosh motion. All characteristic problems can be solved numerically by the Finite Element Method or the Boundary Element Method. It has been demonstrated that the viscosity of the fluid has a dissipative effect at large Reynolds number, while the slosh motion plays a destabilizing role. The non-slosh model of fluid does not affect the stability criterion. Accepted for publication 19 October 1996     

Dynamic stability of weak equations of continuous systems

The stability analysis method is developed for distributed dynamic problems with relaxed assumptions imposed on solutions. The problem is motivated by structural vibrations with external time-dependent parametric excitations which are controlled using surface-mounted or -embedded actuators and sensors. The strong form of equations involves irregularities which lead to computational difficulties for estimation and control problems. In order to avoid irregular terms resulting from differentiation of the force and moment terms the dynamics equations are written in a weak form. The weak form of dynamics equations of linear mechanical structures is obtained using variational calculus. The study of stability of stochastic weak system is based on examining properties of Liapunov functional along a weak solution.     

APPLICATION OF MECHANIZED MATHEMATICS TO ROTOR DYNAMICS

Based on the mechanized mathematics and WU Wen-tsun elimination method, using oil film forces of short-bearing model and Muszynska' s dynamic model, the dynamical behavior of rotor- bearing system and its .stability of motion are investigated . As example , the concept of Wu characteristic set and Maple software , whirl parameters of short- bearing model, which is usually solved by the numerical method, are analyzed. At the same time , stability of zero solution of Jeffcott rotor whirl equation and stability of self-excited vibration are studied. The conditions of stable motion are obtained by using theory of nonlinear vibration .     

Coupled flutter and divergence bifurcation of a double pendulum

The loss of stability of the equilibrium position of a double pendulum with follower force loading and elastic end support is studied. At a special parameter combination the linearized system is characterized by a zero root and a pure imaginary pair of eigenvalues. Therefore, the stability problem is a complicated critical case in the sense of Liapunov and requires a non-linear analysis. A complete post-bifurcation investigation of the coupled divergence and flutter motions is given by means of centre manifold theory, and bifurcation diagrams. Among the different types of motions even the appearance of chaotic behavior is shown.     

Influence of rotatory inertia on dynamic stability of the viscoelastic symmetric cross-ply laminated plates

The dynamic stability problem of the symmetrically laminated cross-ply plates made of the viscoelastic Voigt–Kelvin material, compressed by time-dependent stochastic membrane forces, is investigated. The effect of rotatory inertia is included in the present formulation. It is assumed that all elastic moduli have the same retardation times. By using the direct Liapunov method, bounds of the almost sure stability of cross-ply plates as a function of viscous damping coefficient, retardation time, variances of the stochastic forces, ratio of the principal lamina stiffnesses, number of layers, plate aspect ratio, cross-ply ratio and intensity of the deterministic components of axial loading are obtained. Numerical calculations are performed for the Gaussian process with a zero mean, as well as an harmonic process with random phase.     

Robust adaptive dynamic surface control design for a flexible air-breathing hypersonic vehicle with input constraints and uncertainty

The flight control problem of a flexible air-breathing hypersonic vehicle is presented in the presence of input constraint and aerodynamic uncertainty. A control-oriented model, where aerodynamic uncertainty and the strong couplings between the engine and flight dynamics are included, is derived to reduce the complexity of controller design. The flexible dynamics are viewed as perturbations of the model. They are not taken into consideration at the level of control design, the influence of which is evaluated through simulation. The control-oriented model is decomposed into velocity subsystem and altitude subsystem, which are controlled separately. Then robust adaptive controller is developed for the velocity subsystem, while the controller which combines dynamic surface control and radial basis function neural network is designed for the altitude subsystem. The unknown nonlinear function is approximated by the radial basis function neural network. Minimal-learning parameter technique is utilized to estimate the maximum norm of ideal weight vectors instead of their elements to reduce the computational burden. To handle input constraints, additional systems are constructed to analyze their impact, and the states of the additional systems are employed at the level of control design and stability analysis. Besides, “explosion of terms” problem in the traditional backstepping control is circumvented using a first-order filter at each step. By means of Lyapunov stability theory, it is proved theoretically that the designed control law can assure that tracking error converges to an arbitrarily small neighborhood around zero. Simulations are performed to demonstrate the effectiveness of the presented control scheme in coping with input constraint and aerodynamic uncertainty.     

First integral and stability of motion in the critical cases

In this paper, we indicate that after the Liapunov function by using linear combination of mechanical first integral was suggested by Chetayev in 1946. He and his students solved stability of conservative system by means of this method. But he had trouble to solve the problems by means of cut and try. Moreover, the condition of stability is imperfect. Solution by this method is limited for problems of purely imaginary roots. The cases of zero roots have not been considered. Condition of stability secured is more strict.This paper suggests that the differential equation can be transformed into standard form by method of cancellation of cyclic coordinates (method of lowering degree of order). and condition of stability can be determined by energy integral. By this method not only the computation is clear and concise. But also zero roots can be considered. Therefore the problems of two cyclic coordinates can be transformed into second-order system, and we get new conclusion of the condition of stability simply. As for problems of single cyclic coordinate, infact, Chetayev and his students did not solve the stability of the gyroscope of auter-gimbal with horizontal axis or arbitrary angle. In this paper, it shows that the method suggested here is useful for stability of these problems. The condition of conditional stability and instability were derived.     

Stability of elastic systems under a stochastic parametric excitation

An efficient method to investigate the stability of elastic systems subjected to the parametric force in the form of a random stationary colored noise is suggested. The method is based on the simulation of stochastic processes, numerical solution of differential equations, describing the perturbed motion of the system, and the calculation of top Liapunov exponents. The method results in the estimation of the almost sure stability and the stability with respect to statistical moments of different orders. Since the closed system of equations for moments of desired quantities y _j (t) cannot be obtained, the statistical data processing is applied. The estimation of moments at the instant t _n is obtained by statistical average of derived from the solution of equations for the large number of realizations. This approach allows us to evaluate the influence of different characteristics of random stationary loads on top Liapunov exponents and on the stability of system. The important point is that results found for filtered processes, are principally different from those corresponding to stochastic processes in the form of Gaussian white noises.