作者简介:赵健(1978-),男,副教授,博士.研究方向:汽车地面系统分析与控制.E-mail:zhaojian@jlu.edu.cn
为适应路面附着条件的变化对车辆驱动力控制系统的要求,设计了基于纵滑侧偏组合刷子轮胎模型的路面识别算法,采用卡尔曼滤波方式估算轮胎与路面间的附着力,应用最小二乘法估算路面附着系数,进而确定当前路面对应的最佳车轮滑转率;在此基础上,设计了基于路面自适应的驱动力模糊滑模控制算法,并在Matlab/Simulink与CarSim联合仿真环境下进行了仿真验证。结果表明,本文算法能快速准确地识别路面变化,并根据路面附着情况有效调节车轮滑转率,改善车辆驱动性能。
In order to meet the requirements of variable road conditions for vehicle traction control system, a road identification algorithm was proposed based on combined longitudinal and lateral brushed tire model. First, the tire-road forces were estimated by the Kalman filtering method; the road friction coefficient on each wheel was estimated using least square method, and its optimal slip ratio was obtained subsequently. Then, a fuzzy sliding mode based on traction controller was designed using the identified road information. Finally, the algorithm was verified using Matlab/Simulink and CarSim co-simulation. Results show that the proposed algorithm can identify the variation of road conditions and regulate the slip ratios of the driving wheels quickly and accurately, thus, the vehicle driving pe4rformance can be effectively improved.
汽车驱动力控制系统在汽车起步、加速、爬坡过程中, 通过控制车轮滑转状态, 使车辆驱动轮在恶劣路面或复杂行驶条件下也能产生最佳的纵向驱动力, 并保证一定的转向能力[1]。目前国内外学者对于驱动力控制系统研究已进行了大量的工作, 例如, 赵健等[2]针对四轮驱动汽车, 应用模糊理论设计了牵引力控制系统的油门控制器和制动控制器。张晓兵等[3]提出面向牵引力控制的自动变速器档位干预策略。
这些算法大都采用固定车轮最佳滑转率进行控制, 忽略了路面变化对控制系统的影响。而实际上, 车辆行驶路况复杂, 不但路面附着系数多变, 不同附着系数路面对应的最佳车轮滑转率也有很大差别。因此, 对路面条件的精确识别是保证驱动力控制精度和品质的关键。
对于路面识别算法, 国内外进行了大量的研究。李亮等[4]提出了一种综合路面附着识别方法, 并引入了置信度的概念, 运用模糊逻辑算法实现了附着状态的判断。Takuro等[5]提出了在转弯工况下利用回正力矩及卡尔曼滤波来估算前轮所在路面的附着系数的方法。Edoardo等[6]介绍了倍耐力电子轮胎, 电子轮胎拥有内嵌的传感器和数字计算能力, 能实时地提供地面与轮胎的接触力、摩擦因数和附着系数的最大值。么鸣涛等[7]提出了基于车辆振动加速度响应的路面识别方法, 通过路面不平度的特征信息将路面进行分类。Choi等[8]提出了用线性递归最小二乘法来实时获取轮胎路面附着系数的方法。
虽然目前的路面附着系数估算方法很多, 但大部分都是面向车辆稳定性控制的, 而驱动力控制大多数采用固定的最佳滑转率控制, 难以在各种路面条件下均获得最佳的驱动性能。本文采用卡尔曼滤波和最小二乘法相结合的方式对路面附着系数进行估算, 自适应调节最佳车轮滑转率, 以此为基础, 设计模糊滑模控制器, 对汽车驱动力进行控制, 使之在不同的路面附着系数条件下, 都能获得最佳的驱动性能。
忽略空气阻力, 并将滚动阻力线性化, 汽车的受力分析如图1所示, 车轮的受力分析如图2所示。
车辆的力和力矩平衡方程为:
式中:F=
Hx=
Hy=
式中:Fxi为车轮纵向力(i=fl, fr, rl, rr分别对应左前、右前、左后、右后轮, 下同); Fyi为车轮侧向力; Fyf为前轴侧向力; Fyr为后轴侧向力; δ 为前轮转角; m为整车质量; ax为车辆纵向加速度; ay为车辆侧向加速度; Iz为车辆横摆转动惯量; γ 为车辆横摆角速度; lf为车辆质心到前轴距离; lr为车辆质心到后轴距离; 2t为车辆轮距。
车轮力矩方程为:
式中:Iw为车轮转动惯量; ω i为车轮轮速; Tdi为车轮驱动力矩; Tbi为车轮制动力矩; Fzi为车轮法向力; fri为滚动阻力系数; r为车轮有效滚动半径; G为车轮法向载荷; Fp为车轴作用于车轮的水平力; Tf为滚动阻力矩。
根据式(1)~(4), 建立状态空间方程:
式中:11维系统状态向量x(t)=[Fxfl, Fxfr, Fxrl, Fxrr, Fyf, Fyr, ω fl, ω fr, ω rl, ω rr, γ ]T; 7维测量矩阵z(t)=[ax, ay, ω fl, ω fr, ω rl, ω rr, γ ]T; 系统过程噪声w(t)是11维零均值白噪声向量; 测量噪声v(t)为7维零均值白噪声向量, 且w(t)和v(t)互不相关。
采用零阶保持的方法, 将车辆连续系统状态空间方程进行离散化, 得到相应的离散系统的状态空间方程如下:
进行离散卡尔曼滤波的系统估计过程如下:
状态一步预测:
状态估计:
滤波增益矩阵:
一步预测误差方差阵:
估计误差方差阵:
式中:矩阵W和V分别为过程噪声和测量噪声的协方差矩阵。只要给定初值
由于刷子模型的求解参数少于魔术公式等模型, 且能比较准确地描述如摩擦椭圆影响和饱和轮胎力等轮胎的非线性特征[9, 10], 本文采用纵滑侧偏组合刷子轮胎模型进行路面附着系数识别。
纵滑侧偏组合刷子轮胎模型的计算公式如下:
式中:Cx为轮胎纵滑刚度; Cα 为轮胎侧偏刚度; μ 为路面附着系数; α 为轮胎侧偏角; κ 为实际纵向滑转率; vw为轮心处的纵向速度。
刷子轮胎模型需要辨识的参数包括Cx、Cα 和μ , 根据文献[8]可知, Cx、Cα 与很多因素有关, 如轮胎尺寸、胎面宽度、胎面刚度、充气压力、负载等, 在短时间内变化很小, 可以近似看作趋近于定值, 而路面附着系数μ , 可以采用最小二乘法进行估算。
最小二乘估计使所有偏差的平方和达到最小, 可以有效抑制测量误差[11]。
首先将刷子轮胎模型写成如下非线性形式:
式中:y是测量值, 这里采用由卡尔曼滤波得到的轮胎力的数值, 即y=[Fx, Fy]T; f(k, θ )为刷子轮胎模型的表达式, 其中, θ 是刷子轮胎模型参数, θ =[μ ]; vl是相应的测量噪声。
为满足最小二乘法的需要, 将式(17)线性化, y(k)可以近似地写成:
式中:F
定义变量h(k):
将y(k)的表达式带入上式, 得到:
将非线性的刷子轮胎模型进行线性化处理后, 即可以采用最小二乘法估算方法来进行路面附着系数估算。一般说来, 基于刷子模型的路面附着系数估算算法需要轮胎存在一定的滑转率才能获得较精确的测量结果, 而在低滑转工况下, 估算结果会有一些抖动, 为了避免因此而产生错误的估算结果, 本文算法只在滑转率大于5%时进行估算值的更新, 滑转率小于5%时, 则保持上一个时刻的估算值不变。
根据文献[10], 在驱动工况下, 轮胎的理论纵向滑转率s与实际纵向滑移率κ 存在如下关系:
轮胎与路面间的附着系数随车轮滑转率变化, 由文献[12, 13]可知, 一般情况下, 峰值路面附着系数对应的最佳车轮滑转率在0.2左右, 而附着系数较高的路面对应的最佳车轮滑转率也较大。因此可以根据估算得到的路面附着系数通过查表的方式获得车轮的最佳滑转率s0, 最佳滑转率与路面附着系数的关系如下:μ ≥ 0.4时, s0=0.25; μ ≤ 0.1时, s0=0.15; 0.1< μ < 0.4时, s0=0.2。
选取实际滑转率与最佳车轮滑转率的差值e及其导数作为状态变量:
其切换函数δ slide为:
式中:c为滑模面系数。
考虑到控制的快速性和平顺性的要求, 采用连续函数和模糊调节结合的方式对轮胎的滑转率进行滑模变结构控制, 取:
式中:Ω 为滑模边界值, Ω > 0; ε 为趋近速率, ε > 0。
满足滑模变结构控制的全局达到条件, 则滑动模态存在, 控制方法可行。
为了减少抖振, ε 采用模糊调节的方式取值。模糊调节的输入为e和
求式(21)的一阶和二阶导数, 将
对式(4)两边求导, 得:
假设车轮驱动力矩已知, 将式(26)代入上式, 即可得到所需的驱动轮制动力矩变化率
在Matlab/Simulink和CarSim环境下建立联合仿真平台, 选择普通前驱D级轿车进行仿真, 对路面附着系数的估算结果和驱动力控制的效果进行分析验证。
仿真在附着系数为阶跃变化的路面上进行, 路面附着系数变化情况如下:行驶路程l=0~10 m时, μ =0.5; l=10~25时, μ =0.1; l=25~50时, μ =0.3。路面总长度为50 m, 初始车速为1 km/h, 油门开度在0.1 s内由0上升到100%, 仿真结果如图6~图11所示。
基于卡尔曼滤波算法得到的轮胎力估算值是路面附着系数估算算法的输入量, 因此路面附着系数估算的误差与轮胎力估算的误差有同样的变化趋势, 轮胎纵向力的估算精度越高, 路面附着系数的估算精度也越高。图6为实际的轮胎纵向力与估算的轮胎纵向力对比图, 由图可见, 估算的轮胎纵向力在真实值附近波动, 除路面跃变时有较大的瞬时估算误差外, 其余时刻的估算误差均在10%以下, 算法具有较高的估算精度。图7为实际路面附着系数与估算路面附着系数对比图, 由图可见, 采用本文算法, 在估算的最初阶段, 路面附着系数有一定偏差, 但估算值约在0.3 s内即趋于稳定, 而在路面发生跃变时, 估算值与实际值的最大瞬态偏差只有15%。在路面附着系数为0.5时, 估算结果与实际值之间约有6%的稳态误差, 这是由用于估算的刷子模型和CarSim中的轮胎模型的偏差造成的, 而在附着系数较低时, 估算值与真实值之间几乎没有稳态误差。可见算法可以快速适应路面条件的变化, 实现对轮胎力和路面附着系数的渐近无偏估计。
图8~图10分别为无驱动力控制、基于固定目标滑转率的普通滑模驱动力控制和基于本文设计的路面自适应模糊滑模驱动力控制条件下, 单侧驱动轮滑转率的变化曲线, 由于车辆在均一附着系数路面上行驶, 另一侧驱动轮滑转率曲线与图示曲线基本一致。由图8~图10可见, 无驱动力控制情况下, 驱动车轮迅速发生打滑, 严重影响车辆的起步性能, 而基于固定目标滑转率和固定滑动趋近率的滑模驱动力控制虽然可以较为有效地控制驱动轮的过度滑转, 但并不能将滑转率控制在最佳值, 且轮速存在较大抖动; 基于路面自适应识别的模糊滑模驱动力控制则可以快速有效地响应路面条件的变化, 平滑地将车轮滑转率控制在最佳值附近, 从而得到更佳的控制效果。本文设计的驱动力控制算法能将滑转率稳定控制在目标值附近(一般> 10%), 高于2.2小节中设置的更新估算值的滑转率门限(5%), 因此本文估算算法可以满足本文的驱动力控制算法的需求。
图11为本文设计的控制算法的车轮制动力矩曲线, 由图可见, 车轮制动力矩最大值达到了2350 N· m, 在轮速趋于稳定后, 制动力矩略有波动, 对照图10, 说明本文算法控制效果良好。
建立了离散车辆观测模型, 采用卡尔曼滤波方法进行了轮胎力的估算, 据此基于纵滑侧偏组合刷子轮胎模型, 采用最小二乘法估算得到了路面附着系数, 仿真结果表明, 这一算法可以快速准确地识别路面附着系数的跃变, 实现对路面附着系数的渐近无偏估计。车辆在低附着系数路面上全力加速时, 车轮会发生明显打滑, 而本文所设计的驱动力控制算法可以快速适应路面附着条件的变化, 将驱动轮滑转率始终控制在最佳值附近, 有效提高了车辆驱动性能。同时, 采用模糊规则调节滑模控制趋近速率, 可以减少车轮波动, 提升控制品质。
The authors have declared that no competing interests exist.
[1] |
|
[2] |
|
[3] |
|
[4] |
|
[5] |
|
[6] |
|
[7] |
|
[8] |
|
[9] |
|
[10] |
|
[11] |
|
[12] |
|
[13] |
|