divide power by 8

symbiotic-engineering · Nov 8, 2023 · 7d4dd94 · 7d4dd94
1 parent ffe7f96
commit 7d4dd94
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/mdocean/simulation/modules/dynamics/dynamics.m b/mdocean/simulation/modules/dynamics/dynamics.m
@@ -39,7 +39,7 @@
     X_unsat = get_response(w,m,b,k,Fd);
 
     % confirm unsaturated response doesn't exceed maximum capture width
-    P_unsat = 1/2 * in.B_p * w.^2 .* X_unsat.^2;
+    P_unsat = 1/8 * 1/2 * in.B_p * w.^2 .* X_unsat.^2;
 
     F_ptrain_over_x = sqrt( (in.B_p * w).^2 + (K_p).^2 );
     F_ptrain_unsat = F_ptrain_over_x .* X_unsat;
@@ -54,7 +54,7 @@
     X_sat = get_response(w,m,b_sat,k_sat,Fd);
 
     % calculate power
-    P_matrix = 1/2 * (mult * in.B_p) .* w.^2 .* X_sat.^2;
+    P_matrix = 1/8 * 1/2 * (mult * in.B_p) .* w.^2 .* X_sat.^2;
 
     X_max = max(X_sat,[],'all');
     h_s_extra = (in.h_s - in.T_s - (in.h_f - in.T_f) - X_max) / in.h_s; % extra height on spar after accommodating float displacement