求非线性动力系统周期解的切比雪夫多项式法

周期运动是一种在客观世界中普遍存在的运动形式,它与混沌运动之间存在十分密切的关系,因而具有很重要的研究价值。利用切比雪夫多项式的若干良好性质,对自治非线性动力系统进行分析,将状态矢量在主周期上展开为切比雪夫多项式的形式,从而将原问题转变为非线性代数方程组的求解问题,得出一种可以方便、迅速地获得周期轨道近似多项式表达式的方法。此方法不依赖于小参数假设,可以用于分析强非线性问题,而且对参数激励系统同样有效。在计算机条件允许时,对高维系统也能迅速、精确地得到其周期轨道的近似多项式表达式。以三维Rossler系统和五维非线性磁浮转子系统周期轨道的计算为例,通过与四阶Runge-Kutta数值积分结果比较,说明此方法的精确、高效性。     

Analytical solutions for dynamics of dual-spin spacecraft and gyrostat-satellites under magnetic attitude control in omega-regimes

The attitude dynamics of a dual-spin spacecraft (a gyrostat with one rotor) with magnetic actuators attitude control is considered in the constant external magnetic field at the presence of the spacecraft’s own magnetic dipole moment, which is created proportionally to the angular velocity components (this motion regime can be called as “the omega-regime” or “the omega-maneuver”). The research of the dual-spin spacecraft angular motion under the action of the magnetic restoring torque is fulfilled in the generalized formulation close to the classical mechanics’ task of the heavy body/gyrostat motion in the Lagrange top. Analytical exact solutions of differential equations of the motion are obtained for all parameters in terms of elliptic integrals and the Jacobi functions. New obtained analytical solutions can be classified as results developing the classical fundamental problem of the rigid body and gyrostat motion around the fixed point. The technical application of the omega-regime to the angular reorientation of the spacecraft longitudinal axis along the angular momentum vector is considered.     

Energy harvesting from a magnetic levitation system

Numerical and experimental studies of a magnetic levitation harvester are presented in the paper. The idea is based on the motion of permanent cylinder magnet in a coil exploited for energy harvesting. The novel model is based on a new definition of the coupling coefficient (inductive coefficient) which relates mechanical and an electrical components. The performed static and dynamics experimental tests show that this coefficient is a nonlinear function of the magnet position, and highly depends on the magnet coordinate in the coil, in such a way that the maximum energy is obtained in a coil ends. The comparison between classical – fixed value model – and novel nonlinear model of the inductive coefficient is presented for selected cases. The most essential differences are presented.     

Steady-state analysis of multibody systems with reference to vehicle dynamics

This paper presents a systematic methodology and formulation for determining the steady-state response of multibody systems. The equations of motion for a general multibody system are described in terms of a set of relative joint accelerations. Then, the differential equations of motion are converted to a set of algebraic equations for the steady-state response. These equations are derived based upon a set of conditions that must exist for the steady state. The application of this formulation in determining the steady-state response of a vehicle moving in a circular path is shown. The multibody model of the vehicle for two- or four-wheel steering is presented. The results of the steady-state simulation are compared with those obtained from a transient dynamic analysis.     

Control structure interaction in the nonlinear analysis of flexible mechanical systems

The effect of the control structure interaction on the feedforward control law as well as the dynamics of flexible mechanical systems is examined in this investigation. An inverse dynamics procedure is developed for the analysis of the dynamic motion of interconnected rigid and flexible bodies. This method is used to examine the effect of the elastic deformation on the driving forces in flexible mechanical systems. The driving forces are expressed in terms of the specified motion trajectories and the deformations of the elastic members. The system equations of motion are formulated using Lagrange's equation. A finite element discretization of the flexible bodies is used to define the deformation degrees of freedom. The algebraic constraint equations that describe the motion trajectories and joint constraints between adjacent bodies are adjoined to the system differential equations of motion using the vector of Lagrange multipliers. A unique displacement field is then identified by imposing an appropriate set of reference conditions. The effect of the nonlinear centrifugal and Coriolis forces that depend on the body displacements and velocities are taken into consideration. A direct numerical integration method coupled with a Newton-Raphson algorithm is used to solve the resulting nonlinear differential and algebraic equations of motion. The formulation obtained for the flexible mechanical system is compared with the rigid body dynamic formulation. The effect of the sampling time, number of vibration modes, the viscous damping, and the selection of the constrained modes are examined. The results presented in this numerical study demonstrate that the use of the driving forees obtained using the rigid body analysis can lead to a significant error when these forces are used as the feedforward control law for the flexible mechanical system. The analysis presented in this investigation differs significantly from previously published work in many ways. It includes the effect of the structural flexibility on the centrifugal and Coriolis forces, it accounts for all inertia nonlinearities resulting from the coupling between the rigid body and elastic displacements, it uses a precise definition of the equipollent systems of forces in flexible body dynamics, it demonstrates the use of general purpose multibody computer codes in the feedforward control of flexible mechanical systems, and it demonstrates numerically the effect of the selected set of constrained modes on the feedforward control law.     

Elimination of periodic coefficients from the equations of motion of asymmetric shafts in anisotropic bearings

