首页/文章/ 详情

振动FFT变换

7月前浏览1459

看不懂 FFT,就无法理解振动的本质。很多初学者只会点击软件上的 "Calculate" 按钮,但作为 NVH 工程师,我们需要理解每一个参数背后的物理意义以及它对结果的扭曲。

下面将分三个层次为你拆解:核心逻辑参数设置及其影响、以及结果解读


第一层:核心逻辑——为什么要变换?

在振动试验(如台架测试或整车路试)中,传感器(加速度计)记录的是 时域信号(Time Domain):也就是振动幅值随时间变化的曲线。

  • 时域看到的是现象: “车在抖”、“声音很响”。

  • 频域看到的是根源: “这是由 50Hz 的转子不平衡引起的”、“这是 1200Hz 的齿轮啮合音”。

FFT 的本质就是数学上的棱镜,它将混合在一起的复杂时域波形,拆解成一个个单一频率的正弦波分量。


第二层:关键参数设置及其影响(这是重难点)

在进行 FFT 分析时,无论是使用 Python、Matlab 还是专业的 NVH 软件(如 LMS Test.Lab, Brüel & Kjær, Dewesoft),你都需要配置以下几个核心参数。这几个参数之间存在着**“测不准原理”般的制约关系**。

1. 采样率 (F_s, Sampling Rate)

  • 定义: 每秒钟记录多少个数据点。

  • NVH 铁律(奈奎斯特采样定理): 采样率必须至少是你关注的最高频率(F_{max})的 2 倍。


    F_s > 2*  F_{max}

  • 设置建议: 为了保证波形还原度和防止混叠(Aliasing),工程上通常设置为关注频率的 2.56 倍或 5 倍以上。例如,关注 2000Hz 以内的振动,采样率至少设为 5120Hz 或更高。

2. 频率分辨率 (Δf, Frequency Resolution)

  • 定义: 频谱图上两个相邻数据点之间的频率间隔。它是你能区分两个频率成分的最小能力。

  • 公式:


    Δf=1/T =F_s/N


    (其中 T是采样时长,N 是采样点数/块大小)

  • 参数影响(导师划重点):

    • 想要分得更细(Δf更小),就必须采样更长的时间 (T变大)。

    • 例子: 如果你想区分 50.0Hz 和 50.5Hz 的振动,你的 Δf 必须小于 0.5Hz,这意味着你的采样时间块 T 必须大于 2 秒。

3. 窗函数 (Windowing)

  • 为什么需要: FFT 假设信号是无限周期循环的。如果你截取的一段信号(Frame)的首尾不连续,变换后会产生频谱泄漏(Spectral Leakage),导致尖峰变胖,幅值不准,甚至掩盖微弱信号。

  • 常用窗函数对比:

窗函数类型特点适用场景NVH 常用度
Rectangular (矩形窗)其实就是不加窗,频率分辨率最高,但泄漏最严重。用于瞬态信号(如锤击法测试,信号在采样时间内衰减为0)。⭐⭐
Hanning (汉宁窗)很好的通用性,抑制泄漏能力强,但会损失幅值(需修正)。用于连续随机信号(如路噪、风噪、白噪声振动)。⭐⭐⭐⭐⭐
Flat Top (平顶窗)幅值精度最高,但频率分辨率最差(峰会变宽)。用于校准或主要关注幅值精度的正弦振动。⭐⭐⭐

4. 重叠率 (Overlap) 与 平均 (Averaging)

  • 背景: 振动信号通常有随机性。为了让结果更稳定,我们通常对多帧数据的 FFT 结果取平均。

  • Overlap: 为了充分利用数据(特别是加了 Hanning 窗后,数据头尾被压低,信息丢失),我们会让处理块之间重叠(通常 50% 或 66%)。

  • 平均方式:

    • 线性平均 (Linear/Power Avg): 消除随机噪声,突出周期性信号。

    • 峰值保持 (Peak Hold): 记录测试全过程中出现的每个频率的最大值(常用于寻找最坏工况)。


第三层:FFT 变换后能获取的振动信息

