月永冲路面下电动汽车主动悬架状态反馈H∞控制

李杰 贾长旺 成林海 赵旗

摘要：针对轮毂电机对电动汽车平顺性的负效应，建立轮毂电机电动汽车四自由度振动平面模型，研究被动悬架和主动悬架对电动汽车脉冲平顺性的影响.应用约束状态H_∞控制方法，设计轮毂电机电动汽车主动悬架控制策略，开发了MATLAB/Simulink控制仿真模型.分析无偏心、10%偏心率和20%偏心率对轮毂电机激励的影响，比较轮毂电机电动汽车被动悬架和主动悬架的脉冲路面振动响应的时间历程和4种情况的平顺性评价指标，研究结果表明，轮毂电机偏心对电动汽车产生振动激励，既影响脉冲路面平顺性，也影响状态反馈H_∞控制效果.

关键词：电动汽车；主动悬架系统；鲁棒控制；状态反馈；脉冲路面

中图分类号：U469.72文献标志码：A

State Feedback H_∞Control for Active Suspension of Electric Vehicles on Pulse Road

LI Jie，JIA Changwang，CHENG Linhai，ZHAO Qi

（State Key Laboratory of Automotive Simulation and Control，Jilin University，Changchun 130025，China）

Abstract：Aiming at the negative effect of a hub motor on the ride comfort of an electric vehicle，a four-degree of freedom vibration plane model of hub motor electric vehicle is established，in order to study the influence of passive suspension and active suspension on the pulse road ride comfort of the electric vehicle. An active suspension control strategy for hub motor electric vehicles is designed，and MATLAB/Simulink control simulation model is developed with the constrained state H_∞control method. The effects of no eccentricity，10% eccentricity，and 20% eccentricity on the excitation of the hub motor are analyzed. The time history of pulse road vibration response of passive suspension and active suspension of hub motor electric vehicle，as well as the ride comfort evaluation indexes of four cases，are compared. The results show that the eccentricity of the hub motor can produce vibration excitation on the electric vehicle，which not only affects the pulse road ride comfort but also affects the state feedback H_∞control.

Key words：electric vehicles；active suspension systems；robust control；state feedback；pulse road

與集中式电动汽车相比，轮毂电机电动汽车简化了传动结构，驱动响应快、驱动控制精确和各个车轮独立可控，是底盘优化和控制的理想载体[1-3].然而，轮毂电机也带来平顺性负效应问题[4-6].

针对轮毂电机平顺性负效应问题，已经提出了多种改善方法，涉及轮毂电机轻量化、悬架优化、动力吸振器应用、悬架控制等方面.轮毂电机轻量化是从电机设计方面考虑电机减振[7]，可以减轻非簧载质量，但是难以改变安装电机后汽车发生的变化.悬架优化是通过关键参数优化和结构改进减轻轮毂电机带来的平顺性负效应[6，8]，只能针对特定的路面和车速实现优化，无法全面适应路面和车速的各种变化.动力吸振器可以减少轮毂电机的振动[9]，但会产生在车轮内布置困难和结构复杂化的问题.悬架控制主要围绕悬架控制策略设计和执行器开发等展开，轮毂电机电动汽车可以采用PID[10]、模糊[11]、天棚[12]、地棚[13]、天棚地棚混合[13]、最优[14]和H_∞[15]等控制策略，目前主要针对单轮实现悬架控制，缺乏考虑前后车轮和空间车轮悬架控制的研究，也没有考虑电机偏心的影响.悬架控制执行器开发是通过半主动悬架[16]和主动悬架[17]实现的，半主动悬架执行器目前主要采用磁流变阻尼器[11]，通过控制阻尼力实现，需要外部能量较少；主动悬架主要包括电磁执行器[13]和液压执行器[15]，一般通过输入电能产生主动力，需要外部能量较大.

作为改善轮毂电机平顺性负效应的一种有力措施，主动悬架具有控制更好的优点.然而，主动悬架应用需要解决内部不确定性和外部干扰影响的鲁棒控制问题[18]，轮毂电机电动汽车主动悬架鲁棒控制，即H_∞控制的研究还有待深入开展.

当汽车在道路上行驶时，会遇到脉冲路面，如道路上的凸起或减速带等障碍.虽然脉冲路面的作用时间较短，但会使汽车振动突然加大，立刻降低乘员舒适性，还会对车辆零部件和运载货物造成损伤或破坏.以往的研究较少考虑脉冲路面对轮毂电机电动汽车平顺性的影响，针对轮毂电机电动汽车脉冲路面平顺性开展研究，将使轮毂电机电动汽车平顺性的研究更加全面.

