SciELO - Scientific Electronic Library Online

vol.21 issue1The Lineal Prediction Applied to Distributed Speech Recognition in IP Nets author indexsubject indexsearch form
Home Pagealphabetic serial listing  

Services on Demand




Related links

  • Have no similar articlesSimilars in SciELO


Computación y Sistemas

Print version ISSN 1405-5546

Comp. y Sist. vol.21 n.1 México Jan./Mar. 2017 


Digital Rate MEMS Vibratory Gyroscope Modeling, Tuning and Simulation Results

Valerii V. Chikovani1  * 

Hanna V. Tsiruk1 

1 National Aviation University, Institute of Air Navigation, Kyiv, Ukraine.


This paper is dealing with modeling, tuning and simulation of standing wave control and information processing algorithm for rate MEMS ring-type vibratory gyroscope. Standing wave control algorithm is presented in the form allowing transition from rate to rate integrating and to differential modes of MEMS gyro operation by simple switching command. It enables to implement triple-mode MEMS gyro. Block diagrams of the standing wave control algorithm and resonance frequency tracking subsystem are presented. Drive and compensation signal phases tuning criteria are described and implemented. Detailed explanations to the main operations which are used in information processing algorithm are given. Also, simulation results in case of measurement of a constant and variable angle rates are presented.

Keywords: MEMS gyro; standing wave; control algorithm; model

1 Introduction

MЕМS vibratory gyroscopes as a part of Coriolis vibratory gyroscopes (CVG) are nowadays the most applicable in the modern industry and are intensively developing in the laboratories all over the world. The range of their applications is spreading from automotive to medicine including their traditional applications in aerospace. MEMS CVG technology compared to non-MEMS ones is attractive due to its capability to mass production and as a consequence to low cost.

Digital vibratory gyroscope including MEMS consists of three components: micromechanical component - sensing element, electronic component - dedicated electronics and software component - program implementing information processing and standing wave control algorithms. MEMS gyros differ from others by their sensing element that is manufactured by using surface and bulk micromachining technologies [20, 17, 3, 12], such as film deposition, wet and dry etching, plating and others. MEMS structure design process is very complicated and time consuming, therefore there are some commercial MEMS CADs like Coventor Ware [8], Intellisuite [9] and others which facilitate and accelerate design process and analysis of vibrating microstructure.

This modeling using MEMS CADs can be called modeling at a micro-level. The result of micro-level modeling is a micro vibrating structure with certain geometrical and dynamic parameters. However, vibrating structure is not yet a gyro. In order to predict future gyro output characteristics modeling and simulation on the system level are required based on output parameters obtained at the micro-level modeling. MEMS gyro modeling at the system level consists of two sublevels: (1) electronic schematic and circuit board designs, (2) standing wave control and information processing algorithm modeling.

Digital gyro electronic hardware consists of analog, mixed and digital subcomponents. Analog subcomponent is usually called buffer that receives and amplifies analog signals from sensing element and filters signals before sending them to sensing element. Main elements of the mixed subcomponent are analog-to-digital (ADC) and digital-to-analog (DAC) converters.

Digital subcomponent is digital signal processing (DSP) system [2] or field programmable gate array (FPGA) with an auxiliary surrounding [10, 22]. There are many CADs which facilitate and accelerate designing and simulation different types of electronic schemes and circuit boards. Some of them are SPICE, ProfiCAD, PCB123, EasyEDA, Ngspice, logisim, KTechLab and others. Many of them can be accessed in internet.

In order to model and simulate standing wave control and information processing algorithm mathworks Simulink is mainly used [11, 19].

This paper is dealing with modeling, tuning and simulation of a standing wave control and information processing algorithm to show in detail how can be used the output parameters of micro-level modeling to obtain the output parameters of the future gyro. The considered here control and information processing algorithm will take into account low Q-factor of vibrating structure of preferably ring-type (such as micro fabricated ring, hemisphere, cylinder, and other bodies of rotation). Though considered here algorithms can control a standing wave of any MEMS or non-MEMS vibrating structures. As well-known there are two mode of CVG operation - rate and rate integrating ones [2, 10, 11, 19]. Investigations on differential mode of CVG operation have recently been appeared in the publications [4, 5, 6]. Differential mode of operation can effectively suppress external disturbances. All three modes of operations can be combined in single gyro to implement triple mode gyro [7] in contrast to dual mode gyro designed in [14].

