Modeling and Analysis of Dynamic Systems
|
|
- Allyson Conley
- 5 years ago
- Views:
Transcription
1 Modeling and Analysis of Dynamic Systems Dr. Guillaume Ducard Fall 2017 Institute for Dynamic Systems and Control ETH Zurich, Switzerland G. Ducard c 1 / 54
2 Outline 1 G. Ducard c 2 / 54
3 Outline 1 G. Ducard c 3 / 54
4 Figure: G. Ducard c 4 / 54
5 Geostationary orbit around the Earth in the Equatorial plane, exactly the sideral Earth rotational speed: 23h 56min 4.1s, from Earth, satellite is seen at exact same location in sky. Very useful for communication purposes. G. Ducard c 5 / 54
6 Importance of modeling for: understanding the basic properties and limitations of satellite dynamics, the design of an orbit-stabilizing controller. Main steps followed: 1 derivation of the satellite motion, 2 linearization around the orbit trajectory, 3 study stability, observability, controllability. G. Ducard c 6 / 54
7 Problem Definition R M r(t) m F r (t) ϕ(t) F ϕ(t) Notations: R: radius of the Earth [m] M: mass of the Earth [kg] m: mass of the satellite [kg] r(t): dist. Earth center to satellite [m] ϕ(t): orbit angle of the satellite [rad] F r (t) : radial force [N] F ϕ (t): tangential force [N] G. Ducard c 7 / 54
8 Problem Definition m F r (t) R M r(t) ϕ(t) F ϕ(t) Assumptions: No other celestial bodies considered (only Earth gravitational force considered). M m C.O.G located at center of Earth. satellite always remains in the equatorial plane 2 variables are sufficient to describe its position: r(t) and ϕ(t). the attitude (orientation) of the satellite is kept constant (by an inner control system) F r (t) and F ϕ (t) independant. G. Ducard c 8 / 54
9 Outline 1 G. Ducard c 9 / 54
10 Building Up the Step 1: Inputs & Outputs Inputs: radial force F r (t) tangential force F ϕ (t) Outputs: Earth to satellite distance r(t) orbit angle: ϕ(t) G. Ducard c 10 / 54
11 Building Up the Step 2: Energies involved in the satellite motion kinetic energy T T(r,ṙ, ϕ) = 1 2 m ṙ m (r ϕ)2 potential energy V: is due to the Earth Gravitational field G Earth G Earth(r) = GM r 2 u F grav (r) = mg Earth Universal gravitation constant: G = m 3 /(s 2 kg), Earth mass: M = kg, Earth radius: R = m. Satellite mass: m M G. Ducard c 11 / 54
12 Building Up the Step 2: Energy needed to bring the satellite from altitude R to r V(r) = r R r F(ρ)dρ = G M m R ρ 2 dρ r 1 = GM m R ρ 2 dρ [ = GM m 1 ] r ρ R ( 1 = GM m R 1 ), r > R r Remark: Potential energy of satellite: - only a function of the distance r, V(r). - is zero at Earth surface (arbitrary but convenient choice) G. Ducard c 12 / 54
13 V(r) = GM m ( 1 R 1 ), r > R r V/V r/r The potential energy: is 0 on the surface of the earth (arbitrary but convenient choice) and reaches 90% of its maximum value: V = GMm R, at a distance of only 10 times the radius of the Earth. G. Ducard c 13 / 54
14 Minimum Speed to Reach Orbit Question: Which minimum speed should the rocket reach to overcome the Earth gravitation force to place the satellite in its orbit? Assume a rapid acceleration at the start from standstill (take off) to a speed v 0, the energy balance is ( mv2 0 = GM m R 1 ) = ( V grav ) r R r kinetic energy that a satellite must have to reach a certain orbit height r R. ( V grav ) R r is the change in potential energy from dist. R to r. The minimum speed needed to reach altitude r is thus: ( 1 v 0 (r) = 2GM R 1 ) r G. Ducard c 14 / 54
15 Escape Velocity : v esc or v Definition 1 It is the minimum velocity that an object should have, such that without additional acceleration (no propulsion anymore), this object will move away from massive body, without being pulled back. Definition 2 It is the velocity of an object, such that its kinetic energy at a certain point in space A is equal its gravitational energy at that point A. Definition 3 The velocity that is required to completely leave the influence of the gravitation field of the Earth. G. Ducard c 15 / 54
16 Escape Velocity How to compute the escape velocity? v = lim r v 0 (r) = 2GM ( 1 R 1 ) r v = 2G M R m/s 11.2km/s 40300km/h G. Ducard c 16 / 54
17 Building Up the Step 3: Lagrange formalism The dynamics of the satellite are formulated using Lagrange s method: [ ] d L L dt ṙ r = F r and d dt [ ] L L ϕ ϕ = F ϕ r with the Lagrange function L = T V (difference of the kinetic and potential energies): L = 1 2 m ṙ2 + 1 ( 1 2 m (r ϕ)2 GM m R 1 ) r G. Ducard c 17 / 54
18 Building Up the L = 1 2 m ṙ m (r ϕ)2 GM m ( 1 R 1 ) r d dt ( L ṙ L ṙ = mṙ ) = m r L r = mr ϕ2 GM m 1 r 2 L ϕ = mr2 ϕ d dt ( ) L = mr 2 ϕ+2mr ϕṙ ϕ L ϕ = 0 m r = m r ϕ 2 G M m 1 r 2 +F r m r 2 ϕ = 2m r ϕ ṙ +F ϕ r G. Ducard c 18 / 54
19 Building Up the In order to simplify the notations, let s introduce the control accelerations yields u r = F r m, respectively u ϕ = F ϕ m r = r ϕ 2 G M 1 r 2 +u r ϕ = 2 ϕṙ 1 r + 1 r u ϕ G. Ducard c 19 / 54
20 Building Up the Circular Orbit For geostationary conditions, the periodic reference orbit must be chosen to be: circular with constant angular velocity, i.e., u r = 0, r = 0, ṙ = 0, r = r 0 u ϕ = 0, ϕ = 0, ϕ = ω 0, ϕ = ω 0 t What are the satellite turn rate ω 0, and the turn radius r 0? G. Ducard c 20 / 54
21 Building Up the : Circular Orbit Orbit turn rate: ω 0 The satellite turn rate is ω 0 = 2π/day taking 1 sidereal day = 23h 56min 4.1 s, you find: ω 0 = rad/s G. Ducard c 21 / 54
22 Building Up the : Circular Orbit Geostationary orbit radius: r 0, (F centripetal = F gravitational ( ) G M 1/3 r 0 = m ω 2 0 Remarks: r 0 is approximately 6.2 times the radius of the earth. The energy required is more than 80% of the escape energy. The resulting tangential speed v ϕ = r 0 ω 0 = 3.06m/s 10800km/h G. Ducard c 22 / 54
23 How to Launch a Satellite? What is an energy efficient trajectory to place a geostationary satellite to orbit? G. Ducard c 23 / 54
24 Outline 1 G. Ducard c 24 / 54
25 Building Up the State-Space Formulation In order not to have second-order ODE to work with, the system is put into a state space form, such that only 1st order time-derivatives appear. x 1 (t) = r, x 2 (t) = ṙ, u 1 (t) = u r x 3 (t) = ϕ, x 4 (t) = ϕ, u 2 (t) = u ϕ x 1 (t) [ ] x 2 (t) u1 (t) x(t) =, u(t) = x 3 (t) u 2 (t) x 4 (t) G. Ducard c 25 / 54
26 Nonlinear ODE where f(x(t)) = d x(t) = f(x(t),u(t)) dt x 2 (t) x 1 x 2 4 (t) G M/x2 1 (t)+u 1(t) x 4 (t) 2x 2 (t) x 4 (t)/x 1 (t)+u 2 (t)/x 1 (t) G. Ducard c 26 / 54
27 Measurements Equations y(t) = h(x(t)) where [ x1 (t)/r 0 h(x(t)) = x 3 (t) Notice: y 1 (t) is the normalized radius. ] G. Ducard c 27 / 54
28 Model Linearization The model is now linearized around the nominal orbit r 0 [ ] 0 0 x 0 (t) =, u ω 0 t 0 (t) = 0 Notice that in this case the nominal trajectory is not an equilibrium point but a periodic solution of the system equations. ω 0 G. Ducard c 28 / 54
29 However, introducing the notation x(t) = x 0 (t)+δx(t), u(t) = u 0 (t)+δu(t), y(t) = y 0 (t)+δy(t) yields d dt x 0(t)+ d dt δx(t) = f(x 0(t)+δx(t),u 0 (t)+δu(t)) f(x 0 (t),u 0 (t))+ f δx(t) x + f u x0 (t),u 0 (t) δu(t) x0 (t),u 0 (t) G. Ducard c 29 / 54
30 Since, by construction d dt x 0(t) = f(x 0 (t),u 0 (t)) the linearized system is described by the equation d f δx(t) = δx(t)+ f dt x u x0 (t),u 0 (t) x0 (t),u 0 (t) δu(t) The computation of the Jacobian matrices yields the following results f x = x G M/x x 1 x x 2 x 4 /x 2 1 u 2/x 2 1 2x 4 /x 1 0 2x 2 /x 1 G. Ducard c 30 / 54
31 f u = /x 1 G. Ducard c 31 / 54
32 Along the periodic nominal solution, we get f 3ω r 0 ω 0 x = x0 (t),u 0 (t) ω 0 /r 0 0 0, G. Ducard c 32 / 54
33 f u = x0 (t),u 0 (t) /r 0 Remark: In this special case the resulting system matrices are time invariant. This is not true in general. When linearizing a nonlinear system around a non-constant equilibrium orbit x e (t),u e (t) the resulting system matrices will be known functions of time {A(t),B(t),C(t),D(t)}. G. Ducard c 33 / 54
34 System in standard notation ẋ(t) = A x(t)+b u(t), y(t) = C x(t) with and A = f x x0 (t),u 0 (t), B = f u x0 (t),u 0 (t) [ ] 1/r C = G. Ducard c 34 / 54
35 For the system {A,B,C,0} defined above the following questions will be discussed below: 1 What are the stability properties of the system? 2 Is the system controllable and observable? 3 Which actuators and sensors are most important? 4 What are the main dynamic properties of the system? G. Ducard c 35 / 54
36 System Stability Eigenvalues of the matrix A are the roots of s ω0 2 s 0 2r 0 ω 0 det(s I A) = det 0 0 s 1 0 2ω 0 /r 0 0 s G. Ducard c 36 / 54
37 System Stability det(s I A) = s det s 1 0 3ω 2 0 s 2r 0 ω 0 0 2ω 0 /r 0 s = s [s (s 2 +4ω 2 0) ( 1) ( 3ω 2 0 s) ] = s 2 (s 2 +ω 2 0) G. Ducard c 37 / 54
38 System Stability det(s I A) = s 2 (s 2 +ω 2 0) The roots are {0,0,+jω 0, jω 0 }. The eigenvalue pair ±jω 0 shows that the system includes oscillatory modes with eigenfrequency equal to the angular speed of the satellite: ω 0. The double root in the origin indicates that the system might be unstable. G. Ducard c 38 / 54
39 One way to see that is to analyze the eigenstructure of the system matrix A. The rank of the matrix M = (s I A) s=0 is rank(m)=3. there is only one eigenvector associated with the double eigenvalue s = 0. Conclusion Therefore, the matrix A is cyclic and, as it was shown in the main text, this means that the linearized system is unstable (behavior of two integrators in series). Remark : The same result is obtained when the system s transfer functions will be computed. G. Ducard c 39 / 54
40 Controllability and Observability Recall of the linear system ẋ = ω r 0 ω ω 0 /r x /r 0 u and y = [ 1/r ] x G. Ducard c 40 / 54
41 Controllability matrix R R = [B,A B,A 2 B,A 3 B] = ω /r /r 0 2ω 0 /r This system is completely controllable. In fact the four first columns already are linearly independent (the determinant of that submatrix is 1/r 2 0 ). G. Ducard c 41 / 54
42 Failure Scenario 1: Radial Thruster Failure This yields input matrix b 2 = [0, 0, 0, 1/r 0 ] T. The corresponding controllability matrix has the form 0 0 2ω ω 0 0 2ω 3 0 R 2 = 0 1/r 0 0 4ω 0 /r 0 1/r 0 0 4ω0 2/r 0 0 The determinant of this matrix is 12ω0 4/r2 0, i.e., not zero on the reference orbit. Accordingly, the satellite remains completely controllable even if the radial thruster fails. G. Ducard c 42 / 54
43 Failure Scenario 2: Tangential Thruster Failure This yields b 1 = [0, 1, 0, 0] T and the controllability matrix R 1 = ω ω ω 0 /r ω 0 /r 0 0 2ω 3 0 /r 0 in this case has only three independent columns (the fourth column is equal to ω0 2 times the second column), i.e., In case of tangential-thruster failure the system is no longer completely controllable. G. Ducard c 43 / 54
44 Observability O = C C A C A 2 C A 3 = 1/r /r The system with no faults is completely observable. G. Ducard c 44 / 54
45 Fault Scenario 1: Radial Sensor Failure r meas If the radial sensor fails, the measurement matrix reduces to c 2 = [0, 0, 1, 0] and the resulting observability matrix O 2 has the form O 2 = 0 2ω 0 /r ω0 3/r ω0 2 This matrix is regular (its determinant is 12ω0 4/r2 0 ). Accordingly, the system remains completely observable even if the radial sensor fails. G. Ducard c 45 / 54
46 Fault Scenario 2: Tangential Sensor Failure ϕ meas If the tangential sensor fails, the measurement matrix reduces to c 1 = [1, 0, 0, 0] and the resulting observability matrix O 1 has the form 1/r /r O 1 = 3ω 0 2/r r 0 ω 0 0 ω0 2/r Obviously, this matrix is singular (the third row is zero). Accordingly, the system is not completely observable in this case. G. Ducard c 46 / 54
47 Conclusions on the satellite design This analysis is helpful for the design of the satellite. The main result is: tangential actuator and sensor are more important than radial ones. If radial actuator and sensor fail, the satellite stability is still guaranteed (control reconfiguration). G. Ducard c 47 / 54
48 Transfer Function P(s) = P(s) = [ P11 (s) P 12 (s) P 21 (s) P 22 (s) 1 s 2 +ω 2 0 ] = C(sI A) 1 B 2ω 0 s (s 2 +ω 2 0) 2ω 0 r 0 s (s 2 +ω 2 0 ) s 2 3ω 2 0 r 0 s 2 (s 2 +ω 2 0 ) G. Ducard c 48 / 54
49 Note that even if the system is stabilizable with only the tangential thruster and sensor working, the control problem is difficult, because the corresponding SISO system transfer function P 22 (s) has a non-minimumphase zero at 3 ω 0 (limits the attainable crossover frequency to about 0.85 ω 0. G. Ducard c 49 / 54
50 Remark: The MIMO system has no finite transmission zeros. This can be seen from [ ] (s I A) B Z(s) = C D which turns out to be detz(s) = 1 r 0 Accordingly, no finite s yields a determinant equal to zero and, therefore, no finite zero exists. In other words, the MIMO system P(s) is minimumphase. Compared to the SISO design using only P 22 (s), the design of a suitable controller using both inputs and outputs will yield a much better closed-loop system performance in terms of bandwidth, response times, etc. G. Ducard c 50 / 54
51 Two simulations The main objective is to convey a first impression of the differences between the linear system behavior, as analyzed above, and the nonlinear system behavior. The system is assumed to be in its reference state for all t < 0. At that moment a test signal { d0 0 t < 1 d(t) = 0 1 t is applied. In a first case d 0 is chosen rather small and in a second case very large, in order to clearly show the differences between the linear an the nonlinear system behavior. G. Ducard c 51 / 54
52 r(t) r0 ϕ(t) ω0t Input: u r = d(t), u ϕ = time t/t Input: u r = d(t), u ϕ= time t/t 0 r(t) r0 ϕ(t) ω0t Input: u r = 0, u ϕ = d(t) time t/t 0 Input: u r = 0, u ϕ = d(t) time t/t 0 G. Ducard c 52 / 54
53 r(t) r0 ϕ(t) ω0t Input: u r = d(t), u ϕ = time t/t Input: u r = d(t), u ϕ= time t/t 0 r(t) r0 ϕ(t) ω0t Input: u r = 0, u ϕ = d(t) time t/t Input: u r = 0, u ϕ = d(t) time t/t 0 G. Ducard c 53 / 54
54 Next lecture + Upcoming Exercise Next lecture Zero dynamics Nonlinear systems Thank you for your attention. G. Ducard c 54 / 54
Modeling and Analysis of Dynamic Systems
Modeling and Analysis of Dynamic Systems Dr. Guillaume Ducard Fall 2017 Institute for Dynamic Systems and Control ETH Zurich, Switzerland G. Ducard c 1 / 57 Outline 1 Lecture 13: Linear System - Stability
More informationÜbersetzungshilfe / Translation aid (English) To be returned at the end of the exam!
Prüfung Regelungstechnik I (Control Systems I) Prof. Dr. Lino Guzzella 3. 8. 24 Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam! Do not mark up this translation aid
More informationModeling and Analysis of Dynamic Systems
Modeling and Analysis of Dynamic Systems by Dr. Guillaume Ducard Fall 2016 Institute for Dynamic Systems and Control ETH Zurich, Switzerland based on script from: Prof. Dr. Lino Guzzella 1/33 Outline 1
More information[2] Explain why values of gravitational potential near to an isolated mass are all negative [3]
1 (a) Define gravitational potential.... [2] (b) Explain why values of gravitational potential near to an isolated mass are all negative.... [3] (c) The Earth may be assumed to be an isolated sphere of
More information1 Controllability and Observability
1 Controllability and Observability 1.1 Linear Time-Invariant (LTI) Systems State-space: Dimensions: Notation Transfer function: ẋ = Ax+Bu, x() = x, y = Cx+Du. x R n, u R m, y R p. Note that H(s) is always
More informationModeling and Analysis of Dynamic Systems
Modeling and Analysis of Dynamic Systems by Dr. Guillaume Ducard Fall 2016 Institute for Dynamic Systems and Control ETH Zurich, Switzerland 1/21 Outline 1 Lecture 4: Modeling Tools for Mechanical Systems
More informationGravitation. chapter 9
chapter 9 Gravitation Circular orbits (Section 9.3) 1, 2, and 3 are simple exercises to supplement the quantitative calculations of Examples 4, 5, and 6 in Section 9.3. 1. Satellite near Earth s surface
More informationImplementation of a Communication Satellite Orbit Controller Design Using State Space Techniques
ASEAN J Sci Technol Dev, 29(), 29 49 Implementation of a Communication Satellite Orbit Controller Design Using State Space Techniques M T Hla *, Y M Lae 2, S L Kyaw 3 and M N Zaw 4 Department of Electronic
More informationNewton s Gravitational Law
1 Newton s Gravitational Law Gravity exists because bodies have masses. Newton s Gravitational Law states that the force of attraction between two point masses is directly proportional to the product of
More informationPoincaré Map, Floquet Theory, and Stability of Periodic Orbits
Poincaré Map, Floquet Theory, and Stability of Periodic Orbits CDS140A Lecturer: W.S. Koon Fall, 2006 1 Poincaré Maps Definition (Poincaré Map): Consider ẋ = f(x) with periodic solution x(t). Construct
More informationModeling and Analysis of Dynamic Systems
Modeling and Analysis of Dynamic Systems Dr. Guillaume Ducard Fall 2017 Institute for Dynamic Systems and Control ETH Zurich, Switzerland G. Ducard c 1 / 21 Outline 1 Lecture 4: Modeling Tools for Mechanical
More informationChapter 9. Gravitation
Chapter 9 Gravitation 9.1 The Gravitational Force For two particles that have masses m 1 and m 2 and are separated by a distance r, the force has a magnitude given by the same magnitude of force acts on
More informationRepresent this system in terms of a block diagram consisting only of. g From Newton s law: 2 : θ sin θ 9 θ ` T
Exercise (Block diagram decomposition). Consider a system P that maps each input to the solutions of 9 4 ` 3 9 Represent this system in terms of a block diagram consisting only of integrator systems, represented
More informationWeightlessness and satellites in orbit. Orbital energies
Weightlessness and satellites in orbit Orbital energies Review PE = - GMm R v escape = 2GM E R = 2gR E Keppler s law: R3 = GM s T 2 4π 2 Orbital Motion Orbital velocity escape velocity In orbital motion
More informationPHYSICS CLASS XI CHAPTER 8 GRAVITATION
PHYSICS CLASS XI CHAPTER 8 GRAVITATION Q.1. Can we determine the mass of a satellite by measuring its time period? Ans. No, we cannot determine the mass of a satellite by measuring its time period. Q.2.
More informationSystems and Control Theory Lecture Notes. Laura Giarré
Systems and Control Theory Lecture Notes Laura Giarré L. Giarré 2017-2018 Lesson 5: State Space Systems State Dimension Infinite-Dimensional systems State-space model (nonlinear) LTI State Space model
More informationPhysics Mechanics Lecture 30 Gravitational Energy
Physics 170 - Mechanics Lecture 30 Gravitational Energy Gravitational Potential Energy Gravitational potential energy of an object of mass m a distance r from the Earth s center: Gravitational Potential
More informationApr 14, Calculus with Algebra and Trigonometry II Lecture 20More physics Aprapplications
Calculus with Algebra and Trigonometry II Lecture 20 More physics applications Apr 14, 2015 14, 2015 1 / 14 Motion in two dimensions A particle s motion can be described by specifying how the coordinates
More information10/21/2003 PHY Lecture 14 1
Announcements. Second exam scheduled for Oct. 8 th -- practice exams now available -- http://www.wfu.edu/~natalie/f03phy3/extrapractice/. Thursday review of Chapters 9-4 3. Today s lecture Universal law
More informationPhysics H7A, Fall 2011 Homework 4 Solutions
Physics H7A, Fall 20 Homework 4 Solutions. (K&K Problem 2.) A mass m is connected to a vertical revolving axle by two strings of length l, each making an angle of 45 with the axle, as shown. Both the axle
More informationAssessment Schedule 2011 Physics: Demonstrate understanding of mechanical systems (90521)
NCEA Level 3 Physics (9052) 20 page of 5 Assessment Schedule 20 Physics: Demonstrate understanding of mechanical systems (9052) Evidence Statement Q Evidence Achievement Merit Excellence ONE (a) f = 4500
More information7 - GRAVITATION Page 1 ( Answers at the end of all questions )
7 - GRAVITATION Page 1 1 ) The change in the value of g at a height h above the surface of the earth is the same as at a depth d below the surface of earth. When both d and h are much smaller than the
More informationwhere s is the horizontal range, in this case 50 yards. Note that as long as the initial speed of the bullet is great enough to let it hit a target at
1 PHYS 31 Homework Assignment Due Friday, 13 September 00 1. A monkey is hanging from a tree limbat height h above the ground. A hunter 50 yards away from the base of the tree sees him, raises his gun
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 informationGravitation and Newton s Synthesis
Lecture 10 Chapter 6 Physics I 0.4.014 Gravitation and Newton s Synthesis Course website: http://faculty.uml.edu/andriy_danylov/teaching/physicsi Lecture Capture: http://echo360.uml.edu/danylov013/physics1spring.html
More informationRadial Acceleration. recall, the direction of the instantaneous velocity vector is tangential to the trajectory
Radial Acceleration recall, the direction of the instantaneous velocity vector is tangential to the trajectory 1 Radial Acceleration recall, the direction of the instantaneous velocity vector is tangential
More informationLecture 10. Example: Friction and Motion
Lecture 10 Goals: Exploit Newton s 3 rd Law in problems with friction Employ Newton s Laws in 2D problems with circular motion Assignment: HW5, (Chapter 7, due 2/24, Wednesday) For Tuesday: Finish reading
More informationSolution Set Two. 1 Problem #1: Projectile Motion Cartesian Coordinates Polar Coordinates... 3
: Solution Set Two Northwestern University, Classical Mechanics Classical Mechanics, Third Ed.- Goldstein October 7, 2015 Contents 1 Problem #1: Projectile Motion. 2 1.1 Cartesian Coordinates....................................
More informationPhysics 201, Lecture 23
Physics 201, Lecture 23 Today s Topics n Universal Gravitation (Chapter 13) n Review: Newton s Law of Universal Gravitation n Properties of Gravitational Field (13.4) n Gravitational Potential Energy (13.5)
More informationStatic Equilibrium, Gravitation, Periodic Motion
This test covers static equilibrium, universal gravitation, and simple harmonic motion, with some problems requiring a knowledge of basic calculus. Part I. Multiple Choice 1. 60 A B 10 kg A mass of 10
More information(b) The period T and the angular frequency ω of uniform rotation are related to the cyclic frequency f as. , ω = 2πf =
PHY 302 K. Solutions for problem set #9. Non-textbook problem #1: (a) Rotation frequency of 1 Hz means one revolution per second, or 60 revolutions per minute (RPM). The pre-lp vinyl disks rotated at 78
More informationReview of Linear Momentum And Rotational Motion
Physics 7B-1 (A/B) Professor Cebra Winter 2010 Lecture 7 Review of Linear Momentum And Rotational Motion Slide 1 of 29 Physics 7B Lecture 7 17-Feb-2010 Slide 2 of 29 The Definition of Impulse Recall that
More informationDistance travelled time taken and if the particle is a distance s(t) along the x-axis, then its instantaneous speed is:
Chapter 1 Kinematics 1.1 Basic ideas r(t) is the position of a particle; r = r is the distance to the origin. If r = x i + y j + z k = (x, y, z), then r = r = x 2 + y 2 + z 2. v(t) is the velocity; v =
More informationLecture 6, September 1, 2017
Engineering Mathematics Fall 07 Lecture 6, September, 07 Escape Velocity Suppose we have a planet (or any large near to spherical heavenly body) of radius R and acceleration of gravity at the surface of
More informationChapter 4 Making Sense of the Universe: Understanding Motion, Energy, and Gravity. Copyright 2009 Pearson Education, Inc.
Chapter 4 Making Sense of the Universe: Understanding Motion, Energy, and Gravity How do we describe motion? Precise definitions to describe motion: Speed: Rate at which object moves speed = distance time
More informationLecture 7 (Weeks 13-14)
Lecture 7 (Weeks 13-14) Introduction to Multivariable Control (SP - Chapters 3 & 4) Eugenio Schuster schuster@lehigh.edu Mechanical Engineering and Mechanics Lehigh University Lecture 7 (Weeks 13-14) p.
More information2.10 Saddles, Nodes, Foci and Centers
2.10 Saddles, Nodes, Foci and Centers In Section 1.5, a linear system (1 where x R 2 was said to have a saddle, node, focus or center at the origin if its phase portrait was linearly equivalent to one
More informationTwo-Body Problem. Central Potential. 1D Motion
Two-Body Problem. Central Potential. D Motion The simplest non-trivial dynamical problem is the problem of two particles. The equations of motion read. m r = F 2, () We already know that the center of
More informationVISUAL PHYSICS ONLINE
VISUAL PHYSICS ONLINE EXCEL SIMULATION MOTION OF SATELLITES DOWNLOAD the MS EXCEL program PA50satellite.xlsx and view the worksheet Display as shown in the figure below. One of the most important questions
More information28. Pendulum phase portrait Draw the phase portrait for the pendulum (supported by an inextensible rod)
28. Pendulum phase portrait Draw the phase portrait for the pendulum (supported by an inextensible rod) θ + ω 2 sin θ = 0. Indicate the stable equilibrium points as well as the unstable equilibrium points.
More informationClassical Mechanics Comprehensive Exam Solution
Classical Mechanics Comprehensive Exam Solution January 31, 011, 1:00 pm 5:pm Solve the following six problems. In the following problems, e x, e y, and e z are unit vectors in the x, y, and z directions,
More information6.241 Dynamic Systems and Control
6.241 Dynamic Systems and Control Lecture 24: H2 Synthesis Emilio Frazzoli Aeronautics and Astronautics Massachusetts Institute of Technology May 4, 2011 E. Frazzoli (MIT) Lecture 24: H 2 Synthesis May
More informationis conserved, calculating E both at θ = 0 and θ = π/2 we find that this happens for a value ω = ω given by: 2g
UNIVERSITETET I STAVANGER Institutt for matematikk og naturvitenskap Suggested solutions, FYS 500 Classical Mechanics Theory 2016 fall Set 5 for 23. September 2016 Problem 27: A string can only support
More informationControl, Stabilization and Numerics for Partial Differential Equations
Paris-Sud, Orsay, December 06 Control, Stabilization and Numerics for Partial Differential Equations Enrique Zuazua Universidad Autónoma 28049 Madrid, Spain enrique.zuazua@uam.es http://www.uam.es/enrique.zuazua
More informationÜbersetzungshilfe / Translation aid (English) To be returned at the end of the exam!
Prüfung Regelungstechnik I (Control Systems I) Prof. Dr. Lino Guzzella 9. 8. 2 Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam! Do not mark up this translation aid -
More informationDynamical Systems & Lyapunov Stability
Dynamical Systems & Lyapunov Stability Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University Outline Ordinary Differential Equations Existence & uniqueness Continuous dependence
More informationWelcome back to Physics 215. Review gravity Oscillations Simple harmonic motion
Welcome back to Physics 215 Review gravity Oscillations Simple harmonic motion Physics 215 Spring 2018 Lecture 14-1 1 Final Exam: Friday May 4 th 5:15-7:15pm Exam will be 2 hours long Have an exam buddy
More informationF 12. = G m m 1 2 F 21. = G m 1m 2 = F 12. Review: Newton s Law Of Universal Gravitation. Physics 201, Lecture 23. g As Function of Height
Physics 01, Lecture Today s Topics n Universal Gravitation (Chapter 1 n Review: Newton s Law of Universal Gravitation n Properties of Gravitational Field (1.4 n Gravitational Potential Energy (1.5 n Escape
More informationUniversal Gravitation
Universal Gravitation Newton s Law of Universal Gravitation Every particle in the Universe attracts every other particle with a force that is directly proportional to the product of their masses and inversely
More informationLecture 9: Eigenvalues and Eigenvectors in Classical Mechanics (See Section 3.12 in Boas)
Lecture 9: Eigenvalues and Eigenvectors in Classical Mechanics (See Section 3 in Boas) As suggested in Lecture 8 the formalism of eigenvalues/eigenvectors has many applications in physics, especially in
More informationModeling and Analysis of Dynamic Systems
Modeling and Analysis of Dynamic Systems by Dr. Guillaume Ducard Fall 2016 Institute for Dynamic Systems and Control ETH Zurich, Switzerland 1/22 Outline 1 Lecture 5: Hydraulic Systems Pelton Turbine:
More informationGravitation. One mark questions.
One mark questions. Gravitation. 01. What is the ratio of force of attraction between two bodies kept in air and the same distance apart in water? Ans- 1:1, because gravitational force is independent of
More informationChapter 10 Rotational Kinematics and Energy. Copyright 2010 Pearson Education, Inc.
Chapter 10 Rotational Kinematics and Energy 10-1 Angular Position, Velocity, and Acceleration 10-1 Angular Position, Velocity, and Acceleration Degrees and revolutions: 10-1 Angular Position, Velocity,
More informationv lim a t = d v dt a n = v2 R curvature
PHY 02 K. Solutions for Problem set # 6. Textbook problem 5.27: The acceleration vector a of the particle has two components, the tangential acceleration a t d v dt v lim t 0 t (1) parallel to the velocity
More informationChapter 6 Gravitation and Newton s Synthesis
Chapter 6 Gravitation and Newton s Synthesis If the force of gravity is being exerted on objects on Earth, what is the origin of that force? Newton s realization was that the force must come from the Earth.
More informationSatellite meteorology
GPHS 422 Satellite meteorology GPHS 422 Satellite meteorology Lecture 1 6 July 2012 Course outline 2012 2 Course outline 2012 - continued 10:00 to 12:00 3 Course outline 2012 - continued 4 Some reading
More informationSession 6: Analytical Approximations for Low Thrust Maneuvers
Session 6: Analytical Approximations for Low Thrust Maneuvers As mentioned in the previous lecture, solving non-keplerian problems in general requires the use of perturbation methods and many are only
More informationINTRODUCTION: Ptolemy geo-centric theory Nicolas Copernicus Helio-centric theory TychoBrahe Johannes Kepler
INTRODUCTION: Ptolemy in second century gave geo-centric theory of planetary motion in which the Earth is considered stationary at the centre of the universe and all the stars and the planets including
More informationCircular Motion and Gravitation. Centripetal Acceleration
Circular Motion and Gravitation Centripetal Acceleration Recall linear acceleration! Δv! aavg t 3. Going around urve, at constant speed 1. Speeding up vi vi Δv a ac ac vi ac. Slowing down v velocity and
More informationGravitation & Kepler s Laws
Gravitation & Kepler s Laws What causes YOU to be pulled down to the surface of the earth? THE EARTH.or more specifically the EARTH S MASS. Anything that has MASS has a gravitational pull towards it. F
More informationFall 線性系統 Linear Systems. Chapter 08 State Feedback & State Estimators (SISO) Feng-Li Lian. NTU-EE Sep07 Jan08
Fall 2007 線性系統 Linear Systems Chapter 08 State Feedback & State Estimators (SISO) Feng-Li Lian NTU-EE Sep07 Jan08 Materials used in these lecture notes are adopted from Linear System Theory & Design, 3rd.
More informationPhysics 12. Unit 5 Circular Motion and Gravitation Part 2
Physics 12 Unit 5 Circular Motion and Gravitation Part 2 1. Newton s law of gravitation We have seen in Physics 11 that the force acting on an object due to gravity is given by a well known formula: F
More information9.3 Worked Examples Circular Motion
9.3 Worked Examples Circular Motion Example 9.1 Geosynchronous Orbit A geostationary satellite goes around the earth once every 3 hours 56 minutes and 4 seconds, (a sidereal day, shorter than the noon-to-noon
More informationGravitational Potential Energy and Total Energy *
OpenStax-CNX module: m58347 Gravitational Potential Energy and Total Energy * OpenStax This work is produced by OpenStax-CNX and licensed under the Creative Commons Attribution License 4.0 By the end of
More informationInvariant Manifolds of Dynamical Systems and an application to Space Exploration
Invariant Manifolds of Dynamical Systems and an application to Space Exploration Mateo Wirth January 13, 2014 1 Abstract In this paper we go over the basics of stable and unstable manifolds associated
More informationLecture 16. Gravitation
Lecture 16 Gravitation Today s Topics: The Gravitational Force Satellites in Circular Orbits Apparent Weightlessness lliptical Orbits and angular momentum Kepler s Laws of Orbital Motion Gravitational
More informationLecture 13 REVIEW. Physics 106 Spring What should we know? What should we know? Newton s Laws
Lecture 13 REVIEW Physics 106 Spring 2006 http://web.njit.edu/~sirenko/ What should we know? Vectors addition, subtraction, scalar and vector multiplication Trigonometric functions sinθ, cos θ, tan θ,
More informationPhysics 207 Lecture 10. Lecture 10. Employ Newton s Laws in 2D problems with circular motion
Lecture 10 Goals: Employ Newton s Laws in 2D problems with circular motion Assignment: HW5, (Chapters 8 & 9, due 3/4, Wednesday) For Tuesday: Finish reading Chapter 8, start Chapter 9. Physics 207: Lecture
More informationThe Calculus of Vec- tors
Physics 2460 Electricity and Magnetism I, Fall 2007, Lecture 3 1 The Calculus of Vec- Summary: tors 1. Calculus of Vectors: Limits and Derivatives 2. Parametric representation of Curves r(t) = [x(t), y(t),
More informationGeorgia Institute of Technology Nonlinear Controls Theory Primer ME 6402
Georgia Institute of Technology Nonlinear Controls Theory Primer ME 640 Ajeya Karajgikar April 6, 011 Definition Stability (Lyapunov): The equilibrium state x = 0 is said to be stable if, for any R > 0,
More informationLesson 7. Luis Anchordoqui. Physics 168. Tuesday, October 10, 17
Lesson 7 Physics 168 1 Eruption of a large volcano on Jupiter s moon When volcano erupts speed of effluence exceeds escape speed of Io and so a stream of particles is projected into space Material in stream
More informationAngular Motion, General Notes
Angular Motion, General Notes! When a rigid object rotates about a fixed axis in a given time interval, every portion on the object rotates through the same angle in a given time interval and has the same
More information1. Find the solution of the following uncontrolled linear system. 2 α 1 1
Appendix B Revision Problems 1. Find the solution of the following uncontrolled linear system 0 1 1 ẋ = x, x(0) =. 2 3 1 Class test, August 1998 2. Given the linear system described by 2 α 1 1 ẋ = x +
More informationExplain how it is possible for the gravitational force to cause the satellite to accelerate while its speed remains constant.
YEAR 12 PHYSICS: GRAVITATION PAST EXAM QUESTIONS Name: QUESTION 1 (1995 EXAM) (a) State Newton s Universal Law of Gravitation in words (b) A satellite of mass (m) moves in orbit of a planet with mass (M).
More information9.2 GRAVITATIONAL FIELD, POTENTIAL, AND ENERGY 9.4 ORBITAL MOTION HW/Study Packet
9.2 GRAVITATIONAL FIELD, POTENTIAL, AND ENERGY 9.4 ORBITAL MOTION HW/Study Packet Required: READ Hamper pp 186-191 HL Supplemental: Cutnell and Johnson, pp 144-150, 169-170 Tsokos, pp 143-152 Tsokos Questions
More informationZeros and zero dynamics
CHAPTER 4 Zeros and zero dynamics 41 Zero dynamics for SISO systems Consider a linear system defined by a strictly proper scalar transfer function that does not have any common zero and pole: g(s) =α p(s)
More informationSteve Smith Tuition: Physics Notes
Steve Smith Tuition: Physics Notes E = mc 2 F = GMm sin θ m = mλ d hν = φ + 1 2 mv2 Static Fields IV: Gravity Examples Contents 1 Gravitational Field Equations 3 1.1 adial Gravitational Field Equations.................................
More informationPhysics 2514 Lecture 22
Physics 2514 Lecture 22 P. Gutierrez Department of Physics & Astronomy University of Oklahoma Physics 2514 p. 1/15 Information Information needed for the exam Exam will be in the same format as the practice
More informationMidterm 3 Review (Ch 9-14)
Midterm 3 Review (Ch 9-14) PowerPoint Lectures for University Physics, Twelfth Edition Hugh D. Young and Roger A. Freedman Lectures by James Pazun Copyright 2008 Pearson Education Inc., publishing as Pearson
More informationPHY4116 From Newton to Einstein. Mid-Term Test, 10a.m. Wed. 16 th Nov Duration: 50 minutes. There are 25 marks in Section A and 25 in Section B.
PHY46 From Newton to Einstein Mid-Term Test, 0a.m. Wed. 6 th Nov. 06 Duration: 50 minutes. There are 5 marks in Section A and 5 in Section B. Use g = 0 ms in numerical calculations. You ma use the following
More information2 Lyapunov Stability. x(0) x 0 < δ x(t) x 0 < ɛ
1 2 Lyapunov Stability Whereas I/O stability is concerned with the effect of inputs on outputs, Lyapunov stability deals with unforced systems: ẋ = f(x, t) (1) where x R n, t R +, and f : R n R + R n.
More informationLecture Outline Chapter 10. Physics, 4 th Edition James S. Walker. Copyright 2010 Pearson Education, Inc.
Lecture Outline Chapter 10 Physics, 4 th Edition James S. Walker Chapter 10 Rotational Kinematics and Energy Units of Chapter 10 Angular Position, Velocity, and Acceleration Rotational Kinematics Connections
More informationExam. 135 minutes, 15 minutes reading time
Exam August 6, 208 Control Systems II (5-0590-00) Dr. Jacopo Tani Exam Exam Duration: 35 minutes, 5 minutes reading time Number of Problems: 35 Number of Points: 47 Permitted aids: 0 pages (5 sheets) A4.
More informationSatellite Orbital Maneuvers and Transfers. Dr Ugur GUVEN
Satellite Orbital Maneuvers and Transfers Dr Ugur GUVEN Orbit Maneuvers At some point during the lifetime of most space vehicles or satellites, we must change one or more of the orbital elements. For example,
More informationLecture plan: Control Systems II, IDSC, 2017
Control Systems II MAVT, IDSC, Lecture 8 28/04/2017 G. Ducard Lecture plan: Control Systems II, IDSC, 2017 SISO Control Design 24.02 Lecture 1 Recalls, Introductory case study 03.03 Lecture 2 Cascaded
More informationExam. 135 minutes, 15 minutes reading time
Exam August 15, 2017 Control Systems I (151-0591-00L) Prof Emilio Frazzoli Exam Exam Duration: 135 minutes, 15 minutes reading time Number of Problems: 44 Number of Points: 52 Permitted aids: Important:
More informationAP Physics Multiple Choice Practice Gravitation
AP Physics Multiple Choice Practice Gravitation 1. Each of five satellites makes a circular orbit about an object that is much more massive than any of the satellites. The mass and orbital radius of each
More informationLecture for Week 6 (Secs ) Derivative Miscellany I
Lecture for Week 6 (Secs. 3.6 9) Derivative Miscellany I 1 Implicit differentiation We want to answer questions like this: 1. What is the derivative of tan 1 x? 2. What is dy dx if x 3 + y 3 + xy 2 + x
More informationACM/CMS 107 Linear Analysis & Applications Fall 2016 Assignment 4: Linear ODEs and Control Theory Due: 5th December 2016
ACM/CMS 17 Linear Analysis & Applications Fall 216 Assignment 4: Linear ODEs and Control Theory Due: 5th December 216 Introduction Systems of ordinary differential equations (ODEs) can be used to describe
More informationGetting Some Big Air
Getting Some Big Air Control #10499 February 14, 2011 Abstract In this paper we address the problem of optimizing a ramp for a snowboarder to travel. Our approach is two sided. We first address the forward
More informationChapter 8 - Gravity Tuesday, March 24 th
Chapter 8 - Gravity Tuesday, March 24 th Newton s law of gravitation Gravitational potential energy Escape velocity Kepler s laws Demonstration, iclicker and example problems We are jumping backwards to
More informationMechanics III - Dynamics
Institute of Mechanical Systems Chair in Nonlinear Dynamics Prof. George Haller Mechanics III - Dynamics Institute of Mechanical Systems Swiss Federal Institute of Technology (ETH) Zurich ased on the lecture
More informationTheoretical physics. Deterministic chaos in classical physics. Martin Scholtz
Theoretical physics Deterministic chaos in classical physics Martin Scholtz scholtzzz@gmail.com Fundamental physical theories and role of classical mechanics. Intuitive characteristics of chaos. Newton
More informationRigid Body Dynamics, SG2150 Solutions to Exam,
KTH Mechanics 011 10 Calculational problems Rigid Body Dynamics, SG150 Solutions to Eam, 011 10 Problem 1: A slender homogeneous rod of mass m and length a can rotate in a vertical plane about a fied smooth
More informationGravitation and Newton s Synthesis
Lecture 10 Chapter 6 Physics I 0.4.014 Gravitation and Newton s Synthesis Course website: http://faculty.uml.edu/andriy_danylov/teaching/physicsi Lecture Capture: http://echo360.uml.edu/danylov013/physics1spring.html
More informationBUMPLESS SWITCHING CONTROLLERS. William A. Wolovich and Alan B. Arehart 1. December 27, Abstract
BUMPLESS SWITCHING CONTROLLERS William A. Wolovich and Alan B. Arehart 1 December 7, 1995 Abstract This paper outlines the design of bumpless switching controllers that can be used to stabilize MIMO plants
More informationChapter #4 EEE8086-EEE8115. Robust and Adaptive Control Systems
Chapter #4 Robust and Adaptive Control Systems Nonlinear Dynamics.... Linear Combination.... Equilibrium points... 3 3. Linearisation... 5 4. Limit cycles... 3 5. Bifurcations... 4 6. Stability... 6 7.
More informationCONTROL DESIGN FOR SET POINT TRACKING
Chapter 5 CONTROL DESIGN FOR SET POINT TRACKING In this chapter, we extend the pole placement, observer-based output feedback design to solve tracking problems. By tracking we mean that the output is commanded
More informationLesson 9. Luis Anchordoqui. Physics 168. Tuesday, October 24, 17
Lesson 9 Physics 168 1 Static Equilibrium 2 Conditions for Equilibrium An object with forces acting on it but that is not moving is said to be in equilibrium 3 Conditions for Equilibrium (cont d) First
More informationMSMS Basilio Bona DAUIN PoliTo
MSMS 214-215 Basilio Bona DAUIN PoliTo Problem 2 The planar system illustrated in Figure 1 consists of a bar B and a wheel W moving (no friction, no sliding) along the bar; the bar can rotate around an
More information