Tautochrone curve

 ( Plane Curves ) 

A tautochrone curve or isochrone curve () is the curve for which the time taken by an object sliding without friction in uniform gravity to its lowest point is independent of its starting point on the curve. The curve is a cycloid, and the time is equal to π times the square root of the radius of the circle which generates the cycloid, over the acceleration of gravity. The tautochrone curve is related to the brachistochrone curve, which is also a cycloid.

---

The tautochrone problem The tautochrone problem, the attempt to identify this curve, was solved by Christiaan Huygens in 1659. He proved geometrically in his Horologium Oscillatorium, originally published in 1673, that the curve is a cycloid.

The cycloid is given by a point on a circle of radius $r$ tracing a curve as the circle rolls along the $x$ axis, as: $$\backslash begin\{align\}$$

 x &= r(\theta - \sin \theta) \\
 y &= r(1 - \cos \theta)
\end{align}
     

Note that a tautochrone as typically drawn, where objects slide downwards towards the origin, is instead traced by the cycloid of a circle rolling upside-down along the line $y\; =\; 2r$, such that: $$\backslash begin\{align\}$$

 x &= r(\theta + \sin \theta) \\
 y &= r(1 - \cos \theta)
     

\end{align}

Huygens also proved that the time of descent is equal to the time a body takes to fall vertically the same distance as diameter of the circle that generates the cycloid, multiplied by $\backslash pi\; /\; 2$. In modern terms, this means that the time of descent is $\backslash pi\; \backslash sqrt\{r/g\}$, where $r$ is the radius of the circle which generates the cycloid, and $g$ is the gravity of Earth, or more accurately, the earth's gravitational acceleration.

This solution was later used to solve the problem of the brachistochrone curve. Johann Bernoulli solved the problem in a paper ( Acta Eruditorum, 1697).

The tautochrone problem was studied by Huygens more closely when it was realized that a pendulum, which follows a circular path, was not isochronous and thus his pendulum clock would keep different time depending on how far the pendulum swung. After determining the correct path, Christiaan Huygens attempted to create pendulum clocks that used a string to suspend the bob and curb cheeks near the top of the string to change the path to the tautochrone curve. These attempts proved unhelpful for a number of reasons. First, the bending of the string causes friction, changing the timing. Second, there were much more significant sources of timing errors that overwhelmed any theoretical improvements that traveling on the tautochrone curve helps. Finally, the "circular error" of a pendulum decreases as length of the swing decreases, so better clock could greatly reduce this source of inaccuracy.

Later, the mathematicians Joseph Louis Lagrange and Leonhard Euler provided an analytical solution to the problem.

---

Lagrangian solution For a simple harmonic oscillator released from rest, regardless of its initial displacement, the time it takes to reach the lowest potential energy point is always a quarter of its period, which is independent of its amplitude. Therefore, the Lagrangian of a simple harmonic oscillator is isochronous.

In the tautochrone problem, if the particle's position is parametrized by the arclength from the lowest point, the kinetic energy is then proportional to $\backslash dot\{s\}^2$, and the potential energy is proportional to the height . One way the curve in the tautochrone problem can be an isochrone is if the Lagrangian is mathematically equivalent to a simple harmonic oscillator; that is, the height of the curve must be proportional to the arclength squared:

where the constant of proportionality is $1/(8r)$. Compared to the simple harmonic oscillator's Lagrangian, the equivalent spring constant is $k=mg/(4r)$, and the time of descent is $T/4=\backslash frac\{\backslash pi\}\{2\}\; \backslash sqrt\{\backslash frac\{m\}\{k\}\}=\backslash pi\; \backslash sqrt\{\backslash frac\{r\}\{g\}\}.$ However, the physical meaning of the constant $r$ is not clear until we determine the exact analytical equation of the curve.

To solve for the analytical equation of the curve, note that the differential form of the above relation is

which eliminates , and leaves a differential equation for and . This is the differential equation for a cycloid when the vertical coordinate is counted from its vertex (the point with a horizontal tangent) instead of the cusp.

To find the solution, integrate for in terms of :

{\sqrt{h}}, \\ x &= -4r\int \sqrt{1-u^2} \, du,

where $u\; =\; \backslash sqrt\{h/(2r)\}$, and the height decreases as the particle moves forward . This integral is the area under a circle, which can be done with another substitution $u=\backslash cos\; (t/2)$ and yield:

This is the standard parameterization of a cycloid with $h=2r-y$. It's interesting to note that the arc length squared is equal to the height difference multiplied by the full arch length $8r$.

---

"Virtual gravity" solution The simplest solution to the tautochrone problem is to note a direct relation between the angle of an incline and the gravity felt by a particle on the incline. A particle on a 90° vertical incline undergoes full gravitational acceleration $g$, while a particle on a horizontal plane undergoes zero gravitational acceleration. At intermediate angles, the acceleration due to "virtual gravity" by the particle is $g\backslash sin\backslash theta$. Note that $\backslash theta$ is measured between the tangent to the curve and the horizontal, with angles above the horizontal being treated as positive angles. Thus, $\backslash theta$ varies from $-\backslash pi/2$ to $\backslash pi/2$.

The position of a mass measured along a tautochrone curve, $s(t)$, must obey the following differential equation:

(Based loosely on Proctor, pp. 135–139)

---

