CN117030270A - Method for evaluating time constant of T25 sensor based on airborne measurable parameters - Google Patents
Method for evaluating time constant of T25 sensor based on airborne measurable parameters Download PDFInfo
- Publication number
- CN117030270A CN117030270A CN202311006445.6A CN202311006445A CN117030270A CN 117030270 A CN117030270 A CN 117030270A CN 202311006445 A CN202311006445 A CN 202311006445A CN 117030270 A CN117030270 A CN 117030270A
- Authority
- CN
- China
- Prior art keywords
- time constant
- section
- sensor
- mach number
- measurable
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 24
- 238000012360 testing method Methods 0.000 claims abstract description 15
- BASFCYQUMIYNBI-UHFFFAOYSA-N platinum Chemical compound [Pt] BASFCYQUMIYNBI-UHFFFAOYSA-N 0.000 description 8
- 230000008569 process Effects 0.000 description 8
- 238000012937 correction Methods 0.000 description 6
- 230000007704 transition Effects 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 5
- 230000008859 change Effects 0.000 description 5
- 238000006243 chemical reaction Methods 0.000 description 5
- 230000000875 corresponding effect Effects 0.000 description 5
- 238000010219 correlation analysis Methods 0.000 description 4
- 229910052697 platinum Inorganic materials 0.000 description 4
- 238000010998 test method Methods 0.000 description 4
- 238000007405 data analysis Methods 0.000 description 3
- 238000013461 design Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 230000001052 transient effect Effects 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 239000011540 sensing material Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M15/00—Testing of engines
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M9/00—Aerodynamic testing; Arrangements in or on wind tunnels
- G01M9/02—Wind tunnels
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M9/00—Aerodynamic testing; Arrangements in or on wind tunnels
- G01M9/08—Aerodynamic models
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Fluid Mechanics (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
The application belongs to the technical field of aeroengines, and particularly relates to a method for evaluating a time constant of a T25 sensor based on airborne measurable parameters, which comprises the steps of obtaining a plurality of measurable parameters affecting the Mach number of a T25 section in a full envelope range; calculating the correlation coefficient of each measurable parameter and the Mach number of the T25 section respectively, and taking the measurable parameter corresponding to the maximum correlation coefficient as the most relevant measurable parameter; fitting a fitting formula of the most relevant measurable parameters and the Mach number of the T25 section; based on the relation between the T25 section Mach number and the T25 sensor time constant, a fitting formula of the most relevant measurable parameter and the T25 sensor time constant is fitted, so that the problem that the T25 time constant cannot be obtained in the whole envelope range of the engine because the T25 section Mach number is not measurable is solved, and the problem that the T25 dynamic test data are difficult to analyze due to the influence of multiple factors is essentially avoided.
Description
Technical Field
The application belongs to the technical field of aeroengines, and particularly relates to a method for evaluating a time constant of a T25 sensor based on airborne measurable parameters.
Background
High performance aircraft engines rely on high performance components, and high pressure compressors are one of the important components. In order to achieve excellent performance in a wider range of high-pressure compressors, current engines generally employ high-pressure compressor vanes alpha 2 The adjustable mode enables the compressor to work in a smaller range of inlet attack angles. The conversion rotating speed N2r25 of the high-pressure rotor relative to the inlet temperature of the high-pressure compressor can better reflect the working condition of the compressor, the regulation rule of the guide vane of the high-pressure compressor is designed based on the change of N2r25, and the rotating speed of the high-pressure rotor is converted into a25 section by the N2r25, and the basic formula is as follows:
wherein,the time constant of the acquired value of the T25 sensor has more influence factors in the whole envelope range of the engine, particularly in the acceleration and deceleration process of the engine, and the time constant of the acquired value is more than that of the T25 sensor 25 The time constant can change rapidly, which affects the accuracy of the T25 acquisition during the transition, and thus affects the given control angle of α2.
As shown in fig. 4-5, the T25 sensor adopts a platinum resistor as a temperature sensing element, the platinum resistor changes with the temperature, the time constant T has more influence factors, and in general, when the air flow rate is low, the time constant is larger, and when the air flow rate is high, the time constant is correspondingly smaller. The wind tunnel test adopts the method that the time constant of a sensor is obtained under the measurable air flow rate, and when the engine is used in a full envelope range, the T25 section air flow rate cannot be directly measured, so that the air flow rate is difficult to directly apply.
The method acquires key influence factors of the T25 section air flow velocity based on the multi-data correlation, and finally evaluates the time constant of the T25 sensor based on the airborne measurable parameter when the engine flies within the full envelope range by combining wind tunnel test results, thereby providing an important basis for deep data analysis and sensor correction.
In the prior art, under the measurable steady-state air flow rate, a wind tunnel test method is adopted to measure the time constant of the sensor, and a method for evaluating the time constant of the T25 sensor is lacked when the air flow rate is not measurable in the dynamic change process within the full envelope range.
At present, a wind tunnel test method is adopted to measure the time constant of a sensor under the condition of stable and measurable air flow rate, but in the actual use working condition of an aeroengine, the time constant of the sensor changes along with the air flow rate of a section of the sensor, and the Mach number of a section T25 is not measurable, so that the time constant measured by the wind tunnel test cannot be directly applied, the data measured by the sensor T25 in the transient process is ambiguous as to how much difference is between the data and the real data, and further analysis and correction are difficult, and the acquisition value of the sensor T25 in the transient process directly influences the control of a guide vane of a compressor, thereby influencing the stability margin of the engine.
Disclosure of Invention
In order to solve the above problems, the present application provides a method for estimating a time constant of a T25 sensor based on an onboard measurable parameter, comprising:
acquiring a plurality of measurable parameters affecting the Mach number of the T25 section in the full envelope range;
calculating the correlation coefficient of each measurable parameter and the Mach number of the T25 section respectively, and taking the measurable parameter corresponding to the maximum correlation coefficient as the most relevant measurable parameter;
fitting a fitting formula of the most relevant measurable parameters and the Mach number of the T25 section;
and fitting a fitting formula of the most relevant measurable parameter and the time constant of the T25 sensor based on the relation between the Mach number of the T25 section and the time constant of the T25 sensor.
Preferably, the measurable parameters include: the low-pressure compressor converts the rotating speed N1r, the aircraft altitude H, the engine inlet section air flow speed Mach number Ma and the T25 section air flow speed Mach number.
Preferably, the calculation formula of the correlation coefficient is:
x is the Mach number of the T25 section, and Y is a measurable parameter; the correlation coefficient is 0.0 to +/-0.3, the real correlation is +/-0.3 to +/-0.5, and the high correlation is +/-0.8 to +/-1.
Preferably, the most relevant measurable parameter is preferably the reduced rotation speed N1r of the low-pressure compressor.
Preferably, the relationship between the Mach number of the T25 section and the time constant of the T25 sensor is obtained through wind tunnel test.
Preferably, the corresponding relation between the T25 section air flow velocity and the T25 sensor time constant is obtained through wind tunnel test, and the relation between the Mach number of the T25 section position and the T25 sensor time constant is established based on the relation between the T25 section air flow velocity and the Mach number of the T25 section position
The advantages of the application include: the application creates a method for evaluating the time constant of the T25 sensor based on the airborne measurable parameters, provides important support for further data analysis, T25 sensor correction and T25 airborne model design, and has very good engineering application value.
In the prior art, a sensor time constant is measured under the condition of stable and measurable air flow rate by adopting a wind tunnel test method, but in the actual use working condition, the sensor time constant changes along with the Mach number of the section of the sensor, and the Mach number of the section T25 is not measurable, so that the time constant measured by the wind tunnel test cannot be directly applied, the data measured by the sensor T25 in the transition process is ambiguous in terms of the difference between the data and the real data, the further analysis and correction are difficult, and the time constant of the sensor T25 reflects the response speed of the sensor, and the acquisition value directly influences the control of a compressor guide vane in the transition process, thereby influencing the stability margin of an engine.
According to the application, based on multi-data correlation analysis, key influence factors of the T25 sensor time constant in the whole envelope range can be obtained, strong correlation factors of the T25 sensor time constant in the conversion rotating speed of the low-pressure compressor are extracted, a fitting formula of the corresponding relation between the T25 sensor time constant and the conversion rotating speed N1r of the low-pressure compressor is established, the problem that the T25 time constant cannot be obtained due to the fact that the Mach number of the T25 section of the engine is not measurable in the whole envelope range is solved, and the problem that the T25 dynamic test data is difficult to analyze due to multi-factor influence is essentially avoided.
Drawings
FIG. 1 is a graph showing dependence of parameters on Mach number of section T25 according to a preferred embodiment of the present application;
FIG. 2 is a graph showing the correspondence between Mach numbers Ma25 and N1r in a section T25 according to a preferred embodiment of the present application;
FIG. 3 is a graph showing the time constant versus N1r for a T25 sensor according to a preferred embodiment of the present application;
FIG. 4 is a graph of a T25 sensor transition for different time constants;
FIG. 5T 25 sensor effect schematic for compressor vane.
Detailed Description
In order to make the technical solution of the present application and its advantages more clear, the technical solution of the present application will be further and completely described in detail with reference to the accompanying drawings, it being understood that the specific embodiments described herein are only some of the embodiments of the present application, which are for explanation of the present application and not for limitation of the present application. It should be noted that, for convenience of description, only the part related to the present application is shown in the drawings, and other related parts may refer to the general design, and the embodiments of the present application and the technical features of the embodiments may be combined with each other to obtain new embodiments without conflict.
Furthermore, unless defined otherwise, technical or scientific terms used in the description of the application should be given the ordinary meaning as understood by one of ordinary skill in the art to which the application pertains. The terms "upper," "lower," "left," "right," "center," "vertical," "horizontal," "inner," "outer," and the like as used in the description of the present application are merely used for indicating relative directions or positional relationships, and do not imply that the devices or elements must have a specific orientation, be constructed and operated in a specific orientation, and that the relative positional relationships may be changed when the absolute position of the object to be described is changed, thus not being construed as limiting the application. The terms "first," "second," "third," and the like, as used in the description of the present application, are used for descriptive purposes only and are not to be construed as indicating or implying any particular importance to the various components. The use of the terms "a," "an," or "the" and similar referents in the description of the application are not to be construed as limiting the amount absolutely, but rather as existence of at least one. As used in this description of the application, the terms "comprises," "comprising," or the like are intended to cover an element or article that appears before the term as such, but does not exclude other elements or articles from the list of elements or articles that appear after the term.
Furthermore, unless specifically stated and limited otherwise, the terms "mounted," "connected," and the like in the description of the present application are used in a broad sense, and for example, the connection may be a fixed connection, a removable connection, or an integral connection; can be mechanically or electrically connected; can be directly connected or indirectly connected through an intermediate medium, and can also be communicated with the inside of two elements, and the specific meaning of the two elements can be understood by a person skilled in the art according to specific situations.
1. Mathematical model of T25 sensor
The T25 sensor adopts a platinum resistor as a temperature sensing element, the platinum resistor is utilized to change along with the change resistance value of the temperature, the output resistance value and the temperature are in a linear relation, and the general form of a transfer function is as follows:
the magnitude of the time constant is related to the construction and characteristics of the temperature sensing material, and the external environment, and is a dynamic value. The key influencing parameter is the air flow velocity, i.e. Mach number at the T25 section position.
Table 1T25 sensor time constant wind tunnel test results (example)
Since the Mach number of the T25 section position is not measurable during the flight of the aircraft, the correlation between the T25 section Mach number and the engine measurable parameter must be found.
2. T25 section Mach number key influence factor acquisition
Because the mach number of the T25 section cannot be directly measured, first, the mach number data of the T25 section corresponding to different status points in the full envelope is obtained according to the characteristics of the compressor, for example, as follows:
TABLE 2 Mach number 25 section for different State points (example)
And carrying out correlation analysis on the Mach number of the 25 sections in the full envelope range and a plurality of groups of measurable parameters.
The correlation coefficient is the amount of correlation between study variables, and the expression for the correlation coefficients for variables X and Y is:
typically, the correlation coefficients are 0.0 to.+ -. 0.3 for micro-correlations,.+ -. 0.3 to.+ -. 0.5 for real correlations,.+ -. 0.8 to.+ -. 1 for high correlations.
According to the multi-data correlation analysis, the correlation coefficient between each measurable parameter and the Ma25 can be obtained, and the strong correlation factor can be extracted.
For example, as shown in fig. 1:
corrcoef(n1r,Ma25)=0.9972
corrcoef(T25,Ma25)=0.6872
corrcoef(Ma,Ma25)=0.2765
corrcoef(H,Ma25)=0.260
namely: the n1r and 25 section mach number correlation coefficient is 0.9972, the 25 section temperature and 25 section mach number correlation coefficient is 0.6872, the intake mach number and 25 section mach number correlation coefficient is 0.2765, and the height H and 25 section mach number correlation coefficient is 0.26.
From the above analysis, n1r is highly correlated with the 25 section Mach number, i.e., the T25 sensor time constant over the full envelope only needs to consider the effect of n1r. Further obtain the fitting formula of N1r and T25 section mach number ma_25, as shown in fig. 2:
Ma_25=0.0033*N1r+0.0267
3. estimating T25 sensor time constant based on airborne measurable parameter N1r
Because the wind tunnel test is to measure the sensor time constant under the measurable air flow rate, according to the analysis, the fitting formula of the Mach number of the T25 section and the T25 time constant is obtained, and then the wind tunnel test result is combined, so that the T25 sensor time constant can be estimated based on the airborne measurable parameter N1r within the full envelope range.
Examples: combining the data of table 2, a fitting formula of t_25 and N1r can be obtained, as shown in fig. 3:
T_25=-0.0224*N1r+6.6483
the application creates a method for evaluating the time constant of the T25 sensor based on the airborne measurable parameters, provides important support for further data analysis, T25 sensor correction and T25 airborne model design, and has very good engineering application value.
In the prior art, a sensor time constant is measured under the condition of stable and measurable air flow rate by adopting a wind tunnel test method, but in the actual use working condition, the sensor time constant changes along with the Mach number of the section of the sensor, and the Mach number of the section T25 is not measurable, so that the time constant measured by the wind tunnel test cannot be directly applied, the data measured by the sensor T25 in the transition process is ambiguous in terms of the difference between the data and the real data, the further analysis and correction are difficult, and the time constant of the sensor T25 reflects the response speed of the sensor, and the acquisition value directly influences the control of a compressor guide vane in the transition process, thereby influencing the stability margin of an engine.
According to the application, based on multi-data correlation analysis, key influence factors of the T25 sensor time constant in the whole envelope range can be obtained, strong correlation factors of the T25 sensor time constant in the conversion rotating speed of the low-pressure compressor are extracted, a fitting formula of the corresponding relation between the T25 sensor time constant and the conversion rotating speed N1r of the low-pressure compressor is established, the problem that the T25 time constant cannot be obtained due to the fact that the Mach number of the T25 section of the engine is not measurable in the whole envelope range is solved, and the problem that the T25 dynamic test data is difficult to analyze due to multi-factor influence is essentially avoided.
The foregoing is merely illustrative of the present application, and the present application is not limited thereto, and any changes or substitutions easily contemplated by those skilled in the art within the scope of the present application should be included in the present application. Therefore, the protection scope of the present application shall be subject to the protection scope of the claims.
Claims (6)
1. A method for evaluating a T25 sensor time constant based on an onboard measurable parameter, comprising:
acquiring a plurality of measurable parameters affecting the Mach number of the T25 section in the full envelope range;
calculating the correlation coefficient of each measurable parameter and the Mach number of the T25 section respectively, and taking the measurable parameter corresponding to the maximum correlation coefficient as the most relevant measurable parameter;
fitting a fitting formula of the most relevant measurable parameters and the Mach number of the T25 section;
and fitting a fitting formula of the most relevant measurable parameter and the time constant of the T25 sensor based on the relation between the Mach number of the T25 section and the time constant of the T25 sensor.
2. The method of evaluating a T25 sensor time constant based on an onboard measurable parameter of claim 1, wherein said measurable parameter comprises: the low-pressure compressor converts the rotating speed N1r, the aircraft altitude H, the engine inlet section air flow speed Mach number Ma and the T25 section air flow speed Mach number.
3. The method for estimating a time constant of a T25 sensor based on onboard measurable parameters of claim 1, wherein the correlation coefficient is calculated by the formula:
x is the Mach number of the T25 section, and Y is a measurable parameter; the correlation coefficient is 0.0 to +/-0.3, the real correlation is +/-0.3 to +/-0.5, and the high correlation is +/-0.8 to +/-1.
4. The method of evaluating a T25 sensor time constant based on an on-board measurable parameter of claim 1, said most relevant measurable parameter preferably being low pressure compressor reduced speed N1r.
5. The method for evaluating the time constant of the T25 sensor based on the airborne measurable parameters of claim 1, wherein the relationship between the Mach number of the T25 section and the time constant of the T25 sensor is obtained through wind tunnel test.
6. The method for evaluating the time constant of the T25 sensor based on the airborne measurable parameters of claim 5, wherein the corresponding relation between the T25 section air flow rate and the time constant of the T25 sensor is obtained through wind tunnel test, and the relation between the Mach number of the T25 section position and the time constant of the T25 sensor is established based on the relation between the T25 section air flow rate and the Mach number of the T25 section position.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311006445.6A CN117030270A (en) | 2023-08-10 | 2023-08-10 | Method for evaluating time constant of T25 sensor based on airborne measurable parameters |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311006445.6A CN117030270A (en) | 2023-08-10 | 2023-08-10 | Method for evaluating time constant of T25 sensor based on airborne measurable parameters |
Publications (1)
Publication Number | Publication Date |
---|---|
CN117030270A true CN117030270A (en) | 2023-11-10 |
Family
ID=88624003
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311006445.6A Pending CN117030270A (en) | 2023-08-10 | 2023-08-10 | Method for evaluating time constant of T25 sensor based on airborne measurable parameters |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117030270A (en) |
-
2023
- 2023-08-10 CN CN202311006445.6A patent/CN117030270A/en active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10788399B2 (en) | Apparatus for evaluating turbine engine system stability | |
CN106777554B (en) | State baseline-based health state evaluation method for air circuit unit body of aero-engine | |
CN110717219B (en) | Method and device for acquiring inlet flow of air compressor in complete machine state of aero-engine | |
US20100024536A1 (en) | Methods and systems for estimating operating parameters of an engine | |
EP3171004A1 (en) | Method and system for improving parameter measurement | |
CN108062428B (en) | Turbofan engine online component fault diagnosis method and system | |
CN109376498B (en) | Turbofan engine modeling method | |
Schönweitz et al. | Experimental and numerical examinations of a transonic compressor-stage with casing treatment | |
CN116127863A (en) | Calculation method for determining influence of Reynolds number on engine performance under complete machine condition | |
CN117030270A (en) | Method for evaluating time constant of T25 sensor based on airborne measurable parameters | |
Kong et al. | Study on condition monitoring of 2-Spool turbofan engine using non-linear gas path analysis method and genetic algorithms | |
Bryce et al. | Three-dimensional flow in a highly loaded single-stage transonic fan | |
CN108760329B (en) | Low-pressure turbine noise test method and improvement method thereof | |
CN114139345B (en) | Turboshaft engine output shaft power prediction method based on test data correlation | |
CN114065661B (en) | Method and system for evaluating loss of rear turbine casing under engine complete machine condition | |
Layachi et al. | Effect of the axial spacing between rotor and stator with regard to the indexing in an axial compressor | |
CN115144186A (en) | Gas turbine engine gas path fault continuous high-precision diagnosis method | |
Cherrett et al. | Unsteady three-dimensional flow in a single-stage transonic fan: Part I—Unsteady rotor exit flow field | |
Castner et al. | An engine research program focused on low pressure turbine aerodynamic performance | |
CN114818205B (en) | Online sensing method for blade tip clearance of full life cycle of aero-engine | |
CN116291879A (en) | Stability expansion control method and system for aero-engine subjected to temperature distortion | |
Yu et al. | EGT baseline model of aeroengine based on kernel principal component analysis and deep belief network | |
Kurzke et al. | Model based performance analysis | |
CN112326226B (en) | Resistance calculation method for impeller type centrifugal ventilator | |
Herrera et al. | A comparative analysis of turbine rotor inlet temperature models |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |