# An Inversely Designed Model for Calculating Pull-In Limit and Position of Electrostatic Fixed-Fixed Beam Actuators

Authors: Cevher Ak; Ali Yildiz
Journal: Mathematical Problems in Engineering
 (2014)

---

## Abstract

This study presents an inverse approach to obtain a relation between applied voltage and displacement of the midpoint of fixed-fixed beam actuator. The approach has two main sections. The first one is the inverse design of a model to replace real action of upper beam under electrostatic force. The formula obtained from the first section does not comprise the residual stress and gives very small errors when there is no residual stress on the upper electrode. So, the second part was carried out to add this important system variable into the formula. Likewise, inverse solution was again applied in the later section. The final formula demonstrates that pull-in limit of clamped-clamped actuator is to be at around 40% of original spacing that is in agreement with simulation and previous experimental results. Its percentage errors are within 2% when compared with simulations that are based on finite element method (FEM).  The results are comparable to numerical solutions received from diverse distributed models which require more calculation power in electrostatic and structural domains. On top of that, our formula is valid for all displacements from original position up to pull-in limit.


---

## Body

## 1. Introduction

Electrostatic MEMS based actuators have been used widely as a sensor due to their higher sensitivity, smaller dimensions, low-power consumption, and easy fabrication with new design possibilities. They have been used as a microelectromechanical varactor [1], as a capacitive pressure sensor for measuring blood pressure for a cardiovascular catheter [2], as a mass sensor [3], as an RF microswitch [4, 5], as an energy harvester for MEMS devices [6], and as small force detection [7].An electrostatically controlled MEMS based fixed-fixed actuator is made up of two parallel conductive beams; the ground electrode is fabricated on a substrate and not movable and the upper electrode is suspended above it with an initial gap(
g
) and fixed from both ends. When a voltage difference is applied between bottom and upper electrodes, the middle section of top electrode moves towards bottom electrode due to electrostatic force. The counteract spring force will stop the motion of upper electrode at some equilibrium point. The spring force is a linear function of the movement whereas the electrostatic force increases with the square of distance. After a certain point, the restoring force due to bending cannot balance the electrostatic force any longer. The upper electrode is instable after this point and collides on the bottom electrode. This limit is called pull-in limit and voltage value is named as pull-in voltage. The real behavior of upper electrode can be seen in Figure 1. It obeys two constraints; fixed ends have zero movement and zero-angle.Figure 1
Electrostatic actuator (side view).Calculating the pull-in voltage value accurately is the most crucial issue in MEMS actuators. Nevertheless, behavior is nonlinear because of electrostatic and mechanical coupling. Hence, getting an analytical formula for pull-in limit is very difficult. For more than two decades, researchers have been developing many models and methods for electrostatic MEMS based actuators to calculate pull-in limit [8–16]. However, most of these works have found different pull-in limits. When it is controlled by software (COMSOL) which employs finite element method, pull-in limit appears to be at approximately 40% of initial gap which is also consistent with former experimental measurements [12, 14]. Table 1 points some simulation results out for varying initial gaps. Upper beam was simulated as linear and isotropic poly-Si material in 3D geometry. Its length was chosen as 200 μm for simulation purposes only since length is not a function of pull-in gap. Ratio of pull-in gap to initial gap is constant and independent of material properties of the upper electrode. Other properties of the material in simulation can be seen in Table 1.Table 1
COMSOL Simulation Pull-in Results of a Fixed-Fixed Beam with length = 200μm, width = 50 μm, thickness = 1 μm, Young Modulus = 160 GPa, Poisson Ratio = 0.22.
Initial Gap (μm)
Pull-in Gap (μm)
Pull-in Gap/Initial Gap
2
0.797
0.3985
4
1.594
0.3985
5
1.992
0.3984
10
3.984
0.3984The main intention of this work is to find a relatively simple formula which demonstrates not only the pull-in limit at the actual position, but also a relationship between applied potential difference and displacement of the middle part of a clamped-clamped electrostatic actuator.


## 2. Inversely Designed Model

We have altered the real behavior of upper electrode with a new model for the sake of simplicity. Upper electrode has subdivided into three regions with four pivots. Two of the pivots are located at both fixed ends. They can only rotate around the pivot points whereas they cannot move horizontally and vertically. The other pivot points are located at3
L
/
7 distances from both fixed sides symmetrically. These pivots are free to move vertically though. The new model can be seen clearly in Figure 2. Total length of the electrode is L. The first and third partitions divided as 3
L
/
7 while the second one as L
/
7. These values have been inversely reached by trial and error method since we already know the true results from simulations and previous experimental studies.Figure 2
The movements and capacitances of the partitions introduced in new model.We had studied bisection model for fixed-free actuators earlier [21, 22]. The model was dividing the upper electrode into two partitions. One is not moving at all and fixed to the substrate from one end and the other section can freely move linearly around a pivot which is selected by trial and error method inversely. It was successful as well when compared with simulation results and experimental studies. Later, model for fixed-free actuators was simplified further and named as Pivot Model [23]. In the model, fixed side of the upper electrode is pivoted and the other side is totally free to move linearly around the pivot point. The upper beam was assumed as a rigid body in the model for the sake of simplicity. These were reasonable assumptions since the bending of the beam was very small. It establishes a good relation between applied voltage and displacement span from rest position to pull-in limit.The new model was developed for clamped-clamped actuators and utilizes only one constraint of the fixed ends of top electrode. Fixed ends still have zero displacement while zero-angle constraint is omitted as in Pivot Model. Thus, only one end of the first and third partitions moves linearly around fixed pivot while the second region is totally free to move vertically. The new design model was named as Inverse Pivot Model (IPM) since pivots were used to make the system uncomplicated. In this model, upper electrode is subdivided into three partitions. Locations of the movable pivots were chosen inversely from simulation results and previous experimental measurements. It was also possible to subdivide the beam into more than 3 partitions in order to have a closer representation of the real bending shape of the beam. However, model would get complicated which is not a desired case. In this study, keeping the model as simple as possible is foremost and initiative intention while it gives a good approximation of the real system. After many comparisons with simulation outcomes and previous measurements, length of the upper beam was subdivided into 7 equal segments. One segment was taken as a middle partition, and the other 2 partitions equally shared the rest of the segments due to symmetry. Therefore, side partitions (partitions 1 and 3) consist of 3 segments.The capacitance calculation of the second partition is easy since it establishes a parallel plate capacitance shape during its movement andC
2 can be calculated as
(1)
C
2
=
1
7
ε
0
w
L
(
g
-
δ
)
,
where ε
0 is permittivity of free space. Capacitances of the first and third partitions are equal because of the symmetry. Therefore, calculating one of them will be enough. However, calculation of it is not simple as second partition.Partial capacitance ofC
1 can be considered as a parallel plate capacitance since the inclination of the upper electrode is negligible for this infinitesimally small region. It can be seen in Figure 3 and can be written as
(2)
d
C
1
=
ε
0
d
A
g
-
y
.
d
A is infinitesimal surface area and equal to w
d
x, where w is the width of the upper beam. y can be found from the geometry shown in Figure 3 as
(3)
y
=
δ
x
3
a
,
where δ is the farthest displacement of the upper beam when a potential difference is applied. It equals the displacement of middle point as well. Therefore, (2) can be rewritten as
(4)
d
C
1
=
3
a
ε
0
w
d
x
3
a
g
-
g
x
.
So, capacitances of the first partitions can be calculated as
(5)
C
1
=
∫
d
C
1
=
3
a
ε
0
w
∫
x
=
0
3
a
d
x
3
a
g
-
g
x
.
When value of a is inserted in (5), C
1 and C
3 can be calculated as
(6)
C
1
=
C
3
=
3
7
ε
0
w
L
δ
ln
⁡
(
g
g
-
δ
)
.Figure 3
Capacitance calculation for the first partition.In a real actuator, electrostatic force is nonuniform and distributed along the upper beam. It is extremely difficult to express it in a simple formula. Therefore, it has been exchanged with a single equivalent electrostatic force term. Since the middle region is the closest part of the upper electrode to bottom one, the electrostatic force has the biggest component at this point. Besides, because of symmetry, electrostatic force is placed at the center of the upper electrode as a single equivalent force. Nevertheless, fringing effect of the capacitances was ignored in order to keep the model simple. To compensate the absence of the fringing effect, restoring forces are adjusted accordingly. Instead of using one restoring force term at the center, it has been split into three. Two restoring forces are placed at the one-third of the linearly moving partitions from fixed ends in order to decrease effect of spring force term to make up the missing fringing effect. The final exact point is again reached after many attempts by inverse approach. Optimization algorithms can also be used to find a more precise location in the future studies. The middle spring force is apparently placed at the center. All equivalent forces and their locations can be seen in Figure4.Figure 4
Locating the equivalent forces of new model.The infinitesimal electrostatic force term can be written as(7)
d
F
e
=
1
2
d
C
d
δ
V
2
,
where C is the total capacitance term of the system and V is potential difference applied between electrodes. The total electrostatic force term can be calculated as
(8)
F
e
=
-
ε
0
w
L
14
×
(
5
δ
2
-
6
δ
g
+
(
6
δ
2
+
6
g
2
-
12
δ
g
)
ln
⁡
(
g
/
(
g
-
δ
)
)
δ
2
(
g
-
δ
)
2
)
V
2
and spring constant k can be obtained for a fixed-fixed beam as [24]
(9)
k
=
32
E
^
w
(
t
L
)
3
,
where
(10)
E
^
=
E
(
1
-
ν
2
)
,
where E, ν, and t are Young’s modulus, Poisson ratio of material, and the thickness of the upper electrode, respectively. Instead of plain E, E
^ has been used to keep the formula valid for wide beams w
≥
5
t [17].Since the 3 restoring springs are parallel to each other, spring constant of each partition is equal tok
/
3. The spring displacement of middle partition is a distance δ, while displacements of the other two partitions are δ
/
3. The total restoring force term can be written for 3 partitions as
(11)
F
y
=
δ
(
k
3
)
+
2
(
δ
3
)
(
k
3
)
=
160
9
δ
E
w
(
1
-
ν
2
)
(
t
L
)
3
.The first term in (11) represents the restoring force of the middle partition spring, and the second term represents the total restoring forces of the side partitions springs.Since the electrostatic and restoring forces have to be equal to each other at equilibrium position, (8) and (11) can be equated to each other as
(12)
-
ε
0
w
L
14
(
5
δ
2
-
6
δ
g
+
(
6
δ
2
+
6
g
2
-
12
δ
g
)
ln
⁡
(
g
/
(
g
-
δ
)
)
δ
2
(
g
-
δ
)
2
)
×
V
2
=
160
9
δ
E
w
(
1
-
ν
2
)
(
t
L
)
3
.
Consequently, the relation between displacement and applied voltage can be obtained from (12) as
(13)
V
=
(
h
h
h
k
i
k
k
h
h
×
(
6
δ
2
+
6
g
2
-
12
δ
g
)
+
5
δ
2
-
6
δ
g
ln
⁡
(
g
g
-
δ
)
)
)
-
1
)
2240
9
(
h
h
h
k
i
k
k
h
h
×
(
6
δ
2
+
6
g
2
-
12
δ
g
)
+
5
δ
2
-
6
δ
g
ln
⁡
(
g
g
-
δ
)
)
)
-
1
-
E
t
3
δ
3
(
δ
-
g
)
2
g
k
k
k
k
h
k
k
h
g
×
(
h
h
h
k
i
k
k
h
h
×
(
6
δ
2
+
6
g
2
-
12
δ
g
)
+
5
δ
2
-
6
δ
g
ln
⁡
(
g
g
-
δ
)
)
ε
0
L
4
(
1
-
ν
2
)
y
y
i
y
l
l
l
l
l
h
l
l
l
l
l
l
i
y
y
×
(
ln
⁡
(
g
g
-
δ
)
(
6
δ
2
+
6
g
2
-
12
δ
g
)
h
h
h
k
i
k
k
k
k
k
k
k
k
k
h
h
+
5
δ
2
-
6
δ
g
ln
⁡
(
g
g
-
δ
)
)
)
-
1
)
)
1
/
2
.Equation (13) is valid not only for pull-in limit, but also for all values within the pull-in limit. A comparison of simulation results received from COMSOL package and our IPM results for different displacements up to pull-in limit boundary. All results seem within 1% error level.Many fixed-fixed actuator sensors employ only on/off position like a switch. Therefore, determination of pull-in voltage is very crucial for an actuator before real device fabrication in order to produce a proper sensor which is working in a range of interest. In order to retrieve the pull-in limit, derivative ofV in (13) with respect to δ has to be taken and equate the result to zero and this position is named as pull-in limit [9]. The upper beam will be unstable and collapse towards bottom beam after this critical point. Since the derivative equation equals zero, we can use only the numerator of the derivative since denominator is not affecting the result. So, denominator is dropped from the derivative equation. Subsequently, some common variables are also canceled in order to keep the calculation simple:
(14)
18
ln
⁡
(
g
g
-
δ
)
(
g
-
δ
)
3
-
δ
(
2
g
-
3
δ
)
(
9
g
-
7
δ
)
=
0
.Unfortunately, having an analytical solution to this equation is very cumbersome. Thence, a computational result has been obtained by iteration as(15)
δ
Pull-in
Limit
≅
0.4
g
.This value also can be seen in Table2. When δ
/
g gets closer to 0.4, voltage value starts to saturate. When this important value was inserted into (13) back, crucial value of potential difference can be found as
(16)
V
Pull-in
Limit
=
11.5486
E
g
3
t
3
(
1
-
ν
2
)
ε
0
L
4
.Table 2
Comparison of Voltage Values for Arbitrary Displacements.E
=
169 GPa, L
=
350 μm, t
=
3 μm, ν
=
0.06 and g
=
1 μm.
Displacement (μm)
Voltage (Volt)IPM (9)
Voltage (Volt)COMSOL
Δ%
0.0116
4.9518
5
0.97
0.0493
9.9076
10
0.93
0.0739
11.8935
12
0.90
0.1062
13.8822
14
0.85
0.1496
15.8724
16
0.80
0.1775
16.8617
17
0.82
0.2125
17.8567
18
0.80
0.2594
18.8394
19
0.85
0.3445
19.8053
20
0.98
0.3891
19.9634
20.19
1.14This value is the potential difference value at the pull-in limit boundary. Applying a voltage difference higher than this value causes the upper electrode to be unstable and collapse onto the bottom electrode.


