Movatterモバイル変換


[0]ホーム

URL:


Jump to content
WikipediaThe Free Encyclopedia
Search

Tautochrone curve

From Wikipedia, the free encyclopedia
(Redirected fromTautochrone)
Curve for which the time to roll to the end is equal for all starting points
Four balls slide down a cycloid curve from different positions, but they arrive at the bottom at the same time. The blue arrows show the points' acceleration along the curve. On the top is the time-position diagram.
Objects representing tautochrone curve

Atautochrone curve orisochrone curve (from Ancient Greek ταὐτό (tauto-) 'same' ἴσος (isos-) 'equal' and χρόνος (chronos) 'time') is thecurve for which the time taken by an object sliding withoutfriction in uniformgravity to its lowest point is independent of its starting point on the curve. The curve is acycloid, and the time is equal toπ times thesquare root of the radius of the circle which generates the cycloid, over theacceleration of gravity. The tautochrone curve is related to thebrachistochrone curve, which is also a cycloid.

The tautochrone problem

[edit]
Christiaan Huygens,Horologium oscillatorium sive de motu pendulorum, 1673

It was in the left hand try-pot of the Pequod, with the soapstone diligently circling round me, that I was first indirectly struck by the remarkable fact, that in geometry all bodies gliding along the cycloid, my soapstone for example, will descend from any point in precisely the same time.

Moby Dick byHerman Melville, 1851

The tautochrone problem, the attempt to identify this curve, was solved byChristiaan Huygens in 1659. He proved geometrically in hisHorologium Oscillatorium, originally published in 1673, that the curve is acycloid.

On a cycloid whose axis is erected on the perpendicular and whose vertex is located at the bottom, the times of descent, in which a body arrives at the lowest point at the vertex after having departed from any point on the cycloid, are equal to each other ...[1]

The cycloid is given by a point on a circle of radiusr{\displaystyle r} tracing a curve as the circle rolls along thex{\displaystyle x} axis, as:x=r(θsinθ)y=r(1cosθ),{\displaystyle {\begin{aligned}x&=r(\theta -\sin \theta )\\y&=r(1-\cos \theta ),\end{aligned}}}

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π/2{\displaystyle \pi /2}. In modern terms, this means that the time of descent isπr/g{\textstyle \pi {\sqrt {r/g}}}, wherer{\displaystyle r} is the radius of the circle which generates the cycloid, andg{\displaystyle g} is thegravity of Earth, or more accurately, the earth's gravitational acceleration.

Five isochronous cycloidal pendulums with different amplitudes

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

Schematic of acycloidal pendulum

The tautochrone problem was studied by Huygens more closely when it was realized that a pendulum, which follows a circular path, was notisochronous and thus hispendulum 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 clockescapements could greatly reduce this source of inaccuracy.

Later, the mathematiciansJoseph Louis Lagrange andLeonhard Euler provided an analytical solution to the problem.

Lagrangian solution

[edit]

For asimple 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 isisochronous.

In the tautochrone problem, if the particle's position is parametrized by thearclengths(t) from the lowest point, the kinetic energy is then proportional tos˙2{\displaystyle {\dot {s}}^{2}}, and the potential energy is proportional to the heighth(s). 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:

h(s)=s2/(8r),{\displaystyle h(s)=s^{2}/(8r),}

where the constant of proportionality is1/(8r){\displaystyle 1/(8r)}. Compared to the simple harmonic oscillator'sLagrangian, the equivalent spring constant isk=mg/(4r){\displaystyle k=mg/(4r)}, and the time of descent isT/4=π2mk=πrg.{\displaystyle T/4={\frac {\pi }{2}}{\sqrt {\frac {m}{k}}}=\pi {\sqrt {\frac {r}{g}}}.} However, the physical meaning of the constantr{\displaystyle 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

dh=sds/(4r),dh2=s2ds2/(16r2)=h(dx2+dh2)/(2r),(dxdh)2=2rh1{\displaystyle {\begin{aligned}dh&=s\,ds/(4r),\\dh^{2}&=s^{2}\,ds^{2}/(16r^{2})=h\left(dx^{2}+dh^{2}\right)/(2r),\\\left({\frac {dx}{dh}}\right)^{2}&={\frac {2r}{h}}-1\end{aligned}}}

which eliminatess, and leaves a differential equation fordx anddh. This is the differential equation for acycloid when the vertical coordinateh is counted from its vertex (the point with a horizontal tangent) instead of thecusp.

