Next Article in Journal
Fatigue Crack Initiation Change of Cast AZ91 Magnesium Alloy from Low to Very High Cycle Fatigue Region
Next Article in Special Issue
Comparative Analysis of Temperature Fields in Railway Solid and Ventilated Brake Discs
Previous Article in Journal
High-Temperature Compressive Response of SiCp/6092Al Composites under a Wide Range of Strain Rates
Previous Article in Special Issue
Effect of Frictional Slipping on the Strength of Ribbon-Reinforced Composite
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Effect of Functionally Graded Materials on Temperature during Frictional Heating at Single Braking

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.
Materials 2021, 14(21), 6241; https://doi.org/10.3390/ma14216241
Submission received: 24 September 2021 / Revised: 14 October 2021 / Accepted: 19 October 2021 / Published: 20 October 2021

Abstract

:
A mathematical model for evaluation of the temperature mode of the disc–pad system during single braking is proposed. The model is based on the thermal problem of friction formulated for two semi-infinite bodies, compressed with pressure increasing over time while reducing the sliding velocity from the initial value to zero at the stop. The exact solution to this problem was obtained by means of Duhamel’s theorem. Validation of the solution was performed by achieving in special cases parameters of known solution to this problem with constant pressure and velocity (under uniform sliding). The results of the numerical calculations are presented for a selected friction pair, made of functionally graded materials with titanium alloy (disc) and aluminum alloy (pad) cores coated with ceramics graded toward friction surfaces. For the established values of the parameters such as the rise time in pressure and the FGM gradients, the ability to quickly obtain spatiotemporal temperature distributions in the disc and pad was presented. The influence of the variability of these parameters on the maximum temperature of the brake system was also investigated.

1. Introduction

Friction elements of braking systems are subjected to severe conditions such as high temperature and intensive wear. During braking, performance of these components in terms of efficiency, service life, and dissipation of heat from the contact surface depends on the operating conditions and material properties. It has been shown that the operating characteristics of the entire assembly of the braking system can be significantly improved by introducing a smooth gradient in the microstructure of the friction materials [1]. Such functionally graded materials (FGMs) are a class of heterogeneous materials with continuous variation of properties over their volume. Generally, these materials are composites formed by smooth gradation of two or more constituent phases along certain dimensions of a structure. This gradation can be regulated by changing the volume fraction distribution of component of material from one to another in a controlled manner [2]. As a result, the thermophysical properties of material continuously vary as a function of position along a certain direction. This allows designing a functionally graded material in order to obtain optimized friction characteristics of a brake.
In general, statements of thermal problems of friction contain partial differential equations with variable coefficients. Therefore, the application of analytical methods to their solution is difficult or even impossible. Hence, numerical methods are often used to consider such problems. An FGM disc subjected to thermal load due to frictional heating while taking into account the inertial force due to the rotation of the disc was studied by Afsar and Go [3]. A 2D finite element analysis (FEA) for a circular disc with exponential variations in thermophysical properties in the radial direction was performed. An axisymmetric FEA of a brake disc, with properties distributed according to the power-law function of radial position, was executed by Shahzamanian et al. [4,5]. It was found that the gradation index of the disc material has a crucial influence on the thermomechanical behavior of the entire braking system.
A finite element thermal contact analysis of a functionally graded disc under dry friction was performed by Hosseini and Talebi [6]. The core part of the considered disc was steel, and it gradually changed through the thickness of the disc, according to a power law, approaching pure ceramic at the outer surface. It was shown that the temperature and the corresponding thermal displacements in the FGM disc are much lower than in the conventional steel disc. Furthermore, it was established that the use of an FGM brake disc may eliminate thermal cracking and wear. In particular, functionally graded materials composed of ceramic and metal perform very well in contact problems involving friction, since they combine the advantages of both components [7]. These elements mostly have a metal core in order to maintain strength and rigidity, whereas ceramic is present on the outer surfaces to resist intensive wear and elevated temperature conditions.
Separately, the study of the phenomenon of thermoelastic instability (TEI) of brake systems with FGM should be mentioned. The solution of a 2D thermal contact problem of friction for a functionally graded cermet brake disc was obtained by means of FEA [8]. They investigated TEI caused by the coupled interaction of the mechanical and thermal loads in the sliding system. Generally, this leads to the establishment of localized high-temperature zones on the contact surface, known as hot spots, which are directly attributable to the premature failure of the friction system. This instability is often called the frictionally excited TEI and occurs in tribosystems when the sliding velocity exceeds a certain critical value. It was shown that the value of the critical velocity for a functionally graded brake disc is higher than that for a conventional homogeneous disc. This conclusion was confirmed by further research investigating TEI in an FGM strip sliding uniformly against two homogeneous semi-spaces [9,10]. Assuming an exponential variation of the thermophysical properties along the thickness of an FG strip permitted obtaining an exact solution using the analytical perturbation method. Using the same methodology, the TEI of the brake modeled as an FGM semi-infinite body sliding against a homogeneous semi-space under uniform pressure taking into account the frictional heating and thermal contact resistance was investigated by Mao et al. [11]. As a result, they determined the stability boundaries of thermoelastic instability in the considered sliding system. The effect of the arbitrarily varying thermoelastic properties of the FGM on the TEI was considered by Mao et al. [12]. To simulate the distribution of the FGM properties, a homogeneous multilayered model was employed. This approach is a replacement of the continuous FGM material with a package of homogeneous layers with constant properties. The gradient was simulated by assigning different properties values to each sublayer. It was proven that the results received for the FGM strip, divided into a sufficient number of layers, were close to the results found using the corresponding exact solutions [13]. It should be noted that this conclusion is dependent on the problem under consideration, and the differences between the obtained results may be significant in some cases [14]. This is particularly true for FGM with temperature-dependent properties. A multilayered model was used in [15] in order to establish the coupled effect of the frictional heat and the thermal contact resistance. Since the homogeneous multilayered model deals with the arbitrarily varying properties of FGM, the power-law, exponential, sinusoidal, and cosinusoidal distributions of the brake disc properties were considered. The perturbation and transfer matrix methods were used to deduce the characteristic equation of the TEI problem, to obtain the relationship between the critical sliding velocity and the critical heat flux [15]. The formulated conclusions confirm that the application of ceramic-based FGM in a brake disc, consisting of ceramics at the sliding interface and steel in the middle layer, reduces the susceptibility of braking system toward TEI [9,10,11,15].
However, FGMs are increasingly finding applications in braking systems, in the fabrication of not only discs, but also brake pads. Experimental investigations revealed that FGMs could successfully fulfill the demands for brake pads and improve their characteristics [1,16]. The novel functionally graded ductile iron for brake pads was investigated in a well-controlled model sliding test and a disc-brake machine in [16]. The results of the tribological tests revealed a positive effect of the functional gradient of properties on the wear of the pads and the improved stability of the friction coefficient. Govindaraju et al. [1] developed and investigated Fe-based material on brake pads with graded composition. The FG specimens were subjected to a dry sliding test for studying their tribological behavior. The results were compared with the conventional brake pad specimen. It was found that the wear resistance of the functionally gradient specimen is much greater compared to the conventional pad material [1].
We note that a more comprehensive review of the literature on thermoelastic contact problems with frictional heating for functionally graded materials was provided in our previous article [17]. This article is a continuation of the research cycle started in [17], in which the case of the uniform sliding of an FGM tribosystem was considered. The present article concerns the transient thermal problem of friction during braking, which takes into account the time-dependent specific friction power due to the exponential increase in contact pressure.

