Next Article in Journal
Effect of SiC Abrasive Blasting Parameters on the Quality of the Ceramic and Ni-Cr Dental Alloy Joint
Next Article in Special Issue
Temperature in the Friction Couple Consisting of Functionally Graded and Homogeneous Materials
Previous Article in Journal
Effects of Fixture Configurations and Weld Strength Mismatch on J-Integral Calculation Procedure for SE(B) Specimens
Previous Article in Special Issue
An Effect of a Carbon-Containing Additive in the Structure of a Friction Material on Temperature of the Wet Clutch Disc
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Influence of Thermal Sensitivity of Functionally Graded Materials on Temperature during 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.
Submission received: 23 December 2021 / Revised: 14 January 2022 / Accepted: 24 January 2022 / Published: 26 January 2022

Abstract

:
The model of the frictional heating process during single braking to determine the temperature of the functionally graded friction elements with an account of the thermal sensitivity of materials was proposed. The basis of this model is the exact solution of the one-dimensional thermal problem of friction during braking with constant deceleration. The formulas approximating the experimental data of the temperature dependencies of properties of the functionally graded materials (FGMs) were involved in the model to improve the accuracy of the achieved results. A comparative analysis was performed for data obtained for temperature-dependent FGMs and the corresponding data, calculated without consideration of thermal sensitivity. The results revealed that the assumption of thermal stability of FGMs during braking may cause a significant overestimation of temperature of the friction pair elements.

1. Introduction

During intensive braking, the volume temperature of the disc braking system may be higher than 450 °C [1] and the maximum temperature on the friction surfaces of the pad and the disc during single braking may even reach a level above 1000 °C [2]. In such severe conditions, the thermal and mechanical properties of materials may highly differ from the initial, reached at the ambient temperature. Therefore, in order to improve the theoretical analysis of the thermoelastic behavior of the braking systems, it is necessary to develop mathematical models taking into consideration the thermal sensitivity of friction materials. However, the introduction of the temperature-dependent properties in formulation of the thermal problems of friction leads to nonlinearity, so most of the published analyses have been performed using numerical methods, especially the finite element method [3,4]. One of the alternative techniques used to develop such nonlinear models of frictional heating is linearization by means of the Kirchhoff substitution [5]. This method relies on the reduction of the originally nonlinear heat conduction equation to the linear one. However, it works this way only for materials with simple nonlinearity, which means that their thermal conductivity and specific heat capacity are temperature-dependent, but the thermal diffusivity remains constant [6]. For materials with arbitrary nonlinearity, only the partial linearization by the Kirchhoff substitution of such a problem is possible; as a result, another nonlinear problem is obtained for which the method of solving is known [7]. The Kirchhoff transform has a similar effect in the heat conduction problems formulated for solids with simple nonlinear thermosensitivity under complex heat exchange. Some analytic–numerical methods for the solution of such problems have been proposed in the study [8]. Another technique to take into consideration the thermal sensitivity of materials is the method of successive approximations (iterations), in which the solution of the corresponding linear problem is adopted as the initial approximation, and then the solution found in the previous step is corrected. An iteration algorithm to solve the one-dimensional problem of heat conduction at braking has been proposed in the article [9]. Most models of the frictional heating process taking into account of the thermal sensitivity effect have been developed only for homogeneous materials. Modern friction materials often have a non-uniform, complex internal structure with a changing composition, microstructure, or porosity across the volume of material, such as functionally graded materials (FGMs), which are characterized by smooth variations of properties as a function of position along certain direction. In the case of devices operating at elevated temperatures, including braking systems, FGMs are primarily used in order to obtain high temperature resistance on the friction surface by dissipating heat from it to the inside of the element while maintaining low wear. FGMs of this type are usually two-component. Their friction surface is usually made of metal-ceramic, and the metal opposite surface (core) should have high thermal conductivity. The change of properties in the direction perpendicular to the friction surface is described by continuous functions, usually power or exponential. In the case of the latter, the material gradient parameters are responsible for the speed of transition from one material to another.
The problem of wear of an FGM strip with an account of the heating on sliding contact from friction has been considered in the study [10]. The exact solution of the problem was obtained with the help of the integral Laplace transform technique. It was assumed that the shear modulus is described by means of the function of the vertical coordinate. A comprehensive review of the literature concerning the thermal contact problems of frictional heating for functionally graded materials was provided in our previous article [11]. So far, investigations of the transient heat conduction in FGMs are limited, and most of them have ignored the temperature dependence of the material properties. Therefore, in general, those models are adequate only for relatively low temperatures in an FGM or the materials with insignificant thermal sensitivity. To accurately describe the thermomechanical behavior of FGMs, the temperature dependence of the material properties should be considered. The heat conduction problems formulated for FGMs with non-uniform spatially distributed and temperature-dependent properties are highly nonlinear. Nevertheless, several studies concerning such problems taking into consideration the thermal sensitivity of FGMs can be found, but most of them are solved by means of numerical or semi-analytical methods. The finite element method has been adopted in the paper [12], to perform the nonlinear transient thermal stress analysis of a thick-walled FGM cylinder with temperature-dependent material properties. Another nonlinear transient heat transfer and thermoelastic stress in thermosensitive functionally graded cylinder have been investigated using the Hermitian transfinite element method in the study [13]. The results showed that the effect of thermal sensitivity of materials has a significant influence on the thermal behavior of friction systems.
An analytical approach to solve the one-dimensional transient heat conduction problem for functionally graded materials with temperature-dependent properties has been presented in the article [14]. As for the analytical treatment, the temperature and thermal stress solutions have been obtained in approximate forms for a simplified, homogeneous, multi-layered model of materials. They concluded that the temperature dependence of the material properties is one of the most important factors in the accurate evaluation of temperature and stress distributions [14]. A similar multi-layered model was used to formulate another thermal problem of friction for a thermally sensitive FGM plate in the paper [15]. The authors made an attempt to optimize the functionally graded structure in order to enhance their thermal performance. The proper manufacturing process allows the design of an FGM according to the engineering demands by intentionally setting a specific distribution of the properties. A hybrid genetic algorithm has been developed for the optimization of the FGM composition with temperature-dependent material properties, in order to minimize the thermal stresses under steady-state thermal loads [15]. The optimum composition profile of the functionally graded materials for wide temperature ranges was also studied in the article [16]. The thermoelastic problem for functionally graded material with temperature-dependent properties was considered by means of the perturbation method. Additionally, the crack propagation path was predicted by introducing the fracture mechanics analysis. It was concluded that the proper selection of an FGM gradient can lead to a significant decrease in thermal stresses [16]. A transient thermoelastic behavior of the functionally graded plate with temperature-dependent properties due to a thermal shock was considered in the paper [17]. The temperature and thermal stress distributions in the Cu-W functionally graded composite were found by means of the semi-analytical micromechanical model.
The aim of this study was to investigate the influence of FGMs thermal sensitivity on the distribution of temperature in a disc brake system. This study is a continuation of our previous articles [11,18], which concern the transient thermal problem of friction under uniform sliding and during single braking with an exponential increase in the contact pressure. Due to the appearance of a high temperature level, there is a demand to improve the results by involving the variations of material properties dependent on the actual temperature, since the thermal sensitivity effect is particularly manifested in a high temperature range. In this article, the braking with constant deceleration is considered, when the nominal pressure is reached immediately at the beginning of the process, since the increase in the time of contact pressure growth causes a drop in the achieved temperature.

