PSPICE Hybrid Modeling and Simulation of Capacitive Micro-Gyroscopes.

With an aim to reduce the cost of prototype development, this paper establishes a PSPICE hybrid model for the simulation of capacitive microelectromechanical systems (MEMS) gyroscopes. This is achieved by modeling gyroscopes in different modules, then connecting them in accordance with the corresponding principle diagram. Systematic simulations of this model are implemented along with a consideration of details of MEMS gyroscopes, including a capacitance model without approximation, mechanical thermal noise, and the effect of ambient temperature. The temperature compensation scheme and optimization of interface circuits are achieved based on the hybrid closed-loop simulation of MEMS gyroscopes. The simulation results show that the final output voltage is proportional to the angular rate input, which verifies the validity of this model.

1. Introduction

In recent years, due to a rapid development of microelectromechanical systems (MEMS) technology, MEMS gyroscopes have become an indispensable portion of the inertial navigation, military, and consumer electronics market for the advantages of their small size, light weight, low cost, and high reliability [1]. The MEMS gyroscope is a device for angular rate detection, and its basic principle is based on the Coriolis Effect [1]. The proof mass of a gyroscope oscillates along the drive axis with a stable frequency and stationary amplitude; meanwhile, this gyroscope is subjected to an angular rate Ω in the input axis. Then, it will generate a vibration of the same frequency along the sense axis due to the Coriolis force. The angular rate information is provided by the sense axis, which specifically is proportional to the direct current (DC) output of the sense axis. It is necessary to note that the three axes mentioned above are orthogonal to each other.

MEMS gyroscopes have achieved significant commercial success and have benefited from better and cheaper manufacturing technology. Although the final cost of gyroscopes from mass commercial production is relatively low, it should be noted that prototype design and manufacture and repeat test experiments are costly and consume time. Therefore, building an accurate and complete model is of great significance for reducing the cost of prototype development and guiding the design of peripheral interface electronics. The most commonly used model is the dynamic model, and a good deal of literature [2,3,4,5,6] reveals that the drive mode and sense mode of MEMS gyroscopes can be considered as a second-order Spring-Mass-Damping system. The dynamic model converts the driving force to displacement by the transfer function. The advantage of this model is its intuitive and convenient representation of typical mechanical parameters of MEMS gyroscopes, such as the elastic coefficient of cantilevers and the damping coefficient of the system. Nevertheless, only the displacements of gyroscopes are considered and used to make a system-level simulation, which ignores the actual capacitance changes and the connection to subsequent interface electronics. Considering that micro-Gyroscopes also have characteristic parameters, such as resonant frequency and the quality factor of general resonators, the pure circuit model is also commonly used. In [7], a gyroscope is simplified to a series Resistance-Inductance-Capacitance (RLC) circuit model, where the resonant frequency and quality factor are determined by the size of resistance, inductance, and capacitance in the circuit. This model is easy to implement and it is convenient to combine the RLC model with subsequent interface electronics. However, this model cannot demonstrate some of the mechanical characteristics of MEMS gyroscopes, such as the elastic coupling force and quadrature error. Another novel model is the MIMO T-S fuzzy model, which is built based on nonlinear characteristics of MEMS gyroscopes [8]. The model can strengthen the robustness of the control system with an adaptive sliding controller when external disturbances exist. Nonetheless, this model is non-intuitive and it is challenging to implement the connection between it and subsequent interface electronics.

In this paper, a hybrid model for MEMS gyroscopes is firstly established based on a PSPICE Simulation Model, including two crucial aspects for modeling a micro-gyroscope device, namely different modules in gyroscopes and a principle diagram for them. Then, systematic simulations of the model itself are conducted. At last, hybrid simulations of the model and interface circuits are implemented to achieve closed-loop control of MEMS gyroscopes and optimization of the circuits. The contribution of the present work lies in three aspects.

A detailed model of MEMS gyroscopes is developed. This model considers mechanical thermal noise equivalent disturbance as well as the capacitance model without approximation. The model also considers the effect of temperature on the displacements of both modes.

All output ports of the model are capacitive interfaces, which can be directly connected to conditioning circuits. Therefore, a closed-loop simulation of the model and interface circuits is achieved. A calibration scheme for temperature is developed based on a different Zero Rate Output (ZRO) at different temperatures.

Based on the simulation results, optimization designs of interface circuits are achieved, including the value of the demodulation signal phase and the circuit gain k in the closed-loop detection circuit.

The remainder of this paper is organized as follows. In Section 2, PSPICE models of different modules in gyroscopes are firstly proposed, including a sensitive structure model, a Coriolis force and elastic coupling force model, a mechanical thermal noise equivalent disturbance model, a mechanical model, the temperature models of gyroscope parameters, and a differential capacitance model. Then, the PSPICE device model of capacitive MEMS gyroscopes is established based on these modules together with its principle diagram. Lastly, the interface circuits of the closed-loop simulation are also discussed. In Section 3, a systematic simulation of this model and a closed-loop simulation of MEMS gyroscopes are conducted with an analysis of the effect of ambient temperature and optimization of interface circuits, and simulation results verify the validity of this model. In Section 4, the corresponding discussion is given.

2. Methods

2.1. PSPICE Models of Different Modules in MEMS Gyroscopes

A complete MEMS gyroscope is a system of coupled mechanical and electrical components. Its capacitance interface is used to achieve the electromechanical coupling process. Figure 1 shows the capacitance interface of the modeled gyroscope.

Figure 1

The capacitance interface of the modeled gyroscope.

In Figure 1, the x-axis is the direction of drive mode. The y-axis is the direction of sense mode, and the z-axis is the input direction of the angular rate. The blue movable electrode is the structure (STRC) electrode, and it is biased with a DC voltage U0. In terms of drive mode, the four orange plates are exactly the same and fixed. The upper two plates are actuation electrodes D1+ and D1−. Asinω and −Asinω are differential alternating current (AC) excitation voltages applied to the two electrodes for generating an electrostatic driving force F on the movable plate. The lower two plates D2+ and D2− are detection electrodes of drive mode, which detect the drive-axis displacement of the STRC electrode. The same structure exists in sense mode. The plates S1− and S1+ are detection electrodes of the sense mode, which detect the sense-axis displacement of the STRC electrode. The plates S2− and S2+ are force feedback electrodes of sense mode and are used to compensate for the Coriolis force induced by the input angular rate.

There is some software that can be used in the simulation of MEMS gyroscopes, such as COMSOL, SIMULINK, and PSPICE. COMSOL mainly emulates the mechanical characteristics of MEMS gyroscopes by finite element analysis, so it is incapable of simulating a complete gyroscope system including the interface circuits. Large quantities of simulations on MEMS gyroscopes are conducted in SIMULINK. In [9], a behavioral model of MEMS gyroscopes constructed in SIMULINK is presented. This model simply describes the mechanical model of gyroscopes and behavioral models of interface circuits. However, SIMULINK software is inadequate to simulate actual circuit electronics and changes of environmental parameters, such as ambient temperature. Because of its good convergence, PSPICE software is suitable for system-level and circuit-level simulation. PSPICE has rapid and accurate simulation ability, so it has been successfully used in a wide variety of linear and nonlinear electrical circuit simulations [10]. PSPICE can simulate both the electronic and non-electronic components, which is consistent with the purpose of this paper, namely the modeling of MEMS gyroscopes and the simulation of interface circuits. Other important reasons to choose PSPICE as the simulation tool are that PSPICE uses more accurate and realistic analog electronic component models and the parameters of models are variable in the simulation process. Figure 2 shows the principle diagram of capacitive MEMS gyroscopes, which is based on the diagram of MEMS sensors in [11]. Different modules of MEMS gyroscopes in this principle diagram are established by a PSPICE Simulation Model, and these modules are discussed separately in the following sections. This diagram illustrates that the inputs of gyroscopes include the AC excitation voltage for drive mode, the input angular rate for sense mode, and the mechanical thermal noise existing in both modes. The outputs of gyroscopes are the changes of differential capacitance, which are caused by the vibration displacements of both modes.

Figure 2

The principle diagram of capacitive microelectromechanical systems (MEMS) gyroscopes.

2.1.1. The Sensitive Structure

For drive mode, the electrostatic driving force is induced by the differential AC excitation voltages acting on the sensitive structure. Specifically, Asinω and −Asinω are applied to the actuation electrodes D1+ and D1−. Additionally, the STRC electrode is biased with a DC voltage U0. This electrostatic force F is normal and it is expressed as a gradient of the electrical potential energy in the capacitor.

The energy of a parallel plate capacitor is where C is the capacitance value and V is the voltage applied to the capacitor. The electrostatic driving force F+ between the movable plate and the right actuation plate can be calculated as follows: where d is the initial spacing of the two plates, x is the drive-axis displacement of the STRC electrode, C0 is initial capacitance between the movable plate and actuation plates of drive mode and its value is , ε is the dielectric constant, and a is the plate area.

The same derivation is also applied to the left actuation plate and the movable plate. Additonally, the total electrostatic driving force F for the movable plate is shown below:

The above formula shows that F is related to Asinω and x. In case of x << d, F can be simplified as follows:

Formula (4) indicates that F is proportional to the excitation voltage Asinω when x << . Additionally, the gain factor k is . With the sensitive structure, the AC excitation voltage can be converted into the electrostatic driving force of drive mode. The scale factor k will be used in the closed-loop simulation of MEMS gyroscopes.

2.1.2. The Coriolis Force and Elastic Coupling Force Model

This model describes the mechanism for the generation of Coriolis force and elastic coupling force. The Coriolis force is induced by the Coriolis Effect. It is related to the angular rate input and the displacement velocity of drive mode. The Coriolis force can be calculated according to the specific formula, where is the angular rate of rotation and is along the z-axis, V is the displacement velocity of the STRC electrode in drive mode, and m is the effective mass of drive mode, which also happens to be the proof mass in our gyroscope model. The direction of the Coriolis force is along the y-axis. It should be noted that when the input angular rate is relatively large, the Coriolis force generated by the movement of sense mode will also affect the drive mode. The corresponding Coriolis force is where V is the displacement velocity of the STRC electrode in sense mode. This force will act upon the drive mode of MEMS gyroscopes.

MEMS gyroscopes may be affected by fabrication imperfections which can lead to a coupling error between the drive mode and the sense mode. The most common error is a quadrature error due to spring imbalance, which causes an additional elastic coupling force. This force is orthogonal to the Coriolis force in the phase [12], and it is proportional to coupling stiffness and drive-axis displacement. The size of a coupling force can be calculated as where k is the coupling stiffness, which is used to characterize the effect of spring imbalance. The direction of an elastic coupling force is also along the y-axis.

2.1.3. The Mechanical Thermal Noise Equivalent Disturbance

For the mechanical thermal noise existing in a MEMS gyroscope’s structure, [13] indicates that the effect of the noise can be equivalent to a constant disturbance force upon both modes. The size of the equivalent force is where is Boltzmann’s constant, and its value is 1.38 × 10−23 J/K, is degrees Kelvin, and B is the bandwidth of the gyroscope. For a mode-split gyroscope, its bandwidth is approximately equal to 0.54 times the frequency difference when operating in open-loop detection mode [14]. Additionally, this bandwidth is valid for the simulation of a PSPICE device model of capacitive MEMS gyroscopes. c is the damping coefficient of the gyroscope structure, and it is different for the two modes. The corresponding values of the damping coefficient are shown below: where Q and ω are the quality factor and resonant angle frequency of drive mode, respectively, Q and ω are the quality factor and resonant angle frequency of sense mode, respectively, and m and m are the effective drive mass and effective sense mass, respectively. It is noteworthy that the quality factor and resonant frequency of gyroscopes are greatly affected by ambient temperature. The effect of temperature on parameters will be discussed in the following temperature model.

2.1.4. The Mechanical Model of MEMS Gyroscopes

A mechanical model of MEMS gyroscopes can convert the resultant forces into displacements for both modes. Additionally, a gyroscope can be considered as a second-order Spring-Mass-Damping system. Its characteristics equations are where k and k are the stiffness coefficients of the drive mode and sense mode, respectively, F and F are the mechanical thermal noise equivalent disturbance force of the drive mode and sense mode, respectively, and F and are the driving forces of the two modes, respectively. The transfer functions of the characteristic equations are where s is the variable symbol jω of the transfer function. The transfer functions convert the corresponding resultant forces into displacements for both modes when mechanical model parameters of gyroscopes are known. Additionally, the mechanical model is also related to temperature.

2.1.5. The Temperature Models of Gyroscope Parameters

As discussed earlier, the quality factor and resonant frequency of gyroscopes are greatly impacted by ambient temperature. Therefore, the temperature will affect the displacements of both modes. The authors in [15] indicate that temperature and resonant frequency approximate a linear relationship. A corresponding experiment was carried out to obtain the linear factor based on the gyroscopes in the authors’ research group. The varying resonant frequency measurement is based on the dynamic signal analyzer. At different temperatures, the phase-locked loop controls the drive mode to operate at different working frequencies. Because the Q value of drive mode is very large, the working frequency is approximately equal to the resonant frequency. Therefore, the resonant frequencies of drive mode are obtained by measuring the working frequencies through the “channel” port of the dynamic signal analyzer at different ambient temperatures. For sense mode, the measurement of resonant frequency at different temperatures is achieved through the “source” and “channel” ports of the dynamic signal analyzer. The “source” port outputs a chirp sweep signal to the force feedback electrode. Additionally, there are corresponding response signals with different amplitudes on the detection electrode. These response signals are displayed on the dynamic signal analyzer through the “channel” port. The point with the largest response amplitude corresponds to the resonant frequency of sense mode. Then, the relationship between the resonant frequency and the ambient temperature is shown as follows where and are the resonant frequency of drive mode and sense mode, respectively, and T is degrees centigrade.

In order to obtain a higher quality factor, the microstructure of MEMS gyroscopes is vacuum-packaged. In [16], the quality factor is proportional to the resonant frequency and inversely proportional to the square root of the degrees kelvin. Additionally, the proportional factors of both modes are different. The resonant frequency and quality factor at 27 degrees centigrade will be used to calculate the proportional factor of both modes. Therefore, the Q values of two modes at 27 degrees Celsius need to be measured. The Q value measurement needs the “source” port and the “channel” port of the dynamic signal analyzer. For the drive mode and sense mode, the “source” port outputs a chirp sweep signal to the actuation electrode and the force feedback electrode, respectively. The corresponding response signals are generated on the detection electrodes of both modes and they are displayed on the dynamic signal analyzer through the “channel” port. The Q value is obtained through dividing the resonant frequency by the corresponding −3 dB bandwidth. A corresponding experiment was carried out to obtain the linear factor based on the gyroscopes in the authors’ research group and the formulas are shown as follows:

2.1.6. The Differential Capacitance of MEMS Gyroscopes

In Figure 1, the lower two orange plates and the movable blue electrode form the differential detection capacitance of drive mode, and the capacitance is, respectively:

For sense mode, its differential detection capacitance is formed between the two fixed red plates and the movable blue electrode. The differential capacitance is, respectively: where d is the initial spacing of the capacitance, a is plate area, and y is the sense-axis displacement of the STRC electrode. In our capacitance model, the actual parallel-plate detection mechanism is used. The capacitance change is not simply approximated to a direct proportion of the displacement. This approximation is sometimes used in the simulation of MEMS sensors. The distance of the MEMS capacitor plate is much less than its length and width. Therefore, the edge effect of capacitance can be ignored. The parallel-plate sensing mechanism contributes a nonlinear behavior between the sense capacitance and the sense-axis displacement. This nonlinearity will be eliminated in normal operation because the displacement produced by the Coriolis force is suppressed by the feedback force. In actual gyroscopes, the capacitive nonlinearity is from multiple effects, including microfabrication process errors, parallel plate nonlinearity due to deformation, C-V conversion circuit behavior, and the quadrature error existing in gyroscopes. In our model, these factors are assumed to be ideal because the entire simulation is based on the PSPICE’s own model components and electronics.

2.2. Establishment of Capacitive MEMS Gyroscopes PSPICE Device Model

The parameters of the PSPICE device model of parallel-plate capacitive micro-gyroscopes are shown in Table 1 at a constant temperature of 27 degrees centigrade. All other parameters except are derived from the gyroscopes in the authors’ research group. The data for in Table 1 is not experimental data but simulation data. The authors in [4,17] indicate that the typical quadrature force is much larger than the Coriolis force used for sensing the external angular rate. In the case of other parameters, they are known assuming that = 5 N/m is in accordance with the precondition (assuming the typical angular rate Ω is 0.01°/s). Additionally, = 5 N/m is experimental data according to [4]. Therefore, we take = 5 N/m as our simulation data.

Table 1

Micro-gyroscopes PSPICE model simulation parameters (T = 27 °C).

Parameter	Value
fx	10.157 kHz
fy	10.257 kHz
fd	10.157 kHz
U0	6 V
dx = dy	2 μm
ε	8.85 p
ax = ay	14.256 × 10−9 m2
mx	12.3 μg
Qx	150,000
kyx	5 N/m
my	0.245 μg
Qy	10,000
A	5 V

The PSPICE device model is established according to Figure 2 as shown in Figure 3. The formulas in Figure 3 represent the modules in Figure 2 and they have been discussed in detail in Section 2.1.

Figure 3

The PSPICE device model of capacitive micro-gyroscopes.

This model consists of drive mode and sense mode. In terms of drive mode, the resultant force and mechanical model form a loop to generate a displacement. All of the output ports are differential capacitance interfaces that can be directly connected to a C-V conversion circuit. For the sense mode, the driving force is induced by the Coriolis force and elastic coupling force model. It will be converted into a displacement with the mechanical model. The displacement affects the differential capacitance. The final differential capacitance output is related to the angular rate input. The specific simulation results of a parallel-plate capacitive micro-gyroscopes PSPICE model at a constant temperature of 27 degrees centigrade are shown in Section 3.

The temperature will affect the resonant frequency and quality factor, which changes the gyroscope’s mechanical model. So, the gyroscope will have different motion states at different temperatures. The simulation results at different temperatures with the same excitation voltage are also described in Section 3.

2.3. The PSPICE Closed-Loop Simulation of MEMS Gyroscopes

The temperature will change the motion state of a gyroscope, so it is necessary to implement the closed-loop control of MEMS gyroscopes. For drive mode, the Phase-locked Loop (PLL) method [18] is usually used to control the drive frequency’s stability. The Automatic Gain Control (AGC) method [19] is used to keep the amplitude of the drive-axis displacement constant. For sense mode, the Closed-loop Detection method [20] is used to ensure that the proof mass always vibrates near the equilibrium position. The corresponding block diagram of the closed-loop simulation based on the above three methods and the micro-gyroscopes PSPICE model is shown in Figure 4.

Figure 4

The closed-loop simulation diagram of MEMS gyroscopes.

For drive mode, the control circuits consist of a C-V conversion circuit, a PLL frequency control circuit, and an AGC amplitude control circuit. The gain factors k and k represent the sensitive structure and the C-V conversion circuit, respectively. When the drive mode is stable, the output of the loop filter should be a stabilized value at a fixed temperature, which means that the vibration frequency of displacement reaches stability. It is noteworthy that the stabilized value will change slightly with the ambient temperature due to the change of resonant frequency. Additionally, the output of the amplitude detection module is equal to the reference value when the drive mode is stable, which means that the vibration amplitude of the displacement reaches a constant value. The corresponding simulation results of closed-loop control are shown in Section 3.

For sense mode, the purpose of closed-loop detection is to suppress the Coriolis force completely. This is accomplished by adding differential AC signals to the feedback electrodes of sense mode, namely and in Figure 1. It is noteworthy that the phase of the differential AC signal needs to be opposite to the phase of the displacement velocity of drive mode. A corresponding feedback force is induced to suppress the Coriolis force by adjusting the amplitude of the differential AC signals. The feedback force can be calculated by formula (4), and its value is where C1 is initial capacitance between the movable STRC and each of the four fixed plates of the sense mode, is the AC signal, which has opposite phase of the displacement velocity of drive mode, and is the simplified interface gain of the force feedback structure. In our model, the output of the gain module is the output voltage of MEMS gyroscopes corresponding to different angular rates.

3. Results

3.1. The Simulation Results of the Capacitive MEMS Gyroscopes PSPICE Device Model

The simulation results of the capacitive micro-gyroscopes PSPICE model are described in Figure 5. Figure 5a shows that the drive-axis displacement and sense-axis displacement are stable sine waves with the same frequency as the AC excitation voltage. This means that the Coriolis Effect has been achieved through the established PSPICE model. Figure 5b indicates that the amplitude of the sense-axis displacement changes with the angular rate. It is noteworthy that the sense-axis displacements corresponding to different angular rates will intersect at specific points. At these points, the phases of the sense-axis displacement are some periodic values, which will cause the sense-axis displacement to be zero and the quadrature displacement to be the maximum. Therefore, these points represent the quadrature displacement, which is induced by an elastic coupling force and its amplitude is 4.804 nm. The sense-axis displacement is not proportional to the angular rate input due to the existence of quadrature displacement. Figure 5c shows that the displacements of gyroscopes can be converted into differential capacitance changes. Thus, all output ports of this PSPICE model are capacitance interfaces that can be connected directly to a C-V conversion circuit.

Figure 5

The simulation results of the capacitive micro-gyroscopes PSPICE model: (a) drive-axis displacement and sense-axis displacement, Ω = 30°/s (b) sense-axis displacements at different angular rates (c) the differential detection capacitance of drive mode and sense mode, Ω = 30°/s.

Changes in ambient temperature will cause changes in the resonant frequencies and quality factors of MEMS gyroscopes, which can affect the dynamic characteristics of the gyroscope. Therefore, the same excitation voltage will induce different drive-axis displacement and sense-axis displacement at different ambient temperatures. Figure 6 shows the corresponding simulation results.

Figure 6

Displacements of the capacitive MEMS gyroscopes PSPICE model at different temperatures with the same excitation (): (a) drive-axis displacement (b) sense-axis displacement.

The simulation results show that the drive-axis displacement and sense-axis displacement will have different phases and amplitudes at different temperatures, which indicates that the gyroscope is greatly affected by the ambient temperature. Therefore, subsequent conditioning circuits are necessary to achieve the closed-loop control of MEMS gyroscopes at different ambient temperatures.

3.2. The Simulation Results of the Closed-Loop Control of MEMS Gyroscopes

3.2.1. The Closed-Loop Simulation Results of Drive Mode

Figure 7 shows the simulation results of drive mode at different temperatures.

Figure 7

The closed-loop simulation results of drive mode: (a) the output voltage of amplitude detection at different temperatures (b) the output voltage of the loop filter at different temperatures.

The simulation results show that the outputs of amplitude detection at different temperatures are all 1 V, which is equal to the reference value. This means that the amplitude of the drive-axis displacement is a fixed value as the excitation voltage is adjusted at different ambient temperatures. Additionally, at different temperatures, the output of the loop filter reaches a stabilized value at about 0.5 s. This means that the vibration frequency of the drive-axis displacement has stabilized. Therefore, the drive mode realizes closed-loop control of the vibration amplitude and frequency at different temperatures. It is noteworthy that the output voltages of the loop filter corresponding to different temperatures are different, which indicates that the stable vibration frequency of the drive-axis displacement is different when the ambient temperature changes. This is because the resonant frequency of drive mode is affected by the ambient temperature.

3.2.2. The Closed-Loop Simulation Results of Sense Mode

For sense mode, the closed-loop detection is based on complete compensation for the Coriolis force. This means that the proof mass is only affected by the quadrature force, which causes the proof mass to always vibrate near the equilibrium position. The closed-loop detection method uses the Proportion-Integration (PI) controller algorithm, which will be specifically explained in Section 3.2.3. Figure 8 shows the simulation result of the resultant force of Coriolis force and feedback force.

Figure 8

The resultant force of Coriolis force and feedback force.

The simulation result shows that the resultant force is gradually changed to Zero at 1 s (Ω = −300°/s), which means that the Coriolis force is completely compensated for by the feedback force. Specifically, the closed-loop detection is achieved by setting up feedback electrodes. The electric potential energy is generated between the feedback electrodes and the proof mass by applying external feedback voltages, so a tangential electrostatic feedback force is generated to compensate for the Coriolis force completely. Therefore, the closed-loop detection is implemented based on this method. Figure 9 shows the output voltage corresponding to different input angular rates.

Figure 9

The curve of input angular rate and output voltage.

Different input angular rates can cause different Coriolis force values, so the feedback force induced by the output voltage also changes. Therefore, there is a proportional relationship between the input angular rate and the output voltage. Figure 9 shows that the scale factor is –0.02459 mV/(°/s) with a linear correlation coefficient of 1. The nonlinearity induced by the parallel-plate sensing mechanism is completely eliminated in closed-loop detection because the proof mass always vibrates near the equilibrium position.

Changes in temperature will cause changes in the motion state of a gyroscope. So, the phase of the sense-axis displacement will also be changed. However, the phase of the demodulated signal used in the closed-loop detection circuit cannot be changed in real-time. Therefore, different temperatures will induce different output voltages when the input angular rate is zero. Figure 10 shows different ZRO voltages of gyroscopes at different temperatures. Additionally, a corresponding cubic fitting curve based on the self-compensation method [21] is obtained as where is the curve fitting results for T, and its unit is μV. The final compensated curve is the difference between the raw simulation result and the fitting result. The final temperature result shows that the compensated ZRO voltages are around zero in the temperature range from −20 to 40 °C.

Figure 10

The Zero Rate Output (ZRO) at different temperatures and the temperature compensation curve.

3.2.3. The Optimization Designs of the Closed-Loop Detection Circuit

It is noteworthy that the choice of the demodulated signal phase is important for the aforementioned closed-loop detection scheme. The purpose of closed-loop detection is to control by applying the feedback voltage. This is achieved indirectly by controlling while also requiring a suitable phase . is the controlled variable of the PI controller. The Coriolis displacement and quadrature displacement can be assumed as where is the amplitude of the Coriolis displacement, and this displacement is the result of the interaction between the Coriolis force and the feedback force, is the amplitude of the quadrature displacement, and is the phase of the sense-axis displacement. The orthogonal demodulation process in the closed-loop detection circuit is shown in Figure 11.

Figure 11

The orthogonal demodulation diagram of closed-loop detection.

Therefore, the controlled variable of the PI controller is where φ1 is the phase of the demodulated signal. When , the second term of Equation (27) is zero. This means that the quadrature displacement has no effect on the output voltage. Therefore, and are in a ratio relationship. The loop can control that equals to zero by controlling that is equal to zero. When , the error term due to the quadrature displacement is not equal to zero. This means that the output voltage is affected by the Coriolis displacement and quadrature displacement. Therefore, cannot be controlled to zero by controlling to zero. Additionally, the Coriolis force cannot be completely compensated for in this case. In consequence, the choice of the demodulated signal phase is important for the optimization of Coriolis force cancellation. Figure 12 shows the compensation for the Coriolis force at different demodulation phases.

Figure 12

The resultant force at different demodulation phase angles.

The simulation results show that when , the resultant force of Coriolis force and feedback force is not equal to zero. Therefore, the Coriolis force has not been completely compensated for in this case. The above analysis and simulation show that the optimization of Coriolis force compensation is achieved by choosing a suitable demodulation signal phase .

In addition to the phase of the demodulation signal, the circuit gain of closed-loop detection will also affect the performance of the circuit. For simplicity, we take the product of and as k for analysis. Figure 13 shows different stabilization times at different k.

Figure 13

The stabilization time at different circuit gain k.

The simulation results show that a suitable circuit gain will induce the shortest stabilization time. The feedback force is proportional to the output voltage. So, the feedback force is constantly changing during the process of stabilization. A large k will cause the overshoot of the output voltage to increase, thus increasing the stabilization time of the feedback force. A small k will also slow down the response process of the feedback force. So, the corresponding stabilization time will also be increased.

4. Discussion

In this paper, a PSPICE device model of parallel-plate capacitive micro-gyroscopes was introduced for reducing the cost of prototype development, which consists of a principle diagram and different modules. Systematic simulations of this model and closed-loop control circuits were conducted, and it is demonstrated that this PSPICE device model has two clear advantages:

Some details of MEMS gyroscopes are considered in this model, including the capacitance model without approximation, mechanical thermal noise, and the effect of temperature.

The closed-loop simulation of MEMS gyroscopes is achieved based on the differential capacitance interfaces of this model and subsequent interface circuits.

Some optimization designs of the circuit have been implemented based on this model, including a complete compensation for Coriolis force and the corresponding stabilization time. The test Printed circuit board (PCB) of closed-loop detection is being designed. The corresponding experiment result will be summarized in future work.