Phase Portraits of the Autonomous Duffing Single-Degree-of-Freedom Oscillator with Coulomb Dry Friction
The paper presents phase portraits of the autonomous Duffing single-degree-of-freedom system with Coulomb dry friction in its parameter space. The considered nonlinearities of the cubic stiffness and Coulomb dry friction are widely used throughout the literature. It has been shown that there can be more than one sticking region in the phase plane. It has also been shown that an equilibrium point occurs at the critical combinations of values of the parameters and which gives rise to zero eigenvalue of the linearised system. The unstable limit cycle may appear in the case of negative viscous damping ; .
It is beneficial, if not essential, to know the system behaviour prior to using it in any other application or analysis. Phase portraits of the autonomous Duffing single-degree-of-freedom oscillator with Coulomb dry friction are studied here in a complete parameter space.
In order to understand the phase portraits of the nonlinear system in the next section, it is necessary to review the phase portraits of the linear system at different parameter values.
Without loss of generality, the autonomous linear single-degree-of-freedom (SDOF) system’s equation of motion can be written as follows: where represents viscous damping in the system and represents system’s stiffness. By defining the discriminant the system eigenvalues are as follows: and . If then and the response is . If then and the response is . Finally if then the eigenvalues are complex conjugated, , and the response is . and are constants dependent on the initial conditions.
The equation of motion (1) can be represented as system of first-order differential equations: where and and finally is the phase velocity. Hence, There are only two possibilities regarding the number of the equilibrium points . Either there is only one equilibrium point at the origin of the phase plane or there is infinity number of equilibrium points situated at zero velocity abscissa of the phase plane. Table 1 presents all possibilities of the phase portraits of the linear system.
Some aspects of the viscous and frictional damping acting together are presented in  on the linear system with the Coulomb dry friction.
Negative viscous damping is normally associated with self-exciting systems. They are most commonly found in solid-fluid interactions with the fluid induced vibrations [3, 4], but certainly not exclusively [5, 6]. There are two main reasons for using the self-excited model: either the excitation itself is too difficult to model or only the consequences of the energy-input mechanism are important for the problem under consideration.
Both mechanisms of energy dissipation have been used in many studies concerning the forced system, for example, [7–9]. It was shown that separating the viscous damping and frictional damping from the system response is possible, but difficult [10, 11]. Hence, the interactions between both mechanisms of energy dissipation should be well understood.
2. Nonlinear System
Two most commonly used nonlinearities, the cubic stiffness and the Coulomb dry friction, are added to the linear system, thus adding two terms to (1): where stands for the viscous damping ratio ( is equivalent to ), for the undamped eigenfrequency of the equivalent linear system ( is equivalent to ), for the coefficient of nonlinear stiffness, and for the equivalent friction force.
The Duffing oscillator has been widely studied in the literature; a lot of analyses of the Duffing system behaviour are gathered, ordered, and presented in .
Let us suppose that and throughout the text. The use of Coulomb model of dry friction divides the phase plane into two half-planes: the positive half-plane for the positive velocities and the negative half-plane for the negative velocities. The nonsaddle equilibrium points of the SDOF system are shifted to the left from their previous position on the positive half-plane and toward the right on the negative half-plane and vice versa for the saddle points. Between these pairs of equilibrium points the sticking region is formed. Any trajectory in the phase plane that collides with the sticking region stays there forever. Coulomb dry friction model does not directly change the nature of the equilibrium points, as viscous damping does. It only changes their position . The nature of some types of the equilibrium points is directly changed by the viscous damping parameter (or for the equivalent linear system); see Table 1, column . Hence, the equations of motion are
3. Equilibrium Points
The equilibrium points are obtained at the zero phase velocity: where superscript + stands for the positive phase half-plane and superscript − for the negative one. Index stands for the coordinates of the equilibrium points. Equation (6) gives rise to at most three equilibrium points at each half-plane. The nature of the equilibrium points can be deduced from linearized system, (7), linearized in the proximity of each of the equilibrium points, in the case of its nonzero eigenvalues. Consider the following: It is obvious that coefficients and from (7) directly correspond to the linear system coefficients and , respectively; see (3) and Table 1.
The discriminant , (8), is defined first, when solving the cubic equation (6). The value of the discriminant , (8), determines the number of equilibrium points, for example, number of real solutions of cubic equation (6). Consider the following: where If there is only one equilibrium point on each half-plane; on the other hand if there are three equilibrium points on each half-plane. In the case of there are two equilibrium points on each half-plane. Two of the equilibrium points, which exist at , collide and form an equilibrium point which cannot be analysed by the linearised system. The parameter plane is divided into subregions according to the discriminant’s sign; see Figure 1. The zero discriminant surface in the parameter space is presented in Figure 2. The value of as function of at the given value of is presented in Figure 3. The stability of the equilibrium points and their location as a function of the parameter is presented in Figure 4. If then there are no stable equilibrium points.
3.1. Case A
One has . In this case holds and (6) has only one real solution for each half-plane. Consider the following: The constant can be evaluated from (7) and that gives , which defines both equilibrium points as centres if there is no viscous damping. In the presence of the viscous damping the nature of the equilibrium points changes as described by the column in Table 1. Hence, the stability of the equilibrium points is governed by the viscous damping ratio .
3.2. Case B
One has . There are three subcases depending on the discriminant’s sign. As fast as parameter becomes negative, two unstable equilibrium points appear from , one for each half-plane, and other two from , again one for each half-plane, subcase B1. The absolute value of the for the four equilibrium points grows smaller with every smaller value of the parameter as presented in Figure 4. At certain value of the parameter two of the saddle equilibrium points collide with the stable equilibrium points (at ), which exist also when , subcase B2. At this point only two unstable equilibrium points remain, subcase B3.
3.2.1. Subcase B1
One has . Equation (6) has three real solutions for each half-plane; see Figure 4. By defining the abscises of the equilibrium points can be computed as follows: Constants’ ; signs are estimated with the help of (6) and (7) as follows: The first two equilibrium points on both phase half-planes are saddles. The last equilibrium point on both phase half-planes is centre if there is no viscous damping. In the presence of the viscous damping the nature of these equilibrium points changes as described by the column in Table 1.
3.2.2. Subcase B2
One has . Equation (6) has three real solutions for each half-plane, but two of them are the same; see Figure 4. Consider the following: Constants’ ; signs are estimated with the help of (7) as follows: It is clear from (15) that the first equilibrium point, , is saddle. Other two solutions give a single equilibrium point of higher order, which cannot be properly described by linearised system due to zero eigenvalue . However, the centre manifold theorem [17, 18] allows us to assess the stability of such points. According to the centre manifold theorem the system has one centre manifold and one stable manifold at this point. Hence, the equilibrium points are not unstable.
3.2.3. Subcase B3
One has . There is only one real solution of (6) for each half-plane in this case; see Figure 4. Consider the following: The constant can be evaluated from (7) and that gives , which defines both equilibrium points as saddles.
4. Sticking Regions
The smooth transition of the trajectory from one half-plane to the other one is achieved when and sticking appears when . Hence the sticking regions appear in between two equivalent equilibrium points each from different phase half-plane. The sticking region as function of parameter is presented in Figure 4 as the shaded region.
5. Phase Portraits and Basins of Attraction
The phase portraits are important in understanding the global dynamics of the system, especially in the vicinity of the equilibrium points. Equilibrium points marked with the red circles belong to the positive half-plane and the ones marked with the blue squares to the negative half-plane of the phase plane. The basin of attraction is traced by numerical integration of the system response from the initial conditions. If the trajectory ends up in the attracting set the initial conditions are marked as a point in the phase plane. If the trajectory exhibits unstable motion the point is left unmarked. In the case where the whole phase plane is the basin of attraction no initial point is marked.
5.1. Case A
One has . There are two equilibrium points, see Figure 4, which type is governed by the value of parameter . Hence, four main subcases can be recognized here, defined by the value of the parameter . When parameter reaches value of the nature of the equilibrium points is changed from stable node to stable focus or vice versa. On the contrary, when parameter reaches value of the limit cycle shrinks to the equilibrium points and ceased to exist. It is replaced by two heteroclinic orbits.
5.1.1. Subcase A1
One has . The phase portrait is presented in Figure 5. The equilibrium points are attractors-stable nodes and the sticking region is the attracting set. The basin of attraction is whole phase plane.
5.1.2. Subcase A2
One has . The phase portrait is presented in Figure 6. The equilibrium points are attractors-stable focuses (centres if ) and the sticking region is the attracting set. The basin of attraction is whole phase plane again.
5.1.3. Subcase A3
One has . The phase portrait is presented in Figure 7. The equilibrium points are attractors-unstable focuses and the sticking region is the attracting set. The trajectory that repeats itself where repellency of the equilibrium points is balanced by the attraction of the sticking region can be found. The dry friction energy dissipation mechanism attracts the trajectories closer to the sticking region whereas the negative viscous damping mechanism repels them at the larger distances. Thus the unstable limit cycle forms where these two influences are in balance. The phenomenology is mathematically described in  for a simpler system. The amplitude of the limit cycle as a function of is presented in Figure 8. The basin of attraction lies within the limit cycle. The trajectories that start inside the limit cycle end up in the sticking region. The trajectories that start outside the limit cycle end up in . Hence, the limit cycle divides the phase plane into region of stable and unstable motion; see Figure 7.
5.1.4. Subcase A4
One has . The phase portrait is presented in Figure 9. The equilibrium points are attractors-unstable nodes and the sticking region is the attracting set. Limit cycle from case A3 is shirked towards equilibrium points when approaches and vanishes when . The heteroclinic orbits divide the phase plane into regions of the stable motion and one of the unstable motion in this case, Figure 9. The phenomenology is mathematically described in  for a simpler system.
5.2. Case B
One has . There are three subcases, which depend upon the value of the discriminant . The value of the viscous damping has the largest impact in subcase B1 where and no impact at all in subcase B3 where and all equilibrium points are saddles.
5.2.1. Subcase B1
One has . The phase portraits for values of are qualitatively similar. The phase portrait for the value of is presented in Figure 10 and of in Figure 11. The trajectory coming out from saddle points and ends up in the sticking region between equilibrium points and (blue region).
The first qualitative change happens at , Figure 12, where trajectory coming out from the saddle point ends up in the saddle point and thus becomes the heteroclinic orbit. The same goes for the trajectory from to .
When the value of is lessened further on, Figure 13, the trajectory coming out from saddle point ends up in the sticking region between equilibrium points and (green region). And the trajectory coming out from saddle points ends up in the sticking region between equilibrium points and (red region). The unstable limit cycle is formed between equilibrium points and .
Next qualitative change of the phase portrait happens at the value of , Figure 14. Here the trajectory leaving the saddle point ends up in the saddle point and thus becomes heteroclinic orbit. The same is true for the trajectory between and . The unstable limit cycle is still present. It shrinks in comparison with the limit cycle at larger .
Beyond the basins of attraction split among themselves, Figure 15. The unstable limit cycle shrinks further on.
The limit cycle shrinks toward equilibrium points and ceased to exist at , Figure 16. There are two heteroclinic orbits linking and instead of the limit cycle. The phase portraits with even smaller do not change qualitatively.
5.2.2. Subcase B2
One has . The phase portraits, Figure 17, are qualitatively similar till where the trajectory which ends up in equivalent point starts in and thus becomes a heteroclinic orbit, Figure 18, instead in the second quadrant of the phase space. The phase space does not change quantitatively with decreasing the value of further on, Figure 19.
5.2.3. Subcase B3
One has . The nature of two equilibrium points is saddle. Hence, the variation of parameter does not change the nature of the points and thus the nature of the phase portrait, Figure 20.
The paper presents the phase portraits and the sticking regions of the autonomous nonlinear system in its parameter space. The linear system with single-degree-of-freedom (SDOF) is upgraded with two widely used nonlinearities: the cubic stiffness () and the Coulomb dry friction ().
The phase portrait of the nonlinear system in a way resembles the phase portrait of the Duffing system when . In this case the difference is caused by the Coulomb dry friction model that shifts both phase half-planes in opposite directions.
True difference between Duffing and the system under consideration is shown when . It was shown that there are three distinctive sticking regions for certain set of parameters and no close basin of attraction for . When the sticking regions merge together, at the critical combinations of values of the parameters and , an equilibrium point occurs that gives rise to zero eigenvalue of the linearised system.
Even when it is evident that, although the equilibrium points, which are at the same time attractors of the system, are unstable, the attracting set(s), which is (are) sticking region(s), still attracts considerable amount of trajectories.
The unstable limit cycle appears in the case of , in the case of , and in the case of and .
Conflict of Interests
The author declares that there is no conflict of interests regarding the publication of this paper.
J. M. T. Thompson and H. B. Stewart, Nonlinear Dynamics and Chaos, John Wiley & Sons, 1991.
J. P. den Hartog, “Forced vibrations with combined coulomb and viscous friction,” Transactions of the ASME, vol. 53, pp. 107–115, 1931.View at: Google Scholar
N. Jakić, Parameter identification of the second-order autonomous dynamical system [Ph.D. thesis], University of Ljubljana, Ljubljana, Slovenia, 2002.
J. J. Thomsen, Vibrations and Stability: Order and Chaos, McGraw-Hill, 1997.
A. Panfilov, Non-Linear Dynamical Systems, Utrecht University, Utrecht, The Netherlands, 2001.