汽车操纵稳定性仿真 matlab

1、西安理工大学研究生课程论文/研究报告课程名称： 车辆系统动力学控制与仿真 课程代号： 任课教师： 傅卫平 论文/研究报告题目：汽车操纵稳定性仿真 完成日期： 2010 年 8 月 12 日学 科： 车辆工程 学 号： 0908020287 姓 名： 杨佳丽 成 绩： 汽车操纵稳定性仿真摘要：仿真技术日益广泛地应用于汽车工程领域，操纵稳定性研究越来越多地使用成熟的计算机仿真理论和高性能仿真软件进行分析研究，与辅助设计手段结合可直接指导和参与汽车设计参数的分析、优化与改进。本文在两自由度单轨操纵模型基础上，利用MATLAB对汽车的操纵动力学进行建模、仿真，并以别克轿车和法拉利跑车为例进行性能对比分

2、析，得出了对比图形和分析结果。关键词：MATLAB 操纵动力学 仿真Simulation of Handing Stability of Vehicle Based on MatlabAbstract: Simulation is widely used in vehicle engineering.Computer simulation theory and high-performance simulation software is more and more used for analysis and research in the study of Handing stability

3、 .And it also conbined with CAD are used in car design，such as parameters of analysis，optimization and improvement. In this paper，we can carry out Handing Dynamics modeling and carry on simulation on the basis of two degrees of freedom handing dynamic model.Then take buick and ferrari performance ca

4、rs for example for analysis and comparative，then obtained contrast graphics and analysis results.Key words: MATLAB Handing Dynamics modeling Simulation 1 前轮转向车辆的操纵性能计算机模型 汽车匀速行驶时，只有沿y轴的侧向运动与绕z轴的横摆运动两个自由度。此外，汽车的侧向加速度限定在0.4g以内，轮胎侧偏特性处于线性范围。图1-1所示为一个由前后两个侧向弹性的轮胎支撑于地面、具有侧向及横摆运动的线性二自由度汽车模型，其运动方程为： （1） （2

5、）式中、分别为汽车前、后轴的两个轮胎的侧向力合力。图1-1 前轮转向车辆的线性两自由度操纵模型单轨模型中，前轮的侧向速度为： （3）后轮的侧向速度为： （4）当很小时有tan=，则在后轮为非转向轴的情况下，后轮侧偏角可近似线性地表示为： （5）由于前轮产生一个转向角，且定义顺时针方向为正，得： （6）则前轮侧偏角近似为： （7）根据已知的轮胎的侧偏刚度与侧偏角的大小，可知前后轮侧向力分别为： （8） （9）将式（5）、式（7）、式（8）和式（9）代人式（1）、式（2）中，可得： (10)式中，m为汽车质量；v为汽车侧向速度；r为横摆角速度；为恒定的汽车前进速度；、分别为前后轮胎侧偏刚度。将系统

6、输入（即转向输入）整理到方程右边，并以状态空间方程的形式表示，则得到系统的运动方程为： （11）可将其转化为标准的状态空间方程形式：式中，；状态矢量；输入矢量。2 仿真结果对比分析表2-1给出的两辆不同车型的参数，可对其操纵稳定性进行仿真分析。本文选用别克轿车和法拉利跑车，仿真分析的内容包括：1）时域分析：分析角阶跃输入下横摆角速度的时域响应。2）频域分析：分析不同频率正弦波输入下横摆角速度的频域响应。3）稳定性分析：作出根轨迹图，分析不同车速下系统特征值的变化。表2-1 别克轿车和法拉利跑车模型参数参数符号单位别克法拉利质量横摆转动惯量前轴到质心距离后轴到质心距离前轮侧偏刚度后轮侧偏刚度ma

7、bKgmmkN/radkN/rad204554281.4881.71277.8576.51100810311.2341.022117.44144.932.1 时域分析在MATLAB中，建立前轮转角输入信号时间序列。系统输入U为转向盘转角输入，而非前轮转角。仿真中假设车辆行驶速度和转向系统传动比都恒定，分别为40m/s和45。当车辆以恒定速度直线行驶时，驾驶员突然给转向盘一个很小的角阶跃输入，即： （12）在MATLAB仿真中，根据已给的状态方程矩阵A、B、C、D和系统输入U（即前轮转角）的时间序列信号，应用线性模拟函数语句，即lsim（A，B，C，D，delta，t）（delta变量存储的是车

8、轮转角信号序列），就可方便地对两种车型进行正弦信号输入仿真，得到的横摆角速度时域响应结果如图2所示。图2-1 角阶跃输入下的横摆角速度时域响应由图2-1可知，在同样的转向盘转角输入下，法拉利跑车的瞬时响应比别克轿车的要好，主要体现在较短的稳定时间、较小的超调量以及更好的阻尼特性等。2.2 频域分析一个线性系统传递函数的幅频和相频特性可用伯德（Bode）图来表达。在MATLAB环境下，可直接采用bode（）命令，便可得到伯德图，即：sys=ss（A，B，C，D）bode（sys）当车辆以恒定速度直线行驶时，驾驶员突然给转向盘一个很小的正弦输入，即： （13）则根据转向系传动比=45，得到相应的前

9、轮转角输入为： （14）两种车型的系统频域响应如图3所示。图2-2 横摆角速度频率响应比较（车速）由图2-2可知，在同一行驶车速下，法拉利跑车的响应带宽大于别克轿车的响应带宽，从而也说明法拉利跑车具有更好的频率响应特性。法拉利跑车的系统响应滞后要比别克车的响应滞后少，系统延迟较小。2.3 稳定性分析通常，随着车速的提高，车辆的行驶稳定性下降。对具有过多转向特性的车辆而言，当车速超过其极限车速时，系统将处于不稳定工况，即意味着在很小的干扰输入时，系统将产生很大的响应输出（如高速转向时车辆可能发生侧滑）。而车辆行驶的稳定性与其等效线性系统的特征值有密切的联系。在MATLAB环境中有多种方法可以求系

10、统的特征根（即极点），其中最直接的方法可调用eig（）命令来求一个矩阵的特征值。由控制理论可知，一个可观、可控的线性系统的就是状态矩阵A的特征值，因此，求解车辆在不同车速（1560m/s范围内）的特征根，同时在同一复平面中作图，可采用如下MATLAB程序段：U=15：5：60；for j=1：length(U)A=-(Cf+Cr)/(M*U(j) -(a*Cf-b*Cr)/(M*U(j)-U(j)-(a*Cf-b*Cr)/(Iz*U(j) -(a2*Cf+b2*Cr)/(Iz*U(j)；plot(real(eig(A)，imag(eig(A)，'bx')；hold on；end其中，命令语句real（）和imag（）分别用于系统特征值的实部和虚部求解。不同车速下的根轨迹变化如图2-3所示。图2-3不同车速系统根轨迹图由图4可以得出与时域和频域响应分析相同的结果：随着车速的增加，系统的特征值由复平面的左侧不断向虚轴靠近，因此系统越来越趋向于不稳定。对于两种车型的比较可见，它们的特征值随车速变化而变化的趋势也不同，由于法拉利跑车有相当较大的稳定裕度，因而其特征根位置与别克轿车相比更远离虚轴，而别克轿车的转向特性的稳定裕量较小且容易失稳。3 结语应用计算机建模和仿真分析是现代汽车设计不可或缺的手段。本文以别克轿车和法拉利跑车为例进行性能对比分析，所建模型