# Simple Harmonic Motion

Copyright © David Boettcher 2005 - 2024 all rights reserved.Simple Harmonic Motion occurs in systems involving a mass coupled to a spring. When the mass is displaced from its rest position and released, it oscillates around the rest position at a “natural frequency” determined by the magnitude of the mass and the stiffness of the spring. If the position of the mass is plotted against elapsed time, wave-like curves such as those in the right hand side of the figure below are produced. These are called sine waves because they can be drawn by the trigonometric sine function.

The oscillator at the heart of a mechanical watch comprises a mass, the balance, coupled to a spring, the balance spring. Simple Harmonic Motion illustrates how it is able to keep time despite the amplitude of the oscillation changing, which occurs as the watch is put into different positions, dial up, pendant down etc., and as the main spring runs down. The natural frequency does not change with changes in amplitude. This just what is needed, because it endows a watch with isochronism, the ability to keep the same time whatever the amplitude of the balance.

Note that Simple Harmonic Motion is a simplified ideal that does not really exist. Real balance spring oscillators in watches are not ideal; oscillation of the balance is affected by friction with the air and at the balance staff pivots, and balance springs do not behave like ideal sparing. Nevertheless, the characteristics of Simple Harmonic Motion and sine waves can be used to analyse the behaviour of an balance and spring and give insights into what happens in a real watch oscillator.

To get the basic equations for simple harmonic motion, start with a vector of length (or amplitude) *A* rotating about a point *O* at \(\omega\) radians per second as shown in the figure above. The height \(x\) of the tip of this vector above a horizontal line through the origin *O* plotted against time gives the blue sine wave shown to the right.

Since there are \(2\pi\) radians in a complete circle, the time \(T\) that the vector \(A\) takes to complete one revolution is:

\[ T = \frac{2\pi}{\omega} \quad \text{ and } \quad \omega = \frac{2\pi}{T} \]T is called the period; it is also the time that the sine wave takes to make one complete oscillation, shown by *T* under the sine wave in the figure. Note that one complete oscillation always includes two peaks of amplitude, velocity and acceleration. Frequency *f* is the inverse of the period, so \( f = 1/T \) or:

A watch escapement operating at 18,000 vibrations per hour makes 5 ticks per second. The time between each tick, 0.2 seconds, represents half an oscillation, so the period, the time taken to make one complete oscillation, is 0.4 seconds. The frequency is the inverse of the period, 2.5 oscillations per second or 2.5 Hz.

### Analysis

The displacement *x*, which is the height of the tip of the amplitude vector *A* above the horizontal line passing through the origin *O*, is given by:

The vector *A* is rotating about *O* at \(\omega\) radians per second. At time *t* the angle \(\phi\) is given by \(\phi = \omega t\). The height \(x\) at time \(t\) is then:

Velocity \(v\) in the *x* direction at time *t* is given by the first differential of the equation for *x* with respect to time:

Acceleration \(a\) in the *x* direction at time *t* is given by the second differential of the equation for *x* with respect to time:

Since the second term in this equation, \(A\sin(\omega t)\), is (from the equation above) equal to *x*, this can be substituted to give:

These equations encapsulate the fundamentals of Simple Harmonic Motion. The final equation \(a=-\omega^2 x\) defines the essential condition for simple harmonic motion, which is that acceleration towards the equilibrium or rest position is proportional to displacement, that is \(a \propto x\) . Acceleration is in the opposite direction to displacement, shown by the minus sign.

## Mass Suspended on Spring

A simple system where the restoring force is proportional to displacement from the rest position is a mass suspended on a spring, as shown in the figure. If the mass is moved up or down relative to its rest position and released, it will oscillate about the rest position. To study the motion of this system, laws determined by two of the greatest British scientists are used.

Hooke's Law of springs, *ut tensio, sic vis* or “as the extension, so the force” relates the force exerted by a spring to its extension \(x\) by a spring constant \(k\). So when the mass moves up or down by a distance \(x\), the restoring force exerted on it by the spring is:

Note that the force is given a negative sign because it acts as a restoring force in the opposite direction to the movement of the mass away from its rest position. Note also that the fact that the spring is initially stretched when the mass is in its rest position, by the weight of the mass due to gravity, does not matter; this is a constant that applies all the time. It is the additional force generated in the spring when the mass moves up or down that is relevant in this analysis. If the mass and spring were horizontal, it would make no difference to the analysis.

Newton's Second Law of motion relates force to mass and acceleration:

\[F = ma\]Equating these two forces gives:

\[ma = -kx \quad \text{and therefore} \quad a = - \left ( \frac{k}{m} \right ) x \]This is of the same form as the simple harmonic motion relationship between acceleration and displacement derived above, with the constant \(\frac{k}{m}\) in the place of \(\omega^2\). Equating \(\omega^2\) with \(\frac{k}{m}\) gives:

\[ \omega^2 = \frac{k}{m} \quad \text{and therefore} \quad \omega_n=\sqrt{\frac{k}{m}}\]The frequency of oscillation is determined by the spring constant and the mass and is called the natural frequency of the system, denoted by the subscript n in \(\omega_n\). Substituting \( \sqrt{k/m} \) for \( \omega \) in the equation for period gives the natural period of oscillation:

\[T=\frac{2\pi}{\omega} = 2\pi\sqrt{\frac{m}{k}}\]## Circular Oscillations

Simple harmonic motion can also be used to examine the oscillations of a watch balance and spring.

- Mass is replaced by rotational inertia or moment of inertia \(I\).
- The displacement is angular rather than linear, denoted by \(\theta\) instead of \(x\).
- Linear force is replaced by turning force or torque \(\tau\).
- The spring constant \(S\) gives the torque in proportion to the angular displacement.
- Linear acceleration is replaced by angular acceleration \(\alpha\).

The basic equation for Simple Harmonic Motion can be restated in circular terms by replacing the acceleration \(a\) with angular acceleration \(\alpha\) and the displacement \(x\) with angular displacement \(\theta\) as follows:

\[ \alpha = -\omega^2 \theta \]Hooke's law restated for a spiral or helical spring is that the restoring torque exerted by the spring is proportional to the spring constant multiplied by the angular displacement, the angle that one end of the spring is turned relative to the other, so that:

\[ \tau = -S \theta \]Newton's second law restated for circular motion is that torque is equal to the moment of inertia of the balance multiplied by the angular acceleration, so that:

\[ \tau = I \alpha \]In an oscillating balance and spring assembly, inertia opposes torque, so Hooke's and Newton's laws can be equated as:

\[ \tau = I \alpha = -S \theta \] \[ \therefore \alpha = -\frac{S}{I} \theta \]This is the same form as the basic equation of Simple Harmonic Motion, with the constant \(\frac{S}{I}\) taking the place of \(\omega^2.\) Therefore:

\[ \omega^2 = \frac{S}{I} \quad \text{and} \quad \omega_n = \sqrt{\frac{S}{I}} \]Substituting for \(\omega\) in the equation for period gives:

\[ T = 2 \pi \sqrt{\frac{I}{S}} \]This equation does not contain any term that is dependent on the amplitude of oscillation, which means that the period \(T\) is independent of amplitude and therefore the oscillations are isochronous. A real spring balance oscillator in a watch suffers from additional factors such as air resistance, friction and the actions of unlocking and impulse which detract from perfection, but if these disturbances can be reduced, the isochronism of the oscillator and timekeeping of the watch will be improved.

## Balance Velocity

It is useful to know the velocity of the balance as it swings through the centre line of the escapement, because this is used by timing machines to estimate the amplitude.

The analysis can potentially get confusing; it must be remembered that \(\omega\) is the angular velocity of the SHM vector, which is a constant, whereas \(v\) is the angular velocity of the balance as it is oscillating, which varies from zero at the extremes of its swing, or arc of vibration, to a maximum as it passes through the middle of its swing. From the equations in the analysis section above, the angular velocity of the balance is given by:

\[ v = \omega A \cos(\omega t) \]The value of a cosine function varies between 0 and 1, so velocity is at a maximum when \(cos( \omega t) = 1 \).

\[v_{max} = \omega A \]The equation for period can be used to express \( \omega\ \) in terms of the period \(T\):

\[T = \frac{2\pi}{\omega} \quad \textrm{therefore} \quad \omega = \frac{2\pi}{T} \]Substituting for \(\omega\) into the equation for \(V_{max}\) gives:

\[ v_{max} = \frac{2\pi A}{T} \]### Velocity Example

In watch terms, an oscillation is two vibrations of the balance, or two ticks of a lever escapement. An 18,000 vph (vibrations per hour) lever escapement movement makes 5 ticks per second. Its frequency is half that; 2.5 Hz, and its period \(T=0.4\textrm{sec} \).

For an 18,000 vph balance with an amplitude of 270°:

$$\begin{aligned} v_{max} &= \frac{2\pi\times 270}{0.4} \\ &= 4,241^\circ\textrm{/sec} \\ &= 74.022 \textrm{ radians/sec} \\ &= 11.8 \textrm{ revolutions/sec} \end{aligned}$$Of course the balance isn't actually *revolving* at nearly 12 revolutions per second; this is its instantaneous peak angular velocity as it passes though the central position, which occurs twice per oscillation.

The linear velocity of the rim of the balance is given by angular velocity in radians per second multiplied by the radius of the rim. A 10mm diameter balance will therefore have a peak velocity:

\[ v_{linmax} = 74.022 \times 5.10^{-3} = 0.37 \textrm{ m/sec} \]A Hamilton Model 21 marine chronometer has a balance 35.6mm diameter. The balance oscillates at 14,400 vph to give half second ticks, so the period \( T \) = 0.5 seconds. If the amplitude is 270° or 4.7 radians:

\[v_{max} = \frac{2\pi\times 4.7}{0.5} = 59 \textrm{ radians/sec} \]## Balance Amplitude

The amplitude of the balance is a good indication of the condition of a movement. When the mainspring is fully wound, the amplitude of a clean and well-oiled movement should be at least 270° in the horizontal positions, dial up and down, and 220° when vertical.

It is difficult to measure the balance amplitude directly, so watch timing machines infer it from the timing of the sounds made during escape and impulse.

During escape and impulse, three distinct sounds are recorded by the timing machine. The first is the sound of the impulse pin contacting the lever fork, which is the start of unlocking. Once the escape wheel is unlocked, it takes over driving the lever and the opposite side of the lever fork contacts the impulse pin, which is the start of impulse and causes the second sound. The third and final sound is a combination of a pallet locking the escape wheel and the lever hitting the banking.

Strictly speaking, the lift angle is the angle through which the balance travels during impulse. However, in the context of timing machines, it is conventional to call the lift angle the angle that the balance travels from the initial sound of the impulse pin contacting the lever fork to sound of the lever contacting the banking. The lift angle can be determined physically from the geometry of the escapement. By timing how long it takes the balance to travel through the lift angle, the average velocity of the balance during the period can be calculated, from which the amplitude can be estimated.

At the point of escape and impulse, the balance is travelling at its greatest angular velocity, defined above as \( v_{max} \), which is given by the expression:

\[ v_{max} = \frac{2\pi A}{T} \]This can be rearranged in terms of amplitude:

\[ A = v_{max} \frac{T}{2 \pi} \]The lift time is the length of time it takes the balance to travel through the lift angle. The average angular velocity over the lift angle can be calculated as:

\[ v_{avg} = \frac{\text{lift angle}}{\text{lift time}} \]Although this is not quite the same as \(v_{max}\), which occurs theoretically only for an instant at the centre line position, and is never achieved in a real balance spring oscillator because of the interference of unlocking and impulse, it can be taken as a reasonable approximation. The amplitude can then be estimated as:

\[ A = \frac{\text{lift angle}}{\text{lift time}} \frac{T}{2 \pi} \]### Amplitude Example

An 18,000 vph movement has a lift angle of 52° and the lift time measured on a timing machine of 12 milliseconds (0.012 seconds). The movement makes 5 ticks per second so the ticks are \(\frac{1}{5}\) or 0.2 seconds apart. The period \(T\) of a complete oscillation is two ticks, which is 0.4 seconds. The amplitude is then estimated as:

