79 |
HEV Model Demo Script
80 |
81 |
82 |
83 | Model
84 |
1. Open HEV Model
85 |
2. Run Drive Cycles
86 |
3. Estimate Fuel Economy With Measured Data
87 |
88 |
Simulate
89 |
4. Distribute Simulations Using Parallel Computing
90 |
5. Test and Analyze Test Case Management
91 |
6. Document: View Report, Generate Report
92 |
7. Power Quality Plots: Create Plots
93 |
94 |
Deploy
95 |
8. Deploy to Real-Time Hardware for HIL Testing
96 |
9. Deploy Simulation As a Standalone Executable
97 |
98 |
99 |
100 | Copyright 2011-2024 The MathWorks™, Inc.
101 |
106 |
107 |
138 |
139 |
140 |
--------------------------------------------------------------------------------
/Scripts_Data/HEV_Model_Demo_Script.m:
--------------------------------------------------------------------------------
1 | %% HEV Model Demo Script
2 | %
3 | %
4 | %
5 | %
6 | % Model
7 | %
1. Open HEV Model
8 | %
2. Run Drive Cycles
9 | %
3. Estimate Fuel Economy With Measured Data
10 | %
11 | %
Simulate
12 | %
4. Distribute Simulations Using Parallel Computing
13 | %
5. Test and Analyze Test Case Management
14 | %
6. Document: View Report, Generate Report
15 | %
7. Power Quality Plots: Create Plots
16 | %
17 | %
Deploy
18 | %
8. Deploy to Real-Time Hardware for HIL Testing
19 | %
9. Deploy Simulation As a Standalone Executable
20 | %
21 | %
22 | %
23 | %
24 | % Copyright 2011-2024 The MathWorks(TM), Inc.
25 |
26 |
27 |
--------------------------------------------------------------------------------
/Scripts_Data/HEV_Model_Demo_Script_web.m:
--------------------------------------------------------------------------------
1 | web('HEV_Model_Demo_Script.html')
--------------------------------------------------------------------------------
/Scripts_Data/HEV_Model_PARAM.m:
--------------------------------------------------------------------------------
1 | % Copyright 2011-2024 The MathWorks, Inc.
2 |
3 | %disp('DEFINING PARAMETERS...');
4 | % HEV SERIES PARALLEL PARAMETERS
5 |
6 | % PREDEFINED LI-ION BATTERY PARAMS
7 | HEV_Param.Battery_Det.Nominal_Voltage = 200;
8 | HEV_Param.Battery_Det.Rated_Capacity = 8.1; % Ampere-hours
9 | HEV_Param.Battery_Det.Initial_SOC = 75.75; % Percent
10 | HEV_Param.Battery_Det.Series_Resistance = 0.2/10; % Ohm
11 |
12 |
13 | % FOR ACCEL TEST
14 | %HEV_Param.Battery_Det.Initial_SOC = 97; % Percent
15 |
16 | % GENERIC BATTERY PARAMS
17 | HEV_Param.Battery_Sys.Nominal_Voltage = 217;
18 | HEV_Param.Battery_Sys.Internal_Resistance = 0.24691;% Ohm
19 | HEV_Param.Battery_Sys.Rated_Capacity = 6.9; % Ampere-hours
20 | HEV_Param.Battery_Sys.Initial_Charge = 6.9; % Ampere-hours
21 | HEV_Param.Battery_Sys.Expn_Voltage = 215.0342; % V
22 | HEV_Param.Battery_Sys.Expn_Charge = 2.3438;
23 |
24 | HEV_Param.Battery_Sys.C1.Capacitance = 2500;
25 | HEV_Param.Battery_Sys.C1.Initial_Voltage = 19;
26 | HEV_Param.Battery_Sys.C1.Series_Resistance = 1e-6;
27 | HEV_Param.Battery_Sys.R2 = 0.3;
28 | HEV_Param.Battery_Sys.R1 = 1.8;
29 | HEV_Param.Battery_Sys.Maximum_Capacity = HEV_Param.Battery_Sys.Rated_Capacity; % Ampere-hours
30 |
31 |
32 | % ULTRACAPACITOR PARAMETERS
33 | HEV_Param.UltraCapacitor.Nominal_Capacitance = 1000; % Farad
34 | HEV_Param.UltraCapacitor.Rate_C_V = 0.2; % Farad/Volt
35 | HEV_Param.UltraCapacitor.Series_R = 30/3; % Ohm
36 | HEV_Param.UltraCapacitor.Self_Discharge_R = 500; % Ohm
37 | HEV_Param.UltraCapacitor.Initial_Voltage = 217; % Volt
38 |
39 |
40 | % MOTOR PARAMETERS
41 | HEV_Param.Motor.Stator_Resistance = 0.0065*14; % Ampere-hours
42 | HEV_Param.Motor.Stator_Resistance = 0.0065*14;
43 | HEV_Param.Motor.TorqSpdLUT.SpeedRPM = [0 1200 2000 3000 4000 5000 6000 6500 10000];
44 | HEV_Param.Motor.TorqSpdLUT.TorqueNm = [400 400 225 150 100 80 70 0 0];
45 | HEV_Param.Motor.Damping = 1e-5; %N*m/(rad/s)
46 | HEV_Param.Motor.TorqueControl_TimeConst = 0.02*2/1.5;
47 | HEV_Param.Motor.Shaft_Inertia = 0.2;
48 | HEV_Param.Motor.Series_Resistance = 0.01; %CHG
49 | HEV_Param.Motor.Inductances = [0.001597972349731 0.002057052250467];
50 | HEV_Param.Motor.Efficiency = 91;
51 |
52 |
53 | % GENERATOR PARAMETERS
54 | HEV_Param.Generator.Stator_Resistance = 0.05*0.095;
55 | HEV_Param.Generator.Inductances = [0.000635 0.000635];
56 | HEV_Param.Generator.TorqSpdLUT.SpeedRPM = [ 0 1200 2000 3000 4000 10000 15000];
57 | HEV_Param.Generator.TorqSpdLUT.TorqueNm = [400 400 250 150 110 0 0];
58 | HEV_Param.Generator.Damping = 1e-5; %N*m/(rad/s)
59 | HEV_Param.Generator.TorqueControl_TimeConst = 0.02*2;
60 | HEV_Param.Generator.Shaft_Inertia = 0.2;
61 | HEV_Param.Generator.Series_Resistance = 0.01; %CHG
62 |
63 | % -- trq indep losses adds compliance -- 10W
64 | HEV_Param.Generator.TorqIndep_Elec_Loss = 0; %W
65 |
66 | % CHG -- ADDED COMPLIANCE
67 | HEV_Param.Generator.Shaft_Stiffness = 1e3; %W
68 | HEV_Param.Generator.Shaft_Damping = 1e2; %W
69 |
70 | % DC-DC CONVERTER PARAMETERS
71 | HEV_Param.DCDCConv.Output_Voltage = 500; % Volts
72 | HEV_Param.DCDCConv.Resistance_Losses = 1000/40^2; % Ohm
73 | HEV_Param.DCDCConv.Kp = 0.01;
74 | HEV_Param.DCDCConv.Ki = 10;
75 | HEV_Param.DCDCConv.MinVin = 20;
76 |
77 | HEV_Param.DCDCConv.Mean_Boost.Kp = 0.001;
78 | HEV_Param.DCDCConv.Mean_Boost.Ki = 1;
79 |
80 | HEV_Param.DCDCConv.EPower2Heat = 0.1; % Watts/Watts
81 | HEV_Param.DCDCConv.Thermal_Mass = 0.1*10; % kg
82 | HEV_Param.DCDCConv.Specific_Heat = 100; % J/kg/K
83 | HEV_Param.DCDCConv.Initial_Temperature = 25; % C
84 | HEV_Param.DCDCConv.Air_Temperature = 298; % K
85 | HEV_Param.DCDCConv.Convection.Area = 20; % cm^2
86 | HEV_Param.DCDCConv.Convection.HT_Coefficient = 100; % W/(m^2*K)
87 |
88 |
89 | % CONTROLLER PARAMETERS
90 | HEV_Param.Control.Engine_Start_RPM = 800; % RPM
91 | HEV_Param.Control.Engine_Stop_RPM = 790; % RPM
92 | HEV_Param.Control.Mode_Logic_TS = 0.1;
93 | HEV_Param.Control.ICE.Kp = 0.02;
94 | HEV_Param.Control.ICE.Ki = 0.01;
95 | HEV_Param.Control.Gen.Kp = 10;
96 | HEV_Param.Control.Gen.Ki = 3;
97 | HEV_Param.Control.Mot.Kp = 500;
98 | HEV_Param.Control.Mot.Ki = 300;
99 | HEV_Param.Control.Veh_Spd.Kp = 0.02;
100 | HEV_Param.Control.Veh_Spd.Ki = 0.04;
101 |
102 |
103 | % VEHICLE PARAMETERS
104 | HEV_Param.Vehicle.Mass = 600*2; % kg
105 | HEV_Param.Vehicle.Tire_Radius = 0.3; % m
106 | HEV_Param.Vehicle.Wheel_Inertia = 0.1; % kg*m^2
107 | HEV_Param.Vehicle.Aero_Drag_Coeff = 0.26;
108 | HEV_Param.Vehicle.Incline = 0;
109 | HEV_Param.Vehicle.Engine_Vehicle_Gear_Ratio = 1.3;
110 | HEV_Param.Vehicle.Distance_CG_FrontAxle = 1.35;
111 | HEV_Param.Vehicle.Distance_CG_RearAxle = 1.35;
112 | HEV_Param.Vehicle.Distance_CG_Ground = 0.5;
113 | HEV_Param.Vehicle.Frontal_Area = 2.16;
114 | HEV_Param.Vehicle.Tire.Rated_Vertical_Load = 3000; % N
115 | HEV_Param.Vehicle.Tire.Rated_Peak_Long_Force = 3500; % N
116 | HEV_Param.Vehicle.Tire.Slip_At_Peak_Force = 6; % Percent
117 | %HEV_Param.Vehicle.Tire.Relaxation_Length = 0.2*0.09; % Percent
118 | %HEV_Param.Vehicle.Tire.Relaxation_Length = 0.2;
119 | HEV_Param.Vehicle.Tire.Relaxation_Length = 0.25;
120 | HEV_Param.Vehicle.Trans.Inertia = 0.5;
121 | HEV_Param.Vehicle.Trans.Friction = 0.1*0.001;
122 |
123 |
124 | % POWER SPLIT DEVICE PARAMETERS
125 | HEV_Param.Power_Split.Ratio_Ring2Sun = 2.6;
126 |
127 |
128 | % INTERNAL COMBUSTION ENGINE (ICE) PARAMETERS
129 | HEV_Param.ICE.Shaft_Inertia = 0.25;
130 | HEV_Param.ICE.Max_Power = 57000+57000; % W
131 | HEV_Param.ICE.Speed_At_Max_Power = 5000; % RPM
132 | HEV_Param.ICE.Max_Speed = 6000; % RPM
133 | HEV_Param.ICE.Friction = 0.2079; % N*m*s/rad
134 |
135 | HEV_Param.ICE.sensor_time_constant = 0.005; %0.001
136 |
137 | load FuelConsMap
138 | HEV_Param.ICE.FC.Speed_Vector = FCMap.Speed_Vector;
139 | HEV_Param.ICE.FC.Torque_Vector = FCMap.Torque_Vector;
140 | HEV_Param.ICE.FC.Consumption_Table = FCMap.Consumption_Table;
141 | clear FCMap
142 |
143 |
144 |
145 | % Battery Model Parameters - 50Ah
146 | HEV_Param.Battery_Cell.NominalCap = 50; %(Ah)
147 |
148 | % Battery Block Initial Conditions
149 | % Initial State of Charge (% of full charge)
150 | % Initial Electrolyte Temperature (°C, typically same as ambient temp)
151 | HEV_Param.Battery_Cell.SOC_init = 0.8;
152 | HEV_Param.Battery_Cell.theta_init = 25;
153 |
154 | % Battery Block Thermal Parameters
155 | HEV_Param.Battery_Cell.Ctheta = 400; %(J/°C) Thermal Capacitance
156 | HEV_Param.Battery_Cell.Area = 0.01; % (m^2) Surface area of battery exposed to air
157 | HEV_Param.Battery_Cell.Rtheta = 20; %(W/m^2/K) Convective heat transfer coefficient
158 |
159 | % Battery Block Capacity Parameters
160 | % Charge/discharge cycles at ranges of current/temp
161 | HEV_Param.Battery_Cell.Kc = 1.2; %()
162 | HEV_Param.Battery_Cell.Costar = 1.8e+005; %(As)
163 | HEV_Param.Battery_Cell.Kt_Temps = [25 40 60 75]; % Temperature breakpoints for Kt LUT
164 | HEV_Param.Battery_Cell.Kt = [0.80,1.10,1.20,1.12;]; %() LUT output values
165 | HEV_Param.Battery_Cell.delta = 0.73; %()
166 | HEV_Param.Battery_Cell.Istar = 15; %(A) Nominal Current (=cap/disch_t)
167 | HEV_Param.Battery_Cell.theta_f = -40; %(°C) Electrolyte Freezing Temp
168 |
169 | % Battery Block Parasitic Branch Parameters
170 | % End of charge cycle at ranges of current/temp
171 | HEV_Param.Battery_Cell.Ep = 1.95; %(V) Parasitic emf
172 | HEV_Param.Battery_Cell.Gpo = 2.0e-011; %(s)
173 | HEV_Param.Battery_Cell.Vpo = 0.12; %(V)
174 | HEV_Param.Battery_Cell.Ap = 2.0; %()
175 | HEV_Param.Battery_Cell.Taup = 3; % (s)
176 |
177 | % Battery Block Main Branch Parameters
178 | HEV_Param.Battery_Cell.Emo = 2.18*10; % (V) [max o.c. volts per cell]
179 | HEV_Param.Battery_Cell.Ke = 0.0006; %(V/°C)
180 | HEV_Param.Battery_Cell.Ao = -0.6; % ()
181 | HEV_Param.Battery_Cell.Roo = 0.0042; % (Ohm)
182 | HEV_Param.Battery_Cell.R10 = 0.0010; %(Ohm)
183 | HEV_Param.Battery_Cell.A21 = -10.0; %()
184 | HEV_Param.Battery_Cell.A22 = -8.75; %()
185 | HEV_Param.Battery_Cell.R20 = 0.11; %(Ohm)
186 | HEV_Param.Battery_Cell.Tau1 = 100; %(s)
187 |
188 | % Battery Block Number Of Cells Parameters
189 | HEV_Param.Battery_Cell.ns = 6; %() Number of cells in series in each parallel branch
190 |
191 | % Compute initial extracted charge
192 | HEV_Param.Battery_Cell.Qe_init = (1-HEV_Param.Battery_Cell.SOC_init)*HEV_Param.Battery_Cell.Kc*HEV_Param.Battery_Cell.Costar*interp1([HEV_Param.Battery_Cell.theta_f HEV_Param.Battery_Cell.Kt_Temps],[0 HEV_Param.Battery_Cell.Kt],HEV_Param.Battery_Cell.theta_init,'spline');
193 |
194 |
--------------------------------------------------------------------------------
/Scripts_Data/Select_HEV_Model_Systems.m:
--------------------------------------------------------------------------------
1 | function Select_HEV_Model_Systems(modelname,HEV_Config_Name,HEV_Configs)
2 | % Copyright 2011-2024 The MathWorks, Inc.
3 |
4 | % Find block with mask parameters to set variants
5 | f = Simulink.FindOptions('SearchDepth',1,'regexp',true);
6 | elec_sys = getfullname(Simulink.findBlocks(modelname,'popup_electricalvariant','.*',f));
7 | vehi_sys = [modelname '/Vehicle Dynamics'];
8 |
9 | % Extract configuration
10 | HEV_config_index = find(strcmp({HEV_Configs.Type},HEV_Config_Name));
11 |
12 | elec_config = HEV_Configs(HEV_config_index).Electrical.choice;
13 | batt_config = HEV_Configs(HEV_config_index).Battery.choice;
14 | vehi_config = HEV_Configs(HEV_config_index).Vehicle.choice;
15 |
16 | if(~strcmp(elec_config,'Empty'))
17 | set_param(elec_sys,'popup_electricalvariant',elec_config);
18 | set_param(elec_sys,'popup_batteryvariantsystem',batt_config);
19 | if(~strcmpi(batt_config,'cells'))
20 | set_param(elec_sys,'popup_batteryvariantdetailed',batt_config);
21 | end
22 | end
23 |
24 | if(~strcmp(vehi_config,'Empty'))
25 | set_param(vehi_sys,'LabelModeActiveChoice',vehi_config);
26 | end
27 |
28 |
29 |
30 |
31 |
32 |
33 |
--------------------------------------------------------------------------------
/Scripts_Data/Setup_HEV_Model_Configurations.m:
--------------------------------------------------------------------------------
1 | % Copyright 2011-2024 The MathWorks, Inc.
2 |
3 | expModel = bdroot;
4 |
5 | electrical_path = char(find_system(expModel,'MatchFilter',@Simulink.match.allVariants,'Name','Electrical'));
6 | vehicle_path = char(find_system(expModel,'MatchFilter',@Simulink.match.allVariants,'Name','Vehicle Dynamics'));
7 |
8 | HEV_Electrical_Variants = {'Mean Value' 'System Level' 'Detailed'};
9 |
10 |
11 | HEV_Configs_TBL_COLS = {'Type' 'ELECTRICAL' 'BATTERY' 'VEHICLE' 'SIM TIME'};
12 | HEV_Configs_TABLE{1} = {'Sys BD VS' 'System Level' 'Predefined' 'Simple' 10}; % 15
13 | HEV_Configs_TABLE{2} = {'Sys BC VS' 'System Level' 'Generic' 'Simple' 10}; % 15
14 | HEV_Configs_TABLE{3} = {'Sys BD VF' 'System Level' 'Predefined' 'Full' 10}; % 15
15 | HEV_Configs_TABLE{4} = {'Sys BC VF' 'System Level' 'Generic' 'Full' 10}; % 15
16 | HEV_Configs_TABLE{5} = {'Mvl BD VS' 'Mean Value' 'Predefined' 'Simple' 10}; % 15
17 | HEV_Configs_TABLE{6} = {'Mvl BC VS' 'Mean Value' 'Generic' 'Simple' 10}; % 15
18 | HEV_Configs_TABLE{7} = {'Mvl BD VF' 'Mean Value' 'Predefined' 'Full' 10}; % 15
19 | HEV_Configs_TABLE{8} = {'Mvl BC VF' 'Mean Value' 'Generic' 'Full' 10}; % 15
20 | HEV_Configs_TABLE{9} = {'Det BD VS' 'Detailed' 'Predefined' 'Simple' 10}; % 15
21 | HEV_Configs_TABLE{10} = {'Det BC VS' 'Detailed' 'Generic' 'Simple' 10}; % 15
22 | HEV_Configs_TABLE{11} = {'Det BD VF' 'Detailed' 'Predefined' 'Full' 10}; % 15
23 | HEV_Configs_TABLE{12} = {'Det BC VF' 'Detailed' 'Generic' 'Full' 10}; % 15
24 |
25 |
26 | for i=1:length(HEV_Configs_TABLE)
27 | HEV_Configs(i).Type = HEV_Configs_TABLE{i}{1};
28 | HEV_Configs(i).Electrical.choice = HEV_Configs_TABLE{i}{2};
29 | HEV_Configs(i).Battery.choice = HEV_Configs_TABLE{i}{3};
30 | HEV_Configs(i).Vehicle.choice = HEV_Configs_TABLE{i}{4};
31 | HEV_Configs(i).Sim_Time = HEV_Configs_TABLE{i}{5};
32 | end
33 |
34 | for i=1:length(HEV_Configs)
35 | HEV_Configs(i).Electrical.path = {electrical_path};
36 | HEV_Configs(i).Battery.path = {[electrical_path '/' char(HEV_Configs(i).Electrical.choice) '/Battery']};
37 | HEV_Configs(i).Vehicle.path = {vehicle_path};
38 | end
39 |
--------------------------------------------------------------------------------
/Scripts_Data/Test_HEV_Model_SHORT.m:
--------------------------------------------------------------------------------
1 | % Copyright 2011-2024 The MathWorks, Inc.
2 |
3 | expModel = 'HEV_SeriesParallel';
4 | open_system(expModel);
5 |
6 | ModelVariants = {'System Level' 'Mean Value' 'Detailed'};
7 | BattVariants = {'Predefined' 'Generic' 'Cells'};
8 | VehVariants = {'Simple' 'Full'};
9 |
10 | SimDuration = [max(UrbanCycle1.time) max(UrbanCycle2.time) max(UrbanCycle3.time)];
11 |
12 | MV_testInd = 1;
13 | Batt_testInd = 1;
14 | Veh_testInd = 1;
15 |
16 | set_param([expModel '/Vehicle Dynamics'],'OverrideUsingVariant',VehVariants{Veh_testInd});
17 | set_param([expModel '/Electrical'],'popup_electricalvariant',ModelVariants{MV_testInd});
18 | set_param([expModel '/Electrical'],'popup_batteryvariantsystem',BattVariants{Batt_testInd});
19 |
20 | set_param(bdroot,'FastRestart','on')
21 | for MV_ind = MV_testInd:MV_testInd
22 |
23 | MV_str = char(ModelVariants(MV_ind));
24 | MV_name = strrep(MV_str,' ','_');
25 |
26 | for DC_ind=1:3
27 |
28 | for Veh_ind = 1:1 %length(VehVariants)
29 | Veh_str = char(VehVariants(Veh_ind));
30 | Veh_name = strrep(Veh_str,' ','_');
31 | for Batt_ind=2:2 %length(BattVariants)
32 | Batt_str = char(BattVariants(Batt_ind));
33 | Batt_name = strrep(Batt_str,' ','_');
34 |
35 | %if (strcmp(get_param(bdroot,'FastRestart'),'off'))
36 | set_param(expModel,'StopTime',num2str(SimDuration(DC_ind)));
37 | %end
38 | Drive_Cycle_Num = DC_ind;
39 |
40 | disp(['Simulating UC' num2str(DC_ind) ', ' get_param([expModel '/Electrical'],'ActiveVariant') ', '...
41 | get_param([expModel '/Vehicle Dynamics'],'ActiveVariant') ' Vehicle, ',...
42 | get_param([expModel '/Electrical/' MV_str '/Battery'],'ActiveVariant') ' Battery'] );
43 | sim(expModel);
44 |
45 | if exist('Electricals','var')
46 | eval([MV_name '.Electrical = Electricals;']);
47 | eval([MV_name '.Car = Car;']);
48 | eval([MV_name '.Generator = Generator;']);
49 | eval([MV_name '.Motor = Motor;']);
50 | eval([MV_name '.Control_Logic = Control_Logic;']);
51 | eval([MV_name '.DCDC_Conv = DCDC_Temp;']);
52 | end
53 |
54 | %SaveFolder = [SaveFolderRoot '\UC' num2str(DC_ind) '\' MV_name];
55 | %SaveFileName = [SaveFolder '\' MV_name '_DATA_UC' num2str(DC_ind) '_Veh' num2str(Veh_ind) '_Batt' num2str(Batt_ind)];
56 | %disp(['save ' SaveFileName ' ' MV_name]);
57 | %eval(['save ' SaveFileName ' ' MV_name]);
58 | eval(['clear ' MV_name ' Electricals Car Generator Motor Control_Logic']);
59 | end
60 | end
61 | end
62 | end
63 |
64 | set_param(bdroot,'FastRestart','off');
65 |
66 | open('HEV_Model_Report_SHORT.html');
67 |
--------------------------------------------------------------------------------
/Scripts_Data/UrbanCycle1.mat:
--------------------------------------------------------------------------------
https://raw.githubusercontent.com/mathworks/Simscape-HEV-Series-Parallel/8aa2461df9f8f25aa0939f2f691a06cbf61418d7/Scripts_Data/UrbanCycle1.mat
--------------------------------------------------------------------------------
/Scripts_Data/UrbanCycle2.mat:
--------------------------------------------------------------------------------
https://raw.githubusercontent.com/mathworks/Simscape-HEV-Series-Parallel/8aa2461df9f8f25aa0939f2f691a06cbf61418d7/Scripts_Data/UrbanCycle2.mat
--------------------------------------------------------------------------------
/Scripts_Data/UrbanCycle3.mat:
--------------------------------------------------------------------------------
https://raw.githubusercontent.com/mathworks/Simscape-HEV-Series-Parallel/8aa2461df9f8f25aa0939f2f691a06cbf61418d7/Scripts_Data/UrbanCycle3.mat
--------------------------------------------------------------------------------
/Scripts_Data/UrbanCycle4.mat:
--------------------------------------------------------------------------------
https://raw.githubusercontent.com/mathworks/Simscape-HEV-Series-Parallel/8aa2461df9f8f25aa0939f2f691a06cbf61418d7/Scripts_Data/UrbanCycle4.mat
--------------------------------------------------------------------------------
/Scripts_Data/UrbanCycle5.mat:
--------------------------------------------------------------------------------
https://raw.githubusercontent.com/mathworks/Simscape-HEV-Series-Parallel/8aa2461df9f8f25aa0939f2f691a06cbf61418d7/Scripts_Data/UrbanCycle5.mat
--------------------------------------------------------------------------------
/Scripts_Data/ftpdata.mat:
--------------------------------------------------------------------------------
https://raw.githubusercontent.com/mathworks/Simscape-HEV-Series-Parallel/8aa2461df9f8f25aa0939f2f691a06cbf61418d7/Scripts_Data/ftpdata.mat
--------------------------------------------------------------------------------
/Scripts_Data/plot_HEV_Data.m:
--------------------------------------------------------------------------------
1 | function plot_HEV_Data(FigID,SubPlID,DataStr,SigID,ValID,titlestr,PlotChoice)
2 | % Copyright 2011-2024 The MathWorks, Inc.
3 |
4 | figure(FigID)
5 | temp_colororder = get(gca,'defaultAxesColorOrder');
6 |
7 | SP_color = temp_colororder(1,:);
8 | SP_style = '-';
9 |
10 | SY_color = temp_colororder(1,:);
11 | SY_style = '-';
12 |
13 | MV_color = temp_colororder(2,:);
14 | MV_style = '-';
15 |
16 | DT_color = temp_colororder(3,:);
17 | DT_style = '-';
18 | PlotLineWidth = 0.5;
19 |
20 |
21 | subplot(SubPlID(1),SubPlID(2),SubPlID(3));
22 | %subplot(3,1,1)
23 | for pc_ind = 1:length(PlotChoice)
24 | if PlotChoice(pc_ind) == 'D'
25 | legendstrings{pc_ind} = 'Detailed';
26 | DT_Time = evalin('base', ['Detailed.' DataStr '.time']);
27 | DT_Data = evalin('base', ['Detailed.' DataStr '.signals(' num2str(SigID) ').values(:,' num2str(ValID) ');']);
28 | DT_maxind = length(DT_Time);
29 | plot(DT_Time(1:DT_maxind),DT_Data(1:DT_maxind),...
30 | 'LineStyle',DT_style,'Color',DT_color,'LineWidth',4);%PlotLineWidth);
31 | elseif PlotChoice(pc_ind) == 'S'
32 | legendstrings{pc_ind} = 'System Level';
33 | SY_Time = evalin('base', ['System_Level.' DataStr '.time']);
34 | SY_Data = evalin('base', ['System_Level.' DataStr '.signals(' num2str(SigID) ').values(:,' num2str(ValID) ');']);
35 | SY_maxind = length(SY_Time);
36 | plot(SY_Time(1:SY_maxind),SY_Data(1:SY_maxind),...
37 | 'LineStyle',SY_style,'Color',SY_color,'LineWidth',3)%PlotLineWidth);
38 | elseif PlotChoice(pc_ind) == 'M'
39 | legendstrings{pc_ind} = 'Mean Value';
40 | MV_Time = evalin('base', ['Mean_Value.' DataStr '.time']);
41 | MV_Data = evalin('base', ['Mean_Value.' DataStr '.signals(' num2str(SigID) ').values(:,' num2str(ValID) ');']);
42 | MV_maxind = length(MV_Time);
43 | plot(MV_Time(1:MV_maxind),MV_Data(1:MV_maxind),...
44 | 'LineStyle',MV_style,'Color',MV_color,'LineWidth',1);%,PlotLineWidth);
45 | else
46 | legendstrings{pc_ind} = 'Simplified';
47 | SP_Time = evalin('base', ['Simplified.' DataStr '.time']);
48 | SP_Data = evalin('base', ['Simplified.' DataStr '.signals(' num2str(SigID) ').values(:,' num2str(ValID) ');']);
49 | SP_maxind = length(SP_Time);
50 | plot(SP_Time(1:SP_maxind),SP_Data(1:SP_maxind),...
51 | 'LineStyle',SP_style,'Color',SP_color,SP_style,'LineWidth',PlotLineWidth);
52 | end
53 | hold on
54 | end
55 |
56 | hold off
57 | title(titlestr);
58 | legend(legendstrings);
59 | grid on
60 |
61 | end
62 |
--------------------------------------------------------------------------------
/Scripts_Data/shutdown_HEV_Model.m:
--------------------------------------------------------------------------------
1 | % Shutdown script for HEV_SeriesParallel
2 |
3 | % If parameter sweep plot still open, close it
4 | try close(h4_hev_seriesparallel_pct),end %#ok
5 | clear h4_hev_seriesparallel_pct
--------------------------------------------------------------------------------
/Scripts_Data/spsdrivelogodata.mat:
--------------------------------------------------------------------------------
https://raw.githubusercontent.com/mathworks/Simscape-HEV-Series-Parallel/8aa2461df9f8f25aa0939f2f691a06cbf61418d7/Scripts_Data/spsdrivelogodata.mat
--------------------------------------------------------------------------------
/Scripts_Data/startup_HEV_Model.m:
--------------------------------------------------------------------------------
1 | % Copyright 2011-2024 The MathWorks, Inc.
2 |
3 | HEV_Model_HomeDir = pwd;
4 |
5 | HEV_Model_PARAM
6 |
7 | open_start_content = 1;
8 |
9 | % If running in a parallel pool
10 | % do not open model or demo script
11 | if(~isempty(ver('parallel')))
12 | if(~isempty(getCurrentTask()))
13 | open_start_content = 0;
14 | end
15 | end
16 |
17 | if(open_start_content)
18 | web('HEV_Model_Demo_Script.html')
19 | HEV_SeriesParallel
20 | end
21 |
22 |
23 |
24 |
25 |
--------------------------------------------------------------------------------
/Workflows/Fuel_Consumption/Create_Fuel_Consumption_LUT.m:
--------------------------------------------------------------------------------
1 | %% ESTIMATE FUEL CONSUMPTION WITH MEASURED DATA
2 | % Copyright 2011-2024 The MathWorks, Inc.
3 | % CONSTANTS
4 | gasoline_density_gpm3 = 750*1000; % Grams per cubic meter
5 | gallons_per_m3 = 264.172; % gallons per cubic meter
6 |
7 | %% IMPORT EXPERIMENTAL DATA
8 | % http://www.transportation.anl.gov/pdfs/HV/2.pdf
9 | FC_Measured_Data = xlsread('Fuel_Consumption_DATA.xlsx');
10 | FCData.Speed_RPM = FC_Measured_Data(1:(end-1),1);
11 | FCData.Trq_Nm = FC_Measured_Data(1:(end-1),2);
12 | FCData.FuelUse_gps = FC_Measured_Data(1:(end-1),3)/gallons_per_m3*gasoline_density_gpm3;
13 | Fuel_Consumption_PlotData
14 |
15 | %% OPEN SURFACE FITTING TOOL
16 | cftool(FCData.Speed_RPM, FCData.Trq_Nm, FCData.FuelUse_gps)
17 |
18 | % Polynomial, x=2, y=2
19 | % Fit/Save to Workspace
20 | % sftool('FC_SFit_Session.sfit') % Saved Session
21 |
22 | %% GENERATE LOOKUP TABLE
23 | load FCMap_Poly2
24 | numpts = 10;
25 | FCMap.Speed_Vector = linspace(min(FCData.Speed_RPM),max(FCData.Speed_RPM),numpts);
26 | FCMap.Torque_Vector = linspace(min(FCData.Trq_Nm),max(FCData.Trq_Nm),numpts);
27 | [SpeedArray,TorqueArray] = meshgrid(FCMap.Speed_Vector,FCMap.Torque_Vector);
28 | FCMap.Consumption_Table = FCMap_Poly2(SpeedArray,TorqueArray);
29 |
30 | % PLOT LOOKUP TABLE AS SURFACE
31 | surf(SpeedArray,TorqueArray,FCMap.Consumption_Table);
32 |
33 |
34 |
--------------------------------------------------------------------------------
/Workflows/Fuel_Consumption/FCMap_Poly2.mat:
--------------------------------------------------------------------------------
https://raw.githubusercontent.com/mathworks/Simscape-HEV-Series-Parallel/8aa2461df9f8f25aa0939f2f691a06cbf61418d7/Workflows/Fuel_Consumption/FCMap_Poly2.mat
--------------------------------------------------------------------------------
/Workflows/Fuel_Consumption/FC_SFit_Session.sfit:
--------------------------------------------------------------------------------
https://raw.githubusercontent.com/mathworks/Simscape-HEV-Series-Parallel/8aa2461df9f8f25aa0939f2f691a06cbf61418d7/Workflows/Fuel_Consumption/FC_SFit_Session.sfit
--------------------------------------------------------------------------------
/Workflows/Fuel_Consumption/FuelConsMap.mat:
--------------------------------------------------------------------------------
https://raw.githubusercontent.com/mathworks/Simscape-HEV-Series-Parallel/8aa2461df9f8f25aa0939f2f691a06cbf61418d7/Workflows/Fuel_Consumption/FuelConsMap.mat
--------------------------------------------------------------------------------
/Workflows/Fuel_Consumption/Fuel_Consumption_DATA.xlsx:
--------------------------------------------------------------------------------
https://raw.githubusercontent.com/mathworks/Simscape-HEV-Series-Parallel/8aa2461df9f8f25aa0939f2f691a06cbf61418d7/Workflows/Fuel_Consumption/Fuel_Consumption_DATA.xlsx
--------------------------------------------------------------------------------
/Workflows/Fuel_Consumption/Fuel_Consumption_PlotData.m:
--------------------------------------------------------------------------------
1 | % Copyright 2011-2024 The MathWorks, Inc.
2 |
3 | hndl = plot3(FCData.Speed_RPM, FCData.Trq_Nm, FCData.FuelUse_gps,'o','MarkerFaceColor','blue','MarkerEdgeColor','white');
4 | grid on
5 | box on
6 | xlabel('Engine Speed RPM','FontSize',12,'FontWeight','Bold')
7 | ylabel('Engine Torque (Nm)','FontSize',12,'FontWeight','Bold')
8 | zlabel('Fuel Rate (gram/sec)','FontSize',12,'FontWeight','Bold')
9 | title('Fuel Consumption, Measured Data','FontSize',14,'FontWeight','Bold')
10 | set(gca,'XLim',[1200 4100]);
11 | set(gca,'YLim',[30 99]);
12 | set(gca,'ZLim',[0.25 2.75]);
13 | hold on
14 |
--------------------------------------------------------------------------------
/Workflows/Param_Sweep/HEV_Model_PCT_Sim.m:
--------------------------------------------------------------------------------
1 | %% Script to run HEV_SeriesParallel in parallel
2 | % Copyright 2013-2024 The MathWorks(TM), Inc.
3 |
4 | % Move to folder where script is saved
5 | cd(fileparts(which(mfilename)));
6 |
7 | % Open model and save under another name for test
8 | orig_mdl = 'HEV_SeriesParallel';
9 | open_system(orig_mdl);
10 | mdl = [orig_mdl '_pct_temp' ];
11 | save_system(orig_mdl,mdl);
12 |
13 | %% Configure model for tests
14 |
15 | % Set up parameters
16 | HEV_Vehicle_Mass = HEV_Param.Vehicle.Mass;
17 | HEV_Model_Driver_Ki = 0.04;
18 |
19 | % Model settings
20 | Select_HEV_Model_Systems(mdl,'Sys BC VS',HEV_Configs);
21 | HEVSP_tictoc('off');
22 | set_param([mdl '/Vehicle Dynamics/Simple'],'mass','HEV_Vehicle_Mass');
23 | set_param([mdl '/SLRT Scope'],'Commented','on');
24 | set_param(bdroot,'SimscapeLogType','none');
25 | Drive_Cycle_Num = 1;
26 | set_param(mdl,'StopTime',num2str(eval(['UrbanCycle' num2str(Drive_Cycle_Num) '.time(end)'])));
27 |
28 | save_system(mdl);
29 |
30 | %% Generate parameter sets
31 | Mass_array = [1000:40:1600];
32 |
33 | for i=1:length(Mass_array)
34 | simInput(i) = Simulink.SimulationInput(mdl);
35 | simInput(i) = simInput(i).setVariable('HEV_Vehicle_Mass',Mass_array(i));
36 | end
37 |
38 | %% Run one simulation to see time used
39 | timerVal = tic;
40 | sim(mdl)
41 | Elapsed_Sim_Time_single = toc(timerVal);
42 | disp(['Elapsed Simulation Time Single Run: ' num2str(Elapsed_Sim_Time_single)]);
43 |
44 | %% Run parameter sweep in parallel
45 | timerVal = tic;
46 | simOut = parsim(simInput,'ShowSimulationManager','on',...
47 | 'ShowProgress','on','UseFastRestart','on',...
48 | 'TransferBaseWorkspaceVariables','on');
49 | Elapsed_Time_Time_parallel = toc(timerVal);
50 |
51 | %% Calculate elapsed time less setup of parallel
52 | Elapsed_Time_Sweep = ...
53 | (datenum(simOut(end).SimulationMetadata.TimingInfo.WallClockTimestampStop) - ...
54 | datenum(simOut(1).SimulationMetadata.TimingInfo.WallClockTimestampStart)) * 86400;
55 | disp(['Elapsed Sweep Time Total: ' sprintf('%5.2f',Elapsed_Time_Sweep)]);
56 | disp(['Elapsed Sweep Time/(Num Tests): ' sprintf('%5.2f',Elapsed_Time_Sweep/length(simOut))]);
57 |
58 | %% Plot results
59 | plot_sim_res(simOut,'Parallel Test',Elapsed_Time_Time_parallel)
60 |
61 |
62 | %% Close parallel pool
63 | delete(gcp);
64 |
65 | %% Cleanup directory
66 |
67 | bdclose(mdl);
68 | delete([mdl '.slx']);
69 |
70 |
71 |
72 | %% Plot Function
73 | function plot_sim_res(simOut,annotation_str,elapsed_time)
74 |
75 | % Plot Results
76 | fig_handle_name = 'h4_hev_seriesparallel_pct';
77 |
78 | handle_var = evalin('base',['who(''' fig_handle_name ''')']);
79 | if(isempty(handle_var))
80 | evalin('base',[fig_handle_name ' = figure(''Name'', ''' fig_handle_name ''');']);
81 | elseif ~isgraphics(evalin('base',handle_var{:}))
82 | evalin('base',[fig_handle_name ' = figure(''Name'', ''' fig_handle_name ''');']);
83 | end
84 | figure(evalin('base',fig_handle_name))
85 | clf(evalin('base',fig_handle_name))
86 |
87 | for i=length(simOut):-1:1
88 | data = simOut(i).Motor;
89 | plot(data.time(:,1),data.signals(3).values(:,1),'LineWidth',2)
90 | hold all
91 | end
92 | title('Motor Torque','FontSize',16,'FontWeight','Bold');
93 | xlabel('Time (s)','FontSize',12,'FontWeight','Bold');
94 | ylabel('Motor Torque','FontSize',12,'FontWeight','Bold');
95 | Mass_array = evalin('base','Mass_array');
96 | legend(cellstr(num2str(fliplr(Mass_array(1:1:end))')),'FontSize',10);
97 |
98 | text(0.05,0.05,[annotation_str ', Elapsed Time: ' num2str(elapsed_time)],'Color',[1 1 1]*0.6,'Units','Normalized');
99 | end
100 |
101 |
--------------------------------------------------------------------------------
/Workflows/Power_Quality/DC_Power_Quality_Analysis.m:
--------------------------------------------------------------------------------
1 | %% CONSTANTS
2 | % Copyright 2011-2024 The MathWorks, Inc.
3 |
4 | FigurePos1 = [450 499 574 186];
5 | FigurePos2 = [449 69 574 372];
6 | FigurePos3 = [449 475 574 217];
7 | FigurePos4 = [5 476 437 215];
8 | PurpleColor = [102 0 204]/256;
9 | GenColor = [200 102 100]/256;
10 | SampleTime = (2e-6)*50;
11 |
12 | %% LOAD DATA TEST A
13 | load HEV_PQ_Test_A
14 |
15 | %% PLOT DC VOLTAGE
16 | figure(1);
17 | plot(Time,DC_V,'Color',[102 0 204]/256);
18 | title('DC Bus Voltage','FontSize',16,'FontWeight','Bold');
19 | xlabel('Time','FontSize',12,'FontWeight','bold')
20 | ylabel('Voltage (V)','FontSize',12,'FontWeight','Bold')
21 | set(gcf,'Position',FigurePos1);
22 | set(gca,'YLim',[400 600]);
23 |
24 | %% PLOT SPECTROGRAM
25 | figure(2);
26 | spectrogram(DC_V,1024,512,1024,1/(Time(2)-Time(1)),'yaxis');
27 | %spectrogram(DC_V,256,240,256,1/(SampleTime),'yaxis');
28 | title('Spectrogram, DC Bus','FontSize',16,'FontWeight','Bold')
29 | xlabel('Time','FontSize',12,'FontWeight','bold')
30 | ylabel('Frequency (Hz)','FontSize',12,'FontWeight','Bold')
31 | set(gcf,'Position',FigurePos2);
32 |
33 | %% TEST B
34 | close(1);close(2);
35 | load HEV_PQ_Test_B
36 |
37 | %% PLOT SPECTROGRAM
38 | figure(3);
39 | spectrogram(DC_V,1024,512,1024,1/(Time(2)-Time(1)),'yaxis');
40 | title('Spectrogram, DC Bus','FontSize',16,'FontWeight','Bold')
41 | xlabel('Time','FontSize',12,'FontWeight','bold')
42 | ylabel('Frequency (Hz)','FontSize',12,'FontWeight','Bold')
43 | set(gcf,'Position',FigurePos2);
44 | %set(gca,'YLim',[0 1200]);
45 | %view(0,90)
46 |
47 |
48 | %% PLOT GENERATOR SPEED
49 | figure(4);
50 | plot(Time,Gen_Spd,'LineWidth',3,'Color',GenColor);
51 | title('Generator Speed (RPM)','FontSize',16,'FontWeight','Bold');
52 | xlabel('Time','FontSize',12,'FontWeight','bold')
53 | ylabel('Speed (RPM)','FontSize',12,'FontWeight','Bold')
54 | set(gcf,'Position',FigurePos3);
55 |
56 | %% PLOT MOTOR SPEED
57 | figure(5);
58 | plot(Time,Mot_Spd,'LineWidth',3);
59 | title('Motor Speed (RPM)','FontSize',16,'FontWeight','Bold');
60 | xlabel('Time','FontSize',12,'FontWeight','bold')
61 | ylabel('Speed (RPM)','FontSize',12,'FontWeight','Bold')
62 | set(gcf,'Position',FigurePos4);
63 | set(gca,'YLim',[0 450]);
64 |
65 | %% CLOSE ALL
66 | close all
67 |
--------------------------------------------------------------------------------
/Workflows/Power_Quality/HEV_PQ_Test_A.mat:
--------------------------------------------------------------------------------
https://raw.githubusercontent.com/mathworks/Simscape-HEV-Series-Parallel/8aa2461df9f8f25aa0939f2f691a06cbf61418d7/Workflows/Power_Quality/HEV_PQ_Test_A.mat
--------------------------------------------------------------------------------
/Workflows/Power_Quality/HEV_PQ_Test_B.mat:
--------------------------------------------------------------------------------
https://raw.githubusercontent.com/mathworks/Simscape-HEV-Series-Parallel/8aa2461df9f8f25aa0939f2f691a06cbf61418d7/Workflows/Power_Quality/HEV_PQ_Test_B.mat
--------------------------------------------------------------------------------
/Workflows/RSIM/HEV_Model_RSIM.m:
--------------------------------------------------------------------------------
1 | %% OPEN AND CONFIGURE MODEL
2 | % Copyright 2011-2024 The MathWorks, Inc.
3 |
4 | mdl = 'HEV_SeriesParallel';
5 | open_system(mdl);
6 | Select_HEV_Model_Systems(mdl,'Sys BC VF',HEV_Configs);
7 | HEV_Param.Control.Mode_Logic_TS = 0.1;
8 |
9 | %% BUILD STANDALONE EXECUTABLE
10 | disp(['Building standalone executable for ' mdl]);
11 | set_param(bdroot,'SystemTargetFile','rsim.tlc')
12 | set_param(bdroot,'SimscapeLogType','none');
13 | set_param(bdroot,'DefaultParameterBehavior','Inlined');
14 | set_param([bdroot '/SLRT Scope'],'Commented','on');
15 | slbuild(mdl);
16 | disp(['Finished building standalone executable for ' mdl '.mdl']);
17 |
18 | %% CREATE RSIM PARAMETER FILES
19 | disp(' ');disp('Creating RSIM parameter files...')
20 |
21 | Drive_Cycle_Num = 1;
22 | [rtpstruct]= rsimgetrtp(mdl);
23 | save HEV_Model_Parameters_01 rtpstruct;
24 |
25 | Drive_Cycle_Num = 2;
26 | [rtpstruct]= rsimgetrtp(mdl);
27 | save HEV_Model_Parameters_02 rtpstruct;
28 |
29 | Drive_Cycle_Num = 3;
30 | [rtpstruct]= rsimgetrtp(mdl);
31 | save HEV_Model_Parameters_03 rtpstruct;
32 |
33 | disp('Finished creating parameter data files.')
34 |
35 | %% RUN STANDALONE EXECUTABLE
36 | disp(' ');disp('Running RSIM simulations...')
37 | !HEV_SeriesParallel -p HEV_Model_Parameters_01.mat -o results_01.mat -tf 195
38 | !HEV_SeriesParallel -p HEV_Model_Parameters_02.mat -o results_02.mat -tf 400
39 | !HEV_SeriesParallel -p HEV_Model_Parameters_03.mat -o results_03.mat -tf 400
40 | disp('Finished RSIM simulations...')
41 |
42 | %% PLOT RESULTS
43 | colordef black
44 | figure
45 | set(gcf,'Position',[551 135 449 548]);
46 | YAxLims= [0 150];
47 |
48 | subplot(311)
49 | load results_01.mat; HS_h(1) =plot(rt_Car.time,rt_Car.signals(4).values(:,1),'y'); grid on; set(gca,'YLim',YAxLims);
50 | subplot(312)
51 | load results_02.mat; HS_h(2) =plot(rt_Car.time,rt_Car.signals(4).values(:,1),'m'); grid on; set(gca,'YLim',YAxLims);
52 | subplot(313)
53 | load results_03.mat; HS_h(3) =plot(rt_Car.time,rt_Car.signals(4).values(:,1),'r'); grid on; set(gca,'YLim',YAxLims);
54 | hold off;
55 | set(HS_h,'LineWidth',3);
56 |
57 | subplot(311)
58 | title_h = title('Vehicle Speed (kph)');
59 | subplot(313)
60 | xlabel_h = xlabel('Time (s)');
61 | set(title_h,'FontSize',14);
62 | set(xlabel_h,'FontSize',12);
63 | grid on
64 |
65 | %axis([0 200 0 120]);
66 | colordef white
67 |
68 | %% CLEANUP DIRECTORY
69 | delete HEV_SeriesParallel.exe
70 | delete HEV_Model_Parameters*.mat
71 | delete results*.mat
72 | rtw_dirs = dir('*_rtw');
73 | for i=1:size(rtw_dirs,1)
74 | rmdir(rtw_dirs(i).name,'s');
75 | %eval(['!rmdir ' rtw_dirs(i).name '/S /Q'])
76 | end
77 | !rmdir slprj /S/Q
78 | set_param(bdroot,'SystemTargetFile','grt.tlc')
79 | set_param([bdroot '/SLRT Scope'],'Commented','off');
80 | set_param(bdroot,'DefaultParameterBehavior','Tunable');
81 |
--------------------------------------------------------------------------------
/Workflows/Real_Time/Test_RT_Model.m:
--------------------------------------------------------------------------------
1 | %% OPEN MODEL
2 | mdl = 'HEV_SeriesParallel';
3 | open_system(mdl);
4 | HEVPS_setupRTtest
5 |
6 | %% Get reference results
7 | HEVSP_setdesktop
8 | sim(mdl)
9 | t_ref = Car.time; y_ref = Car.signals(4).values(:,1);
10 | clear tout yout
11 |
12 | %% Create plot for comparing results
13 | figure(1); clf;
14 | temp_colororder = get(gca,'DefaultAxesColorOrder');
15 | set(gcf,'Position',[552 50 472 301]);
16 | plot(t_ref,y_ref,'k','LineWidth',3)
17 | title('Comparing Simulation Results','FontSize',14,'FontWeight','Bold');
18 | xlabel('Time (s)','FontSize',12);ylabel('Vehicle Speed');
19 | legend({'Reference'},'Location','NorthWest')
20 | set(gca,'YLim',[0 20]);
21 |
22 | %% Test with fixed-step, fixed-cost solver
23 | HEVSP_setrealtime
24 | sim(mdl)
25 | t_fs = Car.time; y_fs = Car.signals(4).values(:,1);
26 |
27 | %% Add fixed-step results to plot
28 | figure(1)
29 | hold on
30 | h2=stairs(t_fs, y_fs,'Color',temp_colororder(2,:),'LineWidth',2.5);
31 | hold off
32 | legend({'Reference','Fixed-Step'},'Location','NorthWest')
33 |
34 | %% Enable real-time parameter
35 | tune_bpth = [mdl '/Vehicle Dynamics/Simple/Inertia'];
36 | set_param(tune_bpth,'inertia_conf','runtime');
37 |
38 | % Highlight block
39 | open_system(get_param(tune_bpth,'Parent'),'force')
40 | set_param(tune_bpth,'Selected','on');
41 |
42 | %% Build and download to real-time target
43 | % Choose target
44 | cs = getActiveConfigSet(mdl);
45 | cs.switchTarget('slrealtime.tlc',[]);
46 |
47 | slbuild(mdl);
48 |
49 | %% Download to real-time target
50 | tg = slrealtime;
51 | tg.connect;
52 |
53 | %% Run application
54 | tg.load(mdl)
55 | tg.start('ReloadOnStop',true,'ExportToBaseWorkspace',true)
56 |
57 | open_system(mdl);
58 | disp('Waiting for SLRT to finish...');
59 | pause(1);
60 | while(strcmp(tg.status,'running'))
61 | pause(2);
62 | disp(tg.status);
63 | end
64 | pause(2);
65 |
66 | %% Extract results from logged data in Simulink Data Inspector
67 | y_slrt1 = logsout.getElement('Veh Spd');
68 |
69 | %% Add results from real-time simulation
70 | figure(1)
71 | hold on
72 | h3=stairs(y_slrt1.Values.Time,y_slrt1.Values.Data,'c:','LineWidth',2.5);
73 | hold off
74 | legend({'Reference','Fixed-Step','Real-Time'},'Location','NorthWest');
75 |
76 | %% Change a Simscape block parameter
77 | disp(['Vehicle Mass (current) = ' num2str(getparam(tg,'','HEV_Vehicle_Mass'))]);
78 | setparam(tg,'','HEV_Vehicle_Mass',1800)
79 | disp(['Vehicle Mass (new) = ' num2str(getparam(tg,'','HEV_Vehicle_Mass'))]);
80 |
81 | %% Run simulation with new parameter value
82 | tg.start('ReloadOnStop',true,'ExportToBaseWorkspace',true)
83 |
84 | disp('Waiting for Simulink Real-Time to finish...');
85 | pause(1);
86 | while(strcmp(tg.status,'running'))
87 | pause(2);
88 | disp(tg.status);
89 | end
90 | pause(2);
91 |
92 | %% Extract results from logged data in Simulink Data Inspector
93 | y_slrt2 = logsout.getElement('Veh Spd');
94 |
95 | %% Add results from modified vehicle model
96 | figure(1)
97 | hold on
98 | h4=stairs(y_slrt2.Values.Time,y_slrt2.Values.Data,'Color',temp_colororder(4,:),'LineWidth',2);
99 | hold off
100 | legend({'Reference','Fixed-Step','Real-Time','Modified'},'Location','NorthWest');
101 |
102 | % Copyright 2011-2024 The MathWorks(TM), Inc.
103 | %% Close model and clean up directory
104 | cleanup_rt_dir
105 |
--------------------------------------------------------------------------------
/Workflows/Real_Time/cleanup_rt_dir.m:
--------------------------------------------------------------------------------
1 | % CLEANUP RT TEST DIR
2 | % Copyright 2011-2024 The MathWorks, Inc.
3 |
4 | mdl = 'HEV_SeriesParallel';
5 | delete([mdl '*bio.m'])
6 | delete([mdl '*pt.m'])
7 | delete([mdl '*ref.m'])
8 | delete([mdl '*.dlm'])
9 | delete([mdl '*.xml'])
10 | rtw_dirs = dir('*_rtw');
11 | for i=1:size(rtw_dirs,1)
12 | rmdir(rtw_dirs(i).name,'s');
13 | %eval(['!rmdir ' rtw_dirs(i).name '/S /Q'])
14 | end
15 | bdclose all
16 | delete('*.mex*')
17 |
18 |
--------------------------------------------------------------------------------
/Workflows/Reports/HEV_Model_Report_SHORT.html:
--------------------------------------------------------------------------------
1 |
2 |
3 | HEV Simulation Test ReportHEV Simulation Test Report
Generated by the Simulink Report Generator
Copyright 2024 The MathWorks, Inc.
Kapitel 1. Configuration "Sys BD VS" Tests
% CONFIGURE HEV Model
7 | HEV_Configuration = char(HEV_Config_Set_RPT{HEV_Config_Set_RPT_i});
8 | Select_HEV_Model_Systems(HEV_Configuration, HEV_Configs);
9 | Drive_Cycle_Num = 1;
10 |
11 |
12 |
13 | Model Screenshot
HEV Mode Logic
Kapitel 2. Simulation Time Table
Tabelle 2.1. Simulation Time
| Configuration | Elapsed Time | Simulation Time | Sim Time/Elapsed Time |
|---|
| Sys BD VS | 4.6493 | 195 | 41.9418 |
| Sys BD VS | 6.3675 | 400 | 62.8186 |
| Sys BD VS | 7.5770 | 400 | 52.7915 |
