油气悬架是将弹性元件与减振器合二为一,以蓄能器内部的惰性气体作为弹性介质,以液压油作为传力介质的减振装置[1]。它具有良好的非线性减振特性与较强的环境适应性,因此被广泛应用于运行工况复杂的重载车辆[2]。油气悬架性能的考核与校验主要依靠模拟试验台,其加载方式可显著影响预期载荷或路谱的复现能力。
西华大学司云飞[3]设计了一种采用机械加载且幅频在一定范围可控的汽车悬架试验台,并对试验台的应用进行了初步探讨。东北林业大学赵强[4]和秦玉彬[5]将前馈干扰补偿器和系统最优反馈控制器应用于伺服系统,设计了复合控制式悬架试验台。奇瑞汽车刘刚等[6]采用模糊自整定PID控制方法,完成了高精度电液伺服悬架试验台的设计。杜书等[7]应用RPC模拟迭代等技术,完成道路模拟迭代及耐久试验,验证了悬架道路模拟试验台的相关要求。张辉[8]引入三状态控制策略,应用极点配置的方法对控制器参数进行设计,在保持油气悬架模拟试验系统稳定性的基础上拓展了系统的频宽,改善了动态特性。
分析上述文献不难发现,现有试验台的加载方式主要为机械加载法及电液伺服加载法。其中机械加载法输出信号单一且误差较大,难以满足油气悬架的性能校核要求。电液伺服加载法虽然精度高,但是成本高昂且环境要求苛刻,不利于推广应用。因此有必要设计一种经济可靠、功能丰富且适用于重载车辆油气悬架性能校核的模拟试验台。
重载车辆油气悬架模拟试验台可完成油气悬架的位移、速度、刚度、阻尼、耐久等性能分析试验,以及1/4车辆装车道路模拟试验。其主要性能指标为:最大载荷5 t,最大行程±120 mm,激振频率范围0~30 Hz,活塞杆运动的最大速度0.5 m/s[9-12]。模拟试验台由机械台架、液压加载系统及电液比例控制系统构成,各分系统关系如图1所示。
图1 各分系统功能关系图
Fig.1 Functional relationship diagram of each subsystem
模拟试验台机械台架由底座、光轴立柱、上下横梁、上下滑块、滑块锁紧器、砝码板以及传感器等部件构成,如图2所示。可根据性能分析试验或1/4车辆装车道路模拟试验的不同,将油气悬架、车轮等重载车辆构件安装在不同位置,实现不同模拟试验的功能。
图2 油气悬架模拟试验台三维结构图
Fig.2 Three-dimensional structure diagram of hydro-pneumatic suspension simulation test bench
液压加载系统由液压源、电液比例阀、液压作动缸以及蓄能器等元件组成,主要作用是向试验台台面输出路面激振波形,系统原理图如图3所示。其中,电液比例阀为核心元件,不仅可以控制液压油流向,还可以通过外部电信号成比例地控制液压油流量,实现方向及流量控制元件的合二为一。
1.油箱;2.液压泵及电机;3.溢流阀;4.压力表;5.过滤器;6.单向阀;7.蓄能器;8-1.手控换向阀;8-2.电液比例换向阀;9.液控单向阀;10.液压作动缸
图3 电液比例液压加载系统原理图
Fig.3 Schematic diagram of electro-hydraulic proportional hydraulic loading system
电液比例控制系统由比例放大器、电液比例换向阀、液压作动缸及位移传感器等环节构成,用于控制液压加载系统的波形输出,系统控制框图如图4所示。
图4 电液比例控制系统框图
Fig.4 Block diagram of electro-hydraulic proportional control system
比例放大器用于将输入信号的电压形式转换为电流形式,此过程为比例环节,传递函数为[13]
(1)
式中:I和U分别表示比例放大器的输出电流与输入电压,单位分别为A与V;Ka表示比例放大器增益(A/V)。
电液比例换向阀采用电流信号控制比例电磁铁的位置,改变阀口的大小,从而控制流入与流出作动缸油液的流量,此为电信号到液压信号转换的过程[14]。此过程可视为二阶环节,传递函数为[15]
(2)
式中:QL表示电液比例换向阀的输出流量(m3/s);Ksv表示流量—电流增益[m3/(s·A)];ωsv表示此环节固有频率(rad/s);ζsv表示此环节的阻尼比,取 ζsv =0.6。
液压作动缸原理图如图5所示。图中液压源的供油压力为ps,回流压力为p0,电液比例阀的阀芯位移为xv,作动缸无杆腔压力、流量、有效面积分别为p1、q1、A1,有杆腔压力、流量、有效面积分为p2、q2、A2。作动缸活塞杆位移为y,以活塞杆伸出方向为正方向进行研究。回流压力视为零,即p0=0。
图5 阀控液压作动缸原理图
Fig.5 Schematic diagram of valve-controlled hydraulic actuating cylinder
2.3.1 比例阀线性流量方程
活塞杆伸出时,作动缸活塞杆的输出力为
FL=p1A1-p2A2
(3)
令FL=pLA1,其中pL定义为负载压力,则上式可化为
pLA1=p1A1-p2A2
(4)
即:
(5)
定义面积比n=A2/A1,则上式可化为
pL=p1-np2
(6)
活塞杆伸出时,无杆腔及有杆腔的流量分别为[16]
(7)
(8)
由液压作动缸的工作原理,可知活塞杆运动时两腔的流速相等,因此两腔流量与有效面积成正比为
(9)
联立式(6)与式(9),用ps与pL替换无杆腔与有杆腔的压力p1、p2为
(10)
(11)
由于油液的可压缩性和液压缸外泄漏的影响,液压缸的流出与流入流量不相等,将负载流量定义为两腔流量的均值为
(12)
代入式(7)与式(8),可得:
(13)
代入式(10)与式(11),用ps与pL替换p1、p2,即可得:
(14)
比例阀线性流量方程为
qL=Kqxv-KcpL
(15)
式中,定义流量增益(m2/s)为
(16)
定义流量—压力系数[m5/(Pa·s)]为
(17)
2.3.2 作动缸流量连续方程
两腔流量的连续性方程为[17]
(18)
式中:Cip与Cep分别表示作动缸内泄漏与外泄漏系数[m5/(N·s)];V1、V2分别表示作动缸无杆腔与有杆腔容积(m3);βe表示液压油弹性模量(Pa);y表示活塞的位移(m)。
两腔的工作容积为
(19)
V10、V20分分别为两腔的初始容积(m3)。
当活塞处于液压缸的中间位置时,油液压缩性影响最大,系统稳定性最差,故活塞杆任何其他位置相对于中位都是相对安全的,将中位设定为液压缸的初始位置,此时
(20)
式中:V为液压缸总容积(m3)。
根据液压作动缸特性,易知:
V2=(V-V1)n
(21)
联立方程式(10)—(13)、式(16)—(19),得到液压缸负载流量连续性方程为
(22)
式中,定义如下:
(23)
(24)
等效容积:
(25)
令Am= (A1+A2)/2,则式(22)可变换为
(26)
2.3.3 作动缸负载平衡力方程
当活塞杆伸出时,作动缸输出力与负载的平衡方程为
(27)
式中:mt表示活塞及负载折算到活塞上的总质量(kg);Bp表示活塞及负载的阻尼系数[N/(m·s)-1];Kp表示负载弹簧刚度(m/s);Fex表示作用在活塞上的任意外负载力(N)。
2.3.4 作动缸活塞杆输出位移的表达式
对式(15)、式(26)和式(27)进行拉普拉斯变换,得:
(28)
根据式(28),令Kce=Kc+Cie,称为系统流量系数,消去中间变量PL(s)与QL(s)可得作动缸活塞杆输出位移Y(s)的表达式为
上式上下同除以A1Am,可化为
(30)
式中:Ky=(4βeA1Am)/Vt,称为液压弹簧刚度,是由两腔的液体受到压缩而产生的刚度; Kcem=(A1Am)/Kce。
模拟试验台加载系统液压作动缸的负载以惯性负载为主,弹性负载和阻性负载可忽略不计,因此式(30)可简化为
(31)
式中:ωp为液压固有频率,单位为rad/s,即:
(32)
ζp为液压阻尼比,即:
(33)
2.3.5 液压作动缸环节的传递函数
在式(31)中,分子的第一项为作动缸活塞受电液比例阀移动的影响下输出的位移,以阀芯位移为输入,则传递函数为
(34)
以阀芯输出流量为输入,则传递函数为
(35)
式(31)中,分子的第二项为在外负载力的作用下,活塞输出位移的变化量,其传递函数为
(36)
位移传感器是一种金属感应的线性器件,可看作比例环节[18]。此环节的传递函数为
(37)
式中:Hd为位移传感器的增益系数(V/m);Ud表示位移传感器输出电压(V);Y为液压作动缸活塞杆位移(m)。
根据式(1)、(2),(35)—(37)的电液比例控制系统各环节的传递函数,通过计算可知系统该系统的开环传递函数为:
(38)
系统主要参数取值如表1所示。
表1 电液比例液压加载系统主要参数
Table 1 Main parameters of the electro-hydraulic proportional hydraulic loading system
参数取值参数取值Ka/(V·m-1)0.2βe/Pa7×108Ksv/(m3·(s·A)-1)0.09Kq0.05ζsv0.6Hd1ωsv/(rad·s-1)314mt/kg5 100A1/m23.12×10-3ωp/(rad·s-1)79.53A2/m22.41×10-3ζp0.461
电液比例液压加载系统的传递函数框图如图6所示。利用Matlab绘制电液比例液压加载系统传递函数的开环Bode图,如图7所示。电液比例液压加载系统的幅值裕度为11 dB,幅值穿越频率为69.9 rad/s;相位裕度为74.9°;相位交界频率为16.7 rad/s。根据电液比例控制系统性能判定依据,该幅值裕度大于6 dB;相位裕度大于45°,故动态过程收敛且性能良好。
图6 控制系统传递函数框图
Fig.6 Control system transfer function block diagram
图7 电液比例控制系统开环Bode图
Fig.7 Open loop Bode diagram of electro-hydraulic proportional control system
使用Simulink建立电液比例控制系统模型,以幅值0.1的阶跃信号作为输入激励,响应时间1 s,仿真模型与响应特性曲线分别如图8、图9所示。根据图9可知,电液比例液压加载系统对于阶跃信号的响应平顺性较好,表现为响应曲线较为平滑且无超调;且响应速度较快,完全消除稳态误差的时间为0.35 s,满足油气悬架模拟试验的要求。
图8 电液比例液压加载系统阶跃响应模型
Fig.8 Step response model of electro-hydraulic proportional loading system
图9 阶跃响应特性曲线
Fig.9 Step response characteristic curve
1) 本文中根据重载车辆油气悬架的性能指标要求,运用电液比例控制技术,设计了一套多功能复合式模拟试验台。
2) 通过建立电液比例液压加载系统的数学模型,推导得到了系统的传递函数。
3) 利用Matlab绘制Bode图分析该系统的稳定性;同时,利用Simulink搭建系统仿真模型分析其动态特性。结果表明:电液比例液压加载系统稳定性良好、响应速度快且无超调,可以满足油气悬架模拟试验的加载需求。
4) 本文中油气悬架模拟试验台设计与电液比例液压加载系统特性分析,为下一步重载车辆油气悬架模拟试验台的研制和性能试验的开展提供了参考依据。
[1] 陈开考,庞志康.汽车构造与拆装(下)[M].北京:机械工业出版社,2010:70-71.
CHEN Kaikao,PANG Zhikang.Automotive construction and disassembly (part 2)[M].Beijing:Mechanical Industry Press,2010:70-71.
[2] 汽车工程手册编辑委员会.汽车工程手册-设计篇[M].北京:人民交通出版社,2001.
Automotive engineering handbook editorial committee.automotive engineering handbook-design section[M].Beijing:People’s Communications Press,2001.
[3] 司云飞.悬架匹配试验台设计及应用研究[D].成都:西华大学,2016.
SI Yunfei.Design and application research of suspension matching test bench[D].Chengdu:Xihua University,2016.
[4] 赵强,秦玉彬,孙政等.电液伺服式悬架试验台的模拟仿真研究[J].森林工程,2020,36(3):60-68.
ZHAO Qiang,QIN Yubin,SUN Zheng,et al.Simulation research on the electrohydraulic servo suspension test bench[J].Forest Engineering,2020,36(3):60-68.
[5] 秦玉彬.基于电液伺服系统的路谱复现和悬架控制[D].哈尔滨:东北林业大学,2020.
QIN Yubin.Road spectrum reproduction and suspension control based on electro-hydraulic servo system[D].Harbin:Northeast Forestry University,2020.
[6] 刘刚,高家兵,聂高法,等.基于模糊PID的悬架试验台激振控制系统仿真[J].农业装备与车辆工程,2015,53(2):51-55,67.
LIU Gang,GAO Jiabing,NIE Gaofa,et al.J.Simulation of suspension test bench excitation control system based on fuzzy PID[J].Agricultural Equipment and Vehicle Engineering,2015,53(2):51-55,67.
[7] 杜书,巫洋,刘明,等.乘用车悬架系统道路模拟试验台的研究[J].汽车文摘,2020(6):49-52.
DU Shu,WU Yang,LIU Ming,et al.Research on road simulation test bench for passenger car suspension system[J].Automotive Abstracts,2020 (6):49-52.
[8] 张辉.车辆油气悬架电液伺服加载系统控制方法研究[D].哈尔滨:哈尔滨工业大学,2012.
ZHANG Hui.Research on the control method of electrohydraulic servo loading system for vehicle hydropneumatic suspension[D].Harbin:Harbin Institute of Technology,2012.
[9] GB/T 7031—2005,机械振动 道路路面谱测量数据报告[S].北京:2006.
GB/T 7031—2005,Mechanical vibration road surface spectrum measurement data report[S].Beijing:2006.
[10] JT/T 448—2021,汽车悬架装置检测台[S].北京:2021.
JT/T 448—2021,Automotive suspension device testing platform[S].Beijing:2021.
[11] 汽车减振器性能要求及台架试验方法:QC/T491—2018[S].北京:2018.
Performance requirements and bench test methods for automotive shock absorbers:QC/T491—2018[S].Beijing:2018.
[12] GB1589—2004,道路车辆外廓尺寸、轴荷及质量限值[S].北京:2004.
GB1589—2004,Boundary dimensions,axle load,and mass limits for road vehicles[S].Beijing:2004.
[13] 郑凯锋,陈思忠.比例阀控制非对称液压缸系统的非线性建模与仿真[J].液压与气动,2013(4):25-29.
ZHENG Kaifeng,CHEN Sizhong.Nonlinear modeling and simulation of asymmetric hydraulic cylinder systems controlled by proportional valves[J].Hydraulic and Pneumatic,2013,2013(4):25-29.
[14] WU Y,JIANG T,LU Q,et al.Characteristics of thermally actuated pneumatic proportional pressure valves and their application[J].Journal of The Institution of Engineers (India) Series C,2020,101(55):1-11.
[15] ZHAO W,LI L.A new k-dops collision detection algorithms improved by ga[C]//Lecture Notes of the Institute for Computer Sciences,Social Informatics and Telecommunications Engineering. Springer,Berlin,Heidelberg,2012.
[16] 刘保杰,强宝民,权辉.电液比例位置控制系统建模与仿真[J].液压气动与密封,2011,31(11):45-49.
LIU Baojie,QIANG Baomin,QUAN Hui.Modeling and simulation of electrohydraulic proportional position control system[J].Hydraulic Pneumatic and Sealing,2011,31 (11):45-49.
[17] 汪亮培.采煤机电液比例调高控制仿真研究[D].淮南:安徽理工大学,2019.
WANG Liangpei.Simulation study on hydraulic proportional adjustment control of coal mining machinery[D].Huainan:Anhui University of Science and Technology,2019.
[18] LI M,REN G S.Research on electro-hydraullic proportional self-tuning fuzzy-pid control system based on matlab[J].Applied Mechanics and Materials,2013(423/426):2882-2885.