现代控制理论大作业分析对象：汽车悬架系统
指导老师：***
专业：机械工程
******
学号：S********
1.建模
悬架是车轮或车桥与汽车承载部分之间具有弹性的连接装置的总称，具有传递载荷、缓和冲击、衰减振动以及调节汽车行驶中的车身位置等作用。传统汽车悬驾系统是被动悬驾，其参数不能改变，无法控制其对不同路面激励的响应，因此对不同路面的适应性较差。为提高汽车的行驶平顺性、操纵稳定性和制动性等性能，人们开始用主动悬架系统来代替传统的被动悬架系统。主动悬架系统能根据路面的情况通过一个动力装置改变悬挂架的参数，改善汽车的各方面性能。
对悬驾系统进行仿真计算首先要建立悬驾系统动力学模型，随后对所建立的模型进行仿真分析。为了简化模型，取汽车的一个车轮的悬驾系统进行研究，该模型可简化为一维二自由度的弹簧阻尼质量系统，图1所示为该模型的模拟图。
图1 悬架系统模型的模拟图
其中u 为动力装置的作用力，w 为路面位移，x1为车身位移，x2为悬驾位移，用车身位移来度量车身的振动情况，并视为系统的输出。路面状况以w 为尺度，并视为系统的一个干扰输入。当汽车从平面落入坑时，w 可用一个阶跃信号来模拟。u 为主动悬架的作用力，它是系统的控制量。
进行受力分析，由牛顿第二规律可得车身悬架系统的动力学方程为:
()()()()()
1121212212122s s t m x K x x b x x u
m x K x x b x x u K w x ⎧=-+-+⎪⎨
=-+--+-⎪⎩ 设系统状态变量为:
[]1
21
2x x x x x =
则上面系统动力学方程可改写为状态空间表达式:
x Ax Bu
y Cx Du
=+⎧⎨
=+⎩ 其中:
()1
1
1
1222
200
100001s
s
s t s K K b b A m m m m K K K b b m m m m ⎡⎤⎢⎥⎢⎥⎢⎥--=⎢⎥⎢⎥-+⎢⎥-⎢⎥⎣
⎦
12
200
001
01t B m K m m ⎡⎤
⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥-⎢⎥⎣⎦
[]1000C = []00D = u u w ⎡⎤=⎢⎥⎣⎦
Matlab 系统模型程序代码： m1=800;m2=320;ks=10000;b=30000; kt=10*ks;
a=[0 0 1 0;0 0 0 1;-ks/m1 ks/m1 -b/m1 b/m1; ks/m2 -(ks+kt)/m2 b/m2 -b/m2]; b=[0 0;0 0;1/m1 0;-1/m2 kt/m2]; c=[1 0 0 0]; d=[0 0]; G=ss(a,b,c,d)
[num,den]=ss2tf(a,b,c,d,1);
tfsys=tf(num,den)
eig(a)
figure(1)
impulse(G)
figure(2)
step(G)
系统状态空间模型：:
系统传递函数：
系统矩阵特征值：
系统矩阵特征值都在左半平面，故系统稳定。系统脉冲响应：
系统阶跃响应：
2.系统能控性与能观性分析
Matlab系统能控性与能观性分析程序代码：m1=800;m2=320;ks=10000;b=30000;
kt=10*ks;
a=[0 0 1 0;0 0 0 1;-ks/m1 ks/m1 -b/m1 b/m1;
ks/m2 -(ks+kt)/m2 b/m2 -b/m2];
b=[0 0;0 0;1/m1 0;-1/m2 kt/m2];
c=[1 0 0 0];
d=[0 0];
N=size(a);
n=N(1)
ud=ctrb(a,b);%能控矩阵%
ob=obsv(a,c);%能观矩阵%
if rank(ud)==n%判断能控性%
disp('系统能控')
elseif rank(ud)disp('系统不能控')
end
if rank(ob)==n %判断能观性%
disp('系统能观')
elseif rank(ob)disp('系统不能观')
end
运行结果：
3.系统稳定性分析
Matlab系统稳定性分析程序代码：
m1=800;m2=320;ks=10000;b=30000;
kt=10*ks;
a=[0 0 1 0;0 0 0 1;-ks/m1 ks/m1 -b/m1 b/m1;
ks/m2 -(ks+kt)/m2 b/m2 -b/m2];
Q=eye(4,4);
P=lyap(a,Q);%求解矩阵P
flag=0;
n=length(a);
for i=1:n
det(P(1:i,1:i))
if(det(P(1:i,1:i))<=0)
flag=1;
end
end
if flag==1
disp('系统不稳定');
else
disp('系统稳定');
end
运行结果：
4.系统状态反馈控制器设计
若系统期望性能指标为超调量小于5%，峰值时间小于0.5秒，即：
(
)
22
exp 1505
1%
.p n t w σξπ
ξπ
ξ
=--≤=
≤-
则：
0707
9
.n w ξ≥≥
为计算方便取
0707
10
.n w ξ==
则主导极点为：
2121707707..n n w jw j λλξξ==-±-=-±
