arxiv: v1 [nlin.si] 31 Dec 2015

Transcription

1 Explicit determination of certain periodic motions of a generalized two-field gyrostat P. E. Ryabov a, S.V. Sokolov b,, I.I. Kharlamova c arxiv: v [nlin.si] 3 Dec 205 a Financial University, Leningradsky Avenue, 49, Moscow, 25993, Russia b Institute of Machines Science, Russian Academy of Sciences, 4 Maly Kharitonyevsky Per., Moscow, 0990, Russia c Russian Academy of National Economy and Public Administration,8 Gagarin Street, Volgograd, 4003, Russia Abstract The case of motion of a generalized two-field gyrostat found by V.V. Sokolov and A. V. Tsiganov is known as a Liouville integrable Hamiltonian system with three degrees of freedom. We find a set of points at which the momentum map has rank. This set consists of special periodic motions which correspond to the singular points of a bifurcation diagram on an iso-energetic surface. For such motions the phase variables can be expressed in terms of algebraic functions of a single auxiliary variable. These algebraic functions satisfy a differential equation integrable in elliptic functions of time. It is shown that the corresponding points in the three-dimensional space of the constants of the integrals belong to the intersection of two sheets of the discriminant surface of the Lax curve. Keywords: completely integrable Hamiltonian systems, spectral curve, special periodic solutions 2000 MSC: 70E7, 70G40, 70H06, 70E40 Contents Intrtoduction 2 2 Discriminant surfaces and explicit integration of certain periodic motions 4 3 Conclusion 4 Acknowledgments This work is partially supported by the grants of RFBR and Volgograd Region Authorities No , , and Corresponding author addresses: orelryabov@mail.ru (P. E. Ryabov, sokolovsv72@mail.ru (S. V. Sokolov, irinah@vags.ru (I. I. Kharlamova Preprint submitted to J. of Geometry and Physics August 7, 208

2 . Intrtoduction The motion of a generalized two-field gyrostat is governed by the following system of differential equations: with Hamiltonian function [] Ṁ = M H H H +α +β M α β, α = α H M, β = β H M, H = M 2 +M2 2 +2M2 3 +2λM 3 2ε 2 (α +β 2 +2ε (M 2 α 3 M 3 α 2 +M 3 β M β 3. ( (2 Here M,α and β stand for the total angular momentum and the intensities of the two forces considered in the moving frame formed by the principal axes of inertia of the body. The gyrostatic momentum is directed along the axis of dynamic symmetry and its axial component is denoted by λ. The parameters ε and ε 2 are called deformation parameters since their zero values define important partial cases and establish relations with some previously known integrable cases. Treating R 9 = {(M,α,β} as the Lie coalgebra e(3,2 we obtain the Lie Poisson bracket {M i,m j } = ε ijk M k, {M i,α j } = ε ijk α k, {M i,β j } = ε ijk β k, {α i,α j } = 0, {α i,β j } = 0, {β i,β j } = 0, (3 ε ijk = 2 (i j(j k(k i, i,j,k 3. With respect to this bracket the system ( can be represented in the Hamiltonian form ẋ = {H,x} where x R 9. Note that the Casimir functions of the bracket (3 are α 2, α β and β 2. Therefore we define the phase space P of system ( as a common level of these functions α 2 = a 2, β 2 = b 2, α β = c, (0 < b < a, c < ab. Using the parametric reduction invented by M. P. Kharlamov [2] we can assume that the parameter c is zero. This simplifies calculations significantly. In [], for the system ( with Hamiltonian function (2, V.V.Sokolov and A. V. Tsiganov gave a Lax representation with a spectral parameter and thereby 2

3 proved the Liouville complete integrability of this system. This Lax representation generalizes the L-A pair for the Kowalevski gyrostat in a double field found by A. G. Reyman and M. A. Semenov-Tian-Shansky [3]. For the Hamiltonian function (2, we represent the additional integrals K and G as functions of two deformation parameters ε and ε 2 [4]: K = Z 2 +Z2 2 λ[(m 3 +λ(m 2 +M2 2 +2ε 2(α 3 M +β 3 M 2 ] +λε 2 (α 2 +β 2 M 3 +2λε [α 2 M 2 β M 2 2 (α β 2 M M 2 ] 2λε 2 ω γ, G = ω 2 α +ω2 β +2(M 3 +λω γ 2ε 2 (α 2 β 2 +β 2 α +2ε [β 2 (M 2 α 3 M 3 α 2 α 2 (M β 3 M 3 β ] +2(α β[ε 2 (α 2 +β +ε (α 3 M α M 3 +β 2 M 3 β 3 M 2 ]. Here we use the following notation: Z = 2 (M2 M2 2 +ε 2(α β 2 +ε [M 3 (α 2 +β M 2 α 3 M β 3 ]+ 2 ε2 (β 2 α 2, Z 2 = M M 2 +ε 2 (α 2 +β ε [M 3 (α β 2 +β 3 M 2 α 3 M ] ε 2 (α β, ω α = α M +α 2 M 2 +α 3 M 3, ω β = β M +β 2 M 2 +β 3 M 3, ω γ = M (α 2 β 3 β 2 α 3 +M 2 (α 3 β α β 3 +M 3 (α β 2 α 2 β. In the special case where ε = 0 and ε 2 =, we get the integrals of motion in the problem of the Kowalevski gyrostat subjected to two homogeneous fields [3, 5]. FortheLaxpairduetoSokolovandTsiganov[], theequationofthespectral curve E(z,ζ reads [4] where E(z,ζ : d 4 ζ 4 +d 2 ζ 2 +d 0 = 0, d 4 = z 4 ε 2 (α 2 +β 2 z 2 ε 4 [α 2 β 2 (α β 2 ], d 2 = 2z 6 +[ε 2 (α2 +β 2 h λ 2 ]z 4 +[ε 2 2 (α2 +β 2 ε 2 g]z2 +2ε 2 ε2 2 [α2 β 2 (α β 2 ], d 0 = z 8 +hz 6 +f ε,ε 2 z 4 +ε 2 2gz 2 ε 4 2[α 2 β 2 (α β 2 ]. The most complicated coefficient f ε,ε 2 of z 4 in d 0 can be expressed in terms of the constants of the integrals h, k, and g as follows: f ε,ε 2 = ε 2 g+k ε4 (α β2 4 [h2 +2ε 2 (α2 +β 2 h+ε 4 (α2 β 2 2 ] ε 2 2 (α2 +β 2. 3

4 We define F : P R 3 by F(x = {g = G(x,k = K(x,h = H(x}. The mapping F is called a momentum mapping. By C we denote the set of all critical points of F, i.e., the set of points x such that rankdf(x < 3. The set of critical values Σ = F(C R 3 is called the bifurcation diagram. Normally, Σ is a stratified 2-manifold. In our case consider the bifurcation diagram Σ R 3 as a two dimensional cell complex. Then the singular points form union of skeletons of dimensions 0 and. The determination of -cells of Σ is much more complicated. The corresponding values of the first integrals are generated by periodic trajectories, i.e., closedorbitsforwhichrankdf(x =. Letuscallsuchatrajectoryaspecial periodic motion (SPM [6]. For the classical Kowalevski top in the gravity field all SPMs are permanent rotations around the vertical axis. For a Kowalevski top in two constant fields (λ = 0,ε = 0 the set of SPMs, as shown in [7], consists of three families of pendulum motions pointed out in [2] for an arbitrary rigid body and the families of critical periodic motions of the Bogoyavlensky case [8]. These latter motions were first described in [9] and explicitly integrated in [0]. Note that pendulum motions were first found by Yehia [] with no conditions imposed on the moments of inertia but under some special restrictions on the location of the centers of application of the fields. For the integrable system of Kowalevski Yehia SPMs were presented in [2],[3]. Similar investigations in hydrodynamics were performed in [4], [5] and [6]. The aim ofthe presentarticleisthe constructionofacertainclassofperiodic motions for the generalized two-field gyrostat. To do this we will study the singularities of the spectral curve associated with the Sokolov Tsiganov L-A pair. 2. Discriminant surfaces and explicit integration of certain periodic motions The discriminant surface of the spectral curve E(z, ζ for Sokolov-Tsiganov L A pair was found in [7] and consists of two surfaces Π Π 2 : Π : g(t = ht2 2t 3 ε ε2 2(a 2 b 2 c 2, t k(t = 3t 2 2ht+ ε2 (2t3 ht 2 ε ε2 2 (c2 a 2 b 2 (2ε 2 t+ε2 2 t {h2 +2ε 2 (a 2 +b 2 h+ε 4 [(a 2 b c 2 ]+4ε 2 2(a 2 +b 2 }. 4

5 and Π 2 : g(s = ε4 λ2 [(a 2 b c 2 ] s s 3 {ε2 [ε2 (a2 +b 2 +h+λ 2 ]+2ε 2 2 } 8ε 8 λ2 s 2 + 6λ 2 ε {ε2 [(a 2 b c 2 ]+(a 2 +b 2 (h+λ 2 }, k(s = ε8 λ4 [(a 2 b c 2 ] s 2 + {ε2 [ε2 (a2 +b 2 +h+λ 2 ]+2ε 2 2 } 2ε 4 s 3 s 2 6ε 4 2 λ2[ 2ε 2 (a 2 +b 2 +h+ λ2 ]. 2 Consider a typical example of a discriminant surface on an iso-energetic level h = const. Q 2 Q P P 2 k g Figure : Discriminant surfaces in the section h = const for parameters a =.;c = 0;ε 2 = 0.26;ε = 0.85;b = 0.78;λ =.958;h = 3.4. Q 2 P k g Q Figure 2: Enlarged fragment. The points P,P 2,Q and Q 2 are the self-intersection points of the discriminant leafs Π and Π 2. The coordinates of these points on the iso-energetic level 5

6 h = const are determined by the formulas g = b 2 (h+λ 2 ε P : 2 (a 2 b 2 (ε 2 2 +b 2 ε 4, k = 4 [ε2 (a2 +b 2 +h] 2 + ε2 2 ε 2 (h+λ 2 + ε4 2 ε 4 +ε 2 2 (a2 +b 2 ε 2 b2 λ 2, g = a 2 (h+λ 2 + ε P 2 : 2 (a 2 b 2 (ε 2 2 +a2 ε 4, k = 4 [ε2 (a2 +b 2 +h] 2 + ε2 2 ε 2 (h+λ 2 + ε4 2 ε 4 +ε 2 2 (a2 +b 2 ε 2 a2 λ 2, { g = ±abh, Q,2 : k = 4 (a b2 [2ε 2 h+ε4 (a±b2 +4ε 2 2 ]. Non-degenerate singularities of rank or special periodic motions are the preimages of the points mentioned above. How can one find them? By the definition of such a singularity [8] it is required to find two function g (x and g 2 (x for which at the points x 0 of rank the conditions dg (x 0 = dg 2 (x 0 = 0 are fulfilled. The formulas for the coordinates of the points P k and Q k suggest that ( P : g = G b 2 H, g 2 = K 2 h+ [(a 2 +b 2 ε 4 ε 2 +2ε2 2 H, ] ( P 2 : g = G a 2 H, g 2 = K 2 h+ [(a 2 +b 2 ε 4 ε 2 +2ε2 2 H, ] Q,2 : g = G abh, g 2 = K 2 (a b2 ε 2 H. At the singularities of rank (for which P is the image the following relation must hold: G b 2 H+A (α 2 +α 2 2+α 2 3+B (β 2 +β 2 2+β 2 3+C (α β +α 2 β 2 +α 3 β 3 = 0. (4 Here A, B and C are undetermined multipliers. Equations (4 are equivalent to sgradg = b 2 sgradh. 6

7 In this case the parametrization of special periodic motions looks like b2 ε 4 M = ε2 2 [b(+ξ2 (+η 2 +2aξ( η 2 ] 2λε b(+η 2 ξ ε b( ξ 2 (+η 2, b2 ε 4 M 2 = ε ε2 2 [a(+ξ2 ( η 2 +2bξ(+η 2 ] λε b(+ξ 2 (+η 2 2 ε b, b 2 ε 4 ε2 2 ( ξ2 (+η 2 M 3 = a b 2 ε 4 ε2 2 ( η2 λε b(+η 2 ε b(+η 2, α = 2a ε 2( η 2 ξ bε 2 ( ξ 2 η bε 2 (+η2 (+ξ 2, α 2 = a b2 ε 4 ε2 2 ( η2 bε 2, (+η2 α 3 = a ε 2( η 2 ( ξ 2 +4bε 2 ηξ bε 2 (+η2 (+ξ 2, β = 2 b 2 ε 4 ε2 2 ξ ε 2, β 2 = ε 2 b2 ε 4 (+ξ2 ε 2, β 3 = ε2 2 ( ξ2 ε 2. (+ξ2 (5 Here the auxiliary variables ξ and η satisfy the system of differential equations b2 ε 4 η = ε ε2 2 [a(+ξ2 ( η 2 +2bξ(+η 2 ] λε b(+ξ 2 (+η 2, b2 ε 4 ε2 2 ( ξ2 ξ = a ε 2( η 2 ( ξ 2 +4ε 2 bξη bε (+η 2. Using the equation H = h one can eliminate η and thus obtain a single equation in ξ(t: ( ξ 2 = ε 2 (a 4 ξ 4 +a 3 ξ 3 +a 2 ξ 2 +a ξ +a 0, ε2 2 where a 4 = ε 2 2( ε 2 2ε 2 λ 2 +4q 2 ε q 4 +2rq +ε 2 q 2 h, a 3 = 4ε 2 2 ε λ(qε 2 2 +q3 +r, a 2 = 4ε 2 q 4 h 8q 2 ε r 2 8q 6 8q 4 ε 2 2 2ε 2 q 2 ε 2 2h 2ε 2 ε 4 2λ 2, a = 4ε 2 2 ε λ(qε 2 2 +q3 r, a 0 = ε 2 2( ε 2 2ε 2 λ 2 +4q 2 ε 2 2 2rq +3q 4 +ε 2 q 2 h. Here the parameters q and r can be expressed in terms of a and b by formulas q = b 2 ε 4 ε2 2, r = ε 2 (b2 ε 4 ε2 2 [ε2 b2 (h+2ε 2 b2 +ε 2 2 (a2 +b 2 ] λ 2 b 2 ε 2 ε2 2. 7

8 These formulas can be used for construction of the characteristic polynomial which allows determination of the type of a singularity of rank. In the generic case the characteristic polynomial reads: (µ = µ 4 p 2 µ 2 p 4, where p 2 = 2 tracea2 g, p 4 = 4 [tracea4 g 2 (tracea2 g 2 ]. Here A g denotes a symplectic operator which is linearization of the vector field sgradg at the points of solution (5. Calculation of the coefficients p 2 and p 4 for the function g at point P yields p 2 = 4 ε 2 (a 2 b 2 [a 2 b 2 ε 4 2ε2 2 a2 3b 4 ε 4 ε2 b2 (h+λ 2 ], p 4 = 6 ε 4 (a 2 b 2 3 [(b 2 ε 4 ε2 2 (ε2 b2 (h+2b 2 +ε 2 2 (a2 +b 2 ε 2 ε2 2 λ2 b 2 ]. Now let us describe periodic motions for the point P 2. Given that the skewsymmetric gradients are linearly dependent sgradg = a 2 sgradh one obtains the following parametrization a2 ε 4 M = ε ε2 2 [b(+ξ2 ( η 2 +2aξ(+η 2 ]+λε a(+ξ 2 (+η 2 2 ε a, a 2 ε 4 ε2 2 ( ξ2 (+η 2 a2 ε 4 M 2 = ε2 2 [a(+ξ2 (+η 2 +2bξ( η 2 ]+2λε a(+η 2 ξ ε a( ξ 2 (+η 2, M 3 = b a 2 ε 4 ε2 2 ( η2 +λε a(+η 2 ε a(+η 2, α = ε 2 ε 2, α 2 = 2 a 2 ε 4 ε2 2 ξ a2 ε 4 ε 2, α 3 = ε2 2 ( ξ2 (+ξ2 ε 2, (+ξ2 β = b a 2 ε 4 ε2 2 ( η2 aε 2 (+η2 β 3 = b ε 2( ξ 2 ( η 2 +4aε 2 ξη aε 2 (+ξ2 (+η 2., β 2 = 2b ε 2( η 2 ξ aε 2 ( ξ2 η aε 2 (+ξ2 (+η 2, (6 8

9 Here the auxiliary variables ξ and η satisfy the equations a2 ε 4 η = ε ε2 2 [b(+ξ2 ( η 2 +2aξ(+η 2 ]+λε a(+ξ 2 (+η 2, a2 ε 4 ε2 2 ( ξ2 ξ = b ε 2( η 2 ( ξ 2 +4ε 2 aξη aε (+η 2. Using the equation H = h one can eliminate η and thus obtain a single equation in ξ(t ( ξ 2 = ε 2 ε2 2 m2(b 4ξ 4 +b 3 ξ 3 +b 2 ξ 2 +b ξ +b 0, where b 4 = ε 2 2(4ε 2 2m 2 ε 2 λ 2 ε 2 2 +hε 2 m 2 2nm+3m 4, b 3 = 4ε 2 2 ε λ(ε 2 2 m n+m3, b 2 = 2ε 4 2ε 2 λ 2 4ε 2 hm 4 8ε 4 2m 2 8ε 2 2m 4 +4n 2 8m 6 2ε 2 hε 2 2m 2, b = 4ε 2 2 ε λ(ε 2 2 m+n+m3, b 0 = ε 2 2(4ε 2 2m 2 ε 2 λ 2 ε 2 2 +hε 2 m 2 +3m 4 +2nm. Here the parameters m and n can be expressed in terms of a and b as follows m = a 2 ε 4 ε2 2, n = ε 2 (a2 ε 4 ε2 2 [ε2 a2 (h+2ε 2 a2 +ε 2 2 (a2 +b 2 ] λ 2 a 2 ε 2 ε2 2. Finally the type of a singularity of rank which corresponds to the point P 2 can be determined with the help of the characteristic equation here (µ = µ 4 p 2 µ2 p 4, p 2 = 4 ε 2 (a 2 b 2 [a 2 b 2 ε 4 2ε2 2 b2 3a 4 ε 4 ε2 a2 (h+λ 2 ], p 4 = 6 ε 4 (a 2 b 2 3 [(a 2 ε 4 ε 2 2(ε 2 a 2 (h+2a 2 +ε 2 2(a 2 +b 2 ε 2 ε 2 2λ 2 a 2 ]. At the point Q we obtain a solution of the form M = M 2 = α 3 = β 3 = 0, α = acos(ϕ,α 2 = asin(ϕ, (7 β = bsin(ϕ,β 2 = bcos(ϕ. ϕ = 8ε 2 (a+bsin(ϕ+4ε (a+b[ε (a+bsin(ϕ λ]cos(ϕ. 9

10 Integrating gives ( ϕ 2 = 4[ε 2 (a+b 2 sin 2 ϕ+4ε 2 (a+bcosϕ 2λε (a+bsinϕ+2h+λ 2 ]. Introduce a new variable x(t: Then x = tan ϕ 2. cosϕ = x2 2x 2ẋ +x2, sinϕ = +x2, ϕ = +x 2 The differential equation in x = x(t looks like (ẋ 2 = c 4 x 4 +c 3 x 3 +c 2 x 2 +c x+c 0, here c 4 = 2h+λ 2 4ε 2 (a+b, c 3 = 4λε (a+b, c 2 = 2[2ε 2 (a+b2 +2h+λ 2 ], c = 4λε (a+b, c 0 = λ 2 +2h+4ε 2 (a+b. Therefore x(t can be written in terms of elliptic quadratures. The characteristic equation for determination of the type of the singularity at the point Q reads µ 4 +d 2 µ+d 0 = 0, where d 2 = 2ab[(a b 2 (h+4ε 2 ab 4λ2 ab], d 0 = 8a 3 b 3 [2ε 2 2(a b 4 + ε 2 (a b 2 ((a b 2 (h+2ε 2 ab 4λ 2 ab λ 2 (h(a b 2 2λ 2 ab]. Note that the negative sign of the value d 2 2 4d 0 = 4a 2 b 2 (a b 4 (h 2 6ε 2 2ab indicates that the singularity of rank is a focus. At the point Q 2 the solution is M = M 2 = α 3 = β 3 = 0, α = acos(ϕ,α 2 = asin(ϕ, (8 β = bsin(ϕ,β 2 = bcos(ϕ. 0

11 ( ϕ 2 = 4[ε 2 (a b2 sin 2 ϕ 4ε 2 (a bcosϕ 2λε (a bsinϕ+2h+λ 2 ]. This equation can be solved in terms of elliptic quadratures. The type of the singularity that corresponds to the point Q 2 can be determined from the characteristic equation where µ 4 +d 2 µ+d 0 = 0, d 2 = 2ab[(a+b2 (h 4ε 2 ab+4λ2 ab], d 0 = 8a 3 b 3 [2ε 2 2(a+b 4 + ε 2 (a+b 2 ((a+b 2 (h 2ε 2 ab+4λ 2 ab λ 2 (h(a+b 2 +2λ 2 ab]. Here the value of the expression is always non-negative. d 22 4d 0 = 4a 2 b 2 (a+b 4 (h 2 +6ε 2 2 ab 3. Conclusion The constructed periodic solutions (5, (6, (7 and (8 are singularities of rank of the momentum map. The image of these solutions under the momentum map are singular points in the bifurcation diagram. These solutions play a key role in the construction of the atlas of the bifurcation diagram for the generalized two-field gyrostat. A similar study of the atlas of a bifurcation diagram was carried out by M.P. Kharlamov for the motion of a Kowalevski top in the double field of forces. 4. Acknowledgments The authors are grateful to A. V. Bolsinov for extremely valuable discussions and guidance, National Science Foundation of Poland, Stefan Banach International Mathematical Center, Warsaw Center of Mathematics and Computer Science, Faculty of Physics and Astronomy University of Zielona Gora and the Organizers of FDIS-205 especially V. S. Matveev, S. L. Tabachnikov, A. Maciejewski for hospitality and support. References [] V. V. Sokolov, A. V. Tsiganov, Lax pairs for the deformed Kowalevski and Goryachev Chaplygin tops, Theoretical and Mathematical Physics 3 ( ( doi:0.023/a:

12 [2] M. P. Kharlamov, Bifurcation diagrams of the Kowalevski top in two constant fields, Regular and Chaotic Dynamics 0 (4 ( arxiv: , doi:0.070/rd2005v00n04abeh [3] A. G. Reyman, M. A. Semenov-Tian-Shansky, Lax representation with a spectral parameter for the Kowalewski top and its generalizations, Lett. Math. Phys. 4 ( ( doi:0.007/bf [4] P. E. Ryabov, Phase topology of one irreducible integrable problem in the dynamics of a rigid body, Theoretical and Mathematical Physics 76 (2 ( doi:0.007/s [5] A. I. Bobenko, A. G. Reyman, M. A. Semenov-Tian-Shansky, The Kowalewski top 99 years later: a Lax pair, generalizations and explicit solutions, Commun. Math. Phys. 22 (2 ( doi:0.007/bf [6] M. P. Kharlamov, Periodic motions of the Kowalevski gyrostat in two constant fields, Journal of Physics A: Mathematical and Theoretical 4 (27 ( doi:0.088/75-83/4/27/ [7] M. P. Kharlamov, Singular periodic solutions in the generalized Delone case, Mekh. Tverd. Tela 36 ( [8] O. I. Bogoyavlensky, Euler equations on finite-dimension Lie algebras arising in physical problems, Commun. Math. Phys. 95 (3 ( doi:0.007/bf [9] D. B. Zotev, Fomenko-Zieschang invariant in the Bogoyavlenskyi case, Regular and Chaotic Dynamics 5 (4 ( doi:0.070/rd2000v004n04abeh [0] M. P. Kharlamov, Regions of existence of critical motions for the generalized Kowalevski top and bifurcation diagrams, Mekh. Tverd. Tela 36( [] H. M. Yehia, On certain integrable motions of a rigid body acted upon by gravity and magnetic fields, Int. J. Non-Linear Mech. 36 ( [2] P. V. Kharlamov, One case of integrability of the equations of the motion of a rigid body having a fixed point, Mekh. Tverd. Tela 3 ( [3] E. I. Kharlamova, P. V. Kharlamov, New solution of the differential equations of the motion of a body having a fixed point under the conditions of S.V. Kovalevskaya, Mekh. Tverd. Tela 89 (5 ( [4] S. V. Sokolov, S. M. Ramodanov, Falling motion of a circular cylinder interacting dynamically with a point vortex, Regular and Chaotic Dynamics 8 ( ( doi:0.34/s

13 [5] S. Bezglasnyi, Stabilization of stationary motions of a gyrostat with a cavity filled with viscous fluid, Russian Aeronautics (Iz VUZ 57 (4 ( doi:0.303/s [6] A. V. Borisov, V. G. Lebedev, Dynamics of three vorteces on a plane and a sphere ii. General compact case, Regular and Chaotic Dynamics 3 (2 ( doi:0.070/rd998v003n02abeh [7] P. E. Ryabov, New invariant relations for the generalized twofield gyrostat, Journal of Geometry and Physics 87 ( doi:0.06/j.geomphys [8] A. V. Bolsinov, A. T. Fomenko, Integrable Hamiltonian Systems: Geometry, Topology, Classification, Chapman & Hall/CRC,