Numerical analysis of grasshopper escapement David Ziemkiewicz∗ Institute of Mathematics and Physics, UTP University of Science and Technology, Aleje Prof. S. Kaliskiego 7, 85-789 Bydgoszcz, Poland. The dynamics of driven, damped pendulum as used in mechanical clocks is numerically investi- gated. In addition to the analysis of a well-known mechanisms such as chronometer escapement, the unusual properties of Harrison’s grasshopper escapement are explored, giving some insights regard- ing the dynamics of this system. Both steady state operation and transient effects are discussed, indicating the optimal condition for stable long-term clock accuracy. Possibility of chaotic motion is investigated. I. INTRODUCTION proach to analyse periodic, nonlinear systems [10]. On the other hand, recent advancements in computer tech- A weight-driven pendulum clock is a nonlinear, dy- nology make a direct numerical integration of the partial namic system consisting of a damped, driven pendulum differential equations of pendulum motion with resolution and the so-called escapement mechanism. The role of the much smaller than single period a practical approach, mechanism is twofold. It regulates the speed of the clock even for extended timescales spanning weeks or months by binding it to the period of the pendulum. The second [9]. The advantage of such calculation lies in its flexi- function is to provide the energy to the pendulum, so bility and amount of provided data; the fine structure that a nonzero oscillation amplitude can be maintained of the motion on a time scale smaller than the period is despite the friction. The energy is added by pushing the preserved. Additionally, transient effects can be readily pendulum along some part of its motion, geometrically studied. determined by escapement construction. Thus, the fre- Finally, chaotic dynamics of the whole pendulum- quency of the force is always the same as the pendulum escapement system is explored. The pendulum is an ex- itself, making it a self-excited oscillator [1]. Extensive cellent tool to study chaotic motion; such a behaviour overview of the physics of pendulum clocks is presented may emerge in systems with mass suspended on elas- in [2]. tic string [11], periodically driven [12], parametrically The long, historical struggle to attain highest possi- damped [13] and double [14] pendulums. Strongly ble accuracy (stability of the period) consists of devising damped pendulums may exhibit symmetry breaking [15]. new ways to separate the pendulum from various dis- In a clock, the mechanical oscillations of the mechanism turbances, including the effect of escapement itself; the coupled to the pendulum provide a nonlinear system with same mechanism that keeps the pendulum in motion is chaotic dynamics [16]. also responsible for instability. Accuracy of many popu- lar escapement types has been extensively analysed; some recent studies include dead-beat escapement [3], gravity escapement [3] and anchor escapement [1, 4]. Some re- II. CIRCULAR ERROR search covers mechanisms of historical significance, but relatively poor timekeeping properties such as verge and foliot [5, 6]. The closely related field of dynamics of me- The equation of motion for a pendulum modelled as chanical watches is still dynamically developing [7]. By a point mass m, suspended on a weightless string with moving from an idealized pendulum to physical one, a length L is large number of factors need to be taken into account [8], which can be roughly divided into three categories de- M g mα¨ − γα_ = = sin(α): (1) pending on the source of error: pendulum, environment I L and escapement mechanism. The focus of this paper is the interaction between pendulum circular error present α γ arXiv:2103.14947v2 [physics.class-ph] 7 Apr 2021 where is the pendulum angle, is a damping con- for any nonzero swing angle and escapement error caused stant, M = mgL sin(α) is the total moment of force and by the disturbance of the pendulum motion by the clock I = mL2 is the moment of inertia. In a physical pen- mechanism. dulum, M and I have a more complicated forms. By In [9] the author used an approach that averages the substituting sin(α) ≈ α, one obtains an equation for har- force over period to study the long-term effects on the p monic oscillator with period T0 = 2π L=g and angle pendulum, on a time scales much longer than single pe- α = A sin(!t + φ0), where ! = 2π=T , A is the amplitude riod. The method of averaging is a well-developed ap- and φ0 is the starting phase. Exact solution of Eq. (1) involves a complete elliptic integral [18]. The period of the pendulum T is always larger than T0 and the differ- ∗ [email protected] ence is called the circular error. The value of T can be 2 expressed as a series [19] escapement mechanism of a clock needs to provide the lost energy ∆E in every period by doing work " 2 2 # 1 2 A 1 · 3 4 A T = T 1 + sin + sin + ::: : α2 0 2 2 2 · 4 2 Z ∆E = MF (α)dα (4) (2) α As a first approximation, the error is proportional to 1 the square of amplitude. Another approach, which will where MF is the moment of force (torque) exerted on the be used in this manuscript, is the Arithmetic-Geometric pendulum and hα1; α2i is the range of pendulum angle mean (AGM) [18], which provides particularly accurate where the force is applied. However, any force acting estimation of circular error. Both methods are shown on upon the pendulum will change its rate, introducing an the Fig. 1 along the results of numerical simulation of escapement error EE. pendulum motion. There is an excellent agreement be- T 0 − T EE = ; (5) T where T and T 0 are periods of the free and disturbed pendulum. With the above convention, positive EE indi- cates larger T 0 and slower clock speed. In general, there are two distinct cases of torque adding energy to the pen- dulum: • force acting before the pendulum reaches its lowest point α = 0, directed towards the lowest point, • force acting after the pendulum reaches its lowest point α = 0, directed away from the lowest point. In the first case, the torque adds to the restoring force FIG. 1. The analytically calculated and numerically simu- generated by gravity; therefore, the pendulum acts as lated pendulum period as a function of amplitude though the gravity force was greater and speeds up. Torque applied after α = 0 reduces the "effective gravity" tween the theory and simulation results. Detailed simu- and slows the pendulum down. One can use a symmetric lation description and accuracy estimations are presented range of α to avoid influencing the pendulum speed; in in the Appendix A. Naturally, to minimize the influence other words, the period is constant if the force is an even of circular error one should aim for the smallest practi- function of α [1]. This is the design goal of a chronometer cal swing angle. However, even for A = 2 degrees, the escapement. −5 circular error EC ∼ 7 · 10 which corresponds to a non- Intuitively, to reduce the error EE one needs to min- trivial change of rate on the order of 6 seconds/day. One imise the force MF . This requirement demands low losses of the ways of correcting the circular arc is modifying and high Q factor. Such result is in agreement with obser- the pendulum suspension in such a way that it follows a vations in other systems such as atomic clocks, Q factor cycloid instead of circular arc, as shown by Hyugens in is a measure of stability (precision) of an oscillator [21]. 1673 [20]. Interestingly, as it will be shown, the elimina- However, some amount of friction is necessary for the tion of circular error is not always beneficial to the clock pendulum to stabilise at some amplitude A. Moreover, accuracy. as it will be shown, a pendulum with stronger damping reaches the equilibrium point faster, making it more re- sistant to random external disturbances. Therefore, in III. ESCAPEMENT ERROR the case of mechanical clock, the choice of the optimal Q factor is nontrivial. Any physical pendulum is subject to friction, which causes reduction of amplitude over time. The common IV. CHRONOMETER ESCAPEMENT parameter describing the energy loss of an oscillator is the Q factor Lets assume a L = 1 m, m = 1 kg pendulum driven 2πE m! by a constant torque MF (α1 < α < α2) = const applied Q = = ; (3) ∆E γ when the velocity α_ is positive. The friction is assumed to be proportional to speed, e.g. Mfriction = Lγα_ , where where E is the total energy and ∆E is the energy lost γ is the damping coefficient. Fig. 2 a) shows the typical in one period. In order to sustain pendulum motion, the case where the torque is applied exactly in the middle of 3 One can also notice the lack of data points for A < α2, in- dicating that the force was insufficient to sustain a stable oscillation. Finally, Fig. 2 c) shows a case where the pen- dulum is pushed before reaching the equilibrium point. As expected, the escapement error is negative (period is reduced) and the total error has no local minimum. How- ever, ET (A) is almost linear relation in a large range of A, which could be taken advantage of when designing ad- ditional error correction mechanism. In both asymmetric cases, the absolute value of EE decreases with amplitude. a) This can be explained by the fact that given fixed values of α1; α2, as the A increases, smaller portion of the swing arc is affected by the escapement. This result is consis- tent with findings in [9], where for a short push taking tan α1 place at angle α1, the escapement error is EE ∼ A2 .