A method, system, and storage medium for optimizing the layout of blade vibration monitoring sensors.

Abstract

A method, system, and storage medium for optimizing the layout of blade vibration monitoring sensors include: S100, obtaining geometric constraints and vibration assessment strategy constraints for the sensor layout; S200, generating sensor phase sequences based on geometric constraints; S300, determining multiple vibration orders to be analyzed; S400, for each vibration order, constructing a coefficient matrix for vibration response estimation based on the sensor phase sequence and the total number of cycles to be calculated; S500, for each vibration order, calculating the condition number of the product of the coefficient matrix and its transpose; S600, selecting the maximum value of the condition number as the evaluation parameter of the current sensor phase sequence; S700, adjusting the sensor phase sequence based on the phase adjustment step size under geometric constraints, and re-executing S300-S600; S800, selecting the minimum value among all evaluation parameters and outputting its corresponding sensor phase sequence as the optimization result.

Description

Technical Field

[0001] This invention relates to the field of blade vibration monitoring optimization technology, and in particular to a method, system and storage medium for optimizing the layout of blade vibration monitoring sensors. Background Technology

[0002] As a core component of high-speed rotating machinery such as aero-engines and gas turbines, the vibration state of blades directly affects the safety and reliability of the entire machine. Existing blade vibration measurement methods mainly infer the blade's vibration deformation by analyzing the difference between the actual and theoretical time of arrival at each sensor. Currently, it is generally assumed that the influence of sensor installation position on vibration parameter estimation is negligible, or that sensors are arranged at equal intervals.

[0003] However, in practical applications, blade vibration response fitting is essentially a parameter inversion problem under undersampling conditions, and the sensor placement directly affects the numerical characteristics of the coefficient matrix of the vibration response equation. An inappropriate sensor layout will cause the vibration response equation to exhibit ill-conditioned characteristics, making the estimated results of key parameters such as vibration amplitude and frequency highly sensitive to measurement errors, severely reducing the reliability of vibration assessment. Although there are existing methods that attempt to improve measurement results by increasing the number of sensors, this not only increases system complexity and cost but is also difficult to implement in scenarios with limited structural space.

[0004] Therefore, how to improve the numerical stability of vibration parameter estimation under the condition of limited number of sensors has become a key technical problem that urgently needs to be solved in this field. Summary of the Invention

[0005] To address the aforementioned technical problems, this invention provides a method for optimizing the layout of blade vibration monitoring sensors, comprising the following steps:

[0006] S100. Obtain the geometric constraints of the sensor layout and the strategy constraints of vibration assessment; the geometric constraints include the minimum phase interval between sensors, the total number of sensors, and prohibited installation areas; the strategy constraints include the total number of cycles to be calculated, the order range of vibration fitting, and the phase adjustment step size for local search.

[0007] S200. Generate a sensor phase sequence based on the geometric constraints;

[0008] S300. Based on the total number and order range of the periods to be calculated, determine multiple vibration orders to be analyzed;

[0009] S400. For each vibration order, a coefficient matrix for vibration response estimation is constructed based on the sensor phase sequence and the total number of periods to be calculated.

[0010] S500. For each vibration order, calculate the condition number of the product of the coefficient matrix and its transpose.

[0011] S600. Select the maximum value of the condition number among all vibration orders as the evaluation parameter of the current sensor phase sequence.

[0012] S700: Under the geometric constraints, adjust the sensor phase sequence based on the phase adjustment step size, and re-execute S300-S600 based on each adjusted sensor phase sequence;

[0013] S800: Select the minimum value among all evaluation parameters and output the corresponding sensor phase sequence as the optimization result.

[0014] Optionally, the prohibited installation area is determined based on interference from the housing structure or interference from adjacent components; the minimum phase interval is determined based on the physical dimensions of the sensor and the strength requirements of the housing structure.

[0015] Optionally, the minimum value in the order range is greater than the reciprocal of the total number of periods to be calculated.

[0016] Optionally, the phase adjustment step size is selected based on the following conditions:

[0017] For any natural number, the first product of the natural number and the phase adjustment step size is greater than 0 and less than the second product of the total number of cycles to be calculated and 360; the ratio of the first product value to 360 is not an integer.

[0018] Optionally, S300 may determine multiple vibration orders to be analyzed based on the following first preset formula;

[0019] ;

[0020] in, Let k be the minimum value within the range of orders, k be a positive integer, and C be the total number of periods to be calculated. Let k be the vibration order to be analyzed.

[0021] Optionally, S400 constructs the coefficient matrix for vibration response estimation based on the following second preset formula. :

[0022] ;

[0023] in, Let l be the k-th vibration order to be analyzed, and l be the current calculation period number. Let be the angle of the i-th sensor.

[0024] Optionally, the S500 includes the following steps:

[0025] S510. Calculate the product matrix of the coefficient matrix and its transpose;

[0026] S520. Calculate the condition number of the coefficient matrix.

[0027] Optionally, under the geometric constraints, the sensor phase sequence is adjusted based on the phase adjustment step size, specifically as follows:

[0028] The phase angle of the sensor phase sequence is adjusted sequentially based on the phase adjustment step size, and the adjusted sensor phase sequence is verified to meet the geometric constraints after each adjustment.

[0029] Corresponding to the blade vibration monitoring sensor layout optimization method, the present invention provides a blade vibration monitoring sensor layout optimization system, which includes:

[0030] The constraint acquisition module is used to acquire the geometric constraints of the sensor layout and the strategy constraints of vibration assessment. The geometric constraints include the minimum phase interval between sensors, the total number of sensors, and prohibited installation areas. The strategy constraints include the total number of cycles to be calculated, the order range of vibration fitting, and the phase adjustment step size for local search.

[0031] A sensor phase sequence generation module is used to generate a sensor phase sequence based on the geometric constraints.

[0032] The vibration order generation module is used to determine multiple vibration orders to be analyzed based on the total number and order range of the periods to be calculated.

[0033] The coefficient matrix construction module is used to construct the coefficient matrix for vibration response estimation for each vibration order, based on the sensor phase sequence and the total number of periods to be calculated.

[0034] The calculation module is used to calculate the condition number of the product of the coefficient matrix and its transpose for each vibration order.

[0035] The evaluation parameter selection module is used to select the maximum value of the condition number among all vibration orders as the evaluation parameter of the current sensor phase sequence.

[0036] The dynamic adjustment module is used to adjust the sensor phase sequence based on the phase adjustment step size under the geometric constraints, and to call the vibration order generation module, coefficient matrix construction module, calculation module and evaluation parameter selection module to re-execute the corresponding steps based on each adjusted sensor phase sequence.

[0037] The optimization result output module is used to select the minimum value among all evaluation parameters and output the corresponding sensor phase sequence as the optimization result.

[0038] In addition, to achieve the above objectives, the present invention also provides a computer-readable storage medium storing a blade vibration monitoring sensor layout optimization program, which, when executed by a processor, implements the steps of the blade vibration monitoring sensor layout optimization method described above.

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

[0040] (1) By obtaining the geometric constraints of the sensor layout and the strategy constraints of vibration assessment, the problem of the vibration response equation exhibiting ill-conditioned characteristics due to unreasonable sensor placement in the prior art is solved; by generating the sensor phase sequence based on geometric constraints and determining multiple vibration orders to be analyzed, the coefficient matrix of vibration response estimation is constructed, and the condition number is calculated (to evaluate the numerical stability of parameter estimation under the corresponding vibration order), thereby selecting the maximum value of the condition number as the evaluation parameter, thus realizing the optimization of the sensor layout and improving the numerical stability of vibration parameter estimation.

[0041] (2) By prohibiting installation areas and minimizing phase intervals, the installation position of the sensor is made more reasonable, avoiding interference from the housing structure or adjacent components, and improving the practicality and reliability of the sensor layout.

[0042] (3) By setting the minimum value in the order range to be greater than the reciprocal of the total number of cycles to be calculated, the selection of vibration order is optimized, the accuracy of vibration response estimation is improved, and thus the overall performance of vibration parameter estimation is enhanced.

[0043] (4) By reasonably selecting the phase adjustment step size, it is ensured that the same equivalent phase point will not be checked repeatedly in different periods during local search, thus avoiding "spinning around in place". In this way, search efficiency is guaranteed while also ensuring the accuracy of search results.

[0044] (5) By using the first preset formula, multiple vibration orders to be analyzed are determined, making the selection of vibration orders more scientific and reasonable, which helps to improve the accuracy of vibration response estimation.

[0045] (6) Constructing the coefficient matrix of vibration response estimation through the second preset formula helps to improve the accuracy and reliability of vibration parameter estimation.

[0046] (7) By calculating the condition number of the product of the coefficient matrix and its transpose, a quantitative method is provided for evaluating the sensor phase sequence, which helps to more accurately evaluate and select the optimal sensor layout.

[0047] (8) By adjusting the sensor phase sequence under geometric constraints and verifying whether the adjusted sequence meets the geometric constraints, the optimization process of sensor layout takes into account both performance and actual installation limitations, thus improving the feasibility of the layout scheme. Attached Figure Description

[0048] 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:

[0049] Figure 1 This is a simplified flowchart of an embodiment of the blade vibration monitoring sensor layout optimization method of the present invention;

[0050] Figure 2 This is a schematic diagram showing the coordinate system definition of the sensor in the housing in one embodiment of the blade vibration monitoring sensor layout optimization method of the present invention;

[0051] Figure 3 This is a schematic diagram of blade vibration displacement fitting in an embodiment of the blade vibration monitoring sensor layout optimization method of the present invention;

[0052] Figure 4 This is a framework diagram of an embodiment of the blade vibration monitoring sensor layout optimization system of the present invention. Detailed Implementation

[0053] 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.

[0054] The impeller of the equipment typically has multiple blades. The numbering of each blade on the same cross-section is determined by its position relative to the key phase point (the edge of the groove or protrusion monitored by the key phase sensor) and the rotor's rotation direction. The location of the key phase point is the reference zero phase of the polar coordinate system formed by the projection point of the axis of rotation into the rotor cross-section. If the reverse rotation direction is defined as the positive direction, the incremental order of the blade numbers is determined by the magnitude of the blade's phase in the polar coordinate system. For example, the phase of blade j relative to the reference zero point is... The larger the blade number, the larger the corresponding phase value. Please refer to... Figure 2As shown, the coordinate system takes the impeller rotor rotation direction as the main reference and the housing machining reference as the zero phase direction, which helps to ensure the position accuracy of the sensor.

[0055] In vibration measurement and evaluation, if the torsional deformation of the shaft is ignored, the difference between the arrival time of each blade tip to the sensor and its theoretical value is mainly caused by the deformation of the blades during equipment operation. This can be further divided into two parts: static deformation caused by average load and vibration deformation caused by dynamic load. The theoretical time t from each blade tip to the sensor is determined by the angle between the sensor and the key phase sensor. Phase of the blade relative to the keyway groove and the speed of the impeller rotor Decide.

[0056] The blade deformation P can be estimated from the linear velocity of the section at the blade tip, based on the time difference τ between the time from the blade tip to each sensor and the theoretical value. For a specific blade j, the deformation P passing through the i-th sensor is... Represented as:

[0057] ;

[0058] in, The leaf tip sweeps across the radius of the circle; Let be the time difference between the j-th blade and the i-th sensor.

[0059] From the perspective of blade deformation combination It can also be expressed as:

[0060] ;

[0061] in, For the j-th blade at an excitation frequency of At that time, the estimated value of the blade vibration amplitude; the excitation frequency is the estimated value of the blade vibration frequency; Let be the angle between the i-th sensor and the key phase sensor along the rotor rotation direction (i.e., the sensor's mounting position). 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; The vibration frequency is When, the estimated value of the corresponding static deformation component of the blade.

[0062] When estimating blade vibration Represented as The order component is given by k, where k is the index of the vibration order to be analyzed; assuming the vibration frequency is the m-th order of the rotational frequency, that is... ; Estimates of blade vibration amplitude ,Apart from In addition, it is also affected Influence.

[0063] To facilitate understanding of the present invention, a schematic diagram of blade vibration displacement fitting is also provided; please refer to it. Figure 3 This refers to the deformation measured by four sensors for a specific blade j. And the vibration response waveform obtained by fitting the responses of the four sensors.

[0064] Considering In the vibration response equation, and This is the focus of the evaluation. Related to the phase zero of the timing reference, and It mainly involves static deformation components. To ensure... and To estimate accuracy, this invention provides a method for optimizing the layout of blade vibration monitoring sensors, which includes the following steps:

[0065] S100. Obtain the geometric constraints of the sensor layout and the strategy constraints of vibration assessment; the geometric constraints include the minimum phase interval between sensors, the total number of sensors, and the prohibited installation area; the strategy constraints include the total number of cycles to be calculated, the order range of vibration fitting, and the phase adjustment step size for local search; preferably, the angle range of the prohibited installation area includes the corresponding start angle and end angle, and the calculation needs to cross these unreasonable installation positions.

[0066] S200, Generate sensor phase sequence based on geometric constraints, for example Where N is the total number of sensors. to This represents the first to the Nth sensors;

[0067] S300. Based on the total number of cycles to be calculated and the range of orders (including maximum and minimum values), determine multiple vibration orders to be analyzed.

[0068] S400. For each vibration order, construct the coefficient matrix for vibration response estimation based on the sensor phase sequence and the total number of cycles to be calculated.

[0069] S500. For each vibration order, calculate the condition number of the product of the coefficient matrix and its transpose.

[0070] S600. Select the maximum value of the condition number among all vibration orders as the evaluation parameter of the current sensor phase sequence.

[0071] S700. Under geometric constraints, adjust the sensor phase sequence based on the phase adjustment step size, and re-execute S300-S600 based on each adjusted sensor phase sequence to obtain the corresponding evaluation parameters.

[0072] S800: Select the minimum value among all evaluation parameters and output the corresponding sensor phase sequence as the optimization result.

[0073] This invention addresses the problem of ill-conditioned vibration response equations caused by unreasonable sensor placement in existing technologies by acquiring geometric constraints of sensor layout and strategic constraints of vibration assessment. It generates sensor phase sequences based on geometric constraints, determines multiple vibration orders to be analyzed, constructs a coefficient matrix for vibration response estimation, and calculates the condition number (to evaluate the numerical stability of parameter estimation under the corresponding vibration order). The maximum value of the condition number is then selected as the evaluation parameter, thereby optimizing the sensor layout and improving the numerical stability of vibration parameter estimation.

[0074] In this embodiment, the prohibited installation area is determined based on interference from the housing structure or interference from adjacent components; the minimum phase interval is determined based on the physical dimensions of the sensor and the strength requirements of the housing structure.

[0075] This invention, by prohibiting installation areas and minimizing phase intervals, makes the sensor installation position more reasonable, avoids problems of interference from the housing structure or adjacent components, and improves the practicality and reliability of the sensor layout.

[0076] In this embodiment, the minimum value in the order range is greater than the reciprocal of the total number of cycles to be calculated. That is, under the minimum fitting order, it is ensured that the phase swept by all blade sensors during the analysis covers the vibration response to be fitted for more than one cycle. The maximum and minimum values ​​of the order range are determined according to the design parameters and analysis requirements.

[0077] This invention optimizes the selection of vibration order by setting the minimum value in the order range to be greater than the reciprocal of the total number of cycles to be calculated, thereby improving the accuracy of vibration response estimation and enhancing the overall performance of vibration parameter estimation.

[0078] In this embodiment, the phase adjustment step size is selected based on the following conditions:

[0079] For any natural number, the first product of the natural number and the phase adjustment step size is greater than 0 and less than the second product of the total number of cycles to be calculated and 360; the ratio of the first product value to 360 is not an integer.

[0080] Preferably, the phase adjustment step size satisfies the following formula:

[0081] N×C×S<360 <N×C×(S+1);

[0082] Where C is the total number of cycles to be calculated, and S is the phase adjustment step size.

[0083] By rationally selecting the phase adjustment step size, this invention ensures that the same equivalent phase point is not repeatedly checked in different periods during local search, avoiding "spinning in place" and thus ensuring both search efficiency and the accuracy of search results.

[0084] In this embodiment, S300 specifically determines multiple vibration orders to be analyzed according to the following first preset formula;

[0085] ;

[0086] in, The minimum value in the order range, where k is a positive integer, preferably 0 < k < , Let k be the vibration order to be analyzed.

[0087] This invention determines multiple vibration orders to be analyzed through a first preset formula, making the selection of vibration orders more scientific and reasonable, and helping to improve the accuracy of vibration response estimation.

[0088] In this embodiment, S400 specifically constructs the coefficient matrix for vibration response estimation according to the following second preset formula. :

[0089] ;

[0090] in, Let l represent the k-th vibration order to be analyzed, and l be the current calculation cycle number in the iteration process from 1 to C when forming the coefficient matrix. Let be the angle of the i-th sensor.

[0091] The present invention constructs a coefficient matrix for vibration response estimation using a second preset formula, which helps to improve the accuracy and reliability of vibration parameter estimation.

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

[0093] S510. Calculate the product matrix of the coefficient matrix and its transpose;

[0094] S520, Calculate the condition number of the coefficient matrix.

[0095] This invention provides a quantitative method for evaluating sensor phase sequences by calculating the condition number of the product of the coefficient matrix and its transpose, which helps to more accurately assess and select the optimal sensor layout.

[0096] In this embodiment, under geometric constraints, the sensor phase sequence is adjusted based on the phase adjustment step size, specifically as follows:

[0097] The phase angle of the sensor phase sequence is adjusted by sequentially increasing the phase adjustment step size, and the sensor phase sequence is verified to meet the geometric constraints after each adjustment.

[0098] This invention ensures that the sensor layout optimization process takes into account both performance and practical installation limitations by adjusting the sensor phase sequence under geometric constraints and verifying whether the adjusted sequence meets the geometric constraints, thereby improving the feasibility of the layout scheme.

[0099] like Figure 4 As shown, the present invention also provides a blade vibration monitoring sensor layout optimization system, which includes:

[0100] The constraint acquisition module 10 is used to acquire the geometric constraints of the sensor layout and the strategy constraints of vibration assessment. The geometric constraints include the minimum phase interval between sensors, the total number of sensors, and the prohibited installation area. The strategy constraints include the total number of cycles to be calculated, the order range of vibration fitting, and the phase adjustment step size for local search.

[0101] Sensor phase sequence generation module 20 is used to generate sensor phase sequences based on geometric constraints;

[0102] The vibration order generation module 30 is used to determine multiple vibration orders to be analyzed based on the total number of cycles to be calculated and the order range.

[0103] The coefficient matrix construction module 40 is used to construct the coefficient matrix for vibration response estimation for each vibration order, based on the sensor phase sequence and the total number of cycles to be calculated.

[0104] The calculation module 50 is used to calculate the condition number of the product of the coefficient matrix and its transpose for each vibration order.

[0105] The evaluation parameter selection module 60 is used to select the maximum value of the condition number among all vibration orders as the evaluation parameter of the current sensor phase sequence.

[0106] The dynamic adjustment module 70 is used to adjust the sensor phase sequence based on the phase adjustment step size under geometric constraints, and to call the vibration order generation module, coefficient matrix construction module, calculation module and evaluation parameter selection module to re-execute the corresponding steps based on each adjusted sensor phase sequence.

[0107] The optimization result output module 80 is used to select the minimum value among all evaluation parameters and output the corresponding sensor phase sequence as the optimization result.

[0108] 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 1 The method for optimizing the layout of blade vibration monitoring sensors is shown. The computer-readable storage medium can be a read-only memory, a hard disk, or an optical disk, etc.

[0109] 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.

[0110] 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.

[0111] 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 of optimizing a blade vibration monitoring sensor layout, characterized by, The method comprises the following steps: S100, acquiring geometric constraints of a sensor layout and strategy constraints of vibration evaluation; the geometric constraints comprise a minimum phase interval between sensors, a total number of sensors, and a prohibited installation area; the strategy constraints comprise a total number of to-be-calculated periods, a vibration fitting order range, and a phase adjustment step length for local search; S200, generating a sensor phase sequence based on the geometric constraints; S300, determining a plurality of vibration orders to be analyzed according to the total number of to-be-calculated periods and the order range; The plurality of vibration orders to be analyzed are specifically determined according to the following first preset formula: ; wherein is the minimum value in the order range, k is a positive integer, C is the total number of periods to be calculated, is the kth vibration order to be analyzed; S400, for each vibration order, based on the sensor phase sequence and the total number of periods to be calculated, a coefficient matrix of vibration response estimation is constructed according to the following second preset formula : ; wherein, I is a current calculation cycle number, is an angle of the i-th sensor; S500, for each vibration order, calculating a condition number of a product of the coefficient matrix and its transpose matrix; S600, selecting a maximum value of the condition numbers of all vibration orders as an evaluation parameter of a current sensor phase sequence; S700, under the geometric constraints, adjusting the sensor phase sequence based on the phase adjustment step length, and re-executing S300-S600 based on each adjusted sensor phase sequence; S800, selecting a minimum value of all evaluation parameters, and outputting a sensor phase sequence corresponding to the minimum value as an optimization result.

2. The blade vibration monitoring sensor layout optimization method of claim 1, wherein, The prohibited installation area is determined according to a shell structure interference or an adjacent component interference; and the minimum phase interval is determined according to physical dimensions of the sensors and a shell structure strength requirement.

3. The blade vibration monitoring sensor layout optimization method of claim 1, wherein, A minimum value in the order range is greater than a reciprocal of the total number of to-be-calculated periods.

4. The blade vibration monitoring sensor layout optimization method of claim 1, wherein, The phase adjustment step length is selected according to the following conditions: For any natural number, a first product value of the natural number and the phase adjustment step length is greater than 0 and less than a second product value of the total number of to-be-calculated periods and 360; A ratio of the first product value to 360 is a non-integer.

5. The blade vibration monitoring sensor layout optimization method of claim 1, wherein, S500 comprises the following steps: S510, calculating a product matrix of the coefficient matrix and its transpose matrix; S520, calculating a condition number of the coefficient matrix.

6. The blade vibration monitoring sensor layout optimization method of claim 1, wherein, Under the geometric constraints, the sensor phase sequence is adjusted based on the phase adjustment step length, specifically: The phase angle of the sensor phase sequence is sequentially and incrementally adjusted based on the phase adjustment step length, and whether the adjusted sensor phase sequence satisfies the geometric constraints is verified after each adjustment.

7. A blade vibration monitoring sensor layout optimization system, characterized by, The method comprises the following steps: A constraint condition acquisition module is configured to acquire geometric constraints of a sensor layout and strategy constraints of vibration evaluation; the geometric constraints comprise a minimum phase interval between sensors, a total number of sensors, and a prohibited installation area; the strategy constraints comprise a total number of to-be-calculated periods, a vibration fitting order range, and a phase adjustment step length for local search; A sensor phase sequence generation module is configured to generate a sensor phase sequence based on the geometric constraints; A vibration order generation module is configured to determine a plurality of vibration orders to be analyzed according to the total number of to-be-calculated periods and the order range; the plurality of vibration orders to be analyzed are specifically determined according to the following first preset formula: ; wherein is the minimum value in the order range, k is a positive integer, C is the total number of cycles to be calculated, is the kth vibration order to be analyzed; A coefficient matrix construction module is configured to, for each vibration order, construct a coefficient matrix of vibration response estimation according to a second preset formula based on the sensor phase sequence and the total number of periods to be calculated : ; wherein, I is a current calculation cycle number, is an angle of the i-th sensor; A calculation module is configured to, for each vibration order, calculate a condition number of a product of the coefficient matrix and its transpose matrix; An evaluation parameter selection module is configured to select a maximum value of the condition numbers of all vibration orders as an evaluation parameter of a current sensor phase sequence; The dynamic adjustment module is configured to adjust the sensor phase sequence based on the phase adjustment step under the geometric constraint, and based on each adjusted sensor phase sequence, call the vibration order generation module, the coefficient matrix construction module, the calculation module, and the evaluation parameter selection module to re-execute the corresponding steps. The optimization result output module is configured to select the minimum value among all the evaluation parameters, and output the corresponding sensor phase sequence as the optimization result.

8. A computer-readable storage medium, characterized in that, The computer readable storage medium stores a blade vibration monitoring sensor layout optimization program. When the blade vibration monitoring sensor layout optimization program is executed by the processor, the steps of the blade vibration monitoring sensor layout optimization method according to any one of claims 1 to 6 are implemented.

Images