1. Introduction
In the past few decades, polymer membranes have been widely used to manufacture vibration sensors owing to their mechanical flexibility, high sensitivity, cost-effectiveness, and facile integration into electronic circuits. These sensors, comprising a flexible membrane, detect changes in the surrounding environment based on changes in the vibration characteristics. Owing to the high sensitivity of membrane vibrations to changes in the surrounding environment, these sensors can detect small changes in temperature, pressure, and other physical parameters. As a result, polymer-membrane vibration sensors are applied in numerous fields, including structural health monitoring, the automotive industry, and medical diagnostics [
1,
2,
3,
4,
5].
Understanding the vibration characteristics of membranes is required to extend the applications of membrane sensors. Several studies have documented the natural frequencies of homogeneous membranes [
6,
7]. Generally, exact analytical solutions are limited to simple geometries such as rectangles, circles, sectors, ellipses, and isosceles right triangles [
8,
9,
10]. Various approximate approaches have been developed to solve problems involving polygonal or compound geometries with Dirichlet, Neumann, or mixed boundary conditions [
11,
12,
13,
14,
15]. As alternative practical routes, numerical schemes such as the finite-difference, finite-element (FE), boundary-element (BE), and mesh-free methods have been reported [
16,
17,
18,
19,
20,
21,
22,
23,
24]. Recently, novel numerical methods, such as the Galerkin [
25] and Bezier [
26] methods, have been shown to attain higher stability and accuracy than other numerical methods.
Among previous research topics, homogeneous L-shaped membranes are apparently favored, especially symmetric ones [
27,
28]. Associated eigenvalues are often used as validation criteria or cross-reference sources [
29,
30,
31,
32,
33,
34,
35,
36,
37,
38,
39]. Unfortunately, the theoretical solution for the corresponding free vibration problems is seemingly difficult. The most commonly used method (i.e., the variable separation method) is only feasible for simple geometric membranes whose boundaries fit perfectly into a particular separable coordinate system, such as circular and elliptical membranes [
8]. Although the geometries of L-shaped membranes inherently consist of rectangles, exact analytical solutions for symmetric and asymmetric cases are not readily available. This is because the contours of L-shaped geometries only partially conform to a Cartesian rectangular coordinate system.
As frequency-tunable membrane sensors demonstrate more comprehensive applications (e.g., biomedical sensing, environmental monitoring, and industrial process control), a conceptual model for a tunable microfabricated device comprising asymmetric L-shaped membranes is proposed herein. First, a Fourier series solution is derived, which can be used as an alternative to an exact analytical solution. The region-matching technique combines the domain decomposition method and the variable separation method. When the appropriate auxiliary boundary is selected, the eigenfunctions in each subregion inherently satisfy most parts of the clamped boundary conditions. This semi-analytical nature permits relatively fast convergence and relatively high accuracy of the present results compared to those obtained using commercial FE software Abaqus [
40].
This serves as the cornerstone of our investigation into frequency-tunable sensor design. The operating frequency of these membrane-based sensors can be tuned by either material selection or membrane geometry. The semi-analytical solution presented above is used to determine a suitable material for a membrane vibration sensor with a specific frequency requirement. In addition, frequency tuning is demonstrated by varying the length or width of the membrane segments.
2. Materials and Methods
Consider a tunable microfabricated device comprising a stretched asymmetric L-shaped membrane clamped on all edges.
Figure 1 shows the problem geometry. The thin elastic membrane is assumed to be homogeneous and is characterized by four parameters,
b,
d,
a, and
h, where
b and
d represent the length and width of the left-hand segment, and
a and
h represent those of the right-hand segment, respectively. The uniform tensile force per unit length is
T, and the constant mass per unit area is
ρ. The origin of the Cartesian coordinate system (
x,
y) is the vertical projection point of the concave corner on the underside of the membrane.
Based on the domain decomposition method, a vertical auxiliary boundary
Sa is introduced to divide the entire computational domain into two enclosed regions: regions 1 and 2 (see
Figure 1). The small transverse motions of the membrane, denoted by
uj, must obey the governing Helmholtz equations:
where
is the two-dimensional Laplacian operator in the
x–y plane,
is the wavenumber, and
ω is the angular natural frequency of vibration. The subscript
j, where
j = 1 and 2, represents the region number. The time-harmonic factor is understood throughout this section.
The zero-displacement boundary conditions are imposed along the edge of region 1,
and along the edge of region 2,
By applying the method of eigenfunction expansion, the displacement fields in regions 1 and 2, respectively, can be expressed as follows:
with
where the expansion coefficients
An and
Bn are unknown. Notably, Equations (9) and (10) inherently satisfy the governing Equation (1) and most of the boundary conditions around the membrane edge, except those on
Sa.
Enforcing the displacement continuity condition on
Sa yields
By multiplying Equation (13) by a sequence of sine functions and integrating over the appropriate intervals,
According to Equation (5), the result of integration from
h to
b vanishes. This implies that the upper limit of integration on the left-hand side of Equation (14) can be extended from
h to
b. Therefore, the orthogonal property of sine functions can be applied to Equation (14) directly. Consequently, the following relation holds:
where
Similarly, by considering the slope continuity condition across
Sa,
and applying successive sine functions, and integrating over the range [0,
h],
By exploiting Equation (15) to eliminate the unknown coefficients
An and rearranging the results in a system of linear algebraic equations (with unknown coefficients
Bn), the following matrix form can be obtained:
with
where
δi,j is the Kronecker delta function.
Clearly, Equation (19) constitutes a generalized matrix eigenvalue problem, as expected. Standard techniques can thus be used to evaluate the natural frequencies of the present membrane, which are related to the roots of the determinant equation given by
Once the eigenvalues k are found, the expansion coefficients Bn can be evaluated using partitioned matrices and block multiplication (cf. Equation (19)). Therefore, the expansion coefficients An are determined directly from Equation (15). Eventually, the natural modes (eigenmodes) can be obtained by Equations (9) and (10).
In Equation (20), the summation indices n and the weighting indices i and j are truncated after N terms. Hence, Equation (19) constitutes a system of N equations with N + 1 unknowns. The number of truncation terms in consideration depends only on the accuracy requirement.
5. Conclusions
In this study, a tunable vibration-based sensor comprising asymmetric L-shaped membranes was proposed for use in MEMS applications. A semi-analytical series solution was derived for the corresponding two-dimensional Helmholtz eigenvalue problem to better understand the vibration characteristics of these membranes. The region-matching technique, which combines domain decomposition and variable separation methods, was used to construct the Cartesian displacement fields. The results of this semi-analytical analysis were comparable to those of finite-element simulations over a wide range of asymmetric L-shaped geometries, as well as previous numerical (approximation) solutions for symmetric cases. These results demonstrate the good performance and computational efficiency of the semi-analytical approach.
The asymmetric L-shaped membrane sensor is advantageous as the operating frequency can be tuned by adjusting the geometry, such as the length or width of the membrane segments. The fundamental natural frequency increases with a decrease in the length or width of the membrane segments, providing a means of frequency tuning. In addition, the proposed model can be used to determine the membrane material for membrane sensors with specific frequency requirements and geometries. Based on the changes in resonant frequency due to interactions between the ambient environment and the membrane surface, these membranes could be employed as gas sensors, chemical sensors, or mass sensors. Our numerical results indicate that the performance sensitivity for mass sensing is significantly affected by the geometry of the asymmetric L-shaped membrane. The asymmetry of the membrane can thus be customized for different mass sensing sensitivity requirements.
It is worth noting that the proposed semi-analytical solutions are derived on the basis of the linear, homogeneous, and isotropic material properties of membranes subjected to a uniform tensile stress. In addition, the pre-applied tension and the material properties of the membrane are assumed to remain constant during the vibration, neglecting environmental factors, such as temperature and pressure changes. Further applications are possible if the membrane-based sensors comprise multi-segmented rectangles of different material types. In addition, fabrication issues in micromachined tunable devices will be investigated in the future.