#### Abstract

This paper is concerned with designing a networked controller for a mixed flow two-spool turbofan aeroengine with aging and deterioration. Firstly, the state-space representation of the aeroengine considering aging and deterioration is identified, by which the engine system with aging and deterioration is modeled as an uncertain linear system. Then based on this uncertain linear system, theoretical results from the networked control systems and the regional pole assignment are introduced to formulate the networked engine control design in the form of linear matrix inequalities (LMIs). By solving these LMIs simultaneously, a networked engine controller is obtained which guarantees both the robustness against delay/dropout and the satisfactory dynamic performance. Finally, the proposed method is applied to an aerothermodynamic component-level engine simulator to demonstrate its validity and applicability. The corresponding delay/dropout margin is also calculated, which provides reference for the future development of the distributed engine control system.

#### 1. Introduction

The next revolution in turbine engine control systems would be the physical distribution of control functionality, known as distributed engine control (DEC) [1]. The concept of DEC is borrowed from the distributed control system in automated industry, which is featured by the application of advanced data buses, smart sensors, and smart actuators [2, 3]. The realization of DEC is beneficial to the implementation of the advanced control mode and control algorithm, such as active control, fuzzy control, and model-based control, which bring advantages of weight reduction, modularity, high reliability, high intelligence, etc. [4–10].

However, due to immatureness of the high-temperature electronics [11, 12], transition from the centralized control system to the fully distributed control system will be a gradual process. An intermediate distributed system approach is called the networked engine control (NEC) system, in which the data bus is introduced to replace the analog signal transmission. The A/D conversion functions (sampling, shaping, and quantization) of the centralized controller are moved into smart sensors while the control logic still remains as the central form (single controller).

Due to the introduction of the data bus, transmission delay and packet dropout are inevitable in a NEC system. Depending on how network-induced delays and packet dropouts are handled, a number of methods are developed on stability analysis and control design for networked control systems (NCSs). Among them, time-delay system approach and switched system approach are the two mainstream methods. The time-delay system approach yue 2004 state, Jiang 2008A, e.g., [13, 14], models the packet dropout as input-delay, and therefore, the maximum allowable number of consecutive packet dropouts is calculated. The switched system approach models the NCS as a switched system with arbitrary switching Zhang 2008 modelling, Zhang 2009 robust, Donkers 2011 stability, Kruszewski 2012 switched, e.g., [15–18], or known probability switching [19, 20], which is usually applied in the discrete-time domain.

As a safety-critical power device, it is crucial to figure out to what degree the control system could tolerate the delay and the packet dropout. To settle this problem, different kinds of NEC system models are established by scholars in the aeroengine control field. Constant delay and deterministic dropout are considered by [21, 22] who calculate the worst case bound on the number of consecutive dropouts. The case of stochastic dropout has also been considered. Reference [23] models the packet dropout as a Bernoulli process and calculates the maximum allowable packet dropping probability for the decentralized NEC system. Random delays in the NEC system are considered by [24, 25].

Although a lot of work has been done on the NEC systems, the requirement of the dynamic performance has not been fully considered. Since previous results are developed based on the Lyapunov stability theory, the overall system is Lyapunov stable under the network-induced factors but it may still have a very small decay rate if no constraints are imposed. Besides, the effect of the uncertainty has also not been considered in the design stage; thus, previous results are still far away from practical application. Based on the discussion above, the goal of this paper is to design a NEC system satisfying the following control requirements: (1)Robustness against the transmission delay and packet dropout(2)Satisfactory dynamical performance under different degrees of aging and deterioration(3)Zero tracking error

To satisfy these control requirements, theoretical results from the networked control systems and the regional pole placement [26, 27] which have been studied previously are considered as solutions. In this paper, these theories are further extended to an uncertain case and formulated in the form of LMIs which are numerically tractable. Besides, efforts are also made to establish the uncertain description for aeroengines to make the proposed method applicable.

This paper is organized as follows. Section 2 introduces the modeling process of the aeroengine. Section 3 presents the stability theory and controller design of the NEC system. In Section 4, the robust networked controller and its corresponding delay/dropout margin are calculated. Then simulations are performed to demonstrate the effectiveness of the proposed method. Discussions are also extended on further simulation results. Finally, Section 5 presents conclusions. The symbols used in this paper are listed in Table 1.

#### 2. Modeling of the Aeroengine

In this section, the working principle of the two-spool turbofan aeroengine and its commercial simulator is introduced first. Then by treating the complex simulator as a black box, a group of linear models is identified. These linear models are further redescribed as a nominal model with norm-bounded uncertainty to facilitate the controller design.

##### 2.1. Mixed Flow Two-Spool Turbofan Aeroengine

The cross-section of a mixed flow two-spool turbofan aeroengine is given in Figure 1. In a two-spool aeroengine, the compressors and turbines are linked by concentric shafts which rotate independently, namely, the low-pressure compressor (LPC) is driven by the low-pressure turbine (LPT) and the high-pressure compressor (HPC) is powered by the high-pressure turbine (HPT).

In [28], using GasTurb modeling software, an aerothermodynamic component-level aeroengine simulator is established for the mixed flow two-spool turbofan aeroengine with realistic structural parameter and features of compressors/turbines. The turbine engine modeling software, GasTurb, is a commercial software [29, 30] which is widely used for modeling and simulating gas turbine aeroengines. For the convenience in controller design, the aeroengine simulator developed in [28] is complied under the MATLAB/Simulink environment further. However, due to the limited space and the modeling of simulators being not the main focus of the paper, the modeling process of the simulator is not included in this paper. Readers who are interested could obtain a general knowledge about the modeling process from the gas turbine textbook (see chap. 2.5 of [28, 31]).

The nonlinear state and output equations of the aeroengine simulator are expressed as where and are nonlinear functions of the state vector , the input vector , and the health parameter vector as given in Table 2. The degradation is due to usage and aging results in the higher fuel consumption and exhaust gas temperature, which further yields a shorter component life and a higher cost. Generally, the degradation is reflected as changes in flow characteristics and efficiencies of the rotational components. Therefore, a group of multipliers, which are referred as “health parameters,” are introduced. For instance, denotes a new component while a value being lower than 1 denotes a degraded one. For the engine which needs an overhaul, the health parameters may degrade to a maximum of .

##### 2.2. State-Space Models Obtained by the System Identification

Considering the complexity of the simulator, the system identification technique is utilized to obtain a state-space representation. At the steady-state operating point , the simulator is identified as where symbol denotes the deviation from the steady-state point.

Take the “full health” status as an example. To apply the system identification technology, outputs of the aeroengine simulators and are first driven to their designed value given in Table 2; then two excitation signals are added separately to the input channels to generate and . These two excitation signals are taken, respectively, as where is a swept-frequency cosine signal with amplitude being 1% . As stated in chap. 3.1 of [31], frequencies of the dominant dynamics are generally not higher than 2 Hz for aeroengines. Therefore, the frequency of is set to increase linearly from 0 to Hz within s, which is adequate to extract the spectral characteristics of rotors. Moreover, zero-mean white Gaussian noises are also added to the excitation signals such that the signal-to-noise ratio (SNR) satisfies in which denotes the standard deviation.

The excitation signals and the deviations , , and generated from the aeroengine simulator are depicted in Figure 2, which have already been normalized by their designed value. With these two samples, the system matrices under the “full-health” status are identified by the method proposed in [32] as

For simplicity, it is assumed that all health parameters degrade simultaneously. Then by varying from to evenly and still maintaining at its designed value given in Table 3, a group of matrices , , , and is obtained via the system identification method above, which constitutes the set below

Due to the limited space, the values of these system matrices are not listed here.

##### 2.3. Derivation of the Norm-Bounded Uncertainty

In this paper, degradations of the health parameters (representing aging and deterioration) are viewed as the source of uncertainties for the identified state-space models. The degradation has two adverse effects on the identified linear models: On the one hand, the degradation is reflected in the drift of the steady-state value, and on the other hand, it causes the perturbation of the system matrices.

Since the identified system matrices , , , and vary with , they are supposed to be decomposed as in which , , and are norm-bounded uncertainties satisfying where and are constant matrices of appropriate dimensions. is the uncertain matrix with Lebesgue measurable elements that satisfy .

