14.11.15 Configuring an effective stress analysis for fluid-filled porous media

A coupled pore fluid diffusion/stress analysis allows you to model single phase, partially or fully saturated fluid flow through porous media. For more information, see Coupled pore fluid diffusion and stress analysis, Section 6.7.1 of the ABAQUS Analysis User's Manual.

When you configure this type of procedure, the step editor displays the Basic, Incrementation, and Other tabs. Settings you can configure with these tabbed pages include steady-state or transient pore fluid response and automatic or fixed incrementation.

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. Select a Pore fluid response option:

    Note:  After you have selected a Pore fluid response option, a message appears informing you that ABAQUS/Standard has selected the Default load variation with time option and the Matrix storage option (both located on the Other tabbed page) that correspond to your Pore fluid 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.

      Note:  Fixed incrementation is not generally recommended in this case because the time increments in a typical diffusion analysis can increase over several orders of magnitude during the simulation; automatic incrementation is usually a better choice.

  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 Transient response on the Basic tabbed page, toggle on End step when pore pressure change rate is less than n to enter a minimum value for the pore pressure change rate. The analysis will end if all pore pressures are changing at a rate that is less than the rate that you enter.

  7. If you selected Automatic in Step 2 , do the following:

    1. If you selected Transient response on the Basic tabbed page, enter a value for the Max. pore pressure change per increment. ABAQUS/Standard restricts the time step to ensure that this value is not exceeded at any node (except nodes with boundary conditions) 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 Specifying the tolerance for automatic incrementation” in “Rate-dependent plasticity: creep and swelling, Section 18.2.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.

      Note:  The steady-state coupled equations are strongly unsymmetric; therefore, the unsymmetric matrix solution and storage scheme is selected automatically for steady-state analysis steps (see Procedures: overview, Section 6.1.1 of the ABAQUS Analysis User's Manual).

    • 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 Quasi-Newton to use the quasi-Newton technique for solving nonlinear equilibrium equations. This technique can save substantial computational cost in some cases. Generally it is most successful when the system is large and the stiffness matrix is not changing much from iteration to iteration. You can use this technique only for symmetric systems of equations.

      If you select this technique, enter a value for the Number of iterations allowed before the kernel matrix is reformed. The maximum number of iterations allowed is 25. The default number of iterations is 8.

      For more information, see Quasi-Newton solution technique, Section 2.2.2 of the ABAQUS Theory 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 Pore fluid 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.