4. 海洋结构动力分析基础

海洋结构在波浪、海流、风与地震等环境载荷作用下,本质上是一个持续受迫振动的动力学系统。要判断平台、导管架或系泊浮体是否安全,仅靠静力分析远远不够——结构的惯性力、阻尼耗能与共振风险往往才是控制工况。本节从最基本的单自由度(SDOF)模型出发,逐步过渡到多自由度(MDOF)系统,说明模态分析如何把复杂结构解耦为若干独立振型,并给出阻尼模型与运动方程的建立方法,为后续波浪载荷下的动力响应计算打基础。

一、单自由度系统:动力分析的最小单元

把结构简化为一个集中质量 m、一根等效刚度 k 的弹簧和一个阻尼器 c,就得到单自由度系统。其运动方程由达朗贝尔原理(惯性力 + 阻尼力 + 弹性恢复力 = 外载荷)建立:

m·ü(t) + c·u̇(t) + k·u(t) = F(t)

其中 u 为位移,ü、u̇ 分别为加速度与速度,F(t) 为外激励。约定:位移、速度、加速度均以静平衡位置为原点、向上为正,全文与图示保持一致。对自由振动(F=0),可解出无阻尼圆频率 ωₙ = √(k/m),固有频率 fₙ = ωₙ/(2π),周期 Tₙ = 1/fₙ。例如某海洋平台简化模型 m = 2.0×10⁶ kg、k = 8.0×10⁷ N/m,则 ωₙ ≈ 6.32 rad/s,fₙ ≈ 1.01 Hz,Tₙ ≈ 0.99 s。若波浪主频接近 1 Hz,就可能出现显著放大,这正是动力分析要提前识别的风险。

阻尼比 ζ = c/(2√(km)) 是描述能量耗散的关键无量纲参数。ζ 很小时,共振区放大系数约为 1/(2ζ);ζ = 0.05 时放大倍数约 10 倍,ζ = 0.02 时可达 25 倍。海洋结构常用阻尼比范围:钢质导管架 0.02~0.05,钢筋混凝土结构 0.05~0.07,系泊浮体因附加水动力阻尼可更高。这些取值属于工程经验范围,实际项目应结合规范与实测确定。

二、多自由度系统与运动方程

真实海洋结构有无数自由度,有限元离散后得到 N 个自由度。此时运动方程写成矩阵形式:

[M]·{ü} + [C]·{u̇} + [K]·{u} = {F(t)}

[M]、[C]、[K] 分别为质量、阻尼、刚度矩阵,{u} 为节点位移向量。质量矩阵体现惯性耦合,刚度矩阵体现构件连接与边界约束,阻尼矩阵最难直接获得,通常采用近似模型。对海洋结构,还要注意附加质量:结构在水中振动时带动周围水体一起运动,等效为附加质量矩阵 [Mₐ],总质量取 [M]+[Mₐ]。附加质量与振动方向、浸没深度和水深有关,是海洋结构区别于陆上结构的重要特征。

三、模态分析:把耦合系统解耦

忽略阻尼和外力,令 {u} = {φ}·sin(ωt),代入得广义特征值问题:

([K] − ω²[M])·{φ} = 0

求解得到 N 个特征值 ωᵢ² 和对应特征向量 {φᵢ},即第 i 阶固有频率与振型。振型具有正交性:{φᵢ}ᵀ[M]{φⱼ} = 0(i≠j),{φᵢ}ᵀ[K]{φⱼ} = 0。利用正交性引入模态坐标 {u} = [Φ]{q},可将耦合方程变换为 N 个独立单自由度方程:

q̈ᵢ + 2ζᵢωᵢ·q̇ᵢ + ωᵢ²·qᵢ = Pᵢ(t)/Mᵢ

其中 Mᵢ = {φᵢ}ᵀ[M]{φᵢ} 为模态质量,Pᵢ = {φᵢ}ᵀ{F} 为模态力。工程上通常只取前若干阶低频模态参与叠加,因为波浪能量集中在低频段,高阶模态贡献小。判断截断是否足够,可看累计模态质量参与系数是否达到 90% 以上。下图展示从物理系统到模态解耦、再到响应叠加的完整逻辑关系。

物理系统(MDOF) [M]{ü}+[C]{u̇}+[K]{u}={F} 质量/刚度/阻尼矩阵 含附加质量 [Mₐ] N 个耦合自由度 特征值求解 模态分析 ([K]−ω²[M]){φ}=0 ω₁,ω₂,…,ωₙ 固有频率 {φ₁},{φ₂},… 振型 正交性:φᵢᵀMφⱼ=0 模态叠加 模态坐标解耦 {u}=[Φ]{q} q̈ᵢ+2ζᵢωᵢq̇ᵢ+ωᵢ²qᵢ=Pᵢ/Mᵢ N 个独立 SDOF 方程 取前 k 阶,质量参与≥90% SDOF 基础 mü+c u̇+ku=F ωₙ=√(k/m) ζ=c/(2√(km)) 共振放大≈1/(2ζ) 扩展 阻尼模型 比例阻尼 [C]=α[M]+β[K] ζᵢ=α/(2ωᵢ)+βωᵢ/2 钢导管架 ζ≈0.02~0.05 水动力附加阻尼另计 提供 ζᵢ 响应叠加与判断 {u(t)}=Σ φᵢ·qᵢ(t) 波浪主频 vs fₙ 判断共振 动力放大系数 DAF 疲劳/强度校核依据 关键约定: 位移以静平衡位置为原点、向上为正;附加质量计入总质量;模态截断以累计质量参与系数≥90%为判据。 示例参数(m=2.0×10⁶ kg,k=8.0×10⁷ N/m)为教学假设,实际工程须按规范与实测取值。

四、阻尼模型:最难确定却不可忽略

阻尼矩阵 [C] 无法像 [M]、[K] 那样由几何与材料直接组装,工程上常用瑞利比例阻尼:

[C] = α[M] + β[K]
ζᵢ = α/(2ωᵢ) + βωᵢ/2

给定两个目标模态的阻尼比 ζ₁、ζ₂ 和对应频率 ω₁、ω₂,即可反解 α、β。例如取 ω₁=4.0 rad/s、ω₂=12.0 rad/s,ζ₁=ζ₂=0.03,解得 α≈0.18、β≈0.00375。这样在 ω₁ 与 ω₂ 之间阻尼比接近目标值,两端则偏离,因此选点应覆盖波浪能量集中的频段。对海洋结构,还需注意辐射阻尼(结构振动向水中辐射波能)和粘性阻尼(流体绕流分离耗能),这两部分往往不能简单并入比例阻尼,需通过水动力系数或附加阻尼矩阵单独考虑。

五、运动方程建立的实操要点

建立海洋结构运动方程时,建议按以下顺序推进:第一,明确坐标系与正方向,位移、转角、载荷符号统一;第二,划分有限元网格,保证质量与刚度分布合理,浸没部分施加附加质量;第三,定义边界条件,桩基可简化为固定端或等效弹簧,系泊系统用非线性弹簧模拟;第四,组装 [M]、[K],求解特征值得到频率与振型;第五,根据目标频段确定 α、β 并组装 [C];第六,施加波浪载荷时程或谱,进行模态叠加或直接积分。完成后应检查:前几阶频率是否与同类平台经验值吻合、振型是否合理(如首阶应为整体弯曲)、质量参与系数是否足够。若首阶频率异常偏高,常见原因是边界约束过刚或附加质量漏加。

模态分析得到的频率与振型是后续波浪载荷动力放大、疲劳寿命评估和振动控制设计的共同基础。建议在建模阶段就把附加质量、阻尼取值和模态截断判据记录成参数表,便于不同工况复用与复核。

综上,单自由度系统给出频率、阻尼比与共振放大的基本概念;多自由度系统通过矩阵方程描述真实结构;模态分析利用正交性把耦合系统解耦为独立模态方程;阻尼模型则用比例阻尼近似能量耗散。四者串联起来,才构成海洋结构动力分析的完整起点。