做无损检测方向仿真的人大概率绕不开 COMSOL 里的超声裂纹检测模型。这个模型听起来很标准但真正动手搭的时候坑比想象中多得多。我前前后后用 COMSOL 做过不少压电换能器激励、声波在固体中传播的仿真试件裂纹超声检测算是其中最有代表性、也最容易让新手卡壳的一个案例。这篇就把我从几何建模、物理场选择、网格划分到激励信号设置、结果解读的完整思路和踩坑记录整理出来照着走一遍基本能跑通一个二维的超声裂纹检测模型。先说清楚这个模型是干什么的。试件里有一条裂纹我们在试件表面放一个换能器施加一个短时脉冲激励信号让它向试件内部辐射超声波。声波在传播过程中遇到裂纹会发生反射、散射换能器或另一个接收换能器会收到包含裂纹信息的回波信号。通过分析回波的到达时间、幅度、波形特征就能判断裂纹的位置、尺寸和方向。COMSOL 里做这件事本质上是求解压电耦合、固体弹性波传播以及声-结构边界相互作用的多物理场问题。适合什么人参考一是做超声无损检测课题的学生二是做结构健康监测的工程师三是想快速验证换能器布局和激励参数的人。模型不算复杂但涉及的概念密度很高从压电材料参数到声波吸收边界每个环节都得认真处理。1. 模型整体设计与物理场选型逻辑1.1 为什么选二维模型而不是三维裂纹超声检测仿真第一步就是决定做二维还是三维。COMSOL 里三维模型能直接对应实际试件但代价是网格量爆炸式增长。超声波频率通常在 1 MHz 到 10 MHz在钢里的纵波波速约 5900 m/s算下来波长只有 0.6 mm 到 6 mm。要在每个波长内分 8 到 10 个网格单元三维模型哪怕只是 100 mm × 50 mm × 20 mm 的试件网格数也会轻松超过千万。这个规模不是普通电脑能承受的。我自己的习惯是先用二维模型验证物理逻辑把信号激励方式、边界条件、裂纹特征的影响规律摸清楚再决定要不要扩展成三维。二维模型能把计算量控制在可接受范围内通常几十万到两三百万自由度迭代速度快非常适合做参数扫参。如果你的研究对象本身就适合用二维近似——比如试件在厚度方向尺寸远小于宽度方向、裂纹可以简化为贯穿型裂纹——那么二维模型的结论已经具备很强的工程参考价值。1.2 物理场怎么耦合压电、声学与固体力学超声检测的核心是压电效应。你需要压电换能器把电信号转成机械振动逆压电效应机械振动以应力波形式在试件中传播碰到裂纹后部分能量反射回来再由换能器把机械振动转回电信号正压电效应。COMSOL 里这个完整过程需要三个模块压电设备接口或单独的固体力学接口配合压电材料本构描述压电材料的力-电耦合行为。固体力学接口描述试件中的弹性波传播。试件材料按各向同性或各向异性固体处理阻尼参数要合理设置。静电接口或电流接口只在换能器电极区域激活提供激励电压源。如果你只想分析波在试件中的传播特性暂时不考虑逆压电效应也可以用更简化的做法直接在试件某个表面边界上施加一个随时间变化的位移或力载荷模拟探头的压电激励然后只做固体力学仿真。这种方式计算量更小参数调节也更直接。不过要郑重提醒一点如果做完整压电耦合模型材料参数的准确度直接决定仿真结果的可靠性。压电材料比如 PZT-5H的刚度矩阵、压电常数矩阵、介电常数矩阵、密度任何一个参数偏差过大都会导致仿真得到的谐振频率和实际换能器测试频率对不上。COMSOL 内置材料库里有一些压电材料数据但我遇到过不同版本数据不一致的情况保险起见自己核对原始数据手册。1.3 声波传播区域的物理场选择试件里传播的纵波和横波都属于弹性波描述它们的控制方程是固体力学里的动量方程。这里有一个容易混淆的点COMSOL 里有“声学模块”和“固体力学”两个方向但超声检测中波在固体试件里的传播应该用固体力学接口而不是压力声学接口。压力声学频域或瞬态求解的是声压场适用于流体介质不适用于固体中的应力波。只有你说的是空气耦合超声或者水浸超声换能器通过水耦合层向固体试件发射超声才需要在流体区域使用压力声学在固体区域使用固体力学并在流固边界上用“声-结构边界”耦合。2. 几何建模与裂纹处理方案2.1 试件几何怎么建尺寸选择的三条原则试件几何看似简单一个矩形加一条缝但尺寸设计直接影响到波场演化与裂纹信号的可区分性。第一试件的横向尺寸要足够宽。超声换能器激励出的波束在固体中是发散的波前扩展后如果很快抵达左右边界边界反射会与裂纹信号混叠。仿真模型的边界总不会是无限延伸的只能靠尺寸或吸收边界来减小边界效应。我的经验是试件宽度至少是换能器孔径的 5 到 6 倍如果裂纹位置靠近中心左右边界距离裂纹至少 20 个波长以上。第二试件厚度对于二维模型来说是高度要根据检测对象确定。表面裂纹检测时通常有一个临界问题裂纹太浅反射信号淹没在表面波或直达波里裂纹太深可能有多次反射。厚度取 30 到 60 mm 比较合适这样底面回波和裂纹回波在时间轴上能明显拉开。第三换能器位置和裂纹位置的相对关系决定接收信号的形式。如果要模拟单探头脉冲回波模式就把压电片贴在上表面中心裂纹直接位于探头正下方如果要模拟双探头一发一收模式就分别放置两个换能器裂缝偏置在两者连线的中垂线附近。两种模式的回波特征完全不同建几何前先把检测方案定下来。2.2 裂纹的三种建模方法窄矩形、椭圆孔、无厚度弱界面裂纹建模是这道题的灵魂。CAD 软件里画裂纹很简单直接切一条缝但 COMSOL 里不同的几何表示对力学行为的影响差异很大。最直接的方法是建一个窄矩形区域宽度取裂纹真实开口宽度量级比如 0.1 mm 甚至更小然后在窄矩形和基体之间的边界上设置接触条件或直接设置自由边界。这个方法好处是直观、和实际裂纹形状接近但坏处是窄区域会导致极小网格单元严重影响网格均匀性和计算效率。裂纹尖端附近应力奇异网格稍微密一点计算时间就蹭蹭涨。第二种方法是把裂纹简化成椭圆孔或圆孔。这个方法是模拟内部孔洞类缺陷的常用手段。如果你想模拟的不是严格意义上的“裂纹”而是近圆形的收缩孔、气孔类缺陷这个方法更实用。但严格来说圆孔的反射波特征和窄缝是不同的——圆孔以衍射波为主窄缝以镜面反射波为主这两者的幅值和波形分布有显著差异不能混用。第三种方法是无厚度弱界面方法。在 COMSOL 里设置一条内部边界然后利用“固体力学”中的“对”Pair或“接触”功能定义这条界面的法向和切向刚度。给一个很低的刚度值相当于一个很软的分层区域弹性波经过此处就会发生反射和透射。这个方法的好处是避免了裂纹区域的几何网格问题参数上也能调节“透声程度”适合模拟未完全断开的闭合裂纹或者部分接触裂纹。三种方法怎么选我给的参考矩阵裂纹类型推荐建模方法适用场景表面张开裂纹窄矩形缺口几何直接拉伸切除表面检测、开口裂纹内部孔洞/收缩孔椭圆或圆孔内部缺陷检测闭合裂纹/分层内部边界弱界面刚度复合材料分层、早期裂纹任意轨迹裂纹窄矩形裂纹尖端圆角处理形状复杂的裂纹路径2.3 几何清理与工作平面的正确用法很多人在导入外部 CAD 几何时遇到“实体不规则”“边界丢失”等警告做超声仿真时并不需要多复杂的几何直接在 COMSOL 里用工作平面和草图功能建模反而是最省事的。工作平面的作用在这里就体现出来了二维模型直接选 x-y 平面为工作平面画试件矩形区域然后在矩形上表面画换能器压电片的矩形区域和一个表示裂纹的窄矩形区域。用“布尔分割”或者矩形之间的差值运算就能生成基体和压电片两个域。注意压电片和试件之间要保留一个共同边界供后续设置“接触”——这里说的接触不是力学接触而是模型装配后的连续边界条件。从 SolidWorks 导入 STEP 文件再转换时容易出现警告大多数情况是原有几何带的微小圆角、倒角或短边导致网格退化。对于裂纹检测仿真完全可以直接在 COMSOL 里从头建模不用依赖 CAD 文件。如果确实要导入建议导入前在 CAD 里把裂纹切好、圆角和倒角清干净这样警告能少很多。3. 激励信号设计与参数计算3.1 激励信号为什么用汉宁窗调制的正弦 burst超声检测的激励信号不是随便给一个正弦波那么简单的。实际超声探伤仪里探头接收到的是一个短时脉冲时长通常只有几个周期。这种脉冲的好处是时域上足够短时间分辨率高能把不同界面的反射回波在时间轴上分开频域上也有一定的带宽不至于能量过于集中导致无法激发多种波模式。COMSOL 模型里最常用的激励形式是“汉宁窗调制的正弦信号”也叫做 tone burst。表达式长这样f(t) A * sin(2*pi*f0*t) * 0.5 * (1 - cos(2*pi*t / T))其中 f0 是中心频率T 是脉冲总时长。实际设置时可以定义成A * sin(2*pi*f0*t) * sin(pi*t/T)^2因为 sin(pit/T)^2 和 0.5(1-cos(2pi*t/T)) 等价看哪个顺手。为什么必须加窗因为没有窗的纯正弦波一截断频率成分会出现严重的矩形窗旁瓣等于在频域里引入了很多额外的高频分量。这些高频分量在固体里传播速度不一样会让接收信号出现“拖尾”或者“振铃”与真实的窄带超声脉冲差距很大。汉宁窗旁瓣衰减快激励能量主要集中在中心频率附近既保留了短时脉冲的时间分辨率又控制了频谱的纯净度是超声仿真里最常规也最稳的选择。3.2 中心频率和脉冲周期数怎么定中心频率的选择取决于两个因素检测分辨率和衰减。频率越高波长越短对越小裂纹的分辨能力越强。但频率一高材料对超声的衰减也急剧增大传播相同距离后信号幅度下降更严重。钢中常用 2.25 MHz 或 5 MHz 探头铝件里 1 MHz 到 3 MHz 也常用。仿真中我建议先取一个中间频率比如 2 MHz对应的钢中纵波波长大约 2.95 mm。如果裂纹尺寸和波长相当比如 2 mm这种组合下最能看到散射特征。脉冲周期数也就是 T 中包含多少个正弦周期也值得细调。周期数越少时域脉冲越窄时间分辨率高有利于区分相隔很近的回波但频带变宽传播过程中波形畸变更明显。周期数太多虽然信号更“纯”但两个回波脉冲在时域上容易重叠会增加分辨难度。通常取 3 到 5 个周期。我做这类模型默认给 5 个周期这个数量在大多数几何尺寸下有比较好的平衡。实际操作时建议把激励信号定义为 COMSOL 的解析函数然后把时间范围设置成包含整个脉冲持续时间的区间。给两个实用参数A 1 [V] 或者 1 [um]根据自己的激励形式定 f0 2 [MHz] T 5/f0 2.5 [us]如果激励是以电压形式作用在压电片上A 就是电压幅值如果是以位移形式作用在边界上A 就是位移幅值。单位不同后续接收信号的单位也不同看你想怎么分析。3.3 激励信号在 COMSOL 里的具体设置方法COMSOL 里定义一个解析函数很容易。在“全局定义”里新建“解析函数”把上述表达式写进去自变量为 t单位设置好。但有一个细节很多新手会忽略解析函数默认只在定义时间区间外取零如果你没有给函数设置“平滑过渡”也就是没加窗在起始点和结束点会出现突变等于重新引入了额外的高频成分。所以正确做法是直接使用带汉宁窗的完整表达式这样信号本身在 t0 和 tT 处都是零值、零斜率不需要额外平滑。激励函数定义成F(t) A * sin(2*pi*f0*t) * sin(pi*t/T)^2 * (t T)最后加一个逻辑判断“仅当 t ≤ T 时有效”防止数值计算时在 T 之后还给一个非零值。当然你完全可以通过修改研究的时间步来避免但这个判断会让求解器更安全。3.4 激励频率和裂纹尺寸的匹配关系裂纹检测的本质是声波和缺陷的相互作用。当裂纹尺寸远大于波长时尺寸 ≥ 10λ主要发生镜面反射回波幅度大、波形相对简单。当裂纹尺寸和波长同量级时散射特征明显波形中包含丰富的衍射信息。当裂纹尺寸远小于波长时尺寸 ≤ λ/3几乎无法有效反射只能靠幅度变化和混响特征间接推断这部分已经接近检测极限了。仿真时如果要验证某种检测方法的极限分辨能力可以固定激励频率不变扫描裂纹尺寸从 0.5 mm 到 5 mm观察接收信号幅度的变化趋势。这是 COMSOL 做参数扫描的天然优势一次跑完能画出完整的“缺陷响应-尺寸”曲线。如果你算出来的回波幅度和裂纹尺寸之间没有单调关系也别慌散射过程本来就是振荡的不同尺寸存在干涉增强或相消的区域这正是实际检测中识别裂纹难的原因之一。4. 网格划分、吸收边界与时间步长控制4.1 每波长需要多少个网格单元波传播类仿真的网格划分有一条铁律网格尺寸必须足够分辨最短波长。对纵波来说最短波长是 λ_min c_L / f_max其中 f_max 是激励信号的最高有效频率分量。激励信号虽然是窄带但仍有带宽通常取中心频率的 1.2 到 1.5 倍作为有效上限比较稳妥。网格单元尺寸的经验准则是每波长至少 6 到 10 个单元。用 2 MHz 纵波、钢中 c_L5900 m/s 计算波长 2.95 mm网格最大尺寸取 2.95/8 ≈ 0.37 mm。如果同时存在横波c_S≈3200 m/s横波波长更短约 1.6 mm这时网格尺寸要按横波算取 0.2 mm 左右否则横波的数值色散会很严重。在裂纹尖端附近由于应力场剧烈变化建议额外添加局部细化的加密区域最大网格尺寸可以降到 0.1 mm 量级保证裂纹尖端的散射波计算得准。COMSOL 里可以用“大小表达式”或“区域大小”脚本实现局部细化。4.2 时间步长怎么取才不会发散或失真瞬态求解时间步长由两个因素控制一个是数值稳定性一个是时间精度。稳定性上显式时间积分格式如 COMSOL 的广义 alpha 方法用于波传播问题时需要满足 Courant-Friedrichs-Lewy 条件。粗略理解就是每个时间步内波传播的距离不能超过最小网格单元的尺寸。用公式表达就是dt CFL * h_min / c_max其中 c_max 是介质中最大的波速纵波速度h_min 是最小网格尺寸CFL 取 0.2 到 0.5。举例h_min0.1 mmc_max5900 m/s则 dt ≤ 0.5 * 0.0001 / 5900 ≈ 8.5 ns。这个数值比激励信号的周期0.5 μs小两个数量级所以一个 20 μs 的仿真时长需要 2000 步以上。时间精度上即使 CFL 条件满足时间步长太大也会造成高频成分的相速度误差导致接收波形“拖后腿”或者幅度偏差。这就是为什么大多数人做超声仿真时宁可用更小的 dt 也不省时间数据失真的代价远大于计算时间的代价。研究里的“时间步进”参数我一般设置成“0.1/f0/10 50 ns”这样的固定步长然后在求解器设置里把“最大 BDF 阶次”调整到 2 或 3。这里要注意COMSOL 默认的时间步进可能对超声这类高频瞬态问题过度或不足建议先用固定步长跑一遍再对比自动步长的结果确定最终方案。4.3 吸收边界避免边界反射污染信号超声仿真里最讨厌的物理现象就是试件边界反射。如果模型尺寸不够大边界反射会比裂纹回波还强整个信号本底全是乱七八糟的反射波根本没法分析。解决边界的方案有二一是把试件做得足够大让边界反射在观测时间窗内到达不了接收点二是用吸收边界条件。COMSOL 里推荐的是“完美匹配层”PML——在试件区域外围包裹一层专门吸收波能的域。对瞬态弹性波PML 区域需要设置合适的厚度通常和最大波长同量级取 5 到 10 个单元厚度并且在 PML 外边界固定或设置为低反射。但要注意PML 是频域吸收最有效的方法在瞬态计算中可能不完全理想。如果不需要 PML也可以用“低反射边界条件”Low-Reflecting Boundary作为试件边界条件它本质上是一个阻抗匹配边界能压制大部分反射。实际工程中两者可以叠加使用尤其当试件几何较小、边界离检测区域近时。给一个常用的设置组合试件内部区域正常固体力学试件外部底面和侧边低反射边界如果需要彻底消反射再外加 PML 域我实测下来低反射边界在二维模型中已经能把边界反射压到信号幅度的 1% 以下。前提是边界远离检测点让残余反射波晚进入观测窗口否则还是会干扰。4.4 求解器选择瞬态 vs 频域超声波传播的动态过程必须用瞬态求解器因为我们需要看到波包在时域上的分离。频域求解适合稳态周期激励比如分析换能器的模态和谐振频率可以计算频响曲线。但你不能用频域结果直接看裂纹回波的时域到达时间。COMSOL 的瞬态求解器选项里“隐式 BDF”和“广义 alpha”都可以用。对于波传播推荐广义 alpha因为它对高频成分的数值耗散控制更好波形保真度更高。BDF 在阶次较高时会引入数值膨胀对超声信号这种脉冲型信号尤其不利。设置里还可以开启“自适应时间步长”但建议你先固定步长跑通逐步再看是否需要优化。5. 结果分析与裂纹信号特征提取5.1 A 扫信号怎么看直达波、裂纹回波与底面回波模型算完后最重要的结果是接收点或接收换能器电极处的信号称为 A 扫信号。一个典型脉冲回波序列是这样的初始时刻出现一个很强的信号这是激励信号通过压电耦合直接“串扰”到接收端的也叫电磁串扰或直接耦合波。它不代表任何缺陷。随后可能看到沿着试件表面传播的表面波/爬波它的速度比纵波慢会以一个中间到达时间出现。接着是裂纹回波纵波遇到裂纹后反射回来这个信号是我们要关注的。最后是底面回波纵波穿透整个试件厚度到达底面反射回来。在时间轴上读取裂纹回波和底面回波的时间差乘以纵波速度再除以 2因为要走往返就能算出裂纹深度。这就是超声波检测的基本原理。实际读取信号时可以设置一个时间探针Time Probe或直接在后处理里绘制接收点的时域曲线。接收点可以取试件表面中心线上一个点也可以取整个压电片边界上的平均位移。我推荐用“平均位移”作为接收信号因为压电片大面积的平均响应更接近真实换能器的输出电压比单点信号稳定得多。5.2 波场云图看波前怎么传播、如何被裂纹散射后处理里绘制位移场云图是理解波传播过程最直观的方式。选择“表面”绘图变量选“固体位移大小”solid.disp或“总位移”在不同时间帧下查看波前的位置。你会清楚地看到激励初期压电片下方产生一个圆弧形波前向试件内部扩散。波前遇到裂纹时裂纹尖端附近出现局部应力集中同时产生新的散射波前部分背向散射回探头方向部分绕射到裂纹后方。底面反射波会在底面附近形成另一个强波前可能和裂纹散射波在时间或空间上有重叠。绘制多帧动画可以获得最直观的物理图像。COMSOL 的“动画”功能可以把瞬态结果导出为视频做汇报展示时非常有用。5.3 从回波时间反推裂纹深度一个实操计算示例假设模型几何试件高度 H 50 mm纵波速度 c_L 5900 m/s裂纹位于距上表面 d 25 mm 处底面在 50 mm 处。激励信号开始于 t0。纵波从探头上表面到达裂纹并返回的时间t_crack 2 * d / c_L 2 * 0.025 / 5900 ≈ 8.47 μs纵波到达底面并返回的时间t_bottom 2 * H / c_L 2 * 0.05 / 5900 ≈ 16.95 μs这两个回波之间差了约 8.5 μs刚好对应裂纹到底面往返的时间。在 A 扫曲线上找这两个峰值的位置就能很准确地反推出裂纹深度。如果裂纹不是垂直于表面而是有一定倾角反射信号的方向会发生偏移接收换能器收到的回波幅度会明显下降甚至可能收到来自两个不同路径的多个回波这时信号特征就不再是简单的时间差问题而需要更细致的依赖角度分析。5.4 裂纹长度与角度的灵敏度参数扫描怎么做COMSOL 里可以做“参数扫描”把裂纹长度设为一组值比如 1 mm、2 mm、5 mm、10 mm批量求解后把每个工况下的接收信号汇总到同一个绘图里。对比这组曲线能直观看到回波幅度和裂纹尺寸的关系。同样把裂纹倾角作为扫描参数0°、15°、30°可以看到不同角度的反射效率变化。参数扫描时注意每次计算都可能耗时几分钟到十几分钟可以先降低网格精度快速试跑确认几何和物理场完全正确后再跑高精度版本的正式计算。另外COMSOL 会把扫描结果存成一个大的解集合后处理时用每个参数步的表达式withsol或者直接选择参数值提取曲线别搞混了。6. 常见问题与排查技巧实录6.1 网格太粗导致的高频信号衰减一位用户跑完模型后发现接收信号幅度非常小甚至快到几乎看不见裂纹回波。排查下来最大嫌疑是网格尺寸超过了最小波长的网格分辨率要求。如果网格太粗高频成分在数值上被明显耗散等效于一个低通滤波器把信号里最有用的高频信息滤掉了。排查方法在结果里绘制激励信号加载区域的位移时程曲线看它在第一个微秒内的高频振荡是否清晰。如果曲线看起来很圆滑、缺少高频起伏十有八九是网格分辨率不够。解决把裂纹附近和换能器下方区域的网格最大尺寸缩小到 λ/10 以下同时把全局网格尺寸设置成 λ/8。只缩小局部区域网格能避免全局网格量爆炸。6.2 边界反射信号和裂纹回波重叠如果试件尺寸小或者吸收边界没设置好边界反射信号会叠加在裂纹回波上。识别方法是先跑一个没有裂纹的模型把裂纹几何去掉看接收信号里有哪些峰。这些峰全是几何边界、底部和侧边反射引起的记录下它们的时间位置。然后再跑含裂纹的模型对比前后曲线新增的峰才是裂纹贡献。这个方法我每次做超声仿真都会用相当于做了一次“背景去除”。如果边界反射实在消除不掉有两个调整方向一是加大试件宽度让边界反射信号在时间轴上延后避让裂纹回波窗口二是使用前面提到的低反射边界或 PML把边界反射直接压下去。6.3 计算时间过长或内存不足超声瞬态仿真最让人头疼的就是算不动。内存不足常见原因是网格数量超过了内存上限。对策包括优先用二维模型验证方案。使用“自适应网格细化”或者“网格粗化”在远离裂纹和探头的地方降低密度。调整求解器内存策略使用稀疏矩阵直接求解器的“内存优化”选项。把时间范围缩短——不需要算到所有多次反射都结束只保留第一个底面回波到达后的 5 μs 就够了因为后续信号基本都是多次反射和边界效应没什么信息量。6.4 激励信号的“振铃”和数值振荡有些模型会出现激励结束后信号仍在剧烈高频振荡、衰减很慢的现象。这有两种可能一是数值振荡——因为激励信号的高频成分超出了网格分辨能力导致数值不稳定二是真实的压电“振铃”——压电片自身有谐振激励结束后它继续以固有频率振动。区分方法很简单减小激励周期数或调低中心频率如果振荡大幅减小说明是数值问题如果振荡频率和压电片固有频率一致那就是物理振铃属于真实响应不用特殊处理。如果是因为激励函数没有加窗导致高频旁瓣解决方法是确认激励表达式里已经加了汉宁窗且 tT 之后函数返回零。6.5 从导入 CAD 模型时的警告说起热词里有人提到“SolidWorks 另存为 STEP 后导入 COMSOL 有很多警告”这我也遇到过。多数警告来自原始 CAD 模型中的小面、短边、微小间隙和退化几何体。对超声检测模型来说一个干净的解决办法是直接在 COMSOL 里用工作平面重新绘制二维截面不依赖外部 CAD。如果必须用三维 CAD 几何导入后仔细检查“修复几何”功能生成的修复日志把报警的边和面合并清理干净。压电片的电极区域建议用“选择”手动指定别用容差过大的自动选择否则容易把电极选到邻近边界上。6.6 激励信号和接收信号之间时间延迟异常有时候接收信号出现的时间远超理论计算值。这种情况常见于压电耦合模型换能器的压电片自身会消耗一部分时间把电压转换成振动同时如果采用不同物理场耦合在耦合界面上的计算可能引入额外的相位延迟。排查时可以单独做一个“外部位移激励 接收点”的简化模型对比完整压电耦合模型的结果。不过大多数情况下简化模型比压电耦合模型更适合先验证波的传播物理确定几何参数没问题后再叠加压电效应对最终信号的影响。7. 最后小结与实操建议我再强调一个实操中容易忽略但很重要的细节建议把激励信号和接收信号的时间轴对齐记录下来并在后处理里固定显示窗范围。因为 COMSOL 默认的自动范围有时候会把信号的主要细节截掉或者缩放得很难看。固定 y 轴范围比如设定为零附近的一个区间能让裂纹回波和底面回波更清楚地呈现在一张图里。整个模型从几何、物理场到求解器设置每一步看似独立实际上环环相扣。网格尺寸决定你能分辨多高频的信号时间步长决定信号的时域保真度PML 和低反射边界决定你看到的信号是纯净的物理现象还是被边界污染后的假象。如果你想快速验证一个裂纹检测方案的可行性跑二维模型绝对够了而且效率很高。等二维模型的诊断逻辑清楚了再往三维扩展会顺手很多。