\[ A = \frac{52}{0.12} \frac{0.4}{2 \pi} = 276^\circ \]## Energy

As a balance oscillates, its kinetic energy is converted into potential energy stored in the spring, and vice versa.

At each end of the swing, when the balance is momentarily stationary, its kinetic energy is zero and all the energy is stored in the spring. At the moment the balance passes through the mid point of its swing, it is travelling with maximum velocity and all the energy is in the kinetic energy of the balance.

The total energy in the system at any moment is the sum of the kinetic energy of the balance and the potential energy of the spring.

Kinetic energy of a mass moving in a straight line is given by the equation \( E_k = \frac{1}{2}mv^2 \). The kinetic energy of a rotating object with moment of inertia \(I\) is given by:

\[ E_k = \frac{1}{2}Iv^2 \]The kinetic energy of the balance is greatest when it is travelling with maximum velocity. The instantaneous velocity of the balance is given by:

\[ v=\omega A \cos(\omega t) \]The value of a cosine function varies between 0 and 1, so velocity is at a maximum when \(cos( \omega t) = 1 \).

\[V_{max} = \omega A \]The vector angle \(\omega\) is related to the period \(T\) by the equation:

\[ \omega = \frac{2\pi}{T} \]Substituting these into the first equation gives:

\[ E_{k(max)} = \frac{1}{2}I \left ( \frac{2\pi}{T} A \right ) ^2 \]This can also be written:

\[ E_{k(max)} = I \cdot \left ( \frac{2\pi}{T} \right ) ^2 \cdot \frac{A^2}{2} \]### Energy Example

A Hamilton Model 21 marine chronometer has a balance 35.6mm diameter with a mass of 5.1 gram. Not all of the mass will be in the rim, and the radius of gyration will be less that the radius of the balance, but for a first-order approximation let's take these dimensions. The amplitude is 270°, which is 4.7 radians. The balance oscillates at 14,400 vph, so T = 0.5 seconds.

\[ E_{k(max)} = 5.1\times10^{-3} \cdot (17.8 \times 10^{-3})^2 \cdot \left ( \frac{2\pi}{0.5} \right ) ^2 \cdot \frac{4.7^2}{2} = 2.82 \times 10^{-3} \textrm{joule} \]## Total Energy

As the balance swings the balance spring bends along its length, coiling and uncoiling, which stores and releases energy. When the balance reaches the limits of its swing at, say, +/- 270°, the spring is at its maximum bend and the balance stops moving as its swing reverses. At the this point, all the energy of the system is stored in the spring. The spring accelerates the balance towards the neutral position. As it passes through the neutral position, the velocity of the balance is at its maximum and the spring force becomes zero as it stops accelerating the balance and starts slowing it down.

Kinetic energy is given in angular terms by \( \frac{1}{2}I(v)^2 \). The Hamilton balance has a mass of 5.1 gram and a diameter of 35.6mm giving a radius of 17.8mm. Not all of the mass will be in the rim, and the radius of gyration will be less that the radius of the balance, but for a first-order approximation let's take these dimensions:

\[ \frac{1}{2} \times 5.1 \! \times \! 10^{-3} \times (17.8 \! \times \! 10^{-3})^2 \times 59^2 = 2.81\!\times\!10^{-3} joule \]The maximum **linear velocity** \(L\) of the rim is given by its radius multiplied by the maximum angular velocity:

So although the balance looks like it is moving pretty quickly, its small size means that its speed through the air is quite low. The peak kinetic energy of the balance is then:

\[ KE = \frac{1}{2} mL^2 = \frac{1}{2} 5.1\! \times\! 10^{-3} \times 1.05^2 = 2.81\!\times\!10^{-3} joule\]If you have any comments or questions, please don't hesitate to get in touch via my Contact Me page.

Copyright © David Boettcher 2005 - 2024 all rights reserved. This page updated November 2023. W3CMVS. Back to the top of the page.