2. Statement to the Problem

The frictional heating in a brake disc system during a single braking process is considered. Frictional elements of the system are two identical pads, located symmetrically to the brake disc. At the initial time moment t = 0 , pads are pressed to the friction surfaces of the disc with uniformly distributed on the contact area and time-dependent pressure [18].
p ( t ) = p 0 p * ( t ) , p * ( t ) = 1 e t / t i , 0 t t s ,
where t i 0 is the rise time in contact pressure from zero to the nominal value p 0 , and t s is the time of stop. Due to the interaction of friction forces, the linear velocity of vehicle V is reduced from the initial value V 0 V ( 0 ) to zero at the stop time moment t = t s according to the following law [19,20]:
V ( t ) = V 0 V * ( t ) ,   V * ( t ) = 1 t t s 0 + t i t s 0 p * ( t ) ,   t s 0 = W 0 f p 0 A a V 0 ,   0 t t s ,
where W 0 is the initial kinetic energy of the system, f is the friction coefficient, A a is the nominal contact area between the pad and disc, and t s 0 is the braking time with constant deceleration ( t i 0 ). The braking time, taking into account the temporal profile of the velocity (Equation (2)), is determined from the stop condition V * ( t s ) = 0 . For 0 < t i 0.3 t s 0 it was established [19] that t s t s 0 + 0.99 t i .
The sliding velocity reduction during braking is accompanied by the generation of frictional heat on the contact surface of the friction pair. In order to determine the temperature field generated in this way, the corresponding thermal problem of friction is formulated on the basis of the following assumptions:
  • The materials of the pads and the disc are functionally graded with an exponential decrease in thermal conductivity along their thickness, with invariant specific heat and density;
  • The initial temperature of all elements is the same and equal to the ambient temperature T a ;
  • The whole work of friction goes to heating the bodies, while the wear of the friction surfaces is neglected;
  • The free surfaces of the pads and the disc are adiabatic;
  • The thermal and mechanical properties and coefficient of friction are independent of the temperature T ;
  • Only the change in the temperature gradient in the direction perpendicular to the friction surface is taken into account;
  • The thermal contact of friction between the pads and the disc is perfect; the temperatures of their friction surfaces during braking are the same, and the sum of the intensity of the heat fluxes directed to both elements along the normal to the contact surface is equal to the specific friction power:
    q ( t ) = q 0 q * ( t ) ,   q 0 = f p 0 V 0 ,   q * ( t ) = p * ( t ) V * ( t ) ,   0 t t s ,
    where the temporal profiles of pressure p * ( t ) and velocity V * ( t ) have the forms expressed in Equations (1) and (2), respectively;
  • Due to the symmetry with respect to the center plane of the disc, to establish the temperature of the braking system, it is sufficient to consider the contact scheme of one pad with a disc of half of its thickness.
With such assumptions, a contact scheme of two sliding semi-infinite bodies (semi-spaces) related to the Cartesian system 0 x y z (Figure 1) was adopted to describe the process of frictional heating in the disc–pad system.
The temperature rise Θ = T T a was determined from solution to the following one-dimensional boundary value problem of heat conduction taking into account the generation of heat due to friction:
z [ K 1 ( z ) Θ ( z , t ) z ] = c 1 ρ 1 Θ ( z , t ) t , z > 0 ,   0 < t t s ,
z [ K 2 ( z ) Θ ( z , t ) z ] = c 2 ρ 2 Θ ( z , t ) t , z < 0 ,   0 < t t s ,
K 2 ( z ) Θ ( z , t ) z | z = 0 K 1 ( z ) Θ ( z , t ) z | z = 0 + = q ( t ) , 0 < t t s ,
Θ ( 0 , t ) = Θ ( 0 + , t ) , 0 < t t s ,
Θ ( z , t ) 0 ,   | z | , 0 < t t s ,
Θ ( z , 0 ) = 0 , | z | < ,
where
K l ( z ) = K l , 0 e γ l | z | , | z | < , K l , 0 K l ( 0 ) , γ l 0 ,   l = 1 ,   2 ,
and function q ( t ) has the form expressed in Equation (3). Here and further, the subscript l indicates the parameters and quantities related to the certain element— l = 1 for the disc, and l = 2 for the pad. Taking into account the relations in Equation (10), the problem in Equations (6)–(9) was written in the following form:
2 Θ ( z , t ) z 2 + γ 1 Θ ( z , t ) z = e γ 1 z k 1 , 0 Θ ( z , t ) t , z > 0 ,   0 < t t s ,
2 Θ ( z , t ) z 2 γ 2 Θ ( z , t ) z = e γ 2 z k 2 , 0 Θ ( z , t ) t , z < 0 ,   0 < t t s ,
K 2 , 0 Θ ( z , t ) z | z = 0 K 1 , 0 Θ ( z , t ) z | z = 0 + = q ( t ) ,   0 < t t s ,
Θ ( 0 , t ) = Θ ( 0 + , t ) ,   0 < t t s ,
Θ ( z , t ) 0 ,   | z | , 0 < t t s ,
Θ ( z , 0 ) = 0 , | z | < ,
where
k l , 0 = K l , 0 c l ρ l ,   l = 1 ,   2
are the coefficients of thermal diffusivity of the materials on their contact surfaces; z = 0 .

3. Solution to the Problem

In the case of a uniform slip with a constant specific power of friction q ( t ) = q 0 , t 0 , the solution to the problem in Equations (11)–(16) can be written in the following form [17]:
Θ ^ ( z , t ) = Λ e γ 1 z / 2 [ e γ 1 z / 2 ( 1 + γ ε K ε ) + 4 γ ε n = 1 ϕ 1 ( z , μ n ) Ψ ( μ n ) e p n t ] , z 0 , t 0 ,
Θ ^ ( z , t ) = Λ e γ 2 z / 2 [ e γ 2 z / 2 ( 1 + γ ε K ε ) + 4 γ ε n = 1 ϕ 2 , n ( z , ) Ψ ( μ n ) e p n t ] , z 0 , t 0 ,
where
ϕ 1 ( z , μ n ) = J 1 ( μ n ) J 1 ( γ ε μ n e γ 1 z / 2 ) ,   ϕ 2 ( z , μ n ) = J 1 ( μ n ) J 1 ( γ ε μ n e γ 2 z / 2 ) ,
Ψ ( μ n ) = μ n 2 [ ( 1 + γ ε K ε ) J 0 ( μ n ) J 0 ( γ ε μ n ) ( γ ε + K ε ) J 1 ( μ n ) J 1 ( γ ε μ n ) ] ,
K ε = K 0 * k 0 * ,   γ ε = γ * k 0 * ,   K 0 * = K 1 , 0 K 2 , 0 ,   k 0 * = k 1 , 0 k 2 , 0 ,   γ * = γ 1 γ 2 ,   Λ = q 0 γ 2 K 2 , 0 ,
p n = 0.25 k 1 , 0 γ 1 2 μ n 2 ,
μ n > 0 ,   n = 1 , 2 , 3 , , are the real roots of the following functional equation:
J 0 ( γ ε μ ) J 1 ( μ ) + K ε J 0 ( μ ) J 1 ( γ ε μ ) = 0 ,
where J k ( x ) denotes the Bessel functions of the first kind of the k-th order [21].
The temperature rise Θ ( z , t ) corresponding to the specific friction power q ( t ) in Equation (3) is searched on the basis of Duhamel’s formula [22].
Θ ( z , t ) = t 0 t q * ( t s ) Θ ^ ( z , s ) d s ,   0 < t t s ,
where Θ ^ ( z , t ) is the temperature rise in Equations (18)–(24) for constant specific friction power q ( t ) = q 0 . Taking into account the solutions to Equations (18) and (19) in Duhamel’s integral (Equation (25)), it was achieved that
Θ ( z , t ) = Λ e γ 1 z / 2 [ e γ 1 z / 2 ( 1 + γ ε K ε ) q * ( t ) + 4 γ ε n = 1 ϕ 1 ( z , μ n ) Ψ ( μ n ) G n ( t ) ] , z 0 , 0 t t s ,
Θ ( z , t ) = Λ e γ 2 z / 2 [ e γ 2 z / 2 ( 1 + γ ε K ε ) q * ( t ) + 4 γ ε n = 1 ϕ 2 ( z , μ n ) Ψ ( μ n ) G n ( t ) ] , z 0 , 0 t t s ,
where G n ( t ) is a derivative of the function G n ( t ) , which is determined as
G n ( t ) = 0 t q * ( t s ) e p n t d t , n = 1 , 2 , 3 ,
Substituting the temporal profile of the specific power of friction q * ( t ) in Equation (3) into Equation (28) yielded the following equation:
G n ( t ) = G n , 1 ( t ) 1 t s 0 G n , 2 ( t ) + t i t s 0 G n , 3 ( t ) , n = 1 , 2 , 3 , ,
where
G n , 1 ( t ) = 0 t p * ( t s ) e p n t d t ,   G n , 2 ( t ) = 0 t ( t s ) p * ( t s ) e p n t d t ,   G n , 3 ( t ) = 0 t [ p * ( t s ) ] 2 e p n t d t .
The calculations of integrals in Equation (30) taking into account the time profile of contact pressure p * ( t ) (1), give
G n , 1 ( t ) = p n 1 ( 1 e p n t ) + a n 1 ( e p n t e t / t i ) ,
G n , 2 ( t ) = t ( p n 1 a n 1 e t / t i ) p n 2 ( 1 e p n t ) a n 2 ( e p n t e t / t i ) ,
G n , 3 ( t ) = p n 1 ( 1 e p n t ) + 2 a n 1 ( e p n t e t / t i ) b n 1 ( e p n t e 2 t / t i ) ,
where
a n = p n t i 1 0 ,   b n = p n 2 t i 1 0 n = 1 , 2 , 3 ,
If for any n = k ,   k = 1 , 2 , , the equality p k = t i 1 ( a k = 0 ,   b k = t i 1 ) is true, then the integration of the Equation (30) gives
G k , 1 ( t ) = t i ( 1 e t / t i ) t e t / t i ,
G k , 2 ( t ) = t i [ t t i ( 1 e t / t i ) ] 0.5 t 2 e t / t i ,
G k , 3 ( t ) = t i ( 1 e 2 t / t i ) 2 t e t / t i .
On the other hand, for p k = 2 t i 1 ( a k = t i 1 ,   b k = 0 ) it was obtained that
G k , 1 ( t ) = 0.5 t i ( 1 e t / t i ) 2 ,
G k , 2 ( t ) = 0.5 t i [ t 0.5 t i ( 1 e 2 t / t i ) ] t i [ t t i ( 1 e t / t i ) ] e t / t i ,
G k , 3 ( t ) = 0.5 t i ( 1 e 2 t / t i ) 2 t i ( 1 e t / t i ) e t / t i + t e 2 t / t i .
Substituting the function G n , i ( t ) , i = 1 , 2 , 3 in Equations (31)–(33) into the right side of Equation (29) yields
G n ( t ) = ( 1 + t i t s 0 + 1 t s 0 p n ) ( 1 e p n t ) p n ( 1 + 2 t i t s 0 + 1 t s 0 a n ) ( e t / t i e p n t ) a n + + t i ( e 2 t / t i e p n t ) t s 0 b n t t s 0 ( 1 p n e t / t i a n ) ,   0 t t s ,   n = 1 , 2 , .
The searched derivative of the function G n ( t ) in Equation (41), meeting the conditions in Equation (34), has the following form:
G n ( t ) = ( 1 + t i t s 0 ) e p n t ( 1 e p n t ) t s 0 p n + ( 1 + 2 t i t s 0 + 1 t s 0 a n ) ( t i 1 e t / t i p n e p n t ) a n + + 1 t s 0 a n ( 1 t t i ) e t / t i t i ( 2 t i 1 e 2 t / t i p n e p n t ) t s 0 b n ,   0 t t s ,   n = 1 , 2 , .
Proceeding in a similar manner, from Equations (29) and (35)–(40), the derivative for p k = t i 1 was found.
G k ( t ) = t t s 0 ( 3 + t s 0 t i t 2 t i ) e t / t i + t i t s 0 ( 2 e 2 t / t i e t / t i 1 ) ,
That for p k = 2 t i 1 was also found.
G k ( t ) = ( 1 + 4 t i t s 0 ) ( e t / t i e 2 / t i ) t i 2 t s 0 ( 1 e 2 t / t i ) t t s 0 ( e t / t i + 2 e 2 t / t i ) .
Approaching p n 0 ( a n t i 1 , b n 2 t i 1 ), the limit of Equation (42) was found.
lim p n 0 G n ( t ) = 1 + t i t s 0 t t s 0 ( 1 + t i t s 0 ) e t / t i + t i t s 0 e 2 t / t i t i t s 0 ( 1 t t i ) e t / t i = = ( 1 e t / t i ) [ 1 t t s 0 + t i t s 0 ( 1 e t / t i ) ] = p * ( t ) [ 1 t t s 0 + t i t s 0 p * ( t ) ] = q * ( t ) , 0 t t s ,
where p * ( t ) and q * ( t ) are the dimensionless temporal profiles of pressure (Equation (1)) and specific friction power (Equation (3)), respectively, where the function q * ( t ) occurs beyond the sign of the sum in the solutions in Equations (26)–(28).
It should be noted that, at the initial time moment, from Equation (3), it follows q * ( 0 ) = 0 and, from Equation (42), taking into account Equation (34), it was found that
G n ( 0 ) = 1 + t i t s 0 + ( 1 + 2 t i t s 0 + 1 t s 0 a n ) ( t i 1 p n ) a n t i ( 2 t i 1 p n ) t s 0 b n + 1 t s 0 a n = = 1 + t i t s 0 1 2 t i t s 0 1 t s 0 a n + t i t s 0 + 1 t s 0 a n = 0 ,   p n t i 1 p n 2 t i 1   ,   n = 1 , 2 ,   .
If p k = t i 1 or p k = 2 t i 1 , then, from Equations (43) and (44), it follows that G k ( 0 ) = 0 . In this way, it was shown that the solution in Equations (26) and (27) meets the initial condition in Equation (16).
In the special case for t i 0 , when the pressure p ( t ) in Equation (1) attains the nominal value p 0 immediately, and the velocity V ( t ) in Equation (2) reduces linearly (braking with constant deceleration), the dimensionless temporal profile of the specific friction power q * ( t ) and function G n ( t ) in Equation (42) takes the following form:
q * ( t ) = 1 t t s 0 ,   G n ( t ) = e p n t ( 1 e p n t ) t s 0 p n ,   0 t t s 0 .
From Equations (47) and (48), it follows that q * ( 0 ) = 1 , G n ( 0 ) = 1 . This means that fulfillment of the initial condition in Equation (16) in this case is possible when the following equality is satisfied:
n = 1 ϕ 1 ( z , μ n ) Ψ ( μ n ) = 0.25 γ ε 1 + γ ε K ε e γ 1 z / 2 , z 0 ,   n = 1 ϕ 2 ( z , μ n ) Ψ ( μ n ) = 0.25 γ ε 1 + γ ε K ε e γ 2 z / 2 , z 0 ,
where functions ϕ l ( z , μ n ) , l = 1 , 2 , and Ψ ( μ n ) have the form in Equations (20)–(22). The validation of the summation of functional series in Equation (48) was performed numerically.

