Industrial Revolution involved scientists to pay the attention to the remarkable gyroscope property expressesed in permanent maintaining the axis of a spinning rotor in a space. Since those time brilliant, famous, outstanding, and ordinary scientists studied gyroscopic effects and published fundamental manuscripts that describe interesting gyroscope properties. The applied theory of gyroscope emerged mainly in the twentieth century due to increasing the angular velocities of the spinning components of the movable machines [1M.N. Armenise, C. Ciminelli, F.V. Dell’Olio, and M.N. Passaro, Advances in Gyroscope Technologies., Springer-Verlag Berlin and Heidelberg GmbH & Co. KG: Berlin, .-4J.B. Scarborough, The Gyroscope Theory And Applications, Fairford., GLOS: U. K., .]. Numerous publications have been dedicated to computing of gyroscopic effects in engineering [5B. Neil, Gyroscope., The Charles Stark Draper Laboratory, Inc.: Cambridge, Massachusetts, .-7H. Weinberg, Gyro Mechanical Performance: the most important parameter, Technical Article MS-2158., Analog Devices: Norwood, MA, pp. 1-5.]. All fundamental textbooks of classical mechanics have chapters that represent gyroscope theory [8R.C. Hibbeler, Engineering Mechanics - Statics and Dynamics., 12th edPrentice Hall, Pearson: Singapore, ., 9D.R. Gregory, Classical Mechanics., Cambridge University Press: New York, .
[http://dx.doi.org/10.1017/CBO9780511803789] ]. There are many publications regarding gyroscope theory as well as many approaches and mathematical solutions that describe gyroscopic properties [10J. Syngley, and J.J. Uicker, Theory of Machines and Mechanisms., 3rd edMcGraw-Hill Book Company: New York, ., 11M.D. Aardema, Analytical Dynamics. Theory and Application., Academic/Plenum Publishers: New York, .
[http://dx.doi.org/10.1007/b116020] ]. However, all known manuscripts contain numerous assumptions, simplifications, as well as explaining the gyroscope effects in terms of the change in the angular momentum only [12R.M. Jonsson, "Gyroscope precession in special and general relativity from basic principles", Am. J. Phys., vol. 75, p. 463.
[http://dx.doi.org/10.1119/1.2719202] -13W.C. Liang, and S.C. Lee, "Vorticity, gyroscopic precession, and spin-curvature force", Phys. Rev. D, vol. 87, p. 044024.
[http://dx.doi.org/10.1103/PhysRevD.87.044024] ]. Some researchers intuitively pointed to the action on the spinning rotor of other inertial forces [14E. Butikov, "Inertial rotation of a rigid body", Eur. J. Phys., vol. 27, pp. 913-922.
[http://dx.doi.org/10.1088/0143-0807/27/4/022] , 15O. Tekinalp, T. Elmas, and I. Yavrucuk, "Gimbal angle restricted control moment gyroscope clusters", Proceedings of 4th International Conference on Recent Advances in Space Technologies (RAST), pp. 585-590.
[http://dx.doi.org/10.1109/RAST.2009.5158259] ]. Nevertheless, the action of these inertial forces on gyroscopes is not presented by mathematical models.

Mathematical models for gyroscope properties in known publications do not match practical applications for gyroscopic devices [16T.J. Quinn, and A. Picard, The mass of spinning rotors: No dependence on speed or sense of rotation, Nature 343, France, 1990, pp. 732 – 735., 17L. Zyga, Gyroscope's unexplained acceleration may be due to modified inertia, PhysOrg.com, July 26, 2011.]. From this, researchers have spawned artificial terms, such as gyroscope effects and others that contradict the physics laws [18J. E. Faller, W. J. Hollander, P. G. Nelson, and M. P. McHugh, Gyroscope-weighing experiment with a null result. Physical Review Letters 64, APS Journal, Vol. 64. No. 8, 990, 1990, pp. 825–826., 19Z. Zhang, J. Sun, and K. Wum, "Error analysis and test study of fibre optic gyroscope north-finder", Proceedings of SPIE, vol. 5634, pp. 611-618.Cardiff
[http://dx.doi.org/10.1117/12.575997] ]. This is the reason that gyroscope theory still attracts researchers seeking to find answers on gyroscope challenges [20K. Kodama, M. Adachi, K. Kamimura, and S. Kurosu, "Gyroscopic Force Measuring System", Transactions of SICE, pp. 117-123., 21J.A. Ferrari, "Gyroscope’s Precession and the Principle of Equivalence", Ann. Phys., vol. 501, no. 5, pp. 399-400.
[http://dx.doi.org/10.1002/andp.19895010513] ]. However, the origin of gyroscope effects is more complex than those represented in the known theories. Recent investigations of the physical principles of gyroscope motions demonstrating that the centrifugal, common inertial, Coriolis forces, and the change in the angular momentum manifest the action of gyroscope’s resistance and precession torques [22R. Usubamatov, K.A. Ismail, and J.M. Sah, "Analysis of a Coriolis Acceleration", Journal of Advanced Science and Engineering Research, vol. 4, no. 1, pp. 1-8.-24R. Usubamatov, "Mathematical model for motions of gyroscope suspended from a flexible cord", Cogent Engineering, vol. 3, p. 1245901.
[http://dx.doi.org/10.1080/23311916.2016.1245901] ]. New mathematical models for inertial torques give significant impact to the gyroscope theory and enable for describing all gyroscope properties. This work considers the action of the inertial torques on the aircraft’s propeller with the different number of blades.

The typical design of an aircraft propeller is represented by the several blades mounted on the cylindrical hub that can be considered as the spinning rotor. The inertial forces acting on the propeller are generated by its centre mass of the hub and centre mass of the blades, whose locations and actions are different. The action of inertial forces is considered on a propeller running with a constant angular velocity of ω in a counterclockwise direction when viewed from the tip of axis oz (Fig. 1). The propeller blade’s centre mass m is located on the circle whose radius is r and the mass centre of the hub is m_h. The blades numbered as I, II, III, IV, V, and VI.

Fig. (1) Schematic of the propeller design with two (a), three (b), four (c) and six (d) blades.

The mass moments of inertia for propellers with several blades are defined by the expression J_p = nJ + J_h, where J is the mass moments of inertia of the blade, n is the number of the blades, and J_h is the mass moments of inertia of the hub. The mathematical models for inertial torques generated by the mass of the propeller’s hub are the same as for the spinning disc [24R. Usubamatov, "Mathematical model for motions of gyroscope suspended from a flexible cord", Cogent Engineering, vol. 3, p. 1245901.
[http://dx.doi.org/10.1080/23311916.2016.1245901] ]. The equations of the inertial torques generated by the blades are defined by the following approach. In uniform circular motions, the blade’s mass moves with acceleration that generates the centrifugal forces (Fig. 2). The inclination of the propeller plane on the small angle Δγ around axis ox leads to the change in the direction of the centrifugal force vector f_ct. The vector f_ct, in which the direction coincides with axis ox (i.e. located on 0^o and 180^o along axis ox), does not change. Other vectors f_ct are located on the inclined plane and have the non-identical change in their own directions. The maximal declination of the vector f*_ct is at 90^o and 270^o from axis ox (Fig. 2). These variable directions of the centrifugal force vector generate the change in the vector’s components f_ct.z, whose direction is parallel to the axle oz. The vector of the centrifugal forces f_ct.z and its variable radius of location relative to axis ox generate the resistance torque T_ct acting opposite to the torque T applied on the propeller. The resistance torque is expressed by the following equation:

Fig. (2) Schematic of acting centrifugal forces, torques and motions of the spinning propeller’s blade.

(1)

where T_ct is the torque generated by the centrifugal force of the rotating blade’s mass, f_ct.z is the axial component of the centrifugal force; and y_m = rsinα is the distance of the blade’s mass location in along axis oy, sign (-) means the action of the torque in the clockwise direction.

The equation of f_ct.z is represented by the following expression:

(2)

where f_ct = mrω² is the centrifugal force of the blade’s mass m, ω is the angular velocity of the propeller, α is the angle of the blade’s mass location, Δγ is the angle of the propeller’s plane inclination (sinΔγ = Δγ for the small value of the angle), other parameters are as specified above.

Substituting the defined parameters into Eq. (1) yields the following equation:

(3)

Then replacing sin²α = (1/2)(1 - cos2α) that is the trigonometric identity and substituting into Eq. (3) yields the following equation:

(4)

The rate change in the torque T_ct per time is represented by the following differential equation:

(5)

where t = α/ω is the time has taken relative to the angular velocity of the propeller and other parameters are as expressed above

Then, the differential of time is: while the expression is the angular velocity of the propeller precession around axis ox. Substituting obtained components into Eq. (5), transformation, separating variables and presenting in the integral form with defined limits the following equation emerges:

(6)

Solving of Eq. (6) giving the following rise

(7)

The value of the resistance torque generated by the centrifugal forces of the rotating blade is changed by sinus law. The resistance torque acts on the upper and lower sides of the propeller with several blades. Then the total resistance torque T_ct of the propellers (Figs. 1a-d) is increased according to the number n of blades and expressed by the following equation:

(8)

where J = mr² is the conventional mass moment of inertia of the propeller’s blade, n is the number of blades, other parameters are as specified above.

The change in the resistance torque T_ct generated by the centrifugal forces of the propellers with n blades (Figs. 1a-d) is represented in Figs. (3a-d).

The equations for maximal and minimal values of the resistance torques for the propellers are represented in Table 1. The value of the resistance torque of the propeller depends on the following components:

Fig. (3) The change in the resistance torque that generated by the propeller with two (a), three (b), four (c) and six (d) blades.

Table 1
Equations of resistance torques generated by the centrifugal forces of propeller’s blades

• The vertical, horizontal or angular location of the blades relative axis,
• The mass and the number of the blades,
• The radius of the location of the blade mass,
• The angular velocity of the propeller and
• The angular velocity of the precession.

The uniform circular motion of the propeller experiences the tangential velocity of the blade’s mass. Under the action of the external torque, a propeller turns around axis ox (Fig. 4) that leads to the non-identical change in the direction of the tangential velocity vectors. The maximal changes in direction have the velocity vectors V* of the blade’s mass located on the line of axis ox (Fig. 4). The two vectors V do not have any changes, whose direction is parallel to the line of axis ox, i.e. located on 90^o and 270^o. These variable directions on the tangential velocity vectors generate the change in the vector’s components V_z whose directions are parallel to the axle oz.

Fig. (4) Schematic of acting inertial forces, torques and motions of the spinning propeller.

The change in the velocity vectors refers to the accelerated motions of the rotating blade’s mass that generate their inertial forces f_in and the inertial torque T_in acting around axis oy that is expressed by the following equation:

(9)

where T_in is the torque generated by the inertial force of propeller mass f_in; a_z is the acceleration of the blade’s mass m along axis oz; and x_m = rcosα is the distance to the mass location along axis ox; other components are represented in section 1.

The equation for the acceleration a_z of the blade’s mass element is defined by the first derivative of the change in the tangential velocity, whose value depends on the angle of its location that is variable with time.

(10)

where a_z = dV_z/dt is the acceleration of the mass along axis oz; V_z = VcosαsinΔ γ is the change in the tangential velocity V of the mass; Δγ is the angle of the turn of the propeller’s plane around axis oy (sinΔγ = Δγ for the small values of the angle); V = rω, ω = dα/dt is the angular velocity of the spinning of propeller; α is the angular location of the blade’s mass; other parameters are as specified above.

Substituting the defined parameters into Eq. (10) yields the following equation:

(11)

Then, trigonometric identity is substituted into Eq. (11) that is presented by the following equation:

(12)

Equation (12) is similar to Eq. (4) with the same final expression:

(13)

The value of the precession torque generated by the rotating blade is changed by cosine law and expressed by the following equation:

(14)

Equation (14) is similar to Eq. (8) with the difference in the trigonometric factor that is replaced by cosα. The total and maximal value of the precession torque acting around axis oy is represented as a sum of the precession torques generated by the inertial forces of the blade’s masses and the change in the angular momentum of the propeller’s hub, in which the expression is as follows:

(15)

where T_p is the precession torque, is the torque generated by the centrifugal forces of the hub, is the torque generated by the change in the angular momentum of the hub [24R. Usubamatov, "Mathematical model for motions of gyroscope suspended from a flexible cord", Cogent Engineering, vol. 3, p. 1245901.
[http://dx.doi.org/10.1080/23311916.2016.1245901] ], all parameters are as specified above.

Equation (15) is generic that enable defining the maximal and minimal values of the precession torques for the propellers based on equations of Table 1. The diagrams of the change in the value of the inertial torque T_in of the propeller with different blades are similar to the diagrams for the torque T_ct generated by the centrifugal forces (Fig. 3). The differences are only the diagrams for inertial torques shifted on 90^o along the abscissa.

The action of the Coriolis acceleration and force is revealed on the propeller under the action of the external torque. Fig. (5) depicts the blade’s mass m that travels with the tangential velocity on the circle of the rotating propeller, which turns on the angle Δγ around axis ox. These motions of the mass produce the Coriolis acceleration and force. The turn of the propeller’s plane around axis ox leads to a non-identical change in the directions of the tangential velocity vectors of the mass. The maximal changes in the velocity vectors V* of the blade’s mass located on axis ox (Fig. 5). The vectors V, whose directions are parallel to axis ox, do not have changes. The variable directions of the tangential velocity vectors generate the change in the vector’s components V_z whose directions are parallel to the axle oz. The changes in the values of the velocity are expressed as an acceleration of the blade’s mass and their inertial forces.

Fig. (5) Schematic of acting Coriolis forces, torques and motions of the spinning propeller’s blade.

The resistance torque generated by the Coriolis force of the blade’s mass is expressed by the following equation:

(16)

where T_cr is the torque generated by Coriolis force f_cr of the rotating blade’s mass m; a_z is the acceleration of the mass along axis oz; and y_m = rsinα is the distance to the mass location along axis oy; the sign (-) means the action in the clockwise direction; other components are presented in Section 1.

The equation for Coriolis acceleration a_z of the blade’s mass is defined by the first derivative of the change in the tangential velocity, whose value depends on the angle of its location on the plane yoz that is variable at the time. The expression for a_z is represented by the following equation:

(17)

where V_z = VcosαsinΔγ(t) is the change in the tangential velocity V = rω of the blade’s mass; Δγ is the angle of turn of the propeller’s plane around axis ox (sinΔγ = Δγ for the small values of the angle); α is the angular location of the mass; ω_x = Δγ/dt is the angular velocity of precession around axis ox; other parameters are as specified above.

Substituting defined parameters into Eq. (16) and transformation yields the following equation:

(18)

Fig. 5 depicts the locations of Coriolis forces of the mass m around axes oz and ox. Substituting the defined parameters into Eq. (18) and transformation yields the following equation:

(19)

Substituting the trigonometric identity into Eq. (19) yields the following equation:

(20)

The rate change in the torque T_cr per time is represented by the following differential equation:

(21)

where t = α/ω is the time taken relative to the angular velocity of the rotating blade, and other parameters are as expressed above.

Then, the differential of time is: while the expression is the angular velocity of the rotating blade. Substituting defined components into Eq. (21), transformation, separating variables and presentation in the integral form with limits of the circle’s quarter yield the following equation:

(22)

Solving Eq. (3.104) yields the following result:

(23)

where the change in the limits of integration for the centroid is taken for the quarter of the circle because its location at other quarters is the same.

Solving of Eq. (23):

giving the following rise:

(24)

The value of Coriolis torque generated by the rotating blade is changed by cosine law and expressed by the following equation:

(25)

The equations for maximal and minimal values of Coriolis torques for the propellers are represented in Table 2.

Table 2
Equations of Coriolis torques generated by the mass of propeller’s blades.

The diagrams of the change in the value of Coriolis torque T_cr acting on the propeller with different blades are similar to the diagrams for the centrifugal torque T_ct (Fig. 3). The differences are only the diagrams for Coriolis torques shifted on 90^o along the abscissa and values are less on the factor π according to Eqs. (8) and (26), respectively.

The total value of the resistance torque acting around axis ox is represented as a sum of the resistance torques generated by the centrifugal and Coriolis forces of the blade’s mass and propeller’s hub.

(26)

where T_r is the resistance torque, T_ct.h and T_cr.h is the resistance torques of centrifugal and Coriolis forces, respectively generated by the propeller’s hub, J_ct.h = 2(π/3)²J_h and J_cr.h = (8/9)J_h is their factored mass moment of inertia [24R. Usubamatov, "Mathematical model for motions of gyroscope suspended from a flexible cord", Cogent Engineering, vol. 3, p. 1245901.
[http://dx.doi.org/10.1080/23311916.2016.1245901] ], all other parameters are as specified above.

In engineering propellers of an aircraft and ships manifest gyroscopic properties. The formulated mathematical models for the inertial torques acting on propellers with different designs are distinguished and proportionally depend on the number of blades, their mass moment of inertia and the angular velocity, as well as the angular velocity of its precession. The inertial torques are generated by the centrifugal, Coriolis and common inertial forces as well as the change in the angular momentum of the propeller’s blades and hub, and manifest the fluctuated and variable resistance and precession torques acting around different axes of the propeller. The derived mathematical models of the inertial torques are new and should be used for the computing of forces acting on a propeller aircraft.

The gyroscope theory is one of the most complex and intricate in terms of analytical solutions. The known mathematical models for gyroscope effects do not match adequately the practical tests and for engineering problem the complex numerical modeling is used. The new mathematical models for the torques acting on the aircraft’s propeller consider the several inertial forces of the rotating blade’s mass and propeller’s hub that manifest the gyroscope properties. Derived mathematical models for inertial propeller’s torques enable for solving problems relating to gyroscopic devices and clearly demonstrate the physical principles of the acting forces and motions.