用COMSOL做超声波清洗仿真或者说在一个钢质水槽底部贴几个28kHz压电片、看水里的声压分布能不能把油污和锈皮震下来——这个题目我前前后后啃了一个多月。今天不写教科书就把怎么建模、怎么调参数、怎么把导纳曲线换算成阻抗曲线、以及中途踩过的那些坑全部分享出来。这篇文章适合三类人准备用COMSOL做声学仿真的入门者、在做超声波清洗设备方案预研的工程师、以及手里有28kHz换能器但声场死活不均匀、想在仿真里查原因的人。看完之后你可以照着里面的步骤把“压电片—钢槽—水”这套模型搭起来至少能拿到谐振频率、声压分布、换能器阻抗曲线这三样硬结果。1. 先拆题28kHz、压电片和钢槽这三样凭什么能凑在一起1.1 28kHz的声学物理和清洗逻辑超声波清洗的核心机制是空化效应。换能器在水中激起高频声压负压半周期内将液体撕裂出微小气泡气泡随后在正压半周期急剧崩溃局部产生高温高压和高速微射流把工件表面的油污、锈层、残留研磨膏“炸”下来。这个机制里频率决定了空化气泡的尺寸、崩溃强度和穿透能力。28kHz属于低频段强力清洗。用Minnaert公式粗略估算水中气泡的谐振半径大概是 R≈3.3/f其中频率 f 的单位是HzR的单位是米。代入28kHzR约0.12mm也就是120微米左右如果换成40kHz、68kHz、120kHz气泡半径分别缩到82、48、27微米。气泡越大崩溃时释放的冲击越猛但单个气泡数量偏少、噪声也大。28kHz正好在“冲击力强”和“空化密度尚可”之间特别适合钢质水槽这种表面有厚重氧化皮、油泥、机加工残留物的清洗对象。声波在水里的传播也要考虑。28kHz对应的波长约等于声速除以频率按水中声速1480m/s算λ≈53mm。低频在这个尺度下衰减慢穿透性好声场能覆盖更大的清洗空间。高频清洗在这个频率下的穿透力不足往往只能处理表面微尘放在重工业场景里不太够用。提醒28kHz只是标称工作频率。实际水槽里声波会在液面、槽壁、工件之间来回反射形成驻波场声压分布有节点和波腹。仿真要解决的第一个问题就是摸清这个驻波场放在哪里高、哪里低。一个清洗系统里28kHz能不能“站得住”还取决于换能器与钢槽底板的耦合方式。钢的声阻抗比水高出约30倍换能器直接贴钢槽振动会先经过钢板再辐射到水里。这个路径在模型里必须完整保留否则频率对不上。1.2 别指望一颗裸压电片天生就谐振在28kHz标题里写的是“压电片”这里我必须泼一盆冷水市面上绝大多数28kHz工业换能器并不是一片光杆PZT陶瓷片。PZT裸片的厚度谐振频率大致是 f_t≈N_t/tN_t 对于常用PZT材料大约1900~2000Hz·m。5mm厚的片子厚度谐振在400kHz上下离28kHz差着十万八千里。你就算把电压加到击穿单片在28kHz附近也只能靠“非谐振强迫振动”硬推振幅非常有限。工业上做28kHz换能器通常用夹心式郎之万结构PZT压电环夹在前后金属块之间整体长度约为半个波长靠金属质量块把谐振频率拉下来。仿真时如果你手里有这种换能器的图纸可以直接按实际尺寸建如果只是想验证“钢槽贴压电片能不能洗”那就把压电片、黏结层、钢槽底板当作一个复合振动系统来算重点看复合系统在28kHz附近有没有结构模态。这一步非常关键。我在初版模型里把压电片当独立器件自由边界下扫频怎么扫都找不到28kHz的峰。后来把钢板和水域一起耦合进来在27.8kHz左右出现了一个明显模态这才意识到我们真正激励的从来不是“压电片本身”而是“压电片—钢槽—水负载”这个整机系统。1.3 仿真该回答哪几个核心问题动工之前先列清楚目标不然很容易陷入“模型建得很漂亮但什么都证明不了”的窘境。我这套仿真主要回答三个问题换能器谐振点对不对通过导纳曲线、阻抗曲线确认工作频率附近有没有有效的电声转换峰峰位离28kHz多远。声场分布均不均匀水里哪个区域声压高、哪个区域是死区波腹区是否覆盖目标清洗位置整缸的均匀性指标能不能压到可接受范围。结构模态会不会拖后腿钢槽底板在28kHz附近的振动模态是什么样压电片有没有贴在节点上底板变形模式是活塞式还是弯曲波主导。这三个问题对应三类结果输出阻抗分析、声压云图、模态振型。后面所有建模步骤都是为这三个输出服务的。2. 建模前先定物理模型几何、材料、声学域2.1 几何尺寸和减化原则几何建模别贪大求全。水槽上的倒角、加强筋、排水口、焊瘤在声学尺度下对几十毫米波长的影响有限建模时一律简化或省略。我用的典型尺寸是钢槽600mm×400mm×300mm壁厚2mm水填充到250mm高度压电片直径50mm、厚度5mm贴在槽底外侧。COMSOL里推荐先用二维轴对称模型做“侦察兵”。因为单个换能器贴在圆形或矩形槽底时二维轴对称可以快速验证材料参数、边界条件、PZT极化方向是否正确几分钟就能算完一次扫频。等二维结果合理了再升级成三维完整模型研究多压电片阵列。三维建模时用Work Plane画草图再Extrude拉伸成实体。压电片、钢槽、水三个域要分别建立最后做布尔并集。水槽内壁和压电片表面不需要特意留间隙在物理场设置里通过“接触”和“声-结构界面”把它们关联起来即可。如果你手头真是郎之万换能器几何会更复杂PZT环内径、前后盖板长度、中心螺栓预紧力都不好直接建模。预紧力在COMSOL里可以用“初始应力”近似前后盖板按普通钢材或铝材处理。初学阶段建议先用裸压电片把流程跑通工程阶段再细化。2.2 压电、钢、水的材料参数细节材料参数是仿真精度的天花板。压电材料必须给全三组矩阵恒应变刚度矩阵 cE6×6、压电应力常数矩阵 e3×6、恒应变介电常数矩阵 εS3×3外加密度 ρ。COMSOL材料库里有PZT-5A和PZT-5H可以直接调用但工程上高功率清洗建议用硬性PZT-4/PZT-8如果你手头有厂商提供的数据手册可以按手册参数自建材料。这里要特别盯紧极化方向。COMSOL默认压电材料沿全局坐标z轴极化压电片厚度方向如果是z方向则电极面就是上下两个圆面。如果你从其他软件导入了材料矩阵一定要确认e矩阵的坐标约定和模型坐标一致。极化方向填反了模型表现是换能器在激励电压下纹丝不动或者导纳曲线完全错误。钢槽用结构钢默认参数即可但建议加损耗因子。超声波清洗钢槽在高振幅下会有材料内耗不加阻尼的话谐振峰Q值会高到几十甚至上百仿真出的声压会虚高。我通常给钢加各向同性损耗因子0.01~0.03给PZT加瑞利阻尼α取两个目标频率点的中间值β用阻尼比0.005~0.02换算。水的参数看着简单也有讲究。压力声学模块里需要密度和声速20℃水取998kg/m³和1480m/s就行。如果你要精确匹配实验温度按温度表插值更稳妥。值得说明的是线性压力声学默认忽略水的粘滞吸收短距离内这种忽略是合理的如果水槽很大、传播距离很长可以补一个衰减项。2.3 声学域的边界条件怎么定边界条件是决定声场形态的重要变量同一个模型边界条件改一个边结果能差出好几倍。自由液面水与空气接触的上表面声学上近似为压力释放边界p0。注意如果你把液面当刚性壁声压会高出一大截那是不符合实际的。钢槽内壁水声场与钢结构的交界面用“声-结构界面”绑定。振动速度连续、法向压力连续这是能量从钢板传递到水的通道。钢槽外壁如果槽体悬空或放在橡胶垫上可以直接用自由边界如果外面还有水箱建议在水域外再加完美匹配层PML模拟无限辐射。压电片电极下表面设为接地上表面设为终端电压源激励幅值按功放实际输出取50~100V。这里有个容易忽略的点水槽底部贴换能器后底部钢板本身也在振动底板的模态会直接决定声源分布。如果你把钢板当成刚性边界只让压电片表面辐射声音那就把“板振动辐射”这条主路径给掐断了声压云图自然失真。钢槽内外表面都要设置为声-结构耦合面即便外侧是空气必要时也可以加外部声学域看辐射损耗有多大。3. 物理场链接压电、结构、声场到底怎么接3.1 Piezoelectric多物理场怎么接出来COMSOL里做压电仿真不是手动去耦合两个物理场而是直接用“压电效应”多物理场节点把“固体力学”和“静电”绑在一起。操作路径是添加物理场时选择“固体力学”、“静电”然后在多物理场节点里打开“压电效应Piezoelectric”COMSOL会自动识别压电材料的e矩阵并把机电耦合项加进去。电极设置上记住两点一是压电片的上表面设置“终端Terminal”类型选电压输入V0100V或你需要的驱动电压二是下表面设置“接地Ground”。对于多个压电片每个压电片都可以独立建一组电极想控制相位就分别给不同终端设置复数电压 V_i V0·exp(-j·φ_i)。在频域求解器里这个复数电压会直接参与谐波响应计算非常方便。很多新手会忽略压电域之外的环境。压电片贴在钢槽底部外侧但钢槽内部是水外部是空气。如果不加外部声学域钢槽背面的辐射就不存在这会影响阻尼和真实谐振频率。我们在模型里可以只建模水槽内部水域外部空气用较低精度的网格域配合PML来处理或者干脆忽略外部空气仅用结构阻尼近似能量耗散。3.2 声-结构界面和液体载荷水和钢槽的耦合在COMSOL里有现成边界条件“声-结构交互Acoustic-Structure Interaction”。在钢槽内壁选中所有与水域接触的面添加这个边界条件COMSOL会自动把声压作为载荷施加到钢结构上同时把结构加速度作为声源项反馈给声学方程。说人话就是钢板振一下水被推一下水声压又反过来压钢板。两个方向同时考虑才是真正的“耦合”。如果只加单向耦合比如只把钢板位移当强制声源忽略水对钢板的附加质量谐振频率会往上偏。水的附加质量效应在低频段很明显我实测在28kHz附近带水负载和不带水负载的谐振峰能差出1~2kHz。水与空气的液面用“压力释放”边界前面说过这里再补充一句真实水槽在工作时液面会波动还会形成小喷泉严格说液面不是完全平的。不过在线性声学假设里把液面当平面压力释放边界已经很准了误差小于工程需要。3.3 多压电片的激励方案同相、反相还是分组相移标题里写的是“多个28kHz压电片”那多片之间的相位关系就成了声场均匀性的调节手段。这一步在COMSOL里实现起来很简单只需要给不同终端设置不同复数电压相位但背后的声场机制值得琢磨。激励方式实现方式COMSOL特点与适用场景全同相所有终端电压相位0底板整体做活塞式振动中心声压高适合单一大型工件的定点强洗行/列反相相邻换能器相位相差180°底板产生弯曲波声场更分散均匀性更好适合整槽多工件清洗分组相移不同区域相位按0°、90°、180°、270°递增可以制造声束偏转或螺旋波用于特定方向增强但对一致性要求高随机相位相位在一定范围内随机分布实验中发现能抑制驻波节点但仿真中难复现工程落地也麻烦从我的项目经验看对于600mm×400mm的槽底先试同相方案看中心声压够不够如果中心区域空化太强、边缘太弱就把中间和两边的压电片设为反相声场马上会被“摊开”一些。注意反相会让底板弯曲模态被更充分激发这时要检查钢板的应力是否超过疲劳极限别光顾着均匀性。4. 网格与求解决定仿真能不能用的两道分水岭4.1 声学网格每波长至少6个单元声学有限元最基础的准则是“网格分辨率必须匹配波长”。28kHz在水中波长约53mm常规做法是每个波长至少6个一阶单元COMSOL压力声学默认使用二阶单元网格可以稍微放宽但我习惯仍然控制在λ/6以内算出来稳妥。具体到模型我给水域设置的最大网格尺寸是9mm在压电片正上方、槽底近场区域加密到5mm远离声源的角落放宽到8~10mm。水域体积约60L0.6×0.4×0.25用自由四面体网格大概能产生几十万到上百万个单元三维计算量不小但这也是换取准确性的代价。钢槽底板2mm厚、PZT片5mm厚这里光看波长还不够必须照顾“弯曲波”的网格需求。弯曲波在薄板里的波长远小于水中的压缩波长所以钢板的弯曲网格必须足够细。我的做法是钢板厚度方向分2~3层六面体或棱柱单元面内尺寸控制在5mm以内PZT片厚度方向分4层以上直径方向至少20个单元这样径向模态和厚度模态都能表达。网格太疏的典型表现阻抗曲线出现大量毛刺、谐振峰偏移几百赫兹、声压云图在换能器附近出现非物理的“振荡斑点”。如果你看到这些现象第一步永远是加密水网而不是怀疑物理场设置。4.2 频域扫描与特征频率分析的配合完整的仿真流程我建议分三步走不要上来就猛扫一个频段第一步特征频率分析。求解范围设在20kHz~35kHz找出钢槽-压电-水复合系统的固有模态。这个步骤能帮你快速定位哪些频率附近存在强耦合模态哪些频率附近干净得没有东西。特征频率分析比频域扫描快得多适合用来“探路”。第二步频域扫描。基于特征频率结果选取28kHz附近±3kHz的区间做扫频步长100Hz发现明显谐振峰后在峰位附近加密到20Hz。为什么要加密因为钢结构加水负载后Q值不算特别高但谐振峰依然可能只有几百赫兹宽步长太大可能正好跳过峰顶。第三步在谐振频率点上做更细致的解。如果目的是出漂亮的声压云图和声场动画那就把频点锁定在谐振峰处做一次高精度求解再用“派生值”提取声压、阻抗参数。求解器方面我推荐直接用直接求解器PARDISO。频域声-结构耦合是复数线性方程组迭代求解器容易在共振频率附近找不到收敛解直接法虽然内存吃得多但稳定性好。内存不够时先把网格尺寸稍微放大或考虑用模型对称性缩减规模。4.3 三维全模型和二维模型怎么选单换能器或同心圆阵列用二维轴对称模型在几分钟内就能完成整个频段扫描。但多压电片矩形阵列不具备轴对称性只能上三维。三维模型的计算量要心里有数。我那个600×400×300的槽水域网格9mm、钢板网格1mm整体自由度大概在三百万量级。PARDISO直接求解一个频率点16GB内存的机器勉强能跑但扫几十个频点就要等很久。两个实用技巧如果阵列具备对称性比如6压电片按2×3排布且中心对称可以用半模型或1/4模型加对称边界大幅减算。频域扫描时先射一个粗扫步长500Hz把峰的大致位置摸清楚再缩小范围做细扫。别在无峰区域浪费计算。4.4 瞬态仿真什么时候用很多朋友问要不要做瞬态仿真。我的建议是研究声场分布和谐振特性频域就够研究“声场建立过程”、“气泡动力学”、“驻波形成”这类问题时再考虑瞬态。瞬态仿真的时间步长得满足库朗条件一个声学周期至少采样20步28kHz对应周期约35.7微秒时间步大概取1~2微秒。为了等声场稳定至少要算15~20个周期也就是700微秒左右。这个计算量在二维模型上可行三维全模型就有点吃力了。瞬态模型还能顺便观察非线性和空化初生但那已经不是线性声学能覆盖的范畴了需要额外引入气泡动力学方程。5. 结果分析从声压云图到阻抗曲线5.1 用声压分布评估清洗效果声压云图是超声波清洗仿真里最直观的输出。在COMSOL后处理里选中水域绘制压力声学“声压幅值p”的等值面或切面云图单位用Pa或kPa。色标建议从0拉到1MPa这样波腹区和死区一目了然。怎么判断清洗效果好业内常用的经验是声压幅值大于0.1~0.3MPa的区域才可能发生有效空化低于这个值就属于“弱清洗区”。纯水的空化阈值理论值很高但实际自来水里有大量微气泡和杂质阈值会显著降低所以工程上常用0.1MPa作为粗筛选限。我把后处理做了两个实用动作一是在目标清洗高度比如工件架设高度50mm、100mm、150mm画水平切面分别统计声压超过0.2MPa的面积比例二是用“体积积分”统计整个水域里声压超过0.2MPa的体积占比。这两个指标比单看云图要客观得多写仿真报告或向领导汇报时也更容易说清楚。匀匀性指标可以算目标面上的“变异系数”CV 标准差 / 平均值。CV小于30%已经算不错的声场如果大于50%说明死区太大需要考虑调整压电片排布或激励相位。5.2 导纳曲线换算成阻抗曲线直接抄这个操作谐振器仿真里最常被问到的就是把导纳曲线换算成阻抗曲线。COMSOL默认输出通常是终端电流I不直接给Y和Z换算其实非常简单但网上资料讲得少这里把操作细节写清楚。第一步在压电片上表面激励电极定义一个积分算子intop_elec用于计算电极面的电位移通量。第二步求出电极总电流。频域下电流 I jωQ其中 Q∫D_n·dSD_n是电位移法向分量。在COMSOL全局表达式里定义Itotal j*2*pi*f*intop_elec(Dz)注意Dz要换成你模型里实际的法向电位移分量如果极化方向沿z电极面法向也是z直接用Dz即可。第三步得到导纳 Y Itotal / V0。导纳是复数但COMSOL后处理里可以提取实部G和虚部BYin Itotal / V0 G real(Yin) B imag(Yin)第四步换算阻抗 Z 1/Y。复数取倒数公式如下R G / (G^2 B^2) X -B / (G^2 B^2) Zabs sqrt(R^2 X^2)然后在结果里绘制“频率—R、X”或“频率—Zabs”曲线就能看到串联谐振峰阻抗最小和并联谐振峰阻抗最大了。这个办法和BAW谐振器仿真的流程完全一样本质都是通过谐波分析从电学端口看器件的频率响应。提示导纳曲线上两个明显的频率点一是串联谐振频率fs导纳最大、阻抗接近纯阻二是并联谐振频率fp导纳最小、阻抗很大。两者的间隔与有效机电耦合系数相关。我做完这套模型后在导纳曲线实部上看到28.1kHz附近一个峰与之前实验测到的换能器串联谐振频率非常接近心里的石头就落地了。5.3 模态分析与清洗均匀性的关联特征频率分析不仅能告诉你谐振频率还能给出振型。看钢槽底板在28kHz附近的模态振型特别要关注底板在压电片安装位置处的位移方向。如果压电片正好贴在模态节点上恭喜你白忙一场——节点处位移为零振动完全激励不出来。我遇到过的一个典型情况4片压电片按矩形排布中间两片恰好落在底板(1,1)模态的节线上结果仿真里中间区域的声压始终上不去。后来把中间两片向波腹位置挪了30mm声场立刻改观。这种问题不仿真根本发现不了反正我是不想靠反复焊接换能器来做试验。多压电片阵列的排布原则仿真里可以总结成一句话让每片压电片尽量落在底板模态的波腹附近同时尽量避开与相邻模态节点重叠的位置。初始设计时按“半波长间距”布阵通常是比较稳的28kHz水中半波长约26mm但钢底板中弯曲波长比这短直接把水中的λ/2套用并不合适。稳妥做法是用特征频率分析先看模态再决定间距。6. 坑位实录最容易让仿真翻车的几个地方6.1 谐振峰找不到或者峰位明显偏移这是我最常被问到的问题。排查顺序如下压电材料参数对不对。特别是e矩阵或d矩阵的数量级和单位。COMSOL库里的PZT-5A、PZT-5H可以直接信自定义材料一定要逐项核对。极化方向是否一致。电极面和极化方向不匹配会造成压电耦合系数为零的假象。钢和水域是不是真的耦合上了。如果声-结构界面没加到整个钢槽内壁水负载带来的频率偏移就不会出现。网格够不够细。网格太疏会让谐振峰“矮胖”甚至消失。阻尼设太高。阻尼过大时峰会被抹平看上去像没有谐振。有一个快速验证技巧先跑一次纯结构模态暂时关闭压电和声学看看钢槽底板在25~30kHz之间有没有想象中丰富的模态。如果没有说明钢槽几何或支撑边界设错了如果有再加上水负载观察频率向下偏移这是正常的附加质量效应。6.2 声压云图出现“满缸乱跳”的非物理振荡声压云图如果看起来很热闹布满了高空间频率的亮斑多半不是物理现象而是数值误差。常见原因网格不连续或局部太粗导致声学方程在过渡区产生寄生振荡。钢槽尖角、压电片边缘等几何奇点处应力奇异网格再细也会出现局部高压。可以在结果处理时忽略这些极个别单元把色标上限压到0.5MPa以下就能看清整体分布。PML参数没设好。如果加了完美匹配层PML厚度至少要大于最高关注频率波长的一半也就是28kHz约26mm以上我才敢放心用。遇到振荡首先做网格无关性验证把水域网格从9mm调到6mm看声压峰值变化是否在5%以内。如果变化很大说明原网格不合格继续加密如果变化不大那振荡大概率来自几何奇点可以靠后处理压制。6.3 终端电流异常导纳曲线形状不像换能器正常的换能器导纳曲线在谐振频率附近会出现一个明显的“圆环”或“峰”偏离这个特征多半是电气设置出了问题。检查电极是否只选在压电片的上表面。如果误选了侧面或底面静电容会变得离谱谐振峰会被淹没。检查压电片极化方向和e矩阵符号。某些材料库的符号约定不同导致正压电和逆压电方向不对电流方向反了导纳虚部符号看着就很怪。终端设置为“电压”时要保证V0不是0。检查时可以先扫1kHz、10kHz、20kHz几个频点看在远离谐振时导纳虚部是否接近“C0*jω”即纯容性。如果不是说明静电容路径有问题。6.4 计算时间太长、内存不够三维全模型扫频确实费时间我自己的做法是分层次解决先用二维轴对称模型把所有物理设置跑通确认频率和振型都对再上三维。三维模型先做特征频率分析把关注频率压缩到一个小范围再做细扫。用1/2或1/4对称模型减负。对称边界在声学里就是“硬边界”结构里就是“对称/反对称条件”注意压电片的极化方向在对称面两侧是否对称不能想当然。如果内存不够PARDISO可以用“内核外存储”速度慢一点但能算实在顶不住就换迭代求解器GMRES加多网格预条件但要先在小模型上验证收敛性。7. 仿真结果怎么和真机实验闭环7.1 用阻抗分析仪实测换能器曲线仿真得出导纳和阻抗曲线后最好用阻抗分析仪常见型号有PV520A、PV70A、Keysight 4294A等实测贴在槽底的换能器阻抗特性。测试时把换能器按实际工况贴在钢槽上、槽里装水不要拿裸换能器测否则和仿真对不上。实测结果与仿真对比时允许存在5%~10%的偏差这是正常的。偏差主要来自三处压电片与钢板之间的黏结层环氧树脂胶的弹性没建进去钢槽支撑条件橡胶垫、支架在模型里被简化了实际水温与仿真设定温度不同导致声速变化。如果偏差超过15%优先检查建模细节尤其是黏结层厚度和边界约束。我后来在模型里加了0.2mm的环氧黏结层实测和仿真的串联谐振频率偏差从原来的700Hz缩小到80Hz左右。这个细节很大程度决定了阻抗曲线能不能“对上”。7.2 用铝箔空化实验做声场验证仿真里说哪里声压高不代表那里空化就一定强最好用最便宜的方法验证铝箔空化法。取一片厚度约10~20微米的铝箔剪成合适大小用支架悬在目标清洗位置开机30秒左右取出观察。铝箔上被空化微射流打出密密麻麻的凹坑、小孔的区域就是空化强度高的区域平整光滑的区域则是弱区。我拿仿真预测的声压云图和铝箔实验结果做过对比铝箔受损最重的几个点基本都落在仿真的高声压区说明线性频域模型的趋势判断是准的。不过铝箔坑洞的“密度”和仿真声压幅值并不是严格线性关系空化强度还会受水中气泡含量、温度、气体溶解度影响所以别指望定量完全一致看趋势就够了。7.3 压电片布局调整的经验数值仿真最终要落到工程排布上给出一组可以当初始值的经验参数对于600mm×400mm×300mm的钢槽水量约60L清洗温度40~60℃建议用6~12个28kHz换能片单片功率50~100W总功率密度控制在10~50W/L。初始排布采用3×2或4×2的规则网格间距按仿真模态结果微调优先把高声压区铺到工件清洗面上。每个换能器的激励电压不必完全相同。仿真里如果发现某个位置落在模态节点附近与其硬加电压不如挪位置或改变相位。匹配电路方面各换能器的实际阻抗会因安装位置不同而略有差异量产时用同一套匹配参数可能不是最优至少要在小批阶段逐一测阻抗。最后分享一个我做这个项目后的体会超声波清洗的COMSOL仿真不能保证你一次就能把样机做成功但它能把“瞎试”变成“瞄准再打”。多花两天建模后面省下来的焊接、拆装、烧功放的折腾可远不止两天。仿真给出的换能器布局直接拿去做了铝箔实验声压高的地方铝箔确实坑洞最密这条从模型到样机再回到模型的验证路径是能闭环的。