PIPE16


Elastic Straight Pipe

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

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.

Figure 16.1: PIPE16 Geometry

PIPE16 Geometry

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.NameDescription
1ODPipe outer diameter
2TKWALLWall thickness
3SIFIStress intensification factor (node I)
4SIFJStress intensification factor (node J)
5FLEXFlexibility factor
6DENSFLInternal fluid density
7DENSINExterior insulation density
8TKINInsulation thickness
9TKCORRCorrosion thickness allowance
10AREAINInsulation surface area (replaces program-calculated value)
11MWALLPipe wall mass (replaces program-calculated value)
12STIFFAxial pipe stiffness (replaces program-calculated value)
13SPINRotordynamic spin (required if KEYOPT(7) = 1)

PIPE16 Output Data

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

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.

Figure 16.2: PIPE16 Stress Output

PIPE16 Stress Output


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

NameDefinitionOR
ELElement NumberYY
NODESNodes - I, JYY
MATMaterial numberYY
VOLU:Volume-Y
XC, YC, ZCLocation where results are reportedY6
CORALCorrosion thickness allowance11
TEMPTOUT(I), TIN(I), TOUT(J), TIN(J)22
TEMPTAVG(I), T90(I), T180(I), TAVG(J), T90(J), T180(J)33
PRESPINT, PX, PY, PZ, POUTYY
SFACTI, SFACTJStress intensification factors at nodes I and JYY
STHStress due to maximum thermal gradient through the wall thicknessYY
SPR2Hoop pressure stress for code calculations-Y
SMI, SMJMoment stress at nodes I and J for code calculations-Y
SDIRDirect (axial) stress-Y
SBENDMaximum bending stress at outer surface-Y
STShear stress at outer surface due to torsion-Y
SSFShear 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)YY
S:(AXL, RAD, H, XH)Axial, radial, hoop, and shear stresses44
S:(1, 3, INT, EQV)Maximum principal stress, minimum principal stress, stress intensity, equivalent stress44
EPEL:(AXL, RAD, H, XH)Axial, radial, hoop, and shear strains44
EPTH:(AXL, RAD, H)Axial, radial, and hoop thermal strain44
MFOR:(X, Y, Z)Member forces for nodes I and J (in the element coordinate system)5Y
MMOM:(X, Y, Z)Member moments for nodes I and J (in the element coordinate system)5Y

  1. If the value is greater than 0.

  2. If KEYOPT(1) = 0

  3. If KEYOPT(1) = 1

  4. The item repeats at 0°, 45°, 90°, 135°, 180°, 225°, 270°, 315° at node I, then at node J, all at the outer surface.

  5. If KEYOPT(6) = 2

  6. 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 NameETABLE and ESOL Command Input
ItemECircumferential Location
45°90°135°180°225°270°315°
SAXLLS-1591317212529
SRADLS-26101418222630
SHLS-37111519232731
SXHLS-48121620242832
EPELAXLLEPEL-1591317212529
EPELRADLEPEL-26101418222630
EPELHLEPEL-37111519232731
EPELXHLEPEL-48121620242832
EPTHAXLLEPTH-1591317212529
EPTHRADLEPTH-26101418222630
EPTHHLEPTH-37111519232731
MFORXSMISC1--------
MFORYSMISC2--------
MFORZSMISC3--------
MMOMXSMISC4--------
MMOMYSMISC5--------
MMOMZSMISC6--------
SDIRSMISC13--------
STSMISC14--------
S1NMISC-16111621263136
S3NMISC-38131823283338
SINTNMISC-49141924293439
SEQVNMISC-510152025303540
SBENDNMISC90--------
SSFNMISC91--------
TOUTLBFE-4-1-2-3-
TINLBFE-8-5-6-7-

Table 16.4: PIPE16 Item and Sequence Numbers (Node J)

Output Quantity NameETABLE and ESOL Command Input
ItemECircumferential Location
45°90°135°180°225°270°315°
SAXLLS-3337414549535761
SRADLS-3438424650545862
SHLS-3539434751555963
SXHLS-3640444852566064
EPELAXLLEPEL-3337414549535761
EPELRADLEPEL-3438424650545862
EPELHLEPEL-3539434751555963
EPELXHLEPEL-3640444852566064
EPTHAXLLEPTH-3337414549535761
EPTHRADLEPTH-3438424650545862
EPTHHLEPTH-3539434751555963
MFORXSMISC7--------
MFORYSMISC8--------
MFORZSMISC9--------
MMOMXSMISC10--------
MMOMYSMISC11--------
MMOMZSMISC12--------
SDIRSMISC15--------
STSMISC16--------
S1NMISC-4146515661667176
S3NMISC-4348535863687378
SINTNMISC-4449545964697479
SEQVNMISC-4550556065707580
SBENDNMISC92--------
SSFNMISC93--------
TOUTLBFE-12-9-10-11-
TINLBFE-16-13-14-15-

Table 16.5: PIPE16 Item and Sequence Numbers

Output Quantity NameETABLE and ESOL Command Input
ItemE
STHSMISC17
PINTSMISC18
PXSMISC19
PYSMISC20
PZSMISC21
POUTSMISC22
SFACTINMISC81
SFACTJNMISC82
SPR2NMISC83
SMINMISC84
SMJNMISC85
S1MXNMISC86
S3MNNMISC87
SINTMXNMISC88
SEQVMXNMISC89

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.