Engine Subroutine SUSERnn for Solid Elements

This subroutine calculates the internal forces and moments.

Arguments

The argument list of RUSERnn is:

C--------------------------------------------------------------------------
SUBROUTINE SUSER29(
     1 NEL    ,NUVAR  ,IOUT   ,IPROP  ,IMAT  ,SOLID_ID,TIME  ,TIMESTEP,
     2 EINT   ,VOL    ,UVAR   ,FR_WAVE,OFF    ,RHO    ,SIG    ,
     3 XX1    ,XX2    ,XX3    ,XX4    ,XX5    ,XX6    ,XX7    ,XX8    , 
     4 YY1    ,YY2    ,YY3    ,YY4    ,YY5    ,YY6    ,YY7    ,YY8    ,  
     5 ZZ1    ,ZZ2    ,ZZ3    ,ZZ4    ,ZZ5    ,ZZ6    ,ZZ7    ,ZZ8    ,
     6 UX1    ,UX2    ,UX3    ,UX4    ,UX5    ,UX6    ,UX7    ,UX8    ,
     7 UY1    ,UY2    ,UY3    ,UY4    ,UY5    ,UY6    ,UY7    ,UY8    ,
     8 UZ1    ,UZ2    ,UZ3    ,UZ4    ,UZ5    ,UZ6    ,UZ7    ,UZ8    ,
     9 VX1    ,VX2    ,VX3    ,VX4    ,VX5    ,VX6    ,VX7    ,VX8    ,
     A VY1    ,VY2    ,VY3    ,VY4    ,VY5    ,VY6    ,VY7    ,VY8    ,
     B VZ1    ,VZ2    ,VZ3    ,VZ4    ,VZ5    ,VZ6    ,VZ7    ,VZ8    ,
     C VRX1   ,VRX2   ,VRX3   ,VRX4   ,VRX5   ,VRX6   ,VRX7   ,VRX8   ,
     D VRY1   ,VRY2   ,VRY3   ,VRY4   ,VRY5   ,VRY6   ,VRY7   ,VRY8   ,
     E VRZ1   ,VRZ2   ,VRZ3   ,VRZ4   ,VRZ5   ,VRZ6   ,VRZ7   ,VRZ8   ,
     F FX1    ,FX2    ,FX3    ,FX4    ,FX5    ,FX6    ,FX7    ,FX8    ,
     G FY1    ,FY2    ,FY3    ,FY4    ,FY5    ,FY6    ,FY7    ,FY8    ,
     H FZ1    ,FZ2    ,FZ3    ,FZ4    ,FZ5    ,FZ6    ,FZ7    ,FZ8    ,
     I MX1    ,MX2    ,MX3    ,MX4    ,MX5    ,MX6    ,MX7    ,MX8    ,
     J MY1    ,MY2    ,MY3    ,MY4    ,MY5    ,MY6    ,MY7    ,MY8    ,
     K MZ1    ,MZ2    ,MZ3    ,MZ4    ,MZ5    ,MZ6    ,MZ7    ,MZ8    ,
     L STIFM  ,STIFR  ,VISCM  ,VISCR  )
C-------------------------------------------------------------------------

Definitions

Argument	Format	Description
`NEL`	Integer read only scalar	Number of elements per group. In Radioss Engine subroutines, the element data are treated by groups for vectorization. This argument is machine-dependent and set by Radioss.
`NUVAR`	Integer scalar	Number of extra variables needed for each solid element in `SINInn` or `SUESERnn`
`IOUT`	Integer read only scalar	Output file unit (L01 file)
`IPROP`	Integer read only scalar	Property number.
`IMAT`	Integer read only scalar	Material number.
`SOLID_ID`	Integer read only scalar	Array of size `NEL` containing solid element ID.
`Time`	Float read only	Current time.
`TIMESTEP`	Float read only	Current time step.
`EINT`	Float array read only	Array of size `NEL` containing current total internal energy.
`VOL`	Float array read only	Array of size `NEL` containing initial element volumes. The incoming value is computed with one integration point, which can be recomputed.
`UVAR`	Float array read write	Array of size `NEL`*`NUVAR` containing user element variables.
`OFF`	Float array read write	Array of size `NEL` containing deleted element flags. 0 If the element is OFF 1 If the element is ON
`RHO`	Float array read write	Array of size `NEL` containing current densities.
`SIG`	Float array read write	Array of size 6*`NEL` containing current stress tensor SX, SY, SZ, SXY, SYZ, and SZX.
`XX1`	Float array read only	Array of size `NEL` containing X coordinate at node 1 in the global frame at time TIME.
`UX1`	Float array read only	Array of size `NEL` containing X displacement at node 1 in the global frame at time TIME.
`VX1`	Float array read only	Array of size `NEL` containing initial X velocity at node 1 in the global frame at time TIME.
`VRX1`	Float array read only	Array of size `NEL` containing initial X rotational velocity at node 1 in the global frame at time TIME.
`FX1`	Float array write only	Array of size NEL containing X force at node 1 at time TIME.
`MX1`	Float array write only	Array of size `NEL` containing X moment at node 1 at time TIME.
`STIFM`	Float array write only	Array of size `NEL` containing element translational stiffness in time step computation.
`STIFR`	Float array write only	Array of size `NEL` containing element rotational stiffness in time step computation.
`VISCM`	Float array write only	Array of size `NEL` containing element translational viscosity in time step computation.
`VISCR`	Float array write only	Array of size `NEL` containing element rotational viscosity in time step computation.

Comments

STIFM, STIFR, VISCM, and VISCR compute the nodal or element time step. MASS and XINER compute element time step.
RHO and SIG are used for post-processing in time history and animation. The modification of these variables effects only output. RHO and SIG retain initial values, unless they are modified.