%Series HEV model parameter
%Setting and loading  of Vehicle parameter, control constant and characteristic
clear;

%ECU switch setting
ECU = 0;  %1:Real ECU, 0:ECU model

denormalNumber=10^(-15);

%Loading of Characteristic map
 map_para_s_o;

%Sampling time [s]
stepsize = 0.0005;
%Gravity acceleration [m/s2]
g = 9.80665;

%Vehicle mass[kg]
ms = 12000; 
%Tire radius[m]
rd =0.47;

%Load expression [kgf]   Fr = la + lb*V[km/h] + lc*V[km/h]*V[km/h]
mu = 0.008;  %Rolling resistance coefficient
W = ms * g;  %Vehicle total weight[N]
A = 6.0;  %Frontal projection area [m2]
Cd =0.005*g;  %Air resistance coefficient
la = 96;%(mu*W)/g;
lb = 0.0;
lc = 0.03;%Cd*A/g; %cd=kN*h2/(m2*km2);

%Reduction gear ratio
reduction_gear_ratio = 6.2;

%Differential gear ratio
nf = 6.1;

%Rotaty inertia moment [kgm2]
%j1:Vehicle body and tire, j2:Differential gear, j3:Transmission output side, j4:Transmission input side
j1 = 800;
tire_j=0.001;
j2 = 0.6;
j3 = 0.4;
j4 = 0.0;
%Engine inertia moment
engine_j = 0.13*g;
j5 = engine_j;

%Differential gear efficiency [%]
diff_efficiency = 0.95;

%Reduction gear efficiency [%]
reduce_efficiency =0.95;

gen_gear_eff=0.95;
gen_gear_efficiency=0.95;
gen_gear_ratio=2.23;
Gen2Eg_ratio=gen_gear_ratio;

%Viscosity resistance
cf5 = 0.005;

%Engine revolution limit [rpm]
Rev_limit = 2500.0;
%Dynamo generated torque at engine starting [Nm]
ST_TQ = 700;
%Idle revolution[r/min]
eng_const1=1;
eng_const2=2;

%Stiction at vehicle stop [Nm]
%SEISHI_masatu_tq = la * g * rd * 1.1;

%Capacitor initial state of charge [%]
initial_capa_capacity = 93.0;
%Capacitance [F]
capacity = 35.0;
%Capacitor max voltage [V]
cap_max_vol = 330.0;
%Capacitor max voltage at charge [V]
cap_max_vol_ch = 350.0;
%Capacitor internal resistance [Ħ]
capa_r = 0.095;

%Battery initial state of charge [%]
initial_batt_capacity = 60.0;
%Battery capacity [Ah]
batt_capacity = 6;
%Battery state of charge [%]
SOC = [0 20 40 60 80 100];
%Open circuit voltage [V]
open_voltage = [100 290 300 310 320 330];
%Open circuit voltage (discharge) [V]
disch_open_voltage = [100 290 300 310 320 330];
%Open circuit voltage(charge) [V]
ch_open_voltage = [100 290 300 310 320 330];
%Internal resistance (discharge) [Ħ]
disch_inside_resi = [1 1 1 1 1 1]*0.001;
%Internal resistance (charge) [Ħ]
ch_inside_resi = [1 1 1 1 1 1]*0.001;

%Number of motor
motor_num = 1;
%Motor inertia moment [kgm2] 
motor_j = (0.07) * motor_num;
%Motor normal rotation limit [rpm]
upper_rev = 13000;
%Motor reverse rotation limit [rpm]
lower_rev = -2000.0;

%Generator inertia moment [kgm2]
generator_j = 0.09;

%Time constant
tcm1 = 0.01; %motor
tcm2 = 0.01; %motor
tcm3 = 0.01; %motor
tcm4 = 0.01; %motor
tce1 = 0.005; %engine
tcg1 = 0.01; %generator
tcg2 = 0.01; %generator
tcg3 = 0.01; %generator
tcg4 = 0.01; %generator
tcg5 = 0.01; %generator
tcg6 = 0.005;  %generator

%Engine revolution control constant
adj_EngASR_PID_P = 0.1;
adj_EngASR_PID_I = 0.1;
adj_EngASR_PID_D = 0;
cons_0 = -3000;

%Engine revolution for starter torque addition
starter_rev_limit = 1500;

%Generator revolution control constant
adj_GenASR_PID_P = 1;
adj_GenASR_PID_I = 5;
adj_GenASR_PID_D = 0;