超声相控阵聚焦仿真：COMSOL频域参数化模型实现与光声应用拓展

通常来说，超声相控阵聚焦仿真不算一个特别冷门的领域，但网上能直接跑通的完整参数化模型并不多。前阵子正好有个项目需要评估相控阵在不同焦点深度、不同阵元数下的聚焦声场，我用COMSOL搭了一个频域求解的有限元模型，不需要每次改几何、改边界，所有核心参数都可以在全局参数表里直接调整。模型跑通之后，又从超声扩展到了光声激励的部分应用，整个过程下来积累了不少细节经验。我把这套模型的构建思路和实操过程整理出来，希望能给正在做超声相控阵、声场聚焦仿真或者光声信号分析的同学一些参考。

这个模型解决的核心问题，说白了就是两件事：一是相控阵各阵元的相位延迟到底怎么算、怎么加进边界条件里；二是频域求解下如何得到一个稳定可靠的聚焦声场，而不是被网格、边界反射、求解器设置这些细节坑到怀疑人生。如果你用COMSOL做声学仿真有一段时间了，或者正准备从零开始搭一个相控阵模型，这篇文章可以帮你少走不少弯路。

1. 项目概述与整体设计思路

1.1 超声相控阵聚焦的核心逻辑

超声相控阵的原理如果用一句话概括，就是多个独立小阵元按规则排布，各自发射同频率、可独立控制相位的超声波，在目标位置处因为各阵元到该点的声程不同，通过预先补偿相位差异，让所有阵元在该点产生的声压同相叠加，形成一个局部增强的焦点声场。

打个比方，光学透镜聚焦靠的是曲面形状带来的光程差，让一束平行光汇聚到一个点上；相控阵是反过来的思路——不用透镜形状，而是用电子延迟主动制造出一个虚拟的"声透镜"。只要修改各个阵元上的相位延迟分布，焦点就可以在空间内自由偏转、扫描，甚至同时形成多个焦点。这也是相控阵在无损检测领域取代机械扫查的核心原因，电子扫描速度远快于物理移动。

在本模型中，我采用线性阵列布局：N个阵元等间距排布，每个阵元宽度为w，相邻阵元中心间距为pitch（节距），焦点位置用空间坐标（focus_x, focus_y）定义。第i个阵元所需的相位延迟由下式给出：

Δφ_i = 2πf · (r_i - r_ref) / c

其中f为发射频率，c为介质声速，r_i是第i个阵元中心到焦点的距离，r_ref是参考距离（通常取阵列中心阵元到焦点的距离）。物理含义很直接：把各阵元到焦点的路程差换算成时间差，再换算成相位差，发射时把这个相位差反向补偿掉，就能让所有阵元在焦点处同步到达。

把这条公式写进COMSOL的全局参数和边界条件表达式里，就实现了"参数可任意改变"的核心需求。改频率、改阵元数、改焦点位置，只需要修改参数表里的数值，不用动几何和物理场设置，模型的重用性非常高。

1.2 为什么选择频域而不是时域

这是整个模型选型里最值得展开讲的一点。很多做超声仿真的人上来就用瞬态研究，因为直觉上觉得超声波是脉冲形式的，瞬态更贴近实际。确实，如果关注的是回波信号的时域波形、脉冲宽度、多次反射这些信息，瞬态是绕不开的。

但如果你要评估的是聚焦声场的空间分布、焦点尺寸、焦点增益这些指标，频域方案在效率和可靠性上全面占优。

频域求解的核心假设是声场为稳态简谐变化。在COMSOL压力声学模块中，频域计算求解的是亥姆霍兹方程，时间项被完全消去，最终求解的是空间上的复数声压分布。这么做有两层好处：

第一，求解规模大幅下降。频域模型是稳态问题，不需要在时间轴上迭代，算一次解就用一个频率步进的代价，资源消耗比瞬态低一个量级。尤其在做参数扫描时，频域模型的优势非常明显。

第二，结果解读更直观。复数声压直接给出幅值和相位信息，焦点声压幅值、焦点偏移、旁瓣位置都可以直接量化，不需要做傅里叶变换。

代价也是存在的：频域模型无法直接模拟脉冲激发和宽带信号。如果后续要做脉冲回波检测模拟或宽带优化，一个可行的折中方案是在频域下做多个频点的扫描，再把多个单频结果合成宽带响应。这个思路在工程上被广泛使用，而且可以为后续时域模型提供可靠的初值。

2. 模型构建核心细节与物理场配置

2.1 几何建模与阵列参数化方案

我建的模型是二维的。这个选择是基于实际需求考虑——相控阵聚焦声场的评估，二维模型已经能给出焦点位置、焦斑尺寸、声场分布这些核心指标，计算速度快，调试方便。三维模型的建模逻辑完全一致，只是把阵列从单排变成二维面阵，但计算资源需求会涨很多，建议先把二维跑通再扩展。

几何结构分为三部分：底部是阵列区域，中间是声传播介质域，最外层是完美匹配层（PML）。为了减少不必要的自由度，PML的厚度通常取2~3个波长。矩形计算域的尺寸要根据焦点深度来定，一般要求焦点距离边界至少5个波长，否则边界反射会影响焦点的声压数值。

参数化是这个模型的关键。我在COMSOL的全局参数节点中定义了以下变量：

阵元的几何不再用单独的矩形绘制，而是通过"阵列公式"直接在边界上划分。具体做法是：在几何建模中创建一条长度为 N*pitch 的线段，作为阵列所在的边界。在设置法向加速度边界条件时，利用COMSOL的分段函数或循环表达式，把这条边界按 pitch 分成 N 段，每一段施加各自的相位延迟。

这里有个实操提示：不建议用布尔操作把每个阵元切成单独的小矩形，因为一旦阵元数变化就要重新建模。更灵活的做法是用一个"阵元索引表达式"配合边界选择的分割逻辑，把一条完整边界动态分成N段。COMSOL的"分区"功能可以基于距离范围自动选择对应的边界段，后续修改N值时几何不需要重建。

2.2 材料定义与声学边界条件

材料定义相对简单，水介质（声速约1500 m/s、密度1000 kg/m³）是最常见的默认选择。如果你模拟的是固体中的超声检测，压力声学模块就要换成弹性波模块，这个后面在问题排查部分会提到。这里先聚焦流体内的压力声学模型。

物理场方面，我使用的是COMSOL的压力声学模块（频域），核心因变量是声压p。方程形式为：

∇·(-∇p/ρ_c) - (ω²/ρ_c²)p = 0

其中ω = 2πf，ρ_c是介质密度，c是声速。这个方程是亥姆霍兹方程的标准形式，COMSOL内部已经处理好了。

阵元激励通过"法向加速度"边界条件施加。关键点在于，每个阵元段上的法向加速度表达式为：

a_n = A0 · exp(i·Δφ_i)

即每个阵元施加的激励幅值相同（A0），但相位不同（Δφ_i）。A0是激励幅值，可取1 m/s²作为归一化参考，方便后续对比焦点增益。

其他边界的处理是频域模型的一个重要细节。阵列以外的边界如果不做处理，声波会反射回来，干扰焦点附近的声场。我在计算域外围加了PML完美匹配层，这也是频域声学仿真的标准做法。PML区域在COMSOL中通过选择物理场设置中的"完美匹配层"节点来实现，需要指定PML在几何中的位置和厚度。

边界条件的完整清单如下：

这里有一个很多人容易忽略的点：如果阵列边界和计算域边界共用一条线，PML可能会把阵元的激励也吸收掉。我的做法是让阵列边界位于计算域内部一侧，不与PML重叠，确保激励边界和吸收边界在几何上是分开的两部分。

3. 参数化扫描与聚焦效果实操过程

3.1 相位延迟计算与参数化实现

相位延迟的实现是这个模型最容易出错的地方。直接用全局参数写一个包含循环的表达式并不方便，我采用的是COMSOL的"解析函数"功能来定义相位分布。

具体操作路径是：全局定义 → 解析函数，创建一个函数 phase_delay(x)，输入是阵元位置坐标 x，输出是该位置处的相位延迟值。函数表达式基于分段计算逻辑，先用 floor 函数算出该坐标点对应的阵元索引，再根据焦点坐标计算该阵元中心到焦点的距离：

idx = floor((x - array_start)/pitch) + 1

x_center = array_start + (idx - 0.5)·pitch

r_i = sqrt((x_center - focus_x)² + focus_y²)

r_ref = sqrt(((array_start + (N-1)·pitch/2) - focus_x)² + focus_y²)

Δφ = 2πf·(r_i - r_ref)/c

几个需要注意的细节：

第一，解析函数的参数类型要设为"按坐标"（x），变量名用x即可。COMSOL会自动传入边界上的空间坐标值。

第二，floor 函数在COMSOL中是合法的，但要注意类型转换。建议用 floor(x) + 1 得到从1开始的阵元索引，避免索引从0开始导致分段偏移。

第三，如果阵列中心不在坐标原点，公式里的 array_start 要单独定义一个参数。我的默认模型中，阵列中心在 x=0 处，所以 array_start = -(N-1)·pitch/2。

第四，COMSOL解析函数里的单位问题很容易踩坑。如果几何用的是毫米（mm），而声速单位为 m/s，频率单位为 Hz，那么计算 r_i 时要把几何坐标转换成米。一个简单的方式是在解析函数中显式乘以单位换算系数，或者在几何建模范式里就统一用米。

边界条件中的相位表达式为 A0·exp(i·phase_delay(x))，其中 i 是虚数单位。COMSOL支持复数表达式，直接写即可。这里我建议把解析函数的结果命名为 phi_delay，然后在法向加速度边界中用 ampexp(iphi_delay(x)) 这种形式，代码可读性好，后期排查也方便。

3.2 求解器配置与网格划分

网格划分在声学仿真里直接决定结果可信度。频域压力声学对网格的基本要求是：最大网格尺寸不超过波长的六分之一到十分之一。对于500 kHz、水中声速1500 m/s的情况，波长约为3 mm，因此最大网格尺寸控制在0.3~0.5 mm。

网格剖分策略上，我采用的是自由三角形网格，并为整个计算域设定统一的"最大单元尺寸"限制。虽然均匀网格会增加一些不必要的自由度，但在相控阵聚焦这类问题中，焦点位置相对集中，并不提前确定哪个区域需要加密，所以均匀网格反而是最稳妥的方案。如果你需要精确求解焦点附近的声压梯度，可以在焦点附近用一个圆形域做局部细化，最大网格尺寸再减半。

COMSOL的网格节点中有"大小"功能，选择"预定"参数中的"细化"级别，然后手动修改最大单元尺寸。对于我的默认模型（计算域约80 mm × 60 mm，含PML），设置最大单元尺寸 0.4 mm，可以得到约12万自由度的模型，普通电脑两三分钟就能算完。如果发现内存压力较大，可以把网格放宽到0.6 mm，对焦点声压数值的影响一般在1%以内，但计算时间会明显下降。

求解器选用频域默认的稳态求解器即可。COMSOL通常会自动推荐直接求解器（如MUMPS或PARDISO），对于这种中等规模问题直接求解器稳定性好，不需要手动干预。如果你加了很多频率点做扫描，建议开启"辅助扫描"，COMSOL会自动对每个频率点求解，并把结果保存在一起。

参数扫描有两种方式。一种是用COMSOL的"参数化扫描"研究步骤，扫阵元数量N或焦点深度 focus_y；另一种是做单频求解后手动修改参数重新计算。前者适合批量对比，后者适合快速验证。我的经验是先快速跑一个默认参数组确认设置无误，再开批量扫描，否则一次参数扫描设置错误会浪费大量机时。

3.3 结果后处理与焦点评估

求解完成后，结果评估的核心是焦点处的声压强度和焦斑的空间分布。COMSOL默认生成的声压场图能直接看到波前和聚焦效果，但定量的焦点评估需要自定义一些全局计算表达式。

我习惯定义一个"焦点增益"表达式：

G = abs(p_focus) / abs(p_ref)

其中p_focus是焦点处的复数声压，p_ref是阵元表面附近某参考点处的声压值（也可以定义为无聚焦时阵列前方同一深度处的声压）。这个增益值可以直观反映相控阵聚焦带来的声压提升倍数。

焦斑尺寸的评估推荐使用"3 dB带宽"概念。具体操作时，在焦点深度处（y = focus_y）画一条水平截线，输出声压幅值沿x轴的分布，用"派生值 → 一维绘图组 → 线图"来绘制。然后测量声压幅值下降到峰值一半（-6 dB，即幅度下降一半；或者-3 dB，即强度下降一半，这里取决于你用的是声压幅值还是声强）对应的横向宽度，这就是焦斑直径。

COMSOL的"截线"功能可以直接定义一条从 (x_min, focus_y) 到 (x_max, focus_y) 的直线，求解器会自动计算该线上所有点的声压值，导出数据后可以自己用绘图工具或脚本处理。

一个我经常用的操作是：在后处理中建立一个"表达式" p_abs = abs(p)，然后在派生值中选择"全局计算"，输入 aveop1(p_abs)（其中aveop1是定义在焦点处的平均值算子），可以直接读取焦点声压。或者在焦点位置定义一个"点评估"，计算该点复数声压的模。

如果做参数扫描，可以添加一个"全局评估"节点，勾选焦点声压和焦斑尺寸作为要扫描输出的表达式。这样一次批处理下来，数据表里直接给出每个参数组合对应的焦点声压和焦斑尺寸，非常高效。

4. 从超声到光声的拓展应用

4.1 光声效应与模型扩展思路

标题里提到的"支持超声、光声"是这个模型的一个重要卖点。光声成像的基本原理是：脉冲激光照射生物组织，组织中的吸收体（如血管）吸收光能后迅速升温，产生热弹性膨胀，进而激发出宽带超声波。这些超声波被体表的超声换能器接收后，通过重建算法可以反演吸收体的空间分布。

从数值仿真的角度看，光声问题可以分解为两个步骤：先求解光在介质中的吸收分布，再把吸收分布作为热源项，求解热弹性方程产生超声场。这两个步骤在COMSOL中可以分别实现，也可以耦合在一个模型中。

对于已有的相控阵聚焦模型，扩展到光声方向最简单的方式是在声源阶段做替换：不再用相控阵边界激励产生聚焦声场，而是在吸收体位置放置一个体积源（如单极/偶极声源），模拟光吸收产生的初始声压分布。声源与介质满足相同的波动方程，只不过源项的形式不同。

具体操作时，在压力声学模块中添加"声学边界条件"或"体积源"节点，将源区域定义为吸收体区域。初始声压分布可以通过"压力"边界条件或"法向加速度"来设置。对于脉冲激光激发产生的宽带声信号，频域模型处理起来会有一些局限，因为光声信号是宽带的，包含频率成分很丰富。但这并不妨碍用频域模型做特定频点的光声声场分析——比如评估某个中心频率下的声场聚焦和传播特性。

4.2 光声模型在COMSOL中的具体实现方法

实现光声拓展时，我采用了一种相对简洁的方式：在压力声学模块中定义一个"柱状体热源"（对应搜索词里的"comsol施加柱状体热源"），通过这个热源模拟吸收体受激光激发产生的膨胀。

具体来说，在物理场中找到"压力声学"中的"源项"，在源项表达式中写入一个高斯分布的体积源，中心位于吸收体位置，源强度与光吸收能量密度成正比。频域下这个源项的表达形式为：

Q = Q0 · exp(-((x-x0)² + (y-y0)²)/(2σ²))

其中 (x0, y0) 为吸收体中心，σ 为吸收体半径，Q0 为源强度。这个表达式实现了光声信号从吸收体区域向外传播的激励。

需要注意的是，频域求解时这个源项是复数形式。如果模拟的是脉冲激光，可以用多个频率点的扫描来近似脉冲频谱，然后在后处理中合成峰值声压。扫频点数越多，合成的时域波形越接近真实脉冲。

如果后续想把光声和超声相控阵结合起来（比如光声信号用相控阵接收），可以在同一个模型中同时设置两类源：一类是相控阵发射边界，一类是吸收体体积源。计算结果是两者的叠加，可以分别显示两类源单独作用的声场，也可以通过叠加得到混合激励的结果。

这里要提醒一点：光声模型中的频率范围通常比超声检测更宽。超声检测的中心频率可能只有500 kHz~5 MHz，但光声信号往往覆盖几MHz到几十MHz的频率范围。这意味着在频域扫描时，最高频率对应的波长可能很小，对网格尺寸的要求比纯超声模型严格得多。所以在做光声频域扫描时，建议用最高频率来确定网格尺寸，而不是用中心频率。

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

5.1 模型不出图的排查顺序

COMSOL声学仿真中，最让人头痛的问题就是明明算完了，结果图却是空的或者全零。按照我的排查顺序来，绝大多数情况都能在五分钟内定位：

第一，检查法向加速度表达式的复数形式。这是最容易出问题的地方。如果激励表达式没有包含虚数单位，比如只写了位相（相位）而 amplitude 不是复数，COMSOL默认只解实数域，声场会被错误地衰减或完全抵消。正确写法是 A0 · exp(i·φ)，确保激励为复数。

第二，检查网格和波长的相对关系。如果最大网格尺寸大于波长，物理场就无法正常传播，结果图上会出现严重的不规律模式或大片噪声。这通常是因为改了频率后忘记调整网格尺寸。建议把网格最大尺寸定义为 λ/8（即 c/f/8）这样的参数表达式，频率改变后网格自动更新。

第三，检查PML是否覆盖到激励边界。如果PML层和阵列边界重叠，激励能量会被PML直接吸收，焦点处几乎看不到声场。检查几何中PML区间和阵元边界区间是否有几何交集。

第四，检查焦点深度是否在计算域范围内。如果焦点坐标超出计算域，相位延迟计算得到的是一个发散的波前，焦点不存在。在参数化扫描时要特别注意焦点坐标和计算域尺寸之间的约束关系。

5.2 相位偏转后焦点偏移问题

这个是我在调试中遇到最多的问题。有时候相位延迟明明按公式算好了，但后处理中焦点出现的位置和预设位置有偏差，尤其是偏转角度比较大的时候。

原因是多方面的，最常见的是几何模型使用了毫米单位而声速和频率用了国际单位，导致公式中的 r_i 和 λ 量纲不匹配。COMSOL默认单位处理比较灵活，但如果你在解析函数中手动写了数字，容易出问题。我建议集中在一个表达式里做单位换算，例如读取坐标后统一除以1000转到米。

另一个原因是相控阵远场近似。公式 Δφ_i = 2πf·(r_i - r_ref)/c 是精确的（只要焦点不是特别近），但如果推导时用了平面波近似（比如用 sinθ 近似），在偏转角度超过15°以后误差会明显增大。所以建模时建议用距离公式直接算 r_i，不要用 sinθ 近似。

还有一个原因是阵元边界的位置选择。因为每个阵元是有宽度的，相位延迟应该以阵元中心坐标为准，而不是边界的起点或终点。如果用的坐标偏向阵元边缘，焦点位置会产生在一个阵元宽度以内的偏移。这在小阵元间距微如0.5 mm时会明显影响结果。

5.3 导入外部CAD文件时的常见警告处理

搜索词里有人问"solidworks另存为step后导入comsol有很多警告"，这个我在做三维扩展时也遇到过。一般来说，COMSOL导入STEP文件时的警告主要集中在几何缺陷：窄面、短边、尖角、微小曲面等。这些几何细节在仿真中通常不是必要的，但会让网格生成出现大量退化单元。

最省事的处理方法是导入后用COMSOL的"几何修复"工具，它会自动检测并消除大多数短边和窄面。如果修复后依然有警告，可以手动忽略那些不影响求解区域的微小瑕疵。重点关注的是：导入的几何是否包含多余的面（例如装配体的螺栓孔），如果有，直接用"删除实体"或"忽略面"把它们去掉。

还有一个经验：在CAD软件导出STEP前，尽量把模型简化到只保留必要的求解域。做仿真不是为了展示CAD模型，而是为了得到物理场分布。每个多余的圆角、倒角都会增加几何复杂度，导致网格数量成倍增加，但这些倒角对声场计算没有任何影响。

5.4 其他常见问题速查表

5.5 工作平面与移动网格的使用建议

在做三维模型或需要考虑结构变形时，"工作平面"和"移动网格"这两个功能非常实用。工作平面的作用是在三维模型中定义二维草图面，用来放置相控阵阵元阵列、创建规则几何、或者作为切割后处理的参考平面。我在三维相控阵模型里就是用工作平面来画阵元阵列，然后拉伸或旋转成实体。

移动网格则适用于需要考虑介质运动或结构振动的情况。比如换能器表面在振动时，声场计算要跟着边界一起变化，这时可以把换能器边界设成移动网格条件，让网格跟随边界位移变形。但值得强调的是，频域声学模型中如果只关注小振幅声学波动，通常不需要移动网格——直接用固定网格下的压力声学方程就能获得足够精度。只有在结构大位移或流固耦合强烈的情况下才需要开启移动网格，比如换能器膜片大幅振动的非线性仿真。

如果你只是做标准相控阵聚焦评估，我建议不要动这两个功能，保持模型简单，把精力放在物理场的准确设置上。

6. 模型文件的组织与后续扩展建议

6.1 参数管理的最佳实践

模型跑通之后，参数管理的好坏会直接影响后续使用的体验。我在COMSOL中采用的是一种分层参数管理方式：

第一层是全局参数，只放最常用的核心参数：频率、声速、阵元数、间距、焦点坐标。这些参数在使用界面中就能直接修改，不用进入模型源码。

第二层是派生参数，用全局参数作为输入，经过公式推算出次级参数。比如波长、PML厚度、最大网格尺寸、单阵元宽度等。这些参数一般不直接修改，而是随全局参数自动更新。

第三层是解析函数和边界条件中的表达式。这些实现了模型的核心物理逻辑，尽量保持稳定不变。

这样的分层结构让模型我可以放心地交给别人使用而不用担心被改坏——只要你把第一层参数控制在合理范围内，后面的计算逻辑会自动适应。

6.2 从二维到三维的扩展路径

如果你需要三维模型，建议不要直接重新建模，而是在二维模型的基础上做扩展。具体做法：

先确认二维模型的主要参数和结果指标都是符合预期的；然后在COMSOL中新建一个三维组件，把二维的几何、物理场、网格和求解设置整体复制进去；把阵列几何从一条线扩展成一个面（比如16×16的面阵），把焦点坐标从（focus_x, focus_y）扩展为（focus_x, focus_y, focus_z），相位公式不变。

唯一的难点在于计算资源。三维模型的自由度数通常比二维多2个量级，直接求解可能很吃内存。我建议在三维模型中适当放宽网格要求（最大单元尺寸从λ/8放宽到λ/6），或者利用对称性只计算四分之一区域，在模型对称面上设置对称边界条件即可。

6.3 COMSOL生态工具与数据交互

最近COMSOL也引入了MCP服务（Model Context Protocol），让外部AI工具可以直接调用COMSOL的模型数据和仿真结果，这对自动化批处理和跨工具协同很有价值。如果你经常做参数扫描和数据后处理，可以了解一下这个方向。坦白说，我在这个方面也就刚入门，但感觉它能在重复性工作里省下不少时间。

另一个实用技巧是数据导出。在做完频域扫描之后，我通常会直接把声压场数据导出成文本或CSV格式，再用Python或MATLAB做进一步的信噪比分析和图像重建。COMSOL自带的导出功能可以按网格节点、按截线、按体积分输出数据，用"数据集"选择对应的频率点或参数组合，就能得到特定条件下的声压分布数据。这个流程在做光声图像重建时的价值很大。

结语：实操中的一点个人经验

整套模型从搭建到稳定运行，我花了大概三周左右的时间。回头来看，最花时间的其实不是物理场的设置，而是细节——单位换算、相位表达式的复数写法、PML与激励边界的几何关系、网格尺寸与频率的联动。这些细节任何一个不对，结果都会千奇百怪。

最后再分享一个小技巧：在正式跑大规模参数扫描之前，先拿一个粗网格模型做快速验证——网格尺度放宽到λ/5，用最低的阵元数和最简单的参数组合跑通全流程。这能让你在几秒钟内确认模型设置没有低级错误，再切换到精细网格做正式计算。这么做是效率最高的一种节奏，也能让我这种"粗心型选手"少跑好几轮冤枉路。