Aerodynamic Interference on Trim Characteristics of Quad-Tiltrotor Aircraft

: The aerodynamic interference between the different components of quad-tiltrotor (QTR) aircraft were considered to analyze its inﬂuence on trim characteristics. A comprehensive method with the ﬁxed-wake model was developed for multiple aerodynamic interactions, improving the accuracy of the ﬂight dynamics analysis. Additionally, a more general control strategy was developed to tackle the redundant control issue of the QTR, improving its control efﬁciency by coordinating the authority relationship of various control surfaces across the ﬂight range. Then, the trim features were calculated in the helicopter mode, conversion mode


Introduction
The tiltrotor aircraft has become the research focus in the helicopter industry due to its capability to combine the operational advantages of the fixed-wing airplane and the helicopter, making it capable of both hover and high-speed flight [1,2].However, continuous research for decades has already brought the tiltrotor aircraft to its maximum payload and performance limitations.Therefore, with the aim of further improving its loading capacity and performance requirement, the concept of the quad tiltrotor (QTR) was proposed by Bell Helicopter company [3].The QTR configuration adopts two paralleling nacelle systems, allowing its maximum takeoff weight to be as much as 70,000 kg, similar to the C-130J aircraft.
During the flight of the quad-tiltrotor aircraft, there are significant aerodynamic interactions that would influence its flight dynamic characteristics, including interference inside the multirotor system and the interaction between rotors and wings.Many researchers have implemented different methods to determine the aerodynamic interference, which is usually based on wind tunnel experiments [4][5][6][7][8][9] and CFD calculations [10][11][12].The wind tunnel experiments suffer from extensive time and financial cost, and the CFD method requires high computational cost, making it unsuitable for QTR flight dynamic modelling and analysis.On the other hand, the empirical factor embedded fixed-wake method provides an alternative, more efficient method to quantify the aerodynamic interference of the QTR, which has the capability of synthesizing multiple interactions and has been widely used in the twin-tiltrotor aircraft [13] and tandem helicopter [14,15].As a result, a more accurate aerodynamic model can be established based on this method to investigate the flight dynamics characteristics of QTR.
Aerospace 2022, 9, 262 2 of 21 Furthermore, the control strategy is another critical issue for the QTR aircraft.Reasonable control allocation is the basis for the QTR to convert to safely during flight.The control allocation of the twin-tiltrotor can be directly determined by the tilt angle.The hybrid control equations in the longitudinal plane are established in Reference [16].The attitude angle and the pilot workload are taken into the performance index to obtain the controlled displacement of the stick by using the optimal control method.However, the control redundancy of a quad-tiltrotor aircraft is much more complicated than that of the twin tiltrotor, and its control power is more sensitive to the aerodynamic characteristics in different flight ranges.In Reference [17], the control channels are decoupled to control surfaces, and the speed and height controllers are designed to determine the distribution coefficients of each surface.This method needs to be combined with the control system, which is not suitable for the analysis of flight dynamics.
In light of the preceding discussion, a more accurate aerodynamic model was established, which takes the interference of rotor-rotor and rotor-wing into account based on the synthesis of the empirical factor and fixed-wake method.The developed method is more appropriate for the analysis of flight dynamics.Aiming at the control redundancy problem, a general control strategy for the QTR aircraft is proposed to ensure that the flight dynamics characteristics and safety are satisfied across the flight range.

Flight Dynamics Model and Validation
The flight dynamics model is composed of the aerodynamic models of each component, including the rotor, wing, fuselage, and vertical fin.In this article, the model was developed based on the QTR configuration as shown in Figure 1, which was designed and manufactured by the laboratory at Nanjing University of Aeronautics and Astronautics.It can be seen that the rotors are installed at both ends of the front and rear wings, and the rotors can tilt with the nacelle around the axis.In the four rotors, the rotational directions of the diagonal rotors are the same, while the adjacent is opposite.The rotors are numbered 1, 2, 3, 4, and rotor 1 is right-handed.The flaperons are adopted on each wing surface, and they can be activated collaboratively and differentially to represent the function of the aileron and the elevator.The rudder is installed on the vertical fin to provide heading control.The general parameters of the quad tiltrotor can be seen in Table 1.
and has been widely used in the twin-tiltrotor aircraft [13] and tandem helicopter [14,15].As a result, a more accurate aerodynamic model can be established based on this method to investigate the flight dynamics characteristics of QTR.
Furthermore, the control strategy is another critical issue for the QTR aircraft.Reasonable control allocation is the basis for the QTR to convert to safely during flight.The control allocation of the twin-tiltrotor can be directly determined by the tilt angle.The hybrid control equations in the longitudinal plane are established in Reference [16].The attitude angle and the pilot workload are taken into the performance index to obtain the controlled displacement of the stick by using the optimal control method.However, the control redundancy of a quad-tiltrotor aircraft is much more complicated than that of the twin tiltrotor, and its control power is more sensitive to the aerodynamic characteristics in different flight ranges.In Reference [17], the control channels are decoupled to control surfaces, and the speed and height controllers are designed to determine the distribution coefficients of each surface.This method needs to be combined with the control system, which is not suitable for the analysis of flight dynamics.
In light of the preceding discussion, a more accurate aerodynamic model was established, which takes the interference of rotor-rotor and rotor-wing into account based on the synthesis of the empirical factor and fixed-wake method.The developed method is more appropriate for the analysis of flight dynamics.Aiming at the control redundancy problem, a general control strategy for the QTR aircraft is proposed to ensure that the flight dynamics characteristics and safety are satisfied across the flight range.

Flight Dynamics Model and Validation
The flight dynamics model is composed of the aerodynamic models of each component, including the rotor, wing, fuselage, and vertical fin.In this article, the model was developed based on the QTR configuration as shown in Figure 1, which was designed and manufactured by the laboratory at Nanjing University of Aeronautics and Astronautics.It can be seen that the rotors are installed at both ends of the front and rear wings, and the rotors can tilt with the nacelle around the axis.In the four rotors, the rotational directions of the diagonal rotors are the same, while the adjacent is opposite.The rotors are numbered 1, 2, 3, 4, and rotor 1 is right-handed.The flaperons are adopted on each wing surface, and they can be activated collaboratively and differentially to represent the function of the aileron and the elevator.The rudder is installed on the vertical fin to provide heading control.The general parameters of the quad tiltrotor can be seen in Table 1.In the flight dynamics modeling, the definition of the coordinate systems is shown in Figure 2. In the flight dynamics modeling, the definition of the coordinate systems is shown in Figure 2. where O is the center of gravity of the QTR; V is the inflow velocity; ODXDYDZD is the inertial coordinate system; OXYZ is the body coordinate system; OWXWYWZW is the wind axis; OSXSYSZS is the rotor fixed axis coordinate; OSWXSWYSWZSW is the rotor wind coordinate, in which OSWXSW is the projection of the velocity V onto the OSWXSWYSW plane; βsw, the side-slip angle of the rotor is the angle between OSXS axis and OSWXSW axis.

Rotor Model
The rotor aerodynamics model was developed considering the rotor dynamics characteristics, the induced inflow, and the aerodynamics interaction between the multirotor system.The flapping characteristics of the rotor blade were modeled based on Reference [18].The modified Ferguson model [13] was used to represent the dimensionless induced velocity, as shown in the following equation: where C = CTW/(2B 2 ), B is the blade tip-root loss factor; CTW is the coefficient of the rotor thrust; μ is the advance ratio of the rotor; λ is the inflow ratio of the rotor, which is expressed by Equation (2):

Rotor Model
The rotor aerodynamics model was developed considering the rotor dynamics characteristics, the induced inflow, and the aerodynamics interaction between the multirotor system.The flapping characteristics of the rotor blade were modeled based on Reference [18].The modified Ferguson model [13] was used to represent the dimensionless induced velocity, as shown in the following equation: where C = C TW /(2B 2 ), B is the blade tip-root loss factor; C TW is the coefficient of the rotor thrust; µ is the advance ratio of the rotor; λ is the inflow ratio of the rotor, which is expressed by Equation (2): where w S is the hub vertical velocity in the rotor fixed axis coordinate; Ω is the rotational speed of the rotor; R is the rotor radius.
The blade flapping motion is described in the form of first-order Fourier series, and then the derivation of the tip-path plane dynamic equation can be calculated by the following equations [18].
where a 0 , a 1 , b 1 are the flapping coefficients; D R is the rotor damping matrix; K R is the rotor rigidity matrix; f R is the rotor incentive vector.
According to the induced velocity of the rotor and the flapping motion of the blade, the diagram of solving the rotor model is shown in Figure 3.The details of the procedure can be obtained in the References [18,19].
where wS is the hub vertical velocity in the rotor fixed axis coordinate; Ω is the rotational speed of the rotor; R is the rotor radius.
The blade flapping motion is described in the form of first-order Fourier series, and then the derivation of the tip-path plane dynamic equation can be calculated by the following equations [18].
where a0, a1, b1 are the flapping coefficients; DR is the rotor damping matrix; KR is the rotor rigidity matrix; fR is the rotor incentive vector.
According to the induced velocity of the rotor and the flapping motion of the blade, the diagram of solving the rotor model is shown in Figure 3.The details of the procedure can be obtained in the References [18,19].Then, the corresponding forces and moments acting on the center of gravity of the aircraft in the body coordinate system can be calculated using Equations ( 4) and ( 5), respectively.
where Δ denotes the rotation direction of the rotor, in which Δ = 1 represents right-handed and Δ = −1 represents left-handed; the definition of βsw is shown in Figure 2; ρ is the air density; CHW, CYW, CTW are the rotor coefficients of backward force, sideward force, and thrust respectively, whose calculation formulas can be seen in Reference [18].Q is the rotor shaft torque; LW is the hub moment in the rolling channel; MW is the hub moment in the pitching channel; (xH, yH, zH) T is the vector of the rotor hub relative to the center of gravity.The aerodynamic interference between the multi-rotor system could significantly alter the flight dynamics characteristics across the flight range, which should be considered in the rotor modelling process.Based on References [13,14], this interference can be divided into two aspects, namely, the lateral interaction and the longitudinal interaction.The formation of the lateral interaction shares many similarities with the twin-rotor interaction in the tiltrotor aircraft, and the longitudinal interaction is similar to the rotor interaction in the tandem helicopter.The corresponding interaction between different rotors in the QTR configuration is shown in Table 2.Then, the corresponding forces and moments acting on the center of gravity of the aircraft in the body coordinate system can be calculated using Equations ( 4) and (5), respectively.
where ∆ denotes the rotation direction of the rotor, in which ∆ = 1 represents right-handed and ∆ = −1 represents left-handed; the definition of β sw is shown in Figure 2; ρ is the air density; C HW , C YW , C TW are the rotor coefficients of backward force, sideward force, and thrust respectively, whose calculation formulas can be seen in Reference [18].Q is the rotor shaft torque; L W is the hub moment in the rolling channel; M W is the hub moment in the pitching channel; (x H , y H , z H ) T is the vector of the rotor hub relative to the center of gravity.The aerodynamic interference between the multi-rotor system could significantly alter the flight dynamics characteristics across the flight range, which should be considered in the rotor modelling process.Based on References [13,14], this interference can be divided into two aspects, namely, the lateral interaction and the longitudinal interaction.The formation of the lateral interaction shares many similarities with the twin-rotor interaction in the tiltrotor aircraft, and the longitudinal interaction is similar to the rotor interaction in the tandem helicopter.The corresponding interaction between different rotors in the QTR configuration is shown in Table 2.In Table 2, the long dash "-" means ignoring the interaction of these two rotors.It can be seen that the lateral rotor interaction includes the interference between rotor 1 and rotor 2, and interference between rotor 3 and rotor 4. The longitudinal rotor interaction contains the aerodynamic interference between rotor 1 and rotor 4 and interference between rotor 2 and rotor 3.

Lateral Interaction
The lateral interaction would add additional induced velocity on the rotor disc, which can be calculated using the lateral interaction factor, X S , which is shown in Equation (6).
where X SF is the side-by-side-induced velocity coefficient in the forward flight of the aircraft, determined by the rotor advance ratio µ.X SS is the sideward flight-induced velocity coefficient, decided by the lateral velocity advance ratio.The calculation equations can be written as follows.
if v H,L + v H,R ≥ 0: otherwise: where µ L , µ R are the advanced ratios of the left and right rotors respectively; v H,L , v H,R are the lateral velocity at the hub of the left and right rotors respectively.Therefore, the velocity increment provided by the lateral interaction is calculated based on Equation ( 10) where v i can be obtained from Equation (1).

