Abaqus Verification Guide  


1.4.7 Continuum pore pressure elements

Product: Abaqus/Standard  

I. Plane strain elements with pore pressure

Problem description

Model:

Planar dimension3 × 5
Gravity load vector(1, 1, 0)

Material:

Young's modulus1 × 108
Poisson's ratio0.0
Density1.4142
Permeability1 × 10–5
Specific weight of fluid1.0

Initial conditions:

Initial void ratio1.0
Hydrostatic pressure datum5.0
Hydrostatic pressure elevation0.0
Sink pore pressure14.7

Results and discussion

The calculated reactions are in agreement with the applied loads.

Input files

CPE4P element load tests:


ece4pfde.inp

CENTRIF, BX, BY, GRAV, P1, P2, P3, P4, HP1, HP2, HP3, HP4, Q1, Q2, Q3, Q4, S1, S2, S3, S4.

ece4pfdl.inp

P, HP, Q, S.

CPE4PH element load tests:


ece4phde.inp

CENTRIF, BX, BY, GRAV, P1, P2, P3, P4, HP1, HP2, HP3, HP4, Q1, Q2, Q3, Q4, S1, S2, S3, S4.

ece4phdl.inp

P, HP, Q, S.

CPE4RP element load tests:


ece4prde.inp

CENTRIF, BX, BY, GRAV, P1, P2, P3, P4, HP1, HP2, HP3, HP4, Q1, Q2, Q3, Q4, S1, S2, S3, S4.

ece4prdl.inp

P, HP, Q, S.

CPE4RPH element load tests:


ece4pyde.inp

CENTRIF, BX, BY, GRAV, P1, P2, P3, P4, HP1, HP2, HP3, HP4, Q1, Q2, Q3, Q4, S1, S2, S3, S4.

ece4pydl.inp

P, HP, Q, S.

CPE6MP element load tests:


ece6pkde.inp

CENTRIF, BX, BY, GRAV, P1, P2, P3, HP1, HP2, HP3, Q1, Q2, Q3, S1, S2, S3.

ece6pkdl.inp

P, HP, Q, S.

CPE6MPH element load tests:


ece6plde.inp

CENTRIF, BX, BY, GRAV, P1, P2, P3, HP1, HP2, HP3, Q1, Q2, Q3, S1, S2, S3.

ece6pldl.inp

P, HP, Q, S.

CPE8P element load tests:


ece8pfde.inp

CENTRIF, BX, BY, GRAV, P1, P2, P3, P4, HP1, HP2, HP3, HP4, Q1, Q2, Q3, Q4, Q1D, Q2D, Q3D, Q4D, S1, S2, S3, S4.

ece8pfdl.inp

P, HP, Q, QD, S.

CPE8PH element load tests:


ece8phde.inp

CENTRIF, BX, BY, GRAV, P1, P2, P3, P4, HP1, HP2, HP3, HP4, Q1, Q2, Q3, Q4, Q1D, Q2D, Q3D, Q4D, S1, S2, S3, S4.

ece8phdl.inp

P, HP, Q, QD, S.

CPE8RP element load tests:


ece8prde.inp

CENTRIF, BX, BY, GRAV, P1, P2, P3, P4, HP1, HP2, HP3, HP4, Q1, Q2, Q3, Q4, Q1D, Q2D, Q3D, Q4D, S1, S2, S3, S4.

ece8prdl.inp

P, HP, Q, QD, S.

CPE8RPH element load tests:


ece8pyde.inp

CENTRIF, BX, BY, GRAV, P1, P2, P3, P4, HP1, HP2, HP3, HP4, Q1, Q2, Q3, Q4, Q1D, Q2D, Q3D, Q4D, S1, S2, S3, S4.

ece8pydl.inp

P, HP, Q, QD, S.

II. Axisymmetric elements with pore pressure

Problem description

Model:

Planar dimension3 × 5
Inside radius1.0
Gravity load vector(1, 1, 0)

Material:

Young's modulus1 × 108
Poisson's ratio0.0
Density1.4142
Permeability1 × 10–5
Specific weight of fluid1.0

Initial conditions:

Initial void ratio1.0
Hydrostatic pressure datum5.0
Hydrostatic pressure elevation0.0
Sink pore pressure14.7

Results and discussion

The calculated reactions are in agreement with the applied loads.

Input files

CAX4P element load tests:


eca4pfde.inp

CENTRIF, BR, BZ, GRAV, P1, P2, P3, P4, HP1, HP2, HP3, HP4, Q1, Q2, Q3, Q4, S1, S2, S3, S4.

eca4pfdl.inp

P, HP, Q, S.

CAX4PH element load tests:


eca4phde.inp

CENTRIF, BR, BZ, GRAV, P1, P2, P3, P4, HP1, HP2, HP3, HP4, Q1, Q2, Q3, Q4, S1, S2, S3, S4.

eca4phdl.inp

P, HP, Q, S.

CAX4RP element load tests:


eca4prde.inp

CENTRIF, BR, BZ, GRAV, P1, P2, P3, P4, HP1, HP2, HP3, HP4, Q1, Q2, Q3, Q4, S1, S2, S3, S4.

eca4prdl.inp

P, HP, Q, S.

CAX4RPH element load tests:


eca4pyde.inp

CENTRIF, BR, BZ, GRAV, P1, P2, P3, P4, HP1, HP2, HP3, HP4, Q1, Q2, Q3, Q4, S1, S2, S3, S4.

eca4pydl.inp

P, HP, Q, S.

CAX6MP element load tests:


eca6pkde.inp

CENTRIF, BR, BZ, GRAV, P1, P2, P3, HP1, HP2, HP3, Q1, Q2, Q3, S1, S2, S3.

eca6pkdl.inp

P, HP, Q, S.

CAX6MPH element load tests:


eca6plde.inp

CENTRIF, BR, BZ, GRAV, P1, P2, P3, HP1, HP2, HP3, Q1, Q2, Q3, S1, S2, S3.

eca6pldl.inp

P, HP, Q, S.

CAX8P element load tests:


eca8pfde.inp

CENTRIF, BR, BZ, GRAV, P1, P2, P3, P4, HP1, HP2, HP3, HP4, Q1, Q2, Q3, Q4, Q1D, Q2D, Q3D, Q4D, S1, S2, S3, S4.

eca8pfdl.inp

P, HP, Q, QD, S.

CAX8PH element load tests:


eca8phde.inp

CENTRIF, BR, BZ, GRAV, P1, P2, P3, P4, HP1, HP2, HP3, HP4, Q1, Q2, Q3, Q4, Q1D, Q2D, Q3D, Q4D, S1, S2, S3, S4.

eca8phdl.inp

P, HP, Q, QD, S.

CAX8RP element load tests:


eca8prde.inp

CENTRIF, BR, BZ, GRAV, P1, P2, P3, P4, HP1, HP2, HP3, HP4, Q1, Q2, Q3, Q4, Q1D, Q2D, Q3D, Q4D, S1, S2, S3, S4.

eca8prdl.inp

P, HP, Q, QD, S.

CAX8RPH element load tests:


eca8pyde.inp

CENTRIF, BR, BZ, GRAV, P1, P2, P3, P4, HP1, HP2, HP3, HP4, Q1, Q2, Q3, Q4, Q1D, Q2D, Q3D, Q4D, S1, S2, S3, S4.

eca8pydl.inp

P, HP, Q, QD, S.

III. Three-dimensional elements with pore pressure

Problem description

Model:

Cubic dimension3 × 5 × 1
Gravity load vector(1, 1, 1)

Material:

Young's modulus1 × 108
Poisson's ratio0.0
Density1.7321
Permeability1 × 10–5
Specific weight of fluid1.0

Initial conditions:

Initial void ratio1.0
Hydrostatic pressure datum5.0
Hydrostatic pressure elevation0.0
Sink pore pressure14.7

Results and discussion

The calculated reactions are in agreement with the applied loads.

