CN110232238B - Loader layout method and device, computer equipment and storage medium - Google Patents
Loader layout method and device, computer equipment and storage medium Download PDFInfo
- Publication number
- CN110232238B CN110232238B CN201910502187.8A CN201910502187A CN110232238B CN 110232238 B CN110232238 B CN 110232238B CN 201910502187 A CN201910502187 A CN 201910502187A CN 110232238 B CN110232238 B CN 110232238B
- Authority
- CN
- China
- Prior art keywords
- loader
- parameters
- blade
- range
- wind power
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 58
- 238000003860 storage Methods 0.000 title claims abstract description 13
- 239000011159 matrix material Substances 0.000 claims abstract description 124
- 238000005452 bending Methods 0.000 claims abstract description 53
- 238000012360 testing method Methods 0.000 claims abstract description 29
- 238000004590 computer program Methods 0.000 claims description 15
- 238000010586 diagram Methods 0.000 claims description 15
- 230000009471 action Effects 0.000 claims description 13
- 230000005284 excitation Effects 0.000 claims description 11
- 238000004422 calculation algorithm Methods 0.000 claims description 4
- 230000010365 information processing Effects 0.000 claims description 3
- 230000011218 segmentation Effects 0.000 claims description 2
- 238000009661 fatigue test Methods 0.000 description 14
- 238000009826 distribution Methods 0.000 description 9
- 238000004364 calculation method Methods 0.000 description 8
- 238000012546 transfer Methods 0.000 description 8
- 239000013598 vector Substances 0.000 description 6
- 230000006870 function Effects 0.000 description 5
- 238000013461 design Methods 0.000 description 4
- 238000005457 optimization Methods 0.000 description 4
- 238000013178 mathematical model Methods 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 239000011248 coating agent Substances 0.000 description 1
- 238000000576 coating method Methods 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000004134 energy conservation Methods 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 238000013100 final test Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- 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
- G06F2111/00—Details relating to CAD techniques
- G06F2111/06—Multi-objective optimisation, e.g. Pareto optimisation using simulated annealing [SA], ant colony algorithms or genetic algorithms [GA]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
- Wind Motors (AREA)
Abstract
The invention belongs to the field of computers, and particularly relates to a loader layout method, a loader layout device, computer equipment and a storage medium, wherein the method comprises the following steps: blade parameters of the wind power blade are obtained, and the position range of the loader distributed on the wind power blade and the mass range of the loader are determined; generating a parameter matrix of the loader according to the position range and the quality range, wherein the parameter matrix at least comprises a position matrix and a quality matrix; and calculating modal parameters of a mode formed by the wind power blades and the loader by taking the parameter matrix as a variable, and taking the parameters meeting preset constraint conditions as layout parameters of the loader. The method comprises the steps of firstly determining the position range of a loader on a wind blade and the mass range of the loader through modal analysis of the wind blade and a loader system, then further limiting the position of the loader on the wind blade and the mass of the loader through constraints of conditions such as bending moment, fixed frequency, exciting force and the like, obtaining optimized loader layout parameters, and ensuring the accuracy of wind blade test data.
Description
Technical Field
The invention belongs to the field of computers, and particularly relates to a loader layout method, a loader layout device, computer equipment and a storage medium.
Background
The large-scale wind generating set is a development trend of China and developed countries in the world in the coming years, the large-scale wind generating blade is also more slender, and the reliability problem is gradually concerned. One of the main reliability issues of wind blades is the fatigue loading of the blade, which is now tested experimentally to ensure that the blade can withstand its fatigue loading throughout its service life.
The conventional wind power blade usually uses bending moment load in the spanwise direction of the blade as a fatigue test basis, single-point excitation is applied to complete the fatigue test, the energy required by the large-scale blade fatigue test is larger, the precision is higher, and the problems of precision error, insufficient driving capability and the like exist by adopting a single-point test method. The method for testing the wind power blade by adopting the double-excitation resonance fatigue test has the defects that the full-size fatigue performance of the large blade can be detected, the blade double-excitation resonance test is lack of an effective loader layout method, a tester arranges loaders at different positions of the blade mainly based on experience and observes whether the load of a tested section can reach or approach a target load through strain sensors arranged on different sections, the method is poor in economy and large in matching error, the aim that the loads of multiple sections of the blade meet or approach the designed load is difficult to achieve, the actual service life condition of the blade cannot be accurately obtained, and the like.
Therefore, the existing double-excitation resonance test is not mature enough, and better loader layout setting cannot be provided, so that the technical problem that the actual service life condition of the blade cannot be accurately obtained in the test is caused.
Disclosure of Invention
The embodiment of the invention aims to provide a loader layout method, and aims to solve the technical problem that the actual service life condition of a blade cannot be accurately obtained in a test because the conventional double-excitation resonance test is not mature enough and cannot provide better loader layout setting.
The embodiment of the invention is realized in such a way that a loader layout method comprises the following steps:
blade parameters of the wind power blade are obtained, and the position range of the loader distributed on the wind power blade and the mass range of the loader are determined according to the blade parameters;
generating a parameter matrix of the loader by adopting an exhaustion method according to the position range and the quality range by adopting a preset step length, wherein the parameter matrix at least comprises a position matrix and a quality matrix;
and calculating modal parameters of a mode formed by the wind power blade and the loader by taking the parameter matrix as a variable, and taking parameters in the parameter matrix corresponding to the modal parameters meeting preset constraint conditions as loader layout parameters, wherein the loader layout parameters at least comprise the position of the loader on the wind power blade and the mass of the loader.
Another objective of an embodiment of the present invention is to provide a loader layout apparatus, including:
the information acquisition unit is used for acquiring blade parameters of the wind power blade;
the information processing unit is used for determining the position range of the loader distributed on the wind power blade and the mass range of the loader according to the blade parameters; generating a parameter matrix of the loader by adopting an exhaustion method according to the position range and the quality range by adopting a preset step length, wherein the parameter matrix at least comprises a position matrix and a quality matrix; calculating modal parameters of a mode formed by the wind power blade and the loader by taking the parameter matrix as a variable, and taking parameters in the parameter matrix corresponding to the modal parameters meeting preset constraint conditions as loader layout parameters, wherein the loader layout parameters at least comprise the position of the loader on the wind power blade and the mass of the loader;
and the information output unit is used for outputting the loader layout parameters.
It is another object of an embodiment of the present invention to provide a computer device, including a memory and a processor, where the memory stores a computer program, and the computer program, when executed by the processor, causes the processor to execute the steps of the loader layout method.
It is another object of the embodiments of the present invention to provide a computer-readable storage medium, wherein the computer-readable storage medium stores thereon a computer program, and when the computer program is executed by a processor, the computer program causes the processor to execute the steps of the loader layout method.
The embodiment of the invention provides a loader layout method, a loader layout device, computer equipment and a storage medium thereof.
Drawings
FIG. 1 is a flowchart of a loader layout method according to an embodiment of the present invention;
FIG. 2 is a discrete model diagram of a wind turbine blade double excitation fatigue test system provided in an embodiment of the present invention;
FIG. 3 is a wind turbine blade fatigue damage action coefficient diagram provided in an embodiment of the present invention;
FIG. 4 is a bending moment envelope curve diagram of a wind turbine blade fatigue testing system provided by the embodiment of the invention;
FIG. 5 is a schematic diagram of a loader layout apparatus according to an embodiment of the present invention;
fig. 6 is an internal structural diagram of a computer device according to an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is described in further detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention.
It will be understood that, as used herein, the terms "first," "second," and the like may be used herein to describe various elements, but these elements are not limited by these terms unless otherwise specified. These terms are only used to distinguish one element from another. For example, a first xx script may be referred to as a second xx script, and similarly, a second xx script may be referred to as a first xx script, without departing from the scope of the present application.
Fig. 1 is a flowchart of a loader layout method according to an embodiment of the present invention, which is described in detail below.
In the embodiment of the invention, a loader layout method is provided, and layout parameters of a loader are obtained by analyzing a system formed by wind power blades and the loader, wherein the method comprises the following steps:
in step S102, blade parameters of the wind power blade are obtained, and a position range where the loader is distributed on the wind power blade and a mass range of the loader are determined according to the blade parameters.
In the implementation of the invention, the wind power blade refers to a wind wheel blade on a wind driven generator, and before the wind power blade is put into use, the fatigue damage resistance test needs to be carried out on the blade so as to conveniently know the service life of the wind power blade and the damage resistance capability of the wind power blade; the blade parameters of the wind power blade refer to some inherent characteristics of the wind power blade, and at least comprise the length, the mass, the distribution relation of the mass on the length, the fatigue damage action coefficient of each section of the wind power blade, the modal parameters of the wind power blade and the like, and the modal parameters of the wind power blade at least comprise fixed frequency; the loader is equipment which can be fixed on the wind power blade and provides exciting force; the position range of the loader distributed on the wind power blade refers to the position range of the loader capable of being placed on the wind power blade, and the mass range of the recorder refers to the mass range of the loader applicable to the wind power blade and the whole testing system.
As an embodiment of the invention, blade parameters of the wind power blade are obtained, and the position range of the wind power blade where the loader can be placed is obtained by analyzing the blade parameters, generally, the region of the wind power blade which is easy to fatigue is removed, and the rest parts are the position range of the loader; and then selecting a proper loader mass range according to the field experimental equipment, wherein the proper loader mass range is used as an initial range value of the position range and the mass range in the loader layout parameters.
According to the embodiment of the invention, the blade parameters of the wind power blade are analyzed, the position range of the loader distributed on the wind power blade and the initial mass range of the loader are determined by combining specific experimental equipment, and the configuration parameters of the loader are determined subsequently.
In step S104, a parameter matrix of the loader is generated by an exhaustive method according to the position range and the quality range in a preset step length, where the parameter matrix at least includes a position matrix and a quality matrix.
In the embodiment of the present invention, the parameter matrix of the loader is a matrix formed by layout parameters of the loader on the wind power blade, which facilitates calculation of the layout parameters of the loader subsequently, and the position matrix is a matrix formed by position parameters of the loader distributed on the wind power blade, which is generally preset by a technician, and the position parameters are selected by a certain step length within the position range determined in step S102 to form a position matrix; the quality matrix is a matrix formed by parameters in the quality range determined in step S102, and is also a position matrix formed by selecting quality parameters in a certain step size in the quality range.
According to the embodiment of the invention, the parameters in the position range and the quality range are selected to form the matrix, so that the completeness of the blade parameter test can be ensured, and the calculation of the subsequent loader layout parameters is facilitated.
In step S106, with the parameter matrix as a variable, calculating a modal parameter of a modal composed of the wind turbine blade and the loader, and taking a parameter in the parameter matrix corresponding to the modal parameter that meets a preset constraint condition as a loader layout parameter, where the loader layout parameter at least includes a position of the loader on the wind turbine blade and a mass of the loader.
In the embodiment of the invention, a mode of the wind power blade with the loader is constructed, the mode parameters of the modes of the wind power blade and the loader are calculated by taking the parameter matrix determined in the step S104 as a variable, wherein the mode parameters at least comprise parameters such as fixed frequency, bending moment, exciting force and the like, the parameters are compared with the preset constraint condition, and the parameter group in the parameter matrix corresponding to the mode parameter meeting the preset constraint condition is taken as the layout parameter of the loader. In general, the loader layout parameters meeting the requirements are not only in one group, and the technician can further optimize or select the layout parameters according to the actual situation.
According to the embodiment of the invention, through setting the constraint conditions, the proper loader distribution position and the proper loader quality are selected from the position range and the quality range to be combined, so that the loader layout parameters are obtained.
The embodiment of the invention provides a loader layout method, which comprises the steps of firstly determining the position range of a loader on a wind blade and the mass range of the loader through modal analysis of the wind blade and a loader system, and then further limiting the position of the loader on the wind blade and the mass of the loader through constraints of conditions such as bending moment, fixed frequency, exciting force and the like to obtain optimized loader layout parameters and ensure the accuracy of wind blade test data.
In the loader layout method provided by the embodiment of the present invention, the obtaining blade parameters of a wind turbine blade, and determining a position range of a loader distributed on the wind turbine blade and a mass range of the loader according to the blade parameters includes:
acquiring a blade fatigue damage action coefficient diagram of the wind power blade, dividing a region, with a fatigue damage action coefficient larger than a preset value, of the wind power blade into a vulnerable region according to the blade fatigue damage action coefficient diagram, and taking a part, with the vulnerable region removed, of the wind power blade as a position range where the loader is distributed on the wind power blade;
and determining the mass range according to the mass which can be borne by the clamp of the test excitation equipment and the motor power.
In the embodiment of the invention, as shown in fig. 2, fig. 2 is a discrete model diagram of a wind turbine blade dual-excitation fatigue test system provided in the embodiment of the invention, a wind turbine blade is fixed on a cylindrical support through a bolt flange, a loader is arranged on the blade through a clamp, and it is assumed that the equivalent mass of two loaders is Δ m 1 And Δ m 2 At respective positions of Δ l 1 、Δl 1 +Δl 2 And the blade is divided into a plurality of discrete mass sections along the spanwise direction by adopting a finite element method. As an embodiment of the invention, the length of a wind power blade is 56 meters, the blade is divided into 26 sections, 27 nodes are total, the blade section parameters are shown in table 1, fig. 3 is a fatigue damage action coefficient diagram of the wind power blade provided in the embodiment of the invention, the fatigue test envelope range is determined to be between 0 meter and 38 meters, the suction surface shell at the position of 20 meters is a dangerous interface, the position of the dangerous interface is removed, and other positions in the envelope range are used as the distribution position range of a loader. As an embodiment of the invention, the mass constraint of the loader is determined, and the mass constraint range of each loader is determined to be 400 kg-Delta m according to the mass and the motor power which can be borne by the clamp of the test excitation equipment j ≤8000kg。
TABLE 1 blade parameters
According to the embodiment of the invention, the position range of the wind power blade suitable for placing the loader is obtained through the fatigue damage action coefficient diagram of the wind power blade and the fiber fatigue strength analysis result of each part of the blade main body under the ultimate acting force, the accuracy of test data is ensured, the quality of the loader is determined according to the fixture bearing degree of experimental equipment and the motor power, and the smooth test is ensured.
In the loader layout method provided in an embodiment of the present invention, generating a parameter matrix of the loader by an exhaustive method according to the position range and the quality range in a preset step length, where the parameter matrix at least includes a position matrix and a quality matrix includes:
in the position range, segmenting the position range by taking one end close to the wind power blade as an initial position and a preset length as a step length, and taking the distance between each segmentation point and the initial position as an element in the position matrix;
and in the quality range, constructing the quality matrix by taking the lower limit value of the quality range as an initial value and taking preset quality as a step length.
As an embodiment of the invention, after the position range of the loader distributed on the wind power blade is obtained, the distance in the range is selected by adopting a preset step length by taking one end of the range as a starting point, and then a position matrix is constructed. As a preferred embodiment of the invention, the distribution range of the loader on the wind power blade is 0-38 meters, and then one numerical value is selected for each half meter, and 76 numerical values are selected in total to construct a position matrix; similarly, the mass matrix is constructed by selecting values within the mass range, and as an embodiment of the invention, the mass constraint range of the loader is 400kg ≦ Δ m j Less than or equal to 8000kg, and selecting a quality parameter every 50kg to construct a quality matrix. Of course, in the embodiment of the present invention, the selection of the step size may be selected according to the actual situation of the test, and the smaller the step size is, the more accurate the final test result is.
In the embodiment of the invention, the position parameters and the quality parameters of the loader are listed through an exhaustion method, and the position matrix and the quality matrix are constructed, so that the calculation of the layout parameters of the loader is facilitated.
In the loader layout method provided in the embodiment of the present invention, the calculating a modal parameter of a modality formed by the wind turbine blade and the loader using the parameter matrix as a variable, and using a parameter in the parameter matrix corresponding to the modal parameter that meets a preset constraint condition as a loader layout parameter, where the loader layout parameter at least includes a position of the loader on the wind turbine blade and a mass of the loader, includes:
taking parameters in the quality matrix as the quality of the loader, and taking parameters in the position matrix as the position of the loader on the wind power blade to construct a plurality of different modes;
calculating modal parameters of the modal, and checking whether the modal parameters meet the preset constraint condition; the modal parameters at least comprise a bending moment value of the modal, and the preset constraint condition at least comprises that the bending moment value is within a preset range;
and when the modal parameters meet the preset constraint condition, taking the parameters in the quality matrix corresponding to the modal parameters as the quality parameters of the loader, and taking the parameters in the position matrix as the position parameters of the loader on the wind power blade.
In the embodiment of the invention, as shown in fig. 2, the wind power blade and the loader model are separated, the length and the mass of the wind power blade are divided, and the blade is dispersed into n parts of EI along the spanwise direction i For flexural rigidity, m i Is segment mass,/ i Is the length of the segment; obtaining a discrete equivalent mass matrix M and an equivalent length matrix L of the system along the spanwise direction of the blade:
the bending moment of the wind power blade fatigue test can be divided into two parts: the self-weight of the blade is generated and the loader is generated, the weight of the blade is continuously distributed to enable the bending moment to be distributed nonlinearly, but the mass of the loader is distributed discretely to enable the bending moment to be distributed linearly in a segmented mode. Under the combined action of exciting force and dead weight, the blade needs to ensure that the distribution of bending moment along the spanwise direction is matched with the test requirement, and the bending moment T i Expressed as:
in the formula: t is i The bending moment values of each point of the blade are obtained; m is a group of i The total mass of the discrete mass block of the blade and the loader; l is i For measuring point and mass m i The distance of (d); omega is the excitation frequency; y is i The deflection of the discretized blade is obtained; g is the acceleration of gravity.
As a preferred embodiment of the invention, the target bending moment M of the blade is known target The information satisfies the constraint conditions that the bending moment error is controlled within the range of +/-10% of the target bending moment, and the error of the 7 th to 17 th cross sections of the blade is controlled to be not less than 0 and not more than delta i The range of less than or equal to 10 percent, the error of the 18 th to 22 th cross sections is controlled to be less than or equal to 10 percent and less than or equal to delta i Less than or equal to 10 percent, frequency f c ≥0.4501HZ。
According to the embodiment of the invention, the modal parameters of the system including the bending moment of the system are calculated by analyzing the wind power blades and the loading modes, so that the bending moment of the system can meet the requirements.
In the loader layout method provided in the embodiment of the present invention, the modal parameters further include a fixed frequency of the modal and an excitation force of a loader in the modal, and the preset constraint condition further includes that a mass of the loader is within a specified mass range, a fixed frequency of the modal is within a specified frequency range, and the excitation force recorded therein is within a specified range.
In the embodiment of the invention, a transfer matrix method is adopted, discretization modeling is carried out by utilizing a two-node beam, a series of no-mass units and mass nodes are formed after the span direction of the blade is discretized, and each unitIs represented by 4 vectors, for the blade span i section: deflection y i Angle of rotation theta i Bending moment M i And shear force Q i . Suppose blade loader mass Δ m j In the ith section, the relation between the deformation and the load of the concentrated mass is obtained according to the engineering mechanics theory, and the ith mass block m of the blade is established i With the i-1 th mass block m i-1 The relationship between the state vectors, and is expressed in matrix form:
let the deformation of the cross-section i and the load array be { p i }, corresponding to a concentrated mass matrix, liang Juzhen of [ M i ]And [ S ] i ]The above formula is abbreviated as { P i }=[M i ][S i ]{P i-1 }. The wind power blade is dispersed into N elements and N +1 junction points along the spanwise direction. P ij Representing the bond sites of the ith and jth elements. Each bond site P ij Corresponding to a state vector Z ij The total transfer vector Z can be obtained n,n+1 =H n H n-1 …H 0 Z 0,1 =H·Z 0,1 ,H i And H is a total transfer matrix of the section, and the total transfer matrix of the section is a 4 multiplied by 4 order matrix. The root of the blade is the input end, and the point is the output end, obtains along the total transfer matrix equation of blade spanwise, expresses as:
for a wind power blade with a fixed root, boundary conditions are as follows: blade root y 0 (x,t)=0、θ 0 (x, t) =0, blade tip M N (x,t)=0、Q N (x, t) =0. Substituting the total transmission matrix equation to obtain a solution equation of the wind power blade frequency, wherein the solution equation is as follows:
namely the frequency function of the system, and the natural frequency omega of each order of the system can be obtained by solving the function by using the dichotomy i . Will frequency omega i And the state vector of the initial end can be obtained by a back transfer equation, the state vectors of all points can be obtained by point-by-point transfer, and the main vibration mode corresponding to the frequency is drawn, so that the natural frequency and the corresponding mode of each order of the blade are obtained.
When the blade moves under the action of the exciting force, in a vibration period, the mass of all discrete sections generates the maximum kinetic energy of T, and by utilizing the relation of energy conservation, the two-point exciting force applies work to the blade and the damping energy consumption is equal, so that the blade has the advantages that
Wherein, M i =λ i l i +Δm j ,λ i For the sectional mass linear density, /) i For segment length, Δ m j As mass of the loader, y i Is a mass m i Position shift, m 11 、m 22 Respectively, the moving mass of the loader, Y 1 、Y 2 Respectively, the blade amplitude of the loading node.
During blade fatigue test, the quality and the position of the loader on the surface of the blade are adjusted, and the bending moment errors of each point are ensured to be within an allowable range. Loader layout takes into account the following factors: dangerous section of blade, frequency, bending moment value, load. Selecting the mass and the installation position of the loader as basic variables, and using the basic variables as design variables in bending moment load matching; and establishing a mathematical model by taking the minimum total mass of the loader as an objective function and taking the bending moment error precision and the frequency range of each point along the span direction of the blade as constraint conditions. The loader mass to be added Δ m j And position Δ l j As basic variables, as optimization variables in the bending moment matching, the matrix is expressed as:
an objective function and a constraint condition; minimizing the root mean square value of all section errors to obtain an objective function, and minimizing the error delta of each point along the span direction of the blade i The isoparametric parameters are used as constraint conditions, and the mathematical model is expressed as:
f c ≥f target
m ai ≤Δm i ≤m bi
in the formula, L i The distance between each point of the blade and the root; t is i The target bending moment value of each point along the span direction of the blade is obtained; m is a group of i =Δm j =m i Is the mass of the discrete segment; y is i Is a discrete section deflection; delta is the maximum value of the bending moment error; delta. For the preparation of a coating T Is the maximum value of the energy error; f. of c To calculate a frequency value; f. of target Setting a target lower limit frequency value; u represents the number of excitation points; f j Is an exciting force; y is j The amplitude of the excitation point of the blade; m is ai 、m bi Is the counterweight mass interval limit.
As an embodiment of the invention, the loader quality parameter matrix M of n multiplied by 2 is generated by an exhaustion method exh And a m x 2 position parameter matrix L exh N is determined by the step size of the loader mass and the upper and lower limits, and m is determined by the step size of the loader position on the blade and the upper and lower limits. The section parameters of the blade and the information of the loader are substituted into a transfer matrix calculation formula (2) to obtainThe calculation matrix of 26 segments, and the relation between the first section and the 26 th section by multiplying the 26 matrices in order, further results in the formula (4) for ω. Will balance the weight matrix M exh Each row is assigned to the unknown quantity Δ m in equation (4) 1 、Δm 2 The range ω is set between 0 and 5 and the error is set to 0.001%, and the first order frequency of the corresponding blade loading system is solved by bisection.
According to the design bending moment parameter of the blade root, the initial value M of the bending moment is assigned a1 And substituting the formula (2) to obtain the bending moment and the corresponding deflection of each section. Initially solving a group of minimum quality values, and obtaining a loader layout result by using an optimization algorithm: a mass of 3500kg, 2500kg and a position of 25m, 36m. The theoretical calculation predicted frequency value of the blade loading system is 0.4791Hz, the test frequency is 0.4876Hz, and the maximum error between the bending moment measured by the test and the predicted bending moment is 1.74%.
The design bending moment values versus the theoretical calculated values are shown in table 2 and the available bending moment envelope curves are shown in figure 4. From the comparison of the data in table 2, the blade calculated bending moment value is larger than the equivalent value in the range of 10 meters to 34 meters, and the error is basically within 100% -110%. Within the envelope range of 4 m to 38 m, the error is controlled within 90 percent to 110 percent. The equivalent load of the blade is a concave curve, and the bending moment generated after the loader is added is basically in a linear relation, so that an error is generated. Because the bending moment value of each section of the blade is solved according to the initial value of the bending moment of the blade root, the total error of the bending moment can be effectively controlled by controlling the bending moment value of the root.
TABLE 2 comparison of bending moment results
According to the embodiment of the invention, the resonance fatigue test configuration of the wind power blade influences the system frequency and the driving force required by the loader, extra test time is increased when the performance of the loader is in short, and the driving force required by the loader is reduced near the minimum error of test and target bending moment by optimizing the test configuration, which is particularly important for the fatigue test of the large-size wind power blade. In order to enable the bending moment load distribution during actual work to be met along the blade span direction during fatigue test, a layout optimization mathematical model of the loader of the wind power blade double-excitation fatigue test system is established, the quality and the quantity of the loader are optimized, the layout and the verification are carried out, and the matching of the test bending moment distribution and a design value is ensured.
In the loader layout method provided in the embodiment of the present invention, the method further includes: and constraining all the loader layout parameters by using a quadratic programming algorithm by taking the mean square value of the modal bending moment as a constraint condition, and taking the loader layout parameter which enables the mean square value of the bending moment to be minimum as a final layout parameter of the loader.
In the embodiment of the present invention, generally, more than one set of loader layout parameters is determined by the above embodiment, and a set of loader parameter information satisfying the minimum mean square value of the bending moment envelope error is obtained by optimizing all obtained layout parameters by using a sequential quadratic programming algorithm in an MATLAB optimization toolbox, so as to obtain the optimal layout parameters.
According to the embodiment of the invention, through modal analysis of the wind power blade and the loader system, the position range of the loader on the wind power blade and the mass range of the loader are firstly determined, and then the position of the loader on the wind power blade and the mass of the loader are further limited through the constraints of conditions such as bending moment, fixed frequency, exciting force and the like, so that optimized loader layout parameters are obtained, and the accuracy of wind power blade test data is ensured.
Fig. 5 is a schematic diagram of a loader layout apparatus according to an embodiment of the present invention, which is described in detail below.
An embodiment of the present invention further provides a loader layout apparatus, including:
an information obtaining unit 501, configured to obtain blade parameters of a wind turbine blade;
the information processing unit 502 is configured to determine, according to the blade parameter, a position range where the loader is distributed on the wind turbine blade and a mass range of the loader; generating a parameter matrix of the loader by adopting an exhaustion method according to the position range and the quality range by adopting a preset step length, wherein the parameter matrix at least comprises a position matrix and a quality matrix; calculating modal parameters of a mode formed by the wind power blade and the loader by taking the parameter matrix as a variable, and taking parameters in the parameter matrix corresponding to the modal parameters meeting preset constraint conditions as loader layout parameters, wherein the loader layout parameters at least comprise the position of the loader on the wind power blade and the mass of the loader;
an information output unit 503, configured to output the loader layout parameter.
In the embodiment of the invention, blade parameters of the wind power blade are obtained, and the position range of the loader distributed on the wind power blade and the mass range of the loader are determined according to the blade parameters.
In the implementation of the invention, the wind power blade refers to a fan blade on a wind driven generator, and before the wind power blade is put into use, the fatigue damage resistance test needs to be carried out on the blade so as to conveniently know the service life of the wind power blade and the damage resistance capability of the wind power blade; the blade parameters of the wind power blade refer to some inherent characteristics of the wind power blade, and at least comprise the length, the mass, the distribution relation of the mass on the length, the fatigue damage action coefficient of each section of the wind power blade, the modal parameters of the wind power blade and the like, and the modal parameters of the wind power blade at least comprise fixed frequency; the loader can be fixed on the wind power blade and provides excitation force; the position range of the loader distributed on the wind power blade refers to the position range of the loader capable of being placed on the wind power blade, and the mass range of the recorder refers to the mass range of the loader applicable to the wind power blade and the whole testing system.
As an embodiment of the invention, blade parameters of the wind power blade are obtained, and the position range of the wind power blade where the loader can be placed is obtained by analyzing the blade parameters, generally, a region which is easy to fatigue on the wind power blade is removed, and the rest other parts are the position range of the loader; and then selecting a proper loader mass range according to the field experimental equipment, wherein the proper loader mass range is used as an initial range value of the position range and the mass range in the loader layout parameters.
According to the embodiment of the invention, the blade parameters of the wind power blade are analyzed, the position range of the loader distributed on the wind power blade and the initial mass range of the loader are determined by combining specific experimental equipment, and the configuration parameters of the subsequent loader are determined.
And generating a parameter matrix of the loader by adopting an exhaustion method according to the position range and the quality range by adopting a preset step length, wherein the parameter matrix at least comprises a position matrix and a quality matrix.
In the embodiment of the present invention, the parameter matrix of the loader is a matrix formed by layout parameters of the loader on the wind power blade, which facilitates calculation of the layout parameters of the loader subsequently, and the position matrix is a matrix formed by position parameters of the loader distributed on the wind power blade, which is generally preset by a technician, and the position parameters are selected by a certain step length within the position range determined in step S102 to form a position matrix; the quality matrix is a matrix formed by parameters in the quality range determined in step S102, and similarly, the quality parameters are selected in a certain step size in the quality range to form a position matrix.
According to the embodiment of the invention, the parameters in the position range and the quality range are selected to form the matrix, so that the completeness of the blade parameter test can be ensured, and the calculation of the subsequent loader layout parameters is facilitated.
And calculating modal parameters of a mode formed by the wind power blade and the loader by taking the parameter matrix as a variable, and taking parameters in the parameter matrix corresponding to the modal parameters meeting preset constraint conditions as loader layout parameters, wherein the loader layout parameters at least comprise the position of the loader on the wind power blade and the mass of the loader.
In the embodiment of the invention, a mode of the wind power blade with the loader is constructed, the mode parameters of the modes of the wind power blade and the loader are calculated by taking the parameter matrix determined in the step S104 as a variable, wherein the mode parameters at least comprise parameters such as fixed frequency, bending moment, exciting force and the like, the parameters are compared with the preset constraint condition, and the parameter group in the parameter matrix corresponding to the mode parameter meeting the preset constraint condition is taken as the layout parameter of the loader. In general, there is more than one set of loader layout parameters that meet the requirements, and the technician can further optimize or select the layout parameters according to the actual situation.
According to the embodiment of the invention, through setting the constraint conditions, the proper loader distribution position and the proper loader quality are selected from the position range and the quality range to be combined, so that the loader layout parameters are obtained.
The embodiment of the invention provides a loader layout method, which comprises the steps of firstly determining the position range of a loader on a wind blade and the mass range of the loader through modal analysis of the wind blade and a loader system, and then further limiting the position of the loader on the wind blade and the mass of the loader through constraints of conditions such as bending moment, fixed frequency, exciting force and the like to obtain optimized loader layout parameters and ensure the accuracy of wind blade test data.
The embodiment of the present invention further provides a computer device, which includes a memory and a processor, where the memory stores a computer program, and when the computer program is executed by the processor, the processor executes the steps of the loader layout method according to any one of the above embodiments.
FIG. 6 is a diagram illustrating an internal structure of a computer device in one embodiment. As shown in fig. 6, the computer apparatus includes a processor, a memory, a network interface, an input device, and a display screen connected through a system bus. Wherein the memory includes a non-volatile storage medium and an internal memory. The non-volatile storage medium of the computer device stores an operating system and may also store a computer program that, when executed by the processor, causes the processor to implement the loader layout method. The internal memory may also have a computer program stored therein, which when executed by the processor, causes the processor to perform the loader placement method. The display screen of the computer equipment can be a liquid crystal display screen or an electronic ink display screen, and the input device of the computer equipment can be a touch layer covered on the display screen, a key, a track ball or a touch pad arranged on the shell of the computer equipment, an external keyboard, a touch pad or a mouse and the like.
Those skilled in the art will appreciate that the architecture shown in fig. 6 is merely a block diagram of some of the structures associated with the disclosed aspects and is not intended to limit the computing devices to which the disclosed aspects apply, as particular computing devices may include more or less components than those shown, or may combine certain components, or have a different arrangement of components.
In one embodiment, a computer device is proposed, the computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, the processor implementing the following steps when executing the computer program:
blade parameters of the wind power blade are obtained, and the position range of the loader distributed on the wind power blade and the mass range of the loader are determined according to the blade parameters;
generating a parameter matrix of the loader by adopting an exhaustion method according to the position range and the quality range by adopting a preset step length, wherein the parameter matrix at least comprises a position matrix and a quality matrix;
and calculating modal parameters of a mode formed by the wind power blade and the loader by taking the parameter matrix as a variable, and taking parameters in the parameter matrix corresponding to the modal parameters meeting preset constraint conditions as loader layout parameters, wherein the loader layout parameters at least comprise the position of the loader on the wind power blade and the mass of the loader.
In one embodiment, a computer readable storage medium is provided, having a computer program stored thereon, which, when executed by a processor, causes the processor to perform the steps of: blade parameters of the wind power blade are obtained, and the position range of the loader distributed on the wind power blade and the mass range of the loader are determined according to the blade parameters;
generating a parameter matrix of the loader by adopting an exhaustion method according to the position range and the quality range by adopting a preset step length, wherein the parameter matrix at least comprises a position matrix and a quality matrix;
and calculating modal parameters of a mode formed by the wind power blade and the loader by taking the parameter matrix as a variable, and taking parameters in the parameter matrix corresponding to the modal parameters meeting preset constraint conditions as loader layout parameters, wherein the loader layout parameters at least comprise the position of the loader on the wind power blade and the mass of the loader.
It should be understood that, although the steps in the flowcharts of the embodiments of the present invention are shown in sequence as indicated by the arrows, the steps are not necessarily performed in sequence as indicated by the arrows. The steps are not limited to being performed in the exact order illustrated and, unless explicitly stated herein, may be performed in other orders. Moreover, at least a portion of steps in various embodiments may include multiple sub-steps or multiple stages that are not necessarily performed at the same time, but may be performed at different times, and the order of performance of the sub-steps or stages is not necessarily sequential, but may be performed alternately or alternatingly with other steps or at least a portion of sub-steps or stages of other steps.
It will be understood by those skilled in the art that all or part of the processes of the methods of the embodiments described above can be implemented by a computer program, which can be stored in a non-volatile computer-readable storage medium, and can include the processes of the embodiments of the methods described above when the program is executed. Any reference to memory, storage, database, or other medium used in the embodiments provided herein may include non-volatile and/or volatile memory, among others. Non-volatile memory can include read-only memory (ROM), programmable ROM (PROM), electrically Programmable ROM (EPROM), electrically Erasable Programmable ROM (EEPROM), or flash memory. Volatile memory can include Random Access Memory (RAM) or external cache memory. By way of illustration and not limitation, RAM is available in a variety of forms such as Static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), double Data Rate SDRAM (DDRSDRAM), enhanced SDRAM (ESDRAM), synchronous Link DRAM (SLDRAM), rambus (Rambus) direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM).
The technical features of the embodiments described above may be arbitrarily combined, and for the sake of brevity, all possible combinations of the technical features in the embodiments described above are not described, but should be considered as being within the scope of the present specification as long as there is no contradiction between the combinations of the technical features.
The above-mentioned embodiments only express several embodiments of the present invention, and the description thereof is specific and detailed, but not to be understood as limiting the scope of the present invention. It should be noted that, for a person skilled in the art, several variations and modifications can be made without departing from the inventive concept, which falls within the scope of the present invention. Therefore, the protection scope of the present patent shall be subject to the appended claims.
The above description is only for the purpose of illustrating the preferred embodiments of the present invention and is not to be construed as limiting the invention, and any modifications, equivalents and improvements made within the spirit and principle of the present invention are intended to be included within the scope of the present invention.
Claims (10)
1. A loader layout method, characterized in that the loader layout method comprises:
blade parameters of the wind power blade are obtained, and the position range of the loader distributed on the wind power blade and the mass range of the loader are determined according to the blade parameters;
generating a parameter matrix of the loader by adopting an exhaustion method according to the position range and the quality range by adopting a preset step length, wherein the parameter matrix at least comprises a position matrix and a quality matrix;
and calculating modal parameters of a mode formed by the wind power blade and the loader by taking the parameter matrix as a variable, and taking parameters in the parameter matrix corresponding to the modal parameters meeting preset constraint conditions as loader layout parameters, wherein the loader layout parameters at least comprise the position of the loader on the wind power blade and the mass of the loader.
2. The loader layout method according to claim 1, wherein the obtaining blade parameters of the wind turbine blades and determining the position range of the loader distributed on the wind turbine blades and the mass range of the loader according to the blade parameters comprises:
acquiring a blade fatigue damage action coefficient diagram of the wind power blade, dividing a region, with a fatigue damage action coefficient larger than a preset value, of the wind power blade into a vulnerable region according to the blade fatigue damage action coefficient diagram, and taking a part, with the vulnerable region removed, of the wind power blade as a position range where the loader is distributed on the wind power blade;
and determining the mass range according to the mass which can be borne by the clamp of the test excitation equipment and the motor power.
3. The loader layout method according to claim 1, wherein the generating the parameter matrix of the loader by an exhaustive method according to the position range and the quality range with a preset step size comprises:
in the position range, segmenting the position range by taking one end close to the wind power blade as an initial position and a preset length as a step length, and taking the distance between each segmentation point and the initial position as an element in the position matrix;
and in the quality range, constructing the quality matrix by taking the lower limit value of the quality range as an initial value and taking preset quality as a step length.
4. The loader layout method according to claim 1, wherein the calculating, with the parameter matrix as a variable, modal parameters of a modality formed by the wind turbine blade and the loader, and taking parameters in the parameter matrix corresponding to the modal parameters that satisfy preset constraints as loader layout parameters, where the loader layout parameters at least include a position of the loader on the wind turbine blade and a mass of the loader, includes:
taking parameters in the quality matrix as the quality of the loader, and taking parameters in the position matrix as the position of the loader on the wind power blade to construct a plurality of different modes;
calculating modal parameters of the modal, and checking whether the modal parameters meet the preset constraint condition; the modal parameters at least comprise a bending moment value of the modal, and the preset constraint condition at least comprises that the bending moment value is within a preset range;
and when the modal parameters meet the preset constraint condition, taking the parameters in the quality matrix corresponding to the modal parameters as the quality parameters of the loader, and taking the parameters in the position matrix as the position parameters of the loader on the wind power blade.
5. The loader layout method of claim 4, further comprising:
the modal parameters further include the fixed frequency of the modal and the exciting force of a loader in the modal, and the preset constraint conditions further include that the mass of the loader is within a specified mass range, the fixed frequency of the modal is within a specified frequency range, and the exciting force of the loader is within a specified range.
6. The loader layout method of claim 4, further comprising:
and constraining all the loader layout parameters by using a quadratic programming algorithm by taking the mean square value of the modal bending moment as a constraint condition, and taking the loader layout parameter which enables the mean square value of the bending moment to be minimum as a final layout parameter of the loader.
7. The loader layout method according to claim 1, wherein the modality of the wind blades and the loader comprises at least two loaders.
8. A loader placement apparatus, comprising:
the information acquisition unit is used for acquiring blade parameters of the wind power blade;
the information processing unit is used for determining the position range of the loader distributed on the wind power blade and the mass range of the loader according to the blade parameters; generating a parameter matrix of the loader by adopting an exhaustion method according to the position range and the quality range by adopting a preset step length, wherein the parameter matrix at least comprises a position matrix and a quality matrix; calculating modal parameters of a mode formed by the wind power blade and the loader by taking the parameter matrix as a variable, and taking parameters in the parameter matrix corresponding to the modal parameters meeting preset constraint conditions as loader layout parameters, wherein the loader layout parameters at least comprise the position of the loader on the wind power blade and the mass of the loader;
and the information output unit is used for outputting the loader layout parameters.
9. A computer device comprising a memory and a processor, the memory having stored therein a computer program that, when executed by the processor, causes the processor to perform the steps of the loader layout method of any of claims 1 to 7.
10. A computer-readable storage medium, having stored thereon a computer program which, when executed by a processor, causes the processor to carry out the steps of the loader layout method of any of claims 1 to 7.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910502187.8A CN110232238B (en) | 2019-06-11 | 2019-06-11 | Loader layout method and device, computer equipment and storage medium |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910502187.8A CN110232238B (en) | 2019-06-11 | 2019-06-11 | Loader layout method and device, computer equipment and storage medium |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110232238A CN110232238A (en) | 2019-09-13 |
CN110232238B true CN110232238B (en) | 2023-02-07 |
Family
ID=67858893
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910502187.8A Active CN110232238B (en) | 2019-06-11 | 2019-06-11 | Loader layout method and device, computer equipment and storage medium |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110232238B (en) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007147634A (en) * | 2006-12-22 | 2007-06-14 | Rikogaku Shinkokai | Method and device for vibration analysis and computer-readable recording medium |
CN102607831A (en) * | 2012-02-25 | 2012-07-25 | 中国科学院工程热物理研究所 | Evaluation method for fatigue damage and service life of horizontal axis wind turbine blade |
CN105508147A (en) * | 2015-12-28 | 2016-04-20 | 山东理工大学 | Bending moment matching method for single-point fatigue loading test of wind blade |
CN107038311A (en) * | 2017-04-21 | 2017-08-11 | 中国科学院工程热物理研究所 | A kind of resonant Fatigue Testing Loads collocation method of pneumatic equipment bladess |
-
2019
- 2019-06-11 CN CN201910502187.8A patent/CN110232238B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007147634A (en) * | 2006-12-22 | 2007-06-14 | Rikogaku Shinkokai | Method and device for vibration analysis and computer-readable recording medium |
CN102607831A (en) * | 2012-02-25 | 2012-07-25 | 中国科学院工程热物理研究所 | Evaluation method for fatigue damage and service life of horizontal axis wind turbine blade |
CN105508147A (en) * | 2015-12-28 | 2016-04-20 | 山东理工大学 | Bending moment matching method for single-point fatigue loading test of wind blade |
CN107038311A (en) * | 2017-04-21 | 2017-08-11 | 中国科学院工程热物理研究所 | A kind of resonant Fatigue Testing Loads collocation method of pneumatic equipment bladess |
Also Published As
Publication number | Publication date |
---|---|
CN110232238A (en) | 2019-09-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7006947B2 (en) | Method and apparatus for predicting failure in a system | |
CN103970965B (en) | Test run method for accelerated life test of gas turbine engine | |
CN110274826A (en) | A kind of hard metal material multiaxis high cycle fatigue failure prediction method based on single shaft fatigue S-N curve | |
CN108090963A (en) | A kind of numerical computation method of fibre reinforced composites thermal residual strain at low temperature | |
CN107862170B (en) | Finite element model correction method based on dynamic polycondensation | |
CN109614715A (en) | A kind of lower Field strength method and its application for considering notch effect of multiaxial loading effect | |
CN105260574A (en) | Critical plane method fatigue failure criterion-based high-cycle multi-axial fatigue life prediction method | |
WO2023093874A1 (en) | Building collapse risk assessment method and apparatus, computer device, and storage medium | |
Firat | A numerical analysis of combined bending–torsion fatigue of SAE notched shaft | |
CN101915733A (en) | Frequency change-based structural damage degree assessment method | |
CN110987676A (en) | Full-life prediction method considering crack closure effect under random multi-axis load | |
Sun et al. | A discrete spectral model for intermediate crack debonding in FRP-strengthened RC beams | |
CN111400852A (en) | Method and device for determining turbulence intensity parameters of wind power plant | |
Chen et al. | PSO-driven micromechanical identification of in-situ properties of fiber-reinforced composites | |
CN116720249A (en) | Construction method of stress field of concrete dam | |
Xing et al. | Shear strength reduction of trapezoidal corrugated steel plates with artificial corrosion pits | |
CN110232238B (en) | Loader layout method and device, computer equipment and storage medium | |
Ning et al. | Pipeline element for upheaval buckling analysis of submarine pipelines with geometric-imperfections under high temperature high pressure | |
Jagannathan et al. | Probabilistic strength based matrix crack evolution model in multi-directional composite laminates under fatigue loading | |
Liang et al. | Vibration-based prediction of residual fatigue life for composite laminates through frequency measurements | |
CN110501127A (en) | A kind of uniform beam damnification recognition method based on faulted condition inclination angle slope | |
WO2019200953A1 (en) | Method for measuring and analyzing degree of building strength deterioration | |
Li et al. | A coupled Armstrong-Frederick type plasticity correction methodology for calculating multiaxial notch stresses and strains | |
CN109670230A (en) | A kind of pipe safety appraisal procedure and system | |
Zhang et al. | Numerical identification of meso length-effect and full-field edge-effect of 3D braided composites |
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 |