# Bifurcation Analysis of Stick-Slip Motion of the Vibration-Driven System with Dry Friction

Authors: Peng Li; Ziwang Jiang
Journal: Mathematical Problems in Engineering
 (2018)

---

## Abstract

This paper is concerned with the vibration-driven system which can move due to the periodic motion of the internal mass and the dry friction; the system can be modeled as Filippov system and has the property of stick-slip motion. Different periodic solutions of stick-slip motion can be analyzed through sliding bifurcation, two-parameter numerical continuation for sliding bifurcation is carried out to get the different bifurcation curves, and the bifurcation curves divide the parameters plane into different regions which stand for different stick-slip motion of the periodic solution. Furthermore, continuations with additional conditionv=0 are carried out for the directional control of the vibration-driven system in one period; the curves divide the parameter plane into different progressions.


---

## Body

## 1. Introduction

Recently, mobile mechanisms that can move due to the vibration of the internal mass have been widely researched, and these mechanisms have many advantages over conventional mobile systems (driven by legs, wheels, wings, etc.), for example, easy fabrication, hermetic structure, and locomotion in the narrow environment. Thus they have extensive application in pipeline inspection, life detection in disaster, and medical endoscopy.Chernousko [1] first proposed the horizontal motion of the system driven by the movable internal mass; the friction which acted on the body is anisotropic, which means the coefficient of friction in forward and backward direction is different. The two periodic control modes, velocity-controlled mode and acceleration-controlled mode, are constructed for the relative motion of the internal mass, and optimal parameters of periodic control were decided to realize the maximum mean velocity of the body. Fang et al. [2] used the method of averaging to obtain an approximate expression of the average steady-state velocity when the stick-slip phenomenon was not considered, optimal parameters of the internal controlled mass were determined to maximize the average velocity, and some control strategies were given to control the motion of system under the stick-slip effect. Liu et al. [3] studied the vibroimpact capsule system which has a main body interacting with an internal harmonically driven mass, when the internal mass contact with the plate impact occurs, and the parameters for the maximum mean velocity can be determined through nonlinear dynamics analysis, the energy consumption was also considered, and the parameters for the maximum mean velocity and the minimum energy consumption were not the same. Fang et al. [4, 5] and Zimmermann et al. [6–8] studied the two and more modules vibration-driven systems; the approximate expression of steady-state motion was obtained when the friction is small and the optimal parameters were got to achieve the maximum mean velocity. Bolotnik et al. [9] modeled the system driven by the movable internal mass which can move in the horizontal direction and the vertical direction (change the normal force for anisotropic friction in the different direction). Then the approximate expression of average steady-state velocity was obtained through the method of averaging; optimal parameters (the amplitude and the phase shift of the horizontal and vertical vibration excitation forces) were determined to realize the maximum average velocity and to control the direction of motion. In the paper, we study the model.The dry friction plays an important role in vibration-driven system motion. The systems with dry friction belong to Filippov piecewise-smooth dynamical systems [10, 11]. The Filippov systems may exhibit different types of limit cycles caused by the interaction of a trajectory with the boundary of the sliding regions; the features of Filippov system are called the sliding bifurcation. Kowalczyk et al. [12] investigated a dry friction oscillator through numerical continuation of sliding bifurcation and revealed the codimension two sliding bifurcation points. Marcel Guardia et al. [13] analytically considered sliding bifurcations of periodic orbits in the dry friction oscillator, and the results agreed with the numerical calculation [12]. Fang et al. [14] studied the vibration-driven system through sliding bifurcation, and a two-parameter bifurcation problem was theoretically analyzed. For the numerical continuation of piecewise-smooth system, the software SlideCont [15] and TC-HAT [16] based on ATUO have been developed. Joseph Páez Chávez used the software TC-HAT to study the bifurcation of some mechanical models, the nonsmooth Jeffcott rotor [17], the impact oscillator [18], the piecewise-linear capsule system [19], etc. The continuation toolbox COCO was developed for continuation and bifurcation analysis of smooth and nonsmooth dynamical systems [20, 21]; the soft impact oscillator [22] and the impulsively coupled oscillators [23] were analyzed through the COCO. In this paper, the COCO will be employed to numerically study the sliding bifurcation of vibration-driven system; the sliding bifurcation will help us to understand the stick-slip property of periodic solution and give some instructive ideas to design and control the system.The paper is organized as follows. In Section2, the model of the vibration-driven system is described. The mathematical model of vibration-driven system is studied in detail in order to perform the numerical analysis by the mean of COCO in Section 3. Two-parameter sliding bifurcations are analyzed and the directional control of the vibration-driven system is tackled by the numerical continuation in Section 4, and some conclusions are given in Section 5.


## 2. Modeling of Vibration-Driven System

The vibration-driven system is considered as depicted in Figure1. The vibration-driven system is composed of a rigid body and an internal mass; the rigid body realizes the translational motion along a straight line in the resistive environment. The internal mass can move relative to the rigid body in horizontal and vertical direction. The internal mass is considered a point mass. Dry friction acts between the rigid body and the ground.Figure 1
Vibration-driven system.Two Cartesian reference frames are introduced: the inertial reference frameOxy and the coordinate system O′ξη attached to the rigid body. The x- and ξ-axes are horizontal; the y- and η-axes are vertical. We denote x as the coordinate of the point O′ in the inertial reference frame Oxy, x denotes the displacement of the rigid body, ξ and η denote the coordinates of the internal body in the reference frame O′ξη, m and m1 are the masses of the rigid body and the internal body, respectively, and g is the gravitational acceleration.The dynamics equation of the system alongx-axis can be governed by Newton’s second laws as follows:(1)m+m1x¨+m1ξ¨=Ff,where Ff is the frictional force. The force Ff is described by Coulomb law: (2)Ff=-μFNsgnx˙,ifx˙≠0,-F0,ifx˙=0&F0≤μFN,-μFNsgnF0,ifx˙=0&F0>μFN,where F0 is the resultant force on the body except for the dry friction in the horizontal direction, FN is the normal force exerted on the system by the ground, and μ is the coefficient of dry friction. The asymmetrical friction in forward and backward motions arises from the time-varying normal force FN. The force F0 and force FN can be expressed: (3)F0=-m1ξ¨,(4)FN=m+m1g+m1η¨.The rigid body keeps the contact with the ground; therefore, the force FN must be satisfied the inequality FN≥0. So we have(5)m+m1g+m1η¨≥0.Next, we assume the control acceleration is harmonic function with the same frequency but shifted in the phase; particularly,(6)ξ¨=-Asin⁡ωt,η¨=Bsin⁡ωt+ϕ.HereA and B are the driving-amplitudes, respectively, ω is the driving-frequency, and ϕ is the phase difference between the forces.Substituting (6) into (1)-(5),(7)Mx¨=Fxsin⁡ωt+Ff,(8)Ff=-μFNsgnx˙,ifx˙≠0,-Fxsin⁡ωt,ifx˙=0&Fxsin⁡ωt≤μFN,-μFNsgnsin⁡ωt,ifx˙=0&Fxsin⁡ωt>μFN,(9)FN=Mg+Fysin⁡ωt+ϕ,FN≥0,where M=m+m1, Fx=m1A, and Fy=m1B.To reduce the number of parameters of the system, the nondimensional variablesx∗ and t∗ and the parameters ε and α are introduced:(10)x∗=Mω2Fxx,t∗=ωt,ε=μMgFx,α=FyMgSubstituting these variables above into (7)-(9) (omit the asterisks),(11)x¨=sin⁡t+fc,(12)fc=-εfnsgnx˙,ifx˙≠0,-sin⁡t,ifx˙=0&sin⁡t≤εfn,-εfnsgnsin⁡t,ifx˙=0&sin⁡t>εfn,(13)fn=1+αsin⁡t+ϕ,α≤1.The expressions fc and fn stand for the normalized dry friction and the normal force, respectively. The parameter ε represents the ratio of the possible maximum value of the dry friction force to the amplitude Fx. We assume the value of the parameter ε is in the interval (0,1). Due to periodicity, the phase difference ϕ ranges from 0 to 2π.


## 3. Modeling of the Vibration-Driven System as Filippov System

