汽车平顺建模与分析设计

1、目 录1绪论21.1课题研究的背景和意义 2汽车平顺性研究的背景和意义 31.1.2 汽车平顺性建模的意义 41.2国内外汽车平顺性建模与仿真研究现状 41.2.1 面向结构和面向参数的方法比较 41.2.2 动力学模型的选择 51.2.3 国内汽车平顺性的研究 61.3 MATLAB软件介绍 71.4本文主要研究内容 72路面激励82.1引言82.2随机路面的拟合92.2.1 随机路面拟合的方法 92.2.2 路面激励的拟合92.3路面激励的生成112.4车辆平顺性的评价方法 132.5随机输入平顺性评价指标 143汽车模型建立163.1 模型建立163.2微分方程的建立183.3微分方程的

2、矩阵表达式 194 MATLAB软件下平顺性的计算及分析 214.1利用MATLAB编写改进的欧拉方法后的微分方程的解法 214.2求加速度功率谱密度 244.3平顺性的计算以及评价 264.4各种汽车参数对平顺性的影响 28结论错误!未定义书签。致谢错误！未定义书签。参考文献33毕业设计（论文）原创性声明和使用授权说明原创性声明本人郑重承诺：所呈交的毕业设计（论文），是我个人在指导教 师的指导下进行的研究工作及取得的成果。 尽我所知，除文中特别加 以标注和致谢的地方外，不包含其他人或组织已经发表或公布过的研 究成果，也不包含我为获得及其它教育机构的学位或学历而使用过的材料。对本研究提供过帮助

3、和做出过贡献的个人或集体，均已在文中作了明确的说明并表示了谢意。作者签名：日期：指导教师签名：日期：使用授权说明本人完全了解大学关于收集、保存、使用毕业设计（论文）的规定，即：按照学校要求提交毕业设计（论文）的印刷本和电 子版本；学校有权保存毕业设计（论文）的印刷本和电子版，并提供 目录检索与阅览服务；学校可以采用影印、缩印、数字化或其它复制 手段保存论文；在不以赢利为目的前提下，学校可以公布论文的部分 或全部内容作者签名： 日 期：学位论文原创性声明本人郑重声明：所呈交的论文是本人在导师的指导下独立进行研究所取得的研究成果。除了文中特别加以标注引用的内容外， 本论文 不包含任何其他个人或集体

4、已经发表或撰写的成果作品。 对本文的研 究做出重要贡献的个人和集体，均已在文中以明确方式标明。本人完 全意识到本声明的法律后果由本人承担。作者签名：日期：年月日学位论文版权使用授权书本学位论文作者完全了解学校有关保留、使用学位论文的规定， 同意学校保留并向国家有关部门或机构送交论文的复印件和电子版， 允许论文被查阅和借阅。本人授权大学可以将本学位论文的全部或部分内容编入有关数据库进行检索，可以采用影印、缩印或扫描等复制手段保存和汇编本学位论文。涉密论文按学校规定处理。导师签名:日期： 年 月作者签名：日期：年月日指导教师评阅书指导教师评价：一、撰写（设计）过程1、学生在论文（设计）过程中的治学

5、态度、工作精神优 良 中 及格 不及格2、学生掌握专业知识、技能的扎实程度优良中及格不及格3、学生综合运用所学知识和专业技能分析和解决问题的能力优良中及格不及格4、研究方法的科学性；技术线路的可行性；设计方案的合理性 优 良 中 及格 不及格5、完成毕业论文（设计）期间的出勤情况优良中及格不及格二、论文（设计）质量1、论文（设计）的整体结构是否符合撰写规范？优良中及格不及格2、是否完成指定的论文（设计）任务（包括装订及附件）？优良中及格不及格三、论文（设计）水平1、论文（设计）的理论意义或对解决实际问题的指导意义 优 良 中 及格 不及格2、论文的观念是否有新意？设计是否有创意？优良中及格不及

6、格3、论文（设计说明书）所体现的整体水平优良中及格不及格建议成绩：优良中及格口不及格（在所选等级前的内画“ V）指导教师：（签名）单位：（盖章）年 月日评阅教师评阅书评阅教师评价：一、论文（设计）质量1、论文（设计）的整体结构是否符合撰写规范？优良中及格不及格2、是否完成指定的论文（设计）任务（包括装订及附件）？优良中及格不及格二、论文（设计）水平1、论文（设计）的理论意义或对解决实际问题的指导意义 优 良 中 及格 不及格2、论文的观念是否有新意？设计是否有创意？优良中及格不及格3、论文（设计说明书）所体现的整体水平优良中及格不及格建议成绩：优 良 中 及格 不及格（在所选等级前的内画“ V

