Method and system for processing test data of auxiliary power device of gas turbine

Abstract

The invention belongs to the technical field of device test, and discloses a gas turbine auxiliary power device test data processing method and system, and the method comprises the steps: reading the test data of a gas turbine auxiliary power device, recognizing the test data, calculating the signal characteristics of the test data, and obtaining the test data of the gas turbine auxiliary power device according to the signal characteristics of the test data. Dividing a steady-state interval and a transition-state interval of the test data, representing a switching point of the steady-state interval and the transition-state interval of the test data by using step characteristics of a plurality of signals, calculating interval signal characteristics of the test data according to the test data of each interval, and determining the interval signal characteristics of the test data according to the interval signal characteristics of the test data. And calculating data processing parameters of each interval, smoothing the data curve of the test data of the transition state interval by using a filter, and averaging the test data of the steady state interval to obtain a unique steady state value. According to the method, test original data features are calculated, steady-state and transition-state data processing intervals are divided, data processing parameters of each interval are calculated, and data are processed. Compared with the prior art, the method has the advantages that the test data of the auxiliary power device are efficiently processed in a self-adaptive signal feature extraction mode, the obtained processing result has uniqueness, subjective deviation is avoided, meanwhile, the labor cost is saved, and the test analysis period is shortened.

Description

Technical Field

[0001] The invention belongs to the technical field of device testing, and relates to a test method for a gas turbine auxiliary power device, and in particular to a test data processing method and system for a gas turbine auxiliary power device. Background Art

[0002] Gas turbine auxiliary power units work under harsh conditions of high temperature, high pressure and high speed. There are complex aerodynamic and thermal processes inside them. Testing is the only means to verify and modify the design, which runs through the entire life cycle of auxiliary power unit development, testing, production, use and maintenance. Test data processing is an important part of guiding and verifying the design of auxiliary power units, and is an important criterion for evaluating the performance of auxiliary power unit components and the whole machine.

[0003] The data collected from the auxiliary power unit test has the characteristics of multiple signal types, large quantities, and multiple steady-state transition state intervals. Signal types include temperature, pressure, flow, voltage, current, speed, vibration, etc. The measurement channels range from dozens to thousands, and different types of data often show different signal characteristics. There are multiple different steady-state processes in the same test, and the transition state process during the steady-state switching process is changeable. However, the current bench test can only display the collected raw data and lacks the preliminary data processing process. At present, testers and engineers can only read the collected data at a certain moment, and it is difficult for different personnel to obtain consistent results, which is subjective. When engineers analyze the aerodynamic performance or mechanical performance of the product, they need to spend a lot of energy to process the test data. Summary of the invention

[0004] Purpose of the present invention: In order to solve the above problems, the present invention provides a method and system for processing test data of a gas turbine auxiliary power unit, so as to avoid the subjective problem of test data reading, reduce the data processing workload, and shorten the data processing cycle.

[0005] The technical solution of the present invention is as follows

[0006] A method for processing test data of a gas turbine auxiliary power unit comprises the following steps:

[0007] S1, read the test data of the gas turbine auxiliary power unit;

[0008] S2, identifying the test data and calculating the signal characteristics of the test data;

[0009] S3, dividing the steady-state interval and the transition-state interval of the test data according to the signal characteristics of the test data, and characterizing the switching points of the steady-state interval and the transition-state interval of the test data by using step characteristics of multiple signals;

[0010] S4, calculating the interval signal characteristics of the test data according to the test data of each interval;

[0011] S5, according to the interval signal characteristics of the test data, the data processing parameters of each interval are calculated, the test data of the transition state interval are smoothed by a filter, and the test data of the steady state interval are averaged to obtain a unique steady state value.

[0012] Furthermore, in S2, the signal characteristics of the test data include step signal, variance, range, interquartile range, coefficient of variation, skewness and kurtosis; step signal and coefficient of variation are used to characterize the characteristics of different signals of the test data in a specific time period; variance, range, interquartile range, skewness and kurtosis are used to characterize the degree of dispersion of the test data.

[0013] Furthermore, the type and method of the interval signal feature in S4 are the same as those of the signal feature in S2.

[0014] Furthermore, in S3, the steady-state interval and the transition state interval of the test data appear alternately, the steady-state interval is the stable working state of the gas turbine auxiliary power unit, and the transition state interval is the dynamic adjustment state between the two steady-state intervals.

