最优控制理论在汽车控制系统中运用

最优控制理论在汽车控制系统中运用

董凤鸿1，张皓2

（1 北京科技大学2010级信计2班 41040317）

（2 北京科技大学2012级信计2班 41064044）

摘要：

随着人们生活水平的提高，汽车已经开始走进百姓的生活中。随着人们对汽车消费的增加，越来越多的人开始更多的关注的不仅仅是汽车本身，更多的开始关注汽车的安全性及舒适性。由此，各大汽车厂商更具消费者的需求开始着重研究带有主动控制能力的汽车控制系统。本文引入最优控制理论对当今比较流行的汽车悬挂系统、汽车防抱制动系统（简称ABS 系统）和无级变速器控制系统进行优化。由此达到优化汽车安全性、经济性和舒适性。

关键词：

最优控制理论、悬挂系统、防抱制动系统、无级变速器控制系统

一、引言

汽车防抱制动系统(简称ABS系统) ，实质上是一种制动力的自动调节装置。这种装置使汽车制动系统的结构发生了质的变化,它不仅能充分发挥制动器的制动性能,提高制动减速度和缩短制动距离,而且能有效地提高汽车制动时的方向稳定性,大大改善汽车的行驶安全性。悬挂系统是指车身与车轴之间连接的所有组合体零件的总称，悬挂系统直接影响着汽车的安全性、稳定性和舒适性，是汽车的重要组成部分之一。目前，降低汽车能源消耗和减少废气排放已成为汽车行业最关注的问题，大量试验表明，装有无级变速器(CVT)的汽车比装有传统有级变速器的汽车在改善汽车燃油经济性和排放等方面具有更大的潜力，这是因为CVT连续变化的传动比可以使发动机转速独立于负载和车速的变化，最大限度地发挥发动机的经济性和动力性。

二、正文

（一）、汽车防抱制动系统最优控制

1、方法介绍

最优控制是基于状态空间法的现代控制理论方法。它可以根据车辆一地面系统的数学模型,用状态空间的概念,在时间域内研究汽车防抱制动系统。是一种基于模型分析型的控制系统,它根据防抱系统的各项控制要求,按最优化原理求得控制系统的最优控制指标。我们知道:现代控制理论应用得成功与否,关键在于数学模型是否准确。为此必须首先研究用状态变量表示的防抱系统的数学模型。

2、模型建立

为了便于分析首先作如下假设：

(1)车轮承受的载荷为常数;

(2)不计迎风阻力和滚动阻力;

(3)附着系数Ψ随滑动率变化的规律由图1所示的两条直线所组成 ,

其数学表达式为:

ψ=ψk S T ?S S

式中

Ψh 一最大附着系数

ST 一最大附着系数对应的滑动率

Ψg 一车轮完全抱死时(S = l )的附着系数

S = 1-RW/V 一滑动率

汽车在制动过程中,单个车轮的受力情况如图2所示。根据理论力学的知识和以上假设，可写出车轮作平面运动的运动微分方程式：

MV =F t ?F F =?Z ?ψ S

JW =RF t ?M f ?M b =?M b +Z ?ψ S ?R

图1 Ψ-S 近似曲线图2 制动时车轮受力图

式中

M ———分配到车轮上的汽车质量 J ———车轮的转动惯量 V ———车体的速度W ———车轮的角速度Mf ———地面对车轮的滚动阻力矩

Mb ———制动器制动力矩F =Z Ψ（S ）——地面对车轮的水平作用力

Z ———车轮的法向反力Ψ S ———附着系数，它是滑动率S 的函数Ff ———车体受到的迎风阻力R ———车轮半径

一般地,可把制动力矩表示为如下的时间函数关系

Mb =a ·P （t ）

根据现代控制理论的要求，除需要选取车轮角速度W 和角加速度W 为状态变量外，为了产生闭环控制系统,还应把附着系数Ψ和滑移率S 的关系曲线峰值处的车轮速度V *作为系统的期望值输出,显然它在制动过程中是随时间变化的,因而需要设计跟踪系统,使系统实际输出的是跟踪期望输出值 ,于是可将跟踪输出器设计成二阶积分的型式

,

即 :

I φ1′= ωR ?v ? dt = V ω?V ? t

0dt t 0I φ2′= I φ1dt = V ω?V ? t 10

dt t 0t 0

式中： WR ———车轮的速度

此两式可以写成如下的微分关系：

I φ1I φ2 = ωR ωR + ?V ??V

