Nonlinear dynamical systems exhibiting oscillating limit cycles are found in a large variety of fields including biology, chemistry, mechanics and electronics,   . The historical development of the study of periodic solutions of ordinary differential equations began in a series of articles published in the 19th century. It was Henri Poincare that initiated the qualitative theory of differential equations. This qualitative approach was intended both to guide and bring to light the limits of numerical methods, such as the series solutions, which was quite popular at the time. Most of the qualitative investigations were generally of a local character. Thus, the behaviour of solution is studied in a sufficiently small neighborhood of a given solution, for example in a neighborhood of stationary point or of a periodic solution.
A bead sliding on a wire is a classic problem in theoretical dynamics  . Our interest in this problem is from the point of view of studying the resulting nonlinear differential equations. We seek in particular periodic solutions as well as stability properties, using phase portrait analysis. Furthermore, the shape of our wire is specified by a general curve , and hence we can study different configurations of the given problem.
The temporal evolution of a dynamical system is usually described by points in space and time. The line joining all these points is called the trajectory. A trajectory that comes back upon itself to form a closed loop in phase space is called an orbit  . An orbit for a system usually indicates that the dynamical system under consideration is conservative. We also note that each plotted point along any trajectory has evolved directly from the preceding point. As we plot each successive point in phase space, the plotted points migrate around. Orbits and trajectories therefore reflect the movement or evolution of the dynamical system.
The phase space plot is a world that shows the trajectory and its development. Depending on various factors, different trajectories can evolve for the same system. The phase space plot and such a family of trajectories together are a phase space portrait, phase portrait, or phase diagram.
2. Limit Cycles and Other Closed Paths
The world we live in can be described more accurately by systems of nonlinear differential equations. Furthermore, in nonlinear systems there is particular interest in the existence of periodic solutions, as well as their amplitudes, periods, and phases. If the system is autonomous, and in any solution, so is for any value of , which means that phase is not significant since the solutions map on to the same phase paths. In the phase plane, periodic solutions appear as closed paths  . Conservative systems and Hamiltonian systems are essentially energy systems, and often contain nests of closed paths forming centres, which we might expect since they are generally non-dissipative; meaning that friction is absent.
A limit cycle is an isolated periodic solution of an autonomous system  , represented in the phase plane by an isolated closed path. The neighbouring paths are, by definition, not closed, but spiral into or away from the limit cycle. Furthermore, for a stable limit cycle, the device represented by the system, which might be, for example, an electronic circuit, will spontaneously drift into the corresponding periodic oscillation from a wide range of initial states. The existence of limit cycles is therefore a feature of great practical importance.
3. Statement of the Problem
We derive the differential equation of motion of a bead of mass m, sliding on a wire of given shape in the vertical plane. We consider two cases, namely without friction and with friction. For both cases we perform simulation of the derived differential equations. Our interest is to study the resulting nonlinear ordinary differential equations, and investigate the existence of periodic solutions for this dynamics in the neighbourhood of the critical points. This will be depicted by their corresponding phase portraits.
4. Gradient and Hamiltonian Systems
Let E be an open subset of and let where with . A system of the form
is called a Hamiltonian system with n degrees of freedom on E. For example, the Hamiltonian function
is the energy function for the spherical pendulum
This system is equivalent to the pair of uncoupled harmonic oscillators;
All Hamiltonian systems are conservative in the sense that the Hamiltonian function or the total energy remains constant along trajectories of the system.
5. Theorem (Conservation of Energy)
The total energy of the Hamiltonian system (1) remains constant along its trajectories.
The total derivative of the Hamiltonian function along a trajectory of (1)
Thus, is constant along any solution curve of (1) and the trajectories of (1) lie on the surfaces = constant.
We now establish some very specific results about the nature of the critical points of Hamiltonian systems with one degree of freedom. Note that the equilibrium points or critical points of the system (1) correspond to the critical points of the Hamiltonian function where . We may, without loss of generality, assume that the critical point in question has been translated to the origin.
If the origin is a focus of the Hamiltonian system
then, the origin is not a strict local maximum or minimum of the Hamiltonian function .
A critical point of the system at which (Jacobian) has no zero eigenvalues is called a non-degenerate critical point of the system, otherwise, it is called a degenerate critical point of the system.
We note that any non-degenerate critical point of a planar system is either a hyperbolic critical point of the system or a center of the linearized system.
10. Hamiltonian with One Degree of Freedom
One particular type of Hamiltonian system with one degree of freedom is the Newtonian system with one degree of freedom,
where . This differential equation can be written as a system in :
The total energy for this system where is the kinetic energy and
is the potential energy. With this definition of we see that the Newtonian system (4) can be written as a Hamiltonian system. It is not difficult to establish the following facts for the Newtonian system (4).
The critical points of the Newtonian system (4) all lie on the x-axis. The point is a critical point of the Newtonian system (4) if and only if it is a critical point of the function , i.e., a zero of the function . If is a strict local maximum of the analytic function , it is a saddle for (4). If is a strict local minimum of the analytic function , it is a center for (4). If is a horizontal inflection point of the function , it is a cusp for the system (4). Finally, the phase portrait of (4) is symmetric with respect to the x-axis.
12. Existence of Periodic Solutions
The existence of periodic solutions or otherwise of a system of ODE is settled in the following theorem.
Consider the 2-dimensional autonomous system where . Let be simply connected, such that , we have . We now show that the ODE system has no periodic solutions in Ω.
Towards a contradiction, assume ODE system has a periodic solution in Ω. Let ∂Ω be a boundary on Ω.
is the normal to ∂Ω. Recall then Divergence Theorem:
Let u be a periodic solution with period T, i.e. . Then, a path traversed by a solution starting from to is ∂Ω. Then, ∂Ω is a closed curve.
However, by hypothesis, and . Therefore, , and either or on Ω. Thus, or , a contradiction.
We can now state formally a theorem on the Bendixson criterion for the existence or otherwise of periodic solutions  .
15. Theorem (Bendixson Criterion)
Suppose that the domain is simply connected; is continuously differentiable in D.
Then the equations
can only have periodic solutions if changes sign in D or if in D.
We know that the damped linear oscillator contains no periodic solutions.
Consider now a nonlinear oscillator with nonlinear damping represented by the equation
We assume that and are smooth and that , (damping). The vectorized system is
Thus the divergence of the vector function is which is negative definite. It follows from Bendixson’s criterion that the equation has no periodic solutions.
17. Derivation of the Main Models
Case 1: Smooth wire (no friction)
The bead and wire are shown in Figure 1. The components of the velocity of the bead are given by . The kinetic and potential energies are
Figure 1. Bead sliding on frictionless wire.
Hence the Hamiltonian (total energy) H of the bead is
The equation of motion of the bead can now be derived from the conservation principle of the total energy, i.e. . Hence assuming and are continuous,
Hence the differential equation of motion of the bead is;
Case 2: Non-smooth wire (friction)
Let R be the normal reaction of the bead on the wire and let F be the frictional force opposing the motion as shown in Figure 2. The horizontal and vertical equations of motion are
where θ is the inclination of the tangent of the curve at the bead. Since , then and . Eliminating R between (12) and (13) gives;
where and .
The frictional force F is proportional to the velocity and opposes the motion, hence
Figure 2. Bead sliding on rough wire.
Therefore, the equation of motion is
We observe that we obtain Equation (7) when we set for the frictionless wire.
18. Results: Bead Sliding on a Cubic Shaped Wire
Consider the case without friction, and the initial shape of the wire is a simple cubic given by . The equation of motion of the bead is written;
with , and . We examine two cases, viz;
1) For the frictionless case , and setting , we get
2) For the case with friction , and the equation of motion reads;
Vectorizing Equations (12) and (13), we obtain respectively;
19. Simulation of the Models Using MathCAD 
Case 1 (without friction)
We consider the vectorized systems of differential Equations (14)
Define a function that determines a vector of derivative values at any solution point (t, Y):
Additional arguments for the ODE solver are:
Initial value of independent variable;
End value of independent variable;
Vector of initial function values;
Number of solution values on [t0, t1].
Independent variable values;
First solution function values;
Second solution function values.
We have the following graphical profiles (Table 1).
Case 2 (with friction)
We now consider the vectorized nonlinear systems of differential Equations (15)
The vector of derivatives for this case is:
Employing the same additional arguments for the MathCAD ODE solver as previously, we obtain the solution matrix in Table 2, and the following graphical profiles.
Table 1. Solution matrix for the ODE system (14).
Table 2. Solution matrix for the ODE system (15).
We considered the case of a bead sliding on a smooth wire with shape (Figure 3).
The resulting nonlinear differential equation was computed to give
For the case with friction , and the dynamics is given by the nonlinear ODE;
For the former case, Figure 4 depicts the trajectory of the bead against time, on smooth wire with shape . Figure 5 provides the velocity profile. Both graphs are periodic in nature. Correspondingly, the phase portrait shown in Figure 6 depicts an irregular shaped closed curve, in the phase plane.
For the frictional case analysis of the nonlinear ODE (17) provided us with the trajectory of the bead against time, Figure 7, moving on a rough wire with shape . The corresponding velocity profile is depicted in Figure 8, and for both cases we observe constrained oscillations which die out very fast as a result of friction. The dotted lines placed on the graph shown in Figure 7 is meant to indicate that the bead performs oscillations with decreasing amplitude about the minimum point of the function , with x1-coordinate . The phase portrait for this situation, Figure 9, depicts as expected, a spiral sink which actually spirals to the minimum point with x1-coordinate = 0.577.
Figure 3. Bead sliding on a wire with shape .
Figure 4. Trajectory of bead against time, on smooth wire.
Figure 5. Velocity of bead against time, on smooth wire.
Figure 6. Phase portrait of the dynamics of bead on smooth wire.
Figure 7. Trajectory of bead against time, on rough wire.
Figure 8. Velocity of bead against time, on rough wire.
Figure 9. Phase portrait of the dynamics of bead on rough wire.
The present research study is of practical importance in the field of stability analysis as well as existence of periodic solutions of nonlinear ordinary differential equations. In general, nonlinear dynamical systems exhibiting oscillating limit cycles can be found in a large variety of fields including biology, chemistry, mechanics and electronics. Our contribution to the existing body of knowledge in this field is analyzing highly nonlinear ordinary differential equations, without the possibility of solving them analytically, and obtaining important qualitative properties. We employ instead the phase portrait method which is a numerical method, but with the added advantage of providing a pictorial (graphical) view of the inherent dynamics. Another important aspect of the problem considered, is that the geometry of the curve can be prescribed arbitrarily in the vertical plane by the equation . This then makes it possible to study the dynamics of the bead for any given configuration.