Topics

- Cartesian geometry
- Archimedes trammel
- The cycloid
- Watt's linkage curve

Further reading

- Engine linkage examples
- Does Apollonius deserve some credit for early development of cartesian geometry?
- History of the function concept

Cartesian geometry, also known as analytic geometry, united the subjects of algebra and geometry by interpreting the variables of algebraic equations in terms of a coordinate system to create geometric objects. The methods were developed by rivals Rene Descarte and Pierre Fermat in the early 1600's, with prior art being found Mercator's world map (1569), and in the earlier works by Oresme (circa 1350), Menaechmus (circa 350 BC), among others. The introduction of algebra and coordinates to geometry made it possible to quickly describe many more curves without having to rely on special names and mechanical constructions. Cartesian geometry had a subtle but even more profound effect on our fundamental perception of curves. It broke the symmetry of Euclidean geometry where there was no concept of absolute position and only relative distances mattered. Cartesian geometry lead naturally to the asymmetric relationship of independent and dependent variables so that axes could be *different* units like time or velocity. It set the table for the development of differential calculus within a generation, and remains foundational today.

An important application of Cartesian geometry is the formulation of equations that describe the mobility of mechanical systems. Engineers made huge advances in the design of machines during the industrial revolution (1770-1830), but these advances could only go so far without mathematics. The creation of machines with many interconnected moving parts like the steam engine and the Jacquard loom required a degree of precision in design and construction that had never been needed previously. Cartesian geometry, together with trigonometry, provided a language for this precision engineering.

Many geometric laws can be used to derive mobility equations, but two are of primary importance. The first is the Pythagorean theorem: for any right triangle, the sum of the squares of the length of the legs is equal the square of the length of the hypotenuse. The second is the law of similar triangles: any two triangles are similar if their angles are equal. The ratios of the lengths of corresponding sides of similar triangles are always equal to each other. When applied in the context of Cartesian geometry, it is sometimes convenient to work with the squares of these ratios to avoid complications of signed magnitudes the algebra introduces.

At night in the summer, we sometimes see bicycles ride by with lights upon their wheels. The lights revolve in a circle around the wheel axles, but the path the lights trace appear as cheerful looping bounds down the road as the bicycles move on. If the light is on the rim of the wheel, the curve it traces is called a cycloid. The cycloid's modern study began around 1500, and it became closely connected to many of the interesting mathematical arguments of the next two centuries, including the quadrature, the pendulum clock, optics, and the brachistochrone -- it was nicknamed the "Helen of mathematics" in its time because of the wars it ignited.

A planar curve can be represented in Cartesian geometry implicitly as a function \(f(x,y) = 0\), parametrically as \(x = f(\theta)\), \(y = g(\theta)\), or explicitly as a function \(y = f(x)\). Each has its own preferred uses and limitations, as we will see.

To find an equation for the cycloid, we can make use of parametric coordinates. Let's suppose we start with a wheel of radius \(R\) rolling on a level and flat surface, and we are tracking a bangle a distance \(r\) from the wheel axle. As a wheel rolls, its current state can be represented with a variable \(p\) tracking the position of the axle of the wheel, and a variable \(\theta\) tracking the angle of rotation of the wheel. As long as there is no slippage of the wheel, these two must be related -- each time the wheels angle makes a full revolution of \(2 \pi\), the wheel's position must move forward by a factor of the circumference of the wheel. Since we assumed our wheel has radius \(R\), one rotation of the wheel will move the axle forward \(2 \pi R\). If, initially, \(p=0\) and \(\theta=0\), then \(p=\theta\) as the wheel rolls, whatever the velocity or acceleration.

Now, let \((x,y)\) represent the position of a bangle on our wheel, a distance \(r\) from the axle. Since the axle's position can be expressed in terms of the angle as \((p,1) = (\theta R, 1)\), we can express the position of the bangle with trigonometry as \[x(\theta) = R \theta - r \sin \theta,\quad y(\theta) = R - r \cos \theta,\] assuming the initial position is directly below the axle.

Another example of Cartesian geometry applied to a mechanical system is one you might find in a kid's toy chest -- Archimedes' Trammel, shown above. The handle end is free to be turned, but only to positions that the movement of the two sliders allows. Experimentation with the trammel reveals the handle's motion traces an oval shape taller than it is wide. A natural guess, for those familiar with geometry, is that the motion traces an ellipse. But can we be certain of this?

