PredictionModel of the FusionMicrozone on Sliding Contact Surface

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.


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.

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 [7][8][9], as shown in Figure 1.It consists of two contacting semi-infinite regions, the rough surface with normal load F and the flat one move with the sliding velocity V in the y-direction and the frictional thermal will disappear to these two bodies (without considering the consumption of energy).
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 F = 10 N, test temperature T = 400 • C, with V = 1 m/s, the fusion points are generated (Figure 2).

Thermal Equation of the
(2) Here, Here the thermal boundary is the second boundary condition: Here q is the heat density generated by the sliding contact region's frictional heat, it can be demonstrated by Here μ is friction coefficient; F is pressure of union area; V is sliding velocity; sum i=1 A i is efficient contact area (i 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 t = 0, u x, y, z, t = T 0 . (5)

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), du = q Ai [G] z=0 dx dy , the Green Function can be presented as According to it, the influence coefficient of the temperature rise can be solved: The contact patch is defined as the rectangular region, xdirection from i − 1/2 to i + 1/2, y-direction from j − 1/2 to j +1/2, then the thermal influent coefficient can be presented as ISRN Mechanical Engineering 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 M × N elements of the same size (Figure 3), so when the time of the sliding is t, the frictional heat matrix approach is often used: {U} and {Q}, respectively, are the frictional heat, heat flux vectors with N elements and [T] is the coefficient of frictional heat matrix (M × N) which is only a function of the grid geometry.
[T] matrix needs to be calculated first during the contact simulation.After that, the frictional heat {U} can be obtained simply by multiplying the heat flux {Q} with the matrix [T].The frictional heat matrix method is very effective in calculating the temperature rise.The element T i, j,t of the frictional heat matrix [T] could be physically interpreted as the influence factor of one unit frictional heat at node j to the frictional heat at node i.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 [T] 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).

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 N(0, 1), and then incremental variance  satisfies the equation , and the height of midpoint of time range is given by interpolation; repeat this process n times continuously until the sample variance decreases to σ 2 n = (1/2) 2H σ 2Hn n−1 σ 2 0 .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 [T] 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 {Q}.The overall flow chart for the calculation is given in Figure 4.

Result and Discussion
High-speed steel M50 was chosen as the 3D fractal profile to slide against the ceramic material Si 3 N 4 which was the semiinfinite elastic medium.Their thermal physical parameters are listed in Table 1 [11].
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 Si 3 N 4 was seen as a relative smooth [12,13].
The first calculation condition is as follows: the load F = 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.
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 Si 3 N 4 under the same sliding velocity.And as the speed increases, the temperature rise of M50 increases faster than Si 3 N 4 .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 V = 10 m/s, the load F = 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 i=1 A i , which is related to the roughness of the surface, is in proportional to the frictional heat, which tends to increase the temperature.
The third calculation condition is as follows: the sliding velocity V = 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.
Figures 9 and 10 show a comparison temperature rise contours between M50 (Ra = 6.3 μm, F = 10 N, and V = 10 m/s) and Si 3 N 4 (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 11 shows a temperature rise contours of M50 on the same test condition of Figure 2 (Ra = 3.2 μm, F = 10 N, T = 400 • C, and V = 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.

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.
Input the measure data of the rough surface Using the Brown's fractal function, to simulate the whole Calculate the contact area and the pressure of each node Solve the thermal influent coefficient matrix {T} Solve the heat flux {Q} Solve the temperature rise of the rough surface {U} = {Q}{T}

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

Figure 7 :
Figure 7: Effects of roughness on temperature rise.

Figure 8 :Figure 9 :
Figure 8: Effects of load on temperature rise.

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

Table 1 :
Thermal performance of materials.