Optimization of lubrication characteristics of wind turbine’s transmission system based on Newton Raphson method

Cao, Li; Sun, Wenlei; Gou, Tao

doi:10.21595/jve.2024.23725

Journal of Vibroengineering

Browse Journal

Submit article

Published: 02 April 2024

Check for updates

Optimization of lubrication characteristics of wind turbine’s transmission system based on Newton Raphson method

Li Cao¹

Wenlei Sun²

Tao Gou³

^{1, 2, 3}School of Mechanical Engineering, Xinjiang University, Urumqi, China

Corresponding Author:

Li Cao

Cite the article Download PDF

Downloads 98

Abstract

To improve the reliability and lifespan of wind turbines, this paper takes the two-stage fixed shaft gearbox experimental platform of wind turbines as the research object. Based on Hertz contact theory, the oil film pressure and thickness in the contact area are solved by combining the equations of elastohydrodynamic lubrication and the Newton Raphson method; And the lubrication characteristics of the transmission system were analyzed to verify the correctness of the method; At the same time, in response to partial load phenomenon caused by system coupling deformation, genetic algorithm was selected to modify the gear teeth. The results show that the max unit load on the tooth face and the maximum stress of tooth root decreased by up to 26.48 % and up to 20.35 % respectively after modification which can improve the uneven distribution of oil film and the lubrication performance of the tooth surface.

Optimization of lubrication characteristics of wind turbine’s transmission system based on Newton Raphson method

Highlights

This article presents a computing method to obtain the oil film pressure and thickness in the contact area which combining the equations of elastohydrodynamic lubrication and the Newton Raphson method.
The lubrication characteristics of the transmission system were analyzed to verify the correctness of the method.
In response to partial load phenomenon caused by system coupling deformation, genetic algorithm was selected to modify the gear teeth, improving the uneven distribution of oil film and the lubrication performance of the tooth surface.

1. Introduction

With the continuous progress of industrial level, the size of wind turbines continues to increase, and the relatively harsh working environment of wind turbines also leads to an increase in corresponding faults. Usually, about 20 % of wind turbine faults are caused by gearbox failure, and about 60 % of faults in gearbox come from gears [1]. The main forms of gear failure are pitting, wear, fatigue, bonding, etc. These failure situations can be effectively improved through lubrication to a certain extent. Therefore, understanding, understanding, and mastering the dynamic characteristics of gear lubrication systems and tooth surface lubrication contact performance have certain research significance for studying the vibration and noise of transmission systems and improving transmission performance.

Good lubrication performance can effectively reduce the wear and contact fatigue of interacting tooth surfaces, and improve their durability [2]. Usually, the contact surfaces of high gear pairs and other high gear pairs will form elastohydrodynamic lubrication under high pressure. The lubricating oil film can effectively separate the two contact surfaces, thereby reducing tooth friction and wear, reducing vibration and noise [3].

The transition from classical lubrication theory to elastohydrodynamic lubrication theory has played a driving role in the research progress of gear lubrication to a certain extent. Reasonably designing gear lubrication systems can help improve their service life [4]. Scholars at home and abroad have established a certain theoretical foundation for research. H. P. Evans et al. discussed the impact of oil film thickness between tooth surfaces less than the roughness feature height on lubrication characteristics [5]. S. Li et al. proposed a transient non Newtonian mixed elastohydrodynamic lubrication (EHL) model for involute spur gear tooth contact, which differs from the contact between two cylindrical rollers. The contact of spur gears undergoes many time-varying contact parameters, including normal load, curvature radius, and suction speed. The lubrication behavior under rough conditions under different operating conditions was studied [6]. In 2014, CaichaoZhu and MingyongLiu et al. considered the geometric changes and non Newtonian properties of helical gear pairs and established a thermal finite line contact elastohydrodynamic model for helical gear pairs. They applied the generalized Reynolds equation and multigrid method to solve the model. The analysis showed that the main difference between Ree Eyring non Newtonian fluid and Newtonian fluid was the membrane temperature and shear stress [7]. Mankoun et al. concluded in their analysis that the lubrication effect of the transmission system will directly affect the dynamic characteristics of the transmission system [8]. Guan Hao et al. used Romax software to perform optimization analysis of planetary gear modification, and compared and analyzed the vibration response and transmission error before and after modification [9]. Miryam B. Sánchez et al. determined the modification of spur gears through meshing stiffness, load distribution, and transmission error (including the type of modification curve, modification position, and modification amount), which has higher accuracy compared to the modification results obtained based on the latest meshing stiffness equation [10]. J. Bruyère et al. obtained the modification equations for spur and helical gears through numerical analysis, which can effectively control the transmission error, original stress, etc. of gears [11]. However, high gear tooth surface accuracy, large computational complexity, and relative complexity are required.

