Control of Chaos in Rate-Dependent Friction-Induced Vibration Using Adaptive Sliding Mode Control and Impulse Damper

Two different controlmethods, namely, adaptive slidingmode control and impulse damper, are used to control the chaotic vibration of a block on a belt system due to the rate-dependent friction. In the first method, using the sliding mode control technique and based on the Lyapunov stability theory, a sliding surface is determined, and an adaptive control law is established which stabilizes the chaotic response of the system. In the second controlmethod, the vibration of this system is controlled by an impulse damper. In this method, an impulsive force is applied to the system by expanding and contracting the PZT stack according to efficient control law. Numerical simulations demonstrate the effectiveness of both methods in controlling the chaotic vibration of the system. It is shown that the settling time of the controlled system using impulse damper is less than that one controlled by adaptive slidingmode control; however, it needs more control effort.


Introduction
There exist a lot of works on the theory of friction-driven oscillations in the literature, for example, influence of the belt speed on the system response [1], dynamics of three-block mechanical system with dry friction [2], investigation on the geometry of chaotic attractors for dry friction oscillators [3,4], influence of parametric and external excitations on a dry friction oscillator dynamics [5], and the dynamic behavior of friction-driven oscillator with an impact damper [6].
The characteristics of the friction force between two surfaces are quite complex and depend on many parameters such as, normal pressure, slip velocity, surface, and material properties [7].The friction-actuated oscillation is strongly nonlinear, and discontinuous and has nonsmooth process, which is a source of instabilities generating stick-slip, chatter, squeal, and chaos [8,9].LuGre friction law is one of the most widely used friction laws which models rate dependency of the friction force by one additional inner variable [10].In [11], LuGre friction model is applied to the single degree of freedom friction-induced oscillator, and it is shown that the oscillations of this system turned out to be chaotic for most parameter combinations.
Chatterjee presented a novel method for controlling the vibration, called impulse damper [25].In this method, mass loaded PZT stack is attached to the system, and applying suitable voltage to the PZT stack leads to impulsive force which can control the system by setting control parameters correctly.This technique was used to control a nonlinear friction-driven oscillator [26] and to control the chaotic vibration of friction-driven oscillator with coulomb friction model [27,28].
In the present study, two different control methods, namely, impulse damper and adaptive sliding mode control, are used to control the chaotic vibration of a single degree of freedom system subjected to the rate-dependent friction.To control the chaotic vibration of this system using adaptive sliding mode control, an adaptive control law is established which stabilizes the chaotic response of the system.
Since chaotic systems are sensitive to the initial conditions, applying appropriate impulses at suitable positions causes the system to maintain its path within a specified bound.This was our motivation to investigate the effectiveness of impulse damper to suppress the chaotic behavior of the system.The effectiveness of both methods will be investigated in following sections.

Mathematical Model
In this section, the mathematical modeling of uncontrolled and controlled system, using adaptive sliding mode control and impulse damper is presented.

Friction Mathematical Model.
Experimental investigations, especially at low sliding speeds, have shown that friction laws which give the friction force as an algebraic function of the underlying parameters (i.e., static friction laws) do not capture all of the observable frictional effects.Among the additional effects is the presliding displacement due to lateral contact elasticity, the increase of static friction with time due to diffusion processes on the interface, and frictional lag in sliding, which stands for the effect of the friction force lagging in time behind changes in relative velocity or normal load.To model the effects, the so-called dynamic friction models have been proposed.The following investigation is based on one of the most widespread models, the LuGre model [10], which allows a comparatively simple representation of the dynamic effects.Based on this model, the friction force is given by Governing equation of  is where  represents the average deflection of the bristles, V  is the relative sliding velocity,  is the normal load with a reference value  0 , and (V  ) is the friction force which is formulated as a function of the constant relative sliding velocity V  .Moreover,  0 ,  1 , and  2 parameterize the presliding displacement, an internal viscous frictional damping, and a viscous damping contribution due to the relative velocity [11], respectively.For steady sliding,  2 vanishes and the friction force reduces to Function () is given as follows: where V  is the speed of belt.In this research, sliding block on a rigid belt moving with a constant velocity as shown in Figure 1 is considered.Using LuGre friction, the equations of motion are where the dot denotes differentiation with respect to .
Introducing the new parameter  =   , where   = √/, the following relations will be concluded.
It should be noted that the dimension of ( 6) and ( 7) is length.In these equations, the prime denotes differentiation with respect to , and the values of parameters are set as follows.
Here, a relative damping measure  = ( +  2 )/2 √  has been introduced to combine the effects of linear damping from both structure and friction.
It is shown that for these values of parameters, the system has chaotic vibration [11].

Adaptive Sliding Mode Controller (ASMC) Design.
In this section, the adaptive sliding mode control is used to control the chaotic vibration.The controlled chaotic system can be rewritten as follows: It is desired that  be adjusted to equilibrium state  0 , where  0 is a given constant.The error states are defined as To guarantee the stability of the sliding mode, a switching surface () in the error space is defined as follows: As well known, when the system operates in the sliding mode, it satisfies the following [29]: Therefore, the following sliding mode dynamics can be obtained as From ( 8), it can be seen that the internal parameter  will converge and be stabilized when ( 1 ,  2 ) → (0, 0).The following function is defined: Using ( 13), ( 12) can be rewritten as The controller is designed as where  = |( 1 ,  2 , ,   )| + | 2 | and the adaptive law for parameter  is proposed as The proposed adaptive control scheme will guarantee the global asymptotic stability of the error and is shown in the following.

Stability Analysis.
In the following, we analyze the stability of the sliding mode dynamics based on the Lyapunov stability theory.The Lyapunov function is selected as () = (1/2)( 2 +  2 ), which leads to Solving ( 16) leads to  =  − + .Substituting ( 16) in ( 17), the following is obtained: Taking the limit as  → ∞ on both sides of (19) yields lim Now, Barbalat's lemma is applied as shown in the following.
Thus, the system is in the sliding mode, and using (12), stability of  1 and  2 is surely guaranteed.Therefore, both  1 () and  2 () converge to zero.

Mathematical Model of Impulse Damper.
In this section, for the purpose of controlling and suppressing vibration of the system, an impulse damper is used (Figure 2).The basic principle of the impulse damper is to generate suitable impulsive force by quickly expanding or contracting the PZT stack actuator.The stack is expanded by applying a suitable positive voltage at the instants of the displacement zero crossing with a positive velocity.In addition, a negative voltage is applied to contract the stack when the zero crossing takes place with a negative velocity.In fact, damping is generated by velocity feedback.
Applied force by impulse damper can be written as   +     , where   is the elongation of PZT stack and   is the force produced by the actuator.  can be found from the electromechanical characteristics of the PZT actuator which is given by the following: where   ,   ,   ,   , and C PZT represent charge, elongation, voltage, force, and the effective capacitance of the PZT actuator, respectively. 33 is the piezoelectric constant of each wafer, and  is the total number of wafers in the actuator.  −1 is the compliance of the actuator in short circuit condition and is given by   = /, where  is the elastic modulus,  is the cross-sectional area, and  is the length of the actuator.From ( 22), the following expression is obtained for the actuator force: Employing the previous equation and adding it to the right side of ( 6) as an actuation force which is applied to the system by the impulse damper, one can conclude that The equation of motion of the added mass can be written as [25] where and  2 =  2  2 33   / PZT .Depending upon the problem, the reference displacement  0 can be suitably defined.The following values are chosen for the previous mentioned parameters in this study:   = 0.1,  = 0.7, ℎ  = 25,   = 4,  = 0.003, and  = 100.
The proposed control law is mathematically recast as [25]

Results and Discussion
Equations of this system are solved numerically.For the numerical analysis, Runge-Kutta's integration procedure is employed where the initial condition of system is set as ((0), ẋ (0), (0)) = (0, 0, 0). Figure 3 shows the time history of the mass block displacement, and Figure 4 shows the phase plane of the uncontrolled system.Chaotic dynamics of the system can be clearly seen in these figures.
In the following, the numerical results are given to confirm the validity of the proposed adaptive sliding mode method.In the numerical simulations,  is selected as  = 1.5, and the desired response is set at  0 = 0.0383 (equilibrium point).Synchronization error of the state variables is depicted in Figure 5.
Phase plane of the controlled system has been shown in Figure 6.As can be seen in this figure, the controlled system is not chaotic.
The time response of the switching function () and the control input () are shown in Figures 7 and 8, respectively.
For the purpose of indicating the efficacy of the proposed impulse damper controller, we have examined the system in two different situations.In the first situation, the controller is turned on when the system is at origin, and accordingly, the mass moves towards its equilibrium position at  = 0.0383.In the second situation, controller is off for  < 100; therefore, the mass is vibrating chaotically around its equilibrium position at  = 0.0383.At this moment, the controller is turned on and the effect of the control input on the behavior of the system is evaluated.After that, to show the usefulness of the controller, the controller is turned off, and there is no suppressing signal to damp system vibrations.Figure 9 shows response of the system controlled by the impulse damper when the controller is turned on at  = 0. Figure 10 shows applied force by the impulse damper when the controller is turned on at  = 0. Figures 11 and 12 illustrate the response and control input for the system controlled by the impulse damper when the controller is turned on at  = 100 and is turned off at  = 200.In Figure 9, when the controller is turned on at  = 0, response of the system is absolutely fast.The mass reaches its equilibrium position in a very short time.Although previously mentioned quick response and short settling time demonstrate high performance of the controller, it is achieved at the expense of a high actuating force (shown in Figure 10) which is applied by the impulse damper.The achieved result is rational since the mass is at origin and, thus, is far from its equilibrium position.Therefore, impulse damper controller applies big impulses to the mass for the purpose of moving it towards the equilibrium position.In contrast, when the controller is off for  < 100 and thus the chaotic trajectory enters a small neighborhood of the fixed point  = 0.0383 (m)   (Figure 10), there is no need for such a high force as in the previous case and impulses with a short amplitude as shown in Figure 12 can control the system effectively.In this case, the distance between the current position of the mass and its equilibrium position is not as far as the first case.In other words, the farther (close) the system is to the goal dynamics, the bigger (small) the weight given to the control effort is.Feasibility of using impulse damper for controlling the system is obvious in these figures.As seen in these figures, impulse damper can damp the vibration of the system absolutely in a short time.

Conclusion
In this paper, chaotic vibration control of a single degree of freedom oscillator subjected to a LuGre type friction law was presented.Two different methods were used to control the chaos, ASMC, and impulse damper.
For adaptive control of the chaotic vibration of this system, a switching surface was adopted such that it becomes easy to ensure the stability of the error dynamics in the sliding mode.Then, an ASMC was derived to guarantee the occurrence of the sliding motion.The adaptive laws were derived in the Lyapunov sense to guarantee the stability of the controlled system.Furthermore, impulse damper was also used for controlling the system.The idea is that, since chaotic systems are exponentially sensitive to perturbations, careful choice of even small control perturbations can, after some time, have a large effect on the trajectory location and can be used to guide it.In this method, controlled impulses are generated by expanding and contracting a mass loaded PZT actuator used between the primary system and the secondary mass.Numerical results verified the effectiveness of the proposed methods in controlling the chaotic vibration of the rate-dependent friction-driven oscillator.It was shown that impulse damper decreases the settling time of the controlled system by applying higher control forces in comparison with ASMC.

Figure 1 :
Figure 1: The structure of the block on a rigid belt.

Figure 2 :
Figure 2: PZT stack attached to the system as controller.

Figure 3 :Figure 4 :
Figure 3: Block displacement of the uncontrolled system.

Figure 5 :Figure 6 :
Figure 5: Time response of error states using adaptive sliding mode method.

Figure 7 :
Figure 7: Time response of the switching function using adaptive sliding mode method.

Figure 8 :Figure 9 :Figure 10 :
Figure 8: Time response of the input control using adaptive sliding mode method.

Figure 11 :
Figure 11: Block displacement (controller is turned on at  = 100 and is turned off at  = 200).

Figure 12 :
Figure 12: Applied force (controller is turned on at  = 100 and is turned off at  = 200).