CN113051679B - Load processing method for main bearing of wind power gear box - Google Patents
Load processing method for main bearing of wind power gear box Download PDFInfo
- Publication number
- CN113051679B CN113051679B CN202110291371.XA CN202110291371A CN113051679B CN 113051679 B CN113051679 B CN 113051679B CN 202110291371 A CN202110291371 A CN 202110291371A CN 113051679 B CN113051679 B CN 113051679B
- Authority
- CN
- China
- Prior art keywords
- load
- file
- group
- time
- loads
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/04—Ageing analysis or optimisation against ageing
Abstract
The invention discloses a load processing method of a wind power gear box main bearing, belongs to the field of engineering calculation, and relates to a load processing method technology, wherein a main shaft bearing adopted in the current engineering calculation is calculated based on load spectrums with 6 single degrees of freedom and different rotating speeds, the calculation process is complicated, and the current mainstream fan design software does not support the load input. And the simplified load with 6 degrees of freedom is based on the load at the same rotating speed, and can be conveniently input into the design of the current mainstream commercial fan to carry out rapid fatigue life calculation.
Description
Technical Field
The invention belongs to the field of engineering calculation, relates to a load processing method technology, and particularly relates to a load processing method of a main bearing of a wind power gearbox.
Background
The main shaft bearing adopted in the current engineering calculation is calculated based on load spectrums of 6 single-degree-of-freedom different rotating speeds, the calculation process is complicated, and current mainstream fan design software does not support load input.
The invention provides a time-series-based load processing method, the processed load of a main shaft bearing is a load spectrum containing 64 loads at 6 degrees of freedom and the same rotating speed, and the load spectrum can be conveniently input into current mainstream fan design software to quickly analyze fatigue loads.
Disclosure of Invention
In view of the above, the embodiments of the present disclosure provide a load processing method for a main bearing of a wind turbine gearbox, which at least partially solves the problems in the prior art.
The embodiment of the disclosure provides a load processing method for a main bearing of a wind power gearbox, which comprises the following steps:
generating data, specifically, acquiring a load file at the center of a hub at regular sampling time intervals, and forming a load file after the sampling period is finished;
the data packets, in particular,
integrating all the load files in the time sequence into a new load file;
dividing the new load files obtained in the first step into 64 groups according to a specific arrangement sequence;
the equivalent load of each group of loads is calculated.
Further, before the data is generated, the method further comprises:
establishing data connection with a fan to obtain a load file in the fan time sequence;
verifying the load file and obtaining a verification result, wherein the verification result at least comprises Fx、Fy、Fz、Mx、My、MzTime t, rotational speed r.
Further, the integrating all the payload files in the time series into one payload file includes:
naming the payload file in the time series as n1，n2，n3……nkThe corresponding cycle times are respectively c1，c2，c3…ckAnd the payload files all contain 8 rows of data; specifically, the 1 st column to the 8 th column are sampling times t and F, respectivelyx，Fy，Fz，Mx，My，MzRotating speed r, wherein each row of data is a load;
calculating the sampling time interval Δ t;
replacing the first column of sampling time point data of each load file by delta t;
multiplying a first column of the payload files by a number of cycles for each of the payload filesThen n isiThe first column of each payload file is all Ti＝ci×Δt；
Wherein i is less than k;
and integrating the processed load files into a new load file.
Further, the specific arrangement order is, specifically, the new payload file is arranged according to Fx，Fy，Fz，Mx，My，MzThe sign plus and minus of the numerical values are arranged.
Further, the dividing the new payload files obtained in the first step into 64 groups according to a specific arrangement order includes:
marking the ith row of data from the 1 st column to the 8 th column of the new load file as ti，Fxi，Fyi，Fzi，Mxi，Myi，Mzi，ri；
According to FxThe positive and negative of the first packet are grouped for the first time; according to FyThe positive and negative of the first group are grouped for the second time; according to FzThe positive and negative of the first group are grouped for the third time; according to MxPositive and negative of (4) are grouped for the fourth time; grouping for the fifth time according to the positive and negative of My; according to MzFifth grouping is performed on the positive and negative of (1);
after five groupings are completed, the new payload file is divided into 64 groups.
Further, the calculating the equivalent load of each group of loads includes:
by Peq,i＝(∑k|Pi,k|eci,k|∑kci,k)1/eEquivalent simplification of the loads in group 64;
by passing
Determining Peq,iThe symbol of (a);
the above-mentioned
The calculated interval and the Peq,iThe intervals are the same;
the interval is a positive interval and a negative interval.
Further, before performing the equivalent simplification on the loads in the 64 groups, the method further includes:
through ci,k＝ri,k×ti,kTo ci,kCalculating;
c is mentionedi,kIs the cycle number of the kth load in the ith group of load spectra.
Further, the calculating the equivalent load of each group of loads further includes:
passing through teq,i＝ti1+ti2+ti3+…+tikFor the Peq,iThe total time under load is calculated.
Further, the calculating the equivalent load of each group of loads further includes:
through req,i＝∑kci,k/teq,iFor the Peq,iAnd calculating the equivalent rotating speed under the action of the load.
Compared with the prior art, the invention has the beneficial effects that:
the fan spindle bearing processing method provided by the invention considers the loads in the directions of 6 degrees of freedom, equivalently simplifies the original time sequence load, and separates the load directions in all degrees of freedom during simplification. And the simplified load with 6 degrees of freedom is based on the load at the same rotating speed, and can be conveniently input into the design of the current mainstream commercial fan to carry out rapid fatigue life calculation.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to the drawings without creative efforts.
FIG. 1 is a schematic block diagram of the present invention;
fig. 2 is a schematic view of payload data according to the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention more apparent, the technical solutions of the embodiments of the present invention will be described clearly and completely with reference to the accompanying drawings of the embodiments of the present invention, and it is obvious that the described embodiments are some, but not all embodiments of the present invention. All other embodiments, which can be obtained by a person skilled in the art without any inventive step based on the embodiments of the present invention, are within the scope of the present invention.
Thus, the detailed description of the embodiments of the present invention provided in the following drawings is not intended to limit the scope of the invention as claimed, but is merely representative of selected embodiments of the invention.
As shown in fig. 1-2, a method for processing load of a main bearing of a wind turbine includes:
1) data generation
As shown in fig. 2, the load file is a load file in the time sequence of the wind turbine, the load file includes 8 rows of data, and the 8 rows of data respectively correspond to Fx、Fy、Fz、Mx、My、MzLoad information of time t and rotation speed r.
When sampling, the rotating speed and the rotating speed F at the center of a hub are collected at intervals of certain sampling timex、Fy、Fz、Mx、My、Mz。
And forming a load file after the sampling period is finished.
The sampling period may cycle several times over the design life of the fan.
A complete time series load contains hundreds of load files, each containing tens of thousands of sample points.
The product of the total sampling time of each load file and the corresponding cycle number is added to obtain the design life of the fan, and the design service life of the main shaft bearing is also the design service life of the main shaft bearing, and the design life is generally 20 years.
2) Data packet
The first step is as follows: integrating all the load files in the time sequence into one load file
Naming a payload file in time series as n1，n2，n3…nkThe corresponding cycle times are respectively c1，c2，c3…ckEach payload file contains 8 columns of data, the first to 8 th columns being the sampling times t, F, respectivelyx，Fy，Fz，Mx，My，MzAnd a rotation speed r. Each row of data is a load.
(a) Calculating the sampling time interval Δ t: the sampling time interval for each payload file is the same. And selecting any load file, subtracting the sampling time point of the first row and the first column from the sampling time point of the second row and the first column to obtain a sampling time interval, and counting as delta t.
(b) The first column of sampling point-in-time data of each payload file is entirely replaced with Δ t, i.e. the duration of each payload in each payload file is Δ t.
(c) Multiplying the first column of the payload file by the number of cycles of each payload file, then the nthiThe first column of each payload file is all Ti＝ci×Δt。
(d) And integrating the processed load files into a new load file. That is, the data of the second payload file is connected to the end of the first payload file, the data of the third payload file is connected to the end of the newly formed payload file, and the operations are sequentially performed until all the payload files form a new payload file.
The second step is that: the new load file obtained in the first step is expressed as Fx、Fy、Fz、Mx、MyAnd MzThe signs of the numerical values of (a) are divided into 64 groups.
Marking the ith row of data from the first column to the eighth column of the new load file as ti、Fxi、Fyi、Fzi、Mxi、Myi、MziAnd ri。
(a) According to FxIs first grouped into positive and negative
If F in the ith rowxiIf the value of (1) is positive, the row of data is placed in the first group, if F is in the ith rowxiIs negative or equal to 0, the row of data is placed into the second group.
(b) According to FyIs divided into a second group
The first group in the second step (a) is expressed as FyIs divided into two groups, the second group of data in the second step (a) is divided into FyIs divided into two groups, and the grouping method is the same as that in the second step (a). The payload files are now divided into 4 groups.
(c) According to FzIs grouped for the third time
Each group of data in the second step (b) is expressed as FzIs divided into 2 groups, the grouping method is the same as the second step (a), and the payload file is divided into 8 groups.
(d) According to MxIs divided into a fourth sub-group
Each group of data in the second step (c) is expressed as MxIs divided into 2 groups, the grouping method is the same as the second step (a), and the payload file is divided into 16 groups.
(e) According to MyPositive and negative of (2) are grouped for the fifth time
Each group of data in the second step (d) is expressed as MyIs divided into 2 groups, the grouping method is the same as the second step (a), and the payload file is divided into 32 groups.
(f) According to MzPositive and negative of (2) are grouped for the fifth time
Each group of data in the second step (e) is expressed as MzIs divided into 2 groups, the grouping method is the same as the second step (a), and the payload file is divided into 64 groups.
The third step: calculating the equivalent load of each group of loads
Equivalent simplification of the load in each group according to equation (1)
Peq,i＝(∑k|Pi,k|eci，k/∑kci,k)1/e (1)
For a ball bearing, e is 3, for a roller bearing, e is 10/3, subscript i represents the ith group load obtained after simplification, and k represents that k loads are contained in the ith group load.
In the formula Pi,kFor loads in all directions, i.e. Fxi,k、Fyi,k、Fzi,k、Mxi,k、Myi,k、Mzi,k；
Peq,iIs equivalent load;
ci,kthe number of cycles of the kth load in the ith set of load spectra,
calculated according to the formula (2)
ci,k＝ri,k×ti,k (2)
Load P after treatment by formula (1)eq,iAre all positive values, Peq,iThe symbol of (c) can be determined as in equation (3).
If it is
If the calculated value is positive, Peq,iAlso positive, otherwise, Peq,iThe symbols need to be converted.
teq,iIs Peq,iTotal time under load, teq,iCalculating according to the formula (4):
teq,i＝ti1+ti2+ti3+…+tik (4)
req,iis Peq,iAnd (5) calculating the equivalent rotating speed under the action of the load according to the formula (5).
req,i＝∑k ci,k/teq,i (5)
After the above steps, the original time-series load is simplified to 64 inclusion times t, Fx，Fy，Fz，Mx，My，MzLoad spectrum of the rotational speed r.
The fan spindle bearing processing method provided by the invention considers the loads in the directions of 6 degrees of freedom, equivalently simplifies the original time sequence load, and separates the load directions in all degrees of freedom during simplification. And the simplified load with 6 degrees of freedom is based on the load at the same rotating speed, and can be conveniently input into the design of the current mainstream commercial fan to carry out rapid fatigue life calculation.
An embodiment of the present disclosure further provides an electronic device, including:
at least one processor; and a memory communicatively coupled to the at least one processor; wherein the content of the first and second substances,
the memory stores instructions executable by the at least one processor to enable the at least one processor to perform the method of load handling in the method embodiments described above.
The above formulas are all calculated by taking the numerical value of the dimension-removed data, the formula is one closest to the real situation obtained by collecting a large amount of data and performing software simulation, and the preset parameters in the formula are set by the technical personnel in the field according to the actual situation.
In the embodiments provided by the present invention, it should be understood that the disclosed apparatus, device and method can be implemented in other ways. For example, the above-described apparatus embodiments are merely illustrative, and for example, the division of the modules is only one logical functional division, and there may be other divisions when the actual implementation is performed; the modules described as separate parts may or may not be physically separate, and parts displayed as modules may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the method of the embodiment.
It will also be evident to those skilled in the art that the invention is not limited to the details of the foregoing illustrative embodiments, and that the present invention may be embodied in other specific forms without departing from the spirit or essential attributes thereof.
The present embodiments are therefore to be considered in all respects as illustrative and not restrictive, the scope of the invention being indicated by the appended claims rather than by the foregoing description, and all changes which come within the meaning and range of equivalency of the claims are therefore intended to be embraced therein. Any reference signs in the claims shall not be construed as limiting the claim concerned.
Furthermore, it is obvious that the word "comprising" does not exclude other elements or steps, and the singular does not exclude the plural. A plurality of units or means recited in the system claims may also be implemented by one unit or means in software or hardware. The terms second, etc. are used to denote names, but not any particular order.
Finally, it should be noted that the above examples are only intended to illustrate the technical process of the present invention and not to limit the same, and although the present invention has been described in detail with reference to the preferred embodiments, it will be understood by those skilled in the art that modifications or equivalent substitutions may be made to the technical process of the present invention without departing from the spirit and scope of the technical process of the present invention.
Claims (5)
1. A method of load handling of a main bearing of a wind turbine gearbox, the method comprising:
generating data, specifically, acquiring a load file at the center of a hub at regular sampling time intervals, and forming a load file after the sampling period is finished;
the data packets, in particular,
integrating all the load files in the time sequence into a new load file;
dividing the new load files obtained in the first step into 64 groups according to a specific arrangement sequence;
calculating the equivalent load of each group of loads;
the calculating the equivalent load of each group of loads comprises the following steps:
by passing
Equivalent simplification of the loads in group 64;
for a ball bearing, e is 3, for a roller bearing, e is 10/3, subscript i represents the i-th group load obtained after simplification, k represents that the i-th group load contains k loads, wherein Pi and k are loads in all directions, namely, Fxi, k, Fyi, k, Fzi, k, Mxi, k, Myi, k, Mzi, k, Peq, i are equivalent loads, and ci and k are the cycle number of the k-th load in the i-th group load spectrum;
by passing
Determining Peq,iThe symbol of (a);
the above-mentioned
The calculated interval and the Peq,iThe intervals are the same;
the interval is a positive interval and a negative interval;
before the equivalent simplification of the loads in the 64 groups, the method further comprises the following steps:
by passing
To ci,kCalculating;
c is mentionedi,kThe cycle number of the kth load in the ith group of load spectrums;
the calculating the equivalent load of each group of loads further comprises:
by passing
To the above
Calculating the total time under the action of the load;
the calculating the equivalent load of each group of loads further comprises:
by passing
To the above
And calculating the equivalent rotating speed under the action of the load.
2. The method for processing the load of the main bearing of the wind turbine gearbox according to claim 1, wherein before the data is generated, the method further comprises:
establishing data connection with a fan to obtain a load file in the fan time sequence;
verifying the load file and obtaining a verification result, wherein the verification result at least comprises Fx、Fy、Fz、Mx、My、MzSampling time t and rotating speed r.
3. The method for processing the load of the main bearing of the wind power gearbox according to claim 1, wherein the step of integrating all the load files in the time series into one load file comprises the following steps:
naming the payload file in the time series as n1，n2，n3……nkThe corresponding cycle times are respectively c1，c2，c3…ckAnd the payload files all contain 8 rows of data; specifically, the 1 st column to the 8 th column are sampling times t and F, respectivelyx，Fy，Fz，Mx，My，MzRotating speed r, wherein each row of data is a load;
calculating a sampling time interval delta t;
replacing the first column of sampling time point data of each load file by delta t;
multiplying the first column of the payload file by the number of cycles of each of the payload files, then the nthiFirst column of individual payload filesAll are ti=ci×Δt；
Wherein i is less than k;
and integrating the processed load files into a new load file.
4. The method for processing the load of the main bearing of the wind power gearbox according to claim 1, wherein the specific arrangement sequence is that the new load file is Fx，Fy，Fz，Mx，My，MzThe sign plus and minus of the numerical values are arranged.
5. The method for processing the load of the main bearing of the wind power gearbox according to claim 1, wherein the new load files obtained in the first step are divided into 64 groups according to a specific arrangement sequence, and the method comprises the following steps:
marking the ith row of data from the 1 st column to the 8 th column of the new load file as ti，Fxi，Fyi，Fzi，Mxi，Myi，Mzi，ri；
Grouping for the first time according to the positive and negative of Fx; according to FyThe positive and negative of the first group are grouped for the second time; according to FzThe positive and negative of the first group are grouped for the third time; according to MxPositive and negative of (4) are grouped for the fourth time; according to MyFifth grouping is performed on the positive and negative of (1); according to MzThe positive and negative of (2) are grouped for the sixth time;
after six groupings are completed, the new payload file is divided into 64 groups.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110291371.XA CN113051679B (en) | 2021-03-18 | 2021-03-18 | Load processing method for main bearing of wind power gear box |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110291371.XA CN113051679B (en) | 2021-03-18 | 2021-03-18 | Load processing method for main bearing of wind power gear box |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113051679A CN113051679A (en) | 2021-06-29 |
| CN113051679B true CN113051679B (en) | 2021-10-26 |
Family
ID=76513692
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110291371.XA Active CN113051679B (en) | 2021-03-18 | 2021-03-18 | Load processing method for main bearing of wind power gear box |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113051679B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117592221A (en) * | 2024-01-18 | 2024-02-23 | 洛阳Lyc轴承有限公司 | Load processing method for checking main bearing of fan |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105701337A (en) * | 2015-12-31 | 2016-06-22 | 北京金风科创风电设备有限公司 | Fatigue life prediction method and device for wind turbine generator |
| CN106600066A (en) * | 2016-12-19 | 2017-04-26 | 南京理工大学 | SCADA data-based wind driven generator gearbox fatigue life estimation method |
| CN112016170A (en) * | 2020-08-27 | 2020-12-01 | 国电联合动力技术有限公司 | Wind power main shaft bearing design calculation and model selection method and device |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102385377B (en) * | 2011-09-16 | 2014-05-07 | 上海卡鲁自动化科技有限公司 | Hardware in the loop experiment system used for wind power generator master control system test and method thereof |
| GB201121590D0 (en) * | 2011-12-15 | 2012-01-25 | Lm Wind Power As | A wind turbine blade control method |
| CN109726414A (en) * | 2017-10-30 | 2019-05-07 | 北京万源工业有限公司 | A kind of hub of wind power generator structural fatigue calculation method and analysis of fatigue device |
| CN112883555B (en) * | 2021-01-25 | 2024-03-29 | 运达能源科技集团股份有限公司 | Multi-directional load extraction method for each section of blade of wind generating set |
-
2021
- 2021-03-18 CN CN202110291371.XA patent/CN113051679B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105701337A (en) * | 2015-12-31 | 2016-06-22 | 北京金风科创风电设备有限公司 | Fatigue life prediction method and device for wind turbine generator |
| CN106600066A (en) * | 2016-12-19 | 2017-04-26 | 南京理工大学 | SCADA data-based wind driven generator gearbox fatigue life estimation method |
| CN112016170A (en) * | 2020-08-27 | 2020-12-01 | 国电联合动力技术有限公司 | Wind power main shaft bearing design calculation and model selection method and device |
Non-Patent Citations (2)
| Title |
|---|
| 计算参数对数值修正载荷识别算法的影响;徐菁等;《南京航空航天大学学报》;20160215(第01期);全文 * |
| 起重船海浪扰动及补偿系统研究;王盼;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20180315(第03期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113051679A (en) | 2021-06-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110417011B (en) | Online dynamic security assessment method based on mutual information and iterative random forest | |
| US20160169205A1 (en) | Method for constructing wind power connection system model based on measured data | |
| WO2023063888A2 (en) | Method and apparatus for predicting wind power, and device and storage medium thereof | |
| CN113051679B (en) | Load processing method for main bearing of wind power gear box | |
| CN107483397A (en) | A kind of examinee's self-test assessment method based on B/S framework | |
| CN108052963A (en) | The data screening method, apparatus and wind power generating set of wind power prediction modeling | |
| CN112882899B (en) | Log abnormality detection method and device | |
| CN104574221B (en) | A kind of photovoltaic plant running status discrimination method based on loss electricity characteristic parameter | |
| CN114046903A (en) | Real-time temperature prediction method and system for high-voltage cable core | |
| CN107356431A (en) | Rolling bearing performance degradation assessment method based on ADMM Yu sparse ensemble learning | |
| CN107527093A (en) | A kind of running of wind generating set method for diagnosing status and device | |
| CN116611741A (en) | Construction method and system of service quality index system based on wind power equipment | |
| CN116542030A (en) | Double-fed fan parameter identification method, system and equipment based on track sensitivity | |
| CN115600494A (en) | Low-voltage distribution area topology automatic identification method and device | |
| CN115858613A (en) | Wind power data display method and system based on time sequence information | |
| CN112365039B (en) | Method and system for predicting monthly electricity consumption based on longitudinal and transverse correlation | |
| US3129518A (en) | Method of and apparatus for teaching mathematics | |
| CN105303056B (en) | A kind of wind speed real-time change speed depicting method | |
| CN110796243B (en) | Continuous operation monitoring data simulation generation method and device | |
| CN111950131A (en) | Wind power plant output equivalent aggregation model construction method considering electricity limiting factors | |
| CN113887839A (en) | Clustering-based wind power plant ultra-short-term wind power prediction method and system | |
| CN115392101A (en) | Method and system for generating comprehensive energy random scene | |
| CN112380748A (en) | Wind-induced random vibration analysis method and related device for power transmission tower | |
| CN107133441B (en) | Method and device for determining power data sampling precision in energy Internet | |
| CN105678463B (en) | Method for simulating abnormal data of electric energy metering device for training |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |