wu :: forums (http://www.ocf.berkeley.edu/~wwu/cgi-bin/yabb/YaBB.cgi)
riddles >> general problem-solving / chatting / whatever >> Mechanical Pendulum Gear Clocks, Pendulum
(Message started by: rloginunix on Jul 17th, 2015, 4:38pm)

Title: Mechanical Pendulum Gear Clocks, Pendulum
Post by rloginunix on Jul 17th, 2015, 4:38pm
Mechanical Pendulum Gear Clocks (https://www.ocf.berkeley.edu/~wwu/cgi-bin/yabb/YaBB.cgi?board=riddles_general;action=display;num=1437175174;start=0#0)

Pendulum

1) Law of Motion
2) Solution
3) Small Angles
4) g, http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/pi.gif, Reality
5) Procedure
6) Sample Calculation
7) Huygens
8) More Pendulums


Law of Motion

As it is customary in physics which, as was duly noted, produces exact solutions for approximate problems, we shall at first eliminate some number of practical phenomena from the study of pendulum's motion. We shall discard the friction, the sound waves production, the reactionary forces of the pivot assembly, the mass of the pendulum's rod, its geometry change due to motion, the air resistance and a plethora of atmospheric factors like pressure, temperature, humidity, etc.

Such an arrangement of a point mass m swinging at the end of an absolutely rigid massless rod of length L about a fixed point in a uniform gravitational field is called a mathematical or ideal or simple (gravity) pendulum.

To find a practically useful information - the length of the rod L for a known period T we first need to find the pendulum's law of motion. We shall derive it in three ways.


Newton

Before we jump to the equations let us talk about pendulum's motion at the gut or intuitive level. When the pendulum is at rest it occupies the position along the equilibrium line ZO. The force due to gravity mg is balanced by the rod's tension force N, the pendulum's velocity is zero and so are its acceleration, kinetic and potential energies. Using our muscle power we move the pendulum to its left extreme position, point A, and let it go with zero initial velocity.

At this extreme point A (as well as at B) the magnitude of the vector of m's normal acceleration an which points directly at Z and tracks the change of direction of m's velocity is the smallest and so are the magnitudes of the tension force N and m's kinetic energy. The magnitude of the vector of the tangential acceleration at which is tangent to m's trajectory and tracks the change of magnitude of m's velocity is the largest and so is m's potential energy. After being zero for an instance the magnitude of m's velocity v grows quickly at first and at that stage vectors at and v point in the same direction - towards the equilibrium point O. The magnitude and the direction of the force due to gravity mg remains constant at all times and if a series of snapshots of the oscillating point mass are taken it will show that mg's vectors at all these points remain parallel and are in the same plane in which our pendulum swings (an idealization, see below).

By the time m arrives at O it covers one quarter of its period. The magnitudes of m's an, N, v and kinetic energy reach their maximums while the magnitudes of its at and potential energy are at their lowest. As m moves away from O towards B the magnitudes of an, N and v begin to shrink, at becomes negative as m begins to slow down and as a result vectors at and v point in different directions - at points towards O again while v points towards B.

By the time m arrives at B it covers one half of its period. The magnitudes of m's an, N, v and kinetic energy reach their minimums while the magnitudes of at and potential energy reach their maximums. The velocity of m changes its direction towards O and begins to grow, at points in the same direction and the above A to O scenario plays out again but this time in the opposite direction, from B to O. At O m covers three quarters of its period and by the time it reaches A again the entire period is covered:

http://www.ocf.berkeley.edu/~wwu/YaBBAttachments/rlu_pgc01.png

According to Newton's second law the vector sum of all the forces acting on a body equals the first derivative of mv over time. Let us project these forces on the http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/tau.gif axis that runs perpendicular to ZP where P is the location of the point mass m at an arbitrary time t. Since the mass remains constant we take the derivative of vt only:

d(mvt)/dt = mdvt/dt     (p.1)


where vt is m's velocity. As a vector vt is always tangential to m's trajectory and it is the first derivative of the length of m's path s over time:

vt = ds/dt     (p.2)


Since the rod is rigid m traces a circular arc. Therefore we can express s as a function of an angle http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif if the latter itself is expressed in radians. From the definition of a radian we have:

s = Lhttp://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif


Therefore for very very small portions of the path ds we have (L is const):

ds = Ldhttp://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif


Putting ds back into (p.2) and vt into (p.1) we get:

vt = ds/dt = Ldhttp://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif/dt = Lhttp://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif'

mdvt/dt = md(Lhttp://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif')/dt = mLd2http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif/dt2 = mL http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif''


where ' designates a derivative.


[e]In the drawing above make the tension force T to be N to avoid confusion with the period T[/e]

Title: Re: Mechanical Pendulum Gear Clocks, Pendulum
Post by rloginunix on Jul 17th, 2015, 4:50pm
The projection of the vector N on http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/tau.gif is zero since the two are orthogonal. To project mg on http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/tau.gif we observe that the straight lines ZO and PC are dead verticals and hence are parallel - the length of our pendulum's rod L is way smaller than the radius of Earth whose curvature we ignore in this case. This of course is an idealization since strictly speaking the verticals are the straight lines that pass through the center of the Earth, they are perpendicular to the spherical surface of the Earth and as such are not really parallel. Approximating the shape of the Earth with a perfect sphere is yet another idealization.

The straight lines ZP and CD are both perpendicular to the same straight line http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/tau.gif and hence are parallel. The straight line CD when extended intersects ZO at D' and hence by definition the quadrilateral ZPCD' is a parallelogram in which, according to Euclid (Elements, Book 1 Proposition 34), the opposite angles are equal which in turn means that http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/angle.gif OZP = http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/angle.gif PCD = http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif.

The force Ft that always points towards O is a restoring force - its direction is always opposite to that of m's angular displacement. It means that if we displace m in one direction then Ft would always act in the opposite direction as it would want to bring m to its equilibrium point O. For example, at P the angle http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif is measured clockwise relative to the vertical ZO - the equilibrium line - and hence we take it to be negative while Ft points in the opposite direction, as it wants to move m counterclockwise, and we take it to be positive. If P were on the other side of O, closer to B, then the angle http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif would have been measured counterclockwise, positive, while Ft would be pointing against this displacement and we would take it to be negative. This restoring behavior is captured with the minus sign:

Ft = -mgSin(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif)


Find the analogy here with an ideal spring obeying Hooke's law where the restoring force F = -khttp://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/cdelta.gif x is proportional to a linear measure, http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/cdelta.gif x, while in our case the restoring force is proportional to an angular measure.

Finally we have:

mLhttp://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif'' = -mgSin(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif)

http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif'' + (g/L)Sin(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif) = 0     (p.3)



Torque

Torques and angular momenta make sense only when the reference point is specified. In our case the reference point is Z and ZP is m's radius vector pointing from Z to P with a known and constant length L.

On the one hand, m's torque relative to Z is the cross product of two vectors:

http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/tau.gifZ = L http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/times.gif mg     (vector form)

http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/tau.gifZ = -mgLSin(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif)


On the other hand, it is the first derivative of m's angular momentum M relative to Z over time:

http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/tau.gifZ = dMZ/dt


Since we agreed that the size of the point mass m is much smaller than L for its angular momentum we have:

MZ = L http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/times.gif mv      (vector form)

MZ = Lmv = Lmhttp://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/omega.gif L = mL2dhttp://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif/dt = mL2http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif'

http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/tau.gifZ = d(mL2http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif')/dt = mL2http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif''


since m and L remain constant over time. Equating the two expressions we get:

-mgLSin(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif) = mL2http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif''

http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif'' + (g/L)Sin(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif) = 0     (p.3)

Title: Re: Mechanical Pendulum Gear Clocks, Pendulum
Post by rloginunix on Jul 17th, 2015, 4:57pm
Lagrange

Our simple pendulum can execute an independent motion in only one plane. The nature of this motion is such that just one variable, the angle http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif, is enough to fully describe it. Consequently our pendulum possesses only one degree of freedom and its Lagrangian in terms of http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif is:

http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/clambda.gif = T - U

T = mv2/2 = mL2http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif'2/2

U = mgL(1 - Cos(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif))

http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/clambda.gif = mL2http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif'2/2 - mgL + mgLCos(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif)


Since there is only one degree of freedom the Euler-Lagrangian equations collapse into just one:

d/dt(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/partial.gif http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/clambda.gif/http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/partial.gif http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif') = http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/partial.gif http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/clambda.gif/http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/partial.gif http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif

http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/partial.gif http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/clambda.gif/http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/partial.gif http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif' = mL2http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif'

d(mL2http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif')/dt = mL2http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif''

http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/partial.gif http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/clambda.gif/http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/partial.gif http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif= -mgLSin(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif)


Finally we have:

mL2http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif'' = -mgLSin(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif)

http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif'' + (g/L)Sin(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif) = 0     (p.3)


By solving the above equation we hope to obtain the behavior of the angle http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif as a function of time. Before we do that however let us record the remaining entities as a function of the angle http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif.

We obtain v from the law of conservation of energy. Since the tension force N is always perpendicular to v it does not do any work. The force due to gravity mg is conservative - it does not do any work either. Hence, the entire change in m's potential energy is translated into the change of its kinetic energy:

mghttp://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/cdelta.gif h = mg(h0 - h) = mg(L - LCos(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif0) - L + LCos(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif))

mgL(Cos(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif) - Cos(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif0)) = mv2/2

v(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif) = http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/surd.gif2gL(Cos(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif) - Cos(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif0))


Both accelerations follow:

at(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif) = -gSin(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif)

an(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif) = v2/L = 2g(Cos(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif) - Cos(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif0))


To find N we project it and mg on the normal axis ZP:

man = N - Fn = N - mgCos(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif)

N = man + mgCos(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif) = 2mg(Cos(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif) - Cos(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif0)) + mgCos(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif)

N(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif) = 3mgCos(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif) - 2mgCos(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif0)

Title: Re: Mechanical Pendulum Gear Clocks, Pendulum
Post by rloginunix on Jul 17th, 2015, 5:01pm
Solution

Rearrange (p.3):

http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif'' = - (g/L)Sin(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif)


To save us some typing we designate:

http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/omega.gif2 = g/L

http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif'' = -http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/omega.gif2Sin(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif)     (p.4)


Spell out the definition of a derivative http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif'' on the left hand side of (p.4) and then multiply both sides of (p.4) by dt:

http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif'' = dhttp://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif'/dt = -http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/omega.gif2Sin(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif)

dhttp://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif' = -http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/omega.gif2Sin(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif)dt     (p.5)


Multiply both sides of (p.5) by http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif' and expand its definition on the right hand side:

http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif'dhttp://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif' = -http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/omega.gif2Sin(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif)dthttp://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif' = -http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/omega.gif2Sin(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif)dt(dhttp://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif/dt) = -http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/omega.gif2Sin(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif)dhttp://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif

http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif'dhttp://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif' = -http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/omega.gif2Sin(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif)dhttp://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif    (p.6)


What we managed to do now is to completely separate the variables which means that we can integrate both sides of (p.6):

http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/int.gifhttp://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif'dhttp://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif' = http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif'2/2 + C1

http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/int.gif-http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/omega.gif2Sin(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif)dhttp://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif= http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/omega.gif2Cos(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif) + C2

http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif'2/2 + C1 = http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/omega.gif2Cos(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif) + C2

http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif'2/2 = http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/omega.gif2Cos(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif) + C     (p.7)


Use the boundary conditions to find C:

http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif(t = 0) = http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif0

http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif'(t = 0) = 0

0 = http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/omega.gif2 Cos(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif0) + C

C = -http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/omega.gif2 Cos(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif0)

Title: Re: Mechanical Pendulum Gear Clocks, Pendulum
Post by rloginunix on Jul 17th, 2015, 5:04pm
Put C into (p.7):

http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif'2/2 = http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/omega.gif2Cos(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif) - http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/omega.gif2Cos(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif0)

http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif'2 = 2http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/omega.gif2(Cos(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif) - Cos(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif0))

http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif' = http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/omega.gif(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/surd.gif2)http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/surd.gif(Cos(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif) - Cos(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif0))     (p.8)


Use the power reduction half-angle trigonometric identity:

Cos(x) = 1 - 2Sin2(x/2)


for both cosines under the square root in (p.8):

http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif' = http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/omega.gif(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/surd.gif2)http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/surd.gif(1 - 2Sin2(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif/2) - 1 + 2Sin2(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif0/2))

http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif' = dhttp://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif/dt = 2http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/omega.gifhttp://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/surd.gif(Sin2(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif0/2) - Sin2(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif/2))


Separate the variables one more time:

dt = (1/2http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/omega.gif)dhttp://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif/(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/surd.gif(Sin2(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif0/2) - Sin2(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif/2)))     (p.9)


Integrate the left hand side of (p.9) from 0 to T/4 - one quarter of the period. Integrate the right hand side of (p.9) from 0 to http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif0:

T(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif0) = (2/http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/omega.gif)http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/int.gifdhttp://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif/(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/surd.gif(Sin2(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif0/2) - Sin2(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif/2)))     (p.10)


Within the given integral Sin(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif0/2) is just a constant, designate it as:

k = Sin(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif0/2)     (p.11)


And use the following substitution:

Sin(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif/2) = kSin(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/phi.gif)     (p.12)


Differentiate (p.12):

(dhttp://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif/2)Cos(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif/2) = kCos(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/phi.gif)dhttp://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/phi.gif

dhttp://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif= 2kCos(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/phi.gif)dhttp://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/phi.gif/Cos(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif/2)     (p.13)

Title: Re: Mechanical Pendulum Gear Clocks, Pendulum
Post by rloginunix on Jul 17th, 2015, 5:14pm
Use the basic trigonometric identity Sin2(x) + Cos2(x) = 1 and (p.12) for the cosine in the denominator of (p.13):

Cos(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif/2) = http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/surd.gif(1 - Sin2(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif/2)) = http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/surd.gif(1 - k2Sin2(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/phi.gif))


And put it back into (p.13):

dhttp://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif = 2kCos(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/phi.gif)dhttp://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/phi.gif/(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/surd.gif(1 - k2Sin2(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/phi.gif)))     (p.14)


Use (p.12) to rearrange the difference of squares of sines under the square root in (p.10):

Sin2(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif0/2) - Sin2(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif/2) = k2(1 - Sin2(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/phi.gif)) = k2Cos2(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/phi.gif)     (p.15)


Change the boundaries of integration using (p.12):

http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif = 0, Sin(0/2) = 0 = kSin(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/phi.gif), http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/phi.gif = 0

http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif = http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif0, Sin(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif0/2) = Sin(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif0/2)Sin(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/phi.gif), http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/phi.gif = http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/pi.gif/2


Put (p.14) and (p.15) into (p.10) and reverse substitute http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/omega.gif keeping the boundaries of integration as 0 and http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/pi.gif/2:

T(k) = (4/http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/omega.gif)http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/int.gifkCos(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/phi.gif)/http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/surd.gif(1 - k2Sin2(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/phi.gif)) * dhttp://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/phi.gif/(kCos(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/phi.gif))


T(k) = 4http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/surd.gif(L/g)http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/int.gifdhttp://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/phi.gif/(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/surd.gif(1 - k2Sin2(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/phi.gif)))     (p.16)


Though a French mathematician A. Legendre was working with functions of complex variables z (circa 1825-1837) with a substitution z = Sin(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/phi.gif) he was able to simplify a somewhat different but similarly looking integral to the form you see above. He called it an elliptic integral of the first kind and designated it as F(k). Legendre normalized two other related integrals to what now known as a standard form. All three carry his name, the elliptic integrals of the first, second and third kinds in the Legendre form.

When the upper boundary of the integrals is http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/pi.gif/2 they are called complete, otherwise, if it is an arbitrary angle http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/phi.gif, they are called incomplete and designated as a function of two arguments, http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/phi.gif being the second argument. For example, F(k, http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/phi.gif) stands for an incomplete elliptic integral of the first kind.

These integrals have been studied to a great depth and their numeric values can be estimated to a high degree of accuracy in a multitude of ways. Using Legendre polynomials and the works of J. Wallis, for example, our particular integral of the first kind can be exactly represented as an infinite sum:

F(k) = http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/int.gifdhttp://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/phi.gif/(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/surd.gif(1 - k2Sin2(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/phi.gif))) = http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/pi.gif/2 http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/sum.gif [(2n)!/(2n * n!)2]2 * Sin2n(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif0/2)     (p.17)


where n is a whole number: 0, 1, 2, 3, ... etc. which runs from zero to infinity. Putting it back into (p.16) we get:

T(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif0) = 2http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/pi.gifhttp://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/surd.gif(L/g)http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/sum.gif [(2n)!/(2n * n!)2]2 * Sin2n(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif0/2)     (p.18)

Title: Re: Mechanical Pendulum Gear Clocks, Pendulum
Post by rloginunix on Jul 17th, 2015, 5:17pm
It may be difficult to see (above) but notice that 4 became 2 and http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/pi.gif was added. The first few terms of the above sum are:

T(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif0) = 2http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/pi.gifhttp://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/surd.gif(L/g)(1 + (1/4)Sin2(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif0/2) + (9/64)Sin4(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif0/2) + (25/256)Sin6(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif0/2) + ...)


In clock making, however, the pendulum's period T is known ahead of time while the length of its rod L corresponding to that period must be found. If we decide to use (p.17) for the integral estimate then we can designate the sum in (p.18) as S:

T = 2http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/pi.gifShttp://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/surd.gif(L/g)


and solve it for L:

L(T, http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif0) = (gT2)/(4http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/pi.gif2S2)     (p.19)


where you have to decide on the number of terms in the sum S. In a generic case that sum will be a function of http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif0.

If you find a different way to calculate the integral (p.16) then we can name it as I:

T(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif0) = 4Ihttp://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/surd.gif(L/g)


and solve it for L:

L(T, http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif0) = (gT2)/(16I2)     (p.20)


Keep in mind that http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif0 must be expressed in radians for the above calculations to work. It is highly likely that you will measure or decide on the value of http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif0 in degrees d so here is how to convert them into radians r:

r = (http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/pi.gif/180) * d     (p.21)


Before moving on to small angles I will remark in passing that traditionally, after obtaining the integral in (p.16), we hope to calculate it and then invert that solution to obtain an expression for http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif(t). Unfortunately the analytic or closed form solution for the integral in (p.16) does not exist. However, the behavior of the angle http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif as a function of time can obtained via the elliptic function sn() of Carl Jacobi who, after inverting the integral in (p.16) in particular and further developing the theory of elliptic functions in general, used to famously say "One must always invert!":

http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif(t) = 2Arcsin(k * sn(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/surd.gif(g/L)t, k))


where k is some periodic function. This, however, is way beyond the scope of this tutorial.

Title: Re: Mechanical Pendulum Gear Clocks, Pendulum
Post by rloginunix on Jul 17th, 2015, 5:22pm
Small Angles

There seems to be no end in sight to wonders in math which uncovers things that are true but neither intuitive nor obvious. One such wonder is the fact that functions that have nothing to do with polynomials can be expressed via the very things they have nothing to do with - polynomials of infinite degree.

In 1712 an English mathematician Brook Taylor formalized the work done by a Scottish mathematician James Gregory in a form of a theorem, published in 1715, which states that a function possessing derivatives of all orders at a point x = a can be expressed as a polynomial whose coefficients depend only on the values of derivatives of the function at a. Such representation of a function, in general, is now called Taylor series.

Infinity is a convenient mathematical concept but in practice of course you will not be summing an infinite number of terms and you will not be raising a term to an infinite power. You will pick some reasonable finite number of terms and by doing so you will introduce a certain error. The upper bound of that error is addressed by the Taylor theorem but for the purposes of our discussion it is an overkill to cover it here. We will, however, look at the small angle error in the pendulum's length calculation below.

What makes Taylor series quite remarkable is the fact that we can estimate a value of a function at a certain point (I am glossing over a large amount of gory details here) using only four basic operations - addition, subtraction, multiplication and division - something that computers are good at.

An approximating series of a function in the vicinity of a specific point, zero, is sometimes called a Maclaurin series in memory of a Scottish mathematician Colin Maclaurin, Isaac Newton's disciple. Defending his teacher's analytic methods against the bitter criticism by the contemporaries Maclaurin published one of the first known systematic treatments of the subject, "A Treatise of Fluxions", in 1842. This work was so well and carefully done that it became a sort of an academic standard in Europe for the next almost 80 years. In 1821 a French mathematician Augustin Cauchy published his famous analysis course that introduced a level of mathematical rigor never dreamed of by Newton and Leibniz and, pretty much, put mathematical analysis as we now know and love it on the academic map.

Our prime suspect in (p.3) is Sin(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif). Keeping the angle http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif expressed in radians let us use Taylor (Maclaurin) series in the vicinity of http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif = 0 to represent Sin(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif) as a polynomial (of infinite degree, here n runs from zero to infinity):

Sin(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif) = http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif- http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif3/3! + http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif5/5! - http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif7/7! + ... = http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/sum.gif (-1n)/((2n+1)!) * http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif2n+1     (p.22)


To get a sense for the actual numbers here is what some small and not so small angles look like:

http://www.ocf.berkeley.edu/~wwu/YaBBAttachments/rlu_pgc02.png

Further, here is what (p.22) looks like for 4 degrees expressed in radians:

4 degrees = (3.14159 * 4)/180 = 0.0698131111 rad

Sin(0.0698131111) = 0.0698131111 - 0.0000567100 + 0.0000000138 - ...

Title: Re: Mechanical Pendulum Gear Clocks, Pendulum
Post by rloginunix on Jul 17th, 2015, 5:26pm
We see that starting from as close as the second term the numbers become so small and contribute so little to the grand total that we can throw all of them out and keep the first term only:

Sin(0.0698131111) http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/approx.gif 0.0698131111


or in general for very small angles:

Sin(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif) http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/approx.gif http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif


If we now put the above approximation of Sin(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif) back into (p.3) we get:

http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif'' + (g/L)http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif = 0     (p.23)


I leave it as an exercise for you to solve this equation using the separation of variables approach described above. Eventually you will come to the following integral:

http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/int.gifdhttp://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif/http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/surd.gif(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif02 - http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif2)


which you can calculate with the http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif= http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif0Sin(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/phi.gif) substitution to obtain the period T:

T = 2http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/pi.gifhttp://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/surd.gif(L/g)     (p.24)


which is the limiting case of (p.18) when only the first term of the sum is kept. We now see why for small angles the period of our idealized pendulum does not depend on http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif0 (and m).

Note that if we replace the angular coordinate http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif with a generic coordinate q then purely mathematically we can abstract (p.23) as:

q'' + http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/omega.gif2q = 0


which is applicable to: our simple pendulum, a point mass oscillating on an ideal spring, a current flowing through an LC circuit. For a simple pendulum q means angle and

http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/omega.gif = http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/surd.gif(g/L)


For a point mass on an ideal spring q means linear distance and

http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/omega.gif = http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/surd.gif(k/m)


For a current flowing through an LC circuit q means current and

http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/omega.gif = 1/http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/surd.gif(LC)


All three examples can be covered under one umbrella of the concept known as simple harmonic oscillator for which the behavior of the abstract coordinate as a function of time can be expressed in closed form. When you were solving (p.23) for your substitution for an arbitrary time t you had http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/phi.gif = http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/omega.gif t and for the angle you had http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif(t) = http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif0Cos(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/omega.gif t + http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/delta.gif). Replacing http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif with q we get:

q(t) = q0Cos(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/omega.gif t + http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/delta.gif)


where http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/delta.gif is known as a phase which captures the initial state of motion.

Title: Re: Mechanical Pendulum Gear Clocks, Pendulum
Post by rloginunix on Jul 17th, 2015, 5:29pm
g, http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/pi.gif, Reality

While other factors affect the exact value of g the two major ones are elevation above or submersion below sea level and latitude.

Elevation/submersion. The value of g becomes smaller and smaller as you get farther and farther away from the surface of the Earth in either direction - away from or closer to the center of the Earth. And conversely. It may be intuitively clear why the value of g diminishes as you get away from the Earth's center but why does it happen when you get closer to it? At the risk of straying off topic too much here is an explanation in broad strokes.

Break a solid sphere of uniform density into a very large number of very thin and hollow shells of uniform density. Because of the equality of the angles involved and the inverse proportionality of the gravitational forces to the squares of distances the gravitational pull of such a shell on a point mass located anywhere inside the shell is zero. A gravitational pull of a uniform sphere of mass M on a point mass located anywhere outside the sphere is equivalent to a pull of a material point of mass M located at sphere's center. Apply this model to Earth. As we get deeper and deeper inside our planet we bypass more and more such shells which we can safely throw out - they contribute nothing to the current value of g. The mass of the remaining not yet penetrated sphere diminishes and so is g.

Latitude. Earth is not a perfect sphere - it can be approximated with a high degree of accuracy with a spinning oblate ellipsoid. As such Earth's radius varies with your clock's position on our planet - it is smallest at the poles and largest at the equator. Consequently, the value of g starts out at the highest value at the poles, gp, and diminishes slightly as you get closer and closer to the equator, geq. And conversely:

gp = g(90) = 9.83236 m per s2

g = g(45) = 9.80665 m per s2

geq = g(0) = 9.78033 m per s2


Here degrees stand for the latitude and the unsubscripted g stands for standard gravity whose value is used for definitions and unit conversions. You may experiment with using it in your calculations but you should really use the clock's locale-specific value. For example, you can look it up in a published table - the burden to find one is on you. Keep in mind that when it comes to measurements lone magnitudes like 9.806 by themselves are useless. They become worthy of attention only when an error is given as in 9.806 http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/pm.gif 0.003. Such errors tell us the degree of uncertainty of a measurement. However, I will mostly deviate from this practice in this tutorial.

In US a National Oceanic and Atmospheric Administration  (www.noaa.gov) has a branch called National Geodesic Survey  (www.ngs.noaa.gov) which, at the time of this writing (2015), provided a Surface Gravity Prediction calculator  (www.ngs.noaa.gov/cgi-bin/grav_pdx.prl) which accepted latitude, longtitude and height as input for continental US only. For New York City, NY, US (latitude 40.714167 degrees, longtitude 74.006389 degrees, elevation 4.66 m) it produced the following estimate:

g(NYC, NY, US) = 9.80252 m per s2


If for calculation purposes you need to convert the coordinates from the minutes M and seconds S format into a decimal format d then use the following formula:

d = http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif + M/60 + S/3600


For NYC, NY, US for example, we get:

dlat = 40 + 42'/60 + 51''/3600 = 40.714167 degrees

dlon = 74 + 0'/60 + 23''/3600 = 74.006389 degrees


You can also estimate the local value of g at sea level by using the Helmert's formula:

g(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/phi.gif) = 9.780327(1 + 0.0053024Sin2(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/phi.gif) - 0.0000058Sin2(2http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/phi.gif))     (p.25)


where http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/phi.gif is your latitude in degrees.


[e]Typo fix in g definitions.[/e]

Title: Re: Mechanical Pendulum Gear Clocks, Pendulum
Post by rloginunix on Jul 17th, 2015, 5:33pm
If your clock will be operated above sea level at a considerable height h, measured in meters, then adjust (p.25) by subtracting an extra term from it:

g(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/phi.gif, h) = g(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/phi.gif) - 0.000003086 * h     (p.26)


Using (p.26) to estimate the value of g in NYC, NY, US for comparison purposes, we get:

g = 9.80232 m per s2


If your clock will be operated below sea level at a considerable depth d, measured in meters, then adjust (p.25) by multiplying it by an extra term:

g(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/phi.gif, d) = g(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/phi.gif)(1 - d/Rm)     (p.27)


where Rm = 6,371,009 meters is the mean radius of Earth. If you feel that Rm is not precise enough then you can look up any number of ways to estimate the distance from the center of the Earth to its surface at an arbitrary latitude http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/phi.gif. Wiki is one place to start from.

Overall, the chase after a more precise estimate of a locale-specific value of g is always on - check for the latest developments before running your calculations.

Also, keep in mind that by making it possible for a massive object, bob, to slide along the pendulum's rod you can always fine-tune the period of your clock. If your clock runs fast then you slow it down by sliding the mass m down the rod away from the pivot point effectively increasing the magnitude of L. If you clock runs slow you speed it up by sliding the mass m up the rod towards the pivot point effectively decreasing the magnitude of L.

The value of http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/pi.gif is known with a very high degree of accuracy - way more than we need for our practical purposes. I will use http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/pi.gif http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/approx.gif 3.14159.

Lastly, it goes without saying that our nice mathematical model is not realizable in practice. In reality the period of oscillation of a pendulum is affected by the pivot assembly friction, air resistance, atmospheric pressure, temperature, humidity, etc. I will address the influence of friction and air resistance on T in the escapement chapter but I will ignore the influence of all the other factors on T completely.


Procedure

Here are the steps to find the length of the simple pendulum's rod L:

Decide on the magnitude of T in seconds.
Decide on the magnitude of http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif0 in degrees.
Convert degrees into radians with (p.21).
Estimate the value of the integral (p.16).
Find the locale-specific value of g.
Take http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/pi.gif http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/approx.gif 3.14159.
Put the above values into (p.19) or (p.20) to find L.


Sample Calculation

1). We pick T = 2 seconds. It means that it takes our pendulum 1 second to swing from A to B and 1 second to swing from B to A. In other words this particular pendulum shows up at either extreme point every second. Historically it is known as seconds pendulum and professionals say that it "beats seconds".

2). We pick http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif0 = 2 degrees.

3). From (p.21) we obtain its radian measure:

http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif0(2d) = (3.14159/180) * 2 = 0.03491 rad


4). At first we shall do a small degree approximation by taking only the first term, 1, from S in (p.18) and from (p.19) we get:

S = 1


5). We will use the value of g = 9.80252 specific for NYC, NY, US.

6). Put these numbers into (p.19) to obtain L:

L(2s, 2d) = 9.80252/(4 * (3.14159)2) * 4 = 9.80252/9.86959 = 0.99320 m

L = 99.32 cm


In other words we get 99 centimeters and 3.2 millimeters. On a separate note observe a numeric curiosity - the values of g and http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/pi.gif squared are somewhat close.


[e]changed l to L in 6)[/e]

Title: Re: Mechanical Pendulum Gear Clocks, Pendulum
Post by rloginunix on Jul 17th, 2015, 5:34pm
Now instead of just one let us take three first terms in S in (p.18) still using the Sin() in it (not its small angle approximation):

http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif0(2d)/2 = 0.01745

Sin(0.01745) = 0.00030

Sin2(0.01745) = 0.00000

Sin4(0.01745) = 0.00000


S = 1 + (1/4) * 0.00000 + (9/64) * 0.00000 = 1


In other words with 2 degrees we do not get any perceptible difference for T. What about 4 degrees?

http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif0(4d)/2 = 0.0349065555

Sin(0.0349065555) = 0.0348994671

Sin2(0.0349065555) = 0.0012179728

Sin4(0.0349065555) = 0.0000014834


S = 1 + (1/4) * 0.0012179728 + (9/64) * 0.0000014834 = 1.0003047018

