资源描述：
《用魔术公式建立轮胎模型.pdf》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、汽车系统动力学作业1用魔术公式建立轮胎模型根据资料，通过给定的参数以及对一些比例因子的设定，用魔术公式建立如下轮胎模型。1.纯纵滑情况：%%LongitudinalForce(purelongitudinalslip)%%kisinputvalue,meansslipratio%%Fx0isoutputvalue,meanslongitudinalforcek=linspace(-1,1,300);dfz=(Fz-Fz0)/Fz0;C_x=P_Cx1;%取lam(Cx)=1shapefactoru_x=(P_Dx1+
2、P_Dx2*dfz)*0.9;%取lam(ux*)=0.9compositefrictionscalingfatorD_x=u_x*Fz;%zeta(1-8)都取1S_Hx=(P_Hx1+P_Hx2*dfz);%取lam(Hx)=1horizontalshiftk_x=k+S_Hx;E_x=(P_Ex1+P_Ex2*dfz+P_Ex3*dfz^2)*(1-P_Ex4*sign(k_x));%取lam(Ex)=1curvaturefactorK_xk=Fz*(P_Kx1+P_Kx2*dfz)*exp(P_Kx3*dfz
3、);%取lam(Kxk)=1breakslipstiffnessB_x=K_xk/(C_x*D_x+0.1);%取epsilon=0.1S_vx=Fz*(P_Vx1+P_Vx2*dfz)*(V0/(0.1+V0));%lam(Vx),lam(ux),zeta1=1,epsilon=0.1Fx0=D_x*sin(C_x*atan(B_x*k_x-E_x.*(B_x*k_x-atan(B_x*k_x))))+S_vx;figure(1);plot(k,Fx0,'linewidth',1.5);gridxlabel('纵向
4、滑移率');ylabel('纵向力/（N）');title('纵向力--滑移率(纯纵滑)');汽车系统动力学作业12.纯侧滑情况：%%LateralForceandTorque(puresideslip)%%aisinputvalue,meansslipangle%%Fy0isoutputvalue,meanssideforce%%Mz0isoutputvalue,meansaligningtorquea=linspace(-15,15,300);%侧偏角取-15-15度ta=tan(a*pi/180);sr=sin
5、(r);%r是外倾角，sr表示r*C_y=P_Cy1;%取lam(Cy)=1u_y=(P_Dy1+P_Dy2*dfz)/(1+P_Dy3*sr^2);%取lam(uy*)=1D_y=u_y*Fz;%zeta(1-8)都取1K_ya=P_Ky1*Fz0*sin(P_Ky4*atan(Fz/((P_Ky2+P_Ky5*sr^2)*Fz0)))/(1+P_Ky3*sr^2);%取lam（Kya）=1zeta(3)=1B_y=K_ya/(C_y*D_y);%取epsilony=0；K_yr0=Fz*(P_Ky6+P_Ky7*
6、dfz);%取lam（Kyr)=1S_Vyr=Fz*(P_Vy3+P_Vy4*dfz)*sr;%lam(kyr)=1,camberforcestiffness后面的修正系数全部取1S_Hy=(P_Hy1+P_Hy2*dfz)+(K_yr0*sr-S_Vyr)/(K_ya+0.1);%zeta(0,4)都取1epsilon(k)=0.1S_Vy=Fz*(P_Vy1+P_Vy2*dfz)+S_Vyr;%lam(Vy,uy)=1,zeta(2)=1ay=ta+S_Hy;%tan(sa)表示横向侧偏角a*E_y=(P_Ey1
7、+P_Ey2*dfz)*(1+P_Ey5*sr^2-(P_Ey3+P_Ey4*sr)*sign(ay));%lam(ey)=1Fy0=D_y*sin(C_y*atan(B_y*ay-E_y.*(B_y*ay-atan(B_y*ay))))+S_Vy;figure(2);plot(a,Fy0,'linewidth',1.5);gridxlabel('侧偏角');ylabel('侧向力/（N）');title('侧向力-侧偏角(纯侧滑)');C_r=1;C_t=q_Cz1;B_t=(q_Bz1+q_Bz2*dfz+q_B
8、z3*dfz^2)*(1+q_Bz5*abs(sr)+q_Bz6*sr^2);S_Ht=q_Hz1+q_Hz2*dfz+(q_Hz3+q_Hz4*dfz)*sr;a_t=ta+S_Ht;B_r=(q_Bz9+q_Bz10*B_y*C_y);E_t=(q_Ez1+q_Ez2*dfz+q_Ez3*dfz^2)*(1+(q_Ez4+q_Ez5*sr)
