CN109614640A - A kind of large-scale wind electricity machine unit hub Prediction method for fatigue life and system - Google Patents
A kind of large-scale wind electricity machine unit hub Prediction method for fatigue life and system Download PDFInfo
- Publication number
- CN109614640A CN109614640A CN201811278093.9A CN201811278093A CN109614640A CN 109614640 A CN109614640 A CN 109614640A CN 201811278093 A CN201811278093 A CN 201811278093A CN 109614640 A CN109614640 A CN 109614640A
- Authority
- CN
- China
- Prior art keywords
- wheel hub
- hub
- stress
- fatigue life
- model
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 59
- 230000005611 electricity Effects 0.000 title claims abstract description 24
- 230000006378 damage Effects 0.000 claims abstract description 63
- 230000008569 process Effects 0.000 claims abstract description 27
- 238000004458 analytical method Methods 0.000 claims abstract description 22
- 238000002474 experimental method Methods 0.000 claims abstract description 10
- 238000004088 simulation Methods 0.000 claims description 49
- 238000006073 displacement reaction Methods 0.000 claims description 14
- 229910000831 Steel Inorganic materials 0.000 claims description 12
- 239000010959 steel Substances 0.000 claims description 12
- 238000012216 screening Methods 0.000 claims description 9
- 239000007787 solid Substances 0.000 claims description 4
- 125000004122 cyclic group Chemical group 0.000 claims description 3
- 238000009661 fatigue test Methods 0.000 claims description 3
- 238000012300 Sequence Analysis Methods 0.000 claims 1
- 238000005096 rolling process Methods 0.000 claims 1
- 238000004364 calculation method Methods 0.000 abstract description 7
- 238000004590 computer program Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 238000012545 processing Methods 0.000 description 5
- 238000004422 calculation algorithm Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 3
- 240000002853 Nelumbo nucifera Species 0.000 description 2
- 235000006508 Nelumbo nucifera Nutrition 0.000 description 2
- 235000006510 Nelumbo pentapetala Nutrition 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- 238000012731 temporal analysis Methods 0.000 description 2
- 238000000700 time series analysis Methods 0.000 description 2
- 208000027418 Wounds and injury Diseases 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 230000001771 impaired effect Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 208000014674 injury Diseases 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- 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
A kind of large-scale wind electricity machine unit hub Prediction method for fatigue life and system, comprising: based on the wheel hub finite element model being obtained ahead of time, carry out specific loading load analysis and the analysis of load-time sequence, obtain the stress of wheel hub；Stress based on the wheel hub carries out fatigue experiment and obtains the tired damage degree at each position of the wheel hub, and the tired damage degree is screened by pre-set confidence interval, obtains danger zone；The stress for obtaining the danger zone predicts wind-powered machine unit hub fatigue life by crucial planar process.Calculating process thickness of the invention provided by the invention combines, and the fatigue life of the first entire wheel hub of express delivery rough calculation, obtains danger zone, carries out part to the region later and counts carefully, not only improved the accuracy of calculating, but also shorten and calculate the time, submits working efficiency.
Description
Technical field
The present invention relates to new energy fields, and in particular to a kind of large-scale wind electricity machine unit hub Prediction method for fatigue life and is
System.
Background technique
With the development of Wind turbines enlargement, the strength of parts check of Wind turbines seems more important, and fatigue is strong
It is indispensable that degree analysis, which checks very important a part as strength of parts,.Wheel hub is as connection blade, the weight of main shaft
Component is wanted, its damage fatigue is reasonably analyzed, obtains durability, the reliability data of degree of precision, set as defined in Valuation Standard
Counting the service life is 20 years, to wheel hub meet 20 years Fatigue life designs, cost control and to the accident preventions of entire Wind turbines,
Maintenance and repair has important meaning.Existing forecasting fatigue method is cumbersome, and time-consuming, and calculated result is inaccurate, improves
Life prediction precision has vital effect to the operating status of entire Wind turbines.
Summary of the invention
Cumbersome in order to solve forecasting fatigue method in the presence of the prior art, time-consuming, and calculated result is inaccurate
The problem of, the present invention provides a kind of large-scale wind electricity machine unit hub Prediction method for fatigue life and system.
Present invention provide the technical scheme that
A kind of large-scale wind electricity machine unit hub Prediction method for fatigue life, comprising:
Based on the wheel hub finite element model being obtained ahead of time, specific loading load analysis and the analysis of load-time sequence are carried out,
Obtain the stress of wheel hub；
Stress based on the wheel hub carries out fatigue experiment and obtains the tired damage degree at each position of the wheel hub, will be described tired
Labor damage degree is screened by pre-set confidence interval, obtains danger zone；
The stress for obtaining the danger zone predicts wind-powered machine unit hub fatigue life by crucial planar process.
Preferably, the foundation of the wheel hub finite element model, comprising:
Based on the wheel hub mock-up, is simplified by local feature and handled, obtain wheel hub simulation model；
Based on wheel hub simulation model, by the way of long-range displacement constraint, constraint wheel hub model is obtained；
Based on the constraint wheel hub model, wheel hub finite element model is obtained by 3 D rendering software.
Preferably, described that wheel hub simulation model is obtained by 3 D rendering software based on the wheel hub mock-up, it wraps
It includes:
The wheel hub mock-up is simplified to by local feature and modified process, removes the connection spiral shell between each component of wheel hub
Keyhole obtains initial wheel hub mock-up；
Based on the initial wheel hub mock-up, pitch variable bearings outer ring and hub body binding are connected, by pitch variable bearings
Inner ring and the binding of blade prosthese connect, and main shaft and hub body binding connect, and the steel ball between pitch variable bearings Internal and external cycle uses mould
Quasi- ball is simulated, and the steel ball between main shaft bearing Internal and external cycle is simulated using simulation ball, obtains wheel hub simulation model.
Preferably, described to be based on wheel hub simulation model, by the way of long-range displacement constraint, constraint wheel hub model is obtained,
Include:
Based on wheel hub simulation model, by the long-range displacement constraint, position is arranged in the position other than pitch variable bearings；
Based on the wheel hub model, by the long-range displacement constraint, determine that pitch variable bearings outer ring is active position；
The position is connected with the active position by contacting, limitation pitch variable bearings coordinate system Z-direction turns
It is dynamic, make main shaft Complete Bind, obtains constraint wheel hub model.
Preferably, described to be based on the constraint wheel hub model, wheel hub finite element model is obtained by 3 D rendering software, is wrapped
It includes:
Bolt hole, chamfering are removed by 3 D rendering software based on the constraint wheel hub model, obtaining initial wheel hub has
Limit meta-model；
Based on the initial wheel hub finite element model, hub body is subjected to grid dividing using tetrahedron element, by leaf
Piece prosthese, pitch variable bearings, main shaft, main shaft bearing carry out grid dividing using Three-Dimensional 8-node solid structure unit, will simulate ball
It is built into Three-D limited strain bar unit, obtains wheel hub finite element model.
It is preferably, described that simulation ball is built into Three-D limited strain bar unit, comprising:
The simulation ball of pitch variable bearings is built into the first Three-D limited strain bar unit；
The simulation ball of main shaft bearing is built into the second Three-D limited strain bar unit；
Wherein, the first Three-D limited strain bar unit is equal with second Three-D limited strain bar unit number.
Preferably, described based on the wheel hub finite element model being obtained ahead of time, when carrying out specific loading load analysis and load
Between sequence analyze, obtain the unitstress of wheel hub, comprising:
The unit lotus in X-axis, Y-axis and Z-direction is loaded respectively on three blades of the wheel hub finite element model
It carries, obtains unit load stress；
Stress is obtained in conjunction with load-time sequence based on the unit load stress；
The load includes: moment of flexure and power.
Preferably, the fatigue that the stress progress fatigue experiment based on the wheel hub obtains each position of the wheel hub is impaired
The tired damage degree is screened by pre-set confidence interval, obtains danger zone by degree, comprising:
Using the stress as input quantity, the fatigue test that stress-cyclic fatigue is destroyed is carried out, data matched curve is obtained；
Matched curve obtains the tired damage degree at each position based on the data, and the tired damage degree is set with described
Letter section compares screening, is danger zone by the corresponding position of tired damage degree in confidence interval.
Preferably, the stress for obtaining the danger zone predicts wind-powered machine unit hub fatigue by crucial planar process
Service life, including；
The true stress of the danger zone is acquired, and the true stress is converted into scalar stress；
Based on the scalar stress application into Cyclic Stress fatigue rupture curve, the damage number under load is obtained；
Wind-powered machine unit hub fatigue life is predicted by the fatigue life budget formula based on the damage number.
Preferably, described to predict that wind-powered machine unit hub is tired by the fatigue life budget formula based on the damage number
It in the labor service life, is calculated as follows:
S=1/D；
In formula, D is fatigue damage value, and S is wind-powered machine unit hub fatigue life；
Wherein, D=n/N；
In formula, D is fatigue damage value, and n is the damage number under load, and N is preset cycle-index.
Preferably, the stress for obtaining the danger zone predicts wind-powered machine unit hub fatigue by crucial planar process
Service life, further includes:
According to preset quantity, predict that the wind-powered machine unit hub damages the area of most serious by fatigue life budget formula
Domain, and be ranked up.
A kind of large-scale wind electricity machine unit hub fatigue life predicting system, comprising:
Stress obtains module: based on the wheel hub finite element model being obtained ahead of time, carrying out specific loading load analysis and load
Time series analysis obtains the stress of wheel hub；
Danger zone screening module: the stress based on the wheel hub carries out fatigue experiment and obtains the tired of each position of the wheel hub
The tired damage degree is screened by pre-set confidence interval, obtains danger zone by labor damage degree；
Life prediction module: obtaining the stress of the danger zone, predicts that wind-powered machine unit hub is tired by crucial planar process
The labor service life.
Preferably, it includes wheel hub finite element model setting up submodule in module that the stress, which obtains,；
The finite element model setting up submodule includes:
Wheel hub simulation model establishes unit: being based on the wheel hub mock-up, is simplified by local feature and handled, taken turns
Hub simulation model；
Constraint wheel hub model foundation unit: it is constrained by the way of long-range displacement constraint based on wheel hub simulation model
Wheel hub model；
Finite element model establishes unit: being based on the constraint wheel hub model, it is limited to obtain wheel hub by 3 D rendering software
Meta-model.
Preferably, the finite element model establishes unit, comprising:
Initial wheel hub mock-up establishes subelement: the wheel hub mock-up is simplified and modified by local feature
Journey removes the link bolt hole between each component of wheel hub, obtains initial wheel hub mock-up；
Wheel hub simulation model establishes subelement: the initial wheel hub mock-up is based on, by pitch variable bearings outer ring and wheel hub
Ontology binding connection, pitch variable bearings inner ring and the binding of blade prosthese are connected, and main shaft and hub body binding connect, pitch variable bearings
Steel ball between Internal and external cycle is simulated using simulation ball, and the steel ball between main shaft bearing Internal and external cycle is carried out using simulation ball
Simulation, obtains wheel hub simulation model.
Compared with prior art, the invention has the benefit that
The present invention provides a kind of large-scale wind electricity machine unit hub Prediction method for fatigue life and systems, comprising: based on preparatory
The wheel hub finite element model of acquisition carries out specific loading load analysis and the analysis of load-time sequence, obtains the stress of wheel hub；Base
Fatigue experiment is carried out in the stress of the wheel hub and obtains the tired damage degree at each position of the wheel hub, and the tired damage degree is led to
Pre-set confidence interval screening is crossed, danger zone is obtained；The stress for obtaining the danger zone, it is pre- by crucial planar process
Survey wind-powered machine unit hub fatigue life.Calculating process thickness of the invention combines, the fatigue life of the first entire wheel hub of express delivery rough calculation,
Danger zone is obtained, part is carried out to the region later and is counted carefully, the accuracy of calculating has not only been improved, but also shortens and calculates the time,
Submit working efficiency.Calculating process provided by the invention is suitable for the fatigue mechanisms of wind-powered machine unit hub incessantly, is equally applicable to
The fatigue mechanisms of the mechanical parts such as main shaft, rack.In addition, calculation process is reusable.The present invention establishes fatigue life prediction
After calculation and analysis methods normal process, if Wind turbines are load change, it is only necessary to calculate and analyze in fatigue life prediction
Load-time sequence is replaced in method flow, is not needed to repeat and is established calculation process.
Detailed description of the invention
Fig. 1 is prediction technique flow chart of the invention；
Fig. 2 is wind-powered machine unit hub geometrical model；
Fig. 3 is wind-powered machine unit hub grid model；
Fig. 4 is that wind-powered machine unit hub specific loading deforms calculated result；
Fig. 5 is wind-powered machine unit hub specific loading stress calculating results；
Fig. 6 is wind turbine blade coordinate system.
Specific embodiment
For a better understanding of the present invention, the contents of the present invention are done further with example with reference to the accompanying drawings of the specification
Explanation.
Embodiment 1:
Fig. 1 is the distribution network voltage control method flow chart based on synchronous measure and sensitivity estimation, as shown in Figure 1, this
A kind of large-scale wind electricity machine unit hub Prediction method for fatigue life provided is provided, comprising:
S1, based on the wheel hub finite element model being obtained ahead of time, carry out specific loading load analysis and load-time sequence point
Analysis, obtains the stress of wheel hub
By some local features simplification and modification for not influencing global stress, remove a sequence of all bolts hole of wheel hub.
Pitch variable bearings outer ring and wheel hub binding are connected, pitch variable bearings inner ring and the binding of blade prosthese connect, main shaft and wheel
Hub binding connects, and the steel ball between pitch variable bearings Internal and external cycle is simulated using rod piece, the steel ball between main shaft bearing Internal and external cycle
It is simulated using rod piece.
To prevent pitch variable bearings inner ring relative to the rotation of pitch variable bearings outer ring, by the way of long-range displacement constraint, make
With position of the point other than pitch variable bearings, active position is in pitch variable bearings outer ring, by contact to connecting.It is limited around change
The rotation of paddle coordinate system Z axis, main shaft Complete Bind.
The finite element model for establishing wheel hub is drawn the threedimensional model of component by three-dimensional software, removes the bolt being not concerned with
Hole, chamfering reform into the finite element model for FEM calculation, and wherein wheel hub carries out grid dividing using tetrahedron element,
Blade prosthese, pitch variable bearings, main shaft, main shaft bearing carry out grid dividing using Three-Dimensional 8-node solid structure unit, simulate ball
Bar unit is strained using Three-D limited, for simulating the form of bearing power transmission and the non-linear rigidity of ball.Share 3 pitch axis
Holding has 2 ring raceways with 1 main shaft bearing, each pitch variable bearings, and each raceway is using identical as a main shaft bearing ball
The Three-D limited of group number strains bar unit, for example each raceway has 100 balls, then each raceway just uses 100 groups of three-dimensionals
Finite strain bar unit.Main shaft bearing simulation equally strains bar unit using the Three-D limited of identical as ball number group of number.
Specific loading load is carried out on the finite element model of wheel hub, is individually added on blade 1, blade 2, blade 3
Specific loading Mx, My, Mz, Fx, Fy, Fz are carried, totally 18 operating conditions, Mx: i.e. the moment of flexure of X-direction, Fx: i.e. the power of X-direction obtains 18
A specific loading calculated result file, what is obtained is the calculated stress of specific loading drag, for calculating fatigue damage.
Analysis of fatigue uses load-time sequence, each operating condition be containing 18 sharing part of the load, that is, 1Mx, 1My, 1Mz, 1Fx,
1Fy, 1Fz, 2Mx, 2My, 2Mz, 2Fx, 2Fy, 2Fz, 3Mx, 3My, 3Mz, 3Fx, 3Fy, 3Fz be on three blades moment of flexure and power with
The load history of the variation of time, load-time sequence refer to corresponding load, loadmeter first at the time of emulation occurs
The moment is classified as from 0-600 seconds, behind be classified as the load at corresponding moment.It is mutually tied for the stress under the specific loading with obtained in
It closes, obtains the stress under real load.
S2, the stress based on the wheel hub carry out fatigue experiment and obtain the tired damage degree at each position of the wheel hub, by institute
It states tired damage degree to screen by pre-set confidence interval, obtains danger zone.
The year frequency of occurrence for counting fatigue load operating condition, the true stress obtained by load-time sequence is for S-N
Curve, the number for obtaining the damage generation of the time series obtain final then multiplied by the year frequency of the time series
Damage frequency.
The S-N curve of wheel hub material therefor is calculated, S-N curve, that is, stress-cyclic fatigue damage curve S-N curve is in perseverance
Fatigue test is carried out under width symmetrical loading, passes through drawing data matched curve.
Confidence interval is determined first, and it is integrally tired to carry out quick wheel hub using fatigue life prediction finite element method (fem) analysis method
Labor analysis, determines danger zone, i.e., the region that larger point is concentrated is damaged in calculated result.
S3, the stress for obtaining the danger zone predict wind-powered machine unit hub fatigue life by crucial planar process
The node for extracting danger zone, is assessed in detail again using crucial planar process, final to determine that fatigue damage is maximum
Region and point, crucial planar process (critical plane) is also the computational algorithm of fatigue life a kind of, when compared to rough calculation
Algorithm, time-consuming for the algorithm, calculates more accurate.
Crucial planar process is a kind of Stress superposition method, converts scalar for the stress tensor in LOAD FOR destination file
Method.Calculating process thickness combines, and shortens and calculates the time.
The stress application for being converted to scalar can be obtained damage number under respective loads into SN curve, under all load
The sum of damage frequency n and the ratio of design cycle times N be fatigue damage D.
D=n/N,
Calculated fatigue damage value is sorted, it is damage maximum region that maximum preceding 50 points are hurt in reducible setting loss.
Fatigue life=1/ fatigue damage value.
Embodiment 2:
Based on same inventive concept, the present invention also provides a kind of large-scale wind electricity machine unit hub fatigue life prediction systems
System, comprising:
Stress obtains module: based on the wheel hub finite element model being obtained ahead of time, carrying out specific loading load analysis and load
Time series analysis obtains the stress of wheel hub；
Danger zone screening module: the stress based on the wheel hub carries out fatigue experiment and obtains the tired of each position of the wheel hub
The tired damage degree is screened by pre-set confidence interval, obtains danger zone by labor damage degree；
Life prediction module: obtaining the stress of the danger zone, predicts that wind-powered machine unit hub is tired by crucial planar process
The labor service life.
It includes wheel hub finite element model setting up submodule in module that the stress, which obtains,；
The finite element model setting up submodule includes:
Wheel hub simulation model establishes unit: being based on the wheel hub mock-up, is simplified by local feature and handled, taken turns
Hub simulation model；
Constraint wheel hub model foundation unit: it is constrained by the way of long-range displacement constraint based on wheel hub simulation model
Wheel hub model；
Finite element model establishes unit: being based on the constraint wheel hub model, it is limited to obtain wheel hub by 3 D rendering software
Meta-model.
The finite element model establishes unit, comprising:
Initial wheel hub mock-up establishes subelement: the wheel hub mock-up is simplified and modified by local feature
Journey removes the link bolt hole between each component of wheel hub, obtains initial wheel hub mock-up；
Wheel hub simulation model establishes subelement: the initial wheel hub mock-up is based on, by pitch variable bearings outer ring and wheel hub
Ontology binding connection, pitch variable bearings inner ring and the binding of blade prosthese are connected, and main shaft and hub body binding connect, pitch variable bearings
Steel ball between Internal and external cycle is simulated using simulation ball, and the steel ball between main shaft bearing Internal and external cycle is carried out using simulation ball
Simulation, obtains wheel hub simulation model.
The constraint wheel hub model foundation unit, comprising:
Subelement is arranged in position: wheel hub simulation model is based on, through the long-range displacement constraint, other than pitch variable bearings
Position be arranged position；
Active position determines subelement: determining pitch variable bearings by the long-range displacement constraint based on the wheel hub model
Outer ring is active position；
Constraint wheel hub model obtains subelement: being connected the position with the active position by contacting, limitation becomes
Paddle shaft holds the rotation of coordinate system Z-direction, makes main shaft Complete Bind, obtains constraint wheel hub model.
The finite element model establishes unit, comprising:
Initial finite element model establishes subelement: removing spiral shell by 3 D rendering software based on the constraint wheel hub model
Keyhole, chamfering obtain initial wheel hub finite element model；
Wheel hub finite element model establishes subelement: being based on the initial wheel hub finite element model, hub body is used four
Face body unit carries out grid dividing, and blade prosthese, pitch variable bearings, main shaft, main shaft bearing are used Three-Dimensional 8-node solid structure list
Member carries out grid dividing, and simulation ball is built into Three-D limited strain bar unit, obtains wheel hub finite element model.
The wheel hub finite element model establishes the building of the strain bar unit of the Three-D limited in subelement, sets in the steps below
It sets:
The simulation ball of pitch variable bearings is built into the first Three-D limited strain bar unit；
The simulation ball of main shaft bearing is built into the second Three-D limited strain bar unit；
Wherein, the first Three-D limited strain bar unit is equal with second Three-D limited strain bar unit number.
The stress obtains module, comprising:
Unit load stress acquiring unit: it is loaded respectively on three blades of the wheel hub finite element model along X-axis, Y
Unit load on axis and Z-direction obtains unit load stress；
Stress is obtained in conjunction with load-time sequence based on the unit load stress；
The load includes: moment of flexure and power.
The danger zone screening module, comprising:
Data matched curve acquiring unit: using the stress as input quantity, the fatigue examination that stress-cyclic fatigue is destroyed is carried out
It tests, obtains data matched curve；
Risk zontation unit: matched curve obtains the tired damage degree at each position based on the data, and will be described
Tired damage degree and the confidence interval compare screening, are danger by the corresponding position of tired damage degree in confidence interval
Region.
The life prediction module, including；
Scalar stress transmission unit: the true stress of the danger zone is acquired, and the true stress is converted into mark
Measure stress；
It damages number acquiring unit: based on the scalar stress application into Cyclic Stress fatigue rupture curve, being carried
Damage number under lotus；
Life estimate unit: wind-powered machine unit hub is predicted by the fatigue life budget formula based on the damage number
Fatigue life.
The calculating in Wind turbines hub fatigue service life in the life estimate unit, shown as the following formula:
S=1/D；
In formula, D is fatigue damage value, and S is wind-powered machine unit hub fatigue life；
Wherein, D=n/N；
In formula, D is fatigue damage value, and n is the damage number under load, and N is preset cycle-index.
The life prediction module, further includes:
Degree of injury sorting sub-module: according to preset quantity, the Wind turbines are predicted by fatigue life budget formula
Wheel hub damages the region of most serious, and is ranked up.
Embodiment 3:
By taking the wind-powered machine unit hub of a certain model as an example:
Life prediction is shown in steps are as follows:
Wheel hub geometrical model is cleared up, the geometrical model for clearing up completion is as shown in Figure 2；
Grid dividing is carried out to the geometrical model of Fig. 2, the grid model after division is as shown in Figure 3；
Contiguity constraint carrying out boundary condition constraint, component to the grid model of Fig. 3；
After specific loading load calculates, calculated result is as shown in Figure 4, Figure 5, extracts calculated result file；
Using fatigue life prediction finite element method (fem) analysis method, load history and 18 specific loading operating conditions are added
Static strength results calculated, wind power generation unit blade coordinate system is as shown in Figure 6.
Whole damage is calculated using maximum absolute value slab method first, obtains danger zone, the region is used later
Crucial Plat algorithm carries out part and counts carefully, and obtains final fatigue life.
Obviously, described embodiments are some of the embodiments of the present invention, instead of all the embodiments.Based on the present invention
In embodiment, all other implementation obtained by those of ordinary skill in the art without making creative efforts
Example, shall fall within the protection scope of the present invention.
It should be understood by those skilled in the art that, embodiments herein can provide as method, system or computer program
Product.Therefore, complete hardware embodiment, complete software embodiment or reality combining software and hardware aspects can be used in the application
Apply the form of example.Moreover, it wherein includes the computer of computer usable program code that the application, which can be used in one or more,
The computer program implemented in usable storage medium (including but not limited to magnetic disk storage, CD-ROM, optical memory etc.) produces
The form of product.
The application is referring to method, the process of equipment (system) and computer program product according to the embodiment of the present application
Figure and/or block diagram describe.It should be understood that every one stream in flowchart and/or the block diagram can be realized by computer program instructions
The combination of process and/or box in journey and/or box and flowchart and/or the block diagram.It can provide these computer programs
Instruct the processor of general purpose computer, special purpose computer, Embedded Processor or other programmable data processing devices to produce
A raw machine, so that being generated by the instruction that computer or the processor of other programmable data processing devices execute for real
The device for the function of being specified in present one or more flows of the flowchart and/or one or more blocks of the block diagram.
These computer program instructions, which may also be stored in, is able to guide computer or other programmable data processing devices with spy
Determine in the computer-readable memory that mode works, so that it includes referring to that instruction stored in the computer readable memory, which generates,
Enable the manufacture of device, the command device realize in one box of one or more flows of the flowchart and/or block diagram or
The function of being specified in multiple boxes.
These computer program instructions also can be loaded onto a computer or other programmable data processing device, so that counting
Series of operation steps are executed on calculation machine or other programmable devices to generate computer implemented processing, thus in computer or
The instruction executed on other programmable devices is provided for realizing in one or more flows of the flowchart and/or block diagram one
The step of function of being specified in a box or multiple boxes.
The above is only the embodiment of the present invention, are not intended to restrict the invention, all in the spirit and principles in the present invention
Within, any modification, equivalent substitution, improvement and etc. done, be all contained in apply pending scope of the presently claimed invention it
It is interior.
Claims (14)
1. a kind of large-scale wind electricity machine unit hub Prediction method for fatigue life characterized by comprising
Based on the wheel hub finite element model being obtained ahead of time, specific loading load analysis and the analysis of load-time sequence are carried out, is obtained
The stress of wheel hub；
Stress based on the wheel hub carries out fatigue experiment and obtains the tired damage degree at each position of the wheel hub, by the fatigue by
Damage degree is screened by pre-set confidence interval, obtains danger zone；
The stress for obtaining the danger zone predicts wind-powered machine unit hub fatigue life by crucial planar process.
2. a kind of large-scale wind electricity machine unit hub Prediction method for fatigue life as described in claim 1, which is characterized in that the wheel
The foundation of hub finite element model, comprising:
Based on the wheel hub mock-up, is simplified by local feature and handled, obtain wheel hub simulation model；
Based on wheel hub simulation model, by the way of long-range displacement constraint, constraint wheel hub model is obtained；
Based on the constraint wheel hub model, wheel hub finite element model is obtained by 3 D rendering software.
3. a kind of large-scale wind electricity machine unit hub Prediction method for fatigue life as claimed in claim 2, which is characterized in that the base
Wheel hub simulation model is obtained by 3 D rendering software in the wheel hub mock-up, comprising:
The wheel hub mock-up is simplified to by local feature and modified process, removes the connection bolt between each component of wheel hub
Hole obtains initial wheel hub mock-up；
Based on the initial wheel hub mock-up, pitch variable bearings outer ring and hub body binding are connected, by pitch variable bearings inner ring
It binds and connects with blade prosthese, main shaft and hub body binding connect, and the steel ball between pitch variable bearings Internal and external cycle is using simulation rolling
Pearl is simulated, and the steel ball between main shaft bearing Internal and external cycle is simulated using simulation ball, obtains wheel hub simulation model.
4. a kind of large-scale wind electricity machine unit hub Prediction method for fatigue life as claimed in claim 2, which is characterized in that the base
In wheel hub simulation model, by the way of long-range displacement constraint, constraint wheel hub model is obtained, comprising:
Based on wheel hub simulation model, by the long-range displacement constraint, position is arranged in the position other than pitch variable bearings；
Based on the wheel hub model, by the long-range displacement constraint, determine that pitch variable bearings outer ring is active position；
The position is connected with the active position by contacting, the rotation of pitch variable bearings coordinate system Z-direction is limited, makes
Main shaft Complete Bind obtains constraint wheel hub model.
5. a kind of large-scale wind electricity machine unit hub Prediction method for fatigue life as claimed in claim 2, which is characterized in that the base
In the constraint wheel hub model, wheel hub finite element model is obtained by 3 D rendering software, comprising:
Bolt hole, chamfering are removed by 3 D rendering software based on the constraint wheel hub model, obtain initial wheel hub finite element
Model；
Based on the initial wheel hub finite element model, hub body is subjected to grid dividing using tetrahedron element, by blade vacation
Body, pitch variable bearings, main shaft, main shaft bearing carry out grid dividing using Three-Dimensional 8-node solid structure unit, by simulation ball building
Bar unit is strained at Three-D limited, obtains wheel hub finite element model.
6. a kind of large-scale wind electricity machine unit hub Prediction method for fatigue life as claimed in claim 5, which is characterized in that described to incite somebody to action
Simulation ball is built into Three-D limited strain bar unit, comprising:
The simulation ball of pitch variable bearings is built into the first Three-D limited strain bar unit；
The simulation ball of main shaft bearing is built into the second Three-D limited strain bar unit；
Wherein, the first Three-D limited strain bar unit is equal with second Three-D limited strain bar unit number.
7. a kind of large-scale wind electricity machine unit hub Prediction method for fatigue life as described in claim 1, which is characterized in that the base
In the wheel hub finite element model being obtained ahead of time, specific loading load analysis and the analysis of load-time sequence are carried out, wheel hub is obtained
Unitstress, comprising:
The unit load in X-axis, Y-axis and Z-direction is loaded respectively on three blades of the wheel hub finite element model, is obtained
To unit load stress；
Stress is obtained in conjunction with load-time sequence based on the unit load stress；
The load includes: moment of flexure and power.
8. a kind of large-scale wind electricity machine unit hub Prediction method for fatigue life as described in claim 1, which is characterized in that the base
Fatigue experiment is carried out in the stress of the wheel hub and obtains the tired damage degree at each position of the wheel hub, and the tired damage degree is led to
Pre-set confidence interval screening is crossed, danger zone is obtained, comprising:
Using the stress as input quantity, the fatigue test that stress-cyclic fatigue is destroyed is carried out, data matched curve is obtained；
Matched curve based on the data obtains the tired damage degree at each position, and by the tired damage degree and the confidence area
Between compare screening, by the corresponding position of tired damage degree in confidence interval be danger zone.
9. a kind of large-scale wind electricity machine unit hub Prediction method for fatigue life as described in claim 1, which is characterized in that described to obtain
The stress for taking the danger zone predicts wind-powered machine unit hub fatigue life by crucial planar process, including；
The true stress of the danger zone is acquired, and the true stress is converted into scalar stress；
Based on the scalar stress application into Cyclic Stress fatigue rupture curve, the damage number under load is obtained；
Wind-powered machine unit hub fatigue life is predicted by the fatigue life budget formula based on the damage number.
10. a kind of large-scale wind electricity machine unit hub Prediction method for fatigue life as claimed in claim 9, which is characterized in that described
Wind-powered machine unit hub fatigue life is predicted by the fatigue life budget formula based on the damage number, is calculated as follows:
S=1/D；
In formula, D is fatigue damage value, and S is wind-powered machine unit hub fatigue life；
Wherein, D=n/N；
In formula, D is fatigue damage value, and n is the damage number under load, and N is preset cycle-index.
11. a kind of large-scale wind electricity machine unit hub Prediction method for fatigue life as described in claim 1, which is characterized in that described
The stress for obtaining the danger zone predicts wind-powered machine unit hub fatigue life by crucial planar process, further includes:
According to preset quantity, predict that the wind-powered machine unit hub damages the region of most serious by fatigue life budget formula, and
It is ranked up.
12. a kind of large-scale wind electricity machine unit hub fatigue life predicting system characterized by comprising
Stress obtains module: based on the wheel hub finite element model being obtained ahead of time, carrying out specific loading load analysis and load-time
Sequence analysis, obtains the stress of wheel hub；
Danger zone screening module: stress based on the wheel hub carry out fatigue experiment obtain the fatigue at each position of the wheel hub by
The tired damage degree is screened by pre-set confidence interval, obtains danger zone by damage degree；
Life prediction module: obtaining the stress of the danger zone, predicts the wind-powered machine unit hub tired longevity by crucial planar process
Life.
13. a kind of large-scale wind electricity machine unit hub fatigue life predicting system as claimed in claim 12, which is characterized in that described
It includes wheel hub finite element model setting up submodule in module that stress, which obtains,；
The finite element model setting up submodule includes:
Wheel hub simulation model establishes unit: being based on the wheel hub mock-up, is simplified by local feature and handled, obtain wheel hub mould
Analog model；
It constrains wheel hub model foundation unit: being based on wheel hub simulation model, by the way of long-range displacement constraint, obtain constraint wheel hub
Model；
Finite element model establishes unit: being based on the constraint wheel hub model, obtains wheel hub finite element mould by 3 D rendering software
Type.
14. a kind of large-scale wind electricity machine unit hub fatigue life predicting system as claimed in claim 13, which is characterized in that described
Finite element model establishes unit, comprising:
Initial wheel hub mock-up establishes subelement: the wheel hub mock-up is simplified to by local feature and is modified process,
The link bolt hole between each component of wheel hub is removed, initial wheel hub mock-up is obtained；
Wheel hub simulation model establishes subelement: the initial wheel hub mock-up is based on, by pitch variable bearings outer ring and hub body
Binding connection connects pitch variable bearings inner ring and the binding of blade prosthese, and main shaft and hub body binding connect, inside and outside pitch variable bearings
Steel ball between circle is simulated using simulation ball, and the steel ball between main shaft bearing Internal and external cycle carries out mould using simulation ball
It is quasi-, obtain wheel hub simulation model.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811278093.9A CN109614640A (en) | 2018-10-30 | 2018-10-30 | A kind of large-scale wind electricity machine unit hub Prediction method for fatigue life and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811278093.9A CN109614640A (en) | 2018-10-30 | 2018-10-30 | A kind of large-scale wind electricity machine unit hub Prediction method for fatigue life and system |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109614640A true CN109614640A (en) | 2019-04-12 |
Family
ID=66001722
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811278093.9A Pending CN109614640A (en) | 2018-10-30 | 2018-10-30 | A kind of large-scale wind electricity machine unit hub Prediction method for fatigue life and system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109614640A (en) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110348151A (en) * | 2019-07-17 | 2019-10-18 | 国电联合动力技术有限公司 | Wind turbines tower drum flange connects bolt design method and device |
| CN111625910A (en) * | 2020-05-29 | 2020-09-04 | 华东交通大学 | Static strength and fatigue constraint-based compliant mechanism topology optimization method |
| CN112329305A (en) * | 2020-11-05 | 2021-02-05 | 中车株洲电力机车研究所有限公司 | Method for determining position and polishing amount of soft belt of variable-pitch bearing of wind turbine generator |
| CN113887108A (en) * | 2021-10-20 | 2022-01-04 | 西安热工研究院有限公司 | Fatigue life prediction method for hub of large wind turbine generator |
| CN114021288A (en) * | 2021-11-19 | 2022-02-08 | 西安热工研究院有限公司 | Method for predicting service life of yaw bearing of wind turbine generator |
| CN116423527A (en) * | 2023-06-12 | 2023-07-14 | 之江实验室 | Construction method and device of mechanical arm, storage medium and electronic equipment |
-
2018
- 2018-10-30 CN CN201811278093.9A patent/CN109614640A/en active Pending
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110348151A (en) * | 2019-07-17 | 2019-10-18 | 国电联合动力技术有限公司 | Wind turbines tower drum flange connects bolt design method and device |
| CN111625910A (en) * | 2020-05-29 | 2020-09-04 | 华东交通大学 | Static strength and fatigue constraint-based compliant mechanism topology optimization method |
| CN111625910B (en) * | 2020-05-29 | 2022-06-07 | 华东交通大学 | Static strength and fatigue constraint-based compliant mechanism topology optimization method |
| CN112329305A (en) * | 2020-11-05 | 2021-02-05 | 中车株洲电力机车研究所有限公司 | Method for determining position and polishing amount of soft belt of variable-pitch bearing of wind turbine generator |
| CN112329305B (en) * | 2020-11-05 | 2023-10-20 | 中车株洲电力机车研究所有限公司 | Method for determining soft belt position and polishing amount of variable-pitch bearing of wind turbine generator |
| CN113887108A (en) * | 2021-10-20 | 2022-01-04 | 西安热工研究院有限公司 | Fatigue life prediction method for hub of large wind turbine generator |
| CN114021288A (en) * | 2021-11-19 | 2022-02-08 | 西安热工研究院有限公司 | Method for predicting service life of yaw bearing of wind turbine generator |
| CN114021288B (en) * | 2021-11-19 | 2024-02-20 | 西安热工研究院有限公司 | Method for predicting service life of yaw bearing of wind turbine generator |
| CN116423527A (en) * | 2023-06-12 | 2023-07-14 | 之江实验室 | Construction method and device of mechanical arm, storage medium and electronic equipment |
| CN116423527B (en) * | 2023-06-12 | 2023-09-01 | 之江实验室 | Construction method and device of mechanical arm, storage medium and electronic equipment |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109614640A (en) | A kind of large-scale wind electricity machine unit hub Prediction method for fatigue life and system | |
| Chengwei et al. | Dynamic parametric modeling-based model updating strategy of aeroengine casings | |
| CN110210044A (en) | Load prediction method and device of wind generating set | |
| Tinga et al. | Aligning PHM, SHM and CBM by understanding the physical system failure behaviour | |
| CN105008887A (en) | Turbine blade fatigue life analysis using non-contact measurement and dynamical response reconstruction techniques | |
| CN111339700B (en) | Method and device for evaluating fatigue damage of nuclear turbine blade and storage medium | |
| CN107977516B (en) | It is a kind of to consider that the Notched specimen A LOCAL STRESS-STRAIN of multiaxial loading disproportional degree determines method | |
| CN103646147A (en) | Method for comprehensively evaluating maturity of aerospace component | |
| CN110487578A (en) | Girder construction damnification recognition method based on end reaction and strain | |
| Wu et al. | Development of sensor placement optimization tool and application to large-span cable-stayed bridge | |
| Liu et al. | Crack growth analysis and fatigue life estimation in the piston rod of a Kaplan hydro turbine | |
| Mohanty et al. | An adaptive neuro fuzzy inference system model for studying free in plane and out of plane vibration behavior of curved beams | |
| Hakim et al. | Structural damage detection using soft computing method | |
| Arailopoulos et al. | Integrated reverse engineering strategy for large-scale mechanical systems: application to a steam turbine rotor | |
| Akintunde et al. | Singular value decomposition and unsupervised machine learning for virtual strain sensing: Application to an operational railway bridge | |
| Chondros | Variational formulation of a rod under torsional vibration for crack identification | |
| Kostka et al. | A simulation-based monitoring of a composite plate using an integrated vibration measurement system | |
| CN114021288B (en) | Method for predicting service life of yaw bearing of wind turbine generator | |
| Häckell | A holistic evaluation concept for long-term structural health monitoring | |
| Vettori et al. | Virtual sensing for wind turbine blade full field response estimation in operational modal analysis | |
| Nabiyan et al. | Structural identification for dynamic strain estimation in wind turbine towers | |
| Vettori et al. | A virtual sensing approach to operational modal analysis of wind turbine blades | |
| Flores Terrazas et al. | A streamline approach to multiaxial fatigue monitoring using virtual sensing | |
| Drăgoi et al. | Evaluating the lifetime of hydro-mechanical equipment for low head hydropower plants | |
| Cukor et al. | An inverse approach for load identification of cracked wind turbine components |
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 |