2. Statement to the Problem

To develop an analytical model of frictional heating process in the braking system, the following assumptions were taken into account:
  • The braking process with constant deceleration is considered;
  • At the initial time moment, the temperature of a brake is equal to the ambient temperature Ta;
  • In the heat conduction equation, only the change in temperature gradient in the perpendicular direction to the disc-pad contact surfaces is taken into consideration;
  • The thermal contact on the friction surfaces is perfect, i.e., the temperatures of its contact surfaces are equal, and the sum of frictional heat fluxes intensities, acting along the normal direction to the contact surface to the insides of the elements equal to the specific friction power;
  • Due to the symmetry of the system with respect to the mid plane of the disc, when determining the brake temperature, the contact of one pad and a disc with half of its thickness is considered;
  • The pads and the disc are made of two-component thermally sensitive functionally graded materials, in such a way that their friction surfaces are materials with low thermal conductivity (i.e., cermet), while the core materials are characterized by higher thermal conductivity (titanium alloys, aluminum, etc.);
  • The thermal conductivity of the disc and pads materials increases exponentially with the distance from the contact surface;
  • The whole initial kinetic energy of the vehicle is transformed into heat during braking, neglecting the small part of energy associated with wear on the contact surfaces of the disc and pads;
Based on the assumptions (1)–(5), in order to determine the temperature of the disc-pad system, the scheme of sliding with linearly decreasing velocity of two semi-spaces z 0 (disc) and z 0 (pad) has been adopted. Initiated by the frictional heating temperature field of such a system at a given time instant t 0 depends only on the distance from the friction surface in a perpendicular direction—independent variable z : T = T ( z , t ) .
According to the assumption (6), the thermophysical properties of a friction pair are functions of temperature T :
K l , m = K l , m ( T ) , c l , m = c l , m ( T ) , ρ l , m = ρ l , m ( T ) ,
where K l , m , c l , m and ρ l , m —thermal conductivity, specific heat capacity and density of the first ( m = 1 ) and second ( m = 2 ) component of the materials of the disc ( l = 1 ) and pad ( l = 2 ) , respectively. Corresponding values at the initial system temperature T = T 0 are marked as follows:
K l , m ( 0 ) K l , m ( T 0 ) , c l , m ( 0 ) c l , m ( T 0 ) , ρ l , m ( 0 ) ρ l , m ( T 0 ) .
According to the mixture law, the effective specific heat capacities and densities were also determined:
c l ( 0 ) = c l , 2 ( 0 ) V c + ( 1 V c ) c l , 1 ( 0 ) ,   ρ l ( 0 ) = ρ l , 2 ( 0 ) V ρ + ( 1 V ρ ) ρ l , 1 ( 0 ) ,
where V c , V ρ —volume fractions of the phases c l , m ( 0 ) and ρ l , m ( 0 ) , l = 1 , 2 , m = 1 , 2 , respectively.
Based on the assumption (7), the effective thermal conductivities K l , l = 1 , 2 of the disc and pad materials were established from the equations:
K 1 ( z ) = K 1 , 1 e γ l z , 0 z a ,   K 2 ( z ) = K 2 , 1 e γ 2 z , a z 0 ,
where
γ l = γ l * a ,   γ l * = ln ( K l , 2 ( 0 ) K l , 1 ( 0 ) ) ,
a = max { a 1 , a 2 } , a l = 3 k l ( 0 ) t s ,
k l ( 0 ) = K l , 1 ( 0 ) c l ( 0 ) ρ l ( 0 ) ,
and t s —stop time, and parameters a l , l = 1 , 2 (6) are the thicknesses of the subsurface layers actively participating in heat absorption in the disc and pads, respectively (the so-called effective depth of heat transfer [19]). During braking with constant deceleration, the specific friction power decreases linearly from the nominal value q 0 to zero [20]:
q ( t ) = q 0 q * ( t ) , q 0 = f p 0 V 0 , q * ( t ) = 1 t t s 1 ,   0 t t s ,
t s = W 0 Q 0 1 ,   Q = q 0 A a , A a = 0.5 β ( R e 2 R i 2 ) ,
where A a —nominal area of the contact between the pad and the disc; f —friction coefficient; p 0 —nominal pressure; Q 0 —nominal friction power; 0 β 2 π —nominal friction power; —cover angle of the pad; R i and R e —respectively, the internal and external radii of the pads; V 0 , W 0 —the initial velocity and kinetic energy of the system, respectively. The latter, according to assumption (8), is equal to the total work of friction.
In order to solve the above-formulated nonlinear problem, we will use the idea of adapting an appropriate solution of the linear problem of thermal friction. This approach in the case of homogeneous materials was used in the studies [9,21].

3. Solution with Temperature-Independent FGMs Properties

The key element of the proposed approach is the precise solution of the linear thermal problem of friction during braking with constant deceleration. In the case of FGMs, such a solution for the above-adopted scheme of two sliding semi-spaces for the specific friction power q ( t ) (8) and (9) can be written in the form [18]:
T ( z , t ) = T 0 + Θ ( z , t ) ,   0 t t s ,
Θ ( 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 ,  
Θ ( z , t ) = Λ e γ 1 z / 2 [ e γ 2 z / 2 ( 1 + γ ε K ε ) q * ( t ) + 4 γ ε n = 1 φ 2 ( z , μ n ) Ψ ( μ n ) G n ( t ) ] ,   z 0 ,
G n ( t ) = e p n t ( 1 e p n t ) p n t s , p n = 1 4 k 1 ( γ 1 μ n ) 2 ,
φ 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 * ( k * ) 1 / 2 ,   γ ε = γ * ( k * ) 1 / 2 ,
Λ = q 0 γ 2 K 2 , 1 ( 0 ) ,   K * = K 1 , 1 ( 0 ) K 2 , 1 ( 0 ) ,   k * = k 1 ( 0 ) k 2 ( 0 ) ,   γ * = γ 1 γ 2 ,
where μ n > 0 , n = 1 , 2 , 3 , , are the real roots of the functional equation:
J 0 ( γ ε μ n ) J 1 ( μ n ) + K ε J 0 ( μ n ) J 1 ( γ ε μ n ) = 0 .
J k ( x ) —are the Bessel functions of the first kind of the kth order [22].
The temperature of the friction surfaces of both elements, in accordance with the assumption (4) of their perfect thermal contact of friction, should be the same. Substituting z = 0 in Equations (10)–(12) and (14), the following were obtained:
T ( t ) T ( 0 ± , t ) = T 0 + Θ ( t ) ,   0 t t s ,
Θ ( t ) Θ ( 0 ± , t ) = Λ [ q * ( t ) ( 1 + γ ε K ε ) + 4 γ ε n = 1 φ ^ ( μ n ) Ψ ( μ n ) G n ( t ) ] ,   0 t t s ,
where
φ ^ ( μ n ) = J 1 ( γ ε μ n ) J 1 ( μ n ) .
Dimensionless variables and parameters were introduced:
ζ = z a ,   τ = k 1 t a 2 ,   τ s = k 1 t s a 2 ,   Θ 0 = q 0 a K 1 , 1 ( 0 ) ,   Θ * = Θ Θ 0 ,
where parameters a and q 0 were determined accordingly from Formulas (6) and (8). Taking into account the indications (22) in Formulas (11)–(14), the dimensionless temperature rise of the friction pair elements can be presented in the 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:
φ 1 * ( ζ , μ n ) = J 1 ( γ ε μ n ) J 1 ( μ n e γ 1 * ζ / 2 ) ,   φ 2 ( ζ , μ n ) = J 1 ( μ n ) J 1 ( γ ε μ n e γ 2 * ζ / 2 ) ,
G n ( τ ) = e λ n τ 1 λ n τ s ( 1 e λ n τ ) ,   λ n = 1 4 ( γ 1 * μ n ) 2 ,
and the remaining functions as well as parameters are given by Formulas (15)–(18). Substituting ζ = 0 in Formulas (23)–(25), the dimensionless rise of the temperature on the friction surfaces was obtained:
Θ * ( τ ) Θ * ( 0 ± , τ ) = K 0 * γ 2 * [ q * ( τ ) ( 1 + γ ε K ε ) + 4 γ ε n = 1 φ ^ ( μ n ) Ψ ( μ n ) G n ( τ ) ] ,   0 τ τ s ,
Based on Fourier’s law, the intensity of heat fluxes directed along the normal to the contact surface z = 0 towards the insides of the friction pair elements were defined:
q l ( t ) = ( 1 ) l K l , 1 ( 0 ) T ( z , t ) z | z = 0 ± ,   0 t t s ,   l = 1 , 2 .
Taking into account the indications (22) dimensionless intensities of heat fluxes q l * = q l q 0 1 , l = 1 , 2 were written as:
q 1 * ( τ ) = Θ * ( ζ , τ ) ζ | ζ = 0 + ,   q 2 * ( τ ) = Θ * ( ζ , τ ) K 0 * ζ | ζ = 0 ,   0 τ τ s .
After differentiating the solution (23)–(26) with respect to the variable ζ and subsequent substitution of the found derivatives to the right side of Formula (30), the following was found:
q 1 * ( τ ) = γ ε K ε ( 1 + γ ε K ε ) q * ( τ ) + 2 γ ε K ε n = 1 φ ˜ ( μ n ) Ψ ˜ ( μ n ) G n ( τ ) ,   0 τ τ s ,
q 2 * ( τ ) = 1 ( 1 + γ ε K ε ) q * ( τ ) + 2 n = 1 φ ˜ ( μ n ) Ψ ˜ ( μ n ) G n ( τ ) ,   0 τ τ s ,
where:
φ ˜ ( μ n ) = J 0 ( γ ε μ n ) J 1 ( μ n ) ,   Ψ ˜ ( μ n ) = μ n 1 Ψ ( μ n ) ,
and functions Ψ ( μ n ) and G n ( τ ) can be determined from Equations (15) and (26), respectively.
It should be noted that in the case of homogeneous materials ( γ i 0 , i = 1 , 2 ) of the disc and pads, the dimensionless temperature rise during braking with a constant deceleration has the form [23]:
Θ * ( ζ , τ ) = 2 K * τ ( 1 + K ε ) { ierfc ( ζ 2 τ ) τ τ s [ ( 1 + ζ 2 6 τ ) ierfc ( ζ 2 τ ) e ζ 2 4 τ 3 π ] } , ζ 0 ,   0 τ τ s ,
Θ * ( ζ , τ ) = 2 K * τ ( 1 + K ε ) { ierfc ( ζ 2 k 0 * τ ) τ τ s [ ( 1 + ζ 2 k 0 * 6 τ ) ierfc ( ζ 2 k 0 * τ ) e ζ 2 k 0 * 4 τ 3 π ] } , ζ 0 ,   0 τ τ s ,
where ierfc ( x ) = π 1 / 2 e x 2 x erfc ( x ) , erfc ( x ) = 1 erf ( x ) , erf ( x ) —Gaussian error function. For ζ = 0 from Equations (33) and (34), the known solution of Fazekas was obtained [24]:
Θ * ( τ ) = 2 K * ( 1 + K ε ) τ π ( 1 2 τ 3 τ s ) ,   0 τ τ s .

4. Volume Temperature

With the given input parameters, solutions (19)–(27) make it possible to find the space–time distribution of the temperature inside and its evolution on the friction surfaces of the pad and disc, made of thermally insensitive FGMs. In order to take into account the thermal sensitivity of materials determining the temperature of the braking system using the above-mentioned solutions, the thermal conductivities K l , m ( 0 ) , specific heat capacities c l , m ( 0 ) and densities ρ l , m ( 0 ) at the initial temperature T 0 were replaced with corresponding values K l , m ( ϑ l ) , c l , m ( ϑ l ) and ρ l , m ( ϑ l ) found from Equations (1) and (2) for the volume temperature of the pad and the disc during braking [2,9]:
ϑ l = T 0 + ϑ ^ l ,   l = 1 , 2 ,
where:
ϑ ^ l = 2 α l W 0 3 G l c l ( 0 ) ,
G l = A a a l ρ l ( 0 ) ,
a l —the effective depths of heat penetration (6), α 1 = α , α 2 = 1 α , 0 α 1 —heat partition ratio. Based on Formulas (30)–(32), the heat partition ratio was calculated from the formula:
α q 1 ( t ) q ( t ) γ ε K ε 1 + γ ε K ε .

5. Numerical Analysis

The calculations were performed for the friction pair, one element of which was made of aluminium oxide Al2O3 (friction surface) and cooper Cu (core) [25]. The friction surface and core of the second element are manufactured of zirconium dioxide ZrO2 and titanium alloy Ti-6Al-4V [14]. The temperature-dependent properties of these materials are as follows:
Al2O3 [26,27,28]
K 1 , 1 ( T ) = 39.717 0.130 T + 4.463 10 4 T 2 2.836 10 7 T 3 + 1.941 10 10 T 4 ,
c 1 , 1 ( T ) = 680.72 + 2.432 T 0.53 10 2 T 2 + 0.6 10 5 T 3 0.4 10 8 T 4 + 10 12 T 5 ,
ρ 1 , 1 ( T ) = 3992.2 0.062 T 0.6 10 4 T 2 + 0.4 10 7 T 3 0.9 10 11 T 4 ,
Cu [17,29]
K 1 , 2 ( T ) = 31.985 + 0.0099 T 0.1 10 5 T 2 ,
c 1 , 2 ( T ) = 523.3 + 1.4726 T 0.0024 T 2 + 0.2 10 5 T 3 0.5 10 9 T 4 ,
ρ 1 , 2 ( T ) = 492.45 0.01 T 0.1 10 5 T 2 ,
ZrO2 [27,30,31]
K 2 , 1 ( T ) = 1.9365 + 0.7 10 4 T + 0.5 10 6   T 2 0.2 10 9 T 3 ,
c 2 , 1 ( T ) = 437.96 + 0.7767 T 0.17 10 2 T 2 ,
ρ 2 , 1 ( T ) = 6104.6 0.1212 T 0.4 10 4 T 2 + 0.3 10 7 T 3 0.1 10 10 T 4 ,
Ti-6Al-4V [32,33]
K 2 , 2 ( T ) = 6.6926 + 8.9177 10 3 T + 6.8432 10 6 T 2 ,
c 2 , 2 ( T ) = 529.9316 + 0.4154 T 4.01646 10 4 T 2 + 1.6364 10 7 T 3 ,
ρ 2 , 2 ( T ) = 4434 0.1088 T 0.8 10 4 T 2 + 10 7 T 3 0.6 10 10 T 4 .
Graphs of dimensionless functions K l , m * = K l , m ( T ) / K l , m ( 0 ) , c l , m * = c l , m ( T ) / c l , m ( 0 ) and ρ l , m * = ρ l , m ( T ) / ρ l , m ( 0 ) are illustrated in the Figure 1, Figure 2 and Figure 3.
The calculations were performed according to the following scheme:
(1)
the values of the input parameters were given (Table 1), and then from Equations (8) and (9) the area of the nominal contact was calculated A a = 0.0022 m 2 , specific friction power q 0 = 3.87 MW m 2 , friction power Q 0 = 8510 W and stop time t s = 12.1 s ;
(2)
using the dependencies (40)–(51) the materials properties K l , m ( 0 ) , c l , m ( 0 ) and ρ l , m ( 0 ) , l , m = 1 , 2 at the initial temperature T 0 = 20 C were established (Table 2);
(3)
the effective values of: the specific heat c l ( 0 ) , density ρ l ( 0 ) , thermal diffusivity k l ( 0 ) , the effective depths of heat penetration a l and the dimensionless gradient parameters of materials γ l * , l = 1 , 2 were found from Equations (3) and (5)–(7). Then, the dimensionless parameters K ε and γ ε were determined from the Formulas (16) and (17), and also the weight G l and heat partition ratios α l , l = 1 , 2 were calculated from the Equations (38) and (39) (Table 3);
(4)
the volume temperature values ϑ 1 ( 0 ) = 471.97 of the disc and ϑ 2 ( 0 ) = 260.92 the pad were obtained from the Equations (36) and (37);
(5)
the values of materials properties K l , m ( ϑ l ( 0 ) ) , c l , m ( ϑ l ( 0 ) ) , ρ l , m ( ϑ l ( 0 ) ) , l , m = 1 , 2 , corresponding to the volume temperature ϑ l ( 0 ) were determined from the Formulas (40)–(51);
(6)
the steps (3)–(5) were repeated resulting in the corrected values for the volume temperature ϑ l ( 1 ) = 624.93 , and ϑ 2 ( 1 ) = 292.98 ;
(7)
by means of the formula ϑ l = 0.5 ( ϑ l ( 0 ) + ϑ l ( 1 ) ) , l = 1 , 2 final values of the volume temperature ϑ 1 = 548.45 C , and ϑ 2 = 267.95 C were found;
(8)
based on the dependencies (40)–(51) the values of materials properties K l , m ( ϑ l ) , c l , m ( ϑ l ) , ρ l , m ( ϑ l ) , l , m = 1 , 2 corresponding to the volume temperature ϑ l were established (Table 4) and other parameters necessary to perform the calculations (Table 5);
(9)
the temperature field Θ * ( ζ , τ ) (23)–(26), the temperature evolution Θ * ( τ ) (27), and temporal profiles of heat fluxes intensities q l * ( τ ) , l = 1 , 2 (30)–(32) were determined.
In order to calculate the values of Bessel functions J k ( x ) , k = 0 , 1 the programs BESSJ0 and BESSJ1 from the Numerical Recipes package [34] were used. The roots of the characteristic Equation (18) were searched for by the bisection method with the RTBIS program from this package. In summation of the series in solutions (23), (24), and (30), (31) was performed with an accuracy of 5 × 10 5 . For this accuracy, the minimum number of components was equal to 70.
Changes in the dimensionless temperature rise Θ * ( ζ , τ ) during braking, at few selected distances from the contact surface are presented in Figure 4. The temperature calculated with an account of the thermal sensitivity of the materials (solid lines) is significantly lower in both friction elements compared to the results achieved without taking into account the temperature dependencies of FGMs properties (dashed lines). The maximum dimensionless temperature on the contact surface ζ = 0 without and taking into account the thermal sensitivity of the materials are 0.816 and 0.277, respectively (reduction of about 2.94 times) and are reached at the time moments τ max = 0.37 and τ max = 0.29 (reduction of 21.6%).
Increasing the distance from the contact surface ζ = 0 , the temperature level of both elements drops (Figure 5). The temperature of components made of thermally sensitive materials is lower than their temperature, found for the constant material properties. The greatest difference between these results is on the contact surface.
The conclusions established on the basis of Figure 4 and Figure 5 confirm the results of the calculations, presented in Figure 6. It shows the dimensionless temperature isotherms Θ * ( ζ , τ ) . It can be seen that the effective depth of heat transfer is much greater in the case that material properties remain unchanged under the influence of temperature, than in the case of considering the thermally sensitive FGMs. This effect is most noticeable for the first one ( l = 1 ) , the Al2O3-Cu element. This result is also confirmed by the parameter values a l , l = 1 , 2 presented in Table 3 and Table 5.
The time profiles of the dimensionless intensities of heat fluxes q l * ( τ ) , l = 1 , 2 are shown in Figure 7. They decrease linearly during the braking process from the maximum value at the initial moment to zero at the stop. Most of the frictional heat generated is absorbed by the first element ( l = 1 ) Al2O3-Cu. The linear change in q l * ( τ ) is the result of the specific friction power q * ( τ ) (8), which decreases linearly during braking with a constant deceleration, and the requirement to meet the boundary condition q 1 * ( τ ) + q 2 * ( τ ) = q * ( τ ) , 0 τ τ s . The influence of thermal sensitivity on the intensity of heat fluxes is much smaller than on the temperature. For thermally sensitive materials, the maximum values of the intensity of heat fluxes are q 1 , max * = 0.864 and q 2 , max * = 0.136 , and for constant properties of the materials, we have q 1 , max * = 0.895 and q 2 , max * = 0.105 .

6. Conclusions

A calculation scheme was proposed to determine the temperature field of the friction elements of a disc brake, taking into account the changes in the FGMs properties depending on the actual temperature. The main part of the scheme was the adaptation of a linear solution (with temperature-independent material properties) to the thermal problem of friction during braking to thermally sensitive FGMs. A numerical analysis was performed in the case of braking with constant deceleration of elements made of two-component functionally graded materials with exponential variations in thermal conductivities in the axial direction, across the volume of the materials. It was found that:
  • the influence of thermal sensitivity on the temperature of FGMs may be more significant than in the case of homogeneous materials;
  • for the selected friction pair, taking into account the thermal sensitivity caused an almost threefold reduction in the maximum temperature in comparison to the appropriate temperature values, found with the same properties of the materials;
  • the influence of thermal sensitivity on the intensity of heat fluxes directed from the friction surface to the interior of the friction pair elements is insignificant. This means that to estimate the amount of heat absorbed by the individual elements of the friction pair, appropriate solutions to linear problems can be used.
A verification of the developed theoretical model based on empirical results would be advisable. However, no information on this kind of experimental data has been found in the literature. In particular, it concerns the frictional heating of braking systems with friction elements made of thermally sensitive FGMs. Therefore, the verification of the exact solution was obtained carried out by determining from it, in cases of limit parameters, known solutions of other authors for homogeneous materials, which were verified with appropriate experimental data. A new element, significantly differentiating the results of a given article from those published earlier by us, is the incorporation in the model of the possibility of changing the frictional properties of FGMs under temperature influence. This model includes many new elements, such as determining the intensity of heat fluxes to obtain the form of the heat partition ratio, finding the volume temperature of FGMs, developing a calculation algorithm that takes into account the thermal sensitivity of all materials components, etc. We have shown that taking into consideration the thermal sensitivity of materials can significantly reduce the surface temperature contact of the pad and disc. We proposed a theoretical computational model. We hope that it will be verified with the data obtained from other authors’ research positions. An indirect confirmation of the correctness of our model is also the time profiles of temperature and heat fluxes obtained on its basis, characteristic for braking with a constant deceleration.
It should be noted that all three of our papers constitute a monothematic cycle of interrelated research. We also want to develop a suitable model for braking systems operating in a short-term, repetitive mode. The problem of lowering the temperature level in such systems is up to date.

Author Contributions

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

Funding

This 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 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 , m Specific heat capacity ( J kg 1 K 1 )
f Coefficient of friction (dimensionless)
G l Weight of the friction elements (kg)
J k ( ) The Bessel functions of the first kind of the kth order
k l , m Thermal diffusivity ( m 2 s 1 )
K l , m Thermal conductivity ( W m 1 K 1 )
p Contact pressure ( Pa )
p 0 Nominal value of the contact pressure ( Pa )
R e External radius of the pads (m)
R i Internal radius of the pads (m)
q Specific power of friction ( W m 2 )
q0Nominal value of the specific power of friction ( W m 2 )
Q 0 Nominal friction power ( W )
t Time ( s )
t s Stop time ( s )
T Temperature ( C )
T 0 Initial temperature ( C )
V Velocity ( m s 1 )
V c , V ρ Volume fractions of the material phases
V 0 Initial velocity ( m s 1 )
W 0 Initial kinetic energy of the system ( J )
z Spatial coordinate in axial direction ( m )
lower l Number of the main ( l = 1 ) and frictional ( l = 2 ) elements of the friction pair
lower m Number of the component material m = 1 , 2 of selected friction element
α l Heat partition ratio (dimensionless)
β Cover angle of the pads (rad)
γ 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 , m Density ( kg m 3 )
τ Time (dimensionless)
τ s Time of braking (dimensionless)
ζ Spatial coordinate in axial direction (dimensionless)
ϑ l Volume temperature ( C )

References

  1. Balakin, V.; Sergienko, V. Thermal Calculation of Brakes and Frictional Assemblies; Belyi Metal-Polymer Research Institute of the National Academy of Sciences of Belarus: Gomel, Belarus, 1999. (In Russian) [Google Scholar]
  2. 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]
  3. Adamowicz, A.; Grześ, P. Three-dimensional FE model for calculation of temperature of a thermosensitive disc. Appl. Therm. Eng. 2013, 50, 572–581. [Google Scholar] [CrossRef]
  4. Yevtushenko, A.A.; Grześ, P.; Adamowicz, A. The temperature mode of the carbon-carbon multi-disc brake in the view of the interrelations of its operating characteristics. Materials 2020, 13, 1878. [Google Scholar] [CrossRef] [PubMed]
  5. Kushnir, R.M.; Popovych, V.S. Heat Conduction Problems of Thermosensitive Solids under Complex Heat Exchange; Vyacheslav, V., Ed.; Heat Conduction-Basic Research: Rijeka, Croatia, 2011; ISBN 978-953-307-404-7. [Google Scholar]
  6. Yevtushenko, A.A.; Kuciej, M.; Och, E. Some Methods for Calculating Temperature during the Friction of Thermosensitive Materials. Numer. Heat Transf. Part A Appl. 2015, 67, 696–718. [Google Scholar] [CrossRef]
  7. Yevtushenko, A.; Kuciej, M.; Och, E.; Yevtushenko, O. Effect of the thermal sensitivity in modeling of the frictional heating during braking. Adv. Mech. Eng. 2016, 8, 1–10. [Google Scholar] [CrossRef] [Green Version]
  8. Popovich, V.S.; Garmatii, G.Y. Solution of nonstationary heat conduction problems for thermosensitive bodies under convective heat exchange. J. Math. Sci. 1998, 90, 2037–2041. [Google Scholar] [CrossRef]
  9. 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]
  10. Zelentsov, V.B.; Lapina, P.A.; Mitrin, B.I. Wear of Functionally Graded Coatings under Frictional Heating Conditions. Nanomaterials 2022, 12, 142. [Google Scholar] [CrossRef]
  11. 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]
  12. Abbas, I.A. Nonlinear transient thermal stress analysis of thick-walled FGM cylinder with temperature-dependent material properties. Meccanica 2014, 49, 1697–1708. [Google Scholar] [CrossRef]
  13. Azadi, M.; Azadi, M. Nonlinear transient heat transfer and thermoelastic analysis of thick-walled FGM cylinder with temperature-dependent material properties using Hermitian transfinite element. J. Mech. Sci. Technol. 2009, 23, 2635. [Google Scholar] [CrossRef]
  14. Tanigawa, Y.; Akai, T.; Kawamura, R.; Oka, N. Transient heat conduction and thermal stress problems of a nonhomogeneous plate with temperature-dependent material properties. J. Therm. Stresses 1996, 19, 77–102. [Google Scholar] [CrossRef]
  15. Ding, S.; Wu, C.-P. Optimization of material composition to minimize the thermal stresses induced in FGM plates with temperature-dependent material properties. Int. J. Mech. Mater. Des. 2018, 14, 527–549. [Google Scholar] [CrossRef]
  16. Noda, N. Thermal stresses in functionally graded materials. J. Therm. Stresses 1999, 22, 477–512. [Google Scholar] [CrossRef]
  17. Ueda, S. Thermoelastic analysis of W-Cu functionally graded materials subjected to a thermal shock using a micromechanical model. J. Therm. Stresses 2001, 24, 19–46. [Google Scholar] [CrossRef]
  18. 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]
  19. Chichinadze, A.V. Polymers in Friction Assembles of Machines and Devices: A Handbook; Allerton Press Inc.: New York, NY, USA, 1984. [Google Scholar]
  20. Yevtushenko, A.; Kuciej, M.; Topczewska, K. Frictional Heating during Braking of the C/C Composite Disc. Materials 2020, 13, 2691. [Google Scholar] [CrossRef]
  21. 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]
  22. Watson, G.N. A Treatise of the Theory of Bessel Functions, 2nd ed.; Cambridge University Press: Cambridge, UK, 1995. [Google Scholar]
  23. Yevtushenko, A.A.; Kuciej, M. Two heat conduction problems with frictional heating during braking. J. Theor. Appl. Mech. 2010, 48, 367–380. [Google Scholar]
  24. Fazekas, G.A.G. Temperature gradients and heat stresses in brake drums. SAE Trans. 1953, 61, 279–284. [Google Scholar]
  25. 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]
  26. Burghartz, S.; Schulz, B. Thermophysical properties of sapphire, AlN and MgAl2O4 down to 70 K. J. Nucl. Mater. 1994, 212–215 Pt B, 1065–1068. [Google Scholar] [CrossRef]
  27. Kingery, W.D.; Francl, J.; Coble, R.L.; Vasilos, T. Thermal Conductivity: X, Data for Several Pure Oxide Materials Corrected to Zero Porosity. J. Am. Ceram. Soc. 1954, 37, 107–110. [Google Scholar] [CrossRef]
  28. Archer, D.A. Thermodynamic Properties of Synthetic Sapphire (α-Al2O3), Standard Reference Material 720 and the Effect of Temperature-Scale Differences on Thermodynamic Properties. J. Phys. Chem. Ref. Data 1993, 22, 1441–1453. [Google Scholar] [CrossRef] [Green Version]
  29. Simon, N.J.; Drexler, E.S.; Reed, R.P. Properties of Copper and Copper Alloys at Cryogenic Temperatures; Monograph 177; National Institute of Standards and Technology: New York, NY, USA, 1992.
  30. Pankratz, L.B. Thermodynamic Properties of Elements and Oxides; Bulletin 672; U.S. Bureau of Mines: Washington, DC, USA, 1982. [Google Scholar]
  31. Taylor, D. Thermal expansion data. II: Binary oxides with the fluorite and rutile structures, MO2, and the antifluorite structure, M2O. Trans. J. Br. Ceram. Soc. 1984, 83, 32–37. [Google Scholar]
  32. Deem, H.W.; Wood, W.D.; Lucks, C.F. The relationship between Electrical and Thermal Conductivities of Titanium Alloys. Trans. Metall. Soc. AIME 1958, 212, 520–523. [Google Scholar]
  33. Cezairrliyan, A.; McClure, J.L.; Taylor, R.J. Thermophysical Measurements on 90Ti-6Al-4V Alloy Above 1450 K Using a Transient (Subsecond) Technique. J. Res. Nat. Bur. Stand.-A Phys. Chem. 1977, 81, 251–256. [Google Scholar] [CrossRef]
  34. Press, W.H.; Teukolsky, S.A.; Flannery, B.P.; Vetterling, W.T. Numerical Recipes in Fortran 77: Volume 1, Volume 1 of Fortran Numerical Recipes: The Art of Scientific Computing; Cambridge University Press: Cambridge, UK, 1992. [Google Scholar]
Figure 1. Dependencies of the dimensionless thermal conductivities K l , m * on temperature T .
Figure 1. Dependencies of the dimensionless thermal conductivities K l , m * on temperature T .
Materials 15 00963 g001
Figure 2. Dependencies of the dimensionless specific heat capacities c l , m * on temperature T .
Figure 2. Dependencies of the dimensionless specific heat capacities c l , m * on temperature T .
Materials 15 00963 g002
Figure 3. Dependencies of the dimensionless densities ρ l , m * on temperature T .
Figure 3. Dependencies of the dimensionless densities ρ l , m * on temperature T .
Materials 15 00963 g003
Figure 4. Evolutions of the dimensionless temperature Θ * ( ζ , τ ) during braking at different distances ζ from the surface of friction with (solid lines) and without (dashed lines) taking into account the thermal sensitivity of the materials: (a) Al2O3—Cu; (b) ZrO2—Ti-6Al-4V.
Figure 4. Evolutions of the dimensionless temperature Θ * ( ζ , τ ) during braking at different distances ζ from the surface of friction with (solid lines) and without (dashed lines) taking into account the thermal sensitivity of the materials: (a) Al2O3—Cu; (b) ZrO2—Ti-6Al-4V.
Materials 15 00963 g004
Figure 5. Distribution of the dimensionless temperature Θ max * ( ζ ) = Θ * ( ζ , τ max ) reached at the time moment τ = τ max along the distance ζ from the surface of friction with (solid lines) and without (dashed lines) taking into account the thermal sensitivity of the materials.
Figure 5. Distribution of the dimensionless temperature Θ max * ( ζ ) = Θ * ( ζ , τ max ) reached at the time moment τ = τ max along the distance ζ from the surface of friction with (solid lines) and without (dashed lines) taking into account the thermal sensitivity of the materials.
Materials 15 00963 g005
Figure 6. Isotherms of the dimensionless temperature Θ * ( ζ , τ ) for: (a) thermally sensitivity materials; (b) materials with properties at the initial temperature.
Figure 6. Isotherms of the dimensionless temperature Θ * ( ζ , τ ) for: (a) thermally sensitivity materials; (b) materials with properties at the initial temperature.
Materials 15 00963 g006
Figure 7. Temporal profiles of the dimensionless heat fluxes q l * ( τ ) , l = 1 , 2 during braking with (continuous lines) and without (dashed lines) taking into account the thermal sensitivity of the materials. Dotted lines represent the dimensionless specific power of friction q * .
Figure 7. Temporal profiles of the dimensionless heat fluxes q l * ( τ ) , l = 1 , 2 during braking with (continuous lines) and without (dashed lines) taking into account the thermal sensitivity of the materials. Dotted lines represent the dimensionless specific power of friction q * .
Materials 15 00963 g007
Table 1. Input parameters.
Table 1. Input parameters.
Friction
Coefficient f
Nominal Pressure p 0 , MPa Initial Sliding Speed V 0 , ms 1 Initial Kinetic Energy W 0 , kJ Outer
Radius R e , mm
Inner Radius R i , mm Initial
Temperature T 0 , ° C
0.270.60223.8103.5437.526.520
Table 2. Material properties at the initial temperature T 0 .
Table 2. Material properties at the initial temperature T 0 .
Element IndexMaterial IndexMaterialThermal Conductivity
K l , m ( 0 ) , W m 1 K 1
Specific Heat Capacity
c l , m ( 0 ) , J k g 1 K 1
Density ρ l , m ( 0 ) , kg m 3
l = 1m = 1Al2O337.24727.293990.92
m = 2Cu402.65147.358947.92
l = 2m = 1ZrO21.94452.836102.16
m = 2Ti-6Al-4V6.87538.084431.79
Table 3. Calculated parameters at the initial temperature T 0 .
Table 3. Calculated parameters at the initial temperature T 0 .
Element Indexl = 1l = 2
c l ( 0 ) , J kg 1 K 1 437.3495.5
ρ l ( 0 ) , kgm 3 6469.45267
k l ( 0 ) × 10 6 , m 2 s 1 13.20.743
γ l * 2.3811.266
a l , mm 21.8545.193
G l , kg 0.31270.0605
α l 0.8960.104
Table 4. Material properties at volume temperature ϑ l , l = 1 , 2 .
Table 4. Material properties at volume temperature ϑ l , l = 1 , 2 .
Element IndexMaterial IndexMaterialThermal Conductivity
K l , m ( ϑ l ) , Wm 1 K 1
Specific Heat Capacity
c l , m ( ϑ l ) , J kg 1 K 1
Density ρ l , m ( ϑ l ) , kgm 3
l = 1m = 1Al2O310.191097.933945.59
m = 2Cu367.15401.898690.20
l = 2m = 1ZrO21.99552.676069.84
m = 2Ti-6Al-4V9.57615.444399.06
Table 5. Calculated parameters at volume temperature ϑ l , l = 1 , 2 .
Table 5. Calculated parameters at volume temperature ϑ l , l = 1 , 2 .
Element Indexl = 1l = 2
c l ( ϑ l ) , J kg 1 K 1 749.7584.9
ρ l ( ϑ l ) , kgm 3 6317.25233.8
k l ( ϑ l ) × 10 6 , m 2 s 1 2.150.65
γ l * 3.5851.583
a l , mm 8.8344.854
G l , kg 0.12340.0562
α l 0.8630.137
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. Influence of Thermal Sensitivity of Functionally Graded Materials on Temperature during Braking. Materials 2022, 15, 963. https://0-doi-org.brum.beds.ac.uk/10.3390/ma15030963

AMA Style

Yevtushenko A, Topczewska K, Zamojski P. Influence of Thermal Sensitivity of Functionally Graded Materials on Temperature during Braking. Materials. 2022; 15(3):963. https://0-doi-org.brum.beds.ac.uk/10.3390/ma15030963

Chicago/Turabian Style

Yevtushenko, Aleksander, Katarzyna Topczewska, and Przemysław Zamojski. 2022. "Influence of Thermal Sensitivity of Functionally Graded Materials on Temperature during Braking" Materials 15, no. 3: 963. https://0-doi-org.brum.beds.ac.uk/10.3390/ma15030963

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