表面裂纹大量地存在于工程构件当中，如何快速准确地识别它们已成为无损检测的任务。普通的检测方法不适用于高温、高压、有毒^[1]等恶劣环境，而激光超声无损检测凭借其非破坏、非接触^[2-3]等优点克服了这种局限性，并且具有较高的精度，成为了无损检测的一个重要研究方向^[4]。近些年来，随着计算机技术的飞速发展，计算机数值仿真有了越来越广泛的应用。关于激光超声无损检测的数值模拟已成为无损检测领域的一个热点^[5-6]，而瑞利波进行无损检测的应用也越来越广^[7-8]。通过模拟激光激发超声波的物理过程与超声波的传播过程，可以实现对表面裂纹的检测，从而对实验研究提供一种判定方法。

本文中将激光冲击作用进行简化，将激光的简化模型与激光的参量联系起来，利用ABAQUS软件模拟了激光激发出超声波的物理过程和超声波在材料表面进行传播的过程，并借此模型建立了表面裂纹深度与表面波信号特征量之间的联系，为激光超声检测表面裂纹深度提供了一种有效的方法。

建立一块各向同性且均匀的铝制薄板，铝板的长度为100mm，宽度为30mm，厚度为2mm，裂纹左边界距离板左端75.93mm。现有一束脉冲线源激光沿厚度方向照射到试样表面上，该问题可以简化成平面应变问题来分析, 如图 1所示。

Figure 1.  Schematic diagram of model

为了精确模拟激光激发超声波的物理过程和超声波在表面的传播过程，网格的尺寸应不超过波长的1/4^[9]。又考虑到计算效率，将网格尺寸取为0.040mm×0.013mm。在模型的左边界添加全约束。有限单元取平面应变单元CPE4R。用有限元来模拟超声波，模型必须要有足够大的几何尺寸才能避免来自边界的反射波对计算的影响，但这样会消耗大量的计算资源和计算时间。为了提高计算效率，同时也为了消除边界回波对表面波信号的干扰，在模型的左边界、右边界以及下边界使用无限单元，无限单元为CINPE4。超声波若在物体中能够完整的传播，需要有足够的时间，因此总时间长度取为5×10^-5s，同时保证激光加载的精度，取时间步长为1×10^-9s。

为了使激光冲击作用的简化更接近于实际情况，体现激光在激发过程中的作用，需要将激光冲击作用与力作用相互联系。本文中考虑了激光的主要参量，即脉冲宽度、波长、能量、功率密度，从峰值功率密度、时间分布、空间分布3个角度将激光作用简化成力的作用。

激光波长λ=1064nm，脉冲宽度τ=10ns，激光作用面宽度d=0.6mm，单次脉冲能量E₀=13.5mJ，峰值功率密度${I_0} = \frac{{{E_0}}}{{\tau A}} = 7.5 \times {10^6}{\text{W/c}}{{\text{m}}^{\text{2}}} $。其中，A为面积，线源脉冲激光在作用面上产生的峰值压力为^[10]：

式中，k为与材料有关的量，取k=2.3×10^-7。

激光脉冲的空间分布函数^[11]为：

式中, r₀是激光辐照的半径, r为距激光作用面中心线的距离。

激光脉冲的时间分布函数^[12]为：

式中，α为可由峰值压力确定的修正系数，取α=21，且载荷作用时间t取为脉冲宽度的3倍^[13], 即载荷作用时间为30ns。

由于之后的研究是建立在瑞利表面波基础上的，所以有必要验证在模型上加载的脉冲是否能够产生瑞利表面波。一般瑞利表面波的验证是根据其性质来判断的。

首先在模型表面上离激励处较远的区域取3个点，得到如图 2所示的位移时间曲线。

Figure 2.  Displacement-time curves of nodes 1, 2, 3

图 2中1, 2, 3这3条曲线分别代表了离激励处越来越远的3个接收点的瑞利波信号。从图 2可以看出，在离激励作用位置较远的区域上，各节点的相对幅值的差距变小了，波形在传播的过程中将基本不变，此时产生了瑞利波，此点符合瑞利表面波的性质^[14]。

然后，在模型不同深度的位置取3个点，得到如图 3所示的位移时间曲线。

Figure 3.  Displacement-time curves of nodes 4, 5, 6

图 3中4, 5, 6这3条曲线分别代表了深度越来越大的3个接收点的表面波信号。从图 3可以看出，随着深度的增加，节点振动的位移幅值逐渐衰减，此点符合瑞利表面波的衰减效应。

由此看出，简化后的激光冲击作用的确可以产生瑞利表面波。

在表面上距裂纹左边界10mm处取信号的接收点，同时裂纹深度分别为0μm，230μm，360μm时，接收点处瑞利波信号的位移时间曲线如图 4所示。

Figure 4.  Rayleigh wave signals with different surface crack depths

从图 4得知，当材料表面没有裂纹时，接收点处信号并没有出现瑞利反射回波。当材料表面存在裂纹时，接收点处信号出现了较强的瑞利反射回波。并且随着裂纹深度的增加，瑞利波信号的峰值呈显著增加的状态。因此这可以成为实验判断表面裂纹是否存在的主要依据。

图 5中给出了瑞利反射波的幅值与裂纹深度之间的关系。可以看出，裂纹深度在0μm~260μm时，瑞利波的峰值随着深度的增加而显著增加；裂纹深度在260μm~420μm时，瑞利波的峰值随着深度的增加而缓慢增加；裂纹深度在420μm以上时，瑞利波的峰值随着深度的增加而基本保持不变。这一变化趋势^[15]的原因是瑞利波的能量分布主要集中在表面以下、一个波长的深度范围内，随着裂纹深度的增加，瑞利波被裂纹阻挡，发生反射的能量成分增长，回波的幅值也随深度的增加而增加。当裂纹深度超过一个波长量级时，深度的增加将仅会对反射的能量成分产生微小的作用，此时回波的幅值会随深度增加而产生微小的浮动。这一趋势进一步验证了有限元模型的正确性，但同时也表明了不能简单地根据瑞利波的峰值来判断裂纹深度。

Figure 5.  Peak values variation of Rayleigh wave signals with different surface crack depth

为了使瑞利波的位移时间曲线能够反映与裂纹深度的关系，取反射波的正向峰值与下一次最大正向位移所对应的时间(t₁与t₂如图 6所示)的比值Δ=t₂/t₁为指标。图 7显示了该指标与裂纹深度之间的关系。

Figure 6.  Displacement-time curve of node

Figure 7.  Time ratio Δ versus surface crack depth

图 7中的黑点为通过上述模型计算出来的比值，直线为拟合直线。从图 7可以清晰地看到，时间比值Δ与裂纹深度之间存在着近似线性的关系，并且线性相关系数为98.9%。通过这种关系，可以通过实验测量进而得到表面裂纹的深度。

本文中建立了激光超声无损检测表面裂纹深度的有限元数值模型，考虑了激光工作参量将激光冲击作用等效成脉冲载荷作用，得到了性能良好的瑞利波。与传统的有限元模拟超声波的方法相比，引入无限单元能缩小模型尺寸，提高计算效率，并且也能有效消除边界反射波对计算的影响。模拟结果表明, 材料存在表面裂纹时, 接收点处会有明显的反射回波，借此可以判断材料表面是否存在裂纹。并且瑞利反射波的正向峰值与下一次最大正向位移所对应的时间比值Δ和裂纹深度之间存在线性关系。利用这一关系，可以为激光超声无损检测的实验测量提供一种有效方法。