Research Article | Open Access

Volume 2020 |Article ID 8854136 | https://doi.org/10.1155/2020/8854136

A. I. Ismail, "Applying the Large Parameter Technique for Solving a Slow Rotary Motion of a Disc about a Fixed Point", International Journal of Aerospace Engineering, vol. 2020, Article ID 8854136, 7 pages, 2020. https://doi.org/10.1155/2020/8854136

# Applying the Large Parameter Technique for Solving a Slow Rotary Motion of a Disc about a Fixed Point

Revised20 Aug 2020
Accepted27 Aug 2020
Published16 Sep 2020

#### Abstract

In this paper, the motion of a disk about a fixed point under the influence of a Newtonian force field and gravity one is considered. We modify the large parameter technique which is achieved by giving the body a sufficiently small angular velocity component about the fixed -axis of the disk. The periodic solutions of motion are obtained in the neighborhood tends to 0. This case of study is excluded from the previous works because of the appearance of a singular point in the denominator of the obtained solutions. Euler-Poison equations of motion are obtained with their first integrals. These equations are reduced to a quasilinear autonomous system of two degrees of freedom and one first integral. The periodic solutions for this system are obtained under the new initial conditions. Computerizing the obtained periodic solutions through a numerical technique for validation of results is done. Two types of analytical and numerical solutions in the new domain of the angular velocity are obtained. Geometric interpretations of motion are presented to show the orientation of the body at any instant of time .

#### 1. Introduction

In , the authors considered the limiting case for the motion of a rigid body about a fixed point in the Newtonian force field and gravity one. The small parameter technique is applied to solve this problem. The authors defined this parameter inversely proportional to the third component of the angular velocity which is assumed to be sufficiently high. The periodic solutions and their graphical representations are obtained and illustrated geometrically. In , the authors admitted the KBM technique for solving the problem of a rotating heavy solid about a fixed point under the influence of a gyrostatic moment. They assumed a small parameter as in  and found the analytical and numerical solutions for the body which moves under its gravity and a gyro moment about the minor or the major axis of the ellipsoid of inertia. In , Leshchenko and Ershkov presented a new type of solving procedure for Euler-Poisson equations (rigid body rotation over a fixed point) in the presence of some restricted conditions on the body angular velocity or the applied perturbing torques. The author in  gave the regular precession of an asymmetric rigid body acted upon by a uniform gravity field and magnetic one. He obtained the equations of motion of the body and reduced them to a quasilinear autonomous system. He found the solution to the problem and its geometric interpretation of motion. Nayfeh  presented some perturbation methods such as Poincare’, the KBM, the Multiple scales, and averaging techniques for solving this kind of motion under certain conditions. The authors in  studied the rotating symmetric rigid body about a fixed point in the Newtonian force field in a case analogous to Kovalevskaya’s problem. They described the motion of the body and derived its equations of motion and find the solution to the problem assuming Kovalevskaya’s conditions.

All the previous approximated methods depend on a small parameter achieved inversely proportional to the sufficiently high angular velocity component. This study gives many applications in physics , gyros [8, 9], astronomy, engineering, aerospace, and other sciences.

In our problem, we assume a slow rotation (weak spin) of the body instead of a fast rotation (high spin). We introduce a large parameter inversely proportional to the weak spin about the -axis instead of a small parameter inversely proportional to the high spin about that axis. So, we must apply the large parameter technique instead of the small one used in the previous works. The validation of our results will be given.

#### 2. Formulation of the Problem and Construction of the Periodic Solutions

In this section, we formulate the problem of the motion of the body, deduce the equations of motion, and construct the periodic solutions:

##### 2.1. Formulation of the Problem

Consider a rigid body of mass (M) rotates around a fixed point . Let the two coordinate frames and are fixed in the space and the body, respectively. Initially, let the body spin with a small speed angular velocity component about the -axis. Suppose that and are the generalized coordinates of the reduced quasilinear autonomous system . Let be the gravity acceleration; are the principal moments of inertia of the body for the moving coordinates; the point () is the mass center of the body, and is the position vector of to . Let and be the basic amplitudes for while is the basic amplitude for . Considering , and are the deviations for and from their initial basic amplitudes , , and respectively.

