FLUID81


Axisymmetric-Harmonic Contained Fluid

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

FLUID81 Element Description

Although this archived element is available for use in your analysis, Ansys, Inc. recommends using a current-technology element such as FLUID29, FLUID30, FLUID220, or FLUID221.

FLUID81 is a modification of the axisymmetric structural solid element (PLANE25). The element is used to model fluids contained within vessels having no net flow rate. It is defined by four nodes having three degrees of freedom at each node: translations in the nodal x, y, and z directions. The element is used in a structural analysis as an axisymmetric ring element.

The element is a generalization of the axisymmetric version of FLUID79, the 2D fluid element, in that the loading need not be axisymmetric. Various loading cases are described in Harmonic Axisymmetric Elements with Nonaxisymmetric Loads. The fluid element is particularly well suited for calculating hydrostatic pressures and fluid/solid interactions. Acceleration effects, such as in sloshing problems, as well as temperature effects, may be included. See FLUID81 - Axisymmetric-Harmonic Contained Fluid for more details about this element. Another fluid element (FLUID116) is available to model fluids flowing in pipes and channels.


Note:  This element cannot be used in a modal analysis.


Figure 81.1: FLUID81 Geometry

FLUID81 Geometry

FLUID81 Input Data

The geometry, node locations, and the coordinate system for this element are shown in Figure 81.1: FLUID81 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 isotropic material properties. If MODE = 0 and ISYM = 1, the element behaves similar to the axisymmetric case of FLUID79. The MODE and ISYM parameters are discussed in detail in Harmonic Axisymmetric Elements with Nonaxisymmetric Loads. EX, which is interpreted as the "fluid elastic modulus," should be the bulk modulus of the fluid (approximately 300,000 psi for water). The viscosity property (VISC) is used to compute a damping matrix for dynamic analyses. A typical viscosity value for water is 1.639 x 10-7 lb-sec/in2. Density (DENS) must be input as a positive number.

The use of KEYOPT(2) for gravity springs is discussed in "FLUID80 Input Data". Vertical acceleration (ACELY on the ACEL command) is needed for the gravity springs regardless of the value of MODE. 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 81.1: FLUID81 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.

A summary of the element input is given in "FLUID81 Input Summary". A general description of element input is given in Element Input.

FLUID81 Input Summary

Nodes

I, J, K, L

Degrees of Freedom

UX, UY, UZ

Real Constants

None

Material Properties

MP command: EX, ALPX (or CTEX or THSX), DENS, VISC, ALPD, BETD, DMPR

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

None

KEYOPT(2)

Location of gravity springs:

0 -- 

Place gravity springs on all sides of all elements

1 -- 

Place gravity springs only on face of elements located on Y = 0.0 plane (element must not have positive Y coordinates)

FLUID81 Output Data

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

The pressure and temperature are evaluated at the element centroid. Nodal forces and reaction forces are on a full 360° basis.

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°. Printout for combined loading cases may be obtained from the POST1 routine. 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.

A general description of solution output is given in Solution Output. See the Basic Analysis Guide for ways to view results.

The following notation is used in Table 81.1: FLUID81 Element Output Definitions:

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 81.1: FLUID81 Element Output Definitions

NameDefinitionOR
ELElement NumberYY
NODESNodes - I, J, K, LYY
MATMaterial numberYY
ISYMLoading Key11
MODENumber of waves in loadingYY
VOLU:VolumeYY
XC, YCLocation where results are reportedY2
PRESPressures P1 at nodes J,I; P2 at K,J; P3 at L,K; P4 at I,LYY
TEMPTemperatures T(I), T(J), T(K), T(L)YY
TAVGAverage temperatureY-
PAVGAverage pressureYY

  1. If ISYM is:

    1 - Symmetric loading

    -1 - Antisymmetric loading

  2. Available only at centroid as a *GET item.

Table 81.2: FLUID81 Item and Sequence Numbers lists 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 for more information. The following notation is used in Table 81.2: FLUID81 Item and Sequence Numbers:

Name

output quantity as defined in the Table 81.1: FLUID81 Element Output Definitions

Item

predetermined Item label for ETABLE command

E

sequence number for single-valued or constant element data

I,J,...,L

sequence number for data at nodes I,J,...,L

Table 81.2: FLUID81 Item and Sequence Numbers

Output Quantity NameETABLE and ESOL Command Input
ItemEIJKL
PRESSMISC1----
P1SMISC-32--
P2SMISC--54-
P3SMISC---76
P4SMISC-8--9

FLUID81 Assumptions and Restrictions

  • The area of the element must be positive.

  • The fluid element must lie in an X-Y plane as shown in Figure 81.1: FLUID81 Geometry and the Y-axis must be the axis of symmetry for axisymmetric analyses.

  • An axisymmetric structure should be modeled in the +X quadrants.

  • The Y-axis should be oriented in the vertical direction and the top surface is usually at Y = 0.0.

  • The element temperature is taken to be the average of the nodal temperatures.

  • Temperature-dependent material properties, if any, are evaluated at the reference temperature (TREF).

  • Elements should be rectangular since results are known to be of lower quality for nonrectangular shapes.

  • The nonlinear transient dynamic analysis should be used instead of the linear transient dynamic analysis for this element.

  • A lumped mass matrix may be obtained for this element with the LUMPM command.

  • See FLUID80 for more assumptions and restrictions.

FLUID81 Product Restrictions

There are no product-specific restrictions for this element.