Research Papers

A Numerical Method to Model Dynamic Behavior of Thin Inextensible Elastic Rods in Three Dimensions

[+] Author and Article Information
Stephen Montgomery-Smith

Department of Mathematics,
University of Missouri,
Columbia, MO 65211
e-mail: stephen@missouri.edu

Weijun Huang

Mechanical and Aerospace Engineering,
University of Missouri,
Columbia, MO 65211
e-mail: whfy6@mail.missouri.edu

Some people might argue that a vector of vectors should really be called a tensor.

Contributed by the Design Engineering Division of ASME for publication in the JOURNAL OF COMPUTATIONAL AND NONLINEAR DYNAMICS. Manuscript received October 25, 2012; final manuscript received October 1, 2013; published online November 7, 2013. Assoc. Editor: Johannes Gerstmayr.

J. Comput. Nonlinear Dynam 9(1), 011015 (Nov 07, 2013) (10 pages) Paper No: CND-12-1182; doi: 10.1115/1.4025627 History: Received October 25, 2012; Revised October 01, 2013

Static equations for thin inextensible elastic rods, or elastica as they are sometimes called, have been studied since before the time of Euler. In this paper, we examine how to model the dynamic behavior of elastica. We present a fairly high speed, robust numerical scheme that uses (i) a space discretization that uses cubic splines, and (ii) a time discretization that preserves a discrete version of the Hamiltonian. A good choice of numerical scheme is important because these equations are very stiff; that is, most explicit numerical schemes will become unstable very quickly. The authors conducted this research anticipating describing the dynamic Kirchhoff problem, that is, the behavior of general springs that have natural curvature, and for which the equations take into account torsion of the rod.

Copyright © 2014 by ASME
Your Session has timed out. Please sign back in to continue.


