GENERALIZED WHITTAKER S EQUATIONS FOR HOLONOMIC MECHANICAL SYSTEMS

Similar documents
Tensors, and differential forms - Lecture 2

Physics 235 Chapter 7. Chapter 7 Hamilton's Principle - Lagrangian and Hamiltonian Dynamics

Problem 01 C. 50. What is the length of a line segment with a slope of 4/3, measured from the y-axis to a point (6,4)? B. 25 D. 75 A.

PHYS 705: Classical Mechanics. Euler s Equations

Artificial Intelligence & Neuro Cognitive Systems Fakultät für Informatik. Robot Dynamics. Dr.-Ing. John Nassour J.

1/30. Rigid Body Rotations. Dave Frank

Physics 5153 Classical Mechanics. Canonical Transformations-1

Physical Dynamics (SPA5304) Lecture Plan 2018

Phys 7221 Homework # 8

1. To analyze the deformation of a conical membrane, it is proposed to use a two-dimensional conical-polar coordinate system ( s,

PLANAR KINETICS OF A RIGID BODY FORCE AND ACCELERATION

Review: control, feedback, etc. Today s topic: state-space models of systems; linearization

Differential Kinematics

PHYS 705: Classical Mechanics. Introduction and Derivative of Moment of Inertia Tensor

Part 8: Rigid Body Dynamics

On local normal forms of completely integrable systems

Engineering Physics 1 Dr. B. K. Patra Department of Physics Indian Institute of Technology-Roorkee

3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1

SCHEME OF BE 100 ENGINEERING MECHANICS DEC 2015

Rotational Kinematics and Dynamics. UCVTS AIT Physics

Transcendental cases in stability problem. Hamiltonian systems

CP1 REVISION LECTURE 3 INTRODUCTION TO CLASSICAL MECHANICS. Prof. N. Harnew University of Oxford TT 2017

N mg N Mg N Figure : Forces acting on particle m and inclined plane M. (b) The equations of motion are obtained by applying the momentum principles to

PLANAR KINETIC EQUATIONS OF MOTION (Section 17.2)

Two small balls, each of mass m, with perpendicular bisector of the line joining the two balls as the axis of rotation:

Rigid body dynamics. Basilio Bona. DAUIN - Politecnico di Torino. October 2013

PHY 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

Part II. Classical Dynamics. Year

Physical Dynamics (PHY-304)

Tensors - Lecture 4. cos(β) sin(β) sin(β) cos(β) 0

ON DIVISION ALGEBRAS*

and u and v are orthogonal if and only if u v = 0. u v = x1x2 + y1y2 + z1z2. 1. In R 3 the dot product is defined by

2. A die is rolled 3 times, the probability of getting a number larger than the previous number each time is

Computational Spectroscopy III. Spectroscopic Hamiltonians

A B Ax Bx Ay By Az Bz

Ch 4 Differentiation

Biderivations of the Algebra of Strictly Upper Triangular Matrices over a Commutative Ring

Video 3.1 Vijay Kumar and Ani Hsieh

112 Dynamics. Example 5-3

Input: A set (x i -yy i ) data. Output: Function value at arbitrary point x. What for x = 1.2?

1 Hamiltonian formalism

Math512 PDE Homework 2

Rotational Motion. Chapter 4. P. J. Grandinetti. Sep. 1, Chem P. J. Grandinetti (Chem. 4300) Rotational Motion Sep.

16.2. Line Integrals

Physics 741 Graduate Quantum Mechanics 1 Solutions to Midterm Exam, Fall x i x dx i x i x x i x dx

8 Velocity Kinematics

cc) c, >o,, c, =c7 $"(x\ <o - S ~ Z

16. Rotational Dynamics

or - CHAPTER 7 Applications of Integration Section 7.1 Area of a Region Between Two Curves 1. A= ~2[0- (x :2-6x)] dr=-~2(x 2-6x) dr

Chapter 11. Angular Momentum

Appendix A Review of Vector Algebra and Derivatives of Vectors

Name Class. (a) (b) (c) 4 t4 3 C

Contents. Physical Properties. Scalar, Vector. Second Rank Tensor. Transformation. 5 Representation Quadric. 6 Neumann s Principle

MATH 391 Test 1 Fall, (1) (12 points each)compute the general solution of each of the following differential equations: = 4x 2y.

RAJASTHAN P.E.T. MATHS-1995

Problem 1. Mathematics of rotations

9 Kinetics of 3D rigid bodies - rotating frames

Vector Calculus. A primer

