Aircraft Conﬁguration Parameter Boundaries Based on Closed-Loop Flying Qualities Requirements

: For aircraft employing the ﬂy-by-wire technique, the closed-loop dynamic characteristics are determined by both the conﬁguration design and the ﬂight control system. As the capacity of the control system has certain limitations, the conﬁguration parameters are also constrained by the requirements of the closed-loop ﬂying qualities. This paper presents an aircraft conﬁguration parameter boundaries determination method based on closed-loop ﬂying qualities requirements independent of the actual ﬂight control law design, mainly aiming at the parameters that affect the stability and control characteristics. First, a nonlinear dynamic inversion-based ﬂight control law is adopted to decouple the control law gains from the conﬁguration parameters and to study the relationship between the conﬁguration parameters and closed-loop ﬂying qualities. Second, a ﬂying qualities evaluation scheme is established by selecting the most severe ﬂight conditions and the evaluation criteria that are most sensitive to changes in the parameters. Finally, the parameter boundaries according to the requirements of Level 1 ﬂying qualities are determined by searching for the critical values that lead to degradation of the ﬂying qualities. The proposed method is veriﬁed by an application example of the design ranges of a sample aircraft’s wing position, horizontal tail area, center of gravity, vertical tail vertical tail position.


Introduction
Aircraft configuration design depends on the design requirements of different types of parameters [1,2]. For example, the determination of wing area, fuselage geometry and cross section mainly considers the payload requirement; the airfoil, aspect ratio and sweep angle of the wing are designed according to the aerodynamic characteristics requirements; wing position, horizontal tail area, center of gravity (CG), vertical tail area and position have considerable influences on the dynamic characteristics of aircraft, and the design of these parameters is conducted mainly in consideration of stability and control requirements [3].
The flying qualities evaluation can comprehensively assess the stability and control characteristics, the performance of the flight control system, and the adaptability to the flight missions [4][5][6][7][8]. Thus, the configuration design should take into account the flying qualities requirements. As the developments of digital flight control systems and the fly-by-wire technique allow more advanced flight control laws (FCL) to be implemented on aircraft [9][10][11][12][13][14][15], the configuration design has become more diversified. However, the capacity of FCL has certain limitations, and flying qualities cannot be guaranteed only relying on the control system. Considering that the flying qualities of the closed-loop aircraft are determined by both the configuration design and the flight control system, the configuration parameters are also constrained by the closed-loop flying qualities requirements [16]. Therefore, it is necessary to establish the relationship between the configuration parameters and the closed-loop flying qualities requirements in the conceptual aircraft design.

•
Preliminary Sizing: definition of requirements, selection of an aircraft configuration, selection of a propulsion system, etc. • Layout sizing: fuselage sizing, wing sizing, empennage sizing, design for high lift, etc. • Specific design: stability and control design, flight control system design, landing gear design, structures design, etc. • Aircraft performance analysis: aerodynamic analysis, polar estimation, mass estimation, flying qualities evaluation, mission analysis, etc.
The whole design process has an iterative nature to realize optimization design [27]. Stability and control design and flying qualities evaluation are often in later stages of the design phase. Configuration parameters adjustment due to stability and control requirements or flying qualities requirements may lead to much iterative work. This paper proposes a method to determine the configuration parameter boundaries based on closed- loop flying qualities requirements to simplify the stability and control design. The proposed parameter boundaries aim to transform the closed-loop flying qualities requirements into the design ranges of target configuration parameters. In this way, the flying qualities requirements can be considered at the layout sizing stage during optimization design.
The determination process of the parameter boundaries is as follows. The first step is to construct a high-order closed-loop aircraft model according to the basic configuration design; the model can be divided into an aircraft motion model and the FCL based on the model reference NDI method. The second step is to establish the flying qualities evaluation scheme by selecting the most stringent flight conditions and the evaluation criteria that are most sensitive to the parameter changes. The third step is to establish a reasonable variation series of the target parameters and to search for the critical values that lead to degradation of the flying qualities. Finally, the parameter boundaries according to the requirements of the Level 1 flying qualities (Level 1 FQ) are determined. The method for determining parameter boundaries is shown in Figure 1.
The whole design process has an iterative nature to realize optimization design [27]. Stability and control design and flying qualities evaluation are often in later stages of the design phase. Configuration parameters adjustment due to stability and control requirements or flying qualities requirements may lead to much iterative work. This paper proposes a method to determine the configuration parameter boundaries based on closed-loop flying qualities requirements to simplify the stability and control design. The proposed parameter boundaries aim to transform the closed-loop flying qualities requirements into the design ranges of target configuration parameters. In this way, the flying qualities requirements can be considered at the layout sizing stage during optimization design.
The determination process of the parameter boundaries is as follows. The first step is to construct a high-order closed-loop aircraft model according to the basic configuration design; the model can be divided into an aircraft motion model and the FCL based on the model reference NDI method. The second step is to establish the flying qualities evaluation scheme by selecting the most stringent flight conditions and the evaluation criteria that are most sensitive to the parameter changes. The third step is to establish a reasonable variation series of the target parameters and to search for the critical values that lead to degradation of the flying qualities. Finally, the parameter boundaries according to the requirements of the Level 1 flying qualities (Level 1 FQ) are determined. The method for determining parameter boundaries is shown in Figure 1. It should be emphasized that the parameter boundaries proposed in this paper correspond to specific flying qualities criteria, so they only reflect a part of stability and control requirements. Meanwhile, these configuration parameters should also meet other requirements. For example, the CG position design should meet the trim requirements, the vertical tail design should meet the controllability requirements in the single-engine out condition at take-off, and so forth. Designers should seek to reach the design configuration that satisfactorily meets all requirements in the optimization design process. It should be emphasized that the parameter boundaries proposed in this paper correspond to specific flying qualities criteria, so they only reflect a part of stability and control requirements. Meanwhile, these configuration parameters should also meet other requirements. For example, the CG position design should meet the trim requirements, the vertical tail design should meet the controllability requirements in the single-engine out condition at take-off, and so forth. Designers should seek to reach the design configuration that satisfactorily meets all requirements in the optimization design process.

