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% 以上。下图展示从物理系统到模态解耦、再到响应叠加的完整逻辑关系。
四、阻尼模型:最难确定却不可忽略
阻尼矩阵 [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];第六,施加波浪载荷时程或谱,进行模态叠加或直接积分。完成后应检查:前几阶频率是否与同类平台经验值吻合、振型是否合理(如首阶应为整体弯曲)、质量参与系数是否足够。若首阶频率异常偏高,常见原因是边界约束过刚或附加质量漏加。
综上,单自由度系统给出频率、阻尼比与共振放大的基本概念;多自由度系统通过矩阵方程描述真实结构;模态分析利用正交性把耦合系统解耦为独立模态方程;阻尼模型则用比例阻尼近似能量耗散。四者串联起来,才构成海洋结构动力分析的完整起点。