The mathematical model of vibration-driven system can be defined as a piecewise-smooth system of the Filippov type. We can transform (11) into vector fields, event functions, and reset functions through the approach of multisegment periodic orbits. Let y=(y1,y2)T=(x˙,t)T and p=(ε,α,ϕ)T represent the state variables of the system and the parameters, respectively.The multisegment periodic orbits of the vibration-driven system consist of two or more segments, which can be modeled as follows.Stick: This segment occurs whensin⁡y2<εfn; the motion during this segment is governed by the equation(14)y˙=fs=01.This segment terminates when the resultant force on the body except for the dry friction equals the threshold of the dry friction force. The event functions are detected as follows:hf(y)=sin⁡y2-ε(1+αsin⁡(y2+ϕ))=0 (transition to forward slip)hb(y)=sin⁡y2+ε(1+αsin⁡(y2+ϕ))=0 (transition to backward slip)The next segment initial point is defined by the reset functiongid(y)=y.Forward slip: this segment occurs when the force is larger than the maximum value of the dry friction; that is,siny2>ε(1+αsin⁡(y2+ϕ)); the motion during this segment is governed by the equation(15)y˙=ff=sin⁡y2-ε1+αsin⁡y2+ϕ1.This segment ends when the velocity becomes zero; that is,h(y)=y1=0. The next segment is connected by the reset function gid(y)=y.Backward slip: this segment occurs whensin⁡y2<-ε(1+αsin⁡(y2+ϕ)); the motion of the system during this segments is described by the equation(16)y˙=fb=sin⁡y2+ε1+αsin⁡y2+ϕ1.This segment ends when the velocity becomes zero; the event function is defined:h(y)=y1=0. The next segment is connected by the reset function gid(y)=y.Stick2π: this segment is introduced to keep the variable y2 within the interval [0,2π), when sin⁡y2<εfn, the motion of the system is governed by (14), the segment ends when h2π(y)=2π-y2=0, and the reset function is(17)g2πy=y1y2-2π.Forward slip2π: this segment is introduced to keep the variable y2 within the interval [0,2π), when siny2>ε(1+αsin⁡(y2+ϕ)), the motion of the system is governed by (15), the segment ends when h2π(y)=2π-y2=0, and the reset function is g2π(y).Backward slip2π: this segment is introduced to keep the variable y2 within the interval[0,2π), when sin⁡y2<-ε(1+αsin⁡(y2+ϕ)), the motion of the system is governed by (16), the segment ends when h2π(y)=2π-y2=0, and the reset function is g2π(y).A periodic solution of Filippov system can be described as a sequence of tripletIi=(fi,hi,gi); the segment of system is governed by the vector field fi, terminates at the event function hi, and connects the next segment by the reset function gi. Any periodic trajectory of the system is described by solution signature {Ii}i=1M; M is the length of signature. Therefore the periodic solution of the vibration-driven can be described by combinations of the seven triplets corresponding to above statement:(18)I1=fs,hf,gid,I2=fs,hb,gid,I3=ff,h,gid,I4=fb,h,gid,I5=fs,h2π,g2π,I6=ff,h2π,g2π,I7=fb,h2π,g2π.In the Filippov system, there are four possible sliding bifurcations in the limit cycle because of the interaction of a trajectory with the boundary of a sliding region, including crossing-sliding bifurcation, gazing-sliding bifurcation, switching-sliding bifurcation, and adding-sliding bifurcation. Nondegeneracy conditions for the four sliding bifurcation are given [10, 11]. The sliding bifurcation does not change the number and stability of the system’s solutions, but it will induce the different interaction between the limit cycle and the sliding regions.


## 4. Numerical Bifurcation Analysis

In this section, we will perform the sliding bifurcation and directional control to analyze the dynamics response of the vibration-driven system.


### 4.1. Sliding Bifurcation Analysis

Whenε=0.6,α=0.456,ϕ=2π/3, the periodic trajectory of the motion is shown by Figure 2. The cycle signature is {I1,I3,I4,I5}. We start the numerical continuation of the periodic solution by the method of path-following using the parameter value as an initial value. The additional boundary hb(y)=sin⁡y2+ε(1+αsin⁡(y2+ϕ))=0 may be applied to the start point of the third segment for the crossing-sliding bifurcation continuation. We will use the COCO to carry out the numerical continuation concerning parameters ϕ and ε; the curve cs1 which is the result of this numerical continuation is shown in Figure 3.Figure 2
Periodic solution of vibration-driven system (11) computed for the parameter values ε=0.6,α=0.456,ϕ=2π/3. The trajectory consists of the segments I1,I3,I4,I5.Figure 3
Two-parameter continuation of sliding bifurcation with respect toϕ and ε for α=0.456. adi represent adding-sliding bifurcations, csi represent crossing-sliding bifurcations, and the regions Ei denote the different stick-slip motions of the periodic solutions.Similarly, the numerical continuation is performed in parameter spaceϕ×ε with different cycle signatures and different additional boundary condition for different segment boundary point and the results are depicted in Figure 3. The signature of cycle trajectory in E1 is {I1,I3,I5}, when the additional conditions hb(y)=sin⁡y2+ε(1+αsin⁡(y2+ϕ))=0 and cos⁡y2+εαcos⁡(y2+ϕ)=0 are applied to the point of the segment I5 according to nondegeneracy conditions of adding-sliding bifurcation; hence we can get the adding-sliding bifurcation curve ad1. The adding-sliding bifurcation branch ad2 is got similarly. The curves cs1, cs2, cs3, and cs4 represent crossing-sliding bifurcation branches. The sliding bifurcation curves divide the two parameters plane into eight regions and there are six different stick-slip periodic solutions in the parameter plane. The periodic solution of system in E1 with signature {I1,I3,I5} is depicted in Figure 4(a), which means the velocity of the system is always greater than or equal to 0. The signature of the periodic trajectory in E2 is {I4,I2,I7} as depicted in Figure 4(b), which means the velocity of the system is always lower than or equal to 0. The velocity in the two regions does not change its sign; it is important for the practical application to do some work, such as medical robot for intestinal therapy. Furthermore, when the velocity of system changes its sign, the efficiency will decrease because of more energy dissipated by opposite slip. The limit cycle of system in E3 can be described by the signature {I1,I3,I2,I4,I5} depicted in Figure 4(c). The signature of periodic solution in E4 is {I1,I3,I4,I5} showed in Figure 4(d); the signature of the periodic solution in E5 is {I4,I3,I1,I7} depicted in Figure 4(e). The cyclic signature in Figure 4(f) is {I4,I3,I7}.Figure 4
The velocity of the vibration-driven system forα=0.456, (a) ε=0.8,ϕ=π in E1. (b) ε=0.8,ϕ=0 in E2. (c) ε=0.8,ϕ=2π/3 in E3. (d) ε=0.6,ϕ=2π/3 in E4. (e) ε=0.6,ϕ=5π/3 in E5. (f) ε=0.2,ϕ=2π/3 in E6.
(a) (b) (c) (d) (e) (f)The result of two-parameter continuation for the sliding bifurcation with respect to the parametersϕ and ε by fixing different value α ((a) α=0.2; (b) α=0.45; (c) α=0.7) is shown in Figure 5. The regions E3 and E6 are shrinking and the regions E1, E2, E4, and E5 are expanding as the α increases. Therefore it is easy to realize directional motion in E1 or E2 through changing the parameters when the parameter α increases.Figure 5
Two-parameter continuations of sliding bifurcation with respect toϕ and ε for different parameter α: (a) α=0.2; (b) α=0.45; (c) α=0.7; adi represent adding-sliding bifurcations, csi represent crossing-sliding bifurcations, and the regions Ei denote the different stick-slip motions of the periodic solutions.
(a) (b) (c)We carry out the numerical continuation with respect to the parametersε,α in [0,1]×[0,1] for the different ϕ and the results are depicted in Figure 6 ((a) ϕ=π/3; (b) ϕ=2π/3; (c) ϕ=4π/3). The periodic solutions in region Ei in Figures 3 and 6 have the same stick-slip motion; the curve ss1 represents the switching-sliding bifurcation branch. As shown from Figure 6, there are different regions of stick-slip motion of system in the parameter plane when ϕ is different; the motion in E2 which is always equal to or lower than 0 can be realized by changing the parameter ε and α when ϕ=π/3, but the motion in E1 could not happen no matter the value of ε and α. The velocity of system which is always equal to or larger than zero can be controlled through changing the parameter when ϕ is 2π/3 or 4π/3. Therefore the value of ϕ is important for realizing directional motion.Figure 6
Two-parameter continuation of sliding bifurcation with respect toε and α for different ϕ: (a) ϕ=π/3; (b) ϕ=2π/3; (c) ϕ=4π/3; adi represent adding-sliding bifurcations, csi represent crossing-sliding bifurcations, ssi represent switching-sliding bifurcations, and the regions Ei denote the different stick-slip motions of the periodic solutions.
(a) (b) (c)The numerical continuation for the parametersϕ, α in the [0,2π]×[0,1] by fixing the value of the parameter ε is carried out and the results are depicted in Figure 7. From Figure 7, there are no regions E1 and E2 in the parameter plane when ε=0.4, but the regions E1 and E2 appear in the parameter plane as ε increases to 0.6; the regions E1 and E2 expand when ε=0.8.Figure 7
Two-parameter continuation of sliding bifurcation with respect toϕ and α for different ε: (a) ε=0.4; (b) ε=0.6; (c) ε=0.8; adi represent adding-sliding bifurcations, csi represent crossing-sliding bifurcations, ssi represent switching-sliding bifurcations, and the regions Ei denote the different stick-slip motions of the periodic solutions.
(a) (b) (c)


### 4.2. Directional Control

Based on the above analysis, we can see that the direction of the system progression in one period can be forward (the regionE1 in Figure 3) or backward (the region E2 in Figure 3) owing to different parameters.The average velocity of the system in one period is(19)v=1T∫0Ty1tdt.The regions of different directional progression in one period can be determined by implementing the parameters continuation of the periodic solution with additional condition v=0.Whenα=0.456, there are six different stick-slip period solutions in the parameter plane in Figure 3; two parameters ϕ,ε can be continued numerically based on different period solutions with the additional condition v=0. The result is presented in Figure 8. The curves divide the parameters into three regions: forward drift (v>0), backward drift (v<0), and zero drift (v=0, on the curves). They are shown in Figures 9(a), 9(b), and 9(c), respectively. According to Figure 8, it can be seen that, for any value of the parameter ε, the direction of system progression can be controlled by changing the phase ϕ, the direction of progression is forward when ϕ is in (π/2,π), and the direction of progression is backward when ϕ is in (3π/2,2π). When 0.7239<ε<1, the direction of progression is forward when ϕ is in (π/2,3π/2); and the direction of progression is backward when ϕ is in (0,π/2) and (3π/2,2π).Figure 8
Two-parameter continuations of the periodic orbits for the parameterα=0.456. The average velocity is zero along the curves in the plane.Figure 9
The progression of the vibration-driven system for (a)ε=0.5,α=0.456,ϕ=π/2; (b) ε=0.5,α=0.456,ϕ=3π/2; (c) ε=0.4271,α=0.456,ϕ=1.1121.
(a) (b) (c)Whenϕ=2π/3, there is no solution with additional condition v=0 in the parameter ε-α plane, which indicates that the direction of the progression does not change. The direction of the progression is forward because the average velocity is larger than zero in the region E1 seen from Figure 6(b).Whenε=0.6, the parameter ϕ-α plane for directional continuation is presented in Figure 10, some conclusions are draw from Figure 10 similarly: for any parameter α, the direction of system progression can be controlled by changing the phase ϕ; the direction of progression is forward when ϕ is in (π/2,4.4211); and the direction of progression is backward when ϕ is in (3π/2,2π) and (0,1.2795). When 0<α<0.2278, the direction of progression is forward when ϕ is in (π/2,3π/2); and the direction of progression is backward when ϕ is in (0,π/2) and (3π/2,2π).Figure 10
Two-parameter continuation of the periodic orbits for the parameterε=0.6. The average velocity is zero along the curves in the plane.


## 4.1. Sliding Bifurcation Analysis

Whenε=0.6,α=0.456,ϕ=2π/3, the periodic trajectory of the motion is shown by Figure 2. The cycle signature is {I1,I3,I4,I5}. We start the numerical continuation of the periodic solution by the method of path-following using the parameter value as an initial value. The additional boundary hb(y)=sin⁡y2+ε(1+αsin⁡(y2+ϕ))=0 may be applied to the start point of the third segment for the crossing-sliding bifurcation continuation. We will use the COCO to carry out the numerical continuation concerning parameters ϕ and ε; the curve cs1 which is the result of this numerical continuation is shown in Figure 3.Figure 2
Periodic solution of vibration-driven system (11) computed for the parameter values ε=0.6,α=0.456,ϕ=2π/3. The trajectory consists of the segments I1,I3,I4,I5.Figure 3
Two-parameter continuation of sliding bifurcation with respect toϕ and ε for α=0.456. adi represent adding-sliding bifurcations, csi represent crossing-sliding bifurcations, and the regions Ei denote the different stick-slip motions of the periodic solutions.Similarly, the numerical continuation is performed in parameter spaceϕ×ε with different cycle signatures and different additional boundary condition for different segment boundary point and the results are depicted in Figure 3. The signature of cycle trajectory in E1 is {I1,I3,I5}, when the additional conditions hb(y)=sin⁡y2+ε(1+αsin⁡(y2+ϕ))=0 and cos⁡y2+εαcos⁡(y2+ϕ)=0 are applied to the point of the segment I5 according to nondegeneracy conditions of adding-sliding bifurcation; hence we can get the adding-sliding bifurcation curve ad1. The adding-sliding bifurcation branch ad2 is got similarly. The curves cs1, cs2, cs3, and cs4 represent crossing-sliding bifurcation branches. The sliding bifurcation curves divide the two parameters plane into eight regions and there are six different stick-slip periodic solutions in the parameter plane. The periodic solution of system in E1 with signature {I1,I3,I5} is depicted in Figure 4(a), which means the velocity of the system is always greater than or equal to 0. The signature of the periodic trajectory in E2 is {I4,I2,I7} as depicted in Figure 4(b), which means the velocity of the system is always lower than or equal to 0. The velocity in the two regions does not change its sign; it is important for the practical application to do some work, such as medical robot for intestinal therapy. Furthermore, when the velocity of system changes its sign, the efficiency will decrease because of more energy dissipated by opposite slip. The limit cycle of system in E3 can be described by the signature {I1,I3,I2,I4,I5} depicted in Figure 4(c). The signature of periodic solution in E4 is {I1,I3,I4,I5} showed in Figure 4(d); the signature of the periodic solution in E5 is {I4,I3,I1,I7} depicted in Figure 4(e). The cyclic signature in Figure 4(f) is {I4,I3,I7}.Figure 4
The velocity of the vibration-driven system forα=0.456, (a) ε=0.8,ϕ=π in E1. (b) ε=0.8,ϕ=0 in E2. (c) ε=0.8,ϕ=2π/3 in E3. (d) ε=0.6,ϕ=2π/3 in E4. (e) ε=0.6,ϕ=5π/3 in E5. (f) ε=0.2,ϕ=2π/3 in E6.
(a) (b) (c) (d) (e) (f)The result of two-parameter continuation for the sliding bifurcation with respect to the parametersϕ and ε by fixing different value α ((a) α=0.2; (b) α=0.45; (c) α=0.7) is shown in Figure 5. The regions E3 and E6 are shrinking and the regions E1, E2, E4, and E5 are expanding as the α increases. Therefore it is easy to realize directional motion in E1 or E2 through changing the parameters when the parameter α increases.Figure 5
Two-parameter continuations of sliding bifurcation with respect toϕ and ε for different parameter α: (a) α=0.2; (b) α=0.45; (c) α=0.7; adi represent adding-sliding bifurcations, csi represent crossing-sliding bifurcations, and the regions Ei denote the different stick-slip motions of the periodic solutions.
(a) (b) (c)We carry out the numerical continuation with respect to the parametersε,α in [0,1]×[0,1] for the different ϕ and the results are depicted in Figure 6 ((a) ϕ=π/3; (b) ϕ=2π/3; (c) ϕ=4π/3). The periodic solutions in region Ei in Figures 3 and 6 have the same stick-slip motion; the curve ss1 represents the switching-sliding bifurcation branch. As shown from Figure 6, there are different regions of stick-slip motion of system in the parameter plane when ϕ is different; the motion in E2 which is always equal to or lower than 0 can be realized by changing the parameter ε and α when ϕ=π/3, but the motion in E1 could not happen no matter the value of ε and α. The velocity of system which is always equal to or larger than zero can be controlled through changing the parameter when ϕ is 2π/3 or 4π/3. Therefore the value of ϕ is important for realizing directional motion.Figure 6
Two-parameter continuation of sliding bifurcation with respect toε and α for different ϕ: (a) ϕ=π/3; (b) ϕ=2π/3; (c) ϕ=4π/3; adi represent adding-sliding bifurcations, csi represent crossing-sliding bifurcations, ssi represent switching-sliding bifurcations, and the regions Ei denote the different stick-slip motions of the periodic solutions.
(a) (b) (c)The numerical continuation for the parametersϕ, α in the [0,2π]×[0,1] by fixing the value of the parameter ε is carried out and the results are depicted in Figure 7. From Figure 7, there are no regions E1 and E2 in the parameter plane when ε=0.4, but the regions E1 and E2 appear in the parameter plane as ε increases to 0.6; the regions E1 and E2 expand when ε=0.8.Figure 7
Two-parameter continuation of sliding bifurcation with respect toϕ and α for different ε: (a) ε=0.4; (b) ε=0.6; (c) ε=0.8; adi represent adding-sliding bifurcations, csi represent crossing-sliding bifurcations, ssi represent switching-sliding bifurcations, and the regions Ei denote the different stick-slip motions of the periodic solutions.
(a) (b) (c)


## 4.2. Directional Control

Based on the above analysis, we can see that the direction of the system progression in one period can be forward (the regionE1 in Figure 3) or backward (the region E2 in Figure 3) owing to different parameters.The average velocity of the system in one period is(19)v=1T∫0Ty1tdt.The regions of different directional progression in one period can be determined by implementing the parameters continuation of the periodic solution with additional condition v=0.Whenα=0.456, there are six different stick-slip period solutions in the parameter plane in Figure 3; two parameters ϕ,ε can be continued numerically based on different period solutions with the additional condition v=0. The result is presented in Figure 8. The curves divide the parameters into three regions: forward drift (v>0), backward drift (v<0), and zero drift (v=0, on the curves). They are shown in Figures 9(a), 9(b), and 9(c), respectively. According to Figure 8, it can be seen that, for any value of the parameter ε, the direction of system progression can be controlled by changing the phase ϕ, the direction of progression is forward when ϕ is in (π/2,π), and the direction of progression is backward when ϕ is in (3π/2,2π). When 0.7239<ε<1, the direction of progression is forward when ϕ is in (π/2,3π/2); and the direction of progression is backward when ϕ is in (0,π/2) and (3π/2,2π).Figure 8
Two-parameter continuations of the periodic orbits for the parameterα=0.456. The average velocity is zero along the curves in the plane.Figure 9
The progression of the vibration-driven system for (a)ε=0.5,α=0.456,ϕ=π/2; (b) ε=0.5,α=0.456,ϕ=3π/2; (c) ε=0.4271,α=0.456,ϕ=1.1121.
(a) (b) (c)Whenϕ=2π/3, there is no solution with additional condition v=0 in the parameter ε-α plane, which indicates that the direction of the progression does not change. The direction of the progression is forward because the average velocity is larger than zero in the region E1 seen from Figure 6(b).Whenε=0.6, the parameter ϕ-α plane for directional continuation is presented in Figure 10, some conclusions are draw from Figure 10 similarly: for any parameter α, the direction of system progression can be controlled by changing the phase ϕ; the direction of progression is forward when ϕ is in (π/2,4.4211); and the direction of progression is backward when ϕ is in (3π/2,2π) and (0,1.2795). When 0<α<0.2278, the direction of progression is forward when ϕ is in (π/2,3π/2); and the direction of progression is backward when ϕ is in (0,π/2) and (3π/2,2π).Figure 10
Two-parameter continuation of the periodic orbits for the parameterε=0.6. The average velocity is zero along the curves in the plane.


## 5. Conclusions

This paper studies the dynamical response of the vibration-driven system which is composed of a body with movable internal mass. The asymmetry of friction in forward and backward direction is essential to the motion of system, which arises from the normal force change due to the vertical motion of the internal mass. The vibration-driven system involving dry friction belongs to the Filippov system, the cycle trajectory of the system can be divided into smooth segments, the event functions defined the terminal point of the segments, and the reset functions connected the segments. We take advantage of the software COCO to carry out the bifurcation analysis.Two-parameter sliding bifurcations are carried out by performing the numerical continuation. Different period solutions of stick-slip motion are obtained though the sliding bifurcation curves. For directional control of the vibration-driven system, the drift of the vibration-driven system in one period may change sign, the continuation with additional conditionv=0 is carried out in the parameters plane, and the curves are obtained. The curves divide the parameters plane into three modes of drift (backward, forward, and zero). So the direction of the vibration-driven system progression can be controlled by changing the parameters.The particular contribution of this research is the numerical continuation of the parameters for the vibration-driven system and detailed classifications of parameter space where different system dynamic behaviors can be obtained. The bifurcation analysis improves our understanding of the dynamical behaviors of the vibration-driven system and is of benefit to devise control strategies for the vibration-driven system.



---

*Source: 2305187-2018-08-27.xml*