4.8.3.50. GET_RFLEX_NODEVEL

GET_RFLEX_NODEVEL subroutine returns a velocity vector for a node of RFlex body. This is an auxiliary subroutine for Modal_Force_Ext.

Table 4.145 Function Name

Language type	Subroutine
FORTRAN	call get_rflex_nodevel(ifbody, NodeSeq, MKID, nMK, VEL, ErrFlg)
C/C++	get_rflex_nodevel(ifbody, NodeSeq, MKID, nMK, VEL, &ErrFlg)

Table 4.146 Parameter information

Variable Name	Size	Description
ifbody	int	Sequential id of RFlex body defined in RecurDyn/Solver. This is a related argument with the 5th argument of Modal_Force_Ext subroutine.
NodeSeq	int	Node sequential id defined in RecurDyn/Solver. This is a related argument with the 6th argument of Modal_Force_Ext subroutine.
MKID	int[2]	An array of integer type. Each value should be zero or a marker id. 1st value is defined as base marker. 2nd value is defined as reference marker.
nMK	int	An integer variable for considering base and reference marker. If nMK is 0, then the subroutine returns a global position vector and an orientation matrix of NodeSeq. If nMK is 1, then the subroutine calculates a position vector considering base marker. If nMK is 2, then the subroutine calculates a position vector considering base and reference marker.
Vel	double[6]	An array of double precision type. The array size must be 6. First 3 values mean a relative velocity vector. The last 3 values mean a relative angular velocity vector.
Errflg	int	Error flag. If the result of this argument is -1 (means TRUE in Fortran logical value), there is no error. The others mean that there is an error.

Listing 4.71 C/C++ code for GET_RFLEX_NODEVEL

 #include "stdafx.h"
 #include "DllFunc.h"
 #include <stdio.h>
 FILE* NodeVELwrite;

 RecurDyn_UserSubRoutineWizard9_API void __cdecl modal_force_ext
 (int id, double time, double upar[], int npar, int ifbody, int nodarr[], int nonde, int jflag, int iflag, double result[])
 {
   using namespace rd_syscall;
   // Parameter Information
   //   id     : Modal force sequential identification. (Input)
   //   time   : Simulation time of RD/Solver. (Input)
   //   upar   : Parameters defined by user. (Input)
   //   npar   : Number of user parameters. (Input)
   //   ifbody : RFLEX Body sequential ID. (Input)
   //   nodarr : Node ID array of input node set. (Input)
   //   nonde  : Number of node of node set. (Input)
   //   jflag  : When RD/Solver evaluates a Jacobian, the flag is true. (Input)
   //   iflag  : When RD/Solver initializes arrays, the flag is true. (Input)
   //   result : Returned nodal force vector. Acting point of the nodal force is each center of each node.
   //            Reference frame of each force vector must be Ground.InertiaMarker. (Output, Size: nonde * 6)

   // User Statement
   int errflg =0 , ifinish = 0;
   int nodeseq = 0;
   double vel[6],tvel[3],rvel[3];

   if (iflag)
   {
       for(int k = 0 ; k <nonde ; k++)
       {
         NodeVELwrite=fopen("RFlexNodeVEL_C.txt","w");
         fprintf(NodeVELwrite,"RFlex Node%d Velocity \n",nodarr[k]);

         for (int i = 0; i<6 ;i++) { vel[i]=0.0; }
         for (int i = 0; i<3 ;i++) { tvel[i]=0.0; rvel[i]=0.0; }

         get_rflex_nodeseqid(ifbody,nodarr[k],&nodeseq,&errflg);
         get_rflex_nodevel(ifbody,nodeseq,NULL,0,vel,&errflg);
         get_rflex_nodetvel(ifbody,nodeseq,NULL,0,tvel,&errflg);
         get_rflex_nodervel(ifbody,nodeseq,NULL,0,rvel,&errflg);

         fprintf(NodeVELwrite,"USING GET_RFLEX_NODEVEL \n");
         fprintf(NodeVELwrite,"INITIAL TRANSLATIONAL VELOCITY \n");
         fprintf(NodeVELwrite,"VX = %20.10e\n",vel[0]);
         fprintf(NodeVELwrite,"VY = %20.10e\n",vel[1]);
         fprintf(NodeVELwrite,"VZ = %20.10e\n\n",vel[2]);
         fprintf(NodeVELwrite,"INITIAL ROTATIONAL VELOCITY \n");
         fprintf(NodeVELwrite,"WX = %20.10e\n",vel[3]);
         fprintf(NodeVELwrite,"WY = %20.10e\n",vel[4]);
         fprintf(NodeVELwrite,"WZ = %20.10e\n\n",vel[5]);

         fprintf(NodeVELwrite,"USING GET_RFLEX_NODETVEL \n");
         fprintf(NodeVELwrite,"INITIAL TRANSLATIONAL VELOCITY \n");
         fprintf(NodeVELwrite,"VX = %20.10e\n",tvel[0]);
         fprintf(NodeVELwrite,"VY = %20.10e\n",tvel[1]);
         fprintf(NodeVELwrite,"VZ = %20.10e\n\n",tvel[2]);

         fprintf(NodeVELwrite,"USING GET_RFLEX_NODERVEL \n");
         fprintf(NodeVELwrite,"INITIAL ROTATIONAL VELOCITY \n");
         fprintf(NodeVELwrite,"WX = %20.10e\n",rvel[0]);
         fprintf(NodeVELwrite,"WY = %20.10e\n",rvel[1]);
         fprintf(NodeVELwrite,"WZ = %20.10e\n\n",rvel[2]);

       }
   }
   for(int i=0;i<6*nonde;i++)
   {
       result[i] = 0.0;
   }
 }