# Patent

## Title
Controlled transitory or sustained gliding flight with dihedral angle and trailing flaps

## Abstract

A micro aerial vehicle capable of controlled transitory or sustained gliding flight. The vehicle includes a fuselage. A pair of articulated wings are forward of a center of gravity of the vehicle, the wings being articulated and having trailing edge flaps, and having actuators for controlling the dihedral angles of the wings and the flaps for effective yaw control across the flight envelope. The dihedral angles can be varied symmetrically on both wings to control the aircraft speed independently of the angle of attack and flight-path angle, while an asymmetric dihedral setting can be used to control yaw and the actuators control the dihedral settings of each wing independently. The aircraft lacks a vertical tail or other vertical stabilizer.

## Background
PRIORITY CLAIM AND REFERENCE TO RELATED APPLICATION
The application claims priority under 35 U.S.C. §119 from prior provisional application Ser. No. 61/726,663, which was filed Nov. 15, 2012, and which is incorporated by reference herein.


STATEMENT OF GOVERNMENT INTEREST
This invention was made with government support under contract number FA95500910089 awarded by the Air Force Office of Scientific Research. The government has certain rights in the invention.


FIELD
A field of the invention is flight and flying devices, and particularly articulated wing micro aerial vehicles. Example applications of the invention include bio-inspired flight devices.
BACKGROUND
The dihedral angle is the upward angle measured from horizontal of the wings or tailplane of an aircraft. If there are multiple sections of a wing having different angles from horizontal, then the wing can be referred to as polyhedral. Wing dihedral angle is known to play a role in lateral-directional stability. In standard aircraft including almost all large commercial and private aircraft, the wings are typically set at a fixed positive dihedral angle for stability. Aircraft such as the Lockheed Martin F-104 and AV-8B Harrier use a negative dihedral (anhedral) setting for improved agility. These basic aspects have been well-studied, and are covered in many textbooks related to flight. A flexible wing deforms in a way that the effective dihedral of the wing increases. This only adds to the stabilizing effect of the dihedral. The changes resulting from flexibility have also been documented.
Unlike fixed wing aircraft, flapping wing aircraft possess the ability to change the dihedral angle on demand. An asymmetric wing setting naturally produces a yawing moment, however. Others have investigated the use of articulated wings for roll and yaw control. The approach studied has focused on controllable winglets. An example of this type of control is found in various papers published by Friswell, Bourdin et al. See, e.g., Bourdin et al., “Aircraft Control via Variable Cant-Angle Winglets,” Journal of Aircraft, Vol. 45, No. 2, March-April 2008. In this approach, the canted winglet in the inner wing is deflected upward to induce a turn. This reduces the net lift on the inner wing to induce the desired rolling moment, while the upward deflected winglet acts like a vertical tail and produces a proverse yawing moment in response to the rolling moment. This approach requires a segmented wing, and that only the winglet be deflected. With a monolithic wing, this approach can impair turning because the side force from the deflected wing may adversely dominate the side force required for turn (in a right turn, the deflected winglet produces a leftward force, whereas a rightward force is required to sustain/aid the turn).
Prior work by the present inventors and colleagues has advanced the state of micro aerial craft with articulated wings. Paranjape, A. A., Chung, S.-J., and Selig, M. S., “Flight Mechanics of a Tailless Articulated Wing Aircraft,” Bioinspiration & Biomimetics, Vol. 6, No. 2, (12 Apr. 2011) discloses a micro aerial aircraft with articulated wings and a horizontal tail (no vertical “tail”, and thus the reference to “tailless” aircraft in the title). This paper discloses that the wing dihedral can be varied symmetrically, along with the horizontal tail, to control the flight path angle independently of the fight speed. The micro aerial aircraft in the paper used dihedral angles that could be set asymmetrically for executing rapid zero sideslip turns. This work showed that the sign of the yaw moment derivative with respect to anti-symmetric dihedral depends strongly on the angle of attack and the angular rates. The primary demonstration of this work was to show that the yaw control effectiveness of the anti-symmetric dihedral depends primarily on the angle of attack, and also on the angular rates. When the angular rates are zero, the sign of the effectiveness depends on the sign (xaCL/c+Cm,ac), where xa/c is the non-dimensional distance between the center of gravity and the quarter-chord line. Furthermore, Cm,ac<0 for wings with a positive camber, and therefore, at small angles of attack, the control effectiveness is negative and it is positive at higher angles of attack. For an intermediate range of angles of attack, the sign depends on the angular rates as well. This can cause immense problems for yaw control, particularly when the angle of attack varies across the three regions in the course of a maneuver. There was not a solution proposed for these problems in this paper. This unresolved problem had to do with yaw control effectiveness changing its sign inside the routinely flown flight envelope. The sign was seen to be sensitive to the angle of attack as well as the angular rates. Basically, every control law assumes that the sign of control effectiveness is known; if it changes sign inside the envelope and does so as a function of several variables, it is virtually impossible to design a robust, reliable controller. Thus, this paper did not provide a practical solution for a tailless micro aerial aircraft.
SUMMARY OF THE INVENTION
An embodiment of the invention is a micro aerial vehicle capable of controlled transitory or sustained gliding flight. The vehicle includes a fuselage. A pair of articulated wings are forward of the center of gravity of the vehicle, the wings being articulated and having trailing edge flaps, and having actuators for controlling the dihedral angles of the wings and the flaps for effective yaw control across the flight envelope. The dihedral angles can be varied symmetrically on both wings to control the aircraft speed independently of the angle of attack and flight-path angle, while an asymmetric dihedral setting can be used to control yaw and the actuators control the dihedral settings of each wing independently. The aircraft lacks a vertical tail or other vertical stabilizers.



BRIEF DESCRIPTION OF THE DRAWINGS
FIGS. 1A and 1B show an articulated wing micro aerial vehicle (aircraft) of the invention with respective symmetric and asymmetric dihedral angles;
FIG. 1C shows another articulated wing micro aerial vehicle (aircraft) of the invention with a propulsion mechanism;
FIG. 2 is a block diagram of a closed loop control that is used to control actuated wing portions, elevator and associated flaps for an aircraft of the invention;
FIG. 3 illustrates the forces acting on a wing section;
FIG. 4 illustrates the physics underlying the use of wing dihedral as a control;
FIG. 5 illustrates the flight path angle as a function of the symmetric wing dihedral deflection;
FIG. 6 illustrates the yaw control effectiveness measure sign{Nδasym} is a function of the roll rate p and the yaw rate r;
FIG. 7 illustrates the effect of trailing edge flap deflection over a range of values;
FIG. 8 schematically illustrates the guiding of an aircraft for a change in flight path; and
FIG. 9 is a block diagram for a preferred embodiment inner loop control for the closed loop control of FIG. 2.



DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS
An embodiment of the invention is a vertical tailless micro aerial vehicle with a pair of articulated wings. The aircraft is controlled using wing articulation, together with trailing edge flaps, for effective yaw control across the flight envelope. The dihedral angles can be varied symmetrically on both wings to control the aircraft speed independently of the angle of attack and flight-path angle, while an asymmetric dihedral setting can be used to control yaw in the absence of a vertical tail. Actuators control the dihedral settings of each wing independently, and according to the mathematical formulation of closed loop control in the attachments, to achieve controlled sustained or transitory gliding flight. The wings are set ahead of the vehicle's center of gravity. Wing dihedral alone can be varied asymmetrically to perform rapid turns and regulate sideslip. The wings include trailing edge flaps. The trailing edge flaps are deflected downwards by a predetermined minimum amount.
Preferred embodiments of the invention will now be discussed with respect to the drawings. The drawings may include schematic representations, which will be understood by artisans in view of the general knowledge in the art and the description that follows. Features may be exaggerated in the drawings for emphasis, and features may not be to scale.
FIGS. 1A and 1B show an articulated wing micro aerial vehicle 10 of the invention with respective symmetric and asymmetric dihedral angles. The vehicle 10 is consistent with an experimental device that has been tested for gliding flight. The vehicle 10 is a bird-sized flapping wing aircraft which can provide sustained or transient gliding flight. An airframe 12 supports articulated wings 14 that each have a fixed section 14a closer to the airframe 12 and an actuated (elevator) portion 14b that can be actuated to change its dihedral angle. Each actuated portion 14b includes a trailing flap 16. The actuated portions 14b of the two wings 14 and their flaps 16 are actuated independently from each other. The actuated portions 14b are capable, therefore, of both asymmetric (FIG. 1A) and symmetric dihedral angles (FIG. 1B) accompanied by independent deflection angles for the flaps 16. An on-board controller 18 is mounted to the airframe and conducts closed loop control to maintain flight characteristics by controlling the dihedral angles and the flap deflection angles with actuators. The flaps are set statically because a static deflection suffices for uniform yaw control effectiveness. Dynamic flap deflection can induce undesirable motion in the plane of symmetry of the aircraft and undesirable rolling motion. The controller 18 includes flight sensors to provide data needed for closed loop flight control. A horizontal tail 20 provides stability, and can be a simple fixed horizontal tail. In preferred embodiments, there is no requirement that the tail 20 be actuated.
Other than the flapping wings 14 and trailing flaps 16, no additional yaw control mechanisms (such as a vertical tail, two degree of freedom horizontal tail, or split flaps) are required. The wings 14 shown can be considered rigid, but significantly flexible wings whose bending naturally creates a dihedral effect can be used with the invention while maintaining the control provided by the application of asymmetric and symmetric dihedral angles accompanied by independent deflection angles for the flaps 16. Methods and flight craft provide side force for yaw control in the same direction as that required for turning, which is advantageous compared to the winglet approaches that can create counterproductive side forces. FIG. 1C shows another exemplary micro aerial device 10a that includes a propulsion mechanism. The device 10a includes features of the devices of FIGS. 1A and 1B and these common features are labeled with like reference numerals. In addition, the device 10a includes a propeller 22 and motor 24. The device accomplishes sustained flight.
The following table defines some variable that are used in the following description.
















CL, CR
coefficients of lift and drag



Cm,ac
coefficient of pitching moment abot the




aerodynamic center yawing moment



N
yawing moment



α, β
angle of attack, sideslip



ψ, Θ, φ
Euler angles



γ
flight path angle



p, q, r
body axis roll, pitch and yaw rates



δL, δR
dihedral angle of left and right wing



δa, δs
asymmetric diherdral angle δR =




−δL = δa and symmetric dihedral




angle



δe
Elevator or horizontal tail deflection










FIG. 2 illustrates a closed loop control that is used to control actuated wing portions 14b and associated flaps. The control law has a two-tier hierarchical structure based on time-scale separation, which occurs naturally between the fast rotational dynamics and the slow translational dynamics. An inner loop 30 commands the elevator and the asymmetric components of the wing dihedral. An outer loop 32 commands the angle of attack and turn rate to be tracked by the inner loop based on flight speed and turn rate. The input to the outer loop 32 are position commands and position errors from feedback, and the output includes the alpha and gamma variables, and yaw rate (r). These are passed to the input of the inner loop 30, that then determines and outputs the control surface deflections. The turn rate and the light path angle are computed based on position measurements provided by the aircraft controller 18. Generally, by statically deflecting the trailing edge flaps on the movable sections of the wing, the control loop can ensure that the sign of the control effectiveness remains positive across the entire flight envelope. Flap deflection can be decided a-priori and held constant throughout the flight. Once the sign of the control effectiveness is known and is guaranteed to remain unchanged across the flight envelope, any of a wide range of control designs can be used to achieve guarantees on performance and stability for practical flight conditions. With the control of FIG. 2, the vehicles 10 and 10a are capable of agile flight in the flapping as well as gliding phases.
Control in accordance with FIG. 2 was used to control an experimental aircraft as shown in FIGS. 1A-1B. Control was accomplished by a combination of feedforward and feedback controllers. The angle of attack was controlled by the elevator/tail 20. The tail 20 preferably includes a trailing edge flap to operate as an elevator, but the entire tail can also deflect. In addition, static tails are permitted. The flaps 16 were deflected to a constant angle of 10 deg. The flight path angle was controlled by changing the symmetric component of the wing dihedral (denoted by δ), while the aircraft heading was controlled by anti-symmetric deflection of the two wings (denoted by δasym).

sin γ=cos α cos β sin θ−sin β sin φ cos θ−sin α cos β cos φ cos θ  (1)

sin χ cos γ=cos α cos β cos θ sin ψ+sin β(sin φ sin θ sin ψ+cos φ cos ψ)+sin α cos β(cos φ sin θ sin ψ sin φ cos ψ)  (2)

ω={dot over (χ)}=sign({dot over (χ)})√{square root over (p2+q2+r2)}  (3)

In the above equations, the flight path angle is γ, the wind axis heading angle is χ, and the turn rate is ω.
With a large horizontal tail, and the CG located approximately c/3 behind the wing aerodynamic center, the aircraft was sufficiently stable in pitch. Given the excellent open loop stability characteristics, the angle of attack was not controlled by feedback laws. Rather, the elevator deflection was set as a function of the commanded angle of attack. Specifically:




  
    
      
        
          δ
          
            e
            ,
            c
          
        
        =
        
          
            5
            3
          
          ⁢
          
            (
            
              15
              -
              
                α
                c
              
            
            )
          
        
      
    
    
      
        (
        4
        )
      
    
  




αc and δe,c are given in degrees instead of radians. Yaw control has been often neglected in past studies on perching aircraft, primarily because the aircraft considered in the literature had traditional roll and yaw control surfaces. Lateral directional control is an important concern for aircraft which lack a roll control surface and use a highly unconventional yaw control mechanism.

In an experimental controller, the wind axis heading angle χ was used as feedback for a PI yaw controller. While adding yaw rate feedback would provide for more robust control, the experiments demonstrated that the wind axis heading angle χ was sufficient for control. During test flight of experimental aircraft in as shown in FIGS. 1A and 1B, the dihedral angles of the two wings are controlled independently of each other. A left-right asymmetry in the wing dihedral angles and/or the force distribution on the wings leads to a net side force which, by the virtue of a non-zero distance between the aircraft center of gravity and the line of action of the force, translates into a yawing moment. Due to wing camber, the sign of the yawing moment for a given dihedral configuration depends on the angle of attack and the angular rates of the aircraft. Deflecting the trailing edge flaps on the articulated segments of the wing removes this non-uniformity.
Additional feedback control was also developed for the aerial vehicle of FIGS. 1A-1C. The control theory to be discussed includes use of the wing dihedral angles for control, use of the trailing edge flaps to mitigate control effectiveness problems, stability control, and closed loop control for perching. The relevant forces are discussed first.
FIG. 3 shows the forces acting on a wing section. The lift and drag both act in the plane of the airfoil; the lift is perpendicular to the local wind velocity, while drag acts along the local wind velocity. The quarter-chord pitching moment, Mac, is produced as a result of the pressure distribution on the airfoil and is independent of the angle of attack α. Lift, drag and the quarter-chord pitching moment, all of them per unit wing span, are written as follows:




  
    
      
        
          L
          =
          
            
              1
              2
            
            ⁢
            ρ
            ⁢
            
              
            
            ⁢
            
              V
              
                ∞
                c
              
              2
            
            ⁢
            
              
                C
                L
              
              ⁡
              
                (
                ∝
                )
              
            
          
        
        ,
      
    
    
      
        (
        5
        )
      
    
  
  
    
      
        
          D
          =
          
            
              1
              2
            
            ⁢
            ρ
            ⁢
            
              
            
            ⁢
            
              V
              
                ∞
                c
              
              2
            
            ⁢
            
              
                C
                D
              
              ⁡
              
                (
                ∝
                )
              
            
          
        
        ,
      
    
    
      
        (
        6
        )
      
    
  
  
    
      
        
          
            M
            ac
          
          =
          
            
              1
              2
            
            ⁢
            ρ
            ⁢
            
              
            
            ⁢
            
              V
              
                ∞
                
                  c
                  2
                
              
              2
            
            ⁢
            
              C
              
                m
                ,
                ac
              
            
          
        
        ,
      
    
    
      
        (
        7
        )
      
    
  




where c denotes the chord length, and the non-dimensional numbers CL(∝), CD(∝), and Cm,ac are the coefficients of lift, drag, and quarter chord pitching moment, respectively.

FIG. 4 illustrates the physics underlying the use of wing dihedral as a control. The key point is that changing the wing dihedral re-orients the lift vector with respect to the aircraft z-axis. Increasing the wing dihedral reduces the force acting in the body z-direction, and generates a side force. The reduced z-force manifests in the form of reduction in the net lift acting on the aircraft, which is accompanied by an incommensurate reduction in the drag force. Thus, changing the wing dihedral angle alters the lift-to-drag (L/D) ratio of the aircraft and offers the option of controlling the flight path angle and the aircraft speed independently of each other.
FIG. 5 illustrates the flight path angle as a function of the symmetric wing dihedral deflection. The points shown in the figure are equilibria computed at the same flight speed of 3 m/s.
The side force can be used for providing the centripetal force for turning, and as a source of yawing moment. In particular, if the CG (center of gravity) is located behind the line of action of the side force, then a positive (rightward) side force produces a positive yawing moment and vice-versa. It follows that a positive rolling moment (wherein the lift on the left wing is higher than the right wing) is accompanied by a positive yawing moment if the wings have a positive dihedral deflection. In aircraft which lack a vertical tail, a positive rolling moment is accompanied by a negative yawing moment which leads to an increase in sideslip, with the consequence that the aircraft faces an increased drag and potential reduction in lateral-directional stability. This yawing moment which acts in a sense opposite to the rolling motion is called adverse yaw. Adverse yaw is one of the primary limiting factors for the lateral-directional performance of a vertical-tailless aircraft, and a positive wing dihedral naturally suppresses adverse yaw produced due to rolling
The present inventors and colleagues have previously shown that that yaw control effectiveness of the anti-symmetric dihedral δasym depends not only on the angle of attack, but also on the angular rates. See, A. A. Paranjape, S.-J. Chung, and M. S. Selig, “Flight mechanics of a tailless articulated wing aircraft,” Bioinspiration & Biomimetics, vol. 6, no. 2, 2011, paper number 026005. The non-dimensional control effectiveness is given by

Nδasym=(xaCL+Cm,ac)  (8)

where CL and Cm,ac the coefficients of lift and quarter chord pitching moment, introduced in (1). The term xa denotes the non-dimensional (with respect to chord length) distance between the center of gravity and the quarter-chord line.

Therefore,

sign Nδasym=sign(xaCL+Cm,ac)  (9)

For positively cambered wings, Cm,ac<0 and therefore, at small angles of attack, where CL is small, the control effectiveness is negative. At higher angles of attack, the control effectiveness is positive.

FIG. 6 shows that, for an intermediate range of angles of attack, the sign depends on the angular rates as well. In FIG. 6, the yaw control effectiveness measure sign {Nδasym} is a function of the roll rate p and the yaw rate r. This can cause problems for yaw control, particularly when the angle of attack varies across the three regions in the course of a maneuver.
In the invention, the problem of non-uniform sign of control effectiveness is solved with trailing edge flaps. Trailing flap deflection leads to a greater increase in CL as compared to the reduction in Cm,ac. Using the airfoil theory, it can be shown that the change in CL and Cm,ac due to a flap deflection δf is given by




  
    
      
        
          Δ
          ⁢
          
            
          
          ⁢
          
            C
            L
          
        
        =
        
          (
          
            2
            ⁢
            
              (
              
                π
                -
                
                  θ
                  f
                
                +
                
                  2
                  ⁢
                  
                    
                  
                  ⁢
                  sin
                  ⁢
                  
                    
                  
                  ⁢
                  
                    θ
                    f
                  
                
              
              )
            
            ⁢
            
              δ
              f
            
          
        
      
    
    
      
        (
        10
        )
      
    
  
  
    
      
        
          
            Δ
            ⁢
            
              
            
            ⁢
            
              C
              
                m
                ,
                ac
              
            
          
          =
          
            
              -
              
                δ
                2
              
            
            ⁢
            sin
            ⁢
            
              
            
            ⁢
            
              θ
              f
            
            ⁢
            
              cos
              ⁡
              
                (
                
                  
                    θ
                    f
                  
                  -
                  1
                
                )
              
            
          
        
        ,
      
    
    
      
        (
        11
        )
      
    
  




Where θfε[0,π] depends on the location of the flap (xfc) from the leading edge:

cos θf=1−2xf  (12)

The term θf is defined purely for mathematical convenience in thin airfoil theory. For the example experimental aircraft consistent with FIGS. 1A and 1B, xf≈08.1 and xa=0.25. Thus, θf=2.2143, ΔCL=3.45 δf and ΔCm,ac=−0.14 δf.
FIG. 7 illustrates the effect of using a trailing edge flap from the range of 0-10 degrees of deflection. The flap deflection can be set as a function of α, and can be set to guarantee a certain positive control effectiveness. Consider an example needed effectiveness of at least 0.025 (corresponds to an α of 10 deg in FIG. 7). Then substituting the expression for ΔCL and ΔCm,ac, it follows that




  
    
      
        
          
            
              C
              L
            
            4
          
          +
          
            C
            
              m
              ,
              ac
            
          
          +
          
            0.72
            ⁢
            
              δ
              f
            
          
        
        =
        
          
            0.025
            ∴
            
              
                7.07
                +
                0.5
              
              ∝
              
                
                  -
                  0.1311
                
                +
                
                  0.72
                  ⁢
                  
                    δ
                    f
                  
                
              
            
          
          =
          
            
              0.025
              ⇒
              
                δ
                f
              
            
            =
            
              0.12
              -
              
                0.69
                ⁢
                
                  
                
                ⁢
                α
              
            
          
        
      
    
    
      
        (
        13
        )
      
    
  




Thus, flap deflection of nearly 7 deg is required at α=0, and no flap deflection is required beyond α=10 deg. A static flap deflection is preferred to avoid incurring any adverse effects from a dynamically varying flap deflection. In real terms, the flap deflection is preferably constant while the dynamics “settle down.” When the flap position is changed, the dynamics are perturbed and they settle down after some time has elapsed. The duration that it takes for the dynamics to settle down is called the settling time of the system. Therefore, instead of a statically fixed flap setting, a “piece-wise constant” flap deflection can also be used, where the flap deflection is held fixed for at least the settling time of the dynamics (with some factor of safety). However, the minimum flap deflection angle is still decided by FIG. 7. A larger flap deflection can be used for low speed flight as needed.
The benefit of uniformly positive control effectiveness, however, comes with a trade-off. The aircraft flies in a high-lift (it can be checked that CL>0.64), high draft configuration across the flight envelope. This slows the aircraft. However, the flight path angle can still be controlled effectively using symmetric dihedral deflection.
The rigid flight dynamics, together with the aerodynamics and kinematics, are highly nonlinear. The equations of motion, ignoring terms that arise from the angular velocity of the wing motion (due to flapping), have essentially the following structure:

m({dot over (u)}B+S(ωB)uB+(S({dot over (ω)}B)+S2(ωB))rcg)=Fnet 

J{dot over (ω)}B+S(ωB)JωB+m(S(rcg)uB+S(ωB)S(rcg)uB=Mnet  (14)

where m is the total mass of the aircraft mω is the mass of each wing, J is the moment of inertia tensor for the aircraft S(•) denotes the vector product, and Fnet and Mnet represent the net external (aerodynamic+gravitational) force and moment on the aircraft. Furthermore, ωB=[p,q,r]T is the vector representation of the aircraft angular velocity of the aircraft, with components in the aircraft body axes. The net aerodynamic force depends on the wing orientation. The position of the aircraft center of gravity is denoted by rcg, which is, in turn, approximated closely by




  
    
      
        
          r
          cg
        
        =
        
          
            
              
                
                  m
                  w
                
                ⁢
                b
              
              
                4
                ⁢
                
                  
                
                ⁢
                m
              
            
            ⁡
            
              [
              
                0
                ,
                
                  (
                  
                    
                      sin
                      ⁢
                      
                        
                      
                      ⁢
                      
                        δ
                        L
                      
                    
                    +
                    
                      sin
                      ⁢
                      
                        
                      
                      ⁢
                      
                        δ
                        R
                      
                    
                  
                  )
                
                ,
                
                  (
                  
                    
                      cos
                      ⁢
                      
                        
                      
                      ⁢
                      
                        δ
                        L
                      
                    
                    -
                    
                      cos
                      ⁢
                      
                        
                      
                      ⁢
                      
                        δ
                        R
                      
                    
                  
                  )
                
              
              ]
            
          
          T
        
      
    
    
      
        (
        15
        )
      
    
  




where αL and αR are the dihedral angles of the left and right wings, and b is the total wing span (so that each wing has length b/2). The force and moment vectors (and Fnet and Mnet) depend strongly on the dihedral angles of the wings. This is demonstrated by the yawing moment component of Mnett, which is given by

N=Zw,L(αL)sin δL−Zw,R(αR)sin δR  (16)

where αL and αR are local angles of attack of the left and right wings (which actually vary as a function of y, the spanwise coordinate and Zw,L and Zw,R are the local z-forces on the two wings. The terms αL, and αR themselves depend on δL and δR, respectively. For example, αR at a spanwise coordinate y is given by




  
    
      
        
          α
          R
        
        ≈
        
          
            β
            ⁢
            
              
            
            ⁢
            sin
            ⁢
            
              
            
            ⁢
            
              δ
              R
            
          
          +
          
            α
            ⁢
            
              
            
            ⁢
            cos
            ⁢
            
              
            
            ⁢
            
              δ
              R
            
          
          +
          
            
              py
              +
              
                
                  rx
                  a
                
                ⁢
                c
                ⁢
                
                  
                
                ⁢
                sin
                ⁢
                
                  
                
                ⁢
                
                  δ
                  R
                
              
            
            u
          
          +
          
            ry
            u
          
        
        ∝
        
          
            cos
            ⁢
            
              
            
            ⁢
            
              δ
              R
            
          
          +
          
            
              
                pry
                2
              
              ⁢
              cos
              ⁢
              
                
              
              ⁢
              
                δ
                R
              
            
            
              u
              2
            
          
        
      
    
    
      
        (
        17
        )
      
    
  




where xac is the distance between the center of gravity and the aerodynamic center of the wing. The sideslip is β, roll rate is p, and yaw rate is r. A similar expression can be written for αL.
In the simplest case, Zw,L and Zw,R are linear functions of αL and αR, respectively. Even then, the flight dynamics of aircraft with articulated wings are nonlinear and non-affine-in control. However, recognition of and accounting for the following flight dynamics is used in the invention to simplify control design:
1) The pitch dynamics can be controlled entirely by the elevator, and almost always independently of the lateral directional dynamics. This is true for most aircraft, except those that lack a horizontal tail;
2) For the micro aerial vehicles of the invention, the pitch dynamics (q, ∝) are stable. The only source instability is the lateral-directionals dynamics involving the yaw rate r and sideslip β;
3) The roll dynamics are stable, and they are not controlled directly, since the wing dihedral primarily controls yaw with very little roll control effectiveness.
With these factors in mind, the control laws are determined as follows. Consider the problem of controlling the yaw (β and r) dynamics. For the purpose of controlling the yaw dynamics, can be cast into the following control non-affine form

