Singularity-free augmented Lagrangian algorithms for constrained multibody dynamics

After a general review of the methods currently available for the dynamics of constrained multibody systems in the context of numerical efficiency and ability to solve the differential equations of motion in singular positions, we examine the acceleration based augmented Lagrangian formulations, and propose a new one for holonomic and non-holonomic systems that is based on the canonical equations of Hamilton. This new one proves to be more stable and accurate that the acceleration based counterpart under repetitive singular positions. The proposed algorithms are numerically efficient, can use standard conditionally stable numerical integrators and do not fail in singular positions, as the classical formulations do. The reason for the numerical efficiency and better behavior under singularities relies on the fact that the leading matrix of the resultant system of ODEs is sparse, symmetric, positive definite, and its rank is independent of that of the Jacobian of the constraint equations. The latter fact makes the proposed method particularly suitable for singular configurations.     

约束多体系统的基于离散零空间的隐式龙格库塔法

将离散零空间理论应用于多体系统动力学方程的数值计算,可降低多体系统动力学方程的维数。通过给出离散零空间理论与IRK法相结合的一般数学框架,提出了多体系统动力学的基于离散零空间理论的IRK法。数值算例表明：该算法可获得较满意的数值结果,约束违约程度很小,三种积分算法算例的范数均在10^-16之内。     

约束多体系统动力学正则方程的约束变尺度方法

对约束多体系统动力学的正则方程建立了最优化问题模型,采用辛差分格式,提出相应的最优化问题的约束变尺度法,该方法有效地提高了计算稳定性,算例说明了方法的有效性。     

A new algorithm for solving differential/algebraic equations of multibody system dynamics

I.Intr0ductionTheequationsofmoti0nofconstrainedmultibodysystemdynamicscanbewrittenasfollows11lwherelistimeparameter*q6R'isgeneralizedcoordinatesvectorofsystem,M(q,t):RnXR-R"xnisgeneralizedmassn1atrixofsystem,k6RmisLagrangemultipliervector,fo(q,t):RnXR-R"(…     

Kronecker product and a new matrix form of Lagrangian equations with multipliers for constrained multibody systems

The automatic derivation of motion equations is an important problem of multibody system dynamics. Firstly, an overview of the matrix calculus related to Kronecker product of two matrices is presented. A new matrix form of Lagrangian equations with multipliers for constrained multibody systems is then developed to demonstrate the usefulness of Kronecker product of two matrices in the study of dynamics of multibody systems. Finally, the equations of motion of mechanisms are derived using the proposed matrix form of Lagrangian equations as application examples.     

带约束多体系统动力学方程的隐式算法

研究了带约束多体系统隐式算法,用子矩阵的形式推导出了多体系统正则方程的Ｊａｃｏｂｉ矩阵,它适用于多种隐式算法并给出了隐式Ｒｕｎｇｅ－Ｋｕｔｔａ算法,最后用一算例表明了隐式算法的计算效率和精度明显优于算法。     

A modified constraint force algorithm for flexible multibody dynamics with loop constraints

Nonlinear Dynamics - This paper presents a modified constraint force algorithm (m-CFA) for the dynamics of flexible multibody systems in arbitrary topologies. The m-CFA can efficiently calculate...     

Partition method and experimental validation for impact dynamics of flexible multibody system

The impact problem of a flexible multibody system is a non-smooth, high-transient, and strong-nonlinear dynamic process with variable boundary. How to model the contact/impact process accurately and efficiently is one of the main difficulties in many engineering applications. The numerical approaches being used widely in impact analysis are mainly from two fields: multibody system dynamics (MBS) and computational solid mechanics (CSM). Approaches based on MBS provide a more efficient yet less accurate analysis of the contact/impact problems, while approaches based on CSM are well suited for particularly high accuracy needs, yet require very high computational effort. To bridge the gap between accuracy and efficiency in the dynamic simulation of a flexible multibody system with contacts/impacts, a partition method is presented considering that the contact body is divided into two parts, an impact region and a non-impact region. The impact region is modeled using the finite element method to guarantee the local accuracy, while the non-impact region is modeled using the modal reduction approach to raise the global efficiency. A three-dimensional rod-plate impact experiment is designed and performed to validate the numerical results. The principle for how to partition the contact bodies is proposed: the maximum radius of the impact region can be estimated by an analytical method, and the modal truncation orders of the non-impact region can be estimated by the highest frequency of the signal measured. The simulation results using the presented method are in good agreement with the experimental results. It shows that this method is an effective formulation considering both accuracy and efficiency. Moreover, a more complicated multibody impact problem of a crank slider mechanism is investigated to strengthen this conclusion.     

