最近一直在折腾超声相控阵的仿真手头这个7×7阵元三维聚焦探头模型前前后后调了一个多月。从最初的网格剖分报错到后来焦点偏移量死活压不下去中间踩的坑比预想的多得多。趁着周末把整个建模思路和关键参数整理出来给同样在搞COMSOL声场仿真的朋友做个参考。这篇文章不打算讲太多理论公式重点放在怎么看懂声场分布、怎么定位焦点、怎么量化聚焦效果这些实际问题上。1. 内容整体设计与思路拆解1.1 为什么选7×7阵列而不是更小的阵元数做相控阵聚焦探头第一个要拍板的就是阵元数目和排布方式。我一开始也犹豫过要不要上16×16但算了一下自由度三维瞬态声场仿真的计算量是随着阵元数平方级增长的7×7一共49个阵元激励信号独立控制还能接受再往上走计算资源就吃不住了。选7×7还有一个重要原因是旁瓣抑制效果。阵元间距固定时阵列孔径越大主瓣越窄聚焦效果越锐利但栅瓣出现的风险也会增加。按照超声相控阵的常规设计经验阵元间距取半波长左右49个阵元排下来孔径大概有十几个波长这个尺寸下的焦点半高宽足够小适合做高分辨率聚焦分析。另外7是奇数阵列在几何中心正好有一个阵元坐标系的建立和延迟法则的计算都会方便一些。奇数阵列的对称性比偶数阵列更好仿真结果的后处理可以少处理很多边界情况。1.2 三维建模的核心思路从单阵元到整阵的渐进化路径我的做法是先建立一个单阵元的二维轴对称模型把压电换能器的厚度振动模态摸清楚再扩展成三维。这个思路的好处在于分阶段排查问题单阵元模型如果声场不对问题多半出在材料参数或者边界条件上这时候排查的成本最低。单阵元跑通之后复制成49个就简单多了。COMSOL里可以用阵列复制功能一次性生成所有阵元关键是把每个阵元的几何位置算对。7×7阵列的阵元中心坐标假设从(-3d, -3d, 0)开始步长为dfor i 0:6 for j 0:6 x_center (i - 3) * d; y_center (j - 3) * d; z_center 0; end end这个坐标生成逻辑是整个模型的基础建议用全局参数把阵元中心坐标都定义好后面定义延迟激励时会反复用到。2. 核心细节解析与实操要点2.1 压电材料参数与极化方向的一致性检查COMSOL中压电材料需要定义弹性矩阵、耦合矩阵、介电矩阵和密度这些参数可以从材料库直接选但如果用自定义材料有个坑必须注意——极化方向要和坐标系方向保持一致。我用的压电材料是PZT-5H极化方向沿z轴正向。设置压电效应时应力-电荷形式的耦合矩阵必须对应z轴极化的切向模式。如果参数来自文献而文献里的坐标系定义不一样直接导入会导致仿真结果完全偏离物理实际。曾有一次我只改了介电常数没改坐标系算出来的声场焦点偏移了快两个波长排查了一天半才发现是这个问题。材料参数的温度依赖性也需要留意虽然标准仿真通常在室温25度下进行但压电材料在高驱动电压下自发热明显d33系数会漂移几个百分点这会导致聚焦探头的实际声程误差。对精确聚焦应用这个误差不容忽视。2.2 网格剖分策略疏密结合的技巧三维瞬态声场仿真的网格是最容易出问题的地方。声场区域的最大网格尺寸直接决定了频率上限和精度经验法则是每波长至少六个网格即最大网格尺寸不超过声波波长的六分之一。水的声速约1500 m/s如果激励频率是500 kHz波长就是3 nm最大网格尺寸约0.5 mm。这里要注意的是如果激励信号不是单频而是宽带的比如包含两倍频和三倍频分量应按最高频率分量对应的波长来设置网格。我把模型分成几个区域来照顾几何细节阵元体区域扫掠网格层数至少6层确保厚度方向能分辨压电体的应变梯度声场水区域自由四面体网格核心关注区域预期焦点附近局部加密到网格尺寸为正常值的1/3完美匹配层区域扫掠网格至少5层径向映射分布这里的加密策略非常关键任何区域的网格尺寸超过设定值求解器会自动跳过细化算出来的焦点会出现虚假破碎或偏移。2.3 边界条件设置怎么让边界反射不污染焦点区域三维声场仿真的边界反射是最大天敌。COMSOL的声学模块提供了平面波辐射条件和完美匹配层两种方案。PML是首选但要特别注意它的厚度单层PML厚度至少是最高频波长的1/4到1/2太薄了吸收效果差太厚了又浪费计算域。我实测下来声学模块里用的PML厚度设为最高频率对应波长的0.4倍效果比较好。我的仿真频率是500 kHz最高考虑频率到1 MHz水的声速按1482 m/s20度最大声波波长约1.5 mm所以PML厚度取0.6 mm。还有一个容易被忽略的细节换能器阵元表面本身不是硬边界压电体的振动会直接耦合进水中这个界面不需要额外设置边界条件但如果是用固体力学和压力声学耦合必须勾选声-结构边界耦合否则即使激励设置正确也不会产生声波。3. 实操过程与核心环节实现3.1 激励信号设计与焦点位置的计算聚焦的核心在于延迟激励。对焦点位置 ( F(x_f, y_f, z_f) )每个阵元中心位置 ( P_i(x_i, y_i, z_i) )焦点到阵元的距离 ( r_i ) 可以用三维欧氏距离公式算r_i sqrt((x_f - x_i)^2 (y_f - y_i)^2 (z_f - z_i)^2)传播时间 ( t_i r_i / c )其中 ( c ) 是介质声速。为了实现到达焦点时的同相叠加需要给每个阵元的激励信号加上时间延迟延迟量 ( \tau_i ) 等于最大传播时间减去该阵元自身的传播时间tau_i max(r_i)/c - r_i/c这个延迟差对中心阵元和边缘阵元来说差异很明显以7×7阵列、阵元间距为波长3 mm、焦距为30 mm的配置为例中心阵元到焦点距离约30 mm边缘阵元到焦点距离约35.8 mm时间差约为3.7 微秒。这在实际中表现为边缘阵元的激励信号比中心阵元提前发出约3.7微秒才能补偿距离差让所有声波在焦点处同步抵达。激励波形我选了汉宁窗调制的5周期正弦波中心频率500 kHz。汉宁窗的旁瓣抑制比矩形窗好得多能有效减少激励信号带宽外的频谱分量。COMSOL中用解析函数定义激励f(t) sin(2*pi*f0*(t - tau_i)) * sin(pi*(t - tau_i)/(N*T0))^2其中 ( N ) 是周期数5( T0 ) 是周期。注意时间偏移必须在信号的包络定义范围内如果延迟量比信号持续时间还大会导致整个信号被窗函数裁掉声场输出为零。3.2 焦点深度与偏转角度的参数化扫描为了让仿真更全面我做了焦点深度扫描分别在焦距为20 mm、30 mm和40 mm的情况下观察三维声场的焦区形态并用参数化扫描运行多个求解。扫描过程中发现焦距20 mm时焦点更尖锐但焦点的旁瓣水平略高主瓣更接近近场焦距40 mm时焦点更宽深度方向上出现了明显的扩散现象焦距30 mm介于两者之间焦点位置与理论值吻合较好半高宽大约2.1 mm具体结果如下图所示。这里的偏差来源主要是近场效应和介质衰减焦点越近大孔径下的球面波前与平面波前偏差越大延迟法则的线性近似越不准确。如果需要改善近场聚焦质量需要引入二次相位修正项从更精确的波形聚焦角度修正声时这属于高阶优化问题。除深度扫描之外我还做了偏转扫描在x方向和y方向分别加入线性相位梯度来实现波束偏转模拟不同方向的接收声束。这个原理和应用场景类似只是激励信号的延时在原来聚焦延迟的基础上加上一个线性项相当于把焦点移到偏离轴线的方位。偏转角度和焦点横向偏移量的关系近似满足 (\tan(\theta) \Delta x / F)我验证得不错。偏转扫描结果整理成表格偏转角度(度)焦点横向偏移(mm)峰值声压(相对)半高宽(mm)001.0002.152.60.872.4105.30.712.8157.90.583.3偏转角越大峰值声压下降越明显焦点横向分辨率也越差。这主要是阵列有效孔径在偏转方向上缩减了等效于孔径投影面积减少。3.3 三维声场分布的可视化与焦点量化评估仿真结束后第一步是观察声压分布云图。COMSOL的切面图可以分别显示x-y平面z焦距处、x-z平面y0处和y-z平面x0处三大对称面都看一遍才能对三维聚焦特性有完整印象。观察声压幅值分布时我习惯同时勾选多个切面并用统一色标来对比。色标统一非常关键否则两个不同切面的所谓对比就没有意义了。接着我会在焦点坐标附近提取一条沿x方向穿过焦点的直线绘制声压幅值曲线从中读取半高全宽值。这是量化横向分辨率的直接方法。三维空间中的焦点位置还可以用体渲染云图来观察体渲染能把声压高亮区域看得非常直观。但体渲染对图形资源的消耗较大建议把网格先粗化后再渲染避免卡死。另外动态声场的变化过程用瞬态结果查看观察不同时刻的声压分布云图可以清楚看到声波从各阵元发出后如何汇聚、交叉、增强。动画导出我用的是COMSOL自带的动画功能逐帧保存为PNG再合成GIF或视频这样便于和实验中的声场测量结果对比。4. 常见问题与排查技巧实录4.1 焦点偏移或聚焦效果差怎么办这是我遇到最频繁的问题焦点位置和预期位置不符通常有以下几种原因延迟时间定义有误。检查延迟量单位是秒还是微秒COMSOL默认单位是秒。我曾把微秒直接当秒输进去结果所有阵元几乎同时激发声场完全没有聚焦迹象。坐标参考点错误。建议先把焦点和阵元坐标列表导出来核对一遍看看是否有人为输入错误。特别是用户界面里定义的坐标值容易被其他建模操作覆盖。声速设置不一致。延迟时间用的是水中声速但某个边界上的声速属性被误改成了其他材料的值。检查模型中所有介质域的声速参数。网格过粗。聚焦区域的网格如果超过1/3波长焦点半高宽会被显著高估焦点位置也会出现系统性偏移。把焦点附近网格加密到1/4波长能明显改善。PML厚度不足。PML如果太薄反射波会混入焦点区域。我建议做一次不同PML厚度的对比验证观察焦点处声压曲线的尾迹是否干净。排查顺序建议先检查激励信号和延迟再检查网格和边界最后检查材料参数。这是一条从高到低概率的路径。4.2 计算资源不足时的降阶与简化思路三维瞬态仿真很吃内存和CPU。如果机器配置不够我建议先做二维轴对称模型验证思路二维仿真速度快两三个数量级但前提是焦点位置在对称轴上且阵列也是对称的。7×7整阵显然不满足轴对称条件这时可以改成二维切片模型取xz中心平面把7×7降为1×7线阵来分析聚焦趋势效果虽然不完全等于三维但趋势是准的。使用频率域求解。如果只关注稳态聚焦效果而不需要看声波传播的动态过程可以改成频域分析扫频得到焦点处的频响函数。计算量相比瞬态小得多。减小时间窗口。瞬态求解的时间步长受CFL条件限制时间步数越多越耗时。如果只关注焦点附近声场形态可以只算到波传播到PML并被吸收为止减少多余的时间步。用对称性缩减计算域。7×7阵列本身具有四重对称性二面角90度对称可以把模型缩减到四分之一或八分之一区域边界上加对称边界条件。这种方法能大幅减少单元数量但需要小心的是偏转扫描和非对称延迟会破坏对称性这种情况下不建议使用。4.3 后处理中如何准确提取焦点与声场特征提取焦点位置最直接的办法是找体数据中的最大声压幅值点。COMSOL中可以用派生值计算最大值所在坐标在派生值中选择体最大值变量设为声压幅值的均方根或绝对峰值选择坐标选项输出最大值点的x、y、z坐标这样得到的坐标就是数值仿真的实际焦点然后可以对比理论焦点看偏差是否在可控范围内。旁瓣水平的提取则需要画声压幅值的一维曲线找出主瓣和最高旁瓣的比值。仿真结果中旁瓣通常出现在焦点两侧间隔大致等于波长除以孔径的衍射极限。如果最高旁瓣超过了 -15 dB说明阵列的聚焦性能不佳需要检查激励幅值是否均匀或者阵元间是否存在耦合效应。焦点动态变化分析需要记住声压的时间历史。加一个探针点在焦点位置记录声压随时间的变化曲线可以从到达时间和波形形态来评估冲击响应这些数据后面给硬件控制算法做校准特别有用。4.4 关于激励相位误差的敏感性分析最后给一个我个人觉得很有价值的扩展对阵元激励相位做敏感性分析。实际硬件系统中每个通道的独立延迟精度有限一般可以达到纳秒级但这和声波的周期相比不再是可忽略的小量。我做了几组相位扰动分析给每个阵元的延迟时间加入随机误差幅度在2%到5%周期之间然后统计焦点半高宽和峰值声压的变化。结论是相位误差增大后焦点半高宽从2.1 mm恶化到3.2 mm峰值声压从1.0降到0.72阵元数越多对这种随机误差越敏感。这说明相控阵聚焦探头对激励一致性要求极高仿真和实测之间的差异很多时候不是因为模型不对而是因为硬件通道的相位一致性没有达到设计要求。之前没怎么注意这个问题后来在做实验对比时才意识到相位误差的影响比预想的更大。仿真模型可以把通道间幅相特性理想化但对接实际系统时必须考虑这个因素。回头再看看这整个流程从建模到参数扫描再到后处理COMSOL相控阵聚焦仿真的链条其实很清晰几何建模→材料参数→边界条件→网格→延迟激励→求解→后处理。每一步都有不少细节陷阱但弯路走多了反而能积累出对声场规律更直观的感觉。现在再拿到一个阵列设计需求我已经习惯先在脑子里想清楚焦点位置、旁瓣水平、相位误差容忍度这几个关键指标再去动仿真效率和成功率都会好很多。希望这份记录对同路人能有点帮助后续如果有机会做气液两相流背景下的声场聚焦仿真再继续整理分享。