Longitudinal Interaction
The longitudinal interaction is more complicated than the lateral one as the nacelle incidence angle would play a major effect on the interaction scheme.Therefore, the relevant calculation method should coordinate the interaction phenomena in different modes (helicopter mode, conversion mode, and airplane mode).
For helicopter mode, front rotor interference factor X H,F and rear rotor interference factor X H,B can be written as follows, respectively [14]: where χ F and χ B represent the wake angles of the front and rear rotors, respectively.β SW is the angle of the side-slip of the rotor.η FR is the relative position factor of the front-rear rotor interference, determined by the projection relationship between front and rear rotor, as shown in Equation (13).
where l RRL is the lateral distance between the centre of the front and rear rotors, S FB is the area of the shaded part, as shown in Figure 4a.Thus, the corresponding induced velocity increments can be written as: where v i,B and v i,F represent induced velocities of front and rear rotors, respectively.
Aerospace 2022, 9, x FOR PEER REVIEW 6 of 21 where χF and χB represent the wake angles of the front and rear rotors, respectively.βSW is the angle of the side-slip of the rotor.ηFR is the relative position factor of the front-rear rotor interference, determined by the projection relationship between front and rear rotor, as shown in Equation ( 13).
where lRRL is the lateral distance between the centre of the front and rear rotors, SFB is the area of the shaded part, as shown in Figure 4a.Thus, the corresponding induced velocity increments can be written as: where vi,B and vi,F represent induced velocities of front and rear rotors, respectively.In airplane mode, the rear rotor cannot introduce a significant influence on the aerodynamics of the front rotor, and only the rear rotor can be impacted at this flight range.Therefore, based on the projection relationship in the airplane mode, the induced velocity increment of the rear rotor can be calculated using Equation (15).

,, 2
where RFR is the influential area radius of the front rotor wake at the rear rotor disc, as shown in Figure 4b, and it can be calculated as follows.
( ) where CH, CS, CT are the coefficient of the backward force, sideward force, and thrust of the rotor system, respectively.dn is the longitudinal distance between the centre of the front and rear rotors.In airplane mode, the rear rotor cannot introduce a significant influence on the aerodynamics of the front rotor, and only the rear rotor can be impacted at this flight range.Therefore, based on the projection relationship in the airplane mode, the induced velocity increment of the rear rotor can be calculated using Equation (15).
where R FR is the influential area radius of the front rotor wake at the rear rotor disc, as shown in Figure 4b, and it can be calculated as follows.
where C H , C S , C T are the coefficient of the backward force, sideward force, and thrust of the rotor system, respectively.d n is the longitudinal distance between the centre of the front and rear rotors.
In the conversion mode, the interaction can be calculated using the fitting method to combine the results in the helicopter mode and airplane mode based on the tilt angle β N , which can be expressed as:

Wing Model
In the flight of the quad-tiltrotor aircraft, the downwash of the rotor directly impacts the wing's surface, which will generate an additional aerodynamic load.Reference [20] has indicated that the download to the wing generated from the rotor can reach 10~15% of the rotor's thrust of a tiltrotor aircraft in the hover state.For the quad-tiltrotor aircraft, this downwash effect can be summarized as two types, including the rotor downwash effect on its connected wing (Type 1), and the downwash interference on its unattached wing (Type 2).

Rotor-Wing Interaction: Type 1
To quantify this rotor-wing interaction, the wing can be divided into two parts: the area affected by the rotor wake, i.e., the slipstream area; the area unaffected by the rotor wake, called the freestream area.The resultant aerodynamic force of the wing is the superposition of the two aerodynamic forces.The areas of slipstream and freestream can be calculated as follows [21].
where the subscript WS represents the slipstream area.WF represents the freestream area.µ max is the maximum advanced ratio of the rotor while the wake phasing out the upper connected wing.R i is the influential area radius of the rotor wake at the connected wing.c is the mean chord of the wing.S W is the total area of the wing.The flow velocity of the slipstream can be expressed as: (20) where r WWS is the aerodynamic centre of the slipstream area.

Rotor-Wing Interaction: Type 2
The schematic diagram of the rotor downwash interference on its unattached wing is shown in Figure 5.The cross effect of the rotor-wing interaction (Rotor 1 downwash on the wing connecting Rotor 3) is ignored due to it being minimal compared to other influences.The interaction of the rear rotor with the front wing only exists when the aircraft flies backwards, and its calculation method is the same as the method shown below.The side view of the interference is illustrated in Figure 5a, in which αRW,min is the critical angle of the front rotor that starts to interfere with the rear wing; αRW,max is the angle that the rear wing starts to be entirely immersed in the wake of the front rotor; hWW is the height difference of the front-rear wings; lWW is the longitudinal distance between the front and rear wings; dh is the distance between the rotor hub centre and the tilting axis.Therefore, the relevant angles can be calculated using the following equations.
The top view of the interference is illustrated in Figure 5b, representing the maximum interference area of the front rotor to the rear wing.When βN + χ ≥ αRW,min, the wake of the front rotor alters the aerodynamics of the rear wing, the corresponding interference area SRW can be calculated by the following equations.
where lRW is the lateral distance from the outer end of the rear wing to the hub of the front rotor; RRW is the contraction radius of the front rotor wake at the rear wing.The correction of the flow velocity affected by this interference can be written as: The dynamic pressure of the slipstream qWS and the freestream qWF can be calculated considering this effect.Then, the resultant forces and moments of the wing component can be expressed as: The side view of the interference is illustrated in Figure 5a, in which α RW,min is the critical angle of the front rotor that starts to interfere with the rear wing; α RW,max is the angle that the rear wing starts to be entirely immersed in the wake of the front rotor; h WW is the height difference of the front-rear wings; l WW is the longitudinal distance between the front and rear wings; d h is the distance between the rotor hub centre and the tilting axis.Therefore, the relevant angles can be calculated using the following equations.
The top view of the interference is illustrated in Figure 5b, representing the maximum interference area of the front rotor to the rear wing.When β N + χ ≥ α RW,min , the wake of the front rotor alters the aerodynamics of the rear wing, the corresponding interference area S RW can be calculated by the following equations.
where l RW is the lateral distance from the outer end of the rear wing to the hub of the front rotor; R RW is the contraction radius of the front rotor wake at the rear wing.The correction of the flow velocity affected by this interference can be written as: The dynamic pressure of the slipstream q WS and the freestream q WF can be calculated considering this effect.Then, the resultant forces and moments of the wing component can be expressed as: Aerospace 2022, 9, 262 where L W is the lift of the wing; D W is the drag of the wing; M W is the pitching moment of the wing; C L , C D , C M are the lift, drag, and pitching moment coefficients, respectively, which are related to the Mach number and configuration of the flaperons.Therefore, the forces (X W , Y W , Z W ) T and moments (L W , M W , N W ) T produced by the wing in the body coordinate system can be obtained through coordination transformation.

Model of Other Components
Other aerodynamic components in the model include the fuselage and the vertical fin.The fuselage aerodynamics was obtained by CFD using the fitting method.The detailed introduction of the method is shown in reference [22] and was verified by the experiment.In brief, the ansys-icem is used to mesh grid and the ansys-fluent is used for calculation.The number of boundary layers on the fuselage is 10.The normal height ratio of meshes along the wall is 1.15, and the total meshes are 4 million.The k-ω shear stress transport (SST) turbulence model was employed, and the SIMPLE algorithm was used to solve the model.
Thus, the resultant aerodynamic force coefficients C LF , C DF , C YF , C RF , C MF , C NF , related to the angle of attack and side-slip of the aircraft, can be obtained.Figure 6 indicates the numerical simulation of the fuselage aerodynamics.
where  ̃ is the lift of the wing; DW is the drag of the wing; MW is the pitching moment of the wing; CL, CD, CM are the lift, drag, and pitching moment coefficients, respectively, which are related to the Mach number and configuration of the flaperons.Therefore, the forces (XW, YW, ZW) T and moments (LW, MW, NW) T produced by the wing in the body coordinate system can be obtained through coordination transformation.

Model of Other Components
Other aerodynamic components in the model include the fuselage and the vertical fin.The fuselage aerodynamics was obtained by CFD using the fitting method.The detailed introduction of the method is shown in reference [22] and was verified by the experiment.In brief, the ansys-icem is used to mesh grid and the ansys-fluent is used for calculation.The number of boundary layers on the fuselage is 10.The normal height ratio of meshes along the wall is 1.15, and the total meshes are 4 million.The k-ω shear stress transport (SST) turbulence model was employed, and the SIMPLE algorithm was used to solve the model.
Thus, the resultant aerodynamic force coefficients CLF, CDF, CYF, CRF, CMF, CNF, related to the angle of attack and side-slip of the aircraft, can be obtained.Figure 6  At last, the aerodynamics of the fuselage in the body coordinate can be expressed as follows.At last, the aerodynamics of the fuselage in the body coordinate can be expressed as follows.
where q F∞ is the local dynamic pressure; S F is the flat plate area of the fuselage; l F is the characteristic length of the fuselage; T B W is the transfer matrix from the wind coordinate to the body coordinate; (x F , y F , z F ) T is the distance from the aerodynamic centre of the fuselage to the centre of gravity.
The vertical fin aerodynamics are decided by the lift coefficient C LVT and the drag coefficient C DVT .Thereby, the resultant aerodynamic force can be obtained in the body coordinate.
where q VT∞ is the local dynamic pressure; S VT is the area of the vertical fin; i VT is the incidence of the vertical fin; β VT is the side-slip angle of the vertical fin; α VT is the angle of attack.

Integrated Equations of the Vehicle
The whole model has 16 degrees of freedom, and its general nonlinear equations of motion take the form [23]: in which x is the state vector, x = [u, v, w, p, q, r, ϕ, θ, ψ, v 0,LF , v 0,RF , v 0,RB , v 0,LB , where the subscript LF represents rotor 1, RF represents rotor 2, RB represents rotor 3, LB represents rotor 4; u, v, w, p, q, and r are the rigid body translational velocity components and angular rates; ϕ, θ, and ψ are the rigid body Euler angles; v 0 is the induced velocity.δ S is the control input vector, where δ col is the collective stick displacement; δ lat is the lateral stick displacement; δ lon is the longitudinal stick displacement; δ ped is the pedal displacement.The analysis of δ S is detailed in chapter 3.

Validation
In order to validate the accuracy of the model established in the paper, the experiment of the thrust of the isolated rotor was carried out in the wind tunnel, whose test section size is 3.2 m × 2.4 m and the maximum wind speed is 50 m/s.The equipment of the isolated rotor test is shown in Figure 7a,b.The rotor test involves two aspects: one is the required power changing with the rotor thrust; the other is the rotor thrust changing with the flight velocity.The collective is set to 6 • , 8 • , 10 • , 12 • , and 14 • for the first case, the wind speed is 0 m/s.The collective is fixed at 10 • for the second case, the wind speed changes from 0 m/s to 16 m/s.To avoid the ground effect, the height from the ground of the rotor should not be less than 1.2 times its diameter.
Additionally, the rotor-wing interaction was implemented in the article.However, due to the restriction of the wind tunnel facilities, the interaction test only included rear rotors to the rear wing in the hover state, which can be seen in Figure 7c.The wing's download caused by the rotor thrust could be measured with different collective inputs that involve 6 • , 8 • , 10 • , 12 • , and 14 • .In these two experiments, the rotor rotational speeds were set to be 2100 rpm.The comparisons between test and calculation are shown in Figure 8.Additionally, the rotor-wing interaction was implemented in the article.However, due to the restriction of the wind tunnel facilities, the interaction test only included rear rotors to the rear wing in the hover state, which can be seen in Figure 7c.The wing's download caused by the rotor thrust could be measured with different collective inputs that involve 6 °, 8 °, 10 °, 12 °, and 14 °.In these two experiments, the rotor rotational speeds were set to be 2100 rpm.The comparisons between test and calculation are shown in   Additionally, the rotor-wing interaction was implemented in the article.However, due to the restriction of the wind tunnel facilities, the interaction test only included rear rotors to the rear wing in the hover state, which can be seen in Figure 7c.The wing's download caused by the rotor thrust could be measured with different collective inputs that involve 6 °, 8 °, 10 °, 12 °, and 14 °.In these two experiments, the rotor rotational speeds were set to be 2100 rpm.According to Figure 8, the results from the developed models were in good agreement with the experiment results, indicating the accuracy of the modeling method.The download of the wing accounts for approximately 12~13% of the rotor thrust in the hover state, consistent with that described in Reference [20].To further improve the accuracy of the model, especially the aerodynamic interference during the maneuvring flight, the aircraft flight test was the research focus in the next step.

Control Strategy
The control strategy of the QTR is more complicated than other rotor-powered vehicles due to its redundant control characteristics.The control inputs of the QTR configuration are shown in Table 3.As shown in Table 3, there are 17 control surfaces in the QTR configurations, indicating that this is an underdetermined system with more unknown variables than the number of equations in trimming.
When the quad-tiltrotor aircraft is flying in helicopter mode, the control efficiencies of the flaperon and rudder are relatively small.Therefore, the controllability is mainly provided by the rotor system.When the QTR is in conversion and airplane modes, the control efficiencies of the flaperon and rudder increase with tilting of the nacelle, and the rotor controllability decreases.Furthermore, the change of the tilt angle may lead to additional control coupling problems.Therefore, a mixed control strategy was proposed for the QTR configuration, which is shown below.According to Figure 8, the results from the developed models were in good agreement with the experiment results, indicating the accuracy of the modeling method.The download of the wing accounts for approximately 12~13% of the rotor thrust in the hover state, consistent with that described in Reference [20].To further improve the accuracy of the model, especially the aerodynamic interference during the maneuvring flight, the aircraft flight test was the research focus in the next step.

Control Strategy
The control strategy of the QTR is more complicated than other rotor-powered vehicles due to its redundant control characteristics.The control inputs of the QTR configuration are shown in Table 3.As shown in Table 3, there are 17 control surfaces in the QTR configurations, indicating that this is an underdetermined system with more unknown variables than the number of equations in trimming.
When the quad-tiltrotor aircraft is flying in helicopter mode, the control efficiencies of the flaperon and rudder are relatively small.Therefore, the controllability is mainly provided by the rotor system.When the QTR is in conversion and airplane modes, the control efficiencies of the flaperon and rudder increase with tilting of the nacelle, and the rotor controllability decreases.Furthermore, the change of the tilt angle may lead to additional control coupling problems.Therefore, a mixed control strategy was proposed for the QTR configuration, which is shown below.

Strategy in the Helicopter Mode
In helicopter mode, the vertical channel is controlled by the combined collective with the four rotors.The lateral channel is achieved through the collective differential of left and right rotors.The longitudinal channel is achieved by the collective differential of front and rear rotors.The yawing channel is performed by the differential lateral cyclic pitch of front and rear rotors.Based on Equation (32), the control input in helicopter mode can be written as follows.
The four control input variables are allocated to each rotor so that the corresponding collective, longitudinal cyclic pitch, and lateral cyclic pitch can be determined.Their relationships are shown in Table 4.
It can be seen, in Table 4, the longitudinal motion is controlled by the collective differential of front and rear rotors, resulting in the longitudinal cyclic invariable in this mode.To unify the control inputs in different flight modes, the displacements of the sticks set as the control inputs, and their variation ranges are from 0 to 1, as indicated as follows.
Therefore, the relationship between the control input δ H in helicopter mode and the displacement of the stick can be written as: where k ver , k lat , k lon , k yaw represent collective, lateral, longitudinal, and heading control coefficients in the helicopter mode, respectively.The constants δ col,n , δ lat,n , δ lon,n , and δ ped,n are the neutral displacement of the stick.

Strategy in the Airplane Mode
In the airplane mode, the throttle control is achieved through the combined collective of the four rotors; the longitudinal and lateral control are conducted by the combined and differential flaperon of the wing; the heading is achieved by the rudder.According to Equation (32), the control input in airplane mode can be expressed as: Therefore, the relationship between the control input δ A in the airplane mode and the stick displacement can be written as follows. ) where k a , k e , k t represent lateral, longitudinal, and heading control coefficients in the airplane mode.
To sum up, the control input in each mode is transformed into a function of the stick displacement so the values of rotor collective, rotor cyclic pitch, flaperon, rudder, and attitude are unique, which can be obtained in the trim calculation.Thereby, the control redundant problem of the quad-tiltrotor aircraft is settled.

Results
The aerodynamic interference influence on the trim characteristics is investigated in this part.The relevant calculation points were selected from the conversion corridor which was calculated based on the method in Reference [24], as shown in Figure 9.
stick displacement can be written as follows.
To sum up, the control input in each mode is transformed into a functi displacement so that the values of rotor collective, rotor cyclic pitch, flapero attitude are unique, which can be obtained in the trim calculation.Thereb redundant problem of the quad-tiltrotor aircraft is settled.

Results
The aerodynamic interference influence on the trim characteristics is i this part.The relevant calculation points were selected from the conve which was calculated based on the method in Reference [24], as shown in F At these trimming points, the vehicle is in equilibrium.In this case, the b forces and moments are zero.Therefore, the body translational accelera rates, and angular accelerations of the vehicle are zero.Additionally, the fla of the rotor maintains a steady state, and Equation (3) can be written as foll Combined with Equation ( 29), the Newton-Raphson algorithm is use equilibrium equations.The trim characteristics can be calculated with and ference according to the procedure diagram shown in Figure 10.At these trimming points, the vehicle is in equilibrium.In this case, the body resultant forces and moments are zero.Therefore, the body translational accelerations, angular rates, and angular accelerations of the vehicle are zero.Additionally, the flapping motion of the rotor maintains a steady state, and Equation ( 3) can be written as follows.
Combined with Equation (29), the Newton-Raphson algorithm is used to solve the equilibrium equations.The trim characteristics can be calculated with and without interference according to the procedure diagram shown in Figure 10.In Figure 10, the trim value includes collective stick, lateral stick, longitudinal stick, roll attitude, and pitch attitude.The red parts indicate the aerodynamic interactions, such as rotor lateral interaction, rotor longitudinal interaction, and rotor-wing interaction.
The trim calculation results are shown in Figure 11, which demonstrates that the aerodynamic interaction affects the trim results in the collective and longitudinal sticks and the pitch attitude.The lateral and heading sticks and roll attitudes are close to zero across the flight range.In Figure 10, the trim value includes collective stick, lateral stick, longitudinal stick, roll attitude, and pitch attitude.The red parts indicate the aerodynamic interactions, such as rotor lateral interaction, rotor longitudinal interaction, and rotor-wing interaction.
The trim calculation results are shown in Figure 11, which demonstrates that the aerodynamic interaction affects the trim results in the collective and longitudinal sticks and the pitch attitude.The lateral and heading sticks and roll attitudes are close to zero across the flight range.In Figure 10, the trim value includes collective stick, lateral stick, longitudinal stick, roll attitude, and pitch attitude.The red parts indicate the aerodynamic interactions, such as rotor lateral interaction, rotor longitudinal interaction, and rotor-wing interaction.
The trim calculation results are shown in Figure 11, which demonstrates that the aerodynamic interaction affects the trim results in the collective and longitudinal sticks and the pitch attitude.The lateral and heading sticks and roll attitudes are close to zero across the flight range.

Influence on Displacement of Collective Stick
Figure 11a illustrates that the collective stick would decrease at the beginning, and then increases after 23 m/s.The aerodynamic interference would add the collective pitch in the hover and low-speed flight ranges in the helicopter mode.This is due to the additional downwash load provided by the aerodynamic interaction.
To analyze the impact weight of different aerodynamic interactions on the collective stick, the body forces and rotor thrust are calculated, as shown in Figures 12 and 13.As indicated in Figure 12a, the corresponding direction of the vertical force increment is upward in hover.As the velocity accumulates, its effect increases at first and then decreases.The maximum vertical force increment is about 13 N, which accounts for 2.2% of the weight of the vehicle.In Figure 12b, the corresponding vertical force increment caused by the rotor longitudinal interference decreases at the beginning and then increases.This is because, at low speed, the thrust increment of the front rotor is more significant than the decrement of the rear rotor.As the speed reaches a certain point, the thrust increment of the front rotor becomes smaller than that of the rear rotor, which is illustrated in Figure 13b.The maximum variation of its value is 21 N, accounting for 3.6% of the vehicle's weight.Figure 12c shows that the vertical force increment caused by the interference of the rotor-wing reaches its maximum in the hover state, which is 80 N, accounting for 13.6% of the vehicle's weight.As the flight velocity increases, the rotor wake deviates from the connected wing, decreasing the impact force of the wing.In general, the vertical force increment in the helicopter mode decreases rapidly in the helicopter mode low-speed interval.At this time, the interference of the rotor-wing plays a significant role.When the speed reaches a specific point, the vertical force increment does not change significantly, as shown in Figure 12d.At this time, the interference of longitudinal rotor-rotor plays a major role.

Influence on Displacement of Collective Stick
Figure 11a illustrates that the collective stick would decrease at the beginning, and then increases after 23 m/s.The aerodynamic interference would add the collective pitch in the hover and low-speed flight ranges in the helicopter mode.This is due to the additional downwash load provided by the aerodynamic interaction.
To analyze the impact weight of different aerodynamic interactions on the collective stick, the body forces and rotor thrust are calculated, as shown in Figures 12 and 13.As indicated in Figure 12a, the corresponding direction of the vertical force increment is upward in hover.As the velocity accumulates, its effect increases at first and then decreases.The maximum vertical force increment is about 13 N, which accounts for 2.2% of the weight of the vehicle.In Figure 12b, the corresponding vertical force increment caused by the rotor longitudinal interference decreases at the beginning and then increases.This is because, at low speed, the thrust increment of the front rotor is more significant than the decrement of the rear rotor.As the speed reaches a certain point, the thrust increment of the front rotor becomes smaller than that of the rear rotor, which is illustrated in Figure 13b.The maximum variation of its value is 21 N, accounting for 3.6% of the vehicle's weight.Figure 12c shows that the vertical force increment caused by the interference of the rotor-wing reaches its maximum in the hover state, which is 80 N, accounting for 13.6% of the vehicle's weight.As the flight velocity increases, the rotor wake deviates from the connected wing, decreasing the impact force of the wing.In general, the vertical force increment in the helicopter mode decreases rapidly in the helicopter mode low-speed interval.At this time, the interference of the rotor-wing plays a significant role.When the speed reaches a specific point, the vertical force increment does not change significantly, as shown in Figure 12d.At this time, the interference of longitudinal rotor-rotor plays a major role.In the conversion mode, the vertical force increment caused by the interference of lateral rotor-rotor increases (Figure 12a).In contrast, the vertical force increment caused by the interference of longitudinal rotor-rotor decreases (Figure 12b).The magnitude changes rapidly at the low tilting angle and changes smoothly at the large tilting angle, and then becomes zero approaching the aircraft mode.The vertical force increment caused by the interference of the rotor-wing decreases at first and then increases (Figure 12c).The maximum variation value is about 40 N, accounting for 6.8% of the vehicle's weight, and In the conversion mode, the vertical force increment caused by the interference of lateral rotor-rotor increases (Figure 12a).In contrast, the vertical force increment caused by the interference of longitudinal rotor-rotor decreases (Figure 12b).The magnitude changes rapidly at the low tilting angle and changes smoothly at the large tilting angle, and then becomes zero approaching the aircraft mode.The vertical force increment caused by the interference of the rotor-wing decreases at first and then increases (Figure 12c).The maximum variation value is about 40 N, accounting for 6.8% of the vehicle's weight, and In the conversion mode, the vertical force increment caused by the interference of lateral rotor-rotor increases (Figure 12a).In contrast, the vertical force increment caused by the interference of longitudinal rotor-rotor decreases (Figure 12b).The magnitude changes rapidly at the low tilting angle and changes smoothly at the large tilting angle, and then becomes zero approaching the aircraft mode.The vertical force increment caused by the interference of the rotor-wing decreases at first and then increases (Figure 12c).The maximum variation value is about 40 N, accounting for 6.8% of the vehicle's weight, and tends to be zero in the conversion mode with a higher tilt angle.This change is due to the wake of the front rotor affecting the dynamic pressure of the rear wing, which leads to the change in the lift of the rear wing.In the airplane mode, the vertical force increment of the interferences caused by the lateral and longitudinal rotor-rotor are close to zero.In contrast, the interference caused by the rotor-wing will be reduced slightly.
To sum up, the influence of total interferences on the displacement of the collective stick is significant in hover and low-speed forward flight ranges.At the same time, there is little change in the velocity when the forward speed is greater than 10 m/s.

Influence on Longitudinal Trim Characteristic
The longitudinal trim characteristic of the vehicle includes the displacement of longitudinal stick and pitch angle.The aerodynamic influence on the longitudinal trim characteristic is determined by the pitch moment.
To analyze the impact weight of different aerodynamic interactions on the longitudinal stick and pitching attitude, the body moments and rotor thrust were calculated, as shown in Figures 14 and 15.As can be seen in Figure 14a, the pitch moment increment caused by the interference of the lateral rotor-rotor was relatively small, with a maximum value of 0.6 N•m. Figure 14b illustrates that the interference of the longitudinal rotor-rotor only affected the pitching moment.In the helicopter mode, due to the increase of the front rotor thrust and the decrease of the rear rotor thrust, the nose-up pitching moment of the airframe increased continuously, and its maximum value was 36 N•m.In the conversion mode, by tilting forward the rotor, the moment arm of the rotor thrusts to the centre of gravity decreased.In addition, the rotor thrust increment decreased, so did the nose-up pitching moment.In airplane mode, the nose-up moment increased slowly because of the reduction of the rotor thrust.As shown in Figure 14c, there was an additional nose-down pitching moment to the vehicle in the hover state due to the downwash flow of the rotor.Then, the nose-down pitching moment decreased with forward speed increasing.In conversion mode, the vertical force increment of the rear wing decreased and then increased, as shown in Figure 15, which resulted in an oscillating nose-down pitching moment.In airplane mode, the increased lift of the rear wing led to the nose-down moment rising.Overall, when all interferences were taken into account, as shown in Figure 14d, it would provide a nose-down moment in the hover state.The pilot would need to pull the longitudinal stick backwards.In this case, the displacement of the stick is reduced slightly.By increasing the pitching moment, the displacement of the longitudinal stick also needs to be increased in the helicopter mode.In the conversion mode, the pitching moment decreases rapidly while the tilt angle is small.The displacement of the longitudinal stick also decreases.The influence of the aerodynamic interference on the stick is not significant when β N ≥ 15 deg.In the airplane mode, the increase of nose-down pitching moment results in the reduction of the longitudinal stick displacement.
tends to be zero in the conversion mode with a higher tilt angle.This change is due to the wake of the front rotor affecting the dynamic pressure of the rear wing, which leads to the change in the lift of the rear wing.In the airplane mode, the vertical force increment of the interferences caused by the lateral and longitudinal rotor-rotor are close to zero.In contrast, the interference caused by the rotor-wing will be reduced slightly.
To sum up, the influence of total interferences on the displacement of the collective stick is significant in hover and low-speed forward flight ranges.At the same time, there is little change in the velocity when the forward speed is greater than 10 m/s.