## 3. Inverse Pivot Modified Formula

When it is checked with simulation results obtained by a different software package which utilize finite element method from a previous study, it can be seen that formula gives small errors if there is no residual stress. However, error goes as high as 45% whenever residual stress gets involved since our Inverse Pivot Model does not take residual stress into account. Effect of residual stress can be seen in Table3. Thence, we decided to improve the formula by considering residual stress as a system parameter. In this second section, formula is merely modified inversely since we have a starting point (equation (16)) and the true results from previous researches [12, 15, 17] to compare outcomes of the new formula. Adding residual stress into (16) is easier than starting from scratch to get a new model. From Tables 3 and 4, it can be apparently seen that positive residual stress increases whereas negative one decreases pull-in voltage. Nevertheless, effect of residual stress should increase with the thickness and decrease with the length of upper electrode. We had to move some physical and system parameters from numerator to denominator in order to keep formula dimensionally correct. After many trials, Inverse Pivot Modified Formula (IPMF) has been reached as
(17)
V
=
(
(
σ
t
(
t
/
L
)
1.1
+
(
65
/
9
)
(
E
g
t
3
/
(
1
-
ν
2
)
L
3
)
)
(
61
/
99
)
(
ε
0
L
/
g
2
)
)
1
/
2
,
where σ is the residual stress. Equation (17) has not been attained as an analytical calculation. However, when the residual stress is assumed as zero in this equation, it closely approaches (16). It also establishes very good results for actuators with residual stress at pull-in limit boundary. Errors are now in acceptable ranges which all are smaller than 2%, even for the worst case (see Tables 3 and 4).Table 3
Comparison of IPM and IPMF with previous simulation results.L
=
250 μm, E
=
169 GPa, ν
=
0.06, w
=
50 μm, t
=
3 μm and g
=
1 μm in each case.
Residual Stress,σ, (MPa)
0
100
−25
CoSolve  [17]
40.10
57.60
33.60
V (2D)  [17] (Δ%)
39.50 (1.50%)
56.90 (1.22%)
33.70 (0.30%)
OCI Model  [12] (Δ%)
39.60 (1.25%)
57.40 (0.35%)
33.71 (0.33%)
GDQM  [15] (Δ%)
39.13 (2.41%)
57.62 (0.03%)
33.63 (0.09%)
IPM (Δ%)
39.13 (2.41%)
39.13 (32.07%)
39.13 (16.461%)
IPMF (Δ%)
39.42 (1.70%)
57.03 (0.99%)
33.61 (0.03%)Table 4
Comparison of IPM and IPMF with previous simulation results.L
=
350 μm, E
=
169 GPa, ν
=
0.06, w
=
50
μm, t
=
3
μm and g
=
1
μm in each case.
Residual Stress,σ, (MPa)
0
100
−25
CoSolve  [17]
20.30
35.80
13.70
V (2D)  [17] (Δ%)
20.20 (0.50%)
35.40 (1.11%)
13.80 (0.73%)
OCI Model [12] (Δ%)
20.20 (0.50%)
35.91 (0.31%)
13.71 (0.07%)
GDQM  [15] (Δ%)
20.36 (0.30%)
35.99 (0.53%)
13.60 (0.73%)
IPM (Δ%)
19.96 (1.67%)
19.96 (44.25%)
19.96 (45.69%)
IPMF (Δ%)
20.11 (0.94%)
35.25 (1.54%)
13.97 (1.97%)In Tables3 and 4, Poisson ratios are 0.06 in all cases. So, our model was also checked with a different Poisson ratio which is 0.32. From Table 5, it can be seen that IPM delivers a comparable result with other studies. In particular, IPMF demonstrates very similar error levels when compared with previous methods which depends on numerical distributed models. These kinds of model require more computing power and need more time to get the result. However, our formula is just one step calculation if it is written in an Excel sheet.Table 5
Comparison of IPM and IPMF with previous simulation results [12]. L
=
250 μm, E
=
169 GPa, ν
=
0.32, w
=
50 μm, t
=
3 μm and g
=
1
μm in each case.
Residual Stress,σ, (MPa)
0
CoSolve FEA  [17]
41.20
V (2D)  [17] (Δ%)
41.50 (0.73%)
[18] (Δ%)
42.54 (3.25%)
[19] (Δ%)
41.20 (0.00%)
[20] (Δ%)
41.42 (0.53%)
OCI Model [12]
41.72 (1.26%)
IPM (Δ%)
41.53 (0.80%)
IPMF (Δ%)
41.23 (0.07%)The IPMF is also checked for some other parameters in order to explore the applicable range of the formula. Table6 shows simulation results obtained in COMSOL for a poly-Si beam whose material properties are selected from COMSOL’s library. Mesh numbers are selected automatically for each case in the software.Table 6
Pull-in Voltage Comparison of IPMF to simulation results for Poly-Si material from COMSOL library.E
=
160 GPa, ν
=
0.22 and w
=
50 μm for all cases.
Residual Stress,σ, (MPa)
COMSOL (V)
IPMF (V)
(Δ%)
L
=
100 μm, t
=
2
μm, g
=
1
μm, Number of Meshes = 7533
0
134.70
132.53
1.61
−25
130.30
127.74
1.96
100
150.00
150.18
0.12
L
=
200
μm, t
=
2 μm, g
=
1
μm, Number of Meshes = 4467
0
33.60
33.13
1.40
−25
29.00
28.40
2.07
100
47.50
47.55
0.11
L
=
300
μm, t
=
3
μm, g
=
2
μm, Number of Meshes = 3531
0
78.30
76.52
2.27
−25
67.90
68.49
0.87
100
109.00
102.51
5.95
L
=
400
μm, t
=
3
μm, g
=
2
μm, Number of Meshes = 2817
0
44.01
43.04
2.20
−25
34.00
34.88
2.59
100
72.00
66.30
7.92
L
=
500
μm, t
=
4
μm, g
=
2
μm, Number of Meshes = 3259
0
42.83
42.41
0.98
−25
33.30
32.71
1.77
100
68.10
68.64
0.79
L
=
750
μm, t
=
5
μm, g
=
3
μm, Number of Meshes = 1782
0
49.10
48.39
1.45
−25
32.80
34.99
6.68
100
87.00
82.53
5.14Table6 demonstrates comparison between IPMF and simulation results for wide range of material parameters. IPMF gives very satisfactory results with a maximum error level of 8% for residual stress case. In Table 7, width is changed for 200 μm and 500 μm beam lengths while other parameters are kept constant to show the effect of width. IPMF again shows very satisfactory results for the range of width with a maximum error level of 3.5%. In Table 8, initial gap is changed for again 200 μm and 500 μm beam lengths while other parameters are kept constant to show the effect of the gap between electrodes. It can be seen that error level gets increase as the initial gap heightens for both 200 μm and 500 μm beam lengths. This error stems from ignorance of fringing effect between upper and bottom electrodes. As the initial gap rises, the capacitance value of the system also proliferates due to additional increase of the fringing effect.Table 7
Comparison of IPMF to simulation results for the variety of widths.E
=
160 GPa, ν
=
0.22 and w
=
50 μm for all cases.
Width,w, (μm)
COMSOL (V)
IPMF (V)
(Δ%)
L
=
200 μm, t
=
2 μm, g
=
1 μm
50
33.60
33.13
1.40
20
33.30
33.13
0.51
10
33.10
33.13
0.09
5
32.90
33.13
0.70
2
32.60
33.13
1.63
1
32.20
33.13
2.89
L
=
500 μm, t
=
4 μm, g
=
2
μm
50
42.83
42.41
0.98
20
42.30
42.41
0.26
10
42.10
42.41
0.74
5
41.80
42.41
1.46
2
41.40
42.41
2.44
1
41.00
42.41
3.44Table 8
Comparison of IPMF to simulation results for the variety of initial gaps.E
=
160 GPa, ν
=
0.22 and w
=
50 μm for all cases.
Gap,g, (µm)
COMSOL (V)
IPMF (V)
(Δ%)
L
=
200 μm, w
=
50 μm, t
=
2
μm
1
33.60
33.13
1.40
2
99.50
93.71
5.82
3
196.00
172.17
12.16
5
520.00
370.45
28.76
L
=
500
μm, w
=
50
μm, t
=
4
μm
2
42.83
42.41
0.98
3
79.80
77.91
2.37
4
126.00
119.95
4.81
6
249.00
220.38
11.49


## 4. Conclusions

In the present study, we have proposed a new inversely designed model (IPM) and modified formula (IPMF) which both deliver pull-in limit of a fixed-fixed beam actuator at around 40% of the initial gap. They are consistent with COMSOL simulation results (Tables1 and 2) and previous experimental and distributed models (Tables 3, 4, and 5). Besides, it establishes a good relation between applied voltage and displacements commencing from rest position up to pull-in limit. Particularly, our inversely designed modified formula yields percentage errors less than 2%, even for the worst case. IPMF is a plain formula which does not demand much computing power. Moreover, IPM is successful with gratifying outcomes for applied voltages for given displacements. When compared with previous experimental measurements, IPM and IPMF can be found considerably successful too. However, if the upper electrode has residual stress, our first model IPM cannot deliver good results because model does not take residual stress into account. On the other hand, IPMF still demonstrates very small percentage errors which are comparable to previous distributed models. IPMF also attests successful outcomes for a wide range of top beam geometric parameters (Tables 6 and 7). The formula is valid for both long beams where L
≥
5
w and narrow beams (w
>
5
t) [12].The most eminent benefit of IPMF is that it establishes a good approximation of the system with an unsophisticated plain formula. One can simply calculate required voltage for pull-in limit rather than utilize numerical distributed methods which requires more computing power and is time consuming.Although IPM and IPMF are outstandingly precise especially around pull-in limit, they have some limitations since they neglect some physical parameters of the system for sake of simplification of the final formula. Firstly, both of them do not take the fringing effect into account. It causes the actual capacitance to be higher than what we have calculated in our Inverse Pivot Model. Therefore, when the initial gap gets higher, IPMF also deviates from acceptable error levels (Table8).Secondly, the upper beam was also assumed as a rigid body in the model. Thence, length of the beam was taken as constant even under significant electrostatic force. Therefore, both axial and transverse stresses would not be formed.Lastly, the models assume the actuators are in vacuum since they ignore any atmospheric pressure on the upper electrode.IPMF can be improved more by applying artificial optimization techniques and fringing effect of the capacitance geometry can be taken into account to make the formula deliver smaller error level.



---

*Source: 391942-2014-08-20.xml*