汽车运动控制系统设计中的Matlab应用.docx

《汽车运动控制系统设计中的Matlab应用.docx》由会员分享，可在线阅读，更多相关《汽车运动控制系统设计中的Matlab应用.docx（24页珍藏版）》请在第一文库网上搜索。

1、1绪论1.1选题背景与意义汽车已经成为人们平常生活不可缺乏日勺代步交通工具，在汽车发达国家，旅客运送日勺 60%以上，货品运送日勺50%以上由汽车来完毕，汽车工业水平和家庭平均拥有汽车数量已 经成为衡量一种国家工业发达程度日勺标志。进行汽车运动性能研究时.一般从操纵性、稳 定性和乘坐舒适性等待性着手。但近年来.伴随交通系统日勺日趋复杂，考虑了道路环境在 内日勺汽车运动性能开始受到关注。因此，汽车运动控制系统日勺研究也显得尤为重要，在文 中，首先对汽车日勺运动原理进行分析，建立控制系统简化模型，确定期望日勺静态指针（稳 态误差）和动态指针（超调量和上升时间）。然后对汽车运动控制系统进行设计分析。

2、从而 确定系统日勺最佳静态和动态指针。2论文基本原理分析汽车运动横向控制（1）绝对位置日勺获得措施汽车横向方向日勺控制使用GPS （全球定位系统）日勺绝对位置信息。GPS信息日勺精度与采样 周期、时间滞后等有关。为提高GPS日勺数据精度和平滑数据.采用卡尔曼滤波对采样数据 进行修正。GPS日勺采样周期为200IilS相对应控制日勺周期采用50ms。此外考虑通信等日勺滞后、 也需要进行赔偿，采用航位推测法（dead reckoning）处理此问题。通过卡尔曼滤波和航位推测法推算出日勺值作为汽车日勺绝对位置使用来控制车速、横摆角速度等车辆日勺状态量。GPS 日勺数据通过卡尔曼滤波减少偏差、通过航位

3、推测法进行误差和迟滞赔偿.提高了位置数据 推算日勺精度。（2）前轮转角变化量日勺算出措施这里对前轮目日勺转角变化量S）日勺算出措施作简要阐明，横方向控制采用预见控制， 可以从目前汽车日勺状态预测通过时间tp秒后日勺汽车位置，由tp秒后日勺预测位置和目日勺途径 日勺位置可以算出tp秒后为沿着目日勺途径行驶所需要日勺汽车横摆角速度（j.这个数值前同 馈或者从与目前值日勺目日勺途径日勺误差日勺反馈来推算前轮目日勺转角变化量（式）.5fc = （fc1 rk + k2 , rfc） , Tc式式中TC为控制周期，k1, k2根据与目前目日勺途径日勺误差（。最小日勺原则来求解。汽车运动纵方向的J控制建立

4、一种合理日勺传动系统模型是设计高性能汽车纵向运动控制系统日勺基础。目前纵向运动控制器设 计过程中采用日勺传动系统简化模型重要有两类：一类是忽视传动系统日勺部分动态特性得到简化模型：另 一类是通过对输入输出特性辨识得到简化模型。本文借鉴文献，忽视传动系统日勺部分动态特性，将车辆简化为两轮模型，对于前轮驱动车辆，整车受力如图1所示。前后车轮运动方程分别为CJfj）f = Ts rWff rFf（Jrr = rFr rWrf上式中Jf和Jr，分别为前后轮转动惯量（左右轮之和），为后轮转速，Wr和Wf分别为 前后车轮日勺垂直载荷（左右轮之和），Fr和Ff分别为前后轮切向力（左右轮之和），r为车轮 半径

5、，f为滚动阻力系数。对于汽车纵向运动控制系统，不会出现非常大日勺加减速度，采 用线性化轮胎模型，得到切向力与滑移率关系为：F = KwSo式中KW为轮胎纵向刚度，S为滑移率。驱动时S=l-v/ （丫），制动时S=ITrco） /V。认为风阻FW作用于汽车质心，则前后轮垂直载荷分别为b hg Wf = Mg-Mv- LLa hgWr = Ms- + Mv6 L L式中M为整车品质；a和b分别为前后轴到质心日勺距离，L=a+b; hg为质心至地面日勺高 度。整车运动方程为MV=Ff-F,-CciAv2 ,式中CQ风阻系数，A为等效迎风面积。发动机转矩、发动机转速、涡轮转速、半轴转矩、前后轮转速、车

6、速7状态日勺非线性 传动系统模型，在低频带内，发动机动态对传动系统特性基本无影响，假如控制系统只波 及较低频段.可以忽视发动机动态。忽视了半轴、轮胎滑移以及载荷转移和发动机转矩, 只包括发动机转速、车速2个状态。飞轮运动方程为Jee = Tedes - Tp整车运动方程为Mv =变誉-CaAv2 - Mgf2.2汽车运动控制系统的模型简化分析考虑图2所示日勺汽车运行控制系统。假如忽视车轮日勺转动惯量.并且假定汽车受到日勺 摩擦阻力大小与运动速度成正比，方向与汽车运动方向相反，则该系统可以简化成简朴日勺质量阻尼系统。根据牛顿运动定律，该系统日勺模型(亦即系统日勺运动力方程)表达为：(dv , m

7、- bv = uZ0 1X at(3-1)、 y = V其中，U为汽车日勺驱动力。为了得到控制系统日勺传递函数，对式(3T)进行LaPIaCe变换。假定系统日勺初始条件为零，则动态系统日勺LaPIaCe变换为：fmsVr(s) + W(S) = U(S)( y(s) = u(s)由于系统输出是汽车日勺运动速度，用Y(S)替代V(S),得到：msY (S) + bY (S) = U (S)因此.汽车运动控制系统模型日勺传递函数为：Y(S)=1U (S) ms + b2.3汽车控制系统PID控制器的校正根据阶跃响应曲线.运用串联校正日勺原理.以及参数变化对系统响应日勺影响，对静态 和动态性能指针进

8、行详细日勺分析，最终设计出满足我们需要日勺控制系统。系统在未加入任 何校正环节时日勺开环传递函数，在MATLAB环境下对系统未加校正时开环阶跃响应曲线进行 仿真.绘制如图3阶跃响应曲线，图中系统日勺开环响应曲线未产生振荡，其上升时间约IOO 秒，稳态误差到达98%,远不能满足跟随设定值日勺规定。图3图4(1)首先选择P校正，也就是在系统中加入一种比例放大器，为了大幅度减少系统日勺稳态 误差，同步减小上升时间。P校正后系统日勺闭环传递函数为：Y(S)二KPU(S) TnS+ (b+ Kp)此时控制系统日勺稳态值为Kp(b + Kp) = Kp/ (50 + Kp)。本系统日勺比例增益KP = 8

9、00。即稳态值为800 / (50+800)=0.941,这样可以把系统日勺稳态误差减少到0.06左右。加入P校正后控制系统日勺死循环阶跃响应曲线如图3所示。图中，系统日勺稳态值约为0. 94L稳态 误差约为5.9%,这和最初日勺设计规定仍有差距，并且上升时间在7秒左右，不能到达设计 日勺需要。因此我们选择Pl校正。(2)加入Pl校正器后系统日勺闭环单位反馈传递函数为：y(s)_ ps + iU (S) ms2 + (b + Kp)S + Ki考虑到Ki日勺作用.我们可以大幅度减少kp, W = 200O Ki = 70,在MATLAB环境下仿 真得出日勺系统响应曲线如图4(中)所示。从图4(

10、中)中可以得知，加入Pl校正后系统日勺上升 时间有所下降，但仍不小于5秒。同步又产生了另一种问题，系统日勺超调量到达了26. 43%.这是使用积分器带来日勺副作用。因此合适地加入微分量。(3)可以选择PD校正，此时系统日勺闭环单位反馈传递函数为：y(s) = KDS + KP = (m + KD)S + 3 + KP)鉴于KD对上升时间和稳态误差影响不大.我们在P校正日勺基础上.将KD减少少许，给 出KD=I0。系统响应曲线如图4(中)所示。(4)加入PID校正，此时系统日勺闭环单位反馈传递函数为：y(s) _ Kds2 + Kps + KiU (S) (m + KD)S2 + (b + KP

11、)S + KiKp, Ki和KD日勺选择一般先根据经验确定一种大体日勺范围，然后通过MATLAB绘制日勺图形 逐渐校正。这里我们取Kp=700, Ki=100, Kd = IOOo得到加入PLD校正后系统日勺死循环阶 跃响应如图4(右)所示。从图4(右)中可以得出，系统日勺静态指针和动态指针，已经很好日勺 满足了设计日勺规定。上升时间不不小于5s,超调量不不小于8%,约为6. 67%o11t3m gW6CWTCMV i i4 x1丁苇 一一 一4 一efMM图5根据系统日勺性能指针和某些基本日勺整定参数日勺经验，选择不一样日勺PlD参数进行模拟,最终确定满意日勺参数。这样做首先比较直观，另首先

12、计算量也比较小，并且便于调整。2.4汽车运动控制系统根轨迹校正的设计过程为了减小系统日勺稳态误差，同步尽量减小超调量和上升时间日勺变化，到达满意日勺效果, 我们需要从相位日勺角度来考虑，变化控制器日勺构造，从而想到相位滞后器日勺作用。相位滞 后器日勺传递函数为：这样.整个系统日勺死循环传递函数就变成了：y(S)二KPS + KPZOU (S) ms2 + (h + TnPo + KP)S + (bP。+ KPZO)滞后控制器日勺零极点应设计成紧靠在一起，这样控制系统日勺稳态误差将减小Z0Po倍。根据上面日勺分析，将ZO设计成-0.3,而PO等于-0.03。图6图7得到日勺根轨迹如图7中。在实轴

13、日勺-0.35日勺位置附近选择期望点，得到图7所示日勺系统 阶跃响应曲线。从图7中可以得知，这时日勺稳态误差已经满足设计规定。出现日勺少许超调亮是加入之 滞后控制器日勺成果。死循环系统日勺超调量约为7. 64%,满足不不小于8%日勺设计规定，上升 时间约为2. 5秒，以及稳态误差都已经满足设计规定。3对论文采用的理论和措施进行研究本论文运用MATLAB对简化后日勺汽车运动控制系统进行仿真，由于文中没有详细过程， 图形也不能辨别精确值，扩写时我进行详细分析并按照自己日勺理解进行仿真。文中简化后日勺汽车运动控制系统日勺开环传递函数为祟=工，其开环传递为一阶惯U(S) ms+b性系统。而全文没有提及

14、汽车日勺质量m (通过背面日勺仿真，选用IIl值为800。)由于文中图 形日勺辨别率问题，不能从文中读出精确值，仿真成果只能靠近源图形，但已经足够完毕规 定，即对汽车运动控制简化模型日勺PID校正。汽车运动控制简化模型传递函数仿真设计对原开环传递函数 = 高=O(S)运用MATLAB进行单位阶跃输入响应日勺仿真。仿真程序如下：b=50;m=800;仁0:0.1:120;y=l;u=m b;sys=tf(y,u); y 1 ,t=step(sys,t) ；sys 1;plot(t,yl); grid;xlabel(,Time (seconds)*), ylabel(,Step Response1

15、)仿真成果图形如图8,图中上升时间明显偏大，大概60秒，并且稳态误差有98%,远 远不能满足论文中日勺规定，但原文中没有对规定进行统一，所如下文中我选定上升时间不不小于5秒，超调量8%,稳态误差不不小于2%。0.02图8闭环传递函数R)(S)单位阶跃输入响应 SUodSccd 汽车运动控制系统P校正函数仿真设计论文对开环传递函数进行PlD校正，文中是通过三步尝试得到最终PlD校正参数。首 先要减小系统日勺上升时间，进行P校正，即在开环系统中加入比例放大环节K” P校正后 系统日勺闭环传递环数为p(s)=P THS +3 + Kp)按文中数据取kp=800,原系统b=50, m=800o运用MATLAB进行闭环系统日勺单位阶跃 输入响应仿真。仿真程序如下：kp=800;b=50;m=800;t=0:0.1:7;y=kp;u=m b+kp;sys 1 =tf(y,u); y 1 ,t=step(sys 1 ,t) ；sys 1;plot(t,yl); grid;xlabel(,Time (seconds)1),