Next Article in Journal
Characterization of Co-Cr-W Dental Alloys with Veneering Materials Manufactured via Subtractive Milling and Additive Manufacturing LDED Methods
Next Article in Special Issue
Influence of the Additive of Ceramic and Intermetallic Powders on the Friction Properties and Temperature of the Wet Clutch Disc
Previous Article in Journal
Fundamental Studies on Crystallization and Reaching the Equilibrium Shape in Basic Ammonothermal Method: Growth on a Native Lenticular Seed
Previous Article in Special Issue
Temperature in the Friction Couple Consisting of Functionally Graded and Homogeneous Materials
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Heat Partition Ratio during Braking in a Functionally Graded Friction Couple

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: 1 June 2022 / Revised: 23 June 2022 / Accepted: 28 June 2022 / Published: 30 June 2022

Abstract

:
The theoretical scheme for determining the heat partition ratio (HPR) in a friction couple made of functionally graded materials (FGMs) was proposed. As a result, the formula for the calculation of the HPR was found, which depends on the thermal properties and the parameters of the material’s gradient. In specific cases of these parameters, the known formulas for estimating the HPR for homogeneous materials were obtained. Calculations were carried out for the friction couple consisting of the following two-component FGMs: Al2O3–Cu (first body) and ZrO2–Ti–6Al–4V (second body), under the conditions corresponding to a single braking with a constant deceleration. It was established that the vast majority (almost 90%) of heat that was generated by friction was absorbed by the first body in the selected couple. The possibilities of using the obtained results were discussed herein.

1. Introduction

Functionally graded materials (FGMs) are materials in which, along at least one specific direction, a continuous change in the functional or structural properties has been obtained in a selected technological process. The possibilities of using graded materials in technology seem to be almost unlimited. Proposals of possible FGM applications appear for many industries, such as the following: optical, energy [1,2], aviation [3,4], medical [5], etc. Moreover, they have been used in friction elements, including braking systems [6,7].
The significant influence of temperature on the friction and wear characteristics of friction materials is well known and has been studied by many researchers who are involved in determining the temperature mode of disc brakes [8,9,10]. The basis for establishing the temperature are the solutions (analytical, numerical, or analytical–numerical) to the thermal problems of friction, i.e., the boundary-value problems of heat conduction with two specific boundary conditions on the contact surface of the pad and the disc [11,12]. One of them, the energetic solution, determines the equality of the specific friction power q and the sum of the heat flux intensities q l , l = 1 , 2 that are directed along the normal surface to the contact surface towards the insides of the friction elements q 1 + q 2 = q . Whereas the second solution concerns the type of thermal contact—with (imperfect) or without (perfect) consideration of the thermal resistance on the friction surfaces. In the latter case, the surface temperatures of the friction elements are the same. The coupling of the temperature fields of both of the friction elements through the above-mentioned boundary conditions means that obtaining the analytical [13,14,15,16] or the numerical [17,18,19] solutions of the thermal problems of friction requires the performance of complex mathematical transformations or long-term calculations.
In order to reduce these difficulties, there is also another approach to estimate the temperature of the braking systems, on the basis of solutions to the thermal problems of friction. It is based on a virtual separation of the friction pair elements and the subsequent replacement of the frictional interaction on their working surfaces with the heat fluxes of the following intensities: q l = α l q , l = 1 , 2 , α 1 = α , and α 2 = 1 α , respectively, where α is the heat partition ratio [20,21]. Reviews of the experimental methods, as well as the theoretical methods, for determining the HPR in braking systems have been developed [22,23,24,25]. The theoretical methods primarily rely on the fact that the temperature field, which is first found by means of analytical solutions, contains a priori unknown coefficient α , which is then determined from the condition of the equality of the maximum or the mean temperature on the friction surfaces of the pad and the disc. As a result, the formulas for the calculation of α contain the thermo-physical properties of the friction materials and, in some cases, the operating parameters of the process. Substituting thus determines the value of α back into the analytical or the numerical solutions, and the temperature mode of the braking system is estimated.
It should be noted that the formulas that have been obtained so far for determining the HPR, based on the analytical solutions to the thermal problems of friction, concern homogeneous materials [26,27,28,29,30]. The proportion of the heat distribution by means of the HPR between polymer and steel elements was established experimentally [31]. The aim of this study is to obtain the theoretical formulas for estimating the HPR during the braking of friction pair elements that are made of functionally graded materials (FGMs). One of the proposed applications of the results could be the determination of the bulk temperature of friction elements that are made of thermally sensitive FGM, during a single braking process [32] or a repetitive short-term braking mode [33].

2. Heating of the FGM Semi-Space by the Heat Flux with Constant Intensity

