2025《激光超声表面波检测数值模拟仿真分析案例》1400字.docxVIP

激光超声表面波检测数值模拟仿真分析案例

1.1有限元建模和参数

本文以A60钢轨材料参数为基础，采用有限元方法模拟激光超声波在钢轨表面的传播规律，材料参数如表3.1所示。为了提高运行效率，只截取A60钢轨模型的轨头部分进行模拟，钢轨长100mm。裂纹位于轨距角表面，与列车运行方向呈45°夹角，如图3.2（a）所示。裂纹深度为0.8mm，宽度为0.3mm，长度分别为L1=2.5mm和L2=1.5mm。结合第二章有限元求解理论，对网格进行合理划分，如图3.2（b）所示。网格划分包括网格尺寸、密网格尺寸以及过度网格尺寸，最小网格尺寸为0.08mm，最大网格尺寸为0.5mm。本文只在激光作用区域以及主要声波传播区域采用密集网格，其他区域为宽网格，两者之间则采用过渡网格。仿真时间步长为2ns，总时长为20us，样品材料的初始温度设置为300k。激励源用高斯方程表示，其中光斑直径为0.5mm，距离裂纹10mm处，激光脉冲上升时间为10ns，设置两个接收点，用于接收不同长度裂纹对应的反射波和透射波信号的变化规律，接收点1对应反射波信号，接收点2则对应透射波信号。

此外，对于不同激励位置与裂纹相互作用的规律研究均采用2D截面进行仿真模拟，以节省仿真时间。模型尺寸为22×8mm，裂纹尺寸为0.3×0.8mm，时间步长是2e-9s。设置了三个激励点与接收点，激励与接收点的距离固定为9mm，对应激励点a和接收点a，激励点b和接收点b，激励点c和接收点c。

表3.SEQ表3.\*ARABIC1钢轨特性参数

材料

弹性模量E/Gpa

密度kg/m3

泊松比

热膨胀系数K

A60钢

210

7840

0.29

1.18×10-5

图3.SEQ图3.\*ARABIC2激光超声表面波与轨距角裂纹相互作用规律三维仿真模型

图3.SEQ图3.\*ARABIC3不同激励位置下表面波与裂纹相互作用二维仿真模型

1.2有限元仿真结果分析

激光在热弹效应下，在钢轨表面产生掠面纵波（P）和表面波（R），如图3.3（a）所示，其中P波波速与纵波波速一致，是R波波速的两倍左右，因此传播更快，但能量占比较小，色带较弱。随着时间的推移，表面波传播到裂纹边缘，开始发生散射，部分反射回去，部分沿着裂纹边缘继续向前传播，形成透射波。随着裂纹长度的变化，表面波的传播也发生改变，如图3.4（d）（e）所示，裂纹越长，透射波色带越弱，裂纹越小，透射波色带也越明显，如白色虚线框所示。在图3.2示意的两个接收点位置，获取了L1和L2长度裂纹的反射波以及透射波时域信号，如图3.5所示。与声场图结果一致，裂纹越长，反射波信号越强，透射波信号越弱。

图3.SEQ图3.\*ARABIC4不同时刻，激光超声表面波的传播声场图

（a）接收点1（b）接收点2

图3.SEQ图3.\*ARABIC5不同长度裂纹下反射波与透射波时域图

根据二维仿真模型，进一步分析了不同激励位置下，表面波信号的变化规律。如图3.6所示。当激励点与接收点位于裂纹同一侧时，部分R波与裂纹相互作用发生反射，此时接收到的为反射波(Rr)信号，如图（a）所示；当激励点与接收点分别位于裂纹两侧时，则只能接收到部分透射波（Rt）信号，如图（c）所示；当激励点正好位于裂纹正中心位置时，接收点则能收到部分直达波和部分散射波，如图（b）所示。

图3.7为不同激励位置下对应接收点处的时域信号，a接收点处的R波信号振幅最大，因为包含了直达波以及反射波信号；点c接收到的信号幅值最小，因为大部分的R波信号被裂纹遮挡，只能接收到少量的透射波信号Rt；对于激发点b，激光光斑直径略大于裂纹宽度，小部分光源也被裂纹尖端挡住，大部分光源位于裂纹沟槽处。因此，接收到部分直达波和部分绕射波(Rd)，且绕射波的到达时间滞后于表面直达波。这一结果说明不同激励位置下超声信号存在差异性，为激光超声扫查成像检测提供了理论基础。

（a）激励和接收点在裂纹同侧（b）激励点在裂纹上（c）激励和接收点在裂纹异侧

图3.SEQ图3.\*ARABIC6不同激励位置下，表面波的传播云图

图3.SEQ图3.\*ARABIC7不同激励位置下表面波的时域信号