毕设基于MATLAB的双横臂独立悬架优化设计

1、 基于MATLAB的双横臂独立悬架优化设计 目 录1 绪论2 1.1 引言21.2 国内外对悬架设计的研究概述31.3 本课题研究内容32 双横臂悬架与转向梯形42.1 双横臂独立悬架结构及其应用42.2 双横臂独立悬架及其转向系统52.3 本章小结83 转向梯形断开点位置模型的建立与优化93.1 转向梯形位置模型的建立93.2 转向梯形断开点模型的优化113.3 本章小结144 总结151 绪论1.2 国内外对悬架设计的研究1.3 本课题研究内容本课题主要分析基于matlab条件下对车辆双横臂独立悬架的优化设计，为了综合全面的分析问题，将有关空间运动学的原理和本课题相结合，得到最可行的优化结

2、果。1）分析双横臂式独立悬架其转向梯形在选择不同的位置时对汽车前轮跑偏、振动的影响，通过空间运动学的基本原理，画出各点的运动轨迹，以此来确定转向梯形的位置坐标。2）通过分析与设计，将matlab优化工具包运用到已建立起来的数学模型中，进而来判断建立起来的数学模型是否具有可行性。3）通过对悬架导向机构的优化分析，找到转向梯形的最佳断开点的位置坐标，最小化车轮绕主销的摆动量，降低转向杆系与悬架导向机构之间的不协调性误差。2 双横臂悬架与转向梯形 2.1 双横臂独立悬架结构及其应用双横臂独立悬架有不同的摆臂长度，可以相等，也可以不相等。如图2-1所示，若摆臂长度相等，则出现汽车上下振荡的时候，车轮所

3、在的水平面没有变化，但是两个车轮之间的距离会改变，车轮会向侧向移动。而如果摆臂长度不相同，且处在合适的位置，则上述等长情况下发生的变化则不会太大3。轮胎变形可以在轮距变化不大的情况下轻松地适应，可以允许轮距的改变在小于5mm的情况下不会使车轮出现左右滑动。 图2-1 不等长的双横臂式独立悬架可以在选择适当的参数下，在微小变形的情况下，对汽车的影响微乎其微，这样就可以减小汽车因为路面的不平度而带来的上下振荡。但是双横臂式独立悬架汽车也有一些缺点，比如悬架导向机构及其复杂，悬架制造成本高，悬架占用空间大等。但是因为其在安全性与舒适性方面有着无法比拟的优势，因此自上世纪以来不等长的双横臂式独立悬架便

4、广泛地应用在中小型汽车及微型货车的前轮上。2.2 双横臂独立悬架及转向系统 如图2-2所示为典型的不等长双横臂式独立悬架的构造。 图2-2图示的轿车属于常见的无主销式汽车，但是它设计了一个圆形状的结构代替主销，汽车在发生转弯时，车轮即绕该圆形状物体连心线偏转。路面对于车轮有垂直方向的作用力，然后再通过汽车悬架及其一系列导向机构传递到车身上，所以这就需要汽车的悬架要有强大的刚度，能抵御路面对汽车施加的强大反作用力。汽车转向系统还包括转向器。目前汽车上广泛使用齿轮齿条式转向器。其构造比较简易，且具有较大的刚度，加工简单，简化了转向传动机构，特别适合装配有双横臂式独立悬架的轻型汽车使用。轿车的齿轮齿

5、条式转向器，该转向器配备有较长的横拉杆，这样就可以在路面不平引起汽车震荡的时候4，可以减小因为汽车震荡而引起的横拉杆摆角，避免了双横臂式独立悬架系统与汽车的转向系统之间的运动不协调度，而且该类型的车一般都配有减振器，这样更能减小因为路面不平引起的汽车的上下震动。如下图2-3所示。 图2-3导向机构是导向装置的重要组成部分，对汽车的运动有重要的影响，因此对其设计有严格的要求：1）严格控制轮距的变化，保证不能超过行业规定值，否则对轮胎造成损害。2）要根据特定汽车参数选取标准的汽车前轮定位参数3）保证足够小的车身侧倾角，且车轮与车身倾斜相同。导向机构对汽车的影响主要集中于其对汽车转向机构的影响，由于

6、汽车在转动的时候要保证其前后轴的轴线要交于一点，特别是汽车如果在转向时速度很高，这时汽车的内外侧原件所受到的力度不一样，会有一个明显的倾斜，因为这个倾斜角度对汽车的安全稳定行驶有重要的影响，所以要设法减小这个内外侧倾斜度，这时一般汽车都装配了横向稳定杆，它能平衡汽车内外侧因为转向时所受到的受力不均匀，可以显著提高汽车的安全性与稳定性5。一般汽车为了汽车转向时，转向节与转向轮有相同的角度，大多数都装配有转向梯形机构，我们所要做的就是分析汽车转向机构的空间运动状态，做出其空间运动的方程与图像，找到能保持汽车转向机构运动协调性的最佳转向梯形断开点位置，如下图2-5所示的为汽车转向机构的模型。 图2-

7、4 图中为上摆臂摆动中心，为下摆臂摆动中心，为上摆臂摆动轴线 ，为下摆臂摆动轴线，G为下摆臂球销中心，A为上摆臂球销中心，D为主销与转向节轴线交点，E为车轮中心，F为车轮与地面交点，C为转向节臂球销中心，B为横拉杆断开点。分析转向梯形的空间模型，对汽车悬架的优化有重要影响。不同的悬架有不同的结构与功用。断开时转向梯形的横拉杆是断开的。无论是整体式还是断开式转向梯形，都要确保正确的参数，选择正确的转向梯形断开点位置，降低汽车悬架的运动不协调性6。一般独立悬架都用断开时转向梯形，非独立悬架倾向于用整体式转向梯形。断开时转向梯形能把汽车两侧车轮的运动状态分开，使这两个车轮的运动能够彼此独立，不互相影

8、响。但是断开时转向梯形由于转向杆系球头增多，所以结构很复杂，制造成本也很高，而且在调整前束也比较困难。如下图2-6即为断开时转向梯形的结构图。 图2-5 2.3 本章小结本章首先对汽车双横臂独立悬架的结构进行了介绍，并分析了汽车前轮独立悬架进行设计时候的要求。然后对双横臂式独立悬架汽车的转向系统进行了介绍，分析转向系统的结构以及一些常见的转向系统，并对一些转向系统的设计进行了简单的说明。最后再分析齿轮齿条式转向器，根据其在转向梯形位置中的不同布局，分析其作用。并对断开时转向梯形进行了介绍，分析断开时转向梯形在独立悬架前悬架中的作用，为下一章以转向梯形断开点的位置坐标为优化变量建立转向梯形断开点

9、位置优化模型打下了基础。3 转向梯形断开点位置模型的建立和优化3.1 转向梯形断开点位置模型建立汽车的转向系统、悬架系统与车轮定位之间有密不可分的关系，为了获得舒适安全的驾驶体验，且减小汽车行驶过程中对汽车的磨耗，有必要探讨一下汽车行驶过程中的角度问题，尤其应该重点关注汽车四个车轮之间的角度问题，找到角度之间最理想的数学关系式。众所周知，汽车在转弯时由于路面对汽车的行驶阻力加大，这样会加快路面对轮胎的磨损，危害汽车的行驶安全。转向梯形的设计就是为了减小这种磨损，它要求汽车在转向时能围绕一个瞬时转动中心行驶，保证所有的车轮维持纯滚动运动状态，减少滑动摩擦对汽车轮胎的磨损。显然汽车要想实现这样的运

10、动状态，这对汽车的结构有严格的要求，如图3-1所示，汽车要保证所有车轮在转向时绕着o点转动，o点称为转向中心，要想实现图中的要求，汽车应该满足阿克曼几何学7的要求。1）Ackerman 理论汽车在发生转弯时，想要保证纯滚动的运动状态，则要保证每个车轮都符合自然运动轨迹，即全部车轮绕同一瞬时转向中心o运动，则应有下式成立： (3-1)式中 外侧转向轮转角；内侧转向轮转角；B汽车两端主销轴线和地面交点的距离；L汽车前后轴距，R转弯半径。 图3-1 理想左右轮转向示意图理论上右轮转角可表示为 (3-2) 当汽车发生右轮转向时，转角关系为 (3-3)可得理想的右轮转角，如式所示： (3-4) 2)断开

11、点位置模型的建立转向梯形断开点的最佳位置应使图2-5中点C至点B的空间距离在车轮跳动过程中的变化量最小，并在这个基础上得到优化设计变量、约束条件和目标函数8。梯形断开点坐标： (3-5)约束条件 梯形断开点要受到汽车的尺寸以及汽车传动机构位置的要求，因而得到它所需要满足的条件： ， (3-6) 目标函数 由汽车在不平路面上行驶的时候，悬架会有伸缩的弹性变形，以缓解地面的冲击，这时以汽车轮胎为参照物的话，会上下跳动，这时轮胎就会有个极限跳动量，由可以得到下摆臂的摆角区间为，在这个区间内存在无数个摆动值（i=1,2,3），在这个基础上确定C点的坐标，即 (3-7)目标函数的表达式为 (3-8)式中

