Popular Categories
Control Systems Home
Classical Control
Controller Design
Examples
System Modeling
MATLAB
Control Systems wiki

The online resource for Control Systems




Simplified Lumped Parameter MEMS Tuning Fork Gyroscope Model

From Control Systems at ControlTheoryPro.com, your resource for all things controls.

Jump to: navigation, search
Simplified Lumped Parameter MEMS Tuning Fork Gyroscope Model
Sensors MEMS
To prevent spam, users must register before they can edit/create articles.
There are 195 articles and 323 users of the wiki site.
Authors





Gabe Spradlin, founder of ControlTheoryPro.com, now offers Consulting Services for your Aerospace Systems and Controls needs.
Call me at (720) 210-4078.

Contents

Introduction to Simplified Lumped Parameter MEMS Tuning Fork Gyroscope Model

Figure 1:  Generic 2-axis mass-spring-damper model
Figure 1: Generic 2-axis mass-spring-damper model
See Also




The basics of vibratory MEMS gyro are presented in Vibratory MEMS Gyroscopes. The standard modeling technique is presented in MEMS Gyro Modeling.

Typical is described concisely in Suhas et all[1], quote

 Among various MEMS sensors, a rate gyroscope is one of the most complex 
 sensors from the design point of view.  The gyro normally consists of a 
 proof mass suspended by an elaborate assembly of beams that allow the 
 system to vibrate in two transverse modes.
 In most cases, the FEM analysis becomes prohibitive and one resorts to 
 equivalent electrical circuit simulations...  Here, we present a 
 simplified lumped parameter model of the tuning fork gyro and show how 
 easily it can be implemented using a generic tool like Simulink.  The 
 results obtained are compared with those obtained from more elaborate 
 and intense simulations...  The comparison shows that lumped parameter 
 Simulink model gives equally good results...

The standard proof mass model is shown in Figure 1.

Main Design Parameters for Vibratory Gyroscopes

The main design parameters are stiffness of the structure, damping coefficients and system response to applied rate. The stiffness is key to mode matching between the resonant axes and the driving frequency. Achieving the desired performance is an iterative modeling and simulation process.

Gyroscope Structures

As described in Vibratory MEMS Gyroscopes, mode matching of the resonance is key to maximizing the signal to noise ratio (SNR). The quality (Q) factor is a measure of how large the resonant amplitude will be; when the sense and drive resonances are equal, the output signal is amplified be the Q factor. The larger the Q factor the better the SNR will be. The gyro structure can have a large impact on the Q factor.

Structure with in-plane drive and sense modes

This gyro structure type is symmetric and decoupled having both the drive and sense modes in-plane. In this structure the sensing is done with a pair of comb fingers identical to the fingers used for driving the structure. For more details on this structure and design see Suhas et all[2].

Modeling and Simulation of the Tuning Fork Gyroscope

As seen in Figure 1, the suspended proof mass acts as a single degree of freedom (DOF) system for both the drive and sense motions. The lumped parameters of the structure, mass and stiffness, are obtained using energy methods. The relevant boundary conditions are of the suspended beam model are the anchored end and the free end connected to the proof mass.

The lumped parameters are

LaTeX: M_%7Bd%7D%3DM_%7Bs%7D%3Dm_%7Bp%7D%2Bm_%7Bb%7D%2Bm_%7Bc%7D-m_%7Be%7D


where:

LaTeX: M_%7Bd%7D%20%2F%20M_%7Bs%7D is the mass of the structure,
LaTeX: m_%7Bp%7D is the plate mass,
LaTeX: m_%7Bb%7D is the beam mass,
LaTeX: m_%7Bc%7D is the comb mass, and
LaTeX: m_%7Be%7D is the etch hole mass.
LaTeX: K_%7Bd%7D%3D%5Cfrac%7B4Et%5Comega_%7Bb%7D%5E%7B3%7D%7D%7Bl_%7Bb%7D%5E%7B3%7D%7D


where:

LaTeX: K_%7Bd%7D is the in plane stiffness.
LaTeX: K_%7Bs%7D%3D%5Cfrac%7B4Et%5E%7B3%7D%5Comega_%7Bb%7D%7D%7Bl_%7Bb%7D%5E%7B3%7D%7D


where:

LaTeX: K_%7Bs%7D is the out of plane stiffness.
LaTeX: B_%7Bd%7D%3D%5Cfrac%7B%5Cmu%20A%7D%7B%5Cleft%28%201%2B2K_%7Bn%7D%20%5Cright%29g%7D for rarefied flow


where:

LaTeX: B_%7Bd%7D is the slide film damping.
LaTeX: B_%7Bs%7D%3D%5Cfrac%7B3%5Cmu%20%5Cleft%28%20L_%7Bp%7D%20%5Cright%29%5E4%20n_%7Bh%7D%7D%7B8g%5E3%7D%20%5Cleft%20%5B%204ln%20%5Cleft%28%20%5Ceta%20%5Cright%29%20-%203%2B%5Cfrac%7B4%7D%7B%5Ceta%5E2%7D-%5Cfrac%7B1%7D%7B%5Ceta%5E4%7D%20%5Cright%20%5D%20%2B%20%5Cfrac%7B8%5Cpi%20%5Cmu%20t%20n_%7Bh%7D%20%5Cleft%28%20L_%7Bp%7D%5E%7B2%7D%20-%20L_%7Bh%7D%5E%7B2%7D%20%5Cright%29%5E2%7D%7BL_%7Bp%7D%5E%7B4%7D%7D


where:

LaTeX: B_%7Bs%7D is the squeeze file damping,
LaTeX: L_%7Bp%7D is the pitch,
LaTeX: L_%7Bh%7D is the edge length of hole, and
LaTeX: %5Ceta%3D%5Cfrac%7BL_%7Bp%7D%7D%7BL_%7Bh%7D%7D.

The in plane motion and the push-pull driving circuit is governed by the equation

LaTeX: M_%7Bd%7D%5Cddot%7BY%7D%2BB_%7Bd%7D%5Cdot%7BY%7D%2BK_%7Bd%7DY%3DF_%7Be%7D (1)


where:

LaTeX: F_%7Be%7D is teh electrostatic driving force defined by
LaTeX: F_%7Be%7D%3D%5Cfrac%7B2.28n%5Cepsilon%20t%20V_%7Bdc%7D%20V_%7Bac%7D%20%5Cmbox%7B%20sin%7D%20%5Cleft%28%20%5Comega_%7Bd%7D%20t%20%5Cright%29%7D%7Bg_%7Bc%7D%7D


When a rotation rate Ωr is applied about the rotational axis. The Coriolis force causes motion along the sense direction (normal to the plane of excitation). The motion caused by the Coriolis force is

LaTeX: M_%7Bs%7D%5Cddot%7BZ%7D%2BB_%7Bs%7D%5Cdot%7BZ%7D%2BK_%7Bs%7DZ%3DF_%7Bc%7D (2)


where the Coriolis force is

LaTeX: F_%7Bc%7D%3D2M_%7Bs%7D%5COmega_%7Br%7D%5Cdot%7BY%7D%5Cmbox%7B%20sin%7D%20%5Cleft%28%20%5Comega_%7Bd%7D%20t%20%5Cright%29


The (linearized) change in capacitance is

LaTeX: %5CDelta%20C_%7Bb%7D%3D%5Cfrac%7B%5Cepsilon%20A_%7Be%7D%7D%7Bg%5E2%7DZ (3)


Reduced order models can be developed using the governing equations (Eqn. 1-3). The model can be developed to design and predict the performance in the time and frequency domain.

LaTeX: %5Cfrac%7B%5CDelta%20C_%7Bb%7D%7D%7B%5COmega_%7Br%7D%7D%3D%5Cfrac%7B%5Cleft%28%20G_%7B1%7D%20%5Cbar%7BG_%7B2%7D%7D%20G_%7B3%7D%20%5Cright%29%7D%7B%5Cleft%20%5C%5B%20M_%7Bd%7DM_%7Bs%7Ds%5E4%20%2B%20%5Cleft%28%20M_%7Bd%7DB_%7Bs%7D%20%2B%20M_%7Bs%7DB_%7Bd%7D%20%5Cright%29s%5E3%20%2B%20%5Cleft%28%20M_%7Bd%7DK_%7Bs%7D%20%2B%20M_%7Bs%7DK_%7Bd%7D%20%2B%20B_%7Bd%7DB_%7Bs%7D%20%5Cright%29s%5E2%20%2B%20%5Cleft%28%20B_%7Bd%7DK_%7Bs%7D%20%2B%20B_%7Bs%7DK_%7Bd%7D%20%5Cright%29s%20%2B%20K_%7Bd%7DK_%7Bs%7D%20%5Cright%20%5C%5D%7D (4)


where:

LaTeX: G_%7B1%7D%3D2M_%7Bs%7D,
LaTeX: %5Cbar%7BG_%7B2%7D%7D%3D%5Cfrac%7B2.28n%5Cepsilon%20t%20V_%7Bdc%7D%20V_%7Bac%7D%7D%7Bg_%7Bc%7D%7D and
LaTeX: G_%7B3%7D%3D%5Cfrac%7B%5Cepsilon%20A_%7Be%7D%7D%7Bg%5E2%7D
LaTeX: %5Cfrac%7B%5CDelta%20C_%7Bb%7D%7D%7B%5COmega_%7Br%7D%7D%3D%5Cfrac%7B2YQ_%7Bs%7D%7D%7B%5Comega_%7Bs%7D%7DG_%7B3%7D (5)


where:

LaTeX: Y is the amplitude of the drive mode and
LaTeX: Q_%7Bs%7D%3D%5Cfrac%7B%5Csqrt%7BK_%7Bs%7DM_%7Bs%7D%7D%7D%7BB_%7Bs%7D%7D is the quality factor of the sense mode.



Resources

Notes

  1. Mohite et all, Abstract
  2. Mohite et all, pg. 759
Creative Commons License
Except where otherwise noted, content on this site is licensed under a Creative Commons Attribution 3.0 License
© 2008 ControlTheoryPro.com, LLC
Designed by ControlTheoryPro.com
Powered by MediaWiki