14.11.10 Configuring a fully coupled, simultaneous heat transfer and stress procedure

You must configure a fully coupled temperature-displacement analysis 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. For such cases the thermal and mechanical solutions must be obtained simultaneously rather than sequentially. For more information, see Fully coupled thermal-stress analysis, Section 6.5.4 of the ABAQUS Analysis User's Manual.

When you configure a fully coupled temperature-displacement procedure, the step editor displays the Basic, Incrementation, and Other tabs. Settings you can configure with these tabbed pages include the time period for the step, increment size, and solution technique preferences.

To configure settings on the Basic tabbed page:

  1. In the Edit Step dialog box, display the Basic tabbed page.

    (For information on displaying the Edit Step dialog box, see Creating a step, Section 14.9.2, or Editing a step, Section 14.9.3.)

  2. In the Description field, enter a short description of what occurs during this analysis step. ABAQUS stores the text that you enter in the output database, and the text is displayed in the state block by the Visualization module.

  3. Indicate whether you want Steady-state or Transient response. See the following sections for more information:

    Note:  After you have selected a Response option, a message appears informing you that ABAQUS/Standard has selected the Default load variation with time option (located on the Other tabbed page) that corresponds to your Response selection. Click Dismiss to close the message dialog box.

  4. In the Time period field, enter a time scale for the analysis.

  5. Select an Nlgeom option:

    • Turn Nlgeom Off to perform a geometrically linear analysis during the current step.

    • Turn Nlgeom On to indicate that ABAQUS/Standard should account for geometric nonlinearity during the step. Once you have turned Nlgeom on, it will be active during all subsequent steps in the analysis.

    For more information, see Linear and nonlinear procedures, Section 14.3.2.

  6. Toggle on Use stabilization if you expect the problem to have local instabilities such as surface wrinkling, material instability, or local buckling. ABAQUS/Standard can stabilize this class of problems by applying damping throughout the model. For more information, see the follow sections:

    .

  7. If you toggled on Use stabilization, click the arrow to the right of the Use stabilization field, and select a method for defining the damping factor:

  8. If desired, toggle on Include creep/swelling/viscoelastic behavior. If you leave this option toggled off, you indicate that there is no creep or viscoelastic response occurring during this step even if creep or viscoelastic material properties have been defined.

To configure settings on the Incrementation tabbed page:

  1. In the Edit Step dialog box, display the Incrementation tabbed page.

    (For information on displaying the Edit Step dialog box, see Creating a step, Section 14.9.2, or Editing a step, Section 14.9.3.)

  2. Select a Type option:

    • Select Automatic if you want ABAQUS/Standard to determine suitable time increment sizes.

    • Select Fixed to specify direct user control of the incrementation. ABAQUS/Standard uses an increment size that you specify as the constant increment size throughout the step.

  3. In the Maximum number of increments field, enter the upper limit to the number of increments in the step. The analysis stops if this maximum is exceeded before ABAQUS/Standard arrives at the complete solution for the step.

  4. If you selected Automatic in Step 2, enter values for Increment size:

    • In the Initial field, enter the initial time increment. ABAQUS/Standard modifies this value as required throughout the step.

    • In the Minimum field, enter the minimum time increment allowed. If ABAQUS/Standard needs a smaller time increment than this value, it terminates the analysis.

    • In the Maximum field, enter the maximum time increment allowed.

  5. If you selected Fixed in Step 2, enter a value for the constant time increment in the Increment size field.

  6. If you selected Automatic in Step 2 and if you selected Transient response on the Basic tabbed page, do the following:

    1. Enter a value for the Max. allowable temperature change per increment. ABAQUS/Standard restricts the time step to ensure that this value is not exceeded at any node during any increment of the step.

    2. If you toggled on Include creep/swelling/viscoelastic behavior on the Basic tabbed page, toggle on Creep/swelling/viscoelastic strain error tolerance to enter the maximum difference in the creep strain increment calculated from the creep strain rates at the beginning and at the end of the increment. This value controls the accuracy of the creep integration. For more information, see Automatic incrementation controlled by the creep response” in “Fully coupled thermal-stress analysis, Section 6.5.4 of the ABAQUS Analysis User's Manual.

  7. If you toggled on Include creep/swelling/viscoelastic behavior on the Basic tabbed page, select a Creep/swelling/viscoelastic integration option:

    • Select Explicit/Implicit if you want to allow ABAQUS/Standard to invoke the implicit integration scheme. For most coupled thermal-stress analyses, the unconditional stability of the backward difference operator (implicit method) is desirable.

    • Select Explicit if you want to restrict ABAQUS/Standard to using explicit integration. Explicit integration can be less expensive computationally and simplifies implementation of user-defined creep laws in user subroutine CREEP

    For more information, see Automatic incrementation controlled by the creep response” in “Fully coupled thermal-stress analysis, Section 6.5.4 of the ABAQUS Analysis User's Manual.

