arxiv: v1 [math.oc] 21 Sep 2011
|
|
- Marian Hart
- 6 years ago
- Views:
Transcription
1 Nonlinear Robust Tracking Control of a Quadrotor UAV on SE(3) Taeyoung Lee Melvin Leok and N. Harris McClamroch arxiv:9.4457v [math.oc 2 Sep 2 Abstract This paper provides nonlinear tracking control systems for a quadrotor unmanned aerial vehicle (UAV) that are robust to bounded uncertainties. A mathematical model of a quadrotor UAV is defined on the special Euclidean group and nonlinear output-tracking controllers are developed to follow () an attitude command and (2) a position command for the vehicle center of mass. The controlled system has the desirable properties that the tracking errors are uniformly ultimately bounded and the size of the ultimate bound can be arbitrarily reduced by control system parameters. Numerical eamples illustrating comple maneuvers are provided. I. INTRODUCTION A quadrotor unmanned aerial vehicle (UAV) consists of two pairs of counter-rotating rotors and propellers. It has been envisaged for various applications such as surveillance or mobile sensor networks as well as for educational purposes and several control systems have been studied. Linear control systems have been widely used to enhance the stability properties of an equilibrium of a quadrotor UAV [ [2 [3. In [4 the quadrotor dynamics is modeled as a collection of simplified hybrid dynamic modes and reachability sets are analyzed to guarantees the safety and performance for larger area of operating conditions. Several nonlinear controllers have been developed as well. Backstepping and sliding mode techniques are applied in [5 [6 and a nonlinear H controller is studied in [7. An adaptive neural network based control system is developed in [8. Since all of these controllers are based on Euler angles they ehibit singularities when representing comple rotational maneuvers of a quadrotor UAV thereby significantly restricting their ability to achieve comple flight maneuvers. An attitude control system based on quaternions is applied to a quadrotor UAV [9. Quaternions do not have singularities but they have ambiguities in representing an attitude as the three-sphere S 3 double-covers SO(3). As a result in a quaternion-based attitude control system convergence to a single attitude implies convergence to either of the two disconnected antipodal points on S 3 [. Therefore depending on the particular choice of control inputs a quaternion-based control system may become discontinuous Taeyoung Lee Mechanical and Aerospace Engineering The George Washington University Washington DC 252 tylee@gwu.edu Melvin Leok Mathematics University of California at San Diego La Jolla CA 9293 mleok@math.ucsd.edu N. Harris McClamroch Aerospace Engineering University of Michigan Ann Arbor MI 489 nhm@umich.edu This research has been supported in part by NSF under grants CMMI This research has been supported in part by NSF under grants DMS DMS-52 DMS-687 and CMMI when applied to actual attitude dynamics [ and it may also ehibit unwinding behavior where the controller rotates a rigid body through unnecessarily large angles [2 [3. Attitude control systems also have been developed directly on the special orthogonal group SO(3) to avoid the singularities associated with Euler-angles and the ambiguity of quaternions [4 [5 [6 [7. By following this geometric approach the dynamics of a quadrotor UAV is globally epressed on the special Euclidean group SE(3) and nonlinear control systems are developed to track outputs of several flight modes namely an attitude controlled flight mode a position controlled flight mode and a velocity controlled flight mode [8. Several aggressive maneuvers of a quadrotor UAV are also demonstrated based on a hybrid control architecture. This is particularly desirable since complicated reachability set analysis is not required to guarantee a safe switching between different flight modes as the region of attraction for each flight mode covers the configuration space almost globally. In this paper we etend the results of [8 to construct nonlinear robust tracking control systems on SE(3) for a quadrotor UAV. We assume that there eist unstructured bounded uncertainties with pre-determined bounds on the translational dynamics and the rotation dynamics of a quadrotor UAV. Output tracking control systems are developed to follow an attitude command or a position command for the vehicle center of mass. We show that the tracking errors are uniformly ultimately bounded and the size of the ultimate bound can be arbitrarily reduced. The robustness of the proposed tracking control systems are critical in generating comple maneuvers as the impact of the several aerodynamic effects resulting from the variation in air speed is significant even at moderate velocities [2. The paper is organized as follows. We develop a globally defined model for the translational and rotational dynamics of a quadrotor UAV in Section II. A hybrid control architecture is introduced and a robust attitude tracking control system is developed in Section III. Section IV present results for a robust position tracking followed by numerical eamples in Section V. II. QUADROTOR DYNAMICS MODEL Consider a quadrotor UAV model illustrated in Figure. This is a system of four identical rotors and propellers located at the vertices of a square which generate a thrust and torque normal to the plane of this square. We choose an inertial reference frame e e 2 e 3 and a body-fied frame b b 2 b 3. The origin of the body-fied frame is located at the center of mass of this vehicle. The first and
2 e 3 e e 2 Fig.. f 3 f 2 f 4 b2 Quadrotor model the second aes of the body-fied frame b b 2 lie in the plane defined by the centers of the four rotors as illustrated in Figure. The third body-fied ais b 3 is normal to this plane. Each of the inertial reference frame and the bodyfied reference frame consist of a triad of orthogonal vectors defined according to the right hand rule. Define m R J R 3 3 R SO(3) Ω R 3 R 3 v R 3 d R f i R τ i R the total mass the inertia matri with respect to the body-fied frame the rotation matri from the body-fied frame to the inertial frame the angular velocity in the body-fied frame the position vector of the center of mass in the inertial frame the velocity vector of the center of mass in the inertial frame the distance from the center of mass to the center of each rotor in the b b 2 plane the thrust generated by the i-th propeller along the b 3 ais the torque generated by the i-th propeller about the b 3 ais f R the total thrust magnitude i.e. f = 4 i= f i M R 3 f b3 R b the total moment vector in the bodyfied frame The configuration of this quadrotor UAV is defined by the location of the center of mass and the attitude with respect to the inertial frame. Therefore the configuration manifold is the special Euclidean group SE(3) which is the semidirect product of R 3 and the special orthogonal group SO(3) = R R 3 3 R T R = I det R =. The following conventions are assumed for the rotors and propellers and the thrust and moment that they eert on the quadrotor UAV. We assume that the thrust of each propeller is directly controlled and the direction of the thrust of each propeller is normal to the quadrotor plane. The first and third propellers are assumed to generate a thrust along the direction of b 3 when rotating clockwise; the second and fourth propellers are assumed to generate a thrust along the same direction of b 3 when rotating counterclockwise. Thus the thrust magnitude is f = 4 i= f i and it is positive when the total thrust vector acts along b 3 and it is negative when the total thrust vector acts along b 3. By the definition of the rotation matri R SO(3) the total thrust vector is given by fre 3 R 3 in the inertial frame. We also assume that the torque generated by each propeller is directly proportional to its thrust. Since it is assumed that the first and the third propellers rotate clockwise and the second and the fourth propellers rotate counterclockwise to generate a positive thrust along the direction of b 3 the torque generated by the i-th propeller about b 3 can be written as τ i = ( ) i c τf f i for a fied constant c τf. All of these assumptions are common [3 [9. Under these assumptions the moment vector in the bodyfied frame is given by f f M M 2 = d d f 2 d d f 3. () M 3 c τf c τf c τf c τf f 4 The determinant of the above 4 4 matri is 8c τf d 2 so it is invertible when d and c τf. Therefore for given thrust magnitude f and given moment vector M the thrust of each propeller f f 2 f 3 f 4 can be obtained from (). Using this equation the thrust magnitude f R and the moment vector M R 3 are viewed as control inputs in this paper. The equations of motion of the quadrotor UAV can be written as ẋ = v (2) m v = mge 3 fre 3 + (3) Ṙ = RˆΩ (4) J Ω + Ω JΩ = M + R (5) where the hat map ˆ : R 3 so(3) is defined by the condition that ˆy = y for all y R 3 (see Appendi A). The inverse of the hat map is denoted by the vee map : so(3) R 3. Unstructured uncertainties in the translational dynamics and the rotational dynamics of a quadrotor UAV are denoted by and R R 3 respectively. We assume that uncertainties are bounded: δ R δ R (6) for known positive constants δ δ R R. Throughout this paper λ m ( ) and λ M ( ) denote the minimum eignevalue and the maimum eigenvalue of a matri respectively. III. ATTITUDE CONTROLLED FLIGHT MODE A. Flight Modes Since the quadrotor UAV has four inputs it is possible to achieve asymptotic output tracking for at most four quadrotor UAV outputs. The quadrotor UAV has three translational and three rotational degrees of freedom; it is not possible to achieve asymptotic output tracking of both attitude and position of the quadrotor UAV. This motivates us to introduce
3 several flight modes namely () an attitude controlled flight mode and (2) a position controlled flight mode. A comple flight maneuver can be defined by specifying a concatenation of flight modes together with conditions for switching between them; for each flight mode one also specifies the desired or commanded outputs as functions of time. Unlike a hybrid flight control system that requires reachability analyses [4 the proposed control system is robust to switching conditions since each flight mode has almost global stability properties and it is straightforward to design a comple maneuver of a quadrotor UAV. In this section an attitude controlled flight mode is considered where the outputs are the attitude of the quadrotor UAV and the controller for this flight mode achieves asymptotic attitude tracking. B. Attitude Tracking Errors Suppose that an arbitrary smooth attitude command R d (t) SO(3) is given. The corresponding angular velocity command is obtained by the attitude kinematics equation ˆΩ d = Rd T Ṙd. We first define errors associated with the attitude dynamics of the quadrotor UAV. The attitude error function studied in [4 [9 [2 and several properties are summarized as follows. Proposition : For a given tracking command (R d Ω d ) and the current attitude and angular velocity (R Ω) we define an attitude error function Ψ : SO(3) SO(3) R an attitude error vector e R R 3 and an angular velocity error vector e Ω R 3 as follows: Ψ(R R d ) = 2 tr[ I R T d R (7) e R = 2 (RT d R R T R d ) (8) Then the following statements hold: e Ω = Ω R T R d Ω d (9) (i) Ψ is locally positive-definite about R = R d. (ii) the left-trivialized derivative of Ψ is given by T IŁ R (D R Ψ(R R d )) = e R. () (iii) the critical points of Ψ where e R = are R d R d ep(πŝ) s S 2. (iv) a lower bound of Ψ is given as follows: 2 e R(R R d ) 2 Ψ(R R d ) () (v) Let ψ be a positive constant that is strictly less than 2. If Ψ(R R d ) < ψ < 2 then an upper bound of Ψ is given by Ψ(R R d ) 2 ψ e R(R R d ) 2. (2) Proof: See [2. C. Attitude Tracking Controller We now introduce a nonlinear controller for the attitude controlled flight mode described by an epression for the moment vector: M = k R e R k Ω e Ω + Ω JΩ J(ˆΩR T R d Ω d R T R d Ωd ) + µ R (3) δr 2 µ R = A δ R e A + ɛ R (4) e A = e Ω + c 2 J e R (5) where k R k Ω c 2 ɛ R are positive constants. In this attitude controlled mode it is possible to ignore the translational motion of the quadrotor UAV; consequently the reduced model for the attitude dynamics are given by equations (4) (5) using the controller epression (3)-(5). We now state the result that the tracking errors (e R e Ω ) are uniformly ultimately bounded. Proposition 2: (Robustness of Attitude Controlled Flight Mode) Suppose that the initial attitude error satisfies Ψ(R() R d ()) < ψ 2 < 2 (6) for a constant ψ 2. Consider the control moment M defined in (3)-(5). For positive constants k R k Ω the constants c 2 ɛ R are chosen such that c 2 < min k Ω 4k Ω k R λ m (J) 2 k 2 Ω λ M (J) + 4k R λ m (J) 2 k R λ m (J) (7) ɛ R < λ m(m 2 )λ m (W 2 ) ψ 2 (2 ψ 2 ) (8) λ M (M 22 ) where the matrices M 2 M 22 W 2 R 2 2 are given by M 2 = [ kr c 2 M 2 c 2 λ m (J) 22 = [ 2kR 2 ψ 2 c 2 2 c 2 λ M (J) W 2 = [ c2k R λ M (J) k Ω c 2 Then the attitude tracking errors (e R e Ω ) are uniformly ultimately bounded and the ultimate bound is given by e R 2 + e Ω 2 λ M (M 22 ) λ m (M 2 )λ m (W 2 ) ɛ R. (9) Proof: See Appendi B. From (6) the initial attitude error should be less than 8 in terms of the rotation angle about the eigenais between R and R d. We can further show that the attitude tracking errors eponentially converges to (9) where the size of the ultimate bound can be reduced by the controller parameter ɛ R. It is also possible to achieve eponential attractiveness if the constant ɛ R in (4) is replaced by ɛ R ep( βt) for β >. All of these results can be applied to a nonlinear robust control problem for the attitude dynamics of any rigid body. Asymptotic tracking of the quadrotor attitude does not require specification of the thrust magnitude. As an auiliary problem the thrust magnitude can be chosen in many.
4 different ways to achieve an additional translational motion objective. For eample it can be used to asymptotically track a quadrotor altitude command [8. Since the translational motion of the quadrotor UAV can only be partially controlled; this flight mode is most suitable for short time periods where an attitude maneuver is to be completed. IV. POSITION CONTROLLED FLIGHT MODE We now introduce a nonlinear controller for the position controlled flight mode. This flight mode requires analysis of the coupled translational and rotational equations of motion; hence we make use of the notation and analysis in the prior section to describe the properties of the closed loop system in this flight mode. A. Position Tracking Errors An arbitrary smooth position tracking command d (t) R 3 is chosen. The position tracking errors for the position and the velocity are given by: e = d (2) e v = v ẋ d. (2) Following the prior definition of the attitude error and the angular velocity error we define e R = 2 (RT c R R T R c ) e Ω = Ω R T R c Ω c (22) and the computed attitude R c (t) SO(3) and computed angular velocity Ω c R 3 are given by R c = [b c ; b 3c b c ; b 3c ˆΩc = R T c Ṙc (23) where b 3c S 2 is defined by b 3c = k e k v e v mge 3 + mẍ d + µ k e k v e v mge 3 + mẍ d + µ (24) and b c S 2 is selected to be orthogonal to b 3c thereby guaranteeing that R c SO(3). The constants k k v are positive and the control input term µ is defined later in (29). We assume that k e k v e v mge 3 + mẍ d + µ (25) and the commanded acceleration is uniformly bounded such that for a given positive constant B. B. Position Tracking Controller mge 3 + mẍ d < B (26) The nonlinear controller for the position controlled flight mode described by control epressions for the thrust magnitude and the moment vector are: f = (k e + k v e v + mge 3 mẍ d µ ) Re 3 (27) M = k R e R k Ω e Ω + Ω JΩ J(ˆΩR T R c Ω c R T R c Ωc + µ R ) (28) d (b d ) Fig. 2. Force controller b 3c Controller Moment controller f M v R Ω Quadrotor Dynamics Controller structure for position controlled flight mode δ τ+2 e B e B τ µ = δ τ+ e B τ+ + ɛ τ+ (29) e B = e v + c m e (3) δr 2 µ R = A δ R e A + ɛ R (3) e A = e Ω + c 2 J e R (32) where k k v k R k Ω c c 2 ɛ ɛ R τ are positive constants and τ > 2. The nonlinear controller given by equations (27) (28) can be given a backstepping interpretation. The computed attitude R c given in equation (23) is selected so that the thrust ais b 3 of the quadrotor UAV tracks the computed direction given by b 3c in (24) which is a direction of the thrust vector that achieves position tracking. The moment epression (28) causes the attitude of the quadrotor UAV to asymptotically track R c and the thrust magnitude epression (27) achieves asymptotic position tracking. The closed loop system for this position controlled flight mode is illustrated in Figure 2. The corresponding closed loop control system is described by equations (2)-(5) using the controller epressions (27)-(32). We now state the result that the tracking errors (e e v e R e Ω ) are uniformly ultimately bounded. Proposition 3: (Robustness of Position Controlled Flight Mode) Suppose that the initial conditions satisfy Ψ(R() R c ()) < ψ < (33) e () < e ma (34) for positive constants ψ e ma. Consider the control inputs f M defined in (27)-(32). For positive constants k k v we choose positive constants c c 2 k R k Ω ɛ ɛ R such that c < min k v ( α) c 2 < min k Ω ɛ + ɛ R < 4mk k v ( α) 2 k 2 v( + α) 2 + 4mk ( α) k m 4k Ω k R λ m (J) 2 k 2 Ω λ M (J) + 4k R λ m (J) 2 k R λ m (J) (35) (36) λ m (W 2 ) > W 2 2 4λ m (W ) (37) minλ m (M ) λ m (M 2 ) mine 2 ma ψ (2 ψ ) maλ M (M 2 ) λ M (M 22 ) λ m (W ) (38)
5 where α = ψ (2 ψ ) and the matrices M M 2 M 2 M 22 W W 2 W 2 W R 2 2 are given by M = [ k c M 2 c m 2 = [ k c 2 c m M 2 = [ kr c 2 M 2 c 2 λ m (J) 22 = [ 2kR 2 ψ c 2 2 c 2 λ M (J) [ ck ckv W = m ( α) 2m ( + α) ckv 2m ( + α) k v( α) c [ c m (B + δ W 2 = ) B + δ + k e ma W 2 = W = [ c2k R λ M (J) k Ω c 2 [ λm (W ) 2 W W 2 2 λ m (W 2 ) Then the tracking errors (e e v e R e Ω ) are uniformly ultimately bounded and the ultimate bound is given by e 2 + e v 2 + e R 2 + e Ω 2 < maλ M (M 2 ) λ M (M 22) minλ m (M ) λ m (M 2 )λ m (W ) (ɛ + ɛ R ). (39) Proof: See Appendi C. This proposition shows that the proposed control system is robust to bounded and unstructured uncertainties in the dynamics of a quadrotor UAV. Similar to Proposition 2 the ultimate bound can be arbitrarily reduced by choosing smaller ɛ ɛ R and it is possible to obtain eponential attractiveness. Proposition 3 requires that the initial attitude error is less than 9 in (33). Suppose that this is not satisfied i.e. Ψ(R() R c ()) < 2. We can still apply Proposition 2 which states that the attitude error eponentially decreases until it enters the ultimate bound given by (9). If the constant ɛ R is sufficiently small we can guarantee that the attitude error function decreases to satisfy (33) in a finite time. Therefore by combining the results of Proposition 2 and 3 we can show ultimate boundedness of the tracking errors when Ψ(R() R c ()) < 2. Proposition 4: (Robustness of Position Controlled Flight Mode with a Larger Initial Attitude Error) Suppose that the initial conditions satisfy. Ψ(R() R c ()) < ψ 2 < 2 (4) e () < e ma (4) for a constant ψ 2 e ma. Consider the control inputs f M defined in (27)-(32) where the control parameters k k v k R k Ω c c 2 ɛ ɛ R satisfy (35)-(38) for a positive constant ψ <. If the constant ɛ R is sufficiently small such that ɛ R < λ m(m 2 )λ m (W 2 ) ψ (2 ψ ) (42) λ M (M 22 ) Plane normal to b 3c b3c bd b2c = bc = Proj[ b d = b 2c b 3c b 3c b d b 3c b d Fig. 3. Convergence property of the first body-fied ais: in the proposed control system b 3c is determined by (24). We choose a desired direction of the first body fied ais namely b d that is not parallel to b 3c and project it on to the plane normal to b 3c to obtain b c. This guarantees that the first body-fied ais converges to b c and therefore it asymptotically lies in the plane spanned by b d and b 3c. As b 3c converges to the direction of ge 3 ẍ d in (24) this allows us to specify the direction of the first body-fied ais in the plane normal to ge 3 ẍ d. For all cases the ultimate convergence error is described by (39). then the tracking errors (e e v e R e Ω ) are uniformly ultimately bounded. Proof: See Appendi D. C. Direction of the First Body-Fied Ais As described above the construction of the orthogonal matri R c involves having its third column b 3c specified by a normalized feedback function and its first column b c is chosen to be orthogonal to the third column. The unit vector b c can be arbitrarily chosen in the plane normal to b 3c which corresponds to a one-dimensional degree of choice. This reflects the fact that the quadrotor UAV has four control inputs that are used to track a three-dimensional position command. By choosing b c properly we constrain the asymptotic direction of the first body-fied ais. Here we propose to specify the projection of the first body-fied ais onto the plane normal to b 3c. In particular we choose a desired direction b d S 2 that is not parallel to b 3c and b c is selected as b c = Proj[b d where Proj[ denotes the normalized projection onto the plane perpendicular to b 3c. In this case the first body-fied ais does not converge to b d but it converges to the projection of b d i.e. b b c = Proj[b d as t up to the ultimate bound described by (39). In other words the first body-fied ais converges to a small neighborhood of the intersection of the plane normal to b 3c and the plane spanned by b 3c and b d. This can be used to specify the heading direction of a quadrotor UAV in the horizontal plane (see Figure 3 and [8 for details). V. NUMERICAL EXAMPLES Numerical results are presented to demonstrate the prior approach for performing comple flight maneuvers. The parameters are chosen to match a quadrotor UAV described in [2. J = [ kg m 2 d =.35 m c τf = m. The controller parameters are chosen as follows: m = 4.34 kg k = 59.2 k v = 24.3 k R = 8.8 k Ω =.54 c = 3.6 c 2 =.6 ɛ = ɛ R =.4.
6 (a) Attitude error function Ψ (b) Position error e (m) (a) Attitude error function Ψ (b) Position error e (m) (c) Angular velocity error e Ω (rad/sec) 5 (d) Thrust of each rotor (N) (c) Angular velocity error e Ω (rad/sec) (d) Thrust of each rotor (N) Fig. 4. heli Case I: robust position controlled flight mode to follow an elliptic Fig. 5. Case I: position controlled flight mode to follow an elliptic heli. The robust control input terms are set to zero i.e. µ = µ R = for comparison with Figure 4 We consider a fied disturbance for the translational dynamics and an oscillatory disturbance for the rotational dynamics as follows: = [ T N R (t) = 2 3 [sin(8πt) sin(πt) cos(4πt) T Nm. The corresponding bounds of the disturbances are given by δ = 4.34 and δ R = 2. We consider the following two cases. Case I (elliptic heli): The initial conditions are given by () = [. T m R() = I v() = [ T m/s Ω() = [ T rad/s. The desired position command is an elliptic heli given by d (t) = [.4t.4 sin(πt).6 cos(πt) T m and the desired heading direction is fied as b d = [ T. This corresponds to the position controlled flight mode described in Proposition 3 as the initial attitude error is Ψ() =.4 <. Figure 4 shows simulation results where the position tracking error converges to a small neighborhood of the zero tracking errors and the terminal tracking error is.2 cm. For comparison we set the robust control input terms to zero i.e. µ = µ R = and we repeat numerical simulations to obtain Figure 5. It is observed that the angular velocity tracking error is mostly driven by the disturbance R and the corresponding position tracking error is larger than. m. This illustrates the robustness of the proposed control system for a comple maneuver with larger disturbances. Case II (hovering): The initial conditions are given by () = [. T m R() = ep(.99πê ) v() = [ T m/s Ω() = [ T rad/s where e = [ R 3. The desired position command is given by d (t) = [ T m and the desired heading direction is fied as b d = [ T. This describes a case that a quadrotor UAV should recover from an initially upside-down configuration. The initial attitude error is given by (Ψ() =.9995) < 2 and therefore it corresponds to Proposition 4 that is based on both of the attitude controlled flight mode and the position controlled flight mode. Figure 6 illustrates ecellent convergence properties of the proposed control system for a large initial attitude error where the terminal position tracking error is.2 cm. Figure Figure 7 shows relatively poor tracking performances with a slower convergence when there are no robust control input terms proposed in this paper. A. Properties of the Hat Map APPENDIX The hat map ˆ : R 3 so(3) is defined as ˆ = (43) 2 for = [ ; 2 ; 3 R 3. This identifies the Lie algebra so(3) with R 3 using the vector cross product in R 3. The inverse of the hat map is referred to as the vee map
7 (a) Attitude error function Ψ (b) Position error e (m) (a) Attitude error function Ψ (b) Position error e (m) 5 (c) Angular velocity error e Ω (rad/sec) 5 (d) Thrust of each rotor (N) (c) Angular velocity error e Ω (rad/sec) (d) Thrust of each rotor (N) Fig. 6. Case II: robust position controlled flight mode to recover from an initially upside-down configuration Fig. 7. Case II: robust position controlled flight mode to recover from an initially upside-down configuration. The robust control input terms are set to zero i.e. µ = µ R = for comparison with Figure 6 : so(3) R 3. Several properties of the hat map are summarized as follows. ˆy = y = y = ŷ (44) 2 tr[ˆŷ = T y (45) tr[ˆa = tr[aˆ = 2 tr[ˆ(a A T ) = T (A A T ) (46) ˆA + A T ˆ = (tr[a I 3 3 A ) (47) RˆR T = (R) (48) for any y R 3 A R 3 3 and R SO(3). B. Proof of Proposition 2 We first find the error dynamics for e R e Ω and define a Lyapunov function. Then we find conditions on control parameters to guarantee the boundedness of tracking errors. a) Attitude Error Dynamics: The attitude error dynamics for Ψ e R e Ω are developed in [2 and they are summarized as follows: d dt (Ψ(R R d)) = e R e Ω (49) ė R = E(R R d )e Ω (5) ė Ω = J ( Ω JΩ + u + R ) + ˆΩR T R d Ω d R T R d Ωd (5) where the matri E(R R d ) R 3 3 is given by E(R R d ) = 2 (tr[rt R d I R T R d ). (52) We can show that E(R R d ) to obtain ė R e Ω. (53) Substituting the control moment (3) into (5) Jė Ω = k R e R k Ω e Ω + R + µ R. (54) In short the attitude error dynamics are given by equations (49) (5) (54) and they satisfy (53). b) Lyapunov Candidate: Let a Lyapunov candidate V 2 be V 2 = 2 e Ω Je Ω + k R Ψ(R R d ) + c 2 e R e Ω. (55) We analyzes the properties of V 2 along the solutions of the controlled system in the following domain D 2 : D 2 = (R Ω) SO(3) R 3 Ψ(R R d ) < ψ 2. (56) From () (2) the attitude error function is bounded in D 2 as follows: 2 e R 2 Ψ(R R d ) e R 2 (57) 2 ψ 2 which implies that Ψ is positive-definite and decrescent. It follows that the Lyapunov function V 2 is bounded as z T 2 M 2 z 2 V 2 z T 2 M 22 z 2 (58) where z 2 = [ e R e Ω T R 2 and the matrices M 2 M 22 are given by M 2 = [ kr c 2 M 2 c 2 λ m (J) 22 = [ 2kR 2 ψ 2 c 2. 2 c 2 λ M (J) (59) From equations (49) (5) (54) the time derivative of V 2 along the solution of the controlled system is given by V 2 = k Ω e Ω 2 c 2 k R e R J e R + c 2 E(R R d )e Ω e Ω
8 c 2 k Ω e R J e Ω + (e Ω + c 2 J e R ) ( R + µ R ). (6) Since E(R d R) this is bounded by V 2 z T 2 W 2 z 2 + e A ( R + µ R ) (6) where e A = e Ω +c 2 J e R R 3 and the matri W 2 R 2 2 is given by W 2 = [ c2k R λ M (J) k Ω c 2. (62) Substituting (4) the last term of (6) is bounded by to obtain e A ( R + µ R ) = δ R e A δ2 R e A 2 δ R e A + ɛ R δ R e A = ɛ R ɛ R δ R e A + ɛ R V 2 z T 2 W 2 z 2 + ɛ R (63) c) Boundedness: The condition (7) for the constant c 2 guarantees that the matri W 2 in (63) and the matrices M 2 M 22 in (58) are positive-definite. Therefore we obtain λ m (M 2 ) z 2 2 V 2 λ M (M 22 ) z 2 2 (64) V 2 λ m (W 2 ) z ɛ R. (65) This implies that V2 < when λm (M 22 ) V 2 > λ m (W 2 ) ɛ R d. Consider a sub-level set of the Lyapunov function V 2 defined as S γ = (R Ω) SO(3) R 3 V 2 γ for a positive constant γ. If γ satisfies the following inequality γ < λ m (M 2 )ψ 2 (2 ψ 2 ) d 2 then we can guarantee that S γ is a subset of the domain D 2 defined in (56). In short a sub-level set of the Lyapunov function S γ is a positively invariant set when d < γ < d 2 and any solution starting in S γ eponentially converges to S d. To guarantee the eistence of such a set we require d = λ M (M 22 ) λ m (W 2 ) ɛ R < λ m (M 2 )ψ 2 (2 ψ 2 ) d 2 which can be achieved by (8). Then according to Theorem 5. in [22 the attitude tracking errors are uniformly ultimately bounded and the corresponding ultimate bound is estimated by S d z 2 2 λ M (M 22 ) λ m (M 2 )λ m (W 2 ) ɛ R. C. Proof of Proposition 3 We first derive the tracking error dynamics and a Lyapunov function for the translational dynamics of a quadrotor UAV and later it is combined with the stability analyses of the rotational dynamics in Appendi B to guarantee the boundedness of tracking errors. The subsequent analyses are developed in the domain D D = (e e v R e Ω ) R 3 R 3 SO(3) R 3 e < e ma Ψ < ψ (66) Similar to (57) we can show that 2 e R 2 Ψ(R R c ) e R 2. (67) 2 ψ a) Translational Error Dynamics: The time derivative of the position error is ė = e v. The time-derivative of the velocity error is given by mė v = mẍ mẍ d = mge 3 fre 3 mẍ d +. (68) Consider the quantity e T 3 Rc T Re 3 which represents the cosine of the angle between b 3 = Re 3 and b c3 = R c e 3. Since Ψ(R R c ) represents the cosine of the eigen-ais rotation angle between R c and R we have > e T 3 Rc T Re 3 > Ψ(R R c ) > in D. Therefore the quantity e is T 3 RT c well-defined. To rewrite the error dynamics of e Re3 v in terms f of the attitude error e R we add and subtract e R ce T 3 3 RT c to the right hand side of (68) to obtain Re3 mė v = mge 3 mẍ d where X R 3 is defined by X = f e T 3 RT c Re 3 R c e 3 X + (69) f e T 3 RT c Re 3 ((e T 3 R T c Re 3 )Re 3 R c e 3 ). (7) Let A = k e k v e v mge 3 +mẍ d +µ. Then from (27) (24) we have f = A Re 3 and b 3c = R c e 3 = A/ A i.e. A = A R c e 3. By combining these we obtain f = ( A R c e 3 ) Re 3. Therefore the third term of the right hand side of (69) can be written as f e T R c e 3 = ( A R ce 3 ) Re 3 3 RT c Re 3 e T A 3 RT c Re 3 A = A = k e k v e v mge 3 + mẍ d + µ. Substituting this into (69) the error dynamics of e v can be written as mė v = k e k v e v X + + µ. (7) b) Lyapunov Candidate for Translation Dynamics: Let a Lyapunov candidate V be V = 2 k e m e v 2 + c e e v. (72) The derivative of V along the solution of (7) is given by V = (k v c ) e v 2 c k m e 2 c k v m e e v c + X + + µ m e + e v. (73)
9 From (29) (3) the last part of (73) is bounded by δ τ+2 e B τ+2 e B ( + µ ) δ e B δ τ+ e B τ+ + ɛ τ+ δ e B ɛ τ+ = δ τ+ e B τ+ + ɛ τ+ ɛ. (74) The last inequality is satisfied since if δ e B ɛ δ e B and if δ e B > ɛ δ τ+ e B τ+ δ τ+ e B τ+ + ɛ τ+ ɛ τ+ δ τ+ e B τ+ + ɛ τ+ δ e B ɛ ɛ τ+ ( ) τ (δ e B ) τ ɛ ɛ ɛ. δ e B Now we find a bound of X given by (7). Since f = A (e T 3 R T c Re 3 ) we have X A (e T 3 R T c Re 3 )Re 3 R c e 3 (k e + k v e v + B + δ ) (e T 3 Rc T Re 3 )Re 3 R c e 3. Using (63) and (77) the time-derivative of V is given by The last term (e T 3 Rc T Re 3 )Re 3 R c e 3 represents the sine of the angle between b 3 = Re 3 and b c3 = R c e 3 since (b 3c b 3 )b 3 b 3c = b 3 (b 3 b 3c ). The magnitude of the attitude error vector e R represents the sine of the eigen-ais rotation angle between R c and R (see [8). Therefore we have (e T 3 R T c Re 3 )Re 3 R c e 3 e R. It follows that (e T 3 Rd T Re 3 )Re 3 R d e 3 e R = Ψ(2 Ψ) ψ (2 ψ ) α <. Therefore X is bounded by X (k e + k v e v + B + δ ) e R (75) (k e + k v e v + B + δ )α. (76) Substituting (74) (76) into (73) V (k v ( α) c ) e v 2 c k m ( α) e 2 + c k v m ( + α) e e v + e R (B + δ )( c m e + e v ) + k e e v + ɛ. (77) In the above epression for V there is a third-order error term namely k e R e e v. Using (75) it is possible to choose its upper bound as k α e e v similar to other terms but the corresponding stability analysis becomes complicated and the initial attitude error should be reduced further. Instead we restrict our analysis to the domain D defined in (66) and its upper bound is chosen as k e ma e R e v. c) Lyapunov Candidate for the Complete System:: Let V = V + V 2 be the Lyapunov candidate of the complete system. V = 2 k e m e v 2 + c e e v + 2 e Ω Je Ω + k R Ψ(R R d ) + c 2 e R e Ω. (78) Using (67) the bound of the Lyapunov candidate V can be written as z T M z + z T 2 M 2 z 2 V z T M 2 z + z T 2 M 22z 2 (79) where z = [ e e v T z 2 = [ e R e Ω T R 2 and the matrices M M 2 M 2 M 22 are given by M = [ k c M 2 c m 2 = [ k c 2 c m M 2 = [ kr c 2 M 2 c 2 λ m (J) 22 = [ 2kR 2 ψ c 2. 2 c 2 λ M (J) V z T W z + z T W 2 z 2 z T 2 W 2 z 2 + ɛ + ɛ R (8) where W W 2 W 2 R 2 2 are defined as follows: [ ck ckv W = m ( α) 2m ( + α) ckv 2m ( + α) k (8) v( α) c [ c m (B + δ W 2 = ) (82) B + δ + k e ma W 2 = [ c2k R λ M (J) k Ω c 2. (83) d) Boundedness: Under the given conditions (35) (36) all of the matrices M M 2 M 2 M 22 W and W 2 are positive-definite. Therefore the Lyapunov function V is positive-definite and decrescent to obtain minλ m (M )λ m (M 2 ) z 2 V maλ M (M 2 ) λ M (M 22) z 2 (84) where z = [ z z 2 T R 2 and the time-derivative of V is bounded by V λ m (W ) z 2 + W 2 2 z z 2 λ m (W 2 ) z ɛ + ɛ R = z T W z + ɛ + ɛ R λ m (W ) z 2 + ɛ + ɛ R. (85) where the matri W R 2 2 is given by [ λm (W W = ) 2 W W. (86) 2 2 λ m (W 2 ) Similar to the proof of Proposition 2 we can show that the tracking errors are uniformly ultimately bounded if the constants ɛ ɛ R are sufficiently small as given in (38) and the corresponding ultimate bound is given by (39).
10 D. Proof of Proposition 4 The given assumptions satisfy the assumption of Proposition 2 from which the tracking error z 2 = [ e R e Ω T is guaranteed to eponentially decrease until it satisfies the bound given by (9). But (42) guarantees that the attitude error enters the region defined by (33) in a finite time t. Therefore if we show that the tracking error z = [ e e v T is bounded in t [ t as well then the complete tracking error (z z 2 ) is uniformly ultimately bounded. The boundedness of z is shown as follows. The error dynamics or e v can be written as mė v = mge 3 fre 3 mẍ d +. Let V 3 be a positive-definite function of e and e v : V 3 = 2 e m e v 2. Then we have e 2 2V 3 e v m V 3. The timederivative of V 3 is given by V 3 = e e v + e v (mge 3 fre 3 mẍ d + ) e e v + e v (B + δ ) + e v Re 3 f. From (27) we obtain V 3 e e v + e v (B + δ ) + e v (k e + k v e v + B + δ ) = k v e v 2 + (2(B + δ ) + (k + ) e ) e v d V 3 + d 2 V3 where d = k v 2 m + 2(k + ) m d 2 = 2(B + δ ) 2 m. Suppose that V 3 for a time interval [t a t b [ t. In this time interval we have V 3 V 3. Therefore V 3 (d + d 2 )V 3 V 3 (t) V 3 (t a )e (d+d2)(t ta). Therefore for any time interval in which V 3 V 3 is bounded. This implies that V 3 and therefore z = [ e e v T are bounded for t t. [7 G. Raffo M. Ortega and F. Rubio An integral predictive/nonlinear H control structure for a quadrotor helicopter Automatica vol. 46 pp [8 C. Nicol C. Macnab and A. Ramirez-Serrano Robust neural network control of a quadrotor helicopter in Proceedings of the Canadian Conference on Electrical and Computer Engineering 28 pp [9 A. Tayebi and S. McGilvray Attitude stabilization of a VTOL quadrotor aircraft IEEE Transactions on Control System Technology vol. 4 no. 3 pp [ C. Mayhew R. Sanfelice and A. Teel Quaternion-based hybrid control for robust global attitude tracking IEEE Transactions on Automatic Control 2. [ On the non-robustness of inconsistent quaternion-based attitude control systems using memoryless path-lifting schemes in Proceeding of the American Control Conference 2. [2 S. Bhat and D. Bernstein A topological obstruction to continuous global stabilization of rotational motion and the unwinding phenomenon Systems and Control Letters vol. 39 no. pp [3 C. Mayhew R. Sanfelice and A. Teel On quaternion-based attitude control and the unwinding phenomenon in Proceeding of the American Control Conference 2. [4 F. Bullo and A. Lewis Geometric control of mechanical systems ser. Tets in Applied Mathematics. New York: Springer-Verlag 25 vol. 49 modeling analysis and design for simple mechanical control systems. [5 D. Maithripala J. Berg and W. Dayawansa Almost global tracking of simple mechanical systems on a general class of Lie groups IEEE Transactions on Automatic Control vol. 5 no. pp [6 D. Cabecinhas R. Cunha and C. Silvestre Output-feedback control for almost global stabilization of fully-acuated rigid bodies in Proceedings of IEEE Conference on Decision and Control Ed. 28. [7 N. Chaturvedi A. Sanyal and N. McClamroch Rigid-body attitude control IEEE Control Systems Magazine vol. 3 no. 3 pp [8 T. Lee M. Leok and N. McClamroch Control of comple maneuvers for a quadrotor UAV using geometric methods on SE(3) arxiv. [Online. Available: [9 N. Chaturvedi N. H. McClamroch and D. Bernstein Asymptotic smooth stabilization of the inverted 3-D pendulum IEEE Transactions on Automatic Control vol. 54 no. 6 pp [2 T. Lee Robust adaptive geometric tracking controls on SO(3) with an application to the attitude dynamics of a quadrotor UAV arxiv 2. [Online. Available: [2 P. Pounds R. Mahony and P. Corke Modeling and control of a large quadrotor robot Control Engineering Practice vol. 8 pp [22 H. Khalil Nonlinear Systems 2nd Edition Ed. Prentice Hall 996. REFERENCES [ M. Valenti B. Bethke G. Fiore and J. How Indoor multi-vehicle flight testbed for fault detection indoor multi-vehicle flight testbed for fault detection isolation and recovery in Proceedings of the AIAA Guidance Navigation and Control Conference 26. [2 G. Hoffmann H. Huang S. Waslander and C. Tomlin Quadrotor helicopter flight dynamics and control: Theory and eperiment in Proceedings of the AIAA Guidance Navigation and Control Conference 27 AIAA [3 P. Castillo R. Lozano and A. Dzul Stabilization of a mini rotorcraft with four rotors IEEE Control System Magazine pp [4 J. Gillula G. Hoffmann H. Huang M. Vitus and C. Tomlin Applications of hybrid reachability analysis to robotic aerial vehicles The International Journal of Robotics Research vol. 3 no. 3 pp [5 S. Bouabdalla and R. Siegward Backstepping and sliding-mode techniques applied to an indoor micro quadrotor in Proceedings of the IEEE International Conference on Robotics and Automation 25 pp [6 M. Efe Robust low altitude behavior control of a quadrotor rotorcraft through sliding modes in Proceedings of the Mediterranean Conference on Control and Automation 27.
Nonlinear Robust Tracking Control of a Quadrotor UAV on SE(3)
22 American Control Conference Fairmont Queen Elizabeth Montréal Canada June 27-June 29 22 Nonlinear Robust Tracking Control of a Quadrotor UAV on SE(3) Taeyoung Lee Melvin Leok and N. Harris McClamroch
More informationGeometric Tracking Control of a Quadrotor UAV on SE(3)
49th IEEE Conference on Decision and Control December 5-7, 2 Hilton Atlanta Hotel, Atlanta, GA, USA Geometric Tracking Control of a Quadrotor UAV on SE(3) Taeyoung Lee, Melvin Leok, and N. Harris McClamroch
More informationGeometric Tracking Control of a Quadrotor UAV on SE(3)
Geometric Tracking Control o a Quadrotor UAV on SE(3) Taeyoung Lee, Melvin Leok, and N. Harris McClamroch Abstract This paper provides new results or the tracking control o a quadrotor unmanned aerial
More informationStable Manifolds of Saddle Equilibria for Pendulum Dynamics on S 2 and SO(3)
2011 50th IEEE Conference on Decision and Control and European Control Conference (CDC-ECC) Orlando, FL, USA, December 12-15, 2011 Stable Manifolds of Saddle Equilibria for Pendulum Dynamics on S 2 and
More informationPosition Control for a Class of Vehicles in SE(3)
Position Control for a Class of Vehicles in SE(3) Ashton Roza, Manfredi Maggiore Abstract A hierarchical design framework is presented to control the position of a class of vehicles in SE(3) that are propelled
More informationStabilization of a 3D Rigid Pendulum
25 American Control Conference June 8-, 25. Portland, OR, USA ThC5.6 Stabilization of a 3D Rigid Pendulum Nalin A. Chaturvedi, Fabio Bacconi, Amit K. Sanyal, Dennis Bernstein, N. Harris McClamroch Department
More informationAdaptive Control of a Quadrotor UAV Transporting a Cable-Suspended Load with Unknown Mass
rd IEEE Conference on Decision and Control December -7,. Los Angeles, California, USA Adaptive Control of a Quadrotor UAV Transporting a Cable-Suspended Load with Unknown Mass Shicong Dai, Taeyoung Lee,
More informationLow-Thrust Attitude Control for Nano-Satellite with Micro-Cathode Thrusters
Low-Thrust Attitude Control for Nano-Satellite with Micro-Cathode Thrusters IEPC-3-366 Presented at the 33 rd International Electric Propulsion Conference, The George Washington University, Washington,
More informationAn Intrinsic Robust PID Controller on Lie Groups
53rd IEEE Conference on Decision and Control December 15-17, 2014. Los Angeles, California, USA An Intrinsic Robust PID Controller on Lie Groups D.H.S. Maithripala and J. M. Berg Abstract This paper presents
More informationGlobal Formulations of Lagrangian and Hamiltonian Dynamics on Embedded Manifolds
1 Global Formulations of Lagrangian and Hamiltonian Dynamics on Embedded Manifolds By Taeyoung Lee, Melvin Leok, and N. Harris McClamroch Mechanical and Aerospace Engineering, George Washington University,
More informationNonlinear Landing Control for Quadrotor UAVs
Nonlinear Landing Control for Quadrotor UAVs Holger Voos University of Applied Sciences Ravensburg-Weingarten, Mobile Robotics Lab, D-88241 Weingarten Abstract. Quadrotor UAVs are one of the most preferred
More informationStabilization of a Specified Equilibrium in the Inverted Equilibrium Manifold of the 3D Pendulum
Proceedings of the 27 American Control Conference Marriott Marquis Hotel at Times Square New York City, USA, July 11-13, 27 ThA11.6 Stabilization of a Specified Equilibrium in the Inverted Equilibrium
More informationNonlinear Tracking Control of Underactuated Surface Vessel
American Control Conference June -. Portland OR USA FrB. Nonlinear Tracking Control of Underactuated Surface Vessel Wenjie Dong and Yi Guo Abstract We consider in this paper the tracking control problem
More informationGeometric Mechanics and Global Nonlinear Control for Multi-Body Dynamics
Geometric Mechanics and Global Nonlinear Control for Multi-Body Dynamics Harris McClamroch Aerospace Engineering, University of Michigan Joint work with Taeyoung Lee (George Washington University) Melvin
More informationExperimental Results for Almost Global Asymptotic and Locally Exponential Stabilization of the Natural Equilibria of a 3D Pendulum
Proceedings of the 26 American Control Conference Minneapolis, Minnesota, USA, June 4-6, 26 WeC2. Experimental Results for Almost Global Asymptotic and Locally Exponential Stabilization of the Natural
More informationVelocity-Free Hybrid Attitude Stabilization Using Inertial Vector Measurements
016 American Control Conference (ACC) Boston Marriott Copley Place July 6-8, 016. Boston, MA, USA Velocity-Free Hybrid Attitude Stabilization Using Inertial Vector Measurements Soulaimane Berkane and Abdelhamid
More informationNonlinear Control of a Quadrotor Micro-UAV using Feedback-Linearization
Proceedings of the 2009 IEEE International Conference on Mechatronics. Malaga, Spain, April 2009. Nonlinear Control of a Quadrotor Micro-UAV using Feedback-Linearization Holger Voos University of Applied
More informationGlobal Trajectory Tracking for Underactuated VTOL Aerial Vehicles using a Cascade Control Paradigm
Global Trajectory Tracking for Underactuated VTOL Aerial Vehicles using a Cascade Control Paradigm Roberto Naldi, Michele Furci, Ricardo G. Sanfelice and Lorenzo Marconi Abstract This work proposes a feedback
More informationarxiv: v2 [math.oc] 17 May 2018
Geometric Adaptive Control for a Quadrotor UAV with Wind Disturbance Rejection Mahdis Bisheban, and aeyoung Lee arxiv:803.06363v [math.oc 7 May 08 Abstract his paper presents a geometric adaptive control
More informationUnit quaternion observer based attitude stabilization of a rigid spacecraft without velocity measurement
Proceedings of the 45th IEEE Conference on Decision & Control Manchester Grand Hyatt Hotel San Diego, CA, USA, December 3-5, 6 Unit quaternion observer based attitude stabilization of a rigid spacecraft
More informationAlmost Global Robust Attitude Tracking Control of Spacecraft in Gravity
Almost Global Robust Attitude Tracking Control of Spacecraft in Gravity Amit K. Sanyal Nalin A. Chaturvedi In this paper, we treat the practical problem of tracking the attitude and angular velocity of
More information3D Pendulum Experimental Setup for Earth-based Testing of the Attitude Dynamics of an Orbiting Spacecraft
3D Pendulum Experimental Setup for Earth-based Testing of the Attitude Dynamics of an Orbiting Spacecraft Mario A. Santillo, Nalin A. Chaturvedi, N. Harris McClamroch, Dennis S. Bernstein Department of
More informationQUADROTOR: FULL DYNAMIC MODELING, NONLINEAR SIMULATION AND CONTROL OF ATTITUDES
QUADROTOR: FULL DYNAMIC MODELING, NONLINEAR SIMULATION AND CONTROL OF ATTITUDES Somayeh Norouzi Ghazbi,a, Ali Akbar Akbari 2,a, Mohammad Reza Gharib 3,a Somaye_noroozi@yahoo.com, 2 Akbari@um.ac.ir, 3 mech_gharib@yahoo.com
More informationFurther results on global stabilization of the PVTOL aircraft
Further results on global stabilization of the PVTOL aircraft Ahmad Hably, Farid Kendoul 2, Nicolas Marchand, and Pedro Castillo 2 Laboratoire d Automatique de Grenoble, ENSIEG BP 46, 3842 Saint Martin
More informationGeometric Control and Differential Flatness of a Quadrotor UAV with a Cable-Suspended Load
Geometric Control and Differential Flatness of a Quadrotor UAV with a Cable-Suspended Load Koushil Sreenath, Taeyoung Lee, Vijay Kumar x Q R 3,R SO(3) Abstract A quadrotor with a cable-suspended load with
More informationGeometric Backstepping for Strict Feedback Simple Mechanical Systems
Geometric Backstepping for Strict Feedback Simple Mechanical Systems Guofan Wu and Koushil Sreenath Abstract We propose a geometric control design method for a class of underactuated geometric mechanical
More informationRigid-Body Attitude Control USING ROTATION MATRICES FOR CONTINUOUS, SINGULARITY-FREE CONTROL LAWS
Rigid-Body Attitude Control USING ROTATION MATRICES FOR CONTINUOUS, SINGULARITY-FREE CONTROL LAWS 3 IEEE CONTROL SYSTEMS MAGAZINE» JUNE -33X//$. IEEE SHANNON MASH Rigid-body attitude control is motivated
More informationROBUST SECOND ORDER SLIDING MODE CONTROL
ROBUST SECOND ORDER SLIDING MODE CONTROL FOR A QUADROTOR CONSIDERING MOTOR DYNAMICS Nader Jamali Soufi Amlashi 1, Mohammad Rezaei 2, Hossein Bolandi 2 and Ali Khaki Sedigh 3 1 Department of Control Engineering,
More informationMathematical Modelling and Dynamics Analysis of Flat Multirotor Configurations
Mathematical Modelling and Dynamics Analysis of Flat Multirotor Configurations DENIS KOTARSKI, Department of Mechanical Engineering, Karlovac University of Applied Sciences, J.J. Strossmayera 9, Karlovac,
More informationIntrinsic UDE Control of Mechanical Systems on SO(3)
Margin requirements for first page 72 pt 1 in 25.4 mm Intrinsic UDE Control of Mechanical Systems on SO(3) A. L. M. Sanka Liyanage, Jordan M. Berg, Beibei Ren, and D. H. S. Maithripala Abstract The uncertainty
More informationSemi-Global Attitude Controls Bypassing the Topological Obstruction on SO(3)
Semi-Global Attitude Controls Bypassing the Topological Obstruction on SO(3) Taeyoung Lee, Dong Eui Chang and Yongsoon Eun arxiv:78.7649v [math.oc] 5 Aug 7 Abstract This paper presents global tracking
More informationModeling and Sliding Mode Control of a Quadrotor Unmanned Aerial Vehicle
Modeling and Sliding Mode Control of a Quadrotor Unmanned Aerial Vehicle Nour BEN AMMAR, Soufiene BOUALLÈGUE and Joseph HAGGÈGE Research Laboratory in Automatic Control LA.R.A), National Engineering School
More informationHover Control for Helicopter Using Neural Network-Based Model Reference Adaptive Controller
Vol.13 No.1, 217 مجلد 13 العدد 217 1 Hover Control for Helicopter Using Neural Network-Based Model Reference Adaptive Controller Abdul-Basset A. Al-Hussein Electrical Engineering Department Basrah University
More informationSYMMETRIES OF ROCKET DYNAMICS. Part 1. Executive Summary 1 1. Equations of Motion 1 2. The Problem to Solve 2. Part 2. Introduction 2 3.
SYMMETRIES OF ROCKET DYNAMICS ERNEST YEUNG Abstract. I recap the Euler-Poincaré reduction for the Lagrangian of a free (and non-free) rigid body and discuss the constraints and time steps for the Time
More informationSimulation of Backstepping-based Nonlinear Control for Quadrotor Helicopter
APPLICATIONS OF MODELLING AND SIMULATION http://amsjournal.ams-mss.org eissn 2680-8084 VOL 2, NO. 1, 2018, 34-40 Simulation of Backstepping-based Nonlinear Control for Quadrotor Helicopter M.A.M. Basri*,
More informationQuadrotor Modeling and Control
16-311 Introduction to Robotics Guest Lecture on Aerial Robotics Quadrotor Modeling and Control Nathan Michael February 05, 2014 Lecture Outline Modeling: Dynamic model from first principles Propeller
More informationAdaptive Robust Control (ARC) for an Altitude Control of a Quadrotor Type UAV Carrying an Unknown Payloads
2 th International Conference on Control, Automation and Systems Oct. 26-29, 2 in KINTEX, Gyeonggi-do, Korea Adaptive Robust Control (ARC) for an Altitude Control of a Quadrotor Type UAV Carrying an Unknown
More informationControllability Analysis of A Two Degree of Freedom Nonlinear Attitude Control System
Controllability Analysis of A Two Degree of Freedom Nonlinear Attitude Control System Jinglai Shen, Amit K Sanyal, and N Harris McClamroch Department of Aerospace Engineering The University of Michigan
More informationTrajectory tracking & Path-following control
Cooperative Control of Multiple Robotic Vehicles: Theory and Practice Trajectory tracking & Path-following control EECI Graduate School on Control Supélec, Feb. 21-25, 2011 A word about T Tracking and
More informationIAC-11-C1.5.9 INERTIA-FREE ATTITUDE CONTROL OF SPACECRAFT WITH UNKNOWN TIME-VARYING MASS DISTRIBUTION
6nd International Astronautical Congress, Cape Town, SA. Copyright by the International Astronautical Federation. All rights reserved IAC--C.5.9 INERTIA-FREE ATTITUDE CONTROL OF SPACECRAFT WITH UNKNOWN
More informationA Nonlinear Control Law for Hover to Level Flight for the Quad Tilt-rotor UAV
Preprints of the 19th World Congress The International Federation of Automatic Control A Nonlinear Control Law for Hover to Level Flight for the Quad Tilt-rotor UAV Gerardo R. Flores-Colunga Rogelio Lozano-Leal
More informationPath Following Controller for a Quadrotor Helicopter
Path Following Controller for a Quadrotor Helicopter Ashton Roza, Manfredi Maggiore Abstract A path following controller is presented for a quadrotor helicopter model. The controller relies on input dynamic
More informationarxiv: v2 [cs.sy] 28 Mar 2017
Attitude Tracking Control for Aerobatic Helicopters: A Geometric Approach Nidhish Raj 1, Ravi N. Banavar, Abhishek 3, and Mangal Kothari 3 arxiv:173.88v [cs.sy] 8 Mar 17 Abstract We consider the problem
More informationControl and Navigation Framework for Quadrotor Helicopters
DOI 1.17/s1846-1-9789-z Control and Navigation Framework for Quadrotor Helicopters Amr Nagaty Sajad Saeedi Carl Thibault Mae Seto Howard Li Received: September 1 / Accepted: September 1 Springer Science+Business
More informationNonlinear and Neural Network-based Control of a Small Four-Rotor Aerial Robot
Nonlinear and Neural Network-based Control of a Small Four-Rotor Aerial Robot Holger Voos Abstract Small four-rotor aerial robots, so called quadrotor UAVs, have an enormous potential for all kind of neararea
More informationROBUST NEURAL NETWORK CONTROL OF A QUADROTOR HELICOPTER. Schulich School of Engineering, University of Calgary
ROBUST NEURAL NETWORK CONTROL OF A QUADROTOR HELICOPTER C. Nicol,C.J.B. Macnab, A. Ramirez-Serrano Schulich School of Engineering, University of Calgary Department of Electrical and Computer Engineering
More informationDynamic Modeling and Stabilization Techniques for Tri-Rotor Unmanned Aerial Vehicles
Technical Paper Int l J. of Aeronautical & Space Sci. 11(3), 167 174 (010) DOI:10.5139/IJASS.010.11.3.167 Dynamic Modeling and Stabilization Techniques for Tri-Rotor Unmanned Aerial Vehicles Dong-Wan Yoo*,
More informationThe PVTOL Aircraft. 2.1 Introduction
2 The PVTOL Aircraft 2.1 Introduction We introduce in this chapter the well-known Planar Vertical Take-Off and Landing (PVTOL) aircraft problem. The PVTOL represents a challenging nonlinear systems control
More informationDynamics and Control of a Quadrotor with a Payload Suspended through an Elastic Cable
Dynamics and Control of a Quadrotor with a Payload Suspended through an Elastic Cable Prasanth Kotaru, Guofan Wu and Koushil Sreenath Abstract We study the problem of a quadrotor transporting a payload
More informationAttitude and Earth Velocity Estimation - Part I: Globally Exponentially Stable Observer
Attitude Earth Velocity Estimation - Part I: Globally Eponentially Stable Observer Pedro Batista, Carlos Silvestre, Paulo Oliveira Abstract The main contribution of this paper is the development of a novel
More informationQuadrotors Flight Formation Control Using a Leader-Follower Approach*
23 European Conference (ECC) July 7-9, 23, Zürich, Switzerland. Quadrotors Flight Formation Using a Leader-Follower Approach* D. A. Mercado, R. Castro and R. Lozano 2 Abstract In this paper it is presented
More informationDigital Passive Attitude and Altitude Control Schemes for Quadrotor Aircraft
Institute for Software Integrated Systems Vanderbilt University Nashville, Tennessee, 37235 Digital Passive Attitude and Altitude Control Schemes for Quadrotor Aircraft Nicholas Kottenstette, and Joseph
More informationwith Application to Autonomous Vehicles
Nonlinear with Application to Autonomous Vehicles (Ph.D. Candidate) C. Silvestre (Supervisor) P. Oliveira (Co-supervisor) Institute for s and Robotics Instituto Superior Técnico Portugal January 2010 Presentation
More informationAttitude Regulation About a Fixed Rotation Axis
AIAA Journal of Guidance, Control, & Dynamics Revised Submission, December, 22 Attitude Regulation About a Fixed Rotation Axis Jonathan Lawton Raytheon Systems Inc. Tucson, Arizona 85734 Randal W. Beard
More informationChapter 2 Review of Linear and Nonlinear Controller Designs
Chapter 2 Review of Linear and Nonlinear Controller Designs This Chapter reviews several flight controller designs for unmanned rotorcraft. 1 Flight control systems have been proposed and tested on a wide
More informationAN INTEGRATOR BACKSTEPPING CONTROLLER FOR A STANDARD HELICOPTER YITAO LIU THESIS
AN INEGRAOR BACKSEPPING CONROLLER FOR A SANDARD HELICOPER BY YIAO LIU HESIS Submitted in partial fulfillment of the requirements for the degree of Master of Science in Electrical and Computer Engineering
More informationGeometric path following control of a rigid body based on the stabilization of sets
Preprints of the 19th World Congress The International Federation of Automatic Control Geometric path following control of a rigid body based on the stabilization of sets uri A. Kapitanyuk Sergey A. Chepinskiy
More informationPassivity-based Formation Control for UAVs with a Suspended Load
Passivity-based Formation Control for UAVs with a Suspended Load Chris Meissen Kristian Klausen Murat Arcak Thor I. Fossen Andrew Packard Department of Mechanical Engineering at the University of California,
More informationEnergy-based Swing-up of the Acrobot and Time-optimal Motion
Energy-based Swing-up of the Acrobot and Time-optimal Motion Ravi N. Banavar Systems and Control Engineering Indian Institute of Technology, Bombay Mumbai-476, India Email: banavar@ee.iitb.ac.in Telephone:(91)-(22)
More informationAn Adaptive Block Backstepping Controller for Attitude Stabilization of a Quadrotor Helicopter
WSEAS RANSACIONS on SYSES and CONROL An Adaptive Block Backstepping Controller for Attitude Stabilization of a Quadrotor Helicopter Unmanned Aerial Vehicle Engineering Department echanical Engineering
More informationQuadcopter Dynamics 1
Quadcopter Dynamics 1 Bréguet Richet Gyroplane No. 1 1907 Brothers Louis Bréguet and Jacques Bréguet Guidance of Professor Charles Richet The first flight demonstration of Gyroplane No. 1 with no control
More informationRobust Global Asymptotic Attitude Stabilization of a Rigid Body by Quaternion-based Hybrid Feedback
Robust Global Asymptotic Attitude Stabilization of a Rigid Body by Quaternion-based Hybrid Feedback Christopher G. Mayhew, Ricardo G. Sanfelice, and Andrew R. Teel Abstract Global asymptotic stabilization
More informationGlobal Trajectory Tracking for a Class of Underactuated Vehicles
Global Trajectory Tracking for a Class of Underactuated Vehicles Pedro Casau, Ricardo G. Sanfelice, Rita Cunha, David Cabecinhas, Carlos Silvestre Abstract In this paper, we address the problem of trajectory
More informationModelling of Opposed Lateral and Longitudinal Tilting Dual-Fan Unmanned Aerial Vehicle
Modelling of Opposed Lateral and Longitudinal Tilting Dual-Fan Unmanned Aerial Vehicle N. Amiri A. Ramirez-Serrano R. Davies Electrical Engineering Department, University of Calgary, Canada (e-mail: namiri@ucalgary.ca).
More informationMulti-layer Flight Control Synthesis and Analysis of a Small-scale UAV Helicopter
Multi-layer Flight Control Synthesis and Analysis of a Small-scale UAV Helicopter Ali Karimoddini, Guowei Cai, Ben M. Chen, Hai Lin and Tong H. Lee Graduate School for Integrative Sciences and Engineering,
More informationarxiv: v1 [cs.sy] 29 Aug 2017
Quaternions and Attitude Representation Hardik Parwana, Mangal Kothari Department of Aerospace Engineering, Indian Institute of Technology Kanpur, India 1 Attitude Representation arxiv:178.868v1 [cs.sy]
More informationInvestigation of the Dynamics and Modeling of a Triangular Quadrotor Configuration
Investigation of the Dynamics and Modeling of a Triangular Quadrotor Configuration TONI AXELSSON Master s Thesis at Aerospace Engineering Supervisor: Arne Karlsson Examiner: Arne Karlsson ISSN 1651-7660
More informationENHANCED PROPORTIONAL-DERIVATIVE CONTROL OF A MICRO QUADCOPTER
ENHANCED PROPORTIONAL-DERIVATIVE CONTROL OF A MICRO QUADCOPTER Norman L. Johnson and Kam K. Leang Department of Mechanical Engineering University of Nevada, Reno Reno, Nevada 897-312, USA ABSTRACT This
More informationRobust Hybrid Global Asymptotic Stabilization of Rigid Body Dynamics using Dual Quaternions
Robust Hybrid Global Asymptotic Stabilization of Rigid Body Dynamics using Dual Quaternions Bharani P. Malladi and Eric A. Butcher, University of Arizona, Tucson, AZ, 872, USA Ricardo G. Sanfelice University
More informationAROTORCRAFT-BASED unmanned aerial vehicle
1392 IEEE TRANSACTIONS ON CONTROL SYSTEMS TECHNOLOGY, VOL. 20, NO. 5, SEPTEMBER 2012 Autonomous Flight of the Rotorcraft-Based UAV Using RISE Feedback and NN Feedforward Terms Jongho Shin, H. Jin Kim,
More informationOn Quaternion-Based Attitude Control and the Unwinding Phenomenon
2011 American Control Conference on O'Farrell Street, San Francisco, CA, USA June 29 - July 01, 2011 On Quaternion-Based Attitude Control and the Unwinding Phenomenon Christopher G. Mayhew, Ricardo G.
More informationNomenclature. = the local vertical and local horizontal reference frame. = angular velocity vector of spacecraft (rad/s)
Adaptive Fault-Tolerant Control of Spacecraft Attitude Dynamics with Actuator Failures Yu Han 1 Harbin Institude of Technology, Harbin, Heilongjiang, China, 151 James D. Biggs 2 University of Strathclyde,
More informationRobot Control Basics CS 685
Robot Control Basics CS 685 Control basics Use some concepts from control theory to understand and learn how to control robots Control Theory general field studies control and understanding of behavior
More informationExperiments on Stabilization of the Hanging Equilibrium of a 3D Asymmetric Rigid Pendulum
Proceedings of the 25 IEEE Conference on Control Applications Toronto, Canada, August 28-3, 25 MB4.5 Experiments on Stabilization of the Hanging Equilibrium of a 3D Asymmetric Rigid Pendulum Mario A. Santillo,
More informationKinematics. Chapter Multi-Body Systems
Chapter 2 Kinematics This chapter first introduces multi-body systems in conceptual terms. It then describes the concept of a Euclidean frame in the material world, following the concept of a Euclidean
More informationADAPTIVE SLIDING MODE CONTROL OF UNMANNED FOUR ROTOR FLYING VEHICLE
International Journal of Robotics and Automation, Vol. 30, No. 2, 205 ADAPTIVE SLIDING MODE CONTROL OF UNMANNED FOUR ROTOR FLYING VEHICLE Shafiqul Islam, Xiaoping P. Liu, and Abdulmotaleb El Saddik Abstract
More informationCLF-based Tracking Control for UAV Kinematic Models with Saturation Constraints
CDC3-IEEE45 CLF-based Tracking Control for UAV Kinematic Models with Saturation Constraints Wei Ren Randal W. Beard Department of Electrical and Computer Engineering Brigham Young University Provo, Utah
More informationarxiv:math/ v1 [math.oc] 17 Sep 2006
Global Attitude Estimation using Single Direction Measurements Taeyoung Lee, Melvin Leok, N. Harris McClamroch, and Amit Sanyal arxiv:math/0609481v1 [math.oc] 17 Sep 2006 Abstract A deterministic attitude
More informationControl of a Quadrotor Mini-Helicopter via Full State Backstepping Technique
Proceedings of the 45th IEEE Conference on Decision & Control Manchester Grand Hyatt Hotel San Diego, CA, USA, December 3-5, 006 Control of a Quadrotor Mini-Helicopter via Full State Backstepping Technique
More informationAdaptive Trim and Trajectory Following for a Tilt-Rotor Tricopter Ahmad Ansari, Anna Prach, and Dennis S. Bernstein
7 American Control Conference Sheraton Seattle Hotel May 4 6, 7, Seattle, USA Adaptive Trim and Trajectory Following for a Tilt-Rotor Tricopter Ahmad Ansari, Anna Prach, and Dennis S. Bernstein Abstract
More informationBounded attitude control of rigid bodies: Real-time experimentation to a quadrotor mini-helicopter
Bounded attitude control of rigid bodies: Real-time experimentation to a quadrotor mini-helicopter Fermi Guerrero-Castellanos, Nicolas Marchand, Ahmad Hably, Suzanne Lesecq, Jérôme Delamare To cite this
More informationMini-quadrotor Attitude Control based on Hybrid Backstepping & Frenet-Serret Theory
Mini-quadrotor Attitude Control based on Hybrid Backstepping & Frenet-Serret Theory J. Colorado, A. Barrientos, Senior Member, IEEE, A. Martinez, B. Lafaverges, and J. Valente Abstract This paper is about
More informationReal-time Motion Control of a Nonholonomic Mobile Robot with Unknown Dynamics
Real-time Motion Control of a Nonholonomic Mobile Robot with Unknown Dynamics TIEMIN HU and SIMON X. YANG ARIS (Advanced Robotics & Intelligent Systems) Lab School of Engineering, University of Guelph
More informationCONTROL OF THE NONHOLONOMIC INTEGRATOR
June 6, 25 CONTROL OF THE NONHOLONOMIC INTEGRATOR R. N. Banavar (Work done with V. Sankaranarayanan) Systems & Control Engg. Indian Institute of Technology, Bombay Mumbai -INDIA. banavar@iitb.ac.in Outline
More informationOptimal Fault-Tolerant Configurations of Thrusters
Optimal Fault-Tolerant Configurations of Thrusters By Yasuhiro YOSHIMURA ) and Hirohisa KOJIMA, ) ) Aerospace Engineering, Tokyo Metropolitan University, Hino, Japan (Received June st, 7) Fault tolerance
More informationDifferent Approaches of PID Control UAV Type Quadrotor
Different Approaches of PD Control UAV ype Quadrotor G. Szafranski, R. Czyba Silesian University of echnology, Akademicka St 6, Gliwice, Poland ABSRAC n this paper we focus on the different control strategies
More informationAerial Robotics. Vision-based control for Vertical Take-Off and Landing UAVs. Toulouse, October, 2 nd, Henry de Plinval (Onera - DCSD)
Aerial Robotics Vision-based control for Vertical Take-Off and Landing UAVs Toulouse, October, 2 nd, 2014 Henry de Plinval (Onera - DCSD) collaborations with P. Morin (UPMC-ISIR), P. Mouyon (Onera), T.
More informationMotion control of rigid bodies in SE(3) Ashton Roza
Motion control of rigid bodies in SE(3) by Ashton Roza A thesis submitted in conformity with the requirements for the degree of Master of Applied Science Graduate Department of Electrical and Computer
More informationUAV Rotorcraft in Compliant Contact: Stability Analysis and Simulation
2011 IEEE/RSJ International Conference on Intelligent Robots and Systems September 25-30, 2011. San Francisco, CA, USA UAV Rotorcraft in Compliant Contact: Stability Analysis and Simulation Paul E. I.
More informationDiscontinuous Backstepping for Stabilization of Nonholonomic Mobile Robots
Discontinuous Backstepping for Stabilization of Nonholonomic Mobile Robots Herbert G. Tanner GRASP Laboratory University of Pennsylvania Philadelphia, PA, 94, USA. tanner@grasp.cis.upenn.edu Kostas J.
More informationRigid Body Dynamics, SG2150 Solutions to Exam,
KTH Mechanics 011 10 Calculational problems Rigid Body Dynamics, SG150 Solutions to Eam, 011 10 Problem 1: A slender homogeneous rod of mass m and length a can rotate in a vertical plane about a fied smooth
More informationAsymptotic Smooth Stabilization of the Inverted 3D Pendulum
IEEE TRANSACTIONS ON AUTOMATIC CONTROL 1 Asymptotic Smooth Stabilization of the Inverted 3D Pendulum Nalin A. Chaturvedi, N. Harris McClamroch, Fellow, IEEE, and Dennis S. Bernstein, Fellow, IEEE Abstract
More informationNonlinear Control of a Spacecraft with Multiple Fuel Slosh Modes
11 5th IEEE Conference on Decision and Control and European Control Conference (CDC-ECC) Orlando, FL, USA, December 1-15, 11 onlinear Control of a Spacecraft with Multiple Fuel Slosh Modes Mahmut Reyhanoglu
More informationThe Variational Attitude Estimator in the Presence of Bias in Angular Velocity Measurements
The Variational Attitude Estimator in the Presence of Bias in Angular Velocity Measurements Maziar Izadi 1, Sasi Prabhakaran 1, Amit Sanyal 2,, Carlos Silvestre 3, and Paulo Oliveira 4 Abstract Estimation
More informationA Model-Free Control System Based on the Sliding Mode Control Method with Applications to Multi-Input-Multi-Output Systems
Proceedings of the 4 th International Conference of Control, Dynamic Systems, and Robotics (CDSR'17) Toronto, Canada August 21 23, 2017 Paper No. 119 DOI: 10.11159/cdsr17.119 A Model-Free Control System
More informationAdaptive Robust Tracking Control of Robot Manipulators in the Task-space under Uncertainties
Australian Journal of Basic and Applied Sciences, 3(1): 308-322, 2009 ISSN 1991-8178 Adaptive Robust Tracking Control of Robot Manipulators in the Task-space under Uncertainties M.R.Soltanpour, M.M.Fateh
More informationLoad transportation using rotary-wing UAVs
Load transportation using rotary-wing UAVs Rafael José Figueiras dos Santos rafael.j.f.santos@tecnico.ulisboa.pt Instituto Superior Técnico, Lisboa, Portugal December 5 Abstract The problem of slung load
More informationAdaptive Backstepping Control for Optimal Descent with Embedded Autonomy
Adaptive Backstepping Control for Optimal Descent with Embedded Autonomy Maodeng Li, Wuxing Jing Department of Aerospace Engineering, Harbin Institute of Technology, Harbin, Heilongjiang, 150001, China
More informationTERMINAL ATTITUDE-CONSTRAINED GUIDANCE AND CONTROL FOR LUNAR SOFT LANDING
IAA-AAS-DyCoSS2-14 -02-05 TERMINAL ATTITUDE-CONSTRAINED GUIDANCE AND CONTROL FOR LUNAR SOFT LANDING Zheng-Yu Song, Dang-Jun Zhao, and Xin-Guang Lv This work concentrates on a 3-dimensional guidance and
More informationComputational Geometric Uncertainty Propagation for Hamiltonian Systems on a Lie Group
Computational Geometric Uncertainty Propagation for Hamiltonian Systems on a Lie Group Melvin Leok Mathematics, University of California, San Diego Foundations of Dynamics Session, CDS@20 Workshop Caltech,
More information