/MAT/LAW38 (VISC_TAB)
Block Format Keyword This law describes the viscoelastic foam tabulated material and can only be used with solid elements.
Format
(1)  (2)  (3)  (4)  (5)  (6)  (7)  (8)  (9)  (10) 

/MAT/LAW38/mat_ID/unit_ID or /MAT/VISC_TAB/mat_ID/unit_ID  
mat_title  
${\rho}_{i}$  
E_{0}  ${\nu}_{t}$  ${\nu}_{c}$  ${R}_{\nu}$  Iflag  Itota  
$\beta $  H  R_{D}  K_{R}  K_{D}  $\theta $  
K_{air}  fct_ID_{p}  Fscale_{P}  
P_{0}  R_{P}  P_{max}  $\mathrm{\Phi}$  
fct_ID_{ul}  Fscale_{unload}  ${\dot{\epsilon}}_{\mathrm{unload}}$  a  b  
N_{funct}  CUToff  I_{insta}  
E_{final}  ${\epsilon}_{\mathit{final}}$  $\lambda $  Visc  Tol  
Fscale_{1}  Fscale_{2}  Fscale_{3}  Fscale_{4}  Fscale_{5}  
${\dot{\epsilon}}_{1}$  ${\dot{\epsilon}}_{2}$  ${\dot{\epsilon}}_{3}$  ${\dot{\epsilon}}_{4}$  ${\dot{\epsilon}}_{5}$  
fct_ID_{1L}  fct_ID_{2L}  fct_ID_{3L}  fct_ID_{4L}  fct_ID_{5L}  
fct_ID_{1ul}  fct_ID_{2ul}  fct_ID_{3ul}  fct_ID_{4ul}  fct_ID_{5ul} 
Definitions
Field  Contents  SI Unit Example 

mat_ID  Material identifier (Integer, maximum 10 digits) 

unit_ID  Unit Identifier (Integer, maximum 10 digits) 

mat_title  Material title (Character, maximum 100 characters) 

${\rho}_{i}$  Initial density (Real) 
$\left[\frac{\text{kg}}{{\text{m}}^{\text{3}}}\right]$ 
E_{0}  Minimum tension modulus, used for
interface and time step computation. (Real) 
$\left[\text{Pa}\right]$ 
${\nu}_{t}$  Maximum Poisson's ratio in
tension. Default = 10^{30} (Real) 

${\nu}_{c}$  Maximum Poisson's ratio in
compression. (Real) 

${R}_{\nu}$  Exponent for Poisson's ratio
computation. (Real) 

Iflag  Analysis formulation type flag. 4
(Integer) 

Itota  Incremental formulation flag.
(Integer) 

$\beta $  Relaxation rate for unloading. Default = 10^{30} (Real) 

H  Hysteresis coefficient for
unloading. Default = 1.0 (Real) 

R_{D}  Damping factor on strain rate. Default = 0.5 (Real) 

K_{R}  Recovery model flag for unloading
(hysteresis loop).
Where, ${E}^{\mathrm{int}}$ and ${E}_{\mathrm{max}}^{\mathrm{int}}$ are current internal energy and maximum internal energy, respectively. 6 (Integer) 

K_{D}  Decay model flag, hysteresis type.
(Integer) 

$\theta $  Integration coefficient for
instantaneous module update. Default = 0.67 (Real) 

K_{air}  Air content computation flag. 7
(Integer) 

fct_ID_{p}  Pressure curve identifier (pressure
versus relative volume). (Integer) 

Fscale_{P}  Pressure curve scale
factor. (Real) 
$\left[\text{Pa}\right]$ 
P_{0}  Atmospheric
pressure. (Real) 
$\left[\text{Pa}\right]$ 
R_{P}  Relaxation rate of pressure. Default = 10^{30} (Real) 

P_{max}  Maximum air pressure. Default = 10^{30} (Real) 
$\left[\text{Pa}\right]$ 
$\mathrm{\Phi}$  Porosity (density of foam/density of
polymer). (Real) 

fct_ID_{ul}  Unloading function identifier.
(Integer) 

Fscale_{unload}  Unloading function scale
factor. Default = 1.0 (Real) 
$\left[\text{Pa}\right]$ 
${\dot{\epsilon}}_{\mathit{unload}}$  Unloading strain rate (must be greater
than
${\dot{\epsilon}}_{1}$
). (Real) 
$\left[\frac{\text{1}}{\text{s}}\right]$ 
a  Exponent for stress
interpolation. Default = 1.0 (Real) 

b  Exponent for stress
interpolation. Default = 1.0 (Real) 

N_{funct}  Number of functions defining rate
dependency (five or less). (Integer) 

CUToff  Tension cutoff stress. The element is deleted when one element integration point exceeds the tension cutoff stress value. Default = 10^{30} (Real) 
$\left[\text{Pa}\right]$ 
I_{insta}  Material instability control flag.
(Integer) 

E_{final}  Maximum tension modulus. Default = E_{0} (Real) 
$\left[\text{Pa}\right]$ 
${\epsilon}_{\mathit{final}}$  Absolute value of strain at final
modulus. Default = 1.0 (Real) 

$\lambda $  Modulus interpolation
coefficient. Default = 1.0 (Real) 

Visc  Maximum viscosity. 10 Default = 10^{30} (Real) 
$\left[\text{Pa}\cdot \text{s}\right]$ 
Tol  Tolerance on principal direction
update. Default = 1.0 (Real) 

Fscale_{i}  Scale factor for curve
i. (Real) 
$\left[\text{Pa}\right]$ 
${\dot{\epsilon}}_{i}$  Engineering strain rate for curve
i. (Real) 
$\left[\frac{\text{1}}{\text{s}}\right]$ 
fct_ID_{iL}  Loading function identifier for curve
i. (Integer) 

fct_ID_{iul}  Unloading function identifier for curve
i. (Integer) 
Example (Foam)
#RADIOSS STARTER
#12345678910
/UNIT/1
unit for mat
Mg mm s
#12345678910
# 2. MATERIALS:
#12345678910
/MAT/VISC_TAB/1/1
Foam
# RHO_I
2E10
# E_0 nu_t nu_c R_V Iflag Itota
200 0 0 0 0 0
# Beta H R_D K_R K_D Teta
0 0 0 0 0 0
# K_air fct_ID_p Fscale_P
0 0 1
# P0 Rp Pmax Phi
0 0 0 0
#funID_unl Fscale_unload Eps_._unload a b
0 0 0 0 0
# N_funct CUT_off I_insta
1 0 0
# E_final Eps_final Lambda Visc Tol
0 0 0 0 0
# Fscale_i
1
# Eps_._i
0
# func_ID_iload
4
# func_ID_iunload
0
#12345678910
# 3. FUNCTIONS:
#12345678910
/FUNCT/4
function_4
# X Y
1 200
1 200
#12345678910
#ENDDATA
/END
#12345678910
Comments
 Engineering stress versus engineering strain can be input as functions for different strain rates. The stress and strain is positive in compression and negative in tension. By default (Itota=0), the tension behavior is linear elastic using Young’s modulus, E_{0}. If Itota=1, the engineering stress strain behavior should be input using the functions, fct_ID_{iL}, with the stress strain curve defined both in compression and tension
 When stress strain functions are
defined at different strain rates, the stresses are computed by interpolation from input
functions:
$\sigma =f(\epsilon ,\dot{\epsilon})$
for given $\dot{\epsilon}$ , read two values of function at $\epsilon $ for the two immediately lower and higher strain rates.with
(1) $$\sigma ={\sigma}_{2}+\left({\sigma}_{1}{\sigma}_{2}\right){\left[1{\left(\frac{\dot{\epsilon}{\dot{\epsilon}}_{1}}{{\dot{\epsilon}}_{2}{\dot{\epsilon}}_{1}}\right)}^{a}\right]}^{b}$$The parameters $a$ and $b$ define the shape of the interpolation function within each interval. If $a$ = $b$ = 1, the interpolation is linear.
 A "coupled" set of principal
