Research Article  Open Access
Distributed Modeling and Control of LargeScale Highly Flexible SolarPowered UAV
Abstract
The modeling, stability, and control characteristics of a large scale highly flexible solarpowered UAV with distributed allspan multielevons were presented. A geometrically nonlinear intrinsic beam model was introduced to establish the structural/flight dynamics coupled equation of motion (EOM); based on it, the explicit decoupled linear flight dynamics and structural dynamics EOM were derived through mean axis theory. Undeformed, deformed, and flexible models were compared through trimming and modal analysis. Since the deformation of wing has increased the UAV’s moment of inertia about the pitch axis, the frequency of short period mode has obviously decreased for the deformed model. The serious coupling between short period mode and 1st bending mode also significantly influences the roots of short period mode of flexible model. So flexible model was the only one which is able to accurately estimate the flight dynamics behaviors and was selected as the later control model. Forty distributed elevons and LQG/LTR controller were employed to control the attitude and suppress the aeroelastic deformation of the UAV simultaneously. The dynamics performance, robustness, and simulation results show that they were suitable for large scale highly flexible solarpowered UAV.
1. Introduction
The wing of high aspectratio solarpowered UAV is highly flexible and the surface density of structure is much lower than other kinds of aircraft, which may lead to the huge scale of solarpowered UAV in order to carry sufficient payload. These characteristics will lead to significant geometrically nonlinear aeroelastic deformation during flight and strong coupling between aeroelasticity and flight dynamics. If these issues have not been taken into full account, large error will occur. One of the key recommendations from the flight mishap investigation of the most famous solarpowered UAV “Helios” was to “develop more advanced, multidisciplinary (structures, aeroelastic, aerodynamics, atmospheric, materials, propulsion, controls, etc) “timedomain” analysis methods appropriate to highly flexible, “morphing” vehicles” [1].
In the modeling of “morphing” vehicles, fruitful works have been presented by Hodges, Cesnik, and their coworkers. The analysis method of Hodges et al. underlying NATASHA (Nonlinear Aeroelastic Trim And Stability of HALE Aircraft), a computer program that is based on geometrically exact, fully intrinsic beam equations and a finitestate induced flow model [2–4]. The University of Michigan’s Nonlinear Aeroelastic Simulation Toolbox (UM/NAST) is developed by Cesnik and his coworkers [5, 6], which includes nonlinear strainbased beam model, unsteady aerodynamics with simplified stall models, and a sixdegreesoffreedom flight dynamic formulation.
In the study of flight dynamics of highly flexible UAV, [2] found that the behavior is distinctly different from that of a rigid aircraft while taking into account the flexibility effects. Because the lowfrequency modes of structure are completely coupled to the flight dynamics mode, the phugoid as well as the short period mode is affected by wing flexibility significantly. Reference [3] furthered the study, the paper presented a theoretical basis for the flight dynamic response estimation of a highly flexible flying wing, multiple engines, and multiple control surfaces are taken into account. The results show that the phugoid mode may become unstable during large dihedral angles caused by large mass of payload pod. The pair of complexconjugate short period roots merges to become two real roots. This is because the deformed “U”shape leads to an orderofmagnitude increase in the pitch moment of inertia and so there is corresponding decrease in the frequency. Thus, this highly flexible aircraft does not show a classical shortperiod mode in its deformed state. Similar studies by Raghavan and Patil [4] and Su and Cesnik [5] confirmed this result. All the references paid attention to the nonlinear model and the results analysis, but the modeling of linear dynamics was ignored.
In [7], an optimal Static Output Feedback controller was designed for flutter suppression and gust load alleviation. However, only an elevon on the wing tip was employed in the controller, and the performance was highly dependent on sensor selection and placement.
In [8, 9], dynamic inversion control laws were presented for the trajectory control for Very Flexible Aircraft (VFA). However, in the use of dynamic inversion, there is an assumption that the stability and control derivatives remain fairly constant during the flight. This is not the case with highly flexible UAV.
In [10], a reducedorder controller was employed to control a high speed civil transport aircraft, and the results implied that the reducedorder controller was robust in the presence of unstructured uncertainty arising due to the model reduction error. However, the reduced order controller is weakly robust when employed to control the fullorder model with parametric uncertainty.
LQG/LTR is another popular robust control method. The primary objective of LQG/LTR is to recover the robustness of LQG, an excellent fullstate feedback robust controller. The procedure of LQG/LTR method is as follows: firstly, design a target feedback loop, which satisfies desired frequency performance, such as disturbance rejection and uncertainty deviation of the unmodeled dynamics, and then design a model based compensator to recover the properties of the target feedback loop [11]. LQG/LTR has been used widely in flight control, but only a few references are concentrated on the highly flexible aircraft.
In [12], an adaptive LQG/LTR controller was employed to control a very flexible aircraft, but the dynamics model was greatly simplified to capture only the unstable phugoid mode, and there were only 3 traditional control effectors, throttle, elevator, and aileron, available for the controller.
The large scale highly flexible UAV discussed in this paper is shown in Figure 1. It is equipped with 8 propulsion system nacelles (including secondary batteries, electrical motor, and propeller) and 5 vertical tails. The span is 70 m and the mean aerodynamics chord is 2.44 m. In order to enhance the control ability, trailingedge elevon spans the entire wing. Similar to the Helios [1], the elevon is divided into 40 distributed parts to release the notable aeroelastic deformation caused by large scale highly flexible wing.
For the aeroelastic flight dynamics analysis of such a UAV, most papers ignored the aerodynamics force distribution characteristic along the span and assumed the center of gravity (CG), elastic center (EC), and aerodynamics center (AC) to be coincident for simplification. In this paper, a more practical UAV model is proposed, in which the AC and EC is separated, the 3D quasisteady aerodynamics force is considered, and aerodynamics derivatives in every wing sections are derived for accuracy. In order to get a more accurate characteristic of solarpowered UAV, the mass parameters derived from conceptual design are employed, shown in Table 1.

In this paper, a nonlinear coupled structural and flight dynamics model based on geometrically nonlinear intrinsic beam was established first, and then its explicit linear dynamics model is derived for decoupled dynamics modal analysis by mean axis theory. At last, a LQG/LTR controller is designed with distributed 40 elevons for attitude control and aeroelastic deformation suppression for the largescale highly flexible solarpowered UAV.
2. Modeling of Nonlinear Structural/Flight Dynamics
2.1. Structural Dynamics Model
The geometrically nonlinear structural dynamics and flight dynamics equations are based on the intrinsic beam equation derived by Hodges et al. shown as follows [3, 13]:
The equations above are the core of structural dynamics, but they still cannot be used for the analysis of aerodynamics and flight dynamics without the relationship between displacements and rotations. Hence, the kinematics motion must be established to obtain the relationship between strain and displacement. The position vector of the origin of the deformed beam frame (B) in the root frame (R) and their rotational relationship are [14]:
The spatial derivatives can be replaced by the following equations, namely, central difference method [15]:where is the distance between two points.
The mean axis is applied to simplify and linearize the equation of motion. In the mean axis (M) frame, the linear and angular momentums due to elastic deformation are zero at every instant, which eliminates the coupling between the rigid body modes and the structural mode and results in a reducedorder system with only six degreesoffreedom.
Since the equations in (1) are written in the Bframe, the equations of motion of rigid mode can be derived by premultiplying all terms by to obtain components in the M frame and integrating in it along the beam axis as follows [16]:where the linear and angular velocities are defined asThe relationship of coordinate transformation isThe total moment of inertia of the UAV is solved as follows:
If the relative position of the rigid body does not vary significantly, can be approximated to be zero. If the motion of the UAV is limited to longitudinal only, the last two terms at the right hand of (5) can also be equal to zero.
2.2. Force and Moment Model
2.2.1. Internal Forces and Moments Model
The wing of solarpowered UAV may undergo large deformations during normal flight, while the strain is small, exhibiting geometrically nonlinear behaviors. Therefore, the smallstrain assumption is suitable for this problem. The internal forces and moments are linearly related to the strains and curvatures, shown in (9) [17]:
The wing is assumed to be a prismatic and isotropic beam with the shear center and tension axis coincident with the reference axis, by which the above equation can be derived as follows:where represents the shear modulus; represents Young’s modulus; represents the torsional constant; and represent the area moment of inertia along the and direction.
2.2.2. Aerodynamics Model
In order to analyze the motion of the flexible UAV exactly, the aerodynamic force of vibratingwings should be analyzed. Unsteady aerodynamics models that are useful for aeroelastic analysis are usually Theodorsen unsteady aerodynamic model and finite state induced flow model of Peters and Cao [18]. The Theodorsen model is derived in frequency domain, but it can be converted into time domain by rational function approximation. The finite state induced flow model is derived in time domain, and it can be used in a statespace model easily. These two types of aerodynamics are equivalent and can be derived from one another through Laplace and Fourier transforms.
In this paper, the Theodorsen unsteady aerodynamics model is adopted, which has been used widely for airfoil or large aspect ratio wings, and its quasisteady form is relative to the traditional aerodynamics derivatives used in rigid aircraft.
For the Theodorsen model, the wing section is described in Figure 2, where the definitions of points are as follows: : leading edge, : center of gravity of the wing section, : center of gravity of the UAV, : aerodynamics center, : middle point of chord, : elastic center, the position in direction of frame is , : control point of unsteady aerodynamics, and : aerodynamics center of control surface.
The definitions of variables are as follows: : angle of attack, : length of 1/2 chord, : deformation defined in frame, and frame: the origin is at LE, along the chord and pointing to the trailing edge, perpendicular to and point up.
The formulation of Theodorsen unsteady aerodynamics force and its moment to the elastic center arewhere is the theoretic lift slope of airfoil, the steady aerodynamics center is at 1/4 chord and unsteady aerodynamics center and its control point is at 3/4 chord. All these variables can be replaced by the actual airfoil data. is the Theodorsen function, produced by circulation. is the reduced frequency. If is low enough and , the flow will become quasisteady. As the flight dynamics and control of solarpowered UAV on low frequency is more concerned, quasisteady flow assumption is adopted for analysis in this paper.
By applying the quasisteady flow assumption and ignoring the terms in high orders, (11) becomeswhere the airfoil data have been replaced by the actual ones.
The physical meaning of the above equations is that the lift of wing section is produced by steady and quasisteady flow together, and the steady part is produced by the angle of attack of the section and its increment caused by relative elastic deformation. The quasisteady part is caused by .
It should be noted that, in most cases, the angles of attack in every sections (denoted as ) are not easy to be measured directly. Hence, the relationship between and incident angle at the CG of UAV expressed in frame (notated as and ) can be calculated aswhere is the Euler angle from frame to frame and the subscripts , , and represent the 3 components of Euler angle, which correspond to the roll, pitch, and yaw angle of the aircraft, respectively.
If the wing rotates along the axis only, (13) can be simplified as
The magnitude of drag is small, and its influence on the flight dynamics characteristics is limited, so the drag on every wing section can be divided from the total drag. The expression of total drag is
The aerodynamics force expressed in frame can be derived by the following axis conversion:where
Substituting (17) into (16), the nondimensional force coefficients at the coordinate are derived as follows:
Changing the reference point of moment into CG, the aerodynamics moment to CG of UAV is
It should be noted that the method proposed in this section can be used to calculate either the aerodynamics force or the moment of vertical tails, although the incident angle and Euler angle may vary significantly.
2.2.3. 3D Effect Correction of Aerodynamics
In the analysis of the aeroelastic flight dynamics of solarpowered UAV, most papers ignored the aerodynamics force distribution characteristics along the wing span and considered the flow condition to be same as the airfoils. Actually, the wing span is finite, which will make the flow of each section affect each other; that is to say, 3D effect must be considered for precision. In this paper, 3D effect is corrected based on liftingline theory, and it applies well to wings that have aspect ratio greater than 5 and sweep angle less than 20 degrees.
According to liftingline theory, the lift of airfoil is determined by effective angle of attack , which is formulated as follows:where and represent geometric and induced angle of attack, respectively.
According to the relationship of trailing vortex and circulation, the equation above becomeswhere is an arbitrary spanwise position on the wing and refers to theoretical lift slope of airfoil, which can be substituted by true airfoil value in the actual computation. Since (21) is a nonlinear integraldifferential equation, coordinate transformation method should be applied before solving it for simplification [19].
After solving the circulation distribution from (21), the lift distribution can be obtained from the KuttaJoukowski theorem:
By repeating the above steps for different geometric angles of attack and taking a derivative, one can obtain the lift slope of each cross section of the wing considering 3D effect correction.
3. Model Linearization
3.1. Linearized Aerodynamics
The movement of wing section is separated into two parts: the rigid movement of UAV and the elastic movement of structure. Then, the lift coefficients applied on the wing section can be deployed in the frame as follows:
It should be noted that the variables in the above equations are measured in the frame and frame.
According to the 2nd equation of (12), the aerodynamics moment to the elastic center (EC) is
From the 2nd equation of (18), the difference of lift expressed in the frame and frame is only . Hence, by replacing in (23) with , the lift expressed in the frame will be derived.
Substituting the expression of lift at frame into (19), the pitch moment about the center of gravity of the UAV can be obtained as follows:
3.2. Linearized Flight Dynamics Equation
It can be found that the form of (4) and (5) is similar to that of flight dynamics equations of normal rigid vehicles. Hence, they can also be linearized by small disturbance method and written into two sets of equations: longitudinal and lateraldirectional. Different from the normal rigid vehicle, the external force such as aerodynamics, thrust, and gravity are related to rigid motion and elastic motion for elastic model. The linear equation of flight dynamics for elastic UAV is where is the independent variable of rigid motion, with in longitudinal equations; is the independent variable of elastic motion; is the control variable of rigid motion, ; that is, all the elevons can be controlled independently, but all throttles act conformably again for simplification; is the same as conventional rigid vehicles as shown in [20]; and is similar to the rigid model also, but the column is equal to the length of ; is the influence matrix of elastic motion to rigid motion, for longitudinal motion:
The details of can be found in [21].
3.3. Linearized Structural Dynamics Equation
Ignoring nonlinear coupling terms in (1), and assuming that the translational strain is small and the external moment and the moment of inertia of and direction are equal to zero, (1) can be reduced into
The equations above can be rearranged into a compact form as follows:
Since the focus of this paper is the flight dynamics characteristics of the total UAV, the structural dynamics characteristics are used to derive the increment of aerodynamics force caused by aeroelasticity only and many assumptions can be introduced to linearize the equations above for simplification.
Equation (29) can be simplified to the form as in (30) by rewriting in frame and assuming that the structural elastic vibrations are small perturbations near the static trim state and the elastic motion in direction is ignorable. In the case of geometry, it is assumed that the frame is parallel to frame and the sweep angle of wing is zero. Besides, in mean axis system, the elastic motion is influenced by rigid motion through aerodynamics force only:where represents the angle between axis of frame and the plane of frame. All variables in (30) are in incremental form and the notation has been neglected for simplification.
Replacing the spatial derivatives by (3), and knowing the fact that the wing tip is free, the boundary conditions can be defined that the bending moments and shear forces at the first and last nodes are zero for the 2nd equation in (30); that is [15],Thus, the structural dynamics equation can be derived as follows: where
For the 2nd equation in (30), it is the same as (33); that is,
For the 3rd equation in (30), the boundary conditions are [15]
Then, the structural dynamics equation becomeswhere
In order to further simplify the linear structural dynamics equations, (33), (35), and (37) can be expanded in the modal coordinate system, and the final form can be derived as follows:where the independent variables arewhere is the modal displacements of structural elastic motion, which is the symmetric modes for longitudinal motion; is the transfer matrix from modal coordinate to physical coordinate.
For longitudinal motion, the coefficient matrices in (39) arewhere is related to some derivatives (e.g., ). The exact expressions are the same as the derivatives of general flight dynamics equations. More details about them can be referred to [21].
3.4. Total Linearized Equation
Combining (26) and (39), the final total linear equations can be derived as follows:The independent variable isand the state matrix is And the control matrix is:
The eigenvalues of the above state matrix are the modal eigenvalues of the elastic UAV considering the coupling of aeroelastic and flight dynamics:
4. Model Verification
4.1. Geometrically Nonlinear Deformation
In order to validate the effectiveness of the geometrically nonlinear structural model and central difference method proposed in this paper, a cantilever beam with bending rigidity and length = 1 m is employed. While applying increasing moment at the tip of the beam, the deformation of the beam is shown in Figure 3. It can be seen that while , the angular displacement is , the tip coincides with the root and the beam is deformed into a circle, which is consistent with the analytical results [16]. The effectiveness of structural model is validated.
4.2. Aerodynamics Model
Aerodynamics model is the most complex and important external force model in the problem of flight dynamics and control, so it must be checked carefully before further studies. The lift coefficient distribution along the wing span obtained from two methods are compared in Figure 4; one is derived from the ThinStrip liftingline theory proposed in Section 3.1 in this paper, and the other is the VortexLattice method from AVL software. It can be seen that the precision of ThinStrip liftingline theory model can satisfy the application requirements of this paper.
4.3. Linear Structural Model
In order to investigate the error of different modal truncation, 4 linear structural models of the wing are compared with the geometrically nonlinear one: truncating at the 1st, 2nd, 4th, and 6th bending mode, respectively, and then applying a flapwise impulse force to the tip of the wing. The time domain responses of the wing tip of different models are shown in Figure 5.
It can be seen that, if the deformation is not large enough, for the flapwise deformation, if we employ the first 6 bending modes to constitute the linear structural model, the responses of the wing tip are coincident with the nonlinear one well. Similarly, it can be found from simulation that the first 4 twist modes and 2 chordwise modes are precise enough. So, the first 6 flapwise bending, 4 twist and 2 chordwise bending modes are adopted to constitute the linear structural model of the wing in the following of this paper.
5. Flight Dynamics of Highly Flexible UAV
5.1. Trimming
Two kinetics models are employed for trimming: the “Undeformed” model represents the common rigid flight dynamics model; the “Deformed” model corresponding to (1) represents the flexible flight dynamics model.
In level flight, assuming that all the 40 distributed elevons deflect conformably for trimming for simplification, the highly flexible UAV will deform under the aeroelastic effect. In fact, the deformed UAV in trimming is shown in Figure 6, are the flapwise position and pitch angle of the UAV are shown in Figures 7 and 8.
It can be seen that the maximum wing deformation is 4.58 m at the wing tip, and the pitch angle at wing tip is larger than the root for about 3°. Larger pitch angle brings larger lift, and larger lift at wing tip will bring larger deformation. So, all the distributed elevons deflecting conformably is not the optimal control case while considering aeroelasticity.
5.2. Modal Analysis
The “Underformed” and “Deformed” models are corresponding to conventional rigid vehicles, and the linear “Flexible” model is based on (42). The longitudinal modes of different calculation models are shown in Table 2.

