1. 项目概述与仿真目标

做超声相控阵仿真的人应该都有同感：单阵元或者一维线阵的模型跑起来很轻松，但一旦涉及到面阵、尤其是像7×7这种规模的二维阵列，情况就完全不一样了。这个项目选择COMSOL作为仿真平台，对7×7阵元三维相控阵聚焦探头进行全波仿真，核心任务有两个：一是准确计算焦点在三维空间中的位置，二是观察声场随着阵元激励时序变化的动态过程。

先说结论：COMSOL在处理这类问题上确实有先天优势。它不像传统有限差分法（FDTD）那样需要小心翼翼地控制时间步长和网格尺寸的配比，也不像解析法那样必须把阵元理想化成点源。COMSOL的声学模块里既有频域求解器也有瞬态求解器，配合压电耦合接口还可以直接把电学激励、机械振动和声波传播放在同一个模型里解，这对于研究相控阵聚焦是至关重要的——因为相控阵的本质就是“时序控制”，你必须看到不同阵元在不同时刻发出的声波在空间里怎样叠加、干涉，最后在目标位置形成主瓣和旁瓣结构。

7×7这个数值得单独说一下。它介于稀疏阵列和密集阵列之间，既不会像16×16那样带来极高的计算成本，也不会像4×4那样旁瓣水平高得让人没法看。7×7阵列在无损检测（NDT）和医学超声聚焦中都有实际工程参考价值，仿真的结果可以相对容易地往实际系统迁移。

这个项目适合谁参考？如果你在做超声换能器设计、相控阵聚焦算法验证、声场干涉分析，或者你是COMSOL声学方向的初学者想找一个难度适中的三维阵列案例，这篇文章里记录的建模思路、参数选择和踩坑记录应该能帮你省下不少时间。

2. 阵列设计与物理模型搭建

2.1 7×7阵列的参数确定

阵列设计是整个项目的第一步，也直接决定后面仿真结果的可靠性。我在建模时采用的参数如下：

• 阵元间距（pitch）d = 0.8λ，λ为水中声波波长
• 阵元宽度 a = 0.65λ，阵元高度 b = 1.2λ
• 相邻阵元间隙 = d - a = 0.15λ
• 工作频率 f₀ = 1 MHz，对应的水中波长 λ ≈ 1.5 mm
• 阵元总数 N = 49，阵列整体孔径约 5.6λ × 5.6λ

为什么让阵元间距取0.8λ而不是半波长或者一个波长？这里有个经典的折中。间距大于λ会在空间中产生栅瓣（grating lobes），也就是除了主焦点之外，在周期性的角度位置上还会出现多余的声波汇聚区域，这是相控阵聚焦最要避免的。间距小到半波长以下虽然能彻底消除栅瓣，但阵元之间的互辐射和机械耦合会明显增强，而且阵元尺寸过小意味着单个阵元的辐射阻抗下降，能量传输效率会变差。0.8λ是一个工程上很常用的折中值，它保证在±45°的偏转范围内不出现栅瓣，同时阵元的辐射面积也能维持得不错。

这里有一个COMSOL建模里的细节：阵元几何宽度和激励边界之间的关系。很多人建模时直接给整个矩形区域施加电压边界条件，这会高估阵元的有效辐射面积。我在模型里专门为每个阵元定义了一个压电材料区域（PZT-5H），只有这个区域的上表面才是电学激励面，声波从阵元表面进入水介质。这样设置虽然让网格量上升了一个量级，但得到的声场分布和实际探头的辐射特性更接近。

2.2 几何模型与材料设置

几何上用的是三维全模型——没有用对称边界条件来减半计算域。有人会问，7×7阵列在中心对称的情况下不是可以只建1/4模型吗？理论上确实可以，但实际操作中我不建议这么做，原因是相控阵的聚焦过程是动态的，不同的聚焦点在空间中的位置不同，每一时刻的声场相对于阵列坐标系根本不满足对称性。如果你用对称边界，就相当于强行给声场添加了理想反射边界，动态聚焦仿真的结果会产生偏差。

计算域的形状选的是球形，半径取30λ。这样做的原因有二：第一，球形域配合COMSOL的“完美匹配层”（PML）可以最大程度减少边界反射对焦点区域声场的干扰；第二，球形域可以保证声波从阵元出发到边界各点的传播距离一致，PML的吸收效率均匀。

