多说一句：天宫一号上天，值此举国欢庆的时候，将“CAE自主创新发展战略” 技术科学论坛中发表的“机械系统CAE与计算多体系统动力学”登出来。“天宫一号”空间实验室的关键技术是交会对接，该文就是介绍这方面研究成果的。借花献佛，以飨同行。
刚在电视上见到洪嘉振教授上台当选为上海市先进的报道。28日在上海交通大学闵行校区见到洪嘉振教授，听说又承接了个国家项目，再接再厉取得更大成绩，在此祝贺了！
本来我们之间也有个合作项目要搞，只是不幸夭折，想来有些可惜。
机械系统CAE与计算多体系统动力学
洪嘉振
（上海交通大学工程力学系）
摘要：本文首先简要介绍机械系统动力学仿真软件与计算多体系统动力学的关系与现状，重点介绍计算多体系统动力学的理论、计算方法与软件开发方面的主要成果。最后对我国机械系统动力学仿真软件系统的开发提出一些建议。
关键词：机械系统、动力学仿真、多体系统动力学、CAE软件。
一、 前言
众所周知，由于国内外制造业的工程设计与优化的需求，无论在功能还是用户界面的友好等方面，目前机械系统的CAE大型软件都已经达到相当高的水平。这些软件从功能上可以分为两大类：一类如ProE、UG等以计算机辅助设计(CAD)为目标的工程设计软件；另一类如ADAMS、DADS与RecuDyn等以系统静力学、运动学、动力学与控制性态计算机辅助分析(CAA)与优化(CAO)为目标的系统动力学仿真软件。
多少年来，我国许多高等院校的机械学科以计算机图学为基础，在以CAD为目标的工程设计软件方面取得成果斐然，开发了许多具有我国知识产权的CAD软件，且在工程中得到应用。
从提高工程设计的效率、减少投资风险与增强产品竞争力的角度，系统动力学仿真软件在机械系统CAE的地位日显重要。CAA与CAO软件的特点是将机械系统抽象为具有大范围相对运动物体构成的多体力学模型，利用力学基本原理建立能适用于任意拓扑的通用与程式化的运动学或动力学数学模型，通过数值计算得到系统的静力学、运动学、动力学与控制性态的时间历程，实现对系统设计的分析与优化。这种模式犹如结构分析软件。如同有限元理论在结构分析软件的地位，机械系统动力学仿真软件的理论基础是计算多体系统动力学。以机械系统的构件为刚体假设的多刚体系统动力学，在建模理论、计算方法以及软件实现方面已经非常成熟，能有效解决这类工程对象连续变化的问题。80年代后，这些软件逐渐为国内制造业认可，开始应用解决工程设计问题。
图1对接机构的接触碰撞
工程对象中存在大量碰撞与间隙等复杂的不连续的动力学过程，即在运动过程中系统的自由度会发生突变。例载人航天的关键技术-异体同构对接机构的动力学过程就是有无限多种可能的接触碰撞的工况(见图1)。这种变拓扑的动力学系统仿真技术对于动力学仿真软件必须给予考虑。它将涉及碰撞与间隙等动力学模型的研究、接触点的判断、系统拓扑切换条件实现等等的问题。目前的动力学仿真软件距能稳定、有效地解决不连续的动力学过程尚有一定的距离[1]。
大型、轻质、高速是当前机械系统的特征，构件必须做柔性体假设才能描述系统刚-柔耦合的动力学性态。这类系统的力学模型为柔性多体系统。带柔性体的复杂系统动力学是涉及一般力学，固体力学，计算方法与软件工程多学科相互交叉的新学科。如何建立能精确描述刚-柔耦合的动力学性态的数学模型还有待深入研究，系统动力学模型的数值病态成为动力学仿真的瓶颈[1]。
多体系统动力学的数学模型是微分-代数方程组，其数值求解的方法远不如纯微分方程成熟。探索此类方程可靠稳定的数值方法，仍然有大量问题有待解决[1]。
本文重点介绍能成功解决柔性多体系统动力学问题的单向递推组集的建模方法和柔性多体系统刚柔耦合动力学建模理论的进展。文章最后对我国机械系统动力学仿真软件系统的开发提出一些建议。
二、 计算多体系统动力学的进展
1. 单向递推组集的建模方法[2]
国外的大型动力学仿真软件ADAMS与DADS称它们能解决柔性多体系统动力学问题。采用的方法是将应用于多刚体系统的笛卡儿方法进行拓展，即仍然是最大坐标数的建模方法。其数学模型是微分-代数方程组：
，
(1)
其中p为低频大范围笛卡儿坐标，a为高频微小的变形坐标，两种变量的高度耦合，加上微分-代数方程组数值方法不成熟，很难保证数值的稳定与可靠。
不同推导方法得到的动力学方程尽管形式不同，但是理论上是等价的。然而其数值求解的性态是不同的。基于这样的思想，如果采用多体系统动力学最小坐标数的建模方法，得到的方程为纯微分方程：
，
(2)
其中p为低频大范围拉格朗日坐标，a为高频微小的变形坐标。从数值方法成熟程度显然方程(2)优于方程(1)。但是方程(2)的质量阵A与力阵B中两种坐标变量的耦合程度将更严重，计算复杂。为了解决此矛盾，笔者采取如下的计算策略：(1) 同时采用描述系统的笛卡儿坐标与拉格朗日坐标，质量阵A与力阵B的计算过程实现根据系统的拓扑向外递推，然后进行组集；(2) 通过推导引入可作前处理的不变量，尽可能降低两种坐标变量的耦合。通过上述处理大大地减少了计算量，提高了计算的精度。
这种柔性多体系统的单向递推组集的建模方法不追求方程形式上的优美，而从数值计算的优化上实现动力学模型算法上的构成，故具有较好的数值性态。根据此方法开发的具有自主知识产权的软件系统CADAMB已经运行了十余年。基本上排除了动力学仿真数值病态的困难。在我国多种卫星的总体设计中发挥了作用。例如对我国某对地探测卫星，利用该仿真软件全面深入研究了的大型帆板的展开和撞击动力学过程、帆板展开和撞击对卫星姿态的影响等等。该项研究通过评审，评审意见认为软件系统的数学模型准确、分析结果合理可信，与地面试验结果的趋势吻合，为该卫星的研制提供了重要的依据。
此外，当构件的柔性不考虑时，该建模方法也可应用于多刚体系统的动力学的仿真，由于它的递推组集的特征，仿真计算的效率大大地提高。
2. 柔性多体系统刚柔耦合动力学建模理论
近30多年来，柔性多体系统刚柔耦合动力学建模理论的前期研究大致分为如下三个阶段[1]：
(1) 运动-弹性动力学(KED)方法。该方法首先对系统构件作刚体假设，利用多刚体系统动力学方法进行构件运动分析。然后通过施加惯性力，对构件进行弹性变形和强度分析。这种方法不计构件大范围运动与弹性变形运动的耦合。由于方法简单，对于速度不高的机械系统仿真仍然得到应用。
(2) 混合坐标方法。该方法提出用描述构件大范围运动浮动系的刚体坐标与描述柔性体变形运动的有限元的节点坐标(或模态坐标)建立柔性多体动力学离散的数学模型。这种方法考虑了构件弹性变形与大范围运动的相互耦合。由于对柔性体在线性的范围离散时没有考虑大范围运动对其的影响，得到的柔性多体系统耦合动力学方程的精度实质上处在零次近似范围上。在解决实际问题中长期困惑此方法的问题是上述的数值求解的病态。
(3) 动力刚化问题的研究。1987年Kane发现上述零次近似耦合模型在处理高速旋转的悬臂梁动力学时会得到错误结论。以后的十多年，国内外的学者对该对象，采用不同的假定，引入所谓的“动力刚度项”对零次近似模型加以修正。但是在如何反映物体大范围运动与弹性变形的相互耦合，所采用的假定各不相同[1] 。上述的研究工作说明，柔性多体系统动力学的建模理论有待解决的实质性问题是需揭示刚柔耦合的机理，从力学基本原理建立描述刚柔耦合动力学较精确的动力学模型，应通过实验验证该模型的正确性。
在国家自然科学基金重点项目和教育部博士点基金项目的资助下，笔者的课题组在刚柔耦合动力学建模理论研究过程中，基于这样的新思想[3]：(1)建立较精确的刚柔耦合动力学模型主要目的为有效解决刚-柔耦合多体系统动力学的问题，Kane的非惯性系下的结构动力学问题是其的特殊情况；(2)所谓刚柔耦合“动力刚度项”是客观存在的，不应该通过引入新的假定来“捕捉”。而且刚柔耦合不只影响“刚度”，而且将会对“阻尼”等其他动力学特征产生影响。根据前期研究工作分析认为，造成零次近似耦合动力学方程缺陷的主要原因是应用没有大范围运动的线性变形场理论结果，人为地丢失了其二阶项与大范围运动耦合，而这种耦合仍可达到一阶量级。当与其他一阶项可比时，将会严重影响动力学仿真结果。以中心刚体-悬臂梁系统为对象，以上述建模理论建立刚柔耦合动力学方程[3]：
，
(3)
其中， q为中心刚体定轴转动的转角，a为悬臂梁的变形坐标。该方程的广义质量阵M和广义力阵F中计及了大范围运动和高阶变形运动的耦合项。含这种耦合项的刚柔耦合动力学方程的精度在一次近似的量级上。当这些项不考虑时，方程退化为原零次近似刚柔耦合动力学方程。
定义q为时间的已知函数，方程(3)变为非惯性系下悬臂梁的动力学方程
(4)
在上式的K中可导出没有附加任何条件而理论依据严格的“动力刚度项”。 解释了非惯性系下的悬臂梁产生动力刚化的根本原因和动力学实质。找出了零次模型在高速大范围运动下得到失稳发散错误结论的原因。然而，这是零次近似模型得到错误结果的极端情况；在大范围运动为低速的情况，由于附加耦合项的作用，有时零次近似模型的计算的误差也相当可观。因此需要讨论大范围运动在什么条件下必须应用一次近似耦合模型才能得到正确的结论的条件。
由方程(3)出发对大范围运动为自由的刚体-悬臂梁两体系统进行数值研究，首次发现在一定的初始条件下，传统零次近似模型也会得到系统运动发散错误结果的反例，而利用一次近似耦合模型的计算能得到符合实际情况的结果。此项成果的意义还在于对长期困扰柔性多体系统动力学的数值发散问题做出了新的解释，即计算结果发散不一定是数值病态方面的原因，而应该注意力学模型是否正确与失效。
图2刚柔耦合系统实验
利用自行设计的实验方案，通过单轴气浮台对一柔性梁进行刚-柔耦合动力学试验(见图2)。将实验数据与一次近似模型和零次近似模型的理论结果进行对照[4]，观测到一次近似模型与实验结果吻合相当好。而对于零次近似模型只有大范围运动角速度在一定范围内或在比较小的初始条件范围内才有比较好的对照结果。这项研究从实验的角度验证了一次近似模型在处理非惯性系下的结构动力学问题与多体刚柔耦合动力学问题的正确性和有效性。
三、 国内机械系统动力学仿真软件开发的建议
国内CAE的发展并不平衡。在CAD自主软件的开发上，我国许多院校与研究单位作出巨大的努力，取得了可喜的成果，为我国制造业的工程设计自动化作出了巨大贡献。相对而言，在CAA自主软件的开发上有不小差距。在大型的通用软件方面，上海交通大学的CADAMB的功能已具一定的规模，在数学模型、计算方法与求解器方面具有一定的特点，在工程应用特别在航天器系统姿态耦合动力学的仿真与总体设计中已经发挥重要作用。
尽管如此，国内CAE的发展还存在不少问题：(1)与我国 CAD软件相比， CAA软件的用户界面，前后处理与国外软件有较大的差距。(2)CAD与CAA软件之间没有做到紧密的联系与无缝连接。(3)在CAE软件的产业化方面与国外软件相比问题不少。
据此，特对我国CAE的自主创新的发展建议如下：
(1) 加强我国CAD与CAA软件开发团队的强强联合，实现多学科人员的结合。在已经成熟成果的基础上，实现CAD与CAA软件之间的无缝连接。为我国工程设计提供若干个与国外软件可以相左的自主的CAE软件。
(2) 针对我国国情，利用当前的成果，在行业的层面上提供有特色的CAE软件，为行业新产品的开发提供有效的工具，提高他们在行业中的竞争力
(3) 对于CAA软件，仍需加强支撑软件的计算多体系统动力学的建模理论、计算方法的深入研究。从数值仿真与软件开发的角度，研究具有数值性态好的刚柔耦合的动力学数学模型，仍然是研究的核心。继续发挥国内在这方面的优势，让软件的可靠性与其他性能上具有国际竞争力。
(4) 当前的大量的工程对象的设计必须考虑机械系统与电系统的耦合，机械系统与其他控制系统的耦合等等。如笔者在参与的低压短路器仿真软件系统的开发就是典型的机-电系统的耦合。目前在多物理场工程问题的数值仿真中，尚处在利用各物理场模型的简单迭代层面上，其精度与收敛性方面均存在不少问题。因此当前需注意加紧多物理场工程问题CAE的理论、方法与软件工程的建设。从CAE自主创新角度，这方面的建设将具有国际竞争力。
发表于：
2011-09-30 11:03 阅读(1864)
评论(2) 收藏 好文推荐本博客所有内容，若无特殊声明，皆为博主原创作品，未经博主授权，任何人不得复制、转载、摘编等任何方式进行使用和传播。
作者该类其他博文：