RECEDING HORIZON DRIFT COUNTERACTION AND ITS APPLICATION TO SPACECRAFT ATTITUDE CONTROL
|
|
- Monica Summers
- 5 years ago
- Views:
Transcription
1 AAS 7-46 RECEDING HORIZON DRIFT COUNTERACTION AND ITS APPLICATION TO SPACECRAFT ATTITUDE CONTROL Robert A. E. Zidek, Christopher D. Petersen, Alberto Bemporad, and Ilya V. Kolmanovsky In this paper, a recently developed model predictive control (MPC) approach for drift counteraction optimal control (DCOC) is applied to attitude control of fully actuated and underactuated spacecraft with reaction wheels (RWs). The objective is to maximize the time that prescribed constraints on spacecraft orientation and RW spin rates are satisfied given disturbance torques due to solar radiation pressure. While the MPC/DCOC approach is based on linear programming, all closedloop simulations are performed using the nonlinear model. In case constraints are violated, a control strategy is presented that recovers constraint satisfaction (if possible). We consider the cases where either one, two, or three RWs are operable. The numerical results show that the proposed controller successfully counteracts drift in order to satisfy constraints for as long as possible. INTRODUCTION The objective of drift counteraction optimal control (DCOC) is to satisfy prescribed state and control constraints for as long as possible. Such problems arise in many space flight applications, in particular, when finite resources (fuel/energy) are of concern. For example, control strategies have been developed to extend the operational time of a satellite subject to higher-order orbital perturbations. 4 Similar problems may be formulated for spacecraft attitude control, where tight constraints on spacecraft orientation and/or angular velocity need to be satisfied for as long as possible when the controllability of the spacecraft is limited. This may be the case for underactuated spacecraft,, 6 examples of which are the Kepler spacecraft, which lost two out of its four reaction wheels (RWs), 7 and two out of three RWs of the Hayabusa spacecraft failed. 8 Furthermore, reduced controllability may result from RW saturation. In this paper, a DCOC problem of attitude control of a spacecraft with RWs is considered. The spacecraft orientation is subject to drift caused by solar radiation pressure (SRP). We use a recently developed DCOC approach 9 to maximize the time that prescribed constraints on spacecraft orientation and RW spin rates are satisfied. The approach assumes linear dynamics and solutions to the linear DCOC problem are obtained by either employing linear programming (LP) or mixed integer linear programming (MILP). While the MILP solution is guaranteed to be optimal with respect to the linear DCOC problem, LP efficiently obtains good-quality suboptimal solutions. 9 We present closed-loop simulations based on nonlinear dynamics, where a receding horizon / model predictive control (MPC) implementation of the LP formulation is used to provide state feedback control in Department of Aerospace Engineering, University of Michigan, Ann Arbor, MI 489, USA. IMT Institute for Advanced Studies, Lucca, Italy.
2 order to compensate for unmodeled effects. Furthermore, an additional control strategy is introduced that tries to recover constraint satisfaction in case constraints are violated. We also show that the proposed approach is able to satisfy constraints for an arbitrarily long time when no solution to the DCOC problem exists, i.e., when the spacecraft has enough control authority to satisfy the constraints indefinitely. The paper is structured as follows. First, the spacecraft model is described, including the nonlinear spacecraft kinematics and dynamics, the SRP torque, and the linearized spacecraft model which is used to obtain the DCOC/MPC controller. Then we state the DCOC problem (for linear dynamics) and formulate the LP and MILP to solve it. Moreover, the MPC implementation is described and a control strategy to recover constraint satisfaction is introduced. Numerical results are shown in the subsequent section, where the cases of either one, two, or three operable RWs are treated. For the case of three operable RWs, we also show that the proposed approach may satisfy constraints indefinitely depending on the initial condition and constraints. A conclusion is given in the last section. SPACECRAFT MODEL The spacecraft model, adopted from Petersen et al., 6 is summarized in this section. Notation The following notation is used for the derivation of the spacecraft model. A physical vector is denoted by r and a physical unit vector is denoted by ˆr. Resolving r or ˆr in a given frame F results in a mathematical vector denoted by r F or ˆr F, respectively. Two frames are considered, a given inertial reference frame denoted by I and the body-fixed frame B, which is assumed to be a principal frame. The notation for r or ˆr resolved in the body-fixed frame is r = r B or r = ˆr B, respectively. The skew-symmetric matrix of a mathematical vector r = [r, r, r 3 ] T is given by Kinematics r 3 r S[ r] = r 3 r. () r r The orientation of the spacecraft s body-fixed frame B relative to I is described by Θ = [φ, θ, ψ] T, where ψ (yaw), θ (pitch), and φ (roll) are the 3-- Euler angles. Then the kinematic equations read Θ = cos(θ) cos(θ) sin(φ) sin(θ) cos(φ) sin(θ) cos(φ) cos(θ) sin(φ) cos(θ) sin(φ) cos(θ) where ω is the angular velocity of B relative to I and ω = [ω, ω, ω 3 ] T. Dynamics ω, () The spacecraft is equipped with up to three RWs, where ḡ i denotes the unit vector of the ith RW spin axis resolved in the B frame. Moreover, let W = [ḡ, ḡ, ḡ 3 ] in case of three RWs and W = ḡ or W = [ḡ, ḡ ] in case of one or two operable RWs, respectively. The spin rate of the ith RW is ν i and ν = [ν, ν, ν 3 ] T, [ν, ν ] T, or ν, respectively. We assume that all RWs are identical and thin
3 (moments of inertia about axes transversal to spin axis are approximately zero), where J w denotes the moment of inertia about the RW spin axis. The moment of inertia matrix of the spacecraft bus resolved in the B frame is given by J = diag(j, J, J 3 ) and J = J + J w W W T (3) is the locked inertia. The rotational dynamics of the spacecraft described in the B frame are given by 6 J ω + S[ ω]( J ω + J w W ν) + J w W ν = τ srp, (4) where τ srp is an external torque due to SRP. SRP Torque For modeling the SRP torque in (4), a cuboid spacecraft (six flat panels) is assumed. We define κ = Φ sun,tot c(d SC,sun /d ), β = 4 9 C diff, () where c is the speed of light, Φ sun,tot is the solar flux experienced by the spacecraft, d = AU is the nominal distance between Earth and Sun, d SC,sun is the distance between the spacecraft and Sun (d SC,sun =.99 AU in this paper), and C diff denotes the diffusion coefficient, which is assumed to be the same for all panels. Assuming the SRP acts identically across all points on the jth panel, the SRP acting on panel j is given by6, P j = κ(ˆq j ˆq sun )(ˆq j + β ˆq sun ), (6) where ˆq j is the normal to the surface of the jth panel pointing outward from the spacecraft and ˆq sun represents the direction of the sun. Consequently, the SRP torque due to the jth panel resolved in the B frame is τ srp,j = S[ r j/o r C/O ]A j Pj, (7) where r C/O = [l x, l y, l z ] T denotes the mathematical vector of the position of the spacecraft s center of mass, C, relative to the geometric center, O, of the cuboid. The vector from O to the geometric center of the jth panel is represented by r j/o. With j {x+, x, y+, y, z+, z }, it is r x+/o = r x /O = [L x /,, ] T, r y+/o = r y /O = [, L y /, ] T, r z+/o = r z /O = [,, L z /] T, where the surface areas of the panels are A x+ = A x = L y L z, A y+ = A y = L x L z, and A z+ = A z = L x L y. The total SRP torque is the sum of all panel contributions, τ srp = 6 τ srp,j I j, (8) j= with I j as a function indicating if the jth panel is facing the Sun or not, i.e., I j = if ˆq j ˆq sun > and I j = otherwise. 3
4 Linear Discrete-Time Model The continuous-time nonlinear spacecraft model given by () and (4) is linearized and transformed into a discrete-time model in order to obtain the MPC drift counteraction controller (see subsequent section). We define the state and control input vectors, respectively, by x = [φ, θ, ψ, ω, ω, ω 3, ν,..., ν m ] T, u = [α,..., α m ] T, (9) with m {,, 3} as the number of operable RWs and α i = ν i as the acceleration of the ith RW. The linearization point is given by x lin = [ 6, ν T]T and u lin = m. Thus, the matrices corresponding to the continuous-time state space model, ẋ = Ãx + Bu + d, are I m à = J T J J w S[W ν ] 3 m, m 3 m 3 m m B = 3 m J J w W I m m, d = 3 J τ srp Θ=3 m, () where τ srp Θ=3 is the SRP torque when Θ = 3 and T results from numerically linearizing (8) about Θ = 3, i.e., τ srp T Θ + τ srp Θ=3. The discrete-time model is obtained from the continuous-time model () using Euler s forward method, where t denotes the sampling time and t Z, x t+ = Ax t + Bu t + d, () with A = I (6+m) (6+m) + à t, B = B t, and d = d t. DRIFT COUNTERACTION OPTIMAL CONTROL In this section, a class of DCOC problems for linear discrete-time systems is described. Based on this, we obtain a controller that maximizes the time until constraint violation occurs. The approach, which uses LP and MILP formulations, was recently developed by Zidek et al. 9. Problem Formulation Let the set of admissible control sequences {u t } be defined by where the time-varying control constraints are given by Furthermore, let the first exit-time from a time-varying set, U seq = {{u t } : u t U t for all t}, () U t = {u R m : C c,t u b c,t }. (3) G t = {x R n : C s,t x b s,t }, (4) given an initial state x and a control sequence {u t } U seq, be defined by τ(x, {u t }) = inf{t Z + : x t / G t x G }. () Then the DCOC problem of finding an admissible control sequence {u t } that maximizes the first exit-time is as follows max τ(x, {u t }) {u t} U seq subject to x t+ = A t x t + B t u t + d t, x G. In the following, we denote a solution to (6), if one exists, by {u t } and τ(x, {u t }) is the optimal exit-time. Note that {u t } may not be unique. (6) 4
5 Optimal Solution An optimal solution to problem (6), if one exists, can be obtained by solving the following MILP with time horizon N 9 min N δ t {x t},{u t},{δ τlb,...,δ N } t=τ lb subject to x t+ = A t x t + B t u t + d t, x G δ t δ t C s,t x t b s,t, t =,..., τ lb C s,t x t b s,t + Mδ t, t = τ lb,..., N (7) u t U t δ t {, } Z, t = τ lb,..., N, where M is a large positive constant, denotes the n-dimensional vector of ones, and the binary variable δ t is an indicator variable for the condition x t / G t. Thus, due to the constraint δ t δ t, minimizing the sum of δ t is equivalent to maximizing the first exit-time. The parameter τ lb Z + in (7) is a lower bound on the optimal exit-time, i.e., τ(x, {u t }) τ lb. Note that N τ(x, {u t }) is required when using the MILP (7) to solve problem (6). MILP is in the class of NP-complete problems and the worst-case computation time increases exponentially with the number of integer variables., However, we can relax the MILP and obtain good-quality suboptimal solutions by replacing the binary variables δ τlb,..., δ N with non-negative variables ε t R. 9 The resulting LP formulation in (8) may provide a better balance between computation time and performance of the solution than the MILP (7) and is therefore used in the MPC scheme outlined in Algorithm. The relaxed problem has the following form min N ε t {x t},{u t},{ε τlb,...,ε N } t=τ lb subject to x t+ = A t x t + B t u t + d t, x G ε t ε t C s,t x t b s,t, t =,..., τ lb C s,t x t b s,t + ε t, t = τ lb,..., N u t U t, (8) where τ lb Z + is a lower bound on the optimal exit-time as in (7). In order to obtain a solution that violates the constraints, the time horizon N needs to be greater than the optimal first exittime τ(x, {u t }). In addition, finding a good solution to the DCOC problem (6) requires N to be sufficiently close to τ(x, {u t }). Since τ(x, {u t }) is a priori unknown, a procedure was developed that iteratively adapts N until a proper solution is found. 9 This procedure is outlined in Algorithm. At Step of Algorithm, the lower bound τ lb is initialized using the zero-control solution τ(x, {,,..., }) (assuming U t ; otherwise, any admissible control sequence can be used here). Moreover, the time horizon N is set in Step by adding a constant N add Z + to τ lb. With
6 these parameters, LP (8) is solved. If no constraint violation occurs, N + is used as a new lower bound (Step 6) and the LP is solved for an increased N. This procedure is repeated until constraint violation occurs. A threshold N ub can be specified by the user to terminate the procedure (Step ) when a strict upper bound on the computation time needs to be satisfied (for example, for an MPC implementation). However, the solution may be suboptimal in case of N ub + N add < τ(x, {u t }). Algorithm Iterative procedure based on LP to solve (6) : τ lb min{τ(x, {,,..., }), N ub } : N τ lb + N add, N add Z + 3: {x t }, {u t }, {ε τlb,..., ε N } solution of LP (8) 4: τ max{t N : ε t = } + : if τ = N AND N < N ub then 6: τ lb τ 7: Go to Step 8: end if Strategy to Recover Constraint Satisfaction When constraint violation occurs, i.e., x t / G t for some t, in some cases it may be possible to recover and return to G t. In this regard, we modify the LP in (8) by removing the constraints ε t ε t. Hence, the LP in (9) is obtained, which tries to drive the state vector back into the set G t, see (4), due to the constraints C s,t x t b s,t + ε t and ε t. If x t G t is not possible, the LP solution brings x t as close as possible to G t, assuming that the time horizon N recover is sufficiently small. N recover min {x t},{u t},{ε t} subject to t= ε t x t+ = A t x t + B t u t + d t, x / G (9) ε t C s,t x t b s,t + ε t u t U t. MPC Implementation In order to compensate for unmodeled effects, a DCOC state feedback controller is obtained by employing MPC techniques and recomputing the control, using either LP (8) or (9), over a receding time horizon based on the current state. However, initial numerical results showed that, compared to the linear solution, this MPC strategy may use an excessive amount of control (RW accelerations) when applied to the nonlinear spacecraft model. Hence, in order to minimize the control effort, the absolute values of the control inputs are added to the respective objective functions 6
7 in LP (8) and (9), i.e., LP (8) min N ε t + N u t, t=τ lb LP (9) min Nrecover t= t= ε t + Nrecover t= u t, yielding nonlinear programs (as a consequence of the absolute values). The nonlinear programs are approximated with standard LP by introducing the variables γ t R m as substitutions for u t and adding the constraints γ t u t and γ t u t. Thus, LP (8) and (9), respectively, are extended for the MPC implementation as follows min {x t},{u t},{ε τlb,...,ε N },{γ t} subject to N t=τ lb ε t + w x t+ = A t x t + B t u t + d t, x G N t= T γ t ε t ε t C s,t x t b s,t, t =,..., τ lb C s,t x t b s,t + ε t, t = τ lb,..., N () u t U t γ t u t γ t u t, N recover N recover min ε t + w recover T γ t {x t},{u t},{ε t},{γ t} subject to t= t= x t+ = A t x t + B t u t + d t, x / G ε t () C s,t x t b s,t + ε t u t U t γ t u t γ t u t, where w > and w recover > are weights that we set to w = w recover =. in this paper. The resulting MPC scheme, emphasizing reduced control effort in addition to the DCOC objective, is outlined in Algorithm. At each time instant t sys, the control sequence is recomputed based on the current state. If the current state does not satisfy the constraints (x / G ), LP () is used for control computation to try to recover constraint satisfaction, see Step 8 of Algorithm. Otherwise, the control sequence is computed by the iterative scheme in Algorithm, where LP (8) is replaced by LP (), see Step. The first element, u, of the computed control sequence is applied to the system at Step and the procedure is repeated for the next sampling interval. Note that the moving time horizon approach requires shifting the time-varying constraints and dynamics for control computation. To address this, the original constraints and dynamics, initially stored in Step 3, are used in Steps and 6 to shift the constraints and dynamics based on the current time instant t sys. 7
8 Algorithm MPC scheme based on Algorithm : t sys : t G max{n recover, N ub + N add } 3: {G st,t }, {U st,t }, {A st,t }, {B st,t }, {d st,t } {G t }, {U t }, {A t }, {B t }, {d t } 4: x current state x(t sys ) : for t {,..., t G }: {G t }, {U t } {G st,t+tsys }, {U st,t+tsys } 6: for t {,..., t G }: {A t }, {B t }, {d t } {A st,t+tsys }, {B st,t+tsys }, {d st,t+tsys } 7: if x / G then 8: {u t } solution of () 9: else : {u t } output of Algorithm with LP () instead of LP (8) : end if : Apply u as input u(t sys ) to the system 3: t sys t sys + 4: Go to Step 4 SIMULATION RESULTS Parameter Units Value Parameter Units Value J kg m diag(43,, 3) J w kg m.43 L x m L y m. L z m l x m l y m. l z m ˆq sun I - [, /, / ] T Φ sun W/m 367 c m/sec 99,79,48 C diff -. Table : Model parameters. 6 Numerical results are presented for an example spacecraft with parameters listed in Table. Note that the vector representing the direction of the Sun is resolved in the inertial frame I in Table and needs to be transformed to the B frame continuously using the current attitude of the spacecraft. Throughout this section we assume minimum/maximum angular accelerations of ± rad/sec for each wheel, i.e., U t {u R m : α i [, ] rad/sec, i {,,..., m}}. The receding horizon drift counteraction control scheme in Algorithm is used to maximize the time until the spacecraft violates prescribed constraints on its attitude and RW spin rates. For the Euler angles, we consider constraints of the form φ min φ φ max, θ min θ θ max, ψ min ψ ψ max, () with φ min <, θ min <, ψ min <, φ max >, θ max >, and ψ max >. Similarly, for the spin rate of the ith wheel, ν i,min ν i ν i,max, where ν i,min, ν i,max R and i {,,..., m}. 8
9 The control strategy (Algorithm ) is based on the linear discrete-time dynamics in (). However, all closed-loop simulations are performed on the continuous-time nonlinear model given by () and (4), where the control input is constant in each sampling interval (zero-order hold). Due to the discrete-time formulation of the DCOC algorithm and unmodeled nonlinear effects, slight constraint violation may occur between the discrete time instants when applying the control inputs to the nonlinear continuous-time model. To address this, we tighten the attitude constraints in () and compute the control inputs for the modified constraints. Thus, if possible, we avoid constraint violation between discrete time instants as illustrated in Figure. Linear solution (discrete-time) Nonlinear solution with linear-based control (continuous-time) Constraint violation Linear solution (discrete-time) Tightened constraints Nonlinear solution with linear-based control (continuous-time) Figure : Illustration of effect of constraint tightening when linear-based control (discrete-time) is applied to the nonlinear continuous-time model. Top part: State trajectories without constraint tightening. Bottom part: State trajectories with constraint tightening. With (), the set defining the state constraints is given by G t (η) {x R 6+m : φ [( η)φ min, ( η)φ max ], θ [( η)θ min, ( η)θ max ], ψ [( η)ψ min, ( η)ψ max ], ν i [ν i,min, ν i,max ], i {,,..., m}}, (3) where η is a parameter to tighten the attitude constraints. Note that, for the nonlinear closedloop simulations, η =, and we terminate the simulation and report the corresponding exit-time once the current state vector is not inside G t () anymore. On the other hand, for the control computations according to Algorithm, we set η =. (see Figure ). In addition, the numerical conditioning of each LP is improved by normalizing the state vector x as follows x = Ox + o, (4) where O R 6+m 6+m and o R 6+m transform the state constraints given by G t (.) in (3) to G t = {x : φ, θ, ψ, ν i [, ], i {,,..., m}}. Moreover, ω j = 3 rad/s corresponds to ω j = and ω j = 3 rad/s corresponds to ω j =, where j {,, 3}. 9
10 8, (rad/sec )!! 3 A (rad)!? (rad) 3 (rad) The sampling time for the linear discrete-time model is chosen as t =. sec. The other parameters of the controller are set to N add =, N ub = 3, and N recover =. We use MATLAB a on a laptop with an i-63 processor and 8 GB RAM for computations. The Hybrid Toolbox 3 (lpsol function with default settings) is employed for LP solving. One RW # -3 LP: open-loop linear model (discrete-time) MPC: closed-loop nonlinear model (continuous-time) # # # # Figure : One operable RW: state and control trajectories vs. time. The case of one operable RW (m = ) with spin axis ḡ = [/ 3, / 3, / 3] T is considered. The attitude and spin rate constraints for this case are given by φ min = θ min =.7 rad, ψ min =.7 rad, φ max = θ max =.7 rad, ψ max =.7 rad, ν,min = rad/sec, and ν,max = rad/sec. Note that the lower bound on ν is chosen to avoid zero speed crossings and increase in wear and power consumption at low speeds. Example state and control trajectories are plotted with dotted black lines for an initial condition
11 of φ =. rad, θ =.3 rad, ψ =. rad, ω, = 4 rad/sec, ω, = rad/sec, ω 3, = 4 rad/sec, and ν, = rad/sec in Figure. The state and control constraints are indicated by dashed red lines in Figure. In addition, Figure shows the open-loop solution of LP (8) for the discrete-time linear model with η =, i.e., for G t (), which violates the constraints after 4 discrete time steps corresponding to 4 t = 98.4 sec =.64 min. For comparison, this value is identical to the open-loop solution of MILP (7), which was solved with the milpsol function of the Hybrid Toolbox. 3 It can be seen that the trajectories of the closedloop nonlinear model with the MPC implementation are close to the linear open-loop solution. Constraint violation occurs after 94.4 sec (.7 min), which is as close as 4 % to the optimal exittime of the linear problem. The worst-case computation time for recomputing the control during the closed-loop simulations was.4 sec < t = sec and the average computation time was. sec. Two RWs Now two operable RWs are assumed (m = ), increasing the control authority of the spacecraft. The spin axes of the two wheels are given by ḡ = [,, ] T and ḡ = [/ 3, / 3, / 3] T. The constraints on spacecraft attitude and RW spin rates for this case are set to φ min = θ min =.7 rad, ψ min =.7 rad, φ max = θ max =.7 rad, ψ max =.7 rad, ν,min = rad/sec, ν,max = rad/sec, ν,min = rad/sec, and ν,max = rad/sec. Figure 3 shows example state and control trajectories over time for an initial condition of φ =.7 rad, θ =.9 rad, ψ =.4 rad, ω, =. rad/sec, ω, = rad/sec, ω 3, =, ν, = 4 rad/sec, ν, = 8 rad/sec. In addition to the closed-loop nonlinear continuous-time solution (dotted black lines), the open-loop solution to the linear problem (for η =, i.e., G t ()) is plotted with solid black lines in Figure 3. The optimal exit-time for the linear discrete-time problem is 7 sec = 4. min (8 steps) as obtained by both MILP (7) and LP (8). The MPC implementation generates nonlinear trajectories that are close to the linear solution. Constraint violation occurs after 63.4 sec (4.39 min), which is. % lower than the optimal exit-time of the linear problem. The longest computation times for recomputing the control occur in the beginning when the time until constraint violation is long. In this case, a worst-case computation time of.44 sec < t =. sec was recorded while the average computation time over all sampling instants was.4 sec. Three RWs Finally, the case of three RWs (m = 3) is treated, where ḡ = [,, ] T, ḡ = [,, ] T, and ḡ 3 = [,, ] T. The constraints on the attitude and RW spin rates are φ min = θ min = ψ min =.7 rad, φ max = θ max = ψ max =.7 rad, ν,min = ν,min = ν 3,min = rad/sec, and ν,max = ν,max = ν 3,max = rad/sec. The initial values for the spacecraft attitude and angular velocity are φ =.7 rad, θ = ψ =, ω, =.8 rad/sec, ω, =.9 rad/sec, and ω 3, =.3 rad/sec. For the initial RW spin rates, two cases are considered: ν = [4, 3, 8] T rad/sec and ν = [4, 3, 8] T rad/sec. For the first case, the spin rates are close to the lower boundary which reduces the control authority. Thus, constraint violation occurs in finite time. This can be seen in Figure 4 showing the trajectories of the closed-loop nonlinear continuous-time model with the
12 MPC implementation (dotted black lines) as well as the open-loop solution with linear discretetime dynamics and η = (solid black lines). The optimal exit-time for the linear problem is 34 sec =.67 min (36 time steps), whereas constraint violation occurs with the nonlinear model after 37 sec =.4 min, which is 4 % less than the optimal exit-time of the linear problem. The worst-case computation time was. sec < t =. sec with an average computation time of.6 sec. For the second case, the initial RW spin rates are further away from the boundaries. The increased control authority may enable to satisfy constraints indefinitely. Simulation results for the first sec are plotted in Figure showing that the proposed feedback controller (dotted black lines) can keep the spacecraft within the limits. A worst-case computation time of.66 sec < t =. sec was recorded while the average computation time over all sampling instants was.79 sec. Note that the open-loop solution of the linear problem in Figure (solid black line) was obtained for a time horizon of N = 4 in (8), not considering the system s evolution for t > 4. CONCLUSION This paper considered an optimal control problem of maximizing the time until prescribed constraints are violated with applications to spacecraft attitude control. The proposed state feedback controller is based on linear programming. In addition, a similar controller was developed to recover constraint satisfaction in case constraints are violated. Closed-loop simulations on the nonlinear model showed that drift from solar radiation pressure is successfully counteracted in order to satisfy constraints on the spacecraft attitude and reaction wheel spin rates for as long as possible. ACKNOWLEDGMENT This research is supported by the National Science Foundation Award Number EECS REFERENCES [] P. Romero, J. M. Gambi, and E. Patiño, Stationkeeping manoeuvres for geostationary satellites using feedback control techniques, Aerospace Science and Technology, Vol., No., 7, pp [] B.-S. Lee, Y. Hwang, H.-Y. Kim, and S. Park, East West station-keeping maneuver strategy for COMS satellite using iterative process, Advances in space research, Vol. 47, No.,, pp [3] R. A. E. Zidek and I. V. Kolmanovsky, Deterministic drift counteraction optimal control and its application to satellite life extension, 4th IEEE Conference on Decision and Control (CDC), IEEE,, pp [4] R. A. E. Zidek and I. V. Kolmanovsky, Geostationary satellite station keeping using drift counteraction optimal control, 6th AAS/AIAA Space flight mechanics meeting, 6, pp. AAS 6 7. [] P. Tsiotras and J. Luo, Control of underactuated spacecraft with bounded inputs, Automatica, Vol. 36, No. 8,, pp [6] C. D. Petersen, F. Leve, M. Flynn, and I. Kolmanovsky, Recovering linear controllability of an underactuated spacecraft by exploiting solar radiation pressure, Journal of Guidance, Control, and Dynamics, Vol. 39, No. 4,, pp [7] H. Kuninaka and J. Kawaguchi, Lessons learned from round trip of Hayabusa asteroid explorer in deep space, Aerospace Conference, IEEE, IEEE,, pp. 8. [8] R. Cowen, The wheels come off Kepler, Nature, Vol. 497, 3, pp [9] R. A. E. Zidek, A. Bemporad, and I. V. Kolmanovsky, Optimal and receding horizon drift counteraction control: linear programming approaches, American Control Conference (ACC) 7, accepted. [] R. Linares, M. K. Jah, J. L. Crassidis, F. A. Leve, and T. Kelecy, Astrometric and photometric data fusion for inactive space object mass and area estimation, Acta Astronautica, Vol. 99, 4, pp.. [] A. Bemporad and M. Morari, Control of systems integrating logic, dynamics, and constraints, Automatica, Vol. 3, No. 3, 999, pp
13 , (rad/sec ), (rad/sec ) 8 8!! 3 A (rad)!? (rad) 3 (rad) [] A. Richards and J. How, Mixed-integer programming for control, Proceedings of the, American Control Conference,., IEEE,, pp [3] A. Bemporad, Hybrid Toolbox - User s Guide, 4. bemporad/hybrid/toolbox. # -3 LP: open-loop linear model (discrete-time) MPC: closed-loop nonlinear model (continuous-time) # # # # Figure 3: Two operable RWs: state and control trajectories vs. time. 3
14 , (rad/sec ), 3 (rad/sec ) 8 3, (rad/sec ) 8 8!! 3 A (rad)!? (rad) 3 (rad).. LP: open-loop linear model (discrete-time) MPC: closed-loop nonlinear model (continuous-time) # # # Figure 4: Three operable RWs with ν = [4, 3, 8] T rad/sec: state and control trajectories vs. time. 4
15 .. 3 (rad).? (rad) !. A (rad) # LP: open-loop linear model (discrete-time) MPC: closed-loop nonlinear model (continuous-time)! #-4 - #-4 4! , ,3, Figure : Three operable RWs with ν = [4, 3, 8]T rad/sec, where the optimal exit-time may be indefinite: state and control trajectories vs. time.
Drift Counteraction Optimal Control: Theory and Applications to Autonomous Cars and Spacecraft
Drift Counteraction Optimal Control: Theory and Applications to Autonomous Cars and Spacecraft by Robert A. E. Zidek A dissertation submitted in partial fulfillment of the requirements for the degree of
More informationQuaternion-Based Tracking Control Law Design For Tracking Mode
A. M. Elbeltagy Egyptian Armed forces Conference on small satellites. 2016 Logan, Utah, USA Paper objectives Introduction Presentation Agenda Spacecraft combined nonlinear model Proposed RW nonlinear attitude
More informationIAA-CU A Simulator for Robust Attitude Control of Cubesat Deploying Satellites
A Simulator for Robust Attitude Control of Cubesat Deploying Satellites Giovanni Mattei, George Georgiou, Angelo Pignatelli, Salvatore Monaco Abstract The paper deals with the development and testing of
More informationMixed Control Moment Gyro and Momentum Wheel Attitude Control Strategies
AAS03-558 Mixed Control Moment Gyro and Momentum Wheel Attitude Control Strategies C. Eugene Skelton II and Christopher D. Hall Department of Aerospace & Ocean Engineering Virginia Polytechnic Institute
More informationMAE 142 Homework #5 Due Friday, March 13, 2009
MAE 142 Homework #5 Due Friday, March 13, 2009 Please read through the entire homework set before beginning. Also, please label clearly your answers and summarize your findings as concisely as possible.
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 informationSpacecraft Attitude Control with RWs via LPV Control Theory: Comparison of Two Different Methods in One Framework
Trans. JSASS Aerospace Tech. Japan Vol. 4, No. ists3, pp. Pd_5-Pd_, 6 Spacecraft Attitude Control with RWs via LPV Control Theory: Comparison of Two Different Methods in One Framework y Takahiro SASAKI,),
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 informationSliding Mode Control Strategies for Spacecraft Rendezvous Maneuvers
Osaka University March 15, 2018 Sliding Mode Control Strategies for Spacecraft Rendezvous Maneuvers Elisabetta Punta CNR-IEIIT, Italy Problem Statement First Case Spacecraft Model Position Dynamics Attitude
More informationStabilization of Angular Velocity of Asymmetrical Rigid Body. Using Two Constant Torques
Stabilization of Angular Velocity of Asymmetrical Rigid Body Using Two Constant Torques Hirohisa Kojima Associate Professor Department of Aerospace Engineering Tokyo Metropolitan University 6-6, Asahigaoka,
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 informationStation keeping and momentum management of low-thrust satellites using MPC
MITSUBISHI ELECTRIC RESEARCH LABORATORIES http://www.merl.com Station keeping and momentum management of low-thrust satellites using MPC Weiss, A.; Kalabic, U.; Di Cairano, S. TR2018-048 July 10, 2018
More informationTheory in Model Predictive Control :" Constraint Satisfaction and Stability!
Theory in Model Predictive Control :" Constraint Satisfaction and Stability Colin Jones, Melanie Zeilinger Automatic Control Laboratory, EPFL Example: Cessna Citation Aircraft Linearized continuous-time
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 informationSpacecraft Attitude Control using CMGs: Singularities and Global Controllability
1 / 28 Spacecraft Attitude Control using CMGs: Singularities and Global Controllability Sanjay Bhat TCS Innovation Labs Hyderabad International Workshop on Perspectives in Dynamical Systems and Control
More informationRobust Adaptive Attitude Control of a Spacecraft
Robust Adaptive Attitude Control of a Spacecraft AER1503 Spacecraft Dynamics and Controls II April 24, 2015 Christopher Au Agenda Introduction Model Formulation Controller Designs Simulation Results 2
More informationFeedback Control of Spacecraft Rendezvous Maneuvers using Differential Drag
Feedback Control of Spacecraft Rendezvous Maneuvers using Differential Drag D. Pérez 1 and R. Bevilacqua Rensselaer Polytechnic Institute, Troy, New York, 1180 This work presents a feedback control strategy
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 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 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 informationAttitude Control Strategy for HAUSAT-2 with Pitch Bias Momentum System
SSC06-VII-5 Attitude Control Strategy for HAUSAT-2 with Pitch Bias Momentum System Young-Keun Chang, Seok-Jin Kang, Byung-Hoon Lee, Jung-on Choi, Mi-Yeon Yun and Byoung-Young Moon School of Aerospace and
More informationRobust Control of a 3D Space Robot with an Initial Angular Momentum based on the Nonlinear Model Predictive Control Method
Vol. 9, No. 6, 8 Robust Control of a 3D Space Robot with an Initial Angular Momentum based on the Nonlinear Model Predictive Control Method Tatsuya Kai Department of Applied Electronics Faculty of Industrial
More informationEXAMPLE: MODELING THE PT326 PROCESS TRAINER
CHAPTER 1 By Radu Muresan University of Guelph Page 1 EXAMPLE: MODELING THE PT326 PROCESS TRAINER The PT326 apparatus models common industrial situations in which temperature control is required in the
More informationSpacecraft Rate Damping with Predictive Control Using Magnetic Actuators Only
C. Böhm a M. Merk a W. Fichter b F. Allgöwer a Spacecraft Rate Damping with Predictive Control Using Magnetic Actuators Only Stuttgart, March 2009 a Institute of Systems Theory and Automatic Control, University
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 informationOn the Inherent Robustness of Suboptimal Model Predictive Control
On the Inherent Robustness of Suboptimal Model Predictive Control James B. Rawlings, Gabriele Pannocchia, Stephen J. Wright, and Cuyler N. Bates Department of Chemical & Biological Engineering Computer
More informationVisual Feedback Attitude Control of a Bias Momentum Micro Satellite using Two Wheels
Visual Feedback Attitude Control of a Bias Momentum Micro Satellite using Two Wheels Fuyuto Terui a, Nobutada Sako b, Keisuke Yoshihara c, Toru Yamamoto c, Shinichi Nakasuka b a National Aerospace Laboratory
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 informationATTITUDE CONTROL MECHANIZATION TO DE-ORBIT SATELLITES USING SOLAR SAILS
IAA-AAS-DyCoSS2-14-07-02 ATTITUDE CONTROL MECHANIZATION TO DE-ORBIT SATELLITES USING SOLAR SAILS Ozan Tekinalp, * Omer Atas INTRODUCTION Utilization of solar sails for the de-orbiting of satellites is
More informationEast West GEO Satellite Station-Keeping with Degraded Thruster Response
Aerospace 2015, 2, 581-601; doi:10.3390/aerospace2040581 OPEN ACCESS aerospace ISSN 2226-4310 www.mdpi.com/journal/aerospace Article East West GEO Satellite Station-Keeping with Degraded Thruster Response
More informationWEIGHTING MATRICES DETERMINATION USING POLE PLACEMENT FOR TRACKING MANEUVERS
U.P.B. Sci. Bull., Series D, Vol. 75, Iss. 2, 2013 ISSN 1454-2358 WEIGHTING MATRICES DETERMINATION USING POLE PLACEMENT FOR TRACKING MANEUVERS Raluca M. STEFANESCU 1, Claudiu L. PRIOROC 2, Adrian M. STOICA
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 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 informationHandling Roll Constraints for Path Following of Marine Surface Vessels using Coordinated Rudder and Propulsion Control
2010 American Control Conference Marriott Waterfront, Baltimore, MD, USA June 30-July 02, 2010 FrB15.5 Handling Roll Constraints for Path Following of Marine Surface Vessels using Coordinated Rudder and
More informationFIBER OPTIC GYRO-BASED ATTITUDE DETERMINATION FOR HIGH- PERFORMANCE TARGET TRACKING
FIBER OPTIC GYRO-BASED ATTITUDE DETERMINATION FOR HIGH- PERFORMANCE TARGET TRACKING Elias F. Solorzano University of Toronto (Space Flight Laboratory) Toronto, ON (Canada) August 10 th, 2016 30 th AIAA/USU
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 informationCHAPTER 4 CONVENTIONAL CONTROL FOR SATELLITE ATTITUDE
93 CHAPTER 4 CONVENTIONAL CONTROL FOR SATELLITE ATTITUDE 4.1 INTRODUCTION Attitude control is the process of achieving and maintaining an orientation in space. The orientation control of a rigid body has
More informationAttitude Control of a Bias Momentum Satellite Using Moment of Inertia
I. INTRODUCTION Attitude Control of a Bias Momentum Satellite Using Moment of Inertia HYOCHOONG BANG Korea Advanced Institute of Science and Technology HYUNG DON CHOI Korea Aerospace Research Institute
More informationFeedback Optimal Control of Low-thrust Orbit Transfer in Central Gravity Field
Vol. 4, No. 4, 23 Feedback Optimal Control of Low-thrust Orbit Transfer in Central Gravity Field Ashraf H. Owis Department of Astronomy, Space and Meteorology, Faculty of Science, Cairo University Department
More informationResearch Article Shorter Path Design and Control for an Underactuated Satellite
Hindawi International Journal of Aerospace Engineering Volume 27 Article ID 8536732 9 pages https://doi.org/.55/27/8536732 Research Article Shorter Path Design and Control for an Underactuated Satellite
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 informationAdaptive Dynamic Inversion Control of a Linear Scalar Plant with Constrained Control Inputs
5 American Control Conference June 8-, 5. Portland, OR, USA ThA. Adaptive Dynamic Inversion Control of a Linear Scalar Plant with Constrained Control Inputs Monish D. Tandale and John Valasek Abstract
More informationPRELIMINARY HARDWARE DESIGN OF ATTITUDE CONTROL SUBSYSTEM OF LEONIDAS SPACECRAFT
PRELIMINARY HARDWARE DESIGN OF ATTITUDE CONTROL SUBSYSTEM OF LEONIDAS SPACECRAFT Chak Shing Jackie Chan College of Engineering University of Hawai i at Mānoa Honolulu, HI 96822 ABSTRACT In order to monitor
More informationLesson Rigid Body Dynamics
Lesson 8 Rigid Body Dynamics Lesson 8 Outline Problem definition and motivations Dynamics of rigid bodies The equation of unconstrained motion (ODE) User and time control Demos / tools / libs Rigid Body
More informationTracking Rigid Body Motion Using Thrusters and Momentum. Wheels
JAS 199 Tracking Rigid Body Motion Using Thrusters and Momentum Wheels Christopher D. Hall, Panagiotis Tsiotras, and Haijun Shen Abstract Tracking control laws are developed for a rigid spacecraft using
More informationSATELLITE ATTITUDE CONTROL SYSTEM DESIGN WITH NONLINEAR DYNAMICS AND KINEMTICS OF QUATERNION USING REACTION WHEELS
SATELLITE ATTITUDE CONTROL SYSTEM DESIGN WITH NONLINEAR DYNAMICS AND KINEMTICS OF QUATERNION USING REACTION WHEELS Breno Braga Galvao Maria Cristina Mendes Faustino Luiz Carlos Gadelha de Souza breno.braga.galvao@gmail.com
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 informationA Globally Stabilizing Receding Horizon Controller for Neutrally Stable Linear Systems with Input Constraints 1
A Globally Stabilizing Receding Horizon Controller for Neutrally Stable Linear Systems with Input Constraints 1 Ali Jadbabaie, Claudio De Persis, and Tae-Woong Yoon 2 Department of Electrical Engineering
More informationEVALUATING ORBITS WITH POTENTIAL TO USE SOLAR SAIL FOR STATION-KEEPING MANEUVERS
IAA-AAS-DyCoSS2-14-16-03 EVALUATING ORBITS WITH POTENTIAL TO USE SOLAR SAIL FOR STATION-KEEPING MANEUVERS Thais C. Oliveira, * and Antonio F. B. A. Prado INTRODUCTION The purpose of this paper is to find
More informationDesign of Sliding Mode Attitude Control for Communication Spacecraft
Design of Sliding Mode Attitude Control for Communication Spacecraft Erkan Abdulhamitbilal 1 and Elbrous M. Jafarov 1 ISTAVIA Engineering, Istanbul Aeronautics and Astronautics Engineering, Istanbul Technical
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 informationSimplified Filtering Estimator for Spacecraft Attitude Determination from Phase Information of GPS Signals
WCE 7, July - 4, 7, London, U.K. Simplified Filtering Estimator for Spacecraft Attitude Determination from Phase Information of GPS Signals S. Purivigraipong, Y. Hashida, and M. Unwin Abstract his paper
More informationTrajectory-tracking control of a planar 3-RRR parallel manipulator
Trajectory-tracking control of a planar 3-RRR parallel manipulator Chaman Nasa and Sandipan Bandyopadhyay Department of Engineering Design Indian Institute of Technology Madras Chennai, India Abstract
More informationREACTION WHEEL CONFIGURATIONS FOR HIGH AND MIDDLE INCLINATION ORBITS
REACTION WHEEL CONFIGURATIONS FOR HIGH AND MIDDLE INCLINATION ORBITS Zuliana Ismail and Renuganth Varatharajoo Department of Aerospace Engineering, Universiti Putra Malaysia, Malaysia E-Mail: zuliana.ismail@gmail.com
More informationSatellite Attitude Control System Design Using Reaction Wheels Bhanu Gouda Brian Fast Dan Simon
Satellite Attitude Control System Design Using Reaction Wheels Bhanu Gouda Brian Fast Dan Simon Outline 1. Overview of Attitude Determination and Control system. Problem formulation 3. Control schemes
More informationSemi-Analytical Guidance Algorithm for Fast Retargeting Maneuvers Computation during Planetary Descent and Landing
ASTRA 2013 - ESA/ESTEC, Noordwijk, the Netherlands Semi-Analytical Guidance Algorithm for Fast Retargeting Maneuvers Computation during Planetary Descent and Landing Michèle LAVAGNA, Paolo LUNGHI Politecnico
More informationStation-keeping and momentum-management on halo orbits around L2: Linear-quadratic feedback and model predictive control approaches
MITSUBISHI ELECTRIC RESEARCH LABORATORIES http://www.merl.com Station-keeping and momentum-management on halo orbits around L2: Linear-quadratic feedback and model predictive control approaches Kalabic,
More informationsc Control Systems Design Q.1, Sem.1, Ac. Yr. 2010/11
sc46 - Control Systems Design Q Sem Ac Yr / Mock Exam originally given November 5 9 Notes: Please be reminded that only an A4 paper with formulas may be used during the exam no other material is to be
More informationMulticopter Design Optimization and Validation
Modeling, Identification and Control, Vol 36, No 2, 2015, pp 67 79, ISSN 1890 1328 Multicopter Design Optimization and Validation Øyvind Magnussen Morten Ottestad Geir Hovland University of Agder, Faculty
More informationInertia-Free Spacecraft Attitude Control Using. Reaction Wheels
Inertia-Free Spacecraft Attitude Control Using Reaction Wheels Avishai Weiss, Ilya Kolmanovsky, Dennis S. Bernstein, Department of Aerospace Engineering, The University of Michigan, Ann Arbor, MI 489-4
More informationParametric Integrated Perturbation Analysis - Sequential Quadratic Programming Approach for Minimum-time Model Predictive Control
Preprints of the 9th World Congress The International Federation of Automatic Control Cape Town, South Africa. August 4-9, 4 Parametric Integrated Perturbation Analysis - Sequential Quadratic Programming
More informationMathematical Models for the Triaxial Attitude Control Testbed
Draft Article (For Engineering Journals) Mathematical Models for the Triaxial Attitude Control Testbed Sangbum Cho, Jinglai Shen, N. Harris McClamroch Department of Aerospace Engineering University of
More informationGlobal Trajectory Design for Position and Attitude Control of an Underactuated Satellite *
Trans. Japan Soc. Aero. Space Sci. Vol. 59, No. 3, pp. 7 4, 26 Global Trajectory Design for Position and Attitude Control of an Underactuated Satellite * Yasuhiro YOSHIMURA, ) Takashi MATSUNO, 2) and Shinji
More informationCHAPTER 5 FUZZY LOGIC FOR ATTITUDE CONTROL
104 CHAPTER 5 FUZZY LOGIC FOR ATTITUDE CONTROL 5.1 INTRODUCTION Fuzzy control is one of the most active areas of research in the application of fuzzy set theory, especially in complex control tasks, which
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 informationCOUPLED ORBITAL AND ATTITUDE CONTROL SIMULATION
COUPLED ORBITAL AND ATTITUDE CONTROL SIMULATION Scott E. Lennox AOE 5984: Advanced Attitude Spacecraft Dynamics and Control December 12, 23 INTRODUCTION In the last few years the space industry has started
More informationAN ANALYTICAL SOLUTION TO QUICK-RESPONSE COLLISION AVOIDANCE MANEUVERS IN LOW EARTH ORBIT
AAS 16-366 AN ANALYTICAL SOLUTION TO QUICK-RESPONSE COLLISION AVOIDANCE MANEUVERS IN LOW EARTH ORBIT Jason A. Reiter * and David B. Spencer INTRODUCTION Collision avoidance maneuvers to prevent orbital
More informationPartial Attitude Consensus for Underactuated Satellite Clusters
6 IEEE 55th Conference on Decision and Control (CDC) ARIA Resort & Casino December -4, 6, Las Vegas, USA Partial Attitude Consensus for Underactuated Satellite Clusters John Matthew Brewer and Panagiotis
More informationENAE483: Principles of Space System Design Power Propulsion Thermal System
Power Propulsion Thermal System Team B4: Ben Abresch Jason Burr Kevin Lee Scott Wingate November 8th, 2012 Presentation Overview Mission Guidelines Project Specifications Initial Design Power Thermal Insulation
More informationSERVO TRACKING OF TARGETS AT SEA
Int. J. Appl. Math. Comput. Sci., 2006, Vol. 16, No. 2, 197 207 SERVO TRACKING OF TARGETS AT SEA ANIS AHMED,,MIECZYSŁAW A. BRDYŚ Department of Electronic, Electrical and Computer Engineering School of
More informationSail-Assisted End-of-Life Disposal of High-LEO Satellites
4th International Symposium on Solar Sailing Kyoto, Japan, January 17-20, 2017 Sail-Assisted End-of-Life Disposal of High-LEO Satellites Sergey Trofimov Keldysh Institute of Applied Mathematics Russian
More informationROBUST THREE-AXIS ATTITUDE STABILIZATION FOR INERTIAL POINTING SPACECRAFT USING MAGNETORQUERS *
IAA-AAS-DyCoSS2-4-4-7 AAS 4-52 ROBUST THREE-AXIS ATTITUDE STABILIZATION FOR INERTIAL POINTING SPACECRAFT USING MAGNETORQUERS * INTRODUCTION Fabio Celani In this work feedback control laws are designed
More informationThe ϵ-capacity of a gain matrix and tolerable disturbances: Discrete-time perturbed linear systems
IOSR Journal of Mathematics (IOSR-JM) e-issn: 2278-5728, p-issn: 2319-765X. Volume 11, Issue 3 Ver. IV (May - Jun. 2015), PP 52-62 www.iosrjournals.org The ϵ-capacity of a gain matrix and tolerable disturbances:
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 informationAn Inverse Dynamics Attitude Control System with Autonomous Calibration. Sanny Omar Dr. David Beale Dr. JM Wersinger
An Inverse Dynamics Attitude Control System with Autonomous Calibration Sanny Omar Dr. David Beale Dr. JM Wersinger Outline Attitude Determination and Control Systems (ADACS) Overview Coordinate Frames
More informationAttitude Tracking Control of a Small Satellite in Low Earth Orbit
AIAA Guidance, Navigation, and Control Conference 10-13 August 2009, Chicago, Illinois AIAA 2009-5902 Attitude Tracking Control of a Small Satellite in Low Earth Orbit Amit K. Sanyal Zachary Lee-Ho In
More informationPath Following for Marine Surface Vessels with Rudder and Roll Constraints: an MPC Approach
2009 American Control Conference Hyatt Regency Riverfront, St. Louis, MO, USA June 10-12, 2009 ThC11.6 Path Following for Marine Surface Vessels with Rudder and Roll Constraints: an MPC Approach Zhen Li,
More informationCase Study: The Pelican Prototype Robot
5 Case Study: The Pelican Prototype Robot The purpose of this chapter is twofold: first, to present in detail the model of the experimental robot arm of the Robotics lab. from the CICESE Research Center,
More informationOptElec: an Optimisation Software for Low-Thrust Orbit Transfer Including Satellite and Operation Constraints
OptElec: an Optimisation Software for Low-Thrust Orbit Transfer Including Satellite and Operation Constraints 7th International Conference on Astrodynamics Tools and Techniques, DLR, Oberpfaffenhofen Nov
More informationDesign of Attitude Determination and Control Subsystem
Design of Attitude Determination and Control Subsystem 1) Control Modes and Requirements Control Modes: Control Modes Explanation 1 ) Spin-Up Mode - Acquisition of Stability through spin-up maneuver -
More informationNonlinear and Cooperative Control of Multiple Hovercraft with Input Constraints
Submitted as a Regular Paper to the European Control Conference Revised: December 8, Nonlinear and Cooperative Control of Multiple Hovercraft with Input Constraints William B. Dunbar, Reza Olfati-Saber,
More informationOn Sun-Synchronous Orbits and Associated Constellations
On Sun-Synchronous Orbits and Associated Constellations Daniele Mortari, Matthew P. Wilkins, and Christian Bruccoleri Department of Aerospace Engineering, Texas A&M University, College Station, TX 77843,
More informationSpacecraft Bus / Platform
Spacecraft Bus / Platform Propulsion Thrusters ADCS: Attitude Determination and Control Subsystem Shield CDH: Command and Data Handling Subsystem Payload Communication Thermal Power Structure and Mechanisms
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 informationKINEMATIC EQUATIONS OF NONNOMINAL EULER AXIS/ANGLE ROTATION
IAA-AAS-DyCoSS -14-10-01 KINEMATIC EQUATIONS OF NONNOMINAL EULER AXIS/ANGLE ROTATION Emanuele L. de Angelis and Fabrizio Giulietti INTRODUCTION Euler axis/angle is a useful representation in many attitude
More informationSatellite Attitude Control by Quaternion-Based Backstepping
25 American Control Conference June 8-1, 25. Portland, OR, USA WeB11.4 Satellite Attitude Control by Quaternion-Based Backstepping Raymond Kristiansen* and Per J. Nicklasson** Department of Computer Science,
More informationOPTIMAL SPACECRAF1 ROTATIONAL MANEUVERS
STUDIES IN ASTRONAUTICS 3 OPTIMAL SPACECRAF1 ROTATIONAL MANEUVERS JOHNL.JUNKINS Texas A&M University, College Station, Texas, U.S.A. and JAMES D.TURNER Cambridge Research, Division of PRA, Inc., Cambridge,
More informationEXPERIMENTAL COMPARISON OF TRAJECTORY TRACKERS FOR A CAR WITH TRAILERS
1996 IFAC World Congress San Francisco, July 1996 EXPERIMENTAL COMPARISON OF TRAJECTORY TRACKERS FOR A CAR WITH TRAILERS Francesco Bullo Richard M. Murray Control and Dynamical Systems, California Institute
More informationDesign and Control of Novel Tri-rotor UAV
UKACC International Conference on Control Cardiff, UK, -5 September Design and Control of Novel Tri-rotor UAV Mohamed Kara Mohamed School of Electrical and Electronic Engineering The University of Manchester
More informationEE C128 / ME C134 Feedback Control Systems
EE C128 / ME C134 Feedback Control Systems Lecture Additional Material Introduction to Model Predictive Control Maximilian Balandat Department of Electrical Engineering & Computer Science University of
More informationAttitude Control Method of Spin-Stabilized Satellite Base on Equiangular Precession
International Conference on Artificial Intelligence and Engineering Applications (AIEA 2016) Attitude Control Method of Spin-Stabilized Satellite Base on Equiangular Precession Yonggang Li a, Jiping Ren
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 informationPredictive control of hybrid systems: Input-to-state stability results for sub-optimal solutions
Predictive control of hybrid systems: Input-to-state stability results for sub-optimal solutions M. Lazar, W.P.M.H. Heemels a a Eindhoven Univ. of Technology, P.O. Box 513, 5600 MB Eindhoven, The Netherlands
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 informationRobust Model Predictive Control for Autonomous Vehicle/Self-Driving Cars
Robust Model Predictive Control for Autonomous Vehicle/Self-Driving Cars Che Kun Law, Darshit Dalal, Stephen Shearrow A robust Model Predictive Control (MPC) approach for controlling front steering of
More informationGravity Probe B Data Analysis Challenges, Insights, and Results
Gravity Probe B Data Analysis Challenges, Insights, and Results Mac Keiser April 15, 27 April 15, 27 Jacksonville, FL 1 Topics Data Analysis Strategy and Methods Challenges and Insights Trapped Magnetic
More informationQUATERNION FEEDBACK ATTITUDE CONTROL DESIGN: A NONLINEAR H APPROACH
Asian Journal of Control, Vol. 5, No. 3, pp. 406-4, September 003 406 Brief Paper QUAERNION FEEDBACK AIUDE CONROL DESIGN: A NONLINEAR H APPROACH Long-Life Show, Jyh-Ching Juang, Ying-Wen Jan, and Chen-zung
More informationDISCRETE-TIME TIME-VARYING ROBUST STABILIZATION FOR SYSTEMS IN POWER FORM. Dina Shona Laila and Alessandro Astolfi
DISCRETE-TIME TIME-VARYING ROBUST STABILIZATION FOR SYSTEMS IN POWER FORM Dina Shona Laila and Alessandro Astolfi Electrical and Electronic Engineering Department Imperial College, Exhibition Road, London
More informationEfficient robust optimization for robust control with constraints Paul Goulart, Eric Kerrigan and Danny Ralph
Efficient robust optimization for robust control with constraints p. 1 Efficient robust optimization for robust control with constraints Paul Goulart, Eric Kerrigan and Danny Ralph Efficient robust optimization
More informationLearning Model Predictive Control for Iterative Tasks: A Computationally Efficient Approach for Linear System
Learning Model Predictive Control for Iterative Tasks: A Computationally Efficient Approach for Linear System Ugo Rosolia Francesco Borrelli University of California at Berkeley, Berkeley, CA 94701, USA
More information