Pore pressure submodeling

This problem contains basic test cases for one or more Abaqus elements and features.

This page discusses:

ProductsAbaqus/Standard

Elements tested

  • C3D8PH
  • C3D8PT
  • C3D10MP
  • C3D20P
  • C3D20RP
  • CAX4RP
  • CAX6MP
  • CAX8P
  • CAX8RP
  • CPE4P
  • CPE6MP
  • CPE8P
  • CPE8RP

Features tested

The submodeling capability is applied to two-dimensional, three-dimensional and axisymmetric continuum elements with pore pressure. General steps invoking the soils consolidation procedure are used for both the global and submodel analyses.

Problem description

Model:

All global models have dimensions 3.0 × 5.0 in the xy or rz plane. Each submodel has dimensions 2.05 × 3.45 in the xy or rz plane and occupies the lower right-hand corner of the corresponding global model. In all but the axisymmetric models, the out-of-plane dimension is 1.0. In axisymmetric models the structure analyzed is a hollow cylinder with outer radius 5.0.

Material:

Young's modulus 100 × 106
Poisson's ratio 0.0
Permeability 1 × 10−5
Density 1.4142

Loading:

In all models, a distributed flow of magnitude 0.002 is applied to the right face, where the sink pore pressure is 14.7.

Boundary and initial conditions

In the global model, fixed boundary conditions u1 = 0 and u2 = 0 are prescribed on the left and bottom faces, respectively. In three-dimensional models the additional constraints u3 = 0 are applied to the nodes on the front and back faces. The initial void ratio is unity everywhere and fixed pore pressure boundary conditions are applied on the left face. In the submodel, u2 = 0 is prescribed everywhere on the bottom face, while degrees of freedom 1, 2 and 8 for the nodes on the top and left faces are being driven by the global solution.

Results and discussion

In the global analyses, the pore pressure field predicted by Abaqus varies linearly in the x-direction in nonaxisymmetric models and logarithmically in the r-direction in axisymmetric models. The predicted displacement field is nonuniform in all models. These results are depicted in the pore pressure and x- or r-displacement contour plots shown below. For comparison purposes, the pore pressure and displacement solutions predicted by the submodels are also presented in the same contour plots and excellent agreement between global and submodel results is obtained. Hence, the amplitudes of all driven variables in the submodeled analysis are correctly identified in the global analysis file output and applied at the driven nodes in the submodel analysis.

Global and submodel analyses results for 8-node plane strain elements are shown in Figure 1 and Figure 2.

Global and submodel analyses results for 8-node axisymmetric elements are shown in Figure 3 and Figure 4.

Global and submodel analyses results for 20-node brick elements (front face) are shown in Figure 5 and Figure 6.

Input files

The following input files test the transient soils consolidation procedure. Each test performs a single increment transient consolidation calculation for a time period of one.
pgc38phd.inp

C3D8PH elements; global analysis.

psc38phd.inp

C3D8PH elements; submodel analysis.

pgc3apkd.inp

C3D10MP elements; global analysis.

psc3apkd.inp

C3D10MP elements; submodel analysis.

pgc3kpfd.inp

C3D20P elements; global analysis.

psc3kpfd.inp

C3D20P elements; submodel analysis.

pgc3kprd.inp

C3D20RP elements; global analysis.

psc3kprd.inp

C3D20RP elements; submodel analysis.

pgca4prd.inp

CAX4RP elements; global analysis.

psca4prd.inp

CAX4RP elements; submodel analysis.

pgca6pkd.inp

CAX6MP elements; global analysis.

psca6pkd.inp

CAX6MP elements; submodel analysis.

pgca8pfd.inp

CAX8P elements; global analysis.

psca8pfd.inp

CAX8P elements; submodel analysis.

pgca8prd.inp

CAX8RP elements; global analysis.

psca8prd.inp

CAX8RP elements; submodel analysis.

pgce4pfd.inp

CPE4P elements; global analysis.

psce4pfd.inp

CPE4P elements; submodel analysis.

pgce6pkd.inp

CPE6MP elements; global analysis.

psce6pkd.inp

CPE6MP elements; submodel analysis.

pgce8pfd.inp

CPE8P elements; global analysis.

psce8pfd.inp

CPE8P elements; submodel analysis.

pgce8prd.inp

CPE8RP elements; global analysis.

psce8prd.inp

CPE8RP elements; submodel analysis.

ctp_gbmodel.inp

C3D8PT elements; global analysis.

ctp_sbmodel.inp

C3D8PT elements; submodel analysis.

Figures

Figure 1. Pore pressure contours in global and submodels: 8-node plane strain.

Figure 2. ux contours in global and submodels: 8-node plane strain.

Figure 3. Pore pressure contours in global and submodels: 8-node axisymmetric.

Figure 4. ur contours in global and submodels: 8-node axisymmetric.

Figure 5. Pore pressure contours in global and submodels: 20-node brick.

Figure 6. ux contours in global and submodels: 20-node brick.