Derivation of the analytical solution for torque free motion of a rigid body
|
|
- Godwin Walton
- 6 years ago
- Views:
Transcription
1 Derivation of the analytical solution for torque free motion of a rigid body Dimitar N. Dimitrov Abstract In this report the rotational motion of a single rigid body in the absence of external forces and torques is considered. The dynamic equations governing this motion are integrated in closed form. Only the general case of an unsymmetrical rigid body is considered. A Matlab implementation is available at []. I. Note We could find two textbooks [], [] that contain a detailed derivation of the analytical solution for torque free motion of a rigid body (the one in [] will be adopted here). The solution involves Jacobi elliptic functions, and a brief introduction is included in Appendix. II. List of variables used All the variables used in this report are expressed in a body fixed coordinate frame. This frame is chosen in such a way that just three scalars (principal moments of inertia) are sufficient to describe the inertial characteristics of the body ([5] pp. 488). J,J and J are the principle moments of inertia. It is assumed that, J > J > J. ω,ω and ω are angular velocities and a hat above them will indicate an initial condition (ˆω ), T is rotational kinetic energy, and h is the length of the angular momentum vector. Note that, T and h are constant due to the assumption of zero external torques. In the derivation below we will make use of the ratio D = h T. III. Derivation The rotational equations of motion (with external torques equal to zero) can be written in the following way: J ω (J J )ω ω = () J ω (J J )ω ω = () J ω (J J )ω ω = () where, a dot over a variable expresses a derivative with respect to time. Important for the discussion below will be the dynamical quantities T and h. The rotational kinetic energy is given by: The magnitude of the angular momentum is given by: T = J ω + J ω + J ω (4) h = J ω + J ω + J ω (5) If J = J J the solution derived here is not valid and one should refer to the case for an axi-symmetric body [] pp. 5
2 The ratio of equations (5) and (4) plays an important role in the following derivation (it is denoted by D). The minimum and maximum values of D can be found by consequently multiplying equation (5) by J and J and then subtracting the resulting equations from (4). J (J J )ω + J (J J )ω = T(J D) (6) J (J J )ω + J (J J )ω = T(D J) (7) The left-hand side expressions are non-negative, hence J D J. D = J and D = J correspond to simple rotation around the major and minor axis, respectively. The solutions for these cases are trivial, and will not be included here, hence it will be assumed that J > D > J. Next, the above equations are solved for ω and ω, respectively as a function of ω : where, a and b are non-negative constants and are defined as: ω = J (J J ) J (J J ) (b ω ) (8) ω = J (J J ) J (J J ) (a ω ) (9) a = T(D J ) J (J J ) b = T(J D) J (J J ) Substitution of (8) and (9) in () leads to: J ω + (J J )s J (J J ) J (J J ) (b ω)s J (J J ) J (J J ) (a ω) = () where, s and s are constants that account for the sign of the two square roots (they can be determined from the signs of the initial conditions ˆω and ˆω. s = sign(ˆω ) and s = sign(ˆω ) ). Rearranging the above equation leads to: dω dt (a ω)(b ω) = s (J J )(J J ) s () J J Equation () has to be integrated in order to obtain ω : dω (a ω)(b ω) = s s (J J )(J J ) J J dt + C () where, C is an integration constant. Note that, the expression integrated on the right-hand side is constant, hence after integration one obtains: dω (a ω)(b ω) = s (J J )(J J ) s t + C () J J sign(ˆω i) denotes the sing of the initial condition ˆω i
3 with t being time. Now (for convenience) we set C = s s (J J )(J J ) J J t, where t is a new constant (to be evaluated shortly from the given initial conditions for ω). The above equation becomes: dω (a ω)(b ω) = s (J J )(J J ) s (t t ) (4) J J The integral on the left-hand side of (4) is an elliptic integral of the first kind (see Appendix ), and in order to solve it, three cases have to be distinguished (the reason for this will become apparent later): a) a < b or D < J b) a > b or D > J c) a = b or D = J Cases a) and b) correspond to polhodes which envelop the major and minor axis, respectively. c) is a limiting case, which separates the polhodes from cases a) and b) (for more information see []). Case a) will be considered first. A. Solution for Case a) Equation (4) can be transformed as follows: or dω a (J J )(J J ) = bs ( ω )( a ω s (t t ) (5) a b ) J J a dx ( x )( k x ) = s s P(t t ) (6) with the following abbreviations (b is included in P): x = ω a k = a b P = T(J D)(J J ) J J J (7) The solution to the elliptic integral (of the first kind) in (6) has the form x = s s sn(p(t t ),k) (8) where, sn is a Jacobi elliptic function, with modulus k (it should satisfy k ) and argument P(t t ). The fact that the modulus k should be between and explains why three cases had to be distinguished. Solving (8) for ω gives: ω = as s sn(p(t t ),k) (9) In (9) the constant t is still undetermined. It can be calculated by setting t = and using the initial condition for ω. ˆω = as s sn( Pt,k) = as s sn(pt,k) () In case a) since a < b the restriction for k is satisfied. In case b) k will be defined as b a. In case c) k =, and the solutions for cases a) and b) converge towards the same result.
4 4 In order to solve the above equation for t, we have to compute an elliptic integral (which is an inverse to the elliptic function sn, see Appendix ) of the form: Pt = ul dω ( ω)( k ω) () where, the upper limit of the integral ul = s s ˆω a. The next step in the derivation is to determine ω and ω. This is made by substituting equation (9) into (8) and (9) (The sign determining terms disappear, since ω is needed). ω = J (J J ) J (J J ) [a T(D J ) J (J J ) sn(p(t t ),k) ] () ω = J (J J ) T(D J ) J (J J ) J (J J ) [ sn(p(t t ),k) ] () and the solution for ω : ω = T(D J ) J (J J ) cn(p(t t ),k) (4) ω = J (J J ) J (J J ) [b T(D J ) J (J J ) sn(p(t t ),k) ] (5) ω = J (J J ) T(J D) J (J J ) J (J J ) J (J J ) J (J J ) T(D J ) J (J J ) sn(p(t t ),k) (6) the above equation is multiplied by J (J J ) T(J D), where it was already assumed that D J. ω J (J J ) T(J D) = k sn(p(t t ),k) (7) ω = T(J D) J (J J ) dn(p(t t ),k) (8) The solutions for ω i (i =,,) are (see equations (9), (4), (8)): ω = s T(D J ) J (J J ) cn(p(t t ),k) (9) ω = s T(D J ) J (J J ) sn(p(t t ),k) ()
5 ω = s T(J D) J (J J ) dn(p(t t ),k) () where, s is a constant that determines the sign of the square root in (). From (9) s can be defined as s = s s. The same relation can be found by substituting (9), () and () back in (). Note that, ω is given by: ω = s P T(D J ) J (J J ) cn(p(t t ),k)dn(p(t t ),k) The relation between s, s and s can be expressed as: s s s = () Altogether four combinations of signs satisfying this relation are possible. This is in accordance with the fact that for each value of the parameter D two separate polhodes exist and each of them passes through zero in two different points (see Table I). 5 TABLE I Possible sign combinations S S S Which row form Table has to be chosen can be determined uniquely from the signs of the initial conditions ˆω and ˆω (if the value of one of them is zero its sign is assumed to be plus ). When s and s are known s = s s. The solution for case a) (D < J ) is complete. Following the same patern, the equations for case b) (D > J ) can be derived as well. Here, we give only the result: B. Solution for Case b) ω = s T(D J ) J (J J ) dn(p(t t ),k) () ω = s T(J D) J (J J ) sn(p(t t ),k) (4) ω = s T(J D) J (J J ) cn(p(t t ),k) (5)
6 6 where, the modulus k = b a, and the sign determining terms s, s and s are the same as in the previous case. C. Solution for Case c) P = T(J J )(D J ) J J J It is not necessary to integrate case c) since both solutions for cases a) and b) converge in the limit D J towards the same result. Note that, here k =, and the elliptic functions become hyperbolic functions. Again s, s and s are the same as in the previous cases. ω = s T(J J ) J (J J ) dn(p(t t ),) = s T(J J ) J (J J ) cosh(p(t t )) (6) T T ω = s sn(p(t t ),) = s tanh(p(t t )) (7) J J ω = s T(J J ) J (J J ) cn(p(t t ),) = s T(J J ) J (J J ) cosh(p(t t )) (8) IV. Angular velocity amplitude As noted in [4], the maximum values for the amplitude of ω, ω and ω is equal to: ω max = T(D J ) J (J J ) ω max = T(J D) J (J J ) ω max T(D J = ) ω max = J (J J ) (D < J ) T(J D) J (J J ) (D > J ) ω max = T J (D = J ) these are precisely the terms before the elliptic functions in the solutions for ω i (i =,,). The values for ω max, ω max and ω max can be found alternatively by solving a maximization problem with constraints (see [4] pp. 6). V. Appendix This appendix contains a short introduction to elliptic functions and integrals. For more detailes see [7], [8].
7 7 A. Elliptic integrals Elliptic integrals occur in many applications [8]. Many of them arose from attempts to solve mechanical problems. For example, the period of a simple pendulum was found to be related to an integral which expressed arc length but no form could be found in terms of simple functions. Another simple example is the integral that arises from the circumference of an ellipse. This example will be considered next. Example: The formula for computing arc length of any curve given by its parametric equations x = f(t), y = g(t), over the range c t d is: d The parametric equations of an ellipse are: c ( dx dt ) + ( dy dt ) dt (9) x = asin(t) y = bcos(t) <= t <= π dx dt = acos(t) dy dt = bsin(t) where a and b are the ellipse radii on the major and minor axis, respectively. The eccentricity of an ellipse is defined as: e = b a < (4) Eccentricity is in case of a circle. Substituting the above defined variables in (9) one obtains: C = π a cos(t) + b sin(t) dt (4) The symmetry of an ellipse can be used in order to rewrite the above equation as: C = 4 π/ Using the relation cos(t) = sin(t) one obtains: or C = 4 π/ C = 4 a cos(t) + b sin(t) dt (4) π/ ( b a )sin(t) dt (4) e sin(t) dt (44)
8 8 The above integral is a complete elliptic integral of the second kind. It is called complete, because the limits of integration are and π/. The eccentricity of the ellipse e is called the modulus of the integral in (44). By setting v = sin(t) (dv = cos(t)dt = t dt) one obtains 4 : C = 4 sin(π/) e v dt (45) v Equation (44) is known as Legendre s form of complete elliptic integral of the second kind, (45) is known as Jacobi s form of this integral (they are equivalent). The solution of the problem introduced in the above example, leads to an elliptic integral ( for more examples see [8]). Three basic types of elliptic integrals exist. One of them was already defined. The other two are: Elliptic integral of the first kind: u = F(φ,k) = φ dt ( k sin(t) ) For generality, we consider φ to be the upper limit of the integral (if φ π the integral is called incomplete). The modulus is denoted by k. Transforming (46) to Jacobi s form one obtains (compare with equation (6)): u = F(φ,k) = Elliptic integral of the third kind: sin(φ) dv ( v )( k v ) (46) (47) φ dt ( + nsin(t) ) ( k sin(t) ) (48) sin(φ) dv ( + nv ) ( v )( k v ) (49) where, n is a constant known as the elliptic characteristic. B. Elliptic functions There are three basic elliptic functions. They arise from the inversion of the elliptic integral of the first kind and are denoted by: sn(u,k), cn(u,k) and dn(u,k) (u is defined in (47)). The modulus k is assumed to be in the range k. It is obvious that: sn(u,) = sin(u) cn(u,) = cos(u) dn(u,) = and sn(u,) = tanh(u) cn(u,) = cosh(u) dn(u,) = cosh(u) 4 Note that, when the integration variable is changed the upper limit of integration should be changed as well.
9 9 Hence, for k = the elliptic functions become circular and for k = they become hyperbolic. The upper limit of integration in (46) (also called Jacobi amplitude) φ is defined in terms of sn(u,k) as: The standard Jacobi elliptic functions satisfy the following identities: Derivatives of the Jacobi elliptic functions are: sin(φ) = sn(u, k) (5) sn(u,k) + cn(u,k) = k sn(u,k) + dn(u,k) = d sn(u,k) = cn(u,k)dn(u,k) du d cn(u,k) = sn(u, k)dn(u, k) du d dn(u,k) = k sn(u,k)cn(u,k) du References [] dir/reports/euler analyt.m [] T. Kane, P. Likins, and D. Levinson, Spacecraft Dynamics, McGraw-Hill, 98 [] Wittenburg, Dynamics of Systems of Rigid Bodies,, 977 [4] P. Hughes, Spacecraft Attitude Dynamics, Wiley, 986 [5] J. Wertz, Spacecraft Attitude Determination and Control, Reidel, 978 [6] F. Rimrott, Introductory Attitude Dynamics, Springer-Verlag, 989 [7] A. Greenhill, The applications of elliptic functions, 89 [8] F. Bowman, Introduction to elliptic functions with applications, 96
1/30. Rigid Body Rotations. Dave Frank
. 1/3 Rigid Body Rotations Dave Frank A Point Particle and Fundamental Quantities z 2/3 m v ω r y x Angular Velocity v = dr dt = ω r Kinetic Energy K = 1 2 mv2 Momentum p = mv Rigid Bodies We treat a rigid
More information[ ] ( ) L = L = r p. We can deal with this more conveniently by writing: so that L becomes: α α α. α α N N N N
Lecture 5: Angular Momentum and Energy of a System The total angular momentum is given by the sum of the angular momenta of all particles in the system: We can deal with this more conveniently by writing:
More informationElliptic functions and Elliptic Integrals
Elliptic functions and Elliptic Integrals R. Herman Nonlinear Pendulum We motivate the need for elliptic integrals by looking for the solution of the nonlinear pendulum equation, θ + ω sin θ. () This models
More information27. Euler s Equations
27 Euler s Equations Michael Fowler Introduction We ve just seen that by specifying the rotational direction and the angular phase of a rotating body using Euler s angles, we can write the Lagrangian in
More informationThe Physical Pendulum
The Physical Pendulum Sourendu Gupta TIFR, Mumbai, India Classical Mechanics 2011 August 20, 2011 The physical pendulum θ L T mg We choose to work with a generalized coordinate q which is the angle of
More informationChapter 14: Vector Calculus
Chapter 14: Vector Calculus Introduction to Vector Functions Section 14.1 Limits, Continuity, Vector Derivatives a. Limit of a Vector Function b. Limit Rules c. Component By Component Limits d. Continuity
More informationChapter 11. Angular Momentum
Chapter 11 Angular Momentum Angular Momentum Angular momentum plays a key role in rotational dynamics. There is a principle of conservation of angular momentum. In analogy to the principle of conservation
More informationElliptic Functions. 1 Introduction. 2 Analysis of motion of a Pendulum. Mark Price. Spring 2001
Elliptic Functions Mark Price Spring 00 Introction I have used the following notation in this essay: the set of all complex numbers is denoted by C and the set of all real numbers is denoted by R. Also
More informationPhys 7221, Fall 2006: Midterm exam
Phys 7221, Fall 2006: Midterm exam October 20, 2006 Problem 1 (40 pts) Consider a spherical pendulum, a mass m attached to a rod of length l, as a constrained system with r = l, as shown in the figure.
More informationRigid bodies - general theory
Rigid bodies - general theory Kinetic Energy: based on FW-26 Consider a system on N particles with all their relative separations fixed: it has 3 translational and 3 rotational degrees of freedom. Motion
More informationElliptical Integrals and Functions
Appendix A Elliptical Integrals and Functions Elliptic integrals and functions are mathematical objects, which nowadays are often omitted in the mathematical curricula of universities. One quite trivial
More informationMotion in Space. MATH 311, Calculus III. J. Robert Buchanan. Fall Department of Mathematics. J. Robert Buchanan Motion in Space
Motion in Space MATH 311, Calculus III J. Robert Buchanan Department of Mathematics Fall 2011 Background Suppose the position vector of a moving object is given by r(t) = f (t), g(t), h(t), Background
More informationLecture 35: The Inertia Tensor
Lecture 35: The Inertia Tensor We found last time that the kinetic energy of a rotating obect was: 1 Trot = ωω i Ii where i, ( I m δ x x x i i, k, i, k So the nine numbers represented by the I i tell us
More informationClassical Mechanics. Luis Anchordoqui
1 Rigid Body Motion Inertia Tensor Rotational Kinetic Energy Principal Axes of Rotation Steiner s Theorem Euler s Equations for a Rigid Body Eulerian Angles Review of Fundamental Equations 2 Rigid body
More informationLecture 9 Kinetics of rigid bodies: Impulse and Momentum
Lecture 9 Kinetics of rigid bodies: Impulse and Momentum Momentum of 2-D Rigid Bodies Recall that in lecture 5, we discussed the use of momentum of particles. Given that a particle has a, and is travelling
More informationCIRCULAR MOTION AND ROTATION
1. UNIFORM CIRCULAR MOTION So far we have learned a great deal about linear motion. This section addresses rotational motion. The simplest kind of rotational motion is an object moving in a perfect circle
More informationLecture 38: Equations of Rigid-Body Motion
Lecture 38: Equations of Rigid-Body Motion It s going to be easiest to find the equations of motion for the object in the body frame i.e., the frame where the axes are principal axes In general, we can
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 information16. Rotational Dynamics
6. Rotational Dynamics A Overview In this unit we will address examples that combine both translational and rotational motion. We will find that we will need both Newton s second law and the rotational
More informationLecture 41: Highlights
Lecture 41: Highlights The goal of this lecture is to remind you of some of the key points that we ve covered this semester Note that this is not the complete set of topics that may appear on the final
More informationPHY 5246: Theoretical Dynamics, Fall Assignment # 10, Solutions. (1.a) N = a. we see that a m ar a = 0 and so N = 0. ω 3 ω 2 = 0 ω 2 + I 1 I 3
PHY 54: Theoretical Dynamics, Fall 015 Assignment # 10, Solutions 1 Graded Problems Problem 1 x 3 a ω First we calculate the moments of inertia: ( ) a I 1 = I = m 4 + b, 1 (1.a) I 3 = ma. b/ α The torque
More informationPLANAR KINETICS OF A RIGID BODY FORCE AND ACCELERATION
PLANAR KINETICS OF A RIGID BODY FORCE AND ACCELERATION I. Moment of Inertia: Since a body has a definite size and shape, an applied nonconcurrent force system may cause the body to both translate and rotate.
More informationStructural Dynamics Lecture 4. Outline of Lecture 4. Multi-Degree-of-Freedom Systems. Formulation of Equations of Motions. Undamped Eigenvibrations.
Outline of Multi-Degree-of-Freedom Systems Formulation of Equations of Motions. Newton s 2 nd Law Applied to Free Masses. D Alembert s Principle. Basic Equations of Motion for Forced Vibrations of Linear
More informationChapter 11. Angular Momentum
Chapter 11 Angular Momentum Angular Momentum Angular momentum plays a key role in rotational dynamics. There is a principle of conservation of angular momentum. In analogy to the principle of conservation
More informationLecture 38: Equations of Rigid-Body Motion
Lecture 38: Equations of Rigid-Body Motion It s going to be easiest to find the equations of motion for the object in the body frame i.e., the frame where the axes are principal axes In general, we can
More informationLecture 9 - Rotational Dynamics
Lecture 9 - Rotational Dynamics A Puzzle... Angular momentum is a 3D vector, and changing its direction produces a torque τ = dl. An important application in our daily lives is that bicycles don t fall
More informationAs a starting point of our derivation of the equations of motion for a rigid body, we employ d Alembert s principle:
MEG6007: Advanced Dynamics -Principles and Computational Methods (Fall, 017) Lecture 10: Equations of Motion for a Rigid Body 10.1 D Alembert s Principle for A Free Rigid Body As a starting point of our
More information14. Rotational Kinematics and Moment of Inertia
14. Rotational Kinematics and Moment of nertia A) Overview n this unit we will introduce rotational motion. n particular, we will introduce the angular kinematic variables that are used to describe the
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 informationLecture 32: Dynamics of Rigid Bodies
Lecture 3: Dynamics of Rigid Bodies Our next topic is the study of a special ind of system of particles a rigid body The defining characteristic of such a system is that the distance between any two particles
More information!T = 2# T = 2! " The velocity and acceleration of the object are found by taking the first and second derivative of the position:
A pendulum swinging back and forth or a mass oscillating on a spring are two examples of (SHM.) SHM occurs any time the position of an object as a function of time can be represented by a sine wave. We
More informationPhysics 8, Fall 2011, equation sheet work in progress
1 year 3.16 10 7 s Physics 8, Fall 2011, equation sheet work in progress circumference of earth 40 10 6 m speed of light c = 2.9979 10 8 m/s mass of proton or neutron 1 amu ( atomic mass unit ) = 1 1.66
More informationSOLUTIONS, PROBLEM SET 11
SOLUTIONS, PROBLEM SET 11 1 In this problem we investigate the Lagrangian formulation of dynamics in a rotating frame. Consider a frame of reference which we will consider to be inertial. Suppose that
More informationThe... of a particle is defined as its change in position in some time interval.
Distance is the. of a path followed by a particle. Distance is a quantity. The... of a particle is defined as its change in position in some time interval. Displacement is a.. quantity. The... of a particle
More informationElliptical Integrals and Functions
Epilogue Everything created by humans has its beginning and its end. At the end of this book, we would like to say that this is just a first step towards a much deeper exposition of the topics discussed
More informationexpression that describes these corrections with the accuracy of the order of 4. frame usually connected with extragalactic objects.
RUSSIAN JOURNAL OF EARTH SCIENCES, English Translation, VOL, NO, DECEMBER 998 Russian Edition: JULY 998 On the eects of the inertia ellipsoid triaxiality in the theory of nutation S. M. Molodensky Joint
More informationReducing round-off errors in rigid body dynamics
Reducing round-off errors in rigid body dynamics Gilles Vilmart INRIA Rennes, ENS Cachan Bretagne, avenue Robert Schuman, 357 Bruz, France Université de Genève, Section de mathématiques, -4 rue du Lièvre,
More informationRotational & Rigid-Body Mechanics. Lectures 3+4
Rotational & Rigid-Body Mechanics Lectures 3+4 Rotational Motion So far: point objects moving through a trajectory. Next: moving actual dimensional objects and rotating them. 2 Circular Motion - Definitions
More informationChapter 9 Notes. x cm =
Chapter 9 Notes Chapter 8 begins the discussion of rigid bodies, a system of particles with fixed relative positions. Previously we have dealt with translation of a particle: if a rigid body does not rotate
More informationPhysics 8, Fall 2013, equation sheet work in progress
(Chapter 1: foundations) 1 year 3.16 10 7 s Physics 8, Fall 2013, equation sheet work in progress circumference of earth 40 10 6 m speed of light c = 2.9979 10 8 m/s mass of proton or neutron 1 amu ( atomic
More informationIf the symmetry axes of a uniform symmetric body coincide with the coordinate axes, the products of inertia (Ixy etc.
Prof. O. B. Wright, Autumn 007 Mechanics Lecture 9 More on rigid bodies, coupled vibrations Principal axes of the inertia tensor If the symmetry axes of a uniform symmetric body coincide with the coordinate
More informationB.Sc. (Semester - 5) Subject: Physics Course: US05CPHY01 Classical Mechanics
1 B.Sc. (Semester - 5) Subject: Physics Course: US05CPHY01 Classical Mechanics Question Bank UNIT: I Multiple choice questions: (1) The gravitational force between two masses is (a) Repulsive (b) Attractive
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 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 informationLecture II: Rigid-Body Physics
Rigid-Body Motion Previously: Point dimensionless objects moving through a trajectory. Today: Objects with dimensions, moving as one piece. 2 Rigid-Body Kinematics Objects as sets of points. Relative distances
More informationA Brief History of Elliptic Integral Addition Theorems
Rose-Hulman Undergraduate Mathematics Journal Volume Issue Article A Brief History of Elliptic Integral Addition Theorems Jose Barrios Montclair State University, barriosj@mail.montclair.edu Follow this
More informationDynamics. Dynamics of mechanical particle and particle systems (many body systems)
Dynamics Dynamics of mechanical particle and particle systems (many body systems) Newton`s first law: If no net force acts on a body, it will move on a straight line at constant velocity or will stay at
More informationSECULAR MOTION IN A 2ND DEGREE AND ORDER-GRAVITY FIELD WITH NO ROTATION. 1. Introduction
SECULAR MOTION IN A 2ND DEGREE AND ORDER-GRAVITY FIELD WITH NO ROTATION D. J. SCHEERES and W. HU Department of Aerospace Engineering, The University of Michigan, Ann Arbor, MI 4819-214, e-mail: scheeres@umich.edu
More informationMAT 132 Midterm 1 Spring 2017
MAT Midterm Spring 7 Name: ID: Problem 5 6 7 8 Total ( pts) ( pts) ( pts) ( pts) ( pts) ( pts) (5 pts) (5 pts) ( pts) Score Instructions: () Fill in your name and Stony Brook ID number at the top of this
More informationThe branch of mechanics which deals with the motion of object is called dynamics. It is divided into two branches:
M th KINEMATICS 4 CHAPTER INTRODUCTION The branch of mechanics which deals with the motion of object is called dynamics. It is divided into two branches: Kinematics: (i) Kinematics (ii) Kinetics The branch
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 informationConcepts in Newtonian mechanics
Concepts in ewtonian mechanics For each question, fill in the blank, circle true or false, circle one (or more) of the multiple choice answers, write the definition, or complete the drawing. The percentages
More information5/6/2018. Rolling Without Slipping. Rolling Without Slipping. QuickCheck 12.10
Rolling Without Slipping Rolling is a combination of rotation and translation. For an object that rolls without slipping, the translation of the center of mass is related to the angular velocity by Slide
More informationForced Mechanical Vibrations
Forced Mechanical Vibrations Today we use methods for solving nonhomogeneous second order linear differential equations to study the behavior of mechanical systems.. Forcing: Transient and Steady State
More informationME 391 Mechanical Engineering Analysis
Solve the following differential equations. ME 9 Mechanical Engineering Analysis Eam # Practice Problems Solutions. e u 5sin() with u(=)= We may rearrange this equation to e u 5sin() and note that it is
More informationAP Physics C Mechanics Objectives
AP Physics C Mechanics Objectives I. KINEMATICS A. Motion in One Dimension 1. The relationships among position, velocity and acceleration a. Given a graph of position vs. time, identify or sketch a graph
More informationChapter 15. Oscillations
Chapter 15 Oscillations 15.1 Simple Harmonic Motion Oscillatory Motion: Motion which is periodic in time; motion that repeats itself in time. Examples: SHM: Power line oscillates when the wind blows past.
More information3-D Kinetics of Rigid Bodies
3-D Kinetics of Rigid Bodies Angular Momentum Generalized Newton s second law for the motion for a 3-D mass system Moment eqn for 3-D motion will be different than that obtained for plane motion Consider
More information(You may need to make a sin / cos-type trigonometric substitution.) Solution.
MTHE 7 Problem Set Solutions. As a reminder, a torus with radii a and b is the surface of revolution of the circle (x b) + z = a in the xz-plane about the z-axis (a and b are positive real numbers, with
More informationPHYS 705: Classical Mechanics. Euler s Equations
1 PHYS 705: Classical Mechanics Euler s Equations 2 Euler s Equations (set up) We have seen how to describe the kinematic properties of a rigid body. Now, we would like to get equations of motion for it.
More informationRotational kinematics
Rotational kinematics Suppose you cut a circle out of a piece of paper and then several pieces of string which are just as long as the radius of the paper circle. If you then begin to lay these pieces
More informationPhysics 312, Winter 2007, Practice Final
Physics 312, Winter 2007, Practice Final Time: Two hours Answer one of Question 1 or Question 2 plus one of Question 3 or Question 4 plus one of Question 5 or Question 6. Each question carries equal weight.
More informationAP PHYSICS 1 Learning Objectives Arranged Topically
AP PHYSICS 1 Learning Objectives Arranged Topically with o Big Ideas o Enduring Understandings o Essential Knowledges o Learning Objectives o Science Practices o Correlation to Knight Textbook Chapters
More informationLecture 37: Principal Axes, Translations, and Eulerian Angles
Lecture 37: Principal Axes, Translations, and Eulerian Angles When Can We Find Principal Axes? We can always write down the cubic equation that one must solve to determine the principal moments But if
More information3 Rigid Spacecraft Attitude Control
1 3 Rigid Spacecraft Attitude Control Consider a rigid spacecraft with body-fixed frame F b with origin O at the mass centre. Let ω denote the angular velocity of F b with respect to an inertial frame
More informationMath 53: Worksheet 9 Solutions
Math 5: Worksheet 9 Solutions November 1 1 Find the work done by the force F xy, x + y along (a) the curve y x from ( 1, 1) to (, 9) We first parametrize the given curve by r(t) t, t with 1 t Also note
More information1. The first thing you need to find is the mass of piece three. In order to find it you need to realize that the masses of the three pieces must be
1. The first thing you need to find is the mass of piece three. In order to find it you need to realize that the masses of the three pieces must be equal to the initial mass of the starting rocket. Now
More informationVariation of Parameters
Variation of Parameters James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University April 13, 218 Outline Variation of Parameters Example One We eventually
More informationPhysics 106a, Caltech 4 December, Lecture 18: Examples on Rigid Body Dynamics. Rotating rectangle. Heavy symmetric top
Physics 106a, Caltech 4 December, 2018 Lecture 18: Examples on Rigid Body Dynamics I go through a number of examples illustrating the methods of solving rigid body dynamics. In most cases, the problem
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 informationClassical Mechanics III (8.09) Fall 2014 Assignment 3
Classical Mechanics III (8.09) Fall 2014 Assignment 3 Massachusetts Institute of Technology Physics Department Due September 29, 2014 September 22, 2014 6:00pm Announcements This week we continue our discussion
More informationPHYSICS 210 SOLUTION OF THE NONLINEAR PENDULUM EQUATION USING FDAS
PHYSICS 210 SOLUTION OF THE NONLINEAR PENDULUM EQUATION USING FDAS 1. PHYSICAL & MATHEMATICAL FORMULATION O θ L r T m W 1 1.1 Derivation of the equation of motion O Consider idealized pendulum: Mass of
More informationPhysics 141 Rotational Motion 1 Page 1. Rotational Motion 1. We're going to turn this team around 360 degrees.! Jason Kidd
Physics 141 Rotational Motion 1 Page 1 Rotational Motion 1 We're going to turn this team around 360 degrees.! Jason Kidd Rigid bodies To a good approximation, a solid object behaves like a perfectly rigid
More informationM. van Berkel DCT
Explicit solution of the ODEs describing the 3 DOF Control Moment Gyroscope M. van Berkel DCT 28.14 Traineeship report Coach(es): Supervisor: dr. N. Sakamoto Prof. dr. H. Nijmeijer Technische Universiteit
More informationPHYSICS 110A : CLASSICAL MECHANICS
PHYSICS 110A : CLASSICAL MECHANICS 1. Introduction to Dynamics motion of a mechanical system equations of motion : Newton s second law ordinary differential equations (ODEs) dynamical systems simple 2.
More informationHandout 6: Rotational motion and moment of inertia. Angular velocity and angular acceleration
1 Handout 6: Rotational motion and moment of inertia Angular velocity and angular acceleration In Figure 1, a particle b is rotating about an axis along a circular path with radius r. The radius sweeps
More informationOrbital Mechanics! Space System Design, MAE 342, Princeton University! Robert Stengel
Orbital Mechanics Space System Design, MAE 342, Princeton University Robert Stengel Conic section orbits Equations of motion Momentum and energy Kepler s Equation Position and velocity in orbit Copyright
More informationSections minutes. 5 to 10 problems, similar to homework problems. No calculators, no notes, no books, no phones. No green book needed.
MTH 34 Review for Exam 4 ections 16.1-16.8. 5 minutes. 5 to 1 problems, similar to homework problems. No calculators, no notes, no books, no phones. No green book needed. Review for Exam 4 (16.1) Line
More informationRigid Body Equations of Motion for Modeling and Control of Spacecraft Formations - Part 1: Absolute Equations of Motion.
Rigid ody Equations of Motion for Modeling and Control of Spacecraft Formations - Part 1: Absolute Equations of Motion. Scott R. Ploen, Fred Y. Hadaegh, and Daniel P. Scharf Jet Propulsion Laboratory California
More informationAPPLICATIONS ELLIPTIC FUNCTIONS CLASSICAL AND ALGEBRAIC GEOMETRY. Jamie Snape. Collingwood College, University of Durham
APPLICATIONS OF ELLIPTIC FUNCTIONS IN CLASSICAL AND ALGEBRAIC GEOMETRY Jamie Snape Collingwood College, University of Durham Dissertation submitted for the degree of Master of Mathematics at the University
More informationRotational Motion. Chapter 4. P. J. Grandinetti. Sep. 1, Chem P. J. Grandinetti (Chem. 4300) Rotational Motion Sep.
Rotational Motion Chapter 4 P. J. Grandinetti Chem. 4300 Sep. 1, 2017 P. J. Grandinetti (Chem. 4300) Rotational Motion Sep. 1, 2017 1 / 76 Angular Momentum The angular momentum of a particle with respect
More informationRotation. Rotational Variables
Rotation Rigid Bodies Rotation variables Constant angular acceleration Rotational KE Rotational Inertia Rotational Variables Rotation of a rigid body About a fixed rotation axis. Rigid Body an object that
More informationRotational Kinematics and Dynamics. UCVTS AIT Physics
Rotational Kinematics and Dynamics UCVTS AIT Physics Angular Position Axis of rotation is the center of the disc Choose a fixed reference line Point P is at a fixed distance r from the origin Angular Position,
More informationAPPLIED MATHEMATICS ADVANCED LEVEL
APPLIED MATHEMATICS ADVANCED LEVEL INTRODUCTION This syllabus serves to examine candidates knowledge and skills in introductory mathematical and statistical methods, and their applications. For applications
More informationMath 302 Outcome Statements Winter 2013
Math 302 Outcome Statements Winter 2013 1 Rectangular Space Coordinates; Vectors in the Three-Dimensional Space (a) Cartesian coordinates of a point (b) sphere (c) symmetry about a point, a line, and a
More informationLesson 8. Luis Anchordoqui. Physics 168. Thursday, October 11, 18
Lesson 8 Physics 168 1 Rolling 2 Intuitive Question Why is it that when a body is rolling on a plane without slipping the point of contact with the plane does not move? A simple answer to this question
More informationExercise: concepts from chapter 3
Reading:, Ch 3 1) The natural representation of a curve, c = c(s), satisfies the condition dc/ds = 1, where s is the natural parameter for the curve. a) Describe in words and a sketch what this condition
More informationGeneral Definition of Torque, final. Lever Arm. General Definition of Torque 7/29/2010. Units of Chapter 10
Units of Chapter 10 Determining Moments of Inertia Rotational Kinetic Energy Rotational Plus Translational Motion; Rolling Why Does a Rolling Sphere Slow Down? General Definition of Torque, final Taking
More information3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1
Math Problem a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1 3 6 Solve the initial value problem u ( t) = Au( t) with u (0) =. 3 1 u 1 =, u 1 3 = b- True or false and why 1. if A is
More informationHyperbolic-Type Orbits in the Schwarzschild Metric
Hyperbolic-Type Orbits in the Schwarzschild Metric F.T. Hioe* and David Kuebel Department of Physics, St. John Fisher College, Rochester, NY 468 and Department of Physics & Astronomy, University of Rochester,
More informationNORMAL MODES, WAVE MOTION AND THE WAVE EQUATION. Professor G.G.Ross. Oxford University Hilary Term 2009
NORMAL MODES, WAVE MOTION AND THE WAVE EQUATION Professor G.G.Ross Oxford University Hilary Term 009 This course of twelve lectures covers material for the paper CP4: Differential Equations, Waves and
More informationAstromechanics. 10. The Kepler Problem
Astromechanics 10. The Kepler Problem One of the fundamental problems in astromechanics is the Kepler problem The Kepler problem is stated as follows: Given the current position a velocity vectors and
More informationVariational Discretization of Euler s Elastica Problem
Typeset with jpsj.cls Full Paper Variational Discretization of Euler s Elastica Problem Kiyoshi Sogo Department of Physics, School of Science, Kitasato University, Kanagawa 8-8555, Japan A discrete
More informationVideo 2.1a Vijay Kumar and Ani Hsieh
Video 2.1a Vijay Kumar and Ani Hsieh Robo3x-1.3 1 Introduction to Lagrangian Mechanics Vijay Kumar and Ani Hsieh University of Pennsylvania Robo3x-1.3 2 Analytical Mechanics Aristotle Galileo Bernoulli
More information2007 Problem Topic Comment 1 Kinematics Position-time equation Kinematics 7 2 Kinematics Velocity-time graph Dynamics 6 3 Kinematics Average velocity
2007 Problem Topic Comment 1 Kinematics Position-time equation Kinematics 7 2 Kinematics Velocity-time graph Dynamics 6 3 Kinematics Average velocity Energy 7 4 Kinematics Free fall Collisions 3 5 Dynamics
More informationPhys 7221 Homework # 8
Phys 71 Homework # 8 Gabriela González November 15, 6 Derivation 5-6: Torque free symmetric top In a torque free, symmetric top, with I x = I y = I, the angular velocity vector ω in body coordinates with
More informationInternal Torques and Forces in Gyrostats with Magnetically Suspended Rotors
Internal Torques and Forces in Gyrostats with Magnetically Suspended Rotors by Marcus Carl Pressl Thesis submitted to the Faculty of the Virginia Polytechnic Institute and State University in partial fulfillment
More informationNoether s Theorem. 4.1 Ignorable Coordinates
4 Noether s Theorem 4.1 Ignorable Coordinates A central recurring theme in mathematical physics is the connection between symmetries and conservation laws, in particular the connection between the symmetries
More informationSTATICS Chapter 1 Introductory Concepts
Contents Preface to Adapted Edition... (v) Preface to Third Edition... (vii) List of Symbols and Abbreviations... (xi) PART - I STATICS Chapter 1 Introductory Concepts 1-1 Scope of Mechanics... 1 1-2 Preview
More informationANALYTICAL MECHANICS. LOUIS N. HAND and JANET D. FINCH CAMBRIDGE UNIVERSITY PRESS
ANALYTICAL MECHANICS LOUIS N. HAND and JANET D. FINCH CAMBRIDGE UNIVERSITY PRESS Preface xi 1 LAGRANGIAN MECHANICS l 1.1 Example and Review of Newton's Mechanics: A Block Sliding on an Inclined Plane 1
More information