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

1、重庆大学硕士研究生机械系统动力学课程作业汽车系统的动力学响应简单分析小组成员：李航，王青松，田文昌，谢道坤学号： 20150713058t, 20150702040, 20150702036,20150702051指导教师：何荇兮学院：机械工程学院专业：机械工程班级：传动 3班重庆大学二 O 一五年十一月重庆大学硕士研究生机械系统动力学课程作业汽车系统的动力学响应简单分析一、问题描述汽车系统结构简图以及受力分析简图1、图 2 如示：图 1 汽车系统结构简图图 2 汽车系统受力分析简图系统运动方程如下：mx(c1c2 ) x(l 2 c2l1c1 )(k1 k2 ) x (l 2 k2l1k1

2、)0(1)2( c2 l 2c1l 1 ) x22(k2l 2 k1l1 ) x220(2)mr x(l 2 c2l1 c1 )(l1 k1l 2 k2 )系统结构参数如下：r 20.64m2m 4000kgcc2000N s/m12k1 k220000N/ml 10.9ml 2 1.4m其中， r 为车辆的回转半径。给定的初始条件：xx0。外部冲击力矩： 10 (t) （10 倍单位冲击函数）。试用 MATLAB 里的 ODE45 函数画出 05 秒内的位移 x 与转角的响应。ODE是 MATLAB 专门用于解微分方程的功能函数， 其中 ODE45 求解器属于变步长的一种，采用四阶或五阶 R

3、UNGE KUTTA 法，是解决数值解问题的首选方法。1重庆大学硕士研究生机械系统动力学课程作业二、求解方法已知汽车的运动学方程，在t=0 时，存在一个冲击力矩 10 (t) ，此时运动学方程表示为：mx(c1c2 ) x(l 2c2l1c1 )(k1k2 )x(l 2k2l1k1)0(3)mr 2(c2l2c1l1 ) x(l 22 c2l12c1 )(k2 l2k1l 1) x(l22k2l 12 k1 )10 (t )(4)冲击函数(t) 的定义：本算例中的冲击函数为一个分步函数，由于时间的非负性，冲击函数的定义为：1 ,t0(t )t00 ,根据题目给出的结构参数和初始条件，通过MAT

4、LAB 中的 ODE45函数编写程序求解添加了冲击的运动学方程得出结果，并画出0 5 秒内汽车系统的位移 x 与转角的响应图。三、求解过程MATLAB 中的编程思路及过程如下：首先将运动学方程写成一阶常微分方程组形式：t=0 时：t>0 时：x1x2x2x20.25x410x12.5 x3x3x4x40.0025 0.390625 x22.1640625 x4 3.90625x1 21.640625x3x1x2x2x20.25x410x12.5 x3x3x4x40.390625x22.1640625x43.90625x121.640625x3在 MATLAB 里的 ODE45 函数编写计

5、算主程序，如下所示：clc;clear;t0=0;tf=5;%0<=t<=5x0=0 0 0 0;%定义初始值2重庆大学硕士研究生机械系统动力学课程作业t,x=ode45('f3',t0 tf,x0);%ODE45调用 f3figure;plot(t,x);由 MATLAB 编程可定义函数f 3 ，程序如下所示：function xdot=f3(t,x)%定义函数f3xdot=zeros(4,1);if t<=0xdot=x(2);-(2000+2000)*x(2)+(1.4*2000-0.9*2000)*x(4)+(20000+20000)*x(1)+(1.

6、4*20000-0.9*20000)*x(3)/4000x(4);% 当 t=0 时，作用一个冲击冲击力矩10 (t)elsexdot=x(2);-(2000+2000)*x(2)+(1.4*2000-0.9*2000)*x(4)+(20000+20000)*x(1)+(1.4*20000-0.9*20000)*x(3)/4000x(4);% 当 t>0 时的运动学方程，此时冲击已经结束End最后运行程序即可得到计算结果以及05 秒内汽车系统的位移x 与转角的响应图。3重庆大学硕士研究生机械系统动力学课程作业四、结果与分析经过 MATLAB 运行计算最终得到 05秒内汽车系统的位移 x

7、与转角 的响应图如下：图 3 汽车系统位移x 与转角的响应图由图 3 可看出：在 t=0 时，由于冲击力矩的存在，曲线急剧上升，而冲击力矩对 x 、 x 、的初始状态没有影响，但是使它们在随后的时间产生了波动。同时还可以看出位移x 与转角的响应时间约为4.5s，由于系统存在阻尼，在此之后位移 x 与转角 的响应逐渐趋于平缓，最终趋于无干扰存在时的静止状态， 可见静止状态是一种稳态。4重庆大学硕士研究生机械系统动力学课程作业五、总结与感想在经过这学期的机械系统动力学课程学习，我们收获很多了，学习了不少的内容。首先，我们学习了质点动力学，其中我们主要学习了五大坐标系：笛卡尔坐标系、柱坐标系、极坐标

8、系、球面坐标系还有路径坐标系。通过对这五大坐标系的学习，我们知道了如何在对一个实际的模型进行分析求解时，选择正确的坐标系。然后，我们学习了动力学，其中主要学习了矢量对时间的导数的求法，单位矢量在旋转系统中的导数的求法，同时还学习了刚体的动力学分析。接着，我们学习了相对运动，其中学习了不同坐标系之间的转换，以及不同坐标上单位矢量之间的转换方法，还学习了力与扭矩的平衡问题。最后，我们学习了拉格朗日方程，其中我们学会了求解虚位移、虚功、学会了求解拉格朗日方程的具体步骤。通过，前面内容的学习，我们学会了机械系统动力学的基本知识，学到了简单机械系统运动方程的建立与求解。在课程结束后， 经过老师的指导我们小组四个人， 利用 MATLAB 结合相关机械动力学知识一起完成了课程大作业。我们通过这次汽车系统的动力学的响应简单分析，使用 MATLAB 中的 ODE45 函数求解一阶常微分方程组，通过对比加入初始冲击转矩与存在初始速度时的不同的响应图以及系统的稳定状态，加深了对冲击、阻尼及动力学方程的认识， 这类似于通过 MATLAB 对汽车系统振动进行