本文研究轮毂电机电动汽车状态反馈H_∞控制问题，考虑脉冲路面和轮毂电机实现脉冲路面主动悬架和被动悬架的平顺性对比分析.

1轮毂电机电动汽车振动模型

1.1脉冲路面车轮激励

GB/T 4970—2009规定[19]，脉冲路面车轮激励由三角形凸块确定.脉冲路面前轮激励q_f（t）为：

式中：u为车速；h为凸块高度；l为凸块长度；t₀为汽车以车速u行驶时前轮到达凸块的时间.脉冲路面后轮激励q_r（t）为：

q_r（t）=q_f（t-t_d），t_d=L/u（2）

式中：L为车辆轴距；t_d为后轮滞后前轮的时间.

1.2轮毂电机激励

选取典型的四相8/6极开关磁阻电机作为轮毂电机，其垂向激励为单相转子垂向激励之和[13，20]，即

式中：F_v为电机垂向激励；F_vj分别为6个单相转子a、b、c、a′、b′和c′的垂向激励.

1.3车辆模型

以往研究轮毂电机电动汽车平顺性，多采用汽车二自由度振动单轮模型，具有可以揭示基本概念、基本性能和分析简单明确等优点.然而，二自由度振动单轮模型反映的是汽车一个角的作用，即单个车轮及其上面部分簧载质量的作用，只能用于研究簧载质量和车轴非簧载质量的垂直振动，无法反映簧载质量同时存在的垂直振动和俯仰振动以及两种运动对平顺性的影响，与汽车实际存在差距.而汽车四自由度振动平面模型既能反映车身质量的垂直振动和俯仰振动，也能反映前轴和后轴的非簧载质量的垂直振动，是研究轮毂电机电动汽车平顺性更合适的模型[21].

基于上述分析，建立包含轮毂电机的电动汽车四自由度振动平面模型，如图1所示.

在图1中，m_s和I_sL分别为簧载质量和簧载俯仰转动惯量；m_uf和m_ur分别为包含电机质量的前轴和后轴的非簧载质量；c_sf和c_sr分别为前轴和后轴的悬架阻尼；k_sf和k_sr分别为前轴和后轴的悬架刚度；k_tf和k_tr分别为前轴和后轴的轮胎刚度；L_f和L_r分别为簧载质量质心与前轴和后轴的距离；F_vf和F_vr分別为前轴和后轴的电机垂向激励；F_af和F_ar分别为前轴和后轴的悬架控制力；z_s和φ_s分别为簧载质量的垂向位移和俯仰角位移；z_sf和z_sr分别为前轴和后轴的悬架与簧载质量连接点垂向位移；z_uf和z_ur分别为前轴和后轴的非簧载质量垂向位移.

1.4微分方程

针对z_sf、z_sr、z_uf和z_ur，由Lagrange方程建立4个自

由度的微分方程如下：

1.5状态方程

其中

由可控性定理[18]，上述状态方程可以实现主动悬架控制.此外，当u（t）=0时，上述状态方程也适用于考虑被动悬架的轮毂电机电动汽车.

2约束状态反馈H_∞控制方法

2.1线性矩阵不等式及其求解

线性矩阵不等式F（x）<0，表示对于任意n维非零向量u，u^TF（x）u<0.F（x）的具体表示如下：

式中：x_i=[x₁，x₂，…，x_m]^T为待求的m维向量；F₀和F_i=F_i^T为已知的n×n阶对称矩阵.通过MATLAB的LMI工具箱可以求解得到x^*，以使

F（x^*）<0（11）

成立，或者无解.

2.2约束状态反馈H_∞控制

考虑如下表示：

u（t）=Kx（t）（13）

式中：K为状态反馈增益矩阵.

约束状态反馈H_∞控制问题，可以描述为：对于给定常数γ>0，求使得闭环系统稳定的状态反馈控制方程.对应式（12）第一个方程和第二个方程，有

式中：Q=P^-1；N=KQ；ρ=γ²w_max；w_max=max w（t）.

