一类多体系统非完整运动最优规划的拟牛顿算法

研究带有非完整约束的一类多体系统运动规划问题。多体系统中的非完整约束通常是由不可积的速度约束或不可积的守恒律引起。在系统动量和动量矩守恒情况下,动力学方程降阶为非完整形式约束方程,系统的控制问题可转化为无漂移系统的非完整运动规划问题。文中首先导出具有多体开链系统的非完整运动模型。利用最优控制理论和最优化技术,采用输入参数化的方法将连续的最优控制问题转化为离散的最优控制问题,提出一种非完整多体系统运动规划的拟牛顿算法。最后将该方法用于自由漂浮的空间三连杆机构,仿真结果验证了该方法的有效性。     

Bifurcations and chaos in parametrically excited single-degree-of-freedom systems

The behavior of single-degree-of-freedom systems possessing quadratic and cubic nonlinearities subject to parametric excitation is investigated. Both fundamental and principal parametric resonances are considered. A global bifurcation diagram in the excitation amplitude and excitation frequency domain is presented showing different possible stable steady-state solutions (attractors). Fractal basin maps for fundamental and principal parametric resonances when three attractors coexist are presented in color. An enlargement of one region of the map for principal parametric resonance reveals a Cantor-like set of fractal boundaries. For some cases, both periodic and chaotic attractors coexist.     

柔性多体系统动力学研究现状与展望

对柔性多体系统动力学的研究现状进行了概括和总结,主要从柔性体建模方法、刚柔耦合动力学、接触碰撞问题、多物理场耦合、微分代数方程求解技术、控制方法、设计优化及软件开发和实验研究等几个研究方向进行总结,并对未来的研究方向做了展望.     

Pseudoelastic effect in autoparametric non-ideal vibrating system with SMA spring

In this paper a three degrees of freedom autoparametric system with limited power supply is investigated numerically. The system consists of the body, which is hung on a spring and a damper, and two pendulums connected by shape memory alloy (SMA) spring. Shape memory alloys have ability to change their material properties with temperature. A polynomial constitutive model is assumed to describe the behavior of the SMA spring. The non-ideal source of power adds one degree of freedom, so the system has four degrees of freedom. The equations of motion have been solved numerically and pseudoelastic effects associated with the martensitic phase transformation are studied. It is shown that in this type system one mode of vibrations might excite or damp another mode, and that except different kinds of periodic vibrations there may also appear chaotic vibrations. For the identification of the responses of the system's various techniques, including chaos techniques such as bifurcation diagrams and time histories, power spectral densities, Poincarè maps and exponents of Lyapunov may be used.     

Motion of an inclusion in uniformly and nonuniformly vibrating liquids

An approach to constructing the quantitative nonuniformity characteristics of liquid vibrations is proposed. A new problem of the motion of an inclusion in a vibrating liquid is considered. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 1, pp. 79–85, January–February, 2007.     

Non-existence of Shilnikov chaos in continuous-time systems

In this paper,a non-existence condition for homoclinic and heteroclinic orbits in n-dimensional,continuous-time,and smooth systems is obtained.Based on this result and an elementary example,it can be conjectured that there is a fourth kind of chaos in polynomial ordinary differential equation(ODE) systems characterized by the nonexistence of homoclinic and heteroclinic orbits.     

Bifurcation and chaos of airfoil with multiple strong nonlinearities

The bifurcation and chaos phenomena of two-dimensional airfoils with multiple strong nonlinearities are investigated.First,the strongly nonlinear square and cubic plunging and pitching stiffness terms are considered in the airfoil motion equations,and the fourth-order Runge-Kutta simulation method is used to obtain the numerical solutions to the equations.Then,a post-processing program is developed to calculate the physical parameters such as the amplitude and the frequency based on the discrete numerical solutions.With these parameters,the transition of the airfoil motion from balance,period,and period-doubling bifurcations to chaos is emphatically analyzed.Finally,the critical points of the period-doubling bifurcations and chaos are predicted using the Feigenbaum constant and the first two bifurcation critical values.It is shown that the numerical simulation method with post-processing and the prediction procedure are capable of simulating and predicting the bifurcation and chaos of airfoils with multiple strong nonlinearities.     

Stability chart of parametric vibrating systems using energy-rate method

Based on the integral of energy and numerical integration, we introduce, develop, and apply a general algorithm to predict parameters of a parametric equation to produce a periodic response. Using the new method, called energy-rate, we are able to find not only stability chart of a parametric equation which indicates the boundaries of stable and unstable regions, but also periodic responses that are embedded in stable or unstable regions.There are three main important advantages in energy-rate method. It can be applied not only to linear but also to non-linear parametric equations; most of the perturbation methods cannot. It can be applied to large values of parameters; most of the perturbation methods cannot. Depending on the accuracy of numerical integration method, it can also find the value of parameters for a periodic response more accurate than classical methods, no matter if the periodic response is on the boundary of stability and instability or it is a periodic response within the stable or unstable region.In order to introduce the energy-rate method and indicate its advantages we apply the method to the standard Mathieu's equation,

     

多体系统动力学方程违约修正的数值计算方法

多体系统动力学方程为微分代数方程,一般将其转化成常微分方程组进行数值计算,在数值积分的过程中约束方程的违约会逐渐增大.本文对具有完整、定常约束的多体系统,在修改的带乘子Lagrange正则形式的方程的基础上,根据Baumgarte提出的违约修正的方法,给出了一种多体系统微分代数方程违约修正法和系统的动力学方程的矩阵表达式.通过对曲柄-滑块机构的数值仿真,计算结果表明本文给出的方法在计算精度和计算效率上好于Baumgarte提出的两种违约修正的方法.     

Numerical investigation of a response probability density function of stochastic vibroimpact systems with inelastic impacts

The paper is devoted to numerical calculations of a response probability density function of stochastic vibroimpact systems excited by additive Gaussian white noise. An impact in the considered systems is modeled as a classical, inelastic one against motionless barrier(s) with values of a restitution coefficient both very small and close to unity. Certain empirical formulas are derived based on the results of numerical simulation.     

Parametric open-plus-closed-loop control of chaos in continuous dynamical systems

This paper presents a parametric open-plus-closed-loop control approach to controlling chaos in continuous dynamical systems. As an example, chaos in the Lorenz model is controlled to demonstrate its application. Finally, the relations between the parametric open-plus-closed-loop control and the former control methods, such as the open-plus-closed-loop control and the parametric entrainment control, are discussed.Supported by the Science Foundation of the State Education Commission for Doctorate Program, and the Applied Science Foundation of the State Ministry of Metallurgical Industry.     

Wavelet-Based Analysis of Parametric Vibrations of Flexible Plates

Wavelet-based methods are reviewed, and their advantages in comparison to standard approaches are outlined. Then, a wavelet analysis is performed to investigate parametric vibrations of flexible plates under a sinusoidal load. In particular, a scenario leading from regular motion to chaos is analyzed     

Synchronisation and Chaos in a Parametrically and Self-Excited System with Two Degrees of Freedom

Vibration analysis of a non-linear parametrically andself-excited system of two degrees of freedom was carried out. The modelcontains two van der Pol oscillators coupled by a linear spring with a aperiodically changing stiffness of the Mathieu type. By means of amultiple-scales method, the existence and stability of periodicsolutions close to the main parametric resonances have beeninvestigated. Bifurcations of the system and regions of chaoticsolutions have been found. The possibility of the appearance ofhyperchaos has also been discussed and an example of such solution hasbeen shown.     

Nonlinear flexural waves and chaos behavior in finite-deflection Timoshenko beam

Based on the Timoshenko beam theory, the finite-deflection and the axial inertia are taken into account, and the nonlinear partial differential equations for flexural waves in a beam are derived. Using the traveling wave method and integration skills, the nonlinear partial differential equations can be converted into an ordinary differential equation. The qualitative analysis indicates that the corresponding dynamic system has a heteroclinic orbit under a certain condition. An exact periodic solution of the nonlinear wave equation is obtained using the Jacobi elliptic function expansion. When the modulus of the Jacobi elliptic function tends to one in the degenerate case, a shock wave solution is given. The small perturbations are further introduced, arising from the damping and the external load to an original Hamilton system, and the threshold condition of the existence of the transverse heteroclinic point is obtained using Melnikov＇s method. It is shown that the perturbed system has a chaotic property under the Smale horseshoe transform.     

Experimental study on dynamics of an oblique-impact vibrating system of two degrees of freedom

