PLANE25
Axisymmetric-Harmonic
4-Node Structural Solid
PLANE25 Element Description
Unless you are performing a linear analysis requiring a large number of circumferential modes, consider using a current-technology element such as SOLID272 (KEYOPT(6) = 0). |
PLANE25 is used for 2D modeling of axisymmetric structures with nonaxisymmetric loading. Examples of such loading are bending, shear, or torsion. The element is defined by four nodes having three degrees of freedom per node: translations in the nodal x, y, and z direction. For unrotated nodal coordinates, these directions correspond to the radial, axial, and tangential directions, respectively.
See Harmonic Axisymmetric Elements with Nonaxisymmetric Loads for a description of various loading cases. See PLANE25 in the Mechanical APDL Theory Reference for more details about this element. See PLANE83 for a multi-node version of this element.
PLANE25 Input Data
The geometry, node locations, and the coordinate system for this element are shown in Figure 25.1: PLANE25 Geometry. The element input data includes four nodes, the number of harmonic waves (MODE on the MODE command), the symmetry condition (ISYM on the MODE command) and the orthotropic material properties. The MODE and ISYM parameters are discussed in detail in Harmonic Axisymmetric Elements with Nonaxisymmetric Loads.
The material may be orthotropic, with directions corresponding to the element coordinate directions. The element coordinate system orientation is as described in Coordinate Systems. Harmonically varying nodal forces, if any, should be input on a full 360° basis.
Element loads are described in Element Loading. Harmonically varying pressures may be input as surface loads on the element faces as shown by the circled numbers on Figure 25.1: PLANE25 Geometry. Positive pressures act into the element.
Harmonically varying temperatures may 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 pattern, unspecified temperatures default to TUNIF.
KEYOPT(2) is used to include or suppress the extra displacement shapes.
KEYOPT(3) is used for temperature loading with MODE > 0 and temperature-dependent material properties. Material properties may only be evaluated at a constant (nonharmonically varying) temperature. If MODE = 0, the material properties are always evaluated at the average element temperature. If MODE > 0, TREF must be input as zero.
KEYOPT(4), (5), and (6) provide various element printout options (see Element Solution).
A summary of the element input is given in "PLANE25 Input Summary". Element Input gives a general description of element input.
PLANE25 Input Summary
- Nodes
I, J, K, L
- Degrees of Freedom
UX, UY, UZ
- Real Constants
None
- Material Properties
MP command: EX, EY, EZ, PRXY, PRYZ, PRXZ (or NUXY, NUYZ, NUXZ), ALPX, ALPY, ALPZ (or CTEX, CTEY, CTEZ or THSX, THSY, THSZ), DENS, GXY, GYZ, GXZ, BETD, ALPD, DMPR, DMPS - Surface Loads
- Pressures --
face 1 (J-I), face 2 (K-J), face 3 (L-K), face 4 (I-L)
- Body Loads
- Temperatures --
T(I), T(J), T(K), T(L)
- Mode Number
Number of harmonic waves around the circumference (MODE)
- Loading Condition
Symmetry condition (MODE)
- Special Features
- KEYOPT(1)
Element coordinate system:
- 0 --
Element coordinate system is parallel to the global coordinate system
- 1 --
Element coordinate system is based on the element I-J side.
- KEYOPT(2)
Extra displacement shapes:
- 0 --
Include extra displacement shapes
- 1 --
Suppress extra displacement shapes
- KEYOPT(3)
If MODE is greater than zero, use temperatures for:
- 0 --
Use temperatures only for thermal bending (evaluate material properties at TREF)
- 1 --
Use temperatures only for material property evaluation (thermal strains are not computed)
- KEYOPT(4)
Extra stress output:
- 0 --
Basic element solution
- 1 --
Repeat basic solution for all integration points
- 2 --
Nodal Stress Solution
- KEYOPT(5)
Combined stress output:
- 0 --
No combined stress solution
- 1 --
Combined stress solution at centroid and nodes
- KEYOPT(6)
Include extra surface output (surface solution valid only for isotropic materials):
- 0 --
Basic element solution
- 1 --
Surface solution for face I-J also
- 2 --
Surface solution for both faces I-J and K-L also
PLANE25 Output Data
The solution output associated with the element is in two forms:
Nodal displacements included in the overall nodal solution
Additional element output as shown in Table 25.1: PLANE25 Element Output Definitions.
Several items are illustrated in Figure 25.2: PLANE25 Stress Output.
In the displacement printout, the UZ component is out-of-phase with the UX and UY components. For example, in the MODE = 1, ISYM = 1 loading case, UX and UY are the peak values at θ = 0° and UZ is the peak value at θ = 90°. The same occurs for the reaction forces (FX, FY, etc.). The element stress directions are parallel to the element coordinate system. We recommend that you always use the angle field on the SET command when postprocessing the results. For more information about harmonic elements, see Harmonic Axisymmetric Elements with Nonaxisymmetric Loads
The sign convention on the surface shears is such that for a rectangular element that is lined up parallel to the axes with node J in the positive Y direction from node I, the shear stresses on surfaces I-J and K-L are analogous to the centroidal SYZ in both definition and sign. Stress components which are inherently zero for a load case are printed for clarity. Solution Output gives a general description of 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 letter or number refers to a table footnote that describes when the item is conditionally available, and “-” indicates that the item is not available.
Table 25.1: PLANE25 Element Output Definitions
Name | Definition | O | R |
---|---|---|---|
EL | Element Number | Y | Y |
NODES | Nodes - I, J, K, L | Y | Y |
MAT | Material number | Y | Y |
ISYM | Loading key: 1 = symmetric, -1 = antisymmetric | Y | - |
MODE | Number of waves in loading | Y | - |
VOLU: | Volume | Y | Y |
PRES | Pressure P1 at nodes J,I; P2 at K,J; P3 at L,K; P4 at I,L | Y | Y |
TEMP | Temperatures T(I), T(J), T(K), T(L) | Y | Y |
PK ANG | Angle where component stresses have peak values: 0 and 90/MODE degrees. Blank if MODE = 0. | Y | Y |
XC, YC | Location where results are reported | Y | 3 |
S:X, Y, Z | Direct stresses (radial, axial, hoop) at PK ANG locations | Y | Y |
S:XY, YZ, XZ | Shear stresses (radial-axial, axial-hoop, radial-hoop) at PK ANG locations | Y | Y |
S:1, 2, 3 | Principal stresses at both PK ANG locations as well as where extreme occurs (EXTR); if MODE = 0, only one location is given. | 1 | 1 |
S:INT | Stress intensity at both PK ANG locations as well as where extreme occurs (EXTR); if MODE = 0, only one location is given. | 1 | 1 |
S:EQV | Equivalent stress at both PK ANG locations as well as where extreme occurs (EXTR); if MODE = 0, only one location is given. | 1 | 1 |
EPEL:X, Y, Z, XY | Elastic strain | Y | Y |
EPEL:EQV | Equivalent elastic strain [4] | - | Y |
EPTH:X, Y, Z, XY | Average thermal strains | 1 | 1 |
EPTH:EQV | Equivalent thermal strain [4] | - | 1 |
FACE | Face label | 2 | 2 |
TEMP | Surface average temperature | 2 | 2 |
EPEL(PAR, PER, Z, SH) | Surface strains (parallel, perpendicular, hoop, shear) at PK ANG locations and where extreme occurs (EXTR) | 2 | 2 |
S(PAR, PER, Z, SH) | Surface stresses (parallel, perpendicular, hoop, shear) at PK ANG locations and where extreme occurs (EXTR) | 2 | 2 |
Table 25.2: PLANE25 Item and Sequence Numbers lists output available through the ETABLE command using the Sequence Number method. See The General Postprocessor (POST1) of the Basic Analysis Guide and The Item and Sequence Number Table for more information. The following notation is used in Table 25.2: PLANE25 Item and Sequence Numbers:
- Name
output quantity as defined in the Table 25.1: PLANE25 Element Output Definitions
- Item
predetermined Item label for ETABLE command
- E
sequence number for single-valued or constant element data
- I,J,K,L
sequence number for data at nodes I,J,K,L
Table 25.2: PLANE25 Item and Sequence Numbers
Output Quantity Name | ETABLE and ESOL Command Input | ||||
---|---|---|---|---|---|
Item | I | J | K | L | |
P1 | SMISC | 2 | 1 | - | - |
P2 | SMISC | - | 4 | 3 | - |
P3 | SMISC | - | - | 6 | 5 |
P4 | SMISC | 7 | - | - | 8 |
THETA = 0 | |||||
S1 | NMISC | 1 | 16 | 31 | 46 |
S2 | NMISC | 2 | 17 | 32 | 47 |
S3 | NMISC | 3 | 18 | 33 | 48 |
SINT | NMISC | 4 | 19 | 34 | 49 |
SEQV | NMISC | 5 | 20 | 35 | 50 |
THETA = 90/MODE | |||||
S1 | NMISC | 6 | 21 | 36 | 51 |
S2 | NMISC | 7 | 22 | 37 | 52 |
S3 | NMISC | 8 | 23 | 38 | 53 |
SINT | NMISC | 9 | 24 | 39 | 54 |
SEQV | NMISC | 10 | 25 | 40 | 55 |
EXTR Values | |||||
S1 | NMISC | 11 | 26 | 41 | 56 |
S2 | NMISC | 12 | 27 | 42 | 57 |
S3 | NMISC | 13 | 28 | 43 | 58 |
SINT | NMISC | 14 | 29 | 44 | 59 |
SEQV | NMISC | 15 | 30 | 45 | 60 |
Note: The NMISC items (1 thru 60) in the above table represent the combined stress solution, KEYOPT(5) = 1. If MODE = 0, their values are zero at THETA = 90/MODE and at EXTR.
See Surface Solution for the item and sequence numbers for surface output for the ETABLE command.
PLANE25 Assumptions and Restrictions
The area of the element must be positive.
The element must be defined in the global X-Y plane as shown in Figure 25.1: PLANE25 Geometry and the global X-axis must be the radial direction. Negative X coordinates should not be used.
The element assumes a linear elastic material. Post-analysis superposition of results is valid only with other linear elastic solutions. The element should not be used with the large deflection option.
A triangular element may be formed by defining duplicate K and L node numbers (see Degenerated Shape Elements). The extra shapes are automatically deleted for triangular elements so that a constant strain element results.
Surface stress printout is valid only if the conditions described in Element Solution are met.
You can use only axisymmetric (MODE,0) loads without significant torsional stresses to generate the stress state used for stress stiffened modal analyses using this element.
This element does not support spectrum analysis.
Modeling hint: If shear effects are important in a shell-like structure, you should use at least two elements through the thickness.