2. 单自由度受迫振动系统的物理建模

上一节我们建立了自由振动的概念：系统一旦偏离平衡位置，就会在弹性力与惯性力的交替作用下按固有频率振荡。但真实工程结构很少"自己振给自己看"——电机转子有偏心、路面有起伏、风载荷有脉动、地基有振动，这些外部作用持续向系统注入能量，使振动呈现出与自由振动完全不同的面貌。本节的任务，就是把"弹簧—质量—阻尼器"这个经典原型放到外力作用下，完成从物理系统到数学方程的完整映射。

一、原型系统的物理构成与坐标约定

取一个最简化的单自由度模型：质量为 m 的刚体块置于光滑水平面上，通过刚度为 k 的线性弹簧与固定壁相连，同时连接一个阻尼系数为 c 的粘性阻尼器，质量块上作用一个随时间变化的外力 F(t)。

建模前必须先固定参照约定，否则正负号必然混乱。本节统一规定：以弹簧原长时质量块所在位置为坐标原点 x=0，向右为位移正方向；所有力均按其在正方向上的投影取正负。这一约定将贯穿全文及图示，后续任何符号判断都以此为准。

需要强调的是，这个原型并非"玩具"。机床工作台、建筑隔振层上的设备、汽车悬架上的簧上质量、精密仪器的减振基座，在只关心主导振型时都可以归约为"一个等效质量 + 一个等效弹簧 + 一个等效阻尼 + 一个等效激励"的四要素结构。因此本节推导出的方程具有普遍的工程适用性。

二、四类力的物理意义辨析

对质量块做受力分析，共出现四种力，它们各自扮演不同角色：

• 弹性力（弹簧恢复力）：F_k = −k x。弹簧被拉伸 x>0 时，力指向左，试图把质量拉回平衡位置；被压缩时方向相反。负号体现的是"恢复"本质——它永远与位移反向，是系统能够振荡的根本原因。
• 阻尼力：F_c = −c ẋ。粘性阻尼力与速度成正比且方向相反，它不储存能量，而是把机械能转化为热耗散掉。注意它依赖的是速度而非位移，这是与弹性力的关键区别。
• 惯性力：在惯性系中直接对质量块用牛顿第二定律时，m ẍ 出现在方程左侧，代表"维持加速度所需的合力"；若采用达朗贝尔原理，则可写成 −m ẍ 并称为惯性力，与其余力在形式上平衡。两种写法等价，但物理视角不同，建模时选一种并保持一致即可。
• 激振力：F(t) 是外部施加的已知时间函数，是系统能量的输入源。它可以是简谐的（如偏心转子产生的 F₀sinωt）、周期的、瞬态的或随机的，其具体形式决定了后续求解方法的选择。

三、从物理系统到运动方程的映射

按牛顿第二定律，沿 x 正方向列合力方程：

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

移项整理为标准形式：

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

其中：m —— 等效质量（kg）
      c —— 等效粘性阻尼系数（N·s/m）
      k —— 等效刚度（N/m）
      F(t) —— 等效激振力（N）
      x(t) —— 位移响应（m）

这个二阶常系数线性非齐次微分方程就是单自由度受迫振动的动力方程。它的结构极具规律性：左侧依次是惯性项、阻尼项、弹性项，分别对应 ẍ、ẋ、x 三个状态量；右侧是外部激励。方程两边量纲均为牛顿，这一点可作为建模正确性的快速自检。

为便于分析，通常再引入两个特征参数：固有频率 ω_n = √(k/m) 与阻尼比 ζ = c/(2√(k m))，方程化为 ẍ + 2ζω_n ẋ + ω_n²x = F(t)/m。这一步不是数学游戏——它把四个物理参数压缩为两个无量纲/半量纲特征量，直接决定了系统的响应形态。

四、建模流程与一个教学算例

把上述过程固化为可复用的建模路径：确定等效质量与坐标原点 → 识别等效刚度与等效阻尼 → 明确激振力形式与作用点 → 按统一正方向列力平衡 → 整理为标准二阶方程 → 提取 ω_n、ζ 并校核量纲与符号。

设一教学示例（参数为合理假设，非实测）：m = 2 kg，k = 800 N/m，c = 8 N·s/m，激振力 F(t) = 10 sin(25t) N。代入得：

2·ẍ + 8·ẋ + 800·x = 10·sin(25t)

ωn = √(800/2) = 20 rad/s
ζ  = 8 / (2·√(2·800)) = 8/80 = 0.1
激励频率 ω = 25 rad/s，频率比 r = ω/ωn = 1.25

由此可预判：系统欠阻尼（ζ<1），自由响应会衰减振荡；激励频率高于固有频率（r>1），处于共振峰右侧，稳态振幅将小于静态位移 F₀/k = 0.0125 m，且响应相位滞后于激励。这些判断在下一节求解时会被定量验证，但建模阶段就应能定性说出方向——这正是物理建模的价值所在。

五、建模中的常见陷阱

实际建模时，符号错误远比数学求解错误更常见。三个高频问题值得警惕：其一，坐标原点与正方向未声明，导致弹簧力与阻尼力符号互相矛盾；其二，把重力漏算或重复计算——若以静平衡位置为原点，重力与弹簧静变形产生的力相互抵消，方程中不再出现 mg，这是最容易被忽视的简化；其三，等效参数拼凑错误，例如多个弹簧并联时刚度相加、串联时刚度倒数相加，阻尼器同理，混用会直接改变 ω_n 与 ζ。

建议在列出方程后做两项自检：量纲检查（每一项是否都是牛顿）与极限检查（令 F=0 应退化为自由振动方程 m ẍ+c ẋ+k x=0；令 c=0、k=0 应退化为 m ẍ=F 的纯惯性响应）。两项都通过，方程基本可信。