车辆主动悬架最优控制.doc

《车辆主动悬架最优控制.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、车辆主动悬架的控制研究悬架是汽车的重要装置之一，它对汽车的平顺性、操纵稳定性、通过性等多种使用性能有着很大的影响。设计优良的悬架系统，对提高汽车产品质量有着极其重要的意义。目前，汽车上普遍采用的是弹性元件和减震器组成的常规悬架，从控制力学的角度，将这种悬架称为被动悬架。实践和研究结果都表明，常规悬架受到许多限制，即使采用优化方法来设计也只是将其性能改善到一定程度。为了克服常规悬架对其性能改善的限制，在汽车中采用和发展了新型的主动悬架。主动悬架能够根据路面情况及汽车运行的实际状态进行最优反馈控制，使汽车整体行驶性能达到最佳。主动悬架的主要特点是能够主动提供能量，与

2、传统被动悬架相比，其最大的优点在于具有高度的自适应性。一、车辆主动悬架系统建模主动悬架的分析模型如图3.3所示，图中u为主动悬架执行机构的作用力。主动悬架的运动微分方程为：(1)状态变量、输出向量的选取同被动悬架，且为了便于与被动悬架的比较分析，选取与被动悬架模型相同的输入信号，路面激励仍为选白噪声，根据微分方程组（1），建立如下所示的状态方程和输出方程（2）式中：；；；；汽车悬架可认为是一种连续线性的随机最优控制系统，由最优线性滤波器串接确定性调节器的最优反馈增益系数矩阵组成。这两部分参数可分别加以确定。对于控制要求的性能指标是二次函数积分型的调节器问题，外界

3、干扰是高斯白噪声，综合性能指标为：(3)此处认为汽车主动悬架的最优控制器为一个终端时间无限的线性调节器，问题仍是寻找最优控制，使目标函数J取极小。线性调节器的主要问题之一是如何选择Q、R阵以获得比较满意的控制过程动态响应，计算机仿真可以解决这个问题。在悬架设计中，为提高汽车的操纵稳定性和行驶平顺性，应使簧载质量垂直加速度、悬架动扰度及轮胎动变形较小。此外，从实现控制的角度来看，应使所需的控制能量较小。因此式(3)可写为(4)或写为(5)其中这里，q1——轮胎动变形加权系数q2——悬架动扰度加权系数式(3.9)中第一、二项为误差指标，表示在～∞整个时间系统实际状态

4、与平衡之间的误差总和。这一积分越小，说明控制误差越小，性能越好。积分式中第三项为能量指标，表示在～∞整个时间支付能量的总和。系统状态转移是考控制u(t)来进行的，为要使系统误差很小，则需要支付很大的能量代价。最优反馈增益系数矩阵式可写成(6)式中，增益值k1~k4有明确的物理意义。k1可等效于一放置于簧载和非簧载质量间的弹簧，改变k1则影响簧载质量的固有频率；k2作用于簧载质量的绝对速度上，影响其悬挂阻尼；k3大小涉及轮胎变形，对车轮的垂直弹跳频率产生影响；k4作用于非簧载质量的速度上，影响其非悬挂阻尼。二、主动悬架系统的能控性，能观测性能控性和能观测性是系统的

5、一种特性，是现代控制理论中的两个基本概念。状态完全能控的充分必要条件是能控矩阵满秩；状态完全能观测的充分必要条件是能观测矩阵满秩。主动悬架系统参数的选取如下，即m1=36kg；m2=240kg；kt=160000N/m；将参数值带入矩阵，利用Matlab中的函数Co＝ctrb(A,B)求悬架系统的能控矩阵Co，利用函数rA=rank(Co)得矩阵的秩为rA=4,满秩，故系统是能控的。利用Matlab中的函数Ob＝obsv(A,C)求悬架系统的能观测矩阵Ob，利用函数rB＝rank(Ob)得矩阵的秩为rB=4,满秩，故系统是能观测的。三、主动悬架的频域仿真为了求得

6、主动悬架系统的最优控制u(t),必须先求得反馈增益矩阵K，而K矩阵的求解决定于黎卡提代数方程的解—P矩阵，这可以用计算机来实现求解。程序用Matlab语言编写，给定一组矩阵A、B、Q、R的有关数据，经过计算，便可以最终得到相应的矩阵P、K的数值。下面取三组不同权系数q1，q2进行计算分析；1)取q1=3.35E5,q2=40.5E5时，由程序得k1=2012.5,k2=977.1,k3=－1874.8,k4=－31.3，并求得系统的传递函数及幅频特性，绘制系统的幅频特性图%主动悬架q1=3.35e5;q2=40.5e5时的仿真程序：m1=36;m2=240;kt

7、=160000;q1=3.35e5;q2=40.5e5;A=[010-1;0000;000-1;00kt/m10];B=[0;1/m2;0;-1/m1];D=[0;0;1;0];C=[0000;1000;0010];E=[1/m2;0;0];H=[0;0;0];Q=[q2000;0000;00q10;0000];R=[1];[K,P,F]=lqr(A,B,Q,R)M=A-B*K;N=C-E*K;G=ss(M,D,N,H);G1=tf(G)i=1;fors=0:0.1:80s=s*2*pi*j;G11=(7.811*s^3+580.4*s^2+3.727e004*

8、s+1.422e-)/(