Input files

C3D4P element load tests:


elem_load_c3d4p_1.inp

CENTRIF, BX, BY, BZ, GRAV, P1, P2, P3, P4, HP1, HP2, HP3, HP4, Q1, Q2, Q3, Q4, S1, S2, S3, S4.

C3D6P element load tests:


elem_load_c3d6p_1.inp

CENTRIF, BX, BY, BZ, GRAV, P1, P2, P3, P4, P5, HP1, HP2, HP3, HP4, HP5, Q1, Q2, Q3, Q4, Q5, S1, S2, S3, S4, S5.

C3D8P element load tests:


ec38pfde.inp

CENTRIF, BX, BY, BZ, GRAV, P1, P2, P3, P4, P5, P6, HP1, HP2, HP3, HP4, HP5, HP6, Q1, Q2, Q3, Q4, Q5, Q6, S1, S2, S3, S4, S5, S6.

ec38pfdl.inp

P, HP, Q, S.

C3D8PH element load tests:


ec38phde.inp

CENTRIF, BX, BY, BZ, GRAV, P1, P2, P3, P4, P5, P6, HP1, HP2, HP3, HP4, HP5, HP6, Q1, Q2, Q3, Q4, Q5, Q6, S1, S2, S3, S4, S5, S6.

ec38phdl.inp

P, HP, Q, S.

C3D8RP element load tests:


ec38prde.inp

CENTRIF, BX, BY, BZ, GRAV, P1, P2, P3, P4, P5, P6, HP1, HP2, HP3, HP4, HP5, HP6, Q1, Q2, Q3, Q4, Q5, Q6, S1, S2, S3, S4, S5, S6.

ec38prdl.inp

P, HP, Q, S.

C3D8RPH element load tests:


ec38pyde.inp

CENTRIF, BX, BY, BZ, GRAV, P1, P2, P3, P4, P5, P6, HP1, HP2, HP3, HP4, HP5, HP6, Q1, Q2, Q3, Q4, Q5, Q6, S1, S2, S3, S4, S5, S6.

ec38pydl.inp

P, HP, Q, S.

C3D10MP element load tests:


ec3apkde.inp

CENTRIF, BX, BY, BZ, GRAV, P1, P2, P3, P4, HP1, HP2, HP3, HP4, Q1, Q2, Q3, Q4, S1, S2, S3, S4.

ec3apkdl.inp

P, HP, Q, S.

C3D10MPH element load tests:


ec3aplde.inp

CENTRIF, BX, BY, BZ, GRAV, P1, P2, P3, P4, HP1, HP2, HP3, HP4, Q1, Q2, Q3, Q4, S1, S2, S3, S4.

ec3apldl.inp

P, HP, Q, S.

C3D20P element load tests:


ec3kpfde.inp

CENTRIF, BX, BY, BZ, GRAV, P1, P2, P3, P4, P5, P6, HP1, HP2, HP3, HP4, HP5, HP6, Q1, Q2, Q3, Q4, Q5, Q6, Q1D, Q2D, Q3D, Q4D, Q5D, Q6D, S1, S2, S3, S4, S5, S6.

ec3kpfde.f

User subroutines FLOW and DFLOW used in ec3kpfde.inp.

ec3kpfdl.inp

P, HP, Q, QD, S.

C3D20PH element load tests:


ec3kphde.inp

CENTRIF, BX, BY, BZ, GRAV, P1, P2, P3, P4, P5, P6, HP1, HP2, HP3, HP4, HP5, HP6, Q1, Q2, Q3, Q4, Q5, Q6, Q1D, Q2D, Q3D, Q4D, Q5D, Q6D, S1, S2, S3, S4, S5, S6.

ec3kphdl.inp

P, HP, Q, QD, S.

C3D20RP element load tests:


ec3kprde.inp

CENTRIF, BX, BY, BZ, GRAV, P1, P2, P3, P4, P5, P6, HP1, HP2, HP3, HP4, HP5, HP6, Q1, Q2, Q3, Q4, Q5, Q6, Q1D, Q2D, Q3D, Q4D, Q5D, Q6D, S1, S2, S3, S4, S5, S6.

ec3kprdl.inp

P, HP, Q, QD, S.

C3D20RPH element load tests:


ec3kpyde.inp

CENTRIF, BX, BY, BZ, GRAV, P1, P2, P3, P4, P5, P6, HP1, HP2, HP3, HP4, HP5, HP6, Q1, Q2, Q3, Q4, Q5, Q6, Q1D, Q2D, Q3D, Q4D, Q5D, Q6D, S1, S2, S3, S4, S5, S6.

ec3kpydl.inp

P, HP, Q, QD, S.

IV. CAXA elements with pore pressure

Problem description

Model:

Planar dimension3 × 5
Inside radius1.0
Gravity load vector(1, 1, 0)

Material:

Young's modulus1 × 108
Poisson's ratio0.0
Density1.4142
Permeability1 × 10–5
Specific weight of fluid1.0

Initial conditions:

Initial void ratio1.0
Hydrostatic pressure datum5.0
Hydrostatic pressure elevation0.0
Sink pore pressure14.7

Results and discussion

The calculated reactions are in agreement with the applied loads.

Input files

ecnwpfde.inp

CAXA8P1: BX, BZ, GRAV, P1, P2, P3, P4, HP1, HP2, HP3, HP4, Q1, Q2, Q3, Q4, Q1D, Q2D, Q3D, Q4D, S1, S2, S3, S4.

ecnwprde.inp

CAXA8RP1: BX, BZ, GRAV, P1, P2, P3, P4, HP1, HP2, HP3, HP4, Q1, Q2, Q3, Q4, Q1D, Q2D, Q3D, Q4D, S1, S2, S3, S4.

ecnxpfde.inp

CAXA8P2: BX, BZ, GRAV, P1, P2, P3, P4, HP1, HP2, HP3, HP4, Q1, Q2, Q3, Q4, Q1D, Q2D, Q3D, Q4D, S1, S2, S3, S4.

ecnxprde.inp

CAXA8RP2: BX, BZ, GRAV, P1, P2, P3, P4, HP1, HP2, HP3, HP4, Q1, Q2, Q3, Q4, Q1D, Q2D, Q3D, Q4D, S1, S2, S3, S4.

ecnypfde.inp

CAXA8P3: BX, BZ, GRAV, P1, P2, P3, P4, HP1, HP2, HP3, HP4, Q1, Q2, Q3, Q4, Q1D, Q2D, Q3D, Q4D, S1, S2, S3, S4.

ecnyprde.inp

CAXA8RP3: BX, BZ, GRAV, P1, P2, P3, P4, HP1, HP2, HP3, HP4, Q1, Q2, Q3, Q4, Q1D, Q2D, Q3D, Q4D, S1, S2, S3, S4.

ecnzpfde.inp

CAXA8P4: BX, BZ, GRAV, P1, P2, P3, P4, HP1, HP2, HP3, HP4, Q1, Q2, Q3, Q4, Q1D, Q2D, Q3D, Q4D, S1, S2, S3, S4.

ecnzprde.inp

CAXA8RP4: BX, BZ, GRAV, P1, P2, P3, P4, HP1, HP2, HP3, HP4, Q1, Q2, Q3, Q4, Q1D, Q2D, Q3D, Q4D, S1, S2, S3, S4.

V. Plane strain pore-thermal elements

Problem description

Model:


Planar dimension3 × 5
Gravity direction(1, 1, 0)

Material:


Young's modulus1 × 108
Poisson's ratio0.0
Density1.4142
Permeability1 × 10-5
Thermal conductivity0.1

Initial conditions:


Initial void ratio1.0
Initial temperature0.0
Initial pore pressure0.0

Results and discussion

The calculated reactions are in agreement with the applied loads.

Input files

cpe4pt_mech_loads.inp

CPE4PT: BX, BY, GRAV, CENTRIF, P1, P2, P3, P4, HP1, HP2, HP3, HP4.