? 或写成状态变量Ir1，Ir2的直接关系

I φ1I φ2 = RωI ?φ1 + ?V ?0 在研究中为了便于车体速度V 和峰值车轮速度V* 相比较,将车轮角速度和角加速度两个状态变量 ,用车轮的速度Vw=Rw 和加速度Vw=Rw 来代替,作为系统的状态变量,则可得

R ω=?R 2Z VJ ψn S T ?αi J

p (t )

另外

V ′ω=Rω′

联合写成矩阵形式

V ωV ′ωI φ1′I φ2′ = ?R 2Z ψn T 0 0 010 0 0 0 1 0 0 0 0 1 0

V ′ωV ωI φ1I φ2 + α′J 000 + 00?10 V ?

由此可得汽车防抱制动系统状态方程的规范表达式：

X ′=AX +BU +Bd Y =CX

式中

A = ?R 2Z VJ ψn S T 0 0 010 0 0 0 1 0 0 0 0 1 0

系数矩阵

B = 0 0 ?1 0 T ?误差矩阵B = ?αi J 0 0 0 T ?控制矩阵X = V W ′V W I φ1I φ2 T ?状态变量

C = 0 1 0 0 ?输出矩阵Y =ωR =V W ?输出向量U =p t ?控制向量d f =V ? ?误差向量

用现代控制理论的方法设计汽车防抱制动装置,实质上就是设计一个最优控制系统,使 其在防抱的全过程中能预报出一种控制函数,使防抱系统在防抱过程中以最优的方式工作,

使预先设定的目标函数达到最小值。为了使液压或气压控制系统消耗的能量最小,并使实际输出与期望输出的误差最小,我们选择具有二次型的目标函数，即：

J r = X T θX ?U T RU dt ω

式中θ—状态变量的加权矩阵

R —控制变量的加权矩阵

图3 线性调节器方框图

对于这种线性状态方程的系统,其性能指标为状态变量和控制变量的二次型函数的最优控制问题 ,称为线性自动调节问题。其系统结构可用图3所示的方框图来表示,由最优控制理论不难求出该系统的最优控制规律 :

U dt =?R ?1B T LX =?KX

其中K=R-1BTL 称为反馈控制的线性反馈系数，L 则可Riccati 方程来求得

?LA ?A T T +LBR ?1B T L ?θ=0

代入 K=R-1BTL 则可求出反馈控制系数K=[K1 K2 K3 K4]

（二）、汽车防抱制动系统最优控制

1、方法介绍

本文对主动悬挂系统和被动悬挂系统的特点进行了对比分析，并进一步以主动悬挂系统为研究对象，结合车辆动力学原理，推导了一种基于最优控制理论的主动悬架控制方法。控制方法是主动悬挂系统的核心技术之一，国内外学者提出了自适应控制、预见控制、滑模控制、自校正控制、最优控制理论、模糊控制和神经网络控制等方法。其中最优控制理论基础比较完善，其最大优点是不必根据要求的性能指标确定系统闭环极点的位置，只需根据系统的响应曲线找出合适的状态变量和控制变量的加权矩阵，使系统性能指标函数即目标函数

J

最小。

2、模型建立

主动悬挂系统的状态方程大多具有线性形式：

x′=Ax+Bv

式中：A为，n x n的系统矩阵；B为，n x r控制矩阵；x为n维状态矢量；v为r维控制矢量。加约束后性能指标函数J，为：

J=

1

[x T Qx+v T Rv+2x T Nv]dt t2

t1

式中：Q为状态变量的加权矩阵；R为控制变量的加权矩阵；N为交叉项的权重。这里应注意：要求系统为线性定常系统，且要求系统完全能控；优化后的闭环系统是渐近稳定的。

由于目前的悬挂形式主要是独立悬挂，而二自由度1／4车辆模型能较好地描述汽车独立悬挂系统的实际情况，故取1／4车辆模型作为主动悬挂系统优化控制的研究对象，模型如图4所示。

图4 主动悬挂系统结构示意图

图4中m2为1/4车身质量；m1为车轮质量；u=-Kx表示悬挂对车身或车轮的作用力；K为n 维向量；k1为车轮弹性系数，车轮阻尼由于影响不大而忽略不计，y2、y1、y0分别表示1/4车身质量、车轮质量的位移和路面的激励。

根据牛顿第二定律，可写出系统的微分方程式：

m1y′′=k1y0?y1?u

m2y′′

