《车辆主动悬架最优控制》由会员分享，可在线阅读，更多相关《车辆主动悬架最优控制（15页珍藏版）》请在金锄头文库上搜索。
1、车辆主动悬架的控制研究悬架是汽车的重要装置之一，它对汽车的平顺性、操纵稳定性、 通过性等多种使用性能有着很大的影响。设计优良的悬架系统，对提高汽车产品质量有着极其重要的意义。目前，汽车上普遍采用的是弹性元件和减震器组成的常规悬架，从控制力学的角度，将这种悬架称为被动悬架。 实践和研究结果都表明，常规悬架受到许多限制，即使采用优化方法来设计也只是将其性能改善到一定程度。为了克服常规悬架对其性能改善的限制，在汽车中采用和发展了新型的主动悬架。主动悬架能够根据路面情况及汽车运行的实际状态进行最优反馈控制，使汽车整体行驶性能达到最佳。主动悬架的主要特点是能够主动提供能量，与传统被动悬架相比，其最大的优
2、点在于具有高度的自适应性。一、 车辆主动悬架系统建模主动悬架的分析模型如图3.3 所示，图中u 为主动悬架执行机构的作用力。主动悬架的运动微分方程为：.m2 x2u.m1 x1u kt ( x1x0 )(1)状态变量、输出向量的选取同被动悬架，且为了便于与被动悬架的比较分析，选取与被动悬架模型相同的输入信号，路面激励仍为选白噪声(t) ，根据微分方程组（1），建立如下所示的状态方程和输出方程。xAxBuD(t )yCxEu（ 2）式中：010101100000m2m2AB0000000010E0DC100000kt0110m1；m1 ；0 ；0 0 10 ；汽车悬架可认为是一种连续线性的随机最
3、优控制系统，由最优线性滤波器串接确定性调节器的最优反馈增益系数矩阵组成。这两部分参数可分别加以确定。对于控制要求的性能指标是二次函数积分型的调节器问题，外界干扰是高斯白噪声，综合性能指标为：J (u) X T (t )Q(t) X (t)u T (t ) R(t )u(t )dt(3)0此处认为汽车主动悬架的最优控制器为一个终端时间无限的线性调节器，问题仍是寻找最优控制 u(t) ，使目标函数 J 取极小。线性调节器的主要问题之一是如何选择Q、R 阵以获得比较满意的控制过程动态响应，计算机仿真可以解决这个问题。在悬架设计中， 为提高汽车的操纵稳定性和行驶平顺性， 应使簧载质量垂直加速度、 悬架
4、动扰度及轮胎动变形较小。此外， 从实现控制的角度来看， 应使所需的控制能量较小。因此式 (3) 可写为Jq1 ( x1x0 ) 2q2 ( x2 x1 ) 2Ru 2 dt0(4)JXTQXRu 2 dt(5)或写为0q2000Q000000q10其中0000这里， q1轮胎动变形加权系数q2悬架动扰度加权系数式(3.9) 中第一、 二项为误差指标， 表示在 0 整个时间内系统实际状态与平衡之间的误差总和。这一积分越小，说明控制误差越小，性能越好。积分式中第三项为能量指标，表示在0 整个时间内支付能量的总和。系统状态转移是考控制u(t)来进行的，为要使系统误差很小，则需要支付很大的能量代价。最
5、优反馈增益系数矩阵式可写成。u(t )kXk1 ( x2x1 )k2 x2 k3 ( x1x0 ) k4 x1 (6)式中，增益值k1k 4 有明确的物理意义。k1 可等效于一放置于簧载和非簧载质量间的弹簧，改变 k1 则影响簧载质量的固有频率；k2 作用于簧载质量的绝对速度上，影响其悬挂阻尼；k3 大小涉及轮胎变形，对车轮的垂直弹跳频率产生影响；k4 作用于非簧载质量的速度上，影响其非悬挂阻尼。二、 主动悬架系统的能控性，能观测性能控性和能观测性是系统的一种特性，是现代控制理论中的两个基本概念。状态完全能控的充分必要条件是能控矩阵C o A ABAn 1 B 满秩；状态完全能观测的充分必要条
6、件是能观测矩阵Ob C T ATC T( AT ) n 1C T满秩。主动悬架系统参数的选取如下，即m1=36kg ； m2=240kg ；kt=160000N/m ；将参数值带入矩阵，利用Matlab 中的函数 Cctrb(A,B) 求悬架系统的能控矩阵C ，利oo用函数 rA =rank(C o)得矩阵的秩为rA =4,满秩，故系统是能控的。利用 Matlab 中的函数 Ob obsv(A,C) 求悬架系统的能观测矩阵Ob，利用函数 rB rank(O b)得矩阵的秩为 rB=4, 满秩，故系统是能观测的。三、 主动悬架的频域仿真为了求得主动悬架系统的最优控制u(t), 必须先求得反馈增益
7、矩阵K ，而 K 矩阵的求解决定于黎卡提代数方程的解P 矩阵，这可以用计算机来实现求解。程序用 Matlab 语言编写，给定一组矩阵A 、 B 、Q、 R 的有关数据，经过计算，便可以最终得到相应的矩阵P、 K 的数值。下面取三组不同权系数q1， q2 进行计算分析；1) 取 q1=3.35E5,q 2=40.5E5 时，由程序得 k1 =2012.5,k 2=977.1,k 3= 1874.8,k4= 31.3，并求得系统的传递函数及幅频特性，绘制系统的幅频特性图%主动悬架q1=3.35e5;q2=40.5e5 时的仿真程序：m1=36;m2=240;kt=160000;q1=3.35e5;
8、q2=40.5e5;A=0 1 0 -1;0 0 0 0;0 0 0 -1;0 0 kt/m1 0;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=q2 0 0 0;0 0 0 0;0 0 q1 0;0 0 0 0;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;for s=0:0.1:80s=s*2*pi*j;G11=(7.811 *s3 + 580.4 *s2 + 3.727e004 *s + 1.422e-010)/(s4 +
9、4.942 *s3 + 4457 *s2 + 1.809e004 *s + 3.727e004);G12=(-4385 *s - 1.751e004)/(s4 + 4.942 *s3 + 4457 *s2 + 1.809e004 *s + 3.727e004); G13=(s3 + 4.942 *s2 + 64.29 *s - 2.145e-013)/(s4 + 4.942 *s3 + 4457 *s2 + 1.809e004 *s + 3.727e004);f(i)=abs(G11);h(i)=abs(G12);g(i)=abs(G13);i=i+1;ends=0:0.1:80;figure
10、loglog(s,f,-,s,h,-.,s,g,:)legend(加速度 ,动扰度 ,动载荷 )图 1. q1 3.35E5， q2 40.5E5 的幅频特性图由图 1 可以看出主动悬架的车身加速度、悬架动扰度、 轮胎动载荷幅频特性图同被动悬架相似，同样具有双峰，不同的是在低频固有频率附近，主动悬架的响应幅值明显减小，且变化平缓，主动悬架的减振性能较为突出；在高频固有频率附近，主动悬架的响应幅值变化较大。可知取该组权系数时，主动悬架的减振性能的改善程度不够理想；2) 取 q1 3.35E8，q2 40.5E8 时，由程序得 k1=63640;k 2=4863;k 3= 36146；k4 904
11、；及系统的传递函数和幅频特性，绘制幅频特性图%主动悬架q1=3.35e8;q2=40.5e8 时的仿真程序:m1=36;m2=240;kt=160000;q1=3.35e8;q2=40.5e8;A=0 1 0 -1;0 0 0 0;0 0 0 -1;0 0 kt/m1 0;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=q2 0 0 0;0 0 0 0;0 0 q1 0;0 0 0 0;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;for s=0:0.1:80s=s*2*pi*j;G11=(150.6*s3 + 1.673e004*s2 + 1.179e006*s + 1.653e-008)/(s4 + 45.36*s3 + 5473*s2 + 9.005e004*s + 1.179e006);G12=(-3290*s - 7.332e004
