Products: ABAQUS/Standard ABAQUS/Explicit ABAQUS/CAE
A fully coupled thermal-stress analysis:
is performed when the mechanical and thermal solutions affect each other strongly and, therefore, must be obtained simultaneously;
requires the use of elements with both temperature and displacement degrees of freedom;
can be used to analyze time-dependent material response; and
cannot include cavity radiation effects (Cavity radiation, Section 32.1.1).
neglects inertia effects; and
can be transient or steady-state.
includes inertia effects; and
models transient thermal response.
Fully coupled thermal-stress analysis is needed when the stress analysis is dependent on the temperature distribution and the temperature distribution depends on the stress solution. For example, metalworking problems may include significant heating due to inelastic deformation of the material which, in turn, changes the material properties. In addition, contact conditions exist in some problems where the heat conducted between surfaces may depend strongly on the separation of the surfaces or the pressure transmitted across the surfaces (see Thermal contact properties, Section 30.2.1). For such cases the thermal and mechanical solutions must be obtained simultaneously rather than sequentially. Coupled temperature-displacement elements are provided for this purpose in both ABAQUS/Standard and ABAQUS/Explicit; however, each program uses different algorithms to solve coupled thermal-stress problems.
In ABAQUS/Standard the temperatures are integrated using a backward-difference scheme, and the nonlinear coupled system is solved using Newton's method. ABAQUS/Standard offers an exact as well as an approximate implementation of Newton's method for fully coupled temperature-displacement analysis.
An exact implementation of Newton's method involves a nonsymmetric Jacobian matrix as is illustrated in the following matrix representation of the coupled equations:
Solving this system of equations requires the use of the unsymmetric matrix storage and solution scheme. Furthermore, the mechanical and thermal equations must be solved simultaneously. The method provides quadratic convergence when the solution estimate is within the radius of convergence of the algorithm. The exact implementation is used by default.
Some problems require a fully coupled analysis in the sense that the mechanical and thermal solutions evolve simultaneously, but with a weak coupling between the two solutions. In other words, the components in the off-diagonal submatrices are small compared to the components in the diagonal submatrices . An example of such a situation is the disc brake problem (Thermal-stress analysis of a disc brake, Section 5.1.1 of the ABAQUS Example Problems Manual). For these problems a less costly solution may be obtained by setting the off-diagonal submatrices to zero so that we obtain an approximate set of equations:
As a result of this approximation the thermal and mechanical equations can be solved separately, with fewer equations to consider in each subproblem. The savings due to this approximation, measured as solver time per iteration, will be of the order of a factor of two, with similar significant savings in solver storage of the factored stiffness matrix. Further, in many situations the subproblems may be fully symmetric or approximated as symmetric, so that the less costly symmetric storage and solution scheme can be used. The solver time savings for a symmetric solution is an additional factor of two. Unless you explicitly choose the unsymmetric matrix storage and solution scheme, selection of the scheme will depend on other details of the problem (see Procedures: overview, Section 6.1.1).
This modified form of Newton's method does not affect solution accuracy since the fully coupled effect is considered through the residual vector at each increment in time. However, the rate of convergence is no longer quadratic and depends strongly on the magnitude of the coupling effect, so more iterations are generally needed to achieve equilibrium than with the exact implementation of Newton's method. When the coupling is significant, the convergence rate becomes very slow and may prohibit obtaining a solution. In such cases the exact implementation of Newton's method is required. In cases where it is possible to use this approximation, the convergence in an increment will depend strongly on the quality of the first guess to the incremental solution, which you can control by selecting the extrapolation method used for the step (see Procedures: overview, Section 6.1.1).
|Input File Usage:||Use the following option to specify a separated solution scheme:|
*SOLUTION TECHNIQUE, TYPE=SEPARATED
Step module: Create Step: General: Coupled temp-displacement: Other: Solution technique: Separated
A steady-state coupled temperature-displacement analysis can be performed in ABAQUS/Standard. In steady-state cases you should assign an arbitrary “time” scale to the step: you choose fixed “time” increments and a “time” period. This time scale is convenient for changing loads and boundary conditions through the step and for obtaining solutions to highly nonlinear (but steady-state) cases; however, for the latter purpose, transient analysis often provides a natural way of coping with the nonlinearity.
Frictional slip heat generation is normally neglected in for the steady-state case. However, it can still be accounted for if motions are used to specify translational or rotational nodal velocities in disk brake-type problems or if user subroutine FRIC provides the incremental frictional dissipation through the variable SFD. If frictional heat generation is present, the heat flux into the two contact surfaces depends on the slip rate of the surfaces. The “time” scale in this case cannot be described as arbitrary, and a transient analysis should be performed.
|Input File Usage:|
*COUPLED TEMPERATURE-DISPLACEMENT, STEADY STATE
Step module: Create Step: General: Coupled temp-displacement: Basic: Response: Steady state
Alternatively, you can perform a transient coupled temperature-displacement analysis. You can control the time incrementation in a transient analysis directly, or ABAQUS/Standard can control it automatically. Automatic time incrementation is generally preferred.
The time increments can be selected automatically based on a user-prescribed maximum allowable nodal temperature change in an increment, . ABAQUS/Standard will restrict the time increments to ensure that this value is not exceeded at any node (except nodes with boundary conditions) during any increment of the analysis (see Time integration accuracy in transient problems, Section 7.2.4).
|Input File Usage:|
*COUPLED TEMPERATURE-DISPLACEMENT, DELTMX=
Step module: Create Step: General: Coupled temp-displacement: Basic: Response: Transient; Incrementation: Type: Automatic, Max. allowable temperature change per increment:
If you do not specify , fixed time increments equal to the user-specified initial time increment, , will be used throughout the analysis.
|Input File Usage:|
Step module: Create Step: General: Coupled temp-displacement: Basic: Response: Transient; Incrementation: Type: Fixed: Increment size:
In transient analysis with second-order elements there is a relationship between the minimum usable time increment and the element size. A simple guideline is
There is no upper limit on the time increment size (the integration procedure is unconditionally stable) unless nonlinearities cause convergence problems.
The accuracy of the integration of time-dependent (creep) material behavior is governed by the user-specified accuracy tolerance parameter, . This parameter is used to prescribe the maximum strain rate change allowed at any point during an increment, as described in Rate-dependent plasticity: creep and swelling, Section 18.2.4. The accuracy tolerance parameter can be specified together with the maximum allowable nodal temperature change in an increment, (described above); however, specifying the accuracy tolerance parameter activates automatic incrementation even if is not specified.
|Input File Usage:|
*COUPLED TEMPERATURE-DISPLACEMENT, DELTMX=, CETOL=tolerance
Step module: Create Step: General: Coupled temp-displacement: Basic: Response: Transient, Include creep/swelling/viscoelastic behavior; Incrementation: Type: Automatic, Max. allowable temperature change per increment: , Creep/swelling/viscoelastic strain error tolerance: tolerance
Nonlinear creep problems (Rate-dependent plasticity: creep and swelling, Section 18.2.4) that exhibit no other nonlinearities can be solved efficiently by forward-difference integration of the inelastic strains if the inelastic strain increments are smaller than the elastic strains. This explicit method is efficient computationally because, unlike implicit methods, iteration is not required as long as no other nonlinearities are present. Although this method is only conditionally stable, the numerical stability limit of the explicit operator is in many cases sufficiently large to allow the solution to be developed in a reasonable number of time increments.
For most coupled thermal-stress analyses, however, the unconditional stability of the backward difference operator (implicit method) is desirable. In such cases the implicit integration scheme may be invoked automatically by ABAQUS/Standard.
Explicit integration can be less expensive computationally and simplifies implementation of user-defined creep laws in user subroutine CREEP; you can restrict ABAQUS/Standard to using this method for creep problems (with or without geometric nonlinearity included). See Rate-dependent plasticity: creep and swelling, Section 18.2.4, for further details.
|Input File Usage:|
*COUPLED TEMPERATURE-DISPLACEMENT, CETOL=tolerance, CREEP=EXPLICIT
Step module: Create Step: General: Coupled temp-displacement: Basic: Response: Transient, Include creep/swelling/viscoelastic behavior; Incrementation: Type: Automatic, Creep/swelling/viscoelastic strain error tolerance: tolerance, Creep/swelling/viscoelastic integration: Explicit
You can specify that no creep or viscoelastic response will occur during a step even if creep or viscoelastic material properties have been defined.
|Input File Usage:|
*COUPLED TEMPERATURE-DISPLACEMENT, DELTMX=, CREEP=NONE
Step module: Create Step: General: Coupled temp-displacement: Basic: Response: Transient, toggle off Include creep/swelling/viscoelastic behavior
Some types of analyses may develop local instabilities, such as surface wrinkling, material instability, or local buckling. In such cases it may not be possible to obtain a quasi-static solution, even with the aid of automatic incrementation. ABAQUS/Standard offers a method of stabilizing this class of problems by applying damping throughout the model in such a way that the viscous forces introduced are sufficiently large to prevent instantaneous buckling or collapse, but small enough not to affect the behavior significantly while the problem is stable. ABAQUS/Standard generates an artificial damping matrix by using the mass matrix with a unit density together with a mass-proportional damping factor, as described in Solving nonlinear problems, Section 7.1.1. Whenever possible, the damping factor is chosen such that, based on extrapolation of the results obtained during the first increment, the dissipated energy during the step is a small fraction of the change in strain energy during the step. You control this dissipated energy fraction, which has a default value of 2.0 × 104. If the problem is either unstable or contains rigid body motions during the first increment, an alternative method is used to determine the damping factor; this method is based on making an averaged damping stiffness equal to the dissipated energy fraction times an averaged material stiffness.
|Input File Usage:||Use the following option to activate automatic stabilization with the default dissipated energy fraction:|
*COUPLED TEMPERATURE-DISPLACEMENT, STABILIZE
Use the following option to specify a nondefault dissipated energy fraction:
*COUPLED TEMPERATURE-DISPLACEMENT, STABILIZE=dissipated energy fraction
Use the following option to specify the damping factor directly:
*COUPLED TEMPERATURE-DISPLACEMENT, STABILIZE, FACTOR=damping factor
Step module: Create Step: General: Coupled temp-displacement: Basic: toggle on Use stabilization with, and select dissipated energy fraction or damping factor
In coupled problems where two different fields are active, take care when choosing the units of the problem. If the choice of units is such that the terms generated by the equations for each field are different by many orders of magnitude, the precision on some computers may be insufficient to resolve the numerical ill-conditioning of the coupled equations. Therefore, choose units that avoid ill-conditioned matrices. For example, consider using units of MPascal instead of Pascal for the stress equilibrium equations to reduce the disparity between the magnitudes of the stress equilibrium equations and the heat flux continuity equations.
In ABAQUS/Explicit the heat transfer equations are integrated using the explicit forward-difference time integration rule
The mechanical solution response is obtained using the explicit central-difference integration rule with a lumped mass matrix as described in Explicit dynamic analysis, Section 6.3.3. Since both the forward-difference and central-difference integrations are explicit, the heat transfer and mechanical solutions are obtained simultaneously by an explicit coupling. Therefore, no iterations or tangent stiffness matrices are required.
Explicit integration can be less expensive computationally and simplifies the treatment of contact. For a comparison of explicit and implicit direct-integration procedures, see Dynamic analysis procedures: overview, Section 6.3.1.
The explicit procedure integrates through time by using many small time increments. The central-difference and forward-difference operators are conditionally stable. The stability limit for both operators (with no damping in the mechanical solution response) is obtained by choosing
An approximation to the stability limit for the forward-difference operator in the thermal solution response is given by
In most applications of explicit analysis the mechanical response will govern the stability limit. The thermal response may govern the stability limit when material parameter values are non-physical or a very large amount of mass scaling is used. The calculation of the time increment size for the mechanical solution response is discussed in Explicit dynamic analysis, Section 6.3.3.
ABAQUS/Explicit writes a report to the status (.sta) file during the data check phase of the analysis that contains an estimate of the minimum stable time increment and a listing of the elements with the smallest stable time increments and their values. The initial minimum stable time increment accounts for the stability requirements of both the thermal and mechanical solution responses. The initial stable time increments listed do not include damping (bulk viscosity), mass scaling, or penalty contact effects in the mechanical solution response.
This listing is provided because often a few elements have much smaller stability limits than the rest of the elements in the mesh. The stable time increment can be increased by modifying the mesh to increase the size of the controlling element or by using appropriate mass scaling.
The time increment used in an analysis must be smaller than the stability limits of the central- and forward-difference operators. Failure to use such a time increment will result in an unstable solution. When the solution becomes unstable, the time history response of solution variables, such as displacements, will usually oscillate with increasing amplitudes. The total energy balance will also change significantly.
ABAQUS/Explicit has two strategies for time incrementation control: fully automatic time incrementation (where the code accounts for changes in the stability limit) and fixed time incrementation.
To reduce the chance of a solution going unstable, the stable time increment computed by ABAQUS/Explicit can be adjusted by a constant scaling factor. This factor can be used to scale the default global time estimate, the element-by-element estimate, or the fixed time increment based on the initial element-by-element estimate; it cannot be used to scale a fixed time increment that you specified directly.
Step module: Create Step: General: Dynamic, Temp-disp, Explicit: Incrementation: Time scaling factor: f
The default time incrementation scheme in ABAQUS/Explicit is fully automatic and requires no user intervention. Two types of estimates are used to determine the stability limit: element-by-element for both the thermal and mechanical solution responses and global for the mechanical solution response. An analysis always starts by using the element-by-element estimation method and may switch to the global estimation method under certain circumstances, as explained in Explicit dynamic analysis, Section 6.3.3.
In an analysis ABAQUS/Explicit initially uses a stability limit based on the thermal and mechanical solution responses in the whole model. This element-by-element estimate is determined using the smallest time increment size due to the thermal and mechanical solution responses in each element.
The element-by-element estimate is conservative; it will give a smaller stable time increment than the true stability limit, which is based upon the maximum frequency of the entire model. In general, constraints such as boundary conditions and kinematic contact have the effect of compressing the eigenvalue spectrum, and the element-by-element estimates do not take this into account (see Explicit dynamic analysis, Section 6.3.3)
The stable time increment size due to the mechanical solution response will be determined by the global estimator as the step proceeds unless the element-by-element estimator is chosen, fixed time incrementation is specified, or one of the conditions explained in Explicit dynamic analysis, Section 6.3.3, prevents the use of global estimation. The stable time increment size due to the thermal solution response will always be determined by using an element-by-element estimation method. The switch to the global estimation method in mechanical solution response occurs once the algorithm determines that the accuracy of the global estimation method is acceptable. For details, see Explicit dynamic analysis, Section 6.3.3.
|Input File Usage:||Use the following option to specify the element-by-element estimation method:|
*DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT, ELEMENT BY ELEMENT
Step module: Create Step: General: Dynamic, Temp-disp, Explicit: Incrementation: Type: Automatic, Stable increment estimator: Element-by-element
A fixed time incrementation scheme is also available in ABAQUS/Explicit. The fixed time increment size is determined either by the initial element-by-element stability estimate for the step or by a user-specified time increment.
Fixed time incrementation may be useful when a more accurate representation of the higher mode response of a problem is required. In this case a time increment size smaller than the element-by-element estimates may be used. The element-by-element estimate can be obtained simply by running a data check analysis (see Execution procedure for ABAQUS/Standard and ABAQUS/Explicit, Section 3.2.2).
When fixed time incrementation is used, ABAQUS/Explicit will not check that the computed response is stable during the step. You should ensure that a valid response has been obtained by carefully checking the energy history and other response variables.
If you choose to use time increments the size of the initial element-by-element stability limit throughout a step, the dilatational wave speed and the thermal diffusivity in each element at the beginning of the step are used to compute the fixed time increment size. To reduce the chance of a solution going unstable, the initial stable time increment that ABAQUS/Explicit computes can be adjusted by a constant scaling factor, as described above in “Scaling the time increment.” Alternatively, you can specify a time increment size directly.
|Input File Usage:||Use the following option to request time increments the size of the element-by-element stability limit:|
*DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT, FIXED TIME INCREMENTATION
Use the following option to specify the time increment size directly:
*DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT, DIRECT USER CONTROL
Step module: Create Step: General: Dynamic, Temp-disp, Explicit: Incrementation: Type: Fixed, Use element-by-element time increment estimator or User-defined time increment:
The extreme values defined as the element and nodal variables in a coupled thermal-stress analysis can be monitored exactly as described in Explicit dynamic analysis, Section 6.3.3, for a pure mechanical analysis.
By default, the initial temperature of all nodes is zero. You can specify nonzero initial temperatures. Initial stresses, field variables, etc. can also be defined; Initial conditions, Section 27.2.1, describes all of the initial conditions that are available for a fully coupled thermal-stress analysis.
Boundary conditions can be used to prescribe both temperatures (degree of freedom 11) and displacements/rotations (degrees of freedom 1–6) at nodes in fully coupled thermal-stress analysis (see Boundary conditions, Section 27.3.1). Shell elements in ABAQUS/Standard have additional temperature degrees of freedom 12, 13, etc. through the thickness (see Conventions, Section 1.2.2).
Boundary conditions can be specified as functions of time by referring to amplitude curves (Amplitude curves, Section 27.1.2).
Boundary conditions applied during a dynamic coupled temperature-displacement response step should use appropriate amplitude references (Amplitude curves, Section 27.1.2). If boundary conditions are specified for the step without amplitude references, they are applied instantaneously at the beginning of the step. Since ABAQUS/Explicit does not admit jumps in displacement, the value of a nonzero displacement boundary condition that is specified without an amplitude reference will be ignored, and a zero velocity boundary condition will be enforced.
The following types of thermal loads can be prescribed in a fully coupled thermal-stress analysis, as described in Thermal loads, Section 27.4.4:
Concentrated heat fluxes
Body fluxes and distributed surface fluxes
Node-based film conditions and radiation conditions
Convective film conditions and radiation conditions
Concentrated nodal forces can be applied to the displacement degrees of freedom (1–6); see Concentrated loads, Section 27.4.2.
Predefined temperature fields are not allowed in a fully coupled thermal-stress analysis. Boundary conditions should be used instead to prescribe temperature degree of freedom 11 (and 12, 13, etc. in ABAQUS/Standard shell elements), as described earlier.
Other predefined field variables can be specified in a fully coupled thermal-stress analysis. These values will affect only field-variable-dependent material properties, if any. See Predefined fields, Section 27.6.1.
The materials in a fully coupled thermal-stress analysis must have both thermal properties, such as conductivity, and mechanical properties, such as elasticity, defined. See Part V, Materials,” for details on the material models available in ABAQUS.
In ABAQUS/Standard internal heat generation can be specified; see Uncoupled heat transfer analysis, Section 6.5.2.
Thermal strain will arise if thermal expansion (Thermal expansion, Section 20.1.2) is included in the material property definition.
In ABAQUS/Standard a fully coupled temperature-displacement analysis can be used to analyze static creep and swelling problems, which generally occur over fairly long time periods (Rate-dependent plasticity: creep and swelling, Section 18.2.4); viscoelastic materials (Time domain viscoelasticity, Section 17.7.1); or viscoplastic materials (Rate-dependent yield, Section 18.2.3).
You can specify an inelastic heat fraction in a fully coupled thermal-stress analysis to provide for inelastic energy dissipation as a heat source. Plastic straining gives rise to a heat flux per unit volume of
Inelastic heat fractions are typically used in the simulation of high-speed manufacturing processes involving large amounts of inelastic strain, where the heating of the material caused by its deformation significantly influences temperature-dependent material properties. The generated heat is treated as a volumetric heat flux source term in the heat balance equation.
An inelastic heat fraction can be specified for materials with plastic behavior that use the Mises or Hill yield surface (Inelastic behavior, Section 18.1.1). It cannot be used with the combined isotropic/kinematic hardening model. The inelastic heat fraction can be specified for user-defined material behavior in ABAQUS/Explicit and will be multiplied by the inelastic energy dissipation coded in the user subroutine to obtain the heat flux. In ABAQUS/Standard the inelastic heat fraction cannot be used with user-defined material behavior; in this case the heat flux that must be added to the thermal energy balance is computed directly in the user subroutine.
The default value of the inelastic heat fraction is 0.9. If you do not include the inelastic heat fraction behavior in the material definition, the heat generated by inelastic deformation is not included in the analysis.
|Input File Usage:|
Property module: material editor: Thermal: Inelastic Heat Fraction: Fraction:
Coupled temperature-displacement elements that have both displacements and temperatures as nodal variables are available in both ABAQUS/Standard and ABAQUS/Explicit (see Choosing the appropriate element for an analysis type, Section 21.1.3). In ABAQUS/Standard simultaneous temperature/displacement solution requires the use of such elements; pure displacement elements can be used in part of the model in the fully coupled thermal-stress procedure, but pure heat transfer elements cannot be used. In ABAQUS/Explicit any of the available elements can be used in the fully coupled thermal-stress procedure; however, the thermal solution will be obtained only at nodes where the temperature degree of freedom has been activated (i.e., at nodes attached to coupled temperature-displacement elements).
The first-order coupled temperature-displacement elements in ABAQUS use a constant temperature over the element to calculate thermal expansion. The second-order coupled temperature-displacement elements in ABAQUS/Standard use a lower-order interpolation for temperature than for displacement (parabolic variation of displacements and linear variation of temperature) to obtain a compatible variation of thermal and mechanical strain.
See ABAQUS/Standard output variable identifiers, Section 4.2.1, and ABAQUS/Explicit output variable identifiers, Section 4.2.2, for a complete list of output variables. The types of output available are described in Output, Section 4.1.1.
*HEADING … ** Specify the coupled temperature-displacement element type *ELEMENT, TYPE=CPS4T … ** *STEP *COUPLED TEMPERATURE-DISPLACEMENT or *DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT Data line to define incrementation *BOUNDARY Data lines to define nonzero boundary conditions on displacement or temperature degrees of freedom *CFLUX and/or *CFILM and/or *CRADIATE and/or *DFLUX and/or *DSFLUX and/or *FILM and/or *SFILM and/or *RADIATE and/or *SRADIATE Data lines to define thermal loads *CLOAD and/or *DLOAD and/or *DSLOAD Data lines to define mechanical loads *FIELD Data lines to define field variable values *END STEP