4. Dimensionless Form of Solution

The following denotes are introduced:
ζ = z a ,   τ = k 1 , 0 t a 2 ,   τ s = k 1 , 0 t s a 2 ,   τ s 0 = k 1 , 0 t s 0 a 2 ,   τ i = k 1 , 0 t i a 2 ,   γ l = γ l * a , l = 1 , 2 ,   Θ 0 = q 0 a K 1 , 0 ,   Θ * = Θ Θ 0 ,
where a = max { a 1 ,   a 2 } , a l , and l = 1 , 2 is the thickness of the friction pair element, which actively participates in the absorption of heat. This is the distance from the friction surface, on which the temperature is 5% of maximum values achieved on this surface [23].
a l = 3 k l , 0 t s ,   l = 1 , 2 .
Taking into account the denotes in Equation (49) in Equations (1)–(3), (20), and (42), and the solutions in Equations (26) and (27), the dimensionless temperature rise can be written in the following form:
Θ * ( ζ , τ ) = K 0 * γ 2 * e γ 1 * ζ / 2 [ e γ 1 * ζ / 2 ( 1 + γ ε K ε ) q * ( τ ) + 4 γ ε n = 1 ϕ 1 * ( ζ , μ n ) Ψ ( μ n ) G ( τ , μ n ) ] ,   ζ 0 , 0 τ τ s ,
Θ * ( ζ , τ ) = K 0 * γ 2 * e γ 2 * ζ / 2 [ e γ 2 * ζ / 2 ( 1 + γ ε K ε ) q * ( τ ) + 4 γ ε n = 1 ϕ 2 * ( ζ , μ n ) Ψ ( μ n ) G ( τ , μ n ) ] ,   ζ 0 ,   0 τ τ s ,
where
q * ( τ ) = p * ( τ ) [ 1 τ τ s 0 + τ i τ s 0 p * ( τ ) ] ,   p * ( τ ) = 1 e τ / τ i ,
ϕ 1 * ( ζ , μ n ) = J 1 ( γ ε μ n ) J 1 ( μ n e γ 1 * ζ / 2 ) ,   ϕ 2 * ( ζ , μ n ) = J 1 ( μ n ) J 1 ( γ ε μ n e γ 2 * ζ / 2 ) ,
G n ( τ ) = ( 1 + τ i τ s 0 ) e λ n τ ( 1 e λ n τ ) τ s 0 λ n + ( 1 + 2 τ i λ s 0 + 1 τ s 0 α n ) ( τ i 1 e τ / τ i λ n e λ n τ ) α n τ i ( 2 τ i 1 e 2 τ / τ i λ n e λ n τ ) τ s 0 β n + 1 τ s 0 α n ( 1 τ τ i ) e τ / τ i ,     α n = λ n τ i 1 0 ,   β n = λ n 2 τ i 1 0 ,
λ n = ( 0.5 γ 1 * μ n ) 2 ,   n = 1 , 2 , ,
τ s τ s 0 + 0.99 τ i ,   0 < τ i 0.3 τ s 0 .
Function Ψ ( μ n ) is given by Equation (21), and numbers μ n > 0 are the real roots of the functional Equation (24). From Equations (43) and (44), it follows that
G k ( τ ) = τ τ s 0 ( 3 + τ s 0 τ i τ 2 τ i ) e τ / τ i + τ i τ s 0 ( 2 e 2 τ / τ i e τ / τ i 1 ) ,   λ k = τ i 1 ,
G k ( τ ) = ( 1 + 4 τ i τ s 0 ) ( e t / t i e 2 τ / τ i ) τ i 2 τ s 0 ( 1 e 2 τ / τ i ) τ τ s 0 ( e τ / τ i + 2 e 2 τ / τ i ) ,   λ k = 2 τ i 1 .
Substituting ζ = 0 into Equations (51), (52), and (54), the dimensionless temperature rise on the contact surface can be written in the following form:
G k ( τ ) = ( 1 + 4 τ i τ s 0 ) ( e t / t i e 2 τ / τ i ) τ i 2 τ s 0 ( 1 e 2 τ / τ i ) τ τ s 0 ( e τ / τ i + 2 e 2 τ / τ i ) ,   λ k = 2 τ i 1 ,
where
ϕ * ( μ n ) ϕ 1 * ( 0 , μ n ) = ϕ 2 * ( 0 , μ n ) = J 1 ( γ ε μ n ) J 1 ( μ n ) .
In case of braking with constant deceleration ( τ i 0 ), from Equation (47), it can be obtained that
q * ( τ ) = 1 τ τ s 0 ,   G ( τ , μ n ) = e λ n   τ 1 τ s 0 λ n ( 1 e λ n   τ ) ,   0 τ τ s 0 .
It should be noted that the exact solution to the problem considering the contact scheme of friction for two semi-infinite bodies, made of homogeneous materials ( γ 1 = γ 2 = 0 ), with account of the time of contact pressure increase, was achieved in [19]. A special case of this solution—braking with constant deceleration—was investigated in [24].

