伍治坚提到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里做这个模型的完整链路是：
- 压电片域：用压电材料本构关系，接入固体力学和静电，启用压电效应多物理场耦合；
- 钢槽结构域：普通固体力学，线弹性材料；
- 水槽内液体域