[0015] Furthermore, in S3, the method of characterizing the switching point between the steady-state interval and the transition state interval by the step characteristics of multiple signals is:

[0016] After obtaining the step characteristics of all signals, they are arranged in chronological order, and the time when three or more signal step characteristics are close is taken as the interval division time. The transition state is characterized by the obvious change trend of the window mean. The window width is set to 10 consecutive sampling data. The obvious change trend is characterized by the following formula:

[0017] K = sum(|Mean win,i -Mean|>Mean·p%), i=1,2,3...,

[0018] Among them, K represents the number of signal windows with obvious change trend of a certain parameter, Mean win,i represents the mean of the i-th window, Mean represents the interval mean of the parameter, and p represents the percentage of the mean deviation.

[0019] Furthermore, K>5 indicates that the parameter has an obvious changing trend in the interval; if more than 10 parameters have an obvious changing trend in the same interval, the interval is considered to be a transition state interval; otherwise, it is a steady state interval.

[0020] Furthermore, in S5, the experimental data in the transition state interval is smoothed by a low-pass filter, and the experimental data in the steady state interval is represented by a weighted average.

[0021] A gas turbine auxiliary power unit test data processing system, using the aforementioned gas turbine auxiliary power unit test data processing method, comprises an interaction module, a signal feature extraction module, an interval division module and a data processing module;

[0022] The interactive module reads the test data file according to the data file directory input by the user, draws the characteristic data chart according to the user input signal and the data processing result, and saves it to the specified location. The interactive module also saves the process data file and the result data file according to the data of the interval division module, the feature extraction module and the data processing module; draws the color data characteristic image for the result data file, and saves it to the specified location;

[0023] The signal feature extraction module calculates the signal features based on the input test data. When the input data is all the test data, the calculated signal features are the overall signal features. When the input data is the processing interval data, the calculated signal features are the interval signal features.

[0024] The interval division module divides the steady-state and transition-state intervals according to the signal characteristics;

[0025] The data processing module calculates the data processing parameters of each interval according to the signal characteristics and obtains the data processing results.

[0026] Beneficial effects of the present invention:

[0027] The auxiliary power unit test data processing method provided by the present invention calculates the test original data characteristics, divides the steady-state and transitional state data processing intervals, calculates the data processing parameters of each interval and processes the data. Compared with the prior art, the method realizes efficient processing of the auxiliary power unit test data by adaptively extracting signal characteristics, and the obtained processing results are unique, avoiding subjective deviations, saving labor costs, and shortening the test analysis cycle. BRIEF DESCRIPTION OF THE DRAWINGS

[0028] In order to more clearly illustrate the technical solution of the implementation method of the patent of the present invention, the following is a brief introduction to the drawings required for use in the implementation method. It should be understood that the following drawings only show certain embodiments of the patent of the present invention and should not be regarded as limiting the scope. For ordinary technicians in this field, other relevant drawings can be obtained based on these drawings without paying creative work.

[0029] Figure 1 It is a flow chart of a method for processing test data of a gas turbine auxiliary power unit disclosed in Example 1 of the present application.

[0030] Figure 2 It is a structural schematic block diagram of a gas turbine auxiliary power unit test data processing system disclosed in Example 2 of the present application.

[0031] Figure 3 It is an interactive interface of a gas turbine auxiliary power unit test data processing system disclosed in Example 2 of the present application.

[0032] Figure 4 This is a processing result diagram of a gas turbine auxiliary power unit test data processing system disclosed in Example 2 of the present application. DETAILED DESCRIPTION

[0033] This section is an embodiment of the present invention, which is used to explain and illustrate the technical solution of the present invention. In the absence of conflict, the embodiments of the present invention and the features in the embodiments can be combined with each other.

[0034] In the description of the present invention, it should be understood that the terms "center", "longitudinal", "lateral", "up", "down", "front", "back", "left", "right", "vertical", "horizontal", "top", "bottom", "inside", "outside" and the like indicate directions or positional relationships based on the orientations or positional relationships in the accompanying drawings, and are only for the convenience of describing the present invention and simplifying the description, rather than indicating or implying that the device or case referred to must have a specific orientation, be constructed and operated in a specific orientation, and therefore cannot be understood as limiting the present invention. In addition, the terms "first", "second", etc. are used for descriptive purposes only, and cannot be understood as indicating or implying relative importance or implying the number of technical features indicated. Thus, features defined as "first", "second", etc. may explicitly or implicitly include more of the features. In the description of the present invention, unless otherwise specified, "multiple" means two or more.