To configure settings on the Other tabbed page:

  1. In the Edit Step dialog box, display the Other tabbed page.

    (For information on displaying the Edit Step dialog box, see Creating a step, Section 14.9.2, or Editing a step, Section 14.9.3.)

  2. Select the Matrix storage option of your choice:

    • Select Use solver default to allow ABAQUS/Standard to decide whether a symmetric or unsymmetric matrix storage and solution scheme is needed.

    • Select Unsymmetric to restrict ABAQUS/Standard to the unsymmetric storage and solution scheme. (This is the only matrix storage option available if you select the Full Newton solution technique.)

    • Select Symmetric to restrict ABAQUS/Standard to the symmetric storage and solution scheme.

    For more information on matrix storage, see Matrix storage and solution scheme in ABAQUS/Standard” in “Procedures: overview, Section 6.1.1 of the ABAQUS Analysis User's Manual.

  3. Select the Solution technique of your choice:

    • Select Full Newton to use Newton's method as a numerical technique for solving nonlinear equilibrium equations. For more information, see Nonlinear solution methods in ABAQUS/Standard, Section 2.2.1 of the ABAQUS Theory Manual.

    • Select Separated to specify that linearized equations for the individual fields in the fully coupled procedure are to be decoupled and solved separately for each field. This option provides a less costly solution for an analysis that is fully coupled in the sense that the mechanical and thermal solutions evolve simultaneously, but with a weak coupling between the two solutions. For more information, see Approximate implementation” in “Fully coupled thermal-stress analysis, Section 6.5.4 of the ABAQUS Analysis User's Manual.

    • Select Contact iterations to use contact iterations instead of regular severe discontinuity iterations to speed up computations. Contact iterations are effective for the solution of large, geometrically linear, small-sliding, frictionless static problems with many severe discontinuity iterations.

      If you select this technique, enter the following values:

      • Adjustment factor for the number of solutions in any iteration. This value is a correction factor on the maximum number of right-hand-side solutions during any contact iteration.

      • Maximum number of contact iterations. This value specifies the maximum number of contact iterations allowed before new global matrix assemblage and factorization.

      For more information, see Contact iterations, Section 7.1.2 of the ABAQUS Analysis User's Manual.

  4. Click the arrow to the right of the Convert severe discontinuity iterations field, and select an option for dealing with severe discontinuities during nonlinear analysis:

    • Select Off to force a new iteration if severe discontinuities occur during an iteration.

    • Select On to estimate residual forces associated with severe discontinuities and to check whether the equilibrium tolerances are satisfied. A solution may converge if the severe discontinuities (such as penetrations or tensile contact forces) are small. This option also changes some time incrementation parameters and uses different criteria to determine whether to do another iteration or to make a new attempt with a smaller increment size.

    • Select Propagate from previous step to use the value specified in the previous general analysis step. This value appears in parenthesis to the right of the field.

  5. ABAQUS/Standard automatically selects the Default load variation with time option that corresponds to your Response selection on the Basic tabbed page. It is recommended that you leave the Default load variation with time selection unchanged.

  6. Click the arrow to the right of the Extrapolation of previous state at start of each increment field, and select a method for determining the first guess to the incremental solution:

    • Select Linear to indicate that the process is essentially monotonic and ABAQUS/Standard should use a 100% linear extrapolation, in time, of the previous incremental solution to begin the nonlinear equation solution for the current increment.

    • Select Parabolic to indicate that the process should use a quadratic extrapolation, in time, of the previous two incremental solutions to begin the nonlinear equation solution for the current increment.

    • Select None to suppress any extrapolation.

    For more information, see Extrapolation of the solution” in “Procedures: overview, Section 6.1.1 of the ABAQUS Analysis User's Manual.

Once you have finished configuring settings for the step, click OK to close the Edit Step dialog box.