This article studies the lubrication characteristics of wind power gear transmission systems and conducts micro modification research on the tooth surface to address the phenomenon of eccentric load in spur gear transmission systems. Reasonable modification topology parameters are obtained, and various contact performance parameters before and after modification are compared and analyzed. To a certain extent, it can improve the contact characteristics of the tooth surface, improve the smoothness of the transmission system, reduce vibration noise, and extend the service life of the system.

2. Elastohydrodynamic lubrication model for gear line contact

2.1. Quasi static isothermal elastohydrodynamic lubrication equation

Gear meshing usually operates in an environment containing lubricating oil, so the meshing tooth surface is not in direct contact. There is a thin film of lubricating oil between the two tooth surfaces, which is called “elastohydrodynamic lubrication contact”. According to the Hertz contact theory, gear tooth meshing is considered as the contact between two parallel cylindrical axes, which are compressed against each other under the action of load without lubrication. The contact surface of the gear undergoes elastic deformation, and the actual contact area is not a line but a narrow surface, called the Hertz contact area, with its half width denoted as $B$ [12].

The calculation formula for the half width B of the Hertz contact zone:

1

B = \sqrt{\frac{8 w R}{π E^{'}}},

2

\frac{1}{E^{'}} = \frac{1}{2} ∙ (\frac{1 - μ_{1}^{2}}{E_{1}} + \frac{1 - μ_{2}^{2}}{E_{2}}),

where $w$ is the load per unit length; $R$ is the equivalent radius of curvature; $E$ is the equivalent elastic modulus, and $E_{1}$ , $E_{2}$ , $μ_{1}$ and $μ_{2}$ are the elastic modulus and Poisson’s ratio of the two gears, respectively.

In the Hertz contact area, the contact pressure $p$ is distributed in an elliptical pattern, namely:

3

p = p_{H} \sqrt{1 - {(\frac{x}{B})}^{2}}, p_{H} = \frac{2 w}{π B} = \sqrt{\frac{w E^{'}}{2 π R}} = \frac{E^{'} B}{4 R},

where $p_{H}$ is the maximum Hertz pressure.

By linking Reynolds lubrication theory with Hertz contact theory, the film thickness formula for line contact elastohydrodynamic can be obtained. After considering lubrication, the lubricating oil undergoes a viscous pressure effect under pressure and the gear contact surface undergoes elastic deformation, as shown in Fig. 1.

Fig. 1Deformation diagram of gear oil film thickness

Fig. 1 shows that the coordinate system of oil film in the lubricating oil film. The $x$ -axis of the oil film coordinate system is the flow direction of the lubricating oil film, the $y$ -axis is the oil film thickness direction, and the $z$ -axis is the length direction of the oil film. The length of the oil film along the $z$ -axis of the oil film is $L$ , and the width of the oil film along the $z$ -axis is 2 $B$ . The $y$ -coordinate of the oil film is much smaller than the other two coordinate values, and it can be simplified as a two-dimensional problem, which can be regarded as a line contact elastohydrodynamic lubrication state [13].

The elements of oil film are only affected by the fluid pressure $p$ and viscous force $τ$ . The velocity of the element along the coordinate axis direction are $v_{x}$ , $v_{y}$ , and $v_{z}$ , respectively. And $\frac{\partial v_{z}}{\partial z} =$ 0, $\frac{\partial p}{\partial z} =$ 0. Then the steady-state isothermal Reynolds equation for line contact is:

4

\frac{\partial}{\partial x} [\frac{ρ h^{3}}{12 η} \frac{\partial p}{\partial x}] = (v_{0} + v_{h}) \frac{\partial (ρ h)}{\partial x} + \frac{\partial (ρ h)}{\partial t},

where $x$ is the coordinate of the direction along which the lubricating oil film flows; $ρ$ is the density of lubricating oil; $h$ is the thickness of the lubricating oil; $η$ is the dynamic viscosity of the lubricating oil; $p$ is the oil film pressure; $v_{0}$ and $v_{h}$ are the entrainment velocity of the oil film in the $x$ -direction when $z =$ 0 and $z = h$ , respectively.

When oil film is in a steady state, $\frac{\partial h}{\partial t} =$ 0, $v = v_{0} + v_{h}$ , then:

5

\frac{\partial}{\partial x} [\frac{ρ h^{3}}{η} \frac{\partial p}{\partial x}] = 12 v \frac{\partial (ρ h)}{\partial x} .

In the actual working condition of elastohydrodynamic lubrication, the cylindrical surface often undergoes deformation in the normal direction which causing the shape of the gap is changed. Supposing the sum of the elastic deformations of the two surfaces along the normal direction is $A (x)$ . Therefore, the formula for the oil film thickness at any point $x$ when the elastic cylinder contacts is:

6

h (x) = h_{0} + \frac{x^{2}}{2 R} + A (x),

where $h_{0}$ is the central film thickness without deformation; $R$ is the equivalent radius of curvature; $A (x)$ is the elastic deformation displacement caused by pressure.

According to the theory of elasticity, the elastic displacement $A (x)$ of each point on the surface in the vertical direction is:

7

A (x) = - \frac{2}{π E^{'}} ∙ \int_{e_{1}}^{e_{2}} [p (e) ∙ l n {(e - x)}^{2}] d e + c,

where $e$ represents the distance between any line load $p (e) d e$ and the coordinate origin; $e_{1}$ and $e_{2}$ are the starting and ending coordinates of load $p (x)$ , respectively; $E^{'}$ is the equivalent elastic modulus; $c$ is an undetermined constant.

The density of lubricating oil is a function of pressure and temperature. When the pressure on the lubricating oil increases, its volume decreases, resulting in an increase in density. The change in density with pressure can be represented by the compression coefficient $k$ , which is:

8

k = \frac{1}{ρ} ∙ \frac{d ρ}{d p} = - \frac{1}{V} ∙ \frac{d V}{d p} = - \frac{1}{V} ∙ \frac{V - V_{0}}{p - p_{0}} = \frac{\frac{V_{0}}{V} - 1}{p - p_{0}} = \frac{\frac{ρ}{ρ_{0}} - 1}{p - p_{0}},

where $V_{0}$ , $p_{0}$ , $V$ and $p$ are the volumes and pressures of lubricating oil with a mass of m at pressure $p_{0}$ and $p$ , $ρ_{0}$ and $ρ$ are the densities at pressures $p_{0}$ and $p$ , respectively. Therefore:

9

ρ = ρ_{0} (1 + k (p - p_{0})) .

The formula for the compression coefficient k of lubricating oil is:

10

k = (7.25 - l g η) \times 10^{- 10},

where $η = η_{0} \exp \{(\ln η_{0} + 9.67) [{(1 + 5.1 \times 10^{- 9} p)}^{a} - 1]\}$ , $η_{0}$ and $η$ are the dynamic viscosity under pressure $p_{0}$ and $p$ , respectively, and $ɑ$ is the Reynolds viscosity coefficient.

The isothermal density equation is:

11

ρ = ρ_{0} [1 + \frac{p \times 0.6 \times 10^{- 9}}{1 + p \times 1.7 \times 10^{- 9}}] .

The $p (x)$ is the pressure at the oil film, and the load per unit tooth width can be obtained by integrating, the formula is:

12

w = \int_{x_{i n}}^{x_{o u t}} p (x) d x,

where $x_{i n}$ is the coordinate of the entrance of the lubrication area; $X_{o u t}$ is the coordinate at the outlet of the lubrication area.

2.2. Numerical solution of elastohydrodynamic lubrication equation

In order to improve the stability and increase the convergence speed of the calculation process [14], the number of parameters is reduced by introducing dimensionless parameters. The formulas of dimensionless coordinates, pressure, oil film thickness, density, viscosity, material, velocity, and load are:

13

\bar{x} = \frac{x}{B}, \bar{p} = \frac{p}{p_{H}}, \bar{h} = \frac{h R}{B^{2}}, \bar{ρ} = \frac{ρ}{ρ_{0}}, \bar{η} = \frac{η}{η_{0}},

\bar{G} = α E^{'}, \bar{v} = \frac{v η_{0}}{(E^{'} R)}, \bar{w} = \frac{w}{(E^{'} R)} .

Dimensionless equation of line contact steady-state isothermal Reynolds is:

14

\frac{d}{d \bar{x}} [\bar{ε} \frac{d \bar{p}}{d \bar{x}}] = \frac{d (\bar{ρ h})}{d \bar{x}}, \bar{ε} = \frac{\bar{ρ} {\bar{h}}^{3}}{\bar{η} λ}, λ = \frac{3}{4} π^{2} \frac{\bar{v}}{{\bar{w}}^{2}} .

Dimensionless equations of oil film thickness equation, density equation, viscosity equation, load balance equation are:

15

\bar{h} = {\bar{h}}_{0} + \frac{{\bar{x}}^{2}}{2} - \frac{1}{π} \int_{{\bar{x}}_{im}}^{{\bar{x}}_{out}} p ({\bar{x}}^{'}) \ln |\bar{x} - {\bar{x}}^{'}| d {\bar{x}}^{'},

16

\bar{ρ} = 1 + \frac{0.6 \times 10^{- 9} \bar{p} p_{H}}{1 + 1.7 \times 10^{- 9} \bar{p} p_{H}},

17

\bar{η} = e^{(\ln η_{0} + 9.67) [{(1 + 5.1 \times 10^{- 9} \bar{p} p_{H})}^{z_{0}} - 1]},

18

\bar{w} = \int_{{\bar{x}}_{i n}}^{{\bar{x}}_{o u}} \bar{p} d \bar{x} = \frac{π}{2} .

Dimensionless Equations were resolved using a forward difference format and combine the Newton Raphson numerical calculation method and introducing a “downhill factor” to accelerate the convergence speed of the iterative process [15]. There are 264 nodes are selected. Calculate the domain boundary entry $x_{i n} =$ –2.5, boundary exit $x_{o u t} =$ 1.5. After dispersing, the Reynolds equation, oil film thickness equation, viscous pressure equation, density equation, and load equation as follows:

19

f_{i} = h_{i}^{3} {(\frac{d p}{d x})}_{i} - A {\bar{η}}_{i} (k_{i} - \frac{{\bar{ρ}}_{e} h_{e}}{{\bar{ρ}}_{i}}) = 0,

20

h_{i} = h_{0} + \frac{x_{i}^{2}}{2} + \sum_{j = 1}^{n} k_{i j} p_{j},

21

η = η_{0} \exp \{(\ln η_{0} + 9.67) [{(1 + \frac{p}{p_{0}})}^{z}]\},

22

ρ = ρ_{0} (1 + \frac{0.6 p}{1 + 1.7 p}),

23

w = \sum_{i = 1}^{N} P_{i} Δ x = \frac{π}{2} .

Differentiate the above formulas to gain increment of each derivative matrix. Then these equations were solved by iteration until the convergence accuracy is met, with the oil film thickness convergence accuracy of 0.001 and a pressure convergence accuracy of 0.01.

3. Research on lubrication characteristics of transmission systems

3.1. Oil film thickness and pressure distribution at key meshing points of gear pairs

The two-stage parallel shaft spur gear transmission system of the DDS power transmission fault diagnosis platform which produced by SQI company. The transmission system consists of two-stage gear pairs, three transmission shafts and six deep groove ball bearings. The three-dimensional model and schematic diagram of two-stage gear transmission system is shown in Fig. 2. The parameters of transmission systems is shown in Table 1. During the transmission systems process, gears’ loads varies at different positions along the meshing line. So the load distribution is necessary. The coefficient of load distribution diagram along the meshing line of the gear without considering the transmission error as shown in Fig. 3. Therefore, the contact ratio of gear pair 1 is 1.7667 according Eq. (24). It shows the single tooth meshing time is relatively short. Moreover, the gear in the system has not been modified:

24

ε = \frac{L_{A c t u a l m e s h i n g l i n e}}{p_{b}} = \frac{L_{A c t u a l m e s h i n g l i n e}}{(p c o s α)} .

Fig. 2The three-dimensional model and schematic diagram of two-stage gear transmission system

a) Three-dimensional model

b) Schematic diagram

Table 1The parameters of transmission systems

Gear	The number of teeth	Module $p$	Mesh angle $α$ / °	The width of tooth	Mass / kg	Moment of inertia / kg∙m²
Driving gear of gear pair 1	36	1.5	20	12	0.16	2×10^-4
Driven wheel of gear pair 1	90	1.5	20	12	1.3	3.04×10^-3
Driving gear of gear pair 2	29	1.5	20	12	0.09	1×10^-4
Driven wheel of gear pair 2	100	1.5	20	12	1.6	8.71×10^-3

Fig. 3The coefficient of load distribution along the meshing line

In the Fig. 3, point A and point G are the approach point and mesh-out point of double teeth-meshing, respectively. The point B and point F are the approach point and mesh-out point of single tooth-meshing, respectively. And the point D is the panel point. Based on the elastohydrodynamic lubrication model which was built, the oil film thickness and pressure distribution along the meshing position of the gear can be obtained by solving the elastohydrodynamic lubrication equation along the meshing line using the Newton method. Taking the transmission system’s gear pair 1 as the research object, when the speed of input shaft is 5500 r/min and the torque is 100 N·m, the profile of oil film pressure and oil film thickness at five key points (A, B, D, F, G) as shown in Fig. 4.

Fig. 4(a) is the oil film pressure profile of gear pair 1 along the meshing line. It shows that the second peak pressure of the meshing point F lags behind the point A, and the value is higher than the point A. The oil film pressure distribution of the double teeth-meshing approach point A is relatively stable, and its value is the smallest compared to the other four points. With the exception of point A and point G, every point are in the single tooth-meshing state and the loads are heavy which caused the oil film pressure distribution has significant fluctuations. So the oil film pressure’s value of point B, point D and point F is also higher than the rest.

Fig. 4(b) is the oil film thickness profile along the meshing line of gear pair 1. Point A is the double teeth-meshing approach point which means the process of gearing just start. So the shocking load and the entrainment velocity of point A are relatively small and therefore the oil film thickness at point A is the smallest. The load coefficients are the same of the approach point and the mesh-out point but of the oil film thickness. As the gear teeth gradually starts meshing, the curvature radius gradually increases resulting the oil film thickness gradually increasing along the meshing line because of the oil film thickness is positively correlated with the curvature radius. The amplitude fluctuation of oil film thickness at the single-teeth meshing point is also greater than that at the double-tooth meshing point, and the amplitude fluctuation of oil film thickness at the panel point D is significantly increased which is the key position of oil film fracture. It can also be seen that the maximum thickness of oil film is at the position of point G, which is due to sufficient lubricating oil supply and pile up effect at this point.

Fig. 4The coefficient of load distribution along the meshing line

a) Oil film pressure profile

b) Oil film thickness profile

Examination of Fig. 4 reveals the trend of oil film pressure and thickness by theoretical arithmetic is feasible.

3.2. The influence of speed and load on the oil film thickness of gear pairs

The variation trends of oil film thickness of gear pairs under different speeds and torques are discussed to analyze the impact of different operating conditions on the lubrication characteristics of transmission system.

The variable trend graph of the minimum and central oil film thickness along the meshing line as shown in Fig. 5.

Fig. 5The coefficient of load distribution along the meshing line

a) The state of different speeds

b) The state of different torques

Fig. 5(a) shows the oil film thickness of gear pair 1 along the meshing line under the same torque and different speeds conditions (torque: 100N·m, speed: 5500 r/min, 7000 r/min, 85000 r/min, and 10000 r/min). It can be seen that the minimum oil film thickness is smaller than the central oil film thickness, and oil film thickness not only gradually increase along the meshing line, but also gradually increase with increased speed as a whole. The oil film thickness increases sharply because of the dynamic load increased obviously with the meshing frequency to vibrate in resonance with natural frequency of transmission system when the speed is 10000 r/min.

Fig. 5(b) shows the oil film thickness of gear pair 1 along the meshing line under the same speed and different torques conditions (speed: 5500 r/min, torque:100 N·m, 130 N·m, 160 N·m, and 200 N·m). It shows that the minimum oil film thickness is smaller than the central oil film thickness along the meshing line, and oil film thickness gradually decrease with increased torque as a whole. It can also be seen that the change in torque has little effect on the oil film thickness at the teeth-meshing/tooth-meshing approach point. The principal reason for this regular is extreme pressure (EP) of oil. When the torque increased to a certain value, the distance between the oil film molecules decreases.

3.3. Simulation analysis of lubrication characteristics of transmission system

To verify the oil film thickness results obtained from the numerical analysis, this paper uses Romax designer software to analyze the lubrication performance of the wind turbine’s transmission system. And the rigid-flexible coupling model of the transmission system was built. In order to improve the efficiency of analysis, the model details were reasonably simplified while maintaining the original geometric parameters, material parameters, and spatial positional relationships.

The transmission system box is made of cast iron, with Young’s modulus of 1.25e⁵MPa, density of 7250 kg/m³ and Poisson’s ratio of 0.26. During the meshing process of the transmission system, the mesh size is 10 mm. The meshed model is shown in Fig. 6, including 420154 nodes and 271455 units. And set the bottom of the box as full constraint.

Throughput analysis, the profile of gear pair 1’s oil film thickness under the same torque and different speeds conditions (torque: 100 N·m, speed: 5500 r/min, 7000 r/min, 85000 r/min, and 10000 r/min) were obtained, as shown in Fig. 7.

It can be seen that the oil film thickness gradually increases along the meshing line. And the film thickness at the top of the teeth of the tooth-meshing approach point is a minimum, and the maximum oil film thickness at the root of the teeth of the teeth-meshing/tooth-meshing mesh-out point from Fig. 7. The thickness of the oil film gradually generally increases along the meshing line, and increases with the increase of the speed, but the increase amplitude shows a decreasing trend. Therefore, the results of simulation is consistent with the theoretical analysis.

Fig. 6The rigid-flexible coupling model of the transmission system

Fig. 7The profile of oil film thickness in the state of different speeds

a) 5500 r/min

b) 7000 r/min

c) 8500 r/min

d) 1000 r/min

From Fig. 7, it also can be seen that the minimum oil film thickness is relatively large near the right tooth surface, and the distribution along the tooth surface is not in an equal line pattern. This is due to the alternating occurrence of single and double teeth during the meshing process, the oil film thickness near the right tooth surface is relatively large. So the basic reason for this phenomenon is load deflecting which was caused by the transmission shaft inevitably undergoes bending deformation during the meshing process of the two-staged gear system.

The profile of gear pair 1’s oil film thickness under the same speed and different torques conditions (speed: 5500 r/min, torque: 100 N·m, 130 N·m, 160 N·m and 200 N·m) were obtained, as shown in Fig. 8. In broad terms, the minimum oil film thickness is relatively large near the right tooth surface and increases with torque. But increasing the amplitude of the minimum oil film thickness gradually decreased with torque increased. This is because the lubricating oil itself has extreme pressure characteristics which cause the lubricating oil molecules are difficult to continue to compress when torque reaches a certain value.

Fig. 8The profile of oil film thickness in the state of different torques

a) 100 N·m

b) 130 N·m

