欢迎来到人人文库网! | 帮助中心 人人文库renrendoc.com美如初恋!
人人文库网
首页 人人文库网 > 资源分类 > DOC文档下载

数值波浪水池中船舶水动力系数测试与分析技术.doc

  • 资源大小:405.50KB        全文页数:8页
  • 资源格式: DOC        下载权限:游客/注册会员/VIP会员    下载费用:5
游客快捷下载 游客一键下载
会员登录下载
下载资源需要5

邮箱/手机号:
您支付成功后,系统会自动为您创建此邮箱/手机号的账号,密码跟您输入的邮箱/手机号一致,以方便您下次登录下载和查看订单。注:支付完成后需要自己下载文件,并不会自动发送文件哦!

支付方式: 微信支付    支付宝   
验证码:   换一换

友情提示
2、本站资源不支持迅雷下载,请使用浏览器直接下载(不支持QQ浏览器)
3、本站资源下载后的文档和图纸-无水印,预览文档经过压缩,下载后原文更清晰   

数值波浪水池中船舶水动力系数测试与分析技术.doc

专业文档,值得下载专业文档,值得珍藏数值波浪水池中船舶水动力系数测试与分析技术郭海强,朱仁传,缪国平,余建伟(上海交通大学船舶海洋与建筑工程学院,上海200030)摘要以CFD理论为基础,建立了数值波浪水池,给出了一种基于三维数值波浪水池的船舶水动力系数测试与分析技术。对WIGLEYIII船模在数值波浪水池中受迫振荡进行了数值模拟,计算分析了船体垂荡、纵摇及纵荡运动的附加质量与阻尼,并与三维势流理论计算结果进行了比较,两者吻合良好。此方法能准确给出浮式结构物的水动力系数,细致描述船舶周围的流场,可广泛应用于船舶与海洋工程浮式结构物的水动力性能研究。关键词船舶、舰船工程;数值波浪水池;CFD;水动力系数;船舶运动中图分类号U6617文献标识码A1引言船舶水动力系数如船舶运动时的附加质量与阻尼系数,是准确预报船舶在波浪中运动的基础。有很多船舶在波浪中运动性能预报的理论方法,包括切片理论、细长体理论、三维频域理论及其各种改进方法等等1。这些方法基本上是基于线性势流理论建立起来的。但由于线性势流理论的局限性,很难考虑粘性和非线性的影响,如果不加修正,船舶运动的预报精度较差。因此在工程应用上,一般采用基于实验获得的经验公式或半经验公式对运动计算加以修正2,特别是横摇运动。目前获得船舶附加质量与阻尼的试验方法主要分为两种一种是船舶摇摆的自由衰减试验,但试验结果只能得到在共振频率处的船舶附加质量和阻尼,缺乏频率的相关性;另一种是船舶的强迫振荡运动试验,这虽能给出多个频率下的附加质量、阻尼及水动力,但往往受到试验船模的限制,而且对试验的设备和测试系统要求较高,对此有许多学者仍在致力试验测试系统的研究37。近年来,随着计算机技术和计算技术的飞速发展,计算流体力学(CFD)也有了长足的进步。基于CFD理论的船舶水动力学方面的数值模拟,因为具有费用低、无触点流场测量、无比尺效应、能消除物模中由传感器尺寸及模型变形等因素对流场的影响、可获得较为详细的流场信息等优点而广受关注,且应用的范围越来越广811。本文以CFD理论为基础,建立了数值波浪水池,给出了一种基于三维数值波浪水池的船舶水动力系数的测试与分析技术。文中数值波浪水池的建立,船体在数值波浪水池中受迫振荡的数值模拟实验采用的动网格技术,以及船体运动时的附加质量与阻尼的计算分析,都是基于FLUENT商业软件平台的二次开发实现的。还就WIGLEYIII船模在三维数值波浪水池中受迫振荡进行了数值模拟,船体的附加质量与阻尼系数计算分析结果与势流理论结果进行了比较。对比和研究表明,本文方法能准确给出浮式结构物的水动力系数,能细致描述船舶周围的流场,可广泛用于船舶与海洋工程浮式结构物的水动力性能研究。49卷增刊总第183期郭海强,等数值波浪水池中船舶水动力系数测试与分析技术592数值波浪水池与船体运动模拟船舶水动力系数的测试试验一般是在物理波浪水池中进行的。而笔者的数值波浪水池中船舶水动力系数测试则是一个仿物理试验的数值模拟计算。仿物理试验水池的数值波浪水池也具备了类似于试验中的造波、消波功能(见文献12,鉴于本文未用到水池造波等方面的功能,故对此不予叙述)。数值波浪水池是以多相流理论为基础的,水池形状为一长方体,如图1所示。没有波浪时其上部为空气下部为水,整个流场以连续性方程和NS方程为控制方程01,2,3IIUITX1IJJIIIJJIUUUUUTXXXX1,2,3IIPFIX2式中,IU为流体质点在I方向的速度分量,IF为质量力,P为流体的压力,流体密度定义为21QQQA,其中体积分数AQ表示单元内第Q相流体占的体积与总体积的比例,并且有211QQA,为相体积分数平均的动力粘性系数,与密度定义的形式一致。自由面的波动是采用VOF方法来追踪的,可写为0QIQIAUATX,1,2Q3式中A1、A2分别为空气相、水相的体积分数。船舶水动力系数测试的仿物理实验为数值波浪水池中船模受迫振荡试验,用以测试船体的附加质量与阻尼系数。被测试船模将分别在不同频率下以固定幅值作单一模态的振荡运动。假定船体的某一模态的振荡运动可描述为0SINT4式中为振荡运动的频率,Ξ0和分别为振荡运动的幅值和相位角。由于没有假定运动是微幅的,故式(4)也适用于大振幅的运动,也就是说,测试的结果既可计入粘性影响,又可计入运动、自由面以及船体形状产生的非线性的影响。数值波浪水池中船体受迫振荡运动的模拟采用了移动网格技术见文献13,网格划分均采用结构化网格。如图2所示,船体位于动态网格区域内,动态网格与船体之间相对位置不变,并且整个动态网格区域以给定的振荡运动规律作相应的刚体运动。为了能够准确描述运动船体周围的自由面和准确描述波浪和结构物相互作用的非线性物理现象,在网格划分时自由面附近采用了较为精细的网格。本文采取这样的网格划分意义在于利用了网格的相对运动,无需网格再生、变形,既能保证船体任意幅度的振荡运动,又避免复杂船体周围网格变化,从而提高了计算精度和效率。为了避免水池边壁的波浪反射对振荡辐射运动的船体周围流场的影响,把水池的宽度及长度取得相对较大,并且在离船体远端的区域采用了尺度较大的网格。从数值计算的角度来说这样起到了数值消波的功能,一定程度上可以消去物体振荡运动的辐射波。鉴于船舶水动力系数测试图1三维船模数值波浪水池60中国造船学术论文与分析所需要的监测值的时间序列并不长,笔者未在水池边界区域应用相关的消波技术。图2船体附近网格划分(左为水池顶部,右为对称面)数值波浪水池的边界条件为上面为压力出口条件,船体、水底和侧壁面为壁面边界。船体振荡运动时流场的模拟计算中压力速度采用SIMPLE方法进行迭代求解,动量方程中的瞬态项采用二阶隐格式差分格式,对流项和扩散项的离散都采用二阶迎风差分算法。本文采用VOF方法来追踪自由面波动时,对单元面流量的计算采用的是精确度较高的几何重构法,对两种流体的界面采用分段线性表达。3船舶水动力系数的测试与分析数值波浪水池中船体水动力系数的“测试”(即数值模拟计算)主要借鉴了物理水池中进行船模强迫振荡运动试验的方法,但理论计算不仅比振荡模型物理试验更加容易控制、实现和数据测量,还可以直接采用实船尺度,避免了模型的尺度效应。理论计算中物体作单一模态的强迫振荡运动,且不失一般性,文中只列出船体作垂荡运动时的运动规律,即TTTSINCOSSIN2303303303垂荡加速度垂荡速度垂荡位移5式中Ξ30和分别为垂荡运动的幅值和频率。其它运动,如横荡、纵摇运动等皆可仿此。受迫振荡的船体受力的监测实际上是理论计算中船体物面离散单元受力的积分。若记数值模拟计算中离散的物面单元所受到的压力为P,并假定其中静压力333000SXPGXX,那么物面单元上的动压力就可记为DSPPP。对动压力沿物面S0的积分,船体受到的动力的垂向分量为033DDSFPNS,相应地垂荡运动诱导的船体受到的纵摇力矩就可以写成055DDSFPNS03113DDSPRNRNS,其中N1,N3为物面的法线方向单位分量,指向物体内侧,R1和R3为物面单元的压力中心与转动中心的纵向和垂向距离。船体垂荡运动的线性方程为3333333533533500ABFABF7其中船体所受的力(矩)是振荡辐射运动引起的,无论垂荡力还是纵摇力矩都包含了与振荡运动同相位的力(矩)和与运动相位正交的力(矩),分别表示振荡运动时船体的附加质量力和阻尼力APFPAPFP49卷增刊总第183期郭海强,等数值波浪水池中船舶水动力系数测试与分析技术613033350555SINCOSSINCOSABFFFTFTFFFTFT8式中03F,05F,AF3,BF3,AF5,BF5为三维分解后的幅值。通过对振荡运动的船体周围流场的数值模拟计算,获得力(矩)的时历,再进行相位分解,便可得到船体的附加质量和阻尼系数333230AFA,553230AFA,33330BFB,55330BFB9由此可知,当船体作“纯垂荡”运动时,可以获得垂荡、以及垂荡与纵摇耦合的附加质量和阻尼系数;用同样的方法可以获得“纯纵摇”,“纯纵荡”,“纯横摇”,“纯横荡”及其耦合的附加质量和阻尼系数。4船模水动力学系数的测试与分析实例41测试用船模WIGLEY船型是用数学公式描述的,采用右手坐标系,,,O是纵向坐标,是横向坐标,向右舷为正,是垂向坐标,向下为正。O点在船舯静水面上。无因次化的船型表达式为222428242411111AA10,5050,5050其中WIGLEYIII船模对应的系数2A02,4A0,0,分别对,,各乘以半个船长L/2,半宽B/2和吃水D得到实际的船模尺度,模拟计算用的船模主尺度如表1所示。表1WIGLEYIII模型主尺度船长PPL/M型宽B/M吃水D/M舯剖面系数方形系数300301875066670462242船体水动力系数的计算分析和验证本文就WIGLEYIII船模分别所作的“纯垂荡”、“纯纵摇”及“纯纵荡”运动下的流场进行了数值模拟,以验证本文方法的有效性和准确性。当船体在数值波浪水池中以固定的频率和幅值作单一模态的强迫振荡运动时,通过流场的模拟计算,可以获得相应模态下的力和力矩,以及静力和静力矩。如图3所示,当船模以周期为20S,位移幅值为20MM作“纯垂荡”运动,即以3002SINΠT的位移规律运动时,通过模拟计算获得整个船模的垂向静力和垂荡力的时历曲线。计算中采用了固定的时间步长约为振荡运动周期的六百分之一。图3中的横坐标为时间T,纵坐标为计算所得的垂向静力及垂荡力F。对力和力矩的时历曲线进行相位分解。由于CFD模拟的计算结果都总存在噪声,用最小二乘原理以正弦曲线对平稳的时历曲线拟合,进而获得船体的附加质量和阻尼系数。对应于图3的时历,附加质量与阻尼系数分别为33A=113784,33B=140795。计算发现垂荡和纵摇之间耦合的附加质量和阻尼的量级比垂荡的小67个量级,这与WIGLEY数学船型关于横中剖面对称是一致的。图2WIGLEYIII船型线图

注意事项

本文(数值波浪水池中船舶水动力系数测试与分析技术.doc)为本站会员(dingyx0101)主动上传,人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知人人文库网(发送邮件至[email protected]或直接QQ联系客服),我们立即给予删除!

温馨提示:如果因为网速或其他原因下载失败请重新下载,重复下载不扣分。

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

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

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

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