Abstract
In this paper, we consider the problem of the rotational motion of a rigid body with an irrational value of the frequency . The equations of motion are derived and reduced to a quasilinear autonomous system. Such system is reduced to a generating one. We assume a large parameter proportional inversely with a sufficiently small component of the angular velocity which is assumed around the major or the minor axis of the ellipsoid of inertia. Then, the large parameter technique is used to construct the periodic solutions for such cases. The geometric interpretation of the motion is obtained to describe the orientation of the body in terms of Euler’s angles. Using the digital fourth-order Runge-Kutta method, we determine the digital solutions of the obtained system. The phase diagram procedure is applied to study the stability of the attained solutions. A comparison between the considered numerical and analytical solutions is introduced to show the validity of the presented techniques and solutions.
1. Introduction
The rigid body problem of mass that rotates about a fixed point is classified according to the natural frequency value in either the uniform gravity field of acceleration or a Newtonian force one. The case of the rational value of the natural frequency is studied in [1] for a rotating heavy solid about a fixed point with a small velocity about one of the axes of the ellipsoid of inertia. The solutions obtained contain singular cases of the natural frequencies named when . The singular cases that appeared for the natural frequency values like (the disk case) and are considered in [2, 3], respectively. It remains for us to study four other cases until the solution of the problem is completed up to the third approximation for any natural frequency of movement. Such cases are classified according to the natural frequency values named; the state of irrational values of the natural frequencies which is the subject of this article besides three singular cases will be studied in the future in shaa Allah classified when .
Let the framebe fixed in the body and the framebe fixed in space. Assume that , , and represent the moments of inertia of the body in the moving coordinate system. Suppose () is the mass center of the body. In the case when is irrational, the equations of motion are obtained and reduced to the following system:wherewhere symbols like (ab) mean cyclic permutations and indicate equations which are omitted; and are the -components of the angular velocity vector and the fixed unit vector in space, respectively; and is a constant that depends on the rigid body parameters. The system (1) has the following first integral:where , , and are constants that depend on the rigid body parameters, and the function satisfies the condition when .
2. The Periodic Solutions
In this section, we obtain the generating system of (1), and then, we will solve it under a new condition of motion. We adapt a large parameter method [1] to solve the system (1) in the presence of the first integral (4). By putting into the system (1), we get the generating system as in [4] which has generating periodic solutions of a period .
In this case, from the system (1), the required periodic solutions with a period are assumed in the following forms [5]:where is an arbitrary constant and and are analytic functions of which vanish when . The function is an analytic function of which vanishes when . The functions and are obtained by substituting (5) into (1) and equating coefficients of like powers of
The functions can be written in an expansion series as follows [6]:
Using Equation (3), the periodicity conditions are obtained as follows [7]:
Under condition (7) and series (6), it produces an endless system of equations that designates the coefficients . The functions and are expressed as in [8].
Now, we aim to find the periodic solutions of the generalized system. According to [9], series (6) begins from a term of order not lower than . It follows that the expansion of represents a power series that begins from a term of order not lower than . In this case, the quantities and are obtained.
Introduce the following variables [10]:where and are constants that depend on the rigid body parameters. The vectors and represent the components of the angular velocity and the unit vector of the downward fixed axis in space, respectively, and , satisfies the condition when.
Let be sufficiently small; we define a large parameter and apply the large parameter method [3] to obtain the analytical solutions in the form of power series expansions in terms of as follows:where are constants that depend on the rigid body parameters, and
The correction of the period is obtained in the form
3. The Geometric Interpretation of Motion
In this section, we discuss the problem geometrically to show the orientation of the body at any instant in time. Substituting Equation (9) into Euler’s angles , and in which has been replaced by , using , we obtain the following expressions for the angles [11]:wherewhere is a constant depending on the moments of inertia.
4. Numerical Considerations
In this section, we study the analytical and numerical solutions for our problem as follows.
4.1. The Analytical Solutions
We rewrite the analytical solutions in the following form:
Let the step and , where is the maximum value of the variable . Let us introduce the following data:
Let the body parameters be
Introduce the following computerized symbols:
Using the above data and a computer program, we obtain the analytical solutions in the graphs (see Figures 1 and 2).
4.2. The Numerical Solutions
Using (18), we rewrite the system (1) in the form
Using the fourth-order Runge-Kutta method [12] through a computer program and the data (16) and (17) with the same initial values of the analytical solutions, we obtain the numerical solutions in the graphs (see Figures 3 and 4). For checking the accuracy of both solutions, we draw the graphs (see Figures 5 and 6). We find agreement between the analytical and numerical solutions which satisfied the excellent results for the analytical and numerical techniques. The smooth simple curves obtained (see Figures 7–12) show that the obtained solutions are stable [13, 14].
5. Conclusions
In this section, we conclude that the problem of the motion of a rigid body in a uniform gravity field with an irrational value of the natural frequency which is excluded from the previous works [1–3] is considered. The equations of motion and their first integrals in the presence of new conditions of motion are obtained and reduced to a semilinear autonomous system of two degrees of freedom and one first integral. We assumed a large parameter that is achieved to be inversely proportional to the angular velocity component which is supposed to be sufficiently small. Under this assumption, the well-known Poincare method [15] cannot solve this problem because we cannot achieve the small parameter which must be proportional to a sufficiently high angular velocity component . Therefore, we solve the problem with the large parameter technique. The advantages of this method are as follows: using small energy at the initial moment instead of high energy, obtaining a slow gyroscopic motion instead of fast gyroscopic motion, and giving analytical and numerical solutions in a new domain of the motion . Also, the correction of the period for these solutions is obtained in the new domain. The geometric interpretation of motion is attained. When , we obtain from Equation (17) the cases of regular permutation and pure rotation. Applying the numerical fourth-order Runge-Kutta method [12] through a computer program, we find the numerical solutions for this problem. Another computer program is carried out on the analytical obtained solutions using the large parameter technique and assuming the same initial values and data (16) and (17). We get the numerical results for both solutions and their graphical representations. The agreement of results through the graphical representations shows the advantage of both techniques for obtaining the high accuracy of the obtained solutions. This problem has many applications in aerospace sciences and technologies because of the wide use of gyros in these fields [16]. The procedures used here are useful for solving sophisticated problems such as [17] in a new domain of the considered parameters. This can be done by reflecting the problem parameters. In shaa Allah in next papers, we study the remaining singular cases of the natural frequency values to complete the solutions for the problem up to the third approximation.
Data Availability
Data sharing is not applicable to this article as no datasets were generated or analyzed during the current study.
Conflicts of Interest
The author declares that he has no competing interests.