====== BLZ3T/BWD3T ====== ===== Description ===== 8 node large strain volumetric element with hourglass and locking control for thermo-mechanical analysis. \\ {{ :elements:blz3d.png?300|}} \\ 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