- About this Journal
- Abstracting and Indexing
- Aims and Scope
- Article Processing Charges
- Articles in Press
- Author Guidelines
- Bibliographic Information
- Citations to this Journal
- Contact Information
- Editorial Board
- Editorial Workflow
- Free eTOC Alerts
- Publication Ethics
- Reviewers Acknowledgment
- Submit a Manuscript
- Subscription Information
- Table of Contents
International Journal of Mathematics and Mathematical Sciences
Volume 2011 (2011), Article ID 838924, 22 pages
Exact Solutions of Nonlinear Equation of Rod Deflections Involving the Lauricella Hypergeometric Functions
1Via Negroli, 6, 20136 Milan, Italy
2Dipartimento di Matematica per le Scienze Economiche e Sociali, Viale Filopanti, 5, 40126 Bologna, Italy
Received 4 December 2010; Revised 29 March 2011; Accepted 16 June 2011
Academic Editor: Kenneth Berenhaut
Copyright © 2011 Giovanni Mingari Scarpello and Daniele Ritelli. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
The stress induced in a loaded beam will not exceed some threshold, but also its maximum deflection, as for all the elastic systems, will be controlled. Nevertheless, the linear beam theory fails to describe the large deflections; highly flexible linear elements, namely, rods, typically found in aerospace or oil applications, may experience large displacements—but small strains, for not leaving the field of linear elasticity—so that geometric nonlinearities become significant. In this article, we provide analytical solutions to large deflections problem of a straight, cantilevered rod under different coplanar loadings. Our researches are led by means of the elliptic integrals, but the main achievement concerns the Lauricella hypergeometric functions use for solving elasticity problems. Each of our analytic solutions has been individually validated by comparison with other tools, so that it can be used in turn as a benchmark, that is, for testing other methods based on the finite elements approximation.
1.1. A Literature's Overview
As early as 1691, Jakob Bernoulli proposed to find out the deformed centerline (planar “elastica”) of a thin, homogeneous, straight, and flexible rod under a force applied at its end, and ignoring its weight. Studying the geometry of its deflection, he established, Curvatura laminae elasticae (1694), the differential equation of the centerline, say , of the rod ab appenso pondere curvata to be Not succeeding in solving it, he was led to look for its arclength, arriving at the same elliptic integral of lemniscate rectification. Daniel Bernoulli, Johann's son and Jakob's nephew, obtained in 1738 subspecie integralis the potential energy stored in a bent rod; in October 1742 he wrote to Leonard Euler the rod should attain such a shape as to minimize the function of squared curvature. Accordingly, Euler dealt elastica as an isoperimetric problem of calculus of variations, De curvis elasticis (1744), arriving at elastica's differential equation and identifying nine shapes of the curve, all equilibrium states of minimal energy. Anyway, Colin Maclaurin had realized the elastica's connection to elliptic integrals, see A treatise on fluxions (1742), § 927 entitled “The construction of the elastic curve, and of other Figures, by the rectification of the conic sections”. In 1757, Euler authored a paper, Sur la force des colonnes, concerning the buckling of columns again, where the critical load is approached through a simplified expression to the elastica's curvature, coming again over the problem again in 1770 and 1775; for all this historical concern, see [1, 2]. Some analytic solutions to elastic planar curves through elliptic integrals of I and II kind can be read at [3–5]. Almost half a century ago, the analytic treatment of flexible rods inspired a nice book, , which collects many problems ordered according to the constraint type, and all solved by the elliptic integrals of I and II kind, while the III kind and Theta functions appear there marginally in the 6th Chapter concerning the three dimensional deformations. In order to fix just right from the start, the true concept of rod will be recalled from  as follows:A rod is something of a hybrid, that is, a mathematical curve made up of material points. It has no cross section, yet it has stiffness. It can weigh nothing or it can be heavy. We can twist it, bend it, stretch it and shake it, but we cannot break it, it is completely elastic. Of course, it is a mathematical object, and does not exist in the physical world; yet it finds wide application in structural (and biological) mechanics: columns, struts, cables, thread, and DNA have all been modeled with rod theory.A rod which offers no resistance to bending is called a string: it is completely flexible. The string model finds applications in the study of long slender structural members such as cables, where an external load dominates the mechanics and the rods flexural stiffness is negligible.
The link between elliptic functions and elastica was always understood so close, to induce the elliptic functions to be plotted through suitable curved rubber rods (G. Greenhill: Graphical representation of the elliptic functions by means of a bent elastic beam, 1876), while the paper  is founded on the use of the Weierstrass function. For recent developments, see: [8–15]. Inextensibility and circular shape of undeformed rod are (numerically) analyzed in [16, 17] where the load is a uniform centrally directed force.
Our object is the analytical solutions to large-deflection problems of a planar slender cantilever under coplanar loads. The slender structures are those much longer in one direction than in the other two, and their deflection will be described by elastic rod theory if one is interested in phenomena on length scales much larger than the lateral dimensions of the structure. For Dirichlet problems relevant to large deformations, see [18, 19]. A simplified version of this theory, beam theory, suffices if only small deflections need to be considered; the majority of structural engineering applications are adequately modelled by a beam; bridges and buildings are not designed to suffer large deformations. But in other applications of structural mechanics, deflections may be large, meaning that nonlinear geometric effects must be taken into account; the full (geometrically nonlinear) rod theory is then necessary. Mathematically, a rod is modelled as a curve with effective mechanical properties such as bending and torsional stiffnesses, see . Although such problems belong to statics, they have a strong relation with dynamics; arclength along the rod plays a role similar to time in a dynamical system. In its simplest form, this analogy (Kirchhoff's Dynamic Analogy: Über das Gleichgewicht und die Bewegung eines unendlich dűnnen elastischen Stabes, 1859, see ) is seen between the deformations of an elastic strut and the motions of a swinging pendulum. A less known analogy links a twisted rod and a spinning top, and the precession of a top does correspond to a rod helical deformation. Engineering problems where rod theory is applied include the looping of ocean cables  and the buckling of oil well drilling strings, which nowadays can be over 5 km long. Furthermore, there is a deep reason for a great interest in applying the elastic rod theory to the supercoiling and packing of DNA molecules and proteins (such as collagen): as a matter of fact, due to its double helical nature, DNA is unusually stiff (in both, bending and torsion) for a polymer and, hence, well described by a rod. Finally, the study of strings and rods is of interest to electrodynamic space tethers1, generally credited of great potential for future space missions in order to bring down space debris (derelict spacecraft, etc.). Unlike the majority of tethers, which can be treated as cables, Space European Tether (SET, see ) is a prototype designed to withstand torsional and bending forces so that it needs to be modelled as a spinning rod (and not a cable), namely, a 3D-elastica2. All above provide enough motivations to implement exact solutions within the rod theory, see .The exact governing equations for the finite displacement rod theory become highly nonlinear, and, hence, it is very difficult to solve these equations analytically. Thus, for practical purposes, problems of this kind are usually solved by the methods based on the finite elements approximations. However, the closed form solutions for these problems are still important to the practical point that the accuracy of the approximate methods can be precisely evaluated by these solutions, to say nothing of the mathematical importance.
1.2. Aim of the Paper
It is well known why the elastic curve of a rod is essential, in fact not only the stresses, but also its maximum deflection will not exceed a fixed amount. Furthermore, (see Mattiasson ).When developing finite element methods, there is also a need for checking the results against an exact solution. In the case of finite element methods for geometrically nonlinear beam, frame and shell structures, elliptic integral solutions of large deflection beam and frame problems offer such exact solution.
Nevertheless, his approach is a bit obsolete nowadays because he computes the elliptic integrals numerically through the arithmetic-geometrical means.
We are going to analyze an elastic, thin, flexible rod, clamped at one end (cantilever) bent under different loadings whose plane coincides with that of the strained centerline. Equating the exact curvature in cartesian coordinates to the bending moment divided into the flexural rigidity, one is led to a second-order nonlinear and nonautonomous differential equation. Boundary conditions are given both on the clamped end (Cauchy problem), so that, for example, if the constraint is perfect, displacement and rotation have to be zero there. We discarded the illusion of a general system of loads and rods and restricted ourselves to only one configuration, charged by a single load and treating each single problem by means of the special functions of mathematical physics. The chosen loads are constant couple, evenly distributed load, shear at the tip, hydrostatic load, and sinusoidal load. All our force loads have to be understood as “dead load”, namely, always equal to itself and then not “follower” as are currently named the loads keeping a fixed inclination to the rod before and after the strain.
For one of the above loads, based on the rod's unextensibility, we computed in closed form the free tip coordinates, parametrized to the load itself. We use everywhere nondimensional quantities; the reader will see how advantageous is such a choice in producing “universal” relationships which (in elliptical and/or hypergeometric cloth) keep their formal validity in several situations.
1.3. Special Functions Tools
We assume that the reader knows the elliptic functions (see  for details), but let us recall hereinafter the less known Lauricella hypergeometric ones.
The hypergeometric series first appeared in the Wallis's Arithmetica infinitorum (1656) for and real parameters . The product of factors called Pochhammer symbol (or truncated factorial) allows to write as The first meaningful contributions about various representations are due to Euler3; nevertheless, he does not seem to have known the integral representation traditionally ascribed to him, but really due to A. M. Legendre4. For all this and for the Stirling contributions, see . The above integral relationship is true if and for , even if this bound can be discarded due to analytic continuation.
Several functions have been introduced in the 19th century for generalizing the hypergeometric functions to multiple variables, but we will mention only those introduced and investigated by Lauricella (1893) and Saran (1954). Among them our interest is focused on the hypergeometric function of variables (and parameters), see [27, 28], defined as with the hypergeometric series usual convergence requirements . If , the relevant integral representation theorem (IRT) provides allowing the analytic continuation to deprived of the cartesian -dimensional product of the interval with itself.
2. Bent Slender Rod: A Large Deflections 2D-Model
Let us tackle a -long, thin, and weightless rod whose end A is clamped and B free. We put at B the origin O of a cartesian reference frame with along the undeformed rod, from B to A; and downwards, normal at B to . Our basic assumptions are(A1)the rod is thin, initially straight, homogeneous, with uniform cross-section and uniform flexural stiffness , where is the Young modulus, and the cross-section moment of inertia about a “neutral” axis normal to the plane of bending and passing through the central line;(A2)the slender rod is charged by coplanar dead loads;(A3)linear constitutive load: the induced curvature is proportional via to the bending moment. The rod undergoes large deflections but small strains;(A4)the shear transverse deformation is ignored, and the rod is assumed to keep constant its previous length (inextensibility);(A5)a stationary strain field, by isostatic equilibrium of active loads and reactive forces, takes place, and, due to rest of static equilibrium, no rod element undergoes acceleration.
The key Bernoulli-Euler linear constitutive law connects flexural stiffness and bending moment to the consequent curvature of the deflected centerline . With our (traditional) reference frame , the curvature has always opposite sign to the bending moment where is the unknown elastica's cartesian equation and To (2.1) a useful dress can be given to be solved together to the perfect clamping boundary conditions Throughout all the paper, it will be minded that the Cauchy problem can be solved in such a way. Putting the solution of (2.5) is
In the sequel, we are going to use dimensionless variables , so that the previous relationship becomes whose solution is
2.1. Rod Inflected by a Constant Bending Couple at Its Tip
The clamped section is at , see Figure 1. In such a case, the equation to be integrated is (2.5), with , and (2.9) provides, in nondimensional form, where . The elastica is more readable reverting at a dimensional form half circle of centre and radius , see Figure 1, as it is well known.
2.2. The Horizontal Heavy Rod
The classic elastica theory started by the “lamina”, namely, a thin rod without weight but capable to withstand the flexure and loaded at its free end only. Chapter 5 of  on the specific problem of the horizontal cantilever under uniform load gives (page 173):This problem… has shown considerable resistance to a closed form solution and all investigations dealing with it have employed approximate methods only.
Accordingly, he obtains elastica coordinates as parametric functions of the arclength , stopping their Maclaurin expansions to .
Let us consider at Figure 2 a horizontal, flexible, homogeneous, heavy and cantilever; it will be inflected by the uniform evenly continuous load (weight for unity span length). Its elastica, according to the approximate beam theory, would consist of a fourth-order rational function of . Being we obtain, in nondimensional coordinates, using (2.9) and finally, passing to the variable :
By identification, we obtain the solution to our elastica through a Lauricella function of 3 variables , so that For being compliant with the IRT requirements, it shall be , and then . Otherwise, it will be necessary another method.
We got there three -ranges. The first is the very narrow one, where the small deflections assumption is true, first-order theory. By , we get which at , gives . For instance, if , then(i)if , the first-order approximation can be accepted;(ii)if by (IRT), the hypergeometric function is capable of describing our solution (2.14);(iii)if , the solution (2.14) will be numerically computed.
Formula (2.15), founded on the Lauricella hypergeometric functions, provides a very compact explicit solution to the elastica of a thin clamped heavy rod. We cannot give now the due importance to (2.15) due to the lack of software packages whose quickness and precision of computing are comparable to those available on Mathematica for the Gauss or Appell hypergeometric functions. We did invite Wolfram to consider such a necessity, but no forecasting is possible. As far as we are concerned, the integral (2.14) cannot be led to any other known special functions or to any combination of them.
2.3. A Tip-Sheared Horizontal Rod
We refer to Figure 3, where the force is transversal (in  is studied the inclined end load) but not a follower one; we have , so that, after defining the nondimensional , we obtain being the bending moment always greater than zero, so will be , and then it will be . For this load too, the hypergeometric treatment is possible. Passing from to , with , we obtain
We can see that, recognizing the Lauricella parameters compliant to the IRT requirements if . For this load, we have to the Lauricella's alternative of the elliptic integrals (not applicable to the heavy nor to the submerged rod), so that we will cover with much more details this approach. Investigating the elastica of a thin cantilever, the majority of authors (Landau and Lifchitz , Love , Mathieu , Burgatti , Panayotounakos and Theocaris [11, 12, 14], Chucheepsakul and Phungpaigram , Basoco , Frisch-Fay ) use the slope-arclength formulation they integrate a first time for getting as a function of . But and , so that they get, through a second integration, the elastica's and coordinates parametrized on . Only few develop a different parametric equations, for example, Saalschütz, in . Anyway a treatment leading to explicit sole cartesian coordinates does not seem to be ever established, and it is exactly what we are doing. An old reference, rather fit to Jakob Bernoulli (1694) and Euler (1744) formulae, is that of Kiepert (1912), who5 stops without reducing his elliptic integral to the Legendre standard form. It will be just the elliptic integral (2.17) our startup. First let us change variables in (2.17) passing from to , obtaining For the reader convenience, let us fix some notation. Let . Consider the definite elliptic integral Combining formulae given in : integral 234.16 page 76 and integral 340.01 page 205, introducing we obtain: where and are the I and III kind elliptic integrals. Now assuming that we have the identification which allows to give (2.17) via (2.20) by means of the integration formula where: being Formula (2.25) provides the dimensionless deflection to the position and to the load : its 3D-plot, via Mathematica, as shown in Figure 4.
From it, one could obtain the relevant contour plots cutting by planes of constant. The thing is shown in Figure 5, where the clamped end is on the left. The deflection curves are given for some values of . Each of them is compared with the dot curve coming from the approximate theory6. Mind that is ranged as , and has to be found from (3.12). Accordingly, each single -curve has been plotted between and . From the physical viewpoint, the lines'7 comparison shows at the low loads an agreement which is becoming worst for increasing loads when the approximate theory is progressively failing.
It is worth noting that the power series expansion of (2.25) drives back to the approximate expression for used in literature. Of course, we deeply exploited the computer algebra power for this task. First we compute the Maclaurin third-order expansion of (2.25) with respect to obtaining where ( is the quantity introduced in (2.27)) Now we expand the elliptic integrals of I and III kind with respect to , with initial point arriving at Moreover recalling the Maclaurin first-order expansion with respect to adding (2.31) and (2.32), we obtain eventually that is, the usual approximation.
Going again to formula (2.25), his peculiarities are: it is nondimensional; it is new, in closed form, exact, and explicit. Moreover, for being implemented, no previous transcendental equation has to be solved; it requires only to compute the dimensionless parameter holding all the rod physical data.
On the purpose, it will be observed that in , an eigenvalue problem is solved through a double boundary condition of no rotation at the ends of the rod. The authors get the coordinates of elastica parametrized to the external load and marked by the mode number. For our same load, they show four elastica each relevant to a given load, but referring to the static mode . Other curves are provided referring to the first dynamic mode, . Furthermore, the module of the Jacobi elliptic functions involved there, is on its own depending on the tip load via an equation to be solved previously. Nevertheless the (essential!) relationship keeps unrevealed, because the authors at page 740 write:The modulus of elliptic functions and the parameter play the role of integration constants and they are related to the force and the angle by boundary conditions for each case of rod bending.
In fact, they do not provide the computation at all, marking their plot solutions not in terms of (as it should be done), but of a sequence of values 0.5; 0.6; 0.85; see Figure 2, page 742 of , so that their solution seems to be not easily accessible from .
2.4. The Rod Inflected by Hydrostatic Pressure
We are looking for the shape of a slender and thin cantilever clamped at the bottom of a pool of calm water. The rod weight has been deliberately ignored, being partially balanced by the Archimedes buoyancy. The rod is at full length submerged, being obvious the adjustments for a rod leaning out the water or whose tip is under the sheet of water. We assume, as in Figure 6, a reference with the axis downwards and over the water surface. The hydrostatic head, being the liquid mass density and the gravity, will have a linearly deep-growing profile (Stevin), so that, integrating twice such a head, one arrives at the bending moment
By the usual identification, we obtain the solution to our elastica through the Lauricella function of three variables where , so that For being compliant with the IRT requirements, it will be , and then .
2.5. The Rod Loaded by a Sinusoidal Bending Moment
Assumption A3 implies that our systems are ruled by differential equations affected by a geometrical nonlinearity due to the nonlinear link between the curvature and . But , which is the consequence on the rod, is linked to its cause, namely, the bending moment , by a Linear Constitutive Law , being the flexural stiffness. Such a linearity implies that, by superposition, whenever some moment distribution is repeated periodically, we can make its Fourier expansion, so that the rod bending moment is broken out into the addition of an infinitely many harmonics sine/cosine. In such a way, one will get the elastica by superposition, of course retaining those harmonics deemed useful or necessary.
What above provides enough motivation for analyzing once for all, the effect of the sinusoidal load, whose a half wave only is taken, omitting all the possible variants or subcases. Seeing Figure 7, let be the continuous transverse load (whose highest value is ) over the unit of an -length flexible cantilever Integrating twice such a load, one arrives at the bending moment affecting the -section In such a way, introducing , with8 , following (2.9), we find If , so that By the integration Tables , entries 254.12 page 113 and 337.01 page 201, we can establish an integration formula, embedding (2.40) as a particular case: let and define then where is again the elliptic integral of third kind and of first kind. Being , in (2.40), we do the identification so that our elastica is where Figure 8, performed through Mathematica, provides the overlapping of the integral (2.39) numerically computed and the plot consequent to our compact smart formula (2.44).
3. How to Compute the Tip Position after the Strain
Considering the small deflections, the deformed centerline's curvature is approximately equated to the deflection's second derivative, assuming the rod deformation field to proceed along purely vertical displacements. If it is so, the free end will be mapped to , with , but having the elastica as horizontal projection, its full unstrained length, then the rod will undergo extension.
With large displacements all, this is not true any more, and a further unknown will appear, namely, the abscissa of the new position , where has been mapped to. Rod's behavior under strain can be modeled in several ways. For instance in  after integrating the elastica, the authors write five congruence relationships involving geometric quantities before and after the strain and then leave the rod free of extension.
Following a completely different method, in  some possible alternatives like rotational load, load always perpendicular to the tangent at the tip and so on, are taken into account.
Anyway, the simplest criterion is the inextensibility assumption A4; the rod keeps its originary length9. For the horizontal weightless cantilever inflected by a tip shear (see Figure 3), the elastica has been provided by means of the elliptic integrals, see formula (2.25). Wishing to compute the tip deflection , the approximate beam theory would answer , but with large deflection, we have first to detect what is done assuming the unextensibility We see that is the sole unknown, having been previously found , so that, after , we will compute . let us provide the full treatment for such a case. By plugging in (3.1), after the change , we get Our unknown appears in the upper integration limit; in order to solve (3.3), we need to establish an existence-uniqueness Lemma.
Lemma 3.1. Consider the transcendental -equation where Defining then the equation (3.4) has a unique real root if and only if and such a root is given by where is the Jacobi sine amplitude of module .
Proof. Let . Observe that the function is strictly increasing, and that (3.6) ensures existence of solution of equation (3.4). In order to prove the Lemma, we will refer to , formula 234.00, page 74, to establish that with the and meanings previously introduced. The very last step follows from the inversion: which delivers from (3.9), after equating the right hand side to .
We can then use the Lemma to solve (3.3) via the identification where it has been assumed10 . Minding that and , it should be seen that (3.6) provides an upper limit to the dimensionless load solving numerically a transcendental inequality. To such a limit, a physical meaning can be found as the highest allowable load beyond which the linear constitutive law (curvature versus bending couple) is not true any more.
Solving to , we find out that Therefore, increasing intensities11, and then , the successive positions of the free tip, that is, the components of the displacement will be described as functions of at Figures 9 and 10.
Five elasticae have been analytically solved of a thin cantilever rod bent by different coplanar “dead” loadings12, namely, whose direction remains unchanged during the deformation of the rod. They are concentrated couple, tip shear, evenly continuous load, hydrostatic load, and sinusoidal load. Each of them describes a single physical situation for which we obtain the large deflections in closed form through elliptic and/or hypergeometric functions. Some care has been put (tip sheared rod) in providing analytically, via the isometric assumption, the free end position after the strain, parametrizing its coordinates as a function of the load. Some of our solutions, specially that to the heavy rod by the Lauricella functions, have enriched the literature on the elastica.
The authors are indebted to their friend Aldo Scimone who drew some Figures of this paper; they hereby take opportunity for thanking him warmly. They wish also to thank the anonymous referee whose valuable criticism improved their article.
- A space tether is a long cable used to connect spacecraft to other orbiting bodies such as space stations, boosters, payload, and so forth in order to transfer energy and momentum.
- The problem of 3D-elastica is very old, but out of our purpose. We will say only that it first appeared in Lagrange's Mécanique Analytique 2nd edition (1811). J. Binet (Mémoire sur l'intégration des équations de la courbe élastique à double courbure, 1844) and M. Wantzel (Note sur l'intégration des équations de la courbe élastique à double courbure, 1844) arrived at the same differential system, vaguely formulated, and not affronted. A modern treatment of a particular case can be read in , with its specific bibliography.
- We quote three works: (a) De progressionibus transcendentibus, Op. omnia, S.1, vol. 28; (b) De curva hypergeometrica Op. omnia, S.1, vol. 16; (c) Institutiones Calculi integralis, 1769, vol. II.
- A.M. Legendre, Exercices de calcul intégral, II, quatriéme part, sect. 2, Paris 1811.
- In  putting , he gets the elastica through the integral “Die Kurve welche diese Gleichung entspricht, heißt die elastische Linie weil ein elastischer Stab der an dem Ende befestigt und an dem andern Ende belastet ist, diese Form annimmt.”
- Under transversal load at the cantilever tip, this problem in strength of materials is solved in cartesian coordinates following the linear approximation; the solution is found (e.g., see [37, page 305]) to be in dimensionless form.
- Notice that the plots are coming after Mathematica's computation of our explicit solution (2.25), but adding the map in order to provide the geometrical true configuration of the stressed rod.
- Inserting in some numerical values, for example, , one gets: , so that, unless of assuming unrealistic lengths, is essentially much less than unity.
- Only in some special cases such a tip path description is quick. For instance for a rod loaded by a fixed couple (Section 2.1), one gets .
- Having for instance , we get .
- Nothing to say, apart from the more strict link imposed by the transcendental inequality (linearity range), (3.12) requires .
- The contrary case is that of the follower load which is being unchanged in its relationship with the rod, for example, keeping its direction perpendicular—or inclined—to it during and after the strain.
- C. G. Fraser, “Mathematical technique and physical conception in Euler's investigation of the elastica,” Centaurus, vol. 34, no. 3, pp. 211–246, 1991.
- E. Benvenuto, La Scienza delle Costruzioni e il Suo Sviluppo Storico, Sansoni, Firenze, Italy, 1981.
- W. Schell, Theorie der Bewegung und der Kräfte, vol. 2, B.G. Teubner, Leipzig, Germany, 1880.
- P. Burgatti, Teoria Matematica dell'Elasticità, Zanichelli, Bologna, Italy, 1931.
- F. G. Tricomi, Funzioni Ellittiche, Zanichelli, Bologna, Italy, 1951.
- R. Frisch-Fay, Flexible Bars, Butterworths, London, UK, 1962.
- V. G. A. Goss, “The history of the planar elastica: insights into mechanics and scientific method,” Science and Education, vol. 18, no. 8, pp. 1057–1082, 2009.
- M. A. Basoco, “On the inflexional elastica,” The American Mathematical Monthly, vol. 48, pp. 303–309, 1941.
- S. Chucheepsakul and B. Phungpaigram, “Elliptic integral solutions of variable-arc-length elastica under an inclined follower force,” Zeitschrift für Angewandte Mathematik und Mechanik, vol. 84, no. 1, pp. 29–38, 2004.
- Y. Goto, T. Yoshimitsu, and M. Obata, “Elliptic integral solutions of plane elastica with axial and shear deformations,” International Journal of Solids and Structures, vol. 26, no. 4, pp. 375–390, 1990.
- D. E. Panayotounakos and P. S. Theocaris, “Exact solution of the non-linear differential equation concerning the elastic line of a straight rod due to terminal loading,” International Journal of Non-Linear Mechanics, vol. 17, no. 5-6, pp. 395–402, 1982.
- D. E. Panayotounakos and P. S. Theocaris, “Exact solution for an approximate differential equation of a straight bar under conditions of a non-linear equilibrium,” International Journal of Non-Linear Mechanics, vol. 21, no. 5, pp. 421–429, 1986.
- D. E. Panayotounakos and A. B. Sotiropoulos, “Exact parametric analytic solutions of the elastica ODEs for bars including effects of the transverse deformation,” International Journal of Non-Linear Mechanics, vol. 39, no. 10, pp. 1555–1570, 2004.
- D. E. Panayotounakos and P. S. Theocaris, “Analytic solutions for nonlinear differential equations describing the elastica of straight bars: theory,” Journal of the Franklin Institute, vol. 325, no. 5, pp. 621–633, 1988.
- T. J. Lardner, “A note on the elastica with large loads,” International Journal of Solids and Structures, vol. 21, no. 1, pp. 21–26, 1985.
- J. H. Jung, T. J. Kang, and K. W. Lee, “Elliptic integral solutions of the large deflection of a fiber cantilever with circular wavy crimp,” Textile Research Journal, vol. 73, no. 1, pp. 47–54, 2003.
- R. W. Dickey and J. J. Roseman, “Equilibria of the circular elastica under a uniform central force field,” Quarterly of Applied Mathematics, vol. 51, no. 2, pp. 201–216, 1993.
- G. Mingari Scarpello and D. Ritelli, “Elliptic integrals solution to elastica's boundary value problem of a rod bent by axial compression,” Journal of Analysis and Applications, vol. 5, no. 1, pp. 53–69, 2007.
- G. Mingari Scarpello and D. Ritelli, “Elliptic functions solution to exact curvature elastica of a thin cantilever under terminal loads,” Journal of Geometry and Symmetry in Physics, vol. 12, pp. 75–92, 2008.
- S. S. Antman, Nonlinear Problems of Elasticity, Springer, New York, NY, USA, 1995.
- A. E. H. Love, A Treatise on the Mathematical Theory of Elasticity, Dover, New York, NY, USA, 4th edition, 1944.
- J. Coyne, “Analysis of the formation and elimination of loops in twisted cable,” IEEE Journal of Oceanic Engineering, vol. 15, no. 2, pp. 72–83, 1990.
- V. V. Beletsky and E. M. Levin, Dynamics of Space Tether Systems, Advances in the Astronautical Sciences Volume 83, American Astronautical Society, 1993.
- K. Mattiasson, “Numerical results from large deflection beam and frame problems analysis by means of elliptic integrals,” International Journal for Numerical Methods in Engineering, vol. 17, no. 1, pp. 145–153, 1981.
- P. F. Byrd and M. D. Friedman, Handbook of Elliptic Integrals for Engineers and Scientists, Springer, New York, NY, USA, 2nd edition, 1971.
- J. Dutka, “The early history of the hypergeometric function,” Archive for History of Exact Sciences, vol. 31, no. 1, pp. 16–34, 1984.
- S. Saran, “Hypergeometric functions of three variables,” Ganita, vol. 5, pp. 77–91, 1954.
- G. Lauricella, “Sulle funzioni ipergeometriche a più variabili,” Rendiconti del Circolo Matematico di Palermo, vol. 7, no. 1, pp. 111–158, 1893.
- S. Navaee and R. E. Elling, “Equilibrium configurations of cantilever beams subjected to inclined end loads,” Journal of Applied Mechanics, vol. 59, no. 3, pp. 572–579, 1992.
- L. Landau and E. Lifchitz, Théorie de l'Élasticité, MIR, Moscow, Russia, 1967.
- E. Mathieu, Dynamique Analytique, Gauthier Villars, Paris, Farnce, 1876.
- L. Saalschütz, Der Belastete Stab unter Einwirkung einer Seitlichen Kraft, B.G. Teubner, Leipzig, Germany, 1880.
- Y. V. Zakharov and K. G. Okhotkin, “Nonlinear bending of thin elastic rods,” Journal of Applied Mechanics and Technical Physics, vol. 43, pp. 739–744, 2002.
- B. N. Rao and V. G. Rao, “On the large deflection of cantilever beams with end rotational load,” Zeitschrift für Angewandte Mathematik und Mechanik, vol. 66, no. 10, pp. 507–508, 1986.
- G. Mingari Scarpello and D. Ritelli, “Elliptic integral solutions of spatial elastica of a thin straight rod bent under concentrated terminal forces,” Meccanica, vol. 41, no. 5, pp. 519–527, 2006.
- L. Kiepert, Grundriss der Differential-und Integral-Rechnung II. Teil, Hannover Helwingsche Verlagsbuchhandlung, 1918.
- O. Belluzzi, Scienza delle Costruzioni, vol. 1, Zanichelli, Bologna, Italy, 1970.