轮胎特性图绘制

最新推荐文章于 2024-04-18 19:34:50 发布
最新推荐文章于 2024-04-18 19:34:50 发布
阅读量982 收藏 18
点赞数 3
分类专栏： matlab 文章标签： matlab
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/m0_62299908/article/details/129869101
版权

参考文献《基于“魔术公式”的轮胎动力学仿真分析》。
1、 不同垂直载荷Fz的纵向力Fx图
（1）搭建simulink模型。
在这里插入图片描述
（2）在matlab function 模块中创建魔术公式轮胎。

function Fx0 = pac_tire(lambda,Fz0)
%魔术轮胎公式——期刊《基于魔术公式的轮胎动力学仿真分析》2012
% gamma = 0;%外倾角为0。
%不考虑车轮外倾角和漂移的影响，拟合出的以下参数（即gamma=0;Sh=Sv=0。）
Fz=Fz0/1000;%单位：kN。
%% 纵向力Fx
a0=1.65;
a1=-21.3;
a2=1144;
a3=49.6;
a4=226;
a5=0.069;
a6=-0.006;
a7=0.056;
a8=0.486;
%% 参数Fx
Cx=a0;%曲线形状因子
Dx=(a1*Fz^2+a2*Fz);%峰值因子（曲线的最大值）
BCDx=(a3*Fz^2+a4*Fz)*exp(-a5*Fz);
Bx=BCDx/(Cx*Dx);%刚度因子
Ex=a6*Fz^2+a7*Fz+a8;%曲线曲率因子（曲线最大值附近的形状）
% Shx=a9*Fz+a10;%曲线的水平方向漂移
Shx=0;
Svx=0;%曲线的垂直方向漂移
X=(lambda*100+Shx);%此处lambda*100;lambda∈【0，1】。
% Fx0=(a1*Fz^2+a2*Fz)*sin(a0*atan(((a3*Fz^2+a4*Fz)*exp(-a5*Fz)/(a0*(a1*Fz^2+a2*Fz)))*(x+(a9*Fz+a10))-(a6*Fz^2+a7*Fz+a8)*(((a3*Fz^2+a4*Fz)*exp(-a5*Fz)/(a0*(a1*Fz^2+a2*Fz)))*(x+(a9*Fz+a10))-atan(((a3*Fz^2+a4*Fz)*exp(-a5*Fz)/(a0*(a1*Fz^2+a2*Fz)))*(x+(a9*Fz+a10))))));
Fx0=Dx*sin(Cx*atan(Bx*X-Ex*(Bx*X-atan(Bx*X)))) + Svx;

（3）仿真工况一：Ramp模块的初始值设置为0，仿真时间设置为1s。
在这里插入图片描述
仿真工况二：Ramp模块的初始值设置为-1，仿真时间设置为2s。
在这里插入图片描述
2、 不同垂直载荷Fz的侧向力Fy图。
（1）搭建simulink模型。
在这里插入图片描述
（2）在matlab function 模块中创建魔术公式轮胎。

function Fy0 = pac_tire(alpha,Fz0)
%魔术轮胎公式——期刊《基于魔术公式的轮胎动力学仿真分析》2012
gamma = 0;%外倾角为0。
%不考虑车轮外倾角和漂移的影响，拟合出的以下参数（即gamma=0;Sh=Sv=0。）
Fz=Fz0/1000;%单位：kN。
%% 稳态纯转向工况下侧向力Fy0
%拟合参数Fy0
b0=1.3;
b1=-22.1;
b2=1011;
b3=1078;
b4=1.82;
b5=0.208;
b6=0;
b7=-0.354;
b8=0.707;
b9=0.028;
b10=0;
b11=14.8;
b12=0.022;
%% 参数Fy0
Cy=b0;
Dy=(b1*Fz^2+b2*Fz);
%有点不一样
BCDy=b3*sin(b4*atan(b5*Fz))*(1-b12*abs(gamma));%gamma为外倾角。
By=BCDy/(Cy*Dy);
Ey=b6*Fz^2+b7*Fz+b8;
Shy=b9*gamma;
Svy=(b10*Fz^2+b11*Fz)*gamma;
% 
Y=alpha+Shy;%侧偏角alpha_y
%Fy0——正的测偏角，负的侧偏力。
Fy0=-Dy*sin(Cy*atan(By*Y-Ey*(By*Y-atan(By*Y)))) + Svy;

（3）仿真工况一：Ramp模块的初始值设置为0deg，仿真时间设置为20s。
在这里插入图片描述
仿真工况二：Ramp模块的初始值设置为20deg，仿真时间设置为40s。
在这里插入图片描述
3、侧偏角为5°时，纵向力Fx和侧向力Fy，随滑移率的变化曲线。
（1）搭建simulink模型。
在这里插入图片描述
（2）在matlab function 模块中创建魔术公式轮胎。

function [Fx,Fy] = pac_tire(lambda,alpha,Fz0)
%魔术轮胎公式——期刊《基于魔术公式的轮胎动力学仿真分析》2012
gamma = 0;%外倾角为0。
%不考虑车轮外倾角和漂移的影响，拟合出的以下参数（即gamma=0;Sh=Sv=0。）
% 纵向力中的lambda=lambda*100;%滑移率【0%，100%】。
Fz=Fz0/1000;%单位：kN。
%% 纵向力Fx
a0=1.65;
a1=-21.3;
a2=1144;
a3=49.6;
a4=226;
a5=0.069;
a6=-0.006;
a7=0.056;
a8=0.486;
%% 参数Fx
Cx=a0;%曲线形状因子
Dx=(a1*Fz^2+a2*Fz);%峰值因子（曲线的最大值）
BCDx=(a3*Fz^2+a4*Fz)*exp(-a5*Fz);
Bx=BCDx/(Cx*Dx);%刚度因子
Ex=a6*Fz^2+a7*Fz+a8;%曲线曲率因子（曲线最大值附近的形状）
Shx=0;%曲线的水平方向漂移
Svx=0;%曲线的垂直方向漂移
X=(lambda*100+Shx);%此处lambda*100;lambda∈【0，1】。
% Fx0=(a1*Fz^2+a2*Fz)*sin(a0*atan(((a3*Fz^2+a4*Fz)*exp(-a5*Fz)/(a0*(a1*Fz^2+a2*Fz)))*(x+(a9*Fz+a10))-(a6*Fz^2+a7*Fz+a8)*(((a3*Fz^2+a4*Fz)*exp(-a5*Fz)/(a0*(a1*Fz^2+a2*Fz)))*(x+(a9*Fz+a10))-atan(((a3*Fz^2+a4*Fz)*exp(-a5*Fz)/(a0*(a1*Fz^2+a2*Fz)))*(x+(a9*Fz+a10))))));
Fx0=Dx*sin(Cx*atan(Bx*X-Ex*(Bx*X-atan(Bx*X)))) + Svx;
%% %魔术轮胎公式——期刊《基于魔术公式的轮胎动力学仿真分析》2012
%不考虑车轮外倾角和漂移的影响，拟合出的以下参数（即gamma=0;Sh=Sv=0。）
% 侧偏力中的alpha=alpha*180/pi;%单位换算rad_to_deg。
%% 稳态纯转向工况下侧向力Fy0
%拟合参数Fy0
b0=1.3;
b1=-22.1;
b2=1011;
b3=1078;
b4=1.82;
b5=0.208;
b6=0;
b7=-0.354;
b8=0.707;
b9=0.028;
b10=0;
b11=14.8;
b12=0.022;
%% 参数Fy0
Cy=b0;
Dy=(b1*Fz^2+b2*Fz);
BCDy=b3*sin(b4*atan(b5*Fz))*(1-b12*abs(gamma));%gamma为外倾角。
By=BCDy/(Cy*Dy);
Ey=b6*Fz^2+b7*Fz+b8;
Shy=b9*gamma;
Svy=(b10*Fz^2+b11*Fz)*gamma;
Y=alpha+Shy;%侧偏角alpha_y
%% %Fy0——正的测偏角（deg），负的侧偏力。
Fy0=-Dy*sin(Cy*atan(By*Y-Ey*(By*Y-atan(By*Y)))) + Svy;
%% %联合仿真工况
sigma_x=lambda/(1+abs(lambda));%此处lambda输入【0，1】。
sigma_y=tan(alpha*pi/180)/(1+abs(lambda));%此处alpha的单位为rad。
sigma=sqrt(sigma_x^2+sigma_y^2);
Fx=sign(lambda)*Fx0*sigma_x/sigma;
Fy=sign(alpha)*Fy0*sigma_y/sigma;
% %联合仿真工况
% sigma_x=-lambda/(1+abs(lambda));%此处lambda输入【0，1】。
% sigma_y=-tan(alpha*pi/180)/(1+abs(lambda));%此处alpha的单位为rad。
% sigma=sqrt(sigma_x^2+sigma_y^2);
% Fx=Fx0*sigma_x/sigma;
% Fy=Fy0*sigma_y/sigma;

（3）仿真工况：Ramp模块的初始值设置为0，仿真时间设置为1s。
在这里插入图片描述
4、载荷4kN，外倾角为0，不同侧偏角。
（1）搭建simulink模型。
在这里插入图片描述
（2）在matlab function 模块中创建魔术公式轮胎。

function [Fx,Fy] = pac_tire(lambda,alpha,Fz0)
%魔术轮胎公式——期刊《基于魔术公式的轮胎动力学仿真分析》2012
gamma = 0;%外倾角为0。
%不考虑车轮外倾角和漂移的影响，拟合出的以下参数（即gamma=0;Sh=Sv=0。）
% 纵向力中的lambda=lambda*100;%滑移率【0%，100%】。
Fz=Fz0/1000;%单位：kN。
%% 纵向力Fx
a0=1.65;
a1=-21.3;
a2=1144;
a3=49.6;
a4=226;
a5=0.069;
a6=-0.006;
a7=0.056;
a8=0.486;
%% 参数Fx
Cx=a0;%曲线形状因子
Dx=(a1*Fz^2+a2*Fz);%峰值因子（曲线的最大值）
BCDx=(a3*Fz^2+a4*Fz)*exp(-a5*Fz);
Bx=BCDx/(Cx*Dx);%刚度因子
Ex=a6*Fz^2+a7*Fz+a8;%曲线曲率因子（曲线最大值附近的形状）
Shx=0;%曲线的水平方向漂移
Svx=0;%曲线的垂直方向漂移
X=(lambda*100+Shx);%此处lambda*100;lambda∈【0，1】。
% Fx0=(a1*Fz^2+a2*Fz)*sin(a0*atan(((a3*Fz^2+a4*Fz)*exp(-a5*Fz)/(a0*(a1*Fz^2+a2*Fz)))*(x+(a9*Fz+a10))-(a6*Fz^2+a7*Fz+a8)*(((a3*Fz^2+a4*Fz)*exp(-a5*Fz)/(a0*(a1*Fz^2+a2*Fz)))*(x+(a9*Fz+a10))-atan(((a3*Fz^2+a4*Fz)*exp(-a5*Fz)/(a0*(a1*Fz^2+a2*Fz)))*(x+(a9*Fz+a10))))));
Fx0=Dx*sin(Cx*atan(Bx*X-Ex*(Bx*X-atan(Bx*X)))) + Svx;
%% %魔术轮胎公式——期刊《基于魔术公式的轮胎动力学仿真分析》2012
%不考虑车轮外倾角和漂移的影响，拟合出的以下参数（即gamma=0;Sh=Sv=0。）
% 侧偏力中的alpha=alpha*180/pi;%单位换算rad_to_deg。
%% 稳态纯转向工况下侧向力Fy0
%拟合参数Fy0
b0=1.3;
b1=-22.1;
b2=1011;
b3=1078;
b4=1.82;
b5=0.208;
b6=0;
b7=-0.354;
b8=0.707;
b9=0.028;
b10=0;
b11=14.8;
b12=0.022;
%% 参数Fy0
Cy=b0;
Dy=(b1*Fz^2+b2*Fz);
BCDy=b3*sin(b4*atan(b5*Fz))*(1-b12*abs(gamma));%gamma为外倾角。
By=BCDy/(Cy*Dy);
Ey=b6*Fz^2+b7*Fz+b8;
Shy=b9*gamma;
Svy=(b10*Fz^2+b11*Fz)*gamma;
Y=alpha+Shy;%侧偏角alpha_y
%% %Fy0——正的测偏角（deg），负的侧偏力。
Fy0=-Dy*sin(Cy*atan(By*Y-Ey*(By*Y-atan(By*Y)))) + Svy;
%% %联合仿真工况
sigma_x=lambda/(1+abs(lambda));%此处lambda输入【0，1】。
sigma_y=tan(alpha*pi/180)/(1+abs(lambda));%此处alpha的单位为rad。
sigma=sqrt(sigma_x^2+sigma_y^2);
Fx=sign(lambda)*Fx0*sigma_x/sigma;
Fy=sign(alpha)*Fy0*sigma_y/sigma;
% %联合仿真工况
% sigma_x=-lambda/(1+abs(lambda));%此处lambda输入【0，1】。
% sigma_y=-tan(alpha*pi/180)/(1+abs(lambda));%此处alpha的单位为rad。
% sigma=sqrt(sigma_x^2+sigma_y^2);
% Fx=Fx0*sigma_x/sigma;
% Fy=Fy0*sigma_y/sigma;

（3）仿真工况一：Ramp模块的初始值设置为0deg，仿真时间设置为1s。
在这里插入图片描述
在这里插入图片描述
（4）在matlab中绘制附着椭圆图。
1）代码

%matlab绘图命令
clc
x1=Fx(:,1);
x2=Fx(:,2);
x3=Fx(:,3);
x4=Fx(:,4);
y1=Fy(:,1);
y2=Fy(:,2);
y3=Fy(:,3);
y4=Fy(:,4);
hold on
plot(x1,y1);
plot(x2,y2);
plot(x3,y3);
plot(x4,y4);
legend("1°","3°","5°","7°");
xlabel("纵向力Fx");
ylabel("侧向力Fy" );

2）附着椭圆
在这里插入图片描述
5、不同路面附着系数mu下的侧偏力曲线。
（1）搭建simulink模型。
在这里插入图片描述
（2）在matlab function 模块中创建魔术公式轮胎。

function [Fx,Fy] = pac_tire(lambda,alpha,Fz0,mu)
%魔术轮胎公式——期刊《基于魔术公式的轮胎动力学仿真分析》2012
gamma = 0;%外倾角为0。
%不考虑车轮外倾角和漂移的影响，拟合出的以下参数（即gamma=0;Sh=Sv=0。）
% 纵向力中的lambda=lambda*100;%滑移率【0%，100%】。
Fz=Fz0/1000;%单位：kN。
%% 纵向力Fx
a0=1.65;
a1=-21.3;
a2=1144;
a3=49.6;
a4=226;
a5=0.069;
a6=-0.006;
a7=0.056;
a8=0.486;
%% 参数Fx
Cx=a0;%曲线形状因子
Dx=(a1*Fz^2+a2*Fz) * mu;%峰值因子（曲线的最大值）
BCDx=(a3*Fz^2+a4*Fz)*exp(-a5*Fz);
Bx=BCDx/(Cx*Dx);%刚度因子
Ex=a6*Fz^2+a7*Fz+a8;%曲线曲率因子（曲线最大值附近的形状）
Shx=0;%曲线的水平方向漂移
Svx=0;%曲线的垂直方向漂移
X=(lambda*100+Shx);%此处lambda*100;lambda∈【0，1】。
% Fx0=(a1*Fz^2+a2*Fz)*sin(a0*atan(((a3*Fz^2+a4*Fz)*exp(-a5*Fz)/(a0*(a1*Fz^2+a2*Fz)))*(x+(a9*Fz+a10))-(a6*Fz^2+a7*Fz+a8)*(((a3*Fz^2+a4*Fz)*exp(-a5*Fz)/(a0*(a1*Fz^2+a2*Fz)))*(x+(a9*Fz+a10))-atan(((a3*Fz^2+a4*Fz)*exp(-a5*Fz)/(a0*(a1*Fz^2+a2*Fz)))*(x+(a9*Fz+a10))))));
Fx0=Dx*sin(Cx*atan(Bx*X-Ex*(Bx*X-atan(Bx*X)))) + Svx;
%% %魔术轮胎公式——期刊《基于魔术公式的轮胎动力学仿真分析》2012
%不考虑车轮外倾角和漂移的影响，拟合出的以下参数（即gamma=0;Sh=Sv=0。）
% 侧偏力中的alpha=alpha*180/pi;%单位换算rad_to_deg。
%% 稳态纯转向工况下侧向力Fy0
%拟合参数Fy0
b0=1.3;
b1=-22.1;
b2=1011;
b3=1078;
b4=1.82;
b5=0.208;
b6=0;
b7=-0.354;
b8=0.707;
b9=0.028;
b10=0;
b11=14.8;
b12=0.022;
%% 参数Fy0
Cy=b0;
Dy=(b1*Fz^2+b2*Fz) * mu;
BCDy=b3*sin(b4*atan(b5*Fz))*(1-b12*abs(gamma));%gamma为外倾角。
By=BCDy/(Cy*Dy);
Ey=b6*Fz^2+b7*Fz+b8;
Shy=b9*gamma;
Svy=(b10*Fz^2+b11*Fz)*gamma;
Y=alpha+Shy;%侧偏角alpha_y
%% %Fy0——正的测偏角（deg），负的侧偏力。
Fy0=-Dy*sin(Cy*atan(By*Y-Ey*(By*Y-atan(By*Y)))) + Svy;
%% %联合仿真工况
sigma_x=lambda/(1+abs(lambda));%此处lambda输入【0，1】。
sigma_y=tan(alpha*pi/180)/(1+abs(lambda));%此处alpha的单位为rad。
sigma=sqrt(sigma_x^2+sigma_y^2);
Fx=sign(lambda)*Fx0*sigma_x/sigma;
Fy=sign(alpha)*Fy0*sigma_y/sigma;

（3）仿真工况一：Ramp模块的初始值设置为-20deg，仿真时间设置为40s。
在这里插入图片描述
5、使用carsim软件中的轮胎参数表格绘制三维图。
将carsim中轮胎参数excel表格另存，表格重命名为“轮胎参数”。
在这里插入图片描述
轮胎参数excel表格中的数据命名如下：
在这里插入图片描述

在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
程序：

%% 三维魔术轮胎曲线拟合
% y=d*sin(c*atan(b*x-e*(b*x-atan(b*x))));%魔术公式轮胎
%% 使用工具箱z=f(x,y)——Fx/Fy/Mz=(lambda/alpha,Fz);对应X——；Y——Fz;z——ZZ;
% y=(a1*y^2+a2*y)*y*sin(c*atan(b*x-e*(b*x-atan(b*x))));%曲线拟合公式
% % Fx0=(a1*Fz^2+a2*Fz)*sin(a0*atan(((a3*Fz^2+a4*Fz)*exp(-a5*Fz)/(a0*(a1*Fz^2+a2*Fz)))*lambda*100-(a6*Fz^2+a7*Fz+a8)*(((a3*Fz^2+a4*Fz)*exp(-a5*Fz)/(a0*(a1*Fz^2+a2*Fz)))*lambda*100-atan(((a3*Fz^2+a4*Fz)*exp(-a5*Fz)/(a0*(a1*Fz^2+a2*Fz)))*lambda*100))));
% % Fy0=-(b1*Fz^2+b2*Fz)*sin(b0*atan((b3*sin(b4*atan(b5*Fz))/(b0*(b1*Fz^2+b2*Fz)))*alpha-(b6*Fz^2+b7*Fz+b8)*((b3*sin(b4*atan(b5*Fz))/(b0*(b1*Fz^2+b2*Fz)))*alpha-atan((b3*sin(b4*atan(b5*Fz))/(b0*(b1*Fz^2+b2*Fz)))*alpha))));
% % Mz=(c1*Fz^2+c2*Fz)*sin(c0*atan(((c3*Fz^2+c4*Fz)*exp(-c5*Fz)/(c0*(c1*Fz^2+c2*Fz)))*alpha-(c6*Fz^2+c7*Fz+c8)*(((c3*Fz^2+c4*Fz)*exp(-c5*Fz)/(c0*(c1*Fz^2+c2*Fz)))*alpha-atan(((c3*Fz^2+c4*Fz)*exp(-c5*Fz)/(c0*(c1*Fz^2+c2*Fz)))*alpha))));

%% x——使用工具箱拟合参数
% % Fx0=(a1*Fz^2+a2*Fz)*sin(a0*atan(((a3*Fz^2+a4*Fz)*exp(-a5*Fz)/(a0*(a1*Fz^2+a2*Fz)))*x*100-(a6*Fz^2+a7*Fz+a8)*(((a3*Fz^2+a4*Fz)*exp(-a5*Fz)/(a0*(a1*Fz^2+a2*Fz)))*x*100-atan(((a3*Fz^2+a4*Fz)*exp(-a5*Fz)/(a0*(a1*Fz^2+a2*Fz)))*x*100))));
% % Fy0=-(b1*Fz^2+b2*Fz)*sin(b0*atan((b3*sin(b4*atan(b5*Fz))/(b0*(b1*Fz^2+b2*Fz)))*x-(b6*Fz^2+b7*Fz+b8)*((b3*sin(b4*atan(b5*Fz))/(b0*(b1*Fz^2+b2*Fz)))*x-atan((b3*sin(b4*atan(b5*Fz))/(b0*(b1*Fz^2+b2*Fz)))*x))));
% % Mz=(c1*Fz^2+c2*Fz)*sin(c0*atan(((c3*Fz^2+c4*Fz)*exp(-c5*Fz)/(c0*(c1*Fz^2+c2*Fz)))*x-(c6*Fz^2+c7*Fz+c8)*(((c3*Fz^2+c4*Fz)*exp(-c5*Fz)/(c0*(c1*Fz^2+c2*Fz)))*x-atan(((c3*Fz^2+c4*Fz)*exp(-c5*Fz)/(c0*(c1*Fz^2+c2*Fz)))*x))));
%% Excel表格数据由Carsim得到
clc;close all;clear
%% 轮胎纵向力Fx
lambda=xlsread("轮胎数据","Fx","A2:A52");%第一列
Fz1=xlsread("轮胎数据","Fx","B1:I1")/1000;%单位kN。%第一行
y11=xlsread("轮胎数据","Fx","B2:B52");%列%%y11-y81数据主要用于使用工具箱进行曲线拟合效果观察。
y21=xlsread("轮胎数据","Fx","C2:C52");
y31=xlsread("轮胎数据","Fx","D2:D52");
y41=xlsread("轮胎数据","Fx","E2:E52");
y51=xlsread("轮胎数据","Fx","F2:F52");%列
y61=xlsread("轮胎数据","Fx","G2:G52");
y71=xlsread("轮胎数据","Fx","H2:H52");
y81=xlsread("轮胎数据","Fx","I2:I52");
zz=xlsread("轮胎数据","Fx","B2:I52");%框选Fx\Fy\Mz的取值范围【左上角，右下角】
figure(1)
[x,y]=meshgrid(lambda,Fz1);%二维图%%%[x, y] = meshgrid(x,y);% 快速生成网格所需的数据
z=zz';%转置。
% subplot(1,2,1)
Fx=mesh(x,y,z);%%三维网格图mesh(x,y,z);%%x是n维向量，y是m维向量，z是m*n维向量
% subplot(1,2,2)
% surf(x,y,z) % (X(j), Y(i), Z(i,j))是线框网格线的交点
axis vis3d % 冻结屏幕高宽比，使得一个三维对象的旋转不会改变坐标轴的刻度显示
%% 标注x
xlabel("滑移率(S)");
ylabel("垂直载荷(Fz)");
zlabel("纵向力Fx"); 
%% 轮胎侧向力Fy
alpha=xlsread("轮胎数据","Fy","A2:A52");%第一列
Fz2=xlsread("轮胎数据","Fy","B1:I1")/1000;%单位kN。%第一行
y12=xlsread("轮胎数据","Fy","B2:B52");%列
y22=xlsread("轮胎数据","Fy","C2:C52");
y32=xlsread("轮胎数据","Fy","D2:D52");
y42=xlsread("轮胎数据","Fy","E2:E52");
y52=xlsread("轮胎数据","Fy","F2:F52");%列
y62=xlsread("轮胎数据","Fy","G2:G52");
y72=xlsread("轮胎数据","Fy","H2:H52");
y82=xlsread("轮胎数据","Fy","I2:I52");
zz2=xlsread("轮胎数据","Fy","B2:I52");%框选Fx\Fy\Mz的取值范围【左上角，右下角】
figure(2)
[x2,y2]=meshgrid(alpha,Fz2);%二维图
z2=zz2';%转置。
Fy=mesh(x2,y2,z2);%三维图
surf(x2,y2,z2) % (X(j), Y(i), Z(i,j))是线框网格线的交点
%% 标注y
xlabel("侧偏角alpha/（°）");
ylabel("垂直载荷Fz/（N）");
zlabel("侧偏力Fy/（N）");
%% 轮胎回正力矩
alpha=xlsread("轮胎数据","Mz","A2:A52");%第一列
Fz3=xlsread("轮胎数据","Mz","B1:I1")/1000;%单位kN。%第一行
y13=xlsread("轮胎数据","Mz","B2:B52");%列
y23=xlsread("轮胎数据","Mz","C2:C52");
y33=xlsread("轮胎数据","Mz","D2:D52");
y43=xlsread("轮胎数据","Mz","E2:E52");
y53=xlsread("轮胎数据","Mz","F2:F52");%列
y63=xlsread("轮胎数据","Mz","G2:G52");
y73=xlsread("轮胎数据","Mz","H2:H52");
y83=xlsread("轮胎数据","Mz","I2:I52");
zz3=xlsread("轮胎数据","Mz","B2:I52");%框选Fx\Fy\Mz的取值范围【左上角，右下角】
figure(3)
[x3,y3]=meshgrid(alpha,Fz3);%二维图
z3=zz3';%转置。
subplot(1,2,1)
Mz=mesh(x3,y3,z3);%三维网格图mesh(x,y,z)
subplot(1,2,2)
surf(x3,y3,z3) % (X(j), Y(i), Z(i,j))是线框网格线的交点
%% 标注z
xlabel("侧偏角alpha/（°）");
ylabel("垂直载荷Fz/（kN）");
zlabel("回正力矩Mz（N）");
%% 读取 Microsoft Excel 电子表格文件
%% num = xlsread(filename,sheet,xlRange)
%SSE为方差
%均方根植接近于1是最好的。
%% 标注
% xlabel("滑移率/侧偏角");
% ylabel("垂直载荷Fz");
% zlabel("纵向力Fx/侧偏力Fy/回正力矩Mz");

%% 魔术公式轮胎
%% Fx
% Cx=a0;%曲线形状因子
% Dx=(a1*Fz^2+a2*Fz);%峰值因子（曲线的最大值）
% BCDx=(a3*Fz^2+a4*Fz)*exp(-a5*Fz);
% Bx=BCDx/(Cx*Dx);%刚度因子
% Ex=(a6*Fz^2+a7*Fz+a8);%曲线曲率因子（曲线最大值附近的形状）
% X=lambda*100;%滑移率使用百分比。
% Fx0=Dx*sin(Cx*atan(Bx*X-Ex*(Bx*X-atan(Bx*X))));
% % Bx=((a3*Fz^2+a4*Fz)*exp(-a5*Fz)/(a0*(a1*Fz^2+a2*Fz)));
% % Fx0=(a1*Fz^2+a2*Fz)*sin(a0*atan(((a3*Fz^2+a4*Fz)*exp(-a5*Fz)/(a0*(a1*Fz^2+a2*Fz)))*lambda*100-(a6*Fz^2+a7*Fz+a8)*(((a3*Fz^2+a4*Fz)*exp(-a5*Fz)/(a0*(a1*Fz^2+a2*Fz)))*lambda*100-atan(((a3*Fz^2+a4*Fz)*exp(-a5*Fz)/(a0*(a1*Fz^2+a2*Fz)))*lambda*100))));
%% Fy——忽略外倾角。
% Cy=b0;
% Dy=(b1*Fz^2+b2*Fz);
% % BCDy=b3*sin(b4*atan(b5*Fz))*(1-b12*abs(gamma));%gamma为外倾角。
% BCDy=b3*sin(b4*atan(b5*Fz));
% By=BCDy/(Cy*Dy);
% Ey=(b6*Fz^2+b7*Fz+b8);
% Y=alpha;
% Fy0=-Dy*sin(Cy*atan(By*Y-Ey*(By*Y-atan(By*Y))));%Fy0——正的测偏角（deg），负的侧偏力。
% % By=(b3*sin(b4*atan(b5*Fz))/(b0*(b1*Fz^2+b2*Fz)));
% % Fy0=-(b1*Fz^2+b2*Fz)*sin(b0*atan((b3*sin(b4*atan(b5*Fz))/(b0*(b1*Fz^2+b2*Fz)))*alpha-(b6*Fz^2+b7*Fz+b8)*((b3*sin(b4*atan(b5*Fz))/(b0*(b1*Fz^2+b2*Fz)))*alpha-atan((b3*sin(b4*atan(b5*Fz))/(b0*(b1*Fz^2+b2*Fz)))*alpha))));

结束。

、南城不南 ╯
关注
  • 3
    点赞
  • 踩
  • 18
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
轮胎模型matlab代码-Matlab_vehicle_modelling:MatlabEV仿真模型，用于学士学位论文
05-27
轮胎模型matlab代码Matlab_vehicle_modelling Matlab EV仿真模型，用于学士学位论文 %%读我%% 要运行模拟，必须遵循以下步骤。 确保将Microsoft Office Excel设置为将点用作小数点分隔符，将逗号用作千位分隔符（加载自定义GPC驱动器循环功能需要此功能） 确保已下载并安装了“ MathWorks Automotive的动力总成模块组行驶周期数据”附加组件。 这可以从matlab本身完成，方法是打开“加载项”页面并搜索“动力总成模块组行驶周期”（启用包含正式行驶周期（例如NEDC和EPA行驶周期）的行驶周期块是必需的） 一旦完成了前面的步骤，就可以启动模型。 这将为模型提供所需的车辆布局参数，并指定牵引条件，无论是干停机坪，湿塔玛克面包，雪还是冰。 要初始化模型，请打开initializeModel.m脚本并运行它。 在命令提示符下，指定车辆的参数 一旦运行了initializeModel.m脚本，就可以选择并运行EV_basic_8_1.slx Simulink模型文件。 如果遵循了前面的步骤，则不会出现任何错误代码。 （必须注意，
魔术轮胎matlab建模
wqj5211314的博客
10-13 3553
利用课题组给出的数据搭建魔术轮胎模型，也为整车模型构建做准备
基于Simulink的魔术公式：轮胎动力学模型实现路面附着系数、垂直载荷和车轮侧偏角的输入输出关系
最新发布
WVCaBKKF的博客
04-18 614
总结而言，基于Simulink的轮胎动力学模型是一种有效的工具，用于对车辆在不同路面条件下的行驶性能进行仿真和分析。基于Simulink的轮胎动力学模型，是一种基于魔术公式的模拟工具，用于对车辆在不同路面条件下的行驶性能进行分析和优化。附着椭圆是轮胎与路面之间接触力的分布情况的一种表示方法，它能够直观地反映轮胎的摩擦力分布情况以及车辆的稳定性。因此，在模型中需要考虑如何准确地估计车辆的垂直载荷。通过调整输入参数，我们可以得到不同工况下的纵向力、侧向力和附着椭圆的变化情况，进而评估车辆的操控性能和安全性。
魔术轮胎轮胎侧向力matlab建模
CScodemini的博客
07-13 1438
魔术公式”采用三角函数的组合公式拟合实验轮胎数据，用一套形式相同的公式就可以完整的表达纵向力、侧向力、回正力矩、反转力矩、滚动阻力以及纵向力、侧向力联合作用的工况魔术 公 式 统一性强，需要拟 合 的 参 数 都 有 明 确 的 物 理 意 义，容易确定初值。并且由于魔术公式是基于试验数据，处在试验的范围的高精度外，甚至在极限值以外的一定程度也可使用，它拟合精度高，但计算量大，因此更适用于产品设计、汽车动态模拟及试验对比等要求精确描述轮胎力学特性的领域。
Matlab拟合工具箱拟合轮胎侧偏刚度
YOURUOLI的博客
03-09 859
使用matlab曲线拟合工具对轮胎侧偏刚度进行拟合。
从0开始建立车辆仿真模型 – 车辆轮胎模型的建立
windsnow1991的博客
10-19 2995
从0开始建立车辆仿真模型 – 车辆轮胎模型的建立 本来还想继续把MBD有关代码生成的开发的细节再讲讲，后来我的“产品经理”炸了，开始催更物理模型的文章了，那就先来说说物理模型的话题，后面再去补全代码生成的部分。 那首先来说说为什么要建立物理模型，物理模型一方面可以快速建立被控对象的模型，去支持算法和功能的仿真。另一方面物理模型也可以用来支持Mil / Hil测试，对于真实刷写到XCU中的算法代码进行验证。 本文先从0D/1D的系统仿真的层面介绍物理模型的搭建。 首先从系统层面来看，车辆系统级别的需求是整车
MATLAB绘制万有特性图，完美绘制万有特性图
12-29
MATLAB绘制万有特性图，非最小二乘法（测试网上流传的程序后发现绘制的万有特性完全失真），自己摸索的新程序可完美绘制万有特性
气象领域python色斑图绘制程序
06-09
可绘制降水分布色斑图，温度分布色斑图，对应的图例数据，标题数据，可绘制中心城市名称，添加标注等
Python折线图绘制
04-20
python绘制折线图
Scratch绘制热点图
02-20
Scratch绘制热点图
PHP绘制热点图123
02-20
PHP绘制热点图123
魔术公式.matlab实现代码 纵向力与侧向力实现
06-03
魔术公式
轮胎模型matlab代码-TTC_FSAE:FormulaSAE轮胎测试联盟
05-27
轮胎模型matlab代码Formula SAE轮胎测试联盟- 该项目的目标是利用的数据来协助UFSC方程式开发的原型悬架和转向系统的开发。 关于该项目 必须开发设计良好的汽车悬架和转向系统，并考虑到轮胎在不同情况下的性能。 考虑到这一点，FSAE TTC进行了几次测试，其中轮胎在以下类似机器中进行了测试，可以测量每个轮胎对不同情况的React。 该数据可用于加入该联盟的Formula SAE团队。 在RunDataCorneringTest.m文件中找到的Matlab代码根据exampleData.DAT数据生成5个数字。 图1和2提供了一般的可视化效果，因此用户可以大致了解从数据中得到的期望，而图3、4和5显示了3d图形，其中针对三个不同的外倾角绘制了侧向力与滑移角和垂直载荷的关系。 请注意，FSAE TTC文件仅适用于加入该联盟的团队，因此，对axampleData.DAT文件进行了较小的更改，以便不泄漏有关轮胎模型和测试回合的私人数据信息。
魔术公式轮胎模型
11-17
该文件中包含了魔术公式轮胎模型，代码全部用matlab编写，具有很高的参考价值和实际应用，非常方便！
魔术公式simulink模型
12-14
根据公式搭建了魔术公式的simulink模型，由于魔术公式的公式太多，因此用simulink里面的matlab_fuction编写的，参考相关的论文，考虑了横纵向的耦合。
基于魔术公式的轮胎车辆转向模型
09-20
压缩文件里包括了说明文件与simulink模型，包括了魔术公式的建立，车辆转向的模型，以及相平面绘制分析
铰接式车辆的横向动力学仿真提供车辆模型研究（Matlab代码实现）
Yan_she_He的博客
07-15 566
1. 定义车辆模型：首先建立铰接式车辆的几何模型，并定义车辆的基本参数，如质量、惯性特性、轮胎参数等。2. 车辆动力学模型：根据车辆几何模型和运动学原理，建立车辆的动力学模型。为了获得精确和可靠的仿真结果，建议使用专业的车辆仿真软件，并根据具体的研究目的和问题，选择合适的车辆模型和控制策略进行仿真研究。在每个时间步长内，使用车辆的输入模型、轮胎模型和操纵控制模型，计算车辆的状态和响应。4. 驱动和操纵控制模型：定义车辆的驱动和操纵控制系统模型，包括驱动力和转向力的输入模型。
Simulink根据魔术公式建立轮胎模型
Solititude的博客
12-28 3270
Simulink根据魔术公式建立轮胎模型
adams轮胎和路面的约束_魔术公式轮胎模型建立与仿真
weixin_39710106的博客
01-03 2692
要搭建一个完整的7自由度车辆操稳模型，轮胎模型是关键。本文将以魔术公式为基础，搭建一个能输出纵向力、侧向力的轮胎模型，并进行初步的仿真验证，为下一步的车辆模型搭建做好准备。期望搭建的轮胎模型能实现汽车理论中如下图所示的曲线形状：制动力系数随滑移率先增大后减小、侧偏力随侧偏角先线性增大后减小、纵向力与侧向力满足附着椭圆的约束。下面开始轮胎模型搭建。魔术公式的基本形式如下：Y可以指侧向力Fy、纵向力F...
利用MATLAB绘制万有特性图
05-26
要用MATLAB绘制万有特性图，可以按照以下步骤操作： 1. 收集发动机的性能数据，包括转速、扭矩和功率等参数。 2. 在MATLAB中导入数据，并根据需要进行预处理，如数据清洗、数据插值等。 3. 使用MATLAB绘图工具箱...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值