7、”）评阅教师：（签名）单位：（盖章）年 月 日教研室（或答辩小组）及教学系意见教研室（或答辩小组）评价：一、答辩过程1、毕业论文（设计）的基本要点和见解的叙述情况优 良 中 及格 不及格2、对答辩问题的反应、理解、表达情况优良中及格不及格3、学生答辩过程中的精神状态优良中及格不及格二、论文（设计）质量1、论文（设计）的整体结构是否符合撰写规范？优良中及格不及格2、是否完成指定的论文（设计）任务（包括装订及附件）？优良中及格不及格三、论文（设计）水平1、论文（设计）的理论意义或对解决实际问题的指导意义优 良 中 及格 不及格2、论文的观念是否有新意？设计是否有创意？优良中及格不及格3、论文（设计

8、说明书）所体现的整体水平优良中及格不及格不及格（签名）日评定成绩：优 良 中 及格（在所选等级前的内画“ V）教研室主任（或答辩小组组长）：年 月教学系意见:系主任：（签名）年 月 日1绪论汽车平顺性是指汽车过程中能保证乘员不致因车身振动而引起不舒适和疲乏 感觉，以及保持运载货物完整无损的性能。汽车平顺性的优劣直接影响到乘员的 舒适性和货物的完整性，并制约着汽车动力性和经济性的发挥及零部件的使用寿 命。随着汽车日新月异的发展，汽车平顺性也被越来越加重视，体现着汽车技术 的提高。1.1课题研究的背景和意义汽车平顺性研究的背景和意义随着生活水平的提高和生活节奏的加快，人们对汽车乘坐舒适性的要求越来

9、 越高，正朝着安全、舒适、环保、节能的方向前进，尤其是汽车乘坐的舒适性更 是人们追求高品质生活的重要体现，汽车平顺性是影响汽车乘坐舒适性的重要原 因。因此，汽车平顺性的研究日益受到国内外的重视。对于汽车本身而言，由于各种性能相互影响，因此，作为汽车的性能指标之 一的平顺性也影响着其它性能的发挥。在过程中，如果平顺性很差，强烈振动产 生的冲击会加速零部件的磨损，降低零部件的疲劳寿命。车辆的振动还是车内噪 声的主要来源之一，汽车的强烈振动还会使车轮跳离地面，影响汽车的动力性、 制动性以及操纵稳定性。为了减小振动产生的冲击，必须放慢车速，但却使运输 效率降低。低速会使燃油燃烧不充分，而使燃油经济性变

10、差，排放性也变差。对 于驾驶员而言，好的平顺性能会让人心情愉快，充分享受驾驶带来的快乐。否则， 汽车的不停颠簸，会使驾驶员疲惫，注意力不集中，易引发交通事故，造成人财 损失。目前调查表明，60%以上的人乘车过程中有过难受、晕车甚至呕吐的经历 汽车平顺性的降低不仅使工作效率降低，还严重影响着人的身心健康，长期处于 不舒适的振动环境中，不仅容易引起疲劳、心慌，还容易引发各种心脏疾病。良 好的平顺性会使长途货车保持货物的完整，会使长途客车的旅行者更加舒适，驱 赶长途跋涉带来的烦恼，会使农用运输车不因低级路面而影响劳作。其它行业由 于用途的不同，对汽车平顺性也提出了各自的要求。因此，研究平顺性的意义重

11、 大，必须不断增加研究深度，拓宽研究领域，满足不同行业的需求。因此，平顺 性优劣直接影响到乘员的舒适性，并波及车辆动力性和经济性的发挥，是车辆在 市场竞争中争夺优势的一项重要性能指标。因而如何最大限度地降低汽车在过程 中所产生的振动，甚至更进一步利用振动来为我所用，是一项十分有价值和意义 的工作。随着国民经济连续多年的高速发展，尤其是国家对基础设施建设投入的逐年 加大，使得载货汽车的生产在近年来呈现了爆发式发展。而汽车由于具有运输效 率高、运输成本低的特点，逐渐成为公路运输的宠儿。我国的汽车市场从进入新 世纪以来，取得了突飞猛进的发展，销售业绩基本上是一路走高，增长迅速。年 均增长幅度高达 7

12、0%左右， 2008 年的销售量是 2000 年的销售量的 6.57 倍，汽车 的生产与开发成为国内载货汽车生产厂家竞争的焦点。汽车的工作条件比较恶劣， 的道路标准低、弯道和坡路多，车辆频繁转向与制动，并长期在满载、振动与冲 击载荷下工作，座椅部位的振动水平是一般客车的 9到 16倍，更为严重的是大多 数商用车驾驶员在这种环境中一天工作时间长，这给驾驶员带来较大疲劳和危险。 同时，平顺性对车辆油耗和车体损伤也有重要影响。因此，提高载重车辆的平顺 性十分必要。但较有实力的汽车制造公司基本上都在欧洲和美国。如奔驰公司的 汽车， 2002年的产量约占世界总产量 15.62%，远远超过其它竞争对手。瑞

