COMSOL仿真单一模态兰姆波：超声兰姆波检测电路板缺陷的有限元仿真研究 版本为5.6，低于5.6的版本打不开此模型
在电路板的无损检测中，超声兰姆波的模态选择直接影响缺陷识别精度。最近用COMSOL 5.6复现了一个经典案例：通过单一模态S0波的传播特性检测PCB内部裂纹。这个模型对网格划分要求极其变态——尤其是在薄板结构的边界层处理上。
COMSOL仿真单一模态兰姆波：超声兰姆波检测电路板缺陷的有限元仿真研究 版本为5.6，低于5.6的版本打不开此模型
先看几何建模部分。铝基板厚度设定为1mm，裂纹用0.05mm宽的矩形缺口模拟。这里有个坑：如果直接用布尔运算做减法，会导致相邻网格单元出现畸形。改用"层叠法"构建几何更稳妥：
model.geom.create('comp1', 'Component'); model.geom('comp1').create('r1', 'Rectangle'); model.geom('comp1').feature('r1').set('size', [20 1]); //基板尺寸 model.geom('comp1').create('r2', 'Rectangle'); model.geom('comp1').feature('r2').set('size', [0.05 0.3]); //裂纹尺寸 model.geom('comp1').create('dif1', 'Difference'); model.geom('comp1').feature('dif1').selection('input').set({'r1'}); model.geom('comp1').feature('dif1').selection('input2').set({'r2'});物理场设置阶段，固体力学接口需要特别关注频散方程。这里采用特征频率分析锁定S0模态，关键参数是约化频率kr的取值。实测发现当kr=0.05时，特征频率在1.2MHz附近最稳定：
model.study.create('std1'); model.study('std1').create('freq', 'Frequency'); model.study('std1').feature('freq').set('plist', 'linspace(1e6,2e6,50)'); model.study('std1').feature('freq').set('punit', 'Hz');网格划分是翻车重灾区。薄板结构必须采用边界层网格，特别是裂纹尖端需要二次细化。这里用映射网格控制y方向层数：
model.mesh.create('mesh1', 'geom1'); model.mesh('mesh1').create('map1', 'Map'); model.mesh('mesh1').feature('map1').create('dis1', 'Distribution'); model.mesh('mesh1').feature('map1').feature('dis1').set('numelem', 8); //厚度方向8层 model.mesh('mesh1').feature('map1').set('method', 'specify');后处理阶段，用声压级云图观察波传播时，发现直接输出位移场会丢失相位信息。改用实部虚部分量合成复数位移后，成功捕捉到裂纹引起的波包畸变：
model.result.dataset.create('dset1', 'Solution'); model.result.create('pg1', 'PlotGroup'); model.result('pg1').create('surf1', 'Surface'); model.result('pg1').feature('surf1').set('expr', 'sqrt(emw.ReEx^2 + emw.ImEx^2)'); //合成电场模值这个模型跑起来需要至少32G内存，特别是瞬态分析时建议开启集群计算。遇到过最诡异的bug是求解器自动选择错误导致模态混叠——强制指定PARDISO求解器后问题消失。最后提醒：5.6版本新增的薄层阻尼功能对能量衰减计算帮助很大，别忘记在材料属性里勾选损耗因子选项。
