2. 单自由度系统物理建模

要对弹簧质量阻尼系统做响应仿真，第一步不是写代码，而是把真实机械结构抽象成一个可计算的数学模型。单自由度（SDOF）系统是最基础的建模对象：一个质量块只能沿某一方向平动，两侧分别连接弹簧与阻尼器。建模的核心任务，是把质量块、弹簧、阻尼器三者的力学行为，通过牛顿第二定律统一到一个微分方程中。本节围绕“部件力学特性 → 受力分析 → 力平衡方程 → 标准形式”这条主线展开，并给出可直接用于仿真的参数化结果。

一、三个部件的力学本构关系

建模前必须先明确每个部件“力与运动”的关系，这是列方程的依据。约定如下参照系：取质量块静平衡位置为位移原点，向右为位移 x 的正方向，速度 v = ẋ，加速度 a = ẍ，所有力以向右为正。

• 质量块（惯性元件）：本身不产生恢复力，只储存动能，其力学表现为惯性力 mẍ。在牛顿第二定律中，它出现在方程左侧。
• 弹簧（弹性元件）：力与变形量成正比，F_s = k·Δx，k 为刚度（N/m）。弹簧被拉伸时对质量块施加向左的恢复力，即 −kx。
• 阻尼器（耗能元件）：力与相对速度成正比，F_d = c·Δẋ，c 为阻尼系数（N·s/m）。阻尼器被拉伸（质量块向右运动）时施加向左的阻力，即 −cẋ。

需要注意，弹簧力和阻尼力都是“被动恢复力”，方向总与运动或变形方向相反，因此方程中带负号；而外部激励 F(t) 是主动力，方向按激励定义取正或负。这一符号约定必须在整门课程中保持一致，否则仿真结果会出现相位反转。

二、受力分析与牛顿第二定律

把质量块取为隔离体，沿 x 方向列出全部作用力：弹簧恢复力 −kx、阻尼力 −cẋ、外部激励 F(t)。根据牛顿第二定律 ΣF = ma，得到：

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

移项整理，把所有含位移的项放到左边，即得到单自由度系统的标准运动方程：

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

这个二阶常系数线性微分方程就是整个仿真的“被求解对象”。它的物理含义很清晰：惯性项 mẍ 抵抗加速度变化，阻尼项 cẋ 消耗振动能量，刚度项 kx 提供恢复趋势，三者之和等于外部激励。若 F(t) = 0，则为自由振动；若 F(t) 为简谐力，则为受迫振动。

三、从方程到仿真参数

方程本身还不能直接驱动仿真，需要把它改写为状态空间形式，即用两个一阶方程表示一个二阶方程。令状态变量 x₁ = x，x₂ = ẋ，则：

ẋ₁ = x₂
ẋ₂ = ( F(t) − c·x₂ − k·x₁ ) / m

这样任何数值积分器（如四阶龙格-库塔）都能直接求解。同时，方程中隐含两个决定系统“性格”的组合参数：

• 固有角频率 ω_n = √(k/m)，单位 rad/s，决定自由振动快慢；
• 阻尼比 ζ = c / (2√(km))，无量纲，决定振动衰减形态：ζ < 1 欠阻尼振荡，ζ = 1 临界阻尼，ζ > 1 过阻尼。

把方程两边同除以 m，可写成更通用的标准形式 ẍ + 2ζω_nẋ + ω_n²x = F(t)/m。仿真时通常先由 m、k、c 算出 ω_n 与 ζ，再据此选择时间步长和仿真时长。

四、一个可复现的建模示例

假设某设备上的减振台架可简化为单自由度系统，取一组教学示例参数（非实测）：质量 m = 2 kg，刚度 k = 800 N/m，阻尼系数 c = 8 N·s/m，初始位移 x(0) = 0.05 m，初始速度 ẋ(0) = 0，无外部激励。按上述公式计算：

ωn = √(800/2)      = 20 rad/s
ζ  = 8 / (2·√(800·2)) = 8 / 80 = 0.1
ωd = ωn·√(1−ζ²)    ≈ 19.90 rad/s
Td = 2π/ωd         ≈ 0.316 s

由于 ζ = 0.1 远小于 1，系统为欠阻尼，自由振动呈指数衰减的振荡，相邻峰值按 e^−ζω_nT_d ≈ e^−0.632 ≈ 0.53 的比例递减。仿真时若取时间步长 Δt = 0.001 s，一个振荡周期约含 316 个采样点，足以清晰刻画波形；仿真时长取 3～5 个周期（约 1～1.6 s）即可观察到明显衰减。这些数值可作为校验建模是否正确的基准：若仿真得到的周期明显偏离 0.316 s，或峰值不按约 0.53 倍衰减，就应回查 k、c、m 的取值或方程符号。

五、常见建模偏差

实际建模中，最容易出错的是符号与参照系。例如把弹簧力写成 +kx，会导致系统变成“负刚度”，仿真结果发散；把阻尼力方向与速度同向，则能量被人为注入，振幅越来越大。另一类偏差来自单位：刚度用 N/mm 而位移用 m，会使 k 放大 1000 倍，ω_n 随之失真。此外，若系统存在重力等恒定载荷，应先把静平衡位置设为原点，重力与弹簧静变形相互抵消，方程中不再出现重力项——这一步处理不当，会引入一个恒定的位移偏置，使仿真曲线整体偏移。

完成本节建模后，读者应能独立地把一个“质量块 + 弹簧 + 阻尼器”的物理对象，转化为 m·ẍ + c·ẋ + k·x = F(t) 这一标准方程，并由 m、k、c 计算出 ω_n、ζ 等特征参数，为后续的自由振动、受迫振动与数值求解打下基础。