Next Article in Journal
Trim Tab Flight Stabilisation System Performance Assessment under Degraded Actuator Speeds
Previous Article in Journal
Aerodynamic Exploration for Tandem Wings with Smooth or Corrugated Surfaces at Low Reynolds Number
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Robust LQR Design Method for the Aero-Engine Integral Constant Pressure Drop Control Valve with High Precision

1
School of Energy and Power Engineering, Beijing University of Aeronautics and Astronautics, Beijing 100191, China
2
AECC Guizhou Honglin Aeroengine Control Technology Co., Ltd., Guiyang 550009, China
*
Author to whom correspondence should be addressed.
Aerospace 2023, 10(5), 428; https://doi.org/10.3390/aerospace10050428
Submission received: 9 March 2023 / Revised: 25 April 2023 / Accepted: 26 April 2023 / Published: 30 April 2023
(This article belongs to the Section Aeronautics)

Abstract

:
The closed-loop constant pressure drop control valve is widely used in aero-engine fuel servo metering systems. However, the available constant pressure drop control valve cannot realize servo tracking without static error and, often, a high proportional gain is used to reduce the static error and improve the servo tracking performance, which reduces the stability margin. In this paper, an integral constant pressure drop control valve is designed, which consists of an integral controller and a stabilizing controller. Moreover, a robust LQR design method is proposed to complete the design task. Firstly, the controlled plant’s state–space model is derived, and the augmented model with tracking error is established based on the robust servo system design theory. Secondly, a servo controller with dual functions of integral control and stabilization control is constructed and decoupled, in which the stabilizing controller guarantees the asymptotic stability as well as the anti-disturbance performance, and the integral controller realizes the servo tracking without static error. Finally, based on the robust LQR design method, two key design parameters, including the integral control gain and the stabilization control gain, are designed to complete the design task. The simulation results indicate that, even when suffering 50 mm2 metered flow area step disturbance and 1 MPa inlet pressure step change, the designed integral constant pressure drop control valve can realize the function of servo tracking without static error. The static error is almost 0, the settling time is within 0.01 s, the overshoot is within 10%, and the phase margin is more than 55°.

1. Introduction

The aero-engine fuel servo metering system basically adopts the constant pressure drop control principle to complete fuel metering, except for very few applications of the variable pressure drop control structure, such as the fuel regulator in the Spey MK202 turbofan engine in the Rolls-Royce [1]. However, the available constant pressure drop control valve cannot realize servo tracking without static error and, often, a high proportional gain is used to reduce the static error and improve the servo tracking performance. The control theory explains that if there is no integral part in the closed-loop system, but only a proportional part, then the closed-loop system exhibits static error under any circumstances [2,3]. Although the high proportional gain can reduce the static error, it reduces the system stability and, even, causes aero-engine instability in the acceleration and deceleration process, such as thrust and speed swing [4]. In addition, as the high precision, high stability, and high robustness requirements in the modern aero-engine control system are proposed, it is crucial to improve the servo tracking ability and robustness of the constant pressure drop control valve.
Early studies about the constant pressure drop control valve primarily concentrated on frequency domain modeling and analysis. Specifically, based on the derived frequency domain models, the steady-state and dynamic characteristics of the whole system, as well as the stability variation under different structure parameters, were analyzed [5,6,7]. These studies offer some basic methods, however, due to the limitations to classical control theory and simulation technology, the established models are complex, the proposed methods cannot explain the design theory, and the conclusions cannot be verified because of lacking simulation works. In recent years, simulation technology has developed rapidly, and many studies relying on simulation technology have appeared, mainly concentrating on nonlinear modeling and simulation. For example, when the structure parameter values change, such as the hole diameter and the spring stiffness, the system performance variations are analyzed and simulated, and some guidance on the design structure parameters is offered [8,9,10,11,12]. However, these studies lack theoretical analysis and only rely on nonlinear simulation, causing difficulty in guiding the design process. Moreover, relying on the transfer function models, the system’s stability conditions are analyzed [13]. Indeed, the analysis results should be confirmed, because the model is oversimplified, and the study lacks a simulation comparison. Besides, some detailed studies, concentrating on the system’s characteristic changes with the return orifice profile structure, explain the correlation between the profile structure and the system control gain [14,15,16]. Nevertheless, these studies have limitations, because they only offer guidance for orifice profile design, rather than system design. Furthermore, there are some unique studies. For example, an ideal variable orifice is constructed and tested physically, and the test results indicate that the variable structure can enhance system performance [17]. Despite lacking design theory analysis, the design idea for the variable orifice structure in this study is helpful. In addition, through CFD simulation, the valve balance characteristic related to the flow force are explored [18,19,20,21]. The conclusions indicate that the flow force will affect the system performance, but it is not the most crucial aspect. Moreover, some physical tests have been executed to explore the influence of hysteresis on the speed fluctuation [22,23]. However, the theoretical analysis and results discussion in these studies are unclear, and they do not involve the design process. A new remarkable study, using the linear incremental analysis method, reveals the design theory for a constant pressure drop control valve, and proposes the frequency domain analysis and design methods, realizing a pretty good performance [24]. The study proposes efficient guidance measures and design methods; however, it is still limited to classical control design methods, and cannot realize servo tracking without static error.
These studies concentrate on classical control theory, rather than modern control theory, to study the high proportional gain constant pressure drop control valve, causing complex analysis and design processes. With the design object of servo tracking without static error, this paper uses the modern servo control theory to complete the design of the constant pressure drop control valve with static error-free tracking capability and provides a technical approach that can be used in engineering for high precision fuel metering. This paper makes the following contributions:
  • Firstly, based on the modern servo model design theory, an integral design structure for the constant pressure drop control valve is proposed, and a servo controller with dual functions of integral control and stabilization control is constructed;
  • Secondly, based on the decoupling design theory, the servo controller is decoupled. Where the stabilizing controller guarantees the asymptotic stability, as well as the disturbance rejection performance, and the integral controller realizes the command servo tracking without static error;
  • Finally, a robust LQR design method is proposed to design the control gains of the system, and it completes the design task well. The method is proven to guarantee fine performance and stability, as well as strong robust performance.
The chapters are arranged as follows. In Section 2, an integral design structure is constructed and the relevant design theory is derived. In Section 3, the servo controller is decoupled and realized. In Section 4, a robust LQR design method is proposed, and the design and implementation processes for the system are given. In Section 5, a design example is provided. In Section 6, the conclusions are presented.

2. Theoretical Design