Summary A method is given for eliminating periodic coefficients from rotordynamic systems in which imperfections are accounted for simultaneously in the stationary and rotating components. It is well known that periodic coefficients are unavoidable in such systems whether the equations of motion are formulated in the stationary frame or the rotating frame of reference. Since such equations are more difficult to solve compared to those with constant coefficients, it is desirable to devise a transformation procedure which eliminates periodic coefficients from the governing equations. In this paper, such a procedure is presented. A considerable simplification is thereby obtained of the equations of motion of rotating machinery in which imperfections exist stators and rotors.     

具有广义力场损伤粘弹性梁-柱的广义变分原理及应用

本文从考虑损伤的粘弹性材料的卷积型本构关系出发,建立了在小变形下损伤粘弹性梁-柱的控制方程。提出了以卷积形式表示的梁-柱弯曲问题的泛函,并给出了损伤粘弹性梁-柱的广义变分原理。应用这个广义变分原理,可分别给出梁-柱位移和损伤满足的基本方程,以及相应的初始条件和边界条件。应用Galerkin截断和非线性动力学的数值分析方法,分析了两端简支损伤粘弹性梁柱的动力学行为,给出了不同的材料参数对系统响应的影响。     

超导电磁悬浮力的数值模拟

利用超导体宏观电磁场的Kim模型对高温超导体－永磁铁悬浮系统中二才之间的电磁悬浮力进行了定量计算与模拟。通过与已有实验结果比较表明：对超导悬浮力的数值模拟能较好地反映悬浮力的主要特征。在此基础上，给出了超导体内屏蔽电流，铅直悬浮力和面内电磁力的分布情况，数值仿真结果显示了悬浮力的大小与永磁铁的表磁强度和尺寸有关。这些结果为提高超导悬浮力和改进超导宏观电磁场模型提供了理论依据。     

An index 0 differential-algebraic equation formulation for multibody dynamics: Nonholonomic constraints

A method is presented for formulating and numerically integrating index 0 differential-algebraic equations of motion for multibody systems with holonomic and nonholonomic constraints. Tangent space coordinates are defined in configuration and velocity spaces as independent generalized coordinates that serve as state variables in the formulation. Orthogonal dependent coordinates and velocities are used to enforce position, velocity, and acceleration constraints to within specified error tolerances. Explicit and implicit numerical integration algorithms are presented and used in solution of three examples: one planar and two spatial. Numerical results verify that accurate results are obtained, satisfying all three forms of kinematic constraint to within error tolerances embedded in the formulation.     

Magneto-elastic combination resonances analysis of current-conducting thin plate

Based on the Maxwell equations, the nonlinear magneto-elastic vibration equations of a thin plate and the electrodynamic equations and expressions of electro- magnetic forces are derived. In addition, the magneto-elastic combination resonances and stabilities of the thin beam-plate subjected to mechanical loadings in a constant transverse magnetic filed are studied. Using the Galerkin method, the corresponding nonlinear vibration differential equations are derived. The amplitude frequency response equation of the system in steady motion is obtained with the multiple scales method. The excitation condition of combination resonances is analyzed. Based on the Lyapunov stability theory, stabilities of steady solutions are analyzed, and critical conditions of stability are also obtained. By numerical calculation, curves of resonance-amplitudes changes with detuning parameters, excitation amplitudes and magnetic intensity in the first and the second order modality are obtained. Time history response plots, phase charts, the Poincare mapping charts and spectrum plots of vibrations are obtained. The effect of electro-magnetic and mechanical parameters for the stabilities of solutions and the bifurcation are further analyzed. Some complex dynamic performances such as period- doubling motion and quasi-period motion are discussed.     

Dynamics of vehicles with high gravity centre

The vehicles with high gravity centre are more prone to roll over. The paper deals with a method of dynamics analysis of fire engines which is an example of these types of vehicle. Algorithms for generating the equations of motion have been formulated by homogenous transformations and Lagrange's equation. The model presented in this article consists of a system of rigid bodies connected one with another forming an open kinematic chain. Road maneuvers such as a lane change and negotiating a circular track have been presented as the main simulations when a car loses its stability. The method has been verified by comparing numerical results with results obtained by experimental measurements performed during road tests.     

Changes in the natural frequency of a ferromagnetic rod in a magnetic field due to magnetoelastic interaction

Based on the magnetoelastic generalized variational principle and Hamilton＇s principle, a dynamic theoretical model characterizing the magnetoelastic interaction of a soft ferromagnetic medium in an applied magnetic field is developed in this paper. From the variational manipulation of magnetic scale potential and elastic displacement, all the fundamental equations for the magnetic field and mechanical deformation, as well as the magnetic body force and magnetic traction for describing magnetoelastic interaction are derived. The theoretical model is applied to a ferromagnetic rod vibrating in an applied magnetic field using a perturbation technique and the Galerkin method. The results show that the magnetic field will change the natural frequencies of the ferromagnetic rod by causing a decrease with the bending motion along the applied magnetic field where the magnetoelastic buckling will take place, and by causing an increase when the bending motion of the rod is perpendicular to the field. The prediction by the mode presented in this paper qualitatively agrees with the natural frequency changes of the ferromagnetic rod observed in the experiment.     

Nonlinear vibration of a rotating cantilever beam in a surrounding magnetic field

The nonlinear vibrations of a rotating cantilever beam made of magnetoelastic materials surrounded by a uniform magnetic field are investigated. The kinetic energy, potential energy and work done by the electromagnetic force are obtained. A nonlinear dynamic model, based on the Hamilton principle, which includes the stretching vibration and bending vibration is presented. The Galerkin method is adopted to discretize the dynamic equations. The proposed method is validated by comparison with the literature. The nonlinear behaviors of the responses are studied. Then simulations for different kinds of magnetic field are conducted. The effects of magnetic field parameters, including the amplitude, plane angle, spatial angle and time-varying frequency, on the dynamic behaviors of the stretching motion and bending motion are investigated in detail. The results illustrate that the interaction effects between the rotating cantilever beam and the magnetic field will increase the vibration amplitude and fluctuation of the beam. In particular, we found that: collinear magnetic fields with equal amplitude lead to the same dynamic responses; the amplitude of magnetic field intensity increases the dynamic responses remarkably; the response amplitude changes nonlinearly with the plane angle and spatial angle of the magnetic field; and the increase of time-varying frequency enhances dynamic responses of the rotating cantilever beam.     

大型偏心结构吊装欠约束问题研究

吊装施工过程中被吊模块的水平度是作业要求的重要指标,通常需要增加配重调平。传统有限元方法需要补充约束以消除单元刚体位移,且需要重复计算平衡方程来求解调平载荷,效率不高。将模块的运动分解为随动坐标系的整体运动以及相对该坐标系的弹性变形,可将欠约束问题化为多体系统的静平衡问题。基于虚功率原理推导了吊装平顺时刻的节点力平衡方程以及相应的切线刚度矩阵,并将配重表示为基础配重与载荷系数相乘的形式。通过对节点力平衡方程求导,得到一组以载荷系数为自变量的微分方程,通过求解微分方程并结合水平度判据,可快速搜寻满足水平度要求的载荷系数。数值算例表明,该方法在解决偏心模块吊装欠约束问题方面具有明显的优势,在确定配重载荷方面具有较快的速度和合理的精度。     

Generalized dynamic model for multibodies manipulator

In this paper the general dynamical equations were given for multibodiesmanipulator.The system is a topologic tree structure consisting of arbitrary number ofrigid bodies.The hinges allow the rotational and/or tranlational motion.Inconsideration of influence of friction the dynamic equations are established by meansof Newton-Euler’s method.Further,the equations are separated by way ofconstructing the distribution matrices and a group of force and motion equations areobtained.     

In-plane wave motion in finite element model

The analysis method of lattice dynamics in classical physics is extended to study the properties of in-plane wave motion in the hybrid-mass finite element model in this paper. The dispersion equations of P and SV waves in the discrete model are first obtained by means of separating the characteristic equation of the motion equation, and then used to analyse the properties of P-and SV-homogeneous, inhomogeneous waves and other types of motion in the model. The dispersion characters, cut-off frequencies of P and SV waves, the polarization drift and appendent anisotropic property of wave motion caused by the discretization are finally discussed. The project sponsored by the Earthquake Science Foundation under Contract No. 90141     

Lubricated revolute joints in rigid multibody systems

The main purpose of this work is to present a general methodology for modeling lubricated revolute joints in constrained rigid multibody systems. In the dynamic analysis of journal-bearings, the hydrodynamic forces, which include both squeeze and wedge effects, generated by the lubricant fluid, oppose the journal motion. The hydrodynamic forces are obtained by integrating the pressure distribution evaluated with the aid of Reynolds’ equation, written for the dynamic regime. The hydrodynamic forces built up by the lubricant fluid are evaluated from the system state variables and included into the equations of motion of the multibody system. Numerical examples are presented in order to demonstrate the use of the methodologies and procedures described in this work.     

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

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

Dynamics of nanofibres conveyed by low Reynolds number flow in a microchannel

In this paper we aim to create an experimental and numerical model of nano and micro filaments suspended in a confined Poiseuille flow. The experimental data obtained for short nanofibres will help to elucidate fundamental questions concerning mobility and deformation of biological macromolecules due to hydrodynamic stresses from the surrounding fluid motion. Nanofibres used in the experiments are obtained by electrospinning polymer solutions. Their typical dimensions are 100–1000 μm (length) and 0.1–1 μm (diameter). The nanofibre dynamics is followed experimentally under a fluorescence microscope. A precise multipole expansion method of solving the Stokes equations, and its numerical implementation are used to construct a bead-spring model of a filament moving in a Poiseuille flow between two infinite parallel walls. Simulations show typical behaviour of elongated macromolecules. Depending on the parameters, folding and unfolding sequences of a flexible filament are observed, or a rotational and translation motion of a shape-preserving filament. An important result of our experiments is that nanofibres do not significantly change their shape while interacting with a micro-flow. It appeared that their rotational motion is better reproduced by the shape-preserving Stokesian bead model with all pairs of beads connected by springs, omitting explicit bending forces.