Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Modeling of the Minimum Fluidization Velocity and the Incipient Fluidization Pressure Drop in a Conical Fluidized Bed with Negative Pressure
Next Article in Special Issue
On Invertibility of an Interconnected System Composed of Two Dynamic Subsystems
Previous Article in Journal
Modeling, Simulation, and Cruise Characteristics of Wingtip-Jointed Composite Aircraft
Previous Article in Special Issue
Observer Design for Nonlinear Invertible System from the View of Both Local and Global Levels
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Trajectory/Path-Following Controller Based on Nonlinear Jerk-Level Error Dynamics

Institute of Flight System Dynamics, Technical University of Munich, 85748 Garching, Germany
*
Author to whom correspondence should be addressed.
Appl. Sci. 2020, 10(23), 8760; https://doi.org/10.3390/app10238760
Submission received: 20 November 2020 / Revised: 3 December 2020 / Accepted: 4 December 2020 / Published: 7 December 2020
(This article belongs to the Special Issue Nonlinear Control: Theory and Applications)

Abstract

:
This study proposes a novel, nonlinear trajectory/path-following controller based on jerk-level error dynamics. Therefore, at first the nonlinear acceleration-based kinematic equations of motion of a dynamic system are differentiated with respect to time to obtain a representation connecting the translation jerk with the (specific) force derivative. Furthermore, the path deviation, i.e., the difference between the planned and the actual path, is formulated as nonlinear error dynamics based on the jerk error. Combining the derived equations of motion with the nonlinear error dynamics as well as employing nonlinear dynamic inversion, a control law can be derived that provides force derivative commands, which may be commanded to an inner loop for trajectory control. This command ensures an increased smoothness and faster reaction time compared to traditional approaches based on a force directly. Furthermore, the nonlinear parts in the error dynamic are feedforward components that improve the general performance due to their physical connection with the real dynamics. The validity and performance of the proposed trajectory/path-following controller are shown in an aircraft-related application example.

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:
r A B C = r B C r A C
Here, A and B are the points, while C is the coordinate frame. A velocity is consequently given by
d d t E r A B C = r ˙ D A B C E = v D A B C E
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
d d t F v D A B C E = v ˙ D A B C E F = a D A B C E F
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.
a Linear acceleration vector
Δ a x , T , Δ a y , T , Δ a z , T Trajectory acceleration deviation in three axes
α Rotation acceleration vector
α K Kinematic angle of attack
μ K , γ K , χ K Kinematic trajectory angles
Δ f Specific force vector
Δ f x , Δ f y , Δ f z Specific forces in three axes
G LOES s Actuator transfer function
gGravitational acceleration
𝑗 Jerk vector
K Controller gain matrix
mMass
ν Nonlinear dynamic inversion pseudocontrol
ω Rotational velocity vector
ω 0 Natural frequency
ϕ , θ , ψ Euler angles
r Position vector
Δ x T , Δ y T , Δ z T Trajectory deviation in three axes
x N , y N , z N Aircraft position in navigation frame (generally: North-East-Down)
sLaplace transform variable
u Nonlinear dynamic inversion control vector
Δ u T , Δ v T , Δ w T Trajectory velocity deviation in three axes
V K Kinematic velocity
v Linear velocity vector
ζ Relative damping
Z α Lift force derivative w.r.t. angle of attack
FTrajectory foot point
KKinematic quantity/frame
NNavigation/inertial frame
OAircraft-/Trajectory-fixed North-East-down Frame
RAircraft reference point
TTrajectory frame

Assumptions

