Next Article in Journal
Adaptive Neuro-Fuzzy Technique for Autonomous Ground Vehicle Navigation
Previous Article in Journal
A Review of Camera Viewpoint Automation in Robotic and Laparoscopic Surgery
Previous Article in Special Issue
Neural Networks Integrated Circuit for Biomimetics MEMS Microrobot
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Robust Nonlinear Regulation of Limit Cycle Oscillations in UAVs Using Synthetic Jet Actuators

by
Natalie Ramos Pedroza
1,*,
William MacKunis
1,* and
Vladimir V. Golubev
2
1
Department of Physical Sciences, Embry-Riddle Aeronautical University, 600 S Clyde Morris Blvd, Daytona Beach, FL 32114, USA
2
Department of Aerospace Engineering, Embry-Riddle Aeronautical University, 600 S Clyde Morris Blvd, Daytona Beach, FL 32114, USA
*
Authors to whom correspondence should be addressed.
Robotics 2014, 3(4), 330-348; https://doi.org/10.3390/robotics3040330
Submission received: 31 May 2014 / Revised: 22 July 2014 / Accepted: 15 August 2014 / Published: 26 September 2014
(This article belongs to the Special Issue Advances in Biomimetic Robotics)

Abstract

:
In this paper, a synthetic jet actuators (SJA)-based nonlinear robust controller is developed, which is capable of completely suppressing limit cycle oscillations (LCO) in UAV systems with parametric uncertainty in the SJA dynamics and unmodeled external disturbances. Specifically, the control law compensates for uncertainty in an input gain matrix, which results from the unknown airflow dynamics generated by the SJA. Challenges in the control design include compensation for input-multiplicative parametric uncertainty in the actuator dynamic model. The result was achieved via innovative algebraic manipulation in the error system development, along with a Lyapunov-based robust control law. A rigorous Lyapunov-based stability analysis is utilized to prove asymptotic LCO suppression, considering a detailed dynamic model of the pitching and plunging dynamics. Numerical simulation results are provided to demonstrate the robustness and practical performance of the proposed control law.

1. Introduction

There has recently been a surge of interest in the design and application of unmanned aerial vehicles (UAVs). These UAVs can be used in numerous civilian applications, such as urban reconnaissance, package delivery, and area mapping. Besides, UAVs are utilized in various military applications as well. One of the biggest challenges involved in the autonomous operation of UAVs is in the design of flight tracking controllers for UAVs operating in uncertain and possibly adverse conditions. In particular, suppression of limit cycle oscillations (LCO) (or flutter) is an important concern in UAV tracking control design. This is especially true for applications involving smaller, lightweight UAV systems, where the aircraft wings are more susceptible to LCO. These engineering challenges necessitate the utilization of UAV flight controllers, which achieve accurate flight tracking in the presence of dynamic uncertainty while simultaneously suppressing LCO. Moreover, as practical considerations motivate the implementation of smaller UAVs, there is a growing need for UAV flight control designs that do not require heavy mechanical deflection surfaces.
LCO refer to “flutter” behaviors in UAV wings that manifest themselves as constant-amplitude oscillations [1], which result from nonlinearities inherent in the aeroelastic dynamics of the UAV system [2]. Due to these behaviors, the LCO would surpass the limiting safe flight boundaries of an aircraft [3] and could potentially lead to structural damage and catastrophes. Control applications for LCO suppression are often developed (e.g., see [4,5,6]) using mechanical deflection surfaces (e.g. flaps, ailerons, rudders, and elevators). However, when dealing with small UAVs, practical engineering considerations and physical constraints can preclude the addition of the large, heavy moving parts that are required for the installation of deflection surfaces. To address this challenge, the use of synthetic jet actuators (SJA) in UAV flight control systems is becoming popular as a practical alternative to mechanical deflection surfaces.
The design and application of SJA has recently increased by virtue of their capability to achieve momentum transfer with zero-net mass-flux. This beneficial feature eliminates the need for an external fuel supply, since the working substance is simply the gas (i.e., air) that is already present in the environment of operation [7]. This makes SJA an attractive option in UAV applications, because of the significant reduction in the size of the required equipment. The SJA synthesize the jet flow through the alternating suction and ejection of fluid through an aperture, which is produced via pressure oscillations in a cavity [8] as shown in Figure 1. The pressure oscillations can be generated using various methods, including pistons in the SJA orifices [3] or piezoelectric diaphragms [7,9]. SJA can achieve boundary-layer flow control near the surface of a UAV wing [10], since they can provide instant actuation, unlike conventional mechanical control surfaces. In addition, SJA can expand the usable range of angle of attack, resulting in improved UAV maneuverability [11].
Recently developed nonlinear control methods using SJA typically utilize neural networks and/or complex fluid dynamics computations in the feedback loop (e.g., see [9,10,12,13,14,15,16,17,18,19,20,21]). While such approaches have been shown to yield good SJA-based control performance, they can require increased computational resources, which might not be available in small UAV applications. Adaptive control approaches have been applied to linear time-invariant (LTI) dynamic models to compensate for SJA nonlinearities and external disturbances [13]. Adaptive inverse control schemes are another popularly utilized method to compensate for the actuator nonlinearity inherent in SJA [9,14,15,16,17]. A robust tracking control method is proposed in [7], which builds on results similar to [9,14,15,16,17], to compensate for the SJA nonlinearity at a reduced computational cost. The aforementioned approaches have been shown to achieve good flight tracking performance using SJA; however, nonlinear control approaches have not as often been applied to SJA-based LCO suppression.
Figure 1. Schematic layout of a synthetic jet actuator.
Figure 1. Schematic layout of a synthetic jet actuator.
Robotics 03 00330 g001
In this paper, an SJA-based nonlinear robust controller is developed, which is capable of completely suppressing LCO in UAV systems with uncertain dynamics. This difficulty was handled via innovative algebraic manipulation in the error system development, along with a Lyapunov-based robust control law. A salient feature about this robust control law is that its structure is continuous, and so bounded disturbances can be asymptotically rejected without the need for infinite bandwidth. A rigorous Lyapunov-based stability analysis is utilized to prove asymptotic pitching and plunging regulation, considering a detailed dynamic model of the pitching and plunging dynamics. Numerical simulation results are provided to demonstrate the performance of the proposed control law.

2. Dynamic Model and Properties

The equation describing LCO in an airfoil approximated as a 2-dimensional thin plate can be expressed as
M s p ¨ + C s p ˙ + F ( p ) p = L i f t M o m e n t
where the coefficients M s , C s R 2 × 2 denote the structural mass and damping matrices, F ( p ) R 2 × 2 is a nonlinear stiffness matrix, and p t R 2 denotes the state vector. In Equation (1), p t is explicitly defined as
p = h α
where h t , α t R denote the plunging [meters] and pitching [radians] displacements describing the LCO effects. Also in Equation (1), the structural linear mass matrix M s [22]
M s = m S α S α I α
where the parameters S α , I α R are the static moment and moment of inertia, respectively. The structural linear damping matrix is described as
C s = 2 ζ h k h m 0 0 ζ a k α I α
where the parameters ζ h , ζ α R are the damping logarithmic decrements for plunging and pitching, and m R is the mass of the wing, or in this case, a flat plate. The nonlinear stiffness matrix utilized in this study is
F ( p ) = k h 0 0 k α + k α 3 α 2
where k α , k α 3 R denote structural resistances to pitching (linear and nonlinear) and k h R is the structural resistance to plunging.
In Equation (1), the total lift and moment are explicitly defined as
L i f t M o m e n t = ( L + L v j ) ( M + M v j ) = M a p ¨ + C a p ˙ + K a p + L η η + B u
where L v j t , M v j t R denote the equivalent control force and moment generated by the jth SJA, and L t , M t R are the aerodynamic lift and moment due to the 2 degree-of-freedom motion [10]. In Equation (6), η t R 2 denotes the aerodynamic state vector that relates the moment and lift to the structural modes. Also in Equation (6), u t R 2 denotes the SJA-based control input (e.g., the SJA air velocity or acceleration), and B R 2 × 2 is an uncertain constant input gain matrix that relates the control input u t to the equivalent force and moment generated by the SJA. Also in Equation (6), the aerodynamic and mode matrices M a , C a , K a , L η R 2 × 2 are described as
M a = π ρ b 2 1 b a b a b 2 1 8 + a 2
C a = π ρ b 2 0 U 0 U b 1 2 a + 2 π ρ U b ϕ ( 0 ) 1 b 1 2 a b 1 2 + a b 2 1 2 + a 1 2 a
K a = 2 π ρ U b ϕ ( 0 ) 0 U 0 b 1 2 + a U
L η = 2 π ρ U b a 1 b 1 a 2 b 2 b 1 2 + a a 1 b 1 b 1 2 + a a 2 b 2
where ϕ ( 0 ) is the Wagner solution function at 0, and the parameters a 1 , b 1 , a 2 , b 2 R are the Wagner coefficients. In addition, a , b R denote the relative locations of the rotational axis from the mid-chord and the semi-chord, respectively. The aerodynamic state variables are governed by [22]
η ˙ = C η p ˙ + K η p + S η η
The aerodynamic state matrices in Equation (11), C η , K η , S η R 2 × 2 , are explicitly defined as
C η = U b 1 b 1 2 a 1 b 1 2 a
K η = U b 0 U 0 U
S η = U b b 1 0 0 b 2 .
By substituting Equation (6) into Equation (1), the LCO dynamics can be expressed as
p ¨ = M 1 C p ˙ M 1 K p + M 1 L η η + M 1 B u
where C = C s C a , K = F ( p ) K a , and M = M s M a . By making the definitions x 1 t = h t , x 2 t = α t , x 3 t = h ˙ t , x 4 t = α ˙ t , x 5 t = η 1 t , and x 6 t = η 2 t , the dynamic equation in Equation (15) can be expressed in state form as
x ˙ = A ( x ) x + B ¯ u
where x t R 6 is the state vector, A ( x ) R 6 × 6 is the state matrix (state-dependent). In Equation (16), the input gain matrix B ¯ R 6 × 2 is defined as
B ¯ = 0 2 × 2 M 1 B 0 2 × 2
where 0 2 × 2 denotes a 2 × 2 matrix of zeros. The structure of the input gain matrix in Equation (17) results from the fact that the control input u t only directly affects h ¨ t and α ¨ t .

3. Control Development

The objective is to design the control signal u t to regulate the plunge and pitching dynamics (i.e., h t , α t ) to zero. To facilitate the control design, the expression in Equation (15) is rewritten as
M p ¨ = g ( h , α , η ) + B u
where g ( h , α , η ) is an unknown, unmeasurable auxiliary function.
Remark 1. Based on the open-loop error dynamics in Equation (18), one of the control design challenges is that the control input u t is premultiplied by the uncertain matrix B. In the following control development and stability analysis, it will be assumed that the matrix B is uncertain, and the robust control law will be designed with a constant feedforward estimate of the uncertain matrix. The simulation results demonstrate the capability of the robust control law to compensate for the input matrix uncertainty without the need for online parameter estimation or function approximators.
To quantify the control objective, a regulation error e 1 t R 2 and auxiliary tracking error variables e 2 t , r t R 2 are defined as
e 1 = p p d
e 2 = e ˙ 1 + α 1 e 1
r = e ˙ 2 + α 2 e 2
where α 1 , α 2 > 0 R + are user-defined control gains, and the desired plunging and pitching states p d = h , a T = 0 , 0 T for the plunging and pitching suppression objective. To facilitate the following analysis, Equation (21) is premultiplied by M and the time derivative is calculated as
M r ˙ = M e ¨ 2 + α 2 M e ˙ 2
After using Equations (18)–(21), the open-loop error dynamics are obtained as
M r ˙ = N ˜ + N d + B u ˙ e 2
where the unknown, unmeasurable auxiliary functions N ˜ e 1 , e 2 , r , N d p d , p d R 2 are defined as
N ˜ g ˙ ( p , η ) g ˙ p d , η + α 1 r α 2 e 2 α 1 e 2 + α 1 2 e 1 + α 2 M r α 2 e 2 + e 2
N d p d + g ˙ p d , η
The motivation for defining the auxiliary functions in Equations (24) and (25) is based on the fact that the following inequalities can be developed:
N ˜ ρ 0 z , N d ζ N d , N ˙ d ζ N ˙ d
where ρ 0 , ζ N d , ζ N ˙ d R + are known bounding constants, and z ( t ) R 6 is defined as
z e 1 T e 2 T r T T
Based on the open-loop error dynamics in Equation (23), the control input is designed via
u ˙ = B ^ 1 ( ( k s + I 2 × 2 ) r β s g n ( e 2 ( t ) )
where k s , β R 2 × 2 denote constant, positive definite, diagonal control gain matrices, and I 2 × 2 denotes a 2 × 2 identity matrix. In Equation (28), B ^ R 2 × 2 denotes a constant, feedforward “best guess” estimate of the uncertain input gain matrix B. Note that the control input u t does not depend on the unmeasurable acceleration term r t , since Equation (28) can be directly integrated to show that u t requires measurements of e 1 t and e 2 t only.
To facilitate the following stability proof, the control gain matrix β in Equation (28) is selected to satisfy the sufficient condition
λ min β > ζ N d + 1 α 2 ζ N ˙ d
where λ min · denotes the minimum eigenvalue of the argument. After substituting Equation (28) into Equation (23), the closed-loop error dynamics are obtained as
M r ˙ = N ˜ + N d ( k s + I n × n ) r + β s g n ( e 2 ( t ) ) e 2
To reduce the complexity of the following stability analysis, it is assumed that the product B B ^ 1 is equal to identity. It can be proven that asymptotic regulation can be achieved for the case where the feedforward estimate B ^ is within some prescribed finite range of the actual matrix B. The proof including the uncertainty in B is lengthy and is omitted here for brevity. The complete proof can be found in [23,24]. The following simulation results demonstrate the performance of the controller in the presence of uncertainty in the input gain matrix B .

3.1. Stability Analysis

Theorem 1. The controller given in Equation (28) ensures asymptotic regulation of pitching and plunging displacements in the sense that
e 1 ( t ) 0 a s t
provided the control gain k s is selected sufficiently large, and β is selected according to the sufficient condition in Equation (29).
Lemma 1. To facilitate the following proof, let D R 7 be a domain containing w ( t ) = 0 , where w ( t ) R 7 is defined as
w ( t ) z T P ( t ) T
In Equation (32), the auxiliary function P ( t ) R is the generalized solution to the differential equation
P ˙ ( t ) = L t
P 0 = β e 2 0 N d T 0 e 2 0
where the auxiliary function L ( t ) R is defined as
L ( t ) = r T N d ( t ) β s g n e 2
Provided the sufficient condition in Equation (29) is satisfied, the following inequality can be obtained:
0 t L ( τ ) d τ β e 2 0 N d T 0 e 2 0
Hence, Equation (36) can be used to conclude that P ( t ) 0 .
Proof 1. (See Theorem 1) Let V ( w , t ) : D × [ 0 , ) R be defined as the nonnegative function
V ( w , t ) 1 2 e 1 T e 1 + 1 2 e 2 T e 2 + 1 2 r T M r + P
where e 1 ( t ) , e 2 ( t ) , and r ( t ) are defined in Equations (19)–(21), respectively; and the positive definite function P ( t ) is defined in Equation (33). The function V ( w , t ) satisfies the inequality
U 1 ( w ) V ( w , t ) U 2 ( w )
provided the sufficient condition introduced in Equation (29) is satisfied, where U 1 ( w ) , U 2 ( w ) R denote the positive definite functions
U 1 λ 1 w 2 , U 2 λ 2 w 2
where λ 1 = min 1 2 , λ min M and λ 2 = max 1 , λ max M . After taking the time derivative of Equation (37) and utilizing Equation (20), Equation (21), Equation (30), and Equation (33), V ˙ ( w , t ) can be upper bounded as
V ˙ w , t α 1 1 2 e 1 2 α 2 1 2 e 2 2 r 2 k s r 2 + ρ 0 z r
where the bounds in Equation (26) were used, and the fact that e 1 T e 2 1 2 e 1 2 + 1 2 e 2 2 (i.e., Young’s inequality) was utilized. After completing the squares in Equation (40), the upper bound on V ˙ w , t can be expressed as
V ˙ w , t α 1 1 2 e 1 2 α 2 1 2 e 2 2 r 2 k s r ρ 0 2 k s z 2 + ρ 0 2 4 k s z 2
Since k s > 0 , the upper bound in Equation (41) can be expressed as
V ˙ w , t λ 0 ρ 0 2 4 k s z 2
where λ 0 min α 1 1 2 , α 2 1 2 , 1 .
The following expression can be obtained from Equation (42):
V ˙ ( w , t ) U ( w )
where U ( w ) = c z 2 , for some positive constant c R , is a continuous, positive semi-definite function.
It follows directly from the Lyapunov analysis that e 1 t , e 2 t , r t L . This implies that e ˙ 1 t , e ˙ 2 t L from the definitions given in Equations (20) and (21). Given that e ˙ 1 t , e 2 t , r t L , it follows that e ¨ 1 t L from Equation (21). Thus, Equation (19) can be used to prove that p t , p ˙ t , p ¨ t L . Since p t , p ˙ t , p ¨ t L , Equation (18) can be used to prove that u t L . Since r t , u t L , Equation (28) can be used to show that u ˙ t L . Given that e 1 t , e 2 t , r t , u ˙ t L , Equation (30) can be used along with Equation (26) to prove that r ˙ t L . Since e ˙ 1 ( t ) , e ˙ 2 ( t ) , r ˙ ( t ) L , e 1 t , e 2 t , r t , are uniformly continuous. Equation (27) can then be used to show that z t is uniformly continuous. Given that e 1 t , e 2 t , r t L , Equation (37) and Equation (42) can be used to prove that z ( t ) L L 2 . Barbalat’s lemma [25] can now be invoked to prove that z ( t ) 0 as t . Hence, e 1 ( t ) 0 as t from Equation (27). Further, given that V w , t in Equation (37) is radially unbounded, convergence of e 1 ( t ) is guaranteed, regardless of initial conditions—a global result.

4. Results and Discussion

A numerical simulation was created to demonstrate the performance of the control law developed in Equation (28). In order to develop a realistic stepping stone to high-fidelity numerical simulation results using detailed computational fluid dynamics models, the following simulation results are based on detailed dynamic parameters and specifications. The simulation is based on the dynamic model given in Equation (1) and Equation (11). The dynamic parameters utilized in the simulation are summarized in Table 1 and were obtained from [22].
Table 1. Constant simulation parameters.
Table 1. Constant simulation parameters.
ρ = 1 . 225 k g m 3 a = 0 . 24 U = 18 . 25 m s , 19 . 5 m s , 20 . 5 m s
m = 2 . 55 k g b = 0 . 11 m v = 18 m s
S α = 10 . 4 × 10 3 k g · m a 1 = 0 . 1650 a 2 = 0 . 0455
I α = 2 . 51 × 10 3 k g · m b 1 = 0 . 3350 b 2 = 0 . 3000
k h = 450 N m k α = 9 . 3 N m k α 3 = 55 N m
ζ h = 5 . 5 × 10 3 ζ α = 1 . 8 × 10 2
The following simulation results were achieved using control gains defined as
β = 0 . 001 0 0 25 , k s = 0 . 00001 0 0 0 . 11
α 1 = 1 0 0 35 , α 2 = 1 0 0 35
The control gains given in Equations (44) and (45) were selected based on achieving a desirable response in terms of settling time and required control effort. To test the case where the input gain matrix B is uncertain, it is assumed in the simulation that the actual value of B is the 2 × 2 identity matrix, but the constant feedforward estimate B ^ used in the control law is given by
B ^ = 0 . 9 0 . 1 0 . 1 1 . 1
Figure 2, Figure 3 and Figure 4, Figure 5 show the open-loop pitching and plunging response of the simulated system for the case where the free stream velocity U is 18.25 m / s and 19.5 m / s , respectively. These figures demonstrate that sustained LCO occur in the absence of feedback control. Figure 6 and Figure 7 show the closed-loop pitching and plunging response for the case where U = 18 . 25 m / s , and Figure 8 and Figure 9 show the required control force and moment, respectively, in this case. These simulation results demonstrate the capability of the robust control law to asymptotically suppress LCO using a control force and moment that are within reasonable limits. This demonstrates the practical efficacy of the control design to suppress LCO using the limited control authority that is available from the SJA. It should be noted, however, that in actual application, arrays consisting of several SJAs can be utilized to realize a greater overall control force and moment.
Figure 2. Uncontrolled pitching response for free stream velocity U = 18 . 25 m / s .
Figure 2. Uncontrolled pitching response for free stream velocity U = 18 . 25 m / s .
Robotics 03 00330 g002
Figure 3. Uncontrolled plunging response for free stream velocity U = 18 . 25 m / s .
Figure 3. Uncontrolled plunging response for free stream velocity U = 18 . 25 m / s .
Robotics 03 00330 g003
Figure 10, Figure 11 and Figure 12, Figure 13 show the closed-loop pitching and plunging response and the required control force and moment, respectively, for the case where the free stream velocity U = 19 . 5 m / s . Figure 14, Figure 15 and Figure 16, Figure 17 show the closed-loop pitching and plunging response and the required control moment and force for the case where U = 20 . 5 m / s . As shown in the figures, an increased control effort is required to suppress the LCO in the presence of the increased air velocity.
Figure 18 and Figure 19 show the robustness of the proposed controller in the event of sudden unmodeled disturbances. In these plots, a sudden increase (i.e., step function) in the pitching and plunging rates was programmed to occur at 10 . 9 s . Simultaneously, the free stream velocity U was increased from 19 m / s to 25 m / s for approximately 0.2 seconds to simulate a disturbance due to a wind gust. The pitching and plunging responses in Figure 18 and Figure 19 demonstrate the capability of the control law to compensate for the sudden and unexpected disturbances. Figure 20 and Figure 21 show the commanded control force and moment, respectively, of the closed-loop system in response to the disturbance.
Figure 4. Uncontrolled pitching response for free stream velocity U = 19 . 5 m / s .
Figure 4. Uncontrolled pitching response for free stream velocity U = 19 . 5 m / s .
Robotics 03 00330 g004
Figure 5. Uncontrolled plunging response for free stream velocity U = 19 . 5 m / s .
Figure 5. Uncontrolled plunging response for free stream velocity U = 19 . 5 m / s .
Robotics 03 00330 g005
Figure 6. Closed-loop pitching response for free stream velocity U = 18 . 25 m / s .
Figure 6. Closed-loop pitching response for free stream velocity U = 18 . 25 m / s .
Robotics 03 00330 g006
Figure 7. Closed-loop plunging response for free stream velocity U = 18 . 25 m / s .
Figure 7. Closed-loop plunging response for free stream velocity U = 18 . 25 m / s .
Robotics 03 00330 g007
Figure 8. Actuator force required for free stream velocity U = 18 . 25 m / s .
Figure 8. Actuator force required for free stream velocity U = 18 . 25 m / s .
Robotics 03 00330 g008
Figure 9. Actuator moment required for free stream velocity U = 18 . 25 m / s .
Figure 9. Actuator moment required for free stream velocity U = 18 . 25 m / s .
Robotics 03 00330 g009
Figure 10. Closed-loop pitching response for free stream velocity U = 19 . 5 m / s .
Figure 10. Closed-loop pitching response for free stream velocity U = 19 . 5 m / s .
Robotics 03 00330 g010
Figure 11. Closed-loop plunging response required for free stream velocity U = 19 . 5 m / s .
Figure 11. Closed-loop plunging response required for free stream velocity U = 19 . 5 m / s .
Robotics 03 00330 g011
Figure 12. Actuator force required for free stream velocity U = 19 . 5 m / s .
Figure 12. Actuator force required for free stream velocity U = 19 . 5 m / s .
Robotics 03 00330 g012
Figure 13. Actuator moment required for free stream velocity U = 19 . 5 m / s .
Figure 13. Actuator moment required for free stream velocity U = 19 . 5 m / s .
Robotics 03 00330 g013
Figure 14. Closed-loop pitching response for free stream velocity U = 20 . 5 m / s .
Figure 14. Closed-loop pitching response for free stream velocity U = 20 . 5 m / s .
Robotics 03 00330 g014
Figure 15. Closed-loop plunging response for free stream velocity U = 20 . 5 m / s .
Figure 15. Closed-loop plunging response for free stream velocity U = 20 . 5 m / s .
Robotics 03 00330 g015
Figure 16. Actuator force required for free stream velocity U = 20 . 5 m / s .
Figure 16. Actuator force required for free stream velocity U = 20 . 5 m / s .
Robotics 03 00330 g016
Figure 17. Actuator moment required for free stream velocity U = 20 . 5 m / s .
Figure 17. Actuator moment required for free stream velocity U = 20 . 5 m / s .
Robotics 03 00330 g017
Figure 18. Closed-loop pitching response in the presence of a disturbance.
Figure 18. Closed-loop pitching response in the presence of a disturbance.
Robotics 03 00330 g018
Figure 19. Closed-loop plunging response in the presence of a disturbance.
Figure 19. Closed-loop plunging response in the presence of a disturbance.
Robotics 03 00330 g019
Figure 20. Actuator force use in the presence of a disturbance.
Figure 20. Actuator force use in the presence of a disturbance.
Robotics 03 00330 g020
Figure 21. Actuator moment response in the presence of a disturbance.
Figure 21. Actuator moment response in the presence of a disturbance.
Robotics 03 00330 g021

5. Conclusions

A nonlinear robust control law for SJA-based LCO suppression in UAV wings is presented. The proposed control law is rigorously proven to achieve global asymptotic regulation of the pitching and plunging displacements to zero in the presence of dynamic model uncertainty and parametric actuator uncertainty. Furthermore, the proposed control law is shown via numerical simulation to compensate for unmodeled external disturbances (i.e., due to wind gusts and unmodeled effects). It is further shown that the robust control law can achieve suppression of LCO using minimal control effort. Future work will address LCO suppression control design, which can be proven to achieve simultaneous pitching and plunging regulation using only a scalar control input (i.e., the underactuated control problem). In addition, future work will target controller validation via high-fidelity computational fluid dynamics simulations.

Acknowledgments

The authors would like to thank the Air Force Research Laboratory Munitions Directorate and the American Society for Engineering Education for their support in this research.

Author Contributions

Natalie Ramos Pedroza performed the literature review and prepared the preliminary version of the article; she also prepared the original controller derivation section and assisted with the simulation code-writing. William MacKunis thoroughly checked the mathematical development, and generated most of the numerical simulation results. Vladimir Golubev proofread the article, assisted with numerical simulation results, and performed preliminary testing to determine approximations of the force and moment magnitudes that can be delivered by synthetic jet actuators.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. O’Donnell, K.S.; Marzocca, P.; Milanese, A. Design of a Wind Tunnel Apparatus to Assist Flow and Aeroelastic Control via Zero Net Mass Flow Actuators. In Proceedings of the Structures, Structural Dynamics, and Materials Conference, Fort Worth, TX, USA, 26 April 2007; pp. 1–12.
  2. Satak, N.; Hernandez, E.A.P.; Hurtado, J.E. Rate-Free Control of Nonlinear Wing Section Using Pitch and Plunge Measurements Only. In Proceedings of the 53rd AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics and Materials Conference, Honolulu, HI, USA, 23–26 April 2012; pp. 1–7.
  3. Rubillo, C.; Marzocca, P.; Bollt, E. Active aeroelastic control of lifting surfaces via jet reaction limiter control. Int. J. Bifurc. Chaos 2006, 16, 2559–2574. [Google Scholar] [CrossRef]
  4. Frampton, K.D.; Clark, R.L.; Carolina, N. Experiments on Control of Limit-Cycle Oscillations in a Typical Section Introduction. J. Guid. Control Dyn. 2000, 23, 956–960. [Google Scholar] [CrossRef]
  5. Strganac, T.W.; Kot, J.; Thompson, D.E.; Kurdila, A.J. Identifaction and control of limit cycle oscillations in aeroelastic systems. J. Guid. Control Dyn. 1999, 23, 1127–1133. [Google Scholar] [CrossRef]
  6. Platanitis, G.; Strganac, T.W. Control of a Nonlinear Wing Section Using Leading- and Trailing-Edge Surfaces. J. Guid. Control Dyn. 2004, 27, 52–58. [Google Scholar] [CrossRef]
  7. Mackunis, W.; Subramanian, S.; Mehta, S.; Ton, C.; Curtis, J.W.; Reyhanoglu, M. Robust Nonlinear Aircraft Tracking Control Using Synthetic Jet Actuators. In Proceedings of the IEEE American Control Conference, Firenze, Italy, 10–13 December 2013; pp. 1–8.
  8. Jee, S.; Lopez Mejia, O.D.; Moser, R.D.; Muse, J.A.; Kutay, A.T.; Calise, A.J. Simulation of Rapidly Maneuvering Airfoils with Synthetic Jet Actuators. AIAA J. 2013, 51, 1883–1897. [Google Scholar] [CrossRef]
  9. Deb, D.; Burkholder, J.O.; Smith, D. An Adaptive Inverse Control Scheme for Synthetic Jet Actuator Arrays. In Proceedings of Infotech@Aerospace’ 05, Arlington, VA, USA, 27 September 2005; pp. 1–11.
  10. Milanese, A.; Breuker, R.D.; Marzocca, P.; Abdalla, M. Distributed Synthetic Jet Actuators for the Control of Nonlinear Aeroelastic Systems. In Proceedings of the 19th International Conference on Adaptive Structures and Technologies, Ascona, Switzerland, 6–9 October 2008; p. 12.
  11. Amitay, M.; Smith, D.R.; Kilbens, V.; Parekh, D.E.; Glezer, A. Aerodynamic Flow Control over an Unconventional Airfoil Using Synthetic Jet Actuators. AIAA J. 2001, 39, 361–370. [Google Scholar] [CrossRef]
  12. Tchieu, A.A.; Kutay, A.T.; Muse, J.A.; Calise, A.J.; Leonard, A. Validation of a Low-Order Model for Closed-Loop Flow Control Enable Flight. In Proceedings of the 4th Flow Control Conference, Seattle, WA, USA, 23–26 June 2008; pp. 1–20.
  13. Mondschein, S.T.; Tao, G.; Burkholder, J.O. Adaptive Actuator Nonlinearity Compensation and Disturbance Rejection Applied to Aircraft Models with Synthetic Jet Actuators. In Proceedings of the AIAA Guidance, Navigation, and Control Conference, Toronto, ON, Canada, 2–5 August 2010.
  14. Deb, D.; Tao, G.; Burkholder, J.O.; Smith, D.R. An Adaptive Inverse Control Scheme for A Synthetic Jet Actuator Model. In Proceedings of the 2005 American Control Conference, Portland, OR, USA, 8–10 June 2005; Volume 82071, pp. 2646–2651.
  15. Deb, D.; Burkholder, J.; Smith, D. Adaptive synthetic jet actuator compensation for a nonlinear tailless aircraft model at low angles of attack. In Proceedings of the 2006 American Control Conference, Minneapolis, MN, USA, 14–16 June 2006; p. 6.
  16. Deb, D.; Tao, G.; Burkholder, J.O.; Smith, D.R. Adaptive Compensation Control of Synthetic Jet Actuator Arrays for Airfoil Virtual Shaping. J. Aircr. 2007, 44, 616–626. [Google Scholar] [CrossRef]
  17. Deb, D.; Tao, G.; Burkholder, J.O.; Smith, D.R. Adaptive Synthetic Jet Actuator Compensation for A Nonlinear Aircraft Model at Low Angles of Attack. IEEE Trans. Control Syst. Technol. 2008, 16, 983–995. [Google Scholar] [CrossRef]
  18. Liu, Y.; Ciuryla, M.; Amitay, M.; Chiman, K.; Myatt, J.H.; Zhang, X.; Ren, Z.; Casey, J.P. Integrated flight control and flow control using synthetic jet arrays. In Proceedings of AIAA Guidance, Navigation, and Control Conference and Exhibit, Keystone, CO, USA, 21 August 2006.
  19. Singhal, C.; Tao, G.; Burkholder, J.O. Neural Network-Based Compensation of Synthetic Jet Actuator Nonlinearities for Aircraft Flight Control. In Proceedings of the AIAA Guidance, Navigation, and Control Conference, Chicago, IL, USA, 10–13 August 2009; pp. 1–19.
  20. Tao, G.; Kokotovic, P.V. Adaptive Control of Systems with Actuator and Sensor Nonlinearities, 1st ed.; Wiley: New York, NY, USA, 1996; p. 294. [Google Scholar]
  21. Jee, S.K.; Lopez, O.; Moser, R.D.; Kutay, A.T.; Muse, J.A.; Calise, A.J. Flow Simulation of a Controlled Airfoil with Synthetic Jet Actuators. In Proceedings of the 19th AIAA Computational Fluid Dynamics, San Antonio, TX, USA, 22–25 June 2009; pp. 1–13.
  22. Golubev, V.V.; Dreyer, B.D.; Hollenshade, T.M.; Visbal, M.R. High-Accuracy Viscous Simulations of Gust-Airfoil Nonlinear Aeroelastic Interaction. In Proceedings of the 39th AIAA Fluid Dynamics Conference, San Antonio, TX, USA, 22–25 June 2009; pp. 1–15.
  23. MacKunis, W.; Wilcox, Z.; Kaiser, M.K.; Dixon, W.E. Global Adaptive Output Feedback Tracking Control of an Unmanned Aerial Vehicle. IEEE Trans. Control Syst. Technol. 2010, 18, 1390–1397. [Google Scholar] [CrossRef]
  24. MacKunis, W.; Wilcox, Z.D.; Kaiser, M.K.; Dixon, W.E. Global Adaptive Output Feedback MRAC. In Proceedings of the IEEE Confernece on Decision and Control, Shanghai, China, 15–18 December 2009; pp. 3483–3488.
  25. Khalil, H.K. Nonlinear Systems, 3rd ed.; Prentice Hall: Upper Saddle River, NJ, USA, 2002; p. 750. [Google Scholar]

Share and Cite

MDPI and ACS Style

Pedroza, N.R.; MacKunis, W.; Golubev, V.V. Robust Nonlinear Regulation of Limit Cycle Oscillations in UAVs Using Synthetic Jet Actuators. Robotics 2014, 3, 330-348. https://doi.org/10.3390/robotics3040330

AMA Style

Pedroza NR, MacKunis W, Golubev VV. Robust Nonlinear Regulation of Limit Cycle Oscillations in UAVs Using Synthetic Jet Actuators. Robotics. 2014; 3(4):330-348. https://doi.org/10.3390/robotics3040330

Chicago/Turabian Style

Pedroza, Natalie Ramos, William MacKunis, and Vladimir V. Golubev. 2014. "Robust Nonlinear Regulation of Limit Cycle Oscillations in UAVs Using Synthetic Jet Actuators" Robotics 3, no. 4: 330-348. https://doi.org/10.3390/robotics3040330

APA Style

Pedroza, N. R., MacKunis, W., & Golubev, V. V. (2014). Robust Nonlinear Regulation of Limit Cycle Oscillations in UAVs Using Synthetic Jet Actuators. Robotics, 3(4), 330-348. https://doi.org/10.3390/robotics3040330

Article Metrics

Back to TopTop