The clamped-free rod under inclined end forces and transitions between equilibrium configurations

We investigate the problem of the straight, inextensible and unshearable clamped-free elastic rod subjected to an inclined end force. Exact analytic solutions representing all equilibrium configurations of the deformed rod are presented in elliptic integral form. Those exact solutions, for a given angle of inclination of the end force and number of inflection points, are characterised by two quantities; the end force and the elliptic modulus. Critical points are discussed and analytic conditions for determining their location are presented. Certain critical points where transitions between two equilibrium configurations whose numbers of inflection points differs by one are pointed out. Simple formulae for the total number of equilibrium configurations for a given end force are given. Applying arguments based on the elastic strain energy of the rod, we discuss scenarios where highly inflectional equilibrium configurations can transition to equilibrium configurations with fewer inflection points.


Introduction
Large deformations of clamped-free rods, also known as cantilevers, has been the subject of study by many researchers. For a complete and comprehensive historical review of this subject, we refer the readers to the works of Navaee and Elling [1] and Batista [2] and the references cited therein.
The clamped-free rod under a purely compressive end force succumbs to large post-buckling deformations once the rod has buckled under a critical force. The character of the large post-buckling deformations of the rod will depend on the mode n the rod has buckled into, where n = 1, 2, 3, . . ., see Timoshenko and Gere [3]. Those modes n reflect the number of inflection points present in an exact solution representing an equilibrium configuration for large post-buckling deformations and also provide a means for distinguishing or classifying those exact solutions.
This article revisits the problem of determining exact solutions for all equilibrium configurations of the clampedfree rod under inclined end forces. The inclined end force is provided to the rod in such a manner that even though the force increases, the angle of inclination of the force is kept constant. The angle of inclination of the end force is restricted to an interval which allows the discussion of the clamped-free rod under compression, the pure cantilever, and all other cases in between in a single unified approach. Although that is an old problem and many of the results are well known-see Navaee and Elling [1,4] and Batista [2,5]-this paper draws attention to certain aspects of the mechanics of the problem that have been overlooked. Firstly, transitions and jumps can occur at critical points between equilibrium configurations with different numbers of inflection points when increasing and decreasing the end force, i.e. loading and unloading the rod, respectively. Secondly, the manifestation of those transitions and jumps depends on which parameter is being controlled during loading and unloading the rod; for example, control of the end force (dead loading) leads to different transitions from control of end displacement (rigid loading) [6]. Thirdly, transitions and jumps between equilibrium configurations with different numbers of inflection points involve changes in the elastic strain energy of the rod. That observation offers a means of understanding which equilibrium configuration is selected when a rod undergoes a transition or jump. We apply that elastic strain energy-based argument to explain how a rod in a highly inflectional equilibrium configuration can transition to an equilibrium configuration with fewer inflection points, and ultimately the equilibrium configuration representing the straight rod, by releasing elastic strain energy.
In our analysis of the aforementioned transitions and jumps, we first determine exact solutions for all equilibrium configurations of the clamped-free rod subjected to an inclined end force in terms of elliptic integrals. Those exact solutions are presented in a form that makes their dependency on the number of inflection points explicit. We can then classify exact solutions representing equilibrium configurations by the number of inflection points, as for the case of the rod under pure compression. That is very useful for counting the number of equilibrium configurations for a given end force. We show that the exact solution representing each equilibrium configuration for a given angle of inclination of the end force and number of inflection points can be characterised by two quantities: the end force and the elliptic modulus (a function of the angle at the free end). Consequently, once we have determined those quantities, the shape of the deformed rod is reduced to a simple exercise in plotting points with elliptic integrals.
For each equilibrium configuration of the clamped-free rod under an inclined end force, at any inflection point, we give the slope or angle, the arc length parameter along the rod, and the coordinates. We additionally discuss critical points and present analytic conditions for determining their location. Furthermore, we show that the local force minima (critical points) can be used to define ranges for the values of inclined end force between which the number of inflectional equilibrium configurations can be counted. Simple formulae for the number of equilibrium configurations for a given end force are thereby determined. Finally, we use the results of our analysis to examine how transitions and jumps between equilibrium configurations with different numbers of inflection points can arise.

Mathematical model
We consider an elastic rod parameterised by its non-dimensional (normalised) arc length, 0 ≤ s ≤ 1. In the Cartesian coordinate frame with coordinates (x, y), the rod is clamped at the origin (s = 0) and free at the other end (s = 1). The rod is initially straight, inextensible and unshearable but can deform in the x y-plane when the free end is subjected to a non-dimensional force f inclined at angle α. The angle the tangent at any point (x(s), y(s)) on the deformed rod makes with the x-axis is given by θ(s), see Fig. 1. The (non-dimensional) force f is positive and the angle α is restricted to 0 ≤ α ≤ π . There is no need to consider values of α in π < α ≤ 2π , due to symmetry.
In the Cartesian coordinate frame with coordinates (u, v), which is rotated clockwise through angle α with respect to the (x, y) coordinate frame, as depicted in Fig. 1, we see that Furthermore, the angle the tangent at any point (u(s), v(s)) on the deformed rod makes with the u-axis is ψ(s) given by The (non-dimensional) ordinary differential equations describing the clamped-free rod, which may be derived from force and moment balance arguments applied to an infinitesimal element of deformed rod, are dv ds = sin(ψ).
The boundary conditions are Once ψ, u and v are known, it is easy to determine θ , x and y from Eqs. (1), (2) and (3). Apart from the convenience of not having to carry around all those factors of cos(α) and sin(α) in our calculations, another reason for working in the (u, v) coordinate frame is that the inclined end force f acts along the u-axis.

Exact solutions
When α = 0, we have the rod under pure compression; for α = π/2, we have the pure cantilever; and at α = π , we have the rod under pure tension. We will consider the case α = π first, since the rod is inextensible, and this case is the simplest to deal with. The case α = 0 is dealt with last, as this can be obtained by taking the α → 0 limit in the case 0 < α < π.

α = π
The inextensibility of the rod restricts the exact solutions available. The obvious solution is The above solution represents a rod that remains straight under all f ≥ 0, see Antman [7]. That case was not discussed by Navaee and Elling [1]; however, Batista mentions it briefly [2]. We mention it here for the sake of completeness, but will not discuss it any further.

0 < α < π
Integrating Eq. (4) subject to the boundary condition in Eq. (8), we find where, anticipating elliptic integrals, we introduce the elliptic modulus k given by and γ = θ (1). Note that we have taken f to be a constant during the integration procedure. Equilibrium configurations of the deformed rod may have n inflection points, where n = 1, 2, 3, . . .. At an inflection point dθ/ds = dψ/ds = 0, so from Eqs. (12) and (13), for an equilibrium configuration with n inflection points, the ith inflection point is at where s i is the value of s at the ith inflection point and i = 1, 2, 3, . . . , n. We have s 0 ≤ s 1 < s 2 < s 3 < · · · < s n . Note that s n = 1, and for notational convenience, we set s 0 = 0 and clarify that the latter is not an inflection point. The sign in Eq. (14) is chosen to ensure that the angle at the free end of the rod, where there is an inflection point (see Eq. (8)), is ψ n = ψ(s n ) = ψ(1) = γ + α. Solutions to Eq. (12) exist if k 2 − sin 2 (ψ/2) ≥ 0. That condition places restrictions on the allowed values of γ . Under the inclined end force f the rod deforms with either γ ≥ 0 or γ < 0. As f → ∞, in the former case, we will have γ → π − α, and in the latter case, γ → −(π + α). To explain this, observe when f = 0, Eq. (4) can be re-written as (1/ f )(d 2 ψ/ds 2 ) + sin(ψ) = 0. As f → ∞, the dominant term is sin(ψ). Thus, we get sin(ψ) = 0. For the inflectional solutions considered in this paper, we obtain ψ = ±π . That solution is valid near s = 1 since the boundary condition at s = 1, Eq. (8), holds for ψ = ±π . So we can set s = 1 in that solution and obtain γ + α = ±π as required. Note that ψ = ±π does not satisfy the boundary condition at s = 0, given by Eq. (7). At s = 0, we have a boundary layer, i.e. a concentration of curvature, which changes the orientation of the rod at s = 0 from ψ(0) = α to ψ(0) = ±π . See Audoly and Pomeau [8] for a discussion on the method behind deriving the solution ψ = ±π as f → ∞ and boundary layers. When 0 ≤ γ < π − α, k 2 − sin 2 (ψ/2) ≥ 0 always holds. When −(π + α) < γ < 0, k 2 − sin 2 (ψ/2) ≥ 0 only holds when γ ≤ −2α. Hence, the negative and positive values of γ must be restricted to the intervals −(π + α) < γ ≤ −2α and 0 ≤ γ < π − α, respectively. We refine those restrictions slightly to conform with describing equilibrium configurations by the number of inflection points n. We require Our derivations, based on an examination of plots of sin 2 (ψ/2), will not be elaborated on here since the conditions given by Eqs. (15) and (16) have been presented by Navaee and Elling [4]. The implications of Eqs. (15) and (16) for the elliptic modulus k are In order to completely determine ψ, u and v, introduce and follow the approach similar to the one outlined by Frisch-Fay [9] or Dym [10]. Since ψ/2 always lies between −(γ + α)/2 and (γ + α)/2, cos(ψ/2) > 0. Thus, we can differentiate Eq. (19) with respect to ψ and obtain Also, from Eq. (19), we have where φ is the elliptic argument. The value of φ at ψ = α is denoted by and the value of φ at ψ = ψ i is denoted by for i = 1, 2, 3, . . . , n. Combining Eq. (12) with the correct sign and Eq. (20), we find for for i = 1, 2, 3, . . . , n. It's a straightforward procedure to use Eqs. (5), (6), (19) and (21)-(24) to calculate, for for i = 1, 2, 3, . . . , n. K (k) and E(k) are the complete elliptic integrals of the first and second kind, respectively. F(φ, k) and E(φ, k) are the incomplete elliptic integrals of the first and second kind, respectively. See Gradshteyn and Ryzhik [11] for further details on elliptic integrals. The values of s, u and v at the ith inflection point are where u i = u(s i ), v i = v(s i ) and i = 1, 2, 3, . . . , n. The values of s, u and v at the nth (last) inflection point, where s = s n = 1, are particularly important so we state them explicitly where u n = u(s n ) = u (1) and v n = v(s n ) = v (1). A solution, ψ, u and v, describing an equilibrium configuration of the rod is completely determined once f and k are known. The rod is deformed by providing an inclined force f applied to the free end of the rod, i.e. we input values of f into Eq. (31) and numerically solve for corresponding values of the elliptic modulus k. However, it is much simpler to input the allowed values of k given by Eqs. (17) and (18) into Eq. (31) and obtain the corresponding values for f from which force-elliptic modulus ( f -k) plots can be obtained. Determining the x y-shape of the deformed rod is now a simple exercise in plotting points with elliptic integrals using Eqs. (1), (2) and (25)-(30). It may be noticed that for a given n, solving Eq. (31) for a specified f could result in more than one value for k, since k can be positive (when γ > 0) or negative (when γ < −2α). Thus, for a given n and a specified f , we expect to have non-uniqueness, i.e. multiple solutions. See Sect. 4.1 for a more detailed discussion of non-uniqueness.

α = 0
When α = 0, we have θ = ψ, x = u, y = v and the results given by Eqs. (25)-(33) above reduce to the solutions for deformations of the clamped-free rod from the nth buckled mode. As for the case with 0 < α < π, a solution θ , x and y, describing an equilibrium configuration of the rod is completely determined once f and k are known.
The allowed values of k are given by setting α = 0 in Eqs. (17) and (18). Explicitly, we find Input those values of k into Eq. (31), where we have already set α = 0, and obtain the corresponding values for f .

Properties of the equilibrium configurations
The rod is deformed by providing an inclined force f applied to the free end of the rod. In that case, we can interpret all physical quantities as being functions of f , for example, k = k( f ).
As mentioned above, the case α = π only has the solution given by Eq. (11) for all f ≥ 0. Therefore, we only consider the 0 < α < π and α = 0.

0 < α < π
If we take k = k( f ), the critical points on f -k plots are determined from dk/d f = 0 and d f /dk = 0, which can be easily found from Eq. (31). Taking dk/d f = 0, we find where the corresponding values of f denoted by f (n) α are (36) Those critical points are associated with local minima and maxima for k. Setting d f /dk = 0, gives where sgn(k) denotes the sign of k. Equation (37) can be solved for k = k min . Those critical points are associated with local minima for f . That can be seen by simply rotating the f -k plots in Fig. 2 such that the f -axes point vertically up. Batista presents a version of Eq. (37) in terms of Jacobi elliptic functions, but only solves for n = 1 and −1 < k < − sin(α/2) [5]. In that condition, unlike our condition given by Eq. (37), the dependency on the number of inflection points n is not explicit.
We illustrate our results for α = π/4, π/2 and 3π/4 with f -k plots in Fig. 2. In that figure, the vertical black dashed lines pass through the critical points ( f α , sin(α/2)) for n = 1, 3, 5, . . .. At those critical points, for n ≥ 2, transitions between equilibrium configurations with n and n − 1 inflection points are possible under increasing and decreasing f , i.e. loading and unloading the rod, respectively. For example, at ( f (2) α , − sin( α 2 )), marked as a black dot on each of the plots in Fig. 2, the rod may deform into an equilibrium configuration with either 1 inflection point, if f decreases, or 2 inflection points, if f increases.
As mentioned at the end of Sect. 3.2, we expect multiple solutions for a specified f , i.e. non-uniqueness. Consider solutions for a given n. When we specify f and look for solutions to Eq. (31), we will find solutions where k could be positive or negative, corresponding to deformations of the rod into equilibrium configurations with positive or negative end angle γ , respectively. On an intuitive level, that makes sense since we can deform the rod in two different directions (γ < −2α or γ > 0) and therefore we can reasonably expect multiple solutions for the same f . Even so, there are still two properties which affect non-uniqueness. Firstly, the asymmetry between deforming the rod with γ < −2α (k < − sin(α/2)) and γ > 0 (k > sin(α/2)) results in two different domains for f . To see that, compare plots for a given n with γ < −2α (k < − sin(α/2)) and γ > 0 (k > sin(α/2)) in Fig. 2. When the two domains for f overlap, we have multiple solutions (non-uniqueness), when there is no overlap, we have one solution (uniqueness). Secondly, the presence of local minima for f at ( f (n) min , k (n) min ) means that for f such that f If we denote the number of equilibrium configurations for a given f by N e ( f ), then for n = 1, 2, 3, . . .. An illustration of Eq.hus (39) is provided by Fig. 3.
With multiple solutions representing the equilibrium configurations of the rod for a given f , which solution is selected? According to Love, solutions with internal inflection points are unstable [12]. Therefore, only the solutions for n = 1 are stable. Batista also holds that view and adds that the solution for n = 1 is also unstable for the section of curve between points ( f (2) α , − sin(α/2)) and ( f (1) min , k (1) min ) on any of the plots in Fig. 2 [5]. Everywhere else, the solution for n = 1 is stable. If stability is the criteria for selecting a solution for a given f , then solutions with n = 1 are likely to be selected. However, equilibrium configurations for higher n = 2, 3, . . . although unstable, may be realised in the presence of suitable constraints applied at certain inflection points that suppress equilibrium configurations with lower n. In Bigoni et al, those constraints are referred to as "clamped-guided" [13].  Fig. 3 The x-y shapes for equilibrium configuration with 1 (red), 2 (green) and 3 (blue) inflection points for α = π/4, π/2 and 3π/4 and various values of f . (Color figure online) 4.2 α = 0 The critical force, denoted by f (n) c , at which the rod will buckle into the nth mode can be found by setting k = 0 in Eq. (31) where we have already set α = 0. We find  Let N e (n, f ) denote the number of equilibrium configurations for a given n = 1, 2, 3, . . . and specified f , then and, if we denote the number of equilibrium configurations for a given f by N e ( f ), then for n = 1, 2, 3, . . .. See Fig. 5 for an illustration of Eq. (44).

Transitions between equilibrium configurations
Increasing or decreasing the force f in general will increase or decrease the elastic strain energy of the rod, respectively. Therefore, the elastic strain energy of the rod is important in any discussion where the force f is varied. which after a short calculation becomes for n = 1, 2, 3, . . .. Consider the following scenarios: 1. The rod is in an equilibrium configuration with n = 1 and γ > 0 on the red curve denoted by B in Fig. 6. Increasing f will move to an equilibrium configuration higher up on the red curve and an increase in e, see Fig. 7. Decreasing f moves to an equilibrium configuration lower down on the red curve and a decrease in e.
Continually decreasing f will result in a straight rod at A.
2. The rod is in an equilibrium configuration with n = 1 and γ < −2α on the red curve denoted by C in Fig. 6. Increasing f will move to an equilibrium configuration higher up on the red curve and an increase in e. Decreasing f moves to an equilibrium configuration lower down on the red curve and a decrease in e.
Continually decreasing f will result in reaching the equilibrium configuration with f = f (1) min at D. Once at D, attempting to decrease f further will result in a jump to D 1 , where n = 1 and γ > 0. At D 1 , decreasing f proceeds as described in 1.
3. The rod is in an equilibrium configuration with n = 2 and γ < −2α on the green curve denoted by E in Fig. 6. Increasing f will move to an equilibrium configuration higher up on the green curve and an increase in e. Decreasing f moves to an equilibrium configuration lower down on the green curve and a decrease in e.
We can continue decreasing f to reach F, where f = f (2) α . Decreasing f further the equilibrium configuration transitions from one with n = 2 and γ = −2α to one with n = 1 and γ < −2α. The furthest we can decrease f to is f = f (1) min at D. Once at D, decreasing f proceeds as described in 2. According to Love [12] and Batista [5] the section of curve between D and C is stable. However, we believe that section cannot be reached by starting on the green curve and continuously decreasing f . Our reasoning based on always keeping e as small as possible will result in a path along the green curve that leads to F, then D, and finally proceeds as in 2. 4. The rod is in an equilibrium configuration with n = 2 and γ > 0 on the green curve denoted by G in Fig. 6.
Increasing f will move to an equilibrium configuration higher up on the green curve and an increase in e.
Decreasing f moves to an equilibrium configuration lower down on the green curve and a decrease in e.
Continuously decreasing f will result in reaching the equilibrium configuration with f = f (2) min at H. Once at H, attempting to decrease f further will result in a jump. This time we have 3 choices to jump to: H 1 , H 2 or H 3 . We argue that the process of decreasing f removes work (energy) from the system which must be reflected in a corresponding decrease in e. Since H 2 and H 3 have higher e than the equilibrium configuration at H, the jump will be to H 1 with a lower e than the equilibrium configuration at H, see Fig. 7. Thus, the rod jumps from the equilibrium configuration at H with n = 2 and γ > 0 to the equilibrium configuration at H 1 with n = 1 and γ > 0. From H 1 , decreasing f follows along the lines described in 1. 5. The rod is in an equilibrium configuration with n = 3 and γ > 0 on the blue curve denoted by I in Fig. 6.
Increasing f will move to an equilibrium configuration higher up on the blue curve and an increase in e, see Fig. 7. Decreasing f moves to an equilibrium configuration lower down on the blue curve and a decrease in e.
We can continue decreasing f to reach J, where f = f (3) α . Decreasing f further the equilibrium configuration transitions from one with n = 3 and γ = 0 to one with n = 2 and γ > 0. The furthest we can decrease f to is f = f (2) min at H. Once at H, decreasing f proceeds as described in 4. 6. The rod is in an equilibrium configuration with n = 3 and γ < −2α on the blue curve denoted by K in Fig. 6.
Increasing f will move to an equilibrium configuration higher up on the blue curve and an increase in e, see Fig. 7. Decreasing f moves to an equilibrium configuration lower down on the blue curve and a decrease in e. Continuously decreasing f will result in reaching the equilibrium configuration with f = f (3) min at L. Once at L, attempting to decrease f further will result in a jump. This time we have five choices to jump to: L 1 -L 5 . Since L 1 , L 2 and L 4 have greater e than L, we can discount jumps to them. Both L 3 and L 5 have less e than L. If a jump releases the least amount of e, then the rod jumps to the equilibrium configuration L 5 with n = 1 and γ < −2α. If a jump releases the greatest amount of e, then the rod jumps to an equilibrium configuration L 3 with n = 1 and γ > 0. From L 3 and L 5 , decreasing f follows along the lines described in 1 and 2, respectively. Thus, in this scenario 2 transition paths are possible. 7. Finally consider the rod in an equilibrium configuration at one of the local f minima, f (n) min . Take D for example where f = f (1) min in Fig. 6. Decreasing f will result in a jump to D 1 as described in 2. If f is increased the rod can move to some equilibrium configuration between D and C or D and F. The rod is more likely to choose an equilibrium configuration between D and F because this equilibrium configuration will have a lower e than an equilibrium configuration between D and C, see Fig. 7. The situation is almost like a bifurcation. Indeed, as α → 0, Eq. (42) tells us that we do have a bifurcation since the local f minima f Typical force-displacement ( f -d) plots for 0 < α < π are shown in Fig. 8. During the process of directly increasing and decreasing f (dead loading), all arguments for scenarios 1-7 apply to the f -d plots in Fig. 8. However, f could also be provided indirectly by a device or mechanism that acts on the rod through some end displacement parameter d. That manner of providing f is called rigid loading [6]. In that case, Eq. (47) is interpreted as a boundary condition for u(1) that sets a constraint on the displacement of the free end of the rod in the direction of f . To specify the corresponding boundary value problem, the additional differential equation d f ds = 0 is required, which ensures conservation of f for each displacement d. Under that constraint, the rod can behave differently during loading and unloading. For example, consider scenario 3. The rod is in an equilibrium configuration with n = 1 and γ < −2α on the red curve denoted by C in Fig. 8. Indirectly decreasing f by increasing d moves to an equilibrium configuration higher up on the red curve and a decrease in e. Continually decreasing f will result in reaching the equilibrium configuration with f = f (1) min at D. However, there is still scope for increasing d further until the local maximum (critical point) at M is reached. Once at M, any further attempt to increase d will result in a jump to A, the straight rod

Conclusion
The problem of the straight, inextensible and unshearable clamped-free elastic rod subjected to inclined end forces discussed in this article is well known [1,2,4]. Our contributions to the body of knowledge on this problem are as follows: 1. Exact solutions for all equilibrium configurations of the deformed rod in terms of elliptic integrals are given.
Our version of the exact solutions makes the role of the number of inflection points in an inflectional equilibrium configuration central to their description. We also present explicit formulae for the arc length parameter s and coordinates u, v at those inflection points. 2. A discussion of critical points on f -k plots has been presented and the locations of these critical points are determined. We identify critical points associated with local minima and maxima for k as places where transitions between equilibrium configurations with n and n − 1 inflection points occur for n ≥ 2. For critical points associated with local minima for f , a condition for locating the local force minimum f (n) min on the nth equilibrium configuration of the rod is presented. 3. General formulae are given for the number of equilibrium configurations for a given inclined end force f . 4. Applying arguments based on the elastic strain energy of the rod, we discuss scenarios where highly inflectional equilibrium configurations can transition to equilibrium configurations with fewer inflection points, and ultimately to the equilibrium configuration representing the straight rod. 5. It is pointed out that transitions between equilibrium configurations under rigid loading may differ from those under dead loading.