# Contact Thermal Analysis and Wear Simulation of a Brake Block

Authors: Nándor Békési; Károly Váradi
Journal: Advances in Tribology
 (2013)

---

## Abstract

The present paper describes an experimental test and a coupled contact-thermal-wear analysis of a railway wheel/brake block system through the braking process. During the test, the friction, the generated heat, and the wear were evaluated. It was found that the contact between the brake block and the wheel occurs in relatively small and slowly moving hot spots, caused by the wear and the thermal effects. A coupled simulation method was developed including numerical frictional contact, transient thermal and incremental wear calculations. In the 3D simulation, the effects of the friction, the thermal expansion, the wear, and the temperature-dependent material properties were also considered. A good agreement was found between the results of the test and the calculations, both for the thermal and wear results. The proposed method is suitable for modelling the slowly oscillating wear caused by the thermal expansions in the contact area.


---

## Body

## 1. Introduction

Applications with dry sliding friction between metal parts are common in classical mechanical engineering practice [1, 2]. However, the processes in railway brake systems are quite complex. There are many studies about the steady-state friction of such applications, some of them considering even the wear and the thermal effects [3–6]; yet the real-life braking is typically a transient phenomenon, the vehicle is usually stopped in seconds or in a few minutes. Since most of the kinetic energy of the train is transferred to heat by the friction between the wheels and the brake blocks, a big amount of heat is developing in a short time. It was found earlier that in the beginning of the braking, the edges of the brake blocks contact the wheel, even when the brake block is worn. This can be explained by the thermal expansion of the brake blocks: at the end of the previous braking, the brake block is worn, so the radii of the wheel and the brake block approximately are equal. Then, the system cools down, and since there is a great decrease in temperature, a significant thermal deformation occurs, which causes the brake block to have a shape that fits the wheel only at the edges. There is another interesting feature that is called thermoelastic instability (TEI) in the literature, which is responsible for the appearance and moving of hot spots in the apparent contact area. The TEI was investigated in a number of papers (e.g., in [7–11]), and it was found to be caused by the joint effect of wear and thermal expansion of the high-temperature hot spots.present paper shows experimental and numerical analysis of the thermal and tribological processes of the wheel/brake block system during braking, particularly in the running-in phase.


## 2. Experimental

In order to cast light on the processes during braking, a wear test was performed at the Department of Railway Vehicles and Vehicle System Analysis (Budapest University of Technology and Economics, Budapest, Hungary) on a 1 : 4 scaled test rig. In the experimental setup, a brake block was pressed against a rotating wheel (Figure1).Figure 1
The test rig.


### 2.1. The Wear Test Parameters

Within the frame of this research, our main goal was to understand the phenomena occurring in the running-in phase of the braking process; therefore, the total time of the test was 50 s. The normal and the frictional force was measured by load cells and registered during the test.The brake block was 92 mm high, 20 mm wide, and its initial radius was 119 mm (Figure2). The wheel is 35 mm wide; its radius is 121 mm. The brake block was pressed against the wheel with a normal force of 500 N by a hydraulic cylinder. The wheel was rotating at 875 rpm, so the sliding velocity was 40 km/h. The rotational speed of the wheel and the normal force were kept constant during the test.Figure 2
The main dimensions of the brake block.The heating of the brake block was observed by three thermocouples. The thermocouples were fitted in holes 5 mm under the surface in the brake block at the positions indicated in Figure2. The temperature data were registered during the test online. The temperature of the wheel was measured three times (at 0 s, 25 s, and 50 s) by a laser thermometer (Advanced ThermoLaser TMTL 1400 K, SKF, Sweden). The thermometer was calibrated for the front face of the stationary wheel before the test, then the temperature was measured on the same spot of the rotating wheel.


### 2.2. Results of the Wear Test

Figure3 shows the worn surface of the brake block after the test. It can be seen that the wear of the brake block is limited to the two contact areas at the two edges and furthermore that the wear is not symmetric: the wear on the left side (at the leading edge) is slightly higher than that on the right side (at the trailing edge). The asymmetry of the wear pattern was caused by the friction: the contact pressure distribution was asymmetric due to the friction and so was the wear intensity of the brake block. This phenomenon is well known and is responsible for the pivoting of the brake block on the pin.Figure 3
The worn surface of the brake block after 50 s of sliding (the arrow indicates the sliding direction). The dark areas are worn.The measured coefficient of friction was 0.38.The temperature of the measured points versus time is shown in Figure4. It can be seen that the two edges of the brake block heat up quickly since these are close to the contact zones, while the measured point in the middle is heated less intensively because there is no contact between the central part of the brake block and the wheel. The wavy feature of the upper curve can be attributed to the motion of the hotspots over the apparent contact area. The real contact occurs only in small spots in the apparent contact area, when these hot spots move due to wear and thermal expansion. The change in their distance from the thermocouples affects the measured temperature. The curve of the lower point is less wavy because the generated heat is less in this side. The temperature curve of the middle point shows no waviness because it is relatively far from the contact spots, so the motion of the hot spots has no significant effect.Figure 4
temperature at the measured points of the brake block.Figure5 shows the temperature development of the wheel. It can be seen that the temperature curve is rather linear.Figure 5
The temperature of the wheel versus time.


## 2.1. The Wear Test Parameters

Within the frame of this research, our main goal was to understand the phenomena occurring in the running-in phase of the braking process; therefore, the total time of the test was 50 s. The normal and the frictional force was measured by load cells and registered during the test.The brake block was 92 mm high, 20 mm wide, and its initial radius was 119 mm (Figure2). The wheel is 35 mm wide; its radius is 121 mm. The brake block was pressed against the wheel with a normal force of 500 N by a hydraulic cylinder. The wheel was rotating at 875 rpm, so the sliding velocity was 40 km/h. The rotational speed of the wheel and the normal force were kept constant during the test.Figure 2
The main dimensions of the brake block.The heating of the brake block was observed by three thermocouples. The thermocouples were fitted in holes 5 mm under the surface in the brake block at the positions indicated in Figure2. The temperature data were registered during the test online. The temperature of the wheel was measured three times (at 0 s, 25 s, and 50 s) by a laser thermometer (Advanced ThermoLaser TMTL 1400 K, SKF, Sweden). The thermometer was calibrated for the front face of the stationary wheel before the test, then the temperature was measured on the same spot of the rotating wheel.


## 2.2. Results of the Wear Test

Figure3 shows the worn surface of the brake block after the test. It can be seen that the wear of the brake block is limited to the two contact areas at the two edges and furthermore that the wear is not symmetric: the wear on the left side (at the leading edge) is slightly higher than that on the right side (at the trailing edge). The asymmetry of the wear pattern was caused by the friction: the contact pressure distribution was asymmetric due to the friction and so was the wear intensity of the brake block. This phenomenon is well known and is responsible for the pivoting of the brake block on the pin.Figure 3
The worn surface of the brake block after 50 s of sliding (the arrow indicates the sliding direction). The dark areas are worn.The measured coefficient of friction was 0.38.The temperature of the measured points versus time is shown in Figure4. It can be seen that the two edges of the brake block heat up quickly since these are close to the contact zones, while the measured point in the middle is heated less intensively because there is no contact between the central part of the brake block and the wheel. The wavy feature of the upper curve can be attributed to the motion of the hotspots over the apparent contact area. The real contact occurs only in small spots in the apparent contact area, when these hot spots move due to wear and thermal expansion. The change in their distance from the thermocouples affects the measured temperature. The curve of the lower point is less wavy because the generated heat is less in this side. The temperature curve of the middle point shows no waviness because it is relatively far from the contact spots, so the motion of the hot spots has no significant effect.Figure 4
temperature at the measured points of the brake block.Figure5 shows the temperature development of the wheel. It can be seen that the temperature curve is rather linear.Figure 5
The temperature of the wheel versus time.


## 3. Numerical Modelling

In order to understand the joint effect of wear and the thermal expansions on the contact of the brake block and the wheel, mechanical and thermal numerical simulations were performed.The simulation process is shown in Figure6. One cycle of the simulation can be divided into four main blocks.Figure 6
The flowchart of the simulation process.First, a frictional contact FE calculation is run. In this phase, the brake block is pressed against the wheel, which is rotating at a constant rotational speed according to the measurement. Based on the contact pressure distribution determined by the first simulation block, a transient thermal FE calculation is performed next.The frictional heat increment is determined for each contacting node, and the frictional heat values are applied as point flux type boundary conditions in the thermal FE model of the brake pad.At the beginning of the simulation, room temperature was assumed at the entire model. In the later cycles, the initial temperature distribution was obtained from the result of the previous thermal simulation.In the next block, another FE contact calculation takes place, but now the effect of the thermal expansion is taken into account as well. In the last block of the cycle, we determine the nodal wear increments based on the contact pressure distribution evaluated by (1). Once the nodal wear increments are known, we move the contacting nodes according to them. The Δhi nodal wear of the ith node in the contact area is defined as follows [12]:
(1)Δhi=WspivΔt,
where Ws is the specific wear rate, pi is the nodal contact pressure at node i, v is the sliding velocity, and Δt is the time increment.In order to reduce the effects of the high peaks in the contact pressure distribution caused by the relatively high local temperatures in the contact zone, a smoothing procedure was applied. When a node is in contact, the neighbouring nodes were also moved according to the moving average contact pressure value of the patch around the node.The time increment of the simulation wasΔt=1s, so the total time of the experiment (50s) was modelled in 50 simulation cycles.


### 3.1. The FE Model

The 3D FE model was created using the MSC.Marc software package. The same FE mesh was used in all of the four blocks of the simulation process. The brake block was modelled by 3918 8-node 3D elements of different size (Figure7). Symmetry was considered in the model. In the contact calculations, the wheel was modelled as a rigid body, since the deformation of the wheel would not change the contact conditions significantly. The normal force was applied through a rigid shaft. A coefficient of friction with a value of 0.38 was prescribed between the brake block and the wheel. In the thermal calculations, a thermal face film boundary condition with a value of 20 W/(m2 K) was applied for the surfaces of the brake block except for the symmetry plane, the contact zones, and the hole. On the surface of the hole, another thermal face film boundary condition with a value of 100 W/(m2 K) was applied in order to approximately model the heat flow through the pin.Figure 7
The FE mesh of the brake block with the wheel and the pin considered as rigid bodies.


### 3.2. Frictional Heat Generation

The frictional heat generation was determined analytically as the following. The total heat power equals the work done by the friction force:(2)Q=Ffrv=Fnμv,
where Ffr is the friction force, Fn is the normal force, μ is the coefficient of friction, and v is the sliding velocity.The total heat is distributed between the two contacting bodies. It was found earlier in [13, 14] that the larger part of the heat is conducted by the wheel. The ratio of the heat partition is found to be changing over the time [15]. However, a ratio of heat partition (Khp) with constant value was assumed in the series of calculations. The total heat  is as follows:(3)Q=Qwheel+Qblock=(1-Khp)Q+KhpQ,
where Khp is the ratio of heat partition representing the block side.In the FE calculations, the following material properties were assumed to model the cast iron (P10) material of the brake block:(i)
Poisson’s ratio:υ=0.25,(ii)
density:ρ=7150 kg/m3,(iii)
Young’s modulus:E=-0.00002T2-0.0272T+103,(iv)
linear thermal expansion coefficient:α=0.55E-08T+0.00001014,(v)
specific heat:c=0.0008T2+0.0295T+521.06,(vi)
thermal conductivity:λ=0.00004T2+0.0025T+48.115,(vii)
ratio of heat partition:Khp=0.25.The values of these temperature dependent material properties for different temperatures are shown in Table1.Table 1
Values of the temperature-dependent material properties.
t  (°C)
E  (GPa)
α  (1/K)
c  (J/kgK)
λ  (W/mK)
0
103.00
1.01
E
-
05
521.06
48.12
25
102.31
1.03
E
-
05
522.30
48.03
100
100.08
1.07
E
-
05
532.01
47.47
200
96.76
1.12
E
-
05
558.96
46.02
300
93.04
1.18
E
-
05
601.91
43.77
400
88.92
1.23
E
-
05
660.86
40.72
500
84.40
1.29
E
-
05
735.81
36.87
600
79.48
1.34
E
-
05
826.76
32.22
700
74.16
1.40
E
-
05
933.71
26.77
800
68.44
1.45
E
-
05
1056.66
20.52
900
62.32
1.51
E
-
05
1195.61
13.47
1000
55.80
1.56
E
-
05
1350.56
5.62


## 3.1. The FE Model

The 3D FE model was created using the MSC.Marc software package. The same FE mesh was used in all of the four blocks of the simulation process. The brake block was modelled by 3918 8-node 3D elements of different size (Figure7). Symmetry was considered in the model. In the contact calculations, the wheel was modelled as a rigid body, since the deformation of the wheel would not change the contact conditions significantly. The normal force was applied through a rigid shaft. A coefficient of friction with a value of 0.38 was prescribed between the brake block and the wheel. In the thermal calculations, a thermal face film boundary condition with a value of 20 W/(m2 K) was applied for the surfaces of the brake block except for the symmetry plane, the contact zones, and the hole. On the surface of the hole, another thermal face film boundary condition with a value of 100 W/(m2 K) was applied in order to approximately model the heat flow through the pin.Figure 7
The FE mesh of the brake block with the wheel and the pin considered as rigid bodies.


## 3.2. Frictional Heat Generation

The frictional heat generation was determined analytically as the following. The total heat power equals the work done by the friction force:(2)Q=Ffrv=Fnμv,
where Ffr is the friction force, Fn is the normal force, μ is the coefficient of friction, and v is the sliding velocity.The total heat is distributed between the two contacting bodies. It was found earlier in [13, 14] that the larger part of the heat is conducted by the wheel. The ratio of the heat partition is found to be changing over the time [15]. However, a ratio of heat partition (Khp) with constant value was assumed in the series of calculations. The total heat  is as follows:(3)Q=Qwheel+Qblock=(1-Khp)Q+KhpQ,
where Khp is the ratio of heat partition representing the block side.In the FE calculations, the following material properties were assumed to model the cast iron (P10) material of the brake block:(i)
Poisson’s ratio:υ=0.25,(ii)
density:ρ=7150 kg/m3,(iii)
Young’s modulus:E=-0.00002T2-0.0272T+103,(iv)
linear thermal expansion coefficient:α=0.55E-08T+0.00001014,(v)
specific heat:c=0.0008T2+0.0295T+521.06,(vi)
thermal conductivity:λ=0.00004T2+0.0025T+48.115,(vii)
ratio of heat partition:Khp=0.25.The values of these temperature dependent material properties for different temperatures are shown in Table1.Table 1
Values of the temperature-dependent material properties.
t  (°C)
E  (GPa)
α  (1/K)
c  (J/kgK)
λ  (W/mK)
0
103.00
1.01
E
-
05
521.06
48.12
25
102.31
1.03
E
-
05
522.30
48.03
100
100.08
1.07
E
-
05
532.01
47.47
200
96.76
1.12
E
-
05
558.96
46.02
300
93.04
1.18
E
-
05
601.91
43.77
400
88.92
1.23
E
-
05
660.86
40.72
500
84.40
1.29
E
-
05
735.81
36.87
600
79.48
1.34
E
-
05
826.76
32.22
700
74.16
1.40
E
-
05
933.71
26.77
800
68.44
1.45
E
-
05
1056.66
20.52
900
62.32
1.51
E
-
05
1195.61
13.47
1000
55.80
1.56
E
-
05
1350.56
5.62


## 4. Results of the Simulation

Figure8 shows the calculated wear distribution of the brake block after 50 s of sliding. It can be seen that the wear is asymmetric as it was in the experimental results. The left side of the brake block suffered more wear than the right side. On the other hand, the wear intensity does not differ significantly along the axial direction, which again is in accordance with the experimental test results. The quantitative comparison of the test and the simulation results shows slight differences. The worn areas of the experiment are approximately 14 mm and 10 mm long, while, in the simulation, those are 12 mm and 9 mm, respectively.Figure 8
Calculated wear distribution after 50 s of sliding.The temperature of the nodes in the FE model that were in the position of the thermocouples of the test was obtained from the simulation result files. The calculated temperature curves of the three measured points of the brake block are shown in Figure9. The temperature of the upper point is the highest; the temperature of the lower point is slightly lower. It can be seen that the calculated and the measured curves are in good agreement; however, the calculated curves near the edges show a wavier feature than those of the measurement. The temperature values of the three measured points of the test and the calculation are in good agreement.Figure 9
The calculated and measured temperatures of the three measured points versus time.Figure10 shows the temperature distribution of the brake block at the end of the simulation. One can see that the parts near the edges of the contact area have the highest temperature. On the top, cylindrical face high-temperature spots are visible due to the spot-like contact of the last simulation cycle. The lowest temperature appears around the hole, where the heat flow through the pin was considered.Figure 10
Temperature (in °C) distribution of the brake block at the end of the simulation (50 s).In order to better understand the physical processes, the simulation was extended for an additional 40 s. Figure11 summarizes the contact stress distributions of the contact face at the beginning and after 10, 30, 50, 70, and 90 s of sliding. It can be seen that, in the beginning, the two edges of the brake block were in contact. Later, the contact takes place in spots, which move cycle by cycle as the geometry changes due to wear and thermal expansion.Figure 11
Contact stress distribution at the beginning and after 10, 30, 50, 70, and 90 s of sliding (averaged along the axial direction).Figure12 shows the temperature distribution on the sliding surface of the brake block during the simulation. One can compare the curves of Figures 11 and 12 and see that the locations near the contacting spots have a high peak in the temperature distribution, since the friction heat is proportional to the contact pressure. It can be observed that the contact zones are moving from the initial position on the edges to the middle at first but then moving to the edges again. This can be explained by the following: consider one peak where the brake block contacts the rotating wheel. The contact spot has a peak in the contact pressure distribution, so the evolving frictional heat will cause a relatively high temperature at that location. The thermal expansion increases the contact pressure of this spot because of the high temperature. This area will suffer the highest wear, and the load will be shifted to the areas next to this spot (see the 30 s and the 50 s curves in Figure 11). This motion of the contact zones and the temperature peaks causes the wavy feature of the temperature curves in the measured points because the distance between the heat sources and the measured points is changing.Figure 12
Temperature distribution on the surface at the beginning and after 10, 30, 50, 70 and 90 s of sliding (averaged along the axial direction).


## 5. Conclusions

The braking process of a railway brake block and wheel was investigated experimentally and numerically. The phenomenon known as the motion of hot spots was observed in the experimental test. A coupled numerical contact-thermal-wear simulation method was developed and applied in order to model the physical processes occurring during braking. In the cycles of the simulation, the frictional heat generation was determined based on the frictional contact calculation, and then the thermal expansions were calculated by the thermal calculation. In the last part of the simulation, the wear increments were determined, and the nodes were moved, and the new cycle was started with the new geometry. The same FE mesh was used in the different aspects of the simulation process. The wear and the temperature distribution of the calculation and the experimental test results were mostly similar.The nonuniform wear process caused by the thermal expansion of the contacting “hot spots” was modelled.



---

*Source: 878274-2013-02-27.xml*