[0035] In the description of the present invention, it should be noted that, unless otherwise clearly specified and limited, the terms "installed", "connected" and "connected" should be interpreted in a broad sense, for example, it can be a fixed connection, a detachable connection or an integrated connection; it can be a mechanical connection or a point connection; it can be a direct connection or an indirect connection through an intermediate medium, or it can be the internal communication of two components. For ordinary technicians in this field, the specific meanings of the above terms in the present invention can be understood by specific circumstances.

[0036] Embodiment 1:

[0037] A method for processing test data of a gas turbine auxiliary power unit comprises the following steps:

[0038] S1, read the test data of the gas turbine auxiliary power unit;

[0039] S2, identifying the test data and calculating the signal characteristics of the test data;

[0040] S3, dividing the steady-state interval and the transition-state interval of the test data according to the signal characteristics of the test data, and characterizing the switching points of the steady-state interval and the transition-state interval of the test data by using step characteristics of multiple signals;

[0041] S4, calculating the interval signal characteristics of the test data according to the test data of each interval;

[0042] S5, according to the interval signal characteristics of the test data, the data processing parameters of each interval are calculated, the test data of the transition state interval are smoothed by a filter, and the test data of the steady state interval are averaged to obtain a unique steady state value.

[0043] In S2, the signal characteristics of the test data include step signal, variance, range, interquartile range, coefficient of variation, skewness and kurtosis; step signal and coefficient of variation are used to characterize the characteristics of different signals of the test data in a specific time period; variance, range, interquartile range, skewness and kurtosis are used to characterize the degree of dispersion of the test data.

[0044] The type and method of the interval signal features in S4 are the same as those of the signal features in S2.

[0045] In S3, the steady-state interval and the transition state interval of the test data appear alternately. The steady-state interval is the stable working state of the gas turbine auxiliary power unit, and the transition state interval is the dynamic adjustment state between the two steady-state intervals.

[0046] In S3, the method of characterizing the switching points between the steady-state interval and the transition state interval by the step characteristics of multiple signals is:

[0047] After obtaining the step characteristics of all signals, they are arranged in chronological order, and the time when three or more signal step characteristics are close is taken as the interval division time. The transition state is characterized by the obvious change trend of the window mean. The window width is set to 10 consecutive sampling data. The obvious change trend is characterized by the following formula:

[0048] K = sum(|Mean win,i -Mean|>Mean·p%), i=1,2,3...,

[0049] Among them, K represents the number of signal windows with obvious change trend of a certain parameter, Mean win,i represents the mean of the i-th window, Mean represents the interval mean of the parameter, and p represents the percentage of the mean deviation.

[0050] K>5 means that the parameter has an obvious changing trend in the interval; if more than 10 parameters have an obvious changing trend in the same interval, the interval is considered to be a transition state interval; otherwise, it is a steady state interval.

[0051] In S5, the experimental data in the transition state interval are smoothed by a low-pass filter, and the experimental data in the steady state interval are represented by a weighted average.

[0052] A gas turbine auxiliary power unit test data processing system, using the aforementioned gas turbine auxiliary power unit test data processing method, comprises an interaction module, a signal feature extraction module, an interval division module and a data processing module;

[0053] The interactive module reads the test data file according to the data file directory input by the user, draws the characteristic data chart according to the user input signal and the data processing result, and saves it to the specified location. The interactive module also saves the process data file and the result data file according to the data of the interval division module, the feature extraction module and the data processing module; draws the color data characteristic image for the result data file, and saves it to the specified location;

[0054] The signal feature extraction module calculates the signal features based on the input test data. When the input data is all the test data, the calculated signal features are the overall signal features. When the input data is the processing interval data, the calculated signal features are the interval signal features.

[0055] The interval division module divides the steady-state and transition-state intervals according to the signal characteristics;

[0056] The data processing module calculates the data processing parameters of each interval according to the signal characteristics and obtains the data processing results.

[0057] Embodiment 2:

[0058] like Figure 1 As shown, the present invention provides a flow chart of a method for processing test data of a gas turbine auxiliary power unit, the method comprising:

[0059] Step 1, read the test data;

[0060] Step 2, calculating the overall signal characteristics based on the test data;

