Kinematics of rods

Application of rotations to the dynamics of rods can be traced to the works of Serret [1] and Frenet [2] on space curves in the early 1850s and Kirchhoff’s work [3] on rod theory in the late 1850s. These works feature triads of vectors that rigidly rotate as one moves along a curve in space, and the resulting rate of rotation serves to specify curvatures, torsions, and strains. To explore these topics, we begin with a discussion of space curves and present the Serret-Frenet relations and the Darboux vector; these ideas are applied to a variety of curves for illustration. We then introduce the concept of a material curve and follow with a discussion of Kirchhoff’s celebrated rod theory.

Space curves

The geometric features of space curves, or curves that are embedded in three-dimensional space, are a classic topic in differential geometry. The notion of curvature of a plane curve was extended to a space curve by Serret in 1851 and by Frenet in 1852 (see [1, 2, 4]). Their work forms the basis for most treatments of space curves in modern texts on differential geometry (e.g., [5, 6, 7]). The primary resource for the content we present here is the remarkable text by Kreyszig [5]; our discussion of the forthcoming examples of space curves is adapted from the textbook [8].

The Frenet triad, curvature, and torsion

Consider a space curve \mathcal{C} in Euclidean three-dimensional space \mathbb{E}^3, as illustrated in Figure 1. We assume the curve is parameterized by an arc-length parameter s such that any point along the curve can be located relative to a fixed origin using a set of Cartesian coordinates x_i = x_i(s) (i = 1, \, 2, \, 3) with corresponding orthonormal, right-handed basis \{ \mathbf{E}_i \} by specifying the value of s at the point of interest.

Figure 1. Evolution of the Frenet triad  \{{\bf e}_t, \, {\bf e}_n, \, {\bf e}_b\} for a general space curve  \mathcal{C}.

Therefore, we can express the position vector \mathbf{r} of any point on the curve as follows:

(1)   \begin{equation*} {\bf r} = {\bf r}(s) = x_1(s){\bf E}_1 + x_2(s){\bf E}_2 + x_3(s){\bf E}_3. \end{equation*}

The unit vector {\bf e}_t that is tangent to the curve is given by

(2)   \begin{equation*} {\bf e}_t = \frac{\partial {\bf r}}{\partial s}, \end{equation*}

and the derivative of this vector defines the curvature \kappa and the unit normal vector {\bf e}_n:

(3)   \begin{equation*} \frac{\partial {\bf e}_t}{\partial s} = \kappa {\bf e}_n. \end{equation*}

That is,

(4)   \begin{eqnarray*} && \kappa = \left|\!\left| \frac{\partial {\bf e}_t}{\partial s}\right|\!\right|, \\ \\ \\ && {\bf e}_n = \frac{1}{\kappa}\frac{\partial {\bf e}_t}{\partial s}. \end{eqnarray*}

We then use the unit tangent and normal vectors to construct an orthonormal and right-handed triad known as the Frenet triad:

(5)   \begin{equation*} \left\{ {\bf e}_t, \, {\bf e}_n, \, {\bf e}_b = {\bf e}_t\times{\bf e}_n \right\}, \end{equation*}

where {\bf e}_b is referred to as the unit binormal vector. For a curve whose point locations also vary with time t, i.e., when \mathbf{r} = \mathbf{r}(s, \, t), we can still calculate the Frenet triad, but we must do so at each instant in time. Using the fact that the Frenet triad is orthonormal, we can define the geometric torsion of the space curve, \tau, by the relation

(6)   \begin{equation*} \frac{\partial {\bf e}_b}{\partial s} = -\tau {\bf e}_n. \end{equation*}

In this definition, the minus sign is a convention. A space curve is said to be right-handed if \tau > 0 and left-handed if \tau < 0 (see [5]). Notice that it is possible to define the curvature \kappa and geometric torsion \tau without referring explicitly to the Frenet triad:

(7)   \begin{eqnarray*} && \kappa = \left|\!\left| \frac{\partial^2 {\bf r}}{\partial s^2}\right|\!\right|, \\ \\ \\ && \tau = \left[ \frac{\partial {\bf r}}{\partial s}, \, \, \frac{\partial^2 {\bf r}}{\partial s^2}, \, \, \frac{\partial^3 {\bf r}}{\partial s^3} \right], \end{eqnarray*}

where [{\bf a}, \, {\bf b}, \, {\bf c}] denotes the scalar triple product. Finally, we mention that there are several methods of defining the twist of a space curve. The most popular approach is to use the geometric torsion \tau obtained from the Frenet triad to define the total torsion of a space curve, T_\textrm{w}:

(8)   \begin{equation*} T_\textrm{w} = \frac{1}{2\pi} \int^\ell_0 \tau \textrm{d} s \end{equation*}

for a curve segment of length \ell. Note the conventional division by 2 \pi in this definition.

The Serret-Frenet relations

The Serret-Frenet relations are compact expressions of the rate of change of the Frenet triad basis vectors expressed in the basis \left\{ {\bf e}_t, \, {\bf e}_n, \, {\bf e}_b \right\}. These relations are obtained by using the definitions (3) and (6) and by differentiating {\bf e}_n = {\bf e}_b\times{\bf e}_t:

(9)   \begin{eqnarray*} && \frac{\partial {\bf e}_t}{\partial s} = \kappa {\bf e}_n, \\ \\ \\ && \frac{\partial {\bf e}_n}{\partial s} = - \kappa {\bf e}_t + \tau {\bf e}_b, \\ \\ \\ && \frac{\partial {\bf e}_b}{\partial s} = - \tau {\bf e}_n. \end{eqnarray*}

We can write these relations using the compact notation

(10)   \begin{equation*} \frac{\partial {\bf f}_i}{\partial s} = \bomega_{\textrm{SF}}\times{\bf f}_i \quad \left (i = 1, \, 2 , \, 3 \right ), \end{equation*}

where \left\{ {\bf f}_1, \, {\bf f}_2, \, {\bf f}_3 \right\} = \left\{ {\bf e}_t, \, {\bf e}_n, \, {\bf e}_b \right\} and

(11)   \begin{equation*} \bomega_{\textrm{SF}} = \tau {\bf e}_t + \kappa {\bf e}_b. \end{equation*}

Noting that the Frenet triad is a right-handed basis, the compact form (10) is a statement of the fact that, locally, the triad can be viewed as rigidly rotating along the curve. That is,

(12)   \begin{equation*} {\bf f}_i(s) = {\bf Q}_{\textrm{SF}}(s) {\bf f}_i(0) , \end{equation*}

where {\bf Q}_{\textrm{SF}} is a rotation tensor whose axial vector

(13)   \begin{equation*} \bomega_{\textrm{SF}} = \mbox{ax}\left(\frac{\partial {\bf Q}_{\textrm{SF}}}{\partial s} {\bf Q}^T_{\textrm{SF}}\right). \end{equation*}

The vector \bomega_{\textrm{SF}}, which has the equivalent representations (11) and (13), is often called the Darboux vector. Provided \bomega_{\textrm{SF}}(s) and the initial conditions {\bf e}_t\left(s_0\right), {\bf e}_n\left(s_0\right), and {\bf e}_b\left(s_0\right) for the Frenet triad, we can integrate (9) to determine the evolution of the triad as it moves along the curve. Using the result for {\bf e}_t(s) and integrating (2) with initial condition \mathbf{r}\left(s_0\right) then yields the profile of the space curve, \mathbf{r}(s). Lastly, we note that at points on the curve where the curvature \kappa vanishes, the unit normal vector {\bf e}_n and the unit binormal vector {\bf e}_b typically rotate through 180^\circ about the unit tangent vector {\bf e}_t. At such points, the rotation tensor {\bf Q}_{\textrm{SF}} is not differentiable as a function of s and \bomega_{\textrm{SF}} as given by (11) cannot be obtained from (13).

A plane curve

To illustrate the previous developments related to the Frenet triad, we take as an example a curve in the vertical plane, like the one shown in Figure 2.

Figure 2. Evolution of the Frenet triad  \{{\bf e}_t, \, {\bf e}_n, \, {\bf e}_b\} for a plane curve  x_2 = f(x_1).

For this plane curve, we can take the position vector of a point on the curve, parameterized by the arc-length parameter s, as

(14)   \begin{equation*} {\bf r} = {\bf r}(s) = x_1(s){\bf E}_1 + x_2(s){\bf E}_2. \end{equation*}

From (2), the corresponding unit tangent vector is then

(15)   \begin{equation*} {\bf e}_t = \frac{\partial x_1}{\partial s}{\bf E}_1 + \frac{\partial x_2}{\partial s}{\bf E}_2. \end{equation*}

Suppose we define an angle \beta such that

(16)   \begin{eqnarray*} && \cos(\beta) = \frac{\partial x_1}{\partial s}, \\ \\ \\ && \sin(\beta) = \frac{\partial x_2}{\partial s}, \end{eqnarray*}

which allow us to establish two orthogonal unit vectors:

(17)   \begin{eqnarray*} && {\bf e}_1 = \cos(\beta){\bf E}_1 + \sin(\beta){\bf E}_2, \\ \\ && {\bf e}_2 = \cos(\beta){\bf E}_2 - \sin(\beta){\bf E}_1. \end{eqnarray*}

We can use these vectors and \beta to express {\bf e}_t and its derivative in convenient forms:

(18)   \begin{eqnarray*} && {\bf e}_t = {\bf e}_1, \\ \\[0.10in] && \frac{\partial {\bf e}_t}{\partial s} = \frac{\partial \beta}{\partial s} {\bf e}_2. \end{eqnarray*}

Therefore, from (4), the plane curve’s curvature and the unit normal vector have the representations

(19)   \begin{eqnarray*} && \kappa = \left| \frac{\partial \beta}{\partial s} \right|, \\ \\ \\ && {\bf e}_n = \mbox{sgn}\left( \frac{\partial \beta}{\partial s} \right){\bf e}_2, \end{eqnarray*}

where changes in the sign of \mbox{sgn}\left( \frac{\partial \beta}{\partial s} \right) can be interpreted as rotations about {\bf e}_t through 180^\circ. Taking the cross product of {\bf e}_t with {\bf e}_n yields the unit binormal vector

(20)   \begin{equation*} {\bf e}_b = \mbox{sgn}\left( \frac{\partial \beta}{\partial s} \right) {\bf E}_3. \end{equation*}

Thus, from (6), the geometric torsion \tau for a plane curve is zero. Consequently, (8) reveals that the total torsion T_\textrm{w} is also zero. With some further calculations, we find that the Darboux vector for the plane curve is

(21)   \begin{equation*} \bomega_{\textrm{SF}} = \frac{\partial \beta}{\partial s} {\bf E}_3. \end{equation*}

It is straightforward to generate a representation for the corresponding rotation tensor {\bf Q}_{\textrm{SF}}(s) where neither \frac{\partial \beta(s)}{\partial s} nor its initial value \frac{\partial \beta(0)}{\partial s} are zero:

(22)   \begin{equation*} {\bf Q}_{\textrm{SF}}(s) = {\bf e}_1(s)\otimes{\bf e}_1(0) + \mbox{sgn}\left( \frac{\partial \beta(s)}{\partial s} \right)\mbox{sgn}\left( \frac{\partial \beta(0)}{\partial s} \right)\left({\bf e}_2(s)\otimes{\bf e}_2(0) + {\bf E}_3\otimes{\bf E}_3\right). \hspace{1in} \scalebox{0.001}{\textrm{\textcolor{white}{.}}} \end{equation*}

It is illuminating to consider a few interesting cases of plane curves. For example, for the sinusoidal plane curve x_2 = A\sin\left( \frac{\pi x_1}{\ell}\right), the curvature \kappa = 0 at the points x_1 = \pm n \ell (n = 0, \, 1, \, 2, \, \dots), and thus {\bf Q}_{\textrm{SF}} is not defined at these points. Also, a plane curve with constant non-zero curvature \kappa and geometric torsion \tau = 0 is the arc of a circle. Because this curve is indicative of a constant bending moment, circular arcs feature prominently in design.

A circular helix

We now consider another example of a space curve: a circular helix, such as the ones depicted in Figure 3. A circular helix with radius R can be defined by the position vector

(23)   \begin{equation*} {\bf r} = R{\bf e}_r + R \alpha \theta{\bf E}_3, \end{equation*}

where \theta is a cylindrical polar coordinate such that the unit vectors

(24)   \begin{eqnarray*} && {\bf e}_r = \cos(\theta){\bf E}_1 + \sin(\theta){\bf E}_2, \\ \\ && {\bf e}_\theta = \cos(\theta){\bf E}_2 - \sin(\theta){\bf E}_1. \end{eqnarray*}

It is common to describe a helix in terms of the pitch parameter \alpha, which is related to the helical curve’s pitch angle \gamma:

(25)   \begin{equation*} \alpha = \tan \left( \gamma \right). \end{equation*}

If \alpha > 0, or equivalently \gamma > 0, then the helix is said to be right-handed, whereas the helix is left-handed if \alpha < 0, or equivalently \gamma < 0.

Figure 3. Evolution of the Frenet triad  \{{\bf e}_t, \, {\bf e}_n, \, {\bf e}_b\} for (a) a left-handed circular helix and (b) a right-handed circular helix.

We can determine the Frenet triad for the helix by first differentiating the position vector {\bf r} with respect to the arc-length parameter s and using the chain rule to compute the unit tangent vector {\bf e}_t:

(26)   \begin{equation*} \frac{\partial {\bf r}}{\partial s} = \frac{\textrm{d} \theta}{\textrm{d} s} \left(R{\bf e}_\theta + R \alpha {\bf E}_3 \right) = {\bf e}_t. \end{equation*}

Because {\bf e}_t is a unit vector, we infer that \frac{\textrm{d} \theta}{\textrm{d} s} = \pm \frac{1}{R\sqrt{1 + \alpha^2}}. When \frac{\textrm{d} \theta}{\textrm{d} s} > 0, we find that the Frenet triad basis vectors have the representations

(27)   \begin{eqnarray*} && {\bf e}_t = \frac{1}{\sqrt{1 + \alpha^2}} \left( {\bf e}_\theta + \alpha {\bf E}_3\right), \\ \\[0.10in] && {\bf e}_n = - {\bf e}_r, \\ \\[0.10in] && {\bf e}_b = \frac{1}{\sqrt{1 + \alpha^2}} \left({\bf E}_3 - \alpha {\bf e}_\theta\right). \end{eqnarray*}

Alternatively, when \frac{\textrm{d} \theta}{\textrm{d} s} < 0,

(28)   \begin{eqnarray*} && {\bf e}_t = \frac{- 1}{\sqrt{1 + \alpha^2}} \left( {\bf e}_\theta + \alpha {\bf E}_3\right), \\ \\[0.10in] && {\bf e}_n = - {\bf e}_r, \\ \\[0.10in] && {\bf e}_b = \frac{- 1}{\sqrt{1 + \alpha^2}} \left( {\bf E}_3 - \alpha {\bf e}_\theta\right). \end{eqnarray*}

In terms of the pitch angle \gamma, the relations (27) and (28) simplify to the single prescription

(29)   \begin{eqnarray*} && {\bf e}_t = \cos(\gamma){\bf e}_\theta + \sin(\gamma){\bf E}_3, \\ \\ && {\bf e}_n = - {\bf e}_r, \\ \\ && {\bf e}_b = - \sin(\gamma){\bf e}_\theta + \cos(\gamma){\bf E}_3. \end{eqnarray*}

Some straightforward calculations of the derivatives of the Frenet triad basis vectors and the use of the chain rule reveal that the curvature

(30)   \begin{equation*} \kappa = R \left( \frac{\textrm{d} \theta}{\textrm{d} s} \right)^2 = \frac{1}{R\left(1 + \alpha^2\right)} = \frac{1}{R} \cos^2\left(\gamma\right), \end{equation*}

and the geometric torsion

(31)   \begin{equation*} \tau = \frac{\alpha}{R\left(1 + \alpha^2\right)} = \alpha \kappa = \frac{1}{R} \cos\left(\gamma\right)\sin\left(\gamma\right), \end{equation*}

both of which are constant. A helix is the only curve with constant curvature and constant non-zero torsion. For one segment of the helix of length \ell = \frac{2\pi R}{\cos(\gamma)}, the total torsion

(32)   \begin{equation*} T_\textrm{w} = \frac{1}{2\pi} \int^\ell_0 \tau \textrm{d} s = \frac{\tau \ell}{2 \pi} = \sin(\gamma) = \frac{\alpha}{\sqrt{1 + \alpha^2}} . \end{equation*}

The Darboux vector for the helix has the interesting representations

(33)   \begin{equation*} \bomega_{\textrm{SF}} = \kappa \left( \alpha {\bf e}_t + {\bf e}_b\right) = \frac{\textrm{d} \theta}{\textrm{d} s} {\bf E}_3. \end{equation*}

With the help of (29), the rotation tensor {\bf Q}_{\textrm{SF}}(s) associated with the Frenet triad for a circular helix can be decomposed into the product of a pair of rotations:

(34)   \begin{equation*} {\bf Q}_{\textrm{SF}}(s) = {\bf L}\left(\gamma, \, {\bf e}_r\left(s\right)\right){\bf L}\left( \theta(s) - \theta(0), \, {\bf E}_3\right), \end{equation*}

where the function {\bf L}(\phi, \, {\bf a}) denotes a rotation through a counterclockwise angle \phi about an axis {\bf a}. Finally, we note that a circle of radius R can be obtained from the helix by setting \alpha = 0 in the previous developments. In this case, the curvature \kappa = \frac{1}{R}, and the geometric torsion and total torsion are both zero: \tau = T_\textrm{w} = 0. Additionally, from (27) and (28), the Frenet triad reduces to the basis vectors for cylindrical polar coordinates, as illustrated in Figure 4: {\bf e}_t = \pm {\bf e}_{\theta}, {\bf e}_n = - {\bf e}_r, and {\bf e}_b = \pm {\bf E}_3.

Figure 4. Evolution of the Frenet triad  \{{\bf e}_t, \, {\bf e}_n, \, {\bf e}_b\} for a circle of radius  R. Depending on the sign of  \frac{\textrm{d} \theta}{\textrm{d} s}, the unit binormal vector  {\bf e}_b = \pm {\bf E}_3.

Material curves

In our upcoming discussion of the kinematics of rods, we model the centerline of a rod as a one-dimensional continuum of material points. This continuum is known as a material curve. Based on the works of Green, Naghdi, Wenner, and Laws [9, 10, 11, 12], the concept of a material curve \mathcal{L} embedded in three-dimensional Euclidean space \euclid^3 is described in the following manner. Referring to Figure 5, the current configuration of the material curve, \mathcal{C}, is defined by the vector-valued function {\bf r} = {\bf r} \left(\xi, \, t\right). Here, \xi is a convected coordinate along the material curve that uniquely identifies material points of \mathcal{L}; that is, even as the material curve moves in space, the coordinate \xi associated with a material point remains the same. The vector {\bf r} is the position vector of a material point of \mathcal{L} with respect to a fixed origin. A fixed reference configuration of the material curve, \mathcal{C}_0, is defined by the vector field {\bf R} = {\bf R}(\xi). For convenience, we assume that \xi is the arc-length parameter of the material curve’s fixed reference configuration.

Figure 5. A material curve  \mathcal{L} in its reference and current configurations,  \mathcal{C}_0 and  \mathcal{C}, respectively. The curve has length  \ell in its reference configuration. The material points  \xi = 0 and  \xi = \ell are labeled  \circ and  \bigtriangleup, respectively.

Kirchhoff’s rod theory

To develop a basic model of a rod, its centerline, represented by a material curve, must be capable of resisting bending and torsion. To incorporate these characteristics, additional structure on the concept of a material curve is needed. The added structure we use is a deformable set of vectors associated with each material point of the centerline. These vectors are known as directors, and a material curve with directors is referred to as a directed curve. The terminology “director” was introduced in a seminal paper [13] by Ericksen and Truesdell in 1958. One of the main contributions of Truesdell and Ericksen’s work was a critical examination of the Cosserat brothers’ work on deformable media [14, 15]. The paper [13] helped inspire a renaissance in works on rod and shell theories that continues to this day [9, 16, 17]. From a historical perspective, Kirchhoff (1824-1887) first presented a rod theory capable of modeling bending and torsion in 1859 [3]. Although his treatment differs considerably from the modern treatments that we follow here, the rod theory we present is known as Kirchhoff rod theory in recognition of his contribution. In the early years of the 20th century, the Cosserat brothers presented a formulation of Kirchhoff’s rod theory in [14, 15] using what we now call directors, and so the rod theory we discuss is often considered an example of a Cosserat rod theory.

The directed curve

A directed curve is a material curve for which a set of directors is defined at each material point. The directors represent kinematical quantities pertaining to the cross-sections of a rod-like body, and thus a directed curve serves as a model for such a body. Following the seminal work of the Cosserat brothers in the early 1900s, a directed curve is often known as a Cosserat curve. Referring to Figure 6, a directed curve in its current configuration is defined by the position vector {\bf r} = {\bf r}(\xi, \, t) locating a material point on the curve relative to a fixed origin by the convected coordinate \xi, and by the two directors {\bf d}_\alpha = {\bf d}_\alpha(\xi, \, t) (\alpha = 1, \, 2). A fixed reference configuration of the directed curve is defined by the position vector {\bf R} = {\bf R}(\xi) and by the referential directors {\bf D}_\alpha = {\bf D}_\alpha(\xi), where \xi is assumed to be the arc-length parameter of the curve in its reference configuration. In many treatments of rods, it is common to denote \frac{\partial{\bf r}}{\partial\xi} by {\bf d}_3. This is especially the case with the works of Green, Naghdi, Wenner, and Rubin [10, 11, 18]. Similarly, \frac{\partial{\bf R}}{\partial\xi} is denoted by {\bf D}_3. Alternatively, following Antman [16], one may define {\bf d}_3 = {\bf d}_1\times{\bf d}_2 and {\bf D}_3 = {\bf D}_1\times{\bf D}_2. For Kirchhoff’s rod theory, these two prescriptions are equivalent.

Figure 6. A directed (or Cosserat) curve in its reference and current configurations. In its current configuration, the curve’s directors are denoted by  \mathbf{d}_1 and  \mathbf{d}_2, while the referential directors are  \mathbf{D}_1 and  \mathbf{D}_2.

Kinematics

In Kirchhoff’s rod theory, the directed curve representing a rod-like body’s centerline is assumed to be inextensible, which implies

(35)   \begin{equation*} \frac{\partial{\bf r}}{\partial\xi} = {\bf e}_t = {\bf d}_3. \end{equation*}

Furthermore, we assume the body’s cross-sections remain plane and retain their orientation relative to the centerline. To render this assumption into a mathematically tractable form, we let {\bf D}_1, {\bf D}_2, and {\bf D}_3 = \frac{\partial{\bf R}}{\partial\xi} define a right-handed and orthonormal triad at each material point \xi in the reference configuration. If we also denote a fixed right-handed basis for \euclid^3 by \{ {\bf E}_1, \, {\bf E}_2, \, {\bf E}_3\}, then we can define a rotation tensor {\bf P}_0 that relates the bases \{ {\bf E}_1, \, {\bf E}_2, \, {\bf E}_3\} and \{ {\bf D}_1, \, {\bf D}_2, \, {\bf D}_3\}. That is,

(36)   \begin{equation*} {\bf D}_i = {\bf P}_0{\bf E}_i \quad (i = 1, \, 2, \, 3). \end{equation*}

For many reference configurations, we can conveniently choose {\bf D}_i = {\bf E}_i such that {\bf P}_0 = {\bf I}; some exceptions are for what Love [19] refers to as initially curved rods. Under a motion of the directed curve, the vectors {\bf d}_i retain their relative orientation and magnitude. Consequently,

(37)   \begin{equation*} {\bf d}_i = {\bf P} {\bf D}_i , \end{equation*}

where {\bf P} is a rotation tensor. The most popular method to parameterize {\bf P} is to use the Euler angles. For example, Love [19] used a set of 3-2-3 Euler angles to parameterize {\bf P}. Recently, Maddocks and some of his co-workers used unit quaternions (or the Euler-Rodrigues symmetric parameters) to parameterize this rotation [20, 21]. Regardless of how {\bf P} is parameterized, notice that

(38)   \begin{eqnarray*} && {\bf d}_\alpha = {\bf P}{\bf P}_0 {\bf E}_\alpha, \\ \\[0.10in] && \frac{\partial {\bf r}}{\partial \xi} = {\bf P}\frac{\partial {\bf R}}{\partial \xi}. \end{eqnarray*}

These two equations define the so-called Kirchhoff’s constraints on the rod. Utilizing the identity \bepsilon \left[{\bf Q}{\bf B}{\bf Q}^T\right] = \mbox{det}\left({\bf Q}\right) {\bf Q}\left( \bepsilon \left[{\bf B}\right]\right) for any orthogonal tensor {\bf Q} and any skew-symmetric tensor {\bf B}, the constraints (38) imply that

(39)   \begin{equation*} \frac{\partial {\bf d}_{\alpha}}{\partial \xi} = \left({\bf P}\left(\bnu + \bnu_0\right)\right)\times{\bf d}_\alpha, \end{equation*}

where \bnu and \bnu_0 are axial vectors of skew-symmetric tensors {\bf K} and {\bf K}_0, respectively:

(40)   \begin{eqnarray*} && \bnu = \mbox{ax}\left({\bf K}\right) = \mbox{ax}\left( {\bf P}^T \frac{\partial {\bf P}}{\partial \xi} \right) , \\ \\ \\ && \bnu_0 = \mbox{ax} \left({\bf K}_0\right) = \mbox{ax} \left( \frac{\partial {\bf P}_0}{\partial \xi} {\bf P}_0^T \right) . \end{eqnarray*}