S2 = 1.00030470182 = 1.0006094964

L(2s, 4d) = 9.80252/(4 * (3.14159)2) * 4/1.0006094964) = 0.99260 m

L = 99.26 cm


In other words we get 99 centimeters and 2.6 millimeters. If we take it as a correct value than our previous one, 99.32 cm, which is off by 0.6 mm, introduces about 0.06% error.

In a similar manner for small angles (S = 1) for NYC, NY, US we obtain:

L(0.5s, 4d) = 6.21 cm

L(1s, 4d) = 24.83 cm

L(3s, 4d) = 2 m 23.47 cm

L(4s, 4d) = 3 m 97.28 cm

L(5s, 4d) = 6 m 20.75 cm


Overall keep in mind that the lengths of a small angle simple pendulum relate as squares of its periods. If you double a period (from 2 to 4 seconds) you quadruple the length (from about 1 meter to about 4 meters) and conversely.


[e]Changed l to L for 4 degrees.[/e]

Title: Re: Mechanical Pendulum Gear Clocks, Pendulum
Post by rloginunix on Jul 17th, 2015, 5:36pm
Huygens

By around 1637 Galileo realized that the pendulum can be used as the main part of a time keeping mechanism - a clock. He sketched out a preliminary design of such a clock and started implementing it by the summer of 1641. Unfortunately he never finished it as he passed away in January of 1642. The scientific marathon, however, was not to die with him - on December 25th of that year Isaac Newton was born. In England.

A friar of a Parisian order of Minims, the author of the special kind of numbers of the form 2n - 1, left a well known trace in number theory. It is less well known that while running a hard sciences informational hub out of the rooms of his Franciscan convent Marin Mersenne left a trace in physics as well as he was the one who suggested to a Dutch scientist and mathematician Christiaan Huygens the use of pendulums in time keeping.

The devil is in the details and it took the intellectual power of Huygens to work these details out around 1656-1657 and invent a first known practical rigid rod pendulum clock. A physical functioning copy of this clock, built by someone else, was unveiled by Huygens on December 25th 1657 while its inner workings were described in 1658 in a small book "Horologium". Newton was a tender 15 years of age but a good friar Marin Mersenne did not live to witness this historic moment as he passed away in 1648.

Eventually Huygens realized that a rigid rod pendulum tracing a circular arc is not isochronous (from the Greek "iso" meaning equal and "chrono" meaning time) as its period of oscillation does depend on the size of the swing. In 1673, while in Paris on an invitation from one of the ministers of the French king Louis XIV, in a more voluminous book "Horologium Oscillatorium" he published a solution to that problem by finding a curve tracing which an oscillating point has a period truly independent of the swing size - even for large initial angles. Since Newton's "Principia" and calculus were not to be available for another 14 years (1687) Huygens relied on his personal insight and geometric constructions for the proof.

