3.8.1 Hydrostatic fluid elements

Products: ABAQUS/Standard  ABAQUS/Explicit  

ABAQUS includes a family of elements that can be used to represent fluid-filled cavities under hydrostatic conditions. These elements provide the coupling between the deformation of the fluid-filled structure and the pressure exerted by the contained fluid on the boundary of the cavity. In ABAQUS/Explicit the fluid must be compressible and the pressure is calculated from the cavity volume. In ABAQUS/Standard the fluid inside the cavity can be compressible or incompressible, with the fluid volume given as a function of the fluid pressure, p; the fluid temperature, ; and the fluid mass, m, in the cavity:

We refer to the incompressible case as a “hydraulic” fluid and to the compressible case as a “pneumatic” fluid. The volume, , derived from the fluid pressure and temperature should equal the actual volume, V, of the cavity. In ABAQUS/Standard this is achieved by augmenting the virtual work expression for the structure with the constraint equation

and the virtual work contribution due to the cavity pressure:

where is the augmented virtual work expression and is the virtual work expression for the structure without the cavity. The negative signs imply that an increase in the cavity volume releases energy from the fluid. This represents a mixed formulation in which the structural displacements and fluid pressure are primary variables. The rate of the augmented virtual work expression is obtained as

Here, represents the pressure load stiffness, and is the volume-pressure compliance of the fluid.

Since the pressure is the same for all elements in the cavity, the augmented virtual work expression can be written as the sum of the expressions for the individual elements:

Moreover, since the temperature is the same for all elements in the cavity, the fluid volume can be calculated for each element individually:

where is the element mass. Note that in the solution, the actual volume of the element may be different from the element volume:

The total fluid volume will match the volume of the cavity, however.

Hydraulic fluid with thermal expansion

In ABAQUS/Standard the fluid is incompressible by default and the fluid volume, , is dependent upon temperature but independent of the fluid pressure:

If compressibility is introduced, the fluid volume depends upon both the temperature and pressure:

where K is the fluid bulk modulus and is the reference fluid density at zero pressure and the initial temperature.

The total fluid mass in the cavity is the sum of the fluid masses of the elements making up the cavity:

The mass of a fluid element in the cavity, , is calculated from the initial fluid density, , and the initial element volume, :

where is the initial fluid pressure and is the initial temperature. The initial fluid density follows from the user-defined reference density, :

The fluid density at the current pressure and temperature, , is obtained as

where is the reference temperature for the coefficient of thermal expansion and is the mean (secant) coefficient of thermal expansion, and it is assumed that .

Thus, the fluid volume at the current pressure and temperature is

This volume can be calculated on an element by element basis:

Fluid can be added to or removed from the cavity. The amount of fluid added is given as the change in (fluid) mass . Consequently, the change of the fluid volume at the current cavity temperature is

Ideal gas

In this case the fluid is compressible, and the volume is a function of the pressure and the temperature in the cavity:

where, as before, the total fluid mass in the cavity is the sum of the masses of the elements in the cavity. The fluid is assumed to behave like an ideal gas; hence, the density of the fluid in the cavity can be calculated as

where and are the temperature and pressure at the reference density , is the temperature at absolute zero, and is the ambient pressure. The current fluid volume can again be calculated on an element by element basis:

The corresponding volume-pressure compliance is

Fluid can be added to or removed from the cavity. The amount of fluid added is again given as the change in (fluid) mass. Consequently, the change of the fluid volume at the current cavity temperature is

Volume calculation

The hydrostatic fluid elements appear as surface elements that cover the cavity boundary, but they are actually volume elements when the cavity reference node is accounted for. Figure 3.8.1–1 depicts the 4-node hydrostatic fluid volume element, F3D4. The dashed lines indicate that the element is actually pyramidal in shape.

Figure 3.8.1–1 F3D4 hydrostatic fluid element.

The volume, , of each element must be calculated. The coordinates of any point on the base of the pyramid element can be found by

where are the interpolation functions for the base of the pyramid (Solid isoparametric quadrilaterals and hexahedra, Section 3.2.4), expressed in terms of parametric coordinates g and h; are the nodal coordinates; and the summation extends over all nodes on the base. For three-dimensional elements the Jacobian on the surface is calculated as

The normal to the element face, , multiplied by an infinitesimal area, , of the element face is, hence, obtained as

The infinitesimal volume, , associated with this infinitesimal area is