The following listing summarizes the assumptions used throughout this study.
A. 1
Dynamic model is a rigid-body, i.e., there is no relative movement/velocity between two points on the body.
A. 2
Quasi-stationary mass distribution, i.e., m ˙ 0 .
A. 3
Reference frame is a navigation frame N (e.g., North-East-Down) on a fixed-flat earth (i.e., no influence of earth rotation and ellipsoidal earth shape are considered).
A. 4
Reference point is the center of gravity.
It should be noted that these assumptions are fairly standard in, e.g., aviation or robotics [1], and are only used in this paper for the sake of simplicity. Indeed, the specified assumptions are generally not required for the derivation as well as for the application of the proposed methodology and can thus be relaxed based on the required model fidelity. Thus, they do not imply any loss of generality.

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 R and the trajectory foot point F 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]:
a K R K N K + ω K O K K × v K R K N = Δ f R K
Here, the specific forces Δ f R K = Δ f x Δ f y Δ f z T are normalized by the mass m 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 N in this study (Assumptions A.2 and A.3):
a ˙ K R K N K N 𝑗 K R K N K N + α K O K K N ω ˙ K O K K N × v K R K N + ω K O K K × a K R K N N = Δ f ˙ R K N
Using the Euler differentiation rule [1] and introducing the jerk for the time derivative of the acceleration, Equation (2) results in the following form.
𝑗 K R K N K K + ω K N K K × a K R K N K 𝑗 K R K N K N + α K O K K K + ω K N K K × ω K O K K α K O K K N × v K R K N + ω K O K K × a K R K N K + ω K N K K × v K R K N a K R K N N = Δ f ˙ R K K + ω K N K K × Δ f R K Δ f ˙ R K N
Consequently, Equation (3) provides the desired connection between the jerk 𝑗 K R K N K K and the specific force derivative Δ f ˙ R K K 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]:
r F R T = r R T r F T = Δ x T Δ y T Δ z T T
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 K added inside the bracket) [8,10]:
d d t T r F R T = v K F R T T = Δ u T Δ v T Δ w T T = v K F R T N + ω K T N T × r F R T
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]:
a K F R T T T = Δ a x , T Δ a y , T Δ a z , T T = a K R T N K + ω K T K T × v K R T N a K R T N T a des , K F T N T a K F R T N T + ω K T N T × v K F R T T + α K T N T T × r F R T
Here, the lead acceleration a K F R T N T is split (see Equation (4)) into the a part with respect to the reference point, a K R T N T , and a part with respect to the foot point (the “desired” acceleration), a des , K F T N T . 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:
𝑗 K F R T T T T = 𝑗 K R T N K T + α K T K T T × v K R T N + ω K T K T × a K R T N T 𝑗 des , K F T N T T + α K T N T T × v K F R T T + ω K T N T × a K F R T T T + α ˙ K T N T T T × r F R T + α K T N T T × v K F R T T
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:
𝑗 K F R T T T T = 𝑗 K R T N K K + ω K T K T × a K R T N K 𝑗 K R T N K T + α K T K T T × v K R T N + ω K T K T × a K R T N K + ω K T K T × v K R T N a K R T N T 𝑗 des , K F T N T T + ω K T N T × a K F R T T T + 2 α K T N T T × v K F R T T + α ˙ K T N T T T × r F R T
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 𝑗 K R T N K K (including a transformation from the kinematic K to the trajectory T frame using an orthogonal transformation matrix) and inserted for the first addend in Equation (8):
𝑗 K F R T T T T ν = Δ f ˙ R T K u + ω K N K T × Δ f R T ω K N K T × a K R T N K α K O K T K + ω K N K T × ω K O K T × v K R T N ω K O K T × a K R T N K + ω K N K T × v K R T N + ω K T K T × a K R T N K + α K T K T T × v K R T N + ω K T K T × a K R T N K + ω K T K T × v K R T N 𝑗 des , K F T N T T + ω K T N T × a K F R T T T + 2 α K T N T T × v K F R T T + α ˙ K T N T T T × r F R T
Here, nonlinear dynamic inversion [1,8,10], i.e.,
u = ν g · ,
with g · containing all nonlinear addends from Equation (8) being used to calculate the specific force derivative (command variable u ), which may then be commanded to the inner loop. It should be noted that the values in g · 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 u (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:
ν = 𝑗 K F R T T T T = K p · r F R T K d · v K F R T T K dd · a K F R T T T , K p 0 , K d 0 , K d d 0
Consequently, the pseudocontrol is chosen such that the error between reference and foot point reduces, while following the desired trajectory. The matrices K p , K d , and K d d 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 Δ x T (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 G LOES 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]).
G LOES s = ω 0 2 s 2 + 2 ζ ω 0 s + ω 0 2 , γ K ˙ , cmd G LOES s γ K ˙ , χ K ˙ , cmd G LOES s χ K ˙ , m Δ f z Z α G LOES s α K ˙
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:
ϕ μ K = arctan V K g χ K ˙ , θ α K + γ K , ψ ˙ χ K ˙
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 g 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 180 turn is considered, while the aircraft velocity is set to V K = 22 m / s , 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 180 -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 0.1 m at a turn radius of approximately 5 m , 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 180 -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 9 m / s 2 only yields to a tolerable lateral deviation of 0.1 m .
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 ± 0.02 m occur in climb/descent segments, while the maximum lateral deviation is around ± 0.01 m . 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 ± 0.05 m / s 2 and of ± 0.01 m / s 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.

5. Conclusions

A novel trajectory controller utilizing full nonlinear kinematics on the jerk-level between the reference point of the dynamic model and the planned trajectory foot point was developed in this study. The resulting nonlinear error dynamics was controlled using a nonlinear dynamic inversion control law that connected jerk-based nonlinear error dynamics with kinematic equations of motion. The proposed scheme was applied in an aircraft-related context and shows promising results and significant capabilities in controlling the trajectory/path with a high accuracy.
Further work naturally starts with an incorporation of the trajectory controller into a real flying vehicle to show the capabilities in real flight test. Additionally, surveys on choosing suitable commands other than zero for the desired jerk command should be conducted to achieve, e.g., optimal minimum jerk trajectories. Finally, the gain tuning within the pseudocontrol in Equation (11) should be addressed. To this end, it is especially interesting to consider the nature of the derived nonlinear error dynamics in Equation (9), which are of a linear, parameter-varying type. Here, sophisticated analysis methods can be used to significantly boost performance and already ensure stability as well as desired behavior of the error dynamics in connection with a suitable design of the reference trajectory.

Author Contributions

Conceptualization, P.P., M.M.M. and P.B.; methodology, P.P., M.M.M. and P.B.; software, M.M.M.; validation, P.P. and P.B.; formal analysis, M.M.; investigation, P.P., M.M.M. and P.B.; resources, F.H.; writing—original draft preparation, P.P.; writing—review and editing, M.M.M., P.B. and F.H.; visualization, P.P.; supervision, F.H.; project administration, F.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The authors want to acknowledge the Technical University of Munich (TUM) for supporting the publication in the framework of the TUM Open Access Publication Fund.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Stevens, B.L.; Lewis, F.L.; Johnson, E.N. Aircraft Control and Simulation: Dynamics, Controls Design, and Autonomous Systems, 3rd ed.; Wiley: Hoboken, NJ, USA, 2016. [Google Scholar] [CrossRef]
  2. Elbanhawi, M.; Simic, M. Sampling-Based Robot Motion Planning: A Review. IEEE Access 2014, 2, 56–77. [Google Scholar] [CrossRef]
  3. Claussmann, L.; Revilloud, M.; Gruyer, D.; Glaser, S. A Review of Motion Planning for Highway Autonomous Driving. IEEE Trans. Intell. Transp. Syst. 2020, 21, 1826–1848. [Google Scholar] [CrossRef] [Green Version]
  4. Sujit, P.B.; Saripalli, S.; Sousa, J.B. Unmanned Aerial Vehicle Path Following: A Survey and Analysis of Algorithms for Fixed-Wing Unmanned Aerial Vehicless. IEEE Control Syst. Mag. 2014, 34, 42–59. [Google Scholar] [CrossRef]
  5. Rubí, B.; Pérez, R.; Morcego, B. A Survey of Path Following Control Strategies for UAVs Focused on Quadrotors. J. Intell. Robot. Syst. 2020, 98, 241–265. [Google Scholar] [CrossRef] [Green Version]
  6. Falconi, G.P.; Holzapfel, F. Position tracking of a hexacopter using a geometric backstepping control law—Experimental results. In Proceedings of the 2014 IEEE International Conference on Aerospace Electronics and Remote Sensing Technology, Yogyakarta, Indonesia, 13–14 November 2014; pp. 20–25. [Google Scholar] [CrossRef]
  7. Miao, C.X.; Fang, J.C. An Adaptive Three-Dimensional Nonlinear Path following Method for a Fix-Wing Micro Aerial Vehicle. Int. J. Adv. Robot. Syst. 2012, 9, 206. [Google Scholar] [CrossRef]
  8. Schatz, S.P.; Holzapfel, F. Modular trajectory/path following controller using nonlinear error dynamics. In Proceedings of the ICARES 2014, Yogyakarta, Indonesia, 13–14 November 2014; pp. 157–163. [Google Scholar] [CrossRef]
  9. Schatz, S.P.; Schneider, V.; Karlsson, E.; Holzapfel, F.; Baier, T.; Dorhofer, C.; Hochstrasser, M.; Gabrys, A.; Krause, C.; Lauffs, P.J.; et al. Flightplan flight tests of an experimental DA42 general aviation aircraft. In Proceedings of the 2016 14th International Conference on Control, Automation, Robotics & Vision, Phuket, Thailand, 13–15 November 2016; pp. 1–6. [Google Scholar] [CrossRef]
  10. Schatz, S. Development and Flight-Testing of a Trajectory Controller Employing Full Nonlinear Kinematics. Ph.D. Thesis, Technische Universität München, München, Germany, 2018. [Google Scholar]
  11. Bhardwaj, P.; Raab, S.A.; Zhang, J.; Holzapfel, F. Thrust command based Integrated Reference Model with Envelope Protections for Tilt-rotor VTOL Transition UAV. In Proceedings of the AIAA Aviation 2019 Forum, Dallas, TX, USA, 17–21 June 2019; American Institute of Aeronautics and Astronautics: Reston, Virginia, 2019. [Google Scholar] [CrossRef]
  12. Raab, S.A.; Zhang, J.; Bhardwaj, P.; Holzapfel, F. Consideration of Control Effector Dynamics and Saturations in an Extended INDI Approach. In Proceedings of the AIAA Aviation 2019 Forum, Dallas, TX, USA, 17–21 June 2019; American Institute of Aeronautics and Astronautics: Reston, Virginia, 2019. [Google Scholar] [CrossRef]
  13. Grant, P.R.; Haycock, B. Effect of Jerk and Acceleration on the Perception of Motion Strength. J. Aircr. 2008, 45, 1190–1197. [Google Scholar] [CrossRef]
  14. Kyriakopoulos, K.J.; Saridis, G.N. (Eds.) Minimum jerk path generation. In Proceedings of the 1988 IEEE International Conference on Robotics and Automation, Philadelphia, PA, USA, 24–29 April 1988. [Google Scholar] [CrossRef]
  15. Schneider, V.; Piprek, P.; Schatz, S.P.; Baier, T.; Dorhofer, C.; Hochstrasser, M.; Gabrys, A.; Karlsson, E.; Krause, C.; Lauffs, P.J.; et al. Online trajectory generation using clothoid segments. In Proceedings of the 2016 14th International Conference on Control, Automation, Robotics & Vision, Phuket, Thailand, 13–15 November 2016; pp. 1–6. [Google Scholar] [CrossRef]
  16. Schneider, V. Trajectory Generation for Integrated Flight Guidance. Ph.D. Thesis, Technische Universität München, München, Germeny, 2018. [Google Scholar]
  17. Diaz-Mercado, Y.; Lee, S.G.; Egerstedt, M.; Young, S.Y. Optimal trajectory generation for next generation flight management systems. In Proceedings of the IEEE/AIAA 32nd Digital Avionics Systems Conference (DASC), East Syracuse, NY, USA, 5–10 October 2013; pp. 3–10. [Google Scholar] [CrossRef]
  18. Zheng, H.; Piprek, P.; Hong, H.; Holzapfel, F.; Tang, S. Smooth Sub-Optimal Trajectory Generation for Transition Maneuvers. IEEE Access 2020, 8, 61035–61042. [Google Scholar] [CrossRef]
  19. Piprek, P.; Gros, S.; Holzapfel, F. Rare Event Chance-Constrained Optimal Control Using Polynomial Chaos and Subset Simulation. Processes 2019, 7, 185. [Google Scholar] [CrossRef] [Green Version]
  20. Piprek, P.; Gros, S.; Holzapfel, F. A Distributed Robust Optimal Control Framework Based on Polynomial Chaos. In Proceedings of the 5th CEAS Specialist Conference on Guidance, Navigation & Control, Milano, Italy, 3–5 April 2019. [Google Scholar]
  21. Dicheva, S.; Bestaoui, Y. Three-Dimensional A* Dynamic Mission Planning for an Airborne Launch Vehicle. J. Aerosp. Inf. Syst. 2014, 11, 98–106. [Google Scholar] [CrossRef]
  22. Piprek, P.; Schneider, V.; Fafard, V.; Schatz, S.P.; Dörhöfer, C.; Lauffs, P.J.; Peter, L.; Holzapfel, F. Enhanced Kinematics Calculation for an Online Trajectory Generation Module. Transp. Res. Procedia 2018, 29, 312–322. [Google Scholar] [CrossRef]