5. Numerical Analysis

On the basis of the obtained exact solutions in Equations (51), (52), and (60), the calculations of the temperature generated due to friction in the disc–pad system during single braking were performed. Materials of the friction surfaces of elements were zirconium dioxide ( l = 1 ) and the other ceramic ( l = 2 ). With the distance from these surfaces deeper into the bodies, their thermal conductivity coefficients increased exponentially in accordance with Equation (10), reaching at the effective depths a l , l = 1 , 2 values corresponding to titanium and aluminum alloys, respectively. The thermal properties of the abovementioned materials are listed in Table 1.
Values of the remaining input parameters were as follows: A a = 0.442 10 2 m 2 , f = 0.27 , p 0 = 0.607   MPa , T 0 = 20   C , V 0 = 23.8   m   s 1 , and W 0 = 103.54   kJ [26]. From Equation (2), the braking time with constant deceleration was found t s 0 = 12   s and, next, the stop time t s = 12.49   s . This allowed determining from Equation (50) the effective depths of heat penetration a 1 = 5.556   mm and a 2 = 6.435   mm , as well as the value of the scaling parameter a = a 2 . According to the methodology, described in detail in [17], the dimensionless parameters of the material gradient were also established as γ 1 * = 1.28 and γ 2 * = 4.05 .
Isotherms of the temperature rise Θ ( z , t ) inside the elements of the friction pair are illustrated in Figure 2. The most heated ( Θ = 800 ÷ 943   C ) was a narrow, approximately 0.5   mm thick, near-surface area that appeared 3   s after start of braking. The lifetime of such a high-temperature area is 3   s . The friction surfaces of both elements were cooled down until the stop time moment. At the stop moment, the distance from the friction surface, where the noticeable temperature occurs in the disc was greater than in the pad.
Evolutions of the temperature rise Θ ( z , t ) during braking on the contact surface and inside the friction elements on different depths are presented in Figure 3. At the beginning of braking, the temperature on the friction surfaces z = 0 rapidly increases over time, achieving the maximum value Θ max = 943 C at the moment t max = 5   s . This is followed by a period of cooling of these surfaces until it stops. The temporal profiles of the temperature inside the disc and the pad also have a similar shape. However, the known “delay” effect is visible in the disc, which is that the time to reach the maximum temperature increases with the distance from the contact surface. At the same time, this effect is almost imperceptible. Noteworthy is also the process of rapid cooling on the friction surface of disc after reaching maximum Θ max ; at the stop moment, the temperature inside the disc is higher than on the surface. Again, this effect does not occur in the pad material.
Variations of the temperature during braking on the friction surfaces of disc and pad for different times of contact pressure increase are demonstrated in Figure 4. Extending the time of achieving the nominal value of pressure causes a drop of maximum temperature on the contact surface, while increasing the braking time. The effect of temperature drop with the growth of time of pressure increase is also presented in Figure 5.
The influence of dimensionless parameters of material gradients γ l * , l = 1 , 2 on the dimensionless maximum temperature Θ max * on the contact surface is illustrated in Figure 6. It shows that an increase in the core material volume fraction in selected FGMs (Ti-6Al-4V for disc and aluminum alloy for pad) causes a decrease in the maximum temperature in the brake. The biggest drop in Θ max * occurs when the gradient of the pad material γ 2 * is increased (Figure 6b). However, the highest values of Θ max * are reached for the friction pair in which one of the elements is entirely made of homogeneous material. These are zirconium dioxide ZrO 2 for the disc ( Θ max * = 995   C at γ 1 * = 0 , in Figure 6a) and the other ceramic for the pad ( Θ max * = 1340   C at γ 2 * = 0 , in Figure 6b).

6. Conclusions