Wing Position and Horizontal Tail Area
The longitudinal stability and control characteristics of aircraft are determined mainly by the stability and damping derivatives of the pitch moment, and these derivatives are affected by the wing position, CG position, horizontal tail position, horizontal tail area, and other parameters. The horizontal tail is usually installed at the rear end of the fuselage to ensure the control power of the elevator, and the variation range of the tail position is small. Thus, the longitudinal characteristics are governed predominantly by the wing position and horizontal tail area [1]. These two parameters are usually designed in collaboration according to the requirements of the longitudinal flying qualities.
The stability derivative of the pitch moment C mα can be expressed as [28]: where x ac.WB is the mean aerodynamic center of the wing-body configuration and is affected mainly by the wing position; C Lα.WB is the lift coefficient derivative of the wingbody configuration; x c.g. is the CG position of the aircraft; C Lα.T is the lift coefficient derivative of the horizontal tail; V T is the horizontal tail volume ratio, affected mainly by the horizontal tail position and its area; and C mα.p is the pitch moment derivative of the propulsion system. V T is defined as: where l T is the distance between the aircraft CG and the aerodynamic center of the horizontal tail, S T is the horizontal tail area, c is the mean aerodynamic chord of the aircraft wing, and S is the wing area. Consequently, forward movement of the wing position and a decrease in the horizontal tail area lead to a decrease in |C mα | and a reduction in the static longitudinal stability. Similarly, rearward movement of the wing position and an increase in the horizontal tail area lead to an increase in |C mα | and an increase in the demand on the pitch control moment. The damping derivatives of the pitch moment, namely, C mq and C m . α , are generated mainly by the horizontal tail, whereas the influence of the wing-body configuration accounts for only about 10% of the overall influence [28]. The contributions of the tail to C mq and C m . α are [28]: where ∂ε/∂α is the derivative of the downwash angle of the horizontal tail with respect to the angle of attack. A decrease in horizontal tail area leads to a decrease in both C mq and C m . α . The stability and damping derivatives of the pitch moment mainly affect the longitudinal short-period mode characteristics of aircraft. The natural frequency ω n.sp and damping ratio ξ sp of the short-period mode can be expressed as: α , and Z α are the dimensional dynamic derivatives related to C mα , C mq , C m . α and C Lα , respectively. The detailed derivation process from (3) to (4) is described in the references [29,30].
Combining Equations (3) and (4) demonstrates that a decrease in the horizontal tail area leads to a decrease in ξ sp , slows down the short-period convergence, and decreases the dynamic stability of the aircraft; in contrast, an increase in the tail area leads to an increase in the damping ratio and restrains the pitch response to the pilot control input.
In a general conceptual design process, the horizontal tail is first of all sized from static conditions as done with the "V-diagram" [27,31,32], which corresponds to the basis of the longitudinal stability and control requirements. Then, the horizontal tail design should consider the dynamic stability and control requirements, which are the primary focus of this paper. According to the above analysis, the longitudinal dynamic stability requirements of the aircraft should be considered to determine the forward range of the wing position and the maximum limit of the horizontal tail area. Furthermore, the rearward range of the wing position and the minimum limit of the horizontal tail area should be determined by the longitudinal maneuverability requirements.

Center of Gravity Position
The CG is not a fixed point but moves with the initial loading condition and fuel consumption of the aircraft. For an aircraft with determined aerodynamic shape, a change in the CG position x cg causes C mα to vary and thus affects the longitudinal characteristics. Therefore, the range of the CG location needs to be strictly designed according to the requirements of the longitudinal flying qualities.
A forward CG represents a large |C mα | and requires a larger pitch control moment, while a rearward CG position represents a small |C mα | and poor longitudinal stability. Therefore, in this paper, the forward CG limit mainly considers the longitudinal maneuverability requirements, and the rearward CG limit is determined by the longitudinal stability requirements.
Additionally, the forward CG limit should allow the aircraft to be trimmed at maximum lift coefficient in ground effect. The final forward CG limit will be determined by both the longitudinal maneuverability and trim requirements, which leads to future investigations.

