3.5.3 Euler-Bernoulli beam elements
In these elements it is assumed that the internal virtual work rate is associated with axial strain and torsional shear only. Further, it is assumed that the cross-section does not deform in its plane or warp out of its plane, and that this cross-sectional plane remains normal to the beam axis. These are the classical assumptions of the Euler-Bernoulli beam theory, which provides satisfactory results for slender beams.
Let be material coordinates such that S locates points on the beam axis and measures distance in the cross-section. In addition, let be unit vectors normal to the beam axis in the current configuration: . Then the position of a point of the beam in the current configuration is
where
is the point on the beam axis of the cross-section containing
. Then
and so length on the fiber at
is measured in the current configuration as
Now since the beam is slender, we will neglect terms of second-order in
g and
h, the distance measuring material coordinates in the cross-section. Thus,
Strain measures
The internal virtual work rate associated with axial stress is
where
and
are any material stress and strain measures associated with axial deformation at the point
of the beam, since strains are assumed to be small. For this purpose we will use Green's strain so that
where
, the square of the ratio of current configuration length to reference configuration length in the axial direction on the fiber. From
Equation 3.5.3–1 and its equivalent in the reference configuration, we have
Again, neglecting all but first-order terms in
g and
h because of the slenderness assumption, this becomes
where
This simplification allows us to write the internal virtual work rate associated with axial stress as
where
with
and
Now the cross-sectional base vectors
and
are assumed to remain normal to the beam axis, so
Hence,
So we have
and
This defines the generalized strains associated with axial stretch. For torsional strain the internal virtual work rate is
where
and
is the proportionality constant between shear strains and torsional strains (see
Beam element formulation, Section 3.5.2, for details). For computational simplicity the form of the torsional strains is taken to be
where
and
This assumes a linear interpolation of rotation
along the beam. Thus, the generalized strain measures for these beams are
and
the beam curvature change measures, and
With these measures, the internal virtual work rate can be written
Internal virtual work rate Jacobian
For the Jacobian matrix of the overall Newton method, Equation 2.1.1–3, the variation of this internal work rate with respect to nodal displacement variations must be formed. The constitutive theory is written as
and so
where
is the section stiffness matrix obtained by integration over the cross-section. See
Beam element formulation, Section 3.5.2, for more information on section integration.
First variations of strains
Taking the variations of the above strain definitions gives directly
In these expressions we need
and
as well as
; these are now derived. From the expression for
, namely
another ancillary vector
, normal to the tangent, is defined by
so that
In addition,
so that
Since
form an orthonormal triad,
because
. From the definition of
, it is straightforward to show that
So
and
Thus,
We can also write
and
Combining terms appropriately,
Hence, from
Equation 3.5.3–2 In a similar manner one can show that
The first variations of strain become
so that
In addition,
so that
since
form an orthonormal triad,
because
Second variations of strains
The second variation of the axial strain is simply
To compute the second variations of bending strain, we need expressions for
. These are obtained by approximating
from which
Using these expressions, the second variations of bending strains are written as
and
For the torsional strain contribution to the initial stress matrix, we approximate
This matrix is again unsymmetric.
Interpolation
In the virtual work equation the strains include second derivatives of displacement. For this reason continuity of rotation as well as displacement is needed so that the Hermitian polynomial interpolation functions are the minimum interpolation order needed. These are used here. The Hermite cubic is written in terms of the function value and its derivative at the ends of the interval:
with node 1 at
and node 2 at
.
These functions are used in ABAQUS to interpolate the components of displacement and the initial position vector , so that the elements are basically isoparametric. In addition, rotation of about the beam axis, , is interpolated linearly. This interpolation is unsatisfactory for the user, because the nodal variables are
The last four of these variables are difficult to work with; furthermore, making them the same in all elements sharing the same node causes excessive constraint of axial stretch in these elements, especially if the beam axis is not continuous through the node, as in a frame structure or “T” junction.
To avoid this difficulty, the following procedure is adopted. At a node the tangent to the beam axis is
so
Now suppose we store as degrees of freedom at the node,
where
is the rotation definition introduced above. Since the initial geometry and hence
, the initial direction of the beam axis, is known,
where
is the rotation matrix defined by
, and hence
is defined by these variables and the initial geometry. Furthermore,
is directly available from
and
. Thus, the above set is a satisfactory set of nodal variables. To eliminate the unwanted axial strain constraint, in ABAQUS the stretch
at the node of each such element is taken as an internal variable, local to the element (a third internal node is created for this purpose, and so it is not shared with neighboring elements.)
It should be remarked that the above transformation (Equation 3.5.3–3) is nonlinear. This leads to some complications—for example, the d'Alembert forces no longer have the simple form rather, a matrix replaces where and use the transformation (Equation 3.5.3–3) and its appropriate time derivatives. The resulting Jacobian is nonsymmetric; ABAQUS ignores the nonsymmetric terms.
Integration
The cross-section integration has already been discussed in the context of the shear beams—it is the same for these beams. Along the beam axis, the integration schemes are as described below.
Stiffness and internal forces
Three Gauss points are used. Two Gauss points are not sufficient because the torsional strain is not independent.
Mass and distributed loads
Three Gauss points are used. Rotary inertia is not included in the mass, except for rotation of the section about the beam axis, where it is included to avoid singularity in perfectly straight beams.