第4章ADAMS软件算法基本理论-(陈立平)机械系统动力学分概要文档格式.doc
- 文档编号:13010062
- 上传时间:2022-10-02
- 格式:DOC
- 页数:17
- 大小:411KB
第4章ADAMS软件算法基本理论-(陈立平)机械系统动力学分概要文档格式.doc
《第4章ADAMS软件算法基本理论-(陈立平)机械系统动力学分概要文档格式.doc》由会员分享,可在线阅读,更多相关《第4章ADAMS软件算法基本理论-(陈立平)机械系统动力学分概要文档格式.doc(17页珍藏版)》请在冰豆网上搜索。
对于大多数问题,可以将地球近似为惯性参考标架,虽然地球是绕着太阳旋转而且地球还有自转。
对于每一个刚性体都有一个与之固定的参考标架,称为构件参考标架,刚性体上的各点相对于该构件参考标架是静止的。
4.1.2坐标系的选择
机械系统的坐标系广泛采用直角坐标系,常用的笛卡尔坐标系就是一个采用右手规则的直角坐标系。
运动学和动力学的所有矢量均可以用沿3个单位坐标矢量的分量来表示。
坐标系可以固定在一个参考标架上,也可以相对于参考框架而运动。
合理地设置坐标系可以简化机械系统的运动分析。
在机械系统运动分析过程中,经常使用3种坐标系:
(1)地面坐标系(GroundCoordinateSystem)。
地面坐标系又称为静坐标系,是固定在地面标架上的坐标系。
ADAMS中,所有构件的位置、方向和速度都用地面坐标系表示。
(2)局部构件参考坐标系(LocalPartReferenceFrame,LPRF)。
这个坐标系固定在构件上并随构件运动。
每个构件都有一个局部构件参考坐标系,可以通过确定局部构件参考坐标系在地面坐标系的位置和方向,来确定一个构件的位置和方向。
在ADAMS中,局部构件参考坐标系缺省与地面坐标系重合。
(3)标架坐标系(MarkerSystem)。
标架坐标系又称为标架,是为了简化建模和分析在构件上设立的辅助坐标系,有两种类型的标架坐标系:
固定标架和浮动标架。
固定标架固定在构件上,并随构件运动。
可以通过固定标架在局部构件参考坐标系中的位置和方向,确定固定标架坐标系的位置和方向。
固定标架可以用来定义构件的形状、质心位置、作用力和反作用力的作用点、构件之间的连接位置等。
浮动标记相对于构件运动,在机械系统的运动分析过程中,有些力和约束需要使用浮动标架来定位。
动力学方程的求解速度很大程度上取决于广义坐标的选择。
研究刚体在惯性空间中的一般运动时,可以用它的质心标架坐标系确定位置,用质心标架坐标相对地面坐标系的方向余弦矩阵确定方位。
为了解析地描述方位,必须规定一组转动广义坐标表示方向余弦矩阵。
第一种方法是用方向余弦矩阵本身的元素作为转动广义坐标,但是变量太多,同时还要附加六个约束方程;
第二种方法是用欧拉角或卡尔登角作为转动坐标,它的算法规范,缺点是在逆问题中存在奇点,在奇点位置附近数值计算容易出现困难;
第三种方法是用欧拉参数作为转动广义坐标,它的变量不太多,由方向余弦计算欧拉角时不存在奇点。
ADAMS软件用刚体的质心笛卡尔坐标和反映刚体方位的欧拉角作为广义坐标,即,。
由于采用了不独立的广义坐标,系统动力学方程虽然是最大数量,但却是高度稀疏耦合的微分代数方程,适用于稀疏矩阵的方法高效求解。
4.2ADAMS运动学分析
4.2.1ADAMS运动学方程
利用ADAMS建立机械系统仿真模型时,系统中构件与地面或构件与构件之间存在运动副的联接,这些运动副可以用系统广义坐标表示为代数方程,这里仅考虑完整约束。
设表示运动副的约束方程数为,则用系统广义坐标矢量表示的运动学约束方程组为:
((4.2--1)
考虑运动学分析,为使系统具有确定运动,要使系统实际自由度为零,为系统施加等于自由度()的驱动约束:
…………………………………………………………(4.2--2)
在一般情况下,驱动约束是系统广义坐标和时间的函数。
驱动约束在其集合内部及其与运动学约束合集中必须是独立和相容的,在这种条件下,驱动系统运动学上是确定的,将作确定运动。
由式(4.2-1)表示的系统运动学约束和式(4.2-2)表示的驱动约束组合成系统所受的全部约束:
…………………………….………...(4.2-3)
式(4.2-3)为nc个广义坐标的nc个非线性方程组,其构成了系统位置方程。
对式(4.2-3)求导,得到速度约束方程:
………………………….…(4.2-4)
若令,则速度方程为:
………...………………………..…(4.2-5)
对式(2.24-4)求导,可得加速度方程:
…(4.2-6)
若令,则加速度方程为:
…………………………(4.2-7)
矩阵,为雅可比矩阵,如果的维数为m,q维数为n,那么维数为矩阵,其定义为。
在这里为(nh个运动学约束,nc-nh个驱动约束,nc个广义坐标)的方阵。
4.2.2ADAMS运动学方程的求解算法
在ADAMS仿真软件中,运动学分析研究零自由度系统的位置、速度、加速度和约束反力,因此只需求解系统的约束方程:
(4.2-8)
运动过程中任一时刻位置的确定,可由约束方程的Newton-Raphson迭代法求得:
(4.2-9)
其中,,表示第次迭代。
时刻速度、加速度可以利用线性代数方程的数值方法求解,ADAMS中提供了两种线性代数方程求解方法:
CALAHAN方法(由Michigan大学DonaldCalahan教授提出)与HARWELL方法(由HARWELL的IanDuff教授提出),CALAHAN方法不能处理冗余约束问题,HARWELL方法可以处理冗余约束问题,CALAHAN方法速度较快。
(4.2--10)
(4.2-11)
4.3ADAMS动力学分析
4.3.1ADAMS动力学方程
ADAMS中用刚体B的质心笛卡尔坐标和反映刚体方位的欧拉角作为广义坐标,即,令,,。
构件质心参考坐标系与地面坐标系间的坐标变换矩阵为:
((4.3-12)
定义一个欧拉转轴坐标系,该坐标系的三个单位矢量分别为上面三个欧拉转动的轴,因而三个轴并不相互垂直。
该坐标系到构件质心坐标系的坐标变换矩阵为:
((4.3-132)
构件的角速度可以表达为:
((4.3-3-14)
ADAMS中引入变量为角速度在欧拉转轴坐标系分量:
(4.3-4-15)
考虑约束方程则系统的动力学方程,ADAMS利用带拉格朗日乘子的拉格朗日第一类方程的能量形式得到如下方程:
(4.3--165)
T为系统广义坐标表达的动能,为广义坐标,为在广义坐标方向的广义力,最后一项涉及约束方程和拉格朗日乘子表达了在在广义坐标方向的约束反力。
ADAMS中近一步引入广义动量:
(4.3-6-17)
简化表达约束反力为:
((4.3-7-18)
这样方程(4.3-5)-16)可以简化为:
(4.3-8-19)
动能可以近一步表达为:
(4.3--209)
其中M为构件的质量阵,J为构件在质心坐标系下的惯量阵。
将(4.3-819)分别表达为移动方向与转动方向有:
(4.3-10-21)
(4.3-11-22)
其中,。
(4.3-10-21)式可以简化为:
(4.3-12-23)
,由于中包含欧拉角,为了简化推导,ADAMS中并没有进一步推导,而是将其作一个变量求解。
这样ADAMS中每个构件具有如下15个变量(而非12个)和15个方程(而非12个)。
变量:
(4.3-12-24)
方程:
(4.3-13-25)
集成约束方程ADAMS可自动建立系统的动力学方程――微分-代数方程:
(4.3-1426)
其中,P为系统的广义动量;
H为外力的坐标转换矩阵。
为了更好地说明ADAMS的建模过程下面以一个单摆为例进行建模推导。
图4-1单摆示意图
如图所示,单摆的质量为M、惯量为I,杆长为2L,并在O点以转动副与大地相连接约束在大地的OXY平面内。
在单摆质心处建立单摆的跟随坐标系――局部构件参考坐标系Op-Xp-Yp,其坐标在地面坐标系OXY中为(x,y),单摆的姿态角为θ。
系统的动能表达式:
(4-27)
广义动量表达式:
(4-28)
外力表达式:
(4-29)
约束方程:
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- ADAMS 软件 算法 基本理论 陈立平 机械 系统 动力 学分 概要