数值分析在机械领域的应用.ppt

《机械运动的数值仿真》研究 学校：河北联合大学研究生院—机械工程学院 班级：2011级研究生5班 专业：机械设计及自动化 第21组主讲：杨超 学习的主要内容 一、机械运动相关内容 二、数值分析算法的研究 三、常微分方程初值问题数值求解 四、龙格-库塔算法 五、龙格-库塔方法的实际应用 六、报告总结 一、机械运动相关内容 1、机械运动的相关简单概念　 （1）机械运动 　　 在物理学中，把一个物体相对于另一个物体位置的变化称作为机械运动，简称运动。 　　 （2）参照物 　　 要判断一个物体是否在运动，必须选择另一个物体作为标准，这个作为标准的物体叫做参照物。对于同一个物体的运动，选择的参照物不同，得出的结论也有可能是不同的。 　　 （3）运动和静止的相对性 　　 自然界中一切物体都在运动，因为地球本身在自转，所以绝对静止的物体是不存在的。通常所描述的物体的运动或静止都是相对于某一个参照物而言的。 2、机械运动的前沿科学—导弹的制导 （1）导弹制导的一般原理 在大气层内飞行的导弹，可由改变空气动力获得控制，有翼导弹一般用改变空气动力的方法来改变控制力。 在大气层中或大气层外飞行的导弹，都可以用改变推力的方法获得控制。无翼导弹主要是用改变推力的办法来改变控制力，因无翼导弹在稀薄大气层内飞行时，弹体产生的空气动力很小。 三、常微分方程初值问题数值求解 1、导弹轨迹的描述涉及到常微分方程初值问题数值求解问题 2、《常微分方程初值问题数值求解》方法有多少？ Euler方法、向后Euler方法、梯形方法、改进Euler方法、龙格-库塔方法 3、经典的《常微分方程初值问题数值求解》方法是什么方法？Euler方法 4、《常微分方程初值问题数值求解》的优缺点分析 Euler方法计算简单但精度差；向后Euler方法与Euler方法误差相似；梯形方法比Euler方法精度高但算法复杂、计算量很大；改进Euler方法结合了Euler方法和梯形法的优点； 5、最好的方法是？龙格-库塔方法 四、龙格-库塔算法 1、龙格-库塔算法应用举例 五、龙格-库塔方法的实际应用 如下图4.1所示，假设有一烟花火箭，其初始条件为零。将其放在地方然后点火，该烟花火箭的初始质量为，其中粉末燃料占。经过实验得知，燃料的持续时间为。燃料所产生的恒定推力为。这也说明燃料的消耗率恒定。空气产生的阻力和烟花火箭的速度的平方成正比：。这里，要求选择一种数值方法对其运动过程进行仿真并且其截断误差为或者更高。要求计算出该烟花火箭的最高高度，同时计算出从燃料消耗到该烟花火箭运动到最高点的时间延迟。 该实例要求其截断误差要求大于或等于，这就使得较为简单的欧拉法，中点法不适合本例。龙格-库塔法以其优异的数值特性成为解决本问题的首选，且该数值算法很容易Mmathematic中实现。 很显然，该问题属于变质量的运动学问题，在该运动过程中，其前两秒是在驱动力和阻力的共同作用下加速上升的，而后的时间内，该烟花火箭是在空气的阻力下减速上升的，同时注意到空气的阻力和速度的平方成正比。为了对该运动过程进行数值仿真，那么必须建立相应的微分方程组。分析该运动过程可知，应该将该运动过程分为两部分：加速上升过程和减速上升过程。从而得到相应的微分方程组。 加速上升过程： 式中，为上升的高度，为上升过程的速度，为重力加速度。表示加速上升过程的最终高度，表示加速上升过程的最终速度。 使用龙格-库塔法求解如上的微分方程组。该系统的数值仿真结果已在编程体现。如图是该烟花火箭的上升过程高度的数值仿真，如图是其上升过程速度的数值仿真。同时亦可以得到烟花火箭上升的最大高度和问题中所需的时间延迟： Mmathematic仿真图形 六、报告总结 本文观看结束！！！ 祝各位身体健康！万事如意！！ * * * * 二、数值分析算法的研究 1、数值分析方法意义 数学是一种工具，用于解决日常生活、工业工程上的相关问题。针对于数值分析中的数学方法，我们小组将主要内容概括分解，将使用到的方法进行对比分析。 2、我们要求数值算法的稳定性排除病态 算法的数值稳定性与病态问题： 1）若某算法受初始误差或运算过程中的舍入误差影响较小，则称为数值稳定。 2）若微小的初始误差都会对最终结果产生极大的影响，则称之为病态问题。 3、数值分析主要部分。 1各类插值方法我们讲过拉格朗日插值、牛顿插值、分段插值、样条插值 2函数逼近及拟合3数值积分、欧拉法解常微分方程、龙格-库塔法解常微分方程、方程组。 【1】插值 对于牛顿插值相对于拉氏插值增加一个节点，所有的插值基本多项式要重新取、重新算.2而牛顿插值，节点增加，次数增加，即高次插值函数计算量大，有剧烈震荡，数值稳定性较差（例如龙格现象）