165 lines
4.4 KiB
Matlab
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 |