四桨推进船舶螺旋桨负荷数值计算分析

孙帅，常欣，叶礼裕，王超，张洪雨

(哈尔滨工程大学 船舶工程学院，黑龙江 哈尔滨 150001)

四桨推进船舶螺旋桨负荷数值计算分析

孙帅，常欣，叶礼裕，王超，张洪雨

(哈尔滨工程大学 船舶工程学院，黑龙江 哈尔滨 150001)

针对四桨船舶内外桨负荷分配不均匀的情况，本文采用混合网格技术，对四桨船舶螺旋桨水动力性能进行了计算。分析了附体以及内外桨相互干扰对螺旋桨负荷的影响，从伴流场的角度分析了内外桨负荷不均匀的原因。研究发现：前后桨之间相互影响，内后桨处于外前桨的尾流中，所以内后桨受影响更大；与螺旋桨旋转方向相反的周向预旋伴流，可有效提高螺旋桨负荷，且伴流速度越大，效果越明显。双臂轴支架位于螺旋桨近前方，对桨前伴流场具有重要影响，合理调整轴支架角度，可以达到减小内外桨负荷差的目的。

计算流体力学； 负荷差； 支架安装角度； 周向伴流场；四桨船舶； 螺旋桨

Abstract：This study investigates the uneven load distribution of the internal and external propellers of a four- propeller ship. Hybrid grid technology was applied to calculate the hydrodynamic performances of the propellers. The factors that influence load imbalance between the internal and external propellers were analyzed in detail. The uneven load distribution was analyzed in the context of the wake- flow field. Results showed that the front and rear propellers are mutually affected, and that the internal rear propeller is more affected than the external rear propeller given that it lies in the wake flow of the external front propeller. In addition, the circumferential prewhirling wake- flow opposite that of the rotation direction of propeller can effectively increase the load of the propeller. These effects are more apparent under high wake- flow speed. The double- arm axial bracket, which lies in the near front of propeller, has an important effect on the wake- flow field in front of the propeller. The load difference between the internal and external propeller can be reduced by appropriately adjusting the angle of the axial bracket.

Keywords：CFD; load difference; installation angle of bracket; circumferential wake- flow field; four- propellers ship; propeller

随着航运的发展，出现了两桨一舵、四桨两舵等推进和操纵系统的船舶。船机桨的匹配越来越重要，对于四桨工作的船舶而言，还需考虑不同桨之间的负荷分配问题，只有降低螺旋桨之间的负荷差，才能实现船机桨的合理匹配，获得高效的推进性能、经济性和较长的使用寿命。

CFD方法在船舶水动力方面应用逐渐成熟[1-5]。Visonneau等研究了全附体的船桨干扰的尺度效应问题[6]。 Han等采用数值方法研究了船桨舵干扰问题[7]，Roberto等依据重叠动网格技术模拟螺旋桨的真实旋转，采用非定常 RANS 求解器研究了船桨舵干扰问题。预报的推力、扭矩及速度场与实验结果符合良好[8]。Kawakami等对比六组分布不同的三桨船模试验结果，得出两只边桨内旋所需要的功率最小[9]。国内，王金宝等以哥德堡 2000 会议提供的商船 KCS 为对象，考察了时间步长、自由液面是否固化、湍流模式、网格数量和计算策略对计算结果的影响，得到了一些有意义的结论[10]。傅慧萍比较了 FLUENT的多参考系方法、混合面方法和滑移网格方法在船桨整体计算中的优缺点，并预报了螺旋桨旋转所引起的船体脉动压力[11]。沈海龙等基于滑移网格技术，采用 DES 湍流模型，计算了孤立船体、孤立螺旋桨和船桨干扰的非定常流场，所得计算结果与实验吻合较好[12]。

从以上研究可以看出，对于船桨干扰的研究多基于单桨船，而对于多桨船的干扰问题并没有考虑船体影响。覃新川等通过面元法分析了四桨两舵船舶桨的布局对螺旋桨水动力性能的影响，但是在均匀来流情况下分析的，并未考虑船体伴流场的影响[13]。王展智等以四桨水面舰船为研究对象，研究了螺旋桨不同的纵向和横向位置分布，以及舵的位置变化对螺旋桨水动力性能的影响，为四桨干扰问题的研究提供了方法，但是并未提出减小内外桨负荷差的可行性措施[14]。本文主要针对四桨船内外桨负荷不均匀的问题，采用数值计算的方法计算不同工况下内外桨负荷变化情况，对比分析计算结果，找出影响内外桨负荷差大小的因素，并通过分析计算结果提出减小内外桨负荷差的措施，为工程应用提供参考。

1 计算前步骤

1.1 四桨负荷计算方法和步骤

为系统的研究内外桨负荷差异的原因，提出相应改进措施，采用横向、纵向对比相结合方式，从简单到复杂的设计方案，计算分四个阶段，具体步骤如下：

1)计算无支架状态下，内外桨的负荷差，分析船型对负荷差的影响；

2)计算无支架状态下仅内桨或外桨工作时，内外桨负荷的变化情况，与1)计算结果对比，研究内外桨之间的相互影响；

3)考虑轴支架影响，同样与1)计算结果进行对比，分析轴支架对负荷差的影响；

4)在步骤3)的基础上，改变轴支架的安装角度，研究轴支架角度对内外桨负荷差的影响，并提出减小内外桨负荷差的措施。

内外桨负荷影响因素计算分析流程如图1所示。

1.2 计算模型

参照美国核动力航母公开的资料[15-16]，哈尔滨工程大学舰船总体课题组从减小船体阻力出发设计了一艘四桨高速船舶方案，最大航速32.5 kn。通过Fortran编程生成船体外形数据点，并将其导入ICEM中进行船体几何模型建立，其主参数如表1所示，船模缩尺比为37.5∶1。

在哈尔滨工程大学船模拖曳水池进行了船模系列静水阻力实验。其中，测得2.73 m/s(对应实船32.5 kn)下的船模阻力为170.56 N。根据工程实践要求，内外桨的参数需要完全相同，如果内外桨参数不同，需要准备四个备用桨(两个左旋桨，两个右旋桨),这样会造成资源空间的严重浪费。本文根据相关资料和以往的设计经验引入伴流分数和总推力减额分数，并基于螺旋桨升力线和升力面方法为该型船设计了一个螺旋桨，螺旋桨的主要参数如表1所示。

图1 螺旋桨负荷计算分析流程图Fig.1 Flow chart of the propeller load calculation

船模桨模参数数值参数数值垂线间长/m8.267直径/m0.16吃水/m0.289叶数5舷宽/m1.027毂径比0.183巡航航速/(m·s-1)2.45盘面比0.9最大航速/(m·s-1)2.73旋向外旋

1.3 计算域划分

本文认为兴波对内外桨的影响是相同的，为提高计算效率，未考虑自由面兴波的影响。计算域如图2所示。

为计算船桨一体螺旋桨性能，船桨一体计算还需单独划分艉部流域和内外桨流域。不同流域间通过interface进行数据传递。在桨轴后给内外螺旋桨分别创建一个圆柱形小域，小域与桨毂同轴，其直径等于1.2倍螺旋桨直径，通过MRF方法模拟螺旋桨旋转运动，内外桨的安装需保证桨叶相位一致，船艉流域如图2(b)所示。

1.4 网格划分

根据本次计算对象的复杂性，采用混合网格划分形式，在艉部复杂区域进行非结构化网格划分，其他大部分规整区域进行结构化网格划分。网格划分的全部工作在前处理软件ICEM中完成。

图2 计算域划分Fig.2 The division of computational domain

图3 计算网格划分Fig.3 The division of computational mesh

船艉非结构网格划分比较容易，对于相对船体尺寸较小的附体要注意进行局部加密，以较好地表达其几何形状，另外，为了控制第一层网格高度，必须生成棱柱网格。进行粘性流场计算时，由于SSTk-ω模型对近壁面网格有一定的要求，一般认为将Y+值控制在60左右是合理的。Y+值对应的第一层网格厚度可以参考文献[17-18]进行估算，经估算本文将第一层棱柱网格高度设置为0.6 mm，以1.1倍倍率增长，共计5层。划分好的非结构化网格如图3所示。整体计算网格数目为半船520万。

1.5 湍流模型和参数设置

计算中介质水为不可压缩流体，热交换很小以至于可以忽略不计，可只对质量守恒方程和动量守恒方程进行求解。详细公式可参考文献[19]。计算中采用的湍流模型为进行螺旋桨水动力性能计算时比较常用的SST模型，该模型有效集成了k-ε和k-ω模型的优点，能够较好地模拟存在流动分离和强逆压梯度的复杂流动问题。

入口采用水流速度入口，出口设置为压力出口，流域上边界和船体中纵剖面定义为对称面，船体表面、侧面和底面均定义为无滑移壁面，螺旋桨旋转运动稳态计算采用MRF方法，压力速度耦合迭代采用SIMPLEC方法。

2 结果分析

2.1 无支架计算结果

根据计算方案，本文首先计算了无支架时内外桨的推力、转矩及负荷差。选择4个航速(编号1、2、3、4)进行计算，对应实船航速分别为18、24、29和32.5 kn。4个航速下计算得到内外桨推力、转矩负荷差如图4所示。计算结果显示，无支架时，不同航速下内桨的推力和转矩均要小于外桨，即内桨负荷小于外桨，负荷差最大为4.2%。

图4 无支架螺旋桨性能Fig.4 The propeller performance without shaft bracket

2.2 内外桨相互影响计算结果

为研究内外桨之间的相互影响，在上一节的基础上，重新建立模型，仅有内/外桨，将另一个桨去掉，重新划分网格，分别计算了内/外桨单独工作时桨的性能，并将其与两桨同时工作时进行了对比分析。计算结果如图5所示。从计算结果可以看出，内桨单独工作时，其转矩相比于两桨同时工作增加了10%以上，即认为外桨工作对内桨性能的影响很大。外桨单独工作时，其转矩相比于两桨同时工作时增加了4%左右，即认为内桨工作对外桨性能的影响较小，或者是认为外前桨性能受内后桨影响较小。从流场的角度分析，内后桨位于外前桨的艉流中，螺旋桨的桨前伴流场受干扰明显，因此在外桨工作时，内桨性能变化明显。外前桨位于内后桨的前方，虽然内后桨工作时存在一定的抽吸作用，对外前桨的艉流场有一定影响，但外前桨的桨前伴流场受影响很小，因此在内后桨工作时，外前桨性能变化相对较小。

图5 内外桨之间相互影响Fig.5 The mutual influence between propellers

2.3 轴支架影响计算结果

轴支架安装于螺旋桨的近前方，对船艉伴流有一定的导流或整流作用，对螺旋桨性能也有较大影响。在3.1节的基础上，计算了不同航速有支架情况下内外桨的性能。4个航速下计算得到内外桨推力、转矩以及负荷差如图6所示。计算结果显示，有支架时，不同航速下内桨的推力和转矩均要大于外桨，即内桨负荷大于外桨。这与无支架时计算结果正好相反，可以认为支架对内外桨负荷差影响很大。接下来，分别对比内外桨转矩在有支架和无支架时的变化情况，研究支架对桨性能的影响。

图7为内外桨负荷在有支架和无支架时的变化情况，计算结果显示：在有支架时，内桨的转矩相比于无支架时增加4%左右；而外桨在有支架的情况下，转矩相比于无支架时减小幅度较小，在2%左右。即在有支架的情况下，内桨转矩增加，外桨转矩减小。

图6 有支架螺旋桨性能Fig.6 The propeller performance with shaft bracket

图7 支架对螺旋桨负荷的影响Fig.7 The influence of shaft bracket on propeller load

综上，计算结果得出了如下结论：无支架时，外桨负荷大于内桨；有支架时，内桨负荷大于外桨，即支架对螺旋桨的负荷有很大影响；外桨工作对内桨的影响更大。产生上述结果的原因都可以归结为桨前伴流场的差异。但是上述结果还不能描述伴流场是如何影响内外桨负荷的，即我们还不能根据上述计算结果提出减小内外桨负荷差的改进措施。接下来的计算将尝试解决这一问题。

2.4 轴支架角度影响计算结果

要减小内外桨的负荷差，可以从控制内外桨负荷的大小进行考虑，即如果能够通过控制支架的角度来增加或是减小螺旋桨的负荷，则可以根据当前内外桨负荷的大小来调节桨前支架，从而达到减小内外桨负荷差的目的。

根据船舶原理[14]相关知识可知，如果支架能够产生一个与螺旋桨旋转方向相反的周向预旋伴流，如图8(a)所示，可以提高支架后方螺旋桨的负荷。假设外桨负荷相比内桨要小，那么将外桨桨前支架作图8(b)所示调整，使之产生与螺旋桨旋转方向相反的周向预旋伴流，可以提高外桨负荷，从而达到减小内外桨负荷差的目的。

图8 周向预旋伴流与螺旋桨转速关系Fig.8 The ircumferential preselected wake and rotation direction of the propeller

2.4.1 支架角度调整方案A

如图8所示，当预选流方向与螺旋桨转速相反时，螺旋桨的推力和转矩都应该增大，因此可以对负荷较小的桨前支架按图9方式进行调整，使之产生与螺旋桨旋转方向相反的预选流，称之为方案A。

图9 支架角度调整方案AFig.9 The scheme A of adjusting shaft bracket

图10为航速29 kn时，不同支架角度内外桨负荷计算结果。其中“内支架-内桨”表示调整内桨桨前支架后内桨性能变化规律，“外支架-内桨”表示调整外桨桨前支架后内桨性能变化。按照方案A所示的调整方式，分别调整了内桨桨前支架和外桨桨前支架，调整角度分别为5°和10°。从图中可以看出，按照方案A所示的调整方式进行支架角度调整，内桨负荷随着内桨桨前支架角度的增大而增大，外桨负荷变化很小；外桨负荷随着外桨桨前支架角度的增大而增大，而内桨负荷变化较小。由图9可知，当支架角度变大时，支架的预选或是导流作用越明显，即周向伴流分量Vt越大，因此支架后方螺旋桨的转矩相应增大。因此，我们可以得出如下改变内外桨负荷的措施：调节负荷较小的螺旋桨桨前支架角度，使其能够产生与螺旋桨旋向相反的预选流，从而起到增大螺旋桨负荷的效果。

2.4.2 支架角度调整方案B

根据上述结论，如果支架角度按图11所示的方式调整，称之为方案B，那么该支架后方的螺旋桨负荷应该相应减小。为验证假设的准确性，将外桨桨前支架按方案B进行了调整，调整角度为10°。

图10 支架角度与螺旋桨性能关系Fig.10 The relationship between bracket angle and propeller performance

图11 支架角度调整方案BFig.11 The scheme B of adjusting shaft bracket

表2为原布置方案和方案B内外桨负荷计算结果，按照方案B方式调节外支架后，外桨负荷减小了2%，这与上述假设结果一致，说明上述结论是正确的。图12为螺旋桨负荷随支架角度变化规律，更直观说明本文前文得出的结论的准确性。但是考虑到船艉流线向船舯偏转，如图11所示，如果支架角度按方案B所示方式调整，支架角度与流线存在较大的水动力攻角，当航速较高时，支架靠近船舯一侧会产生较强的空化现象，从而对后方螺旋桨的性能造成较大影响，因此不建议采用方案B所示的支架调整方式。

表2不同支架角度下螺旋桨性能计算结果

Table2Thecomputationalresultswithdifferentbracketangles

方案T/NQ/(N·m)内桨外桨内桨外桨原方案23.8123.775-0.79-0.779方案B23.8423.15-0.7901-0.764

图12 螺旋桨负荷随支架角度变化规律Fig.12 The variation of propeller load with bracket angle

2.4.3 不同支架角度方案船舶阻力的变化

支架角度改变必然会带来船舶阻力的改变，为综合评判轴支架角度改变带来的影响，本节给出了航速29 kn时，不同外支架角度方案的船体阻力的变化，如表3所示。从表中可以看出，支架安装角度为-10°时，船体阻力最大，支架安装角度为5°时，船体阻力最小。轴支架的安装角度会带来预旋效果，但角度过大时，船体阻力增加明显。因此在确定支架角度时，应综合考虑阻力及螺旋桨性能等因素。由上述计算可知，本算例中5°为较为合理的外支架安装角度。

表3 不同支架角度下船体阻力计算结果

2.5 桨前伴流场分析

根据上述得出的结论：桨前周向伴流对螺旋桨的负荷有较大影响，伴流方向与螺旋桨旋向相反，能有效提高螺旋桨负荷。因此，在本节中截取了不同支架角度下桨前伴流场及周向平均伴流速度，结果如图13所示。

图13 伴流场变化Fig.13 The variation of wake flow

横向比较来看，内外桨桨前伴流分布存在较大差别，这是内外桨所处船艉位置不同造成的，这种差别很难通过简单的支架角度调整等措施消除。纵向比较研究支架角度对伴流场的影响，图13(a)、(c)为原布置方案与内桨桨前支架调整5°后(按方案A调整方式)内桨桨前周向伴流场，从图中可以看出，支架角度调整后，支架后方周向伴流明显增大，即产生了与螺旋桨旋转方向相反的周向预旋伴流；图13(b)、(d)为原布置方案与外桨桨前支架调整5°后(按方案A调整方式)外桨桨前周向伴流场，其变化趋势与内桨基本一致，支架后方产生了与螺旋桨旋转方向相反的周向预旋伴流。表4为不同支架角度时，桨盘面周向平均伴流速度：纵向比较来看，随着支架角度的调整，周向伴流速度方向由正变负，数值上逐渐增大，螺旋桨的负荷也逐渐增大，这与上述得出的结论是相吻合的；横向比较来看，也呈现此规律，即与螺旋桨旋向相反的周向伴流能够提高螺旋桨负荷，且周向伴流越大，效果越明显。

表4桨盘面周向平均伴流速度

Table4Thecircumferentialaveragewakeflowvelocityatpropellerplane

角度/(°)内支架-内桨外支架-外桨Vt/(m·s-1)M/(N·m)Vt/(m·s-1)M/(N·m)00.0118955440.790.0519445580.7795-0.0400828050.8131-0.01247430.80810-0.0775002540.8354-0.0359408590.8276

3 结论

1)内后桨处于外前桨的艉流中，因此外前桨工作时对内后桨的性能影响较大，而内后桨工作对外前桨的影响较小；

2)轴支架对船艉伴流有一定的导流作用，轴支架的存在对螺旋桨负荷有较大的影响；

3)船艉周向伴流对螺旋桨的转矩影响较大，如果支架能够产生一个与螺旋桨旋转方向相反的周向预旋伴流，可以有效提高支架后方螺旋桨的负荷，且周向伴流越大，效果越明显。

上述结论都是通过模型尺度计算得到的，若要应用到实船，还存在尺度效应等问题，因此需进一步研究附体尺度效应的影响。

[1] 王超,何苗,郭春雨，等. 一体化节能推进装置定常水动力性能的数值模拟分析[J]. 哈尔滨工程大学学报, 2013, 34(6): 674-679. WANG Chao，HE Miao，GUO Chunyu,et al．Steady hydrodynamic performance simulation of integrative energy- saving propulsion system[J]．Journal of Harbin Engineering University， 2013， 34(6) :674-679.

[2] HE Miao, WANG Chao, CHANG Xin, et al. Numerical analysis of propeller wake flowfield using hybrid mesh technology [J]. Journal of marine science and application, 2012, 3: 295-300.

[3] SENTHIL P M N, BINOD C. Numerical estimation of shallow water resistance of a river- sea ship using CFD [J]. International journal of computer applications, 2013, 71(5): 33-40.

[4] WU Chengsheng, ZHOU Decai, GAO Lei, et al. CFD computation of ship motions and added resistance for a high speed trimaran in regular head waves [J]. International journal of naval architecture and ocean engineering， 2011, 3: 105-110.

[5] WANG Jinbao, YU Hai, ZHANG Yuefeng,et al. CFD- based method of determining form factor k for different ship types and different drafts[J]. Journal of marine science and application, 2016,15(3)： 236-242.

[6] SIMONSEN C D． Rudder propeller and hull interaction by RANS[D]． Copenhagen: Technical University of Denmark， 2000: 12-19．

[7] HAN K J，LARSSON L，REGNSTRÖM B． A numerical study of hull/propeller/rudder interaction[C]//Proceedings of 27th Symposium on Naval Hydrodynamics． Seoul，Korea， 2008: 147-153．

[8] ROBERTO M，ANDREA D M． Simulation of the viscous flow around a propeller using a dynamic overlapping grid approach[C]//Proceedings of First International Symposium on Marine Propulsors． Trondheim Norway， 2009: 32-49．

[9] KAWAKAMI Y. Effect of propulsion devices on the propulsive performance of super tankers : comparison among single, twin and overlapping screw propellers[Z]. Report of ship research institute, 1977, 14:71-94.

[10] 王金宝，蔡荣泉，冯学梅． 计及自由面兴波和螺旋桨非定常旋转效应的集装箱船舶绕流场计算[J]． 水动力学研究与进展: A 辑， 2007， 22(4) :491-500． WANG Jinbao，CAI Rongquan，FENG Xuemei． Practical re- search on viscous flow around container ship taking into account the effect of free surface and propeller rotation[J]． Journal of hydrodynamics: Ser A， 2007， 22(4): 491-500．

[11]付慧萍．船桨整体及螺旋桨诱导的船体表面脉动压力计算[J]．哈尔滨工程大学学报， 2009， 30(7): 728-734． FU Huiping． Numerical study of methods for predicting integral hull and propeller and propeller- induced hull surface pressure fluctuations [J]． Journal of Harbin Engineering University， 2009， 30(7): 728-734．

[12]沈海龙，苏玉民． 船体黏性非均匀伴流场中螺旋桨非定常水动力性能预报研究[J]． 水动力学研究与进展: A 辑， 2009， 24( 2): 232-241． SHEN Hailong，SU Yumin． Study of propeller unsteady performance prediction in ship viscous non- uniform wake[J]．Journal of hydrodynamics: Ser A， 2009， 24(2): 232-241．

[13]覃新川，黄胜，常欣． 四桨两舵推进系统的水动力干扰研究[J]． 中国造船， 2008， 49(1) :112-116． TAN Xinchuan，HUANG Sheng，CHANG Xin． Research on hydrodynamic interaction between four propellers and two rudders[J]． Shipbuilding of China， 2008， 49(1): 112-116．

[14]王展智, 熊鹰, 齐万江,等. 船后桨的布局对螺旋桨水动力性能的影响[J]. 哈尔滨工程大学学报, 2012, 33(4): 427-431. WANG Zhanzhi, XIONG Ying, QI Wangjiang, et al. The effect of propeller arrangement behind the hull on propeller hydrodynamic performance [J]. Journal of Harbin Engineering University， 2012， 33(4): 427-431．

[15]BERNARD Ireland. The illustrated guide to aircraft carriers of the world [M]. London: Hermes house, 2005: 96-188.

[16]BIRKLER J, MATTOCK M, SCHANK J, et al. The US aircraft carrier industrial base force structure, cost, schedule, and technology issues for CVN 77[R]. Rand national defense research instsantamonica, 1998: 18-26.

[17]胡坤, 李振北.ANASYS ICEM CFD 工程实例详解[M]. 北京: 人民邮电出版社, 2014: 151-155. HU Kun, LI Zhenbei. Detailed explanation of engineering examples of ANASYS ICEM CFD [M]. Beijing: Post & Telecom Press, 2014: 151-155.

[18]王福军. 计算流体动力学分析- CFD 软件原理与应用[M]．北京: 清华大学出版社，2004: 124-126． WANG Fujun．The analysis of computation fluid dynamics- the based theory and application of CFD[M]．Beijing: Tsing- hua University Press, 2004: 124-126．

[19]KARIM M. MASHUD P. Numerical simulation of free surface water wave for the flow around NACA 0015 hydrofoil using the volume of fluid (VOF) method[J]. Ocean engineering, 2014(78): 89-94.

Numericalanalysisofthepropellerloadofafour-propellerpropulsionship

SUN Shuai, CHANG Xin, YE Liyu, WANG Chao, ZHANG Hongyu

(College of Shipbuilding Engineering, Harbin Engineering University, Harbin 150001, China)

10.11990/jheu.201607004

http://www.cnki.net/kcms/detail/23.1390.u.20170427.1328.018.html

U661.73

A

1006- 7043(2017)09- 1351- 08

2016-07-01. < class="emphasis_bold">网络出版日期

日期：2017-04-27.

国家自然科学基金项目(51379040，51309061).

孙帅(1989-)，男，博士研究生； 常欣(1978-)，男，副教授.

常欣，E- mail: changxin@hrbeu.edu.cn.

本文引用格式：孙帅，常欣，叶礼裕，等. 四桨推进船舶螺旋桨负荷数值计算分析[J]. 哈尔滨工程大学学报， 2017， 38(9): 1351-1358.

SUN Shuai, CHANG Xin, YE Liyu, et al. Numerical analysis of the propeller load of a four- propeller propulsion ship[J]. Journal of Harbin Engineering University， 2017， 38(9): 1351-1358.