Comsol脉冲涡流无损检测仿真 图一：脉冲涡流仿真，检出电压信号 图二：脉冲涡流模型 图三：磁通密度模 图四：磁通密度模
线圈贴着金属板扫过，检测仪突然滴滴作响——这种非接触式探伤场景在工业检测中司空见惯。今天咱们用COMSOL扒开这个过程的底层物理，看看电磁场如何在毫秒间完成"金属体检"。
打开COMSOL新建模型，先在几何菜单画个直径20mm的圆形线圈，厚度设为0.5mm（毕竟真实线圈都有横截面）。下方的金属板用矩形表示，尺寸建议设为200mm×200mm×5mm，足够大的尺寸能避免边界效应干扰。别忘了在两者之间留0.1mm空气隙，模拟真实探头的提离效应。
材料参数设置有个坑要注意：铝板的电导率别直接填3.5e7 S/m，实战中材料参数往往需要校准。这里咱们用参数化设置：
sigma = 3.5e7 + (defect_depth>0)*1e6; //缺陷区域电导率变化这段脚本实现了缺陷区域的电导率突变，defect_depth参数控制缺陷深度。当检测到缺陷时自动调整材料属性，比GUI操作更灵活。
脉冲激励的设置是核心难点。在电流源设置里，用分段函数构造梯形脉冲：
double[] tPoints = {0, 1e-6, 9e-6, 10e-6}; double[] IValues = {0, 100, 100, 0}; source.set('I0', new Interpolation(tPoints, IValues));这个10微秒脉宽的梯形波，前沿和后沿各留1微秒的过渡，避免数值震荡。注意时间步长要小于脉冲上升时间的1/5，建议设置为0.1微秒。
Comsol脉冲涡流无损检测仿真 图一：脉冲涡流仿真，检出电压信号 图二：脉冲涡流模型 图三：磁通密度模 图四：磁通密度模
求解器配置时记得勾选"初始值一致性检查"，瞬态求解器选BDF并设置最大阶数为2。有个小技巧：在因变量设置里勾选"磁矢势规范变换"，能显著加快计算速度，特别是处理三维涡流问题时效果明显。
图二的磁通密度模分布图（图三、图四）显示，在脉冲下降沿时刻（t=9.5μs），缺陷边缘出现明显的磁场畸变。这是因为涡流在缺陷处被迫改变路径，形成局部的磁场集中。代码层面，可以通过后处理提取特定时刻的场量：
mph.export(model, '磁场数据', 't', [9.5e-6], 'data', 'mfnc')这种磁场畸变会在线圈中感应出特征电压信号（图一）。仔细观察电压曲线的二次衰减过程：第一个峰对应脉冲关断的瞬态，后续的缓变衰减反映缺陷引起的涡流扰动。用峰值比算法可以量化缺陷深度：
Vpeak2/Vpeak1 = 0.32*exp(-defect_depth/0.7)
这个经验公式是前人通过大量仿真总结的，仿真时注意调整探头提离距离参数做验证。
最后提个醒：网格剖分在缺陷区域需要局部加密，建议使用边界层网格。但别超过5层边界层，否则计算量会爆炸。实在拿不准的话，用自适应网格先试算一次，观察哪里需要加密再手动调整。
