2. 模态分析理论基础

海上平台是典型的复杂多自由度结构:甲板、导管架、桩基、设备模块通过焊接与螺栓连接成整体,任意一点的振动都会通过结构路径传递到其他部位。要理解模态识别与损伤诊断,必须先掌握多自由度系统的振动方程、特征值问题、模态正交性以及频响函数与模态参数之间的解析关系。本节从动力学方程出发,逐步推导到频响函数的模态展开式,为后续实测模态参数提取与损伤指标构造奠定理论基础。

一、多自由度系统振动方程

将海上平台离散为 N 个自由度后,其线性振动方程可写为:

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

其中 M、C、K 分别为质量、阻尼和刚度矩阵,均为 N×N 实对称矩阵;x(t) 为位移向量,f(t) 为外激励向量。对海上平台而言,M 主要来自甲板质量、设备质量与附连水质量,K 由导管架杆件轴向刚度、桩土相互作用刚度贡献,C 则包含结构材料阻尼与流体阻尼。工程中常采用比例阻尼假设 C = αM + βK,使方程可在实模态空间解耦。

需要注意坐标约定:本文统一采用“位移向上为正、加速度与位移同向为正”的约定,激励力方向与位移正方向一致时取正。后续频响函数的相位判断均基于此约定,若实测采用相反符号约定,相位需整体平移 180°。

二、特征值问题与模态参数

忽略阻尼和外激励,得到无阻尼自由振动方程 M·ẍ + K·x = 0。设解为 x(t) = φ·ejωt,代入后得到广义特征值问题:

(K − ω²M)·φ = 0

该方程有非零解的条件是 det(K − ω²M) = 0,展开后得到 N 个特征值 ωr²(r = 1, 2, …, N),对应 N 个固有频率 ωr 及特征向量 φr(模态振型)。对海上平台这类大型结构,低阶模态(通常前 10~30 阶)集中了大部分振动能量,是模态识别的主要目标。特征向量的物理意义是:结构按第 r 阶固有频率振动时,各自由度位移的相对比例关系,其绝对幅值需通过质量归一化确定。

三、模态正交性

实模态振型关于质量矩阵和刚度矩阵正交,这是模态分析最核心的数学性质:

φᵣᵀ·M·φₛ = 0  (r ≠ s),  φᵣᵀ·M·φᵣ = mᵣ
φᵣᵀ·K·φₛ = 0  (r ≠ s),  φᵣᵀ·K·φᵣ = kᵣ = ωᵣ²·mᵣ

正交性的工程价值在于:利用模态坐标变换 x = Φ·q(Φ 为振型矩阵),可将耦合的 N 自由度方程解耦为 N 个独立单自由度方程。这意味着每一阶模态可以单独分析、单独识别,实测中通过单点激励多点响应即可分离各阶模态贡献。对损伤诊断而言,正交性还提供了判断振型质量的手段:若实测振型与理论振型关于质量矩阵的正交性明显破坏,往往提示测点布置或识别结果存在问题。

四、频响函数与模态参数关系

对振动方程作傅里叶变换,得到频域方程 (−ω²M + jωC + K)·X(ω) = F(ω),定义位移频响函数矩阵 H(ω) = (−ω²M + jωC + K)−1。在比例阻尼假设下,利用模态正交性可将 H(ω) 展开为模态叠加形式:

Hᵢⱼ(ω) = Σᵣ [ φᵢᵣ·φⱼᵣ / (mᵣ·(ωᵣ² − ω² + j·2ζᵣ·ωᵣ·ω)) ]

式中 Hij(ω) 表示在 j 点激励、i 点测得的位移响应;φir 为第 r 阶振型在 i 点的分量;mr 为模态质量;ζr 为第 r 阶模态阻尼比。该式揭示了频响函数与模态参数的完整对应关系:峰值位置对应固有频率 ωr,半功率带宽对应阻尼比 ζr,峰值幅值与振型分量乘积 φirφjr/mr 相关。这正是实验模态分析(EMA)和运行模态分析(OMA)提取模态参数的物理依据。

