PLANE83
Axisymmetric-Harmonic
8-Node Structural Solid
PLANE83 Element Description
Unless you are performing a linear analysis requiring a large number of circumferential modes, consider using a current-technology element such as SOLID273 (KEYOPT(6) = 0). |
PLANE83 is used for 2D modeling of axisymmetric structures with nonaxisymmetric loading. Examples of such loading are bending, shear, or torsion. The element has three degrees of freedom per node: translations in the nodal x, y, and z directions. For unrotated nodal coordinates, these directions correspond to the radial, axial, and tangential directions, respectively.
This element is a higher order version of the 2D, four-node element (PLANE25). It provides more accurate results for mixed (quadrilateral-triangular) automatic meshes and can tolerate irregular shapes without as much loss of accuracy. The loading need not be axisymmetric. Various loading cases are described in Harmonic Axisymmetric Elements with Nonaxisymmetric Loads.
The 8-node elements have compatible displacement shapes and are well suited to model curved boundaries. See PLANE83 in the Mechanical APDL Theory Reference for more details about this element.
PLANE83 Input Data
The geometry, node locations, and the coordinate system for this element are shown in Figure 83.1: PLANE83 Geometry. Orthotropic material directions correspond to the element coordinate directions. The element coordinate system orientation is as described in Coordinate Systems. Z-direction material properties (EZ, ALPZ, etc.) may be input. MODE and ISYM are used to describe the harmonic loading condition. (See Harmonic Axisymmetric Elements with Nonaxisymmetric Loads for more details.)
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 83.1: PLANE83 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). If all corner node temperatures are specified, each midside node temperature defaults to the average temperature of its adjacent corner nodes. For any other input temperature pattern, unspecified temperatures default to TUNIF.
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 "PLANE83 Input Summary". A general description of element input is given in Element Input.
PLANE83 Input Summary
- Nodes
I, J, K, L, M, N, O, P
- 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, ALPD, BETD, 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), T(M), T(N), T(O), T(P)
- 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(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 (not extra output)
- 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)
Extra surface output (surface solution is valid only for isotropic materials):
- 0 --
Basic element solution (no extra output)
- 1 --
Surface solution for face I-J also
- 2 --
Surface solution for both faces I-J and K-L also
PLANE83 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 83.1: PLANE83 Element Output Definitions
Several items are illustrated in Figure 83.2: PLANE83 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.). 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 element stress directions are parallel to the element coordinate system. 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. 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 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 83.1: PLANE83 Element Output Definitions
Name | Definition | O | R |
---|---|---|---|
EL | Element Number | Y | Y |
NODES | Corner 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 | Y |
VOLU: | Volume | Y | Y |
PRES | Pressures 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), T(M), T(N), T(O), T(P) | Y | Y |
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 | Y | Y |
EPTH:EQV | Equivalent thermal strain [4] | - | Y |
PK ANG | Angle where stresses have peak values: 0 and 90/MODE°. Blank if MODE = 0. | Y | Y |
XC, YC | Location where results are reported | Y | 3 |
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 83.2: PLANE83 Item and Sequence Numbers lists output available through the ETABLE command using the Sequence Number method. See Creating an Element Table in the Basic Analysis Guide and The Item and Sequence Number Table in this reference for more information. The following notation is used in Table 83.2: PLANE83 Item and Sequence Numbers:
- Name
output quantity as defined in the Table 83.1: PLANE83 Element Output Definitions
- Item
predetermined Item label for ETABLE command
- I,J,K,L
sequence number for data at nodes I,J,K,L
Table 83.2: PLANE83 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 |
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.
PLANE83 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 83.1: PLANE83 Geometry and the Y-axis must be the axis of symmetry for axisymmetric analyses. An axisymmetric structure should be modeled in the +X quadrants.
A face with a removed midside node implies that the displacement varies linearly, rather than parabolically, along that face. See Quadratic Elements (Midside Nodes) in the Modeling and Meshing Guide for more information about the use of midside nodes.
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.
The element may not be deactivated with the EKILL command.
The element temperature is taken to be the average of the nodal temperatures.
Surface stress printout is valid only if the conditions described in Element Solution are met.
Modeling hints: —
If shear effects are important in a shell-like structure, at least two elements through the thickness should be used.
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.