A typical fuel metering system structure is shown in Figure 1. In general, the fuel flow metering equation is represented as Q = C q A 2 Δ P / ρ , in which the pressure drop Δ P is designed as a constant, and the required fuel flow Q is metered by the metered flow area A [13,24]. Specifically, the pressure drop control valve is used to guarantee the constant pressure drop Δ P , and the fuel metering valve is used to realize the control of the metered flow area A .
This paper concerns the design problem of the pressure drop control valve. The designed integral constant pressure drop control valve is shown in Figure 2, which is improved from the general pressure drop control valve [12,24].
The parameters are defined as follows: P S is the inlet pressure, P T is the return pressure, P C is the controlled pressure, P Z is the adjusting pressure, P O is the ejection pressure; A J is the metered flow area, A i n is the servo inlet flow area, A o u t is the servo outlet flow area, A Z is the adjusting flow area, A 1 is the fixed inlet flow area, A 2 is the fixed ejection flow area; V C is the controlled chamber volume, V Z is the adjusting chamber volume, V O is the ejection chamber volume; x z is the adjusting valve displacement, x y is the servo valve displacement; A y is the servo valve pressure bearing area, and A z x is the adjusting valve pressure bearing area [24].

2.1. Working Principle

For the integral constant pressure drop control valve, the fuel flow metering equation is correspondingly represented as Q = C q A J 2 ( P S P C ) / ρ , in which the pressure drop is ( P S P C ) , and the metered flow area is A J .
Generally, for the pressure drop control valve, the metered flow area A J is considered as the disturbance input and the inlet pressure P S is considered as the reference input, and the working principle of the system is described as: when the metered flow area A J disturbs or the inlet pressure P S changes, the controlled pressure P C changes, and the pressure drop ( P S P C ) deviates from the designed value. Instantaneously, because the pressure drop ( P Z P C ) changes, the adjusting valve moves and the adjusting flow area A Z changes, realizing a rapid regulation of the P C . Simultaneously, because the pressure drop ( P S P C ) changes, the servo valve also moves and the servo inlet and outlet flow area A i n and A o u t change, realizing a precise regulation of the P C . The dual function restores the pressure drop ( P S P C ) to the design value [24,25,26,27,28,29,30].
Hence, the design objective is to ensure the controlled pressure P C servo tracks the inlet pressure P S without static error, and rejects the disturbance related to the metered flow area A J .

2.2. State–Space Model of the Controlled Plant

The controlled plant is composed of the controlled pressure and the ejection pressure, and their pressure–flow nonlinear dynamic equations are:
d P C d t = B V C ( ( C q 1 A 1 + C q j A J ) 2 ( P S P C ) ρ C q z A Z 2 ( P C P O ) ρ + A y x ˙ y + A z x x ˙ z )
d P O d t = B V O ( C q z A Z 2 ( P C P O ) ρ C q o A 2 2 ( P O P T ) ρ )
where ρ is the oil density, B is the oil bulk modulus, and C q is the flow coefficient.
The calculation formula for C q is:
C q = C q max tanh ( 8 | Δ P | ρ ρ d h N u l a m c )
where C q max is the maximum flow coefficient, d h is the hydraulic diameter, N u is the absolute viscosity, and l a m c is the critical flow number [24].
The relevant linear dynamic differential equations are:
Δ P ˙ C = B V C ( ( K P J + K P Z ) Δ P C + K P Z Δ P O + K P J Δ P S + K A J Δ A J K A Z Δ A Z + A z x Δ x ˙ z + A y Δ x ˙ y )
Δ P ˙ O = B V O ( ( K P Z + K P T ) Δ P O + K P Z Δ P C + K A Z Δ A Z )
where K A J = C q j 2 ( P S P C ) ρ , K A Z = C q z 2 ( P C P O ) ρ , K P Z = C q z A Z 1 2 ρ ( P C P O ) , K P T = C q o A 2 1 2 ρ ( P O P T ) , K P J = ( C q j A J + C q 1 A 1 ) 1 2 ρ ( P S P C ) .
Finally, the state–space model of the controlled plant is:
x ˙ p = A p x p + B p u p + E p w y p = C p x p + D p u p
where x p = [ Δ P C Δ P O ] T , u p = [ Δ A Z Δ x ˙ y Δ x ˙ z ] T , w = [ Δ A J Δ P S ] T , y p = [ Δ P C Δ P O ] T ,
A p = [ B V C ( K P J + K P Z ) B V C K P Z B V O K P Z B V O ( K P Z + K P T ) ] B p = [ B V C K A Z B V C A y B V C A z x B V O K A Z 0 0 ] , E p = [ B V C K A J B V C K P J 0 0 ] C p = [ 1 0 0 1 ] , D p = [ 0 0 0 0 0 0 ]

2.3. Structure Design of the Servo Controller

As mentioned previously, the adjusting flow area A Z is regarded as the main control input. Considering the controlled pressure P C as the output feedback variable, according to the robust servo system design theory, as well as the steady-state deviation design method [2,3], the theoretical design diagram of the control architecture is shown in Figure 3.
Correspondingly, the design objective is to ensure the controlled pressure increment Δ P C servo tracks the inlet pressure increment Δ P S without static error, and rejects the disturbance related to the metered flow area increment Δ A J . The detailed design processes are detailed below.
The tracked input signal Δ P S and the rejected disturbance signal Δ A J are both constant value signals, and their first order differential equations are expressed as Δ A ˙ J = 0 and Δ P ˙ S = 0 , respectively. Since the input signals are constant, an integrator needs to be added to realize servo tracking without static error.
The tracking error signal is defined as:
e = Δ P C Δ P S
Then, the first order differential equation of the tracking error signal is:
e ˙ = Δ P ˙ C
A new state vector is defined as:
z = [ e ξ ] T
The augmented servo system design model is:
z ˙ = A ˜ z + B ˜ μ
where ξ = x ˙ p = [ Δ P ˙ C Δ P ˙ O ] T , μ = u ˙ p = [ Δ A ˙ Z Δ x ¨ y Δ x ¨ z ] T , and
A ˜ = [ 0 1 0 0 B V C ( K P J + K P Z ) B V C K P Z 0 B V O K P Z B V O ( K P Z + K P T ) ] , B ˜ = [ 0 0 0 B V C K A Z B V C A y B V C A z x B V O K A Z 0 0 ]
The output feedback vector is defined as y c = [ e Δ P ˙ C ] T , and the output model is:
y c = C ˜ z + D ˜ μ
where C ˜ = [ 1 0 0 0 1 0 ] , D ˜ = [ 0 0 0 0 0 0 ] .
Then, the output feedback control law of the servo model is:
μ = K c y c = [ k 1 k 2 ] [ e Δ P ˙ C ] = ( k 1 e + k 2 Δ P ˙ C )
where k 1 and k 2 are the elements of the output feedback control gain vector K c .
Then, the control input of the controlled plant can be expressed as:
u = μ   d t = ( k 1 e + k 2 Δ P ˙ C )   d t = k 1 e   d t k 2 Δ P C
Represented by the physical control input Δ A Z , it has:
Δ A Z = k 1 ( Δ P C Δ P S )   d t k 2 Δ P C = k 2 ( k 1 k 2 ( Δ P S Δ P C )   d t + Δ P C )
The specific structure design diagram of the servo controller is shown in Figure 4.
Since the above servo controller is coupled, it cannot be directly realized in a hydraulic control system. The decoupling problem of the servo controller needs to be considered, which is completed in the following chapters.

3. Decoupling Design of the Servo Controller

As mentioned in Section 2.3, the servo controller consists of two parts: an integral link and a stabilizing link. For hydraulic control systems, these two links can only be realized through hydraulic components. The design processes are detailed below.

3.1. Integral Controller

3.1.1. Characteristics Analysis of the Servo Valve

The nonlinear dynamic equation for the servo valve is:
A y P S A y P C F L M y x ¨ y K f 1 x ˙ y K 1 x y = 0
where M y is the mass, K 1 is the spring stiffness, K f 1 is the viscous friction coefficient, and F L is the initial spring force.
The relevant linear dynamic differential equation is:
A y ( Δ P S Δ P C ) M y Δ x ¨ y K f 1 Δ x ˙ y K 1 Δ x y = 0
Its steady-state characteristic is:
Δ x y = A y K 1 ( Δ P C Δ P S )
Obviously, the steady-state value of Δ x y reflects the tracking error.

3.1.2. Characteristics Analysis of the Adjusting Chamber

The pressure–flow nonlinear dynamic equation of the adjusting pressure is:
d P Z d t = B V Z ( C q i n A i n 2 ( P S P Z ) ρ C q o u t A o u t 2 ( P Z P T ) ρ A z x x ˙ z )
The relevant linear dynamic differential equation is:
Δ P ˙ Z = B V Z ( ( K P Y + K P T 2 ) Δ P Z + K A Y Δ A i n K A T 2 Δ A o u t + K P Y Δ P S A z x Δ x ˙ z )
where K A Y = C q i n 2 ( P S P Z ) ρ , K P Y = C q i n A i n 1 2 ρ ( P S P Z ) , K A T 2 = C q o u t 2 ( P Z P T ) ρ , K P T 2 = C q o u t A o u t 1 2 ρ ( P Z P T ) .
The functions A i n = f i n ( x u i n ) and A o u t = f o u t ( x u o u t ) are adopted to represent the geometry relationship of the servo valve inlet orifice and outlet orifice, respectively. According to their linearized gain characteristics:
Δ A i n = d f i n d x u i n Δ x u i n ,   Δ A o u t = d f o u t d x u o u t Δ x u o u t
Since Δ x u i n = Δ x u o u t = Δ x y , then
Δ A i n = d f i n d x u i n Δ x y ,   Δ A o u t = d f o u t d x u o u t Δ x y
If the steady-state servo inlet and outlet flow area are both designed as 0, that is A i n , 0 = A o u t , 0 = 0 , there are:
K P Y = 0 , K P T 2 = 0
Then, it has:
Δ P ˙ Z = B V Z ( ( K A Y d f i n d x u i n + K A T 2 d f o u t d x u o u t ) Δ x y A z x Δ x ˙ z )
where the generalized integral control gain is defined as:
K C = B V Z ( K A Y d f i n d x u i n + K A T 2 d f o u t d x u o u t )
Then, it has:
Δ P Z = Δ P ˙ Z d t = K C Δ x y d t B V Z A z x Δ x z
Because Δ x y reflects the tracking error ( Δ P C Δ P S ) , the integrator for the tracking error is embedded through Equation (27). Obviously, the integral function can be realized, and it ensures the performance of servo tracking without static error.
Finally, the state–space model of the integral controller is:
x ˙ i = A i x i + B i u i y i = C i x i + D i u i
where x i = [ Δ x y Δ x ˙ y Δ P Z ] T , u i = [ e Δ x ˙ z ] T , y i = [ Δ P Z ] ,
A i = [ 0 1 0 K 1 M y K f 1 M y 0 K C 0 0 ] , B i = [ 0 0 A y M y 0 0 B V Z A z x ] C i = [ 0 0 1 ] , D i = [ 0 0 ]

3.2. Stabilizing Controller

The nonlinear dynamic equation of the adjusting valve is:
A z x P Z A z x P C M z x ¨ z K f 2 x ˙ z K 2 x z F L 2 = 0
where M z is the mass, K 2 is the spring stiffness, K f 2 is the viscous friction coefficient, and F L 2 is the initial spring force.
The relevant linear dynamic differential equation is:
A z x ( Δ P Z Δ P C ) M z Δ x ¨ z K f 2 Δ x ˙ z K 2 Δ x z = 0
Its steady-state characteristic is:
Δ x z = A z x K 2 ( Δ P Z Δ P C )
Obviously, the steady-state value of Δ x z reflects the error e u .
The function A Z = f Z ( x u z ) is adopted to represent the geometry relationship of the adjusting valve orifice. Its linearized gain characteristic is:
Δ A Z = d f Z d x u z Δ x u z
Since Δ x u z = Δ x z , then:
Δ A Z = d f Z d x u z Δ x z
Because Δ x z reflects the error ( Δ P Z Δ P C ) , the negative feedback function is embedded through Equation (34), and it ensures the robust asymptotic stability, as well as the disturbance rejection performance.
The generalized stabilizing control gain is defined as:
K Z = d f Z d x u z
Finally, the state–space model of the stabilizing controller is:
x ˙ s = A s x s + B s u s y s = C s x s + D s u s
where x s = [ Δ x z Δ x ˙ z ] T , u s = [ e u ] , y s = [ Δ A z ] ,
A s = [ 0 1 K 2 M z K f 2 M z ] , B s = [ 0 A z x M z ] C s = [ K Z 0 ] , D s = [ 0 ]

3.3. State–Space Model of the Servo Controller

