汽车系统动力学上机指导书

汽车系统动力学上机指导书一、上机目的通过实验使学生了解汽车系统动力学的基本原理和研究方法，掌握系统仿真语言MATLAB的基本编程技术及其SIMULINK工具箱的应用。二、上机要求1熟悉MATLAB软件以及SIMULINK软件包的基本功能；2建立汽车子系统的动力学模型，上机仿真计算，并对结果进行分析；3对上述系统，应用汽车系统动力学理论进行仿真分析，比较仿真结果并得出结论；三、上机环境PC机、MATLAB软件四、实验内容（一）制动系统建模、仿真及ABS控制器设计通过两自由度单轮模型为例，介绍在MATLAB环境下的控制系统力学建模、ABS控制器设计及仿真分析过程。1动力学建模某车辆简化后的制动力学模型如右图所示。其中单轮模型质量M，车轮滚动半径RD，车轮转动惯量为IW，车辆旋转角速度为，车轮轮心前进速度为UW，地面制动力为FXB。作用于车轮的制动力矩为TB。若忽略空气阻力与车轮滚动阻力，则系统的运动方程如下（1）WXBDTIRT（2）XBTUMF式中，地面制动力FXB等于地面作用于车轮的法向反力FX与路面附着系数的乘积，其中为制动滑移率SB的函数。2分段线性的轮胎模型根据第三章中介绍的有关轮胎纵向特性的内容，路面附着系数与车轮滑移率之间存在一定的非线性关系。如果用两段直线近似表示路面附着滑移曲线，可得到分段线性化的附着系数与车辆滑移率SB的关系（即所谓的“DUGOFF”模型），如下图所示。其表达式如下（3）1HBBOOGHGOSS式中，B为峰值附着系数；G为车轮完全抱死时（即SB1）时的路面附着系数；SO为峰值附着系数所对应的滑移率。图1单轮制动力模型图2线性化的路面附着系数与车轮滑移率关系曲线3控制算法这里以门限值控制算法为例，说明ABS控制器设计及制动系统力学的仿真过程。采用门限制控制算法的基本思想是保证车轮滑移率在最理想的范围之内。制动开始后，随着制动压力的升高车轮转速相应减小，车轮出现滑移；当车轮滑移率达到理想范围上限值SMAX时，减小制动压力；随着制动压力的减小，直至减小到滑移率下限值SMIN时再增大制动压力。循环往复这一过程直至车辆停止。因此，在ABS控制器起作用的过程中，滑移率总是保持在理想的范围内，从而保证车辆的最佳制动性能及行驶方向控制的稳定性。4仿真流程及参数输入由上可知，ABS控制器所用到的一些控制参数有1）由路面附着系数与滑移率SB的关系曲线所表示的轮胎模型；2）滑移率控制上限SMAX、下限SMIN；3）车辆模型参数及初始车速O；4）制动器油压增长率KI和减小率KD等。根据分析可知，控制逻辑实现的关键是计算当前车轮滑移率SBT并与预先确定的上限值（SMAX，SMIN）进行比较，来判断对制动液压控制系统的增压或减压操作，控制流程如图3所示。图3仿真流程5实例分析单轮制动动力模型参数由表1给出。设式图2定义的路面附着系数分别为H08，G06以门限值控制算法设计ABS控制器，使车轮滑移率SB保持在最优值（SOPT附近），这里令SMIN018，SMAX022。根据表1给出的模型参数及附着系数，按照图3所示的控制流程采用M语言编制仿真程序。需要指出的是，表1给出的制动系统控制参数仅作为参考，系统设计过程中可根据需要适当调整，已获得满意的结果。表1单轮ABS制动力学模型参数参数符号单位数值车轮质量MKG300车轮动力半径DR025车轮转动惯量WI2KGM12初始车速O/S30初始角速度0/RAD120初始制动力矩BOTNM600制动油压增长率IK/S4500制动油压减小率D/5000采样时间STS0055MATLAB仿真过程图4MATLAB文件编辑调试窗口点击MATLAB指令窗工具条上的NEWFILE图标，打开如图4所示的MATLAB文件编辑调试器，其窗口名为UNTITLED1，我们可在空白窗口中编写程序。输入如下一段程序M300车轮质量RD025车轮动力半径IW12车轮转动惯量VWO30初始车速WO120初始角速度TBO600初始制动力矩KI4700制动油压增长率KD5000制动油压减少率TS005采样时间I1设置数组变量W1WO设置角速度变量V1VWO设置车速变量TB1TBO设置制动压力变量WHILEVI0如果车速大于零，则ABS工作SBIABSVIRDWI/VI计算当前滑移率IFSBI022滑移率大于上限，制动器减压TBI1TBITSKDENDIFSBI018滑移率处于上下限范围内，制动器保压TBI1TBIENDIFSBI018滑移率小于下限，制动器增压TBI1TBITSKIENDV2IWIRD计算车轮线速度V1IFXB/M计算车轮加速度VI1VIV1ITS计算下一采样周期的车轮前进速度WI1WITSFXBRDTBI1/IW计算下一采样周期的车轮角速度II1数组变量增加ENDX0I2绘制ABS控制的滑移率时域结果PLOTX,SB点击编辑调试器工具条的图标，在弹出的“保存为”对话框中，选择保存文件夹，键入新编文件名CARABS，点动保存键，完成文件保存。使CARABSM所在目录成为当前目录或让该目录处在MATLAB的搜索路径上，如该文件放在G盘，则在MATLAB主菜单中应出现点击编辑调试器工具条的DEBUG菜单，选择RUN指令，运行CARABS文件，可得到图5图5ABS控制的滑移率时域仿真结果修改CARABSM文件，输入以下指令V2I0X0I1PLOTX,V,X,V2运行后得到图6图6车轮前进速度与车轮线速度关系曲线（二）行驶动力学计算机建模、仿真及主动悬架控制器设计以单轮车辆模型为例，介绍行驶动力学计算机建模、仿真分析以及利用线性二次最优控制理论进行主动悬架LQG控制器设计过程。1计算机仿真系统模型的建立根据图7所示的主动悬架单轮车辆模型，运用牛顿运动定律，建立系统的运动方程，即（4）BASBWMXUKX（5）WSTGXX这里，采用一个滤波白噪声作为路面输入模型，即（6）002GGXTFXTGUT式中，XG为路面垂向位移（M）；G0为路面不平度系数（M3/CYCLE）；U为车辆前进速度（M/S）；W为数字期望为零的高斯白噪声；F0为下截止频率（HZ）。图7单轮车辆模型结合式（4）、式（5）和式（6），将系统运动方程和路面输入方程写成矩阵形式，即得出系统的空间状态方程（7）XABUFW式中，为系统状态矢量；W（W（T），为高斯白噪TBWBGXXX声输入矩阵；U（UA（T），为输入控制矩阵；000102SSBBSTSAWWKMKAF1BWMB002FGU2LOG控制器设计车辆悬架设计中的主要指标包括代表轮胎接地性的轮胎动载荷；代表轮胎舒适性的车身垂向振动加速度；影响车身姿态且与轮胎布置有关的悬架动行程。因此，LQG控制器设计中的性能指标J即为轮胎动位移、悬架动行程和车身垂向振动加速度的加权平方和在时域T内的积分值，其表达式为（8）222130LIMTWGBWBTJQXTTQXTTQXD式中，Q1、Q2和Q3分别为轮胎动位移、悬架动行程和车身垂向振动加速度的加权系数。加权系数的选取决定了设计者对悬架性能的倾向，如对车身垂向振动加速度项选择较大的权值，则考虑更多的是提高车辆操纵稳定性。为方便起见，这里取车身垂向振动加速度的加权值Q31。将性能指标J的表达式（8）改写成矩阵形式，即（9）01LIM2TTTTJXQURXND式中，；222121000SSBBSSKQQMQQ2BRM00ASNK当车辆参数值和加权系数值确定后，最优控制反馈增益矩阵可有黎卡提（RICCATI）方程求出，其形式如下（10）10TTPABNRPQ最优反馈控制增益矩阵，由车辆参数和加权系数决定。根据K任意时刻的反馈状态变量X（T），就可得到T时刻作动器的最优控制力UA，即（11）AUTT3计算实例这里，以某轿车的后悬架为例，给出一个完整的计算实例，包括车辆模型参数、仿真路面输入参数、控制器的设计参数以及计算结果。此例中车辆以20M/S的速度在某典型路面上行驶，仿真时间T50S。计算中输入的各参数及数值详见表2。表2单轮车辆模型仿真输入参数值车辆模型参数符号单位数值簧载质量MBKG320非簧载质量悬架刚度轮胎刚度悬架工作空间MWKSKTSWSCKGN/MN/MMM402000020000010仿真路面输入参数符号单位数值路面不平度系数车速下截止频率G0UF0M3/CYCLEM/SHZ50X1062001性能指标加权参数符号单位数值轮胎动位移悬架动行程车身加速度Q1Q2Q38000051仿真计算中以式（6）所示的滤波白噪声作为路面输入模型。白噪声的生成可直接调用MATLAB函数WGN（M，N，P）（此函数需要安装信号处理工具箱COMMUNICATIONSTOOLBOX），其中M为生成矩阵的行数，N为列数，P为白噪声的功率（单位为DB）。本例中取M10001，N1，P20。这意味着仿真计算中去一条白噪声，共10001个采集点，噪声强度为20DB。设定采样时间为0005S、车速为20M/S时，相当于仿真路面长度为1000M，仿真时间为50S。根据建立的系统状态方程式（7）及最优化性能指标函数式（9），利用已知的矩阵A、B、Q、R、N，调用MATLAB中的线性二次最优控制器设计函数K，S，ELQR（A，B，Q，R，N），即可完成最优主动悬架控制器的设计。输出的结果中，K为最优控制反馈增益矩阵，S为黎卡提方程的解，E为系统闭环特征根。根据表2给出的仿真输入参数，本例中求得的最优反馈增益矩阵K为K（7118812415192842038520864）同时，还得到了黎卡提方程的解245902847586073969247518630924510376937S在相同的仿真条件下，可将所设计的主动悬架系统与一个被动系统进行对比分析。在被动悬架系统中，取悬架刚度KS22000N/M，阻尼系数CS1000NS/M。除此之外，其他输入参数值均与主动悬架系统完全相同。4MATLAB仿真过程1）生成路面输入模型图8路面输入模型生成右击MATLAB中WORKSPACE栏中的ROAD_FILE，点击SAVESELECTIONAS，保存数据为MATLABMAT文件。图9路面输入模型保存2）参数输入新建一个MFILE，输入如下一段程序LOADMATLABMAT载入路面数据模型KS22000MB320KT200000MW40F001G00000005U20KB20000KS122000CS1000输入仿真有关参数A0,0,KS/MB,KS/MB,0建立主动悬架的状态矩阵0,0,KS/MW,KTKS/MW,KT/MW1,0,0,0,00,1,0,0,00,0,0,0,2PIF0A1CS/MB,CS/MB,KS1/MB,KS1/MB,0建立被动悬架的状态矩阵CS/MW,CS/MW,KS1/MW,KTKS1/MW,KT/MW1,0,0,0,00,1,0,0,00,0,0,0,2PIF0B1/MB,01/MW,00,00,00,2PISQRTG0UB10,00,00,00,00,2PISQRTG0UC1,0,0,0,00,1,0,0,00,0,1,0,00,0,0,1,00,0,0,0,1D0,00,00,00,00,0K71188,12415,19284,20385,20864K10,0,0,0,0点击编辑调试器工具条的DEBUG菜单，选择RUN指令，运行上述文件。3）用SIMULINK创建仿真框图（1）单击MATLAB工具栏上的图标，或在命令窗口中输入SIMULINK，打开SIMULINK模型库，如图10所示图10SIMULINK模型库（2）现在建立一个没有控制环节的仿真模型，来体验一下如何利用SIMULINK来达到目标。单击SIMULINK模型库工具条上的，新建一个后缀名为MDL的SIMULINK模型文件，这里取名DONGLIXUEMDL（名字可以任你选），并单击该模型文件窗口中的工具栏的按钮保存。这个过程如图11所示图11建立SIMULINK模型文件（3）确定输入环节，中间计算环节和输出环节所谓输入环节就是产生一个输入信号发送给系统。根据前面建立的模型，系统的输入是由路面输入和控制力输入组成的列向量。SIMULINK模型库提供了丰富的信号源。在模型库左侧的目录中，单击SIMULINK下的SOURCE，就可以看到各式各样的信号发生器，如图12所示图12SIMULINK提供的信号源在SOURCE中有一个从工作空间取值的信号源SIMIN，用鼠标左键点住它，然后将其拖到模型文件窗口。在模型文件窗口双击这个信号源，一个BLOCKPARAMETERS（块参数）对话框跳出来，将其参数改成图13所示图13块属性对话框至于控制力，由于这里建立的是有控制环节的系统仿真模型，所以单击SIMULINK下的MATHOPERATIONS，选择其中的MATRIXGAIN，将其拖到模型文件窗口。在模型文件窗口双击这个信号源，将其参数改成图14所示图14控制力参数选择接下来是确定中间计算环节。在模型库SIMULINK目录下有一个CONTINUOUS模型库，选择STATESPACE模型（如图15所示），将其拖入模型文件即可。双击模型文件中的STATESPACE模型，在块参数对话框中自定义状态方程的属性。由于我们已经在WORKSPACE中定义了A，B，C，D，所以在对话框中对应填入即可，如图16所示。图15系统的数学模型选择图16系统模型参数设置最后是输出环节。输出的数据和曲线可以存放在模型文件中也可以保存在工作空间供其他程序调用。在模型库SIMULINK目录下点击SINKS，就可以看到很多的输出工具。这里我们先使用示波器SCOPE（图17、18），将其拖入模型文件，双击模型文件窗口的示波器，在弹出的图形框的工具栏中单击，会有一个对话框出现，我们可以在此处订制示波器的属性。选择DATAHISTORY的选项卡，将LIMITDATAPOINTSTOLAST前面的钩取掉即可（如图19所示）。图17示波器图18双击SCOPE后出现的图像框图19示波器参数设置接下来我们将输入模块、输出模块和中间计算模块连接起来。由于输入有两个必须使用信号路由器（SIGNALROUTING）来进行合成。如图20所示，选择MUX，将其拖入模型文件，双击，更改其属性（如图21所示）。其中NUMBERSOFINPUTPORT设置输入端口的数目，DISPLAYOPTION设置MUX的形状。同样系统输出的列向量有5个分量，选用DEMUX可以将其分解成五路输出（如图22所示），从上到下依次对应输出向量的各列分量。由于我们选定的悬架评价指标为轮胎动位移、悬架动行程和车身垂向振动加速度，因此必须对输出各列向量进行微分和加减运算。微分运算在SIMULINK模型库中选择CONTINUOUS中的DERIVATIVE，加减运算在SIMULINK模型库中选择MATHOPERATIONS中的SUM，将其拖入模型文件中，双击，更改加号为减号。（如图23所示）模块之间的连线是很简单的事情。将鼠标靠近模块的输出端口，会出现一个十字光标，按住鼠标左键拖动到另一个端口的输入端口再松开鼠标左键就完成了两个模块的连接。图20选择SIGNALROUTING里的MUX图21更改MUX的属性图22更改SUM的属性图图23更改DEMUX的属性为了与主动悬架性能比较，我们需要建立一个被动悬架的模型。被动悬架模型建立过程与上述主动悬架建立相似，只是参数设置方面有所不同。（如图24，25所示），最后我