Properties of surfaces II: Second moment of area

SOME CHARACTERIZATION

Analytical Dynamics: Lagrange s Equation and its Application A Brief Introduction

I Results in Mathematics

LAGRANGIAN AND HAMILTONIAN

THE UNIVERSITY OF BRITISH COLUMBIA DEPARTMENT OF STATISTICS TECHNICAL REPORT #253 RIZVI-SOBEL SUBSET SELECTION WITH UNEQUAL SAMPLE SIZES

Perturbative QCD. Chul Kim. Seoultech. Part I : Introduction to QCD Structure functions for DIS

Dynamics 12e. Copyright 2010 Pearson Education South Asia Pte Ltd. Chapter 20 3D Kinematics of a Rigid Body

Classical Mechanics. Luis Anchordoqui

THE STRUCTURE OF A RING OF FORMAL SERIES AMS Subject Classification : 13J05, 13J10.

3-D Kinetics of Rigid Bodies

Strain analysis.

Chapter 3 Convolution Representation

Math Review Night: Work and the Dot Product

WBJEE Answer Keys by Aakash Institute, Kolkata Centre

Mechanics. Moment of inertia / Steiner s theorem Dynamics. What you need:

(Received March 30, 1953)

y mx 25m 25 4 circle. Then the perpendicular distance of tangent from the centre (0, 0) is the radius. Since tangent

MATH 2413 TEST ON CHAPTER 4 ANSWER ALL QUESTIONS. TIME 1.5 HRS.

r(j) -::::.- --X U.;,..;...-h_D_Vl_5_ :;;2.. Name: ~s'~o--=-i Class; Date: ID: A

AN INTRODUCTION TO CURVILINEAR ORTHOGONAL COORDINATES

COMMONWEALTH OF AUSTRALIA Copyright Regulations 1969

2.20 Fall 2018 Math Review

Rotational & Rigid-Body Mechanics. Lectures 3+4

Review of Vector Analysis in Cartesian Coordinates

Translational and Rotational Dynamics!

Laws of gyroscopes / cardanic gyroscope

PHYS 705: Classical Mechanics. Rigid Body Motion Introduction + Math Review

Classical Mechanics III (8.09) Fall 2014 Assignment 3

Analytical Mechanics for Relativity and Quantum Mechanics

FIELDS IN SQUARE HELMHOLTZ COILS

This content downloaded from on Mon, 04 Feb :52:49 UTC All use subject to

xi is asymptotically equivalent to multiplication by Xbxi/bs, where S SOME REMARKS CONCERNING SCHRODINGER'S WA VE EQ UA TION

Section 27. The Central Limit Theorem. Po-Ning Chen, Professor. Institute of Communications Engineering. National Chiao Tung University

Power system modelling under the phasor approximation

Supplementary information I Hilbert Space, Dirac Notation, and Matrix Mechanics. EE270 Fall 2017

MATH1013 Calculus I. Introduction to Functions 1

MATH 228: Calculus III (FALL 2016) Sample Problems for FINAL EXAM SOLUTIONS

Extended Application Of Time-Varying Observer For Rigid Rotors Unbalance Estimation During Acceleration

Second quantization: where quantization and particles come from?

Analysis-3 lecture schemes

Lecture Notes #10. The "paradox" of finite strain and preferred orientation pure vs. simple shear

Transcription:

Intnat. J. Mh. Math. S ci. Vol. (1978) 245253 245 GENERALIZED WITTAKER S EQUATIONS FOR OLONOMIC MECANICAL SYSTEMS MUNAWAR USSAIN Department of Mathematics Government College Lahore, Pakistan (Received June 23, 1977) ABSTRACT. In this paper the classical theorem "a conservative holonomic dynamic system is invariantly connected with a certain differential form" is generalized to group variables. This general theorem is then used to reduce the order of a amiltonian system by the use of the integral of energy. Equations of motion of the reduced system so obtained are derived which are the socalled generalized Whittaker s equations. Finally an example is given as an application of the theory. KEY WORDS AND PRASES. Generalized Whlttaker s equations, olonomic Systems. A:M.S(OS) SUBJECT CLASSIFICATION (19,70). 70F20 i. INTRODUCTION It is known [4] that canonical equations for a conservative holonomic system whose amiltonian is are obtained by forming the first Pfaff s system of differential equations of the differential form, Pidqi dr, (i 1,2,,n),

246 MLrNAWAR USSAIN where PI P2 Pn coordinates ql q2 qn are the generalized momenta corresponding to the generalized of the system and summation over a repeated suffix is implied. This result leads to the theorem "A dynamical system is invarlantly connected with the differential form Pidqi dt". Further this theorem has been used to reduce the order of the system by means of the integral of energy. Canonical equations of the reduced system so obtained are known as Whittaker s equations. In what follows in this section we state a few basic results from the theory of group variables in order to generalize the above mentioned results. Consider a conservative holonomlc system having n degrees of freedom and whose position is specified by group variables x I, x2,, Let Xn. n I, 2,,nn be the parameters of real displacement and X I, X 2,, be the corresponding Xn displacement operators expressed by the relations J where i are functions of x I, x 2,, Xn, J Xl" i 8xj (i,j 1,2,,n) then for an arbitrary function () f(x I, x 2,, x n, t) the infinitesimal change df is expressed by df [Xo( f) + nixi(f)]dt (2) where The X s satisfy the relations (Xo, X i) O, (Xi, Xj) CiJk (k,2,,n). (3) Putting f xj in (2), we get J dx.,j.n {i (4) dt Since the operators X i are independent therefore the matrix =o=singular and consequently (4) yields n i Aij xj (5)

GENERALIZED WITTAKER S EQUATIONS 24 7 Let L be the Lagranglan of the system then the canonical equations of the system as obtained in i, 2 are ri i where (i,j,k 1,2,,n), and Yl d Cjlk rj Yk Xi (), (6) L n (Xl, x 2,,Xn; yl,,yn ) nlyi is the amiltonlan of the system and is equal to L (7) the total energy of the system. 2. DERIVATION OF CANONCAL EQUATIONS FROM A CERTAIN DIFERE..LkL, FORM. In order to establish the invarlant relation between the system (6) and a certain differential form we prove the following theorem: TEOREM. The system of equations (6) is equivalent to the first Pfaff s system of differential equations of the dlfferentlal form (nly 2 )dt. PROOF. We put or using (5), we obtain 8d (niyi )dt therefore 6 d yiaij d xj dt (8) O Yl AiJ x 6t (9) where d and denote two independent variations in each of the variables x I, x2,, Xn, t. The bilinear corearlant of YI Yn (8) is given by Od_dO= Yi[Aijdxj AiJ yldt]+xk[yl X_k Aik. t_dy i dxj] Xk dxj AikYi xj + t[d dt] (10)

248 MUNAWAR USSAIN where we have used the relations dx d x (i i i 1,2,,n) dt dt. Equating to zero the coefficients of I n we get the first Pfaff s system of equation in the form dt 0, (i 1,2,,n) Aijdxj Yl A Yl dxj dt dyl Aik Yi BAlk dxj 0 (Ii) (12) d dt, (i i,2,,n) (13) By vrtue of (5), the equations (11) asse the fo (i 1,2,,n). With the help of (4), the equations (12) become k. dy i Ynm Aj Ak J k k mi Ym dq x mi k 8xj i x k (14) which, by means of the relations (I) and (3), finally takes the form dy i rj Yk Cik XI()" (15) dt The relation (13) is a consequence of (14) and (15) and skew symmetric property of Cjik with respect to the first two indices. Since the equations (14) and (15) are identical with (6) the theorem is thus proved. 3. GENERALIZED WITTAKER S EQUATIONS Assume that does not involve the time explicitly and / h O, (16) is the integral of energy of the system. Let the equation (16) be solved for the variable Yl so that it is algebraically equivalent to K(Xl Xn Y2 Yn t, h) + Yl 0. (17)

GENERALIZED WITTAKER S EQUATIONS 249 The differential form associated with the system is (nl Yl + n2 Y2 + + nn Yn + h)dt, where the variables x 1,, x n, Yl Yn h are connected by (17); the differential form can therefore be written as (2 Y2 + n3 Y3 + + rn Yn + h)dt n I where we can regard (Xl,, Xn, Y2 Yn in the phase space. If we express (18) in the form lh_k] I dt [ Y2 + + n Yn + i and put nldt dz then we take T as the new time variable and I 1 as the parameters of real displacement, K dt (18) h, t) as the 2n+l variables n 2 (19) n n i i 2 i n i the corresponding displacement operators and new momenta are respectively X o, X I,, Xn and h, YI Yn" Using the result of section (2), the differential equation corresponding to the form (19) are dyp 8K p %yp dt rl] Yk pk Cj Xp(K) (p 2 3,n) (20) dt 8K dh d r )h Oo The last pair of equations can be separated from the rest of the system since the first (2n2) equations do not involve t and h is a constant. The equations (20) can be further simplified to take the form K[Clpl + nr Crpl] + Yr Clpr + nryq Crpq Xp(K) (21) (p,q,r 2,3,,n).