According to the aforementioned content, the structure of the decoupling controllers is shown in Figure 5.
Combining the integral link and the stabilization link, the state–space model of the servo controller can be expressed as:
x ˙ c = A c x c + B c 1 y + B c 2 r u = C c x c + D c 1 y + D c 2 r
where x c = [ Δ x y Δ x ˙ y Δ P Z Δ x z Δ x ˙ z ] T , y = [ Δ P C ] , r = [ Δ P S ] , u = [ Δ A Z ] ,
A c = [ 0 1 0 0 0 K 1 M y K f 1 M y 0 0 0 K C 0 0 0 B V Z A z x 0 0 0 0 1 0 0 A z x M z K 2 M z K f 2 M z ] , B c 1 = [ 0 A y M y 0 0 A z x M z ] , B c 2 = [ 0 A y M y 0 0 0 ] C c = [ 0 0 0 K Z 0 ] , D c 1 = [ 0 ] , D c 2 = [ 0 ]
Then, the integrated structure of the servo controller is shown in Figure 6.

4. Robust LQR Design Method

As mentioned in Section 3.1.2, the steady-state flow area of the servo inlet and outlet orifices are both designed as 0. It follows that once the flow area variation gradient of the orifice is determined, the generalized integral control gain K C is a constant value at any steady-state working point and cannot be designed arbitrarily.
When combining the servo controller and the controlled plant, embedding the generalized integral control gain K C in the state–space model, and considering the generalized stabilizing control gain K Z as the static output feedback gain, the open-loop augmented state–space model of the whole system can be expressed as:
x ˙ = A x + B u + E w y = C x + D u
where x = [ Δ x y Δ x ˙ y Δ P Z Δ x z Δ x ˙ z Δ P C Δ P O ] T , u = [ Δ A Z ] , y = [ Δ x z ] ,
A = [ 0 1 0 0 0 0 0 K 1 M y K f 1 M y 0 0 0 A y M y 0 K C 0 0 0 B V Z A z x 0 0 0 0 0 0 1 0 0 0 0 A z x M z K 2 M z K f 2 M z A z x M z 0 0 B V C A y 0 0 B V C A z x B V C ( K P J + K P Z ) B V C K P Z 0 0 0 0 0 B V O K P Z B V O ( K P Z + K P T ) ] B = [ 0 0 0 0 0 B V C K A Z B V O K A Z ] , E = [ 0 0 0 A y M y 0 0 0 0 0 0 B V C K A J B V C K P J 0 0 ] C = [ 0 0 0 1 0 0 0 ] , D = [ 0 ]
This is a typical static output feedback design problem, and its output feedback control law is:
u = K Z y
The topology of the output feedback control system is shown in Figure 7.

4.1. Output Feedback Control Gains Design

4.1.1. Integral Control Gain

Defined as:
d f i n d x u i n [ k i n , min , k i n , max ] , d f o u t d x u o u t [ k o u t , min , k o u t , max ]
where k i n , min , k i n , max , k o u t , min and k o u t , max is the extremum values of the flow area variation gradient of the servo inlet orifice and outlet orifice, respectively.
According to Equation (26), it has:
K C [ B V Z ( K A Y k i n , min + K A T 2 k o u t , min ) , B V Z ( K A Y k i n , max + K A T 2 k o u t , max ) ]
According to the structural parameter limitations, appropriate values of the flow area variation gradient can be designed, then the integral control gain can be determined.

4.1.2. Stabilization Control Gain

The quadratic optimization objective for Equation (40) is defined as:
J = 1 2 0 [ y T Q y + u T R u ] d t
where Q = Q T 0 , R = R T > 0 .
Then, the LQR state feedback control law is:
u = K x x = R 1 B T P x
where K x is the state feedback control gain vector, and P is the solution of an ARE algebraic equation, expressed as:
P A + A T P + C T Q C P B R 1 B T P = 0
By designing the weight matrices Q and R to optimize the quadratic objective, a positive definite solution P of the ARE equation can be obtained. Then, the above system yields satisfactory closed-loop reference dynamics, described by:
x ˙ = ( A B K x ) x = F x
Using the output feedback control law expressed by Equation (42), the dominant eigenvalue Λ 1 and its associated eigenvector X 1 of the state feedback design can be retained, and the static output feedback gain is given by [3]:
K Z = K x X 1 ( C X 1 ) 1
where the eigenvector X 1 and the eigenvalue Λ 1 satisfy the eigen equation for the state feedback system, described as:
F X 1 = X 1 Λ 1
Finally, the closed-loop system is calculated as:
x ˙ = ( A B K Z C ) x = A c l x

4.2. Design and Implementation Method of the Servo Controller

Assuming that the performance requirements are:
  • Steady-state requirement: The steady-state pressure drop is designed as P e , and the phase margin should be more than N°;
  • Dynamic requirement: The settling time should be within the t s , and the overshoot should be within the σ .

4.2.1. Control Gains Design

The steady-state flow balance equation of the controlled plant is:
( C q j A J + C q 1 A 1 ) 2 ( P S P C ) ρ = C q z A Z 2 ( P C P O ) ρ = C q o A 2 2 ( P O P T ) ρ
At a steady-state working point, if the inlet pressure is represented as P S , i , and the metered flow area is represented as A J , i , then:
P C , i = P S , i P e
P O , i = ( C q j , i A J , i + C q 1 , i A 1 C q o , i A 2 ) 2 ( P S , i P C , i ) + P T
A Z , i = ( C q j , i A J , i + C q 1 , i A 1 C q z , i ) ( P S , i P C , i ) ( P C , i P O , i )
Subsequently, the parameters K A J , i , K A Z , i , K P Z , i , K P J , i and K P T , i can be calculated.
The steady-state balance equation of the adjusting valve is:
0 = A z x P C , i A z x P Z , i K 2 x s z d , i
Then:
P Z , i = P C , i + K 2 A z x x s z d , i
Subsequently, the parameters K A Y , i and K A T 2 , i can be calculated.
When designing appropriate values of the flow area variation gradient d f i n d x u i n and d f o u t d x u o u t , then the integral control gain K C can be determined according to Equation (26).
When designing optimized weight matrices Q and R that meet the above performance requirements, a positive definite solution P i of the ARE equation can be obtained, and the state feedback control gain vector K x , i can be obtained, then the stabilization control gain is expressed as:
K Z , i = K x , i X 1 ( C X 1 ) 1

4.2.2. Orifice Geometry Relationship Design

