add back min damping plate diameter constraint

symbiotic-engineering · Oct 28, 2024 · 02883db · 02883db
1 parent 28606b1
commit 02883db
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 2 deletions.
diff --git a/mdocean/inputs/parameters.m b/mdocean/inputs/parameters.m
@@ -48,7 +48,8 @@
             'theta_dt', atan(17.5/15),...       % angle from horizontal of damping plate support tubes (rad)
             'FOS_min',  1.5,...                 % minimum FOS (-)	
             'FCR',      0.113,...               % fixed charge rate (-), see RM3 report p63
-            'N_WEC',    100,...                 % number of WECs in array (-)   
+            'N_WEC',    100,...                 % number of WECs in array (-)
+            'D_d_min',  30,...                  % minimum damping plate diameter
             'D_d_over_D_s', 30/6,...            % normalized damping plate diameter (-)
             'T_s_over_D_s', 29/6,... %35/6      % normalized spar draft (-)
             'h_d_over_D_s', 0.1/6,...%in2m/6    % normalized damping plate thickness (-)     

diff --git a/mdocean/simulation/simulation.m b/mdocean/simulation/simulation.m
@@ -63,7 +63,8 @@
 g(8) = FOS3Y / p.FOS_min - 1;           % damping plate survives max force
 g(9) = FOS_buckling / p.FOS_min - 1;    % spar survives max force in buckling
 g(10) = P_elec;                         % positive power
-g(11) = 0; %1 + min(Kp_over_Ks,[],'all');   % spar heave stability (positive effective stiffness)
+g(11) = in.D_d / p.D_d_min - 1;         % damping plate diameter (spar natural freq)
+%1 + min(Kp_over_Ks,[],'all');   % spar heave stability (positive effective stiffness)
 g(12) = h_s_extra(1);                   % prevent float rising above top of spar
 g(13) = h_s_extra(2);                   % prevent float going below bottom of spar
 g(14) = p.LCOE_max/LCOE - 1;            % prevent more expensive than threshold