2

=u

选取状态变量为：

x1=y1?y0

x2=y2?y0

x3=y1′

x4=y2′

由上式得：

x1‘=x3?y0′

x2′=x42?y0

x3′=?k1

m1

x1?

u

m1

x4′=

u m2

由此可得出系统的状态方程为:

x′=Ax+Bu+Eω式中：

A=

00 1 0 10 0 1 k1

m1

0 0 0 0 0 0 0

B=

?1

m1

1

m1

,E=

?1

?1

E为扰动输入矩阵。

在确定目标函数时，应考虑到汽车的平顺性和操纵稳定性，悬挂系统弹簧的动扰度(Y1—y2)会影响到汽车的平顺性，且车轮与路面间的动载荷会影响汽车操控稳定性；从实现控制的角度看，应使所需的控制能量较小，综合以上几种因素，目标函数可写为:

J=[Ru2+q1(y0?y1)+q2(y2?y1)2]dt

∞

或可写成:

J=[x T Qx+Ru2]dt

∞

经验证系统是能控的，其中:

A=q10 1 0 1q2 0 1 0 0 0 0 0 0 0 0

根据二次型性能指标的线性系统最优控制理论，其最优控制律为：

u=?Kx=?R?1B T P

式中：K为最优反馈增益矩阵；P为实对称常值矩

阵，满足黎卡提(Riccati)代数方程，即

?PA?A T P+PBR?1B T P?Q=0

其中矩阵Q的大小与轮胎位移加权系数q1和悬挂动扰度加权系数有关，q1和口2取不同的值就允许对不同的分量加不同的权系数，当某个分量需要特别约束时，可以增大该分量的加权系数。

（三）、汽车防抱制动系统最优控制

1、方法介绍

CVT变速系统是机电液一体化的非线性系统，在升档、降档等不同情况下，其数学模型具有不同的特征，采用固定控制规律和参数不能适应各种不同的工况。由此，本文设计出了电液控制系统的最优控制器。

2、模型建立

cvT电液控制系统可以分为夹紧力控制系统和速比控制系统，控制系统原理（如图5所示）。夹紧力控制系统主要由夹紧力控制阀构成，它控制着液压系统压力的大小，这个压力直接作用于从动轮油缸。在速比一定的条件下，它的大小决定了系统传递转矩的能力。速比控制系

统主要由速比控制阀构成，通过该阀可控制主动轮油缸的充油量，从而改变主动轮油缸内的压力，在金属带的制约下，使主动轮组中可移动带轮的轴向位置发生变化，实现变速比的控制。

根据推导，可以建立速比控制系统的方框图(如图6所示)。因为比例伺服阀的固有频率通常比动力元件固有频率高一倍左右，为简化起见，在此把其视为比例环节。

图5 电液控制系统

图6 速比控制系统传递函数模型方框图

图中：

u为输入电压信号， U f为反馈电压， U e为偏差，I 为力矩马达输入电流，K e为为比例控制器系数，K s为伺服阀的增益，Q(s)为流过电液比例伺服阀的流量，A p为主动轮油缸活塞的有效面积，βe为体积弹性模量，M t为活塞以及与活塞相连的负载总质量，B p为负载的粘性阻尼系数，V1为主动轮油缸的体积，F DR为主动油缸轴向夹紧力，H y为电位计转换系数，Y 为主动轮油缸的输出位移。

为简化系统模型，我通过阅读有关数据资料给出开环传递函数G（s）为：

G s=150

(1)

s1?10?4+1?10?2s+1

3、最优控制设计

设状态空间方程形式为：

X′=AX+BU

Y=CX(2)

式中，A、B、C为空间变量系数矩阵。根据控制理论知识，将（1）式的传递函数转换为状态空间方程：

X1′

X2′X3′=

010

001

0?10000?100

X1

X2

X3

+

1500000 Y= 1 0 0[x1x2x3]T(3)

式中Y和U分别为状态空间的输出和输入变量，Y=[x1x2x3]T为状态向量。

利用最优控制理论，在系统中加入最佳控制信号U(t)=-K（t）+X（t），其中K=R?1B T P,p为Riccati方程的解，并使二次型性能指标取极小值。由此看出，控制系统的优化设计实质上是求最佳调节器K（t）的设计。为了求解K（t），即解Riccati方程：

?PA?A T P+PBR?1B T P?C T QC=0

三、结论

