- About this Journal ·
- Abstracting and Indexing ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Recently Accepted Articles ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents

Journal of Applied Mathematics

Volume 2012 (2012), Article ID 268537, 34 pages

http://dx.doi.org/10.1155/2012/268537

## Towards a Prototype of a Spherical Tippe Top

^{1}Howest, ELIT, University College West Flanders, G. K. De Goedelaan 5, 8500 Kortrijk, Belgium^{2}Department of Mathematical Analysis, Research Group NaM2, University of Ghent, Galglaan 2, 9000 Ghent, Belgium^{3}Department of Architecture, Sint-Lucas Visual Arts, Institute for Higher Education in the Sciences and the Arts, 9000 Ghent, Belgium^{4}Howest, Industrial Design Center, University College West Flanders, Marksesteenweg 58, 8500 Kortrijk, Belgium

Received 14 April 2011; Revised 6 October 2011; Accepted 7 October 2011

Academic Editor: Yuri Sotskov

Copyright © 2012 M. C. Ciocci et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

Among spinning objects, the tippe top exhibits one of the most bizarre and counterintuitive behaviours. The commercially available tippe tops basically consist of a section of a sphere with a rod. After spinning on its rounded body, the top flips over and continues spinning on the stem. The commonly used simplified mathematical model for the tippe top is a sphere whose mass distribution is axially but not spherically symmetric, spinning on a flat surface subject to a small friction force that is due to sliding. Three main different dynamical behaviours are distinguished: *tipping, nontipping, hanging*, that is, the top rises but converges to an intermediate state instead of rising all the way to the vertical state. Subclasses according to the stability of relative equilibria can further be distinguished. Our concern is the degree of confidence in the mathematical model predictions, we applied 3D printing and rapid prototyping to manufacture a “3-in-1 toy” that could catch the three main characteristics defining the three main groups in the classification of spherical tippe tops as mentioned above. We propose three designs. This “*toy*” is suitable to validate the mathematical model qualitatively and quantitatively.

#### 1. Introduction

Spinning toys are among the most ancient toys, and there is a great variety of them. It is quite simple to start spinning objects like a top or a gyroscope, and though it is simple to explain their motion in general, it is challenging to write down the detailed equations of motion. Among spinning objects, the *tippe top* exhibits one of the most bizarre and counterintuitive behaviours. The commercially available tippe tops, patented in Denmark in the 50s, basically consist of a section of a sphere with a rod. After spinning on its rounded body, the top flips over and continues spinning on the stem. It is the friction with the bottom surface and the position of the center of mass below the centre of curvature that cause the tippe top to rise its centre of mass while continuing rotating around its symmetry axis (through the stem). See Figure 1 for an illustration. Remarkably, at the inverted state, the center of mass lies higher than at the initial condition, defying gravity. Experimentally, it is known that such a transition occurs only if the initial spin exceeds a certain critical threshold.

The commonly used simplified mathematical model for the tippe top is a sphere whose mass distribution is axially but not spherically symmetric, spinning on a flat surface subject to a small friction force that is due to sliding. Adopting a bifurcation theory point of view we reach a global geometric understanding of the phase diagram of this dynamical system. According to the eccentricity of the sphere and the Jellett invariant (which includes information on the initial angular velocity) six main classes of tops can be identified within three main groups according to the distinguished dynamical behaviours: *tipping*, *nontipping*, and *hanging*; see Figure 3. Note that objects displaying inversion properties such as the tippe top have been known since the 1800s, see, for example, [1]. After the type of tippe top as in Figure 1 was introduced in Denmark, several theoretical articles have been published since then, see, for example, [2–6] for a survey of the literature. Since it was established that sliding friction was necessary to explain the tippe top inversion [2, 7, 8], many studies have been dedicated to the analysis of models for tippe tops, involving linear stability analysis of the relative equilibria, numerical simulations, and so forth. Some studies have addressed the occurrence of transitions between rolling and sliding during the motion, see [3, 4, 9]. In this paper the presented mathematical results mainly reproduce those in [10–15] but our approach is inspired by the hands-on numerical approach as first attempted by Cohen in [2]. We believe this approach is the best choice in giving a clear view of the role of the different parameters that is necessary during the design process of an actual three-dimensional object that effectively demonstrates the model. We remark that in the mathematical model we stick to the common assumption that the only external force acting on the system consists of a normal reaction force and a frictional force of viscous type opposing the motion of the contact point in the supporting plane. This is the most common assumption in the literature, though in [4] the inclusion of a nonlinear Coulomb-type friction is discussed. It is shown there by numerical simulations that the Coulomb term contributes to the tippe top inversion but the effect is weaker compared to the viscous term. The nonlinear Coulomb term results in algebraic destabilization of the initially spinning top, whereas the viscous friction gives exponential destabilization, see also [5]. This argument motivates our choice of including viscous friction only.

The phase diagram and bifurcation diagrams illustrate the main results that confirm the findings described in [10]; the type of asymptotic dynamics is a function of the Jellett invariant (which includes information on the initial angular velocity) and eccentricity of the sphere. Either the asymptotic state is unique or the system is bistable. Three main different regimes are distinguished: no tippe top phenomenon occurs no matter what the initial spin is, tippe top dynamics may occur if the Jellett invariant (which is proportional to the initial spin) is sufficiently large, or *incomplete* tippe top behaviour occurs, where the top rises but converges to an intermediate state instead of rising all the way to the vertical state.

We underline that though the classification results can be obtained in a less cumbersome way by using the Routhian reduction as in [10], the approach used in this paper is standard and straightforward to implement from a prototyping point of view. Also, it is amenable for extensions to include for example transitions from sliding to rolling and vice versa. Our concern in this paper is the degree of confidence in the mathematical model predictions. We wanted to be sure that the mathematical model as presented here and in [10] reflected the reality. Our goal was to investigate if it was possible to make a “3-in-1 toy” that could catch the three main characteristics *tipping, nontipping, *and* hanging* that define the three main groups in the classification of spherical tippe tops as mentioned above. As far as we know such a toy does not exist yet. We successfully applied the methodology for efficient use of prototyping during the design process as presented in [16]. To the best of our knowledge this is the first time that 3D printing and rapid prototyping is being applied to design and to produce a “toy” suitable to validate the qualitative and quantitative mathematical model describing the behaviour of a dissipative nonlinear dynamical system. From the bifurcation diagram it was clear that it should be possible to hit three out of the 6 classes of tippe tops (one type for each main group) by keeping one of the characterizing parameters of the system fixed and varying the other. We believe that the realization of an actual toy is a more powerful validation of the model than software simulations which are directly affected by the underlying mathematical idealization assumptions.

