摘 要: 为了使悬架减振器在改善汽车行驶平顺性的同时完成振动能量的回收,首先以某品牌减振器为原型设计一种实用新型液压式馈能减振器的原理及结构,依据工作原理构建该减振器的油液流动中流量与压降之间的关系理论模型;然后在AMESim中建立等效参数化仿真模型分析所设计的液压式馈能减振器的外特性及能量回收特性能否满足实际工作需求;最后在Simulink中进行液压式馈能悬架总成AMESim模型与路面时域输入模型的联合仿真,分析此减振器装车后在实际路面激励下对车辆平顺性的影响及能量回收效果. 分析结果表明:此液压式馈能减振器压缩/复原行程阻尼力符合减振器阻尼力允差的国标要求;示功图形状饱满无明显畸变,体现出良好的外特性,满足减振器的设计要求;馈能特性符合最初设计的只在压缩行程回收能量的思想;馈能特性曲线呈现出显著的峰值特性并受到减振器高频响应特性的影响;悬架平顺性满足行驶要求,实际行驶中具有一定的能量回收潜力. 本文所设计的液压式馈能减振器基本能达到预期目标,对节能减排有借鉴价值.
关键词: 液压式馈能减振器;示功特性;速度特性;馈能特性;悬架平顺性
在能源压力日益增加的今天,学者、企业及科研机构从未停止对汽车行驶时能量回收系统方面的研究. 能量吸收及转化形式包括液压能、电能、热能等,将其贮存起来供其他车辆部件或系统应用[1]. 近年来,在悬架振动能量回收系统方面,许多学者研究出多种类型的馈能减振器. 1990年SATOH M等[2]提出了液压式馈能减振器的初始设计思想;2006年陈士安等[3]设计了一款双腔式液压馈能减振器并建立了其力学特性研究模型,然而关于馈能效率的计算方法并没有进行具体阐述;2011年SHIMOKAWA S等[4]构建了两个独立的液压回路将减振器的油液导出,并在各回路中串联蓄能器,调控蓄能器内的油液压力以衰减振动;ZUO团队[5-7]将减振器活塞的直线往复运动通过齿轮齿条机构转换为永磁发电机的旋转运动,以此原理设计制造出了液电馈能式减振器;徐琳、方志刚等[8-10]针对液电馈能减振器进行深入研究,但其还存在复原行程阻尼力低于压缩行程阻尼力、液压回路沿程损失较高、液压整流桥液压传递效率较低等问题.
当下对于馈能减振器的设计方案为液电式和电磁式两种,对液压式的研究极少. 液压式馈能减振器的优点:1)汽车上以液压能为动力的系统较多,便于应用;2)能量转化次数较少,效率高;3)对原有悬架的调整较小. 因此,本文针对液压式馈能减振器进行研究,并对其外特性及能量回收特性进行分析.
1 液压式馈能减振器原理设计
选用某减振器为设计原型,设计液压式馈能减振器原理及结构如图1所示.

(a) 原理图 (b) 结构示意图
1—工作缸; 2—流通阀; 3—活塞; 4—储油缸; 5—补偿阀; 6—活塞杆; 7—复原阀; 8—常通孔; 9—中间缸; 10—高压蓄能器; 11—压缩阀; 12—单向阀; 13—中间缸出油孔; 14—储油缸出油孔; 15—低压蓄能器
图1 液压式馈能减振器原理及结构示意图
Fig.1 Principle diagram and structure diagram of hydraulic energy regenerative damper
活塞3将工作缸1分成上下两腔,传统减振器原有的工作缸1外增加了中间缸9,中间缸9及储油缸4都有出油孔以将油液导出到油管. 因为减振器正常工作时,复原行程中为求快速减振阻尼力应较大,可供回收的能量较少,故设计此液压式馈能减振器只在压缩行程中回收振动能量. 流通阀2、补偿阀5和复原阀7的布置方式与传统减振器相同. 与传统双筒式液压减振器差别之处在于,为在压缩行程中回收能量,将原本安装在底阀处的压缩阀11外置,与高压蓄能器10相连.
减振器在压缩行程中,活塞向下移动,工作缸下腔压力增大,在压力作用下油液经流通阀2进入上腔. 因为活塞杆的存在,上腔增加的容积会小于下腔减小的容积,所以部分油液会从工作缸1下腔上的常通孔8流入中间缸9,经中间缸出油孔13、单向阀12、压缩阀11最终流入高压蓄能器10储存起来,供转向、制动等其他需要液压能的系统使用,以达到回收振动能量的目标. 同时用油管将低压蓄能器15与储油缸出油孔14相连,低压蓄能器15的作用是增补流入高压蓄能器的液压油.
减振器本体外观尺寸没有变化,结构上的改动在于引出了高/低压油管,并需要在车辆底盘合适空间处安装液压阀和管路,以及在发动机舱内布置高/低压蓄能器. 故改动很小,不影响原有悬架的安装和工作空间.
2 液压式馈能减振器仿真模型建立
根据压强作用原理,减振器的阻尼力为

式中:P1为工作缸上腔压力,P2为工作缸下腔压力,Pair为大气压力,Ah为活塞截面积,Ag为活塞杆截面积,f为摩擦力.
本文设计的液压式馈能减振器只在压缩行程回收能量,故复原行程阻尼力与传统双筒式液压减振器基本一致. 而在压缩行程中,压缩阀开启之前,油路如图2中实线部分所示;压缩阀开阀后,油路增加了图2中的虚线部分. 所以压缩行程阻尼力与传统双筒式液压减振器差别较大,本文对油液流经各阀时流量与压降之间的关系进行推导.

图2 压缩行程压缩阀开阀前后油路
Fig.2 The oil route before and after compression valve opening in the compression process
依据流量与压强的串、并联关系及流量计算方法,可知压缩行程工作缸下腔油液流入上腔的过程中,流量与压力差的关系为

式中:Qyh为活塞外圈孔的流量,Qyl为流通阀节流缝隙的流量,QH为流经活塞缝隙的流量,vy为压缩行程活塞相对运动速度,ΔP1为油液流经活塞体时产生的压差,ΔPyh为油液流经活塞外圈孔的压降,ΔPyl为油液流经流通阀节流缝隙的压降.
根据伯努利方程及厚壁小孔流量公式,减振器运动速度较小时,油液流经工作缸常通孔的流量与压降之间的关系为

式中:Qyg为流经工作缸上常通孔的流量,Cv厚壁小孔流量系数,一般取0.82,Ayg为工作缸上常通孔的等效节流面积,ρ为油液密度,ΔPyg为油液流经工作缸上常通孔的压降.
当减振器速度达到压缩阀开启的程度时,油液流经压缩阀的流量与压降之间的关系为

式中:Qyy为流经压缩阀的流量,Ayy为压缩阀等效节流面积,ΔPyy为油液流经压缩阀的压降.
在AMESim中建立液压式馈能悬架的等效参数化仿真模型如图3所示,右侧部分为液压式馈能减振器的仿真模型.

A—橡胶衬套; B—活塞杆与导向座之间泄露; C—摩擦力; D—活塞与工作缸之间泄露; LTF—流通阀; FYF—复原阀; YSF—压缩阀; BCF—补偿阀; DXF—单向阀; CYQ—储油腔; HSYK—活塞圆孔; YSGK—活塞缸常通孔; GYXNQ—高压蓄能器; DYXNQ—低压蓄能器; LTFK—流通阀片缺口非圆孔; YSFK—压缩阀片缺口非圆孔; HSJX—活塞矩形孔; BCFK—底阀矩形孔
图3 液压式馈能悬架的仿真模型
Fig.3 Simulation model of hydraulic energy regenerative suspension
馈能减振器仿真模型中,在储油腔的端口连接了压力传感器,以获得其中的油液压力作为低压蓄能器的开关信号. 因为压缩行程中减振器内的工作油会被压入高压蓄能器并被其他部件消耗,储油腔内的油液压力会持续减小. 当检测到储油腔中的液体压力低于某一阈值时,低压蓄能器就会开启,向储油缸中补充消耗掉的工作油.
蓄能器选用欧洲标准隔膜式蓄能器,如图4所示. 充入氮气体积0.7 L,具体参数如表1所示.

图4 隔膜式蓄能器
Fig.4 Diaphragm accumulator
表1 蓄能器参数
Tab.1 Accumulator parameters

3 外特性仿真分析
减振器的外特性是指减振器伴随(相对)运动的位移或速度,与相应产生的阻尼力之间的关系,一般又体现为示功特性和速度特性[9].
依据减振器原型,设计此液压式馈能减振器的工作缸缸径为30 mm,活塞杆直径为12.5 mm,减振器工作行程为50 mm. 因为汽车在实际正常行驶时,所受激励源的幅值及频率都是随机的,一般在仿真时选择正弦位移为输入信号. 正弦激励的位移与速度满足

式中:A0为所加载位移的幅值,f为谐波激励的频率,v0为活塞速度幅值.
减振器的试验条件可等效为幅值为0.52 m/s,频率为3.31 Hz的正弦速度激励. 减振器各节流孔直径、阀片尺寸及弹簧刚度等参数依据所设计的实际结构设定,仿真得到液压式馈能减振器的阻尼力曲线如图5所示. 根据QC/T 491—1999的要求,活塞直径为30 mm的液压减振器在以0.52 m/s的速度运动时,复原行程阻尼力峰值需要在820~3 288 N之间且压缩行程阻尼力峰值不得大于200 N. AMESim中仿真得到的复原行程最大可提供3 222 N的阻尼力,压缩行程最大可提供1 040 N的阻尼力,满足国标对工作缸直径30 mm的减振器阻尼力允许差值的要求.

图5 活塞速度幅值为0.52 m/s时的阻尼力曲线
Fig.5 Damping force curve when the velocity amplitude of piston is 0.52 m/s
进一步仿真得到此液压式馈能减振器的示功特性曲线和速度特性曲线,如图6所示.

(a) 示功特性曲线

(b) 速度特性曲线
图6 活塞速度幅值为0.52 m/s时的外特性曲线
Fig.6 External characteristic curve when the velocity amplitude of piston is 0.52 m/s
分析仿真结果能够发现复原/压缩行程阻尼力峰值都在减振器行程中点位置处取得,示功图饱满无畸变,减振器可及时提供充足的阻尼力以消减振动. 减振器正常工作时,压缩行程阻尼力较小,以便充分利用弹性元件的弹性来缓和冲击;复原行程阻尼力较大,以求迅速减振. 观察速度特性曲线可以看出,相同速度下压缩行程的阻尼力数值要略小于复原行程,满足设计要求.
为了更全面地剖析本文所设计的液压式馈能减振器的外特性,使用AMESim软件仿真得到不同活塞速度幅值下的示功特性曲线如图7所示,根据正弦激励的位移与速度关系得到的活塞速度幅值及对应的信号频率如表2所示.

图7 不同的活塞速度幅值下的示功特性曲线
Fig.7 Indicator characteristic curve under different velocity amplitude of piston
表2 减振器仿真工况
Tab.2 Damper simulation conditions

观察不同活塞速度幅值下的减振器示功特性曲线能够发现,本文设计的液压式馈能减振器示功图形状饱满,呈现出良好的外特性,仅在高速运动时略有畸变,但并不影响其正常的减振性能,外特性方面满足减振器的设计要求.
4 液压式馈能减振器的馈能特性
4.1 馈能特性定性分析
对馈能特性进行定性分析,以验证本文设计的液压式馈能减振器是否可以完成振动能量的回收. 仿真获得不同活塞速度幅值下高压蓄能器内气体压力随时间变化曲线,如图8所示. 分析图8能够发现,随着减振器工作时间的增加,蓄能器内气体压力从初始的0.6 MPa逐步升高,而且活塞速度幅值越高,高压蓄能器内气体压力上升的幅度越大. 由此可以证明,本文所设计的馈能式液压减振器能够实现能量回收.

图8 不同的活塞速度幅值下的蓄能器内气体的压力曲线
Fig.8 Pressure curve of the gas in the accumulator under different velocity amplitude of piston
同时提取高压蓄能器端口的油液流量曲线,如图9所示. 分析曲线可以发现,在压缩行程时流向高压蓄能器的油液流量呈现正弦型变化规律,在伸张行程时油液流量则为0,呈现显著的周期特性. 这表明本文设计的液压式馈能减振器在复原行程不回收能量而只在压缩行程回收能量,符合最初的设计思想. 而且高压蓄能器端口油液流量随减振器速度的增高而增大,证明减振器速度越大,高压蓄能器所回收到的液压能越多.

图9 高压蓄能器端口油液的流量曲线
Fig.9 Oil flow curve of high pressure accumulator port
4.2 馈能特性定量分析
一般用馈能功率和馈能效率两个评价指标来定量分析液压式馈能减振器的馈能特性[11]. 馈能功率即输入到蓄能器的功率,可以通过在仿真模型中的高压蓄能器端口连接能量传感器获得. 馈能功率的计算方法为
Pk=p·q.
(1)
式中:Pk为蓄能器回收功率,p为蓄能器入口压力,q为油液流量.
蓄能器回收到的能量来源于激励源的激励能,激励源的输出功率可以表示为
P=F·v.
(2)
式中:P为输出功率,F为减振器的实时阻尼力,v为减振器的实时速度.
根据式 (1)、(2),忽略在能量传递过程中系统的泄漏等其他环节的能量损失,馈能效率可表示为

(3)
经仿真得到减振器的馈能功率曲线和馈能效率曲线如图10所示. 观察馈能功率曲线能够发现,每一个周期的功率幅值由10 W逐渐接近12 W,呈现出缓慢增大的趋势. 这是由于蓄能器回收能量的叠加效果会使蓄能器内气体压力不停增大,因而计算出的馈能功率会有缓慢增大的趋势.

(a) 馈能功率曲线

(b) 馈能效率曲线
图10 活塞速度幅值为0.52 m/s时的馈能特性曲线
Fig.10 Feeding characteristic curve when the velocity amplitude of piston is 0.52 m/s
观察馈能功率曲线和馈能效率曲线可以发现,与高压蓄能器端口油液流量曲线类似,液压式馈能减振器的馈能功率和馈能效率也在正弦信号激励下具有周期特性. 并且曲线呈现出明显的峰值特性,即减振器的速度越大,馈能功率和馈能效率也越大. 因而可以推断出馈能减振器的高频响应特性不仅影响系统的阻尼特性,还影响系统的馈能特性.
5 液压式馈能减振器应用性能仿真
为了研究减振器特性与悬架匹配的问题,国内吉林大学的战敏、刘松山[12-13],上海交通大学的ZHANG G等[14],英国剑桥大学的PIRES L等[15]、美国明尼苏达大学的SINGAL K等[16]学者都采用先建立减振器阻尼力模型,再建立整车悬架模型的仿真方法进行的.
由理论分析和具体实践可知,独立悬架系统的振动基本上是独立的,因此可以用二自由度四分之一被动悬架汽车行驶动力学车辆模型来研究车辆在路面激励下的振动情况.
5.1 车辆及悬架模型参数
以哈弗H6为车辆原型,将1/4悬架模型中当量车身质量ms设置为350 kg,非悬挂质量mw设置为50 kg,悬架刚度Ks设置为15 000 N/m,车轮刚度Kw设置为180 000 N/m,悬架阻尼系数Cs设置为1 200 Ns/m.
5.2 悬架的联合仿真模型
路面输入模型分为频域模型和时域模型. 路面时域输入模型使用滤波白噪声来反映低频时路面谱的实际情况,可表示为

式中:q(t)表示路面位移,f0=0.1 Hz表示下截止频率,u表示汽车行驶速度,Sq(n0)表示路面不平度系数,w(t)表示均值为0、强度为1的均匀分布白噪声. B级路面的Sq(n0)上限为128×10-6,下限为32×10-6,几何平均值为64×10-6. 当车辆以100 km/h的速度在B级路面行驶时,得到时域路面不平度曲线如图11所示.

图11 时域路面不平度曲线
Fig.11 Time domain road roughness curve
运用AMESim和Simulink进行联合仿真,设置Visual C++作为编译器,仿真接口采用SimuCosim. AMESim中液压式馈能悬架仿真模型如图3所示,Simulink中搭建的联合仿真模型如图12所示.

图12 Simulink中液压式馈能悬架联合仿真模型
Fig.12 Joint simulation model of hydraulic energy regenerative suspension in Simulink
5.3 悬架平顺性仿真分析
设置路面激励为B级路面,车速为100 km/h,悬架平顺性仿真结果如图13所示.

(a) 车身垂直加速度曲线

(b) 悬架动挠度曲线

(c) 车轮动位移曲线
图13 传统悬架和馈能悬架关于平顺性指标的对比曲线
Fig.13 Contrast curve of conventional suspension and energy regenerative suspension on ride index
由图13能够看出,在数值和曲线趋势上,液压式馈能悬架和传统悬架的车身垂直加速度、悬架动挠度及车轮动位移都十分相近. 经数据处理发现,液压式馈能悬架的车轮动位移指标优于传统悬架,车身垂直加速度和悬架动挠度指标与传统悬架的差值也都不超过5%. 因而此液压式馈能悬架的平顺性基本满足要求.
5.4 悬架馈能特性仿真分析
仿真得到B级路面不同车速条件下的馈能功率和馈能效率如表3所示,车速40 km/h不同路面等级条件下的馈能功率和馈能效率如表4所示.
表3 不同车速下的馈能功率和馈能效率
Tab.3 Feeding power and feeding efficiency at different speeds

综合表3和表4能够看出馈能效率和馈能功率的变化趋势一致. 相同路面等级条件下,车速越高馈能悬架回收的能量越多;相同车速下,路面不平度越高馈能悬架回收的能量越多;而且横向比照能够发现路面不平度对馈能功率的影响更大.
表4 不同路面下的馈能功率和馈能效率
Tab.4 Feeding power and feeding efficiency at different pavement grades

6 结 论
1)设计了一种液压式馈能减振器的结构并分析其工作原理,可实现在车辆行驶的振动过程中,将油液压入高压蓄能器,将悬架的振动能量回收.
2)基于AMESim对液压式馈能减振器进行建模仿真分析,压缩/伸张行程阻尼力峰值均符合国标对相应条件下允差的要求. 速度特性曲线较理想,示功特性曲线形状饱满,呈现出稳定良好的外特性,符合工程实际需求.
3)正弦信号激励下的馈能特性分析表明,液压式馈能减振器的速度越大,馈能功率和馈能效率也越大,馈能特性符合预期的设计目标.
4)对液压式馈能悬架进行工况下联合仿真分析,结果表明,悬架平顺性满足行驶要求,实际行驶中具有一定的能量回收潜力.