Involutes

Let $\alpha ,\beta :(a,b)\to {\mathbb{R}}^3$ be regular curves. We say that $\beta$ is an involute of $\alpha$ when

1. $\beta (t)$ lies on the tangent line to $\alpha$ at $\alpha (t)$, and 
2. the tangents to $\alpha$ and $\beta$ at $\alpha (t)$ and $\beta (t)$ are orthogonal.

In this post, we will explore the basics of what one can say about an involute of a given curve.

Let $\alpha (s)$ be a unit speed curve. Let $\beta$ be an involute of $\alpha$. What can we say about $\beta$?

Evaluate this Sage cell to get an example of a cuve and its involute.

xxxxxxxxxx

 

1

var('s n')

2

a(s) = [cos(s), sin(s)]

3

b(s) = [cos(s) + s*sin(s), sin(s) - s*cos(s)]

4

P = parametric_plot(a(s), [s,0,3*pi/2], color='black')

5

P+= parametric_plot(b(s), [s,0,pi], color='red')

6

string_list = [0,pi/8, pi/2, pi/3, 2*pi/3, pi]

7

for value in string_list:

8

    P+= line( [ a(value), b(value)] )

9

​

10

P.show(aspect_ratio=1, axes=False)

Language: Sage Gap GP HTML Macaulay2 Maxima Octave Python R Singular

---

Messages

In this picture, the original curve is a circle (in black), its involute is in red.

From the definition, we get these conditions:

$$\begin{array}{}\text{(1)}& \beta (s)-\alpha (s)=\lambda (s){\alpha }^{\prime }(s)\end{array}$$

for some function $\lambda (s)$, and 

$$\begin{array}{}\text{(2)}& ⟨{\beta }^{\prime }(s),{\alpha }^{\prime }(s)⟩=0.\end{array}$$

Here the primes denote differentiation with respect to $s$. Remember that $s$ is the arclength parameter for $\alpha$, but probably not for $\beta$.

If we differentiate equation (1), we find

$$\begin{array}{}\text{(3)}& {\beta }^{\prime }(s)-{\alpha }^{\prime }(s)={\lambda }^{\prime }(s){\alpha }^{\prime }(s)+\lambda (s){\alpha }^{″}(s).\end{array}$$

Now, take the inner product of both sides with ${\alpha }^{\prime }=T_{\alpha }$ and you find

$$LHS(3)=⟨{\beta }^{\prime }(s)-{\alpha }^{\prime }(s),{\alpha }^{\prime }(s)⟩=⟨\beta -T_{\alpha },T_{\alpha }⟩=-1$$

and

$$RHS(3)={\lambda }^{\prime }(s)⟨{\alpha }^{\prime }(s),{\alpha }^{\prime }(s)⟩+\lambda (s)⟨{\alpha }^{″}(s),{\alpha }^{\prime }(s)⟩={\lambda }^{\prime }⟨T,T⟩+\lambda ⟨T^{\prime },T⟩={\lambda }^{\prime }$$

If we put these together, we get the equation ${\lambda }^{\prime }(s)=1$, which integrates to 

$$\lambda =C-s$$

for some constant $C\in \mathbb{R}$. This gives us a very reasonable description of the involute $\beta$:

$$\begin{array}{}\text{(involute)}& \beta (s)=\alpha (s)+(C-s)T_{\alpha }(s).\end{array}$$

When is $\beta$ a regular curve? We compute in a straightforward manner:

$${\beta }^{\prime }(s)={\alpha }^{\prime }(s)+(C-s)T_{\alpha }^{\prime }(s)-s{T}_{\alpha }(s)=(C-s){\kappa }_{\alpha }(s)N_{\alpha }(s).$$

Thus $\beta$ is regular at those values of $s$ different from $C$ that have ${\kappa }_{\alpha }(s)\ne 0$.

Observation: 

The quantity $||\beta -\alpha ||=||(C-s)T_{\alpha }(s)||=|C-s|$ serves both as the distance between corresponding points on the curves and as a measure of arclength along $\alpha$. For this reason, involutes are sometimes called string involutes.  (And this is related to the way that Christian Huygens discovered them while trying to build a better pendulum clock.) Imagine the involute as the path taken by an endpoint of a string unwound from the original curve. This is a bit easier to imagine if you view the segment $[\alpha (s)\beta (s)]$ as describing the part of the string which has been pulled off of the curve already. For our picture above, the string segments in are in blue. Unwind the string by picking up the spot where the curves meet on the right, and then go around counter-clockwise.

Some Results about Involutes:

Theorem: An involute of a plane curve lies in the same plane.

Proof: Suppose that $\alpha$ lies in a plane with normal vector $n$. Then by an appropriate translation, we can assume that $⟨\alpha (s),n⟩=0$. Differentiating, we find that $⟨{\alpha }^{\prime },n⟩=0$. By our characterization of involutes above, we see that for an involute $\beta$, there is some constant $C$ so that $\beta =\alpha +(C-s){\alpha }^{\prime }$. Therefore,

$$⟨\beta ,n⟩=⟨\alpha +(C-s){\alpha }^{\prime },n⟩=⟨\alpha ,n⟩+(C-s)⟨{\alpha }^{\prime },n⟩=0$$.

Thus $\beta$ lies in the same plane as $\alpha$.

$◻$

Theorem: An involute of a generalized helix is a planar curve.

Proof: Suppose that the generalized helix $\alpha$ is a unit speed curve which has axis $\ell =P+tv$, where $v$ is a unit vector. By definition of generalized helix, $⟨T_{\alpha },v⟩$ is a constant. Differentiating twice, we find that $⟨N_{\alpha },v⟩=⟨N_{\alpha }^{\prime },v⟩=0$.

Let $\beta =\alpha +(C-s)T_{\alpha }$ be an involute of $\alpha$. We differentiate this twice to find
$${\beta }^{\prime }=(C-s){\kappa }_{\alpha }{N}_{\alpha },{\beta }^{″}=a{N}_{\alpha }+b{N}_{\alpha }^{\prime }$$
for some functions $a$ and $b$. As a consequence, we get that
$$⟨v,{\beta }^{\prime }⟩=0,⟨v,{\beta }^{″}⟩=0.$$
This means that $\beta$ lies in a plane perpendicular to $\ell$.

$◻$

This last theorem calls for a picture (right +Neil Martinsen-Burrell?). So let's show a string unwinding from a standard circular helix. I suggest you use FireFox to see this one, as it requires Java to run the 3D interactive visualization. Hit Evaluate to enjoy the wonders of colored lines giving a perception of space.

xxxxxxxxxx

 

1

a = 2

2

var('s')

3

x(s) = cos(s)

4

y(s) = sin(s)

5

z(s) = a*s

6

x_prime(s) = x.diff(s)

7

y_prime(s) = y.diff(s)

8

z_prime(s) = z.diff(s)

9

X(s) = x(s) - s*x_prime(s)

10

Y(s) = y(s) - s*y_prime(s)

11

Z(s) = z(s) - s*z_prime(s)

12

P = parametric_plot3d([x(s),y(s),z(s)],[s,-2*pi,2*pi],color='red',thickness=4)

13

P+= implicit_plot3d(z==0, [x,-6.3,6.3],[y,-6.3,6.3],[z,-6,6], color='yellow', opacity=.5)

14

Q = parametric_plot3d([X(s), Y(s), Z(s)], [s,0,2*pi], color='black', thickness=2)

15

​

16

@interact

17

def _(t = slider(0,2*pi,pi/10,pi)):

18

    p1 = [x(t),y(t),z(t)]

19

    p2 = [X(t),Y(t),Z(t)]

20

    R = point(p1, size=5) + point(p2, size=5) + line([p1,p2],thickness=5)

21

    (P+Q+R).show()

Language: Sage Gap GP HTML Macaulay2 Maxima Octave Python R Singular

---

Messages