重庆大学硕士研究生机械系统动力学课程作业.doc

重庆大学硕士研究生机械系统动力学课程作业重庆大学硕士研究生机械系统动力学课程作业汽车系统的动力学响应简单分析小组成员： 李航，王青松，田文昌，谢道坤 学 号：20150713058t, 20150702040, 20150702036,20150702051 指导教师： 何荇兮 学 院： 机械工程学院 专 业： 机械工程 班 级： 传动3班 重庆大学二O一五年十一月 1重庆大学硕士研究生机械系统动力学课程作业汽车系统的动力学响应简单分析一、 问题描述汽车系统结构简图以及受力分析简图1、图2如示：5图1 汽车系统结构简图图2 汽车系统受力分析简图系统运动方程如下： (1) (2)系统结构参数如下：其中，为车辆的回转半径。给定的初始条件：。外部冲击力矩：（10倍单位冲击函数）。试用MATLAB里的ODE45函数画出05秒内的位移与转角的响应。ODE是MATLAB专门用于解微分方程的功能函数，其中ODE45求解器属于变步长的一种，采用四阶或五阶RUNGEKUTTA法，是解决数值解问题的首选方法。二、求解方法已知汽车的运动学方程，在t=0时，存在一个冲击力矩，此时运动学方程表示为： (3) (4)冲击函数的定义：本算例中的冲击函数为一个分步函数，由于时间的非负性，冲击函数的定义为：根据题目给出的结构参数和初始条件，通过MATLAB中的函数编写程序求解添加了冲击的运动学方程得出结果，并画出05秒内汽车系统的位移与转角的响应图。三、求解过程MATLAB中的编程思路及过程如下：首先将运动学方程写成一阶常微分方程组形式：t=0时：t0时：在MATLAB里的ODE45函数编写计算主程序，如下所示：clc;clear;t0=0;tf=5; % 0=t=5x0=0 0 0 0; % 定义初始值t,x=ode45(f3,t0 tf,x0); % 调用figure; plot(t,x);由MATLAB编程可定义函数，程序如下所示：function xdot=f3(t,x) %定义函数xdot=zeros(4,1);if t0时的运动学方程，此时冲击已经结束End最后运行程序即可得到计算结果以及05秒内汽车系统的位移与转角的响应图。四、结果与分析经过MATLAB运行计算最终得到05秒内汽车系统的位移与转角的响应图如下：图3 汽车系统位移与转角的响应图由图3可看出：在t=0时，由于冲击力矩的存在，曲线急剧上升，而冲击力矩对、的初始状态没有影响，但是使它们在随后的时间产生了波动。同时还可以看出位移与转角的响应时间约为4.5s，由于系统存在阻尼，在此之后位移与转角的响应逐渐趋于平缓，最终趋于无干扰存在时的静止状态，可见静止状态是一种稳态。五、总结与感想在经过这学期的机械系统动力学课程学习，我们收获很多了，学习了不少的内容。首先，我们学习了质点动力学，其中我们主要学习了五大坐标系：笛卡尔坐标系、柱坐标系、极坐标系、球面坐标系还有路径坐标系。通过对这五大坐标系的学习，我们知道了如何在对一个实际的模型进行分析求解时，选择正确的坐标系。然后，我们学习了动力学，其中主要学习了矢量对时间的导数的求法，单位矢量在旋转系统中的导数的求法，同时还学习了刚体的动力学分析。接着，我们学习了相对运动，其中学习了不同坐标系之间的转换，以及不同坐标上单位矢量之间的转换方法，还学习了力与扭矩的平衡问题。最后，我们学习了拉格朗日方程，其中我们学会了求解虚位移、虚功、学会了求解拉格朗日方程的具体步骤。通过，前面内容的学习，我们学会了机械系统动力学的基本知识，学到了简单机械系统运动方程的建立与求解。在课程结束后，经过老师的指导我们小组四个人，利用MATLAB结合相关机械动力学知识一起完成了课程大作业。我们通过这次汽车系统的动力学的响应简单分析，使用MATLAB中的ODE45函数求解一阶常微分方程组，通过对比加入初始冲击转矩与存在初始速度时的不同的响应图以及系统的稳定状态，加深了对冲击、阻尼及动力学方程的认识，这类似于通过MATLAB对汽车系统振动进行了一次仿真，而通过响应图可以更直观、更准确地分析系统振动的位移、角度、速度以及角速度的变