Project 1

Explanation:

The Matlab program has been downloaded from the above file.

The Engine and the vehicle dynamics are Empty, which we need to get It done.

The detailed Engine section has been explained below.

Engine

Here the input is the Engine speed which is given to the 1D look table which consist of the required torque level in order the move the vehicle. The various points is presented in the matlab section.The saturation Dynamics is used to make sure the torque is as per the requirement. And the output is the torque. In order to Know the Fuel consumption we are having the 2d lookup table which is used to the give the fuel consumption rate of which on integrating we will be getting the fuel consumption. The datas are loaded in the Workspace for the further calculation.

Vehicle Dynamics 

The vehicle dynmaics is designed by considering the various forces.

1.Rolling Resistance force

2.Aerdynamics force

3.Acceralation force

4.Grading force( since there is no slope the grading force is avoided)

This various forces are subracted with the Tractive force which is then divided by the Mass.By which we will be getting the Acceralation on which integrating we will be getting the velocity and on further integrating we will be geting the displacement.

Matlab Program:

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

%% 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

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

 

The simulink file has been pasted below.

Plotting:

The matlab program for the plotting section has been pasted below.

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

%% Skill-Lync -> Hybrid Drives Development

% Project 1 - Conventional Vehicle Model

% By - Ullekh Gambhira

% Description:

% The following file is to used with the Initialization file and Simulink

% results for post-processing of the data.

% The file is incomplete, has some examples of how the data is imported

% into a readable format, data analysis and plotting examples.

% Based on the project task, use commands in a similar format as given

% below, to analyse the data and generate required plots.

% You can change variable names if desired

%% Post-Processing File

% File for post-processing the Conventional 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_Consump = out.Eng_Fuel_Consumption; % 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]

%% Data Analysis

Eng_Pwr = (Eng_Spd.*pi/30).*(Eng_Trq/1000); % Engine brake power [kW]

rho_petrol = 0.7197; % Density of petrol [kg/L]

kmpl = Veh_Dist(end)/(Eng_Fuel_Consump(end)/1000/rho_petrol); % Fuel economy [kmpl]

%% Plotting vehicle characteristics time trace

figure()

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, Trans_In_Spd, 'linewidth', 1, 'Color', 'black');

plot(Time, Trans_Out_Spd, 'linewidth', 1, 'Color', 'red');

ylabel('Speed [rpm]')

xlim([0, Time(end)]);

ylim([0, 6000]);

legend('Trans Input','Trans Output');

subplot(3,1,3)

hold on

grid on

plot(Time, Gear_Cmd, 'linewidth', 1, 'Color', 'black');

xlim([0, Time(end)]);

xlabel('Time [s]');

ylabel('Gear');

%% 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)]);

UDDS CYCLE

for this cycle the output has been pasted below.

In this graph we can see that there is exchange of gears various times by which the speed is within a certain limits

HWFET CYCLE

For this cycle we just need to select the @2 section from the drive cycle programs.