工程提示:海上平台实测中,若某测点恰好位于某阶模态的节点附近(φir ≈ 0),该阶模态在该测点频响中几乎不可见。因此测点布置应避开预期节点,并保证足够数量以完整描述低阶振型。
多自由度系统:从振动方程到频响函数的模态参数链路 时域振动方程 M·ẍ + C·ẋ + K·x = f(t) M:质量(含附连水) K:刚度(含桩土) C:阻尼 傅里叶变换 广义特征值问题 (K − ω²M)·φ = 0 特征值 ωᵣ² → 固有频率 ωᵣ 特征向量 φᵣ → 模态振型 正交性 模态正交性 φᵣᵀMφₛ = 0 (r≠s) φᵣᵀKφₛ = 0 (r≠s) → 方程解耦为 N 个单自由度 频响函数模态展开式 Hᵢⱼ(ω) = Σᵣ [ φᵢᵣ·φⱼᵣ / (mᵣ·(ωᵣ² − ω² + j·2ζᵣ·ωᵣ·ω)) ] i:响应点 j:激励点 r:模态阶次 ωᵣ:固有频率 ζᵣ:阻尼比 mᵣ:模态质量 代入解耦方程 峰值位置 → ωᵣ 共振峰横坐标即固有频率 识别低阶模态的主要依据 半功率带宽 → ζᵣ 幅值降至 0.707 倍处带宽 ζᵣ = Δω / (2ωᵣ) 峰值幅值 → φᵢᵣφⱼᵣ/mᵣ 多点测量可拟合振型 结合归一化得模态质量

五、从理论到实测:一个教学示例

假设某简化两自由度平台模型,质量矩阵 M = diag(1000, 800) kg,刚度矩阵 K = [[2.0×10⁶, −1.0×10⁶], [−1.0×10⁶, 1.5×10⁶]] N/m。求解特征值问题可得两阶固有频率约为 ω₁ ≈ 26.5 rad/s(4.2 Hz)、ω₂ ≈ 62.0 rad/s(9.9 Hz),对应振型 φ₁ ≈ [0.62, 1.00]ᵀ、φ₂ ≈ [−1.00, 0.78]ᵀ(未归一化)。若在自由度 1 处施加 0~15 Hz 扫频激励,在自由度 2 处测量响应,则 H₂₁(ω) 的幅频曲线应在 4.2 Hz 和 9.9 Hz 附近出现两个峰;由半功率带宽可估得各阶阻尼比,由两峰幅值比可反推振型分量比。这一过程正是实验模态分析的基本流程:测量频响函数 → 识别峰值 → 拟合模态参数 → 与理论模型对比。

在损伤诊断中,模态参数的变化提供了敏感指标:刚度退化(如杆件屈曲、连接松动)通常使固有频率下降;局部损伤会改变振型形状,使某些测点的振型分量发生偏移;阻尼比异常增大可能提示裂纹摩擦或连接界面滑移。由于频响函数是模态参数的显式函数,实测频响与基准频响的偏差可直接映射到模态参数偏差,再结合正交性检验排除测量误差,形成“频响偏差 → 模态参数偏差 → 损伤定位”的诊断链条。

注意:上述两自由度算例为教学假设参数,仅用于说明计算流程,不代表任何实际平台结构的实测结果。实际工程中必须基于真实质量、刚度分布及边界条件建模,并通过实测频响函数验证。

掌握本节四个核心关系——振动方程描述系统动力学、特征值问题给出模态参数、正交性实现解耦、频响函数建立测量与模态参数的解析桥梁——是后续开展海上平台模态识别、传感器布置优化与损伤指标构造的前提。下一阶段将在此基础上讨论模态参数识别算法与实测数据处理方法。