Research Papers

Nonlinear Oscillations Induced by Follower Forces in Prestressed Clamped Rods Subjected to Drag

[+] Author and Article Information
Soheil Fatehiboroujeni

Department of Mechanical Engineering,
University of California,
Merced, CA 95343
e-mail: sfatehiboroujeni@ucmerced.edu

Arvind Gopinath

Department of Bioengineering,
Health Science Research Institute,
University of California,
Merced, CA 95343
e-mail: agopinath@ucmerced.edu

Sachin Goyal

Department of Mechanical Engineering,
Health Science Research Institute,
University of California,
Merced, CA 95343
e-mail: sgoyal2@ucmerced.edu

Contributed by the Design Engineering Division of ASME for publication in the JOURNAL OF COMPUTATIONAL AND NONLINEAR DYNAMICS. Manuscript received March 19, 2018; final manuscript received September 30, 2018; published online October 29, 2018. Assoc. Editor: Anindya Chatterjee.

J. Comput. Nonlinear Dynam 13(12), 121005 (Oct 29, 2018) (8 pages) Paper No: CND-18-1111; doi: 10.1115/1.4041681 History: Received March 19, 2018; Revised September 30, 2018

Elastic-driven slender filaments subjected to compressive follower forces provide a synthetic way to mimic the oscillatory beating of biological flagella and cilia. Here, we use a continuum model to study the dynamical, nonlinear buckling instabilities that arise due to the action of nonconservative follower forces on a prestressed slender rod clamped at both ends and allowed to move in a fluid. Stable oscillatory responses are observed as a result of the interplay between the structural elastic instability of the inextensible slender rod, geometric constraints that control the onset of instability, energy pumped into the system by the active follower forces, and motion-driven fluid dissipation. Initial buckling instabilities are initiated by the effect of the follower forces and inertia; fluid drag subsequently allows for the active energy pumped into the system to be dissipated away and results in self-limiting amplitudes. By integrating the equations of equilibrium and compatibility conditions with linear constitutive laws, we compute the critical follower forces for the onset of oscillations, emergent frequencies of these solutions, and the postcritical nonlinear rod shapes for two forms of the drag force, namely linear Stokes drag and quadratic Morrison drag. For a rod with fixed inertia and drag parameters, the minimum (critical) force required to initiate stable oscillations depends on the initial slack and weakly on the nature of the drag force. Emergent frequencies and the amplitudes postonset are determined by the extent of prestress as well as the nature of the fluid drag. Far from onset, for large follower forces, the frequency of the oscillations can be predicted by evoking a power balance between the energy input by the active forces and the dissipation due to fluid drag.

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


