非均匀条件下的船舶锚链静力特性快速计算.pdf

非均匀条件下的船舶锚链静力特性快速计算 2010年5月6日 摘要：为探索锚泊船偏荡时船与锚链相互作用问题，开展了非均匀条件下锚链静力特性 的快速计算研究。采用分段外推-校正方法计算了锚链张力、顶端倾角、悬垂线形状及长度。 通过寻找锚链悬垂线最低点位置以及预估顶端倾角的方法来增加求解速度，实现了快速运 算。 关键词：水路运输；锚泊船；锚链；静力；非均匀条件；快速计算 锚链广泛应用于船舶及海上漂浮结构物的锚泊与定位。锚链形式多种多样，受力情况也 非常复杂。对于最简单的情况：均匀锚链、无潮流作用时，锚链线形状及其受力可以根据悬 垂线理论求出解析解。而在有流速或锚链为非均匀的条件下，通常无法得到锚链线形状及其 受力的解析表达式，需采用数值方法求解。 国内关于锚链系泊性能的研究成果很多。潘斌等计算了均匀流情况下的浮标系泊系统静 力。缪泉明等极限海况下浮标运动及锚链受力进行了时域分析。王冬姣进行了无流速情况下 组合锚链的静力分析。张火明等对深水系泊系统静力特性快速性计算问题做了研究，采用预 先算出的若干拟合曲线的方法，来提高后续计算的速度。滕斌等采用Chebyshev多项式来拟 合锚链受力与顶端位置的关系，节约了计算时间。缪国平介绍了柔性缆索水阻力系数。 所谓非均匀条件，指的是流速分布和锚链单位长度重量皆为非均匀的。采用分段外推- 校正法来计算锚泊船锚链在非均匀条件下的静力特性。这一部分工作是为下一步开展三维锚 链系统以及锚泊船偏荡研究而进行的前期准备。由于问题的复杂性，所需计算时间很长，故 进行快速计算研究。模拟计算表明，对于常规的锚泊船锚泊状态，在普通个人电脑上只需1s 即可完成一个工况的计算。 1锚链静力分析与微分方程建立 锚泊船锚泊状态如图1所示。锚与锚链均受到流的作用。其中流对船的作用力不作专门 1 计算，取经验公式估算。假定该作用力是已知的。锚链受拉力作用而存在拉伸变形。 图1锚泊船与锚链悬垂线示意图 为建立锚链悬垂线微分方程，设如图1所示的坐标系xoy，并从中任取一微段锚链，对 其进行受力分析，如图2所示。由此可建立该微段的静力平衡方程： 图2锚链微元受力分析图 (T+dT)cos(dθ)-T-ωds sinθ+G(1+ε)ds=0 (1) (T+dT)sin(dθ)-ωds cosθ-F(1+ε)ds=0 (2) 式(1)和式(2)中，T为锚链微元张力；θ为张力与水平线夹角；ω为每米锚链重量；ds 为微段长度；F与G分别为水流对锚链微段的法向与切向作用力；dT和dθ分别为张力和角 度的变化量；ε为锚链拉伸变形率： ε=T/AE (3) 2 式(3)中，E为锚链弹性模量；A为锚链横截面积。 空气中每米锚链重量可根据锚链直径估算： ω=2．1462*10·d (4)5 2 式(4)中，d为锚链直径，单位为m；ω的单位为N/m。在水中，其值要乘以系数0.8686。 对于锚泊船，一部分锚链在水面之上，另一部分在水面之下，蹴整个出链长度来讲，每米重 量是不均匀的。 水流作用力由下面两式计算： F=1/2ρCd·Vsinθ (5)2 2 n G=1/2ρC·π·d·Vcosθ (6)t 2 2 式(5)和式(6)中，V为作用在锚链微段上的水流速度，由于流速为非均匀，所以在每个 锚链微段上的流速都不相同；ρ为水的密度；C 和C 分别为法向和切向阻力系数。n t 将式(1)、式(2)展开，忽略高阶小量得：