comsol电磁超声导波检测 磁铁激励静磁场，线圈产生感应涡流，在1mm厚铝板中激励250kHz的lamb，在200mm位置处设置深0.8mm的裂纹缺陷，80mm处铝板表面的点探针接收波形如图4所示，依次为始波，裂纹反射S0模态，裂纹反射A0模态，端面反射S0模态。 模型编号:51＃
material = Model.Material.create('Aluminum_1mm'); material.propertyGroup('Def').set('youngs_modulus', '70e9[Pa]'); material.propertyGroup('Def').set('density', '2700[kg/m^3]'); material.propertyGroup('Def').set('poissons_ratio', 0.33);注意泊松比参数别手滑输成钢材的数值，之前有同行在这里翻过车。板厚直接体现在几何建模里，用拉伸特征做1mm薄板。永磁体建模要特别注意边界条件，这里有个取巧的方法：
magnet = Model.Geom.create('Magnet', 3); magnet.set('pos', [0 0 -5e-3]); //Z轴负方向偏移5mm Model.Physics.create('MagneticField', 'MagneticInsulation');线圈建模时得考虑趋肤效应，250kHz下铝板的趋肤深度大概1.2mm，刚好接近板厚。这里有个关键设置，线圈阻抗矩阵要勾选频域计算，不然瞬态分析会卡死：
coil = Model.Component.create('Coil', true); coil.geom.create('CoilGeometry', 'WorkPlane'); Model.Physics.create('CoilCurrent', 'Coil', coil);裂纹建模用布尔操作切个0.8mm深的V型槽。注意缺陷位置距激励源200mm，这个距离会影响模态分离。网格划分必须用边界层网格，特别是裂纹尖端：
mesh = Model.Mesh.create('CustomMesh'); mesh.automatic(true); mesh.set('physicsControlled', 'on'); mesh.set('elementSizeFactor', 0.7); //太密了算不动求解器设置要开瞬态+频域耦合。时间步长取1/(20*250kHz)=0.2μs，总时长根据声速计算，S0模态在铝板中速度约5km/s，跑完400mm需要80μs，但实际得留余量设到100μs。
接收点信号处理用带通滤波，中心频率250kHz，带宽±20%。特别注意模态分离——S0和A0的到时差是判断缺陷的关键。端面反射出现的位置能用来校准声速，之前有次仿真发现端面反射提前了3mm，结果是材料参数输错了。
波形图里的四个特征峰要对应到时域窗口：
- 始波直达信号（注意线圈耦合干扰）
- 裂纹反射S0（走时约40μs）
- 裂纹反射A0（速度慢，走时约65μs）
- 端面反射S0（走时约80μs）
最后吐槽下后处理，COMSOL导出数据到MATLAB做小波变换更靠谱，内置的FFT工具对瞬态信号不太友好。这个模型跑完大概需要32GB内存，建议关掉其他程序再提交计算。
