操纵运动船舶的水动力计算研究.doc操纵运动船舶的水动力计算研究.doc

收藏 分享

资源预览需要最新版本的Flash Player支持。
您尚未安装或版本过低,建议您

专业文档,值得下载专业文档,值得珍藏操纵运动船舶的水动力计算研究李冬荔,杨亮,聂武摘要以船舶操纵水动力预报为研究背景,通过对商用计算流体力学软件FLUENT的二次开发,采用其动网格技术以及后处理系统,对大型船舶操纵水动力导数进行了数值计算船体按照斜航、不同舵角、纯横荡和纯首摇等状态做运动,得出随船坐标系下作用于其上的水动力及力矩通过进行基于最小二乘法的曲线拟合,最终求得船舶操纵水动力导数计算结果与势流理论计算结果一致,表明了所提出的计算方法适用于复杂船舶运动的水动力导数计算关键词船舶操纵;水动力导数;FLUENT;动网格中图分类号U66133文献标识码A文章编号10006982200902000805COMPUTATIONALINVESTIGATIONOFHYDRODYNAMICFORCESAROUNDMANEUVERINGSHIPLIDONGLI,YANGLIANG,NIEWUSTATEKEYLABORATORYOFAUTONOMOUSUNDERWATERVEHICLE,HARBINENGINEERINGUNIVERSITY,HARBIN150001,CHINAABSTRACTINORDERTOPREDICTTHESHIPMANEUVERABILITYINVISCOUSFLOWS,HYDRODYNAMICDERIVATIVESRELATEDTOSHIPMANEUVERINGWEREANALYZEDUSINGADYNAMICMESHMETHODANDPOSTPROCESSINGSYSTEM,WHICHWASBASEDONTHECOMMERCIALSOFTWAREFLUENTSHIPMANEUVERINGMOTIONSOFTHEDIFFERENTHEADINGANGLES,RUDDERANGLES,SWAYINGANDYAWINGARECONSIDEREDHYDRODYNAMICFORCESANDMOMENTSACTINGONAMANEUVERINGSHIPAREOBTAINEDINTHEBODYFIXEDCOORDINATESYSTEMTHUSHYDRODYNAMICDERIVATIVESOFTHESHIPMANEUVERINGMOTIONSARECALCULATEDBYLEASTSQUARECURVEFITTINGTHECOMPUTATIONALRESULTSARECONFIRMEDBYTHEDATAOFPOTENTIALTHEORYMETHOD,WHICHSHOWTHATTHISMETHODCANBEUSEDTOCALCULATETHEHYDRODYNAMICDERIVATIVESOFSHIPCOMPLICATEDMOTIONSKEYWORDSSHIPMANEUVERING;HYDRODYNAMICDERIVATIVES;FLUENT;DYNAMICMESH0引言操纵运动船体水动力的数值计算方法可分为两大类,一是势流方法,二是粘性流方法其中,势流方法对水动力预报精度能基本满足工程要求,成为日趋成熟且实用的方法而三维船舶粘性流的数值研究是从20世纪60年代后期才开始的在1996年第21届ITTC会议上1,求解船舶粘性流的RANS方法被认为已经成熟随着现代船舶科技水平的不断提高,对船舶水动力方面的研究也提出了更高要求许多新技术,如网格生成技术、并行计算技术、多重网格加速等,被广泛应用于复杂的流场计算中,形成了国际上船舶操纵水动力预报的热门和前沿课题本文通过FLEUNT的非定常动网格的方法,数值求解了船舶操纵水动力导数此方法可应用于一系列和船舶操纵问题相关的粘性流动与水动力计算,得到了令人满意的结果1船舶操纵运动水动力导数的数值计算模型在普通的长条形流场区域内,船体不前进,水以一定速度流向船体,船体就可以作各种规定的运动以及频率很低的振荡测量船体所受的水动力,可求得计算船舶操纵运动所需的各种加速度导数和速度导数船舶操纵运动是低频运动,要求接近于零频率的水动收稿日期20081216;修回日期20090114基金项目国家自然科学基金资助项目(50879014)9力导数船舶操纵水动力导数计算是根据船舶所做各种运动计算水动力,可以分为以下五种工况(各种运动情况如图1所示,各运动参数都在随船坐标系下获得,计算公式见参考文献2)数值计算船体所受水动力需达到一个稳定周期,此时方可将水动力(随船坐标系下船体所受侧向力Y与绕船中的力矩N)与运动状态的函数关系进行非线性最小二乘法拟合,即可得到各个水动力导数此处稳定的含义为当船体周期运动到同一相位时,相邻两次受力大小相差甚微图1船体所做的各种运动2FLUENT的求解原理及应用应用FLUENT软件研究船体操纵运动的非定常水动力性能是本文的主要工作,这项工作的核心方法是基于有限体积法对雷诺平均纳维-斯托克斯方程(RANS方程)及其补充、修正方程的离散的分区网格法将该方法直接应用到基于VOF模型的PISO算法是数值计算的最主要内容321三维非结构网格中动网格技术的实现数值模拟船舶操纵运动中,所采用的动网格技术为弹簧近似法4该方法是将非结构网格的整个区域看作一张具有弹性的网,每条边都看作是一根独立的具有倔强系数的弹簧,于是节点I,J间的弹簧张力为IJIJJIFKXX(1)式中,KIJ为连结节点的位置矢量认为移动后网格点的位置满足各点在初始状态的受力情况,即IJIJFS(2)式中,SI为初始状态时点I所受的合力,对所有与节点I相连的节点进行求和图2为船体表面网格以及对整个流场取剖面后的网格结构示意图,船体表面由型线图画出,其表面网格按外形布置疏密可有效保证运动物体边界附近区域的网格质量此外,为了有效减少三维动网格计算量,将整个流场沿船体中纵剖面(ZX平面)与标准吃水处的自由表面(YX平面)划分为四个子区域,并且对接近船体的中纵剖面处和自由表面处进行了单独的网格划分与加密,从而有效控制整个区间网格的疏密程度,并能达到计算要求图2网格结构示意图22边界条件的设定计算区域内自由液面以下的入口、出口处采用给定速度分布;自由液面以上均为标准大气压,计算区域的侧面以及底面均为不可穿透的固壁面,船体的运动规律由UDF程序给出该处所采用的计算模式与控制参数主要是针对基于PISO算法的流场数值计算方法考虑到非定常运动,采用一阶隐式格式方法离散时间项,此外,考虑到自由液面的VOF处理方法,压力项的离散格式采用将体积力加权的处理方法,体积分数的离散格式采用几何重构的处理方法,湍流模式采用RNGK模式,并对壁面处使用标准壁面函数在RNGK模型中,通过在大尺度运动和修正后的粘度项体现小尺度的影响,使这些小尺度运动有系统地从控制方程中去除3数值计算结果及分析31数值方法验证为验证文中计算方法的有效性,应用前面介绍的求解RANS方程的数值方法对61长椭球体在垂直平面内做定常斜航运动时的流场和水动力进行计算,计算网格结构图如图3所示并将计算结果与WETZEL,和图361长椭球体网格结构图不同舵角斜航运动纯横荡运动纯艏摇运动组合运动舵角0漂角X0Y0YYYX0X0X0Y0Y0Y0YYYYYYYYYYYYYXXYXXXXXXXYXXX漂角固定XYZXYZ10SIMPSON6的试验数据及RHEE和HINO7,8的计算结果进行比较,图4为长椭球体升力系数CL与力矩系数CM计算结果结果表明,本文的数值计算方法能够反映以湍流分离流为主的复杂粘性流动的主要物理特征图4长椭球体升力系数与力矩系数计算结果比较32水动力导数的计算结果数值模拟采用实船大小,设计水线间长286M,型宽367M,型深272M,设计吃水94M,航速18KN1)变舵角运动计算控制导数,每个舵角下采用的时间步长均为001S,数值计算6480时间步后船体所受侧向力与力矩结果稳定图5为不同舵角下(弧度)船体所受无量纲侧向力Y与力矩N曲线,得到无量纲化的控制导数拟合采用如下MATLAB命令流;;,,2;132234,,;,0,0,,1132234;,,,1TYPPOLYFITTYFUNINLINEXTXTXTXXTXLSQCURVEFITFUNTYYXTXTXTXPLOTTYTY(3)图5不同舵角下侧向力与力矩曲线2)斜航运动计算位置导数的方法同上,图6为不同漂角下船体所受无量纲侧向力Y与力矩N曲线(不同的漂角得到不同的侧向速度),线性导数结果有表一给出,非线性导数VVVY030948,VVVN000131图6不同侧向速度下侧向力与力矩曲线3)纯横荡运动时,振幅A33M,横荡频率取002、004、005、008、01,来流速度V926M/S(18KN),将一个数值周期取1600个时间步进行计算,计算4个周期,发现在第3个周期船体所受平均阻力、侧向力与力矩结果均收敛图7为F004HZ计算收敛后一个周期内的无量纲化的受力曲线图7纯横荡运动时船体受侧向力与力矩曲线纯横荡运动时船体在随体坐标系下的水动力表达式的无量纲化形式为222321111222VVYYYVVVVVLLVLL(4)将运动表达式SINYAT、COSVYAT、2SINVYAT代入式(4),可得22SINCOSSINCOSVVALAYYTYTVVATBT(5)其中,22VALAYV,VABYV,对各横荡频率下的受力曲线进行非线性拟合得到各系数项A、B由于步骤2)斜航试验给出的位置导数较准确,所以只用该方法确定线加速度导数YV,即对公式22VALAYV进行第二次拟合,直线斜率即为VY,同理求得VN值,图8给出VY拟合结果0012001000080006000400020000000200040006000800100102030攻角/006005004003002001000001002003004005力矩系数CM升力系数CLRHEE等力矩系数计算值,CMWETZEL等力矩系数计算值,CM本文力矩系数计算值,CMRHEE等升力系数计算值,CLWETZEL等升力系数计算值,CL本文升力系数计算值,CL001500100005000000050010001506040004舵角/侧向力Y′力矩N′船体侧向力数值解,Y′船体力矩数值解,N′Y0000009500023X00000047X200019X3,Y′0015001000050000000500100015020206Y00000036000105X000002X200075X3,N′00300200100000100205040004侧向速度V′侧向力Y′力矩N′船体侧向力数值解,Y′船体力矩数值解,N′Y000009003752X000007X2005158X3,Y′020205Y000000200052X000025X2000022X3,N′0003000200010000000100020301010300100005000000050时间T/S侧向力Y′力矩N′船体侧向力数值解,Y′船体力矩数值解,N′Y003192SINT0002475COST,Y′20Y000186SINT0000525COST,N′100010000500000005Y′N′11图8纯横荡运动水动力导数VY拟合结果4)纯首摇运动时,数值求解过程同步骤3)图9为频率为004HZ时计算稳定后一个周期内受力曲线,图10为加速度导数拟合结果,数值表一给出,非线性导数RRRY452E4,RRRN906E5图9纯首摇运动时船体受力曲线图10纯首摇水动力导数RY、RN拟合结果为验证本文计算结果的准确性,另对上述中的线性水动力导数采用势流理论研究,其中,计算惯性力采用无升力势流理论及两因次切片法,计算粘性力采用细长体理论、短翼理论和分离绕流理论,并且将裸船体和附体所受到的水动力分开来计算,再进行叠加,具体计算方法参见文献9,计算结果如表1所示,二者具有明显的一致性表1本文数值方法和势流方法计算值比较水动力导数104本文数值计算值势流理论计算值Y23266N105124VY375380VN5255VY2622VN09712RY143146RN154127RY3228RN2442295)漂角加首摇组合试验确定交叉耦合导数,由于常数和2220SINRRT在一个周期内区间变号积分之值为零以0SINRRT,20ARV代入水动力表达式无量纲化,并在一个周期内进行变号积分得20236322620236322611,,44111SIN29611,,44111SIN296RVVRRRRRVVRRRRYVRDTYVRDTALALYYYVVNVRDTNVRDTALALNNNVV(6)取频率为002HZ,分别计算了漂角为02、04、06、08、12、14、16、18几种情况的组合运动下的水动力,采用离散化函数积分形式计算等式左边的积分值,即一个周期1600步,步长003125S,得到的一个周期内每个步长对应的水动力值计算区间变号积分,带入已求得的RY,RRRY,可以得到VVRY的值同理求出力矩的交叉耦合导数VVRN扣除斜航运动求得的YV、YVVV,即可求得YVRR表2交叉耦合导数求解结果水动力导数VVRYVVRNVRYVRN本文数值计算值00285001008002770030844结论基于CFD软件FLUENT二次开发,采用RNGK模型,数值计算了实船尺度下船体作操纵运动的水动力导数,计算结果表明1)非结构网格的变形移动,该方法可以保证在物面附近的网格质量,同时适用船体斜航、不同舵角、纯横移和纯摇首等状态大幅度运动0015000500005AL2/V2侧向力系数项AY′数值解2Y00022100026X30010410000300020000000230时间T/S侧向力Y′力矩N′船体侧向力数值解,Y′船体力矩数值解,N′Y00018SINT000144COST0000064SIN3T,Y′50Y000068SINT00014COST00006SIN3T,N′400006000200000004Y′N′00010001000200040015000530A3L2/V3侧向力系数项BY′力矩系数B项N′侧向力系数,BY′力矩系数,BN′Y000032X,BY′50Y0000244X,BN′4000100000001500050010000020100
编号:201311161154464431    类型:共享资源    大小:504.00KB    格式:DOC    上传时间:2013-11-16
  
5
关 键 词:
专业文献 学术论文 精品文档 操纵运动
  人人文库网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
关于本文
本文标题:操纵运动船舶的水动力计算研究.doc
链接地址:http://www.renrendoc.com/p-94431.html

当前资源信息

4.0
 
(2人评价)
浏览:33次
dingyx0101上传于2013-11-16

官方联系方式

客服手机:17625900360   
2:不支持迅雷下载,请使用浏览器下载   
3:不支持QQ浏览器下载,请用其他浏览器   
4:下载后的文档和图纸-无水印   
5:文档经过压缩,下载后原文更清晰   

精品推荐

相关阅读

人人文库
关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服客服 - 联系我们

网站客服QQ:2846424093    人人文库上传用户QQ群:460291265   

[email protected] 2016-2018  renrendoc.com 网站版权所有   南天在线技术支持

经营许可证编号:苏ICP备12009002号-5