13、典的沃 尔沃和斯堪尼亚正在力求在新产品开发方面缩短与奔驰公司产品的距离，以此来 提高自己产品的市场占有份额。目前，我国拥有汽车制造汽车厂家超过 10 家，如 中国第一汽车集团公司、重庆汽车集团公司、中国东风汽车集团公司、陕西汽车 集团公司、北汽福田公司等，各厂家都有自己的特点和能力，但差距各异。所以 提高汽车平顺性已经成为了必然的趋势。1.1.2 汽车平顺性建模的意义 为了保证汽车具有良好的平顺性，汽车性能的预测极其重要。最初是通过试 验评价汽车平顺性的好坏，同时验证汽车建模与仿真结果的正确与否。然而，汽 车平顺性的试验研究固然重要，但是由于采用样车试验的方法，设计周期长，浪 费大量的人力、物

14、力和财力，不能满足现代汽车快速更新、生产周期短和汽车生 产商之间激烈竞争的要求，使得汽车平顺性的建模与仿真研究变得非常重要。汽 车平顺性的建模与仿真研究不受试验条件的限制，可以避免试验中的一些随机因 素的影响，便于比较分析。不仅能够对已有汽车进行评价，提出改进措施。更重 要的是可以预测新产品的平顺性，发现平顺性存在的问题并及时解决，同时还能 优化各参数间的匹配关系，改善平顺性，缩短设计周期，降低生产成本，在设计 初期就能够预测各部件对汽车平顺性的影响。国内外学者在汽车平顺性的建模与 仿真方面进行了大量的研究及其应用，但是，汽车建模与仿真方法等的简单、快 速和精确程度不仅影响着预测结果的真实性，

15、更影响着设计周期的长短、生产成 本的高低。探讨更为简单、快速和精确的建模与仿真方法，对促进汽车平顺性研 究具有重要的意义。1.2 国内外汽车平顺性建模与仿真研究现状1.2.1 面向结构和面向参数的方法比较20世纪 50年代后，仿真技术开始发展，并被逐渐引入到汽车振动研究领域。20 世纪 70 年代，汽车动态仿真技术已在国外得到普及， 产生了不同复杂程度的汽 车型。汽车动态仿真的方法可以分为两大类：多体参数法和集中参数法。但是无 应用哪一种方法进行汽车平顺性分析，都必须建立合理而正确的汽车振动的力模 型。多体参数法，是面向结构的方法，需要给定各部件的详细特征，将汽车的每 一件看作刚体或弹性体，通

16、过各种约束连接建立汽车结构系统振动的拓扑结构， 然后相应的商业化软件，如 ADAM、DADS等进行仿真。商业化的多体参数法的软 件需要使用者有很高的专业水平，同时，应用软件建的汽车模型一般自由度很多， 有些参数难以得到，所以不能从整体上保证系统的准性。而且，复杂的模型在计 算机上求解时间长，且一旦模型出错，很难准确查找问题。集中参数法，是面向参数的方法，不必考虑汽车的具体结构，只要根据汽车 振动分析需要给出描述汽车零部件结构的质量、刚度和阻尼参数，即可进行汽车 平顺性的分析。因此，应用集中参数法建立汽车结构系统振动的力学模型，方法简单。在汽车设计初期，由于不能完全得到汽车结构，应用集中参数法建

17、立汽车 结构系统振动的力学模型，对汽车平顺性的预测与分析是非常有效的。122动力学模型的选择在动力学模型的选择方面，经常使用的基于集中参数法的汽车结构系统振动的力学模型可分为三类：1/4汽车模型、1/2汽车模型和空间整车模型。（1）1/4汽车模型1/4汽车模型认为汽车前后轴上方车身部分的集中质量的运动是相互独立的， 只考虑垂直振动。分析车身垂直振动最简单的模型是单质量系统力学模型，这种 模型旨在将注意力放在悬架设计上，在汽车主动和半主动悬架研究与设计中得到 大量应用。双质量系统力学模型考虑了轮胎弹性和质量，能够反映车轮高频振动 时的动态特性，更接近汽车振动的实际情况，这种模型在汽车的安全性及车

18、轮对 路面的损伤研究方面，能够反映悬架与路面损伤两者之间的关系。考虑了人体和 座椅的三自由度系统力学模型，不仅可以考虑车轮载荷与安全性的关系，还可以 考虑人体加速度与振动舒适性的关系。（2）1/2汽车模型1/2汽车模型（如图1-1所示）考虑了汽车过程中的垂直和俯仰振动，由于 前后车轮对汽车的激励存在时间的滞后关系。因此，在运用1/2汽车模型进行汽车平顺性分析和悬架、轮胎等参数优化时，必须考虑汽车车速的影响。图1-1（3）空间整车模型由于左右车轮对汽车的真实激励并不完全相同，而是具有一定的相关性，使 得汽车产生了侧倾振动。因此，空间整车模型的基本特征是，同时考虑了汽车车 身的垂直、俯仰和侧倾运动

