混频仿真与无损检测：基于Comsol固体力学分析的位移傅立叶变换研究

133混频 comsol 固体力学相关，混频无损检测，两个不同位置不同方向，不同频率混频仿真。 并对位移做了傅立叶变换

混频检测这玩意儿在无损检测里属于那种既骚气又实用的存在。最近折腾了个133混频的COMSOL仿真，核心思路就是让两个不同位置、不同方向的激励源以不同频率干活，看看材料内部怎么"打架"。

建模时直接在固体力学模块里搞事情，两个压电片分别贴在试样两侧。左边那个用131Hz的正弦载荷竖直方向怼，右边玩的是133Hz的水平方向输出——这俩频率差刚好能生成2Hz的低频振动，专门用来钓材料缺陷这条大鱼。

代码部分其实挺有意思，COMSOL的模型树节点用Append方法比图形界面快得多。比如设置第二个压电片的边界条件：

model.component('comp1').physics('solid').feature('bnd1', 2).set('LoadType', 'Harmonic'); model.component('comp1').physics('solid').feature('bnd1', 2).set('F0', '133[Hz]');

这里故意用字符串表达式直接写频率参数，方便后续参数化扫描的时候改数值不费劲。注意载荷方向得用方向矢量来定义，x分量给0.8，y分量0.2这种非正交配置能有效激发复杂模态。

133混频 comsol 固体力学相关，混频无损检测，两个不同位置不同方向，不同频率混频仿真。 并对位移做了傅立叶变换

仿真跑完后处理才是重头戏。位移场的傅里叶变换不能直接用软件自带功能，得把时域数据导出来自己处理。Python脚本这么写：

fft_result = np.fft.fft(displacement_data) freqs = np.fft.fftfreq(len(fft_result), 1/sampling_rate) peak_idx = np.argmax(np.abs(fft_result[1:])) + 1 # 跳过零频 print(f"主频分量在 {freqs[peak_idx]:.2f} Hz")

这里有个坑：采样率必须满足奈奎斯特准则，比如最高激励频率133Hz的话，采样至少得266Hz以上。实际跑仿真的时候建议设个500Hz保平安。

数据可视化时发现个反直觉的现象——在缺陷边缘处，混频产生的低频分量幅值反而比完整区域高。后来想明白这是不同模态的叠加效果，就像两股水流交汇处会形成漩涡一样，材料内部应力波相互干涉产生了"能量洼地"。

这种检测方法最大的优势在于信噪比。单独测131Hz或133Hz的响应可能被环境噪声淹没，但他们的差频信号在频谱上就像深夜里的霓虹灯一样显眼。实际操作时建议扫频范围控制在±5Hz，既能捕捉非线性效应又不至于跑偏。