8 node large strain volumetric element with hourglass and locking control for thermo-mechanical analysis.
Implemented by: Zhu Yongyi, December 1991
Improved by: Lihong Zhang, May 2005 (BWD3T version)
Type: 222
Prepro: BLZ3TA.F
Lagamine: BLZ3TB.F, BWD3TB.F
TITLE (A5) | |
---|---|
TITLE | 'BLZ3T' in columns 1 to 5 |
CONTROL (5I5) | |
NELEM | Number of elements |
INDPP | = 0 if no weight |
= 1 if weight taken into account | |
INSIG | = 0 no initial stresses |
= 1 initial stresses computed from ferrostatic pressure (see below) | |
INSHE | = 0 for automatic calculation of shear locking parameter |
= 1 if shear coefficient taken into account (see below) | |
=-1 for use of element BWD3T (only 1 integration point) | |
ILOAX | = 0 for global axis computation ☛ Objectivity must be verified in the material law (Jaumann correction) ☛ No rotation of material axes |
< 0 for computation with constant and symetrical velocity gradients pseudo local axes : use of local axes on the time step but no evolution of the local axes on the following time step ☛ Objectivity is verified ☛ No rotation of material axes |
|
> 0 for computation with local axes ☛ Objectivity is verified ☛ Rotation of material axes |
|
units: = 1 for rotations included in local tangent matrix ![]() = 2 apply final rotation to local tangent matrix = 3 apply initial rotation to local tangent matrix = 4 compute tangent matrix through global perturbation method |
|
tens (only for ILOAX>0): = 0 for local axes e1, e2, e3 initially parallel to global axes ex, ey, ez = 1 for local axes e1, e2 given (and e3=e1∧e2) = 2 for local axes e1, e2 initially in the plane (ex, ey) forming an angle θ with ex, ey (and e3=e1∧e2) = 3 same as 1 with different local axes for each element = 4 same as 2 with different local axes for each element |
|
Consideration of weight (5G10.0) Only if INDPP = 1 |
|
WSPE(1) | = specific weight in X direction |
WSPE(2) | = specific weight in Y direction |
WSPE(3) | = specific weight in Z direction |
WSPE(4) | = constant heat source |
WSPE(5) | = density |
Consideration of initial stresses from ferrostatic pressure (3G10.0) Only if INSIG = 1 |
|
GAMMA | = specific weight |
TSOL | = solidus temperature |
TLIQ | = liquidus temperature |
Consideration of shear locking (1G10.0) Only if INSHE = 1 |
|
PARSHE | Shear locking coefficient ∈ [0,1] - close to 0: avoid shear locking but higher risk of hourglass modes (use for thin elements in bending) - close to 1: avoid hourglass modes but higher risk of shear locking (use for cubic elements in shear) |
Initial orientation of the local axes (6G10.0) Only if tens of ILOAX = 1 or 3 (only one line if tens of ILOAX = 1, repeated for each element if tens of ILOAX = 3) |
|
e1(x) | coordinate of e1 along ex |
e1(y) | coordinate of e1 along ey |
e1(z) | coordinate of e1 along ez |
e2(x) | coordinate of e2 along ex |
e2(y) | coordinate of e2 along ey |
e2(z) | coordinate of e2 along ez |
Note: These vectors are normalized after reading but should be orthogonal: e1 • e2 = e1(x) * e2(x) + e1(y) * e2(y) + e1(z) * e2(z) = 0 |
|
Initial orientation of the local axes (1G10.0) Only if tens of ILOAX = 2 or 4 (only one line if tens of ILOAX = 2, repeated for each element if tens of ILOAX = 4) |
|
THETA | Angle between e1 and ex in degrees |
Definition of the elements (4I5/14I5) | |
NINTE | Number of integration points (1, 2, 4 or 8) Currently, only NINTE = 1 is available ! |
LMATE1 | Number of the material mechanical law |
LMATE2 | Number of the material thermal law |
LMATE3 | Number of the material metallurgical law (use 0 if not relevant) |
NODES(8) | List of nodes |
$\sigma_x,\sigma_y,\sigma_z,\sigma_{xy},\sigma_{xz},\sigma_{yz},f_x,f_y,f_z,f_{capacitif}$ In global axes