When obtaining the uncertain state-space representation (7), the difficulty lies in how to decompose a group of matrices further into a uniform norm-bounded description , namely, transforming to the set defined below.

To the best of our knowledge, there exists no transformation method such that equals .

Therefore, efforts are made to find a set which could encompass , namely, . This is handled via the concept of *interval matrix* introduced below.

*Definition 1 [33]. *With the given matrices , , and , notation is used to denote the set of matrix which satisfies , namely,
Then is called interval matrix, in which and are the lower- and upper-bound matrices, respectively.

Based on the lemma given below, the interval matrix can be converted equivalently to a matrix with norm-bounded uncertainty. Therefore, the uncertain representation is obtained by seeking a lower-bound matrix and an upper-bound matrix such that .

Lemma 1 [34]. *The interval matrix is equivalent to the set defined below.
in which matrices , , and are taken as
*

*Remark 1. *The equivalence of and has been proved in [34] and thus is omitted here. Readers could verify equalities and for ease of understanding.

With the given matrices , by selecting “tight” bound matrices and such that , matrices , , and are calculated according to Lemma 1 as
Similarly, matrices , , , , and in (7) are calculated as
Therefore, the set encompassing the set is obtained to redescribe a group of linear models (6) as a nominal matrix (9) with norm-bounded uncertainties.

#### 3. Networked Controller Design for the Aeroengine

In this section, the mathematical expression of the NEC system is presented which takes the transmission delay and the packet dropout into account simultaneously. Then a robust networked engine controller is designed with requirements of robustness against delay/dropout and the satisfactory dynamic performance being met with the NCS stability theory and the regional pole placement, respectively. The former will be presented in Section 3.2 while the latter in Section 3.3. Then by combining linear matrix inequalities (LMIs) from these two theories, the controller design is synthesized in Section 3.4.

##### 3.1. Networked Engine Control System

A schematic of the NEC system architecture is shown in Figure 3, where , , and denote the transmission delay between the sensors and the controller, the transmission delay between the controller and the actuators, and the calculation time of the controller, respectively.

To eliminate the potential tracking error due to reasons like modeling biases or plant disturbances, integrators are incorporated in the controller. Based on [35], by putting integrators at the output , format of the controller is taken as where is the state of the integrator satisfying

Supposed that the sensors are clock driven, the controller and actuators are event driven and the data are transmitted with a single packet and the integrated network-induced delay can be approximated as . When packet dropout occurs, the latest packet is utilized for computation and actuation. Then combining (2) with (15) and following the modeling process in [13], the closed-loop system is expressed as in which is the sampling period. is the transmission delay from the sampling instant to the actuating instant . is a nondecreasing sequence satisfying if the data packet disorder is ignored.

More specifically, means that no packet dropout occurs while means a packet loss in the transmission. By the definition of the symbol , it can be seen that (17) models the transmission delay and the packet dropout simultaneously.

Furthermore, by defining an augmented state vector as and considering the output reference to be zero here, the closed-loop system (17) can be expressed compactly as in which . and are system matrices that can still be decomposed into a norm-bounded form as with

The matrices , , , , and can be calculated based on (13) and (14), but they are not listed here due to their high dimensions.

##### 3.2. Stability Theory under Delay/Dropout Condition

In the following, the state feedback controller will be designed. Before we start, a frequently used lemma is introduced below which eliminates the uncertain term in LMIs effectively without introducing conservatism. Thus, by this lemma, existing theorems can be extended for uncertain systems straightforwardly.

Lemma 2 [36]. *With the given matrices , with appropriate dimensions, and a symmetric matrix , the matrix inequality
holds for all satisfying if and only if there exists a scalar such that
**By Lemma 2, the theorem introduced below guarantees the robustness against delay/dropout for the closed-loop system (18).*

Theorem 1 *For the given scalars and , the closed-loop system (18) is exponentially asymptotically stable with if there exist a scalar and matrices , , , , and with appropriate dimensions such that (23) holds.
where
*

*Proof. *It can be known from [14] that for the given scalars and , the closed-loop system (18) is exponentially asymptotically stable with if there exist matrices , , , **,** and with appropriate dimensions such that (25) holds.
Then by substituting into (25) and using Lemma 2 to matrix , it can be found that (25) holds if there exists a scalar satisying (26).
where
Then by Schur complements, (26) is equivalent to (23). The theorem is proved.

*Remark 2. *In Theorem 1, inequality describes the effect of the transmission delay and packet dropout simultaneously; hence, can be viewed as a delay/dropout margin.
(i)If no packet drops in the transmission , inequality holds. The maximum allowable transmission delay is calculated as
(ii)Assuming a constant delay in the transmission, inequality holds. Thus, the maximum allowable number of consecutive packet dropout is calculated as

##### 3.3. -Stability Theory

Though Theorem 1 guarantees that the closed-loop system is exponentially stable, the dynamic performance can still be poor since no constraints are imposed on the rate of decay. To obtain a satisfactory dynamic performance, a simple way is to assign the closed-loop poles of the augmented system (18), namely, the eigenvalues of , to some suitable regions on the left-half complex plane. To achieve this, the below is introduced first which is usually referred as the -stability theory.

Lemma 3 [37]. *With the given LMI region and matrix , the eigenvalues of lie in if and only if there exists a symmetric matrix such that
**Different types of regions, including the half-open plane, vertical strip area, and conic sector, can be covered by the above LMI description. Other complex regions can also be covered as long as they are convex. To guarantee a minimum decay rate and a minimum damping ratio, in this paper, the disk region is chosen as the expected eigenvalue region of which the center is and the radius is . This disk region is described by taking
**Then based on Lemma 3, the theorem below guarantees all the eigenvalues of in the disk region *

Theorem 2. *Given a disk region defined by (31), all eigenvalues of lie in if and only if there exist a symmetric matrix , a matrix with appropriate dimension, and a scalar such that
with the control gain matrix being calculated as .*

*Proof. *Being similar to the proof of Theorem 1, by substituting into and using Lemma 2 to eliminate the uncertain matrix , is equivalent to
Then by Schur complements and taking , (33) is equivalent to (32). The theorem is proved.

##### 3.4. Controller Design Synthesis

Based on the aforementioned work, the controller with strong delay/dropout robustness and satisfactory dynamic performance can be obtained through solving LMIs in Theorem 1 and Theorem 2 simultaneously. These LMIs are written together as in which and are matrices defined in (23) and (32), respectively, and , , , and are predesignated scalars.

Then for the given scalars , , and , the maximum delay/dropout margin can be calculated by a linear search. Correspondingly, the control gain matrix is calculated as .

#### 4. Application on the Turbofan Aeroengine Simulator

In this section, the controller design proposed above is applied on the turbofan aeroengine simulator by selecting appropriate closed-loop poles region. The controller is designed based on the uncertain model established in Section 2 and then verified on the simulator itself to validate the applicability of the proposed method.

##### 4.1. Controller Design

Once the uncertain state-space representation is ready, the center and the radii of the disk pole region are determined to calculate the control gain matrix. In practical turbofan applications, the percentage overshoot (PO) and the setting time of the engine control system are usually required to satisfy

By taking the damping ratio (dampling angle ) and the center , PO and are estimated by the formulae of the second-order system as in which . This fact indicts that the dynamic performance requirements above can be met if and . Based on this criterion, by taking and and searching between , the control gain matrix is calculated for (18) by solving LMIs in (34) as

The corresponding delay/dropout margin is calculated as . Then according to (18), the control gain matrix is implemented in the form of

In the NEC system, the transmission delay is bounded in one sampling period if a time-triggered data bus is utilized [23]. Therefore, taking as its worst value s, indicates a maximum allowable number of consecutive packet dropout based on Remark 2, namely, the NEC system is stable under a consecutive packet dropout being not severer than in every sampling cycles when the time delay s.

##### 4.2. Simulation

For comparison, another controller is calculated under the same condition but without considering the uncertainties. Comparison is performed between these two controllers to verify the effectiveness of our design.

By applying them on the component-level aeroengine simulator and varying the degree of deterioration , unit step response simulations are conducted on the output as shown in Figures 4 and 5. In each simulation, the health parameter is taken as , while the network condition is taken as the extreme situation s and as discussed above. The control input and the output in these figures have already been normalized by their designed values and .

Comparing Figure 4 (with uncertainties) with Figure 5 (without uncertainties), one can see that the settling time in Figure 4 is smaller than that in Figure 5. This shows that the controller design that is taking uncertainties into consideration is more robust and has a smaller settling time than the one without considering the uncertainties. Besides, it can be observed from Figure 4 that (1)The NEC system remains stable under the transmission delay and packet dropout described by the delay/dropout margin (2)The output presents a satisfactory as well as a consistent dynamic performance under different deteriorations. Besides, the response curves do not overshoot significantly since the disk region has a small imaginary part (damping angle )(3)The output could track the reference command without the steady-state error

From the points listed above, the robust networked controller designed in this paper has achieved the expected goal, namely, the three control requirements proposed in Introduction.

##### 4.3. Discussions

In this section, discussions are extended on the relationship between the delay/dropout margin and the disk region . Besides, the conservatism which is introduced when obtaining the norm-bounded uncertainties is also discussed.

###### 4.3.1. Relationship between the Delay/Dropout Margin and the Disk Region

As mentioned above, the expected disk region should satisfy and . Therefore, by varying the center and taking the radii as for a constant damping ratio, different gain matrices are and their delay/dropout margins are calculated in Table 3.

From Table 3, it also can be observed that the delay/dropout margins decrease rapidly while the pole region is being far away from the imaginary axis. Therefore, there should be a trade-off between the robustness against delay/dropout and the dynamical performance when designing the controller.

###### 4.3.2. Conservatism Introduced in the Modeling Process

In this paper, a data-based method is applied to obtain the norm-bounded uncertainties in Section 2.3. By this method, a group of system matrices ( in (6)) is redescribed as a nominal matrix with norm-bounded uncertainties ( in (9)). However, conservatism is introduced in this modeling process since is only a subset of . To demonstrate the conservatism graphically, integrators are incorporated as mentioned in Section 3.1 to augment and to and defined below as and in which is exactly the augmented system matrices in (18).

Based on the gain matrices of various disk regions obtained in the previous section, the closed-loop poles (namely, the eigenvalues of ) are calculated for and . As a result, the convex hulls of these poles are drawn in Figure 6, in which the pole region of is denoted in red while the counterpart of is denoted in blue. It can be seen that certain conservatism is introduced but is still at an acceptable level.

###### 4.3.3. Conservatism Introduced in the Controller Design Process

The conservatism is not just introduced by the modeling but also introduced by the controller. In this subsection, we combine the networked controller design theory [13] with the -stability theorem to design a controller for the uncertain system (18) as a comparison.

By using the same disk region selected in this paper and taking the lower bound of time delay as , the maximum allowable value of is calculated as based on [13]. By contrast, calculated in this paper is . Considering a dropout-free transmission, the gap means a larger maximum allowable delay by ms which is small enough. Therefore, only a little conservatism is introduced in the controller design process.

However, it should be pointed out that the theorem in [13] is in the form of bilinear matrix inequality (BMI), which has to be solved by a complex iterative algorithm. In contrast, the method proposed in this paper is numerically tractable. This makes the controller much easier to be calculated with existing LMI toolbox.

#### 5. Conclusion

In this paper, efforts are made to establish the uncertain description for aeroengines with aging and deterioration. Based on this model, the controller with both strong delay/dropout robustness and satisfactory dynamic performance is designed through solving linear matrix inequalities (LMIs). The validity and applicability of this approach are elaborated through simulations under the deteriorated and networked conditions. The corresponding delay/dropout margin is also calculated, which provides a reference for future development of the aeroengine-distributed control system. Admittedly, with some degree of conservatism, the controller design approach presented in this paper offered a feasible and practical mean to obtain a robust networked controller

#### Data Availability

Data can be provided if needed.

#### Conflicts of Interest

The authors declare that they have no conflicts of interest.

#### Acknowledgments

The authors would like to thank Assoc. Prof. Jiwei Wen in Jiangnan University and Mister Qian Feng in The University of Auckland for their beneficial suggestions.