Compare commits

...

11 Commits

Author SHA1 Message Date
235b1cd14c PT:V1 2024-01-18 23:20:11 +01:00
f7893173d9 V1 2024-01-10 16:03:41 +01:00
fcea576455 Dateien hochladen nach „tire_models“
40375 coenering mit 43164 drivebrake
unloaded radius = 20.574mm
rim = 7inch
normpress = 100000
alle paramenter konvergieren!
2024-01-09 20:18:33 +01:00
c50f370f66 „tire_models/40375mit43164_20.574_7inch_Fx und Mz0 konvergieren nicht.tir“ löschen 2024-01-09 20:16:36 +01:00
93ded23c72 Dateien hochladen nach „tire_models“ 2024-01-09 15:30:57 +01:00
f443fa5e30 „40375mit43164_20.574_7inch_Fx und Mz0 konvergieren nicht.tir“ löschen 2024-01-09 15:30:48 +01:00
0b3786ccdb Dateien hochladen nach „“ 2024-01-09 15:30:28 +01:00
5deb4cb4e4 Merge branch 'main' of https://git.fasttube.de/t.korjakow/ft24_pacejka_tyre 2024-01-05 19:47:45 +01:00
cf7d0d5632 Added mixed model 2024-01-05 19:47:37 +01:00
644bca2082 Merge branch 'main' of ssh://git.fasttube.de:313/t.korjakow/ft24_pacejka_tyre 2024-01-05 19:47:12 +01:00
824f9b1e74 PT:plot to compare. FX, FY, Mx, MZ 2024-01-05 19:46:37 +01:00
9 changed files with 927 additions and 32 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)
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); % SA-MX
FY = smoothdata(FY, "gaussian",100); % figure()
if i < length(subdirectories)/3 % for i=1:length(subdirectories)
plot(SA, FY, 'LineWidth', 1) % fileMeasurement = fullfile("..\sorted_data", subdirectories{i}, "\cornering.mat");
elseif i <= length(subdirectories)/3*2 % parser = tydex.parsers.FSAETTC_SI_ISO_Mat();
plot(SA, FY, '--','LineWidth', 1) % measurements = parser.run(fileMeasurement);
else % measurement = measurements(n);
plot(SA, FY, '-.','LineWidth', 1) %
% [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");
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");
plot(SA, FY, 'LineWidth', 1)
hold on
plot(SA1, FY1, 'LineWidth', 1)
legend("Cornering", "drivebrake")
catch
disp(fileMeasurement)
end
end end
SAA{i} = SA; title("Similarities, n=",n)
FYY{i} = FY; xlabel("SA")
hold on ylabel("FY")
end end
legend(subdirectoryNames{:}) %% without plot
title("Similarities, n=",n) % mkdir('similarityDataFY');
xlabel("Slip Angle") % for n=1:60
ylabel("FY") % for i=1:length(subdirectories)
%% without plot % 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")
@ -81,4 +83,264 @@ title("Mode of the Absolute Residuials")
xlabel("Order of the tires(8 is our traget tire)") 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