通过式（14）和式（15），约束状态反馈H_∞控制问题转化为线性矩阵不等式求解问题.在已知A、B₁、B₂、C₁、D₁、C₂、D₂和γ的条件下，其求解过程为：首先，求解式（14）和式（15）表示的线性矩阵不等式得到Q和N；其次，由Q和N确定K=NQ^-1，将u（t）=Kx（t）代入式（12）第一个方程求解，得到约束状态反馈H_∞控制的状态向量x（t）；最后，由式（12）后两个方程得到控制输出向量z_u（t）和约束输出向量z（t）.

应当说明的是，式（15）是对应于一个约束输出向量z（t）的线性矩阵不等式；当存在多个约束输出向量时，需要依次列出对应于多个约束输出向量的线性矩阵不等式.

3主动悬架状态反馈H_∞控制设计

3.1控制输出向量

3.2约束输出向量

前悬架和后悬架的动行程超过其最大值z_1max和z_2max时，悬架就会撞击限位块而产生振动和噪声.为了避免这种现象，需要

汽车行驶时，只有当车轮动载荷小于车轮静载荷时，才能始终保持车轮接地性和车辆稳定性，即

设z₂（t）=[z_sf-z_uf，z_sr-z_ur，z_uf-q_f，z_ur-q_r]^T，联立式（17）和式（18），得到对应于悬架动行程和车轮载荷的约束输出向量：

其中

为了减少控制过程的能量消耗，应当限制悬架控制力在一定范围内[22]

F_at（t）≤F_afmax，F_ar（t）≤F_armax（20）

z_u（t）=C₁x（t）+D₁u（t）（16）

式中：

式中：F_afmax、和F_armax分别为前悬架和后悬架的最大控制力.

由式（20）得到控制力约束输出向量：

3.3状态反馈H_∞控制策略

联立式（9）、式（16）、式（19）和式（21），轮毂电机电动汽车主动悬架控制问题表示为：

式（22）表示的状态反馈H_∞控制问题，可以理解为：在外界扰动向量w（t）存在最大值w_max的情况下，如果Q、N和ρ=w_maxγ²满足如下线性矩阵不等式

则主动悬架的状态反馈H_∞控制策略为：

u=Kx，K=NQ^-1（26）

4应用分析

4.1状态反馈H_∞控制仿真模型

基于上述理论研究结果开发了电动汽车脉冲路面状态反馈H_∞控制MATALB/Simulink仿真模型，如图2所示.图2中，bump模块用于生成前轮和后轮的脉冲激励；Motor Excitation模块用于生成前轮和后轮的轮毂电机激励；State Function模块用于生成状态方程、控制输出向量、约束输出向量的矩阵，以及获得状态反馈增益矩阵；Output Function用于获得每个车速下各个振动响应量的时间历程和各个振动响应量的绝对值最大值，all_data_bump模块用于以数据和图形两种形式输出所有结果.仿真时，采用某电动汽车的基本参数，如表1所示.约束条件设定为：z_1max=z_2max=0.1m，F_afmax=F_armax=1 500 N.

4.2轮毂电机激励分析

设偏心距e=（R-r）e_r，e_r是以百分率表示的偏心率，分别取e_r=0%、e_r=10%和e_r=20%，前者对应于电机无偏心情况，后者对应于2种电机偏心情况.取某开关磁阻电机参数见文献［13］，3种情况的电机偏心激励如图3所示.

由图3（a）可知，在电机无偏心的情况下，由于电机各处气隙相同，因此电机不存在垂向激励.

由图3（b）和图3（c）可知，当电机转子和定子存在偏心时，由于各处气隙不同，电机产生垂向激励.进一步分析表明，e_r取10%以上的偏心率，偏心率越大，产生的电机垂向激励越大.

由于小偏心率产生的电机垂向激励也小，因此，后续将以偏心率e_r=10%来研究电机垂向激励与脉冲路面激励共同作用，考察电机偏心在脉冲路面下对电动汽车平顺性的影响.以此能更好地说明电机偏心对汽车平顺性的影响，如果小偏心能产生明显的影响，则大偏心会产生更大的影响.

4.3振动响应时间历程

在脉冲路面、车速40 km/h和前后电机偏心率10%的情况下，求取得到两种悬架的簧载质量垂向加速度、簧载质量俯仰角加速度、前悬架动行程、后悬架动行程、前轮相对动载和后轮相对动载的时间历程，如圖4所示.

由图4可知，主动悬架与被动悬架对比，簧载质量垂向加速度和簧载质量俯仰角加速度最大值得到改善；悬架动行程和车轮相对动载最大值有所增加.

