Elastic-visco-plastic constitutive law with thermal effects for solid elements at variable temperature (to check before use, A-M. HABRAKEN, june 91).
Coupled thermo-mechanical analysis of elastic-visco-plastic solids undergoing large strains
Prepro: LIRSTH.F
Lagamine: IRSC2E.F, IRSC2A.F, IRSC2G.F
Plane stress state | NO |
Plane strain state | YES |
Axisymmetric state | YES |
3D state | NO |
Generalized plane state | YES |
Line 1 (2I5, 60A1) | |
---|---|
IL | Law number |
ITYPE | 240 |
COMMENT | Any comment (up to 60 characters) that will be reproduced on the output listing. |
Line 1 (3I5) | |
---|---|
NTEMP | number of temperature at which material data are given |
NINTV | number of sub-steps used to integrate numerically the constitutive equation in a time step |
IENTH | 0 to use the classical formulation of the heat problem 1 to use the enthalpy formulation of the heat problem |
2 lines repeated NTEMP times
Line 1 (7G10.0) | |
---|---|
T | temperature |
E | YOUNG's elastic modulus at temperature T |
ANU | POISSON's ratio at temperature T |
ALPHA | thermal expansion coefficient ($\alpha$) at temperature T |
AN | strain rate exponent (n) at temperature T |
B | strain rate coefficient (B) at temperature T |
AM | hardening exponent (m) at temperature T |
Line 2 (10X, 4G10.0) | |
H1 | hardening coefficient ($H_1$) at temperature T |
AQ | recovery exponent (q) at temperature T |
H2 | recovery coefficient ($H_2$) at temperature T |
AKO | initial yield limit ($K_o$) at temperature T |
6 for the 3D state
4 for the other cases
The stresses are the components of CAUCHY stress tensor in global (X,Y,Z) coordinates.
For the 3-D state:
SIG(1) | $\sigma_{XX}$ |
SIG(2) | $\sigma_{YY}$ |
SIG(3) | $\sigma_{ZZ}$ |
SIG(4) | $\sigma_{XY}$ |
SIG(5) | $\sigma_{XZ}$ |
SIG(6) | $\sigma_{YZ}$ |
For the other cases:
SIG(1) | $\sigma_{XX}$ |
SIG(2) | $\sigma_{YY}$ |
SIG(3) | $\sigma_{XY}$ |
SIG(4) | $\sigma_{ZZ}$ |
4
Q(1) | = element thickness (t) in plane stress state = 1 in plane strain state = circumferential strain rate $\dot{\varepsilon_\theta}$ in axisymmetric state = 0 in 3D state |
Q(2) | current yield limit in tension, its initial value is $K_o$ |
Q(3) | hydrostatic stress ($\sigma_m$) |
Q(4) | difference between the current temperature and the initial temperature |