Investigation of elliptical approximation in the model of the force of dry friction
K. Ragulskis^{1} , A. Bubulis^{2} , P. Paškevičius^{3} , A. Pauliukas^{4} , L. Ragulskis^{5}
^{1}Kaunas University of Technology, K. Donelaičio Str. 73, LT44249, Kaunas, Lithuania
^{2}Kaunas University of Technology, Mechatronics Institute, Studentų Str. 56, LT51424, Kaunas, Lithuania
^{3}Company “Vaivora”, Palemono Str. 2a, LT52191, Kaunas, Lithuania
^{4}Aleksandras Stulginskis University, Studentų Str. 11, LT53361, Akademija, Kaunas District, Lithuania
^{5}Vytautas Magnus University, Vileikos Str. 8, LT44404, Kaunas, Lithuania
^{1}Corresponding author
Mathematical Models in Engineering, Vol. 4, Issue 3, 2018, p. 151156.
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
JVE Conferences
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 ${\dot{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 ${\dot{u}}_{2}$. The average value is calculated as:
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 ${H}_{1}$ 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 ${H}_{2}.$ The average value is calculated as:
The following quantities denoting one half of the difference of the previously defined values are introduced:
The approximation of the transition region is based on the following elliptical equations:
The following elliptic function is introduced:
The force of dry friction is approximated as:
where $h$ denotes the coefficient representing the dry friction, while $\mathrm{\Delta}$ represents the width of the interval for transition between the required values.
The analyzed nonlinear mechanical system is represented by the following numerical model:
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, $\omega $ 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: $\omega =$ 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 $\mathrm{\Delta}=$ 0.8 are shown graphically in Fig. 1.
Some of the results for $\mathrm{\Delta}=\text{0.4}$ are shown graphically in Fig. 2.
Some of the results for $\mathrm{\Delta}=\text{0.2}$ are shown graphically in Fig. 3.
Fig. 1. Motion in steady state regime for wide transition region
a) Time history of displacement
b) Time history of velocity
c) Time history of acceleration
d) Time history of $H$
e) Dependence of $H$ from velocity
f) Dependence of velocity from displacement
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
a) Time history of acceleration
b) Dependence of acceleration from velocity
Fig. 3. Motion in steady state regime for narrow transition region
a) Time history of acceleration
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/\mathrm{\Delta}={h}_{a}/{\mathrm{\Delta}}_{a}$, where ${h}_{a}$ denotes the increase in the approximation of the variation of the coefficient representing the force of dry friction in the investigated region, while ${\mathrm{\Delta}}_{a}$ defines the width of the interval representing the transition connecting the values of the investigated coefficient approximating dry friction $h$ and $h+{h}_{a}.$ Thus ${h}_{a}=h\frac{{\mathrm{\Delta}}_{a}}{\mathrm{\Delta}}.$
In the numerical model the force of dry friction is approximated as:
In the performed investigation values of the parameters of the analyzed nonlinear system were given the quantities: $\omega =\mathrm{}$1, $h=$ 1.6, $\mathrm{\Delta}=$ 0.8, ${\mathrm{\Delta}}_{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
a) Time history of displacement
b) Time history of velocity
c) Time history of acceleration
d) Time history of $H$
e) Dependence of $H$ from velocity
f) Dependence of velocity from displacement
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
 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. 106111. [Publisher]
 Levy S., Wilkinson J. P. D. The Component Element Method in Dynamics with Application to Earthquake and Vehicle Engineering. McGrawHill, New York, 1976. [Search CrossRef]
 Ragulskienė V. VibroShock Systems (Theory and Applications). Mintis, Vilnius, 1974, (in Russian). [Search CrossRef]
 Bolotin V. V. Vibrations in Engineering. Handbook, Vol. 1, Mashinostroienie, Moscow, 1978, (in Russian). [Search CrossRef]
 Inman D. J. Vibration with Control, Measurement, and Stability. PrenticeHall, New Jersey, 1989. [Search CrossRef]
 Lalanne M., Berthier P., Der Hagopian J. Mechanical Vibrations for Engineers. John Wiley and Sons, New York, 1984. [Search CrossRef]
 Thomson W. T. Theory of Vibration with Applications. PrenticeHall, New Jersey, 1981. [Search CrossRef]
 Bathe K. J. Finite Element Procedures in Engineering Analysis. PrenticeHall, New Jersey, 1982. [Search CrossRef]
 Bathe K. J., Wilson E. L. Numerical Methods in Finite Element Analysis. Stroiizdat, Moscow, 1982, (in Russian). [Search CrossRef]
 Zienkiewicz O. C. The Finite Element Method in Engineering Science. Mir, Moscow, 1975, (in Russian). [Search CrossRef]
 Zienkiewicz O. C., Morgan K. Finite Elements and Approximation. Mir, Moscow, 1986, (in Russian). [Search CrossRef]
 Segerlind L. J. Applied Finite Element Analysis. Mir, Moscow, 1979, (in Russian). [Search CrossRef]
 Sumbatov A. S., Yunin Ye. K. Selected Problems of Mechanics of Systems with Dry Friction. Physmathlit, Moscow, 2013, (in Russian). [Search CrossRef]
 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. 528549. [Publisher]
 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. 735741. [Publisher]
 Petean P. G. C., Aguiar M. L. Determination of the adhesion force between particles and rough surfaces. Powder Technology, Vol. 274, 2015, p. 6776. [Publisher]
 Cui Y., Sommerfield M. Forces on micronsized particles randomly distributed on the surface of larger particles and possibility of detachment. International Journal of Multiphase Flow, Vol. 72, 2015, p. 3952. [Publisher]
 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. 119126. [Publisher]
 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. 100112. [Publisher]
Cited By
Mathematical Models in Engineering
K. Ragulskis, A. Bubulis, R. Maskeliūnas, P. Paškevičius, L. Ragulskis

2018