c) 160 N·m

d) 200 N·m

The numerical analysis and simulation analysis of the minimum oil film thickness at the approach point and the mesh-out point of the gear teeth under different working conditions are shown in Fig. 9.

Fig. 9 shows that the minimum oil film thickness increases with increasing the speed whereas decreases with increasing the torque. The results of numerical analysis and simulation analysis have the same trends at different operating conditions and show the least error. This is due to external factors such as node division accuracy and loading time of the finite element model in the simulation analysis can result in certain errors compared to theoretical analysis. However, the variation regular of oil film thickness of the gear pair 1 is consistent with the numerical analysis results under different conditions, verifying the correctness of the theoretical analysis.

Fig. 9The coefficient of load distribution along the meshing line

a) The state of different speeds

b) The state of different torques

The oil film thickness of mesh-out point is much greater than approach point whether the state of different speeds or different torques. The main reason is partial load phenomenon which eventually leads to lifetime of transmission system decreasing due to poor lubrication of the left tooth surface. Therefore, it is necessary to modify the Gear to improve this phenomenon.

4. Simulation analysis of lubrication characteristics of transmission system after modification

Based on the genetic algorithm, optimization variables and goals are regarded as evolving genome and fitness function of biological evolution to perform its flexibility evaluation. Low flexibility individual will be sifted out which rely on mechanisms, and high flexibility individual will generate the next generation which rely on gene recombination and mutation to. After enough algebra, individuals which close to the optimization goals can be obtained. The modification value for drum-shape, tooth inclination of tooth and involute can be obtained, which are 0.3621 μm, –5.96 μm, 13.28 μm, 5.44 μm, respectively. And the topology map topology of tooth surface after modification as shown in Fig. 10.

Fig. 10Topological map of tooth surface after modification

During the meshing process of the transmission system, there will be varying degrees of partial load phenomenon. This phenomenon can cause the stresses to concentrate on one side of the tooth surface. Excessive stress can cause the teeth to break and fail. Perform stress analysis on the tooth root of gear teeth after modification, and obtain the stress distribution diagram. The stress distribution of loaded tooth root before and after modification as shown in Fig. 11 and Fig. 12. It can be seen that the stress distribution at the tooth root after modification is relatively symmetrical and uniform overall. The max stress of the tooth root before and after modification is 481.0 MPa and 383.5 MPa separately when the torque is 700 N·m, which decreased by 20.7 % after modification. Through summarize the decrease of max stress in other three different operating conditions is 3.30 %, 17.81 % and 15.38 % respectively after modification. It can be seen that the modification effectively reduces the maximum stress of the tooth root.

The distribution of unit load on tooth surface directly affects the contact performance of tooth surface. The load gradient sudden change and concentration of the tooth surface will directly impact on the smoothness of gear transmission. In an ideal state, it is better that the load gradient and concentration ratio are more small.

The distribution of unit load on the tooth surface after modification under different working conditions is shown in Fig. 13, and the maximum unit load on the tooth surface of the transmission gear before and after modification is compared, as shown in Table 2. It can be seen that the load of tooth face assumes a uniform distribution after modification. The max unit load on the tooth face before and after modification is 472.0 MPa and 347.0 MPa separately when the torque is 900 N·m, which decreased by 26.48 % after modification. Through summarize the decrease of the max unit load in other three different operating conditions is 9.37 %, 21.62 % and 25.14 % respectively after modification. It can be seen that the modification effectively reduces the peak load on the tooth face.

Fig. 11Stress distribution of loaded tooth root before modification

a) 5500 r/min^-1

b) 7000 r/min^-1

c) 8500 r/min^-1

d) 1000 r/min^-1

Fig. 12Stress distribution of loaded tooth root after modification

a) 100 N·m

b) 130 N·m

c) 160 N·m

d) 200 N·m

Table 2The maximum unit load of tooth surface / N·mm-1

Operating condition	Before modification	After modification
300 N·m	182.6	176.4	3.40 %
500 N·m	334.0	274.5	17.81 %
700 N·m	481.0	383.5	20.27 %
900 N·m	624.0	528.0	15.38 %

The distribution of oil film thickness will directly affect the quality of lubrication performance. After modification, the distribution of oil film thickness in different operating conditions are shown in Fig. 14.

Fig. 13The load distribution of modified tooth surface

a) 300 N·m

b) 500 N·m

c) 700 N·m

d) 900 N·m

Fig. 14The oil film thickness after modification

a) 300 N·m

b) 500 N·m

c) 700 N·m

d) 900 N·m

From Fig. 14, it can be seen that at the top of the tooth, the oil film thickness is the highest in the entire region. The film thickness gradually decreases from the top of the tooth to the root of the tooth which shows a gradient distribution.

This is because the force at the root of the tooth is the largest, and the lubricating oil is squeezed out from the top of the tooth to the bottom of the tooth, resulting in a smaller thickness at the root of the tooth. It also can be seen that the oil film thickness is more uniform along the tooth width direction and the modification effect is better and more obvious with the increase of torque. The results indicate that modification can effectively reduce the partial load phenomenon and teeth wear to improve the lubricating performance and lifetime of gears.

5. Conclusions

1) A linear contact elastohydrodynamic lubrication model for transmission system was established. The oil film thickness and oil film pressure distribution at key points along the meshing line of the gear were obtained by using Newton iteration method. It was found that the oil film thickness gradually increases with the increase of speed and decreases with the increase of load.

2) A rigid flexible coupling model of the transmission system was established based on Romax software, and the variation regulars of oil film thickness on the tooth surface under different speed and torque conditions were analyzed, verifying the correctness of the numerical analysis. It was also found that the coupling deformation of the transmission system caused by partial load phenomenon, resulting in uneven distribution of oil film thickness along the tooth width direction.

3) Taking into account factors such as peak transmission error, load amplitude, and load distribution coefficient, genetic algorithm is used to study the micro modification of gears. Compared to before the modification, the max unit load on the tooth face of the gear transmission system has decreased by up to 26.48 %, and the maximum tooth root stress has decreased by up to 20.35 %, effectively reducing the partial load phenomenon and improving the lubricating performance and lifetime of gears.

References

H. J. Wang, “The research on the strength of planetary gear transmission with flexible pine,” (in Chinese), Chong Qing, China: Chongqing University, 2012.

Search CrossRef
K. A. Irena, H. G. Lemu, and Y. A. Kedir, “Effect of interference size on contact pressure distribution of railway wheel axle press fitting,” Designs, Vol. 7, No. 5, p. 119, Oct. 2023, https://doi.org/10.3390/designs7050119

Publisher
J. Tao, A. Wen, Z. Liu, and S. Yu, “Lean lubrication of ultra large modulus open gear and rack pair: A case study of the gear-rack drive mechanism of the Chinese “Three Gorge Dam” ship lift,” Journal of Cleaner Production, Vol. 282, p. 124450, Feb. 2021, https://doi.org/10.1016/j.jclepro.2020.124450

Publisher
P. Huang, “Numeicial Calculation Methods of Elastohydrodynamic Lubrication,” (in Chinese), Bei Jing: Tsinghua University Press, China, 2013.

Search CrossRef
H. P. Evans, R. W. Snidle, and K. J. Sharif, “Deterministic mixed lubrication modelling using roughness measurements in gear applications,” Tribology International, Vol. 42, No. 10, pp. 1406–1417, Oct. 2009, https://doi.org/10.1016/j.triboint.2009.05.025

Publisher
S. Li and A. Kahraman, “A transient mixed elastohydrodynamic lubrication model for spur gear pairs,” Journal of Tribology, Vol. 132, No. 1, p. 01150, Jan. 2010, https://doi.org/10.1115/1.4000270

Publisher
M. Liu, C. Zhu, H. Liu, H. Ding, and Z. Sun, “Effects of working conditions on TEHL performance of a helical gear pair with non-newtonian fluids,” Journal of Tribology, Vol. 136, No. 2, p. 02150, Apr. 2014, https://doi.org/10.1115/1.4026073

Publisher
M. Ankouni, A. Lubrecht, and P. Velex, “Modelling of damping in lubricated line contacts – Applications to spur gear dynamic simulations,” Proceedings of the Institution of Mechanical Engineers, Part C: Journal of Mechanical Engineering Science, Vol. 230, No. 7-8, pp. 1222–1232, Jan. 2016, https://doi.org/10.1177/0954406216628898

Publisher
H. Guan, J. X. Wang, Z. Shi, and B. Yu, “Tooth profile modification of filtering gear reducer research based on Romax Software,” (in Chinese), Machine Design and Research, Vol. 29, No. 3, pp. 28–31, 2013.

Search CrossRef
M. B. Sánchez, M. Pleguezuelos, and J. I. Pedrero, “Influence of profile modifications on meshing stiffness, load sharing, and trasnsmission error of involute spur gears,” Mechanism and Machine Theory, Vol. 139, pp. 506–525, Sep. 2019, https://doi.org/10.1016/j.mechmachtheory.2019.05.014

Publisher
J. Bruyère and P. Velex, “A simplified multi-objective analysis of optimum profile modifications in spur and helical gears,” Mechanism and Machine Theory, Vol. 80, pp. 70–83, Oct. 2014, https://doi.org/10.1016/j.mechmachtheory.2014.04.015

Publisher
J. Yang, T. Lin, Z. He, and M. Chen, “Novel calculation method for dynamic excitation of modified double-helical gear transmission,” Mechanism and Machine Theory, Vol. 167, p. 104467, Jan. 2022, https://doi.org/10.1016/j.mechmachtheory.2021.104467

Publisher
Z. L. Xiao, “Oil film model and parameters analysis for modified spur gear based on elastohydrodynamic lubrication,” (in Chinese), Chang Sha, China: Hunan University, 2017.

Search CrossRef
Y. Peng, N. Zhao, M. Zhang, W. Li, and R. Zhou, “Non-Newtonian thermal elastohydrodynamic simulation of helical gears considering modification and misalignment,” Tribology International, Vol. 124, pp. 46–60, Aug. 2018, https://doi.org/10.1016/j.triboint.2018.03.025

Publisher
L. P. Yang, “Elastohydrodynamic lubrication and its application in gear transmission,” (in Chinese), Tai Yuan, China: Taiyuan University of Technology, 2002.

Search CrossRef

About this article

Received

23 October 2023

Accepted

23 March 2024

Published

02 April 2024

SUBJECTS

System dynamics in manufacturing system modeling

DOI

https://doi.org/10.21595/jve.2024.23725

Keywords

wind turbine

transmission system

lubrication performance

and oil film

Acknowledgements

This work supported by Research Program supported by the Doctoral Program of Tianchi (tcbs201924), the Talent Youth Doctoral Program of Youth Doctoral (51052300553) and the Xinjiang University scientific research fund (BS190220).

Data Availability

The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.

Author Contributions

Li Cao: conceptualization, data curation, methodology, project administration, resources, supervision, visualization, writing-original draft preparation, writing-review and editing. Wenlei Sun: conceptualization, funding acquisition, project administration, software. Tao Gou: data curation, formal analysis, investigation, validation, visualization, writing-original draft preparation.

Conflict of interest

The authors declare that they have no conflict of interest.

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.

Previous article in issue Previous Next article in issue Next

Research article

2023 12 30

Dynamics model and vibrational response analysis of helical gear-rotor-bearing transmission system

Hongyuan Zhang, Shuo Li, Xin Zhang

Review article

2021 01 27

A review on wind turbines gearbox fault diagnosis methods

H. Gu, W. Y. Liu, Q. W. Gao, Y. Zhang

Research article

2020 02 15

Dynamic performance analysis for wind turbine in complex conditions

Li Cao, Wenlei Sun, Jianxing Zhou, Quanwei Cui

Research article

2020 02 15

Research on gearbox dynamic design methodology based on modal optimization

Zheng Li, Yu Zhou

L. Cao, W. Sun, and T. Gou, “Optimization of lubrication characteristics of wind turbine’s transmission system based on Newton Raphson method,” Journal of Vibroengineering, Vol. 26, No. 5, pp. 1200–1215, Apr. 2024, https://doi.org/10.21595/jve.2024.23725

Copy Extrica

Copied to clipboard!

TY  - JOUR
DO  - 10.21595/jve.2024.23725
UR  - https://doi.org/10.21595/jve.2024.23725
TI  - Optimization of lubrication characteristics of wind turbine’s transmission system based on Newton Raphson method
T2  - Journal of Vibroengineering
AU  - Cao, Li
AU  - Sun, Wenlei
AU  - Gou, Tao
PY  - 2024
DA  - 2024/04/02
PB  - Extrica
SP  - 1200-1215
VL  - 26
IS  - 5
SN  - 1392-8716
SN  - 2538-8460
ER  - 

Copy Ris

Copied to clipboard!

 @article{Cao_2024, title={Optimization of lubrication characteristics of wind turbine’s transmission system based on Newton Raphson method}, volume={26}, ISSN={2538-8460}, url={https://doi.org/10.21595/jve.2024.23725}, DOI={10.21595/jve.2024.23725}, number={5}, journal={Journal of Vibroengineering}, publisher={JVE International Ltd.}, author={Cao, Li and Sun, Wenlei and Gou, Tao}, year={2024}, month=apr, pages={1200–1215} }

Copy Bibtex

Copied to clipboard!

[1]L. Cao, W. Sun, and T. Gou, “Optimization of lubrication characteristics of wind turbine’s transmission system based on Newton Raphson method,” Journal of Vibroengineering, vol. 26, no. 5, pp. 1200–1215, Apr. 2024, doi: 10.21595/jve.2024.23725.

Copy IEEE

Copied to clipboard!

Cao, Li, Wenlei Sun, and Tao Gou. “Optimization of Lubrication Characteristics of Wind Turbine’s Transmission System Based on Newton Raphson Method.” Journal of Vibroengineering 26, no. 5 (April 2, 2024): 1200–1215. https://doi.org/10.21595/jve.2024.23725.

Copy Chicago

Copied to clipboard!

Optimization of lubrication characteristics of wind turbine’s transmission system based on Newton Raphson method

Abstract

Highlights

1. Introduction

2. Elastohydrodynamic lubrication model for gear line contact

2.1. Quasi static isothermal elastohydrodynamic lubrication equation

2.2. Numerical solution of elastohydrodynamic lubrication equation

3. Research on lubrication characteristics of transmission systems

3.1. Oil film thickness and pressure distribution at key meshing points of gear pairs

3.2. The influence of speed and load on the oil film thickness of gear pairs

3.3. Simulation analysis of lubrication characteristics of transmission system

4. Simulation analysis of lubrication characteristics of transmission system after modification

5. Conclusions

References

About this article

Related Articles