To find the solution, integrate forx in terms ofh:

dxdh=2rhh,x=4r1u2du,{\displaystyle {\begin{aligned}{\frac {dx}{dh}}&=-{\frac {\sqrt {2r-h}}{\sqrt {h}}},\\x&=-4r\int {\sqrt {1-u^{2}}}\,du,\end{aligned}}}

whereu=h/(2r){\displaystyle u={\sqrt {h/(2r)}}}, and the height decreases as the particle moves forwarddx/dh<0{\displaystyle dx/dh<0}. This integral is the area under a circle, which can be done with another substitutionu=cos(t/2){\displaystyle u=\cos(t/2)} and yield:

x=r(tsint),h=r(1+cost).{\displaystyle {\begin{aligned}x&=r(t-\sin t),\\h&=r(1+\cos t).\end{aligned}}}

This is the standard parameterization of acycloid withh=2ry{\displaystyle h=2r-y}. It's interesting to note that the arc length squared is equal to the height difference multiplied by the full arch length8r{\displaystyle 8r}.

"Virtual gravity" solution

[edit]

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 accelerationg{\displaystyle g}, while a particle on a horizontal plane undergoes zero gravitational acceleration. At intermediate angles, the acceleration due to "virtual gravity" by the particle isgsinθ{\displaystyle g\sin \theta }. Note thatθ{\displaystyle \theta } is measured between the tangent to the curve and the horizontal, with angles above the horizontal being treated as positive angles. Thus,θ{\displaystyle \theta } varies fromπ/2{\displaystyle -\pi /2} toπ/2{\displaystyle \pi /2}.

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

d2sdt2=ω2s{\displaystyle {\frac {d^{2}s}{{dt}^{2}}}=-\omega ^{2}s}