In order to implement triple or dual mode gyro the algorithm should be represented in the form that allows one to switch from one mode to another without control program reloading. Most MEMS gyros operate in the rate mode. Rate mode algorithm is a simplest of three abovementioned, and the rate-integrating one is the most complex. In this paper the rate mode algorithm will be represented so that all three modes can be transformed from one to other by simple command of switching. Doing so the rate mode algorithm will obviously take more time for one iteration, than minimum possible, but in this case one can take advantage of switching to another mode swiftly without program reloading.

This paper presents operation principle of rate MEMS gyro on the example of ring resonator, detail description of the standing wave control algorithm operating in the rate mode, explanation of resonant frequency tracking principle built on the basis of digital phase lock loop (PLL), tuning of algorithm parameters such as drive and compensation signal’ phases and scale factor calibration bias dependence on resonator’s manufacturing imperfections.

2 MEMS CVG Principle of Operation

Let’s consider rate MEMS CVG operation principle using ring resonator. The primary standing wave is excited on the second resonant mode of the ring oscillations by applying periodical voltage U=Acos (ωrt) with frequency ωr close to the ring resonant frequency on drive electrode. Standing wave on the second resonant mode of the ring is characterized by four antinodes and four nodes of oscillations located along circumferential coordinate of the ring through the equal angles of 45 deg as depicted in the figure 1.

Fig. 1 Standing wave in a ring type MEMS CVG operation principle 

The solid and dashed ellipses in the figure 1 present the ring deformations in the first and second half of vibration period. When the resonator is rotating about its axis of symmetry (perpendicular to the ring’s plane) with angular rate , Coriolis forces F1, F2, F3 and F4 arise, which excite secondary wave called Coriolis wave in the direction of resultant Coriolis force Fc directed under angle 45 deg relative to the direction of the drive force. The resultant Coriolis force is determined from the following relationship:

, (1)

where k is called an angular gain factor (Brian coefficient) depending on resonator geometry, V is linear velocity of ring’s mass points during vibration, αv,Ω is an angle between vector V and vector Ω. As can be seen from the figure 1 αv,Ω=π/2.

Linear velocity V of the ring’s mass points is:

, (2)

where µ1 is transformation coefficient of applied to resonator voltage, U, and responded mechanical deformation.

Substituting (2) into (1) one can obtain:

. (3)

Thus, Coriolis force is proportional to angle rate Ω.

As can be seen from (3) Coriolis force Fc is changing with resonant frequency ωr and its amplitude is proportional to angle rate. This means that sensing element output signal is amplitude modulated one. The resonant frequency is a carrier signal and angle rate Ω is a modulating signal. In order to obtain angle rate demodulation procedure should be applied to signal (3). Mechanical force Fc causes ring’s deformation with amplitude proportional to this force. Sense electrode, located under the angle of 45 deg to drive one, transforms mechanical amplitude into electric signal Uc= µ2Fc, where µ2 is transformation coefficient of mechanical deformation into responded voltage Uc. Then, electric signal Uc is demodulated to obtain signal proportional to angle rate:

, (4)

where SF is a gyro scale factor known as a result of calibration using rotating table.

In order to increase gyro bandwidth that allows one to measure fast changing angle rate force-to-rebalance method of measurement is used [15].

In the force-to-rebalance measurement method the signal Uc is driven to zero with the aid of feedback control system. The signal that compensates to zero for the Coriolis force is proportional to angle rate.

MEMS CVG can also operate as a rate-integrating (whole-angle mode) gyro [15] at which Coriolis force is not compensated for. In a rate-integrating gyro standing wave freely rotates under action of Coriolis force. In this case standing wave angle of rotation is proportional to gyro angle of rotation and coefficient of proportionality is angular gain factor k (Brian coefficient).

Returning to the rate MEMS it should be noted that parameters A and ωr in SF expression (4) should be stabilized to make SF constant. But ωr (resonant frequency) is changing versus temperature during operation and under environment temperature changes, so this parameter should be tracked and measured. Vibration amplitude should be stabilized at prescribed value A0 during angle rate measurements. Besides, quadrature (error) signal that arises due to manufacturing imperfections and Coriolis force should be nulled to implement force-to-rebalance method of measurement. So, four control systems should be implemented to accurately measure angle rate by MEMS gyro.

In order to implement all four above mentioned control circuits, eight diametrically connected electrodes, shown in figure 1, should be used. So, there are only four independent electrodes, two of them are used to apply control signals on sensing element, specifically, to drive primary vibration and to compensate for the Coriolis force, and another two which are used to sense the responded signals.

3 MEMS Gyro Dynamic Equations

In order to model MEMS gyro, first of all sensing element dynamic equations with real errors resulting from manufacturing imperfections should be accepted. There are well-known dynamic equations (5) of two dimensional pendulum [15] that fit our aim. Equation parameters explanations are presented in figure 2.

Fig. 2 Rigidity and damping axes 

In these equations k is Brian coefficient, dxx is the X axis damping coefficient, τ1 is a minimum resonator’s damping time, τ2 is a maximum resonator’s damping time, dxy is a damping cross-coupling coefficient, kxx is a normalized by mass resonator rigidity along the X axis, ω1 and ω2 are maximum and minimum resonant frequencies, kxy is a rigidity cross-coupling coefficient, dyy is the Y axis damping coefficient, kyy is a normalized by mass resonator rigidity along the Y axis, fx, fy are normalized by mass control signals, θω is an angle between minimum frequency axis and standing wave (antinode) axis, θτ is an angle between minimum damping axis and standing wave axis:


Variables x, y, fx and fy in equations (5) correspond to variables Xout, Yout, Xin and Yin presented in figure 2, respectively. In these designations Xout and Yout are sense signals, i.e. displacement of resonator’s mass point from its initial static position expressed in voltages, and Xin and Yin are control signals. System of equations (5) is valid to describe any type of resonator, not only ring type.

In the rate mode of CVG operation standing wave is immovable relative to resonator. They are rotate together with the same angle rate. In this case angle θ between direction of vibration and direction of drive force (Xin axis) is constant and equal to zero. MEMS sensing element as shown in figure 2 can be interpreted as a two-input-two-output plant. The model of the MEMS resonator based on (5) is presented in figure 3 in designations of Simulink blocks.

Fig. 3 MEMS sensing element model based on equations (5)  

Full schematic in the right side of figure 3 (called sensing element) is presented as a subsystem that has four input and two output signals.

Primary standing wave drive signal, fx, Coriolis force compensation signal fy, input angle rate signal, Om, that should be measured and θ (teta in figure 3) signal which is equal to zero for rate mode of MEMS operation.

As discussed above we will present the standing wave control algorithm in such modification that allows one to change the modes of MEMS operation by simple switching without processor’s program reloading. In rate integrating mode of operation standing wave angle θ is changed. This angle θ during simulation should be calculated and sent it to the sensing element model to change the position of standing wave as it is in the real resonator, but in this paper angle θ=0, because rate mode of MEMS gyro operation is only considered.

Upper schematic of figure 3 represents the model of the first equation of system (5) in designations of the Simulink blocks. The second equation of system (5) model has not shown in the figure 3 because it is similar to the first one. The central schematic shows combination of the two equation models with possibility to introduce sensing element parameters like Q-factor, resonant frequency, frequency mismatch, Q-factor mismatch, and angular dispositions of minimums of damping θτ and rigidity θω axes. These parameters mainly determine resonator’s manufacturing imperfections and gyro accuracy.

4 Rate MEMS Digital Control System

MEMS CVG control system block diagram is presented in figure 4. The standing wave is excited at the resonant frequency ωr of the second mode of ring vibrations. Drive wave control is based on two subsystems: phase lock loop (PLL) and automatic gain control (AGC).

Fig. 4 CVG control system block diagram operating in rate mode of operation 

These two subsystems generate drive signal Xin(t) of proper amplitude, frequency and phase to sustain primary vibration in changing environmental conditions. PLL provides resonant frequency tracking when it changes and AGC keeps the primary vibration amplitude at the constant desired amplitude A0.It should be noted that all signals coming from resonator are amplitude modulated ones.

The resonant frequency is a carrier frequency and angle rate Ω is modulating signal. In order to calculate control signals, first of all demodulation should be carried out. In order to apply control signals to the resonator, re-modulation using resonant frequency should be carried out.At gyro rotation there appears at Yout electrode secondary vibration signal caused by Coriolis force together with error signal which is called quadrature, Q. These two signal components are separated by demodulation processes using reference signals generated by PLL.

Sine component amplitude is proportional to angle rate SF*Ω, and cosine component amplitude is a quadrature signal, which is an error caused by resonator manufacturing imperfections.

Then, these two signals are re-modulated, combined and sent to the Yin electrode as a compensation (force-to-rebalance) signal to keep the standing wave in stationary position relative to resonator.

The control system operates as follows. Antinode Xout and node Yout signals after analog-to-digital converters (ADC) are splitting and provided to the block of demodulators where they demodulate with sine and cosine of resonant frequency as the reference signals to obtain four slow changing signal Cx, Sx, Cy, and Sy as shown in figure 4. In essence, Cx, Sx, Cy, and Sy signals can be used as control signals for rate MEMS gyro.

However, to implement so called “universal” algorithm that allows one, as discussed above, to switch from rate to rate-integrating and to differential modes of MEMS gyro operation without program reloading, these four signals Cx, Sx, Cy, and Sy are provided to block of pendulum variables (A0)2=E, Q, θ, φ calculation.

Physical interpretation of pendulum variables is presented in figure 5 [15], where trajectory of resonator point mass during vibration is shown. In the figure 5, A0 is vibration amplitude, Q is quadrature amplitude, θ is a standing wave angle relative to X axis and φ is a vibration phase.

Fig. 5 Pendulum variables 

These pendulum variables can be calculated as follows [16].


For rate MEMS gyro the pendulum variables should be retained at the following values Q=0, E=const, θ=0, φ=0.

These relationships can be implemented with the aid of four, one for each parameter, proportional and integral (PI) controllers. It should be noted that when: θcomm=0 (see figure 4) MEMS gyro operates in the rate mode, when

MEMS gyro operates in the differential mode [4]. When connection between θr and PIθ controller is opened, MEMS gyro operates in rate-integrating mode. Thus, it can be implemented triple mode MEMS gyro [7]. Dual mode CVG has been implemented in hemispherical resonator gyro (HRG) [14].

4.1 Demodulation

In the block of demodulators four identical synchronous demodulators are used to separate envelope from amplitude modulated signals both from Xout and Yout with the aid of two reference signals sinωrt and cosωrt generated by PLL. Each demodulator consists of multiplication unit, low pass filter to suppress high frequency component, arising after multiplication, to separate low frequency envelope and to double it as depicted in figure 6 to obtain the envelop.

Fig. 6 Synchronous demodulator 

For example, let’s take amlitude modulated signal as follows:

, (7)

where Asin(ωt) is an envelope and sinrt) is a carrier. Graph of this function, for fr=2πωr=4 kHz, f=2πω=100 Hz and A=1, is presented in figure 7.

Fig. 7 Amplitude modulated signal 

After multiplication by reference signal sinrt), low pass filtering (LPF) and doubling, in accordance with block diagram shown in figure 6, the following output signal can be obtained:


4.2 Frequency Tracking System

