1. 从聚焦说起超声换能器到底在聚焦什么这两年做超声相关的仿真最常被问到的问题就是换能器怎么模拟聚焦相控阵的相位延迟怎么设焦点为什么不对COMSOL里到底用哪个模块合适这些问题背后其实是同一个需求——把声波能量精准地送到指定位置不管是用于超声检测、治疗还是无损探伤聚焦都是绕不开的关键环节。先说清楚超声聚焦的两种物理路径这决定了你后续所有仿真方案的走向。第一种是几何聚焦。换能器本身就做成凹面形状比如球面凹形晶片或者平面晶片贴一块曲面声透镜。声波从曲面发出之后由于传播路径长短不同会在某一深度自然汇聚形成焦点。这种聚焦不依赖电路控制是靠物理形状把能量憋到一个点上。它适合焦深固定、焦点不动的情况比如很多单探头超声探头、医用HIFU换能器就是典型。第二种是电子聚焦也就是相控阵聚焦。换能器是由几十甚至上百个小阵元组成的阵列每个阵元独立施加激励。通过给不同阵元设置不同的时间延迟频域里就是不同的初始相位可以让每个阵元发出的波前在某个预定位置同相叠加形成焦点。更妙的是只要调整延迟表焦点就能在声束方向上前后移动、甚至横向偏转。这就是相控阵电控聚焦的核心优势——不需要物理改变探头形状就能动态改变聚焦位置。COMSOL对这两种聚焦都能很好地支持。单阵元聚焦用压力声学频域就能跑得很漂亮相控阵聚焦在二维模型下配合参数化扫描也能做得很实在。这个领域仿真最大的价值在于它可以让你在打样之前就看到声场长什么样、焦点落在哪、旁瓣有多高、能量损失有多少——这些东西靠经验公式估算很容易偏差直接上实验又费钱费时间仿真刚好补上这一环。这篇内容就是把我实际做过的COMSOL超声换能器聚焦和相控阵聚焦仿真的完整思路、参数设计、踩坑经验整理出来。适合正在做超声探头设计、相控阵声场分析、或者刚入手COMSOL声学仿真想少走弯路的同学参考。2. 建模型之前物理场选择与几何简化思路2.1 模块与接口怎么选COMSOL里跟声学相关的接口不少但超声换能器聚焦仿真的主力是压力声学模块。做频域稳态分析用压力声学频域物理场接口做瞬态脉冲传播用压力声学瞬态。如果你要模拟压电晶片本身的机电转换过程那还得加固体力学和静电接口通过压电效应多物理场耦合到一起。但在实际工程阶段我大部分时间不会一上来就做完整压电耦合。为什么因为压电模型会引入非常小的特征尺度——压电晶片的厚度可能只有亚毫米级而声场传播区域长达数十厘米。这种跨尺度模型网格量巨大而且压电材料参数弹性矩阵、压电常数矩阵、介电常数矩阵稍微给得不准谐振频率就对不上焦点算出来也差之千里。更实用的做法是把换能器表面当作声源边界来处理。在压力声学接口中通过法向加速度法向速度或者边界声源特征来设置激励把换能器的辐射问题简化成一个给定振幅和相位的边界条件。这样既能模拟聚焦效果又不需要求解压电域计算量小一个数量级而且聚焦声场的核心物理——波前叠加与衍射——一点都不会丢。2.2 几何应该简化到什么程度几何建模最大的原则是能用二维就不用三维能轴对称就轴对称。单凹面换能器的聚焦天然是旋转对称结构。用二维轴对称模型一个截面就能还原完整的三维声场分布计算量只有三维的几十分之一。COMSOL的二维轴对称坐标中x轴是径向z轴是声传播方向换能器凹面用圆弧曲线构建即可。相控阵线阵如果阵元长度方向垂直于扫描平面的尺寸远大于声波波长那可以近似为二维平面应变问题直接在二维直角坐标里建阵列。只有遇到小尺寸矩形阵元、三维聚焦或者偏心聚焦才值得上三维模型。还有一个小技巧模型域只建换能器到焦点附近的主要传播区域外围用完美匹配层PML吸收模拟无限大介质。很多新手把模型建得很大想靠边界隔开反射结果边界反射照样污染声场。PML的有效性远超加大模型范围而且网格开销更低。2.3 材料参数和声学属性聚焦仿真最常见的传播介质是水和空气。水用密度1000 kg/m³、声速约1500 m/s空气用1.25 kg/m³、343 m/s。这些参数在频域声场计算中影响波长进而决定网格尺寸必须先确定。换能器阵元如果做边界声源近似不需要设置材料如果做压电耦合PZT压电陶瓷材料库里有PZT-5A、PZT-5H等常用型号可以直接选用。介质材料设置各向同性声速即可损耗参数如果只关心焦点形态可以先不给算完再评估衰减影响。3. 单换能器几何聚焦仿真凹面换能器全流程3.1 关键参数与几何构建假设我们要仿真一个凹球面聚焦换能器工作频率1 MHz水中声速1500 m/s则波长λ c/f 1.5 mm。换能器有效孔径口径设为20 mm曲率半径R为40 mm。根据球面聚焦公式焦点位置约在曲率中心附近也就是焦距F≈R40 mm处。孔径与波长的比值决定了聚焦的锐度。这里孔径20 mmλ1.5 mm孔径约13.3个波长属于强聚焦换能器焦点声压增益通常会很高。用COMSOL的草图工具在二维轴对称平面中画一条半径为40 mm、对应口径20 mm的圆弧作为换能器辐射面圆弧的几何中心落在对称轴上这就是凹面晶片形状。传播区域用一个矩形域表示径向从0到15 mm轴向从0到80 mm覆盖焦点并留出后续观察区间在域底部设置PML层厚度可取10 mm约6-7个波长。PML外边界固定为零声压或其他吸收边界均可。3.2 激励设置与频域求解在压力声学接口中将凹面边界设为法向加速度边界条件输入值为常数幅值如 a 1 m/s²。频率设为1 MHz。这里如果不加其它说明默认所有边界点同相同幅振动凹面几何自动产生聚焦效果——这正是几何聚焦的仿真原理波前从曲面各点出发曲面中心到焦点的光程最短边缘光程较长但由于各点同时振动波前在焦点处同相汇聚产生声压极大值。网格划分是这一步的重头戏。规则二阶单元时最大网格尺寸不超过λ/5即0.3 mm如果追求焦点附近的声场精度在焦点区和换能器正前方区域局部加密到λ/8到λ/10。PML区域用映射网格分成多层保证吸收效果。求解器用频域直接求解器默认MUMPS或PARDISO都可以。频率单一不用扫频这一步通常几十秒到几分钟就能算完。3.3 声场查看与焦点评价计算完成后在声压级或绝对声压表达式中查看分布在z轴上的声压曲线。焦点位置的判断标准找到轴向声压最大值对应的z坐标如果声压最大值出现在40 mm附近说明几何聚焦设计正确。最大声压值相对换能器表面激励声压的放大倍数就是该换能器的聚焦增益。再看径向分布在焦点平面zF上画一条横向声压分布曲线找到-6 dB宽度即半功率波束宽度这就是焦斑尺寸。对1 MHz、20 mm口径、40 mm焦距的情况焦点声压半宽大概在1-2 mm量级与实际经验吻合。这里有一个非常值得观察的现象如果孔径太小比如几个波长几何聚焦的焦点会明显比理论位置偏近且焦斑变宽。这就是衍射效应压倒了几何聚焦——孔径边缘产生的衍射波与直达波干涉破坏了理想焦点。仿真能直观看到这个过渡过程对设计很有指导意义。3.4 加个透镜怎么仿真透镜聚焦与凹面聚焦原理类似但几何模型要增加一个透镜层。平面换能器贴在透镜平面一侧透镜另一侧做成凸面或凹面。在COMSOL中建一个透镜材料域设置不同于水的声速和密度比如固体透镜用环氧树脂声速约2600 m/s、密度约1200 kg/m³透镜表面与水的界面使用声-声边界连续条件换能器面同样用法向加速度边界。透镜聚焦仿真最容易出问题的点在透镜材料声速和水声速的对比。声速高于水的透镜是汇聚性质低于水的透镜是发散性质设计错了焦点就变成焦散。在COMSOL里把透镜的曲率半径设为变量做一个参数化扫描可以很直观地看到焦距随曲率变化的曲线——这是优化透镜设计最实用的方法。4. 相控阵聚焦仿真相位控制才是核心4.1 阵元设计与相位计算相控阵聚焦与单凹面聚焦最大的区别在于阵列阵元是平面的波前的形状完全靠相位延迟来塑造。仿真中我需要给每一个阵元分配不同的初始相位使得所有阵元的声波在目标焦点处同相叠加。几何构建以线阵为例设计16个阵元中心频率1 MHz阵元中心间距d 0.75 mm刚好λ/2这是医学超声和相控阵探头的通用间距能有效抑制栅瓣。阵元宽度一般取0.5-0.6 mm阵元之间的间隙用绝热或硬边界模拟实际中就是切缝。相位计算是关键。假设焦点位置在坐标(x_f, z_f)第n个阵元中心坐标为(x_n, 0)那么第n个阵元到焦点的距离是d_n sqrt((x_n - x_f)² z_f²)所有阵元发出的声波要在焦点同相就需要让各阵元的初始相位补偿掉距离差造成的相位差。设参考阵元到焦点的距离为d_ref则相位延迟为φ_n 2π(d_n - d_ref) / λ这个相位值以弧度为单位在COMSOL中可以直接在边界条件的复数表达式中使用。推荐的做法是手动计算好16个相位延迟值写成表格然后在COMSOL中建立全局参数每个阵元对应一个参数。阵元与阵元之间的边界是独立的。在压力声学接口中给每个阵元表面施加法向加速度边界幅值相同但相位不同。比如第n个阵元的边界表达式写作a_n A * exp(-1i * φ_n)注意COMSOL中指数函数是exp虚数单位是1i或j需要统一设置。幅值A的单位是m/s²。4.2 焦点固定阵列的中心与偏转聚焦仿真焦点在阵列正前方的情况最简单所有相位关于阵列中心对称。例如焦点在z_f 30 mm、x_f 0处那么阵列中心两侧的阵元相位呈对称分布中心阵元相位最小边缘阵元相位最大。计算完相位后求解频域模型查看声场分布。你会看到在没有相位控制时所有阵元同相声束会保持平面波方向传播没有明显的焦点施加相位延迟后声场在(0, 30 mm)处形成一个清晰的高声压区域焦斑尺寸通常与阵列全孔径对应的衍射极限相当。偏转聚焦更有意思。把焦点放到(x_f 5 mm, z_f 25 mm)相位延迟表就不再对称右侧阵元到焦点距离整体更近所需延迟更大左侧阵元延迟更小。表面上看焦点会从阵列法线方向向左侧偏转一个角度。这也叫电子偏转聚焦。还可以做一个参数化扫描让焦点横向位置x_f从-8 mm变化到8 mm观察声束偏转能力与副瓣的变化。你会发现在偏转角度过大时栅瓣会出现——这就是物理阵元间距过大导致的周期性干涉效应半波长间距准则的重要性从这里看得非常透彻。4.3 用瞬态仿真看聚焦的动态过程频域仿真给出稳态声场但超声波换能器实际工作多以短脉冲模式运行。如果想直观看到声波从阵元发出、互相干涉、在焦点汇聚的过程就需要瞬态仿真。瞬态仿真中相位延迟被替换为时间延迟t_n (d_n - d_ref) / c。激励信号用高斯包络正弦脉冲中心频率1 MHz脉宽2-3个周期。在COMSOL中这个时间延迟通过平滑阶跃函数或在激励边界上使用延迟操作来实现——简单做法是通过事件或分段函数给不同阵元定义不同的激励起始时间也可以直接给每个阵元施加时移后的信号函数。瞬态仿真的计算量比频域大得多。时间步长取激励频率周期的1/20到1/40才能保证精度即0.025 μs左右总时长要覆盖声波从换能器到焦点再加一段约50-80 μs总共数千步。配合二维模型和合理网格现代台式机跑几小时是正常的。瞬态结果里最直观的画面是声波动态演变从阵元表面一个个产生柱面波前各阵元的波前之间发生干涉在焦点处波前汇聚成一个小亮点然后又逐渐散开。这个动画对理解延迟聚焦的物理机制帮助极大也给写论文、做汇报提供了很好的素材。5. 网格、边界与求解器决定仿真成败的三个细节5.1 网格N个波长法则声学仿真有一条铁律一个波长内至少要保证5个二阶单元线性单元要10个以上。低于这个密度数值色散会让声速变慢、相位错乱焦点位移、声压下降都是必然结果。实际操作中我习惯把声场区域网格上限设为λ/6到λ/8这比理论最低要求再保守一些。代价是网格数上升但换来的相位精度非常值得尤其是焦点附近的干涉场对相位误差极其敏感。网格稍粗干涉极大值就可能减弱甚至消失。PML区域不能使用与传播域相同的自由网格。PML需要用映射网格做成规则的层状结构厚度方向上划分至少5-8层并按照PML坐标变换的要求进行拉伸。COMSOL的PML域会自动处理坐标变换但网格层数不足会导致吸收不良边界反射回来的波会叠在焦点场里形成驻波干扰。5.2 边界条件换能器背部与阵元间隙很多初学者只关注换能器辐射面的边界完全忽略阵元之间的间隙和模型外边界。阵元间隙的建模方式会影响阵列的实际有效孔径。当阵元表面振动时相邻阵元间的切缝区域如果设置为硬声场边界法向速度为零模拟的是阵元间切缝处无振动、刚性表面反射的情况这与实际压电阵元的刻槽填充物特性接近。如果间隙很窄小于λ/10可以直接忽略让阵元边界首尾相接也能得到基本准确的聚焦结果但会在相位计算和阵元个数上产生一点歧义。模型外围如果不是PML剩下的边界默认是硬边界声压反射全反射这会直接把边界反射波送回焦点区域形成高旁瓣。所以记住PML只放在外侧末端其他侧边界能放PML也放PML或者至少要设圆柱波辐射或平面波辐射条件做一类近似吸收。5.3 求解器选择的经验之谈COMSOL的默认求解器基于迭代求解在声学频域问题里尤其是含PML和压电耦合时经常不收敛。直接换用直接用MUMPS直接求解器绝大多数情况下都能稳定收敛。原因很简单频域声学矩阵是复对称矩阵直接求解法的鲁棒性明显优于迭代法。模型规模大时直接求解器内存占用会很高。如果矩阵自由度超过几百万可以考虑把模型改成二维周期性、或者用波束包络法如果问题满足慢变包络近似条件。但常规超声聚焦仿真中二维模型自由度基本在几十万量级直接求解器完全能胜任。瞬态仿真建议用BDF向后差分公式或广义α方法COMSOL默认的BDF阶数自适应即可。如果在瞬态中遇到高频振荡不稳定把最大BDF阶数降低到2同时加密时间步通常能解决。6. 我踩过的坑聚焦仿真常见问题排查实录6.1 焦点位置与理论值偏离太大最常遇到的问题是计算得到的焦点位置和理论焦距对不上。排查顺序如下。第一检查激励频率和声速设置。COMSOL中如果介质材料密度或声速设置错误波长就会偏相位计算用的λ与仿真实际的λ不一致焦点就会偏移。第二检查网格密度。如果网格过大导致数值色散等效声速低于真实值焦点会往远处推。这是非常隐蔽的原因因为从声压云图上你不容易察觉色散但焦点位置能精确暴露它。第三孔径不够大导致的衍射聚焦。如果口径与波长之比小于10焦点明显比几何焦距近且更弱。这与理论公式描述的几何焦点脱离有关。这种情况下不是仿真错了而是物理本身如此。6.2 相控阵声场看起来没有聚焦阵列模型算完声场一片平直或者出现多个对称的高压区而不是一个清晰焦点。首先确认相位符号对不对。COMSOL频域中时间因子是e^{iωt}还是e^{-iωt}与相位延迟的加减号相对应。如果用错符号阵元波前在焦点处正好反相叠加后不是增强而是抵消声场会出现分裂或零值。其次检查阵元间距。如果阵元间距大于λ会出现栅瓣声场中出现多个强度相当的伪焦点看起来像聚焦失败。解决方案是加密阵元减小间距或者在允许范围内降低频率。还有容易忽略的一点阵元宽度。阵元如果过宽每个阵元自身的辐射方向性会变窄各阵元的辐射场只在正前方有能量即使相位调整正确旁瓣也低不了。把阵元宽度控制在0.5λ以内获得较宽的单元方向图是相控阵正常工作的重要前提。6.3 PML反射导致焦点区域干涉条纹声场云图中焦点附近出现间隔为半个波长的明暗条纹——这是典型的驻波迹象。检查PML与传播域的界面是否紧密连接PML层数是否足够PML外侧边界条件是否设置正确。还有一个容易忽略的问题是PML厚度如果小于两个波长对斜入射波的吸收效果严重下降。解决方法是加厚PML到5-8个波长并把PML的层状网格数加到10层以上。6.4 瞬态仿真计算时间过长的妥协方案瞬态仿真如果时间步长过小、网格过密计算时间会膨胀到不可接受。我的妥协策略是先用粗网格和较大时间步做快速试跑确认声场形态和焦点位置大致正确然后用细网格、小时间步做最终确认。另一个技巧是只在焦点附近区域局部加密网格其他地方保留相对粗网格这样既保证焦点场精度又控制总体自由度。6.5 快速自查表检查项检查标准频率/波长用水声速1500 m/s重新核算频率对应波长网格密度最大单元尺寸≤λ/5二阶、≤λ/10线性阵元间距≤λ/2防止栅瓣相位符号相位延迟表达式与COMSOL时谐因子一致PML厚度至少2-5个波长厚度方向层状网格≥8层边界激励单位法向加速度单位m/s²不能直接用速度代替求解器频域用MUMPS直接求解不收敛时先排除网格负元7. 后续还能怎么玩从聚焦到更多声学仿真场景前面这些流程跑顺之后你可以在此基础上做很多扩展把压电效应完整加回来用晶片电压激励扫频看阻抗谱和谐振频率附近的声场变化或者把介质换成有衰减的软组织模型在压力声学中加吸收系数研究聚焦深度和热沉积又或者把线阵升级成二维面阵实现三维空间的定点聚焦。我自己做下来的体会是超声聚焦仿真的难点不在于某个具体操作而在于对整个声场干涉机制的直觉建立。COMSOL图形化的声场分布能帮你把抽象公式变成直观画面但前提是你得把模型做得对、边界条件设得准、网格剖得够。多花时间在物理量纲检查和网格独立性验证上远比反复调求解器参数有意义。还有一个我常用的偷懒技巧做参数扫描时先在频域对单个位置跑通再把焦点位置设成参数用COMSOL的参数化扫描一次性算完十几个焦点。这样能得到一张焦点可移动范围图对相控阵系统的偏转能力评估非常直观也方便给团队或客户展示。如果你正准备用COMSOL做超声换能器或相控阵的聚焦仿真建议从最简单的凹面换能器频域模型开始跑通流程后再逐步加阵列、加多物理场耦合、加瞬态分析。这个路径踩坑最少也能最快帮你建立信心。