主动悬架LQG最优控制设计

主动悬架LQG最优控制设计吕福麟113085234379摘要：根据汽车行驶性能的要求，建立二自由度的1/4汽车动力学模型，利用最优控制理论对主动悬架的LQG（LinearQuadraticGaussian）控制器进行设计，运用MATLAB/simulink对模型仿真，对比主动悬架与被动悬架在控制效果上的差别。仿真结果表明，具有LQG控制器的主动悬架可以明显的提高汽车行驶的操稳性能。关键词：主动悬架；LQG控制器；MATLAB/simulink；仿真结果ABSTRACT：Accordingtotherequirementofthevehicledrivingperformance,a2-degree-of-freedoms1/4cardynamicmodelwasbuild.LQG(LinearQuadraticGaussian)controllerforactivesuspensionwasdesignedwiththeoptimalcontrollaw,MATLAB/simulinkwasusedtosimulationthemodel,comparedthedifferenceofcontrolperformancewithactivesuspensionandpassivesuspension.TheresultofsimulationshowthatitcouldimprovethecardrivinghandlingandstabilitywithLQGcontroller.KeyWords:activesuspension,LQGcontroller,MATLAB/simulink,simulationresult1前言传统的悬架主要由弹性元件、减震器和导向装置组成，他们的阻尼和刚度已经确定就不便于调节，而且只能在特定的路面激励和特定的车速下才能达到最优控制，灵活性较差。不能满足人们对驾驶舒适性和操纵稳定性的要求。为了解决这个问题，主动悬架控制渐渐得到发展，他与被动悬架的主要区别在于可以根据不同路面激励和行驶状况，自行调节车辆的动态，从而满足人们对行驶的要求。本问就是在此基础之上研究主动悬架的最优控制。22自由度1/4车辆模型的建立图1单轮车辆模型图中mb为车身质量；mw为车轮质量；xb为车身位移；xw为车轮位移；xg为路面输入；Ks为悬架刚度；Kt为轮胎刚度；U为主控制力。对于图1建立1/4车辆运动微分方程：UxxKxmwbsbb)(UxxKxxKxmwgtwbsww)()(2.1路面模型的建立在分析主动悬架控制过程时，路面输入是一个不可忽略的重要因素，本文利用白噪声信号为路面输入激励，)(2)(2)(000twUGtxftxgg其中，0f为下截止频率，Hz；G0为路面不平度系数，m3/cycle；U0为前进车速，m/sec；w为均值为零的随机输入单位白噪声。上式表明，路面位移可以表示为一随机滤波白噪声信号。这种表示方式来源于试验所测得的路面不平度功率谱密度（PSD）曲线的形状。我们可以将路面输入以状态方程的形式加到模型中：XCYWFXAXroadroadroadroadroad1,2,2,000roadroadroadgroadCUGBfAxX；考虑路面为普通路面，路面不平系数G0=5e-6m3/cycle；车速U0=20m/s；建模中，路面随机白噪声可以用随机数产生（RandomNumber）或者有限带宽白噪声（Band-LimitedWhiteNoise）来生成。本文选择有限带宽白噪声为路面输入激励，运用MATLAB/simulink建立仿真模型如下：图2路面模型仿真构造出的随机路面轮廓如下图所示：图3路面输入仿真结果2.2LQG控制器的设计在汽车悬架设计中，选取的状态变量为：TwbwbxxxxX],,,[则可以根据车辆动力学方程将这些变量写成状态方程空间矩阵的形式：FWAXX其中A为状态矩阵；F为输入矩阵，其值如下：001000010000wstwsbsbsmKKmKmKmKA000wtmKF主要的性能指标有：车身垂直加速度bx；轮胎动变形)(gwxx；悬架动行程)(wbxx；即：wbgwbxxSWSxxDTDxBA:::将这三个性能指标写成矩阵形式：TwbgwbxxxxxY)](),(,[输出矩阵：1100100000bsbsMkMkC转移矩阵：010D由此建立1/4车辆被动悬架动力学模型：图41/4车辆被动悬架动力学模型2.3性能指标的确定LQG控制设计中的目标性能指数J即车身加速度、悬架动行程和位移的加权平方和的积分值，表示如下：dtxqxxqxxqTJTbwbgwT0232221])()([1lim为了据此求解状态反馈增益，必须用状态变量以及输入变量来表示上式：dtNUXRUUQXXTJTTTTT0]21lim其中Q：对应于状态变量的权重矩阵；R：约束输入信号大小的权重矩阵；N：耦合项。对于q1，q2，q3表达的性能函数，可以整理为：dtxqxxqxxqTJTbwbgwT0232221])()([1limdtxxxxxQxxxxxTwbgwbTwbgwbT00][1lim其中2130000000qqqQ由于DUCXY所以CXQDUUDQCXUDQDUXCQCXDUCXQDUCXYQYTTNTTRTTQTTTT000000)()(式中q1，q2，q3分别为轮胎位移，悬架动行程，车身垂直加速度的加权系数。因此可以求出Q、R、N：11122212222222220000000000000000000qqqmKqqmKqmKqmKqCQCQbSbSbSbST；bTmDQDR10；0000ssTKKDQCN；利用MATLAB函数LQR来计算状态反馈K：),,,,(][NRQBAlqrESK其中，K为最优状态反馈矩阵；S为Riccati方程解；E为系统特征值。因此可以得到主动控制力U=-KX。)(54321gwbwbxkxkxkxkxkU加权系数的选取决定悬架的性能，如果车身垂直加速度加权系数较大，则可以提高乘坐的舒适性；若轮胎动位移的加权系数较大，则车辆的操纵稳定性较好。3仿真分析根据ISO标准，设定汽车动力学参数为：车身重量kgmb320；轮胎重量kgmw40；悬架弹簧刚度mNKs/1024；轮胎等效刚度stKK10；悬架等效阻尼mNsCs/1000；路面不平度系数cyclemeG/6530；前进速度hkmU/200；下截止频率1.00f；所选加权系数05.0,1000,100321qqq。求的得反馈增益矩阵]8502.25846.05255.20641.04946.0[K。用simulink建立的主被动悬架的模型如下：图51/4车辆主被动悬架动力学模型仿真结果下图所示：图6车身加速度对比图7轮胎动行程对比图8悬架行程对比图9主动悬架与被动悬架的幅频对比4总结本文首先建立了二自由度1／4车体模型。并运用线性二次最优化理论设计了车辆主动悬架的LQG控制器。借助MATLAB／Simulink软件进行仿真分析。分析结果表明，所设计的最优主动悬架显著地降低了车身的垂向振动加速度。与被动悬架相比，采用LQG控制器的悬架动行程和车轮动位移也都得到了良好的改善。所以基于线性二次最优控制理论进行的主动悬架LQG控制器的设计是行之有效的。