Next Article in Journal
Improvement of Mechanical Properties of Composites with Surface Modified B4C for Precision Machining
Previous Article in Journal
Additive Manufacturing of 3D Anatomical Models—Review of Processes, Materials and Applications
Previous Article in Special Issue
Effect of the Transverse Functional Gradient of the Thin Interfacial Inclusion Material on the Stress Distribution of the Bimaterial under Longitudinal Shear
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Temperature during Repetitive Short-Term Operation of a Brake with Functionally Graded Friction Element

by
Aleksander Yevtushenko
,
Katarzyna Topczewska
* and
Przemysław Zamojski
Faculty of Mechanical Engineering, Bialystok University of Technology (BUT), 45C Wiejska Street, 15-351 Bialystok, Poland
*
Author to whom correspondence should be addressed.
Submission received: 22 November 2022 / Revised: 30 December 2022 / Accepted: 13 January 2023 / Published: 16 January 2023

Abstract

:
The object of study is the temperature of a braking system, operating in repetitive short-term (RST) mode. One element of the considered friction pair is made of a functionally gradient material (FGM), and the other of a homogeneous material. To determine the temperature on the friction surfaces of both elements, the previously obtained, exact solution of the boundary value problem of heat conduction was adopted, with account of the heat generation due to the friction. A calculation scheme was proposed that takes into consideration thermal sensitivity of materials and variations of the friction coefficient under the influence of temperature. Calculations were performed for two-component FGM (ZrO2–Ti-6Al-4V) in combination with gray cast iron (ChNMKh). It was found that for selected friction pair materials, consideration of their thermal sensitivity reduces the time of braking and the value of temperature achieved on the friction surfaces. At the same time, the whole process was characterized by a good stability of braking with a slight decrease in efficiency in each subsequent cycle.

1. Introduction

The operating mode of vehicles may consist of successive cycles of braking and accelerating, which is typical for driving in mountainous or urban terrains. Braking systems of vehicles in such a driving mode also work periodically, so the friction pair elements are frictionally heated during braking applications and are convectively cooled during acceleration stages [1,2,3]. However, the unforced convective heat exchange with the environment is not enough to prevent the problem of overheated friction elements during braking actions, because it has insignificant influence on the temperature distribution [4,5]. Therefore, the characteristic feature of the repetitive short-term (RST) braking mode is that the temperature increases throughout the volumes of the friction couple elements with each subsequent braking cycle. Excessive temperature conditions occurring at repeated braking processes may lead to thermal instability of the friction couple and changes in material properties, and hence a significant reduction in braking effectiveness.
In most of the published studies concerning frictional heating processes in braking systems under the RST mode, numerical methods have been used to solve the thermal problems of friction [1,2,3,4,5]. Partly, this stems from the fact that they allow for direct application of the inhomogeneous temperature state of friction elements that has been found at the end of each cycle as the initial temperature distribution in the next cycle, which is impossible to perform by means of an analytical approach. The thermal behavior of the brake disc system during single as well as RST braking modes have been studied in [4,5], with consideration of the convective heat exchange on the free surfaces of a disc. The numerical solutions to the problem of heat conduction and the corresponding quasistatic thermoelasticity problem were obtained using the finite element method (FEM). The influence of the heat transfer coefficient on the temperature and thermal stress distributions in a brake disc during RST was investigated. It was concluded that free convection cooling of the disc has no pronounced effect on the temperature and thermal stresses during the single braking event, whereas by increasing the number of braking cycles, the heat exchange process has a higher effect on the thermal state of a brake system, particularly in the case of long-lasting cooling periods [5]. The temperature generated during repetitive braking slightly drops with the rise in the heat transfer coefficient, and this relationship has been found to be approximately linear in time. Another computational model for determining the transient temperature field in a brake disc during repeated braking is proposed in [1], with a special emphasis placed on the mutual dependence of velocity and maximum temperature. The calculations were carried out with consideration of the temperature-dependent coefficient of friction and thermal sensitivity of materials. The mean, flash, maximum, and volume temperature of the disc were determined based on the equations of heat dynamics of friction and wear. It was shown that the flash temperature had the highest values at the beginning of process, and gradually decreased with time in each successive braking [1]. This is consistent with the conclusion formulated in [4] that during repetitive braking, an increase in the number of brake cycles affects the local contact conditions, which leads to the growth of the real contact area between the friction elements, where the flash temperature appears. A similar coupled calculation scheme is proposed in [2], which allows one to take into consideration the interdependence of the friction coefficient and the maximum temperature achieved during each stage of the RST braking mode. Such a formulated nonlinear thermal problem of friction was solved using the finite difference method.
Besides numerical models, analytical methods are applied for simulating the frictional heating process during repetitive short-term braking [3,6,7]. Then, in order to establish the initial condition of the thermal problem of friction for the subsequent braking cycle, the volume-averaged temperature of friction components is involved. Comparative analysis of the temperature distribution in a ventilated disc brake system during repeated braking is carried out in [3] by means of the numerical (FEM), analytical, and experimental methods. The proposed numerical model simulated the mutual motion of the stationary pad and the rotating disc by applied moving heat source. The calculated temperature field was compared with the corresponding results obtained from analytical solutions to the problems of heat conduction, as well as with the experimental data achieved from the thermocouple’s measurements. It was concluded that results determined from both theoretical methods are convergent with the experimental data [3]. Another analytical scheme to find the mean and volume temperature during repetitive short-term braking is proposed and successfully verified using experimental data in [6]. In [7], considering the thermal friction problem is considered during repetitive short-term braking mode by means of an analytical approach. A solution to the linear boundary value problem of heat conduction supported by empirical dependencies of material properties was used to calculate the mean temperature. Thermal sensitivity of the friction pair materials was partially taken into account by adjusting constant values of their properties and friction coefficient to the actual thermal state of elements in each braking application.
The above-mentioned studies concern the frictional heating of braking systems with friction elements made of homogeneous materials, or composites with spatially averaged properties. In this paper, an analytical calculation scheme is proposed to determine the temperature during repetitive short-term braking mode, in a tribosystem with a functionally graded friction element.

2. Scheme of Braking and Model Assumptions

Brake system operation during repetitive short-term (RST) mode is based on the successive performance of n cycles. Each of the full cycles k = 1 , 2 , , ( n 1 ) consists of two stages—braking and accelerating—and the last, n-th interrupted cycle has only a braking period [7,8]. In the braking stages, with constant contact pressure p 0 , the velocity of the system V ( k ) changes linearly from the initial value V 0 to the zero at the stop moment t = t s ( k ) [7]:
V ( k ) ( t ) = V 0 V ( k ) ( t ) ,   V ( k ) ( t ) = 1 t t s ( k ) ,   0 t t s ( k ) ,
t s ( k ) = 2 W 0 q 0 ( k ) A a ,   k = 1 , 2 , , n ,
where A a —nominal area of contact between the friction elements, W 0 —initial kinetic energy of the system, q 0 ( k ) —nominal value of the specific friction power. After a stop in each cycle, there is an acceleration stage, which consists of increasing the speed to the initial value, again V 0 , in the time t = t c , as follows [7]:
V ( k ) ( t ) = V 0 V ( k ) ( t ) ,   t s ( k ) t t k = t s ( k ) + t c ,
V ( k ) ( t ) = t t s ( k ) t c ,   k = 1 , 2 , , n 1 .
The full duration of the RST brake mode is equal to:
t b = t s ( 1 ) + t s ( 2 ) + + t s ( n ) + ( n 1 ) t c .
Braking stages are accompanied by intensive frictional heating of the friction elements. Before determining the resulting transient temperature field, the following simplifying assumptions were made:
  • Initial temperature of considered a friction pair at the beginning of the subsequent braking is equal to the volume-averaged temperature of the system;
  • As a result of the friction forces interaction, the heat is generated on the contact area of the elements and absorbed by them along the normal direction to the friction surface. The friction thermal contact of the elements during heating is perfect;
  • Unforced convection cooling of the system during braking stages is omitted.
According to the above assumptions, the friction pair can be formed by two semi-infinite bodies z 0 and z 0 , and the sought transient temperature field is one-dimensional, i.e., T = T ( z , t ) . The scheme of the considered tribosystem is illustrated below, in Figure 1.
Further, all variables and parameters related to the first semispace z 0 are indicated by a subscript l = 1 , and to the second semispace z 0 by a subscript l = 2 . The first semispace z 0 is made of the two-component (base and core) FGM. Suppose that K 1 , m , c 1 , m , ρ 1 , m are the thermal conductivity, specific heat, and density of the materials of the base ( m = 1 ) and core ( m = 2 ), respectively. The thermal conductivity of the FGM increases exponentially (with the gradient parameter γ 0 ) in the direction normal to the working surface of the base. The homogeneous material of the second semispace z 0 has the coefficient of thermal conductivity K 2 , specific heat c 2 , and density ρ 2 .
Additionally, it was assumed that the materials of both elements and the friction coefficient f are thermally sensitive [2]:
K 1 , m ( T ) = K 1 , m ( 0 ) K 1 , m * ( T ) ,   c 1 , m ( T ) = c 1 , m ( 0 ) c 1 , m * ( T ) ,   ρ 1 , m ( T ) = ρ 1 , m ( 0 ) ρ 1 , m * ( T ) ,   m = 1 , 2 ,
K 1 , m ( 0 ) K 1 , m ( T 0 ) ,   c 1 , m ( 0 ) c 1 , m ( T 0 ) ,   ρ 1 , m ( 0 ) ρ 1 , m ( T 0 ) ,   m = 1 , 2 ,
K 2 ( T ) = K 2 ( 0 ) K 2 * ( T ) ,   c 2 ( T ) = c 2 ( 0 ) c 2 * ( T ) ,   ρ 2 ( T ) = ρ 2 ( 0 ) ρ 2 * ( T ) ,
K 2 ( 0 ) K 2 ( T 0 ) ,   c 2 ( 0 ) c 2 ( T 0 ) ,   ρ 2 ( 0 ) ρ 2 ( T 0 ) ,
f ( T ) = f 0 f ( T ) ,   f 0 f ( T 0 ) ,
where T 0 —initial temperature of the system, and the corresponding dimensionless temperature functions are marked with the superscript ‘ ’. Typically, dependencies of the type (6)–(10) are obtained as a result of processing appropriate experimental data [9].

3. Analytical Model

A key element of the proposed mathematical model is the solution to the appropriate thermal problem of friction for the considered system during execution of the individual braking processes. On the basis of the above qualitative assumptions of such a problem, a non-linear boundary value problem of heat conduction can be formulated with account of heat generation due to friction. The considered system consists of two semi-infinite bodies, one of which is made of a FGM and the other of a homogeneous material. Unfortunately, the significant nonlinearity of such a problem, caused by the thermal sensitivity of the friction coefficient as well as mechanical and thermal properties, means that the solutions can be obtained only by the numerical methods [1,2].
Another approach is also known, consisting of adapting an appropriate solution to the linear problem to determine the temperature of the thermally sensitive braking system [8]. It has been carried out for the case of a friction pair made of homogeneous materials [7]. This study presents an algorithm to establish the temperature of a thermally sensitive braking system with a friction couple consisting of a functionally graded element in combination with a homogeneous one. For such a friction pair, an exact solution of the linear, thermal problem of friction was obtained (with the coefficient of friction and thermophysical properties unchanged) during single braking with a constant deceleration [10]. Based on such a solution, the evolution of the temperature on the friction surfaces of the system operating in RST mode during the subsequent k –th braking can be calculated from the formulas [10]:
T ( k ) ( t ) = T ^ ( k ) + Λ ( k ) T ( k ) ( τ ) ,   0 t t s ( k ) ,   k = 1 , 2 , , n ,
where
Λ ( k ) = q 0 ( k ) a ( k ) K 1 , 1 ( k ) ,
T ( k ) ( τ ) = 1 γ ( k ) [ 1 4 π 0 G ( k ) ( x ) P ( k ) ( τ , x ) d x ] , 0 τ τ s ( k ) ,
G ( k ) ( x ) = K ε ( k ) [ J 1 ( x ) ] 2 x 2 { [ J 0 ( x ) ] 2 + [ K ε ( k ) J 1 ( x ) ] 2 } ,
P ( k ) ( τ , x ) = e X k τ ( 1 e X k τ ) X k τ s ( k ) ,   X k = ( γ ( k ) x ) 2 4 ,
τ = k 1 ( k ) t ( a ( k ) ) 2 ,   τ s ( k ) = k 1 ( k ) t s ( k ) ( a ( k ) ) 2 ,   k 1 ( k ) = K 1 , 1 ( k ) ρ 1 ( k ) c 1 ( k ) ,   k 2 ( k ) = K 2 ( k ) ρ 2 ( k ) c 2 ( k ) ,
c 1 ( k ) = v c 1 , 1 ( k ) + ( 1 v ) c 1 , 2 ( k ) ,   ρ 1 ( k ) = v ρ 1 , 1 ( k ) + ( 1 v ) ρ 1 , 2 ( k ) ,   0 v 1 ,
K ε ( k ) = K ( k ) k ( k ) ,   K ( k ) = K 2 ( k ) K 1 , 1 ( k ) ,   k ( k ) = k 2 ( k ) k 1 , 1 ( k ) ,   γ ( k ) = ln ( K 1 , 2 ( k ) K 1 , 1 ( k ) ) ,
K 1 , m ( k ) K 1 , m ( T ^ ( k ) ) ,   c 1 , m ( k ) c 1 , m ( T ^ ( k ) ) ,   ρ 1 , m ( k ) ρ 1 , m ( T ^ ( k ) ) ,   m = 1 , 2 ,
K 2 ( k ) K 2 ( T ^ ( k ) ) ,   c 2 ( k ) c 2 ( T ^ ( k ) ) ,   ρ 2 ( k ) ρ 2 ( T ^ ( k ) ) ,
q 0 ( k ) = f ( k ) p 0 V 0 ,   f ( k ) f ( T ^ ( k ) ) ,
a ( k ) = max ( a 1 ( k ) , a 2 ( k ) ) ,
a l ( k ) = { d l , a l , e f f ( k ) d l a l , e f f ( k ) , a l , e f f ( k ) < d l ,   a l , e f f ( k ) = 3 k l ( k ) t s ( k ) ,   l = 1 , 2 ,
where J n ( x ) , n = 0 , 1 —Bessel functions of the first kind; d l , l = 0 , 1 —thickness of the friction elements (e.g., pad and disc).
The volume temperature T ^ ( k ) of the friction system before the start of the k –th braking was found as [7]:
T ^ ( k ) = 0.5 ( T ^ 0 ( k ) + T ^ 1 ( k ) ) ,   k = 1 , 2 , , n ,
where
T ^ i ( k ) = T 0 + α i ( k ) W 0 2 G 2 , i ( k ) c 2 , i ( k ) ( e β i ( k ) t c e k β i ( k ) t c 1 e β i ( k ) t c ) ,   i = 0 , 1 ,
β i ( k ) = h A vent 2 G 2 , i ( k ) c 2 , i ( k ) ,
G 2 , i ( k ) = A 2 d 2 ρ 2 ( i , k ) ,
c 2 ( 0 , k ) = c 2 ( 0 ) ,   c 2 ( 1 , k ) = c 2 ( T ^ 0 ( k ) ) ,   ρ 2 ( 0 , k ) = ρ 2 ( 0 ) ,   ρ 2 ( 1 , k ) = ρ 2 ( T ^ 0 ( k ) ) ,
where h —coefficient of the convective heat transfer from the surface of the disc with an area A vent during the acceleration stages, α i ( k ) —heat partition ratio (HPR). The methodology for determining HPR for the functionally graded friction couple is proposed in [11]. Based on this methodology, the heat transfer coefficient in formula (25) for the considered friction pair (FGM—homogeneous material) was found in the form [11]:
α i ( k ) = ε i ( k ) 1 + ε i ( k ) ,   ε i ( k ) = 0.625 d * K i ( k ) γ i ( k ) π τ 2 ( i , k ) ,   k = 1 , 2 , , n ,   i = 0 , 1 ,
where
K i ( k ) = K 2 ( i , k ) K 1 , 1 ( i , k ) ,   γ i ( k ) = ln ( K 1 , 2 ( i , k ) K 1 , 1 ( i , k ) ) ,   d * = d 1 d 2 ,
τ 2 ( i , k ) = k 2 ( i , k ) t s ( k ) d 2 2 ,   k 2 ( i , k ) = K 2 ( i , k ) ρ 2 ( i , k ) c 2 ( i , k ) ,
K 1 , m ( 0 , k ) = K 1 , m ( 0 ) ,   K 1 , m ( 1 , k ) = K 1 , m ( T ^ 0 ( k ) ) ,   m = 1 , 2 ,   K 2 ( 0 , k ) = K 2 ( 0 ) ,   K 2 ( 1 , k ) = K 2 ( T ^ 0 ( k ) ) .
It should be noted that from formulas (24) and (25), it follows that before the start of the first braking ( k = 1 ), the volume temperature T ^ ( 1 ) is equal to the initial temperature of the system T 0 . Before the start of subsequent braking, when determining the volume temperature T ^ ( k ) , k = 2 , , n the first component T ^ 0 ( k ) in formula (24) was also established using the properties of the materials at the initial temperature T 0 , while the second component T ^ 1 ( k ) was used to correct the result by taking into consideration the thermal sensitivity of the materials.

4. Numerical Analysis

The following scheme for determination of the temperature evolution on the working surfaces of selected friction pair elements was proposed:
  • Based on experimental data, finding the dependences of material properties and the friction coefficient on temperature in forms (6)–(10). Determining the value of material properties K 1 , m ( 0 ) , c 1 , m ( 0 ) , ρ 1 , m ( 0 ) , m = 1 , 2 , (7), K 2 ( 0 ) , c 2 ( 0 ) , ρ 2 ( 0 ) (9) and the coefficient of friction f 0 (10) at the initial temperature T 0 ;
  • Introduction of the input operational parameters: p 0 , V 0 , T 0 , W 0 , n , A a , A vent , d 1 , d 2 , h , t c , v ;
  • Start of the first braking: k = 1 ;
  • Determination of the volume temperature T ^ ( k ) (24)–(32);
  • Using the dependencies (6)–(10), establishment of the material properties values K 1 , m ( k ) , c 1 , m ( k ) , ρ 1 , m ( k ) , m = 1 , 2 (19), K 2 ( k ) , c 2 ( k ) , ρ 2 ( k ) (20), the friction coefficient f ( k ) , and specific friction power q 0 ( k ) (21) at the volume temperature T ^ ( k ) ;
  • Determination of the stop time t s ( k ) (2) and temporal profile of velocity V ( k ) ( t ) , 0 t t s ( k ) (1);
  • Calculation of the temperature evolution T ( k ) ( t ) , 0 t t s ( k ) (11)–(23);
  • Starting the next k + 1 braking cycle and repeating starting from point (5) or ending the calculation process after reaching the equality k = n .
The above scheme was performed for a selected friction pair, which the first element is made of the two-component FGM: zirconium dioxide ZrO 2 (base, m = 1 )—titanium alloy Ti - 6 Al - 4 V (core, m = 2 ), and the second homogeneous element is made of the gray cast iron ChNMKh.
The properties (7) and (9) of these materials at the initial temperature T 0 = 20 C are as follows [10,11]:
ZrO 2
K 1 , 1 ( 0 ) = 1.94 W m 1 K 1 ,   c 1 , 1 ( 0 ) = 452.83 J kg 1 K 1 ,   ρ 1 , 1 ( 0 ) = 6102.16 kg m 3 ,
Ti - 6 Al - 4 V
K 1 , 2 ( 0 ) = 6.87 W m 1 K 1 ,   c 1 , 2 ( 0 ) = 538.08 J kg 1 K 1 ,   ρ 1 , 2 ( 0 ) = 4431.79 kg m 3 ,
ChNMKh
K 2 ( 0 ) = 52.17 W m 1 K 1 ,   c 2 ( 0 ) = 444.6 J kg 1 K 1 ,   ρ 2 ( 0 ) = 7100 kg m 3 .
Dependencies of material properties on the temperature have the forms:
ZrO 2 [12,13,14,15]
K 1 , 1 ( T ) = 1.9365 + 0.7 10 4 T + 0.5 10 - 6   T 2 0.2 10 9 T 3 ,
c 1 , 1 ( T ) = 437.96 + 0.7767 T 0.17 10 2 T 2 ,
ρ 1 , 1 ( T ) = 6104.6 0.1212 T 0.4 10 4 T 2 + 0.3 10 7 T 3 0.1 10 10 T 4 ,  
Ti 6 Al 4 V [16,17]
K 1 , 2 ( T ) = 6.6926 + 8 . 9177 10 - 3 T + 6.8432 10 6 T 2 ,
c 1 , 2 ( T ) = 529.9316 + 0.4154 T 4.01646 10 4 T 2 + 1.6364 10 7 T 3 ,
ρ 1 , 2 ( T ) = 4434 0.1088 T 0.8 10 4 T 2 + 10 7 T 3 0.6 10 10 T 4 ,
ChNMKh [18]
K 2 ( T ) = 53.24 0.028 T ,
c 2 ( T ) = 432.43 + 0.559 T + 2.712 10 4 T 2 + 1.657 10 6 T 3 + 9.439 10 10 T 4 ,
ρ 2 ( T ) = ρ 2 ( 0 ) .
The dependence of the friction coefficient on temperature for the considered friction pair has the form (10), where [19]
f 0 = 0.27 ,   f ( T ) = e κ ( T 1 ) ,   T = T / T 0 ,   κ = 0.7 10 4 .
Graphs of the corresponding dimensionless functions K 1 , m ( T ) = K 1 , m ( T ) / K 1 , m ( 0 ) , c 1 , m ( T ) = c 1 , m ( T ) / c 1 , m ( 0 ) , ρ 1 , m ( T ) = ρ 1 , m ( T ) / ρ 1 , m ( 0 ) , m = 1 , 2 , K 2 ( T ) = K 2 ( T ) / K 2 ( 0 ) , c 2 ( T ) = c 2 ( T ) / c 2 ( 0 ) , ρ 2 ( T ) = ρ 2 ( T ) / ρ 2 ( 0 ) and f ( T ) are presented in Figure 2.
Calculations were made with the following operating parameters [7]:
p 0 = 1.47 MPa , V 0 = 27.78 m s 1 , W 0 = 392.1 kJ , n = 5 , A a = 4.05 10 2 m 2 , A vent = 4.44 10 2 m 2 , d 1 = 5.5 mm , d 2 = 11 mm , h = 100 Wm 1 K 1 , t c = 5 s , v = 0 . 5 .
Temporal profiles of the velocity V ( k ) ( t ) (1), (2) and specific friction power q ( k ) ( t ) = f ( k ) p 0 V ( k ) ( t ) , 0 t t s ( k ) , k = 1 , 2 , , 5 are illustrated in Figure 3. A noticeable effect is the extension of the braking stage in each subsequent cycle of the process (Figure 3a). During each of the five braking applications, the intensity of the performed friction work (equal to the area under the graph) is the same (Figure 3b). This fact made it possible to compare the relevant temperature evolutions, demonstrated in Figure 4. This figure shows a comparison of friction surface temperature changes during braking T ( k ) ( t ) (11)–(24), found with (solid lines) and without (dotted lines) consideration of the dependencies (37)–(45) of material properties on the temperature. Results corresponding to the dotted curves were obtained for the properties of materials (34)–(36) at the initial temperature T 0 = 20 C . In both variants, temperature changes of the friction coefficient were taken into account in form (46).
With the exception of the first braking, consideration of the materials’ thermal sensitivity resulted in a drop of the temperature on the friction surfaces. This effect is most noticeable in the last, fifth braking. Calculated values of the friction coefficient f ( k ) (21), time of braking t s ( k ) (2), volume temperature T ^ ( k ) (24)–(32), and maximum temperature for each of the five braking actions, obtained with account of thermal sensitivity, are presented in Table 1. Corresponding data found for constant values of material properties are demonstrated in Table 2. Additionally, the data from Table 1 are presented in graphical form in Figure 5. With each successive braking, the coefficient of friction f ( k ) at the volume temperature T ^ ( k ) decreases, while the braking time t s ( k ) , volume temperature T ^ ( k ) , and maximum temperature T max ( k ) increase. Consideration of materials’ thermal sensitivity in the proposed analytical model results in greater stability of the friction coefficient value, shorter time of braking stages, lower values of volume, and maximum temperature, compared to the corresponding data found with unchanged material properties. The differences in the temperature values, obtained with and without account of the thermal sensitivity of the friction pair materials, increase with each successive braking.
Due to the fact that the curve of friction heat resistance (Figure 2d) is monotonically increasing the function of the temperature, the change in the friction coefficient during successive braking has the opposite form to the evolution of temperature (Figure 6). At the beginning of each braking cycle, the coefficient of friction is reduced until the maximum temperature T max ( k ) is reached. In the subsequent period of the temperature drop, which lasts until the stop time t s ( k ) , the coefficient of friction slightly grows. With the increase in the number of braking applications, the minimum value of f ( k ) decreases.
On the basis of the results shown in Figure 6, the parameters characterizing the operation of the braking system during the subsequent cycle, such as the average value of the friction coefficient f m ( k ) , its stability f s ( k ) = f m ( k ) / f max ( k ) , fluctuation f f ( k ) = f min ( k ) / f max ( k ) , and braking efficiency f e f f ( k ) = f s ( k ) / ( t s ( k ) ) 2 , were determined (Table 3). All braking cycles were characterized by good stability, and the first braking turned out to be the most effective for the selected friction pair. With each subsequent braking, the efficiency decreases.

5. Conclusions

An analytical scheme was proposed to determine the temperature during the repeated short-term (RST) operation mode of the braking system, in which one of the friction elements is made of a functionally gradient material (FGM). The proposed approach is a generalization of the authors’ results concerning a single braking process [10,11,12]. Calculations were carried out for a friction pair made of a two-component FGM (base ZrO2, core Ti–6Al–4V) and ChNMKh gray cast iron, for five braking actions. It was found that the braking time, the volume, and maximum temperature values increased almost linearly with the number of braking cycles. Involving the thermal sensitivity of materials into the calculation model causes a decrease in the maximum temperature value in relation to the results obtained for materials with invariant properties under temperature changes. This effect becomes more noticeable with each subsequent braking cycle. The coefficient of friction decreases rapidly at the beginning of each braking to a minimum value, then begins to increase slightly until standstill. The considered friction pair is characterized by good braking stability with sufficient efficiency, slightly decreasing with each successive braking.
It should be noted that the problem of determining the effect of FGM on temperature is currently intensively developed not only for bodies with unidirectional heat extension. An exact solution for transient heat conduction problem in an axisymmetric cylinder made of FGM whose thermal conductivity differs in two (radial and longitudinal) directions has been obtained [20]. Another analytical solution for steady-state heat transfer in a hollow sphere made of functionally graded material has been proposed [21].

Author Contributions

Conceptualization and methodology, A.Y.; software, P.Z.; validation, A.Y. and K.T.; formal analysis, P.Z.; investigation, A.Y., K.T., and P.Z.; writing—original draft preparation, A.Y. and K.T.; writing—review and editing, K.T.; visualization and figure preparation, P.Z.; supervision, A.Y.; project administration, K.T. All authors have read and agreed to the published version of the manuscript.

Funding

This investigation was performed within the framework of research project No. 2017/27/B/ST8/01249, funded by the National Science Centre, Poland, and project financing was through the program of the Minister of Education and Science of Poland named “Regional Initiative of Excellence” in 2019–2022, project No. 011/RID/2018/19; amount financed—PLN 12,000,000.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

a Effective depth of heat penetration of the friction element ( m   )
A a Area of the nominal contact region (m2)
A v e n t Area of the ventilated surface of the brake disc (m2)
c Specific heat ( J kg 1 K 1 )
dThickness of friction elements ( m   )
f Coefficient of friction (dimensionless)
hCoefficient of heat transfer ( W m 2 K 1 )
J k ( ) Bessel functions of the first kind of the kth order
k Thermal diffusivity ( m 2 s 1 )
K Thermal conductivity ( W m 1 K 1 )
K ε Dimensionless coefficient of thermal activity of friction couple
nNumber of brakings in RST brake mode
p Pressure on the contact surface (Pa)
p 0 Nominal value of the contact pressure (Pa)
q Specific friction power ( W m 2 )
q 0 Nominal value of specific friction power ( W m 2 )
t Time ( s )
tbTime of performance of all RST mode of braking (s)
tcCooling time during acceleration stage (s)
t s Stop time ( s )
T Temperature ( C )
T Dimensionless temperature
T ^ Volume temperature ( C )
T 0 Initial temperature ( C )
vVolume fraction of the FGM components (dimensionless)
V Velocity ( m s 1 )
V 0 Initial velocity ( m s 1 )
W0Initial kinetic energy (J)
z Spatial coordinate in axial direction ( m   )
α Heat partition ratio (dimensionless)
γ Parameter of material gradient ( m   - 1 )
Λ Scaling factor of temperature ( C )
ρ Density ( kg m 3 )
τ Dimensionless time
τ s Dimensionless time of braking
superscript kNumber of a braking cycle
subscript lNumber of the friction pair element
subscript mNumber of the component material of functionally graded element

References

  1. Grześ, P. Maximum temperature of the disc during repeated braking applications. Adv. Mech. Eng. 2019, 11, 1–13. [Google Scholar] [CrossRef]
  2. Yevtushenko, A.; Kuciej, M. Calculation of friction characteristics of disc brakes used in repetitive short-term braking mode. J. Frict. Wear 2020, 41, 509–516. [Google Scholar] [CrossRef]
  3. Yevtushenko, A.A.; Kuciej, M.; Grześ, P.; Wasilewski, P. Temperature in the railway disc brake at a repetitive short-term mode of braking. Int. Commun. Heat Mass Transf. 2017, 84, 102–109. [Google Scholar] [CrossRef]
  4. Adamowicz, A.; Grześ, P. Analysis of disc brake temperature distribution during single braking under non-axisymmetric load. Appl. Therm. Eng. 2011, 31, 1003–1012. [Google Scholar] [CrossRef] [Green Version]
  5. Adamowicz, A. Effect of convective cooling on temperature and thermal stresses in disk during repeated intermittent braking. J. Frict. Wear 2016, 37, 107–112. [Google Scholar] [CrossRef]
  6. Dunaevsky, V.V. Prediction of railroad friction braking temperatures: Prediction of average bulk and average surface temperatures of railroad wheels and brake discs. Tribol. Trans. 1991, 34, 343–352. [Google Scholar] [CrossRef]
  7. Yevtushenko, A.; Topczewska, K.; Kuciej, M. Analytical Determination of the Brake Temperature Mode during Repetitive Short-Term Braking. Materials 2021, 14, 1912. [Google Scholar] [CrossRef] [PubMed]
  8. Chichinadze, A.V.; Braun, E.D.; Ginzburg, A.G.; Ignat’eva, Z.V. Calculation, Testing and Selection of Friction Couples; Nauka: Moscow, Russian, 1979. (In Russian) [Google Scholar]
  9. Chichinadze, A.V. Polymers in Friction Assembles of Machines and Devices: A Handbook; Allerton Press Inc.: New York, NY, USA, 1984. [Google Scholar]
  10. Yevtushenko, A.; Kuciej, M.; Topczewska, K.; Zamojski, P. Temperature in the Friction Couple Consisting of Functionally Graded and Homogeneous Materials. Materials 2022, 15, 3600. [Google Scholar] [CrossRef] [PubMed]
  11. Yevtushenko, A.; Topczewska, K.; Zamojski, P. The Heat Partition Ratio during Braking in a Functionally Graded Friction Couple. Materials 2022, 15, 4623. [Google Scholar] [CrossRef] [PubMed]
  12. Yevtushenko, A.; Topczewska, K.; Zamojski, P. Influence of Thermal Sensitivity of Functionally Graded Materials on Temperature during Braking. Materials 2022, 15, 963. [Google Scholar] [CrossRef] [PubMed]
  13. Kingery, W.D.; Francl, J.; Coble, R.L.; Vasilos, T. Thermal Conductivity: X, Data for Several Pure Oxide Materials Corrected to Zero Porosity. J. American Ceramic Soc. 1954, 37, 107–110. [Google Scholar] [CrossRef]
  14. Pankratz, L.B. Thermodynamic Properties of Elements and Oxides-Bulletin 672; U.S. Bureau of Mines: Washington, DC, USA, 1982.
  15. Taylor, D. Thermal expansion data. II: Binary oxides with the fluorite and rutile structures, MO2, and the antifluorite structure, M2O. Trans. J. British Ceramic Soc. 1984, 83, 32–37. [Google Scholar]
  16. Deem, H.W.; Wood, W.D.; Lucks, C.F. The relationship between Electrical and Thermal Conductivities of Titanium Trans. Metall. Soc. AIME 1958, 212, 520–523. [Google Scholar]
  17. Cezairrliyan, A.; McClure, J.L.; Taylor, R.J. Thermophysical Measurements on 90Ti-6Al-4V Alloy Above 1450 K Using a Transient (Subsecond) Technique. J. Res. Nat. Bur. Stand.–A Phys. Chem. 1977, 81, 251–256. [Google Scholar] [CrossRef]
  18. Chichinadze, A.V.; Matveevskii, R.M.; Braun, E.P. Materials in Triboengineering of Unsteady Processes; Nauka: Moscow, Russia, 1986. [Google Scholar]
  19. Fu, P.; Zhao, J.; Zhang, X.; Kang, G.; Wang, P.; Kan, Q. Thermo-mechanically coupled sliding contact shakedown analysis of functionally graded coating-substrate structures. Int. J. Mech. Sci. 2022, 222, 107241. [Google Scholar] [CrossRef]
  20. Delouei, A.A.; Emamian, A.; Karimnejad, S.; Sajjadi, H. A closed-form solution for axisymmetric conduction in a finite functionally graded cylinder. Int. Commun. Heat Mass Transf. 2019, 108, 104280. [Google Scholar] [CrossRef]
  21. Delouei, A.A.; Emamian, A.; Karimnejad, S.; Sajjadi, H.; Jing, D. Two-dimensional analytical solution for temperature distribution in FG hollow spheres: General thermal boundary conditions. Int. Commun. Heat Mass Transf. 2020, 113, 104531. [Google Scholar] [CrossRef]
Figure 1. Scheme of the system.
Figure 1. Scheme of the system.
Materials 16 00881 g001
Figure 2. Temperature dependencies of the dimensionless material properties: (a) coefficient of thermal conductivity; (b) specific heat; (c) density; and (d) friction coefficient of the considered friction pair.
Figure 2. Temperature dependencies of the dimensionless material properties: (a) coefficient of thermal conductivity; (b) specific heat; (c) density; and (d) friction coefficient of the considered friction pair.
Materials 16 00881 g002
Figure 3. Evolutions of the: (a) velocity; (b) specific friction power, during each of the five braking applications.
Figure 3. Evolutions of the: (a) velocity; (b) specific friction power, during each of the five braking applications.
Materials 16 00881 g003
Figure 4. Evolutions of temperature on the friction surface during each of the five braking actions with (solid lines) and without (dotted lines) taking into account the thermal sensitivity of materials.
Figure 4. Evolutions of temperature on the friction surface during each of the five braking actions with (solid lines) and without (dotted lines) taking into account the thermal sensitivity of materials.
Materials 16 00881 g004
Figure 5. Dependencies on the number of braking applications k of: (a) friction coefficient f ( k ) (21); (b) braking time t s ( k ) (2); (c) volume temperature T ^ ( k ) (24); (d) maximum temperature on the friction surface T max ( k ) .
Figure 5. Dependencies on the number of braking applications k of: (a) friction coefficient f ( k ) (21); (b) braking time t s ( k ) (2); (c) volume temperature T ^ ( k ) (24); (d) maximum temperature on the friction surface T max ( k ) .
Materials 16 00881 g005
Figure 6. Evolution of friction coefficient f ( k ) (21) during five successive braking cycles.
Figure 6. Evolution of friction coefficient f ( k ) (21) during five successive braking cycles.
Materials 16 00881 g006
Table 1. Calculation results obtained with consideration of the materials’ thermal sensitivity.
Table 1. Calculation results obtained with consideration of the materials’ thermal sensitivity.
k12345
f ( k ) 0.270.250.230.210.20
t s ( k ) , s 1.771.942.102.262.41
T ^ ( k ) , ° C 2086.85143.21193.48239.74
T max ( k ) , ° C 530.27561.29591.46621.11650.49
Table 2. Calculation results obtained with constant properties of materials.
Table 2. Calculation results obtained with constant properties of materials.
k12345
f ( k ) 0.270.240.220.200.18
t s ( k ) , s 1.771.952.152.362.60
T ^ ( k ) , ° C 2089.83158.79226.89294.12
T max ( k ) , ° C 530.27575.76621.75668.19715.04
Table 3. Parameters of the braking process evaluation.
Table 3. Parameters of the braking process evaluation.
k12345
f m ( k ) 0.1490.1420.1350.1290.123
f s ( k ) 0.5530.5760.5930.6080.620
f f ( k ) 0.4890.5150.5340.5500.563
f e f f ( k ) , s 2 0.1770.1530.1340.1190.107
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Yevtushenko, A.; Topczewska, K.; Zamojski, P. Temperature during Repetitive Short-Term Operation of a Brake with Functionally Graded Friction Element. Materials 2023, 16, 881. https://0-doi-org.brum.beds.ac.uk/10.3390/ma16020881

AMA Style

Yevtushenko A, Topczewska K, Zamojski P. Temperature during Repetitive Short-Term Operation of a Brake with Functionally Graded Friction Element. Materials. 2023; 16(2):881. https://0-doi-org.brum.beds.ac.uk/10.3390/ma16020881

Chicago/Turabian Style

Yevtushenko, Aleksander, Katarzyna Topczewska, and Przemysław Zamojski. 2023. "Temperature during Repetitive Short-Term Operation of a Brake with Functionally Graded Friction Element" Materials 16, no. 2: 881. https://0-doi-org.brum.beds.ac.uk/10.3390/ma16020881

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop