汽车转向梯形优化设计Word文档下载推荐.docx
《汽车转向梯形优化设计Word文档下载推荐.docx》由会员分享,可在线阅读,更多相关《汽车转向梯形优化设计Word文档下载推荐.docx(16页珍藏版)》请在冰豆网上搜索。
转向梯形机构用来保证汽车转弯行驶时所有车轮能绕一个瞬时转向中心,在不同得圆周上做无滑动得纯滚动。
设计转向梯形得主要任务之一就是确定转向梯型得最佳参数与进行强度计算。
一般转向梯形机构布置在前轴之后,但当发动机位置很低或前轴驱动时,也有位于前轴之前得。
转向梯形有整体式与断开式两种,选择整体式或断开式转向梯形方案与悬架采用何种方案有联系。
无论采用哪一种方案,必须正确选择转向梯形参数,做到汽车转弯时,保证全部车轮绕一个瞬时转向中心行驶,使在不同圆周上运动得车轮,作无滑动得纯滚动运动。
同时,为达到总体布置要求得最小转弯直径值,转向轮应有足够大得转角。
二、整体式转向梯形结构方案分析
图5、1整体式转向梯形
1—转向横拉杆2—转向梯形臂3-前轴
整体式转向梯形就是由转向横拉杆1,转向梯形臂2与汽车前轴3组成,如图5、1所示。
其中梯形臂呈收缩状向后延伸。
这种方案得优点就是结构简单,调整前束容易,制造成本低;
主要缺点就是一侧转向轮上、下跳动时,会影响另一侧转向轮。
当汽车前悬架采用非独立悬架时,应当采用整体式转向梯形。
整体式转向梯形得横拉杆可位于前轴后或前轴前(称为前置梯形).对于发动机位置低或前轮驱动汽车,常采用前置梯形。
前置梯形得梯形臂必须向前外侧方向延伸,因而会与车轮或制动底板发生干涉,所以在布置上有困难.为了保护横拉杆免遭路面不平物得损伤,横拉杆得位置应尽可能布置得高些,至少不低于前轴高度。
三、整体式转向梯形机构优化分析
汽车转向行驶时,受弹性轮胎侧偏角得影响,所有车轮不就是绕位于后轴沿长线上得点滚动,而就是绕位于前轴与后轴之间得汽车内侧某一点滚动。
此点位置与前轮与后轮得侧偏角大小有关.因影响轮胎侧偏角得因素很多,且难以精确确定,故下面就是在忽略侧偏角影响得条件下,分析有关两轴汽车得转向问题.此时,两转向前轮轴线得延长线应交在后轴延长线上,如图5—2所示。
设θi、θo分别为内、外转向车轮转角,L为汽车轴距,K为两主销中心线延长线到地面交点之间得距离。
若要保证全部车轮绕一个瞬时转向中心行驶,则梯形机构应保证内、外转向车轮得转角有如下关系:
(1)
图1 理想得内、外车轮转角关系简图
若自变角为θo,则因变角θi得期望值为:
(2)
现有转向梯形机构仅能近似满足上式关系.以图所示得后置梯形机构为例,在图上作辅助用虚线,利用余弦定理可推得转向梯形所给出得实际因变角为
(3)
式中:
m为梯形臂长;
γ为梯形底角。
所设计得转向梯形给出得实际因变角,应尽可能接近理论上得期望值.其偏差在最常使用得中间位置附近小角范围内应尽量小,以减少高速行驶时轮胎得磨损;
而在不经常使用且车速较低得最大转角时,可适当放宽要求。
因此,再引入加权因子,构成评价设计优劣得目标函数为
(4)
由以上可得:
(5)
式中:
x为设计变量,;
θomax为外转向车轮最大转角,由图2得 (6)
式中,Dmin为汽车最小转弯直径;
a为主销偏移距。
考虑到多数使用工况下转角θo小于20°
,且10°
以内得小转角使用得更加频繁,因此取:
(7)
建立约束条件时应考虑到:
设计变量m及γ过小时,会使横拉杆上得转向力过大;
当m过大时,将使梯形布置困难,故对m得上、下限及对γ得下限应设置约束条件。
因γ越大,梯形越接近矩形,值就越大,而优化过程就是求得极小值,故可不必对γ得上限加以限制。
综上所述,各设计变量得取值范围构成得约束条件为:
(8)
梯形臂长度m设计时常取在mmin=0、11K,mmax=0、15K。
梯形底角γmin=70°
此外,由机械原理得知,四连杆机构得传动角δ不宜过小,通常取δ≥δmin=40°
。
如图5-2所示,转向梯形机构在汽车向右转弯至极限位置时达到最小值,故只考虑右转弯时δ≥δmin即可。
利用该图所作得辅助用虚线及余弦定理,可推出最小传动角约束条件为:
(9)
δmin为最小传动角。
δmin=40°
故由式可知,δmin为设计变量m及γ得函数。
由式(6)、式(7)、式(8)与式(9)四项约束条件所形成得可行域,如图3所示得几种情况。
图3b适用于要求δmin较大,而γmin可小些得车型;
图5-3c适用于要求γmin较大,而δmin小些得车型;
图3a适用介于图3b、c之间要求得车型。
图3 转向梯形机构优化设计得可行域
四、整体式转向梯形程序编写
(1)优化编程所需数据:
轴距:
L=2775mm
轮距:
K=1560mm
最小转弯半径:
R=5300mm
转向梯形臂:
m
计算可得底边长:
L—2*a
(2)function fuun。
m编辑过程
在MATLAB窗口新建一个空白M文件
将下式输入
functionc=theatar()
%建立主函数
global optionsLb r aKthetamaxclcrfi0%定义全局变量
K=1638;
%input('
输入主销中心线间距(mm)’);
%依次给予几个变量赋值
L=3308;
%input('输入轴距(mm)');
thetamax=40;
%input('输入外转向轮最大转角(度)’);
x
(1)=175;
%input(’臂长(mm)');
x(2)=74、5;
%input('
底角(度)'
);
b=8;
%input('
内倾角(度)'
);
r=2;
%input('后倾角(度)'
a=1;
%input('
外倾角(度)’);
thetamax=thetamax*pi/180;
%单位转换,弧度与度数转变
lb
(1)=0、11*K;
%设置上下限
lb(2)=1、2217;
%acot(K/(1、2*L));
ub
(1)=0、13*K;
ub
(2)=pi/2;
fil=linspace(0,thetamax,61);
lb=[lb
(1),lb(2)];
ub=[ub
(1),ub
(2)];
x0=[x
(1),x
(2)];
% A=[0、2510、372];
%b=[0、143];
[y,fval]=fmincon(’fuun'
x0,[],[],[],[],lb,ub,[]);
%利用工具箱中得x= fmincon(fun,x0,A,b,Aeq,beq,lb,ub,nonlcon)进行计算
Y=y;
%[y,resnorm]=lsqnonlin(’fuun'
,x0,lb,ub,options) %betae(i) y= fmincon(fuun,x0,[],[],[],[],lb,ub)
fori=1:
61 %设置60个区域
fil=linspace(0,thetamax,61);
%betae(i)=acot(cot(fil(i))—(K/L));
fi=fii(r*pi/180,b*pi/180);
%以下将各公式单位转换,并代入公式
dt=delta(r*pi/180,fii(r*pi/180,b*pi/180));
%=dt d=Di(fii(r*pi/180,b*pi/180),a*pi/180,delta(r*pi/180,fii(r*pi/180,b*pi/180)));
%=d
Mid_w=Ww(a*pi/180,delta(r*pi/180,fii(r*pi/180,b*pi/180)),Di(fii(r*pi/180,b*pi/180),a*pi/180,delta(r*pi/180,fii(r*pi/180,b*pi/180))));
%=w
a1(i)=alfa(fi,d,fil(i),Mid_w);
A(i)=K*cos(b*pi/180)—x(1)*cos(2*(b*pi/180))*cos(x
(2)*pi/180+fil(i));
B(i)=x
(1)*sin(x
(2)*pi/180+fil(i));
C(i)=K*cos(b*pi/180)*cos(x
(2)*pi/180+fil(i))-2*K*cos(b*pi/180)*cos(x
(2)*pi/180)+2*x(1)*(cos(b*pi/180)^2)*(cos(x
(2)*pi/180)^2)—x
(1);
fir(i)=abs(fiir(A(i),B(i),C(i),x(2)*pi/180));
a2(i)=alfa2(fi,d,fir(i),Mid_w);
cl(i)=acos((cos(lamta(d,fil(i)))-sin(a*pi/180)*sin(a1(i)))/(cos(a*pi/180)*cos(a1(i))))*180/pi;
cr(i)=acos((cos(lamta(d,fir(i)))-sin(a*pi/180)*sin(a2(i)))/(cos(a*pi/180)*cos(a2(i))))*180/pi;
betae(i)=acot(cot(cl(i)*pi/180)-(K/L));
Aa(i)=K*cos(b*pi/180)—y
(1)*cos(2*(b*pi/180))*cos(y(2)+fil(i));
Ba(i)=y
(1)*sin(y
(2)+fil(i));
Ca(i)=K*cos(b*pi/180)*cos(y
(2)+fil(i))-2*K*cos(b*pi/180)*cos(y
(2))+2*y
(1)*(cos(b*pi/180)^2)*(cos(y
(2))^2)-y
(1);
fira(i)=abs(fiir(Aa(i),Ba(i),Ca(i),y
(2)));
a2a(i)=alfa2(fi,d,fira(i),Mid_w);
% cl(i)=acos((cos(lamta(d,fil(i)))-sin(a*pi/180)*sin(a1(i)))/(cos(a*pi/180)*cos(a1(i))))*180/pi;
cr1(i)=acos((cos(lamta(d,fira(i)))-sin(a*pi/180)*sin(a2a(i)))/(cos(a*pi/180)*cos(a2a(i))))*180/pi;
end
plot(cl,betae*180/pi,'
r',cl,cr,'
b'
cl,cr1,’-—g'
%plot(cl,cr1,’—-b’);
axis([0,40,0,45]);
xlabel('外转向轮输入角(单位:
度)’);
ylabel('
内转向轮输出角(单位:
度)’);