The results showed that the difference of these models is huge in longitudinal mode. Since the deformation of wing has increased the UAV’s moment of inertia about the pitch axis, the frequency of short period mode has obviously decreased for the deformed model. The serious coupling between short period mode and 1st bending mode also significantly influences the roots of short period mode, which then results in the phugoid changing from stable to unstable. These tendencies are similar to the results presented by Hodges et al. Through these analyses, it is clear that the flexible model is the only one which is able to accurately estimate the longitudinal flight dynamics behaviors.
6. Distributed LQG/LTR Controller Design
6.1. Statement of LQG/LTR Control
For a standard state space system described in (46), the design procedure of LQG/LTR controller in the output port is as follows:(1)Design a LQR regulator , which ensures that the singular value of open loop system on a certain frequency satisfies the robustness requirement.(2)Design a target loop: select appropriate and , so as to the singular value of fullstate feedback loop transfer function on the considering frequency close to that of adequately.
Where, the controlled plant isand the derived LQG/LTR controller plant is
The methodology of LQR is searching a full state feedback controller to minimize the quadratic cost function:subject to the system dynamics described in (46).
The expression of can be calculated bywhere is the solution of the following algebraic Riccati equation:The solution of in (48) is similar to ; one just needs to replace the matrix in (49)–(51) to , respectively, where represents the noise intensity; different noise intensity can be adjusted by ; represents the sensor noise intensity. A welldesigned LQG/LTR controller should be that, when , the singular value of on the considering frequency close to that of adequately. But it should be noted that larger will derive larger gain of .
6.2. Designing LQG/LTR Pitch Angle Controller
Before designing LQG/LTR control law, one shall augment the plant to satisfy the specific performance requirement. For example, adding integrators will assure the zero steadystate errors to step inputs. So, in order to control the pitch angle precisely, the integrator of pitch angle and velocity should be included in the control plant. Then, the control plant in (46) becomes
The critical factor for a practical controller is that not all the state variables are easy to be measured in engineering. So only a part of variables can be used to implement the output feedback control. Because the LTR technology can recover the stable margin of LQR, which is known as an excellent state feedback robust controller, the robustness of LQG/LTR output feedback controller can be guaranteed.
Because the flapwise bending motion of the wing is easy to be measured, so together with , , , , , are selected to be the measured variables for the controller. Then, we can apply the LQG/LTR method described above in (52).
If we select , the singular values of pitch angle loop of denoted by LQR and denoted by LTR are compared in Figure 9.
It can be seen from Figure 9 that, after the loop transfer recovery, the two singular values coincide perfectly on low frequency band, which indicates that the response of the close loop system to common control such as step input is quick and the tracking error is small even with model deviation. On high frequency band, the singular value of LTR decreases more quickly than that of LQR, which will provide excellent robustness to the high frequency disturbance such as unmodeled high frequency structural dynamics and unsteady aerodynamics. So, the LTR satisfies the required performance well.
6.3. Flight Simulation
There are 40 distributed elevons in the large scale highly flexible solarpowered UAV. The details about control actuator are omitted temporally for brief while designing the LQG/LTR controller. But in the stage of flight simulation, two kinds of actuators are employed to verify the robustness of the controller. The model of the elevons actuator isThe model of the throttles isTheir bandwidths are 20 and 2 rad/s, respectively.
Applying the LQG/LTR controller proposed above to the linear structural and flight dynamics model described in (52), when given a desired pitch angle , the longitudinal response of solar aircraft is shown in Figures 10–14. As shown in Figure 10, the adjust time of pitch angle is about 3 seconds; because of the use of integrator, steadystate error is eliminated. Seen from Figure 11, the variation of angle of attack in the whole process does not exceed , so the additional load it brings is very small, indicating that the overall movement of the highlyflexible UAV is relatively stable. Seen from Figure 12, the largest deflection angle of elevons is less than 20 degrees, and it occurs in the initial phase only, and when approaching a steady state, the deflections are no more than 5 degrees, far less than the limitation. It must be noted that the deflection angle of the 40 elevons are all different, because in addition to providing pitch angle control, distributed elevons also need to suppress the elastic deformation movement of the highlyflexible wing in the LQG/LTR controller. Seen from Figure 13, the maximum deformation of the wing occurs in the tip and the maximum magnitude of wing deformation is only about 0.09 m, far less than the deformation of 4.58 m in trimming, reminding that all the elevons deflect conformably in trimming. Seen from Figure 14, the max rate of twist angle is only 0.26 rad/s, corresponding to the reduced frequency of 0.0144. In a word, with distributed elevons, the use of LQG/LTR control method allows faster attitude control of the large scale highly flexible UAV while good aeroelastic deformation suppression is achieved simultaneously.
In addition, reminding that the linear structural model coincide with the nonlinear one perfectly on the condition of small deformation (shown in Figure 5). Since the LQG/LTR controller has restrained the structural deformation and incident angle in a relative small range, it can be deduced that the simulation results in this section are reliable.
7. Conclusions
In this paper, the modeling, stability, and control characteristics of a highly flexible large scale solarpowered UAV with distributed allspan multi elevons were presented.
A geometrically nonlinear intrinsic beam model is introduced to establish the equation of motion of flexible solarpowered UAV. Based on it, the decoupled linear flight dynamics and structural dynamics equations of motions are derived in mean axis frame. Steady aerodynamics model is derived by theoretical airfoil data and corrected for 3D effect by liftingline theory, quasisteady aerodynamics model is derived by Theodorsen unsteady aerodynamics model, and the aerodynamics derivatives of every wing sections are derived and transferred to different frames for the multirole in linear equations of flight dynamics or structural dynamics.
The nonlinear coupled structural and flight dynamics model was employed to trim the large scale flexible UAV, under the assumption that all the 40 distributed elevons deflect conformably. The maximum aeroelastic defection of wingtip is 4.58 m in flapwise, about 13% to the halfspan. The pitch angle at wing tip is larger than the root for about 3°. Larger pitch angle brings larger lift and leads to larger deformation. So all the distributed elevons deflecting conformably are not the optimal control case while considering aeroelastic effects.
Three different linear dynamics models, undeformed, deformed, and elastic, were compared for longitudinal mode analysis. The results showed that the difference of three models is huge. Since the deformation of wing has increased the UAV’s moment of inertia about the pitch axis, the frequency of short period mode has obviously decreased for the deformed model. The serious coupling between short period mode and 1st bending mode also significantly influences the roots of short period mode, which then results in the phugoid changing from stable to unstable. So, it is clear that the flexible model is the only one which is able to accurately estimate the longitudinal flight dynamics behaviors and control the large scale highly flexible UAV exactly.
Forty distributed allspan elevons were employed to control the attitude and suppress the aeroelastic deformation. LQG/LTR methodology was employed to accomplish the dynamics performance and robustness requirement. The flight simulation showed that the designed LQG/LTR controller can provide the adjust time of pitch angle in 3 seconds without steadystate error, and the maximum magnitude of wing deformation is only about 0.09 m, far less than the deformation of 4.58 m in trimming. So, LQG/LTR control method with distributed elevons is suitable for large scale highly flexible solarpowered UAV.
Nomenclature
Linear momentum  
:  Angular momentum 
:  Mass or moment 
:  Moment of inertia 
:  Linear velocity 
:  Angular velocity 
:  Strain 
:  Curvature 
:  External force 
:  External moment 
:  Chord 
:  Wing span 
:  Direction cosine matrix from to 
:  Angle of attack 
:  Angle of side slip 
:  Distance from the origin at , , and directions. 
:  Ground axis frame 
:  Deformed beam axis frame at wing root 
:  Deformed beam axis frame 
:  Undeformed beam axis frame 
:  Mean axis frame 
:  Structure axis frame 
:  Aero axis frame at deformed beam 
:  Aero axis frame at wing root 
:  Aero axis frame at mean axis frame 
:  Vector expressed in frame 
:  part of vector expressed in frame 
:  Spatial derivative of variable 
:  Temporal derivative of variable 
:  Mean value over the beam element 
:  Nodal value. 
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This study is supported by the National Natural Science Foundation of China (Grant no. 11202162) and the Scientific and Technology Foundation of State Key Laboratory on UAV (Grant no. 9140C390103130C39142).
References
 T. E. Noll, S. D. Ishmael, B. Henwood, M. E. PerezDavis, G. C. Tiffany, and J. Madura, “Technical findings, lessons learned, and recommendations resulting from the helios prototype vehicle mishap,” Meeting Proceedings RTOMPAVT145, 2007. View at: Google Scholar
 M. J. Patil, Nonlinear aeroelastic analysis, flight dynamics, and control of a complete aircraft [Ph.D. thesis], Georgia Institute of Technology, Atlanta, Ga, USA, 1999.
 M. J. Patil and D. H. Hodges, “Flight dynamics of highly flexible flying wings,” Journal of Aircraft, vol. 43, no. 6, pp. 1790–1798, 2006. View at: Publisher Site  Google Scholar
 B. Raghavan and M. J. Patil, “Flight dynamics of high aspectratio flying wings: effect of large trim deformation,” AIAA 2007–6386, The American Institute of Aeronautics and Astronautics, 2007. View at: Google Scholar
 W. Su and C. E. S. Cesnik, “Dynamic response of highly flexible flying wings,” AIAA Journal, vol. 49, no. 2, pp. 324–339, 2011. View at: Publisher Site  Google Scholar
 C. M. Shearer and C. E. S. Cesnik, “Nonlinear flight dynamics of very flexible aircraft,” Journal of Aircraft, vol. 44, no. 5, pp. 1528–1545, 2007. View at: Publisher Site  Google Scholar
 M. J. Patil and D. H. Hodges, “Output feedback control of the nonlinear aeroelastic response of a slender wing,” Journal of Guidance, Control, and Dynamics, vol. 25, no. 2, pp. 302–308, 2002. View at: Publisher Site  Google Scholar
 C. M. Shearer and C. E. S. Cesnik, “Trajectory control for very flexible aircraft,” AIAA 20066316, 2006. View at: Google Scholar
 B. Raghavan and M. J. Patil, “Flight control for flexible, highaspectratio flying wings,” Journal of Guidance, Control, and Dynamics, vol. 33, no. 1, pp. 64–74, 2010. View at: Publisher Site  Google Scholar
 X. Li and R. K. Agarwal, “Application of eeducedordermodels to robust control of the dynamics of a flexible aircraft,” AIAA 20035504, 2003. View at: Google Scholar
 M. Athans, “A tutorial on LQG/LTR method,” in Proceedings of the American Control Conference, Seattle, Wash, USA, June 1986. View at: Google Scholar
 T. E. Gibson, A. M. Annaswamy, and E. Lavretsky, “Modeling for control of very flexible aircraft,” in Proceedings of the AIAA Guidance, Navigation, and Control Conference, AIAA 20116202, Portland, Ore, USA, August 2011. View at: Google Scholar
 D. H. Hodges, “A mixed variational formulation based on exact intrinsic equations for dynamics of moving beams,” International Journal of Solids and Structures, vol. 26, no. 11, pp. 1253–1273, 1990. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 R. Palacios, J. Murua, and R. Cook, “Structural and aerodynamic models in nonlinear flight dynamics of very flexible aircraft,” AIAA Journal, vol. 48, no. 11, pp. 2648–2659, 2010. View at: Publisher Site  Google Scholar
 S. S. Rao, Mechanical Vibrations, Pearson Education, Upper Saddle River, NJ, USA, 4th edition, 2004.
 B. Raghavan, Flight dynamics and control of highly flexible flyingwing [Ph.D. thesis], Virginia Polytechnic Institute and State University, 2009.
 D. H. Hodges, “Geometrically exact, intrinsic theory for dynamics of curved and twisted anisotropic beams,” AIAA Journal, vol. 41, no. 6, pp. 1131–1137, 2003. View at: Publisher Site  Google Scholar
 D. A. Peters and W.M. Cao, “Finite state induced flow models part I: twodimensional thin airfoil,” Journal of Aircraft, vol. 32, no. 2, pp. 313–322, 1995. View at: Publisher Site  Google Scholar
 J. D. Anderson, Fundamentals of Aerodynamics, McGrawHill, 4th edition, 2007.
 Z. Fang, W. Chen, and S. Zhang, Aeronautic Vehicle Flight Dynamics, Beihang University Press, 1st edition, 2005 (Chinese).
 R. Wang, “Research on dynamics and control of highly flexible solarpowered UAV,” Postdoctoral Report, Northwestern Polytechnology University, Xi'an, China, 2011 (Chinese). View at: Google Scholar
Copyright
Copyright © 2015 Rui Wang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.