Langthjem, M. , and Sugiyama, Y. , 2000, “ Dynamic Stability of Columns Subjected to Follower Loads: A Survey,” J. Sound Vib., 238(5), pp. 809–851. [CrossRef]
Elishakoff, I. , 2005, “ Controversy Associated With the so-Called “Follower Forces”: Critical Overview,” ASME Appl. Mech. Rev., 58(2), p. 117. [CrossRef]
Bolotin, V. V. , 1999, “ Dynamic Instabilities in Mechanics of Structures,” ASME Appl. Mech. Rev., 55(1), pp. R1–R9. [CrossRef]
Leipholz, H. H. E. , 1980, Stability of Elastic Systems/Horst Leipholz, Sijthoff and Noordhoff/Alphen aan den Rijn, Alphen aan den Rijn, The Netherlands.
Reut, V. I. , 1939, “ About the Theory of Elastic Stability,” Odessa Institute of Civil and Communal Engineering, p. 1.
Pfluger, A. , 1950, Stabilitatsprobleme Der Elastostatik, Springer-Verlag, Berlin.
Beck, M. , 1952, “ Die Knicklast Des Einseitig Eingespannten, Tangential Gedruckten Stabes,” ZAMP Z. Fur Angew. Math. Phys., 3(3), pp. 225–228. [CrossRef]
Païdoussis, M. P. , 2016, Fluid-Structure Interactions: Slender Structures and Axial Flow, 2nd ed., Academic Press, London.
Païdoussis, M. P. , and Li, G. X. , 1993, “ Pipes Conveying Fluid: A Model Dynamical Problem,” J. Fluids Struct., 7(2), pp. 137–204. [CrossRef]
Wood, W. G. , Saw, S. S. , and Saunders, P. M. , 1969, “ The Kinetic Stability of a Tangentially Loaded Strut,” Proc. R. Soc. A: Math., 313(1513), pp. 239–248. [CrossRef]
Païdoussis, M. , 1973, “ Dynamics of Cylindrical Structures Subjected to Axial Flow,” J. Sound Vib., 29(3), pp. 365–385. [CrossRef]
Kang, B. , and Tan, C. A. , 2000, “ Parametric Instability of a Leipholz Column Under Periodic Excitation,” J. Sound Vib., 229(5), pp. 1097–1113. [CrossRef]
Dreyfus, R. , Baudry, J. , Roper, M. L. , Fermigier, M. , Stone, H. A. , and Bibette, J. , 2005, “ Microscopic Artificial Swimmers,” Nature, 437(7060), pp. 862–865. [CrossRef] [PubMed]
Babataheri, A. , Roper, M. , Fermigier, M. , and Du Roure, O. , 2011, “ Tethered Fleximags as Artificial Cilia,” J. Fluid Mech., 678, pp. 5–13. [CrossRef]
Sasaki, Y. , Takikawa, Y. , Jampani, V. S. R. , Hoshikawa, H. , Seto, T. , Bahr, C. , Herminghaus, S. , Hidaka, Y. , and Orihara, H. , 2014, “ Colloidal Caterpillars for Cargo Transportation,” Soft Matter, 10(44), pp. 8813–8820. [CrossRef] [PubMed]
Patteson, A. E. , Gopinath, A. , and Arratia, P. E. , 2016, “ Active Colloids in Complex Fluids,” Curr. Opin. Colloid Interface Sci., 21, pp. 86–96. [CrossRef]
Chelakkot, R. , Gopinath, A. , Mahadevan, L. , and Hagan, M. F. , 2014, “ Flagellar Dynamics of a Connected Chain of Active, Polar, Brownian Particles,” J. R. Soc., Interface, 11(92), p. 20130884. [CrossRef]
Vaziri, A. , and Gopinath, A. , 2008, “ Cell and Biomolecular Mechanics in Silico,” Nat. Mater., 7(1), pp. 15–23. [CrossRef] [PubMed]
Gopinath, A. , and Mahadevan, L. , 2011, “ Elastohydrodynamics of Wet Bristles, Carpets and Brushes,” Proc. R. Soc. A, 467(2130), pp. 1665–1685. [CrossRef]
Vaziri, A. , Gopinath, A. , and Deshpande, V. S. , 2007, “ Continuum-Based Computational Models for Cell and Nuclear Mechanics,” J. Mech. Mater. Struct., 2(6), pp. 1169–1191. [CrossRef]
Qin, B. , Gopinath, A. , Yang, J. , Gollub, J. P. , and Arratia, P. E. , 2015, “ Flagellar Kinematics and Swimming of Algal Cells in Viscoelastic Fluids,” Sci. Rep., 5, p. 9190.
Maghsoodi, A. , Chatterjee, A. , Andricioaei, I. , and Perkins, N. C. , 2017, “ Dynamic Model Exposes the Energetics and Dynamics of the Injection Machinery for Bacteriophage T4,” Biophys. J., 113(1), pp. 195–205. [CrossRef] [PubMed]
Maghsoodi, A. , Chatterjee, A. , Andricioaei, I. , and Perkins, N. C. , 2016, “ A First Model of the Dynamics of the Bacteriophage T4 Injection Machinery,” ASME J. Comput. Nonlinear Dyn., 11(4), p. 041026. [CrossRef]
Goyal, S. , 2006, “ A Dynamic Rod Model to Simulate Mechanics of Cables and DNA,” Ph.D. dissertation, University of Michigan, Ann Arbor, MI.
Goyal, S. , Perkins, N. C. , and Lee, C. L. , 2005, “ Nonlinear Dynamics and Loop Formation in Kirchhoff Rods With Implications to the Mechanics of DNA and Cables,” J. Comput. Phys., 209(1), pp. 371–389. [CrossRef]
De Canio, G. , Lauga, E. , and Goldstein, R. E. , 2017, “ Spontaneous Oscillations of Elastic Filaments Induced by Molecular Motors,” J. R. Soc. Interface, 14(136), p. 20170491.
Herrmann, G. , and Bungay, R. W. , 1964, “ On the Stability of Elastic Systems Subjected to Nonconservative Forces,” ASME J. Appl. Mech., 31(3), pp. 435–440. [CrossRef]
Bayly, P. V. , and Dutcher, S. K. , 2016, “ Steady Dynein Forces Induce Flutter Instability and Propagating Waves in Mathematical Models of Flagella,” J. R. Soc. Interface, 13(123), p. 20160523.
Kirchhoff, G. , 1859, “ Uber Das Gleichgewicht Und Die Bewegung Eines Unendlich Dunnen Elastischen Stabes,” J. Reine Angew. Math., 56(56), pp. 285–343. [CrossRef]
Chung, J. , and Hulbert, G. M. , 1993, “ A Time Integration Algorithm for Structural Dynamics With Improved Numerical Dissipation - the Generalized-Alpha Method,” ASME J. Appl. Mech., 60(2), pp. 371–375. [CrossRef]
Gobat, J. I. , and Grosenbaugh, M. A. , 2001, “ Application of the Generalized-Alpha Method to the Time Integration of the Cable Dynamics Equations,” Comput. Methods Appl. Mech. Eng., 190(37–38), pp. 4817–4829. [CrossRef]
van der Heijden, G. , Neukirch, S. , Goss, V. , and Thompson, J. , 2003, “ Instability and Self-Contact Phenomena in the Writhing of Clamped Rods,” Int. J. Mech. Sci., 45(1), pp. 161–196. [CrossRef]
Fatehiboroujeni, S. , Palanthandalam-Madapusi, H. J. , and Goyal, S. , 2018, “ Computational Rod Model With User-Defined Nonlinear Constitutive Laws,” ASME. J. Comput. Nonlinear Dynam., 13(10), p. 101006. [CrossRef]
Goyal, S. , Perkins, N. , and Lee, C. L. , 2008, “ Non-Linear Dynamic Intertwining of Rods With Self-Contact,” Int. J. Non-Linear Mech., 43(1), pp. 65–73. [CrossRef]
Anwar, Z. , Gopinath, A. , and Armstrong, R. C. , 2012, “ Systems Analysis of Hybrid, Multi-Scale Complex Flow Simulations Using Newton-GMRES,” Rheol. Acta, 51(9), pp. 849–866. [CrossRef]
Gear, C. W. , Kaper, T. J. , Kevrekidis, I. G. , and Zagaris, A. , 2005, “ Projecting to a Slow Manifold: Singularly Perturbed Systems and Legacy Codes,” SIAM J. Appl. Dyn. Syst., 4(3), pp. 711–732. [CrossRef]


Grahic Jump Location
Fig. 1

(a) Schematic representation of a rod of unstressed length L with fixed–fixed boundary condition (clamped at both ends). The end-to-end distance when buckled is Lee < L. (b) The motion of material points comprising the cross section of the rod at arc-length position, s and at time t is determined by tracking the transformations of the body-fixed frame âi(s,t) with respect to the inertial frame of reference êi. (c) The shape (top) and prestress (bottom) in the buckled state for different values of Lee/L. The dashed line corresponds to the unbuckled case Lee/L = 1.0. Prestress here is defined as the component of the internal force in the direction of cross-sectional normal vector â3(s,t) i.e., f3. Here, we use the tension along the filament, f3, to characterize prestress. For the shapes we study in this paper, there is a one-to-one correspondence between the two.

Grahic Jump Location
Fig. 2

(Left) A measure of the prestress may be gained by examining the plot of f3(L) versus 1 − Lee/L. The results of this paper focus on the range 0.05 < 1 − Lee/L < 0.5 shown between the two markers (red circles). (Right) Critical load for onset of oscillations Fcr versus scaled decrease in end-to-end distance 1 − Lee/L for both Stokes [S] drag and Morrison [M] drag. We note that the critical loads are roughly the same for the two forms of fluid drag and deviate little from the nondrag value. This is not surprising for the Morrison drag [M] as it is nonlinear and hence does not enter the linear stability equations at leading order for small perturbations. Given that the Stokes [S] values are close to the no drag values, we surmise that the onset of vibrations and onset of flutter are very close to one another for the parameter range investigated. For 0.05 < 1 − Lee/L < 0.5, the critical force Fcr increases as 1 − Lee/L, or prestress, increases. For the region 1 − Lee/L < 0.05, we are unable to obtain well-defined results since the slack is very small in this region and inextensibility constraint makes the problem very stiff necessitating very small time intervals. Since the focus of our paper is to investigate the role of slack (prestress), we choose to study values of Lee/L < 0.9.

Grahic Jump Location
Fig. 3

Configurations of the oscillating rod when |F|=15 N/m, when the drag force is of the Stokes [S] form. We show kymographs of the curvature in the top column, as well as the shapes over a period (labeled sequentially from 1–10) in the bottom column. We see increasingly sharper shapes (higher values of |κ2|) for smaller values of Lee/L.

Grahic Jump Location
Fig. 4

Configurations of the oscillating rod when |F|=15 N/m, when the drag force is of the Morrison [M] form. We show kymographs of the curvature in the top column, as well as the shapes over a period (labeled sequentially from 1–10) in the bottom column. As in the [S] case, we see increasingly sharper shapes (higher values of |κ2|) for smaller values of Lee/L. Furthermore, the shapes are smoother and the oscillations involve lower amplitudes than the corresponding shapes for the [S] case.

Grahic Jump Location
Fig. 5

Fourier transform (in the time domain) of the shear force at the midspan length of the rod shows that higher harmonics (insets show the raw data) are damped in the case of the nonlinear Morrison drag ((a) and (b)) more effectively than for the linear Stokes drag ((c) and (d)). The right column of the picture corresponds to Lee/L = 0.9 while the left column corresponds to Lee/L = 0.7. The ratio of drag coefficients for both cases is 3.18. Intuitively, we expect this ratio to affect the extent of dampening. The magnitude of |F| for all cases is 14 N/m.

Grahic Jump Location
Fig. 6

Frequency of oscillations for rods as a function of |F| for various values of Lee when the fluid drag is of the Stokes [S] form. The form of the curves suggest possible transitions that may be related to the activation or deactivation of higher order mode shapes (typical shapes are shown alongside the curves). Sudden jumps in the frequency–force curve, for example from point A to point B, are due to the suppression of oscillation modes with smaller wavelength that also have higher energy level. Higher follower force densities yield increasingly steep shapes with localized curvature variations. When plotted on logarithmic axes, we find that the frequencies do not correlate well with the ω∼|F|4/3 power law form for an inertialess cantilever with distributed follower forces, oscillating with a single dominant wavelength and frequency.

Grahic Jump Location
Fig. 7

(a) Frequency for the Morrison [M] drag plotted as a function of the force density |F| plotted in logarithmic scales to illustrate two salient features—(i) as the follower force increases to values much larger than the critical values, the effect of the prestress diminishes, and (ii) the frequencies in the limit |F|≫Fcr scale roughly as ω∼|F|5/6 consistent with our theoretical prediction. (b) Emergent frequency plotted as a function of the scaled end-to-end distance showing nonmonotonic behavior at fixed values of |F|. However, we note that as |F| increases, the effect of the prestress and slack becomes decreasingly important for the range of Lee/L investigated.

Grahic Jump Location
Fig. 8

Fourier transform (in the time domain) of the shear force at the midspan length of the rod subject to Stokes drag shows that higher harmonics are less pronounced as follower load, |F|, increases from |F|=Fcr=11.5 N/m (point A) to |F|=13 N/m (point B) for Lee/L = 0.85. This correlates with the increase of steady-state frequency of the oscillations from point A to pointB.



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