12、为稳定位置时悬架上C点到B点的空间距离，为下横臂转动的角度时候C点到B点的距离。由以上可以看到，我们所要完成的任务就是找到转向梯形断开点的最佳位置坐标，具体的优化程序已附在附件中。3.2 转向梯形断开点模型的优化将上述所建立的优化模型适用于某小型汽车上，运用MATLAB的优化工具箱优化梯形的断开点参数并绘图比较优化前后的转向特性图，进行相关的优化比较，并探讨可行性。该汽车的一些有关悬架的数据如下表1所示。表1 某汽车独立悬架配置参数9前轴载荷（满/空）主销内倾角主销后倾角车轮外倾角前轮距505/410kg8.810.51400mm前束260703452159030017835-3.5645-3

13、009600-42-136600-286刚开始可以利用平面作图的方法把转向梯形断开点的位置坐标确定下来，即 (3-9)可以找到四种不一样的摆臂轴线的投影角（其中X-Y可以定为，），因为要与平面作图法相比较，所以有必要使上述四种方法中的上下摆臂在Y-Z方向的平面内都有不变的投影长度和角度，也就是说使与在Y与Z方向上的坐标不能发生变化。通过查找翻阅资料，我们找到某轻型轿车车轮定位角的相关数据，具体数据可如下表2所示。 表2 四种方案比较方案/（）/（）/（）/（）/mm/mm1000-1.530202.85.8-2.817-5302.85.8017642.82.88.83.81922方案/mm/m

14、m/mm/mm/mm1250-2604.2092-18215-2525.542.23-38215-2465.4842.44-86215-2785.4642.6注意：指的是 的运动不协调误差；指的是 的运动不协调误差。由表2中的四种X,Y,Z的数值，可以得到表2中四种不同的双横臂独立悬架汽车运动的不协调误差值，根据此协调误差值对表格中1、2、3、4方案中的转向梯形断开点进行优化设计，这样便可以得到上述四种方案中的优化后的断开点位置坐标，由此可得出相应的运动不协调度误差。根据表2中的数据前后对比，可以看出早期使用的平面作图方法分析设计转向梯形断开点的结构不能满足汽车悬架优化设计对精确度的要求，因为

15、其忽略了空间角度的影响，所以会和实际结果有较大的偏差。所以根据空间结构运动学原理的优化设计方法可以较好地优化梯形断开点的位置，可以很好地避免因为忽略空间角度而产生的运动不协调误差10。 在上述基础上为了研究X方向的坐标对汽车运动不协调的影响，可以改变X方向的坐标而维持Y,Z方向上坐标不变，然后依次运行上述四种方案。将上述的优化参数输入程序中，运用matlab软件进行优化设计，相关的程序已经附在附件中，运行程序可以得到图3-2所示的的数据。运动不协调误差（）/mm605040302010 0 -200 -100 0 100 20012断a 断开点横坐标（）/mm 5040302010 03-20

16、0 -100 0 100 200运动不协调误差（）/mm4 b 断开点横坐标（）/mm图3-2 关系曲线根据图示可以看出来空间角度的存在与否对的值有很大的影响。这样也可以看出来平面作图法在摆臂轴线存在空间角度的情况下，有很大的不足。尤其当取值不恰当时，这时汽车内外轮的运动不协调度会发生更大的偏差。只有在空间角度不存在的情况下，转向梯形断开点在X方向的坐标值变化才不会对值有很大的影响。这时为了便于对优化结果的分析，我们可以分两种情况来讨论，研究空间角度的存在与否对运动不协调性的影响，并将这两种结果进行对比。可以得到如下图3-3所得的结果：运动不协调误差（）/mm-8 -6 -4 -2 0优化前曲

17、线优化后曲线36353433323130半轮距变化量（L）/mm 图3-3 转向梯形断开点优化前后运动不协调性比较 根据图形分析可以看出来，经过matlab的优化后，转向梯形断开点的位置相对于优化前对轮胎跳动性的影响发生了很大的变化，此时如果轮胎发生跳动将几乎不能够使转向横拉杆发生很大的变化。由此可以确定，优化后的转向梯形设计比优化前的转向梯形性能有很大的变化，这样就可以达到改善汽车运动，减小汽车轮胎磨损的目的，优化效果很明显11。3.3 本章小结本章主要是根据空间结构运动学原理，建立双横臂式独立悬架转向梯形断开点的模型，再运用matlab工具箱对建立的模型进行优化设计。通过优化转向梯形的最佳

