邬 婷,杨 立,梁家健
(1.中国船舶科学研究中心,江苏 无锡 214082;2.江苏省绿色船舶技术重点试验室,江苏 无锡 214082)
小水线面双体船型(SWATH)是一种高性能船型,有优良的耐波性和快速性,且甲板面积大易于布置[1],非常适于用作客船、科考船等。但是小水线面双体船也存在一些问题,比如回转半径大[2]。导致在实船的设计过程中,有些船东会要求增加舵面积来减小回转半径,但是过大的舵面积可能会影响到船体回转时的横倾角。《2008年国际完整稳性规则》[3]A部分第三章规定客船计算回航时的稳定横倾角应不大于10°。由于船舶实际所能达到的最大横倾角要大于定常阶段的稳定横倾角,而这个最大横倾角才是最危险的状态。在2014年SDC1次会议上,波兰提议用瞬态最大横倾角代替稳定回转横倾角作为衡准检验的参数[4];在SLF55次会议上,RINA提议船舶回转瞬态最大横倾角不应超过15°[5];日本等国也对横倾角有相关的提案[6]。针对这一问题,本文将采用CFD数值评估方法计算满舵回转时舵面积对于稳定横倾角的影响,从而为SWATH船型的舵设计提供一定参考。
本文以一艘排水量2 000 t级的小水线面双体船为研究对象。该小水线面双体船为一艘科考船,主要用于深远海科学调查作业,设有双机、双桨、双舵,上层建筑布置于首部作为船员和科考人员起居处所与服务处所,尾部主甲板为作业处所,主要布置科考设备,机舱布置于尾部。2 000 t级的小水线面双体船主参数可见表1,其船体和舵的三维效果图如图 1所示,船尾螺旋桨敞水曲线如图2所示。
表1 2 000 t级小水线面双体船船体和螺旋桨主要参数Table 1 Main parameters of hull and propeller of 2000t SWATH
图1 螺旋桨敞水曲线图Fig.1 Open water curve of propeller
本文的计算所采用的控制方程包含连续性方程和RANS方程:
为了模拟船体周围复杂的流场,本文采用SSTk-ω湍流模型。
湍流动能(k)及湍流脉动频率(ω)的输运方程为:
本文采用VOF方法来处理自由面。VOF方法是在整个流场中定义一个函数C,在每个网格中,这个函数定义为一种流体(我们称之为目标流体)的体积与网格体积的比值。VOF(volume of fluid)函数满足连续性方程:
计算中采用基于非结构网格的有限体积法离散计算区域和控制方程,动量输运方程的对流项采用三阶MUSCL格式离散。
本文粘流数值计算所采用的计算域如图 2所示,计算域的各边界设置如下。1)入口边界:距船首1倍船长,为速度入口边界条件;2)出口边界:距船尾2倍船长,为压力出口边界条件;3)顶部边界:距甲板1倍船长,为速度入口边界条件;4)底部边界:距船底2倍船长,为速度入口边界条件;5)侧面边界:距舷侧1.5倍船长,为速度入口边界条件;6)船体和舵:为无滑移的壁面边界条件。
图2 计算域和边界条件Fig.2 Computational domain and boundary conditions
全局计算域的网格划分如图 3(a)所示,为了能更加清晰地捕捉到自由面的特征,可以看到在自由面附近的网格有所加密。船体和舵的网格如图3(b)所示,为了捕捉船体尾部周围的复杂流动,在船尾包括舵附近网格均加密。
图3 网格划分Fig.3 Meshing
利用所建立的模型对船舶在回转运动定常阶段的横倾角以及船后舵的水动力性能(舵的侧向力)进行计算。在计算过程中,模型采用缩尺比为1∶8,采用虚拟盘体(Virtual Disk)模拟船后螺旋桨旋转。模型总长为 7.875 m,螺旋桨直径为0.475 m,以设计吃水为计算条件。
在舵面积的选择中经常采用的参数为舵面积比x:
式中:L为船长,m;d为夏季载重线吃水,m;A为舵面积,m2。
舵的总面积比选取3.0%,4.0%共2种工况,舵的剖面选用 NACA0020,实船回转速度选取12 kn。数值计算详细工况和参数见表2。
表2 数值计算详细参数Table 2 Detailed parameters of numerical calculation
模型的计算结果如表3所示。
表3 模型计算结果Table 3 Model calculation results
舵的侧向力换算到实船,得到实船的计算结果如表4所示。
表4 实船计算结果Table 4 Calculation results of real ship
从计算结果可以看出,该船舶在回转定常阶段的横倾角会随着舵面积比的增大而增加。当舵面积比从3.0%增加到4.0%的时候,横倾角增加了3.7°。舵面积比为 4.0%的时候,定常回转阶段的横倾角达到了12.2°,如果作为客船,则横倾角超过了客船所要求的最大稳定横倾角10°。
随着舵面积的增加,舵的侧向力也会随之增大。操右舵时,右舵的侧向力大于左舵。当舵面积占比从3.0%增加到4.0%的时候,右舵的侧向力增加了101%,左舵的侧向力增加了87%,从舵机选型的角度考虑,成本会增加,这对于船东来说,是不利的。
通常,在舵设计初期,可以有以下几种方法来进行舵面积的选取[7]:1)完全针对具体一条船的舵面积尺寸选取的最好方式是参考母型船。这种方法在设计初期有母型船的资料,并且母型船实船操纵性优良的情况下,可以为舵的选取提供很好的参照依据。2)按照各类船舶舵面积比的统计资料来选取。这些统计资料仅对于船型有一个大致的分类,不会针对到船舶具体的主尺度和特征,因此仅可作为大致的参考。3)在没有性能优良母型船可参考情况下,通过船级社[8]公式计算或者查找图谱。这些给出的是舵面积的最小值,对舵面积的上限则没有给出相关的计算公式。
从设计的经验来看,由于SWATH船型存在回转半径大的缺点,有时船东会要求增加舵的面积,而从上文的评估可以看出,舵面积过大可能会导致SWATH船型回转横倾角过大,这对于对回转横倾角有规范要求的客船来说是不允许的,对于科考船来说也会降低舒适性和安全性。同时,增大舵面积也会增加舵机的功率,增加成本。因此在进行实船舵设计时,需要通过综合考虑和评估,对舵的面积进行论证。
本文以某2 000 t级小水线面双体船为研究对象,以CFD计算为基础,采用通过缩尺比为1∶8的水动力模型,利用所建立的模型对船舶在2种舵面积比工况下满舵回转运动定常阶段的横倾角以及船后舵的水动力性能(舵的侧向力)进行计算,并结合船舶设计的相关规范,得出如下结论。回转定常阶段的横倾角会随着舵面积比的增大而增加,目标小水线面双体船的舵面积比增加到4%时,横倾角会超过客船所允许的最大稳定横倾角。舵的侧向力也会随着舵面积的增加而增大,会导致舵机的功率增加,不利于成本控制。在舵的设计中,尤其对于SWATH这样比较特殊,没有太多母型船可供参考的船型,不能简单地通过增加舵面积来增加操纵性,而是需要综合考虑操纵性与安全性,对舵的面积进行论证,以保证实船的顺利完成试航。