[2]\fnmZainul \surAbidin

1]\orgdivXI Natural Science, \orgnameDarma Yudha Senior High School, \orgaddress\street189 S.M. Amin Street, \cityPekanbaru, \postcode28292, \stateRiau, \countryIndonesia

[2]Simetri Foundation, Tangerang, Indonesia 15334

Playing Lato-lato is Difficult and This is Why

\fnmFansen Candra \surFunata fansencandra.funata@gmail.com โ€ƒโ€ƒ zainul.abidin@yayasansimetri.or.id [ *
Abstract

Lato-lato, a pendulum-based toy gaining popularity in Indonesian playgrounds, has sparked interest with competitions centered around maintaining its oscillatory motion. While some find it easy to play, the challenge lies in sustaining the oscillation, particularly in maintaining both โ€up and down collisions.โ€ Through a Newtonian dynamics numerical analysis using Python (code by ChatGPT), this study identifies two equilibrium phases - phase 1, characterized by normal pendulum motion, and phase 2, the double collision mode - by using the driven oscillation model. In addition, further analysis and discussion are done using the obtained numeric data. The difficulty in remaining in phase 2 highlights the intricate hand-eye coordination required, shedding light on the toyโ€™s appeal and the skill it demands.

keywords:
lato-lato, driven oscillation, pendulum-based toy, Newtonian dynamics, Python, numerical analysis, ChatGPT

1 Introduction

Lato-lato, or better known as โ€clacker ballsโ€ has existed since the 1960s, originally made of tempered glass. Due to safety issues, the balls are then changed to be made out of plastic [1]. It recently regained its popularity, especially in Indonesia, because it was played by the president of Indonesia [2]. Its appearance in social media, such as TikTok made it even more viral, with videos soaring to millions of views. Furthermore, playing this game positively impacts childrenโ€™s behaviour, increasing the frequency of their social interaction[3].
Previously some reports and articles have been written revolving around lato-lato. However, most of them seem to focus on how educators can implement lato-lato in teaching physics[4][5][6]. As we all know, the physical essence of the toy lato-lato lies in the law of momentum conservation; the collisions that occur between the plastic spheres are such that pi=pfsubscript๐‘๐‘–subscript๐‘๐‘“p_{i}=p_{f}italic_p start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT = italic_p start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT[7][8]. However, the dynamics of the lato-lato itself (2.3), involve gruesome mathematics, as shown by Bartucelli et al[9]. In this paper, we will use the concept of phase diagrams to explain why the lato-lato is such a difficult game. Moreover, A diagram of the Amplitude Conditions111Variable defined to show which phase the pendulum is in will be shown against the initial boundary conditions, ฮธ0subscript๐œƒ0\theta_{0}italic_ฮธ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT and ฮธห™0subscriptห™๐œƒ0\dot{\theta}_{0}overห™ start_ARG italic_ฮธ end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT. Further, we model the lato-lato as a 2 stick pendulum joint at both of its free ends.

Refer to caption
Figure 1: (a) shows the model of the lato-lato used throughout our analysis. We use the approximation r<<lmuch-less-than๐‘Ÿ๐‘™r<<litalic_r < < italic_l, hence the sphere may be considered a point mass. (b) Represents the motion of โ€phase 1โ€, and (c) represents the motion of โ€phase 2โ€. Figure (d) shows a special play style called the โ€tornado styleโ€ (more about it in Subsection 2.6)

2 Formula Derivations

The concept involved in playing with this toy lies mainly in Newtonโ€™s laws and conservation of momentum and energy[10]. Some of the subsections we provide here serve as preliminary materials to aid the readers in grasping the materials as a whole (Section 2.1 and Section 2.2).
The formulas used are as follows

2.1 Single Pendulum Equation of Motion

We begin with a single pendulum case with its energy given as

E=mโขgโขlโข(1โˆ’cosโกฮธ)+mโขl2โขฮธห™22๐ธ๐‘š๐‘”๐‘™1๐œƒ๐‘šsuperscript๐‘™2superscriptห™๐œƒ22E=mgl(1-\cos\theta)+ml^{2}\frac{\dot{\theta}^{2}}{2}italic_E = italic_m italic_g italic_l ( 1 - roman_cos italic_ฮธ ) + italic_m italic_l start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT divide start_ARG overห™ start_ARG italic_ฮธ end_ARG start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 2 end_ARG (1)

Where E๐ธEitalic_E represents the energy of the system, m๐‘šmitalic_m represents the mass of the bob, g๐‘”gitalic_g represents the gravitational acceleration (9.89.89.89.8 m/s2๐‘šsuperscript๐‘ 2m/s^{2}italic_m / italic_s start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT), l๐‘™litalic_l represents the length of the pendulum, ฮธ๐œƒ\thetaitalic_ฮธ represents the deviation angle with respect to the y axis and ฮธห™ห™๐œƒ\dot{\theta}overห™ start_ARG italic_ฮธ end_ARG represents the derivative with respect to time, as shown in Figure 1.

Without energy loss, we consider three separate cases, which are (a) E<<2โขmโขgโขlmuch-less-than๐ธ2๐‘š๐‘”๐‘™E<<2mglitalic_E < < 2 italic_m italic_g italic_l, (b) E=2โขmโขgโขl๐ธ2๐‘š๐‘”๐‘™E=2mglitalic_E = 2 italic_m italic_g italic_l and (c) E>2โขmโขgโขl๐ธ2๐‘š๐‘”๐‘™E>2mglitalic_E > 2 italic_m italic_g italic_l. Here, we ought to find the angular momentum L๐ฟLitalic_L of the system to obtain plots of the systemโ€™s phase diagram

  1. (a)

    E<<2โขmโขgโขlmuch-less-than๐ธ2๐‘š๐‘”๐‘™E<<2mglitalic_E < < 2 italic_m italic_g italic_l

11\displaystyle 11 =\displaystyle== 2โขmโขgโขlโขsin2โกฮธ2E+L22โขmโขl2โขE2๐‘š๐‘”๐‘™superscript2๐œƒ2๐ธsuperscript๐ฟ22๐‘šsuperscript๐‘™2๐ธ\displaystyle\frac{2mgl\sin^{2}\frac{\theta}{2}}{E}+\frac{L^{2}}{2ml^{2}E}divide start_ARG 2 italic_m italic_g italic_l roman_sin start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT divide start_ARG italic_ฮธ end_ARG start_ARG 2 end_ARG end_ARG start_ARG italic_E end_ARG + divide start_ARG italic_L start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 2 italic_m italic_l start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_E end_ARG (2)
11\displaystyle 11 =\displaystyle== mโขgโขlโขฮธ22โขE+L22โขmโขl2โขE๐‘š๐‘”๐‘™superscript๐œƒ22๐ธsuperscript๐ฟ22๐‘šsuperscript๐‘™2๐ธ\displaystyle\frac{mgl\theta^{2}}{2E}+\frac{L^{2}}{2ml^{2}E}divide start_ARG italic_m italic_g italic_l italic_ฮธ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 2 italic_E end_ARG + divide start_ARG italic_L start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 2 italic_m italic_l start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_E end_ARG (3)

It is possible to rewrite equation 2 into 3 using the small angle approximation - sinโกฮธโ‰ˆฮธ๐œƒ๐œƒ\sin\theta\approx\thetaroman_sin italic_ฮธ โ‰ˆ italic_ฮธ. This is possible noting that when E=12โขmโขฯ‰2โขl2โขฮธ2๐ธ12๐‘šsuperscript๐œ”2superscript๐‘™2superscript๐œƒ2E=\frac{1}{2}m\omega^{2}l^{2}\theta^{2}italic_E = divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_m italic_ฯ‰ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_l start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_ฮธ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT (energy of harmonic oscillator) is small, ฮธ๐œƒ\thetaitalic_ฮธ also becomes small. Physically, this condition represents a small oscillation around its stable point.

  1. (b)

    E=2โขmโขgโขl๐ธ2๐‘š๐‘”๐‘™E=2mglitalic_E = 2 italic_m italic_g italic_l

L=mโข2โขgโขl3โขcosโกฮธ2๐ฟ๐‘š2๐‘”superscript๐‘™3๐œƒ2L=m\sqrt{2gl^{3}}\cos\frac{\theta}{2}italic_L = italic_m square-root start_ARG 2 italic_g italic_l start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG roman_cos divide start_ARG italic_ฮธ end_ARG start_ARG 2 end_ARG (4)

Equation 4 may be obtained from equation 1, where we substitute the relation E=2โขmโขgโขl๐ธ2๐‘š๐‘”๐‘™E=2mglitalic_E = 2 italic_m italic_g italic_l. This marks a transition from the first phase to the second phase222Refer to the final sentence in Abstract

  1. (c)

    E>2โขmโขgโขl๐ธ2๐‘š๐‘”๐‘™E>2mglitalic_E > 2 italic_m italic_g italic_l

L2=2โขmโขl2โข[Eโˆ’2โขmโขgโขlโขsin2โกฮธ2]superscript๐ฟ22๐‘šsuperscript๐‘™2delimited-[]๐ธ2๐‘š๐‘”๐‘™superscript2๐œƒ2L^{2}=2ml^{2}[E-2mgl\sin^{2}\frac{\theta}{2}]italic_L start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = 2 italic_m italic_l start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT [ italic_E - 2 italic_m italic_g italic_l roman_sin start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT divide start_ARG italic_ฮธ end_ARG start_ARG 2 end_ARG ] (5)

Equation 5 may be obtained by simply rewriting equation 1 without making any approximations. This equation can be used to draw the phase diagram of phase 2 in the pendulum system. Physically, this condition is reached when the pendulum can perform a full rotation333All phase diagrams can be seen in section 2.7

2.2 Slack Analysis

2.2.1 Conditions

For this part only, we consider a string pendulum. This time, specific boundary conditions need to be fulfilled to perform a complete circular motion.
We consider the following constraint

mโขlโขฮธห™2+mโขgโขcosโกฮธ>0๐‘š๐‘™superscriptห™๐œƒ2๐‘š๐‘”๐œƒ0ml\dot{\theta}^{2}+mg\cos\theta>0italic_m italic_l overห™ start_ARG italic_ฮธ end_ARG start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_m italic_g roman_cos italic_ฮธ > 0

Insertion into the energy equation of the pendulum yields

E>mโขgโขlโˆ’32โขmโขgโขlโขcosโกฮธ๐ธ๐‘š๐‘”๐‘™32๐‘š๐‘”๐‘™๐œƒE>mgl-\frac{3}{2}mgl\cos\thetaitalic_E > italic_m italic_g italic_l - divide start_ARG 3 end_ARG start_ARG 2 end_ARG italic_m italic_g italic_l roman_cos italic_ฮธ

Notice how if E=52โขmโขgโขl๐ธ52๐‘š๐‘”๐‘™E=\frac{5}{2}mglitalic_E = divide start_ARG 5 end_ARG start_ARG 2 end_ARG italic_m italic_g italic_l, slack will never occur. Rewriting in terms of E, yields

cosโกฮธ>23โข[1โˆ’Emโขgโขl]๐œƒ23delimited-[]1๐ธ๐‘š๐‘”๐‘™\cos\theta>\frac{2}{3}[1-\frac{E}{mgl}]roman_cos italic_ฮธ > divide start_ARG 2 end_ARG start_ARG 3 end_ARG [ 1 - divide start_ARG italic_E end_ARG start_ARG italic_m italic_g italic_l end_ARG ] (6)

This means before the lower bound is reached, the strings will remain taut.
We may obtain the angular speed of the pendulum at this instant.

ฮธห™2=2โขg3โขlโข[โˆ’1+Emโขgโขl]superscriptห™๐œƒ22๐‘”3๐‘™delimited-[]1๐ธ๐‘š๐‘”๐‘™\dot{\theta}^{2}=\frac{2g}{3l}[-1+\frac{E}{mgl}]overห™ start_ARG italic_ฮธ end_ARG start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = divide start_ARG 2 italic_g end_ARG start_ARG 3 italic_l end_ARG [ - 1 + divide start_ARG italic_E end_ARG start_ARG italic_m italic_g italic_l end_ARG ] (7)

2.2.2 Slack Time

We find the time required by the string to become taut again. Note that the pendulum will undergo parabolic motion during this time range, and the following equation must be fulfilled for it to become taut again444Due to the complexity in the threaded lato-latoโ€™s motion, we will refrain from using it in our analysis. Instead, the focus will lie on a stick-based lato-lato..

yโ€ฒโฃ2+xโ€ฒโฃ2=l2superscript๐‘ฆโ€ฒ2superscript๐‘ฅโ€ฒ2superscript๐‘™2y^{\prime 2}+x^{\prime 2}=l^{2}italic_y start_POSTSUPERSCRIPT โ€ฒ 2 end_POSTSUPERSCRIPT + italic_x start_POSTSUPERSCRIPT โ€ฒ 2 end_POSTSUPERSCRIPT = italic_l start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT (8)

Where rโ€ฒ=(xโ€ฒ,yโ€ฒ)superscript๐‘Ÿโ€ฒsuperscript๐‘ฅโ€ฒsuperscript๐‘ฆโ€ฒr^{\prime}=(x^{\prime},y^{\prime})italic_r start_POSTSUPERSCRIPT โ€ฒ end_POSTSUPERSCRIPT = ( italic_x start_POSTSUPERSCRIPT โ€ฒ end_POSTSUPERSCRIPT , italic_y start_POSTSUPERSCRIPT โ€ฒ end_POSTSUPERSCRIPT ) is the instantaneous position of the object. We define yโ€ฒ=โˆ’lโขcosโกฮธ+v0โขtโขsinโกฮธโˆ’12โขgโขt2superscript๐‘ฆโ€ฒ๐‘™๐œƒsubscript๐‘ฃ0๐‘ก๐œƒ12๐‘”superscript๐‘ก2y^{\prime}=-l\cos\theta+v_{0}t\sin\theta-\frac{1}{2}gt^{2}italic_y start_POSTSUPERSCRIPT โ€ฒ end_POSTSUPERSCRIPT = - italic_l roman_cos italic_ฮธ + italic_v start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_t roman_sin italic_ฮธ - divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_g italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT and xโ€ฒ=lโขsinโกฮธ+v0โขtโขcosโกฮธsuperscript๐‘ฅโ€ฒ๐‘™๐œƒsubscript๐‘ฃ0๐‘ก๐œƒx^{\prime}=l\sin\theta+v_{0}t\cos\thetaitalic_x start_POSTSUPERSCRIPT โ€ฒ end_POSTSUPERSCRIPT = italic_l roman_sin italic_ฮธ + italic_v start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_t roman_cos italic_ฮธ. Insertion allows us to get the nontrivial equation

14โขgโขt2โˆ’v0โขgโขtโขsinโกฮธ+(v02+gโขlโขcosโกฮธ)=014๐‘”superscript๐‘ก2subscript๐‘ฃ0๐‘”๐‘ก๐œƒsuperscriptsubscript๐‘ฃ02๐‘”๐‘™๐œƒ0\frac{1}{4}gt^{2}-v_{0}gt\sin\theta+(v_{0}^{2}+gl\cos\theta)=0divide start_ARG 1 end_ARG start_ARG 4 end_ARG italic_g italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_v start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_g italic_t roman_sin italic_ฮธ + ( italic_v start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_g italic_l roman_cos italic_ฮธ ) = 0
gโขt=2โขv0โขsinโกฮธยฑโˆ’v02โขcos2โกฮธโˆ’g3โขlโขcosโกฮธ๐‘”๐‘กplus-or-minus2subscript๐‘ฃ0๐œƒsuperscriptsubscript๐‘ฃ02superscript2๐œƒsuperscript๐‘”3๐‘™๐œƒgt=2v_{0}\sin\theta\pm\sqrt{-v_{0}^{2}\cos^{2}\theta-g^{3}l\cos\theta}italic_g italic_t = 2 italic_v start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT roman_sin italic_ฮธ ยฑ square-root start_ARG - italic_v start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT roman_cos start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_ฮธ - italic_g start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT italic_l roman_cos italic_ฮธ end_ARG

inserting the value of cosโกฮธ=โˆ’v02gโขl๐œƒsuperscriptsubscript๐‘ฃ02๐‘”๐‘™\cos\theta=\frac{-v_{0}^{2}}{gl}roman_cos italic_ฮธ = divide start_ARG - italic_v start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_g italic_l end_ARG, will allow us to write

t=4โขv0โขsinโกฮธg๐‘ก4subscript๐‘ฃ0๐œƒ๐‘”t=\frac{4v_{0}\sin\theta}{g}italic_t = divide start_ARG 4 italic_v start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT roman_sin italic_ฮธ end_ARG start_ARG italic_g end_ARG (9)

From this point onward, the velocity component in the direction of the string will be eliminated, leaving the tangential component. This process will result in energy loss.

2.3 Pendulum Dynamics

In this part, we will be deriving the most important formula that is used in the entirety of the paper. This is the equation of motion of the pendulum, with its free end driven by an oscillating force (Fm=โˆ’a0โขฯ‰2โขcosโกฯ‰โขt๐น๐‘šsubscript๐‘Ž0superscript๐œ”2๐œ”๐‘ก\frac{F}{m}=-a_{0}\omega^{2}\cos\omega tdivide start_ARG italic_F end_ARG start_ARG italic_m end_ARG = - italic_a start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_ฯ‰ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT roman_cos italic_ฯ‰ italic_t)555a0subscript๐‘Ž0a_{0}italic_a start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT shows the amplitude of the motion undergone by the pendulumโ€™s free end and ฯ‰๐œ”\omegaitalic_ฯ‰ shows the angular frequency of the oscillatory motion. We consider the frame of the oscillating free end, such that the motion of the spheres is exactly circular.

E๐ธ\displaystyle Eitalic_E =\displaystyle== mโขgโขlโข(1โˆ’cosโกฮธ)+mโขl2โขฮธห™22+โˆซโˆ’mโขa0โขฯ‰2โขcosโกฯ‰โขtโขdโขy๐‘š๐‘”๐‘™1๐œƒ๐‘šsuperscript๐‘™2superscriptห™๐œƒ22๐‘šsubscript๐‘Ž0superscript๐œ”2๐œ”๐‘ก๐‘‘๐‘ฆ\displaystyle mgl(1-\cos\theta)+ml^{2}\frac{\dot{\theta}^{2}}{2}+\int-ma_{0}% \omega^{2}\cos\omega tdyitalic_m italic_g italic_l ( 1 - roman_cos italic_ฮธ ) + italic_m italic_l start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT divide start_ARG overห™ start_ARG italic_ฮธ end_ARG start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 2 end_ARG + โˆซ - italic_m italic_a start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_ฯ‰ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT roman_cos italic_ฯ‰ italic_t italic_d italic_y (10)
dโขEdโขt๐‘‘๐ธ๐‘‘๐‘ก\displaystyle\frac{dE}{dt}divide start_ARG italic_d italic_E end_ARG start_ARG italic_d italic_t end_ARG =\displaystyle== 00\displaystyle 0
00\displaystyle 0 =\displaystyle== mโขgโขlโขsinโกฮธโขฮธห™+mโขl2โขฮธยจโขฮธห™โˆ’mโขa0โขlโขฯ‰2โขcosโกฯ‰โขtโขsinโกฮธโขฮธห™๐‘š๐‘”๐‘™๐œƒห™๐œƒ๐‘šsuperscript๐‘™2ยจ๐œƒห™๐œƒ๐‘šsubscript๐‘Ž0๐‘™superscript๐œ”2๐œ”๐‘ก๐œƒห™๐œƒ\displaystyle mgl\sin\theta\dot{\theta}+ml^{2}\ddot{\theta}\dot{\theta}-ma_{0}% l\omega^{2}\cos\omega t\sin\theta\dot{\theta}italic_m italic_g italic_l roman_sin italic_ฮธ overห™ start_ARG italic_ฮธ end_ARG + italic_m italic_l start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT overยจ start_ARG italic_ฮธ end_ARG overห™ start_ARG italic_ฮธ end_ARG - italic_m italic_a start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_l italic_ฯ‰ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT roman_cos italic_ฯ‰ italic_t roman_sin italic_ฮธ overห™ start_ARG italic_ฮธ end_ARG
โˆ’ฮธยจยจ๐œƒ\displaystyle-\ddot{\theta}- overยจ start_ARG italic_ฮธ end_ARG =\displaystyle== (gโˆ’a0โขฯ‰2โขcosโกฯ‰โขt)โขsinโกฮธ๐‘”subscript๐‘Ž0superscript๐œ”2๐œ”๐‘ก๐œƒ\displaystyle(g-a_{0}\omega^{2}\cos\omega t)\sin\theta( italic_g - italic_a start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_ฯ‰ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT roman_cos italic_ฯ‰ italic_t ) roman_sin italic_ฮธ (11)

Equation 11 is known as the Mathieuโ€™s equations[11], having a general solution of[9]:

uโข(t)=eiโขฮผ0โขtโขp0โข(t)๐‘ข๐‘กsuperscript๐‘’๐‘–subscript๐œ‡0๐‘กsubscript๐‘0๐‘กu(t)=e^{i\mu_{0}t}p_{0}(t)italic_u ( italic_t ) = italic_e start_POSTSUPERSCRIPT italic_i italic_ฮผ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_t end_POSTSUPERSCRIPT italic_p start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_t ) (12)

In general, ฮผ0โˆˆโ„‚subscript๐œ‡0โ„‚\mu_{0}\in\mathbb{C}italic_ฮผ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT โˆˆ blackboard_C. If ฮผ0โˆˆโ„subscript๐œ‡0โ„\mu_{0}\in\mathbb{R}italic_ฮผ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT โˆˆ blackboard_R, the solution is particularly bounded.

2.4 Kapitza Model

To perform code proof-testing later on in section 4, we will consider several constraints to make equation (10) analytically solvable, which are[12]:

  1. 1.

    Small value of a0subscript๐‘Ž0a_{0}italic_a start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT

  2. 2.

    Fast oscillation frequency ฯ‰๐œ”\omegaitalic_ฯ‰

  3. 3.

    This way we may rewrite ฮธ=ฮณ+ฮฒ๐œƒ๐›พ๐›ฝ\theta=\gamma+\betaitalic_ฮธ = italic_ฮณ + italic_ฮฒ, where ฮณ๐›พ\gammaitalic_ฮณ is the slow varying term, with large amplitude and ฮฒ๐›ฝ\betaitalic_ฮฒ is the opposite of ฮณ๐›พ\gammaitalic_ฮณ

We first try to obtain the value of ฮฒ๐›ฝ\betaitalic_ฮฒ. We note that the second derivative of ฮณ๐›พ\gammaitalic_ฮณ is way smaller than that of ฮฒ๐›ฝ\betaitalic_ฮฒ. This way, we may expand the terms, hence ending up with

