三维钢板上SH0模态的超声检测建模手记
comsol压电超声SH波检测3维 利用200kHz的压电片PZT-5H，通过黄铜斜楔在4mm厚钢板中激励SH 0单模态。 图1为波场切面的瞬时图。 压电片为自发自收模式，接收了钢板端面回波如图2所示。 图3为群速度频散曲线。 已在钢板上放置了两个点探针对声波测速，差不多就是SH0模态对应的声速3100m/s左右。 虽然在200kHz时A0模态lamb的群速度与SH0接近，但通过做纵波斜探头产生lamb进行比较，可排除是A0。 本模型运行所需内存约16-32G左右，请知悉 模型编号:20＃
用COMSOL搞超声仿真的人都知道，SH波（水平剪切波）在板材检测里是个难伺候的主儿——模态纯净度、斜楔设计、多物理场耦合，随便一个环节都能让人头大。这次的任务是用PZT-5H压电片在4mm厚钢板上激出SH0单模态，顺带验证端面回波和声速，中间踩过的坑倒是值得唠唠。
斜楔角度：Snell定律的暴力美学
想让SH波在钢板里撒欢跑，黄铜斜楔的角度得算准。这里直接上代码暴力求解：
# 计算斜楔入射角（Snell定律） v_brass = 2100 # 黄铜横波速度(m/s) v_steel_sh = 3100 # 目标SH波速度 theta_rad = np.arcsin(v_brass / v_steel_sh) theta_deg = np.degrees(theta_rad) print(f"斜楔理论角度：{theta_deg:.1f}°") # 输出约43°这角度决定了能量传递效率。COMSOL里画斜楔时，直接在几何模块用拉伸+旋转搞定，边界条件记得标“自由变形”，不然仿真完波型能歪到姥姥家。
压电片：自发自收的骚操作
PZT-5H的参数设置是重头戏。材料库里的压电属性往往不全，得手动怼矩阵：
% COMSOL材料属性输入示例 pzt5h.density = 7500; % kg/m³ pzt5h.epsilon = [1700 1700 1450 0 0 0]' * 8.85e-12; % 介电常数 pzt5h.d = [0 0 0 0 741e-12 0; % 压电应变矩阵 0 0 0 741e-12 0 0; -274e-12 -274e-12 593e-12 0 0 0];重点在压电矩阵d的第3行——厚度方向极化，对应SH波的剪切振动。物理场选压电+固体力学，耦合边界选压电片与斜楔的接触面。激励信号用个200kHz的5周期汉宁窗调幅脉冲，峰峰值电压别超过500V（毕竟PZT-5H的矫顽电场也就个位数kV/mm）。
网格：时间域显式的生存法则
显式时间步求解器对网格尺寸极度敏感。钢板厚度方向至少分4层（λ_SH≈15.5mm），但COMSOL的自动四面体网格容易翻车。手动用扫掠网格划分斜楔和钢板：
// COMSOL扫掠网格设置示例 model.mesh("mesh1").create("sweep1", "Sweep"); model.mesh("mesh1").feature("sweep1").set("type", "manual"); model.mesh("mesh1").feature("sweep1").set("nsweep", 6); // 厚度方向6层 model.mesh("mesh1").run();网格质量检查时重点关注斜楔尖角处，雅可比比低于0.3的直接重画。显式求解器的内存需求跟网格数成正比，32G内存的机器跑起来刚好能喘口气。
模态验证：SH0 vs A0的罗生门
虽说200kHz下SH0和A0的群速度接近（图3频散曲线），但区分它们有损招：
- 波形极性：SH波的质点振动方向平行于板面，A0则是椭圆轨迹。用探针监测质点速度的Z分量，SH0的幅值应该扑街。
- 纵波斜探头对比：换个纵波斜楔激Lamb波，时间域信号里A0会带着它的“拖尾回波”出来嘚瑟，SH0的回波则干净得像刚出道的偶像。
端面回波信号（图2）的时间延迟计算也简单：
理论延时 = 2 * 板长 / 3100m/s 实测值误差在±1%内就妥妥的是SH0，A0的声速这时候应该已经跑到2800m/s以下了。后记：模型编号20#的隐藏剧情
这个模型跑完会生成一个wavefieldslice.png（图1），放大看波前会有轻微震荡——这是斜楔边缘衍射的锅。解决办法？把斜楔横向尺寸加宽到3倍波长以上，或者直接加个完美匹配层（PML）吃衍射波。不过内存可能飙到40G，甲方爸爸的服务器要是扛不住...那就假装没看见吧。
comsol压电超声SH波检测3维 利用200kHz的压电片PZT-5H，通过黄铜斜楔在4mm厚钢板中激励SH 0单模态。 图1为波场切面的瞬时图。 压电片为自发自收模式，接收了钢板端面回波如图2所示。 图3为群速度频散曲线。 已在钢板上放置了两个点探针对声波测速，差不多就是SH0模态对应的声速3100m/s左右。 虽然在200kHz时A0模态lamb的群速度与SH0接近，但通过做纵波斜探头产生lamb进行比较，可排除是A0。 本模型运行所需内存约16-32G左右，请知悉 模型编号:20＃
（完）