Selforganization and evolution in dynamic friction systems
Vilor L. Zakovorotny^{1} , Valery E. Gvindjiliya^{2}
^{1, 2}Don State Technical University (DSTU), RostovonDon, 344000, Russia
^{2}Corresponding author
Journal of Vibroengineering, Vol. 23, Issue 6, 2021, p. 14181432.
https://doi.org/10.21595/jve.2021.22033
Received 6 May 2021; received in revised form 28 July 2021; accepted 5 August 2021; published 30 September 2021
The dynamic features of the friction process of the contacting pairs affect their durability and functional suitability. Therefore, the problems of their dynamics and vibrational activity always attract the attention of researchers in the field of tribology. It is known that selforganization processes take place during friction, such as the formation of selective transfer phenomena during friction, formation of equilibrium roughness. But known studies do not take into account changes in the dynamic properties of the system in the process of its functioning. The primary sources of selforganization are energetic in nature, but the relationship of these phenomena with the dynamic properties of the system has not been explained, and their modeling has not yet been performed. The paper presents the results of experimental and mathematical modeling of evolutionary changes of friction system due to the power of irreversible transformations of energy input to friction zone, which characterizes the goal of research. In order to achieve the goal, the paper deals with the mathematical modeling of the processes of selforganization in a dynamic friction system in combination with the analysis of the equilibrium of the system, its stability, and the manifolds formed in the vicinity of the equilibrium. It is shown that the proposed mathematical models are integrodifferential functionally coupled equations that describe selforganization. Examples are given of dynamic restructuring of the friction system, including bifurcations of the attracting sets of its deformation displacements along the trajectory of the friction forces work. The paper briefly presents the results of the experimental study of the vibration characteristics of the contacting pairs during the development of their wear by the example of two systems: steelbrass and steelglycerinbrass. It is shown that in the system of friction, new properties always establish naturally due to the work of forces. These properties are determined by the evolutionary change of the parameters of the dynamic coupling and the attractive sets of deformation displacements. The study and experiments show that measurement of parameters available for estimation of some coordinates in the dynamic friction system allows one to solve the problem of nondecoupled dynamic monitoring of the state of friction units and also opens new ways of improvement of friction units.
 Dynamic model of cutting forces in the coordinates of the state of the cutting system is developed, taking into account irreversible transformations in the contact zone of interacting surfaces
 Experimental mathematical modeling of evolutionary reorganization in a dynamic friction system is carried out
 The reorganization of the dynamic cutting system in the process of evolution is analyzed using the example of phase trajectories of the dynamic cutting system
Keywords: dynamic system, tribosystem, evolution, selforganization, friction.
1. Introduction
The synergetic paradigm of evolution and selforganization that was developed and widely used in the second half of the 20th century [1][3] has found its way into tribology as well [4][9]. It was established that the process of friction is described in terms of interconnected physical, physicchemical and chemical phenomena that lead to various changes in the transition zone between two interacting surfaces. We will refer to this zone as tribosphere. If we follow the synergetic paradigm, then according to the principle of expansion of statespace dimension, first of all, it is necessary to describe the interactions between mechanical subsystems and the tribosphere [10]. Therefore, the forces that are formed on the surfaces in contact should be represented in the state coordinates of the system. The mathematical description of the forces in the state coordinates defines the properties of the link that connects interacting subsystems. It is important to note that the consideration of the contact interaction forces in the function of the state coordinates fundamentally changes the paradigm of friction analysis. Friction is already considered as a dynamic system in which the forces that depend on the coordinates of the state actually form an intrasystem regulator [11]. Therefore, the parameters of this regulator, i.e., the dynamic coupling formed by the friction zone, determine the asymptotic stability of the tribosystem. When stability is lost in such a system, various attractors (limit cycles, invariant tori, and chaotic attractors) must form in the vicinity of the equilibrium. When studying tribosystems, we gave priority to the analysis of selfoscillations [10][19]. In recent years, several studies have been published that consider the conditions for the appearance of tori and strange (chaotic) attractors [10], [19]. Various ways of selfexcitation are considered, including those associated with the formation of a timedelay argument [20]. In recent years, particularly great interest has been paid to the study of the stability and properties of the attracting sets of deformation displacements in friction systems during cutting [21][28]. At the same time, it is shown that vibrations affect the wear rate of the contacted pairs [29][31]. It was the change in the vibration properties of the friction unit depending on the state of the tribomedium that led to the development of the direction in tribotechnology, which consists in the creation of systems for vibration diagnostics of friction processes [30][33].
However, in all known works on the dynamics of the friction process the dynamic link formed by the process is deemed to be set and constant during the life of the system. In our opinion, the next phase of research should take into account the changes in representation of forces in the state coordinates through modeling how irreversible transformations in the zone of the contact between interacting surfaces influence the tribosphere properties. These irreversible transformations in particular cause evolutionary changes of the properties of the friction system leading to a dynamic reorganization of the system during its functioning.
We emphasize that the power of irreversible transformations stimulates changes in the parameters of the dynamic coupling over time, changes the dynamic properties of the contact and the output characteristics of the friction unit, for example, the wear rate of the contacting pairs. The purpose of the research presented in the article is experimentalmathematical modeling of evolutionary reorganization in a dynamic friction system. It is the study of the evolution of the properties of friction due to the phase trajectory of the power of irreversible transformations in the friction zone that fundamentally distinguishes the materials presented in the article from the known ones.
2. Generalized mathematical model of the evolutionary dynamic system of friction. Method of its analysis
In the study of machine dynamics, it is usually represented in the form of a finitedimensional space structure where separate generalized masses are linked with linear elastic dissipative connections [10], [14]. The principal properties of this system can be seen if we analyze systems with the following features (Fig. 1).
Deformational displacements $\mathbf{X}=\{{X}_{1},{X}_{2}{\}}^{T}$ and $\mathbf{Y}=\{{Y}_{1},{Y}_{2}{\}}^{T}$ of the two interacting subsystems are considered in a plane normal to the axis of rotation of the workpiece. 2. The object of the study is the contact between a spherical indenter of radius $r$ and a round specimen of radius $R$. Note that $R\gg r$. The specimen rotates at constant angular velocity ${\mathrm{\Omega}}_{0}$ represented in $\left[{c}^{1}\right]$. 3. There are no torsional deformations of the specimen and the indenter. Therefore, generated forces of the contact interaction are applied to the centers of mass of the interacting subsystems and deformational displacements occur in the directions shown on Fig. 1. Here forces ${F}^{\left(X\right)}$ and ${F}^{\left(Y\right)}$ are equal and opposite. 4. Width ${\mathrm{\Lambda}}_{0}$ of the tribosphere can be considered constant and its properties  isotropic along the perimeter of rotation. 5. Forces that are formed at the contact of the surfaces depend on the depth of penetration into the tribosphere as well as on speeds of the relative movement of the surfaces.
Thus, forces are represented by dependence (Fig. 1):
where ${\mathbf{\Delta}}_{\mathrm{0,1}}$ is the preset value of the indenter without elastic deformational displacements; ${V}_{0}$ is relative sliding velocity without velocities of elastic deformational displacements; ${p}_{1},...{p}_{n}$ are parameters included in the forces model. If Eq. (1) are defined, then the link that joins subsystems together into one dynamic system of friction is set. The above statements make it possible to define an equation for the friction system dynamics as follows:
where:
${h}^{\left(X\right)}=\left[\begin{array}{ll}{h}_{\mathrm{1,1}}^{\left(X\right)}& {h}_{\mathrm{2,1}}^{\left(X\right)}\\ {h}_{\mathrm{1,2}}^{\left(X\right)}& {h}_{\mathrm{2,2}}^{\left(X\right)}\end{array}\right],{h}^{\left(Y\right)}=\left[\begin{array}{ll}{h}_{\mathrm{1,1}}^{\left(Y\right)}& {h}_{\mathrm{2,1}}^{\left(Y\right)}\\ {h}_{\mathrm{1,2}}^{\left(Y\right)}& {h}_{\mathrm{2,2}}^{\left(Y\right)}\end{array}\right],$
${c}^{\left(X\right)}=\left[\begin{array}{ll}{c}_{\mathrm{1,1}}^{\left(X\right)}& {c}_{\mathrm{2,1}}^{\left(X\right)}\\ {c}_{\mathrm{1,2}}^{\left(X\right)}& {c}_{\mathrm{2,2}}^{\left(X\right)}\end{array}\right],{c}^{\left(Y\right)}=\left[\begin{array}{ll}{c}_{\mathrm{1,1}}^{\left(Y\right)}& {c}_{\mathrm{2,1}}^{\left(Y\right)}\\ {c}_{\mathrm{1,2}}^{\left(Y\right)}& {c}_{\mathrm{2,2}}^{\left(Y\right)}\end{array}\right],$
are matrices of inertia, of velocity, of other coefficients. That is why, if the reaction from the friction process of subsystem is not taken into account, $X$ and $Y$ are, first, asymptomatically stable and, second, autonomous. It is Eq. (1) that change the properties of the system and transform autonomous subsystems into a single dynamic system of friction. Currently existing dynamic models of the friction process differ in the ways Eq. (1) are presented [10], [14]. They are aimed at exposing various physical processes in the friction zone, for example, the decrease of friction forces as the speed of relative sliding increases [14]. At the same time, as a rule, scalar (not special) models of indenter oscillations are being considered. These models enable the explanation some issues with the loss of stability and appearance of autooscillations. Moreover, parameters of the model are deemed to be set and unchanging during the operation. Therefore, the most important question of the synergetic paradigm on the evolutionary reorganization is not examined at all. The basic concept of the evolutionary reorganization of the system means that the changes of tribosphere properties are associated with the phase trajectory of the power of irreversible transformations of input energy in the contact zone for the work done. Then, taking into account Eqs. (1) and (2) we obtain a dynamics evolutionary equation:
where ${w}_{{i}_{i}}(A\varsigma )=\mathrm{e}\mathrm{x}\mathrm{p}[({T}_{i,1}{)}^{1}(A\varsigma )]+{\beta}_{i}\mathrm{e}\mathrm{x}\mathrm{p}[(A\varsigma )\left({T}_{i,2}{)}^{1}\right]$ are kernels of integral operators, with that ${T}_{i.1}$, ${T}_{i.,2}$ are constant works with dimension $kgm$, these parameters describe the evolutionary heredity of trajectories as the work is done; ${\alpha}_{i}$ are parameters that determine the connection of the trajectory of power for the work done; ${p}_{0,i}$ is the initial parameter value. The system Eq. (3) properties are investigated in several stages.
At the first stage stationary evolutionary trajectory (${X}_{1}^{*}\left(A\right)$, ${X}_{2}^{*}\left(A\right)$, ${Y}_{1}^{*}\left(A\right)$, ${Y}_{2}^{*}\left(A\right)$) is analyzed. It constitutes the solution to integral equation:
At the second stage we analyze the stability of trajectories (${X}_{1}^{*}\left(A\right)$, ${X}_{2}^{*}\left(A\right)$, ${Y}_{1}^{*}\left(A\right)$, ${Y}_{2}^{*}\left(A\right)$) corresponding to parameters ${p}_{1}\left(A\right),...{p}_{n}\left(A\right)$. If they are stable, then they will be attractors for system Eq. (3). Otherwise, it will be necessary to analyze attractors in the neighborhood of the stationary evolutionary trajectory, their corresponding variations of contact interaction forces, changes of work and power of irreversible transformations as well as to adjust parameters of the dynamic link, that is to say the very stationary evolutionary trajectories. For this purpose, it is convenient to analyze in advance the properties of system Eq. (3) and to construct system bifurcation diagrams in the zone of identification of evolutionary parameters. However, even if there are no system bifurcations in the course of evolution, it is still possible to evaluate tribosphere properties from evolutionary diagrams of roots of the characteristic polynomial linearized in the neighborhood of the system equilibrium.
Earlier we have put this concept forward as a basis for constructing systems of dynamic diagnostics of the friction site state based on the use of autoregressive spectral analysis of the observed vibrational sequences [10], [34][36].
3. Example of the system and its properties
Let us examine system Eq. (3) where there are no deformations of the specimen present and let us take into account the following main features of modeling of forces in state coordinates [10], [14], [38].
1) The tribosphere exhibits gradient functional services [14], [38]. Therefore, as surfaces come closer, first, the forces of contact interaction increase and, second, normal and tangential forces are redistributed. Their thorough analysis shows that in certain cases there are potential barriers in these characteristics associated with various physical phenomena, for example, the destruction of the servovite film (protective tribofilm) [10], [37], [38]. Let us call changes in contact interaction forces caused by surfaces coming closer together approximation functions. If approximation functions are set, then for every gap we can obtain important characteristics of the friction site. For example, the ratio of tangential force to normal force in a given equilibrium point is usually interpreted as coefficient of friction. From the changes in the gap it is possible to establish force gradients that affect the existence condition of external friction, etc. In the dynamic system these gradients define the dynamic rigidity of the contact. In this paper we will only refer to the most common case where with the narrowing of the gap the forces are increasing steadily. The main properties of these functions can be described in rather simple expressions:
where $\rho $, $\upsilon $ are parameters having meaning of pressure and density respectively. They have dimensions [kg/mm^{2}] and [kg/mm^{3}]. ${X}_{\mathrm{0,1}}$ is the distance between the surface of the tribosphere and the apex of the indentor with the assumption that there are no elastic deformations.
Fig. 1. Diagram of the basic model of the dynamic friction model
2) Research suggests that, firstly, contact interaction forces do not appear instantly. Secondly, there is some delay between variations of normal and tangential components of forces. Thirdly, when examining steady disturbed periodic motions phase shifts between variations of normal and tangential components of the forces are amplitude and frequencydependent [38], [39]. Finally, the trajectory covered by the indenter in relation to the specimen in the speed range (0.13.0) m/s remains practically unchanged once a new force balance is established. Therefore, the amount of delay depends on the relative sliding speed.
3) When the relative sliding speed is rising with the normal force unchanged the tangential component is decreasing which is associated with the temperaturespeed factor. The above specific features allow the presentation of contact interaction forces in the following way:
${F}_{2}({X}_{1},{X}_{2})=\left\{\begin{array}{l}\upsilon \left\{1k\mathrm{e}\mathrm{x}\mathrm{p}\left[a\left(V+\frac{d{X}_{2}}{dt}\right)\right]\right\}({X}_{\mathrm{0,1}}{Y}_{2}{)}^{3},{X}_{\mathrm{0,1}}{X}_{1},\\ 0,\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}{X}_{\mathrm{0,1}}\le {X}_{1},\end{array}\right.$
$\frac{{T}_{1}d{Y}_{1}}{dt}={X}_{1}{Y}_{1},\frac{{T}_{2}d{Y}_{2}}{dt}={X}_{1}{Y}_{2},$
where ${T}_{1},\dots ,{T}_{2}$ are time constants that characterize the delayed transformation of elastic deformational displacement changes into forces; $k$ is a nondimensional parameter $k<1$; $\u0430$ is a parameter that described the steepness of friction force decrease when relative sliding speed is increasing [s/m]. In Eq. (6) relation ${T}_{1}<{T}_{2}$ is fulfilled. And ${T}_{1}={l}_{1}/V$; ${T}_{2}={l}_{2}/V$ is true (${l}_{1}$ is an efficient trajectory needed to establish the normal component and ${l}_{2}$ – the tangential component of forces). State coordinates in Eq. (6) have a limited range of definition $\left(V+\frac{d{X}_{2}}{dt}\right)>0$, ${X}_{1}<{X}_{\mathrm{0,1}}$. Nowadays the paper is published that takes into account delay of friction forces in depend on variation of deformational displacement, but in modeling of dynamic system the author used the model of pure time delay [20].
4. Equilibrium properties in the dynamic friction system
The system equilibrium is defined by conditions: ${X}_{1}^{}={X}_{1}^{*}=const$, ${X}_{2}^{}={X}_{2}^{*}=const$. From Eqs. (2) and (6) we obtain an equation for calculation of the equilibrium. If ${X}_{\mathrm{0,1}}=\u0441onst$ and $V=const$ for ${X}_{1}^{*}$:
where ${\vartheta}_{V}=\upsilon [1k\mathrm{e}\mathrm{x}\mathrm{p}(aV\left)\right]$; ${c}_{\mathrm{\Sigma}}=[{c}_{\mathrm{1,1}}{c}_{\mathrm{2,2}}({c}_{\mathrm{1,2}}{)}^{2}\left]\right({c}_{\mathrm{2,2}}{)}^{1}$. After calculating ${X}_{1}^{*}$ we establish ${X}_{2}^{*}=\left\{\rho \right({X}_{\mathrm{0,1}}{X}_{1}^{*}{)}^{2}{c}_{\mathrm{1,1}}{X}_{1}^{*}\left\}\right({c}_{\mathrm{2,1}}{)}^{1}$. Eq. (7) has three solutions. However, in the range of definition of the parameters it always has only one real solution because $\rho ({X}_{\mathrm{0,1}}{X}_{1}^{*}{)}^{2}\gg {\vartheta}_{V}({X}_{\mathrm{0,1}}{X}_{1}^{*}{)}^{3}{c}_{\mathrm{2,1}}/{c}_{\mathrm{2,2}}$ and ${X}_{\mathrm{0,1}}>{X}_{1}^{*}$. Stability of the equilibrium is of fundamental importance. To analyze it is necessary to examine linearized equation in variations with respect to (${X}_{1}^{*}$, ${X}_{2}^{*}$). After substitution of ${X}_{1}^{}={X}_{1}^{*}+{x}_{1}\left(t\right)$ and ${X}_{2}^{}={X}_{2}^{*}+{x}_{2}\left(t\right)$ from Eqs. (1) and (3) we obtain a linearized equation in variations:
where:
Since ${X}^{*}=\{{X}_{1}^{*},{X}_{2}^{*}\}=const$, then Eq. (8) has constant parameters. Then stability is defined with solutions of characteristic equation:
$\left.\begin{array}{ll}\iff & {h}_{\mathrm{2,1}}p+{c}_{\mathrm{2,1}}\\ \iff & m{p}^{2}+\{{h}_{\mathrm{2,2}}\vartheta ka({X}_{\mathrm{0,1}}{X}_{1}^{*}{)}^{3}\mathrm{e}\mathrm{x}\mathrm{p}(aV)\}p+{c}_{\mathrm{2,2}}\end{array}\right]=0.$
To achieve this the known criteria can be used [40]. To ascertain the mechanisms of stability loss in the system the structure of matrices ${h}^{\left(\mathrm{\Sigma}\right)}$ and ${c}^{\left(\mathrm{\Sigma}\right)}$ can be analyzed; they are not symmetric due to the reaction induced by the friction process [41][43]. They can be represented in the form of symmetric and skewsymmetric components ${h}^{\left(\mathrm{\Sigma}\right)}={h}_{C}^{\left(\mathrm{\Sigma}\right)}+{h}_{K}^{\left(\mathrm{\Sigma}\right)}$ and ${c}^{\left(\mathrm{\Sigma}\right)}={c}_{C}^{\left(\mathrm{\Sigma}\right)}+{c}_{K}^{\left(\mathrm{\Sigma}\right)}$. Therefore, due to the process reaction potential velocity and elastic links are formed (defined by matrices ${h}_{C}^{\left(\mathrm{\Sigma}\right)}$ and ${c}_{C}^{\left(\mathrm{\Sigma}\right)}$), as well as gyroscopic and circulatory forces (defined by matrices ${h}_{K}^{\left(\mathrm{\Sigma}\right)}$ and ${c}_{K}^{\left(\mathrm{\Sigma}\right)}$). If ${h}_{C}^{\left(\mathrm{\Sigma}\right)}$ is positively definite, then it is a dissipation matrix and given the positive definiteness of the symmetric part of the elasticity matrix the system has an asymptomatically stable equilibrium point. In the case when ${h}_{C}^{\left(\mathrm{\Sigma}\right)}$ is negatively definite, then it forms accelerating forces. Then the equilibrium point becomes unstable and gyroscopic forces cannot stabilize the equilibrium. Thus, the first main mechanism of the stability loss is determined by the transformation of the symmetric part of the speed coefficient matrix from positive definite to negative definite. Being positive definite matrix ${c}_{C}^{\left(\mathrm{\Sigma}\right)}$ forms forces directed towards the equilibrium; matrix ${c}_{K}^{\left(\mathrm{\Sigma}\right)}$ forms forces that are orthogonal to the direction of deformational displacements. They cause the formation of circular trajectories in relation to the equilibrium point. If the symmetric part becomes negative definite, then the system equilibrium becomes unstable and periodic variations start to develop in the system in the form of precession oscillations. The above statements complement and generalize the known studies in the field of the friction process stability. First, the system stability depends on control parameters ${X}_{\mathrm{0,1}}$ and $V$. Second, stability loss is crucially affected by delayed arguments ${T}_{1}$ and ${T}_{2}$ and parameter $a$ that influence matrix ${h}_{C}^{\left(\mathrm{\Sigma}\right)}$. Third, stability loss mechanisms can be examined only using a spatial model of the system. Gyroscopic and circulatory forces do not occur in scalar models.
5. Attractors that are formed in the neighborhood of the equilibrium
If the equilibrium is unstable, attractors form in its neighborhood [23]. The system builds a model of two interconnected oscillating circuits that have two main sources of selfexcitation. One is caused by the delay in forces, the other is due to the existence of a region where the increase in speed corresponds to the decrease in the friction force. Therefore, three states can exist in the system: asymptomatically stable equilibrium point, asymptomatically stable limit cycle and twodimensional invariant torus. Let us examine a system for which parameters of the indenter dynamic model are given in Table 1 and parameters of the dynamic link formed by the friction process – in Table 2. The system in question is steelbronze. The bronze indenter is a sphere 8.0 mm in diameter. The specimen is made of steel. Friction with no lubricant. The methodology for parameter identification was detailed in the previous publications [13][15], [23].
Table 1. Matrix parameters of the indenter dynamic model
$m$, [kg⋅s^{2}/mm]

$h$, [kg⋅s/mm]

$c$, [kg/mm]

$\left[\begin{array}{cc}3.0\cdot 1{0}^{3}& 0\\ 0& 3.0\cdot 1{0}^{3}\end{array}\right]$

$\left[\begin{array}{ll}25& 2.0\\ 2.0& 12\end{array}\right]$

$\left[\begin{array}{cc}3000& 50\\ 50& 360\end{array}\right]$

Table 2. Parameters of the dynamic link formed by the friction process
$\rho $, [kg/mm^{2}]

$\upsilon $, [kg/mm^{3}]

$a$, [s/mm]

$k$

${l}_{1}$ [mm]

${l}_{2}$ [mm]

500

600

1200

0.6

5

7

The range of definition of the indenter penetration into the tribosphere in the steady state is $\left({X}_{\mathrm{0,1}}{X}_{1}^{*}\right)\in \left(25;150\right)$mkm. The range of variation of dynamic rigidity of the contact at the site of definition of the indenter penetration into the tribosphere is $\partial {F}_{1}/\partial {X}_{1}\in \left(50;1200\right)$ kg/mm and $\partial {F}_{2}/\partial {X}_{1}\in \left(2.5;150\right)$kg/mm. As the surfaces get closer to each other, for these parameters the coefficient of friction increases from ${k}_{TP.}=$ 0.06 to ${k}_{TP.}=$ 0.18. Without going into too much detail, let us analyze the bifurcation diagram in the plane of these parameters (Fig. 2). The shaded area on Fig. 2 corresponds to an asymptomatically stable equilibrium, the hatched area “DОE” is where a twodimensional invariant torus is formed. Dpartition areas are also shown. They divide the parameter space into “ВОС”, “СОF”, “FOA” areas and a shaded area. They differ in the number of roots of the characteristic polynomial of linearized equations in the right complex semiplane. In areas “ВОС” and “FOA” there is one pair of complex conjugate roots with positive real roots, area “СОF” has two pairs of such roots. However, the boundaries within which the twodimensional torus is formed are displaced on the bifurcation diagram (Fig. 2). Line “СОF” transforms into line “DOE”. Time charts for this case are shown in Fig. 3. It can be clearly seen that there is some pulling of high frequency oscillations towards low frequency ones.
Fig. 2. Bifurcation diagram in plane $(\rho ,\upsilon )$: “BОD” is an area of the orbitally asymptomatically stable limit cycle at high frequency; “EОA” is an area of the forming of the orbitally asymptomatically stable limit cycle at low frequency
The variation of parameters of matrix $c$ results in shifting of roots of the characteristic polynomial Eq. (9). The roots can get closer to each other. In this case the area where the twodimensional torus is formed degenerates.
Fig. 3. An example of transformation of a twofrequency process into a singlefrequency process
6. Reorganization of the dynamic system in the process of evolution.
Since there are no deformations of the specimen in the studies example, then taking into account Eqs. (3) and (6) it is necessary to use the system:
where ${p}_{1}=\rho \left(A\right)$, ${p}_{2}=\upsilon \left(A\right)$, ${p}_{3}={T}_{1}\left(A\right)$,${p}_{4}={T}_{2}\left(A\right)$, ${p}_{5}=k\left(A\right)$,${p}_{6}=a\left(A\right)$. System Eq. (10) is written for area ${X}_{\mathrm{0,1}}>{X}_{1}$. Let us give an example of the evolutionary reorganization of the system, parameters of which are specified in Tables 1 and 2. Parameters of kernels of integral operators for a relative sliding speed of 0.5 m/s are shown in Table 3.
Table 3. Parameters of kernels of integral operators
${p}_{i}$

${p}_{i,0}$

${\alpha}_{i}$

${\beta}_{i}^{}$

${T}_{i}^{\left(1\right)}$

${T}_{i}^{\left(2\right)}$

${p}_{1}=\rho $

500 [kg/(mm^{2})]

2×10^{7} [s/(kg(mm)^{4})]

0.7

2.99×10^{3} [kgm]

5.54×10^{3} [kgm]

${p}_{2}=\upsilon $

600 [kg/(mm^{3})]

1.2×10^{7} [s/(kg(mm)^{5})]

0.6

3.757×10^{3} [kgm]

7.284×10^{3} [kgm]

${p}_{3}=\u0430$

1.0 [s/m]

2.1×10^{12} [s^{2}/((kg)^{2} (mm)^{3})]

0.5

2.84×10^{3} [kgm]

6.156×10^{3} [kgm]

${p}_{4}=k$

0.5

1.1×10^{10} [s/((kg)^{2} (mm)^{2})

0.1

2.84×10^{3} [kgm]

6.156×10^{3} [kgm]

${p}_{5}=$
${T}_{1}$

1.1×10^{2}

1.2×10^{11} [s^{2}/((kg)^{2} (mm)^{2})]

1.5

4.46×10^{3} [kgm]

9.66×10^{3} [kgm]

${p}_{6}={T}_{2}$

1.5×10^{2}

1.7×10^{11} [c^{2}/((kg)^{2} (mm)^{2})]

1.6

3.46×10^{3} [kgm]

9.66×10^{3} [kgm]

Fig. 4 shows examples of trajectories of elastic deformational displacements in plane ${X}_{1}^{}$, ${X}_{2}^{}$ in variations in relation to equilibrium points (trajectories on the left) as well as projections of phase trajectories in plane ${X}_{1}^{}$, $d{X}_{1}^{}/dt$ (trajectories in the middle) and ${X}_{2}^{}$, $d{X}_{2}^{}/dt$ (trajectories on the right). The bifurcation diagram (Fig. 2) shows values of work (black circles) for which trajectories are calculated in Fig. 4. The bifurcation diagram and examples of trajectory changes in the space of elastic deformational displacements demonstrate a large variety of formed attractors that change as friction forces do work. At the same time, as the work is done, certain regions stand out where the topology of the phase space remains unchanged, along with bifurcation points where topology changes. It should be emphasis that the transformation of topology of the phase space of the deformational displacements as a consequence of changes in the interactions of contacted surfaces, changes such properties of the contact as the power of irreversible transformations of the energy introduced into the contact. Therefore, for example, such parameters as wear intensity (increment of wear along the path), transformation of surface roughness, including the formation of equilibrium roughness, etc., do not remain unchanged after bifurcation points in the process of system evolution. These properties affecting the functional suitability of friction pairs during operation can change fundamentally depending on the parameters of the interacting subsystems.
Fig. 4. Examples of trajectories in phase space sections depending on the work done by friction forces
Presented in the article materials allow explaining the experimental investigations carried out during the last thirty years in the direction of evolution of the dynamic friction system. They are described in detail in monographs [10], [14], [36], [38] and articles [15], [22], [29], [30]. Let us briefly consider the results of the experimental study of the vibration characteristics of contacting pairs during the development of their wear. For concreteness, two systems: steelbrass and steelglycerinbrass are considered. The vibrations in the contact interaction force space and the power of irreversible transformations as wear develops were studied. Vibrations were measured with vibration accelerometers, with subsequent integration and elimination of the trend due to uncertainty in initial conditions. The obtained vibration sequences in the form of vibration displacements, velocities and accelerations made it possible to reconstruct the vibration characteristics, their frequency images in the form of auto and mutual spectra as well as in the form of evolution diagrams of root distribution of the characteristic polynomial of autoregressive models of dynamics. Traditional statistical functional analyses of time sequences, as well as sequences based on the Poincare stroboscopic mapping algorithms, have been used in the analysis of oscillations. The following properties were found:
a) as wear progresses, there is a restructuring of the vibrational displacement attractors and associated contact interaction forces, including bifurcations of the attraction sets of deformations;
b) as a rule, there is a critical wear after which chaotic oscillations develop in the system. It is established that the general tendency of evolution of tribosystem from determinism at the initial stage to chaos occurs after some critical work of friction forces;
c) evolution and moments of bifurcations depend on initial parameters of the system, for example, on elements of a matrix of rigidity of indenter suspensions;
d) the intensity of wear depends on the evolution of dynamic contact properties and it can change fundamentally after bifurcation points;
e) in the same external friction conditions, but when changing, for example, the elements of the indenter subsystem stiffness matrix, selective transfer in the steelglycerinbrass system may or may not be formed. Experiments show great sensitivity of selective transfer formation to small variations in friction system parameters as well as to small perturbations such as rotating spindle runout.
We can continue the list of evidence showing the dependence of the output properties of the friction unit on the dynamic characteristics of the friction unit, and it is well known that these properties change in the course of functioning of the friction unit. All previously considered mathematical models of friction process dynamics, were “frozen”. They considered individual aspects of the dynamics, for example, the formation of attracting sets of deformations of the type of limit cycles [14], [15]. They did not allow us to analyze the evolution of the dynamical system, to consider bifurcations of the attracting sets of strain displacements formed during friction as the work is performed by friction forces. Finally, the friction process was considered without emphasis on its systematic character, including the internal connection of the power of irreversible energy transformations, the intensity of wear of contacting surfaces, the transformation of their geometrical topology, etc. But it is the consideration of the system, not the process, the functional connectivity of the system coordinates, their evolution, that are the main distinctive features of the proposed methodology and mathematical modeling, which analytically explain the above mentioned experimentally observed properties of friction units.
Our studies aimed at forming wearfree friction conditions have shown that the selective transfer may not be formed in friction of steelbronze in glycerol in depend on the dynamic parameters of the indenter [37], [38]. There is other, as a rule, experimental studies aimed at revealing the connection of dynamics and wear resistance of friction units [15][19]. In this case at studying of output properties of friction, for example, wear resistance of contacting pairs, there is a fundamental problem of sensitivity of output properties of friction process to variations of dynamic parameters of interacting subsystems. Therefore, by analogy with dynamic systems, it is possible to introduce the concept of robustness on the output properties of the friction system.
The system approaches, the disclosure of functional relationships between coordinates provide an analytical and semantic basis for dynamic monitoring of the state of friction units [10], [21], [3234]. The developed mathematical toolkit and researches of evolution of attractive sets of deformation displacements allow defining reliable information models of diagnostics of knots of a friction knowingly on the basis of mathematical modeling. For example, our analytical and experimental studies have shown that chaotic dynamics always develops in the dynamic system when tool wear approaches its critical value. The dynamic monitoring systems for rotary tribosystems developed on this basis have been implemented at a number of enterprises in Russia. Thus, the developed technique of mathematical modeling of evolution of friction system is a way of estimation of evolutionary transformations of dynamic system of a friction depending on naturally formed phase trajectory of power of irreversible transformations in the friction knot on the produced work.
7. Conclusions
The performed studies allow us to draw a general conclusion about the systemsynergetic properties of the dynamic friction system. The representation of the forces formed in the interface area of the contacting surfaces in the state coordinates leads to the fact that these forces form an internal feedback and begin to play the role of an internal state regulator. Since the force model is nonlinear, such effects as branching of equilibrium points are possible in the system depending on the parameters. It is shown that in the vicinity of these points it is possible to form all the attractive sets known in nonlinear dynamics (limit cycles, invariant tori, chaotic attractors). Since the parameters of the system depend on the power of irreversible energy transformations in the friction zone, during the operation of the friction process, its dynamic restructuring is observed, which changes the output characteristics of friction and wear. It should be emphasized that the naturally released power of input energy in the contact area actually characterizes the factor not considered earlier, which changes the properties of the friction system. It is considered not “frozen”, but evolutionary changing, and the "generator" of these changes is energy release in the friction zone. Depending on the parameters of interacting subsystems, evolutionary trajectories, including bifurcations of attractive sets of strain displacements and their corresponding forces, change the properties of the friction system. The accounting of these changes, the methodology of modeling and studying as a single system fundamentally distinguish the given results from the known ones in the field of friction dynamics. Let us list the following most important provisions that follow from the completed study.
The dynamic friction system is formed as a result of interaction between two mechanical subsystems through a dynamic link created in the zone of the contact of interacting surfaces. The dynamic link that represents a model of forces in the state coordinates of the system is nonlinear in relation to elastic deformational displacements and their speeds. Therefore, the dynamic system of friction can serve as an example of various effects in nonlinear dynamics.
While examining linearized equations in variations in relation to the equilibrium point of elastic deformational displacements, matrices of speed coefficients and elasticity become nonsymmetrical due to the reaction by the process of friction. Therefore, dissipative, acceleration, gyroscopic, potential and circulatory forces are created in the course of the friction process. The main mechanisms of equilibrium stability loss in the system are associated with the transformation of the symmetric part of the speed coefficient matrix from positive definite to negative definite as well as with a significant increase in circulatory forces in relation to their potential components. The formation of circulatory forces makes it possible to explain the always experimentally observed trajectories of periodic deformation displacements of the contacting pairs in the form of circular precessional oscillations.
The principal property of a dynamic friction system is the change in the parameters and properties of the dynamic coupling formed, depending on the phase trajectory of the power of irreversible energy transformations in the contact region according to the work performed. As a result, the dynamic friction system undergoes evolutionary changes in its properties during its operation. They manifest themselves in the branching of equilibrium points (forktype bifurcations) and in the change in the topology of the phase space in the regions surrounding the equilibrium. These changes, depending on the bifurcation points of the deformation displacements and the corresponding forces and powers, cause changes in the physical interactions on the contact surfaces. Moreover, the areas of stationary states alternating with changes in the topology of the phase space lead to the formation of chaotic dynamics in the deformation displacements. Chaotic dynamics usually correspond to the stage of catastrophic wear and tear. Therefore, this opens possibilities for building systems to diagnose the tribosphere state from its vibration data, for example, using models of autoregressive spectral analysis. This method has been tested in diagnostics of bearing units in machines, as well as in machining of materials by cutting. In the latter case, due to the coordination of the trajectories of the machine's executive elements programmed by the CNC system with the evolutionary changes of the system, the possibility of reducing the reduced costs of manufacturing a batch of parts by 1.52.0 times is shown.
It is found that the evolutionary changes in the properties of the dynamic friction system are sensitive to the variations in the parameters and the perturbations that always occur. Therefore, firstly, the wear of the contacting pairs changes with variations in the system parameters, and secondly, it does not remain unchanged at individual stages of evolution. Thus, the results of studying the wear resistance of the contacting pairs obtained on experimental stands can be transferred to a specific machine with great care. Moreover, we can introduce the concept of “roughness” of the system to the evolutionary changes in its properties.
The developed mathematical tools are based on the synergetic paradigm of analysis and synthesis, according to which the mathematical modeling of interacting subsystems is based on the representation of interactions in state coordinates. In this case, the interacting subsystems are transformed into a single interconnected dynamic system. Changes in the properties of the bond, for example, as a result of the release of energy to the contact area, cause a change in the entire system. The given material is limited by mechanical interactions represented by forces, and the coordinates of the state are considered in the form of spatial elastic deformation displacements. It should be taken into account that the release of energy in the contact zone causes many other physical interactions, for example, the power of irreversible transformations is a source of heat production. Then, as the coordinates of the state, it is necessary to consider the thermodynamic parameters and they cannot but affect the mechanical parameters. Thus, the further development of the ideas given in the article should be extended to the system physics of interactions: thermodynamics; tribochemistry, including the problems of selective transport; molecular interactions (adhesion, diffusion), etc.
Acknowledgements
This study was performed with financial support of RFBR grants: “Development of the theory of analysis and synthesis of controlled selforganization in a dynamic cutting system during processing on the example of manufacturing parts on lathes”, No. 190800022; “Development of methods for analysis and synthesis of a dynamic cutting system based on the criterion of minimizing the reduced manufacturing costs”, No. 203890074.
References
 W. R. Ashby, “Principles of the selforganizing dynamic system,” The Journal of General Psychology, Vol. 37, No. 2, pp. 125–128, Oct. 1947, https://doi.org/10.1080/00221309.1947.9918144 [Publisher]
 I. Prigogine, Introduction to Thermodynamics of Irreversible Processes. New York: Interscience, 1961. [Search CrossRef]
 H. Haken, Springer Series in Synergetics. Wiesbaden: Vieweg+Teubner Verlag, 1973, pp. 9–19, https://doi.org/10.1007/9783663015116 [Publisher]
 L. I. Bershadski, “On the selforganization and concepts of wearresistance in tribosystem,” (in Russian), Friction and Wear, Vol. 13, pp. 1077–1094, 1992. [Search CrossRef]
 M. D. Bryant, “Entropy and dissipative processes of friction and wear,” FME Transaction, Vol. 37, No. 2, pp. 55–60, 2009. [Search CrossRef]
 M. Nosonovsky and V. Mortazavi, FrictionInduced Vibrations and SelfOrganization: Mechanics and NonEquilibrium Thermodynamics of Sliding Contact. Florida: CRC Press, 2013. [Publisher]
 A. S. Kuzharov, R. Marchak, J. Guzik, K. Kravchik, and E. G. Zadoshchenko, “Research of the phenomena of the tribological selforganization in the "brassglycerinsteel" system,” (in Russian), Friction and Wear, Vol. 17, No. 1, pp. 113–122, 1996. [Search CrossRef]
 I. Gershman, P. Peretyagin, A. Mironov, and E. Gershman, “Description of tribosystem processes and development of wearresistant materials using nonequilibrium thermodynamics and the theory of selforganization,” (in Russian), VNIIZHT Scientific Journal, Vol. 75, No. 5, pp. 271–277, 2016. [Publisher]
 Y. K. Mashkov, “Nanostructural selforganization and dynamic adaptation of metal–polymer tribosystems,” Technical Physics, Vol. 62, No. 2, pp. 282–286, Feb. 2017, https://doi.org/10.1134/s1063784217020190 [Publisher]
 V. L. Zakovorotny, Dynamics of Tribosystems SelfOrganization. (in Russian), RostovonDon: DSTU, 2003. [Search CrossRef]
 A. A. Kolesnikov, Synergetic and Problems of Control Theory. (in Russian), Moscow: Fizmatlit, 2004. [Search CrossRef]
 O. BenDavid, S. M. Rubinstein, and J. Fineberg, “Slipstick and the evolution of frictional strength,” Nature, Vol. 463, No. 7277, pp. 76–79, Jan. 2010, https://doi.org/10.1038/nature08676 [Publisher]
 Ea. Kuznetsov, “Effect of fluid lubricant on the contact characteristics of rough elastic bodies in compression,” Friction and Wear, Vol. 157, pp. 177–194, 1985. [Publisher]
 I. V. Kragelsky and N. V. Gittis, Frictional Autoscillations. (in Russian), Moscow: Nauka, 1987. [Search CrossRef]
 M. Wiercigroch and A. M. Krivtsov, “Frictional chatter in orthogonal metal cutting,” Philosophical Transactions of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences, Vol. 359, No. 1781, pp. 713–738, Apr. 2001, https://doi.org/10.1098/rsta.2000.0752 [Publisher]
 O. Giannini, A. Akay, and F. Massi, “Experimental analysis of brake squeal noise on a laboratory brake setup,” Journal of Sound and Vibration, Vol. 292, No. 12, pp. 1–20, Apr. 2006, https://doi.org/10.1016/j.jsv.2005.05.032 [Publisher]
 A. Ammar, M. Yousif, and Ab. R. Inzarulfaisham, “Analyzing the disc brake squeal: review and summary,” International Journal of Scientific and Technology Research, Vol. 2, No. 4, pp. 60–72, 2013. [Search CrossRef]
 M. Nosonovsky and G. G. Adams, “Vibration and stability of frictional sliding of two elastic bodies with a wavy contact interface,” Journal of Applied Mechanics, Vol. 71, No. 2, p. 154, 2004, https://doi.org/10.1115/1.1653684 [Publisher]
 J. C. Sprott, “A dynamical system with a strange attractor and invariant tori,” Physics Letters A, Vol. 378, No. 20, pp. 1361–1363, Apr. 2014, https://doi.org/10.1016/j.physleta.2014.03.028 [Publisher]
 A. A. Alifov, “Autoparametric oscillations with delays in elastic and frictional forces,” Journal of Machinery Manufacture and Reliability, Vol. 50, No. 2, pp. 98–104, Mar. 2021, https://doi.org/10.3103/s1052618821020023 [Publisher]
 V. L. Zakovorotny, V. P. Blokhin, and M. I. Alekseichik, Introduction to Dynamics of Tribosystems. (in Russian), RostovonDon: InfoService, 2004. [Search CrossRef]
 R. Rusinek, M. Wiercigroch, and P. Wahi, “Influence of tool flank forces on complex dynamics of cutting process,” International Journal of Bifurcation and Chaos, Vol. 24, No. 9, p. 1450115, Sep. 2014, https://doi.org/10.1142/s0218127414501156 [Publisher]
 R. Rusinek, M. Wiercigroch, and P. Wahi, “Modelling of frictional chatter in metal cutting,” International Journal of Mechanical Sciences, Vol. 89, pp. 167–176, Dec. 2014, https://doi.org/10.1016/j.ijmecsci.2014.08.020 [Publisher]
 V. L. Zakovorotnyi and V. S. Bykador, “Cuttingsystem dynamics,” Russian Engineering Research, Vol. 36, No. 7, pp. 591–598, Jul. 2016, https://doi.org/10.3103/s1068798x16070182 [Publisher]
 A. Masoumi, F. Pellicano, F. S. Samani, and M. Barbieri, “Symmetry breaking and chaosinduced imbalance in planetary gears,” Nonlinear Dynamics, Vol. 80, No. 12, pp. 561–582, Apr. 2015, https://doi.org/10.1007/s1107101418903 [Publisher]
 A. M. Gouskov, M. Guskov, P. Lorong, and G. Panovko, “Influence of flank face on the condition of chatter selfexcitation during turning,” International Journal of Machining and Machinability of Materials, Vol. 19, No. 1, p. 17, 2017, https://doi.org/10.1504/ijmmm.2017.081186 [Publisher]
 V. L. Zakovorotny and V. E. Gvindjiliya, “Bifurcations of attracting sets of cutting tool deformation displacements at the evolution of treatment process properties,” Izvestiya VUZ. Applied Nonlinear Dynamics, Vol. 26, No. 5, pp. 20–38, 2018, https://doi.org/10.18500/0869663220182652038 [Publisher]
 V. L. Zakovorotnyi, A. A. Gubanova, and A. D. Luk’Yanov, “Parametric selfexcitation of a dynamic endmilling machine,” Russian Engineering Research, Vol. 36, No. 12, pp. 1033–1039, Dec. 2016, https://doi.org/10.3103/s1068798x16120194 [Publisher]
 V. L. Zakovorotny and V. E. Gvindjiliya, “Link between the selforganization of dynamic cutting system and tool wear,” (in Russian), Applied Problems of Nonlinear Oscillation and Wave Theory, Vol. 28, No. 1, pp. 46–61, Feb. 2020, https://doi.org/10.18500/0869663220202814661 [Publisher]
 V. Zakovorotny and V. Gvindjiliya, “Influence of cutting dynamic on the selection of the technological regimes to ensure minimal wear of cutting tools,” (in Russian), Metal Working and Material Science, Vol. 22, No. 4, pp. 54–70, Dec. 2020, https://doi.org/10.17212/19946309202022.45470 [Publisher]
 Yu. I. Osenin, Yu. V. Krivosheya, A. V. Chesnokov, and V. K. Antoshkin, “Influence of mutual overlap coefficient on the process of the emergence of the «creak» during the braking of rolling stock by disk brakes,” Friction and Wear, Vol. 42, No. 1, pp. 63–71, 2021, https://doi.org/10.32864/0202497720214216371 [Publisher]
 A. Akay, “Acoustics of friction,” The Journal of the Acoustical Society of America, Vol. 111, No. 4, pp. 1525–1548, Apr. 2002, https://doi.org/10.1121/1.1456514 [Publisher]
 S. I. Builo, V. D. Vereskun, V. I. Kolesnikov, D. S. Manturov, and O. N. Popov, “Determining friction coefficient at runin stage and diagnosing the point of transition to steadystate phase based on acoustic emission signals,” Russian Journal of Nondestructive Testing, Vol. 56, No. 1, pp. 41–48, Jan. 2020, https://doi.org/10.1134/s1061830920010039 [Publisher]
 V. I. Kolesnikov, P. A. Koropets, and E. S. Sinyutin, “Vibrationbased diagnostics of helicopter transmission gear couplings,” Russian Aeronautics, Vol. 63, No. 1, pp. 7–13, Jan. 2020, https://doi.org/10.3103/s106879982001002x [Publisher]
 I. V. Kragelsky, Friction and wear. (in Russian), Moscow: Engineering, 1968. [Search CrossRef]
 V. L. Zakovorotny, M. Marchak, I. V. Usikov, and A. D. Lukyanov, “Interrelation between tribosystem evolution and parameters of dynamic friction system,” Journal of Friction and Wear, Vol. 19, No. 6, pp. 54–64, 1998. [Search CrossRef]
 V. L. Zakovorotny, V. E. Gvindzhiliya, and P. S. Kolodkin, “Effect of dynamic properties of interacting subsystems on evolution of selective transfer formation in friction units,” (in Russian), Vestnik of Don State Technical University, Vol. 19, No. 2, pp. 104–112, Jun. 2019, https://doi.org/10.23947/199259802019192104112 [Publisher]
 V. L. Zakovorotny, NonLinear Tribomechanics. (in Russian), RostovonDon: DSTU, 2000. [Search CrossRef]
 V. L. Zakovorotny, V. A. Shapovalov, and V. V. Shapovalov, “The study of complex friction coefficient,” (in Russian), Friction and Wear, Vol. 8, No. 5, pp. 722–728, 1987. [Search CrossRef]
 A. M. Lyapunov, The General Problem of the Stability of Motion. (in Russian), London: Taylor & Francis, 1992. [Search CrossRef]
 D. R. Merkin, Introduction in the Theory of Motion Stability. (in Russian), Moscow: Nauka, 1971. [Search CrossRef]
 V. M. Likhadanov, “On the influence of force structure on stability of motion,” (in Russian), Applied Mathematics and Mechanics, Vol. 38, pp. 246–253, 1974. [Publisher]
 V. M. Likhadanov, “On stabilizing potential systems,” (in Russian), Applied Mathematics and Mechanics, Vol. 39, pp. 53–58, 1975. [Publisher]