2. 振动系统力学基础

在 Ansys Mechanical 中做振动与模态分析,软件界面上的“求解”按钮只是最后一步;真正决定结果是否可信的,是分析之前对系统力学模型的判断。一个设备、支架或管路,究竟该简化为单自由度(SDOF)还是多自由度(MDOF)?固有频率由哪些参数决定?阻尼在仿真里如何体现?模态振型说明什么?共振判据怎么定?这些问题不解决,网格再细、求解再快,也可能得到没有工程意义的数字。本节把这条力学主线梳理清楚,为后续建模、求解和结果判读提供依据。

一、单自由度系统:振动分析的最小模型

单自由度系统只有一个独立坐标描述运动,典型模型是质量块—弹簧—阻尼器。设质量 m、刚度 k、黏性阻尼系数 c,外力 F(t) 作用在质量上,位移为 x(t),按牛顿第二定律可得运动方程:

m·x''(t) + c·x'(t) + k·x(t) = F(t)

无阻尼自由振动时 c=0、F=0,方程退化为 m·x'' + k·x = 0,其解为简谐振动,圆频率 ωn = √(k/m),固有频率 fn = ωn/(2π)。这个式子说明:固有频率只由质量和刚度决定,与初始位移、初始速度无关;质量增大频率下降,刚度增大频率上升。工程上常把设备简化为“质量—弹簧”系统做快速估算,例如一台质量为 500 kg 的泵组安装在刚度 k = 3.5×10⁶ N/m 的隔振器上,则 ωn = √(3.5×10⁶/500) ≈ 83.7 rad/s,fn ≈ 13.3 Hz。这是示例参数下的估算值,用于说明参数关系,不代表任何具体机组的实测结果。

阻尼的作用是让自由振动幅值随时间衰减。黏性阻尼下,衰减包络为 e^(−ζωn·t),其中 ζ = c/(2√(mk)) 称为阻尼比。ζ 越大衰减越快;ζ=1 为临界阻尼,系统不再振荡。多数金属结构 ζ 在 0.1%~2% 量级,橡胶隔振器可达 5%~15%,土建基础可能更高。Ansys 模态分析中常用恒定阻尼比或材料阻尼输入,若忽略阻尼,模态频率为无阻尼固有频率,而实际共振峰位置和幅值都会受阻尼影响。

二、多自由度系统:耦合、模态与振型

真实设备几乎都是连续体,离散后成为多自由度系统。以两个质量块通过弹簧串联为例,运动方程写成矩阵形式:

[M]{x''} + [C]{x'} + [K]{x} = {F(t)}

其中 [M]、[C]、[K] 分别为质量、阻尼、刚度矩阵。无阻尼自由振动对应广义特征值问题:[K]{φ} = ω²[M]{φ}。求解得到 n 个特征值 ωi² 和 n 个特征向量 {φi},即 n 阶固有频率和对应振型。振型描述的是系统在某一阶固有频率下各自由度位移的相对比例和方向,是一个“形状”,不是绝对幅值。Ansys 模态分析输出的正是这些频率与振型。

多自由度系统与单自由度最大的区别在于耦合:各自由度之间通过刚度或质量矩阵的非对角项相互影响。当两个子系统的固有频率接近时,会出现模态耦合,振型表现为两个方向同时运动,能量在两者之间来回传递。工程上遇到“某方向振动特别大、另一方向也同步放大”时,往往就是耦合模态在起作用。

单自由度与两自由度系统:参数、频率与振型关系 单自由度(SDOF) k c m x(t) ωn = √(k/m) fn = ωn / 2π ζ = c / (2√(mk)) 两自由度(MDOF) k1 m1 k2 m2 振型示意 1阶:同向,低频 2阶:反向,高频

三、固有频率、阻尼与共振判据

固有频率是系统的“身份特征”。Ansys 模态分析求解的就是无阻尼或含阻尼的固有频率与振型。工程上最关心的是:设备在运行转速对应的激励频率附近,是否存在固有频率。若存在,就可能发生共振。共振时即使激励力不大,响应也会被放大,放大倍数近似由阻尼比决定:Q ≈ 1/(2ζ)。ζ=1% 时 Q≈50,ζ=5% 时 Q≈10。这说明阻尼对共振幅值有决定性影响,也解释了为什么模态分析中阻尼设置不能随意忽略。

工程判据通常从频率分离入手。常见做法是要求激励频率与固有频率之间保持一定裕度,例如:

  • 一般旋转设备:|f_exc − fn| / fn ≥ 15%~20%;
  • 重要设备或高振动敏感设备:裕度提高到 20%~30%;
  • 若无法避开,则需评估阻尼、动平衡、隔振或增加刚度/质量以移频。

这里的 f_exc 是激励频率。对转速 n(r/min)的旋转设备,转频 f = n/60;若存在叶片通过频率、啮合频率或两倍转频,还要分别计算。例如转速 1800 r/min,转频 30 Hz;若某阶固有频率为 32 Hz,频率比接近 1,属于高风险;若固有频率为 45 Hz,则裕度约 50%,相对安全。这些数字是示例,用于说明判据的用法。

在 Ansys Mechanical 中,模态分析前先确认单位制、材料密度和弹性模量是否与手算模型一致。密度决定质量矩阵,弹性模量决定刚度矩阵,两者直接决定固有频率。若手算 fn 与仿真结果偏差超过 10%,优先检查约束、连接刚度和材料参数,而不是先怀疑网格。

四、从力学模型到仿真设置

把上述概念落到 Ansys Mechanical 操作上,可以按以下思路组织:

  1. 明确分析对象和边界:设备本体、支架、基础、隔振器分别如何简化;固定约束加在哪里,弹性连接用什么单元模拟。
  2. 建立质量与刚度来源:材料密度、附加质量点、弹簧刚度、接触刚度、螺栓连接刚度。刚度被高估会使频率偏高,质量被漏掉会使频率偏高。
  3. 设置模态提取:通常提取前 6~10 阶,覆盖关注频段;若关心高频局部模态,需提高提取阶数并检查网格是否足够。
  4. 输入阻尼:模态分析可用恒定阻尼比;谐响应分析中可设置材料阻尼或 Rayleigh 阻尼 α、β。
  5. 结果判读:查看各阶频率、振型方向、参与系数和有效质量。若某方向有效质量占比低,说明该方向可能未被充分激发或约束不合理。

以一台风机为例,若模态结果显示第 3 阶频率为 24.5 Hz,振型为电机—叶轮整体俯仰,而风机转频为 25 Hz,则频率裕度仅约 2%,属于共振风险区。此时可考虑增加支架刚度、调整支撑位置或增加阻尼。修改后重新求解,若该阶频率移到 30 Hz 以上,裕度超过 20%,则满足一般工程判据。这里给出的频率和裕度均为示例,实际项目应以现场转速、实测模态或详细仿真为准。

不要把“模态分析结果”直接当作“振动响应结果”。模态分析给出的是固有频率和振型,不包含激励大小和相位;要得到实际振幅、应力和传递率,还需要做谐响应、随机振动或瞬态分析,并正确输入激励谱和阻尼。

总结来看,单自由度模型帮助理解频率与刚度、质量、阻尼的基本关系;多自由度模型揭示模态、振型和耦合;固有频率与阻尼共同决定共振风险;工程判据则把频率裕度、阻尼比和激励频率联系起来。掌握这条主线,才能在 Ansys Mechanical 中做出有依据的简化、设置和判断,而不是只依赖软件默认值。