The presented results are the continuation of an investigation from a previous article of the authors [17], in which, in the dimensionless form, a comparative, qualitative analysis was performed in order to study the influence of gradient of FGMs on the temperature during frictional heating under uniform sliding. However, in this paper the mathematical model was derived to determine the temperature field in a disc–pad system during single braking. An important and unique feature of this model was its taking into account of the time-dependent pressure and velocity for friction elements, made of functionally graded materials with exponentially changing conductivity coefficients with thickness. The proposed model allows for a quick assessment of the brake temperature mode depending on the operational parameters, such as the time of contact pressure increase and the value of the gradient of the friction materials. The analysis was performed in the dimensional form. The friction surfaces of the materials were ceramic, and their cores were titanium alloys (disc) and aluminum alloys (pad). It was established that extending the time of pressure increase causes significant extending of the braking time and, thus, extending of the braking distance. The maximum temperature reached on the friction surfaces drops when the parameters of material gradients are increased.
Application of the proposed model has some limitations, resulting from the simplifying assumptions made, especially the use of only an exponential function to describe the thermal conductivity changes in FGMs. In further research, it is planned to include in the formulation of the boundary value problem of heat conduction, as well as the thermal resistance on the contact surface of the disc and the pad (imperfect thermal contact of friction), and to adapt the obtained exact solution to determine the temperature of the brake during a repeated short-term mode of braking.
As shown in the results of the numerical analysis presented in this article, the maximum temperature achieved even with a single braking is quite high. With such a temperature, the necessary problem is to develop a model that takes into account the thermal sensitivity of the materials. Some steps toward implementing the exact solutions of linear problems for homogeneous materials to take into account their thermal sensitivity have already been made for a single [27] and a repetitive short-term [28] braking modes. On the basis of this methodology, the development of appropriate models for FGM has begun.
Another problem caused by high temperatures is a reduction in the strength of the material, especially when the temperature exceeds the melting point of the aluminum alloy. Investigations of the strength were not the subject of this article, but they should also be considered in the future.

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, P.Z.; writing—review and editing, A.Y. and K.T.; visualization and figure preparation, P.Z.; supervision, A.Y.; project administration, K.T.; funding acquisition, K.T. All authors 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 Center, Poland and project financing through the program of the Minister of Education and Science of Poland “Regional Initiative of Excellence” in 2019–2022, project No. 011/RID/2018/19 (amount of 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 l Effective depth of heat penetration ( m )
A a Area of the nominal contact region ( m 2 )
c l Specific heat ( J   kg 1 K 1 )
f Coefficient of friction (dimensionless)
J k ( ) Bessel functions of the first kind of the k-th order
k l Thermal diffusivity ( m 2 s 1 )
K l Thermal conductivity ( W   m 1 K 1 )
p Contact pressure ( Pa )
p 0 Nominal value of the contact pressure ( Pa )
q Specific power of friction ( W   m 2 )
q 0 Nominal value of the specific power of friction ( W   m 2 )
t Time ( s )
t i Time of the contact pressure increase ( s )
t s 0 Stop time at braking with constant deceleration ( s )
t s Stop time ( s )
T Temperature ( C )
T a Initial (ambient) temperature ( C )
V Velocity ( m   s 1 )
V 0 Initial velocity ( m   s 1 )
W 0 Initial kinetic energy of the system ( J )
x ,   y ,   z Spatial coordinates ( m )
lower l Number of the main ( l = 1 ) and frictional ( l = 2 ) elements of the friction pair
γ l Parameter of material gradient ( m - 1 )
γ l * Parameter of material gradient (dimensionless)
Θ l Temperature rise ( C )
Θ l * Temperature rise (dimensionless)
Θ 0 Temperature scaling factor ( C )
ρ l Density ( kg   m 3 )
τ Time (dimensionless)
τ i Time of contact pressure increase (dimensionless)
τ s 0 Braking time at constant deceleration (dimensionless)
τ s Braking time (dimensionless)
ζ Spatial coordinate in axial direction (dimensionless)

References

  1. Govindaraju, M.; Megalingam, A.; Murugasan, J.; Vignesh, R.V.; Kota, P.K.; Ram, A.S.; Lakshana, P.; Kumar, V.N. Investigations on the tribological behavior of functionally gradient iron-based brake pad material. Proc. Inst. Mech. Eng. Part C J. Mech. Eng. Sci. 2020, 234, 2474–2486. [Google Scholar] [CrossRef]
  2. Swaminathan, K.; Sangeetha, D. Thermal analysis of FGM plates—A critical review of various modeling techniques and solution methods. Compos. Struct. 2017, 160, 43–60. [Google Scholar] [CrossRef]
  3. Afsar, A.; Go, J. Finite element analysis of thermoelastic field in a rotating FGM circular disk. Appl. Math. Model. 2010, 34, 3309–3320. [Google Scholar] [CrossRef]
  4. Shahzamanian, M.; Bin Sahari, B.; Bayat, M.; Mustapha, F.; Ismarrubie, Z.N. Finite element analysis of thermoelastic contact problem in functionally graded axisymmetric brake disks. Compos. Struct. 2010, 92, 1591–1602. [Google Scholar] [CrossRef]
  5. Shahzamanian, M.; Bin Sahari, B.; Bayat, M.; Ismarrubie, Z.N.; Mustapha, F. Transient and thermal contact analysis for the elastic behavior of functionally graded brake disks due to mechanical and thermal loads. Mater. Des. 2010, 31, 4655–4665. [Google Scholar] [CrossRef]
  6. Hosseini, T.P.; Talebi, M. Stress and temperature distribution study in a functionally graded brake disk. Int. J. Automot. Eng. 2012, 2, 172–179. [Google Scholar]
  7. Yaghoobi, M.P.; Ghannad, M. An analytical solution for heat conduction of FGM cylinders with varying thickness subjected to non-uniform heat flux using a first-order temperature theory and perturbation technique. Int. Commun. Heat Mass Transf. 2020, 116, 104684. [Google Scholar] [CrossRef]
  8. Jang, Y.H.; Ahn, S.-H. Frictionally-excited thermoelastic instability in functionally graded material. Wear 2007, 262, 1102–1112. [Google Scholar] [CrossRef]
  9. Lee, S.W.; Jang, Y.H. Frictionally excited thermoelastic instability in a thin layer of functionally graded material sliding between two half-planes. Wear 2009, 267, 1715–1722. [Google Scholar] [CrossRef]
  10. Lee, S.W.; Jang, Y.H. Effect of functionally graded material on frictionally excited thermoelastic instability. Wear 2009, 266, 139–146. [Google Scholar] [CrossRef]
  11. Mao, J.-J.; Ke, L.-L.; Yang, J.; Kitipornchai, S.; Wang, Y.-S. Thermoelastic instability of functionally graded materials with interaction of frictional heat and contact resistance. Mech. Based Des. Struct. Mach. 2017, 46, 139–156. [Google Scholar] [CrossRef]
  12. 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]
  13. Kulchytsky-Zhyhailo, R.; Bajkowski, A.S. Axisymmetrical problem of thermoelasticity for halfspace with gradient coating. Int. J. Mech. Sci. 2016, 106, 62–71. [Google Scholar] [CrossRef]
  14. Martínez-Pañeda, E. On the Finite Element Implementation of Functionally Graded Materials. Materials 2019, 12, 287. [Google Scholar] [CrossRef] [Green Version]
  15. Mao, J.-J.; Ke, L.-L.; Yang, J.; Kitipornchai, S.; Wang, Y.-S. Thermoelastic instability of functionally graded coating with arbitrarily varying properties considering contact resistance and frictional heat. Appl. Math. Model. 2017, 43, 521–537. [Google Scholar] [CrossRef] [Green Version]
  16. Polajnar, M.; Kalin, M.; Thorbjornsson, I.; Thorgrimsson, J.; Valle, N.; Botor-Probierz, A. Friction and wear performance of functionally graded ductile iron for brake pads. Wear 2017, 382-383, 85–94. [Google Scholar] [CrossRef]
  17. 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] [PubMed]
  18. Chichinadze, A.V. Polymers in Friction Assembles of Machines and Devices: A Handbook; Allerton Press Inc: New York, NY, USA, 1984. [Google Scholar]
  19. 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]
  20. Yevtushenko, A.; Kuciej, M.; Topczewska, K. Some theoretical model for determining the temperature field of a multi-disk brake. Adv. Mech. Eng. 2020, 12. [Google Scholar] [CrossRef]
  21. 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, WA, USA, 1964. [Google Scholar]
  22. Özis̨ik, N.M. Heat conduction; John Wiley: New York, NY, USA, 1993. [Google Scholar]
  23. Chichinadze, A.V.; Braun, E.D.; Ginzburg, A.G.; Ignat’eva, E.V. Calculation, Testing and Selection of Friction Couples; Nauka: Moscow, Russia, 1979. (In Russian) [Google Scholar]
  24. Yevtushenko, A.; Kuciej, M. Frictional heating during braking in a three-element tribosystem. Int. J. Heat Mass Transf. 2009, 52, 2942–2948. [Google Scholar] [CrossRef]
  25. Tian, J.; Jiang, K. Heat conduction investigation of the functionally graded materials plates with variable gradient parameters under exponential heat source load. Int. J. Heat Mass Transf. 2018, 122, 22–30. [Google Scholar] [CrossRef]
  26. Chichinadze, A.V.; Kozhemyakina, V.D.; Suvorov, A.V. Method of temperature-field calculation in model ring specimens during bilateral friction in multidisc aircraft brakes with the IM-58-T2 new multipurpose friction machine. J. Frict. Wear 2010, 31, 23–32. [Google Scholar] [CrossRef]
  27. 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]
  28. 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]