Love, A. E. H., 1927, A Treaties on the Mathematical Theory of Elasticity, 4th ed., Dover, New York.
Montgomery-Smith, S. J., and Huang, W., 2012, “Dynamic Equations of Non-Helical Springs,” (to be published).
Levian, R., 2008, “The Elastica: A Mathematical History,” University of California, Berkeley, Technical Report No. UCB/EECS-2008-103.
Falk, R. S., and Xu, J.-M., 1995, “Convergence of a Second-Order Scheme for the Non-Linear Dynamical Equations of Elastic Rods,” SIAM J. Numer. Anal., 32(4), pp. 1185–1209. [CrossRef]
Caflisch, R. E., and Maddocks, J. H., 1984, “Nonlinear Dynamical Theory of the Elastica,” Proc. R. Soc. Edinburgh, 99A, pp. 1–23. [CrossRef]
Coleman, B. D., and Dill, E. H., 1992, “Flexure Waves in Elastic Rods,” J. Acoust. Soc. Am., 91, pp. 2663–2673. [CrossRef]
Dichmann, D. J., Maddocks, J. H., and Pego, R. L., 1996, “Hamiltonian Dynamics of an Elastica and the Stability of Solitary Waves,” Arch. Rational Mech. Anal., 135, pp. 357–396. [CrossRef]
Ito, K., 1998, “Uniform Stabilization of the Dynamic Elastica by Boundary Feedback,” SIAM J. Control Optim., 37, pp. 319–329. [CrossRef]
Ito, K., 2009, “Stabilization of the Dynamic Elastica Only by Damping Torque,” Nonlinear Anal.: Real World Appl., 10(5), pp. 3122–3131. [CrossRef]
Maddocks, J. H., and Dichmann, D. J., 1994, “Conservation Laws in the Dynamics of Rods,” J. Elast., 34, pp. 83–96. [CrossRef]
Lin, Y., and Pisano, A. P., 1987, “General Dynamic Equations of Helical Springs with Static Solution and Experimental Verification,” Trans. ASME, 54, pp. 910–917. [CrossRef]
Pai, F. P., 2007, Highly Flexible Structures: Modeling, Computation, and Experimentation, AIAA, Reston, VA.
Pai, P. F., 2007, Highly Flexible Structures: Modeling, Computation and Experimentation, AIAA, Reston, Virginia. Available at: http://web.missouri.edu/~paip/HFSs
Antman, S. S., 2005, Nonlinear Problems in Elasticity, 2nd ed., Springer, New York.
Goss, V. G. A., 2003, “Snap Buckling, Writhing and Loop Formation in Twisted Rods,” Ph.D. thesis, University College London, London.
van der Heijden, G. H. M., 2012, Home page of Gert van der Heijden. Available at: http://www.ucl.ac.uk/~ucesgvd
Preston, S. C., 2011, “The Geometry of Whips,” http://arxiv.org/abs/1105.1754
Preston, S. C., 2011, “The Motion of Whips and Chains,” J. Differ. Equ., 251(3), pp. 504–550. [CrossRef]
Thess, A., Zikanov, O., and Nepomnyashchy, A., 1999, “Finite-Time Singularity in the Vortex Dynamics of a String,” Phys. Rev. E, 59, pp. 3637–3640. [CrossRef]
Bauchau, O. A., and Bottasso, C. L., 1999, “On the Design of Energy Preserving and Decaying Schemes, for Flexible, Nonlinear Multi-Body Systems,” Comput. Methods Appl. Mech. Eng., 169, pp. 61–79. [CrossRef]
Celledoni, E., Grimm, V., McLachlan, R. I., McLaren, D. I., ONeale, D., Owren, B., and Quispel, G. R. W., 2012, “Preserving Energy Resp. Dissipation in Numerical PDEs Using the ‘Average Vector Field’ Method,” J. Comput. Phys., 231, pp. 6770–6789. [CrossRef]
Futterer, T., Klar, A., and Wegener, R., 2012, “An Energy Conserving Numerical Scheme for the Dynamics of Hyperelastic Rods,” Int. J. Differ. Equ., 2012, Article No. 718308. [CrossRef]
Leimkuhler, B., and Reich, S., 2004, Simulating Hamiltonian Dynamics, Cambridge University Press, Cambridge, UK.
Pace, B., Diele, F., and Marangi, C., 2012, “Energy Preservation in Separable Hamiltonian Systems by Splitting Schemes,” AIP Conf. Proc., 1479, pp. 1204–1207. [CrossRef]
Sanz-Serna, J. M., and Calvo, M. P., 1994, Numerical Hamiltonian Problems, Chapman and Hall, London, UK.
Simo, J. C., Marsden, J. E., and Krishnaprasad, P. S., 1988, “The Hamiltonian Structure of Nonlinear Elasticity: The Material and Convective Representations of Solids, Rods, and Plates,” Arch. Ration. Mech. Anal., 104(2), pp. 125–183. [CrossRef]
Simo, J. C., and Wong, K. K., 1991, “Unconditionally Stable Algorithms for Rigid Body Dynamics that Exactly Preserve Energy and Momentum,” Int. J. Numer. Methods Eng., 31(1), pp. 19–52. [CrossRef]
Simo, J. C., and Tarnow, N., 1992, “The Discrete Energy-Momentum Method. Conserving Algorithms for Nonlinear Elastodynamics,” Z. Math. Phys. (ZAMP), 43(5), pp. 757–792. [CrossRef]
Simo, J. C., Tarnow, N., and Wong, K. K., 1992, “Exact Energy-Momentum Conserving Algorithms and Symplectic Schemes for Nonlinear Dynamics,” Comput. Methods Appl. Mech. Eng., 100(1), pp. 63–116. [CrossRef]
Simo, J. C., Tarnow, N., and Doblare, M., 1995, “Non-Linear Dynamics of Three-Dimensional Rods: Exact Energy and Momentum Conserving Algorithms,” Int. J. Numer. Methods Eng., 38(9), pp. 1431–1473. [CrossRef]
Asmar, N. H., 2004, Partial Differential Equations with Fourier Series and Boundary Value Problems, 2nd ed., Prentice Hall, Englewood Cliffs, NJ.
Dunford, N., and Schwartz, J. T., 1963, Linear Operators. Part II. Spectral Theory. Selfadjoint Operators in Hilbert Space, John Wiley & Sons, Inc., New York.
Burden, R. L., and Faires, J. D., 2005, Numerical Methods, 8th ed., Thomson Brooks/Cole Publishing Co, Pacific Grove, CA.
Hairer, E., and Wanner, G., 2002, Solving Ordinary Differential Equations II: Stiff and Differential-Algebraic Problems, Springer Series in Computational Mathematics, 2nd ed., Springer, New York.
Hall, C. A., and Meyer, W. W., 1976, “Bounds for Cubic Spline Interpolation,” J. Approx. Theory, 16, pp. 105–122. [CrossRef]


Grahic Jump Location
Fig. 1

A heavily damped, coiled, elastic rod unraveling

Grahic Jump Location
Fig. 2

An elastic rod rotating about a point

Grahic Jump Location
Fig. 3

An elastic rod fixed at two ends, one of which is oscillating

Grahic Jump Location
Fig. 4

An elastic rod, clamped at one end, free at the other end, under the influence of gravity

Grahic Jump Location
Fig. 5

A steady state solution

Grahic Jump Location
Fig. 6

Plot of 4rErrr,r+1 for r = 0,1,2

Grahic Jump Location
Fig. 7

Plots of Errr,r+1/Errr+1,r+2 for r = 0,1




Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In