19、，更能全面描述汽车的整体振动特性。国外学者也用 空间整车力学模型进行过汽车平顺性和操纵稳定性研究。由于 1/2 汽车模型考虑了汽车过程中的垂直和俯仰振动， 以及前后车轮对汽车 的激励存在时间的滞后关系。所以此次设计采用的是1/2 汽车模型。1.2.3 国内汽车平顺性的研究从七十年代后期以来，随着随机振动理论、概率论、电子计算机技术在汽车 行业中普及应用，以及一些先进试验测试设备、仪器的引进，我国汽车平顺性的 实验研究工作有了突飞猛进的进展。许多研究所高等院校的研究人员付出了许多 劳动，开展了大量的基础研究工作，制定了一系列近代的试验评价方法。例如八 十年代初，长春汽车研究所、清华大学等首先采用

20、ISO2631 “人体暴露在全振动环 境下的评价指南”国际标准进行汽车在实际道路上平顺性的评价研究。这套标准 主要强调平顺性的不仅与真懂得强度有关，而且与振动的频率和方向有关。再例 如应用数学统计方法成功地设计出了汽车试验跑道。不少研究成果已应用到平顺 性的评价与改进工作中，对我国汽车新产品的开展与老产品的改进起到了重要作 用，大大缩短了我国汽车技术水平在这一领域与国际先进水平的差距。八十年代初，有关改进汽车平顺性的研究工作也取得了相当的进展，如汽车 结构动力学模型及平顺性计算机模拟的研究，汽车座椅振动特性及其改进的研究， 汽车悬架系统的结构分析与改进方面的研究。基中一些研究成果己经应用到了产

21、 品中，使汽车平顺性提到了改善。八十年代对汽车座椅的试验分析工作中，说明 座椅是影响汽车平顺性的重要因素。清华大学、长春汽车研究所进行了许多有关 座椅传递特性、人椅系统动态参数识别、座椅特性与汽车底盘振动特性的合理 匹配等方面的研究。 1994 年重庆大学徐中明等对汽车座椅系统的动力分析进行了 研究。1996年吉林工业大学林逸等对汽车座椅计算机辅助设计方法进行了广泛的 探讨研究，该文通过建立 7 个自由度汽车动力学模型，对汽车座椅平顺性进行了 仿真计算，并对座椅刚度、阻尼进行了优化，提出了汽车座椅参数化绘图软件的 开发方法。1.3 MATLAB 软件介绍由于本文采用的很多方法都要通过 MATL

22、A软件来实现，所以在此介绍MATLAB 软件，方便阅览。MATLA是由MathWork公司与1984年推出的一套计算软件，分为总包和若干个工具箱，可以实现数值分析、优化、统计、偏微分方程数值解、自动控制、信 号处理、图像处理等若干个领域的计算和图形显示功能。它将不同数学分支的算 法以函数的形式分类成库，使用时直接调用这些函数并赋予实际参数就可以解决 问题，快速而且准确。近年来，MATLABt国内的知名度越来越大，并已被广泛的 应用于教学和科研。该软件的特点可以归纳为以下几点：(1) 简单易学MATLAB是一门编程语言，其语法规则与一般的结构化高级编 程语言如 C 语言等大同小异，而且它不需要定

23、义变量和数组，使用更加方便。具 有一般语言基础的用户很快就可以掌握。(2) 代码短小高效由于 MATLAB已经将数学问题的具体算法编成了现成的函 数，用户只要熟悉算法的特点、使用场合、函数的调用格式和参数意义等，使用 调用函数就可以很快的解决问题，而不必花大量的时间纠缠于具体的算法。(3) 计算功能非常强大该软件具有强大的矩阵计算功能，利用一般的符号和 函数就可以对矩阵进行加、减、乘、除的运算以及转置和求逆运算，而且可以处 理稀疏矩阵等特殊的矩阵，非常适合于有限元等大型数值算法的编程。此外，该 软件现有的六十多个工具箱，可以解决数学和工程领域的绝大多数问题。(4) 强大的图形表达功能该软件不仅

24、可以绘制一般的二维、三维图形，如线图、条形图、饼图、散点图、直方图与误差条图等，还可以绘制工程特性较强的 特殊图形，如玫瑰花图、极坐标图、二维、三维等值线图、三维表面图、假彩色 图、二维、三维流线图、三维彩色流锥图、流沙图、流带图、流管图、卷曲图与 切片图等，此外还可以生成快照图和进行动画制作。( 5)可扩展性能可扩展性能是该软件的一大优点，用户可以自己编写 M文件，组成自己的工具箱，以解决本 领域内常见的计算问题。1.4 本文主要研究内容本文以理论力学、统计数学等学科为理论基础，应用 MATLA仿真软件，针对 某商用车，建 9 自由度整车动力学模型并进行平顺性仿真。对于加快新产品开发 进程，

