Nomenclature
The nomenclature of this study is based upon the following principles and introduced here in
concise form. A (relative) position vector is denoted as follows:
Here, A and B are the points, while C is the coordinate frame. A velocity is consequently given by
Here, D denotes the kind of the velocity, e.g., kinematic K. Furthermore, the superscript E describes
the coordinate frame with respect to whom the time derivative is taken. Then, an acceleration is
The superscript F depicts the coordinate frame with respect to whom the second time derivative
is taken.
After this overview on vectors, the following listing summarizes further symbols and indices
used in this study.
| Linear acceleration vector |
| Trajectory acceleration deviation in three axes |
| Rotation acceleration vector |
| Kinematic angle of attack |
| Kinematic trajectory angles |
| Specific force vector |
| Specific forces in three axes |
| Actuator transfer function |
g | Gravitational acceleration |
| Jerk vector |
| Controller gain matrix |
m | Mass |
| Nonlinear dynamic inversion pseudocontrol |
| Rotational velocity vector |
| Natural frequency |
| Euler angles |
| Position vector |
| Trajectory deviation in three axes |
| Aircraft position in navigation frame (generally: North-East-Down) |
s | Laplace transform variable |
| Nonlinear dynamic inversion control vector |
| Trajectory velocity deviation in three axes |
| Kinematic velocity |
| Linear velocity vector |
| Relative damping |
| Lift force derivative w.r.t. angle of attack |
F | Trajectory foot point |
K | Kinematic quantity/frame |
N | Navigation/inertial frame |
O | Aircraft-/Trajectory-fixed North-East-down Frame |
R | Aircraft reference point |
T | Trajectory frame |
1. Introduction
Trajectory control as well as path-following are important topics in the research and industrial communities among others in, e.g., the fields of mobile robotics [
2] and autonomous driving [
3]. Research in these areas is mainly concerned with both a high availability of the path generation as well as a high integrity, stability, and accuracy of the trajectory/path-following controller. Especially, the importance of stable trajectory/path-following control is central as it ensures the safety of, e.g., the vehicle’s and robot’s operation.
As a consequence of the increased onboard computation power and the emergence of small aerial vehicles, which should operate in close proximity to, e.g., housing, high accurate path control became of paramount importance in the context of aviation in the previous years as well. An overview on often used algorithms is given in [
4,
5]: The variety of the applied algorithms ranges from simple geometric approaches to linear and decoupled control law design as well as sophisticated nonlinear control techniques, such as dynamic inversion and backstepping, and even optimal control-based approaches like model predictive control.
Especially, nonlinear control techniques have shown their viability in application as they are able to exploit the full physical capabilities and available control authority of an aircraft, thus ensuring a high availability, integrity, and accuracy of the controller command. Additionally, they provide good interaction with the necessary inner-loop controllers as they can consider their dynamic behavior. Examples comprise backstepping control approaches [
6], (nonlinear) adaptive control [
7], and nonlinear dynamic inversion [
8,
9,
10].
Considering the application to aircraft trajectory/path-following, the approach presented in [
8,
10] has proven to be viable due to its modularity and independence of the dynamic model. Additionally, it allows for an independent control of the rigid-body/rotational dynamics in an inner loop by a stabilizing controller. Still, a drawback of the methodology is the fact that it is based on acceleration-level aircraft dynamics and kinematics, and thus produces force commands that the inner loop must directly follow. While this is theoretically possible, the inner-loop dynamics and actuators do not allow for direct changes of forces. Thus, this study proposes an enhancement of the existing trajectory/path-following approach by using jerk instead of acceleration-based dynamics and kinematics. In doing so, the command to the inner loop is a force derivative instead of the force, which allows for consideration of actuator dynamics [
11,
12] and also yields smoother force commands by construction. A further advantage of the jerk-based control for manned application is that studies suggest an increased passenger comfort by controlling the jerk [
13]. In addition, jerk is also a major contributor to the efficiency and this therefore often a minimization objective in robot manipulator path planning [
14].
The general idea of the control approach of the paper is depicted in
Figure 1: The desired trajectory is generated by a suitable method, e.g., geometric approaches [
15,
16], (robust) optimal control [
17,
18,
19,
20], or shortest path methods [
21]. As the aircraft is generally never able to follow the path exactly, either due to imperfect planning or disturbances (e.g, wind or model errors), a deviation between the aircraft reference point
and the trajectory foot point
arises. This point is either specified by a “virtual” aircraft that the actual aircraft “chases” (this is said to be “trajectory control” in this study) [
5] or by the orthogonal projection of the aircraft reference point on the trajectory (this is said to be “path-following control” in this study), which may be calculated by the methods proposed in [
15,
22]. The trajectory/path-following controller then uses these deviations as well as their derivatives (velocity and acceleration) to calculate a suitable control command to reduce them. It should be noted that already by considering acceleration errors, the controller can react faster to deviation than with the acceleration based formulation in [
8,
10], which only starts acting on velocity errors.
Summarizing, the contributions of this study are as follows.
Derivation of jerk-level dynamic equations, i.e., a connection between the jerk and the (specific) force derivatives.
Derivation of nonlinear, jerk-based error dynamics between the reference point of the dynamic model and the desired point (“foot point”) on the trajectory.
Connection of jerk-based dynamics and kinematics to formulate control law based on nonlinear dynamic inversion.
To achieve this, the study is structured as follows. In
Section 2, the kinematic equations of motion are extended to the jerk-level and force derivative. This is used in
Section 3 together with the jerk-level error dynamics and nonlinear dynamic inversion to formulate a nonlinear control law. Then, the viability of the approach is illustrated by the example in
Section 4. Concluding,
Section 5 gives remarks and an outlook on future work.
2. Jerk-Based Dynamic Equations
The derivation of the jerk-based dynamics starts with the acceleration-based kinematic equation of motion using specific forces on a fixed-flat earth (Assumption A.3) [
1,
10]:
Here, the specific forces are normalized by the mass of the dynamic system (e.g., aircraft or robot). It should be noted that this is possible due to Assumption A.2, i.e., because the model mass is time-invariant.
To get to the jerk-level and introduce the time derivative of the specific forces, Equation (
1) must be differentiated with respect to time in an inertial frame, which is the navigation frame
in this study (Assumptions A.2 and A.3):
Using the Euler differentiation rule [
1] and introducing the jerk for the time derivative of the acceleration, Equation (
2) results in the following form.
Consequently, Equation (
3) provides the desired connection between the jerk
and the specific force derivative
that is utilized in
Section 3 to the dynamic equation with the nonlinear error dynamics.
3. Nonlinear Error Dynamics and Control Strategy
The error dynamics are derived in the following for the difference between the orthogonal projection of the trajectory foot point and the dynamic model reference point (see
Figure 1) [
8,
10]:
Based on Equation (
4), the velocity deviation can be calculated by deriving Equation (
4) with respect to the trajectory frame and applying the Euler differentiation rule (it should be noted that we do not need to derive with respect to an inertial frame in this case, as the position deviation is not based on Newtonian mechanics. Furthermore, the velocity is specified kinematically, and thus has an additional index
added inside the bracket) [
8,
10]:
Then, the acceleration deviation is consequently given by deriving Equation (
5) once more in the trajectory frame as well as applying the Euler differentiation rule to introduce measurable quantities [
8,
10]:
Here, the lead acceleration
is split (see Equation (
4)) into the a part with respect to the reference point,
, and a part with respect to the foot point (the “desired” acceleration),
. This is used in the following to introduce the dynamic equations from
Section 2 into the jerk dynamics derived next. Additionally, the desired acceleration/jerk is an additional design variable in the commanded trajectory specification that can be used to achieve a desired behavior.
Finally, the jerk error can be calculated by deriving Equation (
6) one final time with respect to the trajectory frame:
Now, the Euler differentiation rule may be applied once more to transform the derivatives with respect to the trajectory frame to a more suitable representation for analysis and implementation:
Thus, the nonlinear jerk-based error dynamics have been derived. It is important to note that Equation (
8) shows that the error dynamics are, generally, a linear parameter varying system with respect to the trajectory deviation and its derivatives (last line in Equation (
8)). Additionally, the top two lines in Equation (
8) depict the nonlinear influence of the dynamic model, while the desired jerk of the foot point (second line) is a design variable of the planned trajectory.
As stated in the sections before, the goal of the trajectory controller is to command the specific force derivatives to the inner-loop controller, which stabilizes the rotational dynamics, utilizing the nonlinear error dynamics in Equation (
8). To achieve this, Equation (
3) can be solved for the jerk
(including a transformation from the kinematic
to the trajectory
frame using an orthogonal transformation matrix) and inserted for the first addend in Equation (
8):
Here, nonlinear dynamic inversion [
1,
8,
10], i.e.,
with
containing all nonlinear addends from Equation (
8) being used to calculate the specific force derivative (command variable
), which may then be commanded to the inner loop. It should be noted that the values in
are all available from sensor measurements (e.g., aircraft linear/rotational velocity and acceleration) or can be easily calculated from these (e.g., deviations). For this, all desired (i.e., planned) trajectory related quantities are specified by, e.g., the trajectory generation module and consequently, also known. Thus, Equation (
10) can be used to calculate the desired control
(i.e., the specific force derivative) by solving for it and using the pseudocontrol
as an input variable. For this purpose, the pseudocontrol may, e.g., follow a linear error controller law:
Consequently, the pseudocontrol is chosen such that the error between reference and foot point reduces, while following the desired trajectory. The matrices , , and are positive definite design matrices specified to shape the response characteristics and to ensure the stability of the error dynamics by commanding a stabilizing pseudocontrol. Thus, a proper choice of these matrices ensures not only the stability of the trajectory controller, but also provides disturbance rejection.
It should be noted that the nonlinear error controller in Equations (
9)–(
11) is derived for the case of trajectory control here. Nonetheless, it is straightforward to apply the methodology to path-following control as well by removing the dependence on
(i.e., not controlling this deviation).
Concluding, the control law based on nonlinear error dynamics has been derived in Equations (
9) and (
11) using nonlinear dynamic inversion as specified in Equation (
10). The applicability is shown in an aircraft-related example in the next section.
4. Illustrative Example
This section gives different test results to show the viability of the proposed control scheme in
Section 3. Therefore,
Section 4.1 starts with an introduction of the generic simulation model. Then,
Section 4.2 investigates the effect of the turn radius on the maximal lateral trajectory deviation, i.e., the effects of the model dynamics including actuators and how the error controller deals with it in a simple setup. Finally,
Section 4.3 concludes with an example of the controller in trajectory control mode for a complex trajectory-following task.
4.1. Model Formulation
The model is based on the kinematic equations of motion introduced in Equation (
1) with constants defined in
Table 1. To improve the model’s fidelity for real-world applications, a lower order equivalent system transfer function
is used for course and climb as well as the angle of attack dynamics to mimic a rigid-body behavior, and thus the rotational dynamics. This is defined as follows (in the
Laplace domain [
1]).
Beside the dynamic relations in Equation (
12), the following kinematic relations are used to calculate quantities required to, e.g., evaluate the error dynamics in Equation (
9) or transformation matrices:
Thus, Equations (
12) and (
13) together with the basic kinematic equations of motion in Equation (
1) and a position propagation in navigational frame specify the used generic aircraft model for the simulation assessment. The control inputs, i.e., the specific forces, are calculated by integrating the results obtained when applying Equations (
9) and (
10) (with zero specific force in
x and
y direction as well as negative gravitational acceleration
in
z direction, in order to account for the gravity and ensure a level flight, in the beginning of the simulation). It should be noted that this time-integration step reduces the efficiency of the proposed scheme in the context of this study as the specific force derivatives are not directly used as control variables. However, this formulation will interface with the inner-loop controller described in [
11,
12] that utilizes specific force derivatives as input. Nonetheless, the effectiveness of the approach is still shown in the following examples by this simplified setup.
Furthermore, the gain matrices specified in
Table 1 used for the pseudocontrol calculation in Equation (
11) are not specifically tuned and only designed to provide stability. A specific tuning can further improve the results, while, again, the simplified setup already proves the applicability of the proposed method.
Take into account that the desired trajectory, i.e., the path the simulation aircraft should follow, is specified using command signals for the desired velocity as well as climb and course rate. Finally, the desired jerk in Equation (
9) is set to zero as minimum jerk trajectories are generally most energy-efficient and best perceived by passengers [
13,
14].
4.2. Turn Radius Influence on Maximal Error
This example deals with an analysis of the turn radius influence on the maximal lateral deviation allowed by the controller. Therefore, a turn is considered, while the aircraft velocity is set to , which corresponds to a typical flight velocity of a small unmanned aerial vehicle. Initially, the aircraft’s position and the desired trajectory match.
First,
Figure 2 shows the maximum lateral deviation of the
-turn for different turn radii. It is clear that for a smaller turn radius, i.e., a tighter, more aggressive maneuver, the lateral deviation error is increased. Still, the maximum error is only around
at a turn radius of approximately
, which shows the capabilities of the proposed controller for highly accurate trajectory control. It is furthermore seen that for gentle maneuvers, i.e., an increased turn radius, the lateral deviation asymptotically approaches zero. The remaining magnitude of the error is hereby mainly related to the model lags introduced in Equation (
12). It should be further noted that it is expected that the maximum lateral deviation decreases further when not only the specific force is commanded to the aircraft but also its derivative. Then, the reaction to deviations is improved even further and the introduced lags in Equation (
12) can be directly accounted for.
Continuing,
Figure 3 displays the maximum absolute lateral specific force along the
-turn. As expected, a higher lateral force is required to follow a tighter, i.e., more aggressive, turn maneuver. This strengthens the argumentation used before in
Figure 2 that a more aggressive maneuver, yielding a larger lateral specific force requirement, results in a greater deviation. Still, as also already seen in
Figure 2, even a considerable lateral specific force of around
only yields to a tolerable lateral deviation of
.
Concluding, this application example shows the applicability of the proposed trajectory controller in scenarios with high accuracy demand. Specifically, both gentle as well as aggressive maneuvers are handled by the proposed scheme.
4.3. Four-Dimensional Trajectory Control
The application section concludes with an example of a complex trajectory following task. Here,
Figure 4 shows the desired/commanded (solid red) vs. the actually flown (dashed blue) trajectory for a complex climb and turn maneuver including a holding. Furthermore, the planar projections (dashed black) of the maneuver are shown for convenience and to get a better overview of the maneuver structure.
Then,
Figure 5 shows the position, velocity, and acceleration deviations obtained from following the trajectory in
Figure 4. It can be seen that maximum vertical errors of
occur in climb/descent segments, while the maximum lateral deviation is around
. In general, the horizontal deviation shows the overall best tracking.
This overall shows the high accuracy of the controller. Additionally, the errors are rapidly reduced, which is specifically obtained due to the error controller that already starts acting at the acceleration deviation level. In this context, it can be seen that maximum vertical acceleration deviation of and of for velocity is achieved. This shows both the fast response of the controller to deviations as well as the smoothness, i.e., that small deviations are controlled with small control effort.
Finally,
Figure 6 shows the specific forces and their derivatives commanded to the model when following the trajectory in
Figure 4. It is specifically seen that a high smoothness for the specific forces as well as the decoupling of the axes by nonlinear dynamic inversion is achieved, which results in a more physical command for controlling the aircraft. Still, it is clearly seen that the specific force derivatives directly react to the deviations in
Figure 5. Specifically, it can be observed that the specific force derivatives counteract the acceleration errors in
Figure 5. Furthermore, it is seen that there are generally no high-frequency oscillations in the command signal, which is a specific consequence of the fact that the nonlinear aircraft kinematics are considered by the controller. Due to these, it is known how the aircraft can move, and thus the controller does not try to compensate deviations by large control effort but by reasonable magnitudes and dynamics that the aircraft is able to follow.