CPT212
2D 4-Node Coupled Pore-Pressure-Thermal Mechanical Solid
CPT212 Element Description
CPT212 is a 2D four-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 four nodes and can have the following degrees of freedom at each node:
Translations in the nodal x and y directions
Pore-pressure (PRES)
Temperature (TEMP)
Nonlocal field values (GFV1, GFV2, GFV3)
CPT212 can be used as a plane strain or axisymmetric element. The element has stress stiffening, large deflection, and large strain capabilities. Various output options are also available.
See CPT212 for more details about this element.
CPT212 Input Data
The geometry and node locations for this element are shown in Figure 212.1: CPT212 Geometry.
A degenerated triangular-shaped element can be formed by defining the same node number for nodes K and L. 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 212.1: CPT212 Geometry. 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).
CPT212 surface, body, and nodal-force loads are given in Table 212.1: CPT212 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 212.1: CPT212 Surface, Body, and nodal-force loads
Coupled-Field Analysis | KEYOPT | Load Type | Load | Command Label |
---|---|---|---|---|
Structural-thermal | KEYOPT(11) = 1 | Surface | Structural surface pressure | PRES |
Heat flux | HFLUX | |||
Body | Heat generation | HGEN | ||
Nodal Force | Heat flow | HEAT | ||
Structural-pore-fluid-diffusion | KEYOPT(12) = 1 | Surface | Structural surface pressure | PRES |
Surface flow flux | FFLX | |||
Body | Flow source | FSOU | ||
Temperature | TEMP | |||
Nodal Force | Fluid flow | FLOW | ||
Structural-pore-fluid-diffusion-thermal | KEYOPT(11) = 1 and KEYOPT(12) = 1 | Surface | Structural surface pressure | PRES |
Heat flux | HFLUX | |||
Surface flow flux | FFLX | |||
Body | Heat generation | HGEN | ||
Flow source | FSOU | |||
Nodal Force | Heat flow | HEAT | ||
Fluid flow | FLOW | |||
Structural implicit gradient regularization | KEYOPT(18) = 1, 2, or 3 | Surface | Structural surface pressure | PRES |
Body | Temperature | TEMP |
Input the nodal forces, if any, per unit of depth for a plane analysis and on a full 360° basis for an axisymmetric analysis.
As described in Coordinate Systems, you can use ESYS to orient the material properties and strain/stress output. Issue RSYS 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).
The following table summarizes the element input. For a general description of element input, see Element Input.
CPT212 Input Summary
- Nodes
I, J, K, L
- Degrees of Freedom
UX, UY, 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(3)
Element behavior:
- 1 --
Axisymmetric
- 2 --
Plane strain (Z strain = 0.0) (default)
- 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)
CPT212 Technology
CPT212 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.
CPT212 Output Data
The solution output associated with the element is in two forms:
Nodal degrees of freedom included in the overall nodal solution
Additional element output as shown in Table 212.2: CPT212 Element Output Definitions
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 212.2: CPT212 Stress Output. A general description of solution output is given in Solution Output. See the Basic Analysis Guide for ways to view results.
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 212.2: CPT212 Element Output Definitions
Name | Definition | O | R |
---|---|---|---|
ALL ANALYSES | |||
EL | Element number | - | Y |
NODES | Nodes - I, J, K, L | - | Y |
MAT | Material number | - | Y |
THICK | Thickness | - | Y |
VOLU | Volume | - | Y |
XC, YC | Location where results are reported | Y | 1 |
ALL ANALYSES WITH A STRUCTURAL FIELD | |||
S:X, Y, Z, XY | Stresses | Y | Y |
S:1, 2, 3 | Principal stresses | - | Y |
S:INT | Stress intensity | - | Y |
S:EQV | Equivalent stress | Y | Y |
EPEL:X, Y, Z, XY | Elastic strains | Y | Y |
EPEL:1, 2, 3 | Principal elastic strains | - | Y |
EPEL:EQV | Equivalent elastic strain [2] | Y | Y |
EPTH:X, Y, Z, XY | Thermal strains | Y | Y |
EPTH:EQV | Equivalent thermal strain [2] | - | Y |
EPPL:X, Y, Z, XY | Plastic strains | - | Y |
EPPL:EQV | Equivalent plastic strain [2] | - | Y |
EPTO:X, Y, Z, XY | Total mechanical strains (EPEL + EPPL) | - | Y |
EPTO:EQV | Total equivalent mechanical strain (EPEL + EPPL) | - | Y |
TEMP | Temperatures T(I), T(J), T(K), T(L) | - | Y |
ADDITIONAL OUTPUT FOR ANALYSES WITH A TEMPERATURE FIELD | |||
TG:X, Y | Thermal gradient components | - | Y |
TF:X, Y | Thermal flux components | - | Y |
ADDITIONAL OUTPUT FOR ANALYSES WITH A PORE-PRESSURE FIELD | |||
ESIG:X, Y, Z, XY | Effective stresses | - | Y |
FGRA:X, Y | Fluid pore-pressure gradient components | - | Y |
FFLX:X, Y | Fluid flow flux components | - | Y |
PMSV:VRAT,PPRE,DSAT,RPER | Void volume ratio, pore pressure, degree of saturation, and relative permeability | - | Y |
EPFR | Free strain | - | Y |
ADDITIONAL OUTPUT FOR ANALYSES WITH A NONLOCAL FIELD | |||
MPDP:TOTA,TENS,COMP,RW | Microplane homogenized total, tension, and compression damages (TOTA, TENS, COMP), and split weight factor (RW). | - | Y |
DAMAGE: 1,2,3,MAX | Damage in directions 1, 2, 3 (1, 2, 3) and the maximum damage (MAX). | - | Y |
GMDG | Damage | - | Y |
IDIS | Structural-thermal dissipation rate | - | Y |
For axisymmetric solutions in a global coordinate system, the X, Y, XY, and Z stress and strain outputs correspond to the radial, axial, in-plane shear, and hoop stresses and strains, respectively.
Table 212.3: CPT212 Item and Sequence Numbers lists output available via the ETABLE command using the Sequence Number method. For more information, see Creating an Element Table and The Item and Sequence Number Table in this document. The table uses the following notation:
- Name
output quantity as defined in the Table 212.2: CPT212 Element Output Definitions
- Item
predetermined Item label for ETABLE
- E
sequence number for single-valued or constant element data
- I,J,K,L
sequence number for data at nodes I, J, K, L
CPT212 Assumptions and Restrictions
The area of the element must be positive.
The element must lie in a global X-Y plane as shown in Figure 212.1: CPT212 Geometry and the Y-axis must be the axis of symmetry for axisymmetric analyses. An axisymmetric structure should be modeled in the +X quadrants.
You can form a triangular element by defining duplicate K and L node numbers. (For more information, see Degenerated Shape Elements.)
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.