Edge contact stress three-row cylindrical roller bearing analytical calculation method and device

Abstract

The invention relates to the technical field of wind power generation, and discloses an edge contact stress three-row cylindrical roller bearing analytical calculation method and device, and the method comprises the steps: obtaining design calculation basic data of a three-row cylindrical roller bearing; the three rows of cylindrical roller bearings are sliced according to the design calculation basic data, and k bearing slices are obtained; performing deformation calculation on the k bearing slices to obtain the deformation of the k bearing slices; performing contact force calculation by adopting the deformation amount of the k bearing slices to obtain the real contact stress of the k bearing slices; and carrying out operation life evaluation on the three rows of cylindrical roller bearings by adopting the real contact stresses corresponding to the k bearing slices respectively to obtain the operation life of the three rows of cylindrical roller bearings. According to the edge contact stress three-row cylindrical roller bearing analytical calculation method and device, the problem that edge contact stress cannot be considered in a three-row cylindrical roller bearing engineering algorithm is solved.

Description

Technical Field

[0001] The invention relates to the technical field of wind power generation, and in particular to an analytical calculation method and a device for edge contact stress of a three-row cylindrical roller bearing. Background Art

[0002] Driven by the continuous development of the wind power market, the research and development of wind turbines by domestic main engine manufacturers is developing in the direction of large capacity and high conversion efficiency. As the capacity of wind turbines increases, three-row cylindrical roller bearings have become the main form of variable pitch bearings. Their service life and reliability directly determine whether the entire wind turbine can operate stably during its service life.

[0003] The variable pitch bearing is an important connecting component between the wind turbine blades and the hub. It plays a key role in connecting the upper and lower parts and the front and back of the wind turbine. It bears the wind load transmitted by the blades and the weight of the blades themselves, and transmits these loads to the hub. As a key component of the wind turbine, the variable pitch bearing needs to meet the design requirements related to the extreme working condition strength and fatigue cumulative damage. The design cycle is long and difficult.

[0004] Traditional three-row roller bearing calculations mainly include engineering calculation methods and finite element calculation methods. At present, the engineering calculation method is mainly based on the rigid ring assumption, which cannot consider the impact of edge contact stress and roller modification on raceway contact stress and bearing life; although the finite element calculation method has high accuracy, it has high requirements on mesh quality, computer performance and engineer calculation experience, and the solution is time-consuming.

[0005] Therefore, it is necessary to optimize the engineering calculation method to achieve the effect of significantly improving the accuracy while not increasing the solution time much. Summary of the invention

[0006] In view of the deficiencies in the prior art, the present invention proposes a method and device for analytically calculating edge contact stress of a three-row cylindrical roller bearing, which are used to solve the above technical problems.

[0007] First, an analytical calculation method for edge contact stress of three-row cylindrical roller bearings is provided, including:

[0008] Obtain the basic design calculation data for three-row cylindrical roller bearings;

[0009] Slicing the three-row cylindrical roller bearing according to the design calculation basic data to obtain k bearing slices;

[0010] Calculating the deformation of the k bearing slices respectively to obtain the deformation of the k bearing slices;

[0011] The contact force is calculated by using the deformation of the k bearing slices to obtain the real contact stress of the k bearing slices;

[0012] The operating life of the three-row cylindrical roller bearing is evaluated using the true contact stress corresponding to k bearing slices, and the operating life of the three-row cylindrical roller bearing is obtained.

[0013] Furthermore, the deformation amounts of the k bearing slices are calculated respectively to obtain the deformation amounts of the k bearing slices, including:

[0014] Decompose the axial displacement, radial displacement and rotation angle of the bearing into k bearing slices respectively, and obtain k decomposed bearing slices;

[0015] Correction processing is performed on the k decomposed bearing slices to obtain deformation amounts of the k bearing slices.

[0016] Furthermore, the contact force is calculated based on the deformation of the k bearing slices to obtain the real contact stress of the k bearing slices, including:

[0017] The equilibrium equation is constructed by using the deformation amounts of k bearing slices to obtain the equilibrium equation;

[0018] Iteratively processing the equilibrium equation to obtain contact stresses of k bearing slices;

[0019] The contact stress of the slices corresponding to the m rollers in the three-row cylindrical roller bearing is extracted from the contact stress of the k bearing slices, and the summation is performed to obtain the total contact stress corresponding to the m rollers;

[0020] Construct an influence coefficient matrix according to the total contact stress corresponding to each of the m rollers;

[0021] The contact stress of each bearing slice is updated according to the influence coefficient matrix to obtain the real contact stresses corresponding to the k bearing slices.

[0022] Furthermore, the deformation of k bearing slices is used to construct a balance equation to obtain a balance equation, including:

[0023] According to the rolling element load deformation formula, the deformation of the k bearing slices is respectively calculated to obtain the load of the k bearing slices;

[0024] Decomposing the loads of the k bearing slices respectively to obtain components in n directions;

[0025] The components in the n directions are combined to obtain the equilibrium equation.

[0026] Furthermore, the iterative processing of the equilibrium equation to obtain the contact stress of k bearing slices includes:

[0027] The equilibrium equation is iteratively processed according to the Newton-Raphson method to obtain the roller bearing displacement result;

[0028] If the roller bearing displacement result converges, the roller bearing displacement result is calculated according to the Hertz formula to obtain the contact stress of k bearing slices;

[0029] If the roller bearing displacement results do not converge, adjust the number of slices and iteration tolerance and iterate again.

[0030] Furthermore, the influence coefficient matrix is ​​constructed according to the total contact stresses corresponding to the m rollers, including:

[0031] The total contact stress of the m rollers is calculated separately according to the finite length contact algorithm, and the influence coefficient of the slices in the m rollers is extracted;

[0032] The influence coefficients of the slices in the m rollers are combined to obtain an influence coefficient matrix.

[0033] Furthermore, the three-row cylindrical roller bearing includes three raceways, and the actual contact stress of each slice is calculated to obtain the operating life of the three-row cylindrical roller bearing, including:

[0034] Calculating the actual contact stress of each slice according to the life prediction model to obtain the service life of each raceway;

[0035] The service life of all raceways is combined to obtain the operating life of the three-row cylindrical roller bearing.

[0036] In a second aspect, a device for analyzing and calculating edge contact stress of a three-row cylindrical roller bearing is provided, which is characterized by comprising:

[0037] An acquisition module is used to acquire basic design calculation data of three-row cylindrical roller bearings;

[0038] A slicing processing module, used for slicing the three-row cylindrical roller bearing according to the design calculation basic data to obtain k bearing slices;

[0039] A first calculation module is used to calculate the deformation of the k bearing slices respectively to obtain the deformation of the k bearing slices;

[0040] A second calculation module is used to calculate the contact force using the deformation of the k bearing slices to obtain the real contact stress of the k bearing slices;

[0041] The evaluation module is used to evaluate the operating life of the three-row cylindrical roller bearing by using the real contact stresses corresponding to the k bearing slices, so as to obtain the operating life of the three-row cylindrical roller bearing.

[0042] In a third aspect, a device for analytical calculation of edge contact stress of three-row cylindrical roller bearings is provided, comprising a processor and a memory storing program instructions, wherein the processor is configured to execute the analytical calculation method for edge contact stress of three-row cylindrical roller bearings as described above when running the program instructions.

[0043] In a fourth aspect, an electronic device is provided, characterized in that it includes the edge contact stress analytical calculation device for three-row cylindrical roller bearings as described above.

[0044] The invention adopting the above technical solution has the following advantages:

[0045] 1. The present invention introduces a finite length contact algorithm and considers the roller edge contact stress by constructing a slice influence coefficient matrix, thereby solving the problem that the engineering algorithm of three-row cylindrical roller bearings cannot consider the edge contact stress, and can calculate the bearing edge contact stress efficiently and more accurately.

[0046] 2. Compared with finite element calculation, the analytical method of the present invention greatly reduces the time required for calculation and does not rely on the calculation experience of engineers. BRIEF DESCRIPTION OF THE DRAWINGS

[0047] In order to more clearly illustrate the specific implementation of the present invention, the following will briefly introduce the drawings required for use in the specific implementation. In all the drawings, each element or part is not necessarily drawn according to the actual scale.

[0048] Figure 1 Flow chart of the analytical calculation method of edge contact stress of three-row cylindrical roller bearings of the present invention (I);

[0049] Figure 2 Flow chart (II) of the analytical calculation method for edge contact stress of three-row cylindrical roller bearings of the present invention;

[0050] Figure 3 It is a schematic diagram of a slice model in the analytical calculation method of edge contact stress of a three-row cylindrical roller bearing of the present invention;

[0051] Figure 4 It is a force diagram of a three-row roller bearing in the analytical calculation method of edge contact stress of a three-row cylindrical roller bearing of the present invention;

[0052] Figure 5 A schematic diagram of the influence coefficient in the analytical calculation method of edge contact stress of a three-row cylindrical roller bearing of the present invention;

[0053] Figure 6 A schematic diagram of the contact force of the bearing rollers in the analytical calculation method of the edge contact stress of the three-row cylindrical roller bearing of the present invention;

[0054] Figure 7 This is a flow chart of the analytical calculation device for edge contact stress of three-row cylindrical roller bearings of the present invention. DETAILED DESCRIPTION

[0055] The following embodiments of the technical solution of the present invention are described in detail in conjunction with the accompanying drawings. The following embodiments are only used to more clearly illustrate the technical solution of the present invention, and are therefore only used as examples, and cannot be used to limit the protection scope of the present invention.

[0056] It should be noted that, unless otherwise specified, the technical terms or scientific terms used in this application should be the common meanings understood by technicians in the field to which the present invention belongs. The terms "first", "second", etc. in the specification and claims of the embodiments of the present disclosure and the above-mentioned drawings are used to distinguish similar objects, and are not necessarily used to describe a specific order or sequence. It should be understood that the data used in this way can be interchanged where appropriate, so as to implement the embodiments of the present disclosure described here. In addition, the terms "including" and "having" and any of their variations are intended to cover non-exclusive inclusions. Unless otherwise specified, the term "multiple" means two or more. In the embodiments of the present disclosure, the character " / " indicates that the objects before and after are in an "or" relationship. For example, A / B means: A or B. The term "and / or" is a description of the association relationship of objects, indicating that there can be three relationships. For example, A and / or B means: A or B, or, A and B. The term "corresponding" can refer to an association relationship or a binding relationship, and A and B correspondingly refer to an association relationship or a binding relationship between A and B.

[0057] like Figure 1 to Figure 7 As shown, the present invention provides an analytical calculation method for edge contact stress of a three-row cylindrical roller bearing, comprising:

[0058] Step S01, obtaining basic design calculation data of a three-row cylindrical roller bearing, wherein the basic design calculation data includes a geometric model of the three-row cylindrical roller bearing, roller modification parameters, radial and axial clearances of the bearing, and limit and fatigue loads;

[0059] Step S02, slicing the three-row cylindrical roller bearing according to the design calculation basic data to obtain k bearing slices;

[0060] Step S03, calculating the deformation amounts of the k bearing slices respectively to obtain the deformation amounts of the k bearing slices;

[0061] Step S04, using the deformation of the k bearing slices to calculate the contact force, to obtain the real contact stress of the k bearing slices;

[0062] Step S05, using the real contact stresses corresponding to the k bearing slices to evaluate the service life of the three-row cylindrical roller bearing, to obtain the service life of the three-row cylindrical roller bearing;

[0063] Specifically, the geometric model, roller modification parameters, radial and axial clearances, limit and fatigue loads of the three-row cylindrical roller bearing used for design calculation are obtained;

[0064] The bearing rollers are sliced ​​according to the calculation accuracy, and the roller shaping amount at the center of each slice is calculated. The roller slicing model is shown as follows: Figure 3 shown.

[0065] By introducing the finite length contact algorithm and constructing the slice influence coefficient matrix to consider the roller edge contact stress, the problem that the engineering algorithm of three-row cylindrical roller bearings cannot consider the edge contact stress is solved, and the bearing edge contact stress can be calculated efficiently and more accurately.

[0066] In some embodiments, calculating the deformation amounts of k bearing slices respectively to obtain the deformation amounts of the k bearing slices includes:

[0067] Decompose the axial displacement, radial displacement and rotation angle of the bearing into k bearing slices respectively, and obtain k decomposed bearing slices;

[0068] The k decomposed bearing slices are corrected to obtain the deformation amounts of the k bearing slices.

[0069] Specifically, according to the number of slices, the axial displacement, radial displacement, and rotation angle of the three-row cylindrical roller bearing are decomposed into each slice respectively. The bearing displacement is as follows: Figure 4 As shown; the correction of the slice deformation by axial clearance, radial clearance and roller modification is considered at the same time.

[0070] In some embodiments, the contact force is calculated based on the deformation of k bearing slices to obtain the real contact stress of the k bearing slices, including:

[0071] The equilibrium equation is constructed by using the deformation of k bearing slices to obtain the equilibrium equation;

[0072] The equilibrium equation is iterated to obtain the contact stress of k bearing slices;

[0073] The contact stress of the slices corresponding to the m rollers in the three-row cylindrical roller bearing is extracted from the contact stress of the k bearing slices, and the summation is performed to obtain the total contact stress corresponding to the m rollers;

[0074] Construct an influence coefficient matrix according to the total contact stress corresponding to each of the m rollers;

[0075] The contact stress of each bearing slice is updated according to the influence coefficient matrix to obtain the true contact stress corresponding to k bearing slices.

[0076] In some embodiments, the deformation of k bearing slices is used to construct a balance equation to obtain a balance equation, including:

[0077] According to the rolling element load deformation formula, the deformation of k bearing slices is calculated respectively to obtain the load of k bearing slices;

[0078] The loads of k bearing slices are decomposed separately to obtain components in n directions;

[0079] Combining the components in n directions gives the equilibrium equation.

[0080] Specifically, according to the deformation of each roller slice, the load on each roller slice is calculated using the rolling element load deformation formula, the load on the roller slice is decomposed along the axial, radial and rotational directions of the bearing, the resultant force in the three directions is calculated, and the equilibrium equation is listed.

[0081] Among them, the rolling element load deformation formula is expressed as:

[0082]

[0083] Where: q k is the load on the kth slice of the cylindrical roller; L we is the effective contact length of the cylindrical roller; n s is the number of cylindrical roller slices, which is a cardinality greater than 30; δ k is the elastic deformation of the kth slice of the cylindrical roller.

[0084] In some embodiments, the equilibrium equation is iteratively processed to obtain the contact stress of k bearing slices, including:

[0085] The equilibrium equation is iterated according to the Newton-Raphson method to obtain the displacement result of the roller bearing;

[0086] If the roller bearing displacement results converge, the roller bearing displacement results are calculated according to the Hertz formula to obtain the contact stress of k bearing slices;

[0087] If the roller bearing displacement results do not converge, adjust the number of slices and iteration tolerance and iterate again.

[0088] Specifically, according to the equilibrium equation, the Newton-Raphson method is used for iterative solution. By controlling the reasonable tolerance and the number of roller slices, the bearing displacement result that meets the calculation accuracy is solved. If the result does not converge, the number of slices and the iteration tolerance can be modified to complete the iterative solution. The contact stress on each slice is obtained by solving the Hertz formula.

[0089] Among them, the Hertz formula is expressed as:

[0090]

[0091] Where: σ is the contact stress within the contact range; a is the major semi-axis of the contact area; b is the minor semi-axis of the contact area; Q is the contact force; x and y are the coordinates within the contact area, respectively.

[0092] In some embodiments, an influence coefficient matrix is ​​constructed according to the total contact stresses corresponding to the m rollers, including:

[0093] The total contact stress of the m rollers is calculated separately according to the finite length contact algorithm, and the influence coefficient of the slices in the m rollers is extracted;

[0094] The influence coefficients of the slices in the m rollers are combined to obtain an influence coefficient matrix.

[0095] Specifically, the contact stress of each slice is calculated and summed up to obtain the contact force on each roller;

[0096] According to the contact force on each roller, the influence coefficient between each slice is obtained using the finite length contact algorithm. Figure 5 As shown, the combination becomes the influence coefficient matrix.

[0097] Among them, the finite length contact algorithm formula is expressed as:

[0098]

[0099] Where: Ω is the contact area; Q is the contact force between the roller and the raceway; p(x,y) is the contact stress at the coordinate (x,y); z(x,y) is the initial spacing of the coordinate (x,y); δ is the deformation of the roller and the raceway; E' is the equivalent elastic modulus; υ1 and υ2 are the Poisson's ratios of the roller and raceway materials, respectively; E1 and E2 are the elastic moduli of the roller and raceway materials, respectively; x' and y' are the coordinate values ​​of the reference unit.

[0100] In some embodiments, the three-row cylindrical roller bearing includes three raceways, and the actual contact stress of each slice is calculated to obtain the operating life of the three-row cylindrical roller bearing, including:

[0101] The actual contact stress of each slice is calculated according to the life prediction model to obtain the service life of each raceway;

[0102] The service life of all raceways is combined to obtain the operating life of the three-row cylindrical roller bearing.

[0103] Specifically, according to the influence coefficient matrix and the contact stress on each slice, the stress distribution of each slice of the roller is iteratively solved to obtain the contact stress of each slice. The comparison diagram of whether to consider stress redistribution is shown in the figure below. Figure 6 As shown;

[0104] Based on the calculated contact stress, the service life of each raceway of the bearing is calculated, and the service life of the three-row cylindrical roller bearing is calculated by combining the service lives.

[0105] In some other embodiments, an analytical calculation device for edge contact stress of a three-row cylindrical roller bearing is provided, characterized in that it includes:

[0106] An acquisition module is used to acquire basic design calculation data of three-row cylindrical roller bearings;

[0107] A slice processing module is used to slice the three-row cylindrical roller bearing according to the design calculation basic data to obtain k bearing slices;

[0108] The first calculation module is used to calculate the deformation of k bearing slices respectively to obtain the deformation of the k bearing slices;

[0109] The second calculation module is used to calculate the contact force using the deformation of the k bearing slices to obtain the real contact stress of the k bearing slices;

[0110] The evaluation module is used to evaluate the operating life of the three-row cylindrical roller bearing by using the real contact stresses corresponding to the k bearing slices, so as to obtain the operating life of the three-row cylindrical roller bearing.

[0111] In other embodiments, an edge contact stress analytical calculation device for three-row cylindrical roller bearings is provided, comprising a processor and a memory storing program instructions, wherein the processor is configured to execute the edge contact stress analytical calculation method for three-row cylindrical roller bearings as described above when running the program instructions.

[0112] In other embodiments, an electronic device (eg, a mobile phone, a computer) is provided, comprising the edge contact stress analytical calculation device for a three-row cylindrical roller bearing as described above.

[0113] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit it. Although the present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or replace some or all of the technical features therein by equivalents. These modifications or replacements do not deviate the essence of the corresponding technical solutions from the scope of the technical solutions of the embodiments of the present invention, and they should all be included in the scope of the claims and specification of the present invention.

Claims

1. Analytical calculation method of edge contact stress of three-row cylindrical roller bearings, characterized in that: include: Obtain the basic design calculation data for three-row cylindrical roller bearings; Slicing the three-row cylindrical roller bearing according to the design calculation basic data to obtain k bearing slices; Calculating the deformation of the k bearing slices respectively to obtain the deformation of the k bearing slices; The contact force is calculated by using the deformation of the k bearing slices to obtain the real contact stress of the k bearing slices; The operating life of the three-row cylindrical roller bearing is evaluated using the true contact stress corresponding to k bearing slices, and the operating life of the three-row cylindrical roller bearing is obtained.

2. The edge contact stress analytical calculation method of three-row cylindrical roller bearings according to claim 1 is characterized in that: The step of calculating the deformation amounts of the k bearing slices respectively to obtain the deformation amounts of the k bearing slices comprises: Decompose the axial displacement, radial displacement and rotation angle of the bearing into k bearing slices respectively, and obtain k decomposed bearing slices; Correction processing is performed on the k decomposed bearing slices to obtain deformation amounts of the k bearing slices.

3. The edge contact stress analytical calculation method of three-row cylindrical roller bearings according to claim 1 is characterized in that: The contact force calculation is performed based on the deformation of the k bearing slices to obtain the real contact stress of the k bearing slices, including: The equilibrium equation is constructed by using the deformation amounts of k bearing slices to obtain the equilibrium equation; Iteratively processing the equilibrium equation to obtain contact stresses of k bearing slices; The contact stress of the slices corresponding to the m rollers in the three-row cylindrical roller bearing is extracted from the contact stress of the k bearing slices, and the summation is performed to obtain the total contact stress corresponding to the m rollers; Construct an influence coefficient matrix according to the total contact stress corresponding to each of the m rollers; The contact stress of each bearing slice is updated according to the influence coefficient matrix to obtain the real contact stresses corresponding to the k bearing slices.

4. The analytical calculation method for edge contact stress of three-row cylindrical roller bearings according to claim 3 is characterized in that: The equilibrium equation is constructed by using the deformation amounts of k bearing slices to obtain the equilibrium equation, including: According to the rolling element load deformation formula, the deformation of the k bearing slices is respectively calculated to obtain the load of the k bearing slices; Decomposing the loads of the k bearing slices respectively to obtain components in n directions; The components in the n directions are combined to obtain the equilibrium equation.

5. The edge contact stress analytical calculation method of three-row cylindrical roller bearings according to claim 3 is characterized in that: The iterative processing of the equilibrium equation to obtain the contact stress of k bearing slices includes: The equilibrium equation is iteratively processed according to the Newton-Raphson method to obtain the roller bearing displacement result; If the roller bearing displacement result converges, the roller bearing displacement result is calculated according to the Hertz formula to obtain the contact stress of k bearing slices; If the roller bearing displacement results do not converge, adjust the number of slices and iteration tolerance and iterate again.

6. The edge contact stress analytical calculation method of three-row cylindrical roller bearings according to claim 3 is characterized in that: The influence coefficient matrix is ​​constructed according to the total contact stresses corresponding to the m rollers, including: The total contact stress of the m rollers is calculated separately according to the finite length contact algorithm, and the influence coefficient of the slices in the m rollers is extracted; The influence coefficients of the slices in the m rollers are combined to obtain an influence coefficient matrix.

7. The edge contact stress analytical calculation method of three-row cylindrical roller bearings according to claim 1 is characterized in that: The three-row cylindrical roller bearing includes three raceways, and the actual contact stress of each slice is calculated to obtain the operating life of the three-row cylindrical roller bearing, including: Calculating the actual contact stress of each slice according to the life prediction model to obtain the service life of each raceway; The service life of all raceways is combined to obtain the operating life of the three-row cylindrical roller bearing.

8. An analytical calculation device for edge contact stress of three-row cylindrical roller bearings, characterized in that it includes: An acquisition module is used to acquire basic design calculation data of three-row cylindrical roller bearings; A slicing processing module, used for slicing the three-row cylindrical roller bearing according to the design calculation basic data to obtain k bearing slices; A first calculation module is used to calculate the deformation of the k bearing slices respectively to obtain the deformation of the k bearing slices; A second calculation module is used to calculate the contact force using the deformation of the k bearing slices to obtain the real contact stress of the k bearing slices; The evaluation module is used to evaluate the operating life of the three-row cylindrical roller bearing by using the real contact stresses corresponding to the k bearing slices, so as to obtain the operating life of the three-row cylindrical roller bearing.

9. An analytical calculation device for edge contact stress of three-row cylindrical roller bearings, comprising a processor and a memory storing program instructions, wherein the processor is configured to execute the analytical calculation method for edge contact stress of three-row cylindrical roller bearings as described in any one of claims 1 to 7 when running the program instructions.

10. An electronic device, characterized in that: It comprises the edge contact stress analytical calculation device of three-row cylindrical roller bearings as claimed in claim 8 or 9.

Images