Strains

The components of the vector \bnu with respect to the basis \left\{{\bf D}_1, \, {\bf D}_2, \, {\bf D}_3 \right\} in the reference configuration define three strain measures:

(41)   \begin{equation*} \bnu = \nu_1{\bf D}_1 + \nu_{2} {\bf D}_2 + \nu_{3} {\bf D}_3. \end{equation*}

The strains \nu_1 and \nu_2, commonly referred to as curvatures, correspond to bending resistance of the centerline in two directions, while \nu_3 is known as the torsion. This torsion should not be confused with the geometric torsion \tau defined earlier. Alternative strain measures are given by the components of the vector \bnu + \bnu_0 with respect to \left\{{\bf D}_1, \, {\bf D}_2, \, {\bf D}_3 \right\}:

(42)   \begin{equation*} \bnu + \bnu_0 = \left(\nu_1 + \nu_{01}\right){\bf D}_1 + \left(\nu_2 + \nu_{02}\right) {\bf D}_2 + \left(\nu_3 + \nu_{03}\right){\bf D}_3 , \end{equation*}

where \nu_{0i} are termed intrinsic strains. The strains given by (42) could be considered those relative to a reference configuration in which the rod’s centerline is straight and the directors are constant. Examples of using \bnu as a strain measure can be found in applications such as modeling of double-stranded DNA (e.g., see [22, 23, 24]). The effects of the bending and torsional strains on a rod are illustrated in Figure 7. For the reference configuration shown in this figure, the intrinsic strains \nu_{0i} = 0.

Figure 7. Illustration of the effects of the bending and torsional strains featured in Kirchhoff’s rod theory. The reference configuration of a straight rod is depicted in (a); the rod is subject to pure torsion  (\nu_1 = 0, \ \nu_2 = 0, \ \nu_3 \neq 0 ) in (b); and in (c), the rod is in pure bending  (\nu_1 \neq 0, \ \nu_2 = 0, \ \nu_3 = 0 ) such that the material curve that models the rod forms a circle.

Helical plies

As an illustrative example, we consider the helically plied structures shown in Figure 8. Suppose we wish to model the helical plies and the cylindrical structure about which they are wrapped using a rod theory. We are then faced with the issue of a choice of reference configuration for such a structure and, with that, a choice of the referential directors {\bf D}_1 and {\bf D}_2. Here, we consider one possible choice where these directors “follow” the helical strand. We emphasize that {\bf R}\left(\xi\right) for all three structures shown in Figure 8 is the position vector of material points located on the centerline of the cylindrical body.

Figure 8. Examples of reference configurations of a cylindrical body with helical plies wound on its lateral surface: (a) a single-ply, (b) a two-ply, and (c) a triple-ply. The two-ply is similar to double-stranded DNA. Also shown in this figure are samples of the referential directors  {\bf D}_1 and  {\bf D}_2 associated with the Cosserat rod that is used to model the body; the centerline of the cylindrical body is represented by the material curve  \mathcal{L}. This figure is adapted from [23].

We start by choosing a fixed unit vector {\bf E}_3 that is parallel to the longitudinal axis of the cylindrical body in its reference configuration such that the position vector

(43)   \begin{equation*} {\bf R}(\xi) = \xi {\bf E}_3 = \xi {\bf D}_3. \end{equation*}

Next, we denote the constant radius and pitch angle of the helical strand by R > 0 and \beta \neq 0, respectively. The equation of the helical strand’s centerline may be parameterized using a cylindrical polar angle \phi:

(44)   \begin{equation*} {\bf x} = R \left(\cos(\phi){\bf E}_1 + \sin\left(\phi\right){\bf E}_2\right) + R \tan\left(\beta\right) \phi {\bf E}_3. \end{equation*}

If we choose the material coordinate \xi to be the arc-length parameter for the reference configuration such that \xi = 0 when \phi = 0, then the two coordinates are related by

(45)   \begin{equation*} \xi = R \tan\left(\beta\right) \phi , \end{equation*}

and so the referential directors have the prescription

(46)   \begin{eqnarray*} && {\bf D}_1\left(\xi\right) = \cos\left(\frac{\xi}{R \tan\left(\beta\right)}\right) {\bf E}_1 + \sin\left(\frac{\xi}{R \tan\left(\beta\right)}\right){\bf E}_2, \\ \\ \\ && {\bf D}_2\left(\xi\right) = - \sin\left(\frac{\xi}{R \tan\left(\beta\right)}\right) {\bf E}_1 + \cos\left(\frac{\xi}{R \tan\left(\beta\right)}\right){\bf E}_2. \end{eqnarray*}

Consequently, we can express the representation (44) for \mathbf{x} in the compact form

(47)   \begin{equation*} {\bf x}\left(\xi\right) = R {\bf D}_1\left(\xi\right) +{\bf R}(\xi). \end{equation*}

From (46), it is straightforward to see that the rotation tensor {\bf P}_0 = \sum_{k \, \, = \, 1}^3{\bf D}_k\otimes{\bf E}_k corresponds to a rotation about {\bf E}_3 through an angle \frac{\xi}{R \tan\left(\beta\right)}. As a result, the axial vector of the skew-symmetric tensor \frac{\partial {\bf P}_0}{\partial \xi} {\bf P}^T_0 is given by

(48)   \begin{equation*} \bnu_0 = \nu_0{\bf E}_3 = \frac{\cot\left(\beta\right)}{R} {\bf E}_3. \end{equation*}

Thus, the reference configuration of the cylindrical body has a constant pretwist \nu_0 = \frac{\cot\left(\beta\right)}{R}. This pretwist is often known as the intrinsic twist and can be compared to the geometric torsion \tau = \frac{1}{R} \cos\left(\beta\right)\sin\left(\beta\right) of a helical space curve and the total torsion T_\textrm{w} = \sin\left(\beta\right) of a turn of this helix that we discussed earlier.

Twist, torsion, and tortuosity

The material curve associated with a directed curve can, at each instant, be associated with a space curve. It is natural to question how the strain measure \bnu + \bnu_0 and the rotation tensor {\bf P}{\bf P}_0 for the directed curve are related to the Darboux vector \bomega_\textrm{SF} and its corresponding rotation tensor {\bf Q}_\textrm{SF} for the space curve. We turn to discussing these matters. Referring to Figure 9, recall that the orthonormal and right-handed Frenet triad \left\{ {\bf e}_t, \, {\bf e}_n, \, {\bf e}_b \right\} associated with the space curve satisfies the Serret-Frenet relations (9), yielding the Darboux vector \bomega_{\textrm{SF}} = \tau {\bf e}_t + \kappa {\bf e}_b in terms of the space curve’s geometric torsion \tau and curvature \kappa.

Figure 9. Evolution of the Frenet triad  \{{\bf e}_t, \, {\bf e}_n, \, {\bf e}_b\} and the directors  {\bf d}_1 and  {\bf d}_2 along a directed curve.

Now, the directors associated with the directed curve remain mutually perpendicular to the centerline, establishing an orthonormal, right-handed triad \left\{ {\bf d}_1, \, {\bf d}_2, \, {\bf d}_3 = {\bf e}_t \right\} for the curve. Consequently, the unit normal vector {\bf e}_n and the unit binormal vector {\bf e}_b of the Frenet triad are coplanar with the directors {\bf d}_1 and {\bf d}_2. We can therefore surmise that the two triads are related by a rotation about {\bf d}_3 = {\bf e}_t:

(49)   \begin{eqnarray*} && {\bf d}_1 = {\bf Q}_t{\bf e}_n, \\ \\ && {\bf d}_2 = {\bf Q}_t{\bf e}_b. \end{eqnarray*}

The corresponding rotation tensor {\bf Q}_t has various representations:

(50)   \begin{eqnarray*} {\bf Q}_t \!\!\!\!\! &=& \!\!\!\!\! {\bf d}_1\otimes{\bf e}_n + {\bf d}_2\otimes{\bf e}_b + {\bf d}_3\otimes{\bf e}_t \\[0.075in] &=& \!\!\!\!\! \cos(\phi)\left({\bf d}_1\otimes{\bf d}_1 + {\bf d}_2\otimes{\bf d}_2\right) + \sin(\phi)\left({\bf d}_2\otimes{\bf d}_1 - {\bf d}_1\otimes{\bf d}_2\right) + {\bf d}_3\otimes{\bf e}_t \hspace{1in} \scalebox{0.001}{\textrm{\textcolor{white}{.}}} \\[0.075in] &=& \!\!\!\!\! \cos(\phi)\left({\bf e}_n\otimes{\bf e}_n + {\bf e}_b\otimes{\bf e}_b\right) + \sin(\phi)\left({\bf e}_b\otimes{\bf e}_b - {\bf e}_n\otimes{\bf e}_b\right) + {\bf d}_3\otimes{\bf e}_t, \end{eqnarray*}

where \phi is a counterclockwise angle of rotation. Several results now follow immediately. First, the rotation tensor

(51)   \begin{equation*} {\bf P}{\bf P}_0 = {\bf Q}_t{\bf Q}_\textrm{SF}, \end{equation*}

where

(52)   \begin{equation*} {\bf Q}_\textrm{SF} = {\bf e}_t\otimes{\bf E}_3 + {\bf e}_n\otimes{\bf E}_1 + {\bf e}_b\otimes{\bf E}_2. \end{equation*}

Here, \left\{ {\bf E}_1, \, {\bf E}_2, \, {\bf E}_3 \right\} denotes a fixed orthonormal, right-handed basis. Second,

(53)   \begin{equation*} {\bf P}\left( \bnu + \bnu_0\right) = \bomega_\textrm{SF} + \frac{\partial \phi}{\partial s}{\bf e}_t, \end{equation*}

and therefore the strain measures \nu_i (i = 1, \, 2, \, 3) and \nu_{0i} for the directed curve are related to the geometric torsion \tau and curvature \kappa of the space curve by

(54)   \begin{eqnarray*} && \nu_3 + \nu_{03} = \tau + \frac{\partial \phi}{\partial s}, \\ \\[0.10in] && \nu_1 + \nu_{01} = \kappa \sin(\phi), \\ \\ && \nu_2 + \nu_{02} = \kappa \cos(\phi). \end{eqnarray*}

