To the modeling of processes with repeated damped collisions

Vitaly Krupenin1

1Mechanical Engineering Research Institute of the Russian Academy of Sciences, Moscow, Russia

1Corresponding author

Vibroengineering PROCEDIA, Vol. 20, 2018, p. 36-40.
Received 4 October 2018; accepted 7 October 2018; published 19 October 2018

Copyright © 2018 Vitaly Krupenin. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Creative Commons License
Table of Contents Download PDF References
Cite this article
Views 29
Reads 12
Downloads 614
CrossRef Citations 0

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.

Keywords: collision, bounce, bounce time, contact system, variables “impulse-phase”, dissipation models, systems with relaxation.

1. 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.

2. Modeling approach

A typical example of contact systems of this type is shown in Fig. 1.

Fig. 1. The electronic control relay “PMDsrange” of the company “Pilz International”, Germany

 The electronic control relay “PMDsrange” of the company “Pilz International”, Germany

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 m1 and m2 [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).

Fig. 2. The main working model

The main working model


The main working model


The equations of motion of the system in the absence of constantly acting external excitation forces are obviously of the form:

m 1 x ¨ 1 + c 1 x 1 + Φ x , x ˙ = - B 1 p ; x 1 ,       m 2 x ¨ 2 + c 2 x 2 - Φ x , x ˙ = - B 2 p ; x 2 ,

where x1 and x2 are the coordinates of the colliding bodies; x=x1-x2. Collisions occur at x=Δ. We denote by B1(p;x1) and B2(p;x1) operator functions describing the energy dissipation in systems without considering the energy loss during impacts (pd/dt). 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 Φ(x,x˙) 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 t=tα we have x˙(tα-0)=-Rx˙(tα+0); 0 <R 1 and Φ(x,x˙)=Φ0(x,x˙)+(1-R)Φ1(x,x˙). Here, the elastic [Φ0] and dissipative [Φ1] components of the impact force, written through a generalized function. For each impact [5]: Φ|t=tα=Jδ(t-tα);J=(1+R)|x˙(tα-0)|, where J is impact impulse; δ(t) is Dirac δ-function; tα is the moment of some impact; 0 <R 1 is restitution ratio.

In practice, many commutation systems in which bounce occurs have a regular or quasiregular structure, that is, for our model m1m2, c1c2, B1B2. Let the operators B1, B2 be linear, then from Eq. (1) after the calculations we find for the relative coordinate:

m x ¨ + c x + Φ 0 x , x ˙ + ε r Φ 1 x , x ˙ = - ε B p ; x .

Here 12m=m1=m2, 12 c =c1=c2, 12B=B1=B2; 1-Rr; ε is a small parameter (we consider quasi-conservative models). If m1m2, the model is analogous: we can assume that m 1  and the problem reduces to an analysis of the impact oscillator interacting with a fixed limiter (see Fig. 2(b), Δ< 0).

3. 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:

x t = - J χ t ,       χ t = 1 2 Ω cos Ω t - T 2 sin Ω T 2 ,       0 < t T ,       J = - 2 Ω Δ t g π Ω ω - 1 0 .

Here Ω=cm-1 is natural frequency of a linear oscillator; T=2πω-1 is the process period. The second formula in Eq. (3) is valid only for 0 <tT. If the variable t 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 J is nonnegative. In a conservative system, it is determined from the impact condition x(0)=Δ. In an autonomous system, we can assume that the moment of impact coincides with the beginning of the time count.

We make the transition to the new variables “impulse -phase” (x,x˙)(J,ψ) in accordance with formulas x=-Jχ[ψ,ω(J)],x˙=-Jω(J)χψ[ψ,ω(J)]. In this case, the dependence of ω(J) is determined from the fourth formula Eq. (3). In this case, we choose if Δ> 0, Ω<ω<2Ω, if Δ 0, 2Ω<ω<, if Δ= 0, ω=2Ω (see [5]). After the transformations, the following equations can be obtained for the new variables:

J ˙ = - ε ω J r J δ 2 π ψ + 4 B p ; - J χ ψ , ω J χ ψ ψ , ω J ,
ψ ˙ = ω J 1 + 4 ε J - 1 B p ; - J χ ψ , ω J χ ψ , ω J + J χ ω ψ , ω J ω ' J .

Here δ2π(ψ) is a generalized 2π-periodic δ – function. To determine a slowly changing evolutionary component of the process, it is enough to average Eq. (4) over the fast phase. Let the initial conditions look like this: x(0)=Δ, J=J0> 0. Averaging Eq. (4), after a series of calculations, one can obtain:

J ˙ = ε E 1 J - ε r T - 1 J J + 4 0 T J B p ; - J χ ψ , ω J χ ψ ψ , ω J d ψ .

In the first approximation, Eq. (5) is not used [6]. Now considering various dissipation models it is possible to estimate the regime parameters with repeated damped collisions.

4. Results and discussion

Let's break down several important examples. First, let B= 0 and on the process only affect the energy loss in collisions (R< 1). Using the third relation Eq. (3), we find for Δ>0, (Ω<ω<2Ω)J˙=-12εrJΩπ-arctgJ2ΩΔ-1. The only stationary regime corresponds to the termination of collisions (the impact pulse is zero): J=J0= 0. This regime is asymptotically stable, since for a small perturbation of the stationary regime (δJ0) the equation in variations and its solution take the form: δJ˙=-12εrΩδJ, δJ=δJ0exp-12εrΩt. 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.

For a system with preload (Δ< 0, 2Ω<ω<), Eq. (6) takes the form (B= 0).

J ˙ = - 1 2 ε r J Ω a r c t g J 2 Ω Δ - 1 . In this case, there is no stationary solution, since for J~0, J˙> 0. Assuming the momentum to be small, we linearize the last equation and find its solution: J˙=εrΩ2Δ=-εrΩ2|Δ|, Jt=J0-εrΩ2Δt,J0-ε. Now we can estimate the time of the bounce (tb):tb1=J0(εrΩ2|Δ|)-1. In this case J0:ω(J), T(J) 0. Therefore, an infinite number of collisions must occur in time tb1. R. Nagaev called such a dynamic phenomenon a quasiplastic impact [3]. This is a direct result of idealization, which arises from Newton’s hypothesis.

Consider secondly the case R= 1; Δ> 0, B(p;x)=2bpx2bx˙, b> 0. Using Eqs. (3), (5) after making some calculations (see also [5]) we obtain the equation:

J ˙ = - ε J b ( 4 Ω 2 Δ 2 J - 2 + 1 ) { 1 - Ω T J ] - 1 s i n Ω T J ε E 1 J .

We will estimate the time of the bounce. Let the initial impulse be small enough: J0-ε. Then, estimating the function E1(J), we can show that εE1J=-4εΩ2Δ2bJ-1+ε2 and J(t)=J 02-8εΩ2Δ2bt, tb2=J 02(8εΩ2Δ2b)-1. In this case, when Δ> 0, the period tends to a finite value for J 0, so the number of collisions per time is finite.

The general analysis shows that for Δ< 0, b 0, R< 1, the bounce time is estimated with the help of tb1, and for Δ> 0, b 0, R< 1 we use the estimate for the quantity tb2. 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:

- B p ; x = - t Γ t - s x s d s ,

where Γ(t) 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) (r= 0) is:

ε E 1 J - 4 ε J 0 T J 0 Γ s χ ψ , ω J χ ψ ψ , ω J d s d ψ .

We find that εE1(J)=-4εJω3Γ~s(ω)π-2(Ω2-ω2)-2 ([5]), where Γ~s 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:

x t = - J t χ t , ω J ,         x 1 t = - 1 2 J t χ t , ω J ,         x 2 t = 1 2 J t χ t , ω J ,

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, d shows the law of change in the relative coordinate. A similar graph will correspond to a body colliding with a fixed massive obstacle.

Fig. 3. The laws of motion of colliding bodies

 The laws of motion of colliding bodies


 The laws of motion of colliding bodies


 The laws of motion of colliding bodies


 The laws of motion of colliding bodies


5. 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.


  1. Blum J. Exploring Ardulino. John Wily and Sons, Indianapolis, USA, 2013. [Search CrossRef]
  2. Mandelstam L. I. Lectures on the Theory of Oscillations. Nauka, Moscow,1972, (in Russian). [Search CrossRef]
  3. Nagaev R. F. Mechanical Processes with Repeated Damped Collisions. Nauka, Moscow, 1985, (in Russian). [Search CrossRef]
  4. Kobrinsky A. A. Mechanisms with Elastic Connections and Impact Systems. ILife Books Ltd., London, 1969. [Search CrossRef]
  5. Babitsky V. I., Krupenin V. L. Vibration of Strongly Nonlinear Discontinuous Systems. Springer-Verlag, Berlin, 2001. [Publisher]
  6. Burd Sh V., Krupenin V. L. Averaging in Quasi-conservative Systems. Beliy Veter, Moscow, 2016. [Search CrossRef]
  7. Rabotnov Yu N. Elements of Hereditary Solid Mechanics. Nauka, Moscow, 1977. [Search CrossRef]