(gโˆ’a0โขฯ‰2โขcosโกฯ‰โขt)โขsinโกฮณ=ฮฒยจโขl๐‘”subscript๐‘Ž0superscript๐œ”2๐œ”๐‘ก๐›พยจ๐›ฝ๐‘™(g-a_{0}\omega^{2}\cos\omega t)\sin\gamma=\ddot{\beta}l( italic_g - italic_a start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_ฯ‰ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT roman_cos italic_ฯ‰ italic_t ) roman_sin italic_ฮณ = overยจ start_ARG italic_ฮฒ end_ARG italic_l (13)

noting that ฯ‰2โขa0>>gmuch-greater-thansuperscript๐œ”2subscript๐‘Ž0๐‘”\omega^{2}a_{0}>>gitalic_ฯ‰ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_a start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT > > italic_g, we end up with

ฮฒ=a0lโขsinโกฮณโขcosโกฯ‰โขt๐›ฝsubscript๐‘Ž0๐‘™๐›พ๐œ”๐‘ก\beta=\frac{a_{0}}{l}\sin\gamma\cos\omega titalic_ฮฒ = divide start_ARG italic_a start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG italic_l end_ARG roman_sin italic_ฮณ roman_cos italic_ฯ‰ italic_t (14)

Next, we iterate the obtained ฮฒ๐›ฝ\betaitalic_ฮฒ on the equation of motion to get ฮณ๐›พ\gammaitalic_ฮณ. We will neglect terms of the order y3superscript๐‘ฆ3y^{3}italic_y start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT and so on

gโขsinโกฮณ+gโขcosโกฮณโขฮฒโˆ’a0โขฮฒโขฯ‰2โขcosโกฯ‰โขtโขcosโกฮณ=โˆ’ฮณยจ๐‘”๐›พ๐‘”๐›พ๐›ฝsubscript๐‘Ž0๐›ฝsuperscript๐œ”2๐œ”๐‘ก๐›พยจ๐›พg\sin\gamma+g\cos\gamma\beta-a_{0}\beta\omega^{2}\cos\omega t\cos\gamma=-\ddot% {\gamma}italic_g roman_sin italic_ฮณ + italic_g roman_cos italic_ฮณ italic_ฮฒ - italic_a start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_ฮฒ italic_ฯ‰ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT roman_cos italic_ฯ‰ italic_t roman_cos italic_ฮณ = - overยจ start_ARG italic_ฮณ end_ARG (15)

Noting the 2nd approximation condition, we may average the previous function to obtain:

gโขsinโกฮณl+a02โขฯ‰2l2โขsinโก2โขฮณ4=โˆ’ฮณยจ๐‘”๐›พ๐‘™superscriptsubscript๐‘Ž02superscript๐œ”2superscript๐‘™22๐›พ4ยจ๐›พ\frac{g\sin\gamma}{l}+\frac{a_{0}^{2}\omega^{2}}{l^{2}}\frac{\sin 2\gamma}{4}=% -\ddot{\gamma}divide start_ARG italic_g roman_sin italic_ฮณ end_ARG start_ARG italic_l end_ARG + divide start_ARG italic_a start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_ฯ‰ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_l start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG divide start_ARG roman_sin 2 italic_ฮณ end_ARG start_ARG 4 end_ARG = - overยจ start_ARG italic_ฮณ end_ARG (16)

Next, we will find the average moment of forces acting on the pendulum

ฯ„=mโขฮณยจโขl2๐œ๐‘šยจ๐›พsuperscript๐‘™2\tau=m\ddot{\gamma}l^{2}italic_ฯ„ = italic_m overยจ start_ARG italic_ฮณ end_ARG italic_l start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT
ฯ„=โˆ’(mโขgโขlโขsinโกฮณ+mโขa02โขฯ‰2โขsinโก2โขฮณ4)๐œ๐‘š๐‘”๐‘™๐›พ๐‘šsuperscriptsubscript๐‘Ž02superscript๐œ”22๐›พ4\tau=-(mgl\sin\gamma+ma_{0}^{2}\omega^{2}\frac{\sin 2\gamma}{4})italic_ฯ„ = - ( italic_m italic_g italic_l roman_sin italic_ฮณ + italic_m italic_a start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_ฯ‰ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT divide start_ARG roman_sin 2 italic_ฮณ end_ARG start_ARG 4 end_ARG ) (17)

Then, we define a scalar potential due to this torque

Vโข(ฮณ)=โˆ’โˆซฯ„โข๐‘‘ฮณ๐‘‰๐›พ๐œdifferential-d๐›พV(\gamma)=-\int\tau d\gammaitalic_V ( italic_ฮณ ) = - โˆซ italic_ฯ„ italic_d italic_ฮณ (18)
Vโข(ฮณ)=mโขgโขlโข(1โˆ’cosโกฮณ)+mโขa02โขฯ‰2โขsin2โกฮณ4๐‘‰๐›พ๐‘š๐‘”๐‘™1๐›พ๐‘šsuperscriptsubscript๐‘Ž02superscript๐œ”2superscript2๐›พ4V(\gamma)=mgl(1-\cos\gamma)+ma_{0}^{2}\omega^{2}\frac{\sin^{2}\gamma}{4}italic_V ( italic_ฮณ ) = italic_m italic_g italic_l ( 1 - roman_cos italic_ฮณ ) + italic_m italic_a start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_ฯ‰ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT divide start_ARG roman_sin start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_ฮณ end_ARG start_ARG 4 end_ARG (19)

We continue by analyzing several stability options.

dโขVโข(ฮณ)dโขฮณ=0๐‘‘๐‘‰๐›พ๐‘‘๐›พ0\frac{dV(\gamma)}{d\gamma}=0divide start_ARG italic_d italic_V ( italic_ฮณ ) end_ARG start_ARG italic_d italic_ฮณ end_ARG = 0

from there, we obtain ฮณ1=0subscript๐›พ10\gamma_{1}=0italic_ฮณ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = 0, ฮณ2=ฯ€subscript๐›พ2๐œ‹\gamma_{2}=\piitalic_ฮณ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT = italic_ฯ€, cosโกฮณ3=โˆ’2โขgโขly02โขฯ‰2subscript๐›พ32๐‘”๐‘™superscriptsubscript๐‘ฆ02superscript๐œ”2\cos\gamma_{3}=-\frac{2gl}{y_{0}^{2}\omega^{2}}roman_cos italic_ฮณ start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT = - divide start_ARG 2 italic_g italic_l end_ARG start_ARG italic_y start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_ฯ‰ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG666We define Vmโขaโขxsubscript๐‘‰๐‘š๐‘Ž๐‘ฅV_{max}italic_V start_POSTSUBSCRIPT italic_m italic_a italic_x end_POSTSUBSCRIPT as the value of the scalar potential V๐‘‰Vitalic_V at which ฮณ=ฮณ3๐›พsubscript๐›พ3\gamma=\gamma_{3}italic_ฮณ = italic_ฮณ start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT.. The stability for ฮณ=ฯ€๐›พ๐œ‹\gamma=\piitalic_ฮณ = italic_ฯ€ only works if ฮณ3subscript๐›พ3\gamma_{3}italic_ฮณ start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT has a real solution, such that d2โขVdโขฮณ2superscript๐‘‘2๐‘‰๐‘‘superscript๐›พ2\frac{d^{2}V}{d\gamma^{2}}divide start_ARG italic_d start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_V end_ARG start_ARG italic_d italic_ฮณ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG at that point <0absent0<0< 0.

Here, we can define the energy of the system in the moving frame as

Eโ€ฒ=V+mโขl2โขฮธ2ห™2superscript๐ธโ€ฒ๐‘‰๐‘šsuperscript๐‘™2ห™superscript๐œƒ22E^{\prime}=V+\frac{ml^{2}\dot{\theta^{2}}}{2}italic_E start_POSTSUPERSCRIPT โ€ฒ end_POSTSUPERSCRIPT = italic_V + divide start_ARG italic_m italic_l start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT overห™ start_ARG italic_ฮธ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_ARG start_ARG 2 end_ARG (20)

defining L๐ฟLitalic_L,

L2=2โขmโขl2โข(Eโ€ฒโˆ’mโขgโขlโข(1โˆ’cosโกฮณ)+mโขa02โขฯ‰2โขsin2โกฮณ4)superscript๐ฟ22๐‘šsuperscript๐‘™2superscript๐ธโ€ฒ๐‘š๐‘”๐‘™1๐›พ๐‘šsuperscriptsubscript๐‘Ž02superscript๐œ”2superscript2๐›พ4L^{2}=2ml^{2}(E^{\prime}-mgl(1-\cos\gamma)+ma_{0}^{2}\omega^{2}\frac{\sin^{2}% \gamma}{4})italic_L start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = 2 italic_m italic_l start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_E start_POSTSUPERSCRIPT โ€ฒ end_POSTSUPERSCRIPT - italic_m italic_g italic_l ( 1 - roman_cos italic_ฮณ ) + italic_m italic_a start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_ฯ‰ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT divide start_ARG roman_sin start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_ฮณ end_ARG start_ARG 4 end_ARG ) (21)

2.5 Energy Loss

We consider an energy loss proportional to e2superscript๐‘’2e^{2}italic_e start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT, where e<1๐‘’1e<1italic_e < 1 represents the coefficient of restitution of the two bobs. Due to this, we can write Ensubscript๐ธ๐‘›E_{n}italic_E start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT, that is the energy to the nth collision as En=enโขE0subscript๐ธ๐‘›superscript๐‘’๐‘›subscript๐ธ0E_{n}=e^{n}E_{0}italic_E start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT = italic_e start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT italic_E start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT.

L2=2โขmโขl2โข[Enโˆ’2โขmโขgโขlโขsin2โกฮธ2]superscript๐ฟ22๐‘šsuperscript๐‘™2delimited-[]subscript๐ธ๐‘›2๐‘š๐‘”๐‘™superscript2๐œƒ2L^{2}=2ml^{2}[E_{n}-2mgl\sin^{2}\frac{\theta}{2}]italic_L start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = 2 italic_m italic_l start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT [ italic_E start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT - 2 italic_m italic_g italic_l roman_sin start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT divide start_ARG italic_ฮธ end_ARG start_ARG 2 end_ARG ]

Notice that the system will lose kinetic energy after several collisions, which means that additional energy must be given every time energy is dissipated. For every energy loss, the following ฮดโขW๐›ฟ๐‘Š\delta Witalic_ฮด italic_W must be given.

ฮดโขW=ฮ”โขEn=Enโˆ’1โข[eโˆ’1]๐›ฟ๐‘Šฮ”subscript๐ธ๐‘›subscript๐ธ๐‘›1delimited-[]๐‘’1\delta W=\Delta E_{n}=E_{n-1}[e-1]italic_ฮด italic_W = roman_ฮ” italic_E start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT = italic_E start_POSTSUBSCRIPT italic_n - 1 end_POSTSUBSCRIPT [ italic_e - 1 ] (22)

We can now see the importance of giving additional work to keep the pendulum at its original energy state. We may do this by lifting the pendulum system up and down, with the power defined as:

<P>=<Fโ†’โ‹…vโ†’>expectation๐‘ƒexpectationโ‹…โ†’๐นโ†’๐‘ฃ<P>=<\vec{F}\cdot\vec{v}>< italic_P > = < overโ†’ start_ARG italic_F end_ARG โ‹… overโ†’ start_ARG italic_v end_ARG > (23)

2.6 Tornado Play Style

Apart from the regular pendulum play style, the lato-lato can also be played less conventionally. To model this style, we refer to Figure 1. In this model, we will consider the human finger as a cylindrical wheel having radius R๐‘…Ritalic_R, where (r<<R<<l0much-less-than๐‘Ÿ๐‘…much-less-thansubscript๐‘™0r<<R<<l_{0}italic_r < < italic_R < < italic_l start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT), rotating at a constant angular velocity ฯ‰๐œ”\omegaitalic_ฯ‰.

To analyze the dynamics of the system we shall first consider the movement of the pendulum in the rotating frame ฯ‰๐œ”\omegaitalic_ฯ‰, and then we will transform the kinematic properties back into the inertial lab frame. This will ease the maths involved. The following formula for transformation will be used[13]

vโ†’lโขaโขb=v0โ†’+ฯ‰โ†’ร—rโ†’+ฯห™โ†’rsubscriptโ†’๐‘ฃ๐‘™๐‘Ž๐‘โ†’subscript๐‘ฃ0โ†’๐œ”โ†’๐‘Ÿsubscriptโ†’ห™๐œŒ๐‘Ÿ\vec{v}_{lab}=\vec{v_{0}}+\vec{\omega}\times\vec{r}+\vec{\dot{\rho}}_{r}overโ†’ start_ARG italic_v end_ARG start_POSTSUBSCRIPT italic_l italic_a italic_b end_POSTSUBSCRIPT = overโ†’ start_ARG italic_v start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG + overโ†’ start_ARG italic_ฯ‰ end_ARG ร— overโ†’ start_ARG italic_r end_ARG + overโ†’ start_ARG overห™ start_ARG italic_ฯ end_ARG end_ARG start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT (24)
aโ†’lโขaโขb=a0โ†’+ฯ‰โ†’ร—(ฯ‰โ†’ร—rโ†’)+ฯ‰โ†’ห™ร—rโ†’+2โขฯ‰โ†’ร—ฯห™โ†’r+ฯยจโ†’rsubscriptโ†’๐‘Ž๐‘™๐‘Ž๐‘โ†’subscript๐‘Ž0โ†’๐œ”โ†’๐œ”โ†’๐‘Ÿห™โ†’๐œ”โ†’๐‘Ÿ2โ†’๐œ”subscriptโ†’ห™๐œŒ๐‘Ÿsubscriptโ†’ยจ๐œŒ๐‘Ÿ\vec{a}_{lab}=\vec{a_{0}}+\vec{\omega}\times(\vec{\omega}\times\vec{r})+\dot{% \vec{\omega}}\times\vec{r}+2\vec{\omega}\times\vec{\dot{\rho}}_{r}+\vec{\ddot{% \rho}}_{r}overโ†’ start_ARG italic_a end_ARG start_POSTSUBSCRIPT italic_l italic_a italic_b end_POSTSUBSCRIPT = overโ†’ start_ARG italic_a start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG + overโ†’ start_ARG italic_ฯ‰ end_ARG ร— ( overโ†’ start_ARG italic_ฯ‰ end_ARG ร— overโ†’ start_ARG italic_r end_ARG ) + overห™ start_ARG overโ†’ start_ARG italic_ฯ‰ end_ARG end_ARG ร— overโ†’ start_ARG italic_r end_ARG + 2 overโ†’ start_ARG italic_ฯ‰ end_ARG ร— overโ†’ start_ARG overห™ start_ARG italic_ฯ end_ARG end_ARG start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT + overโ†’ start_ARG overยจ start_ARG italic_ฯ end_ARG end_ARG start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT (25)

where v0โ†’โ†’subscript๐‘ฃ0\vec{v_{0}}overโ†’ start_ARG italic_v start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG defines the velocity of the origin in the lab frame, ฯ‰โ†’โ†’๐œ”\vec{\omega}overโ†’ start_ARG italic_ฯ‰ end_ARG defines the angular velocity of the rotating frame, ฯห™โ†’rsubscriptโ†’ห™๐œŒ๐‘Ÿ\vec{\dot{\rho}}_{r}overโ†’ start_ARG overห™ start_ARG italic_ฯ end_ARG end_ARG start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT defines the velocity of the object in the rotating frame, and lastly ฯยจโ†’rsubscriptโ†’ยจ๐œŒ๐‘Ÿ\vec{\ddot{\rho}}_{r}overโ†’ start_ARG overยจ start_ARG italic_ฯ end_ARG end_ARG start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT defines the acceleration experienced by the object in the rotating frame.

2.6.1 Wrapping

For the initial part, we will derive the time taken for the pendulum to be fully retracted until l=0๐‘™0l=0italic_l = 0. This can be done by giving the sphere an initial momentum such that the string winds around our finger. By assuming the initial speed given is such that v02>>gโขlmuch-greater-thansuperscriptsubscript๐‘ฃ02๐‘”๐‘™v_{0}^{2}>>glitalic_v start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT > > italic_g italic_l, we may ignore the effects of gravity. Therefore, one could write

v0=ฮธห™โขl=cโขoโขnโขsโขtโขaโขnโขtsubscript๐‘ฃ0ห™๐œƒ๐‘™๐‘๐‘œ๐‘›๐‘ ๐‘ก๐‘Ž๐‘›๐‘กv_{0}=\dot{\theta}l=constantitalic_v start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = overห™ start_ARG italic_ฮธ end_ARG italic_l = italic_c italic_o italic_n italic_s italic_t italic_a italic_n italic_t (26)

where ฮธ๐œƒ\thetaitalic_ฮธ is the wrapping angle and l=l0โˆ’Rโขฮธ๐‘™subscript๐‘™0๐‘…๐œƒl=l_{0}-R\thetaitalic_l = italic_l start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT - italic_R italic_ฮธ. Solving the above differential equation yields the following analytical result

t=l0โขฮธโˆ’Rโขฮธ22v0๐‘กsubscript๐‘™0๐œƒ๐‘…superscript๐œƒ22subscript๐‘ฃ0t=\frac{l_{0}\theta-R\frac{\theta^{2}}{2}}{v_{0}}italic_t = divide start_ARG italic_l start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_ฮธ - italic_R divide start_ARG italic_ฮธ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 2 end_ARG end_ARG start_ARG italic_v start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG (27)

inserting ฮธ=l0r๐œƒsubscript๐‘™0๐‘Ÿ\theta=\frac{l_{0}}{r}italic_ฮธ = divide start_ARG italic_l start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG italic_r end_ARG (under the approximation r<<lmuch-less-than๐‘Ÿ๐‘™r<<litalic_r < < italic_l), will yield the final result tR=l022โขv0โขRsubscript๐‘ก๐‘…superscriptsubscript๐‘™022subscript๐‘ฃ0๐‘…t_{R}=\frac{l_{0}^{2}}{2v_{0}R}italic_t start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT = divide start_ARG italic_l start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 2 italic_v start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_R end_ARG.

2.6.2 Unwinding

This part will now use the formulas provided at the beginning of this subsection. The model used here is that we quickly rotate our finger with a constant angular velocity ฯ‰๐œ”\omegaitalic_ฯ‰. We again assume that gravity is negligible. We first notice that considering a frame rotating at angular speed ฮฉ=ฯ‰+ฮธห™ฮฉ๐œ”ห™๐œƒ\Omega=\omega+\dot{\theta}roman_ฮฉ = italic_ฯ‰ + overห™ start_ARG italic_ฮธ end_ARG will be much easier. This is because, by considering this frame, we will be given a system where the pendulumโ€™s string just changes in length, without any rotational motion. In this frame, the kinematic properties are simply :

ฯห™โ†’r=lห™โขl^=Rโขฮธห™โขl^subscriptโ†’ห™๐œŒ๐‘Ÿห™๐‘™^๐‘™๐‘…ห™๐œƒ^๐‘™\vec{\dot{\rho}}_{r}=\dot{l}\hat{l}=R\dot{\theta}\hat{l}overโ†’ start_ARG overห™ start_ARG italic_ฯ end_ARG end_ARG start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT = overห™ start_ARG italic_l end_ARG over^ start_ARG italic_l end_ARG = italic_R overห™ start_ARG italic_ฮธ end_ARG over^ start_ARG italic_l end_ARG (28)
ฯยจโ†’r=lยจโขl^=Rโขฮธห™โขl^subscriptโ†’ยจ๐œŒ๐‘Ÿยจ๐‘™^๐‘™๐‘…ห™๐œƒ^๐‘™\vec{\ddot{\rho}}_{r}=\ddot{l}\hat{l}=R\dot{\theta}\hat{l}overโ†’ start_ARG overยจ start_ARG italic_ฯ end_ARG end_ARG start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT = overยจ start_ARG italic_l end_ARG over^ start_ARG italic_l end_ARG = italic_R overห™ start_ARG italic_ฮธ end_ARG over^ start_ARG italic_l end_ARG (29)

Subsequently, insertion into equation Equation 23

