# Introduction helical gear is an common transmission device and has been widely used in all the fields, especially the machine under high speeds and heavy loads. The load distribution is the foundation of the gear preliminary design and strength rating process. It is known that the load distribution depends on the meshing stiffness of the tooth pair, and the load per unit of length is different at any point in the contact plane. The simple equations in standard ISO [1,2] to describe the load distribution, which is not in good agreement with experimental results. The contact lines of a helical gear are not parallel with the axial line, and the length of contact lines is dynamic changing in the meshing process. Some studies on the meshing stiffness and load distribution of involute gears can be found in technical literature. Z. Chen et al. [3][4][5][6] studied the tooth mesh stiffness and transmission error via the finite element method (FEM). However, this method is feasible but has the problem no generality of the obtained results. Afterward, J.I.Pedreroet al. [7][8][9] proposed a method to calculated non-uniform load distribution along the line of contact from the minimum elastic criterion potential, which depends on the transverse contact ratio. Through this method, the author analyzed the bending strength and pitting load capacity of helical gears. But the balanced load equation was not considered in this method. Thus it can't provide the load distribution of any point on the contact plane, and it is hard to locate the maximum value of the load. The heavy load and high-speed gear generate a lot of heat and temperature rise. High contact temperature of lubricant and tooth surfaces at the instantaneous contact position may lead to the breakdown of the lubricant film at the contact interface. The scuffing failure is unpredictable and fatal for the gears system. Therefor the scuffing load capacity is of great importance in the process of helical gear system preliminary design and strength rating, especially for the heavy load and high-speed gear system. ISO [10,11] provides two methods, namely the flash temperature method and integral temperature method, to evaluate the scuffing load capacity of the gear system, nevertheless the load distribution use the simplified form and the flash temperature calculated based on the Blok flash temperature equation [12], which can't get the accurate flash temperature and need to attach a large safety factor to amend it. The thermal elastohydrodynamic lubrication (EHL) is also a hot research topic. So far, most studies focused on the thermal EHL of spur gears. Wang and Cheng obtained a comprehensive research on a numerical simulation of the contact conditions of straight spur gear pairs [13,14]. L.M. Li proposed an inverse approach to establish the pressure, temperature rise, and apparent viscosity distribution in an EHL line contact [15]. Some methods and beneficial work to study on the EHL of spur gear, and a mass of cases were obtained, but these research mostly stay in the theoretical and ignore the practical application [16][17][18][19]. Besides the spur gears, the EHL research of helical gears is rare, no matter under the isothermal or thermal conditions. Recently, P.Yang and P.R.Yang used the multilevel multi-integration method to study the thermal elastohydrodynamic lubrication of tapered rollers in the opposite orientation; this model can be applied to the helical gear system [20]. These technical literature mentioned above mostly focus on the calculation method and directly offer the load; they can express the thermal EHL characteristic in theory but can't apply to A Global Journal of Researches in Engineering (A ) Volume Xx X Issue I Version I the actual conditions. The scuffing load capacity can be evaluated by the maximum contact temperature and minimum film thickness, which can be solved by the thermal EHL method. The literature of thermal EHL mostly focus on the temperature and film thickness of some single points [16][17][18][19][20]. The literature which makes the thermal EHL theory to design and check gear is absent. This paper proposes a method to study the load per unit of length distribution of all the points on contact plane of helical gears accurately based on the balanced Load equation, transmission error, and meshing stiffness. The feature coordinate to simplify the preliminary design and strength check process of helical gears is established. Based on the load distribution, the thermal EHL model, which corresponds more to actual conditions, is put forward. The hydrodynamic pressure, film thickness, and contact temperature, as well as the flash temperature, to check the scuffing load capacity are calculated via the numerical method. # II. # The Load Distribution Model of Helical Gear a) The contact model of the helical gear under heavy load For operating helical gear pair under heavy load, even though the driving pinion rotates at a constant speed, the gear as well as fall behind the angle ? ? than the theoretical location because of the deformation of driving and driven gear along the action line. As Fig. 1 shows, the meshing condition viewing from helical gear transverse direction, the profile of the solid line is the theoretical location, and the dashed line is the actual location under deformation. N 1 N 2 is the theoretical action line. Thus for any point K at the action line, the deformation along the action line is The analysis model of the helical gear is shown as Fig. 2, the contact plane N 1 N 2 N 3 N 4 is the tangent plane of two gear base circle. K 1 K 2 is one of the contact line. The actual action line is A 1 E 1 .The transverse contact ratio calculated by Equ.2. is the single contact tooth region. In double contact tooth regions, the contact lines always occur double in the same location. The coordinate system is established to describe the actual contact plane A 1 A 2 E 1 E2. The axial direction and action line direction are described by B and Î?" [10]. Î?" is the dimensionless parameter defined as follow: For helical gears, the contact lines at the same time always more than one (depend on the total contact ratio), so setting L denotes the sum of the length of all the contact lines, the sum load of the helical gear is the integral of w k , at any moment, it should be balanced to the extern load, the balanced load equation as follow: C N C N K N 1 1 1 / ) ( ? = Î?"sum b b L L K W r n P dl r l k dl w = = ? = ? ? 1 1 2 0 0 9549 ) ( ? (1) Where C is the pitch point, K the contact point, thus ] , [ 1 1 E A Î?" Î?" ? Î?" , ] , 0 [ b B ? . # Fig. 3: Analysis model of helical real action plane Now assuming A 1 is the gear pair approach point (begin meshing), the contact line will move along with the line A 1 E 2 to the recess action point E 2 . Axial contact ratio is expressed by the ? ? . The length distribution of contact lines with the parameters listed in Table 1 shows as Fig. 5. For different gear pairs, the distribution of load per unit of length is different. Compare the gear pair one and two or three and four, the face width doubled, and the length of the contact line almost double as well. The wider the tooth width, the longer the total length. But the value that the maximum minus the minimum has the upper limit value. We are letting the length ratio max min / L L = ? , the length ratio under different helix angle is show in Fig. 6. It is obvious that when 2 < ? ? , the length ratio is only 0.5, and the fluctuation of length is greatly, the maximum is as double as the minimum length. When 2 > ? ? , the length ratio is close to 1. # Fig. 6: The length ratio varying with the total contact ratio c) The tooth meshing stiffness The gear profile is a complex graphics and can be simply as the combination of a rectangle and a trapezium, as the Fig. 7 shows. The gear deformation includes the bending deformation, shear deformation, and contact deformation. The total is the sum of the rectangle deformation and trapezium deformation. The total of the contact deformation point along the action line can be expressed as the sum of the bending deformation B ? , shear deformation, S ? and lean deformation G ? . ? ? ? ? ? ? + ? = 3 ) ( cos 12 3 3 2 r r x r x f i x Br h h h h h S E w ? ? ? ? ? ? ? ? ? ? ? + + = x i r i r i r f i x S h h h h h h h S E w ln ) ( cos ) 1 ( 2 2 1 ? ? ? 2 2 cos 24 f i x x G S E wh ? ? ? = Where h is the height of the tooth, h x the height of contact point, h r the height of the rectangle, x ? the pressure angle of the contact point, a f a r f i S S S h hS h ? ? = ,v i notes the Poisson's ratio of pinion, E i the elastic modulus of pinion or gear, w the load per unit of length. When a pair of tooth meshing, the sum deformation along the action line can be described as the Equ.6? PV ? ? ? ? + + = ? 2 1 (6) Where E v v w PV ? ? )] 1 ( ) 1 [( 2 2 2 2 1 ? + ? = is the contact deformation of the meshing point, so the stiffness of this point can be expressed as Equ.7. ? = ? / w k (7) Fig. 8 shows the distribution of stiffness along the action line; the maximum occurs close to the pitch point. The regulation of stiffness is similar to the inverse unitary potential ) (? v in Ref. [7-9], which all reflect the capacity that the tooth bears the deformation. ? ? ? = sum L sum dl TE l k W 0 ) ((8) For the convenience of numerical calculation, the line A 1 E 1 is divided into N nodes, and ) 1 /( 1 1 ? = ? N E A l , so the total load can be expressed as Equ.9, the number of the contact line is m. b x ij m j n i ij sum l TE k W ? sin / ) 0 ), max(( 1 1 ? = = ? ? ? = ?? (9) The load distribution calculation process of helical gear is shown as Fig. 9. According to the equations above, we can calculate the load distribution of every point on the contact plane. Firstly, letting L=L 0 , k=k c , and TE 0 will be calculated; secondly, from the Equ.9, the sum load was obtained, then modified the TE according to Equ.10, until the sum load satisfied the Equ.11, then repeat this process until the results of all points were obtained. r i r i x i r i x i r i h h h h S E w ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? = ? ? The Load Distribution with Modification and Misalignment and Thermal Elastohydrodynamic Lubrication Simulation of Helical Gears © 2020 Global Journals Global Journal of Researches in Engineering (A ) Volume Xx X Issue I Version The contact model of the helical gear can be regarded as two tapered rollers [20]. In this model, the actual geometric model is two oval in oxz coordinate systems, and the geometric parameters as follows: According to the Hertz contact theory, the maximum contact press of the helical gear is expressed as Equ.14. # R Ew H ? ? ? 2 4 = (14) Where f) The load distribution calculation of helical gear without modification In the calculation process, the modification, setting ij ? to 0, then the transmission and the unit-line load of contact plane can be obtained. The threedimension unit-line load distribution is shown in Fig. 11. The load distribution on the contact plane is not only depends on the transmission error, but also depends on the stiffness distribution. The value of the load is small in the dedendum and addendum region of helical gear, and large close to pitch point. For the different helix angles and different face width, the load distribution is different. As we know, the sliding speed is maximum in the dedendum or addendum. From the 4 cases, the biggest occurs in the begin meshing and engaging-out point of, the helical contact plane. 2 1 2 1 R R R R R + = the curvature sum, E is equivalent to Young's modulus, ? ? ? ? ? ? ? ? ? + = 2 2 2 1 2 1 1 - 1 2 1 1 E E E ? ? . The transmission errors of helical gear pairs shown in Fig. 12, its distribution law corresponding to the distribution of length. The fluctuation of transmission error is the main indicator of the vibration and dynamic load. It is significant to choose the advisable helix angle and face width to make the fluctuation of transmission error minimum. The contact press distribution of the helical gears is shown in Fig. 13. It is similar to the unit-line load. But the contact stress of foot is greater than that of top of pinion. It is because the sum radius of the curvature of the root is smaller than that of the tooth top. For a helical cylindrical gear without modification, the maximum contact press is located at the engagement of the pinion root. As shown in Fig. 14, it is a fatigue pinion of an electric axle after the loading bench test, fatigue pitting at the meshing point of the root of the pinion. The contact stress with the modification and misalignment In practical application, the tooth surface of helical gear needs to be modified. For one reason, the maximum contact press located at the engagement of the pinion root. For the second reason, the shafts, bearings, and the housing will be deformed under the heavy load. The modification includes the profile modification and helix modification. Profile modification includes profile crowning, pressure angle modification, tip relief, and root relief. The helix modification includes lead crowning, helix angel modification, and end relief. According to the calculation method in Fig. 9, the transmission error(TE), unit-line load distribution, and contact press distribution of helical gear could be obtained. For the four gear pairs in table 1, the modification parameters are shown as table 2. In a transmission system, due to the deformation of the transmission shaft, the machining error of the gearbox, and the changes of the bearing stiffness, the gears pair will operate with misalignment. The contact state with misalignment can be simulated by the helix angle modification. Take gear pair 4 for example, the contact press distribution with different helix modification is shown in Fig. 17. With the increase of tooth inclination deviation, the contact stress inclines to one end, and the maximum stress is also increasing. The feature coordinate system From the above analysis, all the results in the contact plane were calculated by the numerical method. But in the preliminary design and strength rating process of helical gear, our focus is the maximum stress or the highest contact temperature of the contact plane. For the helical gear, the sliding speed is large in the addendum and dedendum region, and most scuffing failure occurs in these regions. The sliding speed approximately zero in the region close to the pitch point, so it is safer than the addendum and dedendum region. The 3-dimensional load distribution covers all the information about the contact plane, but it is timewasting and not intuitionistic. So the feature coordinate system should be established. # Fig. 18: Unit-linear load and transmission error distribution In Fig. 3, the line A 1 E 2 is recommended as the feature coordinate system, marked by ? , defined the same as Î?" . The A 1 is the approach point (begin meshing), and E 2 is the recess point. The parameters in feature coordinate cover both the tooth profile and axial information. The unit-line load distribution along with the feature coordinate is shown in Fig. 18.Compared the feature coordinate and 3-dimension coordinate, the maximum is the same in the addendum and dedendum region, the value in the feature coordinate system is a little lower than the three-dimension coordinate close to the pitch point, because the dangerous region of helical gear is the addendum and dedendum, so that the feature coordinate can satisfy the need of design and strength ratting. Furthermore, the feature coordinate is more succinct than the three-dimension. From the load distribution of gear pair 1 and 3, or gear pair 2 and 4, the helix angle has influenced a lot on the load distribution, in the case, the helix angle from the 9.8 to 20.2, and the maximum load per unit of length from the In the Ref.10, the load distribution of narrow helical gear is given as the Fig. 19, the buttressing effect near the end points A and E of the line of action, compared the results in this paper, the load distribution is similar with the gear pair 1,but is different from the gear pair 2, so the ISO can reflect the characteristic to some extent. Still it can't adapt to all the conditions. In this way, the method proposed by this paper is effective, and the point AU and EU is given by Equ. 15. b EU E A AU ? sin 2 . 0 = Î?" ? Î?" = Î?" ? Î?" (15) In the ISO standard theory, the helical gear with the total contact ratio When gears pair subjected to heavy loads, the lubricating film may not separate the surfaces adequately, leading to localized damage of the tooth surface. This type of failure is known as scuf fing and it is checked by maximum contact temperature and film thickness. Scuffing failure may occur at any stage during the lifetime of a set of gears, and it may lead to failure in a few of hours if the contact temperature is too high. So it is significant to research the scuffing load capacity of the gears. The high contact temperature is the main reason for scuffing failure, and the interfacial contact temperature conceived as the sum of two components, namely the bulk and flash temperature [10]. The bulk temperature fields are constant with time. Nevertheless, the flash temperature varies with time, and only appears in the interface of pinion and gear. In the thermal steady state, the bulk temperature of the tooth is higher than the ambient oil temperature. f M B t t t + = (16) Where t B is the contact temperature, t M the bulk temperature, and t f the flash temperature. The bulk temperature can be calculated by the finite element method (FEM). For the pinion of the pair 1, the bulk temperature is 360K (87?), the ambient oil temperature is 60?. The internal energy of oil film will increase due to the compression and viscous damping, and the temperature of oil film and tooth will rise because of the heat convection and conduction. After some time, the whole system will be on a thermal steady state. The bulk temperature of gear was higher than the ambient oil temperature. Therefore, in the inlet zone, the tooth profile is the heat source to heat the lubrication film, and the lubrication film temperature will rise instantly as high as the bulk temperature because that the lubricating film thickness is so thin. Then the temperature of the lubrication film will get higher due to the compression and viscous damping. In this region, the lubricating film is the heat source to heat the gear by convection and conduction. The ambient temperature was chosen as the inlet temperature in the technical literature [16][17][18][19][20]. By this means, the interface temperatures of the inlet contact region solved by the numerical method will lower than the interface bulk temperature; It is not in good agreement with the fact. In this paper, the interface bulk temperature was taken as the inlet temperature, which is more reasonable. # b) The thermal elastohydrodynamic lubrication equations i. Reynolds equation For the thermal steady state, the Reynolds equation of helical gear can be expressed in the following form as: ? ? ? ? ? ? ? ? ? ? ? ? ? ? = ? ? ? ? ? ? ? ? ? ? 0 1 1 2 2 2 ) ( ) ( F F x u u x h u x p F x ? ?(17) The Load Distribution with Modification and Misalignment and Thermal Elastohydrodynamic Lubrication Simulation of Helical Gears © 2020 Global Journals # Global Journal of Researches in Engineering (A ) Volume Xx X Issue I Version Where dz F h ? = 0 0 / 1 ? , dz z F h ? = 0 1 /? , dz F F z z F h ? ? ? = 0 0 1 2 ) / ( /? ? . The sliding speed of contact point: 30 / 2 , 1 i i i i R n u ? = = ?(18) In this equation, the mass density and the viscosity of lubricant can be described by Dowson and Higginson temperature-pressure-density relationship [21] and Roelands temperature-pressure-viscosity relationship [22] can be expressed as: ? ? ? ? ? ? ? ? ? ? × + × + = ? ? ) (? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? + + ? = ? 0 138 138 ) 1 ( 1 exp ) , ( 0 2 1 0 s z t t p A A t p ? ?(20) Where 67 . 9 ln 0 1 + = ? A , 9 - 2 10 1 . 5 × = A , ) /( 2 1 A A z ? = , )) 138 ( /( 0 1 0 ? = t A s ? . ii. The film thickness is given as ds s x x p E R x h h b a x x ? ? ? + = ln ) ( 4 2 2 0 ? (21) Whereh 0 is film thickness constant depend on the applied load. # iii. The balanced load equation The hydrodynamic pressure of the contact region must be in equilibrium with the input load, which calculated by the numerical method in section 2.Thus the pressure equilibrium is expressed by the following equation: ? = b a x x dx x p w ) ((22) The three-dimensional temperature distribution of oil film was calculated by the energy equation as follows: ? ? ?? ? ?? ?? ? ? ? ? ?? ?? ? ion FluidFrict Expation n Compressio Conduction H p z u x p u t t z t x t u c 2 / eat 2 2 Convection ) ( ) ( ? ? ? ? ? ? ? ? = ? ? ? ? + ? ? ? ? ? ? ? ? ? ? ? (23) The boundary condition for Eq.23 is the thermal interface equation: The equation governing the shear stress in the lubricant film is a momentum equation which for the Newtonian fluid is as follow: ? ? ? ? ? ? ? + ? ? ? = + ? ? ? = ? ? ? ? = ? ? = x b h z x a z t s x ds z t u c h x t t s x ds z t u c x t? ? ? ? ? ? ? ? ? ? = ? ? z u z x p ? (24) iv. Numerical method All equations and boundary conditions mentioned above about thermal were written into dimensionless forms to facilitate the numerical calculation. There are pressure loop and temperature loop in the entire calculation process. In the pressure loop, assuming the temperature of all the nodes is the ambient temperature (360K). In the contact coordinate in the rolling direction ( x direction), the whole contact region (from inlet position to outlet position) was divided into 81 non-equidistant nodes; the interval is large in the inlet region and small in the secondary pressure peak region. The distribution of pressure and the film thickness can be calculated by the Newton-Raphson method depending on the Equ.17-24.Divided the film gap (obtained in the pressure loop, z direction) into 21 equidistant nodes. The temperature of all the nodes can be calculated depend on the equ.20-21 and equ.24-25.Then repeat these two processes until the pressure error less than 5 -10 .The lubricant properties are shown in Table 2. N/m shown in Fig. 14. For the gear transimission system under heavy load, the film center pressure is higher; in this case, the maximum pressure is as high as 1.51GPa, the minimum film thickness is 0.42 m µ . Because of the heavy load, the curve of oil pressure and oil thick has't the typical characteristic of thermal EHL. The secondary pressure peak close to the outlet and inconspicuous. The distribution of the pressure is similar to Hertzian pressure distribution. The maximum occurs the location x=0 and the heavier of the load, the higher of the maximum. The heavier, the wider the contact zone. Furthermore, it has little influence on the minimum oil thickness. The three-dimension and boundaries temperature distribution of the contact region are shown in Fig. 21. The film center temperature is much higher than two boundaries; this is because the heat conducts from the film to all the teeth to maintain the bulk temperature. The distribution of film center temperature is similar to the film center pressure and peaked in the Hertzian maximum, but two boundaries are different from the film center; the temperature increase until close to the outlet and slightly decrease. The reason is in the rear part of the contact zone, the compression effect is no longer generate the heat, so the heat from the film center to the two boundaries via conduction effect and to the inner tooth through convection. # b) Comparison with the Blok's flash temperature The contact and flash temperature is the main reason for the gear scuffing failure, Blok derived the flash temperature equation in 1937 [12] as follow: 0 2 2 2 2 1 1 1 1 2 1 ) ( 7858 . 0 b u c u c u u fw t f ? ? ? ? + ? = (25) Where f=0.06 is coefficient of friction, w the load per unit of length, In the thermal EHL results, the maximum temperature rise of two boundaries is its flash temperature. The flash temperature along the feature coordinate was calculated by the thermal EHL method and Blok equation as show in Fig. 22. In the dedendum region of the pinion (close to approach point), the thermal EHL flash temperature is higher than the Blok result, and in the addendum region of the pinion(close to recess point), the Blok temperature is lower than thermal EHL flash temperature. In theory, the temperature-pressure-density and temperaturepressure-viscosity effect and the compression of oil film were considered in the thermal EHL theory. The Blok flash temperature equation is very concise and corresponding to the experiment to some extent. However, it loses sight of the influence of the bulk temperature, so it is not a good agreement of the experiment under the heavy load. Thus the thermal EHL theory is even close to the actual condition. The safety factor defined as follow [10]: Where t s is the scuffing temperature, t oil the ambient temperature, t cmax the maximum contact temperature. When gear teeth are separated completely by a full fluid film of lubricant, there is no contact between the asperities of tooth surfaces, and usually, there is no scuffing and wear. For the heavy load helical gear, the thickness of oil film is very small, and incidental asperity contact takes place. As the minimum film thickness decreases, the number of contacts increases, abrasive wear, adhesive wear, and scuffing became possible. So the minimum thickness of lubricant film is a property of scuffing load capacity, especially for the gear pair under heavy load. The minimum thickness equation was given by Dowson and Higginson [22] as the Equ.27. The safe factor can be measured by the thickness ratio. When the thickness ratio ) 4 , 1 ) ? ? ,it is considered as mixed friction. When the thickness ratio is less than 1, the scuffing failure probably takes place to a great extent. 3![Load Distribution with Modification and Misalignment and Thermal Elastohydrodynamic Lubrication Simulation of Helical Gears Nomenclature b gear face width, mm ? viscosity of lubricant, Ns/ load per unit of length, N/m 2 b r base radius of gear, m n m normal module, mm 0 b half width of the Hertzian line contact, m ? ? transverse contact ratio H P maximum Hertzian pressure, Pa ? ? axial contact ratio ? Braus viscosity-pressure coeffient,](image-2.png "The 3 mw") ![heat of lubricant and solid, J/kgK 0 t ambient temperature, K bt P transverse circular pitch, m Year 2020 L length of all the contact lines, m TE transmission error, ?m](image-3.png "") 1![Fig. 1: The contact model of helical gear in transverse under heavy load b) The sum length of the contact linesThe analysis model of the helical gear is shown as Fig.2, the contact plane N 1 N 2 N 3 N 4 is the tangent plane of two gear base circle. K 1 K 2 is one of the contact line. The actual action line is A 1 E 1 .The transverse contact ratio calculated by Equ.2.](image-4.png "Fig. 1 :") 2![Fig. 2: The analysis model of the helical gear Viewing from the transverse direction, for](image-5.png "Fig. 2 :") ![divide the A 1 E 1 into three regions, A 1 B 1 , B 1 D 1, and E 1 D 1 , the A1B1 and D1E1 are the double contact tooth regions (the number of tooth contact at the same time is 2), B 1 D 1 is the single contact tooth region. So we can divide the actual contact plane into three regions. It shows as the Fig.3. The A 1 B 1 A 2 B 2 and D 1 D 2 E 1 E 2 are the double contact tooth regions, B 1 B 2 D 1 D 2](image-6.png "") ![Load Distribution with Modification and Misalignment and Thermal Elastohydrodynamic Lubrication Simulation of Helical Gears Global Journal of Researches in Engineering (A ) Volume Xx X Issue I Version I Here assuming the stiffness of the contact point TE is k, then the load per unit of length of point TE is](image-7.png "The.") 5![Fig. 5: The result of the sum length of contact lines](image-8.png "Fig. 5 :") 7![Fig. 7: The stiffness analysis model of gear tooth](image-9.png "Fig. 7 :") 8![Fig. 8: Gear teeth stiffness along the action line of gear pair 1 d) The load distribution calculation process of helical gear In Fig.3, the sum length is shown as follow:](image-10.png "Fig. 8 :") 910![Fig. 9: The load distribution calculation process of helical gear e) The contact stress of the helical gears](image-11.png "Fig. 9 :Fig. 10 :") ![the subscript 1 means the pinion, 2 for gear. The curvature of the contact point of two oval are as follows:](image-12.png "") 11![Fig. 11: Unit-line load distribution of helical gearsThe load distribution solved in this paper is in good agreement with the results which occur in Ref.8 by the finite elements method (FEM). The results show that the contact stress solved by FEM is larger near the pitch point and small in the dedendum and addendum region. Thus the load distribution in this paper is consistent with the actual condition to a great extent without the consideration of the end effect, the machine error, and assembly error.](image-13.png "Fig. 11 :") 12![Fig. 12: The transmission error of the helical gears](image-14.png "Fig. 12 :") 1314![Fig. 13: The contact press of the helical gears](image-15.png "Fig. 13 :Fig. 14 :") 15![Fig. 15: Modification distribution of contact plane The modification distribution of the helical gear contact plane is as shown in Fig.15. And the contact press of helical gear is shown in Fig.16. After modification, the contact press distribution of the tooth surface has been greatly improved. The maximum pressure transferred from the meshing point of the pinion root to the central area of the contact surface.](image-16.png "Fig. 15 :") 16![Fig. 16: The contact press of the helical gears with modification](image-17.png "Fig. 16 :") 17![Fig. 17: The contact press distribution with helix modification](image-18.png "Fig. 17 :") 19![Fig. 19: Load sharing factor for the narrow helical gear and wide helical gear of ISO](image-19.png "Fig. 19 :") ![to spur gears as well as the buttressing effect. And the wide helical gear ( to act near the ends A and E along with the helix teeth over a constant length. From the Fig.6, we know that even though the 2 > ? ? , the fluctuation is considerable for some parameters. So the ISO standards are stereotypical and inaccurate.](image-20.png "") ![Where t a and t b are the bulk interface temperature of the pinion and gear.](image-21.png "") 220![Fig. 20: The pressure and film thickness of oil film Fig.20 shows the oil pressure and film thickness along the x direction of began meshing point of gear pair one and two. The load per unit of length of gear pair one is](image-22.png "? and 2 ?Fig. 20 :") 21![Fig. 21: Three-dimension film temperature and film center and interface temperature of gear pair1](image-23.png "Fig. 21 :") ![the heat conduction coefficient, density, specific heat, sliding speed of the pinion and gear,b 0 the half-width of Hertzian line contact.](image-24.png "") ![with Modification and Misalignment and Thermal Elastohydrodynamic Lubrication Simulation of Helical Gears © 2020 Global Journals Global Journal of Researches in Engineering (A ) Volume Xx X Issue I Version](image-25.png "") 22![Fig.22:The TEHL flash temperature and Blok flash temperature of the helical gear in the feature coordinate c) The minimum thickness method from thermal EHL and Dowson equation When gear teeth are separated completely by a full fluid film of lubricant, there is no contact between the asperities of tooth surfaces, and usually, there is no scuffing and wear. For the heavy load helical gear, the thickness of oil film is very small, and incidental asperity contact takes place. As the minimum film thickness decreases, the number of contacts increases, abrasive wear, adhesive wear, and scuffing became possible. So the minimum thickness of lubricant film is a property of scuffing load capacity, especially for the gear pair under](image-26.png "Fig. 22 :") 23![Fig. 23: The minimum thickness of the film in the feature coordinate system](image-27.png "Fig. 23 :") 1Gear Pair 1Gear Pair 2Gear Pair 3Gear Pair 4Number of teeth (pinion/gear)23/3023/3023/3023/30Normal pressure angle20202020Normal module, mm3333Face width, mm20402040Input power(kw)50505050Modulus of elasticity ?? Helix angle,°9.89.820.220.21 /?? 2 , GPa 206/206 206/206 206/206 206/206 Poisson ratio, ?? 1 /?? 2 0.3/0.3 0.3/0.3 0.3/0.3 0.3/0.3 2Modification formModification/umprofile crowning(barreling) ?? ??10Lead crowning ?? ??10tip relief5End relief0Helix angle modification fH?0Pressure angle modification?fH?0 2Symbol,unitValueAmbient viscosity of lubricant , 0 ? ,Ns/ 2 m0.08Specific heat of lubricant, c, J/kgK2000Specific heats of solids,J/kgK470Thermal conductivity of lubricant, ? , W/mK0.14Thermal conductivities of solids a and b, 1 b r base radius of gear, m ba © 2020 Global Journals x i f x Bt h h h h h h h h h h © 2020 Global Journals g) © 2020 Global Journals h) The minimum film thickness solved by thermal EHL and Dowson equation is shows as Fig. 23. The minimum thickness curves are similar along the feature coordinate of gear pair one and two. The minimum film thick is lowest in the dedendum of the pinion, so the dedendum of the pinion is dangerous. The results are consistent with the maximum flash temperature method. The value from thermal EHL is lower than that from the Dowson equation. The Dowson equation didn't consider the influence of temperature on film thickness. In contrast to the fatigue damage, a single momentary overload may initiate scuffing failure, so the scuffing capacity is crucial for the heavy load and highspeed helical gear system. The minimum film thickness and maximum contact temperature, as well as the flash temperature obtained by the thermal EHL theory, are applied to check the scuffing strength. The scuffing capacity can be checked by the flash temperature or minimum thickness method. V. ## Conclusions In this paper, a model of non-uniform along the contact line of helical gear teeth, obtained from the meshing stiffness, transmission error, and balanced load equation, has been proposed. The feature coordinate system was established to simply preliminary design and strength rating process. The thermal effect of helical gear lubrication will affect the scuffing load capacity. The thermal EHL method is applied to check the scuffing load capacity. The following conclusions can be drawn: 1. The length of contact lines depends on the face width and basic helix angle. And the length is oline is not paralleled with the axial line. The fluctuation is evident when the total contact ratio less than 2.The distribution of transmission error is similar to the sum length. The distribution of load per unit of length is dependent on the meshing stiffness and transmission, it is similar to the meshing stiffness along the contact line and similar to the transmission error along the Baxis. The modification of gears can greatly improve the distribution of contact stress. The maximum contact stress is concentrated in the center of the tooth surface to avoid the contact between the tooth surface. The contact stress inclines to one end, and the maximum stress is also increasing under with misalignment. 2. The feature coordinate system established in this paper can cover bother the tooth profile and axial information. It can satisfy the need for preliminary design and strength check of helical gears. The weakest strength region is the dedendum of the pinion. A load of the dedendum region in the feature coordinate system is as same as the threedimension coordinates system. 3. The thermal EHL theory can provide more information to check the scuffing load capacity. The highest contact temperature and minimum film thickness method were applied to evaluate the scuffing load capacity. The pressure distribution under heavy load is similar to the Hertzian pressure distribution, and the secondary presser peak close the outlet and inconspicuous. The temperature of the film center is much higher than two boundaries. The film center temperature distribution is similar to the pressure distribution. The two boundaries temperature increase until close to the outlet and slightly decrease. The thermal EHL flash temperature is corresponding to the Blok flash temperature but higher in the dedendum region and lower in the addendum region of the pinion. The minimum thickness film is smaller than the Dowson equation. So the scuffing load capacity that solved by thermal EHL is more practical than the traditional method. * Calculation of Load Capacity of Spur and Helical Gears -Part 2: Calculation of Surface Durability (Pitting), International Organization for Standardization ISO Standard 6336-2:1996 1996 Geneva, Switzerland * Calculation of Load Capacity of Spur and Helical Gears -Part 3: Calculation of Tooth Bending Strength ISO Standard 6336-3:1996 International Organization for Standardization 1996 * Mesh stiffness calculation of a spur gear pair with tooth pro file modification and tooth root crack. Mechanism and Machine Theory ZaigangChen YiminShao 2013 62 * A model for the study of meshing stiffness in spur gear transmissions. Mechanism and Machine Theory AFernandez Del Rincon FViadero MIglesias 2013 * Analysis of tooth contact and load distribution of helical gears with crossed axes. Mechanism and Machine Theory YZhang ZFang 1999 34 * Gear tooth stiffness reduction measurement using modal analysis and its use in wear fault severity assessment of spur gears Yesilyurt AndrewDGu Ball NDT&E International 36 2003 * Load distribution model along the line of contact for involute external gears, Mechanism and Machine Theory JIPedrero MPleguezuelos MArtés JAAntona 2010 45 * Critical stress and load conditions for bending calculations of spur and helical gear JIPedrero MPleguezuelos MMuñoz International Journal of Fatigue 48 2013 * Critical stress and load conditions for pitting calculations of spur and helical gear JIPedrero MPleguezuelos MMuñoz Mechanism and Machine Theory 46 2011 * Calculation of scuffing load capacity of cylindrical, bevel and hypoid gears-Part 1: Flash temperature method 13989-1:2000 ISO Standard London The United Kingdom 2000 * Calculation of scuffing load capacity of cylindrical, bevel and hypoid gears-Part 1: Integral temperature method ISO Standard 13989-2:2000 2000 The United Kingdom London * Theoretical study of temperature rise at surfaces of actual contact under oiliness lubricating conditions HBlok Proceedings of the General Discussion on Lubrication and Lubricants[J]. Imech E the General Discussion on Lubrication and Lubricants[J]. Imech E 1937 2 * Anumerical solution to the dynamic load, film thickness and surface temperatures in spur gears, part I analysis KLWang HSCheng Journal of Mechanical Design 103 1981 * A Numerical solution to the dynamic load, film thickness and surface temperatures in spur gears, part II results KLWang HSCheng Journal of Mechanical Design 103 1981 * Jaw-Ren Lin, dtal. Inverse approach for calculating temperature in EHL of line contacts Li-MingChu Hsiang-ChenHsu Tribology International 42 2009 * The Navier-Stokes approach for thermal EHL line contact solutions TAlmqvist RLarsson Tribology International 35 2002 * Effect of lubricant selection on EHL performance of involute spur gears PAnuradha PunitKumar Tribology International 50 2012 * Thermal elastohydrodynamic simulation of involute spur gears incorporating mixed friction LarsBobach RonnyBeilicke DirkBartel Tribology International 2012 * Transient thermo elastohydrodynamic lubrication analysis PWWang HQLi JWTong Tribology International 37 2004 * Analysis on the thermal elastohydrodynamic lubrication of tapered rollers in opposite orientation PingYang PeiranYang Tribology International 40 2007 * Correlation aspects of Viscositytemperature-pressure Relationship of Lubrication oils CjaRoelands 1966 Netherlands Delft Univeristy of Technology * Elastohydrodynamic Lubrication DDowson GRHigginson 1977 Pergamon Press Oxford