CN112179480A - Method and system for identifying timing vibration parameters of harmonic superposition fitting blade tip and electronic equipment - Google Patents
Method and system for identifying timing vibration parameters of harmonic superposition fitting blade tip and electronic equipment Download PDFInfo
- Publication number
- CN112179480A CN112179480A CN202010873955.3A CN202010873955A CN112179480A CN 112179480 A CN112179480 A CN 112179480A CN 202010873955 A CN202010873955 A CN 202010873955A CN 112179480 A CN112179480 A CN 112179480A
- Authority
- CN
- China
- Prior art keywords
- blade
- vibration
- harmonic
- displacement data
- superposition
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H17/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves, not provided for in the preceding groups
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B21/00—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
- G01B21/02—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring length, width, or thickness
Abstract
The invention discloses a method, a system and electronic equipment for identifying timing vibration parameters of a harmonic superposition fitting blade tip, wherein the method comprises the following steps: step 1, measuring to obtain displacement data of a blade top; step 2, obtaining a simplified blade vibration equation; step 3, obtaining a vibration equation set of the corrected harmonic superposition; step 4, calculating to obtain leaf top displacement data; step 5, calculating to obtain a judgment parameter based on the leaf top displacement data obtained by measurement and calculation; judging whether the judgment parameters meet preset requirements or not; step 6, if the difference is not met, calculating to obtain a coefficient correction value based on the judgment parameter, obtaining a coefficient matrix serving as an intermediate iteration value, and skipping to execute the step 4; and 7, substituting the coefficient matrix corresponding to the judgment parameter meeting the preset requirement if the judgment parameter meets the preset requirement, and obtaining a harmonic superposition vibration fitting equation of the actual blade. The invention can better solve the technical problems of frequency aliasing and insufficient dominant frequency in the actual blade vibration.
Description
Technical Field
The invention belongs to the technical field of blade tip timing vibration parameter identification, and particularly relates to a method, a system and electronic equipment for identifying blade tip timing vibration parameters through harmonic superposition fitting.
Background
The rotating blades are core components of a plurality of power systems, and during operation, the blades often bear various complex acting forces and vibrate easily, so that fatigue, breakage and other faults of the blades are caused. Although the vibration parameters of the blades can be obtained through simulation calculation in the design process, the difference between the vibration parameters and the working conditions in a complex actual environment is often larger, and a reliable basis is difficult to provide for real-time safety monitoring and fault diagnosis of the blades, so that the vibration parameters in the actual operation of the blades are very necessary to obtain.
The characteristic of simple structure of the tip timing technology as a non-contact blade vibration measurement technology has become a research hotspot. The technology only needs to install a sensor on a static shell, and calculates the amplitude and the frequency of the vibration of the blade by measuring the time of the blade passing through the sensor; due to the technical characteristics, a single sensor for the timed measurement of the blade tip only obtains one signal for a single blade in each circle, but the vibration frequency of the blade is far higher than the self rotation frequency, so that the data obtained by the timed measurement of the blade tip is undersampled data, and the vibration parameters of the blade need to be identified by a proper method.
At present, a parameter identification method for the timed and synchronous vibration of the blade tip usually adopts a frequency multiplication traversal method to process signals acquired by experiments and fit optimal parameters in an equation. In an actual blade, frequency aliasing caused by severe undersampling of a blade tip vibration sensor exists, the complexity of actual vibration exists in various frequencies, and the main vibration frequency is not particularly prominent, so that a result obtained by the method has a large error, and even obtains a wrong vibration frequency and amplitude.
Disclosure of Invention
The invention aims to provide a method, a system and electronic equipment for identifying timing vibration parameters of a harmonic wave superposition fitting blade tip, so as to solve one or more technical problems. The invention can better solve the technical problems of frequency aliasing and insufficient dominant frequency in the actual blade vibration.
In order to achieve the purpose, the invention adopts the following technical scheme:
the invention discloses a method for identifying timing vibration parameters of a harmonic superposition fitting blade tip, which comprises the following steps of:
and 7, substituting the coefficient matrix corresponding to the judgment parameter meeting the preset requirement into the harmonic wave superposition vibration equation set obtained in the step 3 to obtain a harmonic wave superposition vibration fitting equation of the actual blade.
The invention is further improved in that in step 2, the simplified expression of the blade vibration equation is as follows:
in the formula, AjThe vibration amplitude of the blade under the order of j; p is vibration frequency doubling; phi is the initial phase.
The further improvement of the present invention is that in step 3, the expression of the modified harmonic superimposed vibration equation set is:
in the formula, C is a direct current component brought by the timing measurement of the blade tip; thetaMThe installation angle of the sensor is measured for number M.
The invention is further improved in that, in step 4, the expression of the displacement data of the leaf top obtained by calculation is as follows:
the further improvement of the invention is that in step 5, the simulated annealing algorithm is adopted to calculate and obtain the judgment parameters, and the specific steps comprise:
the given equation expression is:
F(X)＝|Y-BX|，
in the formula, Y is data obtained by timing measurement of an actual blade tip;
given perturbation new solution XnewCalculating
dE＝F(Xnew)-F(Xold)；
Determination of dE<0, then X isnewMaking a current solution; if dE>0, calculating
In the formula, KT is a given value; with P as accepting the new solution XnewThe probability of (d);
wherein, for a given KT, new solutions X are not accepted for a plurality of consecutive timesnewThen, KT is reduced according to the annealing rate, and new solution iteration is continuously generated; and finally iterating until KT meets the requirement.
The further improvement of the invention is that in step 5, the requirement of the judgment coefficient by adopting the simulated annealing algorithm is as follows:
KTmin＜1×10-6。
a further improvement of the present invention is that, in step 6, the expression of the coefficient matrix as the intermediate iteration value is:
Xm＝[A1cos1 m A1sinφ1 m...AqsinφM m Cm]T。
the further improvement of the invention is that in step 7, the expression of the harmonic superposition vibration fitting equation of the actual blade is as follows:
the invention discloses a system for identifying timing vibration parameters of a harmonic superposition fitting blade tip, which comprises:
the displacement data measuring module is used for designing a rotating blade tip timing measuring system with a shell; wherein, M measuring sensors are arranged on the shell, and a synchronous sensor is arranged on the shaft; fixing the rotating speed, and carrying out tip timing measurement on each blade to obtain tip displacement data;
the system comprises an equation set acquisition and coefficient initial value assignment module, a data acquisition module and a parameter calculation module, wherein the equation set acquisition and coefficient initial value assignment module is used for constructing and acquiring a one-dimensional centralized parameter model of the blade based on a rotating blade tip timing measurement system established by the data acquisition module; keeping the front q order of the blade vibration equation to obtain a simplified blade vibration equation; substituting the installation angle of the measuring sensor to obtain a modified vibration equation set with superposed harmonic waves; giving an initial value of a coefficient matrix of a harmonic superposition vibration equation set;
the displacement data calculation module is used for calculating and obtaining the displacement data of the blade top based on the harmonic wave superposition vibration equation set and the coefficient matrix;
the iteration judging module is used for calculating and obtaining judging parameters based on the leaf top displacement data obtained by measurement and the leaf top displacement data obtained by calculation; judging whether the judgment parameters meet preset requirements, if not, skipping to the execution correction value acquisition module, and if so, skipping to the execution parameter identification module;
the correction value acquisition module is used for calculating and obtaining coefficient correction values based on the judgment parameters, obtaining a coefficient matrix serving as an intermediate iteration value, and skipping to execute the displacement data calculation module;
and the parameter identification module is used for substituting the coefficient matrix corresponding to the judgment parameter meeting the preset requirement into the harmonic superposition vibration equation set obtained by the equation set obtaining and coefficient initial value assignment module to obtain a harmonic superposition vibration fitting equation of the actual blade and finish parameter identification.
The invention discloses harmonic superposition fitting blade tip timing vibration parameter identification electronic equipment, which comprises: a memory, a processor, and a computer program stored on the memory and executable on the processor, the computer program being executed by the processor to:
and 7, substituting the coefficient matrix corresponding to the judgment parameter meeting the preset requirement into the harmonic wave superposition vibration equation set obtained in the step 3 to obtain a harmonic wave superposition vibration fitting equation of the actual blade.
Compared with the prior art, the invention has the following beneficial effects:
the invention provides a method for identifying harmonic superposition fitting blade tip timing vibration parameters, which establishes a blade vibration equation of harmonic superposition, processes blade vibration data measured by blade tip timing by adopting a fitting method, and obtains the amplitude and phase of each order vibration of a blade, namely the harmonic superposition vibration data of the blade. The blade vibration equation obtained by the method can be used for analyzing the influence of each order vibration of the blade on the overall vibration; the multi-order vibration fitting equation obtained by the method can better solve the problems of frequency aliasing and insufficient dominant frequency in actual blade vibration; the method considers the influence of each-order vibration of the blade, so that the difference of vibration characteristics of different blades in practice can be better processed, namely, the method can process the blade with prominent main-frequency vibration and the blade with aliasing of each-order vibration, and has higher universality compared with the existing single-order blade tip timing measurement parameter identification method.
The system establishes a harmonic wave superposed blade vibration equation, processes blade vibration data measured by the blade tip at regular time by adopting a fitting method, and obtains the amplitude and the phase of each order vibration of the blade, namely obtains the harmonic wave superposed vibration data of the blade; compared with the existing single-order blade tip timing measurement parameter identification, the method has higher universality.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art are briefly introduced below; it is obvious that the drawings in the following description are some embodiments of the invention, and that for a person skilled in the art, other drawings can be derived from them without inventive effort.
FIG. 1 is a schematic flow chart of a method for identifying harmonic superposition fitting blade tip timing vibration parameters according to an embodiment of the present invention;
FIG. 2 is a schematic structural diagram of a tip timing measurement system according to an embodiment of the present invention;
FIG. 3 is a schematic diagram of a one-dimensional lumped parameter model of a blade according to an embodiment of the present invention.
Detailed Description
In order to make the purpose, technical effect and technical solution of the embodiments of the present invention clearer, the following clearly and completely describes the technical solution of the embodiments of the present invention with reference to the drawings in the embodiments of the present invention; it is to be understood that the described embodiments are only some of the embodiments of the present invention. Other embodiments, which can be derived by one of ordinary skill in the art from the disclosed embodiments without inventive faculty, are intended to be within the scope of the invention.
Referring to fig. 1, a method for identifying a harmonic superposition fitting blade tip timing vibration parameter according to an embodiment of the present invention includes the following steps:
Referring to fig. 3, the vanes are resiliently connected to adjacent disc sectors by springs and dampers, the discs also having interfaces with rotating members to cause the vanes to rotate about their centroidal axes.
For the forced vibration of the actual blade, the differential equation is:
[m]y″+[c]y′+[k]y＝F(t)；
in the formula, m, c and k are respectively a mass matrix, a damping matrix and an elastic matrix, the dimension of a square matrix of the mass matrix, the damping matrix and the elastic matrix is the degree of freedom of the whole system, y is the displacement of the blade unit, and the dimensions of the matrixes are all related to the number of actual units; f (t) is the external excitation force.
For the blade single degree of freedom model shown in fig. 3, for the excitation force f (t), it can be fourier expanded from its dynamic properties as:
in the formula, omega is the fundamental frequency of the exciting force, namely the rotating speed of the corresponding blade; fiIs the ith harmonic amplitude;
is the initial phase of the ith harmonic.
For the single-degree-of-freedom system, the displacement response of the blade obtained through the transformation of the main coordinate and the modal coordinate is as follows:
in the above formula, AiThe vibration amplitude of the blade under the order of i; piFrequency doubling for vibration; phi is aiIs the initial phase.
For the exciting force on the blade, the higher the order is, the smaller the influence on the vibration is, so in the embodiment of the present invention, the exciting force of the first q order is selected, and generally, q is more than or equal to 1 and less than or equal to 20, and then the above formula can be written as:
due to the adoption of a random angle measurement sensor arrangement method. Considering factors such as installation, the number of the sensors is M, generally, M is more than or equal to 4 and less than or equal to 8, and the measured data of any blade meets the following equation set:
and 3, for the equation set which is satisfied by the tip timing data obtained in the step 2, in order to obtain each coefficient in the equation, the equation can be expanded:
order to
X＝[A1cosφ1 A1sinφ1 A2cosφ2 A2sinφ2 ... AqcosφM AqsinφM C]T，
If the angle of the first measuring sensor can be defined as 0, a coefficient matrix can be obtained
Y＝(y1 y2 … yM)T；
The equation can be expressed as Y ═ BX, since P is a natural number as the frequency multiplication of the vibration, the unknown quantities of a, phi and C are now included in X,
X＝(BTB)-1BTY；
at present, the harmonic wave superposition equation is already solved for the equation waiting coefficient, and various algorithms can be adopted for solving the part, in the embodiment of the invention, a simulated annealing algorithm can be adopted, wherein
F(X)＝|Y-BX|；
Given the perturbation, a new solution is generated, yielding F (X)new),dE＝F(Xnew)-F(X)；
Determination of dE<0, then X isnewMaking a current solution; if dE>0, calculating
KT in the above formula is given value. With P as accepting the new solution XnewThe probability of (c).
Continuing the above cycle to continuously receive no new solution X for multiple timesnew。
In the embodiment of the invention, given an annealing rate of 0.98, the KT value is reduced and the continuous circulation is carried out until the KT value is small enoughminThe value is related to the specific accuracy requirement, and is generally 1 × 10-6。
After the iteration is finished, obtaining the final Xfinal，
Then the vibration equation for this model based on harmonic superposition fitting can be expressed as:
and obtaining the vibration parameters of each blade at any moment under the model according to the equation.
The system for identifying the harmonic superposition fitting blade tip timing vibration parameters comprises the following steps:
the displacement data measuring module is used for designing a rotating blade tip timing measuring system with a shell; wherein, M measuring sensors are arranged on the shell, and a synchronous sensor is arranged on the shaft; fixing the rotating speed, and carrying out tip timing measurement on each blade to obtain tip displacement data; the structure of the rotating blade tip timing measurement system is shown in FIG. 2.
The system comprises an equation set acquisition and coefficient initial value assignment module, a data acquisition module and a parameter calculation module, wherein the equation set acquisition and coefficient initial value assignment module is used for constructing and acquiring a one-dimensional centralized parameter model of the blade based on a rotating blade tip timing measurement system established by the data acquisition module; keeping the front q order of the blade vibration equation to obtain a simplified blade vibration equation; substituting the installation angle of the measuring sensor to obtain a modified vibration equation set with superposed harmonic waves; giving an initial value of a coefficient matrix of a harmonic superposition vibration equation set;
the displacement data calculation module is used for calculating and obtaining the displacement data of the blade top based on the harmonic wave superposition vibration equation set and the coefficient matrix;
the iteration judging module is used for calculating and obtaining judging parameters based on the leaf top displacement data obtained by measurement and the leaf top displacement data obtained by calculation; judging whether the judgment parameters meet preset requirements, if not, skipping to the execution correction value acquisition module, and if so, skipping to the execution parameter identification module;
the correction value acquisition module is used for calculating and obtaining coefficient correction values based on the judgment parameters, obtaining a coefficient matrix serving as an intermediate iteration value, and skipping to execute the displacement data calculation module;
and the parameter identification module is used for substituting the coefficient matrix corresponding to the judgment parameter meeting the preset requirement into the harmonic superposition vibration equation set obtained by the equation set obtaining and coefficient initial value assignment module to obtain a harmonic superposition vibration fitting equation of the actual blade and finish parameter identification.
The harmonic superposition fitting blade tip timing vibration parameter identification electronic equipment comprises:
on the hardware level, the electronic device comprises: the processor optionally further comprises an internal bus, a network interface and a memory. The memory may include a memory, such as a high speed random access memory, and may also include a non-volatile memory, such as at least one disk memory. Of course, the electronic device may also include hardware required for other services.
The processor, network interface, memory are interconnected by an internal bus, which may be an industry standard architecture bus, a peripheral component interconnect standard bus, an extended industry standard architecture bus, or the like. The bus may be divided into an address bus, a data bus, a control bus, etc.
The memory is used for storing programs. In particular, the program may include program code, which includes computer operating instructions. The memory may include both memory and non-volatile storage and provides instructions and data to the processor.
The processor reads the corresponding computer program from the nonvolatile memory into the memory and then runs the computer program to form the electronic equipment on the logic level. When the processor executes the program, the following operations are specifically executed:
and 7, substituting the coefficient matrix corresponding to the judgment parameter meeting the preset requirement into the harmonic wave superposition vibration equation set obtained in the step 3 to obtain a harmonic wave superposition vibration fitting equation of the actual blade.
In the method for identifying the multi-harmonic fitting blade tip timing vibration parameters in the embodiment of the invention,
for the excitation force applied to the blade, the higher the order is, the smaller the influence on the vibration is, so in the embodiment of the present invention, the excitation force of the first 20 orders is selected, and the above formula can be written as follows:
due to the adoption of the random angle 4 measurement sensor arrangement method, the measurement data of any blade meets the following equation set:
and 3, for the equation set which is satisfied by the tip timing data obtained in the step 2, in order to obtain each coefficient in the equation, the equation can be expanded:
order to
X＝[A1cosφ1 A1sinφ1 A2cosφ2 A2sinφ2 ... A20cosφ20 A20sinφ20 C]T，
If the angle of the first measuring sensor can be defined as 0, a coefficient matrix can be obtained
Y＝(y1 y2 y3 y4)；
The equation can be expressed as Y ═ BX, since P is a natural number as the frequency multiplication of the vibration, so now 41 unknowns of a, phi and C are included in X,
X＝(BTB)-1BTY；
at present, the harmonic wave superposition equation is already solved for the equation waiting coefficient, and various algorithms can be adopted for solving the part, in the embodiment of the invention, a simulated annealing algorithm can be adopted, wherein
F(X)＝|Y-BX|；
Given the perturbation, a new solution is generated, yielding F (X)new),dE＝F(Xnew)-F(X)；
Determination of dE<0, then X isnewMaking a current solution; if dE>0, calculating
KT in the above formula is given value. With P as accepting the new solution XnewThe probability of (c).
Continuing the above cycle to continuously receive no new solution X for multiple timesnew。
For a given KT, new solutions X are not accepted for a plurality of consecutive timesnewThen, KT is reduced according to the annealing rate, and new solution iteration is continuously generated;
and finally iterating until KT meets the requirement.
In the embodiment of the invention, given an annealing rate of 0.98, the KT value is reduced and the continuous circulation is carried out until the KT value is small enoughminThe value is related to the specific accuracy requirement, and is generally 1 e-6.
After the iteration is finished, obtaining the final Xfinal，
Then the vibration equation for this model based on harmonic superposition fitting can be expressed as:
and obtaining the vibration parameters of each blade at any moment under the model according to the equation.
In engineering application, firstly, a one-dimensional lumped parameter model is established according to the actual number of the blades, and a vibration equation of the blades is obtained; installing reference timing sensors on the shaft, including rotating speed synchronous sensor and angle reference sensor, installing blade tip timing sensors at different positions on the casing, giving numbers and obtaining the installation angle theta of the number k measuring sensork(ii) a Running the whole system, and comparing the time interval sequence t of the signals of the plurality of timing sensors with the rotating speed synchronizing signalijRecording, and obtaining a displacement signal y of blade vibration based on the rotating speedk(ii) a The vibration displacement signals are processed, and the vibration equation of each blade is obtained by adopting the harmonic superposition equation fitting method provided by the invention, so that the vibration state of the blade can be monitored in real time.
As will be appreciated by one skilled in the art, embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
Although the present invention has been described in detail with reference to the above embodiments, those skilled in the art can make modifications and equivalents to the embodiments of the present invention without departing from the spirit and scope of the present invention, which is set forth in the claims of the present application.
Claims (10)
1. A method for identifying timing vibration parameters of a harmonic superposition fitting blade tip is characterized by comprising the following steps:
step 1, designing a rotating blade tip timing measurement system with a shell; wherein, M measuring sensors are arranged on the shell, and a synchronous sensor is arranged on the shaft; fixing the rotating speed, and carrying out tip timing measurement on each blade to obtain tip displacement data;
step 2, constructing a one-dimensional centralized parameter model of the obtained blade based on the rotating blade tip timing measurement system in the step 1; keeping the front q order of the blade vibration equation to obtain a simplified blade vibration equation;
step 3, substituting the installation angle of the measuring sensor to obtain a vibration equation set of the corrected harmonic superposition; giving an initial value of a coefficient matrix of a harmonic superposition vibration equation set;
step 4, calculating to obtain the displacement data of the blade top based on the harmonic wave superposition vibration equation set and the coefficient matrix;
step 5, calculating to obtain a judgment parameter based on the leaf top displacement data obtained by measurement in the step 1 and the leaf top displacement data obtained by calculation in the step 4; judging whether the judgment parameters meet preset requirements, if not, skipping to execute the step 6, and if so, skipping to execute the step 7;
step 6, calculating to obtain a coefficient correction value based on the judgment parameter, obtaining a coefficient matrix as a middle iteration value, and skipping to execute the step 4;
and 7, substituting the coefficient matrix corresponding to the judgment parameter meeting the preset requirement into the harmonic wave superposition vibration equation set obtained in the step 3 to obtain a harmonic wave superposition vibration fitting equation of the actual blade.
2. The method for identifying the timing vibration parameters of the blade tip through harmonic superposition fitting according to claim 1, wherein in the step 2, the expression of the simplified blade vibration equation is as follows:
in the formula, AjThe vibration amplitude of the blade under the order of j; p is vibration frequency doubling; phi is the initial phase.
3. The method for identifying the blade tip timing vibration parameters through harmonic superposition fitting according to claim 2, wherein in the step 3, the expression of the modified harmonic superposition vibration equation set is as follows:
in the formula, C is a direct current component brought by the timing measurement of the blade tip; thetaMThe installation angle of the sensor is measured for number M.
4. The method for identifying the timing vibration parameters of the blade tip through harmonic superposition fitting according to claim 3, wherein in the step 4, the expression of the blade tip displacement data obtained through calculation is as follows:
5. the method for identifying the timing vibration parameters of the harmonic superposition fitting blade tip according to claim 4, wherein in the step 5, the judgment parameters are obtained by calculating through a simulated annealing algorithm, and the specific steps comprise:
the given equation expression is:
F(X)＝|Y-BX|，
in the formula, Y is data obtained by timing measurement of an actual blade tip;
given perturbation new solution XnewCalculating
dE＝F(Xnew)-F(Xold)；
If dE is less than 0, X is determinednewMaking a current solution; if dE > 0, calculate
In the formula, KT is a given value; with P as accepting the new solution XnewThe probability of (d);
for a given KT, new solutions X are not accepted for a predetermined number of consecutive timesnewThen, KT is reduced according to the annealing rate, and new solution iteration is continuously generated; and finally iterating until KT meets the requirement.
6. The method for identifying the timing vibration parameters of the harmonic superposition fitting blade tip according to claim 5, wherein in the step 5, the requirement of the judgment coefficient by adopting the simulated annealing algorithm is as follows:
KTmin＜1×10-6。
7. the method for identifying the harmonic superposition fitting blade tip timing vibration parameters according to claim 4, wherein in the step 6, the expression of the coefficient matrix as the intermediate iteration value is as follows:
Xm＝[A1cosφ1 m A1sinφ1 m … AqsinφM m Cm]T。
8. the method for identifying the harmonic superposition fitting blade tip timing vibration parameters according to claim 7, wherein in the step 7, the expression of the harmonic superposition vibration fitting equation of the actual blade is as follows:
9. a system for identifying timing vibration parameters of a harmonic superposition fitting blade tip is characterized by comprising:
the displacement data measuring module is used for designing a rotating blade tip timing measuring system with a shell; wherein, M measuring sensors are arranged on the shell, and a synchronous sensor is arranged on the shaft; fixing the rotating speed, and carrying out tip timing measurement on each blade to obtain tip displacement data;
the system comprises an equation set acquisition and coefficient initial value assignment module, a data acquisition module and a parameter calculation module, wherein the equation set acquisition and coefficient initial value assignment module is used for constructing a one-dimensional centralized parameter model of the blade based on a rotating blade tip timing measurement system designed by the data acquisition module; keeping the front q order of the blade vibration equation to obtain a simplified blade vibration equation; substituting the installation angle of the measuring sensor to obtain a modified vibration equation set with superposed harmonic waves; giving an initial value of a coefficient matrix of a harmonic superposition vibration equation set;
the displacement data calculation module is used for calculating and obtaining the displacement data of the blade top based on the harmonic wave superposition vibration equation set and the coefficient matrix;
the iteration judging module is used for calculating and obtaining judging parameters based on the leaf top displacement data obtained by measurement and the leaf top displacement data obtained by calculation; judging whether the judgment parameters meet preset requirements, if not, skipping to the execution correction value acquisition module, and if so, skipping to the execution parameter identification module;
the correction value acquisition module is used for calculating and obtaining coefficient correction values based on the judgment parameters, obtaining a coefficient matrix serving as an intermediate iteration value, and skipping to execute the displacement data calculation module;
and the parameter identification module is used for substituting the coefficient matrix corresponding to the judgment parameter meeting the preset requirement into the harmonic superposition vibration equation set obtained by the equation set obtaining and coefficient initial value assignment module to obtain a harmonic superposition vibration fitting equation of the actual blade and finish parameter identification.
10. A harmonic superposition fitting blade tip timing vibration parameter identification electronic device is characterized by comprising: a memory, a processor, and a computer program stored on the memory and executable on the processor, the computer program being executed by the processor to:
step 1, designing a rotating blade tip timing measurement system with a shell; wherein, M measuring sensors are arranged on the shell, and a synchronous sensor is arranged on the shaft; fixing the rotating speed, and carrying out tip timing measurement on each blade to obtain tip displacement data;
step 2, constructing a one-dimensional centralized parameter model of the obtained blade based on the rotating blade tip timing measurement system designed in the step 1; keeping the front q order of the blade vibration equation to obtain a simplified blade vibration equation;
step 3, substituting the installation angle of the measuring sensor to obtain a vibration equation set of the corrected harmonic superposition; giving an initial value of a coefficient matrix of a harmonic superposition vibration equation set;
step 4, calculating to obtain the displacement data of the blade top based on the harmonic wave superposition vibration equation set and the coefficient matrix;
step 5, calculating to obtain a judgment parameter based on the leaf top displacement data obtained by measurement in the step 1 and the leaf top displacement data obtained by calculation in the step 4; judging whether the judgment parameters meet preset requirements, if not, skipping to execute the step 6, and if so, skipping to execute the step 7;
step 6, calculating to obtain a coefficient correction value based on the judgment parameter, obtaining a coefficient matrix as a middle iteration value, and skipping to execute the step 4;
and 7, substituting the coefficient matrix corresponding to the judgment parameter meeting the preset requirement into the harmonic wave superposition vibration equation set obtained in the step 3 to obtain a harmonic wave superposition vibration fitting equation of the actual blade.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010873955.3A CN112179480B (en) | 2020-08-26 | 2020-08-26 | Method and system for identifying timing vibration parameters of harmonic superposition fitting blade tip and electronic equipment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010873955.3A CN112179480B (en) | 2020-08-26 | 2020-08-26 | Method and system for identifying timing vibration parameters of harmonic superposition fitting blade tip and electronic equipment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112179480A true CN112179480A (en) | 2021-01-05 |
| CN112179480B CN112179480B (en) | 2021-11-19 |
Family
ID=73925811
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010873955.3A Active CN112179480B (en) | 2020-08-26 | 2020-08-26 | Method and system for identifying timing vibration parameters of harmonic superposition fitting blade tip and electronic equipment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112179480B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114623994A (en) * | 2022-01-13 | 2022-06-14 | 天津大学 | Blade multi-modal coupling synchronous vibration parameter identification method based on blade tip timing |
| CN114812794A (en) * | 2022-01-12 | 2022-07-29 | 天津大学 | Method for correcting blade tip timing parameter identification error caused by frequency sweep rate |
Citations (18)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0757786A1 (en) * | 1994-04-28 | 1997-02-12 | Electric Power Research Institute | Method and apparatus for isolating and identifying periodic doppler signals |
| EP2237003B1 (en) * | 2008-06-18 | 2014-07-09 | Rolls-Royce plc | Method for analysing vibration in rotor blades |
| CN104697623A (en) * | 2014-10-24 | 2015-06-10 | 南京航空航天大学 | Method of recognizing asynchronous vibration parameters of blades under variable speed excitation |
| CN105300513A (en) * | 2015-11-20 | 2016-02-03 | 天津大学 | Method and device for identifying resonance frequency multiplication number of blade vibration of multiple tip timing sensors |
| CN105352713A (en) * | 2015-11-20 | 2016-02-24 | 天津大学 | Method used for realizing blade fatigue crack detection through static blade displacement detection |
| CN105424160A (en) * | 2015-11-20 | 2016-03-23 | 天津大学 | Method for realizing blade synchronous vibration parameter identification |
| KR20160123606A (en) * | 2015-04-16 | 2016-10-26 | 한양대학교 산학협력단 | System and method for analysing vibration characteristic |
| JP2017156150A (en) * | 2016-02-29 | 2017-09-07 | 三菱重工コンプレッサ株式会社 | Vibration measurement device, vibration measurement system, and vibration measurement method |
| BRPI0502176B1 (en) * | 2004-06-16 | 2017-11-28 | General Electric Company | METHOD AND DEVICE FOR THE INSPECTION OF CONTOURS IN TOUGH SPRINKLER RABBIT |
| CN207420937U (en) * | 2017-05-05 | 2018-05-29 | 河南理工大学 | A kind of axial fan blade embedded radio vibration detecting device |
| CN109313069A (en) * | 2016-04-26 | 2019-02-05 | 比勒陀利亚大学 | A kind of method and system using Tip-Timing (BTT) monitoring turbine rotor blade |
| CN109657397A (en) * | 2018-12-29 | 2019-04-19 | 山东大学 | Turbine blade-rotor system stability prediction technique based on frequency response function |
| CN109883380A (en) * | 2019-03-22 | 2019-06-14 | 西安交通大学 | A kind of rotor blade displacement field measurement method and its system based on blade tip-timing |
| JP2019100210A (en) * | 2017-11-29 | 2019-06-24 | 三菱日立パワーシステムズ株式会社 | Rotor blade monitoring system, rotor blade monitoring device, rotor blade monitoring method, program |
| CN109974849A (en) * | 2019-04-03 | 2019-07-05 | 上海交通大学 | Without under reference signal based on the blade vibration on-line monitoring method of Tip-Timing technology |
| CN110567574A (en) * | 2019-08-02 | 2019-12-13 | 西安交通大学 | Method and system for identifying timing vibration parameters of blade end of rotating blade |
| CN110686764A (en) * | 2019-09-17 | 2020-01-14 | 天津大学 | Method for measuring asynchronous vibration frequency of constant-speed blade based on full-phase difference principle |
| CN111239432A (en) * | 2020-02-10 | 2020-06-05 | 西安交通大学 | Blade arrival time measuring method based on blade end timing and acquisition system thereof |
-
2020
- 2020-08-26 CN CN202010873955.3A patent/CN112179480B/en active Active
Patent Citations (18)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0757786A1 (en) * | 1994-04-28 | 1997-02-12 | Electric Power Research Institute | Method and apparatus for isolating and identifying periodic doppler signals |
| BRPI0502176B1 (en) * | 2004-06-16 | 2017-11-28 | General Electric Company | METHOD AND DEVICE FOR THE INSPECTION OF CONTOURS IN TOUGH SPRINKLER RABBIT |
| EP2237003B1 (en) * | 2008-06-18 | 2014-07-09 | Rolls-Royce plc | Method for analysing vibration in rotor blades |
| CN104697623A (en) * | 2014-10-24 | 2015-06-10 | 南京航空航天大学 | Method of recognizing asynchronous vibration parameters of blades under variable speed excitation |
| KR20160123606A (en) * | 2015-04-16 | 2016-10-26 | 한양대학교 산학협력단 | System and method for analysing vibration characteristic |
| CN105352713A (en) * | 2015-11-20 | 2016-02-24 | 天津大学 | Method used for realizing blade fatigue crack detection through static blade displacement detection |
| CN105424160A (en) * | 2015-11-20 | 2016-03-23 | 天津大学 | Method for realizing blade synchronous vibration parameter identification |
| CN105300513A (en) * | 2015-11-20 | 2016-02-03 | 天津大学 | Method and device for identifying resonance frequency multiplication number of blade vibration of multiple tip timing sensors |
| JP2017156150A (en) * | 2016-02-29 | 2017-09-07 | 三菱重工コンプレッサ株式会社 | Vibration measurement device, vibration measurement system, and vibration measurement method |
| CN109313069A (en) * | 2016-04-26 | 2019-02-05 | 比勒陀利亚大学 | A kind of method and system using Tip-Timing (BTT) monitoring turbine rotor blade |
| CN207420937U (en) * | 2017-05-05 | 2018-05-29 | 河南理工大学 | A kind of axial fan blade embedded radio vibration detecting device |
| JP2019100210A (en) * | 2017-11-29 | 2019-06-24 | 三菱日立パワーシステムズ株式会社 | Rotor blade monitoring system, rotor blade monitoring device, rotor blade monitoring method, program |
| CN109657397A (en) * | 2018-12-29 | 2019-04-19 | 山东大学 | Turbine blade-rotor system stability prediction technique based on frequency response function |
| CN109883380A (en) * | 2019-03-22 | 2019-06-14 | 西安交通大学 | A kind of rotor blade displacement field measurement method and its system based on blade tip-timing |
| CN109974849A (en) * | 2019-04-03 | 2019-07-05 | 上海交通大学 | Without under reference signal based on the blade vibration on-line monitoring method of Tip-Timing technology |
| CN110567574A (en) * | 2019-08-02 | 2019-12-13 | 西安交通大学 | Method and system for identifying timing vibration parameters of blade end of rotating blade |
| CN110686764A (en) * | 2019-09-17 | 2020-01-14 | 天津大学 | Method for measuring asynchronous vibration frequency of constant-speed blade based on full-phase difference principle |
| CN111239432A (en) * | 2020-02-10 | 2020-06-05 | 西安交通大学 | Blade arrival time measuring method based on blade end timing and acquisition system thereof |
Non-Patent Citations (6)
| Title |
|---|
| FAMING YANG: "《An Improved Algorithm for Accelerating the Sparse Reconstruction Speed of Blade Tip-Timing Signal》", 《2019 IEEE 8TH JOINT INTERNATIONAL INFORMATION TECHNOLOGY AND ARTIFICIAL INTELLIGENCE CONFERENCE (ITAIC)》 * |
| JUN LIN: "《Sparse reconstruction of blade tip-timing signals for multi-mode blade vibration monitoring》", 《MECHANICAL SYSTEMS AND SIGNAL PROCESSING》 * |
| WEIMIN WANG: "《A novel none once per revolution blade tip timing based blade vibration parameters identification method》", 《CHINESE JOURNAL OF AERONAUTICS》 * |
| 李孟麟: "《基于叶尖定时的旋转机械叶片振动频率辨识ESPRIT 方法》", 《振动与冲击》 * |
| 许艳芝: "《非接触叶尖振动参数辨识方法验证试验》", 《机械研究与应用》 * |
| 雷杰: "《基于叶尖定时原理的转子叶片同步振动算法研究》", 《机械研究与应用》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114812794A (en) * | 2022-01-12 | 2022-07-29 | 天津大学 | Method for correcting blade tip timing parameter identification error caused by frequency sweep rate |
| CN114812794B (en) * | 2022-01-12 | 2024-03-22 | 天津大学 | Blade tip timing parameter identification error correction method caused by sweep rate |
| CN114623994A (en) * | 2022-01-13 | 2022-06-14 | 天津大学 | Blade multi-modal coupling synchronous vibration parameter identification method based on blade tip timing |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112179480B (en) | 2021-11-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113221280B (en) | Rolling bearing modeling and model updating method and system based on digital twinning | |
| CN112179480B (en) | Method and system for identifying timing vibration parameters of harmonic superposition fitting blade tip and electronic equipment | |
| Vargiu et al. | A reduced order model based on sector mistuning for the dynamic analysis of mistuned bladed disks | |
| US5633800A (en) | Integrated model-based reasoning/expert system diagnosis for rotating machinery | |
| CN105008887B (en) | Use the turbine blade analysis of Fatigue-life of non-cpntact measurement and dynamic response reconfiguration technique | |
| CN109883720B (en) | Arrangement method of blade tip timing sensor for recognizing multi-modal vibration of blade | |
| CN110851963A (en) | Casing circumferential arrangement method of blade end timing sensor | |
| JP5561787B2 (en) | Method and apparatus for correlating a simulation model with a physical device based on a correlation metric | |
| CN114662425B (en) | Water turbine start-stop working condition flow field simulation prediction method and system | |
| CN106289776B (en) | Flexible magnetic suspension bearing rotor rigidity damping identification method | |
| US20150073730A1 (en) | Mechanical strain gauge simulation | |
| Custer et al. | Unsteady simulation of a 1.5 stage turbine using an implicitly coupled nonlinear harmonic balance method | |
| JP6906985B2 (en) | Vibration diagnosis system, vibration diagnosis method and parameter setting method | |
| CN108052963A (en) | The data screening method, apparatus and wind power generating set of wind power prediction modeling | |
| US8839664B2 (en) | Detection and classification of failures of power generating equipment during transient conditions | |
| Fan et al. | An improved mistuning identification and dynamic strain prediction method for rotating blades with application of blade tip timing technology | |
| CN114692388A (en) | Simulation system, method and device for thermal power plant and storage medium | |
| CN110207967B (en) | State identification method and system based on wavelet packet energy characteristics and cross correlation | |
| Dedoussis et al. | Numerical simulation of blade fault signatures from unsteady wall pressure signals | |
| Juul et al. | OMA identification of a wind turbine rotor hub for monitoring purpose | |
| CN117454721B (en) | Wind power plant wake superposition effect evaluation method and medium based on digital simulation experiment | |
| Han et al. | A temperature sensor failure diagnosis approach in frequency | |
| Manzato et al. | Order Based Modal Analysis Versus Standard Techniques to Extract Modal Parameters of Operational Wind Turbine Gearboxes | |
| Kosorus et al. | On the identification of frequencies and damping ratios for structural health monitoring using autoregressive models | |
| CN112539829B (en) | Arrangement method and system of self-adaptive blade tip timing sensors |
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 |