Impeller non-constant speed rotation blade deformation estimation method, system and medium

Abstract

A blade deformation estimation method, system and medium under non-constant speed rotation of an impeller, structural constraint parameters and evaluation strategy parameters of an impeller rotor are obtained; the pulse count quantity and compensation time difference when each rotation period key phase zero point passes are extracted from the key phase sensor response to form a first sequence; the pulse count quantity and compensation time difference when each blade passes the sensor are extracted from the blade sensor response to form a second sequence; based on the first sequence and the second sequence, the key phase zero point passing time of other rotation periods and the measured time when each blade passes each blade sensor are calculated with the key phase time of the starting calculation rotation period as the reference zero time; the instantaneous speed sequence of the impeller is generated based on the key phase zero point passing time sequence; the impeller rotation phase when each blade reaches each blade sensor is calculated based on the structural constraint parameters; the theoretical time and instantaneous speed are solved based on the key phase zero point passing time sequence, the instantaneous speed sequence and the impeller rotation phase sequence; the blade deformation is estimated.

Description

Technical Field

[0001] This invention relates to the field of blade vibration monitoring technology, and in particular to a method, system and medium for estimating blade deformation under non-constant speed rotation of the impeller. Background Technology

[0002] As a core component of high-speed rotating machinery such as aero-engines and compressors, the vibration state of blades directly determines the reliability and safety of equipment operation. Currently, non-contact measurement based on the tip timing principle is the mainstream technology for blade vibration monitoring. This method infers deformation by measuring the difference between the actual time it takes for the blade to pass through the blade sensor and the theoretically calculated time. While this method performs well under constant speed conditions, it faces significant challenges under common non-constant speed conditions such as wind turbine start-up and shutdown, and variable loads. Traditional methods assume a constant speed or use average speed for simplified calculations, leading to severe distortion in the estimation of the theoretical blade passage time, which in turn significantly reduces the accuracy of subsequent deformation extraction.

[0003] Under non-constant rotational speeds, the measurement error of the tip-time method mainly stems from three interdependent key factors: First, the dynamic changes in rotational speed render the theoretical passage time calculated based on the constant rotational speed model meaningless, with systematic errors directly infiltrating the deformation signal; second, the discrete sampling characteristics of the blade sensor lead to asynchrony between the sampling pulse and the actual passage time of the blade, and the resulting "compensation time difference" couples with rotational speed fluctuations under variable speed conditions, making it difficult for traditional methods to accurately compensate; third, deformation calculation relies on instantaneous linear velocity, and insufficient accuracy in estimating instantaneous rotational speed directly amplifies the error in the final result. These problems collectively cause the vibration response equation to exhibit ill-conditioned characteristics, making the estimation results of key parameters such as vibration amplitude and frequency highly sensitive to measurement noise.

[0004] Under non-constant speed conditions, existing technologies suffer from a significant decrease in blade deformation estimation accuracy due to the coupling and inability to accurately decouple the theoretical time estimation distortion, asynchronous sampling bias, and instantaneous speed estimation error. Summary of the Invention

[0005] To address the aforementioned technical problems, this invention provides a method for estimating blade deformation under non-constant speed rotation of an impeller, comprising the following steps:

[0006] S100. Obtain the structural constraint parameters and evaluation strategy parameters of the impeller rotor. The structural constraint parameters include the phase sequence of the blade relative to the key phase zero point, the position sequence of the blade sensor relative to the key phase sensor, the radius of the cross-sectional circle swept by the blade tip, and the number of blades. The evaluation strategy parameters include the number of impeller rotation cycles, the sampling frequency of the key phase sensor, and the sampling frequency of the blade sensor.

[0007] S200: Continuously acquires response data from the key phase sensor and each blade sensor;

[0008] S300. Extract from the key phase sensor response the number of key phase sampling pulse counts relative to the key phase zero trigger time of the previous rotation cycle when the key phase zero point passes through, and the compensation time difference from it to the next key phase sampling pulse, to form a first sequence.

[0009] S400. Extract from the blade sensor response the number of blade sampling pulse counts relative to the current rotation cycle key phase zero-point trigger time of each blade when it passes the blade sensor, and the compensation time difference from the blade to the next blade sampling pulse, to form a second sequence.

[0010] S500. Based on the first sequence and the second sequence, select the bond phase time at the start of the calculation of the rotation period as the reference zero time, and calculate the bond phase zero point passing time of other rotation periods related to the period to be solved, as well as the measured time of each blade passing through different blade sensors.

[0011] S600: Based on the key phase zero-point passage time sequence, generate the instantaneous speed sequence of the impeller;

[0012] S700: Based on structural constraint parameters, calculate the impeller rotation phase when each blade reaches its respective blade sensor;

[0013] S800, based on the key phase zero-point passage time sequence, instantaneous speed sequence, and impeller rotation phase sequence, uses interpolation to solve for the theoretical time and corresponding instantaneous speed of the blade passing through the corresponding blade sensor for different rotation cycles;

[0014] S900. Based on the difference between the measured time and the theoretical time, the radius of the cross-sectional circle swept by the blade tip and the instantaneous rotation speed, the blade deformation is estimated.

[0015] Optionally, the sampling frequency of the blade sensor is an integer multiple of the sampling frequency of the key phase sensor.

[0016] Optionally, the S300 includes the following steps:

[0017] S310. For each rotation cycle, identify the response pulse generated when the key phase sensor captures the key phase zero point;

[0018] S320. Record the number of sampling pulses relative to the key phase zero-point trigger time of the response pulse, based on the sampling frequency of the key phase sensor.

[0019] S330, Record the compensation time difference from the moment the key phase sensor passes through the key phase zero point to the moment the next key phase sensor sampling pulse arrives;

[0020] S340, the pulse count based on S320, and the compensation time difference based on S330 form the first sequence.

[0021] Optionally, the S400 includes the following steps:

[0022] S410. For each blade sensor, sequentially identify the response pulse generated when each blade passes through the blade sensor;

[0023] S420. Record the number of sampled pulses relative to the current rotation cycle key phase zero-point trigger time, based on the blade sensor sampling frequency.

[0024] S430. Record the compensation time difference from the moment the blade passes the blade sensor to the moment the next blade sensor sampling pulse arrives;

[0025] S440, based on the pulse count of S420 and the compensation time difference of S430, forms a second sequence.

[0026] Optionally, in S500, the initial calculation rotation period RS is selected, with its bond phase time as the reference zero time. For any other rotation period after the initial calculation rotation period RS... Its bond phase zero point passes through time The corresponding calculation formula is as follows:

[0027] ;in, For rotation period The number of pulse counts corresponding to the sampling frequency SK of the key phase sensor when it passes through the key phase zero point twice in the next rotation cycle; For rotational period index; This is the compensation time difference between the current period's bond phase sensor passing through the bond phase zero point and the sampling pulse response of the next bond phase sensor; The compensation time difference for the periodic RS bond phase sensor from the bond phase zero point to the sampling pulse response of the next bond phase sensor;

[0028] Alternatively, for any other rotation period prior to the initial calculation of rotation period RS. Its bond phase zero point passes through time The corresponding calculation formula is as follows:

[0029] ;in, For rotation period The number of pulse counts corresponding to the sampling frequency SK of the key phase sensor when it passes through the key phase zero point twice in the next rotation cycle; For rotational period index; For period The key phase sensor compensates for the time difference between the sampling pulse response from the key phase zero point to the next key phase sensor. The compensation time difference is the pulse response from the zero point of the periodic RS bond phase sensor to the sampling pulse response of the next bond phase sensor.

[0030] Optionally, taking the key phase moment of the initial calculation of the rotation period RS as the reference zero moment, the response of the i-th blade sensor in any other rotation period after the initial calculation of the impeller's rotation period RS. The measured time when the p-th blade passes through the i-th blade sensor. The corresponding calculation formula is as follows:

[0031] Where BN represents the number of blades; When the p-th time a blade is detected passing the i-th blade sensor, the blade sensor relative to the i-th blade sensor... The number of pulses collected at the blade sampling frequency SB at the moment of passing through the zero point of the key phase; This is the blade correction time value at the moment when the corresponding blade passes the i-th blade sensor;

[0032] Alternatively, the i-th blade sensor response occurs during any other rotation period before the impeller begins calculating the rotation period RS. The measured time when the p-th blade passes through the i-th blade sensor. The corresponding calculation formula is as follows:

[0033] Where BN represents the number of blades; When the p-th time a blade is detected passing the i-th blade sensor, the blade sensor relative to the i-th blade sensor... The number of pulses collected at the blade sampling frequency SB at the moment of passing through the zero point of the key phase; This is the blade correction time value at the moment when the corresponding blade passes the i-th blade sensor.

[0034] Optionally, in the instantaneous rotational speed sequence, the first and last elements of the sequence are both set to 0;

[0035] For the initial calculated rotation period RS, the instantaneous rotational speed corresponding to the moment the key phase zero point passes through is... The calculation formula is as follows:

[0036] ;in, The time when the bond phase zero passes through the next rotational period corresponding to the initial calculation of the rotational period RS is calculated. The time when the bond phase zero point passes through the previous rotational period corresponding to the initial calculation of the rotational period RS;

[0037] For rotation period Its bond phase zero point passes through time Corresponding instantaneous speed The calculation formula is as follows:

[0038] ;in, For rotation period The bond phase zero point corresponding to the next rotation cycle passes through the time. For rotation period The time when the bond phase zero point corresponding to the previous rotation cycle passes through;

[0039] For rotation period Its bond phase zero point passes through time Corresponding instantaneous speed The calculation formula is as follows:

[0040] ;in, For rotation period The bond phase zero point corresponding to the next rotation cycle passes through the time. For rotation period The time when the bond phase zero point corresponding to the previous rotation cycle passes through.

[0041] Optionally, the S800 includes the following steps:

[0042] S810. Using the time sequence of the bond phase zero point passing through and its corresponding impeller rotation phase sequence as interpolation conditions, and the first derivative of the boundary point of the time sequence of the bond phase zero point passing through as boundary conditions, construct the first interpolation equation and solve for the theoretical time when the blade passes through the corresponding blade sensor.

[0043] S820. Using the time sequence of the key phase zero point and its corresponding instantaneous rotational speed sequence as interpolation conditions, and the first derivative of the boundary point of the instantaneous rotational speed sequence as boundary conditions, a second interpolation equation is constructed to solve for the instantaneous rotational speed of the blade passing through the corresponding blade sensor.

[0044] Corresponding to the blade deformation estimation method under non-constant speed rotation of the impeller, the present invention provides a blade deformation estimation system under non-constant speed rotation of the impeller, which includes:

[0045] The parameter acquisition module is used to acquire the structural constraint parameters and evaluation strategy parameters of the impeller rotor. The structural constraint parameters include the phase sequence of the blade relative to the key phase zero point, the position sequence of the blade sensor relative to the key phase sensor, the radius of the cross-sectional circle swept by the blade tip, and the number of blades. The evaluation strategy parameters include the number of impeller rotation cycles, the sampling frequency of the key phase sensor, and the sampling frequency of the blade sensor.

[0046] The data acquisition module is used to continuously acquire response data from the key phase sensor and each blade sensor;

[0047] The key phase signal processing module extracts from the key phase sensor response the number of key phase sampling pulses when the key phase zero point passes through each rotation cycle relative to the trigger time of the key phase zero point in the previous rotation cycle, as well as the compensation time difference from it to the next key phase sampling pulse, to form a first sequence;

[0048] The blade signal processing module extracts from the blade sensor response the number of blade sampling pulse counts relative to the current rotation cycle key phase zero-point triggering time of each blade when it passes the blade sensor, as well as the compensation time difference from the next blade sampling pulse, to form a second sequence;

[0049] The time calculation module is used to select the bond phase time of the starting calculation of the rotation period as the reference zero time based on the first sequence and the second sequence, and to calculate the bond phase zero point passage time of other rotation periods associated with the period to be solved, as well as the measured time of each blade passing through different blade sensors.

[0050] The instantaneous speed sequence generation module is used to generate the instantaneous speed sequence of the impeller based on the key phase zero-point passage time sequence;

[0051] The impeller rotation phase calculation module is used to calculate the impeller rotation phase when each blade reaches its respective blade sensor based on structural constraint parameters.

[0052] The interpolation solution module is used to solve the theoretical time and corresponding instantaneous speed of the blade passing through the corresponding blade sensor for different rotation cycles based on the key phase zero-point passage time sequence, instantaneous speed sequence, and impeller rotation phase sequence, using the interpolation method.

[0053] The deformation estimation module is used to estimate the blade deformation based on the difference between the measured and theoretical times, the radius of the cross-sectional circle swept by the blade tip, and the instantaneous rotational speed.

[0054] In addition, to achieve the above objectives, the present invention also provides a computer-readable storage medium storing a blade deformation estimation program for impeller non-constant speed rotation, wherein when the blade deformation estimation program for impeller non-constant speed rotation is executed by a processor, the program implements the steps of the blade deformation estimation method for impeller non-constant speed rotation as described above.

[0055] Compared with the prior art, the present invention has the following beneficial effects:

[0056] (1) By using a dual-parameter framework of structural constraint parameters and evaluation strategy parameters, geometric quantities such as the bond phase zero point, blade sensor, radius of the cross section circle swept by the blade tip, and number of blades are simultaneously incorporated into the model along with measurement strategies such as sampling frequency and number of rotation cycles. In each step, the time-domain-sampling domain coupled variables of pulse count and compensation time difference are used to reconstruct the time of the bond phase zero point and the time of each blade passing through each blade sensor in real time for each rotation cycle. Then, the instantaneous speed sequence is generated. Finally, the deformation is estimated by directly linking the difference between the theoretical time and the measured time with the instantaneous linear velocity. Thus, the three major coupled distortion sources of speed fluctuation, sampling asynchrony, and instantaneous speed error are eliminated synchronously in one calculation link, getting rid of the traditional constant speed or average speed assumption and significantly improving the blade deformation estimation accuracy.

[0057] (2) The sampling frequency of the blade sensor is set to an integer multiple of the sampling frequency of the key phase sensor, so that the sampling pulses of the two channels are naturally aligned at integer multiple times, eliminating the cumulative phase slip caused by non-integer multiples of frequency, fundamentally suppressing the secondary coupling error of compensation time difference under variable speed conditions, and further improving the robustness of deformation estimation.

[0058] (3) By establishing a first sequence of pulse count and compensation time difference separately in the key phase channel, the sub-sampling level error between the actual passing time of the key phase zero point and the discrete sampling pulse can be quantified and retained at one time, providing a unified and traceable benchmark for all subsequent time axis calculations.

[0059] (4) A second sequence of pulse counting and time difference compensation is also established in the blade channel to ensure that the actual blade passage time and the reference time of the bond phase maintain the same subsampling accuracy, so as to ensure that the two time chains are aligned at the same error level and eliminate the asynchronous error between the blade and the bond phase across the channel.

[0060] (5) Introduce the blade correction time value into the blade passage time formula, dynamically couple the number of blades, the geometric position of the blade sensor and the key phase reference, and ensure that the output time is aligned with the key phase reference when any blade passes any blade sensor any time, thus completely eliminating the fixed phase error caused by the difference in geometric position.

[0061] (6) The instantaneous speed sequence is set to zero at the beginning and end, and other sequence elements are calculated using specific formulas. This can reduce the sensitivity of speed estimation accuracy to noise without adding extra sensors, while ensuring the smoothness of the speed curve endpoints and avoiding the amplification effect of boundary abrupt changes on deformation estimation.

[0062] (7) By using the three sets of sequences of bond phase, rotational speed and phase as interpolation conditions and adopting the first derivative boundary constraint, the theoretical time and instantaneous rotational speed are solved simultaneously in the same cubic spline framework, ensuring that the two are strictly conjugate mathematically and eliminating the cumulative interpolation error caused by the existing technology of calculating the rotational speed first and then the time. Attached Figure Description

[0063] The accompanying drawings, which are included to provide a further understanding of the invention and form part of this invention, illustrate exemplary embodiments of the invention and are used to explain the invention, but do not constitute an undue limitation of the invention. In the drawings:

[0064] Figure 1 This is a simplified flowchart of an embodiment of the blade deformation estimation method under non-constant speed rotation of the impeller according to the present invention;

[0065] Figure 2 This is a framework diagram of an embodiment of the blade deformation estimation system under non-constant speed rotation of the impeller according to the present invention. Detailed Implementation

[0066] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, not all embodiments. It should be understood that the specific embodiments described herein are only used to explain the present invention and are not intended to limit the present invention. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0067] The impeller of the equipment is usually equipped with multiple blades. In order to correlate a specific blade with its vibration state, each blade needs to be numbered before using the traditional blade tip timing method to measure vibration. The specific number of the blade needs to be determined by comprehensively considering structural parameters such as the positional relationship of the blade relative to the key phase point on the impeller rotor (the edge of the groove or protrusion monitored by the key phase sensor), and the direction of impeller rotation.

[0068] In blade vibration measurement, the key phase sensor is used to capture the response time of the impeller rotor key phase point passing through the key phase sensor, serving as a "reference benchmark" for the blade sensor response time. Simultaneously, the key phase time series can be used to estimate the equipment's rotational speed. The blade sensor, on the other hand, captures the response time of different blade tips passing through the blade sensor. Since the blade's position relative to the key phase point is fixed, by adding the positional relationship between the blade sensor and the key phase sensor, as well as the impeller rotational speed, the theoretical time (under undeformed blade conditions) of the blade passing through different blade sensors can be calculated, using the key phase sensor's capture time of the key phase point as a reference. When the blade deforms, the time measured by the blade sensor will "deviation" from the theoretical time, indicating blade deformation. It can be expressed as follows:

[0069] ;in, The radius of the cross-sectional circle swept by the blade tip. The impeller speed, The theoretical time deviation of the j-th blade passing through the i-th blade sensor is based on the time corresponding to the key phase point captured by the key phase sensor when the blade deforms.

[0070] Blade vibration assessment based on tip timing method, i.e., from the measured deformation vector [ The vibration state of the blade is estimated in [the data]. Specifically, during the estimation, the blade deformation of the j-th blade is measured by the sensor on the i-th blade. It is usually divided into two parts: the blade static deformation component caused by the average load and the blade vibration deformation component caused by the dynamic load. It can be represented as:

[0071] ;

[0072] Where: FN is the total number of vibration orders in the blade vibration assessment, and k represents the sequence number of the vibration order; For the j-th blade at vibration frequency The estimated value of the vibration amplitude; Let be the angle between the i-th blade sensor and the key phase sensor along the rotor rotation direction; The phase of the j-th blade relative to the zero point of the key phase in the opposite direction of the impeller rotation; The vibration frequency is At that time, the estimated value of the initial phase of vibration corresponding to the response time of the key phase sensor; This is an estimate of the static deformation component of the blade.

[0073] In the above-mentioned blade deformation and vibration estimation, the accurate capture of blade passage time and the estimation of theoretical blade passage time under varying rotational speed conditions are the main factors affecting blade deformation calculation, which in turn affect the extraction of blade vibration parameters. Based on this, the present invention provides a method for estimating blade deformation under non-constant speed rotation of the impeller, comprising the following steps:

[0074] S100. Obtain the structural constraint parameters and evaluation strategy parameters of the impeller rotor. The structural constraint parameters include the phase sequence of the blade relative to the key phase zero point, the position sequence of the blade sensor relative to the key phase sensor, the radius of the cross-sectional circle swept by the blade tip, and the number of blades. The evaluation strategy parameters include the number of impeller rotation cycles, the sampling frequency of the key phase sensor, and the sampling frequency of the blade sensor.

[0075] S200: Continuously acquires response data from the key phase sensor and each blade sensor;

[0076] S300. Extract from the key phase sensor response the number of key phase sampling pulse counts relative to the key phase zero trigger time of the previous rotation cycle when the key phase zero point passes through, and the compensation time difference from it to the next key phase sampling pulse, to form a first sequence.

[0077] S400. Extract from the blade sensor response the number of blade sampling pulse counts relative to the current rotation cycle key phase zero-point trigger time of each blade when it passes the blade sensor, and the compensation time difference from the blade to the next blade sampling pulse, to form a second sequence.

[0078] S500. Based on the first sequence and the second sequence, select the bond phase time at the start of the calculation of the rotation period as the reference zero time, and calculate the bond phase zero point passing time of other rotation periods related to the period to be solved, as well as the measured time of each blade passing through different blade sensors.

[0079] S600: Based on the key phase zero-point passage time sequence, generate the instantaneous speed sequence of the impeller;

[0080] S700: Based on structural constraint parameters, calculate the impeller rotation phase when each blade reaches its respective blade sensor;

[0081] S800, based on the key phase zero-point passage time sequence, instantaneous speed sequence, and impeller rotation phase sequence, uses interpolation to solve for the theoretical time and corresponding instantaneous speed of the blade passing through the corresponding blade sensor for different rotation cycles;

[0082] S900. Based on the difference between the measured time and the theoretical time, the radius of the cross-sectional circle swept by the blade tip and the instantaneous rotation speed, the blade deformation is estimated.

[0083] This invention employs a dual-parameter framework of structural constraint parameters and evaluation strategy parameters. It integrates geometric quantities such as the bond phase zero point, blade sensor, radius of the cross-sectional circle swept by the blade tip, and number of blades with measurement strategies such as sampling frequency and number of rotation cycles into the model. At each step, it utilizes a pair of time-domain-sampling-domain coupled variables—pulse count and compensation time difference—to reconstruct in real time the passing time of the bond phase zero point and the passing time of each blade through its respective blade sensor for each rotation cycle. This generates an instantaneous rotational speed sequence. Finally, the difference between the theoretical and measured times is directly linked to the instantaneous linear velocity to estimate the deformation. This eliminates the three major coupled distortion sources—rotational speed fluctuation, sampling asynchrony, and instantaneous speed error—in a single computational link, freeing it from the traditional assumption of constant or average rotational speed and significantly improving the accuracy of blade deformation estimation.

[0084] In this embodiment, the sampling frequency of the blade sensor is an integer multiple of the sampling frequency of the key phase sensor.

[0085] The sampling frequency of the blade sensor in this invention is set to an integer multiple of the sampling frequency of the key phase sensor, so that the sampling pulses of the two channels are naturally aligned at integer multiple times, eliminating the cumulative phase slip caused by non-integer multiples of the frequency, fundamentally suppressing the secondary coupling error of the compensation time difference under variable speed conditions, and further improving the robustness of deformation estimation.

[0086] Preferably, the S300 can use a data acquisition device with a counting function to measure the time it takes for the key phase sensor to pass through the key phase zero point. The counter starts triggering when the key phase sensor passes through the key phase zero point, and the counter stops counting when the SK sampling pulse arrives.

[0087] In this embodiment, S310, for each rotation cycle, the response pulse generated when the key phase sensor captures the key phase zero point is identified;

[0088] S320. Record the number of sampling pulses relative to the key phase zero-point trigger time of the response pulse, based on the sampling frequency of the key phase sensor.

[0089] S330, Record the compensation time difference from the moment the key phase sensor passes through the key phase zero point to the moment the next key phase sensor sampling pulse arrives;

[0090] S340, the pulse count based on S320, and the compensation time difference based on S330 form the first sequence.

[0091] This invention establishes a first sequence of pulse counts and compensation time differences separately within the bond phase channel. This allows for the one-time quantization and retention of the subsampling level error between the actual time of the bond phase zero point and the discrete sampling pulses, providing a unified and traceable benchmark for all subsequent time axis calculations.

[0092] In this embodiment, S400 includes the following steps:

[0093] S410. For each blade sensor, sequentially identify the response pulse generated when each blade passes through the blade sensor;

[0094] S420. Record the number of sampled pulses relative to the current rotation cycle key phase zero-point trigger time, based on the blade sensor sampling frequency.

[0095] S430. Record the compensation time difference from the moment the blade passes the blade sensor to the moment the next blade sensor sampling pulse arrives;

[0096] S440, based on the pulse count of S420 and the compensation time difference of S430, forms a second sequence.

[0097] This invention also establishes a second sequence of pulse counting and time difference compensation within the blade channel, so that the actual blade passage time and the bond phase reference time maintain the same subsampling accuracy, ensuring that the two time chains are aligned at the same error level and eliminating the blade-bond phase cross-channel asynchronous error.

[0098] In this embodiment, in S500, the initial calculation rotation period RS is selected, and its bond phase time is used as the reference zero time. For any other rotation period after the initial calculation rotation period RS... Its bond phase zero point passes through time The corresponding calculation formula is as follows:

[0099] ;in, For rotation period The number of pulse counts corresponding to the sampling frequency SK of the key phase sensor when it passes through the key phase zero point twice in the next rotation cycle; For rotational period index; This is the compensation time difference between the current period's bond phase sensor passing through the bond phase zero point and the sampling pulse response of the next bond phase sensor; The compensation time difference for the periodic RS bond phase sensor from the bond phase zero point to the sampling pulse response of the next bond phase sensor;

[0100] Alternatively, for any other rotation period prior to the initial calculation of rotation period RS. Its bond phase zero point passes through time The corresponding calculation formula is as follows:

[0101] ;in, For rotation period The number of pulse counts corresponding to the sampling frequency SK of the key phase sensor when it passes through the key phase zero point twice in the next rotation cycle; For rotational period index; For period The key phase sensor compensates for the time difference between the sampling pulse response from the key phase zero point to the next key phase sensor. The compensation time difference is the pulse response from the zero point of the periodic RS bond phase sensor to the sampling pulse response of the next bond phase sensor.

[0102] In this embodiment, the key phase moment at which the initial calculation of the rotation period RS is taken as the reference zero moment, and the response of the i-th blade sensor is any other rotation period after the impeller's initial calculation of the rotation period RS. The measured time when the p-th blade passes through the i-th blade sensor. The corresponding calculation formula is as follows:

[0103] Where BN represents the number of blades; When the p-th time a blade is detected passing the i-th blade sensor, the blade sensor relative to the i-th blade sensor... The number of pulses collected at the blade sampling frequency SB at the moment of passing through the zero point of the key phase; This is the blade correction time value at the moment when the corresponding blade passes the i-th blade sensor;

[0104] Alternatively, the i-th blade sensor response occurs during any other rotation period before the impeller begins calculating the rotation period RS. The measured time when the p-th blade passes through the i-th blade sensor. The corresponding calculation formula is as follows:

[0105] Where BN represents the number of blades; When the p-th time a blade is detected passing the i-th blade sensor, the blade sensor relative to the i-th blade sensor... The number of pulses collected at the blade sampling frequency SB at the moment of passing through the zero point of the key phase; This is the blade correction time value at the moment when the corresponding blade passes the i-th blade sensor.

[0106] In this embodiment, the measured time is the time when a certain blade passes through a certain blade sensor minus the number of pulses detected by that blade sensor; specifically, [ ] vector or [ In the vector, for the i-th blade sensor, when the j-th blade first reaches the i-th blade sensor, the number of pulses it has detected is... ,and Where floor represents the integer corresponding to the remainder being greater than zero. In [ ] vector or [ In the vector, the corresponding front end is removed. The pulse is used to obtain the response time vector of each blade in the corresponding rotation cycle, starting from the j-th blade. or response time vector That is, the actual measurement time.

[0107] This invention introduces a blade correction time value into the blade passage time formula, dynamically coupling the number of blades, sensor geometry, and key phase reference. This ensures that the output time is aligned with the key phase reference whenever any blade passes any blade sensor, thus completely eliminating the fixed phase error caused by differences in geometric position.

[0108] In this embodiment, the instantaneous rotational speed sequence is taken , , ..., In the instantaneous rotational speed sequence, the first and last elements of the sequence (i.e., and All are set to 0;

[0109] For the initial calculated rotation period RS, the instantaneous rotational speed corresponding to the moment the key phase zero point passes through is... The calculation formula is as follows:

[0110] ;in, The time when the bond phase zero passes through the next rotational period corresponding to the initial calculation of the rotational period RS is calculated. The time when the bond phase zero point passes through the previous rotational period corresponding to the initial calculation of the rotational period RS;

[0111] For rotation period Its bond phase zero point passes through time Corresponding instantaneous speed The calculation formula is as follows:

[0112] ;in, For rotation period The bond phase zero point corresponding to the next rotation cycle passes through the time. For rotation period The time when the bond phase zero point corresponding to the previous rotation cycle passes through;

[0113] For rotation period Its bond phase zero point passes through time Corresponding instantaneous speed The calculation formula is as follows:

[0114] ;in, For rotation period The bond phase zero point corresponding to the next rotation cycle passes through the time. For rotation period The time when the bond phase zero point corresponding to the previous rotation cycle passes through.

[0115] The instantaneous rotational speed sequence of this invention is set to zero at the beginning and end, and other sequence elements are calculated using specific formulas. This reduces the sensitivity of rotational speed estimation accuracy to noise without adding additional sensors, while ensuring the smoothness of the rotational speed curve endpoints and avoiding the amplification effect of boundary abrupt changes on deformation estimation.

[0116] In this embodiment, S700, based on structural constraint parameters, calculates the impeller rotation phase when each blade reaches its respective blade sensor.

[0117] In this embodiment, S800 includes the following steps:

[0118] S810. Using the key phase zero-point passage time sequence and its corresponding impeller rotation phase sequence as interpolation conditions, and the first derivative of the boundary points of the impeller rotation phase sequence as boundary conditions, construct the first interpolation equation to solve for the theoretical time when the blade passes through the corresponding blade sensor; preferably, the key phase zero-point passage time sequence is selected as... , , ..., Where RN is the number of impeller rotation cycles; in this embodiment, the corresponding boundary condition is: , ;

[0119] S820. Using the key phase zero-point passage time sequence and its corresponding instantaneous rotational speed sequence as interpolation conditions, and the first derivative of the instantaneous rotational speed sequence boundary points as boundary conditions, a second interpolation equation is constructed to solve for the instantaneous rotational speed of the blade passing the corresponding blade sensor; in this embodiment, the corresponding boundary conditions are: , .

[0120] This invention uses three sets of sequences—bond phase, rotational speed, and phase—as interpolation conditions and employs first-order derivative boundary constraints to ensure that theoretical throughput and instantaneous rotational speed are solved simultaneously within the same cubic spline framework. This guarantees that the two are strictly conjugate mathematically and eliminates the cumulative interpolation error caused by calculating rotational speed first and then time in existing technologies.

[0121] like Figure 2 As shown, the present invention also provides a blade deformation estimation system for impeller non-constant speed rotation, which includes:

[0122] The parameter acquisition module 10 is used to acquire the structural constraint parameters and evaluation strategy parameters of the impeller rotor. The structural constraint parameters include the phase sequence of the blade relative to the key phase zero point, the position sequence of the blade sensor relative to the key phase sensor, the radius of the cross-sectional circle swept by the blade tip, and the number of blades. The evaluation strategy parameters include the number of impeller rotation cycles, the sampling frequency of the key phase sensor, and the sampling frequency of the blade sensor.

[0123] The data acquisition module 20 is used to continuously acquire the response data of the key phase sensor and each blade sensor;

[0124] The key phase signal processing module 30 extracts from the key phase sensor response the number of pulses when the key phase zero point passes through each rotation cycle relative to the trigger time of the key phase zero point in the previous rotation cycle, as well as the compensation time difference until the next sampling pulse, to form a first sequence;

[0125] The blade signal processing module 40 extracts from the blade sensor response the number of pulses relative to the current rotation cycle key phase zero trigger time when each blade passes through the blade sensor, as well as the compensation time difference until the next sampling pulse, to form a second sequence;

[0126] The time calculation module 50 is used to select the bond phase time of the starting calculation of the rotation period as the reference zero time based on the first sequence and the second sequence, and to calculate the bond phase zero point passing time of other rotation periods associated with the period to be solved, as well as the measured time of each blade passing through different blade sensors.

[0127] The instantaneous speed sequence generation module 60 is used to generate the instantaneous speed sequence of the impeller based on the key phase zero-point passage time sequence;

[0128] The impeller rotation phase calculation module 70 is used to calculate the impeller rotation phase when each blade reaches its respective blade sensor based on structural constraint parameters.

[0129] The interpolation solution module 80 is used to solve the theoretical time and corresponding instantaneous speed of the blade passing through the corresponding blade sensor for different rotation cycles based on the key phase zero-point passage time sequence, instantaneous speed sequence and impeller rotation phase sequence, using the interpolation method.

[0130] The deformation estimation module 90 is used to estimate the blade deformation based on the difference between the measured time and the theoretical time, the radius of the cross-sectional circle swept by the blade tip, and the instantaneous rotational speed.

[0131] This invention also provides a computer-readable storage medium, which may be a computer-readable storage medium included in the memory described in the above embodiments; or it may be a standalone computer-readable storage medium not assembled into a device. The computer-readable storage medium stores at least one instruction, which is loaded and executed by a processor to implement... Figure 1The method for estimating blade deformation under non-constant speed rotation of the impeller is shown. The computer-readable storage medium can be a read-only memory, a hard disk, or an optical disk, etc.

[0132] It should be noted that the various embodiments in this specification are described in a progressive manner, with each embodiment focusing on the differences from other embodiments. Similar or identical parts between embodiments can be referred to interchangeably. For system embodiments and storage medium embodiments, since they are basically similar to method embodiments, the descriptions are relatively simple; relevant parts can be referred to the descriptions in the method embodiments.

[0133] Furthermore, in this document, the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such a process, method, article, or apparatus. Without further limitation, an element defined by the phrase "comprising one..." does not exclude the presence of other identical elements in the process, method, article, or apparatus that includes said element.

[0134] The foregoing description illustrates and describes preferred embodiments of the present invention. It should be understood that the present invention is not limited to the forms disclosed herein and should not be construed as excluding other embodiments. It can be used in various other combinations, modifications, and environments, and can be altered within the scope of the inventive concept by means of the foregoing teachings or techniques or knowledge in related fields. Any modifications and variations made by those skilled in the art that do not depart from the spirit and scope of the present invention should be within the protection scope of the appended claims.

Claims

1. A method for estimating blade deformation under non-constant speed rotation of an impeller, characterized by, The method comprises the following steps: S100, obtaining structure constraint parameters and evaluation strategy parameters of the impeller rotor; The structure constraint parameters comprise a phase sequence of a blade relative to a key phase zero point, a position sequence of a blade sensor relative to a key phase sensor, a radius of a tip sweep cross section circle, and a number of blades; the evaluation strategy parameters comprise a number of impeller rotation periods, a key phase sensor sampling frequency, and a blade sensor sampling frequency; S200, continuously collecting response data of the key phase sensor and each blade sensor; S300, extracting, from the key phase sensor response, a first sequence comprising a number of key phase sampling pulse counts relative to a key phase zero point passing time of each rotation period, and a compensation time difference from the key phase zero point passing time to a next key phase sampling pulse; S400, extracting, from the blade sensor response, a second sequence comprising a number of blade sampling pulse counts relative to a key phase zero point passing time of a current rotation period, and a compensation time difference from the key phase zero point passing time to a next blade sampling pulse; S500, based on the first sequence and the second sequence, selecting a key phase time of a starting calculation rotation period as a reference zero time, and calculating a key phase zero point passing time of other rotation periods associated with a to-be-solved period and a measured time of each blade passing through different blade sensors; S600, based on the key phase zero point passing time sequence, generating an instantaneous speed sequence of the impeller; S700, based on the structure constraint parameters, calculating an impeller rotation phase of each blade reaching each blade sensor; S800, based on the key phase zero point passing time sequence, the instantaneous speed sequence, and the impeller rotation phase sequence, using an interpolation method to respectively solve a theoretical time of each blade passing through a corresponding blade sensor and a corresponding instantaneous speed in different rotation periods; S900, based on a difference between the measured time and the theoretical time, the radius of the tip sweep cross section circle, and the instantaneous speed, estimating a blade deformation variable.

2. The method according to claim 1, wherein The blade sensor sampling frequency is an integer multiple of the key phase sensor sampling frequency.

3. The method according to claim 1, wherein S300 comprises the following steps: S310, for each rotation period, identifying a response pulse generated when the key phase sensor captures the key phase zero point; S320, recording a number of sampling pulse counts of the response pulse relative to a key phase zero point passing time of a previous rotation period and based on the key phase sensor sampling frequency; S330, recording a compensation time difference from a time when the key phase sensor passes through the key phase zero point to a time when a next key phase sensor sampling pulse arrives; S340, based on the number of pulse counts of S320 and the compensation time difference of S330, forming the first sequence.

4. The method according to claim 1, wherein S400 comprises the following steps: S410, for each blade sensor, sequentially identifying a response pulse generated when each blade passes through the blade sensor; S420, recording a number of sampling pulse counts of the response pulse relative to a key phase zero point passing time of a current rotation period and based on the blade sensor sampling frequency; S430, recording a compensation time difference from a time when the blade passes through the blade sensor to a time when a next blade sensor sampling pulse arrives; S440, based on the number of pulse counts of S420 and the compensation time difference of S430, forming the second sequence.

5. The method of claim 1, wherein In S500, a start calculation rotation period RS is selected, and a reference zero time is set as a key phase time of the start calculation rotation period RS. For any other rotation period after the start calculation rotation period RS The key phase zero point passes through the time The corresponding calculation formula is as follows: ; wherein, is a rotation period The number of pulse counts corresponding to the keyphasor sensor sampling frequency SK is collected by the keyphasor sensor twice when passing through the keyphasor zero point within the next rotation period; is a rotation period index; is the compensation time difference of the current period keyphasor sensor passing through the keyphasor zero point to the next keyphasor sensor sampling pulse response; is the compensation time difference of the period RS keyphasor sensor passing through the keyphasor zero point to the next keyphasor sensor sampling pulse response; Or, for any other rotation period before the start of the calculation of the rotation period RS The key phase zero point passes through the time The corresponding calculation formula is as follows: ; wherein, is a rotation period is the number of pulse counts corresponding to the sampling frequency of the keyphasor sensor SK when the keyphasor sensor passes the keyphasor zero point twice within the next rotation period; is a rotation period index; is a period is the compensated time difference of the keyphasor sensor passing the keyphasor zero point to the next keyphasor sensor sampling pulse response; is the compensated time difference of the keyphasor sensor passing the keyphasor zero point to the next keyphasor sensor sampling pulse response.

6. The method according to claim 5, wherein The key phase moment of the starting calculation rotation period RS is taken as the reference zero moment, and the i-th blade sensor responds in any other rotation period after the starting calculation rotation period RS of the impeller The measured moment of the p-th detection of the blade passing through the i-th blade sensor The corresponding calculation formula is as follows: ; wherein BN is the number of blades; is the time when the pth blade passes through the ith blade sensor, and the ith blade sensor is relative to the first blade sensor; is the number of pulse counts collected by the blade sampling frequency SB when the pth blade passes through the key phase zero point; is the blade correction time value corresponding to the time when the pth blade passes through the ith blade sensor. Or, the i-th vane sensor response in any other rotation period before the impeller start calculation rotation period RS The measured time when the p-th detection of the vane passing through the i-th vane sensor The corresponding calculation formula is as follows: ; wherein BN is the number of blades; is the time when the pth blade passes the ith blade sensor, and the ith blade sensor is relative to the first blade sensor; is the number of pulse counts collected by the blade sampling frequency SB when the pth blade passes the key phase zero point; is the blade correction time value corresponding to the time when the blade passes the ith blade sensor.

7. The method according to claim 6, wherein In the instantaneous speed sequence, the first and last elements are set to 0; For the initial calculation of the rotation period RS, the key phase zero point passes through the moment corresponding to the instantaneous speed The calculation formula is as follows: ; wherein is the time of passage of the key phase zero corresponding to the next rotation period of the start of the calculation of the rotation period RS, is the time of passage of the key phase zero corresponding to the previous rotation period of the start of the calculation of the rotation period RS. For the rotation period The key phase zero point passes through the time The corresponding instantaneous rotation speed The calculation formula is as follows: ; wherein is the time of passage of the key phase zero corresponding to the preceding rotation period of the rotation period ; and is the time of passage of the key phase zero corresponding to the preceding rotation period of the rotation period ; and For the rotation period The key phase zero point passes through the time The corresponding instantaneous rotation speed The calculation formula is as follows: ; wherein is the time of passage of the key phase zero corresponding to the preceding rotation period of the rotation period is the time of passage of the key phase zero corresponding to the preceding rotation period of the rotation period is the time of passage of the key phase zero corresponding to the preceding rotation period of the rotation period is the time of passage of the key phase zero corresponding to the preceding rotation period of the rotation period 8. The method according to claim 1, wherein S800 comprises the following steps: S810, taking the key phase zero point passing time sequence and the corresponding impeller rotation phase sequence as interpolation conditions, and taking the first derivative of the key phase zero point passing time sequence boundary point as a boundary condition, a first interpolation equation is constructed to solve the theoretical time of the blade passing through the corresponding blade sensor; S820, taking the key phase zero point passing time sequence and the corresponding instantaneous speed sequence as interpolation conditions, and taking the first derivative of the instantaneous speed sequence boundary point as a boundary condition, a second interpolation equation is constructed to solve the instantaneous speed of the blade passing through the corresponding blade sensor.

9. A blade deformation estimation system for a non-constant speed rotating impeller, characterized by, Comprise: Parameter acquisition module, used for acquiring structural constraint parameters and evaluation strategy parameters of the impeller rotor; The structural constraint parameters include the phase sequence of the blade relative to the key phase zero point, the position sequence of the blade sensor relative to the key phase sensor, the radius of the blade tip sweeping the cross-sectional circle, and the number of blades; the evaluation strategy parameters include the number of impeller rotation periods, the key phase sensor sampling frequency, and the blade sensor sampling frequency; Data acquisition module, used for continuously collecting response data of the key phase sensor and each blade sensor; The key phase signal processing module extracts the number of key phase sampling pulse counts relative to the key phase zero point trigger time of the previous rotation period and the compensation time difference to the next key phase sampling pulse when the key phase zero point passes through each rotation period from the key phase sensor response, forming a first sequence; The blade signal processing module extracts the number of blade sampling pulse counts relative to the blade passing through the blade sensor from the blade sensor response, and the compensation time difference to the next blade sampling pulse, forming a second sequence; The time calculation module is used for selecting the key phase time of the starting calculation rotation period as the reference zero time based on the first sequence and the second sequence, calculating the key phase zero point passing time of other rotation periods associated with the to-be-solved period and the measured time of each blade passing through different blade sensors; The instantaneous speed sequence generation module is used for generating the instantaneous speed sequence of the impeller based on the key phase zero point passing time sequence; The impeller rotation phase calculation module is used for calculating the impeller rotation phase of each blade reaching each blade sensor based on the structural constraint parameters; The interpolation solving module is used for solving the theoretical time and the corresponding instantaneous speed of the blade passing through the corresponding blade sensor in different rotation periods using interpolation method based on the key phase zero point passing time sequence, the instantaneous speed sequence, and the impeller rotation phase sequence; The deformation estimation module is used for estimating the blade deformation based on the difference between the measured time and the theoretical time, the radius of the blade tip sweeping the cross-sectional circle, and the instantaneous speed.

10. A computer-readable storage medium, characterized in that, The computer readable storage medium stores an impeller non-constant speed rotation blade deformation estimation program, and the impeller non-constant speed rotation blade deformation estimation program is executed by the processor to realize the steps of the impeller non-constant speed rotation blade deformation estimation method according to any one of claims 1 to 8.

Images