/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

Definition


Field	Contents	SI Unit Example
`Eltyp`	Element type. ELEM BRICK Brick elements
`Restype`	AMS Elements using AMS timestep due to /DT/CST_AMS. BFRAC Burn fraction. COLOR 1vfrac1 + 2vfrac2 + .... + last*vfracN (LAW51 and LAW151 only). DAM1 , DAM2 or DAM3 Damage in direction 1, 2 or 3. 9 DAMG Mean damage value over integration points (only for coupled damage models). 11 DENS Density DT Element time step. 10 EINT Internal energy. ENER Specific energy (internal energy divided by the element mass). EPSD Equivalent strain rate in bricks and in shells too (only available in case of strain rate filtering). EPSP Plastic strain $ε_{p}$ . FILL Filling percentage. This option is used only with `Eltyp` =BRICK for /INIBRI/FILL. 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\|. HOURG Hourglass energy per mass unit K Turbulent energy for ALE material law. LAW51 Display results for all Sub-material or specific Sub-material. 8 MACH Mach number (/MAT/LAW151 (MULTIFLUID) only). 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\|. NL_EPSD Non-local plastic strain rate ${\dot{ε}}_{p}^{n l}$ (only if /NONLOCAL/MAT is activated). 12 NL_EPSP Non-local plastic strain $ε_{p}^{n l}$ (only if /NONLOCAL/MAT is activated). 12 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. 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 P Pressure QVIS Displays artificial viscous pressure SCHLIEREN Schlieren image (optical method widely used in CFD field). Only available with ALE material laws. 7 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 SSP Sound speed. Only available with ALE material laws. TDET Detonation times for high explosive JWL EOS. TEMP Temperature TILLOTSON Region identifier for Tillotson Equation of State (/EOS/TILLOTSON). TSAIWU Tsai Wu criterion for /MAT/LAW12 (3D_COMP), /MAT/LAW14 (COMPSO) and /MAT/LAW25 (COMPSH). TVIS Turbulent viscosity for ALE material law. 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. VFRAC All volumetric fractions (for ALE multi-material laws: LAW20, LAW37, LAW51 and LAW151). VONM von Mises stress 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. VOLU Volume of each phase. 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 if 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:
  $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:
  $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:
  $R_{z} (ϕ) = [\begin{matrix} cos ϕ & - sin ϕ & 0 \\ sin ϕ & cos ϕ & 0 \\ 0 & 0 & 1 \end{matrix}]$
- The angles $ψ, θ, and ϕ$ are the Euler angles:
  $R = R_{z} (ϕ) R_{y} (θ) R_{x} (ψ)$
  $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 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 is displayed in the animation only if elementary time step is computed for this element by Radioss. If nodal time step used (/DT/NODA) in the calculation, element time step is not displayed in the animation.
Option DAMG is only used with coupled damage models (/MAT/LAW72 or /FAIL/GURSON) to output damage over integration points. The damage variable is normalized by its critical value.
- For /MAT/LAW72
  $D_{m g} = \frac{D}{D_{C}}$
- For /FAIL/GURSON
  $D_{m g} = \frac{f_{t}}{f_{F}}$
If /NONLOCAL/MAT option is activated, it is possible to output the regularized non-local plastic strain and its rate.