2. 单自由度阻尼系统建模

要计算阻尼振动系统的响应，第一步不是套公式，而是把真实的机械装置抽象成一个可计算的力学模型。单自由度阻尼系统是最基础、也最常用的模型：它把复杂的结构简化为一个集中质量、一个弹簧和一个阻尼器，三者并联（或按实际连接方式组合）承受外力。本节围绕这个模型的建立过程展开，重点讲清三件事——力学模型怎么画、运动方程怎么列、阻尼力与阻尼系数的物理意义是什么。

一、力学模型的三个基本元件

单自由度系统只允许沿一个方向运动，通常取质量块的位移为广义坐标。模型由三类理想元件构成：

• 质量 m：代表系统的惯性，储存动能。它把力转化为加速度，遵循牛顿第二定律。
• 弹簧 k：代表系统的弹性恢复能力，储存势能。它把位移转化为恢复力，力与位移成正比。
• 阻尼器 c：代表系统的耗能能力，把机械能转化为热能耗散。它把速度转化为阻尼力，力与速度成正比。

实际工程中，质量可能是机床的工作台、汽车的车身、建筑楼层；弹簧可能是螺旋弹簧、梁的弯曲刚度、隔振垫；阻尼器可能是液压减振器、材料内阻尼、空气阻力。建模的关键是判断这些实际部件分别对应哪个理想元件，并确定它们的等效参数。

二、坐标约定与受力分析

列方程前必须先约定坐标系，否则正负号必然混乱。本节统一约定：取静平衡位置为位移原点 x=0，向下为正；速度 v=ẋ、加速度 a=ẍ 均以向下为正；所有力以向下为正方向参与求和。

设质量块在静平衡位置时弹簧已有静伸长 δ_st，满足 mg = kδ_st。当质量块从平衡位置向下移动 x 时，弹簧总伸长量为 δ_st+x，弹簧力向上（负方向），大小为 k(δ_st+x)；阻尼力与速度方向相反，向下运动时阻尼力向上，大小为 cẋ；重力向下为 mg。对质量块应用牛顿第二定律：

mẍ = mg − k(δst + x) − cẋ

代入 mg = kδ_st，静伸长项与重力项相互抵消，得到：

mẍ + cẋ + kx = 0        （自由振动）
mẍ + cẋ + kx = F(t)      （受迫振动）

这个消去重力的过程说明一个重要结论：以静平衡位置为原点时，重力和静弹簧力自动抵消，方程中不再出现 mg。这正是建模时选静平衡位置作原点的原因。若改用弹簧原长位置作原点，方程会多出常数项，虽然等价但计算更繁琐。

三、阻尼力与阻尼系数的物理意义

阻尼力写成 F_d = −cẋ，负号表示它始终与速度方向相反，即始终做负功、消耗能量。阻尼系数 c 的单位是 N·s/m，物理含义是单位速度所产生的阻尼力大小，反映系统耗能的强弱。c 越大，振动衰减越快；c=0 时退化为无阻尼简谐振动。

工程中阻尼往往不是单一机构，需要折算成等效阻尼系数。常见来源包括：液压减振器节流阻力、材料内摩擦、结构连接面滑移摩擦、空气或流体粘性阻力。对于粘性阻尼假设，我们只保留与速度一次方成正比的部分，这也是方程保持线性的前提。若实际阻尼以干摩擦为主（库仑阻尼），力与速度符号有关而与大小无关，方程变为非线性，需要分段处理，这已超出线性单自由度模型的范围。

为后续响应计算做准备，通常把方程两边除以 m，化为标准形式：

ẍ + 2ζωn ẋ + ωn² x = f(t)
ωn = √(k/m)          无阻尼固有角频率
ζ = c / (2√(km))     阻尼比
cc = 2√(km)          临界阻尼系数

这里 ζ 是无量纲的阻尼比，它把 c 与系统自身的 m、k 联系起来，是判断系统响应形态的关键参数：ζ<1 为欠阻尼（振荡衰减），ζ=1 为临界阻尼（最快无振荡回到平衡），ζ>1 为过阻尼（缓慢无振荡返回）。同一个 c 值放在不同质量、不同刚度的系统上，效果完全不同，因此比较阻尼强弱应看 ζ 而不是单看 c。

四、建模示例：机床隔振台

设一台设备安装在隔振垫上，设备与台架总质量 m=120 kg，隔振垫等效刚度 k=4.8×10⁴ N/m，等效粘性阻尼系数 c=480 N·s/m。按上述步骤建模：

1. 取静平衡位置为原点，向下为正，确认三个元件并联承受同一位移。
2. 计算固有频率：ω_n=√(48000/120)=√400=20 rad/s，对应 f_n≈3.18 Hz。
3. 计算临界阻尼系数：c_c=2√(km)=2√(48000×120)=2×2400=4800 N·s/m。
4. 计算阻尼比：ζ=480/4800=0.1，属于欠阻尼，系统受扰后会以约 3.18 Hz 的频率振荡并逐步衰减。

这个示例说明：c=480 N·s/m 看似不小，但相对该系统只有 0.1 的阻尼比，衰减较慢。若希望快速稳定，可增大阻尼使 ζ 接近 0.6～0.7；若 ζ 超过 1，则系统不再振荡但回位变慢。建模阶段算出 ζ，就等于预判了后续响应的基本形态。

需要强调的是，等效参数的正确性直接决定计算可信度。刚度可通过静载位移测量估算（k=ΔF/Δx），质量可由称重确定，阻尼系数较难直接测量，工程上常通过自由振动衰减曲线或半功率带宽法反推 ζ，再换算为 c。建模不是一次写方程就结束，而是"抽象—列方程—参数辨识—校验"的循环，本节建立的方程是后续所有响应计算的起点。