重载车辆轮胎模型参数辨识与灵敏度分析

黄通 刘志浩 高钦和 王冬 马栋

摘要：为了探究重载轮胎性能对多轴重型车辆行驶特性的影响，基于轮胎六分力测试试验，对GL073A型重载子午轮胎的力学模型参数进行辨识研究.针对魔术公式轮胎模型参数多，重载子午轮胎垂向载荷范围大所导致的强非线性变化的特点，提出一种基于粒子群算法和Levenberg-Marquardt算法的混合优化参数辨识方法，以轮胎六分力测试试验数据为基础，对该型重载子午轮胎的纵滑和侧偏工况的轮胎模型进行参数辨识和结果验证.结果表明，基于混合优化算法能够提高轮胎模型的参数辨识精度，辨识结果残差控制在5%以内.基于Sobol灵敏度分析方法研究了多特征参数对魔术公式轮胎模型的影响程度，以各特征参数的一阶灵敏度和总阶灵敏度作为评价标准筛选影响轮胎力学性能的主导参数.结果表明，基于Sobol灵敏度分析，从魔术公式轮胎模型的58个特征参数中选择13个主导参数，采用Sobol灵敏度分析所得出的主导参数进行辨识的结果与直接采用混合优化进行辨识的结果相比，残差最大增幅为0.138%，模型收敛速度最大增幅为30.4%.

关键词：轮胎模型；参数辨识；六分力试验；混合优化算法；灵敏度分析

中图分类号：U462文献标志码：A

Parameter Identification and Sensitivity Analysis of Heavy-duty Vehicle Tire Model

HUANG Tong，LIU Zhihao，GAO Qinhe，WANG Dong，MA Dong

（College of Missile Engineering，Rocket Force University of Engineering，Xian 710025，China）

Abstract：In order to explore the influence of heavy-duty tire performance on multi-axle heavy vehicle driving characteristics，the model parameters of the GL073A heavy-duty meridian tire were identified based on a six- component tire test. Aiming at the characteristics of strong nonlinear changes caused by multiple parameters of MF （Magic Formula）model and large vertical load range of heavy radial tire，a hybrid optimization parameter identification method based on Particle Swarm Optimization （PSO）and Levenberg-Marquardt algorithm was proposed. The parameter identification and result verification of the tire model under longitudinal slip and sidesway conditions of the heavy-load radial tire were carried out. The results show that the parameter identification accuracy of the tire model can be improved based on the hybrid optimization algorithm，and the residual of the identification results can be controlled in a range of 5%. Based on the Sobol sensitivity analysis method，the influence of multi-characteristic parameters on the tire MF model was studied. The first-order sensitivity and total order sensitivity of each characteristic parameter were used as evaluation criteria to screen out the leading parameters affecting tire mechanical properties. The results show that based on Sobol sensitivity analysis，13 leading parameters are selected from 58 characteristic parameters of the magic formula tire model. Compared with the results of direct hybrid optimization，the maximum increase of residual error is 0.138%，and the maximum increase of model convergence rate is 30.4%.

Key words：tire model；parameter identification；six component force test；hybrid optimization algorithm；sensitivity analysis

重載子午轮胎作为重载车辆与地面直接接触的部件，其力学性能直接影响车辆动力学和行驶平顺特性.与常见轿车轮胎相比，重载子午轮胎具有气压高、阻尼低、花纹粗大、扁平率大的特点.研究分析重载子午轮胎的力学性能对设计和分析先进的底盘系统和优化先进车辆系统结构至关重要[1].为了研究轮胎力学性能，国内外研究者提出了多种轮胎模型[2-3].

魔术公式（Magic Formula，MF）轮胎模型以试验数据为基础，通过精确拟合试验数据，描述轮胎的力学性能.基于MF模型，相关研究者已发展了许多适应于各种工况的不同改进形MF模型[4-5]. PAC2002模型是Pacejka[6]不断更新和完善的一种广泛应用于车辆动力学仿真和分析的MF模型，能够研究不同载荷、倾角和胎压下的轮胎力学性能，具有较强的工程应用背景.但特征参数繁多，而且高度非线性，给特征参数的辨识带来很大困难.

目前，轮胎模型的特征参数辨识基本采用数值优化算法和智能搜索算法.遗传算法因为具有较强的鲁棒特性，最先开始被应用于轮胎模型的特征参数辨识[7-8].遗传算法虽然能够在全局范围内逼近最优解，但存在局部搜索能力较差，收敛速度慢的问题.为此，相关研究者将全局范围逼近较强的智能搜索算法和局部搜索较强的数值优化算法结合起来，产生了多种轮胎参数辨识的混合优化算法，这些混合算法普遍提升了轮胎模型的辨识精度[9-11].

基于此，本文选择适合在动态和多目标优化环境中寻优的粒子群算法和Levenberg-Marquardt（LM）算法进行纯纵滑和纯侧偏工况下的轮胎模型的特征参数辨识，并对辨识后的轮胎模型进行残差分析和结果验证，最后基于Sobol灵敏度分析对该型轮胎进行参数分析和优化，研究结果可以为该型轮胎的应用提供理论支撑.

1轮胎六分力测试试验

根据轮胎动力学原理，结合现有的试验设备，进行不同垂向载荷作用下纯纵滑和纯侧偏工况的轮胎六分力测试试验.试验轮胎型号为GL073A，试验平台为六自由度轮胎力学特性试验台，试验路面为4 m水泥路面滑台，滑台运行速度为200 mm/s.纯纵滑工况试验中滑移率为-0.018～0；在纯侧偏工况试验中，侧偏角为-10°～13°.轮胎六分力测试试验如图1所示.

在轮胎压力为810 kPa下，当垂向载荷分别为24 800 N、34 700 N和44 600 N时，分别测试轮胎力学性能.测试試验结果分别如图2～图4所示.

根据轮胎力学，侧偏刚度是侧偏角度为0时的侧偏力曲线斜率；回正刚度是侧偏角度为0时的回正力矩曲线斜率；纵滑刚度是滑移率为0时的纵向力曲线斜率.由图2～图4可知：

1）当侧偏角不超过6°时，侧偏力与侧偏角呈线性关系，侧偏刚度从3 828 N/（°）增加到5 320 N/（°）；

2）侧偏角为6°时，回正力矩达到最大值，回正刚度从151 N·m/（°）增加到287 N·m/（°）；

3）纵滑刚度随垂向载荷增大，从301 799增加到375 493.

2魔术公式轮胎模型

魔术公式轮胎模型是基于试验数据通过一个公式即可表示轮胎的纵向力、侧偏力和回正力矩，公式中的特征参数物理意义明确.

基本的魔术公式为纯侧偏和纯纵滑工况下的侧偏力和纵向力：

y（x）=Dsin{C arctan[Bx-E（Bx-arctan（Bx））]}（1）

式中：y为纯纵滑工况的纵向力或纯侧偏工况的侧偏力；x为对应的滑移率或侧偏角；B为刚度因子；C为形状因子；D为峰值因子；E为曲率因子.

对基本魔术公式进行修正，可以得到纯侧偏工况下的侧偏力计算公式：

式中：γ为外倾角；F_wz为轮胎垂向载荷；F_wz0为轮胎名义垂向载荷；df_wz为名义垂向载荷增量；C_y、D_y、E_y，K_y、B_y为侧偏工况侧偏力计算因子；S_Hy、S_Vy分别为侧偏工况横向和垂向漂移量；α为侧偏角；α_y为修正侧偏角；F_y0为侧偏力；p_y为侧偏力特征参数.

df_wz=（F_wz-F_wz0）/F_wz0（4）

在纯侧偏工况中，侧偏力计算公式中待辨识特征参数共18个[12]，如表1所示.

纯侧偏工况下的回正力矩可以用轮胎拖距和残余力矩来计算：

M_z0=-tF_y0+M_zr（5）

式中：t为轮胎拖距；M_zr为残余力矩.均可采用余弦函数表达.

式中：C_t、D_t、E_t、K_t、B_t为侧偏工况回正力矩计算因子；S_Ht、S_Hr分别为侧偏工况拖距和残余力矩的横向漂移量；α_t为拖距修正侧偏角；α_r为残余力矩修正侧偏角；R₀为轮胎半径；q_z为特征参数.

式（6）中各因子均能由类似式（3）的计算公式表达出来，在纯侧偏工况中，回正力矩计算公式中待辨识的特征参数共25个[12]，如表2所示.

对基本魔术公式进行修正，可以得到纯纵滑工况下的纵向力计算公式：

式中：C_x、D_x、E_x、K_x、B_x为侧偏工况回正力矩计算因子；S_Hx、S_Vx分别为纵滑工况横向和垂向漂移量；λ为滑移率；λ_s为修正滑移率；F_x0为纵向力；p_x为特征参数.

在纯纵滑工况中，纵向力计算公式中待辨识的特征参数共15个[12]，如表3所示.

3特征参数辨识

3.1粒子群优化算法

粒子群优化算法的提出来源于鸟类群体行为的规律性，通过模拟鸟类的觅食行为，将寻求优解问题的搜索空间比作鸟类的飞行空间，每只鸟抽象成一个没有质量和体积的粒子，用该粒子表征问题的一个可能解，从搜索空间中的随机解出发，通过迭代运算寻找最优解.粒子群优化算法由于参数少且易于实现，在非线性、多峰问题中具有较强的全局搜索能力，在科学研究与工程实践中均得到广泛关注.

对于魔术公式轮胎模型的特征参数辨识问题，假设在一个D维特征参数向量解搜索空间中，有N个粒子组成一个群，其中第i个粒子在搜索空间中的位置可以表示为：

X_i=（x_i1，x_i2，…，x_iD），i=1，2，…，D（10）

第i个粒子在搜索空间中的“飞行速度”可以表示为：

V_i=（v_i1，v_i2，…，v_iD），i=1，2，…，D（11）

第i个粒子目前搜索到的最优解为：

p_best=（p_i1，p_i2，…，p_iD），i=1，2，…，D（12）

整个粒子群目前搜索到的最优解为：

g_best=（g_i1，g_i2，…，g_iD），i=1，2，…，D（13）

迭代过程中粒子更新后的粒子速度和位置分别为：

v（t+1）=w·v（t）+c₁r₁[p（t）-x（t）]+c₂r₂[g（t）-x（t）]（14）

x（t+1）=x（t）+v（t+1）（15）

式中：v（t）、x（t）分别为目前粒子速度和位置；w为惯性权重；t为迭代次数；c₁、c₂为学习因子；r₁、r₂为[0，1]内的均匀随机数.

式（14）中惯性权重w表示粒子在多大程度保留原来的速度.相关试验研究发现，较大的惯性权重能够使算法全局收敛能力强，而较小的惯性权重能够使算法局部收敛能力强.因此，本文选择由Shi等[13]提出的线性递减动态权值策略，其表达式如下：

式中：T_max为最大进化代数；w_max、w_min分别为最大惯性权重和最小惯性权重.在大多数应用中，通常取

w_max=0.9，w_min=0.4.

3.2Levenberg-Marquardt优化算法

数值优化算法具有相对较强的局部搜索能力. 本文将粒子群优化算法辨识的结果作为数值优化算法的初始值，选择收敛精度较高、非奇异性较好的Levenberg-Marquardt（LM）数值优化算法进行局部最优搜索.

LM算法通常被视为Gauss-Newton算法和最速下降法的结合，其基本原理公式为：

x（t+1）=x（t）-（J^TJ+μl）^-1g（t）（17）

式中：x为特征参数向量；g（t）为残差；μ为参数因子；J为g（t）的Jacobian矩阵.

3.3辨识结果

基于MATLAB编程实现粒子群优化算法和LM优化算法.首先，利用粒子群算法的全局搜索能力辨识出魔术公式轮胎模型特征参数的近似解，并将其作为LM算法的初始值进行局部精确搜索，得到特征参数的最终解.然后，采用式（18）对辨识结果进行残差定量分析与评价.

轮胎侧偏力、回正力矩和纵向力的辨识结果分别如图5～图7所示.