{dot over (η)}(t)=f(t,η(t),κ(t),u(t))  (18)

{dot over (κ)}(t)=ζ(t,η(t),κ(t),u(t),p(t))  (19)

where η(t)=β(t),r(t))ε2 represents the yaw dynamics, while uε is the control input, viz., the asymmetric dihedral deflection (δasym). The term p(t) represents other control inputs, namely the elevator deflection (δe) and symmetric dihedral deflection (δsym) which are used for longitudinal flight control. Finally κε6=[V∞,α,p,q,θ,φ] represents the rolling and pitching motion, as well as translation in the plane of symmetry. The flight dynamic modes corresponding to these six states are to be stable. One of the control objectives is to stabilize η(t)=β(t),r(t) in Equation (19), and ensure that it track a desired trajectory.
In embodiments of the invention, dynamic inversion laws are simplified to provide proportional-integral (PI) or proportional-integral derivative (PID) controllers. Specific example controls were based upon a class of dynamic inversion control laws defined in N. Hovakimyan, E. Lavretsky, and A. Sasane, “Dynamic inversion for nonaffine-in-control systems via time-scale separation. Part I,” Journal of Dynamical and Control Systems, vol. 13, no. 4, pp. 451-465 (2007). Our modification/simplification yields exact gain tuning laws, which allows the control gains to be linked explicitly to the convergence properties of the closed loop system as well as the tracking error bound.
Consider a general system described by equation (19), where η is no longer the yaw dynamics, but represents the state variables of interest for the purpose for control design. For now, we impose the additional condition ηεDη⊂, where Dη is compact. The case ηε2 is addressed later in the section. Let e(t)=η(t)−r(t) be the tracking error signal. Then the open loop error dynamics are given by

{dot over (e)}(t)=f(e(t)+r(t),κ(t),u(t))−{dot over (r)}(t),e(0)=e0 

{dot over (κ)}(t)=ζ(e(t)+r(t),κ(t),u(t)),κ(0)=κ0  (20)

where the unperturbed additional dynamics {dot over (κ)}(t)=ζ(0,κ(t),0) are assumed to be exponentially stabilized by to the control input p(t) (see (19)). We construct the dynamic inversion controller:




  
    
      
        
          
            ε
            ⁢
            
              
            
            ⁢
            
              
                u
                .
              
              ⁡
              
                (
                t
                )
              
            
          
          =
          
            
              -
              
                sign
                ⁡
                
                  (
                  
                    
                      ∂
                      f
                    
                    
                      ∂
                      u
                    
                  
                  )
                
              
            
            ⁢
            
              
                f
                ~
              
              ⁡
              
                (
                
                  t
                  ,
                  η
                  ,
                  κ
                  ,
                  u
                
                )
              
            
          
        
        ,
      
    
    
      
        (
        21
        )
      
    
  




with ε>0 sufficiently, and

{tilde over (f)}(t,η,κ,u)=f(t,e+r,κ,u)−{dot over (r)}(t)+ame(t)  (22)

and am>0 gives the desired rate of convergence of the closed loop dynamics.
Lemma 1. (Theorem 1, [16]) Given the system (10), the controller in (11) ensures that

    
    
        1) the tracking error e(t) ̃(ε), and
        2) the control u(t) converges to the isolated root of (22), i.e. u(t) makes f(η(t),κ(t),u(t))→(ame(t)+{dot over (r)}(t)).
    
    


The proof of this theorem is based on Tikhonov's theorem (H. K. Khalil, Nonlinear Systems, 3rd ed. Pearson Education, Upper Saddle, N.J., 2000, Theorem 11.1), and is also reproduced in N. Hovakimyan, E. Lavretsky, and A. Sasane, “Dynamic inversion for nonaffine-in-control systems via time-scale separation. Part I,” Journal of Dynamical and Control Systems, vol. 13, no. 4, pp. 451-465, 2007.
Remark 1. Although Lemma 1 guarantees that the tracking error  ̃(ε), it does not provide any guarantees on robustness. In particular, using a high gain can be detrimental in the presence of time delays.
Theorem 1. The control law in Equation (31) is equivalent to a PI controller with proportional (kp) and integral (ki) gains tuned to satisfy kp=1/ε and ki=amkp, where am is the desired time constant for the closed loop dynamics.
Proof: Since {dot over (n)}=f(t,η,κ,u) and ė={dot over (η)}−{dot over (r)}, we write the controller as




  
    
      
        
          
            
              
                
                  ε
                  ⁢
                  
                    
                  
                  ⁢
                  
                    
                      u
                      .
                    
                    ⁡
                    
                      (
                      t
                      )
                    
                  
                
                =
                
                ⁢
                
                  
                    -
                    sign
                  
                  ⁢
                  
                    (
                    
                      
                        ∂
                        f
                      
                      
                        δ
                        ⁢
                        
                          
                        
                        ⁢
                        u
                      
                    
                    )
                  
                  ⁢
                  
                    
                      f
                      ~
                    
                    ⁡
                    
                      (
                      
                        t
                        ,
                        η
                        ,
                        κ
                        ,
                        u
                      
                      )
                    
                  
                
              
              ,
            
          
        
        
          
            
              
                =
                
                ⁢
                
                  
                    -
                    sign
                  
                  ⁢
                  
                    (
                    
                      
                        ∂
                        f
                      
                      
                        δ
                        ⁢
                        
                          
                        
                        ⁢
                        u
                      
                    
                    )
                  
                  ⁢
                  
                    (
                    
                      
                        e
                        ⁡
                        
                          (
                          t
                          )
                        
                      
                      -
                      
                        
                          a
                          m
                        
                        ⁢
                        
                          e
                          ⁡
                          
                            (
                            t
                            )
                          
                        
                      
                    
                    )
                  
                
              
              ,
            
          
        
      
    
    
      
        (
        23
        )
      
    
  




Integrating both sides yields a PI controller of the form





  
    
      
        
          u
          ⁡
          
            (
            t
            )
          
        
        =
        
          
            u
            ⁡
            
              (
              0
              )
            
          
          -
          
            
              sign
              ⁡
              
                (
                
                  
                    ∂
                    f
                  
                  
                    ∂
                    u
                  
                
                )
              
            
            ⁢
            
              1
              ε
            
            ⁢
            
              (
              
                
                  e
                  ⁡
                  
                    (
                    t
                    )
                  
                
                -
                
                  e
                  ⁡
                  
                    (
                    0
                    )
                  
                
                -
                
                  
                    a
                    m
                  
                  ⁢
                  
                    
                      ∫
                      0
                      t
                    
                    ⁢
                    
                      
                        e
                        ⁡
                        
                          (
                          t
                          )
                        
                      
                      ⁢
                      
                        
                      
                      ⁢
                      
                        d
                        t
                      
                    
                  
                
              
            
          
        
      
    
    
      
        (
        24
        )
      
    
  




We may choose u(0)=−sign(∂f/∂u)εe(0). If kp and ki denote the proportional and integral gains of the PI controller, the they should be chosen to satisfy

ki=amkp, and kp=1/ε  (25)

so that




  
    u
    ⁡
    
      (
      t
      )
    
  
  =
  
    
      -
      
        sign
        ⁡
        
          (
          
            
              ∂
              f
            
            
              ∂
              u
            
          
          )
        
      
    
    ⁢
    
      (
      
        
          
            k
            p
          
          ⁢
          
            e
            ⁡
            
              (
              t
              )
            
          
        
        +
        
          
            k
            I
          
          ⁢
          
            
              ∫
              0
              t
            
            ⁢
            
              
                e
                ⁡
                
                  (
                  t
                  )
                
              
              ⁢
              
                
              
              ⁢
              
                d
                t
              
            
          
        
      
      )
    
  




This completes the proof of Theorem 1. ▪

Consider a second order system, {umlaut over (η)}(t)=f(t,η,{dot over (η)},κ,u). We can write it in the form:

{dot over (η)}1(t)=η2(t)

{dot over (n)}2(t)=f2(η1,η2,κ,t)+g2(t)u(t)  (26)

With η=η1 as the output, the equation for η2 is affine in u(t). The existence of f2(•) and g2(•) was shown in C. Cao and N. Hovakimyan, “L1 adaptive controller for a class of systems with unknown nonlinearities: Part 1,” in Proceedings of the American Control Conference, Seattle, Wash., 2008, pp. 4093-4098. Note that sign(g2)=sign (∂f/∂{dot over (η)}).
Define the desired value of n2(t) as

η2(t)=a(η1(t)−r(t))+{dot over (r)}(t)  (27)

where r(t) is the reference trajectory for η(t)=η1(t) and a>0. Define the error state for η2(t) as e2(t)=η2(t)−η2,d(t) whose dynamics are given by

ė2(t)=f2(η1,η2,κ,t)+g2(t)u(t)−η2,d(t)  (28)

From Theorem 1, the controller

u(t)=−sign(g2)(kpe2(t)+k1∫0te2(t)dt)  (29)

ensures that η2 tracks η2,d, where the gains kp and kI are chosen as per the guidelines in Theorem 1. The controller still needs further simplification to the PID form. Note that e2(t)=η2(t)+ae1(t)−{dot over (r)}(t)=ė1(t). Substituting into (29), we get

u(t)=−sign(g2)(kpė1(t)+(akp+k1)e1(t)+akI∫0te1(t)dt)  (30)

Which is a PID controller.

Theorem 2. The second-order system (26) can be stabilized using the PID controller (30), and it can be ensured that the tracking error between η(t) and the reference signal r(t) is bounded.
Proof: Theorem 1 guarantees that the control law (29) ensures that the tracking error e2(t) of the bounded (28) is bounded. Thus ∥η2−η2,d∥<O(ε) for some ε>0. Recalling the first equation {dot over (η)}(t)={dot over (η)}1(t)=η2(t), since η2(t)=e2(t)+η2,d(t), we can write

{dot over (η)}(t)=−a(η(t)−r(t)+{dot over (r)}(t)+e2(t)ė1(t){dot over (η)}(t)−{dot over (r)}(t)=ae1(t)+e2(t)  (31)

Since the unperturbed e1 dynamics (obtained by setting e2=0) are potentially stable, it follows from the Comparison Lemma (Lemma 9.1 in H. K. Khalil, Nonlinear Systems, 3rd ed. Pearson Education, Upper Saddle, N.J., 2000) that e1(t) is bounded if e2 is bounded. Moreover, the bound on e1 ̃(ε). This completes the proof. ▪
Remark 2. The following observations summarize the results for the closed loop controller that has been developed.
1) PI and PID controllers can be employed for nonlinear systems of the form (19) provided the additional (κ) dynamics are stable.
2) The dynamic inversion procedure yields a systematic gain tuning procedure (25) and (30).
The objectives of controlled perching can be stated as follows:
1) Apply a control law for the symmetric dihedral deflection (δsym=δL+δR/2) which ensures that the flight path zx tracks the desired profile zd(x). The desired profile zd(x) is a straight line connecting the initial point to the desired final point (xf,zf)
This is illustrated in FIG. 8, which shows the guidance problem. The aircraft is guided in the x-z plane along zd(x), while the y coordinate is regulated separately. Note that the z axis points upwards. Consequently, the x; y; z axes shown do not define a conventional right-handed frame, but are used as a reference purely as a matter of convenience.
2) Apply a control law for the anti-symmetric dihedral deflection (δsym=δL+δR/2) which ensures that y(x)→0 as x→x1 the desired final point. This requires that the yaw (r) dynamics be stabilized.
FIG. 9 shows a preferred controller for the inner loop control 30 and outer loop control of FIG. 2 controller. The control and path guidance algorithms differ from past efforts due to the unconventional control inputs δsym and δasym and in view of the use of trailing flaps. Additionally, the angle of attack is controlled by the elevator. In an example experiment, the flaps are deflected to a constant angle of 10 deg. The stability proofs are presented by rewriting the dynamic equations in the spatial domain rather than as functions of time, thereby permitting the use of dynamic inversion. The flight path angle (γ) is given by A. A. Paranjape, S.-J. Chung, and M. S. Selig, “Flight mechanics of a tailless articulated wing aircraft,” Bioinspiration & Biomimetics, vol. 6, no. 2, 2011, paper number 026005:

.sin γ=cos α cos β sin Θ−sin β sin φ cos Θ−sin α cos β cos φ cos Θ  (22)

while the wind axis heading angle χ is calculated as follows

