A story with twists and turns : how to control the rotation of the notched stick (2025)

Martin Luttmannmartin.luttmann@gmail.comUniversité Paris-Saclay, CEA, LIDYL, 91191 Gif-sur-Yvette, FranceDQML, IMX,École polytechnique fédérale de Lausanne (EPFL) Station 12,CH-1015 Lausanne, Switzerland  Michel Luttmannmichel.luttmann@cea.frCEA, CESTA, F-33114 Le Barp, France

(December 6, 2024)

Abstract

The notched stick, also known as the Gee-Haw-Whammy-Diddle, is a wooden toy able to convert linear vibration into rotational motion, whose behavior has been intriguing both children and physicists for decades. The oldest scientific article one can find on this subject was published 87 years ago in the present journal. Here we derive an analytical model of the system, supported by experimental results. We predict the direction of rotation, and explain why the device is so easy to operate, even without fine control of the various parameters. The potential importance of the vertical displacement of the finger exerting the perturbation force is also highlighted. We finally discuss similarities between the mechanical system described here and the optical effect of birefringence.

I Introduction

What physicist has not marveled, as a child, at a beautiful demonstration involving, for instance, soap bubbles, tippe tops, or Newton’s pendulum? The so-called notched stick is one of those scientific toys that stimulate the curiosity of children [1] and adults alike. It is an old wooden toy, made of a rod with a nail at the end. The nail serves as a shaft for a small propeller. A series of transverse notches is engraved along a portion of one edge of the stick. The player holds the notched stick in one hand, and a small cylindrical wooden dowel in the other. This second bar is used to rub the stick along its notched edge (Fig. 1), causing it to vibrate and emit a sound close to that of a rattle. Experienced players are able, using a method described below, to convert this vibration into a rotation of the propeller. As we shall see, the direction of rotation can be controlled at will by the player. The origin of this toy is sometimes attributed to Native Americans, while it is known today by different names around the world: Gee-Haw-Whammy-Diddle [2], Hooey stick [1] or Ouija Windmill in North America, Hui-Maschine [3] in Germany, Girigiri-Garigari [4] in Japan, etc. In France, the player pronounces the word Bozo-Bozo [5, 6] as a magic spell before the sudden reversal of the rotation, for the greatest delight of his/her audience.

A story with twists and turns : how to control the rotation of the notched stick (1)

In the following, we will refer to the bar with notches as the notched stick, and to the additional stick used to rub the notches as the excitation bar.To control the propeller’s direction of rotation, the player uses the following trick. As they rubs the notched stick with the excitation bar, they brings one of the fingers of the hand holding the latter into contact with the notched stick [7, 2, 1]. For a right-handed player, if the thumb is gently pushed against the top right-hand face of the stick, the propeller is driven in a counter-clockwise motion from their perspective (Fig. 2.a). Conversely, if the tip of the index finger is pressed on the top-left face of the notched stick, the rotation is clockwise (Fig. 2.b). By discreetly pressing on one side and then the other, the player thus reverses the direction of rotation (Supplementary video 1). If the player does not press on either side, the rotation is random, or does not take place at all.

A story with twists and turns : how to control the rotation of the notched stick (2)

A notched stick can be built in minutes, without worrying about the exact dimensions of the stick or the spacing and depth of the notches. With a bit of practice, the device is then surprisingly easy to operate. However, its working principle has been intriguing physicists for almost a century [8], and a dozen of scientific studies have been reported so far. Most of the papers published on the subject aim at describing the coupling between the motions of the nail and that of the propeller, assuming linear or elliptical [9] vibration of the nail. Recently, Marek et al. [10] and Broseghini et al. [11] demonstrated experimentally that a purely linear oscillation of the nail can lead to a rotation of the propeller, amplified by ordinary and/or parametric resonance (a phenomenon resembling the hula-hoop effect [12, 13], or Kapitza pendulums [14]). This is a case of spontaneous symmetry breaking, caused by the mismatch between the position of the nail and the propeller’s center of mass. Although there seem to be no definite consensus on the subject, these studies offer two major lessons. First, in almost all cases, a rotation of the pivot point (the nail) leads to a rotation of the propeller in the same direction, just as we do intuitively when twirling a hula-hoop around the hips, or a ring around the arm. Second, a linear oscillation of the pivot can also trigger a rotation of the propeller, although the direction is then seemingly random and determined by fluctuations or by the exact initial conditions [14, 11]. Hence, the mechanism of control of the propeller’s rotational direction is to be found in the origin of the motion of the nail itself.

A small number of papers have discussed the mechanism of control of the trajectory of the tip of the notched stick [8, 15, 16, 17, 7, 2, 3, 4], often not including a quantitative analysis. The most comprehensive studies of this aspect have been performed by Leonard [8] and Satonobu et al. [4]. In his theoretical work, Leonard considered a different configuration of the device, where the notched stick’s cross section is rectangular and the excitation bar positioned diagonally, and attributed the working mechanism to a difference of resonance frequency between the two cartesian components of the motion. This picture has several limitations. First, later studies [4] showed that the bending of the notched stick is negligible in this problem. Second, we observe that notched sticks with square or circular cross sections are also operational. Finally, the author does not precisely predict the direction of rotation or explains the mechanism for its reversal.Although they could not make quantitative predictions, Satonobu et al. performed valuable experiments and proposed instead that a difference between the waveforms of the force applied in the direction of the perturbation finger and the direction perpendicular to it could explain the behavior of the toy. We will discuss this idea in Section IV of this article. This intuition, however, lacks the support of a rigorous model that provides testable predictions.

In this paper, we clarify and unify previous approaches in a simple, more general mechanical model in order to investigate the origin of the rotation of the nail supporting the propeller. The end of the notched stick is modelled as a point mass moving in a 2D plane under the influence of various forces. Crucially, the point mass dynamics is influenced by two restoring forces corresponding to the action of the hand holding the notched stick and that of the finger pushing on it with adjustable direction [18]. Contrary to previous approaches, our model successfully predicts the direction of rotation, and highlights the role of the energy loss, which, in particular, determines the robustness of the mechanism. The model also accounts for the fact that the only parameter determining the direction of rotation is the pushing angle of the perturbation finger. New behaviors are also predicted, such as the dependence of the trajectory on the driving frequency. Next, we describe an experiment with an idealized version of the notched stick, driving its motion with a sinusoidal force and show that the experimental results are consistent with our model. Finally, we propose that the notched stick undergoes an additional effect, linked to the vertical displacement of the perturbation finger, which may lead to a more pronounced rotation. In the following, we will refer to the device shown in Fig. 1 as the real-life notched stick, to distinguish it from our laboratory version of the device. We also refer to the finger controlling the rotation as the perturbation finger.

II Analytical model

As previous studies demonstrated that the bending of the notched stick is negligible [4], we consider a perfectly rigid notched stick, initially parallel to the z𝑧zitalic_z axis. We assume that one end of the stick is fixed in space, while the other is free to move. Additionally, since the amplitude of the motion of the real-life notched stick’s tip is tiny compared to its length, we assume that the notched stick angle with respect to the z𝑧zitalic_z axis remains small. Therefore, the tip of the stick is approximately confined in the (x,y)𝑥𝑦(x,y)( italic_x , italic_y ) plane (Fig. 3.a), and its position vector in this plane is denoted r.

In this model, the oscillation of the notched stick’s tip is fully analogous to that of a point mass moving in a 2D plane under the influence of four distinct forces (Fig. 3.b). A driving force Fdriving(t)subscriptFdriving𝑡\textbf{F}_{\text{driving}}(t)F start_POSTSUBSCRIPT driving end_POSTSUBSCRIPT ( italic_t ) is directed along the y𝑦yitalic_y axis, and corresponds to the periodic force exerted by the excitation bar as it rubs the notches. A restoring force Frest(r)subscriptFrestr\textbf{F}_{\text{rest}}(\textbf{r})F start_POSTSUBSCRIPT rest end_POSTSUBSCRIPT ( r ) pulls the point mass back to the origin, mimicking the effect of the flexible hand holding the real-life notched stick. We assume Frest(r)=krsubscriptFrestr𝑘r\textbf{F}_{\text{rest}}(\textbf{r})=-k\textbf{r}F start_POSTSUBSCRIPT rest end_POSTSUBSCRIPT ( r ) = - italic_k r, where k𝑘kitalic_k is the stiffness constant of the holding hand. In order to take into account the effect of the perturbation finger, an additional force FpertsubscriptFpert\textbf{F}_{\text{pert}}F start_POSTSUBSCRIPT pert end_POSTSUBSCRIPT is applied on the stick at an angle φ𝜑\varphiitalic_φ away from the y𝑦yitalic_y axis (for instance Fig. 2.a corresponds to φ=π/4𝜑𝜋4\varphi=-\pi/4italic_φ = - italic_π / 4). To take into account the flexibility of the perturbation finger, FpertsubscriptFpert\textbf{F}_{\text{pert}}F start_POSTSUBSCRIPT pert end_POSTSUBSCRIPT is expressed as an elastic force acting along the uφsubscriptu𝜑\textbf{u}_{\varphi}u start_POSTSUBSCRIPT italic_φ end_POSTSUBSCRIPT vector (see Fig. 3.b)

Fpert(r)=q(ruφ)uφ,subscriptFpertr𝑞rsubscriptu𝜑subscriptu𝜑\textbf{F}_{\text{pert}}(\textbf{r})=-q(\textbf{r}\cdot\textbf{u}_{\varphi})%\textbf{u}_{\varphi},F start_POSTSUBSCRIPT pert end_POSTSUBSCRIPT ( r ) = - italic_q ( r ⋅ u start_POSTSUBSCRIPT italic_φ end_POSTSUBSCRIPT ) u start_POSTSUBSCRIPT italic_φ end_POSTSUBSCRIPT ,(1)