In the context of the differential geometry of curves on surfaces, the relations (54) are known as Bonnet’s theorem and Meusnier’s theorem (see Kreyszig [5] and, for additional references, [25]). For Kirchhoff’s rod theory, these relations were first written down by Love in 1892 (see Section 253 of [19]). It is a common error to assume that \tau = \nu_3 + \nu_{03}. We follow Love and call \phi the angle of twist. To interpret \phi in this manner, imagine a straight rod (\tau = 0) and twist it between its ends. Possibly to avoid confusion with the geometric torsion \tau, Love defines the tortuosity \Sigma = \frac{1}{\tau} and uses \Sigma exclusively. We refer the interested reader to the wonderful paper [26] by Antman and Alexander for further discussion on this topic.

References

  1. Serret, J.-A., Sur quelques formules relatives à la théorie des courbes à double courbure, Journal de Mathématiques Pures et Appliquées 16 193-207 (1851).
  2. Frenet, J.-F., Sur quelques propriétés des courbes à double courbure, Journal de Mathématiques Pures et Appliquées 17 437-447 (1852).
  3. Kirchhoff, G., Über des gleichgewicht und die Bewegung eines unendlich dünnen elastichen Stabes, Journal für die reine und angewandte Mathematik 56 285-313 (1859).
  4. Darboux, G., Leçons sur la Théorie Générale des Surfaces et les Applications Géométriques du Calcul Infinitesimal, Parts 1-4, Gauthier-Villars, Paris (1887-1896).
  5. Kreyszig, E., Differential Geometry, revised and reprinted ed., Toronto University Press, Toronto (1964). Reprinted by Dover Publications in 1991.
  6. Oprea, J., Differential Geometry and its Applications, 2nd ed., Pearson/Prentice Hall, Upper Saddle River, NJ (2004).
  7. Spivak, M., A Comprehensive Introduction to Differential Geometry, Volume 2, 2nd ed., Publish or Perish, Berkeley (1979).
  8. O’Reilly, O. M., Engineering Dynamics: A Primer, 2nd ed., Springer-Verlag, New York (2008).
  9. Naghdi, P. M., Finite deformation of elastic rods and shells, in: D. E. Carlson, R. T. Shield (eds.), Proceedings of the IUTAM Symposium on Finite Elasticity, pp. 47-103. Martinus Nijhoff, The Hague (1982).
  10. Green, A. E., Naghdi, P. M., and Wenner, M. L., On the theory of rods. I. Derivations from three-dimensional equations, Proceedings of the Royal Society A 337(1611) 451-483 (1974).
  11. Green, A. E., Naghdi, P. M., and Wenner, M. L., On the theory of rods. II. Developments by direct approach, Proceedings of the Royal Society A 337(1611) 485-507 (1974).
  12. Green, A. E., and Laws, N., Remarks on the theory of rods, Journal of Elasticity 3(3) 179-184 (1973).
  13. Ericksen, J. L., and Truesdell, C., Exact theory of stress and strain in rods and shells, Archive for Rational Mechanics and Analysis 1(1) 295-323 (1958).
  14. Cosserat, E., and Cosserat, F., Sur la statique de la ligne déformable, Compte Rendus de l’Académie des Sciences, Paris 145 1409-1412 (1907).
  15. Cosserat, E., and Cosserat, F., Théorie des Corps Déformables, A. Hermann et Fils, Paris (1909).
  16. Antman, S. S., Nonlinear Problems of Elasticity, 2nd ed., Springer-Verlag, New York (2005).
  17. Naghdi, P. M., The theory of shells and plates, in: S. Flügge, C. A. Truesdell (eds.), Handbuch der Physik, Vol. VIa/2, pp. 425-640, Springer, New York (1972).
  18. Rubin, M. B., Cosserat Theories: Shells, Rods and Points, Kluwer Academic Press, Dordrecht (2000).
  19. Love, A. E. H., A Treatise on the Mathematical Theory of Elasticity, 4th ed., Cambridge University Press, Cambridge (1927).
  20. Dichmann, D. J., Li, Y., and Maddocks, J. H., Hamiltonian formulations and symmetries in rod mechanics, in: J. P. Mesirov, K. Schulten, D. W. Sumners (eds.), Mathematical Approaches to Biomolecular Structure and Dynamics, pp. 71-113, Springer-Verlag, New York (1996).
  21. Kehrbaum, S., and Maddocks, J. H., Elastic rods, rigid bodies, quaternions and the last quadrature, Philosophical Transactions of the Royal Society A 355(1732) 2117-2136 (1997).
  22. Healey, T. J., Material symmetry and chirality in nonlinearly elastic rods, Mathematics and Mechanics of Solids 7(4) 405-420 (2002).
  23. Lauderdale, T. A., and O’Reilly, O. M., On the restrictions imposed by non-affine material symmetry groups for elastic rods: Application to helical substructures, European Journal of Mechanics 26(4) 701-711 (2007).
  24. Marko, J. F., and Siggia, E. D., Bending and twisting elasticity of DNA, Macromolecules 27(4) 981-988 (1994). Errata for this paper were published on pg. 4820 of Vol. 29.
  25. O’Reilly, O. M., On the computation of relative rotations and geometric phases in the motions of rigid bodies, ASME Journal of Applied Mechanics 64(4) 969-974 (1997).
  26. Alexander, J. C., and Antman, S. S., The ambiguous twist of Love, Quarterly of Applied Mathematics 40(1) 83–92 (1982).