25、降低开发费用，提高产品性能和可靠性具有实际的应用价值。 本文具体工作是：第 1 章，阐述课题的背景及其目的以及轿车、汽车的平顺性研究的意义，针对 汽车平顺性仿真的关键问题：路面激励建模、振动模型建模和仿真方法进行了综 述和分析，提出本文研究的内容。第 2 章 路面模型与评价指标研究。汽车振动的主要激励源来源于路面，路面模 型研究是汽车平顺性研究的关键之一。本文重点研究了汽车的时域路面激励模型。 确定了采用谐波叠加原理模拟路面时域模型，针对平稳随机路面进行研究，用不 同形式的三角级数进行模拟，以离散谱逼近目标随机路面模型。根据 B 级标准路 面谱，模拟再现了路面的功率谱密度，验证了模拟路面与标准

26、路面基本吻合。参 考国标GB7031-86车辆振动输入一一路面平度表示方法,采用MATLABS序编写 路面程序，建立了标准 B级随机路面，参考GB5092-86汽车平顺性脉冲输入试 验方法建立了脉冲输入路面。重点介绍国际标准 02631-1(1997)人体承受全 身振动的评价指南及 GB/T4970-1996汽车平顺性随机输入试验方法，在此基 础上确定本文的评价方法，得到的驾驶员座椅处加权加速度均方根值和加权振级， 然后根据该值与人的主观感觉的对应关系进行乘坐舒适度评价。第 3 章，根据 1/2 汽车模型建立 9 自由度动力学模型，模型中包括了驾驶室、 发动机、人体座椅，能够方便对平顺性进行评

27、价。然后利用牛顿第二定律建立9自由动力学模型的微分方程，为方便以后计算求出刚度矩阵、阻尼矩阵、输入刚 度矩阵并将微分方程用矩阵的形式表示。第 4 章，在建立矩阵形式的微分方程的基础上运用改进的欧拉算法，在 MATLAB 软件里实现对垂直加速度进行积分，然后在 0.5-50HZ 内求加速度均方值功率谱密 度，然后根据平顺性评价方法依次求出各种速度下、各种路面等级下的加权加速 度均方根值以及加权振级，然后在特定车速及某种路面等级下分析各种汽车参数 对平顺性的好坏的影响。第 5 章，全文总结。2 路面激励2.1 引言在进行平顺性分析时必须对引起汽车振动的振源进行分析。引起汽车振动的 振源主要有路面不

28、平度激励、发动机激励、轮胎激励等，其中路面不平度激励是 影响汽车平顺性的最主要激励，因此研究汽车平顺性必须研究路面。路面输入激 励一般分为两类：随机输入激励和脉冲输入激励。随机输入激励是指路面长度方 向上的连续激励，如粗糙不平的路面等，它是现代公路的主要形式；脉冲输入激 励是指在较短时间内的离散事件，并且有较强的强度，如平坦道路上明显的凸包 或凹坑等。2.2 随机路面的拟合2.2.1 随机路面拟合的方法路面不平度的时域仿真有如下方法；将任意一条路面轨迹看作一系列离散的 正弦波叠加而成的三角级数法；将路面高程的随机波动抽象为满足一定条件的白 噪声，经一定假定进行适当变换而拟合出路面随机不平度的线

29、性滤波白噪声；根 据实际路面试凑参数而确定形状函数的过滤泊松过程模型；基于频域功率谱密度 采样的傅里叶变换方法。 以上四种方法均属于功率谱密度分析模型，但前三种方 法是针对路面不平度的数值模拟，所得路面不平度对应的功率谱密度与给定的功 率谱密度相比都存在一定的误差，而第四种方法与给定的功率谱密度一致。在实 际中，利用仪器可测得路面不平度的有限数据的特点，建立了路面不平度的统计 模型：AR和ARMA模型。随着分形理论和小波理论的发展，提出了分形分析模型和 小波分析模型，但两者还都处于研究的初始阶段，对于这两种理论在路面不平度 特性的描述以及使用方法还有待进一步研究。应用路面不平度的时域模型进行汽