To answer the question, we can use some geometry and algebra. Let's set the center of the trammel (where the slider grooves cross) to be the origin of our system, and let \((x,y)\) be the position of the handle. In Cartesian coordinates, an ellipse centered at the origin and with principle axes aligned with the coordinate axes has an equation of the form \[a x^2 + b y^2 = c,\] where \(a\), \(b\), and \(c\) are positive constants determining the ellipse's shape. Every point \((x,y)\) on the ellipse will be a solution of this equation. If we can show that there are some constants \(a\), \(b\), and \(c\), such that \(x\) and \(y\) always satisfy this equation, we will know the trammel traces an ellipse.

On close inspection, we see the trammel consists of two sliders trapped in perpendicular grooves. Each slider can move freely, but every time we move one slider, the other slider gets moved in response. The handle's screw-attachments control this -- the screws are always the same distance apart. So when one slider gets moved out, the other sider gets moved in to compensate for the distance change.

There are two parameters in the trammel construction that may change from one toy to another, but are never changed once the trammel is completed -- the distance between the handle and the first screw, and the distance between the first and second screws. Let \(r\) be the distance between the handle and the first screw and let \(s\) be the distance between the two screws.

The things that can change in the trammel are the positions of the slider screws and the position of the handle. Let \(p\) be the horizontal slider screw's position and \(q\) be the vertical slider screw's position while \((x,y)\) is the handle's position. There are 3 equations that constrain the motion of our 4 free variables. By the Pythagorean theorem, we always have \[s^2 = p^2 + q^2.\] And we have two similar triangles in action in the trammel, and since corresponding sides of similar triangles are always in the same proportions, we have \[\frac{q}{y} = \frac{p}{p-x} \quad \text{and} \quad \frac{q}{y} = \frac{s}{s+r}.\] The third equation we obtain from the similar triangles is redundant because it can be recreated from the other two.

It is convenient to write these equations as a system of polynomials, each equal to zero, to which, \[\begin{align*} 0 &= s^2 - p^2 - q^2, \\ 0 &= q \,(p - x) - p y, \\ 0 &= q \,(s +r) - s y, \end{align*}\] in the four variables \((p,q,x,y)\), with \(s\) and \(r\) given as constants. These equations completely specify the mobility of Archimedes's trammel.

The mobility equations of Archimedes's trammel formed a nonlinear system of 3 in 4 variables, so in the local neighborhood of a point on the solution, we expect there to be 4-3 = 1 degree of freedom -- a curve. But globally, there may be a single curve or multiple curves, and each curve might be bounded set or go off to infinity.

For our trammel, we are looking for a curve that relates the handle position to the trammel's built-in parameters, without regard to the slider positions -- something like a function that can be written \(f(x,y; r,s) = 0\) without using \(p\) and \(q\). Well, the proportions can be solved for \(p\) and \(q\). \[q = \frac{ys}{s+r}, \quad p = -\frac{sx}{r}.\] Substituting for \(p\) and \(q\) in the Pythagorian theorem, we find \[s^2 = \left( \frac{sx}{r} \right)^2 + \left(\frac{sy}{s+r}\right)^2\] \[r^2 = x^2 + \left(\frac{1}{\frac{s}{r}+1}\right)^2 y^2\] Thus, the path traced by the trammel's handle will be an ellipse. We can go a step further, and observe that as the ratio \(s/r\) gets smaller, the motion will become less elliptic and more circular.

For the trammel, we were able to obtain an implicit solution, while for the cycloid, we were able to obtain a parametric solution, Implicit representations of a curve can easily be used to test if a point falls on the curve, but does not as easily draw the curve. On the other hand, parametric representations can easily draw the curve, but do not as easily test if a given point is a solution or not. To test if a point is on the cycloid using this parametric formula, we have to search through many different angles \(\theta\), testing if the points match at each step. A quicker method for checking if a point is on the cycloid is to use an inverse cosine function to calculate \(\theta(y)\), and to then compare the given and formula values for \(x\). This will still require a little extra work. Inverse cosine will have infinitely many answers that will have to be sorted through.

