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.

This task shows you how to:

Create or edit a coupled pore fluid diffusion/stress procedure

  1. Display the Edit Step dialog box following the procedure outlined in Creating a step (Procedure type: General; Soils), or Editing a step.

  2. On the Basic, Incrementation, and Other tabbed pages, configure settings such as steady-state or transient pore fluid response and automatic or fixed incrementation as described in the following procedures.

Configure settings on the Basic tabbed page

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

  2. In the Description field, enter a short description of the 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. Choose a Pore fluid response option:

    • Choose Steady-state to specify that there are no transient effects in the wetting liquid continuity equation. The steady-state solution corresponds to constant wetting liquid velocities and constant volume of wetting liquid per unit volume in the continuum. For more information, see Steady-state analysis.

    • Choose Transient consolidation to use the backward difference operator to integrate the continuity equation. This operator provides unconditional stability so that the only concern with respect to time integration is accuracy. For more information, see Transient analysis.

    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 the time period of the step.

  5. Select an Nlgeom option:

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

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

    For more information, see Linear and nonlinear procedures.

  6. Select an automatic stabilization method 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 Unstable problems, and Automatic stabilization of static problems with a constant damping factor

    Click the arrow to the right of Automatic stabilization, and select a method for defining the damping factor:

    • Select Specify dissipated energy fraction to allow Abaqus/Standard to calculate the damping factor from a dissipated energy fraction that you provide. Enter a value for the dissipated energy fraction in the adjacent field (the default is 2.0 × 10−4). For more information, see Calculating the damping factor based on the dissipated energy fraction.

    • Select Specify damping factor to enter the damping factor directly. Enter a value for the damping factor in the adjacent field. For more information, see Directly specifying the damping factor.

    • Select Use damping factors from previous general step to use the damping factors at the end of the previous step as the initial factors in the current step's variable damping scheme. These factors override any initial damping factors that are calculated or specified directly in the current step. If there are no damping factors associated with the previous general step (for example, if the previous step does not use any stabilization or the current step is the first step of the analysis), Abaqus uses adaptive stabilization to determine the required damping factors.

  7. When using automatic stabilization, Abaqus can use the same damping factor over the course of a step, or it can vary the damping factor spatially and temporally during a step based on the convergence history and the ratio of the energy dissipated by damping to the total strain energy. For more information, see Adaptive automatic stabilization scheme. If you selected Specify dissipated energy fraction, adaptive stabilization is optional and turned on by default. If you selected Specify damping factor, adaptive stabilization is optional and turned off by default. If you selected Use damping factors from previous general step, adaptive stabilization is required.

    To use adaptive stabilization, toggle on Use adaptive stabilization with max. ratio of stabilization to strain energy (if necessary), and enter a value in the adjacent field for the allowable accuracy tolerance for the ratio of energy dissipated by damping to total strain energy in each increment. The default value of 0.05 should be suitable in most cases.

  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.

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, or Editing a step.)

  2. Choose a Type option:

    • Choose Automatic if you want Abaqus/Standard to determine suitable time increment sizes.

    • Choose 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:

    1. In the Initial field, enter the initial time increment. Abaqus/Standard modifies this value as required throughout the step.
    2. 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.
    3. 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 the Transient consolidation 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 the Transient consolidation 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.

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, or Editing a step.)

  2. Choose an Equation Solver Method option:

    • Choose Direct to use the default direct sparse solver.

    • Choose Iterative to use the iterative linear equation solver. The iterative solver is typically most useful for blocky structures with millions of degrees of freedom. For more information, see Iterative linear equation solver.

  3. Choose a Matrix storage option:

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

    • Choose 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 Defining an analysis).

    • Choose 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.

  4. Choose a Solution technique:

    • Choose 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.

    • Choose 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 choose 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.

  5. 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, regardless of the magnitude of the penetration and force errors. 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 On to use local convergence criteria to determine whether a new iteration is needed. Abaqus/Standard will determine the maximum penetration and estimated force errors associated with severe discontinuities and check whether these errors are within the tolerances. Hence, a solution may converge if the severe discontinuities are small.

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

    For more information on severe discontinuities, see Severe discontinuities in Abaqus/Standard.

  6. 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.

  7. 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.

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