CN113357099A - Fatigue diagnosis and detection method for fan tower drum based on acceleration sensor - Google Patents

Fatigue diagnosis and detection method for fan tower drum based on acceleration sensor Download PDF

Info

Publication number
CN113357099A
CN113357099A CN202110632610.3A CN202110632610A CN113357099A CN 113357099 A CN113357099 A CN 113357099A CN 202110632610 A CN202110632610 A CN 202110632610A CN 113357099 A CN113357099 A CN 113357099A
Authority
CN
China
Prior art keywords
wind power
power tower
fatigue
tower cylinder
acceleration sensor
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
Application number
CN202110632610.3A
Other languages
Chinese (zh)
Other versions
CN113357099B (en
Inventor
从飞云
洪纵横
童水光
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN202110632610.3A priority Critical patent/CN113357099B/en
Publication of CN113357099A publication Critical patent/CN113357099A/en
Application granted granted Critical
Publication of CN113357099B publication Critical patent/CN113357099B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F03—MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
    • F03D—WIND MOTORS
    • F03D17/00—Monitoring or testing of wind motors, e.g. diagnostics
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F03—MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
    • F03D—WIND MOTORS
    • F03D13/00—Assembly, mounting or commissioning of wind motors; Arrangements specially adapted for transporting wind motor components
    • F03D13/20—Arrangements for mounting or supporting wind motors; Masts or towers for wind motors
    • 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
    • Y02E10/00—Energy generation through renewable energy sources
    • Y02E10/70—Wind energy
    • Y02E10/728—Onshore wind turbines

Abstract

The invention provides a fatigue diagnosis and detection method of a fan tower based on an acceleration sensor, which comprises the following steps: acquiring an acceleration signal of a wind power tower cylinder in the running process of the wind driven generator by using an acceleration sensor and calculating continuous disturbance displacement; simplifying the wind power tower cylinder into a cantilever beam model, and calculating external load acting on the top of the wind power tower cylinder and fatigue load acting on the bottom of the wind power tower cylinder; according to fatigue load taumaxCalculating and counting the stress amplitude of the fatigue load and the corresponding cycle number of the fatigue load by using a time history curve, and drawing a corresponding low-cycle fatigue load amplitude frequency histogram; calculating the total low-cycle fatigue damage of the wind power tower cylinder in the test time and predicting the low-cycle fatigue life of the wind power tower cylinder by using a Miner linear fatigue accumulated damage theory; the method utilizes the acceleration sensor for measurement, can calculate the fatigue damage of the wind power tower drum according to the measurement data, and further estimates the residual fatigue life of the wind power tower drumAnd convenience is provided for the work of maintenance personnel.

Description

Fatigue diagnosis and detection method for fan tower drum based on acceleration sensor
Technical Field
The invention relates to the technical field of wind power generation, in particular to a fatigue detection and fatigue life estimation method for a wind turbine tower.
Background
Nowadays, with the gradual depletion and the continuous rising of prices of traditional energy sources such as petroleum, coal, natural gas and the like, wind energy is regarded as a clean sustainable energy source and is receiving more and more attention of people. From the development trend, the wind energy is favored by various countries due to the advantages of huge storage amount, regeneration, wide distribution, no pollution and the like. Among them, wind power generation has gradually become the most scaled wind power utilization mode with promising commercialization development prospect.
The main structure of the wind driven generator comprises a tower, wind power blades, a cabin, a generator and a gear box in the cabin. The wind power tower cylinder is a tower pole for wind power generation, and mainly plays a supporting role in a wind generating set and absorbs the vibration of the set.
Along with the continuous increase of the single machine capacity of the wind turbine generator, the height of the tower barrel of the wind turbine generator is also continuously increased, so that more variable external loads are brought, and higher requirements are provided for the safe operation of the wind turbine generator. It should be noted that in some extreme weather conditions, the fatigue load experienced by the wind tower is large (exceeding the yield stress of the material), which causes the wind tower to deform plastically during the disturbance, and the actual fatigue life of the wind tower is lower than the design life. Therefore, low cycle fatigue is considered to be one of the important factors affecting the fatigue life of the wind turbine.
At present, most of fatigue diagnosis and detection methods for wind power towers in the market are time-consuming and labor-consuming, and the influence of low-cycle fatigue is not considered, for example, Chinese patent CN109340062B discloses a digital twin fatigue damage prediction method for a low-wind-speed wind turbine, the fatigue life and the fatigue damage condition of the wind turbine are predicted by establishing a wind wheel simulation model and a digital twin model of a virtual wind turbine, and the method is complex and needs to establish a plurality of models; for example, chinese patent CN202010878949.7 discloses a wind turbine blade multi-angle fatigue mechanical property detection device and a use method thereof, and the wind turbine blade multi-angle fatigue mechanical property detection device is arranged to analyze the wind turbine blade fatigue mechanical property when detecting different position states of the wind turbine blade.
Therefore, the low-cycle fatigue detection and diagnosis method for the wind power tower, which is high in applicability and precision, can guide the safe use of the wind driven generator and provide certain reference for the timely maintenance of the wind power tower in the future.
Disclosure of Invention
The invention aims to provide a fatigue diagnosis and detection method of a wind turbine tower based on an acceleration sensor.
In order to achieve the purpose, the invention provides the following technical scheme:
the invention provides a fatigue diagnosis and detection method for a fan tower cylinder of an acceleration sensor, which comprises the following steps:
(1) acquiring an acceleration signal of a wind power tower cylinder in the running process of the wind driven generator by using an acceleration sensor and calculating continuous disturbance displacement;
(2) simplifying the wind power tower cylinder into a cantilever beam model, and calculating external load acting on the top of the wind power tower cylinder and fatigue load acting on the bottom of the wind power tower cylinder;
(3) according to fatigue load taumaxCalculating and counting the stress amplitude of the fatigue load and the corresponding cycle number of the fatigue load by using a time history curve, and drawing a corresponding low-cycle fatigue load amplitude frequency histogram;
(4) and calculating the total low-cycle fatigue damage of the wind power tower cylinder and predicting the low-cycle fatigue life of the wind power tower cylinder by combining the S-N curve of the steel material of the wind power tower cylinder and the Miner linear fatigue accumulated damage theory.
Preferably, said fatigue load τmaxThe time history is plotted against the fatigue load calculated over a certain test time.
Furthermore, the acceleration sensor is a three-axis acceleration sensor and is vertically arranged on the inner wall of the top of the wind power tower; two mutually perpendicular measuring directions of the triaxial acceleration sensor parallel to the horizontal plane are set as an X direction and a Y direction respectively, and the measuring direction perpendicular to the horizontal plane is set as a Z direction.
Further, the specific calculation process of the continuous disturbance displacement in the step (1) is as follows:
(1.1) the three-axis acceleration sensor collects continuous time domain acceleration signals a in the X direction and the Y directionxAnd ay；
(1.2) carrying out low-pass filtering processing on continuous time domain acceleration signals acquired by the triaxial acceleration sensor;
the shake of the wind power tower barrel is mainly a low-frequency signal, and the low-frequency signal needs to be analyzed, so that the acquired acceleration signal needs to be subjected to low-pass filtering; high-frequency signals and other interference signals in the shaking process of the wind power tower can be filtered, and the accuracy of data acquisition is improved.
(1.3) carrying out low-pass filtering on the continuous time domain acceleration signal ax(t) and ay(t) performing integration twice to obtain a calculation formula L of continuous disturbance displacement waveform in X directionxCalculation formula L of continuous disturbance displacement waveform in (t) and Y directionsy(t) is:
Figure BDA0003104286090000041
Figure BDA0003104286090000042
Figure BDA0003104286090000043
Figure BDA0003104286090000044
wherein N is the number of collected data points;
i is one of the acquisition moments;
vx(t) is a continuous time domain velocity waveform in the X direction;
vy(t) is a continuous time domain velocity waveform in the X direction;
Lx(t) is the continuous perturbation displacement waveform in the X direction;
Ly(t) is a continuous disturbance displacement waveform in the Y direction
axiIs the acceleration sampling value at the time i in the X direction;
ayiis the acceleration sampling value at the time i in the Y direction;
vxiis the rate value at time i in the X direction;
vyiis the rate value at time i in the Y direction;
Δ t is the time difference between two samples.
Further, the step (2) specifically includes the following steps:
(2.1) because the tower height of the wind power tower is usually larger than 50m, and the change of the cross section of the wind power tower can be approximately ignored relative to the tower height, the disturbance of the wind power tower can be simplified into a cantilever beam with the equal cross section and subjected to horizontal load, and the deflection line differential equation and the boundary condition of the cantilever beam model are determined by combining relevant parameters of the wind power tower structure, such as the tower height, the outer diameter, the inner diameter and the like; the shearing force F acting on the bottom X direction of the wind power tower barrel is calculated by substituting the boundary condition into a deflection line differential equationxAnd the shearing force F acting on the Y direction of the bottom of the wind power towery：
X-direction deflection line equation:
Figure BDA0003104286090000051
boundary conditions in the X direction: when x is 0, L is 0;
when x is h, L is Lx；
Y-direction deflection line equation:
Figure BDA0003104286090000052
y-direction boundary conditions: when y is 0, L is 0;
when y is h, L is Ly；
Combining the above equations yields:
Figure BDA0003104286090000053
Figure BDA0003104286090000054
wherein x and y are two mutually perpendicular measuring directions of the three-axis acceleration sensor parallel to the horizontal plane respectively;
e is the Young modulus of the steel material of the wind power tower cylinder;
i is the section inertia moment of the wind power tower;
h is the tower height of the wind power tower;
Lxand LyRespectively measuring and calculating continuous disturbance displacement of the wind power tower cylinder in two mutually perpendicular horizontal directions by using an acceleration sensor;
the method for calculating the section moment of inertia I of the wind power tower barrel comprises the following steps:
Figure BDA0003104286090000055
wherein, pi is the circumferential ratio;
d is the outer diameter of the wind power tower;
d is the inner diameter of the wind power tower cylinder.
(2.2) calculating the actual shearing force F at the bottom of the wind power towersCombining the stress condition of the wind power tower cylinder and the actual shearing force F at the bottom of the wind power tower cylindersIs the sum of the two shear force vectors in the horizontal direction, and the actual shear force F at the bottom of the wind power tower cylindersThe calculation formula of (A) is as follows:
Figure BDA0003104286090000061
(2.3) the magnitude of fatigue load acting on the joint of the bottom of the wind power tower and the flange is equal to the maximum shear stress borne by the bottom of the wind power tower; the fatigue load tau acting on the bottom of the wind power tower barrelmaxThe calculation formula of (A) is as follows:
Figure BDA0003104286090000062
wherein, FsThe actual shear force at the bottom of the wind power tower cylinder;
a is the cross section area of the bottom of the cantilever beam model;
preferably, the wind power tower cylinder is simplified into a thin-wall circular ring-shaped section beam with a uniform section, and the maximum shear stress acting on the wind power tower cylinder in the horizontal direction is generated on a neutral axis of the section at the bottom of the wind power tower cylinder, namely the joint of the bottom of the wind power tower cylinder and the flange; for a wind power tower, the joint between the bottom of the wind power tower and a flange is the part which is most likely to cause fatigue damage to the wind power tower due to stress concentration;
(2.4) the external load acting on the top of the wind power tower cylinder is as follows:
according to the stress condition of the wind power tower cylinder, the external load acting on the top of the wind power tower cylinder and the actual shearing force F borne by the bottom of the tower cylindersAre equal.
Further, the specific process of the step (3) is as follows:
(3.1) obtaining the fatigue load tau at the bottom of the wind power tower drum based on the stepsmax-a time history curve, counting fatigue loads by means of a rain flow counting method;
preferably, the rain flow counting method is one of the most commonly used methods for analyzing the fatigue stress spectrum at present. The method is based on a double-parameter method, and can simplify the actually measured load course into a plurality of load cycles for fatigue life estimation.
(3.2) determining a threshold value of the low cycle fatigue load of the wind power tower cylinder according to the yield stress of the wind power tower cylinder material;
preferably, the determination of the threshold value of the low cycle fatigue load refers to the maximum fatigue load that can be theoretically borne when the wind power tower cylinder is subjected to low cycle fatigue.
(3.3) screening the amplitude and frequency of low cycle fatigue load circulation of the wind power tower cylinder in the working process based on the set threshold value of the low cycle fatigue load, and drawing a corresponding amplitude frequency histogram;
and determining the threshold value of the fatigue load when the wind power tower cylinder is subjected to plastic deformation according to the yield stress of the wind power tower cylinder material. When the fatigue load is larger than the threshold value, the bottom of the wind power tower generates plastic deformation. And when the irreversible plastic deformation at the bottom of the wind power tower cylinder is accumulated to a certain degree, the wind power tower cylinder is subjected to low-cycle fatigue damage.
Further, the specific process of the step (4) is as follows:
the amplitude frequency histogram obtained in the step is combined with an S-N curve of the steel material of the wind power tower cylinder and a Miner linear fatigue accumulated damage theory to obtain the low-cycle fatigue total damage D accumulated for the low-cycle fatigue effect of the wind power tower cylinder within a certain operation time ttot：
Figure BDA0003104286090000071
Wherein k is the number of collected data points;
i is one of the acquisition moments;
Diis the low cycle fatigue damage of the wind turbine tower at moment i;
Nfiexpressed in stress amplitude σiFatigue life under action;
nirepresenting stress amplitude σiCorresponding actual cycle times;
the low cycle fatigue life T of the wind power tower cylinder is as follows:
Figure BDA0003104286090000081
the invention has the following beneficial effects:
(1) according to the technical scheme, the acceleration sensor is used as the measuring device, the fatigue damage of the wind power tower can be calculated according to the measured data, the residual fatigue life of the wind power tower is further estimated, convenience is provided for the work of maintenance personnel, and a large amount of manpower and material resources are reduced by timely diagnosing and maintaining the wind power tower in the environment with wide distribution space and complex operating environment of the wind power plant.
(2) According to the invention, the disturbance conditions of the wind power tower in three directions can be measured simultaneously by utilizing the arrangement of the three-axis acceleration sensor, and the corresponding continuous disturbance displacement can be obtained through further calculation. Compared with the combination of multiple sensors, the method has the advantages that the measuring process is simple and quick, the accuracy of results is guaranteed, and meanwhile, the time for processing data of the multiple sensors is effectively saved.
(3) In order to reduce the effect of the wind turbine on the tower, the low-frequency signal of the acceleration sensor is extracted by utilizing the low-pass filtering algorithm according to the influence of the mechanical vibration of the gearbox on the data of the acceleration sensor, so that the influence of other factors on the subsequent calculation process is reduced as much as possible, and the predicted result is more accurate.
Drawings
FIG. 1 is a schematic illustration of a particular installation diagram of a three-axis acceleration sensor of the present invention;
FIG. 2 is a flowchart of a method for diagnosing and detecting low cycle fatigue of a wind power tower based on an acceleration sensor according to an embodiment of the present invention.
FIG. 3 is a schematic diagram of deformation of a wind tower after receiving stress in a horizontal direction according to an embodiment of the invention.
FIG. 4 is a cross-sectional view of a wind tower of the present invention.
Fig. 5 is a schematic diagram of fatigue load-time based on rainflow counting method treatment provided by the embodiment of the invention.
In the figure: the wind power generation device comprises a fan blade fan 1, a wind power cabin 2, an acceleration sensor 3 and a wind power tower 4.
Detailed Description
The following detailed description of the embodiments of the present invention is provided in conjunction with the accompanying drawings, and it should be noted that the embodiments are merely illustrative of the present invention and should not be considered as limiting the invention, and the purpose of the embodiments is to make those skilled in the art better understand and reproduce the technical solutions of the present invention, and the protection scope of the present invention should be subject to the scope defined by the claims.
As shown in fig. 1, the invention provides a fatigue diagnosis and detection method for a wind turbine tower based on an acceleration sensor, and a main structure of a wind turbine includes a fan blade 1, a wind turbine nacelle 2, a wind turbine tower 4, a nacelle, and a generator and a gearbox in the nacelle. The acceleration sensor 3 is a three-axis acceleration sensor and is vertically arranged on the inner wall of the top of the wind power tower barrel 4; two mutually perpendicular measuring directions of the triaxial acceleration sensor parallel to the horizontal plane are set as an X direction and a Y direction respectively, and the measuring direction perpendicular to the horizontal plane is set as a Z direction.
As shown in fig. 2, the detection method includes the following steps:
s1, acquiring an acceleration signal of the wind power tower in the running process of the wind driven generator by using the acceleration sensor and calculating continuous disturbance displacement, wherein the specific calculation process is as follows:
s1.1, the three-axis acceleration sensor collects continuous time domain acceleration signals a in the X direction and the Y directionxAnd ay；
In the operation process of the wind driven generator, the complex operation condition can cause the wind power tower to generate complex vibration, but the main reason causing the fatigue damage of the wind power tower is the transverse vibration of the wind power tower, namely the vibration of the wind power tower on the XY plane. Therefore, the present embodiment only considers selecting the measurement results of the acceleration sensor in the X and Y directions for further processing.
S1.2, carrying out low-pass filtering processing on continuous time domain acceleration signals acquired by the triaxial acceleration sensor;
the shaking of the wind power tower cylinder is mainly a low-frequency signal, and the low-frequency signal needs to be analyzed, so that the collected acceleration signal needs to be subjected to low-pass filtering, the frequency of the low-pass filtering is 10HZ-100HZ, the high-frequency signal and other interference signals in the shaking process of the wind power tower cylinder can be filtered, and the accuracy of data collection is improved.
S1.3, for lowContinuous time domain acceleration signal a after passing filteringx(t) and ay(t) performing integration twice to obtain a calculation formula L of continuous disturbance displacement waveform in X directionxCalculation formula L of continuous disturbance displacement waveform in (t) and Y directionsy(t) is:
Figure BDA0003104286090000101
Figure BDA0003104286090000102
Figure BDA0003104286090000103
Figure BDA0003104286090000104
wherein N is the number of collected data points;
i is one of the acquisition moments;
vx(t) is a continuous time domain velocity waveform in the X direction;
vy(t) is a continuous time domain velocity waveform in the X direction;
Lx(t) is the continuous perturbation displacement waveform in the X direction;
Ly(t) is a continuous disturbance displacement waveform in the Y direction
axiIs the acceleration sampling value at the time i in the X direction;
ayiis the acceleration sampling value at the time i in the Y direction;
vxiis the rate value at time i in the X direction;
vyiis the rate value at time i in the Y direction;
Δ t is the time difference between two samples.
S2, simplifying the wind power tower cylinder into a cantilever beam model, calculating the external load acting on the top of the wind power tower cylinder and the fatigue load acting on the bottom of the wind power tower cylinder, and specifically comprising the following steps:
s2.1, because the tower height of the wind power tower is usually more than 50m, the change of the cross section of the wind power tower can be approximately ignored relative to the tower height, the disturbance of the wind power tower can be simplified into an equi-section cantilever beam which is subjected to horizontal load and a cantilever beam model after the wind power tower is simplified,
as shown in FIG. 3, a high-h wind tower is subjected to a horizontal stress P with an inclination horizontal distance L0；
Determining a deflection line differential equation and boundary conditions of the cantilever beam model by combining related parameters of the wind power tower cylinder structure; calculating the shearing force F acting on the X direction of the bottom of the wind power tower through a connected equation groupxAnd the shearing force F acting on the Y direction of the bottom of the wind power towery：
X-direction deflection line equation:
Figure BDA0003104286090000111
boundary conditions in the X direction: when x is 0, L is 0;
when x is h, L is Lx；
Y-direction deflection line equation:
Figure BDA0003104286090000112
y-direction boundary conditions: when y is 0, L is 0;
when y is h, L is Ly；
Combining the above equations yields:
Figure BDA0003104286090000121
Figure BDA0003104286090000122
wherein x and y are two mutually perpendicular measuring directions of the three-axis acceleration sensor parallel to the horizontal plane respectively;
e is the Young modulus of the steel material of the wind power tower cylinder;
i is the section inertia moment of the wind power tower;
h is the tower height of the wind power tower;
Lxand LyRespectively measuring and calculating continuous disturbance displacement of the wind power tower cylinder in two mutually perpendicular horizontal directions by using an acceleration sensor;
the method for calculating the section moment of inertia I of the wind power tower barrel comprises the following steps:
Figure BDA0003104286090000123
wherein, pi is the circumferential ratio;
d is the outer diameter of the wind tower, as shown in FIG. 4;
d is the inner diameter of the wind tower, as shown in FIG. 4.
S2.2, calculating the actual shearing force F at the bottom of the wind power towersCombining the stress condition of the wind power tower cylinder and the actual shearing force F at the bottom of the wind power tower cylindersIs the sum of the two shear force vectors in the horizontal direction, and the actual shear force F at the bottom of the wind power tower cylindersThe calculation formula of (A) is as follows:
Figure BDA0003104286090000124
s2.3, the fatigue load acting on the connecting part of the bottom of the wind power tower cylinder and the flange is equal to the maximum shear stress borne by the bottom of the wind power tower cylinder; the fatigue load tau acting on the bottom of the wind power tower barrelmaxThe calculation formula of (A) is as follows:
Figure BDA0003104286090000131
wherein, FsThe actual shear force at the bottom of the wind power tower cylinder;
a is the cross section area of the bottom of the cantilever beam model;
preferably, the wind power tower cylinder is simplified into a thin-wall circular ring-shaped section beam with a uniform section, and the maximum shear stress acting on the wind power tower cylinder in the horizontal direction is generated on a neutral axis of the section at the bottom of the wind power tower cylinder, namely the joint of the bottom of the wind power tower cylinder and the flange; for a wind power tower, the joint between the bottom of the wind power tower and a flange is the part which is most likely to cause fatigue damage to the wind power tower due to stress concentration;
s2.4, the calculation formula of the external load acting on the top of the wind power tower cylinder is as follows:
according to the stress condition of the wind power tower cylinder, the external load acting on the top of the wind power tower cylinder and the actual shearing force F borne by the bottom of the tower cylindersAre equal.
S3, according to fatigue load taumaxCalculating and counting the stress amplitude of the fatigue load and the corresponding cycle number of the fatigue load by a time history curve, and drawing a corresponding low-cycle fatigue load amplitude map, wherein the specific process is as follows:
s3.1, obtaining the fatigue load tau based on the stepsmax-a time history curve, counting fatigue loads by means of a rain flow counting method;
preferably, the rain flow counting method is one of the most commonly used methods for analyzing the fatigue stress spectrum at present. The method is based on a double-parameter method, and can simplify the actually measured load course into a plurality of load cycles for fatigue life estimation.
As shown in fig. 5, after the fatigue load spectrum is processed by the rain flow counting method, 2 full cycles 2-3-2 'and 5-6-5' and 2 half cycles 1-2-4 and 4-5-7 can be obtained. After the threshold value of the low cycle fatigue load is set, only one half cycle 1-2-4 of the low cycle fatigue load can be obtained after screening.
S3.2, determining a fatigue load threshold value when the wind power tower cylinder is subjected to low cycle fatigue according to the yield stress of the wind power tower cylinder material;
preferably, the determination of the threshold value of the low cycle fatigue load refers to the maximum fatigue load that can be theoretically borne when the wind power tower cylinder is subjected to low cycle fatigue.
S3.3, screening the amplitude and frequency of low cycle fatigue load circulation of the wind power tower cylinder in the working process based on the set threshold value of the low cycle fatigue load, and drawing a corresponding amplitude frequency histogram;
and determining the threshold value of the fatigue load when the wind power tower cylinder is subjected to plastic deformation according to the yield stress of the wind power tower cylinder material. When the fatigue load is larger than the threshold value, the bottom of the wind power tower generates plastic deformation. And when the irreversible plastic deformation at the bottom of the wind power tower cylinder is accumulated to a certain degree, the wind power tower cylinder is subjected to low-cycle fatigue damage.
S4, calculating the total low cycle fatigue damage of the wind power tower cylinder and predicting the low cycle fatigue life of the wind power tower cylinder by combining the S-N curve of the steel material of the wind power tower cylinder and the Miner linear fatigue accumulated damage theory, wherein the specific process is as follows:
the amplitude frequency histogram obtained in the step is combined with an S-N curve of the steel material of the wind power tower cylinder and a Miner linear fatigue accumulated damage theory to obtain the low-cycle fatigue total damage D accumulated for the low-cycle fatigue effect of the wind power tower cylinder within a certain operation time ttot：
Figure BDA0003104286090000141
Wherein k is the number of collected data points;
i is one of the acquisition moments;
Diis the low cycle fatigue damage of the wind turbine tower at moment i;
Nfiexpressed in stress amplitude σiFatigue life under action;
nirepresenting stress amplitude σiCorresponding actual cycle times;
the low cycle fatigue life T of the wind power tower cylinder is as follows:
Figure BDA0003104286090000151
while the preferred embodiments of the present application have been described, additional variations and modifications in those embodiments may occur to those skilled in the art once they learn of the basic inventive concepts. Therefore, it is intended that the appended claims be interpreted as including preferred embodiments and all alterations and modifications as fall within the scope of the application.

Claims (6)

1. A method for diagnosing and detecting low cycle fatigue of a wind power tower based on an acceleration sensor is characterized by comprising the following steps:
(1) acquiring an acceleration signal of a wind power tower cylinder in the running process of the wind driven generator by using an acceleration sensor and calculating continuous disturbance displacement;
(2) simplifying the wind power tower cylinder into a cantilever beam model, and calculating external load acting on the top of the wind power tower cylinder and fatigue load acting on the bottom of the wind power tower cylinder;
(3) according to fatigue load taumaxCalculating and counting the stress amplitude of the fatigue load and the corresponding cycle number of the fatigue load by using a time history curve, and drawing a corresponding low-cycle fatigue load amplitude frequency histogram;
(4) and calculating the total low-cycle fatigue damage of the wind power tower cylinder within the test time and predicting the low-cycle fatigue life of the wind power tower cylinder by combining the S-N curve of the steel material of the wind power tower cylinder and the Miner linear fatigue accumulated damage theory.
2. The method for diagnosing and detecting the low cycle fatigue of the wind power tower based on the acceleration sensor as claimed in claim 1, wherein the acceleration sensor is a three-axis acceleration sensor vertically mounted on the inner wall of the top of the wind power tower; two mutually perpendicular measuring directions of the triaxial acceleration sensor parallel to the horizontal plane are set as an X direction and a Y direction respectively, and the measuring direction perpendicular to the horizontal plane is set as a Z direction.
3. The method for diagnosing and detecting the low-cycle fatigue of the wind power tower based on the acceleration sensor as claimed in claim 1, wherein the specific calculation process of the continuous disturbance displacement in the step (1) is as follows:
(1.1) the three-axis acceleration sensor collects the X squareContinuous time domain acceleration signal a in direction YxAnd ay；
(1.2) carrying out low-pass filtering processing on continuous time domain acceleration signals acquired by the triaxial acceleration sensor;
(1.3) carrying out low-pass filtering on the continuous time domain acceleration signal ax(t) and ay(t) performing integration twice to obtain a calculation formula L of continuous disturbance displacement waveform in X directionxCalculation formula L of continuous disturbance displacement waveform in (t) and Y directionsy(t) is:
Figure FDA0003104286080000021
Figure FDA0003104286080000022
Figure FDA0003104286080000023
Figure FDA0003104286080000024
wherein N is the number of collected data points;
i is one of the acquisition moments;
vx(t) is a continuous time domain velocity waveform in the X direction;
vy(t) is a continuous time domain velocity waveform in the Y direction;
Lx(t) is the continuous perturbation displacement waveform in the X direction;
Ly(t) is a continuous disturbance displacement waveform in the Y direction
axiIs the acceleration sampling value at the time i in the X direction;
ayiis the acceleration sampling value at the time i in the Y direction;
vxiis the rate value at time i in the X direction;
vyiis the rate value at time i in the Y direction;
Δ t is the time difference between two samples.
4. The method for diagnosing and detecting the low-cycle fatigue of the wind power tower based on the acceleration sensor as claimed in claim 1, wherein the step (2) specifically comprises the following steps:
(2.1) determining a deflection line equation and boundary conditions of the cantilever beam model by combining the simplified cantilever beam model of the wind power tower cylinder with relevant parameters of the wind power tower cylinder structure, and calculating the shearing force F acting on the bottom X direction of the wind power tower cylinder by substituting the boundary conditions into a deflection line differential equationxAnd the shearing force F acting on the Y direction of the bottom of the wind power towery：
X-direction deflection line equation:
Figure FDA0003104286080000031
boundary conditions in the X direction: when x is 0, L is 0;
when x is h, L is Lx；
Y-direction deflection line equation:
Figure FDA0003104286080000032
y-direction boundary conditions: when y is 0, L is 0;
when y is h, L is Ly；
Combining the above equations yields:
Figure FDA0003104286080000033
Figure FDA0003104286080000034
wherein x and y are two mutually perpendicular measuring directions of the three-axis acceleration sensor parallel to the horizontal plane respectively;
e is the Young modulus of the steel material of the wind power tower cylinder;
i is the section inertia moment of the wind power tower;
h is the tower height of the wind power tower;
Lxand LyRespectively measuring and calculating continuous disturbance displacement of the wind power tower cylinder in two mutually perpendicular horizontal directions by using an acceleration sensor;
(2.2) calculating the actual shearing force F at the bottom of the wind power towers：
Figure FDA0003104286080000041
(2.3) the magnitude of fatigue load acting on the joint of the bottom of the wind power tower and the flange is equal to the maximum shear stress borne by the bottom of the wind power tower; the fatigue load tau acting on the bottom of the wind power tower barrelmaxThe calculation formula of (A) is as follows:
Figure FDA0003104286080000042
wherein, FsThe actual shear force at the bottom of the wind power tower cylinder;
a is the cross section area of the bottom of the cantilever beam model;
(2.4) the external load acting on the top of the wind power tower cylinder is as follows:
according to the stress condition of the wind power tower cylinder, the external load acting on the top of the wind power tower cylinder and the actual shearing force F borne by the bottom of the tower cylindersAre equal.
5. The method for diagnosing and detecting the low-cycle fatigue of the wind power tower based on the acceleration sensor as claimed in claim 1, wherein the specific process of the step (3) is as follows:
(3.1) obtaining the fatigue load tau at the bottom of the wind power tower drum based on the stepsmax-a time history curve, counting fatigue loads by means of a rain flow counting method;
(3.2) determining a threshold value of the low cycle fatigue load of the wind power tower cylinder according to the yield stress of the wind power tower cylinder material;
and (3.3) screening the amplitude and frequency of the low cycle fatigue load cycle of the wind power tower cylinder in the working process based on the set threshold value of the low cycle fatigue load, and drawing a corresponding amplitude frequency histogram.
6. The method for diagnosing and detecting the low-cycle fatigue of the wind power tower based on the acceleration sensor as claimed in claim 1, wherein the specific process of the step (4) is as follows:
the amplitude frequency histogram obtained in the step is combined with an S-N curve of the steel material of the wind power tower cylinder and a Miner linear fatigue accumulated damage theory to obtain the low-cycle fatigue total damage D accumulated for the low-cycle fatigue effect of the wind power tower cylinder within a certain operation time ttot：
Figure FDA0003104286080000051
Wherein k is the number of collected data points;
i is one of the acquisition moments;
Diis the low cycle fatigue damage of the wind turbine tower at moment i;
Nfiexpressed in stress amplitude σiFatigue life under action;
nirepresenting stress amplitude σiCorresponding actual cycle times;
the low cycle fatigue life T of the wind power tower cylinder is as follows:
Figure FDA0003104286080000052
CN202110632610.3A 2021-06-07 2021-06-07 Fatigue diagnosis and detection method for fan tower drum based on acceleration sensor Active CN113357099B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110632610.3A CN113357099B (en) 2021-06-07 2021-06-07 Fatigue diagnosis and detection method for fan tower drum based on acceleration sensor

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110632610.3A CN113357099B (en) 2021-06-07 2021-06-07 Fatigue diagnosis and detection method for fan tower drum based on acceleration sensor

Publications (2)

Publication Number Publication Date
CN113357099A true CN113357099A (en) 2021-09-07
CN113357099B CN113357099B (en) 2022-09-06

Family

ID=77532854

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110632610.3A Active CN113357099B (en) 2021-06-07 2021-06-07 Fatigue diagnosis and detection method for fan tower drum based on acceleration sensor

Country Status (1)

Country Link
CN (1) CN113357099B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113719431A (en) * 2021-11-03 2021-11-30 浙江中自庆安新能源技术有限公司 Method and system for measuring residual life of fan tower drum
CN114264484A (en) * 2021-11-30 2022-04-01 江西五十铃汽车有限公司 Method and system for fatigue endurance test of post-processor parts

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101981308A (en) * 2008-08-25 2011-02-23 三菱重工业株式会社 Device, method and program for adjusting restriction on operation of windmill
CN102980651A (en) * 2012-11-02 2013-03-20 华锐风电科技(集团)股份有限公司 Monitoring method and monitoring device and monitoring system of wind turbine generator condition
CN104122013A (en) * 2014-07-15 2014-10-29 西安交通大学 On-line monitoring method for large-scale wind power tower drum structure stress
CN104573172A (en) * 2014-11-19 2015-04-29 新疆金风科技股份有限公司 Fatigue analysis method and fatigue analysis device of structural member in wind generating set
CN107387333A (en) * 2017-09-08 2017-11-24 北京金风科创风电设备有限公司 Method and device for determining fatigue damage of tower component of wind generating set
CN107835899A (en) * 2015-07-03 2018-03-23 远景能源（江苏）有限公司 The method and its wind turbine loaded in prediction and control wind turbine
CN109726411A (en) * 2017-10-31 2019-05-07 北京万源工业有限公司 The calculation method of wind energy conversion system engine room structure fatigue strength
CN110414102A (en) * 2019-07-16 2019-11-05 大连理工大学 Oversea wind power generation structural fatigue analysis method
CN111488678A (en) * 2020-04-01 2020-08-04 浙江运达风电股份有限公司 Wind turbine generator tower accumulated fatigue damage assessment system and method
CN112685939A (en) * 2020-12-28 2021-04-20 大连理工大学 Offshore wind turbine foundation fatigue damage analysis method based on actual measurement

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101981308A (en) * 2008-08-25 2011-02-23 三菱重工业株式会社 Device, method and program for adjusting restriction on operation of windmill
CN102980651A (en) * 2012-11-02 2013-03-20 华锐风电科技(集团)股份有限公司 Monitoring method and monitoring device and monitoring system of wind turbine generator condition
CN104122013A (en) * 2014-07-15 2014-10-29 西安交通大学 On-line monitoring method for large-scale wind power tower drum structure stress
CN104573172A (en) * 2014-11-19 2015-04-29 新疆金风科技股份有限公司 Fatigue analysis method and fatigue analysis device of structural member in wind generating set
CN107835899A (en) * 2015-07-03 2018-03-23 远景能源（江苏）有限公司 The method and its wind turbine loaded in prediction and control wind turbine
CN107387333A (en) * 2017-09-08 2017-11-24 北京金风科创风电设备有限公司 Method and device for determining fatigue damage of tower component of wind generating set
CN109726411A (en) * 2017-10-31 2019-05-07 北京万源工业有限公司 The calculation method of wind energy conversion system engine room structure fatigue strength
CN110414102A (en) * 2019-07-16 2019-11-05 大连理工大学 Oversea wind power generation structural fatigue analysis method
CN111488678A (en) * 2020-04-01 2020-08-04 浙江运达风电股份有限公司 Wind turbine generator tower accumulated fatigue damage assessment system and method
CN112685939A (en) * 2020-12-28 2021-04-20 大连理工大学 Offshore wind turbine foundation fatigue damage analysis method based on actual measurement

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
莫继华: "近海风电机组单桩式支撑结构疲劳分析", 《中国优秀硕士学位论文全文数据库 (工程科技Ⅱ辑)》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113719431A (en) * 2021-11-03 2021-11-30 浙江中自庆安新能源技术有限公司 Method and system for measuring residual life of fan tower drum
CN113719431B (en) * 2021-11-03 2022-02-08 浙江中自庆安新能源技术有限公司 Method and system for measuring residual life of fan tower drum
CN114264484A (en) * 2021-11-30 2022-04-01 江西五十铃汽车有限公司 Method and system for fatigue endurance test of post-processor parts

Also Published As

Publication number Publication date
CN113357099B (en) 2022-09-06

Similar Documents

Publication Publication Date Title
Dong et al. Structural vibration monitoring and operational modal analysis of offshore wind turbine structure
CN113357099B (en) Fatigue diagnosis and detection method for fan tower drum based on acceleration sensor
Gong et al. Imbalance fault detection of direct-drive wind turbines using generator current signals
CN113836762B (en) Wind turbine and wind power plant digital mirror image simulation display system
CN106677996B (en) A kind of tower barrel of wind generating set abnormal vibration detection method
CN111488678B (en) Wind turbine tower accumulated fatigue damage evaluation system and method
Osgood et al. Full-scale modal wind turbine tests: comparing shaker excitation with wind excitation
CN104265577B (en) Wind generating set abnormal detection method based on compressed sensing
CN106815771A (en) A kind of long-term evaluation method of wind power plant load
CN109211390B (en) Safety test and evaluation method for vibration and strength of power transmission tower
CN111504587A (en) Vertical hydraulic generator guide bearing support vibration characteristic evaluation system and method
CN104215323A (en) Method for determining sensitivity of each sensor in mechanical equipment vibrating sensor network
CN114893360B (en) Method and system for identifying abnormal vibration of tower of wind turbine generator and monitoring running state
Devriendt et al. Continuous dynamic monitoring of an offshore wind turbine on a monopile foundation
Lin et al. Field measurement, model updating, and response prediction of a large-scale straight-bladed vertical axis wind turbine structure
CN202453182U (en) Fault diagnosis device of gearbox of wind generation set
Wang et al. EEMD-based videogrammetry and vibration analysis method for rotating wind power blades
CN112177865B (en) Method for solving marking noise and insufficient marks in fan fault detection
CN207600740U (en) A kind of monitoring system for wind-driven generator
CN103630604B (en) The recognition methods of centrifugal compressor half-opened impeller crack fault
CN110469460B (en) Fault detection method and system for wind driven generator
Manzato et al. Wind turbine model validation by full-scale vibration test
CN111173687A (en) On-line monitoring device and method for crack damage of wind power fan blade
CN106338307A (en) Wind turbine generator set state monitoring system having yawing state monitoring function and wind turbine generator set state monitoring method thereof
CN112560916B (en) Wind power tower barrel overturning intelligent diagnosis method based on tilt angle sensor information

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant