The existence and multiplicity of solutions to a quasilinear, elliptic partial differential equation with singular non-linearity is analysed. The partial differential equation is a recently derived variant of a canonical model used in the modelling of micro-electromechanical systems. It is observed that the bifurcation curve of solutions terminates at single dead-end point, beyond which no classical solutions exist. A necessary condition for the existence of solutions is developed, revealing that this dead-end point corresponds to a blow-up in the solution's gradient at a point internal to the domain. By employing a novel asymptotic analysis in terms of two small parameters, an accurate characterization of this dead-end point is obtained. An arc length parameterization of the solution curve can be employed to continue solutions beyond the dead-end point; however, all extra solutions are found to be multi-valued. This analysis therefore suggests that the dead-end is a bifurcation point associated with the onset of multi-valued solutions for the system.

To study electrostatic actuation, researchers commonly use a setup proposed by G. I. Taylor in [Proc. R. Soc. Lond. Ser. A, 306 (1968), pp. 423–434]. It consists of soap film held at a distance h above a rigid plate so that when a voltage difference is applied between the two components, the top film deflects towards the bottom plate. The most striking feature of this system is when the voltage difference exceeds a critical value V*, the electrostatic forces dominate the surface forces and the soap film gets ‘pulled-into’ or collapses onto the bottom plate. This so-called ‘pull-in’ instability is a ubiquitous feature of electrostatic actuation and as a result, has been the subject of many studies. Recently, Siddique et al. [J. Electrostatics, 69 (2011), pp. 1–6] measured the value of V* as a function of the separation distance and found that the standard prediction breaks down as h increases. Here, we continue the work done in [N. D. Brubaker and J. A. Pelesko, European J. Appl. Math., 22 (2011), pp. 455–470] by investigating the cause of this discrepancy. Specifically, we model the effect of gravity on the generalized version of Taylor's model and study whether it provides the proper correction to the predicted value of V*. In doing so, we derive two nonlinear eigenvalue value problems and investigate their solutions sets.

In modelling electrostatically actuated micro- and nano-electromechanical systems, researchers have typically relied on a small-aspect ratio to form a leading-order theory. In doing so, small gradient terms are dropped. Although this approximation has been fruitful, its consequences have not been investigated. Here, this approximation is re-examined, and a new theory which includes often neglected small curvature terms is presented. Furthermore, the solution set of the new theory is explored for the unit disk domain and compared to the standard theory. Also, the analytical results are compared to experimental data.