Such a curve is called a tautochrone from the Greek "tauto" meaning same. It turned out that in a uniform gravitational field a tautochrone is a cycloid. A cycloid is a two-dimensional curve traced by a point on a circumference of a circle rolling along a stationary straight line without slipping. Further, no matter where a material point is initially placed on a properly positioned cycloid (horns up) it will always reach the equilibrium point O in equal amounts of time.

Based on the new find Huygens devised a different pendulum. Instead of using a rigid rod he suspended a point mass at the end of a flexible thread whose swings he limited by the semicycloidal cheeks. Finding the exact shape of these cheeks - a non-trivial task at the time - was yet another discovery made by Huygens. The essence of this discovery is that the evolute of a given cycloid is an equal but shifted cycloid (an evolute is the locus of all the centers of curvature of a given curve).

At this point Huygens was assured that while the swinging thread hugs these cheeks the point mass has no choice but to trace a cycloid (he even went as far as devising a few practical methods of constructing a cycloidal curve). As such we can say that a cycloidal pendulum is isochronous.

I will remark in passing that with calculus we can find the parametric equation of a tautochrone by imposing a requirement that the restoring force mgSin(http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/theta.gif) must be linearly proportional to the length of the path traveled by the point mass along the curve.

However brilliant Huygens' theoretical discovery was unfortunately it did not work as well in practice. For large swings the flexible thread had so much contact with the curbing cheeks that it resulted in a lot of friction which in turn led to a rather rapid energy loss. For small swings the actual contact of the flexible thread with the required trajectory enforcing cheeks was estimated to be just north of one millimeter - not enough to exercise the full benefits of the cycloidal path (A. L. Rawlings, The Science of Clocks and Watches). There were purely technological difficulties as well - at the time it was challenging to manufacture the curbing cheeks in a mathematically precise cycloidal shape. Rigid rod pendulums coupled with the escapement mechanism took the center stage in clock manufacturing.

The book "Horologium Oscillatorium" is highly recommended. If you do not now Latin, the book's original language, either learn it or seek a copy in a tongue you understand. Huygens' style of clarity and conciseness is contagious. I list the translations of the book into several languages at the end of the Hands chapter.


[e]Fixed Mersenne's spelling.[/e]

Title: Re: Mechanical Pendulum Gear Clocks, Pendulum
Post by rloginunix on Jul 17th, 2015, 5:40pm
More Pendulums

Over the years many different types of pendula were invented, each one for a specific purpose. I will mention just a few: O Botafumeiro, military, Kater's, Foucault's, coupled, torsion, chaotic pendula and - quantum pendulum, no less, first studied by Edward Condon (1928) though I am not very ... certain.

After you master a simple mathematical pendulum you may want to take a look at the pendulums called physical or compound. In these pendulums the geometry of the oscillating mass can no longer be ignored - it takes the center stage. I will omit the deduction process and only show the end result of practical use. A period T of a small angle compound pendulum depends on its moment of inertia Iz about the pivot point Z (in our notation) and the distance d from the pendulum's center of mass to the pivot point:

T = 2http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/pi.gifhttp://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/surd.gif(Iz/mgd)     (p.28)


Use the Stiener-Huygens (parallel axis) theorem to calculate Iz.

