汽车平顺性与仿真分析matlab,基于matlab的汽车平顺性的建模与仿真.docx

基于matlab的汽车平顺性的建模与仿真.docx

(1)基于MATLAB的汽车平顺性的建模与仿真车辆工程专硕1601Z1604050李晨1数学建模过程11建立系统微分方程如下图所示，为车身与车轮二自由度振动系统模型图中，M2为悬挂质量(车身质量)；M1为非悬挂质量(车轮质量)；K为弹簧刚度；C为减振器阻尼系数；KT为轮胎刚度；Z1为车轮垂直位移；Z2为车身垂直位移；Q为路面不平度。车轮与车身垂直位移坐标为Z1、Z2，坐标原点选在各自的平衡位置，其运动方程为221210MZCZKZ11110TKZQ(2)(3)(4)(5)(6)12双质量系统的传递特性先求双质量系统的频率响应函数，将有关各复振幅代入，得令23TAMJCK由式(2)得Z2Z1的频率响应函数将式(4)代入式(3)得Z1Q的频率响应函数NKAAAAKAQZTT2212321式中下面综合分析车身与车轮双质量系统的传递特性。车身位移Z2对路面位移Q的频率响应函数，由式(4)及(5)两个环节的频率响应函数相乘得到221ZMJCKZJCK11TTQJ21212AJZKJC2321N21211TTAKZZQN(7)(11)(8)(7)()(9)(10)(12)13车身加速度、悬架弹簧动挠度和车轮相对动载的幅频特性1车身加速度对路面不平度的频率特性2相对动载对路面不平度的频率特性车轮动载荷为车轮静载荷为则车轮与路面相对动载为车轮与路面间相对动载与路面不平度之间的传递函数为3悬架动挠度对路面不平度的频率特性悬架动挠度为222ZQZZHQ12DFMZZ12GG2112121DMZFMZGGG1221/21DDFGQZMQHG2121DFZZQQ(13)悬架动挠度与路面不平度之间的传递函数为2仿真过程通过建模，我们已经得到了各所需的传递函数。下面要利用MATLAB的M文件进行仿真。21公式的进一步推导在公式(7)中，我们需要得到的是传递函数的分子和分母表达式，这样可以通过插值的方法计算传递函数，并以此计算出幅频特性。经进一步推导后我们可得公式(7)的分子为32TJCK分母为同理，对公式(11)、(13)进行推导得公式(11)分子为43222211TMMMJCK分母为21DDFQFZHQ2321NA21MGN公式(13)分子为2TK分母为N22M文件中代码的编写得到了所有传递函数的分子、分母，下面编写代码1一些系统参数的输入2传递函数分子、分母的构建3传递响应函数的构建及频响输出车身加速度对路面不平度响应特性悬架动挠度对路面不平度响应特性相对动载对路面不平度响应特性23图形输出对比汽车理论教材上的内容，作出的曲线基本符合。