材料和物理场的选择按下面这套来：

• 压电阵元域：PZT-5H，极化方向沿Z轴正向，使用压电效应多物理场耦合
• 液体域：水，声速1480 m/s，密度1000 kg/m³
• 结构域：仅压电陶瓷本体，不考虑背衬层和匹配层——这样做的原因是为了先关注声源本身的辐射特性，匹配层在后期可以加进去做对比研究

在压力声学接口下，液体域求解的是亥姆霍兹方程，但为了看动态过程，我使用的是瞬态求解器（时域求解），时间窗口覆盖了从激励开始到声场完全穿过PML的全过程。

2.3 网格划分策略

网格划分是整个三维声场仿真里最折磨人的一步。我的实测经验是：直接使用COMSOL默认的“流体动力学网格”并不理想，它对曲面边界加密过度，反而在阵元表面这种需要精确捕捉振速的地方网格不够密。

我最终采用的策略是分区控制：

• 压电阵元域：最大网格尺寸 λ_piezo/6，λ_piezo是PZT中声波波长约0.4mm，所以网格约67μm
• 阵元表面附近液体域：用边界层网格，第一层厚度25μm，共5层，增长率1.3
• 液体主域：最大网格尺寸 λ_water/6 = 0.25mm
• PML区域：扫掠网格，沿径向厚度方向8层拉伸分布

这套网格设置模型总自由度在450万左右，在64GB内存的工作站上可以流畅求解。如果内存紧张，可以把液体主域的最大网格尺寸放宽到λ/5，计算量能下降约三分之一，但焦点附近声压的峰值精度会有可感知的下降，尤其是对旁瓣电平的评估误差可能超过1dB。

3. 焦点定位与声场动态变化的分析方法

3.1 聚焦延迟法则（Focusing Delay Law）的实现

相控阵聚焦的原理通俗讲就是“分先后、同到达”。所有阵元发出的声波要在焦点处同时到达，那么离焦点远的阵元必须早发力，离焦点近的阵元晚发力。这个时间差就是聚焦延时。公式是：

Δt_n = (√(x_n² + y_n² + z_f²) - d_n) / c

展开说就是：第n个阵元到焦点的距离减去参考距离，再除以声速c，得到该阵元需要提前的时间。这里的x_n和y_n是阵元中心坐标，z_f是焦点深度。

在COMSOL中实现延迟法则有两种路径。第一种是在压电阵元的上表面施加“随时间变化的压力边界条件”或“位移边界条件”，在每个阵元的激励函数里加入一个时间偏移参数。第二种更高效，直接用“瞬态的压力声学边界”里调制每个阵元的法向加速度，让它按照指定的延时曲线启动。

我实际选择的是给每个阵元施加一个加窗的汉宁猝发信号（tone burst），中心频率1MHz，持续5个周期，信号初始相位按延时参数分配。之所以用汉宁窗而不是矩形窗，是因为矩形窗的频率旁瓣太高，会造成声场中出现明显的非期望频率分量，最终表现为焦点周围杂乱的干涉条纹，干扰对主焦点位置的判断。

3.2 焦点定位的计算方式

焦点定位在COMSOL里有两种互为验证的计算方法：一种是直接在求解结果中搜索声压幅值最大的空间点坐标，另一种是根据传播时间差反推焦点。搜索法需要注意一点：声压幅值最大点不一定等于焦点位置，因为瞬态过程初期，阵元近场的干涉会出现局部极大值点，这些点可能比真正的焦点声压还高。所以需要设置一个时间窗口，只在脉冲完全到达预期深度（按理想延时计算）之后的区间内搜索峰值。

我在后处理中使用了一个“域表面最大值”探针组，同时监控三个深度的截面（z = f - 2λ，z = f，z = f + 2λ），以及一个三维体最大值探针。这个做法的好处是可以实时看到焦点在时间维度的建立过程。从结果来看，在激励开始后约8.2μs，焦点位置的声压峰值开始形成，14μs左右达到稳态最大值，随后主瓣逐渐扩展并退出焦点区域。

3.3 声场动态变化的可视化与量化解构

动态声场的分析不应停留在“看动画”的层次。除了经典的声压分布云图动画之外，我做了三组量化提取：

第一组是焦点轴线上（Z轴）的声压幅值包络线。它告诉你聚焦增益的真实效果——焦点处的声压与阵元表面辐射声压的比值，也就是聚焦增益G_F。在本模型中，理论值G_F = N = 49（无衰减理想情况），但由于阵元间隙、有限孔径、水的声吸收等因素，仿真得到的G_F约为31.6。这个差距非常正常，不必为了追求49而强行加密网格或加大阵元占空比，实际工程中也达不到理论极限。

第二组是焦点平面（XY平面）的横向声压分布剖面。通过提取焦点所在截面的声压幅值，采用-6dB波束宽度来衡量焦点的大小。这次仿真得到的-6dB焦斑直径约为0.7λ（约1.05mm），在横向和纵向略有差异——横向稍窄、纵向稍宽，符合圆形活塞聚焦的典型规律。

第三组是时间-声压曲线。在焦点位置放置一个点探针，记录完整的时间波形。你会看到波形呈现明显的“先小后大再衰减”的包络：前半段是边缘阵元先到达的信号，中间是中心阵元后到达的信号叠加后形成的峰值，尾巴则是穿过焦点继续传播的衍射波。这个时间波形其实可以反过来验证延时法则是否正确——如果各阵元到达时间精确对齐，主峰幅度最大且波形对称。

4. 仿真结果对比与实际应用映射

4.1 焦点位置误差与偏转控制

把仿真得到的焦点坐标和理论预设坐标进行对比，本章节值得一提的发现是：当焦点位于Z轴正下方时，定位误差小于0.05λ；但当焦点偏转至轴外（例如偏转30°方向），焦点位置开始出现偏移，最大偏移量达到了0.18λ。原因并不复杂——大角度偏转时，阵元之间的声波传播路径差异变大，介质吸收导致各阵元到达焦点时的幅度不一致，等效于在孔径上施加了一个幅度锥削，这会让焦点位置略微向孔径中心方向回拉。

这在医学超声聚焦和工业无损检测里都是有实际意义的：如果你要在一定区域内进行扫描聚焦，不能想当然地认为“给延时就能指哪打哪”，需要对偏转角度带来的焦点偏移做预校正。

4.2 焦点大小与旁瓣水平的影响因素

我跑了几组对照组仿真来观察参数变化对焦点质量的影响：

这组对比明确反映了几个工程直觉：增大单个阵元面积能提升灵敏度，但也抬高了旁瓣；提高频率能显著收紧焦点，但穿透深度会下降；激励加窗虽然小幅牺牲主瓣幅度，但对旁瓣抑制的贡献非常显著，而且还能缩短脉冲的时间长度，提升距离分辨力。

4.3 从仿真走向工程应用

这套模型能输出的绝不止是几张漂亮的云图。如果你做的是医学超声聚焦（比如高强度聚焦超声HIFU），仿真得到的焦点声压分布可以用来估算焦区温升——把声压幅值作为热源项导入生物热传导方程，粗糙估算就能大致圈定热损伤范围；如果你做的是工业NDT，7×7阵列仿真的声场结果可以直接用于优化全聚焦成像（TFM）的虚拟聚焦算法，因为阵列的辐射声场数据本身就是成像算法的核心先验信息。

从实际应用来讲，这个模型的计算结果和实验室中用激光测振仪或水听器扫描实测的声场分布做对比，能发现一个规律：仿真中-6dB焦斑通常比实际测量小5%~10%。这是因为仿真里没有模拟阵元之间的电学串扰和机械耦合同样，PZT阵元边缘振动的非理想约束也会让实际旁瓣比仿真略高。做设计时把这个误差预留进去就可以。

5. 常见问题与排查技巧实录

首先要说的是求解稳定性问题。三维瞬态声场仿真最容易遇到的就是“计算发散”。现象通常是求解到某个时间步后声压值突然变成10^30量级或者NaN。遇到这种情况，90%的可能是网格尺寸和声波波长不匹配。COMSOL的瞬态声学要求最小网格尺寸不能小于λ/10，但不是说网格越细越好——网格越细，局部Courant数越小，数值色散愈发明显，波前会出现虚假振荡。一个稳健的配置是让液体域的网格尺寸等于λ/8到λ/6之间，时间步长交给COMSOL自动控制，但把“最大时间步长”手动设定约为0.25 / (f₀×10)，也就是每个周期约40步。这样既保证了收敛性，又在消耗适中的前提下维持了波形保真。

其次是PML参数的问题。很多人把PML厚度设为2-3个网格层，结果边界反射依旧明显——反射波从边界折返回计算域后，焦点区域的声压会出现异常的干涉纹路。我在三维域里用的PML厚度是最高声波波长的1.5倍，同时把PML的“缩放因子”从默认值稍作调整（比如设为2），吸收效果明显改善。想确认PML质量，可以放一个在边界附近方向的远场探针，观察反射波到达前与到达后的声压比，通常应低于-40dB。

第三个高频问题出现在压电耦合的瞬态仿真中。COMSOL的“压电效应”多物理场接口在频域求解时非常稳定，但在瞬态求解时容易出现电荷守恒异常。我排查后发现，根源是压电域的电极边界条件没设置好——使用“电浮置”或“零电荷”条件时，压电阵元的电容效应不能被正确模拟，激励电压和振动速度的相位关系因此发生偏移。解决办法是给每个阵元的电极面分配一个“终端”节点，并把终端类型设置为“电压”，给相邻阵列的一个公共电极指定为接地。这样每个阵元的电容、电荷、电流都能独立追踪，仿真结果和实际探头的电声测量数据吻合度显著提升。

最后记录一个看起来玄幻但实际发生过的低级错误：7×7阵列一共49个阵元，延时参数表有49行。有一次我把行号和阵元坐标搞错位了，导致同一行的坐标和延时对应不上，仿真出的焦点完全分裂，声场分布一团乱。检查了很久才发现思路是清醒的、数据是错的。之后我专门在COMSOL里做一个“坐标检查”步骤，把49个阵元中心坐标直接导出，用公式算一遍理论延时，再和Excel里准备的延时表逐行比对，确认无误后才提交求解。这个步骤看起来笨，但真的能救命——特别是阵列规模继续扩大的时候，人工核对坐标几乎不可靠。

针对多个阵元同时激励时计算资源吃紧的问题，我还可以分享一个折短的方案：如果你关注的是稳态聚焦特性（焦点增益、旁瓣电平等），可以改用频域求解器，在每个频率点求解一次复数声场，然后用后处理合成声压幅值分布。频域求解不需要长时间迭代，内存占用比瞬态低一个量级，而且PML的处理在频域下更加稳定。但代价是看不到声波传播的动态干涉过程——所以我的建议是：做方案验证用频域，做动态分析用瞬态，两者交叉验证是最理想的。

提示：在做瞬态声场分析时，千万不要忽略后处理的时间积分。声压云图动画虽然直观，但能量指标（声强、声能密度）才是评价焦点质量的真实标准。在COMSOL里可以从瞬态结果中计算一个周期的RMS声压，用RMS云图来评估焦点的有效区域，比直接看某一时刻的瞬时声压更能说明问题。

6. 经验总结与扩展建议

做这个7×7阵列相控阵聚焦探头的仿真，我个人的核心体会可以浓缩成三句话：阵列参数决定声场的“骨架”，延时法则决定声场的“灵魂”，而网格和求解器设置决定仿真能不能“活着”跑完。三者缺一不可，任何一环的疏忽都不会体现在报错信息里，而是潜伏在看似正常的声场云图之中。

从实用角度看，COMSOL的能力边界并不在于“能不能仿”，而在于仿真效率是否匹配你的迭代需求。表面上看，一次完整的三维瞬态仿真需要数小时，的确比解析模型慢得多，但当你要研究非规则阵列排布、异形阵元结构、声透镜耦合等解析手段搞不定的工程问题时，这个代价完全可以接受。

后续扩展方面，我计划在这套模型上做三件事：一是加入声学匹配层和背衬层，观察它们对带宽和脉冲响应的影响；二是把聚焦点按照一定的路径移动（比如扫描一条线段），记录各焦点位置的声压一致性，这可以直接服务于超声相控阵成像系统设计；三是考虑把生物组织或复合材料的非均匀声学参数加入传播域模型，分析像差的影响——这一步是从“仿真验证”走向“仿真预测”的关键转折，也是相控阵聚焦探头在无损检测和医学超声中真正体现工程价值的场景。