A design formula is provided to deal with the geometry relationship design problem of the flow area increment ( A Z , i A Z , i 1 ) and the orifice underlap increment Δ x u z , i , expressed by:
Δ x u z , i = ( A Z , i A Z , i 1 ) K Z , i 1 , i = 2 , 3 , , N
where N is the number of the designed steady-state working points [24].
Assuming that the adjusting valve steady-state spring compression at the (i − 1)-th steady-state working point is x s z d , i 1 , then the calculation formula of the steady-state spring compression at the i-th steady-state working point is provided as:
x s z d , i = x s z d , i 1 + Δ x u z , i

4.2.3. Valve Initial Parameters Design

(1)
Adjusting valve
The initial underlap and the underlap at the first steady-state working point are designed as x u z , 0 and x u z , 1 , respectively. Since Δ x u z = Δ x s z d , the initial spring compression x s z d , 0 is calculated as:
x s z d , 0 = x s z d , 1 ( x u z , 1 x u z , 0 )
(2)
Servo valve
The steady-state balance equation of the servo valve is:
0 = A y P S , i A y P C , i K x s c d , i
Then, at the first steady-state working point, the steady-state spring compression is:
x s c d , 1 = A y K P e
The initial underlap and the underlap at the first steady-state working point is designed as x u i n , 0 and x u i n , 1 , respectively. Since Δ x u i n = Δ x s c d , the initial spring compression x s c d , 0 is calculated as:
x s c d , 0 = x s c d , 1 ( x u i n , 1 x u i n , 0 )

5. Design Example

The structural parameters of the pressure drop control valve are shown in Table 1.
The input conditions include:
  • The inlet pressure P S is within [3, 9] the MPa;
  • The metered flow area A J is within [10, 240] × 10−6 m2.
The design objectives include:
  • Geometry design of the adjusting orifice A Z = f Z ( x u z ) ;
  • Geometry design of the servo orifices A i n = f i n ( x u i n ) and A o u t = f o u t ( x u o u t ) .
The performance requirements include:
  • The rated pressure drop P e is 0.92 MPa, and the variation range is within 0.01 MPa;
  • The settling time is within 0.01 s, and the overshoot is within 10%;
  • The phase margin should be more than 50°.

5.1. Dynamic Design

5.1.1. First Steady-State Working Point

The input conditions of the first steady-state working point are: the metered flow area A J , 1 is 10 × 10−6 m2, and the inlet pressure P S , 1 is 9 MPa. According to Equations (53)–(55), then:
P C , 1 = P S , 1 Δ P e = 8.08   MPa
P O , 1 = ( C q j , 1 A J , 1 + C q 1 , 1 A 1 C q o , 1 A 2 ) 2 ( P S , 1 P C , 1 ) + P T = 0 . 20265   MPa
A Z , 1 = ( C q j , 1 A J , 1 + C q 1 , 1 A 1 C q z , 1 ) ( P S , 1 P C , 1 ) ( P C , 1 P O , 1 ) = 3 . 48627 10 6   m 2
Then, the parameters K A J , 1 , K A Z , 1 , K P Z , 1 , K P J , 1 and K P T , 1 are calculated as:
K A J , 1 = C q j , 1 2 ( P S , 1 P C , 1 ) ρ = 33 . 99849
K A Z , 1 = C q z , 1 2 ( P C , 1 P O , 1 ) ρ = 99 . 48459
K P Z , 1 = C q z , 1 A Z , 1 1 2 ρ ( P C , 1 P O , 1 ) = 2 . 20144 10 11
K P J , 1 = ( C q j , 1 A J , 1 + C q 1 , 1 A 1 ) 1 2 ρ ( P S , 1 P C , 1 ) = 1 . 88494 10 10
K P T , 1 = C q o , 1 A 2 1 2 ρ ( P O , 1 P T ) = 6 . 54328 10 8
Designing the spring compression x s z d , 1 as 10 mm. According to Equation (57), it has:
P Z , 1 = P C , 1 + K 2 A z x x s z d , 1 = 8 . 22737   MPa
Then, the parameters K A Y , 1 and K A T 2 , 1 are calculated as:
K A Y , 1 = C q i n , 1 2 ( P S , 1 P Z , 1 ) ρ = 31 . 15679
K A T 2 , 1 = C q o u t , 1 2 ( P Z , 1 P T ) ρ = 100 . 42741
When designing the weight matrices Q as 2.041 and R as 2000, the different control gain design schemes are shown in Table 2.

5.1.2. Other Steady-State Working Points

In this paper, the flow area variation gradient of the servo inlet orifice and outlet orifice d f i n d x u i n and d f o u t d x u o u t are both designed as 0.03, and the 4 × 5 steady-state working points are selected as the designed points. Subsequently, when executing the design processes shown in Section 5.1.1, the stabilization control gains that meet the above performance requirements can be obtained, as shown in Table 3.

5.1.3. Valve Initial Parameters Design

(1)
Adjusting valve
The initial underlap x u z , 0 and the underlap at the first steady-state working point x u z , 1 are both designed as 0.1 mm. According to Equation (61), it has:
x s z d , 0 = x s z d , 1 ( x u z , 1 x u z , 0 ) = 10   mm
(2)
Servo valve
According to Equation (63), at the first steady-state working point, the steady-state spring compression is:
x s c d , 1 = A y K P e = 23 . 39155   mm
The initial underlap x u i n , 0 and the underlap at the first steady-state working point x u i n , 1 are both designed as 0 mm. According to Equation (64), it has:
x s c d , 0 = x s c d , 1 ( x u i n , 1 x u i n , 0 ) = 23.39155   mm
Finally, the summary of the designed initial parameters is shown in Table 4.

5.1.4. Orifice Geometry Relationship Design

According to Equations (59) and (60), the geometry relationship value pair of the adjusting orifice is shown in Table 5.
Moreover, the geometry relationship diagram of the adjusting orifice is shown in Figure 8.
In addition, the geometry relationship diagrams of the servo inlet orifice and outlet orifice are shown in Figure 9a,b.

5.2. Simulation and Discussion

A nonlinear model is established based on the AMESim software, and the structure diagram is shown in Figure 10.
According to Table 1, the relevant structure parameters are set, and according to Table 4, the relevant initial parameters are set. In addition, the remaining structure parameters and the simulation parameters are set to the default values of AMESim.

5.2.1. Simulation

Executing the following simulation tasks, the performance of the designed system can be verified. The input conditions and the simulation results are detailed below.
  • The step signal of the inlet pressure P S is shown in Figure 11. Besides, the steady-state working points of the metered flow area A J are designed as 10 mm2, 50 mm2, 100 mm2, 150 mm2, 200 mm2, and 240 mm2, respectively. The simulation results are shown in Figure 12.
  • The step signal of the metered flow area A J is shown in Figure 13. Besides, the steady-state working points of the inlet pressure P S are designed as 3 MPa, 4 MPa, 5 MPa, 6 MPa, 7 MPa, 8 MPa, and 9 MPa, respectively. The simulation results are shown in Figure 14.

5.2.2. Discussion

The simulation results shown in Figure 12 and Figure 14 indicate that:
  • Despite suffering from strong step inputs, the controlled pressure drop is always 0.92 MPa, and the static error is almost 0. Evidently, the designed control device has the ability to perform servo tracking without static error. It follows that the theoretical design architecture is relevant and the LQR design method is effective;
  • During each transient process, the settling time is within 0.01s and the dynamic overshoot is within 10%. Obviously, the dynamic performances match the design requirements. It follows that the derived models are precise and the designed weight matrices, Q and R, are reasonable.
In addition, although the proposed method is based on the linear model, the designed system still has good performance and strong robustness when faced with nonlinear effects, such as leakage. The verification processes are detailed below.
According to the aforementioned content, the adjusting chamber is almost dead at each steady-state working point, hence it is most likely to leak. Assuming that the adjusting chamber leaks, a nonlinear model can be established, and the structure diagram is shown in Figure 15.
The leakage area is represented as A L , and it is designed as 0 mm2, 0.5 mm2, 1.0 mm2, 1.5 mm2, 2.0 mm2, and 2.5 mm2, respectively. Selecting some steady-state working points for verification and executing the following simulation tasks:
  • The step signal of the inlet pressure P S is still shown in Figure 11 and the metered flow area A J is designed as 150 mm2. The simulation results are shown in Figure 16;
  • The step signal of the metered flow area A J is still shown in Figure 13 and the inlet pressure P S is designed as 7 MPa. The simulation results are shown in Figure 17.
The simulation results shown in Figure 16 and Figure 17 indicate that:
  • When the adjusting chamber leaks, the pressure drop in the system deviates from the design value, and as the leakage area increases, the deviation increases; however, even if the leakage area reaches 2.5 mm2, the deviation is still within 0.01 MPa. Thus, the steady-state performance is acceptable;
  • Besides, during each transient process, the settling time is still within 0.01 s and the dynamic overshoot is still within 10%. Obviously, the dynamic performances match the design requirements and the leakage has a small impact on the dynamic performances.
To enhance the reliability of the simulation results, the same input conditions as mentioned in Section 5.2.1 are used and the leakage area is set as 2.5 mm2. When executing corresponding simulation tasks, the simulation results are shown in Figure 18 and Figure 19.
The simulation results indicate that the pressure drop in the system deviates from the design value, but the deviation is all within 0.01 MPa, and the designed system still has good steady-state and dynamic performance.
These simulation results show that, when faced with nonlinear effects, the designed system has strong robustness, and the proposed method is robust and effective.

6. Conclusions

In this paper, with the design object of achieving servo tracking without static error, an integral constant pressure drop control valve is designed, and a robust LQR design method is proposed. The conclusions are as follows:
  • Based on the servo system design theory, a servo control architecture for the pressure drop control valve is constructed and implemented, which can clearly explain the design theory of the system. In addition, the control architecture clearly displays the key structural design parameters, including the generalized stabilization control gain and the generalized servo control gain. Compared with classic design methods, the proposed design architecture and design method are more illustrative;
  • Based on the output feedback design theory, the robust LQR design method can realize the design of the key structural design parameters effectively and obtain optimized structural parameters to guarantee high precision and high robustness in performance. The proposed method provides more accurate guidance for the design of the structural parameters and improves the design efficiency;
  • The simulation results show that the designed integral constant pressure drop control valve has dual control functions of integral control and stabilization control and can realize tracking without static error and pretty good dynamic performance. Besides, when faced with nonlinear effects, the designed system still has good performance and strong robustness. Evidently, the proposed design method is robust and effective, and can also be used in the design process of other fuel system components.
However, nonlinearities also play an essential role in the design process, and it is meaningful to study how nonlinear characteristics affect system performance in a future study. Besides, if executing physical tests, the proposed design methods may not perform as well as the simulation results, which should also be tested in future work.

Author Contributions

Methodology, W.Z.; formal analysis, W.Z.; writing—original draft preparation, W.Z.; writing—review and editing, X.W.; validation, Y.L.; funding acquisition, Z.Z. and L.T. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the National Science and Technology Major Project (J2019-V-0010-0104).

Data Availability Statement

The data used to support the findings in this paper are contained in the text.

Conflicts of Interest

There are no conflict of interest.

References

  1. Wu, D.G. Research on the Performance of Turbofan Engine and Its System; National Defense Industry Press: Beijing, China, 1986; pp. 161–201. [Google Scholar]
  2. Wang, X.; Yang, S.B.; Zhu, M.Y.; Kong, X.X. Aeroengine Control Principles; Science Press: Beijing, China, 2021; pp. 128–257. [Google Scholar]
  3. Lavretsky, E.; Wise, K.A. Robust and Adaptive Control with Aerospace Applications; Springer: London, UK, 2012; pp. 27–72. [Google Scholar] [CrossRef]
  4. Yang, F.; Wang, X.; Cheng, T.; Liu, X. Dynamic characteristics analysis of a pressure differential valve. Aeroengine 2015, 41, 44–50. [Google Scholar] [CrossRef]
  5. Ma, C.Y. The analysis and design of hydraulic pressure-reducing valves. J. Eng. Ind. 1967, 89, 301. [Google Scholar] [CrossRef]
  6. Fan, R.; Zhang, M. The establishment of pilot-operated relief valve’s dynamic mathematic model and the dynamic properties analysis. J. ZhengZhou Text. Inst. 1997, 3, 58–61. [Google Scholar]
  7. Wu, D.; Burton, R.; Schoenau, G.; Bitner, D. Analysis of a pressure—Compensated flow control valve. J. Dyn. Syst. Meas. Control. 2007, 129, 203–211. [Google Scholar] [CrossRef]
  8. Wen, Y.J. Analysis of characteristics of oil return differential pressure valve with EASY5. In Proceedings of the 11th Symposium on Automatic Engine Control of CAA, Beijing, China, 8–10 November 2002; pp. 21–26. [Google Scholar]
  9. Hong, W.; Liu, H.L.; Wang, G.Z.; Qin, J. Research on pressure characteristics of relief valve without pressure overshoot. Chin. Hydraul. Pneum. 2012, 10, 104–106. [Google Scholar] [CrossRef]
  10. Shang, Y.; Guo, Y.Q.; Wang, L. Study of impact of design parameter of differential pressure controller on fuel metering system. Aeronaut. Manuf. Technol. 2013, 6, 89–91. [Google Scholar] [CrossRef]
  11. Wang, B.; Zhao, H.C.; Ye, Z.F. AMESim simulation of afterburning metering unit for fuel system. Aeroengine 2014, 40, 62–66. [Google Scholar] [CrossRef]
  12. Hang, J.; Li, Y.Y.; Yang, L.M.; Li, Y.H. Design and Simulation of Large Flowrate Fuel Metering Valve of Aero engine Based on AMESim. In Proceedings of the 2020 15th IEEE Conference on Industrial Electronics and Applications (ICIEA), Kristiansand, Norway, 9–13 November 2020; IEEE: Piscataway, NJ, USA, 2020. [Google Scholar] [CrossRef]
  13. Wei, Y.Y.; Wang, H.Y.; Miao, W.B. Analysis on modeling of constant pressure difference valve for a turboshaft engine. Aeroengine 2014, 40, 75–78. [Google Scholar] [CrossRef]
  14. Zeng, D.T.; Wang, X. Design and analysis of characteristics of damping hole for a fuel metering valve. In Proceedings of the 2010 International Conference on Mechanical and Electrical Technology, Kyoto, Japan, 1–3 August 2010; IEEE: Piscataway, NJ, USA, 2010. [Google Scholar] [CrossRef]
  15. Zeng, D.T.; Wang, X.; Tan, D.L.; Xu, M. Fuel scavenger contour performance analysis of fuel metering devices. Aeroengine 2010, 36, 38. [Google Scholar] [CrossRef]
  16. Zeng, D.T.; Wang, X.; Tan, D.L. Effects of fuel returned shape on metering devices characteristics. Aeroengine 2012, 38, 46. [Google Scholar] [CrossRef]
  17. Maiti, R.; Pan, S.; Bera, D. Analysis of a load sensing hydraulic flow control valve. Proc. JFPS Int. Symp. Fluid Power 1996, 1996, 307–312. [Google Scholar] [CrossRef]
  18. Amirante, R.; Vescovo, G.D.; Lippolis, A. Flow forces analysis of an open center hydraulic directional control valve sliding spool. Energy Convers. Manag. 2006, 47, 114–131. [Google Scholar] [CrossRef]
  19. Amirante, R.; Vescovo, G.D.; Lippolis, A. Evaluation of the flow forces on an open centre directional control valve by means of a computational fluid dynamic analysis. Energy Convers. Manag. 2006, 47, 1748–1760. [Google Scholar] [CrossRef]
  20. Okungbowa, B.; Stanley, N. CFD Analysis of Steady State Flow Reaction Forces in a Rim Spool Valve. Master’s Thesis, University of Saskatchewan, Saskatoon, SK, Canada, 2006. [Google Scholar]
  21. Valdes, J.R.; Miana, M.J.; Nunez, J.L. Reduced order model for estimation of fluid flow and flow forces in hydraulic proportional valves. Energy Convers. Manag. 2008, 49, 1517–1529. [Google Scholar] [CrossRef]
  22. Li, Z.; Guo, Y.Q.; Liao, G.H. Structure design and performance calculation of a differential pressure measuring device. In Proceedings of the 12th Engine Automatic Control Academic Conference of CAAC, Hong Kong, China, 1–5 November 2004; pp. 158–162. [Google Scholar]
  23. Deng, Z.J.; Guo, L.Y. Analysis of speed fluctuation in bench test of aero-engine numerical control system. China Sci. Technol. Overv. 2020, 15, 63–65. [Google Scholar]
  24. Zhao, W.; Wang, X.; Long, Y.; Zhou, Z.; Tian, L. Frequency Domain Design Method of the Aeroengine Fuel Servo Constant Pressure Difference Control System with High Performance. Aerospace 2022, 9, 775. [Google Scholar] [CrossRef]
  25. Merritt, H.E. Hydraulic Control Systems; John Wiley: New York, NY, USA, 1967; pp. 360–374. [Google Scholar] [CrossRef]
  26. Fitch, E.C.; Hong, I.T. Hydraulic Component Design and Selection; Bardyn Incorporation: Etobicoke, ON, Canada, 2004; pp. 205–213. [Google Scholar]
  27. Li, C.G.; He, Y.M. Modeling and Simulation Analysis of Hydraulic System; Aviation Industry Press: Beijing, China, 2008; pp. 2–11. [Google Scholar]
  28. Wang, Y.; Fan, D.; Zhang, C.; Peng, K.; Shi, D.Y. Design and analysis of the variable pressure-drop fuel metering device. In Proceedings of the 36th Chinese Control Conference, Dalian, China, 26–28 July 2017; pp. 6434–6439. [Google Scholar]
  29. Agh, S.M.; Pirkandi, J.; Mahmoodi, M.; Jahromi, M. Optimum design simulation and test of a new flow control valve with an electronic actuator for turbine engine fuel control system. Flow Meas. Instrum. 2019, 65, 65–77. [Google Scholar] [CrossRef]
  30. Yuan, Y.; Zhang, T.H.; Lin, Z.L.; Zhang, J.M. An investigation into factors determining the metering performance of a fuel control unit in an aero engine. Flow Meas. Instrum. 2020, 71, 101672. [Google Scholar] [CrossRef]
