Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Constant geometric variables #50

Open
wants to merge 2 commits into
base: main
Choose a base branch
from
Open

Constant geometric variables #50

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
3a5451a
hard code geometric variables in least-cost case
MadisonDietrich Apr 25, 2024
99035a0
code for graph that overlays objective values when geometric variable…
MadisonDietrich May 8, 2024
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
23 changes: 22 additions & 1 deletion mdocean/optimization/multiobjective/pareto_curve_heuristics.m
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -250,7 +250,28 @@ function pareto_curve_heuristics()
grid on
xlabel('Percent along the Pareto Curve')
ylabel('Objective Value')
pct = [0 4.6683 16.1186 27.2305 28.6278 30.4717 34.0738 38.4791 40.3735 41.0267 ...
41.9764 43.3829 43.5881 44.4236 44.8667 45.7298 50.1799 51.8743 52.9700 55.1960 ...
56.1828 58.0417 58.8600 60.9711 61.2873 61.5743 62.1110 62.7047 63.3152 64.3656 ...
65.9389 66.9220 67.6261 68.5065 68.7219 69.0767 70.5319 70.8866 71.3303 72.1607 ...
72.2847 73.9855 75.2496 78.0466 78.7258 80.3724 83.5218 84.3179 84.6826 86.4844 ...
87.0652 88.0020 89.1332 89.8590 91.8756 96.5597 97.7855 99.0541]';
lcoe = [10.0000 14.2494 14.9043 16.3479 16.5465 16.8356 17.5015 18.5171 19.0298 ...
19.2161 19.4955 19.9336 20.0000 20.2759 20.4259 20.7251 22.4162 23.1281 23.6075 ...
24.6018 25.0557 25.9355 26.3336 27.3902 27.5518 27.6990 27.9760 28.2847 28.6044 ...
29.1591 30.0000 30.5316 30.9152 31.3981 31.5409 31.7360 32.5223 32.7232 32.9696 ...
33.4381 33.5070 34.4715 35.1936 36.8083 37.2151 38.1557 40.0000 40.4700 40.6851 ...
41.7548 42.1005 42.6605 43.3483 43.7747 44.9983 47.8588 48.6169 49.4073]';
pvar = [158.8509 151.7591 134.3300 117.4241 115.2494 112.3817 106.8146 100.1354 ...
97.3378 96.3835 95.0058 93.0008 92.7122 91.5459 90.9332 89.7518 83.9438 81.8848 ...
80.6024 78.0478 76.9518 74.9637 74.1240 72.0641 71.7683 71.5025 71.0121 70.4791 ...
69.9409 69.0366 67.7310 66.9463 66.3994 65.7339 65.6954 65.4326 64.2794 64.0486 ...
63.7328 63.1846 63.0969 62.0126 61.2449 59.6798 59.4164 58.4211 56.8669 56.5049 ...
56.3351 55.5597 55.3120 54.9383 54.6481 54.2157 53.5928 51.9609 51.5273 51.1327]';
plot(pct,lcoe,'Color',[1 0 0])
plot(pct,pvar,'Color',[0 0 1])
improvePlot
cent = char(0162);
legend(['LCOE (' cent '/kWh)'],'c_v (%)',Location='northeast')
legend(['LCOE (' cent '/kWh)'],'c_v (%)', ...
'LCOE, const geom','Pvar, const geom',Location='northeast')
end
4 changes: 4 additions & 0 deletions mdocean/simulation/simulation.m
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -6,9 +6,13 @@
in = p;

in.D_f = X(1); % outer diameter of float (m)
%in.D_f = 20.8418;
D_s_over_D_f = X(2); % normalized diameter of spar column (-)
%D_s_over_D_f = 0.2879;
h_f_over_D_f = X(3); % normalized vertical thickness of float (-)
%h_f_over_D_f = 0.1;
T_s_over_h_s = X(4); % normalized spar draft below waterline (-)
%T_s_over_h_s = 0.4604;
in.F_max = X(5)*1e6; % max powertrain force (N)
in.B_p = X(6)*1e6; % controller (powertrain) damping (Ns/m)
in.w_n = X(7); % controller (powertrain) natural frequency (rad/s)
Expand Down
Loading