Abel's solution Niels Henrik Abel attacked a generalized version of the tautochrone problem ( Abel's mechanical problem), namely, given a function $T(y)$ that specifies the total time of descent for a given starting height, find an equation of the curve that yields this result. The tautochrone problem is a special case of Abel's mechanical problem when $T(y)$ is a constant.

Abel's solution begins with the principle of conservation of energy – since the particle is frictionless, and thus loses no energy to heat, its kinetic energy at any point is exactly equal to the difference in gravitational potential energy from its starting point. The kinetic energy is $\backslash frac\{1\}\{2\}\; mv^2$, and since the particle is constrained to move along a curve, its velocity is simply $\{d\backslash ell\}/\{dt\}$, where $\backslash ell$ is the distance measured along the curve. Likewise, the gravitational potential energy gained in falling from an initial height $y\_0$ to a height $y$ is $mg(y\_0\; -\; y)$, thus:

\\
     

dt & = - \frac{1}{\sqrt{2g(y_0-y)}} \frac{d\ell}{dy} \,dy \end{align}

In the last equation, we have anticipated writing the distance remaining along the curve as a function of height ($\backslash ell(y))$, recognized that the distance remaining must decrease as time increases (thus the minus sign), and used the chain rule in the form $d\backslash ell\; =\; \backslash frac\{d\backslash ell\}\{dy\}\; dy$.

Now we integrate from $y\; =\; y\_0$ to $y\; =\; 0$ to get the total time required for the particle to fall:

\int_0^{y_0} \frac{1}{\sqrt{y_0-y}} \frac{d\ell}{dy} \, dy
     

This is called Abel's integral equation and allows us to compute the total time required for a particle to fall along a given curve (for which $\{d\backslash ell\}/\{dy\}$ would be easy to calculate). But Abel's mechanical problem requires the converse – given $T(y\_0)\backslash ,$, we wish to find $f(y)\; =\; \{d\backslash ell\}/\{dy\}$, from which an equation for the curve would follow in a straightforward manner. To proceed, we note that the integral on the right is the convolution of $\{d\backslash ell\}/\{dy\}$ with $\{1\}/\{\backslash sqrt\{y\}\}$ and thus take the Laplace transform of both sides with respect to variable $y$:

\mathcal{L} \left [ \frac{1}{\sqrt{y}} \right ]F(s)
     

where $F(s)\; =\; \backslash mathcal\{L\}\; \{\backslash left\}$. Since $\backslash mathcal\{L\}\; \{\backslash left\}\; =\; \backslash sqrt$, we now have an expression for the Laplace transform of $\{d\backslash ell\}/\{dy\}$ in terms of the Laplace transform of $T(y\_0)$:

s^{\frac{1}{2}} \mathcal{L}[T(y_0)]
     

This is as far as we can go without specifying $T(y\_0)$. Once $T(y\_0)$ is known, we can compute its Laplace transform, calculate the Laplace transform of $\{d\backslash ell\}/\{dy\}$ and then take the inverse transform (or try to) to find $\{d\backslash ell\}/\{dy\}$.

For the tautochrone problem, $T(y\_0)\; =\; T\_0\backslash ,$ is constant. Since the Laplace transform of 1 is $\{1\}/\{s\}$, i.e., $\backslash mathcal\{L\}T(y\_0)\; =\; \{T\_0\}/\{s\}$, we find the shape function $f(y)\; =\; \{d\backslash ell\}/\{dy\}$:

s^{\frac{1}{2}} \mathcal{L}[T_0] \\
     

\end{align}

Making use again of the Laplace transform above, we invert the transform and conclude:

( Simmons, Section 54).

Using the fact that $\backslash left(\backslash frac\{d\backslash ell\}\{dy\}\backslash right)^2\; =\; \backslash left(\backslash frac\{dx\}\{dy\}\backslash right)^2\; +\; 1$, we can rewrite the above as:

Let:

So that:

\,dy
     

\end{align}

In order to integrate, substitute:

We can now rewrite the equation for $dx$ as:

{2r \sin^2 \frac{\theta}{2}}} 2r \sin \frac{\theta}{2} \cos \frac{\theta}{2} \,d \theta \\ & = 2r \sqrt{\frac{1 - \sin^2 \frac{\theta}{2}}{\sin^2 \frac{\theta}{2}}} \sin \frac{\theta}{2} \cos \frac{\theta}{2} \,d \theta \\ & = 2r \sqrt{\frac{\cos^2 \frac{\theta}{2}}{\sin^2 \frac{\theta}{2}}} \sin \frac{\theta}{2} \cos \frac{\theta}{2} \,d \theta \\ & = 2r \cos^2 \frac{\theta}{2} \,d \theta \\ & = r(\cos \theta + 1) \,d \theta \end{align}

Integrating both sides arrives at:

And as we parameterized $y$ as:

we now have equations for both $x$ and $y$ parameterized by $\backslash theta$:

which trace a tautochrone sitting on top of the $x$ axis.

(Based on O'neil, pp. 263-264)

---

See also

• Beltrami identity
• Brachistochrone curve
• Calculus of variations
• Catenary
• Cycloid
• Uniformly accelerated motion

---

Bibliography

• George Simmons (1972). 9780070575400, McGraw–Hill. ISBN 9780070575400
• (1983). 9780534011369, Wadsworth Pub. Co. ISBN 9780534011369

---

External links

• Mathworld

Post Comment