2. 安徽理工大学 车辆工程系, 安徽 淮南 232001;
3. 中国第一汽车股份有限公司 汽车振动噪声与安全控制综合技术国家重点实验室, 吉林 长春 130011
2. Department of Vehicle Engineering, Anhui University of Science & Technology, Huainan 232001, China;
3. State Key Laboratory of Comprehensive Technology on Automobile Vibration and Noise & Safety Control, China FAW Group Corporation, Changchun 130011, China
作为主减速器、差速器和半轴的装配基体, 驱动桥壳是汽车的主要承载件和传力件, 其使用寿命直接影响汽车的有效使用寿命.
在实际生产及应用中, 驱动桥壳材料性能的波动、零部件制造误差、载荷的不稳定都会影响桥壳的性能和使用寿命.稳健性设计方法主要包括Taguchi方法、6σ和双响应面法等[1-4].轻量化的方法包括轻量化材料应用、轻量化工艺技术和轻量化设计理念.轻量化材料和工艺技术侧重于生产过程, 轻量化设计理念将轻量化思想融入到设计阶段.目前, 针对驱动桥壳的性能稳健性优化的研究相对较少, 本文通过运用Taguchi方法和综合评价优化方法, 探讨使驱动桥壳具有较高稳健性和较轻质量的最优结构参数组合, 验证Taguchi方法和综合评价方法在驱动桥壳设计开发的实用性.
1 驱动桥壳结构分析 1.1 驱动桥壳疲劳寿命分析驱动桥壳由高强度钢板冲焊而成, 驱动桥壳两侧轴头、法兰盘、板簧座和减速器壳通过焊接固定在桥壳上.
焊接结构疲劳寿命主要由焊缝疲劳强度决定[5].在实验室工况下, 桥壳失效主要出现在驱动桥壳的轴头连接位置和桥壳与法兰盘连接位置, 如图 1所示.
利用ABAQUS软件对实验室工况下的驱动桥壳进行有限元分析计算.在建立好的桥壳有限元模型上施加约束定义, 约束条件为通过建立半轴套管圆环正下方150 mm的约束处的多点约束(multi-point constraint, MPC)单元, 载荷定义为单位载荷, 均匀分配施加在2个板簧座表面单元节点上.
将有限元分析结果文件导入FEMFAT软件, 输入桥壳各个材料的应力-寿命曲线.循环加载方式与设计的台架试验加载方式相同, 在板簧座处施加频率为11 Hz的正弦波加载, 载荷峰值为382.5 kN, 最小载荷取最大载荷的0.1倍.
由桥壳疲劳分析计算结果可知, 驱动桥壳在实验室疲劳寿命台架试验工况下, 最薄弱的部位是法兰盘和桥壳本体的焊接部位以及轴头与桥壳本体的焊接部位, B50下的损伤率计算结果分别是6.543e-6和1.066e-6, 与实验室疲劳寿命台架试验结果一致.
1.2 桥壳性能影响因素分析在生产制造过程中, 桥壳零部件的材料性能波动、桥壳关键部位尺寸的偏差都会影响到桥壳的疲劳寿命.桥壳零部件的材料性能主要包括材料的密度、杨氏模量和泊松比.重卡驱动桥壳壳体厚度是不等厚设计, 桥壳本体结构分为长轴和短轴两部分, 长轴和短轴两部分均为上下对称设计.长短轴根据截面厚度的变化可分为8个区域.桥壳关键部位尺寸包含桥壳长轴和短轴上的8个截面的厚度尺寸, 除了截面1只有1个径向厚度尺寸外, 其余7个截面都具有桥壳水平方向上的厚度、竖直方向上的厚度和圆角半径3个尺寸.此外, 施加在驱动桥壳上载荷的大小也会影响到驱动桥壳的使用寿命.
在驱动桥壳多目标优化模型中, 优化对象为重型卡车驱动桥桥壳, 优化目标有两个, 分别是①桥壳危险部位应力波动降低到最小; ②桥壳的质量降低到最轻, 且质量的稳健性最好.
优化模型包含7个设计变量,这些设计变量分为可控设计变量和不可控设计变量.可控设计变量指的是在设计生产制造过程中可以控制的设计变量, 又称为控制因素, 包括桥壳各截面的水平方向厚度、各截面竖直方向厚度、截面圆角半径和载荷.不可控设计变量是指在设计生产制造过程中存在且无法控制或者消除的变量, 又称为噪声因素, 包括桥壳材料密度、桥壳材料的弹性模量和泊松比.
优化模型以桥壳危险部位B50下的规定寿命值为约束.结合疲劳寿命试验结果和仿真计算结果, 桥壳的危险部位主要有桥壳轴头与桥壳本体焊接部位和法兰盘与桥壳本体焊接等部位.
优化模型具有两种输出:一种是桥壳危险部位的应力; 一种是桥壳的质量.当设计变量发生波动, 危险部位应力也会随之发生变化, 桥壳的疲劳寿命也会发生变化.因此,通过提高输出应力的稳定性来提高桥壳疲劳性能的稳健性.
3 桥壳优化方法 3.1 田口方法运用田口方法中的望目特性和望小特性两种优化策略对驱动桥壳的性能进行优化.
3.1.1 望目特性优化策略运用望目特性对桥壳的稳健性进行优化.不同的设计变量组合导致桥壳危险部位应力值波动.桥壳危险部位应力值的均方差Sm由式(1)计算得到,
(1) |
(2) |
桥壳危险部位应力值平方和ST可以由式(3)计算得到,
(3) |
桥壳稳健性能的信噪比S/N为
(4) |
桥壳稳健性能的敏感度S为
(5) |
运用望小特性对桥壳质量进行优化, 桥壳质量的信噪比如式(6)所示[8-9]:
(6) |
综合评价法用下式表示:
(7) |
式中, U(H1, H2, …, Hn)是系统综合性能评价值.若n个性能评价值之间是相互独立的, 系统综合性能评价值是n个单一性能评价值的线性函数, 可由下式表示[10]:
(8) |
式中, Wi是单一目标性能参数的权重系数, 各权重系数之和为1.
(9) |
式中:Hi是性能参数值; Pi是单个目标性能参数的评价值, 取值范围是0~9.当单个目标性能参数值仅达到最低标准时, Pi取0;当单个目标性能参数值达到最优时, Pi取9.
当Pi=9, 系数A可表示为
(10) |
式中:H*是性能参数最优值; Hi'是性能参数的最低限值.性能参数的权重系数根据产品的性质用途来确定,
(11) |
因此, 系统综合性能评价值函数为
(12) |
根据优化模型的约束条件确定控制因素的取值范围.为了避免桥壳出现应力集中, 以及满足桥壳的生产制造和装配需求, 在优化设计中要求优化前后桥壳截面厚度变化趋势保持一致.每一个控制因素取三个水平的值, 截面3的尺寸参数见表 2.其他截面的尺寸根据截面3的尺寸成比例变化.噪声因素也取三个水平的值, 见表 3.
为了避免桥壳局部位置应力集中导致应力值不准确, 取板簧座处的应力值作为优化系统的输出应力值.
建立正交试验表, 内表为不同控制因素组合的试验组排列, 外表是各试验组在不同噪声因素下的输出应力值和质量值, 见表 4.
运用田口方法计算得到桥壳应力稳健性和轻量化两种单目标的优化结果, 见表 5.
在单一目标优化过程中, 取信噪比最大的控制因素的参数水平, 得到最佳控制因素组合.控制因素的组合为TV3, TH1, FR2, F3时, 桥壳的应力信噪比最高, 桥壳的疲劳性能稳健性能最好.控制因素的组合为TV3, TH3, FR3, F1时, 桥壳的质量最轻, 且质量稳健性最好, 见表 6.
根据驱动桥壳的使用条件和设计要求, 桥壳疲劳性能的稳健性和桥壳的轻量化性能的权重系数均为0.5.运用综合评价方法计算桥壳综合性能的评价值, 相关参数见表 7, 桥壳综合性能的信噪比计算结果见表 8.
控制因素在各个水平下的均值和信噪比值见表 8.多目标下桥壳优化得到的最佳控制因素组合为TV3, TH3, FR2, F3.优化后, 桥壳应力的信噪比增加3.17 dB, 质量的信噪比增加1.21 dB, 质量减轻12 kg, 见表 9.桥壳的疲劳性能稳健性得到了提高, 同时也达到了轻量化设计的目的.
对优化后的桥壳进行疲劳寿命计算, 桥壳本体与法兰盘焊接部位的损伤率为8.143e-6, 桥壳本体与轴头焊接部位的损伤率为3.012e-6.计算结果表明, 虽然优化后的桥壳刚度有所下降, 但疲劳寿命符合要求, 充分地发挥了材料的性能.
5 结论1) 分析影响桥壳疲劳性能稳健性和桥壳质量的因素, 建立桥壳参数化设计模型, 利用正交试验表生成结构参数不同的试验组, 运用田口方法对桥壳的疲劳性能稳健性和桥壳的质量进行单目标优化.
2) 在单目标优化结果的基础上, 对桥壳进行了多目标优化.计算结果表明, 优化后桥壳疲劳性能满足使用要求, 危险部位应力信噪比提高了3.17 dB, 疲劳性能稳定性能得到提高, 桥壳质量减轻12 kg.应用田口方法和综合评价方法节约了桥壳材料, 降低了生产和运营成本.
[1] |
刘为, 薛克敏, 李萍, 等.
汽车驱动桥壳的有限元分析和优化[J]. 汽车工程, 2012, 34(6): 523–527.
( Liu Wei, Xue Ke-min, Li Ping, et al. FE analysis and optimization of vehicle drive axle housing[J]. Automotive Engineering, 2012, 34(6): 523–527. DOI:10.3969/j.issn.1000-680X.2012.06.010 ) |
[2] |
谢晖, 刘行, 洪健程.
基于Taguchi方法的车门结构稳健性优化[J]. 中国机械工程, 2013, 24(12): 1676–1681.
( Xie Hui, Liu Xing, Hong Jian-cheng. Robust optimization of automotive door structure based on Taguchi method[J]. China Mechanical Engineering, 2013, 24(12): 1676–1681. DOI:10.3969/j.issn.1004-132X.2013.12.023 ) |
[3] |
Motta R D, Afonso S M B.
An efficient procedure for structural reliability-based robust design optimization[J]. Structural and Multidisciplinary Optimization, 2016, 54(3): 511–530.
DOI:10.1007/s00158-016-1418-1 |
[4] |
Chiang H L, Yang C B, Hsu C Y.
Combining Taguchi signal-to-noise ratio and grey relational analysis into a multi-objective optimal model for milling inconel 718 superalloy[J]. Journal of the Chinese Society of Mechanical Engineers, 2016, 37(6): 625–633.
|
[5] |
Bhattacharya B, Ellingwood B.
Continuum damage mechanics analysis of fatigue crack initiation[J]. International Journal of Fatigue, 1998, 20(9): 631–639.
DOI:10.1016/S0142-1123(98)00032-2 |
[6] |
Sun G, Fang J, Tian X, et al.
Discrete robust optimization algorithm based on Taguchi method for structural crashworthiness design[J]. Expert Systems with Applications, 2015, 42(9): 4482–4492.
DOI:10.1016/j.eswa.2014.12.054 |
[7] |
Akgul E, Kizilkaya A E, Akay D, et al.
Moora-based Taguchi optimization for selection of murata vortex spinner machine parameters[J]. Tekstil Ve Konfeksiyon, 2016, 26(3): 249–255.
|
[8] |
Abbasi M, Ghafari-Nazari A, Reddy S, et al.
A new approach for optimizing automotive crashworthiness:concurrent usage of ANFIS and Taguchi method[J]. Structural & Multidisciplinary Optimization, 2014, 49(3): 485–499.
|
[9] |
Wilson V H, Udayakumar.
Optimization of diesel engine parameters using Taguchi method and design of evolution[J]. Journal of the Brazilian Society of Mechanical Sciences & Engineering, 2012, 34(4): 423–428.
|
[10] |
Dao T P, Huang S C.
Robust design for a flexible bearing with 1-DOF translation using the Taguchi method and the utility concept[J]. Journal of Mechanical Science & Technology, 2015, 29(8): 3309–3320.
|