Russian Journal of Nonlinear Dynamics, 2019, vol. 15, no. 3, pp. 285-292. Full-texts are available at http://nd.ics.org.ru DOI: 10.20537/nd190307
MATHEMATICAL PROBLEMS OF NONLINEARITY
MSC 2010: 70E17
Precessional Motion of a Rigid Body Acted upon by Three Irreducible Fields
A. M. Hussein
We consider a quit general problem of motion of an asymmetric rigid body about a fixed point, acted upon by an irreducible skew combination of gravitational, electric and magnetic fields. Two of those three fields are uniform and the third has a more complicated structure. The existence of precessional motions about a nonvertical axis is established. Conditions on the parameters of the system are obtained. An alternative physical interpretation is given in the framework of the problem of motion of a rigid body immersed in an incompressible perfect fluid, acted upon by torques due to two uniform fields.
Keywords: rigid body, precessional motion, three irreducible fields
1. Introduction
The motion of a rigid body, fixed at one of its points is said to be a regular precession if the body rotates uniformly about an axis, called the figure axis, fixed in it, while that axis precesses with the same, or another, uniform angular velocity about an axis fixed in space, the precession axis. Routh [12] showed that regular precessions of an asymmetric rigid body, having for the axis of precession the vertical and for the figure axis one of the three principal axes of inertia, are dynamically impossible.
The possibility of the regular precession of an asymmetric heavy rigid body about an axis inclined to the vertical was pointed out by Grioli [3, 4] on a dynamical basis. Later, Gulyaev [5] obtained Grioli's solution by direct integration of the Euler-Poisson equations of motion. The question of the existence of other Grioli-type solutions of the classical problem of the motion of a rigid body has been investigated by many authors, see, e.g., [7, 8]. Grioli's result was generalized by Kharlamova [9] to the case of a gyrostat acted upon by a single gravity field.
Received July 22, 2019 Accepted August 01, 2019
Ashraf M. Hussein ahussein@sci.kfs.edu.eg
Faculty of Science and Arts, University of Bisha Al-Namas 61977, Saudi Arabia
Department of Mathematics, Faculty of Science, Kafrelsheikh University Kafr El-Sheikh 33516, Egypt
The existence of a Grioli-type precession in the problem of the motion of a rigid body in a liquid is shown in [13, 14] when the figure axis is perpendicular to the precession axis. The case when the figure and precession axes are not perpendicular to each other is considered in [10].
The equations of motion of the problem of motion of a rigid body by inertia in a liquid were given in a new form [16] isomorphic to those of a body moving about a fixed point and subject to the effect of a conservative combination of coaxial potential and gyroscopic forces. Physical interpretation of those forces involves a uniform magnetic field, which induces a scalar magnetic potential of magnetized parts of the body, a vector potential representing Lorentz forces on fixed electric charge distribution on the surface of the body and two other fields, gravitational and electric, whose potentials are symmetric around the direction of the magnetic field and involve only first and second harmonics. The analogy of the two problems has made it possible to construct integrable cases and particular solutions in the second problem from all known ones of the first. The second problem turned out to be extremely flexible and has gone a long way in development and generalization in the same framework of axial combinations of two or three fields. For details the reader may be referred to [17] and [18].
Although enormous effort was invested in the problem of motion of a rigid body about a fixed point in a uniform field and in an axial combination of nonuniform fields, the problem involving skew fields acting on the body has escaped attention most of the time.
Conditions for the existence of a regular precession in the problem of the motion of a rigid body (gyrostat) acted upon by an irreducible combination of three fields are obtained in [20]. This generalizes results obtained in [19] for the case of a single rigid body acted upon by two irreducible fields.
The present article considers a rigid body moving in an ideal incompressible fluid while acted upon by an irreducible combination of two uniform fields. We show that the body can perform a Grioli-type precession and obtain some previous results as a special case.
2. Equations of Motion
To describe the attitude dynamics of a rigid body, two coordinate systems OXYZ and Oxyz are introduced: the system OXYZ is fixed in space with unit vectors a, p and y along its axes, while the system Oxyz is rigidly attached to the body with unit vectors i, j and k which, in general, are not directed along the body principal axes at O. So, the inertia tensor is given by
( A -F -e\
-F B -D . (2.1)
y-E -D C j
We consider a dynamical problem defined by the Euler equation
,r dV a dV dV . .
Ill) + u> x (Iu> + K~f + k) = a. x —--1- /3 x —— + 7 x —, 2.2
da d p dY
where u denotes the angular velocity of the body, K = (Kj) is a 3 x 3 symmetric matrix, k is the gyrostatic moment and the dot denotes differentiation with respect to time. Equations (2.2) can be adapted by a proper choice of the system's parameters and the potential V to describe many problems in rigid body dynamics. Let us consider the case when
V = aa + b- (3+ rij7 + sJ -7, (2.3)
where J = (Jj) is a 3 x 3 symmetric matrix. Then the problem under consideration is defined by the equations
Iu + u x (Iu + K7 + k) = a x a + ¡3 x b + 7 x (J7 + s), (2.4a)
a + u x a = 0, ¡3 + u x (3 = 0, 77 + u x 7 = 0, (2.46)
where the dynamical equation (2.4a) is augmented by the Poisson kinematical equations (2.46). A physical interpretation of the problem (2.4) when a = b = 0 is given in [16]. It is shown that the problem describes the inertial motion of a rigid body in an ideal incompressible fluid. In the generic case, the potential is interpreted as representing the action of three crossed fields on a body about a fixed point; one of them is nonlinear (axisymmetric about 7) and the other two are uniform in the a and 3 directions, with their centers at a and b, respectively.
3. Precessional Motion
Without loss of generality, we let 7 and k be the precession and figure axes, respectively. If we denote the constant angle between the two axes by 90, then the regular precession is defined by
7 ■ k = a0, u = Q (7 + k), (3.1)
where a0 = cos 90. The precession (3.1) can be defined in terms of Euler's angles defining the angular orientation of the system Oxyz relative to OXYZ: p denotes the angle of precession around the Z-axis, 9 is the angle of nutation between the z- and Z-axes, and 0 is the angle of rotation of the body around the z-axis. We first recall that (see [11]):
u = [pp sin 9 sin 0 + d cos ^ i + (ip sin d cos 0 — d sin ^ j + [ip cos d + k,
a = (cos p cos 0 — cos d sin p sin 0) i — (cos p sin 0 + cos d sin p cos 0) j + sin d sin p k, (3 2) 3 = (sinp cos 0 + cos d cos p sin 0) i — (sin p sin 0 — cos d cos p cos 0) j — sin d cos pk, 7 = sin 9 sin 0 i + sin 9 cos 0 j + cos 9 k.
The motion (3.1) can be defined by
9 = 9o, p = 0 = u = Q(t — to). (3.3)
Substituting (3.3) into (3.2), we get
u
= (a0Qsinu, aOQcosu, eiQ) ,
a = (e1 cos2 u — a0, —e1 sin ucos u, a'0 sin u) , 3 = (e1 sin u cos u, —61 sin2 u + a0, —a'0 cos u) , 7 = (a0 sin u, a'0 cos u, a0),
(3.4)
where
a0 = v/1 — a0 = sin 90, e1 = 1 + a0 = 2cos2(90/2). (3.5)
For a0 = ±1, the motion degenerates to a permanent rotation about 7 = ±k. Since we are interested here in precessional motions, we will assume in what follows that a0 = ±1.
4. The solution
One can easily see that (3.4) is a solution of the kinematical Poisson equations (2.46). So, we turn our attention now to the dynamical problem. Substituting (3.4) into the dynamical equation (2.4a), we get the set of three scalar equations
- bo {[A + ei (C - B)] Q2 + K3Q + aoJ - J33) - b2 - S3} cosu
- bo [(ei + 1) FQ2 + ao J12 + a^ sinu - e^DQ2 + e^Q - eib3 - ao(ao J23 - b3 + S2) = 0,
- [b0(EQ2 + J13) - e1a3] cos2 u + [b0 (DQ2 + J23) + e1b3] cos u sinu
+ bo {[e1 (A - C) - B] Q2 - K3Q + ao (J33 - J11) - a1 + S3} sinu
+ bo [-(e1 + 1)FQ2 - ao J12 + 61] cos u + (EQ2 + J13)bo - EQ2e2 + QeK
- ao (ao J13 + a3 + S1) = 0,
[2b^(FQ2 + J12) - e1 (a2 - h)] cos2 u + {b2 [(B - A)Q2 + Jn - J22] - e1 (a1 + b2)} cosusinu
+ bo((e1 - 1)EQ2 - k1Q + aoJ13 + s1) cosu + bo [(e1 - 1)DQ2 + k2Q - J23ao - s2] sinu
- b2o(FQ2 + J12) - e1b1 + ao(a2 + b1) = 0, (4.1)
whose solution is a set of conditions on the parameters sufficient for the system to admit the desired solution (3.4).
Since recently, Computer Algebra Systems (CAS) have been playing an increasingly important role in mechanics that goes beyond the numerical approximation methods of computational mechanics. In rigid body dynamics, CAS help, for example, in constructing new integrable cases[18], particular solutions, e.g., polynomial solutions [6, 21] and, generally speaking, in verifying solutions obtained via analytical methods.
Using a Computer Algebra System capable of doing symbolic computations, one can prove that (3.4) is a solution of Eqs. (4.1), and hence of the problem (2.4) if
6 (DQ2 - J23) + £163] cos2 u + [62 (EQ2 + J13) - eia3] cos u sinu
a2 = -61 = -£2Q (2 F Q - K12 ), a3 = £2 (E Q2 - QK23 + J13)
62 = [A - e1 (B - C)] Q2 - (£1K22 - aoK33 - k:i)Q + ao (J22 - J33) - S3
63 = -£2 (DQ2 - QK23 + J23)
S1 = 2aoEQ2 + £2K13Q - J13, S2 = 2aoDQ2 + £2K23Q - J23 K1 = 3aoEQ + (1 - 2ao) K13 - £2Q_1 J13, K2 = 3aoDQ + (1 - 2ao) K23 - £2Q-1 J23,
J11 = 3(A - B)Q2 + 2 (K11 - K22) Q + J22, J12 = -Q (3FQ - 2K12),
(4.2a) (4.26) (4.2c) (4.2d) (4.2e) (4.2/) (4.2g)
in which £2 = 1 - ao, and the angular velocity Q is given by
lA2±y/A^-4AiA3 2 Ai
(4.3)
where
A1 = -ao (2A - 3B + C) - A + B + C, A2 = £2Kn + ao (2K22 - K33) - K3, A3 = ao ( J22 - J33) + a1 - S3.
(4.4)
The potential (2.3) is given now by
V = aa — £2^ (2FQ — K12) (a2 — A) + £2 (EQ2 — QK23 + J13) as + + {[A — £1 (B — C)] Q2 — (£1K22 — aoKss — Ks)Q + ao (J22 — J33) — S3} £2 —
- e-2 (DQ2 - Qii23 + J23) «3 + \ [3(A - B)Q2 + 2 (Kn ~ K22) V + J22] 71 - (4.5)
- Q (3FQ - 2A"i2) 7172 + Ji37i73 + 7^227! + 'hijilz + 7^337! +
+ (2aoEQ2 + £2K13Q — J13) 71 + (2aoDQ2 + £2^3Q — J23) 72 + S373.
The solution (4.2) has the free parameters: a0 G [—1,1], a1, 61, k3, c1 and all the values of Kj and Jj except J11 and J22, in addition to the inertia moments of the body. It is remarkable that all nonfree parameters of the system are given explicitly in terms of the free ones.
5. Special cases
5.1. A gyrostat under three uniform fields
Let us consider the case of a gyrostat, with a gyrostatic moment k acted upon by irreducible three uniform fields derived from the potential
V = a ■ a + b ■ 3 + s ■ 7. (5.1)
The motion is governed by the equations
Iu + u x (Iu + k) = a x a + ¡3 x b + 7 x s,
a + u x a = 0, ¡3 + u x (3 = 0, 7 + u x 7 = 0. (5.2)
The problem is the special version of (2.4) when K = J = 0. Assuming in (4.2) zero values of all Kj and Jj, we get from (4.2g) the conditions B = A, F = 0 on the inertia moments. This means that the z-axis is perpendicular to one of the two circular cross-sections of the ellipsoid of inertia of the body at the fixed point and F vanishes for any other pair of axes in that plane. Thus, we can rotate the x — y axes in their plane to satisfy the additional condition D = 0. We end with the inertia tensor
(a 0 —e\
0 A 0 , (5.3)
—E 0 C
and the following conditions for the existence of the precessional motion
a2 = 61 = 63 = S2 = K2 = 0, a3 = £2EQ2, (5.4a)
62 = (aoA — C) Q2 + K3Q — S3, S1 = 2aoEQ2, m = 3aoEQ2. (5.46)
The values of Xi are
A1 = £1C — aoA, X2 = —K3, X3 = a1 — S3. (5.5)
Substituting (5.5) into (4.3), we obtain the angular velocity in the form
n _ 1 K3 ± v/nj + 4 (q.i - s3) (gpA - eiC)
S 2 eiC-aoA " (5"6)
Using conditions (5.4), the potential (2.3) can be brought to the form
V = ai ai + a3 a3 + 62/2 + «i7i + S3Y3 = , x
(5.7)
= i • (aia + siY) + 62 (j • P) + k • (a3a + S37). The third term in (5.7) can be interpreted as the potential of a uniform gravity field
g = a3a + S37, (5.8)
pointing vertically upwards. Then the axis of precession 7 is inclined to the vertical at an angle
c . -iaa 0 = tan — = tan
_-12
K3 ± Vk3 + 4 (ai - S3) (aoA - eiC)
S3 4s3 (eiC - a0A)2
(5.9)
which depends on the angle d0 between the axis of precession 7 and the axis k of proper rotation.
In the special case of a rigid body (k1 = k3 = 0), we see from (5.46) that the precessional motion is possible only when a0 = 0 (d0 = n/2, a'0 = e1 = e2 = 1). Then we get results established recently by Yehia [20].
5.2. The two-field Kowalewski top
The version of the problem when
a = 0, A = 2C, E = 0 (5.10)
describes a Kowalevski gyrostat in two perpendicular fields. In this case, (5.46) gives
62 = C (2a0 - 1)Q2 + K3Q - S3, si = kI = 0, (5.11)
that is, the gyrostatic momentum is directed along the axis of dynamical symmetry. The velocity Q of precession can be obtained from (5.6) as
2e2C
The potential (5.7) is reduced to
V = 62 (j • P) + S3 (k • 7) = 62/32 + S3a0. (5.13)
The first term in V can be interpreted as the potential of a uniform gravity field in the direction of P and then the vector of the body's center of mass is 62j. The second term in V is the potential of a field in the direction of 7 whose center of application is S3k.
6. Concluding Remarks
We add two uniform fields to what is called the generalized problem of the motion of a body [1, 15]. The main result of the present article is that the rigid body in the generic case, which describes a rigid body moving in an ideal incompressible fluid while acted upon by an irreducible combination of two uniform fields, can perform the Grioli-type precessional motion (3.4) under conditions (4.2), (4.3) and (4.4). As far as we know, only the special case
when the body moves by inertia in the fluid, that is, a = b = 0, is considered in the literature, see, for example, [2, 10, 13, 14].
As a special case, we get a solution of the problem of precessional motions of a gyrostat acted upon by irreducible three uniform fields. In the latter case, it is known that the gyrostat can perform a Grioli-type precessional motion when the axis of proper rotation is perpendicular to the axis of precession [20]. We construct here a new solution of the problem assuming an arbitrary value of the angle between the two axes, from which the results of [20] are obtained under the perpendicularity condition. The results are more specialized to prove the existence of the considered class of motion for a Kowalevski gyrostat in two perpendicular uniform fields.
The author thanks Prof. H. M. Yehia for his interest in this work and for stimulating discussion of the results.
References
[1] Gorr, G. V., A New Solution of a Generalized Problem on the Motion of a Body with a Fixed Point, J. Appl. Math. Mech, 1992, vol. 56, no. 3, pp. 443-446; see also: Prikl. Mat. Mekh, 1992, vol. 56, no. 3, pp. 532-535.
[2] Gorr, G. V., Precessional Motions in Rigid Body Dynamics and the Dynamics of Systems of Coupled Rigid Bodies, J. Appl. Math. Mech, 2003, vol.67, no. 4, pp. 511-523; see also: Prikl. Mat. Mekh., 2003, vol. 67, no. 4, pp. 573-587.
[3] Grioli, G., Esistenza e determinazione delle precessioni regolari dinamicamente possibili per un solido pesante asimmetrico, Ann. Mat. Pura Appl. (4), 1947, vol.26, nos. 3-4, pp. 271-281.
[4] Grioli, G., On the General Theory of Asymmetric Gyros, in Kreiselprobleme Gyrodynamics: IUTAM Symp. Celerina, Berlin: Springer, 1963, pp. 26-31.
[5] Gulyaev, M.P., On a New Particular Solution of the Equations of Motion of a Heavy Rigid Body Having a Fixed Point, Vestn. Mosk. Univ., 1955, vol.9, no. 3, pp. 15-21 (Russian).
[6] Hussein, A.M., On the Motion of a Magnetized Rigid Body, Acta Mech., 2017, vol.228, no. 11, pp.4017-4023.
[7] Kharlamova, E.I., On a Particular Solution of the Euler - Poisson Equations, J. Appl. Math. Mech., 1959, vol. 23, pp. 975-988; see also: Prikl. Mat. Mekh., 1959, vol. 23, pp. 681-690.
[8] Kharlamova, E.I., A Particular Case of the Integrability of the Euler - Poisson Equations, Soviet Phys. Dokl, 1959, vol.4, pp. 313-315; see also: Dokl. AN SSSR, 1959, vol.125, no. 5, pp. 996-997.
[9] Kharlamova, E.I., On a Motion of a Rigid Body Having a Fixed Point, Mekh. Tverd. Tela, 1970, vol.2, pp. 35-37 (Russian).
[10] Khlystunova, N. V., Grioli-Type Precession of a Heavy Rigid Body in a Liquid, J. Appl. Math. Mech., 2000, vol. 64, no. 4, pp. 527-530; see also: Prikl. Mat. Mekh., 2000, vol. 64, no. 4, pp. 551-554.
[11] Leimanis, E., The General Problem of the Motion of Coupled Rigid Bodies about a Fixed Point, Springer Tracts in Natural Philosophy, vol. 7, Berlin: Springer, 1965.
[12] Routh, E. J., On a Theorem of Jacobi in Dynamics, Quart. J. Pure Appl. Math., 1889, vol.23, pp. 34-45.
[13] Rubanovskii, V. N., Some Possible Motions of a Heavy Solid in a Fluid, J. Appl. Math. Mech., 1968, vol.32, no.4, pp. 794-800; see also: Prikl. Mat. Mekh., 1968, vol.32, no.4, pp. 763-768.
[14] Rubanovskii, V. N., On a New Particular Solution of a Heavy Solid in a Fluid, J. Appl. Math.. Mech., 1985, vol.49, no. 2, pp. 160-165; see also: Prikl. Mat. Mekh, 1985, vol.49, no. 2, pp. 212-219.
[15] Verkhovod, Ye. V. and Gorr, G. V., New Cases of Isoconic Motions in the Generalized Problem of the Dynamics of a Rigid Body with a Fixed Point, J. Appl. Math.. Mech., 1993, vol. 57, no. 5, pp. 783-792; see also: Prikl. Mat. Mekh., 1993, vol. 57, no. 5, pp. 25-34.
[16] Yehia, H. M., On the Motion of a Rigid Body Acted upon by Potential and Gyroscopic Forces: 2. A New Form of the Equations of Motion of a Multiconnected Rigid Body in an Ideal Incompressible Fluid, J. Mec. Théor. Appl, 1986, vol. 5, no. 5, pp. 755-762.
[17] Yehia, H. M., New Generalizations of the Integrable Problems in Rigid Body Dynamics, J. Phys. A, 1997, vol. 30, no. 20, pp. 7269-7275.
[18] Yehia, H. M., New Generalizations of All the Known Integrable Problems in Rigid Body Dynamics, J. Phys. A, 1999, vol.32, no.43, pp. 7565-7580.
[19] Yehia, H. M., On the Regular Precession of an Asymmetric Rigid Body Acted upon by Uniform Gravity and Magnetic Fields, EJBAS, 2015, vol. 2, no. 3, pp. 200-205.
[20] Yehia, H. M., Regular Precession of a Rigid Body (Gyrostat) Acted upon by an Irreducible Combination of Three Classical Fields, J. Egypt. Math. Soc, 2017, vol. 25, no. 2, pp. 216-219.
[21] Zyza, A. V., Computer Studies of Polynomial Solutions for Gyrostat Dynamics, Comput. Res. Model., 2018, vol. 10, no. 1, pp. 7-25 (Russian).