4. 微分方程建立与标准化

在弹簧质量阻尼系统中，我们真正关心的往往不是“方程长什么样”，而是“方程里的系数分别代表什么物理能力”。建立微分方程只是第一步，把它整理成标准形式，才能让固有频率 ω_n 与阻尼比 ζ 这两个关键参数浮出水面，进而判断系统是过阻尼、临界阻尼还是欠阻尼。本节围绕“从受力到标准式”的完整推导展开，并给出可直接套用的标准化步骤与示例计算。

一、先定参照系与符号约定

为避免正负号混乱，本节统一约定：取静平衡位置为位移原点 x=0，竖直向下为正方向；速度 v=ẋ，加速度 a=ẍ 均对时间 t 求导。弹簧刚度 k、阻尼系数 c、质量 m 均为正值。外力 f(t) 以向下为正。这一约定在后续所有公式与图示中保持一致。

在静平衡位置，重力 mg 与弹簧静伸长产生的弹力相互抵消，因此以平衡位置为原点列方程时，重力项与静弹力项可以相消，方程中只保留“动态”的弹簧恢复力 −kx 与阻尼力 −cẋ。这是建立方程时最容易出错的地方，务必先做静平衡消项，再写动态方程。

二、从受力分析到二阶微分方程

对质量块做受力分析：向下的外力 f(t)，向上的动态弹簧力 kx，向上的阻尼力 cẋ。由牛顿第二定律：

m·ẍ = f(t) − c·ẋ − k·x

移项整理为“未知量在左、激励在右”的常规形式：

m·ẍ + c·ẋ + k·x = f(t)

这就是弹簧质量阻尼系统的二阶常系数线性常微分方程。它描述了位移 x(t) 在外力作用下的动态响应。注意：若外力为零，则右端为 0，方程描述自由振动；若外力为简谐形式，则描述受迫振动。

三、标准化：除以质量，提取 ω_n 与 ζ

标准形式的目的是让二阶项系数为 1，使两个参数具有明确的物理量纲与几何意义。将上式两边同除以 m：

ẍ + (c/m)·ẋ + (k/m)·x = f(t)/m

对照二阶系统标准形式：

ẍ + 2ζω_n·ẋ + ω_n²·x = f(t)/m

逐项比较系数，得到两个定义式：

• 固有频率：ω_n² = k/m，即 ω_n = √(k/m)，单位 rad/s。它只由质量与刚度决定，反映系统“自身振动的快慢”，与阻尼无关。
• 阻尼比：2ζω_n = c/m，即 ζ = c/(2√(km)) = c/(2mω_n)，无量纲。它反映实际阻尼相对于“临界阻尼”的比例。

其中临界阻尼系数 c_c = 2√(km) = 2mω_n，因此 ζ = c/c_c。当 ζ=1 时系统处于临界阻尼，恰好不发生振荡且最快回到平衡；ζ<1 为欠阻尼，会围绕平衡位置振荡衰减；ζ>1 为过阻尼，缓慢单调返回。这三个区间直接决定后续仿真的响应形态，因此标准化不是形式主义，而是判据的来源。

四、示例：一组可复算的参数

假设某隔振平台简化模型：质量 m = 2 kg，弹簧刚度 k = 200 N/m，阻尼系数 c = 8 N·s/m。这是一组便于手算的示例参数，用于说明标准化流程，并非实测数据。

先算固有频率：ω_n = √(200/2) = √100 = 10 rad/s。再算临界阻尼：c_c = 2√(km) = 2√(200×2) = 2×20 = 40 N·s/m。于是阻尼比 ζ = 8/40 = 0.2。由于 0<ζ<1，系统为欠阻尼，阶跃激励下会出现约 52.7% 的超调（超调量公式 M_p = e^{−πζ/√(1−ζ²)}，代入 ζ=0.2 得 e^−0.642≈0.527），并围绕平衡位置振荡衰减。若把 c 改为 40，则 ζ=1，系统临界阻尼；若 c 改为 60，则 ζ=1.5，过阻尼。可见标准化后，只需看 ζ 一个数就能预判响应形态。

把标准形式写全：ẍ + 2×0.2×10·ẋ + 10²·x = f(t)/2，即 ẍ + 4ẋ + 100x = f(t)/2。这个式子可以直接交给数值求解器（如四阶龙格-库塔）做仿真，状态变量取 x 与 ẋ 即可。

五、标准化操作清单与常见错误

实际建模时可按以下顺序执行，避免返工：

1. 确定坐标原点与正方向，写出静平衡消项后的动态受力式。
2. 按 m·ẍ + c·ẋ + k·x = f(t) 整理，检查各项量纲是否一致（N）。
3. 两边除以 m，使 ẍ 系数为 1。
4. 对照标准式读取 ω_n² = k/m 与 2ζω_n = c/m，反解出 ω_n、ζ。
5. 用 ζ 判定区间，并记录 c_c 作为后续调参基准。

完成本节标准化后，方程已具备统一形式，后续无论是解析求解还是数值仿真，都可以直接以 ω_n、ζ 为输入参数展开，响应曲线的形状也由此被预先“锁定”。