From 7d4dd948ad67df671aca66d026e197787eed4b36 Mon Sep 17 00:00:00 2001 From: rebeccamccabe Date: Wed, 8 Nov 2023 16:18:37 -0500 Subject: [PATCH] divide power by 8 --- mdocean/simulation/modules/dynamics/dynamics.m | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/mdocean/simulation/modules/dynamics/dynamics.m b/mdocean/simulation/modules/dynamics/dynamics.m index d4431c36..599a3ee0 100644 --- 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