CC BY
101
УДК 531.8, 531.44, 528.5-531.7
Диссипация энергии трения в контакте упругих тел под воздействием наложенных нормальных и касательных колебаний
T. Hanisch, I. Richter, Q. Li
Берлинский технический университет, Берлин, 10623, Германия
В статье проведено численное моделирование фрикционного демпфирования в упругом контакте параболического индентора под воздействием комбинации колебаний в нормальном и касательном направлении. Показано, что рассеиваемая энергия сначала линейно растет с увеличением коэффициента трения, затем линейно уменьшается и затем выходит на постоянное значение. Данные участки зависимости соответствуют режимам полного скольжения, частичного скольжения и полного сцепления. Все три асимптотические зависимости могут быть описаны аналитически. Рассеиваемая энергия в безразмерной форме является функцией отношения амплитуды нормальных колебаний к средней глубине вдавливания индентора, отношения изменения площади контакта к площади сцепления, а также величины фазового сдвига между нормальным и тангенциальным колебанием. Предложены обобщающие кривые.
Ключевые слова: трение, диссипация энергии, прерывистое скольжение, двойное колебание, упругий контакт, механика контакта
DOI 10.24411/1683-805X-2020-12007
Frictional energy dissipation in a contact of elastic bodies subjected to superimposed normal and tangential oscillations
T. Hanisch, I. Richter, and Q. Li
Technische Universität Berlin, Berlin, 10623, Germany
Frictional damping in elastic contact of a parabolic indenter subjected to a combination of oscillations in normal and tangential directions is numerically simulated. The dissipated energy first increases linearly with coefficient of friction, then decreases linearly, and finally reaches a constant value. These three regions correspond to the states of complete slip, partial slip and complete stick. All three asymptotical dependencies can be described analytically. The dissipated energy in a dimensionless form is function of the ratio of normal oscillation amplitude and mean indentation depth, the ratio of change in contact area and sticking area, and phase shift between normal and tangential oscillation. Master curves are suggested.
Keywords: friction, energy dissipation, stick-slip, dual oscillation, elastic contact, contact mechanics
1. Introduction
In frictional contacts, energy is dissipated when two contacting bodies have a relative sliding movement. For elastic bodies, it is well known that under the periodic oscillating loading in tangential direction, the microslip appears at the boundary of contact and stick in the middle of contact, which may lead to fretting and initiation of fatigue cracks [1]. This frictional damping occurs very common in the interface of
joints of machine components [2, 3] and plays an important role in many applications of tribology and structure mechanics [4]. The energy dissipation of a spherical indenter subject to tangential oscillation was analyzed early by Mindlin [5]. It was found that the dissipated energy in one cycle is inversely proportional to the coefficient of friction, which indicates that there will be no energy dissipation if the coefficient of friction is infinitely large, because the whole contact is in a state of stick. However, a recent study shows
© Hanisch T., Richter I., Li Q., 2020
that even in the case of infinitely large coefficient, energy dissipation still occurs if the body oscillates in both vertical and tangential directions, because the elastic energy stored at the boundary elements of contact is suddenly relaxed during the composed oscillating process [6]. This kind of energy dissipation is called "relaxation damping". In this paper, we numerically study the frictional damping due to a combination of vertical and tangential oscillation with constant coefficient of friction in contact under the Coulomb's law of friction.
The energy dissipation under varying normal and tangential loading has been studied by many researchers, for example analytically by Davies et al. [7] for smooth two-dimensional indenters and by Putignano et al. [8] for rough surfaces, numerically by Liu and Eriten for two-dimensional wavy surfaces using the finite element method [9], and experimentally early by Johnson [10], Goodman and Brown [11] and recently by Usta et al. [12], where a power-law relation between dissipated energy and maximal applied stress has been intensively discussed and the power-law exponent is argued between 2 and 3. Furthermore, studies have shown that phase shift between normal and tangential oscillation plays an essential role in frictio-nal energy dissipation [13, 14], and the maximal energy dissipation occurs in many cases when phase different is tc/2 [8, 14]. In this paper, we carry out simulations of contact due to a combination of normal and tangential oscillation using the method of dimensionality reduction [15-17]. This is a very effective analytical and numerical tool exactly for this type of contact problems where only the total macroscopic force and displacement are of importance. Both these quantities are determined in the framework of method of dimensionality reduction exactly, provided Coulomb's law of friction is assumed.
2. Mathematical model
We consider a contact between a rigid parabolic in-denter with profile f (r) = r2/(2R) and an elastic half space with elastic modulus E and Poisson' ratio v, where R is radius of indenter and r polar radius in the contact plane. After indentation by d0, indenter is forced to oscillate in vertical and tangential directions with angular frequency ra, phase difference 9 and small amplitudes Auz and Aux according to the following displacement-controlled laws:
Fig. 1. Equivalent contact between a parabolic indenter and an elastic half space in the framework of the method of dimensionality reduction
We consider the case of "no jumping" and small amplitude, therefore Auz << d0.
In the framework of the method of dimensionality reduction, the three-dimensional profile is transformed into a plane profile according to [15, 16]
|x| f'(r) g(x) = |x|f / ( ) dr.
(2)
u(0)(t) = d0 + Au(0)cos(rat), uX0)(t) = AuX0)cos(rat + 9).
(1)
I~2 2
0 v x - r
For the parabolic indenter, its corresponding one-dimensional profile is given by g(x) = x2/R. Furthermore, the elastic half space is replaced by a one-dimensional elastic foundation consisting of an array of independent springs with discrete distance Ax (Fig. 1). The normal and tangential stiffness of springs, Akz and Akx are defined following the rules:
Akz = E*Ax, Akx = G*Ax, (3)
where effective elastic modulus is E* = E/(1 - v2) and shear modulus G = 4G/(2 - v). With profile and elastic foundation defined in Eqs. (2) and (3), one can simply solve the normal or tangential contact problems. The normal and tangential forces on each spring in contact is easily calculated by its stiffness Akz, Akx and displacement uz(x, t) and ux(x, t)
Afz(x, t) = Akzuz (x, t), Afx(x, t) = Akxux(x, t). (4)
The normal displacement is dependent of only the profile g(x) and the given normal oscillation uz(0)(t)
uz (x, t) = u(0)(t) - g(x). (5)
The tangential displacement can be determined by the following Coulomb's law of friction: firstly, we assume that all springs in contact are in stick state and have the same incremental tangential displacement as indenter. If the resulted tangential force on some spring is larger than the production of coefficient of friction ^ and normal force Afz(x, t), then it is in a state of slip, and the tangential force should be corrected according to Coulomb's law. So, for a given small
incremental displacement of indenter dux(0)(t), we have the following rules
dux (x, t) = duf(t), ux (x, t) = ±^Af(xt). (6)
Akx
The details on application of the method of dimensionality reduction to normal and tangential contact can be found in paper [15]. With obtained tangential displacement of springs one can then calculate the force according to Eq. (4) as well as the total tangential force by summing the spring forces
Fx(t) = Z Afx(x, t).
cont
The energy in one period of oscillation T = 2^/ra is then given as
W = J Fx (t )du^(t ).
o
(7)
3. Results
Wœ = 8 G'R^V2 AuZ0)(Auf)2si
sin 9.
(9)
The maximal damping Wœ,max occurs when the phase shift is 9 = rn/2:
Wœ,max = 3 ^R^AuftAuf)2. (10)
Considering another limiting case of very small coefficient of friction where the whole contact area is in a state of slip, then the tangential force is simply following the Coulomb's law of friction: Fx = pFn, and the dissipated energy in one oscillation cycle is equal to
Wc =MFn4Auf. (11)
Substitution of solution of normal load in Hertzian contact, Fn = 4/3E*R1/2d03/2 (we assume Auz << d0, so
the influence of vertical amplitude is neglected) into Eq. (11) provides
WC = 16/3pE*R1/2d03/2AuX0). (12)
Equations (8), (9), (12) are three analytical solutions which will be used in the following analysis.
Now we discuss one important parameter in description of energy dissipation: the ratio of the change in contact area and sticking area Aa/Ac. For the normal contact, it is well known that the contact radius is geometrically related to indentation depth d and sphere radius a = so derivative of contact ra-
dius with respect to indentation depth results in the change in contact radius
Aa =
1 R
Ad.
(13)
3.1. Theoretical analysis
Before presenting numerical results, we introduce the existing important results from the literature and offer a brief discussion on the current study. In the case of only tangential oscillation with finite coefficient of friction, the dissipated energy in one period of oscillation was given by Mindlin [5, 15]
WM = I K-WV^Auf)3 p-1, (8)
where k = E*/G*. It is inversely linear function of coefficient of friction, thus there will be no energy dissipation if coefficient of friction is infinitely large p = <x>. However, in the case of composition of vertical and tangential oscillation according to Eq. (1), the relaxation damping with p = <x> appears and it was given by Popov et al. [6]
In an oscillating contact, the maximum change in indentation depth is Ad = 2Auz(0), then Eq. (13) becomes
Aa = R12d-1/2AuZ0). (14)
For tangential contact, the contact radius of sticking area c is determined by the relation (the normal amplitude of oscillation is still neglected compared with mean indentation depth)
Guf = pE (do -g (c)) = pE
2 A
R
(15)
Similar to normal oscillation, considering the maximal change in tangential displacement (absolute value), it has,
(16)
2 AuX0)j = jpKR-1 • 2cAc I
For small amplitude of oscillation, the radius of slip area is also very small, then we have c « a = 4Rd in comparison with Ac. Following that Eq. (16) has the form
(17)
Ac = Au(.o;R
(0)r>1/2,
-1/2p-1K-1.
From Eqs. (14) and (17) the ratio of change in contact area and sticking area Aa/Ac is then equal to
- Aa Au(0)
-pK, (18)
Ac
Au?->
which is denoted by p.
In this study, we consider dual oscillation but with finite coefficient of friction, the numerically obtained dissipated energy will be normalized by the maximal value Wœ,max in the limiting case of p = <x> in Eq. (9) with 9 = tc/2. Then the normalized solution by Mind-lin in Eq. (8) is
WM
M
W
Au
(0)
Au
(0)
pK
,or Wm = 4. (19) 4P
Fig. 2. Map of slip, stick, and noncontact area in two cycles of oscillation for parameter set: 9 = rc/2, p = 0.7268 and 8 = 0.01 (color online)
The solution for the case of complete slip (12) in the normalized form is
WC
(
W
• = 2
Y
Au
(0)
Au
(0)
Au
(0)
pK, or
WC = 2ô"2p
(20)
with ratio of normal oscillation amplitude and mean indentation depth
- Au (0) 8=^. (21)
d0
Normalized solution for relaxation damping (9) is then phase dependent
W 2
™ - = sin 9. (22)
W =-
W
®,max
From Eqs. (19), (20) and (22), we can see that the normalized dissipated energy in case of Mindlin depends only on p, in the case of relaxation damping by Popov et al. only on phase 9, and in the case of Coulomb on both 5 and p. In the following part, numerical results show that the dimensionless energy dissipation in one oscillation period in a general case is a function of these three parameters:
W ~ f (5, p, 9). (23)
3.2. Numerical results
The frictional contact was numerically simulated using the method of dimensionality reduction as described in Sect. 2. Figure 2 shows an example of contact area, sticking area and area out of contact changing with time in two cycles of oscillation 2T for parameter 9 = tc/2, p = 0.7268 and 5 = 0.01. Vertical axis shows the coordinate in plane ranging from the minimal boundary of stick-slip area cmin to the maximal contact radius amax. One can see that the contact
radius varies with a harmonic-like function which could be also analytically calculated according to a(t) = (Ruz(0)(t))1/2. Focus on only one time moment, for example at time t\ (dashed line in Fig. 2), there is a stick region in the middle (upper part in yellow), slip region at the contact boundary (middle part in orange), and noncontact region (lower part in blue). But from the map it is seen that slip does not exist all the time, for example at time t3 the whole contact is in a state of sticking. Thus, energy dissipation occurs not all the time, but only in the time intervals when slip appears. Interestingly, one can see that in dual oscillation, the slip region appears and spreads gradually, but vanishes suddenly to a state of complete stick (for example at time t2).
A phenomenon should be noted here: in Fig. 2 the first cycle behaves slightly differently than the following one. The reason for that is the initiation of spring locations, therefore only the second cycle is considered below for the calculation of the dissipated energy per cycle.
Figure 3, a shows the dependence of the normalized dissipated energy W on the parameter p for phase 9 = tc/2 and three different values of 5 = 10-4, 10-3 and 10-2. The parameter p was changed by varying the coefficient of friction p. Focus on one single curve, one can see that the dissipated energy increases with coefficient of friction as well as parameter p, then it decreases until reaches to a constant value.
This dependence can be divided into three regions:
- a linear dependence in region I where the complete slip occurs according to Coulomb's law of friction described by Eq. (20),
- an inversely proportional dependence in region II which can be described by the Mindlin's solution (Eq. (19)) for a state of partial sliding, and
- constant value in the region III corresponding to the pure "relaxation damping" for a state of complete sticking described by Popov et al. in Eq. (22) [6].
In this normalized form, three curves overlap at large value of p in .regions II and III where the energy is independent of 5, and the normalized dissipated energy is equal to 1 in the plateau with this example 9 = tc/2: W = 1. A multiplication of the normalized dissipated energy by the ratio 5 and its reciprocal by p, as shown in Fig. 3, b, leads to an overlap of the curves at small ratios of p/5 in regions I and II. These behaviors can be described by Eqs. (19)-(22).
In Fig. 4, a the dependences of dissipated energy W on parameter p for 5 = 10-2 and different phases 9 are shown. A master curve is generated at small values of p in regions I and II, so the dissipated energy
W 103 10210110° 10-1
i-i
10
10"
10"
i-3
10
10"
1-4
cp = 7l/2
— 8= 10"4
✓ \\ / -- ô= 10"3
/ \ ' N... / tA\ / \\ -- s = 10~2
w
104 10 2 10° 102 i
V Region I / \ cp = tz/2
/ \ Region II
/ — § = 10"4 V'" -- 6= 10 3 \ Region III
----- 8 = 10~2 V"
— Coulomb \
Mindlin
— Popov et al. [6]
10
10
10
104
p/ô
Fig. 3. Dependence of dissipated energy on parameter p for different amplitudes of normal oscillation 5 from 10-4 to 10-2: the curves overlap at large values of p (a); the curves overlap at small values of p/5 (b) if the coordinates are multiplied by 5 and its reciprocal (color online)
is independent of phase in this range. With increasing coefficient of friction, they are dispersed because the dissipated energy with very large coefficient of friction in the case of relaxation damping is phase dependent. If the energy is normalized as W/Wœ(9) by taking into account the phase angle, then the curves tends towards W/Wœ(9) = 1 in the region of plateau (Fig. 4, b).
To find an "empirical" equitation describing all three regions, two options are presented below.
The first possibility is to describe three regions separately with already known theories, as discussed above:
P
- 1 0
W = 2=^, region I, for 0 < p< —1=5, 52 v8
W = region II, for 5<p<—(24) 4p V8 4sin2 9
W = sin 9, region III, for p> -
1
4sin 9
The boundary between two regions are simply obtained by equilibrium of two equations. These three relations are shown in Fig. 3, b with solid lines. In Figs. 3
Fig. 4. Dependence of dissipated energy on parameter p for different phases for 5 = 10-2. The energy is normalized by Wœmax (a) and by Wœ(9) (b) in its corresponding phase 9
and 4 numerical results show that there is a smooth transition between regions I and II, as well as between regions II and III. However, Eq. (24) does not cover these transition zones.
Another possibility is to approximate the simulation results using rational functions. For region I, beginning of region II and their transition, a perfect master curve is observed as seen in Figs. 3, b and 4, a, the results are dependent of 8 but independent of phase 9. Taking 8 into account, we have the following rational function for 0 < £ < 4:
W8 = (4.858 • 10-4£5 + 0.2203^4 - 0.2226£3 - 0.1252£2 + 0.1836£1)(£5 -1.932£4 +1.482£3
0.58662,2 + 0.0779£1 + 0.0889)
-1
(25)
with definition of parameter
£1 =|. (26)
8
This approximation is shown in Fig. 5 on the left side with the example of 8 = 10-2, which agrees with numerical results very well. It is noted that the range of £1 in Eq. (25) is evaluated based on the curves in Fig. 3, b. For other cases, for example 8 > 10-2, the range will be reduced. The importance of (25) is the
Fig. 5. Approximation to the numerical result using rational functions
description of transition zone, so for the linear part, more exact solution (24) is suggested.
For the other transition between regions II and III, one can see that shape of the curves in this area are different for different phases 9 (Fig. 4), thus a master curve cannot be generated, therefore we give here only an approximation for a special case of 9 = rc/2 for p > 10-1:
W = (p4 + 0.7169p3 + 0.2668p2 + 0.3342p - 40017-10-6)(p4 + 0.5979p3 + 1.335p2 - 5.079-10-5p)-1.
(27)
4. Conclusion
The frictional contact of a parabolic indenter and an elastic half space, while the indenter is subjected to oscillations in normal and tangential directions, is numerically simulated using the method of dimensionality reduction. The dissipated energy in one oscillation cycle is studied for different coefficients of friction, oscillation amplitudes and phase shifts between vertical and horizontal direction. It is found that the dissipated energy increases linearly with coefficient of friction (region I), then decreases linearly (region II), finally reaches to constant (region III). These three regions correspond to states of complete slip, partial slip and complete stick. The can be described by the known asymptotic solutions based on the Coulomb's law of friction, Mindlin's solution and solution for relaxation damping.
The dissipated energy in dimensionless form occurs to be function of only three dimensionless parameters: normal oscillation amplitude, ratio of change in contact radius and radius of sticking area, and phase shift. Depending on these parameters, two master curves were obtained covering the most part of region, but not the transition between regions II and III where the shape of curves is phase dependent. The depend-
ences in these three regions can be described very well by use of existing analytical solutions, and transitions between them can be defined roughly by their intersection points. The second possibility of approximations is to specify the curve by two rational functions. With that a fairly precise calculation was applies to the case of phase shifts of я/2. For other phase shifts or very large ratios of vertical vibration amplitude to depth of indentation, the numerical simulation should be used as the calculation method.
Acknowledgments
The authors would like to thank V.L. Popov for his advice and contributions to develop the theoretical modelling.
References
1. Popov V.L. Contact Mechanics and Friction. -Springer, 2017.
2. Ahn Y.J. Relaxation damping and friction // Int. J. Me-ch. Sci. - 2017. - V. 128-129. - P. 147-149.
3. Peyret N., Dion J.L., Chevallier G., Argoul P. Microslip induced damping in planar contact under constant and uniform normal stress // Int. J. Appl. Mech. -2010. - V. 2(2). - P. 281-304.
4. Gagnon L., Morandini M., Ghiringhelli G.L. A review of friction damping modeling and testing // Arch. Appl. Mech. - 2020. - V. 90. - P. 107-126.
5. Mindlin R.D., Mason W.P., Osmer J.F., Deresiewicz H. Effects of an Oscillation Tangential Force on the Contact Surfaces of Elastic Spheres // Proceedings of 1st US National Congress of Applied Mechanics, 1952. -P. 203-208.
6. Popov M., Popov V.L. Pohrt R. Relaxation damping in oscillating contacts // Sci. Rep. - 2015. - V. 5. -P. 16189.
7. Davies M., Barber J.R., Hills D.A. Energy dissipation in a frictional incomplete contact with varying normal load // Int. J. Mech. Sci. - 2012. - V. 55. - P. 13-21.
8. Putignano C., Ciavarella M., Barber J.R. Frictional energy dissipation in contact of nominally flat rough surfaces under harmonically varying loads // J. Mech. Phys. Solids. - 2011. - V. 59. - P. 2442-2454.
9. Liu L., Eriten M. Frictional energy dissipation in wavy surfaces // ASME. J. Appl. Mech. - 2016. -V. 83(12). - P. 121001.
10. Johnson K.L. Surface interaction between elastically loaded bodies under tangential forces // Proc. R. Soc. Lond. A. - 1955. - V. 230(1183). - P. 531-548.
11. Goodman L.E., Brown C.B. Energy dissipation in contact friction: Constant normal and cyclic tangential loading // ASME J. Appl. Mech. - 1962. - V. 29(1). -P. 17.
12. Usta A.D., Shinde S., Eriten M. Experimental investigation of energy dissipation in presliding spherical con-
tacts under varying normal and tangential loads // ASME J. Tribol. - 2017. - V. 139(6). - P. 061402.
13. Griffin J.H., Menq C.-H. Friction damping of circular motion and its implications to vibration control // ASME J. Vib. Acoust. - 1991. - V. 113(2). - P. 225229.
14. Jang Y.H., Barber J.R. Effect of phase on the frictional dissipation in systems subjected to harmonically varying loads // Eur. J. Mech. A. Solids. - 2011. -V. 30(3). - P. 269-274.
15. Popov V.L., Heß M. Method of Dimensionality Reduction in Contact Mechanics and Friction. - Berlin-Heidelberg: Springer, 2015.
16. Popov V.L., Hess M. Method of dimensionality reduction in contact mechanics and friction: A users handbook. I. Axially-symmetric contacts // Facta Universi-tat. Mech. Eng. - 2014. - V. 12(1). - P. 1-14.
17. Popov V.L., Heß M, Willert E. Handbook of Contact Mechanics. Exact Solutions of Axisymmetric Contact Problems. - Berlin: Springer, 2019.
Received 16.03.2020, revised 16.03.2020, accepted 24.03.2020
Сведения об авторах
Thilo Hanisch, Master of Science (M.Sci.), Researcher, Technische Universität Berlin, Germany, hanisch@tu-berlin.de Ina Richter, Bachelor of Science (B.Sci.), Technische Universität Berlin, Germany, i.richter@campus.tu-berlin.de Qiang Li, Dr.-Ing., Technische Universität Berlin, Germany, qiang.li@tu-berlin.de