An explicit function representation, which avoids the pitfalls of implicit and parametric representations by using a coordinate as the free parameter, is the most useful because it can easily test if a given point falls on the curve, as well as be used to draw a set of points on the curve. However, as we see above, an explicit form does not always exist.

Sometimes, we can make an implicit solution explicit using local Taylor series approximations. However, these local approximations have to be used carefully to avoid errors that come form accidentally applying them outside their domains of validity. We'll see this in the next example.

In an important sense, the analytic formula's we derive are all parametric solutions because the geometry is characterized by a set of parameters which then appear explicitly in our solution equations. These analytic solutions capture parameter-dependencies, which allows for exploratory sensitivity analysis. This is vital in design, but impossible using non-parametric approaches like Snel's triangulations, where the equations are formulated directly in terms of the data.

A classic application of Cartesian geometry in mechanical engineering is in the characterization of Watt's linkage. James Watt (1736-1819) was a successful Scottish engineer famous for his improvements to steam engine design. One of the primary uses of his engines was to pump water out of coal mines. Early pumps tended to fail regularly, not because of breakdowns in the steam-engine itself, but because of uneven wear on the pump piston. In early designs, piston shafts were linked to the steam engine though a large beam that would drive the piston up and down. But, since the beam-end's path traced part of a circular arc, the piston would rock back and forth in a manner that would systematically wear out the piston seals.

Watt realized that the piston wear could be avoided and reliability increased if a linkage could be found that would move the piston rod up and down in a straight line. He was so proud of his solution, he had it included in his portrait!

In Watt's linkage, two opposite bars in the linkage are of equal length, and the piston rod is connected to the midpoint of connecting rod. Specifically, suppose two pivots are fixed at \((-1,0)\) and \((1,0)\) respectively, and that bars connected to these pivots each have length \(r\). Let \((u_1,v_1)\) be the location of the free end of the first, while \((u_2,v_2)\) be the location of the free end of the second. These two bars are connected by a third bar of length \(\ell\), and the piston rod is connected to the point \((x,y)\) that is the midpoint of this third bar. Based on the opposite convexities of the circles described by the two opposite equal bars, Watt argued that the motion of the midpoint \((x,y)\) would be linear. However, it was almost 100 years before this conjecture was sorted out by Samual Roberts.

By directly translating this description of Watt's linkage into algebra, we end up with 5 equations involving 2 constant parameters (\(r\) and \(\ell\)) and 6 free variables (\(x, y, u_1, u_2, v_1\), and \(v_2\)). \[\begin{align*} \text{Bar 1 length:} && r^2 &= (u_1+1)^2 + v_1^2, \\ \text{Bar 2 length:} && r^2 &= (u_2-1)^2 + v_2^2, \\ \text{Bar 3 length:} && \ell^2 &= (u_2-u_1)^2 + (v_2-v_1)^2, \\ \text{Horizontal midpoint:} && 2 x &= (u_1 + u_2), \\ \text{Vertical midpoint:} && 2 y &= (v_1 + v_2). \end{align*}\] Since there are only 5 equations for 6 variables, we expect the system should have a single degree of freedom. However, in this form, the system is not very useful. There is no easy way to either iterate through points on the curve or to test if a given point falls on the curve. A natural approach to solving this system for midpoint location \((x,y)\) is to parameterize \(u _ 1\) and \(v _ 1\) in polar coordinates, and calculate the midpoint as a function of the angle. However, this approach becomes clumsy because of the repeated need to calculate the intersections of two circles. With some clever algebra (today, we use Groebner bases, a tool you may study elsewhere), one can show that \((x,y)\) solves the single sixth-degree equation

\[0 = 4 y^{2} \left(x^{2} + y^{2} - r^{2}\right) + \left(x^{2} + y^{2}\right) \left(\frac{\ell^{2}}{4} - r^{2} - 1 + x^{2} + y^{2} \right)^{2}.\]

[Show code]

```
from sympy import *
u1,u2,v1,v2,a,r,L = symbols('u1,u2,v1,v2,a,r,L')
sys = [ \
(u1 + u2) - 2*x, \
(v1 + v2) - 2*y, \
(u1 - a)**2 + v1**2 - r**2, \
(u2 + a)**2 + v2**2 - r**2, \
(u2 - u1)**2 + (v2 - v1)**2 - L**2, \
]
ans = (groebner(sys, (v1,v2,u1,u2,x))[-1]*q).collect((q,a),factor).subs(q,1)
pprint(ans)
```

To determine the parameter conditions when Watt's linkage mobility is approximately straight around its midpoint, we can apply a series solution method. We can check by substitution and inspection that \((x,y) = (0,0)\) lies on the curve. Now, using a substitution \[y \approx a_1 x + a_2 x^2 + a_3 x^3 + a_4 x^4 + a_5 x^5,\] and picking the \(a_i\)'s so that coefficients vanish one by one, we can show that when \(\ell = 2 \sqrt{1 - r^2}\), then the motion of the piston link is approximately give by \[y(x) = a_1 x + a_5 x^5 + O(x^6), \quad a_1 = \pm \frac{r}{\sqrt{1-r^2}}, \quad a_5 = \frac{-(1 + a_1^2)^3}{8 a_1 r^2 (1 - r^2)}.\] This shows that Watt's linkage motion is approximately a straight-line motion to fifth-order when the linking bar length \(\ell\) is picked right.

Write a python function to check if a given point \((u,v)\) falls on a cycloid curve with known parameter values \(R\) and \(r\), with an absolute error of no more than \(10^{-2}\). Use your function to test the following 3 points when \(R = 2\) and \(r=3\): \((-25.288,0.073)\), \((0,1.770)\), \((27.5,4)\).

- The implicit 6th degree equation obtained above for Watt's linkage is not very easy to plot. This can be partially solved by choosing a different parameterization and exploiting the linkage-solution's 4-fold symmetry.
- Take \(\alpha = x^2 + y^2\) and find a parametric solution of the form \((\pm x(\alpha), \pm y(\alpha,x))\).
- When you write a program to plot Watt's linkage using this parameterization, how well does it work?

In class we found that Archimedes' trammel traced out an ellipse, based on the separation between the pins and the handle. But it is not the only way to construct a perfect ellipse. Imagine two thumb tacks placed on cardboard a distance \(d\) apart. Take a string of length \(L\), tie it in a circle, loop it over the thumb tacks, and draw the curve of the farthest the string can reach. Show using Cartesian geometry that this curve is also an ellipse, as long as as \(L > 2 d\).

A rotary pump is built by fixing a circular disc to an axle off center, then using the disc to raise and lower a piston as the axle rotates. If the disc has radius \(r\) and the axle is offset a distance \(a < r\) from the center of the disc, find the function \(h(\theta)\) for the height of the piston as a function of the angle of rotation of the axle. Assume \(\theta = 0\) is when the center of the disc is directly below the axle and that \(h(0) = 0\).

- In a flat area of sand on a beach, a piling with circular cross-section of radius \(1\) foot is driven vertically into the ground. A rope has been wrapped tightly around this piling at ground level, with one end free. Grab the free end and walk around the piling to unwind the rope, while keeping the rope taut.
- Find an equation for the path the free end of the rope traces in the sand as you unwind it. (Hint: This curve is
**not**an Archimedean spiral.) - Show that the curve
**is**well-approximated by an Archimedean spiral asymptotically as the length of unwound rope gets larger using a plot.

- Find an equation for the path the free end of the rope traces in the sand as you unwind it. (Hint: This curve is
- The Conchoid of Nicomedes was a curve traced by a linkage mechanism where one pin is fixed to the vertical axis, while the second pin can move freely along the x axis. (see here) This conchoid was designed to help divide an angle into thirds, solving one of the classical Greek puzzles of geometry.
- What are the two physical length parameters that are determined by the linkage construction and control the curves shape. Call these parameters \(a\) and \(b\).
- Using \(x\) to represent the horizontal position of the movable pin, and \((p,q)\) representing the point of the linkage, find a pair of parametric curves for the Conchoid in the forms \((p(x,a,b), q(x,a,b))\).

The Peaucellier-Lipkin consists of 4 bars of length \(b\) forming a rhombus, and an isoceles triangle formed by two longer bars of length \(m\) joined to opposite corners of the linkage. If the peak of the triangle is fixed to the origin, and the inner corner of the rhombus is fixed to move on a circle of radius \(r\) centered at the midpoint between the origin and the inner corner of the rhombus, find a set of equations for the positions of the three other corners.

Find the equation for a torus using a 3-dimensional Cartesian coordinate system.

In class we found that Archimedes' trammel traced out an ellipse, based on the separation between the pins and the handle. But it is not the only way to construct a perfect ellipse. Imagine two thumb tacks placed on cardboard a distance \(d\) apart. Take a string of length \(L\), tie it in a circle, loop it over the thumb tacks, and draw the curve of the farthest the string can reach. Show using Cartesian geometry that this curve is also an ellipse, as long as as \(L > 2 d\)

Show that the intersection between a cylinder and a plane is an ellipse

- A curve traces the intersection point of two lines. The first line is flat, starts at height 1, and moves down with constant speed. The second line is through the origin, starts off vertical, and rotates at a constant angular speed. Both lines reach the x-axis at the same time.
- Find an equation for the curve tracing the intersection.
- Use L'Hopital's rule to find the intersection point of the curve with the x-axis.

(to be completed) A hexagonal trammel is similar to Archimedes's trammel, except that it consists of 3 sliders moving along the 3 axes of a regular hexagon. Model the motion of a point on the handle of this trammel.

(to be completed) Show that Leonardo's ellipsograph does indeed draw ellipses.

- Plot the cycloid path in the following cases.
- Plot the path when R=1 and r = 0.6 for 4 revolutions.
- Imagine that in our cycloid problem, we could extend the position of the bangle beyond the rim of the wheel.

Plot the path when R=1 and r = 1.5 for 3 revolutions.

- Take a 1-foot bar 1 and connect it to the origin. Take a second bar three feet in length and connect the free end of the first bar 1 foot from one end. Fix the free end of the long bar furthest from the pivot point so that it always lies on the negative x-axis, but slides freely. Let \((x,y)\) be the position of the opposite free end of the long bar.
- Find a system of equations that can be used to determine possible positions \((x,y)\).
- (Hard) Reduce your system to a single polynomial equation involving only \(x\) and \(y\).
- Observe that your solution curve is symmetric across the y-axis, even though the true solution does not have any negative values for \(x\). What happened?

- Small changes of the lengths of components in a mechanical linkage can sometimes lead to big changes in the range of motion of the linkage.

Consider the polynomial system \(0 = a - 4 x^2 + 2 v^2 + y^2\), \(0 = x^2 - v\) for the variables \(v,x,y\).- Use substitution to find a single equation relating \(x\) and \(y\) that is independent of \(v\).
- Show that if \(a\) is large, then there are no solutions of this system.
- Show that for any constant value of \(a\), the solution set is bounded.
- Show the solution set contains two unconnected curves when \(a = 1\).
- Show the solution set is a simply connected curve when \(a =-1\).

(Hard) Suppose you have the equations for a plane and a helix, both in 3 dimensions. How many solutions are there to this equation? Discuss their enumeration and calculation.

(Hard) Transform the equation for Watt's curve from Cartesian to polar coordinates. Use this and symmetry to plot the curve when \(r=5/7\) and \(\ell = 2 \sqrt{1-r^2}\).

(Hard) Explain Philon's construction of a cube roots in modern language.

- (Hard) A pair of sophomore students are moving a large mural painting on plaster board into a new apartment, but to get to the room, you have to turn the rectangular mural around a tight corner between two hallways. The mural is about as tall as the hallway, so tilting it only makes the fit harder.
- What 3 variables control whether or not you can get the couch
- Find an equation for the possible width of the first and second hallways for which you will be able to get the couch to your apartment.

There are actually two versions of Watt's linkage. The the second more practical version can be seen here. Draw and label a diagram of this second version of Watt's linkage and construct a system of equations to determine the mobility of the the top of the piston rod.

Show that the intersection between a cylinder and a plane is an ellipse.

Show that the intersection between a cone and a plane is an ellipse.

Find the equation for a torus using a 3-dimensional Cartesian coordinate system.

Design your own linkage of rods and pivots that turns circular motion into parabolic motion.

Find the shape of a curve over which a square could roll without slipping and keeping it's center level as it goes.