MDDP-Cableship/maths/power_model.m
2020-12-01 16:51:10 +00:00

165 lines
4.4 KiB
Matlab

%% power_model.m
%%
%% Vessel power model
close all;clear all;clc;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% Flags
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
CUMULATIVE_ERRORS = false;
ITERATE = true;
SAVE = true;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% Parameters
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
day_to_seconds = 24*60*60;
hours_to_seconds = 60*60;
ITERATIONS = 5;
MIN_P_IN = 0; % W, max power from fuel cells
MAX_P_IN = 8e6; % W, max power from fuel cells
P_IN_LOAD = 0.3; % most efficient load percent
%%%%% DP (SS7)
MAX_P_OUT = 3842e3; % W
MIN_P_OUT = 362e3; % W
TITLE = 'Dyn. Pos. Sea State 7';
SIMULATION_DAYS = 2; % days
cable_drum = get_extra_p(SIMULATION_DAYS*day_to_seconds, 3*hours_to_seconds, 1.5*hours_to_seconds, 946.26e3);
cable_lower = get_extra_p(SIMULATION_DAYS*day_to_seconds, 40*hours_to_seconds, 1.5*hours_to_seconds, 908.41e3);
crane = get_extra_p(SIMULATION_DAYS*day_to_seconds, 12*hours_to_seconds, 15*60, 245.25e3);
rov_launch = get_extra_p(SIMULATION_DAYS*day_to_seconds, 8*hours_to_seconds, 20*60, 454.2e3);
EXTRA_P = [cable_drum ; cable_lower ; crane ; rov_launch];
%%%%% Outbound
% MAX_P_OUT = 1600e3; % W
% MIN_P_OUT = 600e3; % W
% TITLE = 'Outbound Steaming';
% SIMULATION_DAYS = 3; % days
%%%%% Manouvering
% MAX_P_OUT = 800e3; % W
% MIN_P_OUT = 200e3; % W
% TITLE = 'Manouvering';
% SIMULATION_DAYS = 1; % days
%
% %%%%% Home
% MAX_P_OUT = 800e3; % W
% MIN_P_OUT = 200e3; % W
% TITLE = 'Homebound';
% SIMULATION_DAYS = 3; % days
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% Simulate
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
p_av = (MAX_P_OUT + MIN_P_OUT) / 2;
if exist('EXTRA_P')
[power_in,battery_level,power_out,unused_energy,unavailable_energy, batt_capacity] = power_sim(MAX_P_OUT, MIN_P_OUT, MAX_P_IN, MIN_P_IN, SIMULATION_DAYS, p_av, p_av, -1, CUMULATIVE_ERRORS, EXTRA_P);
else
[power_in,battery_level,power_out,unused_energy,unavailable_energy, batt_capacity] = power_sim(MAX_P_OUT, MIN_P_OUT, MAX_P_IN, MIN_P_IN, SIMULATION_DAYS, p_av, p_av, -1, CUMULATIVE_ERRORS);
end
sim_seconds = length(power_in);
x = (1:sim_seconds) / (60 * 60);
x_ticks = (1: sim_seconds / (60 * 60));
if SIMULATION_DAYS > 1
if SIMULATION_DAYS < 4
x_ticks = (1: sim_seconds / (60 * 60));
else
x_ticks = (1: sim_seconds / (60 * 60 * 24));
end
x = x / 24;
end
figure('Renderer', 'painters', 'Position', [10 10 1000 800])
line_width = 1;
subplot(3, 1, 1);
sgtitle(TITLE);
hold on;
grid on;
plot(x, power_in / 1e6, 'g', 'LineWidth', 2);
plot(x, power_out / 1e6, 'r', 'LineWidth', 1);
max_line = yline(MAX_P_OUT / 1e6, '-c', 'LineWidth', line_width * 0.75);
min_line = yline(MIN_P_OUT / 1e6, '-c', 'LineWidth', line_width * 0.75);
max_line.Alpha = 0.5;
min_line.Alpha = 0.5;
yline(p_av / 1e6, '--m', 'LineWidth', line_width * 0.5);
legend('P In', 'P Out', 'Max P Out', 'Min P Out', 'Average P In');
ylabel('Power (MW)')
xlim([0 inf])
% ylim([0 ceil(max(max(power_in/1e6), max(power_out/1e6)))])
xticks(x_ticks)
if SIMULATION_DAYS > 1
xlabel('Time (Days)')
else
xlabel('Time (Hours)')
end
hold off;
% figure(2)
subplot(3, 1, 2);
hold on;
grid on;
plot(x, battery_level * 100 / batt_capacity, 'LineWidth', 2);
legend('Battery Level');
ylabel('Capacity (%)')
xlim([0 inf])
ylim([0 100])
xticks(x_ticks)
if SIMULATION_DAYS > 1
xlabel('Time (Days)')
else
xlabel('Time (Hours)')
end
hold off;
subplot(3, 1, 3);
hold on;
grid on;
plot(x, unused_energy, 'g', 'LineWidth', line_width);
plot(x, unavailable_energy, 'r', 'LineWidth', line_width);
legend('Unused', 'Unavailable');
if CUMULATIVE_ERRORS
ylabel('Energy (J)')
else
ylabel('Power (W)')
end
xlim([0 inf])
ylim([0 inf])
xticks(x_ticks)
if SIMULATION_DAYS > 1
xlabel('Time (Days)')
else
xlabel('Time (Hours)')
end
hold off;
if SAVE
exportgraphics(gcf, sprintf('%s.png', TITLE), 'Resolution', '250', 'ContentType','vector')
end
% FINAL STATS
if CUMULATIVE_ERRORS
fprintf('%.f MJ/day of unused power\n', unused_energy(end) / (1e6 * SIMULATION_DAYS));
fprintf('%.f MJ/day of unavailable power\n\n', unavailable_energy(end) / (1e6 * SIMULATION_DAYS));
fprintf('%.f MJ of unused power\n', unused_energy(end) / 1e6);
fprintf('%.f MJ of unavailable power\n', unavailable_energy(end) / 1e6);
end