横波直探头接收信号示意图](placeholder_waveform.png
comsol压电横波直探头 使用pzt4做的横波直探头,在钢中激励1MHz的横波。 自发自收模式,接收信号如图。
做超声检测的兄弟们都清楚,钢构件探伤最怕遇到倾斜缺陷。这时候传统的纵波探头就有点捉襟见肘了,像我们车间最近遇到的T型焊缝检测难题,最后还是靠这个1MHz横波直探头解决的局。今天咱们就扒一扒这个用PZT-4材料做的横波直探头,看看COMSOL里怎么玩转这种特殊结构。
先说核心配置:直径10mm的PZT-4圆片,厚度刚好卡在横波半波长(钢中波长约3.2mm)。这里有个反常识的点——虽然压电陶瓷本身主要产生纵波,但通过45度斜楔块耦合到钢件时,会发生模式转换。我们在材料属性里要特别注意设置PZT-4的压电矩阵参数:
% COMSOL材料参数片段 material = mphgetprop(model,'mat','PZT-4'); material.def.sevatt.piezoelectricity_type = 'matrix'; material.def.sevatt.d = [ 0 0 0 0 12.3 0; 0 0 0 12.3 0 0; -5.4 -5.4 15.8 0 0 0 ]; % pC/N这个压电应变矩阵d的设置直接关系到振动模态。有意思的是,第三行参数控制着面内剪切振动,正是产生横波的关键。实际仿真时发现,如果d33参数设置偏差超过5%,横波能量会直接腰斩。
边界条件设置是另一个坑点。探头底座必须施加固定约束,但斜楔块与钢件的接触面要玩点花样——这里用了阻抗边界条件模拟实际耦合状态。有个小技巧:把钢件的底面设为低反射边界,能有效避免回波干扰,实测接收信号信噪比提升了40%。
comsol压电横波直探头 使用pzt4做的横波直探头,在钢中激励1MHz的横波。 自发自收模式,接收信号如图。
激励信号方面,1MHz的3周期汉宁窗调幅脉冲是经过实测验证的最佳选择。在COMSOL里这样定义时间步长:
# 时间步设置 t_step = 1/(20e6) # 采样率20MHz t_total = 10e-6 # 总时长10μs当探头自发自收时,接收信号会出现典型的双峰特征(如图)。第一个尖峰对应直达波,大约在3.2μs出现,和钢中横波声速3130m/s的理论值完全吻合。第二个宽峰是模式转换产生的表面波,这个特征峰在检测表面裂纹时特别有用。
网格划分必须讲究——在振动区域采用边界层网格,最薄处仅0.05mm。有个经验公式:网格尺寸≤λ/6,对于1MHz横波来说就是0.5mm左右。但实际操作中发现,当网格长宽比超过1:3时,计算结果会出现明显震荡,这时候就得祭出扫频自适应网格了。
最后说个实战技巧:仿真得到的接收信号建议导出到MATLAB做时频分析。下面这段代码能提取信号特征:
[waveform,Fs] = audioread('received_signal.wav'); cwt(waveform,Fs,'amor'); colorbar('off') title('接收信号时频特征');通过小波变换能清晰看到,除了主频1MHz分量,还存在0.8MHz和1.2MHz的边带,这其实反映了探头阻尼特性的仿真精度。
搞定了这些门道,下次遇到奥氏体不锈钢焊缝检测,或者想玩相控阵横波成像,这套模型框架改改参数就能直接上阵。毕竟在超声检测这个行当,能把横波玩溜了,才算真正入了门道。
