*Article* **Analytical Determination of the Brake Temperature Mode during Repetitive Short-Term Braking**

**Aleksander Yevtushenko, Katarzyna Topczewska and Michal Kuciej \***

Faculty of Mechanical Engineering, Bialystok University of Technology (BUT), 45C Wiejska Street, 15-351 Białystok, Poland; a.yevtushenko@pb.edu.pl (A.Y.); k.topczewska@pb.edu.pl (K.T.) **\*** Correspondence: m.kuciej@pb.edu.pl

**Abstract:** An algorithm to determine the maximum temperature of brake systems during repetitive short-term (RST) braking mode has been proposed. For this purpose, the intermittent mode of braking was given in the form of a few cyclic stages consisting of subsequent braking and acceleration processes. Based on the Chichinadze's hypothesis of temperature summation, the evolutions of the maximum temperature during each cycle were calculated as the sum of the mean temperature on the nominal contact surface of the friction pair elements and temperature attained on the real contact areas (flash temperature). In order to find the first component, the analytical solution to the one-dimensional thermal problem of friction for two semi-spaces taking into account frictional heat generation was adapted. To find the flash temperature, the solution to the problem for the semi-infinite rod sliding with variable velocity against a smooth surface was used. In both solutions, the temperature-dependent coefficient of friction and thermal sensitivity of materials were taken into account. Numerical calculations were carried out for disc and drum brake systems. The obtained temporal variations of sliding velocity, friction power and temperature were investigated on each stage of braking. It was found that the obtained results agree well with the corresponding data established by finite element and finite-difference methods.

**Keywords:** repetitive short-term braking; frictional heating; temperature; thermal sensitivity of materials; friction coefficient

### **1. Introduction**

Repetitive short-term mode of braking (RST) is a sequential performance of a certain number of cycles consisting of two stages: braking (heating) and acceleration (cooling) [1]. In analytical and numerical models, the main characteristics to describe the braking process (temporal profiles of velocity and friction power, time to standstill, maximum temperature, etc.) are found from solutions to the following problems [2,3]:


The schemes of the characteristics' determination, consisting of subsequent solutions to the above-mentioned problems, are the so-called uncoupled models. In this kind of scheme, to determine the sliding velocity profile, time to stop and resulting evolution of specific friction power, first, the initial value problem for the equation of motion of braking system has to be considered. After finding the specific friction power in this manner, its temporal profile is adapted to one of the boundary conditions during formulation of the boundary-value problem of heat conduction. The equality of the specific friction power and the sum of heat flux intensities, directed inside each element of the friction pair perpendicular to the contact surface, is required in this condition. The disadvantage of such an approach is that, when obtaining solutions, a constant, usually averaged with the

**Citation:** Yevtushenko, A.; Topczewska, K.; Kuciej, M. Analytical Determination of the Brake Temperature Mode during Repetitive Short-Term Braking. *Materials* **2021**, *14*, 1912. https://doi.org/10.3390/ ma14081912

Academic Editor: Chuang Dong

Received: 15 February 2021 Accepted: 7 April 2021 Published: 11 April 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Copyright:** © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/).

braking time, the value of the friction coefficient is used. In consequence, the calculations do not take into account the mutual influence of sliding velocity and temperature.

Numerical uncoupled models of braking in the RST mode, based on the finite element method (FEM), were developed for automotive [4–7] and railway [8] brakes. The influence of convective cooling of free surfaces on temperature [4] and thermal stresses were studied in [5–7]. A comparative analysis of volumetric and mean temperature on the friction surface of the disc brake of the diesel multiple unit, obtained by means of: (a) numerical solution to the spatial thermal problem of friction for a disc-pad system, (b) analytical solution to the corresponding one-dimensional problem [9], and (c) experimental data using thermocouples [10] was carried out in [8]. To predict temperature, a computer model of an automotive disc brake operation during short-term multiple brakings was developed [11]. A theoretical and experimental methodology for determining the heat transfer coefficient during the RST braking mode has been proposed [12].

However, if for a selected friction pair the experimental data of frictional thermal stability are known, then the temperature calculations can be made according to the coupled models [13,14]. In this kind of model, the above-mentioned problems are related to the temperature-dependent coefficient of friction, and their solutions at each time step are obtained simultaneously. The implementation of the coupled calculation scheme allows to take into account the interdependence of sliding velocity and temperature during each stage of the RST braking mode. Since both of the above-mentioned problems in the coupled models are non-linear, their solutions were obtained numerically, using FEM [15] or the finite differences method (FDM) [16].

The use of numerical methods allows to take into account the actual dimensions of the braking system, including the contact area of the pad and disc, as well as the cooling of free surfaces, thermal sensitivity of materials, etc. On the other hand, the use of numerical methods requires researcher to have advanced skills when selecting the space–time grid, controlling the stability of calculations, performing verification of the obtained results with appropriate experimental data, etc. Therefore, analytical methods of solving thermal problems of friction during braking are still being investigated, simultaneously with numerical and experimental methods. These are mainly linear boundary value problems of heat conduction [17]. The exact obtained solutions to such problems during single braking cases are in the form of closed engineering formulas and allow for instant assessment of the temperature during braking under light conditions with sufficient accuracy; when the volumetric and the average temperatures of the friction surface do not exceed 100 ◦C and 200 ◦C, respectively [18–20]. Under heavy and medium braking conditions, when the volume temperature reaches values above 250 ◦C, the friction coefficient and properties of the friction pair materials change significantly due to the influence of temperature; therefore, numerical methods are used [21,22]. An attempt was made to use the exact solutions of linear thermal problems of friction to determine temperature during single braking of a system made of thermally sensitive friction materials [23–25].

It should be noted that, in the sense of methodology, the results contained in the monograph by A.V. Chichinadze et al. [9] and in [16] are closest to our study. The algorithm for determining the maximum temperature of a braking system in [9] is based on an analytical solution to the thermal problem of friction during single braking for two layers. This solution was obtained with two very important simplifying assumptions. The first one is the presumption that the temperature of any point on the axis, perpendicular to the friction surface, is directly proportional to the braking time. The second assumption is that the temperature increase, in this case, is equal to the increase in the mean volumetric temperature of the system. With such simplifications, this solution can be defined as approximate in an analytical form. A detailed comparative analysis of the temperature fields during single braking, obtained with the use of an approximate solution [9] and the exact solution, was carried out in [18]. The numerical FDM solution to the thermal problem of friction for the RST operating mode of a disc brake was obtained in [16]. This solution also belongs to the class of approximate solutions; whereas the main purpose of this study

was to obtain exact engineering formulae to this problem and the development, on the basis of them, of a corresponding calculation scheme to find the maximum temperature in a braking system operating in RST mode.

### **2. Statement of the Problem**

A braking system operating in a repetitive short-term (RST) mode is considered. This intermittent braking mode comprised the sequential performance of *n* − 1 full cycles and the final (*n*th) discontinuous cycle (Figure 1). Each of the full cycles consisted of two stages. The first one was braking, where the linear velocity of sliding, *V*(*k*), decreased from the initial value, *V*0, to zero at time moments *t* = *t* (*k*) *<sup>s</sup>* , *k* = 1, 2, ... , *n* − 1. Immediately after the standstill, the second phase of the cycle began, which was acceleration with the released brake. During this stage the velocity increased to initial value *V*<sup>0</sup> at time *t* = *tc*.

**Figure 1.** Scheme of sliding velocity variations during repetitive short-term (RST) braking mode.

Thus, the duration of one full cycle was *tk* = *t* (*k*) *<sup>s</sup>* + *tc*, *k* = 1, 2, ... , *n* − 1, and the overall time of all full braking cycles was *tn*−<sup>1</sup> <sup>=</sup> *<sup>n</sup>*−<sup>1</sup> ∑ *k*=1 *tk*. The last interrupted cycle was performed with the applied brake, so speed *V*(*n*) dropped from initial value *V*<sup>0</sup> to zero at stop moment *t* = *t* (*n*) *<sup>s</sup>* . Thus, the total time of the RST operation mode of braking was *tb* = *tn*−<sup>1</sup> + *t* (*n*) *s* .

Braking phases are accompanied by intense frictional heating of the friction pair elements (disc-pad, drum-shoe, etc.). Because of high temperature gradients, the friction coefficient and thermo-mechanical properties of friction pair elements can differ significantly at the initial phase and after braking. Hence, the assumption is made that, during the *k*th braking, friction coefficient *f* , thermal conductivities *Kl*, specific heat capacities *cl*, densities *ρ<sup>l</sup>* and Brinell hardness *HBl* of friction pair materials *l* = 1, 2 depend on temperature *T* in the form of, Equations (1) and (2):

$$f(T) = f\_0 f^\*(T). \tag{1}$$

$$K\_l(T) = K\_{l,0} K\_l^\*(T), \\ \mathcal{c}\_l(T) = \mathcal{c}\_{l,0} \mathcal{c}\_l^\*(T), \\ \rho\_l(T) = \rho\_{l,0} \rho\_l^\*(T), \\ Hb\_l(T) = HB\_{l,0} B\_l^\*(T), \tag{2}$$

Where Equations (3)–(8),

$$f\_0 = f(T\_0),\ K\_{l,0} = K\_l(T\_0),\ c\_{l,0} = c\_l(T\_0),\ \rho\_{l,0} = \rho\_l(T\_0),\ HB\_{l,0} = HB\_l(T\_0),\tag{3}$$

$$f^\*(T) = f\_1 + \frac{f\_2}{[f\_3(T - f\_4)]^2 + 1} + \frac{f\_5}{[f\_6(T - f\_7)]^2 + 1},\tag{4}$$

$$K\_l^\*(T) = K\_{l,1} + \frac{K\_{l,2}}{[K\_{l,3}(T - K\_{l,4})]^2 + 1} + \frac{K\_{l,5}}{[K\_{l,6}(T - K\_{l,7})]^2 + 1},\tag{5}$$

$$c\_l^\*(T) = c\_{l,1} + \frac{c\_{l,2}}{[c\_{l,3}(T - c\_{l,4})]^2 + 1} + \frac{c\_{l,5}}{[c\_{l,6}(T - c\_{l,7})]^2 + 1},\tag{6}$$

$$\rho\_l^\*(T) = \rho\_{l,1} + \frac{\rho\_{l,2}}{\left[\rho\_{l,3}(T-\rho\_{l,4})\right]^2 + 1} + \frac{\rho\_{l,5}}{\left[\rho\_{l,6}(T-\rho\_{l,7})\right]^2 + 1},\tag{7}$$

$$HB^\*(T) = HB\_1 + \frac{HB\_2}{\left[HB\_3(T - HB\_4)\right]^2 + 1} + \frac{HB\_5}{\left[HB\_6(T - HB\_7)\right]^2 + 1},\tag{8}$$

where *T*0—initial temperature, *fj*, *Kl*,*j*, *cl*,*j*, *HBl*,*j*, *ρl*,*j*, *l* = 1, 2; *j* = 1, 2, ... , 7 are parameters of the experimental data approximations [26]. Here and below, subscripts *l* = 1 and *l* = 2 refer to quantities associated with the primary element (disc, drum, etc.) and friction lining (pad, shoe, etc.), respectively.

The remaining assumptions of the calculation model are as follows:

1. At the initial moment of each braking phase, the friction element is pressed against the primary element contact surface with uniform pressure *p*, which exponentially increases with time *t*, from zero to nominal value *p*0, Equation (9):

$$p(t) = p\_0 p^\*(t), \; p^\*(t) = 1 - e^{-\frac{t}{t\_i}}, \; 0 \le t \le t\_s^{(k)}, k = 1, 2, \dots, n,\tag{9}$$

where *ti*—time of pressure increase.


The applied contact pressure (10) causes a reduction of sliding velocity *V*(*k*) during the *k*th braking, according to the following relation [25], Equation (10):

$$V^{(k)}(t) = V\_0 \, V^{\*(k)}(t), \, V^{\*(k)}(t) = 1 - \frac{t}{t\_{s,0}^{(k)}} + p^\*(t) \frac{t\_i}{t\_{s,0}^{(k)}}, \, 0 \le t \le t\_s^{(k)}, \, k = 1, 2, \ldots, n,\tag{10}$$

Where, Equation (11),

$$\eta\_{s,0}^{(k)} = \frac{2\mathcal{W}\_0}{q\_0^{(k)}A\_a} \,'\tag{11}$$

*W*0—initial kinetic energy of the system, which depends on the mass of the vehicle and its velocity at the moment just before the brake applications, *q* (*k*) <sup>0</sup> —nominal specific friction power, *Aa*—nominal contact surface area, usually determined by the dimensions of the contact surface of friction element. Braking time *t* (*k*) *<sup>s</sup>* is obtained from the stop condition, Equation (12):

*t*

$$V^{\*(k)}(t\_s^{(k)}) = 0.\tag{12}$$

In case of the immediate (*ti* → 0) achievement of nominal pressure value *p*<sup>0</sup> in the Equation (9) from Equation (10) we obtain the linear drop in speed, Equation (13):

$$V^{\*(k)}(t) = 1 - \frac{t}{t\_{s,0}^{(k)}}, \; 0 \le t \le t\_{s,0}^{(k)}.\tag{13}$$

Thus, parameter *t* (*k*) *<sup>s</sup>*,0 (11) is the time of braking with constant deceleration. In the phases of acceleration, the velocity increases linearly with time, Equation (14):

$$V^{\*(k)}(t) = \frac{t - t\_s^{(k)}}{t\_c}, \ t\_s^{(k)} \le t \le t\_{k\prime} \ k = 1, 2, \dots, n. \tag{14}$$

The evolution of maximum temperature *<sup>T</sup>*(*k*) max on the friction surface during *<sup>k</sup>*th braking was sought in the form of [9,15]:

$$T\_{\text{max}}^{(k)}(t) = T\_m^{(k)}(t) + T\_f^{(k)}(t),\ 0 \le t \le t\_s^{(k)},\ k = 1, 2, \dots, n\_\prime. \tag{15}$$

where *<sup>T</sup>*(*k*) *<sup>m</sup>* —mean temperature of the nominal contact surface and *<sup>T</sup>*(*k*) *<sup>f</sup>* —average temperature of the real contact area (flash temperature).

### **3. Solution to the Problem**

### *3.1. Heat Generation on the Nominal Contact Surface*

Taking into account the above assumptions, to find component *<sup>T</sup>*(*k*) *<sup>m</sup>* in the relation (16), we used the a calculation scheme of frictional contact of two different semi-infinite bodies, <sup>0</sup> <sup>≤</sup> *<sup>z</sup>* <sup>&</sup>lt; <sup>∞</sup> (*<sup>l</sup>* <sup>=</sup> 1) and <sup>−</sup><sup>∞</sup> <sup>&</sup>lt; *<sup>z</sup>* <sup>≤</sup> 0 (*<sup>l</sup>* <sup>=</sup> 2), sliding against each other with velocity *<sup>V</sup>*(*k*)(*t*) (10)–(13). Initiated by frictional heating during *k*th braking, we found transient temperature field *T*(*k*)(*z*, *t*) in this system from the solution to the following boundary-value problem of heat conduction, Equations (16)–(20):

$$\frac{\partial^2 T^{(k)}(z,t)}{\partial z^2} = \frac{1}{k\_{1,0}^{(k)}} \frac{\partial T^{(k)}(z,t)}{\partial t}, 0 < z < \infty, 0 < t \le t\_s^{(k)} \,. \tag{16}$$

$$\left. K\_{2,0}^{(k)} \frac{\partial T^{(k)}(z,t)}{\partial z} \right|\_{z=0^-} - \left. K\_{1,0}^{(k)} \frac{\partial T^{(k)}(z,t)}{\partial z} \right|\_{z=0^+} = q^{(k)}(t), 0 < t \le t\_s^{(k)} \tag{17}$$

$$T^{(k)}\left(0^+,t\right) = T^{(k)}\left(0^-,t\right) \equiv T^{(k)}\_{m}\left(t\right), \; 0 < t \le t\_s^{(k)}\;\tag{18}$$

$$T^{(k)}(z,t) \to T\_0^{(k)} \; , \; |z| \to \infty \; , \; 0 < t \le t\_s^{(k)} \tag{19}$$

$$T^{(k)}(z,0) = T\_0^{(k)}, \quad |z| < \infty, \ k = 1, 2, \dots, n \tag{20}$$

where, Equations (21)–(23),

$$q^{(k)}(t) = q\_0^{(k)} q^{\*(k)}(t),\ q\_0^{(k)} = f\_0^{(k)} p\_0 V\_{0\prime} \ q^{\*(k)}(t) = p^\*(t) V^{\*(k)}(t) \tag{21}$$

$$k\_{l,0}^{(k)} = \frac{K\_{l,0}^{(k)}}{\rho\_{l,0}^{(k)} c\_{l,0}^{(k)}}, l = 1, 2,\tag{22}$$

$$f\_0^{(k)} = f(T\_0^{(k)}),\ K\_{l,0}^{(k)} = \mathcal{K}\_l(T\_0^{(k)}),\ \ c\_{l,0}^{(k)} = \mathcal{c}\_l(T\_0^{(k)}),\ \rho\_{l,0}^{(k)} = \rho\_l(T\_0^{(k)}).\tag{23}$$

Averaged volumetric temperature *T*(*k*) <sup>0</sup> of a braking system before *k*th braking was calculated from the formulas, Equations (24)–(26) [9,16]:

$$T\_0^{(k)} = \frac{T\_{0,0}^{(k)} + T\_{0,1}^{(k)}}{2}, k = 1, 2, \dots, n\_\prime. \tag{24}$$

$$T\_{0,i}^{(k)} = T\_0 + \frac{\gamma\_i W\_0}{2G} \left( \frac{e^{-a\_i t\_\mathcal{c}} - e^{-ka\_i t\_\mathcal{c}}}{1 - e^{-a\_i t\_\mathcal{c}}} \right), i = 0, 1,\tag{25}$$

$$\gamma\_i = \frac{\sqrt{K\_{1,i} \rho\_{1,i} c\_{1,i}}}{\sqrt{K\_{1,i} \rho\_{1,i} c\_{1,i}} + \sqrt{K\_{2,i} \rho\_{2,i} c\_{2,i}}} \tag{26}$$

where, Equations (27)–(29),

$$\gamma\_{i} = \frac{\sqrt{\mathcal{K}\_{1,i} \rho\_{1,i} c\_{1,i}}}{\sqrt{\mathcal{K}\_{1,i} \rho\_{1,i} c\_{1,i}} + \sqrt{\mathcal{K}\_{2,i} \rho\_{2,i} c\_{2,i}}} \tag{27}$$

$$\alpha\_{i} = \frac{hA\_{\text{event}}}{\text{G }c\_{1,i}} \tag{28}$$

$$K\_{l,1} = K\_l(T\_{0,0}^{(k)}),\ c\_{l,1} = c\_l(T\_{0,0}^{(k)}),\ \rho\_{l,1} = \rho\_l(T\_{0,0}^{(k)}),\tag{29}$$

*h*—heat transfer coefficient, *G*, *Avent*—weight and total surface area of the primary element of friction pair, respectively; values *Kl*,0, *cl*,0, *ρl*,0, *l* = 1, 2 are defined by Equation (3). For the first braking cycle (*k* = 1)), from relations (25) and (26), it follows that *T*(1) <sup>0</sup> = *T*0.

Exact solution to the boundary-value problem of heat conduction, (17)–(22), during single braking was obtained in [27]. Generalizing this solution to the considered case of repetitive braking, the sought temperature on the friction surface was written in the form, Equation (30):

$$T\_m^{(k)}(t) = T\_0^{(k)} + T\_{m,0}^{(k)} T\_m^{\*(k)}(t), 0 \le t \le t\_s^{(k)}, k = 1, 2, \dots, n \tag{30}$$

where, Equations (31)–(37),

$$T\_m^{\*(k)}(t) = \gamma\_0^{(k)} \sqrt{\frac{t}{t\_{s,0}^{(k)}}} \left[ \left( 1 + \frac{t\_i}{t\_{s,0}^{(k)}} - \frac{2t}{3t\_{s,0}^{(k)}} \right) \frac{2}{\sqrt{\pi}} - \left( 1 + \frac{3t\_i}{2t\_{s,0}^{(k)}} - \frac{t}{t\_{s,0}^{(k)}} \right) \mathcal{F} \left( \sqrt{\frac{t}{t\_{s,0}^{(k)}}} \right) + \frac{t\_i}{t\_{s,0}^{(k)}} \mathcal{F} \left( \sqrt{\frac{2t}{t\_{s,0}^{(k)}}} \right) \right] \tag{31}$$

$$\mathcal{F}(\mathbf{x}) = \frac{2}{\sqrt{\pi}} \sum\_{n=0}^{\infty} \left(-1\right)^{n} \frac{\left(2\mathbf{x}^{2}\right)^{n}}{\left(2n+1\right)!!}, \ 0 \le \mathbf{x} \le \mathbf{3}, \ \mathcal{F}(\mathbf{x}) = \frac{2}{\sqrt{\pi}} \sum\_{n=0}^{\infty} \frac{\left(2n-1\right)!!}{\left(2\mathbf{x}^{2}\right)^{n+1}}, \ \mathbf{x} > \mathbf{3}, \tag{32}$$

