主动悬架LQG最优控制设计

1、主动悬架LQG最优控制设计吕福麟113085234379摘要：根据汽车行驶性能的要求，建立二自由度的1/4汽车动力学模型，利用最优控制理论对主动悬架的LQG(LinearQuadraticGaussian控制器进行设计，运用MATLAB/simulink对模型仿真，对比主动悬架与被动悬架在控制效果上的差别。仿真结果表明，具有LQG控制器的主动悬架可以明显的提高汽车行驶的操稳性能。关键词：主动悬架；LQG控制器；MATLAB/simulink;仿真结果ABSTRACT:Accordingtotherequirementofthevehicledrivingperformance2-degree-

2、of-freedoms1/4cardynamicmodelwasbuild.LQG(LinearQuadraticGaussian)controllerforactivesuspensionwasdesignedwiththeoptimalcontrollaw,MATLAB/simulinkwasusedtosimulationthemodel,comparedthedifferenceofcontrolperformancewithactivesuspensionandpassivesuspension.Theresultofsimulationshowthatitcouldimprovet

3、hecardrivinghandlingandstabilitywithLQGcontroller.KeyWords:activesuspension,LQGcontroller,MATLAB/simulink,simulationresult1 1 前言传统的悬架主要由弹性元件、减震器和导向装置组成，他们的阻尼和刚度已经确定就不便于调节，而且只能在特定的路面激励和特定的车速下才能达到最优控制，灵活性较差。不能满足人们对驾驶舒适性和操纵稳定性的要求。为了解决这个问题，主动悬架控制渐渐得到发展，他与被动悬架的主要区别在于可以根据不同路面激励和行驶状况，自行调节车辆的动态，从而满足人们对行驶的要求

4、。本问就是在此基础之上研究主动悬架的最优控制。22 自由度 1/4 车辆模型的建立图 1 单轮车辆模型图中mb为车身质量；mw为车轮质量；xb为车身位移；xw为车轮位移；xg为路面输入；Ks为悬架刚度；Kt为轮胎刚度；U为主控制力。对于图1建立1/4车辆运动微分方程：mbXb二一 Ksd-Xw)UmwXw=Ks(Xb-Xw)Kt(Xg-Xw)-U2.1 路面模型的建立在分析主动悬架控制过程时，路面输入是一个不可忽略的重要因素，本文利用白噪声信号为路面输入激励，X：(t)=-2 二 foXg(t)2 二.GoUow(t)其中，f0为下截止频率，Hz；Go为路面不平度系数，m3/cycle；Uo为

5、前进车速，m/secw为均值为零的随机输入单位白噪声。上式表明，路面位移可以表Xroad=AoadX+FroadWYroad=CroadXXroad=Xg,Aoad=一2话,B,oad=%3。,%=1；考虑路面为普通路面，路一一一一一，一 3面不平系数Go=5e-6m/cycle；车速Uo=2Om/s；建模中，路面随机白噪严可以用示为一随机滤波白噪声信号这种表示方式来源于试验所测得的路面不平度功率谱密度(PSD)曲线的形状。我们可以将路面输入以状态方程的形式加到模型中:随机数产生(RandomNumber)或者有限带宽白噪严(Band-LimitedWhiteNoise)来生成。本文选择有限带

6、宽白噪声为路面输入激励，运用MATLAB/simulink建立仿真模型如下：图 2 路面模型仿真构造出的随机路面轮廓如下图所示：2.2 LQG 控制器的设计在汽车悬架设计中，选取的状态变量为:X=Xb,Xw,Xb,XwT则可以根据车辆动力学方程将这些变量写成状态方程空间矩阵的形式:X=AXFW其中A为状态矩阵；F为输入矩阵，其值如下:roadroadmodsallBard-LimitedWhiteNoi&e1(tMUBt6oe-d老图 3 路面输入仿真结果KsKsmbmbKs-Kt-Ksmwmw0000_-0KtF=mw0：。主要的性能指标有：车身垂直加速度 Xb；轮胎动变形（Xw-Xg）（X

7、b-Xw）；即:BA：XbDTD:Xw-XgSWS:Xb-Xw将这三个性能指标写成矩阵形式：、，TY=Xb,（Xw-Xg）,（Xb-Xw）输出矩阵:ksksMbMb011-1转移矩阵:一 0D=1-0J由此建立1/4车辆被动悬架动力学模型:00A=0010101；悬架动行程图 41/4 车辆被动悬架动力学模型2.3 性能指标的确定平方和的积分值，表示如下:为了据此求解状态反馈增益，必须用状态变量以及输入变量来表示上式:1TT_T_T_J=lim-oXQXURU2XNUdt其中Q:对应于状态变量的权重矩阵；R约束输入信号大小的权重矩阵；N:耦合项。对于qi,q2,q3表达的性能函数，可以整理为:

8、2,1Tr,、2,、2,J=ljmT(qJXwXg)+q2(XbXw)+q3XbdtXbXw-XgIXb-Xw其中LQG控制设计中的目标性能指数J即车身加速度、悬架动行程和位移的加权1TJ=ljmT0q1(x、2/、2w-Xg)q2(Xb-Xw)2q3XbdtXb-xwQodtQ0q30-0qi0oloq2YTQ0Y=(CXDU)TQ0(CXDU)QRNTTTTTTTTCQoCXUDQoDUXCQoDUUDQO因此可以求出Q、R、N:一 0o-Ks|KsoJ利用MATLAB函数LQR来计算状态反馈K:KSE=lqr(A,B,Q,R,N)其中，K为最优状态反馈矩阵；S为Riccati方程解；E为

9、系统特征值。因此可以得到主动控制力U=-KX0U=-(kiXbk2Xwk3Xbk4Xwksxg)加权系数的选取决定悬架的性能，如果车身垂直加速度加权系数较大，则可以提高乘坐的舒适性；若轮胎动位移的加权系数较大，则车辆的操纵稳定性较好。3 仿真分析根据ISO标准，设定汽车动力学参数为：车身重量mb=320kg；轮胎重量mw=40kg；X由于Y=CXDU所以式中qi,q2,q3分别为轮胎位移,悬架动行程,车身垂直加速度的加权系数。Q=CTQOC=q2ooKSooKS一q2-2晨-q2一-2mb0qimbKSq22mb-qi_qqiR=DTQ0D=-；N=CTQODmb悬架弹簧刚度Ks=2父io4N

10、/m;轮胎等效刚度Kt=iOKs;悬架等效阻尼Cs=iOOONs/m；路面不平度系数Go=5e-6m3/cycle;前进速度Uo=20km/h;下截止频率f0=0.1;所选加权系数q=100,q?=1000,%=0.05。求的得反馈增益矩阵 K=0.4946-0.06412.52550.58462.8502。用simulink建立的主被动悬架的模型如下:图51/4车辆主被动悬架动力学模型仿真结果下图所示:Displayroadmodel1/4vehlicedynamicsToworkspacesTaworkspace4口1/4augmentaynamiesswsM=Ax+Buy=Cx+Du始=Ax+0uy=CK+DUClockkr=AxtBury-CK+DUB-Band-LimitedwniieNoiseDisplay1ToWorkspace1图 6 车身加速度对比图 7 轮胎动行程对比图 8 悬架行程对比图 9 主动悬架与被动悬架的幅频对比4 总结本文首先建立了二自由度1/4车体模型。