Design of shock absorber matching to optimal damping of vehicle suspension
-
摘要: 为了使设计减振器对车辆具有最佳减振效果, 利用悬架最佳阻尼比, 对减振器最佳阻尼系数进行了研究, 建立了减振器最佳速度特性数学模型, 提出了减振器阀系参数设计优化方法, 对设计减振器进行了特性试验和整车振动试验, 并与原车载减振器性能进行了对比。计算结果表明: 减振器特性试验值与最佳阻尼匹配要求值的最大偏差为9%, 而且, 在低频范围内, 设计减振器的整车振动传递函数幅值明显低于原车载减振器的幅值, 有效遏制了簧下质量在13Hz附近的共振, 因此, 减振器速度特性模型和阀系参数优化设计方法是正确的。Abstract: In order to design optimum shock absorber for vehicle, its optimal damping coefficients were studied according to the optimization damping ratio of suspension, a mathematic model of its velocity characteristic was built, an optimal design method of its valves parameters was put forward, its performance and whole vehicle vibration tests were carried out, and its performances were compared with that of original shock absorber. Computation result shows that the maximum warp between testing values and optimal requiring values is 9%, the transfer function values of whole vehicle vibration with designed shock absorbers are lower than that with original shock absorbers in the scope of low frequency, the resonance hereabout 13 Hz is depressed, so the model and method are feasible.
-
0. 引言
汽车减振器是悬架系统中重要部件, 能有效地衰减簧上和簧下质量的振动, 提高车辆行驶平顺性和操纵稳定性[1-2]。在车辆中, 应用最多的是筒式液压减振器, 其阻尼力主要是由减振器阀系参数决定的[3-4]。对于不同车辆类型或同一车辆的前、后轮都需要不同阻尼特性的减振器。
减振器阀系参数设计一直是困扰减振器设计和开发的关键问题[5]。目前, 国内、外对基于车辆参数的减振器阀系参数设计还没有准确、可靠的设计方法[6], 大都是凭经验首先确定一个设计参数, 然后经过反复试验和修改, 最后才确定出参数值。由于减振器多个阀系参数之间相互影响, 对于每个设计参数需要经多次试验才可最终确定其值。虽然曾有学者对基于车辆参数的减振器阀系参数设计进行了大量的研究, 但很少通过建立最佳阻尼匹配减振器特性数学模型, 对阀系参数进行优化设计[7], 大都是利用《机械设计手册》所提供的阀片最大挠度计算公式, 对阀系参数进行近似设计, 参数设计值不可靠, 因此, 目前传统的减振器阀系参数设计方法, 不能满足汽车减振器设计和生产的需要。
本文对车辆悬架系统最佳阻尼匹配减振器速度特性进行研究, 建立了车辆悬架最佳阻尼匹配减振器速度特性数学模型, 对减振器阀系参数进行优化设计, 并对设计减振器进行了特性试验与整车振动特性验证。
1. 筒式减振器工作原理
汽车筒式减振器大都采用双筒式结构, 有4个阀, 分别是复原阀、补偿阀、压缩阀和流通阀, 其中, 复原阀和压缩阀对减振器特性起决定作用, 其结构原理, 见图 1。
减振器处于复原行程时, 复原运动速度为Vf, 活塞缸筒上腔的油液和储油腔的一部分油液, 分别流经复原阀和补偿阀, 产生复原节流压力。减振器处于压缩行程时, 压缩运动速度为Vy, 活塞缸筒下腔中的一部分油液, 分别经过流通阀和压缩阀产生压缩节流压力。在某一速度下, 节流压力与相应面积的乘积, 即为减振器在该速度下的阻尼力。
2. 最佳阻尼匹配减振器阻尼系数
2.1 悬架系统最佳阻尼系数
设单轮总质量为mt, 单轮簧下质量为md, 则单轮簧上质量mu为
mu=mt-md
依照1/4单自由度振动模型, 根据悬架系统最佳阻尼比, 可得悬架系统最佳阻尼系数[8]为
Cs=4πξf0mu (1)
式中: ξ为悬架最佳阻尼比; f0为悬架固有频率。
2.2 减振器最佳阻尼系数
2.2.1 首次开阀最佳阻尼系数
根据悬架系统杠杆比i, 可求得减振器在首次开阀时的阻尼系数为
C1=4πξf0mu/i2 (2)
2.2.2 二次开阀最佳阻尼系数
减振器阻尼是非线性的, 一般将减振器速度特性分段线性化, 并将减振器首次开阀前的速度特性直线斜率k1与二次开阀前的速度特性直线斜率k2的比值, 定义为减振器平安比η。其中, 减振器首次开阀前的速度特性曲线斜率k1为减振器开阀前的阻尼系数C1, 因此, 减振器平安比η可表示为
η=k1/k2=C1/k2 (3)
k2为
k2=C1/η (4)
根据阻尼力与阻尼系数和速度之间的关系, 可得减振器二次开阀阻尼系数C2为
C2=F2V2=F1+k2(V2-V1)V2 (5)
式中: V1、V2分别为减振器首次开阀速度和二次开阀速度; F1和F2分别为减振器在首次和二次开阀时的阻尼力。
将式(4)代入式(5), 得
C2=C1η+C1(1-1η)V1V2 (6)
3. 最佳阻尼减振器特性建模
3.1 复原行程速度特性
当减振器复原行程达到首次开阀速度V1f时, 减振器首次开阀复原阻尼力可表示为
F1f=C1fV1f=C1V1f (7)C1f=C1
式中: C1f为减振器复原首次开阀阻尼系数; V1f为减振器复原行程首次开阀速度。
将式(2)代入式(7), 可得
F1f=4πξf0muV1f/i2 (8)
当减振器复原运动达到二次开阀速度V2f时, 减振器二次开阀阻尼力为
F2f=C2fV2f=C2V2f (9)
C2f=C2
式中: C2f为减振器复原二次开阀阻尼系数; V2f为减振器复原二次开阀速度。
将式(6)代入式(9), 则减振器复原二次开阀阻尼力可表示为
F2f=4πξf0mui2[V1f+1ηf(V2f-V1f)] (10)ηf=η
式中: ηf为减振器复原行程平安比。
3.2 压缩行程速度特性
在减振器压缩行程中, 首次开阀速度和二次开阀速度分别为V1y和V2y。为了增加车辆的乘坐舒适性, 一般情况下, 压缩行程的首次开阀速度V1y小于复原行程的首次开阀速度V1f, 即V1y < V1f。
同理, 当减振器压缩行程达到首次开阀速度V1y时, 减振器首次开阀压缩阻尼力可表示为
F1y=C1yV1y (11)C1y=C1
式中: C1y为减振器压缩行程首次开阀阻尼系数; V1y为减振器压缩行程首次开阀速度。
将式(2)代入式(11), 可得
F1y=4πξf0muV1y/i2 (12)
由式(8)、(12)可知复原和压缩的首次开阀阻尼系数等于减振器首次开阀最佳阻尼系数
C1y=C1f=C1 (13)
考虑车辆行驶的平稳性以及压缩行程中弹簧力的影响, 一般要求减振器压缩行程二次开阀阻尼力小于复原行程二次开阀阻尼力, 即
F2y=F2f/k (2≤k≤3) (14)
式中: k为减振器复原行程二次开阀阻尼力与压缩行程二次开阀阻尼力的比值。
将式(10)代入式(14), 可得
F2y=4πξf0muki2[V1f+1η(V2f-V1f)] (15)
设压缩行程二次开阀阻尼系数为C2y, 且二次开阀速度V2y等于复原行程的二次开阀速度V2f, 则压缩二次开阀阻尼力可表示为
F2y=C2yV2y=C2yV2f (16)
将式(14)代入式(16), 可得减振器压缩行程二次开阀时的阻尼系数C2y为
C2y=C2f/k (17)
可知, 压缩行程二次开阀阻尼系数仅为复原行程的1/k, 因此, 利用减振器复原行程速度特性, 以及减振器复原行程二次开阀阻尼力与压缩行程二次开阀阻尼力的比值k, 可确定减振器压缩行程首次开阀后的速度特性。
利用压缩行程的首次开阀和二次开阀速度, 以及对应的开阀阻尼力, 可以求得减振器压缩行程的平安比为
ηy=kC1y/C2f
结合式(13), 可得压缩行程平安比为
ηy=kη=kηf (18)
即压缩行程平安比是复原行程的k倍。
3.3 最佳阻尼匹配减振器特性曲线
由以上分析可知, 利用车辆参数以及复原行程和压缩行程的开阀速度, 可求得减振器首次开阀阻尼力F1和二次开阀阻尼力F2, 由此可得设计减振器要求的分段线性速度特性曲线, 见图 2。
4. 最佳阻尼匹配减振器设计
筒式减振器由复原阀、压缩阀、流通阀和补偿阀组成, 复原阀和压缩阀的阀系参数决定和影响减振器的速度特性。
4.1 常通节流孔优化设计
根据车辆悬架最佳阻尼减振器在区间[0, V1]的速度特性, 可建立常通节流孔面积优化设计目标函数WA(V)为[9]
WA(V)=∫V10FA(V)dV-∫V10F(V)dV (19)
式中: FA(V)为利用速度设计点V设计常通节流孔面积A所对应减振器开阀前的阻尼力特性, 在减振器结构给定情况下, 可通过仿真得到; F(V)为最佳阻尼匹配减振器速度特性, 见图 2。
当优化设计目标函数WA(V)达到极值时, 所对应的速度设计点为常通节流孔面积最优速度设计点, 利用速度点所设计的常通节流孔面积为复原阀常通节流孔面积A的优化设计值。
4.2 节流阀片厚度优化设计
常通节流孔面积设计值确定之后, 同理可根据减振器在区间[V1, V2]的速度特性, 建立节流阀片厚度优化设计目标函数Wh(V)为
Wh(V)=∫V2V1Fh(V)dV-∫V2V1F(V)dV (20)
式中: Fh(V)为利用速度设计点V设计节流阀片厚度h所对应减振器首次开阀后的阻尼力特性, 在减振器结构给定情况下, Fh(V)可通过仿真得到。
同理, 当阀片厚度优化设计目标函数Wh(V)达到极值时, 所对应的阀片设计厚度为节流阀片厚度h的优化设计值。
4.3 其他参数设计
当阀片厚度h设计确定之后, 根据减振器首次开阀和二次开阀时节流阀片所受节流压力, 利用节流阀片变形系数[10], 可分别对阀片预变形量和最大限位间隙进行设计。减振器节流阀片预变形量的设计值可表示为
fk0=Gkp1h3 (21)
式中: fk0为阀片预变形量, 即阀片在阀口半径位置处的初始变形量; Gk为节流阀片在阀口半径位置的变形系数; p1为首次开阀时节流阀片所受压力。
减振器节流阀片的最大限位间隙设计值为
δmax=Gkp2h3-fk0 (22)
式中: δmax为阀片最大限位间隙, 即减振器达到二次开阀速度时, 节流阀片在阀口半径位置处的开度; p2为二次开阀时节流阀片所受压力。
5. 最佳阻尼匹配减振器设计实例
5.1 车辆参数
某轿车单轮的平均质量mt为450.5 kg, 单轮簧下质量md为35.0 kg, 簧上质量mu为415.5 kg, 悬架系统最佳阻尼比ξ为0.3, 固有频率f0为1.35 Hz; 复原行程首次开阀和二次开阀速度分别为0.3和1.0 m·s-1, 平安比为1.6;压缩行程首次开阀和二次开阀速度分别为0.1和1.0 m·s-1; 减振器复原行程二次开阀阻尼力是压缩行程的2.4倍。
5.2 最佳阻尼匹配减振器特性
根据车辆最佳阻尼匹配减振器速度特性的确定方法, 可得到该车辆最佳阻尼匹配减振器在不同速度下所要求的阻尼力数值, 见表 1。
表 1 速度特性值Table 1. Characteristic values of speeds速度V/(m·s-1) -1.0 -0.1 0 0.3 1.0 阻尼力F/N -656.3 -209.5 0 629.8 1 543.7 5.3 优化设计值
根据单轮簧上质量, 利用最佳阻尼比和平安比系数, 对车辆减振器进行了曲线拟合优化设计, 设计得到的复原阀和压缩阀优化设计具体参数, 见表 2。
表 2 参数优化设计值Table 2. Optimum design values of parameters设计参数 A/mm2 h/mm fk0/mm δmax/mm 复原阀 0.90 0.277 6 0.034 0.06 压缩阀 1.33 0.164 0 0.150 0.21 6. 减振器特性试验与验证
6.1 减振器特性试验
利用日本鹭宝公司生产的多功能液压振动试验台, 对优化设计减振器施加一定幅值和一定频率正弦谐波激励信号, 进行特性试验, 利用拉、压力和位移传感器测量减振器阻尼力和位移, 并利用分析软件对所测得的数据进行分析和处理。试验测得的减振器特性数值见表 3, 速度特性曲线见图 3。
表 3 速度特性试验值Table 3. Testing values of speed characteristic速度V/(m·s-1) -1.0 -0.1 0 0.3 1.0 阻尼力F/N -690 -190 0 608 1 539 可知, 减振器特性试验值与车辆悬架系统最佳阻尼匹配所要求的减振器特性值基本吻合, 最大偏差仅为9%。
6.2 整车振动特性验证
利用液压振动试验台, 对车辆单轮施加一定频率和幅值的激励信号, 在不同激振频率和幅值下, 利用布置在各测量部位的加速度计, 测量整个车辆不同部位的振动加速度, 测点布置见图 4。
为了使设计减振器与原车减振器进行特性对比验证, 在同一位置对设计减振器和原车减振器进行对比试验。对单轮施加激振信号, 并对该轮上的振动信号进行传递函数分析。原车减振器和设计减振器在右前轮位置3的传递函数曲线, 见图 5。
由图 5可知, 在低频范围内, 设计减振器的整车振动传递函数幅值, 明显低于原车载减振器的整车振动传递函数的幅值, 并且有效遏制了簧下质量(车轮)在13 Hz附近的谐波共振, 比原车载减振器具有更好的减振效果。
7. 结语
通过车辆悬架最佳阻尼匹配减振器速度特性的研究, 以及设计减振器的特性试验和整车振动特性试验, 可知根据车辆参数和平安比, 利用车辆悬架最佳阻尼比, 可建立减振器速度特性数学模型, 得到分段线性速度特性曲线。根据最佳阻尼匹配减振器速度特性, 利用阀片变形解析计算式和油液节流压力与流量之间关系, 建立阀系参数优化设计目标函数, 对阀系参数进行优化设计。试验分析结果表明减振器阀系参数设计值准确、可靠, 减振器特性试验值与最佳阻尼匹配所要求的特性值吻合, 与车辆悬架达到最佳阻尼匹配。
-
表 1 速度特性值
Table 1. Characteristic values of speeds
速度V/(m·s-1) -1.0 -0.1 0 0.3 1.0 阻尼力F/N -656.3 -209.5 0 629.8 1 543.7 表 2 参数优化设计值
Table 2. Optimum design values of parameters
设计参数 A/mm2 h/mm fk0/mm δmax/mm 复原阀 0.90 0.277 6 0.034 0.06 压缩阀 1.33 0.164 0 0.150 0.21 表 3 速度特性试验值
Table 3. Testing values of speed characteristic
速度V/(m·s-1) -1.0 -0.1 0 0.3 1.0 阻尼力F/N -690 -190 0 608 1 539 -
[1] 余强, 马建. 半主动悬架系统对汽车侧翻稳定性的改善[J]. 交通运输工程学报, 2005, 5(1): 53-56. doi: 10.3321/j.issn:1671-1637.2005.01.013Yu Qiang, Ma Jian. Stabilityi mprovement of vehicle s antio-verturn withsemi-active suspension system[J]. Journal ofTraffic and Transportation Engineering, 2005, 5(1): 53-56. (in Chinese) doi: 10.3321/j.issn:1671-1637.2005.01.013 [2] 陈涛, 魏朗. 人-车-路互动模式虚拟仿真系统[J]. 长安大学学报: 自然科学版, 2007, 27(1): 67-71. doi: 10.3321/j.issn:1671-8879.2007.01.016Chen Tao, Wei Lang. Driver-vehicle-road virtual simulation withinteraction mode[J]. Journal of Chang an University: Natural Science Edition, 2007, 27(1): 67-71. (in Chinese) doi: 10.3321/j.issn:1671-8879.2007.01.016 [3] 赵成, 胡增荣, 陈大跃. 半主动悬架的滑模变结构控制[J]. 中国公路学报, 2007, 20(3): 109-114. doi: 10.3321/j.issn:1001-7372.2007.03.020Zhao Cheng, Hu Zeng-rong, Chen Da-yue. Sliding mode varying structure control for semi-active suspension[J]. China Journal of Highway and Transport, 2007, 20(3): 109-114. (in Chinese) doi: 10.3321/j.issn:1001-7372.2007.03.020 [4] 吕彭民, 和丽梅, 尤晋闽. 基于舒适性和轮胎动载的车辆悬架参数优化[J]. 中国公路学报, 2007, 20(1): 112-117. doi: 10.3321/j.issn:1001-7372.2007.01.021Lu Peng-min, He Li-mei, You Jin-min. Optimization of vehicle suspension parameters based on comfort and tyre dynamic load[J]. China Journal of Highway and Transport, 2007, 20(1): 112-117. (in Chinese) doi: 10.3321/j.issn:1001-7372.2007.01.021 [5] 吕振华, 李世民, 刘目珍. 筒式液阻减振器工作特性的实验研究[J]. 汽车工程, 2005, 27(2): 203-208. doi: 10.3321/j.issn:1000-680X.2005.02.019Lu Zhen-hua, Li Shi-min, Liu Mu-zhen. An experi mental study on the characteristics of telescopic hydraulic shock absorber[J]. Automotive Engineering, 2005, 27(2): 203-208. (in Chinese) doi: 10.3321/j.issn:1000-680X.2005.02.019 [6] 周长城, 顾亮. 筒式减振器叠加节流阀片开度与特性试验[J]. 机械工程学报, 2007, 43(6): 210-215. doi: 10.3321/j.issn:0577-6686.2007.06.038Zhou Chang-cheng, Gu Liang. Superposition throttle-slices opening size and characteristic test of telescope damper[J]. Chinese Journal of Mechanical Engineering, 2007, 43(6): 210-215. (in Chinese) doi: 10.3321/j.issn:0577-6686.2007.06.038 [7] 陈龙, 江浩斌, 周孔亢, 等. 半主动悬架系统设计及控制[J]. 机械工程学报, 2005, 41(5): 137-141. https://www.cnki.com.cn/Article/CJFDTOTAL-JXXB200505028.htmChen Long, Jiang Hao-bin, Zhou Kong-kang, et al. Control and design for semi-active suspension[J]. Chinese Journal of Mechanical Engineering, 2005, 41(5): 137-141. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-JXXB200505028.htm [8] Zhao He-ping, Huang Hong-cheng, Xi Gang, et al. Dynamic characteristics of vehicle suspension with nonlinear springs[J]. Journal of Mechanical Strength, 2001, 23(2): 165-167. [9] 周长城. 汽车减振器阀系参数解析计算与特性综合仿真[D]. 北京: 北京理工大学, 2006. [10] 周长城, 顾亮, 王丽. 节流阀片弯曲变形与变形系数[J]. 北京理工大学学报, 2006, 26(7): 581-584. https://www.cnki.com.cn/Article/CJFDTOTAL-BJLG200607004.htmZhou Chang-cheng, Gu Liang, Wang Li. Bending deformation and coefficient of throttle-slice[J]. Transactions of Beijing Institute of Technology, 2006, 26(7): 581-584. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-BJLG200607004.htm -