Open access peer-reviewed chapter

Bifurcation Makes a Wave Resonant Solid-State Gyro Be Stable

By Svetlana Pavlovna Nikitenkova and Dmitry Anatolyevich Kovriguine

Submitted: April 26th 2017Reviewed: October 16th 2017Published: April 18th 2018

DOI: 10.5772/intechopen.71726

Downloaded: 278

Abstract

A golf ball having special dimples flies better than an analogous smoothed one. A surprise is that there is a range of Reynolds numbers for which the turbulent drag is somewhat less than that in the laminar case. Analogies always meet together and are accomplished themselves in the physics. We have traced an effect similar to the above mentioned in the area of solid mechanics when the nonlinear system passes through a sequence of bifurcations. In mechanical engineering, the role of such a system can play a solid-state wave gyro entering the family of MEMS/NEMS. It is known that a circular Foucault pendulum can serve as an angular sensor. Standing waves in a thin-walled elastic axisymmetric resonator of a solid-state wave gyro, mounted on a rotating platform, can also detect a rotation rate. Because there are no typical mechanical parts there, such wave sensors have advantages for long-term space missions. However, to maintain the functionality and sensitivity of a wave gyro in practice, the driving of standing waves requires a sophisticated feedback control. Nonetheless, we have demonstrated that such a gyro can operate without any feedback at the expense of the natural nonlinearity of the resonator in a postbifurcation regime.

Keywords

  • wave solid-state gyro
  • Foucault pendulum
  • spring pendulum
  • bifurcation
  • stability

1. Introduction

It is well known that the aerodynamic drag is proportional, first of all, to the kinetic energy of a flow. Moreover, the drag depends on a lot of other almost uncontrolled factors such as the shape, size, inclination of the object in the flow, etc. Implicitly, all these are accumulated all together in the following formula for the hydrodynamic resistance: cdρFV2/2, where ρis the density of a fluid or gas, Vis the velocity of a flow, Fdenotes an effective area of the object, cdstands for the so-called drag coefficient being a dimensionless characteristic function versus the Reynolds number. Recall that the Reynolds number is usually defined as Vl/ν, where lis an effective length, νis the kinematic viscosity coefficient. The function cdis respective for all the uncontrolled points accompanying the interaction between the body and the flow. This function can be recorded experimentally using a wind tube. Experimental results obtained are used to evaluate the drag on other similar objects. In particular, the similarity by the Reynolds number ensures a correct description of viscous motions. Notice that the drag coefficient is nearly a constant in a wide range of Reynolds numbers in most practical cases. As long as the flow velocity would increase, the Reynolds number also increases. It is well known that in a slow flow, the viscosity can be neglected in practice. In this case, one can observe almost an ideal flow with no boundary layer near the surface. Therefore, there is almost no resistance to the motion in accordance with well-known d’Alembert paradox. As the velocity increases further, the drag becomes non-zero as a result of vortex generation. Then, these vortexes oscillate to produce the so-called Kármán vortex street, while the drag increases gradually together with the flow speed. At even higher velocity, the boundary layer turns into a chaotic turbulent flow. It is natural to expect that the turbulent drag would be higher than that in the laminar flow case. Nonetheless, a surprise is that there is a range of Reynolds numbers for which the turbulent drag is somewhat less than that in the laminar case. It is observed experimentally that a roughened cylinder or a ball will pass the turbulent flow at a lower Reynolds number than a smooth cylinder or a ball. At first sight, this results in a rather paradoxical result: there is a small range of Reynolds numbers for which the drag of a roughened body is less than the drag of a smooth body of the same size. For a good example, we can recall a golf ball having special dimples. Indeed, the golf with those dimples flies better than an analogous smoothed one [1].

Physical analogies meet together and are accomplished themselves in most applications. Can we find out an effect similar to the abovementioned dimpled ball that flies better in the area of solid mechanics? We know that the hydrodynamic drag evolves throughout the sequence of bifurcations. But any bifurcation is inherent in nonlinear systems. Therefore, our study should deal with nonlinear mechanical systems. In mechanical engineering, the machines, sensors, and other useful devices are of interest. One of them is a solid-state wave gyro from the family of microelectromechanical system/nanoelectromechanical system (MEMS/NEMS) devices. Let us briefly touch some key points of this topic.

In the nineteenth century, people held no doubt that the Earth is turning about its axis. However, there were no direct proofs. For this reason, Léon Foucault had managed in 1851 his famous experiment with the giant pendulum. The heart of his experiment is the inertia of the pendulum or the resistance to change the motion. This means that the swing plane remains to be fixed in Newton’s absolute frame of references; moreover, this evidence can be traced with the naked eye due to the slow clockwise veering since the Earth rotates anti-clockwise. However, Foucault’s pendulum has one essential flaw, namely the dependence on the latitude and its giant size suppressing the energy dissipation for a while. That is why, in 1852, though an idea, we should note, was not new, Foucault produced a compact and convenient gyroscopic setup based on the immobility of the axis of a rotating mass. It was a prototype of a conventional gyro providing needed horizons at sea and then in the air. The irony of fate, when the twentieth century was ending, the old idea came back to use a Foucault-pendulum-based device as a high precision gyroscope, although, in the other role of a thin axisymmetric solid-state resonator, associated with satellite guidance systems intended for long-term missions extending up to 15 years [6].

In 1890, Bryan demonstrated that the revolving of the standing wave nodes records effectively the rotation of the elastic resonator, as equally as a material point tends to conserve the spatial position in Newton’s space. This evidence was not so new. Yet before Bryan’s experiments, it was well known that the plane of transverse vibration of a straight wire will remain fixed in space instead of turning, though the wire rotates slowly about its axis during oscillations. Bryan noticed that when the vibrating body is such as a bell, rotation about its axis will produce an intermediate effect: the nodal meridians revolve with angular velocity less than that of the resonator that exhibits a new finding which yet nothing is easier than to be verified experimentally. He selected a champagne glass, then struck it to get a pure tone and when the glass turns around, Bryan heard sound beats demonstrating that the nodal meridians do not remain fixed in Newton’s absolute space. He evaluated that the nodal angular velocity is about 3/5 of this almost hemispherical resonator [7]. Although, only in the second half of the twentieth century, Bryan’s effect gets a wide extent turning into a concept of a prototype of an angle sensor that possesses a lot of advantages compared with a conventional gyro because its core is a monolith [8, 9, 10]. Moreover, if the power of the wave gyro is lost for a short time, the resonator conserves the angular rate, so that when the power returns, the gyro need not be reinstalled, etc., unlike, for instance, optical gyroscopes.

Let us now recall some mathematical aspects of the problem along the natural evolution of ideas. The role of the most simple mechanical system, appearing as an abstract oscillatory gyro in vitro, can play a circular Foucault pendulum. The corresponding mathematical model can be given by the following two differential equations:

x¨2Ωsinφẋ+ω2x=0;y¨+2Ωsinφẏ+ω2y=rΩ2sinφcosφ,E1

where φdefines the latitude of geographical place; Ωis the Earth’s angular velocity in absolute value; ris the Earth’s radius plus a distance from the point where the pendulum is suspended. It is supposed that Earth rotates about the axis z. Here, xtand ytdenote the pair of projections of infinitesimal oscillations of the bob, weight mg, in the rotating Cartesian frame of references Oxyz. The origin Obelongs to the Earth’s center while the plane Oxzpasses through both poles. For operating such a gyro in vita, the pendulum should be excited to vibrate. Moreover, the amplitude of oscillations has to be maintained on some sensitive level because of energy dissipation. The presence of damping requires permanently pumping the energy into these oscillations by external forces for a permanent operating.

Almost a half century later, Bryan has derived his equation describing two similar vibration forms which, traveling toward a thin rotating ring, produce phenomena of beats, and which, in the case of high-frequency oscillations like those of sound, can detect effects of slow rotation about the sensitivity axis:

vt,tvθ,θ,t,t4ωvθ,t+ω2μvθ,θ+T2vθ,θ+vθ,θ,θ,θσa2=βvθ,θ,θ,θ,θ,θ+2vθ,θ,t,t+vθ,θσa4.E2

Here, v=vθtand wθtare the tangential and radial components of displacements measured in the rotating polar frame of references; a+wθ, at time t. The symbol ωdenotes the angular rate; ais the radius of the ring; σis mass density; T=σa2ω2μstands for an attractive force μ, times the distance, directed toward the axis. We can consider this model of the circular ring as a prototype of a solid-state wave gyro to trace the precession of the flexural standing wave in vacuo, that is, in the absence of damping. However, the energy dissipation always presents in nature. Therefore, this system requires some external feedback control to supply the wave motion. The principal parametric resonance is a good way to excite such a gyro by driving the axisymmetric tension T. Nonetheless, some feedback is in need since the parametrically excited oscillations are always unstable either in the presence or in the absence of damping [8]. To maintain the functionality as well as the sensitivity of a conventional wave gyro in practice, the driving of standing waves requires a sophisticated feedback control to eliminate the unstable oscillations caused by their parametric excitation. Therefore, such a kind of wave excitation looks as not quite satisfactory from the viewpoint of actuation, driving or stability properties. Since the linear theory is not able to achieve the desired goal, it is necessary to look for other theoretical tools within the nonlinear wave dynamics. Our thought is to use the exclusive property of the axisymmetric mode of the ring to control the amplitude of flexural waves in the gyro in the presence of energy dissipation, just to use this one as a sort of catalysator in the dynamical process. This chapter demonstrates that when both the primary resonant pumping over the axisymmetric mode and the principal parametric resonant excitation are combined, such a gyro can operate without any feedback, just at the expense of the natural nonlinearity of the resonator in a post-bifurcation regime.

2. Spring pendulum

Before all the theoretical approaches related to a wave resonant solid-state gyro, it may be of place to consider a most simple mechanical system which behaves analogously and contains no excess detail, in order to explain from first principles why the nonlinearity stabilizes unstable parametric oscillations, at the expense of resonance experienced in the mechanical system. The best candidate seems to be a spring pendulum.

Therefore, let us consider a pendulum swinging in a plane and consisting of a bob of mass mattached to a spring with the stiffness kand of the natural length l. Under the Earth gravity, the equilibrium length of the pendulum becomes l+Δ, where Δ=mg/kdenotes the spring elongation. After introducing the time-dependent radial coordinate rt, perpendicularly to the pivot of the spring, and the angle θtbetween the spring and the vertical, the equations governing the motion can be derived easily by using the following Lagrangian function:

L=mṙ2+l+r2θ̇2ṙ2+l+r2θ̇2/2+kr2/2+mgl+r1cosθr.E3

Alternatively, we may follow the Hamiltonian approach:

H=p12m+p22ml+q1212mp12m2+p22l+q12m2+kq122+mgl+q11cosq2q1,E4

and then writing the following set of Hamiltonian equations:

q̇1=p1m;q̇2=p2ml+Δ+q12;ṗ1=p22ml+Δ+q13kq1+mgcosq21;ṗ2=mgl+Δ+q1sinq2,E5

where q1=r, q2=θ, p1=mṙ, p2=ml+r2θ̇are the canonical coordinates. For further preparations, it is convenient to pass from the physical coordinates to dimensionless variables: t=τ/Ω1; p1=lmΩ1z3; p2=l2mΩ12+Ω222z4/Ω13; q1=lz1; q2=z2, where Ω1=k/mand Ω2=g/ldenote the natural radial and angular frequencies of the pendulum with an inextensible thread at small oscillations. Notice that the angular frequency would be somewhat lower in the case of extensible thread, namely ν=Ω2/Ω12+Ω22. In the new dimensionless variables, equations of infinitesimal oscillations would read

dz1=z3;dz2=z4;dz3=z1;dz4=ν2z2.E6

Let us now consider small-but-finite oscillations of the pendulum. In this case, the general solution to the set (6) can be represented as follows:

z1τ=μA1Tcosφ1T+τ;z2τ=μA2Tcosντ+φ2T/ν;z3τ=μA1Tsinφ1T+τ;z4τ=μA2Tsinντ+φ2TE7

to look for solutions to Eq. (5) with unknown amplitudes Aiand phases φislowly evolving in the new slow time T=μτ, where μis arbitrary small-scaling parameter. The substitution from Eq. (7) into Eq. (5), and ordering in μ, yields the following set of first-order approximation evolution equations:

dA1dT=A222sinΦ;dA2dT=9A2A132sinΦ;dΦdT=9A1162A22A1cosΦ,E8

under the one assumption that the following integer-valued ratio between the frequencies, manifesting the primary parametric resonance, that is, 1:ν=1:2, takes place. Notice that the set (8) contains the so-called generalized phase ΦT=φ1T2φ2T, as the unknown variable. The set (8) is of Hamiltonian structure. So, the related average Hamiltonian function reads

A22TA1TcosΦT=A220A10.E9

Moreover, Eq. (8) posses one more additional integral of motion:

E=9A1232+A222,E10

where E=9A120/32+A220/2is the average kinetic energy of the pendulum defined at the initial instant of time. These two integrals, Eqs. (9) and (10), allow us to integrate the set (8) analytically in terms of Jacobi elliptic functions [3].

Let us now refer to the damped forced motion of the pendulum. In this case, equations governing the motion (5) can be easily modified to the following form:

q̇1=p1m;q̇2=p2ml+Δ+q12;ṗ1=p22ml+Δ+q13kq1+mgcosq212μδp1μFcosΩ1t+ψ1;ṗ2=mgl+Δ+q1sinq22μδp2,E11

where δis the viscous drag coefficient; Fis the amplitude of an external harmonic force acting at the resonant frequency of free radial oscillations of the pendulum. We can utilize the same anzats (7) to derive the truncated set of modulation equation for slowly varying amplitudes and phases. Small parameter μemphasizes that both the damping and forcing are small but finite. After the substitution from Eq. (7) into Eq. (11), and ordering in μ, the first-order approximation evolution equations describing the resonant excitation of the pendulum over the radial mode can be written as

dA1dT=A22sinΦ/2+GΛA1;dA2dT=9A2A1sinΦ32ΛA2;dΦdT=A222A1+9A116cosΦ.E12

These equations take in a natural manner into account the so-called triad-angle locking phenomenon described in [2], when both the phase of external force and the phase of radial oscillations have to be coupled, that is, φ1T+ψ1=π/2, accordingly the phase matching conditions. Notice that the parameter G=F/6mg, entering Eq. (12), is interpreted as the dimensionless force, while Λ=δ/Ω1stands for the dimensionless drag coefficient.

In the case of stationary motion, when the amplitudes are constants, Eq. (12) produces the following algebraic set of equations:

A22/2+GΛA1=0;9A2A132ΛA2=0;ΦT=π/2,E13

having the following solutions:

A1=GΛ,A2=0,whenK32Λ29;E14

and

A1=32Λ9,A2=18G64Λ2/3,asK>32Λ29.E15

These solutions, Eqs. (14) and (15), are plotted in Figure 1. As we can see, the dimensionless parameter Kplays the role of control parameter governing the bifurcation of the system. The stationary solution (14) near the bifurcation point K¯=32Λ2/9becomes unstable and gives way to the new stationary steady state (15).

Figure 1.

Bifurcation diagram—Amplitudes of radial and angular modes versus control parameter K. Solid circles refer to stable angular oscillations, circles mark stable radial oscillations. Crosses denote unstable radial mode. Change in the stability takes place near the bifurcation point K¯ (parameter Λ=3/82 has been scaled arbitrary for the best view).

To examine the stability properties without any excess mathematical preparations, let us analyze the power of drag forces using the energy balance resulted from the set (12):

ddT9A1216+A22=Λ9A1216+A22+9GA116.E16

For the stationary orbits, the left-hand term of this equation is zero. Since the right-hand term, representing changes of the average kinetic energy in time, is also zero, then the power of drag forces should be, in turn, proportional to the kinetic energy with a negative sign. The diagram of the average kinetic energy as the function of control parameter Kis shown in Figure 2. As we can see, the power of drag forces reduces after the bifurcation under the same external harmonic excitation. This means stability.

Figure 2.

Bifurcation diagram—Average kinetic energy versus control parameter K. Solid lines refer to stable oscillations while dotted line marks unstable regime of oscillations. Parameter Λ is the same as in Figure 1.

Let us now refer to the original Eq. (11) to verify results yielded by the modulation theory. To plot the numerical results, we have used the following transform: t=τm/k; p1t=lmk/mz3τ; p2t=mk/mmg+lk2z4τ/k2; q1t=lz1τ; q2t=z2τ, from the physical to dimensionless variables. The numerical parameters of the system are chosen by us as follows: F=1.0N; β=0.01kg/s; δ=0.01s1; g=9.8m/s2; k=29.40kg/s2; l=1.0m; m=1.0kg.

Since we are studying a damped forced motion, the initial conditions to the governing equations rewritten in terms of either the physical or the dimensionless variables ziτcannot play any role at large times, though we have decided to start in our calculations from the static equilibrium point. Figure 3(a) displays the time history of the radial coordinate of the system governed by Eq. (11), rewritten in new dimensionless variables ziτ, while the evolution of the angular oscillations is shown in Figure 3(b). Figure 4 displays a typical Lissajous curve that appears on z1z2cross section of the four-dimensional phase space for the fully developed stationary state. Figure 5 shows individual cross sections related to the radial and angular vibrations in the developed stationary regime of oscillations. Finally, Figure 6 exhibits the evolution of the Hamiltonian function (4) or the energy of the system (11) in comparison with that one where the angular coordinate is fixed to zero but the force excitation remains the same. Notice that if the angular coordinate is fixed, then the pendulum just performs linear oscillations in the radial direction.

Figure 3.

The evolution of the damped forced radial coordinate of the pendulum in time τ (a) and the time history of the resonantly excited angular coordinate of the pendulum in time τ (b).

Figure 4.

The evolution of the damped forced radial coordinate of the pendulum in time τ.

Figure 5.

Cross sections of the phase space demonstrating individual dynamical behavior over the radial (a) and angular coordinates (b) at the damped forced oscillations in a spring pendulum in time τ.

Figure 6.

Energy of the pendulum becomes some lower in comparison with the same system with fixed angular coordinate. Dashed line refers to a spring pendulum with constrained to zero angular coordinate.

In new dimensionless variables ziτ, the energy of the oscillating pendulum is given by the following explicit expression:

Eτ=118z1+1254z1+13cosz2+81z14+162z13+81z32+135z12+162z32+108z1+81z32+256z42+54.E17

Mathematically, the above results represent direct proofs that in our case the modulation theory has more advantages in comparison with numerical investigations. At small and even moderate oscillations, this asymptotic approach allows us to make up a complete parametric analysis of the system under the study. Pragmatically, stability properties demonstrated using the above example of a spring oscillator in post-bifurcation regime can be used when designing an idea of a high-precision angular sensor which can appear in the form of a circular Foucault pendulum, in most simple case, or as a solid-state wave gyro with axisymmetric resonator in more general case (for instance, see [4, 5] and references therein).

Gradually, one may always construct an appropriate theory by further developing the model given by Eq. (1). Nonetheless, here we decided to invoke more general case in our consideration, namely a thin-ring resonator of which is described in detail in [11].

3. Thin-ring resonator

In order to define the in-plane position in a thin-ring resonator, thickness hand radius R, rotating with angular velocity Θabout its sensitive axis, we introduce the two frames of polar coordinates in the absolute Newton’s space R+wϕ+v/Rand in the rotating plane R+wφ+v/Rwhere both have the same center as pole. Here, wand vare radial and tangential displacements of the ring, respectively. Therefore, we shall have

ϕ=φ+Θt,

if the ring rotates uniformly in time t. Also, we may suppose that the angular rate as a slowly varying function of time. In this case, the latter expression is slightly modified by the following integral:

ϕ=φ+t0tΘτ.

Equations of motion are derived from the theory of thin-walled shells that uses Kirchhoff-Love hypotheses. So, the field of displacements in the ring is expressed as us=vζwsv/Rand uζ=w, where v=vstand w=wstare the components of displacements, rewritten as functions of the circumferential coordinate s; ζis the distance from the mid-line along the radius.

In the rotating frame of references, the Lagrangian density of the system reads

L=ρF2vt+Θw+RΘ2+wtΘv2EFκ22R2v212h/2h/2Eϵss2,

where Eis Young’s modulus; Fdenotes cross-section square; ρis mass density; κcharacterizes the stiffness of the rotating platform of the ring; ϵss=vs+w/R+ξws,svs/R+ws2/2is the circumferential component of the deformation tensor. According to paradigms of the variational analysis, equations governing the motion in the ring have the following form:

Lv̇t+LvssLv=QvRv̇;Lẇt+LwssLws,ss,sLw=QwRẇ,E18

where Qvand Qware introduced as generalized forces; the function R=ηK, expressed through the kinetic energy K, is respective for the energy dissipation model linearly scaled by the coefficient η. For further preparations, it seems to be continent to rearrange these equations to dimensionless notation:

v¨+2μΩẇ+μΩ̇w+Ω̇/ϵμ2Ω2vVφ+ϵ2Wφ,φ+κ2v=ϵμ2wφ2φ+Qv;w¨2μΩv̇μΩ̇vμ2Ω2wμΩ2/ϵ+V+ϵ2Wφ,φ,φ=ϵμVwφφ+wφ22+ϵ2μ22wφ3φ+Qw,E19

where ϵ=h/12R1is the relative thickness of the ring; small parameter μ=a/R1is intended for making up further procedures of the perturbation analysis; the functions V=vφ+wand W=wφvare written for brevity. The transform to these dimensionless variables reads: vτφvst/μa; wτφwst/μa; φ=s/R; τ=tc/R; Ωτ=ΘtR/μc. Here, c=E/ρdenotes the wave propagation velocity; a=h/12.

Finally, we should provide the set (19) by the periodicity conditions

vφτ=vφ+2πτ;wφτ=wφ+2πτ.E20

3.1. Dispersion relation

For studying the wave propagation in the rotating frame of references, it is convenient to define linear modes of oscillations in the ring on the fixed platform. If the rotation is absent, then the normal modes of vibrations represent standing waves being a superposition of two waves traveling toward identical wave numbers, frequencies as well as amplitudes. First of all, we should understand how the spectrum of these oscillatory modes would change in the uniformly rotating ring. Are there the standing modes taking place? And if these cannot be detected, then, what thing should we invoke instead into our study? We know that the precession of waves appears as a reaction on the rotation. Therefore, there will be expected some asymmetry in the polarization vectors of traveling waves. One more finding is we know that the precession rate of waves has to be proportional to the difference between the frequencies in the wave pair as if those compose a standing wave in the ring being at rest. Let the angular rate in the ring be constant and then the wave precession should appear as some sort of kinematic reactions on the rotation of the ring platform. But if the platform rotates with acceleration, then we can expect a kind of some dynamical response that may appear as essentially differ from kinematic one. Fortunately, experiments with revolving axisymmetric bodies tell us that the expression for the uniformly rotating ring, which will be done subsequently, is equally valid in both cases of the uniformly rotating ring and in the motion with arbitrary, but moderate, angular acceleration [8].

In the case of uniform rotating, a simple solution of the linearized set (19) is given by

vτφ=μΩ2ϵμ2Ω21+Bexpi+ωτ;wτφ=Aexpi+ωτ,

Here, we can trace the appearance of the constant extension in radial direction caused by the centrifuge force. Notice that the amplitudes, Aand B, entering therein, are linearly interrelated, that is, B=pA, through the interrelation coefficients, p, defined for both the high- and low-frequency sets of normal modes:

pk,n=i2μΩωk,n+n1+ϵ2n2n21+ϵ2+κ2ωk,n2+μ2Ω2.E21

The high- and low-frequency branches are indexed by k=1and k=2, correspondingly. Recall from the linear algebra that these coefficients should satisfy the orthogonality condition, that is, p1,np2,n=1, for arbitrary wave number n, excluding the one case of linearly decoupled oscillations, taking place at n=0, that is, at axisymmetric radial oscillations. The natural frequencies of waves, ωk,n, are defined by the dispersion relation

n21+ϵ2+κ2ω2+μ2Ω21+ϵ2n4ω2+μ2Ω22μΩωn1+ϵ2n22=0.E22

Let the angular velocity be zero and then the dispersion relation turns into the simple algebraic equation

n6ϵ2ϵ2ω22+κ2n4+ω21ϵ2+ω2n2ω4+1+κ2ω2κ2=0,E23

Two identical ones of all the four roots of Eq. (23) are shown in Figure 7. The low-frequency branch refers in general to the bending modes while the high-frequency ones reply mainly to the circumferential modes of extension.

Figure 7.

The high- and low-frequency dispersion branches for waves on the fixed platform (ϵ=0.01, κ=0.45, Ω=0).

In the case of a small angular rate, that is, Ωωk,n, the roots of the dispersion relation (22) are approximately represented as follows:

ωk,nΩωk,n0+2μΩn1+ϵ2n21+ϵ2n4+1+ϵ2n2+κ22ωk,n02.E24

It may not, perhaps, be out of place to note that these roots possess by asymmetry because of the angular rate. Nonetheless, each frequency, ωk,n, can be decomposed on the antisymmetric term as well as on the symmetric part. The antisymmetric term is respective for the wave precession caused by the ring rotation while the symmetric one describes a small inessential correction.

4. Nonlinear resonant coupling between the axisymmetric radial oscillation and two similar bending wave forms being in phase

This section invokes briefly some general points from the common theory of waves, including understanding phenomena of resonance experienced in nonlinear mechanical systems and some mathematical preparations which are necessary to investigate the given problem from the viewpoint of the perturbation analysis, using methods of slowly varying amplitudes and phases. Then, we generalize the problem by considering the case of damped forced oscillations in a thin circular ring. This generalization of the problem leads to a thought on how to excite stable wave precession regimes which need no feedback to repair unwanted motions in the solid-state wave gyro.

4.1. Triple-mode resonant coupling between waves in the ring

Truncated equations or evolution equations describing the modulation phenomena of amplitudes and phases can be directly obtained using the following anzats [12, 13, 14, 15]:

vτφ=k=13pkAkTexpiϕk;wτφ=k=13AkTexpiϕk,E25

which represent an approximated solution to Eq. (19). Here, AkT(k=1..3) are the slowly varying complex wave amplitudes, evolving at slow time T=μτ; ϕ1=+ω1τ, ϕ2=+ω1τ, ϕ3=ω3τare the fast rotating phases; denote the complex conjugation. Notice that the amplitude and phase number three refer to the axisymmetric radial oscillation of the ring.

We assume that the angular rate and acceleration of the ring are small enough, that is, Θτμand dΘ/μ2. Also, we should suppose that spectral parameters of waves in the solution (25) are matched by the following phase synchronicity conditions:

ω3=ω1+ω2+μΔω,E26

where Δωdenotes a small discrepancy between the phases of the high-frequency axisymmetric radial oscillation and two low-frequency quasi-harmonic waves with amplitudes, numbers one and two. Certainly, all the frequencies ω1, ω2, ω3and the wave numbers ±nshould satisfy the dispersion relation (23). In the theory of nonlinear waves, such a trio is used to name a resonant triplet while the synchronicity (26) is associated with the so-called phase-matching conditions.

One more convenient way of obtaining the evolution equation is walking down along the Hamiltonian formalism. Let us consider the average Lagrangian of the system

L=12π302π02π02πLdϕ1dϕ2dϕ3,E27

in order than to expand this one in a formal series in small parameter μas follows:

L=L0+μL1+μ2L2+,

Here, zero-order approximation term, that is, L0=0, would coincide exactly with the dispersion relation, accordingly to findings of Whitham [16]. The second term L1, entering this expansion, describes effects arising in the first-order nonlinear approximation analysis which we need to study here. Then, the nontrivial average Hamiltonian of the system would read

H=k=13Ak,TL1Ak,T+A¯k,TL1A¯k,TL1,

In this case of resonant coupling between the axisymmetric radial oscillation and two low-frequency bending waves, being in phase (26), this average Hamiltonian is rearranged in more pragmatic form

H=2iΩj=12ωjpj¯pjAj2+ϵn2A¯1A¯2A3expiΔωT+A1A2A¯3expiΔωT.E28

Recall that the bending wave frequency ω=ω1=ω2should be close to the frequency ω3/2that denotes the frequency of the axisymmetric radial oscillation, that is, ω3=1. Note that the case under consideration is associated with the principal parametric resonance. His Hamiltonian produces a set of following evolution equations:

dA1dT=2Ωp1¯p11+p2A1iϵn2ω1+p2A¯2A3eiΔωTdA2dT=2Ωp2¯p21+p2A1iϵn2ω1+p2A¯1A3eiΔωT;dA3dT=iϵn2A1A2eiΔωT,E29

where pk=1kin1+ϵ2n2/1+ϵ2n2ω2are the wave polarization coefficients which satisfy the equality p=p1=p2, due to the symmetry of the problem.

Now, we use the following transform of variables:

A1=a1exp2p1¯p11+p20TΩς;A2=a2exp2p2¯p21+p20TΩς;A3=a3,E30

that allows us to rewrite the set (29) by getting rid of Coriolis terms:

da1dT=iϵn2ω1+p2a¯2a3eiΔωT;da2dT=iϵn2ω1+p2a¯1a3eiΔωT;da3dT=iϵn2a1a2eiΔωT,E31

Eq. (32) have now obtained a standard form similar to Euler’s kinematic equations describing rotation of a rigid body with a fixed point. This set can be integrated exactly in terms of the Jacobi elliptic functions [17]. Also, here is a place to note that in the theory of nonlinear waves similar equations describe break-up instability phenomena when the high-frequency mode becomes unstable triad with respect to small low-frequency perturbations [12].

4.2. Resonant excitation of the gyro

Let us invite in our consideration the generalized forces describing the damping and forcing of the ring resonator over the axisymmetric form of oscillation. By substituting the terms Qv=2μηv̇; and Qw=2μηẇQcosϖτinto Eq. (19), we can obtain, after the exchange of variables (30), the following set of evolution equations:

da1dT=ηa1iϵn2ω1+p2a¯2a3eiΔωT;da2dT=ηa2iϵn2ω1+p2a¯1a3eiΔωT;da3dT=ηa3+iQeiδT/2iϵn2a1a2eiΔωT,E32

that can be rewritten, after one more exchange of variables; ajT=bjTeiφjT, in terms of real-valued amplitudes and phases:

db1dT=ηb1+ϵn2b2b3sinψω1+p2;db2dT=ηb2+ϵn2b1b3sinψω1+p2;db3dT=ηb3ϵn2b1b2sinψ+QsinδTφ3/2;dT=Δω+2n2ϵb12ω1+p2b22b32b22b32cosψ+ωb1b2QcosδTφ31+p22ω1+p2b1b2b3;dφ3dT=ϵn2b1b2cosψ/b3+QcosδTφ3/2b3,E33

where ψT=φ3Tφ2Tφ1T+ΔωT. These equations are convenient in the study of stationary damped forced motions performed in the system (32). There are two subsets of such motions; the first can be written as

b1=b2=0;b3=K1+p2ηωϵn2,E34

while the second would be

b1=b2=2ηω1+p2K1ϵn2;b3=ηω1+p2ϵn2.E35

Here, we have introduced a new notation: K=Qn2ϵ/41+p2η2ω2for the characteristic number by analogy with Reynolds numbers in the hydrodynamics [18, 19].

Now, we may return to our preliminary study over the spring pendulum to be acquainted in similarities and analogies. Indeed, the first stationary solution (34) appears as stable only within the range 0KK. Here, K=1denotes the critical value of the control parameter governing the bifurcations. This stationary solution coincides exactly with the solution if neglecting all the nonlinear terms in the system (19). Although near the point Kthis stationary solution loses its stability to give its place for the new stationary state (35). The new stationary state would be stable at K>1. It is not now surprising that as if the bifurcation parameter Kgrows even further, then the energy would be pumped only into the low-frequency bending modes while the amplitude of the high-frequency axisymmetric radial mode remains to be fixed, that is, b3=ηω1+p2/ϵn2. Probably, this mechanism represents an effective way in the problem of installation of the solid-state wave gyro.

When returning to the old notation, the solution of our problem can be expressed as it follows:

vτφ=4pb1cosψτcosω+μδ+Δω2τ+Oμ2τ;wτφ=4b1sinψτcosω+μδ+Δω2τ2b3sinϖτ+Oμ2τ,E36

where

ψτ=4p1+p20τΩς

is the precession rate. Typical patterns of the wave precession are shown in Figure 8.

Figure 8.

The wave precession. Modal number n=7. Solid lines refer to flexural mode, while the dashes to the circumferential one (η=0.1,μ=ϵ≈0.01,κ=0.45,Ω=0.085, Q = 0.1). The dotted radius indicates the rotation angle Ωτ. The stiffness of the gyro platform κ is tuned to the seventh bending modes.

5. Conclusion

Léon Foucault has made a sensation with his famous giant pendulum experiment in the Panthéon in Paris in 1851. That time, people could perceive the inertia of the pendulum with his naked eye to have the first direct proof that the Earth is turning anti-clockwise. Moreover, this experiment has demonstrated that a pendulum could be used as a vibratory gyro. However, a flaw of the Foucault pendulum is the dependence on the latitude. A device was needed that was unaffected by the latitude. In 1852, Foucault proposed an angular sensor based on the immobility of the axis of a rotating mass. Unlike the pendulum, this sensor could detect a fixed direction in Newton’s absolute space. Recall that a torque directing the axis of rotation of the rotor of the gyro is known as the precession. For the precision, however, this gyro required an exquisite construction, since the gyro’s rotor should be balanced as well as possible.

In 1890, Bryan investigated the nature of the beats which may be heard when a vibrating cylinder or other axisymmetric thin-walled shell of revolution has been involved in a rotatory motion about its axis [7]. Perhaps, this study had been inspired as a reply to A. E. H. Love, famous for his work on the mathematical theory of elasticity. In the paper [20], Love supposed that unless a body was revolving with the angular velocity comparable with the frequencies of the vibrations, the latter would not be almost affected. The only important effect of rotation would be an extension because of the centrifugal force. Bryan has corrected that in the axisymmetric shells at high-frequency vibrations phenomena of beats may be observed, which appears as the most noticeable effects of the rotation. It is probable that Bryan had no idea how to utilize his finding. Though those times people were in controversy on how far theory of thin-walled shells elaborated by Lord Rayleigh is capable of practice, is this perhaps just a sort of one more abstract theory? Lo and behold, initially conceived in 1890 through the observation of beats from a ringing wine glass, the concept was lost until uncovered in 1965. Due to technologies of the second half of the twentieth century, perhaps in times of Vietnam War, this physical effect has inspired a concept of a solid-state wave gyro [6]. Because there are no typical mechanical parts, these wave sensors have a lot of advantages for long-term space missions. However, to maintain the functionality as well as the sensitivity of a conventional wave gyro in practice, the driving of standing waves requires somewhat sophisticated feedback control. This chapter demonstrates that when both the primary resonant pumping over the axisymmetric mode of oscillations and advantages of the principal parametric resonance are combined, such a gyro can operate without any feedback at the expense of the natural nonlinearity of the resonator in a post-bifurcation regime.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Svetlana Pavlovna Nikitenkova and Dmitry Anatolyevich Kovriguine (April 18th 2018). Bifurcation Makes a Wave Resonant Solid-State Gyro Be Stable, Vibration Analysis and Control in Mechanical Structures and Wind Energy Conversion Systems, Francisco Beltran-Carbajal, IntechOpen, DOI: 10.5772/intechopen.71726. Available from:

chapter statistics

278total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Wind-Induced Vibrations to Tall Buildings and Wind Turbines

By Abdel Rahman Elbakheit

Related Book

First chapter

Adaptive Tuned Vibration Absorbers: Design Principles, Concepts and Physical Implementation

By Philip Bonello

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us