comsol静电梳状谐振器 MEMS加速度传感器仿真

梳齿结构在MEMS加速度计里就像精密钟表的齿轮，微小位移就能引发电容变化。最近在COMSOL里复现静电驱动谐振器时，发现梳齿间距的设置直接能让仿真结果"翻车"——不信邪的可以试试把初始间隙设成5μm以上，看谐振频率会不会像过山车一样乱跳。

建模第一步先画个丐版梳齿结构：

gap = 2e-6  # 梳齿间距2μm
finger_width = 4e-6
model.geom.create("geom1", 2)
model.geom("geom1").feature().create("r1", "Rectangle")
model.geom("geom1").feature("r1").set("size", [finger_width, 100e-6])

这段代码生成的单根梳齿宽度4μm、长度100μm，重点在于尺寸参数必须用科学计数法。有次手滑写成0.000004，结果边界条件直接报错——COMSOL对单位制的执念堪比强迫症。

静电耦合模块要特别注意非线性设置：

![静电物理场设置示意图]

comsol静电梳状谐振器 MEMS加速度传感器仿真

移动梳齿和固定梳齿需要分别指定电势，边界条件里的浮动电位选项就像给结构加了隐形的弹簧。这里有个坑：当加速度加载方向与梳齿运动方向存在夹角时，记得勾选空间电荷密度选项，否则静电力计算结果会漏掉关键分量。

调试时发现网格剖分是玄学重灾区。试过用自由四面体网格，结果电场分布像毕加索的抽象画。改成边界层网格后，梳齿边缘的场强计算终于正常：

# 边界层网格参数
boundary_layer_thickness = 0.1e-6
model.mesh("mesh1").feature().create("bl1", "BoundaryLayer")
model.mesh("mesh1").feature("bl1").set("thickness", boundary_layer_thickness)

这种设置能在不爆内存的前提下捕捉边缘效应，就像给梳齿套了层纳米级紧身衣。仿真跑完后别急着看位移云图，先检查电容变化曲线是否光滑——有次发现曲线出现锯齿状波动，排查半天原来是求解器步长设得比梳齿间距还大。

当施加1g加速度时，谐振频率偏移量达到312Hz。这个数值和工艺误差的关系堪比硬币的两面：实测中发现若锚点支撑梁的刚度偏差5%，频率偏移会直接缩水到200Hz以下。所以仿真时材料属性栏里的杨氏模量最好别用默认值，老老实实填单晶硅的169GPa才是正道。

参数化扫描时加速度从-10g扫到+10g，得到的灵敏度曲线居然不是直线——仔细看才发现大加速度下梳齿发生了pull-in效应。这时候该祭出非线性求解器了，就像给狂奔的野马套上缰绳。修改求解器配置后，曲线终于呈现出优美的线性段，斜率刚好对应设计的0.3mV/g灵敏度。

最后说个冷知识：COMSOL的后处理模块能直接导出电容梯度矩阵，这对标定交叉轴干扰比手动计算快十倍。不过要小心坐标系方向，有次忘记旋转坐标系导致交叉灵敏度计算全错，这教训值三杯浓缩咖啡的清醒时间。