# External Periodic Force Control of a Single-Degree-of-Freedom Vibroimpact System

Authors: Jingyue Wang; Haotian Wang; Tie Wang
Journal: Journal of Control Science and Engineering
 (2013)

---

## Abstract

A single-degree-of-freedom mechanical model of vibro-impact system is established. Bifurcation and chaos in the system are revealed with the time history diagram, phase trajectory map, and Poincaré map. According to the bifurcation and chaos of the actual vibro-impact system, the paper puts forward external periodic force control strategy. The method of controlling chaos by external periodic force feedback controller is developed to guide chaotic motions towards regular motions. The stability of the control system is also analyzed especially by theory. By selecting appropriate feedback coefficients, the unstable periodic orbits of the original chaotic orbit can be stabilized to the stable periodic orbits. The effectiveness of this control method is verified by numerical simulation.


---

## Body

## 1. Introduction

In the field of nonlinear, along with people understanding the nature of chaos, how to control chaos and chaos synchronization has been a hot topic studied by researchers. In the 1990s, Ott et al. proposed the OGY chaos control method [1]. Scholars put forward a lot of modified control methods [2, 3] based on the OGY method. Although small perturbation of the system parameters for the control of chaos has been confirmed, the method requires detailed information about the target trajectory and brings a lot of inconvenience to the practical application. Therefore the domestic and foreign researchers put the traditional control theory and chaotic motion characteristics used in chaos control and present a lot of chaos control method, such as delayed feedback control [4, 5], periodic parameter perturbation control [6], continuous feedback control [7], pulse feedback control [8], and adaptive control [9], and so forth. There are also examples of control on the vibroimpact system [10, 11].Vibroimpact system as a typical nonsmooth dynamical system generally exists in practical engineering. Because of the frequent collision, the system has strong nonlinearity and discontinuities compared with a smooth nonlinear system, presents more complex nonlinear phenomena, and causes hazards on the safe operation of the system. Because of the collision interface differential discontinuities, the original method applied to continuous system can not be used for such system.This paper puts forward a sine periodic force feedback controller based on the periodic external force feedback control strategy and analyzes the stability of control theory. When selecting the appropriate feedback coefficients, the chaotic orbits can be controlled onto the stable periodic orbits. A single-degree-of-freedom vibroimpact system is transformed into a form of Poincaré map for numerical simulation. The results of numerical simulation show that the method is effective in practical engineering, so it has certain practical significance.


## 2. Mechanics Model of the Vibroimpact System

Figure1 shows a single-degree-of-freedom mechanical model of vibroimpact system. Oscillator M is connected to the left side support by the spring with stiffness K and the damper with damping C. In the harmonic excitation Fsin(ΩT+θ), motion of the oscillator is in the horizontal direction. X is the displacement of motion. When the oscillator is in the equilibrium position, the gap between the rigid constraints on the right side is D. Considering the collision as the rigid collision and R as the coefficient of restitution, the differential equation of motion of system is
(1)MX¨+CX˙+KX=Fsin(ΩT+θ)(X<B).Figure 1
Mechanics model of a single-degree-of-freedom vibroimpact system.The shock equation of system is(2)X˙+=-RX˙-(X=B),
where X˙- and X˙+ represent the impacting mass velocities of approach and departure at the instant of impacting, respectively.After the dimensionless transformation, whenx<d, the differential equations of motion of the system between two collisions are as follows:
(3)x¨+2ζx˙+x=sin(ωt+θ).
In which, the nondimensional quantities are
(4)x=XKF,ζ=C2MK,ω=ΩMK,t=TKM,d=DKF.Whenx=d, shock equation of the system at collision transient is given by
(5)x˙+=-Rx˙-,
where x˙- and x˙+ represent the impacting mass velocities of approach and departure at the instant of impacting, respectively.By (1), the general solution of the system between two collisions between is
(6)x=e-ζt(acosωdt+bsinωdt)+Asin(ωt+θ)+Bcos(ωt+θ),
where ωd=1-ζ2. A and B are amplitude constants. By the initial conditions of the system: x(t0)=x0 and x˙(t0)=x˙0, let t0=0; the integral constants A and B can determined
(7)a=x0-Asinθ-Bcosθ,b=(x˙0+ζx0+(Bω-Aζ)sinθ-(Aω+Bζ)cosθ)ωd,A=1-ω2(1-ω2)2+(2ζω)2,B=-2ζω(1-ω2)2+(2ζω)2.Periodic motion of the system under certain parameter conditions can be expressed asn-p, n represents a force cycle number, and p represents the number of collisions. Considering the n-1 periodic motion, the collision instantaneous dimensionless time t=0, and next collision time tz=2nπ/ω, then boundary conditions of the n-1 periodic motion are x(0)=x(tz)=d and x˙(0)=-Rx(tz). With (6) applying to the boundary conditions, the existence conditions of periodic motion system are
(8)|-Htanθ--1+tan2θ--H21+tan2θ-|≤1,
where θ-=arctan(-H2/H1), c=cosωdtz, e=e-ζtz, s=sinωdtz,
(9)H=dH1((ζ+Re(ωds+ζc))es+(ωd+Re(ωdc-ζs))(ec-1)),H1=A(ζ+Re(ωds+ζc))-Bω(1+R)es+A(ωd+Re(ωdc-ζs))(ec-1),H2=B(ζ+Re(ωds+ζc))-Aω(1+R)es+B(ωd+Re(ωdc-ζs))(ec-1).At the same time, theory fixed point of then-1 periodic motion is given by
(10)θd=arccos(-Htanθ--1+tan2θ--H21+tan2θ-),xd=d,x˙d=(d-Asinθd-Bcosθd)((1-ec)ωd-esζ)es+Aωcosθd-Bωsinθd.Define the following section:σ={(θ,x,x˙)∈R2×S1|x=d}. In the paper, we choose the section σ to establish the Poincaré map
(11)θk+1=ωt+θk,(12)xk+1=e-ζt(a^cosωdt+b^sinωdt)+Asin(ωt+θk)+Bcos(ωt+θk),(13)x˙k+1=-R(e-ζt((b^ωd-a^ζ)cosωdt-(b^ζ+a^ωd)sinωdt)+Aωcos(ωt+θk)-Bωsin(ωt+θk)b^ωd)),
where a^=xk-Asinθk-Bcosθk,   and  b^=(x˙k+ζxk+(Bω-Aζ)sinθk-(Aω+Bζ)cosθk)/ωd. In which, initial iteration value of the Poincaré map of the n-1 periodic motion is
(14)θ0=θd,xk+1=xk=x0=xd=d,x˙0=x˙d.


## 3. Chaos and Bifurcation

The single-degree-of-freedom mechanical model of vibroimpact system, with system parametersζ=0.01, d=0.05, and R=0.8, has been chosen to be analyzed. The system parameter ω is taken as the bifurcation parameter. The global bifurcation diagram can be obtained with the ω changing in the range of [0,6] as shown in Figure 2. As you can see from Figure 2, the system has stable n-1 periodic motion in a certain range. But in different period of single touch movement, the system will produce double periodic bifurcation, from periodic motion to chaos in the process of transition. But with the increase in ω, the system will degenerate to the periodic motion. The Figures 2(b) and 2(c) are a partial enlargement of Figure 2(a). It can be seen that, when ω=2.5585, the vibroimpact system will be from periodic 1 motion to periodic 2 motion by bifurcation and when ω=2.624, the vibroimpact system will be from periodic 2 motion to periodic 4 motion by double periodic bifurcation. When ω=2.641, the vibroimpact system will be from periodic 4 motion to periodic 8 motion by double periodic bifurcation. When ω=2.644, the vibroimpact system will be from periodic 8 motion to periodic 16 motion by double periodic bifurcation. With the vibration frequency ω increasing further, the system leads to chaotic motion.Bifurcation diagram of the system.
(a)
(b)
(c)The excitation frequencyω takes 2.55, 2.6, 2.635, 2.642, 2.644, and 2.652. Poincaré map, phase portrait, and time course diagram of the system are shown in Figures 3, 4, 5, 6, 7, and 8. The system is, respectively, periodic 1 motion, periodic 2 motion, periodic 4 motion, periodic 8 motion, periodic 16 motion, and chaotic motion.Poincaré map, phase portrait, and time course diagram of the system withω=2.55.
(a)
(b)
(c)Poincaré map, phase portrait, and time course diagram of the system withω=2.6.
(a)
(b)
(c)Poincaré map, phase portrait, and time course diagram of the system withω=2.635.
(a)
(b)
(c)Poincaré map, phase portrait, and time course diagram of the system withω=2.42.
(a)
(b)
(c)Poincaré map, phase portrait, and time course diagram of the system withω=2.644.
(a)
(b)
(c)Poincaré map, phase portrait, and time course diagram of the system withω=2.52.
(a)
(b)
(c)


## 4. Chaos Control

The paper chooses the sine driving force for the periodic force excitation. Periodic force is easy to produce and control the external driving force in the actual project, so the sine driving force is used to suppress the bifurcation and chaotic motion of single-degree-of-freedom vibroimpact system. Based on the principle of parameter perturbation, periodic force excitation method can stabilize the chaotic motion by applying disturbance directly into the system. An unstable periodic motion of the chaotic system can produce resonance with external periodic force. The system can be from its unstable limit cycle to a stable limit cycle by resonating with external driving signal. So the chaos is controlled.Periodic force excitation can be expressed as(15)f=hsin(γωθ).
In the formula, h and γ are adjustable parameters.As the previous analysis, whenω=2.652, the system is in chaotic motion. After the introduction of the external periodic force feedback controller, using the fourth-order Runge-Kutta method of numerical simulation, the Poincaré map, phase portrait, and time course diagram of the controlled system are obtained, as shown in Figures 9 and 10. When γ=0.05 and h=4.6, the system is controlled to periodic 1 motion. When γ=0.05 and h=5.6, the system is controlled to periodic 2 motion. The simulation results show the effectiveness of the method. Because the method does not change the original system parameters, it can be applied to a 2-DOF and multiple-DOF vibroimpact mechanical system.Poincaré map, phase portrait, and time course diagram of the controlled system withγ=0.05 and h=4.6.
(a)
(b)
(c)Poincaré map, phase portrait, and time course diagram of the controlled system withγ=0.05 and h=5.6.
(a)
(b)
(c)


## 5. Conclusion

Based on a single-degree-of-freedom vibroimpact system as the research object, bifurcation and chaos have been researched with the system parameters changing. By adopting an external periodic force excitation method for suppressing its chaotic behavior, it delayed the occurrence of fault. Because this method does not change the original system parameters, it is easy to implement in engineering. This method is not limited to this kind of mechanical system with clearance collision and can also be used in other similar nonlinear system.



---

*Source: 570137-2013-11-21.xml*