柔体动力学介绍

《柔体动力学介绍》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、柔体动力学介绍一、KED（Kineto-Elastodynamics）法KED法，即运动弹性动力学，由美国学者Erdman和Sandor提出。该方法的研究始于上个世纪60年代，早期研究者仅把部件（一般是一个，如四杆机构的连杆）看作是柔性的，并且只考虑其一种变形（如杆件的弯曲变形），方程中也引入较多假设。70年代初期，Erdman和Sandor将结构动力学中的有限元方法移植到机构分析中来，克服了模型过于简单的缺陷。我国自80年代初开始研究机构弹性力学，学者张策对KED法做了大量研究。KED法在分析机构的真实运动时，均假设：Ø与采用刚性机构的运动分析法的到的机构名义运动的位移相比，由

2、构件变形引起的弹性位移很小；Ø这种弹性位移不会影响机构的名义运动。依据上述假设，机构真实运动的位移可以看作是名义运动的位移和弹性位移的叠加。名义运动可以用刚体机构运动和动力学分析方法求出，弹性位移则用弹性动力学分析方法求出。为了使所建模型较准确反应原机构系统的特性，现在普遍采用“子结构分析方法”，即把系统按结构划分为子结构单元，然后建立单元和子结构的运动方程，最后将单元和子结构的运动方程组合成系统的运动方程。对于连续体的离散，有1）集中参数模型2）有限元模型两种建模方法。以一个简单例子为例：一般弹性动力学方程为：其中，第一个方程描述的是机构的刚体动力学方程，第二个方程描述的是机构

3、的结构振动方程。表示机构广义刚体位移，表示机构广义弹性位移，表示机构所受外力，表示机构的科氏力和离心力。对于KED方法，变形对刚体运动的影响忽略不计，因此，忽略耦合项，上述方程变为：从上式可以看出，由于KED方法的假设，使方程得到很大的化简，提高了计算效率，此方法对于作大范围刚体运动，机构刚度大（即弹性变形小的系统）适用。但随着轻质、高速运动、大尺寸机构的发展，KED方法计算结果的精确度不再令人满意。在这些系统中，刚体运动和弹性变形的惯性耦合非常重要，在动力学分析中不能被忽略，因此，KED方法在这些机构的动力学分析中不再适用。二、浮动坐标法（FloatingFrameofRefe

4、rence）浮动坐标法是目前进行计算机柔体动力学仿真时最广泛运用的方法，这种方法已经被应用在几种商用动力学分析软件中。在浮动坐标法中，共使用两种坐标系来描述变形体的构型，一种是用来描述变形体连体坐标系的位置和方向，另一种是用来描述变形体相对于其连体坐标系的变形。如下图所示：图一1）运动分析变形体上任意一点在全局坐标系中的位置为：(1.1)其中，为连体坐标系相对于全局坐标系的方向余弦矩阵。为点在未变形时在连体坐标系中的位置，为变形位移。对于一个具体问题，如何选择合适的连体坐标系是难点。在多刚体动力学中，选取通过质心的主轴坐标系为跟随坐标系而使动力学方程中平移与转动惯性解耦。柔体中各

5、质点的位置时刻都在变换，其质心相对于其内部的质点也一直在不停地变化，因而不存在一个固定的连体坐标系，选取不同的连体坐标意味着选取了不同的柔性体变形。但研究表明这并不影响最终的位移分析结构。系统广义坐标为。其中，为描述变形体相对位置的笛卡尔坐标，为描述变形体方向的角度坐标，为变形体上任意点变形坐标，与变形位移的关系为，其中为型函数。之后便可以对系统进行运动学分析，对公式（1.3）求一次导数，便可得到任一点的运动速度：(1.2)其中，，为系统广义速度阵，为系数矩阵。对1.4再求一次导数，可得到点的运动加速度：(1.3)2）质量矩阵：（1）系统的动能为：(1.4)（2）系统质量阵为(1

6、.5)是一个非线性的对称矩阵。变系数，随位形变化3）系统广义力利用虚功原理，求解弹性力和外力所产生的关于广义坐标的广义力（1）系统广义弹性力(1.6)（2）系统广义外力(1.7)其中，和为关于移动和转动坐标的广义力。4）运动约束方程图二系统的约束方程可写为向量的形式：（1.10）其中，为系统的广义坐标，为时间，为独立的约束方程。例如，如图二所示，如果点和点相连，则有(1.8)当时，则表示两个点始终相连。为了将约束方程（1.10）引入动力学方程中，对于广义坐标取微小变化，公式（1.10）可写为：(1.9)其中，为系统约束压雅可比矩阵。5）系统动力学方程将以上求解各式带入到第一类拉格

7、朗日方程中，即可得到系统动力学方程(1.10)其中，为拉格朗日乘子。从以上分析可以看出，在浮动坐标法中，质量矩阵为一个非线性的对称矩阵，刚度矩阵为一个常量矩阵。计算广义力时，需考虑系统的科氏力和离心力。浮动坐标法适用于作大范围移动小变形的系统，对于作大范围移动大变形的系统，此方法的求解不够精确，不再适用。三、绝对节点坐标法（AbsoluteNodalCoordinateFormulation）该方法由AhmedA.Shabana于1996年提出，其理论基础主要是有限元与连续介质力