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.
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
where the coefficients
,
denote the structural mass and damping matrices,
is a nonlinear stiffness matrix, and
denotes the state vector. In Equation (
1),
is explicitly defined as
where
,
denote the plunging [meters] and pitching [radians] displacements describing the LCO effects. Also in Equation (
1), the structural linear mass matrix
[
22]
where the parameters
are the static moment and moment of inertia, respectively. The structural linear damping matrix is described as
where the parameters
,
are the damping logarithmic decrements for plunging and pitching, and
is the mass of the wing, or in this case, a flat plate. The nonlinear stiffness matrix utilized in this study is
where
denote structural resistances to pitching (linear and nonlinear) and
is the structural resistance to plunging.
In Equation (
1), the total lift and moment are explicitly defined as
where
denote the equivalent control force and moment generated by the
jth SJA, and
,
are the aerodynamic lift and moment due to the 2 degree-of-freedom motion [
10]. In Equation (
6),
denotes the aerodynamic state vector that relates the moment and lift to the structural modes. Also in Equation (
6),
denotes the SJA-based control input (e.g., the SJA air velocity or acceleration), and
is an uncertain constant input gain matrix that relates the control input
to the equivalent force and moment generated by the SJA. Also in Equation (
6), the aerodynamic and mode matrices
,
,
,
are described as
where
is the Wagner solution function at 0, and the parameters
,
,
,
are the Wagner coefficients. In addition,
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]
The aerodynamic state matrices in Equation (11),
,
,
, are explicitly defined as
By substituting Equation (
6) into Equation (
1), the LCO dynamics can be expressed as
where
,
, and
. By making the definitions
,
,
,
,
, and
, the dynamic equation in Equation (15) can be expressed in state form as
where
is the state vector,
is the state matrix (state-dependent). In Equation (16), the input gain matrix
is defined as
where
denotes a
matrix of zeros. The structure of the input gain matrix in Equation (17) results from the fact that the control input
only directly affects
and
.
3. Control Development
The objective is to design the control signal
to regulate the plunge and pitching dynamics (
i.e.,
) to zero. To facilitate the control design, the expression in Equation (15) is rewritten as
where
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 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
and auxiliary tracking error variables
are defined as
where
are user-defined control gains, and the desired plunging and pitching states
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
After using Equations (18)–(21), the open-loop error dynamics are obtained as
where the unknown, unmeasurable auxiliary functions
,
are defined as
The motivation for defining the auxiliary functions in Equations (24) and (25) is based on the fact that the following inequalities can be developed:
where
,
,
are known bounding constants, and
is defined as
Based on the open-loop error dynamics in Equation (23), the control input is designed via
where
,
denote constant, positive definite, diagonal control gain matrices, and
denotes a
identity matrix. In Equation (28),
denotes a constant, feedforward “best guess” estimate of the uncertain input gain matrix
B. Note that the control input
does not depend on the unmeasurable acceleration term
, since Equation (28) can be directly integrated to show that
requires measurements of
and
only.
To facilitate the following stability proof, the control gain matrix
β in Equation (28) is selected to satisfy the sufficient condition
where
denotes the minimum eigenvalue of the argument. After substituting Equation (28) into Equation (23), the closed-loop error dynamics are obtained as
To reduce the complexity of the following stability analysis, it is assumed that the product
is equal to identity. It can be proven that asymptotic regulation can be achieved for the case where the feedforward estimate
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
3.1. Stability Analysis
Theorem 1. The controller given in Equation (28) ensures asymptotic regulation of pitching and plunging displacements in the sense thatprovided the control gain is selected sufficiently large, and β is selected according to the sufficient condition in Equation (29).
Lemma 1. To facilitate the following proof, let be a domain containing , where is defined asIn Equation (32), the auxiliary function is the generalized solution to the differential equationwhere the auxiliary function is defined asProvided the sufficient condition in Equation (29) is satisfied, the following inequality can be obtained:Hence, Equation (36) can be used to conclude that . Proof 1. (See Theorem 1)
Let be defined as the nonnegative functionwhere , , and are defined in Equations (19)–(21), respectively; and the positive definite function is defined in Equation (33). The function satisfies the inequalityprovided the sufficient condition introduced in Equation (29) is satisfied, where , denote the positive definite functionswhere and . After taking the time derivative of Equation (37) and utilizing Equation (20), Equation (21), Equation (30), and Equation (33), can be upper bounded aswhere the bounds in Equation (26) were used, and the fact that (i.e., Young’s inequality) was utilized. After completing the squares in Equation (40), the upper bound on can be expressed asSince , the upper bound in Equation (41) can be expressed aswhere . The following expression can be obtained from Equation (42):where , for some positive constant , is a continuous, positive semi-definite function. It follows directly from the Lyapunov analysis that . This implies that from the definitions given in Equations (20) and (21). Given that , , it follows that from Equation (21). Thus, Equation (19) can be used to prove that , . Since , , Equation (18) can be used to prove that . Since , , Equation (28) can be used to show that . Given that , , Equation (30) can be used along with Equation (26) to prove that . Since , , are uniformly continuous. Equation (27) can then be used to show that is uniformly continuous. Given that , Equation (37) and Equation (42) can be used to prove that . Barbalat’s lemma [25] can now be invoked to prove that as . Hence, as from Equation (27). Further, given that in Equation (37) is radially unbounded, convergence of 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.
| | |
---|
| | |
| | |
| | |
| | |
| | |
The following simulation results were achieved using control gains defined as
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
identity matrix, but the constant feedforward estimate
used in the control law is given by
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
and 19.5
, 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
, 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 .
Figure 2.
Uncontrolled pitching response for free stream velocity .
Figure 3.
Uncontrolled plunging response for free stream velocity .
Figure 3.
Uncontrolled plunging response for free stream velocity .
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
.
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
. 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
. Simultaneously, the free stream velocity U was increased from
to
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 .
Figure 4.
Uncontrolled pitching response for free stream velocity .
Figure 5.
Uncontrolled plunging response for free stream velocity .
Figure 5.
Uncontrolled plunging response for free stream velocity .
Figure 6.
Closed-loop pitching response for free stream velocity .
Figure 6.
Closed-loop pitching response for free stream velocity .
Figure 7.
Closed-loop plunging response for free stream velocity .
Figure 7.
Closed-loop plunging response for free stream velocity .
Figure 8.
Actuator force required for free stream velocity .
Figure 8.
Actuator force required for free stream velocity .
Figure 9.
Actuator moment required for free stream velocity .
Figure 9.
Actuator moment required for free stream velocity .
Figure 10.
Closed-loop pitching response for free stream velocity .
Figure 10.
Closed-loop pitching response for free stream velocity .
Figure 11.
Closed-loop plunging response required for free stream velocity .
Figure 11.
Closed-loop plunging response required for free stream velocity .
Figure 12.
Actuator force required for free stream velocity .
Figure 12.
Actuator force required for free stream velocity .
Figure 13.
Actuator moment required for free stream velocity .
Figure 13.
Actuator moment required for free stream velocity .
Figure 14.
Closed-loop pitching response for free stream velocity .
Figure 14.
Closed-loop pitching response for free stream velocity .
Figure 15.
Closed-loop plunging response for free stream velocity .
Figure 15.
Closed-loop plunging response for free stream velocity .
Figure 16.
Actuator force required for free stream velocity .
Figure 16.
Actuator force required for free stream velocity .
Figure 17.
Actuator moment required for free stream velocity .
Figure 17.
Actuator moment required for free stream velocity .
Figure 18.
Closed-loop pitching response in the presence of a disturbance.
Figure 18.
Closed-loop pitching response in the presence of a disturbance.
Figure 19.
Closed-loop plunging response in the presence of a disturbance.
Figure 19.
Closed-loop plunging response in the presence of a disturbance.
Figure 20.
Actuator force use in the presence of a disturbance.
Figure 20.
Actuator force use in the presence of a disturbance.
Figure 21.
Actuator moment response in the presence of a disturbance.
Figure 21.
Actuator moment response in the presence of a disturbance.