通常来说超声相控阵聚焦仿真不算一个特别冷门的领域但网上能直接跑通的完整参数化模型并不多。前阵子正好有个项目需要评估相控阵在不同焦点深度、不同阵元数下的聚焦声场我用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的全局参数节点中定义了以下变量参数名含义默认值f超声频率500 kHzc介质声速1500 m/sN阵元数量16pitch阵元中心间距0.8 mmw阵元宽度0.6 mmfocus_x焦点x坐标0 mmfocus_y焦点y坐标深度40 mmlambda波长c/fpml_tPML厚度3*lambda阵元的几何不再用单独的矩形绘制而是通过阵列公式直接在边界上划分。具体做法是在几何建模中创建一条长度为 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在几何中的位置和厚度。边界条件的完整清单如下边界位置边界条件类型参数阵元表面法向加速度a_n A0·exp(i·Δφ_i)以阵元索引分段阵列以外的刚性边界硬声场边界默认法向速度为零计算域外圈PML厚度约3个波长对称面如果有对称/理想声场可根据模型对称性简化这里有一个很多人容易忽略的点如果阵列边界和计算域边界共用一条线PML可能会把阵元的激励也吸收掉。我的做法是让阵列边界位于计算域内部一侧不与PML重叠确保激励边界和吸收边界在几何上是分开的两部分。3. 参数化扫描与聚焦效果实操过程3.1 相位延迟计算与参数化实现相位延迟的实现是这个模型最容易出错的地方。直接用全局参数写一个包含循环的表达式并不方便我采用的是COMSOL的解析函数功能来定义相位分布。具体操作路径是全局定义 → 解析函数创建一个函数 phase_delay(x)输入是阵元位置坐标 x输出是该位置处的相位延迟值。函数表达式基于分段计算逻辑先用 floor 函数算出该坐标点对应的阵元索引再根据焦点坐标计算该阵元中心到焦点的距离idx floor((x - array_start)/pitch) 1x_center array_start (idx - 0.5)·pitchr_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 要单独定义一个参数。我的默认模型中阵列中心在 x0 处所以 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 其他常见问题速查表现象可能原因解决办法计算结果显示声压全为零激励表达式没加虚数单位检查法向加速度表达式确保含 exp(iφ)焦点位置偏移单位不匹配或相位用平面波近似用统一国际单位用距离公式算相位PML区域出现数值波动PML厚度设置过薄增加PML厚度至3个波长以上网格剖分失败几何中存在极小边或窄面使用几何修复工具清除微瑕疵参数扫描结果趋势异常某参数改变了物理域尺寸但网格未更新把网格尺寸定义为频率/波长的函数表达式求解不收敛或内存不足网格过于细密适当放宽网格或改用更大内存/分批扫描导出后处理数据为空没有选择正确的派生值类型使用体/面/线选择而非点选择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用最低的阵元数和最简单的参数组合跑通全流程。这能让你在几秒钟内确认模型设置没有低级错误再切换到精细网格做正式计算。这么做是效率最高的一种节奏也能让我这种粗心型选手少跑好几轮冤枉路。