Engine state parameter extraction method corresponding to crystal temperature measurement

Abstract

The invention provides an engine state parameter extraction method corresponding to crystal temperature measurement, and belongs to the field of engine test data processing, and the method comprises the steps: setting a measurement point on an engine provided with a crystal to measure an engine state parameter, operating the engine to a temperature measurement state, and collecting the measurement point data of the engine state parameter; selecting one parameter from the engine state parameters as a discrimination parameter, and screening the collected measurement point data of the engine state parameters and the corresponding transient time according to the discrimination parameter to form a measurement point data set; arranging the measuring point data set according to the collected transient time, calculating the time used by each measuring point, sorting all the measuring points of the discrimination parameters, and calculating the accumulated time of the measuring points; selecting each measurement point of the discrimination parameters with the accumulation time within the target time, and calculating a statistical value of each measurement point; and comparing and de-weighting to obtain a sequencing set, extracting all the other measuring points except the discriminant parameters in the sequencing set, and carrying out statistics to obtain each engine state parameter corresponding to the crystal temperature measurement.

Description

Technical Field

[0001] The present application belongs to the field of engine test data processing, and in particular relates to a method for extracting engine state parameters corresponding to crystal temperature measurement. Background Art

[0002] Crystal temperature measurement technology is a temperature measurement technology based on the temperature memory effect and the thermal stability of irradiation defects in crystal materials. Crystal temperature measurement technology is an ideal measurement method for measuring surface or airflow temperatures at typical locations on aircraft engines.

[0003] When using crystal to measure the engine wall temperature, it is necessary to measure for 3 to 7 minutes at the highest wall temperature test state. However, due to the limitations of test tasks and procedures, the engine state parameters at the highest wall temperature test state, such as gas or air conditioning temperature, pressure, etc., are often non-continuous, intermittent, and show an oscillating trend.

[0004] The existing method for extracting engine state parameters corresponding to the crystal temperature measurement state generally directly selects the engine state parameters within a large state temperature measurement period, and obtains the required data after simple averaging. However, this processing method does not take into account the discontinuity and oscillation of the engine state parameters, and has a large error, which is not conducive to judging the engine state based on the engine state parameters.

[0005] Therefore, a method is needed to extract the engine state parameters corresponding to the crystal temperature measurement. Summary of the invention

[0006] The purpose of the present application is to provide a method for extracting engine state parameters corresponding to crystal temperature measurement to solve or alleviate at least one problem in the background technology.

[0007] The technical solution of the present application is: a method for extracting engine state parameters corresponding to crystal temperature measurement, comprising:

[0008] Setting measuring points on the engine with crystals to measure engine status parameters, running the engine to a temperature measurement state, and collecting measuring point data of the engine status parameters;

[0009] A parameter is selected from the engine state parameters as a discrimination parameter, and according to a predetermined range of the discrimination parameter, the collected measurement point data of the engine state parameter and its corresponding transient time are screened to form a measurement point data set, and the measurement point data set is smoothed;

[0010] Arrange the measurement point data set according to the collected transient time, set the transient time of the first measurement point to zero, and the time used by the remaining measurement points is the difference between the transient time of the next measurement point and the transient time of the previous measurement point, so as to obtain the time used by each measurement point, sort all the measurement points of the discrimination parameters in order from large to small according to the temperature value, and calculate the cumulative time of the measurement points;

[0011] Select each measuring point of the discriminant parameter whose cumulative time is within the target time, extract and rearrange it, record the original order of the measuring point, and calculate the statistics of each selected measuring point;

[0012] The original sorting numbers of the measurement points of the selected discrimination parameters are compared and deduplicated to obtain a sorted set, and all the remaining measurement points except the discrimination parameters in the sorted set are extracted, and statistics are performed on all the remaining measurement points to obtain the engine state parameters corresponding to the crystal temperature measurement.

[0013] Preferably, the engine state parameters include the gas inlet and outlet total temperature, total pressure, static pressure and the cold gas inlet and outlet cavity temperature and cavity pressure.

[0014] Preferably, the discrimination parameter is the pre-turbine temperature or the engine exhaust temperature.

[0015] Preferably, in the screening and collecting data, the predetermined range includes the maximum value of the discrimination parameter, which is used to characterize the highest state of the engine.

[0016] Preferably, the method for smoothing the screened measurement point data set includes the Lowess method.

