摘要
轮胎与路面的接触关系的模拟影响着车辆悬架响应和控制的精度.在运用分形理论构建的三维路面模型基础上,建立了轮胎与路面的三维动态接触关系.针对1/4车辆悬架模型,分别搭建了被动悬架、模糊PID控制和天棚阻尼控制算法,在两种接触条件下对比了被动悬架和两种控制算法的车辆悬架垂向加速度、悬架动挠度及轮胎力的时域和频域响应情况.计算结果表明,采用模糊PID控制算法在车辆平顺性方面优于被动悬架和天棚阻尼控制算法;车辆在考虑轮胎和路面三维接触时,在面接触工况下,车身加速度、悬架动挠度及轮胎力响应的峰值都低于点接触工况.
轮胎与道路表面在接触印迹上的面积,虽然只有成人手掌大小,但车辆动力学性能的发挥主要受轮胎动态特性的影响,即通过轮胎/路面的接触面产生的三向力、回正力矩等实现载荷传递和操纵功
本文针对车辆悬架的平顺性优化问题,建立了如

图1 1/4车辆模型
Fig.1 1/4 vehicle model
(a)被动悬架 (b)半主动悬架
(a)Passive suspension (b) Semi-active suspension
(1) |
半主动悬架动力学方程为
(2) |
半主动悬架系统总的阻尼力为
(3) |
可变阻尼力为
(4) |
式中,为簧下质量,为簧上质量, 为轮胎刚度,为悬架弹簧刚度,为阻尼系数, 为可调阻尼力,为路面不平度, 为可调阻尼.本文选取二自由度1/4汽车悬架模型参数,如
作为车辆振动输入路面不平度,主要采用功率谱密度来描述其统计特性.可以将路面速度功率谱密度幅值在整个频率范围为一常数,即为一白噪声.计算公式如下:
(5) |
式中,
——下截止空间频率,;
——路面不平度系数();
——均值为零的白噪声;
——路面随机高程位移(m);
——车辆的行驶速度().
根据上式路面不平度时域的表达,当车辆行驶速度v=20m/s时,搭建二维随机路面模型,如

图2 二维随机路面(C级)
Fig.2 Random road surface (Grade C)
运用分形插值原理,从空间上将二维路面拓展到三维路面,构建了对照于国标中C级路面激励的三维模型路面谱,如

图3 三维路面谱
Fig.3 3D pavement spectrum

图4 三维接触模型
Fig.4 3D contact model
为了探究路面和轮胎接触情况,本文进行了轮胎接地面积测试,以重载子午线轮胎10.00R20为例,在标准胎压830 kPa和载荷F=30 000 N下,轮胎与地面的接触面宽度和长度分别为200 mm和230 mm,路面与轮胎接触点数为725个,纵向和横向分别25个点和29个点,具体过程参见文献[
悬架是汽车的重要组成部分,对汽车的平顺性与行驶安全性具有重要意义.本文采用天棚阻尼控制和模糊PID控制两种策略,并与被动悬架进行对比,分别在轮胎和路面的点接触和面接触工况条件下,进行了时域和频域分析,对比了两种工况下的车辆的垂向动态性能.
天棚阻尼控制策略是将阻尼器安装在车体与某固定的“天棚”之间,当阻尼达到一定值时可以起到减振的作用,在汽车悬架中采用广泛.其模型如

图5 天棚结构原理图
Fig.5 Schematic diagram of the shed structure
(6) |
比例积分微分控制,简称PID(Proportional I- ntegral Derivative),PID控制以其算法整体简便、可靠性高、应用范围广、控制效果较好,在工业控制中被广泛应用.PID控制的数学表达式为:
(7) |
式中,e(t)是系统误差;Kp是比例系数,表示对系统误差信号比例的加权;Ki是积分系数,表示对系统误差信号积分的加权;Kd是微分系数,表示对系统误差信号微分的加权.
而模糊PID控制是在PID控制规律基础上,以车身垂直加速度与理想加速度值(设为0)之间的偏差(e)和其变化率(ec)作为控制器的输入,半主动悬架系统的可调阻尼力作为控制器的输出,满足不同时刻e和ec来调整PID控制器的三个参数值,结合1/4车辆模型,模糊PID控制模型如

图6 模糊PID控制仿真结构图
Fig.6 Fuzzy PID control simulation structure
在点接触工况下,重载车辆分别采用被动悬架、半主动悬架天棚阻尼振动控制和半主动悬架模糊PID控制三种情况,对车体振动加速度、悬架动挠度及轮胎力进行分析,结果如

(a) 时域图
(a) Time domain diagram

(b) 频域图
(b) Frequency domain map
图7 点接触下车身加速度
Fig. 7 Body acceleration under point contact

(a) 时域图
(a) Time domain diagram

(b) 频域图
(b) Frequency domain map
图8 点接触下悬架动挠度
Fig. 8 Suspension deflection under point contact

(a) 时域图
(a) Time domain diagram

(b) 频域图
(b) Frequency domain map
图9 点接触下轮胎动载荷
Fig.9 Tire dynamic load under point contact
(1)从时域图中可知,天棚阻尼控制和模糊PID控制算法下车身加速度峰值小于被动悬架车身加速度峰值,其控制算法下均方根值比被动悬架低19.9%和37.16%.
(2)从频域图中可知,两种控制算法和被动悬架下的主频率相近,但被动悬架对应峰值明显高于其他两种控制算法.
(1)从时域图中可知,天棚阻尼控制和模糊PID控制算法下悬架动挠度峰值小于被动悬架车身加速度峰值,其控制算法下均方根值比被动悬架低23.08%和42.95%.
(2)从频域图中可知,两种控制算法和被动悬架下的主频率相近,但被动悬架对应峰值明显高于其他两种控制算法.
(1)从时域图中可知,天棚阻尼控制算法下轮胎力峰值高于被动悬架峰值,其控制算法下均方根值比被动悬架高13.1%,该控制算法对 轮胎力控制不能达到理想效果;模糊PID控制算法下轮胎力峰值低于被动悬架峰值,其控制算法下均方根值比被动悬架低17.23%.
(2)从频域图中可知,在低频下,两种控制算法和被动悬架下的主频率相近,被动悬架对应峰值明显高于其他两种控制算法;在高频下,模糊PID控制算法和被动悬架下的主频率几乎相同,天棚阻尼控制算法对应峰值明显大于被动悬架.
在面接触工况下,同样分别采用半主动悬架天棚阻尼控制和半主动悬架模糊PID控制算法,并与被动悬架对比,对车体振动加速度、悬架动挠度和轮胎力进行分析,结果如

(a) 时域图
(a) Time domain diagram

(b) 频域图
(b) Frequency domain map
图10 面接触下车身加速度
Fig.10 Body acceleration under surface contact

(a) 时域图
(a) Time domain diagram

(b) 频域图
(b) Frequency domain map
图11 面接触下悬架动挠度
Fig.11 Suspension deflection under surface contact

(a) 时域图
(a) Time domain diagram

(b) 频域图
(b) Frequency domain map
图12 面接触下轮胎动载荷
Fig. 12 Tire dynamic load under surface contact
将面接触下车辆的各动力学响应情况汇总为
(1)从时域图中可知,模糊PID控制算法下加速度的峰值低于被动悬架车身加速度峰值,其控制算法均方根值低于被动悬架37.66%;天棚阻尼控制加速度的峰值高于被动悬架车身加速度峰值,但其控制算法均方根值低于被动悬架5.54%,天棚阻尼控制算法对车身加速度的控制达不到预期效果.
(2)从频域图中可知,在低频下,两种控制算法和被动悬架下的主频率相近,但被动悬架对应峰值高于其他两种控制算法,在高频60Hz左右,天棚阻尼控制算法下曲线出现明显波动.
(1)从时域图中可知,天棚阻尼控制和模糊PID控制算法下悬架动挠度峰值小于被动悬架峰值,其控制算法下均方根值比被动悬架低14.96%和29.92%.
(2)从频域图中可知,两种控制算法和被动悬架下的主频率相近,但天棚阻尼控制算法对应峰值高于被动悬架对应峰值,模糊PID控制算法对应峰值低于被动悬架对应峰值.
(1)从时域图中可知,模糊PID控制算法下轮胎力的峰值低于被动悬架轮胎力峰值,其控制算法均方根值低于被动悬架36.71%;天棚阻尼控制下轮胎力高于被动悬架轮胎力,但其控制算法均方根值高于被动悬架123%,天棚阻尼控制算法不能有效对轮胎力进行控制.
(2)从频域图中可知,在低频下,两种控制算法和被动悬架下的主频率相近,但被动悬架对应峰值高于其他两种控制算法,在高频下,天棚阻尼控制算法下曲线对应峰值远大于被动悬架.
通过上一节分析,半主动悬架天棚阻尼控制在两种工况下对轮胎力和车身加速度控制效果都不好.为了更加清晰直观看出点接触和面接触两种工况对车辆平顺性影响,通过仿真,将两种工况的反应曲线分别在被动悬架和模糊PID控制算法动力学响应图像中显现.结果如

(a) 车身加速度
(a) Body acceleration

(b) 轮胎力
(b) Tire force

(c) 悬架动挠度
(c) Suspension deflection
图13 被动控制仿真结果
Fig. 13 Passive control simulation results

(a) 车身加速度
(a) Body acceleration

(b) 轮胎力
(b) Tire force

(c) 悬架动挠度
(c) Suspension deflection
图14 模糊PID控制仿真结果
Fig. 14 Fuzzy PID control simulation results
从
由
在被动控制算法下,把轮胎和路面看成多点接触工况,各参数的均方根值均比点接触工况都低.车体加速度、悬架位移和轮胎力的均方根值分别小18.38%、18.59%、39.73%.
由
本文具体研究了在模糊PID控制、天棚阻尼控制和被动悬架状态下,轮胎与路面单点接触和多点面接触两种不同接触模型下车辆悬架的响应区别.天棚阻尼控制算法对悬架动挠度和轮胎力不能达到期望效果;模糊PID控制算法优于天棚阻尼控制算法和被动悬架;车身振动、悬架动挠度和轮胎力方面,面接触工况都小于点接触工况响应.
参考文献
张建莉,陈榕峰,姚岢.公路路面不平整度对汽车振动的影响分析.公路交通科技,2019,36(11):129~133 [百度学术]
Zhang J L,Chen R F, Yao K.Analysis on Influence of road pavement roughness on vehicle vibration. Journal of Highway and Transportation Research and Development,2019,36(11):129~133 (in Chinese) [百度学术]
段虎明,石峰,谢飞,等.路面不平度研究综述.振动与冲击,2009,28(9):95~101,216 [百度学术]
Duan H M, Shi F,Xie F,et al.A survey of road roughness study. Journal of Vibration and Shock,2009,28 (9):95~101,216(in Chinese) [百度学术]
刘永臣,王国林,孙丽,等. 面向用户的汽车道路谱重构方法. 中国机械工程,2014,25(15): 2112~2116 [百度学术]
Liu Y C, Wang G L, Sun Li, et al. Reconstruction methods of vehicle road spectrum faced users. China Mechanical Engineering,2014,25(15): 2112~2116(in Chinese) [百度学术]
Dodds C J,Robson J D.The description of road surface roughness. Journal of Sound and Vibration,1973,31(2):175~183 [百度学术]
全国机械振动与冲击标准化技术委员会. GB/T7031-2005, 机械振动—道路路面谱测量数据报告. 北京: 中国标准出版社, 2006(National Mechanical Vibration and Shock Standardization Technical Committee. GB/T7031-2005, Mechanical vibration-road surface profiles-reporting of measured data.Beijing: China Standard Press, 2006(in Chinese)) [百度学术]
赵济海,王哲人,关朝雳. 路面不平度的测量、分析与应用.北京:北京理工大学出版社,2000(Zhao J H, Wang Z R, Guan C L. Measurement, analysis and application of road unevenness. Beijing:Beijing Institute of Technology Press,2000(in Chinese)) [百度学术]
Michelberger P,Palkovic L,Bokor J. Robust design of active suspension system. International Journal of Vehicle Design,1993,14(2 /3):145~165 [百度学术]
王亚,陈思忠,郑凯锋.时空相关路面不平度时域模型仿真研究.振动与冲击,2013,32(5):70~74 [百度学术]
Wang Y, Chen S Z, Zhen K F. Simulation research on time domain model of road roughness with time-space correlation. Journal of Shock and Vibration,2013,32 (5):70~74 (in Chinese) [百度学术]
王汉平,张哲,李倩.路面时域模拟的谐波叠加组分中相位角的相干性研究.北京理工大学学报,2019,39(10):1034~1038 [百度学术]
Wang H P,Zhang Z, Li Q.Coherence of phase angle of harmonic superposition components in time domain simulation of road roughness for right and left wheel tracks. Transactions of Beijing Institute of Technology, 2019,39(10):1034~1038 (in Chinese) [百度学术]
鲁植雄,徐浩,刘奕贯,等.基于分形插值的三维路面重构与分析.农业工程学报,2014,30(22):188~194 [百度学术]
Lu Z X,Xu H,Liu Y G,et al. Reconstruction and analysis of 3D pavement based on fractal interpolation. Journal of Agricultural Engineering, 2014, 30(22): 188~194 (in Chinese) [百度学术]
徐中源. 基于实测路面信息的三维虚拟路面重构研究[硕士学位论文].青岛:青岛大学,2018 [百度学术]
Xu Z Y. Research on 3D virtual pavement(reconstruction based on actual pavement information[Master Thesis]. Qingdao :Qingdao University, 2018 (in Chinese)) [百度学术]
路永婕,李振宇,怀文青,等.基于分形理论的三维路面谱重构及在多体动力学软件中的应用.图学学报,2019,40(2):328~334 [百度学术]
Lu Y J, Li Z Y, Huai W Q,et al.Reconstruction of three-dimensional pavement spectrum based on fractal theory and its application in Multi-Body dynamics software. Journal of Graphics,2019,40(2):328 ~334(in Chinese) [百度学术]
Lu Y J, Huai W Q, Zhang J N. Construction of three-dimensional road surface and application on interaction between vehicle and road.Shock and Vibration,2018,2535409 [百度学术]