250 MUNAWAR USSAIN The original differential equations can therefore be replaced by the reduced system (21) which has only ni degrees of freedom. The equations (21) are the desired Whittaker s equations. 4. AN EXAMPLE Consider a rigid body which is moving about one of its fixed points 0 under the action of gravity. We introduce a fixed frame of reference Oxyz such that Oz is vertically upwards and a moving frame Ox y z which coincides with the principal axes of inertia of the body at O. Let us choose the Eulerian angles 8,, (8 is the angle of nutation, the angle of precession and the angle of proper rotation) as the group variables which specify the position of the body at time t. Obviously the dynamical system under conslderatlon is a conservative one and it has three degrees of freedom. Choosing the parameters of real displacement as the components of angular velocity. os + along the moving axes, we have the relations sin sin n 2 Sln + Sin 8 Cos + $ co,, Consequently the displacement operators X X are given by 2 3 X Cos Cosec 1 + e sin Cot e sin X Sin Cosec 2 e Cos Cot e Cos + x 3 which satisfy the commutation relations (X 1, X2) X XlX2X2X1 3 X 1 (X 2, X3) X2X3 X3X2 (22) (23) (X 3, X1) X3X1 XlX3 X 2

GENERALIZED WITTAKER S EQUATIONS 251 The nonvanlshlng C s are therefore expressed by the relations C123 C213 i, C231 C321 1, C312 C132 i. Let T and U denote the kinetic and potential energies of (24) the system respectively, then T (A n + Bn 2 + Cn3), (25) U Mg(x Sln B Sin + y Sin B Cos + z Cos B where A, B, C are te principal moments of inertia at O; x, y, z are the coordinates of the centre of gravity of the body wth respect to the moving axls and M Is the mass of the body. Using (25), we have the Lagranglan L and momenta expressed by the relations: Yl Y2 Y3 L T U (An + Bn + Cn ) + Mg ( Sin 0 Sin + y Sin O Cos +z Cos O), Yl A 1 Y2 B C 2 Y3 n (26) 3. In view of (26) the amiltonlan is given by = + + Mg(x Sin e sin + y Sin e Cos + z Cos e) (27) "A B Using (6), (22), (24), (26), and (27), canonlcal equations of the system are Y Y 73 n 1 1, n 2 2, n 3 A B C dyl BC y2y3 dt BC + g (" Cos e " Sin y e Cos ), dy2 CA y3y I + Mg ( Cos O + Sin e Sin ) dt CA dy 3 A B yly 2 + Mg Sin e(x Cos Sin ). (28)

252 MUNAWAR USSAIN Now the relation (16) gives + and consequently + 2Ms( Sin 6 Sin + y Sin 6 Cos + z Cos 6) + 2h 0, Yl A[2S( Sin 0 Sin, + Sin 0 Cos + Cos 0) 2 Y 2hi, Comparing this relation wth (17), we get K A[2ME(" Sin e Sin + y Sin e Cos $ + z Cos e).. 2h] ;. B C Therefore by the application of (21) the canonical equations of the system reduce to 3K 2 Y2 3K r3 3 dy2 Y3 n3 Yl X2(K), d r dy3 Y2 + n2 Yl X3(K). (29) Now dy2 dy2 K A dy2 d dt h K dt dy3 dy3 K A dy3 de dt 9h K dt X2(K) =. Mg(x Cos 0 + z Sin e sla ), x3 ) " t s cos s ) s e.

Therefore equation (29) GENERALIZED WITTAKER S EUATIONS 25 3 assume the form A A r BY2 rl y3 _._.dy2 KA_._y 3 + Ms( Cos e dt CA + Sin e Sin,), dy3 dt K(BA) Y2 + S; Sin e(" Cos Sin ). " These are the Whittaker s equations for the system under consideration. REFERENCES 1. Cetaev, N. G. On the Equations of Poincare, Prlkl. Mg.t...Me.h.(1941) 253262. 2. ussain, M. amlltonjacobi Theorem in Group Variables, Journal of Applied Mathematics and Physics (ZAMP), Vol 27, (1976) 285287. 3. Poincare,. On a New Form of the Equations of Mechanics, C. R. Acad. Sci. 132 (1901) 369371. 4. Whittaker, E. T. Analytic.al Dynamics of Particle.s and Rigid Bodies, Cambridge University Press, 1961.