第四节 单自由度系统自由振动：运动方程、固有频率推导、无阻尼与有阻尼情况
隔振系统的核心分析起点，是将其简化为单自由度（SDOF）弹簧-质量-阻尼模型。无论是被动隔振器中的橡胶垫、空气弹簧，还是主动隔振系统中的作动器与传感器配合，其动力学行为都首先通过这一基础模型来理解。本节将系统推导单自由度系统的自由振动方程，明确固有频率的物理意义与计算方法，并对比无阻尼与有阻尼两种情形下的响应特征。这是后续理解隔振效率、共振传递率以及主动控制策略的数学与物理基石。
上图展示了单自由度系统的物理模型、运动方程与响应特征之间的逻辑关系。质量块 m 通过刚度为 k 的弹簧和阻尼系数为 c 的阻尼器连接在基础上。当系统受到初始扰动（初位移或初速度）后，不再施加外部激励，系统将依靠自身弹性恢复力和阻尼力完成运动，这就是自由振动。理解这一过程，是分析隔振系统在瞬态扰动（如设备启停、路面冲击）下表现的前提。
4.1 运动方程的建立与固有频率推导
4.1.1 核心原理：牛顿第二定律与达朗贝尔原理
取质量块 m 为隔离体，设其相对于静平衡位置的位移为 x(t)，向下为正。在自由振动过程中，质量块受到三个力：弹簧恢复力 F_k = -kx（方向始终指向平衡位置）、阻尼力 F_c = -c·ẋ（与速度方向相反）、惯性力 F_i = -m·ẍ。根据牛顿第二定律或达朗贝尔原理，可得：
m·ẍ + c·ẋ + k·x = 0这是单自由度系统自由振动的通用运动方程。若将方程两边同时除以 m，并引入两个关键参数：
- 固有角频率：ωₙ = √(k/m)，单位 rad/s，反映系统固有的振动快慢属性，仅与质量和刚度有关，与初始条件无关。
- 阻尼比：ζ = c / (2√(km)) = c / (2mωₙ)，无量纲，反映阻尼相对于临界阻尼的大小。
则方程可改写为标准形式：
ẍ + 2ζωₙ·ẋ + ωₙ²·x = 0固有频率的物理意义在于：如果系统无阻尼且受到初始扰动，它将以此频率持续振荡。在隔振工程中，固有频率决定了隔振系统的起始隔振频率——通常当激励频率超过 √2 倍的固有频率时，隔振系统才开始发挥隔振作用。
4.1.2 判断方法：从物理参数到固有频率的计算路径
实际工程中，计算固有频率通常有两种路径：
- 理论计算法：已知设备质量 m 和隔振器总刚度 k（多个隔振器并联时 k = k₁ + k₂ + ...），直接代入 ωₙ = √(k/m)，再换算为频率 fₙ = ωₙ/(2π) Hz。
- 静态变形法：测量隔振器在设备自重下的静态压缩量 δ_st（单位 m），则 fₙ = (1/2π)·√(g/δ_st)。这是因为 k·δ_st = mg，代入 ωₙ = √(k/m) = √(g/δ_st)。该方法在工程现场非常实用，无需精确知道刚度和质量的具体数值。
4.1.3 案例拆解：精密仪器隔振平台的固有频率计算
背景与约束：某半导体检测设备质量 m = 800 kg，拟安装在 4 个相同的空气弹簧隔振器上，每个隔振器的垂直刚度 k_i = 40 N/mm。要求系统固有频率低于 5 Hz，以保证对 10 Hz 以上的地面微振动有足够衰减。
分析过程：4 个隔振器并联，总刚度 k = 4 × 40 N/mm = 160 N/mm = 160,000 N/m。代入公式：
ωₙ = √(k/m) = √(160000 / 800) = √200 ≈ 14.14 rad/s
fₙ = ωₙ / (2π) ≈ 14.14 / 6.283 ≈ 2.25 Hz结果与验证：计算得到固有频率约 2.25 Hz，满足低于 5 Hz 的设计要求。进一步验证静态变形法：每个隔振器承受 200 kg 重量，压缩量 δ_st = F/k_i = (200×9.8)/40000 ≈ 0.049 m = 49 mm，fₙ = (1/2π)√(9.8/0.049) ≈ 2.25 Hz，两种方法结果一致。
复盘经验：空气弹簧的刚度会随内部气压变化，实际安装后应通过激振测试或环境振动测试验证固有频率。若实测频率偏高，可适当降低气压或增加配重来调整。
4.2 无阻尼自由振动：等幅简谐运动
4.2.1 核心原理：能量守恒与简谐解
当 c = 0 时，运动方程简化为 ẍ + ωₙ²x = 0。这是一个二阶常系数齐次线性微分方程，其通解为：
x(t) = A·cos(ωₙt) + B·sin(ωₙt)或合并为振幅-相位形式：
x(t) = X·cos(ωₙt - φ)其中振幅 X = √(A² + B²)，相位角 φ = arctan(B/A)。A 和 B 由初始条件确定：若 t=0 时初位移 x₀、初速度 v₀，则 A = x₀，B = v₀/ωₙ。无阻尼自由振动的特点是：振幅恒定不变，系统总机械能（动能 + 弹性势能）守恒，振动永不衰减。这在物理上是一种理想化情形，实际系统中阻尼总是存在的。
4.2.2 判断方法：无阻尼模型的适用边界
在隔振工程中，无阻尼模型适用于以下场景：
- 系统阻尼比极小（ζ < 0.05），如金属弹簧隔振器、空气弹簧（无附加阻尼时）。
- 仅需估算固有频率和振型，不关心振幅衰减速率。
- 分析系统在共振点附近的峰值响应上限（无阻尼时共振振幅理论上趋于无穷，实际由阻尼限制）。
但需注意：无阻尼模型不能用于预测共振时的实际振幅，也不能用于评估系统在瞬态扰动后的稳定时间。对于隔振设计，通常仍需考虑至少 2%~5% 的阻尼比。
4.2.3 案例拆解：无阻尼假设下的共振风险评估
场景：某风机转速 1200 rpm（对应激励频率 20 Hz），安装在四个金属螺旋弹簧隔振器上，系统固有频率计算为 8 Hz。试评估启动过程中通过共振区时的风险。
分析：风机启动时转速从 0 逐渐升至 1200 rpm，激励频率连续扫过 8 Hz 固有频率。若按无阻尼模型，共振时振幅理论上无穷大，实际中金属弹簧阻尼比约 0.01~0.03，共振放大因子 Q = 1/(2ζ) ≈ 17~50。这意味着共振时传递的力幅是静态力的 17~50 倍，可能导致隔振器损坏或设备剧烈晃动。
措施：在弹簧隔振器上增加粘滞阻尼器（将 ζ 提升至 0.15~0.2），使共振放大因子降至 2.5~3.3；或采用变频启动控制，使风机快速通过共振转速区间（如 3 秒内从 0 升至额定转速），减少共振持续时间。
4.3 有阻尼自由振动：衰减振荡与临界阻尼
4.3.1 核心原理：特征根与三种阻尼状态
对于有阻尼系统，设解 x = e^(λt)，代入运动方程得到特征方程：
λ² + 2ζωₙλ + ωₙ² = 0解得特征根 λ = -ζωₙ ± ωₙ√(ζ² - 1)。根据阻尼比 ζ 的大小，系统呈现三种截然不同的响应状态：
| 阻尼状态 | 条件 | 特征根性质 | 响应形式 | 工程应用 |
|---|---|---|---|---|
| 欠阻尼 | 0 < ζ < 1 | 共轭复根 | x(t) = X·e^(-ζωₙt)·cos(ω_d·t - φ) | 绝大多数隔振系统 |
| 临界阻尼 | ζ = 1 | 重根 | x(t) = (C₁ + C₂t)·e^(-ωₙt) | 仪表指针、快速回零机构 |
| 过阻尼 | ζ > 1 | 两个负实根 | x(t) = C₁e^(λ₁t) + C₂e^(λ₂t) | 避免振荡的慢速回零机构 |
其中 ω_d = ωₙ√(1-ζ²) 称为有阻尼固有频率，略低于无阻尼固有频率。当 ζ 较小时（如 <0.2），ω_d 与 ωₙ 差异很小（差异 <2%），工程上常近似取 ω_d ≈ ωₙ。
4.3.2 判断方法：阻尼比估算与衰减率计算
工程中常用对数衰减率法从实测自由振动波形中估算阻尼比。设相邻两个正峰值分别为 x₁ 和 x₂，则对数衰减率 δ = ln(x₁/x₂)，阻尼比 ζ = δ / √(4π² + δ²)。当 ζ 较小时，可近似 ζ ≈ δ/(2π)。
此外，阻尼比也可通过半功率带宽法从频响曲线中获取，或通过阶跃响应超调量估算：ζ = -ln(M_p) / √(π² + ln²(M_p))，其中 M_p 为最大超调量百分比。
4.3.3 案例拆解：实测自由振动波形评估隔振器阻尼性能
背景：某光学平台采用橡胶隔振器，受到脉冲激励后，用位移传感器记录到自由振动波形。相邻三个正峰值分别为：x₁ = 4.2 mm，x₂ = 2.8 mm，x₃ = 1.9 mm。试评估该隔振器的阻尼比和系统固有频率。
分析过程：
- 计算对数衰减率：δ₁₂ = ln(4.2/2.8) = 0.405，δ₂₃ = ln(2.8/1.9) = 0.388，取平均 δ ≈ 0.397。
- 计算阻尼比：ζ = 0.397 / √(4π² + 0.397²) ≈ 0.063。
- 从波形时间轴读取振动周期 T ≈ 0.25 s，则有阻尼频率 f_d = 1/T = 4 Hz，无阻尼固有频率 fₙ = f_d / √(1-ζ²) ≈ 4 / √(1-0.004) ≈ 4.01 Hz。
结果与决策：该橡胶隔振器阻尼比约 6.3%，固有频率约 4 Hz，满足光学平台对 8 Hz 以上微振动的衰减要求。但阻尼比偏低，若地面存在 4~6 Hz 的持续扰动，可能引起明显共振放大。建议在平台与隔振器之间增加粘滞阻尼器，将阻尼比提升至 0.15 左右。
复盘经验：实测波形分析是验证隔振系统建模准确性的最直接手段。建议在设备安装完成后进行锤击试验或释放初位移试验，记录自由振动波形，提取实际固有频率和阻尼比，与设计值对比。若偏差超过 15%，应检查隔振器安装状态、附加刚度（如管道连接、电缆拖链）或质量分布是否与假设一致。
4.4 常见误区与总结复盘
4.4.1 误区清单
- 误区一：将固有频率与激励频率混淆。固有频率是系统固有属性，与外部激励无关；激励频率是外部扰动的频率。只有当两者接近时才会发生共振。
- 误区二：认为阻尼越大隔振效果越好。阻尼在抑制共振的同时，也会为高频振动提供传递路径。隔振设计需在共振抑制和高频衰减之间权衡。
- 误区三：忽略附加刚度的影响。设备与隔振器之间的管道、电缆、接地线等会引入附加刚度，使实际固有频率高于理论计算值。某案例中，一根粗电缆使系统固有频率从设计的 5 Hz 升至 8 Hz，导致隔振效果大幅下降。
- 误区四：将无阻尼固有频率直接用于有阻尼系统。当阻尼比大于 0.2 时，有阻尼固有频率与无阻尼值的差异不可忽略，应使用 ω_d = ωₙ√(1-ζ²) 进行修正。
4.4.2 总结复盘
单自由度自由振动分析为隔振系统提供了三个核心参数：固有频率 ωₙ（决定隔振起始频率）、阻尼比 ζ（决定共振放大倍数和衰减速率）、以及响应形式（决定系统在瞬态扰动下的行为）。在实际工程中，应遵循以下流程：
- 建立物理模型，明确质量、刚度、阻尼的等效值。
- 计算固有频率，验证是否满足隔振设计要求（通常要求 fₙ < f_excitation / √2）。
- 评估阻尼比，确保共振峰值可控且高频衰减不被过度牺牲。
- 通过实测自由振动或频响测试验证设计参数，必要时修正模型。