Models - Mems.biased Resonator 2d Basic
Models - Mems.biased Resonator 2d Basic
3a
Model Definition
The model consists of a poly-silicon resonator, which is manufactured through a surface micromachining process. Initially a silicon wafer is coated with 0.75 m of oxide and 0.15 m of silicon nitride, to isolate the micromachined parts from the wafer ground plane. Polysilicon electrodes with a thickness of 0.3 m are deposited next. A sacrificial layer of oxide is then deposited to a thickness of 1.985 m. Note that in Ref. 1 the sacrificial oxide is actually 1.3 m, but the gap thickness was adjusted to this value for the purposes of simulation, to account for the depletion layer in the silicon. In this model we use the same adjustment to enable the simulations to be directly compared with those presented in the paper. Holes are etched in the sacrificial layer (to provide anchor points for the resonator) and the structural polysilicon is deposited with a thickness of 1.9 m. The geometry is modeled in COMSOL in two dimensions, with a plane strain assumption made for the mechanical deformation. The structure has a plane of symmetry, so it is possible to model only half of the geometry explicitly. However, since we are interested in the normal modes, modeling only half of the structure and assuming symmetry will eliminate all the anti-symmetric modes from the analysis. In this case the full structure is therefore analyzed in COMSOL. Figure 1 shows the geometry.
2012 COMSOL
1 |
Figure 1: The model geometry. The wafer itself is not drawn explicitly, but is represented in the model by a ground plane on the underside of the geometry. The layers of deposited material from the ground plane up are: silicon oxide, silicon nitride (too thin to see clearly), polysilicon electrodes/air gap (etched sacrificial oxide), polysilicon resonator, air. In operation both the silicon resonator and the underlying wafer are grounded and an electric voltage is applied to the driving electrode, which is bisected by the symmetry plane. Typically a DC bias of 35 V is applied in normal operation of the device. The assumption is made that the polysilicon is a perfect conductor, so the bias voltage is applied on all exterior surfaces of the resonator and its anchor as a potential boundary condition. In this model the deformation of the structure is computed with the applied DC bias. Note that the silicon oxide and nitride are assumed to be rigid for the solid mechanics simulations, so the structure is anchored at the base of its electrode, and these domains are not included in the solid mechanics equations.
ELECTROMECHANICAL FORCES
Within a vacuum or other medium, forces between charged bodies can be computed on the assumption that a fictitious state of stress exists within the field. The Electromagnetic or Maxwell stress tensor can be used to compute the induced stresses in a material as a result of an electric field as well as surface forces acting on bodies in air or vacuum. Within a material, COMSOL Multiphysics uses the following form of the stress tensor TEM,S, which is appropriate for isotropic materials (Ref. 2): 1 - ( E D + a 2 E E ) I + ED T + 1 -- ( a a 1 ) EE T T EM, S = -2 2 2 where E is the electric field, D is the electric displacement field, I is the identity tensor, and 0 is the permittivity of free space and a1 and a2 are material parameters that specify the electrostrictive properties of the material (for this device, assume a1 = a2 = 0 because the field is in any case very low within the material). This additional stress is applied to the material by the electromechanical solid node. Note that mechanical stresses are usually induced in the material as a result of the net forces acting on the surfaces, in addition to the stress induced by the electric field. The forces on the surfaces of a solid body can be computed by applying a similar stress term within the vacuum of the form:
2 |
1 - ( E D ) I + ED T T EM, V = -2 A net force on the surface typically results from the discontinuity of the stress tensor at the interface. However, since it is undesirable to apply a stress term throughout the vacuum, the force is only available on the surface of solid bodies, via the electromechanical interface node. The surface force is given by: 1 -n 1 T EM, V = 2 E D n 1 + ( n 1 E ) D where n1 is the surface normal, pointing out from the mechanical body.
2012 COMSOL
3 |
Figure 2: The y-displacement of the resonator as a function of position. The maximum displacement occurs in the center of the resonator, immediately over the biasing electrode.
Figure 3: Electric potential contours in the gap between the grounded resonator and the biased driving electrode.
4 |
References
1. F.D. Bannon III, J.R. Clark and C.T.-C. Nguyen, High-Q HF Microelectromechanical Filters, IEEE Journal of Solid State Circuits, vol. 35, no. 4, pp. 512526, 2000. 2. J.A. Stratton, Electromagnetic Theory, Cambridge (MA), 1941.
Modeling Instructions
MODEL WIZARD
1 Go to the Model Wizard window. 2 Click the 2D button. 3 Click Next. 4 In the Add physics tree, select Structural Mechanics>Electromechanics (emi). 5 Click Next. 6 Find the Studies subsection. In the tree, select Preset Studies>Stationary. 7 Click Finish.
GEOMETRY 1
1 In the Model Builder window, under Model 1 click Geometry 1. 2 In the Geometry settings window, locate the Units section. 3 From the Length unit list, choose m.
Import 1
1 Right-click Model 1>Geometry 1 and choose Import. 2 In the Import settings window, locate the Import section. 3 Click the Browse button. 4 Browse to the models Model Library folder and double-click the file
biased_resonator_2d.mphbin.
2012 COMSOL
5 |
Parameters
1 In the Model Builder window, right-click Global Definitions and choose Parameters. 2 In the Parameters settings window, locate the Parameters section. 3 In the table, enter the following settings:
Name Vdc Expression 35[V] Description DC bias voltage
Explicit 1
1 In the Model Builder window, under Model 1 right-click Definitions and choose Selections>Explicit. 2 Select Domains 5, 6, and 9 only. 3 Right-click Model 1>Definitions>Explicit 1 and choose Rename. 4 Go to the Rename Explicit dialog box and type Resonator in the New name edit field. 5 Click OK.
Adjacent 1
1 Right-click Definitions and choose Selections>Adjacent. 2 In the Adjacent settings window, locate the Input Entities section. 3 Under Input selections, click Add. 4 Go to the Add dialog box. 5 In the Input selections list, select Resonator. 6 Click the OK button. 7 Right-click Model 1>Definitions>Adjacent 1 and choose Rename. 8 Go to the Rename Adjacent dialog box and type Resonator Boundaries in the New name edit field.
6 |
9 Click OK.
Explicit 2
1 Right-click Definitions and choose Selections>Explicit. 2 Select Domain 8 only. 3 Right-click Model 1>Definitions>Explicit 2 and choose Rename. 4 Go to the Rename Explicit dialog box and type Drive Electrode in the New name
edit field.
5 Click OK.
Adjacent 2
1 Right-click Definitions and choose Selections>Adjacent. 2 In the Adjacent settings window, locate the Input Entities section. 3 Under Input selections, click Add. 4 Go to the Add dialog box. 5 In the Input selections list, select Drive Electrode. 6 Click the OK button. 7 Right-click Model 1>Definitions>Adjacent 2 and choose Rename. 8 Go to the Rename Adjacent dialog box and type Drive Electrode Boundaries in
Material Browser
1 In the Model Builder window, under Model 1 right-click Materials and choose Open Material Browser. 2 In the Material Browser window, locate the Materials section. 3 In the tree, select MEMS>Semiconductors>Poly-Si. 4 Right-click and choose Add Material to Model from the menu. 5 In the Model Builder window, right-click Materials and choose Open Material Browser. 6 In the Material Browser window, locate the Materials section.
2012 COMSOL
7 |
7 In the tree, select MEMS>Insulators>SiO2. 8 Right-click and choose Add Material to Model from the menu.
SiO2
1 In the Model Builder window, under Model 1>Materials click SiO2. 2 Select Domain 1 only.
Material Browser
1 In the Model Builder window, right-click Materials and choose Open Material Browser. 2 In the Material Browser window, locate the Materials section. 3 In the tree, select MEMS>Insulators>Si3N4. 4 Right-click and choose Add Material to Model from the menu.
Si3N4
1 In the Model Builder window, under Model 1>Materials click Si3N4. 2 Select Domain 2 only.
Material Browser
1 In the Model Builder window, right-click Materials and choose Open Material Browser. 2 In the Material Browser window, locate the Materials section. 3 In the tree, select Built-In>Air. 4 Right-click and choose Add Material to Model from the menu.
Air
1 In the Model Builder window, under Model 1>Materials click Air. 2 Select Domains 3, 4, 7, and 10 only.
Set the air material property to be non-solid, to ensure the interface solves the electrostics equations in the spatial frame.
3 In the Material settings window, click to expand the Material Properties section. 4 From the Domain type list, choose Non-solid.
Set up the solid mechanics and electrostatics boundary conditions. Add a selection to the default linear elastic dielectric feature.
ELECTROMECHANICS
8 |
2 In the Linear Elastic Dielectric settings window, locate the Domain Selection section. 3 From the Selection list, choose Resonator.
Fixed Constraint 1
1 In the Model Builder window, right-click Electromechanics and choose the boundary
Use a fixed mesh in the solid dielectric domains that are assumed to be rigid.
Fixed Mesh 2
1 Right-click Electromechanics and choose the domain setting Deformed Mesh>Fixed Mesh. 2 Select Domains 1, 2, and 8 only.
The default mesh boundary conditions are suitable for this model, so none need to be added. Apply electric potential boundary conditions on grounded surfaces.
Electric Potential 1
1 Right-click Electromechanics and choose the boundary condition Electrical>Electric Potential. 2 Select Boundaries 2 and 42 only.
Electric Potential 2
1 Right-click Electromechanics and choose the boundary condition Electrical>Electric Potential. 2 In the Electric Potential settings window, locate the Boundary Selection section. 3 From the Selection list, choose Resonator Boundaries.
Electric Potential 3
1 Right-click Electromechanics and choose the boundary condition Electrical>Electric Potential. 2 In the Electric Potential settings window, locate the Boundary Selection section. 3 From the Selection list, choose Drive Electrode Boundaries.
2012 COMSOL
9 |
4 Locate the Electric Potential section. In the V0 edit field, type Vdc.
Modify the default mesh settings to improve the mesh resolution in the narrow gaps.
MESH 1
Size
1 In the Model Builder window, under Model 1 right-click Mesh 1 and choose Edit Physics-Induced Sequence. 2 In the Model Builder window, under Model 1>Mesh 1 click Size. 3 In the Size settings window, locate the Element Size section. 4 Click the Custom button. 5 Locate the Element Size Parameters section. In the Resolution of narrow regions edit
1 In the Model Builder window, right-click Study 1 and choose Rename. 2 Go to the Rename Study dialog box and type Stationary in the New name edit field. 3 Click OK. 4 Right-click Study 1 and choose Compute.
Check that the renamed study node is referenced correctly by the X-file.
RESULTS
Data Sets
Plot the y-displacement of the structure.
Displacement (emi)
1 In the Model Builder window, under Results right-click Displacement (emi) and choose Rename. 2 Go to the Rename 2D Plot Group dialog box and type Biased Displacement in the New name edit field. 3 Click OK.
Biased Displacement
1 In the Model Builder window, expand the Results>Biased Displacement node, then
click Surface 1.
10 |
2 In the Surface settings window, locate the Expression section. 3 In the Expression edit field, type v. 4 Click the Plot button. 5 Click the Zoom Extents button on the Graphics toolbar.
Compare the resulting plot with that in Figure 2. Create a plot to show the electric potential contours in the gap.
2D Plot Group 3
1 In the Model Builder window, right-click Results and choose 2D Plot Group. 2 Right-click 2D Plot Group 3 and choose Contour. 3 In the Contour settings window, locate the Expression section. 4 In the Expression edit field, type V. 5 Locate the Levels section. In the Total levels edit field, type 10. 6 Click the Plot button. 7 In the Model Builder window, right-click 2D Plot Group 3 and choose Rename. 8 Go to the Rename 2D Plot Group dialog box and type Electric Potential
Contours in the New name edit field.
9 Click OK.
View 2
1 In the Model Builder window, under Model 1 right-click Definitions and choose View. 2 In the View settings window, locate the View section. 3 Select the Lock axis check box.
Axis
1 In the Model Builder window, expand the View 2 node, then click Axis. 2 In the Axis settings window, locate the Axis section. 3 In the x minimum edit field, type -12. 4 In the x maximum edit field, type 12. 5 In the y minimum edit field, type -2. 6 In the y maximum edit field, type 4. 7 Click the Apply button.
2012 COMSOL
11 |
12 |