欢迎来到天天文库
浏览记录
ID:41764847
大小:140.50 KB
页数:6页
时间:2019-09-01
《三种汽车振动模型的MATLAB分析比较》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库。
1、三种汽车振动模型的MATLAB分析比较赖永生工程力学专业学号:0743058043一、背景及意义1.没有振动的汽车是不存在的,不同汽车的振动大小也是不一样的。2.汽车的振动问题具有强烈的系统依赖性、时间依赖性以及耦合关联性。3.为了避免成本更高、耗时更长而效果有限的生产后的整改,在设计开发的早期就进行系统设计与控制,则能够达到良好的减振效果。4.实现“早期就进行系统设计与控制”也是一件不容易的工作。尤其在以前模拟和仿真系统欠缺与不完善的情况下更是如此,进行汽车振动分析“成本更高、耗时更长而效果有限”。5.MATLAB的出现完全改变了这种窘况,用MATLAB进行汽车振动分析“成本更低、
2、耗分析比较。时更短而效果显著”,而且更加方便快捷。6.但是不管进行什么分析,都要先建立模型,不同的模型,得到的结果有一定差异。所以有必要分别对每个模型分别进行分析,以确立最佳模型。本文就是对三种汽车振动模型的MATLAB分析比较。二、具体内容本文首先给岀每一种具体的汽车振动模型,包括图示、振动方程、初始条件、边界条件和相关参数,充分了解并利用相应的MATLAB的库函数,然后编制程序求解并画出振动波形图,最后对三种模型进行分析比较,得出结论。图1单自由度系统的振动模型三、技术路线及技术难点和可行性分析技术路线:tny+cy+ky=Q化为其中»+阿+©J=o是固有频率5是阻尼比(1)单个
3、自由度的振动模型单自由度阻尼系统的振动方程如下:方程解为:X0=4幺缽sin(©f+9)其中(v°+曲)尸+g©)2,0=0(弩),©=%冋y()-初始位置心-初始速度这里假设^=0.1到1,公共参数^=2O,xo=l,心=0。计算的终点时间tf=3。(2)双自由度的振动模型(车体上下振动图2和绕质心的旋转图3)图3二自由度系统的振动cmx=-k、(x+l})-k2(x-l23)
4、^
5、)x+伙]/]"+)0—0令心+灯3k7kJ?—kJ、kl—,b=—mmmmk°—k
6、J、3k』kJ、+k2】「9k—9d—Jcml得rx+ax-b0=0〔0-ex+dO—0(3)二自由度可解耦系统的振动模型,其微分方程如下所示:(I)(X)//7q////C2//图4二自由度可解耦系统的振动模型mxx{+(q+c2)Xj一c2x2++k2)x{-k2x2=0m{x2+c2(x2-xl)+k2(x2-xj=0可写成矩阵形式Mx+Cx+Kx=0其中,Mm]0K=W,X士叮tv27、a2_②在原始方程加+Kx二0两端各作乘以U及在屮间乘以对角矩阵,X.=U*1=Ux得*M*(/*/*壬+(/*K*U*(/*x=0_Z]作变量置换Z—7也2z+L*z=0这是一个对角矩阵方程,即可把它分成两个方程求解:Zj+右Z]=0z9+A2z2=0这就是说两种振动模态可解耦,令©2=&,①是第i个模式的固有频率(i二1,2)。③由上述的解耦模态中,给出初始条件兀0,兀〃0,化为z°,乙必,即可分别求出其分量z2(此处阻尼比为0),再按兀=inv(U)*z变换为X。假设位置和速度的初始条件分别为_[IT_[1T2x(t)=^ui([ui]TZ=1X0=lX0]XQ29Xd0=l8、XdOXd02,则这三个步骤得到最后公式为1丁mx0cosco^H[wz]Mxd0sin(4)编写程序上机调试求解并画图,对三种模型的结果进行比较分析,得出和应结论。技术难点:(1)当阻尼不等于零时,用MATLAB无法求解,需要借助数值方法求解。(2)如果选取不同的广义坐标,则建立的微分方程的形式也不同,方程的耦合情况和输出相应也将发生相应变化。(3)对于汽车振动的的两种模态,当两种模态叠加时,初始条件变化,则模态的幅值也将发生相应变化。(4)可能会遇到程序的编制困难、调试困难和图像显示困难。可行性分析:只要建立了合乎实际的的模型,并列出了汽车振动的微分方程、初始条件和边界条件,9、再利用MATLAB的相关函数求解微分方程和画出振动波形图,问题也就得到了解决。需要说明的是,这里选取的汽车振动模型已是经过实践检验的和具体实验结果吻合得比较好的模型,所以这里省略了和具体实验数据的比较,仅就三种不同汽车模型进行比较。四、时间安排14周查阅相关的文献。15周完成程序编制并上机调试运行,保存结果。16周完成比较分析和最后作业报告。五、成果形式我将以程序的形式提交课题成果,并配以一定量分析结果的图片。下面是部分成果展示:10X因匡-2012345
7、a2_②在原始方程加+Kx二0两端各作乘以U及在屮间乘以对角矩阵,X.=U*1=Ux得*M*(/*/*壬+(/*K*U*(/*x=0_Z]作变量置换Z—7也2z+L*z=0这是一个对角矩阵方程,即可把它分成两个方程求解:Zj+右Z]=0z9+A2z2=0这就是说两种振动模态可解耦,令©2=&,①是第i个模式的固有频率(i二1,2)。③由上述的解耦模态中,给出初始条件兀0,兀〃0,化为z°,乙必,即可分别求出其分量z2(此处阻尼比为0),再按兀=inv(U)*z变换为X。假设位置和速度的初始条件分别为_[IT_[1T2x(t)=^ui([ui]TZ=1X0=lX0]XQ29Xd0=l
8、XdOXd02,则这三个步骤得到最后公式为1丁mx0cosco^H[wz]Mxd0sin(4)编写程序上机调试求解并画图,对三种模型的结果进行比较分析,得出和应结论。技术难点:(1)当阻尼不等于零时,用MATLAB无法求解,需要借助数值方法求解。(2)如果选取不同的广义坐标,则建立的微分方程的形式也不同,方程的耦合情况和输出相应也将发生相应变化。(3)对于汽车振动的的两种模态,当两种模态叠加时,初始条件变化,则模态的幅值也将发生相应变化。(4)可能会遇到程序的编制困难、调试困难和图像显示困难。可行性分析:只要建立了合乎实际的的模型,并列出了汽车振动的微分方程、初始条件和边界条件,
9、再利用MATLAB的相关函数求解微分方程和画出振动波形图,问题也就得到了解决。需要说明的是,这里选取的汽车振动模型已是经过实践检验的和具体实验结果吻合得比较好的模型,所以这里省略了和具体实验数据的比较,仅就三种不同汽车模型进行比较。四、时间安排14周查阅相关的文献。15周完成程序编制并上机调试运行,保存结果。16周完成比较分析和最后作业报告。五、成果形式我将以程序的形式提交课题成果,并配以一定量分析结果的图片。下面是部分成果展示:10X因匡-2012345
此文档下载收益归作者所有