6.5.3 Pressure loadings on elbow elements

Elbow elements are often used to model pipelines in which the curvature of the pipe can change significantly while the pipe is subjected to uniform or hydrostatic pressure. Therefore, pressure loadings that include large geometry changes are developed for these elements, as described in this section.
The virtual work contribution of pressure on the lateral surface of the elbow is
where

is the pressure magnitude,

is the variational displacement at a point on the midsurface of the lateral wall of the elbow,

is the normal to the lateral wall midsurface, and

is the area of space occupied by the lateral surface in the current configuration.
The product of the surface normal and the differential area can be rewritten in terms of material coordinates
along the pipe and
around the pipe section:
where

is the initial pipe radius, so that
where

is the area of space occupied by the lateral surface in the reference configuration.
For hydrostatic pressure, the pressure magnitude is a function of position:
where

is the reference pressure magnitude,

is the zero pressure height,

is the reference height, and

is
, a unit vector in the vertical direction.
In the elbow elements position on the lateral surface,
, is interpolated as
where

is the position of a point on the pipe axis,

is the radial displacement,

is the tangential displacement,

is
, and
The first variation of the position can now be expressed as
and the derivatives of the position with respect to the parametrization are
and
Assuming that (1) terms in

and

can be ignored due to negligible twist in the pipe, (2) terms in

and its derivatives can be ignored due to negligible warping in the pipe, (3)

, (4) the stretch

is unity, and (5)

and

are small compared to

, we arrive at the following expression for the integrand of

:
For closed-end loading the virtual work contribution of pressure on the end-caps of the elbow is
where

and

are the material coordinates in a two-dimensional cylindrical coordinate system of points on the end-caps of the elbow element,

represents position on the end-caps, and

is the area of space occupied by the end-caps of the elbow element. We assume the following deformation for the end-caps:
where

is the parameter that identifies the end-cap being considered. The assumed deformation arises naturally on considering a deformation of the end-caps in which radial rays of the reference end-cap configuration remain straight lines under deformation. It can be shown easily that the assumed deformation of the end-caps is differentiable as long as the deformations of the circumferential curves of the end-caps are differentiable. For the end-cap boundary shapes that arise in applications (primarily ovalized modes), the assumed deformation will be locally invertible so that integration of functions over the deformed surface is not likely to be a problem.
Ignoring the terms due to warping in the expression for position on the lateral surface, the first variation of position on an end-cap is
and the derivatives of

with respect to

and

are
and
The integrand of the expression for the virtual work of pressure on the end-caps can now be expressed as
where

is

if the center of the end-cap is node 1 of the element and

if the center is node 2 or 3 of the element.
The load stiffness for the pressure loading, which by definition is the first variation of the virtual work of the pressure load, is given by
. The following expressions are required for its calculation:
and
In the above

is nonzero only in the case of hydrostatic pressure, when it is given in the first case by
and in the second case by