The Triskele Curve

loopspace

2nd August 2018

Creative Commons License

Contents

  1. Home

  2. 1. Introduction

  3. 2. Defining the Curve in 3D

  4. 3. Defining the Curve in 2D

  5. 4. Special Case: Ellipse

  6. 5. Finding Functions

1 Introduction

I recently came across Triskele Globes. It was sufficiently near to Christmas that I thought I'd have a go at making one as they look quite festive. I was right: they do look quite festive. However, looking very closely at them then I wasn't quite sure that the fold was along the right curve. There seemed to be a couple of points of deformation, particularly at the corners. So I thought I'd have a look at figuring out what the curve ought to be.

2 Defining the Curve in 3D

Let's generalise. Suppose we wrinkle a piece of paper along its length so that its cross section is the same all through.

Let's assume that we start with a piece of paper lying in the positive quadrant of the x-y–plane up against the axes. We then wrinkle it in the y–direction, meaning that as you travel in the y–direction then the paper changes height but along the x–direction then it stays constant. We keep it above the positive quadrant and so that its edges stay above the axes (though we give ourselves leave to move it vertically).

This means that the paper is described by a parametrisation of the form (s,f(t),g(t)).

We want to fold it straight down so that the folded part has constant cross section in the z–direction. Once it has been folded, the bottom edge of the paper will be contained in a horizontal plane since we can only flex the paper in one direction at a time. We can therefore arrange it so that the bottom edge (once folded) lies in the x-y–plane.

Pre-folding, this edge was in the y-z–plane and post-folding it is in the x-y–plane. Therefore the distance along the paper from the fold to the y-z–plane prior to folding is the same as the distance along the paper from the fold to the x-y–plane after folding. That is to say, the fold lies on the plane x=z.

On the paper, that translates to s=g(t). Thus the fold curve is (g(t),f(t),g(t)).

After folding, this then continues down to the x-y–plane and the parametrisation of the folded part is therefore (g(t),f(t),s). A full parametrisation of the paper is therefore:

(max(s,g(t)),f(t),min(s,g(t)))

3 Defining the Curve in 2D

The previous section defines the curve in space, but what we want to do is to draw it on the paper at the outset. To do that, we need to draw it on the paper and then "unwrinkle" the paper. If we let (u,v) be the paper's natural coordinates then in the pre-folded parametrisation, we have u=s and v is the arc length along (f(t),g(t)). That is,

v(t)=0tf'(r)2+g'(r)2dr

Hence the curve (g(t),f(t),g(t)) unwrinkles to the curve:

(g(t),0tf'(r)2+g'(r)2dr)
1

When the cross-section is given by a height function – that is, f(t)=t – the fold curve is (g(t),t,g(t)) in 3D space and on the paper is given by:

(g(t),0t1+g'(r)2dr)
2

4 Special Case: Ellipse

When making the Triskele Globes, each piece of paper has the following cross section prior to the folding.

This consists of two straight parts (which form the inside of the globe) and two elliptical ends (which will get folded down). The part we are most interested in is the elliptical end.

If we scale our measurements so that the inner cube has length 2, then the parametrisation of the ellipse is:

(sin(ϕ),bcos(ϕ)),-π/2ϕπ/2

with 0<b1. The parametrisation of the paper is therefore:

(s,sin(ϕ),bcos(ϕ))

To get the fold curve, we need to figure out the curve:

(bcos(ϕ),-π/2ϕcos2(ψ)+bsin2(ψ)dψ)

The integral simplifies (slightly) to:

-π/2ϕcos2(ψ)+b2sin2(ψ)dψ=0ϕ1-(1-b2)sin2(ψ)dψ=-π/2ϕ1-msin2(ψ)dψ

where m2=1-b2.

This is an incomplete elliptic integral of the 2nd kind. Unless b=1, it has no closed form solution. Fortunately, there are numerical solutions and the scipy module for python can compute them. The program curve.py uses this to produce a table of coordinates for the paper curve which can be read in by the LaTeX document triskele.tex to produce the template; the version for b=.8 is triskele.pdf.

Below this is plotted with the approximation of an arc for comparison. It should be noted that they are very close.

If b=1 (whence the cross section is circular) then the integral simplifies to

-π/2ϕcos2(ψ)+sin2(ψ)dψ=-π/2ϕdψ=ϕ+π/2

and so the fold curve has equation

(cos(ϕ),ϕ+π/2)

which we can plot directly to obtain the following curve, overlaid again with an arc so show how the separation has increased.

Note that the case where b=1 means that the ends are circular and that the two folds meet at the centre.

5 Finding Functions

In Equation 1 there are three functions: f(t), g(t), and the arc length function, say h(t). These are related via a Pythagorean relationship on the derivatives that states that:

h'(t)2=f'(t)2+g'(t)2
3

Thus the problem of finding such triples (f,g,h) which have a closed form expression (for some suitable definition thereof) is closely related to a functional version of finding Pythagorean triples.

The fact that we are working with derivatives is what makes this a non-trivial question. As we have already seen, this problem is complicated even for an ellipse.

Trigonometric and hyperbolic functions form a natural source of these triples:

f(t) g(t) h(t)
sin(t) cos(t) t
ln|cos(t)| t ln|sec(t)+tan(t)|
cosh(t) t sinh(t)
tan-1(sinh(t)) ln(cosh(t)) t
11+x 11-x sin-1(t)