整体式转向梯形优化设计MATLAB程序

kk123

整体式转向梯形优化设计MATLAB程序
第一个m文件 objfun22.m
function f=objfun22(x)L=4430;M=2150;M1=2360;gam=x(1)*pi/180;m=x(2);for i=1:35;    afa(i)=i*pi/180;    N=sqrt(M^2+m^2-2*M*m*sin(gam+afa(i)));    deta1=asin(m*cos(gam+afa(i))/N);    b=M-2*m*sin(gam);    beta=gam+deta1-asin((N^2+m^2-b^2)/(2*N*m));    beta0=atan(1/(M1/L+1/tan(afa(i))));    c(i)=abs(beta-beta0)*180/pi;endf=max(c);12345678910111213141516
解释：
L为轴距
M 为两转向主销的轴线延长线与地面交点之间的距离
afa(i) 外侧转角
beta 实际转角
beta0 理论转角

第二个m文件 confun22.m
function[c,ceq]=confun22(x)L=4430;M=2150;M1=2360;gam=x(1)*pi/180;m=x(2);afa35=35*pi/180;N=sqrt(M^2+m^2-2*M*m*sin(gam+afa35)); deta1=asin(m*cos(gam+afa35)/N);    b=M-2*m*sin(gam);    beta35=gam+deta1-asin((N^2+m^2-b^2)/(2*N*m));    dlt=180-acos((b^2-M^2+2*M*x(2)*sin(gam-beta35))/(2*b*x(2)))*180/pi;    c=[x(1)-40;(90-atan(1.2*L/M)*180/pi)-x(1);x(2)-0.22*M;0.11*M-x(2);30-dlt];    ceq=[];1234567891011121314
解释：
afa35 外侧转角35度
beta35 实际转角
dlt 外转角35度时的传动角
主程序
clear allx0=[15,210];options=optimset('largescale','off');[x,fval]=fmincon('objfun22',x0,[],[],[],[],[],[],'confun22',options)L=4430;M=2150;M1=2360;gam=x(1)*pi/180;%将足球比分yaoji1真钱，扎金花角转化为弧度m=x(2);for i=1:35;     afa(i)=i*pi/180;%外侧转角    N=sqrt(M^2+m^2-2*M*m*sin(gam+afa(i)));    deta1=asin(m*cos(gam+afa(i))/N);    b=M-2*m*sin(gam);    beta=gam+deta1-asin((N^2+m^2-b^2)/(2*N*m));%实际内转角    beta0=atan(1/(M1/L+1/tan(afa(i))));%理论内转角    c(i)=abs(beta-beta0)*180/pi;   bt0(i,:)=beta0.*180/pi;%建立理论外侧转角数组    bt(i,:)=beta.*180/pi;%建立实际外侧转角数组   ai(i,:)=afa(i).*180/pi;endplot(ai,bt0,ai,bt)%画出实际与理论转角曲线figure(2)plot(ai,c)123456789101112131415161718192021222324
按照这个给定条件得到结果为：
x =

21.9551 236.4940

转向梯形角为 21.9551度 梯形臂长度236.4940mm
fval =

0.5773外侧理论与实际转角误差

注意：直接运行主程序




外侧车轮转角




内测理论与实际转角偏差