Generalized reciprocity theorems for infinitesimal deformations superimposed upon finite deformations of rods: the plane problem

The present paper is concernedwith bending of extensible and shear-deformable rods,wherewe treat the case of plane deformations. Our emphasis is laid upon infinitesimal deformations superimposed upon finite deformations. In this framework, we focus on deriving generalized reciprocity theorems that must hold with respect to the superimposed infinitesimal deformations, when produced by different types of loadings. Besides imposed distributed and concentrated force and moment loadings, we consider the action of eigenstrains. Moreover, we also take into account the presence of jumps in kinematic entities, particularly concerning jumps of axial rod positions and of cross-sectional rotations. We first derive local and global universal reciprocity relations, which must hold for rods made of any material. The global universal relation then is specialized to the case of hyperelastic rods. From the latter relation, we prove that reciprocity theorems known from the linear theory, such as the theorems of Betti, Maxwell and Land, as well as Maysel’s formula, formally do hold for infinitesimal deformations superimposed upon finite deformations. As a formulation that appears to be novel also with respect to the purely linear theory, we present a combination of Maysel’s formula with the theory of Land, the former dealing with eigenstrains, the latter involving jumps in kinematic entities.


Introduction
The present paper is concerned with bending of extensible and shear-deformable rods. For a comprehensive presentation of the theory of rods, see the book on the mechanics of deformable solids [1] by V.V. Eliseev, to whom our contribution is devoted in grateful memory.
Subsequently, we treat the case of plane deformations of rods, the finite strain theory of which dates back to Reissner [2]. Our emphasis is laid upon infinitesimal deformations superimposed upon finite deformations. For fundamentals of this type of superposition, see Knops and Wilkes [3]. In the present contribution, we focus on deriving generalized reciprocity theorems that must hold with respect to the superimposed infinitesimal deformations, when produced by different loadings. Besides imposed distributed and concentrated force and moment loadings, we also consider the action of eigenstrains. Moreover, we also take into account the presence of jumps in kinematic quantities, particularly concerning jumps of axial rod positions and of cross-sectional rotations.
Our contribution is organized as follows. In Sect. 2, we present the local relations of equilibrium for the finite strain case, as well as the local incremental equilibrium relations that hold for superposed infinitesimal deformations. Following and extending the procedure for rods rigid in shear presented by DaDeppo [4], the incremental equilibrium relations are derived as time rate forms of the equilibrium relations for the finite strain case. In Sect. 3, we derive a local virtual work relation by multiplying the local incremental equilibrium relations with proper infinitesimal virtual deformation rates. This relation is further simplified by introducing H. Irschik (B) Institute of Technical Mechanics, Johannes Kepler University of Linz, Linz, Austria E-mail: hans.irschik@jku.at the rates of generalized strain measures that were introduced by Reissner [2]. Local and global universal reciprocity relations are eventually derived in Sect. 4. We talk about a universal relation, when it does hold for any constitutive relation of the rod under consideration. In our derivation, we consider the virtual deformation quantities to be produced by a set of virtual internal and external forces and moments, which forms an equilibrium system, likewise to the original incremental problem. We are thus allowed to interchange original and virtual quantities in the local virtual work relation derived before in Sect. 3. This yields an adjoint virtual work relation. Subtracting the two virtual work relations provides the local universal reciprocity relation. The global form of the latter follows by integration along the rod segment under consideration in the reference configuration. In doing so, we also take into account the presence of concentrated external forces and external moments, which correspond to equilibrium relations of jump of the internal forces and moments. Besides of these static jump relations, we also consider kinematic jumps to be present in the infinitesimal superimposed deformations, namely jumps of the position vector or of the cross-sectional rotation at a certain location. Jumps may be present in both, the original and the virtual problem context. The static and kinematic relations of jump enter the universal global reciprocity relation, since integration must be performed piecewise, due to the presence of jumps in the local virtual work relations.
In Sect. 5, we specialize the universal global reciprocity relation derived in Sect. 4 to the case of hyperelastic rods, where we involve the generalized strain measures of Reissner [2], as well as proper integrability relations. We particularly allow the presence of incremental superimposed eigenstrains in the original problem. The generalized reciprocity theorem so obtained includes not only distributed forces and moments as well as jumps in forces and moments, but also eigenstrains and jumps in the kinematic relations. As special forms of this generalized formulation, various reciprocity relations known from the purely linear theory are proved to hold also for infinitesimal deformations superimposed upon finite deformations of shear-deformable and extensible rods, such as the Betti and Maxwell reciprocity theorems, see DaDeppo [4] for rods rigid in shear, and see Ziegler [5] for the linear theory. We further deal with an extension of a reciprocity theorem that involves jumps in the kinematic entities, and which in the purely linear context dates back to Land, see Kurrer [6] for a contemporary exposition. Concerning eigenstrains, we derive an extension of Maysel's formula of the linear theory of thermoelasticity, see Ziegler and Irschik [7] for the linear theory. As a topic which appears to be novel also with respect to the purely linear theory, we derive a reciprocity formulation that connects superimposed jumps of the virtual kinematic entities with the presence of superimposed eigenstrains in the original problem, and which thus allows to compute eigenstrain-induced internal forces and moments, when the solution due to the superimposed jumps in the kinematic entities is known. In Sect. 6, we finally present numerical and analytic examples for the correctness of the latter formulation.

Relations of equilibrium and incremental equilibrium relations
In the course of the present paper, we assume that the deformation of the rod takes place in the (x, z)-plane of a global Cartesian coordinate system with unit base vectors e x , e y , and e z , such that e y is perpendicular to the plane of deformation. The curved axial coordinate of the rod in the undeformed reference configuration is denoted as s 0 . The rod may have a curved axis already in the reference configuration. The axial coordinate of the rod in the deformed configuration is called s, see Fig. 1. Studying equilibrium of a rod element of differential axial length ds, but formulating any entities as functions of the reference coordinate s 0 , one obtains the local relations of equilibrium for the deformed configuration of the rod as, see [1,2], The derivative of some entity f with respect to s 0 is denoted by a prime: The force vector R represents the resultant of the inner forces acting upon a cross section of the rod in the deformed configuration, the position vector of the axis point in this cross section being denoted as r , see Fig. 1.
The moment vector M is the resultant couple of the inner forces with respect to the axis point. Note that the deformed configuration may represent a large deformation, possibly accompanied by large displacements from the reference configuration. For the sake of clearness, we will denote this deformed configuration as   In the intermediate configuration, we introduce a local cross-sectional coordinate system with unit base vectors e 1 and e 2 , the origin of which is situated at the axis point, and which stems from a positive rotation about the y axis of the global system by the angle ϕ, such that e 3 = e x sin ϕ + e z cos ϕ.
The corresponding decomposition of the resultants of the inner forces in the local system is written as Following Reissner [2], the e 1 -direction is taken perpendicular to the cross section in the deformed configuration, such that e 3 is tangential to the latter, where the angle ϕ is interpreted as the cross-sectional rotation, see Fig. 1. Hence, we call N and Q normal and shear force, respectively, and M is the bending moment. An analogous decomposition for the imposed forces and couples reads: We now formulate the local equilibrium conditions, Eqs. (1) and (2), in the local coordinate system. From Eqs.
(4) and (5), the derivatives of the local base vectors with respect to s 0 follow as The axial deformation gradient vector r introduced in Eq. (2) is tangential to the axis, The axial stretch is denoted as Λ = ds/ds 0 , and χ is the shear angle, see (1) and (2): We now consider the situation of a quasi-static loading that slowly evolves with time t. Taking into account an infinitesimal deformation with infinitesimal displacements superimposed upon the intermediate configuration, we are seeking for the rate form of Eqs. (13) and (14). Recall that such a superposition makes practical sense, since we have required the intermediate configuration to be infinitesimally superstable [3]. Following considerations by Da Deppo for rods rigid in shear [4], the rate form can be obtained from the equilibrium conditions in Eqs. (13) and (14) by differentiation with respect to the parameter t. Subsequently, the derivative of some entity f with respect to t, holding the place s 0 fixed, is denoted by a superimposed dot: By tacitly skipping dt, one then is allowed to interpret this rateḟ directly as the increment of f that corresponds to the superimposed infinitesimal deformation under consideration, being expressed as a function of the place s 0 in the reference configuration. Analogous to Eqs. (10) and (11), we now havė Using the latter relations, the rate forms of Eqs. (13) and (14) can be formulated aṡ For subsequent use, we note the increments of the resultant internal force and its space-wise derivative aṡ  18) and (19), in this special case indeed yields the well-known equilibrium conditions of the linear theory of rods, see [1] and Ziegler [5].

Local virtual work relations
Since we intend to derive reciprocity relations later on, we here provide local virtual work relations fitting to the rate equilibrium conditions stated in Eqs. (18) and (19). This is performed by a scalar multiplication of Eq.
(18) with a virtual displacement increment. In the following, virtual increments are indicated by underlined symbols, and they are understood as functions of the axial coordinate s 0 in the reference configuration again. These underlined virtual entities are assumed to correspond to an infinitesimal test deformation superimposed upon the same intermediate configuration introduced above, such that we are allowed to use superimposed dots also for the virtual entities. Hence, we obtain the local virtual work relation for the forces in the forṁ Now assume, for the moment being, that the entities are smooth enough at the place s 0 , such that the product rule of differentiation can be applied to Eq. (21): The increment of the virtual axial deformation gradient vector r follows by time-wise differentiation of Eq.
Substituting Eqs. (20) and (23) into (22) yields Similarly, multiplying the incremental equilibrium relation for the moment, Eq. (19), with a virtual rate of the rotation angleφ yields Adding Eqs. (24) and (25), we obtain the local relation of virtual work of the incremental external and internal forces and moments done with respect to the virtual incremental test deformations as This relation can be further simplified by utilizing generalized strain measures that were introduced by Reissner [2]. These measures consist of a generalized extensional strain, a generalized shear strain, and a generalized bending strain, Substituting the deformation rates in Eqs. (27) and (28) We now consider the underlined virtual deformation quantities to correspond to a set of virtual internal and external forces and moments. The virtual forces and moments are designated by corresponding underlined symbols, and they are considered to produce a state of equilibrium, too, likewise to the original forces and moments. We are thus allowed to interchange the roles of the underlined and non-underlined quantities in Eq.
The relation in Eq. (31) represents the local virtual work of the underlined virtual forces and moments done upon the non-underlined incremental deformations, which we denote as original ones in the following, while underlined quantities will be called virtual ones. Subtracting Eq. (31) from Eq. (30) yields the following universal local reciprocity relation: We call this relation as universal, since it must hold for any constitutive relation of the rod under consideration. Note that Eq. (31) also nicely reflects the fact that Reissner's generalized strain measures, see Eqs. (27)-(29), are work conjugate to the internal forces and moments. We now integrate over the rod segment under consideration, over the region 0 ≤ s 0 ≤ l 0 in the reference configuration. Assuming that the original incremental forces, moments, and deformations may suffer jumps at the place s 0 = a 0 , while the virtual entities may jump at s 0 = a 0 , integration of the bracketed terms at both sides of Eq. (32), the differentials of which are indicated by primes, must be performed piecewise. Eventually, the following global universal reciprocity relation including jumps of static and/or kinematic quantities is obtained: Qγ +Ṅε +Ṁφ ds 0 .
Jumps are indicated by double brackets, In deriving Eq. (33), it has been assumed with little loss of generality, but for the sake of brevity that This relation is trivially satisfied for classical homogeneous boundary conditions at the rod ends, e.g. at a clamped end, kinematic quantities do vanish. Moreover, note that the amounts of the jumps in Eq. (33) usually can be taken as known, because they correspond to the following static and kinematic relations of jump, respectively: ˙ r The index imp in Eqs. (36)-(39) refers to an imposed quantity.

