Evolutes

An Evolute of a Circular Helix

In a previous post, we studied the involute of a curve. Now we will reverse the direction of the problem: Given a regular curve $\alpha$, another curve $\beta$ is called an evolute of $\alpha$ when $\alpha$ is an involute of $\beta$.

This means that

1. $\beta (s)-\alpha (s)$ lies along the tangent to $\beta$ at $\beta (s)$.
2. ${\beta }^{\prime }(s)$ is orthogonal to ${\alpha }^{\prime }(s)$.

Our goal is to understand what we can say about evolutes.

Basic Description of Evolutes

So, suppose that $\alpha$ is a unit speed curve with Frenet-Serret apparatus $\{\kappa ,\tau ,T,N,B\}$ and arclength parameter $s$. Let $\beta$ be an evolute of $\alpha$. Let's find as precise a description of $\beta$ as possible.

Begin by writing $\beta (s)=\alpha (s)+a(s)T(s)+b(s)N(s)+c(s)B(s)$. By our defining conditions, we must have that $a=0$. Thus $\beta (s)=\alpha (s)+b(s)N(s)+c(s)B(s)$. Differentiating, we find
$${\beta }^{\prime }=T+b^{\prime }N+b{N}^{\prime }+c^{\prime }B+c{B}^{\prime }=(1-b\kappa )T+(b^{\prime }-c\tau )N+(b\tau +c^{\prime })B.$$
Again, by our defining conditions, $⟨{\beta }^{\prime },T⟩=0$, so we deduce that $b=1/\kappa$, which is the radius of curvature of $\alpha$ at $\alpha (s)$. It is traditional to denote this by $\rho$, to remind us that it is a radius---in fact, it is the radius of the osculating circle. We shall do so in what follows.

We will have to work a little harder to find a reasonable form for $c$. So far, we have
$${\beta }^{\prime }=({\rho }^{\prime }-c\tau )N+(\rho \tau +c^{\prime })B.$$
Note that ${\beta }^{\prime }$ and $\beta -\alpha$ are parallel by definition of an evolute! Since these vectors lie in the normal plane $\mathrm{span}\{N,B\}$ for $\alpha$, we deduce that the ratios of the normal to the binormal coordinates must be equal. Hence,
$${\displaystyle \frac{{\rho }^{\prime }-c\tau }{\rho }}={\displaystyle \frac{\rho \tau -c^{\prime }}{c}}.$$
Some algebraic manipulation eventually yields
$$\tau ={\displaystyle \frac{c{\rho }^{\prime }-c{\rho }^{\prime }}{{\rho }^2+c^2}}=-{\displaystyle \frac{\left(c/\rho \right)}{1+{\left(c/\rho \right)}^2}}.$$
Integrating, we find that
$${\int }_{s_0}^s\tau ds+C=\mathrm{arccot}\left(c/\rho \right),$$
which means that
$$c=\rho \cot ({\int }_{s_0}^s\tau ds+C)$$
for some constant $C$.

So, there we have it. If $\beta$ is an evolute for the unit-speed curve $\alpha$, we can find some constant $C$ so that
$$\begin{array}{}\text{(evolute)}& \beta (s)=\alpha (s)+\rho (s)N(s)+\rho (s)\cot ({\int }_{s_0}^s\tau (u)du+C)B(s).\end{array}$$
This means that the set of all possible evolutes forms a one-parameter family, and for any particular choice of $C$, we can compute the corresponding evolute with the above formula.

Let's make a plot of a piece of a helix and part of an evolute for that helix. Note that I chose $C=\pi /2$ here. That makes the picture reasonable. Other choices lead to evolutes which are very long compared to the original curve, and the plot is hard to read. If you don't have Java working on your machine, add the code viewer='tachyon' inside the last set of parentheses. You will at least get a plot, though you won't be able to move it around.

xxxxxxxxxx

 

1

var('s')

2

x(s) = cos(s)/sqrt(2); y(s) = sin(s)/sqrt(2); z(s) = s/sqrt(2)

3

kappa = 1/sqrt(2); tau = 1/sqrt(2); rho = sqrt(2)

4

n1(s) = -cos(s); n2(s) = -sin(s); n3(s) = 0

5

b1(s) = sin(s)/2; b2(s) = -cos(s)/2; b3(s) = 1/2

6

alpha(s) = [x(s),y(s),z(s)]

7

P = parametric_plot3d(alpha(s), [s,-pi,pi], color='black',thickness=2, aspect_ratio=1)

8

​

9

mu(s) = rho * cot( tau*s + pi/2)

10

beta1(s) = x(s) + rho*n1(s) + mu(s)*b1(s)

11

beta2(s) = y(s) + rho*n2(s) + mu(s)*b2(s)

12

beta3(s) = z(s) + rho*n3(s) + mu(s)*b3(s)

13

beta(s) = [beta1(s), beta2(s), beta3(s)]

14

​

15

tmin = -pi/2; tmax=pi/2

16

times = [tmin,3*tmin/4, tmin/2,tmin/4, 0, tmax/4, tmax/2,3*tmax/4, tmax]

17

​

18

Q = parametric_plot3d([beta1(s),beta2(s),beta3(s)], [s,tmin,tmax],color='red',thickness=2)

19

for time in times:

20

    Q+= line([alpha(time), beta(time)], color='blue')

21

(P+Q).show()

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

---

Messages

Evolutes of Plane Curves

Consider the case where $\alpha$ is a plane curve. Then $\tau =0$, and by Equation (evolute) above, we can write
$$\begin{array}{}\text{(1)}& \beta =\alpha +\rho (N+\cot (C)B)\end{array}$$
for some choice of constant $C$. To understand this curve better, let us compute its curvature, ${\kappa }_{\alpha }$, and torsion, ${\tau }_{\beta }$. Recall that $s$ is the arclength parameter for $\alpha$, not for $\beta$, so we will have to use the formulae for general curves. In what happens, the new feature to our work will be that $B^{\prime }=-\tau N=0$.

Preparation: Derivatives of $\beta$

To get ready to use our crazy formulae, we must compute the first three derivatives of $\beta$ with respect to $s$. We start from Equation (1) and follow our noses (we are all looking in the same direction, right?), and eventually, one finds these:

$${\beta }^{\prime }={\rho }^{\prime }(N+\cot (C)B).$$

$${\beta }^{″}=-{\rho }^{\prime }\kappa T+{\rho }^{″}(N+\cot (C)B)$$

$${\beta }^{‴}=(-2{\rho }^{″}\kappa -{\rho }^{\prime }{\kappa }^{\prime })T-{\rho }^{\prime }{\kappa }^{2}N+{\rho }^{‴}(N+\cot (C)B).$$
Now, it may seem odd at first that we have chosen to express things this way, instead of sticking to our basis $\{T,N,B\}$. That normally serves us so well! But these expressions will make our work below a bit easier.

Preparation: funny expressions we need

Now it is not too hard to compute the different pieces of our crazy formulae.

$${\beta }^{\prime }\times {\beta }^{″}={\rho }^{\prime }(N+\cot (C)B)\times [-{\rho }^{\prime }\kappa T+{\rho }^{″}(N+\cot (C)B)]$$

Of course, for any vector $v$, one has $v\times v=0$, so we only need to worry about the first term on the right hand side. Using that the Frenet framing is a right handed orthonormal basis, we obtain
$${\beta }^{\prime }\times {\beta }^{″}=-({\rho }^{\prime }{)}^2\kappa (-B+\cot (C)N).$$
So now
$$⟨{\beta }^{\prime }\times {\beta }^{″},{\beta }^{\prime }\times {\beta }^{″}⟩=({\rho }^{\prime }{)}^4{\kappa }^2(1+{\cot }^2(C)).$$
Also,
$$⟨{\beta }^{\prime },{\beta }^{\prime }⟩=({\rho }^{\prime }{)}^2{(1+{\cot }^2(C))}^2$$
Finally, we need the triple product $({\beta }^{\prime }{\beta }^{″}{\beta }^{‴})$.  Since this is a determinant, the expressions above come to the rescue. A determinant is unchanged if we add a multiple of one row to another, so we get
$$({\beta }^{\prime }{\beta }^{″}{\beta }^{‴})=det\left(\begin{array}{c}{\beta }^{\prime }\\ {\beta }^{″}\\ {\beta }^{‴}\end{array}\right)=det\left(\begin{array}{c}{\rho }^{\prime }(N+\cot (C)B)\\ -{\rho }^{\prime }\kappa T+{\rho }^{″}(N+\cot (C)B)\\ (-2{\rho }^{″}\kappa -{\rho }^{\prime }{\kappa }^{\prime })T-{\rho }^{\prime }{\kappa }^{2}N+{\rho }^{‴}(N+\cot (C)B)\end{array}\right)$$
Removing multiples of the first row from the others, and then a multiple of the second row from the third, we obtain
$$({\beta }^{\prime }{\beta }^{″}{\beta }^{‴})=det\left(\begin{array}{c}{\rho }^{\prime }(N+\cot (C)B)\\ -{\rho }^{\prime }\kappa T\\ (-2{\rho }^{″}\kappa -{\rho }^{\prime }{\kappa }^{\prime })T-{\rho }^{\prime }{\kappa }^{2}N\end{array}\right)=det\left(\begin{array}{c}{\rho }^{\prime }(N+\cot (C)B)\\ -{\rho }^{\prime }\kappa T\\ -{\rho }^{\prime }{\kappa }^{2}N\end{array}\right)$$
Now, let's pull out constants and then remove a multiple of the third row from the first row to obtain
$$({\beta }^{\prime }{\beta }^{″}{\beta }^{‴})=({\rho }^{\prime }{)}^3{\kappa }^{3}det\left(\begin{array}{c}(N+\cot (C)B)\\ T\\ N\end{array}\right)=({\rho }^{\prime }{)}^3{\kappa }^{3}det\left(\begin{array}{c}\cot (C)B\\ T\\ N\end{array}\right)=({\rho }^{\prime }{)}^3{\kappa }^3\cot (C).$$

Computation of Curvature and Torsion for $\beta$:

Finally, we apply our formulae for curves which are not parametrized by arc-length. We find that the curvature is

$${\kappa }_{\beta }={\left[{\displaystyle \frac{⟨{\beta }^{\prime }\times {\beta }^{″},{\beta }^{\prime }\times {\beta }^{″}⟩}{⟨{\beta }^{\prime },{\beta }^{\prime }{⟩}^3}}\right]}^{1/2}={\left[{\displaystyle \frac{({\rho }^{\prime }{)}^4{\kappa }_{\alpha }^2(1+{\cot }^2(C))}{({\rho }^{\prime }{)}^6(1+{\cot }^2(C))}}\right]}^2={\displaystyle \frac{{\kappa }_{\alpha }}{{\rho }^{\prime }(1+{\cot }^2(C))}}.$$

Now we can compute the torsion
$${\tau }_{\beta }={\displaystyle \frac{({\beta }^{\prime }{\beta }^{″}{\beta }^{‴})}{⟨{\beta }^{\prime }\times {\beta }^{″},{\beta }^{\prime }\times {\beta }^{″}⟩}}={\displaystyle \frac{({\rho }^{\prime }{)}^3{\kappa }^3\cot (C)}{({\rho }^{\prime }{)}^4{\kappa }^4(1+{\cot }^2(C))}}={\displaystyle \frac{{\kappa }_{\alpha }\cot (C)}{{\rho }^{\prime }(1+{\cot }^2(C))}}.$$

Hey, Wait A Minute!

Those expressions have a lot in common. In particular, if we look at the ratio of curvature to torsion we see

$${\kappa }_{\beta }/{\tau }_{\beta }=\tan (C)$$

is a constant! So by Lancret's Theorem, the evolute must be a curve of constant slope, i.e. a generalized helix. And examining the proof of Lancret's theorem, we see that the angle the helix takes against its axis is $C$.

Theorem: An evolute of a plane curve is a (generalized) helix.

[[ That needs a picture: the evolute of an ellipse with a reasonably big angle would be good ]]

Special Choice

Suppose we choose $C=\pi /2$, then $\beta =\alpha +\rho N$ lies in the same plane as $\alpha$. In fact, $\beta$ is the locus of the centers of curvature of $\alpha$. We can be more clear.

Theorem: Suppose that $\beta$ is an evolute of the planar curve $\alpha$. Then $\beta$ lies in the same plane as $\alpha$ if, and only if, $\beta$ is the locus of centers of $\alpha$. Moreover, such $\beta$ is a regular curve exactly at those points where ${\kappa }_{\alpha }^{\prime }\ne 0$.

Proof: We already have one direction. Next, suppose that $\beta$ is an evolute of $\alpha$ which happens to lie in the same plane as $\alpha$. Then for some functions $a=a(s)$ and $b=b(s)$ and some constant $C$ we have 

$$\beta =\alpha +aT+bN=\alpha +\rho N+\cot (C)B.$$

Since the vector $\beta -\alpha$ has a unique expression as a linear combination of the orthonormal basis $\{T,N,B\}$, we deduce that $a=0,b=\rho ,C=\pi /2$. Therefore $\beta$ is the locus of centers of curvature of $\alpha$.

Now, is $\beta$ a regular curve? We can easily compute the speed from our work above:

$$||{\beta }^{\prime }||=⟨{\beta }^{\prime },{\beta }^{\prime }{⟩}^{1/2}={\rho }^{\prime }\cdot (1+0)={\left(1/{\kappa }_{\alpha }\right)}^{\prime }={\displaystyle \frac{-{\kappa }_{\alpha }^{\prime }}{{\kappa }_{\alpha }^2}},$$

which clearly vanishes exactly at those points where ${\kappa }_{\alpha }^{\prime }=0$.

$◻$

[[ I need a picture here, too: an ellipse with the planar evolute showing off the four vertices clearly.]]

The condition ${\kappa }^{\prime }=0$ is now obviously an interesting one. Such a point is an example of a vertex of a plane curve. We might have more to say about those later. For now, the vertices of the plane curve $\alpha$ are those points which correspond to singular points on its coplanar evolute, its locus of centers of curvature.

There is a paper about evolutes in the issue of the American Mathematical Monthly that just appeared, so this serves as a small introduction. Perhaps we'll have more to say about the topic from a more geometric vein later.