Oloid

Download as pdf or txt
Download as pdf or txt
You are on page 1of 14

Journal for Geometry and Graphics

Volume 1 (1997), No. 2, 105–118

The Development of the Oloid

Hans Dirnböck∗ and Hellmuth Stachel×


∗ Institute
of Mathematics, University Klagenfurt
Universitätsstr. 65, A-9020 Klagenfurt, Austria

× Institute
of Geometry, Vienna University of Technology
Wiedner Hauptstr. 8-10/113, A-1040 Wien, Austria
email: [email protected]

Abstract. Let two unit circles kA , kB in perpendicular planes be given such that
each circle contains the center of the other. Then the convex hull of these circles is
called Oloid. In the following some geometric properties of the Oloid are treated
analytically. It is proved that the development of the bounding torse Ψ leads to
elementary functions only. Therefore it is possible to express the rolling of the
Oloid on a fixed tangent plane τ explicitly. Under this staggering motion, which is
related to the well-known spatial Turbula-motion, also an ellipsoid Φ of revolution
inscribed in the Oloid is rolling on τ . We give parameter equations of the curve
of contact in τ as well as of its counterpart on Φ.
The surface area of the Oloid is proved to equal the area of the unit sphere. Also
the volume of the Oloid is computed.
Keywords: Oloid, Turbula-motion, development of torses.
MSC 1994: 51N05, 53A05, 53A17

1. Introduction
Let kA , kB be two unit circles in perpendicular planes Π1 , Π2 such that kA passes through
the center MB of kB and kB passes through the center MA of kA (see Fig. 1)1 , The torse
(developable) Ψ connecting kA and kB is the enveloping surface of all planes τ that touch kA
and kB simultaneously. If any tangent plane τ contacts kA at A and kB at B, then the line
AB is a generator of Ψ. In this case the tangent line of kA at A must intersect the tangent
line of kB at B in a finite or infinite point T on the line 12 of intersection between Π1 and Π2
(see Fig. 2; the triangle ABT can also be found in Fig. 5 and Fig. 6).
1
All figures in this paper are orthogonal views. But only in Fig. 5 and Fig. 6 the superscript “n” is used
to indicate that geometric objects have been projected orthogonally into a plane.
c 1997 Heldermann Verlag
ISSN 1433-8157/$ 2.50 °
106 H. Dirnböck, H. Stachel: The Development of the Oloid

Figure 1: Circles kA , kB defining the Oloid

We choose the planes Π1 , Π2 as coordinate planes and the midpoint O of MA MB as the


origin of a cartesian coordinate system. Then we may set up the equations of kA , kB as
kA : x2 + (y + 12 )2 = 1 and z=0
(1)
kB : (y − 21 )2 + z 2 = 1 and x = 0.
We parametrize the torse Ψ by the arc-length t of kA with the starting point t = 0 at U on

Figure 2: Coordinate system and notation

the negative y-axis. Then we obtain the coordinates


µ ¶
1
A = sin t , − − cos t , 0 . (2)
2
H. Dirnböck, H. Stachel: The Development of the Oloid 107
Since the point T on the y-axis is conjugate to A with respect to kA , we get
µ ¶
2 + cos t
T = 0, − , 0 . (3)
2 cos t

In the same way conjugacy between T and B with respect to kB implies


µ √ ¶
1 cos t 1 + 2 cos t
B = 0, − , ± . (4)
2 1 + cos t 1 + cos t

The upper sign of the z-coordinate corresponds to the upper half of Ψ.2
From (2) and (4) we compute the squared length of the line segment AB as
µ ¶2
2 2 cos t 1 + 2 cos t
AB = sin t + 1 + cos t − + = sin2 t + (1 + cos t)2 − 2 cos t + 1 ,
1 + cos t (1 + cos t)2

which results in
Theorem 1: All line segments AB of the torse Ψ are of equal length

AB = 3 . (5)

This surprising result has already been proved in [7]. But probably also P. Schatz was aware
of this result when he took out a patent for the Oloid (cf. [8]) in 1933 (see also [9], Figures
155, 156 and p. 122).

Let u denote the arc-length of kB , starting on the positive y-axis. Then A ∈ kA and
B ∈ kB are points of the same generator of Ψ if and only if the parameters t of A and u of B
obey the involutive relation
cos t t u 1
cos u = − or cos2 cos2 = . (6)
1 + cos t 2 2 4
For real generators of Ψ the condition 1 + 2 cos t ≥ 0 is necessary. By the restriction
2π 2π 2π 2π
− <t< and − <u< (7)
3 3 3 3
we avoid vanishing denominators. It has to be noted that for Ψ the parametrization by t
becomes singular at t = ±2π/3.
In the following we restrict each generator of Ψ to the line segment AB. Thus we obtain
just the boundary of the convex hull of kA and kB .

2. Development of the Torse Ψ


When Ψ is developed into a plane τ , then the circles kA , kB are isometrically transformed into
planar curves kAd , kBd , respectively. It is well-known from Differential Geometry (see e.g. [11],
p. 209 or [12], p. 72) that at corresponding points A ∈ kA ⊂ Ψ and Ad ∈ kAd ⊂ τ the geodesic
curvatures are equal. This can be expressed in a more geometric way as follows (cf. [4], p. 295):
2
In the generalization presented in [5] the circles kA , kB are replaced by congruent ellipses with a common
axis.
108 H. Dirnböck, H. Stachel: The Development of the Oloid

When τ is specified as the tangent plane of Ψ along the generator AB, then the curvature
center K of kAd at Ad = A is located on the curvature axis of kA at A, which is the axis of
revolution of (the curvature circle) kA (see Fig. 2, compare Fig. 3). Since K = (− 21 , 0, ±k) 3
is aligned with T and B, we get for the squared curvature radius
2 2 + 2 cos t
ρ2 = AK = 1 + k 2 = .
1 + 2 cos t
Hence the curvature κ of kAd reads
s
1 1 + 2 cos t
= κ(t) = . (8)
ρ 2(1 + cos t)

This is the so-called natural equation of kAd with arc-length t .4


In order to deduce an explicit representation of kAd , we choose τ as the tangent plane at
the point U ∈ kA with minimal y-coordinate. In τ we introduce a cartesian coordinate system
with origin U d = U and axes I and II (see Fig. 3). We define the first coordinate-axis I parallel
to the tangent vector of kA at U . Then due to Riccati’s formula (see e.g. [10], p. 44) we get
Z t
IA (t) = I0 + cos α(t) dt Z t
0
Z t for α(t) := α0 + κ(t) dt (9)
0
IIA (t) = II0 + sin α(t) dt
0

with the specifications α0 = I0 = II0 = 0 . By integration of (8) we obtain


√ √
6 sin t 3 tan 2t
α(t) = 2 arcsin √ − arcsin . (10)
3 1 + cos t 3

Theorem 2: In the cartesian coordinate system (I, II ) (see Fig. 4) the arc length parame-
trization of the development kAd of the circle kA reads
√ " √ #
2 3 p 2 cos t
IA (t) = 2(1 + 2 cos t)(1 − cos t) + arccos √
9 1 + cos t
√ · ¸ (11)
3 2
IIA (t) = 4(1 − cos t) + ln .
9 1 + cos t

Proof: The integrals in the left column of (9) could not be immediately solved with the use
of common computer-algebra-systems. We succeeded as follows: The integral for IIA (t) can
be transformed into
Z t à √ √ !
6 sin t 3 tan 2t
IIA (t) = sin 2 arcsin √ − arcsin dt =
0 3 1 + cos t 3
Z t" √ √ #
4 6 sin 2t (1 + 2 cos t) 3 tan 2t (4 cos t − 1)
= √ − dt ,
0 9 1 + cos t 9
3
The sign of the z-coordinate √ is equal to that of B in (4).
4 1
Note ρ̇(0) = 0 , but ρ̈(0) = 18 3 6= 0. This proves that at U d there is exactly a four-point contact between
d
kA and its curvature circle (see Fig. 4 or Fig. 5).
H. Dirnböck, H. Stachel: The Development of the Oloid
Figure 3: Axonometric view of the Oloid and its development into τ

109
110 H. Dirnböck, H. Stachel: The Development of the Oloid

and this gives rise to the second equation in (11). From


√ µ ¶
dIIA 3 1
= sin t 4 + = sin α (12)
dt 9 1 + cos t

due to (9) we obtain


s µ ¶2 √ 3
dIA dIIA 6(1 + 2 cos t) 2
= cos α = 1− = √ . (13)
dt dt 9 1 + cos t

Then the integration can be carried out using the substitution t := tan 2t . The first quarter
of the developed curve kAd ends at
à √ √ !
¡ ¡ 2π ¢ ¡ 2π ¢¢ 2π 3 2 3
IA 3 , IIA 3 = , (3 + ln 2 ) ≈ (1.2092 , 1.4215 ).
9 9

There is an analogous representation of the developed image kBd of the circle kB in terms
of its arc-length u. The curves kAd and kBd are congruent since halfturns about the axes
x ± z = y = 0 interchange kA and kB while the Oloid is transformed into itself. However,
based on (11) and due to (5) the curve kBd can also be parametrized in the form

IB (t) = IA (t) + √3 cos eAB ,
(14)
IIB (t) = IIA (t) + 3 sin eAB .

Here the angle eAB = α + γ (see Fig. 4) defines the direction of the developed generator Ad B d .
Angle α has already been computed in (12) and (13). γ is the angle made by the generator
AB of Ψ and the tangent vector

vA = (cos t, sin t, 0) (15)


−→
of kA at A. The dot product of vA and the vector AB according to (2) and (4) gives
√ −→ sin t cos t sin t
3 cos γ = vA · AB = − sin t cos t + sin t + sin t cos t − = .
1 + cos t 1 + cos t
Elementary trigonometry leads to
s
2(1 + 2 cos t)
sin γ = (16)
3(1 + cos t)

and finally to

7 + 7 cos t + 4 cos2 t
sin eAB =
9(1 + cos t)
√ p (17)
2 2(2 + cos t) (1 − cos t)(1 + 2 cos t)
cos eAB = − .
9(1 + cos t)

We substitute these formulas in (14). Then due to (11) we obtain


H. Dirnböck, H. Stachel: The Development of the Oloid
Figure 4: The development of the Oloid with the images kAd of kA and kBd of kB together with the evolute ekAd of kAd

111
112 H. Dirnböck, H. Stachel: The Development of the Oloid

n
Figure 5: Detail of Fig. 4 with the image kO of the center curve kO
under orthogonal projection into τ

Theorem 3: In the cartesian coordinate system (I, II ) of τ (see Fig. 4 or Fig. 3) the devel-
opment kBd of the circle kB has the parametrization with respect to the arc-length t of kA as
follows:
√ " √ p #
2 3 2 cos t 2(1 − cos t)(1 + 2 cos t)
IB (t) = arccos √ −
9 1 + cos t (1 + cos t)
√ · ¸ (18)
3 2 11 + 7 cos t
IIB (t) = ln + .
9 1 + cos t 1 + cos t

In a similar way also the evolute ekAd of kAd (see Fig. 4 or Fig. 5) can be computed. The
parameter representation
IK (t) = IA (t) − ρ sin α
IIK (t) = IIA (t) + ρ cos α
H. Dirnböck, H. Stachel: The Development of the Oloid 113
of ekAd makes use of the curvature radius ρ according to (8). From (12) and (13) we obtain
p
(5 + 4 cos t) 6(1 − cos t)
ρ sin α = √
√ 9 1 + 2 cos t
2 3
ρ cos α = (1 + 2 cos t)
9
and finally as parametrization of the evolute ekAd of kAd
√ √ p
2 3 2 cos t 2(1 − cos t)
IK (t) = arccos √ − p
9 1 + cos t 3(1 + 2 cos t)
√ · ¸ (19)
3 2
IIK (t) = 6 + ln .
9 1 + cos t

The evolute ekAd obviously (see Fig. 4) does not pass through the cusps of kAd ; the curvature
radius ρ tends to infinity. This reveals that these cuspidal points are not ordinary. For kBd
the Taylor-series expansion of the parameter representation (18) at t = 0 is

1 5 7
√ 3 2
IB (t) = t + O(t ), IIB (t) = 3 + t + O(t4 ).
360 12
Therefore the singularities of kAd and kBd are of order 2 and class 3 (German: Rückkehrflach-
punkte).

3. Motions Related to the Oloid


According to Fig. 3 we assume that the Oloid is rolling on the upper side of τ . In the following
we therefore choose for point B in (4) the negative z-coordinate. For the sake of brevity we
substitute

s := sin t and c := cos t with − 1


2
< c ≤ 1, −1 ≤ s ≤ 1 , s 2 + c2 = 1 . (20)

In order to describe the rolling of Ψ on τ we introduce a moving frame of Ψ with origin


A ∈ kA . The first vector of this frame is the tangent vector vA according to (15). The second
vector wA perpendicular to vA is specified in the tangent plane τ . We define
µ ¶ ³ √ ´
1 1 −→ 1 √
wA := √ AB − vA cos γ = p −s 1 + 2c , c 1 + 2c , −1 .
sin γ 3 2(1 + c) (21)

The vector
1 ³ √ ´
nA := vA × wA = p −s , c , 1 + 2c (22)
2(1 + c)

perpendicular to τ completes this cartesian frame. nA is pointing to the interior of Ψ .


While the Oloid is rolling on the fixed plane τ , the frame (A; vA , wA , nA ) shall be moving
along Ψ in such a way, that A is the running point of contact between kA and τ . This
implies that A ∈ kA is always coincident with the corresponding point Ad ∈ kAd . Therefore
114 H. Dirnböck, H. Stachel: The Development of the Oloid

the elements of the moving frame get the following coordinates with respect to the cartesian
coordinate system (U d ; I, II, III ) attached to τ :

A = (IA (t), IIA (t), 0) , vA = (cos α, sin α, 0) ,


(23)
wA = (− sin α, cos α, 0) , nA = (0, 0, 1) .

Let (x, y, z) denote the coordinates of any point P attached to the Oloid. The required
representation of the motion consists of a matrix equation which allows to compute the in-
stantaneous coordinates (I, II, III ) of point P with respect to the fixed plane τ , in dependence
of the motion parameter t . In order to obtain this equation we firstly compute the coordinates
(ξ, η, ζ) of P with respect to the moving frame (A; vA , wA , nA ). Though P is attached to Ψ,
these coordinates are dependent on t. From (2), (15), (21) and (22) we get
     p √  
x s c p2(1 + c) −s 1 + 2c −s ξ

 y  =  −1 − c  + p 1  s 2(1 + c) c 1 + 2c c   η .
2
2(1 + c) √
z 0 0 −1 1 + 2c ζ (24)

Secondly, according to (23) the motion of the moving frame with respect to τ (see Fig. 3)
reads       
I IA (t) cos α − sin α 0 ξ
 II  =  IIA (t)  +  sin α cos α 0   η .
III 0 0 0 1 ζ
Now we eliminate (ξ, η, ζ) from these two matrix equations with orthogonal 3 × 3-matrices.
After substituting (11), (12) and (13) we obtain by straight-forward calculation
Theorem 4: Based on the cartesian coordinate systems (x, y, z) in the moving space and
(I, II, III ) in the fixed space, the rolling of the Oloid on the tangent plane τ can be represented
as
 √ √ 
cs 1 + 2c c 2
 2(1 + c)p2(1 + c) + 2 arccos √1 + c 
     
I √   ³ ´ x
 2 
 II  = 3  15 + 13c − c
+ ln
2  + aij  y  , where
9  2(1 + c) 1 + c 
III √  z
 
 3 3 (2 + c) 
p
2 2(1 + c)
 √ √ 
(5 + c) 1 + 2c (2 + c)s 1 + 2c (5 + 4c)s
 p2(1 + c) p p
(1 + c) 2(1 + c) (1 + c) 2(1 + c) 
 
³ ´ √3   (c − 1)s 5 + 5c − c 2
(1 + 2c)

1 + 2c


aij =  − .
9  1+c 1+c 1+c 

 √ √ p 
 3s 3 3c 3 3 3(1 + 2c) 
−p p p
2(1 + c) 2(1 + c) 2(1 + c)
Here c and s stand for cos t and sin t, respectively, while the motion parameter t obeys (7).