30、 车平顺性仿真较多是线性滤波白噪声法、三角级数法和傅里叶变换方法。尽管每 一种方法都给出了理论推导，并且都有广泛的应用，但由于路面不平度的时域模 型均是在统计特征一致的情况下构造的，每种模型只能是路面不平度这一随机过 程的一个样本。同时，上述方法均把路面不平度看成是各态历经的平稳随机过程。2.2.2 路面激励的拟合不管是标准道路谱，还是实测道路谱，其功率谱密度函数（PSD）是路面不平 度的一个统计量。因此，对应于测量范围内某一种确定的路面不平度，其PSD是唯一的；但对于给定的 PSD，其模拟设计的路面不平度并不唯一，也就是说频域 模型和时域模型并非一对一的映射，因此从频域模型所得的路面不平度的

31、时域模型只能看成是满足给定路谱的全部可能的路面不平度中的一个样本函数。为从已 知的路面谱，获得路面激励时域模型随机路面不平度，通用的方法是将路面不 平度定性为平稳的随机过程。大量的测量分析结果表明，路面不平度具有随机、 平稳和各态历经的特性，可以用平稳随机过程理论来分析描述。通常把道路垂直 纵断面与道路表面的交线作为路面不平度的样本，通过样本的数学特征方差或 功率谱密度函数来描述路面。均值为零时，方差可以反映路面不平度大小的总体 情况；功率谱密度函数能够表示路面不平度能量在空间频域的分布，它说明了路 面不平度或者说路面波的结构。从功率谱密度函数不仅能了解路面不平度的结构， 还能反映出路面的总体

32、特征。上文中提到的谐波叠加法拟合不平路面的原理是：设路面高程为平稳的、 遍历的均值为0的Gaussian过程，则可以用不同形式的三角级数进行模拟。本文 以正弦波进行模拟，随机正弦波叠加法采用以离散谱逼近目标随机过程的模型， 是一种离散化数值模拟路面的方法。随机信号可以通过离散 Fourier 分析变换分解 为一系列具有不同频率和幅值的正弦波。谱密度就等于由带宽划分的这些正弦波 幅值的平方。路面的不平对汽车产生的激励具有随机性，其统计特性通常主要采 用路面功率谱密度来描述。按照国标GB7031-86车辆振动输入一一路面平度表示方法，路面不平度位移功率谱密度拟合表达式采用下式：G q(n) G q

33、(n0) (n / n0 ) w式中no是参考空间频率；n是空间频率，它表示波长的倒数，即每米长度中包括几 个波长；Gq(n。)为参考空间频率no下的路面功率谱密度值，称为路面不平度系数， 单位为米的立方；w为频率指数，为双对数坐标上斜线的斜率，它决定路面功 率谱密度的频率结构。已知在空间频率n! n n2内的路面位移谱密度为 Gq( n)，利用平稳随机过程的平均功率的频谱展开性质，路面不平度的方差2为Gq (n)dn将区间ni，n2划分为n个小区间，取每个小区间的中心频率nmid(i)(i=l , 2, 3n)处的谱密度值Gq(nmid i代替Gd(n)在整个小区间内的值，则离散化后近视为n

34、zGq (nmid i )nii 1对应每个小区间，现在要找到具有频率nmid(i)( i=1，2，3. ，n)且其标准差为niGq(“mid i)正弦波函数，这样的正弦波函数可为.2Gq(nmid i) ni sin(2 nmid iX i)将对应于各个小区间的正弦波函数叠加起来，就得到时域路面随机位移输入n q(x). 2Gq(nmid i) ni sin(2 nmid iX i)i 1式中 为(0，2 )上的均匀分布的随机数；X为时域路面的X方向同样将空间频率功率谱密度化为时间功率谱密度后可以得到：nq(t). 2Gq (f mid i) f Sin (2 % it i)i 1于是便可以

35、得到路面垂直位移随时间的变化关系。在GB7031文件里规定，按功率谱密度把路面分成八个等级并规定每种路面等级下 不平度系数Gq(n)的取值范围和集合平均值，如表 2-1所示：表2-1路面等级Gq (n0 )/(10 6m3)n00.1m 1V10 3m)0.011m 11n 2.83m几何平均值几何平均值A163.81B647.61C25615.23D102430.45E409660.90F16384121.80G65536243.61H262144487.222.3路面激励的生成根据某一等级路面不平度系数 Gq(n。)的取值，可计算出一定空间频率范围内 的路面不平度功率谱密度Gq(n)数据曲

36、线，将路面不平度功率谱密度Gq(n)离散化, 作为已知输入带入式，根据谐波叠加原理，应用MATLAB件编制路面谱生成程序， 可以生成这种等级下的时域路面文件。根据上面的理论公式，根据国标GB7031-86车辆振动输入一一路面平度表示 方法中规定的B级路面数据，拟合出B级随机路面，B级随机路面的空间谱如图 2-1所示。图2-1把按照谐波叠加法理论拟合出的 B级随机路面的时间谱，再作谱分析拟合出B级随机路面的频域谱(具体求法将在第四章提到)，将此频域谱与按国标GB7031-86车辆振动输入一一路面平度表示方法提供的B级路面理论上的频域谱进行对比，如图2-2所示所示。从图中可得，拟合路面于标准路面基

