1. 课程导入与仿真目标

弹簧质量阻尼系统是机械振动课程中最基础、也最具工程穿透力的模型。它由三个理想元件构成：质量为 m 的刚体、刚度为 k 的线性弹簧、阻尼系数为 c 的黏性阻尼器。三者并联后，在外力 F(t) 作用下产生位移 x(t)。按牛顿第二定律，取静平衡位置为位移零点、向右为正方向，可得运动方程：

m·x''(t) + c·x'(t) + k·x(t) = F(t)
初始条件：x(0) = x0，x'(0) = v0

这个二阶常系数线性微分方程看似简单，却几乎覆盖了振动工程的全部核心概念：固有频率、阻尼比、共振、瞬态与稳态、幅频与相频特性。理解它，就理解了后续多自由度系统、连续体振动乃至振动控制的共同语言。

工程意义与典型应用场景

在工程中，弹簧质量阻尼系统很少以“纯模型”的形式出现，而是作为大量真实装置的等效描述。汽车悬架把车身视为质量、螺旋弹簧视为刚度、减振器视为阻尼，路面激励视为输入；精密仪器的隔振平台、建筑结构的调谐质量阻尼器（TMD）、电机安装座的橡胶隔振垫、包装缓冲结构，都可以用同一组 m、c、k 参数来近似。仿真这类系统的价值在于：在物理样机加工之前，就能预测不同参数组合下的响应幅值、稳定时间与共振风险，从而指导参数选型。

举一个教学示例（以下数值为合理假设，用于说明方法，非实测）：某设备质量 m = 10 kg，安装于刚度 k = 4000 N/m 的隔振垫上，阻尼系数 c = 40 N·s/m。可先手算特征量：固有角频率 ω_n = √(k/m) = √400 = 20 rad/s，对应固有频率约 3.18 Hz；阻尼比 ζ = c/(2√(km)) = 40/(2×200) = 0.1，属于欠阻尼。若初始位移 0.01 m、初速为零，则响应为衰减振荡，包络按 e^−ζω_nt = e^−2t 衰减，振荡频率 ω_d = ω_n√(1−ζ²) ≈ 19.9 rad/s。这些手算结果正是后续仿真曲线应当对照的“基准答案”。

课程学习路径

本课程按“建模—求解—分析—设计”的主线展开，共 30 节。前期建立数学模型与参数体系，明确 m、c、k 与 ω_n、ζ 的换算关系；中期分别用解析法、数值积分法（如四阶龙格–库塔）和频域法求解响应，并对比自由振动、简谐激励、阶跃激励与任意激励下的差异；后期进入参数影响分析、共振识别、隔振设计与仿真结果校核。学习路径不是线性背诵，而是反复在“公式—曲线—物理直觉”之间往返。

仿真任务总览

本课程的仿真任务可归纳为四类，它们共享同一套参数与方程，但观察重点不同：

• 自由振动仿真：给定初始位移或初速，观察欠阻尼、临界阻尼、过阻尼三种情形下位移曲线的形态差异，验证衰减包络与对数衰减率。
• 简谐激励响应：施加 F(t) = F₀sin(ωt)，扫描激励频率，得到幅频曲线，定位共振峰并观察相位滞后。
• 阶跃与脉冲响应：模拟突加负载或冲击，关注最大超调量、上升时间与稳定时间。
• 参数灵敏度分析：固定激励，改变 c 或 k，观察响应如何迁移，为隔振与减振设计提供依据。

下图给出本课程的整体任务结构与各任务之间的输入输出关系，可作为后续每节仿真练习的路线参照。

仿真前必须明确的约定

为避免后续各节出现符号混乱，本课程统一约定：位移 x 以静平衡位置为零点，向右（或向上）为正；速度与加速度同向为正；外力 F(t) 与位移正方向一致时取正。阻尼力方向始终与速度相反，故方程中写作 +cx'。若某节需要改用其他坐标系，会在该节开头明确说明并保持图示一致。

需要强调的是，仿真不是“跑出漂亮曲线”就结束。一个合格的仿真结果应当能回答三个问题：曲线形态是否符合该阻尼状态的理论预期？关键指标（峰值、周期、稳定时间）与手算基准是否在合理误差内？参数微小变化时响应是否按物理规律迁移？后续各节将围绕这三问逐步展开，从单自由度自由振动一直走到带控制的隔振系统设计。