The first vector on the right side of this matrix equation represents the path kO of the
Oloid’s center O under this rolling motion. In particular, the third coordinate of this vector
H. Dirnböck, H. Stachel: The Development of the Oloid 115
gives the oriented distance
2+c
r := p (25)
2 2(1 + c)
between O and the tangent plane for each t. Due to the introduced moving frame, this
−→
distance r equals the dot product nA · AO . One can verify that for each t the velocity vector
of the center curve kO is perpendicular to the axis Ad B d of the instantaneous rotation5 (see
n
orthogonal view kO of kO in Fig. 5 or Fig. 6).
The rolling of the Oloid is truly staggering. It is related to the Turbula motion (see [13] or
[7]6 and the references there) which is used for shaking liquids. It turns out that the Turbula
motion is inverse to the motion of the moving frame in (24).
The circles kA and kB can also be seen as singular surfaces b kA , b
kB of 2nd class. The
coordinates (u0 : u1 : u2 : u3 ) of their tangent planes
u0 + u 1 x + u 2 y + u 3 z = 0
match the “tangential equations”
b
kA : 4u20 − 4u0 u2 − 4u21 − 3u22 = 0 , b
kB : 4u20 + 4u0 u2 − 3u22 − 4u23 = 0 .
Then due to a standard theorem of Projective Geometry the torse Ψ is not only tangent to
kA and kB but to all surfaces of 2nd class included in the range which is spanned by b
kA and
b 7
kB . Among these surfaces there is an ellipsoid Φ of revolution obeying the equation
1 b
b:
Φ : 6x2 + 4y 2 + 6z 2 = 3 or Φ 2
(k A + b
kB ) = 4u20 − 2u21 − 3u22 − 2u23 = 0 (26)

3
with focal points MA , MB and semi-axes 2
and √12 (cf.
[13], p. 31). The curve l of contact
between Φ and the torse Ψ is located on cylinders which are the images of kA and kB ,
respectively, in the polarity with respect to Φ. Therefore this curve has the representations
¡ ¢2 ¡ ¢2
l : 3x2 + y − 21 = 3z 2 + y + 21 = 1 or
√ (27)
s 3c ± 1 + 2c
x= , y=− , z= .
2+c 2(2 + c) 2+c
Together with the Oloid also the inscribed ellipsoid Φ is rolling on τ . In the fixed plane τ the
point of contact with the rolling ellipsoid traces a curve l d . The parameter representation
√ √ √ · ¸
d 2 3 c 2 3 2 3(5 + c)
l : I= arccos √ , II = ln + , III = 0 (28)
9 1+c 9 1+c 2+c
of this isometric image of l ⊂ Ψ is obtained by transforming the coordinates of l given in (27)
(negative sign) under the matrix equation of Theorem 4 .
Fig. 6 shows not only the fixed tangent plane τ with the developed curves kAd , kBd and ld
in true shape. In this figure also an orthogonal view of the Oloid with the inscribed ellipsoid
Φ and the curve l of tangency is displayed.
5
In general the instantaneous motion is a helical motion. However when a torse is rolling on a plane, the
helical parameter must vanish (cf. [3], p. 161 or [6]).
6
In this paper a very particular plane-symmetric six-bar loop is studied which is also displayed in [1],
Figure 1. In each position of this loop and for each two opposite links Σ, Σ0 there is a plane τ of symmetry.
It turns out that relatively to Σ these planes τ are tangent to a torse of type Ψ. The Turbula motion is the
motion of Σ relative to τ , when in τ the generator of the torse is kept fixed.
7
In the cases treated in [5] kA and kB are ellipses, but Φ is a sphere. This implies that the center O of
gravity has a constant distance to τ during the rolling motion.
116
H. Dirnböck, H. Stachel: The Development of the Oloid
Figure 6: The Oloid and the inscribed ellipsoid Φ are rolling on τ while the curve l of contact between Ψ and Φ traces l d
H. Dirnböck, H. Stachel: The Development of the Oloid 117
4. Surface Area and Volume of the Oloid
As the development of a torse is locally an isometry, the computation of the area of Ψ can
be carried out either in the 3-space or after the development into the plane τ . We prefer the
latter and use a formula given in [2], p. 118, eq. (5): The area swept out by the line segment
AB under a planar motion for t0 ≤ t ≤ t1 can be computed according to
Z t1 ° °
°1 −→ °
S= ° 2 (vA + vB ) × AB ° dt , (29)
t0

where vA , vB are the velocity vectors of the endpoints. For vectors in R2 the norm in this
formula can be cancelled which gives rise to an even oriented area.
In the coordinate system (I, II ) of τ we obtain due to (14) and (11)
µ ¶
−→ ¡√ √ ¢ dIA dIIA
AB = 3 cos eAB , 3 sin eAB , vA = , ,
dt dt
µ ¶
dIA √ d eAB dII √ d eAB
vB = − 3 sin eAB , + 3 cos eAB
dt dt dt dt
and according to (12) and (13)
µ ¶
dS √ dIA dIIA 3 d eAB √ 3 d eAB
= 3 sin eAB − cos eAB − = 3 sin γ − .
dt dt dt 2 dt 2 dt
Eq. (16) and the derivation of (17) lead to
p √
dS 2(1 + 2 cos t) 3 2 cos t
= √ − p
dt 1 + cos t 2 (1 + cos t)(1 + 2 cos t)
which finally results in
dS 2 + cos t
=p . (30)
dt 2(1 + cos t)(1 + 2 cos t)
After integration we obtain up to a constant k
· ¸
1 1 − 4 cos t 1 + 5 cos t
S(t) = arcsin − arcsin +k
2 3 3(1 + cos t)
and for the complete torse Ψ
h ³π ´ i · µ ¶ ¸

SΨ = 8 S − S(0) = 4 S − S(0) = 4π . (31)
2 3
Theorem 5: The surface area of the Oloid equals that of the unit sphere.
The computation of the Oloid’s volume starts from (30): Each surface element of Ψ is
the base of a volume element forming a pyramid with apex O. Its altitude r has already been
computed in (25) as it equals the distance between O and the corresponding tangent plane.
Thus we obtain
r (2 + cos t)2
dV = dS = √ dt . (32)
3 12(1 + cos t) 1 + 2 cos t
A numerical integration gives
h ³π ´ i
VΨ = 8 V − V (0) ≈ 3.05241 . (33)
2
118 H. Dirnböck, H. Stachel: The Development of the Oloid

Acknowledgements
The authors express their gratitude to Christine Nowak (Klagenfurt) for her successful help
in solving the integrals. Sincere thanks also for permanent contributions given by Jakob
Kofler (Klagenfurt). The authors are indebted to Gerhard Hainscho (Wolfsberg) for
fruitful discussions and to Sabine Wildberger (Klagenfurt) for manufacturing a model.
Finally the authors thank Manfred Husty (Leoben) for useful comments and Gunter Weiss
(Dresden) for pointing their attention to the singularities of kAd and kBd .

References
[1] J.E. Baker: The Single Screw Reciprocal to the General Plane-Symmetric Six-Screw
Linkage. JGG 1, 5–12 (1997).
[2] W. Blaschke, H.R. Müller: Ebene Kinematik. Verlag von R. Oldenbourg, München
1956.
[3] O. Bottema, B. Roth: Theoretical Kinematics. North-Holland Publ. Comp., Ams-
terdam 1979.
[4] H. Brauner: Lehrbuch der Konstruktiven Geometrie. Springer-Verlag, Wien 1986.
[5] C. Engelhardt, C. Ucke: Zwei-Scheiben-Roller. MNU, Math. Naturw. Unterr. 48/5,
259–263 (1995).
[6] M. Husty, A. Karger, H. Sachs, W. Steinhilper: Kinematik und Robotik.
Springer, Berlin 1997.
[7] S. Kunze, H. Stachel: Über ein sechsgliedriges räumliches Getriebe. Elem. Math. 29,
25–32 (1974).
[8] P. Schatz: Deutsches Reichspatent Nr. 589 452 (1933) in der allgemeinen Getriebe-
klasse.
[9] P. Schatz: Rhythmusforschung und Technik. Verlag Freies Geistesleben, Stuttgart 1975.
[10] K. Strubecker: Differentialgeometrie I. 2. Aufl., Sammlung Göschen, Bd. 1113/1113a,
Walter de Gruyter & Co, Berlin 1964.
[11] K. Strubecker: Differentialgeometrie II. 2. Aufl., Sammlung Göschen, Bd. 1179/
1179a, Walter de Gruyter & Co, Berlin 1969.
[12] T.J. Willmore: An Introduction to Differential Geometry. Oxford University Press
1969.
[13] W. Wunderlich: Umwendung einer regelmäßigen sechsgliedrigen Würfelkette. Pro-
ceedings IFToMM Symposium Mostar, May 1980, 23–33.

Received July 31, 1997; final form November 12, 1997

You might also like