可展桁架展开过程模拟的违约校正

在 Yoon等学者的基础上 ,利用泰勒级数和 Moore-Penrose广义逆对主动校正法进行了深入、清晰的阐述 ,并针对可展桁架结构展开模拟所遇到的完整定常约束 ,发展了一种简单实用且精度较高的能量和速度违约校正方法。文中给出的算例说明了本文的违约校正算法的有效性     

多体系统指标2运动方程HHT方法违约校正1）

采用Cartesian绝对坐标建模方法,完整约束多体系统运动方程是指标3的微分--代数方程（differentialalgebraic equations,DAEs）,数值求解指标3的DAEs属于高指标问题,通过对位置约束方程求导,可使运动方程的指标降为2.位置约束方程求导得到的是速度约束方程.直接求解指标3的运动方程,速度约束方程得不到满足,而且高指标DAEs的数值求解存在一些问题.论文首先采用HHT（Hilber--Hughes--Taylor）直接积分方法求解降指标得到的指标2运动方程,此时速度约束方程参与离散计算,从机器精度上讲速度约束自然得到满足,而位置约束方程没有参与计算,存在“违约”.针对违约问题,采用基于Moore--Penrose广义逆理论的违约校正方法,消除位置约束方程的违约.指标2运动方程HHT方法违约校正,将HHT方法和违约校正方法很好地结合,在数值求解指标2运动方程的过程中,位置约束方程和速度约束方程都不存在违约问题,而且新方法没有引入新的未知数向量,离散得到的非线性方程组的方程数量与原指标2运动方程的方程数量相同,求解规模没有扩大.新方法的实用和有效性通过算例的数值实验得到验证,数值实验也说明新方法保持了HHT方法本身具有的数值阻尼可以控制和二阶精度的特性.最后从非线性方程组的求解规模和计算速度上与其他方法进行了比较分析,说明新方法的优势所在.     

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

多体系统动力学方程为3阶微分代数方程,已有的约束违约稳定法存在位移违约问题,数值仿真准确性和稳定性不足。本文将求解高阶微分代数方程的降阶理论、ε嵌入处理方式与隐式龙格库塔法相结合,提出了直接满足位移约束条件的多体系统动力学方程的无违约算法,避免了约束违约问题。该方法先将多体动力学方程转化为2阶微分代数方程,并与位移约束方程联立;再应用ε嵌入隐式龙格库塔法进行数值求解。应用两种方法分别对单摆机构进行数值仿真,结果表明本文的方法不仅能适应较大步长,且准确性和稳定性均优于约束违约稳定法。     

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

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

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

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

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

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

黎曼位形空间中约束多体系统的动力学分析

在黎曼位形空间中研究了约束多体系统的动力学问题．通过对约束矩阵的奇异值分解和修正的Ｇｒａｍ Ｓｃｈｍｉｄｔ过程构造了系统流形的法向和切向子空间的正交归一化基，将系统的动力学方程沿这双基进行投影，得到了求系统动力学响应的新型公式，给出了其数值分析的一种方法，并举了算例．     

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

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

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

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

含非理想约束多柔体系统递推建模方法

基于多体系统中邻接物体运动学递推关系,可以证明树状多体系统中末端物体的作用体现为传递给其内接物体的惯性和外力. 由于闭环系统切断铰约束反力和非理想约束反力可看作为系统外力,任何复杂系统都可以转化为等效的树系统,并且系统约束方程中所涉及的广义加速度可以系统化地用描述约束反力的拉氏乘子替换. 基于以上结果,提出了针对含非理想约束多柔体系统递推建模方法. 利用该方法可以将复杂多体系统动态减缩为单个物体,从而在求解系统加速度时不需对整个系统的质量矩阵进行求逆运算,同时大幅度地降低了非理想约束反力方程的维数. 通过一个算例具体说明了所提方法的求解过程,算例结果与现有商业软件所得结果一致.      

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

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

一种建立多体系统计算模型的新方法

本文根据约束的数学力学性质提出了一种新的多体系统建模方法。同传统的缩并法相比，本文所提出的方法具有更好的数值稳定性和较高的计算效率。     

柔性多体系统的计算策略

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

多刚体系统动力学的子系统递推法

本文采用递归推导并定义铰基矢阵得到以相对速度或伪速度表示的动力学方程。对非树系统利用回路转换矩阵及铰基矢阵自动建立约束方程,并用奇异值分解缩聚法得到以独立伪速度为变量的动力学方程。并导出递推计算公式便于程序执行。相应的通用程序SNERM已产生。     

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时间进行了比较分析．     

含摩擦与碰撞平面多刚体系统动力学线性互补算法

基于接触力学理论和线性互补问题的算法, 给出了一种含接触、碰撞以及库伦干摩擦, 同时具有理想定常约束(铰链约束) 和非定常约束(驱动约束) 的平面多刚体系统动力学的建模与数值计算方法. 将系统中的每个物体视为刚体, 但考虑物体接触点的局部变形, 将物体间的法向接触力表示成嵌入量与嵌入速度的非线性函数,其切向摩擦力采用库伦干摩擦模型. 利用摩擦余量和接触点的切向加速度等概念, 给出了摩擦定律的互补关系式; 并利用事件驱动法, 将接触点的黏滞-滑移状态切换的判断及黏滞状态下摩擦力的计算问题转化成线性互补问题的求解. 利用第一类拉格朗日方程和鲍姆加藤约束稳定化方法建立了系统的动力学方程, 由此可降低约束的漂移, 并可求解该系统的运动、法向接触力和切向摩擦力, 还可以求解理想铰链约束力和驱动约束力. 最后以一个类似夯机的平面多刚体系统为例, 分析了其动力学特性, 并说明了相关算法的有效性.      

有限转动张量的保角参量及在多柔体系统中的应用

采用保角转动参数描述了多体系统中的大转动张量．该方法消除了传统的欧拉参数描述所必需的约束方程，并且适于大变形部件的建模需要．利用以上结果建立了含大变形梁状部件的多体系统的力学模型．