cpe4rpt_mech_loads.inp

CPE4RPT: BX, BY, GRAV, CENTRIF, P1, P2, P3, P4, HP1, HP2, HP3, HP4.

cpe4pht_mech_loads.inp

CPE4PHT: BX, BY, GRAV, CENTRIF, P1, P2, P3, P4, HP1, HP2, HP3, HP4.

cpe4rpht_mech_loads.inp

CPE4RPHT: BX, BY, GRAV, CENTRIF, P1, P2, P3, P4, HP1, HP2, HP3, HP4.

VI. Plane strain pore-thermal elements with flow loads

Problem description

Model:


Planar dimension3 × 5

Material:


Young's modulus1 × 108
Poisson's ratio0.0
Density1.4142
Permeability1 × 10-5
Thermal conductivity0.1
Specific heat0.39

Initial conditions:


Initial void ratio1.0
Initial temperature0.0
Initial pore pressure0.0

Results and discussion

The calculated reactions are in agreement with the applied loads.

Input files

cpe4pt_flow_loads.inp

CPE4PT: Q1, Q2, Q3, Q4, S1, S2, S3, S4.

cpe4rpt_flow_loads.inp

CPE4RPT: Q1, Q2, Q3, Q4, S1, S2, S3, S4.

cpe4pht_flow_loads.inp

CPE4PHT: Q1, Q2, Q3, Q4, S1, S2, S3, S4.

cpe4rpht_flow_loads.inp

CPE4RPHT: Q1, Q2, Q3, Q4, S1, S2, S3, S4.

VII. Plane strain pore-thermal elements with heat loads

Problem description

Model:


Planar dimension7 × 7

Material:


Young's modulus30 × 106
Poisson's ratio0.3
Permeability1 × 10-5
Coefficient of thermal expansion0.0
Thermal conductivity3.77 × 10-5
Density82.9
Specific heat0.39
Coefficient of thermal expansion, pore fluid0.0
Thermal conductivity, pore fluid3.77 × 10-5
Density, pore fluid82.9
Specific heat, pore fluid0.39

Initial conditions:


Initial void ratio1.0
Initial temperature0.0
Initial pore pressure0.0

Results and discussion

The calculated reactions are in agreement with the applied loads.

Input files

cpe4pt_heat_loads.inp

CPE4PT: BF, S1, S2, S3, S4, F1, F2, F3, F4, R1, R2, R3, R4.

cpe4rpt_heat_loads.inp

CPE4RPT: BF, S1, S2, S3, S4, F1, F2, F3, F4, R1, R2, R3, R4.

cpe4pht_heat_loads.inp

CPE4PHT: BF, S1, S2, S3, S4, F1, F2, F3, F4, R1, R2, R3, R4.

cpe4rpht_heat_loads.inp

CPE4RPHT: BF, S1, S2, S3, S4, F1, F2, F3, F4, R1, R2, R3, R4.

VIII. Three-dimensional pore-thermal elements

Problem description

Model:


Cubic dimension7 × 7 × 7
Gravity direction(1, 0, 0)

Material:


Modulus3 × 106
Density10.0
Expansion0.0001
Specific heat1.0
Conductivity0.1
Density, pore fluid10.0
Expansion, pore fluid0.0001
Specific heat, pore fluid1.0
Conductivity, pore fluid0.1
Permeability0.01
Specific weight of fluid1.0

Initial conditions:


Initial void ratio1.0
Initial temperature0.0
Initial pore pressure0.0

Results and discussion

The calculated reactions are in agreement with the applied loads.

Input files

c3d4pt_loads.inp

C3D4PT: BX, BY, BZ, GRAV, CENT, CENTRIF, P1, P2, P3, P4, P5, P6, HP1, HP2, HP3, HP4, HP5, HP6.

c3d8pt_loads.inp

C3D8PT: BX, BY, BZ, GRAV, CENT, CENTRIF, P1, P2, P3, P4, P5, P6, HP1, HP2, HP3, HP4, HP5, HP6.

c3d8rpt_loads.inp

C3D8RPT: BX, BY, BZ, GRAV, CENT, CENTRIF, P1, P2, P3, P4, P5, P6, HP1, HP2, HP3, HP4, HP5, HP6.

c3d8rpht_loads.inp

C3D8RPHT: BX, BY, BZ, GRAV, CENT, CENTRIF, P1, P2, P3, P4, P5, P6, HP1, HP2, HP3, HP4, HP5, HP6.

c3d10pt_loads.inp

C3D10PT: BX, BY, BZ, GRAV, CENT, CENTRIF, P1, P2, P3, P4, HP1, HP2, HP3, HP4.

c3d10mpt_loads.inp

C3D10MPT: BX, BY, BZ, GRAV, CENT, CENTRIF, P1, P2, P3, P4, HP1, HP2, HP3, HP4.

IX. Three-dimensional pore-thermal elements with flow loads

Problem description

Model:


Cubic dimension5 × 3 × 1
Gravity direction(1, 1, 1)

Material:


Modulus1 × 108
Density1.7321
Expansion0.0
Specific heat10.0
Conductivity1.0
Density, pore fluid1.7321
Expansion, pore fluid0.0
Specific heat, pore fluid10.0
Conductivity, pore fluid1.0
Permeability1 × 10–5
Specific weight of fluid1.0

Initial conditions:


Initial void's ratio1.0
Initial temperature0.0
Initial pore pressure0.0

Results and discussion

The calculated reactions are in agreement with the applied loads.

Input files

c3d4ph_dsflow_loads.inp

C3D4PH: P, HP, Q, S.

c3d4pt_dsflow_loads.inp

C3D4PT: P, HP, Q, S.

c3d8pt_dflow_loads.inp

C3D8PT: BX, BY, BZ, GRAV, CENT, CENTRIF, P1, P2, P3, P4, P5, P6, HP1, HP2, HP3, HP4, HP5, HP6, Q1, Q2, Q3, Q4, Q5, Q6, S1, S2, S3, S4, S5, S6.

c3d8pt_sflow_loads.inp

C3D8PT: P, HP, Q, S.

c3d8rpt_dflow_loads.inp

C3D8RPT: BX, BY, BZ, GRAV, CENT, CENTRIF, P1, P2, P3, P4, P5, P6, HP1, HP2, HP3, HP4, HP5, HP6, Q1, Q2, Q3, Q4, Q5, Q6, S1, S2, S3, S4, S5, S6.

c3d8rpht_dflow_loads.inp

C3D8RPHT: BX, BY, BZ, GRAV, CENT, CENTRIF, P1, P2, P3, P4, P5, P6, HP1, HP2, HP3, HP4, HP5, HP6, Q1, Q2, Q3, Q4, Q5, Q6, S1, S2, S3, S4, S5, S6.

c3d10mpt_flow_loads.inp

C3D10MPT: BX, BY, BZ, GRAV, CENT, CENTRIF, P1, P2, P3, P4, HP1, HP2, HP3, HP4, Q1, Q2, Q3, Q4, S1, S2, S3, S4.

c3d10p_dsflow_loads.inp

C3D10P: P, HPQ, S.

c3d10pt_dsflow_loads.inp

C3D10PT: P, HPQ, S.

c3d10mpt_dsflow_loads.inp

C3D10MPT: P, HPQ, S.

X. Three-dimensional pore-thermal elements with heat loads

Problem description

Model:

For this set of verification problems both the solid and the pore fluid used identical heat transfer properties so that results could be compared with conventional heat transfer elements.


Cubic dimension7 × 7 × 7

Material:


Modulus3 × 106
Density82.9
Expansion0.0
Specific heat0.39
Conductivity3.77 × 10–5
Density, pore fluid82.9
Expansion, pore fluid0.0
Specific heat, pore fluid0.39
Conductivity, pore fluid3.77 × 10–5
Permeability0.001
Specific weight of fluid10.0

Initial conditions:


Initial void's ratio1.0
Initial temperature0.0
Initial pore pressure0.0

Results and discussion

The calculated reactions are in agreement with the applied loads.

Input files

c3d4pt_heat_loads.inp

C3D4PT: S, F, R.

c3d8pt_heat_loads.inp

C3D8PT: S, F, R.

c3d8rpt_heat_loads.inp

C3D8RPT: S, F, R.

c3d8rpht_heat_loads.inp

C3D8RPHT: S, F, R.

c3d10pt_heat_loads.inp

C3D10PT: S, F, R.

c3d10mpt_heat_loads.inp

C3D10MPT: S, F, R.

XI. Axisymmetric pore-thermal elements

Problem description

Model:


Cubic dimension3 × 3

Material:


Modulus3 × 106
Density5 × 10–5
Expansion0.0001
Specific heat1.0
Conductivity0.1
Density, pore fluid5 × 105
Expansion, pore fluid0.0001
Specific heat, pore fluid1.0
Conductivity, pore fluid0.1
Permeability0.01
Specific weight of fluid1.0

Initial conditions:


Initial void's ratio1.0
Initial temperature0.0
Initial pore pressure0.0

Results and discussion

The calculated reactions are in agreement with the applied loads.

Input files

cax4pt_loads.inp

CAX4PT: CENT, CENTRIF, BR, GRAV,HP1, HP2, HP3, HP4, P1, P2, P3, P4.

cax4rpt_loads.inp

CAX4RPT: CENT, CENTRIF, BR, GRAV,HP1, HP2, HP3, HP4, P1, P2, P3, P4.

cax4rpht_loads.inp

CAX4RPHT: CENT, CENTRIF, BR, GRAV,HP1, HP2, HP3, HP4, P1, P2, P3, P4.

XII. Axisymmetric pore-thermal elements with flow loads

Problem description

Model:


Cubic dimension3 × 5

Material:


Modulus1 × 108
Density1.4142
Expansion0.0
Specific heat10.0
Conductivity1.0
Density, pore fluid1.4142
Expansion, pore fluid0.0
Specific heat, pore fluid10.0
Conductivity, pore fluid1.0
Permeability1 × 10–5
Specific weight of fluid1.0

Initial conditions:


Initial void's ratio1.0
Initial temperature0.0
Initial pore pressure0.0

Results and discussion

The calculated reactions are in agreement with the applied loads.

Input files

cax4pt_dflow_loads.inp

CAX4PT: CENTRIF, BR,HP1, HP2, HP3, HP4, Q1, Q2, Q3, Q4, S1, S2, S3, S4.

cax4rpt_dflow_loads.inp

CAX4RPT: CENTRIF, BR,HP1, HP2, HP3, HP4, Q1, Q2, Q3, Q4, S1, S2, S3, S4.

cax4rpht_dflow_loads.inp

CAX4RPHT: CENTRIF, BR,HP1, HP2, HP3, HP4, Q1, Q2, Q3, Q4, S1, S2, S3, S4.

XIII. Axisymmetric pore-thermal elements with heat loads

Problem description

Model:

For this set of verification problems both the solid and the pore fluid used identical heat transfer properties so that results could be compared with conventional heat transfer elements.


Cubic dimension7 × 7

Material:


Modulus30 × 106
Density82.9
Expansion0.0
Specific heat0.39
Conductivity3.77 × 10–5
Density, pore fluid82.9
Expansion, pore fluid0.0
Specific heat, pore fluid0.39
Conductivity, pore fluid3.77 × 10–5
Permeability0.001
Specific weight of fluid10.0

Initial conditions:


Initial void's ratio1.0
Initial temperature0.0
Initial pore pressure0.0

Results and discussion

The calculated reactions are in agreement with the applied loads.

Input files

cax4pt_heat_loads.inp

CAX4PT: S, F, R.

cax4rpt_heat_loads.inp

CAX4RPT: S, F, R.

cax4rpht_heat_loads.inp

CAX4RPHT: S, F, R.