For compound pendulums the geometric parameter of interest to us will be packed into Iz and d and the above equation will have to be solved anew for each particular case. Here are a few simple practical ones. All the bodies below are absolutely rigid, solid, thin and uniform. Keep in mind that in comparison with the simple pendulum the characteristic sizes of the compound pendulums come out quite large and will require more vertical and horizontal clearance:

http://www.ocf.berkeley.edu/~wwu/YaBBAttachments/rlu_pgc03.png

A rod

A rod of length L and mass m swinging about a point at one of its extremities. The distance d = 0.5L and the rod's moment of inertia about Z is:

Iz = mL2/12 + mL2/4 = mL2/3


From (p.28) we get:

T = 2http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/pi.gif http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/surd.gif((2mL2)/(3mgL)) = 2http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/pi.gif http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/surd.gif((2L)/(3g))


Solve for L:

L(T) = 3gT2/(8http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/pi.gif2)


For a 2-second NYC, NY, US pendulum for L we get a rough estimate of 148.98 cm.


A Hoop

A hoop of radius R and mass m swinging about a point on its circumference. The distance d = R and the hoop's moment of inertia about Z is:

Iz = mR2 + mR2 = 2mR2


From (p.28) we get:

T = 2http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/pi.gifhttp://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/surd.gif((2mR2)/(mgR)) = 2http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/pi.gifhttp://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/surd.gif((2R)/(g))


Solve for R:

R(T) = gT2/(8http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/pi.gif2)


For a 2-second pendulum for R we get a rough estimate of 49.66 cm. The hoop's characteristic size then is about one meter.


A Disk

A disk of radius R and mass m swinging about a point on its circumference. The distance d = R and the moment of inertia of a disk about Z is:

Iz = mR2/2 + mR2 = 3mR2/2


From (p.28) we get:

T = 2http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/pi.gifhttp://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/surd.gif((3mR2)/(2mgR)) = 2http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/pi.gifhttp://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/surd.gif((3R)/(2g))


Solve for R:

R(T) = gT2/(6http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/pi.gif2)


For a 2-second pendulum for R we get a rough estimate of 66.21 cm. In comparison with a hoop we need more disk for the same period.


A Triangular Boundary

Three rods each having the length L and mass m form a rigid equilateral triangle which is swinging about one of its vertexes. The moment of inertia of the perimeter of the equilateral triangle (the boundary only) about Z is:

Iz = 2 * mL2/3 + mL2/12 + m3L2/4 = 3mL2/2


To find the location of the center of mass of a triangular boundary known as Spieker point or Spieker center we replace the massive rods with the equivalent point masses located at the rod's centers. We construct a medial triangle using these mid points. The intersection of the inner angles' bisectors of the medial triangle locates the center of mass. In our case of an equilateral triangle we get:

d = (2/3) * (Lhttp://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/surd.gif3)/2 = (Lhttp://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/surd.gif3)/3


From (p.28) we get:

T = 2http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/pi.gifhttp://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/surd.gif((3mL2)/(2 (3m)g) * (3)/(Lhttp://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/surd.gif3)) = 2http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/pi.gifhttp://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/surd.gif((Lhttp://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/surd.gif(3))/(2g))


Solve for L:

L(T) = gT2/(2http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/surd.gif(3)http://www.ocf.berkeley.edu/~wwu/YaBBImages/symbols/pi.gif2)


For a 2-second pendulum for L we get a rough estimate of 114.68 cm.


[e]Fixed superscript typos in L(T).[/e]

Title: Re: Mechanical Pendulum Gear Clocks, Pendulum
Post by rloginunix on Jul 23rd, 2015, 5:31pm
The Escapement chapter is here. (https://www.ocf.berkeley.edu/~wwu/cgi-bin/yabb/YaBB.cgi?board=riddles_general;action=display;num=1437253052)



Powered by YaBB 1 Gold - SP 1.4!
Forum software copyright © 2000-2004 Yet another Bulletin Board