拿到频谱图后,NVH 工程师需要从中解读出什么?

1. 主要频率成分(Identify Dominant Frequencies)

  • 旋转阶次 (Orders): 如果振动频率随着转速变化(例如 $f = \frac{RPM}{60} \times k$),这是由旋转部件引起的(不平衡、不对中、叶片通过频率)。

  • 固定频率 (Resonance): 如果不管转速怎么变,某个频率的峰值一直存在且很高,这通常是结构的固有频率(共振)

2. 幅值信息 (Amplitude)

  • 峰值 (Peak): 该频率下的最大振动量。

  • 有效值 (RMS): 这是一个能量概念,代表该频段内的总能量。在宽带随机振动测试中非常重要。

  • PSD (功率谱密度): 当处理随机振动(如路谱模拟)时,单纯的 FFT 幅值会受 Δf影响而波动。必须转换为 PSD (g^2/Hz),它是标准化的能量密度,与频率分辨率无关,便于不同测试标准间的对比。

3. 谐波失真 (Harmonics)

  • 除了基频(比如电机转速频率),是否还有 2倍频、3倍频?丰富的谐波通常意味着机械系统中存在“非线性”,如间隙(Looseness)或摩擦。


第三部分:FFT 之后我们能获取什么信息?

拿到频谱图后,NVH 工程师主要看三点:

  1. 主要频率(Dominant Frequency): 谁是罪魁祸首?

    • 如果是转速的整数倍(1X, 2X, 3X...),通常与旋转部件有关(不平衡、不对中)。

    • 如果是电源频率(50Hz/60Hz)及其倍频,可能是电气噪声或磁致伸缩。

    • 如果是固定频率且不随转速变化,通常是结构的共振频率(Natural Frequency)

  2. 幅值(Amplitude): 问题有多严重?

    • 对应各个频率点的振动能量大小。

  3. 边带与谐波(Sidebands & Harmonics): 细节诊断。

    • 例如齿轮箱故障,通常表现为啮合频率两侧出现很多小的边带。


第四部分:实际案例分析——电机异常振动诊断

背景: 某款电动汽车的空调压缩机在 3000 RPM (50 Hz) 运行时,车内感到明显抖动。

1. 时域数据观察:

数据看起来就是一团杂乱的波浪,幅值大约在 $\pm 2g$ 之间波动。肉眼看不出规律。

2. FFT 参数设置:

  • 采样率 (F_s): 5120 Hz (保证能看到 2500Hz 以内的信号)。

  • 点数 (N): 4096 点。

  • 分辨率 (Δf):5120 / 4096 = 1.25 Hz (足够区分主要频率)。

  • 窗函数: Hanning。

3. 频谱分析结果:

  • 现象 A: 在 50 Hz 处有一个高达 1.5g 的尖峰。

    • 分析: 50 Hz 正好等于电机转速 (3000 RPM / 60)。

    • 结论: 这是一阶(1st Order)振动,通常直接指向转子动不平衡(Unbalance)

  • 现象 B: 在 500 Hz 处有一簇较小的峰值,且伴有边带。

    • 分析: 压缩机内部有 10 个叶片,50 * 10 = 500 Hz

    • 结论: 这是叶片通过频率(Blade Pass Frequency),属于流体脉动引起的振动,幅值较小,可忽略。

最终对策: 对压缩机转子进行动平衡校正,削减 50Hz 处的能量。



来源:ABAQUS仿真世界
振动非线性电源通用汽车MATLABpythonLMS电机NVHAlias试验电气
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2026-01-22
最近编辑:7月前
yunduan082
硕士 | 仿真主任工程... Abaqus仿真世界
获赞 203粉丝 392文章 486课程 0
点赞
收藏
作者推荐

基于ABAQUS的PET瓶结构优化与轻量化设计