4.4平顺性评价指标的车速特性

振动响应时间历程用于分析振动响应随着时间变化的过程，而平顺性需要一个反映振动响应的整体指标.GB/T 4970—2009规定，采用振动响应时间历程的绝对值最大值作为脉冲路面平顺性评价指标，并且通过其与车速的关系评价汽车平顺性.

取最低车速10 km/h、车速增量2 km/h和最高60 km/h，在脉冲路面上仿真1 s，设定4种悬架情况：1）无偏心被动，对应于前后电机无偏心的被动悬架；2）无偏心主动，对应于前后电机无偏心的主动悬架；3）偏心被动，对应于前后电机偏心10%的被动悬架；4）偏心主动，对应于前后电机偏心10%的主动悬架.

针对1种情况，首先，求解得到1个车速对应的振动响应时间历程；其次，求解得到该车速对应的振动响应绝对值最大值；再次，求解得到各个车速对应的振动响应的绝对值最大值；最后，求解得到4种情况的全部结果，如图5所示.

由图5可知，当电机无偏心时，主动悬架与被动悬架对比结果如下：1）簧载质量垂向加速度对比改善44.11%，平均改善1.132 9 m/s²；2）簧载质量俯仰角加速度对比改善28.50%，平均改善0.622 rad/s²；3）前悬架动行程对比增加0.09%，平均增加3.06 mm，后悬架动行程对比增加0.07%，平均增加1.99 mm；4）前轮相对动载对比增加7.77%，平均增加0.133 3，后悬架相对动载对比增加5.09%，平均增加0.081 0.当电机偏心10%时，主动悬架与被动悬架对比结果如下：1）簧载质量垂向加速度对比改善42.82%，平均改善0.969 6 m/s²；2）簧载质量俯仰角加速度对比改善22.41%，平均改善0.449 9 rad/s2；3）前悬架动行程对比增加14.72%，平均增加口4.596 5 mm；后悬架动行程对比增加17.22%，平均增加4.402 5 mm；4）前輪相对动载对比增加14.72%，平均增加0.119 9；后轮相对动载对比增加14.72%，平均增加0.118 4.

5结论

1）考虑脉冲路面激励和轮毂电机激励，建立了轮毂电机电动汽车四自由度振动平面模型，既可以用于分析主动悬架的控制效果，也可以分析被动悬架的作用，适用于研究各种控制方法.

2）总结和说明了可以应用MATLAB工具箱LMI 求解的状态反馈H_∞控制问题，实现了轮毂电机电动汽车悬架状态反馈H_∞控制设计，开发了相应的MATLAB/Simulink仿真模型，通过参数设置可以研究单独脉冲路面车轮激励作用、单独电机激励作用、脉冲路面车轮和电机两种激励共同作用的效果.

3）在脉冲路面上，轮毂电机电动汽车主动悬架H_∞控制改善了簧载质量垂向加速度和簧载质量俯仰角加速度，增加了前悬架动行程、后悬架动行程、前轮相对动载和后轮相对动载，增加程度都在约束范围内；电机无偏心和偏心的主动悬架对比说明，偏心影响了主动悬架的改善能力，主动悬架设计需要考虑电机偏心情况.

参考文献

[1] KULKARNI A，RANJHA S A，KAPOOR A. A quarter-car suspension model for dynamic evaluations of an in-wheel electric vehicle [J]. Journal of Automobile Engineering，2018，232（9）：1139-1148.

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

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

[3]袁小芳，陈秋伊，黄国明，等.基于FNN的电动汽车自适应横向稳定性控制[J].湖南大学学报（自然科学版），2019，46（8）：98-104.

YUAN X F，CHEN Q Y，HUANG GM，et al. Adaptive lateral stability control of electric vehicle based on FNN[J]. Journal of Hunan University（Natural Sciences），2019，46（8）：98-104.（In Chinese）

[4] BlCEK M，KUNC R，ZUPAN S. Mechanical impact on in-wheel motors performance [J]. Journal of Mechanics，2017，33 （5）：607-618.

[5] OKSUZTEPE E. In-wheel switched reluctance motor design for electric vehicles by using a Pareto-based multiobjective differential evolution algorithm [J]. IEEE Transactions on Vehicular Technology，2017，66（6）：4706-4715.

[6]陈龙，董红亮，李利明.适合轮毂电机驱动的新型悬架系统设计[J].振动与冲击，2015，34（8）：174-180.

CHEN L，DONG H L，LI L M. A new type suspension design suitable for an in-wheel motor driving system[J]. Journal of Vibration and Shock，2015，34（8）：174-180. （In Chinese）

[7] CVETKOVSKI G，PETKOVSKA L. Genetic algorithm as a tool for multi-objective optimization of permanent magnet disc motor [J]. Archives of Electrical Engineering，2016，65（2）：285-294.

[8] YANG F X，ZHAO L L，YU Y W，et al. Analytical description of ride comfort and optimal damping of cushion-suspension for wheel-drive electric vehicles[J]. International Journal of Automotive Technology，2017，18（6）：1121-1129.

[9] LUO Y T，TAN D. Study on the dynamics of the in-wheel motor system[J]. IEEE Transactions on Vehicular Technology，2012，61 （8）：3510-3518.

[10] JIN L Q，YU Y J，FU Y. Study on the ride comfort of vehicles driven by in-wheel motors [J]. Advances in Mechanical Engineering，2016，8（3）：168781401663362.

[11]杨蔚华，李友荣，方子帆，等.多功能电动轮的轮内MR减振器设计与仿真[J].机械设计与制造，2015（6）：183-187.

YANG W H，LI Y R，FANG Z F，et al. Design and simulation of in-wheel MR damper for multifunctional motorized wheel [J]. Machinery Design & Manufacture，2015（6）：183-187. （In Chinese）

[12] KATSUYAMA E，OMAE A. Improvement of ride comfort by unsprung negative skyhook damper control using in-wheel motors[J]. SAE International Journal of Alternative Powertrains，2016，5 （1）：214-221.

[13]陳少帅.基于磁流变减振器逆模型的轮毂电机式电动汽车悬架控制研究[D].长春：吉林大学，2019：7-9.

CHEN S S. Research on suspension control of electric vehicle with in-wheel motor based on magnetorheological damper inverse model [D]. Changchun：Jilin University，2019：7-9. （In Chinese）

[14]胡一明，李以农，李哲，等.主被动一体悬架构型的多目标粒子群最优控制[J].控制理论与应用，2020，37（3）：574-583.

HU Y M，LIY N，LIZ，et al. Multi-objective particle swarm optimization linear quadratic regulator controller base on integrated suspension[J]. Control Theory & Applications，2020，37（3）：574-583.（In Chinese）

[15] SHAO X X，NAGHDYF，DU H P，et al. Output feedback H_∞control for active suspension of in-wheel motor driven electric vehicle with control faults and input delay[J]. ISA Transactions，2019，92：94-108.

[16] MA F W，WANG J W，WANG Y，et al. Optimization design ofsemi-active controller for in-wheel motors suspension[J]. Journal of Vibroengineering，2018，20（8）：2908-2924.

[17] LI Z，ZHENG L，REN Y，et al. Multi-objective optimization of active suspension system in electric vehicle with in-wheel-motor against the negative electromechanical coupling effects[J]. Mechanical Systems and Signal Processing，2019，116：545-565.

[18]姜长生，孙隆和，吴庆宪，等.系统理论与鲁棒控制[M].北京：航空工业出版社，1998：372-395.

JIANG C S，SUN L H，WU Q X，et al. System theory and robust control [M]. Beijing：Aviation Industry Press，1998：372-395. （In Chinese）

[19]汽车平顺性试验方法：GB/T4970-2009[S].北京：中国标准出版社，2009：1-10.

Method of running test-automotive ride comfort：GB/T 4970-2009 [S]. Beijing：China Standard Press，2009：1-10.（In Chinese）

[20]李杰，高雄，王培德，等.路面和开关磁阻电机作用下电动汽车振动分析[J].湖南大学学报（自然科学版），2018，45（8）：7-14.

LI J，GAO X，WANG P D，et al.Vibration analysis of electric vehicle under the action of road and switched reluctance motor[J]. Journal of Hunan University（Natural Sciences），2018，45（8）：7- 14.（In Chinese）

[21]林逸，俞凡.汽车系统动力学[M].北京：机械工业出版社，2016：196-198.

LIN Y，YU F. Vehicle system dynamics[M]. Beijing：China Machine Press，2016：196-198.（In Chinese）

[22] SHAO X X，NAGHDY F，DU H P. Reliable fuzzy H_∞control for active suspension of in-wheel motor driven electric vehicles with dynamic damping[J]. Mechanical Systems and Signal Processing，2017，87：365-383.