where q𝑞qitalic_q is the stiffness constant of the perturbation finger. This additional force breaks the mirror symmetry of the system, thus allowing for clockwise or counter-clockwise oscillations of the stick. From an energy perspective, the combined forces FrestsubscriptFrest\textbf{F}_{\text{rest}}F start_POSTSUBSCRIPT rest end_POSTSUBSCRIPT and FpertsubscriptFpert\textbf{F}_{\text{pert}}F start_POSTSUBSCRIPT pert end_POSTSUBSCRIPT correspond to an anisotropic energy potential, as mentioned in Ref. [17] which compares the notched stick to a ball rolling in a distorted watch glass. Finally, a damping force, that we write for simplicity Fdamp=αdrdtsubscriptFdamp𝛼drd𝑡\textbf{F}_{\text{damp}}=-\alpha\frac{\text{d}\textbf{r}}{\text{d}t}F start_POSTSUBSCRIPT damp end_POSTSUBSCRIPT = - italic_α divide start_ARG roman_d bold_r end_ARG start_ARG d italic_t end_ARG, accounts for the energy loss.

A story with twists and turns : how to control the rotation of the notched stick (3)

The point mass then obeys Newton’s second law

md2rdt2=Fdriving+Frest+Fpert+Fdamp,𝑚superscriptd2rdsuperscript𝑡2subscriptFdrivingsubscriptFrestsubscriptFpertsubscriptFdampm\frac{\text{d}^{2}\textbf{r}}{\text{d}t^{2}}=\textbf{F}_{\text{driving}}+%\textbf{F}_{\text{rest}}+\textbf{F}_{\text{pert}}+\textbf{F}_{\text{damp}},italic_m divide start_ARG d start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT r end_ARG start_ARG d italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG = F start_POSTSUBSCRIPT driving end_POSTSUBSCRIPT + F start_POSTSUBSCRIPT rest end_POSTSUBSCRIPT + F start_POSTSUBSCRIPT pert end_POSTSUBSCRIPT + F start_POSTSUBSCRIPT damp end_POSTSUBSCRIPT ,(2)

where m𝑚mitalic_m is the equivalent mass associated to the notched stick’s. Let us now consider the simplest possible driving force, a vertical sinusoidal force of frequency ω𝜔\omegaitalic_ω: Fdriving(t)=F0eiωtuysubscriptFdriving𝑡subscript𝐹0superscript𝑒𝑖𝜔𝑡subscriptu𝑦\textbf{F}_{\text{driving}}(t)=F_{0}e^{i\omega t}\,\textbf{u}_{y}F start_POSTSUBSCRIPT driving end_POSTSUBSCRIPT ( italic_t ) = italic_F start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT italic_i italic_ω italic_t end_POSTSUPERSCRIPT u start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT, where uysubscriptu𝑦\textbf{u}_{y}u start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT is the unit vector along y𝑦yitalic_y. Note that the back and forth motion of the excitation bar actually creates a variable torque on the notched stick, resulting in an additional low frequency component of the driving force. This additional component is neglected here for simplicity. To compute the steady-state solution, we express the cartesian coordinates of the point mass as x(t)=x0eiωt𝑥𝑡subscript𝑥0superscript𝑒𝑖𝜔𝑡x(t)=x_{0}\,e^{i\omega t}italic_x ( italic_t ) = italic_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT italic_i italic_ω italic_t end_POSTSUPERSCRIPT and y(t)=y0eiωt𝑦𝑡subscript𝑦0superscript𝑒𝑖𝜔𝑡y(t)=y_{0}\,e^{i\omega t}italic_y ( italic_t ) = italic_y start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT italic_i italic_ω italic_t end_POSTSUPERSCRIPT, where x0subscript𝑥0x_{0}italic_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT and y0subscript𝑦0y_{0}italic_y start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT are complex numbers. Plugging these and the expression of the four forces in Eq. 2, we obtain

mω2x(t)=kx(t)+q(x(t)sin2φy(t)cosφsinφ)+iωαx(t)mω2y(t)=ky(t)+q(y(t)cos2φx(t)cosφsinφ)+iωαy(t)+F0eiωt.𝑚superscript𝜔2𝑥𝑡𝑘𝑥𝑡𝑞𝑥𝑡superscript2𝜑𝑦𝑡𝜑𝜑𝑖𝜔𝛼𝑥𝑡𝑚superscript𝜔2𝑦𝑡𝑘𝑦𝑡𝑞𝑦𝑡superscript2𝜑𝑥𝑡𝜑𝜑𝑖𝜔𝛼𝑦𝑡subscript𝐹0superscript𝑒𝑖𝜔𝑡\begin{split}m\omega^{2}x(t)&=\\kx(t)&+q\big{(}x(t)\sin^{2}\varphi-y(t)\cos\varphi\sin\varphi\big{)}+i\omega%\alpha x(t)\\m\omega^{2}y(t)&=\\ky(t)&+q\big{(}y(t)\cos^{2}\varphi-x(t)\cos\varphi\sin\varphi\big{)}+i\omega%\alpha y(t)\\&+F_{0}e^{i\omega t}.\end{split}start_ROW start_CELL italic_m italic_ω start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_x ( italic_t ) end_CELL start_CELL = end_CELL end_ROW start_ROW start_CELL italic_k italic_x ( italic_t ) end_CELL start_CELL + italic_q ( italic_x ( italic_t ) roman_sin start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_φ - italic_y ( italic_t ) roman_cos italic_φ roman_sin italic_φ ) + italic_i italic_ω italic_α italic_x ( italic_t ) end_CELL end_ROW start_ROW start_CELL italic_m italic_ω start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_y ( italic_t ) end_CELL start_CELL = end_CELL end_ROW start_ROW start_CELL italic_k italic_y ( italic_t ) end_CELL start_CELL + italic_q ( italic_y ( italic_t ) roman_cos start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_φ - italic_x ( italic_t ) roman_cos italic_φ roman_sin italic_φ ) + italic_i italic_ω italic_α italic_y ( italic_t ) end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL + italic_F start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT italic_i italic_ω italic_t end_POSTSUPERSCRIPT . end_CELL end_ROW(3)

Dropping the t𝑡titalic_t dependence, we rewrite the first of the two above equations as

x=ξy,𝑥𝜉𝑦x=\xi\cdot y,italic_x = italic_ξ ⋅ italic_y ,(4)

with

ξ=qcosφsinφm(ω02ω2)+iαω,𝜉𝑞𝜑𝜑𝑚superscriptsubscript𝜔02superscript𝜔2𝑖𝛼𝜔\xi=\frac{q\cos\varphi\sin\varphi}{m(\omega_{0}^{2}-\omega^{2})+i\alpha\omega},italic_ξ = divide start_ARG italic_q roman_cos italic_φ roman_sin italic_φ end_ARG start_ARG italic_m ( italic_ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_ω start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) + italic_i italic_α italic_ω end_ARG ,(5)

and where

ω0=(k+qsin2φ)/m.subscript𝜔0𝑘𝑞superscript2𝜑𝑚\omega_{0}=\sqrt{(k+q\sin^{2}\varphi)/m}.italic_ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = square-root start_ARG ( italic_k + italic_q roman_sin start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_φ ) / italic_m end_ARG .(6)

The modulus of ξ𝜉\xiitalic_ξ corresponds to the ratio of the amplitudes of the x𝑥xitalic_x and y𝑦yitalic_y components of the motion, while its complex argument provides us with the dephasing between them. Since the argument of ξ𝜉\xiitalic_ξ is not zero in general, the trajectory of the notched stick’s tip is generally an ellipse. Eq. 4 alone is sufficient to compute the shape of the trajectory, while the full solution of Eq. 2 is obtained by plugging Eq. 4 in the second row of Eq. 3 (Supplementary Materials S1).

First, we note in Eq. 5 that a π/2𝜋2\pi/2italic_π / 2 variation of the perturbation angle φ𝜑\varphiitalic_φ changes ξ𝜉\xiitalic_ξ to ξ𝜉-\xi- italic_ξ, which corresponds to a π𝜋\piitalic_π phase shift between the two cartesian components of the motion. It thus yields a reversal of the ellipse’s direction of travel. Hence, the angle φ𝜑\varphiitalic_φ determines the direction of rotation of the stick’s tip, as observed with the real-life notched stick. For instance, φ=π/4𝜑𝜋4\varphi=\pi/4italic_φ = italic_π / 4 corresponds to counter-clockwise rotation, in agreement with the practical observation.

Secondly, we note in Eq. 5 that values of φ𝜑\varphiitalic_φ that are integer multiples of π/2𝜋2\pi/2italic_π / 2 yield ξ=0𝜉0\xi=0italic_ξ = 0. In that case, we have x(t)=0𝑥𝑡0x(t)=0italic_x ( italic_t ) = 0, thus the motion is linear vertical. On the other hand, if φ𝜑\varphiitalic_φ is an odd multiple of π/4𝜋4\pi/4italic_π / 4, the modulus of ξ𝜉\xiitalic_ξ is maximized, allowing non-degenerate elliptical motion. Note that the square cross-section of the notched stick (Fig. 2b,c) results in a non-zero component of the perturbation force pointing in the direction φ=±π/4𝜑plus-or-minus𝜋4\varphi=\pm\pi/4italic_φ = ± italic_π / 4, even if the player’s finger pushes horizontally.

Equation. 5 also indicates that at the resonance frequency, ω=ω0𝜔subscript𝜔0\omega=\omega_{0}italic_ω = italic_ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, ξ𝜉\xiitalic_ξ is a purely imaginary number. As such, at the resonance frequency, the components x(t)𝑥𝑡x(t)italic_x ( italic_t ) and y(t)𝑦𝑡y(t)italic_y ( italic_t ) are ±π/2plus-or-minus𝜋2\pm\pi/2± italic_π / 2 out of phase, the sign of the dephasing (determining the rotation direction) being that of cosφsinφ𝜑𝜑\cos\varphi\sin\varphiroman_cos italic_φ roman_sin italic_φ. The direction of rotation thus switches every π/2𝜋2\pi/2italic_π / 2 variation of φ𝜑\varphiitalic_φ, as observed with the real-life device. If, in addition, the modulus of ξ𝜉\xiitalic_ξ is close to 1111, then the trajectory of the notched stick’s tip is close to a circle.

A story with twists and turns : how to control the rotation of the notched stick (4)
A story with twists and turns : how to control the rotation of the notched stick (5)

The aforementioned features can be observed in Fig. 4, showing |ξ|𝜉|\xi|| italic_ξ | and the dephasing arg(ξ)arg𝜉\text{arg}(\xi)arg ( italic_ξ ) versus the perturbation angle and the driving frequency. The range of φ𝜑\varphiitalic_φ is limited to [π/2,π/2]𝜋2𝜋2[-\pi/2,\pi/2][ - italic_π / 2 , italic_π / 2 ] for clarity. The parameters k𝑘kitalic_k, q𝑞qitalic_q, m𝑚mitalic_m and α𝛼\alphaitalic_α of this plot correspond to the experiment reported later in the text. As expected, |ξ|𝜉|\xi|| italic_ξ | peaks for φ=±π/4𝜑plus-or-minus𝜋4\varphi=\pm\pi/4italic_φ = ± italic_π / 4 and ω𝜔\omegaitalic_ω close to ω0subscript𝜔0\omega_{0}italic_ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, while the dephasing is equal to ±π/2plus-or-minus𝜋2\pm\pi/2± italic_π / 2 for the resonance frequency ω=ω0𝜔subscript𝜔0\omega=\omega_{0}italic_ω = italic_ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT (blue curve). Note however that |ξ|𝜉|\xi|| italic_ξ | remains significant, and the dephasing close to ±π/2plus-or-minus𝜋2\pm\pi/2± italic_π / 2, for a broad range of frequency centered on ω0subscript𝜔0\omega_{0}italic_ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, of width approximately given by α/m𝛼𝑚\alpha/mitalic_α / italic_m [19]. This means that driving frequencies sufficiently close to ω0subscript𝜔0\omega_{0}italic_ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT will still excite a rotation of the device. Similarly, |ξ|𝜉|\xi|| italic_ξ | and the dephasing vary rather slowly with φ𝜑\varphiitalic_φ. As we will discuss later in the text, this weak dependence of ξ𝜉\xiitalic_ξ on the variables ω𝜔\omegaitalic_ω and φ𝜑\varphiitalic_φ explains why the device continues to work even if the driving frequency and the perturbation force direction are not precisely adjusted.

III Experiment with an idealized version of the device

The simple model detailed above indicates that the shape of the trajectory of the notched stick’s tip should depend on the driving frequency and the perturbation angle φ𝜑\varphiitalic_φ, while its direction of travel should only depend on φ𝜑\varphiitalic_φ.

The amplitude of the oscillation of the notched stick’s tip is small compared to its length, typically less than a millimeter for a 15 cm long stick. As such, recording this motion requires special techniques. Some previous experiments have used an accelerometer mounted on the tip of the notched stick, or high-speed cameras [11, 10]. Here, we propose a different, simpler method, that provides a visualization of the trajectory in real time.

The experimental setup is sketched in Fig. 5.a. The idea is to magnify the oscillation using a continuous laser beam, reflected at an angle on a mirror mounted on the tip of the notched stick. The other end of the stick is screwed onto a fixed vertical panel. The flexibility of the screw allows small oscillations of the stick about its equilibrium position. The total mass of our stick is 15 grams, thus yielding an equivalent mass of m=5×103𝑚5superscript103m=5\times 10^{-3}italic_m = 5 × 10 start_POSTSUPERSCRIPT - 3 end_POSTSUPERSCRIPT kg (this stems from the moment of inertia of a thin stick of length l𝑙litalic_l and mass m𝑚mitalic_m, equal to ml2/3𝑚superscript𝑙23ml^{2}/3italic_m italic_l start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / 3). We use a stretched rubber band to exert the perturbation force corresponding to the action of the finger in the real-life device. The direction of this force, i.e. the angle φ𝜑\varphiitalic_φ, can be varied by attaching the rubber band in various positions on the edge of a circular piece centered on the stick. The stiffness constant of the rubber band is estimated to be q50𝑞50q\approx 50italic_q ≈ 50 N/m. A vertical sinusoidal driving force with adjustable frequency is exerted on the device by the mean of a magnetic shaker (model Ling V203). A screen is placed 2m away from the mirror in order to observe the trajectory of the laser spot. The stick being 20 cm long, the trajectory magnifying factor is about 10. Since the persistence of vision is typically longer than the period of the oscillation, one can directly see the magnified trajectory on the screen.

A story with twists and turns : how to control the rotation of the notched stick (6)

Fig. 5.b shows the trajectories photographed by a camera (0.2-s exposure time), for φ=π/4𝜑𝜋4\varphi=\pi/4italic_φ = italic_π / 4 and various driving frequencies ω𝜔\omegaitalic_ω. Close to the resonance frequency, located around 250 rad/s, the trajectory takes the shape of a near-perfect ellipse. The inclination angle and eccentricity of this ellipse vary strongly with the frequency. At lower frequencies (ω<200𝜔200\omega<200italic_ω < 200 rad/s) and higher frequencies (ω>300𝜔300\omega>300italic_ω > 300 rad/s), the trajectory is almost purely linear vertical (which corresponds to a low value of |ξ|𝜉|\xi|| italic_ξ |) with an amplitude that decreases the further we are from the resonance frequency. Fig. S 5 of the Supplementary Materials shows a frequency sweep performed in the same conditions for φ=π/4𝜑𝜋4\varphi=-\pi/4italic_φ = - italic_π / 4, which exhibits a similar behavior, except for a reversal of the ellipse orientation angle and direction of rotation. These observations are in full agreement with our analytical model. However, tangled, non-elliptical trajectories are observed at certain low frequencies (for instance at ω=157𝜔157\omega=157italic_ω = 157 rad/s in Fig. 5.b). Although we have not yet identified the origin of these imperfections, we believe they may be caused by deviations from a perfectly sinusoidal excitation force.

To further understand the fine variations of the trajectory, we extract from these images the inclination angle ΨΨ\Psiroman_Ψ and the ratio b/a𝑏𝑎b/aitalic_b / italic_a of the minor axis to the great axis of the ellipses [20]. These are shown in Fig. 6 for φ=±π/4𝜑plus-or-minus𝜋4\varphi=\pm\pi/4italic_φ = ± italic_π / 4, with respect to the driving frequency, and fitted with the values expected from Eq. 5 (Supplementary Materials S3, S4). Importantly, the resonance frequency ω0subscript𝜔0\omega_{0}italic_ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT and the damping constant α𝛼\alphaitalic_α are adjustable parameters of the fits. Ideally, the curve of Fig. 6.a should be the opposite of that of Fig. 6.c, while the curves of Fig. 6.b,d should be identical. We attribute the discrepancies to experimental errors in the setting of the perturbation angle φ𝜑\varphiitalic_φ, as well as anisotropies in the restoring force of the screw holding the notched stick to the vertical panel. Nevertheless, as the match between the experimental points and the fitted curves is satisfying, the shapes of the trajectories provide us with an estimation of the parameters α𝛼\alphaitalic_α and ω0subscript𝜔0\omega_{0}italic_ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT in our experimental conditions. Taking the mean value of the best parameters of each fit, we find ω0=246.7±7.3subscript𝜔0plus-or-minus246.77.3\omega_{0}=246.7\pm 7.3italic_ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 246.7 ± 7.3 rad/s and α=0.297±0.138𝛼plus-or-minus0.2970.138\alpha=0.297\pm 0.138italic_α = 0.297 ± 0.138 kg/s (we provide here the 2σ2𝜎2\sigma2 italic_σ uncertainty of the set of four values). These values were plugged in Eq. 5 in order to plot the graphs of Fig. 4. Note that a direct measurement of α𝛼\alphaitalic_α would have been difficult. In particular, the dynamic contact between the magnetic shaker and the stick has an influence on the energy loss of the system, hence α𝛼\alphaitalic_α cannot be rigorously estimated by, for instance, measuring the damping time of the stick in the absence of the magnetic shaker.

IV Discussion

The model derived here, as well as the experiment reported, considers an idealized version of the notched stick, which is quite far from the real-life device. In the case of the real-life version, the perturbation finger and the excitation bar move together along the length of the notched stick. In addition, the vertical driving force is not a sine function, but rather a pulsed excitation resulting from the periodic impacts of the excitation bar on the notches. Furthermore, the real-life notched stick is heavily damped, being held by a tight hand. For these reasons, we can wonder whether the conclusions derived above still apply to the real toy manipulated by hand.

A simple way to test the effect of the displacement of the perturbation finger is to rub the notched stick with the excitation bar, while asking another person to push laterally with his finger on a given point of the notched stick (Supplementary video 2). While it is still possible to control the direction of rotation of the propeller in this way [4], it is less easy than with the standard method, shown in Fig. 2.

A possible explanation of this fact goes as follows [4]. The perturbation finger, pressing on either side of the notched stick, is integral with the excitation bar, and therefore moves synchronously with it (Fig. 7.a). As the excitation bar oscillates vertically as it passes each notch, the perturbation finger also moves up and down periodically (Fig. 7.b). The force FpertsubscriptFpert\textbf{F}_{\text{pert}}F start_POSTSUBSCRIPT pert end_POSTSUBSCRIPT exerted by the perturbation finger must vary as a result of this vertical displacement. Note that this effect was not present in the laboratory device of Fig. 5, since the attachment point of the rubber band exerting the perturbation force was static. To quantify the importance of this additional effect, we propose to modify the perturbation force FpertsubscriptFpert\textbf{F}_{\text{pert}}F start_POSTSUBSCRIPT pert end_POSTSUBSCRIPT of Eq. 1 to

Fpert(r)=q((rheiωtuy)uφ)uφ.subscriptFpertr𝑞rsuperscript𝑒𝑖𝜔𝑡subscriptu𝑦subscriptu𝜑subscriptu𝜑\textbf{F}_{\text{pert}}(\textbf{r})=-q\left((\textbf{r}-he^{i\omega t}\textbf%{u}_{y})\cdot\textbf{u}_{\varphi}\right)\textbf{u}_{\varphi}.F start_POSTSUBSCRIPT pert end_POSTSUBSCRIPT ( r ) = - italic_q ( ( r - italic_h italic_e start_POSTSUPERSCRIPT italic_i italic_ω italic_t end_POSTSUPERSCRIPT u start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ) ⋅ u start_POSTSUBSCRIPT italic_φ end_POSTSUBSCRIPT ) u start_POSTSUBSCRIPT italic_φ end_POSTSUBSCRIPT .(7)

hhitalic_h is a complex number whose modulus equals half the depth of the notches (about 1 mm), and whose argument corresponds to the dephasing between the vertical displacement of the finger and the driving force Fdriving(t)subscriptFdriving𝑡\textbf{F}_{\text{driving}}(t)F start_POSTSUBSCRIPT driving end_POSTSUBSCRIPT ( italic_t ). Supplementary Materials S5 details the computation of x(t)𝑥𝑡x(t)italic_x ( italic_t ) and y(t)𝑦𝑡y(t)italic_y ( italic_t ) given this new form for the perturbation force.

In addition, we note that our analytical model allows to compute the steady-state trajectory for a single-frequency excitation. Thus, to refine the model, we compute a sum of these trajectories for different driving frequencies, corresponding to the Fourier series of a sawtooth function. Such temporal driving force is a better model of the periodic impacts of the excitation bar on the notches. In addition, we consider a more realistic situation, with a spring constant q𝑞qitalic_q equal to k𝑘kitalic_k, which stems from the fact that the flexibility of the perturbation finger’s is close to that of the holding hand. Moreover, we consider a higher value of the damping factor α=2.97𝛼2.97\alpha=2.97italic_α = 2.97 kg/s, equal to ten times that of our experiment. This value is greater than the critical value α=2mω0𝛼2𝑚subscript𝜔0\alpha=2m\omega_{0}italic_α = 2 italic_m italic_ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT below which the freely moving notched stick oscillates periodically upon relaxation [21] (which does not occur with the real-life device).

A story with twists and turns : how to control the rotation of the notched stick (7)

Fig. 7.c shows the trajectory of the notched stick’s tip in this more realistic configuration, for h=00h=0italic_h = 0 and for various values of the driving frequency. The vertical driving force is given by the sum of the first 25 terms of the sawtooth Fourier series. Importantly, the trajectories are still close to ellipses, and are all travelled in the same (clockwise) direction. We also notice that the shape of the trajectory does not show a strong dependence over the value of ω𝜔\omegaitalic_ω, which stems from the width of the resonance scaling as α/msimilar-toabsent𝛼𝑚\sim\alpha/m∼ italic_α / italic_m. For the sake of comparison, Fig. S3 of the Supplementary Materials displays the trajectories in the same conditions, for α=2.97×101𝛼2.97superscript101\alpha=2.97\times 10^{-1}italic_α = 2.97 × 10 start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT kg/s, showing way more diverse trajectories. We conclude that the robustness of the system is conditioned by a sufficiently high value of the damping term, the real-life notched stick being heavily damped.

Secondly, in Fig. 7.d, we demonstrate that a vertical oscillation of the perturbation finger (h00h\neq 0italic_h ≠ 0) can reinforce the rotation mechanism described in detail above (resulting in a more eccentric elliptical trajectory). These trajectories were obtained for |h|=22|h|=2| italic_h | = 2 mm and arg(h)=π/2arg𝜋2\text{arg}(h)=\pi/2arg ( italic_h ) = italic_π / 2. We argue that close to the resonance frequency, the vertical displacement of the excitation bar is in phase with y(t)𝑦𝑡y(t)italic_y ( italic_t ), but in phase quadrature with the vertical driving force Fdriving(t)subscriptFdriving𝑡\textbf{F}_{\text{driving}}(t)F start_POSTSUBSCRIPT driving end_POSTSUBSCRIPT ( italic_t ) [22], effectively yielding arg(h)=±π/2argplus-or-minus𝜋2\text{arg}(h)=\pm\pi/2arg ( italic_h ) = ± italic_π / 2. We therefore conclude that a vertical displacement of the perturbation finger in parallel with the excitation rod might explain why the real-life notched stick is so efficient compared to the laboratory-built device. This effect needs to be studied in greater detail in future experiments, in order to quantify its importance in the operation of the real-life toy.

V Conclusion

Despite the simplicity of its design, the underlying physics of the notched stick is surprisingly rich. Here we focused on the causes of the motion of the notched stick’s tip. Our main conclusions are the following. First, we demonstrated that the control of the direction of rotation is made possible by an anisotropic restoring force acting on the tip of the notched stick, caused by applying a perturbation force at an angle of 45°  with respect to the driving force. Remarkably, the direction of the rotation only depends on the angle of the perturbation force, and in particular is independent from the driving frequency, the temporal shape of the driving force, or the value of the damping factor. Secondly, we revealed a rich frequency behavior that allowed us to validate our model and measure its parameters with an experiment. We also found a strong dependence of the shape of the notched stick’s trajectory on the value of the damping factor. More precisely, a significant energy loss in the system is necessary in order for the device to work outside of laboratory conditions, where the driving frequency cannot be controlled precisely. Finally, we identified an additional effect caused by the vertical oscillation of the perturbation finger, which could play a significant role in the operation of the notched stick.

The behavior of the real-life notched stick is a complex problem, and the rotation of the stick’s tip probably results from several effects combined. Two mechanisms are identified in this paper. The second mechanism was only studied numerically in Sec. IV, and it would be of great interest to conduct dedicated experiments on this aspect. Other subtle effects would be worth tackling, such that the role of the shape of the cross-section of the notched stick. Indeed, a rectangular cross-section seems to induce a rotation without the need for a lateral perturbation force [8, 7].

Since it is able to convert linear vibration into a controllable rotational motion, the notched stick is a type of vibrot [11, 23, 24]. Applications of such systems have been identified in engineering, in the context of screw or bolts loosening or tightening [25], or in robotics, for instance to build ultrasound-driven motors [26], or for the study of collective behaviors [27]. Similar vibration-induced rotations have been observed in various other systems with broken symmetry [28, 29, 30, 31].

Finally, armed with our new understanding, we can draw an analogy between the mechanism of the notched stick and birefringent materials in optics. In classical electromagnetism, it is common to assume that the negative electric charges in the matter are elastically bound to atomic nuclei. Birefringent materials exhibit an anisotropic electrical susceptibility, which corresponds to an anisotropic restoring force acting on the charge [32]. Thus, the electric charge is analogous to our point mass model of the notched stick’s tip. In this analogy, the driving force corresponds to the Lorentz’s force exerted on the charge by an incident linearly polarized light field. Note that since the oscillating electric charge radiates an electromagnetic wave, it must also lose kinetic energy, i.e. there is a damping force. In this picture, the change of polarization state caused by the propagation of light in the medium originates in the elliptical motion of the electric charges, which is similar to the elliptical motion’s of the notched stick’s tip. With this analogy, we explain why the notched stick has the same symmetries as waveplates. Indeed, rotating a quarter-waveplate by π/2𝜋2\pi/2italic_π / 2 in front of an incident linearly polarized beam of light reverses the helicity of the output field.

Acknowledgements.

The authors thank François Nadal, Laurent Potier and Robert Bédoret for helpful discussions.

AUTHOR DECLARATIONS

The authors have no conflicts to disclose.

References

A story with twists and turns : how to control the rotation of the notched stick (2025)

References

Top Articles
Latest Posts
Recommended Articles
Article information

Author: Errol Quitzon

Last Updated:

Views: 6208

Rating: 4.9 / 5 (59 voted)

Reviews: 82% of readers found this page helpful

Author information

Name: Errol Quitzon

Birthday: 1993-04-02

Address: 70604 Haley Lane, Port Weldonside, TN 99233-0942

Phone: +9665282866296

Job: Product Retail Agent

Hobby: Computer programming, Horseback riding, Hooping, Dance, Ice skating, Backpacking, Rafting

Introduction: My name is Errol Quitzon, I am a fair, cute, fancy, clean, attractive, sparkling, kind person who loves writing and wants to share my knowledge and understanding with you.