Consider a temperature of T of a semi-limited body z 0 , which is uniform and equal to T 0 at the initial point in time of t = 0 . The body is made of a functionally graded material (FGM) with the thermal conductivity of K , increasing exponentially in the positive direction of the axis z as follows:
K ( z ) = K 0 e γ z , z 0 , K 0 K ( 0 ) ,
where γ > 0 is the gradient of the material. Next, the semi-space is heated on the outer surface as z = 0 by the heat flux with a constant intensity of q = q 0 in time as t > 0 (Figure 1).
The transient temperature field of the semi-space was searched in following the form:
T ( z , t ) = T 0 + Θ ( z , t ) ,   z 0 ,   t 0 ,
where the temperature rise of Θ ( z , t ) was determined from the solution to the following boundary-value heat conduction problem:
z [ K ( z ) Θ ( z , t ) z ] = ρ c Θ ( z , t ) t ,   z > 0 ,   t > 0 ,
K 0 Θ ( z , t ) z | z = 0 + = q 0 ,   t > 0 ,
Θ ( z , t ) 0 ,   z ,   t > 0 ,
Θ ( z , 0 ) = 0 ,   z 0 ,
and ρ , c are the density and the specific heat capacity of the material, respectively.
The following dimensionless variables and parameters were introduced:
ζ = z a ,   γ = γ a ,   τ = k t a 2 ,   k = K 0 c ρ ,   q = q q 0 ,   Θ = Θ Θ 0 ,   Θ 0 = q 0 a K 0 ,
where a is the effective heating depth, i.e., the distance from the heated surface, at which the temperature is equal to 5% of the maximum temperature on the heated surface [2]. Taking into account the designations (7), the problem (3–6) was written in following the form:
2 Θ * ( ζ , τ ) ζ 2 + γ * Θ * ( ζ , τ ) ζ e γ ζ Θ * ( ζ , τ ) τ = 0 ,   ζ > 0 ,   τ > 0 ,
Θ * ( ζ , τ ) ζ | ζ = 0 = 1 ,   τ > 0 ,
Θ * ( ζ , τ ) 0 ,   ζ ,   τ > 0 ,
Θ * ( ζ , 0 ) = 0 ,   ζ 0 .
Applying the following Laplace integral transform [34]:
Θ ¯ ( ζ , p ) L [ Θ ( ζ , τ ) ; p ] = 0 Θ ( ζ , τ ) e p τ d τ ,   Re p 0 ,
to the boundary-value problem (8)–(11), the following were obtained:
d 2 Θ ¯ * ( ζ , p ) d ζ 2 + γ d Θ ¯ * ( ζ , p ) d ζ p e γ ζ Θ ¯ * ( ζ , p ) = 0 ,   ζ > 0 ,
d Θ ¯ * ( ζ , p ) d ζ | ζ = 0 + = 1 p ,
Θ ¯ ( ζ , p ) 0 ,   ζ .
The solution to the boundary problem (13)–(15) takes the following form:
Θ ¯ * ( ζ , p ) = e 1 2 γ 1 ζ φ ( ζ , p ) Φ ( p ) ,   ζ 0 ,
where
φ ( ζ , p ) = I 1 ( 2 γ p e 1 2 γ ζ ) ,
Φ ( p ) = p p I 0 ( 2 γ p ) ,
I k ( x ) , k = 0 , 1 are the modified Bessel functions of the first kind [35]. Differentiating the function (18) with consideration of the relation of I 0 ( x ) = I 1 ( x ) , the following was found:
Φ ( p ) = 3 2 p I 0 ( 2 γ p ) + p γ I 1 ( 2 γ p ) .
The transition from the space of the transforms to the originals in the solution (16)–(19) was carried out on the basis of the Vashchenko-Zakharchenko theorem [36,37] as follows:
Θ ( ζ , τ ) = e 1 2 γ ζ [ lim p 0 φ ( ζ , p ) p Φ ( p ) + n = 1 φ ( ζ , p n ) Φ ( p n ) e p n τ ] ,   ζ 0 ,   τ 0 ,
where
I 0 ( 2 γ 1 p n ) 0 ,   n = 1 , 2 ,
By using the expressions of the modified Bessel functions [35] as follows:
I 0 ( x ) 1 + x 2 4 + x 4 64 + ,   I 1 ( x ) x 2 + x 3 16 +
and by limiting them to the first components, the following representations of functions (17) and (18) for the small values of the parameter p were obtained:
φ ( ζ , p ) 1 γ e 1 2 γ ζ p ,   Φ ( p ) p p .
Taking into account the Formulas (23) and relationships [35] as follows:
I 0 ( x ) = J 0 ( i x ) ,   I 1 ( x ) = i J 1 ( i x ) ,   i 1 ,
where J k ( x ) , k = 0 , 1 are the Bessel functions of the first kind, solutions (20) and (21) were written in the following form:
Θ ( ζ , τ ) = 1 γ e 1 2 γ ζ [ e 1 2 γ ζ 4 n = 1 J 1 ( μ n e 1 2 γ ζ ) μ n 2 J 1 ( μ n ) e λ n τ ] ,   ζ 0 ,   τ 0 ,
where
λ n = ( 1 2 γ μ n ) 2 ,
J 0 ( μ n ) 0 ,   n = 1 , 2 ,
By adopting ζ = 0 in the solution (25)–(27), the dimensionless temperature rise on the heated surface was found in the following form:
Θ ( τ ) Θ * ( 0 , τ ) = 1 γ * ( 1 4 n = 1 e λ n τ μ n 2 ) ,   τ 0 .
Verification of the developed model was carried out by checking the boundaries (9), (10), and the initial (11) conditions. For this, by differentiating the solution (25) with respect to the spatial variable ζ , and taking into account the relationship [35] as follows:
J 1 ( x ) = J 0 ( x ) x 1 J 1 ( x )
the following was found:
Θ ( ζ , τ ) ζ = e γ ζ [ 1 2 n = 1 J 0 ( μ n e 1 2 γ ζ ) μ n J 1 ( μ n ) e λ n τ ] ,   ζ 0 ,   τ 0 .
Approaching Formula (30) to the limit of ζ 0 , with consideration of the fact that μ n are the roots of Equation (27), the following was achieved:
Θ ( ζ , τ ) ζ | ζ = 0 + = 1 ,
which confirms that the boundary condition (9) was satisfied. However, by approaching to the limit ζ in the solution (25), we have confirmed the fulfillment of the condition of the temperature fade (10). In addition, it was numerically established that the initial condition (11) was met. In particular, according to the Formula (25), on the heated surface ζ = 0 at the initial point in time, the dimensionless temperature rise is zero if the following occurs:
n = 1 1 μ n 2 = 1 4 .
Based on the calculations, it was found that the sum of the first 102 components of the series on the left-hand side of Equation (32) is equal to 0.248985.
In addition to the exact (25), the appropriate asymptotic solutions to the problem at small and large values of the Fourier number τ (dimensionless time) were also found.
Small values τ (large values of the parameter p ). Taking into account that, in Formulas (17) and (18), the first components of the asymptotic of modified Bessel functions at large values of the argument [35] were as follows:
I 0 ( x ) e x 2 π x ( 1 + 1 8 x + ) ,   I 1 ( x ) e x 2 π x ( 1 + 3 8 x + ) ,
the transformed solution (16) was written as follows:
Θ ¯ * ( ζ , p ) = e 1 4 γ ζ e b p p p ,   b = 2 γ * ( 1 e 1 2 γ ζ ) ,   ζ 0 .
Using the relation [38] as follows:
L 1 [ e b p p p ; τ ] = 2 τ ierfc ( b 2 τ ) ,
from the Formula (34), the following form of dimensionless temperature rise at the initial moments of heating was obtained:
Θ ( ζ , τ ) 2 e 1 4 γ ζ τ ierfc ( b 2 τ ) ,   ζ 0 ,   0 τ < < 1 ,
where ierfc ( x ) = π 1 2 e x 2 x erfc ( x ) , erfc ( x ) = 1 erf ( x ) , erf ( x ) is the Gauss error function [35].
At ζ = 0 from the solution (36), a known result for the evolution of the temperature on the heated surface of a homogeneous semi-space was obtained [18] as follows:
Θ ( τ ) 2 τ π ,   0 τ < < 1 ,
Large values τ (small values of the parameter p ). By including the first two components in the distributions (22), from the transformed solution (16)–(18) it follows that:
Θ ¯ * ( ζ , p ) e γ * ζ γ * [ 1 p ( 1 1 2 e γ * ζ ) 1 ( γ * 2 + p ) ] ,   ζ 0 .
Considering the relations [38] as follows:
L 1 [ p 1 ; τ ] = 1 ,   L 1 [ ( γ * 2 + p ) 1 ; τ ] = e γ * 2 τ ,
from the solution (38), the following asymptotic representation for the dimensionless temperature rise at large values of time was obtained:
Θ * ( ζ , τ ) e γ * ζ γ * [ 1 ( 1 1 2 e γ * ζ ) e γ * τ ] ,   ζ 0 ,   τ > > 1 .
The temperature rise of the heated surface was found, substituting to the formula (40) ζ = 0 as follows:
Θ ( τ ) 1 γ * ( 1 1 2 e γ * τ ) ,   τ > > 1 .

