轮胎魔术公式.zip

里面为魔术轮胎公式的仿真程序，为m文件，助于对轮胎模型的理解

%-------------轮胎魔术公式--------------
%% 参数设置
a0 = 1.65; a1 = -34; a2 = 1250; a3 = 3036; a4 = 12.8; a5 = 0.00501; a6 = -0.02103; 
a7 = 0.77394; a8 = 0.0022890; a10 = 0.00370; a11 = 19.1656;a12 =1.21356; a13 = 6.26206;
b0 = 2.37272; b1 =-9.46; b2 = 1490; b3 = 130; b4 = 276; b5 = 0.0886; b6 = 0.00402;
b7 = -0.0615; b8 = 1.2; b9 = 0.0299; b10 = -0.176;
c0 = 2.340; c1 = 1.495; c2 = 6.416654; c3 = -3.57403; c4 = -0.087737;c5 = 0.098410;
c6 = 0.0027699; c7 = -0.0001151; c8 = 0.1; c9 = -1.33329; c10 = 0.025501; c11 = -0.02357; 
c12 = 0.03027; c13 = -0.0647; c14 = 0.0211329; c15 = 0.89469; c16 = -0.099443; c17 = -3.336941;
%% 纵向力计算
Fz = 5;                       % 垂直载荷：kN
s = linspace（-202040）;
Sh = b9*Fz+b10;               % 水平方向漂移
Sv = 0;                       % 垂直方向漂移
C = b0;                       % 曲线形状因子
D = b1*Fz^2+b2*Fz;            % 曲线巅因子
E = b6*Fz^2+b7*Fz+b8;         % 曲线曲率因子
B = （b3*Fz^2+b4*Fz）*exp（-b5*Fz）/（C*D）;
x = s + Sh;
F1 = D*sin（C*atan（B*x-E*（B*x - atan（B*x））））+Sv;
figure 
plot（sF1‘r‘‘LineWidth‘2）
xlabel（‘滑移率s‘）;
ylabel（‘纵向力F1‘）;
title（‘滑移率与纵向力的关系‘）
grid minor

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2019-01-12 20:17  轮胎魔术公式\
     文件        1215  2019-01-12 09:38  轮胎魔术公式\tyreMagic.m