Table of Contents
BLZ3T/BWD3T
Description
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
Input file
| 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  Not available = 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 |
Results
$\sigma_x,\sigma_y,\sigma_z,\sigma_{xy},\sigma_{xz},\sigma_{yz},f_x,f_y,f_z,f_{capacitif}$ In global axes