由图5～图7可以看出，通过粒子群优化算法的辨识结果与试验数据大致吻合，说明该算法能够有效获得魔术公式轮胎模型特征参数的近似解；通过混合优化算法的辨識结果与试验数据吻合度较高，辨识结果残差相对减小，说明混合优化算法能够获得更为精确的解.

受试验数据采集量多的影响，垂直载荷为44 100 N时的辨识残差相对于垂直载荷分别为24 500 N和34 300 N时的辨识残差要低，辨识结果与试验数据的吻合度相对较高.由图7可知，在纵向力辨识结果中，由于试验数据的类线性特征，混合优化算法的辨识结果对粒子群算法的辨识结果提升相对较弱.

4多特征参数敏感性分析

轮胎模型特征参数灵敏度能够反映特征参数变化对轮胎模型响应的影响程度，灵敏度较大的特征参数对轮胎模型拥有更强的调整能力，通常可以当作主导特征参数以反映轮胎模型的变化情况，特别是魔术公式这类特征参数较多的轮胎模型，确定主导特征参数能够有效提高模型调整的便捷性.

Sobol灵敏度分析方法是基于多重积分的分解方法，可以对非线性的复杂模型进行分析，研究各参数同时变化下各参数之间的耦合作用.本文选择Sobol灵敏度分析方法对魔术公式轮胎模型特征参数的影响进行分析，并以各特征参数的一阶灵敏度和总阶灵敏度作为评价标准.

4.1侧偏力特征参数

在采样数为5 000时，侧偏力特征参数的一阶和总阶灵敏度分析结果如图8所示

由图8可知，垂直漂移因子特征参数p_vy1的一阶灵敏度最大，约为0.35，总阶灵敏度稍大，约为0.18，说明p_vy1不仅对模型响应影响较大，而且与其他参数存在耦合关系，这是因为在式（2）和式（3）中，p_vy1作为一个叠加项独立存在，对模型影响比较直接.重点分析剩余一阶灵敏度较大的刚度因子特征参数p_ky1和p_ky2和水平漂移因子特征参数p_hy2.

在p_ky1和p_ky2以及p_hy2不同值下，各主要特征参数总阶灵敏度如图9所示

由图9可知：

1）峰值因子特征参数p_dy1受p_ky1、p_ky2和p_hy2的影响较大，这与图8所示的p_dy1一阶灵敏度较小，总阶灵敏度较大的结果是一致的.其灵敏度随刚度因子特征参数p_ky1和p_ky2的增大而减小，随水平漂移因子特征参数p_hy2的增大而增大.

2）垂直漂移因子特征参数p_vy1灵敏度随刚度因子特征参数p_ky1和p_ky2的增大而增大，随水平漂移因子特征参数p_hy2的增大而减小.

3）刚度因子特征参数p_ky1和p_ky2以及水平漂移因子特征参数p_hy2之间的耦合灵敏度较高，p_ky1和p_ky2相互耦合灵敏度呈正比关系.当p_hy2为14时，p_ky1和p_ky2的灵敏度均达到最高.

在侧偏力18个特征参数中，选择刚度因子特征参数p_ky1和p_ky2、水平漂移因子特征参数p_hy2和垂直漂移因子特征参数p_vy1为主导特征参数.通过调整这4 个参数来反映魔术公式轮胎模型侧偏力的变化规律，并选择峰值因子特征参数p_dy1为主要特征参数，可以通过调整峰值因子特征参数，使主导特征参数更为准确.

4.2回正力矩特征参数

在采样数为5 000时，回正力矩特征参数的一阶和总阶灵敏度分析结果如图10所示.

由图10可知，一阶灵敏度大于0.10的特征参数有q_ez1、q_ez4、q_bz1和q_hz1其中，曲率因子特征参数q_ez1的一阶灵敏度最大，约为0.38，但总阶灵敏度相对较小，说明q_ez1对模型影响较大，但与其他参数耦合相对较小；曲率因子特征参数q_ez4和水平漂移因子特征参数q_hz1也存在相同的规律.而峰值因子特征参数q_dz1和形状因子特征参数q_cz1虽然一阶灵敏度较小，但总阶灵敏度相对较大，q_dz1的总阶灵敏度最大，约为3.75.因此，重点分析总阶灵敏度较大的q_dz1和q_cz1，如图11（a）和图11（b）所示.

由图11可知，峰值因子特征参数q_dz1主要与形状因子特征参数q_cz1耦合相关，二者相互耦合灵敏度较高，随着q_dz1的增加，q_cz1灵敏度增加幅度较小，基本控制在0.4～0.45，其他特征参数灵敏度也相对稳定. 随着q_cz1的增加，q_dz1灵敏度增长幅度较大，当q_cz1为- 0.5时，q_dz1灵敏度达到最大，约为0.95；同时随着q_cz1的增大，其他特征参数的灵敏度均逐渐减小.

在回正力矩的25个特征参数中，选择曲率因子特征参数q_ez1和q_ez4、刚度因子特征参数q_bz1、水平漂移因子特征参数q_hz1以及峰值因子特征参数q_dz1为主导特征参数；选择形状因子特征参数q_cz1为主要特征参数.

4.3纵向力特征参数

在采样数为5 000时，縱向力特征参数的一阶和总阶灵敏度分析结果如图12所示.

由图12可知，刚度因子特征参数p_kx1的一阶灵敏度和总阶灵敏度最大，分别约为0.37和0.25，说明p_kx1对模型响应和其他参数的耦合关系均有较大影响.剩余一阶灵敏度较大的有：水平漂移因子特征参数p_hx1、垂直漂移因子特征参数p_vx1和峰值因子特征参数p_dx1.其中p_vx1为模型叠加项且总阶灵敏度较小，p_dx1和p_hx1的总阶灵敏度都相对较大.因此，重点分析峰值因子特征参数p_dx1和水平漂移因子特征参数p_hx1变化下的总阶灵敏度分布规律.

p_dx1和p_hx1变化下的总阶灵敏度如图13所示.由图13可知，随着p_dx1的增大，刚度因子特征参数p_kx1和水平漂移因子特征参数p_hx1均表现出较高的灵敏度，并随之增加；其他特征参数灵敏度随着p_dx1的增大而逐渐减小.随着水平漂移因子p_hx1的逐渐增大，p_dx1和p_kx1灵敏度虽然占据较大分量，但灵敏度总体呈减小趋势.

因此，在纵向力15个特征参数中，选择刚度因子特征参数p_kx1、水平漂移因子特征参数p_hx1、垂直漂移因子特征参数p_vx1以及峰值因子特征参数p_dx1为主导特征参数.

将通过Sobol灵敏度分析选择出的主导特征参数代入轮胎模型进行计算，结果如表4所示.

由表4可知，采用Sobol灵敏度分析所得的主导参数进行辨识的结果与直接采用混合优化进行辨识的结果相比，辨识残差均相对增大，但增大的幅值较小.其中垂向载荷为34 300 kN的回正力矩残差增幅相对较大，残差增幅为0.138%，这是由于主导参数辨识省略了其他特征参数影响，使残差增大，但仍小于单独粒子群算法的辨识残差.而稳定迭代次数相对于直接采用混合优化进行辨识的结果减小明显，模型收敛速度增加，最大增幅为30.4%，对轮胎模型拥有较强的调整能力，采用Sobol灵敏度分析所得出的主导参数能够有效提高模型调整的便捷性.

5结论

本文以某型轮胎六分力测试试验为基础，对该型重载子午轮胎的魔术公式轮胎模型进行参数辨识和灵敏度分析.得出以下结论：

1）基于粒子群算法能够实现对魔术公式轮胎模型特征参数的辨识，且辨识结果残差分布在5%左右.混合优化算法能够提高特征参数辨识的精确性，将辨识结果残差控制在5%以内.

2）基于Sobol灵敏度分析方法可以从魔术公式轮胎模型的58个特征参数中选择出13个特征参数作为主导特征参数，可以优先调整这13个主导参数来控制魔术公式轮胎模型的变化规律.

参考文献

[1]刘志浩，高钦和，于传强，等.重载轮胎面内欧拉梁模型及耦合振动特性研究[J].振动工程学报，2019，32（1）：107-119.

LIU Z H，GAO Q H，YU C Q，et al. In-plane Euler beam tire model and coupling vibration analysis of heavy loaded radial tire [J]. Journal of Vibration Engineering，2019，32（1）：107-119. （In Chinese）

[2]王军年，于田雨，孙娜娜，等.改善后轮独立驱动汽车过弯效率的转矩分配控制[J].湖南大学学报（自然科学版），2020，47 （12）：9-17.

WANG JN，YUTY，SUN NN，et al. Torque vectoring control of rear-wheel-independent-drive vehicle for cornering efficiency im- provement[J]. Journal of Hunan University（Natural Sciences），2020，47（12）：9-17.（In Chinese）

[3]许男，周健锋，郭孔辉，等.胎压载荷耦合效应下复合工况UniTire轮胎模型[J].机械工程学报，2020，56（16）：193-203.

XU N，ZHOU J F，GUO K H，et al. UniTire model under combined slip conditions with the coupling effect of inflation pressure and vertical load [J]. Journal of Mechanical Engineering，2020，56 （16）：193-203.（In Chinese）

[4]周衛琪，齐翔，陈龙，等.基于无迹卡尔曼滤波与遗传算法相结合的车辆状态估计[J].汽车工程，2019，41（2）：198-205.

ZHOU W Q，QI X，CHEN L，et al. Vehicle state estimation basedon the combination of unscented Kalman filtering and genetic algorithm[J]. Automotive Engineering，2019，41（2）：198-205.（In Chinese）

[5]李玲，马力，牟宇，等.模型车轮胎侧偏刚度的参数辨识方法[J].汽车工程，2016，38（12）：1508-1514.

LI L，MA L，MU Y，et al. Parameter identification method for the tire cornering stiffness of model vehicle[J].Automotive Engineer- ing，2016，38（12）：1508-1514.（In Chinese）

[6]PACEJKA H B. Tyre and vehicle dynamic [M]. 3rd ed. Oxford：Butterworth-Heinemann，2012：157-163.

[7]边伟，龚佳慧，文爱民，等.基于遗传算法的魔术公式轮胎模型参数两级辨识[J].重庆交通大学学报（自然科学版），2017，36 （5）：115-120.

BIAN W，GONG J H，WEN A M，et al. Two levels of parameter identification of Magic Formula tire model based on genetic algorithm[J].Journal of Chongqing Jiaotong University（Natural Science），2017，36（5）：115-120.（In Chinese）

[8]田晶晶，阳冬波，李枭.基于遗传算法的魔术公式轮胎模型参数辨识方法研究[J].交通节能与环保，2014，10（4）：13-16.

TIAN J J，YANG D B，LI X. Study on the parameters identification method of Magic Formula tyre model based on genetic algorithm [J]. Transport Energy Conservation & Environmental Protection，2014，10（4）：13-16.（In Chinese）

[9]张云清，陈伟，陈立平，等.Magic Formula轮胎模型参数辨识的一种混合优化方法[J].汽车工程，2007，29（3）：250-253.

ZHANG Y Q，CHEN W，CHEN L P，et al. A hybrid optimization on parameters identification for Magic Formula tire model[J]. Automotive Engineering，2007，29（3）：250-253.（In Chinese）

[10]张家旭，李静.基于混合优化方法的UniTire轮胎模型参数辨识[J].吉林大学学报（工学版），2017，47（1）：15-20.

ZHANG J X，LI J. Parameter identification for UniTire model based on hybrid optimization method[J].Journal of Jilin University（Engineering and Technology Edition），2017，47（1）：15-20. （In Chinese）

[11]田洪刚，芮强，王红岩.基于试验数据的魔术轮胎模型参数估計[J].机械工程师，2014（5）：19-22.

TIAN H G，RUI Q，WANG H Y. Identification of Magic tyre mode parameters based on experimental data[J].Mechanical Engineer，2014（5）：19-22.（In Chinese）

[12] PACEJKA H B，BESSELINK I J M. Magic formula tyre model with transient properties[J].Vehicle System Dynamics，1997，27 （1）：234-249.

[13] SHI Y，EBERHART R. A modified particle swarm optimize[C]//Proceeding IEEE International Conference on Evolutionary Computation. Anchorage：IEEE，1998：69-73.