Project 2

QUESTION : DEVELOPMENT OF FORWARD ENERGY BASED FUEL CONSUMPTION MODEL OF A P1 HYBRID VEHICLE

MATLAB SCRIPT FOR INITIATION OF HYBRID VEHICLE :-

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%% Skill-Lync -> Hybrid Drives Development

% Project 1 - Conventional Vehicle Model

% By - Ullekh Gambhira

% 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

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 = 2; % 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);

MATLAB COMMAND WINDOW OF HYBRID VEHICLE :-

>> Project2_InitFile
>> 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 DriveCycles/UDDS
>> load DriveCycles/HWFET
>>

• Drive cycle which one we required for analysis is selected by loading that cycle has provided in question link

SIMULINK MODELLING :-

Initiation file script MATLAB Run :-

Hybrid power train simulink model :-

Driver Model :

Controller Model :-

Engine Control Unit :

Idle Speed Controller :-

Electric Motor Control Unit :-

Transmission Control Block :

Brake Control Unit :

Powertrain Model :-

Battery Model :-

SOC Estimation Model :

Engine Model :-

Electric Motor Model :-

Crank Hub :-

Transmission Model :-

Torque Converter :-

Differential Model :-

Wheel Model :-

Vehicle Dynamics :

SIMULATION RESULTS :-

1. Drive Cycle

UDDS :

HWFET :

3. The Time trace for engine, Transmission input,transmission output and wheel speed

UDDS :

HWFET :

4. Engine troque, Engine power, fuel rate and fuel consumption

UDDS :

HWFET :

5. Accelerator and brake pedal, gear command, vehicle desired and actual speed, vehicle distance

UDDS:

HWFET:

6. The engine operating points on the engine fuel map

7. The Motor speed, torque, motor input power.

UDDS:

HWFET:

8. The Pack voltage, current and Soc

UDDS:

HWFET:

9. The motor operating points on the motor input power map, with the motor max motoring and generating
   torque curves.