expression that describes these corrections with the accuracy of the order of 4. frame usually connected with extragalactic objects.
|
|
- Vivian Merritt
- 6 years ago
- Views:
Transcription
1 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 Institute of the Physics of the Earth, Bolshaja Gruzinskaja, Moscow 38 Abstract. Modern nutation theories of the Earth planets use the assumption, that the planet s inertia ellipsoid is symmetrical about the axes of rotation. Below we estimate the eects of the inertia ellipsoid triaxiality accurate to ((B ; A)=(C ; A)) 4, A B C are equatorial polar principal moments of inertia, respectively.. Introduction As known, the nutational motion of the Earth planets in space is fully determined by Euler angles # E E E, which describe the position of the body- xed (Tesser's) axes (x y z) with respect to an immobile reference frame. At the same time, dynamic equations of motion of an actual planet model dene only the components (! x,! y,! z ) of the angular velocity vector! with respect to a moving reference frame. In the case of an axially symmetrical planet (equatorial moments of inertia satisfy the condition A = B) the Euler angles are connected with these components through well-known Poinsot's formulae (see, for example, [Lau Lifshitz, 964]), but in the case A 6= B Euler's relations between # E E E (! x! y! z ) present a rather complicated (nonlinear) system of ordinary differential equations which has not been investigated in all details as yet. As was shown by Zharkov et al. [996], in the case of small violation of the axial symmetry (if the parameter =(A ; B)=(C ; A) satises the condition j j) Euler's equations can be solved by using the perturbation method in powers of this parameter in such an approach, rst-order (in ) corrections describe only short-period (semidiurnal) perturbations of the values # E, E, E. Nevertheless, this result cannot be considered as a nal solution of the problem under consideration, because second-order corrections of the order of can contain also longperiod terms, which perturb not only instant values of the nutational amplitudes, but also their mean values. c998 Russian Journal of Earth Sciences. Paper No. TJE983. Online version of this paper was published on July 998. URL: Modern VLBI measurements are most sensible just to these terms, their analysis is most interesting for any practical purposes. Below we present a simple analytical expression that describes these corrections with the accuracy of the order of 4.. Statement of the Problem As mentioned above, modern astrometric measurements determine the motion of xed points of a planetary surface R with respect to an immobile reference frame usually connected with extragalactic objects. At the same time, well-known Liouville's equation M M = L () (! is the angular velocity of the Tesser's reference frame (x y z), M is the angular momentum of the planet, L is the torque of external forces, the dot above symbol denotes the time derivative) determines only the motion of the vector with respect to the moving system of coordinates (x y z). To express the components of the vector R in terms of the known values (! x! y! z ), it is convenient to use the well-known Poinsot's kinematic relations (see, for example, [Zharkov et al., 996]) which connect the trajectory of the vector! with respect to the moving reference frame (x y z) (polhode) the trajectory of the same vector in space (herpolhode). In the case of an axially symmetrical rigid planet ( = the products of inertia I ik vanish), the components M! are connected by the well-known relation M x = A! x M y = B! y M z = C! z () the torque L may be represented in the form L = l(i sin(t) +j cos(t)) (3) 5
2 5 molodensky: effects of the inertia ellipsoid triaxiality i, j are unit vectors oriented along the axes x y, respectively, l =(C ; A)v t =a (4) v t is the amplitude of the nearly diurnal tidal potential, V i = v t a (xz cos(t)+yz sin(t)) = v t a r sin # cos # cos(t ; ) (5) a is the mean radius of the planet, (r # ) are spherical coordinates, is tidal frequency. In the most general case, Liouville equation () represents a system of three nonlinear ordinary dierential equations with respect to three unknown functions! x (t)! y (t)! z (t). In the case =, it has a very simple exact solution. In fact, substituting ()-(4) into (), we obtain! =!(i" cos(t) +j" sin(t) +k) (6) k is a unit vector that coincides with the z axis " = v t(a ; C)(A ;!(A ; C)) a A! +(A ; C)! (7) is the amplitude of nutational motion with respect to the moving reference frame (x y z). The fact that the amplitudes the x y components of vector! in (6) are equivalent is a direct consequence of the axial symmetry of the problem under consideration. In view of this, the surfaces of polhode herpolhode are exactly conical in accordance with Poinsot's kinematic relations, the ratio of apex angles of these cones is dened by the relation = arcsin(;"=( )) (8) the amplitude of nutational motion of a vector R in the body-xed reference frame (x y z) is equal to + " in the case of prograde nutational motion (if the directions of the! motion nutation coincide) to in the opposite case of retrograde motion. At 6= the solutions of dynamic kinematic nonlinear equations can be rather complicated. Because long-period perturbations are most important in practical applications, below we use the method of its approximate solution based on the following procedure:. Dynamic equations () provide expressions for the components! x (t),! y (t),! z (t) of the vector! with respect to the moving reference frame, in which the eects of the inertia ellipsoid triaxiality are taken into account. Inverting Euler's kinematic relations between the known components! x (t)! y (t)! z (t) of the vector! in the moving reference frame, the components! (t),! (t),! 3 (t) of the same vector in space, Euler's angles # E E E, these angles can be expressed in terms of! x (t)! y (t)! z (t). Taking into account, that the kinematic relations are nonlinear, we use the rst-order method of perturbations with respect to the small parameter. 3. To connect short-period rst-order perturbations of the functions! (t)! (t)! 3 (t) long-period second-order perturbations of the same functions, we use the well-known Poinsot's theorem stating that the rotation of an arbitrary triaxial rigid body may be represented as the rolling of the polhode over the herpolhode without sliding. In accordance with this theorem, the ratio of trajectory lengths of the vector! in the moving xed reference frames L L is equal to the ratio of the corresponding periods. The perturbation of the trajectory L in turn consists of () the known terms which describe semidiurnal oscillations () unknown terms of higher orders of smallness which describe the mean radius of this trajectory. As is shown below, this permits the mean radius of the trajectory L to be found with sucient accuracy (of the order of 4 ). 3. Tidal Torque for a Triaxial Body Using expression (5) for the components of tidal potential in Cartesian coordinates, it is easy to nd the tidal torques which act on an axially nonsymmetrical planet L = Z r rv t d = = t [i(c ; B) sin(t) +j(a ; C) cos(t)] a (9) is the total volume of the Earth. Replacing the relation between! y M y in equations () by the formula M y = B! y then substituting () (9) in Liouville's equation (), we obtain A _! x +(C ; B)! y! z = t (C ; B) sin(t) a (a) B _! y +(A ; C)! x! z = t a (A ; C) cos(t) C _! z +(B ; A)! x! y =: (b) (c) Using relations (6) (7), it is easy to see that the
3 molodensky: effects of the inertia ellipsoid triaxiality 53 ratio of terms C _! z (B ; A)! x! y in (c) is of the order of ev t =(ga), e is the attening of the planet g is the gravity on its surface. This ratio is negligibly small for all planets of the solar system, the approximation! z =! = const is valid with a very high accuracy. In this approximation, the solution of equations (a) (b) describing forced nutation precession is # (if this amplitude satises the condition j# j) 4. First-order Approximation for a Slightly Triaxial Body (3d)! =!(i" cos(t) +j" sin(t) +k) (a) " = v t (B ; C)(B ;!(A ; C)) a AB! +(A ; C)(B ; C)! 3 (b) To estimate the eects of triaxiality of the planet in the rst approximation, we rst perform some simple transformations. Multiplying (a) (b) by cos E sin E then subtracting the results, we express # E in terms of " " as follows: _# E =!(" cos E ; " sin E ): (4a) " = v t (A ; C)(A ;!(B ; C)) a AB! +(A ; C)(B ; C)! 3 (c) Analogously, multiplying (a) (b) by sin E cos E then summing the results, we have Substituting relations () in the well-known Euler relations between Euler angles (# E E E ) (! x! y! z )[Lau Lifshitz, 964]! x = _ E sin # E sin E + _ # E cos E (a)! y = _ E sin # E cos E + _ # E sin E (b)! z = _ E cos # E cos E + _ E (c) we obtain a system of three nonlinear ordinary dierential equations with respect to three unknown functions # E E E. In the case of an axially symmetrical body (A = B " = " = ", " is dened by (7)) these equations have an elementary (Poinsot) solution in the form _ E sin # E =!(" sin E + " cos E ) (4b) To linearize these equations, their solutions are represented by the superposition E = + E = + # E = # + # (5) # are the solutions of axially symmetrical problem (3a)-(3c), # are small corrections due to the triaxiality eects of the inertia ellipsoid. Substituting (5) (a) into (4a), (4b), (c) neglecting the terms of the order of,we obtain _ sin # + _ sin # + # _ cos # =!" cos(t)(sin + cos )+!" sin(t)(cos ; sin ) (6a) # = arctan (3a) _# + _ # =!" cos(t)(cos ; sin );!" sin(t)(sin + cos ) (6b) = ; t (3b) = sin # t + const (3c) which show, that the angular frequency amplitude of nutational motion in space are _ = sin # _ cos # + _ cos # ; # _ sin # + + _ + _ =: (6c) Transforming, in accordance with (3b), the main terms in the right-h sides of these equations,! " + "!" cos(t) sin " sin(t) cos = sin( + t) " ; "! " + " " ; " sin( ; t) = cos(t)!" cos(t)cos ;!" sin(t) sin =
4 54 molodensky: effects of the inertia ellipsoid triaxiality! " + " cos( + t) " ; "! " ; " sin(t) cos( ; t) = taking into account, that the zero-order values of # are connected with "! by relations (3), we obtain, after dierentiation of these relations with respect to the time, Since the system of equations (9) is of the third order, it has three linearly independent homogeneous solutions. It is easy to show, however, that these solutions are not interesting for the problem under consideration, because they describe the small variations in Euler angles caused by an arbitrary innitesimal rotation of the xed reference frame. The inhomogeneous solution of (9) has the form _ sin # =!" =! " + " (7a) = ; = a sin(t) # = b cos(t) () _# = (7b) _ cos # + _ =: (7c) Using these relations, equations (6) can be rewritten as follows: _ sin # + # _ cos # ;! " + " =! " ; " cos(t) sin(t) = (8a) a = ; (" ; " )(3 )( )) "(! ; 3 ) b = ;a : _# +! " + " =! " ; " _ cos # ; # _ sin # + _ : sin(t) (8b) (8c) Eliminating # _ from these relations with the help of (3a) (7c), we nally obtain the system of three inhomogeneous linear ordinary dierential equations with the coecients c _ + c # = cos(t) _# + c 3 = sin(t) _ + c _ 4 ; c 3 # =: c =sin# = (() +( ) ) = c =!" cot # = c 4 =cos# = c 3 =!" (() +( ) ) = =!(" ; " )=: (9) These expressions fully determine rst-order shortperiod perturbations of the nutational precessional motion caused by the triaxiality of the planet s gure. To nd short-period motion perturbations of the vector!, we introduce a xed Cartesian reference frame (e e e 3 )insuchaway that the axis e 3 coincides with the axis of the unperturbed cone of herpolhode, E is equal to the angle between the nodal line vector e. Then the relations between components (!!! 3 ) Euler angles (# E E E ) are completely analogous to (a)-(c) after replacement of E by E E by E :! = _ E sin # E sin E + _ # E cos E (a)! = _ E sin # E cos E ; _ # E sin E (b)! 3 = _ E cos # E sin E + _ E : (c) In the case of nutational motion of any actual planet, its spatial amplitude =( ) does not exceed a few seconds of arc, we can use, instead of (3a), a more simple expression # = : () Substituting (5), (3b), (3c), (), (7) into (), we have
5 molodensky: effects of the inertia ellipsoid triaxiality 55! =(; ; a cos(t))( + +b cos(t)) sin(( )t + a sin(t)); ;b sin(t) cos(( )t + a sin(t))! =(; ; a cos(t))( + +b cos(t)) cos(( )t + a sin(t))+ +b sin(t) cos(( )t + a sin(t)) (3a) (3b) part of this function f (s) (B;A) which is connected with short-period (semidiurnal) perturbations of herpolhode. Using this value, the total perturbation of herpolhode length can be represented as a sum of the known shortperiod perturbation L () f (s) (B ;A) unknown longperiod perturbation L () f () (B ;A) that describes vari- ation in the \mean" herpolhode radius: = L () (s) (f L () f (B ; A) = () (B ; A)+f (B ; A)): (5a)! 3 =!: (3c) Like the rst-order corrections to Euler angles (), these expressions contain only short-period (semidiurnal) perturbations. 5. Long-period terms in the second-order approximation Now, we will use the well-known Poinsot's theorem stating that the rolling of polhode over herpolhode takes place without sliding. On the strength of this theorem, a general method for the calculation of long-period herpolhode perturbations may be formulated as follows. Let the total length of the vector! trajectory in the moving reference frame (x y z) is dened by the expression L = L () ( + f (B ; A)) (4a) L () is the length of the vector in the case of an axially symmetrical planet (A = B) f (B ; A) isa small correction due to +eects of triaxiality. Because this trajectory coincides with circle (6) if A = B with ellipse (a) if A 6= B, the function f actually describes the known dierence between the length of ellipse (a) with the semiaxes!"!" the length!" of circle (6). Likewise, using relations (3), the total length of the! trajectory in the xed reference frame (e e e 3 )can be represented in the form L = L () ( + f (B ; A)) (4b) L () = j=( )jl () is the total length of the herpolhode in the biaxial case f (B ; A) describes the dierence between the lengths in the triaxial biaxial cases. As mentioned before, the absence of exact expressions for the Euler angles in terms of known components! x (t)! y (t) precludes direct determination of this function. However, expressions (3) readily provide the In accordance with Poinsot's theorem, the condition of the polhode rolling over herpolhode without sliding is that the trajectories described by the angular velocity vector during equal time intervals have the same lengths in moving xed reference frames: L =T = L () ( + f (B ; A))=T = = L () ( + f (B ; A))=T = = j=( )jl () (s) ( + f (B ; A)+ +f () (B ; A))=T T T are full revolution periods of the vector! in the moving xed reference frames, respectively. Since these values are connected through the relation T = j=( )jt the unknown value f () (B ; A) can be expressed in terms of the known functions f (B ;A) f (s) (B ;A) in the following, very simple form: f () (B ; A) =f (B ; A) ; f (s) (B ; A) (5b) Then, variation in the mean radius of herpolhode (which coincides with variation in the nutation amplitude A multiplied by!) may be represented in the form A = A f () (B ; A) = (6) = A (f (B ; A) ; f (s) (B ; A)) A is the unperturbed nutation amplitude corresponding to the axially symmetrical approximation. In accordance with (), the complete trajectory length of the vector! in the moving reference frame (x y z) is L =! =! Z Z (" cos + " sin ) = d = ( " + " + " ; " cos()) = d:
6 56 molodensky: effects of the inertia ellipsoid triaxiality Developing the integr function as the Taylor s series in small parameter =(" ; + ) with " )=(" " the help of (b) (c) neglecting terms of the order of 4,wehave Z _! = k cos( + )+k sin( + ) (3a) _! = k sin( + )+k cos( + ) (3b) L =! (( + " " ) = ; ; (" ; " ) 8 p (" + " )3= cos ())d = =!(" + (" ; " ) 8" ; (" ; " ) 64" 3 )= =!(" + (" ; " ) ): 6" (7a) Comparing this expression relation (4a) yields f (A ; B) = (" ; " ) L = lim N! N 6" (7b) the nutation amplitudes " " " are dened by expressions (b), (c), (7), respectively. Strictly speaking, the denition of the complete trajectory length of the vector! in the xed reference frame (e e e 3 needs some renement, because this trajectory is generally nonclosed. In fact, it is obvious that, since the period of short-period perturbations is equal to = the period of nutational motion is equal to =j j, this trajectory is closed only if = k( ), k is any integer number. Of course, nutation frequencies of actual planets do not satisfy this condition, the herpolhode is not a closed curve. In order to apply Poinsot's theorem to the case of nonclosed trajectory of the angular velocity vector in space, we dene the \mean length" of this trajectory as follows: Z N ) ) ) = dt: (8) Substituting (3) into (8) exping, as before, the integr function in Taylor's series with the accuracy of the order of 4,we obtain lim N! N L = j Z j N ( + ; )dt (9a) _! 3 = (3c) k = ; ; (b( ))+ +a ( cos(t); ;(ab( ) +4a ) cos (t) k =(4a +b( )) sin(t) = = ;a"! sin(t): (3a) (3b) Substituting these expressions into (9), (3a), (3b) neglecting the terms of the fourth higher orders, we have (_! ) +(_! ) +(_! 3 ) = k + k = =() ( ; a! ;! + cos(t)+ +a (! ; ) ; 4( ) ( ) cos (t)+ +4a! ( ) sin (t)) j j lim N! N Z N ( + )dt = L () ( + a 5! ; 3 ; 6! 4( ) ) j j lim N! N Z N 8 dt = (3a) (3b) (3c) (! ; L () ) (a 4( ) ) Z N 3 dt =: (33) = () [( _! ) +(_! ) +(_! 3 ) ] ; : Substituting (a-c) into (6a, b), we have (9b) Summing these expressions, we nally obtain L = j=( )jl () (s) ( + f (B ; A)) (34a)
7 molodensky: effects of the inertia ellipsoid triaxiality 57 f (s) (B ; A) =a! ;! ; ( ) = = (" ; " ) (3 ) (! ;! ; ) 4" (! ; 3 ) = = (" ; " ) (! ;! ; ) 4" ( ;!) (34b) Substituting (34b) (6b) into (3), we obtain the nal value of the correction to the mean nutation amplitudes A = A + A = A ( + f () (B ; A)) = = A ( + (" ; " ) 4" ( 4 ; (! ;! ; ) ( ;!) )) = = A ( + (" ; " ) 6" ( )(5 ; 3!) ( ;!) ): (35) respectively, T ==! is one sidereal day). Nutational motion with the period T s = T = takes place only in the case of synchronised diurnal orbital rotation (if at the points of equinoxes the orientation of the planet's principal equatorial moments of inertia is invariable in space). In this case, the nutation amplitudes are obviously determined not by the mean values of A, B averaged over the total period, but by their \eective" values in the vicinity ofpoints at which the tidal torque is maximum. In this specic case, the corrections to the mean nutation amplitudes are proportional to (" " )=" rather than (" " ) =". This qualitative consideration explains, why second-order correction (35) tends to in- nity in the limiting case!!. Acknowledgement. This work was supported by the Russian Foundation for Basic Research, project No Note that this correction vanishes if = ;! or = 3!=5. The rst root corresponds to the precession frequency. Vanishing of this expression for the precession frequency implies that the precessional constant depends only on the value ; (A + B)= but not on the dierence of equatorial moments of inertia B ; A. The second root, as well as the denominator root = ;!, corresponds to the case of high nutation frequencies in space as compared with the period of diurnal rotation (these frequencies correspond to the spatial periods of nutational motion T s =5T =8T s = T =, References Lau, L. D. E. M. Lifshitz, Mechanics, Moscow, Nauka, 964 (in Russian). Molodensky, S.M.,Tides, Nutation, Internal Structure of the Earth, Moscow, Nauka, 98, 4 pp. (in Russian). Zharkov V. N., S. M. Molodensky, E. Groten, A. Brzezinski, P. Varga. The Earth its rotation. Low-frequency geodynamics, Heidelberg, Wichman Verlag, 996, 53 pp. (Received May )
13. Rigid Body Dynamics II
University of Rhode Island DigitalCommons@URI Classical Dynamics Physics Course Materials 2015 13. Rigid Body Dynamics II Gerhard Müller University of Rhode Island, gmuller@uri.edu Creative Commons License
More informationProblem 1. Mathematics of rotations
Problem 1. Mathematics of rotations (a) Show by algebraic means (i.e. no pictures) that the relationship between ω and is: φ, ψ, θ Feel free to use computer algebra. ω X = φ sin θ sin ψ + θ cos ψ (1) ω
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 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 informationPerturbed Earth rotation. L. D. Akulenko, S. A. Kumakshev, and A. M. Shmatkov
Perturbed Earth rotation L. D. Akulenko, S. A. Kumakshev, and A. M. Shmatkov 1 Introduction Numerous astrometric studies are based on the dynamic theory of the Earth s rotation with respect to the center
More information2 so that the (time-dependent) torque, as measured in space-xed coordinates, is ~N = ma2! 2 (^x cos!t +^y sin!t) 10 = a p 5 2 ^z ~ F 2 ~F 1 = a p 5^z
1 PHYS 321 Homework Assignment #9 Due: Friday, 22 November 2002 (4 probs) 1. Problem 7-14, B&O p. 280. The situation is as shown above. We take the body-xed axes ^x 0 and ^y 0 to lie in the plane of the
More informationPhysical Dynamics (SPA5304) Lecture Plan 2018
Physical Dynamics (SPA5304) Lecture Plan 2018 The numbers on the left margin are approximate lecture numbers. Items in gray are not covered this year 1 Advanced Review of Newtonian Mechanics 1.1 One Particle
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 informationTorque and Rotation Lecture 7
Torque and Rotation Lecture 7 ˆ In this lecture we finally move beyond a simple particle in our mechanical analysis of motion. ˆ Now we consider the so-called rigid body. Essentially, a particle with extension
More informationRevision of the physical backgrounds of Earth s rotation
Revision of the physical backgrounds of Earth s rotation L. Völgyesi Budapest University of Technology and Economics, H-151 Budapest, Müegyetem rkp. 3. Hungary e-mail: lvolgyesi@epito.bme.hu bstract. Rotation
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 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 informationPhysics for Scientists and Engineers 4th Edition, 2017
A Correlation of Physics for Scientists and Engineers 4th Edition, 2017 To the AP Physics C: Mechanics Course Descriptions AP is a trademark registered and/or owned by the College Board, which was not
More information06. Lagrangian Mechanics II
University of Rhode Island DigitalCommons@URI Classical Dynamics Physics Course Materials 2015 06. Lagrangian Mechanics II Gerhard Müller University of Rhode Island, gmuller@uri.edu Creative Commons License
More informationSatellite Communications
Satellite Communications Lecture (3) Chapter 2.1 1 Gravitational Force Newton s 2nd Law: r r F = m a Newton s Law Of Universal Gravitation (assuming point masses or spheres): Putting these together: r
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 informationProblem Solving Session 11 Three Dimensional Rotation and Gyroscopes Solutions
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Physics Physics 8.01 Problem Solving Session 11 Three Dimensional Rotation and Gyroscopes Solutions W14D3-1 Rotating Skew Rod Solution Consider a simple
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 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 informationOur goal is to solve a general constant coecient linear second order. this way but that will not always happen). Once we have y 1, it will always
October 5 Relevant reading: Section 2.1, 2.2, 2.3 and 2.4 Our goal is to solve a general constant coecient linear second order ODE a d2 y dt + bdy + cy = g (t) 2 dt where a, b, c are constants and a 0.
More informationMotion under the Influence of a Central Force
Copyright 004 5 Motion under the Influence of a Central Force The fundamental forces of nature depend only on the distance from the source. All the complex interactions that occur in the real world arise
More informationSTATICS & DYNAMICS. Engineering Mechanics. Gary L. Gray. Francesco Costanzo. Michael E. Plesha. University of Wisconsin-Madison
Engineering Mechanics STATICS & DYNAMICS SECOND EDITION Francesco Costanzo Department of Engineering Science and Mechanics Penn State University Michael E. Plesha Department of Engineering Physics University
More informationFINAL EXAM GROUND RULES
PHYSICS 507 Fall 2011 FINAL EXAM Room: ARC-108 Time: Wednesday, December 21, 10am-1pm GROUND RULES There are four problems based on the above-listed material. Closed book Closed notes Partial credit will
More informationLECTURE 15 + C+F. = A 11 x 1x1 +2A 12 x 1x2 + A 22 x 2x2 + B 1 x 1 + B 2 x 2. xi y 2 = ~y 2 (x 1 ;x 2 ) x 2 = ~x 2 (y 1 ;y 2 1
LECTURE 5 Characteristics and the Classication of Second Order Linear PDEs Let us now consider the case of a general second order linear PDE in two variables; (5.) where (5.) 0 P i;j A ij xix j + P i,
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 information7. The gyroscope. 7.1 Introduction. 7.2 Theory. a) The gyroscope
K 7. The gyroscope 7.1 Introduction This experiment concerns a special type of motion of a gyroscope, called precession. From the angular frequency of the precession, the moment of inertia of the spinning
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 information*School of Aeronautics and Astronautics, Purdue University, West Lafayette, IN
NEW CONTROL LAWS FOR THE ATTITUDE STABILIZATION OF RIGID BODIES PANAGIOTIS TSIOTRAS *School of Aeronautics and Astronautics, Purdue University, West Lafayette, IN 7907. Abstract. This paper introduces
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 informationROTATIONAL STABILITY OF A CHARGED DIELEC- TRIC RIGID BODY IN A UNIFORM MAGNETIC FIELD
Progress In Electromagnetics Research Letters, Vol. 11, 103 11, 009 ROTATIONAL STABILITY OF A CHARGED DIELEC- TRIC RIGID BODY IN A UNIFORM MAGNETIC FIELD G.-Q. Zhou Department of Physics Wuhan University
More informationChapter 5. Sound Waves and Vortices. 5.1 Sound waves
Chapter 5 Sound Waves and Vortices In this chapter we explore a set of characteristic solutions to the uid equations with the goal of familiarizing the reader with typical behaviors in uid dynamics. Sound
More information1/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 informationCircular Membranes. Farlow, Lesson 30. November 21, Department of Mathematical Sciences Florida Atlantic University Boca Raton, FL 33431
The Problem Polar coordinates Solving the Problem by Separation of Variables Circular Membranes Farlow, Lesson 30 Department of Mathematical Sciences Florida Atlantic University Boca Raton, FL 33431 November
More information9 Kinetics of 3D rigid bodies - rotating frames
9 Kinetics of 3D rigid bodies - rotating frames 9. Consider the two gears depicted in the figure. The gear B of radius R B is fixed to the ground, while the gear A of mass m A and radius R A turns freely
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 information8.012 Physics I: Classical Mechanics Fall 2008
IT OpenCourseWare http://ocw.mit.edu 8.012 Physics I: Classical echanics Fall 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. ASSACHUSETTS INSTITUTE
More informationMODELLING OF FLEXIBLE MECHANICAL SYSTEMS THROUGH APPROXIMATED EIGENFUNCTIONS L. Menini A. Tornambe L. Zaccarian Dip. Informatica, Sistemi e Produzione
MODELLING OF FLEXIBLE MECHANICAL SYSTEMS THROUGH APPROXIMATED EIGENFUNCTIONS L. Menini A. Tornambe L. Zaccarian Dip. Informatica, Sistemi e Produzione, Univ. di Roma Tor Vergata, via di Tor Vergata 11,
More informationWork - kinetic energy theorem for rotational motion *
OpenStax-CNX module: m14307 1 Work - kinetic energy theorem for rotational motion * Sunil Kumar Singh This work is produced by OpenStax-CNX and licensed under the Creative Commons Attribution License 2.0
More informationNewton's second law of motion
OpenStax-CNX module: m14042 1 Newton's second law of motion Sunil Kumar Singh This work is produced by OpenStax-CNX and licensed under the Creative Commons Attribution License 2.0 Abstract Second law of
More informationANNEX 1. DEFINITION OF ORBITAL PARAMETERS AND IMPORTANT CONCEPTS OF CELESTIAL MECHANICS
ANNEX 1. DEFINITION OF ORBITAL PARAMETERS AND IMPORTANT CONCEPTS OF CELESTIAL MECHANICS A1.1. Kepler s laws Johannes Kepler (1571-1630) discovered the laws of orbital motion, now called Kepler's laws.
More informationAy 1 Lecture 2. Starting the Exploration
Ay 1 Lecture 2 Starting the Exploration 2.1 Distances and Scales Some Commonly Used Units Distance: Astronomical unit: the distance from the Earth to the Sun, 1 au = 1.496 10 13 cm ~ 1.5 10 13 cm Light
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 informationUnderstanding Precession
University of Rochester PHY35 Term Paper Understanding Precession Author: Peter Heuer Professor: Dr. Douglas Cline December 1th 01 1 Introduction Figure 1: Bicycle wheel gyroscope demonstration used at
More informationPhysics 70007, Fall 2009 Answers to Final Exam
Physics 70007, Fall 009 Answers to Final Exam December 17, 009 1. Quantum mechanical pictures a Demonstrate that if the commutation relation [A, B] ic is valid in any of the three Schrodinger, Heisenberg,
More informationaction? Pawe l Wegrzyn z Jagellonian University, Institute of Physics Reymonta 4, 30{059 Krakow, Poland Abstract
Boundary terms in Nambu{Goto string action? Leszek Hadasz y Pawe l Wegrzyn z Jagellonian University, Institute of Physics Reymonta 4, 30{059 Krakow, Poland Abstract We investigate classical strings dened
More informationORBITS WRITTEN Q.E. (June 2012) Each of the five problems is valued at 20 points. (Total for exam: 100 points)
ORBITS WRITTEN Q.E. (June 2012) Each of the five problems is valued at 20 points. (Total for exam: 100 points) PROBLEM 1 A) Summarize the content of the three Kepler s Laws. B) Derive any two of the Kepler
More informationCONTINUED FRACTIONS AND THE SECOND KEPLER LAW
CONTINUED FRACTIONS AND THE SECOND KEPLER LAW OLEG KARPENKOV Abstract. In this paper we introduce a link between geometry of ordinary continued fractions and trajectories of points that moves according
More informationMECH 576 Geometry in Mechanics September 16, 2009 Using Line Geometry
MECH 576 Geometry in Mechanics September 16, 2009 Using Line Geometry 1 Deriving Equations in Line Coordinates Four exercises in deriving fundamental geometric equations with line coordinates will be conducted.
More informationGeneral Theoretical Concepts Related to Multibody Dynamics
General Theoretical Concepts Related to Multibody Dynamics Before Getting Started Material draws on two main sources Ed Haug s book, available online: http://sbel.wisc.edu/courses/me751/2010/bookhaugpointers.htm
More informationPhysical Dynamics (PHY-304)
Physical Dynamics (PHY-304) Gabriele Travaglini March 31, 2012 1 Review of Newtonian Mechanics 1.1 One particle Lectures 1-2. Frame, velocity, acceleration, number of degrees of freedom, generalised coordinates.
More information8.012 Physics I: Classical Mechanics Fall 2008
MIT OpenCourseWare http://ocw.mit.edu 8.012 Physics I: Classical Mechanics Fall 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. MASSACHUSETTS INSTITUTE
More information8 Rotational motion of solid objects
8 Rotational motion of solid objects Kinematics of rotations PHY166 Fall 005 In this Lecture we call solid objects such extended objects that are rigid (nondeformable) and thus retain their shape. In contrast
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 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 informationPhysics 106b/196b Problem Set 9 Due Jan 19, 2007
Physics 06b/96b Problem Set 9 Due Jan 9, 2007 Version 3: January 8, 2007 This problem set focuses on dynamics in rotating coordinate systems (Section 5.2), with some additional early material on dynamics
More informationLab #4 - Gyroscopic Motion of a Rigid Body
Lab #4 - Gyroscopic Motion of a Rigid Body Last Updated: April 6, 2007 INTRODUCTION Gyroscope is a word used to describe a rigid body, usually with symmetry about an axis, that has a comparatively large
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 informationthe precession-nutation parameters and Greenwich sidereal time, which Ephemeris Pole (CEP) has to be extended in order to be consistent with
Denition of the Celestial Ephemeris Pole and the Celestial Ephemeris Origin Nicole Capitaine DANOF/UMR8630-CNRS, Observatoire deparis 61, avenue de l'observatoire, 75014 - Paris, France, capitain@danof.obspm.fr
More informationRotational motion of a rigid body spinning around a rotational axis ˆn;
Physics 106a, Caltech 15 November, 2018 Lecture 14: Rotations The motion of solid bodies So far, we have been studying the motion of point particles, which are essentially just translational. Bodies with
More informationDerivation of dual forces in robot manipulators
PERGAMON Mechanism and Machine Theory (1998) 141±148 Derivation of dual forces in robot manipulators V. Brodsky, M. Shoham Department of Mechanical Engineering, Technion Israel Institute of Technology,
More informationOn a WKB-theoretic approach to. Ablowitz-Segur's connection problem for. the second Painleve equation. Yoshitsugu TAKEI
On a WKB-theoretic approach to Ablowitz-Segur's connection problem for the second Painleve equation Yoshitsugu TAKEI Research Institute for Mathematical Sciences Kyoto University Kyoto, 606-8502, JAPAN
More informationCALCULATION OF POSITION AND VELOCITY OF GLONASS SATELLITE BASED ON ANALYTICAL THEORY OF MOTION
ARTIFICIAL SATELLITES, Vol. 50, No. 3 2015 DOI: 10.1515/arsa-2015-0008 CALCULATION OF POSITION AND VELOCITY OF GLONASS SATELLITE BASED ON ANALYTICAL THEORY OF MOTION W. Góral, B. Skorupa AGH University
More informationPhysics 351, Spring 2015, Final Exam.
Physics 351, Spring 2015, Final Exam. This closed-book exam has (only) 25% weight in your course grade. You can use one sheet of your own hand-written notes. Please show your work on these pages. The back
More informationThe Rocket Equation. Lukas Lundin. 26th January 2016
The Rocket Equation Lukas Lundin 26th January 2016 Abstract In this project we study the basics of rocket propulsion and rocket motion in the vicinity of the Earth. Furthermore we will compare dierent
More informationANALYTISK MEKANIK I HT 2016
Karlstads Universitet Fysik ANALYTISK MEKANIK I HT 2016 Kursens kod: FYGB08 Kursansvarig lärare: Jürgen Fuchs rum 21F 316 tel. 054-700 1817 el.mail: juerfuch@kau.se FYGB08 HT 2016 Exercises 1 2016-12-14
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 informationUNIVERSITY OF MARYLAND Department of Physics College Park, Maryland. PHYSICS Ph.D. QUALIFYING EXAMINATION PART II
UNIVERSITY OF MARYLAND Department of Physics College Park, Maryland PHYSICS Ph.D. QUALIFYING EXAMINATION PART II January 22, 2016 9:00 a.m. 1:00 p.m. Do any four problems. Each problem is worth 25 points.
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 informationIntroduction Fundamental definitions Motivation
1 Introduction Fundamental definitions Motivation 1.1 Rotation and global shape of the Earth At a very elementary level, the Earth is considered to be an axially symmetric ellipsoid, rotating with uniform
More information1. Introduction As is well known, the bosonic string can be described by the two-dimensional quantum gravity coupled with D scalar elds, where D denot
RIMS-1161 Proof of the Gauge Independence of the Conformal Anomaly of Bosonic String in the Sense of Kraemmer and Rebhan Mitsuo Abe a; 1 and Noboru Nakanishi b; 2 a Research Institute for Mathematical
More informationGravitation fundamentals. By: Sunil Kumar Singh
Gravitation fundamentals By: Sunil Kumar Singh Gravitation fundamentals By: Sunil Kumar Singh Online: < http://cnx.org/content/col10459/1.2/ > C O N N E X I O N S Rice University, Houston, Texas This
More information10-6 Angular Momentum and Its Conservation [with Concept Coach]
OpenStax-CNX module: m50810 1 10-6 Angular Momentum and Its Conservation [with Concept Coach] OpenStax Tutor Based on Angular Momentum and Its Conservation by OpenStax College This work is produced by
More informationHomework 1. Due whatever day you decide to have the homework session.
Homework 1. Due whatever day you decide to have the homework session. Problem 1. Rising Snake A snake of length L and linear mass density ρ rises from the table. It s head is moving straight up with the
More information1 Matrices and Systems of Linear Equations
Linear Algebra (part ) : Matrices and Systems of Linear Equations (by Evan Dummit, 207, v 260) Contents Matrices and Systems of Linear Equations Systems of Linear Equations Elimination, Matrix Formulation
More informationRotational Motion About a Fixed Axis
Rotational Motion About a Fixed Axis Vocabulary rigid body axis of rotation radian average angular velocity instantaneous angular average angular Instantaneous angular frequency velocity acceleration acceleration
More informationUniform circular motion *
OpenStax-CNX module: m13871 1 * Sunil Kumar Singh This work is produced by OpenStax-CNX and licensed under the Creative Commons Attribution License 2.0 Abstract (UCM) is the basic unit of rotational kinematics
More information! 4 4! o! +! h 4 o=0! ±= ± p i And back-substituting into the linear equations gave us the ratios of the amplitudes of oscillation:.»» = A p e i! +t»»
Topic 6: Coupled Oscillators and Normal Modes Reading assignment: Hand and Finch Chapter 9 We are going to be considering the general case of a system with N degrees of freedome close to one of its stable
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 informationAPPENDIX B SUMMARY OF ORBITAL MECHANICS RELEVANT TO REMOTE SENSING
APPENDIX B SUMMARY OF ORBITAL MECHANICS RELEVANT TO REMOTE SENSING Orbit selection and sensor characteristics are closely related to the strategy required to achieve the desired results. Different types
More informationContents. Objectives Torque on an Object Rotational Kinetic Energy Yo yo Rolling on an Incline Physical Pendulum Angular Momentum and Torque Recap
Physics 121 for Majors Class 21 Rotating Objects Last Class We learned to find angular momentum and torques of point masses and objects. We learned how to use torques and forces to solve problems with
More informationGeneralization of the Euler Angles
The Journal of the Astronautical Sciences, Vol. 51, No. 2, AprilJune 2003, pp. 123123 Generalization of the Euler Angles Malcolm D. Shuster 1 and F. Landis Markley 2 Abstract It is shown that the Euler
More informationResearch Article Variation of the Equator due to a Highly Inclined and Eccentric Disturber
Hindawi Publishing Corporation Mathematical Problems in Engineering Volume 009, Article ID 467865, 10 pages doi:10.1155/009/467865 Research Article Variation of the Equator due to a Highly Inclined and
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 informationRELATIVE MOTION ANALYSIS: VELOCITY (Section 16.5)
RELATIVE MOTION ANALYSIS: VELOCITY (Section 16.5) Today s Objectives: Students will be able to: a) Describe the velocity of a rigid body in terms of translation and rotation components. b) Perform a relative-motion
More informationAppendix. Vectors, Systems of Equations
ppendix Vectors, Systems of Equations Vectors, Systems of Equations.1.1 Vectors Scalar physical quantities (e.g., time, mass, density) possess only magnitude. Vectors are physical quantities (e.g., force,
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 informationHM NAUTICAL ALMANAC OFFICE NAO TECHNICAL N OTE
Council for the Central Laboratory of the Research Councils HM NAUTICAL ALMANAC OFFICE NAO TECHNICAL N OTE No. 73 2003 June Approximate ephemerides and phenomena of natural satellites by D.B. Taylor and
More informationAstronomy 6570 Physics of the Planets. Precession: Free and Forced
Astronomy 6570 Physics of the Planets Precession: Free and Forced Planetary Precession We have seen above how information concerning the distribution of density within a planet (in particular, the polar
More informationFineman Honors Physics Final Study Guide
All Science Tests are on Wednesday, June 17 th. Students who take more than one Science class will take their second science final on Thursday, June 18 from 8:00-10:00 AM in the Library. The Honors Physics
More informationDynamics 12e. Copyright 2010 Pearson Education South Asia Pte Ltd. Chapter 20 3D Kinematics of a Rigid Body
Engineering Mechanics: Dynamics 12e Chapter 20 3D Kinematics of a Rigid Body Chapter Objectives Kinematics of a body subjected to rotation about a fixed axis and general plane motion. Relative-motion analysis
More informationProblem 1 : Solution/marking scheme Two Problems in Mechanics (10 points)
Problem 1 : Solution/marking scheme Two Problems in Mechanics (10 points) Part A. The Hidden Disk (3.5 points) A1 (0.8 pt) Find an expression for b as a function of the quantities (1), the angle φ and
More informationLie Groups for 2D and 3D Transformations
Lie Groups for 2D and 3D Transformations Ethan Eade Updated May 20, 2017 * 1 Introduction This document derives useful formulae for working with the Lie groups that represent transformations in 2D and
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 informationinertia of a body, principal axes of inertia, invariants of an inertia tensor, and inertia triangle inequalities are illustrated and discussed.
Preface This book belongs to a series of three books written simultaneously (the remaining two are titled Classical Mechanics: Dynamics and Classical Mechanics: Applied MechanicsandMechatronics). This
More informationAnswers to Problem Set Number MIT (Fall 2005).
Answers to Problem Set Number 6. 18.305 MIT (Fall 2005). D. Margetis and R. Rosales (MIT, Math. Dept., Cambridge, MA 02139). December 12, 2005. Course TA: Nikos Savva, MIT, Dept. of Mathematics, Cambridge,
More informationMechanics Lecture Notes
Mechanics Lecture Notes Lectures 0 and : Motion in a circle. Introduction The important result in this lecture concerns the force required to keep a particle moving on a circular path: if the radius of
More informationAstronomy 6570 Physics of the Planets
Astronomy 6570 Physics of the Planets Planetary Rotation, Figures, and Gravity Fields Topics to be covered: 1. Rotational distortion & oblateness 2. Gravity field of an oblate planet 3. Free & forced planetary
More informationChapter 4. Motion in Two Dimensions
Chapter 4 Motion in Two Dimensions Kinematics in Two Dimensions Will study the vector nature of position, velocity and acceleration in greater detail Will treat projectile motion and uniform circular motion
More informationProperties of surfaces II: Second moment of area
Properties of surfaces II: Second moment of area Just as we have discussing first moment of an area and its relation with problems in mechanics, we will now describe second moment and product of area of
More informationSolutions to Math 53 Math 53 Practice Final
Solutions to Math 5 Math 5 Practice Final 20 points Consider the initial value problem y t 4yt = te t with y 0 = and y0 = 0 a 8 points Find the Laplace transform of the solution of this IVP b 8 points
More information