CN104573172A - Fatigue analysis method and fatigue analysis device of structural member in wind generating set - Google Patents
Fatigue analysis method and fatigue analysis device of structural member in wind generating set Download PDFInfo
- Publication number
- CN104573172A CN104573172A CN201410664081.5A CN201410664081A CN104573172A CN 104573172 A CN104573172 A CN 104573172A CN 201410664081 A CN201410664081 A CN 201410664081A CN 104573172 A CN104573172 A CN 104573172A
- Authority
- CN
- China
- Prior art keywords
- load
- structural member
- fatigue
- finite element
- stress
- 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.)
- Granted
Links
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Wind Motors (AREA)
Abstract
The invention provides a fatigue analysis method and a fatigue analysis device of a structural member in a wind generating set. The method comprises the following steps: establishing a finite element model for a to-be-measured structural member involving bearing connection and related structural members thereof; applying load to the finite element model according to various preset load working conditions respectively; acquiring each node stress result of the to-be-measured structural member under each load working condition through finite element analysis; scaling down the acquired stress results by taking a fatigue load sequence which is pretreated through regularization as a coefficient; synthesizing the stress results under the different load components after scaling down to obtain a stress spectrum of each node of the finite element model of the to-be-measured structural member; calculating a fatigue damage value of the to-be-measured structural member according to the stress spectrum of each node and the stress-life curve of the to-be-measured structural member. By applying the method and the device, the fatigue analysis result of the structural member in the wind generating set is more accurate.
Description
Technical field
The present invention relates to fatigue analysis method, particularly relate to a kind of method and analysis of fatigue device for carrying out analysis of fatigue to the structural member relating to bearing connection in wind power generating set.
Background technology
Wind generating set structure part mainly bears or the dynamic loading of transfer complex, and therefore fatigure failure is one of dominant failure mode of wind generating set structure part.The accuracy of computation of structure stress result to the fatigue damage position of accurate computation structure part and lesion size most important.
Bearing is link important in wind power generating set, has very strong non-linear load transfer characteristic.And this characteristic of bearing, to the calculating of structural member stress, and then produces important impact to the accuracy of its Analysis of Fatigue exactly.
When existing wind generating set structure part fatigue analysis method carries out finite element analysis to structural member, each sharing part of the load applies separately a ultimate load, calculate the node stress result of structural member, node stress directly combines with the tired sequential load after process, carries out Linear Amplifer and the damage of synthesizing to come calculating unit to node stress.
When the method calculates fatigue damage, a kind of ultimate stress result is only calculated for each sharing part of the load, does not consider the load stress difference that structural member causes because of bearing non-linear factor when positive and negative two different directions effects.But due to the nonlinear characteristic of bearing, to each sharing part of the load, even if practical structures part is under positive and negative two load effects of identical numerical value, the dynamic changes process feature of bearing also can be significantly different, cause same node location to have different stress states.Therefore, differentiation is not added to sharing part of the load action direction, fatigue calculations can be caused to occur very large error, even can cause the change of tired hotspot location.
Summary of the invention
The object of the present invention is to provide a kind of fatigue analysis method for the structural member in wind power generating set and device, consider bearing nonlinear influencing factors, by finite element method, Analysis of Fatigue is carried out to the structural member in wind power generating set, thus make fatigue calculations more accurate.
According to an aspect of the present invention, a kind of fatigue analysis method for the structural member in wind power generating set is provided.Described method comprises: for the structural member to be measured and relative structural member that relate to bearing connection set up finite element model; According to the multiple load working conditions pre-set, difference imposed load on described finite element model, each described multiple load working condition comprises the fatigue load statistics ultimate value of the single sharing part of the load in positive dirction or negative direction, the described sharing part of the load is power or the moment of X, Y or Z-direction, and described multiple load working condition comprises the set of the fatigue load statistics ultimate value of each sharing part of the load in positive dirction and negative direction; By finite element analysis, obtain each node stress result of structural member to be measured under each load working condition effect; Using pretreated for regularization fatigue load sequence as coefficient, scaled down is carried out to the stress result obtained, then the stress result under the different loads component after comparative example reduction synthesizes, and obtains the stress spectra of each node of described structural member finite element model to be measured; According to the stress spectra of described each node and the S-N curve of described structural member to be measured, calculate the fatigue damage value of described structural member to be measured.
According to a further aspect in the invention, a kind of analysis of fatigue device for the structural member in wind power generating set is provided.Described device comprises: modeling unit, sets up finite element model for structural member to be measured and relative structural member for relating to bearing connection; Load applying unit, for multiple load working conditions that basis pre-sets, difference imposed load on described finite element model, each described multiple load working condition comprises the fatigue load statistics ultimate value of the single sharing part of the load in positive dirction or negative direction, the described sharing part of the load is power or the moment of X, Y or Z-direction, and described multiple load working condition comprises the set of the fatigue load statistics ultimate value of each sharing part of the load in positive dirction and negative direction; Stress acquiring unit, for by finite element analysis, obtains each node stress result of structural member to be measured under each load working condition effect; Reduction and synthesis unit, for using pretreated for regularization fatigue load sequence as coefficient, scaled down is carried out to the stress result that described stress acquiring unit obtains, then the stress result under the different loads component after comparative example reduction synthesizes, and obtains the stress spectra of each node of described structural member finite element model to be measured; Fatigue mechanisms unit, for according to reduction and synthesis unit obtain the stress spectra of described each node and the S-N curve of described structural member to be measured, calculate the fatigue damage value of described structural member to be measured.
The present invention propose the fatigue analysis method for the structural member in wind power generating set and use the device of described method to have following beneficial effect:
First, the main components of wind power generating set is being carried out to finite element modeling and in calculating, considering the non-linear effects of bearing, stress calculating results can be made more accurate.
Secondly, to the finite element model built, apply respectively to add up ultimate value with the forward of each sharing part of the load and oppositely corresponding respectively fatigue load, and extract the stress result of each sharing part of the load in positive and negative both direction respectively for analysis of fatigue, to the situation of positive and negative two different directions effects of each sharing part of the load, the component stress state change that the non-linear force-transfer characteristic of bearing causes can be reflected exactly.
Again, by each sharing part of the load Regularization is used for analysis of fatigue, instead of the stress result of direct unit of account load, avoid the error of calculation because bearing rigidity causes in low load section nonlinear characteristic this feature the most obvious, result of calculation is more accurate.
Again, by point column processing to tired time series load, each sharing part of the load is divided into positive and negative two passages, thus combines with corresponding stress result, make the actual loading feature of Fatigue Simulation analysis and parts more identical, result of calculation is more accurate.
Accompanying drawing explanation
Fig. 1 illustrates the process flow diagram of the fatigue analysis method for the structural member in wind power generating set according to exemplary embodiment of the present invention;
Fig. 2 is the logic diagram of the analysis of fatigue device for the structural member in wind power generating set that exemplary embodiment of the present invention is shown;
Fig. 3 illustrates the example of the tired time series load of the emulation before performing regularization pre-service;
Fig. 4 illustrates the example of the fatigue load sequence that exemplary embodiment according to the present invention obtains after regularization pre-service.
Embodiment
Below in conjunction with accompanying drawing detailed description according to the fatigue analysis method for the structural member in wind power generating set of exemplary embodiment of the present invention and the device using described method.
It may be noted that the fatigue analysis method for structural member of the present invention is applicable to any analysis of fatigue of blower fan structure part relating to bearing and connect, and be not limited only to the concrete structure part example that provides in this application.In addition, the fatigue analysis method for structural member of the present invention is applicable to the analysis of fatigue of the structural member in the blower fan of any type, has good versatility.
Fig. 1 illustrates the process flow diagram of the fatigue analysis method for the structural member in wind power generating set according to exemplary embodiment of the present invention.
In order to describe the fatigue analysis method for the structural member in wind power generating set of the present invention intuitively, described exemplary embodiment is described using base plate as the example of structural member at this.
With reference to Fig. 1, in step S110, for the structural member to be measured and relative structural member relating to bearing connection sets up finite element model.
Particularly, in step S110, set up described finite element model according to the geometry of structural member to be measured, physics and material properties, and in described finite element model, adopt the bearer properties of non-linear unit simulation bearing only pressurized.
For base plate, solid element grid division can be adopted to prostheses such as base body, main shaft and pylons; Solid element grid division is adopted to the inside and outside circle of driftage bearing, and the characteristic of bearing only pressurized is simulated to the Link unit in the roller employing ANSYS software of driftage bearing.
After the finite element model loading aforementioned foundation, in step S120, according to the multiple load working conditions pre-set, difference imposed load on described finite element model, each described multiple load working condition comprises the fatigue load statistics ultimate value of the single sharing part of the load in positive dirction or negative direction, the described sharing part of the load is power or the moment of X, Y or Z-direction, and described multiple load working condition comprises the set of the fatigue load statistics ultimate value of each sharing part of the load in positive dirction and negative direction.
Particularly, can online under formulate described multiple load working condition by following process: according to physics during structural member actual motion and mechanical state, judge to cause structure member to occur the sharing part of the load of fatigue damage, define the load working condition of finite element analysis for these sharing parts of the load; After this, the sharing part of the load described in tired time series load is added up, determines the ultimate value of this sharing part of the load.
Multiple load working conditions that exemplary embodiment according to the present invention is formulated are listed in following table:
After this, in step S130, by finite element analysis, obtain each node stress result of structural member to be measured under each load working condition effect.Such as, for described base plate, the stress result under each sharing part of the load effect of hub centre can be calculated.
In step S140, using pretreated for regularization fatigue load sequence as coefficient, carry out scaled down to the stress result obtained, the stress result under the different loads component then after comparative example reduction synthesizes, and obtains the stress spectra of each node of described structural member finite element model to be measured.The stress result of described each sharing part of the load is corresponding with the sharing part of the load corresponding in tired time series load.The exemplary process obtaining the pretreated fatigue load sequence of described regularization is illustrated after a while with reference to Fig. 3 and Fig. 4.
After this, in step S150, according to the stress spectra of described each node and the S-N curve (i.e. S-N curve) of described structural member to be measured, calculate the fatigue damage value of described structural member to be measured.
Usually, the S-N curve of described structural member to be measured is the S-N curve of the S-N curve amendment utilizing material.
Particularly, step S150 can comprise: by the time dependent stress spectra obtained from step S140, add up according to rain flow method, in conjunction with the S-N curve of described structural member, and utilize Palmgren-Miner linear cumulative damage criterion, calculate in fan design lifetime, the fatigue damage value of described structural member under described fatigue load.It may be noted that, method for calculating fatigue damage value according to the stress spectra of described each node and the S-N curve of described structural member to be measured in the present invention is not limited to above-mentioned use rain flow method and Palmgren-Miner linear cumulative damage criterion, and any stress spectra according to each node of structural member to be measured and S-N curve in prior art can be used to calculate the algorithm of fatigue damage value.
The exemplary process obtaining described pretreated fatigue load sequence is illustrated hereinafter with reference to Fig. 3 and Fig. 4.
Fig. 3 illustrates the example of the tired time series load of the emulation before performing regularization pre-service.Fig. 4 illustrates the example of the fatigue load sequence that exemplary embodiment according to the present invention obtains after regularization pre-service.Shown in Fig. 3 and Fig. 4 be respectively hub centre the tired time series load of the emulation before regularization pre-service and through the pretreated fatigue load sequence data of regularization.
With reference to Fig. 3, the tired time series load data of described emulation comprises moment and six sharing part of the load data (as shown in Figure 3, totally 7 row), described six sharing part of the load data are respectively the moment sharing part of the load in the X, Y, Z direction (Mx, My, Mz) and the power sharing part of the load in the X, Y, Z direction (Fx, Fy, Fz).Such as, 0.26 second moment, torque load is in the X direction 216093.328Nm, torque load component is in the Y direction 177349.906Nm, torque load component is in z-direction-79174.828Nm, the power sharing part of the load is in the X direction 148906.063N, and the power sharing part of the load is in the Y direction 274.755N, and the power sharing part of the load is in z-direction-922222N.
According to an alternative embodiment of the invention, obtain described regularization pretreated fatigue load sequence by following process: obtain the tired time series load of emulation, the tired time series load of described emulation comprises multiple dummy load component values in each moment in simulation time; Further, to each dummy load component value in each moment in described emulation tired time series load, Regularization is performed.
Particularly, described Regularization comprises: by each sharing part of the load of original tired time series load in each moment, be divided into positive and negative two passages according to loading direction; Being the absolute value ratio that fatigue load that described fatigue load value and described finite element apply adds up ultimate value by value assignment corresponding for the passage consistent with described dummy load component direction, is zero by value assignment corresponding for the passage contrary with described dummy load component direction.
Through the example of pretreated fatigue load sequence data shown in Fig. 4.Wherein, the fatigue load data in each moment divide two row displays, comprise 12 row altogether: sharing part of the load Mx, the sharing part of the load-Mx, sharing part of the load My, the sharing part of the load-My, sharing part of the load Mz, the sharing part of the load-Mz, sharing part of the load Fx, the sharing part of the load-Fx, sharing part of the load Fy, the sharing part of the load-Fy, sharing part of the load Fz and the sharing part of the load-Fz.Can find out, through pre-service, for each sharing part of the load, each dummy load is split as the fractionation sharing part of the load value of positive and negative both direction, and the fractionation sharing part of the load value of described positive and negative both direction be counted as to apply with the absolute value ratio of the corresponding ultimate load of the described fractionation sharing part of the load and its direction.Such as, 0.26 second moment, the sharing part of the load Mz shown in Fig. 3 is split as sharing part of the load Mz and-Mz, is counted as 0.0000000000E+0 and 0.7917482800E-01 respectively in the diagram in corresponding row.
If structural member to be measured is also subject to permanent load effect, then the load working condition of described finite element analysis also needs the operating mode comprising permanent load effect.
According to a preferred embodiment of the invention, the described ultimate load of applying also comprises the permanent load of the pretightning force of such as quality load and bolt etc.Step S120 also comprises: apply described permanent load respectively to described finite element model, and step S130 also comprises: obtain the stress result of each node of described finite element model under described permanent load effect respectively.In addition, in described fatigue load sequence, also comprise the permanent load coefficient in each moment, described permanent load coefficient is 1.In step S150, be that the fatigue load sequence of the permanent load coefficient of 1 is as coefficient using being included in the forward of each moment and each sharing part of the load and oppositely corresponding respectively load ratio and value, scaled down and synthesis are carried out to the stress result obtained, obtains the stress spectra of described each node.
Fig. 2 is the logic diagram of the analysis of fatigue device for the structural member in wind power generating set that exemplary embodiment of the present invention is shown.
With reference to Fig. 2, the analysis of fatigue device for the structural member in wind power generating set of exemplary embodiment of the present invention comprises modeling unit 210, load applying unit 220, stress acquiring unit 230, reduction and synthesis unit 240 and Fatigue Damage Calculation unit 250.
Modeling unit 210 sets up finite element model for the structural member to be measured and relative structural member for relating to bearing connection.Particularly, modeling unit 210 can set up described finite element model according to the physics of the parts in structural member to be measured and relative structural member, geometry and material properties, and in described finite element model, adopt the bearer properties of non-linear unit simulation bearing only pressurized.
Multiple load working conditions that load applying unit 220 pre-sets for basis, difference imposed load on the described finite element model that modeling unit 210 is set up, wherein, each described multiple load working condition comprises the fatigue load statistics ultimate value of the single sharing part of the load in positive dirction or negative direction, the described sharing part of the load is power or the moment of X, Y or Z-direction, and described multiple load working condition comprises the set of the fatigue load statistics ultimate value of each sharing part of the load in positive dirction and negative direction.
Stress acquiring unit 230, for by finite element analysis, obtains each node stress result of structural member to be measured under each load working condition effect.
Reduction and synthesis unit 240 for using pretreated for regularization fatigue load sequence as coefficient, scaled down is carried out to the stress result that described stress acquiring unit obtains, then the stress result under the different loads component after comparative example reduction synthesizes, and obtains the stress spectra of each node of described structural member finite element model to be measured.The stress result of described each sharing part of the load is corresponding with the sharing part of the load corresponding in tired time series load.
Alternatively, reduction and synthesis unit 240 obtain described pretreated fatigue load sequence from fatigue load pretreatment unit (not shown).Wherein, described fatigue load pretreatment unit comprises in said device or is included in another device, or alternatively, described device also comprises described fatigue load pretreatment unit.
Particularly, described fatigue load pretreatment unit is for obtaining the tired time series load of emulation, described emulation tired time series load comprises multiple dummy load component values in each moment in simulation time, and to each dummy load component value in each moment in described emulation tired time series load, perform Regularization.Described Regularization can comprise: by each sharing part of the load of original tired time series load in each moment, positive and negative two passages are divided into according to loading direction, being the absolute value ratio of the ultimate load value that described fatigue load value and described finite element apply by value assignment corresponding for the passage consistent with described dummy load component direction, is zero by value assignment corresponding for the passage contrary with described dummy load component direction.
Fatigue damage value computing unit 250, for according to reduction and the stress spectra of described each node that obtains of synthesis unit 240 and the S-N curve of described structural member to be measured, calculates the fatigue damage value of described structural member to be measured.
Alternatively, if structural member to be measured is also subject to permanent load effect, then the load working condition of described finite element analysis also comprises the operating mode of permanent load effect.In the case, load applying unit 120 is also for applying described permanent load respectively to described finite element model, and described stress acquiring unit 230 is also for obtaining the stress result of each node of described finite element model under described permanent load effect respectively; In described fatigue load sequence, also comprise the permanent load coefficient in each moment, described permanent load coefficient is 1.
The present invention propose the fatigue analysis method for the structural member in wind power generating set and use the device of described method to have following beneficial effect:
First, the main components of wind power generating set is being carried out to finite element modeling and in calculating, considering the non-linear effects of bearing, stress calculating results can be made more accurate.
Secondly, to the finite element model built, apply respectively to add up ultimate value with the forward of each sharing part of the load and oppositely corresponding respectively fatigue load, and extract the stress result of each sharing part of the load in positive and negative both direction respectively for analysis of fatigue, to the situation of positive and negative two different directions effects of each sharing part of the load, the component stress state change that the non-linear force-transfer characteristic of bearing causes can be reflected exactly.
Again, by each sharing part of the load Regularization is used for analysis of fatigue, instead of the stress result of direct unit of account load, avoid the error of calculation because bearing rigidity causes in low load section nonlinear characteristic this feature the most obvious, result of calculation is more accurate.
Again, by point column processing to tired time series load, each sharing part of the load is divided into positive and negative two passages, thus combines with corresponding stress result, make the actual loading feature of Fatigue Simulation analysis and parts more identical, result of calculation is more accurate.
In addition, the fatigue analysis method that the present invention proposes has versatility, be applicable to the analysis of fatigue relating to all blower fan main components that bearing connects, and be not limited to some blower fan type or certain concrete structural member, to the design of blower fan structure part, there is directive significance widely.
It may be noted that the needs according to implementing, each step described can be split as more multi-step, also the part operation of two or more step or step can be combined into new step, to realize object of the present invention in the application.
Above-mentioned can at hardware according to method of the present invention, realize in firmware, or be implemented as and can be stored in recording medium (such as CD ROM, RAM, floppy disk, hard disk or magneto-optic disk) in software or computer code, or be implemented and will be stored in the computer code in local recording medium by the original storage of web download in remote logging medium or nonvolatile machine readable media, thus method described here can be stored in use multi-purpose computer, such software process on the recording medium of application specific processor or able to programme or specialized hardware (such as ASIC or FPGA).Be appreciated that, computing machine, processor, microprocessor controller or programmable hardware comprise and can store or receive the memory module of software or computer code (such as, RAM, ROM, flash memory etc.), when described software or computer code by computing machine, processor or hardware access and perform time, realize disposal route described here.In addition, when the code for realizing the process shown in this accessed by multi-purpose computer, multi-purpose computer is converted to the special purpose computer for performing the process shown in this by the execution of code.
The above; be only the specific embodiment of the present invention, but protection scope of the present invention is not limited thereto, is anyly familiar with those skilled in the art in the technical scope that the present invention discloses; change can be expected easily or replace, all should be encompassed within protection scope of the present invention.Therefore, protection scope of the present invention should be as the criterion with the protection domain of described claim.
Claims (15)
1. for a fatigue analysis method for the structural member in wind power generating set, it is characterized in that, described method comprises:
For the structural member to be measured and relative structural member that relate to bearing connection set up finite element model;
According to the multiple load working conditions pre-set, difference imposed load on described finite element model, each described multiple load working condition comprises the fatigue load statistics ultimate value of the single sharing part of the load in positive dirction or negative direction, the described sharing part of the load is power or the moment of X, Y or Z-direction, and described multiple load working condition comprises the set of the fatigue load statistics ultimate value of each sharing part of the load in positive dirction and negative direction;
By finite element analysis, obtain each node stress result of structural member to be measured under each load working condition effect;
Using pretreated for regularization fatigue load sequence as coefficient, scaled down is carried out to the stress result obtained, then the stress result under the different loads component after comparative example reduction synthesizes, and obtains the stress spectra of each node of described structural member finite element model to be measured;
According to the stress spectra of described each node and the S-N curve of described structural member to be measured, calculate the fatigue damage value of described structural member to be measured.
2. method according to claim 1, is characterized in that, the stress result of described each sharing part of the load is corresponding with the sharing part of the load corresponding in tired time series load.
3. method according to claim 2, is characterized in that, the process that described structural member to be measured and relative structural member for relating to bearing connection sets up finite element model comprises:
Set up described finite element model according to the geometry of structural member to be measured, physics and material properties, and in described finite element model, adopt the bearer properties of non-linear unit simulation bearing only pressurized.
4. method according to claim 3, is characterized in that, described regularization pretreated fatigue load sequence is obtained by following process:
Obtain the tired time series load of emulation, the tired time series load of described emulation comprises multiple dummy load component values in each moment in simulation time,
To each dummy load component value in each moment in described emulation tired time series load, perform Regularization:
By each sharing part of the load of original tired time series load in each moment, be divided into positive and negative two passages according to loading direction,
Being the absolute value ratio that fatigue load that described fatigue load value and described finite element apply adds up ultimate value by value assignment corresponding for the passage consistent with described dummy load component direction, is zero by value assignment corresponding for the passage contrary with described dummy load component direction.
5. the method according to any one of Claims 1 to 4, is characterized in that, if structural member to be measured is also subject to permanent load effect, then the load working condition of described finite element analysis also comprises the operating mode of permanent load effect.
6. method according to claim 5, is characterized in that, described method also comprises: apply described permanent load respectively to described finite element model, and obtains the stress result of each node of described finite element model under described permanent load effect respectively.
7. method according to claim 6, is characterized in that, in described fatigue load sequence, also comprise the permanent load coefficient in each moment, described permanent load coefficient is 1.
8. for an analysis of fatigue device for the structural member in wind power generating set, it is characterized in that, described device comprises:
Modeling unit, sets up finite element model for structural member to be measured and relative structural member for relating to bearing connection;
Load applying unit, for multiple load working conditions that basis pre-sets, difference imposed load on described finite element model, each described multiple load working condition comprises the fatigue load statistics ultimate value of the single sharing part of the load in positive dirction or negative direction, the described sharing part of the load is power or the moment of X, Y or Z-direction, and described multiple load working condition comprises the set of the fatigue load statistics ultimate value of each sharing part of the load in positive dirction and negative direction;
Stress acquiring unit, for by finite element analysis, obtains each node stress result of structural member to be measured under each load working condition effect;
Reduction and synthesis unit, for using pretreated for regularization fatigue load sequence as coefficient, scaled down is carried out to the stress result that described stress acquiring unit obtains, then the stress result under the different loads component after comparative example reduction synthesizes, and obtains the stress spectra of each node of described structural member finite element model to be measured;
Fatigue mechanisms unit, for according to reduction and synthesis unit obtain the stress spectra of described each node and the S-N curve of described structural member to be measured, calculate the fatigue damage value of described structural member to be measured.
9. device according to claim 8, is characterized in that, the stress result of described each sharing part of the load is corresponding with the sharing part of the load corresponding in tired time series load.
10. device according to claim 9, is characterized in that, geometry, physics and material properties that described modeling unit is used for according to structural member to be measured set up described finite element model.
11. devices according to claim 10, it is characterized in that, described reduction and synthesis unit obtain described pretreated fatigue load sequence from fatigue load pretreatment unit, and wherein, described fatigue load pretreatment unit comprises in said device or is included in another device.
12. devices according to claim 11, it is characterized in that, described fatigue load pretreatment unit is for obtaining the tired time series load of emulation, described emulation tired time series load comprises multiple dummy load component values in each moment in simulation time, and to each dummy load component value in each moment in described emulation tired time series load, perform Regularization:
By each sharing part of the load of original tired time series load in each moment, be divided into positive and negative two passages according to loading direction,
Being the absolute value ratio of the fatigue load ultimate value that described fatigue load value and described finite element apply by value assignment corresponding for the passage consistent with described dummy load component direction, is zero by value assignment corresponding for the passage contrary with described dummy load component direction.
Device described in 13. according to Claim 8 ~ 12 any one, is characterized in that, if structural member to be measured is also subject to permanent load effect, then the load working condition of described finite element analysis also comprises the operating mode of permanent load effect.
14. devices according to claim 13, it is characterized in that, described load applying unit is also for applying described permanent load respectively to described finite element model, and described stress acquiring unit is also for obtaining the stress result of each node of described finite element model under described permanent load effect respectively.
15. devices according to claim 14, is characterized in that, in described fatigue load sequence, also comprise the permanent load coefficient in each moment, described permanent load coefficient is 1.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410664081.5A CN104573172B (en) | 2014-11-19 | 2014-11-19 | The fatigue analysis method and analysis of fatigue device of structural member in wind power generating set |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410664081.5A CN104573172B (en) | 2014-11-19 | 2014-11-19 | The fatigue analysis method and analysis of fatigue device of structural member in wind power generating set |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104573172A true CN104573172A (en) | 2015-04-29 |
| CN104573172B CN104573172B (en) | 2017-10-31 |
Family
ID=53089226
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410664081.5A Active CN104573172B (en) | 2014-11-19 | 2014-11-19 | The fatigue analysis method and analysis of fatigue device of structural member in wind power generating set |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104573172B (en) |
Cited By (29)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104809311A (en) * | 2015-05-14 | 2015-07-29 | 东南大学 | Structural part remaining life predicting method based on multi-factor fusion correction |
| CN107220449A (en) * | 2017-06-07 | 2017-09-29 | 重庆大学 | A kind of power of fan curve modeling method based on piecewise nonlinear function |
| CN107590356A (en) * | 2017-10-31 | 2018-01-16 | 国电联合动力技术有限公司 | A kind of wind driven generation set main shaft bearing automatic Model Selection method and storage device |
| WO2018050010A1 (en) * | 2016-09-19 | 2018-03-22 | 舍弗勒技术股份两合公司 | Evaluation method and apparatus for remaining life of bearing |
| CN108153949A (en) * | 2017-12-14 | 2018-06-12 | 新疆金风科技股份有限公司 | The force calculation method and device of transmission belt in the drive component of unsymmetric structure |
| CN108241761A (en) * | 2016-12-26 | 2018-07-03 | 北京金风科创风电设备有限公司 | Method and device for determining fatigue damage of generator component |
| CN109299559A (en) * | 2018-10-08 | 2019-02-01 | 重庆大学 | A kind of Surface hardened layer gear wear and fatigue failure competition mechanism analysis method |
| CN109426685A (en) * | 2017-08-28 | 2019-03-05 | 长城汽车股份有限公司 | Parts analysis method |
| CN109614640A (en) * | 2018-10-30 | 2019-04-12 | 中国电力科学研究院有限公司 | A kind of large-scale wind electricity machine unit hub Prediction method for fatigue life and system |
| 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 |
| CN109977441A (en) * | 2017-12-28 | 2019-07-05 | 北京金风科创风电设备有限公司 | Method and device for determining stress of bearing bolt hole of wind generating set |
| CN109977432A (en) * | 2017-12-27 | 2019-07-05 | 江苏金风科技有限公司 | Pylon door opening fatigue safety coefficient calculation method and computing device |
| CN110390146A (en) * | 2019-07-04 | 2019-10-29 | 山东中车风电有限公司 | Wind turbines tower weld fatigue damage measurement method and product based on sector load |
| CN110874515A (en) * | 2018-08-31 | 2020-03-10 | 福建金风科技有限公司 | Method and equipment for determining fatigue damage of three-pile foundation |
| CN110879925A (en) * | 2018-09-05 | 2020-03-13 | 江苏金风科技有限公司 | Fatigue damage determination method and device for top flange of wind generating set |
| CN110895621A (en) * | 2018-08-23 | 2020-03-20 | 江苏金风软件技术有限公司 | Method and device for determining fatigue damage of tower circumferential weld of wind turbine generator |
| CN111080015A (en) * | 2019-12-19 | 2020-04-28 | 中国铁建重工集团股份有限公司 | Shield equipment and real-time service life prediction system and prediction method thereof |
| CN112231842A (en) * | 2020-10-26 | 2021-01-15 | 哈尔滨工程大学 | Unmanned aerial vehicle structure design method based on fatigue strength |
| CN112395727A (en) * | 2019-07-31 | 2021-02-23 | 比亚迪股份有限公司 | CAE-based vehicle endurance condition creating method and device |
| CN112883555A (en) * | 2021-01-25 | 2021-06-01 | 浙江运达风电股份有限公司 | Method for extracting multidirectional loads of all sections of blade of wind generating set |
| CN113090473A (en) * | 2019-12-23 | 2021-07-09 | 北京金风科创风电设备有限公司 | Wind generating set and fatigue damage calculation method and device thereof |
| CN113139250A (en) * | 2021-04-19 | 2021-07-20 | 珠海格力智能装备有限公司 | Golllin column fatigue life prediction method, device, equipment and computer readable medium |
| CN113357099A (en) * | 2021-06-07 | 2021-09-07 | 浙江大学 | Fatigue diagnosis and detection method for fan tower drum based on acceleration sensor |
| CN113468711A (en) * | 2020-03-30 | 2021-10-01 | 新疆金风科技股份有限公司 | Method and device for reducing load condition of wind turbine generator |
| CN113569347A (en) * | 2021-06-24 | 2021-10-29 | 南方科技大学 | Fan component damage analysis method and device and storage medium |
| CN113673023A (en) * | 2020-05-15 | 2021-11-19 | 中航西飞民用飞机有限责任公司 | Method for determining fatigue load spectrum of complex block structure of airplane |
| CN113864137A (en) * | 2021-12-06 | 2021-12-31 | 天津发现技术有限公司 | Fatigue life monitoring method and system for whole wind turbine generator |
| CN113887108A (en) * | 2021-10-20 | 2022-01-04 | 西安热工研究院有限公司 | Fatigue life prediction method for hub of large wind turbine generator |
| CN117828954A (en) * | 2024-03-04 | 2024-04-05 | 质子汽车科技有限公司 | Swing arm fatigue analysis method and system considering contact state and electronic equipment |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102508970A (en) * | 2011-11-09 | 2012-06-20 | 保定天威风电科技有限公司 | Finite element modeling method for hub of megawatt wind generating set |
| CN102831260A (en) * | 2012-07-27 | 2012-12-19 | 南车株洲电力机车研究所有限公司 | Hub strength calculation method of wind power generating unit |
| US20140336954A1 (en) * | 2013-05-13 | 2014-11-13 | Michael Bruyneel | Calculating Fatigue and Fatigue Failure of Structures |
-
2014
- 2014-11-19 CN CN201410664081.5A patent/CN104573172B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102508970A (en) * | 2011-11-09 | 2012-06-20 | 保定天威风电科技有限公司 | Finite element modeling method for hub of megawatt wind generating set |
| CN102831260A (en) * | 2012-07-27 | 2012-12-19 | 南车株洲电力机车研究所有限公司 | Hub strength calculation method of wind power generating unit |
| US20140336954A1 (en) * | 2013-05-13 | 2014-11-13 | Michael Bruyneel | Calculating Fatigue and Fatigue Failure of Structures |
Non-Patent Citations (4)
| Title |
|---|
| J.M. JONKMAN等: "Dynamics of offshore floating wind turbines-analysis of three concepts", 《WIND ENERGY》 * |
| WENBIN DONG等: "Time domain-based gear contact fatigue analysis of a wind turbine drivetrain under dynamic conditions", 《INTERNATIONAL JOURNAL OF FATIGUE》 * |
| 杜静等: "兆瓦级风力发电机组主轴疲劳分析", 《太阳能学报》 * |
| 高俊云等: "基于Ansys/FE-SAFE的风力机主轴疲劳计算分析", 《机械强度》 * |
Cited By (46)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104809311B (en) * | 2015-05-14 | 2017-11-03 | 东南大学 | A kind of structural member method for predicting residual useful life based on multifactor fusion amendment |
| CN104809311A (en) * | 2015-05-14 | 2015-07-29 | 东南大学 | Structural part remaining life predicting method based on multi-factor fusion correction |
| WO2018050010A1 (en) * | 2016-09-19 | 2018-03-22 | 舍弗勒技术股份两合公司 | Evaluation method and apparatus for remaining life of bearing |
| CN107843427A (en) * | 2016-09-19 | 2018-03-27 | 舍弗勒技术股份两合公司 | Method and device for evaluating residual life of bearing |
| CN108241761A (en) * | 2016-12-26 | 2018-07-03 | 北京金风科创风电设备有限公司 | Method and device for determining fatigue damage of generator component |
| CN107220449A (en) * | 2017-06-07 | 2017-09-29 | 重庆大学 | A kind of power of fan curve modeling method based on piecewise nonlinear function |
| CN109426685A (en) * | 2017-08-28 | 2019-03-05 | 长城汽车股份有限公司 | Parts analysis method |
| CN109426685B (en) * | 2017-08-28 | 2023-09-22 | 长城汽车股份有限公司 | Component fatigue analysis 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 |
| CN107590356A (en) * | 2017-10-31 | 2018-01-16 | 国电联合动力技术有限公司 | A kind of wind driven generation set main shaft bearing automatic Model Selection method and storage device |
| CN107590356B (en) * | 2017-10-31 | 2020-04-28 | 国电联合动力技术有限公司 | Automatic model selection method and storage device for main shaft bearing of wind generating set |
| CN108153949A (en) * | 2017-12-14 | 2018-06-12 | 新疆金风科技股份有限公司 | The force calculation method and device of transmission belt in the drive component of unsymmetric structure |
| CN109977432A (en) * | 2017-12-27 | 2019-07-05 | 江苏金风科技有限公司 | Pylon door opening fatigue safety coefficient calculation method and computing device |
| CN109977441B (en) * | 2017-12-28 | 2023-02-03 | 北京金风科创风电设备有限公司 | Method and device for determining stress of bearing bolt hole of wind generating set |
| CN109977441A (en) * | 2017-12-28 | 2019-07-05 | 北京金风科创风电设备有限公司 | Method and device for determining stress of bearing bolt hole of wind generating set |
| CN110895621A (en) * | 2018-08-23 | 2020-03-20 | 江苏金风软件技术有限公司 | Method and device for determining fatigue damage of tower circumferential weld of wind turbine generator |
| CN110895621B (en) * | 2018-08-23 | 2023-07-25 | 江苏金风软件技术有限公司 | Method and device for determining fatigue damage of tower girth weld of wind turbine generator |
| CN110874515A (en) * | 2018-08-31 | 2020-03-10 | 福建金风科技有限公司 | Method and equipment for determining fatigue damage of three-pile foundation |
| CN110874515B (en) * | 2018-08-31 | 2023-04-14 | 福建金风科技有限公司 | Method and equipment for determining fatigue damage of three-pile foundation |
| CN110879925A (en) * | 2018-09-05 | 2020-03-13 | 江苏金风科技有限公司 | Fatigue damage determination method and device for top flange of wind generating set |
| CN110879925B (en) * | 2018-09-05 | 2023-07-11 | 江苏金风科技有限公司 | Method and device for determining fatigue damage of top flange of wind generating set |
| CN109299559A (en) * | 2018-10-08 | 2019-02-01 | 重庆大学 | A kind of Surface hardened layer gear wear and fatigue failure competition mechanism analysis method |
| CN109614640A (en) * | 2018-10-30 | 2019-04-12 | 中国电力科学研究院有限公司 | A kind of large-scale wind electricity machine unit hub Prediction method for fatigue life and system |
| CN110390146A (en) * | 2019-07-04 | 2019-10-29 | 山东中车风电有限公司 | Wind turbines tower weld fatigue damage measurement method and product based on sector load |
| CN110390146B (en) * | 2019-07-04 | 2023-04-18 | 山东中车风电有限公司 | Method for calculating fatigue damage of welding seam of tower cylinder of wind turbine generator based on sector load and product |
| CN112395727A (en) * | 2019-07-31 | 2021-02-23 | 比亚迪股份有限公司 | CAE-based vehicle endurance condition creating method and device |
| CN111080015B (en) * | 2019-12-19 | 2023-07-25 | 中国铁建重工集团股份有限公司 | Shield equipment and real-time life prediction system and prediction method thereof |
| CN111080015A (en) * | 2019-12-19 | 2020-04-28 | 中国铁建重工集团股份有限公司 | Shield equipment and real-time service life prediction system and prediction method thereof |
| CN113090473A (en) * | 2019-12-23 | 2021-07-09 | 北京金风科创风电设备有限公司 | Wind generating set and fatigue damage calculation method and device thereof |
| CN113468711A (en) * | 2020-03-30 | 2021-10-01 | 新疆金风科技股份有限公司 | Method and device for reducing load condition of wind turbine generator |
| CN113468711B (en) * | 2020-03-30 | 2024-02-06 | 金风科技股份有限公司 | Method and equipment for reducing load working condition of wind turbine generator |
| CN113673023A (en) * | 2020-05-15 | 2021-11-19 | 中航西飞民用飞机有限责任公司 | Method for determining fatigue load spectrum of complex block structure of airplane |
| CN113673023B (en) * | 2020-05-15 | 2024-05-24 | 中航西飞民用飞机有限责任公司 | Method for determining fatigue load spectrum of complex block structure of aircraft |
| CN112231842B (en) * | 2020-10-26 | 2022-12-23 | 北京氢源智能科技有限公司 | Unmanned aerial vehicle structure design method based on fatigue strength |
| CN112231842A (en) * | 2020-10-26 | 2021-01-15 | 哈尔滨工程大学 | Unmanned aerial vehicle structure design method based on fatigue strength |
| CN112883555A (en) * | 2021-01-25 | 2021-06-01 | 浙江运达风电股份有限公司 | Method for extracting multidirectional loads of all sections of blade of wind generating set |
| CN112883555B (en) * | 2021-01-25 | 2024-03-29 | 运达能源科技集团股份有限公司 | Multi-directional load extraction method for each section of blade of wind generating set |
| CN113139250B (en) * | 2021-04-19 | 2023-08-11 | 珠海格力智能装备有限公司 | Method, device, equipment and computer readable medium for predicting fatigue life of tie bar |
| CN113139250A (en) * | 2021-04-19 | 2021-07-20 | 珠海格力智能装备有限公司 | Golllin column fatigue life prediction method, device, equipment and computer readable medium |
| CN113357099B (en) * | 2021-06-07 | 2022-09-06 | 浙江大学 | Fatigue diagnosis and detection method for fan tower drum based on acceleration sensor |
| CN113357099A (en) * | 2021-06-07 | 2021-09-07 | 浙江大学 | Fatigue diagnosis and detection method for fan tower drum based on acceleration sensor |
| CN113569347A (en) * | 2021-06-24 | 2021-10-29 | 南方科技大学 | Fan component damage analysis method and device and storage medium |
| CN113887108A (en) * | 2021-10-20 | 2022-01-04 | 西安热工研究院有限公司 | Fatigue life prediction method for hub of large wind turbine generator |
| CN113864137A (en) * | 2021-12-06 | 2021-12-31 | 天津发现技术有限公司 | Fatigue life monitoring method and system for whole wind turbine generator |
| CN117828954A (en) * | 2024-03-04 | 2024-04-05 | 质子汽车科技有限公司 | Swing arm fatigue analysis method and system considering contact state and electronic equipment |
| CN117828954B (en) * | 2024-03-04 | 2024-06-07 | 质子汽车科技有限公司 | Swing arm fatigue analysis method and system considering contact state and electronic equipment |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104573172B (en) | 2017-10-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104573172A (en) | Fatigue analysis method and fatigue analysis device of structural member in wind generating set | |
| CN110210044A (en) | Load prediction method and device of wind generating set | |
| Dykes et al. | Sensitivity analysis of wind plant performance to key turbine design parameters: a systems engineering approach | |
| Nathi et al. | Coupled structural/thermal analysis of disc brake | |
| Bediz et al. | Three dimensional dynamics of rotating structures under mixed boundary conditions | |
| CN103471854A (en) | Aero-engine body vibration characteristic analysis method | |
| Babu et al. | Stress analysis on steering knuckle of the automobile steering system | |
| CN106021768B (en) | The power distribution network of the access containing distributed generation resource simplifies modeling method | |
| CN103942091A (en) | MATLAB user-defined model and PSASP joint simulation excitation system simulation method and system | |
| US20200091764A1 (en) | System and method for a fast power network simulator | |
| Fei et al. | Transient probabilistic analysis for turbine blade-tip radial clearance with multi-component and multi-physics fields based on DCERSM | |
| Adegas et al. | Reduced-order LPV model of flexible wind turbines from high fidelity aeroelastic codes | |
| Annoni et al. | Wind farm flow modeling using input-output dynamic mode decomposition | |
| Erzan Topçu | PC‐based control and simulation of an electro‐hydraulic system | |
| Berg et al. | Force estimation via kalman filtering for wind turbine blade control | |
| CN115270351A (en) | Method and device for calculating wall thickness of high-temperature superheater outlet header and electronic equipment | |
| Song et al. | Parametric fluid-structural interaction reduced order models in continuous time domain for aeroelastic analysis of high-speed vehicles | |
| Solcia et al. | Efficient multirate simulation of complex multibody systems based on free software | |
| CN114117847A (en) | Method for predicting fatigue life of hub tapered roller bearing | |
| CN106960077A (en) | A kind of modeling and simulation method of gear case of blower kinetic model | |
| Cai et al. | A comparison and simulation study of nonlinearly designed robust controllers for power system transient stability | |
| CN105334749A (en) | Value resolving method of fault correlation voting system reliability mathematic model | |
| CN108393928B (en) | Multi-rigid-finite-element hybrid analysis method for contact collision of mechanical arm of flexible robot | |
| Fan et al. | Power system equilibrium tracing and bifurcation detection based on the continuation of the recursive projection method | |
| CN107220459B (en) | Finite element analysis method for push rod of bulldozer |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CP01 | Change in the name or title of a patent holder |
Address after: 830026 No. 107, Shanghai Road, Urumqi economic and Technological Development Zone, the Xinjiang Uygur Autonomous Region Patentee after: Jinfeng Technology Co.,Ltd. Address before: 830026 No. 107, Shanghai Road, Urumqi economic and Technological Development Zone, the Xinjiang Uygur Autonomous Region Patentee before: XINJIANG GOLDWIND SCIENCE & TECHNOLOGY Co.,Ltd. |
|
| CP01 | Change in the name or title of a patent holder |