Optimization Plug Mode of External Fuel and Weapons for Less Changing Aircraft Center of Gravity

To achieve the minimum change of aircraft center of gravity, plug mode optimization was carried out. Plug mode, including weight and position of external fuel and weapon load, was selected as design mode. Numerical model for calculating aircraft center of gravity was established and used to generate mode samples. Random search algorithm (RSA) was designed and applied for solving the optimization plug mode, and the precise solution of plug mode was determined by full array method (FAM). Compared with typical plug mode, the optimization mode can reach a minimum change in aircraft center of gravity with fewer samples than the solution by FAM. RSA coupled with FAM is an effective scheme for the optimization of aircraft plug mode.


Introduction
Demands for more fuel and weapons plug and smaller change of aircraft center of gravity have caused aircraft designers and staff to develop system or model to determine aircraft center of gravity.Han et al. developed an adaptive nonsingular fast terminal sliding mode control for aircraft with center of gravity variations to improve the dynamic and the quality characteristics of the aircraft [1].Yang et al. established aircraft center of gravity estimation using Gaussian process regression models [2].An aircraft weight and center of gravity indicator system was developed [3].Gupta et al. proposed a finite element based structural model of a small flexible aircraft, focusing on obtaining a simple model suitable for control design based on a two-step procedure [4].Li et al. analyzed the position of center of gravity for an unmanned aerial vehicle based on load task [5].
Much work has been done on the influence of aircraft center of gravity on flight performance and aerodynamic characteristics.Li et al. did research on effects of partial fuel pump failure on center of gravity control for high-speed aircraft to gain potential performance benefits.Besides, the influence of partial fuel transfer line leak on active center of gravity control was analyzed [6,7].Jungo et al. investigated the development of the CEASIOM aircraft design environment for novel aircraft configurations and analyzed static and dynamic stability of the aircraft [8].Mardanpour et al. did research on engine placement effect on nonlinear trim and stability of flying wing aircraft [9].Zhen probed into the CG position of minimum drag to check stability boundary and front control boundary for J-7 aircraft [10].An analytical expression was developed which showed that wing/tail interference drag was determined by wing downwash at downstream infinity [11,12].Wang et al. did optimization design of unmanned aerial vehicle with considering the uncertainty of the center of gravity and the influence of the center of gravity on the navigation time is greatly reduced [13].Diekmann did trim analysis of nonlinear flight dynamics for a civil aircraft with active high-lift system [14].Tang et al. put forward a calculation and analysis method for aircraft center of gravity to improve flight performance and safety [15].
Developing an efficient design and optimization method for plug systems is attracting more and more attention from researchers.Dancila et al. researched fuel burn prediction algorithm for cruise, constant speed, and level flight segments and considered the continuous fuel burn rate variation with time caused by the gross weight modification due to the fuel burn process itself [16].Patron et al. put forward low calculation time interpolation method on the altitude optimization algorithm for the FMS CMA-9000 improvement on the A310 and L-1011 aircraft [17].lavatory partition for an aircraft [18].Jia et al. developed airborne weapon CAPP system and met the application of an aviation enterprise [19].Xia et al. investigated the research of release safety analysis method for external stores of weapon [20][21][22].Luo developed 1553B bus error correction technique in the management system of aircraft weapon to increase the adaptability between the aircraft and external plug [23].
Concentrating on the plug mode based on various missions for fighters with much heavy fuel load or quite a few weapons, a method based on RSA is designed and used to solve the optimization plug mode for minimum changing aircraft center of gravity.
In this paper, the design mode and objective function are presented in Section 2; and then the method for solving the optimization problem is performed in Section 3; the optimization results are analyzed in detail in Section 4; finally, the manuscript is concluded in Section 5.

Design Mode and Objective Function
The plug mode is defined as the corresponding relations between external loads and plug positions where the external loads include fuel and weapons loads such as common size auxiliary fuel tanks, various missiles, and bombs.The schematic of plug and spare load for a fighter on active service in current research is shown in Figure 1.European EF-2000 is selected as the fighter support, with empty weight 11.150 ton, 13 plug positions, and two types of external fuel tank.The coordinate of aircraft center of gravity can be calculated as where   is the coordinate of aircraft center of gravity,   is the aircraft empty weight,   is the coordinate of aircraft center of gravity with empty weight,  if is the total weight of internal fuel,  if is the equivalent coordinate of internal fuel,  ef is the external fuel weight,  ef is the coordinate of the external fuel plug,  ew is the weight of external weapon,  ew is the coordinate of external weapon,   is the total weight of the current aircraft.Generally, the position of aircraft center of gravity is defined as a percentage of mean aerodynamic chord length: where  is the position of aircraft center of gravity,   is the coordinate of the leading edge of mean aerodynamic chord,  mac is the length of mean aerodynamic chord, and   is the trailing edge of mean aerodynamic chord.By using symmetrical loading, there are 7 plug positions to carry weapons or fuel tanks as shown in Figure 1.With specific task, external fuel and weapons can be determined.Then plug mode can be designed as The plug mode, namely, the one-to-one correspondence between   and   ,  = 1, 2, . . ., 7, can be solved by determining the array of   and the center of gravity can be calculated by the combination of  and  based on formula (1).

Methods for Solving Optimization Problem
Because of various types of missiles and different missions, the group of  and  can reach hundreds of thousands or more, from the theoretical perspective regardless of plug position limit, including the fixed position of heavy external fuel tank or some anti-ship missile.Supposing the weight of internal fuel equals zero and all the plug positions mounted with weapons or fuel tanks, the optimization plug mode can be solved by the following methods.

Full Array Method.
For giving mission and fixed type of external fuel load and weapon, the number of the group of  and  is an arrangement number, from the theoretical perspective.And a fixed plug position can be considered as cutting the length of the array of .Then the matrix of all the group of  and  can be expressed as where FAM is a function to list all the array of ,  is  1 × 7 matrix,  0 is the number of plug modes in different missions,  1 is the number of plug modes in giving missions,  2 is the number of plug modes by RSA,  is the number of active plugs, and  is the number of all the types of external fuel tanks and weapons.According to (1)-( 3), the center of gravity and minimum change of center of gravity can be calculated as where  0 is the aircraft center of gravity with a typical plug mode under the giving mission.

Random Search Algorithm.
Taking into account EF-2000 fighter, the value of  can reach 17; then it is formidable for FAM to couple with the huge number  0 .Therefore, RSA is designed to search the optimization plug mode in a small samples.For a giving mission and optional types of external fuel load and weapons, RSA is designed and used to solve the optimization plug mode as shown in Figure 2.
Here  2 is a random number or a giving number,  is critical error of center of gravity change, and  is  2 × 7 random matrix by arrangement of  as described in (7).Similarly, the center of gravity and minimum change of center of gravity can be calculated as The RSA operator can proceed as follows.
Step 1. Obtain or update  2 value.
Step 2. For each  = 1, 2, . . .,  and  is a true random operator to generate a set of nonrepeating permutations to simulate the various plug mode.
Step 3. Calculate  and  for each   .
Step 4. Determine the optimal plug mode with the minimum absolute value of  among .
Concerning that the valve of  2 is not large enough to find the optimization solution,   2 is used to update  2 to search one more time until meeting the error test.Even if the above 4 steps still cannot get the answer, resetting  can be available to continue the algorithm.The updating and resetting operator for  2 can be described as follows: The optimization plug mode can be determined based on (9) with the termination of the random search algorithm.In addition, the ultimatum solution of the optimal plug mode should pass the results test such as flight qualities and operational and wing structural constrains as shown in Figure 2.

The Results of Typical Plug Mode.
For comparison of the solutions between RSA and FAM, a typical plug load including 10 missiles and 3 external fuel tanks is given as shown in Table 1.There are 10 missiles and 3 fuel tanks for the fighter to possess plenty of missiles and fuel.By using symmetrical loading, the typical plug mode is presented in Table 2.

The Results and Analysis of FAM.
To obtain the precise solution of the typical plug mode, FAM is used to produce all the array of  matrix.Concerning the fact that there is a unique position ( 7 ) for fuel tank 2 as shown in Figure 1 and Table 2, the value of  1 is 720.With this in mind, the precise plug mode can be determined as presented in Figure 3.The results of center of gravity are obtained as shown in Figure 4 with a range of 13.1343%-27.8909%.The change of  is presented in Figure 5 in the scope [−6.3657, 8.3909].To observe the range of , the results of center of gravity is sorted as shown in Figure 6.Besides, there are 164 plug modes with less changing of the value of  than the typical plug mode.Then the minimum change of  can be determined as 0.0594% where  = 19.5594%with the precise plug mode.Because  matrix is a full arrangement of , the precise solution can be considered reliable.

The Results and Analysis of RSA.
Figure 7 shows that it is formidable to search the robust results of plug mode when the initial size of  2 is small.Thanks to the action   of updating and resetting operator of  2 , the relative error of the optimal  is significantly reduced when the value of  2 increases to 20 as shown in Table 3. Figure 8 presents the fact that there are more solutions close to the value  = 19.5594%by FAM.In addition, the error is significantly improved when  2 ≥ 20 from Table 3 and the proposed algorithm is theoretically and practically convergent.Figure 9 supports the fact that a new plug mode can be applied to finish the external load distribution when it is suitable for  6 plug position to carry Paveway missile, including the strength of the plug pylons and the space around the plug point.Unlike the other random search method for producing some single values, the RSA operator concentrates on generating a certain number of nonrepeating combinations to simulate the various and promising plug modes.It can be noted from Figure 10 that the center of gravity ranges from 12.4865% to 22.4288% close to the result by FAM. Figure 11 shows that the value of Δ is between −7.0135% and 2.9288%.It could be remarkable from Figure 12 that there are 59 plug modes which changes center of gravity less than the typical plug mode by comparing the absolute value of Δ between RSA and the typical plug mode.The minimum value of center of gravity position can reach 19.5605% with Δ = 0.0605% which is close to the solution by FAM.Furthermore, it just costs RSA 100 initial plug modes to get the optimization solution by searching only one time.Compared with FAM, RSA is really fast and efficient with enough precision.Figure 13 supports the fact that there are 13 plugs occupying all the plug positions not only to simulate aircraft high payload but also to test the proposed algorithm.By pouring attention to the weight of Paveway missile heavier than that of fuel tank 1 in spite of the former volume smaller than the later,  2 or  4 plug position is strong enough to carry the other plugs lighter than fuel tank 1 from the perspective of wing structure and this is operable.In addition, the high weight of Paveway plug at  2 position is not good for rolling manipulation, thus affecting the lateral flight qualities.Figure 14 indicates that a precise plug mode with the minimum changing center of gravity is obtained by FAM.Although the longitudinal flight qualities are improved because of the minimum changing aircraft center of gravity, the heavy weight of Paveway missile at  1 position is extremely detrimental to the lateral flight qualities as shown in FAM results.In addition, high load often occupying the wing tip is unfavorable for aircraft wing structure and easily leads to structural fatigue along with the aircraft taking off, climbing, and maneuvering.Thanks to the optimal solution by RSA, the disadvantages above are well improved as shown in RSA results.As previously demonstrated,  2 is strong enough to carry fuel tank 1 lighter than the Paveway missile and the auxiliary tank could be preferentially exhausted and abandoned; thus the flight qualities could be improved significantly.The remaining operation is placing the Paveway missile at  6 position and this is convenient and feasible.

Conclusions
RSA for solving the optimization plug mode was carried out.A typical plug mode was used to be the comparison of reference between RSA and FAM.FAM is precise and comprehensive to calculate the precise solution with computing a large number of the arrangements of .RSA is accurate and efficient for generating the optimization plug mode with searching among a small number of samples of the arrangement of .From this research, it can be found that RSA coupled with FAM is an efficient tool for optimization of aircraft plug mode.

Figure 2 :
Figure 2: The flow chart of the random search algorithm and results test.

Figure 5 :
Figure 5: Change of  by FAM.

Figure 6 :
Figure 6: Sorting results of center of gravity by FAM.

Figure 13 :Figure 14 :
Figure 13: The schematic of aircraft, plug, plug position, and typical plug mode.
Ivester et al. worked on compact aircraft galley and lavatory arrangement and articulating

Table 1 :
A typical plug load.

Table 2 :
The distribution of the typical plug load.

Table 3 :
Center of gravity generations of the optimization plug mode by RSA.