《1 前言》

1 前言

四轮转向 (four-wheel steering, 4WS) 是提高车辆操纵稳定性和主动安全性的有效控制策略, 已成为目前的研究热点之一。4WS控制系统的研究大多以经典控制理论为基础[1], 依赖于一个确定性模型。考虑汽车是包含各种不确定因素的非线性、模型不确定性的复杂系统, 车辆在行驶过程中受到很多不确定性因素的影响, 诸如载荷变化、外部干扰、路面不平及刹车/加速时轮胎侧偏特性的变化等, 因而针对此类情况的鲁棒控制设计受到重视。设计车辆鲁棒控制器有可能解决使车辆控制系统性能恶化的不确定性, 常见的有H2/H控制等[1]

常规的H控制由于实质是将各类不确定性总体叠加, 导致对鲁棒性而言常常过分保守。而μ综合鲁棒控制方法是鲁棒控制分析与综合的一种有效手段。该方法对不同不确定性分别建模, 可同时兼顾系统的稳定鲁棒性和性能鲁棒性, 通过定义结构奇异值作为度量, 采用矩阵不等式方法 (LMI) 确定其上下界, 引入定标矩阵实现上界的迭代逼近[2]。此方法大大改善了H控制设计的保守性。

本文根据4WS车辆控制器设计对稳定鲁棒性和性能鲁棒性的严格要求, 提出横摆角速度跟踪控制算法框架, 根据系统性能的要求选取不同环节的权函数, 采用μ综合方法实现4WS车辆控制器设计。仿真结果表明, 4WS控制系统对模型参数的摄动具有良好的动态性能、稳定鲁棒性和性能鲁棒性, 有效地提高了车辆操纵稳定性和安全性。

《2 车辆四轮转向系统动力学模型》

2 车辆四轮转向系统动力学模型

研究四轮转向汽车操纵稳定性可以采用两自由度转向模型, 其简化数学表达式为

mv(β˚+γ)=Fy1+Fy2,Ιzγ˙=LfFy1-LrFy2(1)Fy1=Κf(δf-β-Lfγv),Fy2=Κr(δr-β+Lrγv)(2)

式中δfδr为前、后轮转角;Fy1Fy2为前、后轮所受的侧向力, 线性轮胎模型是假定汽车转向时前、后轮所受侧向力大小与前后轮侧偏角成线性关系;m为汽车质量;Iz为横摆惯性力矩; β为质心点的侧偏角;γ为横摆角速度;v为行驶速度;Lf, Lr为质心至前、后轴的距离;Kf, Kr为前、后轮的侧偏刚度。

根据汽车转向响应的横摆角速度γ和车辆质心的侧偏角β, 可得到侧向加速度

ay=d2ydt2=v(γ+dβdt)(3)

式中y为侧向位移。

由式 (1) 和式 (2) 可得到如下状态方程形式的线性汽车模型:

[β˙γ˙]=[-(Κf+Κr)mv-1-(LfΚf-LrΚr)mv2(-LfΚf+LrΚr)Ιz)-(Lf2Κf+Lr2ΚrΙzv][βγ]+[ΚfmvΚrmvLfΚfΙz-LrΚrΙz][δfδr](4)

考虑到方向盘转角与前轮转角一般都是比例关系, 仿真计算以前轮转角为输入。由于横摆角速度便于测量, 而质心侧偏角难以测量, 控制系统只取横摆角速度作为反馈变量来设计方向盘转角跟踪横摆角速度的μ综合鲁棒控制, 如图1所示4WS车辆闭环控制系统。前、后轮转向角到车辆质心横摆角速度的传递函数为

Gfγ=b1s+b2s2+a1s+a2,Grγ=b3s+b4s2+a1s+a2(5)

其中

b1=LfΚfΙz;b2=ΚfΚr(Lf+Lr)mvΙz;b3=-LrΚrΙz,b4=-b2;a1=Lf2Κf+Lr2ΚrΙzv+Κf+Κrmv;a2=ΚfΚr(Lf+Lr)2mv2Ιz+LrΚr-LfΚfΙz

理想横摆角速度增益定义为γfγ(s)=1+ks1+τsb2a2kτ为理想增益系数, 经反复试算得到。

《图1》

图1 4WS车辆闭环控制系统

图1 4WS车辆闭环控制系统  

Fig.1 4WS vehicle closed-loop control system

《3 μ综合鲁棒控制器设计》

3 μ综合鲁棒控制器设计

《3.1μ综合鲁棒性能定理[2]》

3.1μ综合鲁棒性能定理[2]

用结构奇异值μ定量地表征结构化不确定性对线性动态系统稳定性及性能的影响。如图2所示, 对于传递函数矩阵M表示的广义标称系统, MCn×n, 块对角阵集合Δ (s) 表示结构化不确定性摄动, w表示外部输入, z表示外部误差输出, μ可定义为

μΔ¯(Μ)={(minΔΔ¯¯{σ¯(Δ):ΔΔ¯,det(Ι-ΜΔ)=0})-10det(Ι-ΜΔ)0,ΔΔ¯

《图2》

图2 鲁棒性能分析的M-Δ结构

图2 鲁棒性能分析的M-Δ结构  

Fig.2 System performance analysis

Dolye等通过引入一个假想不确定性Δp, Δ是实际摄动, 定义:

Δq={[Δ00Δp]|ΔCq1×p1,ΔpCq2×p2}

M (s) 有q1+q2个输入和p1+p2个输出, M (s) 是P (s) 和K (s) 的线性分式变换 (LFT) ,

Μ(s)=F1(Ρ(s),Κ(s))=[Μ11Μ12Μ21Μ22](6)

为无模型摄动情况下的系统闭环阵。

如图3所示闭环系统对所有摄动Δ稳定的充要条件为μΔ (M) <1, ∀ω, ω表示频率域。

如图2, 图3所示的鲁棒性能分析的M-Δ结构框图, 鲁棒性能定理可描述为

{μ1(Μ11)<1maxμ2(FuΔBΔ¯(Μ,Δ))<1

基于鲁棒性能的4WS车辆广义控制逻辑框图如图4所示。

《图3》

图3 鲁棒稳定性要求鲁棒性能的M-Δ结构

图3 鲁棒稳定性要求鲁棒性能的M-Δ结构  

Fig.3 Stability structure for robust performance

《图4》

图4 4WS车辆广义系统框图

图4 4WS车辆广义系统框图  

Fig.4 Vehicle generalized system

《3.2μ综合鲁棒控制器设计》

3.2μ综合鲁棒控制器设计

P (s) 为广义4WS车辆系统, K (s) 为μ综合反馈控制传递矩阵, γ为传感器输出的横摆角速度, u为控制输出的控制量;dis1, dis2和 zf, zr为前后轮车辆系统不确定性模块的输入及输出;以及实际系统中不确定性的环节ΔfW11, ΔrW12, 建模中所有不确定性都归入归一化Δf, Δr, 它描述了标称模型P (s) 和实际系统模型M (s) 之间的差别。另外Δf, Δr是稳定的且仅知道有范数边界条件‖Δf<1, ‖Δr<1及其结构。W11, W12为模型不确定性频率加权函数, W2为系统性能加权函数反映在指定频率范围内系统所要求达到的性能权重。δf为前轮转角输入, e^为四轮转向车辆系统跟踪横摆角速度的控制误差信号输出。对式 (6) 表达的四轮转向闭环控制系统用LFT形式描述:{Fu (P, Δ) :ΔΔ¯,Δ<1}。四轮转向车辆μ综合控制闭环系统如图5所示。

μ综合算法是在所有稳定的控制器K中最小化闭环系统的F1 (P (s) , K (s) ) 的μΔ¯()峰值。即minΚmaxωμΔ[F1(Ρ,Κ)(jω)], 引入尺度变换矩阵D (s) , D (s) 为稳定可逆的最小相位矩阵, Dω是尺度矩阵集合中的元素, 且在每个ω处相互独立, 上式可写成minΚminDωD¯[DωF1(Ρ,Κ)(jω)Dω-1]

《图5》

图5 4WS车辆μ综合控制闭环系统

图5 4WS车辆μ综合控制闭环系统  

Fig.5 4WS vehicle closed-loop control system based on μ synthesis

对稳定的车辆系统摄动Δf, Δr, 且‖Δf<1, ‖Δr<1, 找寻一个稳定的控制器K (s) , 使得闭环车辆系统在摄动情况下仍能保持稳定, 方向盘转角跟踪横摆角速度的瞬态响应, 加权灵敏度传递函数矩阵

S(Δ)=W2(γfr-(W11Δf+Ι)Gfγ)((W12Δr+Ι)ΚGrγ+Ι)-1

有‖S (Δ) ‖<1成立。

4WS车辆系统μ控制器可描述为

x˙μ=Aμxμ+Bμu,yμ=Cμxμ+Dμu(7)

μ控制器K (s) 的设计采用D-K迭代法, 在D-K迭代中, 设计μ综合鲁棒控制器转化为寻求稳定的控制阵K (s) 和尺度矩阵D (s) , 首先保持D (s) 不变, 则上式实际上是H优化控制问题, 通过求解Riccati方程得到;保持K (s) 不变, 即在已知一稳定的控制器K (s) 前提下, 对D (s) 进行凸优化。

《4 四轮转向车辆控制系统仿真》

4 四轮转向车辆控制系统仿真

定义:Gf0=Gfγ (v0) , Gr0=Grγ (v0) , 车辆系统仿真参数如表1所示。

表1 四轮转向车辆参数

Table 1 4WS vehicle data

《表1》


m
/kg
Lf
/m
Lr
/m
Kf/
N·rad-1
Kr/
N·rad-1
v/
m·s-1
Iz/
kg·m2

1 740
1.121.5237 50075 000352 045

基于多变量不确定性的四轮转向车辆系统模型, 定义[3]

Gf=(1+ΔfW11)Gf0,Gr=(1+ΔrW12)Gr0(8)

对于车辆摄动Δf, Δr, 考虑‖Δf<1, ‖Δr<1, 把式 (8) 改写为

|Gf(v)-Gf0Gf0||W11(jω)|;|Gr(v)-Gr0Gr0||W12(jω)|(9)

式 (9) 也可表达成另外一种形式:

W11(jω){maxGf(v)-Gf0Gf0vV},W12(jω){maxGr(v)-Gr0Gr0vV}V[050](10)

根据式 (10) 及4WS车辆系统要求横摆角速度跟踪与尽量控制质心侧偏角为零的控制目标, W1, W2经过不断的试算, 现取模型不确定性加权函数W1=diag[W11, W12], Gfγ, Grγ对应的加权函数W11=W12=0.06(s+1)s/200+1;加权函数矩阵W1取对角, 且对角元素W11, W12相等, 说明车辆系统不确定性描述在某种意义上关于标称模型P是“球形”的。加权函数表明在低频处存在6 %的建模误差, 在高频处模型不确定性高达100 %, 在更高频处不确定性比重会更大。

系统的跟踪误差e^要求:e^<1/W2, 性能加权函数取W2=0.4(s+3)s+0.03, 其稳态跟踪误差<0.03/1.2=2.5 %。

模型不确定性加权函数W11, W12为乘法摄动, 如图5所示在simulink中建模, 方向盘转角输入为单位角阶跃输入, 现假设前轮转角为4°。考虑实际车辆系统乘法摄动外扰以白噪声为理论近似, 假设白噪声功率谱密度图 (PSD) 的高度为0.02 rad/s, 采样时间为0.01 s。

基于以上加权函数, 通过三次D-K迭代得到20阶控制器, 由于阶次过高, 应用MAILAB/simulink控制工具箱[4], 基于截断的均衡模型降阶, 控制器降为4阶, 经分析系统仍能维持原有的鲁棒性。鲁棒性能的满足只有在每一个频率点处闭环系统频率响应μΔ<1才行, μ分析法在ω的大部分频段上改善了σ¯ (·) 的值, 但在峰值处没有改变, 通过下一步的迭代来降低μ的峰值。这里用三阶的多项式拟合尺度矩阵D (s) , 经过三次迭代就可以得到一个性能较好的Kμ控制器, 闭环系统的鲁棒性不断改善, 直至μ<1。从表2中可以看出μ的峰值最终是0.456, 这说明用Kμ控制器的闭环系统达到了鲁棒性能的要求。

表2 D-K迭代结果

Table 2 Summary of D-K Iteration

《表2》


迭代
控制器阶数尺度矩阵D阶数γμ峰值

1
8 0 12.6880.525

2
20120.4910.475

3
20120.4570.456

图6为横摆角速度操纵稳定性响应, 可以看出, 给方向盘一个角阶跃输入, 前轮转角为4°的工况下, 4WS横摆角速度峰值为0.06 rad/s;在时间1s时基本稳定;控制器设计的目标即为横摆角速度跟踪控制。在白噪声的外扰下, 横摆角速度基本上跟踪方向盘的转角, 稳定时间也比较快, 说明γKμ控制器对系统的鲁棒性能有改善。图7为质心侧偏角响应曲线, 如图4所示, 系统反馈量为γ, 设计的是理想横摆角速度增益, 但4WS车辆质心侧偏角只是在较短时间内有小的波动, 最大峰值不到0.008 rad, 基本保持为零, 说明设计的4WS μ控制器Kμ对运动姿态也得到较好的控制。图8为侧向加速度响应曲线, 汽车侧向加速度在常速下的最大幅度一般不高于3~4 m/s2, 4WS由于后轮转向的加入, 使侧向加速度的过渡过程变得复杂;侧向加速度阶跃响应的最大幅值升高, 性能指标值的变化规律变得也较为复杂, 设计的控制器使得侧向加速度ay最大幅值为3.4 m/s2, 说明控制器很好地控制了侧向加速度性能, 响应曲线也充分证明侧向加速度反应时间迅速, 转向运动平稳, 在扰动工况下车辆仍能保持良好的动态性能。

《图6》

图6 横摆角速度响应曲线

图6 横摆角速度响应曲线  

Fig.6 Yaw rate response

《图7》

图7 质心侧偏角阶跃响应曲线

图7 质心侧偏角阶跃响应曲线  

Fig.7 Slip angle response

《图8》

图8 侧向加速度响应曲线

图8 侧向加速度响应曲线  

Fig.8 Lateral acceleration response

仿真结果证明, 通过μ综合鲁棒控制设计, 不仅可以使4WS车辆系统具有良好的操纵稳定性和鲁棒性, 而且对外界干扰具有较好的抑制性能。

《5 结语》

5 结语

笔者应用μ综合鲁棒控制设计方法, 以模型不确定性的四轮转向车辆为对象, 基于横摆角速度跟踪理论对4WS车辆控制系统进行分析, 同时考虑了系统鲁棒稳定性、鲁棒性能及系统的动态性能。结果表明, 系统具有良好的跟踪性能, 对外扰有很强的抑制作用, 改善了以往控制器设计的保守性及车辆转向响应性能, 有效地提高了车辆操纵稳定性和安全性。