[0061] In this embodiment, the overall signal characteristics include step signal, variance, range, interquartile range, coefficient of variation, skewness and kurtosis.

[0062] Among them, the step signal and the coefficient of variation are used to characterize the characteristics of different signals in a specific time period. For example, the coefficient of variation of the bleed air mass flow rate and the coefficient of variation of the bleed air temperature can be used to characterize different steady-state processes, and the starting fuel flow step signal and the anti-surge valve opening signal are used to mark the determined steady-state and transition state change moments.

[0063] Variance, range, interquartile range, skewness and kurtosis are used to characterize the degree of data dispersion. For example, variance can be used to characterize the degree of dispersion of atmospheric temperature. The atmospheric temperature changes slowly, and the sensor measurement deviation is small. The measurement deviation of the bleed air mass flow sensor is large, and its dispersion degree needs to be characterized by five signal characteristic parameters: variance, range, interquartile range, skewness and kurtosis.

[0064] Step 3: Divide the steady state and transition state intervals according to signal characteristics.

[0065] In this embodiment, the steady state and transition state intervals of the test data appear alternately, the steady state refers to the stable working state of the gas turbine auxiliary power device, and the transition state refers to the dynamic adjustment process between the two steady states.

[0066] Optionally, the steady-state and transition state switching points are characterized by the step characteristics of multiple signals. After obtaining the step characteristics of all signals, they are arranged in chronological order, and the time when the step characteristics of 3 or more signals are close is taken as the interval division time. In particular, 3 consecutive sampling data moments are considered to be close. The transition state is characterized by the presence of an obvious change trend in the window mean, and the window width is set to 10 consecutive sampling data. The obvious change trend is characterized by the following formula:

[0067] K = sum(|Mean win,i -Mean|>Mean·p%), i=1,2,3...,

[0068] Among them, K represents the number of signal windows with obvious change trend of a certain parameter, Mean win,i represents the mean of the i-th window, Mean represents the interval mean of the parameter, and p represents the percentage of the mean deviation.

[0069] Optionally, it is considered that K>5 indicates that the parameter has an obvious changing trend in the interval; if more than 10 parameters have an obvious changing trend in the same interval, the interval is considered to be a transition state interval; otherwise, it is a steady state interval.

[0070] Step 4: Calculate the interval signal feature based on each interval data. The feature type is the same as step 2.

[0071] Step 5, calculate the data processing parameters of each interval according to the interval signal characteristics, use a filter to smooth the data curve in the transition state, and use the average value in the steady state to obtain a unique steady-state value.

[0072] This embodiment takes the rotation speed as an example to provide a general method for setting filter parameters used for different signals.

[0073] Optionally, the filter parameters include window width, cutoff frequency, weight parameters, etc. depending on the filter type.

[0074] In this embodiment, all transition states use low-pass filters to smooth data curves. For the rotation speed, a low-pass Chebyshev filter is used. The filter parameters include the cutoff frequency and the filter order. The cutoff frequency calculation formula can be:

[0075]

[0076] Where n is the maximum speed measured by the auxiliary power unit, k is the multiple of the cutoff frequency and the parameter natural frequency;

[0077] The filter order calculation formula is:

[0078]

[0079] Among them, Mean win,i represents the mean of the i-th window, Range win,i It represents the range of the i-th window, round is the rounding function, max is the maximum value function, and the window width is set to 10 consecutive sampling data.

[0080] In this embodiment, all steady states are represented by weighted mean. The speed steady state calculation formula is:

[0081]

[0082] Among them, a, b, c are weight constants, is the average speed, n q1 , n q2 , n q3 are the first quartile, the second quartile (median), and the third quartile of the interquartile range, and S is the square root of the variance (standard deviation).

[0083] As can be seen from the above embodiments of the present invention, the embodiments of the present invention first calculate different signal features based on the original test data, divide multiple processing intervals within the full test cycle, calculate signal features for each data interval, set filter parameters in the transition state to process the transition state interval data, and use weighted average values ​​in the steady state to obtain a unique steady state value. Compared with the prior art, this embodiment realizes efficient processing of auxiliary power unit test data by adaptively extracting signal features, and the obtained processing results are unique, avoiding subjective bias, while saving labor costs and shortening the test analysis cycle.

[0084] Embodiment 3:

[0085] like Figure 2 The present invention discloses a gas turbine auxiliary power unit test data processing system, comprising:

[0086] The interactive module is used to read the test data file according to the data file directory input by the user, and to draw the characteristic data chart according to the user input signal and data processing results and save it to the specified location.

[0087] In this embodiment, the test data file is acquired by the software and hardware of the test bench, and contains test data and other test information, wherein the test data is the operation object of this embodiment. The interactive module is used to read the test data file and extract the test data into a table for use by other modules.

[0088] like Figure 3 As shown, Figure 3 This is an interactive module operation interface disclosed in Example 2 of the present application.

[0089] Furthermore, a test data file represents the test data recorded from the last start to the first stop of the auxiliary power unit; other test information includes the test name, test time, operator, test description, etc.

[0090] In this embodiment, the interactive module can also be used to save process data files and result data files according to the data of the interval division module, feature extraction module and data processing module; and to draw a color data feature image for the result data file and save it to a specified location.

[0091] Furthermore, the data feature images include the test data time variation curve, the original test data box plot, the time proportion graph of different steady-state intervals and the bubble graph. Figure 4 As shown, Figure 4 This is the time variation curve of the processed test data disclosed in Example 2 of this application.

[0092] A signal feature extraction module is used to calculate signal features based on test data;

[0093] In this embodiment, the signal feature extraction module is used to calculate the signal features according to the input test data, including step signal, variance, range, interquartile range, coefficient of variation, skewness and kurtosis. When the input data is all test data, the overall signal features are calculated, and when the input data is processing interval data, the interval signal features are calculated.

[0094] An interval division module is used to divide the steady-state and transition-state intervals according to signal characteristics;

[0095] In this embodiment, the steady state and transition state intervals of the test data appear alternately, the steady state refers to the stable working state of the gas turbine auxiliary power device, and the transition state refers to the dynamic adjustment process between the two steady states.

[0096] Optionally, the steady-state and transition state switching points are characterized by the step characteristics of multiple signals. After obtaining the step characteristics of all signals, they are arranged in chronological order, and the time when the step characteristics of 3 or more signals are close is taken as the interval division time. In particular, 3 consecutive sampling data moments are considered to be close. The transition state is characterized by the presence of an obvious change trend in the window mean, and the window width is set to 10 consecutive sampling data. The obvious change trend is characterized by the following formula:

[0097] K = sum(|Mean win,i -Mean|>Mean·p%), i=1,2,3...,

[0098] Among them, K represents the number of signal windows with obvious change trend of a certain parameter, Mean win,i represents the mean of the i-th window, Mean represents the interval mean of the parameter, and p represents the percentage of the mean deviation.

[0099] Optionally, it is considered that K>5 indicates that the parameter has an obvious changing trend in the interval; if more than 10 parameters have an obvious changing trend in the same interval, the interval is considered to be a transition state interval; otherwise, it is a steady state interval.

[0100] A data processing module is used to calculate the data processing parameters of each interval according to the signal characteristics to obtain the data processing results;

[0101] This embodiment takes the rotation speed as an example to provide a general method for setting filter parameters used for different signals.

[0102] Optionally, the filter parameters include window width, cutoff frequency, weight parameters, etc. depending on the filter type.

[0103] In this embodiment, all transition states use low-pass filters to smooth data curves. For the rotation speed, a low-pass Chebyshev filter is used. The filter parameters include the cutoff frequency and the filter order. The cutoff frequency calculation formula can be:

[0104]

[0105] Where n is the maximum speed measured by the auxiliary power unit, k is the multiple of the cutoff frequency and the parameter natural frequency;

[0106] The filter order calculation formula is:

[0107]

[0108] Among them, Mean win,i represents the mean of the i-th window, Range win,iIt represents the range of the i-th window, round is the rounding function, max is the maximum value function, and the window width is set to 10 consecutive sampling data.

[0109] In this embodiment, all steady states are represented by weighted mean. The speed steady state calculation formula is:

[0110]

[0111] Among them, a, b, c are weight constants, is the average speed, n q1 , n q2 , n q3 are the first quartile, the second quartile (median), and the third quartile of the interquartile range, and S is the square root of the variance (standard deviation).

[0112] Through a gas turbine auxiliary power unit test data processing system of this embodiment, the corresponding aviation integrated electromechanical gas turbine auxiliary power unit test data processing method in the aforementioned method embodiment can be implemented, and has the beneficial effects of the corresponding method embodiment. The specific beneficial effects will not be repeated here.

[0113] Embodiment 4:

[0114] Step 1, read the original data file, remove the test name, test time, operator, test description and other test related information, and save the original test data table;

[0115] Step 2, calculate the overall mean of the measured parameters one by one, and then use the sliding window method to calculate the window mean. If a certain number of window means exceed or fall below the overall mean, it means that the parameter has changed significantly and is marked as a valid parameter. Difference is calculated for the time series of the test records to obtain the auxiliary power unit working time series, and the auxiliary power unit characteristic data table is obtained.

[0116] Step 3: Get the step moment of the position feedback parameter and mark it as the starting point of the auxiliary power unit working range; get the speed stabilization moment and mark it as the end point of the auxiliary power unit working range. The speed stability in the range is marked as the steady-state range, and the speed fluctuation is marked as the dynamic range. The test data is divided into multiple different steady-state mode and dynamic mode working ranges.

[0117] Step 4: Extract the data of each working interval into a separate characteristic data table and mark them in order of time to facilitate drawing a complete test curve after data statistics are completed.

[0118] Step 5: The dynamic interval data is filtered using a sliding mean filter, and the steady-state interval data is segmentedly calculated using the sliding mean method. For parameters with constant steady-state performance, the only data output that can characterize the steady-state performance is selected, and for parameters with continuously changing steady-state performance, the sliding mean data is retained.

[0119] Step 6, rearrange the interval data into a complete data table in chronological order, and calculate the auxiliary power unit performance related parameters, including converted flow, total output shaft power, total fuel consumption, fuel consumption rate, thermal efficiency, intake flow rate, exhaust flow rate, average exhaust temperature, average bleed air temperature, average bleed air pressure, bleed air equivalent power, dynamic time, steady-state time, speed acceleration, speed overshoot; use the auxiliary power unit simulation model to calculate the performance value under steady-state conditions, and calculate the relative deviation between the test data and the model simulation calculation results.

[0120] Step 7, draw the auxiliary power unit characteristic data graph and save it in the form of a picture to a designated location.

[0121] It should be noted that, when collecting test data once, a corresponding characteristic data table of the auxiliary power device is obtained, and data processing is performed according to the obtained original test data table.

[0122] This embodiment takes the original data obtained from the performance test of the auxiliary power unit ground sample as an example.

[0123] Read the original data file, remove the test name, test time, operator, test description and other test related information, and save the original test data table;

[0124] The overall mean of the measured parameters is calculated one by one, and then the sliding window method is used to calculate the window mean. In this embodiment, due to the large amount of data, the window width is selected as 20. If there are a certain number of window means that exceed or fall below the overall mean, it means that the parameters have changed significantly and are marked as valid parameters. The valid parameters of this embodiment include bleed air shut-off valve pressure, bleed air mass flow, bleed air valve position feedback, hydraulic pump oil supply flow, hydraulic valve return oil flow, bleed air total pressure, bleed air static pressure, bleed air temperature, hydraulic pump outlet pressure, fuel flow, compressor inlet pressure, inlet atmospheric temperature, exhaust total pressure, exhaust static pressure, exhaust temperature, speed monitoring, shut-off valve inlet pressure, pressure after pressure reducing valve, generator excitation voltage, permanent magnet machine voltage, lubricating oil radiator outlet oil pressure, lubricating oil radiator inlet oil pressure, lubricating oil flow, lubricating oil radiator outlet water temperature, hydraulic pump input shaft power, fuel mass flow, exhaust temperature, bleed air coarse adjustment valve adjustment, bleed air fine adjustment valve adjustment, hydraulic pump I oil temperature, hydraulic pump I oil supply flow, hydraulic pump I return flow, hydraulic pump II oil temperature, hydraulic pump II oil supply flow, hydraulic pump II return flow, lubricating oil supply temperature, total lubricating oil return flow, generator load voltage, generator load current, fuel centrifugal pump regulating valve position feedback, fuel centrifugal pump oil supply flow, fuel supply pressure, bleed air shut-off valve temperature, plant inlet pressure, plant inlet temperature, plant inlet regulating valve position feedback, hydraulic pump I unloading valve current, hydraulic pump II unloading valve current, generator lubricating oil flow, generator lubricating oil inlet temperature, generator lubricating oil outlet temperature, motor end temperature, motor slot temperature, starting motor torque, starting motor speed, starting power supply current, etc. The auxiliary power unit working time sequence is obtained by subtracting the time series of the test records, the sequence start time is 0 hours, 0 minutes, 0 seconds, and the year, month, and day information is ignored to obtain the auxiliary power unit characteristic data table;

[0125] The position feedback parameters include the position feedback of the fuel centrifugal pump regulating valve and the position feedback of the air induction valve. The step moment of the position feedback parameters is obtained and marked together with the step moment of the hydraulic pump input shaft power as the starting point of the auxiliary power unit working range; the speed stabilization moment is obtained and the end point of the auxiliary power unit working range is marked. The speed stability within the interval is marked as the steady-state interval, and the speed fluctuation is marked as the dynamic interval. The test data is divided into multiple different steady-state mode and dynamic mode working intervals. In this embodiment, 14 interval moments are obtained, and a total of 3 steady-state working intervals are divided.

[0126] The data of each working interval is extracted into a separate characteristic data table and marked in time sequence to facilitate drawing a complete test curve after data statistics are completed.

[0127] The dynamic interval data is filtered by sliding mean, and the steady-state interval data is calculated segmentally by sliding mean method. For parameters with constant steady-state performance, the only data output that can characterize the steady-state performance is selected, and the sliding mean data is retained for parameters with continuously changing steady-state performance.

[0128] The interval data are rearranged into a complete data table in chronological order, and the auxiliary power unit performance-related parameters are calculated, including converted flow, total output shaft power, total fuel consumption, fuel consumption rate, thermal efficiency, intake flow rate, exhaust flow rate, average exhaust temperature, average bleed air temperature, average bleed air pressure, bleed air equivalent power, dynamic time, steady-state time, speed acceleration, and speed overshoot; the auxiliary power unit simulation model is used to calculate the performance values ​​under steady-state conditions, and the relative deviation between the test data and the model simulation calculation results is calculated.

[0129] Figure 3 This is a screenshot of the interactive interface of the present invention, where the required parameters are selected to draw the auxiliary power unit characteristic data image. The drawn image includes the following: Figure 4 The speed, total output shaft power, bleed air mass flow rate, and fuel mass flow rate of an embodiment of the present invention are normalized and changed over time, and also include a curve graph of total fuel consumption changing over time, a box graph of original data of steady-state bleed air mass flow rate, a box graph of original data of steady-state bleed air pressure, a box graph of original data of steady-state bleed air temperature, a box graph of original data of steady-state total output shaft power, a fingerprint graph of deviations between the average value of total bleed air pressure, the average value of bleed air temperature, the fuel mass flow rate, and the exhaust temperature MAX and the calculation results of the simulation model, a graph of exhaust temperature based on time proportion, and a graph of fuel mass flow rate based on time proportion.

[0130] like Figure 2 As shown, the present invention also provides a gas turbine auxiliary power unit test data processing system, the system comprising:

[0131] Interactive module: used to read the original test data, remove the test related information, and save the original test data in the form of a table.

[0132] Feature extraction module: Calculate the overall mean of the measured parameters one by one, and then use the sliding window method to calculate the window mean. If a certain number of window means exceed or fall below the overall mean, it means that the parameters have changed significantly and are marked as valid parameters. Difference is calculated for the time series of the test records to obtain the auxiliary power unit working time series and obtain the auxiliary power unit characteristic data table.

[0133] Interval division module: obtain the step moment of the position feedback parameter and mark it as the starting point of the auxiliary power unit working interval; obtain the speed stabilization moment and mark the end point of the auxiliary power unit working interval. The speed stability in the interval is marked as the steady-state interval, and the speed fluctuation is marked as the dynamic interval. The test data is divided into multiple different steady-state mode and dynamic mode working intervals; extract the data of each working interval into a separate feature data table and mark it in time order.

[0134] Data processing module: Dynamic interval data adopts sliding mean filtering, and steady-state interval data adopts sliding mean method to calculate steady-state performance in sections. For parameters with basically constant steady-state performance, the only data output of the interval that can characterize the steady-state performance is selected, and the sliding mean data is retained for parameters with continuously changing steady-state performance;

[0135] Interactive module: Rearrange the interval data into a complete data table in chronological order, and calculate the auxiliary power unit performance related parameters, including converted flow, total output shaft power, fuel consumption rate, thermal efficiency, intake flow rate, exhaust flow rate, average exhaust temperature, average bleed air temperature, average bleed air pressure, dynamic time, steady-state time, speed acceleration, and speed overshoot.

[0136] The interactive module receives the data obtained by other modules and the user input signal, draws the auxiliary power unit characteristic data image and saves it in the form of a picture to a specified location.

[0137] The interactive module uses the convertTDMS program module to read the test data file, remove the test related information, and save the original test data in the form of a table.

[0138] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, but the protection scope of the present invention is not limited thereto. Any technician familiar with the field can easily think of various equivalent modifications or substitutions within the technical scope disclosed by the present invention, and these modifications or substitutions should all be covered within the protection scope of the present invention.

Claims

1. A method for processing test data of a gas turbine auxiliary power unit, characterized in that: The following steps are involved: S1, read the test data of the gas turbine auxiliary power unit; S2, identifying the test data and calculating the signal characteristics of the test data; S3, dividing the steady-state interval and the transition-state interval of the test data according to the signal characteristics of the test data, and characterizing the switching points of the steady-state interval and the transition-state interval of the test data by using step characteristics of multiple signals; S4, calculating the interval signal characteristics of the test data according to the test data of each interval; S5, according to the interval signal characteristics of the test data, the data processing parameters of each interval are calculated, the test data of the transition state interval are smoothed by a filter, and the test data of the steady state interval are averaged to obtain a unique steady state value.

2. A method for processing test data of a gas turbine auxiliary power unit according to claim 1, characterized in that: In S2, the signal characteristics of the test data include step signal, variance, range, interquartile range, coefficient of variation, skewness and kurtosis; step signal and coefficient of variation are used to characterize the characteristics of different signals of the test data in a specific time period; Variance, range, interquartile range, skewness and kurtosis are used to characterize the degree of dispersion of experimental data.

3. A gas turbine auxiliary power unit test data processing method according to claim 2, characterized in that: The type and method of the interval signal features in S4 are the same as those of the signal features in S2.

4. A method for processing test data of a gas turbine auxiliary power unit according to claim 1, characterized in that: In S3, the steady-state interval and the transition state interval of the test data appear alternately. The steady-state interval is the stable working state of the gas turbine auxiliary power unit, and the transition state interval is the dynamic adjustment state between the two steady-state intervals.

5. A method for processing test data of a gas turbine auxiliary power unit according to claim 4, characterized in that: In S3, the method of characterizing the switching points between the steady-state interval and the transition state interval by the step characteristics of multiple signals is: After obtaining the step characteristics of all signals, they are arranged in chronological order, and the time when three or more signal step characteristics are close is taken as the interval division time. The transition state is characterized by the obvious change trend of the window mean. The window width is set to 10 consecutive sampling data. The obvious change trend is characterized by the following formula: K=sum([Mean win i -Mean|>Mean·p%),i=1,2,3... Among them, K represents the number of signal windows with obvious change trend of a certain parameter, Mean win,i represents the mean of the i-th window, Mean represents the interval mean of the parameter, and p represents the percentage of the mean deviation.

6. A method for processing test data of a gas turbine auxiliary power unit according to claim 5, characterized in that: K>5 means that the parameter has an obvious changing trend in the interval; if more than 10 parameters have an obvious changing trend in the same interval, the interval is considered to be a transition state interval; otherwise, it is a steady state interval.

7. A method for processing test data of a gas turbine auxiliary power unit according to claim 1, characterized in that: In S5, the experimental data in the transition state interval are smoothed by a low-pass filter, and the experimental data in the steady state interval are represented by a weighted average.

8. A gas turbine auxiliary power unit test data processing system, using a gas turbine auxiliary power unit test data processing method as claimed in any one of claims 1 to 7, comprising an interaction module, a signal feature extraction module, an interval division module and a data processing module; The interactive module reads the test data file according to the data file directory input by the user, draws the characteristic data chart according to the user input signal and the data processing result, and saves it to the specified location. The interactive module also saves the process data file and the result data file according to the data of the interval division module, the feature extraction module and the data processing module; draws the color data characteristic image for the result data file, and saves it to the specified location; The signal feature extraction module calculates the signal features based on the input test data. When the input data is all the test data, the calculated signal features are the overall signal features. When the input data is the processing interval data, the calculated signal features are the interval signal features. The interval division module divides the steady-state and transition-state intervals according to the signal characteristics; The data processing module calculates the data processing parameters of each interval according to the signal characteristics and obtains the data processing results.

Images