船舶操纵运动波浪力计算

船舶操纵运动波浪力计算2.1 不规则波入射力计算模型依据概率统计理论，不规则波的波面可以看作是由一系列具有不同的频率、 波数、波幅、传播方向以及随机分布初相位角的规则波叠加而成在实际应用中 寻求海浪的统计特性，通常采用“波能谱”的概念来描述海浪海浪形成的过程是风把能量传递给水的过程这一过程大致可分为两个阶段 第一阶段为波浪生长阶段，当风最初作用于海面上时，海面开始出现较小的波， 随着时间的增长，风不断地把能量传递给水，波浪越来越大，显然这一阶段海浪 是比较复杂，其统计特性随时间不断变化，这一阶段的海浪描述描述相当复杂 但是，当波浪渐趋稳定时，波的能量达到一定值，其统计特征基本上不随时间变 化，为了这一阶段海浪的数学描述，应用波谱密度函数，从大量观察分析结果表 明海浪以及船舶在波浪中的运动等均属于狭带谱的正态随机过程，因此基于以下 假设：1．波浪为弱平稳的、各态历经的、均值为零的正态(高斯)随机过程2．波谱的密度函数为窄带3．波峰(最大值)为统计上独立的 由波的方向性谱密度，不规则波的波面可用下列随机积分表示来描述：g(g，耳,t) = J2 ncos[k(g cosesin0) -wt + e(w,0)]、'2S (®,0)dwd0 (2-1)严o g2其中， Sg (w,0 )为波谱密度函数，表示了不规则波浪中各种频率波的能量在总能 g量中所占的份量。

仅考虑波沿主浪向运动的情况，并将式(2-1)转化为随船坐标系下表示为：g(x,y,t) = cos[k(xcos p-ysin p) — ① t + e(①)]、.2S (①)d① (2-2)o e g为了方便计算，将波能谱密度函数进行离散，用求和形式代替上式的积分如 下：g (x, y, t) = 2 ,;'2S (①)Arn cos[k (xcos p - y sin p) 一① t + & ] (2-3)7 g i i ei ii=1其中，相位角£可视为均匀分布在(0, 2兀)区间内的随机变量i 由于不规则波可看作是多个规则谐波分量叠加的结果，因而航行于不规则波 浪中的船舶所受到的主干扰力仍然依据傅汝德-克雷洛夫(Froude-Krylov )假设类比规则波主干扰力的推导过程，深水中不规则波浪对船体的主干扰力(力 矩)仍然是对压力差沿船体表面进行的积分，同样将船体简化成箱体，经推广可 得不规则波对船体的主干扰力和力矩的数学模型表达如下：x =—2 4 pg 3 w 0 4 sin zi=1 iY 丄 4 Pg 3 e w 0 4 cos z i i=1 iZ =—2 4pg 3w 04 sinzcoszi=1 iK 一区他3w 04 coszi=l i、0 2 dig、0.2 dig02Lcosz 02Bsinzsin t sin t sin2 g 2 g02Lcosz 02Bsinzsin t sin — sin2 g 2 g02Lcosz 02Bsinzsin— sin^2g 2g、^2dig02Lcosz 02Bsinzsin— sin^2g 2g 0 2V cos z 、0 i kk i g 丿 0 2V cos z 、0 i kk i g 丿cos 0ikksin 0ikk0 2V cosz、> —i 'g丿0 2V cos z、 ―i g丿、t — s zi丿 02Bsinz 02Bsinz 、2g2 si Bgco 2g 2gE 1—e 02cosz i ik+£ 2Rg2i=lMw 02 sinz i口 iE 1—en =—2n 遊w 02 ii=1 iE 1—e02Lcoszsin—O2d-ig2gsin 0ikk02V cos^、、t —s 丿丿'丿04sin2zi02sinzi02Bsinzsin-^2gsin 002Vcosz、—i i g 丿丿、 t —s 丿丿 i丿 02Lcosz2g2 si 2g04 cos zi02Bcosz 、Lgco 2g 02coszi02Bsinzsin~^2gcos 0kk i®2V cosz、、 t —s 丿 「°.02Lcosz 2g 04cos2ziC c • 92g2sin^02Bcosz、Lgco 2g02cosz 丿i2-4)其中，E「卑讦为各离散规则波的单幅值。

为简化问题，通常假定波浪是二因次的，即波浪只沿一个固定方向传播，而 且波峰线是无限长彼此平行的，它与平面行进波不同的是波浪周期、波高是随即 变化的，通常称这类不规则波浪为长峰不规则波线性水波理论已证明，二因次不规则波波浪是由无限多个不同的波幅和波长 的单元规则线性叠加而成（各单元规则波相位是随机的）这样，长峰不规则波 的数学表达式可写成2-5)z (t) = 2 a cos (kx —01 + s )i i 0 i ii=1式中s 为随机相位同理可得不规则波波浪压力的表达式：2-6)P = £ —pga e-kz cos {k [xcos仲 + 卩)—y sin(屮 + p)]+w t + £ }i i ei ii=1式中 a 可根据波浪谱求的： ia = 2 （①丿A®i Y z0 i当船舶在大地坐标系中运动时，随船坐标系也随之运动假设船舶的航向角为屮，即随船坐标系与大地坐标系X方向的夹角当随船坐标系静止，且原点0与大地坐标系原点重合时，两坐标系之间的转化关系为x = x cos 屮一 y sin 屮2-7)< y = x sin屮 + y cos屮z = z0考虑船的航速，则上式可写为x =(x + ut )cos 屮—(y + vt)sin 屮2-8)2-9)< y = (x + ut )sin屮 +(y + vt )cos屮 0z = z0若在初始时，船艏与波浪有夹角 p ，上式可写为x = (x + ut )cos3+P) — (y + vt )sin(V + P)< y = (x + ut)sin(^+P) + (y + vt)cos(y + P) 0z = z0在本文建立的坐标系中，浪向角的定义见如图：-90将上式代入式（2-6）可得波浪中任意一点的压力分布在计算规则波作用于船舶上的波浪干扰力时，应用著名的傅汝德—克雷洛夫(Froude-Krylov )假设：波浪中船舶的存在不影响波浪的压力分布。

作用在船体上的波浪扰动力与力矩即为波浪动压力沿船体湿表面的积分'F = -ff APndS2-10)wM =-ffAP(r x n )dSwS式中，S 船体的湿表面积n——s的单位外法线矢量，方向指向船体外部 r——动压力作用点相对于随船坐标系的位置向量 写成坐标轴上的投影形式wave= ff p n dsI1wave=ff p n dsI2p n dsI32-11)船体面元网格K=ffp (:yn - zn)dswaveI32M=fsfp(zn - xn)dswaveI13Nwave=ffsp(In - yn21)dsZwaves在数值计算时，采用将船舶表面划分成多个网格，分别计算每个网格上的波浪力再求和的方法2.2 不规则波辐射力计算模型2.2.1 横摇水动力及力矩模型对于船舶在波浪中的横摇运动，由于粘性的影响较大，但粘性系数除了试验 之外尚无准确的理论计算方法所以本文采用较简单的非耦合的横摇运动模型K 二—九 p-2K p-A-GM • sinQ (2-12)H 44 p横摇阻尼系数如下：K =卩 I + 九)-A-GM (2-13)p 申 xx 44式中一般卩在0.055〜0.07之间。

Q 此外由于船舶在回转时会由于横向流体动力的作用而产生横倾，加上横倾力 矩后，式(2-12)变为：K 二一九 p - 2 K p-A-GM - sin Q-Y - z (2-14)H 44 p H H式中，z为横向流体动力Y作用点的z向坐标，HH2-15)z为船舶重心距基线的高度，a为横向力作用点高度系数，当2.5 < BT < 8.0时, g可用下式计算：a = 4 - B T + 0.02(B T - 5.35)3 (2-16)2.3 不规则波漂移力和力矩计算模型波浪的二阶漂移力则会改变船舶航行的航向和航迹，由于理论计算波浪漂移 力较复杂，但已知二阶波浪干扰力与波幅的平方成比例，因此，将波浪漂移力的 数学模型表示如下：XwD一一pgLa 2 cos 咒C2 XwDwD=—pgLa 2 sin 咒C2 YwD2-17)NwD=—pgL a 2 sin 咒C2 NwD其中， C 、C 、C 是关于波浪波长与船长比的系数， Daidola 根据 EnglishXwD YwD NwD的船模试验结果回归得到f 九 九 九C = 0.05 - 0.2(—) + 0.75(—)2 - 0.51(—)3xwd L L L丸 丸 丸< C = 0.46 + 6.83(—)-15.65(—)2 + 8.44(—)3 (2-18)YwD L L L九 九 九C =-0.11 + 0.68(—) - 0.79(—)2 + 0.21(—)3NwD L L L不规则波浪漂移力可以看成各种频率的规则波浪漂移力的叠加，由式（2-17） 可推导不规则波浪漂移力：X = pgL cos 咒》CwD XwDi=1Y = pgL sin 咒工 CwD YwDi=1N = pgL2 sin 咒工C wDNwDi=12兀①2] / \ IS边丛①I gL J G ， 叱]s 6皿I gL 丿 G ，/ 2兀①2 ) /、 1 IS边丛①I gL 丿 G i2-19)。