Development of Unified High-Fidelity Flight Dynamic Modeling Technique for Unmanned Compound Aircraft

This study presents the unified high-fidelity flight dynamic modeling technique for compound aircraft. The existing flight dynamic modeling technique is absolutely depended on the experimental data measured by wind tunnel. It means that the existing flight dynamic model cannot be used for analyzing a new configuration aircraft. The flight dynamic modeling has to be implemented when a performance analysis has to be performed for new type aircraft. This technique is not effective for analyzing the performance of the new configuration aircraft because the shapes of compound aircraft are very various. The unified highfidelity flight dynamic modeling technique is developed in this study to overcome the limitation of the existing modeling technique. First, the unified rotor and wing models are developed to calculate the aerodynamic forces generated by rotors and wings. The revolutions per minute (RPM) and pitch change with rotation direction are addressed by rotor models. The unified wing model calculates the induced velocity by using the vortex lattice method (VLM) and the Biot–Savart law. The aerodynamic forces and moments for wings and rotors are computed by strip theory in each model. Second, the performance analysis such as propeller performance and trim for compound aircraft is implemented to check the accuracy between the proposed modeling technique and the helicopter trim, linearization, and simulation (HETLAS) program which is validated. It is judged that this study raises the efficiency of aircraft performance analysis and the airworthiness evaluation.


Introduction
Unmanned aerial vehicles (UAVs) are used for various purposes in military and civilian applications owing to their enormous advantages such as wide range of operational missions and flight altitudes [1]. It is necessary to develop a suitable configuration for the operational purpose and accurately analyze the flight performance as the demand for UAVs in the aviation industry. Today, the various maneuvers such as hovering, low-speed flight, and high-speed flight are simultaneously required to UAVs. However, the fixed-wing aircrafts cannot easily perform hovering and low-speed flight mission. Moreover, the rotorcrafts cannot implement the high-speed flight mission due to the problem of compressibility and flow separation in the rotor. For this reason, many aircraft manu-facturers and researchers have proposed the compound UAVs. They conducted performance analysis and verification for compound UAVs. The most famous configurations are the vertical/short takeoff and landing (V/STOL) wheel [2]. Figure 1 shows the V/STOL aircraft and propulsion concepts.
Basset et al. [3] tried to define the five categories of rotorcraft shapes as shown in Figure 2. The categories mean that all of compound aircraft can be represented by combinations of rotating components such as rotors, propellers, and fixed components such as wings and engines. Thus, the flight dynamic modeling techniques are the most important section when the designer is developing the compound aircraft.
The traditional modeling techniques are considerably depending on the wind tunnel data [4,5]. This method is not efficiency because the experimental using wind tunnel has to be perform for each aircraft. It is the limitation to develop the aircraft system since the development cost and period are highly raised. Also, it is difficult to modify the flight dynamic model (FDM) according to the purpose of analysis or to add functions since the most FDM is distributed as commercial program [6]. The flight dynamic modeling techniques have been developed by several researchers to overcome this disadvantage. Kim et al. [7] studied the aerodynamic and inertial modeling of propeller. Cook et al. [4] and Etkins et al. [5] researched the flight dynamic principle. (1) Cook MV and Etkins used the traditional modeling techniques using wind tunnel data. Leishman et al. [8] tried to develop the mathematical helicopter aerodynamic model. Chaffin et al. [9] defined the guide of the use for pressure disk rotor model. Taamallah et al. [10] studied the flight dynamic modeling for a smallscale helicopter UAV. Howlett et al. [11] developed the black hawk engineering simulation program using mathematical model. Pearson et al. [12] researched the aerodynamic characteristics of tapered wings. Talbot et al. [13] studied the mathematical model for single main rotor helicopter. The researchers in Ref [8][9][10][11][12][13] had tried to develop the mathematical modeling technique for helicopter. However, (2) these researches can be used for certain aircraft. If the configurations of aircraft are changed, the flight dynamic modeling procedure has to be newly started. This is the limitation for developing compound aircraft because it is difficult to develop the flight dynamic model for each compound aircraft. It raises the development cost and period. So, most companies designing the aircraft system cannot try to develop the compound aircraft.
(3) Therefore, the flight dynamic modeling technique has to be developed to easily express the characteristics of all type compound aircraft.
The contribution of this studied is the development of unified flight dynamic modeling technique to overcome the disadvantages of existing modeling techniques and to easily improve the fidelity of FDM while responding to changes in aircraft geometry. Especially, a unified rotor and wing models are developed to selectively calculate the aerodynamic forces of the rotors with wings and to analyze the flight performance of any aircraft without modifying the mathematical model. This technique can easily modify and add the rotor or wing models according to the analysis purpose.
The rest of this paper is organized as follows. Section 2 introduces the unified aerodynamic component modeling technique, and Section 3 deals with the application to validate the accuracy of proposed modeling method by comparing with the HETLAS program. Section 4 is the conclusion of this research.    3 International Journal of Aerospace Engineering 2.1. Unified Rotor Modeling. The rotors can be classified depending on whether they produce the thrust or propulsive forces. The mathematical aerodynamic model of the rotors can be expressed by the blade element method (BEM). However, there are differences between a rotor and a propeller in terms of the mounting position, orientation, and rotor dynamics such as flapping motion. Therefore, various shapes of the rotor and propeller are classified according to the mounting position, orientation, rotor dynamics, and control mechanism as presented in Table 1. Table 1 presents some standardized rotor types. The reference position and direction should be selected to indicate the mounting position and direction of the rotor using the unified rotor model. The aerodynamic forces and moments can be expressed on the body coordinate because all the forces and moments of the aircraft act on the center of gravity (CG). Therefore, the mounting position of the rotor is expressed using the relative distance vector from the CG of the aircraft.
2.1.1. Mounting Position and Orientation. The main rotor of a conventional helicopter is located at the top center of the aircraft and mounted perpendicular to the fuselage or with small forward tilt angles. The tail rotor is generally located on the rear of the aircraft and mounted vertically on the lateral axis. The propeller is mounted on the front (tractor type) or the rear (pusher type) perpendicular to the axis of the fuselage. Here, the azimuth changes for main rotor (MR), propeller, and tail rotor (TR) are expressed as shown in Figure 4.

