汽车半悬挂系统建模与分析现代控制理论大作业.pdf
《汽车半悬挂系统建模与分析现代控制理论大作业.pdf》由会员分享,可在线阅读,更多相关《汽车半悬挂系统建模与分析现代控制理论大作业.pdf(19页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、 XX 大学 现代控制理论 汽车半主动悬架系统的建模与分析 姓名:XXX 学号:XXXX 专业:XXXX 一 课题背景 汽车的振动控制是汽车设计的一个重要研究内容,涉及到汽车的平顺性和操纵稳定性。悬架系统是汽车振动系统的一个重要子系统,其振动传递特性对汽车性能有很大影响。因此设计性能良好的悬架系统以减少路面激励的振动传递,从而提高汽车的平顺性和操纵稳定性是汽车振动控制研究的重要课题。悬架系统是汽车车身与轮胎间的弹簧和避震器组成整个支撑系统,用于支撑车身,改善乘坐舒适度。而半主动悬架是悬架弹性元件的刚度和减振器的阻尼系数之一可以根据需要进行调节控制的悬架。目前,半主动悬架研究主要集中在调节减振器
2、的阻尼系数方面,即将阻尼可调减振器作为执行机构,通过传感器检测到汽车行驶状况和道路条件的变化以及车身的加速度,由 ECU根据控制策略发出脉冲控制信号实现对减振器阻尼系数的有级可调和无级可调。二 系统建模与分析 1.1 半主动悬架系统的力学模型 以二自由度 1/4 半主动悬架模型为例,并对系统作如下假设:(1)悬挂质量与非悬挂质量均为刚体;(2)悬架系统具有线性刚度和阻尼;(3)悬架在工作过程中不与缓冲块碰撞;(4)轮胎具有线性刚度,且在汽车行驶过程中始终与地面接触。综上,我们将该系统等效为两个质量块 M,m;两个弹簧系统 Ks,Kt;一个可调阻尼器(包含一个常规阻尼器 Cs 和一个变化阻尼力
3、F),如图 1 所示。图 1 系统力学模型 1.2 半主动悬架系统的数学模型 由减振器的简化模型得:NS FC VF 对 m 进行分析:211201122()tssd zdzdzmK zzKzzCFdtdtdt 即:1011212()tssmzK zzKzzCzzF 对 M 进行分析:2212122()ssd zdzdzMKzzCFdtdtdt 即:21212ssMzKzzCzzF 选取状态变量:1102213142xzzxzzxzxz,输入变量:uF 输出变量:1122yxyx,综上可得,系统状态空间表达式为:11032214331234423411tssssssxzzxxzzxxKKCCx
4、xxxxFmmmmmKCCxxxxFMMMM 整理得:0010000110110tssssssKKCCmmmmmKCCMMMMxxu 10000100yx 三 数值化分析 选取系统参数为:M=391 kg,m=50.7 kg,Ks=60KN/m,Kt=362 KN/m,Cs 取 1 KNs/m。状态空间表达式变为:001000011071401183.4319.7219.72-0.020-153.452.56-2.560.0026xxu 10000100yx 四 能控性与能观性分析 00100001101000,71401183.4319.7219.72-0.0201000-153.452.5
5、6-2.560.0026Abc 4.1 能控性分析 能控性矩阵:23(,)Mb Ab A b A b 通过 matlab 计算得:Rank(M)=4,满秩,故系统可控。4.2 能观性分析 能观性矩阵:23,TNC CA CACA 通过 matlab 计算得:Rank(N)=4,满秩,故系统可观。五 稳定性分析 存在唯一平衡点 x=0,对矩阵 A 进行特征值计算:通过 MATLAB 计算,我们得到特征值为:-10.2018+90.5683i,-10.2018-90.5683i,-0.9382+11.4463i,-0.9382-11.4463i。由于矩阵 A 的特征值均有负实部,所以系统是大范围渐
6、近稳定的。六 状态观测器设计 因为系统完全能观,所以可以设计状态观测器。6.1 全维观测器 将系统极点配置为:-1,-2,-3,-4.MATLAB 程序:A=0,0,1,0;0,0,-1,1;-7140,1183.43,-19.72,19.72;0,-153.45,2.56,-2.56;b=0;0;-0.02;0.0026;c=1,0,0,0;0,1,0,0;opt=-1,-2,-3,-4;G=(place(A,c,opt);输出结果为:7877.00131.08140.09822.03395.101621.13873.789405.1G 所以,全维观测器方程为:001001.940578.3
7、873001101.162110.3395()71401183.4319.7219.72-0.020.98220.81400-153.452.56-2.560.00260.01310.7877yyxxu+6.2 降维观测器 由于 rank(c)=2,n=4,所以将系统极点配置为-1,-2.构造变换阵作线性变换,设10 0 1 00 0 1 00 0 0 10 0 0 1,1 0 0 01 0 0 00 1 0 00 1 0 0TT。则,119.719.771401183.42.62.60153.410001100ATAT,10.020.00260010,000010BT BCCT。MATLAB
8、 程序:opt2=-1,-2;T=0,0,1,0;0,0,0,1;1,0,0,0;0,1,0,0;Tni=inv(T);A_2=Tni*A*T;B_2=Tni*B;C_2=C*T;A_11=A_2(1:2,1:2);A_21=A_2(3:4,1:2);G2=(place(A_11,A_21,opt2);输出结果为:56.0272.1912G。所以,降维观测器方程为:110.980.0271401183.40.020.042.040153.40.0026119.7220.56wxyuxwy 七 最优控制 对于半主动悬架系统,最优控制器的设计目的就是寻找最优控制 F,使实现控制所需的能量为最小:d
9、txxqxqJ024222211,其中,1q,2q分别为轮胎动变形加权系数,悬架动挠度加权系数,为车身加速度加权系数。将目标性能泛函改写成二次型性能指标形式:dtRuuQxxJTT0)(,这里,0000000000000021qqQ,为半正定常数矩阵;21MR,为正定常数矩阵。所以,最优控制存在,且唯一:)()(1tPxBRtuT 式中,P 为44维正定常数矩阵,满足黎卡提矩阵代数方程:01QPBPBRPAPATT 采用试探法取三组不同权系数1q、2q,运用 MATLAB 进行计算分析:(1)q1=3.35e5,q2=40.5e5;(2)q1=3.35e8,q2=40.5e8;(3)q1=3.
10、35e9,q2=40.5e9;Matlab 程序:%最优控制 clc;clear;M=391;A=0,0,1,0;0,0,-1,1;-7140,1183.43,-19.72,19.72;0,-153.45,2.56,-2.56;B=0;0;-0.02;0.0026;C=1,0,0,0;0,1,0,0;D=0;R=1/M2;%求不同 Q、R 下的状态反馈阵 K Q1=3.35e5;Q2=40.5e5;Q=Q1,0,0,0;0,Q2,0,0;0,0,0,0;0,0,0,0;K P e=lqr(A,B,Q,R)Ac=(A-B*K);Bc=B;Cc=C;Dc=D;T=0:0.05:5;U=0.2*on
11、es(size(T);Y,X1=lsim(Ac,Bc,Cc,Dc,U,T);Q1=3.35e8;Q2=40.5e8;Q=Q1,0,0,0;0,Q2,0,0;0,0,0,0;0,0,0,0;K P e=lqr(A,B,Q,R)Ac=(A-B*K);Bc=B;Cc=C;Dc=D;T=0:0.05:5;U=0.2*ones(size(T);Y,X2=lsim(Ac,Bc,Cc,Dc,U,T);Q1=3.35e9;Q2=40.5e9;Q=Q1,0,0,0;0,Q2,0,0;0,0,0,0;0,0,0,0;K P e=lqr(A,B,Q,R)Ac=(A-B*K);Bc=B;Cc=C;Dc=D;T=0:0
12、.05:5;U=0.2*ones(size(T);Y,X3=lsim(Ac,Bc,Cc,Dc,U,T);figure;hold on;plot(T,X1(:,1),-,color,black);plot(T,X2(:,1),-,color,green);plot(T,X3(:,1),-.,color,red);xlabel(时间(s));ylabel(轮胎动变形(m));hold off;legend(q1=3.35e5,q2=40.5e5,q1=3.35e8,q2=40.5e8,q1=3.35e9,q2=40.5e9);figure;hold on;plot(T,X1(:,2),-,colo
13、r,black);plot(T,X2(:,2),-,color,green);plot(T,X3(:,2),-.,color,red);xlabel(时间(s));ylabel(悬架动挠度(m));hold off;legend(q1=3.35e5,q2=40.5e5,q1=3.35e8,q2=40.5e8,q1=3.35e9,q2=40.5e9);figure;hold on;plot(T,X1(:,3),-,color,black);plot(T,X2(:,3),-,color,green);plot(T,X3(:,3),-.,color,red);xlabel(时间(s));ylabel
14、(悬架动载荷(N));hold off;legend(q1=3.35e5,q2=40.5e5,q1=3.35e8,q2=40.5e8,q1=3.35e9,q2=40.5e9);figure;hold on;plot(T,X1(:,4),-,color,black);plot(T,X2(:,4),-,color,green);plot(T,X3(:,4),-.,color,red);xlabel(时间(s));ylabel(车身加速度(m/s2));hold off;legend(q1=3.35e5,q2=40.5e5,q1=3.35e8,q2=40.5e8,q1=3.35e9,q2=40.5e
15、9);matlab 仿真结果如下:00.511.522.533.544.55-6-4-2024681012x 10-8时 间(s)轮胎动变形(m)q1=3.35e5,q2=40.5e5q1=3.35e8,q2=40.5e8q1=3.35e9,q2=40.5e9 图 2 轮胎动变形变化趋势 00.511.522.533.544.5500.511.522.53x 10-7时 间(s)悬架动挠度(m)q1=3.35e5,q2=40.5e5q1=3.35e8,q2=40.5e8q1=3.35e9,q2=40.5e9 图 3 悬架动挠度的变化趋势 00.511.522.533.544.55-3-2-10
16、12345x 10-6时 间(s)悬架动载荷(N)q1=3.35e5,q2=40.5e5q1=3.35e8,q2=40.5e8q1=3.35e9,q2=40.5e9 图 4 悬架动载荷的变化趋势 00.511.522.533.544.55-3-2-1012345x 10-6时 间(s)车身加速度(m/s2)q1=3.35e5,q2=40.5e5q1=3.35e8,q2=40.5e8q1=3.35e9,q2=40.5e9 图 5 车身加速度的变化趋势 通过 MATLAB 仿真得到,加权系数对悬架性能有较大的影响,当1q、2q取得较大值时,车身加速度,悬架动挠度及轮胎动变形的波动很小。当 q1=3
17、.35e9,q2=40.5e9 时,由图1,2 可以看出,悬架动挠度和轮胎动变形几乎为 0,可视为最优状态。八 总结 本次大作业主要完成了对汽车半主动悬架系统的建模与分析。在这次过程中,首先,建立系统状态空间表达式,然后对系统进行能观能控性及稳定性分析;其次,通过学习也对系统观测器进行设计,了解全维观测器和降维观测器的应用和区别;最后,对最优控制有了一定的了解,通过设置 Q 矩阵的参数,可以使系统最后的误差和过程中的能量损耗达到一个设计者预想的一个结果,即使汽车平顺性和操纵稳定性达到最优状态,此外由于加权系数的选取存在随机性,仿真结果仍存有误差,希望在今后的理论和实践学习中进一步完善和改进该模
18、型。参考文献 朱明.汽车半主动悬架系统的研究.重庆大学硕士学位论文.2004附录资料:不需要的可以自行删除 预应力锚索桩板墙施工工法 一、前言 在山岭陡峭、地形复杂、山高谷深的地区,高等级公路通过的地段造成大量的高填深挖,高桥及隧道处处可见。在山谷深、地面横坡陡峭的地段,路基难于填筑,旱桥跨越在经济和技术上造成较大的浪费,同时也给路基稳定及桥梁的桥桩、墩柱带来隐患。采用新型高挡墙跨越不仅开挖面小,也可消耗废方,起到安全、经济和环保的作用。个旧至冷墩二级公路预应力锚索桩板墙工程是采用 40 米高预应力锚索桩板墙进行边坡治理的项目,稳定了高填方路基,减少了陡坡旱桥,预应力锚索结构由于其合理的受力机
19、理以及在软弱岩体中能更有效的发挥土体承载力而提供了较大锚固力,通过施工经总结形成本工法。二、工法特点 1.采用 MG-50A 型潜孔冲击钻跟套管无水干钻,能有效的预防塌孔,保证水泥浆与孔壁岩体的粘结强度。2.锚索材料选用低松弛环氧喷涂无粘结钢绞线(ASTMT416-88a 标准 270 级,强度 Rby=1860Kpa,松弛率为 3.5%,j=15.24mm),配套 OVM15 型锚具,钢绞线强度高,性能好,可以在张拉结束后有效的进行放张或补偿张拉且弥补了钢绞线在特殊环境下中长久防腐的问题。3.该体系能主动提供抗滑力,有效的控制岩体的位移,在锚索的锚固范围内产生亚应力带,从而从根本上改善岩体的
20、力学性能。4.根据现场实际地质情况,大吨位锚索主要锚于碎石土、亚粘土中,鉴于土体破碎,抗剪强度低,在锚索结构上,通过对拉力型锚索与分散压缩型锚索工作性能的比较,采用分散压缩型锚索结构有突出优点。拉力型锚索与分散压缩型锚索工作性能的比较见表 1。表 1 拉力型锚索与分散压缩型锚索工作性能的比较 项 目 拉力型锚索 分散压缩型锚索 岩体水泥浆体间的粘结摩阻应力分布状况 沿锚固体长度分布极不均匀,应力集中严重,易发生渐进性破坏 沿锚固长度分布较均匀 岩体水泥浆体间的粘结摩阻应力值 总拉力大,粘结摩阻应力值大 总拉力可分散成几个较小的压力,粘结摩阻应力值显著减小 粘结摩阻强度 灌浆体受拉不会引起水泥浆
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 汽车 悬挂 系统 建模 分析 现代 控制 理论 作业
限制150内