3. 随机振动系统建模

随机振动分析的第一步，是把真实结构抽象成能写出运动方程的力学模型。模型一旦确定，后续的激励谱、响应谱、均方值与疲劳损伤估计才有统一的数学基础。本节围绕单自由度（SDOF）与多自由度（MDOF）两条主线，说明运动方程如何建立、频响函数如何推导、模态叠加如何降阶，以及状态空间表示为何适合数值计算与随机响应求解。

一、单自由度系统：从受力到运动方程

先约定坐标与符号：取质量块静平衡位置为位移 x 的原点，向右为正；速度、加速度同向为正。设质量 m、阻尼系数 c、刚度 k，基础或外力激励为 F(t)。由牛顿第二定律，惯性力、阻尼力、弹性力与外力平衡，得到

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

两边除以 m，引入无阻尼固有频率 ω_n = √(k/m) 与阻尼比 ζ = c/(2√(km))，可写成标准形式

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

若激励来自基础运动 y(t)，令相对位移 z = x − y，则方程变为 m·z̈ + c·ż + k·z = −m·ÿ，右端是基础加速度的等效惯性激励。这一形式在车载、机载与地震激励问题中非常常见，务必区分“绝对位移响应”与“相对位移响应”的物理含义。

二、频响函数：随机激励到随机响应的桥梁

对运动方程做傅里叶变换（零初始条件），令 X(ω)、F(ω) 为对应频谱，可得

H(ω) = X(ω)/F(ω) = 1 / (k − mω² + jωc)
     = 1 / [m(ω_n² − ω² + j·2ζω_nω)]

这就是位移频响函数。它的模在 ω ≈ ω_n 处出现峰值，峰值高度约为 1/(2ζk)，阻尼越小峰越尖锐。随机振动中真正起作用的是功率谱关系：若激励为平稳随机过程，其自功率谱为 S_FF(ω)，则响应自功率谱为

S_xx(ω) = |H(ω)|² · S_FF(ω)

响应均方值 σ_x² = ∫S_xx(ω)dω。这说明随机振动分析的核心不是单频放大倍数，而是“频响函数模平方 × 激励谱”在整个频带上的加权积分。若激励是加速度谱 S_aa，则需用加速度到位移的频响 H_xa(ω) = −1/(ω_n² − ω² + j2ζω_nω)，避免量纲与参考量混淆。

三、多自由度系统：矩阵方程与模态解耦

把结构离散为 n 个自由度，位移向量 x ∈ Rⁿ，运动方程为

M·ẍ(t) + C·ẋ(t) + K·x(t) = F(t)

其中 M、C、K 分别为质量、阻尼、刚度矩阵，通常对称正定或半正定。对无阻尼系统求解广义特征值问题 Kφ = ω_n²Mφ，得到固有频率 ω_nr 与模态向量 φ_r。将物理坐标变换到模态坐标 x = Φq，并左乘 Φ^T，若阻尼满足比例阻尼假设 C = αM + βK，则方程解耦为 n 个独立单自由度方程：

q̈_r + 2ζ_r ω_r q̇_r + ω_r² q_r = φ_rᵀ F / m_r
m_r = φ_rᵀ M φ_r   （模态质量）

这一步是随机振动工程计算的关键：不必直接求解耦合的大矩阵方程，而是对每一阶模态分别计算频响与响应谱，再按模态叠加得到物理响应。实际中常只保留前若干阶对响应贡献大的模态，例如激励频带只覆盖前 3～5 阶时，高阶模态贡献可忽略。

四、状态空间表示：面向数值求解的统一形式

把二阶方程降为一阶方程组，令状态向量 z = [x; ẋ]，则

ż(t) = A·z(t) + B·F(t)
y(t) = C·z(t) + D·F(t)

A = [ 0    I ]
    [ -M⁻¹K  -M⁻¹C ]
B = [ 0 ; M⁻¹ ]

状态空间形式便于用龙格–库塔等数值方法做时域仿真，也便于用李雅普诺夫方程求解随机响应协方差。对白噪声激励，稳态协方差 P 满足 AP + PA^T + Q = 0，其中 Q 由激励强度与 B 决定。这样，随机振动的均方响应就可以从代数方程直接得到，而不必逐频积分。

五、建模实践中的关键判断

实际建模时，先判断激励类型与频带。若激励谱在某一阶固有频率附近能量集中，该阶模态必须保留；若激励为宽带白噪声且覆盖多阶模态，则需按模态贡献排序截断。阻尼处理上，若无法获得完整的阻尼矩阵，比例阻尼是最常用的工程近似；对局部阻尼集中或非比例阻尼明显的结构，则需回到状态空间或复模态方法。

以一个三自由度示例说明：设 m₁=m₂=m₃=1 kg，k₁=k₂=k₃=1000 N/m，比例阻尼 ζ=0.02。求解广义特征值可得前三阶固有频率约为 8.0 Hz、22.4 Hz、32.0 Hz（示例值，实际以数值求解为准）。若基础加速度激励为 0.01 g²/Hz 的平直谱，且只覆盖 5～40 Hz，则三阶模态均需保留；若激励带宽仅 5～15 Hz，则主要保留第一阶，第二、三阶贡献可忽略。判断依据是各阶模态的模态参与质量与频响峰值位置是否落在激励频带内。

综上，单自由度模型给出随机振动分析的基本频响关系，多自由度模型通过模态叠加把耦合问题化为若干独立单自由度问题，状态空间表示则为数值求解与协方差分析提供统一框架。三者不是互相替代，而是同一物理问题在不同计算阶段的不同表达。掌握它们之间的转换关系，才能在面对实际随机振动与宽带噪声联合分析时，选对模型、算对谱、判对结果。