A、在汽车防抱制动系统问题的研究中，可知ABS是一种简单实用的系统,其质量和功能在不断地完善。本文所介绍的最优控制法在理论上很成熟,它将车轮的角速度和角加速度作为状态变量对系统进行优化控制,能达到很好的防抱制动效果,但它要求建立较高质量的数学模型。

B、在汽车悬挂系统问题中，最优控制理论知识同汽车控制系统结合较好并可知主动悬挂系统很大程度上降低了车身的加速度，增加通过性。对二自由度1／4主动悬挂系统的分析与速度，提高了乘坐的舒适性。由于主动控制力的存在，悬挂系统的动行程变化平缓，且比被动悬挂系统行程小，许用的悬挂空间得到了充分利用。

C、通过对无级变速器控制系统研究，可以知道最优控制理论的可运用广空间的广泛性。

综上，最优控制理论在汽车控制系统中的运用可以大大提高汽车的经济性、舒适性和安全性。它对于未来汽车控制系统研究是一种科学、实用的工具理论，有很大的运用发展空间。

四、致谢

感谢廖福成老师一学期辛苦的付出，以及孜孜不倦的教导。由此，让我们学会了很多论文研究知识，同时也是廖老师的现代控制理论课程给本次最优化控制理论在汽车控制系统中的研究论文提供了理论基础。在此，特表感谢！

参考文献：

[1] 何渝生.汽车控制理论基础及应用.重庆：重庆大学出版社，1995

[2] 张孝祖.车辆控制理论基础及应用.北京：化学工业出版社，2007

[3]贾启芬，刘习军．机械与结构振动[M]．天津：天津大学出社，2006·

[4]王家春、董申、李旦.超紧密机床的主动隔振系统研究[J].振动与冲击，

2000,19(3):54-56·

[5]李海波、何天明．汽车半主动悬架系统的研究现状及趋势[J].北京汽车，2007（3）：22-25·

[6]胡建军、秦大同等.金属带式无极变速传动比变化特性研究[J].汽车工程，2003,25（1）：25-29·

[7]李泳龙、张伯英、董秀国等.金属带式无极变速器液压控制系统[J]吉林大学学报，2000,30（3）：14-19·

[8] 王春行.液压伺服控制系统[M].北京：机械工业出版社

[9] 王占林.近代液压控制[M]北京：机械工业出版社

[10]魏巍.Matlab控制工具箱技术手册[M].北京：国防工业出版社

[11]吴光强，杨伟斌，秦大同．双离合器式自动变速器控制系统的关键技术[J].机械工程学报，2007，43(2)：13-21．

[12]牛铭奎，程秀生，等．双离合器式自动变速器换挡特性研究[J]．汽车工程，2004，26(4)：453—457．

[13]金伦，程秀生，等．双离合器自动变速器仿真研究[J]．汽车技术，2005(8)：4-7

[14]PietroDolcini，Hubert B’echart．Observer-based Optimal Control ofdry Clutch Engagement[C]．Proc．of 44“IEEE Conference onDecision and Control，and the European Control Conference 2005：440—445．

[15]Luigi Glielmo，Francesco Vasca．Optimal Control of Dry Clutch

Engagement[c]．SAE Paper 2000—01-0837．

[16] 邢继祥．最优控制应用基础[M]．北京：科学技术出版社，2003．

[17] 柴山、孙焕纯．离散变量结构优化设计的组合算法[J]．应用力学和数学，1997，(9)

[18]朱浩鹈、李为吉．结构多目标优化非劣迹遗传算法．西北工业大学学报，2001，19(1)

[19]王馥．智能优化算法及其应用[M]．清华大学出版社，200l

[20]王小平，曹立明．遗传算法一理论、应用与软件实现[M]．西安文通大学出版社，2002

[21]陈翠彪，刘正士，彭敏，等．双轴汽车受路面作用的动载荷分析[J]．农业机械学报，2008，39(2)：185～187．

[22]曹源文，梁乃兴，于清，等．路面不平整引起的车辆动载计算方法[J]．交通运输工程学报，2008，8(2)：69--73．

[23] 陈潇凯，林逸，施国标．等．多体系统动力学软件在汽车工程中应用的新进展[J]．计算机仿真，2005，22(6)：201--204．

[24] 缪炳荣，方向华，傅秀通．SIMPACK动力学分析基础教程[M]．成都：西南交通大学出版社，2008．

[26] 喻凡，林逸．汽车系统动力学[M]．北京：机械工业出版社，2005：34～149．