6.5.5 Adiabatic analysis

**Products: **ABAQUS/Standard ABAQUS/Explicit ABAQUS/CAE

An adiabatic stress analysis:

is used in cases where mechanical deformation causes heating but the event is so rapid that this heat has no time to diffuse through the material—for example, a very high-speed forming process;

can be conducted as part of a dynamic analysis (“Implicit dynamic analysis using direct integration,” Section 6.3.2, or “Explicit dynamic analysis,” Section 6.3.3) or as part of a static analysis (“Static stress analysis,” Section 6.2.2);

is available only for the isotropic hardening metal plasticity models with a Mises yield surface (“Classical metal plasticity,” Section 18.2.1);

can be conducted if parts of the model are elastic only—no change in temperature occurs in the elastic regions; and

requires that a material's density, specific heat, and inelastic heat fraction (fraction of inelastic dissipation rate that appears as heat flux) be specified.

Adiabatic analysis

Adiabatic thermal-stress analysis is typically used to simulate high-speed manufacturing processes involving large amounts of inelastic strain, where the heating of the material caused by its deformation is an important effect because of temperature-dependent material properties. The temperature increase is calculated directly at the material integration points according to the adiabatic thermal energy increases caused by inelastic deformation; temperature is not a degree of freedom in the problem. No allowance is made for conduction of heat in an adiabatic analysis. For problems where both inelastic heating and conduction of the heat are important, a fully coupled temperature-displacement analysis must be performed (“Fully coupled thermal-stress analysis,” Section 6.5.4).

In an adiabatic analysis plastic straining gives rise to a heat flux per unit volume of

where is the heat flux that is added into the thermal energy balance, is the user-specified inelastic heat fraction (assumed constant; discussed below), is the stress, and is the rate of plastic straining. The heat equation solved at each integration point is

where is the material density and is the specific heat (see “Density,” Section 16.2.1, and “Specific heat,” Section 20.2.3).

ABAQUS/CAE Usage: | Use any of the following procedures to perform an adiabatic analysis: |

Step module: |

Subsequent thermal diffusion analysis in ABAQUS/Standard

In ABAQUS/Standard thermal diffusion analysis can be performed after the adiabatic calculation (for example, to study the cool-down of a component after sudden deformation). In this case the temperatures at the end of the adiabatic analysis must be written to the ABAQUS/Standard results file as element variables averaged at the nodes. Since temperature values in an adiabatic analysis can be written to the results file as element quantities only by using the TEMP output variable identifier, they cannot be read directly into a subsequent thermal diffusion analysis as initial conditions. However, if you postprocess the results file to produce a second results file in which the temperature data are provided as nodal quantities, a subsequent heat transfer analysis can be performed with these temperatures as initial conditions. See “Predefined fields,” Section 27.6.1, and “Accessing the results file information,” Section 5.1.3, for details. Alternatively, you could postprocess the results file to produce a data list containing data pairs consisting of nodes and temperatures.

The temperatures, NT, obtained from the heat transfer analysis can then be used to drive a continuation of the previous stress analysis. This stress analysis should be restarted from the end of the adiabatic analysis and will provide the response to the change of the temperature field obtained during the heat transfer analysis. In this case ABAQUS/Standard will automatically read the temperatures from the results file that was obtained from the heat transfer analysis and apply them in the restarted analysis.

The following input options could be used to perform a heat transfer analysis using the temperatures from an adiabatic analysis and then continue the stress analysis:

**Static adiabatic analysis … *STEP *STATIC, ADIABATIC … **Write the temperatures to the results file as element **variables averaged at the nodes *EL FILE, POSITION=AVERAGED AT NODES TEMP *END STEP **Heat transfer analysis using the temperatures from the **static analysis as initial conditions … *INITIAL CONDITIONS, TYPE=TEMPERATURE, FILE=new results file, STEP=step, INC=increment*STEP *HEAT TRANSFER … *NODE FILE NT *END STEP **Restart from the adiabatic analysis using temperatures **obtained from the heat transfer analysis *RESTART, WRITE, READ, STEP=k, INC=i, END STEP … *STEP *STATIC … *TEMPERATURE, FILE=heat_transfer_results_file… *END STEP

If the continuation of the analysis into thermal diffusion requires a fully coupled temperature-displacement analysis (see “Fully coupled thermal-stress analysis,” Section 6.5.4), the simplest (but more expensive) approach is to use coupled temperature-displacement elements throughout the adiabatic analysis. At the end of the static or the dynamic adiabatic calculations, the temperatures must be written to the results file as element variables averaged at the nodes. In addition, you must constrain all temperature degrees of freedom since they are not used in the adiabatic analysis. The adiabatic analysis can then be restarted to apply the correct temperature distribution obtained from the adiabatic analysis to the temperature degree of freedom of each node in the model. To create the input for the boundary conditions, you must postprocess the results file obtained from the adiabatic analysis and extract the value of TEMP at each node in the model (see “Accessing the results file information,” Section 5.1.3). The temperature boundary conditions can be released as needed in subsequent coupled temperature-displacement analysis steps.

The following input options could be used to perform a coupled temperature-displacement analysis using the temperatures from an adiabatic analysis:

**Static adiabatic analysis, coupled temperature-displacement **plane stress elements … *ELEMENT, TYPE=CPS4T, ELSET=EALL … *BOUNDARYnodes, 11, 11, 0.0 *STEP *STATIC, ADIABATIC … **Write the temperatures to the results file as element **variables averaged at the nodes *EL FILE, POSITION=AVERAGED AT NODES TEMP *END STEP **Restart from the adiabatic analysis *RESTART, WRITE, READ, STEP=k, INC=i, END STEP … *STEP *STATIC **Dummy step to associate the temperature variable TEMP with **the temperature degree of freedom at each node 1.0, 1.0 … *BOUNDARY, OP=NEWnode, 11, 11,temperature… *END STEP **Coupled temperature displacement run for cool down of **structure: continuation of the restart analysis … *STEP *COUPLED TEMPERATURE-DISPLACEMENT 0.1, 1.0 … *BOUNDARY, OP=NEW **no temperature boundary condition specified *END STEP

Initial temperatures can be prescribed at nodes as initial conditions. Initial values of stresses, field variables, solution-dependent state variables, etc. can also be specified (see “Initial conditions,” Section 27.2.1).

Boundary conditions can be applied to displacement degrees of freedom in an adiabatic analysis in the same way that they are applied in nonadiabatic dynamic, explicit dynamic, or static analysis steps (see “Boundary conditions,” Section 27.3.1). Temperature is not a degree of freedom in an adiabatic analysis.

The loading options available for an adiabatic analysis are the same as those available for nonadiabatic dynamic, explicit dynamic, or static analysis steps (see “Applying loads: overview,” Section 27.4.1).

The following types of mechanical loads can be prescribed:

Concentrated nodal forces can be applied to the displacement degrees of freedom (1–6); see “Concentrated loads,” Section 27.4.2.

Distributed pressure forces or body forces can be applied; see “Distributed loads,” Section 27.4.3. The distributed load types available with particular elements are described in Part VI, “Elements.”

Predefined temperature fields cannot be used during an adiabatic analysis step.

The values of user-defined field variables can be specified; these values affect only field-variable-dependent material properties, if any. See “Predefined fields,” Section 27.6.1.

Only Mises plasticity with isotropic elasticity and isotropic hardening is allowed in adiabatic stress analysis (“Inelastic behavior,” Section 18.1.1). Kinematic or combined hardening is not available, but rate effects can be included. However, portions of the model can include only elastic material; no change in temperature occurs in the elastic regions, since there is no source of heat generation.

You must specify the density, the inelastic heat fraction, and the specific heat as part of the material definition for the material in which heat will be generated by plastic dissipation. You can also specify latent heat if necessary (“Latent heat,” Section 20.2.4).

The inelastic heat fraction is the amount of inelastic dissipation used to calculate the increase in temperature. The default value of the inelastic heat fraction is 0.9. If the inelastic heat fraction is not included in the material definition, the heat generated by inelastic deformation is not included in the analysis.

In ABAQUS/Standard adiabatic analyses can also be carried out with user subroutine `UMAT`. In this case the temperature must be defined as a solution-dependent state variable, and all coupling terms must be included in the user subroutine. If conductivity (“Conductivity,” Section 20.2.2) is defined for the material, it will be ignored during adiabatic analysis steps.

Input File Usage: | All of the following options must be included in the material definition: |

*DENSITY *INELASTIC HEAT FRACTION *SPECIFIC HEAT The following option can be included if latent heat effects are important: *LATENT HEAT |

ABAQUS/CAE Usage: | All of the following must be included in the material definition: |

Property module: The following can be included if latent heat effects are important: Property module: material editor: |

Material properties can be temperature dependent. Since the only source of temperature change in adiabatic analysis is inelastic deformation, the temperature can only rise. This temperature rise may cause thermal expansion (usually a small effect) and localization of the deformation if the flow stress is reduced by the temperature rise. Since the adiabatic assumption applies only in rapid events and inelastic deformation usually causes significant temperature rises only if the deformation is substantial, the strain rates are often large in adiabatic analysis. The softening of the material caused by the temperature rise may, thus, be offset somewhat by strengthening associated with rate dependence if the material is rate sensitive.

Any of the stress/displacement or coupled temperature-displacement elements in ABAQUS can be used in an adiabatic analysis (see “Choosing the appropriate element for an analysis type,” Section 21.1.3). Mass or spring elements will not contribute to the heating of the material since they cannot generate plastic strains.

If coupled temperature-displacement elements are used in an adiabatic analysis, the temperature degrees of freedom will be ignored.

Since temperatures are updated at the material calculation points, output of temperature is available with output variable TEMP, not with output variable NT.

The element output available for an adiabatic analysis includes stress; strain; energies; the values of state, field, and user-defined variables; and composite failure measures. The nodal output available includes displacements, reaction forces, and coordinates. All of the output variable identifiers are outlined in “ABAQUS/Standard output variable identifiers,” Section 4.2.1, and “ABAQUS/Explicit output variable identifiers,” Section 4.2.2.

*HEADING … *MATERIAL, NAME=name*ELASTIC, TYPE=ISOTROPICData lines to define isotropic linear elasticity*PLASTICData lines to define metal plasticity*DENSITYData lines to define density*INELASTIC HEAT FRACTIONData line to define inelastic heat fraction*SPECIFIC HEATData lines to define specific heat… *BOUNDARYData lines to specify zero-valued boundary conditions*INITIAL CONDITIONS, TYPE=typeData lines to specify initial conditions*AMPLITUDE, NAME=nameData lines to define amplitude variations** *STEP, NLGEOMThe NLGEOM parameter is used in ABAQUS/Standard to include geometric nonlinearity*DYNAMIC, ADIABATICor*DYNAMIC, EXPLICIT, ADIABATICor*STATIC, ADIABATICData line to control time incrementation or to specify the time period of the step*BOUNDARY, AMPLITUDE=nameData lines to describe nonzero or zero-valued boundary conditions*CLOAD and/or *DLOAD and/or *DSLOADData lines to specify loads*FIELDData lines to specify field variable values*END STEP