The paper presents a detailed experimental study of an oblique-impact vibration system of two degrees of freedom. The primary objective of the study is to verify the hypothesis of instantaneous impact in the oblique-impact process of two elastic bodies such that the incremental impulse method works for computing the nonlinear dynamics of the oblique-impact vibrating systems. The experimental setup designed for the objective consists of a harmonically excited oscillator and a pendulum, which obliquely impacts the oscillator. In the study, the dynamic equation of the experimental setup was established first, and then the system dynamics was numerically simulated by virtue of the incremental impulse method. Afterwards, rich dynamic phenomena, such as the periodic vibro-impacts, chaotic vibro-impacts and typical bifurcations, were observed in a series of experiments. The comparison between the experimental results and the numerical simulations indicates that the incremental impulse method is reasonable and successful to describe the dynamics during an oblique-impact process of two elastic bodies. The study also shows the limitation of the hypothesis of instantaneous impact in an oblique-impact process. That is, the hypothesis only holds true in the case when the impact angle is not too large and the relative approaching velocity in the normal direction is not too low. Furthermore, the paper gives the analysis of the tangential rigid-body slip on the contact surface in the case of a large impact angle, and explains why there exist some discrepancies between the numerical simulations and the experimental results.     

Impact phenomena of rotor-casing dynamical systems

Rubbing and impacting between a rotor and adjacent motion-constraining structures is a serious malfunction in rotating machinery. A shaver rotor-casing system with clearance and mass imbalance is modelled with two second-order ordinary differential equations and inelastic impact conditions. The dynamics is investigated analytically, as well as by numerical simulation. A Lyapunov exponent technique is developed to characterize the topologically different behavior as the parameters are varied. The dry friction coefficient and the eccentricity of the rotor imbalance were chosen to be the two variable parameters, the effect of which on the system dynamics is illustrated through phase plots, bifurcation diagrams, as well as Poincaré maps. The results demonstrate the existence of both rubbing and impacting behavior. Depending on values of the parameters, rubbing motion in both the clockwise and counter-clockwise directions may occur. Within the impact regime, the impact behavior could be periodic, quasi-periodic or chaotic, as confirmed by the calculation of Lyapunov exponents.     

弹性杆-刚性体系统动力学分析的广义逆矩阵方法

将多刚体系统的广义逆矩阵方法推广到含弹性杆与刚性体的混合系统的动力学分析中,建立了以节点坐标表示的基于全局惯性坐标系的刚体-柔性体混合系统动力学方程.首先以两端节点坐标为变量推导了弹性杆的动力学方程,以刚性体节点坐标为变量推导了刚性体节点速度约束方程和刚性体动力学方程,最后得到弹性杆与刚性体混合系统的动力学方程和速度约束方程.本方法在同一个惯性坐标系对刚柔多体系统进行描述,具有方法简洁、便于计算建模的特点.论文最后给出两个数值算例,检验了方法的有效性.     

Vibrations and self-heating of a layered elastic-viscoelastic rectangular prism under a vibrating punch

The vibrations and self-heating of a layered metal-polymer rectangular prism kinematically excited by a normally vibrating punch are analyzed numerically. The effect of reinforcement layers on the rate of variation in and the spatial distribution of the temperature field is examined. It is established that thin metallic layers in a polymer matrix change the dissipation mechanism from shear to bulk. The effect of contact stress concentration on the localization of temperature field is studied. It is established that the energy flux through the load application area decreases due to softening of the material. Several reinforcement effects are revealed. It is shown that they are determined by the value of the load parameter. The numerical results are in good agreement with experimental data on self-heating kinetics __________ Translated from Prikladnaya Mekhanika, Vol. 43, No. 8, pp. 71–79, August 2007.     

Piecewise analysis of oblique vibro-impacting systems

A great number of studies have shown the complex nonlinear dynamics of mechanical systems with repeated normal impacts. An oblique frictional impact introduces even more complicated dynamics such as stick-slip motions to those systems. Hence, the dynamics of oblique vibro-impacting systems with possible sliding motion is an open problem. Based on a hybrid analysis of vibro-impact dynamics, kinematics and complementary conditions, a piecewise analysis method is developed in the paper to describe the sliding motion during an oblique impact. Thereby, a parametrically excited planar pendulum between two parallel rigid walls is studied as an illustrative example. The example, together with the corresponding numerical results, shows that the sliding impacts occur in such a system with a set of properly selected parameters. The project supported by the National Natural Science Foundation of China (59905010) and by the Scientific Research Foundation for the Returned Overseas Chinese Scholars, Ministry of Education, China