/MAT/LAW64 (UGINE_ALZ)

Block Format Keyword This law describes the Ugine & Alz trip steel material. This material law can be used only with shell elements.

Format


(1)	(2)	(3)	(5)	(7)	(9)
/MAT/LAW64/`mat_ID`/`unit_ID` or /MAT/UGINE_ALZ/`mat_ID`/`unit_ID`
`mat_title`
$ρ_{i}$
`E`		$ν$	`C`_p
`D`		`n`	`M`_d	`V`₀	`V`_m
`fct_ID`₀	`fct_ID`₁	`Fscale`₀	`Fscale`₁	`T`₀

Definition


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)
$ρ_{i}$	Initial density. (Real)	$[\frac{kg}{m^{3}}]$
`E`	Initial Young's modulus. (Real)	$[Pa]$
$ν$	Poisson's ratio. (Real)
`C`_p	Specific heat capacity. Default = 10³⁰ (Real)	$[\frac{J}{kg \cdot K}]$
`D`	Material parameter 1. (Real)
`n`	Material parameter 2. (Real)
`M`_d	Limit martensite transformation temperature. (Real)	$[K]$
`V`₀	Material parameter. (Real)
`V`_m	Constant martensite fraction for second yield stress function 0 < `V`_m ≤ 1. (Real)
`fct_ID`₀	Yield stress function identifier for 0 martensite fraction. (Integer)
`fct_ID`₁	Yield stress function identifier for `V`_m martensite fraction. (Integer)
`Fscale`₀	Scale factor for yield function for `fct_ID`₀. (Real)	$[Pa]$
`Fscale`₁	Scale factor for yield function for `fct_ID`₁. (Real)	$[Pa]$
`T`₀	Initial temperature. (Real)	$[K]$

Example (Steel)

#RADIOSS STARTER
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/UNIT/1
unit for mat
                  Mg                  mm                   s
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#-  2. MATERIALS:
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/MAT/LAW64/1/1
Steel
#              RHO_I
              7.8E-9                   
#                  E                  Nu                  Cp
              210000                  .3           460000000
#                  D                   n                  Md                  V0                  Vm
                   4                 3.5                 356                  .2                  .6
# func_ID0  func_ID1             Fscale0             Fscale1                  T0
         1         2                   1                   1                 323
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#-  3. FUNCTIONS:
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/FUNCT/1
function_1
#                  X                   Y
                   0                 250                                                            
                .001                 350                                                            
                  .5                1100                                                            
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/FUNCT/2
function_2
#                  X                   Y
                   0                 930                                                            
                .001                1000                                                            
                  .5                1500                                                            
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#ENDDATA
/END
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|

Comments

Martensite fraction:
$V_{m} (ε_{p}, T) = V_{m}^{max} (T) \cdot (1 - e^{- (D \cdot ε_{p}) n})$
$V_{m}^{max} (T) = V_{0} Ln (M_{d} - T + 1)$
if
$T < M_{d}$
$V_{m}^{max} (T) = 0$
if
$T > M_{d}$
Mechanical behavior:
The yield plastic stress is computed by linear interpolation between two curves fct_ID₁ and fct_ID₀.
The temperature is computed assuming the adiabatic condition (by default the condition is isothermal with C_p = 10³⁰):
$T = T_{0} + \frac{E_{int}}{ρ C_{p} (Volume)}$
Where, E_int is the internal energy of the element.