where is the position of the cavity reference node. The volume of the element, , is then obtained by integration. For a quadrilateral base this yields

The integration boundaries will be different for a triangular base. Introducing the relative position, , this becomes

The variation in the element volume is readily obtained as

This expression includes contributions to the volume change due to variation in the “sides” of the pyramid element. Hence, the equivalent forces will also include the effect of the pressure on these sides. The pressure on the sides will be balanced by the pressure on the side of the adjacent pyramid element; hence, we only need to calculate the contribution from the pressure on the “base” of the pyramid. This can be done by separating the contributions by using partial integration, which yields

The last two integrals represent the contributions on the sides of the pyramid; hence, the resulting expression is

This equation can readily be obtained directly (see Pressure load stiffness, Section 6.5.1).

The second variation of the expression for the volume is

Similar expressions can be obtained for planar and axisymmetric fluid elements. For a planar element the volume is readily obtained as

where is the unit vector in the direction perpendicular to the plane. Similarly, for axisymmetric elements

The first and second variations are obtained in the same way as for three-dimensional elements.

The integrations can be carried out analytically. For instance, for element type F3D4 the above expressions yield, after some manipulation,

where , denote the relative nodal coordinates. The first variation (involving the base only) is equal to

and the second variation of the volume is

Fluid link element

In addition to the fluid cavity elements ABAQUS also offers a 2-node fluid link element that can be used to model fluid flow between two cavities or between a cavity and the outside world. This is typically used when the fluid has to flow through a narrow orifice. It is assumed that the mass flow rate, q, through the link is a homogeneous function of the pressure difference, . In addition, it is assumed that the flow rate may depend on the average temperature——and, for a compressible fluid, on the average pressure—:

The flow rate needs to be integrated over a finite increment. We assume that the dependence on the average pressure, , is weak. Hence, we use a semi-implicit method: we use at the end of the increment and at the start of the increment. For the temperature we choose to be the average of at the start and end of the increment, because this is likely to be the most accurate. Hence, we obtain the mass flow through the link:

This mass change needs to be converted to a volume change in each cavity. It is assumed that the fluid in both cavities is the same, but the pressures (and possibly the temperatures) may be different. With use of the pressure and temperature-dependent density for each cavity i, the relations become

and

Observe that, except for isothermal incompressible fluids, the resulting matrix is nonsymmetric. For an incompressible fluid with thermal expansion, the nonsymmetric contribution is small and the equations can be symmetrized without significant degradation in convergence rate. For an ideal gas, the nonsymmetric terms can have a significant effect if the drop in pressure is large.

Many applications of fluid flow through a fluid link involve dynamic loading in the form of steady-state vibration; and often in such cases the dissipative losses in the fluid link must be modeled to obtain useful results. In most problems of this class the fluids on either side of the fluid link are first pressurized statically. In the implementation in ABAQUS/Standard, it is assumed that the vibration amplitude is sufficiently small that the fluid link response in the dynamic phase of the problem can be treated as linear perturbations about the prepressurized state.

For small vibrations about a prepressurized state we linearize Equation 3.8.1–1 to give

The mass flow through the link can, therefore, be derived as follows:

Substituting the above expression for mass flow into Equation 3.8.1–2 and noting that yields

In ABAQUS this model is provided only for the direct-solution steady-state dynamic analysis procedure.

Negative eigenvalues

It is possible that negative eigenvalues will be encountered in the solution of certain hydrostatic fluid element problems. With standard elements this can indicate that a bifurcation or buckling load has been exceeded. However, this is not necessarily true with hydrostatic fluid elements; negative eigenvalues can result solely from the numerical implementation.

Consider the simple hydrostatic fluid model depicted in Figure 3.8.1–2.

Figure 3.8.1–2 Simple hydrostatic fluid model.

If the fluid is considered incompressible, application of the downward force causes the fluid to compress vertically and expand horizontally, while maintaining the original fluid volume. Thus, the model can be adequately discretized as a three degree of freedom system: the horizontal displacement of the right platen, u; the vertical displacement of the top platen, v; and the fluid pressure, p. The corresponding system of equations in matrix form is

The equation solver processes the equations in sequence. Hence, it will process the submatrix relating the displacements to the forces prior to processing the constraints. This leads to a negative eigenvalue if . However since the mode associated with the negative eigenvalue is subsequently constrained by the continuity equation, no instability occurs.

Reference