资源描述:
振动与冲击 第38 卷第23 期JOURNAL OF VIBRATION AND SHOCKVol.38 No.23 2019 基金项目2018 年度交通运输战略规划政策项目计划 2018- 22- 3 ; 2019 年中央高校基本科研业务费 自然科学类 300102229108 收稿日期2019 -04 -11修改稿收到日期2019 -07 -08 第一作者 李彬 男, 博士, 讲师, 1982 年生 基于试验验证的商用车驾驶室悬置参数优化研究 李彬1,曹阳2,运伟国3 1. 长安大学 汽车学院, 西安710064; 2. 吉林大学 交通学院, 长春 130012; 3. 浙江吉利新能源商用车集团有限公司, 杭州310052 摘要针对目前中国商用车驾驶室舒适性问题, 以整车刚柔耦合多体动力学理论基础, 采用 ADAMS 软件建立某典 型商用车驾驶室及其悬置系统动力学模型, 并对该模型进行了准确性验证。设计了 DOE 优化试验, 经过仿真分析, 对驾驶 室悬置的刚度和阻尼进行优化, 并对试验车辆按照优化参数进行改进。最后, 通过道路试验表明, 优化后的商用车行驶平顺 性较优化前提升了9.6, 改善了驾驶室的隔振性能。 关键词驾驶室悬置系统; 试验; 刚柔耦合模型; DOE 分析; ADAMS 中图分类号U491.254文献标志码ADOI 10.13465/j. cnki. jvs.2019.23.019 Optimizing of commercial vehicle cab suspension paraments by means of experiments LI Bin1,CAO Yang2,YUN Weiguo3 1. School of Automobile,Chang’ an University,Xi’ an 710064,China; 2. Transportation College,Jilin University, Changchun 130012,China; 3. Zhejiang Geely New Energy Commercial Vehicle Group Co. ,LTD,Hangzhou 310052,China Abstract Aiming at the comfort problem of commercial vehicle cabs,the dynamic model of cab and its mounting system of a typical commercial vehicle was established by using ADAMS software, based on the rigid- flexible coupling multi- body dynamics theory. The accuracy of the model was verified. The DOE optimization test was designed,the stiffness and damping of the cab mount are optimized by simulation analysis,and the test vehicle was improved according to the optimization parameters. Finally,the road test shows that the ride comfort of the optimized commercial vehicle is 9. 6 higher than that before optimization,and the vibration isolation perance of the cab is improved. Key wordscab mounting system;test;rigid- flexible coupling model;DOE analysis;ADAMS 随着生活水平的不断提高和科学技术的不断进步以 及公路运输在我国经济发展中重要性的日益增加, 人们 不仅对乘用车的平顺性有较高要求, 而且对商用车平顺 性的要求也越来越高 [ 1 ]。据不完全统计从事专业驾驶 员工作的人, 患骨骼肌肉方面的疾病的概率要高很多 [ 2 ] 。 因此, 驾驶室悬置隔振技术近年来在国内外受到高度重 视, 开展了大量理论和试验研究, 用来改善商用车驾驶室 平顺性。在传统的设计中, 商用车驾驶室直接安装在车 架上, 车辆行驶的振动直接传递到驾驶室, 驾驶室的平顺 性较差, 同时还可能影响行驶安全性 [ 3 ]。而驾驶室悬置 隔振技术是通过弹簧阻尼元件组成的悬置系统将驾驶室 悬置在车架上, 因而能极大的提高驾驶室的平顺性。国 内对于商用车驾驶室悬置隔振系统的应用和研究起步较 晚, 相关的理论分析和应用研究均相对滞后, 特别是综合 利用多体系统动力学和连续弹性理论在整车范围内对基 于刚柔多体参数化仿真模型的驾驶室悬置隔振仿真研究 还很少 [ 4- 5 ]。国内大多数企业主要以模仿设计和经验设 计为主, 通过参照国外实际样车来进行悬置隔振系统的 设计。但驾驶室悬置系统的隔振效果不仅取决于其结构 形式、 几何位置以及弹性元件的刚度、 阻尼等参数的选 择, 而且与底盘悬架系统性能参数的匹配有重要关系, 还 与车辆主要承受载荷的部件如车架等的柔性有较大关 系 [ 6- 8 ]。国内企业在进行模仿设计时由于对国内外车辆 使用情况以及车辆总体结构和结构参数的差异分析不 足, 因而常常使得悬置隔振系统性能参数选取不当, 或与 主悬架性能参数匹配不好, 或者车架设计不好, 导致隔振 效果较差 1整车刚柔耦合多体动力学理论基础 根据 ISO 26312003 标准, 车辆平顺性研究的主要 频率范围为0.5 ~80 Hz, 其中人体对0.5 ~12.5 Hz 内的 振动较为敏感。以往工程师在进行驾驶室悬置匹配与设 计过程中, 一般采用多刚体动力学的研究方法 [ 9 ]。即将 整车部件中除弹性、 阻尼元件之外的其他部件均视为刚 ChaoXing 体, 这种简化作为驾驶室悬置的初步选型在一定程度上 是合理的, 但是在实际车辆结构中, 车架、 车桥和车轮等 主要零件均是具有弹性的 [ 10 ]。根据商用车模态分布情 况, 对于全悬浮式驾驶室悬置结构, 其刚体模态频率分布 在1 ~8 Hz 之间, 发动机和车桥的模态在 80 Hz 之上, 远 远高于驾驶室悬置系统刚体模态频率, 而车架的第一阶 固有频率一般在 10 Hz 左右 与车架的长度有直接关 系 , 根据模态理论, 任意激励的动态响应频率可视为系 统各阶固有模态频率的线性组合, 一般来说, 固有模态频 率与响应频率相隔越近则该阶模态对响应的贡献越大。 因此可以判定发动机和车桥的弹性对整车平顺性的影响 很小, 而车架弹性的影响是较大的 [ 11- 13 ]。 在本文研究中, 假设系统中车架、 及驾驶室悬置支 架为柔性体; 驾驶室、 驾驶员座椅、 车桥等均为刚性体; 悬 架系统中的弹簧和阻尼力为力学单元, 导向机构视为 约束。 根据随机振动理论, 一个具有 m 个输入 Xi t i 1, 2, , m 和 n 个输出 Yk t k 1, 2, , n 的常参数线 性系统可用如图 1 所示的系统表示, 其中 Xi ω 和 Yk ω 分别是Xi t 和Yk t 的傅里叶变换。Hki ω 为第 i 个输入和第 k 个输出之间的频率响应函数。 经过傅里叶变换, 并将输入输出写成矩阵形式, 有 X ω [ X1 ω X2 ω Xm ω ] T 1 Y ω [ Y1 ω Y2 ω Yn ω ] T 2 则系统的频率响应矩阵 H ω为 H ω H11 ω H12 ω H1m ω H21 ω H22 ω H2m ω Hn1 ω Hn2 ω Hnm ω 若设系统 m 个输入的自谱与互谱构成的 m n 阶输入功 率谱矩阵 SX ω 为 SX ω SX1X1 ω SX1X2 ω SX1Xn ω SX2X1 ω SX2X2 ω SX2Xn ω SXmX1 ω SXmX2 ω SXmXn ω n 个输出的自谱与互谱构成的 n n 阶输出功率谱矩阵 SY ω 为 SY ω SY1Y1 ω SY1Y2 ω SY1Yn ω SY2Y1 ω SY2Y2 ω SY2Yn ω SYnY1 ω SYnY2 ω SYnYn ω 则输出功率谱与输入功率谱的关系可写成 SY ω H* ω SX ω HT ω 3 其中H* ω 为H ω 的共轭矩阵, HT ω 为H ω 的转置矩阵。 根据随机振动理论, 设整车系统在 l 个车轮处激励 相互独立, 其相应的激励和功率谱分别为 Xii t 和 Sii f i 1, 2, , l, 那么在仅考虑座椅处垂直方向加速度时 其输出功率谱可表示为 Szz f H* f Sxx f HT f∑ l i 1 Hiz f2Sii f 式中 Szz f 为整车的输入功率谱矩阵, 该矩阵为对角 阵, 可以表示为 Sxx f S11 S22 Sl-1l-1 S ll H f 为系统的传递函数矩阵, 可以表示为 H f [ H1z f H2z f H1l f ] 4 相应的驾驶员座椅垂直方向加速度均方根值, 表 示为 σz ∫ ∞ 0 Szz f d 槡 f 5 2驾驶室悬置系统 ADAMS 建模与验证 2.1刚柔耦合整车模型的建立 根据企业提供的 CAD 模型, 应用多体动力学软件 ADAMS/CAR 模块建立整车的虚拟样机, 整车模型包括 驾驶室及悬置系统、 车架、 前后悬架系统、 前后桥等子系 统。首先建立车架、 驾驶室、 驾驶室悬置系统、 前后悬架 等的模板,这其中包括构件的拓扑结构、 运动副、 衬套以 及输入、 输出通讯器; 其次, 根据模板建立子系统; 最后由 各个子系统组装成整车模型。整车各硬点坐标均来自 CAD 模型, 对关键硬点相对位置进行了实车测量校核。 驾驶室质量参惯性参数及各级悬置、 衬套的刚度及阻尼 参数均来自于台架试验结果, 驾驶室质量及转动惯量测 试如图 1 所示, 各参数见表 1 所示, 虚拟样机如图 2 所示。 2.2模型准确性验证 在进行整车模型的振动仿真分析之前, 需要对整车 模型进行验证及确认所建立模型的准确性, 本次验证包 531第23 期李彬等基于试验验证的商用车驾驶室悬置参数优化研究 ChaoXing 图1驾驶室质量及转动惯量测试 Fig.1Cab quality and rotational inertia test 图2刚柔耦合多体动力学 Adams 模型 Fig.2Rigid and Flexible Multi- Body Adams Model 表1驾驶室质量参数 Tab.1Quality Parameters of Cab 参数数值 驾驶室质量/kg 892 驾驶室惯性矩/ kgm2 Ixx 1024.18 Iyy 808.11 Izz 1015.51 Ixy -26.76 Ixz -77.50 Iyz -9.32 驾驶室前悬置刚度/ Nmm -1 26 驾驶室后悬置刚度/ Nmm -1 34 驾驶室前悬置阻尼/ Nmms -1 11.3 驾驶室后悬置阻尼/ Nmms -1 11.8 括两个方面的内容 1自由度与静平衡验证。该验证在模型建立完成 后即可进行, 目的是检验模型是否欠约束或是过约束。 另外样机模型能否静平衡是能否顺利进行仿真计算的首 要条件。系统静平衡的成功表明所建样机模型的约束和 力学参数可以支撑起整车系统, 从而保证整车系统正常 的仿真计算。 2载荷谱输入输出验证。本整车模型采用六立柱 试验台进行驱动, 通过试验台给汽车轴端加载位移信号, 该位移信号是通过对试验采集的轴端振动加速度信号积 分得到。同时获取模型中驾驶室座椅导轨位置的响应加 速度与试验采集的相同位置的振动加速度进行对比, 如 果加速度均方根值误差控制在 15之内即认为模型准 确, 可以进行下一步的优化工作。通过该方法不仅可以 验证模型的准确性还可以检验各级悬置及衬套刚度试验 结果的准确性。 3优化分析与结果 3.1DOE 优化试验设计 采用虚拟 DOE 正交试验技术对驾驶室悬置进行系 统参数优化, 以驾驶室悬置前后刚度, 前后阻尼作为优化 设计因数, K1为前悬置弹簧刚度; K2为后悬置弹簧刚度; C1为前悬置弹簧压缩阻尼; C2为后悬置弹簧压缩阻尼。 每个因数分为三个水平。以驾驶员座椅导轨位置垂向加 权加速度均方根值作为输出, 计算影响驾驶员乘坐舒适 性的每个因素的最佳水平。 表2不同因素的水平值 Tab.2Level value of different factors 因素 K1/ Nmm -1 K2/ Nmm -1 C1/ Nmms -1 C2/ Nmms -1 水平113 178.4758.85 水平219.5 25.514.12514.75 水平332.5 42.516.9517.7 优化时试验台的激励信号是选择车辆在 B 级路面 上以70 km/h 匀速行驶时通过测试得到的各轴头位移信 号, 该信号同样是通过对加速度信号积分得到的, 并通过 频谱分析进一步得到位移信号的自功率谱。 3.2优化结果 经过仿真分析, 对每种方案驾驶员座椅导轨位置的 加速度均方根值进行计算, 最终挑选出 7 种较为合理的 组合方案, 见表3 所示。 表3优化结果 Tab.3Optimized results 因素 K1/ N mm -1 K2/ N mm -1 C1/ Nmm s -1 C2/ Nmm s -1 加权均 方根值/ ms -2 1131716.9514.1250.345 9 2131716.9517.70.348 9 313178.47517.70.351 0 4131714.12517.70.355 4 532.51716.958.850.356 5 619.51714.12514.1250.357 7 732.51716.9514.1250.358 9 通过将上述 7 个方案分别在整车下通过虚拟 B 级 路面进行路试。计算时分别以 60 km/h、 70 km/h、 80 km/h 的车速进行匀速直线行驶, 检测驾驶员座椅导轨处 加速度, 并通过式 5 对加速度数据分别进行加权均方 根值计算。7 组方案中方案 5、 6、 7 对应的加权加速度均 631振 动 与 冲 击2019 年第38 卷 ChaoXing 方根值较小。计算结果如表4 所示。优化方案较原方案 性能提升百分比见表5 所示。 表4三种较优方案计算结果 Tab.4Three preferred optimized resultsm/s2 车速/ kmh -1 原方案方案5方案6方案7 601.049 10.787 90.806 90.793 5 701.375 71.200 41.141 01.125 0 801.337 81.189 41.235 01.154 7 表5三种较优方案性能提升率 Tab.5The improvement rate of three preferred optimized results 车速/ kmh -1 方案5 方案6 方案7 60252324 70131718 8011814 平均161619 根据以上分析方案 7 为最优方案, 是考虑交互影响 的前提下对垂直振动衰减最好的选择, 使驾驶员座椅导 轨处加权加速度均方根值改善了19。 4试验验证 按照上优化方案对车辆进行了改进, 对改进后的车 辆进行试验验证。根据 ISO 26312003 的规定, 主要以 驾驶室座椅导轨位置的加权加速度均方根值 aw来评价 平顺性。aw的计算公式为 aw ∑ 20 j 1 ω2j∫ fuj fij Ga f d[]f 1 2 6 式中 fij、 fuj分别为1/3 倍频带的中心频带为 f j的下、 上限 频率; Ga f 为等带宽的加速度自功率谱密度; a ω为单轴 向加权加速度均方根值; ωj为第 j 个 1/3 倍频带的加权 系数。 测试系统为 LMS SCM09 多通道振动测试系统, 加速 度传感器选择 PCB ICP 型加速度传感器, 所用的 Test. lab 软件可以记录加速度时间历程, 计算出测试过程中的 aω, 部分测试设备与传感器布置位置如图3 所示。 本测试在西安绕城高速进行, 测试工况均为匀速行 驶工况, 车速依次从 40 km/h 按照 10 km/h 递增至 80 km/h, 测试结果如图4 所示。 从图4 得知, 优化后的驾驶室悬置对在各车速下整 车平顺性均有所提升, 平均提升 9. 6。其中, 60 km/h 到70 km/h 车速时加权加速度均方根值下降, 主要原因 是60 km/h 时轮胎不平衡引起的一阶激励与驾驶室刚体 模态重复。由于试验道路与仿真所加载的 B 级路面激 励的不一致, 其优化提升幅度低于仿真结果, 但在一定程 度上证明了仿真优化的有效性和科学性。 aLMS 数据采集仪 bPCB 三向加速传感器 cLMS 数据采集仪 d传感器布置位置 图3座椅导轨位置传感器布置图 Fig.3Sensor layout on Seat rail position 图4改进前后驾驶员座椅导轨处加速度均方根值对比 Fig. 4Comparison of the root mean square value vibration acceptationofdriverseatrailsbeforeand after improvement 5结论 1以多体系统动力学理论为基础, 以机械动力学 仿真软件 ADAMS 为工具, 建立了基于整车的全浮式驾 驶室悬置系统刚柔耦合的参数化仿真模型, 利用 DOE 试 验方法以驾驶员座椅导轨位置的加速度均方根值为目标 进行了优化。通过实车道路试验证明优化后整车行驶平 顺性较优化前提升了9. 6, 说明该优化方案是正确的, 科学的。 2改优化方案充分考虑了车架刚度及驾驶室悬置 支架刚度等对驾驶室悬置匹配的影响, 把其作为弹性体 建模, 该方法可以作为车架等刚度设计目标的制定依据, 另外还可以利用该模型对车架刚度, 驾驶室悬置支架刚 度对车辆平顺性的影响进行研究。 参 考 文 献 [1] 高菲. 基于骨肌生物力学的驾驶员疲劳机理与舒适性评价 研究[ D] . 长春 吉林大学, 2017. [2] PORTER J M,GYI D E. The prence of musculoskeletal troubles among car drivers [J] .Occupational Medicine Lond , 2002, 52 1 4 - 12. [3] DE TEMMERMAN J,DEPREZ K,ANTHONIS J,et al. 731第23 期李彬等基于试验验证的商用车驾驶室悬置参数优化研究 ChaoXing Conceptualcabsuspensionsystemforaself- propelled agricultural m- achine, part1 developmentofalinear mathematical model[J] . Biosystems Engineering,2004,89 4 409- 416. [4] 王楷焱,史文库,杨昌海,等. 基于 ADAMS 的商用车驾驶 室悬置系统振动模态和传递特性[ J] . 吉林大学学报 工学 版 , 2010, 40 2 330- 334. WANG Kaiyan, SHIWenku, YANGChanghai, etal. Commercial vehicle cab suspension system vibration modes and transmission characteristics by means of ADAMS[ J] . Journal of Jilin University Engineering and Technology Edition , 2010, 40 2 330- 334. [5]苑雷. 商用车驾驶室空气悬置匹配优化[ D] . 长春 吉林大 学, 2008. [6] SOLIMAN A M A,GAZALY N M,KADRY F S. Parameters affecting truck ride comfort[ C] . SAE Papers, 2014. [7] 王国林,李凯强,杨建, 等. 商用车驾驶室悬置仿真与隔振 性能优化[ J] . 汽车工程, 2017, 39 9 1081- 1086. WANG Guolin,LI Kaiqiang,YANG Jian,et al. Simulation and vibration isolation optimization on the cab mount of a commercial vehicle [J] . Automotive Engineering,2017,39 9 1081- 1086. [8] WARD C H. The application of a new cab mounting to address cab shake on the 2003 Chevrolet Kodiak and GMC TopKick [ C] . SAE Paper, 2002. [9]WANG Tao, WANGLiangmo, WANGChenzhi, etal. Development in the acquisition of vehicle loads integrated with a rigid and flexible multi- body model [ J] . SAE International Journal of Commercial Vehicle, 2014, 7 1 163- 169. [ 10] FENG Chunliang,HE Feng,FENG Xusheng. Model analysis for a heavy truck cab based on ANSYS[ J] . Machinery Design & Manufacture, 2013 4 73- 74. [ 11] LIU Jingyi,DONG Qiangqiang,YANG Jiansen. Commercial vehicle cab fatigue load decomposition and validation[J] . Machinery Design Manufacture, 2017 8 187- 193. [ 12] 王正华,喻凡,庄德军. 汽车座椅舒适性的主观和客观评 价研究[ J] . 汽车工程, 2006, 28 9 817- 819. WANG Zhenghua,YU Fan,ZHUANG Dejun. A Study on objective and subjective uations for ride comfort of vehicle seats[ J] . Automotive Engineering, 2006, 28 9 817- 819. [ 13] WANG Tao,WANG Liangmo,LI Tan,et al. Fatigue failure assessment based on loading reproduction for commercial vehicle cab[ J] . Journal of Huazhong University of Science and Technology Natural Science Edition , 2017, 45 5 61- 66. [ 12] 范政武,王铁,王永红, 等. 基于改进多目标粒子群算法的 商用车悬架系统优化[ J] . 振动与冲击,2018,37 15 216- 224. FAN Zhengwu, WANGTie, WANGYonghong, etal. Commercial vehicle suspension system optimization based on improved multi- objective particle swarm optimization algorithm [ J] .Journal of Vibration and Shock,2018,37 15 檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿 216- 224. 上接第109 页 [ 14] SEIREG A,HOUSER D R. uation of dynamic factors for spur and helical gears [J] .Journal of Engineering for Industry, 1970, 92 504. [ 15] 方宗德. 修形斜齿轮的轮齿接触分析[ J] . 航空动力学报, 1997, 12 3 247- 250. FANG Zongde. Tooth contact analysis of modified helical gears [ J] . Journal of Aerospace Power, 1997, 12 3 247- 250. [ 16] 方宗德. 修形斜齿轮的承载接触分析[ J] . 航空动力学报, 1997, 12 3 252- 254. FANG Zongde.Loaded tooth contact analysis of modified helical gears[ J] . Journal of Aerospace Power, 1997, 12 3 252- 254. [ 17] 铁摩辛柯,古地尔,徐芝纶. 弹性理论Theory of elasticity [ M] . 北京高等教育出版社, 2013. [ 18] AHMED F,DEB K. Multi- objective optimal path planning using elitist non- dominated sorting genetic algorithms[ J] . Soft Computing, 2013, 17 7 1283- 1299. 831振 动 与 冲 击2019 年第38 卷 ChaoXing
展开阅读全文