Global reciprocity relations for hyperelastic rod
Following Reissner [2], the constitutive relations of an elastic rod can be formulated as functions of the generalized strains stated in Eqs. (27)-(29), from which their rate forms follow by differentiation with respect to time: In Eqs. (40)-(42), we also have introduced generalized eigenstrains ε e , γ e , and ϕ e , as well as their rates, these eigenstrains being taken as known in the following. For the sake of simplicity, we assume forces and moments to depend on the respective generalized eigenstrains only, see the last terms in Eqs. (40)-(42). Now, for hyperelastic rods, certain necessary and sufficient integrability conditions must be satisfied in order that the existence of a potential function can be guaranteed. These conditions read: For rods rigid in shear, where γ = 0 andγ = 0, see DaDeppo [4]. Relations analogous to Eqs. (40)-(43) also hold for the virtual (underlined) quantities; however, we subsequently consider eigenstrains to be absent in the virtual problem: Using Eqs. (40)-(44), one obtainṡ Substituting Eq. (45) into the universal global relation of reciprocity, Eq. (33), and considering Eq. (35), yields the following generalized reciprocity relation for hyperelastic rods in the presence of rates of eigenstrain in the original problem: As will be discussed subsequently, this relation includes and extends reciprocity relations known from the linear theory of rods to the case of infinitesimal deformations superimposed upon possibly large, but infinitesimally superstable deformations of rods. A main result of the present paper is that the form of these reciprocity relations is the same in both, the linear theory of rods (small deformations from an undeformed reference configuration) and in the theory of small deformations from an intermediate state with a possibly large deformation with respect to the reference state. In the latter case, the static and kinematic rates of course do depend on the deformation of the intermediate state.
To the best knowledge of the present author, the combined presence of kinematic jumps and eigenstrains in reciprocity relations yields a novel formulation also with respect to the linear theory of hyperelastic rods.

Betti's reciprocal theorem
Assuming eigenstrains as well as any jumps to be absent, from Eq. (46) we obtain the celebrated reciprocity theorem of Betti, see [5] for the linear theory, and Da Deppo [4] for infinitesimal deformations superimposed on finite deformations: The relation stated in Eq. (46) extends the result of DaDeppo [4] for rods rigid in shear to shear-deformable rods.

Reciprocity theorem dating back to Land
We now study the presence of jumps of the static and kinematic entities in the absence of eigenstrains and distributed forces and moments in Eq. (46). One possible case is the presence of static entities suffering jumps in the original problem and kinematic entities that jump in the virtual one. This yields: In the linear theory, such a relation is known as the reciprocity theorem of Land. We refer to the book of Kurrent [6], who deserves the merit of bringing back this theorem to the attention of a wider audience.

Extension of Maysel's formula
Consider the case of jumps being present in the static entities of the virtual problem, and of eigenstrains being present in the original problem in Eq. (46). One then obtains: Practically, this means that specific kinematic entities, located at s 0 = a 0 and produced by eigenstrains in the original problem, can be obtained by integration along the rod, once the solution of the virtual problem is known due to jumps of the work-conjugated static quantities, also applied at s 0 = a 0 . This technique is known as Maysel's formula in the linear theory of thermal eigenstrains, see Ziegler and Irschik [7]. See Irschik et al. [8] for the finite strain context, and Irschik and Pichler [9] for small vibrations superimposed upon large static deformations of piezoelastic bodies.
This allows computing internal forces and moments due to the eigenstrains, when the solution of the virtual problem due to the jumps of the kinetic entities is known. To the best knowledge of the present author, Eq.
(51) represents a novel extension of Maysel's procedure, not only with respect to the case of infinitesimal deformations superimposed upon a finite strain, but also in the framework of the purely linear theory.

Example problem
For an exemplary confirmation of the validity of our above results, we consider the case of an initially straight beam, the axial coordinate in the undeformed configuration being denoted as s 0 = x, where 0 ≤ x ≤ l 0 = l. The rod is taken as simply supported and axially movable at its left end, in x = 0, while it is clamped and fixed at the right end, in x = l. For the intermediate configuration, we consider an external moment M 0 and an external axial force R x0 to be applied at x = 0, see Fig. 2. We now assume imposed distributed forces and moments to be absent, where we first formulate relations for the intermediate configuration. The local equilibrium relations, Eqs. (13) and (14), now read We need to introduce constitutive hyperelastic relations, which, for the sake of simplicity, are taken as linear functions of Reissner's generalized strains, Eqs. (27)-(29), see also [1,2], In accordance with a frequent notation in structural mechanics, EI denotes the bending stiffness, EA is the tensile stiffness, and kGA stands for the shear stiffness, k being the shear factor, and E denoting Young's modulus. The cross-sectional area is A, and I denotes the cross-sectional moment of inertia.
The position vector of a point of the axis in the deformed configuration can be written as r = (x + u) e x + w e z , the displacement gradient vector following to r = 1 + u e x + w e z . The axial displacement is denoted by u, and w is the deflection. In order to close the problem, the following kinematic relations are considered, which follow by representing the deformation gradient vector in the form of Eq. (12) in the global coordinate system, using Eqs. (4) and (5), and resolving for the displacement gradients: In order to eliminate the shear angle χ from the relations in Eqs. (54), (58) and (59), we utilize Eqs. (55) and (56), which yields We thus have arrived at a system of six first-order nonlinear ordinary differential equations, Eqs.
The results enter the corresponding linear boundary-value problem for the superimposed infinitesimal deformations as known parameters, which however generally are non-constant.
For the original problem, the corresponding linear equilibrium relations follow in the absence of distributed imposed incremental forces and moments as, see Eqs. (13), (14), and (59): Having eliminated the shear angle and its infinitesimal increment, the kinematic relations become, see also Eqs. (61) and (62):u We particularly wish to obtain exemplary evidence for the correctness of the extended Maysel's formula, Eq. (51). Subsequently, for the sake of brevity, we consider the presence of a bending-type generalized incremental eigenstrainφ e only. The corresponding linear constitutive bending relation is taken aṡ The relations in Eqs. (65)-(70) represent a system of six first-order ordinary differential equations for the six incremental quantitiesu,ẇ,φ,Ṅ ,Q andṀ that are superimposed upon the intermediate configuration in the original problem, due to the generalized eigenstrainφ e , sometimes also denoted as an imposed curvature.
In order to exemplarily confirm the validity of Eq. (51), we consider the clamped end of the rod, i.e. we let a 0 → l, at which place we wish to compute the incremental clamping momentṀ x=l due toφ e . This means that we have to apply an imposed virtual negative cross-sectional rotation just left to the clamping end, in order to make φ Note that, despite the virtual problem (the underlined quantities) is considered to represent an infinitesimal deformation, a unit jump of the cross-sectional rotation angle has been introduced. This is admissible due to the linearity of the problem at hand and is customary in linear structural mechanics, see e.g. [6]. The extended Maysel's formula, Eq. (51), now specializes to Assuming that the constitutive behaviour is linear also for the intermediate configuration, and that, for the moment being, bending eigenstrains would be present, we find were used, e.g. concerning a selection of integration procedures. It turned out that, as long as the procedures did converge, results were obtained with a small computational effort, "just at fingertips", despite the problems at hand appear to be comparatively involved. In our numerical computations, differences between the numerical outcomes of the original problem and of Eq. (77), the latter involving the virtual problem, remained comparatively small. For instance, we considered a thick rod of circular cross section with diameter d = 0.2 m, length l = 1 m, and Young's modulus E = 100 GPa. Bending stiffness then is E I = 2.5 πMPa, tensile stiffness is E A = πMPa, and we set kG A = 0.5 πMPa. For the boundary conditions of the intermediate problem, Eq. (63), the end moment was taken as M 0 = 10π 2 KNm, and an axial compressive force R x0 was applied, which was 95% of the classical Euler buckling force for the problem at hand, not taking into account shear and extension, see e.g. [5]. When the diameter of the beam was decreased to d = 0.01 m, this difference increased to 0.4203%. The following explanation has been provided: As already mentioned, the above Land-type extension of Maysel's formula in the context of an intermediate state with a large deformation, Eq. (51), appears to represent a novel contribution to the eigenstrain theory. Extended numerical studies, concerning initially curved beams and eigenstrain loadings that are more complex than the constant bending-type one acting on an initially straight beam, which has been treated in the above exemplary numerical study with respect to Eq. (77), will be published in a forthcoming paper.
Note, however, that analytic checks of the more general Eq. (76) are possible via exemplary results from the literature on initially straight beams with non-constant eigenstrain loadings, when the undeformed configuration and the intermediate configuration do coincide, i.e. in the purely linear case, assuming M 0 = 0 and R x0 = 0. Consider again the case of a simply supported-clamped beam. Assume that the imposed curvature in the original problem is a box-type temperature difference between the upper and lower side of the beam,