摘要采用Abaqus/Explicit程序进行了显式非线性分析,计算了PET瓶的屈曲载荷和应力分布。根据数值模拟结果,对PET瓶的结构进行了优化,以提高其屈曲载荷。同时,利用Abaqus获得的PET瓶应力云图,对PET瓶的塑料分布进行了优化,提高了PET材料的利用效率,降低了瓶子重量。这为PET瓶的轻量化设计提供了新思路和新方法。关键词:PET瓶;有限元分析;ABAQUS;结构优化;轻量化设计;屈曲载荷引言PET瓶是一种广泛使用的饮料包装容器,具有质量轻、强度高、透明度好等优点。然而,在运输、储存和使用过程中,PET瓶需要承受轴向压力、内部气压等载荷,容易发生屈曲失稳或局部破坏。为了满足强度要求的同时实现轻量化,减少材料消耗和环境影响,优化PET瓶结构已成为重要研究方向。传统设计主要依赖物理原型实验,周期长、成本高。本文利用有限元软件ABAQUS,建立PET瓶的三维模型,进行非线性屈曲分析和应力计算。通过模拟结果指导结构优化和材料分布调整,实现瓶子轻量化。有限元模型的建立几何模型本文以典型600ml碳酸饮料PET瓶为研究对象。瓶子几何包括瓶口、瓶颈、瓶肩、瓶身和瓶底。瓶壁厚度不均匀:瓶肩和瓶底较厚(约0.5-0.8mm),瓶身较薄(约0.2-0.3mm)。模型忽略螺纹等细小特征,以简化计算。材料模型PET材料为高分子聚合物,具有非线性弹塑性行为。模拟采用弹性-塑性本构模型,材料参数如下:密度:1.37g/cm³弹性模量:约2.5-3.5GPa(考虑温度影响)泊松比:0.4考虑大变形几何非线性。载荷与边界条件主要载荷:轴向顶载(模拟堆垛压力),典型值为4-8kN。边界条件:瓶底固定支撑,瓶口施加向下位移或集中力。内部压力:对于碳酸饮料瓶,考虑0.4-0.6MPa内压(可选)。分析类型:ABAQUS/Explicit显式动态分析,适合处理高度非线性的屈曲和接触问题,避免隐式分析收敛困难。网格划分采用壳单元(S4R)或三维实体单元(C3D8R)划分网格。应力集中区(如瓶肩过渡、瓶底瓣状结构)加密网格,总单元数约5-10万。屈曲分析结果与结构优化通过显式非线性分析,计算原设计PET瓶的轴向屈曲临界载荷。结果显示:原瓶屈曲载荷约为50N。屈曲模式:瓶身或瓶肩局部凹陷,伴随整体失稳。最大VonMises应力出现在瓶肩和瓶底支撑环区域。基于分析结果,进行结构优化:增加瓶肩圆角半径,减少应力集中。在瓶身添加纵向或环向筋条(ribs),增强抗屈曲能力。优化瓶底“瓣状”或“香槟底”结构,提高底部刚度。优化后模拟结果:屈曲载荷提高15%-30%。失稳模式延迟,瓶子在标准堆垛载荷(例如8层堆码)下保持稳定。应力分布优化与轻量化设计ABAQUS输出详细的应力云图(VonMises应力、principalstress等),清晰显示瓶子各区域应力水平:高应力区:瓶肩过渡区、瓶底支撑点(应力>30-50MPa)。低应力区:瓶身上部和中部(应力<10MPa)。根据应力分布,进行塑料厚度优化(变量厚度设计):保持或略增高应力区厚度。显著减薄低应力区壁厚(例如从0.3mm减至0.2mm)。轻量化结果:瓶子重量减少10%-20%(例如从25g减至20g以下)。材料利用率提高,应力分布更均匀。优化后瓶子仍满足顶载、跌落和内压等标准要求。该方法实现了“等强度轻量化”,减少原材料消耗,降低生产成本和碳足迹。结论本文利用ABAQUS/Explicit对PET瓶进行非线性有限元分析,准确预测了屈曲行为和应力分布。通过模拟指导的结构优化和材料分布调整,显著提高了瓶子强度并实现了轻量化。该数值方法高效、经济,避免了大量物理实验,为PET瓶设计提供了科学工具。来源:ABAQUS仿真世界

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习计划 福利任务
下载APP
联系我们
帮助与反馈