作者简介:胡满江(1988-),男,助理研究员,博士研究生.研究方向:整车技术集成与安全应用.E-mail:manjiang_h@vip.163.com
为了能够实现对不同车辆质量的通用型估计,提出了一种基于纵向频响特性的整车质量估计方法。该方法通过智能网联汽车普遍获取的纵向加速度和车轮转速信息,基于车辆纵向动力学模型,推导出加速度与轮速信号间的幅频函数关系,并求取在不同频率下车辆质心处加速度与轮速之间的幅值比,进而通过最小二乘法拟合得到整车质量。CarSim-Simulink联合仿真结果表明:本文的整车质量估计方法精度在93.9%以上,可满足多车协同避撞系统车辆控制的要求。
A vehicle mass estimation method based on longitudinal frequency response characteristics is proposed for wide application for different types of vehicles. First the longitudinal acceleration and wheel speed information obtained from intelligent connected vehicle are used to derive the amplitude-frequency function between acceleration and wheel speed signal based the vehicle longitudinal dynamics model. Then, the amplitude ratio of the acceleration of the center of mass to the wheel speed at different frequencies is obtained. Finally the vehicle mass is calculated by least square fitting. The CarSim-Simulink co-simulation results show that the estimation accuracy of the proposed method is above 93.9%, which can meet the requirement of vehicle control for multi-vehicle coordinated collision avoidance system.
车辆状态数据是实现多车协同避撞控制逻辑和下层控制方法的基础和前提[1], 而整车质量作为控制期望输出力与车辆加速度间的惯性连接参数, 其估计精度的准确与否影响着车辆的控制响应[2]。整车质量和道路坡度对于车辆纵向运动的影响是同时存在的, 相关研究通常是实现两者的联合估计; 或者先估计坡度, 再估计质量。坡度的估计可以借助运动学方法实现, 如利用车辆纵向速度的变化率, 纵向加速度传感器测量值和道路坡度间的运动学关系, 还可以利用GPS测量行驶的距离, 结合大气压力传感器测量出的海拔变化, 得到车辆行驶道路的坡度角[3]。Eriksson[4]基于整车纵向动力学方程, 假定坡度已知, 通过测量发动机扭矩, 利用卡尔曼滤波观测器[5]实现质量估计, 如果道路坡度变化, 则该方法会出现较大误差。Vahidi等[6, 7, 8]利用最小二乘法实现质量和坡度的同时估计。McIntyre等[9]尝试利用模型预测控制方法估计质量。Winstead等[10]应用李雅普诺夫方法, 设计了卡尔曼滤波观测器与模型预测控制器联合估计的自适应滤波器, 理论上保证了观测器的稳定性。同济大学冯源等[11]首先利用轴荷分配规律排除道路坡度的影响, 进而结合轮毂电机驱动转矩精确可知的特点, 以分阶段式的最小二乘法实现质量估计。北京航空航天大学王健等[12]在线性二自由度单轨车辆模型的基础上, 建立了滑膜观测器估计车辆的质心侧偏角。清华大学Chu等[13]采用高频信息提取的方法, 实现了质量与坡度间的解耦估计。上述辨识控制对象关键参数的观测方法均需要车辆处于较大激励, 而智能网联汽车通常会在智能底盘集成控制[14]下处于稳态区域, 仅依靠传统时域信息的观测方法难以在近稳态激励下剥离其他参数对待估计参数的影响。
为此, 本文提出了基于频响特性的车辆关键参数辨识方法来实现整车质量估计。利用智能车辆上广泛装配的加速度传感器和轮速传感器采集纵向加速度信号和车轮转速信号; 基于整车纵向动力学方程, 推导出加速度与轮速信号间的幅频函数, 求取在不同频率下车辆质心处加速度与轮速之间的幅值比, 进而通过最小二乘法拟合得到整车质量最佳估计值表达式; 最后基于CarSim-Simulink联合仿真验证了该方法的准确性。
车辆纵向加速度是由驱动力、空气阻力、坡道阻力、滚动阻力等内外力综合作用引起的, 如图1所示, 其中
车辆的纵向动力学模型可简化为如式(1)所示。
式中:
将式(1)左、右两端同时对车辆纵向加速度
从较长距离看坡度一般为定值(除特殊情况外), 且与车辆纵向加速度不直接相关, 而坡度的瞬时变化是随机的, 可以看作为噪声项, 故道路坡度对纵向加速度的微分为零:
动力传动系统的驱动力矩一般有较高的变化频率, 相比之下, 纵向加速度的积分值— — 纵向车速的变化率要远小于加速度变化率。与纵向加速度的微分信号相比, 可以认为纵向速度的微分信号是很小的量, 即:
综合式(3)(4), 式(3)可以简化为式(5):
由式(5)可知, 当驱动力的变化率较大时, 空气阻力、坡道阻力和滚动阻力对质量估计的影响可忽略不计, 即完成了道路纵向坡度与整车质量两者估计过程的解耦。
在车辆智能安全系统中普遍运用的车辆状态信号主要包含纵向加速度、轮速等信号。若整车质量可由这两种信号估计得出, 这种方法将获得在不同智能车辆架构下的普遍适用性。为此, 本部分将结合式(5)进一步推导纵向加速度与轮速的频域关系。
式(5)右侧的分母和分子同时对时间进行微分, 可得:
进一步整理可得轮胎瞬态纵向力微分式为:
将式(7)左、右两侧同时乘以时间延迟系数
式(8)左、右两端同时对车辆纵向加速度
结合式(3)(4), 式(9)可以简化为式(10):
式中:
式(10)左、右两侧的分母和分子同时对时间进行微分, 可得:
对式(11)进一步整理可得:
结合驱动工况下, 轮胎滑移率的定义:
式中:
可将式(12)写成式(14):
由于智能群体中单车通过智能驾驶系统将自身控制在稳定区域, 即轮胎侧偏角或者滑移率很小。那么可假设
对式(15)两侧进行拉布拉斯变化, 可得:
对(16)整理, 可得两个状态量间的传递函数为:
式(17)为车辆质心处加速度与车轮轮速之间的频率关系表达式。
通过式(17)可知, 车辆质心处加速度与轮速之间的频率关系是与整车质量有关的。假设车速为已知信号, 轮胎半径为已知参数, 式(17)右侧仅整车质量为未知。为此, 首先利用周期图法(Welch's averaged periodogram method)[15]求取在不同频率下, 车辆质心处加速度与轮速之间的幅值比, 进而通过最小二乘法拟合得到整车质量。
(1)幅值比的求取方法采用下式:
式中:
通过式(18)可以得到在不同频率下的幅值比, 设
(2)以式(17)为基础, 令
$\frac{A_{\omega}(f)}{A_{a_{x}}(f)}$=$\frac{m_{v_{x}}}{K_{s}R}\sqrt{(2\pi \tau f)^{2}}+1$(19)
(3)令 $x_{i}$=$\frac{v_{x}}{K_{s}R}\sqrt{(2\pi \tau f_{i})^{2}+1}$, $y_{i}=A_{i}$。采用最小二乘法估计参数时, 要求观测值
式中:
式(20)对
求解式(21), 即可得到整车质量的最佳估计值如下:
利用CarSim-Simulink联合仿真平台验证基于整车纵向频响特性的整车质量估计方法, 仿真车辆参数设置如表1所示。
验证轮速与纵向加速度之间的传递函数, 即式(19)。设置质量
从图2中可以看出, 轮速的稳态值为27.8 rad/s, 扫频幅值为2.78 rad/s。扫频幅值相比稳态值较小且频率较高, 可以推断车辆应近似匀速行驶。受此信号驱动下, 得到的纵向速度信号如图3所示。
从图3可知, 在20 s之前, 由于轮速扫频的频率较低, 车速能够跟随轮速进行微小变化。随着频率的增大, 车速趋于稳定, 收敛于9.5 m/s。这样就符合了质量估计公式推导过程中的假设条件之一(见式(4)), 即车速变化速率远小于轮速的变化率。因此仅关注20 s以后的加速度信号, 如图4所示。
从图4可以看出, 加速度小于0.5 m/s2, 符合质量估计公式推导过程中的另一假设条件(见式(14)(15)), 即轮胎处于小滑移率范围。图4中, 虚线和划线是结合扫频信号频率与时间的关系, 由式(19)计算得到的加速度幅值-时间曲线。对比该曲线和仿真加速度信号的峰值点可以发现, 两者较好地重合在一起, 表明式(19)可以精确地描述在不同频率下的轮速与加速度间的幅值关系。
为验证基于纵向频响特性整车质量估计方法的准确性和适应性, 考虑到与商用车及其他车辆相比, 乘用车具有较高的智能化装配率, 为此, 对市场上在售的不同类型乘用车[16]整车装备质量进行调研, 将每种车辆类型中销量前10名的车型整车装备质量进行统计, 根据调研结果, 质量分布如图5所示, 统计结果见表2。
在质量估计方法理论核心, 即式(19)成立的基础上, 基于周期图的方法, 实现质量估计。在正常匀速行驶时, 轮速可以看作以随机信号。在仿真平台中构造随机的轮速信号, 如图6所示。根据上述整车质量分布与统计情况表, 分别设置仿真整车质量为1000、1500、2000和2500 kg。在图6信号作用下, 获得在不同整车质量设置下的加速度信号, 如图7所示。
对比图7(a)~(d), 可以看出质量增加导致了加速度在同等随机输入下振荡幅值减小。基于加速度信号和轮速信号数据, 采用周期图法得到不同频率下的轮速与加速度的幅值比,
如图8所示。由图8可知, 不同质量下的轮速与加速度的幅值比原始数据点均与频率间存在近似线性关系, 这与式(19)是相吻合的。 基于以上不同频率下的幅值比信息, 采用最小二乘法, 可得到整车质量估计值如表3所示。
将利用式(22)得到的数值解带入到式(17), 并以此画出轮速与加速度间的幅值比-频率曲线, 与利用功率图法得到的参考值进行趋势对比, 如图8所示。可以看出:原始数据点均匀分布于拟合曲线周围。将质量估计值与对应参考值进行对比可得, 估计精度在93.9%以上。由此可见本文的质量估计方法精度较高且具有广泛的适应性。
(1)利用纵向加速度信号和车轮转速信号, 基于整车纵向动力学方程, 推导出加速度与轮速信号间的幅频函数, 该函数与整车质量存在正相关特性。通过求取在不同频率下辆质心处加速度与轮速之间的幅值比, 采用最小二乘法拟合求得整车质量的最佳估计值表达式。
(2)采用CarSim-Simulink联合仿真验证了不同频率下的轮速与加速度间的幅值关系, 在此基础上, 根据实际整车装备质量统计结果, 对比整车质量估计值和理论值, 可得估计精度在93.9%以上, 表明本文的整车质量估计方法具有较高的准确性和适应性。
(3)本文着重从理论上探讨了该方法的可行性, 并利用仿真对方法进行了准确性验证。该方法在工程实践中的应用, 如目标车辆车轮转速信息的获取, 将在车辆智能化和网联化技术的工程应用中, 进行进一步的实践。
The authors have declared that no competing interests exist.
[1] |
|
[2] |
|
[3] |
|
[4] |
|
[5] |
|
[6] |
|
[7] |
|
[8] |
|
[9] |
|
[10] |
|
[11] |
|
[12] |
|
[13] |
|
[14] |
|
[15] |
|
[16] |
|