# Attractor and Vector Structure Analyses of Bursting Oscillation with Sliding Bifurcation in Filippov Systems

Authors: Rui Qu; Shaolong Li
Journal: Shock and Vibration
 (2019)

---

## Abstract

The main purpose of this paper is to investigate the mechanism of sliding phenomenon in Filippov (nonsmooth) dynamical systems by attractor analysis and vector analysis. A corresponding simple model based on Chua’s circuit with periodic excitation was introduced as an example. The attractor analysis proposed in our previous work is used to discuss the complicated oscillations of the Filippov system. However, it failed to perfectly explain the sliding phenomena and establish an analytical method of constant voltage control. Therefore, the geometric structure and analytic conditions of sliding bifurcations in the generaln-dimensional piecewise smooth system are discussed in detail by vector structure analysis. The prospects of practical application of this method are also discussed in the end.


---

## Body

## 1. Introduction

In practical engineering applications, there are many nonsmooth factors, such as switching [1] and impulse control [2]. The nonsmooth models describing the special dynamic behavior are generally divided into three categories: continuous nonsmooth systems [3], Filippov systems [4], and nonsmooth impulsive systems [5]. The Filippov-type (piecewise smooth) systems are characterized by discontinuous vector fields and continuous Jacobian matrices. The introduction of nonsmooth terms will result in a series of special dynamic behaviors near the interface, such as sliding, traversing, grazing, and jumping [6–9]. Due to the discontinuity of the vector field, the traditional attractor analysis method [10] employed to deal with smooth systems cannot be directly used to study the special dynamic behaviors occurring in nonsmooth region. Therefore, new theoretical methods are needed to analyze such problems.The sliding motions [11, 12] with profound physical backgrounds exist widely in engineering practice. For example, in mechanical systems, the introduction of dry friction factor will result in the sliding phenomenon of mechanical oscillator [13]. Similarly, in the circuit system, sliding bifurcation indicates the conservation of voltage, which has certain guiding significance for the design of a constant voltage transformer. Therefore, the analytical conditions of sliding bifurcation [14] are of vital importance.Taking Chua’s circuit system with bilateral diodes as an example, the bursting oscillation and sliding bifurcation are studied by equilibrium curve, bifurcation diagram, and phase portrait. The validity of the theoretical results is verified by calculation of the sliding time and analysis of simulation data near special points. The application prospects of methods proposed are discussed in the end.


## 2. Mathematical Model

The current-driven Chua’s circuit is shown in the Figure1. The circuit consists of capacitors C1andC2, inductor L, resistors RandRN, and current source is, in which RN is a nonlinear resistor. According to Kirchhoff’s current law and voltage law, the state equations of the circuit can be formed as(1)dvc1dt=1RC1vc1−vc2−1C1fvc1+1C1is,dvc2dt=1RC2vc1−vc2+1C1il,dildt=−1Lvc2,where vc1 and vc2 are the voltages of C1 and C2, respectively; il is the current of L; and fvc1 is the volt-ampere characteristic function of a nonlinear resistor. The excitation current [15, 16] is denoted by is=FsinΩt+Φ in which F,Ω,andΦ represent the amplitude, frequency, and initial phase of the excitation current, respectively. By the coordinate and dimensionless transformation,(2)x=vc1,y=vc2,z=Ril,t=RC2τ,α=C2C1,β=R2C2L,A=RC2FC1,ω=RC2Ω.Figure 1
The current-driven Chua’s circuit.The state equations of Chua’s circuit excited by current can be expressed as(3)x˙=αy−x−fx+Asinωt,y˙=x−y+z,z˙=−βy.where the volt-ampere characteristic function fx of the nonlinear resistance can be expressed in the form of cubic polynomial proposed by Hartley and Mossayebi [17] in the study of deformed Chua’s circuit as fx=mx+nx3, where m=−8/7andn=2/7. It is noteworthy that the nonlinearity is weak due to the relatively small nonlinear coefficient.Taking the simple Chua’s circuit system as an example, dynamic behaviors of Filippov systems with multiscale effect can be studied. By adding a bilateral diode as a nonlinear element, the dimensionless mathematical model of the modified circuit can be expressed as follows:(4)x˙=αy+ax3+bx+w+csgnx,y˙=x−y+z,z˙=−βy,where w=AsinΩt is the external excitation, and when the frequency of external excitation is far less than the natural frequency ω, there will be two time scale effects. Nonsmooth interface x=0 divides the original system into two subsystems located in different phase spaces.


## 3. Bifurcation Analysis of Two Subsystems

There are magnitude differences between external excitation frequency and natural frequency when0<Ω≪1. In any natural frequency period Tω=2π/ω,t∈t0,t0+Tω, the external excitation will change between WA=AsinΩt0 and WB=AsinΩt0+2πΩ/ω. Because Ω/ω≪1, WA≈WB which means that in any natural frequency period, the external excitation is almost a constant. Therefore, the whole external excitation term can be regarded as a slowly varying parameter and two generalized autonomous subsystems can be derived as follows.Whenx<0,(5)S1:x˙=αy+ax3+bx+w−c,y˙=x−y+z,z˙=−βy.Whenx>0,(6)S2:x˙=αy+ax3+bx+w+c,y˙=x−y+z,z˙=−βy.The nonsmooth interfaceΣΣ:x,y,zx=0 divides the phase space into two regions. The corresponding regions of x>0 and x<0 are represented as D+ and D−, respectively. When the trajectory is located in D−, it is restricted by subsystem S1. Once the trajectory reaches D+ through the interface, the governing equation will be converted to S2 immediately. The equilibrium point of two subsystems can be expressed as E0x,y,z=X0,0,−X0, where X0 satisfying(7)αaX03+bX0+w+csgnX0=0.The properties of the equilibrium point depend on the following characteristic equation:(8)λ3+1−αb−3αaX02λ2+β−α−αb−3αaX02−αβb+3aX02=0.Therefore, the conditions for the stability of equilibrium pointE0 are as follows:(9)1−αb−3αaX02>0,αβb+3aX02<0,1−αb−3αaX02β−α−αb−3αaX02+αβb+3aX02>0.Two types of codimension-one bifurcations can be observed, in which fold bifurcations occur under the following conditions:(10)FB:b+3aX02=0,αaX03+bX0+w+csgnX0=0.With1−αb−3αaX02>0, 1−αb−3αaX02β−α−αb−3αaX03+αβb+3aX02>0.Based on the Routh–Hurwitz stability criterion [2, 3], for a cubic polynomial,(11)P3z=z3+a1z2+a2z3+a3.Its Routh–Hurwitz conditions can be formed as(12)a1>0,a11a3a2>0,a3a1a2−a3>0,⇕a1>0,a3>0,a1a2>a3.When conditiona1a2>a3 degenerates to a1a2=a3, a pair of pure imaginary roots appear in the solution of characteristic polynomial (8) indicating that Hopf bifurcation occurs. Therefore the conditions of Hopf bifurcation are as follows:(13)HB:1−αb−3αaX02β−α−αb−3αaX02+αβb+3aX02=0,αaX03+bX0+w+csgnX0=0.With1−αb−3αaX02>0,αβb+3aX02<0.The bifurcation and bursting oscillation behaviors of this system have been studied in our previous work [14]. Therefore, No further elaboration will be made here. The two-parameter bifurcation diagram and the equilibrium branches with bifurcations are shown in Figure 2 as a supplement to our previous work.Figure 2
(a) Two-parameter bifurcation diagram and (b) the equilibrium branches with bifurcations in different regions.
(a)
(b)In Figure2, FBi±and  HBi±i=1,2 corresponds to the different bifurcation points in different segments of equilibrium branches, respectively. The bifurcation lines divided the parametric plane into different regions marked by number ii=1,2,…,7, and the stability of the equilibrium points corresponding to the adjacent regions is opposite. The results of Figure 2 are consistent with those of our previous study.


## 4. Attractor Structure Analysis of Sliding Motion

It is noteworthy that all the bifurcation results obtained are related to slow-varying parameters [18–22]. Therefore, in order to reveal the complete bursting oscillation process and special dynamic behaviors of the system (4), the concept of transformed phase portrait (TPP) [23, 24] should be introduced.


### 4.1. TPP for the TrajectoryΠT:xt,yt,zt,t∈R

The traditional phase portrait reflects the relationship among state variables in the phase spacex,y,z or the projection planes x,y and y,z. If the slowly varying parameter is regarded as a generalized state variable, the generalized phase trajectory can be obtained ΠGT:xt,yt,zt,wt,t∈R, where wt=AsinΩt. The generalized phase trajectories in the generalized phase space x,y,z,w and its projected phase plane are called transformed phase portraits, which are used to characterize the relationship between different state variables and slowly varying parameters.Variable parametersb were taken b=0.14 and b=0.70, respectively, while other parameters were fixed as(14)α=6.9,β=10.0,a=−27,c=−2,A=6.0.Therefore,w varies periodically with frequency 0.01 in the interval −6.0,6.0. Taking A = 0, the natural frequency of the system can be calculated. For b=0.14, the equilibrium point of the two subsystems is approximately E±=±1.1662,0.0,∓1.1662, with its corresponding eigenvalue as follows:(15)λ1=−7.920,λ±=−0.0784±2.988I,indicating that the natural frequency is about 2.988. Similarly, when b=0.70, the coordinates of the three equilibrium points and their corresponding eigenvalues are(16)E1±=±0.4517,0.0,∓0.4517,λ1=4.5159,λ±=−0.9463±2.6698I,E2±=±1.2897,0.0,∓1.2897,λ1=−6.0386,λ±=0.01534±2.8797I,E3±=±1.7414,0.0,∓1.7417,λ1=−13.622,λ±=−0.2417±3.0923I.Therefore, when the trajectory oscillates near different equilibrium pointsEj±j=1,2,3, its natural frequencies are about 2.6698, 2.8797, and 3.0923.


### 4.2. Scenario 1:b=0.14

The 3D phase portrait and bifurcation diagram with parameterb are shown in Figure 3. It can be seen that the attractor structure keeps the symmetry of the original system. The alternation of large-amplitude spiking states and small-amplitude quiescent states can be observed. Specially, the sliding phenomenon appears on the nonsmooth boundary indicating that the voltage of the capacitor remained a constant during this period of time. It should be pointed out that the phenomenon in Figure 3(b) does not mean the changes of periodicity (similar to period doubling bifurcation) but that the position of equilibrium points migrates due to the occurrence of sliding bifurcation which will be discussed in detail.Figure 3
(a) 3D phase portrait withb = 0.14 and (b) bifurcation diagram with parameter b.
(a)
(b)In order to reveal the mechanism of bursting oscillation, the overlap of the transformed phase portrait [14] with the equilibrium branches can be observed in Figure 4. EBjii=1,2;j=1,2,…,6 corresponds to the equilibrium curve segments of subsystem Sii=1,2, respectively. The solid part of the equilibrium line represents the stable equilibrium state, and the dashed line represents the unstable state. The black lines represent the equilibrium points that can be found in the corresponding region, while the red parts denote the equilibrium state that exists but is not located in the phase space of the corresponding system in control.Figure 4
Overlap of the transformed phase portrait with the equilibrium branches.


#### 4.2.1. Bursting Oscillation Mechanism

As shown in Figure4, the trajectory dominated by subsystem S1 starts from point A1 w=−6.0 and runs almost strictly along its stable equilibrium state EB11. Despite the existence of a stable limit cycle LC11 between the two supercritical Hopf bifurcation points HB11 and HB21, the trajectory still runs along the stable equilibrium, showing the quiescent state. The stability of limit cycles produced by these two Hopf bifurcations is opposite. They disappear after contact suggesting the occurrence of LPC bifurcation. When trajectory reaches the fold bifurcation point FB11, it jumps directly to the stable equilibrium. As time t continues to increase, the trajectory is attracted by the only stable equilibrium segment EB61 of subsystem S1 until it intersects at A2 with the nonsmooth boundary Σ.


#### 4.2.2. Supplementary Notes

It should be pointed out that although there are Hopf bifurcation pointsHB11 and HB21 on the equilibrium as shown in Figure 4, the trajectory still runs along the equilibrium due to the slow path effect which causes the bifurcation lag [25]. The phenomenon that the trajectory contacts the nonsmooth boundary and runs along Σ is called sliding motion. The sliding process is still in quiescent state until it reaches point A3. Besides, it can be clearly seen from Figure 3(a) that when trajectory is located on the nonsmooth boundary (sliding motions), the original system will degenerate into a two-dimensional system:(17)y˙=x−y+z,z˙=−βy,indicating that there will be a new equilibrium point (stable focus equilibrium point) on the interface x=0, named pseudoequilibrium point. With the variation of parameter β=R2C2/L, the pseudoequilibrium will also bifurcate.


#### 4.2.3. Attractor Structure When the Sliding Phenomenon Occurs

When the phase trajectory approaches the boundary from sideD−, it is attracted by EB61 located in D+ on subsystem S1, so it has the tendency to move across the boundary to D+. However, once the trajectory reaches D+, the control system will be switched to S2 immediately. The corresponding stable equilibrium EB62 is located at D−, which makes the trajectory return to D−. The alternating control of two subsystems is the intrinsic cause of the sliding phenomenon.


#### 4.2.4. Sliding Time Calculation

In our previous work [14], the sliding time has been estimated by the time required from FB11 to A3 with insufficient accuracy. In addition, this estimation method is not universal. New method are ought to be proposed, which will be discussed in the next section.


#### 4.2.5. Discussion on Physical Meaning

In general, fold, Hopf, and sliding bifurcation, respectively, correspond to mutation, oscillation, and conservation of circuit component voltage in the circuit system. A constant voltage transformer is widely used in engineering practice. The above results show that constant voltage control can be achieved by choosing appropriate nonlinear resistors and the output time of constant voltage can accurately be controlled by the calculation method of sliding time. The key point is how to deduce the analytical expression of sliding motions (constant voltage) which will be discussed in Section5.


### 4.3. Scenario 2:b=0.70

In this case, a problem which was not well resolved in [14] should be pointed out. The overlap of the transformed phase portrait with the equilibrium branches is shown in Figure 5(a), and the local enlargement of time history is shown in Figure 5(b).Figure 5
(a) Overlap of the transformed phase portrait with the equilibrium branches and (b) local enlargement of time history.
(a)
(b)Although the distribution of attractor structure in case 2 is different from that in case 1, most of its dynamic behavior can be explained on the basis of previous studies, so it is not necessary to elaborate on this for article length consideration. It is noteworthy that as shown in local enlargement at the lower right corner of Figure5(a), when the trajectory slides from point A2, a small jump occurs immediately and then returns to the interface to continue sliding. This seemingly insignificant phenomenon cannot be well explained by attractor structure analysis due to the coexistence of multiple equilibrium points. In order to further study the special dynamic behavior of nonsmooth systems at nonsmooth boundary, it is obviously not enough to use attractor structure analysis alone. Therefore, based on differential inclusion theory and focusing on the local dynamics near the interface, another method to analyze such problems named vector structure analysis is proposed.


## 4.1. TPP for the TrajectoryΠT:xt,yt,zt,t∈R

The traditional phase portrait reflects the relationship among state variables in the phase spacex,y,z or the projection planes x,y and y,z. If the slowly varying parameter is regarded as a generalized state variable, the generalized phase trajectory can be obtained ΠGT:xt,yt,zt,wt,t∈R, where wt=AsinΩt. The generalized phase trajectories in the generalized phase space x,y,z,w and its projected phase plane are called transformed phase portraits, which are used to characterize the relationship between different state variables and slowly varying parameters.Variable parametersb were taken b=0.14 and b=0.70, respectively, while other parameters were fixed as(14)α=6.9,β=10.0,a=−27,c=−2,A=6.0.Therefore,w varies periodically with frequency 0.01 in the interval −6.0,6.0. Taking A = 0, the natural frequency of the system can be calculated. For b=0.14, the equilibrium point of the two subsystems is approximately E±=±1.1662,0.0,∓1.1662, with its corresponding eigenvalue as follows:(15)λ1=−7.920,λ±=−0.0784±2.988I,indicating that the natural frequency is about 2.988. Similarly, when b=0.70, the coordinates of the three equilibrium points and their corresponding eigenvalues are(16)E1±=±0.4517,0.0,∓0.4517,λ1=4.5159,λ±=−0.9463±2.6698I,E2±=±1.2897,0.0,∓1.2897,λ1=−6.0386,λ±=0.01534±2.8797I,E3±=±1.7414,0.0,∓1.7417,λ1=−13.622,λ±=−0.2417±3.0923I.Therefore, when the trajectory oscillates near different equilibrium pointsEj±j=1,2,3, its natural frequencies are about 2.6698, 2.8797, and 3.0923.


## 4.2. Scenario 1:b=0.14

The 3D phase portrait and bifurcation diagram with parameterb are shown in Figure 3. It can be seen that the attractor structure keeps the symmetry of the original system. The alternation of large-amplitude spiking states and small-amplitude quiescent states can be observed. Specially, the sliding phenomenon appears on the nonsmooth boundary indicating that the voltage of the capacitor remained a constant during this period of time. It should be pointed out that the phenomenon in Figure 3(b) does not mean the changes of periodicity (similar to period doubling bifurcation) but that the position of equilibrium points migrates due to the occurrence of sliding bifurcation which will be discussed in detail.Figure 3
(a) 3D phase portrait withb = 0.14 and (b) bifurcation diagram with parameter b.
(a)
(b)In order to reveal the mechanism of bursting oscillation, the overlap of the transformed phase portrait [14] with the equilibrium branches can be observed in Figure 4. EBjii=1,2;j=1,2,…,6 corresponds to the equilibrium curve segments of subsystem Sii=1,2, respectively. The solid part of the equilibrium line represents the stable equilibrium state, and the dashed line represents the unstable state. The black lines represent the equilibrium points that can be found in the corresponding region, while the red parts denote the equilibrium state that exists but is not located in the phase space of the corresponding system in control.Figure 4
Overlap of the transformed phase portrait with the equilibrium branches.


### 4.2.1. Bursting Oscillation Mechanism

As shown in Figure4, the trajectory dominated by subsystem S1 starts from point A1 w=−6.0 and runs almost strictly along its stable equilibrium state EB11. Despite the existence of a stable limit cycle LC11 between the two supercritical Hopf bifurcation points HB11 and HB21, the trajectory still runs along the stable equilibrium, showing the quiescent state. The stability of limit cycles produced by these two Hopf bifurcations is opposite. They disappear after contact suggesting the occurrence of LPC bifurcation. When trajectory reaches the fold bifurcation point FB11, it jumps directly to the stable equilibrium. As time t continues to increase, the trajectory is attracted by the only stable equilibrium segment EB61 of subsystem S1 until it intersects at A2 with the nonsmooth boundary Σ.


### 4.2.2. Supplementary Notes

It should be pointed out that although there are Hopf bifurcation pointsHB11 and HB21 on the equilibrium as shown in Figure 4, the trajectory still runs along the equilibrium due to the slow path effect which causes the bifurcation lag [25]. The phenomenon that the trajectory contacts the nonsmooth boundary and runs along Σ is called sliding motion. The sliding process is still in quiescent state until it reaches point A3. Besides, it can be clearly seen from Figure 3(a) that when trajectory is located on the nonsmooth boundary (sliding motions), the original system will degenerate into a two-dimensional system:(17)y˙=x−y+z,z˙=−βy,indicating that there will be a new equilibrium point (stable focus equilibrium point) on the interface x=0, named pseudoequilibrium point. With the variation of parameter β=R2C2/L, the pseudoequilibrium will also bifurcate.


### 4.2.3. Attractor Structure When the Sliding Phenomenon Occurs

When the phase trajectory approaches the boundary from sideD−, it is attracted by EB61 located in D+ on subsystem S1, so it has the tendency to move across the boundary to D+. However, once the trajectory reaches D+, the control system will be switched to S2 immediately. The corresponding stable equilibrium EB62 is located at D−, which makes the trajectory return to D−. The alternating control of two subsystems is the intrinsic cause of the sliding phenomenon.


### 4.2.4. Sliding Time Calculation

In our previous work [14], the sliding time has been estimated by the time required from FB11 to A3 with insufficient accuracy. In addition, this estimation method is not universal. New method are ought to be proposed, which will be discussed in the next section.


### 4.2.5. Discussion on Physical Meaning

In general, fold, Hopf, and sliding bifurcation, respectively, correspond to mutation, oscillation, and conservation of circuit component voltage in the circuit system. A constant voltage transformer is widely used in engineering practice. The above results show that constant voltage control can be achieved by choosing appropriate nonlinear resistors and the output time of constant voltage can accurately be controlled by the calculation method of sliding time. The key point is how to deduce the analytical expression of sliding motions (constant voltage) which will be discussed in Section5.


## 4.2.1. Bursting Oscillation Mechanism

As shown in Figure4, the trajectory dominated by subsystem S1 starts from point A1 w=−6.0 and runs almost strictly along its stable equilibrium state EB11. Despite the existence of a stable limit cycle LC11 between the two supercritical Hopf bifurcation points HB11 and HB21, the trajectory still runs along the stable equilibrium, showing the quiescent state. The stability of limit cycles produced by these two Hopf bifurcations is opposite. They disappear after contact suggesting the occurrence of LPC bifurcation. When trajectory reaches the fold bifurcation point FB11, it jumps directly to the stable equilibrium. As time t continues to increase, the trajectory is attracted by the only stable equilibrium segment EB61 of subsystem S1 until it intersects at A2 with the nonsmooth boundary Σ.


## 4.2.2. Supplementary Notes

It should be pointed out that although there are Hopf bifurcation pointsHB11 and HB21 on the equilibrium as shown in Figure 4, the trajectory still runs along the equilibrium due to the slow path effect which causes the bifurcation lag [25]. The phenomenon that the trajectory contacts the nonsmooth boundary and runs along Σ is called sliding motion. The sliding process is still in quiescent state until it reaches point A3. Besides, it can be clearly seen from Figure 3(a) that when trajectory is located on the nonsmooth boundary (sliding motions), the original system will degenerate into a two-dimensional system:(17)y˙=x−y+z,z˙=−βy,indicating that there will be a new equilibrium point (stable focus equilibrium point) on the interface x=0, named pseudoequilibrium point. With the variation of parameter β=R2C2/L, the pseudoequilibrium will also bifurcate.


## 4.2.3. Attractor Structure When the Sliding Phenomenon Occurs

When the phase trajectory approaches the boundary from sideD−, it is attracted by EB61 located in D+ on subsystem S1, so it has the tendency to move across the boundary to D+. However, once the trajectory reaches D+, the control system will be switched to S2 immediately. The corresponding stable equilibrium EB62 is located at D−, which makes the trajectory return to D−. The alternating control of two subsystems is the intrinsic cause of the sliding phenomenon.


## 4.2.4. Sliding Time Calculation

In our previous work [14], the sliding time has been estimated by the time required from FB11 to A3 with insufficient accuracy. In addition, this estimation method is not universal. New method are ought to be proposed, which will be discussed in the next section.


## 4.2.5. Discussion on Physical Meaning

In general, fold, Hopf, and sliding bifurcation, respectively, correspond to mutation, oscillation, and conservation of circuit component voltage in the circuit system. A constant voltage transformer is widely used in engineering practice. The above results show that constant voltage control can be achieved by choosing appropriate nonlinear resistors and the output time of constant voltage can accurately be controlled by the calculation method of sliding time. The key point is how to deduce the analytical expression of sliding motions (constant voltage) which will be discussed in Section5.


## 4.3. Scenario 2:b=0.70

In this case, a problem which was not well resolved in [14] should be pointed out. The overlap of the transformed phase portrait with the equilibrium branches is shown in Figure 5(a), and the local enlargement of time history is shown in Figure 5(b).Figure 5
(a) Overlap of the transformed phase portrait with the equilibrium branches and (b) local enlargement of time history.
(a)
(b)Although the distribution of attractor structure in case 2 is different from that in case 1, most of its dynamic behavior can be explained on the basis of previous studies, so it is not necessary to elaborate on this for article length consideration. It is noteworthy that as shown in local enlargement at the lower right corner of Figure5(a), when the trajectory slides from point A2, a small jump occurs immediately and then returns to the interface to continue sliding. This seemingly insignificant phenomenon cannot be well explained by attractor structure analysis due to the coexistence of multiple equilibrium points. In order to further study the special dynamic behavior of nonsmooth systems at nonsmooth boundary, it is obviously not enough to use attractor structure analysis alone. Therefore, based on differential inclusion theory and focusing on the local dynamics near the interface, another method to analyze such problems named vector structure analysis is proposed.


## 5. Structural Characteristics of Vector Field for Sliding Bifurcation

Consider a generaln-dimensional piecewise smooth system:(18)Y˙=F1Y,hY>0,F2Y,hY<0.The nonsmooth boundary can be defined asΣ=Y∈Rn+1hY=0. The vector field of the sliding region can be represented by convex closures of F1Y and F2Y:(19)FsY=F1Y+F2Y2+F1Y−F2Y2αY,where −1≤αY≤1. Assuming that hY=∂h/∂y1,∂h/∂y2,…T=∇h, so hYFSY=0 (normal vector hY is perpendicular to the nonsmooth interface, while hYFSY=0 is parallel to the interface). Therefore, αY can be expressed as(20)αY=−hYF1+F2hYF2−F1.The sliding region and its boundary can be expressed as(21)Σ⌢=Y∈Σ−1≤αY≤1,(22)∂Σ⌢±=Y∈ΣαY=±1,where ∂Σ⌢− is the tangent point of F1Y on nonsmooth boundary Σ. So, hYF1Y=0 and hYF2Y=0. By deriving the two sides of formula (20) separately,(23)αY=−hYF1+F2YhYF2−F1hYF2−F12+hYF1+F2hYF2−F1YhYF2−F12=−hYF1Y+hYYF1+hYF2Y+hYYF2hYF2−hYF1hYF2−hYF12+hYF2Y+hYYF2−hYF1Y−hYYF1hYF1+hYF2hYF2−hYF12=−hYF1Y+hYF2YhYF2−hYF1hYF2−hYF12+hYF2Y−hYF1YhYF1+hYF2hYF2−hYF12,∂Σ⌢− satisfying FSY=F1Y, indicating that hYF1=0. Combining with (23), normal vectors on ∂Σ⌢− can be expressed as(24)αY=−2hYF2hYF1Y.The physical meaning of differential inclusion theory can be understood simply through dry friction: sliding friction corresponds to the boundary of sliding region, static friction corresponds to the interior of sliding region, and the numerical value of static friction depends on the equilibrium condition.


### 5.1. Geometric Classification of Sliding Bifurcations

