/ANIM/BRICK/Restype

Engine Keyword Generates animation files containing brick element data for the specified result. Options used for BRICK element type.

Format

/ANIM/Eltyp/Restype

Definitions

Field	Content	SI Unit Example
`Eltyp`	Element type. ELEM BRICK Brick elements
`Restype`	AMS Elements using AMS timestep due to /DT/CST_AMS. COLOR 1vfrac1 + 2vfrac2 + .... + last*vfracN (LAW51 and LAW151 only). DENS Density DT Element time step. 10 ENER Specific energy (internal energy divided by the element mass). EINT Internal energy. EPSP Plastic strain $ε_{p}$ . EPSD Equivalent strain rate in bricks and in shells too (only available in case of strain rate filtering). HOURG Hourglass energy P Pressure VONM von Mises stress SIGEQ Equivalent stress, based on the yield criteria used for the corresponding material. 3 SIGX Stress XX 5 SIGY Stress YY SIGZ Stress ZZ SIGXY Shear stress XY SIGYZ Shear stress YZ SIGZX Shear stress ZX TEMP Temperature THIC Thickness ORTHDIR/I Euler angles ( $ψ, θ, and ϕ$ in unit $[\deg]$ in layer I). ORTHDIR/ALL Euler angles ( $ψ, θ, and ϕ$ in unit $[\deg]$ in all layers) Defining orthotropic reference system with global reference system - only for orthotropic solid properties. 2 FLAY Number of failed layers for solids. DAM1 , DAM2 or DAM3 Damage in direction 1, 2 or 3. 9 USRi User law variable ( $i$ =1 to 18). USRII User law variable on each integration points (II =1 to 99 number of the variable). 6 Tip: Since the number of user variables that can be output is limited, use /H3D/SOLID/USER/ instead. FILL Filling percentage. This option is used only with `Eltyp` =BRICK for /INIBRI/FILL. K Turbulent energy for ALE material law. TVIS Turbulent viscosity for ALE material law. VORTX Vorticity around x for QUAD (2D) and SOLID (3D) with ALE material law. VORTY Vorticity around Y for SOLID with ALE material law. VORTZ Vorticity around Z for SOLID with ALE material law. VORT Vorticity around for SOLID with ALE material law. QVIS Displays artificial viscous pressure VFRAC Volumetric fractions (for ALE multi-material laws: LAW20, LAW37 and LAW51). BFRAC Burn fraction (for JWL EOS: LAW5 and LAW51). SSP Sound speed. Only available with ALE material laws. SCHLIEREN Schlieren image (optical method widely used in CFD field). Only available with ALE material laws. 7 VOLU Volume of each phase. LAW51 Display results for all Sub-material or specific Sub-material. 8 TDET Detonation times for high explosive JWL EOS. VEL Cell velocity for FVM with /INTER/TYPE22 for all components and magnitude. To request each separately use: VELX, VELY, VELZ, VELXY, VELYZ, VELXZ, \|VEL\|. FIN Cell internal Forces for FVM with /INTER/TYPE22 for all components and magnitude. To request each separately use: FINX, FINY, FINZ, FINXY, FINYZ, FINXZ, \|FIN\|. MOM Cell momentum Density (MOM) for FVM with /INTER/TYPE22 for all components and magnitude. To request each separately use: MOMX, MOMY, MOMZ, MOMXY, MOMYZ, MOMXZ, \|MOM\|. MACH Mach number (/MAT/LAW151 (MULTIFLUID) only). OFF Element status. Where, = -1 Element is not active (it is defined in an activated rigid body). = 0 Deleted element. Between 0 and 1 Under failure process. = 1 Active element.
	WPLA Plastic work for /MAT/LAW12 (3D_COMP) and /MAT/LAW25 (COMPSH)

Comments

For brick elements, the stresses are output in the elemental (corotational) coordinate system for /PROP/TYPE14 (SOLID) elements, and in the orthotropic material coordinate system for /PROP/TYPE6 (SOL_ORTH) elements, when corotational formulation is used. For all other cases the stresses are output in the global coordinate system.
Using /ANIM/BRICK/ORTHDIR with properties 6, 21, and 22 will output 3 real values of angle $ψ, θ, and ϕ$ in unit [deg]. Defining rotation matrix $R$ , to go from global reference system to orthotropic reference system:
- Rotation matrices
- A rotation of $ψ$ radians about the x-axis is defined as:(1)
  $R_{x} (ψ) = [\begin{matrix} 1 & 0 & 0 \\ 0 & cos ψ & - sin ψ \\ 0 & sin ψ & cos ψ \end{matrix}]$
- Similarly, a rotation of $θ$ radians about the y-axis is defined as:(2)
  $R_{y} (θ) = [\begin{matrix} cos θ & 0 & sin θ \\ 0 & 1 & 0 \\ - sin θ & 0 & cos θ \end{matrix}]$
- Finally, a rotation of $ϕ$ radians about the z-axis is defined as:(3)
  $R_{z} (ϕ) = [\begin{matrix} cos ϕ & - sin ϕ & 0 \\ sin ϕ & cos ϕ & 0 \\ 0 & 0 & 1 \end{matrix}]$
- The angles $ψ, θ, and ϕ$ are the Euler angles:(4)
  $R = R_{z} (ϕ) R_{y} (θ) R_{x} (ψ)$
  (5)
  $R = [\begin{matrix} cos θ cos ϕ & sin ψ sin θ cos ϕ - cos ψ sin ϕ & cos ψ sin θ cos ϕ + sin ψ sin ϕ \\ cos θ sin ϕ & sin ψ sin θ sin ϕ + cos ψ cos ϕ & cos ψ sin θ sin ϕ - sin ψ cos ϕ \\ - sin θ & sin ψ cos θ & cos ψ cos θ \end{matrix}]$
The option SIGEQ is available with Eltyp = SHELL or BRICK only (/ANIM/SHELL/SIGEQ). Each material law, in Radioss has its own yield criterion to calculate the equivalent stress. For some it is von Mises; for others, it is Hill or Barlat or something else. For any non-von Mises criterion, the corresponded equivalent stress (or criterion) is computed within all the integration points of the element. Therefore, the output field /ANIM/BRICK/SIGEQ is computed as a mean value over the all integration points.
For brick elements, the stresses are output in the elemental (corotational) coordinate system for /PROP/SOLID elements, and in the orthotropic material coordinate system for /PROP/SOL_ORTH elements, when corotational formulation is used. For all other cases the stresses are output in the global coordinate system.
The option /ANIM/ELEM/SIGX is only applied for shell elements. For brick elements, /ANIM/BRICK/TENS must be used.
User variables are only available for shell and brick elements. When an integration point is not explicitly described, returned integration point means the integration point is superior; computed as [(number of integration points in thickness + 1) / 2]. The result is then rounded up to the superior value.
- Example:
  For two integration points in thickness, second integration point from bottom (top of thickness) is returned.
  
  For three integration points in thickness, second integration point from bottom (middle one) is returned.
  
  For four integration points in thickness, third integration point from bottom is returned.
The Schlieren contour value is computed using $ξ = e^{- c \nabla ρ}$ . Radioss outputs $ξ$ using c=1. The constant $c$ can be updated using HyperView by creating a derived result $ξ^{c_{u s e r}}$ .
For solid elements using LAW51, it is possible to display results for a given sub-material number (1 to 4) using:
- /ANIM/ELEM/LAW51/ALL: results for all sub-material, or
  - /ANIM/ELEM/LAW51/1: results for sub-material 1
  - /ANIM/ELEM/LAW51/2: results for sub-material 2
  - /ANIM/ELEM/LAW51/3: results for sub-material 3
  - /ANIM/ELEM/LAW51/4: results for sub-material 4
  In this case, the following options can be displayed per phase:
  - /ANIM/ELEM/P
  - /ANIM/ELEM/DENS
  - /ANIM/ELEM/ENER
  - /ANIM/ELEM/SSP
  - /ANIM/ELEM/EPSP
  - /ANIM/ELEM/TEMP
  - /ANIM/ELEM/VOLU
  - /ANIM/MASS
For quad or brick elements, /ANIM/ELEM/DAM1, /DAM2, and /DAM3 are available for material LAW24. These values are the principal values of the damage (values in the local cracking skew).
Element time step shows in animation only if elementary time step is computed for this element by Radioss. If nodal time step used (/DT/NODA) in computation, then no element time step shows in animation.