基于转矩协调分配的分布式驱动 电动汽车稳定性控制

张细政 郑 亮

China Mechanical Engineering - - 中国机械工程 -

引言车辆稳定性控制要求汽车行驶在一理想轨迹的同时,车辆状态始终跟随理想状态。近年来,将

收稿日期: 2017-11-17

基金项目:国家自然科学基金资助项目( 61673164,61203019);湖南省教育厅科学研究重点资助项目( 14A032);湖南省高校科技创新团队支持计划资助项目( 2014-2016) Vehicle Stability Control of Distributed-driven Electric Vehicles Based on Optimal Torque Allocation

ZHANG Xizheng ZHENG Liang

School of Electrical and Information,Hunan Institute of Engineering,Xiangtan,Hunan,411104 Abstract: A hierarchical vehicle ⁃ stability ⁃ control method was presented based on the longitudinal force distribution optimization for the handling and stability control of the distributed⁃driven electric vehi⁃ cles. The eight ⁃ degree ⁃ freedom vehicle models and the three ⁃ layer control systems were developed. By selecting the sideslip angle and the yaw rate as the state variables and introducting the virtual control to decouple two control variables,the integral 2⁃DOF vehicle models were adopted to calculate the equiva⁃ lent yaw moments for the vehicle stability in upper controllers. The linear quadratic regulator ( LQR) method was utilized to optimize the distribution of the front and rear steering angles and the tire longitudi⁃ nal forces in middle controllers. The sliding ⁃ mode ⁃ based slip controller in the lower layer was also de⁃ signed to reallocate the wheel torques. Simulation results show that the control system may make full use of the adhesion potential of the tire under high speed and extreme conditions, realize the coordinated dis⁃ tribution of wheel torques and improve the steering stability of the vehicles. When the actuators fail, the system may reconstruct effectively and realize the reallocation of control inputs to improve the safety of the vehicles.

Key words: distributed ⁃ driven electric vehicle;vehicle stability control;control allocation;wheel slip control

电机置于车轮内,车轮采用独立电机作为驱动源的新型分布式驱动电动汽车( distributed ⁃ driven electric vehicle,DEV)得到了广泛关注。在分布式驱动形式下,各轮驱动/制动转矩独立可控,只需结合汽车运行状态直接对电机施加控制,便可以精确地将转矩控制分配到每个车轮 。因此,

[] 1⁃2获得可物理实现的纵向力与转矩,并将其协调地

湖南工程学院电气信息学院,湘潭, 411104

摘要:提出一种基于车轮转矩优化分配的层次化车辆稳定性控制方法,用于分布式驱动电动汽车的操纵稳定性控制。建立八自由度车辆模型,分三层设计控制系统,上层控制器以质心侧偏角和横摆角速度为状态变量,采用积分二自由度控制模型,引入虚拟控制解耦两控制变量,计算车辆稳定的等效横摆力矩;中层采用线性二次型方法,优化分配前后轮转向角和轮胎纵向力;下层控制器设计滑模滑移率控制器,完成定滑移率下的车轮转矩再分配。仿真结果表明,该控制系统在高速极限工况下能充分利用轮胎的附着潜力,实现车轮转矩的协调分配,提高车辆的操纵稳定性;当执行机构出现故障时,系统能有效重构并实现控制量再分配,提高车辆的安全性。

关键词:分布式驱动电动汽车;稳定性控制;控制分配;滑移率控制

中图分类号: U461.6

DOI:10.3969/j.issn.1004⁃132X.2018.15.003 开放科学(资源服务)标识码(OSID) :

分配到4个车轮上,是实现DEV操纵稳定的关键性问题。

从已有研究文献来看, DEV稳定性控制多采用层次化的结构,控制上层为横摆力矩控制器设计,下层为横摆力矩的控制分配 。横摆力矩控

[] 3⁃9制器完成横摆力矩的决策和计算,根据车辆当前状态与参考值的误差计算出所需的期望横摆力矩。控制策略开始采用一些现代控制方法,如滑模控制 、模糊与神经网络控制 、模型预测

[] 3⁃4 [] 5⁃6 [] 7和最优化方法 等。目前,大部分DEV的稳定性