Let the moving -axis make an angle with the downward fixed -axis. This case of a slow spinning rotary body when the natural frequency gives the rotary motion of a disc which is excluded from the previous case . We achieve a large parameter inversely proportional to the sufficiently small value of . Assume that the disc rotates about the -axis in the presence of a Newtonian force field  and a gravitational one. We reduce the equations of motion for this case to the following autonomous system of two degrees of freedom : where where where the symbols like (), (), and () denote omitted equations and , , , and are the components of the angular velocity vector and the cosines direction of the unit vector .

The system (2) has the following integral : where where the quantities and are the initial values of the corresponding variables.

##### 2.2. Construction of the Periodic Solutions

In this subsection, we use the large parameter technique  to construct the periodic solutions of system (2) when the disc rotates with slow velocity about the minor axis of the ellipsoid of inertia. Assume that

Using (8) and the definition of , we get where and and are constants.

Substituting (9) into (3), we get where

Let

Using (15), (11), and (14), we get

Using (16) and McLaren’s expansion in the form : when , we get where

Using (18), the independent periodicity conditions give  where

The expansions are put in power series form as follows:

Substituting (22) into (19) and equating the coefficients of in both sides, we get

From  we get

The following quantities are derived:

Making use of (18), (23), (8), (24), (3), and (25), we obtain the following periodic solutions forms:

The correction of the period is:

#### 3. Geometric Interpretation of the Motion

In this section, we search the description of the motion of the disk at any instant of the time using Euler’s angles [11, 12]: where

These formulas depend on four arbitrary constants and , where is sufficiently small. The obtained analytical periodic solutions are considered a generalization of the ones obtained in [13, 14].

#### 4. Numerical Solutions

In this section, we use a smooth and suitable numerical method for obtaining the approximated numerical solutions for the autonomous system (2). Such a method is named the fourth-order Runge-Kutta method  which is used through a computerized program to find the numerical solutions of the considered problem. In the other side, we computerize the obtained analytical solutions and compare them with the numerical ones aiming to find the errors between them:

##### 4.1. The Analytical Solutions

Rewriting the resulted analytical solutions in the form:

Let the disk parameters are

We obtain the following parameters of the motion:

Making use of (30), (31), and (32) through a computer program, we obtain the values of the analytical solutions (see Table 1).

 0 7.048539 1 -2.028946 0 10 4.52171 0.8102261 -5.775176 -0.5861173 20 2.79E-01 0.3129329 -7.329452 -0.9497752 30 -4.070127 -0.3031333 -6.101851 -0.9529482 40 -6.874124 -0.804146 -2.558307 -0.5944318 50 -7.069064 -0.999947 1.956236 -0.01030037 60 -4.580957 -0.8162205 5.728294 0.5777406 70 -3.54E-01 -0.3226997 7.326192 0.9465014 80 4.007059 0.2933016 6.143451 0.95602 90 6.847407 0.7979805 2.628977 0.6026833 100 7.088839 0.9997878 -1.883318 0.02059964 110 4.639718 0.8221281 -5.680805 -0.5693026 120 4.30E-01 0.3324316 -7.322155 -0.9431273 130 -3.943567 -0.2834384 -6.1844 -0.9589905 140 -6.819966 -0.7917299 -2.699373 -0.6108714 150 -7.107862 -0.9995226 1.810199 -0.03089697 160 -4.697986 -0.8279484 5.632714 0.5608044 170 -5.05E-01 -0.3421288 7.317342 0.9396531 180 3.879655 0.2735457 6.224691 0.961859 190 6.791799 0.7853958 2.769477 0.6189939 200 7.12613 0.9991513 -1.736891 0.04119054 210 4.755757 0.8336813 -5.584021 -0.5522458 220 5.80E-01 0.3517892 -7.311752 -0.9360793 230 -3.815332 -0.263623 -6.264326 -0.9646258 240 -6.762911 -0.7789784 -2.839286 -0.6270508 250 -7.143644 -0.998674 1.663392 -0.05148069 260 -4.813022 -0.8393253 5.534739 0.5436295 270 -6.56E-01 -0.3614123 7.305387 0.9324061 280 3.750605 0.2536723 6.303295 0.9672902 290 6.733307 0.7724778 2.908801 0.6350418 300 7.160398 0.9980907 -1.589723 0.06176443
##### 4.2. The Numerical Solutions

The system (30) is rewritten in the form: where .

Using (33), (31), (32), the fourth-order Runge-Kutta method, and the same initial values in the Table 1, we find the results of the numerical solutions in Table 2. Table 1 is in good agreement with Table 2; that is, the analytical solutions are in full approximation with the numerical ones.

 0 7.048539 1 -2.028946 0 10 4.521935 0.8102441 -5.774878 -0.5860706 20 0.2793855 0.3130168 -7.329234 -0.9497204 30 -4.069078 -0.3029831 -6.102208 -0.9529552 40 -6.87327 -0.8039892 -2.559511 -0.5945567 50 -7.06908 -0.9998797 1.954393 -0.01054169 60 -4.582269 -0.8163247 5.726515 0.5774587 70 -0.3566143 -0.3229907 7.325408 0.9463065 80 4.004261 0.292901 6.144369 0.9560346 90 6.845462 0.7976249 2.63166 0.6029608 100 7.088834 0.9996485 -1.879637 0.02108086 110 4.642086 0.8223141 -5.677522 -0.5687839 120 0.4337946 0.3329276 -7.320768 -0.9427877 130 -3.939007 -0.2827875 -6.185843 -0.9590072 140 -6.816898 -0.7911727 -2.703507 -0.6112964 150 -7.107798 -0.9993057 1.804683 -0.0316163 160 -4.701379 -0.8282109 5.627904 0.5600467 170 -0.5109161 -0.3428261 7.315316 0.9391646 180 3.873325 0.2726439 6.226625 0.9618729 190 6.78758 0.7846335 2.775045 0.6195633 200 7.125968 0.9988525 -1.729537 0.04214693 210 4.760143 0.8340154 -5.577667 -0.5512487 220 0.5879716 0.3526854 -7.309052 -0.9354378 230 -3.807219 -0.2624713 -6.266709 -0.9646314 240 -6.75751 -0.7780079 -2.846264 -0.6277602 250 -7.143345 -0.998288 1.654209 -0.05267154 260 -4.818369 -0.8397263 5.526817 0.5423905 270 -0.6649504 -0.3625042 7.301974 0.9316074 280 3.7407 0.2522707 6.306089 0.9672824 290 6.726693 0.7712965 2.917158 0.6358864 300 7.159923 0.9976127 -1.578707 0.063189

#### 5. Conclusions

In this paper, we studied the singular value of the natural frequency problem in a limiting case when . In this case, the body rotates about the -axis which should coincide with the minor axis of the ellipsoid of inertia. This problem corresponds to the disc motion . The motion takes a slow spin rotation about the symmetric -axis of the disc. The equations of motion and their first integrals are derived and reduced to two quasilinear differential equations of the second order and one first integral . The periodic solutions for this problem are constructed applying the periodicity conditions and assuming a large parameter  proportional to . We used here the large parameter technique instead of the small one well-known in . The advantage of this technique comes from the saving of the high initial energy which is given for the body to start the motion, and the solving of the problem in a new domain of the motion (, , ) and under new considerations. A geometric interpretation using Euler’s angles of motion of the body as a function of time is presented. This interpretation shows the orientation of the disc at any instant of time . There are numerical considerations of the problem and the solutions using the fourth-order Runge-Kutta method which is a specialized and smooth method for finding the approximated periodic solutions of the nonlinear differential equations. Computerized programs are given as a validation of the technique and the results. These programs are carried in a closed interval of time and showed that a full agreement between the analytical solutions and the numerical ones. The agreement of the numerical results with the analytical ones proves that the accuracy (validation) of the resulted solutions and studied techniques (see Tables 1 and 2). From the above tables, we find that the errors between the analytical and the numerical solutions are very small and can be neglected. To study the behavior of the body, we investigate a geometric interpretation of the motion using Euler’s angles. We obtain an arbitrary initial angle of nutation , precession , and pure rotation . We note also that the expressions for the Euler’s angles depend on four arbitrary constants and (where is sufficiently small). Moreover, we note that the disc spins slowly about the minor axis of the ellipsoid of inertia (that is a case of weak oscillations is obtained). In the first approximation, the case of a pseudoregular precession about the vertical axis is attained. As an example, the case of a regular Precession of a slowly spinning Lagrange gyroscope (, ) is obtained as a special case of this motion. There are many applications of these results in both military and civil life. This study is important for the satellite motion which has the correspondence of inertia moments, the antennas, the navigations, the solar collectors, and aerospace dynamics.

#### 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.

#### Authors’ Contributions

I am the individual author of the manuscript.

1. F. A. El-Barki and A. I. Ismail, “Limiting case for the motion of a rigid body about a fixed point in the Newtonian force field,” ZAMM, vol. 75, no. 12, pp. 821–829, 1995. View at: Publisher Site | Google Scholar
2. T. S. Amer, A. I. Ismail, and W. S. Amer, “Application of the Krylov-Bogoliubov-Mitropolski technique for a rotating heavy solid under the influence of a gyrostatic moment,” Journal of Aerospace Engineering, vol. 25, no. 3, pp. 421–430, 2012. View at: Publisher Site | Google Scholar
3. S. V. Ershkov and D. Leshchenko, “On a new type of solving procedure for Euler–Poisson equations (rigid body rotation over the fixed point),” Acta Mechanica, vol. 230, no. 3, pp. 871–883, 2019. View at: Publisher Site | Google Scholar
4. H. M. Yehia, “On the regular precession of an asymmetric rigid body acted upon by uniform gravity and magnetic fields,” Egyptian Journal of Basic and Applied Sciences, vol. 2, pp. 200–205, 2015. View at: Publisher Site | Google Scholar
5. A. H. Nayfeh, Perturbations Methods, WILEY-VCH Verlag GmbH & co. KGaA, Weinheim, 2004.
6. M. Iñarrea, V. Lanchares, A. I. Pascual, and A. Elipe, “Stability of the permanent rotations of an asymmetric gyrostat in a uniform Newtonian field,” Applied Mathematics and Computation, vol. 293, pp. 404–415, 2017. View at: Publisher Site | Google Scholar
7. M. M. Bhatti, R. Ellahi, A. Zeeshan, M. Marin, and N. Ijaz, “Numerical study of heat transfer and Hall current impact on peristaltic propulsion of particle-fluid suspension with compliant wall properties,” Modern Physics Letters B, vol. 33, no. 35, article 1950439, 2019. View at: Publisher Site | Google Scholar
8. L. Zhang, M. B. Arain, M. M. Bhatti, A. Zeeshan, and H. Hal-Sulami, “Effects of magnetic Reynolds number on swimming of gyrotactic microorganisms between rotating circular plates filled with nanofluids,” Applied Mathematics and Mechanics, vol. 41, no. 4, pp. 637–654, 2020. View at: Publisher Site | Google Scholar
9. M. M. Bhatti, A. Shahid, T. Abbas, S. Z. Alamri, and R. Ellahi, “Study of activation energy on the movement of gyrotactic microorganism in a magnetized nanofluids past a porous plate,” Processes, vol. 8, no. 3, p. 328, 2020. View at: Publisher Site | Google Scholar
10. W. S. Amer, “On the motion of a flywheel in the presence of attracting center,” Results in Physics, vol. 7, pp. 1214–1220, 2017. View at: Publisher Site | Google Scholar
11. W. S. Amer, “The dynamical motion of a gyroscope subjected to applied moments,” Results in Physics, vol. 12, pp. 1429–1435, 2019. View at: Publisher Site | Google Scholar
12. T. S. Amer and W. S. Amer, “The rotational motion of a symmetric rigid body similar to Kovalevskaya’s case,” Iranian Journal of Science and Technology, Transactions A: Science, vol. 42, no. 3, pp. 1427–1438, 2018. View at: Publisher Site | Google Scholar
13. T. S. Amer and W. S. Amer, “Substantial condition for the fourth first integral of the rigid body problem,” Mathematics and Mechanics of Solids, vol. 23, no. 8, pp. 1237–1246, 2017. View at: Publisher Site | Google Scholar
14. A. Gilat, “Numerical methods for engineers and scientists,” Wiley, 2013. View at: Google Scholar
15. T. S. Amer, “The Dynamical Behavior of a Rigid Body Relative Equilibrium Position,” Advances in Mathematical Physics, vol. 2017, Article ID 8070525, 13 pages, 2017. View at: Publisher Site | Google Scholar
16. T. S. Amer, “On the dynamical motion of a gyro in the presence of external forces,” Advances in Mechanical Engineering, vol. 9, no. 2, pp. 1–13, 2017. View at: Google Scholar
17. T. S. Amer, M. A. Bek, and M. K. Abouhmr, “On the vibrational analysis for the motion of a harmonically damped rigid body pendulum,” Nonlinear Dynamics, vol. 91, no. 4, pp. 2485–2502, 2018. View at: Publisher Site | Google Scholar
18. I. M. Abady and T. S. Amer, “On the application of KBM method for the 3-D motion of an asymmetric rigid body,” Nonlinear Dynamics, vol. 89, pp. 1591–1609, 2017. View at: Publisher Site | Google Scholar
19. T. S. Amer, “The rotational motion of the electromagnetic symmetric rigid body,” Applied Mathematics & Information Sciences, vol. 10, no. 4, pp. 1453–1464, 2016. View at: Publisher Site | Google Scholar
20. A. I. Ismail, “Solving a problem of rotary motion for a heavy solid using the large parameter method,” Advances in Astronomy, vol. 2020, Article ID 2764867, 7 pages, 2020. View at: Publisher Site | Google Scholar
21. F. L. Chernousko, L. D. Akulenko, and D. D. Leshchenko, Evolution of motions of a rigid body about its center of mass, Springer, 2017.
22. L. D. Akulenko, Y. S. Zinkevich, T. A. Kozachenko, and D. D. Leshchenko, “The evolution of the motions of a rigid body close to the Lagrange case under the action of an unsteady torque,” Journal of Applied Mathematics and Mechanics, vol. 81, no. 2, pp. 79–84, 2017. View at: Publisher Site | Google Scholar
23. A. I. Ismail and T. S. Amer, “The fast spinning motion of a rigid body in the presence of a gyrostatic momentum ,” Acta Mechanica, vol. 154, no. 1-4, pp. 31–46, 2002. View at: Publisher Site | Google Scholar
24. T. S. Amer, “On the motion of a gyrostat similar to Lagrange’s gyroscope under the influence of a gyrostatic moment vector,” Nonlinear Dynamics, vol. 54, no. 3, pp. 249–262, 2008. View at: Publisher Site | Google Scholar
25. A. I. Ismail, T. S. Amer, S. A. el Banna, and M. A. el-Ameen, “Electromagnetic gyroscopic motion,” Journal of Applied Mathematics, vol. 2012, Article ID 812645, 14 pages, 2012. View at: Publisher Site | Google Scholar
26. L. D. Akulenko, D. D. Leshchenko, and F. L. Chernous'ko, “Perturbed motions of a rigid body, close to the Lagrange case,” Journal of Applied Mathematics and Mechanics, vol. 43, no. 5, pp. 829–837, 1979. View at: Publisher Site | Google Scholar

#### More related articles

Article of the Year Award: Outstanding research contributions of 2020, as selected by our Chief Editors. Read the winning articles.