According to the geometric structure of the periodic solution trajectory intersecting with the sliding region and its boundary, the sliding bifurcation can be roughly divided into the following four categories, as shown in Figure6.Figure 6
(a) Crossing, (b) grazing, (c) switching, and (d) adding sliding bifurcations.
(a)
(b)
(c)
(d)The characteristic of crossing sliding bifurcation is that the trajectory crosses the boundary of the sliding region at the bifurcation point and leaves (as shown in Figure6(a) trajectory b). Then the sliding phenomena occur when the contact points enter the sliding area with the change of parameters (as shown in the trajectory c of Figure 6(a)).The grazing sliding bifurcation, as the name implies, grazes along the boundary of the sliding region at the bifurcation point and leaves (as shown in Figure6(b) track b). Then, the sliding phenomena occur when the contact points enter the sliding area with the change of parameters (as shown in the trajectory c of Figure 6(b)).The switching sliding bifurcation is similar to crossing sliding bifurcation, but the difference is that a section of track is still in the sliding region after traversing the boundary of the sliding region at the bifurcation point (as shown in Figure6(c) trajectory b). Then, a segment of the trajectory remains in the sliding region as the parameters change (as shown in Figure 6(c) trajectory c).Unlike the other three bifurcations, the bifurcated trajectories of the adding sliding bifurcation are completely located in the sliding region and tangent to the boundary of the sliding region at the bifurcation point (as shown in track b of Figure6(d)). As the parameters change, a section of the track is still located in the sliding region (as shown in track c of Figure 6(d)).


### 5.2. Analytical Conditions for Sliding Bifurcation

The above four kinds of sliding bifurcations can accurately be described by analytic conditions.


#### 5.2.1. Crossing Sliding Bifurcation and Grazing Sliding Bifurcation

Because the bifurcation pointY∗ is located on the nonsmooth boundary,(25)hY∗=0,hYY∗≠0.Due toY∗∈∂Σ⌢−,(26)ddthϕY∗,tt=0=hYF1Y∗=0,where ϕ is the corresponding smooth solution manifold. In addition, after crossing the interface, the trajectory moves away from the sliding area, so(27)ddtαϕY∗,tt=0=αYF1Y∗<0.By the law of derivation of the dot product and the combination of formulas (24), (26), and (27),(28)d2dt2hϕY∗,tt=0=hYF1YF1Y∗>0.


#### 5.2.2. Switching Sliding Bifurcation

This situation also satisfies the two basic conditions (25) and (26), but the track crosses the sliding area at the bifurcation point and remains in the sliding area; therefore,(29)ddtαϕY∗,tt=0=αYF1Y∗>0,(30)d2dt2hϕY∗,tt=0=hYF1YF1Y∗<0.


#### 5.2.3. Adding Sliding Bifurcation

The two basic conditions (25) and (26) are still valid, but in this case, the bifurcation trajectory is completely located within the sliding region and tangent to the boundary of the sliding region at the bifurcation point, so the conditions (28) and (30) for the above three cases degenerate:(31)d2dt2hϕY∗,tt=0=hYF1YF1Y∗=0.Because the bifurcation trajectory has a minimum value at∂Σ⌢−, there is a nondegeneration condition:(32)d2dt2αϕY∗,tt=0=αYF1YF1Y∗>0.By the law of derivation of the dot product and the combination of formulas (24), (26), and (32),(33)d3dt3hϕY∗,tt=0=hYF1Y2F1Y∗<0.


### 5.3. Vector Field Structure Characteristics of Sliding Bifurcation

Based on the previous analysis, the vector field structure of sliding bifurcation can be summarized.Regardless of the type of sliding bifurcation, as long as the trajectory is located in the sliding region, its vector field must have the structural characteristics as shown in Figure7(a).Figure 7
The structure of the vector field in (a) sliding region and (b) traversing region.
(a)
(b)Sliding phenomena occur when there are opposite vector fields on both sides of the nonsmooth interface. The direct reason is that the trajectory is forced by the opposite vector field near the boundary and thus can only oscillate along the boundary. When the limit of step size approaches zero, it is a strict sliding motion. With opposite vector field in Figure7(a), if trajectory locates in D−, it will be pushed to nonsmooth boundary. For the same reason, if it locates in D+, it will also approach the boundary. Therefore, it can be understood as that the trajectory is squeezed by the vector field and can only strictly slide along the nonsmooth boundary. This result is consistent with the conclusions of attractor structure analysis. In contrast, when vector fields with the same sign exist on both sides of the nonsmooth boundary, the trajectory traverses the interface directly, as shown in Figure 7(b).In scenario 2, the transient jumping phenomenon after sliding is caused by the change of external excitation in the sliding region, which leads to the change of vector field, until the vector field on one side weakens gradually and eventually reverses. It should be pointed out that although there is only one stable attractor on both sides of the interface, the complete sliding phenomenon may not occur. This is because the stability of the attractor does not directly determine the direction of the normal component of the vector field at the interface, as shown in diagrammatic sketch Figure8.Figure 8
The vector field of a stable focus.For instance, supposing that subsystem defined inD+ has a stable focus (or limit cycle) P+ in D−. The direction of the vector field in its domain of attraction is shown by the blue arrows. It can be clearly seen that although P+ is attractive to the phase trajectory in D+ in the long run, its vector field on the side of D+ does not point to the nonsmooth boundary (as shown in the dotted part). Therefore, although the attractor structure analysis methods (such as fast-slow analysis method) have some advantages in studying the global dynamic behaviors of nonsmooth periodic excitation system, it also has some limitations in dealing with some special local dynamic details such as sliding bifurcation.


### 5.4. Numerical Verification

