Compare commits

...

13 Commits

11 changed files with 1250 additions and 129 deletions

BIN
scripts/FX.fig Normal file

Binary file not shown.

BIN
scripts/FY.fig Normal file

Binary file not shown.

BIN
scripts/MX.fig Normal file

Binary file not shown.

BIN
scripts/MZ.fig Normal file

Binary file not shown.

Binary file not shown.

View File

@ -1,7 +1,8 @@
%% please open Magic Formular APP first %% please open Magic Formular APP first
clear clear
clc
n = 20; n = 20;
mkdir('similarityData');
%% %%
directoryPath = "..\sorted_data"; directoryPath = "..\sorted_data";
subdirectories = getSubdirectories(directoryPath); subdirectories = getSubdirectories(directoryPath);
@ -15,32 +16,207 @@ subdirectories_list = contents([contents.isdir] & ~ismember({contents.name}, {'.
% Extract and display the names of the subdirectories % Extract and display the names of the subdirectories
subdirectoryNames = {subdirectories_list.name}; subdirectoryNames = {subdirectories_list.name};
%% %%
figure() % figure()
% for i=1:length(subdirectories)
% fileMeasurement = fullfile("..\sorted_data", subdirectories{i}, "\cornering.mat");
% parser = tydex.parsers.FSAETTC_SI_ISO_Mat();
% measurements = parser.run(fileMeasurement);
% measurement = measurements(n);
%
% [SX,SA,FZ,IP,IA,VX,FX,FY,MZ,MY,MX,W,T] = unpack(measurement);
% FY = smoothdata(FY, "gaussian",100);
% if i < length(subdirectories)/3
% plot(SA, FY, 'LineWidth', 1)
% elseif i <= length(subdirectories)/3*2
% plot(SA, FY, '--','LineWidth', 1)
% else
% plot(SA, FY, '-.','LineWidth', 1)
% end
% SAA{i} = SA;
% FYY{i} = FY;
% hold on
% end
% legend(subdirectoryNames{:})
% title("Similarities, n=",n)
% xlabel("Slip Angle")
% ylabel("FY")
for i=1:length(subdirectories)
% SA-MX
% figure()
% for i=1:length(subdirectories)
% fileMeasurement = fullfile("..\sorted_data", subdirectories{i}, "\cornering.mat");
% parser = tydex.parsers.FSAETTC_SI_ISO_Mat();
% measurements = parser.run(fileMeasurement);
% measurement = measurements(n);
%
% [SX,SA,FZ,IP,IA,VX,FX,FY,MZ,MY,MX,W,T] = unpack(measurement);
% MX = smoothdata(MX, "gaussian",100);
% if i < length(subdirectories)/3
% plot(SA, MX, 'LineWidth', 1)
% elseif i <= length(subdirectories)/3*2
% plot(SA, MX,'--','LineWidth', 1)
% else
% plot(SA, MX, '-.','LineWidth', 1)
% end
% hold on
% end
% legend(subdirectoryNames{:})
% title("Similarities, n=",n)
% xlabel("Slip Angle")
% ylabel("MX")
% SA - MZ
% figure()
% for i=1:length(subdirectories)
% parser = tydex.parsers.FSAETTC_SI_ISO_Mat();
% fileMeasurement = fullfile("..\sorted_data", subdirectories{i}, "\cornering.mat");
% measurements = parser.run(fileMeasurement);
% measurement = measurements(n);
%
%
% [SX,SA,FZ,IP,IA,VX,FX,FY,MZ,MY,MX,W,T] = unpack(measurement);
% MZ = smoothdata(MZ, "gaussian",100);
% if i < length(subdirectories)/3
% plot(SA, MZ, 'LineWidth', 1)
% elseif i <= length(subdirectories)/3*2
% plot(SA, MZ,'--','LineWidth', 1)
% else
% plot(SA, MZ, '-.','LineWidth', 1)
% end
% SAA{i} = SA;
% MZZ{i} = MZ;
% hold on
% end
% legend(subdirectoryNames{:})
% title("Similarities, n=",n)
% xlabel("Slip Angle")
% ylabel("MZ")
% for n=1:1
% figure()
% j = 1;
% for i=8:8%length(subdirectories)
% fileMeasurement = fullfile("..\sorted_data", subdirectories{i}, "\cornering.mat");
% fileMeasurement1 = fullfile("..\sorted_data", subdirectories{i}, "\drivebrake.mat");
% parser = tydex.parsers.FSAETTC_SI_ISO_Mat();
% try
% measurements = parser.run(fileMeasurement);
% measurement = measurements(n);
%
% measurements1 = parser.run(fileMeasurement1);
% measurement1 = measurements1(n);
%
% [SX,SA,FZ,IP,IA,VX,FX,FY,MZ,MY,MX,W,T] = unpack(measurement);
% [SX1,SA1,FZ1,IP1,IA1,VX1,FX1,FY1,MZ1,MY1,MX1,W1,T1] = unpack(measurement1);
%
% MX = smoothdata(MX, "gaussian");
% if i < length(subdirectories)/3
% plot(SA, FY, 'LineWidth', 1)
% elseif i <= length(subdirectories)/3*2
% plot(SA, FY, '--','LineWidth', 1)
% else
% plot(SA, FY, '-.','LineWidth', 1)
% end
% lengends{j} = subdirectoryNames{i};
% j = j+1;
% hold on
% catch
% disp(fileMeasurement)
% end
% end
% legend(lengends{:})
% title("Similarities, n=",n)
% xlabel("SA")
% ylabel("MZ")
% end
%% prove FY are the same
for n=1:10:60
figure()
for i=9:9%length(subdirectories)
fileMeasurement = fullfile("..\sorted_data", subdirectories{i}, "\cornering.mat"); fileMeasurement = fullfile("..\sorted_data", subdirectories{i}, "\cornering.mat");
fileMeasurement1 = fullfile("..\sorted_data", subdirectories{i}, "\drivebrake.mat");
parser = tydex.parsers.FSAETTC_SI_ISO_Mat(); parser = tydex.parsers.FSAETTC_SI_ISO_Mat();
try
measurements = parser.run(fileMeasurement); measurements = parser.run(fileMeasurement);
measurement = measurements(n); measurement = measurements(n);
measurements1 = parser.run(fileMeasurement1);
measurement1 = measurements1(n);
[SX,SA,FZ,IP,IA,VX,FX,FY,MZ,MY,MX,W,T] = unpack(measurement); [SX,SA,FZ,IP,IA,VX,FX,FY,MZ,MY,MX,W,T] = unpack(measurement);
FY = smoothdata(FY, "gaussian",100); [SX1,SA1,FZ1,IP1,IA1,VX1,FX1,FY1,MZ1,MY1,MX1,W1,T1] = unpack(measurement1);
if i < length(subdirectories)/3
MX = smoothdata(MX, "gaussian");
plot(SA, FY, 'LineWidth', 1) plot(SA, FY, 'LineWidth', 1)
elseif i <= length(subdirectories)/3*2
plot(SA, FY, '--','LineWidth', 1)
else
plot(SA, FY, '-.','LineWidth', 1)
end
SAA{i} = SA;
FYY{i} = FY;
hold on hold on
plot(SA1, FY1, 'LineWidth', 1)
legend("Cornering", "drivebrake")
catch
disp(fileMeasurement)
end
end end
legend(subdirectoryNames{:})
title("Similarities, n=",n) title("Similarities, n=",n)
xlabel("Slip Angle") xlabel("SA")
ylabel("FY") ylabel("FY")
end
%% without plot %% without plot
% mkdir('similarityDataFY');
% for n=1:60
% for i=1:length(subdirectories)
% fileMeasurement = fullfile("..\sorted_data", subdirectories{i}, "\cornering.mat");
% parser = tydex.parsers.FSAETTC_SI_ISO_Mat();
% measurements = parser.run(fileMeasurement);
% measurement = measurements(n);
%
% [SX,SA,FZ,IP,IA,VX,FX,FY,MZ,MY,MX,W,T] = unpack(measurement);
% FY = smoothdata(FY, "gaussian",100);
% SAA{i} = SA;
% FYY{i} = FY;
% end
% file_name = ['similarityDataFY\nEqu' num2str(n) '.mat'];
% save(file_name, 'SAA', 'FYY');
% end
% mkdir('similarityDataFX');
% for n=1:60
% for i=1:length(subdirectories)
% fileMeasurement = fullfile("..\sorted_data", subdirectories{i}, "\cornering.mat");
% parser = tydex.parsers.FSAETTC_SI_ISO_Mat();
% measurements = parser.run(fileMeasurement);
% measurement = measurements(n);
%
% [SX,SA,FZ,IP,IA,VX,FX,FY,MZ,MY,MX,W,T] = unpack(measurement);
% FX = smoothdata(FX, "gaussian",100);
% SAA{i} = SA;
% FXX{i} = FX;
% end
% file_name = ['similarityDataFX\nEqu' num2str(n) '.mat'];
% save(file_name, 'SAA', 'FXX');
% end
%
% mkdir("similarityDataMZ")
% for n=1:60
% for i=1:length(subdirectories)
% fileMeasurement = fullfile("..\sorted_data", subdirectories{i}, "\cornering.mat");
% parser = tydex.parsers.FSAETTC_SI_ISO_Mat();
% measurements = parser.run(fileMeasurement);
% measurement = measurements(n);
%
% [SX,SA,FZ,IP,IA,VX,FX,FY,MZ,MY,MX,W,T] = unpack(measurement);
% MZ = smoothdata(MZ, "gaussian",100);
% SAA{i} = SA;
% MZZ{i} = MZ;
% end
% file_name = ['similarityDataMZ\nEqu' num2str(n) '.mat'];
% save(file_name, 'SAA', 'MZZ');
% end
mkdir("similarityDataMX")
for n=1:60 for n=1:60
for i=1:length(subdirectories) for i=1:length(subdirectories)
fileMeasurement = fullfile("..\sorted_data", subdirectories{i}, "\cornering.mat"); fileMeasurement = fullfile("..\sorted_data", subdirectories{i}, "\cornering.mat");
@ -49,10 +225,10 @@ for n=1:60
measurement = measurements(n); measurement = measurements(n);
[SX,SA,FZ,IP,IA,VX,FX,FY,MZ,MY,MX,W,T] = unpack(measurement); [SX,SA,FZ,IP,IA,VX,FX,FY,MZ,MY,MX,W,T] = unpack(measurement);
FY = smoothdata(FY, "gaussian",100); MX = smoothdata(MX, "gaussian",100);
SAA{i} = SA; SAA{i} = SA;
FYY{i} = FY; MXX{i} = MX;
end end
file_name = ['similarityData\nEqu' num2str(n) '.mat']; file_name = ['similarityDataMX\nEqu' num2str(n) '.mat'];
save(file_name, 'SAA', 'FYY'); save(file_name, 'SAA', 'MXX');
end end

View File

@ -1,6 +1,7 @@
clear; clear;
clc; clc;
% Use the dir function to get information about files and directories
% % FY
directoryPath = ".\similarityData"; directoryPath = ".\similarityData";
% Get a list of all files in the directory % Get a list of all files in the directory
@ -29,7 +30,7 @@ for tire=1:length(SAA)
end end
residuals_tire{tire} = residuals_one_tire; residuals_tire{tire} = residuals_one_tire;
end end
%% mean % mean
for tire=1:length(SAA) for tire=1:length(SAA)
num = 0; num = 0;
% each tire under 60 measurement conditions % each tire under 60 measurement conditions
@ -42,11 +43,11 @@ for tire=1:length(SAA)
end end
residuals_mean{tire} = sum(residuals_tire{tire})/num; residuals_mean{tire} = sum(residuals_tire{tire})/num;
end end
%% median % median
for tire=1:length(SAA) for tire=1:length(SAA)
residuals_median{tire} = median(residuals_tire{tire}); residuals_median{tire} = median(residuals_tire{tire});
end end
%% mode % mode
for tire=1:length(SAA) for tire=1:length(SAA)
residuals_mode{tire} = mode(residuals_tire{tire}); residuals_mode{tire} = mode(residuals_tire{tire});
end end
@ -64,6 +65,7 @@ ylabel("in 1")
xticks(idx) xticks(idx)
grid on grid on
subplot(312) subplot(312)
y = cell2mat(residuals_median); y = cell2mat(residuals_median);
plot(idx, y, "-o") plot(idx, y, "-o")
@ -82,3 +84,263 @@ xlabel("Order of the tires(8 is our traget tire)")
ylabel("in 1") ylabel("in 1")
xticks(idx) xticks(idx)
grid on grid on
sgtitle("FY")
%% FX
directoryPath = ".\similarityDataFX";
% Get a list of all files in the directory
files = dir(fullfile(directoryPath, '*.mat')); % You can adjust the file extension as needed
% Extract file names
fileNames = {files.name};
for n=1:length(fileNames)
% fit the curve of 43075 , 8 inch
filePath = fullfile(directoryPath, fileNames(n));
load(filePath);
[fitresult, gof] = createFit(SAA{8},FXX{8});
% Evaluate the fitted curves
for tire=1:length(SAA)
y_fit = fitresult(SAA{tire});
% Calculate residuals
residuals{tire}{n} = abs(FXX{tire} - fitresult(SAA{tire}));
end
end
for tire=1:length(SAA)
residuals_one_tire = [];
for n=1:length(fileNames)
residuals_one_tire = [residuals_one_tire; residuals{tire}{n}];
end
residuals_tire{tire} = residuals_one_tire;
end
% mean
for tire=1:length(SAA)
num = 0;
% each tire under 60 measurement conditions
for i=1:length(residuals_tire{tire})
if not(isnan(residuals_tire{tire}(i)))
num = num + 1;
else
residuals_tire{tire}(i) = 0;
end
end
residuals_mean{tire} = sum(residuals_tire{tire})/num;
end
% median
for tire=1:length(SAA)
residuals_median{tire} = median(residuals_tire{tire});
end
% mode
for tire=1:length(SAA)
residuals_mode{tire} = mode(residuals_tire{tire});
end
% plots
idx = linspace(1,length(SAA),length(SAA));
figure()
subplot(311)
y = cell2mat(residuals_mean);
plot(idx, y, "-*")
title("Average Absolute Residuials")
xlabel("Order of the tires(8 is our traget tire)")
ylabel("in 1")
xticks(idx)
grid on
subplot(312)
y = cell2mat(residuals_median);
plot(idx, y, "-o")
title("Median of the Absolute Residuials")
xlabel("Order of the tires(8 is our traget tire)")
ylabel("in 1")
xticks(idx)
grid on
subplot(313)
y = cell2mat(residuals_mode);
plot(idx, y, "-o")
title("Mode of the Absolute Residuials")
xlabel("Order of the tires(8 is our traget tire)")
ylabel("in 1")
xticks(idx)
grid on
sgtitle("FX")
%% Mx
directoryPath = ".\similarityDataMX";
% Get a list of all files in the directory
files = dir(fullfile(directoryPath, '*.mat')); % You can adjust the file extension as needed
% Extract file names
fileNames = {files.name};
for n=1:length(fileNames)
% fit the curve of 43075 , 8 inch
filePath = fullfile(directoryPath, fileNames(n));
load(filePath);
[fitresult, gof] = createFit(SAA{8},MXX{8});
% Evaluate the fitted curves
for tire=1:length(SAA)
y_fit = fitresult(SAA{tire});
% Calculate residuals
residuals{tire}{n} = abs(MXX{tire} - fitresult(SAA{tire}));
end
end
for tire=1:length(SAA)
residuals_one_tire = [];
for n=1:length(fileNames)
residuals_one_tire = [residuals_one_tire; residuals{tire}{n}];
end
residuals_tire{tire} = residuals_one_tire;
end
% mean
for tire=1:length(SAA)
num = 0;
% each tire under 60 measurement conditions
for i=1:length(residuals_tire{tire})
if not(isnan(residuals_tire{tire}(i)))
num = num + 1;
else
residuals_tire{tire}(i) = 0;
end
end
residuals_mean{tire} = sum(residuals_tire{tire})/num;
end
% median
for tire=1:length(SAA)
residuals_median{tire} = median(residuals_tire{tire});
end
% mode
for tire=1:length(SAA)
residuals_mode{tire} = mode(residuals_tire{tire});
end
% plots
idx = linspace(1,length(SAA),length(SAA));
figure()
subplot(311)
y = cell2mat(residuals_mean);
plot(idx, y, "-*")
title("Average Absolute Residuials")
xlabel("Order of the tires(8 is our traget tire)")
ylabel("in 1")
xticks(idx)
grid on
subplot(312)
y = cell2mat(residuals_median);
plot(idx, y, "-o")
title("Median of the Absolute Residuials")
xlabel("Order of the tires(8 is our traget tire)")
ylabel("in 1")
xticks(idx)
grid on
subplot(313)
y = cell2mat(residuals_mode);
plot(idx, y, "-o")
title("Mode of the Absolute Residuials")
xlabel("Order of the tires(8 is our traget tire)")
ylabel("in 1")
xticks(idx)
grid on
sgtitle("MX")
%% Mz
directoryPath = ".\similarityDataMZ";
% Get a list of all files in the directory
files = dir(fullfile(directoryPath, '*.mat')); % You can adjust the file extension as needed
% Extract file names
fileNames = {files.name};
for n=1:length(fileNames)
% fit the curve of 43075 , 8 inch
filePath = fullfile(directoryPath, fileNames(n));
load(filePath);
[fitresult, gof] = createFit(SAA{8},MZZ{8});
% Evaluate the fitted curves
for tire=1:length(SAA)
y_fit = fitresult(SAA{tire});
% Calculate residuals
residuals{tire}{n} = abs(MZZ{tire} - fitresult(SAA{tire}));
end
end
for tire=1:length(SAA)
residuals_one_tire = [];
for n=1:length(fileNames)
residuals_one_tire = [residuals_one_tire; residuals{tire}{n}];
end
residuals_tire{tire} = residuals_one_tire;
end
% mean
for tire=1:length(SAA)
num = 0;
% each tire under 60 measurement conditions
for i=1:length(residuals_tire{tire})
if not(isnan(residuals_tire{tire}(i)))
num = num + 1;
else
residuals_tire{tire}(i) = 0;
end
end
residuals_mean{tire} = sum(residuals_tire{tire})/num;
end
% median
for tire=1:length(SAA)
residuals_median{tire} = median(residuals_tire{tire});
end
% mode
for tire=1:length(SAA)
residuals_mode{tire} = mode(residuals_tire{tire});
end
% plots
idx = linspace(1,length(SAA),length(SAA));
figure()
subplot(311)
y = cell2mat(residuals_mean);
plot(idx, y, "-*")
title("Average Absolute Residuials")
xlabel("Order of the tires(8 is our traget tire)")
ylabel("in 1")
grid on
xticks(idx)
subplot(312)
y = cell2mat(residuals_median);
plot(idx, y, "-o")
title("Median of the Absolute Residuials")
xlabel("Order of the tires(8 is our traget tire)")
ylabel("in 1")
xticks(idx)
grid on
subplot(313)
y = cell2mat(residuals_mode);
plot(idx, y, "-o")
title("Mode of the Absolute Residuials")
xlabel("Order of the tires(8 is our traget tire)")
ylabel("in 1")
xticks(idx)
grid on
sgtitle("MZ")

View File

@ -0,0 +1,229 @@
[MDI_HEADER]
FILE_TYPE = 'tir'
FILE_VERSION = 3
FILE_FORMAT = 'ASCII'
$-----------------------------------------------------------------units
[UNITS]
LENGTH = 'meter'
FORCE = 'newton'
ANGLE = 'radians'
MASS = 'kg'
TIME = 'second'
$-----------------------------------------------------------------model
[MODEL]
FITTYP = 61
TYRESIDE = 'LEFT'
LONGVL = 10
VXLOW = 1
$-------------------------------------------------------------dimension
[DIMENSION]
UNLOADED_RADIUS = 0.20574
WIDTH = 0.1905
ASPECT_RATIO = 0.8267
RIM_RADIUS = 0.127
RIM_WIDTH = 0.2032
$--------------------------------------------------operating conditions
[OPERATING_CONDITIONS]
NOMPRES = 100000
$---------------------------------------------------------------inertia
[INERTIA]
MASS = kg
$--------------------------------------------------------------vertical
[VERTICAL]
FNOMIN = 1500
$------------------------------------------------------------structural
[STRUCTURAL]
$---------------------------------------------------------contact patch
[CONTACT_PATCH]
$----------------------------------------------inflation pressure range
[INFLATION_PRESSURE_RANGE]
$--------------------------------------------------vertical force range
[VERTICAL_FORCE_RANGE]
$-------------------------------------------------------long slip range
[LONG_SLIP_RANGE]
$------------------------------------------------------slip angle range
[SLIP_ANGLE_RANGE]
$-----------------------------------------------inclination angle range
[INCLINATION_ANGLE_RANGE]
$--------------------------------------------------scaling coefficients
[SCALING_COEFFICIENTS]
LFZO = 1
LCX = 1
LMUX = 1
LEX = 1
LKX = 1
LHX = 1
LVX = 1
LCY = 1
LMUY = 1
LEY = 1
LKY = 1
LHY = 1
LVY = 1
LTR = 1
LRES = 1
LXAL = 1
LYKA = 1
LVYKA = 1
LS = 1
LKYC = 1
LKZC = 1
LVMX = 1
LMX = 1
LMY = 1
LMP = 1
$---------------------------------------------longitudinal coefficients
[LONGITUDINAL_COEFFICIENTS]
PCX1 = 1.5
PDX1 = 2.1186
PDX2 = -0.13379
PDX3 = 11.9723
PEX1 = -0.27146
PEX2 = 0.40109
PEX3 = 0.6
PEX4 = -0.9
PKX1 = 35.4142
PKX2 = 20.0246
PKX3 = -0.97661
PHX1 = -0.0055756
PHX2 = 0.0002412
PVX1 = 0.058717
PVX2 = -0.033067
PPX1 = -2.4576
PPX2 = -4.043
PPX3 = -0.28527
PPX4 = 0.36685
RBX1 = 14.6533
RBX2 = 15.583
RBX3 = 48.2717
RCX1 = 1.2
REX1 = 0.50883
REX2 = 0.24823
RHX1 = 0.012901
$----------------------------------------------overturning coefficients
[OVERTURNING_COEFFICIENTS]
QSX1 = -0.023509
QSX2 = 2.0014
QSX3 = 0.11274
QSX4 = 0.071804
QSX5 = 0
QSX6 = 0
QSX7 = 0.024372
QSX8 = 0
QSX9 = 0
QSX10 = 0
QSX11 = 0
QSX12 = 0
QSX13 = 0
QSX14 = 0
PPMX1 = 0.46712
$--------------------------------------------------lateral coefficients
[LATERAL_COEFFICIENTS]
PCY1 = 1.5
PDY1 = 2.2152
PDY2 = -0.24006
PDY3 = 8.7883
PEY1 = 0.22956
PEY2 = 0.23261
PEY3 = 1.2231
PEY4 = -73.217
PEY5 = -368.3453
PKY1 = -27.4287
PKY2 = 1.3665
PKY3 = 1
PKY4 = 2
PKY5 = 16.4211
PKY6 = -5.0355
PKY7 = -2.3904
PHY1 = 0.0029497
PHY2 = 0.0030317
PVY1 = 0.097172
PVY2 = 0.1046
PVY3 = -1.9398
PVY4 = -0.060559
PPY1 = 0.62414
PPY2 = 1.4562
PPY3 = -0.21694
PPY4 = 0.20207
PPY5 = -0.2294
RBY1 = 15.2511
RBY2 = 4.6851
RBY3 = -0.15
RBY4 = 5.524
RCY1 = 0.85866
REY1 = -0.19774
REY2 = 5
RHY1 = 0.007696
RHY2 = 0.019597
RVY1 = 0
RVY2 = 0
RVY3 = 0
RVY4 = 0
RVY5 = 0
RVY6 = 0
$--------------------------------------------------rolling coefficients
[ROLLING_COEFFICIENTS]
QSY1 = 0
QSY2 = 0
QSY3 = 0
QSY4 = 0
QSY5 = 0
QSY6 = 0
QSY7 = 0
QSY8 = 0
$-------------------------------------------------aligning coefficients
[ALIGNING_COEFFICIENTS]
QBZ1 = 2.2492
QBZ2 = -2.0829
QBZ3 = -1.6133
QBZ4 = 113.4182
QBZ5 = -113.1977
QBZ9 = -2.2113
QBZ10 = -0.10713
QCZ1 = 3.1732
QDZ1 = 0.2302
QDZ2 = -0.089693
QDZ3 = 2.4601
QDZ4 = -34.0481
QDZ6 = 0.015409
QDZ7 = -0.0076995
QDZ8 = -0.084439
QDZ9 = 9.0558
QDZ10 = -8.1707
QDZ11 = -26.6613
QEZ1 = 0.58479
QEZ2 = -2.5778
QEZ3 = -3.8586
QEZ4 = 0.015591
QEZ5 = -0.090861
QHZ1 = 0.0019082
QHZ2 = -0.0033306
QHZ3 = -0.02611
QHZ4 = 1.3539
PPZ1 = -2.2988
PPZ2 = 0.24851
SSZ1 = -0.031619
SSZ2 = -0.078744
SSZ3 = 0
SSZ4 = 0
$-------------------------------------------------turnslip coefficients
[TURNSLIP_COEFFICIENTS]
PDXP1 = 0
PDXP2 = 0
PDXP3 = 0
PKYP1 = 0
PDYP1 = 0
PDYP2 = 0
PDYP3 = 0
PDYP4 = 0
PHYP1 = 0
PHYP2 = 0
PHYP3 = 0
PHYP4 = 0
PECP1 = 0
PECP2 = 0
QDTP1 = 0
QCRP1 = 0
QCRP2 = 0
QBRP1 = 0
QDRP1 = 0

View File

@ -0,0 +1,228 @@
[MDI_HEADER]
FILE_TYPE = 'tir'
FILE_VERSION = 3
FILE_FORMAT = 'ASCII'
$-----------------------------------------------------------------units
[UNITS]
LENGTH = 'meter'
FORCE = 'newton'
ANGLE = 'radians'
MASS = 'kg'
TIME = 'second'
$-----------------------------------------------------------------model
[MODEL]
FITTYP = 61
TYRESIDE = 'LEFT'
LONGVL = 10
VXLOW = 1
$-------------------------------------------------------------dimension
[DIMENSION]
UNLOADED_RADIUS = 0.20574
WIDTH = 0.2032
RIM_RADIUS = 0.254
RIM_WIDTH = 0.2032
$--------------------------------------------------operating conditions
[OPERATING_CONDITIONS]
NOMPRES = 97000
$---------------------------------------------------------------inertia
[INERTIA]
MASS = 3.6287
$--------------------------------------------------------------vertical
[VERTICAL]
FNOMIN = 1140
$------------------------------------------------------------structural
[STRUCTURAL]
$---------------------------------------------------------contact patch
[CONTACT_PATCH]
$----------------------------------------------inflation pressure range
[INFLATION_PRESSURE_RANGE]
$--------------------------------------------------vertical force range
[VERTICAL_FORCE_RANGE]
$-------------------------------------------------------long slip range
[LONG_SLIP_RANGE]
$------------------------------------------------------slip angle range
[SLIP_ANGLE_RANGE]
$-----------------------------------------------inclination angle range
[INCLINATION_ANGLE_RANGE]
$--------------------------------------------------scaling coefficients
[SCALING_COEFFICIENTS]
LFZO = 1
LCX = 1
LMUX = 1
LEX = 1
LKX = 1
LHX = 1
LVX = 1
LCY = 1
LMUY = 1
LEY = 1
LKY = 1
LHY = 1
LVY = 1
LTR = 1
LRES = 1
LXAL = 1
LYKA = 1
LVYKA = 1
LS = 1
LKYC = 1
LKZC = 1
LVMX = 1
LMX = 1
LMY = 1
LMP = 1
$---------------------------------------------longitudinal coefficients
[LONGITUDINAL_COEFFICIENTS]
PCX1 = 1.5
PDX1 = 2.0056
PDX2 = -0.14219
PDX3 = 3.9302
PEX1 = -0.44146
PEX2 = -0.60999
PEX3 = 0.6
PEX4 = -0.9
PKX1 = 30.6895
PKX2 = 22.1344
PKX3 = -1.2545
PHX1 = -0.0039672
PHX2 = -0.00099913
PVX1 = -0.0019729
PVX2 = -0.085842
PPX1 = -1.9846
PPX2 = -2.39
PPX3 = -0.6612
PPX4 = -0.70183
RBX1 = 20.8165
RBX2 = 16.5244
RBX3 = 719.7481
RCX1 = 1.04
REX1 = 0.67305
REX2 = 0.62374
RHX1 = -0.0025848
$----------------------------------------------overturning coefficients
[OVERTURNING_COEFFICIENTS]
QSX1 = -0.016561
QSX2 = 2.1865
QSX3 = 0.081194
QSX4 = 0.44488
QSX5 = 0
QSX6 = 0
QSX7 = 0.019222
QSX8 = 0
QSX9 = 0
QSX10 = 0
QSX11 = 0
QSX12 = 0
QSX13 = 0
QSX14 = 0
PPMX1 = 0.37295
$--------------------------------------------------lateral coefficients
[LATERAL_COEFFICIENTS]
PCY1 = 1.8586
PDY1 = 2.258
PDY2 = -0.32829
PDY3 = 9.9601
PEY1 = 0.64364
PEY2 = -0.55542
PEY3 = -0.10996
PEY4 = -5.3879
PEY5 = -45.5915
PKY1 = -41.2181
PKY2 = 2.0528
PKY3 = 0.72559
PKY4 = 2
PKY5 = 36.8632
PKY6 = -4.0914
PKY7 = -2.0202
PHY1 = 0.00027403
PHY2 = -2.2872e-05
PVY1 = 0.048602
PVY2 = 0.030584
PVY3 = -0.26546
PVY4 = 1.9087
PPY1 = 0.50301
PPY2 = 1.2958
PPY3 = -0.15219
PPY4 = 0.093692
PPY5 = -1.1309
RBY1 = 38.03
RBY2 = 14.4955
RBY3 = 0.039486
RBY4 = 49.0927
RCY1 = 0.91003
REY1 = 0.74082
REY2 = 0.1359
RHY1 = 0.0025767
RHY2 = 0.0059086
RVY1 = 0
RVY2 = 0
RVY3 = 0
RVY4 = 0
RVY5 = 0
RVY6 = 0
$--------------------------------------------------rolling coefficients
[ROLLING_COEFFICIENTS]
QSY1 = 0
QSY2 = 0
QSY3 = 0
QSY4 = 0
QSY5 = 0
QSY6 = 0
QSY7 = 0
QSY8 = 0
$-------------------------------------------------aligning coefficients
[ALIGNING_COEFFICIENTS]
QBZ1 = 0.42286
QBZ2 = -0.14241
QBZ3 = -0.14767
QBZ4 = 27.0482
QBZ5 = -27.2817
QBZ9 = 0.058045
QBZ10 = 0.0026325
QCZ1 = 17.3533
QDZ1 = 0.14691
QDZ2 = -0.059137
QDZ3 = 1.3754
QDZ4 = -27.248
QDZ6 = 0.01906
QDZ7 = -0.009968
QDZ8 = -2.7485
QDZ9 = 6.0439
QDZ10 = 5.4202
QDZ11 = -24.0105
QEZ1 = -1.6389
QEZ2 = -9.0585
QEZ3 = -11.2643
QEZ4 = 2.2186
QEZ5 = 0.22602
QHZ1 = 0.0028312
QHZ2 = -0.013895
QHZ3 = -0.26076
QHZ4 = 1.0186
PPZ1 = -2.7245
PPZ2 = -0.74267
SSZ1 = -0.048316
SSZ2 = -0.091767
SSZ3 = 0
SSZ4 = 0
$-------------------------------------------------turnslip coefficients
[TURNSLIP_COEFFICIENTS]
PDXP1 = 0
PDXP2 = 0
PDXP3 = 0
PKYP1 = 0
PDYP1 = 0
PDYP2 = 0
PDYP3 = 0
PDYP4 = 0
PHYP1 = 0
PHYP2 = 0
PHYP3 = 0
PHYP4 = 0
PECP1 = 0
PECP2 = 0
QDTP1 = 0
QCRP1 = 0
QCRP2 = 0
QBRP1 = 0
QDRP1 = 0

View File

@ -17,15 +17,16 @@ LONGVL = 10
VXLOW = 1 VXLOW = 1
$-------------------------------------------------------------dimension $-------------------------------------------------------------dimension
[DIMENSION] [DIMENSION]
UNLOADED_RADIUS = 0.2286 UNLOADED_RADIUS = 0.20574
$--------------------------------------------------operating conditions $--------------------------------------------------operating conditions
[OPERATING_CONDITIONS] [OPERATING_CONDITIONS]
NOMPRES = 100000 NOMPRES = 99000
$---------------------------------------------------------------inertia $---------------------------------------------------------------inertia
[INERTIA] [INERTIA]
MASS = kg
$--------------------------------------------------------------vertical $--------------------------------------------------------------vertical
[VERTICAL] [VERTICAL]
FNOMIN = 1500 FNOMIN = 1560
$------------------------------------------------------------structural $------------------------------------------------------------structural
[STRUCTURAL] [STRUCTURAL]
$---------------------------------------------------------contact patch $---------------------------------------------------------contact patch
@ -70,40 +71,40 @@ LMP = 1
$---------------------------------------------longitudinal coefficients $---------------------------------------------longitudinal coefficients
[LONGITUDINAL_COEFFICIENTS] [LONGITUDINAL_COEFFICIENTS]
PCX1 = 1.5 PCX1 = 1.5
PDX1 = 1.9208 PDX1 = 2.0648
PDX2 = -0.17852 PDX2 = -0.15342
PDX3 = 3.9292 PDX3 = 12.1289
PEX1 = -0.6473 PEX1 = -0.17784
PEX2 = -0.67795 PEX2 = 0.563
PEX3 = 0.6 PEX3 = 0.6
PEX4 = -0.9 PEX4 = -0.9
PKX1 = 23.6691 PKX1 = 35.1456
PKX2 = 18.5248 PKX2 = 24.3339
PKX3 = -1.6735 PKX3 = -1.2273
PHX1 = -0.0044404 PHX1 = -0.0050602
PHX2 = -0.0017213 PHX2 = 0.0013983
PVX1 = -0.024192 PVX1 = 0.034156
PVX2 = -0.1 PVX2 = -0.1
PPX1 = -2.3463 PPX1 = -1.7796
PPX2 = -2.7112 PPX2 = -2.9784
PPX3 = -0.74204 PPX3 = -0.3492
PPX4 = -0.7617 PPX4 = 0.38929
RBX1 = 5.99 RBX1 = 11.1213
RBX2 = 5.99 RBX2 = 12.6753
RBX3 = 0 RBX3 = 5.0211
RCX1 = 1 RCX1 = 1.1591
REX1 = -1 REX1 = -0.06655
REX2 = -0.1 REX2 = 0.055878
RHX1 = 0 RHX1 = 0.032551
$----------------------------------------------overturning coefficients $----------------------------------------------overturning coefficients
[OVERTURNING_COEFFICIENTS] [OVERTURNING_COEFFICIENTS]
QSX1 = 0 QSX1 = -0.011071
QSX2 = 1.5 QSX2 = 1.9673
QSX3 = 0.1 QSX3 = 0.11823
QSX4 = 0.1 QSX4 = 0.066517
QSX5 = 0 QSX5 = 0
QSX6 = 0 QSX6 = 0
QSX7 = 0.1 QSX7 = 0.05706
QSX8 = 0 QSX8 = 0
QSX9 = 0 QSX9 = 0
QSX10 = 0 QSX10 = 0
@ -111,45 +112,45 @@ QSX11 = 0
QSX12 = 0 QSX12 = 0
QSX13 = 0 QSX13 = 0
QSX14 = 0 QSX14 = 0
PPMX1 = 0 PPMX1 = 0.36777
$--------------------------------------------------lateral coefficients $--------------------------------------------------lateral coefficients
[LATERAL_COEFFICIENTS] [LATERAL_COEFFICIENTS]
PCY1 = 1.5 PCY1 = 1.6675
PDY1 = 2.4789 PDY1 = 2.1768
PDY2 = -0.5 PDY2 = -0.48644
PDY3 = 10 PDY3 = 5.9709
PEY1 = 1.2322 PEY1 = 0.51141
PEY2 = 0.47308 PEY2 = -0.012972
PEY3 = -0.3052 PEY3 = -0.050898
PEY4 = -7.8434 PEY4 = -4.5261
PEY5 = -115.0383 PEY5 = -70.188
PKY1 = -21.5408 PKY1 = -58.0824
PKY2 = 1.2661 PKY2 = 2.1273
PKY3 = -0.59067 PKY3 = 1
PKY4 = 2 PKY4 = 2
PKY5 = -25.9435 PKY5 = 34.5112
PKY6 = -2.5098 PKY6 = -5.7422
PKY7 = 0.61093 PKY7 = -3.8244
PHY1 = -0.0083677 PHY1 = 0.0019699
PHY2 = -0.0052997 PHY2 = 0.00075738
PVY1 = -0.16592 PVY1 = 0.051696
PVY2 = 0.071875 PVY2 = 0.014416
PVY3 = -0.74422 PVY3 = 1.0271
PVY4 = -3.5858 PVY4 = 1.2347
PPY1 = 0.32063 PPY1 = 0.82385
PPY2 = 1.3726 PPY2 = 1.4074
PPY3 = -0.2143 PPY3 = -0.17124
PPY4 = -0.032422 PPY4 = -0.11678
PPY5 = -1.5161 PPY5 = -0.60425
RBY1 = 5 RBY1 = 27.8314
RBY2 = 2.99 RBY2 = 18.3875
RBY3 = 0.02 RBY3 = 0.030113
RBY4 = 0 RBY4 = -90
RCY1 = 1 RCY1 = 0.98399
REY1 = -0.1 REY1 = 0.66603
REY2 = 0.1 REY2 = 0.42965
RHY1 = 0 RHY1 = 0.0027678
RHY2 = 0 RHY2 = 0.0092365
RVY1 = 0 RVY1 = 0
RVY2 = 0 RVY2 = 0
RVY3 = 0 RVY3 = 0
@ -168,37 +169,37 @@ QSY7 = 0
QSY8 = 0 QSY8 = 0
$-------------------------------------------------aligning coefficients $-------------------------------------------------aligning coefficients
[ALIGNING_COEFFICIENTS] [ALIGNING_COEFFICIENTS]
QBZ1 = 0.10008 QBZ1 = 1.167
QBZ2 = 0.00023246 QBZ2 = -0.12363
QBZ3 = 0.00029086 QBZ3 = -0.024724
QBZ4 = -0.0031101 QBZ4 = -178.5697
QBZ5 = -0.0031102 QBZ5 = 178.4563
QBZ9 = 1.5088 QBZ9 = -15.5984
QBZ10 = -0.19344 QBZ10 = -0.42485
QCZ1 = 0.11987 QCZ1 = 7.1867
QDZ1 = 0.1152 QDZ1 = 0.19446
QDZ2 = 0.0047181 QDZ2 = -0.073114
QDZ3 = 0.96498 QDZ3 = 1.4741
QDZ4 = -0.73585 QDZ4 = -33.4633
QDZ6 = -0.0036203 QDZ6 = 0.018614
QDZ7 = -0.050338 QDZ7 = -0.0071083
QDZ8 = -0.084844 QDZ8 = -1.8286
QDZ9 = 5.8235 QDZ9 = 2.4371
QDZ10 = 1.612 QDZ10 = 0.30258
QDZ11 = -0.47189 QDZ11 = -2.7828
QEZ1 = -0.025293 QEZ1 = 0.2928
QEZ2 = -0.096812 QEZ2 = -10.1916
QEZ3 = -0.090171 QEZ3 = -31.1995
QEZ4 = -0.0044756 QEZ4 = -0.30557
QEZ5 = -0.0044755 QEZ5 = 27.0423
QHZ1 = 0.01593 QHZ1 = 0.0039846
QHZ2 = -0.0020194 QHZ2 = -0.0045528
QHZ3 = -0.002358 QHZ3 = 0.085773
QHZ4 = 0.0013057 QHZ4 = -0.02284
PPZ1 = 0 PPZ1 = -3.0722
PPZ2 = 0 PPZ2 = -2.7578
SSZ1 = 0 SSZ1 = -0.032766
SSZ2 = 0 SSZ2 = -0.066404
SSZ3 = 0 SSZ3 = 0
SSZ4 = 0 SSZ4 = 0
$-------------------------------------------------turnslip coefficients $-------------------------------------------------turnslip coefficients

View File

@ -0,0 +1,225 @@
[MDI_HEADER]
FILE_TYPE = 'tir'
FILE_VERSION = 3
FILE_FORMAT = 'ASCII'
$-----------------------------------------------------------------units
[UNITS]
LENGTH = 'meter'
FORCE = 'newton'
ANGLE = 'radians'
MASS = 'kg'
TIME = 'second'
$-----------------------------------------------------------------model
[MODEL]
FITTYP = 61
TYRESIDE = 'LEFT'
LONGVL = 10
VXLOW = 1
$-------------------------------------------------------------dimension
[DIMENSION]
UNLOADED_RADIUS = 0.2667
$--------------------------------------------------operating conditions
[OPERATING_CONDITIONS]
NOMPRES = 99000
$---------------------------------------------------------------inertia
[INERTIA]
MASS = kg
$--------------------------------------------------------------vertical
[VERTICAL]
FNOMIN = 1560
$------------------------------------------------------------structural
[STRUCTURAL]
$---------------------------------------------------------contact patch
[CONTACT_PATCH]
$----------------------------------------------inflation pressure range
[INFLATION_PRESSURE_RANGE]
$--------------------------------------------------vertical force range
[VERTICAL_FORCE_RANGE]
$-------------------------------------------------------long slip range
[LONG_SLIP_RANGE]
$------------------------------------------------------slip angle range
[SLIP_ANGLE_RANGE]
$-----------------------------------------------inclination angle range
[INCLINATION_ANGLE_RANGE]
$--------------------------------------------------scaling coefficients
[SCALING_COEFFICIENTS]
LFZO = 1
LCX = 1
LMUX = 1
LEX = 1
LKX = 1
LHX = 1
LVX = 1
LCY = 1
LMUY = 1
LEY = 1
LKY = 1
LHY = 1
LVY = 1
LTR = 1
LRES = 1
LXAL = 1
LYKA = 1
LVYKA = 1
LS = 1
LKYC = 1
LKZC = 1
LVMX = 1
LMX = 1
LMY = 1
LMP = 1
$---------------------------------------------longitudinal coefficients
[LONGITUDINAL_COEFFICIENTS]
PCX1 = 1.5
PDX1 = 2.0648
PDX2 = -0.15342
PDX3 = 12.1289
PEX1 = -0.17784
PEX2 = 0.563
PEX3 = 0.6
PEX4 = -0.9
PKX1 = 35.1456
PKX2 = 24.3339
PKX3 = -1.2273
PHX1 = -0.0050602
PHX2 = 0.0013983
PVX1 = 0.034156
PVX2 = -0.1
PPX1 = -1.7796
PPX2 = -2.9784
PPX3 = -0.3492
PPX4 = 0.38929
RBX1 = 11.1331
RBX2 = 12.6807
RBX3 = 4.153
RCX1 = 1.1596
REX1 = -0.068608
REX2 = 0.052784
RHX1 = 0.032035
$----------------------------------------------overturning coefficients
[OVERTURNING_COEFFICIENTS]
QSX1 = -0.0078067
QSX2 = 0.51737
QSX3 = 0.091207
QSX4 = -0.10015
QSX5 = 0
QSX6 = 0
QSX7 = 10.9752
QSX8 = 0
QSX9 = 0
QSX10 = 0
QSX11 = 0
QSX12 = 0
QSX13 = 0
QSX14 = 0
PPMX1 = 1.0792
$--------------------------------------------------lateral coefficients
[LATERAL_COEFFICIENTS]
PCY1 = 1.6675
PDY1 = 2.1768
PDY2 = -0.48644
PDY3 = 5.9709
PEY1 = 0.51141
PEY2 = -0.012972
PEY3 = -0.050898
PEY4 = -4.5261
PEY5 = -70.188
PKY1 = -58.0824
PKY2 = 2.1273
PKY3 = 1
PKY4 = 2
PKY5 = 34.5112
PKY6 = -5.7422
PKY7 = -3.8244
PHY1 = 0.0019699
PHY2 = 0.00075738
PVY1 = 0.051696
PVY2 = 0.014416
PVY3 = 1.0271
PVY4 = 1.2347
PPY1 = 0.82385
PPY2 = 1.4074
PPY3 = -0.17124
PPY4 = -0.11678
PPY5 = -0.60425
RBY1 = 27.8314
RBY2 = 18.3875
RBY3 = 0.030113
RBY4 = -90
RCY1 = 0.98399
REY1 = 0.66603
REY2 = 0.42965
RHY1 = 0.0027678
RHY2 = 0.0092365
RVY1 = 0
RVY2 = 0
RVY3 = 0
RVY4 = 0
RVY5 = 0
RVY6 = 0
$--------------------------------------------------rolling coefficients
[ROLLING_COEFFICIENTS]
QSY1 = 0
QSY2 = 0
QSY3 = 0
QSY4 = 0
QSY5 = 0
QSY6 = 0
QSY7 = 0
QSY8 = 0
$-------------------------------------------------aligning coefficients
[ALIGNING_COEFFICIENTS]
QBZ1 = 0.16633
QBZ2 = 0.0056748
QBZ3 = 0.067521
QBZ4 = 50.8558
QBZ5 = -50.9946
QBZ9 = 15.4176
QBZ10 = 0.42498
QCZ1 = 49.5583
QDZ1 = 0.14773
QDZ2 = -0.059219
QDZ3 = 1.6287
QDZ4 = -35.6742
QDZ6 = 0.01414
QDZ7 = -0.0050238
QDZ8 = -1.3791
QDZ9 = 1.8212
QDZ10 = -0.034803
QDZ11 = -1.5176
QEZ1 = -0.21511
QEZ2 = 1.0726
QEZ3 = 1.5408
QEZ4 = -1.1908
QEZ5 = 16.7096
QHZ1 = 0.0042517
QHZ2 = -0.0030017
QHZ3 = 0.05872
QHZ4 = -0.17859
PPZ1 = -3.0638
PPZ2 = -2.8093
SSZ1 = -0.025289
SSZ2 = -0.051196
SSZ3 = 0
SSZ4 = 0
$-------------------------------------------------turnslip coefficients
[TURNSLIP_COEFFICIENTS]
PDXP1 = 0
PDXP2 = 0
PDXP3 = 0
PKYP1 = 0
PDYP1 = 0
PDYP2 = 0
PDYP3 = 0
PDYP4 = 0
PHYP1 = 0
PHYP2 = 0
PHYP3 = 0
PHYP4 = 0
PECP1 = 0
PECP2 = 0
QDTP1 = 0
QCRP1 = 0
QCRP2 = 0
QBRP1 = 0
QDRP1 = 0