Iform = 0

Block Format Keyword

/MAT/LAW51 - Iform = 0: Multi-material Law with up to 3 Elasto Materials (Solid, Liquid, or Gas)

Description

Law 51 is able to handle up to three elastic materials (solid, liquid, or gas). The material boundaries inside an element are not explicitly defined, but an anti-diffusive technique is used to avoid expansion of transition zone (Comment 1).

Compatible only with 3D analysis and Euler or ALE formulation.

LAW51 is based on equilibrium between each material present inside the element. RADIOSS computes and outputs a relative pressure . At each cycle:

= 1 = 2 = 3

Total pressure can be calculated with external pressure:

P = + Pext

Where, P is positive for a compression and negative for traction.

Hydrostatic stresses are computed from Polynomial EOS:

dEint = -( + Pext ) dV

Where, which means that the EOS is linear for an expansion and cubic for a compression.

Deviatoric stresses are computed with shear modulus:

Format

(1)	(2)	(3)	(4)	(5)	(6)	(7)	(8)	(9)	(10)
/MAT/LAW51/mat_ID
mat_title
Blank
Iform

#Global Parameters

(1)	(2)	(3)	(4)	(5)	(6)	(7)	(8)	(9)	(10)
Pext

#Material1 Parameters

(1)	(3)	(5)	(7)	(9)
		E0mat_1		C0mat_1
C1mat_1	C2mat_1	C3mat_1	C4mat_1	C5mat_1
G1mat_1

#Material2 Parameters

(1)	(3)	(5)	(7)	(9)
		E0mat_2		C0mat_2
C1mat_2	C2mat_2	C3mat_2	C4mat_2	C5mat_2
G1mat_2

#Material3 Parameters

(1)	(3)	(5)	(7)	(9)
		E0mat_3		C0mat_3
C1mat_3	C2mat_3	C3mat_3	C4mat_3	C5mat_3
G1mat_3

Field	Contents	SI Unit Example
mat_ID	Material identifier (Integer, maximum 10 digits)
mat_title	Material title (Character, maximum 100 characters)
Iform	Formulation flag (Integer)
Pext	External pressure (Comment 2) Default = 0 (Real)
	Kinematic viscosity shear (Comment 3) Default = 0 (Real)
	Kinematic viscosity (volumetric), which corresponds to Stokes Hypothesis (Comment 3) Default = 0 (Real)
	Initial volumetric fraction (Comment 4) (Real)
	Initial density (Real)
E0mat_i	Initial energy per unit volume (Real)
	Hydrodynamic cavitation pressure (Comment 5) If fluid material (), then default = -Pext If solid material (), then default = -1e30 (Real)
C0mat_i	Initial pressure (Real)
C1mat_i	Hydrodynamic coefficient (Real)
C2mat_i	Hydrodynamic coefficient (Real)
C3mat_i	Hydrodynamic coefficient (Real)
C4mat_i	Hydrodynamic coefficient (Real)
C5mat_i	Hydrodynamic coefficient (Real)
G1mat_i	Elasticity shear modulus Default = 0 (Fluid material) (Real)

#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|

/MAT/LAW51/97

100%Water-MULTIMAT:AIR+WATER(total pressure formulation:Pext=0),units {kg,m,s,Pa}

#--------------------------------------------------------------------------------------------------#

# Material Law No 51. MULTI-MATERIAL SOLID LIQUID GAS -ALE-CFD-SPH

#--------------------------------------------------------------------------------------------------#

# Blank format

# IFLG

#---Global parameters------------------------------------------------------------------------------#

# P_EXT NU LAMDA

0 0 0

#---Material#1:AIR(PerfectGas)---------------------------------------------------------------------#

# ALPHA_1 RHO_0_1 E_0_1 P_MIN_1 C_0_1

0.0 1.2 2.5E+05 -1e30 0

# C_1_1 C_2_1 C_3_1 C_4_1 C_5_1

0 0 0 0.4 0.4

# G_1

#---Material#2:WATER(Linear_Incompressible)--------------------------------------------------------#

# ALPHA_2 RHO_0_2 E_0_2 P_MIN_2 C_0_2

1.0 1000.0 2.5E+05 -1e30

# C_1_2 C_2_2 C_3_2 C_4_2 C_5_2

0 0 0 0 0

# G_2

#---Material#3:not defined-------------------------------------------------------------------------#

# ALPHA_3 RHO_0_3 E_0_3 P_MIN_3 C_0_3

0.0 0 0 0 0

# C_1_3 C_2_3 C_3_3 C_4_3 C_5_3

0 0 0 0 0

# G_3

#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|

/MAT/LAW51/98

100%Water-MULTIMAT:AIR+WATER(relative pressure formulation:Pext>0),units {kg,m,s,Pa}

#--------------------------------------------------------------------------------------------------#

# Material Law No 51. MULTI-MATERIAL SOLID LIQUID GAS -ALE-CFD-SPH

#--------------------------------------------------------------------------------------------------#

# Blank format

# IFLG

#---Global parameters------------------------------------------------------------------------------#

# P_EXT NU LAMDA

1E+5 0 0

#---Material#1:AIR(PerfectGas)---------------------------------------------------------------------#

# ALPHA_1 RHO_0_1 E_0_1 P_MIN_1 C_0_1

0.0 1.2 2.5E+05 -1E+5 -1E+5

# C_1_1 C_2_1 C_3_1 C_4_1 C_5_1

0 0 0 0.4 0.4

# G_1

#---Material#2:WATER(Linear_Incompressible)--------------------------------------------------------#

# ALPHA_2 RHO_0_2 E_0_2 P_MIN_2 C_0_2

1.0 1000.0 0 -1E+5 0

# C_1_2 C_2_2 C_3_2 C_4_2 C_5_2

0 0 0 0 0

# G_2

#---Material#3:not defined-------------------------------------------------------------------------#

# ALPHA_3 RHO_0_3 E_0_3 P_MIN_3 C_0_3

0.0 0 0 0 0

# C_1_3 C_2_3 C_3_3 C_4_3 C_5_3

0 0 0 0 0

# G_3

#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|

1.	The anti-diffusive technique can be adjusted with /UPWIND from the Starter Input. The 3 flag is the upwind coefficient for a damp area:

3 = 1: Full Upwind (default, recommended value)

3 = 1e-30: Zero Upwind (less diffusive, but potentially unstable)

3 = -1: Full Downwind (Anti-Diffusive technique, potentially unstable)

2.	RADIOSS computes and outputs a relative pressure .

However, total pressure is essential for energy integration (dEint = -PdV). It can be computed with the external pressure flag Pext.

P = + Pext leads to dEint = -(Pext + ) dV.

This means that if Pext = 0, the computed pressure is also the total pressure: = P

3.	Kinematic viscosities are global and is not specific to each material. It allows computing viscous stress tensor:

Where, is the cinematic shear viscosity flag, and

is the cinematic volumetric viscosity flag.

4.	Volumetric fractions enable the sharing of elementary volume within the three different materials.

For each material must be defined between 0 and 1.

Sum of initial volumetric fractions must be equal to 1.

For automatic initial fractioning of the volume, refer to the /INIVOL card.

5.	flag is the minimum value for the computed pressure . It means that total pressure is also bounded to:

For fluid materials and detonation products must remain positive to avoid any tensile strength so must be set to -Pext.

For solid materials, default value = 1e-30 is suitable but may be modified.

6.	Material tracking is possible through animation files:

/ANIM/BRIC/VFRAC (All material volumetric fractions)

7.	The following global outputs are available for animation files:

/ANIM/BRICK/PLA51 (global plasticity)

/ANIM/BRICK/TEM51 (global temperature)