# Friction coefficient and meshing efficiency of spiral bevel gears

The mixed lubrication model is composed of basic equations, Reynolds equation is composed of oil film equation, density equation, viscosity equation and load balance equation, which is solved by composite iterative method, and the local elastic deformation integral is solved by fast Fourier transform method. Elastic modulus of gear material e = 219.78 GPA, viscosity of lubricating oil η 0 = 0.09 Pa.s, viscosity pressure coefficient α= 12.5 GPa–1。 The surface roughness of spiral bevel gear is 0.2 μ m. Root mean square roughness σ= zero point two eight μ m。

The spiral bevel gear used in a helicopter gearbox is used as the research object. The gear blank parameters and machining parameters are listed in Table 1 and table 2 The input speed of the small wheel is 2933 R / min and the torque of the small wheel is 190 n · M. the load distribution calculation of spiral bevel gear adopts the method proposed in document, and its load distribution and maximum Hertz contact pressure are shown in Figure 1 It can be seen from Figure 1 that the maximum Hertz contact pressure decreases gradually from meshing in to meshing out (small wheel angle), and the maximum Hertz contact pressure is as high as 1.8 GPa According to the contact analysis model of spiral bevel gear, the relative sliding speed and coiling speed vector between the tooth surfaces of spiral bevel gear can be obtained by formula, and the results are shown in Fig. 2 It can be seen from Fig. 2 (a) that the entrainment speed increases gradually from meshing in to meshing out, and the entrainment angle decreases gradually The relative sliding velocity first decreases and then increases, the relative sliding velocity of the node is the smallest, and the angle direction of the sliding velocity changes suddenly near the node Through the contact analysis of spiral bevel gear, the contact load, velocity vector and tooth surface contact geometry are obtained. These parameters are input into the hybrid lubrication model, the lubrication analysis of spiral bevel gear can be carried out, and then the friction characteristics of meshing pair under lubrication can be obtained.

The mixed lubrication model of spiral bevel gear comprehensively considers the real contact situation of spiral bevel gear. Compared with the simplified algorithm and empirical formula, the solution time of friction coefficient is longer In order to compare the prediction of friction coefficient by empirical formula and accurate algorithm, Fig. 3 shows the calculation results in this paper and the calculation results of friction coefficient by Benedict and Kelley, Xu and Kahraman It can be seen from Fig. 3 that the friction coefficient of a pair of meshing pairs increases first and then decreases by using the mixed lubrication analysis model from meshing in to meshing out, and reaches the maximum at the node position It can be seen from Fig. 2 (b) that the relative sliding speed first decreases and then increases, while the change trend of friction coefficient is opposite When the contact lubrication zone of spiral bevel gear shows the characteristics of non-Newtonian fluid, the changes of contact load, relative sliding speed, entrainment angle, contact geometry (contact ellipse, radius of curvature, etc.) and oil film thickness will affect the shear rate and ultimate shear stress between tooth surfaces, so as to affect the friction characteristics between tooth surfaces The friction coefficient between the tooth surfaces of spiral bevel gears is affected by many complex parameters. Therefore, it is difficult to comprehensively consider the influence of various parameters when using empirical formula to predict the friction coefficient of spiral bevel gears in the whole meshing cycle.

The friction coefficient obtained by Benedict and Kelley, Xu and Kahraman is relatively small, and the friction coefficient results are quite different near the node This is because in Benedict and Kelley’s empirical formula, the relative sliding speed appears on the denominator, while the relative sliding speed near the node is small, resulting in a sudden increase in the friction coefficient at the node If the relative sliding speed is 0, the friction coefficient obtained by Benedict and Kelley’s empirical formula is infinite, which is not in line with the actual situation The friction coefficient obtained by Xu and Kahraman suddenly decreases near the node because Xu and Kahraman adopt the line contact model, that is, it is considered that the direction of entrainment velocity is consistent with the main direction of contact ellipse (short axis of contact ellipse) Zhang et al. Studied the friction characteristics of spiral bevel gear based on the point contact model in which the coiling speed coincides with the short axis of the contact ellipse. The prediction results of the friction coefficient are shown in Fig. 4 (a), which is similar to the prediction results obtained by Xu and Kahraman using the simplified algorithm, that is, there is a sudden change in the friction coefficient near the node Mohammadpour et al. Established a point contact lubrication model considering the included angle of entrainment speed, and analyzed the change of friction coefficient of spiral bevel gear meshing pair in one meshing cycle. The change of friction coefficient is shown in Fig. 4 (b) The model can simulate the contact of spiral bevel gear relatively truly, but the entrainment angle is assumed to be constant It can be seen from Fig. 4 (b) that the variation trend of friction coefficient of mohammadpour et al. Is similar to the results in this paper It can be seen from Fig. 3 that the empirical formula and simplified algorithm of friction coefficient have poor prediction accuracy near the node of spiral bevel gear, while the solution results of friction coefficient at other meshing positions are similar to the calculation results in this paper.