The periodic rotary motions of a rigid body in a new domain of angular velocity

Abstract

In the previous works, the limiting case for the motion of a rigid body about a fixed point in a Newtonian force field, which comes from a gravity center lies on Z-axis, is solved. The authors apply the small parameter technique which is achieved giving the body a sufficiently large angular velocity component r_o about the fixed z-axis of the body. The periodic solutions of motion are obtained in neighborhood r_o tends to $\infty$. In our work, we aim to find periodic solutions to the problem of motion in the neighborhood of r₀ tends to $0$. So, we give a new assumption that: r_o is sufficiently small. Under this assumption, we must achieve a large parameter and search for another technique for solving this problem. This technique is named; a large parameter technique instead of the small one well known previously. We see the advantage of the new technique which appears in saving high energy used to begin the motion and give the solution of the problem in another domain. The obtained solutions by the new technique depend on r_o. We consider that the center of mass of this body does not necessarily coincide with the fixed point O. We reduce the six nonlinear differential equations of the body and their three first integrals to a quasilinear autonomous system of two degrees of freedom and one first integral. We solve the rational case when the frequencies of the generating system are rational except $(\,\omega = \,1,\,2,1/2,3,1/3, \ldots )$ under the condition $\gamma^{\prime\prime}_{0} = \cos \theta_{o} \approx 0$. We use the fourth-order Runge–Kutta method to find the periodic solutions in the closed interval of the time t and to compare the analytical method with the numerical one.

Introduction

Some asymptotic perturbed techniques [1,2,3] are widely used by many authors for solving the ordinary linear and nonlinear systems for differential equations in different problems of engineering, mathematical physics, and astronomy. As an extension of this type of problem, we use some perturbation and numerical techniques in the movement of coherent bodies around a fixed point in the presence of new conditions to obtain periodic solutions of different scales to those obtained before.

As in [4], let the fixed z-axis of the ellipsoid of inertia of the body makes an angle $\theta_{o} \approx {\raise0.7ex\hbox{$\pi $} \!\mathord{\left/ {\vphantom {\pi 2}}\right.\kern-0pt} \!\lower0.7ex\hbox{$2$}}$ with the fixed Z-axis in space. We assume the new value of r_o, which is sufficiently small instead of sufficiently large value in [4]. We define a large parameter $\mu$ proportional to 1/r_o instead of the small one in [4]. Consider A, B, C represent the moments of inertia of the body, p, q, r are the components of the angular velocity vector and $\gamma ,\gamma^{\prime},\gamma^{\prime\prime}$ are the direction cosines of the unit vector in direction of the Z-axis. So, the equations of motion and their three first integrals are derived in the form:

$$\begin{aligned} & \dot{p}_{1} + A_{1} q_{1} r_{1} = \mu^{ - 1} a^{ - 1} (\gamma^{\prime\prime}y^{\prime}_{0} - \gamma^{\prime}z^{\prime}_{0} + kaA_{1} \gamma^{\prime}\gamma^{\prime\prime}), \\ & \dot{q}_{1} + B_{1} p_{1} r_{1} = \mu^{ - 1} b^{ - 1} (\gamma z^{\prime}_{0} - \gamma^{\prime\prime}x^{\prime}_{0} + kbB_{1} \gamma^{\prime\prime}\gamma ), \\ & \dot{r}_{1} = \mu^{ - 2} ( - C_{1} p_{1} q_{1} + \gamma^{\prime}x^{\prime}_{0} - \gamma y^{\prime}_{0} + kC_{1} \gamma \gamma^{\prime}), \\ & \dot{\gamma } = r_{1} \gamma^{\prime} - \mu^{ - 1} q_{1} \gamma^{\prime\prime},\quad \dot{\gamma }^{\prime} = \,\mu^{ - 1} p_{1} \gamma^{\prime\prime} - \,r_{1} \gamma ,\quad \dot{\gamma}^{\prime\prime} = \mu^{ - 1} (q_{1} \gamma - p_{1} \gamma^{\prime}).\\ \end{aligned}$$

(1)

$$r_{1}^{2} = 1 + \mu^{ - 2} s_{1} ,\quad r_{1} \gamma^{\prime\prime} = \gamma^{\prime\prime}_{0} + \mu^{ - 1} s_{2} ,\quad \gamma^{2} + \gamma^{\prime 2} + \gamma^{\prime \prime 2} = 1.$$

(2)

Such that:

$$\begin{aligned} & s_{1} = s_{3} - 2z^{\prime}_{0} (\gamma^{\prime\prime}_{0} - \gamma^{\prime\prime}) + k(\gamma_{0}^{\prime \prime 2} - \gamma^{\prime \prime 2} ), \\ & s_{2} = a(p_{10} \gamma_{0} - p_{1} \gamma ) + b(q_{10} \gamma^{\prime}_{0} - q_{1} \gamma^{\prime}) \\ & s_{3} = a(p_{10}^{2} - p_{1}^{2} ) + b(q_{10}^{2} - q_{1}^{2} ) - 2[x^{\prime}_{0} (\gamma_{0} - \gamma ) + y^{\prime}_{0} (\gamma^{\prime}_{0} - \gamma^{\prime})] \\ & \quad + k[a(\gamma_{0}^{2} - \gamma^{2} ) + b(\gamma_{0}^{\prime 2} - \gamma^{\prime 2} )]. \\ \end{aligned}$$

(3)

$$p_{1} = {p \mathord{\left/ {\vphantom {p c}} \right. \kern-0pt} c},\quad (pqr),\quad \gamma_{0} > 0,\quad \mu = {c \mathord{\left/ {\vphantom {c {r_{0} ,}}} \right. \kern-0pt} {r_{0} ,}}\quad l^{2} = x_{0}^{2} + y_{0}^{2} + z_{0}^{2} ,$$

(4)

$$k = {N \mathord{\left/ {\vphantom {N {c^{2} ,}}} \right. \kern-0pt} {c^{2} ,}}\quad A_{1} = (C - B)/A,\quad (ABC).$$

(5)

Reduction of the equations of motion to a quasilinear autonomous system

Solving the first and the second equations of (2) for $r_{1}$ and $\gamma^{\prime\prime}$, we get:

$$r_{1} = 1 + \frac{1}{2}\mu^{ - 2} s_{3} + \ldots ,\quad \gamma^{\prime\prime} = \gamma^{\prime\prime}_{0} + \mu^{ - 1} s_{2} - \frac{1}{2}\mu^{ - 2} s_{3} \gamma^{\prime\prime}_{0} + \ldots .$$

(6)

Differentiating the first and the fourth equations of (1) and using (6) for reduction in the four remaining equations into two differential equations of the second order, we get:

$$\begin{aligned} & \ddot{p}_{1} + \omega^{2} p_{1} = \mu^{ - 1} [z^{\prime}_{0} (a^{ - 1} - A_{1} b^{ - 1} )\gamma + A_{1} b^{ - 1} x^{\prime}_{0} \gamma^{\prime\prime}_{0} - k(A_{1} - \omega^{2} )\gamma^{\prime\prime}_{0} \gamma ] \\ & \quad + \,\mu^{ - 2} \{ A_{1} x^{\prime}_{0} (b^{ - 1} s_{2} - q_{1} \gamma^{\prime}) - \omega^{2} p_{1} s_{3} + p_{1} (A_{1} C_{1} q_{1}^{2} - a^{ - 1} z^{\prime}_{0} \gamma^{\prime\prime}_{0} ) \\ & \quad + \,y^{\prime}_{0} [(2A_{1} + a^{ - 1} b)q_{1} \gamma - a^{ - 1} p_{1} \gamma^{\prime}] + kA_{1} [p_{1} (\gamma_{0}^{\prime \prime 2} - \gamma^{\prime 2} ) \\ & \quad + \,q_{1} (1 - C_{1} )\gamma \gamma^{\prime} - s_{2} (1 + B_{1} )\gamma ]\} + \mu^{ - 3} \{ 0.5(a^{ - 1} - A_{1} b^{ - 1} )z^{\prime}_{0} s_{3} \gamma \\ & \quad - \,s_{2} p_{1} [z^{\prime}_{0} (2\omega^{2} + a^{ - 1} ) - 2k(A_{1} + \omega^{2} )\gamma^{\prime\prime}_{0} ]\} + \ldots , \\ & \ddot{\gamma } + \gamma = \mu^{ - 1} \gamma^{\prime\prime}_{0} p_{1} (1 + B_{1} ) + \mu^{ - 2} \{ p_{1} [(1 + B_{1} )s_{2} + (1 - C_{1} )q_{1} \gamma^{\prime}] \\ & \quad - \,\gamma [s_{3} + y^{\prime}_{0} \gamma^{\prime} + z^{\prime}_{0} b^{ - 1} \gamma^{\prime\prime}_{0} + q_{1}^{2} + k(C_{1} \gamma^{\prime 2} - B_{1} \gamma_{0}^{\prime \prime 2} )] \\ & \quad + \,x^{\prime}_{0} (b^{ - 1} \gamma_{0}^{\prime \prime 2} + \gamma^{\prime 2} )\} + \mu^{ - 3} \{ 2b^{ - 1} x^{\prime}_{0} \gamma^{\prime\prime}_{0} s_{2} \\ & \quad + \,[2k(1 - B_{1} )\gamma^{\prime\prime}_{0} - z^{\prime}_{0} (2 + b^{ - 1} )]s_{2} \gamma_{2} \} + \ldots . \\ \end{aligned}$$

(7)

Here:

$$\omega^{2} = - A_{1} B_{1} = (A - C)(B - C)/AB = (a - 1)(b - 1)/ab.$$

Solving the first and the fourth equations of the system (1) and using (6), we get:

$$\begin{aligned} & q_{1} = - \dot{p}_{1} A_{1}^{ - 1} + \mu^{ - 1} \,a^{ - 1} A_{1}^{ - 1} (y^{\prime}_{0} \gamma^{\prime\prime}_{0} - z^{\prime}_{0} \gamma^{\prime} + a\,k\,A_{1} \gamma^{\prime}\gamma^{\prime\prime}_{0} ) + \ldots , \\ & \gamma^{\prime} = \dot{\gamma } - \mu^{ - 1} A_{1}^{ - 1} \gamma^{\prime\prime}_{0} \dot{p}_{1} + \ldots . \\ \end{aligned}$$

(8)

Introducing the new variables $p_{2}$ and $\gamma_{2}$ such that:

$$\begin{aligned} & p_{1} = p_{2} + \mu^{ - 1} \chi_{1} \gamma_{2} + \mu^{ - 1} \frac{{x^{\prime}_{0} }}{1 - a}\gamma^{\prime\prime}_{0} , \\ & \gamma = \gamma_{2} + \mu^{ - 1} a\gamma^{\prime\prime}_{0} p_{2} , \\ \end{aligned}$$

(9)

where:

$$\chi_{1} = (1 - \omega^{2} )^{ - 1} [ - z^{\prime}_{0} (a^{ - 1} - A_{1} \,b^{ - 1} ) + k\gamma^{\prime\prime}_{0} (A_{1} - \omega^{2} )]\,.$$

Substituting from (9) into (8), we get:

$$\begin{aligned} & q_{1} = - A_{1}^{ - 1} \dot{p}_{2} + \mu^{ - 1} A_{1}^{ - 1} (y^{\prime}_{0} a^{ - 1} \gamma^{\prime\prime}_{0} - \chi_{2} \dot{\gamma }_{2} ) + \ldots , \\ & \gamma^{\prime} = \dot{\gamma }_{2} - \mu^{ - 1} \gamma^{\prime\prime}_{0} A_{1}^{ - 1} b\dot{p}_{2} + \ldots ,\quad \chi_{2} = \chi_{1} + z^{\prime}_{0} a^{ - 1} - kA_{1} \gamma^{\prime\prime}_{0} . \\ \end{aligned}$$

(10)

Substituting (10) and (9) into (3), we get:

$$s_{i} = s_{i1} + 2^{2 - i} \mu^{ - 1} s_{i2} + \ldots ,\quad i = 1,2,$$

(11)

where:

$$\begin{aligned} & {{s_{11} = a(p_{20}^{2} - p_{2}^{2} ) + b(\dot{p}_{20}^{2} - \dot{p}_{2}^{2} )} \mathord{\left/ {\vphantom {{s_{11} = a(p_{20}^{2} - p_{2}^{2} ) + b(\dot{p}_{20}^{2} - \dot{p}_{2}^{2} )} {A_{1}^{2} - 2[x^{\prime}_{0} (\gamma_{20} - }}} \right. \kern-0pt} {A_{1}^{2} - 2[x^{\prime}_{0} (\gamma_{20} - }}\gamma_{2} ) \\ & \quad + \,y^{\prime}_{0} (\dot{\gamma }_{20} - \dot{\gamma }_{2} )] + k[a(\gamma_{20}^{2} - \gamma_{2}^{2} ) + b(\dot{\gamma }_{20}^{2} - \dot{\gamma }_{2}^{2} )], \\ & s_{12} = a\left[ {\frac{{x^{\prime}_{0} }}{1 - a}\gamma^{\prime\prime}_{0} (p_{20} - p_{2} ) + \chi_{1} (p_{20} \gamma_{20} - p_{2} \gamma_{2} )} \right] \\ & \quad - \,bA_{1}^{ - 2} [y^{\prime}_{0} a^{ - 1} \gamma^{\prime\prime}_{0} (\dot{p}_{20} - \dot{p}_{2} ) - \chi_{2} (\dot{p}_{20} \dot{\gamma }_{20} - \dot{p}_{2} \dot{\gamma }_{2} )] \\ & \quad - \,ax^{\prime}_{0} \gamma^{\prime\prime}_{0} (p_{20} - p_{2} ) + y^{\prime}_{0} bA_{1}^{ - 1} \gamma^{\prime\prime}_{0} (\dot{p}_{20} - \dot{p}_{2} ) + (z^{\prime}_{0} - k\gamma^{\prime\prime}_{0} )s_{21} \\ & \quad + \,k\gamma^{\prime\prime}_{0} [a^{2} (p_{20} \gamma_{20} - p_{2} \gamma_{2} ) - b^{2} A_{1}^{ - 1} (\dot{p}_{20} \dot{\gamma }_{20} - \dot{p}_{2} \dot{\gamma }_{2} )], \\ & s_{21} = a(p_{20} \gamma_{20} - p_{2} \gamma_{2} ) - bA_{1}^{ - 1} (\dot{p}_{20} \dot{\gamma }_{20} - \dot{p}_{2} \dot{\gamma }_{2} ), \\ & s_{22} = a\left[ {a\gamma^{\prime\prime}_{0} (p_{20}^{2} - p_{2}^{2} ) + \frac{{x^{\prime}_{0} }}{1 - a}\gamma^{\prime\prime}_{0} (\gamma_{20} - \gamma_{2} ) + \chi_{1} (\gamma_{20}^{2} - \gamma_{2}^{2} )} \right] \\ & \quad + \,bA_{1}^{ - 1} [bA_{1}^{ - 1} \gamma^{\prime\prime}_{0} (\dot{p}_{20}^{2} - \dot{p}_{2}^{2} ) + y^{\prime}_{0} a^{ - 1} \gamma^{\prime\prime}_{0} (\dot{\gamma }_{20} - \dot{\gamma }_{2} ) - \chi_{2} (\dot{\gamma }_{20}^{2} - \dot{\gamma }_{2}^{2} )], \\ & s_{3} = s_{11} + 2\mu^{ - 1} (s_{12} - z^{\prime}_{0} s_{21} + ks_{21} \gamma^{\prime\prime}_{0} ) + \ldots . \\ \end{aligned}$$

(12)

Substituting (11), and (12) into (6) yields:

$$r_{1} = 1 + \frac{1}{2}\mu^{ - 2} s_{11} + \ldots ,\quad \gamma^{\prime\prime} = \gamma^{\prime\prime}_{0} + \mu^{ - 1} s_{21} + \mu^{ - 2} s_{22} - \frac{1}{2}\mu^{ - 2} s_{11} \gamma^{\prime\prime}_{0} + \ldots .$$

(13)

Substituting (9), (10), (11), and (12) into (7), we obtain a quasilinear autonomous system [5]:

$$\ddot{p}_{2} + \omega^{2} p_{2} = \mu^{ - 2} F(p_{2} ,\dot{p}_{2} ,\gamma_{2} ,\dot{\gamma }_{2} ,\varepsilon ),\quad \ddot{\gamma }_{2} + \gamma_{2} = \mu^{ - 2} \Phi (p_{2} ,\dot{p}_{2} ,\gamma_{2} ,\dot{\gamma }_{2} ,\varepsilon ),$$

(14)

where:

$$\begin{array}{*{20}l} {F = F_{2} + \mu^{ - 1} F_{3} + \ldots ,} \hfill & {\varPhi = \varPhi_{2} + \mu^{ - 1} \varPhi_{3} + \ldots ,} \hfill \\ {F_{2} = f_{2} - a\gamma^{\prime\prime}_{0} \chi_{1} (1 - \omega^{2} )p_{2} ,} \hfill & {\varPhi_{2} = \varphi_{2} + a\gamma^{\prime\prime}_{0} (1 - \omega^{2} )[\chi_{1} \gamma_{2} + \gamma^{\prime\prime}_{0} x^{\prime}_{0} (1 - a)^{ - 1} ],} \hfill \\ \end{array}$$

$$\begin{aligned} & F_{3} = f_{3} - \chi_{1} \{ \varphi_{2} + a\gamma^{\prime\prime}_{0} (1 - \omega^{2} )[\chi_{1} \gamma_{2} + \gamma^{\prime\prime}_{0} x^{\prime}_{0} (1 - a)^{ - 1} ]\} , \\ &\Phi _{3} = \varphi_{3} + a\gamma^{\prime\prime}_{0} [a\gamma^{\prime\prime}_{0} \chi_{1} (1 - \omega^{2} )p_{2} - f_{2} ], \\ & f_{2} = x^{\prime}_{0} (A_{1} b^{ - 1} s_{21} + \dot{p}_{2} \dot{\gamma }_{2} ) + A_{1}^{ - 1} \dot{p}_{2} [C_{1} p_{2} \dot{p}_{2} - y^{\prime}_{0} (2A_{1} + ba^{ - 1} )\gamma_{2} ] - p_{2} (\omega^{2} s_{11} \\ & \quad + \,y^{\prime}_{0} a^{ - 1} \dot{\gamma }_{2} + z^{\prime}_{0} a^{ - 1} \gamma^{\prime\prime}_{0} ) - k\{ \gamma_{2} [(1 - C_{1} )\dot{p}_{2} \dot{\gamma }_{2} + A_{1} (1 + B_{1} )s_{21} ] - A_{1} p_{2} (\gamma_{0}^{\prime \prime 2} - \dot{\gamma }_{2}^{2} )\} , \\ \end{aligned}$$

$$\begin{aligned} & \varphi_{2} = (1 + B_{1} )p_{2} s_{21} + x^{\prime}_{0} (\dot{\gamma }_{2}^{2} + b^{ - 1} \gamma^{\prime\prime}_{0} ) - A_{1}^{ - 1} \dot{p}_{2} [A_{1}^{ - 1} \gamma_{2} \dot{p}_{2} + (1 - C_{1} )p_{2} \dot{\gamma }_{2} ] \\ & \quad - \,\gamma_{2} [s_{11} + z^{\prime}_{0} b^{ - 1} \gamma^{\prime\prime}_{0} + y^{\prime}_{0} \dot{\gamma }_{2} - k(C_{1} \dot{\gamma }_{2}^{2} - B_{1} \gamma_{0}^{\prime \prime 2} )], \\ & f_{3} = - \omega^{2} \left[ {s_{11} \left( {\chi_{1} \gamma_{2} + \frac{{x^{\prime}_{0} }}{1 - a}\gamma^{\prime\prime}_{0} } \right) + 2p_{2} s_{12} } \right] + C_{1} A_{1}^{ - 1} \dot{p}_{2} \\ & \quad \times \,\left[ {\dot{p}_{2} \left( {\chi_{1} \gamma_{2} + \frac{{x^{\prime}_{0} }}{1 - a}\gamma^{\prime\prime}_{0} } \right) + 2p_{2} \dot{p}_{2} (\chi_{2} \dot{\gamma }_{2} - y^{\prime}_{0} a^{ - 1} \gamma^{\prime\prime}_{0} )} \right] \\ & \quad + \,x^{\prime}_{0} [A_{1} b^{ - 1} s_{22} - bA_{1}^{ - 1} \gamma^{\prime\prime}_{0} \dot{p}_{2}^{2} - \gamma^{\prime\prime}_{0} \dot{\gamma }_{2} (y^{\prime}_{0} a^{ - 1} - \chi_{2} \dot{\gamma }_{2} )] \\ & \quad - y^{\prime}_{0} a^{ - 1} \gamma^{\prime\prime}_{0} \left[ {\dot{\gamma }_{2} \left( {\chi_{1} \gamma_{2} + \frac{{x^{\prime}_{0} }}{1 - a}} \right) - bA_{1}^{ - 1} p_{2} \dot{p}_{2} } \right] \\ & \quad + \,z^{\prime}_{0} a^{ - 1} \left[ {\frac{1}{2}b^{ - 1} (2b - 1)s_{11} \gamma_{2} - \gamma^{\prime\prime}_{0} \left( {\chi_{1} \gamma_{2} + \frac{{x^{\prime}_{0} }}{1 - a}} \right) - p_{2} s_{21} } \right] \\ & \quad + \,y^{\prime}_{0} (2 - b)(1 - b)^{ - 1} [\gamma_{2} (y^{\prime}_{0} a^{ - 1} \gamma^{\prime\prime}_{0} - \,\chi_{2} \dot{\gamma }_{2} ) - a\gamma^{\prime\prime}_{0} p_{2} \dot{p}_{2} ] \\ & \quad + \,k\left\{ {(1 - C_{1} )[(y^{\prime}_{0} a^{ - 1} \gamma^{\prime\prime}_{0} - \chi_{2} \dot{\gamma }_{2} )\gamma_{2} \dot{\gamma }_{2} - a\gamma^{\prime\prime}_{0} p_{2} \dot{p}_{2} \dot{\gamma }_{2} + bA_{1}^{ - 1} \gamma^{\prime\prime}_{0} \gamma_{2} \dot{p}_{2}^{2} ]} \right. \\ & \quad + \,p_{2} \gamma^{\prime\prime}_{0} [2b\dot{p}_{2} \dot{\gamma }_{2} - aA_{1} (1 + B_{1} )s_{21} + 2A_{1} s_{21} ] \\ & \left. {\quad + \,A_{1} \left[ {\left( {\chi_{1} \gamma_{2} + \gamma^{\prime\prime}_{0} \frac{{x^{\prime}_{0} }}{1 - a}} \right)(\gamma_{0}^{\prime \prime 2} - \dot{\gamma }_{2}^{2} ) - (1 + B_{1} )\gamma_{2} s_{22} } \right]} \right\}, \\ & \varphi_{3} = - 2s_{12} \gamma_{2} - a\gamma^{\prime\prime}_{0} p_{2} s_{11} + (1 + B_{1} )\left[ {p_{2} s_{22} + \left( {\frac{{x^{\prime}_{0} }}{1 - a}\gamma^{\prime\prime}_{0} + \chi_{1} \gamma_{2} } \right)s_{21} } \right] \\ & \quad + \,2x^{\prime}_{0} \gamma^{\prime\prime}_{0} (b^{ - 1} s_{21} - \,bA_{1}^{ - 1} \dot{p}_{2} \dot{\gamma }_{2} ) - z^{\prime}_{0} b^{ - 1} (a\gamma_{0}^{\prime \prime 2} p_{2} + s_{21} \gamma_{2} ) + A_{1}^{ - 1} (1 - C_{1} ) \\ & \quad \times \,\left[ {bA_{1}^{ - 1} \gamma^{\prime\prime}_{0} p_{2} \dot{p}_{2}^{2} - \left( {\frac{{x^{\prime}_{0} }}{1 - a}\gamma^{\prime\prime}_{0} + \chi_{1} \gamma_{2} } \right)\dot{p}_{2} \dot{\gamma }_{2} + \,b^{ - 1} \gamma_{0}^{\prime \prime 2} x^{\prime}_{0} + p_{2} \dot{\gamma }_{2} (y^{\prime}_{0} a^{ - 1} - \chi_{2} \dot{\gamma }_{2} )} \right] \\ & \quad - \,y^{\prime}_{0} \gamma^{\prime\prime}_{0} (ap_{2} \dot{\gamma }_{2} - bA_{1}^{ - 1} \gamma_{2} \dot{p}_{2} ) + A_{1}^{ - 2} [2\dot{p}_{2} (y^{\prime}_{0} a^{ - 1} \gamma^{\prime\prime}_{0} - \,\chi_{2} \dot{\gamma }_{2} )\gamma_{2} - a\gamma^{\prime\prime}_{0} p_{2} \dot{p}_{2}^{2} ] \\ & \quad + \,k\gamma^{\prime\prime}_{0} [ap_{2} (C_{1} \dot{\gamma }_{2}^{2} - B_{1} \gamma_{0}^{\prime \prime 2} ) - 2\gamma_{2} (bA_{1}^{ - 1} C_{1} \dot{p}_{2} \dot{\gamma }_{2} + B_{1} s_{21} )]. \\ \end{aligned}$$

The last equation of (2) gives the first integral of the system (14) as follows [6]:

$$\gamma_{2}^{2} + \dot{\gamma }_{2}^{2} + 2\mu^{ - 1} \gamma^{\prime\prime}_{0} (ap_{2} \gamma_{2} - bA_{1}^{ - 1} \dot{p}_{2} \dot{\gamma }_{2} + s_{21} ) + \ldots = 1 - \gamma_{0}^{\prime \prime 2} .$$

(15)

In the next, we will look for the periodic solutions of the system (14) under the conditions $A > B > C$ or $A < B < C$ ($\omega^{2}$ is positive). The first condition [7] gives the slow rotation of the body about the major axis of the ellipsoid of inertia and the second gives a slow rotation of the body about the minor axis of the ellipsoid of inertia. We apply the large parameter method [8] to solve the autonomous system (14).

The formal construction of the periodic solutions for a rational value of the natural frequency $\omega$

We achieve the periodic solutions $p_{2} (\tau ,\mu^{ - 1} ),\,\,\dot{p}_{2} (\tau ,\mu^{ - 1} ),\,\,\gamma_{2} (\tau ,\mu^{ - 1} ),\,\,\dot{\gamma }_{2} (\tau ,\mu^{ - 1} )$ of system (14) when:

$$p_{2} (0,0) = \dot{p}_{2} (0,0) = \dot{\gamma }_{2} (0,\mu^{ - 1} ) = 0\,.$$

(16)

The generating system of (14) is obtained when $\mu \to \infty$ in the form:

$$\ddot{p}_{2}^{(0)} + \omega^{2} p_{2}^{(0)} = 0,\quad \ddot{\gamma }_{2}^{(0)} + \gamma_{2}^{(0)} = 0.$$

(17)

So the periodic solutions of system (17), when the period $T_{0} = 2\pi n$, become:

$$p_{2}^{(0)} = M_{1} \cos \,\omega \tau + M_{2} \sin \,\omega \tau ,\quad \gamma_{2}^{(0)} = M_{3} \cos \,\tau ,$$

(18)

where $M_{i} ,\,\,\,i = 1,2,3$ are constants to be determined.

Assuming the following solutions of system (14) [9]:

$$\begin{aligned} & p_{2} (\tau ,\mu ) = \tilde{M}_{1} \cos \,\omega \tau + \tilde{M}_{2} \sin \,\omega \tau + \sum\limits_{k = 2}^{\infty } {\mu^{ - k} G_{k} (\tau )} , \\ & \gamma_{2} (\tau ,\mu ) = \tilde{M}_{3} \cos \,\tau + \sum\limits_{k = 2}^{\infty } {\mu^{ - k} H_{k} (\tau )} \,, \\ \end{aligned}$$

(19)

with a period $T(\,\mu^{ - 1} ) = T_{0} + \alpha (\,\mu^{ - 1} )$ which reduces to (18) at $\mu \to \infty$. Let us define the quantities $\tilde{M}_{i} ,\,\,i = 1,2,3$ as follow:

$$\tilde{M}_{i} = M_{i} + \beta_{i} (\mu^{ - 1} ),\quad i = 1,2,3,$$

(20)

where $\beta_{i}$ are functions of $\mu^{ - 1}$ which represent the deviations of the initial values of $p_{2} ,\,\,\dot{p}_{2} ,\,\,\gamma_{2}$ for system (14) from their initial values of generating system (17) such that $\beta_{i} (\,0) = 0$.

Let us express the initial conditions (16) by the relations:

$$p_{2} (0,\mu^{ - 1} ) = \tilde{M}_{1} ,\quad \dot{p}_{2} (0,\mu^{ - 1} ) = \omega \tilde{M}_{2} ,\quad \gamma_{2} (0,\mu^{ - 1} ) = \tilde{M}_{3} ,\quad \dot{\gamma }_{2} (0,\mu^{ - 1} ) = 0.$$

(21)

We rewrite the periodic solutions (18) in the form:

$$p_{2}^{(0)} = E\cos \,(\omega \tau - \varepsilon ),\quad \gamma_{2}^{(0)} = M_{3} \cos \,\tau ,$$

(22)

where $E = \sqrt {M_{1}^{2} + M_{2}^{2} }$ and $\varepsilon = \tan^{ - 1} {{M_{2} } \mathord{\left/ {\vphantom {{M_{2} } {M_{1} }}} \right. \kern-0pt} {M_{1} }}$. Using (22) and (12), we get:

$$\begin{aligned} & s_{11}^{(0)} = E^{2} [(a\cos^{2} \varepsilon - 0.5\,) + b\omega^{2} A_{1}^{ - 2} (\sin^{2} \varepsilon - 0.5) + 0.5(b\omega^{2} A_{1}^{ - 2} - a)\cos 2(\omega \tau - \varepsilon )] \\ & \quad - \,2M_{3} [x^{\prime}_{0} (1 - \cos \tau ) + y^{\prime}_{0} \sin \tau ] - 0.5kM_{3}^{2} C_{1} (1 - \cos 2\tau ), \\ & s_{12}^{(0)} = \text{a}^{\text{2}} E\gamma^{\prime\prime}_{0} x^{\prime}_{0} (1 - a)^{ - 1} [\cos \varepsilon - \cos (\omega \tau - \varepsilon )] \\ & \quad + \,Ey^{\prime}_{0} b\gamma^{\prime\prime}_{0} A_{1}^{ - 1} \omega (1 - a^{ - 1} A_{1}^{ - 1} )[\sin \varepsilon + \sin (\omega \tau - \varepsilon )] \\ & \quad + \,aEM_{3} (\chi_{1} + a\gamma^{\prime\prime}_{0} k)\{ \cos \varepsilon - 0.5\cos [(\omega + 1)\tau - \varepsilon ] - 0.5\cos [(\omega - 1)\tau - \varepsilon ]\} \\ & \quad + \,\omega bA_{1}^{ - 1} EM_{3} (\chi_{2} A_{1}^{ - 1} - bk\gamma^{\prime\prime}_{0} )\{ 0.5\cos [(\omega + 1)\tau - \varepsilon ] - 0.5\cos [(\omega - 1)\tau - \varepsilon ]\} \\ & \quad + \,EM_{3} (z^{\prime}_{0} - \gamma^{\prime\prime}_{0} k)\{ a\cos \varepsilon + 0.5(b\omega A_{1}^{ - 1} - a)\cos [(\omega - 1)\tau - \varepsilon ] \\ & \quad - \,0.5(b\omega A_{1}^{ - 1} + a)\cos [(\omega + 1)\tau - \varepsilon ]\} , \\ & s_{21}^{(0)} = EM_{3} \{ a\cos \varepsilon + 0.5(b\omega A_{1}^{ - 1} - a)\cos [(\omega - 1)\tau - \varepsilon ] \\ & \quad - \,0.5(b\omega A_{1}^{ - 1} + a)\cos [(\omega + 1)\tau - \varepsilon ]\} , \\ & s_{22}^{(0)} = E^{2} \gamma^{\prime\prime}_{0} [a^{2} (\cos^{2} \varepsilon - 0.5) + b^{2} \omega^{2} A_{1}^{ - 2} (\sin^{2} \varepsilon - 0.5) \\ & \quad - \,0.5(a^{2} - b^{2} \omega^{2} A_{1}^{ - 2} )\cos 2(\omega \tau - \varepsilon )] + 0.5M_{3}^{2} (a\chi_{1} + bA_{1}^{ - 1} \chi_{2} )(1 - \cos 2\tau ) \\ & \quad + \,\gamma^{\prime\prime}_{0} M_{3} [a\,x^{\prime}_{0} (1 - a)^{ - 1} (1 - \cos \tau ) + by^{\prime}_{0} \,a^{ - 1} A_{1}^{ - 1} \sin \tau ]. \\ \end{aligned}$$

(23)

Substituting (22) and (23) into (16), we get:

$$\begin{aligned} & F_{2}^{(0)} = M_{1} L(\omega )\cos \omega \tau + M_{2} L(\omega )\sin \omega \tau + \ldots , \\ & \varPhi_{2}^{(0)} = M_{3} N(\omega )\cos \tau + \ldots , \\ & F_{3}^{(0)} = M_{1} K(\omega )\cos \omega \tau + M_{2} K(\omega )\sin \omega \tau + \ldots , \\ \end{aligned}$$

(24)

where:

$$\begin{aligned} & L(\omega ) = \omega^{2} [ - (aM_{1}^{2} + b\omega^{2} A_{1}^{ - 2} M_{2}^{2} ) + 0.25(M_{1}^{2} + M_{2}^{2} )(C_{1} A_{1}^{ - 1} + 3a + b\omega^{2} A_{1}^{ - 2} )] \\ & \quad + \,2\omega^{2} M_{3} x^{\prime}_{0} - \gamma^{\prime\prime}_{0} [z^{\prime}_{0} a^{ - 1} + a\chi_{1} (1 - \omega^{2} )] + k\{ A_{1} (\gamma_{0}^{\prime \prime 2} - 0.5M_{3}^{2} ) \\ & \quad + \,0.5M_{3}^{2} [a(A_{1} - 2\omega^{2} ) + \omega^{2} b]\} , \\ & N(\omega ) = - (aM_{1}^{2} + b\omega^{2} A_{1}^{ - 2} M_{2}^{2} ) - 0.5(M_{1}^{2} + M_{2}^{2} )[aB_{1} + \omega^{2} A_{1}^{ - 2} (1 - b)] \\ & \quad + \,2M_{3} x^{\prime}_{0} - \gamma^{\prime\prime}_{0} \,[z^{\prime}_{0} b^{ - 1} - a\chi_{1} (1 - \omega^{2} )] + k[M_{3}^{2} (b - a) - B_{1} \gamma_{0}^{\prime \prime 2} ], \\ & K(\omega ) = - 2\omega^{2} \gamma^{\prime\prime}_{0} \,[a^{2} x^{\prime}_{0} M_{1} (1 - a)^{ - 1} + \omega y^{\prime}_{0} bM_{2} (1 - a^{ - 1} A_{1}^{ - 1} )A_{1}^{ - 1} \\ & \quad + \,0.25a^{ - 1} A_{1}^{ - 1} C_{1} y^{\prime}_{0} (M_{1}^{2} + M_{2}^{2} )] - 2\omega^{2} aM_{1} M_{3} [\chi_{1} + a\gamma^{\prime\prime}_{0} k + z^{\prime}_{0} - \gamma^{\prime\prime}_{0} k] \\ & \quad + \,aM_{1} M_{3} [2A_{1} k\gamma^{\prime\prime}_{0} - z^{\prime}_{0} a^{ - 1} - aA_{1} k\gamma^{\prime\prime}_{0} (1 + B_{1} ) - \chi_{1} (1 + B_{1} )]. \\ \end{aligned}$$

(25)

Using (24) and (25), the following functions are obtained:

$$\begin{array}{*{20}l} {g_{2} (T_{0} ) = - \pi n\omega^{ - 1} M_{2} L(\omega ),} \hfill & {\dot{g}_{2} (T_{0} ) = \pi nM_{1} L(\omega ),} \hfill \\ {h_{2} (T_{0} ) = 0,} \hfill & {\dot{h}_{2} (T_{0} ) = \pi nM_{3} N(\omega ),} \hfill \\ {g_{3} (T_{0} ) = - \pi n\omega^{ - 1} M_{2} K(\omega ),} \hfill & {\dot{g}_{3} (T_{0} ) = \pi nM_{1} K(\omega ).} \hfill \\ \end{array}$$

(26)

Substituting by the initial conditions (21) into the first integration (17) when $\tau = 0$, we get:

$$M_{3}^{2} + 2M_{3} \beta_{3} + \beta_{3}^{2} + 2\mu^{ - 1} a\gamma^{\prime\prime}_{0} M_{3} (M_{1} + \beta_{1} ) = 1 - \gamma_{0}^{\prime \prime 2} .$$

(27)

Let γ ^″₀ depends on μ⁻¹, we get:

$$\gamma^{\prime\prime}_{0} = \mu^{ - 1} \varGamma ,\quad 0 < \varGamma < 1.$$

(28)

Taking into consideration, Eqs. (27) and (28), we get $M_{3} ,\,\,\beta_{3}$ as follows:

$$M_{3} = 1,\quad \beta_{3} = - a\varGamma \mu^{ - 2} \tilde{M}_{1} - \frac{1}{2}\mu^{ - 2} \varGamma^{2} + \ldots .$$

(29)

The independent conditions for periodicity are:

$$\begin{aligned} & - (L_{1} (\omega ) - \omega^{2} N_{1} (\omega ))\pi n\omega^{ - 1} \tilde{M}_{2} + \mu^{ - 1} \,G_{3} (T_{0} ) + \ldots = 0, \\ & (L_{1} (\omega ) - \omega^{2} N_{1} (\omega ))\pi n\tilde{M}_{1} + \mu^{ - 1} \,\dot{G}_{3} (T_{0} ) + \ldots = 0, \\ & \mu^{ - 2} (\dot{H}_{2} (T_{0} ) + \mu^{ - 1} \dot{H}_{3} (T_{0} ))\tilde{M}_{3}^{ - 1} + \ldots = \alpha (\mu^{ - 1} ), \\ \end{aligned}$$

(30)

where $L_{1} (\omega ),\,\,N_{1} (\omega )$ are obtained from $L(\omega ),\,\,N(\omega )$ replacing $M_{i}$ by $(M_{i} + \beta_{i} ),\,\,i = 1,2,3$ to get:

$$L_{1} (\omega ) - \omega^{2} N_{1} (\omega ) = W_{0} (\omega )(\tilde{M}_{1}^{2} + \tilde{M}_{2}^{2} ) - \gamma^{\prime\prime}_{0} [z^{\prime}_{0} \,W_{1} (\omega ) + k\gamma^{\prime\prime}_{0} \,W_{2} (\omega )] - kW_{3} (\omega )\tilde{M}_{3}^{2} ,$$

(31)

where:

$$\begin{array}{*{20}l} {W_{0} (\omega ) = (a - 1)(a + b - 2)/2b,} \hfill & {W_{1} (\omega ) = [3(a + b) - 2(2ab + 1)]/ab,} \hfill \\ {W_{2} (\omega ) = 2\omega^{2} [1 - (a + b)],} \hfill & {W_{3} (\omega ) = \omega^{2} b.} \hfill \\ \end{array}$$

(32)

For zeros approximation for power series of $1/\mu$, Eq. (30) give:

$$M_{1} = M_{2} = 0\,.$$

(33)

Since the z-axis is directed along with the major or the minor axis of the ellipsoid of inertia of the body, we get: $W_{0} (\omega ) > 0$ for all $\omega$ under consideration.

Assume that:

$$\gamma^{\prime\prime}_{0} \,[z^{\prime}_{0} \,W_{1} + k\gamma^{\prime\prime}_{0} \,W_{2} ] + kW_{3} (\omega )M_{3}^{2} \ne 0\,.$$

(34)

Using (30), we get $\beta_{1} ,\,\,\beta_{2}$ in power series expansions of powers less than $\mu^{ - 2}$. Then for the rational values of the natural frequency $\omega$ does not equal to $(1,2,1/2,3,1/3, \ldots )$, we get the required periodic solutions and the correction of the period $\alpha (\mu^{ - 1} )$ as:

$$\begin{aligned} & p_{1} (\tau ,\mu^{ - 1} ) = \mu^{ - 1} [x^{\prime}_{0} (1 - a)^{ - 1} \gamma^{\prime\prime}_{0} + \chi_{1} M_{3} \cos \tau ] + \ldots , \\ & q_{1} (\tau ,\mu^{ - 1} ) = \mu^{ - 1} [y^{\prime}_{0} (1 - b)^{ - 1} \gamma^{\prime\prime}_{0} + A_{1}^{ - 1} M_{3} \chi_{2} \sin \tau ] + \ldots , \\ & r_{1} (\tau ,\mu^{ - 1} ) = 1 - 0.25\mu^{ - 2} M_{3} [kM_{3} C_{1} + 4x^{\prime}_{0} (1 - \cos \tau ) + y^{\prime}_{0} \sin \tau - kM_{3} C_{1} \cos 2\tau ] + \ldots , \\ & \gamma (\tau ,\mu^{ - 1} ) = M_{3} \cos \tau - 0.5\mu^{ - 2} \varGamma^{2} \cos \tau + \ldots , \\ & \gamma^{\prime}(\tau ,\mu^{ - 1} ) = - M_{3} \sin \tau + 0.5\mu^{ - 2} \varGamma^{2} \sin \tau + \ldots , \\ & \gamma^{\prime\prime}(\tau ,\mu^{ - 1} ) = \gamma^{\prime\prime}_{0} + \mu^{ - 2} M_{3} [x^{\prime}_{0} (1 - a)^{ - 1} \gamma^{\prime\prime}_{0} - 0.5M_{3} C_{1} \left( {\frac{{z^{\prime}_{0} }}{a + b - 1} + 0.5k\gamma^{\prime\prime}_{0} } \right) \\ & \quad - x^{\prime}_{0} (1 - a)^{ - 1} \gamma^{\prime\prime}_{0} \cos \tau + y^{\prime}_{0} (1 - b)^{ - 1} \gamma^{\prime\prime}_{0} \sin \tau + 0.25M_{3} C_{1} (\frac{{2z^{\prime}_{0} }}{a + b - 1} + k\gamma^{\prime\prime}_{0} )\cos 2\tau ] + \ldots , \\ \end{aligned}$$

(35)

$$\alpha (\mu^{ - 1} ) = 2\mu^{ - 2} \pi n\left\{ {M_{3} x^{\prime}_{0} - z^{\prime}_{0} \gamma^{\prime\prime}_{0} - 0.5k} \left[ {\gamma_{0}^{\prime \prime 2} (bB_{1} - aA_{1} ) + \,B_{1} \gamma^{\prime\prime}_{0} (1 - \gamma^{\prime\prime}_{0} )} - 0.125C_{1} M_{3}^{2} \right] \right\} + \ldots \,.$$

(36)

The obtained solutions (35) and (36) are considered as the generalization of the corresponding problem in gravity field which studied in previous works [10] (when k = 0), the deviations between them are given by:

$$\begin{aligned} & \Delta p_{1} = \mu^{ - 1} (1 - \omega^{2} )^{ - 1} kM_{3} \gamma^{\prime\prime}_{0} (A_{1} - \omega^{2} )\cos \tau + \ldots , \\ & \Delta q_{1} = \mu^{ - 1} A_{1}^{ - 1} M_{3} \left\{ {\left. {(1 - \omega^{2} )^{ - 1} [k(A_{1} - \omega^{2} )\gamma^{\prime\prime}_{0} ] - kA_{1} \gamma^{\prime\prime}_{0} } \right\}} \right.\sin \tau + \ldots , \\ & \Delta r_{1} = - 0.25\mu^{ - 2} M_{3} kC_{1} (1 - \cos 2\tau ) + \ldots , \\ & \Delta \gamma = \mu^{ - 2} [0] + \ldots ,\quad \Delta \gamma^{\prime} = \mu^{ - 2} [0] + \ldots , \\ & \Delta \gamma^{\prime\prime} = - 0.25\mu^{ - 2} M_{3}^{2} C_{1} k\gamma^{\prime\prime}_{0} (1 - \cos 2\tau ) + \ldots , \\ \end{aligned}$$

$$\Delta \alpha (\mu^{ - 1} ) = - \mu^{ - 2} \pi nk\left[ {\gamma_{0}^{\prime \prime 2} (bB_{1} - aA_{1} ) - 0.5M_{3}^{2} (b - a) + B_{1} \gamma^{\prime\prime}_{0} (1 - \gamma^{\prime\prime}_{0} )} \right] + \ldots \,\,.$$

(37)

Geometric interpretation of the motion

The geometric interpretation for the motion of the body at any instant of time to Euler’s angles definitions $\theta ,\,\,\psi ,\,\,\phi$ is given by [11]:

$$\begin{aligned} & \theta = \theta_{0} + \mu^{ - 2} {\text{cosec}}\,\theta_{0} [\theta_{2} (t + t_{0} ) - \theta_{2} (t_{0} )] + \ldots , \\ & \psi = \psi_{0} + 0.5Mg\ell\,C^{ - 1} r_{0}^{ - 1} (\chi_{1} - \chi_{2} A_{1}^{ - 1} )t + \mu^{ - 1} \sqrt {Mg\ell /C} \,[\psi_{1} (t + t_{0} ) - \psi_{1} (t_{0} )] + \ldots , \\ & \phi = \phi_{0} + [r_{0}^{ - 1} - 0.5Mg\ell C^{ - 1} r_{0}^{ - 1} \cos \theta_{0} (\chi_{1} - \chi_{2} A_{1}^{ - 1} )]t + \ldots , \\ \end{aligned}$$

(38)

where:

$$\begin{aligned} \theta_{2} (t) &= M_{3} \left[ {y^{\prime}_{0} \,\gamma^{\prime\prime}_{0} (1 - b)^{ - 1} \sin r_{0}^{ - 1} t - x^{\prime}_{0} \gamma^{\prime\prime}_{0} (1 - a)^{ - 1} \cos r_{0}^{ - 1} t} \right. \\ & \left. {\quad+ \,0.5C_{1} \left( {\left. {\frac{{z^{\prime}_{0} }}{a + b - 1} - 0.667k\gamma^{\prime\prime}_{0} } \right)\cos 2r_{0}^{ - 1} t} \right.} \right], \\ & {\psi_{1}(t) = 0.25 (\chi_{1} + \chi_{2}{A_{1}^{-1}}) + \sin 2 r_{0}^{ - 1}t + y_{0}^{\prime} {\gamma_{0}^{\prime\prime}}(1-b)^{-1} \cos r_{0}^{-1}t + x_{0}^{\prime}(1-a)^{-1}\gamma_{0}^{\prime\prime}\sin r_{0}^{-1}t.}\\ \end{aligned}$$

(39)

Numerical solutions

In this section, we use a computer program to determine the obtained solutions (19) and their derivatives for the time in the interval t $\in$ [0, 300]. On the other hand, we use the fourth-order Runge–Kutta method [12] through another program to obtain numerical solutions for the autonomous system (14). In the end, we compare both solutions to check the accuracy of the method of solutions. These results are obtained through Tables 1 and 2. From these Tables, we deduce that the numerical solutions are in agreement with the analytical ones which prove the accuracy of considered methods.

Table 1 Represents the values of the obtained analytical solutions using the large parameter method in the interval t $\in$ [0, 300]

Table 2 represents the values of the numerical solutions using fourth-order Runge–Kutta method in the interval t $\in$ [0, 300]

Conclusions

We conclude that the problem of the motion of a rigid body about a fixed point is studied in many works [13,14,15,16,17,18] in both the uniform and gravity fields. We study our problem in case of a right angle of nutation $\theta_{0}$ when its center of mass does not necessarily coincide with the fixed point. The equations of motions of the problem are obtained and reduced to a quasilinear autonomous system. The obtained system is solved by assuming a large parameter achieved from an angular velocity component tends to zero. The obtained solutions are treated through computer programs in a bounded interval of time. The autonomous system is treated with the Runge–Kutta method in the same interval of time to obtain the numerical solutions of the motion. Both obtained solutions are in full agreement with others which prove the accuracy of both numerical and analytical techniques used in solving the problem. For the geometric interpretation obtained, we note that:

The precession angle $\psi$

From (38) when $\mu \to \infty$, we deduce that the precession angle $\psi$ is sufficiently large because r₀ is sufficiently small, that is, we obtain a case of large precession $\psi = \psi_{0} + 0.5Mg\ell\,C^{ - 1} r_{0}^{ - 1} (\chi_{1} - \chi_{2} A_{1}^{ - 1} )\,t$.

The nutation angle $\theta$

We obtain a case of steady regular permutation: $\theta = \theta_{0}$.

The pure rotation angle $\phi$

The case of a large pure rotation is obtained which depends on 1/r₀ in the form:

$$\phi = \phi_{0} + [r_{0} - 0.5Mg\ell C^{ - 1} r_{0}^{ - 1} \cos \theta_{0} (\chi_{1} - \chi_{2} A_{1}^{ - 1} )]t.$$

The large parameter technique used here is considered as the only one suitable for this problem in the origin domain of r_o tends to zero. Poincaré–Lindstedt method or Krylov Boboliubov Mitropolski one is failed to solve this problem because they depend on achieving a small parameter in domain r₀ tends to infinity. We conduct a comparison of the results of this manuscript with the results of the previous work. The results were obtained in [19] deals with the disk problem which satisfies the symmetry moments of inertia about two principal axes of the ellipsoid of inertia but our results here treat the general rigid problem in a limiting value of the Euler’s angle $\theta_{o} \approx {\raise0.7ex\hbox{$\pi $} \!\mathord{\left/ {\vphantom {\pi 2}}\right.\kern-0pt} \!\lower0.7ex\hbox{$2$}}$. The advantage of our used technique [20] depends on a large parameter $\mu \to \infty$. The obtained solutions are checked using two programs to assert their accuracy through Tables 1 and 2. The main results in our work are the obtained analytical solutions in Eq. (35) which is represented through computerized digital data in Table 1. The secondary results are proving the validity of these solutions which are given through the Runge–Kutta method in Table 2. So we confirm that the presented numerical results would correlate to results obtained by other distinguished numerical techniques (except the Runge–Kutta method) will reveal the same results obtained in Table 1.