PIPE16
Elastic
Straight Pipe
PIPE16 Element Description
Although this archived element is available for use in your analysis, Ansys, Inc. recommends using a current-technology element such as PIPE288 with KEYOPT(3) = 2. |
PIPE16 is a uniaxial element with tension-compression, torsion, and bending capabilities. The element has six degrees of freedom at two nodes: translations in the nodal x, y, and z directions and rotations about the nodal x, y, and z axes. See PIPE16 - Elastic Straight Pipe for more details about this element.
PIPE16 Input Data
The geometry, node locations, and the coordinate system for this element are shown in Figure 16.1: PIPE16 Geometry. The element input data include two or three nodes, the pipe outer diameter and wall thickness, stress intensification and flexibility factors, internal fluid density, exterior insulation density and thickness, corrosion thickness allowance, insulation surface area, pipe wall mass, axial pipe stiffness, rotordynamic spin, and the isotropic material properties.
The element X-axis is oriented from node I toward node J. For the two-node option, the element Y-axis is automatically calculated to be parallel to the global X-Y plane. Several orientations are shown in Figure 16.1: PIPE16 Geometry. For the case where the element is parallel to the global Z-axis (or within a 0.01 percent slope of it), the element Y-axis is oriented parallel to the global Y-axis (as shown). For user control of the element orientation about the element X-axis, use the third node option. The third node (K), if used, defines a plane (with I and J) containing the element X and Z axes (as shown). Input and output locations around the pipe circumference identified as being at 0° are located along the element Y-axis, and similarly 90° is along the element Z-axis.
The stress-intensification factor modifies the bending stress. Stress intensification factors can be input at end I (SIFI) and end J (SIFJ), if KEYOPT(2) = 0, or determined by the program using a tee-joint calculation if KEYOPT(2) = 1, 2, or 3. Stress-intensification factor values less than 1.0 are set equal to 1.0. The flexibility factor (FLEX) is divided into the cross-sectional moment of inertia to produce a modified moment of inertia for the bending stiffness calculation. FLEX defaults to 1.0 but may be input as any positive value.
The element mass is calculated from the pipe wall material, the external insulation, and the internal fluid. The insulation and the fluid contribute only to the element mass matrix. The corrosion thickness allowance contributes only to the stress calculations. A positive wall mass real constant overrides the pipe wall mass calculation. A nonzero insulation area real constant overrides the insulation surface area calculation (from the pipe outer diameter and length). A nonzero stiffness real constant overrides the calculated axial pipe stiffness.
Element loads are described in Element Loading. Pressures may be input as surface loads on the element faces as shown by the circled numbers on Figure 16.1: PIPE16 Geometry. Internal pressure (PINT) and external pressure (POUT) are input as positive values. The internal and external pressure loads are designed for closed-loop static pressure environments and therefore include pressure loads on fictitious "end caps" so that the pressure loads induce an axial stress and/or reaction in the pipe system. If a dynamic situation needs to be represented, such as a pipe venting to a lower pressure area or the internal flow is past a constriction in the pipe, these end cap loads may need to be modified by applying a nodal force normal to the cross-section of the pipe with the magnitude representing the change in pressure. Alternatively, the precomputed end cap loads can be removed using KEYOPT(8) = 1 and the appropriate end cap loads added by the user. The transverse pressures (PX, PY, and PZ) may represent wind or drag loads (per unit length of the pipe) and are defined in the global Cartesian directions. Positive transverse pressures act in the positive coordinate directions. The normal component or the projected full pressure may be used (KEYOPT(5)). Tapered pressures are not recognized. Only constant pressures are supported for this element. See PIPE16 - Elastic Straight Pipe for more information.
Temperatures may be input as element body loads at the nodes. Temperatures may have wall gradients or diametral gradients (KEYOPT(1)). The average wall temperature at θ = 0° is calculated as 2 * TAVG - T(180) and the average wall temperature at θ = -90° is calculated as 2 * TAVG - T(90). The element temperatures are assumed to be linear along the length. The first temperature at node I (TOUT(I) or TAVG(I)) defaults to TUNIF. If all temperatures after the first are unspecified, they default to the first. If all temperatures at node I are input, and all temperatures at node J are unspecified, the node J temperatures default to the corresponding node I temperatures. For any other pattern of input temperatures, unspecified temperatures default to TUNIF.
For piping analyses, the PIPE module of PREP7 may be used to generate the input for this element. KEYOPT(4) is used to identify the element type for output labeling and for postprocessing operations.
KEYOPT(7) is used to compute an unsymmetric gyroscopic damping matrix (often used for rotordynamic analyses). The rotational frequency is input with the SPIN real constant (radians/time, positive in the positive element x direction).
A summary of the element input is given in "PIPE16 Input Summary". A general description of element input is given in Element Input.
PIPE16 Input Summary
- Nodes
I, J, K (K, the orientation node, is optional)
- Degrees of Freedom
UX, UY, UZ, ROTX, ROTY, ROTZ
- Real Constants
OD, TKWALL, SIFI, SIFJ, FLEX, DENSFL, DENSIN, TKIN, TKCORR, AREAIN, MWALL, STIFF, SPIN See Table 16.1: PIPE16 Real Constants for a description of the real constants - Material Properties
EX, ALPX (or CTEX or THSX),
PRXY (or NUXY), DENS, GXY, BETD, ALPD, DMPR
- Surface Loads
- Pressures --
1-PINT, 2-PX, 3-PY, 4-PZ, 5-POUT
- Body Loads
- Temperatures --
TOUT(I), TIN(I), TOUT(J), TIN(J) if KEYOPT (1) = 0, or TAVG(I), T90(I), T180(I), TAVG(J), T90(J), T180(J) if KEYOPT (1) = 1
- Special Features
Stress stiffening Large deflection Birth and death - KEYOPT(1)
Temperatures represent:
- 0 --
The through-wall gradient
- 1 --
The diametral gradient
- KEYOPT(2)
Stress intensification factors:
- 0 --
Stress intensity factors from SIFI and SIFJ
- 1 --
Stress intensity factors at node I from tee joint calculation
- 2 --
Stress intensity factors at node J from tee joint calculation
- 3 --
Stress intensity factors at both nodes from tee joint calculation
- KEYOPT(4)
Element identification (for output and postprocessing):
- 0 --
Straight pipe
- 1 --
Valve
- 2 --
Reducer
- 3 --
Flange
- 4 --
Expansion joint
- 5 --
Mitered bend
- 6 --
Tee branch
- KEYOPT(5)
PX, PY, and PZ transverse pressures:
- 0 --
Use only the normal component of pressure
- 1 --
Use the full pressure (normal and shear components)
- KEYOPT(6)
Member force and moment output:
- 0 --
Do not print member forces or moments
- 2 --
Print member forces and moments in the element coordinate system
- KEYOPT(7)
Gyroscopic damping matrix:
- 0 --
No gyroscopic damping matrix
- 1 --
Compute gyroscopic damping matrix. Real constant SPIN must be greater than zero. DENSFL and DENSIN must be zero.
Note: The real constant MWALL is not used to compute the gyroscopic damping matrix.
- KEYOPT(8)
End cap loads:
- 0 --
Internal and external pressures cause loads on end caps
- 1 --
Internal and external pressures do not cause loads on end caps
Table 16.1: PIPE16 Real Constants
No. | Name | Description |
---|---|---|
1 | OD | Pipe outer diameter |
2 | TKWALL | Wall thickness |
3 | SIFI | Stress intensification factor (node I) |
4 | SIFJ | Stress intensification factor (node J) |
5 | FLEX | Flexibility factor |
6 | DENSFL | Internal fluid density |
7 | DENSIN | Exterior insulation density |
8 | TKIN | Insulation thickness |
9 | TKCORR | Corrosion thickness allowance |
10 | AREAIN | Insulation surface area (replaces program-calculated value) |
11 | MWALL | Pipe wall mass (replaces program-calculated value) |
12 | STIFF | Axial pipe stiffness (replaces program-calculated value) |
13 | SPIN | Rotordynamic spin (required if KEYOPT(7) = 1) |
PIPE16 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 16.2: PIPE16 Element Output Definitions
Several items are illustrated in Figure 16.2: PIPE16 Stress Output.
The direct stress (SAXL) includes the internal pressure (closed end) effect. The direct stress does not include the axial component of the transverse thermal stress (STH). The principal stresses and the stress intensity include the shear force stress component, and are based on the stresses at the two extreme points on opposite sides of the neutral axis. These quantities are calculated at the outer surface and might not occur at the same location around the pipe circumference. Angles listed in the output are measured as shown (θ) in Figure 16.2: PIPE16 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 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 16.2: PIPE16 Element Output Definitions
Name | Definition | O | R |
---|---|---|---|
EL | Element Number | Y | Y |
NODES | Nodes - I, J | Y | Y |
MAT | Material number | Y | Y |
VOLU: | Volume | - | Y |
XC, YC, ZC | Location where results are reported | Y | 6 |
CORAL | Corrosion thickness allowance | 1 | 1 |
TEMP | TOUT(I), TIN(I), TOUT(J), TIN(J) | 2 | 2 |
TEMP | TAVG(I), T90(I), T180(I), TAVG(J), T90(J), T180(J) | 3 | 3 |
PRES | PINT, PX, PY, PZ, POUT | Y | Y |
SFACTI, SFACTJ | Stress intensification factors at nodes I and J | Y | Y |
STH | Stress due to maximum thermal gradient through the wall thickness | Y | Y |
SPR2 | Hoop pressure stress for code calculations | - | Y |
SMI, SMJ | Moment stress at nodes I and J for code calculations | - | Y |
SDIR | Direct (axial) stress | - | Y |
SBEND | Maximum bending stress at outer surface | - | Y |
ST | Shear stress at outer surface due to torsion | - | Y |
SSF | Shear stress due to shear force | - | Y |
S:(1MX, 3MN, INTMX, EQVMX) | Maximum principal stress, minimum principal stress, maximum stress intensity, maximum equivalent stress (over eight points on the outside surface at both ends of the element) | Y | Y |
S:(AXL, RAD, H, XH) | Axial, radial, hoop, and shear stresses | 4 | 4 |
S:(1, 3, INT, EQV) | Maximum principal stress, minimum principal stress, stress intensity, equivalent stress | 4 | 4 |
EPEL:(AXL, RAD, H, XH) | Axial, radial, hoop, and shear strains | 4 | 4 |
EPTH:(AXL, RAD, H) | Axial, radial, and hoop thermal strain | 4 | 4 |
MFOR:(X, Y, Z) | Member forces for nodes I and J (in the element coordinate system) | 5 | Y |
MMOM:(X, Y, Z) | Member moments for nodes I and J (in the element coordinate system) | 5 | Y |
The item repeats at 0°, 45°, 90°, 135°, 180°, 225°, 270°, 315° at node I, then at node J, all at the outer surface.
Available only at centroid as a *GET item.
The following tables list output available through the ETABLE command using the Sequence Number method. See The General Postprocessor (POST1) in the Basic Analysis Guide and The Item and Sequence Number Table of this manual for more information. The following notation is used in Table 16.3: PIPE16 Item and Sequence Numbers (Node I) through Table 16.5: PIPE16 Item and Sequence Numbers:
- Name
output quantity as defined in the Table 16.2: PIPE16 Element Output Definitions
- Item
predetermined Item label for ETABLE command
- E
sequence number for single-valued or constant element data
- I, J
sequence number for data at nodes I and J
Table 16.3: PIPE16 Item and Sequence Numbers (Node I)
Output Quantity Name | ETABLE and ESOL Command Input | |||||||||
---|---|---|---|---|---|---|---|---|---|---|
Item | E | Circumferential Location | ||||||||
0° | 45° | 90° | 135° | 180° | 225° | 270° | 315° | |||
SAXL | LS | - | 1 | 5 | 9 | 13 | 17 | 21 | 25 | 29 |
SRAD | LS | - | 2 | 6 | 10 | 14 | 18 | 22 | 26 | 30 |
SH | LS | - | 3 | 7 | 11 | 15 | 19 | 23 | 27 | 31 |
SXH | LS | - | 4 | 8 | 12 | 16 | 20 | 24 | 28 | 32 |
EPELAXL | LEPEL | - | 1 | 5 | 9 | 13 | 17 | 21 | 25 | 29 |
EPELRAD | LEPEL | - | 2 | 6 | 10 | 14 | 18 | 22 | 26 | 30 |
EPELH | LEPEL | - | 3 | 7 | 11 | 15 | 19 | 23 | 27 | 31 |
EPELXH | LEPEL | - | 4 | 8 | 12 | 16 | 20 | 24 | 28 | 32 |
EPTHAXL | LEPTH | - | 1 | 5 | 9 | 13 | 17 | 21 | 25 | 29 |
EPTHRAD | LEPTH | - | 2 | 6 | 10 | 14 | 18 | 22 | 26 | 30 |
EPTHH | LEPTH | - | 3 | 7 | 11 | 15 | 19 | 23 | 27 | 31 |
MFORX | SMISC | 1 | - | - | - | - | - | - | - | - |
MFORY | SMISC | 2 | - | - | - | - | - | - | - | - |
MFORZ | SMISC | 3 | - | - | - | - | - | - | - | - |
MMOMX | SMISC | 4 | - | - | - | - | - | - | - | - |
MMOMY | SMISC | 5 | - | - | - | - | - | - | - | - |
MMOMZ | SMISC | 6 | - | - | - | - | - | - | - | - |
SDIR | SMISC | 13 | - | - | - | - | - | - | - | - |
ST | SMISC | 14 | - | - | - | - | - | - | - | - |
S1 | NMISC | - | 1 | 6 | 11 | 16 | 21 | 26 | 31 | 36 |
S3 | NMISC | - | 3 | 8 | 13 | 18 | 23 | 28 | 33 | 38 |
SINT | NMISC | - | 4 | 9 | 14 | 19 | 24 | 29 | 34 | 39 |
SEQV | NMISC | - | 5 | 10 | 15 | 20 | 25 | 30 | 35 | 40 |
SBEND | NMISC | 90 | - | - | - | - | - | - | - | - |
SSF | NMISC | 91 | - | - | - | - | - | - | - | - |
TOUT | LBFE | - | 4 | - | 1 | - | 2 | - | 3 | - |
TIN | LBFE | - | 8 | - | 5 | - | 6 | - | 7 | - |
Table 16.4: PIPE16 Item and Sequence Numbers (Node J)
Output Quantity Name | ETABLE and ESOL Command Input | |||||||||
---|---|---|---|---|---|---|---|---|---|---|
Item | E | Circumferential Location | ||||||||
0° | 45° | 90° | 135° | 180° | 225° | 270° | 315° | |||
SAXL | LS | - | 33 | 37 | 41 | 45 | 49 | 53 | 57 | 61 |
SRAD | LS | - | 34 | 38 | 42 | 46 | 50 | 54 | 58 | 62 |
SH | LS | - | 35 | 39 | 43 | 47 | 51 | 55 | 59 | 63 |
SXH | LS | - | 36 | 40 | 44 | 48 | 52 | 56 | 60 | 64 |
EPELAXL | LEPEL | - | 33 | 37 | 41 | 45 | 49 | 53 | 57 | 61 |
EPELRAD | LEPEL | - | 34 | 38 | 42 | 46 | 50 | 54 | 58 | 62 |
EPELH | LEPEL | - | 35 | 39 | 43 | 47 | 51 | 55 | 59 | 63 |
EPELXH | LEPEL | - | 36 | 40 | 44 | 48 | 52 | 56 | 60 | 64 |
EPTHAXL | LEPTH | - | 33 | 37 | 41 | 45 | 49 | 53 | 57 | 61 |
EPTHRAD | LEPTH | - | 34 | 38 | 42 | 46 | 50 | 54 | 58 | 62 |
EPTHH | LEPTH | - | 35 | 39 | 43 | 47 | 51 | 55 | 59 | 63 |
MFORX | SMISC | 7 | - | - | - | - | - | - | - | - |
MFORY | SMISC | 8 | - | - | - | - | - | - | - | - |
MFORZ | SMISC | 9 | - | - | - | - | - | - | - | - |
MMOMX | SMISC | 10 | - | - | - | - | - | - | - | - |
MMOMY | SMISC | 11 | - | - | - | - | - | - | - | - |
MMOMZ | SMISC | 12 | - | - | - | - | - | - | - | - |
SDIR | SMISC | 15 | - | - | - | - | - | - | - | - |
ST | SMISC | 16 | - | - | - | - | - | - | - | - |
S1 | NMISC | - | 41 | 46 | 51 | 56 | 61 | 66 | 71 | 76 |
S3 | NMISC | - | 43 | 48 | 53 | 58 | 63 | 68 | 73 | 78 |
SINT | NMISC | - | 44 | 49 | 54 | 59 | 64 | 69 | 74 | 79 |
SEQV | NMISC | - | 45 | 50 | 55 | 60 | 65 | 70 | 75 | 80 |
SBEND | NMISC | 92 | - | - | - | - | - | - | - | - |
SSF | NMISC | 93 | - | - | - | - | - | - | - | - |
TOUT | LBFE | - | 12 | - | 9 | - | 10 | - | 11 | - |
TIN | LBFE | - | 16 | - | 13 | - | 14 | - | 15 | - |
PIPE16 Assumptions and Restrictions
The pipe must not have a zero length or wall thickness. In addition, the OD must not be less than or equal to zero, the ID must not be less than zero, and the corrosion thickness allowance must be less than the wall thickness.
The element temperatures are assumed to vary linearly along the length.
The element may be used for both thin and thick-walled situations; however, some of the stress calculations are based on thin-wall theory.
The pipe element is assumed to have "closed ends" so that the axial pressure effect is included.
Shear deflection capability is also included in the element formulation.
Eigenvalues calculated in a gyroscopic modal analysis can be very sensitive to changes in the initial shift value, leading to potential error in either the real or imaginary (or both) parts of the eigenvalues.
PIPE16 Product Restrictions
When used in the product(s) listed below, the stated product-specific restrictions apply to this element in addition to the general assumptions and restrictions given in the previous section.
Ansys Professional —
The SPIN real constant (R13) is not available.
The only special features allowed are stress stiffening and large deflections.
KEYOPT(7) (gyroscopic damping) is not allowed.