资源描述：
《汽车五自由度建模.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、汽车振动大作业..一、汽车悬架系统振动模型汽车是一个复杂的振动系统，在振动分析的建模过程当中，要根据所分析的问题对汽车进行简化，建立相应的模型。现在考虑汽车车身悬架的五自由度模型，如下图1所示，该模型主要考虑左右车辙的不平度差异和较小的轮胎阻尼而得到的，该模型中主要有车身的垂直、俯仰两个自由度和前后车轴质量两个垂直自由度，汽车座椅一个垂直自由度，系统共五个自由度，其中车身质量的垂直、俯仰两个自由度的振动对系统平顺性的影响较大，假设车身是具有垂直和俯仰两个自由度的刚体，其车身的质量和转动惯量分别为：，前后车
2、轮质量、悬架参数和轮胎刚度的符合前加入了分别表示前(front)和后(rear)的下标“f”和“r”，如图1示：图1五自由度汽车悬架系统图1中：表示前轮转动位移自由度；表示车体垂直位移自由度；表示后轮转动位移自由度；俯仰转动位移自由度；表示驾驶员座椅垂向自由度；表示驾驶员座椅质量；表示车体质量；表示前轮质量；表示后轮质量；表示座椅弹簧刚度;悬架弹簧刚度；表示座椅弹簧阻尼；表示悬架弹簧阻尼；表示车身质心至前轴距离；车身质心至后轴距离，分别为前后轮随机激励力。..二、运动微分方程由图1可得到下述理论值：(1)
3、系统的动能为：(1-2)(2)系统的势能为：(1-3)(3)系统阻尼耗散的能量：(1-4)由拉格朗日运动方程：可得到多自由度的运动微分方程：式中：..表一汽车结构参数汽车结构参数数值—驾驶员座椅质量—车体质量—右前、左前轮胎质量—左后、右后轮胎质量—转动惯量—座椅弹簧刚度—右前、左前悬架弹簧刚度—左后、右后悬架弹簧刚度—座椅弹簧阻尼—悬架弹簧阻尼—车身质心至前轴距离—车身质心至后轴距离-座椅到质心距离取汽车结构参数如表一所示，则可求得系统的质量矩阵，阻尼矩阵，刚度矩阵分别为：..由特征方程求得固有频率与振
4、型。根据系统的模型方程，用MATLAB得到系统的固有频率与振型，固有频率如表2所示，固有振型如图2所示.表2各阶固有频率数值阶数12345单位（Hz）387.94.388.61872.81856.6固有振型为：..图2系统固有振型三、自由振动分析当系统的初始条件确定时，可以求得系统的自由振动，假设初始条件为：..初始位移：初始速度：1.无阻尼自由振动....2.有阻尼自由振动....3.频响函数阻尼系统的频响函数矩阵为：..将式3-2左乘，右乘得频响函数矩阵的模态展开式：计算了有阻尼的频响特性，如下图；.
5、...相关程序：1.固有阵型：clearallclcm=[650000070800000106000000800000080];k=[23071-20292-432600-2307162689-39882-20292-19326-4326-398825570915219289890-202921521914905200-19326289890148086];[v,d]=eig(k,m)[omeg,w_order]=sort(sqrt(diag(d)));df=omeg./(2*pi)plot(v(:,w_
6、order(1)),'-rs','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','g','MarkerSize',7,title'')%第一阶振型%plot(v(:,w_order(2)),'-rs','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','g','MarkerSize',7)%第二阶振型..%plot(v(:,w_order(3)),'-rs','LineWidth',2,'
7、MarkerEdgeColor','k','MarkerFaceColor','g','MarkerSize',7)%第三阶振型%plot(v(:,w_order(4)),'-rs','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','g','MarkerSize',7)%第四阶振型%plot(v(:,w_order(5)),'-rs','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','g
8、','MarkerSize',7)%第五阶振型1.无阻尼自由振动：clcM=[650000070800000106000000800000080];K=[23071-20292-432600-2307162689-39882-20292-19326-4326-398825570915219289890-202921521914905200-19326289890148086];[E,F]=eig(K,M);W=diag(sqr
