2. 结构动力学基础回顾

模态试验与有限元模型修正都建立在结构动力学的基本方程之上。如果把这些方程理解清楚，后面看到的频响函数曲线、模态频率、阻尼比和振型，就不再是孤立的“试验结果”，而是同一物理系统在不同表达方式下的呈现。本节从单自由度系统出发，逐步过渡到多自由度系统，并说明固有频率、阻尼比、模态振型、频响函数、脉冲响应函数以及模态叠加原理之间的内在联系。

一、单自由度系统：最简动力学模型

单自由度（SDOF）系统只有一个独立的位移坐标，典型模型是质量块—弹簧—阻尼器系统。设质量块质量为 m，弹簧刚度为 k，粘性阻尼系数为 c，外力为 f(t)，位移为 x(t)。根据牛顿第二定律，运动方程为：

m·x''(t) + c·x'(t) + k·x(t) = f(t)

两边除以 m，可写成标准形式：

x''(t) + 2ζω_n·x'(t) + ω_n²·x(t) = f(t)/m

其中 ω_n = √(k/m) 为无阻尼固有圆频率，单位是 rad/s；ζ = c/(2√(mk)) 为阻尼比，是无量纲量。固有频率 f_n = ω_n/(2π)，单位是 Hz。固有频率由系统的质量和刚度决定，与外部激励无关；阻尼比则决定振动衰减的快慢。当 ζ < 1 时为欠阻尼，系统会以有阻尼固有频率 ω_d = ω_n√(1-ζ²) 振荡衰减；工程结构通常处于欠阻尼状态，ζ 常在 0.1%～10% 之间。

自由振动时，若给质量块一个初始位移 x₀ 后释放，其响应为：

x(t) = x₀·e^(-ζω_n t)·[cos(ω_d t) + (ζ/√(1-ζ²))·sin(ω_d t)]

从这条曲线可以同时读出三个关键信息：振荡周期对应 ω_d，包络衰减速率对应 ζω_n，而初始幅值和相位则与初始条件有关。这也是后面用自由衰减响应识别模态参数的物理基础。

二、多自由度系统：从单坐标到模态坐标

实际结构有无限多个自由度，有限元离散后得到有限个自由度。设系统有 n 个自由度，位移向量为 x，外力向量为 f，则多自由度（MDOF）运动方程为：

M·x''(t) + C·x'(t) + K·x(t) = f(t)

其中 M、C、K 分别为质量、阻尼和刚度矩阵，通常为 n×n 对称矩阵。与单自由度不同，多自由度系统的各坐标之间存在耦合，一个点受力会引起多个点响应。直接求解耦合方程很困难，而模态分析的核心思想就是通过坐标变换把耦合方程解耦。

忽略阻尼时，无阻尼自由振动方程为 M·x'' + K·x = 0。设解为 x = φ·e^(jωt)，代入后得到广义特征值问题：

(K - ω²M)·φ = 0

该方程有非零解的条件是系数矩阵行列式为零，即 det(K - ω²M) = 0。解此特征方程可得到 n 个特征值 ω_r²（r = 1, 2, …, n），对应 n 个固有频率；每个特征值对应一个特征向量 φ_r，即第 r 阶模态振型。振型描述的是结构在该阶固有频率下各自由度位移的相对比例，是一个形状而非绝对幅值。振型具有正交性：关于质量矩阵和刚度矩阵均正交，即 φ_rᵀMφ_s = 0、φ_rᵀKφ_s = 0（r ≠ s）。正是这一性质使方程能够解耦。

三、频响函数与脉冲响应函数

对运动方程做傅里叶变换，令 X(ω) 和 F(ω) 分别为位移和力的频谱，可得：

(-ω²M + jωC + K)·X(ω) = F(ω)

定义位移频响函数矩阵为 H(ω) = (-ω²M + jωC + K)⁻¹，则 X(ω) = H(ω)·F(ω)。矩阵元素 H_ij(ω) 表示在第 j 点施加单位简谐力时，第 i 点的位移响应，是模态试验中最常测量的量。对单自由度系统，频响函数可写成：

H(ω) = 1 / (k - mω² + jωc)

其幅值在 ω ≈ ω_n 附近出现峰值，峰值高度与阻尼成反比，相位在共振区快速变化约 180°。脉冲响应函数 h(t) 则是频响函数的傅里叶逆变换，二者构成一对傅里叶变换对：

H(ω) = ∫ h(t)·e^(-jωt) dt
h(t) = (1/2π)·∫ H(ω)·e^(jωt) dω

物理上，h(t) 是系统受到单位脉冲激励后的响应，表现为衰减振荡；H(ω) 是系统对不同频率简谐激励的稳态响应。二者从时域和频域两个角度描述同一系统，因此模态试验既可以在时域用脉冲响应识别参数，也可以在频域用频响函数识别参数。

四、模态叠加原理

模态叠加原理是连接多自由度系统与单自由度系统的桥梁。由于振型关于质量矩阵和刚度矩阵正交，可将物理坐标 x(t) 表示为各阶振型的线性组合：

x(t) = Σ_r φ_r · q_r(t) = Φ · q(t)

其中 Φ 为振型矩阵，q(t) 为模态坐标向量。将上式代入运动方程，并左乘 Φᵀ，利用正交性可得一组相互独立的单自由度方程：

m_r·q_r''(t) + c_r·q_r'(t) + k_r·q_r(t) = F_r(t)

式中 m_r = φ_rᵀMφ_r 为模态质量，k_r = φ_rᵀKφ_r 为模态刚度，F_r = φ_rᵀf 为模态力。这样，一个 n 自由度耦合系统就被分解为 n 个独立的单自由度系统。频响函数也可写成各阶模态贡献之和：

H(ω) = Σ_r [ φ_r·φ_rᵀ / (m_r(ω_r² - ω² + j2ζ_rω_rω)) ]

这一表达式说明：实测频响函数上的每一个共振峰对应一阶模态，峰的位置给出 ω_r，峰的宽度给出 ζ_r，而各测点在该峰处的幅值和相位关系则给出振型 φ_r。这正是模态参数识别的基本依据。工程中通常只需前若干阶模态即可较好近似结构响应，因为高阶模态贡献较小，且受测量带宽和噪声限制。

五、从方程到试验的对应关系

把上述理论与模态试验对照，可以建立清晰的映射：试验中在结构上布置加速度传感器，用力锤或激振器施加激励，测量得到的是频响函数 H_ij(ω) 或脉冲响应 h_ij(t)；通过曲线拟合识别出各阶 ω_r、ζ_r 和振型 φ_r；有限元模型则通过求解 (K - ω²M)φ = 0 给出计算模态。模型修正的任务，就是调整有限元模型中的质量、刚度、边界条件等参数，使计算模态与试验模态一致。因此，本节这些方程和概念，是后续所有试验设计、参数识别和模型修正工作的共同语言。