Moving Mass Actuated Missile Control Using Convex Optimization Techniques
|
|
- Elaine Andrews
- 5 years ago
- Views:
Transcription
1 Moving Mass Actuated Missile Control Using Convex Optimization echniques Sesha Sai Vaddi Optimal Synthesis, 868 San-Antonio Road, Palo Alto, CA, 9433, USA his paper deals with the control of a moving mass actuated missile in the presence of state and control constraints Quadratic programming which is a convex optimization formulation is employed for control computation Planar equations of motion of a moving mass actuated missile are obtained in closed form using Lagrangian approach hese equations are analytically linearized and used for control purposes Control is computed at each instant of time using the time-to-go information Results obtained from closed loop nonlinear simulations demonstrate the effectiveness of the controller Nomenclature I Moment of inertia of the missile M b Mass of the missile m Moving mass Position of the moving mass along its axis Kinetic energy F Axial rocket motor thrust g Acceleration due to gravity X North position of the missile in an inertial NED frame Z Altitude of the missile θ Attitude angle of the missile M t otal mass M Mass matrix u Actuator force X State vector of the dynamic system Z State vector of the optimization problem Subscript b Body t otal c Continuous d discrete I Introduction Moving mass actuation system is a novel concept introduced recently for the control of missiles hese actuators are completely enclosed within the geometric envelope of the vehicle and are equally effective inside and outside the atmosphere he objective of this research is to design a control system for a moving mass actuated missile he dynamic model of the system is nonlinear and has both actuator and state constraints Control of nonlinear dynamical systems is more difficult than the control of their linear counterparts he presence of Research Scientist, Optimal Synthesis, 868 San-Antonio Road, Palo Alto, CA, of 13
2 state and actuator constraints adds up to the complexity of the problem No general framework exists for control design of nonlinear systems with state and actuator constraints Current approach is based on obtaining a convex relaxation Current controller design approach is based on quadratic programming A linear discretised version of the nonlinear system is used for this purpose he approach is particularly attractive for its ability to explicitly accommodate the before mentioned constraints A Background Moving mass actuation system was recently discussed by Menon 1 etal An integrated guidance-control system for target interception was designed in that paper Numerical feedback linearization technique which is a feature of commercially distributed software Nonlinear Synthesis oolbox is used for the controller design State and actuator constraints arising from limited mass displacement and peak actuator force are not included in the design process he controller is implemented on a simulation where the states and controls are saturated when these constraints are violated Similar setting has been used for testing a controller in a recent publication 2 that is designed using a combination of feedback linearization and finiteinterval optimal control Even though the controller works well in a wide range of scenarios it cannot prevent saturation he current research is primarily motivated by the need to avoid saturation and thereby prevent any undesirable consequences associated with it A controller formulation that can accommodate the state and actuator constraints is essential for achieving this Convex optimization 3 formulations and solution methodologies have been gaining popularity in the recent pasta fundamental property of convex optimization problems is that any locally minima is also a global minima Moreover, unlike general nonlinear optimization problems numerical techniques for solving these problems are well developed and guarantee convergence Some of the convex optimization formulations are quadratic programming, second-order cone programming, geometric programming, semi-definite programming etc illerson etal 4 have used a convex optimization formulation for the control of spacecraft formation Controller design by posing them as convex optimization problems is similar to predictive control approach discussed in Ref5 Equations of motion governing the full three dimensional translation, rotation and actuator dynamics are derived in Ref 1 he current paper derives a planar version of these equations, however, in closed form and also analytically obtains a linearized version of the same Moreover, equations of motion are obtained directly in terms of inertial co-ordinates of the missile in contrast to the body-frame velocity components his facilitates designing an integrated guidance-controller with this model he linear model is then used in the controller design for which quadratic programming methodology is adopted his paper is organized as follows: Section II deals with the dynamic model of the moving mass actuated missile, Section III deals with the quadratic programming formulation of the controller, closed loop simulation results are presented in Section IV II Model he following assumptions are made in modeling the dynamics of the moving mass actuated missile Motion of the missile is restricted to the vertical plane Gravity vector is assumed to be act along the inertial Z direction he missile is flying outside the atmosphere therefore no aerodynamic forces act on the vehicle Shown in Fig1 is a diagram of the moving mass actuated missile External rocket motor thrust acts along the axis of the missile he position of the mass is manipulated by the controller to control the attitude of the missile By changing the attitude of the missile the orientation of the rocket-motor thrust is changed in an inertial frame, thereby changing the path of the missile Body-frame to inertial frame transformation matrix is given by:- C b = c θ s θ s θ c θ (1) 2 of 13
3 Figure 1 Schematic of the Moving Mass Actuated Missile Velocity of the body in the body-frame is given by:- Ẋ Ṙ b = C b Ż = c θ Ẋ s θ Ż s θ Ẋ + c θ Ż Inertial velocity of the mass in the body-frame is given by:- c θ Ẋ s θ Ż θ Ṙ m = + = s θ Ẋ + c θ Ż + l θ he total kinetic energy of the system can be written as:- (2) c θ Ẋ s θ Ż + θ s θ Ẋ + c θ Ż + + l θ (3) = 1 2 M bẋ M bż I θ m(c θẋ s θż + θ) m(s θẋ + c θż + + l θ) 2 (4) = 1 2 Mt(Ẋ2 + Ż2 ) (I + m2 + ml 2 ) θ m 2 +mc θ Ẋ θ ms θ Ż θ + ms θ Ẋ + ms θ Ẋl θ + mc θ Ż + mc θ Żl θ + m l θ Equations of motion corresponding to the four degrees of freedom X = [X Z θ ] can be written in Lagrangian 6 form as shown below:- d dt Ẋ X = Q X (6) d dt Ż Z = Q Z (7) d dt d dt θ θ = Q θ (8) = Q (9) (5) 3 of 13
4 A detailed evaluation of the above equations is given in the appendix he exact form of the system dynamics is finally obtained as:- where, M = M(X)Ẍ + C(X, Ẋ) = Q (1) M t mc θ + ms θ l ms θ M t mc θ l ms θ mc θ mc θ + ms θ l ms θ + mc θ l (I + m 2 + ml 2 ) ml ms θ mc θ ml m ms θ θ 2 + 2mc θ θ + mcθ l θ 2 mc C = θ θ 2 2ms θ θ msθ l θ 2 2m θ m θ 2 F c θ F Q = s θ + M t g mgs θ + mgc θ l u + mgc θ (11) (12) (13) Ẍ = M(X) 1 [Q C(X, Ẋ)] (14) Inverting the mass matrix using Maple software and simplifying the equations further, the following are obtained:- Ẍ = c θ M t m F c θ Im (M t I + (M t m)m 2 ) F s θ M t m u + mlc θ (M t I + (M t m)m 2 ) u + (15) M t Imlc θ θ2 + 2M t Imc θ θ 2Im 2 c θ θ Im 2 c θ l θ 2 (M t m)(m t I + (M t m)m 2 ) Z = g + c θ u mls θ M t I + (M t m)m 2 u s θ M t m F s θ mi (m M t )(M t I + (M t m)m 2 ) F + (16) m 2 Ils θ θ2 + 2m 2 Is θ θ mlsθ θ2 M t I 2ms θ θmt I (M t m)(m t I + (M t m)m 2 ) m θ = M t I + (M t m)m 2 F M t l M t I + (M t m)m 2 u (17) 2 m2 θ m 2 l θ M t m θ + M t lm θ 2 M t I + (M t m)m 2 = ml M t I + (M t m)m 2 F M t + m(m t m) u + Mt 2 l 2 (M t I + (M t m)m 2 )(M t m) u + (18) + 2ml θ(m t m) (M t I + (M t m)m 2 ) + θ 2 (M t m)ml 2 + (M t I + (M t m)m 2 ) θ 2 he above equations are the exact nonlinear equations governing the translational, rotational and actuator dynamics of the missile hese can be used to design guidance and control laws for the missile Linearizing the above equations about θ =, θ =, =, = the following are obtained:- Ẍ = 1 M t F θ M t m u + ml M t I u (19) 4 of 13
5 Z = g u M t m θ M t m F + θm M t (M t m) F = g u M t m F M t θ (2) θ = m M t I F l I u (21) = ml M t I F M t + m(m t m) u + M t l 2 I(M t m) u (22) III Controller Design Using Quadratic Programming Controller design using quadratic programming will be discussed in this section for a general linear system of the form:- Ẋ = A c X + B c u (23) Discretized version of the above equations obtained using the zero-order-hold are as follows:- X[k + 1] = A d X[k] + B d u[k] (24) where X is a n 1 vector, A d is a n n matrix and B d is a n 1 vector he general form of a quadratic programming problem can be written as follows:- min Z SZ (25) subject to FZ = g ll <= Z <= ul Define Z = [X [2] X [3] X [N] u[1] u[2] u[n 1]] a vector of length (n+1)(n 1) as the optimization state vector Both the state variables and control variables over a N steps are treated as optimization variables he controller objective is posed as the objective function of the optimization problem he constraint matrix is formed by invoking the discrete-time propagation equation over N time steps X[2] X[3] X[N] = A d A d A d X[1] X[2] X[N 1] + B d B d B d u[1] u[2] u[n 1] (26) FZ = g (27) where F is a n(n 1) (n + 1)(N 1) matrix given by the following equation: F = [ ] I n(n 1) n(n 1) n(n 1) (N 1) and g is a n(n 1) 1 vector: B d A d B d A d B d (28) 5 of 13
6 g = A d X[1] State and control constraints such as limits can be equivalently imposed as limits on the optimization state vector (29) Z min <= Z <= Z max (3) IV Results Even though the controller design was done using a discrete-time linear model it was tested on continuous time nonlinear closed loop simulation Equations 1-14 are used for this purpose he optimization process is conducted using the quadprog function of MOSEK 7 software he following values have been used for the missile parameters,m = 17slugs, m = 1775slugs, I = 6slug ft 2, F = 5lbf, max = 5ft, u max = 75lbf Normalized position and control are used for plotting in the following sub section hey are both normalized by their respective maximum values A Body Rate Stabilization θ dot(rad/s) Figure 2 Angular Velocity ime History he state vector corresponding to attitude dynamics is represented by X = associated with these states can be written as:- θ mft M ti θ l I = 1 + mlf t M ti M t M b + Mtl2 IM b [ θ ] Linear dynamics u (31) Continuous to discrete(c2d) time conversion utility in Matlab with zero-order-hold is used to numerically compute discretized version of the system matrices A d and B d 6 of 13
7 Figure 3 Normalized Moving Mass Position Control(lbf) Figure 4 Normalized Control ime History he objective of the controller is to regulate the body-rate of the missile Origin is an equilibrium(unstable) point for the attitude actuator dynamics of the system herefore, the objective function for the optimization problem is chosen to minimize the sum of squares of the body-rate as well as the actuator states Σ N k=1( θ[k] 2 + [k] 2 + [k] 2 ) (32) In other words all entries of the S matrix corresponding to body-rate, actuator position and actuator speed are set to 1 Entries corresponding to the control terms are set to zero A receding horizon controller with constant value of N = 99 is used he lower and upper limits on the optimization state vector are chosen as follows:- max <= [k] <= max, k = 1N 1 (33) u max <= u[k] <= u max, k = 1N 1 (34) No limits have been placed on θ[k] and [k] Initial condition on the angular velocity has been chosen as 1rad/s he primary objective of regulating the body-rate is achieved as seen in Fig2 within 2 seconds he controller avoids both control and state saturation as shown in Fig3 and Fig4 respectively he controller is guaranteed to avoid control saturation in the nonlinear simulation However, it is not theoretically guaranteed to avoid state saturation in the nonlinear simulation his is due to the fact that a linear model is used in the optimization process instead of the actual nonlinear model Closed loop results conducted indicate that 7 of 13
8 the nonlinear simulation also avoids state saturation Primary reason for this is the fact that the control is computed in a feedback manner at each instant of time herefore, to avoid state saturation it is only required that the linear and nonlinear model predictions be close to each other for only time step which is 1 sec in this case B Attitude Command racking Attitude command tracking problem is solved using the optimization approach over a fixed horizon of 3 seconds he sample time of 1 seconds has been retained he control is computed at each instant of time to go time using the state at the current time as the initial condition he parameter N = sample time is also re-computed at each instant of time as a function of the time-to-go where the final time is chosen as 3 seconds he objective function for the optimization problem is chosen as follows:- Σ N k=1(1(θ[k] θ c [k]) 2 + θ[k] 2 + [k] 2 ) (35) Limits on the actuator mass position and actuator force are retained from the previous example Attitude command of θ c [k] = 1degrees has been chosen for this example Initial condition on the attitude angle has been set to 45degrees and the initial condition on body-rate has been set to 5rad/s Shown in figures 5, 6, and 7 are the attitude angle time history, normalized mass position and normalized control histories Clearly, the controller works within the limits of the state and constraints while tracking the attitude command within 3 seconds θ(deg) Figure 5 Attitude Angle ime History Figure 6 Normalized Moving Mass Position 8 of 13
9 Control(lbf) Figure 7 Normalized Control ime History C Integrated Guidance and Control with Seeker Constraints Missile guidance and control have been traditionally addressed in a decoupled fashion An integrated guidance-control system involving the position states, attitude states and actuator states can be designed using approach discussed in this paper he above approach can be easily extended to include guidance component as well Linear model corresponding to the inertial position has been obtained in this paper Inertial frame terminal position co-ordinate specification can be included in the optimization function hus an integrated guidance-controller can be designed using the same optimization approach A free falling target model is assumed Vertical translational dynamics of the missile are augmented to the control model he effect of gravity is eliminated by considering the relative motion dynamics with respect to the target Ẍ target =, Ztarget = g (36) X = X target X, Z = Z target Z (37) he resulting linear dynamical system is of sixth order he linear system matrices are:- Z θ Z θ Ft M t mft M ti = mlf t M ti Z θ Z θ + 1 M t m l I M t M b + Mtl2 IM b u (38) An objective function including the relative position states, missile attitude states and the moving mass actuator states is formed as shown in Eq39 Whereas the attitude and actuator states are penalized over the entire course of interception, the relative position states are penalized only over the terminal phase of interception Σ N k= 2N 3 (1 Z[k] 2 + Ż[k]2 ) + Σ N k=1( θ[k] 2 + 1[k] [k] 2 ) (39) o generate maximum lateral acceleration the missile has a tendency to rotate large angles Off bore sight angle is the angle between the line of sight vector and the longitudinal axis of the missile When the off bore sight angle assumes large values the target leaves the field of view of the seeker herefore, to always keep the target within the field of view of the seeker it is necessary to impose limits on the attitude of the missile hese limits are not constant unlike the position limits on the moving mass hey are dependent on the relative geometry between the missile and the target However, in this example the line of sight angle is 9 of 13
10 so small that limits on off bore sight angle could be approximated as limits on the attitude angle of missile A 45 degree field of view seeler is chosen in this example herefore, the lower and upper limits on the optimization state vector are set to:- 45 degrees <= θ[k] <= 45 degrees, k = 1N 1 (4) max <= [k] <= max, k = 1N 1 (41) u max <= u[k] <= u max, k = 1N 1 (42) Initial conditions on the missile and target are chosen as follows, [X, Z, θ,, Ẋ, Ż, θ, ] = [, 1ft,,, 6ft/s, 1ft/s,, ], [X t, Z t, Ẋ t, Ż t ] = [1ft, 11ft, 5ft/s, ] he value of N is updated at each instant of time based on the instantaneous range and range rate of the missile target separation ime to go is computed as time to go = range range rate and the value of N is computed as N = time to go sample time Shown in Fig8 are the trajectories of the missile and the target he missile successfully intercepts the target from behind Attitude angle time history is shown in Fig9 he controller keeps the attitude of the missile within 45 degrees even in the nonlinear simulation, thereby keeping the target within the field of view of seeker during the entire course of interception Z(Kft) Missile arget X(Kft) Figure 8 arget Interception θ(deg) Figure 9 Attitude Angle ime History 1 of 13
11 Figure 1 Normalized Moving Mass Position Control(lbf) Figure 11 Normalized Control ime History V Conclusion Planar equations of motion governing the dynamics of a moving mass actuated missile have been derived in closed form A linear model was then obtained analytically Different control problems such attituderate regulation, attitude command tracking and target interception are posed as a quadratic programming problems Limitations on the moving mass position, actuator force and seeker field of view are explicitly accounted as optimization state constraints he resulting controller has been effectively tested in nonlinear closed loop simulations 11 of 13
12 Appendix Ṙ b = F b = c θ Ẋ s θ Ż s θ Ẋ + c θ Ż Ẋ = M tẋ + mc θ θ + ms θ + msθ l θ, Ż = M tż ms θ θ + mc θ + mcθ l θ, X = (43) Z = (44) θ = (I + m2 + ml 2 ) θ + mc θ Ẋ ms θ Ż + ms θ Ẋl + mc θ Żl + m l (45) θ = ms θẋ θ mc θ Ż θ + mc θ Ẋ + mc θ Ẋl θ ms θ Ż ms θ Żl θ (46) M b gs θ M b gc θ, Ṙ 1 = = m + ms θ Ẋ + mc θ Ż + ml θ (47) = m θ 2 + mc θ Ẋ θ ms θ Ż θ (48) F F x F x mgs θ, F 1 =, F 2 =, F m = (49) u u + mgc θ c θ Ẋ s θ Ż l θ + s θ Ẋ + c θ Ż, Ṙ 2 = c θ Ẋ s θ Ż + θ s θ Ẋ + c θ Ż + l θ, Ṙ m = c θ Ẋ s θ Ż + + θ s θ Ẋ + c θ Ż + l θ (5) Q X = F b Ṙ b Ẋ + F m Ṙ m Ẋ + F 1 Ṙ 1 Ẋ + F 2 Ṙ 2 Ẋ Q Z = F b Ṙ b Ż + F m Ṙ m Ż + F 1 Ṙ 1 Ż + F 2 Ṙ 2 Ż Q θ = F b Ṙ b θ + F m Ṙ m θ + F 1 Ṙ 1 θ + F 2 Ṙ 2 θ (51) (52) (53) Q X = { Q = F b Ṙ b M b gs θ M b gc θ + F m Ṙ m + F 1 Ṙ 1 + F 2 Ṙ 2 Ṙ b Ẋ = Ṙ 1 Ẋ = Ṙ 2 Ẋ = Ṙ m Ẋ = + F + F x u c θ s θ + (54) (55) F x mgs θ u + mgc θ c θ } s θ (56) Q X = F c θ (57) Ṙ b Ż = Ṙ 1 Ż = Ṙ 2 Ż = Ṙ m Ż = s θ c θ (58) 12 of 13
13 Q Z = { M b gs θ M b gc θ + F + F x u + F x mgs θ u + mgc θ } s θ c θ (59) Q θ = Ṙ b θ =, Ṙ 1 θ = F l Q Z = F s θ + M t g (6) l +, F x u Ṙ 2 θ = l l, + Ṙ m θ = F x mgs θ u + mgc θ l (61) l (62) Q θ = mgs θ + mgc θ l (63) Ṙ b =, Ṙ 1 =, Ṙ 2 =, Ṙ m = [ ] Q = F x mgs θ u + mgc θ 1 1 (64) (65) Q = u + mgc θ (66) References 1 Menon, PK, Sweriduk, G, Ohlmeyer E J, and Malyevac DS, Integrated Guidance and Control of Moving Mass Actuated Kinetic Warheads, Journal of Guidance Control and Dynamics, Vol 27, No 1, Jan-Feb 24, pp Vaddi, S S, Menon, PK, Sweriduk, and Ohlmeyer E J, Multi-Stepping Approach to Finite-Interval Missile Integrated Control, Journal of Guidance Control and Dynamics, Vol 29, No 4, Jul-Aug 26, pp Boyd, S, and Vandenberghe, L, Convex Optimization, Cambridge University Press, 24 4 illerson, M, Inalhan, G, and How J P, Co-Ordination and Control of Distributed Spacecraft Systems Using Convex Optimization echniques, International Journal of Robust and Nonlinear Control, Vol 12, No 2, 22, pp Maciejowski J M, Predictive Control with Constraints, Prentice Hall, 22 6 Likins, P W, Elements of Engineering Mechanics, McGraw-Hill Book Company 7 wwwmosekcom 13 of 13
Linear Feedback Control Using Quasi Velocities
Linear Feedback Control Using Quasi Velocities Andrew J Sinclair Auburn University, Auburn, Alabama 36849 John E Hurtado and John L Junkins Texas A&M University, College Station, Texas 77843 A novel approach
More informationCHAPTER 1. Introduction
CHAPTER 1 Introduction Linear geometric control theory was initiated in the beginning of the 1970 s, see for example, [1, 7]. A good summary of the subject is the book by Wonham [17]. The term geometric
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 informationIntegrator Backstepping using Barrier Functions for Systems with Multiple State Constraints
Integrator Backstepping using Barrier Functions for Systems with Multiple State Constraints Khoi Ngo Dep. Engineering, Australian National University, Australia Robert Mahony Dep. Engineering, Australian
More informationDistributed Coordination and Control of Formation Flying Spacecraft
Distributed Coordination and Control of Formation Flying Spacecraft Michael Tillerson, Louis Breger, and Jonathan P. How MIT Department of Aeronautics and Astronautics {mike t, lbreger, jhow}@mit.edu Abstract
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 information4F3 - Predictive Control
4F3 Predictive Control - Lecture 3 p 1/21 4F3 - Predictive Control Lecture 3 - Predictive Control with Constraints Jan Maciejowski jmm@engcamacuk 4F3 Predictive Control - Lecture 3 p 2/21 Constraints on
More informationPierre Bigot 2 and Luiz C. G. de Souza 3
INTERNATIONAL JOURNAL OF SYSTEMS APPLICATIONS, ENGINEERING & DEVELOPMENT Volume 8, 2014 Investigation of the State Dependent Riccati Equation (SDRE) adaptive control advantages for controlling non-linear
More informationINTEGRATED INTERCEPT MISSILE GUIDANCE AND CONTROL WITH TERMINAL ANGLE CONSTRAINT
6th INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES INTEGRATED INTERCEPT MISSILE GUIDANCE AND CONTROL WITH TERMINAL ANGLE CONSTRAINT H. S. Shin*, T. W. Hwang**, A. Tsourdos***, B. A. White***, M. J.
More informationDesign and modelling of an airship station holding controller for low cost satellite operations
AIAA Guidance, Navigation, and Control Conference and Exhibit 15-18 August 25, San Francisco, California AIAA 25-62 Design and modelling of an airship station holding controller for low cost satellite
More informationEXTREMAL ANALYTICAL CONTROL AND GUIDANCE SOLUTIONS FOR POWERED DESCENT AND PRECISION LANDING. Dilmurat Azimov
EXTREMAL ANALYTICAL CONTROL AND GUIDANCE SOLUTIONS FOR POWERED DESCENT AND PRECISION LANDING Dilmurat Azimov University of Hawaii at Manoa 254 Dole Street, Holmes 22A Phone: (88)-956-2863, E-mail: azimov@hawaii.edu
More informationLecture Schedule Week Date Lecture (M: 2:05p-3:50, 50-N202)
J = x θ τ = J T F 2018 School of Information Technology and Electrical Engineering at the University of Queensland Lecture Schedule Week Date Lecture (M: 2:05p-3:50, 50-N202) 1 23-Jul Introduction + Representing
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 informationFINITE HORIZON ROBUST MODEL PREDICTIVE CONTROL USING LINEAR MATRIX INEQUALITIES. Danlei Chu, Tongwen Chen, Horacio J. Marquez
FINITE HORIZON ROBUST MODEL PREDICTIVE CONTROL USING LINEAR MATRIX INEQUALITIES Danlei Chu Tongwen Chen Horacio J Marquez Department of Electrical and Computer Engineering University of Alberta Edmonton
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 information(W: 12:05-1:50, 50-N202)
2016 School of Information Technology and Electrical Engineering at the University of Queensland Schedule of Events Week Date Lecture (W: 12:05-1:50, 50-N202) 1 27-Jul Introduction 2 Representing Position
More informationAdaptive Control of Hypersonic Vehicles in Presence of Aerodynamic and Center of Gravity Uncertainties
Control of Hypersonic Vehicles in Presence of Aerodynamic and Center of Gravity Uncertainties Amith Somanath and Anuradha Annaswamy Abstract The paper proposes a new class of adaptive controllers that
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 informationThe Inverted Pendulum
Lab 1 The Inverted Pendulum Lab Objective: We will set up the LQR optimal control problem for the inverted pendulum and compute the solution numerically. Think back to your childhood days when, for entertainment
More informationAddis Ababa University Addis Ababa Institute of Technology School Of Mechanical and Industrial Engineering Extension Division` Assignment 1
Assignment 1 1. Vehicle B is stopped at a traffic light, as shown in the figure. At the instant that the light turns green, vehicle B starts to accelerate at 0.9144m/s 2. At this time vehicle A is 91.44m
More informationReal-Time Implementation of a LQR-Based Controller for the Stabilization of a Double Inverted Pendulum
Proceedings of the International MultiConference of Engineers and Computer Scientists 017 Vol I,, March 15-17, 017, Hong Kong Real-Time Implementation of a LQR-Based Controller for the Stabilization of
More informationLecture AC-1. Aircraft Dynamics. Copy right 2003 by Jon at h an H ow
Lecture AC-1 Aircraft Dynamics Copy right 23 by Jon at h an H ow 1 Spring 23 16.61 AC 1 2 Aircraft Dynamics First note that it is possible to develop a very good approximation of a key motion of an aircraft
More informationLine following of a mobile robot
Line following of a mobile robot May 18, 004 1 In brief... The project is about controlling a differential steering mobile robot so that it follows a specified track. Steering is achieved by setting different
More informationLinear Quadratic Integrated vs. Separated Autopilot-Guidance Design
Proceedings of the EuroGNC 2013, 2nd CEAS Specialist Conference on Guidance, Navigation & Control, Delft University of Technology, Delft, The Netherlands, April 10-12, 2013 WeCT2.1 Linear Quadratic Integrated
More informationA Stable Block Model Predictive Control with Variable Implementation Horizon
American Control Conference June 8-,. Portland, OR, USA WeB9. A Stable Block Model Predictive Control with Variable Implementation Horizon Jing Sun, Shuhao Chen, Ilya Kolmanovsky Abstract In this paper,
More informationNONLINEAR BACKSTEPPING DESIGN OF ANTI-LOCK BRAKING SYSTEMS WITH ASSISTANCE OF ACTIVE SUSPENSIONS
NONLINEA BACKSTEPPING DESIGN OF ANTI-LOCK BAKING SYSTEMS WITH ASSISTANCE OF ACTIVE SUSPENSIONS Wei-En Ting and Jung-Shan Lin 1 Department of Electrical Engineering National Chi Nan University 31 University
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 informationECEN 420 LINEAR CONTROL SYSTEMS. Lecture 6 Mathematical Representation of Physical Systems II 1/67
1/67 ECEN 420 LINEAR CONTROL SYSTEMS Lecture 6 Mathematical Representation of Physical Systems II State Variable Models for Dynamic Systems u 1 u 2 u ṙ. Internal Variables x 1, x 2 x n y 1 y 2. y m Figure
More informationTopic # Feedback Control Systems
Topic #1 16.31 Feedback Control Systems Motivation Basic Linear System Response Fall 2007 16.31 1 1 16.31: Introduction r(t) e(t) d(t) y(t) G c (s) G(s) u(t) Goal: Design a controller G c (s) so that the
More informationAdaptive Nonlinear Control Allocation of. Non-minimum Phase Uncertain Systems
2009 American Control Conference Hyatt Regency Riverfront, St. Louis, MO, USA June 10-12, 2009 ThA18.3 Adaptive Nonlinear Control Allocation of Non-minimum Phase Uncertain Systems Fang Liao, Kai-Yew Lum,
More informationNew implicit method for analysis of problems in nonlinear structural dynamics
Applied and Computational Mechanics 5 (2011) 15 20 New implicit method for analysis of problems in nonlinear structural dynamics A. A. Gholampour a,, M. Ghassemieh a a School of Civil Engineering, University
More informationProblem 1 Problem 2 Problem 3 Problem 4 Total
Name Section THE PENNSYLVANIA STATE UNIVERSITY Department of Engineering Science and Mechanics Engineering Mechanics 12 Final Exam May 5, 2003 8:00 9:50 am (110 minutes) Problem 1 Problem 2 Problem 3 Problem
More informationInverted Pendulum. Objectives
Inverted Pendulum Objectives The objective of this lab is to experiment with the stabilization of an unstable system. The inverted pendulum problem is taken as an example and the animation program gives
More informationIn this section of notes, we look at the calculation of forces and torques for a manipulator in two settings:
Introduction Up to this point we have considered only the kinematics of a manipulator. That is, only the specification of motion without regard to the forces and torques required to cause motion In this
More informationAutomatic Control II Computer exercise 3. LQG Design
Uppsala University Information Technology Systems and Control HN,FS,KN 2000-10 Last revised by HR August 16, 2017 Automatic Control II Computer exercise 3 LQG Design Preparations: Read Chapters 5 and 9
More informationMultibody simulation
Multibody simulation Dynamics of a multibody system (Newton-Euler formulation) Dimitar Dimitrov Örebro University June 8, 2012 Main points covered Newton-Euler formulation forward dynamics inverse dynamics
More informationSystem simulation using Matlab, state plane plots
System simulation using Matlab, state plane plots his lab is mainly concerned with making state plane (also referred to as phase plane ) plots for various linear and nonlinear systems with two states he
More informationEQUIVALENT SINGLE-DEGREE-OF-FREEDOM SYSTEM AND FREE VIBRATION
1 EQUIVALENT SINGLE-DEGREE-OF-FREEDOM SYSTEM AND FREE VIBRATION The course on Mechanical Vibration is an important part of the Mechanical Engineering undergraduate curriculum. It is necessary for the development
More information4F3 - Predictive Control
4F3 Predictive Control - Lecture 2 p 1/23 4F3 - Predictive Control Lecture 2 - Unconstrained Predictive Control Jan Maciejowski jmm@engcamacuk 4F3 Predictive Control - Lecture 2 p 2/23 References Predictive
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 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 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 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 informationRobot Dynamics - Rotary Wing UAS: Control of a Quadrotor
Robot Dynamics Rotary Wing AS: Control of a Quadrotor 5-85- V Marco Hutter, Roland Siegwart and Thomas Stastny Robot Dynamics - Rotary Wing AS: Control of a Quadrotor 7..6 Contents Rotary Wing AS. Introduction
More informationWhat is flight dynamics? AE540: Flight Dynamics and Control I. What is flight control? Is the study of aircraft motion and its characteristics.
KING FAHD UNIVERSITY Department of Aerospace Engineering AE540: Flight Dynamics and Control I Instructor Dr. Ayman Hamdy Kassem What is flight dynamics? Is the study of aircraft motion and its characteristics.
More informationAnalytical Dynamics: Lagrange s Equation and its Application A Brief Introduction
Analytical Dynamics: Lagrange s Equation and its Application A Brief Introduction D. S. Stutts, Ph.D. Associate Professor of Mechanical Engineering Missouri University of Science and Technology Rolla,
More informationAgile Missile Controller Based on Adaptive Nonlinear Backstepping Control
Agile Missile Controller Based on Adaptive Nonlinear Backstepping Control Chang-Hun Lee, Tae-Hun Kim and Min-Jea Tahk 3 Korea Advanced Institute of Science and Technology(KAIST), Daejeon, 305-70, Korea
More informationIn the presence of viscous damping, a more generalized form of the Lagrange s equation of motion can be written as
2 MODELING Once the control target is identified, which includes the state variable to be controlled (ex. speed, position, temperature, flow rate, etc), and once the system drives are identified (ex. force,
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 informationDYNAMICS OF SERIAL ROBOTIC MANIPULATORS
DYNAMICS OF SERIAL ROBOTIC MANIPULATORS NOMENCLATURE AND BASIC DEFINITION We consider here a mechanical system composed of r rigid bodies and denote: M i 6x6 inertia dyads of the ith body. Wi 6 x 6 angular-velocity
More informationMultibody simulation
Multibody simulation Dynamics of a multibody system (Euler-Lagrange formulation) Dimitar Dimitrov Örebro University June 16, 2012 Main points covered Euler-Lagrange formulation manipulator inertia matrix
More informationPLANAR KINETIC EQUATIONS OF MOTION (Section 17.2)
PLANAR KINETIC EQUATIONS OF MOTION (Section 17.2) We will limit our study of planar kinetics to rigid bodies that are symmetric with respect to a fixed reference plane. As discussed in Chapter 16, when
More informationNUMERICAL SEARCH OF BOUNDED RELATIVE SATELLITE MOTION
NUMERICAL SEARCH OF BOUNDED RELATIVE SATELLITE MOTION Marco Sabatini 1 Riccardo Bevilacqua 2 Mauro Pantaleoni 3 Dario Izzo 4 1 Ph. D. Candidate, University of Rome La Sapienza, Department of Aerospace
More informationChapter 4 The Equations of Motion
Chapter 4 The Equations of Motion Flight Mechanics and Control AEM 4303 Bérénice Mettler University of Minnesota Feb. 20-27, 2013 (v. 2/26/13) Bérénice Mettler (University of Minnesota) Chapter 4 The Equations
More information1 x(k +1)=(Φ LH) x(k) = T 1 x 2 (k) x1 (0) 1 T x 2(0) T x 1 (0) x 2 (0) x(1) = x(2) = x(3) =
567 This is often referred to as Þnite settling time or deadbeat design because the dynamics will settle in a Þnite number of sample periods. This estimator always drives the error to zero in time 2T or
More informationAutonomous Helicopter Landing A Nonlinear Output Regulation Perspective
Autonomous Helicopter Landing A Nonlinear Output Regulation Perspective Andrea Serrani Department of Electrical and Computer Engineering Collaborative Center for Control Sciences The Ohio State University
More informationStable Limit Cycle Generation for Underactuated Mechanical Systems, Application: Inertia Wheel Inverted Pendulum
Stable Limit Cycle Generation for Underactuated Mechanical Systems, Application: Inertia Wheel Inverted Pendulum Sébastien Andary Ahmed Chemori Sébastien Krut LIRMM, Univ. Montpellier - CNRS, 6, rue Ada
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 informationA Comparison of Closed-Loop Performance of Multirotor Configurations Using Non-Linear Dynamic Inversion Control
Aerospace 2015, 2, 325-352; doi:10.3390/aerospace2020325 OPEN ACCESS aerospace ISSN 2226-4310 www.mdpi.com/journal/aerospace Article A Comparison of Closed-Loop Performance of Multirotor Configurations
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 informationGravity Turn Concept. Curvilinear Coordinate System Gravity Turn Manoeuvre concept Solutions for Constant Pitch Rate
Gravity Turn Concept Curvilinear Coordinate System Gravity Turn Manoeuvre concept Solutions for Constant Pitch Rate Inclined Motion Concept In reality, vertical motion is used only for a very small part
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 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 information1. Type your solutions. This homework is mainly a programming assignment.
THE UNIVERSITY OF TEXAS AT SAN ANTONIO EE 5243 INTRODUCTION TO CYBER-PHYSICAL SYSTEMS H O M E W O R K S # 6 + 7 Ahmad F. Taha October 22, 2015 READ Homework Instructions: 1. Type your solutions. This homework
More informationDynamics. 1 Copyright c 2015 Roderic Grupen
Dynamics The branch of physics that treats the action of force on bodies in motion or at rest; kinetics, kinematics, and statics, collectively. Websters dictionary Outline Conservation of Momentum Inertia
More informationRobust Nonlinear Design of Three Axes Missile Autopilot via Feedback Linearization
Robust Nonlinear Design of Three Axes Missile Autopilot via Feedback Linearization Abhijit Das, Ranajit Das and Siddhartha Mukhopadhyay, Amit Patra 1 1 Abstract The nonlinearity and coupling of the missile
More informationChapter 3 Numerical Methods
Chapter 3 Numerical Methods Part 3 3.4 Differential Algebraic Systems 3.5 Integration of Differential Equations 1 Outline 3.4 Differential Algebraic Systems 3.4.1 Constrained Dynamics 3.4.2 First and Second
More informationIntroduction to Flight Dynamics
Chapter 1 Introduction to Flight Dynamics Flight dynamics deals principally with the response of aerospace vehicles to perturbations in their flight environments and to control inputs. In order to understand
More informationPartial J 2 -Invariance for Spacecraft Formations
Partial J 2 -Invariance for Spacecraft Formations Louis Breger and Jonathan P. How MIT Department of Aeronautics and Astronautics Kyle T. Alfriend Texas A&M University Department of Aerospace Engineering
More information5.12 The Aerodynamic Assist Trajectories of Vehicles Propelled by Solar Radiation Pressure References...
1 The Two-Body Problem... 1 1.1 Position of the Problem... 1 1.2 The Conic Sections and Their Geometrical Properties... 12 1.3 The Elliptic Orbits... 20 1.4 The Hyperbolic and Parabolic Trajectories...
More informationSwinging-Up and Stabilization Control Based on Natural Frequency for Pendulum Systems
9 American Control Conference Hyatt Regency Riverfront, St. Louis, MO, USA June -, 9 FrC. Swinging-Up and Stabilization Control Based on Natural Frequency for Pendulum Systems Noriko Matsuda, Masaki Izutsu,
More informationAS3010: Introduction to Space Technology
AS3010: Introduction to Space Technology L E C T U R E 22 Part B, Lecture 22 19 April, 2017 C O N T E N T S Attitude stabilization passive and active. Actuators for three axis or active stabilization.
More informationA RECEDING HORIZON CONTROL FOR LIFTING PING-PONG BALL. Sumiko Majima and Keii Chou
A RECEDING HORIZON CONROL FOR LIFING PING-PONG BALL Sumiko Majima and Keii Chou Graduate School of Systems and Information Engineering, University of sukuba, sukuba, Ibaraki, Japan Abstract: his paper
More informationAircraft Stability & Control
Aircraft Stability & Control Textbook Automatic control of Aircraft and missiles 2 nd Edition by John H Blakelock References Aircraft Dynamics and Automatic Control - McRuler & Ashkenas Aerodynamics, Aeronautics
More informationAppendix A Solving Linear Matrix Inequality (LMI) Problems
Appendix A Solving Linear Matrix Inequality (LMI) Problems In this section, we present a brief introduction about linear matrix inequalities which have been used extensively to solve the FDI problems described
More informationSpeed Profile Optimization for Optimal Path Tracking
Speed Profile Optimization for Optimal Path Tracking Yiming Zhao and Panagiotis Tsiotras Abstract In this paper, we study the problem of minimumtime, and minimum-energy speed profile optimization along
More information6.3 MODERN GUIDANCE LAWS
Keywords. Modern guidance law, Augmented proportional navigation, Zero-effort-miss 6.3 MODERN GUIDANCE LAWS Essentially the modern guidance laws attempt to take into account target accelerations (i.e.,
More informationArtificial Intelligence & Neuro Cognitive Systems Fakultät für Informatik. Robot Dynamics. Dr.-Ing. John Nassour J.
Artificial Intelligence & Neuro Cognitive Systems Fakultät für Informatik Robot Dynamics Dr.-Ing. John Nassour 25.1.218 J.Nassour 1 Introduction Dynamics concerns the motion of bodies Includes Kinematics
More informationOptimization of Orbital Transfer of Electrodynamic Tether Satellite by Nonlinear Programming
Optimization of Orbital Transfer of Electrodynamic Tether Satellite by Nonlinear Programming IEPC-2015-299 /ISTS-2015-b-299 Presented at Joint Conference of 30th International Symposium on Space Technology
More informationVariable-gain output feedback control
7. Variable-gain output feedback control 7.1. Introduction PUC-Rio - Certificação Digital Nº 611865/CA In designing control laws, the usual first step is to describe the plant at a given operating point
More informationThe dynamics of a Mobile Inverted Pendulum (MIP)
The dynamics of a Mobile Inverted Pendulum (MIP) 1 Introduction Saam Ostovari, Nick Morozovsky, Thomas Bewley UCSD Coordinated Robotics Lab In this document, a Mobile Inverted Pendulum (MIP) is a robotic
More informationM2A2 Problem Sheet 3 - Hamiltonian Mechanics
MA Problem Sheet 3 - Hamiltonian Mechanics. The particle in a cone. A particle slides under gravity, inside a smooth circular cone with a vertical axis, z = k x + y. Write down its Lagrangian in a) Cartesian,
More informationEE Homework 3 Due Date: 03 / 30 / Spring 2015
EE 476 - Homework 3 Due Date: 03 / 30 / 2015 Spring 2015 Exercise 1 (10 points). Consider the problem of two pulleys and a mass discussed in class. We solved a version of the problem where the mass was
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 informationAutonomous navigation of unicycle robots using MPC
Autonomous navigation of unicycle robots using MPC M. Farina marcello.farina@polimi.it Dipartimento di Elettronica e Informazione Politecnico di Milano 7 June 26 Outline Model and feedback linearization
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 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 informationIntegrated Guidance and Control of Missiles with Θ-D Method
Missouri University of Science and Technology Scholars' Mine Mechanical and Aerospace Engineering Faculty Research & Creative Works Mechanical and Aerospace Engineering 11-1-2006 Integrated Guidance and
More informationApplication of the state-dependent riccati equation and kalman filter techniques to the design of a satellite control system
Shock and Vibration 19 (2012) 939 946 939 DOI 10.3233/SAV-2012-0701 IOS Press Application of the state-dependent riccati equation and kalman filter techniques to the design of a satellite control system
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 informationQUALIFYING EXAMINATION, Part 1. Solutions. Problem 1: Mathematical Methods. r r r 2 r r2 = 0 r 2. d 3 r. t 0 e t dt = e t
QUALIFYING EXAMINATION, Part 1 Solutions Problem 1: Mathematical Methods (a) For r > we find 2 ( 1 r ) = 1 ( ) 1 r 2 r r2 = 1 ( 1 ) r r r 2 r r2 = r 2 However for r = we get 1 because of the factor in
More information3 Space curvilinear motion, motion in non-inertial frames
3 Space curvilinear motion, motion in non-inertial frames 3.1 In-class problem A rocket of initial mass m i is fired vertically up from earth and accelerates until its fuel is exhausted. The residual mass
More informationAn Introduction to Model-based Predictive Control (MPC) by
ECE 680 Fall 2017 An Introduction to Model-based Predictive Control (MPC) by Stanislaw H Żak 1 Introduction The model-based predictive control (MPC) methodology is also referred to as the moving horizon
More informationOptimized Trajectory Shaping Guidance for an Air-to-Ground Missile Launched from a Gunship. Craig Phillips Ernie Ohlmeyer Shane Sorenson
Optimized Trajectory Shaping Guidance for an Air-to-Ground Missile Launched from a Gunship Craig Phillips Ernie Ohlmeyer Shane Sorenson Overview Mission Scenario Notional Munition Concept Guidance Laws
More informationSUM-OF-SQUARES BASED STABILITY ANALYSIS FOR SKID-TO-TURN MISSILES WITH THREE-LOOP AUTOPILOT
SUM-OF-SQUARES BASED STABILITY ANALYSIS FOR SKID-TO-TURN MISSILES WITH THREE-LOOP AUTOPILOT Hyuck-Hoon Kwon*, Min-Won Seo*, Dae-Sung Jang*, and Han-Lim Choi *Department of Aerospace Engineering, KAIST,
More informationControlling the Inverted Pendulum
Controlling the Inverted Pendulum Steven A. P. Quintero Department of Electrical and Computer Engineering University of California, Santa Barbara Email: squintero@umail.ucsb.edu Abstract The strategies
More informationInfinite-dimensional nonlinear predictive controller design for open-channel hydraulic systems
Infinite-dimensional nonlinear predictive controller design for open-channel hydraulic systems D. Georges, Control Systems Dept - Gipsa-lab, Grenoble INP Workshop on Irrigation Channels and Related Problems,
More informationON THE INTEGRATION OF EQUATIONS OF MOTION: FEM AND MOLECULAR DYNAMICS PROBLEMS
8th International Congress on Computational Mechanics, Volos, 1-15 July 015 ON THE INTEGRATION OF EQUATIONS OF MOTION: FEM AND MOLECULAR DYNAMICS PROBLEMS E.G. Kakouris, V.K. Koumousis Institute of Structural
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 informationNonholonomic Constraints Examples
Nonholonomic Constraints Examples Basilio Bona DAUIN Politecnico di Torino July 2009 B. Bona (DAUIN) Examples July 2009 1 / 34 Example 1 Given q T = [ x y ] T check that the constraint φ(q) = (2x + siny
More information