# A Novel Method for Investigating the Casimir Effect on Pull-In Instability of Electrostatically Actuated Fully Clamped Rectangular Nano/Microplates

Authors: Arman Mohsenzadeh; Masoud Tahani; Amir R. Askari
Journal: Journal of Nanoscience
 (2015)

---

## Abstract

The objective of the present paper is to represent a novel method to investigate the stable and unstable behaviors of fully clamped rectangular nano/microplates under the effects of electrostatic and Casimir pressures. To this end, the governing partial differential equation of equilibrium is considered and reduced to an algebraic equation using a simple and computationally efficient single degree of freedom (SDOF) model through the Galerkin weighted residual method. The linear and undamped mode-shapes of the plate are used in the Galerkin procedure as the weight function which is obtained by the extended Kantorovich method (EKM). The present findings are compared and validated by available empirical and theoretical results in the literature as well as those obtained by finite element (FE) simulation carried out using COMSOL Multiphysics commercial software and excellent agreements between them are observed.


---

## Body

## 1. Introduction

Micron and submicron scale structures are frequently used in different applications such as nano/microelectromechanical systems (N/MEMS) nowadays. Because of their small size, low power consumption, and the reliability of batch fabrications, there are lots of potential applications in engineering [1, 2]. N/MEMS are also widely spread in biological purposes. Drug delivery, chemical analysis, molecular separation, DNA analysis, and microfluidic systems are some of biological usages of N/MEMS [3, 4].Some of N/MEMS devices consist of a parallel-plate capacitor as shown in Figure1. In this figure, the upper nano/microplate is considered as movable electrode which is modelled as a fully clamped nano/microplate and the other one is stationary electrode. Both mechanical and electrical domains play crucial role in N/MEMS. The nano/microplate as mechanical part of the N/MEMS deforms under the electrical loading. This deformation changes the capacitance and also the attracting electrostatic force which acts on the movable electrode. By increasing the input voltage, the microplate deflection may be increased. The induced electrostatic force has an upper limit that overcomes the elastic restoring force and causes the sudden collapse conditions in the nano/microstructure. This unstable behaviour of electrically actuated nano/microstructures is called pull-in instability and the associated voltage and associated displacement are known as pull-in voltage and pull-in displacement, respectively.Figure 1
A schematic view of an electrostatically actuated nano/microplate.Decreasing the dimension of structures from microscale to nanoscale makes the intermolecular surface forces more considerable. The Casimir and the van der Waals attractions are the most important forces which act at the dimensions of N/MEMS [5]. Emission of retarded electromagnetic waves in the separations less than few micrometers makes an intermolecular force, known as the Casimir force [6]. The van der Waals force acts in the shorter distance range compared to the Casimir force which is more considerable in nanoscale [6]. The Casimir force plays a vital role in N/MEMS design [7] inasmuch as neglecting this effect in designing N/MEMS may lead to the collapse condition of the structure without applying any external voltage [8]. Therefore, investigating the Casimir effect on N/MEMS is inescapable.The Casimir effect was mentioned and presented theoretically in 1948 for the first time by the Dutch physicist Hendrik Brugt Gerhard Casimir (1909–2000), who investigates the attraction of two parallel plates in the vacuum [9]. Several researches investigated the role of the Casimir effect on the pull-in phenomenon in N/MEMS [10–13]. Lin and Zhao studied the Casimir effect on the static pull-in and dynamic behaviour of the nanoscale electrostatic actuator using a lumped model [14]. Ramezani et al. investigated pull-in instability of cantilever type NEMS considering both electrostatic and intermolecular surface forces using a distributed parameter model [15]. They used Green’s function to transform the nonlinear differential equation of the model into the integral form. Koochi et al. investigated the influence of surface effects including residual surface stress and surface elasticity as well as the effect of Casimir attraction on the pull-in instability of a cantilever nanoactuator [16]. Beni et al. proposed a model to consider the influence of surface effect and molecular forces on pull-in voltage of rotational nano/micromirror using two degrees of freedom model [17]. They also investigated the effect of the Casimir attraction on the torsion/bending coupled pull-in instability of an electrostatic nanoactuator [18]. Farrokhabadi et al. modelled the Casimir force-induced instability in freestanding nanowires with circular cross-section analytically [19]. Zhao et al. suggested a reduced-order model for electrically actuated rectangular microplates. They found the linear mode-shapes numerically using the hierarchical finite element method [20]. Batra et al. presented a review of modelling electrostatically actuated MEMS [21]. They also considered the Casimir force to develop reduced-order models for prestressed clamped rectangular and circular electrostatically actuated microplates [22]. They used their reduced-order models to study the static pull-in instability as well as free vibrations of the plate about its predeformed configuration [22]. They found that using only one mode for transverse deflection in the Galerkin procedure can provide high accurate results [22].One of the most important issues to model the fully clamped nano/microplates is that the governing equation of the model is strongly nonlinear, so there is no analytical solution for such systems in the literature. In addition, the available numerical solutions suffered from long run time. The main idea of present paper is to represent an accurate solution for such problems which rapidly converges and reduces the time and the cost of the calculation. To this end, the method of extended Kantorovich is employed to extract the accurate mode-shapes of the fully clamped microplate. Using the extracted mode-shapes, the nonlinear governing boundary value equation of equilibrium is reduced to a nonlinear algebraic equation. The resulting nonlinear algebraic equation is also solved through an iterative procedure developed based on Taylor’s series expansion. It is shown that the results of present study are in excellent agreement with the finite element simulations carried out using the COMSOL Multiphysics as well as those observed experimentally in the literature.


## 2. Problem Formulations

As shown in Figure1, a flexible clamped nano/microplate is considered with length a, width b, thickness h, and density ρ under the combined effects of electrostatic and Casimir pressures. The governing equation of motion for deflection of thin nano/microplates can be written as [23]
(1)ρh∂2W∂t2+D∂4W∂X4+2∂4W∂X2∂Y2+∂4W∂Y4=FElectrostaticW,V+FCasimirW,
where W and D are the deflection and flexural rigidity of the nano/microplate, respectively. Also FElectrostatic and FCasimir, respectively, denote the distributed electrostatic and the Casimir pressures. These quantities can be defined as [22]
(2)D=Eh312(1-v2),(3)FElectrostatic=ε0εrV22G-WX,Y2,(4)FCasimir=π2ħc240G-WX,Y4,
where v is Poisson’s ratio, V is the applied DC voltage, G is the initial gap between two electrodes, εr is the relative permittivity, ε0 is the permittivity of the air, ħ=1.055×10-34 is Planck’s constant divided by 2π, and c=2.998×108 is the speed of light in vacuum. Neglecting the inertia term in (1) leads to the governing static equilibrium equation as follows:
(5)D∂4W∂X4+2∂4W∂X2∂Y2+∂4W∂Y4=ε0V22G-WX,Y2+π2ħc240G-WX,Y4.
To nondimensionalize the governing equation of equilibrium, the nondimensional parameters x=X/a,  y=Y/b, and w=W/G are considered. Substituting these quantities into (5) leads to
(6)∂4w∂x4+2α2∂4w∂x2∂y2+α4∂4w∂y4=β1-w2+γ1-w4,
where α, β, and γ are the aspect ratio of the plate (α=a/b), the nondimensional parameter of voltage, and the Casimir force parameter, respectively. β and γ can also be obtained as
(7)β=a4ε0εrV22DG3,(8)γ=a4π2ħc240DG5.The corresponding nondimensional boundary conditions for fully clamped microplate are(9a)wx,0.5=wx,-0.5=0,(9b)w0.5,y=w-0.5,y=0,(9c)∂w∂xx=0.5,-0.5=0,(9d)∂w∂yy=0.5,-0.5=0.


## 3. Solution Procedure

Due to the high nonlinearity involved in (6), a closed-form solution for this equation cannot be found. Hence, an approximate solution will be developed here through the Galerkin weighted residual method. Based on this procedure, the nano/microplate deflection can be expressed as a linear combination of a complete set of linearly independent basis functions [24]. It is noted that these functions must satisfy all kinematic boundary conditions [24]. Therefore, linear and undamped mode-shapes of the undeformed nano/microplate can be used as these basis functions. It is proved that using only the first mode for approximating the transverse deflection in static, dynamic, vibrational, and pull-in analysis of electrically actuated nano/microbeams and plates may be very accurate [2, 4, 22, 25, 26]. Hence, the deflection of the plate based on one-mode solution can be expressed as
(10)wx,y=uφ11x,y,
where u is an unknown parameter determined through the Galerkin procedure and φ11x,y is the first linear and undamped mode-shape of the undeformed nano/microplate which will be determined in the next section through the extended Kantorovich method (EKM).According to the Galerkin procedure, (6) is multiplied by φ11x,y, (10) is substituted into the resulting equation, and the outcome is integrated over the dimensionless domain, obtaining
(11)u∫-1/21/2∫-1/21/2∂4φ11∂x4+2α2∂4φ11∂x2∂y2+α4∂4φ11∂y4φ11dxdy=∫-1/21/2∫-1/21/2βφ111-uφ112+γφ111-uφ114dxdy.
Equation (11) is a nonlinear algebraic equation which relates the nondimensional input voltage (i.e., β) and the Casimir parameter (i.e., γ) to the microplate deflection. To solve this equation, the iterative method presented by Tahani and Askari [4] is used and (11) is linearized by expanding its nonlinear terms about an unknown static deflection u0 using zeroth-order Tailor’s series expansion as
(12)u∫-1/21/2∫-1/21/2∂4φ11∂x4+2α2∂4φ11∂x2∂y2+α4∂4φ11∂y4φ11dxdy=∫-1/21/2∫-1/21/2βφ111-u0φ112+γφ111-u0φ114dxdy.
Therefore, the solution takes the form
(13)u=∫-1/21/2∫-1/21/2βφ111-u0φ112+γφ111-u0φ114dxdy·∫-1/21/2∫-1/21/2∂4φ11∂x4+2α2∂4φ11∂x2∂y2+α4∂4φ11∂y4φ11dxdy-1.
Hence, the solution can be found through the following iterative formula:
(14)uN+1=∫-1/21/2∫-1/21/2βφ111-uNφ112+γφ111-uNφ114dxdy·∫-1/21/2∫-1/21/2∂4φ11∂x4+2α2∂4φ11∂x2∂y2+α4∂4φ11∂y4φ11dxdy-1.At the first step, the value ofuN is set to zero and the iterative procedure is continued till the convergence is achieved or pull-in happened. The convergence criterion is defined as
(15)uN+1-uNuN≤10-6,
and the pull-in will happen if
(16)wmid-point=φ110,0uN≥1.
It is to be noted that although we use wmid-point≥1 for pull-in condition in our solution, this instability can be observed when the slope of deflection-voltage graph reaches infinity [4].


## 4. Application of the EKM for Obtaining Linear and Undamped Mode-Shapes

The governing eigenvalue equation for the free vibrations of a microplate can easily be obtained as(17)∂4φmn∂x4+2α2∂4φmn∂x2∂y2+α4∂4φmn∂y4-ωmn2φmn=0,
where ωmn represents the natural frequency of the microplate. According to the general procedure of the EKM, the microplate mode-shape should be considered as a multiplication of single term separable functions as
(18)φmnx,y=fxgy,
where fx and gy are unknown functions in terms of x and y to be determined. The dimensionless boundary conditions presented in (9a), (9b), (9c), and (9d) in terms of the separable functions can be written as(19a)f-12=f12=dfdxx=-1/2=dfdxx=1/2=0,(19b)g-12=g12=dgdxx=-1/2=dgdxx=1/2=0.Based on the general procedure of the Galerkin weighted residual method (GWRM) [24], the variational form of (17) can be obtained as
(20)∫-1/21/2∫-1/21/2∂4φmn∂x4+2α2∂4φmn∂x2∂y2+α4∂4φmn∂y4-ωmn2φmnδφmndxdy=0.
According to the iterative nature of the EKM, at the first step, it is assumed that g is a prescribed known function. Hence, the variation of φmn would be just because of the variation of the unknown function fx. Substituting (18) into (20) leads to
(21)∫-1/21/2∫-1/21/2gd4fdx4+2α2d2fdx2d2gdy2+α4fd4gdy4-ωmn2fggdyδfdx=0.
Using the fact that δf is an arbitrary value except at the boundaries and employing the fundamental lemma of variational calculus, we obtain
(22)∫-1/21/2g2dyd4fdx4+2α2∫-1/21/2gd2gdy2dyd2fdx2+α4∫-1/21/2gd4gdy4dy-ωmn2∫-1/21/2g2dyf=0.
Next, using the integration by parts, one can write
(23)d4fdx4-I1d2fdx2+I2-ωmn2f=0,
where(24a)I1=2α2∫-1/21/2dg/dy2dy∫-1/21/2g2dy,(24b)I2=α4∫-1/21/2d2g/dy22dy∫-1/21/2g2dy.At the next step, using the prescribed known functionf obtained from (23), the unknown function g can be determined through the following equation:
(25)d4gdy4-I1′d2gdy2+I2′-ωmn2g=0,
where
(26a)I1′=2∫-1/21/2df/dx2dxα2∫-1/21/2f2dx,(26b)I2′=∫-1/21/2d2f/dx22dxα4∫-1/21/2f2dx.
By choosing an initial guess function gy=4y2-12, the solution of (23) takes the form
(27)fx=C1sinh⁡ϑx+C2cosh⁡ϑx+C3sinυx+C4cos⁡υx,
where(28a)ϑ=I12+ωmn2+I122-I2,(28b)υ=-I12+ωmn2+I122-I2.Imposing the boundary conditions presented in (19a) yields the following eigenvalue, eigenfunction equations:(29a)C1sinh⁡ϑ2-C2cosh⁡ϑ2+C3sinυ2-C4cos⁡υ2=0,(29b)C1sinh⁡ϑ2+C2cosh⁡ϑ2+C3sinυ2+C4cos⁡υ2=0,(29c)ϑC1cosh⁡ϑ2-ϑC2sinh⁡ϑ2+υC3cos⁡υ2+υC4sinυ2=0,(29d)ϑC1cosh⁡ϑ2+ϑC2sinh⁡ϑ2+υC3cos⁡υ2-υC4sinυ2=0.This system of equations can be simplified to(30a)C2cosh⁡ϑ2+C4cos⁡υ2=0,(30b)ϑC2sinh⁡ϑ2-υC4sinυ2=0,(31a)C1sinh⁡ϑ2+C3sinυ2=0,(31b)ϑC1cosh⁡ϑ2+υC3cos⁡υ2=0.It is worth noting that (30a), (30b), (31a), and (31b) are related to the symmetric and antisymmetric mode-shapes, respectively. Equating the determinant of the coefficients of these two systems of equations would lead to the following frequency equations for symmetric and antisymmetric modes, respectively:
(32)υtanυ2=-ϑtanhϑ2,ϑtanυ2=υtanhϑ2.
After finding the values of ϑ and υ, by normalizing the mode-shapes (i.e., setting C4=1 for symmetric modes and C3=1 for antisymmetric modes), the value of C2 and C1 can be obtained using (30a) and (31a), respectively. Next, by continuing the iterative procedure, one can utilize the solution of (23) (i.e., (27)) as a prescribed function and find the unknown function gy using (25). The procedure should be continued till the convergence occurs.


## 5. Finite Element Analysis

In this section, a fully clamped nano/microplate is modelled using commercial FE software COMSOL Multiphysics. The model is designed to investigate pull-in instability of a nano/microplate under the effect of electrostatic field induced by applying DC voltage. Geometrical and material properties of the model are also given in Table1. To find pull-in voltage using COMSOL Multiphysics software, for the first step, the stationary problem should be solved [27]. To do so, a 3D tetrahedral free mesh is employed. Also the electrical field is applied by assigning zero and stationary electric potentials to the movable and fixed electrodes, respectively. The stationary electric potential is also set to vary from 0 V to the value in which pull-in instability occurs. Pull-in voltage would be found by solving an inverse problem in COMSOL [27]. In this procedure, COMSOL changes the applied voltage to satisfy the equation of the inverse problem. In this study, the FE results obtained by COMSOL are utilized to validate the findings of the present approach.Table 1
Specifications of the system for FE simulation.
Specification
Value
Width (µm)
250
Length (µm)
250
Thickness (µm)
3
Young’s modulus (GPa)
169
Poisson’s ratio
0.28
Density (Kg/m3)
2331
Initial gap (µm)
1
Dielectric of air (F/m)
8.85
×
1
0
-
12In order to apply the Casimir effect to the FE model, COMSOL functions are employed. The Casimir force is applied to the bottom surface of the flexible electrode (i.e., the microplate) as a distributed variable boundary load per unit area and its value is obtained through (4). It is to be noted that the Casimir force is a function of deflection, so it is varying in a range that should be set by the user based on the usual range of the pull-in displacements.


## 6. Results and Discussions

### 6.1. Comparison and Validation

To find the number of iteration which must be carried out in the EKM, a convergence study is performed in Table2 for square microplate first mode-shape parameters as well as its fundamental natural frequency. From Table 2, it is apparent that the convergence is so rapid and that in all cases it occurs after the second iteration. It should be noted that the calculated present fundamental natural frequency of the system is in an excellent agreement with ω11=36.1087 obtained by Arenas [28]. To illustrate the accuracy of the present analysis, the first mode-shape of the microplate is also depicted in Figure 2.Table 2
Convergence study for parameters of the first mode-shape as well as dimensionless fundamental natural frequencies of a square microplate.
i
ϑ
i
,
i
-
1
ϑ
i
,
i
υ
i
,
i
-
1
υ
i
,
i
C
2
i
,
i
-
1
C
2
i
,
i
ω
11
i
,
i
-
1
ω
11
i
,
i
1
—
6.52637011
—
4.31201887
—
0.04221306
35.99948133
35.99948133
2
6.52610255
6.52610827
4.31206780
4.31206675
0.04222025
0.04222010
35.99896276
35.99896276
3
6.52610815
6.52610815
4.31206678
4.31206678
0.04222010
0.04222010
35.99896276
35.99896276
4
6.52610815
6.52610815
4.31206678
4.31206678
0.04222010
0.04222010
35.99896276
35.99896276Figure 2
The first mode-shape of a fully clamped square microplate.To validate the accuracy of present pull-in results, variation of dimensionless mid-point deflection versus the nondimensional applied voltage for a square microplate without the effect of Casimir force is depicted in Figure3. The results are compared and validated by available experimental and other existing results in the literature [20, 29]. As one can observe from this figure, the present findings are in excellent agreement with those reported in previous studies [20, 29]. As Figure 3 shows, the slope of mid-point deflection-voltage graph reaches infinity when the input voltage exceeds a critical value called pull-in voltage. It is noteworthy that this fact was used by some researchers to determine the pull-in condition in electrically actuate microbeams [14, 30].Figure 3
Comparison between present findings and available empirical [29] and numerical [20] results for a square microplate.To emphasize the accuracy of the present analysis, FE simulation is also employed in the present paper. Figure4 shows the deformed configuration of a microsystem modelled in COMSOL Multiphysics software. It should be noted that the properties of the model were also presented in Table 1. As it is seen from Figure 4, the pull-in displacement (the maximum displacement occurred just before bechancing pull-in instability) is 41.61% of the initial gap counterpart. To compare the results of the present iterative approach with those of FE simulation, variation of the mid-point deflection versus the applied voltage is depicted in Figure 5. According to the results shown in Figure 5, there exists an excellent agreement between the findings of present iterative approach and those of 3D FE analysis carried out by COMSOL.Figure 4
Deformed configuration of the system with properties given in Table1 at pull-in condition.Figure 5
Comparison between present findings and FE results for a square microplate with properties listed in Table1.


### 6.2. Effect of the Casimir Force

By decreasing the dimensions of electromechanical systems from microscale to nanoscale, the effect of dispersion forces may be considerable. To investigate the effect of Casimir force on the stable and unstable responses of the system, variations of the nondimensional mid-point deflection versus the dimensionless applied voltage for microplates with some different aspect ratios are depicted in Figure6. According to this figure, by increasing the values of Casimir parameter, the mid-point deflection is increased which results in earlier pull-in instability of the system. The values of nondimensional pull-in parameter (i.e., βPI) versus the Casimir parameter of the system (i.e., γ) at some different aspect ratios of the plate are also given in Table 3. As it is seen from this table, decreasing the aspect ratio (i.e., α) reduces the nondimensional parameter of pull-in voltage and the percentage of the reduction increases as the value of γ is increased. Furthermore, according to Table 3, the system may be collapsed without applying the external voltage, if the Casimir parameter of the system reaches its critical value. It is to be noted that this condition is undesirable in designing N/MEMS and should be avoided.Table 3
Nondimensional pull-in voltage versus the Casimir parameter.
γ
β
Pl (α = 1)
β
Pl (α = 0.75)
β
Pl (α = 0.5)
β
Pl (α = 0.25)
0
176.7
113.7
84.6
73.2
20
134.7
73.2
45.3
34.5
40
97.5
38.4
12.3
0
60
63
6.9
0
0
80
30.6
0
0
0
100
0
0
0
0Effect of the Casimir force on the variation of nondimensional mid-point deflection versus the dimensionless pull-in voltage for nano/microplates with different aspect ratios: (a)α=1, (b) α=0.75, (c) α=0.5, and (d) α=0.25.
(a)
(b)
(c)
(d)The variation of critical Casimir parameter versus the variation of the plate’s aspect ratio is depicted in Figure7. According to the results of this figure, the critical value of Casimir parameter increases with an increase of the plate’s aspect ratio. According to (8), γ is a function of geometrical and material properties of the system such as the initial gap between the two electrodes and the dimensions of the movable electrode. Based on this equation, when the initial gap between two electrodes is less than its minimum allowable size or the dimensions of the nano/microplate become greater than their allowable sizes, γ reaches its critical value. In such conditions, the movable electrode collapses without applying the external voltage. Hence, such a system cannot work appropriately. The maximum allowable dimensions of the movable electrode as well as the minimum allowable initial gap can be easily determined using (8) as
(33)amax⁡=240DG5γcrπ2ħc4,bmax⁡=αa,Gmin⁡=a4π2ħc240Dγcr5.
It is noteworthy that the initial gap between the two electrodes and the dimensions of the movable electrode should be controlled to avoid the collapse condition in N/MEMS.Figure 7
Variation of the critical Casimir parameter versus the plate’s aspect ratio.


## 6.1. Comparison and Validation

To find the number of iteration which must be carried out in the EKM, a convergence study is performed in Table2 for square microplate first mode-shape parameters as well as its fundamental natural frequency. From Table 2, it is apparent that the convergence is so rapid and that in all cases it occurs after the second iteration. It should be noted that the calculated present fundamental natural frequency of the system is in an excellent agreement with ω11=36.1087 obtained by Arenas [28]. To illustrate the accuracy of the present analysis, the first mode-shape of the microplate is also depicted in Figure 2.Table 2
Convergence study for parameters of the first mode-shape as well as dimensionless fundamental natural frequencies of a square microplate.
i
ϑ
i
,
i
-
1
ϑ
i
,
i
υ
i
,
i
-
1
υ
i
,
i
C
2
i
,
i
-
1
C
2
i
,
i
ω
11
i
,
i
-
1
ω
11
i
,
i
1
—
6.52637011
—
4.31201887
—
0.04221306
35.99948133
35.99948133
2
6.52610255
6.52610827
4.31206780
4.31206675
0.04222025
0.04222010
35.99896276
35.99896276
3
6.52610815
6.52610815
4.31206678
4.31206678
0.04222010
0.04222010
35.99896276
35.99896276
4
6.52610815
6.52610815
4.31206678
4.31206678
0.04222010
0.04222010
35.99896276
35.99896276Figure 2
The first mode-shape of a fully clamped square microplate.To validate the accuracy of present pull-in results, variation of dimensionless mid-point deflection versus the nondimensional applied voltage for a square microplate without the effect of Casimir force is depicted in Figure3. The results are compared and validated by available experimental and other existing results in the literature [20, 29]. As one can observe from this figure, the present findings are in excellent agreement with those reported in previous studies [20, 29]. As Figure 3 shows, the slope of mid-point deflection-voltage graph reaches infinity when the input voltage exceeds a critical value called pull-in voltage. It is noteworthy that this fact was used by some researchers to determine the pull-in condition in electrically actuate microbeams [14, 30].Figure 3
Comparison between present findings and available empirical [29] and numerical [20] results for a square microplate.To emphasize the accuracy of the present analysis, FE simulation is also employed in the present paper. Figure4 shows the deformed configuration of a microsystem modelled in COMSOL Multiphysics software. It should be noted that the properties of the model were also presented in Table 1. As it is seen from Figure 4, the pull-in displacement (the maximum displacement occurred just before bechancing pull-in instability) is 41.61% of the initial gap counterpart. To compare the results of the present iterative approach with those of FE simulation, variation of the mid-point deflection versus the applied voltage is depicted in Figure 5. According to the results shown in Figure 5, there exists an excellent agreement between the findings of present iterative approach and those of 3D FE analysis carried out by COMSOL.Figure 4
Deformed configuration of the system with properties given in Table1 at pull-in condition.Figure 5
Comparison between present findings and FE results for a square microplate with properties listed in Table1.


## 6.2. Effect of the Casimir Force

By decreasing the dimensions of electromechanical systems from microscale to nanoscale, the effect of dispersion forces may be considerable. To investigate the effect of Casimir force on the stable and unstable responses of the system, variations of the nondimensional mid-point deflection versus the dimensionless applied voltage for microplates with some different aspect ratios are depicted in Figure6. According to this figure, by increasing the values of Casimir parameter, the mid-point deflection is increased which results in earlier pull-in instability of the system. The values of nondimensional pull-in parameter (i.e., βPI) versus the Casimir parameter of the system (i.e., γ) at some different aspect ratios of the plate are also given in Table 3. As it is seen from this table, decreasing the aspect ratio (i.e., α) reduces the nondimensional parameter of pull-in voltage and the percentage of the reduction increases as the value of γ is increased. Furthermore, according to Table 3, the system may be collapsed without applying the external voltage, if the Casimir parameter of the system reaches its critical value. It is to be noted that this condition is undesirable in designing N/MEMS and should be avoided.Table 3
Nondimensional pull-in voltage versus the Casimir parameter.
γ
β
Pl (α = 1)
β
Pl (α = 0.75)
β
Pl (α = 0.5)
β
Pl (α = 0.25)
0
176.7
113.7
84.6
73.2
20
134.7
73.2
45.3
34.5
40
97.5
38.4
12.3
0
60
63
6.9
0
0
80
30.6
0
0
0
100
0
0
0
0Effect of the Casimir force on the variation of nondimensional mid-point deflection versus the dimensionless pull-in voltage for nano/microplates with different aspect ratios: (a)α=1, (b) α=0.75, (c) α=0.5, and (d) α=0.25.
(a)
(b)
(c)
(d)The variation of critical Casimir parameter versus the variation of the plate’s aspect ratio is depicted in Figure7. According to the results of this figure, the critical value of Casimir parameter increases with an increase of the plate’s aspect ratio. According to (8), γ is a function of geometrical and material properties of the system such as the initial gap between the two electrodes and the dimensions of the movable electrode. Based on this equation, when the initial gap between two electrodes is less than its minimum allowable size or the dimensions of the nano/microplate become greater than their allowable sizes, γ reaches its critical value. In such conditions, the movable electrode collapses without applying the external voltage. Hence, such a system cannot work appropriately. The maximum allowable dimensions of the movable electrode as well as the minimum allowable initial gap can be easily determined using (8) as
(33)amax⁡=240DG5γcrπ2ħc4,bmax⁡=αa,Gmin⁡=a4π2ħc240Dγcr5.
It is noteworthy that the initial gap between the two electrodes and the dimensions of the movable electrode should be controlled to avoid the collapse condition in N/MEMS.Figure 7
Variation of the critical Casimir parameter versus the plate’s aspect ratio.


## 7. Conclusions

The main idea of this study was to present an accurate novel method to investigate the behaviour of the parallel-plate nano/microcapacitors. To this end, the extended Kantorovich method tighter with the iterative procedure proposed by Tahani and Askari was used for obtaining both stable and unstable behaviors of the system. The results were compared and validated by previous theoretical and experimental findings as well as those of FE simulation carried out by COMSOL Multiphysics software. It was found that accounting for the Casimir effect decreases the pull-in voltage of the system especially for the cases with large Casimir parameter. Also, it was observed that the movable electrode may be collapsed without applying the external voltage for systems in which the Casimir parameter reaches its critical value. The latest finding was also utilized to obtain the critical dimensions of the movable electrode as well as the minimum allowable initial gap of the system.



---

*Source: 328742-2015-02-09.xml*