Next Article in Journal
Non-Linear Inner Structure of Topological Vector Spaces
Next Article in Special Issue
Modelling an Industrial Robot and Its Impact on Productivity
Previous Article in Journal
Existence and Concentration Behavior of Solutions of the Critical Schrödinger–Poisson Equation in R3
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A 2D Membrane MEMS Device Model with Fringing Field: Curvature-Dependent Electrostatic Field and Optimal Control

1
Dipartimento di Ingegneria Industriale e dell’Informazione, University of Pavia, Via A. Ferrata 5, 27100 Pavia, Italy
2
Dipartimento di Ingegneria dell’Informazione Infrastrutture Energia Sostenibile, “Mediterranea” University, Via Graziella Feo di Vito, 89124 Reggio Calabria, Italy
3
Dipartimento di Ingegneria Civile Energia Ambiente e Materiali, “Mediterranea” University, Via Graziella Feo di Vito, 89124 Reggio Calabria, Italy
*
Author to whom correspondence should be addressed.
Mathematics 2021, 9(5), 465; https://doi.org/10.3390/math9050465
Submission received: 10 February 2021 / Revised: 19 February 2021 / Accepted: 22 February 2021 / Published: 25 February 2021
(This article belongs to the Special Issue Mathematical Problems in Mechanical Engineering)

Abstract

:
An important problem in membrane micro-electric-mechanical-system (MEMS) modeling is the fringing-field phenomenon, of which the main effect consists of force-line deformation of electrostatic field E near the edges of the plates, producing the anomalous deformation of the membrane when external voltage V is applied. In the framework of a 2D circular membrane MEMS, representing the fringing-field effect depending on | u | 2 with the u profile of the membrane, and since strong E produces strong deformation of the membrane, we consider | E | proportional to the mean curvature of the membrane, obtaining a new nonlinear second-order differential model without explicit singularities. In this paper, the main purpose was the analytical study of this model, obtaining an algebraic condition ensuring the existence of at least one solution for it that depends on both the electromechanical properties of the material constituting the membrane and the positive parameter δ that weighs the terms | u | 2 . However, even if the the study of the model did not ensure the uniqueness of the solution, it made it possible to achieve the goal of finding a stable equilibrium position. Moreover, a range of admissible values of V were obtained in order, on the one hand, to win the mechanical inertia of the membrane and, on the other hand, to ensure that the membrane did not touch the upper disk of the device. Lastly, some optimal control conditions based on the variation of potential energy are presented and discussed.

1. Introduction

The remarkable development of embedded technologies in recent years is, in large part, due to the small size of the devices used that manage the link between the physical nature of the problem and the logic of the machine language [1,2]. In this context, interest in micro-electro-mechanical system (MEMS) devices has been very high since they represent a good approximation of the human–machine interface [3]. MEMS are devices that have been recognized as one of the most promising of the 21st century, capable of revolutionizing both the industrial and the consumer-product worlds. MEMS are devices integrated in miniaturized form on the same substrate as that of semiconductor material (silicon), combining the electrical properties of the integrated semiconductor with optomechanical properties [1]. These are “intelligent” systems that combine electronic, optical, biological, chemical, and mechanical fluid management functions in a very small space, integrating sensor and actuator technology, and the most diverse process-management functions [1,3]. MEMS technologies are adopted in the most varied application fields, many of which are based on microscopic mirrors [1,2] or oscillating lenses in a single or array version that are used to create complex optoelectronic devices [2]. In microwave electronics, MEMS devices are often used as single switches for complex applications such as matching networks, resonant filters, power-supply networks for array antennas, and generally reconfigurable systems [2]. MEMS are also used for new solutions in chemistry and bioengineering [2]. Since the first batch device was produced [4], technological development has strongly influenced the production of physicomathematical models, describing increasingly complex multiphysics [5,6,7]. However, such models, although theoretically valid, rarely provide explicit solutions, so that the conditions of existence and uniqueness of the solution need to be obtained with reasonable computational costs [8,9,10]. Furthermore, if the solution is not analytically obtainable, one can rely on numerical techniques that provide approximate solutions that, if they satisfy the aforementioned conditions of existence and uniqueness [11,12,13], do not represent ghost solutions [14,15,16,17,18]. The scientific community is currently working hard both on the analysis and synthesis of multiphysical models, and on technology transfer [19,20,21,22,23,24]. In such contexts, it is preferred to study MEMS devices equipped with symmetries in order to obtain models that can be studied more easily, both mathematically and physically [12]. Thus, for application reasons, the authors focus their attention on a 2D circular membrane MEMS device used in many industrial and biomedical applications [1,12,25,26,27,28]. Particularly, the authors consider a bounded circular smooth domain Ω , where a 2D circular membrane MEMS device is represented by two parallel disks placed d apart [12]. Furthermore, an electrically conductive elastic membrane (which has the same size as the disks at rest) was clumped on the edge of the lower disk (which acted as a support to the membrane) deforming toward the upper disk when an external electric voltage V is applied between the disks. This deformation, in this work, is described by the profile of the membrane, u ( r ) , r [ 0 , R ] , with r being the radial coordinate [1,5,6,11]. The model took the following well-known form [1,5,6,11]:
Δ u ( r ) = λ 2 ( 1 u ( r ) ) 2 u ( R ) = 0 , u ( 0 ) = 0 , 0 < u ( r ) < d
in which λ 2 is a parameter depending on V. However, electrostatically, if d R , the fringing-field phenomenon occurs [29,30,31]. In other words, the lines of force of electric field E away from the edges of the device appear parallel and uniform, while as one approaches the edges of the disks, they become a curve [29,30,31]. Thus, with these premises, in Model (1), λ 2 becomes λ 2 ( 1 + δ | u ( r ) | 2 ) where term δ | u ( r ) | 2 , with δ R + , considers the fringing-field effects (when δ = 0 , no fringing field phenomenon is observed). Since the device presents radial symmetry with respect to vertical axes r = 0 , in Δ u ( r ) , only the radial component is taken into account, and u ( r ) = u ( r ) . Upon deformation of the membrane, the electrostatic capacitance of the device varies, as the distance between the membrane and upper disk is variable. Furthermore, λ 2 ( 1 u ( r ) ) 2 | E | 2 , so that λ 2 ( 1 u ( r ) ) 2 = θ | E | 2 , θ R + [11]. Moreover, physically, E on the membrane is locally orthogonal to the straight-line tangent to the profile of the membrane. Thus, it makes sense to consider | E | proportional to mean curvature K ( r , u ( r ) ) of the membrane [12]. So, as detailed in Section 3, the model becomes the following 2D second-order semilinear elliptic model:
u ( r ) = 1 r u ( r ) 4 ( 1 u ( r ) d * ) 2 θ λ 2 ( 1 + δ | u ( r ) | 2 ) u ( R ) = 0 , u ( 0 ) = 0 , 0 < u ( r ) < d ,
in which d * is the critical security distance ensuring that the membrane does not touch the upper disk. In recent years, there has been an interesting surge in the application of machine learning and statistical frameworks to solve similar problems as those in this paper, as highlighted in the studies shown in [32,33,34]. Model (2) is very difficult to solve and, in any cases, does not give explicit solutions, so that we focus on if any conditions ensuring both the existence and uniqueness of the solution are obtainable. In particular, considering Model (2) in a more general formulation, an interesting result of the existence of at least one solution is obtained if an algebraic condition, depending on the electromechanical properties of the material constituting the membrane and fringing-field parameter δ , is satisfied. However, although the uniqueness of the solution is not ensured, the unique admitted equilibrium position results in being stable. Moreover, the range of admissible values of V, to both win the mechanical inertia of the membrane ( V m i n ) i n e r t i a and to avoid the membrane touching the upper disk ( V m a x ) p e r m i s s i b l e , is achieved and discussed. Lastly, some optimal control conditions, formulated in terms of variation of potential energy, are discussed. In order to improve the readability of the paper, Figure 1 shows a schematic flowchart of the entire work.

2. Circular-Membrane MEMS Devices: Overview

The MEMS device studied in this paper consisted of two parallel metal disks of which the radius is R, placed at a distance d between them [12]. A circular membrane of radius R in its rest condition was anchored to the edges of the lower disk (which acted as a support). The membrane, under the effect of an external V, deformed towards the upper plate without touching it. Then, E in the device generated electrostatic pressure p e l = 0.5 ϵ 0 | E | 2 ( ϵ 0 , permittivity of free space) that deflected the membrane. Electrostatically, when the membrane deformed, E , which depended on the distance between membrane and upper disk, was locally orthogonal to the line tangent to the membrane at the point. Furthermore, the electrostatic capacitance of the device, C e l , was also variable, as the distance between membrane and upper disk locally varied. The higher | E | was, the greater the curvature of the membrane, so that | E | was proportional to curvature K of the membrane.

2.1. Circular-Membrane MEMS as Transducer

For our purposes, we exploited some similarities with the model of a MEMS circular plate transducer subjected to mechanical pressure p.
Remark 1.
Usually, the terms “sensor” and “transducer” are used interchangeably. However, a distinction must be made, as a sensor is a sensitive element that converts an input into a physical output, which is acquired into an electrical signal. On the other hand, a transducer consists of an input interface, a sensor, and an output interface. In this paper, we refer to a MEMS transducer in order to consider both interfaces.
In this section, we report some important results concerning the mechanics, as already discussed in [12], of when a metal plate is subjected to mechanical pressure p, deforms, and its deflection u satisfies the known equation [1,35]
ρ h u t t T Δ u + D Δ 2 u = 0
with ρ being the density of the material constituting the plate, with thickness h, T, and [35]
D = Y h 3 12 ( 1 ν 2 ) ,
with Y (Young’s modulus) and ν (Poisson ratio) being the mechanical tension and flexural stiffness of the plate, respectively. Moreover, if the plate is circular, u only depends on radial coordinate r (assuming radial symmetry). With these premises, and in steady-state conditions, u is valuable as [1,35]
u ( r ) = R 4 64 D 1 r R 2 2 p
with 0 r R , for z-directed u [35]. Of course, if r = 0 , displacement at the center of the plate u 0 becomes [1,35]
u 0 = R 4 64 D p
so that (5) becomes [1]
u ( r ) = u 0 1 r R 2 2 p .
Then, unlike the actuator, the device works like a transducer. In fact, as verified in [12], p generates u ( r ) , so that C e l becomes [1,35]
C e l ( u 0 ) = 0 R 2 ϵ 0 π r d 1 u ( r ) d d r , if | u 0 | d .
Both h and D are limited values. Then, u ( r ) becomes inconspicuous, so that the distance between the two plates remains d. Using the Taylor series (up to the third term) and considering the electrostatic capacitance at equilibrium, C 0 = ϵ 0 π R 2 d for p = 0 , Equation (8) becomes [35]
C e l ( u 0 ) C 0 1 + u 0 3 d + u 0 2 5 d 2 ,
which can be exploited to achieve the coenergy of the system, the charge in the membrane, and the electrical force [11]. Furthermore, in [12], C e l ( u 0 ) , as computed in (9), is a nonlinear function of u 0 . Moreover, d C e l ( u 0 ) d u 0 p by electrostatic force f e l . In other words, d C e l ( u 0 ) d u 0 2 f e l V 2 , and, if | u 0 | d
| E ( r ) | V d u 0 1 r R 2 2 .
All involved physical quantities depend on d because the circular plate had a significant value of D; then, u ( r ) was extremely limited, so that any dependencies on d u ( r ) could be replaced by the simpler dependence on d. If, instead of the deformable plate, a membrane was considered (see [12]), h is negligible, so that D, as formulated in (4) significantly decreases compared to the case in which a deformable plate is present. Considering that the lower the D value, the more flexible the membrane, u 0 became higher with the risk that the membrane touched the upper plate. In this, it is imperative to consider d u ( r ) in the denominator of (9). By subjecting the membrane to p, u ( r ) becomes [1]
u ( r ) = u 0 1 r R 2
with
u 0 = p R 2 4 T .
In this case, f e l became [1,35]
f e l = 0.5 ϵ 0 π R 2 V 2 ( d u ( r ) ) 2
from which
p e l f e l π R 2 = 0.5 ϵ 0 V 2 ( d u ( r ) ) 2 .
Remark 2.
In the calculation of f e l and p e l , the surface of the membrane was approximated as π R 2 , even in the presence of membrane deformation. This approximation is justifiable because d R , and the surface of the deformed membrane could then be approximated to the surface of the membrane in the rest position.

2.2. p and p e l : An Interesting Relationship

The relationship between p and p e l arises from, when V is applied, E is generated inside the device, also generating p e l , which deforms the membrane, obtaining u ( r ) and making it clear that there is a link between p and p e l . Let us first observe that, from (12), u 0 depends on p and, setting
k 1 = R 2 4 T ,
we can write:
u 0 = p R 2 4 T = k 1 p .
If there are no further stresses, p originates exclusively from p e l , which, in turn, is generated by | E | . Then, the following chain of equalities makes sense:
u 0 = k 1 p = k 1 k 2 p e l = k p e l ,
with both k 2 and k being constant.
Remark 3.
Considering (14), u 0 becomes
u 0 = k p e l = k ϵ 0 V 2 2 ( d u ( r ) ) 2
where d u ( r ) is locally the distance between membrane and upper disk. The profile of the membrane must also absolutely not touch the upper disk. Then, once the membrane is deformed,
u ( r ) d d * ,
must occur, so that
1 ( d u ( r ) ) 2 1 d * 2 .
Graphical details are shown in Figure 2. Thus, Equation (11), considering Equation (18) and Remark 3, can be written as
u ( r ) u ¯ ( r ) = k ϵ 0 V 2 2 d * 2 1 r R 2 .
The relationship between p and p e l constitutes a dual transducer–actuator model. So, the behavior of the transducer helps us to understand the function of an actuator and vice versa.
Remark 4.
Equation (21) is a very important relationship because it represents an upper solution that is useful to exploit a lemma known in the literature to obtain an algebraic condition governing the existence of the solution for (2).

3. Problem Formulation

As introduced above, Model (1), in the presence of a fringing field, becomes
Δ u ( r ) = λ 2 ( 1 + δ | u ( r ) | 2 ) ( 1 u ( r ) ) 2 u ( R ) = 0 , u ( 0 ) = 0 , 0 < u ( r ) < d ,
highlighting a radial symmetry with respect to vertical axes r = 0 , Δ u ( r ) = 1 r u ( r ) + u ( r ) . Moreover, with u ( r ) = u ( r ) , Equation (22) can be written as
1 r u ( r ) + u ( r ) = λ 2 ( 1 + δ | u ( r ) | 2 ) ( 1 u ( r ) ) 2 u ( R ) = 0 , u ( 0 ) = 0 , 0 < u ( r ) < d .
Furthermore, when the membrane deformed, the electrostatic capacitance of the device varied, as the distance between membrane and upper disk was variable. Moreover, in (23), λ 2 ( 1 u ( r ) ) 2 | E | 2 , so that it makes sense to write λ 2 ( 1 u ( r ) ) 2 = θ | E | 2 , θ R + [11]. Thus, Equation (23) becomes
u ( r ) + 1 r u ( r ) = θ | E | 2 ( 1 + δ | u ( r ) | 2 ) u ( R ) = 0 , u ( 0 ) = 0 θ R + , 0 < u ( r ) < d .
In addition, E on the membrane was physically locally orthogonal to the straight-line tangent to the profile of the membrane, so that | E | could be considered proportional to mean curvature K ( r , u ( r ) ) of the membrane [12]
K ( r , u ( r ) ) = 1 2 u ( r ) + 1 r u ( r )
The larger | E | is, the more the membrane deforms. Then, it makes sense to assume that | E | K ( r , u ( r ) ) (as detailed in (25)). Thus,
| E | = μ ( r , u ( r ) , λ ) K ( r , u ( r ) )
where μ ( r , u ( r ) , λ ) , the function of proportionality, can be written as [11,12]
μ ( r , u ( r ) , λ ) = λ ( 1 u ( r ) d * )
with μ ( r , u ( r ) ) C 0 ( A ) and A = [ R , R ] × [ 0 , 1 ) , and, in the radial symmetry framework, K ( r , u ( r ) ) represents the mean curvature, so that (26), considering both (27) and (25), can be written as
| E | 2 = 1 4 λ 2 ( 1 u ( r ) d * ) 2 u ( r ) + 1 r u ( r ) 2
and Model (24) becomes
u ( r ) + 1 r u ( r ) = θ λ 2 4 ( 1 u ( r ) d * ) 2 u ( r ) + 1 r u ( r ) 2 ( 1 + δ | u ( r ) | 2 ) u ( R ) = 0 , u ( 0 ) = 0 θ R + , 0 < u ( r ) < d .
Therefore, from (29), Equation (2) follows.
Remark 5.
Equation (2) was achieved from (29) because it was easy to prove that u ( r ) + 1 r u ( r ) 0 . For details, see [12].

4. General Problem Formulation

Model (2) is a special case of the following general problem:
u ( r ) + F ( r , u ( r ) , u ( r ) ) = 0 u ( b ) = B , u ( a ) = m ,
in which F C 0 ( ( a , b ] × R × R ) and B , m R . In fact, setting
F ( r , u ( r ) , u ( r ) ) = 1 r u ( r ) + 4 ( 1 u ( r ) d * ) 2 θ λ 2 ( 1 + δ | u ( r ) | 2 ) ,
B = m = 0
and
b = R , a = 0
we obtain (2). Now, we obtain conditions ensuring both the existence and uniqueness of the solution for (2). We introduce two preliminary lemmas.

5. Preliminary Lemmas

For further details, see [36] ). Particularly, we exploited Lemma 1 to prove that (2) admits at least one solution. On the other hand, Lemma 2 was used to prove that the solution for (2) is not unique.
Lemma 1.
Let us consider Model (30) and two twice continuously differentiable functions u 1 ( r ) and u 2 ( r ) , such that
u 1 ( r ) < u 2 ( r ) , r ( a , b ) .
Let us also consider
u 1 ( r ) + F ( r , u 1 ( r ) , u 1 ( r ) ) > 0
and
u 2 ( r ) + F ( r , u 2 ( r ) , u 2 ( r ) ) < 0
for r ( a , b ) . Furthermore, let F ( r , y ( r ) , y ( r ) ) be a continuous function satisfying the following Lipschitz condition:
K 1 ( r ) ( u ( r ) v ( r ) ) + L 2 ( r ) ( u ( r ) v ( r ) ) F ( r , u ( r ) , u ( r ) ) F ( r , v ( r ) , v ( r ) ) K 2 ( r ) ( u ( r ) v ( r ) ) + L 1 ( r ) ( u ( r ) v ( r ) )
in U × ( , + ) in which
U = { ( r , u ) : a < r < b and u 1 ( r ) u ( r ) u 2 ( r ) }
and K 1 ( r ) , K 2 ( r ) , L 1 ( r ) and L 2 ( r ) are continuous functions in ( a , b ] . If u 1 ( a ) u 2 ( a ) , with u 1 ( b ) = u 2 ( b ) = B , Equation (30) has at least one solution, u ( r ) , such that u 1 ( r ) u ( r ) u 2 ( r ) in [ a , b ] .
Lemma 2.
Let us suppose that the conditions of Lemma 1 are satisfied, and u 1 ( r ) and u 2 ( r ) satisfy the boundary conditions. Let us also suppose that
u ( r ) + K 2 ( r ) u ( r ) + L 1 ( r ) u ( r ) = 0
has no nontrivial solution that satisfies zero-boundary conditions on any subinterval of [ a , b ] ; then, Equation (30) has only one solution, u ( r ) , such that u 1 ( r ) u ( r ) u 2 ( r ) .

6. An Interesting Result of the Existence of at Least One Solution

In this section, we prove the existence of at least one solution for (2), referring to Section 7, the discussion on uniqueness.
Theorem 1.
Let us consider Problem (2). Let us also consider u 1 ( r ) and u 2 ( r ) as two functions twice continuously differentiable and both defined on [ 0 , R ] , such that u 1 ( r ) < u 2 ( r ) . Moreover, let us suppose that
u 1 ( r ) + 1 r u 1 ( r ) + 4 ( 1 u 1 ( r ) d * ) 2 θ λ 2 ( 1 + δ | u 1 ( r ) | 2 ) > 0
and
u 2 ( r ) + 1 r u 2 ( r ) + 4 ( 1 u 2 ( r ) d * ) 2 θ λ 2 ( 1 + δ | u 2 ( r ) | 2 ) < 0
for r ( 0 , R ) . Furthermore, if 1 r u ( r ) + 4 ( 1 u ( r ) d * ) 2 θ λ 2 ( 1 + δ | u ( r ) | 2 ) is a continuous function (obviously, except for r = 0 ), which satisfies the Lipschitz condition in U × ( , + ) , with U = { ( r , u ) : 0 < r < R and u 1 ( r ) u ( r ) u 2 ( r ) } , and if u 1 ( 0 ) u 2 ( 0 ) , with u 1 ( R ) = u 2 ( R ) = 0 with
θ λ 2 > 2 d * 2 R 2 k ϵ 0 V 2 1 + δ k ϵ 0 V 2 r d * 2 R 2 2 ;
thus, Problem (2) admits at least one solution.
Proof of Theorem 1.
To prove this theorem, we need to exploit Lemma 1. Particularly, we assume that
u 1 ( r ) = 0 , r [ 0 , R ]
and
u 2 ( r ) = u ¯ ( r ) = k ϵ 0 V 2 2 d * 2 1 r R 2 .
as derived from (21). Figure 2 displays both u 1 ( r ) and u 2 ( r ) , a possible recovery of the membrane. By construction, it is clear that u 1 ( r ) < u 2 ( r ) ; moreover, they are twice continuously differentiable functions. We easily observe that u 1 ( r ) = u 1 ( r ) = 0 , so that
u 1 ( r ) + 1 r u 1 ( r ) + 4 ( 1 u 1 ( r ) d * ) 2 θ λ 2 ( 1 + δ | u 1 ( r ) | 2 ) = ( 1 d * ) 2 θ λ 2 > 0 .
Thus, (40) is verified. Moreover,
u 2 ( r ) = k ϵ 0 V 2 r d * 2 R 2
and
u 2 ( r ) = k ϵ 0 V 2 d * 2 R 2
so that we need to verify (41), which becomes
u 2 ( r ) + 1 r u 2 ( r ) + 4 ( 1 u 2 ( r ) d * ) 2 θ λ 2 ( 1 + δ | u 2 ( r ) | 2 ) = 2 k ϵ 0 V 2 d * 2 R 2 + 4 1 k ϵ 0 V 2 2 d * 2 1 r R 2 d * 2 θ λ 2 1 + δ k ϵ 0 V 2 r d * 2 R 2 2 < 0
that is
4 1 k ϵ 0 V 2 2 d * 2 1 r R 2 d * 2 θ λ 2 1 + δ k ϵ 0 V 2 r d * 2 R 2 2 < 2 k ϵ 0 V 2 d * 2 R 2 .
For this purpose, since (49) 0 1 k ϵ 0 V 2 2 d * 2 1 r R 2 d * 2 < 1 , to verify (41) and then (48), it is sufficient to impose
1 θ λ 2 1 + δ k ϵ 0 V 2 r d * 2 R 2 2 < k ϵ 0 V 2 2 d * 2 R 2 ;
moreover, from (43), (44) and (46), we can easily verify u 1 ( 0 ) u 2 ( 0 ) and u 1 ( R ) = u 2 ( R ) = 0 .
Lastly, as Lemma 1 requires, we need to prove that F given by (31) satisfies the Lipschitz condition (37). Thus, starting from (31), ( u ( r ) , v ( r ) U :
F ( r , u ( r ) , u ( r ) ) F ( r , v ( r ) , v ( r ) ) = 1 r u ( r ) + 4 ( 1 u ( r ) d * ) 2 θ λ 2 ( 1 + δ | u ( r ) | 2 ) 1 r v ( r ) 4 ( 1 v ( r ) d * ) 2 θ λ 2 ( 1 + δ | v ( r ) | 2 ) = 1 r ( u ( r ) v ( r ) ) + 1 θ λ 2 4 ( 1 u ( r ) d * ) 2 1 + δ | u ( r ) | 2 4 ( 1 v ( r ) d * ) 2 1 + δ | v ( r ) | 2 = 1 r ( u ( r ) v ( r ) ) + 4 θ λ 2 ( ( 1 u ( r ) d * ) 2 ( 1 + δ | u ( r ) | 2 δ | u ( r ) | 2 ) 1 + δ | u ( r ) | 2 ( 1 v ( r ) d * ) 2 ( 1 + δ | v ( r ) | 2 δ | v ( r ) | 2 ) 1 + δ | v ( r ) | 2 ) = 1 r ( u ( r ) v ( r ) ) + 4 θ λ 2 ( ( 1 u ( r ) d * ) 2 1 δ | u ( r ) | 2 1 + δ | u ( r ) | 2 ( 1 v ( r ) d * ) 2 1 δ | v ( r ) | 2 1 + δ | v ( r ) | 2 ) = 1 r ( u ( r ) v ( r ) ) + 4 θ λ 2 ( ( 1 u ( r ) d * ) 2 ( 1 v ( r ) d * ) 2 ( 1 u ( r ) d * ) 2 δ | u ( r ) | 2 1 + δ | u ( r ) | 2 + ( 1 v ( r ) d * ) 2 δ | v ( r ) | 2 1 + δ | v ( r ) | 2 ) .
Observing that ( 1 u ( r ) d * ) 2 δ | u ( r ) | 2 1 + δ | u ( r ) | 2 > 0 , Equation (51) becomes:
F ( r , u ( r ) , u ( r ) ) F ( r , v ( r ) , v ( r ) ) < 1 r ( u ( r ) v ( r ) ) + 4 θ λ 2 ( ( 1 u ( r ) d * ) 2 ( 1 v ( r ) d * ) 2 + ( 1 v ( r ) d * ) 2 δ | v ( r ) | 2 1 + δ | v ( r ) | 2 ) .
Moreover, since δ | v ( r ) | 2 1 + δ | v ( r ) | 2 1 , Equation (52) can be written as
F ( r , u ( r ) , u ( r ) ) F ( r , v ( r ) , v ( r ) ) < 1 r ( u ( r ) v ( r ) ) + 4 θ λ 2 ( ( 1 u ( r ) d * ) 2 ( 1 v ( r ) d * ) 2 + ( 1 v ( r ) d * ) 2 ) .
Furthermore, with ( 1 v ( r ) d * ) 2 > 0 from (53)
F ( r , u ( r ) , u ( r ) ) F ( r , v ( r ) , v ( r ) ) < 1 r ( u ( r ) v ( r ) ) + 4 θ λ 2 ( ( 1 u ( r ) d * ) 2 + ( 1 v ( r ) d * ) 2 ) < 1 r ( u ( r ) v ( r ) ) + 4 θ λ 2 ( ( 1 u ( r ) ) 2 + ( 1 v ( r ) ) 2 ) < 1 r ( u ( r ) v ( r ) ) + 4 θ λ 2 ( 2 + u 2 ( r ) + v 2 ( r ) 2 u ( r ) 2 v ( r ) ) < 1 r ( u ( r ) v ( r ) ) + 4 θ λ 2 ( 2 + u ( r ) + v ( r ) 2 u ( r ) 2 v ( r ) ) = 1 r ( u ( r ) v ( r ) ) + 4 θ λ 2 ( 2 u ( r ) v ( r ) ) < 1 r ( u ( r ) v ( r ) ) + 4 θ λ 2 ( 2 + u ( r ) v ( r ) ) .
As it is safe to assume that u ( r ) > v ( r ) , and with u ( r ) < 1 and v ( r ) < 1 ,
2 < K ¯ ( u ( r ) v ( r ) )
with K ¯ = K ¯ ( r ) = sup { 2 / ( u ( r ) v ( r ) ) } | ( r ( r ) , v ( r ) ) U , from which K ¯ 2 u ( r ) v ( r ) so that (54) becomes
F ( r , u ( r ) , u ( r ) ) F ( r , v ( r ) , v ( r ) ) < 1 r L 1 ( r ) ( u ( r ) v ( r ) ) + 4 ( K ¯ ( r ) + 1 ) θ λ 2 K 2 ( r ) ( u ( r ) v ( r ) ) .
Moreover,
F ( r , u ( r ) , u ( r ) ) F ( r , v ( r ) , v ( r ) ) = 1 r ( u ( r ) v ( r ) ) + 4 ( 1 u ( r ) d * ) 2 θ λ 2 ( 1 + δ | u ( r ) | 2 ) 4 ( 1 v ( r ) d * ) 2 θ λ 2 ( 1 + δ | v ( r ) | 2 ) = 1 r ( u ( r ) v ( r ) ) + 4 θ λ 2 ( 1 u ( r ) d * ) 2 1 + δ | u ( r ) | 2 ( 1 v ( r ) d * ) 2 1 + δ | v ( r ) | 2 1 r ( u ( r ) v ( r ) ) + 4 θ λ 2 ( 1 u ( r ) d * ) 2 1 + δ | u ( r ) | 2 ( 1 v ( r ) d * ) 2 1 r ( u ( r ) v ( r ) ) + 4 θ λ 2 ( 1 u ( r ) d * ) 2 1 + δ | u ( r ) | 2 ( 1 v ( r ) d * ) = 1 r ( u ( r ) v ( r ) ) + 4 θ λ 2 ( ( 1 u ( r ) d * ) 2 1 + δ | u ( r ) | 2 ( 1 + δ | u ( r ) | 2 δ | u ( r ) | 2 ) ( 1 v ( r ) d * ) ) = 1 r ( u ( r ) v ( r ) ) + 4 θ λ 2 ( ( 1 u ( r ) d * ) 2 ( 1 u ( r ) d * ) 2 δ | u ( r ) | 2 1 + δ | u ( r ) | 2 ( 1 v ( r ) d * ) ) > 1 r ( u ( r ) v ( r ) ) + 4 θ λ 2 ( ( 1 u ( r ) d * ) 2 ( 1 u ( r ) d * ) 2 ( 1 v ( r ) d * ) ) = 1 r ( u ( r ) v ( r ) ) + 4 θ λ 2 ( ( 1 v ( r ) d * ) ) > 1 r ( u ( r ) v ( r ) ) + 4 θ λ 2 ( ( 1 v ( r ) d * ) ( 1 u ( r ) d * ) ) = 1 r ( u ( r ) v ( r ) ) + 4 θ λ 2 ( 1 + v ( r ) + d * 1 + u ( r ) + d * ) = 1 r ( u ( r ) v ( r ) ) + 4 θ λ 2 ( 2 + 2 d * + u ( r ) + v ( r ) ) 1 r ( u ( r ) v ( r ) ) + 4 θ λ 2 ( 2 + 2 d * + u ( r ) v ( r ) ) > 1 r ( u ( r ) v ( r ) ) + 4 θ λ 2 ( 2 + u ( r ) v ( r ) ) = 1 r ( u ( r ) v ( r ) ) 8 θ λ 2 + 4 ( u ( r ) v ( r ) ) θ λ 2 .
Supposing again that u ( r ) > v ( r ) , Equation (55) holds, from which
8 > 4 K ¯ ( r ) ( u ( r ) v ( r ) ) ;
therefore, Equation (57) becomes
F ( r , u ( r ) v ( r ) , u ( r ) v ( r ) ) > 1 r ( u ( r ) v ( r ) ) 8 θ λ 2 + 4 ( u ( r ) v ( r ) ) θ λ 2 = 1 r ( u ( r ) v ( r ) ) + 1 θ λ 2 ( 8 + 4 ( u ( r ) v ( r ) ) > 1 r ( u ( r ) v ( r ) ) + 1 θ λ 2 ( 4 K ¯ ( r ) ( u ( r ) v ( r ) ) + 4 ( u ( r ) v ( r ) ) = 1 r L 2 ( r ) ( u ( r ) v ( r ) ) + 4 4 K ¯ ( r ) θ λ 2 K 1 ( r ) ( u ( r ) v ( r ) ) .
Lastly, it is required that u 1 ( a ) u 2 ( a ) . For this goal, since a = r = 0 , we achieve u 1 ( a ) = u 1 ( 0 ) . Moreover, u 2 ( a ) = u 2 ( 0 ) = 0 . Furthermore, u 1 ( R ) = u 2 ( R ) = 0 . Thus, the proof of the theorem is complete. □
Remark 6.
Physically, as T increases, and V is fixed (if the intended use of the device was chosen), the membrane lifts from its rest position to a lesser extent than that when T is reduced. On the other hand, by increasing p e l , the deformation of the membrane is more accentuated. This is confirmed in (42). In fact, from (16),
R 2 = 4 k 1 T
and from (17), we obtain
k 1 k = 1 k 2 .
Thus, Equation (42) is writable as
θ λ 2 > 8 T d * 2 k 2 ϵ 0 V 2 1 + δ k 2 ϵ 0 V 2 r 4 T d * 2 2 .
Analyzing (62), increasing T also increases θ λ 2 , so that, observing (2), | u ( r ) | significantly decreases with the consequent reduction in the concavity of the membrane. Furthermore, increasing p e l and considering (17) increases k 2 , so that θ λ 2 decreases. Then, | u ( r ) | in (2) increases. We lastly observe that (42), when δ = 0 (i.e., without fringing field), we obtain the well-known algebraic condition studied in [12].

7. Solution Uniqueness

As verified above, Equation (2) admits at least one solution, u ( r ) , such that u 1 ( r ) < u ( r ) < u 2 ( r ) , where both u 1 ( r ) and u 2 ( r ) satisfy the hypotheses of Theorem 1. However, the uniqueness of the solution is not ensured, as proved in the following theorem.
Theorem 2.
If the hypotheses of Theorem 1 regarding (2) are satisfied, and u 1 ( r ) and u 2 ( r ) together satisfy the assigned boundary conditions, then the uniqueness of solution u ( r ) , such that u 1 ( r ) u ( r ) u 2 ( r ) , is not guaranteed.
Proof of Theorem 2.
As specified in (54),
L 1 ( r ) u ( r ) + K 2 ( r ) u ( r ) = 1 r u ( r ) + 4 ( K ¯ ( r ) + 1 ) θ λ 2 u ( r ) .
Therefore, considering Lemma 2, we can take into account the following ODE:
u ( r ) + L 1 ( r ) u ( r ) + K 2 ( r ) u ( r ) = 0 ,
which, in our case, since L 1 ( r ) = 1 r and K 2 ( r ) = 4 ( K ¯ ( r ) + 1 ) θ λ 2 , becomes
u ( r ) + 1 r u ( r ) + 4 ( K ¯ ( r ) + 1 ) θ λ 2 u ( r ) = 0 ,
which represents a special case of the following Bessel equation:
u ( r ) + 1 r u ( r ) + α 2 u ( r ) = 0 ,
with α = 4 ( K ¯ ( r ) + 1 ) θ λ 2 R + . From the Bessel theory of ordinary differential equations, the general solution for (65) can be written in terms of the linear combination of two linearly independent Bessel functions of the first and second kind of the zero-th order, J 0 4 ( K ¯ ( r ) + 1 ) θ λ 2 r and Y 0 4 ( K ¯ ( r ) + 1 ) θ λ 2 r , respectively. In other words:
u ( r ) = c 1 J 0 4 ( K ¯ ( r ) + 1 ) θ λ 2 r + c 2 Y 0 4 ( K ¯ ( r ) + 1 ) θ λ 2 r
where c 1 and c 2 are arbitrary constants, and
J 0 4 ( K ¯ ( r ) + 1 ) θ λ 2 r = 1 + m = 1 + ( 1 ) m 4 ( K ¯ ( r ) + 1 ) θ λ 2 2 m 2 2 m ( m ! ) 2
Y 0 4 ( K ¯ ( r ) + 1 ) θ λ 2 r = 2 π [ γ + ln 0.5 4 ( K ¯ ( r ) + 1 ) θ λ 2 r J 0 4 ( K ¯ ( r ) + 1 ) θ λ 2 r + m = 1 + ( 1 ) m + 1 H m 4 ( K ¯ ( r ) + 1 ) θ λ 2 2 m 2 2 m ( m ! ) 2 ]
in which γ = 0.5772 is the Euler–Mascheroni constant, and H m = 1 + 1 2 + 1 3 + + 1 m . We observe that, for r 0 , we obtain that J 0 1 . Y 0 also has logarithmic singularity when r = 0 . Therefore, considering a linear combination with c 2 = 0 and c 1 0 , we obtain the following general solution:
u ( r ) = c 1 1 + m = 1 + ( 1 ) m 4 ( K ¯ ( r ) + 1 ) θ λ 2 2 m 2 2 m ( m ! ) 2
Equation (70) represents a nontrivial solution for (65). Thus, by Lemma 2, it follows that the uniqueness of the solution to Problem (2) is not guaranteed because a subinterval of [ 0 , R ] exists (in our case, itself) on which (64) admits no trivial solution. □

8. On Research of Critical Points and Stability

8.1. A More Suitable Writing of the Differential Model

For searching any critical points, we need to rewrite (2) as a system of two ordinary differential equations of the first order in normal form [37]. For this purpose, let us consider two functions, u 1 ( r ) and u 2 ( r ) , such that
u 1 ( r ) = u ( r ) , u 2 ( r ) = u ( r ) .
Then, Equation (71) allows for rewriting (2) as a system of differential equations of the first order in which the unknown functions are the profile of the membrane u ( r ) and its speed of variation u ( r ) . In fact, from (2) and considering (71), we can write
u 1 ( r ) = u 2 ( r ) ; u 2 ( r ) = 1 r u 2 ( r ) ( 1 u 1 ( r ) d * ) 2 θ λ 2 ( 1 + δ | u 2 ( r ) | 2 ) u 1 ( R ) = u 2 ( 0 ) = 0 .

8.2. Critical Points and Stability

System (72) is a special case of the following general formulation:
u 1 ( r ) = f ¯ ( u 1 ( r ) , u 2 ( r ) ) ; u 2 ( r ) = g ¯ ( u 1 ( r ) , u 2 ( r ) ) ,
where, in our case,
f ¯ ( u 1 ( r ) , u 2 ( r ) ) = u 2 ( r )
and
g ¯ ( u 1 ( r ) , u 2 ( r ) ) = 1 r u 2 ( r ) ( 1 u 1 ( r ) d * ) 2 θ λ 2 ( 1 + δ | u 2 ( r ) | 2 ) .
Therefore, for achieving the critical points, we set u 1 ( r ) = u 2 ( r ) = 0 , which in our case, since θ λ 2 0 , we are given the following unique critical point:
( u 1 0 , u 2 0 ) = ( 1 d * , 0 ) .
To evaluate the stability of (76), we exploit the first Lyapunov criterion [37] that is based on the linearization of System (72) in the neighborhood of the critical point. For this purpose, we consider the following change of variable:
u 1 ( r ) = u 1 0 + ϵ ξ ( r ) ; u 2 ( r ) = u 2 0 + ϵ η ( r )
with ϵ being a small-enough quantity. Therefore, considering (73)–(75) and (77),
u 1 ( r ) = ϵ ξ ( r ) = f ¯ ( u 1 ( r ) , u 2 ( r ) ) u 2 ( r ) = ϵ η ( r ) = g ¯ ( u 1 ( r ) , u 2 ( r ) )
from which, developing in Taylor series both f ¯ ( u 1 ( r ) , u 2 ( r ) ) and g ¯ ( u 1 ( r ) , u 2 ( r ) ) , and neglecting the terms of an order higher than the linear one and setting τ = ξ 2 + η 2 , it follows that
ϵ ξ ( r ) = f ¯ ( u 1 0 + ϵ ( r ) , u 2 0 + ϵ η ( r ) ) f ¯ ( u 1 0 , u 2 0 ) + ϵ f ¯ u 1 ( u 1 0 , u 2 0 ) ξ ( r ) + ϵ f ¯ u 2 ( u 1 0 , u 2 0 ) η ( r ) + o ( τ ) ϵ η ( r ) = g ¯ ( u 1 0 + ϵ ( r ) , u 2 0 + ϵ η ( r ) ) g ¯ ( u 1 0 , u 2 0 ) + ϵ g ¯ u 1 ( u 1 0 , u 2 0 ) ξ ( r ) + ϵ g ¯ u 2 ( u 1 0 , u 2 0 ) η ( r ) + o ( τ ) .
Remark 7.
Equation (79) makes sense because both u 1 ( r ) and u 2 ( r ) are analytical functions allowing for the linearization procedure by means of computing f ¯ u 1 , f ¯ u 2 , g ¯ u 1 and g ¯ u 2 .
In (79), f ¯ ( u 1 0 , u 2 0 ) = g ¯ ( u 1 0 , u 2 0 ) = 0 , so that it becomes
ξ ( r ) = f ¯ u 1 ( u 1 0 , u 2 0 ) ξ ( r ) + f ¯ u 2 ( u 1 0 , u 2 0 ) η ( r ) η ( r ) = g ¯ u 1 ( u 1 0 , u 2 0 ) ξ ( r ) + g ¯ u 2 ( u 1 0 , u 2 0 ) η ( r ) .
Moreover,
f ¯ u 1 ( u 1 0 , u 2 0 ) = 0 , f ¯ u 2 ( u 1 0 , u 2 0 ) = 1 , g ¯ u 2 ( u 1 0 , u 2 0 ) = 1 r
and
g ¯ u 1 ( u 1 0 , u 2 0 ) = 2 ( 1 ( 1 d * ) d * ) = 0 ,
so that System (79) can be written as
ξ ( r ) = η ( r ) η ( r ) = η ( r ) r
that, when solved, gives us
ξ + ln η C 1 = C 2
with both C 1 and C 2 being constant. Equation (82), indicated by
z = ξ ( r ) η ( r ) , z ˙ = ξ ( r ) η ( r ) , A = 0 1 0 1 r ,
can easily be written as
z ˙ = A z .
Let us consider the following definitions [38].
Definition 1.
If A is a square matrix with order n, with r n distinct eigenvalues λ 1 , λ 2 , … λ r , i j , the characteristic polynomial of A is
P ( s ) = ( s λ 1 ) ν 1 ( s λ 2 ) ν 2 ( s λ r ) ν r
such that i = 1 r = n . ν i N + defines the algebraic multiplicity of generic eigenvalue λ i . Moreover, the geometric multiplicity of λ i is defined by number μ i of the linearly independent eigenvectors corresponding to it.
Definition 2.
If λ is an eigenvalue of A with algebraic multiplicity ν linked to Jordan form J , the index of λ, π, is the order of the largest Jordan block associated with λ in J .
Now, we present the following result evaluating the stability of (85) exploiting a criterion on the basis of A eigenvalues [38].
Theorem 3.
Dynamic systems in form (85) admit at least one stable equilibrium position if, and only if, its matrix, A, does not have eigenvalues with a positive real part and if any eigenvalues with a real part zero have a unit index. Furthermore, if z 0 = [ z 0 , 1 , z 0 , 2 ] T , then
z ˙ ( r ) = e A r z ( 0 ) = e A r z 0 .
Thus, the following result holds.
Proposition 1.
Equation (82) admits a stable equilibrium position.
Proof. 
From A, two eigenvalues, λ 1 and λ 2 , are computed:
λ 1 = 0 ; λ 2 = 1 r ,
so that, by Theorem 3, it follows that System (82) is stable, and (76) is a stable equilibrium position. □
Remark 8.
The number of the eigenvalues of A, counted with their algebraic multiplicity, is equal to the order of A. Moreover, the geometric multiplicity of each eigenvalue is equal to the algebraic multiplicity. Then, it follows that A is diagonalizable [38].
Therefore, since A is diagonalizable (see Remark 8), e A r can be written as [38]:
e A r = k = 1 n t k · s k T e λ k r = t 1 · s 1 + t 2 · s 2 e 1
in which t k and s k are, respectively, the left and right eigenvectors corresponding to λ k . Thus,
t 1 = [ 1 0 ] T , t 2 = [ 1 r 1 ] T , s 1 = [ 1 r ] , s 2 = [ 0 1 ] ,
so that (89) becomes
e A r = 1 r r + 1 e 0 1 e r .
Moreover, Equation (91) for r 0 is limited in norm, so that (87) becomes
ξ ( r ) = ξ 0 r + ( r + e 1 ) η 0 η ( r ) = η 0 r e
from which, eliminating r, we achieve
ξ ( r ) η ( r ) = ξ 0 e η 2 ( r ) η 0 η 0 2 e + η 0 e η ( r ) .
Indicated by
H = ξ 0 e η 0 1 2 1 1 2 0 0 1 0 η 0 2 e
and
H 1 = ξ 0 e η 0 1 2 1 2 0
one easily achieves that | H | 0 and | H 1 | < 0 , so that (93) on ξ η plane represents a hyperbole. Furthermore, from (92),
η ( r ) = ξ ( r ) r ( e r + 1 ) + ξ 0 r 2 ( e r + 1 ) ,
in which r ( 0 , R ] represents a straight line of which the slope, as r varies on ( 0 , R ] , changes from 1 r ( e r + 1 ) to 1 R ( e R + 1 ) , intercepting point ( x i 0 r , 0 ) on axis η ( r ) = 0 and point ( 0 , x i 0 r 2 ( 1 + e r ) ) on axis ξ ( r ) = 0 . Therefore, the straight line, when its slope changes, intercepts an arc of hyperbole representing points ( ξ ( r ) , η ( r ) ) . This arc of hyperbole represents the place of points ξ ( r ) η ( r ) , which admit stability for (85).
Concerning the stability of System (72), let us introduce the following lemma [37].
Lemma 3.
If linearized System (85) is stable, the critical point of nonlinear System (72) is stable itself.
Therefore, we can introduce the following important result.
Proposition 2.
Critical point (76) is an equilibrium position characterized by stability for System (72).
Proof of Proposition 2.
It is an immediate consequence of Lemma 3. □
Remark 9.
Point ( 1 d * , 0 ) identifies a profile of the membrane when u 0 is very close to the upper disk. This poses a risk, as the membrane could touch the upper disk. Electrostatically, considering that
1 ( 1 u ( r ) ) 2 1 d * ,
p e l becomes
p e l = 1 2 ϵ 0 V 2 ( 1 u ( r ) ) 2 ϵ 0 V 2 2 d * .
Thus, once V is fixed, p e l does not fluctuate, so that any variations of p e l are not appreciable. Then, from (17), not even any fluctuations of p are appreciable, so that if the membrane reaches the unstable equilibrium position, the risk of touching the upper disk is minimal.
Remark 10.
The only equilibrium position that we obtained was fixed by quantity 1 d * . It was a constant amount once safety distance d * was set. The greater d * is, the lower the deformation of the membrane under stable conditions. Obviously, a lower bound is obtained if d * = 1 . This lower-bound condition, however, requires that the membrane does not deform despite the application of an electrical voltage V, even of strong amplitude. We, therefore, deduce that the condition of a stable lower bound, even if it has full meaning from a mathematical point of view, from a physical point of view, it has no relevance because it does not admit membrane deformations even if electrostatic stress is relevant.

9. On Admissible Values of V

9.1. Minimal Value of V to Win the Mechanical Inertia of the Membrane

Proposition 3.
Let us consider Model (2) and its condition of existence (42). Thus, the minimal value of V to win membrane mechanical inertia ( V m i n ) i n e r t i a satisfies the following inequality:
( V m i n ) i n e r t i a = T d * 6 d 3 ϵ 0 2 + 4 δ k 2 p e l 2 4 .
Proof of Proposition 3.
From (42),
λ 2 > 2 R 2 d * 2 θ V 2 ϵ 0 k 1 + δ k ϵ 0 V 2 r d * 2 R 2 2 2 R 2 d * 2 θ V 2 ϵ 0 k 1 + δ k ϵ 0 V 2 d * 2 2 = 2 R 2 d * 6 θ V 2 ϵ 0 k ( d * 4 + δ ( k ϵ 0 V 2 ) 2 ) 2 R 2 d * 6 θ V 2 ϵ 0 k ( d * 4 + δ ( k V 2 ) 2 )
from which, considering that [1]
λ 2 = 2 ϵ 0 V 2 R 2 d 3 T ,
we achieve
2 ϵ 0 V 2 R 2 d 3 T 2 R 2 d * 6 θ V 2 ϵ 0 k ( d * 4 + δ ( k V 2 ) 2 ) ,
and again
V 4 ( d * 4 + δ ( k V 2 ) 2 ) T d * 6 d 3 θ k ϵ 0 2 .
However, from (18),
k V 2 = 2 k p e l ( d u ( r ) ) 2 ϵ 0
from which
1 + δ ( k V 2 ) 2 = 1 + 4 δ k 2 p e l 2 ( d u ( r ) ) 4 ϵ 0 2 ,
so that (103) becomes
V 4 ( 1 + δ ( k V 2 ) 2 ) = V 4 1 + 4 δ k 2 p e l 2 ( d u ( r ) ) 4 ϵ 0 2 T d * 6 d 3 θ k ϵ 0 2
from which, with θ k < 1 and d u ( r ) < 1 , we can write
V 4 T d * 6 d 3 θ k ( ϵ 0 2 + 4 δ k 2 p e l 2 ( d u ( r ) ) 4 ) T d * 2 d 3 ϵ 0 2 + 4 δ k 2 p e l 2
and
V T d * 6 d 3 ϵ 0 2 + 4 δ k 2 p e l 2 4
obtaining the (99). □
Remark 11.
Equation (99) makes sense because T appears in the numerator of its right side, i.e., the mechanical tension of the membrane at rest. Then, the greater T is, the greater V must be to overcome the mechanical inertia of the membrane.

9.2. Maximal Value of V so the Membrane Does Not Touch the Upper Disk

Remark 12.
Equation (21) was obtained by exploiting the theory of elasticity of circular membranes [1] when external V is applied. Since the considered device is circular, Equation (21) retains its validity and the geometry has axial symmetry (with axis r = 0 ). This symmetry was maintained even in the presence of a fringing field because this phenomenon (more evident at the edges of the device) was also symmetrical with respect to the same vertical axis. Therefore, the presence of the fringing field would seem not to invalidate the validity of (21). However, in (21), there is no trace of C e l which represents the most influenced electrostatic parameter by the fringing-field effect, so the presence of terms due in the fringing field is not explicitly evident. On the other hand, studying the existence of the solution for (2), made (21) provide an algebraic condition in which the fringing field (presence of δ) was evident. It follows that, even in the presence of a fringing field, Equation (21) is still valid. Therefore, the achieved results in [12] and concerning the maximal value of V are still valid.
In particular, the following results hold [12].
Proposition 4.
Indicating by ( V m a x ) p e r m i s s i b l e the maximum value of V in order that the membrane does not touch the upper desk, for Model (2), the following inequality holds (for details, see [12]):
( V m a x ) p e r m i s s i b l e < 2 d * ( 1 d * ) k ϵ 0
Thus, considering both (108) and (109), the range of the admissible values for V is
T d * 6 d 3 ϵ 0 2 + 4 δ k 2 p e l 2 4 V < 2 d * ( 1 d * ) k ϵ 0 .
Remark 13.
Then, Equation (110) makes sense because
T d * 6 d 3 ϵ 0 2 + 4 δ k 2 p e l 2 4 < 2 d * ( 1 d * ) k ϵ 0 .
In fact, if absurdly
T d * 6 d 3 ϵ 0 2 + 4 δ k 2 p e l 2 4 > 2 d * ( 1 d * ) k ϵ 0 ,
one would obtain
δ < ϵ 0 2 T d 3 k 2 d * 4 4 ( 1 d * ) 2 < 1 1 < 0 1 4 k 2 p e l 2
so that δ < 0 , which represents an impossible condition because δ R + . Therefore, Equation (111) makes sense, so that (110) is true.

10. Interesting Optimal Control Conditions

If the membrane is at rest, its distance from the upper disk is d, and the C e l of the device along any plane whose support is the straight line r = 0 , in presence of fringing field, is [30]
( C e l ) c u r v e = 2 ϵ 0 R d 1 + d 2 π R ln 2 π R d ,
where d 1 represents the distance so that the total C e l , ( C e l ) t o t a l , becomes:
( C e l ) t o t a l = 0 π Z ( ϕ ) ( C e l ) c u r v e d ϕ = ( C e l ) c u r v e 0 π Z ( ϕ ) d ϕ .
in which Z ( ϕ ) is a bounded and a continuous electrostatic function depending on angular coordinate ϕ [39]. Thus, Z ( ϕ ) being both continuous and bounded,
0 π Z ( ϕ ) d ϕ = D < + ,
so that (115) becomes
( C e l ) t o t a l = ( C e l ) c u r v e D = 2 ϵ 0 R D d 1 + d 2 π R ln 2 π R d .
Thus, the potential energy (indicated by W i n i t i a l ) when the device is at rest is
W i n i t i a l = 1 2 ( C e l ) t o t a l V 2 = ϵ 0 R D V 2 d 1 + d 2 π R ln 2 π R d .
If the membrane deforms, and in the presence of a fringing field, C e l d e f , along with any curve obtained by intersecting the profile of the deformed membrane with any plane of which the support is line r = 0 , would become [30]
( C e l d e f ) c u r v e = 2 ϵ 0 R R + R d r d u ( r ) 1 + 1 2 π R R + R d r d u ( r ) ln 2 π R R + R d r d u ( r ) .
However,
1 d u ( r ) 1 d *
from which
R + R d r d u ( r ) R + R d r d * = 2 R d * .
Thus,
( C e l d e f ) c u r v e 4 ϵ 0 R 2 d * 1 + 1 2 π R R + R d r d u ( r ) ln 4 π R 2 d * .
Moreover,
d u ( r ) d ,
from which
1 d u ( r ) 1 d
and
R + R d r d u ( r ) R + R d r d
Thus, Equation (122), considering (125), becomes
( C e l d e f ) c u r v e 4 ϵ 0 R 2 d * 1 + d 4 π R 2 ln 4 π R 2 d * .
Therefore, as (117),
( C e l d e f ) t o t a l 4 ϵ 0 R 2 D d * 1 + d 4 π R 2 ln 4 π R 2 d * .
From (127), considering (109), it follows that
W f i n a l = 1 2 ( C e l d e f ) t o t a l V 2 2 ϵ 0 R 2 D d * 2 d * ( 1 d * ) k ϵ 0 1 + d 4 π R 2 ln 4 π R 2 d * .
from which
( W f i n a l ) max = 2 ϵ 0 R 2 D d * 2 d * ( 1 d * ) k ϵ 0 1 + d 4 π R 2 ln 4 π R 2 d * .
Thus, considering both (118) and (129),
Δ W = ( W f i n a l ) max W i n i t i a l 2 ϵ 0 R 2 D d * 2 d * ( 1 d * ) k ϵ 0 1 + d 4 π R 2 ln 4 π R 2 d * ϵ 0 R D V 2 d 1 + d 2 π R ln 2 π R d
On the other hand, from (110),
V 2 > T d * 6 d 3 ϵ 0 2 + 4 δ k 2 p e l 2 ;
thus,
Δ W > ϵ 0 R D T d * 6 d 3 ϵ 0 2 + 4 δ k 2 p e l 2 { 2 R d * 1 + d 4 π R 2 ln 4 π R 2 d * 1 d 1 + d 2 π R ln 2 π R d } .

10.1. V Maximizing Δ W

From (130), it is clear that we obtain an upper limitation for the value of V, which maximizes Δ W is ( V m a x ) p e r m i s s i b l e .

10.2. An Interesting Limitation for Δ W Starting from | E |

From the equations of Models (22) and (2):
λ 2 ( 1 u ( r ) ) 2 ( 1 + δ | u ( r ) | 2 ) = u ( r ) r + 1 θ λ 2 ( 1 u ( r ) d * ) 2 ( 1 + δ | u ( r ) | 2 ) .
Cconsidering that λ 2 ( 1 u ( r ) ) 2 = θ | E | 2 , Equation (133) becomes:
| E | 2 ( 1 + δ | u ( r ) | 2 ) = u ( r ) r θ + 1 θ 2 λ 2 ( 1 u ( r ) d * ) 2 ( 1 + δ | u ( r ) | 2 ) ,
so that
W f i n a l = 1 2 ϵ 0 | E | 2 = ϵ 0 2 u ( r ) r θ ( 1 + δ | u ( r ) | 2 ) + 1 θ 2 λ 2 ( 1 u ( r ) d * ) 2 ( 1 + δ | u ( r ) | 2 ) 2 .
Therefore, considering both (118) and (135),
Δ W = W f i n a l W i n i t i a l = ϵ 0 2 u ( r ) r θ ( 1 + δ | u ( r ) | 2 ) + 1 θ 2 λ 2 ( 1 u ( r ) d * ) 2 ( 1 + δ | u ( r ) | 2 ) 2 ϵ 0 R D V 2 d 1 + d 2 π R ln 2 π R d .
Moreover, with 1 u ( r ) d * < 1 d * and u ( r ) < H R + (see, [11,12]), Equation (136) becomes:
Δ W = W f i n a l W i n i t i a l < ϵ 0 2 H r θ ( 1 + δ | u ( r ) | 2 ) + 1 θ 2 λ 2 ( 1 d * ) 2 ( 1 + δ | u ( r ) | 2 ) 2 ϵ 0 R D V 2 d 1 + d 2 π R ln 2 π R d .
Lastly, with
1 1 + δ | u ( r ) | 2 1
Equation (137) becomes
Δ W < ϵ 0 2 H r θ + ( 1 d * ) 2 θ 2 λ 2 ϵ 0 R D V 2 d 1 + d 2 π R ln 2 π R d .
In (139), the right side is non-negative. In fact, if it were negative, we would easily achieve
H < r θ 2 R D V 2 d 1 + d 2 π R ln 2 π R d ( 1 d * ) 2 θ 2 λ 2
from which, substituting the usual values for each parameter, we would find that H was increased by a very small non-negative quantity. This means that the slope of the membrane at the edges is very small, as if the deformation of the membrane was extremely small even for high values of V. It is evident that this condition is physically impossible. Hence, it follows that the right side of (139) is always a non-negative quantity.

11. Conclusions

In this study, a new nonlinear second-order differential 2D model for a membrane MEMS device in which fringing-field effects occur was presented and discussed. Once the possibility of formulating | E | in terms of average membrane curvature was justified, the mathematical model was detailed from the point of view of the actuator, highlighting the link between p and p e l , highlighting the actuator–transducer duality. Once two important lemmas were presented, the existence of at least one solution for the proposed model was ensured if an algebraic condition depending on fringing-field parameter δ was verified. However, even if the uniqueness of the solution was not assured, the only permitted equilibrium position of the model is stable. Furthermore, this equilibrium position is associated with the profile of the membrane that was closest to the upper disk; therefore, even if the membrane was very close to the upper disk, the stability of this equilibrium configuration reduced the risk of touching the upper disk. The range of possible values for V was obtained by taking into account the need, on the one hand, to overcome the mechanical inertia of the membrane and, on the other, to prevent the membrane from touching the upper disk, highlighting that the fringing-field effect appeared explicitly only in the extreme of V, delegated to overcome the mechanical inertia of the membrane. Lastly, some conditions for optimal membrane control were obtained. In particular, an increase was obtained for Δ W when the membrane passed from the rest position to a generic deformed configuration. This is very interesting because it contains all the parameters that came into play when the fringing-field phenomenon occurs. Sophisticated mathematical models describing the behavior of the membrane under the effect of V and in the presence of fringing fields, while offering excellent food for thought, are poorly suited to real-time industrial applications. Then, some simplifications of the model and, in our case, of the formulation of the mean curvature, appear necessary. Locally, this is proportional to | E | , but determines the effective deformation of the membrane inside the device. However, the study of the model, although it does not provide results that can be superimposed with experimental data, certainly provides qualitative indications on the electromechanical behavior of a membrane MEMS device. Lastly, the mathematical models adhering to the physical reality of MEMS modeling are extremely complex and obviously do not allow for in-depth analytical studies. Then, some simplifications in the geometry of the devices are necessary to obtain simplified analytical models so that they can be easily studied. Obviously, the results obtained by studying the model proposed in this paper poorly agree with any experiment results, but provide interesting qualitative indications of the behavior of MEMS device membranes characterized by simplified geometries.

Author Contributions

Conceptualization, M.V. and L.F.; methodology, M.V. and L.F.; validation, L.F. and P.D.B.; formal analysis, M.V. and L.F.; investigation, M.V.; resources, L.F.; writing—original draft preparation, M.V.; writing—review and editing, L.F. and P.D.B.; supervision, L.F. and P.D.B. Luisa Fattorusso is Member of the National Group for Mathematical Analysis, Probability, and their Applications (GNAMPA-INdAM). All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
rradial coordinate
Rradius of the membrane
u ( r ) profile of the membrane
Vexternal voltage
λ 2 parameter depending on V
ddistance between the parallel disks
E electrostatic field
θ coefficient of proportionality between u ( r ) and | E | 2
δ parameter concerning the fringing field effect
K ( r , u ( r ) ) mean curvature of the membrane
d * critical security distance
( V m i n ) i n e r t i a V to win the mechanical inertia of the membrane
( V m a x ) p e r m i s s i b l e V to avoid that the membrane touches the upper disk
p e l electrostatic pressure
ϵ 0 permittivity of the free space
C e l electrostatic capacitance
pmechanical pressure
ρ density
hthickness of the plate
Tmechanical tension of the membrane at rest
YYoung modulus
ν Poisson ratio
u 0 displacement at the center of the membrane
f e l electrostatic force
k 1 coefficient of proportionality between u 0 and p
k 2 coefficient of proportionality between p e l and p
kcoefficient of proportionality between u 0 and p e l
u 1 ( r ) , u 2 ( r ) twice continuously differentiable functions (upper and lower solutions)
K 1 ( r ) , K 2 ( r ) , L 1 ( r ) , L 2 ( r ) continuous functions
Z ( ϕ ) bounded and continuous electrostatic function
( C e l ) c u r v e total electrostatic capacitance when the membrane is at rest
Dconstant
( C e l d e f ) c u r v e electrostatic capacitance when the membrane is deformed
( C e l d e f ) t o t a l total electrostatic capacitance when the membrane is deformed

References

  1. Pelesko, J.A.; Bernstein, D.H. Modeling MEMS and NEMS; Chapman and Hall, CRC Press Company: Boca Raton, FL, USA; London, UK; New York, NY, USA; Washington, DC, USA, 2003. [Google Scholar]
  2. Gad-el-Hak, M. MEMS: Design and Fabrication; Chapman and Hall, CRC Taylor and Francis: Boca Raton, FL, USA, 2006. [Google Scholar]
  3. Cauchi, M.; Grech, I.; Mallia, B.; Mollicone, P.; Sammut, N. Analytical, Numerical and Experimental Study of a Horizontal Electrothermal MEMS Microgripper fo the Deformability Characterization of Human Red Blood Cells. Micromachines 2018, 9, 108. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Nathanson, H.; Newell, W.; Wickstrom, R.; Lewis, J. The Resonant Gate Transistor. IEEE Trans. Electron Devices 1964, 14, 117–133. [Google Scholar] [CrossRef]
  5. Cassani, D.; Tarsia, A. Periodic Solutions to Nonlocal MEMS Equations. Discret. Contin. Dyn. Syst. Ser. S 2016, 12, 188–199. [Google Scholar] [CrossRef]
  6. Cassani, D.; do Ó, J.M.; Ghoussoub, N. On a Fourth Order Elliptic Problem with a Singular Nonlinearity. Nonlinear Stud. 2009, 9, 189–209. [Google Scholar] [CrossRef]
  7. Zhang, Y.; Wang, T.; Luo, A.; Hu, Y.; Li, X.; Wang, F. Micro Electrostatic Enerhy Harvester with both Broad Bandwidth and High Normalized Power Density. Appl. Energy 2018, 212, 363–371. [Google Scholar] [CrossRef]
  8. Mohammadi, A.; Ali, N. Effect of High Electrostatic Actuation on Thermoelastic Dameing in Thin Rectangular Microplate Resonators. J. Theor. Appl. Mech. 2015, 53, 317–329. [Google Scholar] [CrossRef] [Green Version]
  9. Vinyas, M.; Kattimani, S. Investigation of the Effect of BaTiO3-CoFe24 Particle Arrangement on the Static Response of Magneto-Electro-Thermo-Elastic Plates. Compos. Struct. 2018, 185, 51–56. [Google Scholar] [CrossRef]
  10. Feng, J.; Liu, C.; Zhang, W.; Hao, S. Static and Dynamic Mechanical Behaviors of Electrostatic MEMS Resonator witgh Surface Processing Error. Micromachines 2018, 9, 34. [Google Scholar] [CrossRef] [Green Version]
  11. Di Barba, P.; Fattorusso, L.; Versaci, M. Electrostatic Feld in Terms of Geometric Curvature in Membrane MEMS Devices. Commun. Appl. Ind. Math. 2017, 8, 165–184. [Google Scholar]
  12. Di Barba, P.; Fattorusso, L.; Versaci, M. A 2D Non-Linear Second-Order Differential Model for Electrostatic Circular Membrane MEMS Devices: A Result of Existence and Uniqueness. Mathematics 2019, 7, 1193. [Google Scholar] [CrossRef] [Green Version]
  13. Pelesko, J.A.; Driscoll, T.A. The Effect of the Small-Aspect-Ratio Approximation on Canonical Electrostatic MEMS Models. J. Eng. Math. 2005, 53, 129–152. [Google Scholar] [CrossRef] [Green Version]
  14. Angiulli, G.; Jannelli, A.; Morabito, F.C.; Versaci, M. Reconstructing the Membrane Detection of a 1D Electrostatic-Driven MEMS Device by the Shooting Method: Convergence Analysis and Ghost Solutions Identification. Comp. Appl. Math. 2018, 37, 4484–4498. [Google Scholar] [CrossRef]
  15. Versaci, M.; Angiulli, G.; Fattorusso, L.; Jannelli, A. On the Uniqueness of the Solution for a Semi-Linear Elliptic Boundary Value Problem of the Membrane MEMS Device for Reconstructing the Membrane Profile in Absence of Ghost Solutions. Int. J. Non-Linear Mech. 2019, 109, 24–31. [Google Scholar] [CrossRef]
  16. Versaci, M.; Jannelli, A.; Angiulli, G. Electrostatic Micro-Electro-Mechanical-Systems (MEMS) Devices: A Comparison Among Numerical Techniques for Recovering the Membrane Profile. IEEE Access 2020, 8, 125874–125886. [Google Scholar] [CrossRef]
  17. Versaci, M.; Di Barba, P.; Morabito, F.C. Curvature-Dependent Electrostatic Field as a Principle for Modelling Membrane-Based MEMS Devices. A Review. Membranes 2020, 10, 361. [Google Scholar] [CrossRef] [PubMed]
  18. Versaci, M.; Morabito, F.C. Membrane Micro Electro-Mechanical Systems for Industrial Applications. In Handbook of Research on Advanced Mechatronic Systems and Intelligent Robotics; IGI Global: Hershey, PA, USA, 2003; pp. 139–175. [Google Scholar]
  19. Zozulya, V.V.; Saez, A. A High-Order Theory of a Thermoelastic Beams and Its Application to the MEMS/NEMS Analysis and Simulations. A Review. Arch. Appl. Mech. 2016, 86, 1255–1273. [Google Scholar] [CrossRef]
  20. Zega, V.; Frangi, A.; Guercilena, A.; Gattere, G. Analysis of Frequency Stability and Thermoelastic Effects for Slotted Tuning Fork MEMS Resonators. Sensors 2018, 8, 2157. [Google Scholar] [CrossRef] [Green Version]
  21. Sravani, K.G.; Narayana, T.L.; Guha, K.; Rao, K.S. Role of Dielectric Layer and Beam Membrane in Improving the Performance of Capacitive RF MEMS Switches for Ka-Band Applications. Math. Technol. 2018, 9, 145–156. [Google Scholar] [CrossRef]
  22. Neff, B.; Casset, F.; Millet, A.; Agache, V.; Verplanck, N.; Boizot, F.; Fanget, S. Development and Characterization of MEMS Membrane Based on Thin-Film PZT Actuators for Microfluidic Applications. In Proceedings of the 19th International Conference on Thermal, Mechanical and Multi-Physics Simulation and Experiments in Microelectronics and Microsystems (EuroSimE), Toulouse, France, 15–18 April 2018. [Google Scholar]
  23. Javaher, H.; Ghanati, P.P.; Azizi, S. A Case Study on the Numerical Solution and Reduced Order Model of MEMS. Sensors 2018, 19, 35–48. [Google Scholar] [CrossRef]
  24. de Oliveira Hansen, R.; Mátéfi-Tempfli, M.; Safonovs, R.; Adam, J.; Chemnitz, S.; Reimer, T.; Wagner, B.; Benecke, W.; Mátéfi-Tempfli, S. Magnetic Films for Electromagnetic Actuation in MEMS Switches. Microsyst. Technol. 2018, 24, 1987–1994. [Google Scholar] [CrossRef]
  25. Di Barba, P.; Wiak, S. MEMS: Field Models and Optimal Design; Springer International Publishing: Berlin/Heidelberg, Germany, 2020. [Google Scholar]
  26. Di Barba, P.; Gotszalk, T.; Majastrzyk, W.; Mognasti, M.; Orlowska, K. Optimal Design of Electromagnetically Actuated MEMS Cantilevers. Sensors 2018, 18, 2533. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Velosa-Moncada, L.A.; Aguilera-Cortés, L.A.; González-Palacios, M.A.; Raskin, J.P.; Herrera-May, A.L. Design of a Novel MEMS Microgripper with Rotatory Electrostatic Combdrive Actuators for Biomedical Applications. Sensors 2018, 18, 1664. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Nisar, A.; Afzulpurkar, N.; Mahaisavariya, B.; Tuantranont, A. MEMS-Based Micropumps in Drug Delivery and Biomedical Applications. Sens. Actuators B Chem. 2008, 2, 917–942. [Google Scholar] [CrossRef]
  29. Wei, J.; Ye, D. On MEMS Equation with Fringing Field. Proc. Am. Math. Soc. 2010, 138, 1693–1699. [Google Scholar] [CrossRef] [Green Version]
  30. Leus, V.; Elata, D. Fringing Field Effect in Electrostatic Actuator. In Technical Report ETR-2004-2; Israel Institute of Technology, Faculty of Mechanical Engineering: Haifa, Israel, 2004. [Google Scholar]
  31. Hosseini, M.; Zhu, G.; Peter, Y.A. A New Formulation of Fringing Capacitance and its Application to the Control of Parallel-Plate Electrostatic Micro Actuators. Analog. Integr. Circ. Sig. Process. 2007, 53, 119–128. [Google Scholar] [CrossRef]
  32. Li, Y.; Li, H.; Pickard, F.C., IV; Narayanan, B.; Sen, F.G.; Chan, M.K.; Sankaranarayanan, S.K.; Brooks, B.R.; Roux, B. Machine Learning Force Field Parameters from Ab Initio Data. J. Chem. Theory Comput. 2017, 13, 4492–4503. [Google Scholar] [CrossRef] [Green Version]
  33. Imani, M.; Ghoreishi, S.F. Scalable Inverse Reinforcement Learning Through Multi-Fidelity Bayesian Optimization. IEEE Trans. Meural Netw. Learn. Syst. 2021, 1–8. [Google Scholar] [CrossRef]
  34. Roudneshin, M.; Sayrafian-Pour, K.; Aghdam, A. A Machine Learning Approach to the Estimation of Near-Optimal Electrostatic Force in Micro Energy-Harvesters. In Proceedings of the 2019 IEEE International Conference on Wireless for Space and Extreme Environments, Ottawa, ON, Canada, 16–18 October 2019. [Google Scholar]
  35. Timoshenko, S.; Woinowsly-Krieger, S. Theory of Plates and Shells; McGraw Hill: New York, NY, USA, 1959. [Google Scholar]
  36. Bayley, P.B.; Shampine, L.F.; Waltman, P.E. Nonlinear Two Point Boundary Value Problems; Academic Press: Cambridge, MA, USA, 1968. [Google Scholar]
  37. Di Barba, P.; Fattorusso, L.; Versaci, M. Curvature Dependent Electrostatic Field in the Deformable MEMS Device: Stability and Optimal Control. Commun. Appl. Ind. Math. 2020, 11, 35–54. [Google Scholar]
  38. Barreira, L.; Valls, C. Dynamical Systems; Springer: London, UK, 2020. [Google Scholar]
  39. Fujimoto, M. Physics of Classical Electromagnetism; Springer: New York, NY, USA, 2007. [Google Scholar]
Figure 1. Work flowchart.
Figure 1. Work flowchart.
Mathematics 09 00465 g001
Figure 2. u 1 ( r ) and u 2 ( r ) for model (2).
Figure 2. u 1 ( r ) and u 2 ( r ) for model (2).
Mathematics 09 00465 g002
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Di Barba, P.; Fattorusso, L.; Versaci, M. A 2D Membrane MEMS Device Model with Fringing Field: Curvature-Dependent Electrostatic Field and Optimal Control. Mathematics 2021, 9, 465. https://doi.org/10.3390/math9050465

AMA Style

Di Barba P, Fattorusso L, Versaci M. A 2D Membrane MEMS Device Model with Fringing Field: Curvature-Dependent Electrostatic Field and Optimal Control. Mathematics. 2021; 9(5):465. https://doi.org/10.3390/math9050465

Chicago/Turabian Style

Di Barba, Paolo, Luisa Fattorusso, and Mario Versaci. 2021. "A 2D Membrane MEMS Device Model with Fringing Field: Curvature-Dependent Electrostatic Field and Optimal Control" Mathematics 9, no. 5: 465. https://doi.org/10.3390/math9050465

APA Style

Di Barba, P., Fattorusso, L., & Versaci, M. (2021). A 2D Membrane MEMS Device Model with Fringing Field: Curvature-Dependent Electrostatic Field and Optimal Control. Mathematics, 9(5), 465. https://doi.org/10.3390/math9050465

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop