Mathematical Problems in EngineeringVolume 2011 (2011), Article ID 351269, 26 pagesdoi:10.1155/2011/351269

Research Article

Vector Rotators of Rigid Body Dynamics with Coupled Rotations around Axes without Intersection

Katica R. (Stevanović) Hedrih¹ and Ljiljana Veljović²

¹Mathematical Institute, SANU, 11001 Belgrade, Serbia²Faculty of Mechanical Engineering, University of Kragujevac, 34000 Kragujevac, Serbia

Received 27 January 2011; Revised 30 May 2011; Accepted 1 June 2011

Academic Editor: Massimo Scalia

Copyright © 2011 Katica R. (Stevanović) Hedrih and Ljiljana Veljović. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Vector method based on mass moment vectors and vector rotators coupled for pole and oriented axes is used for obtaining vector expressions for kinetic pressures on the shaft bearings of a rigid body dynamics with coupled rotations around axes without intersection. Mass inertia moment vectors and corresponding deviational vector components for pole and oriented axis are defined by K. Hedrih in 1991. These kinematical vectors rotators are defined for a system with two degrees of freedom as well as for rheonomic system with two degrees of mobility and one degree of freedom and coupled rotations around two coupled axes without intersection as well as their angular velocities and intensity. As an example of defined dynamics, we take into consideration a heavy gyrorotor disk with one degree of freedom and coupled rotations when one component of rotation is programmed by constant angular velocity. For this system with nonlinear dynamics, a series of tree parametric transformations of system nonlinear dynamics are presented. Some graphical visualization of vector rotators properties are presented too.

1. Introduction

Well-known toy top or a tern is just a simple toy for many that has the unusual property that when it rotates sufficiently by high angular velocity about its axis of symmetry and it keeps in the state of stationary rotation around this axis. This feature has attracted scientists around the world and as a result of years of research created many devices and instruments, from simple to very complex structures, which operate on the principle of a spinning top that plays an important role in stabilizing the movement. Ability gyroscope that keeps the line was used in many fields of mechanical engineering, mining, aviation, navigation, military industry, and in celestial mechanics.

Gyroscopes’ name comes from the Greek words γυρο (turn) and σκοπεω (observed) and is related to the experiments that the 1852nd were painted by Jean Bernard Leon Foucault. The principle of gyroscope based on the principle of precession pseudoregular.

Gyroscopes are very responsible parts of instruments for aircraft, rockets, missiles, transport vehicles, and many weapons. This gives them a very important role, and they need to be under the strict control of the design and inner workings because in case of damage they could lead to catastrophic consequences. Gyroscope (gyro, top) is a homogeneous, axis-symmetric rotating body that rotates by large angular velocity about its axis of symmetry and is now one of the most inertial sensors that measure angular velocity and small angular disturbances angular displacement around the reference axis.

Properties of gyroscopes possess heavenly bodies in motion, artillery projectiles in motion, rotors of turbines, different mobile installations on ships, aircraft propeller rotating, and so forth. The modern technique of gyroscope is an essential element of powerful gyroscopic devices and accessories, which are used to automatically control movement of aircraft, missiles, ships, torpedoes, and so on. They are used in navigation to stabilize the movement of ships in a seaway, to change direction, and direction of angular and translator velocity projectiles, and in many other special purposes.

There are many devices that are applied to the military, and their design is based on the principles of gyroscopes. Technical applications gyros today are so manifold and diverse that there is a need to get out of the general theory of gyroscopes allocates a separate discipline, called “applied theory of gyroscopes.”

An overview of optical gyroscopes theory with practical aspects, applications, and future trends is presented in [1] written by Adi in 2006.

The original research results of dynamics and stability of gyrostats were given in 1979 by Ančev and Rumjancev [2].

Three of papers [3–5] written by Rumjancev related to stability of rotation of a heavy rigid body with one fixed point in S. V. Kovalevskaya's case, on the stability of motion of gyrostats and Stability of rotation of a heavy gyrostat on a horizontal plane pointed out important research results in this area.

Subjects of series of published papers (see [3–15]) are construction models, dynamics, and applications of gyroscopes as well as special phenomena of nonlinear vibration properties of the gyroscope, analysis of gyroscope dynamics for a satellites, analytical research results on a synchronous gyroscopic vibration absorber, inertial rotation sensing in three dimensions using coupled electromagnetic ring-gyroscopes, gyroscopes for orientation and inertial navigation, and others.

By Cavalca et al. [10] published in 2005 an investigation result on the influence of the supporting structure on the dynamics of the rotor system is presented.

Each mechanical gyroscope is based on coupled rotations around more axes with one point intersection. Most of the old equipment was based on rotation of complex and coupled component rotations which resulting in rotation about fixed point gyroscopes.

The classical book [16] by Andonov et al. contains a classical and very important elementary dynamical model of heavy mass particle relative motion along rotate circle around vertical axis through its centre, whose nonlinear dynamics and singularities are primitive model of the simple case of the gyrorotor, and present an analogous and useful dynamical and mathematical model of nonlinear dynamics.

No precisions and errors appear in the functions of gyroscopes caused by eccentricity and unbalanced gyrorotor body as well by distance between axis of rotations are reason to investigate determined task as in the title of our paper.

This vector approach proposed by us is very suitable to obtain new view to the properties of dynamics of pure classical task, investigated by numerous generations of the researchers and serious scientists around the world.

Using Hedrih’s (see [17–22]) mass moment vectors and vector rotators, some characteristics members of the vector expressions of derivatives of linear momentum and angular momentum for the gyrorotor coupled rotations around two axes without intersection obtain physical and dynamical visible properties of the complex system dynamics.

Between them there are vector terms that present deviational couple effect containing vector rotators whose directions are the same as kinetic pressure components on corresponding gyrorotor shaft bearings.

Also, we can conclude that the impact of different possibilities to establish the phenomenological analogy of different physical vector models (see [17, 20]) expressed by vectors connected to the pole and the axis and the influence of such possibilities to applications allows researchers and scientists to obtain larger views within their specialization fields. This is the reason for introducing mass moment vectors to the rotor dynamics, as well as vector rotators.

The primary-main vector is vector of the body mass inertia moment at the point for the axis oriented by the unit vector and there is a corresponding vector of the rigid body mass deviational moment for the axis through the point (see [17, 20]).

Also, there are a number of the vector rotators, pure kinematics vectors depending on angular velocity and angular acceleration of the body rotation as well as of the mass center position or deviational plane of the body in relation to the axis.

For the case of a rigid body simple rotation about one axis there are two orthogonal vector rotators with same intensity depending on angular acceleration and angular velocity. Directions of these vector rotators are the same as components of kinetic pressures to shaft bearings. The vector rotators correspond to the rotation axis and one in the deviational plane through the axis and second orthogonal to the deviational plane and both with intensity . In the listed papers [17–22] as well as in others, written by first author of this paper, no listed heir, many applications of the discovered vector method by using mass moment vectors are presented for to express kinetic parameters of heavy rotors dynamics as well as of coupled multistep rotors dynamics and for gyrorotors dynamics.

Organizations of this paper based on the vector method applications with use of the mass moment vectors and vector rotators for obtaining vector expressions for linear momentum and angular momentum and their derivatives of the rigid body coupled rotations around two axes without intersections. These obtained expressions are analyzed and series of conclusions are pointed out, all useful for analysis of the rigid body coupled rotations around two axes without intersections when system dynamics is with two degrees of mobility as well as with two degrees of freedom, or for constrained by programmed rheonomic constraint and with one degree of freedom.

By using two vector equations of dynamic equilibrium of rigid body dynamics with coupled rotations around two axes without intersection for two degrees of freedom it is possible to obtain two nonlinear differential equations in scalar form for rotations about each axes and also corresponding kinetic pressures in vector form bearing of both shafts.

2. Mass Moment Vectors for the Axis to the Pole

The monograph [20], IUTAM extended abstract [17], and monograph paper [21] contain definitions of three mass moment vectors coupled to an axis passing through a certain point as a reference pole. Now, we start with necessary definitions of mass momentum vectors.

Definitions of selected mass moment vectors for the axis and the pole, which are used in this paper are as follows.(1)Vector of the body mass linear moment for the axis, oriented by the unit vector , through the point—pole , in the following form (see Figure 1): where is the position vector of the elementary body mass particle in point , between pole and mass particle position .(2)Vector of the body mass inertia moment for the axis, oriented by the unit vector , through the point—pole , in the following form

Figure 1: Arbitrary position of rigid body coupled rotations around two axes without intersection. System is with two degrees of mobility (two freedom or one degree of freedom and one rheonomic constraint) where and are generalized coordinates fixed coordinate system and two moveable coordinate systems and that are rotating with component angular velocities of rigid body coupled rotations: independent generalized (/or rheonomic) coordinates are coordinate of precession rotation and coordinate around self rotation axis. Vector rotators , , and are presented.

For special cases, the details can be seen in [17–22]. In the previously cited references, the spherical and deviational parts of the mass inertia moment vector and the inertia tensor are analysed. In monograph [20] knowledge about the change (rate) in time and the derivatives of the mass moment vectors of the body mass linear moment, the body mass inertia moment for the pole, and a corresponding axis for different properties of the body, is shown, on the basis of results from the first author’s reference [22].

This expression is the vector form of the theorem for the relation of material body mass inertia moment vectors, and , for two parallel axes through two corresponding points, pole   and pole . We can see that all the terms in the last expression have the same structure. These structures are , and .

In the case when the pole is the centre of the body mass, the vector (the position vector of the mass centre with respect to the pole ) is equal to zero whereas the vector turns into so that the last expression (2.3) can be written in the following form:

This expression (2.4) represents the vector form of the theorem of the rate change of the mass inertia moment vector for the axis and the pole, when the axis is translated from the pole at the mass centre C to the arbitrary point, pole O.

The Huygens-Steiner theorems (see [20, 21]) for the body mass axial inertia moments, as well as for the mass deviational moments, emerged from this theorem (2.4) on the change of the vector of the body mass inertia moment at point for the axis oriented by the unit vector passing through the mass center , and when the axis is moved by translate to the other point .

Mass inertia moment vector for the axis to the pole is possible to decompose in two parts: first collinear with axis and second normal to the axis. So we can write

Collinear component to the axis corresponds to the axial mass inertia moment of the body. Second component, , orthogonal to the axis, we denote by the , and it is possible to obtain by both side double vector products by unit vector with mass moment vector in the following form: In case when rigid body is balanced with respect to the axis the mass inertia moment vector is collinear to the axis and there is no deviational part. In this case axis of rotation is main axis of body inertia. When axis of rotation is not main axis then mass inertial moment vector for the axis contains deviation part . That is case of rotation unbalanced rotor according to axis and bodies skew positioned to the axis of rotation.

3. Linear Momentum and Angular Momentum Vector Expressions for Rigid Body Dynamic with Coupled Rotation around Axes without Intersection

3.1. Model of a Rigid Body Rotation around Two Axes without Intersection

Let us to consider rigid body rotation around two axes first oriented by unit vector with fixed position and second oriented by unit vector which is rotating around fixed axis with angular velocity . Axes of rotation are without intersection. Rigid body is positioned on the moving rotating axis oriented by unit vector and rotate around self-rotating axis with angular velocity and around fixed axis oriented by unit vector with angular velocity . Then, axes of rigid body coupled rotations are without intersection. The shortest orthogonal distance between axes is defined by length and it are perpendicular to both axes that is to the direction of angular velocities and . This vector is (see Figure 1): and it can be seen on Figure 1.

When any of three main central axes of rigid body mass inertia moment is not in direction of self rotation axis, then we can see that rigid body is skew positioned. The angles , are angles of skew position of rigid body to the self rotation axis. When center of the mass of rigid body is not on self rotation axis of rigid body rotation, we can say that rigid body is skew. Eccentricity of position is normal distance between mass center and axis of self rotation and it is defined by . Here is vector position of mass center with origin in point , and position vector of mass center with fixed origin in point is .

A plane in which lies the shortest distance, lenght , that is perpendicular to fixed axis of precession rotation by angular velocity is denoted as . A plane that is formed by the shortest distance and fixed axis of component (transmission) rotation gyrorotor system is denoted as and in referent position with we denote axis of fixed coordinate system, with we denote axis in line with axis of component rotation by angular velocity while third axis is perpendicular to it. Lets choose a moveable axis in line to vector , axis that rotates by angular velocity around the moveable coordinate system is rotating as it can be seen on Figure 2.

(a)

(b)

Figure 2: Vector rotators (a) and (b) in relations to corresponding mass moment vectors and , and their corresponding deviational components and as well as to corresponding deviational planes.

In the rigid body, an elementary mass around point we denote with position vector , and with origin in the point on the movable self rotation axis and with vector positions of the same body elementary mass with origin in the point where point is fixed on the axis oriented by unit and is on self-axis rotation oriented by unit and both points are on the end of shortest orthogonal distance betwen axis of body coupled rotations. Position vector of elementary mass with origin in pole is , and velocity of mass particle is: .

3.2. Linear Momentum and Angular Momentum of a Rigid Body Coupled Rotations around Two Axes without Intersection

By using basic definition of linear momentum and angular momentum as well as expresson for velocity of rotation elementary body mass , we can write the following vector expressions: (a)for linear momentum in the following vector form (see [20, 23]): where and are correspond body mass linear moment of the rigid body for the axes oriented by direction of component angular velocities of coupled rotations through the movable pole on self-rotating axis;(b)for angular momentum in the following vector form (see [20, 23]): where and are corresponding rigid body mass inertia moment vectors for the axes oriented by directions of component rotations through the pole on self-rotating axes.

First term in expression (2.6) presents transmission part of linear momentum as if all rigid body mass is concentrate in pole on self-rotating axis and rotate around fixed axes with angular velocity . This part is equal to zero in case when axes are with intersection. Second and third terms in expression for linear momentum present linear momentum of pure rotation, as relative motion around two axes with intersection in the pole on self rotation axes. This two parts are different from zero in all case.

Term is corresponding linear mass moment vector as if all rigid body mass is concentrate in pole on the self rotation axis for the axis oriented by direction of precision rotation, threw the pole .

First term in expression (3.3) presents transmission part of angular momentum as if all rigid body is concentrate in pole on self-rotating axes and rotate arround fixed axis with angular velocity . This part is equal to zero in case when axes are with intersection. First, second, third, and fourth members present transmission parts and fifth and sixth parts present relative angular momentum with respect to pole of pure rotation by two axes as they are with intersection in pole on self axis rotation. In case when axes are with intersection first four members in expression for angular moment are equal to zero.

3.3. Derivatives of Linear Momentum and Angular Momentum of Rigid Body Coupled Rotations around Two Axes without Intersection

By using expressions for linear momentum (3.2) after taking in account derivatives of parts, the derivative of linear momentum of rigid body coupled rotations around two axes without intersection, we can write the following vector expression:

After analysis structure of linear momentum derivative terms, we can see that it is possible to introduce pure kinematic vectors depending on component angular velocitie and component angular accelerations of component coupled rotations that are useful to express derivatives of linear moment in following form

By using vector expressions for angular momentum (4.1) after taking in account derivatives of parts, the derivative of angular momentum of rigid body coupled rotations around two axes without intersection, we can write the folowing expression:

After analysis structure of angular momentum terms, we can see, as in previous chapter for the derivatives of linear momentum, that it is possible to introduce pure kinematic vectors rotators depending on angular velocities and angular accelerations of component coupled rotations and that is used to express derivatives of angular momentum in the following shorter form: where the following denotation is used:

4. Vector Rotators of Rigid Body Coupled Rotations around Two Axes without Intersection

We can see that in previous expression (3.5) for derivative of linear momentum the following three vectors are introduced:

The first two vector rotators and are orthogonal to the direction of the first fixed axis and third vector rotator is orthogonal to the self rotation axis. But, first vector rotator is coupled for pole on the fixed axis and second and third vector rotators, and , are coupled for the pole at self rotation axis and for corresponding direction oriented by directions of component angular velocities of coupled rotations. Intensity of two first rotators is equal and is expressed by angular velocity and angular acceleration of the first component rotation, and intensity of third vector rotators is expressed by angular velocity and angular acceleration of the second component rotation, and are in the following forms:

Lets introduce notation , , and denote difference between corresponding component angles of rotation and of the rigid body component rotations and corresponding absolute angles of pure kinematics vector rotators about axes oriented by unit vectors and . These angles are determined by the following relations:

Angular velocity of relative kinematics vectors rotators , and which rotate about corresponding axes in relation to the component angular velocities of the rigid body component rotations are

In Figure 1. Vector rotators , and are presented.

Fourth vector rotator is in the following vector form and with intensity :

This vector rotator depends on both components of coupled rotations.

We can see that in previous vector expression (3.6) or (3.7) for derivative of angular momentum are introduced following two vectors rotators: and in the following vector form:

The first is orthogonal to the fixed axis oriented by unit vector and second is orthogonal to the self rotation axis oriented by unit vector . Intensity of first rotator is equal to intensity of previous defined rotator and intensity of second rotator is equal to intensity of previous defined rotator defined by expressions (3.7). Their intensities are

In Figure 2 vector rotators (in Figure 2(a)) and (in Figure 2(b)) in relations to corresponding mass moment vectors and , and their corresponding deviational components and as well as to corresponding deviational planes are presented.

Vector rotators and are pure kinematical vectors first presented in [20, 21] as a function on angular velocity and angular acceleration in a form . Also from Section 3.3 expressions (3.5) and (3.6) or (3.7) for derivatives for linear and angular momentum contain members with in tree types of different pure kinematical vectors rotators which rotate around first and second axis in corresponding directions of coupled rotation components, but with pole in or in . These vector rotators are possible to separate by following criteria: (1) intensity of vector rotator is expressed by angular velocity and angular acceleration in the form or angular velocity and angular acceleration in the form and ; (2) intensity of the vector rotators is expressed by both angular velocity components and , and no contain angular accelerations and ; (3) vector rotators are coupled by pole in or in ; (4) type of angular velocities components of vector rotators.

Rotators from first set are rotated around through pole axis in direction of first component rotation angular velocity and depend of angular velocity and angular acceleration . There are two vectors of such type and all trees have equal intensity. Rotators from second set are rotated around axis in direction of second component rotation and depend of angular velocity and angular acceleration . There are two vectors of such type and they have equal intensity.

Let us introduce notation, and denote difference between corresponding component angles of rotation and of the rigid body component rotations and corresponding absolute angles of pure kinematics vector rotators about axes oriented by unit vectors and through pole . These angles are determined by following relations:

Angular velocity of relative kinematics vectors rotators and which rotate about axes in corresponding directions in relation to the component angular velocities of the rigid body component rotations through pole are

Also, it is possible to separate a few numbers of rotators and between the following: where unit vector orthogonal to the axis oriented by unit vector and mass moment vector for the axis oriented by unit vector through pole , and intensity equal twice multiplication of product of intensities of component angular velocities and of rigid body coupled rotations around exes without intersection.

5. Vector Rotators of Rigid Body-Disk Dynamics with Coupled Rotations around Two Orthogonal Axes without Intersection

Let us consider vector rotators for the special case when rigid body-disk rotate around two orthogonal axes without intersection.

Vector of relative mass center position in relation to the pole and self rotation axis oriented by unit vector , we can express in the movable coordinate systems with axes oriented by basic unit vectors: and which rotate around self rotation axis with angular velocity in the form , as well as by basic unit vectors and which rotate around fixed axis oriented by unit vector with angular velocity in the following form: . is angle between mass center vector position and self rotation axis oriented by unit vector . Vector of the orthogonal distance between orthogonal axes without intersection is .

For this case unit vectors and are orthogonal, and after taking into account this orthogonality and corresponding formulas (4.1), (4.5), (4.6), and (4.10) for vector rotators we obtain the following vector expressions:

Previous expressions for vectors rotators are derived with supposition that rigid body is disk and that unit vectors in different deviation planes are:

In Figure 3. four schematic presentations of deviational planes and component directions of the vector rotators of rigid body-disk dynamics with coupled rotation around two orthogonal axes without intersection are presented. In Figure 3(a) deviation plane containing body mass center , vector of relative mass center position in relation to the pole and self rotation axis oriented by unit vector is visible. In Figure 3(b) deviation plane containing self rotation axis oriented by unit vector and body mass inertia moment vector and its deviational component vector of mass deviational moment for self rotation axis and pole is visible. In Figure 3(c) two deviational planes through pole : deviation plane containing self rotation axis oriented by unit vector and body mass inertia moment vector and its deviational component vector of mass deviational moment for self rotation axis and pole and deviation plane containing axis parallel to fixed axis oriented by unit vector and body mass inertia moment vector and its deviational component vector of mass deviational moment for axis oriented by unit vector and through pole are visible. In Figure 3(d) schematic presentation of the rigid body-disk skew and eccentrically positioned on the self rotation axis with corresponding mass moment vectors and deviation plane as a detail of the rigid body-disk coupled rotation around two orthogonal axes without intersection is visible. In all form of the parts in Figure 3. the component directions of the vector rotators components are visible.

(a)

(b)

(c)

(d)

Figure 3: Schematic presentation of deviational planes and component directions of the vector rotators of rigid body-disk dynamics with coupled rotation around two orthogonal axes without intersection. (a) Deviation plane containing body mass center , vector of relative mass center position in relation to the pole , and self rotation axis oriented by unit vector . (b) Deviation plane containing self rotation axis oriented by unit vector and body mass inertia moment vector and its deviational component vector of mass deviational moment for self rotation axis and pole . (c) Two deviational planes through pole : deviation plane containing self rotation axis oriented by unit vector and body mass inertia moment vector and its deviational component vector of mass deviational moment for self rotation axis and pole and deviation plane containing axis parallel to fixed axis oriented by unit vector and body mass inertia moment vector and its deviational component vector of mass deviational moment for axis oriented by unit vector and through pole . (d) Schematic presentation of the rigid body-disk skew and eccentrically positioned on the self rotation axis with corresponding mass moment vectors and deviation plane as a detail of the rigid body-disk coupled rotation around two orthogonal axes without intersection.

By use derived vector expressions of the vector rotators we can obtain some angles between corresponding vector rotator and basic vectors of corresponding movable coordinate systems coupled with corresponding compounding axis of component coupled rotations in the following form:

For the case that , where is relative angle of rotation in comparison with angle of rotation , when is absolute angle of rotor rotation about axis oriented by unit vector , taking into account its rotation about axis oriented by unit vector .

6. Dynamic of Rigid Body Coupled Rotation around Two Orthogonal Axes without Intersection and with One Degree of Freedom

6.1. Model Description of a Gyrorotor Coupled Rotations around Two Orthogonal Axes without Intersection and with One Degree of Freedom

We are going to take into consideration special case of the considered heavy rigid body with coupled rotations about two axes without intersection with one degree of freedom, and in the gravitation field. For this case generalized coordinate is independent, and coordinate is programmed. In that case, we say that coordinate is rheonomic coordinate and system is with kinematical excitation, programmed by forced support rotation by constant angular velocity. When the angular velocity of shaft support axis is constant, , we have that rheonomic coordinate is linear function of time, , and angular acceleration around fixed axis is equal to zero .

Special case is when the support shaft axis is vertical and the gyrorotor shaft axis is horizontal, and all time in horizontal plane, and when axes are without intersection at normal distance . So we are going to consider that example presented in Figure 5. The normal distance between axes is . The angle of self rotation around moveable self rotation axis oriented by the unit vector is and the angular velocity is . The angle of rotation around the shaft support axis oriented by the unit vector is and the angular velocity is . The angular velocity of rotor is . The angle is generalized coordinates in case when we investigate system with one degree of freedom, but system has two degrees of mobility. Also, without loss of generality we take that rigid body is a disk, eccentrically positioned on the self rotation shaft axis with eccentricity , and that angle of skew inclined position between one of main axes of disk and self rotation axis is , as it is visible in Figure 4.

Figure 4: Model of heavy gyrorotor with two component coupled rotations around orthogonal axes without intersections.

(a)

(b)

Figure 5: Intensity of vector rotators and , for different disk eccentricity (a) and for different initial conditions (b).

For that example, differential equation of the heavy gyrorotor-disk self rotation of reviewed model in Figure 4, for the case coupled rotations about two orthogonal axes, we can obtain in the following form: where

Here it is considered an eccentric disc (eccentricity is ), with mass and radius , which is inclined to the axis of its own self rotation by the angle (see Figure 5.), so that previous constants (6.11) in differential equation (6.10) become the following forms:

6.2. Phase Portrait of the Heavy Gyrorotor Disk Coupled Rotations About Two Axes without Intersection and Their Three Parameter Transformations

Relative nonlinear dynamics of the heavy gyrorotor-disk around self rotation shaft axis is possible to present by means of phase portrait method. Forms of phase trajectories and their transformations by changes of initial conditions, and for different cases of disk eccentricity and angle of its skew, as well as for different values of orthogonal distance between axes of component rotations may present character of nonlinear oscillations.

For that reason it is necessary to find first integral of the differential (6.10). After integration of the differential (6.3), the nonlinear equation of the phase trajectories of the heavy gyrorotor disk dynamics with the initial conditions , we obtain in the following for As the analyzed system is conservative it is the energy integral.

6.3. Kinematical Vector Rotators of the Heavy Gyrorotor Disk Coupled Rotations about Two Axes without Intersection and Their Three Parameter Transformations

In the considered case for the heavy gyrorotor-disk nonlinear dynamics in the gravitational field with one degree of freedom and with constant angular velocity about fixed axis, we have three sets of vector rotators.

Three of these vector rotators , , and , from first set, are with same constant intensity and rotate with constant angular velocity and equal to the angular velocity of rigid body precession rotation about fixed axis, but two of these three vector rotators, and are connected to the pole on the self rotation axis, and are orthogonal to the axis parallel direction as direction of the fixed axis. All these three vector rotators , and are in different directions (see Figures 3(a), 3(b), 3(c), and 4). Two of these vector rotators, and , from second set, are with same intensity equal to , and connecter to the pole and orthogonal to the self rotation axis oriented by unit vector and rotate about this axis with relative angular velocity defined by second expression (4.9), , in respect to the self rotation angular velocity . These two of these vector rotators, and are oriented in the following directions:

By use expressions (5.1) we can list following series of vector rotators of the gyrorotor-disk with coupled rotation around orthogonal axes without intersection and with :

One of the vectors rotators from the third set is with intensity and direction: . This vector rotator is connecter to the pole and orthogonal to the axis oriented by unit vector and relative rotate about this axis. Intensity of this vector rotator expressed by generalized coordinate , angle of self rotation of heavy disk, taking into account first integral (6.4) of the differential equation (6.1) obtain the following form: where denotes .

Intensity of two of these vector rotators, and , from second set, depends on angular velocity and angular acceleration . For the considered system of the heavy gyrorotor-disk dynamics, for obtaining expression of intensity of vector rotators, and , from second set, in the function of the generalized coordinate , angle of self rotation of heavy disk self rotation, we take into account a first integral (6.4) of nonlinear differential equation (6.1), and by using these result and previous expressions (6.6) of vector rotator we can write the following.

(i) Intensity of the vectors rotators, and , connected for the pole and rotate around self rotation axis, in the following form:

(ii) Vector rotators orthogonal to the self rotation axes are in the following vector forms:

Parametric equations of the trajectory of the vector rotators and are in the following same forms: but it is necessary to take into consideration that is not in same directions, but is in the same plane orthogonal to the axis oriented by unit vector and through pole .

Relative angular velocity of both vector rotators and in plane orthogonal to the axis oriented by unit vector and through pole . in relation on angular velocity of self rotation, is possible to express by using second expression (4.9), , and we can write the following:

By using previous derived expression (6.8) for intensity of the vectors rotators, and , connected for the pole and rotate around self rotation axis, oriented by unit vector in the orthogonal plane through pole and by changing some parameters of heavy gyrorotor structure, as it is eccentricity , angle of disk inclination , orthogonal distance between axes , as well as parameter contained in the coefficients of the nonlinear differential equation (6.1) and presented by expressions (6.5), we obtain series of the graphical presentation, and some of these are presented in Figure 5.

By using parametric equations, in the form (6.11), of the trajectory of the vector rotators and connected for the pole and rotate around self rotation axis, oriented by unit vector in the orthogonal plane through pole and by changing some parameters of heavy gyrorotor-disk, as it is eccentricity , angle of disk inclination , orthogonal distance between axes , as well as parameter contained in the coefficients of the nonlinear differential equation (6.1) and presented by expressions (6.5), we obtain series of the graphical presentation, and some of these are presented in Figure 6.

(a)

(b)

Figure 6: Transformation of the trajectory of the vector rotator (and ) in the plane through pole and orthogonal to the self rotation axis for different values of parameter .

In Figure 6. transformation of the trajectory (hodograph) of the vector rotator (and ) in the plane through pole and orthogonal to the self rotation axis for different values of parameter is presented.

In Figure 9. transformation of the trajectory (hodograph) of the vector rotator (and ) in the plane through pole and orthogonal to the self rotation axis for different values of parameter is presented.

In Figure 7. transformation of the trajectory (hodograph) of the vector rotator (and ) in the plane through pole and orthogonal to the self rotation axis, for different values of parameter , orthogonal distance between axes of gyrorotor-disk coupled component rotations is presented.

(a)

(b)

Figure 7: Transformation of the trajectory of the vector rotator (and ) in the plane through pole and orthogonal to the self rotation axis, for different values of parameter , orthogonal distance between axes of gyrorotor-disk coupled component rotations.

By using expression, in the form (6.12), of relative angular velocity of the vector rotator (and ) rotation in the plane through pole and orthogonal to the self rotation axis, oriented by unit vector and by changing some parameters of heavy gyrorotor structure, as it is eccentricity , angle of disk inclination , orthogonal distance between axes , as well as parameter contained in the coefficients of the nonlinear differential (6.1) and presented by expressions (6.5), we obtain series of the graphical presentation, and some of these are presented in Figure 8.

(a)

(b)

Figure 8: Relative angular velocity of the vector rotator (and ) in the plane through pole and orthogonal to the self rotation axis, for different values of parameter , orthogonal distance between axes of gyrorotor-disk coupled component rotations.

(a)

(b)

Figure 9: Transformation of the trajectory of the vector rotator (and ) in the plane through pole and orthogonal to the self rotation axis for different values of parameter .

In Figure 8. relative angular velocity of the vector rotator (and ) in the plane through pole and orthogonal to the self rotation axis, for different values of parameter , orthogonal distance between axes of gyrorotor-disk coupled component rotations is presented.

7. Concluding Remarks

First main result presented is successful application the vector method by use mass moment vectors for investigation of the rigid body coupled rotation around two axes without cross-sections and vector decomposition of the dynamic structure into series of the vector parameters useful for analysis of the coupled rotation kinetic properties.

By introducing mass moment vectors and vector rotators we expressed linear momentum and angular momentum, as well as their derivatives with respect to time for the case of the rigid body coupled rotations around two axes without intersections. By applications of the new vector approach for the investigations of the kinetic properties of the nonlinear dynamics of the rigid body coupled rotations around two axes without intersections, we show that vector method, as well as applications of the mass moment vectors and vector rotators simples way show characteristic vector structures of coupled rotation kinetic properties.

Appearance, as it is visible, of the vector rotators, their intensity, and their directions as well as their relative angular velocity of rotation around component directions parallel to components of the coupled rotations, is very important for understanding mechanisms of coupled rotations as well as kinetic pressures on shaft bearings of both shafts.

Special attentions are focused to the vector rotators, as well as to the absolute and relative angular velocities of their rotations. These kinematical vector rotators of the heavy gyrorotor disk coupled rotations about two axes without intersection and their three parameter transformations are done as a second main result of this pepar.

A complete analysis of obtained vector expressions for derivatives of linear momentum and angular momentum give us a series of the kinematical vectors rotators around both directions determined by axes of the rigid body coupled rotations around axes without intersection. These kinematical vectors rotators are defined for a system with two degrees of freedom as well as for rheonomic system with two degrees of mobility and one degree of freedom and coupled rotations around two coupled axes without intersection as well as their angular velocities and intensity.

Acknowledgments

Parts of this research were supported by the Ministry of Sciences and Technology of Republic of Serbia through Mathematical Institute, SANU, Belgrade Grant ON174001 “Dynamics of hybrid systems with complex structures. Mechanics of materials”, supported by the Faculty of Mechanical Engineering University of Niš and Faculty of Mechanical Engineering University of Kragujevac.

References

S. Adi, An Overview of Optical Gyroscopes, Theory, Practical Aspects, Applications and Future Trends, 2006.

A. Ančev and V. V. Rumjancev, “On the dynamics and stability of gyrostats,” Advances in Mechanics, vol. 2, no. 3, pp. 3–45, 1979 (Russian).

V. V. Rumjancev, “On stability of rotation of a heavy rigid body with one fixed point in S. V. Kovalevskaya's case,” Akademia Nauk SSSR. Prikladnaya Matematika i Mekhanika, vol. 18, pp. 457–458, 1954 (Russian).

V. V. Rumjancev, “On the stability of motion of gyrostats,” Journal of Applied Mathematics and Mechanics, vol. 25, pp. 9–19, 1961 (Russian).

V. V. Rumjancev, “Stability of rotation of a heavy gyrostat on a horizontal plane,” Mechanics of Solids, vol. 15, no. 4, pp. 11–22, 1980.

F. Ayazi and K. Najafi, “A HARPSS polysilicon vibrating ring gyroscope,” Journal of Microelectromechanical Systems, vol. 10, no. 2, pp. 169–179, 2001.

A. Banshchikov, Analysis of Dynamics for a Satellite with Gyros with the Aid of to Software Lin Model, Institute for System Dynamics and Control Theory, Siberian Branch of the Russian Academy of the Sciences, Russia.

A. Burg, A. Meruani, B. Sandheinrich, and M. Wickmann, MEMS Gyroscopes and Their Applications, A Study of the Advancements in the Form, Function, and Use of MEMS Gyroscopes, ME-38/ Introduction to Microelectromechanical System.

E. Butikov, “Inertial rotation of a rigid body,” European Journal of Physics, vol. 27, no. 4, pp. 913–922, 2006.

K. L. Cavalca, P. F. Cavalcante, and E. P. Okabe, “An investigation on the influence of the supporting structure on the dynamics of the rotor system,” Mechanical Systems and Signal Processing, vol. 19, no. 1, pp. 157–174, 2005.

W. Flannely and J. Wilson, Analytical Research on a Synchronous Gyroscopic Vibration Absorber, Nasa CR-338, National Aeronautic and Space Administration, 1965.

P. W. Forder, “Inertial rotation sensing in three dimensions using coupled electromagnetic ring-gyroscopes,” Measurement of Science and Technology, vol. 6, no. 12, pp. 1662–1670, 1995.

Y. Z. Liu and Y. Xue, “Drift motion of free-rotor gyroscope with radial mass-unbalance,” Applied Mathematics and Mechanics (English Edition), vol. 25, no. 7, pp. 786–791, 2004.

Y. A. Karpachev and D. G. Korenevskii, “Single-rotor aperiodic gyropendulum,” International Applied Mechanics, vol. 15, no. 1, pp. 74–76.

R. M. Kavanah, “Gyroscopes for orientation and Inertial navigation systems,” Cartography and Geoinformation, vol. 6, pp. 255–271, 2007, KIG 2007, Special issue

A. A. Andonov, A. A. Vitt, and S. E. Haykin, Teoriya Kolebaniy, Nauka, Moskva, Russia, 1981.

K. Hedrih (Stevanović), “On some interpretations of the rigid bodies kinetic parameters,” in Proceedings of the 18th International Congress of Theoretical and Applied Mechanics (ICTAM '92), pp. 73–74, Haifa, Israel, August 1992.

K. Hedrih (Stevanovic), “Same vectorial interpretations of the kinetic parameters of solid material lines,” Zeitschrift für Angewandte Mathematik und Mechanik, vol. 73, no. 4-5, pp. T153–T156, 1993.

K. Hedrih (Stevanović), “The mass moment vectors at an -dimensional coordinate system,” Tensor, vol. 54, pp. 83–87, 1993.

K. Hedrih (Stevanović), Vector Method of the Heavy Rotor Kinetic Parameter Analysis and Nonlinear Dynamics, Monograph, University of Niš, 2001.

K. Hedrih (Stevanović), “Vectors of the body mass moments,” in Topics from Mathematics and Mechanics, vol. 8, pp. 45–104, Zbornik radova, Mathematical institute SANU, Belgrade, Serbia, 1998.

K. Hedrih (Stevanović), “Derivatives of the mass moment vectors at the dimensional coordinate system , dedicated to memory of Professor D. Mitrinović,” Facta Universitatis. Series: Mathematics and Informatics, no. 13, pp. 139–150, 1998.

D. Raškovic, Mehanika III—Dinamika (Mechanics III—Dynamics), Naučna knjiga, 1972.