From 56c63e86cf54922a8dca39fa4cd8e19b9a9a24d1 Mon Sep 17 00:00:00 2001 From: West Date: Wed, 23 Nov 2022 15:56:48 -0700 Subject: [PATCH] Fixed a bug on line 1265 of MoorDyn_Line.f90 that lead to bending moment for rod EndB to be computed incorrectly. This fix makes the fortran version of MoorDyn v2 consistent with the C++ version. --- modules/moordyn/src/MoorDyn_Rod.f90 | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/modules/moordyn/src/MoorDyn_Rod.f90 b/modules/moordyn/src/MoorDyn_Rod.f90 index b9ac986806..dca892ffff 100644 --- a/modules/moordyn/src/MoorDyn_Rod.f90 +++ b/modules/moordyn/src/MoorDyn_Rod.f90 @@ -827,7 +827,7 @@ SUBROUTINE Rod_DoRHS(Rod, m, p) Rod%Mext = Rod%Mext + (/ Mtemp*sinBeta, -Mtemp*cosBeta, 0.0_DbKi /) ! axial drag - Rod%Dq(:,I) = Rod%Dq(:,I) + VOF * 0.25* Pi*Rod%d*Rod%d * p%rhoW*Rod%CdEnd * MagVq * Vq + Rod%Dq(:,I) = Rod%Dq(:,I) + 0.5 * VOF * 0.25* Pi*Rod%d*Rod%d * p%rhoW*Rod%CdEnd * MagVq * Vq ! >>> what about rotational drag?? <<< eqn will be Pi* Rod%d**4/16.0 omega_rel?^2... *0.5 * Cd... @@ -857,7 +857,7 @@ SUBROUTINE Rod_DoRHS(Rod, m, p) Rod%Mext = Rod%Mext + (/ Mtemp*sinBeta, -Mtemp*cosBeta, 0.0_DbKi /) ! axial drag - Rod%Dq(:,I) = Rod%Dq(:,I) + VOF * 0.25* Pi*Rod%d*Rod%d * p%rhoW*Rod%CdEnd * MagVq * Vq + Rod%Dq(:,I) = Rod%Dq(:,I) + 0.5 * VOF * 0.25* Pi*Rod%d*Rod%d * p%rhoW*Rod%CdEnd * MagVq * Vq ! long-wave diffraction force Rod%Aq(:,I) = Rod%Aq(:,I) + VOF * p%rhoW* Rod%CaEnd * (2.0/3.0*Pi*Rod%d**3 /8.0) * aq