Vertical Tail Area and Position
Directional stability, also known as weathercock stability, is provided primarily by the vertical tail. Meanwhile, a large vertical tail is accessible spiral stability. However, a large vertical tail surface may produce an excessive aerodynamic drag force during high-speed flight at high altitudes in cruise, and such a large surface would require a higher structural strength and mass of the airframe. Therefore, it is necessary to restrict the vertical tail area to optimize both the aerodynamic characteristics in cruise and the airframe structure. To ensure the directional stability after the vertical tail area is reduced, it is necessary to increase the vertical tail arm by adjusting the position of the vertical tail. Similar to the horizontal tail, the vertical tail is also installed at the end of the fuselage to ensure the control power of the rudder. Thus, the vertical tail position should be adjusted by changing the chord length, sweep angle and other parameters without moving the root position. It should be noted that a large sweep angle may lead to the decrease of the vertical tail lift coefficient. Thus, the vertical tail design should avoid a large sweep angle [29,30].
The contributions of the vertical tail to the stability and damping derivatives of the directional moment C nβ.V and C nr.V can be expressed as [28]: where C Yβ.V is the lift coefficient derivative of the vertical tail; V V is the vertical tail volume ratio; S V is the vertical tail area; l V is the force arm of the vertical tail; and b is the wingspan of the aircraft. The stability and damping derivatives of the directional moment affect mainly the Dutch roll mode characteristics of the aircraft. The natural frequency ω n.dr and damping ratio ξ dr of the Dutch roll mode can be expressed as: where N β , N r , Y β , and Y r are the dimensional dynamic derivatives related to C nβ , C nr , C Yβ and C Yr , respectively. The detailed derivation process from (5) to (6) is also described in the references [29,30]. Combining Equations (5) and (6) reveals that a reduced vertical tail arm and a decrease in the tail area cause a reduction in both C nβ and ξ dr and decreases the directional static and dynamic stability. In a general conceptual design process, vertical tail sizing for static stability has to compensate for the destabilizing effect of the fuselage; vertical tail sizing for static control has to compensate the yaw moment from a one-engine out flight case [29,30]; vertical tail sizing for roll response control requires the proper balance of the rolling time constant and the roll sensitivity; there are also spiral stability and directional dynamic stability requirements. Directional dynamic stability mainly refers to the Dutch roll mode stability, which is the primary focus of this paper. The eventual vertical tail area will be the minimum area that meets all these requirements.

Aircraft Modelling
The aircraft motion model is based on the six-degree-of-freedom equations of motion.
To determine the configuration parameter boundaries based on closed-loop flying qualities requirements, the motion model should be constantly updated according to the changes in the configuration parameters.
To easily adapt the control laws to continuously change the aircraft configuration during the design process, NDI controllers offer considerable advantages. For instance, the model reference NDI control law [33,34] uses transfer function-based reference models to describe the desired aircraft dynamics. The reference models can be expressed as [33]: where p c , q c , and r c are the rate commands of the pilot; p r , q r , and r r are the corresponding reference rate commands generated by the reference models; K q.r and K p.r are the control gains; ξ sp.r , ω sp.r , and T θ2.r are the desired damping ratio, natural frequency and time constant of the short-period mode, respectively; and ω p.r and ω r.r are the desired roll and directional frequencies. The reference commands are achieved by the actual dynamic inversion at the core of the control law, which computes the necessary surface positions. A proportional-integral (PI) compensator is added to drive down the error between the desired dynamics and the actual dynamics. Filters are also needed to attenuate undesirable noise from the feedback sensors. The structure of the high-order closed-loop aircraft model with the model reference NDI control law is shown in Figure 2.
The parameters of the model reference NDI control law are designed according to the desired dynamic characteristics and do not need to be adjusted according to changes in the configuration parameters [25]. As a result, the FCL achieves the desired closed-loop flying qualities of the varying aircraft configuration and eliminates the disturbance attributable to variation in the control system. Therefore, the model reference NDI method is employed to design the control law in this paper, which is convenient for studying the relationship between the configuration parameters and the closed-loop flying qualities.  It should be emphasized that the design ranges of the configuration parameters derived from the NDI control law represent the capability of control laws that makes the closed-loop aircraft meet the flight quality requirements. For aircraft adopting other control laws, such as SCAS, the conclusions are still valid.

