电磁超声检测（EMAT）在工业无损检测领域一直是个热门方向，最近在COMSOL 5.6上折腾了个L型铝板裂纹检测模型，整个过程像极了在虚拟实验室里玩"找茬"游戏。模型文件必须用5.6版本打开这点要注意——之前用5.5死活加载不了，后来发现是新增的磁致伸缩效应模块搞的鬼。
先说说模型的核心配置：铝板厚度3mm，L型转角处预设了0.5mm深的表面裂纹。激发线圈采用8字型缠绕结构，这种设计能让涡流场在裂纹位置产生明显的扰动。材料参数设置里有个小技巧，在铝的磁导率参数栏直接输入if(emw.normE>1e3, 1.05, 1)这种条件表达式，可以模拟高频磁场下的非线性效应。
% 线圈激励参数设置示例
f0 = 2e6; % 2MHz激励频率
tau = 0.5/f0; % 脉宽调节
model.param.set('f0', [num2str(f0) '[Hz]']);
model.component('comp1').physics('emw').feature('ew1').set('WaveType', 'Harmonic');
物理场耦合要特别注意多物理场接口的顺序——先电磁场后固体力学场的设置顺序能让声波传播计算更稳定。网格划分阶段建议在裂纹尖端采用边界层网格，这里有个取巧的方法：用自适应网格生成后，手动在裂纹区域添加圆柱形网格细化。
求解器配置经常让人抓狂，推荐尝试分离式求解器搭配GMRES迭代法。特别是当计算发散时，试试把相对容差从1e-6调到1e-5，往往能救命。后处理阶段有个好玩的功能：在声压场图上叠加涡流密度流线图，裂纹位置的流线扭曲会非常明显。
时域信号的提取要注意采样间隔设置，这里有个实测有效的代码片段：
t = linspace(0, 50e-6, 1000);
sensor_data = mphinterp(model, 'espr.V', 'coord', [0.02;0.03;0.001], 't', t);
plot(t*1e6, real(sensor_data));
xlabel('Time (\mu s)');
ylabel('Signal (V)');
最后给个实用建议：当声波传播路径存在L型拐角时，在转角处添加人工阻尼层能有效抑制数值震荡。这个模型的妙处在于可以通过修改几何参数快速适配不同角度的折弯件检测，下次准备试试把铝板换成碳纤维复合材料看看效果如何。
