Abstract
The lift force acting on a purely translational galloping oscillator can be well approximated by using the quasi-steady theory, which states that the flow around a galloping body in motion is very similar to the known flow around a fixed body provided a minimum free stream velocity threshold and a similarity principle are satisfied. However, for oscillators undergoing coupled translational–rotational galloping, the rotation of the bluff body breaks the similarity principle, and therewith the quasi-steady assumption. Traditionally, the effects of rotation are accounted for by including explicit time-dependent terms in the lift force. However, we argue that the time-explicit representation of the lift force is unnecessary because the phenomenon of galloping is not time-explicit and acts only as a function of the body motion. Thus, we propose a modified time-implicit lift force representation, which has the same form of the quasi-steady theory, but with an additional dependence on the free stream velocity. The modified lift force is obtained by studying the transient growth of the amplitude of the bluff body oscillation. The proposed approach is used to model the lift force for square and trapezoidal prisms undergoing coupled translational–rotational galloping showing excellent prediction capabilities.
Similar content being viewed by others
Notes
The reader should bear in mind that perfectly symmetric bodies like a square are prone to galloping. This is because although the production of lift requires asymmetry, the asymmetry does not need to be in the geometry of the bluff body, but can also occur due to the angle of incidence between the moving body and the flow.
The reader can refer to the book by Païdoussis for a detailed review on the galloping instability [6].
References
Den Hartog, J.P.: Mechanical Vibrations. Courier Corporation, North Chelmsford (1985)
Parkinson, G.: Phenomena and modelling of flow-induced vibrations of bluff bodies. Prog. Aerosp. Sci. 26(2), 169–224 (1989)
Parkinson, G.V., Smith, J.D.: The square prism as an aeroelastic non-linear oscillator. Q. J. Mech. Appl. Math. 17(2), 225–239 (1964)
Nakamura, Y., Tomonari, Y.: Galloping of rectangular prisms in a smooth and in a turbulent flow. J. Sound Vib. 52(2), 233–241 (1977)
Nakamura, Y., Hirata, K., Urabe, T.: Galloping of rectangular cylinders in the presence of a splitter plate. J. Fluids Struct. 5(5), 521–549 (1991)
Païdoussis, M.P., Price, S.J., De Langre, E.: Fluid–Structure Interactions: Cross-Flow-Induced Instabilities. Cambridge University Press, Cambridge (2010)
Fung, Y.C.: An Introduction to the Theory of Elasticity. Wiley, New York (1955)
Blevins, R.D.: Flow-Induced Vibration. Van Nostrand Reinhold Co., New York (1977)
Blevins, R.D.: Flow-Induced Vibration, 2nd edn. Van Nostrand Reinhold Co., New York (1990)
Bearman, P.W., Gartshore, I.S., Maull, D.J., Parkinson, G.V.: Experiments on flow-induced vibration of a square-section cylinder. J. Fluids Struct. 1(1), 19–34 (1987)
Daqaq, M.F., Bibo, A., Akhtar, I., Alhadidi, A.H., Panyam, M., Caldwell, B., Noel, J.: Micropower generation using cross-flow instabilities: a review of the literature and its implications. J. Vib. Acoust. 141(3), 030801 (2019)
Ryabinin, A.N., Lyusin, V.D.: Galloping of small aspect ratio square cylinder. ARPN J. Eng. Appl. Sci. 10(1), 134–138 (2015)
Ewere, F., Wang, G.: Performance of galloping piezoelectric energy harvesters. J. Intell. Mater. Syst. Struct. 25(14), 1693–1704 (2013)
Nakamura, Y.: On the aerodynamic mechanism of torsional flutter of bluff structures. J. Sound Vib. 67(2), 163–177 (1979)
Van Oudheusden, B.W.: Rotational one-degree-of-freedom galloping in the presence of viscous and frictional damping. J. Fluids Struct. 10(7), 673–689 (1996)
Van Oudheusden, B.W.: Aerodynamic stiffness and damping effects in the rotational galloping of a rectangular cross-section. J. Fluids Struct. 14(8), 1119–1144 (2000)
Van Oudheusden, B.W.: Aerodynamic stiffness effects in rotational galloping at high wind speeds. J. Wind Eng. Ind. Aerodyn. 64(1), 31–46 (1996)
Kluger, J.M., Moon, F.C., Rand, R.H.: Shape optimization of a blunt body vibro-wind galloping oscillator. J. Fluids Struct. 40, 185–200 (2013)
Scanlan, R.H.: Problematics in formulation of wind-force models for bridge decks. J. Eng. Mech. 119(7), 1353–1375 (1993)
Dickinson, M.: The effects of wing rotation on unsteady aerodynamic performance at low Reynolds numbers. J. Exp. Biol. 192(1), 179–206 (1994)
Dickinson, M.H., Lehmann, F.-O., Sane, S.P.: Wing rotation and the aerodynamic basis of insect flight. Science 284(5422), 1954–1960 (1999)
Sane, S.P., Dickinson, M.H.: The aerodynamic effects of wing rotation and a revised quasi-steady model of flapping flight. J. Exp. Biol. 205(8), 1087–1096 (2002)
Sane, S.P.: The aerodynamics of insect flight. J. Exp. Biol. 206(23), 4191–4208 (2003)
Nayfeh, A.H., Mook, D.T.: Nonlinear Oscillations. Wiley, New York (2008)
Nayfeh, A.H.: Introduction to Perturbation Techniques. Wiley, New York (2011)
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Conflict of interest
The authors declare that they have no conflict of interest.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Appendices
Appendix A: Equation of motion
Following linear Euler–Bernoulli beam theory for which the thickness, \(t_b\), and deflection w(x, t) of the beam are assumed to be much smaller than its length L, we can write the equation governing the motion of the beam as
where t represents the temporal variable, and x represents the spatial variable along the length of the beam, the dot represents a derivative with respect to t and the prime represents a derivative with respect to x, \(M_b\) is the mass of the beam, \(c_v\) is the linear damping coefficient, \(E_b\) and \(I_b\) are, respectively, the modulus of elasticity, and mass moment of inertia of the beam, \(F_N\) is the vertical component of the aerodynamic force, and \(\delta \) is the Dirac delta function. Note that the aerodynamic force is assumed to be concentrated at the tip, which implicitly assumes that the width of the bluff body, D, is much smaller than the length of the beam, L.
The proper boundary conditions pertinent to Fig. 2 are given by
where \(M_t\) and \(I_t\) are, respectively, the mass and mass moment of inertia of the bluff body.
As shown in Fig. 2, the aerodynamic force acting on the moving tip body is produced by a relative wind velocity \(V_{\mathrm{rel}}\) at an effective angle of attack \(\alpha \), which can be written as
where \(\alpha _t=\tan ^{-1}\left[ \frac{{\dot{w}}(L,t)}{V}\right] \). The corresponding drag and lift force components are given by
where \(\rho _a\) is the air density, \(C_{F_D}\) and \(C_{F_L}\) are, respectively, the drag and lift coefficients. The normal force acting in the y-direction can be, therefore, written as
When \({\dot{w}}(L,t)/V \ll 1\), the previous equation becomes
Furthermore, under the assumption that \({\dot{w}}(L,t)/V \ll 1\), \(V_{\mathrm{rel}}\) can be approximated by V, which yields
Assuming that the beam is free of internal resonances and that the fluid forces are incapable of directly exciting the higher vibration modes, a single-mode representation of the beam dynamics can be written as
where \(\phi (x)\) is a spatial representation of the beam shape at a given time, chosen here to be the first mode shape of a cantilever beam with a tip mass. The function \({\bar{y}}(t)\) is a temporal function representing the time variation of the mode shape. The function \(\phi (x)\) can be written as
where
Substituting Eq. (A.8) into Eq. (A.1), multiplying by \(\phi (x)\) and integrating over the length of the beam yields
where
Normalizing the mode shape, \(\phi (x)\), by choosing A such that \(\phi (L)=1\), yields \(f_0=1\) and \(\alpha = \dot{{\bar{y}}}/V+ \phi ^{\prime }(L) {\bar{y}}\). Furthermore, for small deflections \({\dot{w}}(L,t)/L < 0.3\), the function \(\phi ^{\prime }(L)\) can be further approximated by \(\phi ^{\prime }(L) \approx -\frac{3}{2 L} \phi (L)\) . Thus, Eq. (A.10) can be further simplified as
which is the same as Eq. (2) used to study the dynamics of the galloping oscillator.
Appendix B: Approximate analytical solution of the equation of motion
An approximate solution of Eq. (6) can be obtained using the method of multiple scales; a method widely utilized to analyze the response of weakly nonlinear systems containing mechanisms that act on distinct time scales [25]. In the case of a galloping body, a “true” time scale describes the instantaneous position of the body, while a slow time scale tracks the weakly nonlinear evolution of oscillation amplitude.
The scaling is accomplished using an arbitrarily small bookkeeping parameter \(\epsilon \) such that “true” time is \(T_0 = \tau \) and “slow” time is \(T_1 = \epsilon \tau \). For the present problem, higher orders of \(\epsilon \) (even slower time scales) are neglected. The derivative operator is adapted such that
where
Additionally, the response \(y(\tau )\) is described into multiple scales such that
It can be reasoned that the mechanism of amplitude growth occurs slowly because the energy input to the system is on the order of the energy dissipated in the context of underdamped vibrations, energy dissipation must be small (controlled by \(\zeta _m\)). As such, a small net energy input requires small input flow energy (controlled by \(\mu \)). Therefore, the appropriate parameters are scaled by \(\epsilon \) as
When Eqs. (B.1), (B.3), and (B.4) are substituted into Eq. (6), the resulting system can be divided into subsystems by collecting the zeroth and first power of \(\epsilon \) as following:
For \(\epsilon ^0\):
admits a zeroth-order solution
where a is the amplitude of oscillation and \(\gamma \) is a phase shift function.
For \(\epsilon ^1\):
Letting \(\phi = T_0 + \gamma \) and substituting the zeroth-order solutions, as well as the expression for \(C_N\) from Eq. (3) into Eq. (B.7) yields
Because the physical system does not have oscillations which grow to infinity, the “forcing” terms on the right-hand side of Eq. (B.8) must not resonate with \(y_1\). This condition is enforced by ensuring the right-hand side of Eq. (B.8) is orthogonal to its homogeneous solution throughout one cycle of oscillation. That is, the right-hand side of the equation should be orthogonal to both \(\sin \phi \) and \(\cos \phi \).
For orthogonality with \(\cos \phi \), we obtain
Noting that \(\alpha =\frac{y^{\prime }}{U}-\xi U\), and substituting Eq. (B.6) into the expression for \(\alpha \) we obtain
which implies that the angle of attack changes sign when
Thus, Eq. (B.9) becomes
where \(\theta _0=\tan ^{-1}\left( \xi U\right) \). Upon carrying the integration in Eq. (B.12), we obtain
where the \(C_j\) are given by
Following a similar procedure to enforce the orthogonality with \(\sin \phi \) yields
At steady state, \(\frac{\text {d} \gamma }{\text {d} \tau } = \frac{\text {d} a}{\text {d} \tau } = 0\). Root-solving algorithms can then be used to solve Eqs. (B.13) and (B.15) for the steady-state phase shift and amplitude of the response, respectively.
Rights and permissions
About this article
Cite this article
Noel, J., Alhadidi, A.H., Alhussien, H. et al. A time-implicit representation of the lift force for coupled translational–rotational galloping. Nonlinear Dyn 103, 2183–2196 (2021). https://doi.org/10.1007/s11071-021-06232-6
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s11071-021-06232-6