18、断开点位置，可以得到最佳断开点的位置明显能改变汽车的转向特性。同时再分析具体实例，给出优化的可行性。 4总结 汽车的悬架系统是一个复杂的系统，对整个汽车的安全舒适行驶有着最重要的影响。双横臂式独立悬架由于其特殊的机构，对对改善汽车的平顺性，降低汽车的振动频率及汽车的安全行驶，起到了重要的作用，越来越多的适用于轿车和轻型货车上，因此对双横臂式独立悬架的优化设计就有着十分重要的现实意义。 本设计在汽车悬架优化设计方面主要有以下可取的地方。 （1）本设计并没有使用早期经常使用的平面作图分析方法，而是根据双横臂式独立悬架传动机构的空间运动轨迹，运用空间运动学分析法，优化转向梯形断开点的位置坐标，并进行

19、实际验证，这种方法能很好的反映转向梯形传动机构的运动特性。 （2）通过运用MATLAB软件进行导向机构的优化设计，再通过模拟曲线，得出经过优化后的转向梯形传动机构的性能明显优于优化前的，进而验证了优化的可行性。 （3）本设计所建立起来的研究方法具有可延展性与广泛性，能通用于目前大多数轿车的双横臂式独立悬架的优化设计中，具有深远的现实意义。附录（1）转向梯形断开点坐标主程序：clearh=-10:1:10;G3=QQ1(-10),QQ1(-9),QQ1(-8),QQ1(-7),QQ1(-6),QQ1(-5),QQ1(-4),QQ1(-3),QQ1(-2),QQ1(-1),QQ1(0),QQ1(1

20、),QQ1(2),QQ1(3),QQ1(4),QQ1(5),QQ1(6),QQ1(7),QQ1(8),QQ1(9),QQ1(10);I=cc(-10),cc(-9),cc(-8),cc(-7),cc(-6),cc(-5),cc(-4),cc(-3),cc(-2),cc(-1),cc(0),cc(1),cc(2),cc(3),cc(4),cc(5),cc(6),cc(7),cc(8),cc(9),cc(10);Gn3=new(-10),new(-9),new(-8),new(-7),new(-6),new(-5),new(-4),new(-3),new(-2),new(-1),new(0),ne

21、w(1),new(2),new(3),new(4),new(5),new(6),new(7),new(8),new(9),new(10);In=ee(-10),ee(-9),ee(-8),ee(-7),ee(-6),ee(-5),ee(-4),ee(-3),ee(-2),ee(-1),ee(0),ee(1),ee(2),ee(3),ee(4),ee(5),ee(6),ee(7),ee(8),ee(9),ee(10);J=hh(-10),hh(-9),hh(-8),hh(-7),hh(-6),hh(-5),hh(-4),hh(-3),hh(-2),hh(-1),hh(0),hh(1),hh(2)

22、,hh(3),hh(4),hh(5),hh(6),hh(7),hh(8),hh(9),hh(10);Jn=hhn(-10),hhn(-9),hhn(-8),hhn(-7),hhn(-6),hhn(-5),hhn(-4),hhn(-3),hhn(-2),hhn(-1),hhn(0),hhn(1),hhn(2),hhn(3),hhn(4),hhn(5),hhn(6),hhn(7),hhn(8),hhn(9),hhn(10);subplot(1,2,1),plot(G3,I,-),hold on;plot(G3,In,-),box off,hold offxlabel(Deltaz/mm)ylabe

23、l(Deltatau /(circ)subplot(1,2,2),plot(G3,J,-),hold on;plot(G3,Jn,-),box off,hold offxlabel(Deltaz/mm)ylabel(Deltatheta /(circ)white(2)梯形断开点坐标调用程序：function G3=QQ1(h)u=h*pi/180;q0=cos(u/2);q1=1.*sin(u/2); %这几个值将来再进行输入q2=0.*sin(u/2); %这几个值将来再进行输入q3=0.*sin(u/2); %u相当于以前的阿尔法B0=-1.5;617;-324.8; %这几个值将来再进行

24、输入O1=-1.5;257;-281; %这几个值将来再进行输入B=Q1*(B0-O1)+O1;B1=B(1:1,1);B2=B(2:2,1);B3=B(3:3,1);R1=363; %这几个值将来再进行输入R2=252; %这几个值将来再进行输入R3=220; %这几个值将来再进行输入R4=260; %这几个值将来再进行输入t=4.31;M=(R3.2-(R1.2+R2.2+R4.2+2*R1.*R4.*cos(u+6.1588-t)/(2*R2.*(R1.*cos(u+6.1588)+R4.*cos(t);%w相当于以前的噶玛，t相当于以前的依浦西龙N=(R1.*sin(u+6.1588)

25、+R4.*sin(t)/(R1.*cos(u+6.1588)+R4.*cos(t);w=asin(M.*N+sqrt(N.2-M.2+1)/(1+N.2);c=asin(R1.*sin(u+6.1588)+R2.*sin(w)+R4.*sin(t)/R3);v=c+0.1629; p0=cos(v/2);p1=1.*sin(v/2); %这几个值将来再进行输入p2=0.*sin(v/2); %这几个值将来再进行输入p3=0.*sin(v/2); %v相当于以前的贝塔D0=6.9;578;-76; %这几个值将来再进行输入O2=6.9;360.7;-43; D=Q2*(D0-O2)+O2;D1=

26、D(1:1,1);D2=D(2:2,1);D3=D(3:3,1);a=0.5975;E=(1/(1+a)*B+(a/(1+a)*D; %a相当于以前的蓝姆达E1=E(1:1,1);E2=E(2:2,1);E3=E(3:3,1);v1=asin(D1-B1)/(R2.*cos(0.1798); %计算贝塔一v2=asin(B2-D2)/(R2.*cos(0.1798).*sqrt(1+cos(v1).2)-atan(tan(0.1798)/cos(v1);%计算贝塔二H=(-0.0665-v2)/pi*180; %前轮外倾角F=E+0;95.*cos(v2);95.*sin(v2);F1=F(1

27、:1,1);F2=F(2:2,1);F3=F(3:3,1);I=(atan(B1-D1)/(B3-D3)-0.0337)/pi*180; %主销后倾角G=F+0;334.*sin(v2);-334.*cos(v2)-1.6000;675.1321; -573.5527;G1=G(1:1,1);G2=G(2:2,1); %前轮半轮距变化量G3=G(3:3,1) %轮跳变化量 K=1/2.*B+1/2.*D;K1=K(1:1,1);K2=K(2:2,1);K3=K(3:3,1);J=360.*(atan(-(K1-F1)/(F2-K2)/pi*1340+1.6908e+003; %前轮前束L=(a

28、tan(B2-D2)/(D3-B3)-0.1555)/pi*180;%主销内倾角（3）生成图3-2程序：function wwh=-10:1:10;J=func(-10),func(-9),func(-8),func(-7),func(-6),func(-5),func(-4),func(-3),func(-2),func(-1),func(0),func(1),func(2),func(3),func(4),func(5),func(6),func(7),func(8),func(9),func(10);L=ff1(-10),ff1(-9),ff1(-8),ff1(-7),ff1(-6),f

29、f1(-5),ff1(-4),ff1(-3),ff1(-2),ff1(-1),ff1(0),ff1(1),ff1(2),ff1(3),ff1(4),ff1(5),ff1(6),ff1(7),ff1(8),ff1(9),ff1(10);subplot(1,2,1),plot(h,J,r*),axis(-10,10,-5,20),box off,title(子图(1)subplot(1,2,2),plot(h,L,g-),axis(-10,10,-2,2),title(子图(2)（4）生成图3-3程序：function G3=yhlt(h)u=h*pi/180;q0=cos(u/2);q1=1.*

30、sin(u/2); %这几个值将来再进行输入q2=0.*sin(u/2); %这几个值将来再进行输入q3=0.*sin(u/2); %u相当于以前的阿尔法B0=-1.62;617.30;-322; %这几个值将来再进行输入O1=-1.62;257.0;-281.0; %这几个值将来再进行输入B=Q1*(B0-O1)+O1;B1=B(1:1,1);B2=B(2:2,1);B3=B(3:3,1);R1=363; %这几个值将来再进行输入R2=252; %这几个值将来再进行输入R3=220; %这几个值将来再进行输入R4=260; %这几个值将来再进行输入t=16.1*pi/180;M=(R3.2-(R1.2+R2.2+R4.2+2*R1.*R4.*cos(u-t)/(2*R