COOMSOL超声仿真——复合材料空气耦合超声单侧检测仿真研究（论文仿真） 模型介绍：本模型采用压力声学、固体力学物理场，通过给定一个声源，描述声音在复合材料中的传播。
搞超声仿真最头疼的就是多物理场耦合，尤其是复合材料这种各向异性的主儿。今儿拿COMSOL折腾空气耦合超声单侧检测模型，上来就得面对声波在固体和空气中反复横跳的物理现象。
先说说模型架构：左边是压电陶瓷片发射超声波，右边是8层碳纤维板，最右端接了个空气层当接收介质。重点在于让声波从固体传入空气还能被检测到——这玩意儿在实际检测中就是个玄学问题，毕竟99%的声波能量在固-气界面就被反射了。
建模时先得处理几何结构。复合材料层用参数化堆叠，代码里这么搞：
for i = 1:8 geometry.create('layer'+i, 'Block'); geometry('layer'+i).set('size', [0.1 0.1 0.002]); geometry('layer'+i).set('pos', [0 0 0.002*(i-1)]); end每层0.2毫米厚度，手动循环堆八层比GUI里拖拽快得多。这里有个坑：各向异性材料需要单独定义材料坐标系，不然纤维方向对不上实际工况。
物理场设置才是重头戏。压力声学模块和固体力学模块必须勾选多物理场耦合里的“声-结构边界”。关键参数在代码里长这样：
% 空气域边界条件 model.physics('ac').feature('abc1').set('V', {'0' '0' '0'}); % 硬声场边界 model.physics('ac').feature('pr1').set('p0', '1e-3*exp(-((t-1e-6)/0.5e-6)^2)'); % 高斯脉冲声源 % 固体域材料参数 model.material('comp').propertyGroup('def').set('youngs_modulus', {'30e9' '8e9' '8e9'}, 'c11', 'c22', 'c33'); model.material('comp').propertyGroup('def').set('density', '1800');注意各向异性参数得按C11、C22、C33顺序填，搞反了会导致声速计算结果抽风。声源时间参数控制着脉冲宽度，0.5e-6秒的脉宽对应1MHz中心频率——这是经过实测验证过的最佳参数。
求解器配置最容易翻车。建议先用频域求解找共振点，再用时域模块抓波形：
model.study('std1').feature('freq').set('plist', 'linspace(0.5e6,2e6,20)'); model.study('std2').feature('time').set('tlist', '0:1e-7:5e-6');时域计算必须开阻尼牛顿法，不然数值振荡能让你怀疑人生。网格剖分要特别注意固-气交界面，这里需要手动加密到波长1/5以下：
mesh.fun('custom').set('expr', 'sqrt((x-0.05)^2+(y-0.05)^2)'); mesh.automatic(true); mesh.elementSize('elementSize').set('custom', 'on');跑完仿真别急着看云图，先在接收点位置设个探针：
model.result.dataset('pt1').set('data', 'cpt1'); model.result.table('tbl1').set('data', 'pt1');时域信号得做FFT才能看出有效频率成分。实测发现层间脱粘缺陷会导致2MHz分量增强，这跟实际探伤中的频移现象吻合。不过要注意，空气耦合的接收信号幅值通常只有发射信号的0.1%不到，后处理时记得把纵坐标调成对数刻度。
最后说个血泪教训：千万别在固体域和空气域之间用自动生成的接触对，COMSOL会默认创建连续性边界条件，但实际物理过程应该用阻抗边界条件手动处理声压耦合。之前因为这个错误导致仿真结果比实测信号强了三个数量级，排查了整整两天才发现问题所在...
