会员注册 | 登录 | 微信快捷登录 支付宝快捷登录 QQ登录 微博登录 | 帮助中心 人人文库renrendoc.com美如初恋!
站内搜索 百度文库

热门搜索: 直缝焊接机 矿井提升机 循环球式转向器图纸 机器人手爪发展史 管道机器人dwg 动平衡试验台设计

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

宽屏显示 收藏 分享

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

专业文档,值得下载专业文档,值得珍藏操纵运动船舶的水动力计算研究李冬荔,杨亮,聂武摘要以船舶操纵水动力预报为研究背景,通过对商用计算流体力学软件FLUENT的二次开发,采用其动网格技术以及后处理系统,对大型船舶操纵水动力导数进行了数值计算.船体按照斜航、不同舵角、纯横荡和纯首摇等状态做运动,得出随船坐标系下作用于其上的水动力及力矩.通过进行基于最小二乘法的曲线拟合,最终求得船舶操纵水动力导数.计算结果与势流理论计算结果一致,表明了所提出的计算方法适用于复杂船舶运动的水动力导数计算.关键词船舶操纵水动力导数FLUENT动网格中图分类号U661.33文献标识码A文章编号10006982200902000805ComputationalinvestigationofhydrodynamicforcesaroundmaneuveringshipLIDongli,YANGLiang,NIEWuStateKeyLaboratoryofAutonomousUnderwaterVehicle,HarbinEngineeringUniversity,Harbin150001,ChinaAbstractInordertopredicttheshipmaneuverabilityinviscousflows,hydrodynamicderivativesrelatedtoshipmaneuveringwereanalyzedusingadynamicmeshmethodandpostprocessingsystem,whichwasbasedonthecommercialsoftwareFLUENT.Shipmaneuveringmotionsofthedifferentheadingangles,rudderangles,swayingandyawingareconsidered.Hydrodynamicforcesandmomentsactingonamaneuveringshipareobtainedinthebodyfixedcoordinatesystem.Thushydrodynamicderivativesoftheshipmaneuveringmotionsarecalculatedbyleastsquarecurvefitting.Thecomputationalresultsareconfirmedbythedataofpotentialtheorymethod,whichshowthatthismethodcanbeusedtocalculatethehydrodynamicderivativesofshipcomplicatedmotions.KeywordsshipmaneuveringhydrodynamicderivativesFLUENTdynamicmesh0引言操纵运动船体水动力的数值计算方法可分为两大类,一是势流方法,二是粘性流方法.其中,势流方法对水动力预报精度能基本满足工程要求,成为日趋成熟且实用的方法.而三维船舶粘性流的数值研究是从20世纪60年代后期才开始的.在1996年第21届ITTC会议上1,求解船舶粘性流的RANS方法被认为已经成熟.随着现代船舶科技水平的不断提高,对船舶水动力方面的研究也提出了更高要求.许多新技术,如网格生成技术、并行计算技术、多重网格加速等,被广泛应用于复杂的流场计算中,形成了国际上船舶操纵水动力预报的热门和前沿课题.本文通过FLEUNT的非定常动网格的方法,数值求解了船舶操纵水动力导数.此方法可应用于一系列和船舶操纵问题相关的粘性流动与水动力计算,得到了令人满意的结果.1船舶操纵运动水动力导数的数值计算模型在普通的长条形流场区域内,船体不前进,水以一定速度流向船体,船体就可以作各种规定的运动以及频率很低的振荡.测量船体所受的水动力,可求得计算船舶操纵运动所需的各种加速度导数和速度导数.船舶操纵运动是低频运动,要求接近于零频率的水动收稿日期20081216修回日期20090114基金项目国家自然科学基金资助项目(50879014)9力导数.船舶操纵水动力导数计算是根据船舶所做各种运动计算水动力,可以分为以下五种工况(各种运动情况如图1所示,各运动参数都在随船坐标系下获得,计算公式见参考文献2).数值计算船体所受水动力需达到一个稳定周期,此时方可将水动力(随船坐标系下船体所受侧向力Y与绕船中的力矩N)与运动状态的函数关系进行非线性最小二乘法拟合,即可得到各个水动力导数.此处稳定的含义为当船体周期运动到同一相位时,相邻两次受力大小相差甚微.图1船体所做的各种运动2FLUENT的求解原理及应用应用FLUENT软件研究船体操纵运动的非定常水动力性能是本文的主要工作,这项工作的核心方法是基于有限体积法对雷诺平均纳维-斯托克斯方程(RANS方程)及其补充、修正方程的离散的分区网格法.将该方法直接应用到基于VOF模型的PISO算法是数值计算的最主要内容3.2.1三维非结构网格中动网格技术的实现数值模拟船舶操纵运动中,所采用的动网格技术为弹簧近似法4.该方法是将非结构网格的整个区域看作一张具有弹性的网,每条边都看作是一根独立的具有倔强系数的弹簧,于是节点i,j间的弹簧张力为ijijjiFKxx(1)式中,Kij为连结节点的位置矢量.认为移动后网格点的位置满足各点在初始状态的受力情况,即ijijFS(2)式中,Si为初始状态时点i所受的合力,对所有与节点i相连的节点进行求和.图2为船体表面网格以及对整个流场取剖面后的网格结构示意图,船体表面由型线图画出,其表面网格按外形布置疏密可有效保证运动物体边界附近区域的网格质量.此外,为了有效减少三维动网格计算量,将整个流场沿船体中纵剖面(ZX平面)与标准吃水处的自由表面(YX平面)划分为四个子区域,并且对接近船体的中纵剖面处和自由表面处进行了单独的网格划分与加密,从而有效控制整个区间网格的疏密程度,并能达到计算要求.图2网格结构示意图2.2边界条件的设定计算区域内自由液面以下的入口、出口处采用给定速度分布自由液面以上均为标准大气压,计算区域的侧面以及底面均为不可穿透的固壁面,船体的运动规律由UDF程序给出.该处所采用的计算模式与控制参数主要是针对基于PISO算法的流场数值计算方法.考虑到非定常运动,采用一阶隐式格式方法离散时间项,此外,考虑到自由液面的VOF处理方法,压力项的离散格式采用将体积力加权的处理方法,体积分数的离散格式采用几何重构的处理方法,湍流模式采用RNGk模式,并对壁面处使用标准壁面函数.在RNGk模型中,通过在大尺度运动和修正后的粘度项体现小尺度的影响,使这些小尺度运动有系统地从控制方程中去除.3数值计算结果及分析3.1数值方法验证为验证文中计算方法的有效性,应用前面介绍的求解RANS方程的数值方法对61长椭球体在垂直平面内做定常斜航运动时的流场和水动力进行计算,计算网格结构图如图3所示.并将计算结果与Wetzel,和图361长椭球体网格结构图不同舵角斜航运动纯横荡运动纯艏摇运动组合运动舵角0漂角x0y0yyyx0x0x0y0y0y0yyyyyyyyyyyyyxxyxxxxxxxyxxx漂角固定xyzxyz10Simpson6的试验数据及Rhee和Hino7,8的计算结果进行比较,图4为长椭球体升力系数Cl与力矩系数Cm计算结果.结果表明,本文的数值计算方法能够反映以湍流分离流为主的复杂粘性流动的主要物理特征.图4长椭球体升力系数与力矩系数计算结果比较3.2水动力导数的计算结果数值模拟采用实船大小,设计水线间长286m,型宽36.7m,型深27.2m,设计吃水9.4m,航速18kn.1)变舵角运动计算控制导数,每个舵角下采用的时间步长均为0.01s,数值计算6480时间步后船体所受侧向力与力矩结果稳定.图5为不同舵角下(弧度)船体所受无量纲侧向力Y与力矩N曲线,得到无量纲化的控制导数.拟合采用如下matlab命令流,,21.32234,,,0,0,,1132234,,,1typpolyfittyfuninlinextxtxtxxtxlsqcurvefitfuntyyxtxtxtxplottyty(3)图5不同舵角下侧向力与力矩曲线2)斜航运动计算位置导数的方法同上,图6为不同漂角下船体所受无量纲侧向力Y与力矩N曲线(不同的漂角得到不同的侧向速度),线性导数结果有表一给出,非线性导数vvvY0.30948,vvvN0.00131.图6不同侧向速度下侧向力与力矩曲线3)纯横荡运动时,振幅a3.3m,横荡频率取0.02、0.04、0.05、0.08、0.1,来流速度V9.26m/s(18kn),将一个数值周期取1600个时间步进行计算,计算4个周期,发现在第3个周期船体所受平均阻力、侧向力与力矩结果均收敛.图7为f0.04Hz计算收敛后一个周期内的无量纲化的受力曲线.图7纯横荡运动时船体受侧向力与力矩曲线纯横荡运动时船体在随体坐标系下的水动力表达式的无量纲化形式为222321..111222VVYYYvvVVVLLVLL(4)将运动表达式sinyat、cosvyat、2sinvyat代入式(4),可得22sincossincosVVaLaYYtYtVVAtBt(5)其中,22VaLAYV,VaBYV,对各横荡频率下的受力曲线进行非线性拟合得到各系数项A、B.由于步骤2)斜航试验给出的位置导数较准确,所以只用该方法确定线加速度导数Yv,即对公式22VaLAYV进行第二次拟合,直线斜率即为VY,同理求得VN值,图8给出vY拟合结果.0.0120.0100.0080.0060.0040.0020.0000.0020.0040.0060.0080.0100102030攻角/°0.060.050.040.030.020.010.000.010.020.030.040.05力矩系数Cm升力系数ClRhee等力矩系数计算值,CmWetzel等力矩系数计算值,Cm本文力矩系数计算值,CmRhee等升力系数计算值,ClWetzel等升力系数计算值,Cl本文升力系数计算值,Cl0.0150.0100.0050.0000.0050.0100.0150.60.40.00.4舵角/°侧向力Y′力矩N′船体侧向力数值解,Y′船体力矩数值解,N′y0.00000950.0023x0.0000047x20.0019x3,Y′0.0150.0100.0050.0000.0050.0100.0150.20.20.6y0.00000360.00105x0.00002x20.0075x3,N′0.030.020.010.000.010.020.50.40.00.4侧向速度V′侧向力Y′力矩N′船体侧向力数值解,Y′船体力矩数值解,N′y0.000090.03752x0.00007x20.05158x3,Y′0.20.20.5y0.0000020.0052x0.00025x20.00022x3,N′0.0030.0020.0010.0000.0010.0020.30.10.10.30.0100.0050.0000.0050时间t/s侧向力Y′力矩N′船体侧向力数值解,Y′船体力矩数值解,N′y0.03192sint0.002475cost,Y′20y0.00186sint0.000525cost,N′100.0100.0050.0000.005Y′N′11图8纯横荡运动水动力导数vY拟合结果4)纯首摇运动时,数值求解过程同步骤3).图9为频率为0.04Hz时计算稳定后一个周期内受力曲线,图10为加速度导数拟合结果,数值表一给出,非线性导数rrrY4.52E4,rrrN9.06E5.图9纯首摇运动时船体受力曲线图10纯首摇水动力导数rY、rN拟合结果为验证本文计算结果的准确性,另对上述中的线性水动力导数采用势流理论研究,其中,计算惯性力采用无升力势流理论及两因次切片法,计算粘性力采用细长体理论、短翼理论和分离绕流理论,并且将裸船体和附体所受到的水动力分开来计算,再进行叠加,具体计算方法参见文献9,计算结果如表1所示,二者具有明显的一致性.表1本文数值方法和势流方法计算值比较水动力导数104本文数值计算值势流理论计算值Y23266N10.512.4vY375380vN5255vY2622vN0.971.2rY14.314.6rN15.412.7rY3.22.8rN2.442.295)漂角加首摇组合试验确定交叉耦合导数,由于常数和2220sinrrt在一个周期内区间变号积分之值为零.以0sinrrt,20arV代入水动力表达式无量纲化,并在一个周期内进行变号积分得20236322620236322611,,44111sin29611,,44111sin296rvvrrrrrvvrrrrYvrdtYvrdtaLaLYYYVVNvrdtNvrdtaLaLNNNVV(6)取频率为0.02Hz,分别计算了漂角为0.2°、0.4°、0.6°、0.8°、1.2°、1.4°、1.6°、1.8°几种情况的组合运动下的水动力,采用离散化函数积分形式计算等式左边的积分值,即一个周期1600步,步长0.03125s,得到的一个周期内每个步长对应的水动力值计算区间变号积分,带入已求得的rY,rrrY,可以得到vvrY的值.同理求出力矩的交叉耦合导数vvrN.扣除斜航运动求得的Yv、Yvvv,即可求得Yvrr.表2交叉耦合导数求解结果水动力导数vvrYvvrNvrYvrN本文数值计算值0.02850.010080.02770.030844结论基于CFD软件FLUENT二次开发,采用RNGk模型,数值计算了实船尺度下船体作操纵运动的水动力导数,计算结果表明1)非结构网格的变形移动,该方法可以保证在物面附近的网格质量,同时适用船体斜航、不同舵角、纯横移和纯摇首等状态大幅度运动.0.0150.0050.0005aL2/V2侧向力系数项AY′数值解2y0.002210.0026x30.0104100.0030.0020.0000.00230时间t/s侧向力Y′力矩N′船体侧向力数值解,Y′船体力矩数值解,N′y0.0018sint0.00144cost0.000064sin3t,Y′50y0.00068sint0.0014cost0.0006sin3t,N′400.0060.0020.0000.004Y′N′0.0010.0010.0020.0040.0150.00530a3L2/V3侧向力系数项BY′力矩系数B项N′侧向力系数,BY′力矩系数,BN′y0.00032x,BY′50y0.000244x,BN′400.0100.0000.0150.0050.0100.00020100
编号:201311161154464431    大小:504.00KB    格式:DOC    上传时间:2013-11-16
  【编辑】
5
关 键 词:
专业文献 学术论文 精品文档 操纵运动
温馨提示:
1: 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
2: 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
3.本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
5. 人人文库网仅提供交流平台,并不能对任何下载内容负责。
6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
  人人文库网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
0条评论

还可以输入200字符

暂无评论,赶快抢占沙发吧。

当前资源信息

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

官方联系方式

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

相关资源

相关资源

相关搜索

专业文献   学术论文   精品文档   操纵运动  
关于我们 - 网站声明 - 网站地图 - 友情链接 - 网站客服客服 - 联系我们
copyright@ 2015-2017 人人文库网网站版权所有
苏ICP备12009002号-5