Partition method for impact dynamics of flexible multibody systems based on contact constraint

The impact dynamics of a flexible multibody system is investigated. By using a partition method, the system is divided into two parts, the local impact region and the region away from the impact. The two parts are connected by specific boundary conditions, and the system after partition is equivalent to the original system. According to the rigid-flexible coupling dynamic theory of multibody system, system＇s rigid-flexible coupling dynamic equations without impact are derived. A local impulse method for establishing the initial impact conditions is proposed. It satisfies the compatibility con- ditions for contact constraints and the actual physical situation of the impact process of flexible bodies. Based on the contact constraint method, system＇s impact dynamic equa- tions are derived in a differential-algebraic form. The contact/separation criterion and the algorithm are given. An impact dynamic simulation is given. The results show that system＇s dynamic behaviors including the energy, the deformations, the displacements, and the impact force during the impact process change dramatically. The impact makes great effects on the global dynamics of the system during and after impact.     

A METHOD FOR SOLVING THE DYNAMICS OF MULTIBODY SYSTEMS WITH RHEONOMIC AND NONHOLONOMIC CONSTRAINTS

ＡＭＥＴＨＯＤＦＯＲＳＯＬＶＩＮＧＴＨＥＤＹＮＡＭＩＣＳＯＦＭＵＬＴＩＢＯＤＹＳＹＳＴＥＭＳＷＩＴＨＲＨＥＯＮＯＭＩＣＡＮＤＮＯＮＨＯＬＯＮＯＭＩＣＣＯＮＳＴＲＡＩＮＴＳ￥ＳｈｕｉＸｉａｏｐｉｎｇ（水小平）ＺｈａｎｇＹｏｎｇｆａ（张永发）（Ｄｅｐａｒ...     

A projective criterion to the coordinate partitioning method for multibody dynamics

Summary The paper discusses some developments in the coordinate partitioning method for the dynamic analysis of constrained/closed-loop multibody systems. First, the method is reformulated to a more compact form. Then, a simple and reliable projective criterion for choosing the best coordinates from the redundant ones is proposed, and some advantages are pointed out that may arise in the method by applying inverse kinematics algorithms. Finally, the problem of determination ofphysical reactions of constraints and closing conditions is discussed. A four-bar linkage mechanism serves for an illustration of some aspects of the paper.

---

Ein Projektionskriterium zur Koordinatenwahl bei Mehrkörpersystemen mit kinematischen Schleifen

Übersicht Bei der Simulation der Bewegunsgleichungen von Mehrkörpersystemen mit kinematischen Schleifen in Minimalform stellt sich die Frage nach der Wahl günstiger verallgemeinerter Koordinaten. Hierfür wird ein Projektionskriterium vorgeschlagen, welches die Trennung der systembeschreibenden redundanten Koordinaten in die verallgemeinerten und davon abhängige Koordinaten gestattet. Durch Anwendung eines Verfahrens zur Rückwärtstransformation der kinematischen Beschreibung lassen sich in diesen explizite Schließbedingungen formulieren, was sich bei der Simulation vorteilhaft auswirkt; diese kann ohne Verletzung der Schließbedingungen erfolgen. Es wird auch gezeigt, wie auch bei Verwendung verallgemeinerter Koordinatenphysikalisch interpretierbare Reaktionen ermittelt werden können. Als Anwendungsbeispiel dient ein ebenes Viergelenk.

---

---

On leave from the Department of Mechanics, Technical University of Radom, ul. Malczewskiego 29, 26-600 Radom, Poland     

Direct sensitivity analysis of multibody systems with holonomic and nonholonomic constraints via an index-3 augmented Lagrangian formulation with projections

Optimizing the dynamic response of mechanical systems is often a necessary step during the early stages of product development cycle. This is a complex problem that requires to carry out the sensitivity analysis of the system dynamics equations if gradient-based optimization tools are used. These dynamics equations are often expressed as a highly nonlinear system of ordinary differential equations or differential-algebraic equations, if a dependent set of generalized coordinates with its corresponding kinematic constraints is used to describe the motion. Two main techniques are currently available to perform the sensitivity analysis of a multibody system, namely the direct differentiation and the adjoint variable methods. In this paper, we derive the equations that correspond to the direct sensitivity analysis of the index-3 augmented Lagrangian formulation with velocity and acceleration projections. Mechanical systems with both holonomic and nonholonomic constraints are considered. The evaluation of the system sensitivities requires the solution of a tangent linear model that corresponds to the Newton–Raphson iterative solution of the dynamics at configuration level, plus two additional nonlinear systems of equations for the velocity and acceleration projections. The method was validated in the sensitivity analysis of a set of examples, including a five-bar linkage with spring elements, which had been used in the literature as benchmark problem for similar multibody dynamics formulations, a point-mass system subjected to nonholonomic constraints, and a full-scale vehicle model.     

