5.2.2 Pressure and fluid flow in pore pressure contact

Product: ABAQUS/Standard  

ABAQUS/Standard provides a surface-based capability for modeling fully saturated porous media. The surface-based capability can be used for small or finite sliding. Only flow normal to the surface is considered; tangential flow cannot be modeled.

The pore fluid constraints on the contact interface

Let and be the pore pressures at the two sides of the interface. It is at all times required that the pore pressures on opposite sides be equal:

Similarly, let and be the volume flow rate densities normal to the interface at the two sides, and let be the relative velocity of the two sides in the direction of the interface normal.

It is assumed that the interface is filled with fluid at all times. Hence, continuity requires that

whereas the difference

is undetermined and is to be treated as an independent variable. Inversion of these equations yields

The transient equations

The contribution of the interfacial virtual work equation and its linearized form are first obtained in the general form including finite sliding. The equations are then specialized to the various formulations implemented in ABAQUS.

Since we want to achieve force and volume flow rate equilibrium at each side of the interface, as well as obtain continuity in the pore pressures, we add the following integral to the virtual work equation:

where is an arbitrary Lagrange multiplier and is the interface area. Eliminating and and using a suitable choice for ,

we obtain

where

Since ABAQUS uses displacements (not velocities) and fluxes integrated over , the equation can be multiplied by to obtain

where is the incremental change in in the direction of the interface normal. Linearization yields

Closed contact

If the two sides are locally in contact, and ; therefore, the virtual work simplifies to

Similarly, the linearized form simplifies to

The steady-state equations

For steady-state analysis the transient terms can be omitted, and the terms involving fluid flow are written in rate form. In this case we can assume that the interface displacements vanish, which leads to the simplified virtual work contribution

and the linearized form

Small sliding

When the small-sliding contact formulation is used, the terms , , , and in the linearized virtual work equation will vanish.

Reference