1. 项目背景与SH0模态机理1.1 为什么在钢板检测里盯上SH0模态钢材在桥梁、管道、储罐、船体里用量极大而这类结构最常见的失效模式就是腐蚀减薄和裂纹扩展。传统超声检测靠纵波直探头、横波斜探头一点一点扫效率低不说大面积钢板检测时漏检风险相当高。导波检测的好处是单点激励就能覆盖很长的传播路径但通常说的Lamb波兰姆波有两个让人头疼的问题频散和多模态。频散意味着不同频率成分传播速度不一样脉冲信号越走越散定位和定量都受影响多模态则意味着信号里可能同时存在好几个波包稍不注意就认错目标。SH波是水平剪切波质点振动方向平行于板面且垂直于传播方向。SH0模态是SH波族里最低阶的那一支它有几个天然优势零阶模态几乎不频散群速度和相速度在很宽的频厚积范围内保持稳定它的位移场是面内剪切型不容易和不需要的模态耦合再有就是SH0对沿板厚的贯穿型缺陷、表面腐蚀这类损伤的敏感度比较符合工程检测的需求。我在这个项目里做的就是把“三维钢板上激励和接收SH0模态”这件事用建模手段完整复现一遍把激发条件、传播行为、缺陷响应都算清楚为后续做定量检测做铺垫。1.2 SH0的质点振动与传播特性回顾这里先铺垫一个关键背景后面建模参数都会用到。SH波在自由板中的频散方程可以写成[ \sin(qd)0 ]其中(q\omega\sqrt{\rho/G - 1/c^2})(d)是半板厚(\omega)是角频率(G)是剪切模量(c)是相速度。这个方程解出来的最低阶就是SH0它的特点是频率越低、板越薄越接近纯剪切波速[ c_{SH0} \approx \sqrt{G/\rho} ]钢结构里(G)大约是80GPa上下(\rho)取7850kg/m³算出来剪切波速约3200m/s。竖向看钢板的色散曲线在SH0这一支上基本是条水平线这就是为什么SH0在工业导波检测里这么受欢迎——信号稳定、好判读。需要强调的是SH0的质点位移方向是平行于板面的横向剪切不产生板厚方向的位移。这个特性在建模时很要命如果激励方式给的是面外力z向力根本激不出SH0只会激出对称或反对称的Lamb波。真实检测里常用电磁超声换能器EMAT来激发SH波因为EMAT的洛伦兹力和磁致伸缩效应可以产生面内剪切力建模时如果也要仿真就得在力边界条件上做文章而不是简单加一个z向载荷。2. 三维建模方案设计思路2.1 建模对象与整体方案选型我这次建模的目标是一块尺寸为1000mm×600mm×10mm的Q345钢板。板厚10mm这意味在1MHz附近的频厚积是10MHz·mm在这个范围内SH0频散很小适合用窄带激励。我选择的建模工具是COMSOL的固体力学模块做时域显式求解同时也用了一个简单的MATLAB解析模型作对照。选择COMSOL的原因很直接它处理压电、电磁、结构多物理场耦合方便后处理里提取节点位移、做波场动画都很顺手MATLAB解析模型则用来在半分钟内快速估算频散曲线和理论波速两个工具形成互验。三维模型比二维模型麻烦不少但必须做三维。原因有三条第一钢板上换能器有实际尺寸不是理想点源激励区域的面积效应只有三维能模拟第二缺陷形状比如圆孔、刻槽的散射波具有方向性二维模型只能计算单一截面的情况无法反映真实的三维散射场第三SH0的质点运动方向在板面内如果只建一个二维截面难以正确表现位移分量的空间分布。2.2 几何建模与边界条件处理几何建模时我先用了一个“半模型”思路利用板中面的对称性只建上半部分对称面设置成低反射边界或者完全对称边界来缩减计算量。前提是激励和结构几何要关于板厚中面对称。SH0模态的位移分布是沿板厚均匀的如果激励力是中面对称的面内力半模型加对称边界可以精确复现全板响应计算量直接砍掉一半。边界条件的处理是这次项目里最需要小心的地方。钢板是有限尺寸的波传到板的边缘会发生反射如果不处理反射波会和缺陷回波混在一起后处理时根本无法分辨。我的处理方式是在板的四周设置完美匹配层也就是PML。PML的厚度取激励中心频率对应波长的1.52倍材料属性逐渐变为复值从而吸收波能。要注意PML内部的网格不能太粗否则吸收效果急剧恶化。2.3 材料参数与阻尼设置Q345钢的弹性参数我取的是杨氏模量E206GPa泊松比ν0.3密度ρ7850kg/m³。开了时间显式求解之后还有一个很重要的问题是材料阻尼。实际检测中SH0在钢板里传播几米之后振幅会衰减但数值模型里如果不加阻尼波会一直保持良好的形状导致仿真结果过于理想。我加了一个各向同性损耗因子用瑞利阻尼的等效方式实现把阻尼比控制在0.1%左右。这个数值不能给太大否则回波幅度衰减过猛缺陷散射信号会被淹没。注意瑞利阻尼的系数不能直接套用结构动力学里常用的5%阻尼比那对导波检测完全不适用。导波问题里通常只加很小的材料衰减实际的衰减系数最好参照你手头实验数据来定。3. 从网格到求解核心实操要点3.1 网格尺寸怎么选才稳显式有限元计算中网格尺寸必须足以分辨最短波长。SH0在1MHz时波长约3.2mm为保险起见我按最短波长伽马20个单元来加密即单元尺寸0.5mm左右。如果只按常规的“每波长510个单元”去做虽然算得动但波前会出现明显震荡后续做振幅定量分析时误差偏大。网格类型选择上我用的是六面体单元。自由网格剖分虽然省事但六面体在波传播方向上形状规则数值频散小。板的厚度方向我分了6层单元这个层数在10mm板上可以保证高频成分在厚度方向也能被分辨不至于出现伪模态。网格尺寸和时间步长是有联动关系的。显式求解为了保证稳定时间步必须满足Courant条件即[ \Delta t \le \frac{h_{min}}{c_{max}} ]其中(h_{min})是最小网格边长(c_{max})是模型中最大波速。钢结构中纵波波速约5900m/sh取0.5mm时临界时间步大概在8.5e-8秒。我实际取的是一个安全系数大约0.8倍也就是6.8e-8秒左右。这个数值如果设置太大求解直接发散设置太小计算时间成倍增加。3.2 激励信号的构造与施加方式激励信号我选择汉宁窗调制的5周期正弦脉冲。中心频率1MHz信号时长约5μs。选5个周期是权衡后的结果周期太少频带太宽会激发出不需要的模态周期太多时间窗太长波包变宽在短板上缺陷回波和边界回波容易重叠。关键在力的方向。SH0的质点位移方向是面内横向激励力必须也在这个方向上。我是在板的表面施加一个沿着y轴方向的面内集中力传播方向是x方向这样激发的刚好是沿x传播、在y方向振动的SH0波。这个看起来简单的细节很多新手在这里栽跟头——用了z向力或者面外力结果屏幕上跑出来的是A0或S0还在纳闷为什么看不到SH0。3.3 后处理如何证明你激励出来的是SH0后处理中判断模态的金标准是看位移分量。SH0在钢板上传播时主导位移分量是y方向位移u_y如果激励位置附近同时出现了明显的u_z面外位移那就说明Lamb波也被激发出来了。我在模型中设置了三个探针分别记录u_x、u_y、u_z随时间的变化。结果显示u_y的振幅比u_z大两个数量级以上这说明SH0占绝对主导整个激励方案是成功的。另外我还在沿传播方向的路径上提取了整条线的位移场并做成动画可以看到一个沿x方向平滑移动的波包质点做y方向的横向振动没有出现厚度方向的弯曲运动。这比单看某一点的时域信号直观得多也更能说服自己模型建的没问题。3.4 缺陷模型的典型处理方法作为检测建模缺陷是不可少的。我在板的中心位置设置了一个贯穿厚度的圆形通孔直径8mm。通孔建模很简单几何上用差集操作把圆柱体减掉就行。但要注意网格上的处理——孔周围的应力集中区必须局部加密否则散射场的细节会被粗网格抹掉。我做了两层局部细化孔边网格尺寸从0.5mm加密到0.25mm保证了孔周围位移梯度的分辨率。更复杂的缺陷比如表面刻槽需要在钢板上表面做一个1mm宽、3mm深的矩形槽网格也需要在槽底和槽壁做局部细化。这类缺陷对SH0的散射特性往往比圆孔更敏感因为槽的方向性和深度直接影响反射系数的大小。我建议在建模初期就把缺陷几何参数做成变量方便后续做参数扫描。4. 关键参数计算与典型结果4.1 频率与波长的换算逻辑建模中所有几何尺寸、网格密度、PML厚度都要围绕波长来定。已知道钢板的SH0波速约3200m/s中心频率1MHz对应的波长为3.2mm。板厚10mm频厚积10MHz·mm。若查阅频散曲线SH0在该频厚积下基本无频散群速度和相速度几乎一致。这个结论很重要因为如果频率提高到2MHz频厚积变成20MHz·mm虽然SH0仍然无频散但波长为1.6mm网格尺寸就得缩到0.25mm模型自由度成倍上升。实际操作中我还要预估激励频带。5周期汉宁窗信号的频带宽度大约是中心频率的20%左右也就是0.91.1MHz。在这个范围内SH0的速度变化不超过0.1%所以波形传播后不会明显畸变。这是SH0最可爱的地方——用来做定位时间-距离换算极其稳定。4.2 模拟结果中读取波速与衰减做完时域仿真后我在传播方向上距离激励点200mm和400mm处各取一个探针得到两段时域信号。通过互相关函数计算两个信号之间的时间延迟我得到波速约3185m/s和理论值3200m/s偏差不到0.5%这从侧面验证了模型和网格设置的合理性。如果要做衰减分析就取同一传播路径上多个探针的峰值振幅做指数拟合。注意在三维板中SH0虽然在板平面上是柱面波扩散但由于板的几何约束振幅随距离的衰减并不像体波那样是1/r关系而更接近(1/\sqrt{r})。如果看到衰减过大先检查是不是PML吸收不够反射回来了再看是不是网格太粗导致数值耗散。4.3 缺陷散射信号的辨识有孔和无孔两种模型分别跑完后我对比了接收点信号的差异。可以看到在激励后的第一波包之后无孔模型只有边界反射回波有孔模型在中间时间段出现了一个明显的额外波包就是孔造成的散射信号。这个散射波包的到达时间对应激励点到孔再到接收点的路径和带宽和激励信号基本一致说明散射过程没有引入严重的频散效应。需要特别提醒的是孔散射波的幅值通常比直达波小一个数量级以上如果激励信号或者边界反射处理不好很可能把散射信号淹没掉。如果仿真中接收点的信噪比不够优先改进边界吸收而不是增强激励幅值——过大的激励幅值可能让模型进入非线性区结果反而失真。5. 常见问题与排查技巧实录5.1 为什么激不出来SH0我在调试中遇到过几次激励方向正确但SH0幅值依然很低的情况。排查后发现几何模型里板的法向方向虽然是z但力加载时使用的坐标系是全局坐标系如果几何体做过旋转局部坐标和全局坐标就会不一致导致看似给了y方向力实际分解到板面上后既有面内分量又有法向分量。解决办法是加载时务必先确认力分量的方向相对几何表面是平行的。另一个常见问题是加载区域过大。如果面内力的加载区域尺寸和波长相当激发的波束方向性会变强侧向泄漏增加沿主传播方向的SH0能量反而减少。我会建议在看模态纯度时先用一个点力做初步验证等确认波型正确后再改成模拟实际换能器尺寸的面力。5.2 边界反射和缺陷回波分不清这是导波建模里最烦的问题。如果板太短或者PML效果不好边界反射会混进缺陷回波的时间窗。我的排查顺序是第一步跑一个无缺陷模型把所有边界反射波峰的位置和振幅标出来第二步跑含缺陷模型把新增的波包挑出来第三步如果新增波包刚好和某个边界反射波混在一起就调整PML厚度或者把接收点位置稍微移开让两者在时间上错开。PML效果不好的原因一般是厚度不够或者PML内网格太粗。有一次我把PML厚度从1.5倍波长增加到3倍波长背景噪声立刻降了10dB效果立竿见影。如果计算资源有限也可以考虑在板边设置“低反射边界条件”尽管不完全等于PML但对散斑噪声的压制已经够用。5.3 计算内存和耗时失控显式计算最怕模型规模失控。1000mm×600mm×10mm的板如果全部用0.5mm六面体网格单元数超过2400万一般工作站根本吃不消。我的做法是先用半模型减一半自由度再把激励和接收区域附近的网格细化、远离关注区域的网格适当放大。但要注意放大网格时必须保证最大网格尺寸仍然小于最短波长的1/10否则波在粗网格区域会产生明显反射。此外时间步长不需要全模型统一。COMSOL的显式求解器可以设置局部时间步但大多数情况下还是靠调整最小网格尺寸来加速。如果只是为了看远场传播和缺陷回波可以先跑一个30μs的短时程等验证算法没问题再延长到时程。5.4 结果中出现高频振荡伪影有时候探针信号尾部会出现高频毛刺这通常是网格尺寸不均匀造成的数值反射。一个典型的错误是把板的大面积区域网格设为0.8mm而孔边加密到0.25mm两者之间没有过渡。波从细网格进入粗网格时会因为数值阻抗不匹配产生伪反射。解决方法是设置合理的网格过渡比相邻网格尺寸变化不超过1.5倍或者使用边界层的渐变过渡。6. 一些个人体会和扩展建议把这个建模做完之后我的最大感受是SH0模态的建模难点不在波动方程本身而在工程细节的把控。激励方向、边界吸收、网格一致性、时间步长的匹配每一项单独拿出来都不算高深但串在一起就决定了模型能不能真实反映物理过程。我建议刚接触导波建模的朋友按照“先解析公式算波速再二维模型验波型再三维模型做定量”的顺序推进这样可以避免一开始就被三维模型的计算量和后处理复杂度拖住。后续扩展方面我打算往两个方向走。一是把电磁超声换能器的电磁场和结构场耦合起来做全物理场仿真这样就能直接模拟EMAT线圈的激励过程而不是简单给一个力边界条件。二是把缺陷参数扫描和机器学习结合起来建立SH0散射信号与缺陷尺寸、方向之间的映射关系这部分做好了会对实际检测中的缺陷定量评价很有帮助。如果你也在做类似的工作欢迎从这些方向继续往下挖模型里的坑踩过了后面就能少走很多弯路。