Closed Form Representations of Some Series in Darling ’ s Model for Squeeze Film Damping with a Rectangular Plate

Compressible squeeze film damping is a phenomenon of great importance for micromachines. For example, for the optimal design of an electrostatically actuated micro-cantilever mass sensor that operates in air, it is essential to have a model for the system behavior that can be evaluated efficiently. An analytical model that is based upon a solution of the linearized Reynolds equation has been given by R.B. Darling. In this paper we explain how some infinite sums that appear in Darling’s model can be evaluated analytically. As an example of applications of these closed form representations, we compute an approximation for the critical frequency where the spring component of the reaction force on the microplate, due to the motion through the air, is equal to a certain given multiple of the damping component. We also show how some double series that appear in the model can be reduced to a single infinite series that can be approximated efficiently.


Introduction
An analytical model for squeeze film damping is presented by Darling et al. in [2] for different venting conditions.Based upon the solution of the linearized Reynolds equations, series representations of the resulting reaction forces are presented for different air venting boundary OPEN ACCESS conditions; in particular, for ideally vented rectangular vibrating plates.In [1], based upon the model of Darling et al. [2], an optimal sensor design problem is studied.
For optimal sensor design problems, the aim is to find system parameters for which a certain real valued function that measures the sensor's performance is maximized.In order to determine such parameters that solve the optimal design problems, it is useful to have a compact analytical model for the system behavior.
In [1], two kinds of boundary conditions appear for rectangular plates: Plates with two opposite edges closed and two opposite edges venting and plates with two adjacent edges venting and two adjacent edges closed.For the former, we give a closed form solution of the corresponding spring and damping constants that are given as infinite series in [1] and [2].This model has also been used in [3] to model squeeze film effects.For the latter case, we present an approximation of the resulting reaction forces solution that can be evaluated efficiently.We also give a bound for the approximation error.
Let us start with a short presentation of the spring and damping constants for a rectangular plate with two opposite edges venting.Let b denote the plate dimension in the direction of venting and let A denote the plate area.P A is the ambient pressure and g 0 is the nominal gap.Let μ denote the viscosity of the trapped gas.Define k n = nπ/b and the constant For the vibration frequency ω, according to Darling's model we have the spring constant In the sequel, we will give compact representations for k and β that do not require the evaluation of infinite series.In order to do this, we introduce some auxiliary functions in the next section.

Closed form Representations for Some Infinite Series
In this section we define some auxiliary functions that appear in Darling's model as infinite series.In Section 4, we will show how stiffness and damping can be expressed using these series.For these series we present closed form representations that can be evaluated efficiently and thus are useful for the solution of the optimal design problem.For a real number x let F(x) = 1/(1 + x 4 ) + 1/(3 4 + x 4 ) + 1/(5 4 + x 4 ) + 1/(7 4 + x 4 ) + ...
then for all x unequal to zero we have the representation (with w(x)= π x/2 1/2 ) and F(0)= π 4 /96, which is also the maximal value attained by F.
As a second real function, define then we have for all x unequal to zero and B(0) = π⁴/960, which is also the maximal value attained by B. Note that Equation (4) implies that x⁴B(x) converges to the limit π²/8 as x tends to infinity.
For a complex number z with −z 2 not equal to the square of any odd integer, let then for all z unequal to zero we have the representation and H(0) = F(0) = π⁴/96.

Efficient Approximation of a Double Series
For certain boundary conditions, in Darling's model infinite double series appear.We give approximations for these double series that can be evaluated efficiently.We also provide bounds for the approximation error.Let the parameter λ greater than or equal to 1 be given.
For a complex number z, let D denote the double series For k = 7, we obtain |R(z,7)| < 2 × 10 −5 /λ 3 .Thus for k = 7 and purely imaginary z 2 , A approximates D with a uniform bound that is less than 2 × 10 −5 .

Computation of the Spring and Damping Constants: Two Opposite Edges Venting
The spring and damping constants for a plate with two opposite edges venting can be expressed in terms of the functions F and B that have been introduced in Section 2 in Equations ( 2) and ( 4).According to Darling's model, we have (see Li and Miller [1]) the spring constant k = (8ω²AP A /(π²g 0 ))(α⁴b⁴/π 4 )B(αbω 1/2 /π) (12) with the notation introduced in Section 1.The dimension in Equation ( 12) is N/m.The damping constant β is given by The dimension in Equation ( 13) is N s/m.Using Equations ( 4) and ( 2), these functions can easily be implemented for numerical evaluation.
For a sensor design, Li and Miller [1] consider a square plate of side length p with a square hole of side length s in the middle.This plate is divided into eight elements.Four plates with opposite edges venting appear as side elements.The other four elements are square corner plates with adjacent edges venting.The elements with opposite venting are described by Equations ( 12) and ( 13), where the value of b is given by (p−s)/2.The optimal design problem consists in determining an optimal value for s.For the corner elements with adjacent edges venting the side length is (p−s)/2.The corresponding spring and damping constants are given by double series that are more difficult to evaluate analytically.
We will give a suitable approximation in Section 5.For the elements with opposite venting described by Equations ( 12) and (13) according to the Equation (46) in reference [3], the reaction force on the plate due to motion through the air is where C 0 is a suitably chosen real constant, j 2 = −1 and H can be computed using Equation (6).

An Application: Approximation of the Critical Frequency where k = γ β
Our representations allow the approximate computation of the frequency for which k = γ β, where γ is a given constant with the appropriate dimensions.At this frequency, the ratio k/ β of viscous force and spring force is equal to the given value γ.According to Equations ( 12) and ( 13), this problem is equivalent to the solution of the equation with x = αbω 1/2 /π and the dimensionless constant C = γα²b²/π².Equation ( 15) is equivalent to

Conclusions
The air between two parallel microplates can compress to store energy or vent to dampen energy.This compressible squeeze film damping is a phenomenon of great importance for micromachines.In reference [2], Darling et al have presented an analytical model for this phenomenon that is based upon expansions in series of eigenfunctions.To apply this model for optimization purposes, it is useful to have closed form representations of the corresponding functions.In this paper we give such a representation for the case that two opposite edges of the plate are closed.The case where two adjacent edges of the plate are closed, leads to functions given by double series.For this double series we provide an approximation that can be computed easily.Moreover, we give a bound for the corresponding approximation error.