现代接触动力学1

《现代接触动力学1》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、现代接触动力学其中为刚体的惯性张量，代表所有作用在刚体上的力之和，包括主动力矩之和和理想约束力矩之和。动量力矩方程是相对一个确定的参考点而建立的。通常这个参考点取为刚体的质心。如果参考点变动，式（2．4）中所有的量必须相应地改变。虽然理想约束力和力矩由于运动方向正交而不影响物体的运动，但是他们却出现在动量和动量守恒方程中。虚功原理则指出理想约束力和力矩的虚功为零。即（2．25）在静平衡状态下主动力的虚功也为零，，这个关系式能被用于确定静平衡的条件。拉格朗日形式的达朗伯原理相应于动力学问题的虚功原理，利用式（2．23）和（2．24）求解和，

2、并将它们带入式（2．25）中得（2．26）用广义坐标上式可转化为(2．27)由的任意性可导出完整约束多刚体系统的运动微分方程（2．28）其中为对称正定的质量矩阵，代表广义的陀螺力、离心力和科氏力，代表主动作用力。利用广义坐标，动量和动量矩守恒定律也可用分块矩阵形式加以描述为（2．29）其中——质量和转动惯量分块对角矩阵；——全局雅克比矩阵；——陀螺力，离心力，科氏力；现代接触动力学——主动作用力和力矩；——约束作用力和力矩。牛顿-欧拉方程式（2．29）中出现的约束作用力和力矩可通过左乘全局雅克比矩阵而消除，因为根据达朗伯原理可导出。左乘矩

3、阵后，牛顿-欧拉方程式（2．29）变成运动方程式（2．28）。简洁的运动方程式以及上述的推导只适用于完整约束力多刚体系统。对于一般的约束方程与速度有关的非完整约束多刚体系统运动学微分方程可写为（2．30）动力学微分方程类似于（2．28），即（2．31）2．3多刚体系统计算软件和求解过程多刚体系统的建模通常是借助于计算机来实现的。这里只简单地介绍一些斯图加特大学B所开发的计算机软件。在软件包NEWMOD（参见EberhardNeerpasch1996）中，简单的基本元素诸如坐标系统，支座和刚体等能通过交互界面直观地描述，见图2．3。系统的若

4、干基本元素在软件包NEWMOD中被装配在一起，然后输入到自动生成运动微分方程的符号推导软件包NEWEUL（参见Eberhard1979）中，为便于应用求解常微分方程组的标准软件，个二阶微分方程组被转化为个一阶微分方程组（2．32）通过对时间积分可获得广义位置和速度在离散时间点上的数值，动画显示软件包如NEWANIM(参见[Kberhard1994])可形象地表示这些数值结果，展示系统的运动过程。运动微分方程亦可进一步地用于分析，比如导出的带有符号的运动微分方程可运用计算机代数系统软件包MAPLE进行线性化。用软件系统MATLAB还可为线性

5、化的系统配置控制器。2．4简单的多刚体系统示例运动方程的推导可通过示例加以说明。一个质量为现代接触动力学，标号为1的刚体在一个倾斜角为的光滑斜面上向下滑动，见图2．24；另一个摆长为，质量为，标号为2的单摆围绕着刚体1的质心以角速度匀速转动，该系统只有一个自由度，广义坐标可选为，矢径列阵为（2．33）图2.4简单的多刚体系统示例雅克比矩阵为（2．34）由式和可得相对速度和相对加速度（2．35）主动作用力只有重力（2．36）根据式（2．27），运动方程中式（2．28）出现的各项一次变为（2．37）现代接触动力学（2．38）（2．39）运动方

6、程为（2．40）设初始条件，参数通过求解微分方程式（2．40），图（2．5）展示广义位移和广义速度随时间的变化规律。图2．6动画显示该系统的运动过程。在图2．6中不同的时间的运动状态重叠在一幅图中，以便于对整体的运动有更清晰的直观认识。尽管这个例子和第5章中5．6节所举例的带有接触的多体系统看上去十分类似，然而约束的本质却大相径庭。这里刚体1和斜面的约束是双面约束，不仅刚体1不可能脱离斜面，而且斜面也允许承受拉应力。而在5．6节中所讨论的接触问题则涉及单面约束，刚体1允许脱离斜面，并且在刚体1和斜面的接触过程中，斜面不允许承受拉应力。3线

7、性有限元法有限元法可看作为求解偏微分方程的数值计算方法。其基本原理是将几何形状复杂的求解域离散为一组有限个且按一定方式相互联结在一起的单元的组合体。尽管现代有限元法起源于结构力学，但如今它已在许多计算领域内如流体力学、热传导、空气动力学中得到广泛的应用。接触问题源于固体力学。计算接触问题时通常是先引进固体力学中线弹性理论的基本方程，然而讨论求解这些基本方程的近似解法，最后常用加权余量法和积分弱解形式建立可解的线性方程组。由于几何形体常离散成几种基本类型的单元的组合体，人们只需要对这几种基本的单元进行分析，建立单元刚度矩阵和质量矩阵，再考虑

8、单元的应力、应变以及边界条件，然后将所有的结构单元组装在一起。接触问题总是非线性的，然而线性有限元的许多部分可以直接地或通过修正运用于非线性有限元计算。这一章介绍的线性有限元法将作为非线性有限