做超声仿真这些年，我总结出一个经验：不是模型越复杂越显水平，而是能用最简单的模型把物理讲清楚才是真功夫。压电纵波直探头水耦模型就是我反复向新人推荐的一个入门案例，它把压电换能器、声场传播、耦合介质这三样超声模拟的核心要素全部囊括，计算量又小到普通笔记本就能跑完，非常适合用来开启超声模拟之旅。

这个模型到底解决了什么问题？一句话：让你弄明白超声探头是怎么发出纵波、纵波又是怎么通过水进入被测介质的。在无损检测、医学超声、水声工程的仿真里，这套物理过程几乎无处不在。适合的人群很广——打算用多物理场仿真软件做压电超声的新手、做超声检测课题需要先验证原理的学生、以及想快速搭一个振动-声学耦合Demo的研发人员。我接下来说的这套流程，基于常见实践做了整理，具体参数不一定完美，但足够你第一次就跑通。

1. 模型设计：先想清楚物理场再动手

1.1 压电直探头的工作原理

先说压电效应。压电材料有个独特的性质：给它施加机械应力，内部会出现电荷，这叫正压电效应；反过来，在两个电极之间加一个交变电压，材料会发生形变，这是逆压电效应。超声换能器靠的就是逆压电效应——电压驱动晶片振动，把电能变成机械能，所以本质上它是一个机电转换器件。

直探头名字里的“直”，指的是声波传播方向和探头表面垂直。晶片在厚度方向振动时，辐射出来的声波主要是纵波，也就是质点振动方向和波的传播方向一致。这个纵波的产生机制并不复杂：当交变电压作用在厚度极化的压电晶片上时，晶片沿厚度方向伸缩，把振动传递给前方的介质。探头的结构一般分三层：最上面是背衬（阻尼块），中间是压电晶片，最下面是保护层或耦合层。背衬的作用是吸收晶片向后辐射的声波、抑制杂散振动，让发射脉冲更干净；保护层则用于直接接触被测介质。

在COMSOL Multiphysics里模拟这个过程，就是把这三层结构画出来，给压电晶片加上电压激励，让它厚向振动，从而在水里激发纵波。需要特别强调的是，压电材料的方向和极化方向决定了哪个方向会产生应变、哪个方向会产生电场，这一点在后面设置阶段非常关键，一旦搞反，整个模型就会按照一个不存在的振动模式运行，结果自然完全错误。

1.2 水耦层的价值与替代方案

为什么要在探头和被测介质之间加一层水？核心原因是声阻抗匹配。空气的声阻抗大约只有 430 瑞利，而压电陶瓷的声阻抗高达 3000 万瑞利量级，两者相差五六个数量级，声波一碰到空气就直接被强反射回去，根本进不了介质。水的声阻抗大约是 150 万瑞利，比空气高得多，能作为中间媒介让声波顺利进入被测固体，这就是水耦存在的物理意义。

在水耦模型里，这段水既是声学窗口，也是声传播介质。在COMSOL里我们可以先不关心被测固体，就让声波在水中自由传播、被边界反射，重点看脉冲在时域上的行为。这么做还有一个额外好处：水的声速是已知且稳定的，大约 1480m/s，你可以用这个数值验证波形的到达时间、波长计算，确认仿真的物理量是否对得上。如果以后想模拟接触式探头，可以把水替换成耦合剂或者直接换成固体介质，但作为入门的第一个模型，水是最直观、最不容易出错的。

1.3 为什么用简单模型入门

我带过不少新人，发现一上来就上复杂3D模型的人多半会在调试阶段崩溃。复杂模型的问题在于，一旦结果不对，你根本分不清是几何错了、材料错了，还是耦合设置错了。而简单模型里每个环节都能快速验证：网格密度对不对，看波形平不平滑；边界条件对不对，看反射时间准不准；材料方向对不对，看振幅有没有变化。每一个物理量之间都能靠手算验证，排查问题非常高效。

同时，简单模型的计算量小。2D轴对称模型只求解一个矩形截面，普通笔记本几分钟到十几分钟就能跑完，足够你一天内反复修改几十次参数。这比在大型模型上每次等两小时调试要高效得多。所以，我常说的一句话是：先跑通简单的，复杂的是后来才加进去的。很多研发人员做换能器设计时，也是从这种简化模型起步，先验证振动模态和声场分布趋势，再逐渐加入细节，这是一种非常成熟的工程降阶思路。

2. 从零开始建模：几何、材料和物理场配置

2.1 模块选择与几何搭建

在COMSOL Multiphysics里，做压电超声需要两个物理场接口：一个是压电设备（Piezoelectric Devices），另一个是瞬态压力声学（Pressure Acoustics, Transient）。它们通过“声-结构边界”多物理场节点自动耦合。如果不用压电模块，也可以分别加固体力学和静电，再手动添加压电耦合节点，但那只是曲线救国，没必要。

几何建议用2D轴对称建模。为什么不用3D？因为直探头和水柱在几何上可以看作绕轴向的旋转体，轴对称模型既保留了三维的物理效果，计算量又比3D小一个数量级。在轴对称模型里，画的是一个轴向截面——下方是背衬块，中间是压电晶片，上方是水柱。尺寸可以参考常用探头规格：晶片厚度2mm，直径20mm；背衬厚度10mm左右；水柱高度35mm。水柱高度别太短，不然边界反射会过早进入观察窗口，干扰对直达波的判断。水柱半径可以略大于晶片半径，比如取12到15mm，以减少侧面边界对轴向声场的影响。

2.2 材料参数与边界条件

材料这一步是个重点。压电材料要选择PZT-5A或者PZT-5H这类常用压电陶瓷，在材料库里有现成参数，也可以从文献里找到弹性常数矩阵、压电常数矩阵、介电常数矩阵。对于新手，直接选材料库里的PZT-5A最省事，这些参数是经过验证的，不容易出错。这里特别提一下，COMSOL里压电材料参数可能以“应力-电荷”或“应变-电荷”形式给出，两种形式之间有数学转换关系，选择时注意保持一致性，否则后面可能出现难以置信的应力结果。

水域的参数只有两个：密度1000kg/m³，声速1480m/s。在压力声学里，声速直接影响波长和波数，尽量按照实际使用场景来设置。如果以后想模拟几十兆赫的医疗超声，可以把水换成具体耦合胶的参数，但原理是一样的。背衬材料的选取比较自由，常用钨粉加环氧树脂的复合材料，仿真时可以用一个声阻抗匹配的固体材料近似，背衬域也可以设成吸收性材料，以模拟真实阻尼块的消声效果。

边界条件有几个坑要提前避开。晶片的电极面要设为终端（Terminal）或电压边界，在终端上施加电压激励。背衬的外表面以及水的侧边界、上边界建议使用低反射边界或匹配边界，避免模拟出来的反射波干扰直达波观察。如果你故意要看边界反射，才把它们设成硬声场边界。对于瞬态声压模拟，最稳妥的办法是加完美匹配层，也就是PML，把边界反射压到最低，让波形看起来干干净净。

2.3 压电材料的极化方向设置

这个环节是最容易翻车的地方。压电材料必须设置极化方向，默认坐标系轴方向并不一定等于实际极化方向。如果极化方向选错，模型要么不振动，要么振幅完全不对。在COMSOL里，压电材料属性节点中可以指定极化轴向，典型设置是让极化方向沿z轴，而z轴要和晶片的厚度方向重合。由于我们用的是2D轴对称模型，物理域主要在r-z平面里，要特别注意轴对称坐标下应力-电荷方程的分量映射。

一个稳妥的自检办法是：先只给压电晶片加一个z向电压，不做声学耦合，观察z方向的应变方向。如果应变方向和电压方向的关系符合压电常数矩阵的符号约定，说明设置正确。如果应变方向反了，说明坐标轴设置有问题。我见过有人反复调网格、调求解器，最后才发现就是极化方向反了，白白浪费了一整天的调试时间。这种基础问题在简单模型里暴露得越早，后续做复杂模型时就越省心。

3. 网格划分与求解器配置

3.1 新手必看：网格尺寸怎么定才不出错

压力声学仿真的网格规则很简单：一个波长至少要有6到10个网格单元，太少会把波形的细节磨平，太多则浪费计算量。以1MHz为例，水中纵波的波长 λ = 1480m/s ÷ 1e6 Hz ≈ 1.48mm。如果按照每个波长8个单元来划，最大网格尺寸在0.18mm左右。压电陶瓷里纵波声速大约4000m/s，对应波长约4mm，网格可以放宽到0.5mm。背衬里声速和密度更复杂，但机械振动频率相近，网格取0.5到1mm通常安全。

在划分网格时，水区是波传播的主区域，网格必须满足波长分辨率要求；压电区次之；背衬可以更粗。为了让压电片和水之间的声-结构耦合更稳定，两个域的边界处可以加边界层网格，把边界附近的网格加密到水体网格尺寸的一半以下。网格划分完成后，建议先跑一个小时间范围，比如5微秒，看看波形有没有出现锯齿状的高频振荡，有就说明网格还不够细，需要加密，没有就继续往下走。判断网格收敛性的标准很简单：把网格细化一倍，如果波形关键峰谷的幅值和到达时间变化在可接受范围内，说明当前网格可以信任。

3.2 激励信号与时间步长怎么搭配

激励电压信号常用窄带脉冲，比如一个中心频率1MHz、包含3到5个周期的汉宁窗调制的正弦波。这样一个波包在几个微秒内发射完，频带集中，不容易激发高频数值振荡。如果直接用理想阶跃电压，上升沿太硬，会激励出远超网格分辨率的高频分量，让波前出现非物理的振荡，这是新手最常见的问题之一。

时间步长设置要满足两个条件。一是保证一个波动周期内至少采样20到40个点，1MHz对应的时间步长约为2.5e-8到5e-8秒。二是满足CFL条件，即时间步长乘以声速再除以网格尺寸要小于一个安全值，避免数值误差在介质中扩散。简单起见，直接用1e-8秒做时间步长，对1MHz、网格0.2mm的情况相当安全。总时长建议取20到40微秒，具体看水柱高度。声波往返一次大约需要2倍水柱高度除以声速，35mm的水柱对应约47微秒，所以如果研究反射波，总时长至少要覆盖一次完整往返；如果只看直达波，可以缩短到15到20微秒，节省计算时间。

3.3 求解器设置与计算效率

这种瞬态声-压电耦合模型，推荐使用全耦合求解器，把压电结构振动和水声场放在一起求解，比分开迭代更稳定，也不会出现耦合界面上的“时间差”问题。如果在大型3D模型里觉得太慢，可以把轴对称模型进一步压缩为四分之一模型或半模型，但这种对称简化要确保物理场也具备相应对称性。

求解器的时间步进方式，COMSOL默认会自动选时间步长，但为了波形一致性，我通常手动把时间步长设为固定值，等于上面根据分辨率和CFL条件算出的值。相对容差保持默认的1e-4即可；如果波形出现明显的数值抖动，可以把时间步长再缩小，或者把网格局部加密。记住一个经验：先用大网格加大时间步长快速跑通，再细化网格复核结果的收敛性。这样整个调试循环会快很多，也不容易因为一次等待过久而失去耐心。

4. 结果后处理与波形解读

4.1 如何读取波形数据

跑完之后，不要只盯着最后的云图看，要先把时域波形读出来。在“派生值”或“探针”功能里添加一个点探针，位置选在轴线上距换能器表面10mm处。探针变量用压力p或者质点位移都可以，但最直观的还是压力。

在“一维绘图组”里选择一个点，把p随时间变化画出来。你会看到类似这样的一段波形：先是在t靠近0的地方有一个小幅信号，这往往是压电片表面电极的瞬时影响或数值扰动；紧接着主脉冲以大振幅到达；然后可能出现若干次反射回声。如果水边界设置了低反射条件，那么主脉冲之后还有一些小的残余振荡是正常的，它们来自数值截断和边界不完全吸收，不影响对直达波的判断。我第一次跑通时，看到那个和理论到达时间吻合到小数点后一位的主峰，瞬间觉得之前折腾的几何和网格都值了。

4.2 直达波与界面反射波的区别

在一维绘图里，横轴是时间，纵轴是压力，每条波峰都能和具体传播路径对应。直达波到达时间 t1 = 距离 ÷ 声速。比如探针离探头表面10mm，水中声速1480m/s，那么大约6.8微秒后在探针处出现主信号。如果水柱上方边界是全反射硬界面，那么还会有一列反射波：声波到达界面后反射回来，再经过探针，到达时间大约是探针到边界距离的两倍加上初始距离，再除以声速。用这个时间差验证边界条件的正确性非常直观。

当你在波形图里看到间隔相等、振幅逐渐递减的多个波包时，它们很可能是水柱上边界和侧边界的多次反射。如果这些反射不是你要研究的对象，就应该在边界上使用吸收边界或PML，只保留需要的那部分信号。这一步看似只是调试，但它其实是所有超声类模拟最关键的基本功：学会从时域信号里反推声波的传播路径。我经常在带人时强调一个理念：波形不是一团乱码，它是声场在时间和空间上的投影，每个峰都有一段旅程。

4.3 声场分布云图与频响分析

时域波形能回答“信号长什么样”，要回答“声场在空间里怎么分布”，就要画云图和动画。通过二维绘图组可以查看任意时刻的声压分布，用颜色表示压力正负。你会看到从压电晶片出发的波前随着时间向外扩展，进入水区后逐渐扩散，形成一个典型的由直探头辐射的声场图案。主轴上压力最大，偏离轴线后压力下降，这是因为直探头具有方向性。

把时间步长设得足够小，创建动画，可以看到纵波波前在介质中的完整传播过程。这个过程特别有助于建立直觉，后续分析复杂几何时就有参照系了。如果想做频响分析，可以把时域波形做FFT，观察发射频谱的中心频率和带宽。对于窄带激励，FFT主峰应该落在1MHz附近。如果中心频率偏移很多，要么是压电材料参数不对，要么是晶片几何尺寸对应的实际谐振频率和激励频率不匹配。晶片厚度和谐振频率之间有个近似关系 f ≈ v ÷ (2d)，用这个公式可以快速估算。

5. 常见错误与调试实录

5.1 压电材料参数报错

最典型的报错是“材料参数矩阵不正定”或者“刚度为负”。这通常说明压电材料坐标系和数据表格对不上。解决办法是检查极化方向设置，并确保弹性常数、压电耦合常数和介电常数的单位一致。在COMSOL的材料库中，某些PZT材料的参数是以“应变-电荷”形式给出的，如果你选了“应力-电荷”形式，注意矩阵之间的转换关系，否则求解器会给出明显异常的结果，比如应力达到十几兆帕。

还有一个容易踩的坑：材料参数矩阵里的压电常数符号弄反。压电常数矩阵各项的符号决定了拉伸和电压激励的对应关系，整体反转会导致晶片“不振动”或“振动方向反向”。这时可以通过极化方向的反向设置或者坐标映射校正。简单模型里这种错误很容易发现，因为只要电压激励方向对，晶片就应该在z向有可观察的应变；如果没有，优先检查符号和坐标方向，而不是去调网格。

5.2 波形失真或“炸场”

我见过不少新手的第一个波形是噪声形态——振幅巨大、毫无规律、正负剧烈震荡。这几乎都是网格太粗或时间步长太大导致的数值不稳定性，而不是物理场设置错了。先缩小时间步长看看波形是否改善，再检查网格是否满足每个波长至少6到10个单元的要求。另一个常见原因是电压激励的上升沿太硬，改用汉宁窗或高斯包络后往往立刻好转。

如果你用了PML但仍然出现莫名其妙的边界反射，可能是PML层的层数和厚度不够。PML应该覆盖目标频带内的主要波长，厚度一般取中心频率对应波长的1到2倍，层数建议不低于8层。在调试PML时，可以跑一个只有水介质、没有压电片的纯声学模型，用点声源激励，验证不同入射角下的吸收效果。这个方法虽然绕了一点，但能快速定位反射源。

5.3 内存不足与计算时间过长

3D模型很容易出现内存不足。我的建议很简单：先做2D轴对称，不要一上来就贪3D。轴对称模型已经能覆盖大部分直探头的研究需求。如果确实需要3D，可以把水柱高度缩短，只保留感兴趣的区域，远离轴线的区域网格可以变粗，因为那里声压梯度小，这些做法都能显著减少内存占用。时间总长也别开太长，够观察到目标反射即可，不要漫无目的地跑几百微秒。

计算超时还有一个容易被忽略的原因：求解器设成了自动时间步长，并且载荷是一个很窄的脉冲，导致系统为了解析脉冲而自动把步长压得非常小。这种情况下，手动设定一个合理的时间步长，并关闭“基于局部误差估计”的自适应步进，往往能大幅缩短计算时间，同时不损失波形精度。我在实际项目中常用这个技巧，把一次仿真从半小时压到四五分钟。

6. 模型扩展方向

6.1 加入缺陷反射体：从发射模型变成探伤模型

模型跑通后，往水里放一个刚性小球或者一个平底孔，就变成了一个简单的探伤模型。声波遇到反射体后产生回波，回波到达探针的时间可以通过路程除以声速精确算出来。比如在轴线上、距离探头30mm处放一个直径3mm的小球，那么回波经过的路程大约是 30mm × 2 ÷ 1480m/s，约40.5微秒后到达探针。在探针波形里，你会看到缺陷回波和边界反射同时出现，通过时间窗口可以清晰辨认它们。

这个扩展特别适合当练习题：改变反射体位置、尺寸、材料，观察回波幅值和到达时间的变化，然后把仿真结果和理论公式对照。整个过程只需要在几何里加一个球或圆柱，再设置相应的声学边界条件，比如硬边界代表水中刚性缺陷，阻抗边界代表非刚性缺陷，其余设置完全不需要改动。这是从“会跑模型”到“会用模型”非常自然的一步。

6.2 换能器优化：背衬和匹配层

想提高探头的带宽，可以降低背衬的声阻抗，比如把实心背衬换成钨粉环氧复合背衬，或者在背衬里设置渐变的声阻抗结构。想提高发射效率，可以在晶片和水之间加一层四分之一波长匹配层，厚度约等于该频率下匹配层材料中声波波长的四分之一。例如匹配层材料声速取2000m/s，1MHz时四分之一波长约0.5mm，加这么一层后，晶片发射到水里的声压幅值通常会有明显提升。

这些改动在COMSOL里实现起来很简单：几何里增加一个薄层，给薄层指定材料参数和厚度，重新求解。但物理结果却非常直观，你可以对比有匹配层和没有匹配层的发射波形，看主脉冲幅值、脉宽、中心频率如何变化。这类参数化扫参是最能增强仿真理解的练习，同时也能给实际探头设计提供方向性参考。

6.3 双探头收发模式：还原真实A扫描

如果你想更接近工程实际，可以设置两个探头：一个发射，一个接收。接收探头不需要加激励，只把压电晶片设置为“传感器”模式，读取接收面上的电压。这时你能看到完整的发射脉冲、界面回波、缺陷回波，和实验系统测得的A扫描非常接近。这种收发分离的模式可以消除发射脉冲本身的干扰，更容易识别小缺陷信号。

双探头模型在几何上多了一个压电晶片域和一段水，物理场设置和单探头大同小异，关键是要把接收探头设置成开路或高阻抗状态，让它作为传感元件工作而不是辐射元件。这个扩展做完之后，你对换能器收发特性的理解会深入很多，比如发射探头和接收探头在频带、灵敏度上的差异，以及水中的扩散衰减对接收信号幅值的影响。

我个人的体会是，这个简单水耦模型真正的价值在于“手感”。跑通它，你对压电模块怎么设置、声学与结构的耦合在哪里、网格和时间步长怎样配合，都会有实实在在的把握。之后再加匹配层、加反射体、换复杂耦合结构，都是顺手的事。很多学员问我要一个万能模板，其实没有；但我一定会建议他们先自己搭一次这个模型，因为只有这样，后面遇到再花哨的仿真，你心里始终有一根接地气的基线。

最后再分享一个小技巧：跑成功后，把不同参数下的结果分别保存，然后把波形画在同一张图里对比。比如控制变量调整背衬材料或水柱高度，你会发现很多理论书上的曲线，在仿真图里变得比看十遍文献都更容易理解。仿真这种事，很多时候就是先粗糙、再细致、最后精炼的三步走，急不得。