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

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

Projected Runge-Kutta methods for constrained Hamiltonian systems

Projected Runge-Kutta (R-K) methods for constrained Hamiltonian systems are proposed. Dynamic equations of the systems, which are index-3 differential-algebraic equations (DAEs) in the Heisenberg form, are established under the framework of Lagrangian multipliers. R-K methods combined with the technique of projections are then used to solve the DAEs. The basic idea of projections is to eliminate the constraint violations at the position, velocity, and acceleration levels, and to preserve the total energy of constrained Hamiltonian systems by correcting variables of the position, velocity, acceleration, and energy. Numerical results confirm the validity and show the high precision of the proposed method in preserving three levels of constraints and total energy compared with results reported in the literature.     

Lyapunov stable penalty methods for imposing holonomic constraints in multibody system dynamics

This paper presents stability and convergence results on a novel approach for imposing holonomic constraints for a class of multibody system dynamics. As opposed to some recent techniques that employ a penalty functional to approximate the Lagrange multipliers, the method herein defines a penalized dynamical system using penalty-augmented kinetic and potential energies, as well as a penalty dependent constraint violation dissipation function. In as much as the governing equations are not typically cocreive, the usual convergence criteria for linear variational boundary value problems are not directly applicable. Still numerical simulations by various researchers suggest that the method is convergent and stable. Despite the fact that the governing equations are nonlinear, the theoretical convergence of the formulation is guaranteed if the multibody system is natural and conservative. Likewise, stability and asymptotic stability results for the penalty formulation are derived from well-known stability results available from classical mechanics. Unfortunately, the convergence theorem is not directly applicable to dissipative multibody systems, such as those encountered in control applications. However, it is shown that the approximate solutions of a typical dissipative system converge to a nearby collection of trajectories that can be characterized precisely using a Lyapunov/Invariance Principle analysis. In short, the approach has many advantages as an alternative to other computational techniques:

Stability and vibration of a helical rod constrained by a cylinder

The stability and vibration of an elastic rod with a circular cross section under the constraint of a cylinder is discussed. The differential equations of dynamics of the constrained rod are established with Euler’s angles as variables describing the attitude of the cross section. The existence conditions of helical equilibrium under constraint are discussed as a special configuration of the rod. The stability of the helical equilibrium is discussed in the realms of statics and dynamics, respectively. Necessary conditions for the stability of helical rod are derived in space domain and time domain, and the difference and relationship between Lyapunov’s and Euler’s stability concepts are discussed. The free frequency of flexural vibration of the helical rod with cylinder constraint is obtained in analytical form. The project supported by the National Natural Science Foundation of China (10472067). The English text was polished by Yunming Chen.     

多体动力学程序切断铰的处理方法

多体系统动力学从非树系统派生出树系统的计算需要进行切断铰的处理。切断铰约束方程的形成是进行多体系统程序编写时的重要部分,其处理过程复杂,需要一定的技巧。本文引入了约束正交补轴的概念,详细介绍了几种典型(旋转铰、万向节、棱柱铰、旋转棱柱组合铰)切断铰位移约束方程、速度约束方程、加速度约束方程的形成方法,并给出了详细的程序化过程,该方法适用于任何类型的切断铰。最后给出相应算例,结果表明本文的方法能快速、正确地形成切断铰约束方程。     

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

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

具有非线性homologous变形约束的桁架结构形态分析

研究了具有非线性homologous变形约束条件的桁架结构形态分析问题。在已有的线性homologous变形约束桁架形态分析的基础上,将结构的节点分成三类:homologous变形约束节点,形状可变节点和边界点。运用Moore-Penrose广义逆矩阵性质,将基础方程组解的存在条件表示为包含形状可变节点未知坐标的非线性方程组,为采用Newton-Raphson方法求解非线性方程组,对AA (A为任意矩阵,A 为A的Moore-Penrose广义逆矩阵)求偏导数,找到了满足保型要求的形态,给出的桁架算例说明了本文方法的有效性。     

关于一类电-力振动模型的建模与分析

从电学、力学的基本原理出发,通过数学方法建立和分析了一种电-力振动模型。这是一种较复杂的电力混合作用的线性振动系统,由模型的特殊结构(类扬声器结构)作者确立了两种磁场,即感生磁场和外磁场两者正交独立,并规定了电学、力学两种不同物理量的坐标取向关系。该模型需要求得三阶正系数常微分方程的收敛解,再求得包含暂态、稳态项的完整解。另外,文章从能量和做功的角度,通过对电压电流间的相位差分析,对所建模型的正确性作了论证,同时也为这类建模引荐了一种论证手段。     

基体开裂的正交复合材料对称层合板的粘弹性力学行为

本文讨论了含有均匀基体裂纹的正交复合材料对称层合板的线性粘弹性力学行为.采用二维剪切滞后模型并对其层间剪应力在厚度方向进行线性假设分布,求得层合板的平均应力应变的线弹性解,利用等效约束模型和经典层合理论可得到层合板因为含有基体裂纹而所引起的刚度退化现象.在弹性-粘弹性对应原理的基础上对其层合板的线粘弹性进行了讨论和研究.结果表明:层合板的松弛模量和蠕变泊松比随着时间的增加而减少,到达稳态后其值基本上是恒值.并跟Zocher的解析解和有限元数值解作了比较,发现结果非常吻合.     

板材冲压成形有限元数值模拟界面摩擦约束处理方法

基于弹塑性大变形板材冲压成形增量有限元法,采用线性弹簧单元提出一种空间三维板材冲压成形等效界面摩擦力的处理方法.这种方法的优点是可以反映界面摩擦力的同步、被动产生效果,同时它还可以保证约束后刚度矩阵的对称性.将提出的摩擦约束处理算法引入自主开发的板材成形模拟软件FASTAMP,计算了汽车油底壳的成形过程,计算结果的厚度分布与实际冲压件吻合比较好,验证了等效界面摩擦力约束处理方法的有效性.