Page:EB1911 - Volume 17.djvu/1009

From Wikisource
Jump to navigation Jump to search
This page has been proofread, but needs to be validated.
990 
MECHANICS
[KINETICS


For a real solution we must have k2 < 1, which is equivalent to 2BT > Γ2. If the initial conditions are such as to make 2BT < Γ2, we must interchange the forms of p and r in (7). In the present case the instantaneous axis returns to its initial position in the body whenever φ increases by 2π, i.e. whenever t increases by 4K/σ, when K is the “complete” elliptic integral of the first kind with respect to the modulus k.

The elliptic functions degenerate into simpler forms when k2 = 0 or k2 = 1. The former case arises when two of the principal moments are equal; this has been sufficiently dealt with in § 19. If k2 = 1, we must have 2BT = Γ2. We have seen that the alternative 2BT ≷ Γ2 determines whether the polhode cone surrounds the principal axis of least or greatest moment. The case of 2BT = Γ2, exactly, is therefore a critical case; it may be shown that the instantaneous axis either coincides permanently with the axis of mean moment or approaches it asymptotically.

When the origin of the moving axes is also in motion with a velocity whose components are u, v, w, the dynamical equations are

dξ rη + qζ = X,   dη pζ + rχ = Y,   dζ qχ + pη = Z,
dt dt dt
(11)


dλ rμ + qνwη + vζ = L,   dμ pν + rλ- uζ + wξ = M,   dν qλ + pμvξ + uη = N.
dt dtdt
(12)


To prove these, we may take fixed axes O′x′, O′y′, O′z′ coincident with the moving axes at time t, and compare the linear and angular momenta ξ + δξ, η + δη, ζ + δζ, λ + δλ, μ + δμ, ν + δν relative to the new position of the axes, Ox, Oy, Oz at time t + δt with the original momenta ξ, η, ζ, λ, μ, ν relative to O′x′, O′y′, O′z′ at time t. As in the case of (2), the equations are applicable to any dynamical system whatever. If the moving origin coincide always with the mass-centre, we have ξ, η, ζ = M0u, M0v, M0w, where M0 is the total mass, and the equations simplify.

When, in any problem, the values of u, v, w, p, q, r have been determined as functions of t, it still remains to connect the moving axes with some fixed frame of reference. It will be sufficient to take the case of motion about a fixed point O; the angular co-ordinates θ, φ, ψ of Euler may then be used for the purpose. Referring to fig. 36 we see that the angular velocities p, q, r of the moving lines, OA, OB, OC about their instantaneous positions are

p = θ.  sin φ − sin θ cos φψ. ,   q = θ.  cos φ + sin θ sin φψ. , r = φ.  + cos θψ. ,
(13)

by § 7 (3), (4). If OA, OB, OC be principal axes of inertia of a solid, and if A, B, C denote the corresponding moments of inertia, the kinetic energy is given by

2T = A (θ.  sin φ − sin θ cos φψ. )2 + B (θ.  cos φ + sin θ sin θψ)2 + C (φ.  + cos θψ. )2.
(14)

If A = B this reduces to

2T = A (θ.2 + sin2 θ ψ.2) + C (φ.  + cos θ ψ. )2;
(15)

cf. § 20 (1).

§ 22. Equations of Motion in Generalized Co-ordinates.—Suppose we have a dynamical system composed of a finite number of material particles or rigid bodies, whether free or constrained in any way, which are subject to mutual forces and also to the action of any given extraneous forces. The configuration of such a system can be completely specified by means of a certain number (n) of independent quantities, called the generalized co-ordinates of the system. These co-ordinates may be chosen in an endless variety of ways, but their number is determinate, and expresses the number of degrees of freedom of the system. We denote these co-ordinates by q1, q2, . . . qn. It is implied in the above description of the system that the Cartesian co-ordinates x, y, z of any particle of the system are known functions of the q’s, varying in form (of course) from particle to particle. Hence the kinetic energy T is given by

2T = Σ {m (2 + 2 + ż2) }
= a1112 + a2222 + . . . + 2a1212 + . . .,
(1)

where

arr = Σ [ m { ( x 2 + y 2 + z 2 } ],
qr qr qr
ars = Σ { m x   x + y   y + z   z ) } = asr.
qr qs qr qs qr qs
(2)


Thus T is expressed as a homogeneous quadratic function of the quantities 1, 2, . . . n, which are called the generalized components of velocity. The coefficients arr, ars are called the coefficients of inertia; they are not in general constants, being functions of the q’s and so variable with the configuration. Again, If (X, Y, Z) be the force on m, the work done in an infinitesimal change of configuration is

Σ (Xδx + Yδy + Zδz) = Q1δq1 + Q2δq2 + . . . + Qnδqn,
(3)

where

Qr = Σ X x + Y y + Z z .
qr qr qr
(4)

The quantities Qr are called the generalized components of force.

The equations of motion of m being

mẍ = X,   mÿ = Y,   mz̈ = Z,
(5)

we have

Σ { m x + ÿ y + z ) } = Qr.
qr qr qr
(6)

Now

= x 1 + x 2 + . . . + x n,
q1 q2 qn
(7)

whence

/r = x/qr   (8)

Also

d x = 2x 1 + 2x 2 + . . . + 2x r = .
dt qr q1qr q2qr qnqr qr
(9)


Hence

x = d x d x = d .
qr dt qr dt qr dt r qr
(10)


By these and the similar transformations relating to y and z the equation (6) takes the form

d ∂T ∂T = Qr.
dt r qr
(11)

If we put r = 1, 2, . . . n in succession, we get the n independent equations of motion of the system. These equations are due to Lagrange, with whom indeed the first conception, as well as the establishment, of a general dynamical method applicable to all systems whatever appears to have originated. The above proof was given by Sir W. R. Hamilton (1835). Lagrange’s own proof will be found under Dynamics, § Analytical. In a conservative system free from extraneous force we have

Σ (X δx + Y δy + Z δz) = −δV,
(12)

where V is the potential energy. Hence

Qr = − ∂V ,
qr
(13)

and

d ∂T ∂T = − ∂V .
dt r qr qr
(14)


If we imagine any given state of motion (1, 2 . . . n) through the configuration (q1, q2, . . . qn) to be generated instantaneously from rest by the action of suitable impulsive forces, we find on integrating (11) with respect to t over the infinitely short duration of the impulse

∂T/r = Qr′,   (15)

where Qr′ is the time integral of Qr and so represents a generalized component of impulse. By an obvious analogy, the expressions ∂T/∂r may be called the generalized components of momentum; they are usually denoted by pr thus

pr = ∂T / ∂r = a1r1 + a2r2 + . . . + anrn.
(16)

Since T is a homogeneous quadratic function of the velocities 1, 2, . . . n, we have

2T = ∂T 1 + ∂T 2 + . . . + ∂T n = p12 + p22 + . . . + pnn.
1 2 n
(17)

Hence

2dT/dt= 11 + 22 + . . . + nn + 11 + 22 + . . . + nn
= ∂T + Q1 1 + ∂T + Q2 2 + . . . + ∂T + Qn n + ∂T 1 + ∂T 2 + . . . + ∂T n
1 2 n 1 2 n
=dT/dt + Q11 + Q22 + . . . + Qnn,
(18)



or

dT/dt = Q11 + Q22 + . . . + Qnn.   (19)