伍治坚提到COMSOL超声清洗仿真，很多人的第一反应是"压电片加上水槽、求解出水里声压不就行了"。真跑起来就发现不是这么回事，尤其当模型里出现多个28kHz压电片、还要讲究"按一定方式分布"贴在2mm钢质水槽外侧的时候，问题会一下子跳到多物理场耦合、模态匹配、阵列相位关系这些更深的层面。这篇东西我打算从一个实际工程项目拆起：为什么要做多压电片分布仿真、单片激励和阵列激励差在哪、多物理场怎么搭、边界条件怎么给、扫频时看什么，以及最后怎么把仿真结果拿回到工程上判断清洗效果。
开篇能看出来，这个项目想解决的已经是"槽内声场均匀性"问题，不是简单演示"能不能激发出声波"。这个阶段跳过去之后，真正花时间的是压电片怎么排、激励相位怎么设置、声场为什么总是不均匀。我尽量把每一步推导和踩过的坑都说清楚。
1. 从"超声清洗仿真"到"多压电片阵列建模"，问题发生了什么变化
1.1 单片压电片仿真为什么让你越算越心虚
先做一个实验。建一个直径20mm、厚2mm的PZT压电片，贴在2mm厚不锈钢水槽底板外侧，槽内是水，给压电片加28kHz正弦电压，然后扫频看结果。你会发现槽内声压分布极不均匀，离压电片近的地方声压高得吓人，槽的另一端几乎没什么响应。把压电片挪到不同位置，声压云图立刻变一个样。
不是模型错了，是物理上本来如此。一个压电片相当于一个点激励源，激发的板振动在钢槽里传播时会形成复杂弯曲波场，再通过板-水界面辐射声波。因为2mm钢板很薄，弯曲波在板内传播时损耗小、反射强，驻波效应特别明显。单片换能器在工程上做小功率清洗、局部清洗可以，想覆盖整个清洗槽，很难均匀。
所以一旦项目标题里出现"多个28kHz压电片"，本质上已经从元器件仿真切换成了阵列设计问题。你要回答的不再是"压电片能不能振动"，而是"振出来的声场能不能在全槽范围内达到足够且相对均匀的声压幅值"。
1.2 工程上到底想从仿真里拿到什么结果
很多做清洗设备的朋友第一次提需求时会说："你帮我仿真一下，看看超声能不能洗干净。"这个目标没法直接仿真，因为清洗效果跟空化泡动力学、液体温度、含气量、表面张力都有关系，有限元模型不可能把所有气泡都建出来。但工程上需要知道的三件事，仿真完全可以给：
- 哪些区域声压幅值低于空化阈值，意味着清洗盲区；
- 压电片贴在不同位置时，声压空间分布是否均匀；
- 驱动频率偏移（比如从26kHz到30kHz扫频）时，声场是否稳定，还是只在某个狭窄频点冒尖。
说到底，超声清洗仿真的评价指标就是声场分布和空化阈值的关系。纯水在常压下的空化阈值大约在0.1MPa到0.3MPa之间，实际超声清洗槽由于水中溶解气体、微小杂质的存在，空化起始的负压幅值往往比理论低不少。所以在后处理里把你的声压结果换算成声压级，再叠加空化阈值参考线，就比只给一张漂亮的彩色云图有用得多。
1.3 标题里没写完的半句话，决定整个建模边界
这个项目标题的结尾是"向内……"，被截断了。按照常见设计思路补全，应该是"向内辐射超声，实现槽内液体清洗"。这半句话意味着模型至少要包含三段物理区域：压电片和钢槽组成的结构域、槽内水组成的声学域、以及水面上方的空气或自由边界。
空气通常不建模，直接在水面位置给一个声学边界条件。但"给什么条件"是有讲究的。水面如果是开放槽体，近似看成声压为零的软边界；如果清洗槽上方还有封闭盖板、或者液体中插入了工件，边界条件就要改。最稳妥的做法是，水面以上留一段薄薄的空气层并设置完美匹配层，代价是网格和自由度上去了，但能避免因为边界反射造成声场假驻波。
还有更隐蔽的决定因素：整个模型的范围。你有没有把清洗槽的四壁全部建模？槽壁是悬挂还是支撑在台面上？进水口、溢流口这些结构要不要忽略？这些在低频（28kHz不算高）结构振动问题里影响很大，板一薄越容易被边界条件带偏。
需要模型API调用？ 免费领10W Token，多模型网关一键接入 Claude、DeepSeek 等主流模型。
2. 压电-结构-声学耦合建模：不能省的中间环节
2.1 三个域分别用什么物理场，为什么结构场往往被低估
在COMSOL里做这个模型的完整链路是：
- 压电片域：用压电材料本构关系，接入固体力学和静电，启用压电效应多物理场耦合；
- 钢槽结构域：普通固体力学，线弹性材料；
- 水槽内液体域：压力声学（频域）；
- 多物理场耦合节点：压电效应负责压电片上的力-电耦合；声-结构边界负责水与钢槽内壁的力-声耦合。
结构场一定是关键角色，因为这个模型里振动从压电片传到水里，完全要靠钢槽的弹性波传播。如果只把压电片底面的位移当成一个给定的简谐位移边界，直接加在水域上，你相当于假设钢槽是刚体、不产生任何板振动模式，计算结果会丢掉所有槽壁共振细节。真实清洗设备恰恰很依赖槽壁的弯曲振动模式来形成声场，这不是可以随意简化掉的物理过程。
钢槽建模不要偷懒用壳单元。2mm厚相对整个水槽尺寸确实很薄，但在28kHz激励下，沿壁厚方向的应力梯度对弯曲波传播很重要。用实体单元，厚度方向至少剖分两到三层，哪怕为此全局网格变密也值得。
2.2 几何细节：压电片、胶层、水槽怎么逼近真实结构
先说这个项目最关键的尺寸搭配问题。28kHz在水中波长大约53mm，在钢板中的纵波波长大约200多毫米，但2mm薄板在28kHz条件下的主导波动是弯曲波，波长通常在几十毫米量级，与压电片尺寸同阶甚至可能小于几个压电片尺寸。这个量级关系导致的结果是：压电片不是一个点，它在板表面的一小块区域内施加应变，板的局部弯曲响应，而不是简单厚度方向的打桩式激励。
具体建模时给定：
- 不锈钢槽：2mm壁厚，尺寸按实际槽体定，比如200mm宽、300mm长、150mm深，做一个参照模型；
- 压电片：选用PZT-8类大功率压电陶瓷，直径20mm、厚度2mm的薄圆片是一种常见教学设定。但这里必须提醒，真实28kHz功率超声换能器通常不是裸陶瓷片直接贴槽壁，而是通过前后金属块组成夹心式换能器再与槽壁连接。如果项目里确实就是裸片，那你必须在建模时意识到，裸陶瓷片在28kHz附近没有天然强谐振，实际响应非常弱，往往要借助钢槽的复合谐振才能把能量放大。
如果你只是做教学性质的模型验证，用一个厚度方向极化的PZT薄片完全可以，但要在结论里注明这里模拟的是"压电陶瓷片+钢槽体系"的综合响应，不是商用换能器的原样照搬。真要贴近工程，可以在压电片下方再建一块金属匹配块，把整套东西贴到钢槽外壁，模型会复杂很多，但模态特征才真正接近实物。
胶层是一个容易被忽略但影响巨大的中间层。真实贴装时压电片用环氧树脂胶粘结在钢槽外壁，胶层厚度可能只有0.05到0.2mm，但胶的弹性模量比陶瓷和钢低一两个数量级。如果仿真里假设压电片与钢槽完美共节点连接，计算出的振动传递效率会偏高，共振频率也会略有偏移。严谨一些的话在压电片与槽壁之间加一层薄胶域，开始时用0.1mm厚试算。很多仿真与实验对不上，最后排查发现是胶层没建。
2.3 材料参数和损耗：为什么共振峰会变成无限大
结构钢的参数好办。水的参数更简单，声速约1500m/s，密度1000kg/m3。真正麻烦的是压电材料矩阵和系统损耗。
COMSOL材料库自带一些PZT参数，但要注意区分材料型号。大功率超声清洗一般用PZT-8，损耗小、抗去极化能力强；PZT-5H压电常数大但发热严重，不适合连续大功率工况。材料库里如果没有PZT-8，可以自建材料，填入密度、弹性刚度矩阵、压电应力常数矩阵、相对介电常数矩阵。这个矩阵方向容易填错，尤其压电常数矩阵的上下标对应关系，建议先做一块自由压电片的特征频率验证，确认模态频率在合理范围再继续。
损耗项必须给。仿真物理模型默认没有阻尼，结构在谐振频率处的位移响应会趋向无穷大。工程上真实结构总有机电损耗、胶层损耗、结构阻尼，所以扫频曲线上看到的共振峰有明显宽度。给材料加损耗因子的方法很多，最省事的是给压电陶瓷和各结构域指定各向同性损耗因子，取值参考：PZT在0.01到0.03，不锈钢约0.001到0.006。水的声衰减在这个频率很低，可以忽略。
没有阻尼时，你会发现28kHz附近某个结构共振峰的声压幅值高到离谱，这种结果拿到工程评审会上一眼就会被问倒。而且没有阻尼，某些对称模态可能退化成同一个频率，求解时矩阵出现数值病态，甚至导致直接求解器报错。
3. 多个压电片到底怎么排：从波长、相位到电极极性
3.1 先把28kHz的"折合尺度"算出来，才有布置依据
布置压电片之前，先算几个尺度，能帮你避开后面大部分雷区：
水中声速约1500m/s，28kHz对应波长约54mm。槽内声场想要均匀，不能所有压电片间距都远大于一个水中波长，否则它们各自辐射的声波基本不叠加，声场会碎成几块，各管一片区域的感觉。
钢中弯曲波的波长不好直接给死数字，因为它与板厚、频率相关。按2mm厚钢板估算，在28kHz下弯曲波相速度大概在1000m/s级别，波长约40mm上下，和水中波长在同一个量级。这个巧合其实对设计有利：布置一个压电片阵列时，主要考虑的间距范围大概落在30mm到70mm之间是一个较合理的扫描区间。
模型中具体的压电片间距不要拍脑袋定，强烈建议做参数化扫描。在COMSOL里把相邻压电片中心距设成一个参数，比如d，扫描范围为0.5倍到1.5倍水中波长，输出水中平均声压级和声压标准差。你会发现当d取到某个区间时，声场均匀性最优；超出后，声场迅速分裂。这是"按一定方式分布"的最好量化依据。
3.2 贴槽底还是贴侧壁：这是两个完全不同的物理场景
压电片贴槽底外壁和贴侧壁外壁，在仿真里差别极大。槽底通常面积大、平整、下方无遮挡，是优先选择；但槽底本身像一个平板，圆周边界被四个侧壁约束住，它的弯曲模态呈现明显的中心-周边几何形态。如果你在槽底正中央只贴一个压电片，激发出的通常是(1,1)弯曲模态，那么对应在槽内水中的声压分布会有中心高、四周低的强梯度。反过来，贴槽底四周边缘区域，可以激发更多不对称模态，声场整体均匀性反而好。
侧壁的情况完全不同。侧壁板更窄更高，边界条件一边连底，一边可能与顶框相连，压电片贴侧壁外时声波要先经过壁板的弯曲传输再辐射给水，路径距离更长。对于清洗槽尺寸比较扁的情况，贴侧壁的压电片能量主要覆盖邻近区域，对远端的传播效率低。
工程实践里较好的做法是"槽底为主、侧壁为辅"的混合分布。但混合分布有一个建模难点：不同面上的压电片如果驱动相位完全相同，它们辐射的声波会因为在空间中传播路径不同而产生干涉增强或减弱，干涉条纹在槽内乱七八糟。这时候想要均匀声场，单靠同相恒定频率几乎不可能做到，除非整个系统有特殊设计。
3.3 相位关系：同相、反相和交替反相，分别是什么效果
多只压电片的电学连接方式直接决定相位关系。仿真中给每片压电片加电压时，可以在"终端"边界条件里指定复数电压幅值。比如第1片给20V幅值、相位0度，第2片给20V幅值、相位180度，实际效果就是这个压电片的极化方向相对另一片是反的，或者电极反接。
你可以快速做一个实验对比，三组设置：
- 所有压电片同相驱动；
- 所有压电片交替反相（相邻相差180度）；
- 中央一片相位0度、四角各片相位90度梯度设置。
第一组最常见的特征是会形成对称驻波场，声压极大极小值位置比较固定。第二组会把板面分成振动方向相反的若干区域，但激励频率不变时，节点线会清晰可见，这些节点线附近水中的声压往往很弱。第三组相当于旋转激励，能在一定程度上把能量分散到多个板模态，但结构复杂。
需要强调一个结论：相位设置不能消除驻波，只能改变驻波形态。真实超声清洗设备为了让声场均匀，普遍采用扫频或频率调制，让驻波位置随时间不断移动，时间平均后声场才显得均匀。仿真里如果你只算单频28kHz，声压标准差大概率不好看。更贴近真实工况的方法是算一段频率范围，比如26kHz到30kHz，步长100Hz，然后对结果做频率平均，再看平均后的声压分布。这往往比单纯调相位更接近用户实际体验。
3.4 压电片焊脚极性和极化方向的检查清单
仿真中压电片材料默认厚度方向极化。建模时要注意极化轴和全局坐标的关系。圆形薄片如果厚度沿全局坐标z方向，压电矩阵里的d33分量对应的就是z向。很多刚上手的人把材料矩阵拷进去后发现声压为零，查半天，最后发现极化方向设成了沿径向。
电极面也容易搞错。压电片上下两个表面分别设置"电势"和"接地"边界。比如下表面接地，上表面给交流电压。如果你要给多个压电片配置不同相位，每个压电片的上表面终端电压幅值要写成复数，形式如20[V]、20[V]exp(ipi)，分别对应0度和180度。
一个特别影响结果的小细节：极性接反不会导致声压为零，而是导致该片在钢槽中产生的振动位移与预期相反，最终声场干涉图样变化。如果两片物理贴装方向相反，但仿真都按同一极化方向处理，算出来的结果就和实物对不上。装夹方向必须作为几何输入，而不能靠后期调电压相位补偿。
4. 从几何到结果：可复现的COMSOL仿真流程
4.1 建模顺序和分步调试策略
直接一股脑建完整模型很容易翻车。我建议按下面几个阶段递进：
第一步，先不建压电片，先做纯结构特征频率分析。把不锈钢槽的底部壁板或整槽作为自由振动结构，找到25kHz到30kHz范围内的结构模态和振型，记录模态频率、节线位置。这一步的信息量极大。你能看到这个尺寸的槽在这个频段附近有多少可激励模态，以及各模态的振动形态。
第二步，建一片压电片贴到槽底特征振型中应变较大的区域，做频域分析。固定频率范围扫描，看钢槽外壁振动位移和水中声压响应。先不追求绝对声压值，先看共振峰出现在哪里，与第一步算的特征频率是否接近。
第三步，复制该压电片成阵列，加上电压相位关系，做参数化扫描。这个阶段再做布置间距、片数、贴装位置的优化。前面的分步方式，能让问题定位时知道是结构模态算错了、还是声学边界给错了、又或是阵列摆放不当。
4.2 网格怎么画才能保证声场计算可信
网格尺度是整个模型中最多人含糊的地方。压力声学域要求每个波长至少有5到6个二阶单元。28kHz水中波长54mm，按5节点波长算，最大单元尺寸大约10mm，这是一个很粗的量级，普通电脑能轻松处理。但如果你想捕捉钢槽壁的弯曲波传播，结构域网格要细得多，因为钢板弯曲波波长约40mm，可是板厚才2mm，厚度方向至少要剖2到3层，也就是单元尺寸在0.7到1mm量级。这一下子就把表面网格限制到了毫米量级。水声场网格反而被结构网格带动变密。
对压电片区域建议用扫描网格，沿厚度方向指定至少3层。因为压电分析中电场、应力在厚度方向变化快，单元太厚会明显高估压电耦合刚度，算出来的谐振频率偏高。普通四面体网格容易在薄片中出现扭曲单元，影响频率精度。
内存吃紧时可以这样处理：钢槽和压电片用高阶实体单元细画，水域网格最大尺寸按声波波长隔点数放宽，让水域和结构域的网格保持独立，只靠声-结构边界一致网格传递。COMSOL会自动处理交界面的网格映射，但两种网格在边界上的几何尺寸差太大会导致耦合矩阵产生大量小面，计算效率下降。建议边界处结构表面网格尺寸控制在2到5mm左右，避免过细到0.5mm以下。
4.3 研究类型选频域还是特征频率，求解器怎么选
主要研究用频域（频率响应）直接给28kHz正弦激励并扫频。扫频范围的选择，要看你是把28kHz当成中心频率还是标称频率。压电换能器和槽体系统组装后，实际共振频率经常与单片标称频率偏离几千赫兹，所以建议扫频范围设为25kHz到31kHz。频率步长有讲究：共振峰越尖锐，步长越要小。当无阻尼或阻尼很小时，100Hz步长可能跨过峰顶导致漏掉强响应，可以先用500Hz步长粗看峰位，决定是否做自适应细扫。
求解器方面，三维多压电片+结构+声学耦合模型涉及复数代数方程，直接法求解器比迭代法稳定得多，不容易出现声学域收敛残差大的问题。缺点是内存占用高。如果模型自由度超过数百万，可以改选MUMPS直接求解器配合稀疏矩阵存储优化。模型自由度通过网格控制，不用一上来就追求10mm全域细网格，试算阶段放粗到25mm声学网格完全够看趋势。
如果只要找结构共振频率，可以先做特征频率研究，模态数设置为比如20，频率搜索范围15kHz到40kHz。有了特征频率结果，你就能把频域扫描范围缩小到感兴趣的几个峰附近，减小计算成本。
4.4 后处理提取哪些量：声压、声压级、空化区域，一个都不能少
仿真完成后，槽内声压数据要换算出三个层面：
第一层，声压级云图。计算式是20*log10(p/(20e-6))，单位dB。做切片图时，通常在液体中部取一个与槽底平行的水平面切片，用来表示工作区域声场分布。没有声压级而只保留线性声压时，客户看到的是带有负值压力的云图，难以直接理解。
第二层，水中声压的负压峰值。空化泡主要在声波负压相被拉大，随后在正压相猛烈崩溃。压力声学频域解给出的是复声压p，真实瞬态声压可写成p(t)=|p|cos(ωt+φ)。当|p|>1个大气压环境压力时，负压相最低值低于-1atm，液体就容易空化。工程上把0.1MPa到0.3MPa作为粗略参考阈值，看哪里低于阈值。
第三层，槽底结构和压电片的位移。如果压电片输出很大，但槽壁整体位移都很小，说明振动被压电片附近局部区域吸收了，能量没有有效辐射到水中。这种情况要先检查是不是压电片贴装区域的刚度太大，把振动短路掉了，或者胶层过软导致能量耗散。
用全局计算探针可以获得水域平均声压和标准差。平均声压反映总体功率注入强度，标准差反映均匀性。优化压电片分布时，把标准差最小化作为目标，比只盯着某个点上声压最大化要合理得多。
5. 仿真与实槽效果对不上的根因排查
5.1 你扫出了"漂亮的模态"，但水里根本没有声压
这个现象在刚搭建的模型里非常常见。频域扫频结果一看，槽壁在28kHz附近振动位移很大，可是水里的声压却弱得离谱。问题往往出在声-结构边界没有正确启用。检查一下物理场界面里是否真的添加了"声-结构边界"多物理场耦合，并且耦合边界是钢槽内壁与水全部交界面。如果只给了一半交界面，另一半被默认当成声学硬边界，那声压自然只剩一半。
另一个原因是结构振动模式是面内剪切型或扭转型，这种模态法向位移很小，难以向流体辐射声波。所以光看槽壁总位移不一定可靠，需要单独看结构边界法向位移幅值。法向位移大，辐射效率才高；法向几乎为零的模态，哪怕振幅再大也只是做无用功。这也是为什么必须通过声学结果而不是纯结构结果判断超声清洗性能。
5.2 所有压电片同相驱动，声场成了一条"腊肠"
有一个实际案例我印象很深。四条压电片全部贴在槽底，同相驱动，仿真得到的水中声压云图呈现出一条高强声压带，像腊肠一样横在槽内，其余区域声压很差。客户想通过多加几个压电片解决，方向是对的，但同相激励把所有片子的能量集中叠加到了某些驻波腹上，反而强化了不均匀。
解决思路不是单纯加片，而是让阵列的相位、位置错开。你可以用模具胶实贴方案尝试把相邻压电片的激励极性做成反相，或把阵列中心距拉开，让每个片子的近场主瓣不要互相重叠，再看整体声场是否被"抹平"。还有一种手段，把几片压电片的驱动信号频率略微错开，比如26.8kHz、27.2kHz、27.8kHz等错频方案，但这种仿真要用瞬态研究或含多个频率源的谐波叠加，计算量明显增大，需要权衡。
5.3 仿真声压级不低，但实际清洗总是有盲区
听多了会发现一个共性规律：仿真云图上某个局部声压级达到很高，不等于整个槽都能洗干净。实际清洗时，工件本身会遮蔽声波，液体中的气泡、污物也会吸收衰减，而且槽内液体不是纯水，温度、表面活性剂都影响空化阈值。所以把一个静态声压云图当作"清洗能力地图"是危险的。仿真真正可以指导的是趋势对比，用来判断"改A方案是否优于B方案"，而不是预测绝对清洗效果。
实践中我的做法是，把仿真给出的负压峰值和空化阈值对比，划出可能的空化区域范围，再和实际清洗实验的铝箔腐蚀分布对照。如果铝箔腐蚀强的区域确实落在仿真声压高值区，这个模型就可以作为后续优化的标定模型。如果完全不相关，优先检查模型漏了什么边界条件或负载条件。
5.4 频率响应峰的偏移：仿真和阻抗分析仪读数对不上
最后说一个很实用的小技巧。当仿真样机做出来了，用阻抗分析仪实测压电片电阻抗，通常能看到阻抗模值最低或相位翻转的频率点，这就是系统的机电共振频率。把这个实测频率和仿真频响曲线的峰值频率对比，如果差得较多，多半不是网格问题，而是模型里刚度偏大或偏小。刚度偏大表现为仿真共振频率偏高，需要检查胶层、边界约束、压电片附近是否存在过度固定；刚度偏小则反过来，检查不锈钢材料杨氏模量是否错成软钢。
每次调完模型参数，我只改一个变量，重新计算共振峰偏移方向。这比一次改一堆参数再去猜原因高效得多。积累几次标定之后，你会发现自己对这类模型可判断的信心会强很多，也能在迭代中总结出针对新结构的设计直觉。
仿真不是目的，理解声场怎么形成、如何控制和改善，才是做这类项目真正磨练出来的能力。一套COMSOL模型从搭好到可信，靠的是不断和试验现象校准，哪次对不上就往回找遗漏的物理过程。做超声清洗槽的设计，不妨把这篇的方法作为起点，先在模型里把不同分布方案的声场差异算明白，再去做实槽验证，工程效率会高很多。