Discrete time transfer matrix method for dynamics of multibody system with flexible beams moving in space

In this paper,by defining new state vectors and developing new transfer matrices of various elements moving in space,the discrete time transfer matrix method of multi-rigid-flexible-body system is expanded to study the dynamics of multibody system with flexible beams moving in space.Formulations and numerical example of a rigidflexible-body three pendulums system moving in space are given to validate the method.Using the new method to study the dynamics of multi-rigid-flexible-body system moving in space,the global dynamics equations of system are not needed,the orders of involved matrices of the system are very low and the computational speed is high,irrespective of the size of the system.The new method is simple,straightforward,practical,and provides a powerful tool for multi-rigid-flexible-body system dynamics.     

单轴旋转惯性导航系统的水平姿态角误差修正方法

在单轴旋转惯性导航系统中,旋转轴不正交误差对于水平方向上的姿态角会产生一定的影响。针对传统最小二乘拟合估计旋转轴不正交角存在拟合近似误差大、未考虑到数据的统计特性等缺陷,提出一种基于Kalman滤波的水平姿态角误差修正方法。在传统方案的基础上,增加以加速度计的输出作为观测量,并采用带有确定性控制项的Kalman滤波方法估计旋转轴不正交角,进而修正转轴不正交误差。仿真实验结果表明,转轴不正交角较大时,该方案将水平姿态角误差峰峰值从传统的最小二乘拟合修正方案的50'进一步降低至10'以内,精度提升了约80%。实际试验结果表明该方案相比于传统的最小二乘拟合修正方案,提升了水平姿态角精度。     

开放式多体系统动力学仿真算法软件研发(I)DAEs求解算法构架设计

基于开放式工程与科学计算集成化软件平台SiPESC,研发了用于多体系统动力学时程分析的一类通用求解算法构架。该构架的核心思想是算法与数据相分离,整个构架由五个基本类及子类组成。本文重点阐述基本类的抽象过程,利用插件技术设计求解器的构架,进一步应用该构架实现了Newmark方法,HHT(Hilber-Hughes-Taylor)方法,Generalized α方法,Bathe方法及祖冲之类Symplectic方法等微分-代数方程组(DAEs)求解器的开发。研究工作表明,本文所提出的DAEs求解算法构架对多体系统动力学的时程分析具有良好的开放性和通用性,可方便进行各种新的DAEs求解算法的动态扩展。     

A constrained particle dynamics for continuum-particle hybrid method in micro- and nano-fluidics

A hybrid method of continuum and particle dynamics is developed for micro- and nano-fluidics, where fluids are described by a molecular dynamics (MD) in one domain and by the Navier–Stokes (NS) equations in another domain. In order to ensure the continuity of momentum flux, the continuum and molecular dynamics in the overlap domain are coupled through a constrained particle dynamics. The constrained particle dynamics is constructed with a virtual damping force and a virtual added mass force. The sudden-start Couette flows with either non-slip or slip boundary condition are used to test the hybrid method. It is shown that the results obtained are quantitatively in agreement with the analytical solutions under the non-slip boundary conditions and the full MD simulations under the slip boundary conditions.The project supported by Chinese Academy of Sciences under the innovative project “Multi-scale modelling and simulation in complex system” (KJCX-SW-L08) and National Natural Science Foundation of China (10325211).     

多体系统动力学子系统求解算法

为了降低系统求解规模,提高复杂多体系统动力学求解效率,研究了多体动力学子系统求解算法.通过分析系统拓扑构型,基于雅克比矩阵分块特性,提出子系统综合算法.算法采用半递归法建立含开环、闭环的子系统运动方程,将子系统质量和外力等效作用在内接基体上,解耦系统运动方程组装和加速度求解.此外,针对松散耦合子系统,提出了予系统分治算法.算法针对子系统间的耦合特点,分解复杂多体系统模型,实现子系统的独立求解.最后以机构动力学数值求解实例,验证上述算法的正确性和求解效率.