基于某Matlab的汽车主动悬架控制器设计与仿真.doc
word目录1、研究背景32、仿真系统模型的建立42.1被动悬架模型的建立42.2主动悬架模型的建立53、LQG控制器设计64、仿真输出与分析7795、总结10附录：MATLAB程序源代码11一主动悬架车辆模型11二被动悬架车辆模型12三均方根函数131、研究背景汽车悬架系统由弹性元件、导向元件和减振器组成,是车身与车轴之间连接的所有组合体零件的总称,也是车架(或承载式车身)与车桥(或车轮)之间一切力传递装置的总称,其主要功能是使车轮与地面有很好的附着性,使车轮动载变化较小,以保证车辆有良好的安全性,缓和路面不平的冲击,使汽车行驶平顺,乘坐舒适,在车轮跳动时,使车轮定位参数变化较小,保证车辆具有良好的操纵稳定性。(a)被动悬架系统(b)半主动悬架系统(c)主动悬架系统图1 悬架系统汽车的悬架种类从控制力学的角度大致可以分为被动悬架、半主动悬架、主动悬架3种如图1所示。目前,大局部汽车使用被动悬架 ,这种悬架在路面不平或汽车转弯时,都会受到冲击,从而引起变形,这时弹簧起到了减缓冲击的作用,同时弹簧释放能量时,产生振动。为了衰减这种振动 ,在悬架上采用了减振器 ,这种悬架作用是外力引起的,所以称为被动悬架。半主动悬架由可控的阻尼与弹性元件组成 ,悬架的参数在一定围可以任意调节。主动悬架是在控制环节中安装了能够产生上下移动力的装置,执行元件针对外力的作用产生一个力来主动控制车身的移动和车轮受到的载荷,即路面的反作用力。随着电控技术的开展，微处理器在车辆中的应用已经日趋普遍,再加上作动器、可调减振器和变刚度弹簧等重大技术的突破,使人们更加注对主动悬架系统的研究。车辆悬架的特性可以从车身垂直加速度,悬架动行程以与轮胎动位移来研究。本文对主动悬架采用LQG最优设计谋略,利用MATLAB/Simulink软件进展仿真,分别对被动悬架与主动悬架建立动力学模型,并对两种悬架的仿真结果做了详细的比拟分析与说明。2、仿真系统模型的建立根据牛顿运动定律，利用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为被动悬架单轮车辆模型，其微分方程为：12图2 1/4车辆被动悬架模型引用路面输入模型为：3式中：f0为下截止频率Hz；G0为路面不平度系数m3/cycle，v0为前进车速m/s；w为数字期望为零的高斯白噪声。选取状态变量为：，结合式1、式2、式3，将系统运动方程与路面激励写成矩阵形式，得出系统空间状态方程：4式中，A为状态矩阵；F为输入矩阵；W=wt，为高斯白噪声输入矩阵。其值如下：；。将车身加速度、悬架动行程、轮胎动位移作为性能指标：Y=T。如此可将性能指标写成状态变量与输入信号的线性组合形式：5式中，C为输入矩阵：。同理，运用牛顿运动定律，利用1/4车辆模型特性，建立一个具有主动悬架1/4车辆的动力学模型。其中Us为作动力控制力。图3为主动悬架单轮车辆模型。图3 1/4车辆主动悬架模型其微分方程为：67整理为状态方程：89式中：U为作动器控制矩阵；B，D为新增输入矩阵。其值为：;。3、LQG控制器设计车辆悬架设计中的主要性能指标包括：代表轮胎接地性的轮胎动载荷；代表乘坐舒适性的车身垂向振动加速度；影响车身姿态且与结构设计和布置有关的悬架动行程。因此，LQG控制器设计中的性能指标J即为轮胎动位移、悬架动行程和车身垂向振动加速度的加权平方和在时域T的积分值，其表达式为：10式中，q1、q2和q3分别为轮胎动位移、悬架动行程和车身垂向振动加速度的加权系数。将性能指标J的表达式8改写成矩阵形式：11式中，Q对应于状态变量的权重矩阵；R为约束输入信号大小的权重矩阵；N为耦合项。；。当车辆参数值和加权系数值确定后，最优控制反应增益矩阵可由黎卡提Riccati方程求出，其形式如下：12最优控制反应增益矩阵，由车辆参数和加权系数决定。根据任意时刻的反应变量X(t),就可得出t时刻作动器的最优控制力，即：134、仿真输出与分析选择某轿车的后悬架作为相关计算参数：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，主动悬架Ks=20 kN/m，Cs=0 N·s/m，被动悬架Ks=22 kN/m，Cs=1 kN·s/m。仿真计算中以式3所示的滤波白噪声作为路面输入模型。白噪声的生成可直接调用MATLAB函数WGNM，N，P，其中M为生成矩阵的行数，N为列数，P为白噪声的功率单位dB。仿真计算中取一条白噪声，共10001个采样点，噪声强度为20dBM=10001，N=1，P=20。设定采样时间为0.005s、车速为20m/s时，仿真路面长度为1km，仿真时间为50s。根据所建立的系统状态方程式4、8与最优性能指标函数式11，利用的矩阵A、B、Q、R、N，调用MATLAB中的线性二次最优控制设计函数K，S，E=LQRA，B，Q，R，N，即可完成最优主动悬架控制器的设计。输出的结果中，K为最优控制反应增益矩阵，S为黎卡提方程的解，E为系统闭环特征根。带入仿真输入参数，可求得最优反应增益矩阵K为：K=711.88-1241.5-19284-2038.520864黎卡提方程的解为：。在Simulink环境下建立的最优主动悬架车辆仿真模型框图如图4所示。LQG主动悬架系统和被动悬架系统的时域仿真结果分别如如图5、图6所示，包括路面位移输入xgt、悬架动行程SWSt、轮胎动位移DTDt与车身加速度BAt。图4 Simulink环境系统仿真框图BA/(m/s2)SWS/mDTD/mxg/m时间/s图5 在某路面输入下的最优主动悬架仿真结果BA/(m/s2)SWS/mDTD/mxg/m时间/s图6 在某路面输入下的最优被动悬架仿真结果计算得出的两个不同系统的性能指标均方根值如表1所示。表1 主动悬架与被动悬架性能指标均方根值的比拟性能指标单位主动悬架均方根值被动悬架均方根值车身加速度 BAm/s2悬架动行程 SWS轮胎动位移 DTD比拟图5、图6可得，车身垂直加速度在12s，主动悬架车辆和被动悬架车辆的垂直加速度相差不大；随着时间的变化，被动悬架的车身垂直加速度比主动悬架的加速度越来越大，垂直加速度越大，车辆的舒适性越差。主动悬架的轮胎动位移始终控制在-22mm之间，而被动悬架的动位移很大，动位移越大，车辆的舒适性越差。悬架动行程在01s，车辆悬架的动行程相差不大，但是随着时间的推移，主动悬架的车辆的动行程波动较小，比被动悬架车辆的舒适性好。同时由表1可见，在轮胎动位移根本一样的情况下，所设计的最优主动悬架显著的降低了车身的垂向振动加速度，与被动悬架系统相比，其均方根值减少了18%。主动悬架系统的仿真结果还明确，其悬架动行程同时也被很好的控制在设计要求的围±100mm，意味着许用的悬架工作空间得到了更充分的利用。5、总结主动悬架与被动悬架的主要区别是有无阻尼，本文结合现代控制理论，通过建立车辆的主动悬架与被动悬架的数学模型，利用MATLAB/Simulink软件进展仿真输出，对车身垂直加速度、悬架位移与轮胎变形进展比拟，可以明显看出主动悬架要优于被动悬架的车辆，说明主动悬架车辆比被动悬架车辆舒适性好，与实际情况一致。通过教师的指导和同学们的帮助，对学习过程中遇到的问题逐一解决，将MATLAB与现代控制理论合理的结合一起，提高了问题的解决效率，加深了对现代控制理论的理解和对MATLAB软件的熟悉程度。本文以1/4车辆模型为研究容，相对整车模型仍有不足，在今后的学习过程中有待进一步研究。附录：MATLAB程序源代码一主动悬架车辆模型clcclear all%主动悬架车辆模型参数mb=320; %簧载质量kgmw=40; %非簧载质量kgks=20000; %主动悬架刚度N/mkt=200000; %轮胎刚度N/mswsc=+-100; %悬架工作空间mm%仿真路面输入参数g0=5*10-6; %路面不平度系数m3/cycleu=20; %车速m/sf0=0.1; %下截止频率Hz%性能指标加权系数q1=80000; %轮胎动位移q2=5; %悬架动行程q3=1; %车身加速度%系统输入x0=wgn(10001,1,20); %生成高斯白噪声采样点数 列数 功率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)/(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; %簧载质量kgmw=40; %非簧载质量kgks=22000; %被动悬架刚度N/mcs=1000; %阻尼系数N.s/mkt=200000; %轮胎刚度N/mswsc=+-100; %悬架工作空间mm%仿真路面输入参数g0=5*10-6; %路面不平度系数m3/cycleu=20; %车速m/sf0=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-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为最优控制反应增益矩阵三均方根函数%性能指标均方根fBA=sqrt(sum(BA(:,2).2)/6136) %车身加速度均方根fSWS=sqrt(sum(SWS(:,2).2)/6136) %悬架动行程均方根fDTD=sqrt(sum(DTD(:,2).2)/6136) %轮胎动位移均方根12 / 12