[] 8控制研究集中在横摆力矩控制器设计,研究比较成熟,基本都能实现优良的控制效果。这些方案在计算期望横摆力矩时,使用的控制变量是横摆角速度和质心侧偏角,而由于横摆运动和侧偏运动间的耦合关系,控制两变量所需的横摆力矩往往相反,单独依靠直接横摆力矩控制( direct yaw⁃ moment,control,DYC)很难实现两者的有效解耦,在极限工况下两者的冲突更加突出,车辆的稳定性受到较大影响 。由此,需要设计简单

[] 9⁃10可靠的横摆与侧偏运动解耦方法以提高极限工况下的车辆稳定性。

期望横摆力矩由下层控制器精确地分配到各执行器上,形成车轮的纵向力和转向角。文献[ 11 ]基于前轮转角的前馈控制和车辆状态误差的反馈控制,仅采用简单的转矩平均分配法,稳定效果不佳。文献[ 12 ]以轮胎利用率为指标进行力的优化分配,但要求所有车轮可以独立制动、独立转向。文献[ 13 ]提出将纵向力分配问题等价为带有不等式约束的非线性规划问题,但所设计的逐点二次规划算法过于复杂。文献[ 14 ]在解决了纵向合力与横摆转矩可行域的实时估计基础上,通过合力计算与优化分配,构造了内环、外环结构的车辆稳定控制系统,但忽视了车辆的转向控制能力。文献[ 15 ]采用由计算总横摆力矩的运动控制器和优化各车轮驱动转矩的控制分配器组成的分层模块化控制框架,增强了汽车的侧向横摆和侧翻稳定性能,但该方法需要进行多次在线优化,运算复杂。分配的实质是一个有约束优化问题的求解,如何充分利用协调优化的技术优势,综合考虑执行器约束和车轮滑移率等当前工作状况,设计开发更为有效的稳定性控制系统是当前DEV动力学控制研究的重要内容。

本文在八自由度车辆模型基础上构造了层次化的DEV稳定性控制系统。

1 车辆稳定控制模型

1.1 车体仿真模型

本文的车体模型假设车轮始终接触路面,侧 倾轴线在汽车中心平面上,忽略汽车的俯仰及垂直运动,共有纵向、侧向、横摆和侧倾4个自由度,再加上4个车轮自身的转动,采用八自由度非线性车辆模型进行稳定性分析与设计 ,见图1。

[] 16

FyrL FyfL

( b)悬架切面

图1 八自由度车辆模型

Fig.1 Eight degrees of freedom vehicle model车身纵向动力学方程为

ΣFx = (- v̇ m γvy )=

x

( FxfL + FxfR ) cosδf - ( FyfL + FyfR ) sinδf +

( FxrL + FxrR ) cosδr - ( FyrL + FyrR ) sinδr车身侧向动力学方程为

ΣFy = (+ v̇ γvx ) = m

y

(+ FyfL FyfR ) cosδf + ( FxfL + FxfR ) sinδf +

( FyrL + FyrR ) cosδr + ( FxrL + FxrR ) sinδr车身横摆动力学方程为

ΣMz == γ̇

I z a [( FyfL + FyfR ) cosδf - ( FxfL +

FxfR ) sinδf ] - 0.5tf [ ( FyfL - FyfR ) sinδf -

( FxfL - FxfR ) cosδf ] -b [( FyrL + FyrR ) cosδr +

( FxrL + FxrR ) sinδr ] - 0.5t r [ ( FyrL - FyrR ) sinδr -

( FxrL - FxrR ) cosδr ] ( 3)车身侧倾动力学方程为

ΣMϕ = =-

I x ϕ̈ ms hs cosϕ (+ v̇ γvx ) -

y

I xz γ̇ - Bϕ ϕ̇ - Kϕϕ车轮运动学方程为

Jω ωij = Td, - Tb, - Fx, Rw

ij ij ij

由式( 2)得车身侧偏动力学方程: β̇ = - v̇ tanβ/ v - γ + ( v̇ + γvx ) /( vcosβ )

y

式中,、v vx、vy分别为质心、纵向和横向车速; δf、δr分别为前轮、后轮转向角; γ为横摆角速度; β为质心侧偏角; Iz、Ix分别为横摆转动惯量和侧倾转动惯量; ϕ 为车身侧倾角; FxfL、FxfR、FxrL、FxrR、FyfL、FyfR、FyrL、FyrR分别为沿纵向、侧向的前左、前右、后左和后右轮胎力分量;下标fr、分别表示前后轮, LR、 分别表示左右轮; Td,、Tb, i= f,r; j= L,R)

ij( ij分别为车轮驱动力矩和制动力矩; m为整车质量; ms为悬挂质量; ab、分别为前轴、后轴与质心间的水平距离; tf、tr分别为前轮、后轮距; Bϕ、Kϕ分别为悬架侧倾角阻尼和单侧轮胎侧偏刚度; Jω、ωij、Rw分别为车轮滚动惯量、角速度和滚动半径; hs为簧载质量的侧倾力臂。

1.2 非线性轮胎模型

采用“魔术公式”轮胎模型模拟轮胎的非线性力学特性 ,为简单起见,不考虑纵向力和轮胎回

[] 16正力矩的作用,忽略曲线的垂直和水平方向漂移,轮胎侧向力

Fy = -D sin ( Carctan ( Bα - E ( Bα - arctan ( Bα )))) ( 7) D = a1 F + a2 Fz C = a0 E = a6 Fz + a7

2 z

B = a3 ( CD )- ( 1 - a5 | γ | ) sin ( 2arctan ( Fz a- ))

1 1

4

式中, D为峰值因子; C为曲线形状因子; α为侧偏角; E为曲率因子; B为刚度因子;系数a ( i= 1,2,…,7)的值采用SAE⁃870421中的轮胎模型数据 。

[] 17轮胎附着力大小由近似摩擦圆限制,最大轮胎纵向力

Fx, max =

式中, μ为路面附着系数。

由式( 7)、式( 8)知,在车辆行驶过程中,轮胎受到的垂直载荷对纵向力和侧向力有着重要影响,考虑汽车加速、减速及转向时汽车的载荷转移影响,前轴、后轴垂直载荷差为

ΔFzf = t ( Bϕf ϕ̇ + Kϕf ϕ + ms hf ( v̇ sinβ + ü

-1 f

ï vβ̇cosβ + γvcosβ ))

ý

ΔFzr = t ( Bϕr ϕ̇ + Kϕr ϕ + ms hr ( v̇ sinβ + ï

-1 r vβ̇cosβ + γvcosβ ))

þ轮胎垂直载荷为

FzrL = Fzr - ΔFzr ü

FzfL = Fzf - ΔFzf

ý

FzrR = Fzr + ΔFzr

FzfR = Fzf + ΔFzf þ轮胎的侧偏角计算公式为

vsinβ + aγ ü αfL = - δf + arctan ï vcosβ + 0.5tf γ ï

ï vsinβ + aγ αfR = - δf + arctan vcosβ - 0.5tf γ ï

ý vsinβ - bγ αrL = - δr + arctan ïï vcosβ + 0.5tr γ ï

vsinβ - bγ αrR = - δr + arctan ïï

vcosβ - 0.5tr γ þ车辆稳定的分层控制器设计基于分层控制思想,分三层设计控制系统:上· · 层为横摆力矩控制器,中层为线性二次型调节器( LQR)优化分配器,下层为基于滑模变结构的滑移率控制器。

2.1 稳定性控制参考模型

为得到车辆横摆运动的期望值,本文采用二自由度车辆模型设计理想质心侧偏角βd和理想横摆角速度γd。对于质心侧偏角,应该尽量保持βd为零;而理想横摆角速度是要保证在某一车速下,驾驶员感觉其转向输入与车辆横摆角速度间成近似线性关系,以减小操纵难度。不考虑各轮轮心的纵向速度差异,参考模型的稳态响应为

[] 17 β* =0 ( 12)

d γd | γd | < μgv -

1 x γ* = ( 13) d 1 1 μgv - sgn γd | γd |≥ μgv -

x x v2 vx γd = δf v ch = gL/Kus

ch

L v2 + v2

ch x mg b a

Kus = ( - )

L Cαf Cαr

式中, vch为稳定性参数; Cα f、Cα r为前轴、后轴轮胎的侧偏刚度。

2.2 横摆力矩LQR控制器

当车辆的实际横摆角速度和质心侧偏角与理想的横摆角速度和质心侧偏角出现偏差时,上层控制器在汽车上施加一给定的附加横摆力矩来完成对车辆行驶状态的纠正,使其始终跟随参考模型所产生的理想值。考虑如下积分二自由度车辆控制模型:

é C0 é β̇ ù êê -

mv ê ú ê γ= C1

- ë û êê Iz

ë 0

( 14) C0 =+ Cαf Cαr C1 = aCαf - bCαr C2 = a2 Cαf + b2 Cαr式中, ψ为横摆角; uu ˉ 、ˉ 分别为跟踪参考模型的质心侧偏角

β γ

和横摆角速度所需的等效控制量。

采用LQR方法设计状态反馈控制律:

é ù é u ˉ ù β

=- K ė u ˉ γ

ë eγ û

∫ eγ = ė ė =- γ γd

γ γ

式中, ė 为跟踪误差; K为LQR增益。

γ

给定闭环控制系统式( 14)的期望极点,利用极点配置方法,即可求出闭环反馈增益K。

不同于其他基于二自由度车辆模型的稳定性系统 ,式( 14)中控制矩阵的结构特点表明等

[ 1,4 ,7]效控制两分量已经解耦,可以独立设计;同时,式( 15)的控制量设计引入了积分模型,能有效改进横摆角速度跟踪时的稳态误差。

i

2.3 控制分配器

计算出的等效控制量需要在轮间进行分配,即根据车辆行驶状态的识别、驾驶员转向盘的检测等,判断车辆在哪个轮上或哪几个轮上分配横摆力矩,通过控制分配器向4个轮毂电机发出驱动或制动指令,对车轮进行转矩调节,最终转化为轮胎与地面的相互作用力,以此来达到保持车辆稳定的目的。对于DEV电动汽车,车辆由4个轮毂电机进行驱动,前后轴均可独立转向,整个系统可通过多种灵活的输入组合形式来完成等效横摆力矩的轮间分配。考虑如下分配:

Bu =[ˉ u u ˉ ]

T β γ

式中, B为控制输入矩阵; u为实际控制输入量。

本文采用前后轮独立转向、转矩独立控制的DEV,则有é Cαf Cαr ê ê mv mv 0 0 0

B =

Cαf Cαr tf tf tr a - b - - ë Iz Iz 2I z 2I z 2I z u =[ δf δr FxfR FxfL FxrR FxrL ]

T

式( 16)中存在执行器的冗余。控制分配器的作用是在同时考虑地面附着条件限制和执行器位置/速率限制的情况下,完成各控制输入量的最优分配,其本质上是一个有约束的多变量优化问题。等式( 16)的控制量分配可转化为有约束的二次优化问题,即min

u s.t. uT Qu + CT u

ý

Bu =[ u ˉ u ˉ]

T β γ u lb, ≤≤ ui u ub, [ ui ] = u þ

i i

式中, Q为对称正定矩阵; C为线性系数矩阵; u lb, 、u ub, 为执

i i行器约束的上下界。

式( 17)的优化问题既包含了车辆稳定性控制中对状态误差的要求,又体现了对独立制动产生的汽车横摆力矩的要求。观察控制矩阵B的结构,发现控制输入[ δf δr ]直接作用并决定了质心侧偏角和横摆角速度,具有更大的控制效果;而四轮驱动力Fx 仅直接作用在横摆角速度状态方程上,经

,ij由两个状态间的耦合作用来间接对质心侧偏角施加影响。这种耦合可能使优化过程中控制目标间冲突,尤其是在执行器故障的情况下,需要优先保证横摆角速度的稳定性控制。在式( 16)中引入虚拟控制量τ,新的状态方程为

é C0 C1 ù é β̇ ù êê mv +1 úú

mv2 é β ù ê ú =C1 +[

ê Bκu ]ˉ γ̇ C2 γ ê úë ë Iz Iz v û u ˉ =[ uτ ] κ =[ κβ κγ ]

T T

其中, κ为控制矩阵。虚拟输入u ˉ中通过引入τ,增 加了控制自由度,给定实际控制输入u,在实现横摆角速度跟踪优先控制的同时,可以放宽对质心侧偏角的等式约束,从而达到两个状态间解耦目的,因而, κ值选择时,要求γ → τ的直流增益为零。由式( 18)求得输入τ对状态γ的传输函数:

C0 C1 ü sL ( β )=- Lβ ( )-( 1 + ) Lγ ( )+ κβ L ( τ ) mv mv2

ý

C1 C2 sL ( γ )=- Lβ ( )- Lγ ( )+ κγ L ( τ )

Iz Iz v þ

( 19)

式中, L表示拉普拉斯变换。

由式( 19)可求得: Lγ ( )

=

L ( τ )

κγ mIz vs + κγ Iz C0 - κβ C1 mv

( 20) mIz vs2 +( Iz C0 + mC2) s - C1 mv2 -+ C C0 C2 Iz

2

1

令传输函数式( 20)的直流增益为零,则有κγ Iz C0 -

mvC1 κβ C1 mv = 0,求得κβ = 1, κγ = 。则优化问

Iz C0

题式( 17)转化为如下模型min u ˉ TQu ˉ+ CT u ˉ

u ˉ s.t. [ Bκu ] ˉ =[ u ˉ

β u ˉ ≤ˉ u ≤ˉ u

lb, i ub, i

2.4 车轮滑模滑移率控制器

车轮纵向力与轮胎滑移率直接相关,在给定期望滑移率下,设计滑移率控制器( slip rate con⁃ troller,SRC)来调节各车轮纵向力大小。在车辆转向和横摆运动的影响下,沿车轮坐标系的轮心纵向车速和滑移率分别计算如下: v fL = ( vx - 0.5tf γ ) cosδf + ( vy + aγ ) sinδf ü vfR = ( vx + 0.5tf γ ) cosδf + ( vy + aγ ) sinδf

ý vrL = ( vx - 0.5tr γ ) cosδr + ( vy - aγ ) sinδr vrR = ( vx + 0.5tr γ ) cosδr + ( vy - aγ ) sinδr þ

Rw ωij / v ij -1 车轮制动时σ ij = ( 23)

1- v ij / ( Rw ωij ) 车轮加速时在定侧偏角下,轮胎纵向力会随着纵向滑移率的增加而先快速增大,在临界滑移率σ∗处达到峰值,滑移率再继续增加后轮胎纵向力反而减小。当实际滑移率较小时( |σ | < | σ ∗ |),控制车轮电机转矩Tij与目标纵向力Fx, 成比例,即

ij

Tij = Rw, Fx, ( 24)

ij ij然而,在一定的垂直载荷下,由于车轮附着力的限制,当目标轮胎纵向力Fx, 过大时,按式( 24)控制

ij车轮转矩会使得实际滑移率超出临界滑移率,从而造成车轮抱死和纵向力减小,因此,当实际滑移率较大时( |σ |≥ | σ ∗ |),需要采取滑移率控制策略来调节电机转矩,从而限制车轮滑移率。

制动工况下,构造以滑移率误差为变量的滑模函数为

Sσ =- σ σ∗ σ∗ <0

选取Lyapunov稳定性函数

V = 0.5S ≥ 0 ( 26)轮胎纵向力可表示为线性部分和非线性增量之和,即

Fx, = Cσ σ + ΔFx, ( 27)

∗ ij ij ij

式中, Cσ为轮胎工作在线性区域时的纵向刚度; ΔFx, 为非

ij线性增量,且满足0 ≤ ΔFx, ≤ ij ;为非线性增量最

ij x, x, ij

大值。

结合车轮运动方程式( 5)和式( 27),可得

Rw R2 w ΔFx, v̇ σ̇ = - ( Tb, + Rw Cσ σ∗ ij ) - -( 1 + σ ij )

ij x, ij ij ij

Jω vx, Jω vx, vx, ij ij ij

( 28)控制器稳定的必要条件是滑动模态到达条件需要得到满足,即

= Sσ σ̇ ≤- ησ | Sσ |

ij

将式( 28)代入式( 29),有

Rw R2 w ΔFx,

S- [ ( Tb, + Rw Cσ σ )- -

∗ ij ij ij

Jω vx, Jω vx, ij ij

(1+ σ ij ) ]≤- ησ | Sσ | ( 30)

x, ij vx,

ij当车辆制动行驶时,设计车轮转矩控制量

J ω v̇ Iω vx,

Tb, = Rw Cσ σ + (1+ σ ij )- kσ sgnSσ ( 31)

∗ x, ij ij ij

Rw Rw ij

R2 kσ > + ησ

W

x, ij

Iω vx,

ij

式中, kσ滑模增益。

由式( 30)和式( 31)知,当Sσ ≠0时, < 0,而正V定,则此时整个滑移率误差动态是一致渐近稳定的,零点是系统的一个稳定点。同上,加速工况下,设计车轮转矩控制量

I ω v̇ Iω vx,

Td, = Rw Cσ σ + - kσ sgnSσ

∗ x, ij ij ij

Rw (1- σ ij ) Rw (1- σ ij )

2 ij

( 32)设计的滑模滑移率控制器的工作原理见图2。

vx ω

图2 滑模滑移率控制器工作原理

Fig.2 Principle of sliding mode slip controller仿真结果与分析

在 MATLAB/Simulink软件环境下构造车体和轮胎模型,并对控制系统进行仿真。选取DEV 轿车参数为: a = 1.013 m, b = 1.702 m,

1784 · Tf = 1.554 m, Tr = 1.534 m, m = 1.401 kg, Ix= 550.2 kg · m2, Iy = 2 587.1 kg · m2 , Iz = 2 677.2 kg · m2, hs = 0.582 2 m, hf = hr = 0.127 m, Kϕf = 4 297.2 N · m/rad, Kϕr = 3 724.2 N · m/rad, Bϕf = 900 N · m · s /rad, Bϕr = 850 N · m · s / rad, Cαf = 113.2 kN · m/rad, Cαr = 90.2 kN · m/rad。

闭环系统期望极点设置为[ - 7 - 2 + 5i -2-5i ],求得反馈增益K= [ 0.466 7 - 0.926 8 - 0.008 3 0.042 8 3.993 7 29.000 3 ]。分层控制系统见图3。δend β*

d

-+ ven d γ* - + d

∫(

图 3 分层控制系统框图

Fig.3 Hierarchical control system block diagram分别在转向角阶跃输入、正弦输入和转向器故障3种工况下进行仿真分析。为便于比较,仿真时给出了本文设计的优化分配( O)、DYC 控制( DYC)和无控制( w/o)下的稳定性能曲线。

V ̇

3.1 阶跃转向工况

图4~图8所示为车速80 km / h、路面摩擦因数 0.8,转向盘阶跃输入(转角2π/3)时的仿真结果。仿真结果显示,无控制时车辆的横摆角速度γw/ o和质心侧偏角βw/ o很快达到危险的临界值,表明此时车辆已经接近失控,需要采取相应的控制策略才能保证车辆的期望转向。 图4 质心侧偏角响应曲线(阶跃转向工况)

Fig.4 The responses of centroid side-angle

( step steering)

从图4、图5中可以看出,在采取DYC控制和优化分配后,两种控制器的横摆角速度和质心侧偏角都能较好地跟踪理想模型,使车辆稳定转向。 - F

图5 横摆角速度响应曲线(阶跃转向工况) Fig.5 The yaw rate responses( step steering)

图6 四轮轮胎纵向力(阶跃转向工况) Fig.6 Four-wheel tire longitudinal forces

( step steering)

图7 四轮轮胎滑移率(阶跃转向工况) Fig.7 Four-wheel tire slip rates( step steering)

图8 车辆行驶轨迹(阶跃转向工况)

Fig.8 The vehicle desired and actual rajectories

( step steering)

相较于DYC控制,本文设计的优化分配方法通过在四轮施加适当的制动力矩,车辆的横摆角速度响应γo能更好地跟踪理想横摆角速度γd ,并稳定在 0.226 rad/s左右;同时车辆质心侧偏角βo也得到很好的抑制,并稳定在理想的零侧偏角范围内,表明该控制系统有很好的稳定性。图6给出了四轮轮胎的纵向力。图7为四轮滑移率绝对值,在转 向开始后,滑移率开始上升,但在滑移率控制器的作用下始终小于0.2;同时,由图可见四轮中左后轮滑移率相对较大,而右前轮的滑移率相对较小,反映了车辆左转时车速不变情况下左后轮轮速较小、右前轮轮速较大的实际情况。图8给出了本文优化分配方法下车辆质心的期望轨迹、实际轨迹和航向,可见车辆能很好地跟踪理想行驶轨迹。

3.2 正弦转向工况

图9~图13所示为车速80 km/h、路面摩擦因数0.8,转向盘正弦输入(转角±π/2,频率0.25 Hz)时的仿真结果。结果显示,无控制时车辆的横摆角速度γw/ o已不能跟踪理想单移线行驶,质心侧偏角 βw/ o则很快达到极危险的临界值,车辆此时已经失控,需要进行必要控制才能保证车辆的单移线变道行驶。 图9 质心侧偏角响应曲线(正弦转向工况) Fig.9 The centroid lateral responses

( sinusoidal steering)

图10 横摆角速度响应曲线(正弦转向工况) Fig.10 The yaw rate responses( sinusoidal steering)

图11 四轮轮胎纵向力(正弦转向工况) Fig.11 Four-wheel tire longitudinal forces

( sinusoidal steering)

图12 四轮轮胎滑移率(正弦转向工况)

Fig.12 Four-wheel tire slip rates( sinusoidal steering)

图13 车辆行驶轨迹(正弦转向工况)

Fig.13 The vehicle desired and actual rajectories

( sinusoidal steering)

从图9、图10中可以看出,在采取DYC控制和优化分配后,两种控制器的横摆角速度和质心侧偏角都能较好地跟踪理想模型。相较于DYC控制,本文优化分配下的车辆横摆角速度响应γo 能更好地跟踪理想横摆角速度γd,车辆质心侧偏角βo也得到了很好抑制,显著小于DYC控制质心侧偏角 βDYC ,车辆能够很好地稳定变道。图11给出了四轮轮胎纵向力。图12为轮胎滑移率,开始变道后,滑移率迅速上升,但在滑移率控制器的作用下始终小于0.2。图13给出了车辆质心的期望轨迹、实际轨迹和航向,表明车辆能按理想轨迹完成变道行驶。

3.3 转向故障工况

正弦工况的仿真过程中, 4s≤ to < 5s时设置前轮转向失效, 5s≤ to ≤6s时后轮发生转向失效。图14~图 17给出了转向故障工况的仿真结

图14 质心侧偏角响应曲线(转向故障工况) Fig.14 The centroid angle response( steering failure) 图15 横摆角速度响应曲线(转向故障工况) Fig.15 The yaw rate response( steering failure) 图16 前轮、后轮转向角曲线(转向故障工况) Fig.16 Front and rear wheel steering angles

( steering failure)

图17 车辆行驶轨迹(转向故障工况)

Fig.17 The vehicle desired and actual rajectories

( steering failure)

果。由图14可以看出,在转向失效的时间段,虽然质心侧偏角有较大波动但小于10°;图15中横摆角速度有一定的波动,但仍在车辆可控范围内。事实上在转向器发生故障后,控制策略需要在余下可用的执行器中重新计算并分配控制量,以适应新的控制环境。图16给出了前轮、后轮转向角,在图中的故障时段内,转向角调节失效并锁定。图17给出了车辆行驶轨迹。

3种工况下的仿真结果表明,所设计的控制系统能够很好地跟踪理想模型,保持车辆在高速转向、变道和转向故障等情况下的操纵稳定性,且系统响应迅速,具有很好的控制效果。

4 结论

本文建立了DEV车辆稳定性模型,根据各车轮转矩可以直接控制的特点,以质心侧偏角和横摆角速度为控制变量,设计了基于转矩协调的分层车辆稳定性控制方案。采用引入横摆角后的积分二自由度控制模型来计算期望横摆力矩,实现两控制变量的解耦,减小稳态跟踪误差;基于LQR

方法进行控制量的优化分配,引入虚拟控制量,在不影响横摆角速度跟踪的情况下,放宽对质心侧偏角的约束。设计滑模滑移率控制器,实现定滑移率的车轮转矩控制。对汽车在高速转向和移线行驶的不同工况进行了仿真,结果表明控制系统能很好地跟踪横摆角速度的期望值,并保持质心侧偏角最小;当转向机构出现故障时,系统能有效重构并实现控制量再分配。该系统充分发挥了DEV的优点,有效地改善了汽车行驶的稳定性,提高了车辆安全性。

参考文献:

[] 1 TERASHIMA M,ASHIKAGA T,MIZUNO T, et al. Novel Motors and Controllers for High⁃perfor⁃ mance Electric Vehicle with Four In ⁃ wheel ⁃ motors [] J . IEEE Transactions on Industrial Electronics, 1997,44(1):28⁃38.

[ 2 ] 张立国,宁国宝. 国内电动汽车发展综述[ J. ]汽车工业研究, 2006,13(12):30⁃33.

ZHANG Liguo,NING Guobao. Overview of Do⁃ mestic Electric Vehicle Development [] J . Automo⁃ tive Industry Research,2006,13(12):30⁃33.

WANG J M,LONGORIA R G. Coordinated and Reconfigurable Vehicle Dynamics Control [] J . IEEE Transactions on Control Systems Technology, 2009,17(3):723⁃732.

[ 4 ] 赵艳娥,张建武.基于滑模控制的四轮驱动电动汽车稳定性控制[ J. ] 上海交通大学学报, 2009,(10): 1526⁃1530.

ZHAO Yan e,ZHANG Jianwu. Stability Control for a Four ⁃ motor ⁃ wheel Drive Electric Vehicle Based on Sliding Mode Control [] J . Journal of Shanghai Ji⁃ aotong University,2009,(10):1526⁃1530.

[] 5 TAHAMI F,FARHANGI S,KAZEMI R. A Fuz⁃ zy Logic Direct Yaw Moment Control System for All ⁃ wheel ⁃ drive Electric Vehicles [] J . Vehicle Sys⁃ tem Dynamics,2004,41(3):203⁃221.

[ 6 ] 沈勇,吴新文. 基于复合神经网络模型的四轮独立驱动电动车控制[ J. ]汽车工程, 2004,26(4):458⁃460. SHEN Yong,WU Xinwen. Multiplex Neural Net⁃ work Control of Four⁃wheel Drived Electric Vehicle [] J . Automotive Engineering,2004,26(4):458 ⁃ 460.

[ 7 ] 杨建森,李飞,丁海涛,等.基于广义预测控制的汽车横摆稳定性控制[]农业机械学报, J. 2012,43(1):1⁃5. YANG Jiansen,LI Fei,DING Haitao,et al. Vehi⁃ cle YawStability Control Based on Generalized Pre⁃ dictive Control [] J . Chinese Journal of Agricultural Machinery,2012,43(1):1⁃5. [] 8 MOKHIAMAR O,ABE M. How the Four Wheels Should Share Forces in an Optimum Cooperative Chassis Control [] J . Control Engineering Practice, 2006,14(3):295⁃304.

KATSUYAMA E. Decoupled 3D Moment Control Using In⁃weel⁃motors [] J . Vehicle System Dynam⁃ ics,2012,51(1):18⁃31.

[] 10 NAGASE H,INOUE T,HORI Y. Decoupling Control of P and Y for High Performance AFS and DYC of 4 Wheel Motor Electric Vehicle [] C // Proceedings of the International Symposium on Ad⁃ vanced Vehicle Control ( AVEC). Hiroshima, 2002:1⁃7.

SHINO M,NAGAI M. Yaw Moment Control of Electric Vehicle for Improving Handling and Stabil⁃ ity [] J . JSAE Review,2001,22(4):473⁃480. MOKHIAMAR O,ABE M. Simultaneous Opti⁃ mal Distribution of Lateral Tire Forces for the Model Following Control [] J . Journal of Dynamic System, Measurement and Control, 2004, 126 ( 4):753⁃763.

[] 13 ONE E,HATTORI Y. Vehicle Dynamics Control Based on Tire Grip Margin [] C //The 7th Interna⁃ tional Symposium on Advanced Vehicle Control. Netherlands:Swets and Zeitlinger,2004:531⁃536. [ 14 ] 刘跃,方敏,汪洪波.车辆稳定控制中的合力计算与分配[] J. 控制理论与应用, 2013,30(9):1122⁃1130. LIU Yue,FANG Min, WANG Hongbo. Force Calculation and Distribution i n Vehicle Stability Control [] J . Control Theory and Application, 2013,30(9):1122⁃1130.

KANG J Y,YOO J H,YI K. Driving Control Al⁃ gorithm for Maneuverability,Lateral Stability and Rollover Prevention of 4WD Electric Vehicles with Independently ⁃ driven Front and Rear Wheels [] J. IEEE Transactions on Vehicular Technology, 2011,60(7):2987⁃3001.

[] 16 RAJAMANI R. 车辆动力学及控制[ M. ] 北京:机械工业出版社, 2011.

RAJAMANI R. Vehicle Dynamics and Control [] M . Beijing:China,Machine Press,2011.

[] 17 WANKI C,JANGYEOL Y,JEONGTAE K,et al. An Investigation i nto Unified Chassis Control Scheme for Optimised Vehicle Stability and Ma⁃ noeuvrability [] J Vehicle System Dynamics,2008, 46:(S1),87⁃105.

(编辑 袁兴玲)

作者简介:张细政,男, 1979年生,教授,博士。研究方向为智能控制技术及应用。获省部级二等奖2项。出版专著1部,发表论文30余篇。E⁃ mail:[email protected]

Newspapers in Chinese (Simplified)

Newspapers from China

© PressReader. All rights reserved.