To the modeling of processes with repeated damped collisions

The problem of calculating and analyzing vibro-impact processes with bounce in two colliding subsystems is considered. The nature of these processes depends on the adopted models of energy dissipation and the design features of the colliding objects. A method for constructing the laws of motion is presented, based on the technique of averaging in vibro-impact systems. Examples are given. The determining role of energy dissipation models is shown. The results can be used, for example, in the calculation and analysis of specific phenomena that arise in the elements of switching devices, as well as in many other systems whose operation is accompanied by repeated damped collisions.


Introduction
In the general case, regimes with repeated collisions (modes with bounce or chatter) can be called movements accompanied by many collisions of contacting solids after the first impact occurred.In the simplest case, this phenomenon consists in the known effect of the bobbing of a coin, which at a certain height was released above a massive metal or stone slab.From a mechanical point of view, a bounce occurs when the colliding bodies are in a field of forces that return them to the point of contact, resulting in repeated collisions.
It is known that in several cases, the bounce of contacts can cause great damage to technical objects.It can occur and occurs in numerous electromechanical switching devices, as well as in the components of various systems -buttons, relays, reed switches, contactors, magnetic starters and so on.The effect is manifested in the fact that after the first connection numerous repeated disconnects and closures of the contacts are observed.
This significantly reduces the reliability of devices, introduces undesirable distortions into their work and can lead to premature failures and loss of operability of equipment elements.For example, the commutation of powerful electric circuits can be accompanied by ignition and extinction of the electric arc and the appearance of sparks, which causes rapid contact wear.
The most important characteristics of the bounce are impulses of impacts, the number of collisions from the beginning of the first contact to the last, and the duration (time) of the bounce.These parameters depend on the design parameters of the systems.Bounce times range from order values of 1 to 2 microseconds from modern miniature reed switches to hundreds of milliseconds in powerful contactors.
The methods used to combat bounce are based on hardware and algorithmic techniques.Their description is widely presented on the Internet, due to the very wide distribution of devices related to the Arduino platform [1].Basically, these methods are palliative, do not consider the mechanics of processes and cannot be completely effective.Mathematical modeling of electromechanical contacts or similar systems, apparently, was first conducted by L. Mandelstam [2].It should be noted the work of R. Nagaev [3], who gave complete and detailed descriptions of regimes with repeated collisions for several nontrivial mechanical problems.
In this paper, we consider some model problems and methods for describing the chatter by means of the averaging technique.

Modeling approach
A typical example of contact systems of this type is shown in Fig. 1.As the main working model, we take the system shown in Fig. 2(a).At the ends of elastic weightless beams (springs), colliding point bodies with masses  and  [4] are fixed.The assumption that bodies are pointlike means that all their geometric parameters can be neglected in comparison with the gap size Δ ≥ 0 (negative values Δ correspond to the preload).
where  and  are the coordinates of the colliding bodies;  =  −  .Collisions occur at  = Δ .We denote by  (;  ) and  (;  ) operator functions describing the energy dissipation in systems without considering the energy loss during impacts ( ≡ /).Here, not only the standard factors leading to the use of viscous friction models can be taken into account, but also the circumstances associated with internal imperfections and peculiarities of the material of the beams (springs).For example, they can be made from materials with relaxation.In Eq. ( 1) the function Φ(,  ) denotes the impact force, which depends on the relative coordinate and velocity.Let the impact be modeled according to Newton.Then at the instant of impact  =  we have  ( − 0) = − ( + 0); 0 <  ≤ 1 and Φ(,  ) = Φ (,  ) + (1 − )Φ (,  ).Here, the elastic [Φ ] and dissipative [Φ ] components of the impact force, written through a generalized function.For each impact [5]: Φ| = ( −  );  = (1 + )| ( − 0)|, where  is impact impulse; () is Dirac -function;  is the moment of some impact; 0 <  ≤ 1 is restitution ratio.

Method of analysis
We analyze Eq. ( 2) using the averaging technique in form, which in respect of vibro-impact systems was developed in books [5].The solution of the problem for the conservative case  = 0 is well known and studied.It is described by expressions: Here Ω = √ is natural frequency of a linear oscillator;  = 2 is the process period.The second formula in Eq. ( 3) is valid only for 0 <  ≤ .If the variable  lies outside this interval, then it is necessary to use the conditions of periodicity or the representation of the solution in the form of a Fourier series [5].The impact impulse  is nonnegative.In a conservative system, it is determined from the impact condition (0) = Δ.In an autonomous system, we can assume that the moment of impact coincides with the beginning of the time count.

Results and discussion
Let's break down several important examples.First, let  = 0 and on the process only affect the energy loss in collisions (  < 1).Using the third relation Eq. ( 3), we find for Δ > 0, (Ω <  < 2Ω)  = − Ω  − arctg .The only stationary regime corresponds to the termination of collisions (the impact pulse is zero):  =  = 0.This regime is asymptotically stable, since for a small perturbation of the stationary regime ( ) the equation in variations and its solution take the form:  = − Ω,  =  exp − Ω .So, the time of returning to the undamped free sinusoidal oscillations (the time of termination of the vibro-impact process, with a bounce) is infinite in such a model.
The general analysis shows that for Δ < 0,  ≠ 0,  < 1, the bounce time is estimated with the help of  , and for Δ > 0,  ≠ 0,  < 1 we use the estimate for the quantity  .For systems with preload, we have a process with an infinite number of collisions (quasiplastic impact).
It can easily be shown that for Δ = 0 both mechanisms of dissipation have the same effect on the process.The time for the termination of the regime with collisions is infinite.
Let's consider one more example.Let the energy losses occur due to the properties of the material.For example, the beams (springs) (Fig. 2) have a polymeric braid.Then: where Γ() is the relaxation kernel [7], whose parameters are determined by the characteristics of the material of a beam with a polymer braid.The right side of the Eq. ( 6) ( = 0) is: We find that  () = −4 Γ () (Ω −  ) ( [5]), where Γ is the sine of the Fourier transformation of the relaxation kernel.The dissipative properties of the material are determined precisely by this quantity [7].Further analysis depends on the specific type of relaxation nucleus.A detailed examination here is difficult.However, qualitative differences with the cases considered are small.The influence of material properties on the bounce should be considered only when carrying out high-precision calculations.
Thus, using the introduced change of variables, one can obtain a representation for the laws of variation of the relative and absolute coordinates: where the function is determined with the help of the second relation Eq. ( 3) or in other.Fig. 3 shows the laws of motion of colliding bodies and the times of rebound.For the case Δ = 0, the decay time of the vibro-impact process is infinite, and the model needs to be refined.Fig. 3,  shows the law of change in the relative coordinate.A similar graph will correspond to a body colliding with a fixed massive obstacle.

Conclusions
The problems of analyzing the motions of systems with repeated damped collisions of elements are of engineering value, since they make it possible to estimate the parameters of the bounces accompanying these movements.
It is shown that the averaging technique in vibro-impact systems makes it possible to obtain important relationships for estimating the motion parameters of such systems.This is important for carrying out practical calculations.
It is shown that the character of the bounce modes in a special way depends on the models of energy dissipation and structural parameters.

Fig. 3 .
The laws of motion of colliding bodies