In order to better understand the mechanism of sliding bifurcation, the sliding process in case 2 is numerically analyzed as Table1.Table 1
Simulation data near points of special behaviors.
w
x
y
x
˙
A
2 (jumping)
1.0085
≈0, >0
0.1413
≈0, <0
1.0091
≈0, >0
0.1436
≈0, >0
…
…
…
…
…<0
A
2
⟶
A
3 (sliding)
1.5001
≈0, >0
0.00285
−0.4769 < 0
1.5006
≈0, <0
0.00281
+3.5230 > 0It is noteworthy that in sliding regionA2⟶A3, the signs of vector field x˙ on both sides of nonsmooth boundary x=0 are opposite and the difference is exactly equal to 2c=2×2.0=4.0 (c is the coefficient of the symbolic function shown in formula (14)). Based on the analytical expressions (20)–(22) of the sliding region of the upper section, it can be found that the numerical results are in good agreement with the theoretical results.


### 5.5. Discussion on Physical Meaning

The proposed analytical conditions of sliding bifurcation imply that the constant voltage can be controlled by changing the external excitation (AC source) rather than the properties of the circuit elements themselves. Besides, the sliding phenomenon also exists in mechanical systems such as spring oscillator considering dry friction factor, etc. Actually, in engineering practice, the sliding (sticking) phenomenon is not always expected. For instance, the occurrence of sliding phenomenon of the lathes for producing models with spring oscillator means the production of defective products. Therefore, finding out a way to eliminate the sliding motions is also of great concern. Controlling the frequency of external excitation instead of changing the property of the material itself is a relatively simple and effective solution because the change of external excitation invalidates the necessary condition of sliding bifurcation mentioned in the manuscript. This part of the content will be further explored in our next work.


## 5.1. Geometric Classification of Sliding Bifurcations

According to the geometric structure of the periodic solution trajectory intersecting with the sliding region and its boundary, the sliding bifurcation can be roughly divided into the following four categories, as shown in Figure6.Figure 6
(a) Crossing, (b) grazing, (c) switching, and (d) adding sliding bifurcations.
(a)
(b)
(c)
(d)The characteristic of crossing sliding bifurcation is that the trajectory crosses the boundary of the sliding region at the bifurcation point and leaves (as shown in Figure6(a) trajectory b). Then the sliding phenomena occur when the contact points enter the sliding area with the change of parameters (as shown in the trajectory c of Figure 6(a)).The grazing sliding bifurcation, as the name implies, grazes along the boundary of the sliding region at the bifurcation point and leaves (as shown in Figure6(b) track b). Then, the sliding phenomena occur when the contact points enter the sliding area with the change of parameters (as shown in the trajectory c of Figure 6(b)).The switching sliding bifurcation is similar to crossing sliding bifurcation, but the difference is that a section of track is still in the sliding region after traversing the boundary of the sliding region at the bifurcation point (as shown in Figure6(c) trajectory b). Then, a segment of the trajectory remains in the sliding region as the parameters change (as shown in Figure 6(c) trajectory c).Unlike the other three bifurcations, the bifurcated trajectories of the adding sliding bifurcation are completely located in the sliding region and tangent to the boundary of the sliding region at the bifurcation point (as shown in track b of Figure6(d)). As the parameters change, a section of the track is still located in the sliding region (as shown in track c of Figure 6(d)).


## 5.2. Analytical Conditions for Sliding Bifurcation

The above four kinds of sliding bifurcations can accurately be described by analytic conditions.


### 5.2.1. Crossing Sliding Bifurcation and Grazing Sliding Bifurcation

Because the bifurcation pointY∗ is located on the nonsmooth boundary,(25)hY∗=0,hYY∗≠0.Due toY∗∈∂Σ⌢−,(26)ddthϕY∗,tt=0=hYF1Y∗=0,where ϕ is the corresponding smooth solution manifold. In addition, after crossing the interface, the trajectory moves away from the sliding area, so(27)ddtαϕY∗,tt=0=αYF1Y∗<0.By the law of derivation of the dot product and the combination of formulas (24), (26), and (27),(28)d2dt2hϕY∗,tt=0=hYF1YF1Y∗>0.


### 5.2.2. Switching Sliding Bifurcation

This situation also satisfies the two basic conditions (25) and (26), but the track crosses the sliding area at the bifurcation point and remains in the sliding area; therefore,(29)ddtαϕY∗,tt=0=αYF1Y∗>0,(30)d2dt2hϕY∗,tt=0=hYF1YF1Y∗<0.


### 5.2.3. Adding Sliding Bifurcation

The two basic conditions (25) and (26) are still valid, but in this case, the bifurcation trajectory is completely located within the sliding region and tangent to the boundary of the sliding region at the bifurcation point, so the conditions (28) and (30) for the above three cases degenerate:(31)d2dt2hϕY∗,tt=0=hYF1YF1Y∗=0.Because the bifurcation trajectory has a minimum value at∂Σ⌢−, there is a nondegeneration condition:(32)d2dt2αϕY∗,tt=0=αYF1YF1Y∗>0.By the law of derivation of the dot product and the combination of formulas (24), (26), and (32),(33)d3dt3hϕY∗,tt=0=hYF1Y2F1Y∗<0.


## 5.2.1. Crossing Sliding Bifurcation and Grazing Sliding Bifurcation

Because the bifurcation pointY∗ is located on the nonsmooth boundary,(25)hY∗=0,hYY∗≠0.Due toY∗∈∂Σ⌢−,(26)ddthϕY∗,tt=0=hYF1Y∗=0,where ϕ is the corresponding smooth solution manifold. In addition, after crossing the interface, the trajectory moves away from the sliding area, so(27)ddtαϕY∗,tt=0=αYF1Y∗<0.By the law of derivation of the dot product and the combination of formulas (24), (26), and (27),(28)d2dt2hϕY∗,tt=0=hYF1YF1Y∗>0.


## 5.2.2. Switching Sliding Bifurcation

This situation also satisfies the two basic conditions (25) and (26), but the track crosses the sliding area at the bifurcation point and remains in the sliding area; therefore,(29)ddtαϕY∗,tt=0=αYF1Y∗>0,(30)d2dt2hϕY∗,tt=0=hYF1YF1Y∗<0.


## 5.2.3. Adding Sliding Bifurcation

The two basic conditions (25) and (26) are still valid, but in this case, the bifurcation trajectory is completely located within the sliding region and tangent to the boundary of the sliding region at the bifurcation point, so the conditions (28) and (30) for the above three cases degenerate:(31)d2dt2hϕY∗,tt=0=hYF1YF1Y∗=0.Because the bifurcation trajectory has a minimum value at∂Σ⌢−, there is a nondegeneration condition:(32)d2dt2αϕY∗,tt=0=αYF1YF1Y∗>0.By the law of derivation of the dot product and the combination of formulas (24), (26), and (32),(33)d3dt3hϕY∗,tt=0=hYF1Y2F1Y∗<0.


## 5.3. Vector Field Structure Characteristics of Sliding Bifurcation

Based on the previous analysis, the vector field structure of sliding bifurcation can be summarized.Regardless of the type of sliding bifurcation, as long as the trajectory is located in the sliding region, its vector field must have the structural characteristics as shown in Figure7(a).Figure 7
The structure of the vector field in (a) sliding region and (b) traversing region.
(a)
(b)Sliding phenomena occur when there are opposite vector fields on both sides of the nonsmooth interface. The direct reason is that the trajectory is forced by the opposite vector field near the boundary and thus can only oscillate along the boundary. When the limit of step size approaches zero, it is a strict sliding motion. With opposite vector field in Figure7(a), if trajectory locates in D−, it will be pushed to nonsmooth boundary. For the same reason, if it locates in D+, it will also approach the boundary. Therefore, it can be understood as that the trajectory is squeezed by the vector field and can only strictly slide along the nonsmooth boundary. This result is consistent with the conclusions of attractor structure analysis. In contrast, when vector fields with the same sign exist on both sides of the nonsmooth boundary, the trajectory traverses the interface directly, as shown in Figure 7(b).In scenario 2, the transient jumping phenomenon after sliding is caused by the change of external excitation in the sliding region, which leads to the change of vector field, until the vector field on one side weakens gradually and eventually reverses. It should be pointed out that although there is only one stable attractor on both sides of the interface, the complete sliding phenomenon may not occur. This is because the stability of the attractor does not directly determine the direction of the normal component of the vector field at the interface, as shown in diagrammatic sketch Figure8.Figure 8
The vector field of a stable focus.For instance, supposing that subsystem defined inD+ has a stable focus (or limit cycle) P+ in D−. The direction of the vector field in its domain of attraction is shown by the blue arrows. It can be clearly seen that although P+ is attractive to the phase trajectory in D+ in the long run, its vector field on the side of D+ does not point to the nonsmooth boundary (as shown in the dotted part). Therefore, although the attractor structure analysis methods (such as fast-slow analysis method) have some advantages in studying the global dynamic behaviors of nonsmooth periodic excitation system, it also has some limitations in dealing with some special local dynamic details such as sliding bifurcation.


## 5.4. Numerical Verification

In order to better understand the mechanism of sliding bifurcation, the sliding process in case 2 is numerically analyzed as Table1.Table 1
Simulation data near points of special behaviors.
w
x
y
x
˙
A
2 (jumping)
1.0085
≈0, >0
0.1413
≈0, <0
1.0091
≈0, >0
0.1436
≈0, >0
…
…
…
…
…<0
A
2
⟶
A
3 (sliding)
1.5001
≈0, >0
0.00285
−0.4769 < 0
1.5006
≈0, <0
0.00281
+3.5230 > 0It is noteworthy that in sliding regionA2⟶A3, the signs of vector field x˙ on both sides of nonsmooth boundary x=0 are opposite and the difference is exactly equal to 2c=2×2.0=4.0 (c is the coefficient of the symbolic function shown in formula (14)). Based on the analytical expressions (20)–(22) of the sliding region of the upper section, it can be found that the numerical results are in good agreement with the theoretical results.


## 5.5. Discussion on Physical Meaning

The proposed analytical conditions of sliding bifurcation imply that the constant voltage can be controlled by changing the external excitation (AC source) rather than the properties of the circuit elements themselves. Besides, the sliding phenomenon also exists in mechanical systems such as spring oscillator considering dry friction factor, etc. Actually, in engineering practice, the sliding (sticking) phenomenon is not always expected. For instance, the occurrence of sliding phenomenon of the lathes for producing models with spring oscillator means the production of defective products. Therefore, finding out a way to eliminate the sliding motions is also of great concern. Controlling the frequency of external excitation instead of changing the property of the material itself is a relatively simple and effective solution because the change of external excitation invalidates the necessary condition of sliding bifurcation mentioned in the manuscript. This part of the content will be further explored in our next work.


## 6. Conclusions

The introduction of bidirectional diodes in the circuit system will result in the occurrence of sliding bifurcation indicating the conservation of voltage. The selection of different nonlinear resistors corresponds to different sliding time and modes. Sliding motions can be divided into four categories according to the corresponding geometric characteristics. The specific expression of sliding area and its boundary along with the special structure of vector field are presented. Vector fields on both sides of nonsmooth boundary are opposite in sliding region while remain the same in traversing region.In addition, this paper mainly focuses on the single variable of nonlinear resistance. Coupling effects between multiple factors will be discussed separately in other literature.



---

*Source: 8213808-2019-07-18.xml*