$$\gamma\_0^{(k)} = \frac{\sqrt{K\_{1,0}^{(k)} \rho\_{1,0}^{(k)} c\_{1,0}^{(k)}}}{\sqrt{K\_{1,0}^{(k)} \rho\_{1,0}^{(k)} c\_{1,0}^{(k)}} + \sqrt{K\_{2,0}^{(k)} \rho\_{2,0}^{(k)} c\_{2,0}^{(k)}}},\tag{33}$$

$$T\_{m,0}^{(k)} = \frac{q\_0^{(k)} a^{(k)}}{K\_{1,0}^{(k)}} \tag{34}$$

$$a^{(k)} = \max\left\{ a\_l^{(k)} \right\} \tag{35}$$

$$a\_l^{(k)} = \begin{cases} \quad d\_{l\prime} \quad a\_{l,eff\prime}^{(k)} \ge d\_{l\prime},\\ \quad a\_{l,eff\prime}^{(k)} \quad a\_{l,eff\prime}^{(k)} < d\_{l\prime} \end{cases} \tag{36}$$

$$a\_{l,eff}^{(k)} = \sqrt{\mathfrak{B}\_{l,0}^{(k)} t\_{s,0}^{(k)}}\,\tag{37}$$

*dl*, *l* = 1, 2 thicknesses of the friction pair components.

The algorithm to find the mean temperature for the selected friction pair consisted in the subsequent performance of the following steps:


### *3.2. Temperature of the Real Contact Region*

Real friction surfaces of braking system elements are not perfectly smooth, they are characterized by significant roughness and waviness. Therefore, these elements are not in contact on an entire nominal friction area, but only in some parts of one, and consists of roughness waves. These waves form the so-called contour surface of contact, the area of which (*A*(*k*) *<sup>c</sup>* ) changes during *k*th braking, according to the law, Equations (38) and (39) [28]:

$$A\_c^{(k)}(t) = A\_d \left[ \frac{p(t)b\_0^{\nu-1}}{HB\_{\text{min}}^{(k)}(t)} \right]^{\frac{1}{\nu+1}}, \ 0 \le t \le t\_s^{(k)},\tag{38}$$

$$HB\_{\min}^{(k)}(t) = \min \left\{ HB\_1[T\_m^{(k)}(t)], HB\_2[T\_m^{(k)}(t)] \right\} \tag{39}$$

where *<sup>p</sup>*, *<sup>T</sup>*(*k*) *<sup>m</sup>* —pressure (9) and mean temperature (30)–(37) on the nominal contact surface with area *Aa*, *b*0, *ν*—parameters of the reference curve for the harder element of the friction pair materials. Usually this is the material of a primary element (disc, drum etc.). Pressure *p* (*k*) *<sup>c</sup>* , on contour area *<sup>A</sup>*(*k*) *<sup>c</sup>* (38), (39) was found from the relation, Equation (40):

$$p\_c^{(k)}(t) = p(t) \frac{A\_d}{A\_c^{(k)}(t)}, \ 0 \le t \le t\_s^{(k)} \tag{40}$$

In a typical braking system, the friction linings are made of more deformable and less durable material than the primary element material. Thus, we assume that:


3. Before coming into contact with the friction lining, the temperature of asperity does not change along its height and is equal to the mean temperature of nominal contact area *<sup>T</sup>*(*k*) *<sup>m</sup>* .

Contact of a single asperity with a smooth friction lining surface creates the real region of contact with a diameter as follows Equation (41) [29]:

$$d\_r^{(k)}(t) = \left(\frac{8r\_{av}h\_{\max}}{\nu}\right)^{\frac{1}{2}} \left[\frac{p\_c^{(k)}(t)}{HB\_{\min}^{(k)}(t)b\_0}\right]^{\frac{1}{2v}}, \ 0 \le t \le t\_s^{(k)}\tag{41}$$

where *rav*—averaged radius of the asperities rounding and *h*max—maximum height of the roughness. Total area *A*(*k*) *<sup>r</sup>* of the real contact region is related to area *Aa* of the nominal contact surface according to the formula, Equation (42):

$$A\_r^{(k)}(t) = A\_a \frac{p(t)}{HB\_{\text{min}}^{(k)}(t)}, \ 0 \le t \le t\_s^{(k)}.\tag{42}$$

Evolution of temperature *T*(*k*) *<sup>f</sup>* on the real contact region with area *<sup>A</sup>*(*k*) *<sup>r</sup>* (42) (flash temperature) was determined from the following Equation (43) [26,30]:

$$T\_f^{(k)}(t) = \left(1 + \frac{1}{\sqrt{2}}\right) \frac{f\_m^{(k)}(t)p(t)V^{(k)}(t)A\_d d\_r^{(k)}(t)}{A\_r^{(k)}(t)[4K\_{1,m}^{(k)}(t) + \sqrt{\pi V^{(k)}(t)\ d\_r^{(k)}(t)}K\_{2,m}^{(k)}(t)\ c\_{2,m}^{(k)}(t)\ \rho\_{2,m}^{(k)}(t)]}, \ 0 \le t \le t\_s^{(k)}, \tag{43}$$

Where Equation (44):

$$f\_{\mathfrak{m}}^{(k)}(t) \equiv f\left[T\_{\mathfrak{m}}^{(k)}(t)\right],\ \mathbf{K}\_{l,\mathfrak{m}}^{(k)}(t) \equiv \mathbf{K}\_{l}\left[T\_{\mathfrak{m}}^{(k)}(t)\right],\ c\_{l,\mathfrak{m}}^{(k)}(t) \equiv c\_{l}\left[T\_{\mathfrak{m}}^{(k)}(t)\right],\ \rho\_{l,\mathfrak{m}}^{(k)}(t) \equiv \rho\_{l}\left[T\_{\mathfrak{m}}^{(k)}(t)\right],\ l=1,2.\tag{44}$$

Summarizing, another scheme was proposed for determining the flash temperature evolution during *k*th braking:


### **4. Numerical Analysis**

Based on the analytical model proposed above, the temperature generated during RST braking mode was studied for two tribosystems: disc and drum brakes. The temperature field of the disc brake system was analyzed using the finite differences method in [16]. Additionally, the corresponding evolutions of temperature in the drum brake were presented in [9], where the calculations were performed on the basis of approximate solutions to the thermal problem of friction for two strips and experimental data. Results of calculations presented below were obtained for the same friction materials and operational parameters as the above-mentioned studies. In both systems, the RST braking mode was considered, which includes 4 cycles (*n* = 4): three full cycles (braking-acceleration) and the last interrupted cycle (braking). The brake disc is made of cast iron ChNMKh,

and brake pads—were manufactured from cermet FMC-11. In the second system, the drum brake was made of 30KhHSA steel and the brake shoes were manufactured with retinax FC-16L. The sintered cermet friction material FMC-11 contains 64% Fe, 15% Cu, 3% SiO2, 6% BaSO4, 3% asbestos and 9% graphite. Retinax FC-16L is a composite based on phenol-formaldehyde resins and reinforced with brass shavings [21]. The thermophysical and mechanical properties of materials and friction coefficients of selected pairs at the initial temperature are included in Table 1. The methodology of calculations adopted in the present study is shown in Figure 2.

**Table 1.** Properties of friction pair materials at the initial temperature *T*<sup>0</sup> = 20 ◦C [26].


**Figure 2.** Flowchart of the calculation procedure.

Experimental data of frictional thermal stability and property variations under different temperature conditions of friction materials ChNMKh/FMC-11 and 30KhHSA/FC-16L are included in [9,26]. The values of coefficients in the functions (4)–(8), which approximate these data, are presented in Table 2. In this paper, it was assumed that density *ρ<sup>l</sup>* of all selected materials changes slightly under temperature variations, so Equation (7) takes the form *ρ*∗ *<sup>l</sup>* (*T*) = 1, *l* = 1, 2. Graphs of functions *f*(*T*) (4), *Kl*(*T*) (5), *cl*(*T*) (6) and *HBl*(*T*), *l* = 1, 2 (8) are presented in Figure 3.


**Table 2.** Coefficients in approximation functions (4)–(6) and (8) for considered materials [31].

**Figure 3.** Graphs of functions approximating the experimental data of the dependence on temperature: (**a**) coefficient of friction *f* ; (**b**) thermal conductivity *Kl*; (**c**) specific heat capacity *cl*; (**d**) Brinell hardness *HBl*, *l* = 1, 2.

### *4.1. Disc Brake System*

The input values of operating parameters [15,16]: *p*<sup>0</sup> = 1.47 MPa, *V*<sup>0</sup> = 27.78 m s<sup>−</sup>1, *<sup>W</sup>*<sup>0</sup> <sup>=</sup> 392.1 kJ, *Aa* <sup>=</sup> 4.047 · <sup>10</sup>−<sup>2</sup> m2, *Avent* <sup>=</sup> 4.44 · <sup>10</sup>−<sup>2</sup> m2, *<sup>G</sup>* <sup>=</sup> 1.58 kg, *<sup>h</sup>* <sup>=</sup> <sup>100</sup> Wm<sup>−</sup>2K−1, *d*<sup>1</sup> = 5.5 mm, *d*<sup>2</sup> = 10 mm, *ti* = 0.5 s, *tc* = 5 s, *rav* = 450 μm, *h*max = 2.5 μm, *b*<sup>0</sup> = 1.0, *ν* = 2.1. Due to the geometric symmetry of the disc brake system with respect to the center plane of the disc, to determine the temperature field we consider the tribosystem, which consists of disc replacement and the pad with thicknesses of *dl*, *l* = 1, 2, respectively. The thickness of the disc replacement element *d*<sup>1</sup> is equal to the half of the real disc dimension.

The temporal profiles of velocity *V*(*k*) (10) and specific friction power *q*(*k*) (22) during four considered braking cycles are presented in Figure 4. Except for the short initial period of braking, when the pressure grows, the sliding speed linearly decreases with time and the stop time is longer with every subsequent cycle (Figure 4a). The specific friction power increases at the beginning of each braking cycle until reaching the maximum value, and then decreases to zero at the stop moment (Figure 4b), which is characteristic for the socalled rational braking mode [18,19]. It should be noted that the friction work performed at each braking cycle (the area under each of the four curves in Figure 4b) was the same and equal to the initial kinetic energy of the system *W*<sup>0</sup> = 392.1 kJ.

**Figure 4.** Changes in: (**a**) velocity *V*(*k*) and (**b**) specific friction power *q*(*k*) during *k* = 1, 2, 3, 4 braking applications of a disc brake system.

The coefficient of friction *f* (*k*) <sup>0</sup> (24) decreases with every subsequent braking from 0.45 during the first cycle to 0.28 during the fourth, last braking (Figure 5a). The drop of the friction coefficient value causes the elongation of braking stage *t* (*k*) *<sup>s</sup>* —its values are equal to 1.54 s, 1.73 s, 1.96 s and 2.22 s, respectively for *k* = 1, 2, 3, 4 (Figure 5b).

Variations of mean temperature *<sup>T</sup>*(*k*) *<sup>m</sup>* (30), flash temperature *<sup>T</sup>*(*k*) *<sup>f</sup>* (43) and maximum temperature *<sup>T</sup>*(*k*) max (15) during each braking cycle of the process are shown in Figure 6.

Evolutions of mean temperature on the nominal contact area during every braking stage correspond to the temporal profiles of specific friction power, presented in Figure 4b. At the beginning of each braking cycle, the mean temperature *<sup>T</sup>*(*k*) *<sup>m</sup>* growth takes place until maximum value is reached, and then the temperature drop to the standstill occurs. According to the adapted plastic roughness deformation mechanism, the flash temperature reaches significant values at the initial moments of the braking stages, when the contact area of the pad with the disc is relatively "cold". Subsequent heating of this area causes a rapid decrease in the flash temperature *T*(*k*) *<sup>f</sup>* . The mean temperature of the friction surface has a decisive influence on the time profile and the maximum temperature value. With each subsequent braking, the highest values of *<sup>T</sup>*(*k*) *<sup>m</sup>* and *<sup>T</sup>*(*k*) max increase, while the *<sup>T</sup>*(*k*) *<sup>f</sup>* decreases.

**Figure 5.** The values of: (**a**) friction coefficient *f* (*k*) <sup>0</sup> ; (**b**) stop time *t* (*k*) *<sup>s</sup>* for disc brake in each braking cycle *k* = 1, 2, 3, 4.

**Figure 6.** Evolutions of the mean *<sup>T</sup>*(*k*) *<sup>m</sup>* , flash *<sup>T</sup>*(*k*) *<sup>f</sup>* and maximum *<sup>T</sup>*(*k*) max temperatures for disc brake during each cycle: (**a**) *k* = 1; (**b**) *k* = 2; (**c**) *k* = 3; (**d**) *k* = 4.

Temporal profiles of *<sup>T</sup>*(*k*) *<sup>m</sup>* , *<sup>T</sup>*(*k*) *<sup>f</sup>* and *<sup>T</sup>*(*k*) max during the entire RST braking mode are presented in Figure 7. The combined time of all four braking cycles duration is 7.45 s and with account of the three stages of acceleration 15 s, gives the total time of the whole RST braking mode in the disc brake system *tb* = 22.45 s.

**Figure 7.** Variations of the maximum *<sup>T</sup>*(*k*) max (solid lines), the mean *<sup>T</sup>*(*k*) *<sup>m</sup>* (dashed lines) and the flash *<sup>T</sup>*(*k*) *<sup>f</sup>* (dotted lines) temperatures during disc brake RST mode.

### *4.2. Drum Brake System*

The values of input parameters of braking in the RST mode by friction pair consisting of the drum (*l* = 1) and brake shoe (*l* = 2), were adapted [9]: *p*<sup>0</sup> = 0.44 MPa, *<sup>V</sup>*<sup>0</sup> = 11.5 m s<sup>−</sup>1, *<sup>W</sup>*<sup>0</sup> = 215.7 kJ, *Aa* = *Ac* = 3.85 · <sup>10</sup>−<sup>2</sup> m2, *Avent* = <sup>15</sup> · <sup>10</sup>−<sup>2</sup> <sup>m</sup>2, *G* = 5.5 kg, *h* = 80 Wm<sup>−</sup>2K−1, *d*<sup>1</sup> = 10 mm, *d*<sup>2</sup> = 18 mm, *ti* = 0.5 s, *tc* = 25 s, *rav* = 500 μm, *h*max = 4.5 μm, *b*<sup>0</sup> = 1.0, *ν* = 2.2.

The evolutions of sliding velocity and specific friction power for the drum brake are presented in Figure 8. There are some similarities with the corresponding profiles for a disc brake, which are shown in Figure 4. Sliding velocities also have short time of nonlinear drop (Figure 8a), and the profiles of specific friction power have a local maximum, shifted closer to the beginning of the braking period in this case (Figure 8b). Unlike the disc system, however, the duration of the braking stage in the drum brake is shortened and the maximum value of the specific friction power decreases with each subsequent braking.

In Figure 8a the shortening of the time of subsequent braking can be noticed, which is related to the corresponding changes in the friction coefficient: in the considered drum system its value during the first braking is 0.39 and increases linearly to 0.42 during the fourth braking (Figure 9a). The respective stop times are equal to 6.17 s and 5.77 s (Figure 9b).

Much smaller changes in the specific friction power time profiles during each braking (Figure 8b) than in the case of a disc brake result in the fact that the temperature evolution in the drum brake during each braking stage also differs much less (Figure 10). On the other hand, the contribution of the mean temperature *<sup>T</sup>*(*k*) *<sup>m</sup>* and the flash temperature *<sup>T</sup>*(*k*) *f*

to the value of the maximum temperature *<sup>T</sup>*(*k*) max in the drum brake are completely different than in the disc brake. In the latter, as noted above, the time profile of the maximum temperature was primarily shaped by the mean temperature on the friction surface. In the drum brake, both components *<sup>T</sup>*(*k*) *<sup>m</sup>* and *<sup>T</sup>*(*k*) *<sup>f</sup>* show a significant influence on evolution and

values of *<sup>T</sup>*(*k*) max. The decisive factor in the initial braking phase is the flash temperature, and in the final stage—the mean temperature of the drum-brake shoe contact area.

**Figure 8.** Changes in time during RST mode of the drum brake: (**a**) velocity *V*(*k*); (**b**) specific friction power *q*(*k*), *k* = 1, 2, 3, 4.

**Figure 9.** The values of: (**a**) coefficient of friction *f* (*k*); (**b**) stopping time *t* (*k*) *<sup>s</sup>* for the drum brake in each braking cycle *k* = 1, 2, 3, 4.

The changes in the temperature with time during RST braking in the drum brake are presented in Figure 11. Duration of four braking cycles was equal to 23.81 s, three accelerations 75 s, so the total RST operation mode lasted 98.81 s. As in the case of the disc brake, with each subsequent braking, the highest values of the mean and the maximum temperatures increase, whereas the flash temperatures decrease. However, these relations are not as noticeable as in the disc brake. One of the reasons for such relatively small changes in temperature can result from much longer (five times) cooling phase during vehicle acceleration.

**Figure 10.** Evolutions of the mean *<sup>T</sup>*(*k*) *<sup>m</sup>* , flash *<sup>T</sup>*(*k*) *<sup>f</sup>* and maximum *<sup>T</sup>*(*k*) max temperatures for drum brake during each cycle: (**a**) *k* = 1; (**b**) *k* = 2; (**c**) *k* = 3; (**d**) *k* = 4.

**Figure 11.** Evolutions of the maximum *<sup>T</sup>*(*k*) max (solid lines), mean *<sup>T</sup>*(*k*) *<sup>m</sup>* (dashed lines) and flash *<sup>T</sup>*(*k*) *f* (dotted lines) temperatures during drum brake RST mode.

Values of the friction coefficient *f* (*k*), braking time *t* (*k*) *<sup>s</sup>* , volumetric temperature *<sup>T</sup>*(*k*) 0 and the highest values of the mean temperature *<sup>T</sup>*(*k*) *<sup>m</sup>* , flash temperature *<sup>T</sup>*(*k*) *<sup>f</sup>* and maximum temperature *<sup>T</sup>*(*k*) max during performing each of the four cycles are included in Table 3.


**Table 3.** Calculated values of some characteristics at each braking cycle.

### **5. Summary of the Results and Discussion**

The analytical model to determine the maximum temperature reached as a result of friction during the repetitive short-term operating mode of the braking systems, consisting of the *n* braking and acceleration cycles was proposed. The thermal sensitivity of the materials of the friction pairs and the temperature dependence of the friction coefficient were taken into account. According to the Chichinadze's summation hypothesis, the maximum temperature *<sup>T</sup>*(*k*) max during each braking (*<sup>k</sup>* <sup>=</sup> 1, 2, ... , *<sup>n</sup>*) was searched for as the sum of the mean temperature *<sup>T</sup>*(*k*) *<sup>m</sup>* on the nominal contact area and the average temperature *T*(*k*) *<sup>f</sup>* of the real contact area (flash temperature). The numerical analysis was carried out for two systems: disc brake (cast iron/cermet) and drum brake (steel/retinax) with *n* = 4:


◦C, respectively. The greatest relative percentage difference of the results occurred in the fourth stage and is equal to ≈ 11%. In the drum brake the maximum temperatures *<sup>T</sup>*(*k*) max determined by means of the proposed model, are equal to 353 ◦C, 387 ◦C, 416 ◦C and 443 ◦C (Table 3), and corresponding results presented in monograph [9] are 295 ◦C, 330 ◦C, 400 ◦C and 440 ◦C. The highest relative difference in outcomes occurs in the stage one and is equal to ≈ 21%. It should be noted that the mean temperature *<sup>T</sup>*(*k*) *<sup>m</sup>* of the selected friction pair drum-brake shoe with similar input parameters during single braking (*k* = 1) was analyzed with the use of the finite difference method in the article [32]. The highest value of the mean temperature on the nominal contact area of the drum-brake shoe obtained was equal to 210 ◦C [32], which is in good agreement with the value 280 ◦C presented in Table 3.

### **6. Conclusions**

The main advantage of the proposed approach is the demonstration of the possibility of adapting numerous exact solutions to the linear thermal problems of friction existing in the scientific literature, to determine the maximum temperature of the brake not only during a single, but also during a repetitive operation mode. This allows for express estimation with sufficient accuracy, not only of the maximum temperature, but also of important braking characteristics such as variations of the sliding velocity and specific friction power, stopping time and braking distance during each stage of RST braking mode. It should be emphasized that the proposed calculation model is coupled—it allows to take into account the mutual influence of all of the above-mentioned characteristics in the braking process, by the temperature-dependent coefficient of friction. From the point of view of application possibilities, it is important that this model can be used for calculations not only for materials with stable thermo-physical properties, but also for thermally sensitive materials.

**Author Contributions:** Conceptualization and methodology, A.Y and M.K.; software, K.T.; validation, A.Y.; formal analysis, K.T.; investigation, A.Y. and K.T.; writing—original draft preparation, A.Y. and K.T.; writing–review and editing, A.Y. and K.T.; visualization and figures preparation, K.T.; supervision, A.Y. and M.K.; project administration, M.K.; funding acquisition, M.K. All authors have read and agreed to the published version of the manuscript.

**Funding:** This study was performed within the framework of research project No. 2017/27/B/ST8/ 01249, funded by the National Science Centre, Poland and project financing 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 financing 12,000,000 PLN.

**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**


### **References**


*Article*
