具有可积微分约束的力学系统的Lie对称性

研究具有可积微分约束的力学系统的Ｌｉｅ对称性与守恒量。采用两种方法：一是用不可积微分约束系统的方法;另一是用积分后降阶系统的方法,研究两种方法之间的关系。     

On the brachistochronic motion of mechanical systems with unilateral constraints

The problem of brachistochronic motion of a two-degree of freedom mechanical system is considered. The ideal bilateral constraints and one unilateral constraint with Coulomb friction are imposed on the system. It has been assumed that the components of the metric tensor are constant. After introducing the transformation of generalized coordinates, whereby the quadratic form of kinetic energy is reduced to a canonical form, the considered problem is geometrically interpreted as the motion of a representative point in a two-dimensional Euclidean space. Thus, the analogy with the brachistochronic motion of a particle along a curve which is treated as a unilateral constraint is achieved. The results are illustrated via an example.     

On the Noether symmetry and Lie symmetry of mechanical systems

The Noether symmetry is an invariance of Hamilton action under infinitesimal transformations of time and the coordinates. The Lie symmetry is an invariance of the differential equations of motion under the transformations. In this paper, the relation between these two symmetries is proved definitely and firstly for mechanical systems. The results indicate that all the Noether symmetries are Lie symmetries for Lagrangian systems meanwhile a Noether symmetry is a Lie symmetry for the general holonomic or nonholonomic systems provided that some conditions hold. The project supported by the National Natural Science Foundation of China (19972010)     

Form invariance and Lie symmetry of variable mass nonholonomic mechanical system

IntroductionThereisacloserelationbetweenthesymmetryandtheconservedquantityinamechanicalsystem .ModernmethodstofindconservedquantityofamechanicalsystemaremainlyNoethersymmetrymethod[1]andLiesymmetrymethod[2 ].NoethersymmetryisaninvarianceoftheHamiltonactionundertheinfinitesimaltransformations.Liesymmetryisaninvarianceofthedifferentialequationsundertheinfinitesimaltransformations.Inthepasttenyears,aseriesofimportresultshavebeenobtainedonthestudyoftheNoethersymmetryandLiesymmetry[3～12 ].Thefo…     

动力学与控制论力学系统的自由度

力学系统的自由度定义源自描述系统位形的独立坐标数.在分析力学发展过程中,人们通过对非完整约束的研究,将其拓展为独立的坐标 变分数.本文指出,对于含非完整约束的力学系统,该定义存在不妥之处,给出的自由度会过度限制系统的力学行为.文中研究力学系统在状态空间中的可达流形,指出可达流形维数与描述系统动力学的一阶常微分方程组的最少未知函数个数一致,例如Gibbs-Appell方程与广义速度方程联立的未知函数个数,进而将可达流形维数的一半定义为系统自由度.通过含黏弹性支承的振动系统、在倾斜平面上运动的冰橇等案例,讨论了单个非完整约束导致的半自由度概念,指出其力学意义和与相邻整数自由度的关系.此外,文中还给出两个非完整约束导致系统减少一个自由度的案例,讨论了系统的切丛和余切丛维数.      

An Augmented Formulation for Mechanical Systems with Non-Generalized Coordinates: Application to Rigid Body Contact Problems

In computational multibody algorithms, the kinematic constraintequations that describe mechanical joints and specified motiontrajectories must be satisfied at the position, velocity andacceleration levels. For most commonly used constraint equations, onlyfirst and second partial derivatives of position vectors with respect tothe generalized coordinates are required in order to define theconstraint Jacobian matrix and the first and second derivatives of theconstraints with respect to time. When the kinematic and dynamicequations of the multibody systems are formulated in terms of a mixedset of generalized and non-generalized coordinates, higher partialderivatives with respect to these non-generalized coordinates arerequired, and the neglect of these derivatives can lead to significanterrors. In this paper, the implementation of a contact model in generalmultibody algorithms is presented as an example of mechanical systemswith non-generalized coordinates. The kinematic equations that describethe contact between two surfaces of two bodies in the multibody systemare formulated in terms of the system generalized coordinates and thesurface parameters. Each contact surface is defined using twoindependent parameters that completely define the tangent and normalvectors at an arbitrary point on the body surface. In the contact modeldeveloped in this study, the points of contact are searched for on lineduring the dynamic simulation by solving the nonlinear differential andalgebraic equations of the constrained multibody system. It isdemonstrated in this paper that in the case of a point contact andregular surfaces, there is only one independent generalized contactconstraint force despite the fact that five constraint equations areused to enforce the contact conditions.     

