Table of Contents
ISRN Mechanical Engineering
Volume 2011 (2011), Article ID 801732, 6 pages
Research Article

Prediction Model of the Fusion Microzone on Sliding Contact Surface

Institute of Tribology, Wuhan University of Technology, Wuhan 430070, China

Received 12 August 2011; Accepted 5 September 2011

Academic Editor: T. Basak

Copyright © 2011 Yan Lu and Zuomin Liu. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


The current paper is motivated by the need to understand the factors in generating the fusion microzone in sliding systems. The objectives are to analyze the different elements' varied influence on the engineering surface's temperature rise. The current paper developed the prediction model based on the thermal conduct theory. A solution based on the Green's function method is combined with the grid method for calculating the temperature rise and distribution. The research indicates that: frictional heat is closely related to the sliding velocity, its value is in proportion to the sliding velocity; the thermal properties of the material are one of the key points to decide the temperature rise; the load is another main factor in increasing the temperature rise; comparing with other elements, the roughness may be the least effective to the temperature rise.

1. Introduction

In the sliding contact process, a temperature variation in the contacting bodies that resulted from frictional heating will cause many related engineering problems, such as wear, fusion, and deformation. Since Blok [1] and Jaeger and Carslaw [2] pioneered the studies of temperature rise at the contact surfaces due to moving heat surfaces, there are more researches aimed at studying the temperature rise and distribution. Vick and Furey [3] deduced the base theory of the frictional heat and calculated a set of regularly arranged contact areas’ temperature distribution; Chen et al. [4] used the classic theory to calculate the local temperature rise that greatly influences the glue, wear, and local plastic deformation. Chen and Li [5] established a wear model to simulate sliding processes in order to analyze the frictional heating. Bogdanovich and Tkachuk [6] presented some experiments to understand thermal phenomena in a sliding contact, revealing that asperities in contact were found to undergo the pulse effect of the temperature approaching the melting point of one of the friction members. But few researches are focused on the temperature distribution of the realistic engineering surface which was constructed by a set of patches whose sizes and shapes are irregularly. The major objective of the present paper is to study the real contact area’s frictional heat. The developed sliding thermal model can deduce the equation of thermal conduction. The solution is used to demonstrate the effects of the contact area’s roughness, the parameters of materials, velocity of the sliding, and the ability of load.

2. Mathematic Model of Friction Heat

When the true engineering surfaces slide over each other, the frictional heat is generated by the contact area which is a set of discontinuous asperities with different sizes, varied shapes. Then the system of heat generation on frictional surfaces can be modeled as finite heat sources moving on half-infinite homogeneous body [79], as shown in Figure 1. It consists of two contacting semi-infinite regions, the rough surface with normal load 𝐹 and the flat one move with the sliding velocity 𝑉 in the 𝑦-direction and the frictional thermal will disappear to these two bodies (without considering the consumption of energy).

Figure 1: Theory figure of the rough surface’s sliding.

It can be seen clearly from the Figure 1 that the frictional heat is generated by part of the hole. It is the main factor in generating the fusion point. As the test shows that when M50 (Ra = 0.24–0.38 μm) sliding under the load 𝐹=10 N, test temperature 𝑇=400°C, with 𝑉=1 m/s, the fusion points are generated (Figure 2).

Figure 2: Fusion micro-zone on the M50 engineering surface.
2.1. Thermal Equation of the Contact Region

According to the theory of the thermal conduction:heat𝑡2heat𝑡1=theboundaryheatinow+heatgeneratedbytheheatsource.(1) Set 𝑢(𝑥,𝑦,𝑧,𝑡) to represent the temperature in which the body moves to the position (𝑥,𝑦,𝑧) at time 𝑡, 𝐶 is the specific heat capacity, 𝑉 is the relative velocity, 𝜌 is density, and 𝑘 is the coefficient of heat conducting, then𝐷𝑢𝐶𝜌𝑥,𝑦,𝑧,𝑡2𝑢𝑥,𝑦,𝑧,𝑡1=𝑑𝑥𝑑𝑦𝑑𝑧𝑡2𝑡1𝑑𝑡𝜕𝐷𝑘𝜕𝑢𝜕𝑛𝑑𝑆𝑡2𝑡1𝑑𝑡𝐷𝐶𝜌𝜕𝑢+𝜕𝑡𝑑𝑥𝑑𝑦𝑑𝑧𝑡2𝑡1𝑑𝑡𝐷𝐶𝜌𝜕𝑢𝜕𝑦𝜕𝑦=𝜕𝑡𝑑𝑥𝑑𝑦𝑑𝑧𝑡2𝑡1𝑑𝑡𝐷(𝑘𝑢)𝐶𝜌𝜕𝑢𝜕𝑡+𝐶𝜌𝑉𝜕𝑢𝜕𝑦=(𝑘𝑢)𝜕𝑢𝜕𝑡+𝑉𝜕𝑢𝜕𝑦=𝑎2𝜕2𝑢𝜕𝑥2+𝜕2𝑢𝜕𝑦2+𝜕2𝑢𝜕𝑧2𝑢𝑡=𝑎2Δ𝑢𝑉𝑢𝑦.(2) Here, Δ=𝜕2/𝜕𝑥2+𝜕2/𝜕𝑦2+𝜕2/𝜕𝑧2, 𝑎2=𝑘/𝐶𝜌 is heat diffusion coefficient.

Here the thermal boundary is the second boundary condition: 𝑘𝜕𝑢=𝜕𝑛𝑞,contactregion(𝑧=0)0,uncontactregion(𝑧0).(3) Here 𝑞 is the heat density generated by the sliding contact region’s frictional heat, it can be demonstrated by𝐹𝑞=𝜇sum𝑖=1𝐴𝑖𝑉.(4) Here 𝜇 is friction coefficient; 𝐹 is pressure of union area; 𝑉 is sliding velocity; sum𝑖=1𝐴𝑖 is efficient contact area (𝑖 is the number of the contact area). Assume that there is no energy loss during the sliding process. Initial condition at the beginning of the sliding 𝑡=0,𝑢(𝑥,𝑦,𝑧,𝑡)=𝑇0.(5)

2.2. Solving the Thermal Control Equation

The frictional heat control equation of the sliding contact region is constructed by (2), (3), and (5), then according to Green’s method of the partial differential equation (3), 𝑑𝑢=𝑞𝐴𝑖[𝐺]𝑧=0𝑑𝑥𝑑𝑦, the Green Function can be presented as1𝐺=4𝑎𝜋(𝑡𝑡)3exp𝑥𝑥+𝑦𝑦𝑉𝑡𝑡24𝑎2(𝑡𝑡).(6) According to it, the influence coefficient of the temperature rise can be solved: 𝑇=𝑡𝑡=0𝐴𝑖[𝐺]𝑧=0𝑑𝑥𝑑𝑦𝑑𝑡.(7) The contact patch is defined as the rectangular region, 𝑥-direction from 𝑖1/2 to 𝑖+1/2, 𝑦-direction from 𝑗1/2 to 𝑗+1/2, then the thermal influent coefficient can be presented as1𝑇=4𝜋𝑡𝑡=0𝑇𝑥𝑇𝑦𝑑𝑡,𝑇(8)𝑥=𝑥erf𝑎𝑥𝑖1/24𝑎2(𝑡𝑡)𝑥erf𝑎𝑥𝑖+1/24𝑎2(𝑡𝑡),𝑇𝑦=𝑦erf𝑏𝑦𝑗1/2𝑉𝑡𝑡4𝑎2(𝑡𝑡)𝑦erf𝑏𝑦𝑗+1/2𝑉𝑡𝑡4𝑎2(𝑡𝑡).(9)

The frictional heat in each node is influenced by the surrounding friction force over the entire nominal area. To calculate the temperature in just one node, the heat influence function of all the nodes must be integrated. To overcome this difficulty, the contact region Ω has been divided into 𝑀×𝑁 elements of the same size (Figure 3), so when the time of the sliding is 𝑡, the frictional heat matrix approach is often used:[𝑇]{𝑈}={𝑄}or𝑢=𝑀𝑁𝑚=0𝑛=0𝑞𝑖,𝑗,𝑡×𝑇𝑖,𝑗,𝑡.(10){𝑈} and {𝑄}, respectively, are the frictional heat, heat flux vectors with 𝑁 elements and [𝑇] is the coefficient of frictional heat matrix (𝑀×𝑁) which is only a function of the grid geometry. [𝑇] matrix needs to be calculated first during the contact simulation. After that, the frictional heat {𝑈} can be obtained simply by multiplying the heat flux {𝑄} with the matrix [𝑇]. The frictional heat matrix method is very effective in calculating the temperature rise.

Figure 3: Grid geometry of rough surface’s contact area.

The element 𝑇𝑖,𝑗,𝑡 of the frictional heat matrix [𝑇] could be physically interpreted as the influence factor of one unit frictional heat at node 𝑗 to the frictional heat at node 𝑖. The validity of this statement depends on the types of the interpolating functions used for the frictional heat but the overall concept is still valid. According to the laws of thermodynamics, this influence factor is only a function of the distance between the two nodes. This means that the elements of [𝑇] are determined only by the grid geometry.

The frictional heat coefficient of each element can be calculated by (9), then substitute its result to (8), and the temperature rise coefficient which was generated by the moving heat source can be worked out. Next couple it to the heat flux of every node, and finally, the temperature rise and distribution of the sliding rough surface can be obtained by (10).

2.3. Algorithm and Program

Engineering rough surface can be like a fractal phenomenon, appearing random, chaotic, self-similarity, self-affinity, and multiscale nature [8]. The fractal surface generated by square grid can be used to simulate the engineer surface efficiently. The main idea of this algorithm is that first to get a few points of the surface, then to increase their height randomly according to the normal distribution 𝑁(0,1), and then incremental variance satisfies the equation 𝑉(𝑡)=𝐸[𝑥(𝑡)𝑥(0)]2=|𝑡|2𝐻𝜎20, and the height of midpoint of time range is given by interpolation; repeat this process 𝑛 times continuously until the sample variance decreases to 𝜎2𝑛=(1/2)2𝐻𝜎2𝐻𝑛𝑛1𝜎20. According to it the coordinate value of the rough surface’s nodes can be obtained. And the contact area has been solved by the previous research in [10]; then the thermal influent coefficient of total nodes can be worked out; the frictional heat coefficient matrix [𝑇] can be solved by this results easily. At the same time, the total contact area of whole patches can be obtained as well then to calculate the heat flux {𝑄}. The overall flow chart for the calculation is given in Figure 4.

Figure 4: Overall flow chart for solving true rough surface’s temperature rise.

3. Result and Discussion

High-speed steel M50 was chosen as the 3D fractal profile to slide against the ceramic material Si3N4 which was the semi-infinite elastic medium. Their thermal physical parameters are listed in Table 1 [11].

Table 1: Thermal performance of materials.

M50: it is a square with the length of side 0.1 μm. The true surface was established by using the fractal theory (Figure 5). Its roughness Ra = 6.3 μm, the total of the nodes is 10×10; the surface of material Si3N4 was seen as a relative smooth [12, 13].

Figure 5: Fractal model of engineering surface.

The first calculation condition is as follows: the load 𝐹=10 N, the sliding velocity was varied between 1 and 100 m/s. So the calculated temperature could be changed from room temperature to 1300°C. The relationship between the sliding velocity and the temperature rise has been plotted as two curves in Figure 6.

Figure 6: Effects of sliding speed on temperature rise.

Each curve in Figure 6 shows a significant temperature rise following the sliding velocity increase. This indicates that the sliding velocity is one of the most important factors which generate the frictional heat.

The temperature of M50 is higher than the Si3N4 under the same sliding velocity. And as the speed increases, the temperature rise of M50 increases faster than Si3N4. The result shows that the frictional heat of the material which thermal characteristic performances better will increase less. It indicates that heat diffusion coefficient, which decides the thermal property, may be one of the key points to temperature rise.

The second calculation conditions is as follows: the sliding velocity 𝑉=10 m/s, the load 𝐹=10 N, and the roughness, respectively, was 1.2 μm, 1.6 μm, 3.2 μm, and 6.3 μm. So the calculated temperature could be changed from 160 to 230°C. The relationship between the roughness and the temperature rise has been plotted as two curves in Figure 7. The result shown here reflects effects of roughness on temperature. The total contact area sum𝑖=1𝐴𝑖, which is related to the roughness of the surface, is in proportional to the frictional heat, which tends to increase the temperature.

Figure 7: Effects of roughness on temperature rise.

The third calculation condition is as follows: the sliding velocity 𝑉=10 m/s, the roughness Ra = 6.3 μm, and the load was varied between 10 and 100 N. So the calculation temperature could be changed from 200 to 1800°C. The relationship between the roughness and the temperature rise has been plotted as two curves in Figure 8. The result shown here reflects one major effect of load on temperature. The amount of frictional energy deposited on the surface is directly proportional to the load, which tends to increase temperature.

Figure 8: Effects of load on temperature rise.

Figures 9 and 10 show a comparison temperature rise contours between M50 (Ra = 6.3 μm, 𝐹=10 N, and 𝑉=10 m/s) and Si3N4 (relative smooth). It can be observed that frictional heat distribution is changed by the contact area, especially that below asperity contacts in the rough surface where thermoelastic deformation promoted the interaction of the thermal conduct of the neighboring asperity. These analytical results are supported by FEM simulations and it shows a strong effect of the contact area on the location of the frictional heating.

Figure 9: Temperature rise contours of M50.
Figure 10: Temperature rise contours of Si3N4.

Figure 11 shows a temperature rise contours of M50 on the same test condition of Figure 2 (Ra = 3.2 μm, 𝐹=10 N, 𝑇=400°C, and 𝑉=1 m/s). The result can totally interpret why there are some fusion points appearing on the M50 friction surface (the melting point of material M50 is 550°C). From this method, the fusion micro-zone can be decided. It helps us to predict the fusion situation of each material’s engineering surface in extreme friction condition.

Figure 11: Temperature rise contours of M50 on the same condition as Figure 2.

4. Conclusions

The focus of this paper is to deduce prediction model which can calculate the fusion micro-zone on the sliding contact surface from the viewpoint of the surface temperature rise and its distribution. It can be seen clearly that high local temperature decides the fusion zone on the engineering surface; the elements which contribute to the friction heat are the elements effect on the fusion zone. Case studies reveal five important trends.(1)The maximum frictional heat is closely related to the sliding velocity, and it increases as the sliding velocity increases.(2)The thermal properties of the material are one of the key points to temperature rise. The temperature rise decreases as the heat diffusion coefficient increases.(3)The shape of the temperature distribution looks like an eccentric ring. The maximum hot point lies in its center.(4)The roughness of the engineering surface has relative smaller influence on the temperature rise compared with the sliding velocity and thermal properties.(5)The load is another main factor in increasing the temperature rise.

The important conclusion from this study is that it predicts the true engineer surface’s fusion micro-zone, in the sliding process. The results may be useful for selecting material and designing the surface accuracy in the extreme working condition.


The authors would like to thank National Natural Science Foundation of China for financial support (ID 51075311).


  1. H. Blok, “Theoretical study of temperature rise at surfaces of actual contact under oiliness lubricating conditions,” in Proceedings of the General Discussion on Lubrication and Lubricants, vol. 2, pp. 222–235, Institution of Mechanical Engineers, London, UK.
  2. J. C. Jaeger and H. S. Carslaw, Conduction of Heat in Solids, Clarendon Press, Oxford, UK, 1959.
  3. B. Vick and M. J. Furey, “A basic theoretical study of the temperature rise in sliding contact with multiple contacts,” Tribology International, vol. 34, no. 12, pp. 823–829, 2001. View at Publisher · View at Google Scholar · View at Scopus
  4. H. Chen, W. Hu, H. Wang, and W. Wang, “Calculation of temperature fields of bodies in sliding con tact without lubrication,” Journal of Tsinghua University(Science and Technology), vol. 47, no. 11, pp. 1962–1964, 2007. View at Google Scholar
  5. Q. Chen and D. Y. Li, “A computational study of frictional heating and energy conversion during sliding processes,” Wear, vol. 259, no. 7–12, pp. 1382–1391, 2005. View at Publisher · View at Google Scholar · View at Scopus
  6. P. N. Bogdanovich and D. V. Tkachuk, “Temperature distribution over contact area and “hot spots” in rubbing solid contact,” Tribology International, vol. 39, no. 11, pp. 1355–1360, 2006. View at Publisher · View at Google Scholar · View at Scopus
  7. G. Liu, Q. Wang, and Y. Ao, “Convenient formulas for modeling three-dimensional thermo-mechanical asperity contacts,” Tribology International, vol. 35, no. 7, pp. 411–423, 2002. View at Publisher · View at Google Scholar · View at Scopus
  8. B. B. Mandelbrot, “Multifractal and fractal,” Physics Today, vol. 39, no. 9, pp. 11–12, 1986. View at Publisher · View at Google Scholar
  9. A. Majumdar and B. Bhushan, “Fractal model of elastic-plastic contact between rough surfaces,” ASME—Journal of Tribology, vol. 113, no. 11, pp. 1–11, 1991. View at Publisher · View at Google Scholar
  10. Y. Lu and Z. Liu, “Friction surface contact deformation model based onmaterial thermal parameters,” Chinese Journal of Mechanical Engineering, vol. 46, no. 9, pp. 28–33, 2010. View at Google Scholar
  11. Z. Liu, “Study of critical friction and wear temperature of M50 high speed steel,” Materials for Mechanical Engineering, vol. 21, no. 2, pp. 21–24, 1997. View at Google Scholar
  12. J. Yang and K. Komvopoulos, “A mechanics approach to static friction of elastic-plastic fractal surfaces,” ASME—Journal of Tribology, vol. 127, no. 2, pp. 315–324, 2005. View at Publisher · View at Google Scholar · View at Scopus
  13. H. Sofuoglu and A. Ozer, “Thermomechanical analysis of elastoplastic medium in sliding contact with fractal surface,” Tribology International, vol. 41, no. 8, pp. 783–796, 2008. View at Publisher · View at Google Scholar · View at Scopus