[0017] Preferably, the target time is 3 minutes to 7 minutes.

[0018] Preferably, the statistical values ​​include an average temperature value, a maximum temperature value, and a minimum temperature value.

[0019] The engine state parameter extraction method corresponding to the crystal temperature measurement of the present application can accurately obtain the corresponding engine state parameters in the crystal measurement process by screening and smoothing the measurement point parameter data of the test state, sorting the discrimination parameters, and calculating the time. The factors that the parameters are non-continuous, intermittent, and present an oscillating trend in the highest state are taken into account, and the calculation results are more accurate. BRIEF DESCRIPTION OF THE DRAWINGS

[0020] In order to more clearly illustrate the technical solution provided by the present application, the following is a brief introduction to the accompanying drawings. Obviously, the accompanying drawings described below are only some embodiments of the present application.

[0021] Figure 1 This is a schematic diagram of the engine state parameter extraction method of the present application. DETAILED DESCRIPTION

[0022] In order to make the purpose, technical solutions and advantages of the implementation of this application clearer, the technical solutions in the embodiments of this application will be described in more detail below in conjunction with the drawings in the embodiments of this application.

[0023] The present application proposes a method for extracting engine status parameters corresponding to the temperature measurement of a certain part of the engine using a crystal. The method accurately extracts the engine status parameters corresponding to the crystal temperature measurement through data screening, smoothing, sorting of discrimination parameters, and time calculation. The method takes into account the fact that the engine status parameters are non-continuous, intermittent, and show an oscillating trend in the highest temperature state, and the calculation result is more accurate.

[0024] like Figure 1 As shown, the method for extracting engine state parameters corresponding to crystal temperature measurement provided by the present application includes the following process:

[0025] Step S10: data measurement.

[0026] A crystal is set on a certain part or component of the engine for temperature measurement. For example, when the maximum wall temperature of the engine needs to be measured, the crystal can be embedded in the outer wall of the turbine cooling blade to measure the maximum wall temperature. At the same time, measurement points for measuring engine state parameters are arranged as needed, and the engine is operated to collect engine state parameters during engine operation. These engine state parameters include the total temperature, total pressure, static pressure of the gas inlet and outlet, and the cavity temperature and cavity pressure of the cold air inlet and outlet. The above measurement point data can be measured by setting temperature measurement points and pressure measurement points at the turbine inlet and outlet and the cold air flow path inlet and outlet.

[0027] Step S20: data preprocessing.

[0028] S21, data screening:

[0029] A parameter is selected from the engine parameters as a discrimination parameter, and according to a predetermined range of the discrimination parameter, the collected engine state parameter data and its corresponding transient time measurement point data are screened to form a measurement point data set. In a preferred embodiment of the present application, the turbine pre-temperature or the engine exhaust temperature (i.e., the total temperature of the gas outlet) is selected as a discrimination parameter for engine state determination and monitoring during the crystal temperature measurement process.

[0030] The screening includes the required engine state parameters for measuring the large state area (i.e., the engine high speed state) and the transient time corresponding to the engine state parameter measurement. The large state area can be judged by discriminating parameters such as the turbine pre-temperature or the engine exhaust temperature. When screening, a discriminating parameter is selected and a certain range is given. When the discriminating parameter meets the range requirements, all types of parameter data of the corresponding time are extracted.

[0031] For example, in the embodiment of the present application, the pre-turbine temperature is selected as the discrimination parameter, and considering the uneven distribution of the combustion chamber outlet and the measurement deviation, (T4*-1000, T4*max) is selected as the screening range, wherein T4* is the pre-turbine temperature of the corresponding state of the crystal temperature measurement, and T4*max is the maximum value of all T4* measurement points.

[0032] S22, data smoothing: smoothing the filtered data.

[0033] Affected by engine operation and test accuracy, normal measurement point data will also show up and down jumps within a small range, so the filtered measurement point data set is smoothed.

[0034] For example, in this embodiment of the present application, Lowess (local weighted regression algorithm) is used to smooth the temperature before the turbine, and local regression is performed using weighted linear least squares based on a first-order polynomial model.

[0035] Step S30, calculating the time used for measuring the points.

[0036] Time calculation: Arrange the measurement point data set according to the transient time collected during the test, and give a sorting number n. Set the time used by the first measurement point to 0, and the time used by the remaining measurement points is the transient time of the next measurement point minus the transient time of the previous measurement point, so as to calculate the time used by each measurement point. Since the measurement point data collected during the test is very dense, the time is magnified.

[0037] Table 1 shows an example of time magnification of the measurement point data in this embodiment of the present application. In this embodiment, the time is magnified by 10 7 times.

[0038] Table 1 Time calculation (partial data)

[0039] Sort number Raw instant time - double precision format Zoomed in time 1 0.445548379629630 0 2 0.445549074074074 6.94444444504505 3 0.445549884259259 8.10185185162737 4 0.445550578703704 6.94444444448994 5 0.445552546296296 19.6759259257773 6 0.445553356481482 8.10185185162737 7 0.445553935185185 5.78703703679739 8 0.445555324074074 13.8888888895350 9 0.445556365740741 10.4166666659022 10 0.445557060185185 6.94444444504505 11 0.445558564814815 15.0462962961173 12 0.445559490740741 9.25925925931992 13 0.445560069444444 5.78703703679739 14 0.445561805555556 17.3611111109473 15 0.445562384259259 5.78703703735251 16 0.445563773148148 13.8888888889799 17 0.445565277777778 15.0462962966724 18 0.445565972222222 6.94444444393483 19 0.445566782407407 8.10185185218248 20 0.445567592592593 8.10185185162737 …… …… …… 340 0.445813657407407 6.94444444448994 341 0.445814351851852 6.94444444448994 342 0.445815162037037 8.10185185162737 343 0.445815856481482 6.94444444504505 344 0.445816666666667 8.10185185162737 …… …… ……

[0040] Afterwards, all the measurement points of the discrimination parameters are sorted in descending order according to the temperature value, and the cumulative time of the measurement points is calculated.

[0041] Table 2 shows an example of sorting some measurement points of the pre-turbine temperature T4* in this embodiment of the present application.

[0042] Table 2 Sorted T4* measurement point 1 (partial data)

[0043]

[0044]

[0045]

[0046] Step S40, extraction and calculation of discrimination parameters.

[0047] Select each measuring point of the discrimination parameter whose cumulative time is within the target time (including the target time), and record the original sorting number n_i of each. In this application, the target time is selected within the range of 3-7 minutes.

[0048] For example, in one embodiment of the present application, the target time is 5 minutes, and the judgment parameters of measuring point 1, measuring point 4, and measuring point 16 extract 4238 groups, 3554 groups, and 4113 groups of data respectively, as shown in Table 3.

[0049] Table 3 Original order n_i of T4* measuring points within target time (partial data)

[0050]

[0051]

[0052] Perform discriminant parameter analysis and calculate the average temperature value Tta_i, maximum temperature value Ttmax_i, minimum temperature value Ttmin_i, distribution law, etc. of each selected measuring point.

[0053] As shown in Table 4, there are the average temperature value Tta_i, the maximum temperature value Ttmax_i and the minimum temperature value Ttmin_i of the discrimination parameters of this embodiment of the present application.

[0054] Table 4 Average temperature, maximum value, minimum value (partial parameters) of T4* measuring point

[0055] Average temperature(K) Maximum temperature(K) Minimum temperature (K) T4* measuring point 1 1329.413 1350.338 1310.154 T4*Measurement point 4 1439.361 1462.034 1418.196 T4* measuring point 16 1547.851 1557.766 1539.195

[0056] Step S50: extracting engine status parameters.

[0057] The original sorting numbers n_i of the selected discriminant parameters are compared and removed to obtain a sorting set N. All other measuring points except the discriminant parameters are extracted from the sorting set N. Statistical analysis is performed on the remaining measuring points, such as obtaining the average value, maximum value, minimum value, distribution law, etc., to obtain the engine state parameters corresponding to the crystal temperature measurement.

[0058] In the embodiment shown in Table 5, the original sorting n_i of the 24 turbine pre-temperature T4* measuring points is compared, and duplicate data is removed to obtain the sorting set N. The remaining measuring point data after extraction are shown in Table 6. As shown in Table 7, according to the above method, the engine state parameters such as the total pressure and static pressure of the gas inlet and outlet, the cavity temperature and cavity pressure of the cold gas inlet and outlet are extracted, and the average value, maximum value, minimum value, distribution law, etc. of each engine state parameter are calculated. At this point, the engine state parameters corresponding to the crystal temperature measurement are obtained.

[0059] Table 5 Repeated T4* measurement point numbers (partial data)

[0060] Duplicate measurement point numbers Repetitions 1 2537 11 2 2538 11 3 2539 11 4 2540 14 5 2541 14 6 2542 21 7 2543 18 8 2544 19 9 2545 18 10 2546 21 11 2547 21 12 2548 21 13 2549 16 14 2550 16 15 2551 15 16 2552 19 17 2553 20 18 2554 16 19 2555 15 20 2556 15 …… …… 340 2876 13 341 2877 15 342 2878 14 343 2879 12 344 2880 11 …… …… 3551 11844 10 3552 11845 11 3553 11846 18 3554 11847 12 3555 11848 12 …… …… ……

[0061] Table 6 Data of other measuring points (partial parameters)

[0062]

[0063]

[0064] Table 7 Average temperature, maximum value, minimum value (partial parameters) of other measurement points

[0065] average value Maximum Minimum Gas inlet total pressure measuring point 2-3 (KPa) 2688.69 2721.96 2622.2 Gas inlet static pressure measuring point 2 (KPa) 1269.32 1298.42 1233.36 Air conditioning inlet temperature measurement point 2 (K) 569.85 580.00 549.84 Air inlet pressure measurement point 1 (KPa) 2762.36 2795.19 2694.73 …… …… …… ……

[0066] The engine state parameter extraction method corresponding to the crystal temperature measurement of the present application can accurately obtain the corresponding engine state parameters in the crystal measurement process by screening and smoothing the measurement point parameter data of the test state, sorting the discrimination parameters, and calculating the time. The factor that the parameters are non-continuous, intermittent, and show an oscillating trend in the highest state is taken into account, and the calculation result is more accurate.

[0067] The above is only a specific implementation of the present application, but the protection scope of the present application is not limited thereto. Any changes or substitutions that can be easily thought of by a person skilled in the art within the technical scope disclosed in the present application should be included in the protection scope of the present application. Therefore, the protection scope of the present application shall be based on the protection scope of the claims.

Claims

1. A method for extracting engine state parameters corresponding to crystal temperature measurement, characterized in that: include: Setting measuring points on the engine with crystals to measure engine status parameters, running the engine to a temperature measurement state, and collecting measuring point data of the engine status parameters; A parameter is selected from the engine state parameters as a discrimination parameter, and according to a predetermined range of the discrimination parameter, the collected measurement point data of the engine state parameter and its corresponding transient time are screened to form a measurement point data set, and the measurement point data set is smoothed; Arrange the measurement point data set according to the collected transient time, set the time of the first measurement point to zero, and the time of the remaining measurement points to the difference between the transient time of the next measurement point and the transient time of the previous measurement point, so as to obtain the time used by each measurement point, sort all the measurement points of the discrimination parameters in descending order according to the temperature value, and calculate the cumulative time of the measurement points; Select each measuring point of the discriminant parameter whose cumulative time is within the target time, record the original sorting sequence number of the measuring point, and calculate the statistical value of each selected measuring point; The original sorting numbers of the measurement points of the discrimination parameters are compared and deduplicated to obtain a sorting set, and all the remaining measurement points except the discrimination parameters in the sorting set are extracted. All the remaining measurement points are counted to obtain the engine state parameters corresponding to the crystal temperature measurement.

2. The method for extracting engine state parameters corresponding to crystal temperature measurement according to claim 1, characterized in that: The engine status parameters include the total temperature, total pressure, static pressure of the gas inlet and outlet, and the cavity temperature and cavity pressure of the cold gas inlet and outlet.

3. The method for extracting engine state parameters corresponding to crystal temperature measurement according to claim 2, characterized in that: The discrimination parameter is the pre-turbine temperature or the engine exhaust temperature.

4. The method for extracting engine state parameters corresponding to crystal temperature measurement according to claim 3, characterized in that: In the screening and collecting data, the predetermined range includes the maximum value of the discrimination parameter, which is used to characterize the highest state of the engine.

5. The method for extracting engine state parameters corresponding to crystal temperature measurement according to claim 1, characterized in that: Methods for smoothing the selected measurement point data set include the Lowess method.

6. The method for extracting engine state parameters corresponding to crystal temperature measurement according to claim 1, characterized in that: The target time is 3 minutes to 7 minutes.

7. The method for extracting engine state parameters corresponding to crystal temperature measurement according to claim 1, characterized in that: The statistical values ​​include an average value, a maximum value, and a minimum value.

Images