3. 阻尼类型与数学模型

在阻尼振动系统的响应计算中，“阻尼”并不是一个单一的量，而是一类耗能机制的统称。不同物理机制对应不同的数学描述，进而改变运动方程的形式、求解方法和响应特征。工程上最常见的四类描述是：粘性阻尼、库仑阻尼、结构阻尼和比例阻尼。理解它们的数学表达与适用场景，是正确建立模型、选择求解策略并解释计算结果的前提。

一、粘性阻尼：与速度成正比的线性模型

粘性阻尼假设阻尼力与相对速度成正比，方向与速度相反。设位移为 x(t)，则阻尼力为

F_d = -c · ẋ

其中 c 为粘性阻尼系数（N·s/m）。代入单自由度系统得到线性常系数微分方程：

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

两边除以 m，引入固有频率 ω_n=√(k/m) 与阻尼比 ζ=c/(2√(km))，得到标准形式 ẍ + 2ζω_nẋ + ω_n²x = F(t)/m。这一形式的最大优点是线性：自由振动解为指数衰减正弦，ζ<1 时 x(t)=Xe^-ζω_ntsin(ω_dt+φ)，其中 ω_d=ω_n√(1-ζ²)；强迫振动可用频响函数 H(ω)=1/(k-mω²+jcω) 直接计算幅值与相位。粘性阻尼适用于流体粘性摩擦、油压减振器、空气动力阻尼等场景，也是模态分析与频响计算中最常用的等效模型。

二、库仑阻尼：恒定摩擦力的非线性模型

库仑阻尼来自干摩擦，其大小近似恒定，方向始终与运动方向相反：

F_d = -μN · sgn(ẋ)

其中 μ 为摩擦系数，N 为正压力。由于 sgn(ẋ) 的存在，方程 m·ẍ + k·x = -μN·sgn(ẋ) 是分段线性的，不能直接用叠加原理。其自由振动最显著的特征是振幅按等差递减：每个半周期振幅减少 2μN/k，而不是按指数比例衰减。当振幅降到约 μN/k 以下时，摩擦力足以平衡弹性恢复力，系统停在某个非零位置而不再回到原点，形成“死区”。库仑阻尼适用于干摩擦导轨、螺栓连接滑移、制动器接触等场景。工程上常把它等效为一个“等效粘性阻尼比” ζ_eq=4μN/(πkX) 以便粗略估算，但需注意该等效值随振幅 X 变化，只在小振幅范围内近似有效。

三、结构阻尼：与振幅同相、滞后 90° 的复刚度模型

金属材料内部晶格摩擦、复合材料层间摩擦等引起的能量耗散，通常与频率关系不大，而与振幅相关。结构阻尼（又称迟滞阻尼）常用复刚度描述：

F_d = -j·η·k·x        （频域）
m·ẍ + k(1+jη)x = F(t)

其中 η 为损耗因子（loss factor），表示每弧度能量耗散与最大弹性能之比。其等效粘性阻尼比为 ζ=η/2。结构阻尼的核心特点是：每周期耗散能量与频率无关、与振幅平方成正比，因此频响曲线在共振峰附近呈现更“宽而平”的特征。它适用于金属结构、焊接件、复合材料构件等在共振区的响应估算。需要注意，结构阻尼在时域中对应非因果的卷积形式，一般只在频域分析中使用；若必须做时域积分，应将其近似转换为等效粘性阻尼。

四、比例阻尼：面向多自由度系统的工程假设

对于多自由度系统 M·ẍ + C·ẋ + K·x = F，阻尼矩阵 C 一般难以直接测定。比例阻尼（Rayleigh 阻尼）假设 C 是 M 与 K 的线性组合：

C = α·M + β·K

其中 α（质量比例系数，单位 s⁻¹）与 β（刚度比例系数，单位 s）由实验或经验确定。该假设的关键价值在于：使 C 与模态振型关于 M、K 正交，从而把耦合方程组解耦为 n 个独立的单自由度方程，第 i 阶模态阻尼比为 ζ_i=(α+βω_i²)/(2ω_i)。给定两阶实测阻尼比 ζ₁、ζ₂ 及对应频率 ω₁、ω₂，即可反解 α、β：

α = 2ω₁ω₂(ζ₁ω₂ - ζ₂ω₁)/(ω₂² - ω₁²)
β = 2(ζ₂ω₂ - ζ₁ω₁)/(ω₂² - ω₁²)

例如某两自由度结构实测 ω₁=20 rad/s、ζ₁=0.02，ω₂=60 rad/s、ζ₂=0.05（示例数据），代入可得 α≈0.53 s⁻¹、β≈1.33×10⁻³ s。比例阻尼适用于阻尼分布较均匀、模态耦合较弱的钢结构、机架等，是模态叠加法能顺利实施的重要前提；若阻尼明显局部化（如集中安装的液压减振器），比例阻尼假设会失真，需改用非比例阻尼方法。

五、如何选择与配合使用

四类模型并非互斥，而是对应不同层次的建模需求。实际工程中可按以下思路判断：若耗能主要来自流体或专用减振元件，优先用粘性阻尼，直接进入线性频响或时域积分；若存在明显的干摩擦接触面且振幅较大，用库仑阻尼做分段线性分析，并关注振幅等差衰减与死区；若关注金属结构在共振峰附近的频响形状，用结构阻尼的复刚度形式更贴近实测；若系统为多自由度且需要模态叠加，则用比例阻尼构造 C 矩阵，并用两阶实测阻尼比反解 α、β。

需要特别提醒的是，等效转换必须说明前提。例如把库仑阻尼等效为粘性阻尼时，等效阻尼比随振幅变化，只适合在给定振幅附近做粗略估算；把结构阻尼用于时域积分时，需先转换为等效粘性阻尼比 ζ=η/2，否则会引入非因果误差。建模报告中应写明所选阻尼类型、参数来源（实测、经验或假设）以及适用范围，避免把示例参数当作实测结论使用。