基于Matlab的汽车主动悬架控制器设计与仿真

1、现代控制理论及其应用课程小论文基于Matlab的汽车主动悬架控制器设计与仿真学院：机 械 工 程 学 院班级：（XX）姓名： X X X 2015年6月3号河 北 工 业 大 学目录1、研究背景32、仿真系统模型的建立42.1被动悬架模型的建立42.2主动悬架模型的建立53、LQG控制器设计64、仿真输出与分析74.1仿真的输出74.2仿真结果分析95、总结10附录：MATLAB程序源代码11（一）主动悬架车辆模型11（二）被动悬架车辆模型12（三）均方根函数131、研究背景汽车悬架系统由弹性元件、导向元件和减振器组成,是车身与车轴之间连接的所有组合体零件的总称,也是车架(或承载式车身)与车桥

2、(或车轮)之间一切力传递装置的总称,其主要功能是使车轮与地面有很好的附着性,使车轮动载变化较小,以保证车辆有良好的安全性,缓和路面不平的冲击,使汽车行驶平顺,乘坐舒适,在车轮跳动时,使车轮定位参数变化较小,保证车辆具有良好的操纵稳定性。(a)被动悬架系统(b)半主动悬架系统(c)主动悬架系统图1 悬架系统汽车的悬架种类从控制力学的角度大致可以分为被动悬架、半主动悬架、主动悬架3种（如图1所示）。目前,大部分汽车使用被动悬架 ,这种悬架在路面不平或汽车转弯时,都会受到冲击,从而引起变形,这时弹簧起到了减缓冲击的作用,同时弹簧释放能量时,产生振动。为了衰减这种振动 ,在悬架上采用了减振器 ,这种悬

3、架作用是外力引起的,所以称为被动悬架。半主动悬架由可控的阻尼及弹性元件组成 ,悬架的参数在一定范围内可以任意调节。主动悬架是在控制环节中安装了能够产生上下移动力的装置,执行元件针对外力的作用产生一个力来主动控制车身的移动和车轮受到的载荷,即路面的反作用力。随着电控技术的发展，微处理器在车辆中的应用已经日趋普遍,再加上作动器、可调减振器和变刚度弹簧等重大技术的突破,使人们更加注对主动悬架系统的研究。车辆悬架的特性可以从车身垂直加速度,悬架动行程以及轮胎动位移来研究。本文对主动悬架采用LQG最优设计策略,利用MATLAB/Simulink软件进行仿真,分别对被动悬架与主动悬架建立动力学模型,并对两

4、种悬架的仿真结果做了详细的比较分析与说明。2、仿真系统模型的建立2.1被动悬架模型的建立根据牛顿运动定律，利用1/4车辆模型特性，建立被动悬架1/4车辆的动力学模型。其中mb代表车身质量（kg），mw代表车轮质量（kg），xb代表车身位移（m），xw代表车轮位移（m），Ks代表悬架弹簧刚度（N/m），Kt代表轮胎刚度（N/m），Cs代表悬架阻尼（N·s/m），xg代表路面位移（m），代表车身加速度（m/s2），代表车轮速度（m/s），代表车轮加速度（m/s）。图2为被动悬架单轮车辆模型，其微分方程为：（1）（2）图2 1/4车辆被动悬架模型引用路面输入模型为：（3）式中：f0为下截止

5、频率（Hz）；G0为路面不平度系数（m3/cycle），v0为前进车速（m/s）；w为数字期望为零的高斯白噪声。选取状态变量为：，结合式（1）、式（2）、式（3），将系统运动方程及路面激励写成矩阵形式，得出系统空间状态方程：（4）式中，A为状态矩阵；F为输入矩阵；W=（w（t），为高斯白噪声输入矩阵。其值如下：；。将车身加速度、悬架动行程、轮胎动位移作为性能指标：Y= T。则可将性能指标写成状态变量及输入信号的线性组合形式：（5）式中，C为输入矩阵：。2.2主动悬架模型的建立同理，运用牛顿运动定律，利用1/4车辆模型特性，建立一个具有主动悬架1/4车辆的动力学模型。其中Us为作动力控制力。图3

