制动系统建模仿真及ABS控制器设计

1、目录1. 动力学建模-.0-2. 分段线性的轮胎模型-.0-3. 控制算法-.1-4. 仿真流程及参数输入-.1-5. 实例分析-.2-6. MATLAB仿真过程-.2-6.1. 逻辑门限值控制器-.2-6.2. 模糊控制器-.6-6.2.1模糊控制器设计-.6-622模糊控制matlab仿真过程7-6.3. PID控制器-11-通过两自由度单轮模型为例，介绍在MATLAB环境下的控制系统力学建模、ABS控制器设计及仿真分析过程。1.动力学建模某车辆简化后的制动力学模型如右图所示。其中单轮模型质量转动惯量为IW,车辆旋转角速度为3,车轮轮心前进速度为UW，地面制动力为Fxb。作用于车轮的制动力

2、矩为Tb。若忽略空气阻力与车轮滚动阻力，则系统的运动方程如下：Tb(1)(2)Fxb等于地面作用于车轮的法向反力Fx的乘积，其中为制动滑移率Sb的函数。IWddtdum尸加dt式中，地面制动力与路面附着系数2.分段线性的轮胎模型根据第三章中介绍的有关轮胎纵向特性的内容，路面附着系数与车轮滑移率之间存在一定的非线性关系。如果用两段直线近似表示路面附着-滑移曲线，可得到分段线性化的附着系数卩与车辆滑移率Sb的关系(即所谓的“Dugoff”模型),如下图所示。其表达式如下：Shgohg1So1So式中，yb为峰值附着系数；yg为车轮完全抱死时(即Sb=1)时的路面附着系数；So为峰值附着系数所对应的

3、滑移率。SbSo(3)也援柳晤轲SE图2线性化的路面附着系数与车轮滑移率关系曲线3. 控制算法这里以门限值控制算法为例，说明ABS控制器设计及制动系统力学的仿真过程。采用门限制控制算法的基本思想是保证车轮滑移率在最理想的范围之内。制动开始后，随着制动压力的升高车轮转速co相应减小，车轮出现滑移；当车轮滑移率达到理想范围上限值SmaxB寸，减小制动压力；随着制动压力的减小，直至减小到滑移率下限值Smin时再增大制动压力。循环往复这一过程直至车辆停止。因此，在ABS控制器起作用的过程中，滑移率总是保持在理想的范围内，从而保证车辆的最佳制动性能及行驶方向控制的稳定性。4. 仿真流程及参数输入由上可知

4、，ABS控制器所用到的一些控制参数有：1）由路面附着系数卩与滑移率Sb的关系曲线所表示的轮胎模型；2）滑移率控制上限Smax、下限Smin；3）车辆模型参数及初始车速卩so；4）制动器油压增长率ki和减小率kd等。根据分析可知，控制逻辑实现的关键是计算当前车轮滑移率Sb（t）并与预先确定的上限值（Smax，Smin）进行比较，来判断对制动液压控制系统的增压或减压操作，控制流程如图3所示。图3仿真流程5. 实例分析单轮制动动力模型参数由表1给出。设式图2定义的路面附着系数分别为曲=0.8，旧=06以门限值控制算法设计ABS控制器，使车轮滑移率Sb保持在最优值（Sopt附近）,这里令Smin=0.

5、18,Smax=0.22。根据表1给出的模型参数及附着系数，按照图3所示的控制流程采用m语言编制仿真程序。需要指出的是，表1给出的制动系统控制参数仅作为参考，系统设计过程中可根据需要适当调整，已获得满意的结果。表1单轮ABS制动力学模型参数参数符号单位数值车轮质量mkg300车轮动力半径rdm0.25车轮转动惯量Iwkgm212初始车速wom/s30初始角速度0rad/s120初始制动力矩TboNm600制动油压增长率kiNm/s4500制动油压减小率kdNm/s5000采样时间s0.056. MATLAB仿真过程6.1.逻辑门限值控制器Matlab代码如下：clearclc%输入初始参数m=

6、300;%车轮质量rd=0.25;%车轮动力半径逻辑门控制Iw=12;%车轮转动惯量u0=30;%初始车速w0=120;%初始角速度Tb0=600;%初始制动力矩ki=4500;%制动油压增长率kd=5000;%制动油压减少率dt=0.05;%采样时间u(1)=u0;%设置车速数组变量w(1)=w0;%设置轮速数组变量Tb(1)=Tb0;%设置制动压力数组变量s(1)=0;%设置滑移率数组变量Fxb(1)=0;%设置地面制动力数组变量%滑移率逻辑门控制i=1;while(u(i)>0)s(i+1)=(u(i)-w(i)*rd)/u(i);%计算滑移率ifs(i+1)<0.2mu=4

7、*s(i+1);elsemu=0.85-0.25*s(i+1);endFxb(i+1)=m*9.8*mu;ifs(i+1)<0.18Tb(i+1)=Tb(i)+ki*dt;elseifs(i+1)>0.22Tb(i+1)=Tb(i)-kd*dt;elseTb(i+1)=Tb(i);enddw=(Fxb(i+1)*rd-Tb(i+1)/Iw;du=-Fxb(i+1)/m;u(i+1)=u(i)+du*dt;w(i+1)=w(i)+dw*dt;i=i+1;endt=0:dt:(i-1)*dt;plot(t,u'w'*rd)figure()%计算纵向附着系数%计算地面制动

8、力%与滑移率门限值比较，确定制动力矩%计算车轮角加速度%计算车身加速度%计算下一时刻车速%计算下一时刻轮速plot(t,Fxb,t,Tb)figure。plot(t,s)运行结果如下：1)制动器制动力矩与地面制动力随时间变化曲线力动制-50000.511.533.54制动力变化曲线0020015000O005O-地面制动力N-制动器制动力矩N*m2)滑移率随时间变化曲线0.3522.5时间t/s4.50.30.250.2率移滑滑0.150.10.0500.51滑移率变化曲线1.533.5422.5时间t/s03）车速与轮速随时间变化曲线30车速与轮速变化曲线Sm度速25201510-50.51

9、.52.53.5时间t/s4.5V、:速:速弋、轮、7、*04可以看出：基于滑移率的逻辑门限值控制器，可以限制滑移率在给定的最佳滑移率附近波动，汽车完全停止时，车速为零，滑移率趋于无穷大，所以在仿真时间最后一段，滑移率有较大波动。制动器制动力矩在滑移率控制器的控制下，可以完成增压、保压、减压动作，防止车轮过早抱死，减小制动距离。62模糊控制器621模糊控制器设计1）输入量模糊化为简化，仅选取滑移率偏差e作为输入量，当最佳滑移率选为0.2时，滑移率偏差范围为【-0.2,0.8】，经过适当变化，使e的变化范围为【-6,6】，以便于控制。为简化，仅选取三个语义词【N,Z，P】，即正、零、负。分别定义

10、各自的隶属度函数，进行模糊化，如下图：4X0inputa"2）输出量反模糊化选取制动器制动力矩改变量为输出量，取值范围设定为【-1，+1】。同样，为简化，选取三个语义词【I,K，D】，即增压、保压、减压。并定义相应隶属度函数，如下图：3）模糊规则由于仅有一个输入量，且仅选取三个语义词，只需定义三条模糊规则，如下图:4）控制曲线观察制动器制动力矩改变量与滑移率偏差的关系0.3622模糊控制matlab仿真过程程序如下：clearclc%输入初始参数m=300;%车轮质量rd=0.25;%车轮动力半径lw=12;%车轮转动惯量u0=30;%初始车速w0=120;%初始角速度Tb0=600

11、;%初始制动力矩ki=4500;%制动油压增长率kd=5000;%制动油压减少率dt=0.05;%采样时间u(1)=u0;%设置车速数组变量w(1)=w0;%设置轮速数组变量Tb(1)=Tb0;%设置制动压力数组变量s(1)=0;%设置滑移率数组变量Fxb(1)=0;%设置地面制动力数组变量dtp(1)=0;%设置地面制动力变化量数组变量%模糊控制i=1;while(u(i)>0)s(i+1)=(u(i)-w(i)*rd)/u(i);ifs(i+1)<0.2%计算滑移率%计算纵向附着系数模糊控制mu=4*s(i+1);elsemu=0.85-0.25*s(i+1);endFxb(i

12、+1)=m*9.8*mu;e=12*(s(i+1)-0.5);abs=readfis('abs_fzy2.fis'dtp(i+1)=evalfis(e,abs);Tb(i+1)=Tb(i)+2*5000*dtp(i+1)*dt;dw=(Fxb(i+1)*rd-Tb(i+1)/Iw;du=-Fxb(i+1)/m;u(i+1)=u(i)+du*dt;w(i+1)=w(i)+dw*dt;i=i+1;end%绘图t=0:dt:(i-1)*dt;plot(t,u'w'*rd)figure()plot(t,Fxb,t,Tb)figure()plot(t,s)figure()

13、plot(t,dtp)%计算地面制动力%将滑移率偏差范围转化到【);%调用设计好的模糊控制器%输出制动器制动力矩变化量%计算制动器制动力矩运行结果：1）模糊控制器输出，即Tb的变化量模糊控制器输出（Tb变化量）2）地面制动力与制动器制动力矩制动力变化曲线Fb(N)-Tb(N*m)00O50050T00T力动制00-544.533.522.5时间t/s3)滑移率滑移率变化曲线4)车速与轮速车速与轮速变化曲线302520S15m度速1050时间t/s-56.3.PID控制器PID控制器的设计选取滑移率s的偏差e作为输入，若选定最佳滑移率为0.2,则e=s-0.2;输出量为制动器制动力矩的变化量dt

14、p。则有，dtpKp*eedtKd*詈选取一组合适的系数，Kp、Ki、Kd，即可完成PID控制MATLAB仿真过程MATLAB程序代码如下：clearclc%输入初始参数m=300;%车轮质量rd=0.25;%车轮动力半径Iw=12;%车轮转动惯量u0=30;%初始车速w0=120;%初始角速度Tb0=600;%初始制动力矩ki=4500;%制动油压增长率kd=5000;%制动油压减少率dt=0.05;%采样时间u(1)=u0;%设置车速数组变量w(1)=w0;%设置轮速数组变量Tb(1)=Tb0;%设置制动压力数组变量s(1)=0;%设置滑移率数组变量Fxb(1)=0;%设置地面制动力数组变

15、量dtp(1)=0;%设置地面制动力变化量数组变量%PID控制器i=1;while(u(i)>0)s(i+1)=(u(i)-w(i)*rd)/u(i);ifs(i+1)<0.2mu=4*s(i+1);elsemu=0.85-0.25*s(i+1);endFxb(i+1)=m*9.8*mu;PIDe=s(i+1)-0.2;%计算滑移率偏差控ec=(s(i+1)-s(i)/dt;%计算滑移率偏差变化率制器es=sum(s);%计算滑移率偏差的累积dtp(i+1)=Kp*e+Kd*ec+Ki*es;%计算PID控制输出量，即Tb变化量Tb(i+1)=Tb(i)-800*dtp(i+1)*dt;%计算Tbdw=(Fxb(i+1)*rd-Tb(i+1)/Iw;du=-Fxb(i+1)/m;u(i+1)=u(i)+du*dt;w(i+1)=w(i)+dw*dt;i=i