基于无量纲大挠度控制方程的MEMS电容式压力传感器归一化灵敏度理论与数值研究

《Micromachines》：Theoretical and Numerical Investigation of Normalized Sensitivity in MEMS Capacitive Pressure Sensors Based on Dimensionless Large-Deflection Governing Equations

　　

MEMS电容式压力传感器（MEMS capacitive pressure sensors）是压力传感器的主要类别之一。然而，在非接触工作模式下，其灵敏度通常随施加压力非线性变化，这增加了设计阶段多压力点灵敏度评估的难度和工作量。在本研究中，基于大挠度控制方程（large-deflection governing equations）建立了MEMS电容式压力传感器的归一化灵敏度分析框架。在该框架内，归一化灵敏度与归一化压力之间的非线性关系由五个无量纲参数表征。结果表明，当这些无量纲参数相同时，具有不同结构和材料参数组合的传感器会收敛到几乎相同的归一化灵敏度-压力特性曲线，验证案例中的平均相对偏差不超过0.05%。此外，进一步分析了每个无量纲参数对特性曲线非线性行为的影响。所提出的框架为MEMS电容式压力传感器的灵敏度比较、性能评估和结构优化提供了统一且系统的分析基础。

论文解读：基于无量纲大挠度控制方程的MEMS电容式压力传感器归一化灵敏度研究

**研究背景与问题**

MEMS（微机电系统）压力传感器是广泛应用的微型化压力测量器件，其中电容式压力传感器（capacitive pressure sensors）因其低功耗、高灵敏度、良好的温度稳定性及适合批量制造等优点，在工业控制、汽车电子和航空航天等领域得到大量研究与应用。其基本工作原理为：压力敏感膜片在压力载荷下弯曲，改变电极间隙，从而将压力信号转换为电容信号。灵敏度是这类传感器最重要的性能指标之一，需要在设计阶段进行量化以指导结构尺寸和材料参数的选择。

然而，现有研究在评估灵敏度时存在显著问题：首先，传感器在非接触工作模式下，灵敏度随施加压力非线性变化，这导致在设计阶段需要在多个压力点进行逐点电容和灵敏度计算，增加了迭代设计的工作量；其次，对于大挠度工况及膜片几何形状（如矩形、圆形）等因素对归一化灵敏度的影响，尚缺乏系统分析和统一表征。先前的研究（如Rochus等、Jindal等、Xu等）主要基于小挠度假设（small-deflection assumption）的圆形膜片，但未能涵盖大挠度效应及膜片几何形状的系统影响。

因此，研究人员开展本研究，旨在建立基于大挠度控制方程的归一化灵敏度分析框架，通过无量纲参数统一表征不同设计参数组合下传感器的灵敏度特性，以简化比较、评估和优化过程。该论文发表在《Micromachines》。

**主要关键技术方法**

研究人员采用以下关键技术方法：（1）大挠度控制方程的无量纲化处理：将膜片大挠度控制方程（包括残余应力项）通过选取特征尺度（坐标、挠度、压力、残余应力、应力函数）进行无量纲化，得到控制无量纲参数（ν、AR、σ*、λ、γ）的方程。（2）归一化流程：选取接触压力（contact pressure, pc）和初始灵敏度（initial sensitivity, S0）作为基准，定义归一化压力pn和归一化灵敏度Sn，推导出Sn-pn关系由无量纲参数决定。（3）间接计算方法：基于无量纲参数一致性，通过有限元仿真（FEA）计算代表性传感器的Sn-pn曲线，避免逐点计算。（4）有限元仿真设置：采用COMSOL固体力学有限元公式，考虑几何非线性，对膜片进行网格划分，通过网格独立性验证（接触压力偏差<0.5%），并对电容进行面积积分。所有仿真基于膜片自身，未涉及具体试剂或质粒构建，亦无样本队列来源。

**研究结果**

**3.1. 归一化灵敏度特性相似性验证**

通过选择两对传感器（圆形膜片传感器1和2，γ=0.9647, λ=0.7, ν=0.28, σ*=-0.2；矩形膜片传感器3和4，γ=0.9474, AR=2, λ=1.8, ν=0.28, σ*=-0.5），采用间接计算方法计算Sn-pn曲线。结果表明：圆形膜片传感器的平均相对偏差仅为0.01%，最大相对偏差0.04%（pn=0.99）；矩形膜片传感器的平均相对偏差为0.05%，最大相对偏差0.13%（pn=0.99）。这验证了当无量纲参数相同时，不同尺寸参数传感器具有相似的归一化灵敏度特性曲线。

**3.2. 关键无量纲参数的影响**

**3.2.1. 间隙厚度比λ（Gap-to-Thickness Ratio）**

对圆形膜片传感器，λ从0.1到2.0变化，固定ν=0.28, σ*=0, γ=0.96。结果显示：λ=0.1和0.2时的Sn-pn曲线几乎重合（最大相对偏差3.5%），属于小挠度工况；λ>0.2时，曲线逐渐偏离小挠度参考曲线，且随pn增大偏离加剧。λ=0.1时pn=0.99的Sn约为λ=2.0时的4.9倍。这表明λ越大，Sn-pn曲线的非线性程度越低，因为大挠度下w₀*与p*呈三次函数关系，单位压力增量引起的挠度增量减小。

**3.2.2. 纵横比AR（Aspect Ratio）**

对矩形膜片（AR=1,2,3）和圆形膜片，在γ=0.96, ν=0.28, σ*=0条件下，分别计算λ=0.1和λ=2.0时的Sn-pn曲线。结果表明：AR=1的矩形膜片与圆形膜片曲线几乎重叠（最大相对偏差<2.0%）；随着AR增大，曲线偏离并高于AR=1的曲线，且pn=0.99时AR=3的Sn约为AR=1的1.8倍。这表明增大AR增强Sn-pn曲线的非线性，原因在于大AR的膜片形状函数W*(x*, y*)使更多区域挠度接近最大挠度，导致局部间隙倒数变化更显著。

**3.2.3. 泊松比ν（Poisson’s Ratio）**

比较圆形膜片传感器在ν=0.22和ν=0.28下的Sn-pn曲线（σ*=0, γ=0.96, λ=0.1,0.5,2.0）。结果显示：两条曲线差异不显著，最大偏差发生在λ=2.0、pn=0.99处，仅3.7%。这是因为ν以1-ν²形式进入方程，ν从0.28变为0.22时，1-ν²仅变化约3.3%，因此对Sn-pn曲线非线性影响较小。

**3.2.4. 无量纲残余应力σ*（Dimensionless Residual Stress）**

对圆形膜片传感器，在γ=0.96, ν=0.28下，比较σ*=0.50, 0.25, 0, -0.25, -0.50时的Sn-pn曲线。λ=0.1时各曲线接近重叠，但pn=0.99时σ*=-0.50和σ*=0.50相对于σ*=0的偏差分别为4.0%和3.3%；λ=0.5时，最大相对偏差均超过10%（pn=0.99）。总体而言，压应力（负σ*）降低非线性，拉应力（正σ*）增加非线性，且绝对值越大影响越显著。原因在于压应力使初始灵敏度S₀增大，归一化后Sn降低。

**3.2.5. 等效间隙系数γ（Equivalent Gap Coefficient）**

对圆形膜片传感器，在ν=0.28, σ*=0下，比较γ=0.92, 0.94, 0.96, 0.98时不同λ（0.1,0.5,1.0,2.0）的Sn-pn曲线。结果表明：给定pn下，Sn随γ增大而增大，且γ越大曲线非线性越强；随着λ增大，不同γ曲线之间的相对偏差提前出现（如λ=0.1时pn=0.7超过10%，λ=2.0时pn=0.4超过10%）。这可通过局部归一化灵敏度表达式解释，γ增大导致局部灵敏度因子增大，而λ增大放大了这一效应。

**讨论与结论**

讨论部分指出，与先前仅适用于小挠度圆形膜片的模型相比，本框架覆盖小挠度和大挠度工况，明确纳入λ表征膜片变形状态，并考虑了AR、ν、σ*和γ的影响。γ的影响在不同λ下得到系统分析。该框架使设计人员能通过无量纲参数快速评估和比较不同传感器非线性趋势，避免逐点计算。此外，可通过建立归一化Sn-pn曲线族，仅需确定pc和S0即可通过查表获取全范围灵敏度，提高设计效率。

研究结论（翻译）：
在本研究中，基于大挠度控制方程系统地研究了非接触工作模式下电容式压力传感器的非线性灵敏度特性。建立了电容式压力传感器输出的归一化分析模型。证明当无量纲参数λ、AR、ν、σ*和γ相同时，具有不同结构和材料参数组合的电容式压力传感器表现出相似的Sn-pn曲线。通过间接计算方法进一步验证，具有不同尺寸参数但相同无量纲参数的传感器的Sn-pn曲线平均相对偏差不超过0.05%。还揭示了无量纲参数对归一化灵敏度非线性特性的影响。其中，增大γ或AR倾向于增加Sn-pn曲线的非线性程度，而增大λ则倾向于降低之。在所研究的预屈曲范围内，压应力降低非线性程度，而拉应力增加之；两种情况下，影响随无量纲残余应力幅值增大而更显著。泊松比ν的影响相对有限。总体而言，所提出的无量纲框架为比较、预测和定制非接触电容式压力传感器的非线性灵敏度特性提供了一种系统且统一的方法，从而促进传感器设计评估与优化。