三种汽车振动模型的MATLAB分析比较.docxVIP

三种汽车振动模型的MATLAB分析比较 赖永生 工程力学专业 学号：0743058043 一、 背景及意义 没有振动的汽车是不存在的，不同汽车的振动大小也是不一样的。 汽车的振动问题具有强烈的系统依赖性、时间依赖性以及耦合关联性。 为了避免成本更高、耗时更长而效果有限的生产后的整改，在设计开发的早期就进 行系统设计与控制，则能够达到良好的减振效果。 实现“早期就进行系统设计与控制”也是一件不容易的工作。尤其在以前模拟和仿 真系统欠缺与不完善的情况下更是如此，进行汽车振动分析“成本更高、耗时更长而效果有 限”。 MATLAB的出现完全改变了这种窘况，用MATLAB进行汽车振动分析“成本更低、耗 分析比较。时更短而效果显著”，而且更加方便快捷。 但是不管进行什么分析，都要先建立模型，不同的模型，得到的结果有一定差异。 所以有必要分别对每个模型分别进行分析，以确立最佳模型。本文就是对三种汽车振动模型 的MATLAB分析比较。 二、 具体内容 本文首先给岀每一种具体的汽车振动模型，包括图示、振动方程、初始 条件、边界条件和相关参数，充分了解并利用相应的MATLAB的库函数，然 后编制程序求解并画出振动波形图，最后对三种模型进行分析比较，得出结 论。 图1单自由度系统的振动模型三、技术路线及技术难点和可行性分析 图1单自由度系统的振动模型 技术路线： tny + cy + ky = Q化为其中? +阿+ ?J = o是固有频率5是阻尼比单个自由度的振 tny + cy + ky = Q 化为 其中 ? +阿+ ?J = o 是固有频率 5 是阻尼比 方程解为：X0 = 4幺缽sin(?f+9) 其中（v° + 曲）尸 +g?）2,0=0（弩）,? =%冋 其中 y（）-初始位置 y（）-初始位置 心-初始速度 这里假设^=0.1到1,公共参数^=2O,xo=l,心=0。计算的终点时间tf=3。 （2）双自由度的振动模型（车体上下振动图2和绕质心的旋转图3） 图3二自由度系统的振动 图3二自由度系统的振动 cmx = -k、（x + l｝）- k2（x- l23） |^/C0 = —kJ】（X + 厶） + （兀—‘2） 整理得 诚+(々 +k2)x-0(k2l2 一紹])=0 JcO _ (心厶 _ k、Z| )x+伙]/] + )0 — 0 令 心+灯 3 k 7 kJ? — kJ、 kl TOC \o 1-5 \h \z — ，b = — m m m m k ° — k J、 3 k 』k J、+ k 2】「 9 k — 9 d — Jc ml 得 r x + ax-b0 = 0 〔0 - 〔0 - ex + dO — 0 (3)二自由度可解耦系统的振动模型，其微分方程如下所示: (I) ( X ) / / 7 q / / / / C2 / / 图4 二自由度可解耦系统的振动模型 mxx{ + (q + c2)Xj 一c2x2 + + k2)x{ -k2x2 = 0 m{x2 + c2(x2 -xl) + k2(x2 -xj = 0 可写成矩阵形式 Mx + Cx+ Kx =0 其中，Mm ]0 其中，M m ] 0 K= W , X士叮 tv 2 v 2 设80,即无阻尼情况，则系统可解耦为两种独立的振动模态。 这个方程可由下列步骤来解 用eig函数求出矩阵K-AM的特征值L和特征向量U, U和L满足 u *M*U = I 「人 01 L = U *K*t/二勺 _0 a2_ 在原始方程加+Kx二0两端各作乘以U及在屮间乘以对角矩阵 , X.=U * 1=U x得*M *(/*/*壬 , X. =U * 1 =U x _ Z] 作变量置换Z— 7 也2 z + L*z = 0 这是一个对角矩阵方程，即可把它分成两个方程求解： Zj + 右 Z] = 0 z 9 + A 2 z 2 = 0 这就是说两种振动模态可解耦，令?2 =，①是第i个模式的固有频率 (i 二 1,2)。 由上述的解耦模态中，给出初始条件兀0，兀〃0,化为z°,乙必，即可分别 求出其分量z2 (此处阻尼比为0),再按兀=inv(U ) * z变换为X。 假设位置和速度的初始条件分别为 _[ IT _[ 1T 2x(t) = ^ui([ui]TZ=1X0 = lX0] XQ2\ 9 Xd0 = 2 x(t) = ^ui([ui]T Z=1 1 丁 mx0 cos co^ H [wz] Mxd0 sin (4)编写程序上机调试求解并画图，对三种模型的结果进行比较分析，得出和 应结论。 技术难点： 当阻尼不等于零时，用MATLAB无法求解，需要借助数值方法求解。 如果选取不同的广义坐标，则建立的微分方程的形式也不同，方程的耦 合情况和输出相应也将发生相应变化