Products: ABAQUS/Standard ABAQUS/Explicit
ABAQUS allows the introduction of the inelastic heat fraction, , which defines heat generation caused by mechanical dissipation associated with plastic straining. This term can be introduced as a source of coupling for thermal-mechanical analysis. Such coupling might be important in a simulation in which extensive inelastic deformation is occurring fairly rapidly in a material whose mechanical properties are temperature dependent. If the process is very slow, the heat generated by the plastic deformation has time to dissipate; and uncoupled, isothermal, analysis is sufficient to model the process. If the process is extremely rapid, the heat has no time to diffuse; and uncoupled, adiabatic, analysis (in which each integration point is treated as if it is thermally insulated from its neighbors) is sufficient. Fully coupled thermal-stress analysis is required for cases that lie far enough from both extremes. This section defines the heat generation term caused by inelastic straining and describes how this term contributes to the overall Newton solution scheme.
The model assumes that plastic straining gives rise to a heat flux per unit volume of
ABAQUS generally uses a backward Euler scheme to integrate the plastic strain, so at the end of the increment is approximated as
When Newton's method is used to solve the nonlinear equations, the coupling term gives rise to three contributions to the Jacobian matrix for the Newton method:
The mechanical constitutive model has the following general form. The elasticity defines the stresses by
where is the strain energy density potential and is the mechanical elastic strain. We implicitly assume that the elasticity is not fully incompressible, although the derivation is not significantly different if this is not the case, since the pressure stress will do no work in a fully incompressible material and so makes no contribution to the terms under discussion.We assume that there is an additive strain rate decomposition that can be integrated to give
where is the total strain and is the strain caused by thermal expansion. In the constitutive models in ABAQUS only. This form of decomposition of the deformation depends on being measured as the integrated rate of deformation and on the elastic and thermal strains being small: this is true for the standard plasticity models provided in the program.The plastic flow definition is integrated by the backward Euler method to give
Finally, assuming there is a single active yield surface or a single active flow surface, rate-independent models introduce a yield surface constraint, while rate-dependent models provide an integrated flow rate constraint, both of which are incorporated in the general form
where is a scalar stress function (for example, the Mises or Hill stress function for simple metal plasticity models) and is the yield stress for a rate-independent model, while for a rate-dependent model, whereEquation 4.3.10–1 to Equation 4.3.10–4 are a general definition of all of the standard isotropic hardening plasticity models integrated by the backward Euler method.
We now take variations of these equations with respect to all quantities at the end of the increment:
For simplicity of notation we now define
These expressions allow us to write