Project 2
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% Skill-Lync -> Hybrid Drives Development
% Project 2 - Conventional Vehicle Model
% By - pavanraj
% Description:
%% Initialization File
% Initialization file for loading the variables for the Conventional vehicle simulation model
%% Vehicle Parameters
M = 1500; % Vehicle weight [kg]
Cr = 0.015; % Rolling resistance coefficient
Rw = 0.35; % Rolling radius [m]
Cd = 0.33; % Drag coefficient
Af = 2; % Frontal area [m2]
rhoa = 1.225; % Air density [kg/m3]
Vwind = 0; % Wind velocity [kmph]
g = 9.81; % Gravitational acceleration [m/s2]
theta = 0; % Slope [deg]
myu_road_tire = 1; % Coefficient of friction between tire - road interface
%% Powertrain Parameters
FDR = 3.73; % Final Drive Ratio / Diff Ratio
DiffEff = 0.99; % Differential efficiency (assumed independent of operating conditions)
NumGear = 4; % Number of gears in transmission
TransGear = 1:NumGear; % Transmission Gears
TransRatio = [3.06, 1.62, 1, 0.7]; % Corresponding Transmission ratios
TransEff = 0.90; % Transmission efficiency (assumed independent of operating conditions)
TC_SR = [0, .5, .6, .7, .8, .87, .92, .94, .96, .97, 1]; % Torque Converter Speed Ratio (Matlab default) -> TurbSpd/PumpSpd
TC_TR = [2.232, 1.5462, 1.4058, 1.2746, 1.1528, 1.0732, 1.0192, 1, 1, 1, 1]; % Torque Converter Torque Ratio (Matlab default) -> TurbTrq/PumpTrq
TC_Cf = 1e-4 * [.72558, .66322, .63463, .59042, .51331, .4144, .29287, .22444, .12186, .04386, 0]; % Torque Converter Capacity Factor (Matlab default) -> PumpTrq/PumpSpd^2 [Nm/rpm^2]
load Vehicle_Data/Engine_data_1_9L_Saturn_95kW
load Vehicle_Data/Trans_Shift_map
load Vehicle_Data/EM_Honda_Insight_10kW
load Vehicle_Data/ESS_NiMH6
load Pack_Current_test.mat
ESS_init_soc = 0.8;
ESS_soc_max = 0.9;
ESS_soc_min = 0.6;
%% Driver Parameters
Kp = 10; % P gain for driver PI model
Ki = 1; % I gain for driver PI model
Drive_Cycle_Select = 1; % Drive cycle selection
% Drive cycle values in m/s
if Drive_Cycle_Select == 1
load DriveCycles/UDDS;
elseif Drive_Cycle_Select == 2
load DriveCycles/HWFET;
elseif Drive_Cycle_Select == 3
load DriveCycles/FTP;
elseif Drive_Cycle_Select == 4
load DriveCycles/US06;
elseif Drive_Cycle_Select == 5
load DriveCycles/NEDC;
elseif Drive_Cycle_Select == 6
load DriveCycles/WLTP;
end
time_end = cyc_time(end);
vel_end = cyc_vel(end);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% Skill-Lync -> Hybrid Drives Development
% Project 2 - p1 Vehicle Model
% By - PAVANRAJ
%% Post-Processing File
% File for post-processing the p1 vehicle simulation model data
%% Renaming variables
Time = out.tout;
Acc_Pdl = out.Acc_Pdl; % Accelerator pedal [%]
Brk_Pdl = out.Brk_Pdl; % Brake pedal [%]
Gear_Cmd = out.Gear_Cmd; % Transmission gear
Trans_In_Spd = out.Trans_In_Spd*30/pi; % Transmission input speed [rpm]
Trans_Out_Spd = out.Trans_Out_Spd*30/pi;% Transmission output speed [rpm]
Eng_Trq = out.Eng_Trq; % Engine brake torque [Nm]
Eng_Spd = out.Eng_Spd*30/pi; % Engine speed [rpm]
Eng_fuel_rate = out.Eng_Fuel_Rate; % Engine fuel rate [gkw]
Eng_Fuel_Consump = out.Eng_Fuel_Consump; % Engine fuel consumption [g]
Veh_Dist = out.Veh_Dist/1000; % Distance traveled by vehicle [km]
Veh_Spd = out.Veh_Spd*18/5; % Vehicle speed [kmph]
Whl_Spd = out.Wheel_Spd*30/pi; % wheel speed [rpm]
%% Data Analysis
Eng_Pwr = Eng_Spd*pi/30.*Eng_Trq/1000; % Engine brake power [kW]
Rhoa_petrol = 0.7197; % Density of petrol [kg/L]
kmpl = Veh_Dist(end)/(Eng_Fuel_Consump(end)/1000/Rhoa_petrol); % Fuel economy [kmpl]
%display Fuel Economy of the cycle
fprintf('Fuel Economy = %2.2fKmplnn',kmpl);
%% Plotting vehicle characteristics time trace
figure(1)
subplot(3,1,1)
hold on
grid on
plot(Time, Acc_Pdl, 'linewidth', 1, 'Color', 'black');
plot(Time, Brk_Pdl, 'linewidth', 1, 'Color', 'red');
ylabel('Pedal [%]')
xlim([0, Time(end)]);
legend('Accel','Brake');
subplot(3,1,2)
hold on
grid on
plot(Time,Gear_Cmd,'LineWidth',1,'Color','black');
xlim([0, Time(end)]);
xlabel('Time [S]')
ylabel('Gear');
subplot(3,1,3)
hold on
grid on
plot(Time,Veh_Spd,'LineWidth',1,'Color','red');
plot(Time,Veh_Dist,'LineWidth',1,'Color','blue');
xlim([0, Time(end)]);
xlabel('Time [S]')
ylabel('speed,distance')
legend('speed','distance [km]');
figure(2)
subplot(4,1,1)
hold on
grid on
plot(Time,Eng_Trq, 'linewidth', 1, 'Color', 'black');
ylabel('Torque [Nm]')
xlim([0, Time(end)]);
subplot(4,1,2)
hold on
grid on
plot(Time,Eng_Pwr, 'linewidth', 1, 'Color', 'red');
ylabel('Power [Kw]')
xlim([0, Time(end)]);
subplot(4,1,3)
hold on
grid on
plot(Time,Eng_fuel_rate, 'linewidth', 1, 'Color', 'red');
xlabel('Time [S]');
ylabel('Engine Fuel Rate');
subplot(4,1,4)
hold on
grid on
plot(Time,Eng_Fuel_Consump,'linewidth',1,'Color','blue');
xlim([0,Time(end)]);
xlabel('Time [sec]')
ylabel('Fuel Consump [g]')
figure(3)
hold on
grid on
plot(cyc_time,cyc_vel,'LineWidth',1,'Color','red');
xlabel('Time [sec]')
ylabel('velocity [Mph]')
title('drive cycle')
%% ploting batery properties
% ploting the all out put required graphs
figure(4)
subplot(5,1,1)
plot(I_pack,'b');
grid on;
xlabel('TIME [S]');
ylabel('Pack Current [A]');
legend('Pack_Current');
ylim([-60 60]);
subplot(5,1,2)
plot(out.SOC,'r');
grid on;
xlabel('Time [S]');
ylabel('SOC');
legend('SOC');
ylim([-5 5]);
subplot(5,1,3)
plot(out.R_graph,'m-');
grid on;
xlabel('Time [S]');
ylabel('Term Volt [v]');
legend('Terminal voltage');
ylim([120 190]);
subplot(5,1,4)
plot(out.Internal_resistance,'g');
grid on;
xlabel('Time [S]');
ylabel('Int res [ohm]');
legend('Cell internal resistance');
ylim([0.015 0.035]);
figure(5)
plot(I_pack,'b');
grid on;
xlabel('TIME [S]');
ylabel('Pack Current [A]');
legend('Pack_Current');
ylim([-60 60])
% Test Results
Avg_Speed = sum(Veh_Spd)/length(Veh_Spd);
Max_Speed = maxk(Veh_Spd,1);
table(Avg_Speed,Max_Speed,kmpl,'VariableNames',{'Average speed [kmph]','maximum speed [kmph]','Fuel economey [kmph]'});
%% Plotting engine operating points on engine fuel map
figure('Name','Engine Operating Points');
hold on
grid on
title('Engine Operating Points','fontsize',12);
p1=plot(Eng_map_spd*30/pi,Eng_max_trq,'linewidth',2,'color','r');
p2=line([Eng_map_spd(end)*30/pi Eng_map_spd(end)*30/pi],[0 Eng_max_trq(end)]);
set(p2,'linewidth',1,'linestyle','--','color','k');
p3=line([Eng_map_spd(1)*30/pi Eng_map_spd(1)*30/pi],[0 Eng_max_trq(1)]);
set(p3,'linewidth',1,'linestyle','--','color','k');
[c,h]=contour(Eng_map_spd'*30/pi,Eng_map_trq,Eng_fuel_map');
clabel(c,h);
xlabel('Engine Speed [rpm]');
ylabel('Engine Torque [Nm]');
title('Engine Fuel Consumption Map [g/s]');
p4 = plot(Eng_Spd,Eng_Trq,'.','color','k');
legend([p1 p4],{'Max Curve','Operating Points'});
axis([Eng_map_spd(1)*30/pi Eng_map_spd(end)*30/pi 0 Eng_map_trq(end)]);