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

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

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

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

4、真系统模型的建立2.1被动悬架模型的建立根据牛顿运动定律，利用1/4车辆模型特性，建立被动悬架1/4车辆的动力学模型。其 中mb代表车身质量（kg），mw代表车轮质量（kg），Xb代表车身位移（m），Xw代表车轮位移（m）, K.代表悬架弹簧刚度（N/m），Kt代表轮胎刚度（N/m），CS代表悬架阻尼（Ns/m），Xg代表路面位移（m）， Rb代表车身加速度（m/s?）， *w代表车轮速度（m/s）， Xw代表车轮加速度（m/s）图2为被动悬架单轮车辆模型，其微分方程为:mbK s (xb xw ) Cs (Rb *w )(1)mwRW(2)Ks(XbXw)Kt(Xg-5乍身位移Rb车身这度”

5、*车轮位移 *车轮速度轮胎刚度K.路面位移图2 1/4车辆被动悬架模型引用路面输入模型为:(3)Xg2 foXg(t)2、._GoUw(t)式中：fo为下截止频率（Hz）； G为路面不平度系数（m/cycle ），Vo为前进车速（m/s）; w为 数字期望为零的高斯白噪声。选取状态变量为：X Xb Xw Xb Xw XgT，结合式（1）、式（2）、式（3），将系统运动方程及路面激励写成矩阵形式，得出系统空间状态方程:(4)AX BU FW式中，A为状态矩阵；F为输入矩阵；W=(w(t),为高斯白噪声输入矩阵。其值如下:CsmbCsCl mbCsKsmbKsKsmbKt Ks0Kt00Amwmw

6、mwmwmw ； F0100000010002 , G0U00002 f。将车身加速度、悬架动行程、轮胎动位移作为性能指标:Y=Rb (Xb xw) (Xw Xg)T。则可将性能指标写成状态变量及输入信号的线性组合形式:(5)1000001000式中，C为输入矩阵：C 0010000010Y CX0 0 0 0 12.2主动悬架模型的建立同理，运用牛顿运动定律，利用1/4车辆模型特性，建立一个具有主动悬架1/4车辆的动力学模型。其中US为作动力控制力。图车轮悬架弹戡刚度仏乍身3为主动悬架单轮车辆模型图3 1/4车辆主动悬架模型其微分方程为:Ks(Xb Xw) Cs(*b *w)mwXW(6)U

7、aKs(Xb Xw) Kt(Xg Xw) Cs(Xb Xw)整理为状态方程：X AX BU FW(7)(8)Y CX DU(9)式中：U为作动器控制矩阵；B, D为新增输入矩阵。其值为：1 1B mbmw0 0 0 ; D 0 0 0 0 0。3、LQG空制器设计车辆悬架设计中的主要性能指标包括：代表轮胎接地性的轮胎动载荷；代表乘坐舒适性的车身垂向振动加速度；影响车身姿态且与结构设计和布置有关的悬架动行程。因此，LQG空制器设计中的性能指标J即为轮胎动位移、悬架动行程和车身垂向振动加速度的加权平方和在时域T的积分值，其表达式为:TlimT T 0qXw(t)2 2 2Xg(t)q2Xb(t)

8、Xw(t)qsXb(t) dt(10)式中，q1、q2和q3分别为轮胎动位移、悬架动行程和车身垂向振动加速度的加权系数将性能指标J的表达式(8)改写成矩阵形式:1 T TTTJ lim X QX U RU 2X NU dt0式中，Q对应于状态变量的权重矩阵；R为约束输入信号大小的权重矩阵;(11)N为耦合项。000000Q00q2K：2 mb00q2K： m：0 0 00000q2Ks2mb0q1 q2Ks2 m；q1qq112mbKs。Ks0R A ； Nmb00当车辆参数值和加权系数值确定后， 最优控制反馈增益矩阵可由黎卡提 (Riccati )方程求出, 其形式如下：PA ATP (PB

9、 N)R 1(BTP NT ) Q 0(12)最优控制反馈增益矩阵 K BTP NT ，由车辆参数和加权系数决定。根据任意时刻的反馈变 量 X(t), 就可得出 t 时刻作动器的最优控制力 Ua ，即：Ua(t)KX (t)(13)4、仿真输出与分析4.1 仿真的输出选择某轿车的后悬架作为相关计算参数： mb=320 kg，mw=40 kg， Kt=200 kN/m ，悬架工 作空间 SWSC=100 mm G0=5cm3/cycle , u=20 m/s, fO=O.1Hz , q仁80000, q2=5, q3=1，主 动悬架 Ks=20 kN/m, Cs=0 N - s/m，被动悬架 K

10、s=22 kN/m, Cs=1 kN s/m。仿真计算中以式( 3)所示的滤波白噪声作为路面输入模型。白噪声的生成可直接调用MATLA函数WGNM 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、QR、N,调用MATLA中的线性二次最优控制设计函数K , S , E=LQR(

11、 A,B,QR,N) 即可完成最优主动悬架控制器的设计。输出的结果中，K为最优控制反馈增益矩阵，S为黎卡提方程的解，E为系统闭环特征根。带入仿真输入参数可求得最优反馈增益矩阵K为：K=(711.88-1241.5-19284-2038.520864)黎卡提方程的解为：2.45590.02892.47458.66077.3090.02890.48860.02987.52627.2364S2.47450.02984.97448.67545.1033 。8.66077.52628.67542710.12700.47.3097.23645.10332700.42693.7在Simulink环境下建立的

12、最优主动悬架车辆仿真模型框图如图 4所示。LQG主动悬架系 统和被动悬架系统的时域仿真结果分别如如图 5、图 6所示 包括路面位移输入 xg(t )、悬架 动行程SW( t )、轮胎动位移DTD( t)及车身加速度BA( t )。road -diEplscement inpul图4 Simulink环境系统仿真框图5C5051015202530354045 5QJV2T mABn.020 0100.01510152D253035404550e/hdII02nm0.2510Z54050丄35时间/s图5在某路面输入下的最优主动悬架仿真结果m)r mABmo型50454035305220O4O2D

13、 贬04 a a 00.时间/s图6在某路面输入下的最优被动悬架仿真结果计算得出的两个不同系统的性能指标均方根值如表1所示。表1主动悬架与被动悬架性能指标均方根值的比较性能指标单位主动悬架均方根值被动悬架均方根值车身加速度BAm/s21.41281.7828悬架动行程SWSmm32.317.7轮胎动位移DTDmm5.86.04.2仿真结果分析比较图5、图6可得，车身垂直加速度在12s,主动悬架车辆和被动悬架车辆的垂直加 速度相差不大；随着时间的变化，被动悬架的车身垂直加速度比主动悬架的加速度越来越大，垂直加速度越大，车辆的舒适性越差。主动悬架的轮胎动位移始终控制在-22mm之间，而被动悬架的动

14、位移很大，动位移越大，车辆的舒适性越差。悬架动行程在01s，车辆悬架的动行程相差不大，但是随着时间的推移，主动悬架的车辆的动行程波动较小，比被动悬架车辆 的舒适性好。同时由表 1 可见，在轮胎动位移基本相同的情况下，所设计的最优主动悬架显著的降低 了车身的垂向振动加速度， 与被动悬架系统相比， 其均方根值减少了 18%。主动悬架系统的仿 真结果还表明，其悬架动行程同时也被很好的控制在设计要求的围（土100mm,意味着许用的悬架工作空间得到了更充分的利用。5、总结主动悬架与被动悬架的主要区别是有无阻尼，本文结合现代控制理论，通过建立车辆的 主动悬架与被动悬架的数学模型，利用 MATLAB/Sim

15、ulink 软件进行仿真输出，对车身垂直加 速度、悬架位移及轮胎变形进行比较，可以明显看出主动悬架要优于被动悬架的车辆，说明 主动悬架车辆比被动悬架车辆舒适性好，与实际情况一致。通过老师的指导和同学们的帮助，对学习过程中遇到的问题逐一解决，将MATLAB!现代控制理论合理的结合一起， 提高了问题的解决效率， 加深了对现代控制理论的理解和对 MATLAB 软件的熟悉程度。本文以 1/4 车辆模型为研究容，相对整车模型仍有不足，在今后的学习过 程中有待进一步研究。附录：MATLA程序源代码（一）主动悬架车辆模型clcclear all%主动悬架车辆模型参数 mb=320; % 簧载质量（ kg）

16、mw=40; % 非簧载质量（ kg） ks=20000; % 主动悬架刚度（ N/m） kt=200000; % 轮胎刚度（ N/m） swsc=+-100; % 悬架工作空间（ mm） %仿真路面输入参数g0=5*10A-6; % 路面不平度系数（mT/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;% 系统输入序列 仿

17、真时间序列 白噪声序列%系统状态方程矩阵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=O;O;O;O;2*pi*(gO*u)95C=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 O O O%二次最优控制器设计Q=O O O O O;O O O O O;O O q2+(ksA2)/(mbA2) -q2-(ksA2)/(mbA2)O;O O -q2-(ksA2)/(mbA

18、2)q1+q2+(ksA2)/(mbA2) -q1;O O O -q1 q1;R=1/(mbA2);N=O;O;-ks/(mbA2);ks/(mbA2);O;K,S,E=LQR(A,B,Q,R,N); %K 为最优控制反馈增益矩阵(二) 被动悬架车辆模型clc clear all%被动悬架车辆模型参数mb=32O;%簧载质量( kg)mw=4O; % 非簧载质量( kg) ks=22OOO; % 被动悬架刚度( N/m) cs=1OOO; % 阻尼系数( N.s/m) kt=2OOOOO; % 轮胎刚度( N/m) swsc=+-1OO; % 悬架工作空间( mm)%仿真路面输入参数gO=5*1OA-6;% 路面不平度系数( mA3/cycle )u=2O;%车速( m/s)fO=O.1;%下截止频率( Hz)%性能指标加权系数q1=8OOOO; % 轮胎动位移q2=5;%悬架动行程q3=1;%车身加速度%系统输