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)提取模态参数的物理依据。
五、从理论到实测:一个教学示例
假设某简化两自由度平台模型,质量矩阵 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 附近出现两个峰;由半功率带宽可估得各阶阻尼比,由两峰幅值比可反推振型分量比。这一过程正是实验模态分析的基本流程:测量频响函数 → 识别峰值 → 拟合模态参数 → 与理论模型对比。
在损伤诊断中,模态参数的变化提供了敏感指标:刚度退化(如杆件屈曲、连接松动)通常使固有频率下降;局部损伤会改变振型形状,使某些测点的振型分量发生偏移;阻尼比异常增大可能提示裂纹摩擦或连接界面滑移。由于频响函数是模态参数的显式函数,实测频响与基准频响的偏差可直接映射到模态参数偏差,再结合正交性检验排除测量误差,形成“频响偏差 → 模态参数偏差 → 损伤定位”的诊断链条。
掌握本节四个核心关系——振动方程描述系统动力学、特征值问题给出模态参数、正交性实现解耦、频响函数建立测量与模态参数的解析桥梁——是后续开展海上平台模态识别、传感器布置优化与损伤指标构造的前提。下一阶段将在此基础上讨论模态参数识别算法与实测数据处理方法。