混合对转推进系统水动力性征数值研究.pdf
中 国 舰 船 研 究 第14卷扫码阅读全文第14卷 第2期2019年4月中 国 舰 船 研 究Chinese Journal of Ship ResearchVol.14 No.2Apr. 2019期刊网址www.ship-混合对转推进系统水动力性征数值研究曾旺1,李晔1,秦江涛*21上海交通大学海洋工程国家重点实验室,上海2002402武汉理工大学高性能舰船技术教育部重点实验室,湖北武汉430063摘 要[目的]为研究混合对转推进系统的推进性能,[方法]采用滑移网格方法,分别对螺旋桨、吊舱推进器、混合对转推进系统的水动力性征进行数值模拟,分析混合对转推进系统中各推进成分相对于其单独工作时的水动力性征差异,研究混合对转推进系统中各部分之间的相互影响规律。[结果]结果表明在混合对转推进系统中,后桨对前桨的抽吸作用导致前桨推力系数、转矩系数、敞水效率随进速系数的减小而减小,与单独螺旋桨的变化趋势不同;吊舱和桨毂的阻塞效应以及前桨非定常尾流对后桨的干扰作用,导致后桨推力系数、转矩系数及敞水效率与单独吊舱推进器的结果不同;混合对转推进系统效率低于单独螺旋桨,未能体现后桨对前桨尾流的吸能作用,因此在混合对转推进系统的设计研究工作中,需要进一步考虑前、后桨几何参数的匹配问题。[结论]研究结果可为混合对转推进系统的优化设计和工程应用提供参考。关键词混合对转推进;对转推进器;吊舱推进器;水动力性征;计算流体力学中图分类号U661.1 文献标志码A DOI10.19693/j.issn.1673-3185. 01022Numerical study on hydrodynamic characteristics ofhybrid contra-rotating propellerZeng Wang1,Li Ye1,Qin Jiangtao*21 State Key Laboratory of Ocean Engineering,Shanghai Jiao Tong University,Shanghai 200240,China2 Key Laboratory of High Perance Ship Technology of Ministry of Education,Wuhan University of Technology,Wuhan 430063,ChinaAbstract[Objectives]In order to study the propulsion perance of the hybrid contra-rotatingpropeller,[s]the hydrodynamic characteristics of the propeller,the podded propeller and thehybrid contra-rotating propeller(CRP) are numerically simulated by using sliding mesh respectively,the differences in the hydrodynamic characteristics of each propulsion component in thehybrid CRP compared to those of the entire hybrid CRP are analyzed in detail,and the interaction betweenvarious parts of the hybrid CRP is studied.[Results]The results show that the thrusting coefficient,thetorque coefficient and the open water efficiency of the forward propeller decrease as the speed coefficientdecreases due to the suction effect of the rear propeller on the forward propeller in the hybrid CRP,inwhich the trend is different from that in the single propeller. The thrusting coefficient,the torquecoefficient and the open water efficiency of the rear propeller are different from those of the single poddedpropeller due to the blocking effect of pod and hub and the interference effect of unsteady wake of theforward propeller on the rear propeller. The calculated efficiency of the hybrid CRP is lower than that of thesingle propeller,and the energy absorption effect of the rear propeller on the forward propeller wake is notreflected. Therefore,the matching of the geometry parameters of the forward and rear propellers should befurther considered in the design and study of the hybrid CRP.[Conclusions]This study results canprovide reference for the optimal design and engineering application of the hybrid CRP.Key wordshybrid contra-rotating propeller(CRP);contra-rotating propeller;podded propulsion;hydrodynamic characteristics;computational fluid dynamics收稿日期2017 - 06 - 14 网络首发时间2019-1-18 837基金项目国家自然科学基金资助项目(51479114);海洋能综合配套服务平台建设资助项目(GHME2014ZC01);MIITC浮动支撑平台资助项目(201622)作者简介曾旺,男,1994年生,硕士生。研究方向船舶水动力学。E-李晔,男,1977年生,博士,教授。研究方向流固耦合,水动力学。E-秦江涛,男,1982年生,博士,讲师。研究方向船舶水动力学。E-*通信作者秦江涛网络首发地址http// W,Li Y,Qin J T. Numerical study on hydrodynamic characteristics of hybrid contra-rotating propeller[J].Chinese Journal of Ship Research,2019,14(2)30-35,62.downloaded from www.ship-第2期0 引 言船舶推进器是船舶设计领域的重要研究对象,推进器的种类和性能直接影响船舶的航行性能和营运经济性。随着船舶大型化、高速化和绿色化的发展,传统螺旋桨推进方式在部分船舶中往往不能满足操纵、效率、振动、噪声以及布置等方面的特殊使用要求[1]。部分特种形式的推进器由于在效率和操纵性等方面的优势,受到了业内的广泛关注,如对转螺旋桨和吊舱推进器等。对转螺旋桨的后桨通过回收前桨旋转尾流的部分能量来提高推进系统的效率[2-3];吊舱推进器的吊舱桨由吊舱模块中的驱动电机驱动,且整个吊舱可以全方位转动,从而改善了船舶的操纵性和机动性[4-5]。以对转螺旋桨和吊舱推进器为基础,ABB公司于2000年提出了一种混合对转推进系统(CRP)(图1)。该推进系统的前桨为传统螺旋桨,后桨为可水平全方位旋转的吊舱推进器。CRP综合了对转螺旋桨和吊舱推进器的优点,具有较高的推进效率、良好的操纵性能和空泡性能[6]。国外关于CRP的相关研究较多,日本的三菱重工水池、韩国的三星拖曳水池[7]、荷兰的MARIN水池、芬兰的VTT[8]等机构均针对CRP开展了研究,主要研究了CRP的功率分配比、螺旋桨参数等对推进系统效率与空化性能的影响。国内学者盛立等[9-10]对CRP进行了数值研究和试验研究,对其敞水性征进行了数值预报,并借助空泡水筒开展了相关模型试验与数值结果验证工作。CRP目前已应用于实船,日本三菱重工公司采用CRP的2艘客滚船Hamanasu与Akashia号已经投入营运,并表现出了优良的推进性能。当前针对CRP的研究大多针对给定几何及布置参数的CRP的水动力性征进行数值模拟或试验,而关于CRP中螺旋桨和吊舱推进器各自的推进性征,以及各推进成分之间的相互影响规律等方面的系统研究较少。本文将通过求解单独螺旋桨、单独吊舱推进器和CRP这3种推进器的水动力性能,探讨CRP中螺旋桨和吊舱推进器各自的水动力性能,进一步研究CRP中螺旋桨和吊舱推进器相对于它们单独工作时水动力性能的差异,并将分析前后桨相互作用对各桨推进性能的影响规律,以期为CRP的优化设计和工程应用提供参考。1 数值方法与数值模型1.1 数值模拟对象为了对比分析CRP与单独螺旋桨、单独吊舱推进器的敞水性征差异,分别对3种推进器的敞水性能进行数值模拟,其中螺旋桨和吊舱桨的主要参数如表1所示。CRP由上述螺旋桨与吊舱推进器组合,即CRP的前桨与单独螺旋桨一致,后桨采用吊舱桨,前、后桨间距为0.4倍前桨直径,3种推进器桨叶均采用MAU图谱桨型。3种推进器的几何模型如图2所示。(a)对转螺旋桨 (b)吊舱推进器(c)CRP图1 特种形式推进器Fig.1 Special type propellers(a)单独螺旋桨 (b)吊舱推进器图2 推进器几何模型Fig.2 3D-model of propellers(c)CRP表1 螺旋桨和吊舱桨主要参数Table 1 Main parameters of propeller and pod桨型螺旋桨吊舱桨桨叶45直径/m0.250.22毂径比0.180.24侧斜角/(°)1010螺距比0.801.04盘面比0.5500.614旋转方向右左转速/(r·min-1)600800曾旺等混合对转推进系统水动力性征数值研究 31downloaded from www.ship-中 国 舰 船 研 究 第14卷1.2 流动控制方程以推进器的粘性湍流场为数值模拟对象,数值模型的控制方程为RANS方程¶ρ¶t¶ ρui¶xi0 (1)¶ ρui¶t¶ ρuiuj¶xj-¶p¶xi¶¶xjéëêêùûúúμtæèçöø÷¶ui¶xj¶uj¶xi-23δij¶ui¶xi¶ -ρ- ---ui uj ¶xj(2)式中ρ为流体密度;ui,uj分别为绝对速度矢量的 i 分量和j分量;p 为静压;μt为湍流粘度;-ρ- ---ui uj 为雷诺应力项;当ij时,δij1,当i≠j时,δij0。1.3 计算域、边界条件与离散网格由于存在前、后桨桨叶、舱体、吊舱支架之间的相互干扰,因此采用滑移网格法(Sliding mesh)来处理各部分之间的相对旋转问题。该方法将计算域划分成2个或多个区域,每个区域与其相邻区域存在交界面,通过保证交界面两侧的通量相等来进行数据交换,从而模拟流场中相对旋转物体间强烈的相互作用,该方法具有捕捉流场非定常特性的能力[11]。本文采用滑移网格方法,分别对单独螺旋桨、单独吊舱推进器和CRP的敞水性征进行了数值模拟。单独吊舱推进器的流场计算域分为2个部分吊舱桨旋转区域,以及包含舱体与支架的静止区域。CRP的流场计算域分为3个部分前桨右旋区域、后桨左旋区域,以及包含舱体与支架的静止区域。数值模型边界条件中,桨叶、螺旋桨轴与桨毂、舱体与支架均为无滑移壁面,入口为速度入口,出口为压力出口,存在相对滑移运动的界面设置为交界面。采用切割体网格(Cutcell)对计算域进行离散,其中在桨叶、舱体、支架等近壁区域采用棱柱层网格,以满足近壁流场捕捉的要求;通过面控制和特征线控制,对曲率较大的螺旋面附近的网格进行加密,通过体控制,对近壁面流域进行局部加密,以提高流场捕捉精度。单独螺旋桨计算域网格总数约为110万,单独吊舱推进器计算域网格总数约为132万,CRP计算域网格总数约为178万,各推进器计算域典型剖面的网格分布如图3所示。1.4 湍流模型与近壁处理方法选用Realizable k- ε 湍流模型确保方程组封闭。该模型由标准k- ε 模型改进而来,湍流黏度系数 μt与应变率相关,避免在时均应变率较大时产生负的正应力。Realizable k- ε 湍流模型适用于Re数较大的湍流核心区,在黏性力占主导地位的近壁区,湍流的脉动影响不如分子黏性的影响大,需要对湍流模型进行修正。采用壁面函数法对近壁面区域进行处理,网格布置上需将距壁面法向的第1个节点布置于对数律区域。2 数值计算结果与分析2.1 水动力性能数据表达数值计算结果以无因次化参数表述,其中进速系数J统一表述为(a)螺旋桨横向切面(b)吊舱桨纵向切面(c)CRP水平切面图3 典型剖面的离散网格Fig.3 Discrete grids of propeller,pod and CRP32downloaded from www.ship-第2期JVAnADA(3)式中VA为来流速度(进速);nA为吊舱桨转速;DA为吊舱桨直径。单独螺旋桨与吊舱桨的推进系数采用常规表达方式,此处不再赘述。CRP的敞水效率η0C为η0C(TFTA-R)VA2πnFQF2πnAQA(4)式中nF为螺旋桨转速;TF为螺旋桨的推力;QF为螺旋桨的转矩;TA为吊舱桨的推力;QA为吊舱桨的转矩;R为舱体与支架的力。2.2 数值结果将文献[12]的MAU型系列图谱的试验结果与单独螺旋桨敞水性征数值结果和试验结果进行对比,结果如图4所示,吊舱推进器敞水性征数值结果如图5所示。图中KTF,KTA分别为单独螺旋桨、吊舱桨的推力系数;KQF,KQA分别为单独螺旋桨、吊舱桨的转矩系数;η0F,η0A分别为单独螺旋桨、吊舱桨的敞水效率;KT,KQ,η0为文献[12]的试验结果。由图 4 可以看出,当进速系数 J0.7 时,KTF,KQF,η0F的数值结果与试验结果间误差在7.35以内;当 J0.7时,误差变大,最大误差为9.89。经过分析,在数值计算时,桨叶推力的剪切力部分直接由各个单元分量叠加而得,压力项部分由叶面、叶背压力作差求得。随着进速系数的增大,压力项逐渐占主导地位,压力项部分两较大值作差时会产生较大的相对误差。该较大的进速系数通常偏离了螺旋桨设计工况,计算结果基本满足工程应用要求。CRP的敞水性征数值结果如图6所示。图中KTFC,KTAC分别为CRP前、后桨推力系数,KQFC,KQAC分别为 CRP 前、后桨转矩系数,η0FC,η0AC分别为CRP前、后桨敞水效率。由图6可以看出,CRP的 KTFC,KQFC,η0FC均随J的增大而缓慢增大;CRP的 KTAC,KQAC随J的增大而减小,敞水效率在J0.65附近达到峰值;在不同的J值下,CRP的总效率均大于后桨效率而小于前桨效率,在J0.734附近达到峰值。2.3 水动力性能对比为了对比CRP中螺旋桨和吊舱推进器相对于单独工作时的水动力性征差异,分别绘制了不同推进器推力系数、转矩系数和敞水效率的对比图(图7图9),其中,图7为CRP前桨与单独螺旋桨的对比,图8为CRP后桨与单独吊舱推进器的对比,图9为3种推进器的敞水效率对比。0.2 0.3 0.4 0.5 0.6 0.7 0.8J0.70.60.50.40.30.20.1KT,10KQ,η0KTF10KQFη0FKT10KQη00.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9J0.60.50.40.30.10.1KTA,10KQA,η0AKTA10KQAη0A图5 吊舱推进器敞水性征曲线Fig.5 Open water characteristic curves of pod图4 螺旋桨敞水性征曲线Fig.4 Open water characteristic curves of propeller0.50 0.55 0.60 0.65 0.70 0.75 0.80 0.85 0.90J0.70.60.50.40.30.20.10KT,10KQ,η0KTFC10KQFCη0FCKTAC10KQACη0η0AC图6 CRP敞水性征曲线Fig.6 Open water characteristic curves of contra-rotatingpropeller图7 螺旋桨敞水性能对比Fig.7 Comparison of propeller open water perance0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9J0.80.70.60.50.40.30.20.10KT,10KQ,η0KTFCη0FCKTFη0F10KQF10KQFC曾旺等混合对转推进系统水动力性征数值研究 33downloaded from www.ship-中 国 舰 船 研 究 第14卷由图7可知,CRP的 KTFC,KQFC,η0FC均随进速系数的增大而增大。由于J较小(后桨转速较大)时,后桨对前桨流场产生了较强的抽吸作用而使CRP中前桨效率较低;随着J的增大(后桨转速减小),后桨对前桨流场的抽吸作用逐渐减弱,导致CRP中 KTFC增大,η0FC也随J的增大而增大;和单独敞水桨相比,前桨效率较高的工况对应的J增大,表现出与单独螺旋桨不同的变化趋势。由图8可知,当J在0.5附近时,CRP的 KTAC,KQAC,η0AC基本与单独吊舱推进器相等,随着J的增大,CRP的 KTAC,KQAC,η0AC均小于单独吊舱推进器,且差值越来越大。这是由于随着J的增大(后桨转速减小),KTAC减小,由吊舱和桨毂的阻塞效应产生的阻力在推力成分中的比例增大,且CRP后桨处于前桨尾流中,前桨导致的非定常诱导速度会使 KTAC减小,导致 CRP 的 KTAC,KQAC,η0AC相对于单独吊舱推进器均呈现不同程度的降低。由图9可知,在所计算的工况范围内,η0C始终低于单独螺旋桨;当J0.7时,η0C低于单独吊舱推进器。经分析,认为本文选取的前、后桨参数在CRP中并非较优的参数组合,该CRP未能体现出利用前桨尾流旋转能量提高效率的优势。2.4 流场分析图10为J0.6时,CRP中的前桨与单独螺旋桨叶面(左)、叶背(右)压力系数Cp分布图。图中高压区与低压区的分布规律和变化趋势较相似,但CRP前桨的最大压力系数低于单独螺旋桨,最小压力系数高于单独螺旋桨。分析认为,在该进速系数工况下,CRP后桨对前桨叶面附近流体的抽吸作用较强,导致CRP螺旋桨叶面的最大压力系数低于单独螺旋桨;CRP前方的流体向后流动以补充前桨附近区域,使前桨叶背压力系数有所增大,导致CRP螺旋桨叶背最小压力系数高于单独螺旋桨。图11为J0.6时,CRP后桨与单独吊舱桨叶面、叶背压力系数Cp分布图。由图可以看出,叶面与叶背的压力分布规律仍大致相同,但CRP后桨的最大压力系数和最小压力系数均高于单独螺旋桨。这是由于CRP的后桨处于前桨旋转尾流内且旋转方向相反,前桨的右旋尾流产生向右的周向诱导速度使左旋的吊舱桨表面受压,CRP后桨压力系数整体高于单独吊舱桨。实际上,前桨尾流产生的向右的周向诱导速度会使后桨水动力螺距角减小,攻角增大,轴向推力分量增大,这也是后桨吸收前桨旋转尾流从而提高CRP推进效率的根本原因。同时也可以看出,前桨产生的周向诱导速度对攻角的改变量会受到前桨转速、前桨直径0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9J0.70.60.50.40.30.20.10KTACη0ACKTAη0A10KQA10KQACKT,10KQ,η0图8 吊舱推进器敞水性能对比Fig.8 Comparison of pod open water perance0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9J0.80.70.60.50.40.30.20.1η0Fη0η0Aη0图9 各推进器的敞水效率对比Fig.9 Comparison of open water efficiency of different kinds ofpropellers(b)CRP前桨图10 螺旋桨叶面与叶背的压力系数Cp对比Fig.10 Comparison of propeller on Cp of front and back surface(a)单独螺旋桨-5.386 -3.272 -1.159 0.9548 3.068 5.182Cp-3.766 -2.239 -0.7116 0.8156 2.343 3.870Cp34downloaded from www.ship-第2期(下转第62页)等众多因素的影响,并且后桨攻角的大小和攻角变化量会直接影响后桨的升力系数以及推进效率,这种影响关系并非为正相关。因此,CRP通过前桨尾流能量的部分回收来提高CRP推进效率的目的在实际应用过程中并非一定能达到,需要进一步考虑前、后桨几何参数的匹配问题。对CRP提高效率的机理将在以后的研究中给出。3 结 论通过对CRP水动力性征数值进行研究分析,得出以下结论1)采用Cutcell网格和滑移网格方法所得的螺旋桨推力系数、转矩系数和敞水效率相对于试验结果的误差均在较小的范围内,表明了本文计算方法的合理性。2)CRP后桨对前桨的抽吸作用随着J的增大(后桨转速的减小)而减弱,导致CRP的 KTFC,KQFC,η0FC随J的增大而增大,效率较高的工况对应的J增大,表现出与单独螺旋桨不同的变化趋势;吊舱与桨毂的阻塞效应以及前桨非定常诱导速度对后桨的干扰作用,导致 KTAC,KQAC,η0AC相对于单独吊舱推进器均呈现出不同程度的降低。3)本文所计算的CRP的效率低于单独螺旋桨,未能体现后桨对前桨尾流的能量吸收作用,认为该CRP的参数并非推进性能较优的参数组合,故在CRP的设计研究工作中,需要进一步考虑前后桨几何参数的匹配问题,否则无法发挥CRP的高效节能优势。参考文献[1] 张涛,杨晨俊,宋保维. 基于MRF模型的对转桨敞水性能数值模拟方法探讨[J].船舶力学,2010,14(8)847-853.Zhang T,Yang C J,Song B W. Investigations on thenumerical simulation forthe open-water perfor⁃mance of contra-rotating propellers based on the MRFmodel[J]. Journal of Ship Mechanics,2010,14(8)847-853(in Chinese).[2] RICHARDS J,ENGELSKIRCHEN J,PRAEFKE E.Counter rotating propellers without complex shaftingfor a fast monohull ferry[C]//Proceedings of the 6thIn⁃ternational Conference on Fast Sea Transportation.Southampton,UKSpringer-verlag,2001.[3] 辛公正,刘登成,施小勇,等.对转吊舱推进器水动力性能预报与模型试验[C]//第二十三届全国水动力学研讨会暨第十届全国水动力学学术会议论文集.西安中国力学学会,中国造船工程学会,201112.Xin G Z,Liu D C,Shi X Y. Numerical prediction andmodel test of hydrodynamic perance for contra-ro⁃tating POD propulsion[C]. Proceedings of the 23rd Na⁃tional Conference on Hydrodynamics and 10th Nation⁃al Congress on Hydrodynamics. Xi an Chinese Societyof Mechanical Engineering, China Shipbuilding Engi⁃neering Society,201112.[4] 沈兴荣,蔡跃进,蔡荣泉,等.吊舱式推进器水动力性能研究综述[J]. 船舶力学,2011,15(1/2)189-197.Shen X R,Cai Y J,Cai R Q,et al. Summarization ofhydrodynamics perance study of podded propulsor[J]. Ship Mechanics,2011,15(1/2)189-197(in Chi⁃nese).[5] 范露.吊舱推进器的水动力性能优化设计研究[D].武汉武汉理工大学,2008.Fan L. Study on optimization design of hydrodynamiccharacteristics of pod propeller[D]. Wuhan WuhanUniversity of Technology,2008.[6] 张庆文. 吊舱式CRP推进系统发展及应用前景[J].船海工程,2007,36(2)57-60.Zhang Q W. Prospect and utilization of the poddedCRP propulsion system[J]. Shipbuilding Engineering,2007,36(2)57-60(in Chinese).(a)单独吊舱桨(b)CRP后桨图11 吊舱桨叶面与叶背的压力系数Cp对比Fig.11 Comparison of pod on Cp of front and back surface-9.676 -6.220 -2.763 0.6931 4.149 7.606Cp-7.639 -4.403 -1.168 2.068 5.304 8.539Cp曾旺等混合对转推进系统水动力性征数值研究 35downloaded from www.ship-中 国 舰 船 研 究 第14卷形的区域,其von Mises等效应力随计算点距裂尖距离的增大而减小,两者之间基本上呈线性或双线性关系;在未发生塑性变形的区域,基本上是裂尖距离的三次函数。2)对于裂纹横断面,裂纹总是沿裂纹边缘形状的法向方向扩展,应力最大值总是出现在裂纹边缘。在裂纹内部,应力值随着距裂纹中心的距离而呈现高度的非线性分布;在裂纹外部,应力值随着距裂纹中心距离的增大而减小。当距裂纹中心的距离大于裂纹半径1.67倍时,应力集中消失。3)同样,对于裂纹横断面,当载荷为材料屈服应力的2/3倍时,裂纹边缘最大应力值总是为载荷的1.45倍;当载荷为材料屈服应力的1/3倍时,裂纹边缘最大应力值总是为载荷的2倍。参考文献[1] Cui C,Yan G P,Xia T,et al. Assessment of residualultimate strength of cracked steel plates under longitu⁃dinal compression[J]. Ocean Engineering,2016,121174-183.[2] Bayatfar A,Khedmati M R,Rigo P. Residual ultimatestrength of cracked steel unstiffened and stiffenedplates under longitudinal compression[J]. Thin-Walled Structures,2014,84378-392.[3] 闫小顺,周心桃,黄小平. 球柱结合壳焊趾表面裂纹坐标变换有限元建模方法及应力强度因子计算[J]. 中国舰船研究,2016,11(4)59-66.Yan X S,Zhou X T,Huang X P. Stress intensity factorcalculation for surface cracks at weld toe of sphere-cyl⁃inder shells[J]. Chinese Journal of Ship Research,2016,11(4)59-66(in Chinese).[4] 吴建国,王奇志,张行,等. 基于能量释放率理论的三维应力强度因子显式闭合解[J]. 中国舰船研究,2014,9(1)81-92.Wu J G,Wang Q Z,Zhang X,et al. Explicitclosed- solutions for three-dimensional stress in⁃tensity factors based on the energy release rate theory[J]. Chinese Journal of Ship Research,2014,9(1)81-92(in Chinese).[5] Rybicki E F,Kanninen M F. A finite element calcula⁃tion of stress intensity factors by a modified crack clo⁃sure integral[J]. Engineering Fracture Mechanics,1977,9(4)931-938.[6] Shivakumar K N,Tan P W,Newman J C,Jr. A virtualcrack-closure technique for calculating stress intensityfactors for cracked three dimensional bodies[J]. Inter⁃national Journal of Fracture,1988,36(3)R43-R50.[7] Narayana K B,Dattaguru B,Ramamurthy T S,et al.A general procedure for modified crack closure integralin 3D problems with cracks[J]. Engineering FractureMechanics,1994,48(2)167-176.[8] Xie D,Biggers S B,Jr. Progressive crack growth analy⁃sis using interface element based on the virtual crackclosure technique[J]. Finite Elements in Analysis andDesign,2006,42(1)977-984.(上接第35页)[7] 汪小翔.混合式CRP吊舱推进器的水动力性能研究[D].哈尔滨哈尔滨工程大学,2013.Wang X X. Research on the hydrodynamic perfor⁃mance of hybrid CRP podded propulsion[J]. HarbinHarbin Engineering University,2013(in Chinese).[8] 廖慧清,魏海波,李祝清,等.混合电力吊舱式推进系统[J].船电技术,2006,26(1)1-3.Liao H Q,Wei H B,Li Z Q,et al. Hybrid podded pro⁃pulsion with contra-rotating propellor[J]. Marine Elec⁃tric Electronic Engineering,2006,26(1)1-3(inChinese).[9] 盛立,熊鹰.混合式CRP吊舱推进器水动力性能数值模拟及试验[J]. 南京航空航天大学学报,2012,44(2)184-190.Sheng L,Xiong Y. Numerical simulation and experi⁃mentalinvestigation on hydrodynamics perance ofhybrid CRP podded propulsion[J]. Journal of NanjingUniversity of Aeronautics Astronautics,2012,44(2)184-190(in Chinese).[10] 张可,熊鹰,王展智.前后桨间距对吊舱式CRP水动力性能的影响[C]//中国造船工程学会2013年船舶水动力学学术会议论文集.西安中国造船工程学会,20139.[11] 郭春雨,汪小翔,赵大刚,等.吊舱式CRP推进器的定常水动力性能模拟[J]. 船海工程,2013,42(3)110-114,118.Guo C Y,Wang X X,Zhao D G,et al. Simulation ofthe steady hydrodynamics perance on poddedCRP propulsion[J]. Ship Ocean Engineering,2013,42(3)110-114,118(in Chinese).[12] 盛振邦,刘应中.船舶原理[M].上海上海交通大学出版社,2004264.Sheng Z B,Liu Y Z. Ship principles[M]. ShanghaiShanghai Jiao Tong University Press,2004264.62downloaded from www.ship-