基于重叠网格方法的中型邮轮减摇鳍数值和试验分析.pdf

收稿日期:修回日期:录用日期:作者简介:姚汝林(),博士生,主要从事船舶与海洋工程先进制造研究通信作者:余龙,副教授,电话(T e l):;E m a i l:y u l o n e s j t u e d u c n 文章编号:()S D O I:/j c n k i j s j t u S 基于重叠网格方法的中型邮轮减摇鳍数值和试验分析姚汝林,樊奇东,余龙,汪学锋,(上海交通大学 海洋工程国家重点实验室,上海 ;上海交通大学 船舶海洋与建筑工程学院,上海 )摘要:减摇鳍是解决船舶横摇的有效办法之一,可提升邮轮的舒适性该装置一般由一对或两对可伸缩的鳍组成,对称地安装在船舶的两侧舭部,通过控制算法依靠鳍的水动力升力提供运动阻尼然而,不同航速和不同攻角下减摇鳍与船舶的耦合分析还较少,尤其是船体流场对鳍阻力的影响尚缺乏研究利用模型试验和数值分析方法研究某邮轮的阻力性能和鳍的运动影响,比较数值和模型试验结果,对减摇鳍引起的阻力和船鳍耦合运动进行研究,为减摇鳍的设计和选择、邮轮设计和性能预报提供了参考关键词:减摇鳍;模型试验;邮轮;不同攻角;不同航速中图分类号:U 文献标志码:AN u m e r i c a l a n dT e s t i n gA n a l y s i so f F i nS t a b i l i z e r so fAM e d i u mS i z e dC r u i s eS h i pw i t hO v e r s e tG r i d sY A OR u l i n,F ANQ i d o n g,Y UL o n g,WANGX u e f e n g,(S t a t eK e yL a b o r a t o r yo fO c e a nE n g i n e e r i n g,S h a n g h a i J i a oT o n gU n i v e r s i t y,S h a n g h a i ,C h i n a;S c h o o l o fN a v a lA r c h i t e c t u r e,O c e a na n dC i v i lE n g i n e e r i n g,S h a n g h a i J i a oT o n gU n i v e r s i t y,S h a n g h a i ,C h i n a)A b s t r a c t:F i ns t a b i l i z e r s r e p r e s e n t a ne f f e c t i v es o l u t i o nt oa d d r e s s t h e r o l lm o t i o no f s h i p sa n d i m p r o v e t h ec o m f o r to fp a s s e n g e r so nc r u i s es h i p s T h e s ed e v i c e s t y p i c a l l yc o m p r i s eo fo n eo r t w op a i r so f r e t r a c t a b l ef i n s,s y mm e t r i c a l l ym o u n t e do ne i t h e rs i d eo f t h es h i p,w h i c hu t i l i z eh y d r o d y n a m i c l i f t t od a m p e nm o t i o nt h r o u g hac o n t r o l a l g o r i t h mH o w e v e r,c o u p l i n ga n a l y s i so f f i ns t a b i l i z e r sa n ds h i p sa tv a r i o u ss p e e d sa n da n g l e so fa t t a c kr e m a i n sl i m i t e d,p a r t i c u l a r l y w i t hr e g a r dt ot h ei m p a c to ft h eh u l lf l o wf i e l do nf i nr e s i s t a n c e T h i sp a p e r i n v e s t i g a t e st h ed r a gp e r f o r m a n c ea n dt o w i n gm o t i o no fac r u i s es h i pu s i n gm o d e lt e s t sa n dn u m e r i c a l a n a l y s i sm e t h o d s,a n dc o m p a r e st h er e s u l t so f t h en u m e r i c a l a n dm o d e l t e s t s I ta l s oe x a m i n e st h ed r a g r e s u l t i n g f r o mf i ns t a b i l i z e r s a n d t h e c o u p l i n gm o t i o no f t h e s h i p,o f f e r i n g i n s i g h t f o r t h ed e s i g na n ds e l e c t i o no f f i ns t a b i l i z e r s,c r u i s es h i pd e s i g n,a n dp e r f o r m a n c ep r e d i c t i o n K e yw o r d s:f i ns t a b i l i z e r;m o d e l t e s t;r e s i s t a n c e;c r u i s es h i p;d i f f e r e n t f i na n g l e s;d i f f e r e n t s h i ps p e e d第 卷 增刊 年 月上 海 交 通 大 学 学 报J OUR NA LO FS HAN GHA I J I AOT ON GUN I V E R S I T YV o l S u p O c t 增刊姚汝林,等:基于重叠网格方法的中型邮轮减摇鳍数值和试验分析 减摇鳍是解决船舶横摇的有效办法之一,可提升乘客的舒适感该装置一般由一对或两对可伸缩的鳍组成,对称地安装在船舶的两侧舭部,利用控制算法依靠鳍的水动力升力提供运动阻尼,削减船舶横摇许多学者对减摇装置进行了研究,比如减摇水舱、移动重量、减摇陀螺以及舭龙骨,可以在较低航速时平衡船舶然而,应用这些技术的难点在于船舶的尺寸、功率和成本典型的减摇鳍结构如图所示一套减摇鳍通常有一对分布在两舷的鳍片,由液压单元和本地控制单元(L C U)驱动,由一个驾驶台控制单元(B C P)、一个主控制单元(MC U)和一个横摇传感器单元(RM S U)构成关于鳍的选择与其他减摇装置的比较参见文献图可收放减摇鳍装置F i g R e t r a c t a b l e f i nd e v i c eG a i l l a r d e和R a m等针对多种不同的带减摇鳍船型进行了模型和数值模拟试验,研究减摇鳍升阻力性能但是在不同角度、不同速度下鳍片的升力和阻力特性研究还较为欠缺 P e r e z等的研究表明,有效攻角应被限制以避免出现动态失速有学者对一艘鳍翼型为NA C A 的水翼巡逻艇进行数值分析,发现在、和 攻角,船速为 k n(k n m/h)时,鳍的攻角增加引起阻力增 大攻 角 时,船 体 阻 力 增 加 总 量 可 达 K i m等对安装减摇鳍的邮轮在遭遇波浪条件下的运动进行了数值分析当船速为 k n(中等速 度)时,减 摇 鳍 控 制 系 统 将 横 摇 振 幅 降 低 了 当船速为k n(低速)时,横摇振幅降低了 为进一步研究船鳍耦合,以一艘中型邮轮为对象,使用重叠网格计算模拟方法对不同攻角下带减摇鳍的邮轮总阻力进行分析,并通过模型试验对计算流体动力学(C F D)结果进行验证同时也研究了在静水中船鳍耦合作用以及带减摇鳍船舶的升沉和横摇运动,为计算带减摇鳍船舶的设计提供参考 带减摇鳍邮轮船型介绍以一艘中型邮轮为研究对象,该船 的模型主要参数如表所示表船模型主要参数T a b M a i np a r a m e t e r so f t h e s h i pm o d e l参数数值垂线间长,LP P/m 型宽,B/m 吃水,T/m 排水量,/t 湿表面积,Sw/m 重心,LC G/m 重心,KG/m 鳍切片翼型NA C A 鳍弦长,C/m 长宽比 鳍长,LF/m 三维模型如图所示,船中带有一对减摇鳍,鳍安装在水线之下,通过沿自身轴的旋转实现不同角度图带有一对减摇鳍的船舶模型F i g S h i pm o d e lw i t hap a i ro f f i ns t a b i l i z e r s 上海交通大学学报第 卷在使用过程中,船两侧的鳍片并不总是对称地旋转当横摇达到最大程度时,两侧的鳍会发生不对称旋转,以便获得更大抗摇力矩;比如,左舷的鳍向上旋转(顺时针旋转),而右舷的鳍则向下旋转(逆时针旋转),可旋转幅度为 ,如图所示图鳍片的不同的旋转方向F i g F i n sw i t hd i f f e r e n t r o t a t i o nd i r e c t i o n s数值方法由于带减摇鳍邮轮航行时强非线性与黏性的影响,为了获得高精度的模拟,应用雷诺平均纳维斯托克斯(R AN S)求解器和重叠网格进行数值分析 控制方程采用通用C F D软件S T A R C CM开展带减摇鳍邮轮的黏性流动模拟,使用有限体积法方法捕捉自由液面 R AN S方法被证实在船舶的阻力计算中拥有较高精准度 连续性和动量守恒方程如下:uixi()(ui)txj(ujui)pxiFixj(j iRj i)()式中:ui、uj分别为雷诺平均速度分量;xi、xj分别为位置分量;为流体密度;t为时间;p为压力;Fi为混合函数;j i为应力分量;Rj iu ju i为雷诺应力项或湍流应力项;在雷诺平均方程中,这个未知量使方程不封闭为了使方程封闭,加入S S TK Om e g a湍流模型 重叠网格技术重叠网格又称O v e r s e t,是将模型中流动区域划分成几何边界比较简单的子区域形成网格,然后嵌套进背景网格中计算通过在重叠网格区域相互的插值,使得子区域的网格在重叠边界进行数据交换,从而完成整个流场的求解使用重叠网格技术能在特定位置进行局部加密,使物体不受与网格之间的拓扑结构上的约束,能够让物体在计算域内自由运动减摇鳍作为附体被船舶安装后,会对船体的受力和运动产生一定影响利用重叠网格处理这种耦合流场问题下的船舶阻力和运动,在过往研究中已经被证明能够获得较好结果 计算域及网格划分计算针对带鳍船,建立两套网格第一套为背景网格,尺寸为 LP P LP PLP P,网格数量为 ,背景网格的入口、底部和顶部的边界条件设置为速度入口,出口设置为压力出口第二套为船体带鳍附近网格,尺寸为 LP P B D,D为船的型深,网格量为 ,计算域如图所示计算不同攻角减摇鳍使用相同的背景网格,不同攻角的船体带鳍网格根据几何作变化,如图所示为顺利捕捉减摇鳍周围流场,重叠网格下相同船速下鳍片表面的局部流线、Y的比较分别如图和图所示Y为研究流体在近壁面区域的流动状态而引入的无量纲距离图计算域设置F i g C a l c u l a t i o nd o m a i n增刊姚汝林,等:基于重叠网格方法的中型邮轮减摇鳍数值和试验分析 图船体带鳍周围网格F i g M e s ha r o u n dt h eh u l lw i t hf i n s图带鳍的船舶模型F i g S h i pm o d e lw i t hf i n s图 减摇鳍流线信息F i g S t r e a m l i n e i n f o r m a t i o no f f i ns t a b i l i z e r图船和 鳍的YF i g Y o fv e s s e l a n d f i n s相关计算表明,在Y小于 时,能捕捉到减摇鳍附近的流线细节,可获得阻力和运动较好的模拟结果 试验设置在上海交通大学多功能船模拖曳水池进行阻力系列试验拖曳水池的尺寸 m m m,拖车最大速度为 m/s 图显示了整个船模的几何和减摇鳍的安装位置:在船中部安装了一对减摇鳍,纵向位置距离船尾 上海交通大学学报第 卷边界约 Lp p,在左舷和右舷各有一只,垂向位置比基线高出 D,倾斜度为 鳍片的水翼采用NA C A系列,工作角度范围是 鳍片的几何参数见表表中列出各种试验状态,主要研究船鳍耦合作用下的阻力试验工况涵盖了不同鳍片攻角和不同船速下有或无鳍片情况表模型试验工况T a b C o n d i t i o n so fm o d e l t e s t工况编号航速/k n弗劳德数(F r)模型速度/(ms)鳍攻角/()数值计算与试验结果比较分析 C F D结果验证实验流体力学(E F D)和数值结果如图所示,该图比较了没有鳍片和有鳍片攻角为 和 的实验结果 C F D的结果与试验数据也较为吻合,最大误差范围在 以内图带鳍和不带鳍船的阻力曲线比较F i g D r a gc u r v e so f s h i pw i t ho rw i t h o u t f i n s表为在 k n即服务航速时,利用重叠网格模拟 计 算 阻 力 与 动 态 流 体 相 互 作 用(D F B I)模拟、模型试验结果的比较 D F B I模拟与模型试验在 与 时的误差分别为 和 ;重叠网格计算结果与模型试验在 与 时的误差分别为 和 计算结果与试验结果吻合度较高,误差均控制在以内不同攻角的鳍平均阻力增量为 ,而在最大鳍片角度下阻力增量高达 船鳍耦合的C F D计算误差低于,而仅计算均匀流场中鳍的阻力,其结果会高估约 表阻力试验结果对比T a b C o m p a r i s o no f r e s u l t so f r e s i s t a n c e t e s t减摇鳍攻角/()阻力/N模型试验 模型试验 D F B I模拟 D F B I模拟 重叠网格模拟 重叠网格模拟 鳍的阻力评估在指定的雷诺数(R e)下,Uf为自由流速度,A为鳍片面积,当FL和FD分别代表在一定攻角下鳍片所引起的提升力和阻力,则升力系数CL和阻力系数CD可以定义如下:CLFL A Uf()CDFD A Uf()由于局部速度的不稳定,引入了平均R e的评估:R eUfL()式中:为动态黏性系数;L为鳍片长度根据鳍片的局部雷诺数,计算鳍片的阻力,并与实验结果进行比较鳍片的阻力如表所示表模型状态下减摇鳍 时阻力T a b D r a go f f i n i nt h em o d e l s t a t eF rR e减摇鳍阻力/N试验获得值C F D计算值不考虑船体计算值 为了详细分析鳍片效应,采用重叠网格计算不同船速下鳍片表面的局部压力和摩擦力并进行比较,如图 所示分析鳍片的摩擦因数可知,随着船航行速度增加,在 的鳍片上有两个集中的高摩擦区域,其中相对更高的区域靠近船体然而,在 鳍片上,摩擦因数高的区域集中在鳍的导边从压力分布分析,随着鳍片攻角较大,高压区集中在导边,以提供升力增刊姚汝林,等:基于重叠网格方法的中型邮轮减摇鳍数值和试验分析 图 在不同的速度下船鳍的压力和表面摩擦因数的详细分析F i g D e t a i l e da n a l y s i so f f i np r e s s u r ea n ds k i nf r i c t i o nf a c t o ra td i f f e r e n t s p e e d s 静水中带鳍船运动分析在静水中 k n服务航速下,对不同攻角带鳍船体的运动进行模拟在速度较低时,船舶的升沉和纵倾不会影响船舶的阻力结果减摇鳍在船中布置时主要影响船舶的横摇,故在模拟时开放了船舶的升沉和横摇两个自由度,计算结果参见表由表可知,减摇鳍为 时,C F D计算与模型试验在船舶运动误差较小,平均横摇误差为 ,说明使用重叠网格能够较好计算出带减摇鳍船舶在航行时的运动状态表 k n航速时船舶运动的试验与数值结果对比T a b C o m p a r i s o no ft e s t i n ga n dn u m e r i c a lr e s u l t so fs h i pm o t i o n sa t k n类型升沉/c m横摇/()减摇鳍攻角/()模型试验 重叠网格 结论减摇鳍是一种有效的提升船舶舒适性的装置,在邮轮上应用广泛减摇鳍的效果以及船鳍耦合作用的影响一直是工程界研究的热点问题以一艘中型邮轮及其减摇鳍为研究对象,通过模型试验和使用重叠网格进行数值计算,分析对比结果,得出以下结论:()试验提供了可靠的船鳍作用下阻力结果在不同船速下,不同攻角的鳍平均阻力增量超过,而在最大鳍片角度下阻力增量高达 由试验得出,鳍片阻力在 攻角比 下增加倍试验结果表明鳍片诱导阻力影响显著,在设计中应该予以重视()采用重叠网格C F D进行非定常分析计算了船鳍耦合和非耦合情况的鳍片阻力船鳍耦合的C F D计算误差低于,而仅计算均匀流场中鳍的阻力,其结果会高估约,船鳍耦合效应不可忽略根据详细的流场信息和计算结果,发现随着速度提高,减摇鳍的高压区向前缘移动,而摩擦区向后缘移动()对静水中减摇鳍的数值计算方法进行了研究,发现计算时减摇鳍部分Y需小于 ,才能精确捕捉船鳍耦合运动的流场,获得精度较高的船舶运动数据使用重叠网格对不同攻角减摇鳍进行分析,有利于理解船鳍耦合作用下鳍的阻力并提出高保真数值分析使用重叠网格的数值计算与试验结果吻合 上海交通大学学报第 卷较好,对实验中船体运动预报误差较小,为后续进一步计算波浪中 的带减摇 鳍船体 运 动 分 析 提 供 了基础参考文献:L I ANGLH,Z HAOP,Z HAN GST,e t a l S i m u l a t i o na n da n a l y s i so fM a g n u sr o t a t i n gr o l ls t a b i l i z e ra t l o ws p e e dJ O c e a nE n g i n e e r i n g,:G A I L L A R D EGD y n a m i cb e h a v i o u ra n do p e r a t i o n a ll i m i t so f s t a b i l i z e r f i n sC t hC o n g r e s so f I n t e r n a t i o n a l M a r i t i m e A s s o c i a t i o n o ft h e M e d i t e r r a n e a n C r e t a,G r e e c e R AMBRR,S UR E N D R ANS,L E ESK C o m p u t e ra n de x p e r i m e n t a l s i m u l a t i o n so nt h e f i ne f f e c to ns h i pr e s i s t a n c eJ S h i p s&O f f s h o r eS t r u c t u r e s,():P E R E ZT,B L ANK E MS h i pr o l ld a m p i n gc o n t r o lJ A n n u a lR e v i e w si nC o n t r o l,():KU S D I ANA W,S EM I N,MA D EA R I ANAI,e t a lL o c a t i o na n a l y s i so fp a t r o l b o a t f i ns t a b i l i z e rb a s e do nn u m e r i c a lm e t h o dJ I O PC o n f e r e n c eS e r i e s:E a r t h&E n v i r o n m e n t a l S c i e n c e,():K I MJH,K I M Y HM o t i o nc o n t r o l o f ac r u i s e s h i pb yu s i n ga c t i v es t a b i l i z i n gf i n sJ P r o c e e d i n g so f t h eI n s t i t u t i o no fM e c h a n i c a lE n g i n e e r s,P a r tM:J o u r n a lo fE n g i n e e r i n gf o rt h eM a r i t i m eE n v i r o n m e n t,():吴乘胜,陈雄,孙立宪,等静水中自由船模拖曳C F D模拟 方 法研 究J船 舶力 学,():WUC h e n g s h e n g,CHE NX i o n g,S UNL i x i a n,e t a lN u m e r i c a l s i m u l a t i o no f f r e e s h i pm o d e l t o w e d i ns t i l lw a t e rJ J o u r n a lo fS h i p M e c h a n i c s,():L E ES,RHE EKP,CHO IJWD e s i g no ft h er o l ls t a b i l i z a t i o nc o n t r o l l e r,u s i n gf i ns t a b i l i z e r sa n dp o dp r o p e l l e r sJ A p p l i e dO c e a nR e s e a r c h,():曹阳,朱仁传,蒋银,等波浪中KV L C C 运动与阻力增加的C F D计算及分析J哈尔滨工程大学学报,():C A O Y a n g,Z HU R e n c h u a n,J I ANG Y i n,e ta lC F Dv e r i f i c a t i o na n da n a l y s i so fa d d e dr e s i s t a n c ea n dm o t i o n so fKV L C C i nh e a d w a v e sJJ o u r n a lo fH a r b i nE n g i n e e r i n gU n i v e r s i t y,():V E R S T E E G H K,MA L A L A S E K E R A WA ni n t r o d u c t i o nt oc o m p u t a t i o n a lf l u i dd y n a m i c s:T h ef i n i t ev o l u m em e t h o dM n de d H a r l o w,E n g l a n d;N e wY o r k,U S A:P e r s o nE d u c a t i o nL t d,ME N T E RFRT w o e q u a t i o ne d d y v i s c o s i t yt u r b u l e n c em o d e l sf o re n g i n e e r i n ga p p l i c a t i o n sJA I A AJ o u r n a l,():王慧,朱仁传,杨云涛基于R AN S方程的高速滑行艇阻力计算C第三十届全国水动力学研讨会暨第十五届全国水动力学学术会议,北京:海洋出版社,:WAN G H u i,Z HU R e n c h u a n,YAN G Y u n t a o T h er e s i s t a n c ec a l c u l a t i o no f t h eh i g h s p e e dp l a n i n gb a s e do nR AN Se q u a t i o n sCT h e t hN a t i o n a lS y m p o s i u mo nH y d r o d y n a m i c s&t hN a t i o n a lC o n f e r e n c eo nH y d r o d y n a m i c s,B e i j i n g,C h i n a:C h i n aO c e a nP r e s s,王建华基于重叠网格技术的船舶操纵运动直接数值模拟D上海:上海交通大学,WAN GJ i a n h u a D i r e c t s i m u l a t i o n so fs h i pm a n e u v e ru s i n go v e r s e tg r i dt e c h n i q u eD S h a n g h a i:S h a n g h a iJ i a oT o n gU n i v e r s i t y,YA ORL,YU L,F AN Q D,e ta lE x p e r i m e n t a la n dn u m e r i c a l r e s i s t a n c ea n a l y s i s f o r ac r u i s e s h i pW/Of i ns t a b i l i z e r sJ J o u r n a l o fM a r i n eS c i e n c e&E n g i n e e r i n g,():