开放式多体系统动力学仿真算法软件研发(Ⅱ)DAEs求解算法对比

研究求解微分-代数方程组(DAEs)的高效率、高精度和高稳定性数值积分方法一直是多体系统动力学领域的热点问题之一。本文将求解结构动力学方程的Bathe数值积分策略应用于DAEs的求解,并基于SiPESC平台开发了开放式多体系统动力学仿真算法软件,综合比较研究了Newmark法、HHT-I3法、Generalizedα方法、Bathe方法和祖冲之类Symplectic方法。通过复摆、刚-柔耦合双摆和对称陀螺三个数值算例研究了算法参数与数值阻尼的关系。数值实验表明,Newmark方法在特定参数下引入的数值阻尼通常不可控,HHT-I3方法、Generalizedα方法和Bathe方法通过选择特定步长和参数可引入可控的数值阻尼,祖冲之类Symplectic方法无数值阻尼。在求解真实高频和低频耦合问题以及高速旋转的陀螺问题时,采用祖冲之类Symplectic方法或者无耗散的Newmar方法能够对系统的高频成分进行准确模拟。     

考虑弹簧阻尼作动器解析雅可比矩阵的多刚体动力学分析

弹簧-阻尼-作动器(spring-damper-actuator,SDA)是多体系统中常见的力元,在工程领域中有着广泛的应用.采用绝对坐标方法建立的多体系统动力学控制方程通常是复杂的非线性微分-代数方程组.为了保证数值解的精度和稳定性,通常需要采用隐式算法求解动力学方程,而雅可比矩阵的计算在隐式数值求解过程中至关重要.对于含有SDA的多体系统,SDA造成的附加雅可比矩阵是与广义坐标和广义速度相关的高度非线性函数.目前的很多研究工作专注于广义力向量的计算,然而对附加雅克比矩阵的计算则少有关注.针对含SDA的多刚体系统进行动力学分析,首先基于Newmark算法研究其在动力学方程求解中的雅可比矩阵的构成形式;然后推导SDA的广义力向量对应的附加雅可比矩阵,其中包括广义力向量对广义坐标和对广义速度的偏导数矩阵.最后通过两个数值算例研究附加雅可比矩阵对动力学分析收敛性的影响;数值分析表明:当SDA的刚度、阻尼和作动力数值较大时,SDA导致的附加雅可比矩阵对数值解的收敛性有重要影响;当考虑SDA对应的附加雅可比矩阵时,动力学分析可以以较少的迭代步实现收敛,从而减少分析时间.     

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

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