The idea of the resonant frequency tracking is based on the fact that the drive signal (Xin) and responded signal (Xout) in resonator have π/2 phase difference when drive signal frequency is equal to resonant one. When drive frequency defers from resonant one there appears addition ±φ to π/2 as shown in figure 8. This means that Xout signal has quadrature component which amplitude is proportional to φ and its sign indicates whether the drive frequency is more or less of resonant one. The basic component of Xout signal is used to control (stabilize) the vibration amplitude, and quadrature component is used to track the resonant frequency.

Fig. 8 Amplitude and phase frequency characteristics of resonator 

Really, let’s consider the typical resonator normalized amplitude-frequency (AFC) and phase-frequency (PFC) characteristics. From figure 8 follows that when the drive frequency is not equal to resonant frequency the phase of the Xout signal acquires additional phase ±φ.

For example, let the drive signal is not at the resonant frequency Xin=Bcos((ωr+∆ω)t), then the responded signal, in accordance with figure 8, is Xout=A0sinrt-φ).

This means that to the sine component of Xout signal the cosine component is added. Cosine component of the Xout signal is extracted by the block of demodulators (see figure 4) which is proportional to the phase difference φ between reference signal cos(ωrt) and Xout signal:


The signal proportional to phase difference φ, which is calculated by the last expression of (6), as an error signal is provided to proportional and integral PIPLL controller which, after multiplication by sample time T0 and amplification Kϕ times changes the output frequency of numerically controlled oscillator (NCO) so that to null the error signal φ as shown in figure 9.

Fig. 9 Block diagram of the resonant frequency tracking system based on PLL 

For rate mode of MEMS operation (φ=0) is equivalent to Cx=0.

4.3 Vibration Energy Stabilization

Block diagram of figure 4 shows stabilization of vibration energy E (first equation of (6)). In the rate mode E=(Sx)2=(Ax)2 is X channel vibration amplitude squared, because Cx, Sy, and Cy are driving to null. Thus, sine component of the X channel signal Sx is extracted by the corresponding demodulator, then squares it and sends to discriminator (see figure 4).

The error signal, difference between real square of vibration amplitude and desired one (A0)2 is provided to the PIE controller, which gives the signal on its output so that to drive the error signal to null, that is [(A0)2 -- (Ax)2] → 0.

Thus, vibration energy stabilization results in vibration amplitude stabilization at desired value A0. The controller output is DC signal which should be re-modulated by multiplying by reference signal cos(ωrt) to send it to DAC and then, in analog form (voltage), to driving electrode.

4.4 Standing Wave Angle Control

When θcomm=0, PIθ controller (see figure 4) generates such an output signal that drives to null the difference (θr - θcomm) → 0, where θr is an actual standing wave angle and θcomm is a commanded value of standing wave angle. When θr=0, drive Xin and compensation Yin signals are:


where Qout is an output signal of the PIQ controller, and Zy is a signal proportional to angle rate measured. Thus, X is a drive channel and Y is a measurement channel.

When: θcomm=θ0=const≠0,πk/4, k=1,2,…, i.e. standing wave angular position does not coincide with any of eight electrodes, there arise two X and Y measurement channels with opposite angle rates. The biases and scale factors of these two channels are dependent on angle θ0. One can choose angle θ0 such that biases of X and Y channels will be equal to each other, then the difference of the two channels increases angle rate and compensates for the bias in difference channel.

There is another option of choosing angle θ0 such that scale factors of the two channels are equal to each other. In this case external disturbances acting on the gyro can be effectively compensated for in difference channel [4, 5]. There are also other advantages of differential mode of operation [6].

When connection between θr signal and PIθ controller is opened MEMS gyro operates in rate-integrating (whole angle) mode. In this case standing wave freely rotates under action of Coriolis force and θr is output signal proportional to angle of gyro rotation.

5. Drive and Compensation Signal Phases Tuning

In the gyro model like in real gyro it is important to tune drive and compensation signal phases. Tuning criteria can be different. In this work to tune drive phase, criterion of minimum of drive signal amplitude providing desired vibration amplitude A0 is chosen.

Figure 10 shows dependence of drive signal amplitude on its phase. From this figure one can see that this minimum is not sharp and permits variation in the vicinity of minimum point without significant change in amplitude. Thus, minimum amplitude is reached in the phase range [1-1.1] rad.

Fig. 10 Drive signal phase tuning 

In spite of that model output is dimensionless and in real digital MEMS gyro output signals are represented in digital code, it is advisable to normalize output signal of ADC such that output code would be numerically equal to the voltage on its input. It will be very convenient to visually analyze such digital signals. Therefore, output signals of the gyro model in this paper will be represented in volts.

In order to tune phase of compensation signal, minimum of module of angle rate in percent which is visible in the quadrature signal is chosen as a tuning criterion. Figure 11 shows the dependence of angle rate portion in percent which is visible in quadrature signal versus compensation signal phase.

Fig. 11 Compensation signal phase tuning 

As can be seen from this figure minimum is reached at the phase value of 1.05 rad which is in the range of minimum of drive signal phase. It is interesting to note, that in real gyro drive and compensation signal phases are also close to each other. When compensation signal phase is more than 1.2 rad or less than 0.2 rad, then feedback control subsystem providing compensation for the Coriolis and quadrature signals becomes unstable. When these tunings are not properly accomplished, then the higher portion of input angle rate will be visible in the drive and quadrature signals that results in accuracy reduction or even instability.

6 Simulation Results

In order to start simulation let’s introduce a noise in the model to make the model close to real gyro. The noise can be introduced on the output of ADC with such standard deviation that the output angle rate signal noise would be close to known MEMS gyro, for example, one of the gyro from the block of sensors ADIS16367 (Analog Devices, USA).

Gyro model operates at sampling frequency 50 kHz and output LPF has cutting frequency 300 Hz. Besides, gyro model scale factor calibration should be made before the beginning of simulation, like it occurs in real gyro. After that, the output signal proportional to angle rate measured should be divided by scale factor to obtain angle rate in deg/s.

Usually, scale factor (SF) calibration is conducted for some different angle rates from the gyro measurement range, and then, these scale factor values are averaged to reduce errors from SF nonlinearity. We simplify this procedure and conduct gyro model SF calibration at one value of angle rate equal to 100 deg/s. For this purpose, in accordance with international IEEE standard [15], it is necessary to rotate the gyro clock-wise and counter-clock-wise with equal angle rates of Ωr=±100 deg/s, then, mean values of measurements Np and Nm are calculated for each direction of rotation, and then, SF is calculated by the formula:

. (11)

Figure 12 shows noisy output signal of the gyro model under clock-wise and counter-clock-wise rotation with Ωr=±100 deg/s with indication of mean values Np and Nm. Calculation shows that SF=0.03241 V/(deg/s).

Fig. 12 Scale factor and bias measurement during rotation 

Because gyro model (and real digital gyro) output signal is represented in digital code numerically equal to voltage on the input of ADC, as we discussed above, SF dimension can also be represented in 1/(deg/s). Using the same measurements, Np and Nm, gyro bias, B, during rotation can be calculated as follows:

. (12)

After scale factor determination, rate MEMS gyro model output signal will be presented in deg/s.

It should be noted that rate MEMS gyro bias depends on resonator parameters imperfections, and SF linearly depends on resonant frequency and vibration amplitude (4). Bias and SF have been obtained at the following resonator parameters:

Figure 13 shows dependence of bias on frequency mismatch expressed in percent (∆f/f)*100%. Negative percent means that axes of maximum and minimum resonant frequency exchanged their positions.

Fig. 13 Bias versus frequency mismatch 

As can be seen frequency mismatch significantly affects bias and this dependence is close to parabola. This graph says that bias is very sensitive to frequency mismatch. This sensitivity is about 130 deg/s/%. So, frequency mismatch should be made close to zero. In order to make frequency mismatch close to zero in real MEMS gyro electrostatic stiffness tuning [21] or laser trimming [22] can be used.

Figure 14 shows dependence of bias on Q-factor mismatch expressed in percent, (∆Q/Q)*100%.

Fig. 14 Bias versus Q-factor mismatch, ∆Q/Q 

In contrast to figure 13 Q-factor mismatch influences on bias non-symmetrically and has much lower sensitivity of about 0.03 deg/s/%, than that of frequency mismatch.

Figure 15 shows bias dependence on θω angle.

Fig. 15 Bias versus θω angle 

This dependence is close to harmonic with two periods on full revolution.

Closer to harmonic one is the dependence of bias on θτ angle shown in figure 16.

Fig. 16 Bias versus θτ angle 

This figure shows that bias sensitivity to θτ angle is much lower with the maximum value of about 1.4×10-3 (deg/s)/deg in contrast to bias sensitivity to θω with the maximum value of about 3×10-2 (deg/s)/deg.

From the data presented in figures 15 and 16 follows that by forced rotation of a standing wave, through the angle 90 deg, change of a gyro bias, caused by frequency and Q-factor mismatches, can be compensated for.

Figure 17 shows gyro model output signal when input angle rate changes its frequency from DC to 200 Hz. This figure shows bandwidth of the gyro model on the level of 70 Hz, and non-uniformity in band pass of no more than 5%. SF does not depend on the considered here parameters. Gyro model and real gyro SF linearly depends on amplitude of vibration and resonant frequency. SF sensitivity to resonant frequency is much less than that of vibration amplitude.

Fig. 17 Model gyro bandwidth determination 

7 Conclusion

This paper presents MEMS gyro model based on “universal” algorithm enables implementing any of three modes of MEMS operation - rate, rate-integrating and differential ones. Also, this gyro model gives opportunity to implement new triple mode MEMS gyro.

Using this MEMS gyro model tuning procedure has been shown. Dependences of bias on orientation angle θω of rigidity and angle θτ of damping axes have been obtained. Based on these dependences maximum value of bias sensitivities to such resonator imperfections as frequency, Q-factor mismatches, θω and θτ angles have been determined. It was found out that gyro bias can be compensated for by forced rotation of standing wave through the angle 90 deg.

This rate MEMS gyro model can be used to investigate influence of shock and vibration on the output signal in order to estimate gyro sensitivity to them.

Moreover, presented in this paper gyro model can help developer to implement rate-integrating (whole angle) and differential modes of MEMS operation, investigate their sensitivity to different types of internal and external disturbances.

In this paper PI controllers are used to provide necessary bandwidths for the MEMS gyro control subsystems. In order to improve MEMS gyro control in harsh environment robust PID tuning [13], adaptive PID, fuzzy logic, sliding mode [1] and their combinations can be used.


1. Aguilar-Mejía, O., Tapia-Olvera, R., Rivas-Cambero, I., & Minor-Popocatl, H. (2016). Diseño de un controlador de velocidad adaptativo para un MSIP utilizando inteligencia artificial. Computación y Sistemas, Vol. 20, No. 1, pp. 41-54, DOI: 10.13053/CyS-20-1-2361. [ Links ]

2. Bowles, S. R. (2015). Design and implementation of a control scheme for a MEMS rate integrating gyroscope. PhD thesis, Newcastle University, Oct., pp. 148. [ Links ]

3. Bustillo, J. M., Howe, R. T., & Muller, R. S. (1998). Surface micromachining for microelectromechanical systems. Proceedings of the IEEE, Vol. 86, No. 8, pp. 1552-1574. [ Links ]

4. Chikovani, V. V., & Suschenko, O. A. (2014). Differential mode of operation for ring-like resonator CVG. IEEE International Conference on Electronics and Nanotechnology (ELNANO), NTUU-KPI, Kyiv, Ukraine, pp. 451-455. [ Links ]

5. Chikovani, V. V., & Tsiruk, H. V. (2014). Shock resistance of differential type ring-like resonator vibratory gyroscope. Electronics and control systems, Vol. 4, No. 42, pp. 79-83. [ Links ]

6. Chikovani, V. V., Suschenko, O. H., & Tsiruk, H. V. (2016). Redundant information processing techniques comparison for differential vibratory gyroscope. Eastern-European Journal of enterprise technologies 4/7 (82), pp. 45-52, DOI: 10.15587/1729-4061.2016.75206. [ Links ]

7. Chikovani, V. V., & Tsiruk, H. V. (2015). Differential mode of operation for multimode vibratory gyroscope. IEEE International Conference on actual problems of unmanned aerial vehicles development (APUAVD), NAU, Kyiv, Ukraine, pp. 87-90. [ Links ]

8. CoventorWare (2016). http://www.coventor.comLinks ]

9. IntelliSuite (2016). Woburn, MA: IntelliSense Corp. http://www.intellisense.comLinks ]

10. Gregory, J., Cho, J., & Najafi, K. (2012). Characterization and control of a high-Q MEMS inertial sensor using low-cost hardware. IEEE Conf. Proc., Myrtle Beach, S.C., pp. 239-247, DOI: 10.1109/plans.2012.6236886. [ Links ]

11. Hui, L., & Quan, H. (2013). Research of ring MEMS rate integrating gyroscopes. World Academy of Science, Engineering and Technology, Vol. 7. [ Links ]

12. Kovacs, G. T. A., Maluf, N. I., & Petersen, K. E. (1998). Bulk micromachining of silicon. Proceedings of the IEEE, Vol. 86, No. 8, pp. 1536-1551. [ Links ]

13. Lagunas-Jiménez, J. R., Moo-Yam, V., & Ortíz-Moctezuma, B. (2014). Two-degrees-of-freedom robust PID controllers tuning via a multiobjective genetic algorithm. Computación y Sistemas, Vol. 18, No. 2, pp. 259-273, DOI: 10.13053/CyS-18-2-2014-031. [ Links ]

14. Lynch, D. D., & Matthews, A. (1996) Dual mode hemispherical resonator gyro operating characteristics. 3-rd S. Petersburg int. conf. on integrated navigation systems, S. Petersburg, Russia, pp. 37-44. [ Links ]

15. Lynch, D. D. (2004). Coriolis vibratory gyroscope. IEEE standard specification format guide and test procedure for Coriolis vibratory gyros. IEEE std.1431TM, annex B, pp. 56-66. [ Links ]

16. Lynch, D. D. (1995). Vibratory gyro analysis by the method of averaging. Proc. of the 2nd S.Petersburg International Conference on gyroscopic technology and navigation, St. Petersburg, Russia, part 1, pp. 26-34. [ Links ]

17. Quabili, A., & Dang, T. (2009). Bulk and Surface Micromachining. Group Activity #3 Report, 10 May, P. 5. [ Links ]

18. Sonmezoglu, S., Alper, S. E., & Akin, T. (2014). An automatically mode-matched MEMS gyroscope with wide and tunable bandwidth. Journal of microelectromechanical systems, Vol. 23, No. 2, pp. 284-297. [ Links ]

19. Su, Z., Liu, N., Li, Q., Fu, M., Liu, H., & Fan, J. (2014). Research on the signal process of a bell-shaped vibratory angular rate gyro. Sensors, Vol. 14, pp. 5254-5277, DOI: 10.3390/s140305254. [ Links ]

20. Tai-Ran, H. S. U. (2008). MEMS and Microsystems: Design, Manufacture, and Nanoscale Engineering. John Wiley & Sons Inc., 2nd ed., pp. 527. [ Links ]

21. Xia, D., Yu, C., & Kong, L. (2014). The development of micromachined gyroscope structure and circuitry technology. Sensors, Vol. 14, pp. 1394-1473, DOI: 10.3390/s140101394. [ Links ]

22. Zhao, X., Zhou, B., Zhang, R., & al. (2011). Research on laser trimming of silicon MEMS vibratory gyroscopes. Vibratory gyroscopes. Integrated ferroelectrics Jan., Vol. 129, No. 1, pp. 37-44. [ Links ]

Received: December 08, 2016; Accepted: January 20, 2017

* Corresponding author: Valerii V. Chikovani, e-mail:;

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License