Figure 1. Scheme of the frictional heating in the disc–pad system.
Figure 1. Scheme of the frictional heating in the disc–pad system.
Materials 14 06241 g001
Figure 2. Isotherms of the temperature rise Θ ( z , t ) in the disc and the pad at t i = 0.5   s .
Figure 2. Isotherms of the temperature rise Θ ( z , t ) in the disc and the pad at t i = 0.5   s .
Materials 14 06241 g002
Figure 3. Evolutions of the temperature rise Θ ( z , t ) during braking at t i = 0.5   s for different distances from the friction surface: (a) the disc; (b) the pad.
Figure 3. Evolutions of the temperature rise Θ ( z , t ) during braking at t i = 0.5   s for different distances from the friction surface: (a) the disc; (b) the pad.
Materials 14 06241 g003
Figure 4. Evolutions of the temperature rise Θ ( 0 , t ) during braking for different values of the time t i of contact pressure increase.
Figure 4. Evolutions of the temperature rise Θ ( 0 , t ) during braking for different values of the time t i of contact pressure increase.
Materials 14 06241 g004
Figure 5. Dependence of the maximum temperature rise Θ max on the time t i of contact pressure increase.
Figure 5. Dependence of the maximum temperature rise Θ max on the time t i of contact pressure increase.
Materials 14 06241 g005
Figure 6. Dependence of the maximum temperature rise Θ max at t i = 0.5   s on the dimensionless gradient of material: (a) γ 1 * for γ 2 * = 4.05 ; (b) γ 2 * for γ 1 * = 1.28 .
Figure 6. Dependence of the maximum temperature rise Θ max at t i = 0.5   s on the dimensionless gradient of material: (a) γ 1 * for γ 2 * = 4.05 ; (b) γ 2 * for γ 1 * = 1.28 .
Materials 14 06241 g006
Table 1. Thermophysical properties of the FGM components [15,25].
Table 1. Thermophysical properties of the FGM components [15,25].
Element SubscriptMaterial Thermal   Conductivity   K   [ Wm 1 K 1 ] Thermal   Diffusivity   k   ×   10 6   [ m 2 s 1 ]
l = 1 ZrO22.090.86
Ti-6Al-4V7.53.16
l = 2 ceramic31.15
aluminum alloy17367.16
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 Effect of Functionally Graded Materials on Temperature during Frictional Heating at Single Braking. Materials 2021, 14, 6241. https://doi.org/10.3390/ma14216241

AMA Style

Yevtushenko A, Topczewska K, Zamojski P. The Effect of Functionally Graded Materials on Temperature during Frictional Heating at Single Braking. Materials. 2021; 14(21):6241. https://doi.org/10.3390/ma14216241

Chicago/Turabian Style

Yevtushenko, Aleksander, Katarzyna Topczewska, and Przemysław Zamojski. 2021. "The Effect of Functionally Graded Materials on Temperature during Frictional Heating at Single Braking" Materials 14, no. 21: 6241. https://doi.org/10.3390/ma14216241

APA Style

Yevtushenko, A., Topczewska, K., & Zamojski, P. (2021). The Effect of Functionally Graded Materials on Temperature during Frictional Heating at Single Braking. Materials, 14(21), 6241. https://doi.org/10.3390/ma14216241

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