Figure 1. Structure diagram of the fuel metering system.
Figure 1. Structure diagram of the fuel metering system.
Aerospace 10 00428 g001
Figure 2. Structure diagram of the integral constant pressure drop control valve.
Figure 2. Structure diagram of the integral constant pressure drop control valve.
Aerospace 10 00428 g002
Figure 3. Theoretical design diagram of the control architecture.
Figure 3. Theoretical design diagram of the control architecture.
Aerospace 10 00428 g003
Figure 4. Structure design diagram of the servo controller.
Figure 4. Structure design diagram of the servo controller.
Aerospace 10 00428 g004
Figure 5. Structure diagram of the decoupling controllers.
Figure 5. Structure diagram of the decoupling controllers.
Aerospace 10 00428 g005
Figure 6. Integrated structure diagram of the servo controller.
Figure 6. Integrated structure diagram of the servo controller.
Aerospace 10 00428 g006
Figure 7. Topology of the output feedback control system.
Figure 7. Topology of the output feedback control system.
Aerospace 10 00428 g007
Figure 8. Geometry relationship diagram of the adjusting orifice.
Figure 8. Geometry relationship diagram of the adjusting orifice.
Aerospace 10 00428 g008
Figure 9. (a) Geometry relationship diagram of the servo inlet orifice; (b) geometry relationship diagram of the servo outlet orifice.
Figure 9. (a) Geometry relationship diagram of the servo inlet orifice; (b) geometry relationship diagram of the servo outlet orifice.
Aerospace 10 00428 g009
Figure 10. Structure diagram of the nonlinear model.
Figure 10. Structure diagram of the nonlinear model.
Aerospace 10 00428 g010
Figure 11. Step signal of the inlet pressure.
Figure 11. Step signal of the inlet pressure.
Aerospace 10 00428 g011
Figure 12. Response curves when the inlet pressure changes.
Figure 12. Response curves when the inlet pressure changes.
Aerospace 10 00428 g012
Figure 13. Step signal of the metered flow area.
Figure 13. Step signal of the metered flow area.
Aerospace 10 00428 g013
Figure 14. Response curves when the metered flow area disturbs.
Figure 14. Response curves when the metered flow area disturbs.
Aerospace 10 00428 g014
Figure 15. Structure diagram of the nonlinear model when the adjusting chamber leaks.
Figure 15. Structure diagram of the nonlinear model when the adjusting chamber leaks.
Aerospace 10 00428 g015
Figure 16. Response curves of the step in inlet pressure when the adjusting chamber leaks.
Figure 16. Response curves of the step in inlet pressure when the adjusting chamber leaks.
Aerospace 10 00428 g016
Figure 17. Response curves of the step in metered flow area when the adjusting chamber leaks.
Figure 17. Response curves of the step in metered flow area when the adjusting chamber leaks.
Aerospace 10 00428 g017
Figure 18. Response curves when the inlet pressure changes in the 2.5 mm2 leakage area.
Figure 18. Response curves when the inlet pressure changes in the 2.5 mm2 leakage area.
Aerospace 10 00428 g018
Figure 19. Response curves when the metered flow area is disturbed in the 2.5 mm2 leakage area.
Figure 19. Response curves when the metered flow area is disturbed in the 2.5 mm2 leakage area.
Aerospace 10 00428 g019
Table 1. Structural parameters of the pressure drop control valve.
Table 1. Structural parameters of the pressure drop control valve.
Parameter/UnitValueParameter/UnitValue
M y /Kg0.08 K 1 /(N/m)4 × 104
M z /Kg0.05 K f 1 /(N/(m/s))200
d y /m0.036 K 2 /(N/m)1.5 × 104
d z /m0.036 K f 2 /(N/(m/s))200
V C /m32 × 10−6 ρ /(Kg/m3)780
V Z /m32 × 10−6 B /MPa1.7 × 103
V O /m34.9087 × 10−4 C q max 0.7
A 1 /m22.8274 × 10−7 N u /Pas0.051
A 2 /m21.9007 × 10−4 l a m c 1 × 103
P T /MPa0.2
Table 2. Different control gain design schemes at the first steady-state working point.
Table 2. Different control gain design schemes at the first steady-state working point.
d f i n d x u i n / d f o u t d x u o u t K Z , 1 Phase/°
0.0050.021360.4
0.0100.030759.8
0.0150.030359.0
0.0200.030258.1
0.0250.028957.2
0.0300.031156.4
0.0350.031855.6
0.0400.038555.0
Table 3. Stabilization control gains at each steady-state working point.
Table 3. Stabilization control gains at each steady-state working point.
P S MPa A J mm2 A Z mm2 K Z Phase/°
3107.14131780.030569.2
3021.2589050.025187.6
8058.7209180.016189.4
160138.753300.029983.1
240359.532860.030589.8
5104.96916550.030261.5
3014.7505470.028884.5
8039.9049310.015184.5
16085.5461940.020176.0
240148.389490.025970.4
7104.03608390.029858.2
3011.9699120.030182.6
8032.1752690.015281.1
16067.2138190.014173.7
240109.703000.015968.5
9103.48626890.031156.4
3010.3347030.028181.3
8027.6932580.022176.6
16057.1608840.012170.7
24090.9952310.014962.8
Table 4. Summary table of the designed initial parameters.
Table 4. Summary table of the designed initial parameters.
Parameter/UnitValue
x s z d , 0 /mm10
x s c d , 0 /mm23.39155
x u z , 0 /mm0.1
x u i n , 0 /mm0
x u o u t , 0 /mm0
Table 5. Geometry relationship value pair of the adjusting orifice.
Table 5. Geometry relationship value pair of the adjusting orifice.
x u z mm A Z mm2
00
0.13.4862689
0.1176789394.0360839
0.1489904024.9691655
0.2209159757.1413178
0.32561712910.334703
0.3838096211.969912
0.4761895214.750547
0.70217417321.258905
0.95852289827.693258
1.16132882632.175269
1.66985922139.904931
2.81263756557.160884
2.94156599558.720918
3.46907537467.213819
4.76924381385.546194
5.04034018190.995231
6.295895148109.70300
8.122958041138.75330
8.445238643148.38949
16.59749231359.53286
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhao, W.; Wang, X.; Long, Y.; Zhou, Z.; Tian, L. Robust LQR Design Method for the Aero-Engine Integral Constant Pressure Drop Control Valve with High Precision. Aerospace 2023, 10, 428. https://doi.org/10.3390/aerospace10050428

AMA Style

Zhao W, Wang X, Long Y, Zhou Z, Tian L. Robust LQR Design Method for the Aero-Engine Integral Constant Pressure Drop Control Valve with High Precision. Aerospace. 2023; 10(5):428. https://doi.org/10.3390/aerospace10050428

Chicago/Turabian Style

Zhao, Wenshuai, Xi Wang, Yifu Long, Zhenhua Zhou, and Linhang Tian. 2023. "Robust LQR Design Method for the Aero-Engine Integral Constant Pressure Drop Control Valve with High Precision" Aerospace 10, no. 5: 428. https://doi.org/10.3390/aerospace10050428

APA Style

Zhao, W., Wang, X., Long, Y., Zhou, Z., & Tian, L. (2023). Robust LQR Design Method for the Aero-Engine Integral Constant Pressure Drop Control Valve with High Precision. Aerospace, 10(5), 428. https://doi.org/10.3390/aerospace10050428

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop