Update beams self-weight (only sensor outputs)

Author: RBergua

modules/subdyn/src/SubDyn_Output.f90

@@ -247,10 +247,9 @@ SUBROUTINE ConfigOutputNode_MKF_ID(pLst, iElem, iiNode, iStore, NodeID2)
       CALL ElemM(p%ElemProps(iElem),         pLst%Me(:,:,iiNode,iStore))
       CALL ElemK(p%ElemProps(iElem),         pLst%Ke(:,:,iiNode,iStore))
       CALL ElemF(p%ElemProps(iElem), Init%g, pLst%Fg(:,iiNode,iStore), FCe)
-      ! NOTE: Removing this force contribution for now 
-      ! The output of subdyn will just be the "Kx" part for now
-      !pLst%Fg(:,iiNode,iStore) = pLst%Fg(:,iiNode,iStore) + FCe(1:12) ! Adding cable element force 
-      pLst%Fg(:,iiNode,iStore) = FCe(1:12) ! Adding cable element force 
+      ! Apply superposition correction: include both the element gravity force vector 
+      ! and the cable pretension nodal force when evaluating internal member loads.
+      pLst%Fg(:,iiNode,iStore) = pLst%Fg(:,iiNode,iStore) + FCe(1:12) ! gravity force + cable element force
    END SUBROUTINE ConfigOutputNode_MKF_ID