CN116877212A - Method and system for timing monitoring of blade tips of rotational speed-free reference turbine blade - Google Patents

Method and system for timing monitoring of blade tips of rotational speed-free reference turbine blade Download PDF

Info

Publication number
CN116877212A
CN116877212A CN202310835523.7A CN202310835523A CN116877212A CN 116877212 A CN116877212 A CN 116877212A CN 202310835523 A CN202310835523 A CN 202310835523A CN 116877212 A CN116877212 A CN 116877212A
Authority
CN
China
Prior art keywords
blade
vibration
tip timing
tip
sensors
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202310835523.7A
Other languages
Chinese (zh)
Inventor
王亚南
乔百杰
李贵龙
朱昱达
杨志勃
陈雪峰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xian Jiaotong University
Original Assignee
Xian Jiaotong University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xian Jiaotong University filed Critical Xian Jiaotong University
Priority to CN202310835523.7A priority Critical patent/CN116877212A/en
Publication of CN116877212A publication Critical patent/CN116877212A/en
Pending legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F01—MACHINES OR ENGINES IN GENERAL; ENGINE PLANTS IN GENERAL; STEAM ENGINES
    • F01D—NON-POSITIVE DISPLACEMENT MACHINES OR ENGINES, e.g. STEAM TURBINES
    • F01D21/00—Shutting-down of machines or engines, e.g. in emergency; Regulating, controlling, or safety means not otherwise provided for
    • F01D21/003—Arrangements for testing or measuring
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F01—MACHINES OR ENGINES IN GENERAL; ENGINE PLANTS IN GENERAL; STEAM ENGINES
    • F01D—NON-POSITIVE DISPLACEMENT MACHINES OR ENGINES, e.g. STEAM TURBINES
    • F01D5/00—Blades; Blade-carrying members; Heating, heat-insulating, cooling or antivibration means on the blades or the members
    • F01D5/12—Blades
    • F01D5/26—Antivibration means not restricted to blade form or construction or to blade-to-blade connections or to the use of particular materials
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • 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)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mechanical Engineering (AREA)
  • Data Mining & Analysis (AREA)
  • Evolutionary Computation (AREA)
  • Operations Research (AREA)
  • Algebra (AREA)
  • Computer Hardware Design (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Geometry (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

The method comprises the steps of establishing a finite element model of a turbine blade, carrying out modal analysis to obtain a vibration mode of blade vibration, determining the installation number and angle of sensors, collecting the actual arrival time of each blade at the sensors by using a blade tip timing sensor, calculating the vibration acceleration of the blade according to the arrival time collected by the sensors, calculating the vibration acceleration parameters of the blade by using the vibration acceleration of the blade and the angle layout of the sensors, and calculating the vibration displacement parameters by using the vibration acceleration parameters.

Description

Method and system for timing monitoring of blade tips of rotational speed-free reference turbine blade
Technical Field
The application relates to the technical field of non-contact measurement of turbine blades, in particular to a method and a system for timing monitoring of blade tips of a rotational speed-free reference turbine blade.
Background
Steam turbines are still important as one of the main devices in the traditional energy field. Steam turbine blades are an important component of a steam turbine, and the state of the steam turbine blades directly affects the safety and economy of the steam turbine. However, the conventional turbine blade detection method generally adopts offline detection, and needs to stop and overhaul the turbine, which consumes time and has high labor cost. In order to improve the operation efficiency and reliability of a steam turbine, a method and a system capable of monitoring the state of the blade end of a blade of the steam turbine on line in real time need to be developed. Since the working environment of the steam turbine is extremely bad and usually operates in an environment which is not designed to be standard, the vibration signals suffered by the blades contain a great deal of noise and interference, so that the identification and judgment of the states of the blades become very complex and difficult. Thus, for such a complex system to operate at high speed, blade monitoring needs to be highly real-time and scalable. In this regard, tip timing (BTT) technology may meet this requirement. The tip timing technology is a measurement technology which is rising in recent years, adopts non-contact measurement, utilizes a tip sensor to determine the arrival time of a blade, utilizes a rotating speed sensor to determine ideal time, and obtains the vibration displacement of the blade through the time difference of the tip sensor and the rotating speed sensor. However, this approach may result in large displacement calculation errors due to calculation bias of the ideal arrival time. In order to calculate the blade vibration displacement more accurately, it is necessary to eliminate the error of the ideal arrival time. By calculating the tip vibration velocity to calculate the blade vibration displacement, deviations in the ideal arrival time can be eliminated, thereby obtaining more accurate results. However, calculating the tip vibration speed to calculate the high-frequency vibration is not ideal in effect, and a new tip timing method is required in order to realize accurate identification of the high-frequency vibration.
The above information disclosed in the background section is only for enhancement of understanding of the background of the application and therefore may contain information that does not form the prior art that is already known to a person of ordinary skill in the art.
Disclosure of Invention
Aiming at the defects in the prior art, the application aims to provide a method and a system for identifying the leaf-end timing vibration parameters of a full frequency domain, which omits a rotating speed sensor, eliminates the influence caused by inaccurate ideal arrival time and calculates displacement more accurately; the installation requirement of the sensor is reduced, and the engineering application capacity of the blade tip timing is greatly improved; the vibration acceleration signal has the characteristics of high precision, high sensitivity, lower noise and the like, so that the blade tip timing acceleration algorithm can realize high-precision blade vibration displacement calculation; the identification effect is better in high frequency.
In order to achieve the above object, the present application provides the following technical solutions:
the application relates to a method for monitoring the blade tip timing of a rotating speed-free reference turbine blade, which comprises the following steps:
step S1, a finite element model of a turbine blade is established, modal analysis is carried out on the finite element model, and the vibration mode of the turbine at a preset working speed is calculated to obtain the blade vibration mode parameter f= [ f ] 1 f 2 …f m ]；
Step S2, circumferentially installing k groups of blade tip timing sensors on a rotor blade casing, and acquiring a time sequence of the rotor blade reaching the blade tip timing sensors;
step S3, calculating to obtain the vibration acceleration of the blade by using the time sequence and the installation angle of the blade tip timing sensor;
and a fourth step S4, obtaining the vibration frequency, amplitude and phase of the turbine blade at a constant rotating speed through a circumferential Fourier algorithm based on the blade tip vibration acceleration.
In the method, in a first step S1, the number m of the blade modes of the steam turbine is obtained through mode analysis, and the number k of the blade tip timing sensor groups is determined according to the number m of the blade modes.
In the method, the relationship between the number k of the blade tip timing sensor groups and the number m of the blade modes is k=3m.
In the method, in the first step S1, the number of the tip timing sensors in each group is three, and the angular interval of the tip timing sensors in each group is less than 5 degrees.
In the method, the preset working rotating speed of the steam turbine is 3000r/min.
In the method, an angle matrix of an nth-circle blade tip timing sensor is constructed according to the vibration mode parameters of the blade and the time sequence
wherein fm The vibration frequency of m-order modes, m is the mode number, t k,n Average of blade arrival times acquired for the kth group of tip timing sensors at the nth turn, where For the time acquired in the nth turn for the smaller angle of the kth group tip timing sensor, +.>Time acquired in the nth turn for the middle of the kth group tip timing sensor angle, +.>The time acquired at the nth turn is the one with the greatest angle of the kth group of tip timing sensors.
In the third step S3, the vibration acceleration of the blade tip in the nth turn is calculated by using the angle between the kth group of tip timing sensors and the acquired arrival time wherein ,/>R is the radius of the steam turbine, theta k,1 、θ k,2 、θ k,3 Respectively the angles of the kth group of sensors after the sensor is sequenced according to the angles, and the rotating speed f of the steam turbine Ω Constant.
In the method, in the fourth step S4, tip timing vibration parameters are obtained through a circumferential fourier algorithm based on tip vibration acceleration, where the tip timing vibration parameters include:
a n ＝Φ n S n, wherein wherein />Vibration acceleration values obtained for the kth group of blade tip timing sensors at the nth turn; s is S n ＝[i 1 j 1 i 2 j 2 …i m j m ] T ，S n The vibration parameter matrix is the nth circle;
vibration acceleration amplitude values under different blade vibration modesVibration velocity phaseVibration speed amplitude in different blade vibration modes>Vibration velocity phaseVibration displacement amplitude in different blade vibration modes>Vibration velocity phase
A system for implementing the rotational speed reference-free turbine blade tip timing monitoring method includes,
the blade modeling and finite element analysis module is configured to establish a finite element model of a turbine blade and perform modal analysis to obtain a blade vibration modal parameter f= [ f ] 1 f 2 …f m ]；
A sensor layout module connected to the blade modeling and finite element analysis module, the sensor layout module determining a number of tip timing sensors based on the vibration mode parameters and constructing a sensor angle matrix based on the blade vibration mode parameters and the arrival time
A time of arrival acquisition module configured to acquire a time of arrival sequence of the blade through the tip timing sensor, the individual blade time of arrival sequences acquired by the set of sensors being
The vibration acceleration calculation module is connected with the arrival time acquisition module and calculates the vibration speed by using the arrival time, and the vibration speed sequence of the single blade is a= [ a ] 1 a 2 …a n ]；
The vibration parameter acquisition module is connected with the sensor layout module and the vibration acceleration calculation module, and the vibration parameter acquisition module solves the vibration parameter S by using a circumferential Fourier according to the sensor angle matrix and the vibration speed n ＝[i 1 j 1 i 2 j 2 …i m j m ] T 。
Advantageous effects
The application has no rotation speed sensor, simplifies the installation structure and improves the application space in a large-scale device. The rotational speed-free design can eliminate measurement errors caused by inaccurate ideal arrival time, and the vibration displacement domain of the blade is transferred to the acceleration domain for analysis, so that the identification accuracy of the blade vibration in high frequency can be improved. The method for calculating the vibration acceleration amplitude of the blade by adopting the circumferential Fourier method can ensure the accuracy of a calculation result to the greatest extent, and deduce the vibration parameter of the blade by utilizing the relation among the acceleration, the speed and the displacement.
Drawings
In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the drawings required for the embodiments will be briefly described below, and it is apparent that the drawings in the following description are only some embodiments described in the present application, and other drawings may be obtained according to these drawings for a person having ordinary skill in the art.
Various other advantages and benefits of the present application will become apparent to those of ordinary skill in the art upon reading the following detailed description of the preferred embodiments. The drawings are only for purposes of illustrating the preferred embodiments and are not to be construed as limiting the application. It is evident that the figures described below are only some embodiments of the application, from which other figures can be obtained without inventive effort for a person skilled in the art. Also, like reference numerals are used to designate like parts throughout the figures.
In the drawings:
FIG. 1 is a schematic flow diagram of a turbine blade tip timing on-line monitoring system according to one embodiment of the present application;
FIG. 2 is a diagram of a turbine blade monitoring system in accordance with one embodiment of the present application;
FIG. 3 is a graph of simulated acceleration results according to one embodiment of the present application;
FIG. 4 is a graph of simulated acceleration phase results according to one embodiment of the present application.
The application is further explained below with reference to the drawings and examples.
Detailed Description
For the purpose of making the objects, technical solutions and advantages of the embodiments of the present application more apparent, the technical solutions in the embodiments of the present application will be clearly and completely described in conjunction with the embodiments of the present application, and it is apparent that the described embodiments are some embodiments of the present application, but not all embodiments. All other embodiments, based on the embodiments of the application, which are apparent to those of ordinary skill in the art without inventive faculty, are intended to be within the scope of the application.
Thus, the following detailed description of the embodiments of the application, as presented in the figures, is not intended to limit the scope of the application, as claimed, but is merely representative of selected embodiments of the application. All other embodiments, based on the embodiments of the application, which are apparent to those of ordinary skill in the art without inventive faculty, are intended to be within the scope of the application.
It should be noted that: like reference numerals and letters denote like items in the following figures, and thus once an item is defined in one figure, it need not be further defined and explained in the following figures.
In the description of the present application, it should be understood that the terms "center", "longitudinal", "lateral", "length", "width", "thickness", "upper", "lower", "front", "rear", "left", "right", "vertical", "horizontal", "top", "bottom", "inner", "outer", "clockwise", "counterclockwise", etc. indicate orientations or positional relationships based on the orientations or positional relationships shown in the drawings are merely for convenience in describing the present application and simplifying the description, and do not indicate or imply that the apparatus or elements referred to must have a specific orientation, be configured and operated in a specific orientation, and thus should not be construed as limiting the present application.
Furthermore, the terms "first," "second," and the like, are used for descriptive purposes only and are not to be construed as indicating or implying a relative importance or implicitly indicating the number of technical features indicated. Thus, a feature defining "a first" or "a second" may explicitly or implicitly include one or more such feature. In the description of the present application, the meaning of "a plurality" is two or more, unless explicitly defined otherwise.
In the present application, unless explicitly specified and limited otherwise, the terms "mounted," "connected," "secured," and the like are to be construed broadly, and may be, for example, fixedly connected, detachably connected, or integrally formed; can be directly connected or indirectly connected through an intermediate medium, and can be communicated with the inside of two elements or the interaction relationship of the two elements. The specific meaning of the above terms in the present application can be understood by those of ordinary skill in the art according to the specific circumstances.
In the present application, unless expressly stated or limited otherwise, a first feature "above" or "below" a second feature may include both the first and second features being in direct contact, as well as the first and second features not being in direct contact but being in contact with each other through additional features therebetween. Moreover, a first feature being "above," "over" and "on" a second feature includes the first feature being directly above and obliquely above the second feature, or simply indicating that the first feature is higher in level than the second feature. The first feature being "under", "below" and "beneath" the second feature includes the first feature being directly under and obliquely below the second feature, or simply means that the first feature is less level than the second feature.
In order to make the technical scheme of the present application better understood by those skilled in the art, the present application will be further described in detail with reference to the accompanying drawings. As shown in fig. 1 to 4, a method for monitoring the tip timing of a non-rotational speed reference turbine blade includes,
in the first step S1, a finite element model of a turbine blade is established, modal analysis is carried out on the finite element model, and a vibration mode of the working speed of 3000r/min is calculated to obtain a blade vibration mode parameter f= [ f ] 1 f 2 …f m ]；
In the second step S2, k groups of blade tip timing sensors are circumferentially arranged on a rotor blade casing, a test is carried out, and a time sequence of the rotor blade reaching the sensors is obtained;
in the third step S3, based on the identification requirement of the vibration parameters of the turbine blade, respectively utilizing the arrival time and the installation angle of the sensor to convert the measured physical quantity into the blade tip vibration acceleration;
in the fourth step S4, parameter identification of tip timing is completed by using the obtained different tip vibration accelerations and using a circumferential fourier algorithm, so as to obtain information such as vibration frequency, amplitude, phase and the like of the turbine rotor blade at a constant rotation speed.
In a preferred embodiment of the method, in a first step S1, finite element modal analysis is performed on the rotor blade to obtain the number m of modes of the turbine rotor blade under a given working condition; in a second step S2, a tip timing sensor group number k is determined from the number m of steam turbine rotor blade modes, wherein k=3m; the number of sensors in each group is three, and the angular interval between the sensors in each group should be less than 5 degrees, and the smaller the condition of meeting the practical situation, the better.
In a preferred embodiment of the method, in a first step S1, an nth turn sensor angle matrix is constructed according to the frequencies of different modal vibrations and the arrival times of different sensors
wherein fm The vibration frequency of m-order modes, m is the mode number, t k,n Average value of blade arrival time acquired at nth turn for kth group of sensors, wherein For the time acquired in the nth turn for the smaller sensor angle of the kth group,/->For the time acquired in the nth turn, which is the middle of the sensor angle of the kth group,/is>Is the k groupThe time acquired at the nth turn is the time when the sensor angle is the largest.
In a preferred embodiment of the method, in a second step S2, a tip timing sensor is mounted on a turbine casing, a tip timing test is performed, and the time for the blade to reach the sensor is measured;
s501, the tip timing system is not provided with an OPR rotation speed sensor;
in a preferred embodiment of the method, in a third step S3, according to vibration parameter identification requirements of different modes, blade tip vibration acceleration is calculated by using arrival time and a sensor installation angle;
s601, calculating the vibration acceleration of the blade tip in the nth turn by using the angles among the kth group of sensors and the arrival time acquired by the kth group of sensors according to the arrival time of the blade at the sensors and the angles among the sensors wherein ,r is the radius of the steam turbine;
s6011, turbine speed f Ω Constant;
in a fourth step S4, according to the obtained blade tip vibration speed, blade tip timing vibration parameters are obtained through a circumferential Fourier solution; the vibration parameter is calculated by the following formula:
a n ＝Φ n S n, wherein wherein />Vibration acceleration values obtained for the kth group of sensors at the nth turn; s is S n ＝[i 1 j 1 i 2 j 2 …i m j m ] T ，S n The vibration parameter matrix is the nth circle;
s701, according to the conversion relation between the vibration displacement and the acceleration of the blade tip, finally obtaining the vibration displacement parameter of the blade;
s7011 vibration acceleration amplitude values under different blade vibration modesVibration velocity phase
S7012 vibration speed amplitude under different blade vibration modesVibration velocity phase
S7012 vibration displacement amplitude values under different blade vibration modesVibration velocity phase
The measurement system comprises a measurement device which is arranged to measure,
blade modeling and finite element analysis module: establishing a finite element model of a turbine blade, and performing modal analysis to obtain a blade vibration modal parameter f= [ f ] 1 f 2 …f m ]；
Sensor layout module: the grouped turbine blade tip timing sensors are arranged on the casing of the rotary blade;
constructing a sensor angle matrix from modal frequencies and arrival times
Arrival time acquisition module: the arrival time sequence of the single blade acquired by the group of sensors is
Vibration acceleration calculation module: calculating vibration velocity by using arrival time, wherein the vibration velocity sequence of the single blade is a= [ a ] 1 a 2 …a n ]；
Vibration parameter acquisition module: solving parameter S by using circumferential Fourier according to sensor angle matrix and vibration speed n ＝[i 1 j 1 i 2 j 2 …i m j m ] T ；
The application will be further described with reference to fig. 1 and 2, with a specific simulation embodiment instead of the experimental example, to clearly illustrate the steps of the application in use on a steam turbine.
FIG. 1 is a flow chart of a turbine blade tip timing on-line monitoring method of the application, FIG. 2 is a diagram of a turbine blade monitoring system in one embodiment of the application, finite element mode analysis is firstly carried out on blades to obtain the mode of blade vibration, in the simulation, the blade vibration mode and the blade vibration function are set, the number of required sensors and the installation angle are further determined, after the sensors 2 are installed around a casing 1, a turbine is operated, a blade 4 is driven by a turbine main shaft 3 to rotate, the sensors obtain the blade arrival time, then a vibration acceleration conversion formula is utilized to obtain the blade vibration acceleration corresponding to each pair of sensors, and then the blade vibration parameters are solved according to a sensor angle matrix determined by the installation number of the blade tip timing sensors. The method comprises the following specific steps:
1) The number of vibration modes was set to 1, and the vibration function was set to 0.001sin (2000 pi t+30).
2) The number of sensors was determined to be 6 and their angles were 40, 42, 44, 83, 85, 87, respectively.
3) The rotating speed of the turbine is 3000r/min, the blades rotate for 1s, the blade vibration is converted into an angle domain signal, the arrival time obtained by a 50-circle sensor of blade rotation is obtained, and the arrival time is obtained
4) The acceleration of the blade tip vibration can be obtained by using the arrival time
5) Averaging t to obtain wherein />
6) Constructing a sensor angle matrix
7) Solving for
8) Vibration acceleration amplitudeVibration acceleration phase->
9) Vibration velocity amplitude
10 Amplitude of vibration displacementThe results are shown in fig. 4, and the simulation results are substantially the same as the simulation parameter settings.
Finally, it should be noted that: the described embodiments are intended to be illustrative of only some, but not all, of the embodiments of the present application and, based on the embodiments herein, all other embodiments that may be made by those skilled in the art without the benefit of the present disclosure are intended to be within the scope of the present application.
While certain exemplary embodiments of the present application have been described above by way of illustration only, it will be apparent to those of ordinary skill in the art that modifications may be made to the described embodiments in various different ways without departing from the spirit and scope of the application. Accordingly, the drawings and description are to be regarded as illustrative in nature and not as restrictive of the scope of the application, which is defined by the appended claims.