3. Heating of the FGM Semi-Space by Heat Flux with the Intensity Linearly Decreasing in Time

As presented above, the exact (25), asymptotic (36), and (40) solutions were obtained at a constant intensity of heat flux of q = q 0 . This section concerns the heating of the surface of the semi-space by heat flux with the following time profile of intensity:
q ( t ) = q 0 q ( t ) ,   q ( t ) = 1 t t s ,   0 t t s ,
where t s is the stop moment of heating. The evolution of the heat flux intensity (42) corresponds to the temporal profile of the specific friction power during braking with constant deceleration [39]. The corresponding dimensionless temperature rise of Θ ^ * was searched based on Duhamel’s theorem [40] as follows:
Θ ^ * ( ζ , τ ) = τ 0 τ q ( τ s ) Θ * ( ζ , s ) d s ,   ζ 0 , 0 τ τ s ,
where Θ * is the dimensionless temperature rise (25)–(27), q ( τ ) is the function (42), and τ s is the dimensionless stop time as follows:
τ s = k t s a 2 .
After the integration with the next differentiation on the right-hand side of Formula (42), the following was obtained:
Θ ^ * ( ζ , τ ) = 1 γ * e 1 2 γ ζ [ q * ( τ ) e 1 2 γ ζ 4 n = 1 J 1 ( μ n e 1 2 γ ζ ) μ n 2 J 1 ( μ n ) G n ( τ ) ] , ζ 0 , 0 τ τ s ,
where
G n ( τ ) = e λ n τ ( 1 e λ n τ ) λ n τ s ,
the coefficients λ n were determined from Formula (26), and the numbers μ n , n = 1 , 2 , , are the single positive roots of Equation (27).
By putting ζ = 0 in the solution (45), the following formula was used to determine the evolution of the dimensionless temperature rise of the heated surface that was obtained:
Θ ^ * ( τ ) = 1 γ * [ q * ( τ ) 4 n = 1 G n ( τ ) μ n 2 ] ,   0 τ τ s .

4. The Heat Partition Ratio

We will generalize the results that were obtained above to the case of two ( l = 1 , 2 ) FGM semi-spaces that exponentially increased with the distance from the surface z = 0 of thermal conductivity as follows:
K l ( z ) = K l , 0 e γ l z ,   z 0 ,   K l , 0 K l ( 0 ) ,   γ l 0 ,   l = 1 , 2 .
The surfaces of z = 0 of each semi-space were heated with the heat flux at the following intensities:
q l ( t ) = α l q ( t ) ,   0 t t s ,   l = 1 , 2 ,   α 1 = α ,   α 2 = 1 α ,
where q ( t ) is the function (42) and α is the unknown heat partition ratio (HPR) (Figure 2).
Taking into account the form of the solutions (45) and (46), the temperature field in each semi-space heated by heat fluxes of intensity (49) respectively, was written in the following form:
T l ( z , t ) = T 0 + Θ l , 0 Θ ^ l ( ζ l , τ ) ,   z 0 ,   0 t t s ,   l = 1 , 2 ,
where
Θ ^ l ( ζ l , τ ) = 1 γ l e 1 2 γ l ζ l [ ( 1 τ ) e 1 2 γ l ζ l 4 n = 1 J 1 ( μ n e 1 2 γ l ζ l ) μ n 2 J 1 ( μ n ) G l , n ( τ ) ] ,   ζ l 0 ,   0 τ 1 ,
G l , n ( t ) = e λ l , n τ ( 1 e λ l , n τ ) λ l , n τ l , s ,
λ l , n = ( 1 2 γ l μ n ) 2 τ l , s ,
τ = t t s ,   ζ l = z a l ,   τ l , s = k l t s a l 2 ,   γ l = γ l a l ,   k l = K l , 0 c l ρ l ,   Θ ^ l = Θ ^ l Θ l , 0 ,   Θ l , 0 = α l q 0 a l K l , 0 , l = 1 , 2 ,
c l , ρ l are the specific heat capacity and the density of the materials, respectively, μ n > 0 , n = 1 , 2 , are the single roots of Equation (27).
Substituting z = 0 ( ζ l = 0 ) to the solution (51)–(54), the change in the temperature on the heated surface was obtained in the following form:
T l ( t ) = T 0 + Θ l , 0 Θ ^ l ( τ ) ,   0 t t s ,   l = 1 , 2 ,
where
Θ ^ l ( τ ) = 1 γ l [ 1 τ 4 n = 1 G l , n ( τ ) μ n 2 ] ,   0 τ 1 .
The coefficients of Θ l , 0 (54) in Formulas (50) and (55) contain an unknown heat partition ratio of α . In order to determine it, the equality condition of the mean temperature in time on the heated surfaces was used as follows:
T ˜ 1 = T ˜ 2 ,
where
T ˜ l = 1 t s 0 t s T l ( t ) d t ,   l = 1 , 2 ,
Considering the temperature on the heated surfaces (55) and (56) in the Formula (58), after integration was found as follows:
T ˜ l = T 0 + Θ l , 0 Θ ˜ l ,   l = 1 , 2 ,
where
Θ ˜ l = 1 2 γ l ( 1 8 n = 1 G ˜ l , n μ n 2 ) ,  
G ˜ l , n = ( 1 e λ l , n ) λ l , n ( 1 + 1 λ l , n τ l , s ) 1 λ l , n τ l , s ,   l = 1 , 2 ,
and the coefficients of λ l , n were determined from the formula (53).
Substituting the mean temperature (59)–(61) to Equation (57), the HPR was obtained in the following form:
α = K a Θ ˜ + K ,
where
a = a 1 a 2 ,   K = K 1 , 0 K 2 , 0 ,   Θ ˜ = Θ ˜ 1 Θ ˜ 2 .
If the effective heating depths are determined from the Formula [9] as follows:
a l = 3 k l t s ,   l = 1 , 2 ,
then from Formula (54) it follows that the dimensionless heating time of each semi-space is equal to τ l , s = 3 1 0.333 , and from Formulas (62) and (63) it follows that:
α = K ε Θ ˜ + K ε ,
where
K ε = K k ,   k = k 1 k 2 .
Omitting the component containing a series, on the right-hand side of Formula (60), i.e., assuming the following:
Θ ˜ l 1 2 γ l ,   l = 1 , 2 ,
where
α K γ a + K γ ,
and Formula (65) is as follows:
α K ε γ 1 + K ε γ ,
where
γ = γ 1 γ 2 .
With the same effective heating depths ( a = 1 ) and dimensionless gradients of materials ( γ = 1 ) from Formula (68), the known Block’ result [20] as follows:
α K 1 + K ,
and from Formula (69), the classic Charron’s Formula [41] was achieved as follows:
α K ε 1 + K ε ,
which is often used in analytical and numerical modeling of the frictional heating of homogeneous materials during braking [42,43,44]. For the same materials of both semi-spaces ( K = k = γ = a = 1 ), all of the above obtained formulas give the value of the heat partition ratio equal to α = 0.5 .
It should be noted that the received formulas for the determination of HPR need to be verified in the future with appropriate experimental data. Due to the lack of such opportunity at that moment, the authors would be grateful for the provision of such data or for carrying out the cooperative research. This would also allow us to establish the limits of the applicability of the formulas. However, already at this stage of research, we can assume the possibility of the practical use of the obtained formulas to determine HPR, due to the usage of the classic methodology that has been approved for homogeneous materials. Furthermore, the experimentally confirmed formulas for homogeneous materials can be obtained from the proposed solution for FGMs, as a results of the proper limit approach.
An important element of the most general Formula (62) for determining the HPR are the dimensionless time-averaged temperature rises Formulas (60) and (61) found in the case of the heat flux intensity linearly decreasing in time. This case is often considered when calculating the temperature of the brake systems that are operating in the mode with a sudden increase in the contact pressure to the nominal value at the beginning of braking. A classification of the remaining heat flux intensity temporal profiles has been proposed for homogeneous materials, without and with consideration of the contact pressure rise [45,46]. Obtaining the appropriate solutions in the case of FGM is also one of the directions of our research in the future.

5. Example of Calculation of the Heat Partition Ratio for an FGM Couple

Calculations were performed for two semi-limited bodies made of two-component FGMs. The first ( l = 1 ) element forms aluminum oxide Al2O3 (base, m = 0 ) and copper Cu (core, m = 1 ), and the other ( l = 2 ) contains zircon dioxide ZrO2 (base, m = 0 ) and titanium alloy Ti-6Al-4V (core, m = 1 ). The thermo-physical properties of these materials are demonstrated in Table 1.
The other operating parameters are the nominal heat flux intensity of q 0 = 3.78 MW m 2 , the braking time of t s = 121 s, and the temperature at the initial time moment of T 0 = 20   ° C [32]. The dimensionless material gradient parameters of each element were calculated from the relation [47] as follows:
γ l = ln ( K l , 1 K l , 0 ) ,   l = 1 , 2 .
Obtaining the values of γ 1 * = 2.38 , γ 2 * = 1.26 . The specific heat capacity and density of materials of the heated elements were determined according to the following mixture law:
c l = c l , 1 v + ( 1 v ) c l , 0 ,   ρ l = ρ l , 1 v + ( 1 v ) ρ l , 0 ,   l = 1 , 2 ,
where v is the volume fractions of the base and core components. For their equal participation ( v = 0.5 ), the following properties of FGM were found:
c 1 = 437.32 J kg 1 K 1 ,   ρ 1 = 6469.42 kg m 3 ,   k 1 = 1.32 · 10 5 m 2 s 1 ,
c 2 = 495.46 J kg 1 K 1 ,   ρ 2 = 5266.96 kg m 3 ,   k 2 = 7.43 · 10 7 m 2 s 1 .
From Equation (64) the effective depths of the heat penetration were determined as a 1 = 21.8 mm , a 2 = 5.2 mm . Next, based on Formulas (63), (66), and (70), the values of the dimensionless parameters were calculated as follows: a = 4.208 , K * = 19.196 , k * = 17.706 , K ε = 4.562 , γ * = 1.883 . This allowed for an estimation by means of Equations (60), (61) dimensionless, time-averaged temperature rises of Θ ˜ 1 = 0.329 , Θ ˜ 2 = 0.703 , and their ratio of Θ ˜ = 0.468 (63).
By substituting the found parameter values successively to the right-hand sides of Formulas (65) and (69), the proper values of the HPR were calculated as α = 0.907 and α = 0.896 . A slight (1.2%) difference in the obtained results allowed us to analyze the influence of the dimensionless parameters of the thermal activity K ε (66) and the relative gradient of the FGMs γ (70) on the HPR value (Figure 3), on the basis of Formula (69) only. For a fixed value of the thermal activity coefficient K ε , a rise of the parameter γ causes an increase in the amount of heat that is directed to the first element of the heated couple. Conversely, by increasing the thermal activity of the friction pair at a predetermined value γ , the amount of heat that is directed to the first element increases.
For the HPR values α 1 = 0.896 , α 2 = 0.104 (49), from Equation (54), the scaling factors of the temperature rises were determined to be Θ 1 , 0 = 2034 ° C and Θ 2 , 0 = 1080 ° C . Then, based on the relations Formulas (55) and (56) the evolutions of the temperature on the heated surfaces of each element were found (Figure 4).
It has been established that the temporal profiles of the temperature of the elements are different, in particular in the final stage of the heating process. The time course of the temperature of the first element ( l = 1 , Al2O3–Cu) is typical for the evolution of the friction surface temperature during braking with a constant deceleration—a rapid increase in the temperature at the beginning of braking, reaching its maximum value in the middle of the process, followed by a temperature reduction until the standstill. However, in the second element ( l = 2 , ZrO2–Ti–6Al–4V) a rise of temperature on the heated surface is monotonic during the whole process. Such temperature behavior is decisively influenced by the thermo-physical properties of the component materials of each element. In the functionally graded friction couple under consideration, the materials of both of the components of the first element have a significantly greater ability to dissipate the heat from the heated surface than the materials of the second element (Table 1). Moreover, this is confirmed, by the values of the coefficients α 1 = 0.896 , α 2 = 0.104 , which prove that the first element absorbs almost 90 % , and the second only slightly more than 10 % of the entire heat flux intensity q (42). Due to the low thermal conductivity of the component materials of the second element, in particular the zircon dioxide, the temperature of the heated surface of this element continues to rise during heating, even with a linearly decreasing intensity of the heat flux.
The results that are demonstrated in Figure 3 and Figure 4 were obtained by using the HPR value α = 0.896 that was found from the dependency (69). The approximated formula for determining the α (69), as well as the exact Equation (65), was obtained by assuming that the effective depths of the heating of the elements a l , l = 1 , 2 were determined from the empirical Formula (64), so the corresponding values of the Fourier numbers τ l , s , l = 1 , 2 (54) were the same and equal to 1/3. The introduction of Formula (64) to the model is justified when the thickness of the heated element is greater than the appropriate value al, l = 1, 2 [2]. In general, for the determination of the temperature mode at the design stage of the brake, the Fourier numbers τ l , s ,   l = 1 , 2 are given, and the effective heating depths based on the relation (54) are calculated in the following form:
a l = k l t s τ l , s ,   l = 1 , 2 .
With consideration of Equation (75) for the selected friction couple and the established time of heating (braking) t s = 12.1 s the ratio of the depths of heating a (63) is determined from the following relation:
a = k τ s ,   τ s = τ 1 , s τ 2 , s ,
where k * is the ratio of the thermal diffusivities of the FGMs (66). By treating the Fourier numbers τ l , s , l = 1 , 2 (54) as independent variables in Formula (62), their influence on the HPR α was investigated with the previously determined values of the dimensionless parameters, K * = 19.196 , k * = 17.706 , γ 1 * = 2.38 , and γ 2 * = 1.26 for a given pair of elements. The dimensionless parameter Θ ˜ (63) included in Formula (62) was determined from Formulas (60) and (61), and the finding of a was made by using Formula (76). The results of the calculations are presented in Figure 5.
At a fixed value τ 2 , s , the HPR α quickly increases with increasing τ 1 , s , reaching its nominal value at τ 1 , s 1.5 (Figure 5a). The largest nominal value α is achieved for τ 1 , s = 1 / 3 , i.e., when determining the effective heating depth using Formula (64). The results obtained at τ 2 , s = 1 and τ 2 , s = 2 practically coincide. A different nature of the HPR α change occurs with a fixed value of the Fourier number τ 1 , s and increasing the values of τ 2 , s (Figure 5b). By increasing τ 2 , s , the amount of heat that is absorbed by the first element decreases, reaching a minimum value at τ 2 , s 1.5 . It should be noted that the Fourier numbers τ l , s , l = 1 , 2 (54) are dimensionless input parameters that play an important role in estimating the temperature mode of the brakes. The designers calculate their values after determining the geometric parameters a l , l = 1 , 2 (63). In the quality of these parameters, they choose the thickness of the friction elements, or, if they are greater than the effective depth of the heat penetration (64), then the latter. Then, having previously found the maps of those, as shown in Figure 5, the designer can easily estimate the HPR of a given braking system.

6. Conclusions

A methodology was proposed in order to determine the heat partition ratio (HPR) in the braking systems with friction pair elements that were made of functionally graded materials (FGMs). The basis of this methodology is an exact solution to the initial-boundary problem of heat conduction that is formulated for the FGM semi-space that is heated on its surface by a heat flux with an intensity that decreases linearly with time and includes an unknown a priori HPR. For its determination, the condition of the equality of the time-averaged temperature on the heated surfaces of the two different semi-spaces was used.
As an example, the heat partition ratio was found for the FGM friction couple that consisted of Al2O3–Cu and ZrO2–Ti–6Al–4V. It was established that element Al2O3–Cu absorbs most of the heat that is generated due to friction (almost 90%). The maximum temperature on the friction surface of the Al2O3–Cu element is about 820 °C and is achieved in the middle of the braking time. However, the highest temperature of the friction surface of the ZrO2–Ti–6Al–4V element is achieved at the stop moment and amounts to 970 °C. Thus, the crucial influence on the evolution of the temperature in the functionally graded friction couple have the thermo-physical properties of the component materials.
Then, the simulation of the HPR dependency on the dimensionless parameters, such as the thermal activity of the friction couple, the gradient materials, and the Fourier numbers was carried out. It was found that increasing the values of the friction couple thermal activity or the gradient materials ratio causes an increase in the HPR. The scheme of the map development for a given functionally graded friction pair was proposed, which allows for a quick estimation of the HPR, depending on such input parameters as their Fourier numbers. The importance of having such maps for the designer is indicated by the results that have been obtained in the case of the homogeneous materials [9,28].

Author Contributions

Conceptualization and methodology, A.Y.; software, K.T. and P.Z.; validation, P.Z.; formal analysis, K.T. and P.Z.; investigation, A.Y., K.T. and P.Z.; writing—original draft preparation, K.T. and P.Z.; writing–review and editing, A.Y., K.T. and P.Z.; visualization and figures preparation, K.T.; supervision, A.Y.; project administration, K.T.; funding acquisition, 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 obtained 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

a Effective depth of heat penetration ( m   )
c Specific heat capacity ( J kg 1 K 1 )
f 0 Coefficient of friction
I k ( ) Modified Bessel functions of the first kind of the kth order
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
p Dimensionless parameter of the Laplace integral transform
q Intensity of heat flux ( W m 2 )
q 0 Nominal intensity of the heat flux ( W m 2 )
t Time ( s )
t s Braking time ( s )
T Temperature ( C )
T 0 Initial temperature ( C )
v Volume fraction of the material phases (dimensionless)
z Spatial coordinate in axial direction ( m   )
α Heat partition ratio
γ Parameter of material gradient ( m 1 )
γ FGMs gradient ratio
Θ Temperature rise ( C )
Θ Dimensionless temperature rise
Θ 0 Scaling factor of temperature rise ( C )
ρ Density ( kg m 3 )
τ * Dimensionless time
τ s Dimensionless braking time
ζ Dimensionless spatial coordinate in axial direction

References

  1. Padture, N.P.; Gell, M.; Jordan, E.H. Thermal Barrier Coatings for Gas-Turbine Engine Applications. Science 2002, 296, 280–284. [Google Scholar] [CrossRef] [PubMed]
  2. Mahamood, R.M.; Akinlabi, E.T. Functionally Graded Materials; Springer: Cham, Switzerland, 2017. [Google Scholar]
  3. Zhou, W.; Ai, S.; Chen, M.; Zhang, R.; He, R.; Pei, Y.; Fang, D. Preparation and thermodynamic analysis of the porous ZrO2/(ZrO2 + Ni) functionally graded bolted joint. Compos. Part B Eng. 2015, 82, 13–22. [Google Scholar] [CrossRef]
  4. Zhou, W.; Zhang, R.; Ai, S.; He, R.; Pei, Y.; Fang, D. Load distribution in trends of porous metal-ceramic functionally graded composite joints subjected to thermomechanical loading. Compos. Struct. 2015, 134, 680–688. [Google Scholar] [CrossRef] [Green Version]
  5. Mahmoud, D.; Elbestawi, M. Lattice Structures and Functionally Graded Materials. Applications in Additive Manufacturing of Orthopedic Implants: A Review. J. Manuf. Mater. Process. 2017, 1, 13. [Google Scholar] [CrossRef]
  6. Hovorun, T.P.; Berladir, K.V.; Pererva, V.I.; Rudenko, S.G.; Martynov, A.I. Modern Materials for Automotive Industry. J. Eng. Sci. 2017, 4, 8–18. [Google Scholar] [CrossRef] [Green Version]
  7. Strojny-Nędza, A.; Pietrzak, K.; Gili, F.; Chmielewski, M. FGM based on copper-alumina composites for brake disc applications. Arch. Civ. Mech. Eng. 2020, 20, 83. [Google Scholar] [CrossRef]
  8. Newcomb, T.P. Temperatures reached in disc brakes. J. Mech. Eng. Sci. 1960, 2, 167–177. [Google Scholar] [CrossRef]
  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. Day, A.J. Braking of Road Vehicles; Butterworth-Heinemann: Oxford, UK, 2014. [Google Scholar]
  11. Ling, F.F. Surface Mechanics; John Wiley & Sons: New York, NY, USA, 1973. [Google Scholar]
  12. Barber, J.R.; Comninou, M. Thermoelastic contact problems. In Thermal Stresses, 3rd ed.; Hetnarsky, R.R., Ed.; Elsevier: Amsterdam, The Netherlands, 1989. [Google Scholar]
  13. Awrejcewicz, J.; Pyryev, Y. Nonsmooth Dynamics of Contacting Thermoelastic Bodies; Springer: New York, NY, USA, 2009. [Google Scholar]
  14. Yevtushenko, A.A.; Kuciej, M. One-dimensional thermal problem of friction during braking: The history of development and actual state. Int. J. Heat Mass Transf. 2012, 55, 4148–4153. [Google Scholar] [CrossRef]
  15. Evtushenko, O.; Kuciej, M.; Topczewska, K. Determination of the maximal temperature of a pad-disk tribosystem during one-time braking. Mater. Sci. 2020, 56, 152–159. [Google Scholar] [CrossRef]
  16. Pyryev, Y.; Yevtushenko, A. The influence of the brakes friction elements thickness on the contact temperature and wear. Heat Mass Transf. 2000, 36, 319–323. [Google Scholar] [CrossRef]
  17. Yevtushenko, A.A.; Grześ, P. The FEM-modeling of the frictional heating phenomen in the pad/disc tribosystem (a review). Num. Heat Trans. Part A-Appl. 2010, 58, 207–226. [Google Scholar] [CrossRef]
  18. Wasilewski, P. Frictional heating in railway brakes: A review of numerical models. Arch. Computat. Methods Eng. 2020, 27, 45–58. [Google Scholar] [CrossRef]
  19. Deressa, K.T.; Ambie, D.A. Thermal load simulations in railway disc brake: A systematic review of modelling temperature, stress and fatigue. Arch. Computat. Methods Eng. 2022, 29, 2271–2283. [Google Scholar] [CrossRef]
  20. Blok, H. Theoretical field study of temperature rise at surfaces of actual contact under oiliness lubricating conditions. Proc. Instn. Mech. Engrs. (Gen. Discuss. Lubr. Lubr.) 1937, 45, 222–235. [Google Scholar]
  21. Day, A.J.; Newcomb, T.P. Dissipation of frictional energy from the interface of an annular disc brake. Proc. Inst. Mech. Eng. Part D J. Automob. Eng. 1984, 198, 201–209. [Google Scholar] [CrossRef]
  22. Grześ, P. Partition of heat in 2D finite element model of a disc brake. Acta Mech. Autom. 2011, 5, 35–41. [Google Scholar]
  23. Yevtushenko, A.A.; Grześ, P. Finite element analysis of heat partition ratio in a pad/disc brake system. Num. Heat Trans. Part A-Appl. 2011, 58, 521–542. [Google Scholar] [CrossRef]
  24. Qiu, L.Q.; Qi, H.S.; Wood, A. Two-dimensional finite element analysis of the heat partition ratio of a friction brake. Proc. Inst. Mech. Eng. Part J J. Eng. Trib. 2018, 232, 1489–1501. [Google Scholar] [CrossRef] [Green Version]
  25. Afzal, A.; Abdul Mujeebu, M. Thermo-Mechanical and Structural Performances of Automobile Disc Brakes: A Review of Numerical and Experimental Studies. Arch. Comput. Methods Eng. 2019, 26, 1489–1513. [Google Scholar] [CrossRef]
  26. Carslaw, H.C.; Jaeger, J.C. Conduction of Heat in Solids, 2nd ed.; Clarendon Press: Oxford, UK, 1959. [Google Scholar]
  27. Hasselgruber, H. Der Schaltvorgang einer Trockenreibung Kupplung bei kleinster Erwärmung. Konstruction 1963, 15, 41–45. [Google Scholar]
  28. Ginzburg, A.G. Coefficient of distribution of the heat flows during braking. In Calculation and Testing of Friction Couples; Mechanical Engineering: Moscow, Russia, 1974. (In Russian) [Google Scholar]
  29. Yevtushenko, A.; Kuciej, M.; Topczewska, K. Effect of the Temporal Profile of the Friction Power on Temperature of a Pad-Disc Brake System. J. Theoret. Appl. Mech. 2019, 57, 461–473. [Google Scholar] [CrossRef]
  30. Xiong, C.; Chen, M.; Yu, L. Analytical model and material equivalent methods for steady state heat partition coefficient between two contact discs in multi-disc clutch. Proc. J. Mech. Eng. J. Automob. Eng. 2019, 234, 857–871. [Google Scholar] [CrossRef]
  31. Xia, Y.; Yano, A.; Hayashi, N.; Noraguchi, N.; Xie, G. Analysis of temperature and heat partitioning coefficient during friction between polymer and steel. Tribol. Int. 2022, 171, 107561. [Google Scholar] [CrossRef]
  32. 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]
  33. 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]
  34. Sneddon, I.N. The Use of Integral Transforms; McGraw-Hill: New York, NY, USA, 1972. [Google Scholar]
  35. Abramowitz, M.; Stegun, I. Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables; United States Department of Commerce, National Bureau of Standards (NBS): Washington, DC, USA, 1964. [Google Scholar]
  36. Vashchenko-Zakharchenko, M.E. Symbolic Calculus and Its Application to Integration of Linear Differential Equations; University Press: Kiev, Ukraine, 1862. [Google Scholar]
  37. Yevtushenko, A.; Topczewska, K.; Zamojski, P. The Effect of Functionally Graded Materials on Temperature during Frictional Heating: Under Uniform Sliding. Materials 2021, 14, 4285. [Google Scholar] [CrossRef]
  38. Bateman, H.; Erdélyi, A. Tables of Integral Transforms; McGraw-Hill: New York, NY, USA, 1954; Volume 1. [Google Scholar]
  39. Yevtushenko, A.; Topczewska, K.; Zamojski, P. The Effect of Functionally Graded Materials on Temperature during Frictional Heating at Single Braking. Materials 2021, 14, 6241. [Google Scholar] [CrossRef]
  40. Özişik, N.M. Heat Conduction; John Wiley & Sons: New York, NY, USA, 1993. [Google Scholar]
  41. Charron, F. Partage de la Chaleur Entre Deux Corps Frottants; Publications Scientifiques et Techniques; Blondel La Rougery: Paris, France, 1943; 182p. [Google Scholar]
  42. Talati, F.; Jalalifar, S. Analysis of heat conduction in a disk brake system. Heat Mass Transf. 2009, 45, 1047–1059. [Google Scholar] [CrossRef]
  43. 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]
  44. Loizou, A.; Qi, H.S.; Day, A.J. A fundamental study on the heat partition ratio of vehicle disk brakes. J. Heat Transf. 2013, 135, 121302. [Google Scholar] [CrossRef]
  45. Yevtushenko, A.; Kuciej, M.; Topczewska, K. Analytical model for investigation of the effect of friction power on temperature in the disk brake. Adv. Mech. Eng. 2017, 9, 1–12. [Google Scholar] [CrossRef]
  46. Topczewska, K. Influence of the Time of Increase in Contact Pressure in the Course of Braking on the Temperature of a Pad–Disc Tribosystem. Mater. Sci. 2018, 54, 250–259. [Google Scholar] [CrossRef]
  47. Mao, J.J.; Ke, L.L.; Yang, J.; Kitipornchai, S.; Wang, Y.S. The coupled thermoelastic instability of FGM coatings with arbitrarily varying properties: In-plane sliding. Acta Mech. 2018, 229, 2979–2995. [Google Scholar] [CrossRef]
Figure 1. Scheme of the heating of the FGM semi-space.
Figure 1. Scheme of the heating of the FGM semi-space.
Materials 15 04623 g001
Figure 2. Scheme of separation of the friction pair elements.
Figure 2. Scheme of separation of the friction pair elements.
Materials 15 04623 g002
Figure 3. Isolines of the heat partition ratio α (69) in the coordinate system ( γ , K ε ).
Figure 3. Isolines of the heat partition ratio α (69) in the coordinate system ( γ , K ε ).
Materials 15 04623 g003
Figure 4. Evolutions of temperature T l , l = 1 , 2 (55) and (56) on the heated surfaces z = 0 of the considered semi-spaces made of FGMs.
Figure 4. Evolutions of temperature T l , l = 1 , 2 (55) and (56) on the heated surfaces z = 0 of the considered semi-spaces made of FGMs.
Materials 15 04623 g004
Figure 5. Dependences of the heat partition ratio α (62) on: (a) Fourier number τ 1 , s for different values of τ 2 , s ; (b) Fourier number τ 2 , s for different values of τ 1 , s .
Figure 5. Dependences of the heat partition ratio α (62) on: (a) Fourier number τ 1 , s for different values of τ 2 , s ; (b) Fourier number τ 2 , s for different values of τ 1 , s .
Materials 15 04623 g005
Table 1. Thermo-physical properties of FGM components [32].
Table 1. Thermo-physical properties of FGM components [32].
Element SubscriptMaterial Subscript, mMaterialThermal Conductivity
K l , m ( 0 ) , Wm 1 K 1
Specific Heat Capacity
c l , m ( 0 ) , J kg 1 K 1
Density   ρ l , m ( 0 ) , kgm 3
l = 1 base, 0Al2O337.24727.293990.92
core, 1Cu402.65147.358947.92
l = 2 base, 0ZrO21.94452.836102.16
core, 1Ti-6Al-4V6.87538.084431.79
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yevtushenko, A.; Topczewska, K.; Zamojski, P. The Heat Partition Ratio during Braking in a Functionally Graded Friction Couple. Materials 2022, 15, 4623. https://0-doi-org.brum.beds.ac.uk/10.3390/ma15134623

AMA Style

Yevtushenko A, Topczewska K, Zamojski P. The Heat Partition Ratio during Braking in a Functionally Graded Friction Couple. Materials. 2022; 15(13):4623. https://0-doi-org.brum.beds.ac.uk/10.3390/ma15134623

Chicago/Turabian Style

Yevtushenko, Aleksander, Katarzyna Topczewska, and Przemysław Zamojski. 2022. "The Heat Partition Ratio during Braking in a Functionally Graded Friction Couple" Materials 15, no. 13: 4623. https://0-doi-org.brum.beds.ac.uk/10.3390/ma15134623

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