aโ†’lโขaโขbsubscriptโ†’๐‘Ž๐‘™๐‘Ž๐‘\displaystyle\vec{a}_{lab}overโ†’ start_ARG italic_a end_ARG start_POSTSUBSCRIPT italic_l italic_a italic_b end_POSTSUBSCRIPT =\displaystyle== โˆ’ฮฉ2โขRโ†’โˆ’ฮฉ2โขlโ†’+lยจโขlRโขet^+2โขฯ‰โขlห™โขet^+lยจโขl^superscriptฮฉ2โ†’๐‘…superscriptฮฉ2โ†’๐‘™ยจ๐‘™๐‘™๐‘…^subscript๐‘’๐‘ก2๐œ”ห™๐‘™^subscript๐‘’๐‘กยจ๐‘™^๐‘™\displaystyle-\Omega^{2}\vec{R}-\Omega^{2}\vec{l}+\ddot{l}\frac{l}{R}\hat{e_{t% }}+2\omega\dot{l}\hat{e_{t}}+\ddot{l}\hat{l}- roman_ฮฉ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT overโ†’ start_ARG italic_R end_ARG - roman_ฮฉ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT overโ†’ start_ARG italic_l end_ARG + overยจ start_ARG italic_l end_ARG divide start_ARG italic_l end_ARG start_ARG italic_R end_ARG over^ start_ARG italic_e start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT end_ARG + 2 italic_ฯ‰ overห™ start_ARG italic_l end_ARG over^ start_ARG italic_e start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT end_ARG + overยจ start_ARG italic_l end_ARG over^ start_ARG italic_l end_ARG
(aโ†’lโขaโขb)et^subscriptsubscriptโ†’๐‘Ž๐‘™๐‘Ž๐‘^subscript๐‘’๐‘ก\displaystyle(\vec{a}_{lab})_{\hat{e_{t}}}( overโ†’ start_ARG italic_a end_ARG start_POSTSUBSCRIPT italic_l italic_a italic_b end_POSTSUBSCRIPT ) start_POSTSUBSCRIPT over^ start_ARG italic_e start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT end_ARG end_POSTSUBSCRIPT =\displaystyle== โˆ’(ฯ‰+lห™R)2โขR+lยจโขlR+2โขฮฉโขlห™superscript๐œ”ห™๐‘™๐‘…2๐‘…ยจ๐‘™๐‘™๐‘…2ฮฉห™๐‘™\displaystyle-(\omega+\frac{\dot{l}}{R})^{2}{R}+\ddot{l}\frac{l}{R}+2\Omega% \dot{l}- ( italic_ฯ‰ + divide start_ARG overห™ start_ARG italic_l end_ARG end_ARG start_ARG italic_R end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_R + overยจ start_ARG italic_l end_ARG divide start_ARG italic_l end_ARG start_ARG italic_R end_ARG + 2 roman_ฮฉ overห™ start_ARG italic_l end_ARG
(aโ†’lโขaโขb)et^subscriptsubscriptโ†’๐‘Ž๐‘™๐‘Ž๐‘^subscript๐‘’๐‘ก\displaystyle(\vec{a}_{lab})_{\hat{e_{t}}}( overโ†’ start_ARG italic_a end_ARG start_POSTSUBSCRIPT italic_l italic_a italic_b end_POSTSUBSCRIPT ) start_POSTSUBSCRIPT over^ start_ARG italic_e start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT end_ARG end_POSTSUBSCRIPT =\displaystyle== l2ห™R+lยจโขlRโˆ’ฯ‰2โขRห™superscript๐‘™2๐‘…ยจ๐‘™๐‘™๐‘…superscript๐œ”2๐‘…\displaystyle\frac{\dot{l^{2}}}{R}+\frac{\ddot{l}l}{R}-\omega^{2}Rdivide start_ARG overห™ start_ARG italic_l start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_ARG start_ARG italic_R end_ARG + divide start_ARG overยจ start_ARG italic_l end_ARG italic_l end_ARG start_ARG italic_R end_ARG - italic_ฯ‰ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_R (30)

The vector et^=etโ†’et^subscript๐‘’๐‘กโ†’subscript๐‘’๐‘กsubscript๐‘’๐‘ก\hat{e_{t}}=\frac{\vec{e_{t}}}{e_{t}}over^ start_ARG italic_e start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT end_ARG = divide start_ARG overโ†’ start_ARG italic_e start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT end_ARG end_ARG start_ARG italic_e start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT end_ARG, represents the direction perpendicular to the string l^^๐‘™\hat{l}over^ start_ARG italic_l end_ARG. Notice that the force in this direction is negligible (gravity), therefore we may immediately set it to 0. The equation can be turned into a perfect integral

0=ddโขtโข(lโขlห™)โˆ’ฯ‰2โขR20๐‘‘๐‘‘๐‘ก๐‘™ห™๐‘™superscript๐œ”2superscript๐‘…20=\frac{d}{dt}(l\dot{l})-\omega^{2}R^{2}0 = divide start_ARG italic_d end_ARG start_ARG italic_d italic_t end_ARG ( italic_l overห™ start_ARG italic_l end_ARG ) - italic_ฯ‰ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_R start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT (31)

integrating both sides by dโขt๐‘‘๐‘กdtitalic_d italic_t, and setting the initial condition l0=0subscript๐‘™00l_{0}=0italic_l start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 0, allows us to write

lโขlห™๐‘™ห™๐‘™\displaystyle l\dot{l}italic_l overห™ start_ARG italic_l end_ARG =\displaystyle== ฯ‰2โขR2โขtsuperscript๐œ”2superscript๐‘…2๐‘ก\displaystyle\omega^{2}R^{2}titalic_ฯ‰ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_R start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_t
โˆซ0llโข๐‘‘lsuperscriptsubscript0๐‘™๐‘™differential-d๐‘™\displaystyle\int_{0}^{l}l\,dlโˆซ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_l end_POSTSUPERSCRIPT italic_l italic_d italic_l =\displaystyle== 12โขฯ‰2โขR2โขt212superscript๐œ”2superscript๐‘…2superscript๐‘ก2\displaystyle\frac{1}{2}\omega^{2}R^{2}t^{2}divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_ฯ‰ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_R start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT
l๐‘™\displaystyle litalic_l =\displaystyle== ฯ‰โขRโขt๐œ”๐‘…๐‘ก\displaystyle\omega Rtitalic_ฯ‰ italic_R italic_t (32)

It turns out, the lโข(t)๐‘™๐‘กl(t)italic_l ( italic_t ) function is linear, hence we may write unwinding time tusubscript๐‘ก๐‘ขt_{u}italic_t start_POSTSUBSCRIPT italic_u end_POSTSUBSCRIPT as

tu=l0ฯ‰โขRsubscript๐‘ก๐‘ขsubscript๐‘™0๐œ”๐‘…t_{u}=\frac{l_{0}}{\omega R}italic_t start_POSTSUBSCRIPT italic_u end_POSTSUBSCRIPT = divide start_ARG italic_l start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG italic_ฯ‰ italic_R end_ARG (33)

2.6.3 Total time

We assume motion starts from unwinding, and when l๐‘™litalic_l reaches l0subscript๐‘™0l_{0}italic_l start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, the ฯ‰๐œ”\omegaitalic_ฯ‰ value is immediately set to 0. Hence, by combining the total time T=tR+tu+tt๐‘‡subscript๐‘ก๐‘…subscript๐‘ก๐‘ขsubscript๐‘ก๐‘กT=t_{R}+t_{u}+t_{t}italic_T = italic_t start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT + italic_t start_POSTSUBSCRIPT italic_u end_POSTSUBSCRIPT + italic_t start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT (third term = transition time), we can approximate the period of each โ€tornadoโ€ motion as

T=l0Rโข(1+ฯ‰โขl02โขv0)+ฯ€vet^๐‘‡subscript๐‘™0๐‘…1๐œ”subscript๐‘™02subscript๐‘ฃ0๐œ‹subscript๐‘ฃ^subscript๐‘’๐‘กT=\frac{l_{0}}{R}(1+\frac{\omega l_{0}}{2v_{0}})+\frac{\pi}{v_{\hat{e_{t}}}}italic_T = divide start_ARG italic_l start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG italic_R end_ARG ( 1 + divide start_ARG italic_ฯ‰ italic_l start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG 2 italic_v start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG ) + divide start_ARG italic_ฯ€ end_ARG start_ARG italic_v start_POSTSUBSCRIPT over^ start_ARG italic_e start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT end_ARG end_POSTSUBSCRIPT end_ARG (34)

Noting that vโ†’=โˆ’(ฯ‰โขR)โขl^+ฯ‰โขlโขet^โ†’๐‘ฃ๐œ”๐‘…^๐‘™๐œ”๐‘™^subscript๐‘’๐‘ก\vec{v}=-(\omega{R})\hat{l}+\omega l\hat{e_{t}}overโ†’ start_ARG italic_v end_ARG = - ( italic_ฯ‰ italic_R ) over^ start_ARG italic_l end_ARG + italic_ฯ‰ italic_l over^ start_ARG italic_e start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT end_ARG, we can find the velocity of the sphere when the length has reached l0subscript๐‘™0l_{0}italic_l start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, vet^=ฯ‰โขl0subscript๐‘ฃ^subscript๐‘’๐‘ก๐œ”subscript๐‘™0v_{\hat{e_{t}}}=\omega l_{0}italic_v start_POSTSUBSCRIPT over^ start_ARG italic_e start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT end_ARG end_POSTSUBSCRIPT = italic_ฯ‰ italic_l start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT. The radial velocity can be ignored because when the string quickly goes back to being taut - noting (R<<lmuch-less-than๐‘…๐‘™R<<litalic_R < < italic_l) - the radial component just vanishes. v0subscript๐‘ฃ0v_{0}italic_v start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT defined previously is also equal to vet^subscript๐‘ฃ^subscript๐‘’๐‘กv_{\hat{e_{t}}}italic_v start_POSTSUBSCRIPT over^ start_ARG italic_e start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT end_ARG end_POSTSUBSCRIPT, due to the periodicity defined. Our expression simplifies into:

T=3โขl02โขR+ฯ€ฯ‰๐‘‡3subscript๐‘™02๐‘…๐œ‹๐œ”T=\frac{3l_{0}}{2R}+\frac{\pi}{\omega}italic_T = divide start_ARG 3 italic_l start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG 2 italic_R end_ARG + divide start_ARG italic_ฯ€ end_ARG start_ARG italic_ฯ‰ end_ARG (35)

2.7 Preliminary Figures

Refer to caption
Figure 2: The orange contours in the plots show the phase of the pendulum system, whereas the green contour shows the potential barrier experienced by the pendulum in the kapitza model. Figures (a),(b), and (c) show the phase diagram of the special cases, where its equations are given in subsection 1, representing a simple mathematical pendulum. Meanwhile, figures (d), (e),(f), and (g), represent the phase diagram for the kapitza model, that analyzes the conditions 0<Eโ€ฒ<2โขmโขgโขl0superscript๐ธโ€ฒ2๐‘š๐‘”๐‘™0<E^{\prime}<2mgl0 < italic_E start_POSTSUPERSCRIPT โ€ฒ end_POSTSUPERSCRIPT < 2 italic_m italic_g italic_l, 2โขmโขgโขl<Eโ€ฒ<Vmโขaโขx2๐‘š๐‘”๐‘™superscript๐ธโ€ฒsubscript๐‘‰๐‘š๐‘Ž๐‘ฅ2mgl<E^{\prime}<V_{max}2 italic_m italic_g italic_l < italic_E start_POSTSUPERSCRIPT โ€ฒ end_POSTSUPERSCRIPT < italic_V start_POSTSUBSCRIPT italic_m italic_a italic_x end_POSTSUBSCRIPT, Eโ€ฒ=Vmโขaโขxsuperscript๐ธโ€ฒsubscript๐‘‰๐‘š๐‘Ž๐‘ฅE^{\prime}=V_{max}italic_E start_POSTSUPERSCRIPT โ€ฒ end_POSTSUPERSCRIPT = italic_V start_POSTSUBSCRIPT italic_m italic_a italic_x end_POSTSUBSCRIPT, and Eโ€ฒ>Vmโขaโขxsuperscript๐ธโ€ฒsubscript๐‘‰๐‘š๐‘Ž๐‘ฅE^{\prime}>V_{max}italic_E start_POSTSUPERSCRIPT โ€ฒ end_POSTSUPERSCRIPT > italic_V start_POSTSUBSCRIPT italic_m italic_a italic_x end_POSTSUBSCRIPT respectively. The orange contour in the four figures represents a plot of the systemโ€™s potential barrier, whereas the red contour shows the phase diagram Eโ€ฒsuperscript๐ธโ€ฒE^{\prime}italic_E start_POSTSUPERSCRIPT โ€ฒ end_POSTSUPERSCRIPT vs ฮธ๐œƒ\thetaitalic_ฮธ of the kapitza model. Finally, (h) and (i) show the phase transitions that occur for a simple pendulum and the kapitza model. Both are done at the most extreme case of each respective condition.

In this part, we have provided the plots of trivial single pendulum phase diagrams, as well as other figures that might aid in illustrating Section 2.

3 Application in Lato-lato

Having derived all the necessary equations we may now solve the equation of motion (equation 10) numerically. Notice how this equation shows how the lato-lato is usually moved around by the player. Since equation 10 only shows the motion of a single pendulum, it is required that we add specific constraints, characterizing the geometry of a double pendulum system. The constraints are as follows: We apply 2 boundary conditions at ฮธ=0๐œƒ0\theta=0italic_ฮธ = 0 and ฮธ=ฯ€๐œƒ๐œ‹\theta=\piitalic_ฮธ = italic_ฯ€. Forcing, ฮธห™(0)subscriptห™๐œƒ0\dot{\theta}_{(0)}overห™ start_ARG italic_ฮธ end_ARG start_POSTSUBSCRIPT ( 0 ) end_POSTSUBSCRIPT and ฮธห™(ฯ€)subscriptห™๐œƒ๐œ‹\dot{\theta}_{(\pi)}overห™ start_ARG italic_ฮธ end_ARG start_POSTSUBSCRIPT ( italic_ฯ€ ) end_POSTSUBSCRIPT to always be reversed.

(ฮธห™0)i=โˆ’(ฮธห™0)i+1subscriptsubscriptห™๐œƒ0๐‘–subscriptsubscriptห™๐œƒ0๐‘–1(\dot{\theta}_{0})_{i}=-(\dot{\theta}_{0})_{i+1}( overห™ start_ARG italic_ฮธ end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT = - ( overห™ start_ARG italic_ฮธ end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) start_POSTSUBSCRIPT italic_i + 1 end_POSTSUBSCRIPT (36)
(ฮธห™ฯ€)i=โˆ’(ฮธห™ฯ€)i+1subscriptsubscriptห™๐œƒ๐œ‹๐‘–subscriptsubscriptห™๐œƒ๐œ‹๐‘–1(\dot{\theta}_{\pi})_{i}=-(\dot{\theta}_{\pi})_{i+1}( overห™ start_ARG italic_ฮธ end_ARG start_POSTSUBSCRIPT italic_ฯ€ end_POSTSUBSCRIPT ) start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT = - ( overห™ start_ARG italic_ฮธ end_ARG start_POSTSUBSCRIPT italic_ฯ€ end_POSTSUBSCRIPT ) start_POSTSUBSCRIPT italic_i + 1 end_POSTSUBSCRIPT (37)

This represents the almost elastic collisions between the plastic spheres, which immediately models the momentum conservation law between the plastic spheres. The codes needed to solve the equation numerically and to obtain plots of the system are all provided by chatGPT. These codes can be seen in the Appendices.

4 Code Proof-testing in Kapitza Model

To test the performance of the code, we simulate a special case called the Kapitza model, as in section 2.4 To do this, we just need to use code (b) (ฮธ๐œƒ\thetaitalic_ฮธ as a function of time) and insert some numeric constants which are in agreement with the approximations in section 2.4. It is also important to note that the equation of motion - equation 11 - of the pendulum is strictly independent of the bobโ€™s mass, hence for convenience, we may as well set m=0.020๐‘š0.020m=0.020italic_m = 0.020kg. In addition, to cover all ranges of mass, we define the specific energy Eโ€ฒsuperscript๐ธโ€ฒE^{\prime}italic_E start_POSTSUPERSCRIPT โ€ฒ end_POSTSUPERSCRIPT as the energy per unit mass (kg). We use the following numeric values:

  1. (a)

    a0=0.01subscript๐‘Ž00.01a_{0}=0.01italic_a start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 0.01 m

  2. (b)

    l=0.25๐‘™0.25l=0.25italic_l = 0.25 m

  3. (c)

    g=9.81๐‘”9.81g=9.81italic_g = 9.81 m/s2๐‘šsuperscript๐‘ 2m/s^{2}italic_m / italic_s start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT

  4. (d)

    ฯ‰=100๐œ”100\omega=100italic_ฯ‰ = 100 rad/s

  5. (e)

    ฮธ(0)=1.5subscript๐œƒ01.5\theta_{(0)}=1.5italic_ฮธ start_POSTSUBSCRIPT ( 0 ) end_POSTSUBSCRIPT = 1.5 rad

  6. (f)

    ฮธห™(0)=0subscriptห™๐œƒ00\dot{\theta}_{(0)}=0overห™ start_ARG italic_ฮธ end_ARG start_POSTSUBSCRIPT ( 0 ) end_POSTSUBSCRIPT = 0 rad/s

Refer to caption
Figure 3: We notice the ฮธ๐œƒ\thetaitalic_ฮธ may be decomposed just as our approximation in section 2.4 (ฮธโ‰ˆฮณ+ฮฒ๐œƒ๐›พ๐›ฝ\theta\approx\gamma+\betaitalic_ฮธ โ‰ˆ italic_ฮณ + italic_ฮฒ)

From here, it is safe to say that the codes provided should work just fine

5 Results in Lato-lato

Now that we have the codes that represent the lato-latoโ€™s mechanics, we can analyze the reason why playing the lato-lato is difficult. In general, we have observed that there are 2 main phases in the systemโ€™s motion. To analyze we will use the 4 codes we have attached in GitHub[14].

5.1 Maximum Amplitude

The code provided in [14] will plot a graph of the maximum amplitude vs ฯ‰๐œ”\omegaitalic_ฯ‰. From this graph, we can see that there are 2 phases. The first phase represents the normal pendulum motion and the second phase represents the condition at which the spheres collide at both ฮธ=0๐œƒ0\theta=0italic_ฮธ = 0 and ฮธ=ฯ€๐œƒ๐œ‹\theta=\piitalic_ฮธ = italic_ฯ€ rads.

Refer to caption
Figure 4: The maximum amplitude is only considered within the 3-second range after the motion starts. This is done to imitate the way people usually play the lato-lato. The figures are plotted for a0=0.05subscript๐‘Ž00.05a_{0}=0.05italic_a start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 0.05 m, a0=0.1subscript๐‘Ž00.1a_{0}=0.1italic_a start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 0.1 m, a0=0.5subscript๐‘Ž00.5a_{0}=0.5italic_a start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 0.5 m and a0=1subscript๐‘Ž01a_{0}=1italic_a start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 1 m.

The three graphs represent the plots for different a0subscript๐‘Ž0a_{0}italic_a start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT. The graph in the upper left corner shows the plot for a0=0.5subscript๐‘Ž00.5a_{0}=0.5italic_a start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 0.5 m, the one in the upper right corner is the zoomed version of the latter, the one in the lower left shows for a0=0.1subscript๐‘Ž00.1a_{0}=0.1italic_a start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 0.1 m, and in the lower right corner shows for a0=1subscript๐‘Ž01a_{0}=1italic_a start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 1 m. Notice how different a0subscript๐‘Ž0a_{0}italic_a start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT yield different maximum ฯ‰๐œ”\omegaitalic_ฯ‰ graphs.
Generally, one may observe that when the driven amplitude a0subscript๐‘Ž0a_{0}italic_a start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT is increased, the minimum driven angular frequency ฯ‰๐œ”\omegaitalic_ฯ‰ needed to reach phase 2 is lower. For lower amplitudes on the other hand, specific ฯ‰๐œ”\omegaitalic_ฯ‰ should be maintained to reach phase 2. This strengthens the argument as to why the lato-lato is a difficult game.

5.2 ฮธ๐œƒ\thetaitalic_ฮธ as a function of time

In this part, we can see the graph that represents the polar coordinates of the lato-lato as a function of time. Through these plots, we should be able to get a general feel of the systemโ€™s motion.

   Refer to caption โ€… Refer to caption

(a)                                        (b)

Figure 5: Figure on the left shows ฮธ๐œƒ\thetaitalic_ฮธ as a time function for ฯ‰=1๐œ”1\omega=1italic_ฯ‰ = 1 rad/s and the one on the right shows for ฯ‰=10๐œ”10\omega=10italic_ฯ‰ = 10 rad/s. For the smaller ฯ‰๐œ”\omegaitalic_ฯ‰, we have a sinusoidal amplitude function and a less periodic wave function for ฯ‰>10๐œ”10\omega>10italic_ฯ‰ > 10 rad/s

Figure 4 shows the ฮธ๐œƒ\thetaitalic_ฮธ vs time function for a0=0.5subscript๐‘Ž00.5a_{0}=0.5italic_a start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 0.5 m. In general, we can see that increasing the ฯ‰๐œ”\omegaitalic_ฯ‰ whilst keeping the a0subscript๐‘Ž0a_{0}italic_a start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT constant will help the lato-lato reach the 2nd phase. However, this will not be the case if a0subscript๐‘Ž0a_{0}italic_a start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT is too small.
At lower angular frequencies, we may observe an amplitude undergoing sinusoidal change. Furthermore, at higher angular frequencies, the second phase is easily reached, though the motion does not seem as periodic as the latter.

5.3 Energy as a function of time and ฮธ๐œƒ\thetaitalic_ฮธ

Refer to caption
Figure 6: Figures on the left show one of the pendulumโ€™s energy as a time function and ฮธ๐œƒ\thetaitalic_ฮธ for ฯ‰=1๐œ”1\omega=1italic_ฯ‰ = 1 rad/s and the on the right shows for ฯ‰=10๐œ”10\omega=10italic_ฯ‰ = 10 rad/s. For the smaller ฯ‰๐œ”\omegaitalic_ฯ‰, we have a periodic graph, whereas the larger one shows a rather irregular progression. All the graphs are plotted at a0=0.5subscript๐‘Ž00.5a_{0}=0.5italic_a start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 0.5 m

Through the plots in this section, we should be able to understand the energy transfer777Energy transfer here, refers to how much work is being transferred from the personโ€™s hand to the lato-lato system. that occurs between the system and the driving force. We will yield plots that show how the energy transfer varies as a function of the driven ฯ‰๐œ”\omegaitalic_ฯ‰. Consequently, one may compare how the energy of the system relates to the lato-latoโ€™s position.

First, it is important to notice that at lower angular frequency ฯ‰๐œ”\omegaitalic_ฯ‰, the energy indeed has an amplitude varying periodically, similar to the ฮธ๐œƒ\thetaitalic_ฮธ plot we have obtained. Likewise for higher ฯ‰๐œ”\omegaitalic_ฯ‰, a fluctuating pattern can be seen. Moreover, notice how the energy transfer is larger for a higher ฯ‰๐œ”\omegaitalic_ฯ‰. The energy reached at ฮธ=0๐œƒ0\theta=0italic_ฮธ = 0 may yield relatively high values, reaching a maximum of Emโขaโขx=17.60subscript๐ธ๐‘š๐‘Ž๐‘ฅ17.60E_{max}=17.60italic_E start_POSTSUBSCRIPT italic_m italic_a italic_x end_POSTSUBSCRIPT = 17.60 J/kg and minimum Emโขiโขn=2.23subscript๐ธ๐‘š๐‘–๐‘›2.23E_{min}=2.23italic_E start_POSTSUBSCRIPT italic_m italic_i italic_n end_POSTSUBSCRIPT = 2.23 J/kg for ฯ‰=1๐œ”1\omega=1italic_ฯ‰ = 1 rad/s. A less orderly graph can be seen, with a maximum energy of Emโขaโขx=86.88subscript๐ธ๐‘š๐‘Ž๐‘ฅ86.88E_{max}=86.88italic_E start_POSTSUBSCRIPT italic_m italic_a italic_x end_POSTSUBSCRIPT = 86.88 J/kg and Emโขiโขn=0.01subscript๐ธ๐‘š๐‘–๐‘›0.01E_{min}=0.01italic_E start_POSTSUBSCRIPT italic_m italic_i italic_n end_POSTSUBSCRIPT = 0.01 J/kg for the higher ฯ‰=10๐œ”10\omega=10italic_ฯ‰ = 10 rad/s. This explains why it is easier for the system to reach the 2nd phase when given a higher angular frequency, recalling that the energy transfer is higher in the high angular frequency case.

5.4 Systemโ€™s stability

The mesh plots provided in this section will directly show us the phase conditions in the lato-latoโ€™s phase diagram, that is whether it reaches the 2nd phase or not. The mesh plots will give a diagram representing the amplitude conditions vs the initial boundary conditions ฮธ0subscript๐œƒ0\theta_{0}italic_ฮธ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT and ฮธห™0subscriptห™๐œƒ0\dot{\theta}_{0}overห™ start_ARG italic_ฮธ end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT. In this section, we will be able to analyze which sort of amplitudes will give the easiest condition required to reach the 2nd phase of the motion. We will compare three driven amplitudes, which are (1) a0=0.1โขmsubscript๐‘Ž00.1๐‘ša_{0}=0.1mitalic_a start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 0.1 italic_m, (2) a0=0.5โขmsubscript๐‘Ž00.5๐‘ša_{0}=0.5mitalic_a start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 0.5 italic_m, and (3) a0=1โขmsubscript๐‘Ž01๐‘ša_{0}=1mitalic_a start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 1 italic_m.

5.4.1 a0=0.05โขmsubscript๐‘Ž00.05๐‘ša_{0}=0.05mitalic_a start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 0.05 italic_m

We can see the yellow dot reaches the 2nd phase at around ฯ‰=13โˆ’16๐œ”1316\omega=13-16italic_ฯ‰ = 13 - 16 rad/s888People usually play the lato-lato at this angular frequency, equivalent to a frequency of about 2-3 Hz). From the mesh plots, we can also see more unstable999โ€unstableโ€ refers to not being able to reach phase 2 (colored red). regions (there is a larger blue area in between the red ones), especially at around 15-18 rad/s. Due to the irregularity in the region distribution, we may infer that relatively small changes in the angular frequency can cause huge changes in the systemโ€™s mesh plots. For the purpose of our analysis, we will focus on the transition experienced by the yellow dot at around ฯ‰=16.25๐œ”16.25\omega=16.25italic_ฯ‰ = 16.25 rad/s, and ฯ‰=13๐œ”13\omega=13italic_ฯ‰ = 13 rad/s.
Taking the average of the lower (13 rad/s) and upper limit (16.25 rad/s) of the typical ฯ‰๐œ”\omegaitalic_ฯ‰ values, yields ฯ‰=14.625๐œ”14.625\omega=14.625italic_ฯ‰ = 14.625 rad/s. This gives a ฮดlโขoโขwsubscript๐›ฟ๐‘™๐‘œ๐‘ค\delta_{low}italic_ฮด start_POSTSUBSCRIPT italic_l italic_o italic_w end_POSTSUBSCRIPT101010Defined as ฮ”โขฯ‰ฯ‰ฮ”๐œ”๐œ”\frac{\Delta\omega}{\omega}divide start_ARG roman_ฮ” italic_ฯ‰ end_ARG start_ARG italic_ฯ‰ end_ARG. ฮ”โขฯ‰ฮ”๐œ”\Delta\omegaroman_ฮ” italic_ฯ‰ is calculated as ฯ‰fโˆ’ฯ‰0subscript๐œ”๐‘“subscript๐œ”0\omega_{f}-\omega_{0}italic_ฯ‰ start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT - italic_ฯ‰ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, where ฯ‰fsubscript๐œ”๐‘“\omega_{f}italic_ฯ‰ start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT represents the ฯ‰๐œ”\omegaitalic_ฯ‰ at which the point experiences a transition in phase. A lower value of ฮด๐›ฟ\deltaitalic_ฮด shows a larger sensitivity of 12.5%percent\%% and ฮดhโขiโขgโขhsubscript๐›ฟโ„Ž๐‘–๐‘”โ„Ž\delta_{high}italic_ฮด start_POSTSUBSCRIPT italic_h italic_i italic_g italic_h end_POSTSUBSCRIPT of about 11.1%percent\%%. At this ฯ‰๐œ”\omegaitalic_ฯ‰, we get a relatively stable condition. Therefore, to reach phase 2 in the lato-latoโ€™s motion, one has to maintain a specific angular frequency and initial conditions.

Here, we have made the assumption that the changes are done at the yellow dot. This means if the error is done at random initial conditions ฮธ0subscript๐œƒ0\theta_{0}italic_ฮธ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT and ฮธห™0subscriptห™๐œƒ0\dot{\theta}_{0}overห™ start_ARG italic_ฮธ end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, for instance at (1111 rad, โˆ’1.51.5-1.5- 1.5 rad/s ), the system will no longer be able to reach the desired phase 2, in 3 seconds.

Refer to caption
Figure 7: Mesh plot progression for a0=0.05โขmsubscript๐‘Ž00.05๐‘ša_{0}=0.05mitalic_a start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 0.05 italic_m, a0=0.5โขmsubscript๐‘Ž00.5๐‘ša_{0}=0.5mitalic_a start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 0.5 italic_m and a0=1โขmsubscript๐‘Ž01๐‘ša_{0}=1mitalic_a start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 1 italic_m showing Amplitude Conditions for different ฯ‰๐œ”\omegaitalic_ฯ‰. Phase 1 in the mesh plots is denoted by the color blue and phase 2 by the color red. The yellow dot shows the observed condition, that is at ฮธ0=2subscript๐œƒ02\theta_{0}=2italic_ฮธ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 2 rad and ฮธ0ห™=โˆ’0.5ห™subscript๐œƒ00.5\dot{\theta_{0}}=-0.5overห™ start_ARG italic_ฮธ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG = - 0.5 rad/s

5.4.2 a0=0.5โขmsubscript๐‘Ž00.5๐‘ša_{0}=0.5mitalic_a start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 0.5 italic_m

Notice at a very high ฯ‰๐œ”\omegaitalic_ฯ‰, the system will reach the 2nd phase regardless of the boundary conditions. Focusing on the yellow dot, it can be seen the system starts to enter the 2nd phase at ฯ‰=6.1๐œ”6.1\omega=6.1italic_ฯ‰ = 6.1rad/s and ฯ‰=6.5๐œ”6.5\omega=6.5italic_ฯ‰ = 6.5 rad/s. Subsequently, it enters the 2nd phase again at ฯ‰=8๐œ”8\omega=8italic_ฯ‰ = 8rad/s, from which it stays there for any ฯ‰๐œ”\omegaitalic_ฯ‰ larger than the former. If we consider the ฯ‰=6.3๐œ”6.3\omega=6.3italic_ฯ‰ = 6.3 rad/s, the ฮด๐›ฟ\deltaitalic_ฮด towards ฯ‰=6.1๐œ”6.1\omega=6.1italic_ฯ‰ = 6.1 rad/s is about 3.8%percent3.83.8\%3.8 % and 3.17%percent3.173.17\%3.17 % towards ฯ‰=6.5๐œ”6.5\omega=6.5italic_ฯ‰ = 6.5 rad/s. This shows that at the hill-like contour, the system is not very stable when compared with the case (a0,ฯ‰)=(0.05(a_{0},\omega)=(0.05( italic_a start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , italic_ฯ‰ ) = ( 0.05 m, 14.6314.6314.6314.63 rad/s). However, specific initial conditions, such as (3,1.5), allow a very stable condition, as shown in the mesh plot progression.

5.4.3 a0=1โขmsubscript๐‘Ž01๐‘ša_{0}=1mitalic_a start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 1 italic_m

Even though the numeric value for this consideration may be considered unrealistic, we may still get a decent analysis by considering this case. Notice how the small budge in Figure 3 now reaches the maximum ฮธ=3.14๐œƒ3.14\theta=3.14italic_ฮธ = 3.14 rad. This means that when we increase the driven amplitude a0subscript๐‘Ž0a_{0}italic_a start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, it becomes easier to reach the 2nd phase of the lato-latoโ€™s motion. The yellow dot quickly reaches the 2nd phase at ฯ‰=5๐œ”5\omega=5italic_ฯ‰ = 5 rad/s, and for all ฯ‰๐œ”\omegaitalic_ฯ‰ bigger than 5, it stays in the 2nd phase. Moreover, it becomes clearer that the considered initial condition reaches the โ€œall-redโ€ part of the diagram more quickly, as expected when the amplitude is increased.

Therefore, it is better to just play at a high frequency, because that way we can guarantee a larger energy transfer111111Energy transfer here, refers to the energy transfer from the playerโ€™s hand to the lato-lato system (stable at phase 2).

6 Conclusion

Having analyzed the conditions for reaching phase 2 of the lato-latoโ€™s motion, we may finally infer that there are certain driven angular frequencies and amplitudes required to maintain its motion. These characteristics may be sensitive as observed in Section 5.4.1, for a0=0.05โขmsubscript๐‘Ž00.05๐‘ša_{0}=0.05mitalic_a start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 0.05 italic_m. Even at a very specific ฯ‰=14.63๐œ”14.63\omega=14.63italic_ฯ‰ = 14.63 rad/s, we only get a ฮด๐›ฟ\deltaitalic_ฮด of as much as 12%, which can be considered extremely sensitive. As a result, subtle changes at the incorrect moment may prevent continuous collisions between the bobs. Through the analysis of phase diagrams, we have qualitatively proven that playing the lato-lato is indeed difficult. Furthermore, throughout this paper, we have only considered the dynamics of a stick-based pendulum, which is easier to control than the one made of string. Therefore, trained muscle memory and experience are some of the necessities required if we wish to master this game.

\bmhead

Acknowledgments We would like to deliver our deepest appreciation to Sandy Adhitia Ekahana for the spark of ideas[15]. Further, this paper would not have been made possible without the codes required for the numerical analysis which are provided by the AI ChatGPT.

References