Anhedral tip sweep
Airfoil #3 Figure 5: General information for rotor model [15].  The mounting angles are defined to represent the mounting orientation, namely, the starboard side tilt angle θ Right , which rotates in the CCW direction on the x-axis, the forward tilt angle θ FWD , which rotates in the CCW (counter clockwise) direction on the y-axis, and finally, the azimuth rotation angle θ AZIM , which rotates in the CCW direction on the z-axis. The passive transformation (alias transformation) of the coordinate system applied to the individual rotor is defined as given in Equation (1) for each axis. The matrices L 1 , L 2 , L 3 are the transformation matrices around x-, y-, and z-axis, respectively.   Engine torque transferred to rotor to overcome main rotor torque The reference for the mounting orientation is the orientation of the conventional main rotor. The coordinate transformation matrix from CG to the rotor hub L HC is defined as given in Equation (2).
2.1.2. Rotor Dynamics. The rotor and the propeller are separated based on the rotor dynamics. The rotor dynamics is not applied when the propeller and ABC (advanced blade concept) rotor developed by Sikorsky are used [14]. Figures 5  and 6 show the general information for rotor dynamics. The unified rotor dynamics can select the rotor dynamics for each rotor. Table 2 is the descriptions for variables of rotor dynamics. The aerodynamic forces and moments of the blade elements are calculated using BEM. The aerodynamic coefficients corresponding to the angle of attack (AOA) and Mach number of the blade elements are estimated from the nonlinear aerodynamic coefficient table con-  International Journal of Aerospace Engineering structed in C81 form using the table lookup method. A uniform inflow model or pitt-peter (dynamic inflow) models is selectively applied to the rotor inflow model as depending on the rotor type, analysis efficiency, and various applications [15]. The kinematics at a blade element has to be defined to express the rotor dynamics. The coordinate of the rotorcraft is defined as shown in Figure 7, and the basic velocity, acceleration, and angular velocity are expressed as Eqs. (3) and (4). The superscript CI is used to define the vectors at the rotorcraft CG relative to the center of the inertial reference frame, whereas the subscript represents the frame (C for CG-fixed frame) used to define vector components. For instance, the vector ω CI C is the angular velocity of the CGfixed frame relative to inertia frame, which is expressed on the CG-fixed frame. Figure 7 shows the front view of the rotor fixed nonrotating coordinates ðx H , y H , z H Þ with rotational speed Ω in the CW direction.
The linear velocity and acceleration at the rotor center can be expressed by using the definition of the rotorcraft's CG position r CI C and the center of the rotor r HC C as shown in Equation (5). The subscript R in Equation (5) means the rotating hub frame of rotor. The equations of the rotor dynamics have to reflect the order of the hinge changes, adding, and subtraction. It means that the rotor dynamics cannot be expressed by one formula without constraints. This study utilizes the generalized vector kinematics [7]. The generalized vector kinematics is derived by relative motion from several coordinate systems as shown in Figure 7. Also, the generalized vector kinematics is utilized to derive the implicit formulation of the rotor motion equations. It can be applied regardless of the shape and arrangement order of various types of rotor. The accuracy of the derived rotor motion equation can be easily verified. The linear velocity and acceleration are described at the lead-lag hinge because the feathering motion is not applied to rotor dynamics in this study. Moreover, the lead-lag hinge is the outer most hinge as shown in Figure 5. It is decided to denote the variables related in feathering motion in Equations (6) and (7) for convenience. The linear velocity and acceleration at the lead-lag hinge for blade section can be expressed as Equations (6) and (7), respectively, where the subscripts L and T mean the lead-lag hinge and the blade section frame, respectively. Furthermore, the linear velocity and   International Journal of Aerospace Engineering acceleration are defined on CG and AC of the blade section to reflect the rotor dynamics and to calculate the aerodynamic forces, respectively.
The velocity _ r TI L is composed by three parts such as linear velocity, angular velocity from lead-lag hinge, and effect of the angular rate from lead-lag hinge to the blade section. Although the acceleration as Equation (7) can be derived by using similar way to define Equation (6), the only difference between Equation (7) and Equation (6) is that Equation (6) considers the sweep length y S and an hedral length z A . The variable σ denotes the rotating direction which is CW rotation (σ = −1) and CCW rotation (σ = 1).
The flap, lead-lag, and RPM dynamics associated with rotor motion are defined as Equations (8)-(10), respectively. Each dynamics can be selectively used by rotor type. If the rotor type is defined as the propeller, the rotor motion does not apply as noticed in Table 1, whereas if the general rotor is selected, the rotor motion is used. Moreover, the relative angles such as flapping angles are computed with each rotor motion. The flapping motion occurs of the y-axis associated with rotating rotor hub frame because of the flap hinge. The unit vector of y-axis and z-axis on rotor hub frame is defined as e Y and e Z , respectively. The lead-lag and RPM dynamics move on the z-axis as the same principle of flapping motion. The flapping motion and RPM dynamics are needed to consider rotating direction additionally. The   Figure 23: Concept of strip theory [16]. Figure 27: AOA with V c < 0, V n < 0. 9 International Journal of Aerospace Engineering RPM dynamics does not consider the damping ratio and spring constant because the rotor shaft linked directly to the gear box and engine. So, the inertia forces in Equation (10) are divided as two forces generated by acceleration at the blade section and shaft I DRV _ Ω, respectively. The variables which use in the following equations are listed in Table 3.
The aerodynamic forces dF Aero Sec and moments dM Aero Sec for the unit section on the blade have to be calculated to compute the flap, lead-lag, and RPM dynamics. Figure 8 shows the general information for generating aerodynamic forces with induced drag at the blade section. The aerodynamic forces generated by blade section can be defined as Equation (11).
where dR, dL, dD, and dM are the spanwise axial forces, lift, drag, and moment on the unit section for blade,respectively, to calculate the aerodynamic forces. These four variables are expressed as Equation (13).
Here, c is the distance of the camber for blade. C L , C D , and C R are the aerodynamic coefficients for lift, drag, and spanwise axial force, respectively. The coefficients are estimated using C81 tables by v TI L , v tur L , and v ind L mentioned in Equation (14) are _ r TI L , turbulence, and induced inflow, respectively.     (1) are applied for the rotation of the wing. Thus, the wing coordinate rotation matrix L WING using Equation (1) is defined as given in Equation (15).
The unified wing model can handle the angle of elevation, attachment, and rotation angle for wing. The elevation angle θ Elev and the attachment angle θ Attach for wing can be defined as shown in Figures 9 and 10. The fuselage line and the starboard side of the wing are the standard line to define each angle. Figures 11-14 show the results of wing orientation for elevation angle θ Elev . The blue color expresses the port wing, and the red color means the starboard wing in Figures 11-14. Figure 11 shows the wing orientations on the X-Y-Z coordinate when the elevation angle for starboard side wing and port side wing are defined as 45 (deg) and 0 (deg), respectively. Moreover, the attachment angles are same as 0 (deg) for port side and starboard side. Figures 12-14 are expressing the wing orientation results on the X-Y, Y-Z, and X-Z coordinate, respectively. This is the very useful function to make the geometry for various compound aircraft. Figures 15-18 show the wing geometry when the attachment angles of starboard side and port side wing are defined as 45 (deg) and 0 (deg). The difference between the starboard side and the port side is expressed on the X-Y-Z coordinate in Figure 15     Finally, the unified wing model considers the upsidedown configuration of airfoil. This configuration is mostly applied to horizontal stabilizer. The horizontal stabilizer is attached to generate a positive pitching moment. It is the same as when the direction of the normal vector of the airfoil is reversed. Figures 19 and 20 show the normal vector u n and tangential vector u c on the X-Z airfoil coordinate.
u n and u c on the wing coordinate as shown in Figures 9 and 10 can be defined as Equation (16), where e x , e y , e z are unit vectors of each axis on the wing model.
The upside-down airfoil configuration changes a direction of normal vector u n as shown in Equation (17). The aerodynamic forces and moments can be easily computed by Equation (17).

12
International Journal of Aerospace Engineering F C means an origin point of the body coordinate. r P and r v are a distance vector from the origin point of the body coordinate to control point of each lattice and to unit vortex filament dl, respectively. r c is a core radius for vortex filament, and h is a perpendicular distance from vertex filament to control point P. If the number of lattice is increased, the number of control point is expended. The lift of each lattice can be computed by Equation (19).
where ρ and V are the density of air and the true air speed, respectively. Γ c is the vortex strength. Kutta-Joukowski theory as shown in Equation (20) is used to compute Γ c .
C L and s are the lift coefficient and the area of unit section. Finally, the velocity V has to be computed, and it is a kV ac k. V ac is a velocity vector that the Y-component is zero. The velocity vector is calculated after computing v k . v k is a velocity vector including the Y-component on each lattice. The subscript of v k is a number of lattice in the main wing and fin stabilizer. v k is defined as Equation (21).
where v CI C is a true air speed vector associated with CG, and ω CI C is a angular velocity vector of aircraft. r PC C is a distance vector between the CG of aircraft and the control point e y is defined as shown in Equation (23).
The calculated aerodynamic forces and moments for each lattice have to be expressed on the body coordinate. This study use the strip theory to incorporate the aerodynamic forces and moments on the origin point of the body coordinate. The CG is the same with the origin point of the body coordinate. Figure 23 shows the concept of strip theory.
Thus, the computed aerodynamic forces and moments of each lattice can be expressed on CG as shown in Equation (24). Here, α and β are the AOA and the sideslip angle (SSA), respectively. 14 International Journal of Aerospace Engineering The aerodynamic forces and moments of main wing and fin stabilizers can be calculated by applying Equation (24) to all of lattice as shown in Equation (25).
where the subscript j of the P is order of lattice for main wing and fin stabilizer. So, jth control point P is described as P j .

