您的当前位置:首页正文

ADAMS在求解微分方程组中的应用

2022-04-09 来源:伴沃教育
ADAMS在求解微分方程组中的应用

众所周知ADAMS具有强大的结算功能,在求解动力学问题方面可谓得心应手。在此我想介绍一下它在求解非线性微分方程组方面的应用。

在工程分析中我们经常用一组微分方程来对一个过程进行描述。我所熟悉的内外弹道现象最后都归纳为一个微分方程组。这样最后对问题的解决也就转移到对方程组的求解上了。而这一类微分方程绝大部分不能给出解析解,只能采用数值解法来求解近似解,常用的一种便是龙格库塔法了。这又牵扯到编程,界面的组织及结果的处理,虽然解决方法明了,但实际操作是很麻烦的。我以一个简单的例子说明用ADAMS来解决类似的方程组。在这个过程中,我们无需编程,不需要了解什么数值解法,只需把ADAMS的解算器看成一个黑匣子就可以了。

假设一个物体在地面上运动,它由于受到空气阻力会做减速运动。而空气阻力是与物体运动速度相关联的,它不像摩擦系数那样是一个稳定的值或系数。也就是说空气阻力是速度的一个函数。假设关系为:F=f(v)。在此不考虑摩擦力,那么对物体的运动过程可以用以下方程组来描述:

v=

dx

dtdvFf(v)

a==−=−

dtmm

上述的方程组有解,并且解是唯一的。现在我们用ADAMS来解此问题。首先建一个

任意形状的物体,比如一个块体,修改其质量,使其为m,赋予其初速v0,这都为初始条件。然后给其加一与运动方向相反的力F=ma=f(v)。然后一切OK了,接着进行仿真,ADAMS便会自动求解该运动过程了。在后处理中可以看到各运动诸元了,自己就无需为结

果的处理而犯愁了。

上面只是举一个简单的例子。有的虽然复杂,但解决思路是一样的。象我在求解外弹道时,弹丸受的空气阻力要分解成两个方向,而且还有一个跟速度有关的阻力系数,还有跟弹丸高度有关的函数式,虽然有些复杂,但只要理顺好了关系,在弹丸上施加一个力的表达式,而这个表达式应该是弹丸运动过程中的一个可测量的值,比如速度,加速度或质心位置等就可以了。最后我还设计了一个界面,可以解任意的弹丸外弹道问题。界面如下所示:

输入初始条件,按“确定”键就可进行计算了,比起自己编程,省事多了,而且绝对准确。

另外由于ADAMS不支持汉字,因此我们在写自己的对话框时,为了使各参数更明确,在定义Label时在“Label text”里不要填,而要把汉字说明保存成bmp图片形式,再在“Icon file”里写出关联的图片地址,如下所示,虽然有些麻烦,但明了直观。

希望对诸位有益。

Nijinf于南京理工大学 2004-11-1

因篇幅问题不能全部显示,请点此查看更多更全内容