nominal stresses is computed with anisotropic Poisson's ratios:
${\nu}_{ij}={\nu}_{c}+\left({\nu}_{t}{\nu}_{c}\right)\left(1\mathrm{exp}\left({R}_{v}\left{\epsilon}_{ij}\right\right)\right)$ in tension ( ${\epsilon}_{ij}\ge 0$ )
${\nu}_{ij}={\nu}_{c}$ in compression.
Where, ${\epsilon}_{ij}=\frac{({\epsilon}_{i}+{\epsilon}_{j})}{2}$
 ${\epsilon}_{ij}\ge 0$
 Analysis formulation type Iflag.
Iflag=0: corresponds to the viscoelastic foam tabulated material (viscoelasticity is computed in each principal stress direction).
Iflag=1: behavior will be linear in both tension and compression, following Hook's relations.
For compression, Young's modulus ${E}_{0}$ and Poisson's ratio ${\nu}_{c}$ are used.
In tension, the instantaneous Young's modulus ratio ${E}_{t}$ is used.
The instantaneous Young’s modulus is updated using:(2) $${E}_{t}={E}_{final}+({E}_{0}{E}_{final})\left[1{e}^{\lambda ({V}_{R}1+{\epsilon}_{final})}\right]$$with(3) $${E}_{0}<E<{E}_{f}{}_{inal}$$Where, ${E}_{0}$
 Minimum tension modulus
 ${E}_{final}$
 Maximum tension modulus
 ${V}_{R}$
 Relative volume computed in Radioss
 ${\epsilon}_{final}$
 Absolute value of the strain corresponding to the maximum compression modulus.
The instantaneous modulus is only used for tension.
 For stability,
$\dot{\epsilon}$
is filtered using:
(4) $${\dot{\epsilon}}_{\mathit{filt}}^{n}={\dot{\epsilon}}_{\mathit{filt}}^{n1}+{R}_{D}\left({\dot{\epsilon}}^{n}{\dot{\epsilon}}_{\mathit{filt}}^{n1}\right)$$  Hysteresis is applied in linear
tension case.
If K_{R}=1, Hysteresis is only applied in compression.
If K_{R}=2, Hysteresis is applied both in compression and in tension.
 For air pressure
${P}_{air}$
(when
K_{air}=1)If fct_ID_{p}≠0:
(5) $${P}_{\mathit{air}}={\mathit{Fscale}}_{p}.\mathrm{f}\left(\frac{V}{{V}_{0}}\right)$$Where, $f$ refers to function number fct_ID_{p}.
If fct_ID_{p}=0:(6) $${P}_{\mathit{air}}={P}_{0}\frac{\left(1\frac{V}{{V}_{0}}\right)}{\left(\frac{V}{{V}_{0}}\mathrm{\Phi}\right)}$$Relaxation is applied as:(7) $${P}_{\mathit{air}}=\mathrm{min}\left({P}_{\mathit{air}},{P}_{\mathrm{max}}\right)\mathrm{exp}\left({R}_{p}t\right)$$Where, ${R}_{p}$ is the relaxation rate of pressure and t is the time.
 During unloading, without an
unloading curve defined fct_ID_{iul} =
fct_ID_{ul}=0,
$\sigma $
is computed from the first loading curve,
fct_ID_{1L}.
If the unloading curve is defined, $\sigma $ is interpolated between the first loading curve fct_ID_{1L} and the defined unloading curve fct_ID_{ul} or fct_ID_{iul}. In this case, fct_ID_{1L} must correspond to a quasistatic state.
 Unloading functions fct_ID_{iul} (Line 12) are used only if the unloading curve fct_ID_{ul} is not defined.
 If Visc
is input, interpolated stress will be limited by this value to have a larger
timestep:
(8) $$\sigma \le {\sigma}_{1}+\mathit{Visc}\left(\dot{\epsilon}{\dot{\epsilon}}_{1}\right)$$  The behavior is strain rate independent when $\dot{\epsilon}\le {\dot{\epsilon}}_{1}$ .
 /VISC/PRONY can be used with this material law to include viscous effects.