6、为主动悬架单轮车辆模型。图3 1/4车辆主动悬架模型其微分方程为：（6）（7）整理为状态方程：（8）（9）式中：U为作动器控制内矩阵；B，D为新增输入矩阵。其值为：;。3、LQG控制器设计车辆悬架设计中的主要性能指标包括：代表轮胎接地性的轮胎动载荷；代表乘坐舒适性的车身垂向振动加速度；影响车身姿态且与结构设计和布置有关的悬架动行程。因此，LQG控制器设计中的性能指标J即为轮胎动位移、悬架动行程和车身垂向振动加速度的加权平方和在时域T内的积分值，其表达式为：（10）式中，q1、q2和q3分别为轮胎动位移、悬架动行程和车身垂向振动加速度的加权系数。将性能指标J的表达式（8）改写成矩阵形式：（11）

7、式中，Q对应于状态变量的权重矩阵；R为约束输入信号大小的权重矩阵；N为耦合项。；。当车辆参数值和加权系数值确定后，最优控制反馈增益矩阵可由黎卡提（Riccati）方程求出，其形式如下：（12）最优控制反馈增益矩阵，由车辆参数和加权系数决定。根据任意时刻的反馈变量X(t),就可得出t时刻作动器的最优控制力，即：（13）4、仿真输出与分析4.1仿真的输出选择某轿车的后悬架作为相关计算参数：mb=320 kg，mw=40 kg，Kt=200 kN/m，悬架工作空间SWSC=±100 mm，G0=5 cm3/cycle，u=20 m/s，f0=0.1Hz，q1=80000，q2=5，q3=1

8、，主动悬架Ks=20 kN/m，Cs=0 N·s/m，被动悬架Ks=22 kN/m，Cs=1 kN·s/m。仿真计算中以式（3）所示的滤波白噪声作为路面输入模型。白噪声的生成可直接调用MATLAB函数WGN（M，N，P），其中M为生成矩阵的行数，N为列数，P为白噪声的功率（单位dB）。仿真计算中取一条白噪声，共10001个采样点，噪声强度为20dB（M=10001，N=1，P=20）。设定采样时间为0.005s、车速为20m/s时，仿真路面长度为1km，仿真时间为50s。根据所建立的系统状态方程式（4）、（8）及最优性能指标函数式（11），利用已知的矩阵A、B、Q、R、N，

9、调用MATLAB中的线性二次最优控制设计函数K，S，E=LQR（A，B，Q，R，N），即可完成最优主动悬架控制器的设计。输出的结果中，K为最优控制反馈增益矩阵，S为黎卡提方程的解，E为系统闭环特征根。带入仿真输入参数，可求得最优反馈增益矩阵K为：K=（711.88-1241.5-19284-2038.520864）黎卡提方程的解为：。在Simulink环境下建立的最优主动悬架车辆仿真模型框图如图4所示。LQG主动悬架系统和被动悬架系统的时域仿真结果分别如如图5、图6所示，包括路面位移输入xg（t）、悬架动行程SWS（t）、轮胎动位移DTD（t）及车身加速度BA（t）。图4 Simulink环境

10、系统仿真框图BA/(m/s2)SWS/mDTD/mxg/m时间/s图5 在某路面输入下的最优主动悬架仿真结果BA/(m/s2)SWS/mDTD/mxg/m时间/s图6 在某路面输入下的最优被动悬架仿真结果计算得出的两个不同系统的性能指标均方根值如表1所示。表1 主动悬架与被动悬架性能指标均方根值的比较性能指标单位主动悬架均方根值被动悬架均方根值车身加速度 BAm/s21.41281.7828悬架动行程 SWSmm32.317.7轮胎动位移 DTDmm5.86.04.2仿真结果分析比较图5、图6可得，车身垂直加速度在12s内，主动悬架车辆和被动悬架车辆的垂直加速度相差不大；随着时间的变化，被动悬

11、架的车身垂直加速度比主动悬架的加速度越来越大，垂直加速度越大，车辆的舒适性越差。主动悬架的轮胎动位移始终控制在-22mm之间，而被动悬架的动位移很大，动位移越大，车辆的舒适性越差。悬架动行程在01s内，车辆悬架的动行程相差不大，但是随着时间的推移，主动悬架的车辆的动行程波动较小，比被动悬架车辆的舒适性好。同时由表1可见，在轮胎动位移基本相同的情况下，所设计的最优主动悬架显著的降低了车身的垂向振动加速度，与被动悬架系统相比，其均方根值减少了18%。主动悬架系统的仿真结果还表明，其悬架动行程同时也被很好的控制在设计要求的范围内（±100mm），意味着许用的悬架工作空间得到了更充分的利用。

12、5、总结主动悬架与被动悬架的主要区别是有无阻尼，本文结合现代控制理论，通过建立车辆的主动悬架与被动悬架的数学模型，利用MATLAB/Simulink软件进行仿真输出，对车身垂直加速度、悬架位移及轮胎变形进行比较，可以明显看出主动悬架要优于被动悬架的车辆，说明主动悬架车辆比被动悬架车辆舒适性好，与实际情况一致。通过老师的指导和同学们的帮助，对学习过程中遇到的问题逐一解决，将MATLAB与现代控制理论合理的结合一起，提高了问题的解决效率，加深了对现代控制理论的理解和对MATLAB软件的熟悉程度。本文以1/4车辆模型为研究内容，相对整车模型仍有不足，在今后的学习过程中有待进一步研究。附录：MATLA

13、B程序源代码（一）主动悬架车辆模型clcclear all%主动悬架车辆模型参数mb=320; %簧载质量（kg）mw=40; %非簧载质量（kg）ks=20000; %主动悬架刚度（N/m）kt=200000; %轮胎刚度（N/m）swsc=+-100; %悬架工作空间（mm）%仿真路面输入参数g0=5*10-6; %路面不平度系数（m3/cycle）u=20; %车速（m/s）f0=0.1; %下截止频率（Hz）%性能指标加权系数q1=80000; %轮胎动位移q2=5; %悬架动行程q3=1; %车身加速度%系统输入x0=wgn(10001,1,20); %生成高斯白噪声（采样点数 列数

14、 功率）t=0:0.005:50; %仿真时间Tw=t' x0; %系统输入序列仿真时间序列 白噪声序列%系统状态方程矩阵A=0 0 -ks/mb ks/mb 0;0 0 ks/mw -(kt+ks)/mw kt/mw;1 0 0 0 0;0 1 0 0 0;0 0 0 0 -2*pi*f0B=1/mb;-1/mw;0;0;0F=0;0;0;0;2*pi*(g0*u)0.5C=1 0 0 0 0;0 1 0 0 0;0 0 1 0 0;0 0 0 1 0;0 0 0 0 1D=1/mb 0 0 0%二次最优控制器设计Q=0 0 0 0 0;0 0 0 0 0;0 0 q2+(ks2)/

15、(mb2) -q2-(ks2)/(mb2) 0;0 0 -q2-(ks2)/(mb2) q1+q2+(ks2)/(mb2) -q1;0 0 0 -q1 q1;R=1/(mb2);N=0;0;-ks/(mb2);ks/(mb2);0;K,S,E=LQR(A,B,Q,R,N); %K为最优控制反馈增益矩阵（二）被动悬架车辆模型clcclear all%被动悬架车辆模型参数mb=320; %簧载质量（kg）mw=40; %非簧载质量（kg）ks=22000; %被动悬架刚度（N/m）cs=1000; %阻尼系数（N.s/m）kt=200000; %轮胎刚度（N/m）swsc=+-100; %悬架工作

16、空间（mm）%仿真路面输入参数g0=5*10-6; %路面不平度系数（m3/cycle）u=20; %车速（m/s）f0=0.1; %下截止频率（Hz）%性能指标加权系数q1=80000; %轮胎动位移q2=5; %悬架动行程q3=1; %车身加速度%系统输入x0=wgn(10001,1,20); %生成高斯白噪声（采样点数 列数 功率）t=0:0.005:50; %仿真时间Tw=t' x0; %系统输入序列仿真时间序列 白噪声序列%系统状态方程矩阵A=-cs/mb cs/mb -ks/mb ks/mb 0;cs/mw -cs/mw ks/mw -(kt+ks)/mw kt/mw;1 0 0 0 0;0 1 0 0 0;0 0 0 0 -2*pi*f0B=0;0;0;0;2*pi*(g0*u)0.5C=1 0 0 0 0;0 1 0 0 0;0 0 1 0 0;0 0 0 1 0;0 0 0 0 1D=0 0 0 0%二次最优控制器设计Q=0 0 0 0 0;0 0 0 0 0;0 0 q2+ks2/mb2 -q2