A velocity transformation method for the nonlinear dynamic simulation of railroad vehicle systems

The solution of the constrained multibody system equations of motion using the generalized coordinate partitioning method requires the identification of the dependent and independent coordinates. Using this approach, only the independent accelerations are integrated forward in time in order to determine the independent coordinates and velocities. Dependent coordinates are determined by solving the nonlinear constraint equations at the position level. If the constraint equations are highly nonlinear, numerical difficulties can be encountered or more Newton–Raphson iterations may be required in order to achieve convergence for the dependent variables. In this paper, a velocity transformation method is proposed for railroad vehicle systems in order to deal with the nonlinearity of the constraint equations when the vehicles negotiate curved tracks. In this formulation, two different sets of coordinates are simultaneously used. The first set is the absolute Cartesian coordinates which are widely used in general multibody system computer formulations. These coordinates lead to a simple form of the equations of motion which has a sparse matrix structure. The second set is the trajectory coordinates which are widely used in specialized railroad vehicle system formulations. The trajectory coordinates can be used to obtain simple formulations of the specified motion trajectory constraint equations in the case of railroad vehicle systems. While the equations of motion are formulated in terms of the absolute Cartesian coordinates, the trajectory accelerations are the ones which are integrated forward in time. The problems associated with the higher degree of differentiability required when the trajectory coordinates are used are discussed. Numerical examples are presented in order to examine the performance of the hybrid coordinate formulation proposed in this paper in the analysis of multibody railroad vehicle systems.     

Group classification for path equation describing minimum drag work and symmetry reductions

The path equation describing the minimum drag work first proposed by Pakdemirli is reconsidered （Pakdemirli, M. The drag work minimization path for a fly- ing object with altitude-dependent drag parameters. Proceedings of the Institution of Mechanical Engineers, Part C, Journal of Mechanical Engineering Science 223（5）, 1113- 1116 （2009））. The Lie group theory is applied to the general equation. The group classi- fication with respect to an altitude-dependent arbitrary function is presented. Using the symmetries, the group-invariant solutions are determined, and the reduction of order is performed by the canonical coordinates.     

Generating functions for volume-preserving transformations

A general implicit solution for determining volume-preserving transformations in the n-dimensional Euclidean space is obtained in terms of a set of 2n generating functions in mixed coordinates. For n=2, the proposed representation corresponds to the classical definition of a potential stream function in a canonical transformation. For n=3, the given solution defines a more general class of isochoric transformations, when compared to existing methods based on multiple potentials. Illustrative examples are discussed both in rectangular and in cylindrical coordinates for applications in mechanical problems of incompressible continua. Solving exactly the incompressibility constraint, the proposed representation method is suitable for determining three-dimensional isochoric perturbations to be used in bifurcation theory. Applications in non-linear elasticity are envisaged for determining the occurrence of complex instability patterns for soft elastic materials.     

A matrix formulation for the dynamic analysis of planar mechanisms using point coordinates and velocity transformation

This paper presents a matrix formulation for the dynamic analysis of planar mechanisms consisting of interconnected rigid bodies. The formulation initially uses the rectangular Cartesian coordinates of an equivalent constrained system of particles to define the configuration of the mechanical system. This results in a simple and straightforward procedure for generating the equations of motion. The equations of motion are then derived in terms of relative joint coordinates through the use of a velocity transformation matrix. The velocity transformation matrix relates the relative joint velocities to the Cartesian velocities. For the open-loop case, this process automatically eliminates all of the non-working constraint forces and leads to an efficient integration of the equations of motion. For the closed-loop case, suitable joints should be cut and few cut-joints constraint equations should be included for each closed loop. Two examples are used to demonstrate the generality and efficiency of the proposed method.     

General dynamical systems described by first order quasilinear PDEs reducible to homogeneous and autonomous form

Within the framework of Lie group analysis of differential equations, a theorem is determined stating necessary and sufficient conditions allowing one to recover an invertible point transformation mapping a general dynamical system described by nonhomogeneous and nonautonomous first order quasilinear partial differential equations to homogeneous and autonomous form. The proof of the theorem is constructive and the new independent and dependent variables are obtained by determining the canonical variables associated to a suitable subalgebra of the Lie algebra of point symmetries admitted by the source system. The theorem is applied by considering some examples of physical interest arising from different contexts.     

Impact dynamics of multibody systems with frictional contact using joint coordinates and canonical equations of motion

This paper presents a methodology in computational dynamics for the analysis of mechanical systems that undergo intermittent motion. A canonical form of the equations of motion is derived with a minimal set of coordinates. These equations are used in a procedure for balancing the momenta of the system over the period of impact, calculating the jump in the body momentum, velocity discontinuities and rebounds. The effect of dry friction is discussed and a contact law is proposed. The present formulation is extended to open and closed-loop mechanical systems where the jumps in the constraints' momenta are also solved. The application of this methodology is illustrated with the study of impact of open-loop and closed-loop examples.     

Contact-impact formulation for multi-body systems using component mode synthesis

The efficiency and accuracy are two most concerned issues in the modeling and simulation of multi-body systems involving contact and impact. This paper proposed a formulation based on the component mode synthesis method for planar contact problems of flexible multi-body systems. A flexible body is divided into two parts： a contact zone and an un-contact zone. For the un-contact zone, by using the fixed-interface substructure method as reference, a few low-order modal coordinates are used to replace the nodal coordinates of the nodes, and meanwhile, the nodal coordinates of the local impact region are kept unchanged, therefore the total degrees of freedom （DOFs） are greatly cut down and the computational cost of the simulation is significantly reduced. By using additional constraint method, the impact constraint equations and kinematic constraint equations are derived, and the Lagrange equations of the first kind of flexible multi-body system are obtained. The impact of an elastic beam with a fixed half disk is simulated to verify the efficiency and accuracy of this method.     

Weak formulation and first order form of the equations of motion for a class of constrained mechanical systems

Some new theoretical results are presented on modeling the dynamic response of a class of discrete mechanical systems subject to equality motion constraints. Both the development and presentation are facilitated by employing some fundamental concepts of differential geometry. At the beginning, the equations of motion of the corresponding unconstrained system are presented on a configuration manifold with general properties, first in strong and then in a primal weak form, using Newton׳s law of motion as a foundation. Next, the final weak form is obtained by performing a crucial integration by parts step, involving a covariant derivative. This step required the clarification and enhancement of some concepts related to the variations employed in generating the weak form. The second part of this work is devoted to systems involving holonomic and non-holonomic scleronomic constraints. The equations of motion derived in a recent study of the authors are utilized as a basis. The novel characteristic of these equations is that they form a set of second order ordinary differential equations (ODEs) in both the coordinates and the Lagrange multipliers associated to the constraint action. Based on these equations, the corresponding weak form is first obtained, leading eventually to a consistent first order ODE form of the equations of motion. These equations are found to appear in a form resembling the form obtained after application of the classical Hamilton׳s canonical equations. Finally, the new theoretical findings are illustrated by three representative examples.     

Lie group integration for constrained generalized Hamiltonian system with dissipation by projection method

For the constrained generalized Hamiltonian system with dissipation, by introducing Lagrange multiplier and using projection technique, the Lie group integration method was presented, which can preserve the inherent structure of dynamic system and the constraintinvariant. Firstly, the constrained generalized Hamiltonian system with dissipative was converted to the non-constraint generalized Hamiltonian system, then Lie group integration algorithm for the non-constraint generalized Hamiltonian system was discussed, finally the projection method for generalized Hamiltonian system with constraint was given. It is found that the constraint invariant is ensured by projection technique, and after introducing Lagrange multiplier the Lie group character of the dynamic system can‘ t be destroyed while projecting to the constraint manifold. The discussion is restricted to the case of holonomic constraint. A presented numerical example shows the effectiveness of the method.     

柔性机械臂碰撞系统的一种建模方法及精细积分数值模拟

对一单连杆柔性机械臂与外界发生碰撞时的动力学特性进行了研究。在建立动力学模型时，利用有限元思想，用有限个刚性单元和抗扭连接弹簧来模拟机械臂及其柔性性质，得出了较为简单的非线性动力学方程，根据几何关系，建立了系统在发生碰撞时的约束条件。本文在系统的动力学方程中引入对偶变量，并将方程导向Hamilton体系，用精细积分法对系统的动力学特性进行了模拟。     

GAUSS PRINCIPLE OF LEAST CONSTRAINT IN GENERALIZED COORDINATES AND ITS GENERALIZATION

当采用广义坐标描述系统的运动时,相比质点形式的高斯最小拘束原理,通过广义坐标形式的高斯最小拘束原理来建立动力学优化模型,计算效率更高. 从高斯原理的变分形式出发推导了广义坐标形式的高斯最小拘束原理,并研究了非理想约束、单边约束及刚体碰撞情形下的高斯最小拘束原理的形式. 研究认为:对刚体碰撞情形下,高斯最小拘束原理不能取代碰撞恢复定律,碰撞恢复定律以碰撞后广义速度的约束方程形式起作用.     

Exact solutions to magnetogasdynamics using Lie point symmetries

In the present work, we find some exact solutions to the first order quasilinear hyperbolic system of partial differential equations (PDEs), governing the one dimensional unsteady flow of inviscid and perfectly conducting compressible fluid, subjected to a transverse magnetic field. For this, Lie group analysis is used to identify a finite number of generators that leave the given system of PDEs invariant. Out of these generators, two commuting generators are constructed involving some arbitrary constants. With the help of canonical variables associated with these two generators, the assigned system of PDEs is reduced to an autonomous system whose simple solutions provide nontrivial solutions of the original system. Using this exact solution, we discuss the evolutionary behavior of weak discontinuities.     

Forward and inverse dynamics of nonholonomic mechanical systems

This paper deals with the forward and the inverse dynamic problems of mechanical systems subjected to nonholonomic constraints. The intrinsically dual nature of these two problems is identified and utilised to develop a systematic approach to formulate and solve them according to an unified framework. The proposed methodology is based on the fundamental equations of constrained motion which derive from Gauss’s principle of least constraint. The main advantage arising from using the fundamental equations of constrained motion is that they represent an effective method capable to derive the generalised acceleration of a mechanical system, constrained in general by a set of nonholonomic constraints, together with the generalized constraint forces (forward dynamics). When the constraint equations are used to represent the desired behaviour of the mechanical system under study, the generalised constraint forces deriving from the fundamental equations of constrained motion provide the control actions which reproduce the specified motion for the system (inverse dynamics). This approach is systematically extended to underactuated mechanical systems introducing a new method named underactuation equivalence principle. The underactuation equivalence principle is founded on the key idea that the underactuation property of a mechanical system can be mathematically represented using a particular set of nonholonomic constraint equations. Two simple case-studies are reported to exemplify the proposed methodology. In the first case-study the computation of the generalised constraint forces relative to the revolute joint constraints of a physical pendulum is illustrated. In the second case-study the calculation of the control action which solves the swing-up problem for an inverted pendulum is described.