37、本吻合。综上所述，按照上面“谐波叠 加法来构建随即路面不平度的时域模型”理论建立的随机露面是准确的。2-.2图2.4车辆平顺性的评价方法汽车平顺性的评价是一个极为复杂的过程，它包括人、车、路三个环节，其 中人是最活跃的因素，当前对汽车平顺性进行评价主要分两类：主观评价和客观 评价。主观评价是人对汽车平顺性最直接的评价方法。主观评价方法主要是根据 有一定经验的试车人员对汽车振动的直观感受进行统计分析并对车辆进行评价。 经过对主观评价的研究分析，设计人员根据经验可以简单地改变如汽车悬架参数 来提高汽车的平顺性。但是，由于车辆的动态特性和人体对振动响应的复杂性， 主观评价只能是对汽车平顺特性的一个比

38、较模糊的描述。此外，由于人与人之间 存在的差异，以及人体自身复杂的心理、生理特性，即使对同样的汽车振动的感 觉也会不一致，这样，就会导致难以对汽车平顺性进行定量、准确的评价分析， 因此需要专门评价人员进行。而客观评价方法由于排除了人的个人差异，以具有量值的概念对汽车平顺性 进行评价，从而可以比较精确合理地评价、分析汽车的平顺性。它主要考虑车辆 的隔振性能，以机械振动的各个物理量（如振幅、频率、加速度等）作为评价指 标，通过测试传递到人体的振动量的大小，来确定影响人体舒适性的程度，以此 来评价汽车的平顺性，因此，这是一种较为合适的评价方法。国际标准ISO2631法用于评价振动得到多数人的认同，我

39、国也参照该标准制定了相应的汽车振动评 价方面的国家标准。2.5随机输入平顺性评价指标如图2-3所示ISO2631-1(1997)35标准规定在进行舒适性评价时，它除了考 虑座椅支撑面处输入点3个方向的线振动，还考虑该点三个方向的角振动，以及 座椅靠背和脚支撑面两个输入点各 3个方向的线振动，共3个输入12个轴向振 动。此标准认为人体对不同频率振动的敏感程度不同，不仅给出了各轴向0.580的频率加权函数，又考虑不同输入点、不同轴向的振动对人体影响的差异，还给 出了各轴向振动的轴加权系数。图2-3ISO2631-1(1997)标准规定，当振动波形峰值系数9(峰值系数是加权加速度 时间历程aw(t)

40、的峰值与加权加速度均方根值aw的比值)时，用基本的评价方法一 加权加速度均方根值来评价振动对人体舒适和健康的影响。由于本次设计只是考虑垂直路面位移对平顺性的影响，所以只计算垂直方向上的加权加速度均方根值。计算垂直加权加速度均方根值aw有两种计算方法1)对纪录的加速度时间历程a (t)，通过相应频率加权函数 w(f)的滤波网络 得到加权加速度时间历程aw(t)，按下式计算加权加速度均方根值。TO2aw 1/ T aw dt0式中：T为振动的分析时间，取100s.其中频率加权函数w(f)用下列的公式表示:0.50.5f/42Wk(f )1412.5/ f12.510.5wd (f)2/ f210.

41、5wc(f)8/f8 f10.5we(f)2412.5802 ；808；801 ；802)对记录的加速度时间历程a ( t)进行频谱分析得到功率谱密度函数Ga ( f)计算 aw。先计算1 /3倍频带加速度均方根谱值:1/2fuiai Ga(f),fli式中ai是中心频率为fci的第i个1/3倍频带加速度均方根值，m/S2 ； fui、fli 分别是第i个1/3倍频带的上下限频率，H乙fd第i个1/3倍频带的中心频率,Hz； Ga( f)为awMai ri加速度自功率谱密度函数，m/s3;然后计算式中W为第i个1/3倍频带的加权系数;为了表示方便，可以将得到的垂直加权加速度均方根值 aw换算成

42、对数形式，用加权振级进行比较。Law 20 lg(aw/a)；式中a0为参考加速度均方根值，a010 6m/s2。为了方便对客观数据进行评价，表2-2给出了加权振级Law和加权加速度均方根 值aw与人的主观感受之间的关系表2-2加权加速度均方根值aw2(m/s )加权振级Law(db)人的主观感受2.0126极不舒服3汽车模型建立3.1模型建立动力学建模是一个提取结构力学特征的过程，需要适当的简化。整车动力学 模型应包括驾驶员、驾驶室、发动机、车架及载荷、前轮簧下质量、后轮簧下质 量。1/2汽车九自由度结构系统不仅考虑了车身、驾驶室和发动机的俯仰，还考虑 了人体、驾驶室和发动机的垂直振动，能够

43、反映载重汽车的平顺性，建立的动力 学模型如图3-1。图中参数意义如表3-1和表3-2所示。LL J |11 1.TFT8 1IT“斗哄览啄tI图3-1.汽车9自由度模型表3-1质量kgM9驾驶员51M7驾驶室730M5发动机1,653M3车架及载荷19,702M前轮簧下775M2后轮簧下2,140惯性矩kgm 2I 8驾驶室380I 6发动机270I 4车架及载荷238,000刚度系数 N/mK9人体座椅系统318.75冷驾驶室前悬架1.96E+4K8r驾驶室后悬架1.96E+4K5发动机前悬置3.14E+6K6发动机后悬置2.45E+6K前悬架3.82E+5K4后悬架1.78E+6Ki前轮胎

44、1.96E+6K2P后轮胎7.84E+6阻尼系数N/(m/s)C9人体座椅系统63.75C驾驶室前悬架7.50E+3C8P驾驶室后悬架7.30E+3C5发动机前悬置5.00E+3C6P发动机后悬置3.64E+3C3前悬架3.18E+4C4后悬架1.17E+5布置尺寸mXcf驾驶室前悬架-1.18XcrP驾驶室后悬架0.64Xc驾驶室质心-0.54Xef发动机前悬置-0.32Xerr发动机后悬置0.65Xe发动机质心0.21Xf前悬架0Xr后悬架6.53Xs车架及载荷质心5.423.2微分方程的建立利用牛顿第二定律可以得到九个微分方程，如下:ID1Z1C3(Z1Z1)k3 (Z1 z1)k1Z1

45、k1q1；m2Z2C4(Z2Zz)k4(Z2z2)k2Z2k2 q2 ；m3Z3k3(z1Z1)C3(z1Z1)1k 4 1(Z2Z2)C4(z2 Z2)k7(Z3 z7)C7(z31Z7) k5 (z4z1 1：9)C5 (Z 4z9)1k6 (z 51Z10 )C6(z51Z 10 )k8 (Z16Z8)C8 (Z16z8)0；打4(k3(z1Z1)C3(z1Z1)11(k4(Z2Z2)C4 (Z2Z2)12 (k7(Z3z7)C7(Z31Z7) 13 (k5(Z1 14Z9JC5(Z14Z:9) 14(k6 (z5乙0)C5 (Z5z10)15(k8 (z6Z8)C8 (Z6Z8)160

46、；m5Z5k5(z9Z4)C5(z9z4)k6(Z10Z5)1C6 (z 10；Z；)0；166(k5(Z9Z4)C5(z9z4)C(k61(z 10 z：5)C6 (Z110 z5) d05m7Z7k7(z7Z3)C7(z7z3)k8(z81 1 1Z6 ) C8 (z 8 Z 6)k9 (Z7Z9)C9(z1 88(k7(z7Z3)C7(z7z3)a(k81 1(Z8Z61)C8 (z81Z 6)0；mgZ9kg(Z9Z7)C9(Z9Z7)0；Z9)0；3.3 微分方程的矩阵表达式为了方便计算，要将方程改进成矩阵的表达式mz c z k zkt QZZ1,Z2,Z3, 4,Z5, 6 ,Z

47、7, 8,Z9z11z3 l1 4;z2 z3 l24;z3z3l31 14; z4 z3 l4 4 ;z5z3 l5 4 ;z6z3 l6 4; z7 z7 a8;z8z7b1 18 ;z9 z5 c 6; z10z5 d 6;m,k , c分别表示质量矩阵，刚度矩阵，阻尼矩阵，具体表示为mdiag (m1 , m2 ,m3 , I 4,m5,I 6,m7,I8,m9)k11 k12 k13 k14 k15k16k17k18k19k22 k23 k24 k25k26k27k28k29k33 k34 k35k36k37k38k39k44 k45k46k47k48k49K=k55k56k57k5

48、8k59k66k67k68k69对称k77k78k79k88k89k99式中 K11=K1+K3;K13=-K3;K14=l 1*K3;K22=K2+K4;K23=-K4;K24=-l 2*K4K33=K3+K4+K5+K6+K7+K8;K34=-l 1*K3+l 2*K4-l 3*K7-l 4*K5-l 5*K6-l 6*K8;K35=-K5-K6;K36=c*K5-d*K 6;K37=-K7-K8;K38=a*K7-b*K 8;K44 = l 1A2*K3+l 22*K4 + | 32*K7+| M2*K5 + | 52*K6+| 62*K8；K45 = l 4*K5 + | 5*K6；K46=-C*l 4*K5 + d*| 5*K6；K47 = | 3*K7 + | 6*K8；K48=-a*| 3*K7 + b*| 6*K8；K55我+；K56=-c*K