热源强度分布函数
激光激发超声检测板材中缺陷 COMSOL有限元仿真模型，板状材料激光激发超声模拟，包括板状材料带有缺陷和无缺陷两种情况。
激光激发超声检测板材中缺陷，这活儿干起来真挺考验耐性的。以前总觉得做个有限元模型就是把参数填进去，等到跑 COMSOL 的时候才发现，光是网格和时间步长就能把人绕进去。这次重点是想看看板子里藏着的缺陷到底怎么在声波里显形，无缺陷和有缺陷两种情况对比着看，效果才直观。
物理机制上其实就是热弹性效应，激光能量打上去瞬间升温，材料膨胀产生应力波。在 COMSOL 里，这需要耦合热传导和结构力学。最麻烦的是那个热源怎么定义，毕竟激光不是均匀照在板子上的，是个高斯分布。
# t 为时间变量，用于模拟激光脉冲宽度 # x, y 为平面坐标，x0, y0 为光斑中心 # w 为光斑半径 power_density(t, x, y) = (2 * P / (pi * w^2)) * ... exp(-2 * ((x-x0)^2 + (y-y0)^2) / w^2) * ... sin(pi * t / pulse_width)这段逻辑对应到 COMSOL 的表达式编辑器里，其实就是个分段函数。关键是那个时间项，sin 波模拟短脉冲，如果pulse_width设得太宽，频谱就窄了，分辨率不够；设得太窄，高频噪声又太多，收敛困难。之前有一次忘了乘时间因子，直接把连续波当脉冲算，结果后面全是驻波，根本看不出缺陷反射。
几何建模上，无缺陷就是一个完整的长方体。要想模拟缺陷，通常是在板中间挖个坑。这里不建议用复杂的曲面，尽量用简单的几何特征，比如圆柱形气孔或者扁平的椭圆夹杂。因为后续网格要剖分得特别细，形状越怪，网格质量越难控制，尤其是缺陷周围，必须加密。
% 这里的脚本逻辑是为了自动生成缺陷模型 % 避免手动操作导致的位置偏差 defect_geom = boolean_difference('plate', 'void_hole'); mesh.param.set('size', 'maxh', 'lambda_min/3'); % lambda_min 是最高频对应的波长别看这几行伪代码，实际在 LiveLink for MATLAB 里跑的时候，每次改个参数都得重新生成网格，耗时很长。所以网格策略一般是：整体用较粗的网格保证计算效率，在激光照射点和缺陷周围手动划定区域进行局部加密。特别是厚度方向，因为是薄板，至少要保证每个波长上有 6-8 个单元，不然相位误差太大了。
边界条件处理也有讲究。四周要是自由的，上下表面也是自由的，模拟无限大板子的话，还得加完美匹配层（PML）吸收边界，防止波在边界反射回来干扰结果。不过对于这种短时探测，只要记录的时间窗口够短，波还没跑到边界就被截断了，这时候不加 PML 也能凑合用，省点内存。
两组模型跑完，后处理的数据对比很明显。无缺陷的情况下，接收探头收到的信号主要是直达波，后面跟着上下表面的多次反射，波形比较单一。
激光激发超声检测板材中缺陷 COMSOL有限元仿真模型，板状材料激光激发超声模拟，包括板状材料带有缺陷和无缺陷两种情况。
但是加上缺陷之后，情况就变了。
# 提取特定时间窗口的信号进行分析 signal_with_defect = data.time_series[t_start:t_end] reflection_peak = find_local_max(signal_with_defect, window=10ns) print(f"缺陷回波时间差：{reflection_peak - direct_wave_time} ns")代码里这种简单的峰值查找逻辑，用来定位缺陷回波位置很方便。你看那个reflection_peak，它出现在直达波之后，底面回波之前。这个时间差乘以声速，除以 2，大概就是缺陷的深度。
有意思的是，如果缺陷比较小，反射波的能量会明显减弱。在 COMSOL 的 2D 绘图里，你可以把位移场做成动画，能看到波阵面扫过缺陷区域时，能量发生了散射，一部分绕过去了，一部分折返回来了。这就是为什么有时候仿真里看不到明显的回波峰——可能是缺陷太小，散射成了衍射波，能量分散了。
还有一个坑是数值色散。如果你的网格不够密，或者时间步长太大，算出来的波速会比理论值快。这时候即使加了缺陷，反射波的位置也会漂移。所以跑完结果第一件事，就是先看无缺陷模型的波速对不对得上材料的声速参数。如果不对，那就是离散误差惹的祸，得回去加密网格。
另外，激光激发源的位置也很敏感。如果打在缺陷正上方，反射信号最强；稍微偏一点，信号就弱下去了。这在仿真里很容易通过扫描光源位置来观察，虽然实际检测不一定这么扫，但仿真可以帮你找到最佳的检测布局。
总的来说，这套模型的核心不在于方程有多复杂，而在于参数的敏感度分析。从热源功率到网格密度，再到缺陷的尺寸，每一个变量都会影响最终波形。有时候为了验证算法，故意把缺陷设得很大，波形一目了然；有时候为了贴近工程，把缺陷设得很隐蔽，信噪比就很考验滤波算法了。
这玩意儿跑通了也就那样，但想到以后能靠这个优化无损检测的参数，还是觉得值得折腾。下次准备试试不同材料属性下的波传播，或者引入各向异性，估计又有新的坑等着填。