VEHICLE TRACTION FORCE CONTROL BASED ON THE ROAD ADHESION COEFFICIENT ESTIMATION BY FFUKF
-
摘要: 以后驱牵引车为研究对象, 设计了基于路面附着系数估计的牵引力控制系统(TCS). 在路面附着系数估计方面, 针对传统卡尔曼滤波难以跟踪时变非线性系统的问题, 本文将模糊控制理论和衰减记忆滤波思想引入无迹卡尔曼滤波, 设计一种基于模糊遗忘因子的无迹卡尔曼滤波路面附着系数估计方法, 提高了算法的跟踪性能. 牵引力控制包括扭矩控制和制动控制. 在TCS扭矩控制方面, 分别利用路面附着系数和驱动轮滑转率在目标滑转率附近时的车辆加速度计算目标基础扭矩, 根据车辆行驶状态和抖振度参量, 基于可拓控制理论划分经典域、可拓域和非域, 通过可拓集的关联函数得到动态权重系数, 将上述两种方法计算得到的目标基础扭矩进行可拓融合设计出基础扭矩. 之后, 以实际滑转率和目标滑转率之间的误差作为输入, 采用模糊自整定PI控制器得到目标反馈扭矩. 在制动控制方面, 针对两种典型路面分别设计了PI控制压力和附着差压力. 实车试验结果表明, 基于模糊遗忘因子的无迹卡尔曼滤波算法能够更加快速地跟踪路面附着系数的变化, 同时基于路面附着系数估计的TCS控制策略能够有效抑制驱动轮过度滑转, 将驱动轮滑转率控制在最佳范围内, 显著提高了车辆的动力性.Abstract: Traction control system (TCS) based on road adhesion coefficient estimation is designed for the rear drive tractor. Firstly , to track the time-varying nonlinear system in the classical Kalman filter, a fuzzy control theory and an attenuated memory filter are introduced into untraced Kalman filter, and the fuzzy forgetting factor unscented Kalman filter (FFUKF) is proposed to estimate road adhesion coefficient, so as to improve the tracking performance of the algorithm. Traction control includes torque control and braking control. Then, for the TCS torque control, the road adhesion coefficient and the vehicle acceleration when the slip rate of the driving wheel is near the target slip rate are used to calculate the target base torque, respectively. According to the vehicle state and chattering parameters, the extension set is divided into the classical domain, extension domain and non-domain with the extension control theory, and the dynamic weight coefficient is obtained by the correlation function of extension set. The base torque is designed by extension fusion of the target base torques calculated by the above methods. Next, the error between the actual slip rate and the target slip rate is used as the input, and the target feedback torque is obtained by using the fuzzy self-tuning PI controller. In terms of braking control, the PI pressure control and adhesion difference pressure are designed on two typical road surfaces, respectively. Finally, the test results show that the proposed FFUKF algorithm can track the tire-road friction coefficient more quickly, and the proposed TCS control strategy can effectively restrain the excessive slip of the driving wheels to control the slip rate of drive wheel in the best range, so that improving vehicle dynamic performance significantly.
-
引 言
汽车牵引力控制系统 (traction control system, TCS)是汽车主动安全控制关键技术之一, TCS系统通过发动机输出转矩控制和主动制动控制调节驱动轮滑转率在最优范围内, 在实现驱动稳定性前提下提高汽车的加速性能[1]. 汽车—驾驶员—道路组成的系统具有瞬变的特点, 若想达到较好的控制效果, 对其关键状态参量的估算尤其是对路面状态的估计至关重要[2].
国内外诸多学者对路面附着系数的估计做了广泛的研究工作[3-5]. 王博等[6]指出目前路面附着系数估计主要分为Cause-based和Effect-based两类. Cause-based主要使用传感器(光学、声学、微波等)测量影响路面附着系数的主要因素, 通过经验模型估计附着系数[7-8], 但此方法受工作环境等影响较大, 并且成本较高, 难以商用. Effect-based主要通过建立观测器分析附着系数, 对车轮或车身运动响应进行估计, 目前多数研究属于Effect-based方法. 基于Effect-based的附着系数估计算法主要有基于
$ \mu - s $ 曲线斜率的估计方法[9-10]、卡尔曼滤波[11-13]等. 与扩展卡尔曼滤波(extended Kalman filter, EKF)相比, 无迹卡尔曼滤波(unscented Kalman filter, UKF)采用真实的非线性模型, 用随机变量逼近状态分布, 而不是将非线性函数线性化, 可以处理不可导的非线性函数, 保证了轮胎在非线性情况下的滤波精度[14], 适用于轮胎−路面附着系数这种强非线性系统. Chen等[15]基于改进的Dugoff轮胎模型, 利用无迹卡尔曼滤波和MSE加权融合获得道路附着系数, 但传统的UKF具有记忆性, 如果过滤器的初始化不够接近实际状态, 它可能会发散, 同时旧数据的累积会减小新获得的数据对滤波的更新作用, 因此导致UKF无法快速跟踪系统参数的变化[16-18].在TCS控制器方面, 对控制方法研究集中在 PID、逻辑门限、模糊逻辑、神经网络以及滑模控制等[19-20], 而以路面附着系数估计为基础进行TCS控制能够实现更好的控制效果. Ding等[21]提出了一种基于模糊逻辑理论的轮胎−路面附着系数估计方法, 在此基础上综合PID控制器对车轮滑转率进行控制, 实现了良好的控制效果. 在多数研究中TCS控制依赖路面附着系数估计, 一旦路面附着系数估计不准确, 控制效果会急剧下降.
高精度、快收敛的路面附着系数估计算法能够有效提升TCS控制效果, 本文将模糊控制和衰减记忆滤波思想引入UKF, 设计基于模糊遗忘因子的无迹卡尔曼滤波估计算法, 以提高滤波算法的跟踪性. 在TCS设计方面, 本文根据车辆行驶状态和路面状态, 应用可拓理论, 将基于附着系数和基于车辆加速度计算的基础扭矩进行可拓融合得到目标基础扭矩, 之后以实际滑转率和目标滑转率之间的误差作为输入, 设计模糊自整定PI控制器得到反馈扭矩, 控制驱动轮滑转率接近其目标滑转率, 以避免因附着系数估计不准确而导致TCS控制变差的问题.
1. 车辆和轮胎模型
1.1 五自由度整车动力学模型
本文研究的牵引力控制只考虑车辆水平直线行驶的情况, 不涉及车辆横向运动和横摆运动, 试验车型为4 × 2后驱牵引车. 因此, 考虑模型精度和计算量, 本文采用五自由度整车动力学模型来描述车辆运动状况, 分别为车身的纵向运动和四个车轮的转动. 五自由度整车动力学模型简图如图1所示, OX正方向为车辆行驶方向, 其中O为车辆质心位置, 车轮顺时针旋转方向定为正方向, 下标fl, fr, rl, rr分别表示车辆的左前轮、右前轮、左后轮和右后轮.
整车的纵向运动方程
$$ m{a_x} = {F_{x,fl}} + {F_{x,fr}} + {F_{x,rl}} + {F_{x,rr}} + {F_w} $$ (1) 式中, m为车辆质量, 单位为kg,
$ {a_x} $ 为纵向加速度, 单位为$ {\text{m/}}{{\text{s}}^{\text{2}}} $ ,$ {F_w} $ 为空气阻力;$ {F_{x,ij}} $ 为各轮处纵向轮胎力, 单位为N.对于空气阻力
$ {F_w} $ , 计算公式如下$$ {F_w} = \frac{{{C_D}Av_x^2}}{{21.15}} $$ (2) 式中,
$ {C_D} $ 为空气阻力系数, A为车辆迎风面积, 单位为$ {{\text{m}}^{\text{2}}} $ ,$ {v_x} $ 为纵向车速, 单位为km/h.车轮旋转运动微分方程
$$ \left. {\begin{array}{*{20}{l}} {{I_{w,fl}}{{\dot \omega }_{fl}} = {T_{b,fl}} - {F_{x,fl}}{R_e}} \\ {{I_{w,fr}}{{\dot \omega }_{fr}} = {T_{b,fr}} - {F_{x,fr}}{R_e}} \\ {{I_{w,rl}}{{\dot \omega }_{rl}} = {T_{d,rl}} + {T_{b,rl}} - {F_{x,rl}}{R_e}} \\ {{I_{w,rr}}{{\dot \omega }_{rr}} = {T_{d,rr}} + {T_{b,rr}} - {F_{x,rr}}{R_e}} \end{array}} \right\} $$ (3) 式中,
$ {I_w} $ 为车轮转动惯量, 单位为$ {\text{kg}} \cdot {{\text{m}}^{\text{2}}} $ ,$ \omega $ 为车轮转速, 单位为$ {\text{rad/s}} $ ,$ {T_b} $ 为作用于车轮的制动力矩, 单位为$ {\text{N}} \cdot {\text{m}} $ ,$ {T_d} $ 为作用于车轮的驱动力矩,$ {R_e} $ 为车轮有效转动半径, 单位为m.进一步地, 当车轮未打滑或滑转率较小时, 有
$$ {F_t} = {F_j} + {F_f} + {F_w} $$ (4) 式中,
$ {F_t} $ 为驱动力,$ {F_j} $ 为加速阻力,$ {F_f} $ 为滚动阻力.对于
$ {F_t} $ ,$ {F_j} $ ,$ {F_f} $ , 计算公式如下$$ {F_t} = \frac{{{T_e}{i_g}{i_0}{\eta _T}}}{{{R_e}}} $$ (5) $$ {F_j} = \delta m{a_x} $$ (6) $$ {F_f} = mgf $$ (7) 式中,
$ {T_e} $ 为发动机扭矩,$ {i_g} $ 为变速器当前传动比,$ {i_0} $ 为主减速器传动比,$ {\eta _T} $ 为传动效率,$ \delta $ 为车辆旋转质量换算系数, f为滚动阻力系数.考虑加速、制动过程中的纵向载荷转移, 各个车轮所受到的地面法向反作用力近似计算公式如下
$$ \left. {\begin{array}{*{20}{l}} {{F_{z,fl}} = \dfrac{{mg}}{2}\left( {\dfrac{{{L_r}}}{L} - \dfrac{{{h_g}{a_x}}}{{gL}}} \right)} \\ {{F_{z,fr}} = \dfrac{{mg}}{2}{{\left( {\dfrac{{{L_r}}}{L} - \dfrac{{{h_g}{a_x}}}{{gL}}} \right)} }} \\ {{F_{z,rl}} = \dfrac{{mg}}{2}\left( {\dfrac{{{L_f}}}{L} + \dfrac{{{h_g}{a_x}}}{{gL}}} \right)} \\ {{F_{z,rr}} = \dfrac{{mg}}{2}\left( {\dfrac{{{L_f}}}{L} + \dfrac{{{h_g}{a_x}}}{{gL}}} \right)} \end{array}} \right\} $$ (8) 式中, L为轴距,
$ {L_f} $ 为车辆质心到前轴距离,$ {L_r} $ 为车辆质心到后轴的距离,$ {h_g} $ 为质心高度.1.2 轮胎模型
Dugoff轮胎模型[22]是一种应用广泛的半经验轮胎公式, 具有近似高精度的特点. 为了准确地表达轮胎的非线性力学特性, 该轮胎模型引入了边界值L对模型进行了修正, 从而可以得到描述轮胎纵向力和侧向力与路面附着系数之间的函数表达式. Dugoff轮胎模型的输入输出如图2所示.
对于单个轮胎受力分析, 在轮胎坐标系下经过受力分解得到纵向力与侧向力的表达式如下
$$ {F_x} = \mu {F_{z,kN}}{C_x}\frac{\lambda }{{1 - \lambda }}f(L) $$ (9) $$ {F_y} = \mu {F_{z,kN}}{C_y}\frac{{\tan \alpha }}{{1 - \lambda }}f(L) $$ (10) 其中
$$ f(L)= \begin{cases}L(2-L) & L < 1 \\ 1 & L \geqslant 1\end{cases} $$ (11) $$ L = \frac{{(1 - \lambda )\left( {1 - \varepsilon {v_x}\sqrt {C_x^2{\lambda ^2} + C_y^2{{\tan }^2}\alpha } } \right)}}{{2\sqrt {C_x^2{\lambda ^2} + C_y^2{{\tan }^2}\alpha } }} $$ (12) 式中,
${F_{z,kN}} = \dfrac{{{F_z}}}{{1000}}$ ,$ \mu $ 为路面附着系数,$ \varepsilon $ 为速度影响系数, 可通过仿真和实车试验标定得到, L为边界值, 主要用于轮胎非线性区间轮胎力的拟合;$ {C_x} $ 和$ {C_y} $ 分别为轮胎的纵向刚度和侧偏刚度, 其与轮胎结构、材料和轮胎垂直载荷有关. 通过仿真与实车试验, 本文得到轮胎参数$ {C_x} = 76.8\;{\text{kN/rad}} $ ,$ {C_y} = 103\;{\text{kN/rad}} $ ,$ \varepsilon = 1.16. $ 由Dugoff轮胎模型公式可知, 路面附着系数与公式中其他参数没有关联, 因此可将Dugoff轮胎模型中的路面附着系数分离出来, 式(9)和式(10)修改如下
$$ {F_x} = \mu F_x^0 = \mu {F_{z,kN}}{C_x}\frac{\lambda }{{1 - \lambda }}f(L) $$ (13) $$ {F_y} = {\mu _r}F_y^0 = \mu {F_{z,kN}}{C_y}\frac{{\tan \alpha }}{{1 - \lambda }}f(L) $$ (14) 式中,
$ F_x^0 $ 和$ F_y^0 $ 分别为轮胎在纵向和侧向受力的归一化表示, 其与路面附着系数无关. 将路面附着系数分离出来, 采用归一化的表达形式有利于路面附着系数估计算法中各系数或表达式的确定, 减小计算量, 也提高计算的准确度.2. 路面附着系数估计
2.1 基于模糊遗忘因子的无迹卡尔曼滤波
传统路面附着系数算法常常选取需要估计的状态变量为四个车轮处路面附着系数, 而对于本文所研究的后驱商用车TCS来说, 其控制对象为驱动轮, 同时车辆处在均一和对开路面需要采取不同的控制策略, 所以状态向量选取左右两侧驱动轮下的路面附着系数
$ {\mu _l} $ 和$ {\mu _r} $ , 这样可以简化模型, 减小计算量, 同时也满足TCS控制的设计要求. 车辆的纵向加速度和左、右驱动轮轮速对路面附着系数比较敏感, 因此, 量测变量选取$ {a_x} $ ,$ {\omega _{rl}} $ 和$ {\omega _{rr}} $ . 本文采用的基于离散非线性系统的无迹卡尔曼滤波, 其实现步骤较长, 由于篇幅限制, 在此不再赘述.对于时不变系统, 传统无迹卡尔曼滤波能够实现较高的估计精度[23], 但随着无迹卡尔曼滤波过程的进行, 获得的数据不断增加, 旧数据的累积对协方差矩阵P的影响逐渐增大, 使之趋于稳定, 这导致无迹卡尔曼滤波无法快速准确地跟踪参数变化. 路面附着系数是时变参数, 路面附着系数随时可能发生突变, 而在这时新数据比老数据更能反应参数变化的情况. 在估计算法中引入遗忘因子可以用来确定数据更新的权重, 算法的收敛速度和跟踪速度受遗忘因子的影响[24]. 对于常遗忘因子f来说, 其主要用于决定新老数据的权重, f越小, 加强新数据弱化老数据的程度就越强, 系统的跟踪能力越强, 但是对噪音越敏感, 收敛性较差; f越大, 加强新数据弱化老数据的程度相对较弱, 系统跟踪能力就越弱, 对噪声越不敏感, 收敛性较好[25]. 因此常遗忘因子难以快速且准确地追踪路面附着系数这种时变参数. 如果能够综合f在较大时的收敛性和f在较小时的跟踪性, 那算法的性能就能得到很大的提升. 于是本文引入了模糊遗忘因子的概念, 通过模糊控制动态地调整遗忘因子的大小, 使估计算法在不同情况下都能有较好的收敛性和跟踪能力.
在路面附着系数发生突变时, 实际量测向量
$ {\boldsymbol{z}} = [{a_x}\; {\omega _{rl}} \;{\omega _{rr}}] $ 会产生明显变化, 实际残差协方差矩阵$ {{\boldsymbol{r}}_k} $ 也会产生明显变化, 而旧数据的累积使得无迹卡尔曼滤波理论差协方差矩阵$ {{\boldsymbol{P}}_{{z_k}{z_k}}} $ 趋于稳定, 这种情况导致的算法跟踪能力变差, 能够通过e (理论差协方差矩阵与残差协方差矩阵之差矩阵的迹)很好地体现出来. 但是三个量测变量$ {a_x} $ ,$ {\omega _{rl}} $ ,$ {\omega _{rr}} $ 的量纲不一样, 不能直接将方差的偏差相加得到e, 因此采用归一化处理对e进行修改. 本文采用最大−最小归一化方法对方差进行处理, 将方差归一化到区间[0, 1]上.通过仿真和实车试验数据分析可知方差
$ D\left( {{a_x}} \right) $ ,$ D\left( {{\omega _{rl}}} \right) $ ,$ D\left( {{\omega _{rr}}} \right) $ 最小值为0, 最大值分别为$ D{\left( {{a_x}} \right)_{max}} = 9 $ ,$D{\left( {{\omega _{rl}}} \right)_{max}} = 400,$ $ D{\left( {{\omega _{rr}}} \right)_{max}} = 400 $ , 从而得到归一化矩阵N, 表示如下$$ {\boldsymbol{N}} = \left[ {\begin{array}{*{20}{c}} {\dfrac{1}{{D{{\left( {{a_x}} \right)}_{max}}}}}&{}&{} \\ {}&{\dfrac{1}{{D{{\left( {{\omega _{rl}}} \right)}_{max}}}}}&{} \\ {}&{}&{\dfrac{1}{{D{{\left( {{\omega _{rr}}} \right)}_{max}}}}} \end{array}} \right] $$ (15) 得到归一化后的e, 表达式如下
$$ {e_k} = {\rm{tr}}\left[ {\left( {{{\boldsymbol{P}}_{{z_k}{z_k}}} - {{\boldsymbol{r}}_k}} \right){\boldsymbol{N}}} \right] $$ (16) 对遗忘因子f模糊控制调节, 具体设计如下.
输入量e的模糊论域设定为[−3, 3], ec的模糊论域为[−10,10], e和ec的模糊子集均为{NB, NM, ZO, PM, PB}. 输出量f的模糊论域为[0.95,1], 模糊子集为{LL, LM, M, BM, BB}. 模糊控制的隶属度函数如图3所示.
本文以e及其变化率ec作为模糊控制的输入, 输出为遗忘因子f. 具体步骤如下.
①输入/输出量的模糊化.
②模糊推理规则设计
当
$ |e| $ 较大时, 应取较小的f, 以提高滤波算法对系统参数变化的跟踪性能.当
$ |e| $ 较小时, 应取较大的f, 以提高滤波算法的收敛性能和抗噪性能.当
$ec \cdot e < 0$ 时,$ |e| $ 有减小的趋势, 故此时$ |ec| $ 越大应取较大的f.当
$ec \cdot e > 0$ 时,$ |e| $ 有增大的趋势, 故此时$ |ec| $ 越大应取较小的f.根据以上结论设计模糊规则表如表1所示.
表 1 模糊推理则表Table 1. Fuzzy inference tableNB NM ZO PM PB NB LL LM BM BM M NM LM LM BM M M ZO LM M BB M LM PM M M BM LM LM PB M BM BM LM LL ③解模糊
利用MATLAB的fuzzy工具, 选择重心法进行解模糊. 常用的解模糊法有最大隶属度法、加权平均法和重心法, 由于重心法具有更加平滑的输出推理控制[26], 因此采用所示重心法[27]进行解模糊得到遗忘因子f, 从而获得修正的协方差矩阵为
$$ \begin{split} &{{\boldsymbol{P}}_{k + 1\mid k}} = \hfill \\ & \frac{1}{f}\left\{ {\left[ {\sum\limits_{i = 1}^{2n} {\omega _c^{(i)}} \left( {{{\widehat x}_{k + 1\mid k}} - x_{k + 1\mid k}^{(i)}} \right){{\left( {{{\widehat x}_{k + 1\mid k}} - x_{k + 1\mid k}^{(i)}} \right)}^{\rm{T}}}} \right] + {\boldsymbol{Q}}} \right\} \end{split} $$ (17) 式中,
$ {\omega _c} $ 为协方差权值,$ \widehat x $ 为状态预测值; Q为W的方差矩阵,$ {\boldsymbol{Q}} = {\boldsymbol{E}}({{\boldsymbol{W}}_k}{{\boldsymbol{W}}_k}^{\rm{T}}) $ , W为系统噪声.2.2 路面附着系数估计仿真分析
为验证路面附着系数估计算法的有效性, 本研究以传统UKF路面附着系数估计算法为对照组, 采用Simulink/TruckSim联合仿真平台进行对接路面和对开路面仿真试验, 对所设计的路面附着系数估计算法的有效性进行验证.
工况一: 对接路面工况
设置路面模型为对接路面, 车辆从静止开始加速, 车辆先进入附着系数为0.75的高附着路面然后进入附着系数为0.3的低附着路面, 路面附着系数估计结果如下图所示.
由图4(a)和图4(b)可以看出, 在估计的初始阶段传统UKF算法收敛速度较慢且估计偏差较大, 在高附着路面上估计最终偏差超过10%, 而FFUKF算法在初始阶段能够在0.5 s内快速收敛到真实值附近, 在高附路面上最终偏差不超过5%. 当路面附着系数突然发生变化时, 传统UKF算法不能快速跟踪真实值, 在路面附着发生突变近3 s后估计偏差仍超过50%. FFUKF算法通过引入模糊遗忘因子能够根据滤波器输出的误差值动态地调整协方差矩阵, 估计算法能够自适应地跟踪路面附着系数的变化, 结果显示在路面附着系数发生变化的0.5 s左右FFUKF算法收敛到真实值附近, 在低附路面最终估计偏差不超过5%.
工况二: 对开路面工况
设置路面模型为对开路面, 左侧为高附侧, 路面附着系数为0.75, 右侧为低附侧, 路面附着系数为0.3, 车辆从静止开始加速, 路面附着系数估计结果如图5所示.
由图5(a)和图5(b)可以看出, 传统UKF算法需要2.5 s左右才能收敛到真实值附近, 且估计误差较大, 左侧附着系数最终误差为10%左右, 右侧附着系数最终误差为20%左右. 而FFUKF估计算法在0.5 s内就能够收敛到真实值附近, 左、右侧附着系数最终估计误差均小于5%.
3. TCS控制器
3.1 TCS基础扭矩
根据上文得到的左右路面附着系数
$ {\mu _l} $ 和$ {\mu _r} $ , 计算路面附着所能提供的最大驱动力$ {F_{tmax }} $ 和所需的发动机净输出扭矩$ {T_{B1}} $ 为$$ \qquad\qquad{F_{tmax }} = {\mu _l}{F_{z,rl}} + {\mu _r}{F_{z,rr}} $$ (18) $$ \qquad\qquad{T_{B1}} = \frac{{{F_{tmax}}{R_e}}}{{{i_0}{i_g}{\eta _T}}} $$ (19) 式中,
$ {i_0} $ 为主减速器传动比,$ {i_g} $ 为变速器传动比,$ {\eta _T} $ 为传动系统传动效率.在一般路面上通过路面附着系数
$ {\mu _l} $ 和$ {\mu _r} $ 计算得到的基础扭矩$ {T_{B1}} $ 较为准确, 但当车辆在不平道路上加速行驶时, 车轮将会产生振动和轮速波动, 这将影响路面附着系数估计效果和基础扭矩$ {T_{B1}} $ 的计算准确度, 从而降低了TCS的控制效果. 因此针对不平路面驱动工况, 本文提出了路面不平参量B的概念, 用以表征车辆在当前路面上行驶时路面附着系数估计的准确度.研究表明路面不平会对轮加速度信号的振动幅值和周期产生影响[28], 因此可以设定加速度门限值, 记录单位时间内驱动轮加速度穿越此门限的次数, 如图6所示, 该方法称为“穿越计次法”[29], 单位时间内的穿越次数即为路面不平参量B, 具体实现步骤如下.
(1) 设定m组穿越门限
$ {a_{ui}} $ 和$ {a_{li}} $ ($ i = 1,2, \cdots ,m $ ,$ {a_u} $ 表示上限,$ {a_l} $ 表示下限).(2) 设定m个不平度累加值
$ {p_i} $ ($ i = 1,2, \cdots ,m $ ), 驱动轮加速度波动的幅度越大, 其能够穿越的门限值越大, 每次穿越门限值所加的不平度累加值$ {p_i} $ 越大.(3) 驱动轮的加速度每穿越一次门限值累加一次
$ {p_i} $ , 记录一定时间T内总的累加值, 记为路面不平参量B, B每过T时间更新一次.由以上可知, 在不平路面上路面附着系数估计不准确, 这导致基于路面附着系数计算得到的基础扭矩
$ {T_{B1}} $ 不准确. 为解决以上问题, 可以驱动轮滑转率在目标滑转率附近时记录此时车辆加速度作为车辆在该路面上的目标加速度, 以此加速度计算目标基础扭矩, 计算公式如下$$ {T_{B2}} = \frac{{\left( {m{a_r} + {F_f} + {F_w}} \right){R_e}}}{{{i_0}{i_g}{\eta _T}}} $$ (20) 式中,
$ {a_r} $ 为目标加速度.此方法计算的基础扭矩
$ {T_{B2}} $ 受路面不平影响较小, 计算较为准确. 但在路面附着系数发生突变, 驱动轮滑转率可能不在目标滑转率附近时, 此时目标加速度保持上一次的值, 没有更新, 导致基础扭矩$ {T_{B2}} $ 计算不准确. 针对以上问题, 本文基于可拓控制理论划分经典域、可拓域和非域, 通过可拓关联函数得到动态权重系数, 将上述两种方法计算得到的目标基础扭矩$ {T_{B1}} $ 和$ {T_{B2}} $ 进行可拓融合, 以提高目标扭矩的准确度. 具体实现步骤如下.(1) 特征量的选取
基于上文分析可将路面不平参量B作为一个特征量. 另外, 当车辆起步, 车速较低时, 路面附着系数估计算法还未收敛, 估计不准确. 同时, 在低速时轮速传感器误差影响较大, 这也会使估计不准确, 导致
$ {T_{B1}} $ 的计算也不准确. 而基于轮速加速度计算得到的$ {T_{B2}} $ 在低速时较为准确, 且反应迅速. 故同时选择车速v作为另一特征量, 为方便可拓集的划分, 将特征量重新定义为$ {v'} = \dfrac{1}{v} $ . 基于以上考虑本文选取的特征量为$ S(R,{v'}) $ .(2) 可拓集合划分
可拓集合划分如图7所示. 根据仿真和实车调试选取路面不平度的经典域和可拓域的边界值
$ {B_1} = 20 $ , 可拓域与非域的边界值$ {B_2} = 40 $ .$ {v'} $ 的经典域和可拓域的边界值$ v_1' = 0.4 $ , 可拓域与非域的边界值$ v_2' = 0.7 $ .(3) 关联函数计算
$ {P_0} $ 表示当前特征量的值. 将$ {P_0} $ 与原点O相连并延长, 直线$ O{P_0} $ 与经典域边界交于$ {P_1} $ , 与可拓域边界交于$ {P_2} $ . 根据可拓理论定义,$ {P_0} $ 到经典域与可拓域的可拓距为$ d\left( {{P_0}, \langle O,{P_1} \rangle } \right) $ 和$ d\left( {{P_0}, \langle O,{P_2} \rangle } \right) $ . 可拓距和关联函数的计算公式如下, 其中i = 1, 2$$ \begin{gathered} {\text{ }}d\left( {{P_0}, \langle O,{P_i} \rangle } \right) = \left\{ {\begin{array}{*{20}{l}} { - \left| {O{P_0}} \right|},&{{P_0} \in \langle O,{P_i}m \rangle } \\ { - \left| {{P_0}{P_i}} \right|},&{{P_0} \in \langle {P_1}m,{P_i} \rangle } \\ {\left| {{P_0}{P_i}} \right|},&{{P_0} \in \langle {P_i}, + \infty \rangle } \end{array}} \right. \end{gathered} $$ (21) $$ K(S) = \frac{{d \left( {P_0},\langle {O,{P_2}} \rangle \right) }}{{D\left( {{P_0}, \langle {P_O},{P_2} \rangle ,{P_0}, \langle O,{P_1} \rangle } \right)}} $$ (22) $$ \begin{split}& D\left( {{P_0}, \langle {P_1},{P_2} \rangle ,{P_0}, \langle O,{P_1} \rangle } \right) = \hfill \\ & \qquad d\left( {{P_0}, \langle O,{P_2} \rangle } \right) - \left( {{P_0}, \langle O,{P_1} \rangle } \right) \end{split} $$ (23) 式中,
$ {P_1}m $ 为线段$ O{P_1} $ 的中点,$ {P_2}m $ 为$ O{P_2} $ 的中点.(4) 测度模式判别
测度模式
$ {M_1} $ : 当$ K(S) \geqslant 1 $ 时, 特征量$ S(R,v') $ 在经典域内, 路面附着系数估计较为准确,$ {T_{B1}} $ 计算准确, 所以权重系数$ \xi = 1 $ .测度模式
$ {M_2} $ : 当$ 0 < K(S) < 1 $ 时, 特征量$ S(R,v') $ 在可拓域内, 权重系数$ \xi = K(S) $ .测度模式
$ {M_3} $ : 当$ K(S) \leqslant 0 $ 时, 特征量$ S(R,v') $ 在非域内, 路面附着系数估计有较大误差,$ {T_{B1}} $ 计算不准确, 此时基础扭矩采用$ {T_{B2}} $ , 权重系数$ \xi = 0 $ .(5) 基础扭矩计算
根据可拓融合理论, 得到的基础扭矩为
$$ {T_B} = \xi {T_{B1}} + (1 - \xi ){T_{B2}} $$ (24) 3.2 TCS反馈扭矩
TCS的设计目标一般是将驱动轮的滑转率控制在0.1~0.2之间, 然而当车辆起步时, 由于此时车速很小, 轻微的轮速波动会导致滑转率的值出现很大的变化. 本文目标驱动轮轮速的设计思路为: 车速较低时维持一定滑转量, 随着车速增长不断减小滑转量, 车速达到一定程度时, 再将驱动轮滑转率维持在10%~20%范围内. 基于上述分析得到目标驱动轮轮速计算公式如下
$$ \left. {\begin{array}{*{20}{l}} {{v_{tar{\text{ }}}} = 5 + {v_x},\quad {\text{ }}{v_x} < 5} \\ {{v_{tar{\text{ }}}} = 8 + 0.4{v_x},\quad {\text{ }}5 \leqslant {v_x} < 10} \\ {{v_{tar{\text{ }}}} = 2 + {v_x},\quad {\text{ }}10 \leqslant {v_x} < 20} \\ {{v_{tar{\text{ }}}} = 1.1{v_x},\quad {\text{ }}{v_x} \geqslant 20} \end{array}} \right\}$$ (25) 式中,
$ {v_{tar{\text{ }}}} $ 为目标驱动轮轮速.在PID控制中, 微分控制项需要误差的微分, 而在实际控制系统中误差的微分往往波动很大, 并且会有突变的现象, 如果对误差的微分做滤波处理又会产生延迟, 失去了微分控制项能够进行“预测控制”的效果. 因此在实际工程中有时会摒弃微分控制项, 转而采用前馈 + PI控制器, 在上节已经介绍了基于前馈思想的基础扭矩设计, 因此, PI控制器更适合本文中反馈扭矩的设计. 传统PI算法参数固定, 难以适应复杂多变的情况, 模糊自整定控制技术是具有自学习功能的模糊控制类技术[30], 模糊自整定PID控制能通过使用模糊推理实时对PID控制器的参数进行自整定[31], 算法的适应性强.
模糊自整定PI控制算法中模糊控制以驱动轮轮速误差的绝对值
$ |E| $ 为输入, PI控制的两个参数$ {k_p} $ ,$ {k_i} $ 为输出. E为驱动轮目标轮速与实际轮速之差, 左、右驱动轮轮速的平均值作为实际驱动轮轮速$ {v_{act}} $ ,$ {v_{act}} $ 和E的计算如下$$ \qquad\qquad{v_{act}} = \frac{{{\omega _{rl}}{R_e} + {\omega _{rr}}{R_e}}}{2} $$ (26) $$ \qquad\qquad E = {v_{tar{\text{ }}}} - {v_{act}} $$ (27) 模糊自整定PI控制具体实现步骤如下.
(1)输入/输出量的模糊化
输入量
$ |E| $ 的模糊论域设定为[0,10], 模糊子集为{LL, LM, M, BM, BB}. 输出量$ {k_p} $ 的论域为[400, 600], 模糊子集为{LL, LM, M, BM, BB}. 输出量$ {k_i} $ 的论域为[100, 200], 模糊子集为{LL, LM, M, BM, BB}. 模糊控制的隶属度函数如图8所示.(2) 模糊推理规则设计
根据参数
$ {k_p} $ 和$ {k_i} $ 的作用特点, 设计模糊规则如表2所示.表 2 模糊自整定PI的模糊规则Table 2. Fuzzy rules of fuzzy adaptive PI$ |E| $ LL LM M BM BB $ {k_p} $ LL LM M BM BB $ {k_i} $ BB BM M LM LL (3) 解模糊
采用重心法进行解模糊得到
$ {k_p} $ 和$ {k_i} $ .(4) 反馈扭矩计算
根据以上得到的PI控制参数
$ {k_p} $ ,$ {k_i} $ 计算反馈扭矩, 计算公式如下$$ \qquad\qquad{T_{w,FB}} = {k_p}E(t) + {k_i}\int_0^t E (t){\rm{d}}t $$ (28) $$ \qquad\qquad{T_{e,FB}} = \frac{{{T_{w,FB}}}}{{{i_0}{i_g}{\eta _T}}} $$ (29) 式中,
$ {T_{w,FB}} $ 为轮端反馈扭矩,$ {T_{e,FB}} $ 为发动机端反馈扭矩.3.3 制动控制器
本文所设计的制动控制结构框图如图9所示. TCS制动压力控制在均一和对开路面上采取不同控制策略. 在均一路面上, 制动控制的主要作用是保持左右驱动轮轮速一致, 因此在均一路面上以左右驱动轮轮速差为输入采用PI压力控制得到
$ {P_P} $ , 适时的压力控制保证左右驱动轮轮速保持一致. 而在对开路面上, 由于平衡式差速器平分扭矩到两驱动轮上的特点, 这会导致低附侧驱动轮发生较大滑转, 而高附侧车轮几乎不滑转, 无法充分利用高附侧附着力, 因此需要对低附侧驱动轮施加克服两侧附着差的附着差压力$ {P_\mu } $ , 保证高附着车轮得到充足动力性, 同时以左右驱动轮轮速差为输入采用PI压力控制得到的$ {P_P} $ 作为补偿制动压力.(1)压力PI控制
在所设计的PI控制器中, 输入量为两驱动转速度差, 以两驱动转速度差等于零为控制目标, 计算公式如下
$$ \qquad\qquad{e_v} = {\omega _{rr}}{R_e} - {\omega _{rl}}{R_e} $$ (30) $$\qquad\qquad {P_P} = {K_p}{e_\omega }(t) + {K_i}\int_0^t {{e_\omega }} (t){\rm{d}}t $$ (31) 式中,
$ {K_p} $ 和$ {K_i} $ 为PI控制器的参数.当
$ {P_\omega } > 0 $ 时, 说明右侧驱动轮速比左侧快, 故对右侧驱动轮施加制动压力$ {P_\omega } $ , 当$ {P_\omega } < 0 $ 时, 说明左侧驱动轮速比右侧快, 故对左侧驱动轮施加制动压力$ - {P_\omega } $ .(2)附着差压力
对开路面驱动轮受力情况如图10所示, 假设左侧为高附侧, 右侧为低附侧. 当TCS扭矩控制介入后, 由于降扭控制的作用, 高附侧驱动轮不会出现过度滑转的情况, 而低附侧驱动轮可能还在滑转, 此时高附侧和低附侧驱动轮能够达到最大的轮端驱动力矩, 分别为
$$ {T_{w,rl}} = \frac{{{T_{e,act}}{i_0}{i_g}{\eta _T}}}{2} $$ (32) $$ \begin{split} &\\ & {T_{w,rr}} = {\mu _r}{F_{z,rr}}{R_e} \end{split}$$ (33) 式中,
$ {T_{e,act}} $ 为当前发动机实际扭矩.从图10可以看出, 要想高附侧驱动轮达到最大驱动力, 充分利用高附侧驱动力, 就必须对低附侧驱动轮施加
$ \Delta {T_w} $ 的制动力矩, 以克服平衡式差速器的作用, 以此制动力矩计算得到的制动压力即为附着差压力, 再通过两种典型路面判断即可知道制动力需要作用于哪侧驱动轮上. 当然, 施加$ \Delta {T_w} $ 的制动力矩是理论值, 在实际情况中可能会因为路面附着估计误差以及发动机扭矩信号误差等导致$ \Delta {T_w} $ 计算不准确. 因此, 为了防止过度制动引起的驱动能力和舒适性降低, 将在车轮上施加$ 0.7 \Delta {T_w} $ 的制动力矩, 剩下的制动力可通过PI控制压力补偿. 附着差压力计算公式如下$$ \Delta {T_w} = {T_{w,rl}} - {T_{w,rr}} $$ (34) $$ {P_\mu } = \frac{{0.7 \Delta {T_w}}}{{{K_b}}}\quad\quad $$ (35) 式中,
$ \Delta {T_w} $ 为附着差引起的两侧驱动轮扭矩差,$ {P_\mu } $ 为附着差压力,$ {K_b} $ 为制动器系数.需要说明的是, 为了防止制动控制的过度介入, 在均一路面上制动控制介入的门限设置为15 kPa. 另外, 制动控制的精度为10 kPa.
4. 实车试验
本次实车试验场地为黑龙江省的北安市极地汽车高寒试验场, 设计了均一低附路面、对接路面和对开路面这三种路面试验工况, 实车试验所选择的车辆为解放JH6牵引车头(长宽高为6.165 m × 2.495 m × 3.91 m), 轴数为2, 后驱4 × 2, 车辆相关配置参数如表3所示, 其中
$ {P_m} $ 为发动机最大输出功率,$ {T_m} $ 为发动机最大扭矩. 试验过程中开启TCS, 驾驶员快速踩下加速踏板至100%开度并保持, 车辆由静止开始直线加速行驶, 换挡逻辑为车辆自身变速器换挡逻辑, 从一挡开始自动换挡.表 3 试验车辆参数配置Table 3. Vehicle parametersParameter Value m/kg 6800 L/m 3.9 $ {L_f} $/m 1.49 $ {L_r} $/m 2.485 $ {h_g} $/m 1.123 $ {P_m} $/kW 370 $ {T_m} $/(N·m) 2300 4.1 均一路面试验
均一低附路面实车试验结果如图11所示. 从图11(a)可以看出, 在车辆速度大于0后, 估计算法开始进行估计, 本文所设计的FFUKF路面附着系数估计算法的估计值能够在1 s内快速收敛到0.2左右并保持稳定, 其估计值在冰面路面附着系数参考值为0.1~0.25内, 说明本文设计的算法能够快速和准确地对路面附着系数进行估计, 为TCS控制算法提供了有力的支持, 提高了TCS控制效果.
结合图11(b)~图11(d)可以看出, 在5 s左右车辆开始起步, 此时驱动轮还未出现过度滑转, TCS未介入, TCS目标扭矩为驾驶员期望扭矩, 2 s后右后驱动轮出现过度滑转, 7 s左右TCS系统开始介入, 立即进行了降扭控制, 驱动轮滑转程度迅速减小. 同时, 制动控制适时进行单侧驱动轮制动, 使左右两侧驱动轮轮速保持一致. 在10 s到20 s, 30 s, 40 s和50 s左右驱动轮轮速并不能很好地收敛到目标值附近, 这是由于在起步加速阶段换挡导致的, 商用车发动机扭矩大, 离合器分离时驱动轮失去动力导致驱动轮轮速下降至车速附近, 离合器结合时产生冲击造成驱动轮轮速突然增大和波动, 但是在整个控制过程中驱动轮并未出现过度滑转和过大的轮速波动. 除去换挡阶段, TCS能够将驱动轮轮速控制在目标值附近, 整个过程中平均加速度为0.4
$ {\text{m/}}{{\text{s}}^{\text{2}}} $ 左右, 车辆有较强的加速能力, TCS控制效果良好.4.2 对接路面试验
对接路面实车试验结果如图12所示. 从图12(a)可以看出, 在车辆速度大于0后, FFUKF估计算法开始进行路面附着系数估计, 在1 s内左右两侧附着系数估计值达到0.5左右, 7.5 s时收敛到了0.7左右, 由于车辆在高附路面时驱动轮滑转率非常小且车辆进行换挡导致识别速度较慢, 准确度不够高. 但是本文路面附着系数算法主要用于TCS控制, 在驱动轮滑转率非常小时TCS未介入, 另外TCS基础扭矩通过可拓融合得到, 并非只依赖于路面附着系数, 因此这对TCS控制影响较小. 在10 s左右车辆驶入低附路面, FFUKF算法能够快速准确地跟踪路面附着系数的变化, 2 s内快速收敛到冰面附着系数参考值附近且保持稳定.
结合图12(b)~图12(d)可以看出, 在2 s左右车辆开始起步, 驾驶员期望扭矩快速上升, 10 s之前车辆行驶在高附路面上, 驱动轮并未出现过度滑转, TCS未介入, TCS目标扭矩为驾驶员期望扭矩. 在10 s左右车辆由高附路面行驶至低附路面, 驱动轮出现过度滑转, 此时TCS迅速介入, 立即进行了降扭控制, 在13 s左右驱动轮轮速被控制在目标轮速附近. 同样地, 在10 s, 17 s和18 s左右由于换挡的原因出现了轮速的突然增大, 在TCS的控制下驱动轮轮速很快降了下去, 但是由于过度降扭导致17~19 s左右驱动轮轮速并不能很好地收敛到目标值附近, 这是需要优化的方向之一. 总的来说在对接路面上TCS能够在驱动轮出现过度滑转时迅速介入抑制驱动轮的滑转, 在冰面加速阶段车辆的平均加速度为0.37
$ {\text{m/}}{{\text{s}}^{\text{2}}} $ 左右, 车辆加速能力较强, TCS控制效果较好.4.3 对开路面试验
对开路面实车试验结果如图13所示. 从图13(a)可以看出, 在车辆速度大于0后, 估计算法开始进行估计, 左侧车轮在水泥石板路面上, 本文所设计的FFUKF估计算法在1 s内达到0.4左右, 在之后10 s内收敛到石板路面附着系数参考值0.75~0.85左右, 高附侧估计值收敛速度较慢, 原因是高附侧驱动轮滑转率非常小且车辆进行换挡导致估计速度较慢, 准确度不够高, 这是算法优化方向之一. 右侧车轮在冰面上, 估计值在2 s内收敛到冰面附着系数参考值0.1~0.25附近, 估计快速且准确.
结合图13(a)~图13(c)可以看出, 在2 s左右车辆开始起步, 驾驶员期望扭矩快速上升, 此时驱动轮还未出现过度滑转, TCS未介入, TCS目标扭矩为驾驶员期望扭矩, 之后2 s内实际扭矩增大, 右后驱动轮出现过度滑转, 4 s左右TCS开始介入, 立即进行了降扭控制, 驱动轮滑转程度迅速减小, TCS制动控制持续对右侧车轮施加制动压力, 以保证高附侧驱动轮能够充分利用路面附着, 当进行换挡时, 两侧驱动轮失去了动力, 故此时不进行制动. 综合来看, TCS通过扭矩控制和对低附侧车轮进行制动控制抑制了驱动轮的过度滑转, 整个过程中车辆的平均加速度为0.6
$ {\text{m/}}{{\text{s}}^{\text{2}}} $ 左右, 提高了车辆的动力性, 控制效果较好.5. 结 论
本文以后驱牵引车为研究对象, 设计了基于路面附着系数估计的牵引力控制系统.
(1)在影响牵引力控制的路面附着系数估计方面, 设计了基于遗忘因子的无迹卡尔曼滤波估计算法, 通过模糊控制动态地调整遗忘因子的大小, 解决了传统无迹卡尔曼滤波无法快速跟踪时变非线性系统的问题.
(2) TCS扭矩控制器采用基础扭矩加反馈扭矩控制的设计思路. 在TCS基础扭矩设计方面, 通过可拓融合将两种计算基础扭矩的方法进行加权融合计算, 提高了基础扭矩计算的准确性.
(3) TCS制动控制器在均一和对开路面上采取不同的压力控制策略. 在均一路面上, 以两侧驱动轮轮速差为输入利用PI控制算法得到均一路面目标制动压力, 控制左右两侧驱动轮轮速保持一致; 在对开路面上, 利用路面附着系数和发动机当前扭矩计算附着差压力, 以两侧驱动轮轮速差为输入利用PI控制计算补偿压力, 两者相加得到对开路面目标制动压力. 该方法提高了制动控制的准确性, 防止制动控制过渡介入, 提高了舒适性.
(4) 实车试验结果表明路面附着系数估计算法在不同路面状态下收敛迅速, 估计结果较准确, 满足TCS控制使用要求, 提高了TCS控制效果; 在三种路面状态下TCS控制都能够有效地提高车辆的动力性和横向稳定性, 控制效果良好.
-
表 1 模糊推理则表
Table 1 Fuzzy inference table
NB NM ZO PM PB NB LL LM BM BM M NM LM LM BM M M ZO LM M BB M LM PM M M BM LM LM PB M BM BM LM LL 表 2 模糊自整定PI的模糊规则
Table 2 Fuzzy rules of fuzzy adaptive PI
$ |E| $ LL LM M BM BB $ {k_p} $ LL LM M BM BB $ {k_i} $ BB BM M LM LL 表 3 试验车辆参数配置
Table 3 Vehicle parameters
Parameter Value m/kg 6800 L/m 3.9 $ {L_f} $/m 1.49 $ {L_r} $/m 2.485 $ {h_g} $/m 1.123 $ {P_m} $/kW 370 $ {T_m} $/(N·m) 2300 -
[1] 李贵强. 基于自抗扰算法的电动汽车牵引力系统优化控制技术研究. [博士论文]. 北京: 中国科学院大学, 2020 Li Guiqiang. Study on optimization of electric vehicle traction control system base on ADRC algorithm. [PhD Thesis]. Beijing: University of Chinese Academy of Sciences, 2020 (in Chinese)
[2] 冉旭. 牵引力控制系统中压力与扭矩的协调控制. [硕士论文]. 秦皇岛: 燕山大学, 2013 Ran Xu. Coordinated control based on brake pressure and engine torque in traction control system. [Master Thesis]. Qinhuangdao: Yanshan University, 2013 (in Chinese)
[3] 林棻, 王少博, 赵又群等. 基于改进Keras模型的路面附着系数估计. 机械工程学报, 2021, 57(12): 74-86 (Lin Fen, Wang Shaobo, Zhao Youqun, et al. Road friction coefficient estimation based on improved Keras model. Journal of Mechanical Engineering, 2021, 57(12): 74-86 (in Chinese) Lin Fen, Wang Shaobo, Zhao Youqun et al. Road friction coefficient estimation based on improved Keras model. Journal of Mechanical Engineering, 2021, 57(12): 74-86 (in Chinese)
[4] 李韶华, 王桂洋, 杨泽坤等. 基于DRBF-EKF算法的车辆质心侧偏角与路面附着系数动态联合估计. 力学学报, 2022, 54(7): 1-13 (Li Shaohua, Wang Guiyang, Yang Zekun, et al. Dynamic joint estimation of vehicle sideslip angle and road adhesion coefficient based on drbf-ekf algorithm. Chinese Journal of Theoretical and Applied Mechanics, 2022, 54(7): 1-13 (in Chinese) Li Shaohua, Wang Guiyang, Yang Zekun et al. Dynamic joint estimation of vehicle sideslip angle and road adhesion coefficient based on drbf-ekf algorithm. Chinese Journal of Theoretical and Applied Mechanics, 2022, 54(7): 1-13 (in Chinese)).
[5] 胡敬宇, 汪, 严永俊等. 基于限定记忆随机加权扩展卡尔曼滤波的车辆状态估计. 东南大学学报, 2022, 52(2): 387-393 (Hu Jingyu, Wang Yan, Yan Yongjun, et al. Vehicle state estimation based on limited memory random weighted extended Kalman filter. Journal of Southeast University, 2022, 52(2): 387-393 (in Chinese) Hu Jingyu, Wang Yan, Yan Yongjun et al. Vehicle state estimation based on limited memory random weighted extended Kalman filter. Journal of Southeast University, 2022, 52(2): 387-393 (in Chinese)
[6] 王博, 卢萍萍, 管欣等. 路面附着系数识别方法发展现状综述. 汽车技术, 2014, 8: 1-7 (Wang Bo, Lu Pingping, Guan Xin, et al. A review on the development status of road adhesion coefficient identification approach. Automobile Technology, 2014, 8: 1-7 (in Chinese) doi: 10.3969/j.issn.1000-3703.2014.08.001 (Wang Bo, Lu Pingping, Guan Xin et al. A review on the development status of road adhesion coefficient identification approach. Automobile Technology, 2014 (8): 1-7 (in Chinese)). doi: 10.3969/j.issn.1000-3703.2014.08.001
[7] Sato Y, Kobayashi D, Watanabe K, et al. Study on recognition method for road friction condition. Transactions of the Society of Automotive Engineers of Japan, 2007, 38(2): 51-56
[8] Koskinen S. Sensor data fusion based estimation of tyre-road friction to enhance collision avoidance. [PhD Thesis]. Tampere University of Technology, 2010
[9] Cho W, Choi J, Kim C, et al. Unified chassis control for the improvement of agility, maneuverability, and lateral stability. IEEE Transactions on vehicular Technology, 2012, 61(3): 1008-1020 doi: 10.1109/TVT.2012.2183152
[10] Sharifzadeh M, Akbari A, Timpone F, et al. Vehicle tyre/road interaction modeling and identification of its parameters using real-time trust-region methods. IFAC-Papers On Line, 2016, 49(3): 111-116
[11] 顾凯峰. 驾驶员行为特性及路面附着系数的车辆纵向主动安全控制研究. [硕士论文]. 重庆: 重庆交通大学, 2021 Gu Kaifeng. Research on vehicle longitudinal active safety control considering driver behavior characteristics and road adhesion coefficient. [Master Thesis]. Chongqing: Chongqing Jiaotong University, 2021 (in Chinese)
[12] Wang Y, Yin G, Dong H. A novel approach for tire-road friction coefficient estimation using adaptive cubature kalman filter//2020 4th CAA International Conference on Vehicular Control and Intelligence (CVCI). Hangzhou, 18-20 Dec. 2020
[13] Liu YH, Li T, Yang YY, et al. Estimation of tire-road friction coefficient based on combined APF-IEKF and iteration algorithm. Mechanical Systems and Signal Processing, 2017, 88: 25-35 doi: 10.1016/j.ymssp.2016.07.024
[14] Pei X, Hu X, Liu W et al. State estimation of vehicle’s dynamic stability based on the nonlinear kalman filter. Automot. Innov., 2018, 1: 281-289 doi: 10.1007/s42154-018-0028-6
[15] Chen L, Bian M, Luo Y, et al. Estimation of road-tire friction with unscented Kalman filter and MSE-weighted fusion based on a modified Dugoff tire model//SAE Technical Paper, 2015-01-1601, 2015
[16] 钱华明, 葛磊, 彭宇. 多渐消因子卡尔曼滤波及其在 SINS 初始对准中的应用. 中国惯性技术学报, 2012, 20(3): 287-291 (Qian Huaming, Ge Lei, Peng Yu. Multiple fading factors Kalman filter and its application in SINS initial alignment. Journal of Chinese Inertial Technology, 2012, 20(3): 287-291 (in Chinese) doi: 10.3969/j.issn.1005-6734.2012.03.008 Qian Huaming, Ge Lei, Peng Yu. Multiple fading factors Kalman filter and its application in SINS initial alignment. Journal of Chinese Inertial Technology, 2012, 20(3): 287-291(in Chinese)). doi: 10.3969/j.issn.1005-6734.2012.03.008
[17] 张浩. 基于强跟踪无迹卡尔曼滤波的结构时变参数识别. [硕士论文]. 兰州: 兰州理工大学, 2016 Zhang Hao. Structural time-varying parameter identification based on strong tracking unscented Kalman filter. [Master Thesis], Lanzhou: Lanzhou University of Technology, 2016 (in Chinese)
[18] 李兴佳, 李建芬, 朱敏等. 基于无迹卡尔曼滤波的定位融合与校验算法研究. 汽车工程, 2021, 43(6): 825-832 (Li Xingjia, Li Jianfen, Zhu Min, et al. Research on positioning fusion and verification algorithm based on UKF. Automotive Engineering, 2021, 43(6): 825-832 (in Chinese) Li Xingjia, Li Jianfen, Zhu Min, et al. Research on positioning fusion and verification algorithm based on UKF. Automotive Engineering, 2021, 43(6): 825-832 (in Chinese)
[19] 张成宝, 吴光强, 丁玉兰等. 汽车驱动防滑的控制方法研究. 汽车工程, 2000, 22(5): 324-328 (Zhang Chengbao, Wu Guangqiang, Ding Yulan, et al. A research on the control method of vehicle anti-slip regulation. Automotive Engineering, 2000, 22(5): 324-328 (in Chinese) (Zhang Chengbao, Wu Guangqiang, Ding Yulan et al. A research on the control method of vehicle anti-slip regulation. Automotive Engineering, 2000, 22(05): 324-328 (in Chinese)).
[20] 田浩. 轮毂电机汽车驱动防滑控制研究. [硕士论文]. 吉林: 吉林大学, 2022 (Tian Hao. Research on acceleration slip regulation system for electric vehicles with in-wheel motors. [Master Thesis]. Jilin: Jilin University, 2022 (in Chinese)
[21] Ding X, Zhang L, Wang Z, et al. Acceleration slip regulation for four-wheel-independently-actuated electric vehicles based on road identification through the fuzzy logic. IFAC-Papers On Line, 2018, 51(31): 943-948 doi: 10.1016/j.ifacol.2018.10.059
[22] Dugoff H, Fancher PS, Segel L. An analysis of tire traction properties and their influence on vehicle dynamic performance. SAE transactions, 1970, 79: 1219-1243
[23] 王震坡, 薛雪, 王亚超. 基于自适应无迹卡尔曼滤波的分布式驱动电动汽车车辆状态参数估计. 北京理工大学学报, 2018, 38(7): 698-702 (Wang Zhenbo, Xue Xue, Wang Yachao. State parameter estimation of distributed drive electric vehicle based on adaptive unscented Kalman filter. Transactions of Beijing Institute of Technology, 2018, 38(7): 698-702 (in Chinese) (Wang Zhenbo, Xue Xue, Wang Yachao. State parameter estimation of distributed drive electric vehicle based on adaptive unscented Kalman filter. Transactions of Beijing Institute of Technology, 2018, 38(7): 698-702 (in Chinese)).
[24] 谢文超, 赵延明, 方紫微等. 带可变遗忘因子递推最小二乘法的超级电容模组等效模型参数辨识方法. 电工技术学报, 2021, 36(5): 996-1005 (Xie Wenchao, Zhao Yanming, Fang Ziwei, et al. Variable forgetting factor recursive least squales based parameter identification method for the equivalent circuit model of the supercapacitor cell module. Transactions of China Electrotechnical Society, 2021, 36(5): 996-1005 (in Chinese) Xie Wenchao, Zhao Yanming, Fang Ziwei et al. Variable forgetting factor recursive least squales based parameter identification method for the equivalent circuit model of the supercapacitor cell module. Transactions of China Electrotechnical Society, 2021, 36(5): 996-1005 (in Chinese)).
[25] 韩伟, 王大志, 刘震. 基于可变遗忘因子RLS算法的谐波电流检测方法. 电工技术学报, 2013, 28(12): 70-74 (Han Wei, Wang Dazhi, Liu Zhen. A harmonic current detection method based on variable forgetting factor RLS algorithm. Transactions of China Electrotechnical Society, 2013, 28(12): 70-74 (in Chinese) Han Wei, Wang Dazhi, Liu Zhen. A harmonic current detection method based on variable forgetting factor RLS algorithm. Transactions of China Electrotechnical Society, 2013, 28(12): 70-74 (in Chinese
[26] 潘冬, 郭林, 李玉善. 高速公路隧道智能照明控制系统研究. 山东交通科技, 2021, 3: 89-92 (Pan Dong, Guo Lin, Li Yushan. Study on intelligent lighting control system of expressway tunnel. Shandong Communications Technology, 2021, 3: 89-92 (in Chinese) Pan Dong, Guo Lin, Li Yushan, Study on intelligent lighting control system of expressway tunnel. Shandong Communications Technology, 2021(3): 89-92 (in Chinese)
[27] 刘佳佳, 左兴权. 交叉口交通信号灯的模糊控制及优化研究. 系统仿真学报, 2020, 32(12): 2401-2408 (Liu Jiajia, Zuo Xingquan. Research on fuzzy control and optimization for traffic lights at single intersection. Journal of System Simulation, 2020, 32(12): 2401-2408 (in Chinese) Liu Jiajia, Zuo Xingquan. Research on fuzzy control and optimization for traffic lights at single intersection. Journal of System Simulation, 2020, 32(12): 2401-2408 (in Chinese)
[28] 苑绍志, 李静, 李幼德. 考虑路面不平的牵引力控制系统. 吉林大学学报(工学版), 2007, 37(5): 990-993 (Yuan Shaozhi, Li Jing, Li Youde. Traction control system with consideration of road surf ace roughness. Journal of Jilin University (Engineering and Technology Edition) , 2007, 37(5): 990-993 (in Chinese) (Yuan Shaozhi, Li Jing, Li Youde. Traction control system with consideration of road surf ace roughness. Journal of Jilin University (Engineering and Technology Edition), 2007, 37(5): 990-993 (in Chinese)).
[29] 孔磊, 宋健, 严忆泉等. 用于防抱制动系统的路面不平特征识别算法. 机械工程学报, 2007, 43(11): 229-234 (Kong Lei, Song Jian, Yan Yiquan, et al. Recognition algorithm of road roughness for anti-lock braking system. Chinese Journal of Mechanical Engineering, 2007, 43(11): 229-234 (in Chinese) doi: 10.3321/j.issn:0577-6686.2007.11.037 (Kong Lei, Song Jian, Yan Yiquan et al. Recognition algorithm of road roughness for anti-lock braking system. Chinese Journal of Mechanical Engineering, 2007, 43(11): 229-234 (in Chinese)). doi: 10.3321/j.issn:0577-6686.2007.11.037
[30] 王焕清, 陈明, 刘晓平. 一类非线性系统模糊自适应固定时间量化反馈控制. 自动化学报, 2021, 47(12): 2823-2830 (Wang Huanqing, Chen Ming, Liu Xiaoping. Fuzzy adaptive fixed-time quantized feedback control for a class of nonlinear systems. ACTA Automatica Sinica, 2021, 47(12): 2823-2830 (in Chinese) (Wang Huanqing, Chen Ming, Liu Xiaoping. Fuzzy adaptive fixed-time quantized feedback control for a class of nonlinear systems. ACTA Automatica Sinica, 2021, 47(12): 2823-2830 (in Chinese)).
[31] 邱灿文. 无人赛车纵向控制系统设计与研究. [硕士论文]. 广州: 广东工业大学, 2019 Qiu Canwen. Design and research of longitudinal control system of unmanned racing car. [Master Thesis]. Guangzhou: Guangzhou University of Technology, 2019 (in Chinese)
-
期刊类型引用(1)
1. 王德莉,李晨莹,吴炳增,焦一宇,裴海清,徐伟. 联合噪声驱动下耦合记忆阻尼Rayleigh振子族的节律模式过渡分析. 力学学报. 2024(08): 2381-2396 . 本站查看
其他类型引用(2)