Figure 1. Visualization of basic idea for nonlinear trajectory/path-following controller based on trajectory reference point F , i.e., the foot point along the desired trajectory, as well as the definition of the path deviations and their derivatives (adapted and extended from the work in [8]).
Figure 1. Visualization of basic idea for nonlinear trajectory/path-following controller based on trajectory reference point F , i.e., the foot point along the desired trajectory, as well as the definition of the path deviations and their derivatives (adapted and extended from the work in [8]).
Applsci 10 08760 g001
Figure 2. Maximal absolute lateral deviation over turn radius for 180 -turn maneuver.
Figure 2. Maximal absolute lateral deviation over turn radius for 180 -turn maneuver.
Applsci 10 08760 g002
Figure 3. Maximal absolute lateral specific force over turn radius for 180 turn maneuver.
Figure 3. Maximal absolute lateral specific force over turn radius for 180 turn maneuver.
Applsci 10 08760 g003
Figure 4. Commanded vs. actually flown trajectory for four-dimensional trajectory control case.
Figure 4. Commanded vs. actually flown trajectory for four-dimensional trajectory control case.
Applsci 10 08760 g004
Figure 5. Trajectory position, velocity, and acceleration deviations in all three axes along the trajectory.
Figure 5. Trajectory position, velocity, and acceleration deviations in all three axes along the trajectory.
Applsci 10 08760 g005
Figure 6. Specific force and its derivative along the trajectory.
Figure 6. Specific force and its derivative along the trajectory.
Applsci 10 08760 g006
Table 1. Constants used in the simulation model.
Table 1. Constants used in the simulation model.
NameSymbolValueUnit
Actuator natural frequency ω 0 20 rad s
Actuator damping ζ 1
Gravitational accelerationg 9.81 m s 2
Linear force derivative Z α 3 1 s
Position error gain matrix K p I 3 × 3
Velocity error gain matrix K d 3 · I 3 × 3
Acceleration error gain matrix K dd 9 · I 3 × 3
Aircraft massm100 kg
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Piprek, P.; Marb, M.M.; Bhardwaj, P.; Holzapfel, F. Trajectory/Path-Following Controller Based on Nonlinear Jerk-Level Error Dynamics. Appl. Sci. 2020, 10, 8760. https://doi.org/10.3390/app10238760

AMA Style

Piprek P, Marb MM, Bhardwaj P, Holzapfel F. Trajectory/Path-Following Controller Based on Nonlinear Jerk-Level Error Dynamics. Applied Sciences. 2020; 10(23):8760. https://doi.org/10.3390/app10238760

Chicago/Turabian Style

Piprek, Patrick, Michael M. Marb, Pranav Bhardwaj, and Florian Holzapfel. 2020. "Trajectory/Path-Following Controller Based on Nonlinear Jerk-Level Error Dynamics" Applied Sciences 10, no. 23: 8760. https://doi.org/10.3390/app10238760

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop