《汽车运动控制系统设计中的Matlab应用.docx》由会员分享,可在线阅读,更多相关《汽车运动控制系统设计中的Matlab应用.docx(24页珍藏版)》请在优知文库上搜索。
1、1绪论1.1选题背景与意义汽车已经成为人们平常生活不可缺乏时代步交通工具,在汽车发达国家,旅客运送的60%以上,货品运送05O%以上由汽车来完毕,汽车工业水平和家庭平均拥有汽车数量已经成为衡量一种国家工业发达程度的标志。进行汽车运动性能研究时.一般从操纵性、稳定性和乘坐舒适性等待性着手。但近年来.伴随交通系统的日趋复杂,考虑了道路环境在内的汽车运动性能开始受到关注。因此,汽车运动控制系统的研究也显得尤为重要,在文中,首先对汽车B运动原理进行分析,建立控制系统简化模型,确定期望的静态指针(稳态误差)和动态指针(超调量和上升时间)。然后对汽车运动控制系统进行设计分析。从而确定系统0最佳静态和动态指
2、针。2论文基本原理分析汽车运动横向控制(1)绝对位置的获得措施汽车横向方向g控制使用GPS(全球定位系统)0绝对位置信息。GPS信息B精度与采样周期、时间滞后等有关。为提高GPSB数据精度和平滑数据.采用卡尔曼滤波对采样数据进行修正。GPS的采样周期为200InS相对应控制的周期采用50ms。此外考虑通信等时滞后、也需要进行赔偿,采用航位推测法(deadreckoning)处理此问题。通过卡尔曼滤波和航位推测法推算出时值作为汽车时绝对位置使用来控制车速、横摆角速度等车辆的状态量。GPS的数据通过卡尔曼滤波减少偏差、通过航位推测法进行误差和迟滞赔偿.提高了位置数据推算的精度。(2)前轮转角变化量
3、的算出措施这里对前轮目的转角变化量(Ab)的算出措施作简要阐明,横方向控制采用预见控制,可以从目前汽车B状态预测通过时间tp秒后0汽车位置,由tp秒后的预测位置和日日勺途径的位置可以算出tp秒后为沿着目日勺途径行驶所需要B汽车横摆角速度(3。这个数值前回馈或者从与目前值0目的途径0误差的反馈来推算前轮目的转角变化量(式(1).k=(Z1rk+k2rc)Tc式(1)式中TC为控制周期,k1,k2根据与目前目的途径的误差()最小的原则来求解。汽车运动纵方向的控制建立一种合理的传动系统模型是设计高性能汽车纵向运动控制系统的基础。目前纵向运动控制器设计过程中采用的传动系统简化模型重要有两类:一类是忽视
4、传动系统口勺部分动态特性得到简化模型:另一类是通过对输入输出特性辨识得到简化模型。本文借鉴文献,忽视传动系统的部分动态特性,将车辆简化为两轮模型,对于前轮驱动车辆,整车受力如图1所示。前后车轮运动方程分别为(Jff=Ts-rWff-rFfIJrr=rFr-rWrf上式中Jf和Jr,分别为前后轮转动惯量(左右轮之和),再为后轮转速,Wr和Wf分别为前后车轮日勺垂直载荷(左右轮之和),Fr和Ff分别为前后轮切向力(左右轮之和),r为车轮半径,f为滚动阻力系数。对于汽车纵向运动控制系统,不会出现非常大的加减速度,采用线性化轮胎模型,得到切向力与滑移率关系为:F=KwSo式中KW为轮胎纵向刚度,S为滑
5、移率。驱动时s=lr/(r),制动时S=I-(r)/V。认为风阻FW作用于汽车质心,则前后轮垂直载荷分别为bhgWf=Mg-Mv-LLahWr=M2-+Mv1L式中M为整车品质;a和b分别为前后轴到质心的距离,L=a+b;hg为质心至地面的高度。整车运动方程为MV=Ff-Fr-CQAV2,式中Ca风阻系数,A为等效迎风面积。发动机转矩、发动机转速、涡轮转速、半轴转矩、前后轮转速、车速7状态的非线性传动系统模型,在低频带内,发动机动态对传动系统特性基本无影响,假如控制系统只波及较低频段.可以忽视发动机动态。忽视了半轴、轮胎滑移以及载荷转移和发动机转矩,只包括发动机转速、车速2个状态。飞轮运动方程
6、为Jee=Tedes-Tp整车运动方程为Mv=RgRd-CAv2-Mgf2.2汽车运动控制系统的模型简化分析考虑图2所示的汽车运行控制系统。假如忽视车轮时转动惯量.并且假定汽车受到的摩擦阻力大小与运动速度成正比,方向与汽车运动方向相反,则该系统可以简化成简朴的质量阻尼系统。根据牛顿运动定律,该系统0模型(亦即系统B运动力方程)表达为:(m三+hv=u(3-1)Iy=v其中,U为汽车0驱动力。为了得到控制系统的传递函数,对式(3T)进行LaPIaCe变换。假定系统的初始条件为零,则动态系统的LaPlaCe变换为:CmsV(s)+bP(s)=U(三)Iy(s)=U(三)由于系统输出是汽车的运动速度
7、,用Y(三)替代V(三),得到:msY(s)+bY(三)=U(三)因此.汽车运动控制系统模型的传递函数为:V(三)_1Z=.U(三)ms+b2.3汽车控制系统PID控制器的校正根据阶跃响应曲线.运用串联校正的原理.以及参数变化对系统响应0影响,对静态和动态性能指针进行详细的分析,最终设计出满足我们需要的控制系统。系统在未加入任何校正环节时的开环传递函数,在MATLAB环境下对系统未加校正时开环阶跃响应曲线进行仿真.绘制如图3阶跃响应曲线,图中系统的开环响应曲线未产生振荡,其上升时间约100秒,稳态误差到达98%,远不能满足跟随设定值的规定。图4(1)首先选择P校正,也就是在系统中加入一种比例放
8、大器,为了大幅度减少系统日勺稳态误差,同步减小上升时间。P校正后系统的闭环传递函数为:Xs)二KPU(三)ms+(b+Kp)此时控制系统的稳态值为Kp(b+KP)=Kp/(50+Kp)。本系统B比例增益KP=800。即稳态值为800/(50+800)=0.941,这样可以把系统0稳态误差减少到0.06左右。加入P校正后控制系统0死循环阶跃响应曲线如图3所示。图中,系统的稳态值约为0.941.稳态误差约为5.9%,这和最初B设计规定仍有差距,并且上升时间在7秒左右,不能到达设计的需要。因此我们选择Pl校正。(2)加入Pl校正器后系统的闭环单位反馈传递函数为:Y(三)_KPS+KiU(三)ms2+
9、(b+KP)S+Ki考虑到Ki日勺作用.我们可以大幅度减少kp,取kp=2000Ki=70,在MATLAB环境下仿真得出的系统响应曲线如图4(中)所示。从图4(中)中可以得知,加入Pl校正后系统的上升时间有所下降,但仍不小于5秒。同步又产生了另一种问题,系统0超调量到达了26.43%.这是使用积分器带来的副作用。因此合适地加入微分量。(3)可以选择PD校正,此时系统的闭环单位反馈传递函数为:丫_KDS+pi(s)=(m+Kd)s+(h+Kp)鉴于KD对上升时间和稳态误差影响不大.我们在P校正的基础上.将KD减少少许,给出KD=I0。系统响应曲线如图4(中)所示。(4)加入PlD校正,此时系统的
10、闭环单位反馈传递函数为:y(s)=KDS2+KPS+KtU(三)(n+KD)S2+(b+KP)S+KiKp,Ki和KD的选择一般先根据经验确定一种大体的范围,然后通过MATLAB绘制的图形逐渐校正。这里我们取Kp=700,Ki=100,Kd=100o得到加入PLD校正后系统的死循环阶跃响应如图4(右)所示。从图4(右)中可以得出,系统的静态指针和动态指针,已经很好日勺满足了设计的规定。上升时间不不小于5s,超调量不不小于8%,约为6.67%0图5根据系统0性能指针和某些基本的整定参数的经验,选择不一样BPID参数进行模拟,最终确定满意的参数。这样做首先比较直观,另首先计算量也比较小,并且便于调
11、整。2.4汽车运动控制系统根轨迹校正的设计过程为了减小系统的稳态误差,同步尽量减小超调量和上升时间B变化,到达满意的效果,我们需要从相位的角度来考虑,变化控制器的构造,从而想到相位滞后器的作用。相位滞后器的传递函数为:GC(S)=s + Zo s+ Po这样.整个系统的死循环传递函数就变成了:Y(三)KPS+KPZoU(三)ms2+(b+mP0+KP)S+(bP0+KPZO)滞后控制器0零极点应设计成紧靠在一起,这样控制系统0稳态误差将减小Z0Po倍。根据上面的分析,将ZO设计成-0.3,而PO等于-0.03。图6得到的根轨迹如图7中。在实轴的-0.35的位置附近选择期望点,得到图7所示的系统
12、阶跃响应曲线。从图7中可以得知,这时的稳态误差已经满足设计规定。出现的少许超调亮是加入之滞后控制器的成果。死循环系统的超调量约为7.64%,满足不不小于8%0设计规定,上升时间约为2.5秒,以及稳态误差都已经满足设计规定。3对论文采用的理论和措施进行研究本论文运用MATLAB对简化后的汽车运动控制系统进行仿真,由于文中没有详细过程,图形也不能辨别精确值,扩写时我进行详细分析并按照自己的理解进行仿真。文中简化后B汽车运动控制系统的开环传递函数为供=一三,其开环传递为一阶惯U(三)ms+b性系统。而全文没有提及汽车的质量m(通过背面的仿真,选用In值为800。)由于文中图形0辨别率问题,不能从文中
13、读出精确值,仿真成果只能靠近源图形,但已经足够完毕规定,即对汽车运动控制简化模型的PlD校正。汽车运动控制简化模型传递函数仿真设计对原开环传递函数=高=中O(三)运用MATLAB进行单位阶跃输入响应的仿真。仿真程序如下:b=50;m=800;t=0:0.1:120;y=l;u=mb;sysO=tf(y,u);y1,t=step(sysO,t);sysl;plot(t,yl);grid;xlabel(Time(seconds)1),ylabel(,StepResponse1)仿真成果图形如图8,图中上升时间明显偏大,大概60秒,并且稳态误差有98%,远远不能满足论文中的规定,但原文中没有对规定进
14、行统一,所如下文中我选定上升时间不不小于5秒,超调量8%,稳态误差不不小于2%。0.02O1020406080100120Time (seconds).O6.O O.0.O.0.00.0.0 O O O O SUodS H d图8闭环传递函数(s)单位阶跃输入响应汽车运动控制系统P校正函数仿真设计论文对开环传递函数进行PID校正,文中是通过三步尝试得到最终PID校正参数。首先要减小系统的上升时间,进行P校正,即在开环系统中加入比例放大环节Kp,P校正后系统的闭环传递环数为p(s)=Pns+(b+KP)按文中数据取kp=800,原系统b=50,m=800运用MATLAB进行闭环系统I向单位阶跃输
15、入响应仿真。仿真程序如下:kp=800;b=50;m=800;t=0:0.1:7;y=kp;u=mb+kp;sys1=tf(y,u);yl,t=ste(sysl,t)jsys1;lot(t,yl);grid;xlabel(,Time(seconds),ylabel(,StepResponse1)仿真成果图形如下图90.90.80.70.60.50.40.30.20.1234Time (seconds)67esuods B d图9闭环传递函数p(s)单位阶跃输入响应详细分析:令(三)=Pk7ms+(b+Kp)Ts+1比较系数得T=1617,一阶系统的阶跃响应是一种按指数规律单调上升的过程,其动态性能指标中不存在超调量、峰值时间、上升时间等项。按一阶系统的过渡过程时间定义:ts=3T,计算得ts=2.82s,