Since the two parameters on which the whole classification is based, *inertia ratio* and *eccentricity*, are not independent, the challenge was to come up with a feasible prototype which could be easily mechanically driven. After detailed mathematical calculations and the development of 3D animations (see electronic attachments or http://cage.ugent.be/~bm/tippetop/tippetop.html) we used 3D printing to create a functional model giving us a quick and easy hands-on demonstration capability.

#### 2. Mathematical Results

In this paper we consider a sphere whose mass distribution is axially but not spherically symmetric, spinning on a flat surface subject to a small friction force that is due to sliding. In [13], Ebenfeld and Scheck presented a detailed analysis of the dynamics of the eccentric spinning sphere on a flat surface where friction is assumed to be only due to sliding, see also [14]. Without making any other assumptions, we show that their results imply a full qualitative understanding of the asymptotic long-term dynamics. Whereas the treatment of [13, 14] is mainly analytical, here we adopt a bifurcation theory point of view leading to a global geometric understanding of the phase diagram. The phase diagram in Figure 2 and bifurcation diagrams in Figure 3 illustrate our main results. Recall that an -limit set of a dynamical system is a closed invariant set that is accumulated by a (forward) trajectory [17]. Our main result is summarized in the following theorem.

Theorem 2.1. *A spinning eccentric sphere on a flat surface with small slipping friction admits three types of (asymptotically stable) -limit sets:*(i)*vertically spinning top (), which has its center of mass straight below its geometric center,*(ii)*vertically spinning top (), which has its center of mass straight above its geometric center,*(iii)*intermediate spinning top (), whose center of mass is neither straight below nor straight above its geometric center.**These are solutions of constant energy that are purely rolling due to the assumption on sliding friction (i.e., they display no slipping). The vertical states are periodic, whereas, in general, the intermediate states are quasiperiodic.**At most two of the above types of solutions can be stable at the same time. In case the stable solution is unique its basin of attraction consists of almost the entire phase space (subset of full measure). If the system is bistable, the separatrix between the two different domains of attraction for the asymptotically stable states is expected to be formed by the stable manifold of an unstable intermediate spinning top solution.*

All the analytical results needed to arrive at this conclusion can in principle be found in [13, 14]. However, these papers stop short of drawing the full global conclusions as formulated in the above theorem, and also crucially they did not present the phase diagram and bifurcation diagrams that we present in Figures 2 and 3.

We note that for the eccentric sphere in regime I, the state is always asymptotically stable and thus does not display the tippe top phenomenon. Similarly, no such dynamics arises in regime III since there the inverted state is always unstable. Tippe top dynamics may occur in regime II if the Jellett invariant (which is proportional to the initial spin) is sufficiently large, corresponding to the empirical observation that tippe top dynamics requires a sufficiently large initial spin. In the subregimes IIb and III it is also possible to observe *incomplete* tippe top behaviour, where the top rises but converges to an intermediate state instead of rising all the way to the vertical state . Note that, in regime I, tipping might occur if the top is initially spun sufficiently fast under an angle not close to .

It is important to recognize the existence of symmetries. Recall that symmetries are transformations of the phase space that map solutions of a system to other solutions. In our model of the eccentric sphere on a flat surface, symmetries arise due to the homogeneity of the surface on which the sphere moves and the rotational symmetry of the eccentric sphere. The combined symmetries are thus the Euclidean group (acting as translations and rotations in the plane) and the rotation group acting as rotation of the sphere around its axis of symmetry.

It turns out that the -limit sets mentioned above are all relative equilibria with respect to the symmetry group , see Section 4. Recall that relative equilibria with respect to a group are equilibria for the associated flow on a reduced phase space that is obtained from the original phase space by taking the quotient with respect to the action of . The existence and type of such relative equilibria depend solely on the inertia ratio, the eccentricity of the sphere, and the Jellett integral of motion. We identify a number of regimes characterizing the relative equilibria as a function of the Jellett invariant (which is proportional to the initial angular velocity). The vertical states and are always relative equilibria and their stability depends on the inertia ratio, the eccentricity of the sphere, and the Jellett integral of motion. In addition, intermediate states may exist, which branch off from the vertically spinning solutions. We sketch the phase diagram in Figure 2. For the labeled regions in this phase diagram, the corresponding bifurcation diagrams for the relative equilibria are presented in Figure 3.

The proof of Theorem 2.1, which builds upon the results by [13, 14], can be found in the appendix. In order to present our point of view clearly and in a self-contained way, in Section 3 we also present a derivation of the equations of motion of the eccentric sphere model of the tippe top, including a discussion of the symmetries and their consequences. Here also, one finds a precise description of the assumed nature of the friction and a definition of all the relevant variables that appear as parameters in Figures 2 and 3. In Sections 4 and 5 the relative equilibria of the system and their stability are discussed. The readers who are acquainted with the topic can start reading from Section 6.

We would like to point out that the strategy of proof used here may well be applicable to a large number of similar examples of mechanical systems under the influence of some kind of friction, such as the *Rattleback* [18] or *Hycaro* tippe top of Tokieda [19]. The key observation is that for mechanical systems under the influence of friction, in a natural way the energy becomes a *Lyapunov function* since friction causes energy loss. The next observation is that orbits which do not dissipate energy need to lie entirely on a subvariety of the phase space that is defined by the condition that friction is absent. Equilibria naturally lie on this subset since they have zero velocity and friction is absent at zero velocity. However, one would expect that typically no solution lies entirely on this subvariety, unless the solution lies on the orbit of a symmetry group that leaves the zero friction subvariety invariant. (One can make this precise by constructing a small local perturbation that moves a solution off the zero friction subvariety.) In many cases, the set of such relative equilibria can be accurately analyzed, either analytically or numerically, and the local stability properties can be deduced from a dissipation-induced instability point of view, based on the local stability properties of the relative equilibria in the absence of friction. We note in this respect that the set of -limit sets on the zero friction subvariety is independent of the form and size of the friction. The final step is to draw global conclusions from this local information. The latter is within reach if one has a good understanding of the -limit sets. We would like to stress that Theorem 2.1 concerns the asymptotic dynamics. In an experiment with small friction, the observation may well be dominated by *transient* dynamics which bears strong resemblance (on short time scales) to the dynamics of the spherical top without friction. The dynamics of the latter is rather complicated as it is a nonintegrable Hamiltonian system. In Section 7 we present some results from numerical simulations demonstrating explicitly some examples where the transient dynamics does not appear to prevent fast convergence to the asymptotic states (although of course for sufficiently small friction coefficient the transient dynamics would dominate on finite time intervals).

#### 3. The Equations of Motion

We consider an *eccentric sphere* as in Figure 4, where denotes the center of mass and the center of the sphere. The line joining the center of mass and the geometrical center is an axis of inertial symmetry: in the plane perpendicular to this axis the moment of inertia tensor of the sphere has two equal principal moments of inertia.

We describe the motion of the sphere using three reference frames:(I)an inertial (laboratory) frame , where is some point on the table and the -axis is the vertical.(II)a (noninertial) rotating frame whose origin is in the center of mass and whose 3rd axis is always parallel to the vertical. The and axes are specified below.(III)a principal axis system , whose -axis is the symmetry axis of the sphere (Note that in [5] the origin of the reference frames (II) and (III) is at the center of the sphere, and not at the center of mass.)

The reference frames (II) and (III) are indicated in Figure 4. The eccentricity is the distance between the center of mass and the geometric center of the sphere, with , where denotes the radius of the sphere. We denote the moments of inertia and . The point of contact with the plane of support is denoted by .

Let be the Euler angles of the body with respect to , see Figure 5 for an illustration. The -plane contains the vector which joins the center of mass to the point of contact. The plane is inclined at angle to the vertical -plane and processes with angular velocity around the vertical . We choose the horizontal in so that is perpendicular to . For the rotating frame, is in and perpendicular to the symmetry axis and the axis coincides with . (Note that the axes and are principle axes, but they are not body fixed.) The angle between the vertical and the axis of the top is denoted by . The angular velocity describes the nutation of the body in the vertical plane . The angle describes the orientation of the body with respect to the frame and is the spin of the sphere around its symmetry axis. We denote by the unit vectors along and by those along , and . Note that . Because of the inherent translational symmetry (of the body on the plane), it is convenient to describe the body in terms of the relative (moving) reference frames (II) and (III), rather than in the absolute reference frame (I). By doing so we thus ignore the translational motion on the plane and focus on the relative motion of the body, which captures the tippe top behaviour.

The (relative) position vector of the body is
or . Note that the coordinates of the reference frames (II) and (III) are related by the relations
The reference frames (II) and (III) rotate with respective angular velocities
The angular velocity of the body involves, in addition, the angular velocity :
where denotes the component of about , better known as the *spin*. For later use we introduce the notation . Consequently, with denoting the inertia tensor of the sphere, the angular momentum of the sphere is given by
The point of contact has coordinates
The velocity of the point of contact is
where denotes the vector from the center of mass to the point of contact and is the velocity of the center of mass. We set and use the fact that to obtain
Hence, we obtain
The fact that the sphere remains in contact with the table is expressed by the (holonomic) constraint
where denotes the height of the center of mass above the table, compare Figure 4. From this constraint it follows that the component of is
so that the -component of vanishes, consistent with the constraint.

We note that the physical interpretation of concerns the phenomenon of *slipping*. In case the body slips on the surface. In contrast, a rolling motion of the body is characterized by the fact that .

The equations of motion will be derived, in Newton’s spirit, as a consequence of the action of external forces. We distinguish the following forces acting on the sphere(i)The gravitational force: , where is the total mass of the sphere.(ii)A force acting on the point of contact , where , is the normal reaction force at (due to the stiffness of the surface) and is a friction force. For completeness we mention [20] where a mathematical model for the tippe top is proposed taking elasticity properties of the table and tippe top into account.

Friction is the resistive force acting between bodies that tends to oppose and damp out motion. Friction is usually distinguished as being either static friction (the frictional force opposing placing a body at rest into motion) or kinetic friction (the frictional force tending to slow a body in motion). Importantly, we assume that the friction force is entirely due to the slipping of the sphere on the surface and neglect all other sources of friction. Friction forces can be complicated, and there are various models in circulation. We adopt the assumption of *viscous friction* [9, 15, 21] and assume the friction force to be given by
where is the coefficient of *sliding* friction with the dimension of (velocity)^{−1}. is proportional to the size of the normal reaction force and vanishes smoothly when . (An alternative model for the friction force is the so-called Coulomb friction . This model is not appropriate when due to the singular nature of this force when .) Euler’s equations of motion for the sphere,
govern the evolution of the angular momentum in a noninertial reference frame, rotating with frequency , due to the influence of the external torque . The equation of motion for the center of mass in the rotating frame is
In terms of the coordinates in reference frame (III) the equations of motion (3.13) yield
where and . From the equation for the motion of the center of mass (3.14), in terms of reference frame (II), we obtain
Recalling that , from the last of the latter equations we may derive an expression for :
The equations of motion (3.15) and (3.16) can be written as a system of six coupled first-order nonlinear ordinary differential equations in the variables , where , , , and .

Setting for simplicity, these may be arranged in the standard form (when ) It should be remembered that , , and are still functions of the other variables. For instance, from (3.17) and (3.18) one finds Recall that we require that . If this condition fails, the sphere loses contact with the surface. Expressions for and follow similarly from (3.12).

It is important to recognize that some of the structure of the equations of motion (3.18) is due to symmetry. We recall that the symmetries are the Euclidean group (acting as translations, rotations, and reflections in the plane) and the rotation group acting as rotation of the sphere around its axis of symmetry. The effect of the Euclidean symmetry is that the right-hand side of the equations of motion contain no reference to the position of the sphere on the surface. In a similar way, due to the rotational symmetry the equations of motion do not depend explicitly on . The system can be viewed as three coupled systems, where the coupling is of *skew product* type: the evolution of , and does not depend on , and , and the evolution of , and does not depend on . Moreover, note that the position of the center of mass relative to the surface (in coordinates) could in principle be obtained by integrating the velocities and over time. Because of the fact that we take friction into account, Noether’s theorem does not apply, so the continuous symmetries we observe need not (and do not) give rise to conserved quantities. However, it was discovered by Jellett [8] by an approximate argument, and later proved by Routh [6], that the system (3.18) has the following conserved quantity:
Indeed, it follows from (3.13) that , so that
Note that the Jellett invariant can be written as ([10, 15, 21])

#### 4. -Limit Sets Are Relative Equilibria

Our aim is to describe the asymptotic dynamics of the eccentric sphere. Recall that a subset of the phase space is an -limit set if this set is accumulated by (forward) orbits. While the friction force destroys the Hamiltonian nature of the dynamics, it greatly simplifies the asymptotic dynamics. This follows from the fact that in the presence of friction the energy, which is conserved in the absence of friction, is almost always decreasing along solutions.

The energy is given by , where is the kinetic energy with its rotational and translational part and is the potential energy. With our choice of variables we may write and , where .

Lemma 4.1 (see [13]). *The energy is a Lyapunov function (recall that a Lyapunov function is nonincreasing along orbits) for (3.18). In particular,
*

As is parallel (and opposite) to , vanishes if and only if vanishes. Observe that decreases monotonically and hence is a suitable Lyapunov function. Moreover, is analytic and therefore along orbits it is either strictly decreasing or constant. The energy is constant only if , that is in the absence of friction. Thus, the -limit sets must consist of orbits which do not experience friction. We show that such orbits are necessarily *relative equilibria*.

Proposition 4.2. *Solutions have constant energy only if they are relative equilibria with respect to the action of .*

*Proof. *We already concluded that needs to be equal to 0 along any orbit in an -limit set. A straightforward calculation shows that indeed implies that , so that such a solution must be an relative equilibrium.

This observation is in fact what one would generically expect to find. If is a submanifold of the phase space that corresponds to the absence of friction, in general it would be quite unexpected to find a nonequilibrium solution that lies entirely inside .

#### 5. Stability and Bifurcations of Relative Equilibria

Having determined that the relative equilibria are the only possible asymptotic states in the presence of friction, we derive in this section these solutions of constant energy using the explicit equations of motion (3.15)-(3.16), see also [5, 13, 15]. With and , the equations of motion yield , , , , , and These equations have the following three types of solutions. The linear stability analysis differs from [13, 14] in methodology.

*Vertical States*

(1)Vertical state :
The top is spinning about its axle with center of mass straight below the geometric center.(2)Vertical state :
The top is spinning about its axle with center of mass straight above the geometric center.

*Intermediate States*

For these solutions we have , , and are related by
Elimination of from the above yields
Hence, the condition for the existence of intermediate states is

It is natural to divide the solutions into three groups, according to regimes of the parameters and [21].Group I: . Intermediate states exist with Group II: . Intermediate states exist with any .Group III: . Intermediate states exist with

As in [10] we further refine this classification. Note that the intermediate states discussed here correspond to the tumbling solutions discussed in [13].

The intermediate states are completely determined by (5.4), (5.5), and the Jellett invariant . More precisely, combining the square of (3.22) with (5.4) and (5.6), they are obtained by solving The next theorem summarizes the linear stability and local bifurcation results, as depicted in Figure 3 (cf. [21]), and the proof is sketched in the Appendix. We identify six different groups according to how the value of is related to the eccentricity . In the literature, results have previously been expressed in terms of variables and , referring to the spin of an initial condition in a vertical state. We define , which is the value of the spin at for motion with Jellett invariant . Similarly, denotes the spin of the solution with Jellett invariant at . Note that, for a fixed value of the Jellett invariant , these spins are related by . We further define and . Furthermore we define and .

Theorem 5.1. *The bifurcation diagrams of the eccentric sphere spinning on a flat surface with small friction fall in one of the following six categories (Figure 3).**Group I: .*

(i)*The vertical state is stable for any value of .*(ii)*The vertical state is stable if and unstable otherwise.*(iii)*Intermediate states exist for all values of satisfying .** Group Ia: . The entire branch of intermediate states is unstable.** Group Ib: . The branch of intermediate states has a fold point at . The branch with is stable, while the branch with is unstable.**Group II: .*

(i)*The vertical state is stable if and unstable otherwise.*(ii)*The vertical state is stable if and unstable otherwise.*(iii)*Intermediate states exist for all . We distinguish the following three subgroups.**Group IIa: and . A fold bifurcation of intermediate states occurs.*

*Group IIb: or . The entire branch of intermediate states is stable.*

*Group IIc: and . The entire branch of intermediate states is unstable.*

*Group III: .*

(i)

*The vertical state is stable if .*(ii)

*The vertical state is unstable for all .*(iii)

*Intermediate states exist for and are all stable.*

The proof of Theorem 5.1 can be mainly recovered from [10]; for completeness we provide the calculations based on a direct approach in the appendix.

#### 6. Prototype of a Spherical Tippe Top

Rapid prototyping (RP) technologies enable solid models to be obtained from designs generated with CAD applications. Their increasing popularity in industry is due to the reduction in cost and time associated with the use of these models when verifying product development stages and improvements in end quality. These technologies can also be applied to verify the correctness and/or accuracy of mathematical models and last but not least to enhance students’ active learning in the frame of a learning-by-doing approach. Students can bring their designs to fruition and develop a deeper insight in abstract concepts. We made a prototype of a spherical tippe top for educational use in the Product Development Laboratory of Howest.

As pointed out earlier, to realize a 3-in-1 toy an axially symmetric sphere where one has control over and is needed. We considered three possible designs:(1)a solid sphere with a cylindrical hole through the center where a setscrew can move;(2)a hollow sphere with a cylindrical rod on which a weight can move;(3)a hollow sphere with a toroidal band at the equator fitted with a cylindrical rod on which a weight can be screwed.

From the bifurcation diagram in Figure 2, it is clear that it is possible to hit the three main groups by fixing and changing . Therefore, it is important to understand for the three designs how and vary with respect to each other when the weight is moved. We set up a Maple worksheet based on the given mathematical description and calculated and in function of the position of the midpoint of the moving weight with respect to the center of the sphere; this will be further on denoted by . We took into account the physical parameters: dimensions of the different parts (radii, heights, and thickness) and the density of the materials.

From this we realized that for the solid sphere the goal of the three types is within reach, whereas for the hollow sphere the design has to be modified. Our modifications resulted in the third design as given above. We now discuss our findings for the realized prototypes. Our realizations were all printed with the commercial available Dimension SST1200es with printing technology based on the FDM principle (fused deposition modeling) in ABSplus.

##### 6.1. Sphere with Cylindrical Hole and Setscrew

For the first design we realized three different tops, varying the geometrical dimensions. This was done because the calculation showed that for the given materials some zones are hard to achieve or are very narrow, see Figure 7. The prototype consists of a sphere with a cylindrical hole through the center, together with a piece of adjustable cylindrical iron wire (setscrew), see Figure 6. With a caliper, it can be checked how deep the setscrew is set in the hole. The position of the midpoint of the setscrew with respect to the center of the sphere is denoted by . The hole is suitable for a setscrew M12. The dimensions of the toy were chosen based on the mathematical calculations derived from the model. The diameter of the sphere was chosen so that one can comfortably spin the toy by hand. With a sphere of diameter 50 mm, good values for the chosen design are a hole of radius 5.5 mm, filled with the setscrew of height 15 mm, or a hole of radius 1.5 mm, filled with the setscrew of height 3 mm, see Figure 7. The densities are 1.08 g/cm^{3} for ABSplus and 7.87 g/cm^{3} for the setscrew.

The prototype is axially symmetric; therefore, only the eccentricity and the moment of inertia are functions of , they are easily calculated; remains constant when moving the setscrew up and down. In Figure 7 the quantities , , and are plotted as functions of , (a) for the prototype with a M12 setscrew and (b) for the prototype with a M3 setscrew. The printed prototype in Figure 6, is of the first type and according to the mathematical calculations will exhibit the predicted behaviour as follows: for between mm and mm the toy does not show tippe top dynamic no matter what the initial spin is (type I), for between mm and mm complete tippe top dynamic is observed (type IIc). For above mm the top is of type IIb (incomplete tipping is observed if initial spin is not sufficiently high).

For the different positions of the weight, we launched by hand the toy plenty (≥50) of times and registered each time tipping, nontipping, or hanging. In Table 1 we report the typical results for 5 launches. Note that the tippe top is hand spun, so there will be a deviation from the starting position . Tipping and nontipping were mostly observed in the setup of type IIa and I, respectively. For the setup of type IIb, the expected behaviour (tipping) was not observed; we always observed the hanging behaviour. This is because we were not able to launch the toy fast enough by hand and also because the setscrew sticking out the toy does not allow ideal launching position. Our observations indicate that the prototype behaves as predicted by the model. In details, tipping and hanging at can be explained by the presence of a branch of intermediate states and stable position that one could hit if the toy is not launched exactly from the position, see Figure 3. The hanging at is due to the fact that we did not launch the top fast enough.

For a prototype fit for a M3 screw, the intervals for are as follows: for between mm and mm the toy does not show tippe top dynamic no matter what the initial spin is (type I); for between mm and mm complete tippe top dynamic is observed (type II). For above mm the top is of type III (incomplete tipping is observed). Also, this prototype was spun ≥50 times, and we registered similar observations as for the previous one.

We conclude that this prototype can give a working 3-into-1 toy but has some disadvantages.(i)The setscrew can come loose after intensive use.(ii)When there are three zones present, at least one of the zones is small.(iii)Using a caliper to know if the setscrew is in the center is not practical.

Several attempts were done in the computations to improve the design, for example, by adding holes into the solid sphere. These attempts were not successful, so no other prototype was printed. Instead, we concentrated on the sphere with a cylindrical rod.

##### 6.2. Sphere with Cylindrical Rod

The second design consists of a spherical shell with a cylindrical solid rod through the center along which a symmetric bead is spun; this bead can be put at different heights along the rod. In this design, the user must open the sphere and change the position of the weight by screwing it up or down, after which the sphere can be closed and spun. See Figure 8.

The advantage of this design is that different ABS colors can easily be used for both sides making the tipping more visible and that the rod can be marked at the critical positions. Where in the first design setscrews of different length can be used, in this design weights of different lengths and different widths can be considered within one toy. Also, different materials for the rod can be considered.

However, computations show that a 3-into-1 toy was difficult to obtain with the chosen materials: ABS, iron, and nylon. We briefly summarize our findings that form the basis of further improvements that lead to the third prototype. We tried three different possibilities for the rod: iron, nylon, and ABS. Nylon and ABS seem to work best. With an iron rod and physical parameters that allow easy playing with the toy we did not succeed in catching all the three zones. See below for the specific values of the parameters. As Figure 9(a) shows, only zone three is hit, which means that this tippe top never shows complete inversion but may tip up to a certain angle.

The prototype is axially symmetric; therefore, only the eccentricity and the moment of inertia are functions of , they are easily calculated; remains constant when moving the weight up and down. The physical parameters for the construction are radius of the spherical shell of 25 mm, thickness of the shell of mm, radius of the cylindrical rod of mm, radius of the weight of 10 mm, and height of the weight of mm.

In the case of a nylon or ABS rod, it was possible to obtain a tipping top. As illustrated in Figure 9(b), both zones II and III are hit, which means that complete tipping and incomplete tipping may be observed according to the position of the weight. Note that the section of the blue curve in zone III is very small, which makes it very difficult to observe the hanging phenomenon. A similar observation holds for the section of the blue curve in zone I; this is so small that nontipping behaviour cannot in practice be observed. Many launches of this prototype only confirmed the observations above. This seems to be a good type II tippe top, but we were not able to observe the other behaviours.

##### 6.3. Sphere with a Toroidal Band and a Cylindrical Rod

The third design consists of a spherical shell with a toroidal band around the equator and a cylindrical solid rod through the center along which a symmetric bead is spun. Also in this design, the user must open the sphere and change the position of the weight by screwing it up or down, after which the sphere can be closed and spun. See Figure 10.

Adding a toroidal band was a way to find a compromise between a spherical shell and a solid sphere. The prototype is still axially symmetric, and the band provides a better click system to open and close the toy. The physical parameters for the construction are radius of the spherical shell of mm, thickness of the shell of mm, radius of the cylindrical rod of mm, radius of the weight of mm, and height of the weight of mm. The band has the form of a solid of revolution generated by an ellipse rotating around the rod, the semiaxes of the ellipse measure, respectively, mm and mm. The rod is made of iron. According to the mathematical model, the top should behave as follows: for between mm and mm the toy does not show tippe top dynamic (type I); for between mm and mm complete tippe top dynamic is observed (type II). For above mm the top is of type III (incomplete tipping is observed).

The toy was launched many times, allowing to observe without problems zones I and II predicted by the model. Zone III is, however, difficult to observe at first, the maximum that can be obtained with the constructed prototype was only 15 mm (due to the fastening system of the rod), and moreover which is not so easy to see. Different weights can, however, be used so that zone III becomes visible; for example, this is the case when using an iron weight with radius mm and height mm. We conclude that this last prototype is a good candidate for the 3-in-1 toy, although some further optimizations of the physical parameters can be considered (Figure 11).

#### 7. Numerical Illustrations

##### 7.1. System Trajectories

In this section we present some simulations of (3.15)-(3.16). We focus on the parameter regime of Group II since these are the tops exhibiting “tipping” behaviour. Indeed, if the initial spin , then tipping occurs. The trajectories lie in the (reduced) 6-dimensional phase space (we ignore the equation); here we show their projections in the 3-dimensional subspace of the variables .

Figure 12 shows a number of trajectories for a tippe top of Group IIb starting from initial conditions and . Other input parameters are , the friction coefficient , the eccentricity , and inertia ratio . Points with are stable whereas those for which are unstable. Let be the value of the initial spin calculated for the angle at the Jellett where the change in stability for the inverted position () occurs. Trajectories originating near an unstable noninverted position are attracted either to one of the intermediate states at an angle when or when to a steady state for which ; in this case, the ball rises fully to a stable inverted vertical (tippe top) position with a final spin determined by the Jellett. The blowups in the insets show oscillations in the immediate neighborhood of the fixed points; these depend on the precise initial conditions. We note also that changing does not affect the final destination of the trajectories but it does affect the time needed to follow such trajectories in phase space, this at least within our parameter range of computations.

Figure 14 shows a number of trajectories for a tippe top of Group IIa with friction coefficient , starting now under an angle close to the inverted position . Physical parameters are , , and .

Recall from Figure 3 that trajectories starting near the inverted position will, depending on the initial spin, remain in the neighborhood of go all the way down to the noninverted position , or reach a stable intermediate state. For a clearer overview of the possible behaviours we use the symmetry and sketch in Figure 13 the curve of intermediate states (5.10) in -plane also indicating the essential ’s at which changes in stability type occur. Recall that for a given the relations (5.6) and (5.4) determine and .

The depicted trajectories have been obtained with the initial conditions , is arbitrary, and is one of the following: , , or . These choices were made to reduce to a minimum in the drawings the initial oscillations in the -direction. We start near the equilibrium on the eigenvector of the positive eigenvalue(s) so that the motion quickly evolves in the unstable direction away from the initial position. Trajectories originating near an unstable inverted position will either reach a stable intermediate state at when or fall in the noninverted vertical position when . Here is the value of the initial spin corresponding to the angle at which intermediate states change stability type. Points on with are stable and trajectories starting near the inverted state will be attracted to it.

##### 7.2. 3D Animations

In this section we comment on the 3D animations illustrating the phenomena of “tipping” or “hanging on an intermediate state” for an eccentric sphere (http://cage.ugent.be/~bm/tippetop/tippetop.html). In the films, the eccentric sphere is drawn as a transparent ball with a top in it. We focus on eccentric spheres belonging to Groups IIb and IIa, see Figure 3 for the corresponding bifurcation diagrams. The films have been made using Maple to solve the ODE system and feeding the results to Povray, Imagemagick and ffmpeg. The films are in 5x slow motion for the sake of clarity, so 1 second takes 5 seconds in the animation, with 30 frames for every second. For the sake of clarity the evolution of the nutation angle is shown in each film.

The films for a top of Group IIb show a complete flip (See http://cage.ugent.be/~bm/tippetop/tippetop_IIb_flip.mpg), and the rising to a stable intermediate state (http://cage.ugent.be/~bm/tippetop/tippetop_IIb_IntSt_Comp.mpg) for a top of Group IIa the film shows how the top launched upside-down migrates to a stable intermediate state (http://cage.ugent.be/~bm/tippetop/tippetop_IIaComb.mpg).

In the first film one sees a complete flip (tippe top effect) of the sphere; the physical parameters used are gram, cm, , mg/m^{2}, and mg/m^{2}. The friction coefficient is 0.3. We show 90 seconds, which is presented as 7.5 minutes of the film.

The second film shows the motion towards a stable intermediate state from the unstable noninverted or unstable inverted position. The initial data is chosen so that the Jellett coincides with that of a stable intermediate state. It corresponds to 45 seconds of the tippe top movement. The physical parameters are here as in the first film except for the friction coefficient which is now . The initial conditions around and are chosen to correspond to a Jellett of approximately . The intermediate state is at degrees.

The animations for a top out of Group IIa are meant to illustrate how, depending on the initial conditions, the top started at the inverted position () can fall either to an intermediate state or to the noninverted position (), see Figure 3. The physical parameters are gram, cm, eccentricity , mg/(m^{2}), and mg/(m^{2}). The friction coefficient is 0.08. In the left of the animation we see the behaviour in which the tippe top flips towards the noninverted state. However, oscillation of theta occurs. To the right, we see the movement for a slightly different initial state, with motion towards an intermediate state. Also here oscillation of theta occurs.

#### 8. Further Remarks

*Global Dynamics*

Concerning Theorem 2.1 we wish to stress the importance of the local bifurcation diagrams for the global dynamics. Clearly, if we have a unique asymptotically stable -limit set it is clear that the basin of attraction for this set is equal to nearly the full measure set of the phase space defined by the complement to the stable manifolds of all coexisting (unstable) -limit sets.

It thus remains to analyze the situation when we have two coexisting asymptotically stable -limit sets. From the bifurcation analysis we know that in such case the coexisting stable -limit sets are the vertical states.

Theorem 8.1. *The -limit sets of the eccentric sphere on a flat surface with small friction are asymptotically stable relative equilibria (with respect to ). The -limit set is a unique relative equilibrium, in which case the basin of attraction is the complement of the stable manifolds of the unstable relative equilibria (and hence dense, and of full measure, in the phase space). Otherwise, there are at most two stable relative equilibria (the vertical states) and the system is bistable. In this case, the union of the basins of attraction of the two vertical states is the complement of the stable manifold of the unstable relative equilibria, which is an intermediate state. This union is dense and of full measure in the phase space. The separatrix between the two basins of attraction (inside a level set of the Jellett invariant) consists of the stable manifold of an unstable intermediate state.*

*Proof. *Most of the above statement is a direct consequence of the existence of the energy as a Lyapunov function (through La Salle’s principle). One readily verifies (from the local bifurcation analysis) that the stable manifold of the intermediate state coexisting with two asymptotically stable vertical states has codimension one (inside the level set of the Jellett invariant) and divides the phase space into two parts.

The regions of bistability (as a function of the Jellett invariant ) follow from the local bifurcation diagrams discussed in Theorem 5.1.

*Remark 8.2. *Note that the specific ’s where the changes in stability type of the steady states occur do not depend on . The viscous friction influences the time needed for an orbit to reach such a point. This fact was already clear in [13] and could be proved in advance also in our setup. The result remains true also for more general forms of friction laws proportional to such as those proposed in [4]. This suggests that the study of the asymptotic dynamics of other (mechanical) problems, as for example the rattleback, might notably simplify by the introduction of viscous friction in the model.

*Rolling Model*

A “rolling” eccentric ball does not tip. In this section we give a simple argument showing that, if pure rolling is assumed, then the tippe top phenomenon cannot occur.

Solving the tippe top under the constraint of pure rolling (i.e., when the nonholonomic constraint is satisfied) allows for complete reduction of the equations of motion to a second-order . See [3] for a discussion of this approach. In the pure rolling regime the system is not anymore dissipative and admits three conserved quantities: the energy, , the Jellett as before, and the Routhian, *Routh*, given by [3]

“Tipping” in the rolling model would violate the conservation of . Indeed, from the Jellett invariant we know that the sign of has to change in a complete inversion since and . But this is not allowed if *Routh* = constant has to hold.

The motion in the rolling model is governed by a functional relation of the type . Indeed, the conserved quantities give three relations for the components of the angular velocity . In details, for a given , the Routhian (8.1) fixes , then the Jellett fixes , and finally the energy fixes the tipping rate (cf. (3.4)), yielding a functional relation of the type . Note further that the constraint gives , from (3.9). In this approach, as mentioned in [3] one has to check whether a found solution is physically possible, that is, one has to take into account that rolling cannot be sustained if , where is the coefficient of static friction. In [3] it is remarked that only a few pure rolling precessional solutions satisfy this condition. The analysis, however, leaves open the possibility to have pure rolling periodic motions around the intermediate states as we discuss below.

*Sliding versus Rolling*

A debatable issue is whether transitions between sliding and rolling are possible during the motion of the top. As it was pointed out in [9], such transitions must also be considered when setting up a realistic model to describe the dynamics of the tippe top. To the different regimes there correspond different sets of equations. A switch between sliding and rolling occurs as the absolute contact velocity vanishes. In the Coulomb-friction model, a switch from rolling to sliding occurs when the tangent reaction force required to maintain rolling exceeds . We refer to [4] for considerations and simulations on this topic and to [3] for a detailed analysis of the pure rolling model. We consider two test cases, the pendulum motion and the behaviour of the tippe top around a stable intermediate state.

#### Appendices

#### A. Proof of Theorem 5.1

For the interested reader, the following sections contain the calculations needed for a straightforward linear stability analysis of the steady states. These form the proof of Theorem 5.1.

##### A.1. Stability of the Vertical State

With the Taylor expansions in linearizing in and noting that and , the linearization of the equations of motion (3.16) and (3.15) at yields Introducing the complex coordinates(A.3)–(A.6) can be reduced to two complex equations. The addition (A.3) (A.4) yields whereas (A.5) + (A.6) leads to

These equations admit a solution of the form when satisfies the determinant equation When the roots of (A.10) are where we set In the absence of friction, that is, when , the vertical state is marginally stable as are purely imaginary since .

We now analyze the effect of small friction () by examining how the roots (A.11) are perturbed to first order in : As , the vertical state is stable if . This is the case when It follows that, in Group I (), the vertical state is always stable, while for Group II and Group III () stability requires that .

It remains to be shown how yields the relation (A.14). We focus on the inequality , the arguments are similar for . The inequality yields Using (A.12), this gives Note that if the above condition is satisfied for all . If on the other hand the inequality holds, then ; squaring both sides yields where . Since we are in the case , we can rewrite this last condition as (A.14).

*Remark A.1. *Ignoring translational effects, that is, throwing everything in the variable away (cf. [5]), one is left with
where we set for simplicity , . Equation (A.19) is of Maxwell-Bloch type [5] and allows us to recover the analysis carried out in [5]. An analogous result holds when linearizing around .

##### A.2. Stability of the Vertical State

The stability of the vertical state is studied in a similar way as in Section A.1. From the equation of motion (3.15)-(3.16), introducing complex coordinates we obtain the coupled complex equations The corresponding determinant equation for eigenvalue is given by When , the roots of (A.22) are Thus, at , the vertical state is at most marginally stable, when is purely imaginary, that is, when This is of course also a necessary condition for stability if is small. Provided , corresponding to a resonance at , the roots of (A.22) are perturbed at order to Thus, for stability we have to require , which yields Condition (A.26) is never satisfied for Group III, so is unstable. In the case of Groups I and II, when , the condition for stability is

Note that for tippe tops of Group I and II , with . The equality holds when .

##### A.3. Stability of Intermediate States

In this section we consider intermediate asymptotic states, which exist if the condition (5.7) is satisfied. Such a state, if it exists, is of the form , with constant and related by (5.4). In order to study the stability properties, we study the eigenvalues of the -reduced equations of motion, obtained from (3.18) by omitting the equation.

With denoting the corresponding Jacobian of this reduced equation, eigenvalues satisfy the determinant equation where is the (six-by-six) identity matrix and is a polynomial of degree 5 in . So, 0 is always a solution of (A.28).

*Remark A.2. *It is not possible to reduce the system to a “Maxwell-Bloch” form around an intermediate state as on the contrary was the case around the two vertical spin states.

When the six roots of (A.28) are where and is given by (5.6). Note that from (5.7). All the eigenvalues , are on the imaginary axis, and the intermediate states at are marginally stable.

*Remark A.3. *It is worth noting that for the resonance occurs when . Using the expressions for and given before, one checks that this equality is satisfied when the equation
admits a real root between −1 and 1. The resonance disappears when higher-order terms in are added.

When is small, , we write the first-order perturbation of the roots as From (A.28) it follows that and are as before and The coefficients and are given by with where The coefficients are calculated with the help of Maple For tippe tops of Group II they have a fixed sign for varying in . Note that and . The first inequality is obvious, and the second is less straightforward and is proved below. The friction is stabilizing at when .

The sign of depends on only, since the two terms , are never zero for . The zero of is the bifurcation point where the change in stability type happens. When the intermediate states are not anymore stable. For the sake of brevity we refer to [10] for the details on the behaviour of , where the same crucial function is encountered in the stability analysis via the Routhian reduction. This analysis completes the proof of Theorem 5.1.

*Proof of . *To prove that it is sufficient to prove that , see (A.33), since and . The sign is determined by sign. Considering as a square polynomial in , , we have that the discriminant of is given by
which is negative for all , . Hence, the sign of remains fixed. It is easily verified that for , or , sign is always negative (e.g., , hence for all , and . We conclude that .

#### References

- J. Perry,
*Spinning Tops and Gyroscopic Motions*, Dover, New York, NY, USA, 1957. - C. M. Cohen, “The tippe top revisited,”
*American Journal of Physics*, vol. 45, pp. 12–17, 1977. View at Publisher · View at Google Scholar - C. G. Gray and B. G. Nickel, “Constants of the motion for nonslipping tippe tops and other tops with round pegs,”
*American Journal of Physics*, vol. 68, no. 9, pp. 821–828, 2000. View at Publisher · View at Google Scholar - A. C. Or, “The dynamics of a Tippe top,”
*SIAM Journal on Applied Mathematics*, vol. 54, no. 3, pp. 597–609, 1994. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - N. M. Bou-Rabee, J. E. Marsden, and L. A. Romero, “Tippe top inversion as a dissipation-induced instability,”
*SIAM Journal on Applied Dynamical Systems*, vol. 3, no. 3, pp. 352–377, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - E. J. Routh,
*Dynamics of a System of Rigid Bodies*, MacMillan, New York, NY, USA, 1905. - C. M. Braams, “On the influence of friction on the motion of a top,”
*Physica*, vol. 18, pp. 503–514, 1952. View at Zentralblatt MATH - J. H. Jellett,
*A Treatise on the Theory of Friction*, MacMillan, London, UK, 1872. - T. R. Kane and D. Levinson, “A realistic solution of the symmetric top problem,”
*Journal of Applied Mechanics*, vol. 45, pp. 903–909, 1978. View at Publisher · View at Google Scholar - M. C. Ciocci and B. Langerock, “Dynamics of the tippe top via Routhian reduction,”
*International Journal of Bifurcation and Chaos*, vol. 12, no. 6, pp. 602–614, 2007. View at Publisher · View at Google Scholar - B. Y. M. Branicki, H. K. Moffatt, and Y. Shimomura, “Dynamics of an axisymmetric body spinning on a horizontal surface. III. Geometry of steady state structures for convex bodies,”
*Proceedings of the Royal Society of London. Series A*, vol. 462, no. 2066, pp. 371–390, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - B. Y. M. Branicki and Y. Shimomura, “Dynamics of an axisymmetric body spinning on a horizontal surface. IV. Stability of steady spin states and the `rising egg' phenomenon for convex axisymmetric bodies,”
*Proceedings of the Royal Society of London. Series A*, vol. 462, no. 2075, pp. 3253–3275, 2006. View at Publisher · View at Google Scholar · View at MathSciNet - S. Ebenfeld and F. Scheck, “A new analysis of the tippe top: asymptotic states and Liapunov stability,”
*Annals of Physics*, vol. 243, no. 2, pp. 195–217, 1995. View at Publisher · View at Google Scholar · View at MathSciNet - S. Rauch-Wojciechowski, M. Sköldstam, and T. Glad, “Mathematical analysis of the tippe top,”
*Regular & Chaotic Dynamics*, vol. 10, no. 4, pp. 333–362, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - H. K. Moffatt, Y. Shimomura, and M. Branicki, “Dynamics of an axisymmetric body spinning on a horizontal surface. I. Stability and the gyroscopic approximation,”
*Proceedings of The Royal Society of London. Series A*, vol. 460, no. 2052, pp. 3643–3672, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - R. Bastiaens, J. Detand, O. Rysman, and T. Defloo, “Efficient use of traditional-, rapid- and virtual prototyping in the industrial product development process,” in
*Proceedings of the 3rd International Conference on Advanced Research in Virtual and Rapid Prototyping (VRAP '09)*, Lleira, Portugal, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - V. I. Arnol’d,
*Dynamical Systems. Encyclopedia of Mathematical Sciences*, vol. 3, Springer, New York, NY, USA, 1988. View at Zentralblatt MATH - H. K. Moffatt and T. Tokieda, “Celt reversals: A prototype of chiral dynamics,”
*Proceedings of the Royal Society of Edinburgh Section A*, vol. 138, no. 2, pp. 361–368, 2008. View at Publisher · View at Google Scholar - T. Tokieda, “Private communications,” in
*Proceedings of the Geometric Mechanics and its Applications (MASIE)*, Lausanne, Switzerland, July 2004. - C. Friedl, Der Stehoufkreisel, Zulassungsarbeit zum 1. Staatsexamen, Universität Augsburg, http://www.physik.uniaugsburg.de/%18wobsta/tippetop/index.shtml.en.
- T. Ueda, K. Sasaki, and S. Watanabe, “Motion of the tippe top: gyroscopic balance condition and stability,”
*SIAM Journal on Applied Dynamical Systems*, vol. 4, no. 4, pp. 1159–1194, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet