行驶动力学建模、仿真与主动悬架控制器设计说明.doc

《行驶动力学建模、仿真与主动悬架控制器设计说明.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、......目录1.计算机仿真系统模型的建立22.LOG控制器设计33.计算实例44.MATLAB仿真过程55.半车模型建模及仿真95.1随机线性最优控制105.2预瞄控制125.3结果比较13专业.专注.......以单轮车辆模型为例，介绍行驶动力学计算机建模、仿真分析以及利用线性二次最优控制理论进行主动悬架LQG控制器设计过程。1.计算机仿真系统模型的建立根据图7所示的主动悬架单轮车辆模型，运用牛顿运动定律，建立系统的运动方程，即：（4）（5）这里，采用一个滤波白噪声作为路面输入模型，即：（6）式中，xg为路面垂向位移（m）；

2、G0为路面不平度系数（m3/cycle）；u为车辆前进速度（m/s）；w为数字期望为零的高斯白噪声；f0为下截止频率（Hz）。图7单轮车辆模型结合式（4）、式（5）和式（6），将系统运动方程和路面输入方程写成矩阵形式，即得出系统的空间状态方程：（7）专业.专注.......式中，，为系统状态矢量；W=（w（t）），为高斯白噪声输入矩阵；U=（Ua（t）），为输入控制矩阵；；；1.LOG控制器设计车辆悬架设计中的主要指标包括：①代表轮胎接地性的轮胎动载荷；②代表轮胎舒适性的车身垂向振动加速度；③影响车身姿态且与轮胎布置有关的悬架动行

3、程。因此，LQG控制器设计中的性能指标J即为轮胎动位移、悬架动行程和车身垂向振动加速度的加权平方和在时域T内的积分值，其表达式为：（8）式中，q1、q2和q3分别为轮胎动位移、悬架动行程和车身垂向振动加速度的加权系数。加权系数的选取决定了设计者对悬架性能的倾向，如对车身垂向振动加速度项选择较大的权值，则考虑更多的是提高车辆操纵稳定性。为方便起见，这里取车身垂向振动加速度的加权值q3=1。将性能指标J的表达式（8）改写成矩阵形式，即：（9）式中，专业.专注.......；；当车辆参数值和加权系数值确定后，最优控制反馈增益矩阵可有黎卡

4、提（Riccati）方程求出，其形式如下：（10）最优反馈控制增益矩阵，由车辆参数和加权系数决定。根据任意时刻的反馈状态变量X（t），就可得到t时刻作动器的最优控制力Ua，即：（11）1.计算实例这里，以某轿车的后悬架为例，给出一个完整的计算实例，包括车辆模型参数、仿真路面输入参数、控制器的设计参数以及计算结果。此例中车辆以20m/s的速度在某典型路面上行驶，仿真时间T=50s。计算中输入的各参数及数值详见表2。表2单轮车辆模型仿真输入参数值车辆模型参数符号单位数值簧载质量非簧载质量悬架刚度轮胎刚度悬架工作空间mbmwKsKtSW

5、ScKgKgN/mN/mmm3204020000200000仿真路面输入参数符号单位数值专业.专注.......路面不平度系数车速下截止频率G0Uf0m3/cyclem/sHz5.0x10-6200.1性能指标加权参数符号单位数值轮胎动位移悬架动行程车身加速度q1q2q3———8000051仿真计算中以式（6）所示的滤波白噪声作为路面输入模型。白噪声的生成可直接调用MATLAB函数WGN（M，N，P）（此函数需要安装信号处理工具箱Communicationstoolbox），其中M为生成矩阵的行数，N为列数，P为白噪声的功率（单位

6、为dB）。本例中取M=10001，N=1，P=20。这意味着仿真计算中去一条白噪声，共10001个采集点，噪声强度为20dB。设定采样时间为0.005s、车速为20m/s时，相当于仿真路面长度为1000m，仿真时间为50s。根据建立的系统状态方程式（7）及最优化性能指标函数式（9），利用已知的矩阵A、B、Q、R、N，调用MATLAB中的线性二次最优控制器设计函数[K，S，E]=LQR（A，B，Q，R，N），即可完成最优主动悬架控制器的设计。输出的结果中，K为最优控制反馈增益矩阵，S为黎卡提方程的解，E为系统闭环特征根。根据表2给出

7、的仿真输入参数，本例中求得的最优反馈增益矩阵K为：K=（711.88-1241.5-19284-2038.520864）同时，还得到了黎卡提方程的解：专业.专注.......在相同的仿真条件下，可将所设计的主动悬架系统与一个被动系统进行对比分析。在被动悬架系统中，取悬架刚度Ks=22000N/m，阻尼系数Cs=1000NS/m。除此之外，其他输入参数值均与主动悬架系统完全相同。1.MATLAB仿真过程1）生成路面输入模型代码如下：a=wgn(10001,1,20);t=0:0.005:50;road_file(:,1)=t';ro

8、ad_file(:,2)=a;saveroad_fileroad_file2）参数输入代码如下：loadroad_file.mat%载入路面数据模型Ks=22000;mb=320;Kt=200000;mw=40;f0=0.1;G0=0.000005