Claims (9)

1. The method for monitoring the blade tip timing of the turbine blade without the rotating speed reference is characterized by comprising the following steps of:
step S1, a finite element model of a turbine blade is established, modal analysis is carried out on the finite element model, and the vibration mode of the turbine at a preset working speed is calculated to obtain the blade vibration mode parameter f= [ f ] 1 f 2 … f m ]；
Step S2, circumferentially installing k groups of blade tip timing sensors on a rotor blade casing, and acquiring a time sequence of the rotor blade reaching the blade tip timing sensors;
step S3, calculating to obtain the vibration acceleration of the blade by using the time sequence and the installation angle of the blade tip timing sensor;
and a fourth step S4, obtaining the vibration frequency, amplitude and phase of the turbine blade at a constant rotating speed through a circumferential Fourier algorithm based on the blade tip vibration acceleration.
2. The method according to claim 1, characterized in that preferably in the first step S1, the number m of blade modes of the steam turbine is obtained by a mode analysis, and the number k of blade tip timing sensor groups is determined according to the number m of blade modes.
3. The method of claim 2 wherein the number of tip timing sensor groups k and the number of blade modes m are related by k = 3m.
4. The method according to claim 1, wherein in the first step S1, the number of tip timing sensors per group is three and the angular interval of the tip timing sensors per group is less than 5 degrees.
5. The method of claim 1, wherein the predetermined operating speed of the turbine is 3000r/min.
6. The method of claim 1, wherein an nth turn tip timing sensor angle matrix is constructed from blade vibration modal parameters and time series
wherein fm The vibration frequency of m-order modes, m is the mode number, t k,n Average of blade arrival times acquired for the kth group of tip timing sensors at the nth turn, where For the time acquired in the nth turn for the smaller angle of the kth group tip timing sensor, +.>Time acquired in the nth turn for the middle of the kth group tip timing sensor angle, +.>The time acquired at the nth turn is the one with the greatest angle of the kth group of tip timing sensors.
7. The method according to claim 6, wherein in the third step S3, the vibration acceleration of the blade tip in the nth turn is calculated using the angle between the kth group of tip timing sensors and the acquired arrival time wherein ,r is the radius of the steam turbine, theta k,1 、θ k,2 、θ k,3 Respectively the angles of the kth group of sensors after the sensor is sequenced according to the angles, and the rotating speed f of the steam turbine Ω Constant.
8. The method according to claim 7, wherein in the fourth step S4, the tip-timing vibration parameter is obtained via a circumferential fourier algorithm based on the tip-vibration acceleration, the tip-timing vibration parameter including:
a n ＝Φ n S n, wherein wherein />Vibration acceleration values obtained for the kth group of blade tip timing sensors at the nth turn; s is S n ＝[i 1 j 1 i 2 j 2 … i m j m ] T ，S n The vibration parameter matrix is the nth circle;
vibration acceleration amplitude values under different blade vibration modesVibration velocity phase
Vibration velocity amplitude under different blade vibration modesVibration velocity phase->
Vibration displacement amplitude values under different blade vibration modesVibration velocity phase
9. A system for implementing the method for monitoring the blade tip timing of a rotational speed-free reference turbine blade according to any one of claims 1 to 8, comprising,
the blade modeling and finite element analysis module is configured to establish a finite element model of a turbine blade and perform modal analysis to obtain a blade vibration modal parameter f= [ f ] 1 f 2 … f m ]；
A sensor layout module connected to the blade modeling and finite element analysis module, the sensor layout module determining a number of tip timing sensors based on the vibration mode parameters and constructing a sensor angle matrix based on the blade vibration mode parameters and the arrival time
A time of arrival acquisition module configured to acquire a time of arrival sequence of the blade through the tip timing sensor, the individual blade time of arrival sequences acquired by the set of sensors being
The vibration acceleration calculation module is connected with the arrival time acquisition module and calculates the vibration speed by using the arrival time, and the vibration speed sequence of the single blade is a= [ a ] 1 a 2 … a n ]；
The vibration parameter acquisition module is connected with the sensor layout module and the vibration acceleration calculation module, and the vibration parameter acquisition module solves the vibration parameter S by using a circumferential Fourier according to the sensor angle matrix and the vibration speed n ＝[i 1 j 1 i 2 j 2 … i m j m ] T 。
CN202310835523.7A 2023-07-10 2023-07-10 Method and system for timing monitoring of blade tips of rotational speed-free reference turbine blade Pending CN116877212A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310835523.7A CN116877212A (en) 2023-07-10 2023-07-10 Method and system for timing monitoring of blade tips of rotational speed-free reference turbine blade

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310835523.7A CN116877212A (en) 2023-07-10 2023-07-10 Method and system for timing monitoring of blade tips of rotational speed-free reference turbine blade

Publications (1)

Publication Number Publication Date
CN116877212A true CN116877212A (en) 2023-10-13

Family

ID=88263733

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310835523.7A Pending CN116877212A (en) 2023-07-10 2023-07-10 Method and system for timing monitoring of blade tips of rotational speed-free reference turbine blade

Country Status (1)

Country Link
CN (1) CN116877212A (en)

Similar Documents

Publication Publication Date Title
CN101625260B (en) Method for detecting high speed rotating blade synchronous vibration parameters under speed change
CN110608710B (en) Rotor blade dynamic strain field measuring method and system based on blade end timing
CN101629846B (en) Method for detecting synchronous vibration parameters of high-speed rotary blade under constant speed
CN109883380B (en) Rotor blade displacement field measuring method and system based on blade end timing
CN101929917B (en) Fault diagnosis method for rotary machine
CN110851963A (en) Casing circumferential arrangement method of blade end timing sensor
CN102721462B (en) Method for quickly computing Bode plot and Nyquist plot of rotary mechanical vehicle starting and parking processes
CN105234745B (en) Machine tool spindle thermal error, deviation from circular from are separated and processing method with turn error
CN111397877B (en) Rotary machine beat vibration fault detection and diagnosis method
CN111507043A (en) Rotor blade dynamic stress field measuring method and system based on blade end timing
CN109540482B (en) Method and device for analyzing keyless synchronous vibration parameters of turbine blade
CN109115408B (en) Dynamic balance test method of large hydroelectric generating set based on centrifugal force equation
CN110553844A (en) Method and system for detecting misalignment fault of rotary machine
CN108871543A (en) The anharmonic Fourier's analysis method of blade asynchronous vibration frequency under constant speed
Fan et al. Shifting straight-line fitting method to calculate blade vibration based on blade tip timing
CN116577050A (en) Method for testing dynamic damping ratio of aero-engine compressor rotor blade
CN113190786B (en) Vibration prediction method for large-scale rotating equipment by utilizing multidimensional assembly parameters
CN106092534A (en) Blade modal damping assay method
CN116877212A (en) Method and system for timing monitoring of blade tips of rotational speed-free reference turbine blade
CN116953279A (en) Timing monitoring method and system for high-frequency vibration blade tip of turbine blade of aeroengine
CN115081271B (en) Leaf end timing system checking method and checking system based on digital simulator
CN117030004A (en) Turbine blade tip timing on-line monitoring method and system
CN117147057B (en) Opposite rotating propeller fan axis track reconstruction method based on Mems triaxial vibration acceleration sensor
CN112903271A (en) Non-contact asynchronous vibration parameter identification method for rotor blade
Liu et al. Tip timing based non-contact vibration measurement of aero-engine turbine blades

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