sinx cosγ=cos α cos β cos φ sin ψ+sin β(sin φ sin Θ sin ψ+cos Θ cos ψ+sin α cos β(cos φ sin Θ sin ψ sin φ cos ψ  (32)

Yaw control has been often neglected in the literature on perching. On the other hand, lateral-directional control is an important concern for aircraft which lack a roll control surface, such as an aileron, and use a highly unconventional yaw control mechanism.
With a large horizontal tail, and the CG located approximately c/3 behind the wing aerodynamic center, the aircraft was sufficiently stable in pitch. Given the excellent open loop stability characteristics, the angle of attack was not controlled by feedback laws. Rather, on the basis of open loop glide tests, the elevator deflection was set as a function of the commanded angle of attack:




  
    
      
        
          
            δ
            
              e
              ,
              c
            
          
          =
          
            
              5
              3
            
            ⁢
            
              (
              
                15
                -
                
                  α
                  c
                
              
              )
            
          
        
        ,
        
          ...
          ⁢
          
            
          
          ⁢
          
            (
            
              α
              ⁢
              
                
              
              ⁢
              in
              ⁢
              
                
              
              ⁢
              degrees
            
            )
          
        
      
    
    
      
        (
        34
        )
      
    
  




where δe,c and α have been specified in degrees instead of the more conventional units of radians. Equation (34) is identical to equation (4), but is repeated here for clarity of explanation and ease of reference.
Consider the equation




  
    
      
        
          γ
          .
        
        =
        
          
            
              
                ρ
                ⁢
                
                  
                
                ⁢
                
                  SC
                  L
                
              
              
                2
                ⁢
                
                  
                
                ⁢
                m
              
            
            ⁢
            
              V
              ∞
            
            ⁢
            cos
            ⁢
            
              
            
            ⁢
            
              δ
              sym
            
          
          -
          
            
              g
              
                V
                ∞
              
            
            ⁢
            cos
            ⁢
            
              
            
            ⁢
            γ
          
        
      
    
    
      
        (
        35
        )
      
    
  




where m is the aircraft mass, S us the wing area and δsym is the symmetric dihedral deflection (δsym=δL+δR)/2). Note that




  
    
      
        
          
            d
            γ
          
          
            d
            x
          
        
        =
        
          
            
              γ
              .
            
            
              x
              .
            
          
          =
          
            
              γ
              .
            
            
              
                V
                ∞
              
              ⁢
              cos
              ⁢
              
                
              
              ⁢
              γ
            
          
        
      
    
    
      
        (
        36
        )
      
    
  




If we further assume that cos γ≈1, then the equations of translational motion simplify to the following form:




  
    
      
        
          
            
              d
              z
            
            
              d
              x
            
          
          =
          
            tan
            ⁡
            
              (
              
                γ
                ⁡
                
                  (
                  x
                  )
                
              
              )
            
          
        
        ⁢
        
          

        
        ⁢
        
          
            
              d
              γ
            
            
              d
              x
            
          
          =
          
            
              
                
                  ρ
                  ⁢
                  
                    
                  
                  ⁢
                  
                    SC
                    L
                  
                
                sm
              
              ⁢
              
                cos
                2
              
              ⁢
              
                
                  δ
                  sym
                
                ⁡
                
                  (
                  x
                  )
                
              
            
            -
            
              g
              
                V
                ∞
                2
              
            
          
        
      
    
    
      
        (
        37
        )
      
    
  




where δsym is the control input.

The system in (37) is in the form (26) except that derivatives and functions are defined with respect to x, not t. Hence, by replacing t and dt in (30) with x and dx, the following controller is designed (see Theorem 2):




  
    
      
        
          
            δ
            sym
          
          ⁡
          
            (
            x
            )
          
        
        =
        
          -
          
            (
            
              
                
                  k
                  p
                
                ⁢
                
                  
                    d
                    
                      e
                      z
                    
                  
                  
                    d
                    x
                  
                
              
              +
              
                
                  (
                  
                    
                      ak
                      p
                    
                    +
                    
                      k
                      I
                    
                  
                  )
                
                ⁢
                
                  
                    e
                    z
                  
                  ⁡
                  
                    (
                    x
                    )
                  
                
              
              +
              
                
                  ak
                  I
                
                ⁢
                
                  
                    ∫
                    0
                    x
                  
                  ⁢
                  
                    
                      
                        e
                        z
                      
                      ⁢
                      
                        
                      
                      (
                      x
                      )
                    
                    ⁢
                    
                      d
                      x
                    
                  
                
              
            
            )
          
        
      
    
    
      
        (
        38
        )
      
    
  




where ez(x)=z(x)−zd(x), and a is the desired rate of convergence of z(x) to the desired trajectory zd(x). The stability of the controller is guaranteed by Theorem 2. Note that the independent variable used in the controller is x, and not t. This is consistent with the recasting of the equations of motion in the spatial domain.
The motion in the x-y plane is given by

{dot over (x)}=V∞ cos γ cos X,{dot over (y)}=V∞ cos γ sin X  (39)

We assume that the yaw rate r settles rapidly to the commanded yaw rate rc, so that rc≈r={dot over (X)}/cos θ cos φ. The inner loop yaw controller in FIG. 9 ensures that r converges rapidly to rc. This is an expression of the time scale separation between the fast yaw (r) dynamics and the slower heading (y,X) dynamics. Thus we get




  
    
      
        
          
            
              d
              y
            
            
              d
              x
            
          
          =
          
            tan
            ⁢
            
              
            
            ⁢
            
              χ
              ⁡
              
                (
                x
                )
              
            
          
        
        ⁢
        
          

        
        ⁢
        
          
            d
            χ
          
          
            d
            x
          
        
        =
        
          
            
              V
              ∞
            
            ⁢
            γ
            ⁢
            
              
            
            ⁢
            
              
                χ
                .
              
              ⁡
              
                (
                t
                )
              
            
          
          =
          
            
              
                
                  V
                  ∞
                
                ⁢
                cos
                ⁢
                
                  
                
                ⁢
                γ
              
              
                cos
                ⁢
                
                  
                
                ⁢
                Θcos
                ⁢
                
                  
                
                ⁢
                φ
              
            
            ⁢
            
              
                r
                c
              
              ⁡
              
                (
                x
                )
              
            
          
        
      
    
    
      
        (
        40
        )
      
    
  




where φ is the bank angle of fuselage. Note that cos γ, and cos Θ, and cos φ are all positive since they lie in [−π, π] (very generally), which implies that the control coefficient V∞ cos γ/cos Θ cos φ is uniformly positive. The control problem is very similar to what is encountered for the flight path angle, show in (37), and we derive a similar controller as (38):

rc(x)=−(kp tan(X(x))+(akp+kI)y(x)+akI∫0xy(x)dx)  (41)

The stability of this controller is guaranteed by Theorem 2. Further, as in the case of flight path control, the use of dynamic inversion was considerably facilitated by rewriting the dynamical equations in the spatial domain which rendered them in a strict feedback form. Next, we describe the design of the inner yaw control loop in FIG. 9.
Remark 3. The choice of y=0 as the desired path can be replaced by any suitable path yd (x), such as a straight line connecting the initial point and the desired final point, as for the flight path angle guidance law.
The objective of the inner yaw control loop in FIG. 9 is to command the anti-symmetric wing dihedral, δasym, so that the rate r tracks the yaw rate commanded by the outer loop, rc. Yaw control is achieved using a PI controller, motivated by Theorem 1, given by

δasym=kp(rc(t)−r(t))+kI∫0trc(t)−r(t))dt  (42)

The proportional and integral gains were set to kp=kI=2 during experiments.
Remark 4.
In practice, the heading and flight path angles are computed and commanded preferably every 0:05 s. This prevents undesirable oscillatory behavior which arises due to continuous changes in the wing dihedral, and the underlying coupling between the wing dihedral angle and the pitching dynamics of the aircraft.
The experimental air craft consistent with FIGS. 1A and 1B have the following characteristics.









TABLE 2





Property
Metric Measurement
Units

















Mass (m)
44.0
g


Wing span (S)
41.8
cm


Wing Chord (c)
9.5
cm


Wing incidence angle (left and right)
6.0
deg


Wing dihedral
controlled-variable


MAV length
35.0
Cm


Elevator area
39.12
cm2









The articulated portion was hinged and included the outboard 60% of the wing. It has five control surfaces:
1) An elevator, which is a movable flap attached to the horizontal tail, and whose deflection is denoted by δe.
2) The dihedral angles (δR and δL) of the outboard segments of the right and left wing can be changed independently of each other. The actuators for changing the wing dihedral angle were attached on the lower surface of the center (non-rotating) wing section, along with the radio receiver.
3) The outboard segments were equipped with flaps which are capable of being actuated independently. In the experiments, both flaps were deflected by the same amount δf. Flaps are used for ensuring uniform yaw control effectiveness. These flaps, when deflected in an anti-symmetric manner, can also act as the traditional ailerons.
Both wings in the experimental aircraft can rotate from a maximum 45 deg dihedral to minimum −15 deg for a total arc range of 60 deg. The actuators for wing dihedral, are controlled independently on both wings. Digital actuators with a torque rating of 0.29 kg-cm are used to maneuver the wings. The time required for the wings to rotate from the minimum −15 deg to maximum 45 deg is about 0.05 s, and the actuators have a time-delay of 200 ms.
While not bound by this theory and not needed to demonstrate patentability, the experiments demonstrate that the invention provides what is believed to be the first perching demonstration on a laterally unstable aerial robot. In experiments, the lateral-directional motion of the aircraft was controlled actively for the entire duration of the maneuver. This complements and advances the state of the art experiments reported in the literature which were concerned almost exclusively with the longitudinal motion of stable aircraft.
The experimental aircraft used variable, asymmetric wing dihedral effectively to control the flight path as well as the heading of the aerial robot. Trailing edge flaps ensure that the wing dihedral provided uniform yaw control effectiveness across the flight envelope. The closed loop guidance laws of the invention provide perching through control developed by rewriting the equations of motion in the spatial domain and applying dynamic inversion-motivated PID control. Experiments demonstrate the ability for the aircraft to conduct controlled guided flight and perch on hand, which is a significant advance in developing micro aerial vehicles capable of close operational interaction with humans.
The above example aircraft has been flight-tested. During flight, the dihedral angles of the two wings are controlled independently of each other. A left-right asymmetry in the wing dihedral angles and/or the force distribution on the wings leads to a net side force which, by the virtue of a non-zero distance between the aircraft center of gravity and the line of action of the force, translates into a yawing moment. Due to wing camber, the sign of the yawing moment for a given dihedral configuration depends on the angle of attack and the angular rates of the aircraft. Deflecting the trailing edge flaps on the articulated segments of the wing removes this non-uniformity.
The deflection of trailing edge flaps symmetrically in methods and aerial vehicle of the invention, achieves an increase in the sectional lift coefficient to the extent that the yaw control effectiveness is uniformly positive across the flight envelope. It is well-known that trailing edge flaps increase the sectional lift coefficient, and calculations from thin airfoil theory show that this increase in lift is much more than the increase in the magnitude quarter chord pitching moment coefficient. However this use of trailing edge flaps runs counter to intuition for improving the yaw control authority, and would have been generally viewed: flap deflection increases the wing camber which is, in principle, responsible for a non-zero quarter-chord pitching moment in the first place.
While specific embodiments of the present invention have been shown and described, it should be understood that other modifications, substitutions and alternatives are apparent to one of ordinary skill in the art. Such modifications, substitutions and alternatives can be made without departing from the spirit and scope of the invention, which should be determined from the appended claims.
Various features of the invention are set forth in the appended claims.

## Claims

1. A micro aerial vehicle capable of controlled transitory or sustained gliding flight, comprising:
a fuselage; and
a pair of articulated wings forward of a center of gravity of the vehicle, the wings being articulated and having trailing edge flaps, and having actuators for controlling the dihedral angles of the wings and the flaps for effective yaw control across the flight envelope, wherein the dihedral angles can be varied symmetrically on both wings to control the aircraft speed independently of the angle of attack and flight-path angle, while an asymmetric dihedral setting can be used to control yaw and the actuators control the dihedral settings of each wing independently;
wherein the aircraft lacks a vertical tail or other vertical stabilizer.

2. The vehicle of claim 1, having closed loop control means for controlling the actuators for sustained or transitory gliding flight.

3. The vehicle of claim 2, wherein the closed loop control means further controls the actuators for perching.

4. The vehicle of claim 3, further comprising an elevator and associated actuator, and wherein said closed loop control means further control the associated actuator for sustained or transitory gliding flight and for perching.

5. The vehicle of claim 1, further comprising:
an elevator and associated actuator; and
a closed loop controller, wherein the controller controls the actuators and associated actuator such that dihedral angles of the two wings are controlled independently of each other; a left-right asymmetry in the wing dihedral angles is used to produce yawing moment, wherein the sign of the yawing moment for a given dihedral configuration depends on the angle of attack and the angular rates of the aircraft and wherein the trailing edge flaps are set to guarantee uniform yaw control effectiveness across the flight envelope.

6. The vehicle of claim 1, further comprising:
an elevator and associated actuator; and
a controller, the controller that applies a two-tier hierarchical structure based on time-scale separation, with an inner loop that commands the elevator and the asymmetric components of the wing dihedral and an outer loop commands the angle of attack and turn rate to be tracked by the inner loop based on flight speed and turn rate, and the controller deflects the trailing edge flaps to ensure that the sign of the control effectiveness remains positive across the entire flight envelope.

7. The vehicle of claim 1, further comprising:
an elevator and associated actuator; and
a controller that controls flight path angle by changing a symmetric component of the wing dihedral, controls heading by anti-symmetric component of the wing dihedral, and controls angle of attack by the elevator deflection.

8. The vehicle of claim 7, wherein the controller controls the elevator deflection as a function of a commanded angle of attach according to




  
    
      δ
      
        e
        ,
        c
      
    
    =
    
      
        5
        3
      
      ⁢
      
        (
        
          15
          -
          
            α
            c
          
        
        )
      
    
  
  ,




where αc and δe,c are given in degrees and δe,c is the elevator deflection and αc is the commanded angle of attack.

9. The vehicle of claim 7, wherein the controller uses a wind axis heading angle as feedback for a PI (proportional integral) yaw controller.

10. The vehicle of claim 1, wherein said pair of articulated wings are rigid.

11. The vehicle of claim 1, wherein said pair of articulated wings and the trailing edge flaps are actuated by the actuators independently from each other.