做超声检测和声学设计的朋友一定经历过这种尴尬新的相控阵探头还在论证阶段阵元数量、间距、频率都悬着打样试错的钱花得心疼实验台又排不上队。这时候COMSOL有限元超声相控阵聚焦仿真模型就派上用场了——把阵元参数全部参数化在频域里一次算出稳态声场顺手把焦点位置、焦点尺寸、聚焦增益这些指标拉出来整个方案可行性判断一两天就能落地。这篇文章不是把COMSOL官网案例抄一遍而是把实际建模过程中那些真正影响结果的细节、命令和坑整理出来。内容包括相控阵聚焦的频域求解逻辑、模型搭建里物理场和边界的设置思路、延时法则在COMSOL里的工程化实现、参数化扫描和后处理分析方法最后把常见的边界反射、内存爆炸、CAD拓扑报错这些问题一并排掉。适合正在做超声聚焦、无损检测换能器设计、水声工程仿真的工程师参考哪怕你刚接触COMSOL跟着这套思路也能跑通一套能用的聚焦仿真模型。1. 相控阵聚焦仿真的底层逻辑1.1 相控阵到底在做什么相控阵换能器的核心是一排独立的压电阵元每个阵元可以单独控制激励信号的延时——也就是相位。生活中可以这样理解一排人喊口号如果所有人同时张嘴声音是朝四面八方散的如果让最远的人先喊、最近的人后喊这个“错峰喊话”能让声音在一个预先指定的点上汇聚恰好像聚光灯。超声相控阵聚焦就是这个原理的工程版。设焦点坐标为F(x0, y0)第i个阵元的中心坐标为(xi, yi)声速为c。声波从第i个阵元传到焦点需要的时间是ti sqrt((xi - x0)^2 (yi - y0)^2) / c要让所有阵元发出的波同时到达焦点激励信号就需要按时间差做补偿。通常以最远阵元为基准取Δt_i t_max - ti也就是距离焦点越近的阵元越晚激励。把Δt_i换算成频域相位偏移就得到每个阵元的激励相位。这一整套延时分布业内叫“聚焦法则”是相控阵聚焦的灵魂。1.2 频域仿真为什么是工程主力COMSOL里做声学仿真最常用的有两类研究时域研究和频域研究。时域能捕捉脉冲波传播的完整瞬态过程看得到波包如何运动、如何反射物理画面很直观。但时域求解有严格的时间步长约束网格越细、域越大步数越多一个二维模型跑几个小时稀松平常三维模型更是灾难。频域仿真假设系统处于稳态正弦激励下求解的是线性方程组不存在时间步进问题。相控阵聚焦的典型工作状态是连续波或窄带脉冲激励频率基本固定在中心频率附近用频域求解在物理上完全合理。更重要的是频域天然适合做参数化扫描频率、阵元间距、焦点深度、阵元数量这些东西一变重新解一次线性系统就行不需要从头跑时间历程。工程选型时我的建议很明确如果只看聚焦效果、旁瓣水平、焦点尺寸这类稳态指标直接上频域如果非要研究脉冲波形畸变、非线性效应或者瞬态聚焦过程再考虑时域。1.3 有限元方法在这个场景里的不可替代性超声换能器声场也可以用解析方法算最经典的是瑞利-索末菲积分通过把每个阵元等效为点源叠加得到空间任一点的声压。解析方法的优点是快几十秒就能出来一幅声场图缺点是边界条件简化得厉害碰到复杂几何形状、近场区域、非均匀介质、多层结构就不那么好用。有限元方法把计算域离散成网格单元在每个单元内用形函数逼近声压场能处理任意几何形状和材料分布近场精度比解析法高得多。代价是网格要划得足够密矩阵求解要吃掉不少内存。现实中做相控阵聚焦仿真通常就是几何还算规整、但阵元排布和边界情况比较复杂有限元这种“一刀切”的通用求解器反而最省心。2. 模型搭建几何、物理场与边界设定2.1 几何建模2D还是3D对称性该用就用我一开始直接烧三维模型结果网格数量大得离谱一台32G内存的工作站跑一个8阵元阵列就到极限了。说到底相控阵聚焦仿真不一定非得3D。如果阵列是一维排布也就是阵元沿x轴排列、声场在xy平面内传播那用2D模型完全够用。COMSOL里新建“二维”组件计算域设成一个矩形水域阵列边界放在上边沿或者左边沿焦点落在域内部。2D模型网格量小两个数量级参数扫描时可以一口气跑几十组。二维模型还能进一步利用对称性。阵列排列和焦点都关于某条轴线对称时只建一半模型在对称轴上施加对称边界条件计算量立刻减半甚至减更多。COMSOL的“对称”边界条件在压力声学里对应的是硬边界条件也就是法向振速为零物理上等同镜像反射。三维模型也不是不能做适合阵列是面阵、焦点在三维空间里、需要考察横向和纵向聚焦分辨率的场景。实际经验是先用2D把方案跑通把参数趋势摸清再上3D做最终验证这条路径效率最高。CAD模型的导入偶尔会遇到“转换为CAD内核时不支持的拓扑”这个报错通常是STEP文件里存在重叠面、微小边或者退化几何。我的习惯是直接用COMSOL内置的几何工具重建声学域——反正就是一个矩形或圆柱水域阵列边界用矩形阵列操作生成几十秒就能搞定比在CAD软件里反复修补导出再导入省心得多。2.2 物理场接口压力声学、频域是默认选项COMSOL里做常规声学仿真物理场接口选择“声学模块 压力声学 压力声学频域”。这个接口求解的是亥姆霍兹方程也就是频域波动方程变量是声压p。在均匀介质中方程形式是∇·(−1/ρ0 ∇p) − ω^2/(ρ0 c^2) p 0其中ρ0是介质密度c是声速ω2πf是角频率。这个接口会求解复声压场实部代表瞬时声压在某一相位时刻的值虚部加入了相位信息。声学域的材料设置为水密度998 kg/m³声速1480 m/s。如果目标介质是其他液体或固体替换材料参数即可。这里要注意材料参数里的声速必须与你设置的频率、波长配套否则网格尺寸和PML厚度即使看起来合理实际却可能差出一个数量级。换能器的激励边界有几种处理方式。最简化的方法是把每个阵元边界设成“法向加速度”边界条件或者直接给阵元边界赋一个复声压边界条件。如果你不关心压电陶瓷内部电场和机械振动的耦合这就是最快的路子。如果要做完整的压电超声换能器仿真那必须加上“压电设备”物理场接口把PZT材料的压电耦合矩阵、电极边界全建出来网格和计算量会相应增大。2.3 边界条件别让声波“碰壁”无限大水域在有限元模型里只能截取有限区域截断边界如果处理不好声波会在边界上反射回来聚焦声场里就会混入一圈圈的假干涉条纹。这个干扰早期特别容易忽略因为图像上看起来像漂亮的波纹实际上全是假信号。COMSOL里处理开放声场的最标准做法是加完美匹配层也就是PML。在声学域的外围再包上一圈环状区域把这个区域设为“完美匹配层”的域条件。PML的原理是在内部对方程做坐标拉伸变换使进入该区域的波被指数衰减吸收理论上不产生反射。PML参数设置有两个关键点。第一是厚度建议至少一个波长稳妥一点取1.5到2个波长。第二是网格PML区域内网格应该尽量规整通常用映射网格或扫掠网格生成规则矩形网格这样可以避免不规则网格引起的局部数值误差。还有一个容易被忽略的点就是PML外面不需要再设任何边界条件程序默认处理为吸收。物理域和PML交界处的网格连续性要保证否则出现数值反射焦点位置的声压分布可能已经受了污染。2.4 网格尺寸怎么定频域声学仿真的网格规则可以用一句话概括每个波长内至少5到6个二阶单元。声速c、频率f和波长λ的关系是λc/f。以水声速1480 m/s为例不同频率对应的波长和推荐网格尺寸如下频率波长推荐网格每波长6单元250 kHz5.92 mm约1.0 mm500 kHz2.96 mm约0.5 mm1 MHz1.48 mm约0.25 mm在COMSOL网格节点里可以直接选择“声学”预定义网格大小软件会根据物理场接口自动估算波长并给出推荐尺寸。我自己的流程是从较粗网格开始试算逐步细化做一次标准的网格收敛性检查。如果两次加密后焦点声压变化小于1%网格量就是合适的。PML区域网格要单独处理。我的做法是在PML区域用映射网格划定结构化网格层内沿厚度方向划分4到6层即可太多太密只会徒增计算量。3. 聚焦延时法则在COMSOL里的工程化实现3.1 延时计算先算传播时间再转相位聚焦法则的工程实现在COMSOL里有几种套路先说最朴素的。假设有N个阵元沿x轴排列阵元间距为pitch第i个阵元的中心坐标为xi。焦点坐标是(x_f, y_f)。第i个阵元到焦点的距离为ri sqrt((xi - x_f)^2 y_f^2)对应的传播时间为ti ri/c。取所有阵元中传播时间最大值t_max第i个阵元的延时为Δt_i t_max − ti。在频域谐波分析中时间延迟Δτ对应的相位偏移是−2πfΔτ符号约定取决于COMSOL的时谐因子约定工程计算中你只需要保证各阵元间的相对相位差是对的。这样第i个阵元的激励就可以写成复值边界条件p_i p0 * exp(−j * 2 * pi * f * Δt_i)p0是单个阵元的激励声压幅值设置为1 Pa便于后处理归一化。3.2 用分段函数把延时映射到阵列边界每个阵元的坐标不同延时也逐元变化。在COMSOL里不需要为每个阵元单独建边界条件更聪明的办法是用一个分段函数把延时映射到空间坐标上。在“全局定义 函数”中新建一个“分段函数”变量设为x区间按阵元边界划分。比如阵元1的x范围是0到0.6mm此区间内函数值是Δt_1阵元2的x范围是0.6到1.2mm函数值是Δt_2依此类推。然后把换能器整个阵列边界设成一个统一的压力边界条件赋值为p0 * exp(−j * 2 * pi * f * delay_func(x))由于分段函数只在相应区间内取相应延时值就实现了各阵元独立相位控制。这个思路最大的好处是参数扫描时只要更新延时函数整个模型自动响应不需要动几何结构。阵元数量很多时也可以用外部CSV表格导入延时数据再通过“插值函数”映射到坐标上。听说有些项目里还有把MATLAB算好的聚焦法则直接灌进COMSOL的做法本质上都是一样。3.3 参数化扫描让所有关键参数活起来COMSOL的“全局参数”是参数化仿真的心脏。我把f0工作频率、pitch阵元间距、elem_count阵元数、focus_x和focus_y焦点坐标全部定义成全局参数。几何构建时阵列位置用参数驱动。二维建模时先用一个“矩形”表示单个阵元边界再用“阵列”操作x方向间距设为pitch个数设为elem_count。这样pitch和elem_count一变几何自动重建网格也跟着自适应。研究设置中在“研究1”里选择“频域”把频率设为f0。然后在研究步骤前添加“参数化扫描”把需要扫描的参数加进去比如扫f0从300kHz到700kHz步长50kHz同时pitch从0.4mm扫到1.0mm步长0.2mm。COMSOL会自动计算所有参数组合每个组合解一次频域方程不需要你手动干预。这里给一个参数表参考方便你直接抄作业参数名含义典型值扫描范围f0工作频率500 kHz250k–1 MHzpitch阵元间距0.6 mm0.3–1.2 mmelem_count阵元数量168–32focus_x焦点横向坐标0 mm−10–10 mmfocus_y焦点纵向深度30 mm10–60 mmp0阵元激励幅值1 Pa固定1参数化扫描跑完后所有结果都存放在不同数据集里。后处理时可以通过“数据集”下拉列表切换不同参数组合的结果非常方便。4. 聚焦效果分析怎么证明“聚焦成功”了4.1 第一眼判定声压幅值分布图最直观的分析是在“二维绘图组”里画声压幅值分布。这里的第一个坑是变量选择COMSOL频域求解得到的是复声压acpr.p直接用实部或虚部绘图会看到明显的不对称图案好像聚焦非常怪异。必须用幅值也就是abs(acpr.p)或者声压级Lp才能得到真实的声场包络。画出来后判断聚焦是否成功的三个视觉标准焦点位置是否有明显的亮斑也就是声压幅值极大值。焦点周围是否存在过大的旁瓣主瓣和旁瓣的亮度差是否清晰。声束是否沿着预期的轴线传播有没有明显偏斜。如果亮斑不在预设焦点位置先检查延时法则的符号和坐标定义是不是反了这是我之前踩过最多次的坑。4.2 量化指标这样算视觉判断只能用来粗筛写报告还需要量化指标。COMSOL的“派生值”菜单里可以直接计算这些第一焦点位置。使用“派生值 体最大值”或“面最大值”在“表达式”里填abs(acpr.p)让软件帮你定位声压幅值最大的坐标。如果这个坐标和预设焦点坐标偏差在可接受范围内说明聚焦法则设置正确。第二焦点尺寸。一般是提取通过焦点的直线上声压分布然后找主瓣下降6dB两个位置的距离。在COMSOL里用“一维绘图组”和“截线”数据集在焦点处画一条横向截线再把表达式调成20*log10(abs(acpr.p)/max)读取−6dB对应的横坐标差。横向分辨率就是这个值。第三聚焦增益。把相控阵聚焦时焦点的峰值声压除以同等条件下单个阵元在焦点位置的声压得到的就是聚焦增益通常用dB表示。这个指标直接反映相控阵相比单阵元的增强能力工程上很关键。4.3 参数化扫描结果怎么“看”参数化扫描后“一维绘图组”里选择不同参数组合的数据集可以叠加多个焦点压力曲线。我通常把频率或阵元间距作为横轴画出焦点声压随参数变化曲线用来判断最优工作点。举个例子扫描阵元间距时如果发现焦点声压在某个间距之后明显衰减同时旁瓣开始升高说明可能出现栅瓣效应这个间距就该避开。反之如果焦点声压随阵元数线性上升说明聚焦在正常增长阶段。也可以把不同参数下的声场图横竖排列在一起做一个“报告”用COMSOL自带的“报告工具”生成带图带的PDF沟通方案、写专利都够用。4.4 经验公式焦点尺寸可以提前算个底做参数扫描之前用近似公式估算一下预期的焦点尺寸能帮你判断仿真结果靠不靠谱。空间分辨率的近似公式是焦点横向宽度 ≈ λ * F / D其中λ是波长F是焦距D是阵列总孔径宽度。比如f500 kHz水中声速1480 m/sλ≈2.96mm阵元数16间距0.6mm孔径D9.6mm焦距F30mm。代入得到焦点宽度约9.25mm。如果仿真结果和这个数量级对不上就该回头检查模型了。另一方面阵元间距有一个硬限制为了避免栅瓣干扰阵元间距一般要求小于等于半个波长。上面例子里半个波长是1.48mmpitch取0.6mm是安全的如果扫描到1.5mm以上图像上激光瓣就会冒出来。5. 实操踩坑记录与排错技巧5.1 边界反射带来的假干涉条纹这是新手上路最容易撞的墙。现象是声场图中不仅焦点处有亮斑域边缘还出现一圈圈规则的弧形条纹越靠近边界越明显。原因就是截断边界没有加PML或者PML厚度不够、网格不匹配。解决办法很简单给计算域外层包一圈厚度至少一个波长的PML区域确保PML和内部物理域的网格交界处节点连续。如果你用了映射网格注意PML区域的网格不要跨层扭曲。曾经有人图省事把PML厚度设成0.5个波长结果边界反射还是明显加厚到1.5个波长后干净多了。5.2 内存不足导致求解失败COMSOL求解器报“内存不足”在三维模型中很常见。我的排查顺序是优先利用对称性能切一半就切一半。检查网格是否局部过密比如阵元边缘处是否出现了极小尺寸的网格可以用网格统计里的“最小网格质量”来看。求解器上二维模型直接用MUMPS或者PARDISO直接求解器最稳三维大模型才考虑迭代求解器加预处理器。迭代求解器对声学问题收敛不一定快实际经验是三维模型中用“GMRES多网格”有时能压住内存但参数不合适会发散。如果你不是特别熟还是老实靠网格控制和对称性降低规模。5.3 单位换算与复声压的显示误区频域计算得到的声压是一个复数很多新手画出来图后觉得声场不对称其实只是显示了实部。把表达式改成abs(acpr.p)之后瞬间就正常了。另一个单位坑是分贝。COMSOL里“声压级”默认参考声压是20微帕但在水声和超声领域参考声压标准并不统一。如果要自定义基准需要手动写表达式20*log10(abs(acpr.p)/p_ref)p_ref按你自己的标准设置别直接套默认结果。5.4 CAD导入报错“不支持的拓扑”这个报错通常出现在使用外部CAD几何导入COMSOL时尤其是STEP格式文件含有微小倒角、重叠面、破面等。解决方法有几个在CAD软件里简化模型去掉对声场影响可以忽略的小特征。在COMSOL“几何”里利用“修复”工具自动删除短边、合并薄片。最省心的是直接用COMSOL内建几何图元重建声学域毕竟仿真域就是一个水域加阵列边界没必要非得从复杂CAD模型导进来。5.5 模型验证的“金标准”参数化模型跑出来再漂亮也得有一个正确的参照物来验证。我常用的验证方法是把阵列激励设定为所有阵元同相位计算远场指向性和理论公式对比。单阵元、线阵远场指向性有解析解误差在1dB以内说明模型基础设置没问题。再进阶一点就是用4.4节的经验公式核对焦点尺寸和焦点增益的数量级。焦点宽度、旁瓣水平这些指标和理论趋势一致模型就可信。如果仿真结果和理论差了一倍以上别急着调网格先回头检查介质声速、频率单位和延时符号。5.6 那些看起来很吓人的热词其实与本模型无关搜索COMSOL相关资料时总会看到“塑性变形”“欧拉角”“摩擦角”“粘度随温度变化”这些词。这些大部分是结构力学、流体力学或者电磁场的仿真议题和超声相控阵聚焦声场模型没有直接关系。你如果在声学域里面设置了塑性参数或者欧拉角变量反而会让模型变得混乱。我个人的经验是仿真之前先想清楚这一版要回答什么问题。超声相控阵聚焦核心就是频域声场、聚焦法则、参数扫描这三个点。把边界反射挡在外面把延时法则算对把网格按波长划好这个模型就能踏踏实实给设计提供依据。之后就算要扩展比如加入压电耦合、时域脉冲分析、热效应分析也是在现有框架上加模块而不是返工重来。
COMSOL超声相控阵聚焦仿真：频域建模与参数化实现全解析
做超声检测和声学设计的朋友一定经历过这种尴尬新的相控阵探头还在论证阶段阵元数量、间距、频率都悬着打样试错的钱花得心疼实验台又排不上队。这时候COMSOL有限元超声相控阵聚焦仿真模型就派上用场了——把阵元参数全部参数化在频域里一次算出稳态声场顺手把焦点位置、焦点尺寸、聚焦增益这些指标拉出来整个方案可行性判断一两天就能落地。这篇文章不是把COMSOL官网案例抄一遍而是把实际建模过程中那些真正影响结果的细节、命令和坑整理出来。内容包括相控阵聚焦的频域求解逻辑、模型搭建里物理场和边界的设置思路、延时法则在COMSOL里的工程化实现、参数化扫描和后处理分析方法最后把常见的边界反射、内存爆炸、CAD拓扑报错这些问题一并排掉。适合正在做超声聚焦、无损检测换能器设计、水声工程仿真的工程师参考哪怕你刚接触COMSOL跟着这套思路也能跑通一套能用的聚焦仿真模型。1. 相控阵聚焦仿真的底层逻辑1.1 相控阵到底在做什么相控阵换能器的核心是一排独立的压电阵元每个阵元可以单独控制激励信号的延时——也就是相位。生活中可以这样理解一排人喊口号如果所有人同时张嘴声音是朝四面八方散的如果让最远的人先喊、最近的人后喊这个“错峰喊话”能让声音在一个预先指定的点上汇聚恰好像聚光灯。超声相控阵聚焦就是这个原理的工程版。设焦点坐标为F(x0, y0)第i个阵元的中心坐标为(xi, yi)声速为c。声波从第i个阵元传到焦点需要的时间是ti sqrt((xi - x0)^2 (yi - y0)^2) / c要让所有阵元发出的波同时到达焦点激励信号就需要按时间差做补偿。通常以最远阵元为基准取Δt_i t_max - ti也就是距离焦点越近的阵元越晚激励。把Δt_i换算成频域相位偏移就得到每个阵元的激励相位。这一整套延时分布业内叫“聚焦法则”是相控阵聚焦的灵魂。1.2 频域仿真为什么是工程主力COMSOL里做声学仿真最常用的有两类研究时域研究和频域研究。时域能捕捉脉冲波传播的完整瞬态过程看得到波包如何运动、如何反射物理画面很直观。但时域求解有严格的时间步长约束网格越细、域越大步数越多一个二维模型跑几个小时稀松平常三维模型更是灾难。频域仿真假设系统处于稳态正弦激励下求解的是线性方程组不存在时间步进问题。相控阵聚焦的典型工作状态是连续波或窄带脉冲激励频率基本固定在中心频率附近用频域求解在物理上完全合理。更重要的是频域天然适合做参数化扫描频率、阵元间距、焦点深度、阵元数量这些东西一变重新解一次线性系统就行不需要从头跑时间历程。工程选型时我的建议很明确如果只看聚焦效果、旁瓣水平、焦点尺寸这类稳态指标直接上频域如果非要研究脉冲波形畸变、非线性效应或者瞬态聚焦过程再考虑时域。1.3 有限元方法在这个场景里的不可替代性超声换能器声场也可以用解析方法算最经典的是瑞利-索末菲积分通过把每个阵元等效为点源叠加得到空间任一点的声压。解析方法的优点是快几十秒就能出来一幅声场图缺点是边界条件简化得厉害碰到复杂几何形状、近场区域、非均匀介质、多层结构就不那么好用。有限元方法把计算域离散成网格单元在每个单元内用形函数逼近声压场能处理任意几何形状和材料分布近场精度比解析法高得多。代价是网格要划得足够密矩阵求解要吃掉不少内存。现实中做相控阵聚焦仿真通常就是几何还算规整、但阵元排布和边界情况比较复杂有限元这种“一刀切”的通用求解器反而最省心。2. 模型搭建几何、物理场与边界设定2.1 几何建模2D还是3D对称性该用就用我一开始直接烧三维模型结果网格数量大得离谱一台32G内存的工作站跑一个8阵元阵列就到极限了。说到底相控阵聚焦仿真不一定非得3D。如果阵列是一维排布也就是阵元沿x轴排列、声场在xy平面内传播那用2D模型完全够用。COMSOL里新建“二维”组件计算域设成一个矩形水域阵列边界放在上边沿或者左边沿焦点落在域内部。2D模型网格量小两个数量级参数扫描时可以一口气跑几十组。二维模型还能进一步利用对称性。阵列排列和焦点都关于某条轴线对称时只建一半模型在对称轴上施加对称边界条件计算量立刻减半甚至减更多。COMSOL的“对称”边界条件在压力声学里对应的是硬边界条件也就是法向振速为零物理上等同镜像反射。三维模型也不是不能做适合阵列是面阵、焦点在三维空间里、需要考察横向和纵向聚焦分辨率的场景。实际经验是先用2D把方案跑通把参数趋势摸清再上3D做最终验证这条路径效率最高。CAD模型的导入偶尔会遇到“转换为CAD内核时不支持的拓扑”这个报错通常是STEP文件里存在重叠面、微小边或者退化几何。我的习惯是直接用COMSOL内置的几何工具重建声学域——反正就是一个矩形或圆柱水域阵列边界用矩形阵列操作生成几十秒就能搞定比在CAD软件里反复修补导出再导入省心得多。2.2 物理场接口压力声学、频域是默认选项COMSOL里做常规声学仿真物理场接口选择“声学模块 压力声学 压力声学频域”。这个接口求解的是亥姆霍兹方程也就是频域波动方程变量是声压p。在均匀介质中方程形式是∇·(−1/ρ0 ∇p) − ω^2/(ρ0 c^2) p 0其中ρ0是介质密度c是声速ω2πf是角频率。这个接口会求解复声压场实部代表瞬时声压在某一相位时刻的值虚部加入了相位信息。声学域的材料设置为水密度998 kg/m³声速1480 m/s。如果目标介质是其他液体或固体替换材料参数即可。这里要注意材料参数里的声速必须与你设置的频率、波长配套否则网格尺寸和PML厚度即使看起来合理实际却可能差出一个数量级。换能器的激励边界有几种处理方式。最简化的方法是把每个阵元边界设成“法向加速度”边界条件或者直接给阵元边界赋一个复声压边界条件。如果你不关心压电陶瓷内部电场和机械振动的耦合这就是最快的路子。如果要做完整的压电超声换能器仿真那必须加上“压电设备”物理场接口把PZT材料的压电耦合矩阵、电极边界全建出来网格和计算量会相应增大。2.3 边界条件别让声波“碰壁”无限大水域在有限元模型里只能截取有限区域截断边界如果处理不好声波会在边界上反射回来聚焦声场里就会混入一圈圈的假干涉条纹。这个干扰早期特别容易忽略因为图像上看起来像漂亮的波纹实际上全是假信号。COMSOL里处理开放声场的最标准做法是加完美匹配层也就是PML。在声学域的外围再包上一圈环状区域把这个区域设为“完美匹配层”的域条件。PML的原理是在内部对方程做坐标拉伸变换使进入该区域的波被指数衰减吸收理论上不产生反射。PML参数设置有两个关键点。第一是厚度建议至少一个波长稳妥一点取1.5到2个波长。第二是网格PML区域内网格应该尽量规整通常用映射网格或扫掠网格生成规则矩形网格这样可以避免不规则网格引起的局部数值误差。还有一个容易被忽略的点就是PML外面不需要再设任何边界条件程序默认处理为吸收。物理域和PML交界处的网格连续性要保证否则出现数值反射焦点位置的声压分布可能已经受了污染。2.4 网格尺寸怎么定频域声学仿真的网格规则可以用一句话概括每个波长内至少5到6个二阶单元。声速c、频率f和波长λ的关系是λc/f。以水声速1480 m/s为例不同频率对应的波长和推荐网格尺寸如下频率波长推荐网格每波长6单元250 kHz5.92 mm约1.0 mm500 kHz2.96 mm约0.5 mm1 MHz1.48 mm约0.25 mm在COMSOL网格节点里可以直接选择“声学”预定义网格大小软件会根据物理场接口自动估算波长并给出推荐尺寸。我自己的流程是从较粗网格开始试算逐步细化做一次标准的网格收敛性检查。如果两次加密后焦点声压变化小于1%网格量就是合适的。PML区域网格要单独处理。我的做法是在PML区域用映射网格划定结构化网格层内沿厚度方向划分4到6层即可太多太密只会徒增计算量。3. 聚焦延时法则在COMSOL里的工程化实现3.1 延时计算先算传播时间再转相位聚焦法则的工程实现在COMSOL里有几种套路先说最朴素的。假设有N个阵元沿x轴排列阵元间距为pitch第i个阵元的中心坐标为xi。焦点坐标是(x_f, y_f)。第i个阵元到焦点的距离为ri sqrt((xi - x_f)^2 y_f^2)对应的传播时间为ti ri/c。取所有阵元中传播时间最大值t_max第i个阵元的延时为Δt_i t_max − ti。在频域谐波分析中时间延迟Δτ对应的相位偏移是−2πfΔτ符号约定取决于COMSOL的时谐因子约定工程计算中你只需要保证各阵元间的相对相位差是对的。这样第i个阵元的激励就可以写成复值边界条件p_i p0 * exp(−j * 2 * pi * f * Δt_i)p0是单个阵元的激励声压幅值设置为1 Pa便于后处理归一化。3.2 用分段函数把延时映射到阵列边界每个阵元的坐标不同延时也逐元变化。在COMSOL里不需要为每个阵元单独建边界条件更聪明的办法是用一个分段函数把延时映射到空间坐标上。在“全局定义 函数”中新建一个“分段函数”变量设为x区间按阵元边界划分。比如阵元1的x范围是0到0.6mm此区间内函数值是Δt_1阵元2的x范围是0.6到1.2mm函数值是Δt_2依此类推。然后把换能器整个阵列边界设成一个统一的压力边界条件赋值为p0 * exp(−j * 2 * pi * f * delay_func(x))由于分段函数只在相应区间内取相应延时值就实现了各阵元独立相位控制。这个思路最大的好处是参数扫描时只要更新延时函数整个模型自动响应不需要动几何结构。阵元数量很多时也可以用外部CSV表格导入延时数据再通过“插值函数”映射到坐标上。听说有些项目里还有把MATLAB算好的聚焦法则直接灌进COMSOL的做法本质上都是一样。3.3 参数化扫描让所有关键参数活起来COMSOL的“全局参数”是参数化仿真的心脏。我把f0工作频率、pitch阵元间距、elem_count阵元数、focus_x和focus_y焦点坐标全部定义成全局参数。几何构建时阵列位置用参数驱动。二维建模时先用一个“矩形”表示单个阵元边界再用“阵列”操作x方向间距设为pitch个数设为elem_count。这样pitch和elem_count一变几何自动重建网格也跟着自适应。研究设置中在“研究1”里选择“频域”把频率设为f0。然后在研究步骤前添加“参数化扫描”把需要扫描的参数加进去比如扫f0从300kHz到700kHz步长50kHz同时pitch从0.4mm扫到1.0mm步长0.2mm。COMSOL会自动计算所有参数组合每个组合解一次频域方程不需要你手动干预。这里给一个参数表参考方便你直接抄作业参数名含义典型值扫描范围f0工作频率500 kHz250k–1 MHzpitch阵元间距0.6 mm0.3–1.2 mmelem_count阵元数量168–32focus_x焦点横向坐标0 mm−10–10 mmfocus_y焦点纵向深度30 mm10–60 mmp0阵元激励幅值1 Pa固定1参数化扫描跑完后所有结果都存放在不同数据集里。后处理时可以通过“数据集”下拉列表切换不同参数组合的结果非常方便。4. 聚焦效果分析怎么证明“聚焦成功”了4.1 第一眼判定声压幅值分布图最直观的分析是在“二维绘图组”里画声压幅值分布。这里的第一个坑是变量选择COMSOL频域求解得到的是复声压acpr.p直接用实部或虚部绘图会看到明显的不对称图案好像聚焦非常怪异。必须用幅值也就是abs(acpr.p)或者声压级Lp才能得到真实的声场包络。画出来后判断聚焦是否成功的三个视觉标准焦点位置是否有明显的亮斑也就是声压幅值极大值。焦点周围是否存在过大的旁瓣主瓣和旁瓣的亮度差是否清晰。声束是否沿着预期的轴线传播有没有明显偏斜。如果亮斑不在预设焦点位置先检查延时法则的符号和坐标定义是不是反了这是我之前踩过最多次的坑。4.2 量化指标这样算视觉判断只能用来粗筛写报告还需要量化指标。COMSOL的“派生值”菜单里可以直接计算这些第一焦点位置。使用“派生值 体最大值”或“面最大值”在“表达式”里填abs(acpr.p)让软件帮你定位声压幅值最大的坐标。如果这个坐标和预设焦点坐标偏差在可接受范围内说明聚焦法则设置正确。第二焦点尺寸。一般是提取通过焦点的直线上声压分布然后找主瓣下降6dB两个位置的距离。在COMSOL里用“一维绘图组”和“截线”数据集在焦点处画一条横向截线再把表达式调成20*log10(abs(acpr.p)/max)读取−6dB对应的横坐标差。横向分辨率就是这个值。第三聚焦增益。把相控阵聚焦时焦点的峰值声压除以同等条件下单个阵元在焦点位置的声压得到的就是聚焦增益通常用dB表示。这个指标直接反映相控阵相比单阵元的增强能力工程上很关键。4.3 参数化扫描结果怎么“看”参数化扫描后“一维绘图组”里选择不同参数组合的数据集可以叠加多个焦点压力曲线。我通常把频率或阵元间距作为横轴画出焦点声压随参数变化曲线用来判断最优工作点。举个例子扫描阵元间距时如果发现焦点声压在某个间距之后明显衰减同时旁瓣开始升高说明可能出现栅瓣效应这个间距就该避开。反之如果焦点声压随阵元数线性上升说明聚焦在正常增长阶段。也可以把不同参数下的声场图横竖排列在一起做一个“报告”用COMSOL自带的“报告工具”生成带图带的PDF沟通方案、写专利都够用。4.4 经验公式焦点尺寸可以提前算个底做参数扫描之前用近似公式估算一下预期的焦点尺寸能帮你判断仿真结果靠不靠谱。空间分辨率的近似公式是焦点横向宽度 ≈ λ * F / D其中λ是波长F是焦距D是阵列总孔径宽度。比如f500 kHz水中声速1480 m/sλ≈2.96mm阵元数16间距0.6mm孔径D9.6mm焦距F30mm。代入得到焦点宽度约9.25mm。如果仿真结果和这个数量级对不上就该回头检查模型了。另一方面阵元间距有一个硬限制为了避免栅瓣干扰阵元间距一般要求小于等于半个波长。上面例子里半个波长是1.48mmpitch取0.6mm是安全的如果扫描到1.5mm以上图像上激光瓣就会冒出来。5. 实操踩坑记录与排错技巧5.1 边界反射带来的假干涉条纹这是新手上路最容易撞的墙。现象是声场图中不仅焦点处有亮斑域边缘还出现一圈圈规则的弧形条纹越靠近边界越明显。原因就是截断边界没有加PML或者PML厚度不够、网格不匹配。解决办法很简单给计算域外层包一圈厚度至少一个波长的PML区域确保PML和内部物理域的网格交界处节点连续。如果你用了映射网格注意PML区域的网格不要跨层扭曲。曾经有人图省事把PML厚度设成0.5个波长结果边界反射还是明显加厚到1.5个波长后干净多了。5.2 内存不足导致求解失败COMSOL求解器报“内存不足”在三维模型中很常见。我的排查顺序是优先利用对称性能切一半就切一半。检查网格是否局部过密比如阵元边缘处是否出现了极小尺寸的网格可以用网格统计里的“最小网格质量”来看。求解器上二维模型直接用MUMPS或者PARDISO直接求解器最稳三维大模型才考虑迭代求解器加预处理器。迭代求解器对声学问题收敛不一定快实际经验是三维模型中用“GMRES多网格”有时能压住内存但参数不合适会发散。如果你不是特别熟还是老实靠网格控制和对称性降低规模。5.3 单位换算与复声压的显示误区频域计算得到的声压是一个复数很多新手画出来图后觉得声场不对称其实只是显示了实部。把表达式改成abs(acpr.p)之后瞬间就正常了。另一个单位坑是分贝。COMSOL里“声压级”默认参考声压是20微帕但在水声和超声领域参考声压标准并不统一。如果要自定义基准需要手动写表达式20*log10(abs(acpr.p)/p_ref)p_ref按你自己的标准设置别直接套默认结果。5.4 CAD导入报错“不支持的拓扑”这个报错通常出现在使用外部CAD几何导入COMSOL时尤其是STEP格式文件含有微小倒角、重叠面、破面等。解决方法有几个在CAD软件里简化模型去掉对声场影响可以忽略的小特征。在COMSOL“几何”里利用“修复”工具自动删除短边、合并薄片。最省心的是直接用COMSOL内建几何图元重建声学域毕竟仿真域就是一个水域加阵列边界没必要非得从复杂CAD模型导进来。5.5 模型验证的“金标准”参数化模型跑出来再漂亮也得有一个正确的参照物来验证。我常用的验证方法是把阵列激励设定为所有阵元同相位计算远场指向性和理论公式对比。单阵元、线阵远场指向性有解析解误差在1dB以内说明模型基础设置没问题。再进阶一点就是用4.4节的经验公式核对焦点尺寸和焦点增益的数量级。焦点宽度、旁瓣水平这些指标和理论趋势一致模型就可信。如果仿真结果和理论差了一倍以上别急着调网格先回头检查介质声速、频率单位和延时符号。5.6 那些看起来很吓人的热词其实与本模型无关搜索COMSOL相关资料时总会看到“塑性变形”“欧拉角”“摩擦角”“粘度随温度变化”这些词。这些大部分是结构力学、流体力学或者电磁场的仿真议题和超声相控阵聚焦声场模型没有直接关系。你如果在声学域里面设置了塑性参数或者欧拉角变量反而会让模型变得混乱。我个人的经验是仿真之前先想清楚这一版要回答什么问题。超声相控阵聚焦核心就是频域声场、聚焦法则、参数扫描这三个点。把边界反射挡在外面把延时法则算对把网格按波长划好这个模型就能踏踏实实给设计提供依据。之后就算要扩展比如加入压电耦合、时域脉冲分析、热效应分析也是在现有框架上加模块而不是返工重来。