Selection of Flight Conditions for Evaluation
Any variation in the flight conditions will lead to a change in the flying qualities evaluation results, so the requirements for the configuration parameters differ under different flight conditions. Consequently, to determine configuration parameter boundaries based on closed-loop flying qualities requirements, it is necessary to determine the flight conditions with the most severe requirements for the configuration parameters [1,3].
The previous discussion shows that the forward range of the wing position, the minimum limit of the horizontal tail area, and the rearward CG limit are the main factors influencing the longitudinal dynamic stability requirements of the aircraft. In contrast, the longitudinal maneuverability requirements are governed primarily by the rearward range of the wing position, the maximum limit of the horizontal tail area and the forward CG limit. For aircraft whose load and fuel are arranged mainly in the front region of the airframe, the forward CG limit appears during the take-off phase [1]. In this case, the stability margin of the aircraft is large during the take-off phase, and the longitudinal maneuverability is poor. Considering the low flight speed during the take-off phase, the aerodynamic efficiency of the control surface is low, which worsens the maneuverability. Therefore, the rearward range of the wing position, the maximum limit of the horizontal tail area and the forward CG limit should be determined in the take-off phase. The flight speed, flight altitude, weight distribution, flaps, landing gear and other settings should be consistent with the most forward CG position. With loading and fuel consumption, the CG position reaches the rearward limit in landing phase, during which the stability margin of the aircraft decreases and the most stringent requirements for the longitudinal stability are generated. Therefore, the forward range of the wing position, the minimum limit of the horizontal tail area and the rearward CG limit should be determined in landing phase, and the exact flight condition should be set to the state with the most rearward CG position. The selected flight conditions for evaluating an aircraft whose load and fuel are arranged mainly at the rear of the airframe are the opposite. For an aircraft whose load and fuel cannot be judged clearly whether they are arranged in the front or rear position (for example, the fuel is in the wings and the wingbox), the flight conditions for longitudinal flying qualities evaluation can be selected directly according to the change of the CG position during flight.
The vertical tail design for the directional dynamic stability should be determined under the flight condition with the most stringent requirements. Due to the decrease in damping coefficients during high-speed and high-altitude flight, the Dutch roll mode of the aircraft will deteriorate, thereby increasing the dynamic stability requirements. Therefore, the cruise phase at high speed and high altitude is selected to determine the

Selection of the Flying Qualities Evaluation Criteria
Current flying qualities standards provide multiple evaluation criteria [35][36][37][38]. Each evaluation criterion has its own inspection focus, and thus, the same configuration design may yield different evaluation results according to different criteria. To determine configuration parameter boundaries based on closed-loop flying qualities requirements, it is necessary to determine the evaluation criteria that are most sensitive to changes in the configuration parameters to obtain the strictest parameter design ranges.

CAP Criterion
The control anticipation parameter (CAP) is defined as the ratio of the initial pitch acceleration to the steady-state normal acceleration [38]. The CAP value can reflect the short-period mode flying qualities of an aircraft, including the longitudinal stability and maneuverability. To apply the CAP criterion to a high-order closed-loop aircraft, the loworder equivalent system (LOES) theorem is adopted to match the characteristics of the high-order system. The longitudinal LOES of the aircraft can be expressed as: where α is the aircraft angle of attack; q is the pitch rate; F e is the pitch control input; 1/T θ 2 is the zero of the pitch attitude transfer function; K α and K q are the equivalent system gains of the transfer functions; and τ α and τ q are the equivalent time delays.
The CAP of the LOES is defined as: where V is the flight speed and g is the acceleration due to gravity. The Level 1 FQ of the CAP criterion are defined by the short-period damping ratio ζ sp of the equivalent system and the CAP value. The upper boundary of the CAP value and the lower boundary of ζ sp reflect the stability requirements of an aircraft. A large CAP value and small ζ sp indicate a lack of longitudinal stability, and thus, a pilot's small control input may lead to a violent pitch response. Therefore, the upper boundary of the CAP value and the lower boundary of ζ sp can be used to determine the forward range of the wing position, the minimum limit of the horizontal tail area and the rearward CG limit.
In contrast, the lower boundary of the CAP value and the upper boundary of ζ sp reflect the maneuverability requirements. A small CAP value and large ζ sp represent a large longitudinal stability margin, which may increase the pilot's workload to achieve pitch motion. Thus, the rearward range of the wing position, the maximum limit of the horizontal tail area and the forward CG limit should be designed according to the lower boundary of the CAP value and the upper boundary of ζ sp .

Chalk Criterion
The Chalk criterion is a time-domain criterion that focuses on the characteristics of the pitch rate response and defines three parameters for evaluation: the effective time delay t 1 , the transient peak ratio ∆q max /∆q min and the effective rise time ∆t, as shown in Figure 3 [38].

Chalk Criterion
The Chalk criterion is a time-domain criterion that focuses on the characteristics of the pitch rate response and defines three parameters for evaluation: the effective time delay 1 t , the transient peak ratio max min / q q   and the effective rise time t  , as shown in Figure 3 [38].  In Figure 3, q ss is the steady-state pitch rate; q max is the maximum pitch rate; q min is the first minimum pitch rate; the instant t 1 is defined by the intersection between the tangent at the maximum slope and the time axis; the instant t 2 is defined by the intersection between the tangent at the maximum slope and the horizontal line of the steady-state pitch rate; and ∆t is the difference between t 1 and t 2 . t 1 reflects the time delay between the real pitch rate response and an ideal pitch rate response with the maximum pitch acceleration, ∆q max /∆q min reflects the short-period damping ratio of the pitch response, and ∆t reflects the rise time to accelerate to the steady-state pitch rate of the pitch rate ideal response.
The Chalk criterion is based on the time-domain pitch response of the aircraft without identifying the frequency-domain characteristics. Thus, the Chalk criterion can be used to avoid the situation in which a frequency-domain criterion is not applicable to aircraft with an unconventional response. The design ranges of the configuration parameters should be determined according to the more stringent requirements of the CAP and Chalk criteria to enhance the universality of the design method based on closed-loop flying qualities requirements.

Dutch Roll Criterion
The flying qualities with respect to the yaw axis can be evaluated by the characteristics of the lateral-directional oscillatory (Dutch roll) response to the yaw controller. The Dutch roll frequency and damping ratio can be defined by matching the higher-order sideslip response to the yaw control input to the following lower-order form [38]: where β is the sideslip angle; F r is the yaw control input; K β is the equivalent system gain of the transfer function; and τ eβ is the equivalent time delay of the transfer function. The flying qualities requirements of the Dutch roll criterion are specified in terms of the minimum values of ω n.dr , ξ dr and ξ dr ω n.dr , and the evaluation results can reflect the lateral-directional stability of the aircraft. This criterion can be used to determine the vertical tail area and position boundaries.

Application of the Proposed Parameter Boundaries Design Method
This section takes a medium-range transport aircraft as an application example of the configuration parameter boundaries design method. The definitions of the target configuration parameters of the sample aircraft are illustrated in

Application of the Proposed Parameter Boundaries Design Method
This section takes a medium-range transport aircraft as an application example of the configuration parameter boundaries design method. The definitions of the target configuration parameters of the sample aircraft are illustrated in  The configuration parameters of the basic aircraft design are depicted in Table 1. The aircraft motion model and the inner loop of the NDI control law are based on the basic aircraft design, and the reference model parameters of the NDI control law are designed according to the flying qualities requirements, as shown in Table 2. By combining the motion model with the FCL, the high-order closed-loop aircraft model is established.

Boundaries of the Wing Position and Horizontal Tail Area
For the sample aircraft, the forward CG limit appears in the take-off phase (sea-level, low-speed, forward CG position, maximum take-off weight, maximum thrust, flaps and landing gear retracted), while the rearward CG limit appears in the landing phase (sealevel, low-speed, rearward CG position, minimum landing weight, landing engine setting, flaps and landing gear down). According to the previous discussion, the forward range of x W and the minimum limit of S T should be designed by evaluating the flying qualities in landing phase according to the most stringent requirement among the upper boundary of the CAP value, the lower boundary of ξ sp , and the Chalk criterion. Likewise, the rearward range of x W and the maximum limit of S T should be designed in take-off phase according to the most stringent requirement among the lower boundary of the CAP value, the upper boundary of ξ sp and the Chalk criterion.
Preliminary sizing requires that the CG position is on the MAC of the wing; thus, the first estimate of x W ranges from 10 m to 16 m. Then, according to the requirements of the longitudinal Level 1 FQ of closed-loop aircraft, the design range of S T corresponding to different x W is determined. Since the rearward range of x W and the maximum limit of S T both reflect the maneuverability requirements, the upper boundary of S T corresponding to different x W can also be regarded as the rearward boundary x W corresponding to different S T . Similarly, the lower boundary of S T corresponding to different x W can also be regarded as the forward boundary of x W corresponding to different S T . The design range of x W is determined by the design range of S T .

•
Maximum limit of the horizontal tail area The variation series of x W is established as follows: within the preliminary range from 10 m to 16 m, one point is selected every 1 m, yielding seven points in total. For each X W , a series of S T with an interval of 2 m 2 is selected. Then, the CAP criterion is applied to evaluate the flying qualities of aircraft with different x W -S T combinations in the take-off phase. Partial evaluation results with x W values of 10 m, 12 m, 14 m, and 16 m are shown in Figure 5. The CAP criterion characteristics with x W values of 12 m and 14 m are shown in Table 3.  Table 3.    For a certain x W , the CAP value decreases gradually with increasing S T , and the damping ratio ξ sp increases. The maximum limit of S T defined by the CAP criterion is the critical value that leads to the degradation of the flying qualities to Level 2. When x W is set to 12 m and 14 m, the corresponding maximum limit of S T is 36 m 2 and 32 m 2 , respectively. Hence, as the wing position moves rearward, the maximum limit of S T gradually decreases.
According to all the above evaluation results, the level boundaries of the x W -S T combination defined by the CAP criterion in take-off phase are presented in Figure 6.
According to this figure, to ensure the Level 1 FQ of the CAP criterion, when x W is set as 10 m, the corresponding design range of S T is from 24 m 2 to 38 m 2 ; when x W is set as 16 m, the S T design range is from 14 m 2 to 26 m 2 . The right boundary of Level 1 in Figure 6 determines the maximum limit of S T corresponding to different x W , and the rearward range of x W corresponding to different S T . According to all the above evaluation results, the level boundaries of the W x -T S combination defined by the CAP criterion in take-off phase are presented in Figure 6. According to this figure, to ensure the Level 1 FQ of the CAP criterion, when W x is set as 10 m, the corresponding design range of T S is from 24 m 2 to 38 m 2 ; when W x is set as 16 m, the T S design range is from 14 m 2 to 26 m 2 . The right boundary of Level 1 in Figure 6 determines the maximum limit of T S corresponding to different W x , and the rearward range of W x corresponding to different T S . The Chalk criterion is also applied to evaluate the flying qualities in the take-off phase. Figure 7 shows the influence of changing T S on the pitch rate response with W x = 12 m. The variations in the tail area exert the main influences on the response peak and response speed of the pitch rate. As T S increases, the response peak decreases, and the response time delay increases. The influence of the rearward displacement of the wing is similar to that of changing T S , which means that decreasing W x will reduce both the response peak and the response speed.  Table 4.
The transient peak ratio The Chalk criterion is also applied to evaluate the flying qualities in the take-off phase. Figure 7 shows the influence of changing S T on the pitch rate response with x W = 12 m. The variations in the tail area exert the main influences on the response peak and response speed of the pitch rate. As S T increases, the response peak decreases, and the response time delay increases. The influence of the rearward displacement of the wing is similar to that of changing S T , which means that decreasing x W will reduce both the response peak and the response speed.
According to all the above evaluation results, the level boundaries of the W x -T S combination defined by the CAP criterion in take-off phase are presented in Figure 6. According to this figure, to ensure the Level 1 FQ of the CAP criterion, when W x is set as 10 m, the corresponding design range of T S is from 24 m 2 to 38 m 2 ; when W x is set as 16 m, the T S design range is from 14 m 2 to 26 m 2 . The right boundary of Level 1 in Figure 6 determines the maximum limit of T S corresponding to different W x , and the rearward range of W x corresponding to different T S . The Chalk criterion is also applied to evaluate the flying qualities in the take-off phase. Figure 7 shows the influence of changing T S on the pitch rate response with W x = 12 m. The variations in the tail area exert the main influences on the response peak and response speed of the pitch rate. As T S increases, the response peak decreases, and the response time delay increases. The influence of the rearward displacement of the wing is similar to that of changing T S , which means that decreasing W x will reduce both the response peak and the response speed.  Table 4.
The transient peak ratio Among all Chalk criterion parameters, the effective time delay t 1 is the most affected by the variations in and S T . The evaluation results of t 1 are shown in Table 4. The transient peak ratio ∆q max /∆q min and the effective rise time ∆t are also affected, but the FQ Level is not degraded. With the rearward movement of the wing and an increase in the tail area, the t 1 of the pitch response increases, which indicates poor maneuverability. The FQ Level is degraded to Level 2 (0.12 s < t 1 < 0.17 s) or even Level 3 (t 1 > 0.17 s). In the chosen ranges of x W and S T , an increase in x W and a decrease in S T will not degrade the flying qualities of the Chalk criterion. The level boundaries of the x W -S T combination defined by the Chalk criterion in the take-off phase are shown in Figure 8. According to this figure, to ensure the Level 1 FQ of the Chalk criterion, the maximum limit of S T is 38 m 2 when x W = 10 m, and the maximum limit is 26 m 2 when x W = 16 m. Compared with the parameter boundaries defined by the CAP criterion shown in Figure 6, the Chalk criterion defines a more relaxed Level 1 boundary of the rearward range of x W and the maximum limit of S T . In contrast, the Level 2 boundary defined by the Chalk criterion is stricter and can be used as a supplement to the CAP criterion boundaries. > 0.17 s). In the chosen ranges of W x and T S , an increase in W x and a decrease in T S will not degrade the flying qualities of the Chalk criterion. The level boundaries of the W x -T S combination defined by the Chalk criterion in the take-off phase are shown in Figure 8. According to this figure, to ensure the Level 1 FQ of the Chalk criterion, the maximum limit of T S is 38 m 2 when W x = 10 m, and the maximum limit is 26 m 2 when W x = 16 m. Compared with the parameter boundaries defined by the CAP criterion shown in Figure 6, the Chalk criterion defines a more relaxed Level 1 boundary of the rearward range of W x and the maximum limit of T S . In contrast, the Level 2 boundary defined by the Chalk criterion is stricter and can be used as a supplement to the CAP criterion boundaries. • Minimum limit of the horizontal tail area The selection of the x W and S T series is the same as the design of the maximum S T range. The CAP criterion is applied to evaluate the flying qualities in the landing phase. Partial evaluation results with x W values of 10 m, 12 m, 14 m, and 16 m are shown in Figure 9, where the corresponding minimum limits of S T are 26 m 2 , 22 m 2 , 20 m 2 , and 16 m 2 , respectively. As the wing position moves forward, the minimum limit of S T gradually increases.
The level boundaries of the x W -S T combination defined by the CAP criterion in landing phase are shown in Figure 10. According to this figure, to ensure the Level 1 FQ of the CAP criterion, the design range of S T is from 26 m 2 to 40 m 2 when x W = 10 m and from 16 m 2 to 32 m 2 when x W = 16 m. The left boundary of Level 1 in Figure 10 can determine the minimum limit of S T corresponding to different x W , and the forward range of x W corresponding to different S T .   According to the evaluation results of the Chalk criterion in the landing phase, the effective rise time and the response time delay decrease as the wing moves forward and the tail area decreases, which represents a better response speed of the pitch rate. Meanwhile, the pitch damping characteristics become worse as the transient peak ratio   According to the evaluation results of the Chalk criterion in the landing phase, the effective rise time and the response time delay decrease as the wing moves forward and the tail area decreases, which represents a better response speed of the pitch rate. Meanwhile, the pitch damping characteristics become worse as the transient peak ratio According to the evaluation results of the Chalk criterion in the landing phase, the effective rise time and the response time delay decrease as the wing moves forward and the tail area decreases, which represents a better response speed of the pitch rate. Meanwhile, the pitch damping characteristics become worse as the transient peak ratio ∆q max /∆q min increases, but the FQ Level is not degraded. Therefore, the forward boundary of x W and the minimum limit of S T should be determined only by the CAP criterion.

•
Level boundaries of the wing position and horizontal tail area To determine the design ranges of x W and S T according to the Level 1 FQ requirements in the entire flight envelope, the most stringent boundaries in Figures 6, 8 and 10 are combined to form the final parameter level boundary diagram, as shown in Figure 11.
According to Figure 11, the Level 1 boundary of the rearward x W and the maximum S T (right side in Figure 11) is determined by the evaluation results in the take-off phase; the upper half is defined by the Chalk criterion and the lower half is defined by the CAP criterion. The Level 2 boundary (top right in Figure 11) is determined in take-off phase by the Chalk criterion. The Level 1 boundary of the forward x W and the minimum S T (left side in Figure 11) is defined in the landing phase by the CAP criterion. The design range of S T determined by the Level 1 FQ requirements is from 26 m 2 to 38 m 2 when x W = 10 m and from 16 m 2 to 26 m 2 when x W = 16 m. In the configuration design process, the design range of S T can be defined as the horizontal interval within the Level 1 zone (shaded part in Figure 11) corresponding to the designed x W . Similarly, the design range of x W is defined as the vertical interval within the Level 1 zone according to the designed S T .
To determine the design ranges of W x and T S according to the Level 1 FQ requirements in the entire flight envelope, the most stringent boundaries in Figures 6, 8 and 10 are combined to form the final parameter level boundary diagram, as shown in Figure  11.

Boundaries of the Center of Gravity Position
The The series for the forward cg x limit design is set as follows: from 6% to 36%, a point is selected every 3% of the wing MAC, yielding ten points in total. Likewise, the series for the rearward cg x limit design is set as follows: from 38% to 65%, a point is selected every 3% of the wing MAC, also yielding ten points in total. Then, the CAP criterion is applied to evaluate the flying qualities of the aircraft with the forward cg x in the take-off phase, while the rearward cg x designs are evaluated in the landing phase. The evaluation results are presented in Figure 12. Compared with the basic design of the sample aircraft in Figure 11. Complete level boundaries of x W and S T .

Boundaries of the Center of Gravity Position
The forward limit of x cg should be designed in the take-off phase according to the most stringent requirement among the lower boundary of the CAP value, the upper boundary of ξ sp , and the Chalk criterion. Likewise, the rearward limit of x cg should be designed by evaluating the flying qualities in the landing phase according to the most stringent requirement among the upper boundary of the CAP value, the lower boundary of ξ sp , and the Chalk criterion.
The series for the forward x cg limit design is set as follows: from 6% to 36%, a point is selected every 3% of the wing MAC, yielding ten points in total. Likewise, the series for the rearward x cg limit design is set as follows: from 38% to 65%, a point is selected every 3% of the wing MAC, also yielding ten points in total. Then, the CAP criterion is applied to evaluate the flying qualities of the aircraft with the forward x cg in the take-off phase, while the rearward x cg designs are evaluated in the landing phase. The evaluation results are presented in Figure 12. Compared with the basic design of the sample aircraft in the take-off phase, as the CG position moves forward (TO phase, x cg -), the CAP value decreases gradually, the damping ratio ξ sp increases, and the FQ Level is degraded from Level 1 to Level 2. In the landing phase, as the CG position moves rearward (LD phase, x cg +), the CAP value increases, and the damping ratio ξ sp decreases.
The forward x cg designs are also evaluated by the Chalk criterion in the take-off phase. The effective time delay t 1 is the parameter most affected by the variation in x cg ; the evaluation results of t 1 are shown in Table 5. As CG moves forward, the t 1 value of the pitch response increases, which indicates poor maneuverability. The FQ Level is degraded to Level 2 or even Level 3. According to the evaluation results of the Chalk criterion in the landing phase, the pitch damping characteristics of the aircraft become worse as the CG moves rearward, but the FQ Level is not degraded. Therefore, the Chalk criterion is only applied to design the forward x cg limit.
The level boundaries of x cg are defined as in Figure 13. According to this figure, the Level 1 boundary of the forward CG limit is 21% of the wing MAC, determined by the upper boundary of ξ sp in the CAP criterion in the take-off phase. The Level 2 boundary of the forward CG limit is 12% of the wing MAC, determined by the Chalk criterion. Finally, the Level 1 boundary of the rearward CG limit is 59% of the wing MAC, determined by the upper boundary of the CAP value in the landing phase. The x cg range should be kept within the Level 1 zone (shaded part in Figure 13).  The forward cg x designs are also evaluated by the Chalk criterion in the take-off phase. The effective time delay 1 t is the parameter most affected by the variation in cg x ; the evaluation results of 1 t are shown in Table 5. As CG moves forward, the 1 t value of the pitch response increases, which indicates poor maneuverability. The FQ Level is degraded to Level 2 or even Level 3. According to the evaluation results of the Chalk criterion in the landing phase, the pitch damping characteristics of the aircraft become worse as the CG moves rearward, but the FQ Level is not degraded. Therefore, the Chalk criterion is only applied to design the forward cg x limit.
The level boundaries of cg x are defined as in Figure 13. According to this figure, the Level 1 boundary of the forward CG limit is 21% of the wing MAC, determined by the upper boundary of sp  in the CAP criterion in the take-off phase. The Level 2 boundary of the forward CG limit is 12% of the wing MAC, determined by the Chalk criterion. Finally, the Level 1 boundary of the rearward CG limit is 59% of the wing MAC, determined by the upper boundary of the CAP value in the landing phase. The cg x range should be kept within the Level 1 zone (shaded part in Figure 13).

Boundaries of the Vertical Tail Area and Position
The

Boundaries of the Vertical Tail Area and Position
The vertical tail area S V and position x V should be designed in the cruise phase (high speed, high altitude, cruise engine setting and configuration) according to the requirements of the Dutch roll criterion. The variation series of S V is set as follows: below the basic value of S V , a point is selected every 1 m 2 , and for each S V , a series of x V with an interval of 1 m is selected. As x V is adjusted by changing the quarter chord point of the vertical tail without moving the root position, the design range of x V is limited from 26 m to 31 m.
Then, the Dutch roll criterion is applied to evaluate the flying qualities of different S V -x V combinations. Partial evaluation results with S V values of 10 m 2 , 12 m 2 , 14 m 2 , and 16 m 2 are shown in Figure 14, and the Dutch roll characteristics with S V values of 12 m 2 and 14 m 2 are shown in Table 6.
(high speed, high altitude, cruise engine setting and configuration) according to the requirements of the Dutch roll criterion. The variation series of V S is set as follows: below the basic value of V S , a point is selected every 1 m 2 , and for each V S , a series of V x with an interval of 1 m is selected. As  Table 6.  For a certain S V , both ω n.dr and ξ dr of the Dutch roll mode decrease gradually with the forward movement of the tail. The forward range of x V is defined as the critical value that leads to the degradation of directional flying qualities. When S V is set as 12 m 2 and 14 m 2 , the corresponding forward range of x V is 30 m and 28 m, respectively. With decreasing S V , the design range of x V gradually decreases. The minimum limit of S V is defined as the critical value without the corresponding x V that meets the requirements of the Level 1 FQ.
According to all evaluation results, the level boundaries of the S V -x V combination are defined as in Figure 15.
Considering that the drag force due to the vertical tail is largely related to the vertical tail area, it is necessary to restrict the tail area while satisfying the flying qualities requirements. According to Figure 15, x V within the design range can ensure Level 1 FQ of the closed-loop aircraft when S V is not less than 15 m 2 . The minimum limit of S V is 11 m 2 and the corresponding x V is 31 m. In the configuration design process, the design range of S V can be defined as the horizontal interval within the shaded Level 1 zone in Figure 15, corresponding to the designed x V . Similarly, the design range of x V is defined as the vertical interval within the Level 1 zone according to the designed S V . m 2 and 14 m 2 , the corresponding forward range of V x is 30 m and 28 m, respectively.
With decreasing V S , the design range of V x gradually decreases. The minimum limit of V S is defined as the critical value without the corresponding V x that meets the requirements of the Level 1 FQ. According to all evaluation results, the level boundaries of the V S -V x combination are defined as in Figure 15. Considering that the drag force due to the vertical tail is largely related to the vertical tail area, it is necessary to restrict the tail area while satisfying the flying qualities requirements. According to Figure 15

Conclusions
In this study, the influences of typical configuration parameters on flying qualities are analyzed, and a configuration parameter boundaries design method based on closed-loop flying qualities requirements is proposed and verified by an application example.

Conclusions
In this study, the influences of typical configuration parameters on flying qualities are analyzed, and a configuration parameter boundaries design method based on closed-loop flying qualities requirements is proposed and verified by an application example.
(1) Five typical configuration parameters that should be designed based on flying qualities requirements are determined: wing position, horizontal tail area, center of gravity position, vertical tail area and position.
(2) The proposed configuration parameter boundaries design method based on closedloop flying qualities requirements is achieved as follows. First, the flight control law is designed based on the model reference nonlinear dynamic inversion method, and a highorder closed-loop aircraft model is established. Second, the most stringent flight conditions and the evaluation criteria that are most sensitive to the target parameters are determined. Third, by searching for the critical values that lead to the degradation of flying qualities, the closed-loop flying qualities requirements are transformed into the design ranges of target configuration parameters.
(3) An application example involving the typical parameter boundaries design of a sample aircraft is provided. The wing position of the sample aircraft ranges from 10 m to 16 m, the design range of the horizontal tail area corresponding to the forward limit of wing position is from 26 m 2 to 38 m 2 , and the horizontal tail area range corresponding to the rear limit of wing position is from 16 m 2 to 26 m 2 . The design ranges of the wing position and horizontal tail area are defined by the level boundary diagram. The forward limit of the center of gravity is 21% of the average aerodynamic chord, while the rearward limit is 59%. The minimum limit of the vertical tail area is 11 m 2 , and the design ranges of the vertical tail area and position are defined by the level boundary diagram.

Data Availability Statement:
The high-order closed-loop aircraft model ( Figure 2) and data presented in this study can be obtained from the corresponding author on request.