HD: Add outputs of Morison mesh nodal loads when OutAll is true

Author: luwang00

docs/source/user/hydrodyn/input_files.rst

@@ -899,9 +899,10 @@ HydroDyn standalone, or by the OpenFAST program when running a coupled
 simulation. See :numref:`hd-summary-file` for summary file details.
 
 If **OutAll** is set to TRUE, HydroDyn will output the total strip-theory
-forces and moments on each user-defined member and joint. These additional
-output channels are inserted as additional columns in the output file
-independent of any user-requested outputs. See :ref:`hd-outall-option`
+forces and moments on each user-defined member and joint, followed by the
+forces and moments on each computational node of the strip-theory mesh.
+These additional output channels are inserted as additional columns in the
+output file independent of any user-requested outputs. See :ref:`hd-outall-option`
 for more information.
 
 If **OutSwtch** is set to 1, outputs are sent to a file with the name

docs/source/user/hydrodyn/output_files.rst

@@ -174,8 +174,8 @@ data for a single 6\ **NBody**-by-6\ **NBody** matrix. If **NBodyMod** > 1
 If **OutAll** is set to TRUE, HydroDyn will output the total strip-theory
 forces and moments on each user-defined member and joint. These are
 included as additional columns in the output file independent of any
-user-requested outputs. The forces and moments on the members (loads on
-the side wall) will be written first. For example, the 6 load components
+user-requested outputs. The forces and moments on the members (integrated loads across
+all side walls) will be written first. For example, the 6 load components
 on the first member in the **MEMBERS** table (the first row of the table)
 are given by **M1TotFxi**, **M1TotFyi**, **M1TotFzi**, **M1TotMxi**,
 **M1TotMyi**, and **M1TotMzi**. After the member loads, the total lumped
@@ -197,7 +197,27 @@ force and moment components are resolved in the earth-fixed inertial frame
 of reference, and the moments are computed about the instantaneous principal
 reference point (PRP), same as the output channels **HydroFxi**,
 **HydroFyi**, etc. As a reminder, the PRP is a body-fixed point located at
-the earth-fixed origin when the HydroDyn structure is undisplaced.
+the earth-fixed origin when the HydroDyn structure is undisplaced. Summing
+all member and joint loads gives the total strip-theory loads on the entire
+structure.
+
+After the member and joint loads, HydroDyn also outputs the total forces
+and moments on each computational node of the HydroDyn strip-theory (Morison)
+mesh. This internal mesh is used to map the loads to other structural
+modules, such as SubDyn, and contains joint nodes at the user-defined joint
+locations followed by member internal nodes created from member discretization.
+The force and moment components are again resolved in the earth-fixed inertial
+frame of reference. However, the moment on each node is about the node itself,
+not about the PRP as with the member and joint load outputs above. Additionally,
+the joint mesh nodes can have load contributions from both member side walls
+and from the joint/endplates. This is because part of the side-wall loads on
+the first and last element of a member can be assignd to the joint nodes. As
+a result, the load outputs at the joint nodes do not necessarily match the joint
+load outputs above, which do not contain contributions from member side walls.
+The nodal load output column names indicate the node number, e.g., **N1TotFxi**,
+**N1TotFyi**, **N1TotFzi**, **N1TotMxi**, **N1TotMyi**, and **N1TotMzi** for
+the first node. The node numbering follows the **Nodes** table in the HydroDyn
+summary file.
 
 Results File
 ~~~~~~~~~~~~

modules/hydrodyn/src/Morison_Output.f90

@@ -7936,6 +7936,11 @@ SUBROUTINE MrsnOut_MapOutputs( y, p, u, m )
       y%WriteOutput(I+1:I+6) = m%F_tot_End(:,im)
       I = I + 6_IntKi
    END DO
+   DO im = 1,p%NNodes
+      y%WriteOutput(I+1:I+3) = y%Mesh%Force(:,im)
+      y%WriteOutput(I+4:I+6) = y%Mesh%Moment(:,im)
+      I = I + 6_IntKi
+   END DO
 
 END SUBROUTINE MrsnOut_MapOutputs
 
@@ -8957,7 +8962,7 @@ SUBROUTINE SetOutParam(OutList, p, ErrStat, ErrMsg )
    ! If a selected output channel is not available in this module, set error flag.
    !-------------------------------------------------------------------------------------------------
    if (p%OutAll) then
-      ALLOCATE ( p%OutParam(1:(p%NumOuts + 6_IntKi*p%NMembers + 6_IntKi*p%NJoints)) , STAT=ErrStat2 )
+      ALLOCATE ( p%OutParam(1:(p%NumOuts + 6_IntKi*p%NMembers + 6_IntKi*p%NJoints + 6_IntKi*p%NNodes)) , STAT=ErrStat2 )
    else
       ALLOCATE ( p%OutParam(1:p%NumOuts) , STAT=ErrStat2 )
    endif
@@ -9042,6 +9047,25 @@ SUBROUTINE SetOutParam(OutList, p, ErrStat, ErrMsg )
             p%OutParam(k)%SignM =  1_IntKi
          end do
       end do
+      do i = 1,p%NNodes
+         p%OutParam(k+1)%Name  = "N"//trim(num2lstr(i))//"TotFxi"
+         p%OutParam(k+1)%Units = "(N)"
+         p%OutParam(k+2)%Name  = "N"//trim(num2lstr(i))//"TotFyi"
+         p%OutParam(k+2)%Units = "(N)"
+         p%OutParam(k+3)%Name  = "N"//trim(num2lstr(i))//"TotFzi"
+         p%OutParam(k+3)%Units = "(N)"
+         p%OutParam(k+4)%Name  = "N"//trim(num2lstr(i))//"TotMxi"
+         p%OutParam(k+4)%Units = "(N-m)"
+         p%OutParam(k+5)%Name  = "N"//trim(num2lstr(i))//"TotMyi"
+         p%OutParam(k+5)%Units = "(N-m)"
+         p%OutParam(k+6)%Name  = "N"//trim(num2lstr(i))//"TotMzi"
+         p%OutParam(k+6)%Units = "(N-m)"
+         do j = 1,6
+            k = k + 1_IntKi
+            p%OutParam(k)%Indx  = -1_IntKi
+            p%OutParam(k)%SignM =  1_IntKi
+         end do
+      end do
    end if
 
    RETURN