which, along with the initial conditionss(0)=s0{\displaystyle s(0)=s_{0}} ands(0)=0{\displaystyle s'(0)=0}, has solution:

s(t)=s0cosωt{\displaystyle s(t)=s_{0}\cos \omega t}

It can be easily verified both that this solution solves the differential equation and that a particle will reachs=0{\displaystyle s=0} at timeπ/2ω{\displaystyle \pi /2\omega } from any starting positions0{\displaystyle s_{0}}. The problem is now to construct a curve that will cause the mass to obey the above motion.Newton's second law shows that the force of gravity and the acceleration of the mass are related by:

gsinθ=d2sdt2=ω2s{\displaystyle {\begin{aligned}-g\sin \theta &={\frac {d^{2}s}{{dt}^{2}}}\\&=-\omega ^{2}s\,\end{aligned}}}

The explicit appearance of the distance,s{\displaystyle s}, is troublesome, but we candifferentiate to obtain a more manageable form:

gcosθdθ=ω2dsds=gω2cosθdθ{\displaystyle {\begin{aligned}g\cos \theta \,d\theta &=\omega ^{2}\,ds\\\Longrightarrow ds&={\frac {g}{\omega ^{2}}}\cos \theta \,d\theta \end{aligned}}}

This equation relates the change in the curve's angle to the change in the distance along the curve. We now usetrigonometry to relate the angleθ{\displaystyle \theta } to the differential lengthsdx{\displaystyle dx},dy{\displaystyle dy} andds{\displaystyle ds}:

ds=dxcosθds=dysinθ{\displaystyle {\begin{aligned}ds={\frac {dx}{\cos \theta }}\\ds={\frac {dy}{\sin \theta }}\end{aligned}}}

Replacingds{\displaystyle ds} withdx{\displaystyle dx} in the above equation lets us solve forx{\displaystyle x} in terms ofθ{\displaystyle \theta }:

ds=gω2cosθdθdxcosθ=gω2cosθdθdx=gω2cos2θdθ=g2ω2(cos2θ+1)dθx=g4ω2(sin2θ+2θ)+Cx{\displaystyle {\begin{aligned}ds&={\frac {g}{\omega ^{2}}}\cos \theta \,d\theta \\{\frac {dx}{\cos \theta }}&={\frac {g}{\omega ^{2}}}\cos \theta \,d\theta \\dx&={\frac {g}{\omega ^{2}}}\cos ^{2}\theta \,d\theta \\&={\frac {g}{2\omega ^{2}}}\left(\cos 2\theta +1\right)d\theta \\x&={\frac {g}{4\omega ^{2}}}\left(\sin 2\theta +2\theta \right)+C_{x}\end{aligned}}}

Likewise, we can also expressds{\displaystyle ds} in terms ofdy{\displaystyle dy} and solve fory{\displaystyle y} in terms ofθ{\displaystyle \theta }:

ds=gω2cosθdθdysinθ=gω2cosθdθdy=gω2sinθcosθdθ=g2ω2sin2θdθy=g4ω2cos2θ+Cy{\displaystyle {\begin{aligned}ds&={\frac {g}{\omega ^{2}}}\cos \theta \,d\theta \\{\frac {dy}{\sin \theta }}&={\frac {g}{\omega ^{2}}}\cos \theta \,d\theta \\dy&={\frac {g}{\omega ^{2}}}\sin \theta \cos \theta \,d\theta \\&={\frac {g}{2\omega ^{2}}}\sin 2\theta \,d\theta \\y&=-{\frac {g}{4\omega ^{2}}}\cos 2\theta +C_{y}\end{aligned}}}

Substitutingϕ=2θ{\displaystyle \phi =2\theta } andr=g4ω2{\textstyle r={\frac {g}{4\omega ^{2}}}\,}, we see that theseparametric equations forx{\displaystyle x} andy{\displaystyle y} are those of a point on a circle of radiusr{\displaystyle r} rolling along a horizontal line (acycloid), with the circle center at the coordinates(Cx+rϕ,Cy){\displaystyle (C_{x}+r\phi ,C_{y})}:

x=r(sinϕ+ϕ)+Cxy=rcosϕ+Cy{\displaystyle {\begin{aligned}x&=r\left(\sin \phi +\phi \right)+C_{x}\\y&=-r\cos \phi +C_{y}\end{aligned}}}

Note thatϕ{\displaystyle \phi } ranges fromπϕπ{\displaystyle -\pi \leq \phi \leq \pi }. It is typical to setCx=0{\displaystyle C_{x}=0} andCy=r{\displaystyle C_{y}=r} so that the lowest point on the curve coincides with the origin. Therefore:

x=r(ϕ+sinϕ)y=r(1cosϕ){\displaystyle {\begin{aligned}x&=r\left(\phi +\sin \phi \right)\\y&=r\left(1-\cos \phi \right)\\\end{aligned}}}

Solving forω{\displaystyle \omega } and remembering thatT=π2ω{\displaystyle T={\frac {\pi }{2\omega }}} is the time required for descent, being a quarter of a whole cycle, we find the descent time in terms of the radiusr{\displaystyle r}:

r=g4ω2ω=12grT=πrg{\displaystyle {\begin{aligned}r&={\frac {g}{4\omega ^{2}}}\\\omega &={\frac {1}{2}}{\sqrt {\frac {g}{r}}}\\T&=\pi {\sqrt {\frac {r}{g}}}\\\end{aligned}}}

(Based loosely onProctor, pp. 135–139)

Abel's solution

[edit]

Niels Henrik Abel attacked a generalized version of the tautochrone problem (Abel's mechanical problem), namely, given a functionT(y){\displaystyle 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 whenT(y){\displaystyle T(y)} is a constant.

Abel's solution begins with the principle ofconservation of energy – since the particle is frictionless, and thus loses no energy toheat, itskinetic energy at any point is exactly equal to the difference ingravitational potential energy from its starting point. The kinetic energy is12mv2{\textstyle {\frac {1}{2}}mv^{2}}, and since the particle is constrained to move along a curve, its velocity is simplyd/dt{\displaystyle {d\ell }/{dt}}, where{\displaystyle \ell } is the distance measured along the curve. Likewise, the gravitational potential energy gained in falling from an initial heighty0{\displaystyle y_{0}} to a heighty{\displaystyle y} ismg(y0y){\displaystyle mg(y_{0}-y)}, thus:

12m(ddt)2=mg(y0y)ddt=±2g(y0y)dt=±d2g(y0y)dt=12g(y0y)ddydy{\displaystyle {\begin{aligned}{\frac {1}{2}}m\left({\frac {d\ell }{dt}}\right)^{2}&=mg(y_{0}-y)\\{\frac {d\ell }{dt}}&=\pm {\sqrt {2g(y_{0}-y)}}\\dt&=\pm {\frac {d\ell }{\sqrt {2g(y_{0}-y)}}}\\dt&=-{\frac {1}{\sqrt {2g(y_{0}-y)}}}{\frac {d\ell }{dy}}\,dy\end{aligned}}}

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

Now we integrate fromy=y0{\displaystyle y=y_{0}} toy=0{\displaystyle y=0} to get the total time required for the particle to fall:

T(y0)=y=y0y=0dt=12g0y01y0yddydy{\displaystyle T(y_{0})=\int _{y=y_{0}}^{y=0}\,dt={\frac {1}{\sqrt {2g}}}\int _{0}^{y_{0}}{\frac {1}{\sqrt {y_{0}-y}}}{\frac {d\ell }{dy}}\,dy}

This is calledAbel's integral equation and allows us to compute the total time required for a particle to fall along a given curve (for whichd/dy{\displaystyle {d\ell }/{dy}} would be easy to calculate). But Abel's mechanical problem requires the converse – givenT(y0){\displaystyle T(y_{0})\,}, we wish to findf(y)=d/dy{\displaystyle f(y)={d\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 theconvolution ofd/dy{\displaystyle {d\ell }/{dy}} with1/y{\displaystyle {1}/{\sqrt {y}}} and thus take theLaplace transform of both sides with respect to variabley{\displaystyle y}:

L[T(y0)]=12gL[1y]F(s){\displaystyle {\mathcal {L}}[T(y_{0})]={\frac {1}{\sqrt {2g}}}{\mathcal {L}}\left[{\frac {1}{\sqrt {y}}}\right]F(s)}

whereF(s)=L[d/dy]{\displaystyle F(s)={\mathcal {L}}{\left[{d\ell }/{dy}\right]}}. SinceL[1/y]=π/s{\textstyle {\mathcal {L}}{\left[{1}/{\sqrt {y}}\right]}={\sqrt {{\pi }/{s}}}}, we now have an expression for the Laplace transform ofd/dy{\displaystyle {d\ell }/{dy}} in terms of the Laplace transform ofT(y0){\displaystyle T(y_{0})}:

L[ddy]=2gπs12L[T(y0)]{\displaystyle {\mathcal {L}}\left[{\frac {d\ell }{dy}}\right]={\sqrt {\frac {2g}{\pi }}}s^{\frac {1}{2}}{\mathcal {L}}[T(y_{0})]}

This is as far as we can go without specifyingT(y0){\displaystyle T(y_{0})}. OnceT(y0){\displaystyle T(y_{0})} is known, we can compute its Laplace transform, calculate the Laplace transform ofd/dy{\displaystyle {d\ell }/{dy}} and then take the inverse transform (or try to) to findd/dy{\displaystyle {d\ell }/{dy}}.

For the tautochrone problem,T(y0)=T0{\displaystyle T(y_{0})=T_{0}\,} is constant. Since the Laplace transform of 1 is1/s{\displaystyle {1}/{s}}, i.e.,L[T(y0)]=T0/s{\textstyle {\mathcal {L}}[T(y_{0})]={T_{0}}/{s}}, we find the shape functionf(y)=d/dy{\textstyle f(y)={d\ell }/{dy}}:

F(s)=L[ddy]=2gπs12L[T0]=2gπT0s12{\displaystyle {\begin{aligned}F(s)={\mathcal {L}}{\left[{\frac {d\ell }{dy}}\right]}&={\sqrt {\frac {2g}{\pi }}}s^{\frac {1}{2}}{\mathcal {L}}[T_{0}]\\&={\sqrt {\frac {2g}{\pi }}}T_{0}s^{-{\frac {1}{2}}}\end{aligned}}}

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

ddy=T02gπ1y{\displaystyle {\frac {d\ell }{dy}}=T_{0}{\frac {\sqrt {2g}}{\pi }}{\frac {1}{\sqrt {y}}}}

It can be shown that the cycloid obeys this equation. It needs one step further to do the integral with respect toy{\displaystyle y} to obtain the expression of the path shape.

(Simmons, Section 54).

See also

[edit]

References

[edit]
  1. ^Blackwell, Richard J. (1986).Christiaan Huygens' The Pendulum Clock. Ames, Iowa: Iowa State University Press. Part II, Proposition XXV, p. 69.ISBN 0-8138-0933-9.

Bibliography

[edit]

External links

[edit]
Retrieved from "https://en.wikipedia.org/w/index.php?title=Tautochrone_curve&oldid=1276375923"
Categories:
Hidden categories:

[8]ページ先頭

©2009-2025 Movatter.jp