CPT215


3D 8-Node Coupled Pore-Pressure-Thermal Mechanical Solid

Valid Products: Pro | Premium | Enterprise | PrepPost | Solver | AS add-on

CPT215 Element Description

CPT215 is a 3D eight-node coupled physics solid element capable of modeling coupled physics phenomena such as structural-pore-fluid-diffusion-thermal analysis and structural implicit gradient regularization using a nonlocal field. The element is defined by eight nodes and can have the following degrees of freedom at each node:

  • Translations in the nodal x, y, and z directions

  • Pore-pressure (PRES)

  • Temperature (TEMP)

  • Nonlocal field values (GFV1, GFV2, GFV3)

CPT215 has elasticity, stress stiffening, large deflection, and large strain capabilities. Various printout options are available.

For more details about this element, see CPT215 .

Figure 215.1: CPT215 Structural Solid Geometry

CPT215 Structural Solid Geometry

A higher-order version of this element is CPT216.

CPT215 Input Data

The geometry and node locations for this element are shown in Figure 215.1: CPT215 Structural Solid Geometry. A prism-shaped element can be formed by defining the same node numbers for nodes K and L, and nodes O and P. A tetrahedral-shaped element and a pyramid-shaped element can also be formed, as shown in the illustration. (CPT217 is a similar element, but is a 10-node tetrahedron.)

In addition to the nodes, for structural-pore-fluid-diffusion-thermal analysis, the element input data includes the orthotropic material properties. Orthotropic material directions correspond to the element coordinate directions. The element coordinate system orientation is described in Coordinate Systems.

Element loads are described in Element Loading. Loads can be input (SF and SFE) on the element faces indicated by the circled numbers in Figure 215.1: CPT215 Structural Solid Geometry. Positive pressures act into the element. Positive pressures act into the element. Body loads may be input (BF and BFE) at the element nodes or as a single element value. Nodal forces can be applied to the nodes directly (F).

CPT215 surface, body, and nodal-force loads are given in Table 215.1: CPT215 Surface, Body, and nodal-force loads. Also see Loading Types in the Coupled-Field Analysis Guide.

Most loads can be defined as a function of primary variables by using tabular input. For more information, see Applying Loads Using Tabular Input in the Basic Analysis Guide and the descriptions of individual loading commands in the Command Reference.

Table 215.1: CPT215 Surface, Body, and nodal-force loads

Coupled-Field AnalysisKEYOPTLoad TypeLoadCommand Label
Structural-thermalKEYOPT(11) = 1SurfaceStructural surface pressurePRES
Heat fluxHFLUX
BodyHeat generationHGEN
Nodal ForceHeat flowHEAT
Structural-pore-fluid-diffusionKEYOPT(12) = 1SurfaceStructural surface pressurePRES
Surface flow fluxFFLX
BodyFlow sourceFSOU
TemperatureTEMP
Nodal ForceFluid flowFLOW
Structural-pore-fluid-diffusion-thermalKEYOPT(11) = 1 and KEYOPT(12) = 1SurfaceStructural surface pressurePRES
Heat fluxHFLUX
Surface flow fluxFFLX
BodyHeat generationHGEN
Flow sourceFSOU
Nodal ForceHeat flowHEAT
Fluid flowFLOW
Structural implicit gradient regularizationKEYOPT(18) = 1, 2, or 3SurfaceStructural surface pressurePRES
BodyTemperatureTEMP

For problems that do not consider the optional temperature degrees of freedom, temperatures can be input as element body loads at the nodes. The node I temperature T(I) defaults to TUNIF. If all other temperatures are unspecified, they default to T(I). For any other input temperature pattern, unspecified temperatures default to TUNIF.

As described in Coordinate Systems, you can use the ESYS command to orient the material properties and strain/stress output. Use the RSYS command to choose output that follows the material coordinate system or the global coordinate system.

The element generally produces an unsymmetric matrix. To avoid convergence difficulty, use the unsymmetric solver (NROPT,UNSYM).

"CPT215 Input Summary" contains a summary of element input. For a general description of element input, see Element Input.

CPT215 Input Summary

Nodes

I, J, K, L, M, N, O, P

Degrees of Freedom

UX, UY, UZ, PRES, TEMP, GFV1, GFV2, GFV3

Real Constants
None
Material Properties
TB command: See Element Support for Material Models for this element.
MP command: EX, EY, EZ, ALPX, ALPY, ALPZ (or CTEX, CTEY, CTEZ or THSX, THSY, THSZ), PRXY, PRYZ, PRXZ (or NUXY, NUYZ, NUXZ), DENS, GXY, GYZ, GXZ, ALPD, BETD
Surface Loads
Body Loads
Special Features
KEYOPT(6)

Element formulation in coupled-field analyses with structural degrees of freedom:

0 -- 

Pure displacement formulation (default)

1 -- 

Mixed u-P formulation

KEYOPT(11)

Temperature degree of freedom:

0 -- 

Disabled (default)

1 -- 

Enabled

KEYOPT(12)

Pressure degree of freedom:

0 -- 

Disabled (default)

1 -- 

Enabled

KEYOPT(18)

Nonlocal degree of freedom:

0 -- 

Disabled (default)

1 -- 

Enabled (adds one extra degree of freedom per node)

2 -- 

Enabled (adds two extra degrees of freedom per node)

3 -- 

Enabled (adds three extra degrees of freedom per node)

CPT215 Technology

CPT215 uses the method (also known as selective reduced integration). This approach helps to prevent volumetric mesh locking in nearly incompressible cases. It replaces volumetric strain at the Gauss integration point with the average volumetric strain of the elements.

CPT215 Output Data

The solution output associated with the element is in two forms:

By default, the integration point results are copied to the nodes (ERESX).

The element stress directions are parallel to the element coordinate system, as shown in Figure 215.2: CPT215 Stress Output. A general description of solution output is given in The Item and Sequence Number Table. See the Basic Analysis Guide for ways to view results.

Figure 215.2: CPT215 Stress Output

CPT215 Stress Output

The element stress directions are parallel to the global coordinate system.


The Element Output Definitions table uses the following notation:

A colon (:) in the Name column indicates that the item can be accessed by the Component Name method (ETABLE, ESOL). The O column indicates the availability of the items in the file jobname.out. The R column indicates the availability of the items in the results file.

In either the O or R columns, Y indicates that the item is always available, a number refers to a table footnote that describes when the item is conditionally available, and - indicates that the item is not available. All output is available only if calculated (based on input values).

Table 215.2: CPT215 Element Output Definitions

NameDefinitionOR
ALL ANALYSES
ELElement number-Y
NODESNodes - I, J, K, L -Y
MATMaterial number-Y
THICKThickness-Y
VOLUVolume-Y
XC, YCLocation where results are reportedY1
ALL ANALYSES WITH A STRUCTURAL FIELD
S:X, Y, Z, XYStressesYY
S:1, 2, 3Principal stresses-Y
S:INTStress intensity-Y
S:EQVEquivalent stressYY
EPEL:X, Y, Z, XYElastic strainsYY
EPEL:1, 2, 3Principal elastic strains-Y
EPEL:EQVEquivalent elastic strain [2]YY
EPTH:X, Y, Z, XYThermal strainsYY
EPTH:EQVEquivalent thermal strain [2]-Y
EPPL:X, Y, Z, XYPlastic strains-Y
EPPL:EQVEquivalent plastic strain [2]-Y
EPTO:X, Y, Z, XYTotal mechanical strains (EPEL + EPPL)-Y
EPTO:EQVTotal equivalent mechanical strain (EPEL + EPPL)-Y
TEMPTemperatures T(I), T(J), T(K), T(L)-Y
ADDITIONAL OUTPUT FOR ANALYSES WITH A TEMPERATURE FIELD
TG:X, YThermal gradient components-Y
TF:X, YThermal flux components-Y
ADDITIONAL OUTPUT FOR ANALYSES WITH A PORE-PRESSURE FIELD
ESIG:X, Y, Z, XYEffective stresses-Y
FGRA:X, YFluid pore-pressure gradient components-Y
FFLX:X, YFluid flow flux components-Y
PMSV:VRAT,PPRE,DSAT,RPERVoid volume ratio, pore pressure, degree of saturation, and relative permeability-Y
EPFRFree strain-Y
ADDITIONAL OUTPUT FOR ANALYSES WITH A NONLOCAL FIELD
MPDP:TOTA,TENS,COMP,RWMicroplane homogenized total, tension, and compression damages (TOTA, TENS, COMP), and split weight factor (RW).-Y
DAMAGE: 1,2,3,MAXDamage in directions 1, 2, 3 (1, 2, 3) and the maximum damage (MAX).-Y
GMDGDamage-Y
IDISStructural-thermal dissipation rate-Y

  1. Available only at centroid as a *GET item

  2. The equivalent strains use an effective Poisson's ratio: for elastic and thermal this value is set by the user (MP,PRXY); for plastic and creep this value is set at 0.5.

Table 215.3: CPT215 Item and Sequence Numbers lists output available via ETABLE using the Sequence Number method. See Element Table for Variables Identified By Sequence Number and The Item and Sequence Number Table in this document for more information. The following notation is used in Table 215.3: CPT215 Item and Sequence Numbers:

Name

output quantity as defined in the Table 215.2: CPT215 Element Output Definitions

Item

predetermined Item label for ETABLE command

I,J,...,P

sequence number for data at nodes I, J, ..., P

Table 215.3: CPT215 Item and Sequence Numbers

Output Quantity Name ETABLE and ESOL Command Input
ItemIJKLMNOP
P1SMISC2143----
P2SMISC56--87--
P3SMISC-910--1211-
P4SMISC--1314--1615
P5SMISC18--1719--20
P6SMISC----21222324

CPT215 Assumptions and Restrictions

  • The element must not have a zero volume. Also, the element may not be twisted such that the element has two separate volumes (which occurs most frequently when the element is numbered improperly). Elements may be numbered either as shown in Figure 215.1: CPT215 Structural Solid Geometry or may have the planes IJKL and MNOP interchanged.

  • When degenerated into a tetrahedron, wedge, or pyramid element shape (described in Degenerated Shape Elements), the corresponding degenerated shape functions are used. Degeneration to a pyramidal form should be used with caution. The element sizes, when degenerated, should be small to minimize the stress gradients. Pyramid elements are best used as filler elements or in meshing transition zones.

  • Stress stiffening is always included in geometrically nonlinear analyses (NLGEOM,ON). It is ignored in geometrically linear analyses (NLGEOM,OFF). Prestress effects can be activated by the PSTRES command.

CPT215 Product Restrictions

There are no product-specific restrictions for this element.