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