# I. Introduction ith the increasingly complex battlefield environment, the war has put forward higher requirements for the mission capability and viability of UAVs. The cost of large strategic deterrence UAVs continues to rise, and the cluster UAVs have gradually become a new focus [1] . However, due to the cost and volume constraints, cluster UAVs are inherently short of range and endurance. At present, there are three common methods to improve the range and endurance of the aircraft: First, the aerial refuelling technology is adopted, but for a large number of cluster UAVs, the efficiency of aerial refuelling cannot meet the actual needs. The second is to use a large aspect ratio design combined with solar power module, but for clustered UAVs that need high flexibility and mobility, the aspect ratio will be greatly limited, while the small aspect ratio combined with solar power module will have very low income. Third, like the "Gremlins" in the United States, large transport planes are used for long-distance delivery, but this only solves the problem of long distance, not the endurance time. In recent years, with the rise of clustered UAVs, a technology that once appeared has been re-valued, namely wingtip connection technology. Scientists hope that the fixed-wing cluster UAVs can form a combination with a large aspect ratio through wingtip connection. At this time, if the UAVs combined with solar technology, the range and endurance time can be greatly extended [2] . Wingtip connection was first proposed by Richard Vogt, a German scientist who immigrated to the United States after World War II. Based on this concept, the United States carried out "MX-1018" (Tip-Tow) and "Tom Tom" projects. However, due to the aerodynamic interference at the wingtip, this research has caused great potential safety hazards and even casualties. In addition, the impact of the gradual maturity of aerial refuelling technology at that time temporarily caused the flight test of wingtip connection stop [3] . Since the 21st century, scientists in the United States, Germany and China have seen the great value of wingtip connection, see Figure 1, for it can improve the ceiling and endurance of cluster UAVs, and gradually resumed and accelerated the study on wingtip connection. Through computational fluid dynamics analysis and wind tunnel test, American scientists analysed and studied aerodynamics and aerodynamic control in the process of wingtip connection [4][5][6][7] , and verified the flight control system through flight test [8][9][10] ; German scientists mainly explored and practiced the wingtip connection technology in engineering [11] , and focused on designing the control system of the combination [12] ; Since 2018, China has started the exploration and study of relevant technologies, demonstrated that wingtip connection can improve range and endurance [13,14] , and conducted relevant flight tests [15,16] . In this paper, through the study on the aerodynamic theory of the UAV combination, simulation and example verification are carried out on the improvement of the ceiling and endurance of the wingtip connection, so as to verify that the clustered UAVs can improve the ceiling and endurance by wingtip connection. # II. Basic Theoretical Research on Flight Performance of Wingtip Connection Taking the conventional aircraft with a wing, without offset, dihedral and twist, as an example, the aerodynamic performance of the combination formed by the wingtip connection of the flying units with identical n frames is analysed theoretically. # a) Effect of Wingtip Connection on Lift-drag Ratio Lift-drag ratio is the most intuitive embodiment of aircraft aerodynamic performance. The larger the lift-drag ratio, the better the aerodynamic performance. In this section, it is assumed that the external environment before and after wingtip connection is unchanged, and the flight altitude and speed of the aircraft are unchanged. When the flight state is approximately at the minimum resistance, the endurance can reach the maximum. This is an ideal state, which can predict the trend and simplify the formula. For this state, there is always a waste resistance coefficient, C D0n , equal to the induced resistance coefficient, C Din , because the induced resistance is inversely proportional to the square of the velocity, and the waste resistance is proportional to the square of the velocity. When the two are equal, the total resistance is the smallest, see Figure 2, that is The lift-drag ratio, K, can be derived from equations ( 1) and (2). 0 0 2 Dn Din D n D n C C C C = + =(1)0 0 4 4 Ln n n n n Dn D n D n C A e A K ne C C C ? ? = = = ?(3) Under the assumptions in this section, the parasitic resistance coefficient of the combination is equal to the coefficient of the flight unit, which is a constant, that is 0 0 D n D C C = 0 4 D n A C ? . n n K ne ?(4) That is, the lift-drag ratio of the assembly K n is proportional to n ne . Here, the Oswald efficiency factor can be calculated according to the formula mentioned in literature [2], that is ( ) 0 1 1 n d e nAkC s ? µ = +(5) Where, according to the Oswald efficiency coefficient of the high aspect ratio wing, k is assumed to be 0.18 [17] , C d0 is the resistance coefficient of the airfoil.? is the theoretical Oswald coefficient without considering the influence of the fuselage, s indicating the influence of the fuselage on the induced resistance, which is a function of the fuselage diameter d and the wingspan b. ? ? = -? ? ? ?(6) For the specified flight unit, in equation ( 5), all are constants except e n and n, so it is advisable to set ( ) 1 0 2 1 d k AkC k s ? µ = ? ? = ? (7) Then equation ( 5) can be reduced to 1 2 1 n e k n k = +(8) It can be seen that Oswald efficiency factor decreases with n increasing. # b) Effect of Wingtip Connection on Cruise Speed This section assumes that the maximum liftdrag ratio point of the flight unit is taken at 0°, so that the cruise speed of the flight unit is the speed of horizontal flight at 0° angle of attack. It is assumed that the force balance of flight units before and after the wingtip connection is always satisfied, and the flight altitude remains unchanged. Define that in the two-dimensional plane, the horizontal flight direction of the aircraft is the x-axis positive direction, and the y-axis is perpendicular to the x-axis and upward. Carry out force analysis on the X and Y directions of the aircraft respectively, and then equation ( 9) is on X axis, equation ( 10) is on Y axis. 2 1 2 n n n n Dn n F D V C S ? = = (9) 2 1 2 n n n n Ln n G L V C S ? = =(10) It can be seen from equation ( 9) and (10), n n n F G K = (11) That is, when the maximum lift-drag ratio is reached, the cruising pull is the smallest, so that the endurance time is the longest. The lift coefficient calculation formula given by reference [15] under the constant speed assumption. 1 1 2 2 2 1 1 = 2 2 2 i i Ln n n n n nL L L L C V V nS S nS ? ? - - -? ? ? ? ? ? ? ? ? - ? ? ? ? ? ? ? ? ? ? ? ?(12) It can be seen that the lift coefficient increases with the increasing of n. Where L i is the ideal lift of a single wing calculated according to the two-dimensional infinite length model, and L ? is the lift difference for one side of a single wing. Under the assumption of this section, when G n unchanged, V n gradually decreases with the increase of C Ln , and equation ( 10) can be reduced to equation (13). 1 n Ln V C ?(13) That is, the cruising speed of the combination is inversely proportional to Ln C . # III. Simulation Study on Flight Performance of Wingtip Connection The vortex lattice method is used to simulate the aerodynamic performance of the wingtip connection. Under the constant lift mode, the simplified full USA-35b airfoil model of a rectangular wing with 3 m wingspan and 0.4 m chord length is calculated, as shown in Figure 3, assuming that the wingtip is rigidly connected and there is no gap. The mass of the flight unit is 60 kg, when the angle of attack is ? = 0° and the torque coefficient C m = 0 at approximately the maximum liftdrag ratio. © 2023 Global Journ als The approximate equivalent method of direct size increase is used in the calculation of the combination, which will lead to some errors in the simulation results. This method is used to simulate the changes of lift-drag ratio and flight speed with different number of units. The simulation results of aerodynamic characteristics at an altitude of 2000 m are shown in Figure 4. 1 for the sorted data, and draw Figure 5 for the simulation lift-drag ratio at 0° and the calculated theoretical lift-drag ratio. It can be seen that the law of lift-drag ratio approximately meets n times, and the failure to reach n times is caused by e n < 1. This result is consistent with the conclusion in Section 2.1. Shown in Figure 6, two parameters of lift coefficient and flight speed with different n are extracted from the simulation results. It can be seen that the demand for cruise speed gradually decreases with the increase of n, thus reducing the consumption of power and energy. The relationship between C Ln ?V n 2 and n is shown in Figure 7, verifying equation ( 13), which is consistent with the conclusion in Section 2.2. In this section, it is assumed that the flight unit or combination will fly at a constant speed and level under the condition of force balance. When n is constant, both G n and S n are constant, too. Then let equation ( 10) be reduced to 2 3 n n Ln k V C ? = . 3 1 2 n n n G L S k = =(14) It can be seen that G n , S n and k 3 remain unchanged when the height changes. Then equation ( 9) can be reduced to 3 1 1 2 n n n n F D S k K = =(15) That is, when the flight altitude changes, the pull required for the level flight of the combination is only related to the lift-drag ratio, which mainly depends on the flight speed and angle of attack [18], and is not sensitive to the altitude change, so the demand for pull at different altitudes is basically unchanged. Using the model of Section 3, since the influence of the fuselage is not considered, the induced resistance is used to estimate the resistance of the whole aircraft. The induced resistance of the flight unit is calculated according to the following formula. 2 1 1 1 1 1 1 2 i Di D V C S ? =(16) In the formula, the lower corner mark 1 represents the flight unit, and the induced resistance of the flight unit is calculated, D i1 =14.91N. The induced resistance of the small fixed wing UAV accounts for 25% ~ 40% of the total resistance. Here, the waste resistance of the flight unit is calculated according to the induced resistance accounting for 30% of the total resistance, D p1 =41.74N. Similarly, the induced resistance of the combination, D in , can be obtained, while the waste resistance of the combination, D pn , can be calculated using the following formula. n pn p V D nD V =(17) Then calculate the tensile demand of a single propeller through the following formula. in pn n D D F n + =(18) The relationship between tensile force F n and n is shown in Figure 8. It can be seen that the required tension F n decreases with the increase of n, and the change rate decreases. When estimating the resistance according to Section 4.1 using the model of section 3, for a motor drive, the 1.5:1 reduction mechanism and 16 × 8 APC propeller is calculated under 32 V power supply. The maximum pull of single propeller at different altitudes is shown in Table 3, and Figure 9 is drawn accordingly. The fitting formula of the curve is obtained. The relationship between n and the ceiling is shown in Figure 10. It can be seen that increasing the number of flight units can effectively improve the ceiling of the aircraft. The ceiling increases with the increase of n, while the change rate decreases. Assuming a 128 Ah 32 V power supply is used and the combination flies horizontally at an altitude of 2 km, the working current and flight time of the combination with different n are calculated as shown in Table 5. The relationship between n and flight time is shown in Figure 11. # V. Results and Discussion Conclusion With the rise of cluster UAVs, wingtip connection technology has been paid attention again. Wingtip connection can make the flight unit fly higher, farther and more durable. Firstly, the theoretical analysis of wingtip connection is carried out, and two conclusions are obtained through the theoretical analysis. First, the liftdrag ratio of the combination K n is proportional to n ne , and the Oswald efficiency factor e n decreases with the increase of n; The second is that the cruise speed of the combination V n is inversely proportional to Ln C . Then, the correctness of the above two conclusions is verified by simulation. Based on this, the influence of wingtip connection on endurance performance is analysed, and the following conclusions are obtained. First, the flight altitude has little effect on the tensile demand, but tensile demand decreases with the increase of n, while the change rate decreases; The second is that the ceiling increases with the increase of n, while the rate of change decreases; Third, the endurance time increases with the increase of n, while the change rate decreases. Theoretically, increasing the number of units n can continuously improve the flight performance of aircraft. # c) Wingtip Connection can Increase Endurance Time However, as the number of assembly units increases, the ceiling and endurance benefits brought by connection gradually decrease. In addition, largescale wingtip connections have high requirements for control systems and connecting mechanisms. Too many units will lead to more complex control and connecting mechanisms with greater loads. Therefore, the number of wingtip connections is not the more the better. 1![Figure 1: Cluster UAVs with wingtip connection](image-2.png "Figure 1 :") 2![Figure 2: Relationship between resistance and speed](image-3.png "Figure 2 :D") 3![Figure 3: Simulation Model of flight unit](image-4.png "Figure 3 :") 4![Figure 4: Variation of aerodynamic parameters vary with respect to the number of units in the assembly. (a) Variation of lift coefficient, resistance coefficient, torque coefficient and lift-drag ratio with respect to angle of attack. (b) Variation of lift-drag ratio with respect to angle of attack See Table1for the sorted data, and draw Figure5for the simulation lift-drag ratio at 0° and the calculated theoretical lift-drag ratio. It can be seen that the law of lift-drag ratio approximately meets n times,](image-5.png "Figure 4 :") 5![Figure 5: Relationship between lift-drag ratio and 1n](image-6.png "Figure 5 :") 6![Figure 6: Relationship between lift coefficient, flight speed and n](image-7.png "Figure 6 :D") 7![Figure 7: Relationship between C Ln ?V n 2 and n IV. Effect of Wingtip Connection on Cruise Performance](image-8.png "Figure 7 :") 8![Figure 8: Relationship and trend between tensile force F n and n](image-9.png "Figure 8 :") 9![Figure 9: Relationship between altitude and maximum pulling force of single propeller](image-10.png "Figure 9 :") 10![Figure 10: Relationship between n and the ceiling](image-11.png "Figure 10 :D") 11![Figure 11: Relationship between n and flight time It can be seen that increasing the number of flight units can effectively improve the flight time. The endurance time increases with the increase of n, while the change rate decreases.](image-12.png "Figure 11 :") ![](image-13.png "") 1bYear 20235Volume Xx XIII Issue I V ersion I( )n 1 2 3 4 5 6 7Simulated lift-drag ratio 24.161 34.397 41.050 45.870 50.199 53.539 56.095Theoretical lift-drag ratio 24.161 34.164 41.847 48.322 54.024 59.170 63.930Speed[m/s] 42.00 39.38 38.40 37.88 37.56 37.34 37.18Lift coefficient 0.553 0.629 0.661 0.679 0.691 0.699 0.705C Ln ?V n 2 975.49 975.44 974.68 974.29 974.60 974.56 974.83Global Journal of Researches in Engineering© 2023 Global Journ als D Study on the Influence of UAV Wingtip Connection on Ceiling and Endurance 2nTensile demand of single propeller F n [N]156.64245.59341.48439.15537.81636.78736.05 3Year 20238Volume Xx XIII Issue I V ersion IAltitude b) Wingtip Connection can Raise the Ceiling Maximum pulling forceAltitudeMaximum pulling forceH[km]F max [N]H [km]F max [N]Global Journal of Researches in Engineering0 1 2 3 4 5 660.22 58.61 56.91 55.10 53.15 51.05 48.767 8 9 10 11 1246.30 43.61 40.65 37.43 33.83 29.84© 2023 Global Journ als 44max F0.1037 = -H2-1.0283H+61.071(19)nTensile demand of single propeller F n [N]Ceiling H max [km]156.643.25245.598.22341.489.65439.1510.40537.8110.82636.7811.13736.0511.35 5nCurrent [A]Flight time [min]Endurance improvement129.93256.60.00%221.78352.737.45%318.97404.857.76%417.45440.271.55%516.59463.080.44%615.93482.087.83%715.48496.293.37% ## Acknowledgments The authors are grateful for the financial support from the China Aviation Science Foundation Program (201829U8007). They would also like to thank the Aeromodelling Team of RFUE for providing support. ## Data Availability The aerodynamic parameters in this paper are calculated by xflr5 v6.48, and the file has been uploaded. The file name is model.xfl. ## Conflicts of Interest The authors declare that there is no conflict of interest regarding the publication of this paper. * Miniature fixed wing UAV swarms: review and outlook XWang ZLiu YCong Acta Aeronautica et Astronautica Sinica 41 4 2020 * Effect of wingtip connection on the energy and flight endurance performance of solar aircraft MWu ZShi TXiao Aerospace Science and Technology 108 106404 2021 * Dangerous experiments: wingtip coupling at 15,000 feet CAnderson Flight Journal 5 6 2000 * Compound aircraft transport study: wingtip-docking compared to formation flight SAMagill 2002 Blacksburg, Virginia, USA Virginia Polytechnic Institute and State University * Meta aircraft flight dynamics and controls CMontalvo 2014 Atlanta, Georgia, US Georgia Institute of Technology * Meta aircraft flight dynamics CMontalvo MCostello Journal of Aircraft 52 1 2015 * Meta aircraft connection dynamics CMontalvo MCostello AIAA guidance, navigation, and control conference Minneapolis, Minnesota, USA 2012 * Experimental control of two connected fixed wing aircraft CCollin CMontalvo Aerospace 113 5 2018 * Takeoff and Landing of a Wing-Tip-Connected Meta Aircraft with Feedback Control MCobar CMontalvo Journal of Aircraft 58 4 2021 * Dynamics and control of in-flight wing tip docking JRCooper PMRothhaar Journal of Guidance Control & Dynamics 41 11 2018 * ELAHAelastic aircraft for high altitudes concept and current development state of an unconventional stratospheric UAV SWlach GBalmer MHermann 23rd ESA Symposium on European Rocket and Balloon Programmes and Related Research Visby, Sweden 2017 * Closedloop flight tests with an unmanned experimental multi-body aircraft AKöthe ABehrens AHamann 17th International Forum on Aeroelasticity and Structural Dynamics Como, Italy 2017 * Overview of UAV link wing system technology RZhou WZhou DWang the 5th Joint Conference of Aerospace Propulsion (JCAP) Technical Conference Nanjing 2020 * Flexible modular swarming UAV: innovative, opportunities, and technical challenges YYang ZZhang PYing Flight Dynamics 39 2 2021 * Research progress of combined fixed-wing unmanned aerial vehicle based on wingtip chained WZhou PMa ZGuo Acta Aeronautica et Astronautica Sinica 43 9 2022.2022 * Flight mechanical analysis and test of unmanned multi-body aircraft CAn CXie YMeng International Forum on Aeroelasticity and Structural Dynamics Savannah, Georgia, USA 2019 * Induced drag of high-aspect ratio wings TChen JKatz 42nd AIAA Aerospace Sciences Meeting and Exhibit Reno, Nevada, USA 2004 * ZFang WChen SZhang Flight dynamics of aircraft Beijing Beihang University Press 2005