多体系统动力学微分／代数方程组数值方法

多体系统动力学微分／代数混合方程组又称Ｅｕｌｅｒ－ｌａｇｒａｎｇｅ方程,是近十年来动力学和计算数学领域研究的热点之一．本文介绍这两个领域中引入的传统的数值积分方法与新的理论．      

多体系统Lagrange方程数值算法的研究进展

Lagrange方法是建立多体系统动力学方程的普遍方法之一, 其方程的形式为常微分方程组或微分-代数方程组,数值计算与数 值分析是研究多体系统动力学特性的重要方法。本文简要介绍了多 体系统动 力学方程的第一、二类Lagrange方程和修正的Lagrange方 程的基本形式及这些方程的正则形式,着重介绍了正则方程在数值 计算中的特点,就多体系统Lagrange方程的隐式算法、辛算法和多 体系统动力学特性的数值分析方法(包括数值仿真、 Poincarè映射 和Lyapunov指数的计算方法)的研究现状进行了综述。     

多体系统Lagrange方程数值算法的研究进展

Lagrange方法是建立多体系统动力学方程的普遍方法之一,其方程的形式为常微分方程组或微分 - 代数方程组,数值计算与数值分析是研究多体系统动力学特性的重要方法.本文简要介绍了多体系统动力学方程的第一、二类Lagrange方程和修正的Lagrange方程的基本形式及这些方程的正则形式,着重介绍了正则方程在数值计算中的特点,就多体系统Lagrange方程的隐式算法、辛算法和多体系统动力学特性的数值分析方法(包括数值仿真、Poincar'e映射和Lyapunov指数的计算方法)的研究现状进行了综述.     

微分几何与向后差分相结合的多柔体系统动力学数值方法

将处理微分－代数混合方程的微分几何方法与向后差分相结合，建立了一种新的求解多柔体系统的力学方程的数值方法，通过典型算例验证了方法的正确性和有效性。     

多体系统动力学方程符号推导

符号演算又名计算机代数,是近二十年来蓬勃发展的计算机科学。用符号演算导出复杂的多体系统动力学非线性微分方程组,可以提高计算效率,计算精度,减轻人的劳动並为设计方案的选择提供快速分析的依据。本文介绍用FORTRAN-77在VAX/11—780上符号推导以6自由度机器人为例的多体系统动力学方程。文中方法可应用于航天器、机构等动力学方程推导。     

树形多体系统非线性动力学的数值分析方法

研究了树形多体系统大线性动力学分析的数值方法,利用多体系统的正则方程及其线性化程,给出了多体系统Ｌｙａｐｕｎｏｖ指数和Ｐｏｉｎｃａｒｅ映射的计算方法,该算法具有较好的计算精度和通用性,既适用于说明该算法的有效性,并对该系统的动力学行为进行分析,最后用算例说明该算法的有效性,并对该系统的动力学特征（周期解、准周期解、分岔、混沌以及通往混沌的道路等）进行了分析。     

含随机参数的多体系统动力学分析

基于Lagrange方程建立了含随机参数的多体系统的动力学 模型,利用广义坐标分离法将随机微分代数方程转化为随机纯微分方程,利用Newmark法进行数值解算. 应用随机因子法求解系 统随机响应的数字特征,获得统计意义下的解. 以旋转杆滑块系统为例,考虑系统中载荷、物理和几何参数的随机性,通过与Monte Carlo法结果的对比验证了文中方法的正确性和有效性. 计算结果表明,部分随机参数的分散性对多体系统动力响应的影响不可忽略,利用随机参数的动力学模型将能客观地反映出系统的动力学行为.     

动力刚化多体系统动力学

本文利用几何非线性的应变－－位移关系，在小变形假设条件下，得到了一般柔性构件弹性有的广义价值标二阶小量表达式。在此基础上，利用Ｋａｎｅ方程的Ｈｕｓｔｏｎ方法，在推导偏（角）速度表达式后，作适当的线性化处理，以使动力刚度项得以保留，从而建立了动力刚化多体系统的动力学方程，仿真算例证明了该理论的正确性和有效性。     

多柔体系统动力学方程一种数值求解方法

本文提出用含有二阶自校正功能的增广法和半隐式Runge-Kutta法相结合求解多柔体系统动力学中的微分-代数型方程。并给出了计算实例和实验验证。     

柔性多体系统动力学通用算法研究

柔性多本系统运动学动力学仿真通用软件可广泛服务于工程领域,基于Ｋａｎｅ方程建立一种开发该类软件的一般算法,通过定义的辨识函数解决束识别问题,利用振型描述构件的弹性变形,建立了递推格式的运动学模型,提出了计算偏带度、偏角速度和运动微分方程系数矩阵的“０－１”法,最后给出的算例表明所建立的算法是可行的。     

NEW GENERALIZED-a METHOD FOR OVER-DETERMINED MOTION EQUATIONS IN MULTIBODY DYNAMICS

完整约束多体系统第一类Lagrange方程建模得到的运动方程是指标-3形式的微分-代数方程（differental-algebraicequations,DAEs）．如果同时考虑速度约束,将得到超定运动方程,该方程是指标-2的超定微分-代数方程（over．determineddifferential-algebraicequations,ODAEsl．基于结构动力学中常用的广义-OZ方法,将其拓展,求解包含速度约束的超定运动方程,相对于其他求解指标-2ODAEs的算法,新的算法没有增加离散得到的非线性方程组方程的数目．通过数值实验验证算法,并说明其求解ODAEs不存在精度降阶的现象,仍然具有二阶精度,同时算法的数值耗散也是可以控制的．最后新方法与其他求解多体系统0DAEs形式运动方程算法的CPU时间进行了比较分析．     

多体系统Euler-Lagrange方程的最小二乘法与违约修正

针对受完整约束的多体系统动力学Euler-Lagrange方程,在其传统的直接增广算法和零空间方法基础上提出了当系统存在冗余约束情形下的最小二乘法.同时,对应于最小二乘法提出了改进的约束违约修正方法.本文还针对Euler-Lagrange方程的计算过程给出了相应Jacobi矩阵的QR分解和零空间连续正交基的算法.最后,以平行五连杆机构给出了数值结果并与部分现有方法进行比较.     

柔性多体系统动力学研究及存在的问题

本文对柔性多体系统动力学研究领域中动力学模型的建立、求解以及该领域中存在的问题进行了简要综述。     

非光滑多体系统动力学数值算法的研究进展

非光滑多体系统动力学数值计算方法是多体系统动力学研究的重要内容之一. 本文介绍了近年来含摩擦与碰撞的非光滑多体系统动力学数值算法方面的研究进展. 首先, 讨论了库仑摩擦模型和修正的库仑摩擦模型, 以及具有单边和双边约束的多体系统中法向约束力的特点. 其次, 回顾了基于连续模型和非连续模型的多体系统动力学方程的数值计算方法, 详细介绍了基于互补概念的非光滑多体系统动力学的事件驱动法和时间步进法, 分析比较了相关的数值算法. 最后, 指出了一些需要进一步研究的问题.     

柔性多体系统的计算策略

对柔性多体系统计算建模的研究现状和近期进展进行了总结. 重点讨论了柔性多体 动力学的以下内容: 柔性构件的建模, 约束建模, 求解技术, 控制策略, 耦合问 题, 设计和实验的研究. 对柔性多体系统建模的浮动坐标系,转动坐标系和惯性系 等3种坐标系的特点进行了对比. 指出了未来的研究方向, 包括柔性多体系统的新 的应用,如微观力学系统和超微观力学系统等; 提高这些模型的计算精度和效率的 技巧和策略; 以及可以用于改善柔性多体系统的工具. 本综述文章引用了877篇参 考文献.     

多体动系统动力学方程在流形上的辛算法

多体系统动力学方程的数值方法一直是数学与力学家们的热门研 究课题.特别是多体系统动力学微分/代数方程组形式的数学模型，是所 谓的指标-3问题，它的求解是一个难题.目前流行的关于它的数值方法都 有不尽人意的地方，主要是对出现的所谓的违约问题和刚性问题未很好地 解决.多体系统动力学方程在流形上的辛算法是近几年出现的新的数值方 法，它将闭环型多体系统动力学的方程的约束部分和常微分方程部分利 用所谓的辛方法同时进行处理，其中的一些方法已证明是有效的，所以， 用它求解多体系统动力学方程前景看好.本文介绍了这些新的理论，并提 出了一些有待解决的问题.