Investigation of elliptical approximation in the model of the force of dry friction

K. Ragulskis1 , A. Bubulis2 , P. Paškevičius3 , A. Pauliukas4 , L. Ragulskis5

1Kaunas University of Technology, K. Donelaičio Str. 73, LT-44249, Kaunas, Lithuania

2Kaunas University of Technology, Mechatronics Institute, Studentų Str. 56, LT-51424, Kaunas, Lithuania

3Company “Vaivora”, Palemono Str. 2a, LT-52191, Kaunas, Lithuania

4Aleksandras Stulginskis University, Studentų Str. 11, LT-53361, Akademija, Kaunas District, Lithuania

5Vytautas Magnus University, Vileikos Str. 8, LT-44404, Kaunas, Lithuania

1Corresponding author

Mathematical Models in Engineering, Vol. 4, Issue 3, 2018, p. 151-156. https://doi.org/10.21595/mme.2018.20001
Received 30 April 2018; received in revised form 22 May 2018; accepted 30 May 2018; published 30 September 2018

Copyright © 2018 K. Ragulskis, et al. 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 Acknowledgements References
Cite this article
Views 34
Reads 10
Downloads 866
Abstract.

In the numerical procedures involving the force of dry friction some type of numerical approximation for this investigated phenomenon is used in the process of calculations. Elliptic approximations for the local transition region are proposed, and numerical results are presented in this paper. The smoothness and continuity of this approximation is demonstrated by several computational experiments.

Keywords: force of dry friction, elliptic representation, results of calculations, vibrations, visual representation.

1. Introduction

Engineering problems and investigations of vibrating systems using various numerical procedures involving the force of dry friction are important in practical applications. Applications of various types of numerical approximations for this phenomenon are common. In the presented investigation elliptic approximation of the local transition region is proposed. Computational experiments are performed for typical problems and numerical results are presented. This approximation has good smoothness and continuity properties. Thus, it is advantageous in various numerical procedures.

This paper is a continuation of investigations for numerical representation of the similar model for the force of dry friction presented in [1]. Piecewise linear approximations of the model for the force representing dry friction as well as numerical investigations of various nonlinear problems of vibrations of mechanical systems are presented in [2]. Essentially nonlinear dynamical systems and their analytical as well as computational investigations are presented in [3]. Basic models of vibrating systems and their investigations are presented in [4]. Contemporary methods of investigation of vibrating systems are described in [5]. Basic models of vibrating systems and their analysis are described in [6]. Vibrating dynamical systems and their contemporary applications are presented in [7].

Numerical procedures of integration of vibrating mechanical systems are presented in [8]. Related presentation of numerical procedures for integration of vibrating systems is available in [9]. Method of integration using finite elements in time is described in [10]. More detailed description of the numerical methods based on finite elements in time is available in [11]. Numerical procedures for integration of differential equations of first order are presented in [12].

Solutions of various basic problems involving the force of dry friction are given in [13].

Nonlinear effects are important in problems of surface cleaning and they are described in [14]. The process of cleaning is also investigated in [15]. Interactions of particles with surfaces are investigated in [16]. Interactions of particles with surfaces of larger particles are analysed in [17]. Measurements of adhesion of particles are described in [18]. Related experimental investigations of particles using atomic force microscopy are presented in [19]. The investigations of numerical representation of the investigated model of the force representing dry friction are important for their possible use in the models of surface cleaning.

2. Investigation of elliptic representation for the analyzed force representing dry friction

Further u denotes displacement of the mechanical system, H represents the approximate model for the investigated force representing dry friction, while the dot over the variable is used to denote the operation of differentiation performed with respect to the time variable t.

The numerical value of the variable representing velocity at the beginning of the local transition interval between the constant values of the force representing the investigated dry friction is denoted as u˙1 and the numerical value of the variable representing velocity at the end of the local transition interval between the constant values of the force representing the investigated dry friction is denoted as u˙2. The average value is calculated as:

(1)
u ˙ c = u ˙ 1 + u ˙ 2 2 .

The value of the approximation for the force representing dry friction at the beginning of the local transition interval interconnecting the constant values for the force representing dry friction is denoted as H1 and the value of the force representing dry friction at the end of the local transition interval interconnecting the constant values for the force representing dry friction is denoted as H2. The average value is calculated as:

(2)
H c = H 1 + H 2 2 .

The following quantities denoting one half of the difference of the previously defined values are introduced:

(3)
a = H 2 - H 1 2 ,
(4)
b = u ˙ 2 - u ˙ 1 2 .

The approximation of the transition region is based on the following elliptical equations:

(5)
H - H c 2 a 2 + u ˙ - u ˙ c - b 2 b 2 = 1 ,
(6)
H - H c 2 a 2 + u ˙ - u ˙ c + b 2 b 2 = 1 .

The following elliptic function is introduced:

(7)
E u ˙ 1 , u ˙ 2 , H 1 , H 2 , u ˙ = H c - a 1 - u ˙ - u ˙ c - b b 2 ,         u ˙ u ˙ c , H c + a 1 - u ˙ - u ˙ c + b b 2 ,           u ˙ > u ˙ c .

The force of dry friction is approximated as:

(8)
H = E - Δ , Δ , - h , h , u ˙ , u ˙ < Δ , h , u ˙ Δ ,         u ˙ > 0 , - h , u ˙ Δ ,         u ˙ < 0 ,

where h denotes the coefficient representing the dry friction, while Δ represents the width of the interval for transition between the required values.

The analyzed nonlinear mechanical system is represented by the following numerical model:

(9)
m u ¨ + c u ˙ + k u = P s i n ω t - H ,

where m is the notation used for the mass of the investigated system, c denotes the coefficient representing the force of viscous friction, k denotes the coefficient representing the stiffness of the investigated mechanical system, P denotes the amplitude of the force which excites the vibrating motion, ω denotes the frequency of the force performing excitation.

In the performed investigation values of the parameters of the analyzed nonlinear system were given the quantities: ω= 1, h= 3.2, P= 4, m= 1, c= 0.1, k= 1. In the numerical investigations zero initial conditions were assumed. Steady state motion was investigated. Two periods were represented.

Results for Δ= 0.8 are shown graphically in Fig. 1.

Some of the results for Δ=0.4 are shown graphically in Fig. 2.

Some of the results for Δ=0.2 are shown graphically in Fig. 3.

Fig. 1. Motion in steady state regime for wide transition region

Motion in steady state regime for wide transition region

a) Time history of displacement

Motion in steady state regime for wide transition region

b) Time history of velocity

Motion in steady state regime for wide transition region

c) Time history of acceleration

Motion in steady state regime for wide transition region

d) Time history of H

Motion in steady state regime for wide transition region

e) Dependence of H from velocity

Motion in steady state regime for wide transition region

f) Dependence of velocity from displacement

Motion in steady state regime for wide transition region

g) Dependence of acceleration from velocity

The presented graphical relationships for three typical widths of the transition region show the influence of the width of the interval representing transition to the obtained numerical results. The effect of the width of the interval representing transition is especially clearly seen in the graphical relationships involving the accelerations.

Fig. 2. Motion in steady state regime for transition region of medium width

Motion in steady state regime for transition region of medium width

a) Time history of acceleration

Motion in steady state regime for transition region of medium width

b) Dependence of acceleration from velocity

Fig. 3. Motion in steady state regime for narrow transition region

Motion in steady state regime for narrow transition region

a) Time history of acceleration

Motion in steady state regime for narrow transition region

b) Dependence of acceleration from velocity

3. Investigation of more complicated model for the analyzed force representing dry friction

In this model the investigated coefficient representing dry friction increases in the vicinity of the value of velocity equal to zero. This fact is taken into account in a more complicated model of the representation for the force approximating dry friction. Here it is usually considered that it is valid h/Δ=ha/Δa, where ha denotes the increase in the approximation of the variation of the coefficient representing the force of dry friction in the investigated region, while Δa defines the width of the interval representing the transition connecting the values of the investigated coefficient approximating dry friction h and h+ha. Thus ha=hΔaΔ.

In the numerical model the force of dry friction is approximated as:

(10)
H = E - Δ - Δ a , Δ + Δ a , - h - h Δ a Δ , h + h Δ a Δ , u ˙ , u ˙ < Δ + Δ a , E Δ + Δ a , Δ + 2 Δ a , h + h Δ a Δ , h , u ˙ , Δ + Δ a u ˙ < Δ + 2 Δ a ,         u ˙ > 0 , E - Δ - 2 Δ a , - Δ - Δ a , - h , - h - h Δ a Δ , u ˙ , Δ + Δ a u ˙ < Δ + 2 Δ a ,         u ˙ < 0 , h , u ˙ Δ + 2 Δ a ,         u ˙ > 0 , - h , u ˙ Δ + 2 Δ a ,         u ˙ < 0 .

In the performed investigation values of the parameters of the analyzed nonlinear system were given the quantities: ω= 1, h= 1.6, Δ= 0.8, Δa= 0.8, P= 4, m= 1, c= 0.1, k= 1. In the numerical investigations zero initial conditions were assumed. Steady state motion was investigated. Two periods were represented.

Results of calculations are shown graphically in Fig. 4.

Typical features of dynamic behavior of the nonlinear vibrating system with more complicated model of dry friction are reproduced by this elliptic representation for approximation of the effect of dry friction.

Fig. 4. Motion in steady state regime for more complicated model of the force approximating dry friction

Motion in steady state regime for more complicated model of the force approximating dry friction

a) Time history of displacement

Motion in steady state regime for more complicated model of the force approximating dry friction

b) Time history of velocity

Motion in steady state regime for more complicated model of the force approximating dry friction

c) Time history of acceleration

Motion in steady state regime for more complicated model of the force approximating dry friction

d) Time history of H

Motion in steady state regime for more complicated model of the force approximating dry friction

e) Dependence of H from velocity

Motion in steady state regime for more complicated model of the force approximating dry friction

f) Dependence of velocity from displacement

Motion in steady state regime for more complicated model of the force approximating dry friction

g) Dependence of acceleration from velocity

4. Conclusions

In engineering investigations of vibrating systems various types of nonlinear behavior are to be used in numerical procedures. It is important to use numerically acceptable approximations of this behaviour. Typical often used nonlinearity is dry friction. Usually some type of approximation for this phenomenon is used. Here elliptic approximation of the local transition region is proposed. Numerical results are presented in detail. This approximation can be noted for its high smoothness and continuity. The influence of the width of the interval interconnecting the investigated values to the results obtained from calculations is shown in the obtained graphical relationships.

Numerical analysis of the mechanical model with more advanced model of the approximation for the investigated force representing dry friction has been presented. In this model the increase of the investigated coefficient representing the force approximating the effect of dry friction in the vicinity of the value of velocity equal to zero is taken into account. Typical features of dynamic behavior of the nonlinear vibrating system with more complicated model of dry friction are reproduced by this elliptic representation for the force approximating dry friction. Based on the presented numerical results the elliptic approximation is recommended for engineering calculations involving the action of force representing the effect of dry friction.

Acknowledgements

The authors thank the anonymous reviewers for providing new ideas and comments. On the basis of them it was possible to substantially improve the material of the paper.

References

  1. Ragulskis K., Paškevičius P., Bubulis A., Pauliukas A., Ragulskis L. Improved numerical approximation of dry friction phenomena. Mathematical Models in Engineering, Vol. 3, Issue 2, 2017, p. 106-111. [Publisher]
  2. Levy S., Wilkinson J. P. D. The Component Element Method in Dynamics with Application to Earthquake and Vehicle Engineering. McGraw-Hill, New York, 1976. [Search CrossRef]
  3. Ragulskienė V. Vibro-Shock Systems (Theory and Applications). Mintis, Vilnius, 1974, (in Russian). [Search CrossRef]
  4. Bolotin V. V. Vibrations in Engineering. Handbook, Vol. 1, Mashinostroienie, Moscow, 1978, (in Russian). [Search CrossRef]
  5. Inman D. J. Vibration with Control, Measurement, and Stability. Prentice-Hall, New Jersey, 1989. [Search CrossRef]
  6. Lalanne M., Berthier P., Der Hagopian J. Mechanical Vibrations for Engineers. John Wiley and Sons, New York, 1984. [Search CrossRef]
  7. Thomson W. T. Theory of Vibration with Applications. Prentice-Hall, New Jersey, 1981. [Search CrossRef]
  8. Bathe K. J. Finite Element Procedures in Engineering Analysis. Prentice-Hall, New Jersey, 1982. [Search CrossRef]
  9. Bathe K. J., Wilson E. L. Numerical Methods in Finite Element Analysis. Stroiizdat, Moscow, 1982, (in Russian). [Search CrossRef]
  10. Zienkiewicz O. C. The Finite Element Method in Engineering Science. Mir, Moscow, 1975, (in Russian). [Search CrossRef]
  11. Zienkiewicz O. C., Morgan K. Finite Elements and Approximation. Mir, Moscow, 1986, (in Russian). [Search CrossRef]
  12. Segerlind L. J. Applied Finite Element Analysis. Mir, Moscow, 1979, (in Russian). [Search CrossRef]
  13. Sumbatov A. S., Yunin Ye. K. Selected Problems of Mechanics of Systems with Dry Friction. Physmathlit, Moscow, 2013, (in Russian). [Search CrossRef]
  14. Chahine G. L., Kapahi A., Choi J.-K., Hsiao Ch.-T. Modeling of surface cleaning by cavitation bubble dynamics and collapse. Ultrasonics Sonochemistry, Vol. 29, 2016, p. 528-549. [Publisher]
  15. Witte A. K., Bobal M., David R., Blättler B., Schoder D., Rossmanith P. Investigation of the potential of dry ice blasting for cleaning and disinfection in the food production environment. LWT – Food Science and Technology, Vol. 75, 2016, p. 735-741. [Publisher]
  16. Petean P. G. C., Aguiar M. L. Determination of the adhesion force between particles and rough surfaces. Powder Technology, Vol. 274, 2015, p. 67-76. [Publisher]
  17. Cui Y., Sommerfield M. Forces on micron-sized particles randomly distributed on the surface of larger particles and possibility of detachment. International Journal of Multiphase Flow, Vol. 72, 2015, p. 39-52. [Publisher]
  18. Jiang Y., Turner K. T. Measurement of the strength and range of adhesion using atomic force microscopy. Extreme Mechanics Letters, Vol. 9, Issue 1, 2016, p. 119-126. [Publisher]
  19. Kumar N., Zhao C., Klaassen A., Van Den Ende D., Mugele F., Siretanu I. Characterisation of the surface charge distribution on kaolinite particles using high resolution atomic force microscopy. Geochimica et Cosmochimica Acta, Vol. 175, 2016, p. 100-112. [Publisher]

Cited By

Mathematical Models in Engineering
K. Ragulskis, A. Bubulis, R. Maskeliūnas, P. Paškevičius, L. Ragulskis
2018