Influence on Longitudinal Trim Characteristic
The longitudinal trim characteristic of the vehicle includes the displacement of longitudinal stick and pitch angle.The aerodynamic influence on the longitudinal trim characteristic is determined by the pitch moment.
To analyze the impact weight of different aerodynamic interactions on the longitudinal stick and pitching attitude, the body moments and rotor thrust were calculated, as shown in Figures 14 and 15.As can be seen in Figure 14a, the pitch moment increment caused by the interference of the lateral rotor-rotor was relatively small, with a maximum value of 0.6 N•m. Figure 14b illustrates that the interference of the longitudinal rotor-rotor only affected the pitching moment.In the helicopter mode, due to the increase of the front rotor thrust and the decrease of the rear rotor thrust, the nose-up pitching moment of the airframe increased continuously, and its maximum value was 36 N•m.In the conversion mode, by tilting forward the rotor, the moment arm of the rotor thrusts to the centre of gravity decreased.In addition, the rotor thrust increment decreased, so did the nose-up pitching moment.In airplane mode, the nose-up moment increased slowly because of the reduction of the rotor thrust.As shown in Figure 14c, there was an additional nose-down pitching moment to the vehicle in the hover state due to the downwash flow of the rotor.Then, the nose-down pitching moment decreased with forward speed increasing.In conversion mode, the vertical force increment of the rear wing decreased and then increased, as shown in Figure 15, which resulted in an oscillating nose-down pitching moment.In airplane mode, the increased lift of the rear wing led to the nose-down moment rising.Overall, when all interferences were taken into account, as shown in Figure 14d, it would provide a nose-down moment in the hover state.The pilot would need to pull the longitudinal stick backwards.In this case, the displacement of the stick is reduced slightly.By increasing the pitching moment, the displacement of the longitudinal stick also needs to be increased in the helicopter mode.In the conversion mode, the pitching moment decreases rapidly while the tilt angle is small.The displacement of the longitudinal stick also decreases.The influence of the aerodynamic interference on the stick is not significant when βN ≥ 15 deg.In the airplane mode, the increase of nose-down pitching moment results in the reduction of the longitudinal stick displacement.Figure 11f shows that the pitching attitude decreased in the helicopter mode.The vehicle needs to change the pitching attitude to make gravity create the longitudinal force to balance the aircraft.The additional drag caused by the interference increased at the beginning and then decreased, as shown in Figure 12d.Therefore, the vehicle will first pitch down faster, and then the trend of variation will slow down compared to non-interaction.In the conversion mode, the pitching attitude first increased and decreased later.The reason is that the forward longitudinal force generated by the rotor is more remarkable as the tilt angle increases.The aircraft needs to raise its nose so that the longitudinal component of gravity transforms to the drag to balance the body.When the tilt angle increases to a particular value, the thrust of the rotor and pitching attitude decreases rapidly with velocity increases.Since the interference will increase the drag of the vehicle, the pitching attitude is smaller than the case where the interference is not considered.In the airplane mode, the aircraft is in the state of nose-down, and its pitching attitude keeps decreasing as the speed increases.It is because the wing has an initial incidence.As the dynamic pressure increases, the vehicle needs to reduce the wing angle of attack to balance gravity.Figure 11f shows that the pitching attitude decreased in the helicopter mode.The vehicle needs to change the pitching attitude to make gravity create the longitudinal force to balance the aircraft.The additional drag caused by the interference increased at the beginning and then decreased, as shown in Figure 12d.Therefore, the vehicle will first pitch down faster, and then the trend of variation will slow down compared to non-interaction.In the conversion mode, the pitching attitude first increased and decreased later.The reason is that the forward longitudinal force generated by the rotor is more remarkable as the tilt angle increases.The aircraft needs to raise its nose so that the longitudinal component of gravity transforms to the drag to balance the body.When the tilt angle increases to a particular value, the thrust of the rotor and pitching attitude decreases rapidly with velocity increases.Since the interference will increase the drag of the vehicle, the pitching attitude is smaller than the case where the interference is not considered.In the airplane mode, the aircraft is in the state of nose-down, and its pitching attitude keeps decreasing as the speed increases.It is because the wing has an initial incidence.As the dynamic pressure increases, the vehicle needs to reduce the wing angle of attack to balance gravity.Figure 11f shows that the pitching attitude decreased in the helicopter mode.The vehicle needs to change the pitching attitude to make gravity create the longitudinal force to balance the aircraft.The additional drag caused by the interference increased at the beginning and then decreased, as shown in Figure 12d.Therefore, the vehicle will first pitch down faster, and then the trend of variation will slow down compared to non-interaction.In the conversion mode, the pitching attitude first increased and decreased later.The reason is that the forward longitudinal force generated by the rotor is more remarkable as the tilt angle increases.The aircraft needs to raise its nose so that the longitudinal component of gravity transforms to the drag to balance the body.When the tilt angle increases to a particular value, the thrust of the rotor and pitching attitude decreases rapidly with velocity increases.Since the interference will increase the drag of the vehicle, the pitching attitude is smaller than the case where the interference is not considered.In the airplane mode, the aircraft is in the state of nose-down, and its pitching attitude keeps decreasing as the speed increases.It is because the wing has an initial incidence.As the dynamic pressure increases, the vehicle needs to reduce the wing angle of attack to balance gravity.

Conclusions
In this paper, the calculation methods of aerodynamic interaction in the QTR are developed.Then an integrated modelling approach was developed for corresponding flight dynamics analysis.A mixed control strategy was proposed to combine multiple redundant control inputs.Based on this strategy, the trim characteristics of the QTR were investigated.Additionally, the influence of different types of aerodynamic interference on the trim characteristic was analyzed.The following conclusions can be obtained: 1.
The interference of lateral rotor-rotor has little influence on the body forces and moments.The maximum variation of the vertical force increment is 13 N, which accounts for 2.2% of the vehicle's weight.On the other hand, the interference of the longitudinal rotor-rotor mainly affects the vertical force and pitching moment.The maximum vertical force increment is 21 N, which accounts for 3.6% of the overall weight.In addition, the nose-up pitching moment varies greatly, with a maximum value of 36 N•m; 2.
The interference of the rotor-wing mainly affects the vertical force and pitching moment.The variation of the vertical force increment reaches its maximum in the hover state, accounting for 13.6% of the vehicle's weight.At this time, a nose-down pitching moment of about 5 N•m will be generated.The pitching moment will be changed significantly in the conversion mode, and the corresponding maximum pitching moment is 15 N•m.

3.
The trim values of the longitudinal stick and the pitching angle are affected significantly by rotor-rotor and rotor-wing interference in the helicopter mode and conversion mode with a lower tilt angle.The interaction also affects the longitudinal stick in the airplane mode and the collective stick at low-speed range in the helicopter mode.

Figure 1 .
Figure 1.Configuration of the QTR.(a) Three-dimensional model of QTR.(b) Hover state of QTR flight test.Figure 1. Configuration of the QTR.(a) Three-dimensional model of QTR.(b) Hover state of QTR flight test.

Figure 1 .
Figure 1.Configuration of the QTR.(a) Three-dimensional model of QTR.(b) Hover state of QTR flight test.Figure 1. Configuration of the QTR.(a) Three-dimensional model of QTR.(b) Hover state of QTR flight test.

Figure 2 .
Figure 2. Definition of coordinate systems.whereO is the center of gravity of the QTR; V is the inflow velocity; ODXDYDZD is the inertial coordinate system; OXYZ is the body coordinate system; OWXWYWZW is the wind axis; OSXSYSZS is the rotor fixed axis coordinate; OSWXSWYSWZSW is the rotor wind coordinate, in which OSWXSW is the projection of the velocity V onto the OSWXSWYSW plane; βsw, the side-slip angle of the rotor is the angle between OSXS axis and OSWXSW axis.

Figure 2 .
Figure 2. Definition of coordinate systems.whereO is the center of gravity of the QTR; V is the inflow velocity; O D X D Y D Z D is the inertial coordinate system; OXYZ is the body coordinate system; O W X W Y W Z W is the wind axis; O S X S Y S Z S is the rotor fixed axis coordinate; O SW X SW Y SW Z SW is the rotor wind coordinate, in which O SW X SW is the projection of the velocity V onto the O SW X SW Y SW plane; β sw , the side-slip angle of the rotor is the angle between O S X S axis and O SW X SW axis.

Figure 3 .
Figure 3. Diagram of solving the rotor model.

Figure 3 .
Figure 3. Diagram of solving the rotor model.

Figure 4 .
Figure 4. Interactional diagram of front rotor wake to the rear rotor.(a) Top view of longitudinal interaction (helicopter mode).(b) Top view of longitudinal interaction (airplane mode).where V is inflow velocity.

Figure 4 .
Figure 4. Interactional diagram of front rotor wake to the rear rotor.(a) Top view of longitudinal interaction (helicopter mode).(b) Top view of longitudinal interaction (airplane mode).where V is inflow velocity.

Figure 5 .
Figure 5. Interactional diagram of front rotor wake to the rear wing.(a) Side view of interference.(b) Top view of interference.

Figure 5 .
Figure 5. Interactional diagram of front rotor wake to the rear wing.(a) Side view of interference.(b) Top view of interference.

Figure 7 .
Figure 7. Experiment of the rotor thrust and the rotor-wing interaction.(a) The wind tunnel.(b) Experiment of the rotor thrust.(c) Experiment of the rotor-wing interaction.

Figure 7 .
Figure 7. Experiment of the rotor thrust and the rotor-wing interaction.(a) The wind tunnel.(b) Experiment of the rotor thrust.(c) Experiment of the rotor-wing interaction.

Figure 7 .
Figure 7. Experiment of the rotor thrust and the rotor-wing interaction.(a) The wind tunnel.(b) Experiment of the rotor thrust.(c) Experiment of the rotor-wing interaction.

Figure 8 .
Figure 8. Comparisons between tests and calculations.(a) Required power of the isolated rotor.(b) Rotor thrust changes with velocity.(c) Percentages of wing's download to rotor thrust.

Figure 9 .
Figure 9. Conversion corridor of quad tiltrotor and selected trimming points.

Figure 9 .
Figure 9. Conversion corridor of quad tiltrotor and selected trimming points.

Figure 10 .
Figure 10.Diagram of the trim calculation procedure.

Figure 10 .
Figure 10.Diagram of the trim calculation procedure.

Figure 11 .
Figure 11.Comparison of trim results with and without total interference.(a) Comparison of collective stick displacement with and without interference.(b) Comparison of lateral stick displacement with and without interference.(c) Comparison of longitudinal stick displacement with and without interference.(d) Comparison of heading stick displacement with and without interference.(e) Comparison of rolling angle with and without interference.(f) Comparison of pitching angle with and without interference.

Figure 12 .Figure 13 .
Figure 12.Influence on body forces with different interaction.(a) Body forces increment caused by rotor-lateral interaction.(b) Body forces increment caused by rotor-longitudinal interaction.(c) Body forces increment caused by rotor-wing interaction.(d) Body forces increment caused by total interaction.where ΔX, ΔY, ΔZ represent the body forces increment along body axis.

Figure 12 .Figure 12 .Figure 13 .
Figure 12.Influence on body forces with different interaction.(a) Body forces increment caused by rotor-lateral interaction.(b) Body forces increment caused by rotor-longitudinal interaction.(c) Body forces increment caused by rotor-wing interaction.(d) Body forces increment caused by total interaction.where ∆X, ∆Y, ∆Z represent the body forces increment along body axis.

Figure 13 .
Figure 13.Influence on rotor thrust with different rotor-rotor interaction.(a) Rotor thrust increment caused by lateral interaction.(b) Rotor thrust increment caused by longitudinal interaction.

Figure 15 .
Figure 15.Vertical forces increment of front and rear wings caused by rotor-wing interaction.

Figure 14 .Figure 14 .
Figure 14.Influence on body moments with different interaction.(a) Body moments increment caused by rotor lateral interaction.(b) Body moments increment caused by rotor longitudinal interaction.(c) Body moments increment caused by rotor-wing interaction.(d) Body moments increment caused by total interaction.where ∆L, ∆M, ∆N represent the body moments increment along body axis.

Figure 15 .
Figure 15.Vertical forces increment of front and rear wings caused by rotor-wing interaction.

Figure 15 .
Figure 15.Vertical forces increment of front and rear wings caused by rotor-wing interaction.

Table 1 .
General parameters of the quad tiltrotor.

Table 1 .
General parameters of the quad tiltrotor.

Table 2 .
Categories of rotor-rotor interaction.

Table 3 .
Forms and numbers of control inputs.

Table 3 .
Forms and numbers of control inputs.

Table 4 .
Relationship between control value for each rotor and control input.