Unified Wing
Modeling: Aerodynamic Coefficients. The AOA is calculated by 4 cases by using the absolute values of velocities. Each case is described as shown in Figures 24-27. The influence of the wing lift coefficient on the control surface deflection is expressed as given in Equation (26). The control surface is considered as a plain flap. Here, C l,surface and δ surface are the lift coefficient associated with control surface and control surface deflection angle, respectively.
The main coefficient in Equation (26) is a C l,surface because other coefficients such as C l,0 and C l α are the constant or the functions associated with AOA, SSA, and Mach number. Thus, the calculation accuracy of C l,surface is the most important for the wing model. Several tools such as chester, torenbeek, and DATCOM have suggested the theory to calculate the coefficient C l,surface . A comparison of various empirical formulas calculating C l,surface indicates that the theory proposed in DATCOM is the most accurate as shown in Figure 28 [17]. The reference data in Figure 28 is the wind tunnel test data for NACA23012 airfoil. The specification of NACA23012 airfoil is described in Table 4. This study selects the theory proposed in DATCOM to compute C l,surface .

Application
This section tries to validate the FDM built by the unified rotor and wing model. The single rotor helicopter and compound helicopter using coaxial rotor are used for checking fidelity of FDM. The helicopter models are used to the BO-105 helicopter and the reference rotorcraft which is built by HETLAS. A reference rotorcraft is a Surion helicopter. HETLAS is a program for implementing trim, linearization, and simulation of the helicopter. Especially, this program had been validated by the Surion helicopter. Moreover, it has been already used to develop the FBW (fly-by-wire) system for Surion helicopter. Figure 29 shows the structure of the HETLAS program. The detail information of HETLAS is described in Ref. [15,18].
The configuration data for reference helicopter are described as Table 5, and Figure 30 shows the shape of reference helicopter.
The configuration data and shape for the BO-105 helicopter are described as Table 6 and Figure 31, respectively.
The validation of proposed FDM would be conducted by comparing with analysis results computed from the HETLAS program. Trim analysis of the reference helicopter and simulation for the BO-105 rotorcraft in the time domain are used as the validation method for proposed FDM.
The trim analysis of the reference helicopter would be implemented by the HETLAS program and proposed FDM, where the validation of fidelity for proposed FDM is conducted by checking the difference of trim results between the HETLAS program and proposed FDM since the almost 15 International Journal of Aerospace Engineering analysis and experimental data of reference helicopter are the confidential data. The simulation of the BO-105 helicopter would be performed by proposed FDM because the flight test data of the BO-105 helicopter has been opened. Especially, the AC-120 (advisory circular 120) criteria are applied to the simulation process to compare thoroughly.

Unified Rotor
Model. This study shows that the rotor model employing the unified technique can be used to analyze various types of rotor dynamics. This study uses the rotor model validated by the HETLAS program. Figures 32  and 33 show the aerodynamic analysis results of the rotor and propeller using the unified rotor model.
The rotor of a typical helicopter is used for analyzing rotor performance as shown in Figure 32. The main rotor is used for analysis, and the collective inputs are defined as 5 degrees under the normal working states with increased forward speed. All results are compared with no rotor dynamics, flap dynamics, and flap and lag dynamics. The induced velocity has been decreased when the flight speed is raised. The induced velocity is inversely proportional to forward speed as v i = T/2ρAV ∞ as understanding the momentum theory. The induced velocity is significantly decreased as forward speed increased as shown in Figures 32 and 33. Although the induced velocity is similar responses, the torques represent different behaviors. Torque with no dynamics is steadily increased, but the others have maintained those values. This is because the AOA of the blade section affects the drag force of the blade section. The drag coefficient at the blade section is reduced with a small effective AOA due to the increasing forward speed U T and the decreasing induced velocity U P as shown in Figure 8. This is the main reason to maintain the torque level. Furthermore, the flapping angles which can be considered as damping effect make the significant bias among all speed regions. The lag dynamics with negative lead-lag angles show less effect to torques increasing as a slightly bigger effective AOA with changing U T .
The tractor-type propeller is used for analyzing propeller performance as shown in Figure 33 [7]. It shows the result of the analysis associated with advanced ratio and collective change. The thrusts and the required powers of the propeller have been increased when the advance ratio and collective input are raised, because the increments of flow velocity for propeller make the effective AOA of the blade section falls in the stall region. It means that the maximum thrust level is almost constant when the effective AOA is in the stall region. Furthermore, the propeller efficiency can be reached to maximum point. Figures 32 and 33 have showed that the rotor dynamics can be selectively applied. Therefore, it is judged that the unified rotor model is highly effective to calculate the aerodynamic forces of various types of rotor or propeller.

Unified Wing
Model. This research tries to validate the fidelity of the wing model. The tapered wing is used for checking the fidelity of the wing model. The specification of tapered wing without control surface is defined as Figure 34. The airfoil is a NACA 23012, and the wind velocity kVk is a 21.34(m/s).
Figures 35 is a lift coefficient and drag coefficient. These two aerodynamic coefficients computed by the wing model are very similar with reference data. The reference data is a wind tunnel test data for tapered wing.
The coefficients computed by the unified wing model for lift and drag are very similar with reference data. Thus, it is judged that the fidelity of the proposed wing model is appropriate. The specification of tapered wing with control surface is defined as Figure 36 [22]. The validation is performed by the BO-105 helicopter and its flight test data. This validation is checked to meet helicopter simulator qualification [23]. Table 7 shows the tolerances with control inputs, attitude angles, and handling qualities in each axis.

18
International Journal of Aerospace Engineering Therefore, the aircraft performance can be accurately analyzed by proposed FDM. Furthermore, if the aircraftspecific adjustments are applied to FDM, a more precise prediction of performance can be made. Figures 40-43 indicate the response property of aircraft for collective pitch input. The input signal is given as 3-2-1-1 for each axis. Also, it shows the handling qualities criteria of collective pitch input for longitudinal, lat-eral, and directional axis. The simulation has been conducted in forward flight with 80 knots for 10 seconds. The related responses, for example, pitch rate in longitudinal handling qualities with lateral cyclic inputs as described in Table 7, show the very similar results with flight test data even though they have small errors. Therefore, it can be also considered that the FDM can precisely predict the flight performances.  The compound aircraft using coaxial rotor is based on the single rotor helicopter used in HETLAS. The coaxial rotors were assumed to be ABC rotors, and the blade chords were adjusted to have the same solidity. In addition, the wings were attached to the center at both sides of the fuselage, and a pusher propeller was attached to the rear fuselage. The control inputs were collective, longitudinal and lateral cyclic pitches, and differential collective pitch as given in Equation (27). Furthermore, the collective control of the pusher was applied through thrust scheduling as shown in Figure 45 to create appropriate propulsive force according to the forward flight velocity.
Figures 46 shows the trim results of compound aircraft. The compound aircraft has a high pitch attitude angle with small roll attitude angle and lateral cyclic input due to torque from the pusher at the hovering status.
Then, the collective input is getting decreased with the pitch angle until 80 knots because of the lift of wing    21 International Journal of Aerospace Engineering increased. After 80 knots, the pitch angle decreased, and the pusher collective input increased more steeply. This makes lateral cyclic inputs increased due to high torque from the pusher. Although there is the complicated interplay between the longitudinal axis with collective and longitudinal cyclic input, the lateral cyclic input is working for balancing the rolling moments. Therefore, it is judged that this simulation results are physically appropriate, and the proposed FDM can be widely used to analyze the flight performance of compound aircraft.

Conclusion
This paper developed the unified high-fidelity flight dynamic modeling technique for compound aircraft. The unified rotor model calculated the aerodynamic forces and moments by using BEM with enhancement of the coefficient estimation. Furthermore, the coefficient interpolation technique has been applied to the rotor sweep, its orientation, and rotor dynamics. The proposed FDM uses the aerodynamic coefficients of airfoil for rotor and wing. The aerodynamic coefficients are the necessary data to build the high-fidelity FDM. This is the same with the traditional FDM technique. The analysis performance of unified FDM do not depend on the configuration of aircraft such as wing span and rotor diameter. Especially, although the aircraft has not manufactured, the performance analysis can be performed by unified FDM. This is the main advantage of this research, whereas the traditional FDM using experimental data has the dependency of configuration for aircraft since the experimental results is different. This is the limitation of development for compound aircraft as mentioned in section of introduction. Thus, once the aerodynamic coefficients for airfoil are given to proposed FDM, most shape of rotor and wing will be accurately analyzed.
The unified rotor model was defined as general rotor and propeller to analyze its performance as shown in Figures 32 and 33. The analysis results for rotor and propeller showed the physically appropriate rotor performance as mentioned in section of application. The unified wing model computed its aerodynamic forces and moments using strip theory and VLM using Biot-Savart laws to estimate induced velocity. Moreover, it included the wing orientation as well as control surface influences when the induce velocity was calculated. The unified wing model was defined as the tapered wing model to check the property of lift and drag. The lift and drag coefficients calculated by the unified wing model are highly accurate. It is judged that the unified rotor and wing model proposed by this research are very precise and useful.
The unified rotor and wing models were used to build FDM. The validation of fidelity for the proposed modeling technique was performed by single rotor helicopter and compound aircraft using coaxial rotors. The applications were defined as trim and simulation with flight test data. A reference helicopter including confidential data and the BO-105 single rotor helicopter was used in section of application. First, the results of trim analysis for the reference helicopter between HETLAS and unified FDM were practically the same since the small errors were induced by enhancing the estimation accuracy of aerodynamic coefficients. Second, the correlated simulation responses with control inputs which were 3-2-1-1 control to collective, lateral cyclic, longitudinal cyclic and pedal inputs. The results of simulation computed by proposed FDM for BO-105 with the flight test data showed the simulation results of proposed FDM for BO-105 that are satisfying the qualification for helicopter simulator with small error. Third, the flight performance of compound rotorcraft using coaxial rotors was analyzed using the unified FDM as shown in Figures 45 and 46. Figures 45  and 46 show the physically appropriate results as mention in application section.
Therefore, all type aircraft can be accurately analyzed by this research. In addition, the analysis efficiency of the compound aircraft system can be dramatically raised since the flight performance based on airworthiness standards can be interpreted by nonlinear optimal control programming (NOCP) or nonlinear programing (NLP) method. It is judged that this study can affect to develop all types of aircraft and can reduce the development cost and period by providing high-fidelity analysis results for flight performance.
Data Availability [1] The data used to support the finding of this study are available from the first and corresponding authors upon request.