張 露,肖昌潤(rùn),焦玉超
(海軍工程大學(xué) 艦船工程系,湖北 武漢 430033)
十字舵與X舵潛艇的水動(dòng)力性能數(shù)值比較
張 露,肖昌潤(rùn),焦玉超
(海軍工程大學(xué) 艦船工程系,湖北 武漢 430033)
現(xiàn)代潛艇尾操縱面建筑形式主要選擇十字舵與X舵,這2種形式的舵各有其優(yōu)缺點(diǎn)。在以往關(guān)于這2種舵形潛艇操縱性水動(dòng)力的研究中,有研究者通過(guò)編制潛艇六自由度運(yùn)動(dòng)仿真程序,對(duì)十字舵與X舵在水平面和垂直面的水動(dòng)力性能進(jìn)行了綜合比較,得出X舵的水動(dòng)力性能優(yōu)于十字舵。本文以十字舵和X舵Suboff為研究對(duì)象,通過(guò)CFD數(shù)值方法模擬了潛艇直航運(yùn)動(dòng)和垂直面變攻角運(yùn)動(dòng)這2種具有典型代表意義的運(yùn)動(dòng)情況,分別計(jì)算了2種舵形潛艇的操縱性水動(dòng)力,通過(guò)分析計(jì)算結(jié)果,定量地比較了在舵面積相等的情況下十字舵與X舵潛艇的水動(dòng)力性能,得出與編程仿真相同的結(jié)論:X舵的水動(dòng)力性能優(yōu)于十字舵。
CFD數(shù)值模擬;十字舵潛艇;X舵潛艇;Suboff;水動(dòng)力性能
潛艇尾操縱面發(fā)展至今已40多年,十字舵仍是尾操縱面的基本形式,十字舵與早期的雙槳潛艇共同發(fā)展起來(lái),成為應(yīng)用最早、使用最多的尾操縱面基本形式,如中國(guó)、美國(guó)等國(guó)家的戰(zhàn)斗潛艇基本上采用十字舵形,十字舵有利于對(duì)潛艇航向和深度的直觀操縱,簡(jiǎn)化了潛艇的操縱性分析,其結(jié)構(gòu)和布置相對(duì)簡(jiǎn)單,應(yīng)急情況或戰(zhàn)斗狀態(tài)一般人工駕駛,此時(shí)十字舵的操縱顯得簡(jiǎn)捷、方便、安全。近年來(lái),歐洲一些潛艇強(qiáng)國(guó)如瑞典、德國(guó)、挪威、丹麥等在其新建的潛艇上陸續(xù)采用X舵,并得到了很好的效果。由于X舵的特殊結(jié)構(gòu),相對(duì)于傳統(tǒng)十字舵而言X舵有許多顯著的優(yōu)點(diǎn):X舵在總體布置上可以做到橫向尺度不超寬、垂向尺度不突出基線;X舵可以對(duì)潛艇橫傾進(jìn)行控制;采用X舵后減小了操舵后回轉(zhuǎn)中的艇重、尾重現(xiàn)象;X舵可大大減少尾舵卡住時(shí)造成的嚴(yán)重后果,提高了潛艇的安全性和水下動(dòng)力抗沉能力;X舵的舵裝置容易標(biāo)準(zhǔn)化和系列化等。X舵的缺點(diǎn)是由其操縱特性帶來(lái)的,由于X舵的每個(gè)舵板都有潛浮功能和轉(zhuǎn)首功能,因此X舵操縱不直觀,人工操縱時(shí)不如十字舵方便,X舵所產(chǎn)生的垂向和水平方向的力都對(duì)稱,因此很難調(diào)節(jié)垂直面與水平面對(duì)操縱性的不同要求。因此,選擇十字舵或X舵需要根據(jù)對(duì)潛艇性能的不同需求來(lái)決定。本文通過(guò)CFD數(shù)值計(jì)算方法,定量地比較了舵面積相同的情況下十字舵與X舵的水動(dòng)力性能,為今后潛艇尾操縱面的設(shè)計(jì)提供了理論參考。
潛艇的操縱性水動(dòng)力計(jì)算主要依靠工程估算和母型修正方法,操縱性預(yù)報(bào)主要依靠模型試驗(yàn)方法,由于國(guó)內(nèi)現(xiàn)有的潛艇操縱性水動(dòng)力估算公式形式簡(jiǎn)單,不能充分準(zhǔn)確地反映艇型的變化,而且其適用性和預(yù)報(bào)精度有限,達(dá)不到工程實(shí)用精度,模型試驗(yàn)周期長(zhǎng)耗資大,往往難以滿足進(jìn)度要求,而當(dāng)今計(jì)算流體力學(xué)(CFD)的飛速發(fā)展及其在潛艇操縱性水動(dòng)力預(yù)報(bào)方面的應(yīng)用為解決以上困難開(kāi)辟了有效的新途徑,CFD是以流體力學(xué)和數(shù)值計(jì)算方法的基本理論為基礎(chǔ)而建立的,本文將通過(guò)流體計(jì)算軟件Fluent對(duì)十字舵與X舵Suboff這2種模型的直航運(yùn)動(dòng)和垂直面變攻角(α)運(yùn)動(dòng)進(jìn)行數(shù)值計(jì)算,并比較這2種舵形Suboff的水動(dòng)力性能,同時(shí)對(duì)比文獻(xiàn)中Suboff模型(十字舵)的實(shí)驗(yàn)數(shù)據(jù),驗(yàn)證了CFD數(shù)值計(jì)算方法的可靠性。
1.1 計(jì)算模型
Suboff標(biāo)準(zhǔn)模型艇長(zhǎng)L=4.356 m,最大直徑D=0.508 m,模型包含指揮臺(tái)圍殼、螺旋槳、尾舵3種附體。本文利用Solidworks進(jìn)行幾何建模,得到十字舵Suboff與X舵Suboff的幾何模型,十字舵潛艇的尾操縱面可分為升降舵和方向舵,且升降舵與方向舵呈正交布置,升降舵控制潛艇的深度和縱傾,方向舵用于改變或保持潛艇航向,然而,X舵潛艇不存在升降舵與方向舵,其尾操縱面的4個(gè)舵板呈正交布置,舵軸中心線與潛艇的中線面成±45°角,每個(gè)舵板的偏轉(zhuǎn)都能引起潛艇的潛浮與轉(zhuǎn)向,因此,必須采用多個(gè)(一般是2個(gè)或4個(gè))操舵控制系統(tǒng)綜合控制,以改變潛艇航向或深度。其中,舵面積相等的十字舵與X舵的具體布局方式分別如圖1和圖2所示。
1.2 控制方程及數(shù)值方法
潛艇在粘性流場(chǎng)中運(yùn)動(dòng),從牛頓力學(xué)角度出發(fā),由于水通常被認(rèn)為是不可壓縮粘性流體,故流動(dòng)必遵循質(zhì)量守恒和動(dòng)量守恒定律,用張量的形式表述定常、不可壓流動(dòng)作時(shí)間平均處理的控制方程為:
文中采用RNGk-ε湍流模型封閉雷諾時(shí)均方程,RNGk-ε湍流模型的輸運(yùn)方程為:
近壁面流體采用壁面函數(shù)法,F(xiàn)luent建議Y+集中于(5,25),使用SIMPLE算法對(duì)壓力速度耦合方程組進(jìn)行求解,以有限體積法離散流體運(yùn)動(dòng)的控制方程和湍流運(yùn)動(dòng)方程,其中擴(kuò)散項(xiàng)以中心差分格式差分,對(duì)流項(xiàng)采用2階迎風(fēng)差分格式。
1.3 計(jì)算域設(shè)置、邊界條件設(shè)定及網(wǎng)格劃分
計(jì)算域設(shè)置為一長(zhǎng)方體區(qū)域:距離潛艇前部有1L距離,尾部有2L距離,上下左右均有1L距離,如圖3所示。
邊界條件設(shè)定為:進(jìn)口處為速度入口邊界條件,出口處為壓力出口邊界條件,上下面、左右面均為壁面邊界條件,艇體表面為壁面無(wú)滑移邊界條件。
采用ICEM CFD軟件劃分網(wǎng)格,網(wǎng)格劃分大多采用結(jié)構(gòu)化網(wǎng)格,對(duì)于指揮臺(tái)圍殼和尾舵等這種帶有小角度夾角的端面進(jìn)行網(wǎng)格加密,保證非加密與附近加密區(qū)網(wǎng)格尺寸相當(dāng),有利于計(jì)算插值和物理量的交換,加快收斂速度。十字舵Suboff模型網(wǎng)格總數(shù)為207萬(wàn),X舵Suboff模型網(wǎng)格總數(shù)為212萬(wàn),部分網(wǎng)格細(xì)節(jié)如圖4和圖5所示。
1.4 計(jì)算結(jié)果及分析
1.4.1 兩種舵形Suboff直航阻力數(shù)值模擬
試驗(yàn)的計(jì)算結(jié)果以不同速度下阻力值的形式表示出來(lái),如表1所示。
結(jié)果分析:
表 1 直航運(yùn)動(dòng)艇體阻力計(jì)算結(jié)果Tab. 1 The calculation results of the straight motion
1)從十字舵計(jì)算值與試驗(yàn)值的對(duì)比可知,相對(duì)誤差均在4%內(nèi),說(shuō)明此CFD數(shù)值方法可靠,選擇的RNG k-ε湍流模型、邊界條件設(shè)定及網(wǎng)格劃分均合適。
1.4.2 計(jì)算過(guò)程中的表面壓力分布
結(jié)果分析:
通過(guò)分析v=3.045 m/s,2種舵型Suboff的表面壓力分布,易發(fā)現(xiàn)X舵Suboff的表面壓力要小于十字舵Suboff。
2.1 計(jì)算模型
計(jì)算模型仍取圖1和圖2中的十字舵與X舵Suboff幾何模型。
2.2 控制方程及數(shù)值方法
控制方程仍選擇用張量的形式表述的定常、不可壓流動(dòng)作時(shí)間平均處理的控制方程:
當(dāng)攻角|α|<8°時(shí),采用RNGk-ε湍流模型封閉雷諾時(shí)均方程,當(dāng)攻角|α|≥8°時(shí),采用SSTk-ω湍流模型封閉雷諾時(shí)均方程,SSTk-ω湍流模型的輸運(yùn)方程為:
同樣以有限體積法離散流體運(yùn)動(dòng)的控制方程和湍流運(yùn)動(dòng)方程,其中擴(kuò)散項(xiàng)以中心差分格式差分,對(duì)流項(xiàng)采用2階迎風(fēng)差分格式,使用SIMPLE算法對(duì)壓力速度耦合方程組進(jìn)行求解。
2.3 計(jì)算域設(shè)置、邊界條件設(shè)定及網(wǎng)格劃分
計(jì)算域設(shè)置為一長(zhǎng)方體區(qū)域:距離潛艇前部有1L距離,尾部有2L距離,上下有1L距離,左右有1L距離。
邊界條件設(shè)定為:進(jìn)口處為速度入口邊界條件,出口處為壓力出口邊界條件,上下面、左右面均為壁面邊界條件,艇體表面為壁面無(wú)滑移邊界條件。
同樣地,采用ICEM CFD軟件劃分網(wǎng)格,網(wǎng)格劃分大多采用結(jié)構(gòu)化網(wǎng)格,對(duì)于指揮臺(tái)圍殼和尾舵等這種帶有小角度夾角的端面進(jìn)行網(wǎng)格加密,保證非加密與附近加密區(qū)網(wǎng)格尺寸相當(dāng),有利于計(jì)算插值和物理量的交換,加快收斂速度。由于本文主要目的是比較2種舵形的水動(dòng)力性能,所以對(duì)2種舵形Suboff均在RNG k-ε與SST k-ω這2種湍流模型中取相同網(wǎng)格,其中十字舵Suboff的網(wǎng)格總數(shù)為224萬(wàn),X舵Suboff的網(wǎng)格總數(shù)為232萬(wàn),部分網(wǎng)格細(xì)節(jié)如圖8~圖10所示。
2.4 計(jì)算結(jié)果及分析
2.4.1 兩種舵形Suboff垂直面變攻角數(shù)值模擬
試驗(yàn)的計(jì)算結(jié)果以不同攻角下升力值的形式表示出來(lái),如表2所示。
結(jié)果分析:
表 2 垂直面變攻角試驗(yàn)艇體升力計(jì)算結(jié)果Tab. 2 Test results of vertical drift angle experiment
1)從表中易看出,對(duì)于2種舵形潛艇,升力值均隨著攻角的增大而增大。
2)當(dāng)攻角|α|<8°時(shí),十字舵計(jì)算結(jié)果的相對(duì)誤差在8%之內(nèi),滿足工程需要;當(dāng)攻角|α|≥8°時(shí),十字舵計(jì)算結(jié)果的相對(duì)誤差超過(guò)了10%,不滿足工程需要,原因是當(dāng)攻角|α|≥8°時(shí),選擇SSTk-ω湍流模型采用了與RNGk-ε湍流模型相同的網(wǎng)格,實(shí)際上,SSTk-ω湍流模型的網(wǎng)格總數(shù)要大于RNGk-ε湍流模型的網(wǎng)格總數(shù),但是本文主要目的是比較2種舵形的水動(dòng)力性能,所以可以認(rèn)為此CFD數(shù)值方法滿足精度要求。
2.4.2 計(jì)算過(guò)程中的流場(chǎng)分布
結(jié)果分析:從2種舵型Suboff的流場(chǎng)分布可以發(fā)現(xiàn),X舵Suboff的尾部流場(chǎng)要比十字舵尾部流場(chǎng)均勻,因此可以推斷X舵比十字舵更有利于螺旋槳的推進(jìn)效率。
1)通過(guò)將十字舵SUBOFF模型的CFD數(shù)值計(jì)算結(jié)果與試驗(yàn)值進(jìn)行比較,可以發(fā)現(xiàn);潛艇直航運(yùn)動(dòng)時(shí),CFD計(jì)算精度均在4%絕對(duì)誤差范圍之內(nèi);垂直面變攻角運(yùn)動(dòng)時(shí),雖然當(dāng)攻角|α|≥12°時(shí),由于湍流模型計(jì)算不能準(zhǔn)確捕捉到大攻角時(shí)的流動(dòng)分離,而且大攻角時(shí)指揮臺(tái)圍殼引起背流區(qū)分離,影響舵效及伴流場(chǎng),從而導(dǎo)致計(jì)算值偏離實(shí)驗(yàn)值超過(guò)10%,但是當(dāng)攻角|α|<12°時(shí),CFD計(jì)算精度均在8%絕對(duì)誤差范圍之內(nèi),而且本文中CFD數(shù)值模擬垂直面變攻角運(yùn)動(dòng)的攻角選擇范圍只需滿足|α|<10°,所以可以認(rèn)為本文采用的CFD數(shù)值方法可靠。
2)在潛艇直航運(yùn)動(dòng)結(jié)果中,通過(guò)對(duì)比十字舵Suboff與X舵Suboff阻力系數(shù)計(jì)算值,發(fā)現(xiàn)X舵S u boff的阻力性能優(yōu)于十字舵Suboff,由于2個(gè)Suboff模型只存在尾操縱面建筑形式的差別,且舵面積相等,所以可以認(rèn)為X舵潛艇比十字舵潛艇的阻力性能好。
3)在垂直面變攻角運(yùn)動(dòng)結(jié)果中,通過(guò)對(duì)比十字舵SUBOFF與X舵SUBOFF升力系數(shù)計(jì)算值,發(fā)現(xiàn)X舵SUBOFF的升力性能優(yōu)于十字舵SUBOFF,同理,可以認(rèn)為X舵潛艇比十字舵潛艇的升力性能好;另外,通過(guò)對(duì)比流場(chǎng)的計(jì)算結(jié)果,發(fā)現(xiàn)X舵潛艇的尾部流場(chǎng)要好于十字舵潛艇,更有利于螺旋槳的推進(jìn)效率。
4)本文運(yùn)用CFD數(shù)值模擬方法比較了十字舵與X舵潛艇的部分水動(dòng)力性能,得出了X舵潛艇的阻力性能、升力性能均優(yōu)于十字舵潛艇,但是在以后的工作中還可以進(jìn)一步比較其他水動(dòng)力性能,比如比較十字舵與X舵這2種舵型對(duì)潛艇操縱性的不同影響。
5)潛艇尾舵對(duì)艇的操縱性和安全性起著關(guān)鍵作用,在研制發(fā)展新型潛艇中占有重要地位,隨著新型潛艇線性設(shè)計(jì)的不斷發(fā)展,現(xiàn)多為單槳尖回轉(zhuǎn)體尾型,這種尾型給潛艇尾舵的設(shè)計(jì)提供了更加充分的想象空間,所以今后可以充分運(yùn)用CFD數(shù)值方法,來(lái)研究新舵型(如“米”字舵、“木”字舵、“干”字舵等)的水動(dòng)力性能,推動(dòng)潛艇新舵型的研發(fā)。
[ 1 ]施生達(dá), 王京齊. 潛艇操縱面的建筑形式評(píng)述[J]. 海軍工程學(xué)院學(xué)報(bào), 1995(1): 78-86.
[ 2 ]JEANS T, WATT G, GERBER A, et al. High-resolution Reynolds-averaged navier-stokes flow predictions over axisymmetric bodies with tapered tails[J]. AIAA Journal, 2009, 47(1): 19–32.
[ 3 ]PHILLIPS A, TURNOCK S, FURLONG, et al. Influence of turbulence closure models on the vertical flow field around a submarine body undergoing steady drift [J]. Journal of Marine Science and Technology, 2010, 15(3): 201–217.
[ 4 ]SEN D. A study on sensitivity of maneuverability performance on the hydrodynamic coefficients for submerged bodies[J]. Journal of Ship Research, 2000, 44(3): 186–196.
[ 5 ]馮德生. 潛艇“X”舵發(fā)展概況和瑞典“X”舵的研究應(yīng)用簡(jiǎn)述[Z]. 出國(guó)考察論文集. 九江: 707所分部, 1991.
[ 6 ]ITTC 1996 report of the maneuverability committee [C]// Proceedings, 20th International Towing Tank Conference. Trondheim Norway, 1996. 347–398.
[ 7 ]KIM S, RHEE B, MILLER R. Anatomy of turbulent flow around DARPA SUBOFF body in a turning maneuver using high-fidelity RANS computations[J]. International shipbuilding Progress, 2013, 60(1): 207–231.
[ 8 ]SUNG C, JIANG M, RHEE B, et al. Validation of the flow around a turning submarine[C]// Fukouka: The Twenty-Fourth Sympo. on Naval Hydrodynamics, 2002.
[ 9 ]ZHANG J, MAXWELL J, GERBER A, et al. Simulation of the flow over axisymmetric submarine hulls in steady turning[J]. Ocean Engineering, 2013, 57(1): 180–196.
[10]PAN Yu-chun, ZHANG Huai-xin. Numerical prediction of submarine hydrodynamic coefficients using CFD simulation [J]. Journal of Hydrodynamic, 2012, 24 (6): 840–847.
[11]恩蓋伊恩TD S, 高英D, 博欽斯基. 調(diào)整片輔助操縱面在船舶航行中的應(yīng)用[C]// 第六屆國(guó)際潛艇譯文集, 中國(guó)船舶科學(xué)研究中心, 2001. 7.
Numerical comparison on hydrodynamic performance of cross rudder and X rudder submarine
ZHANG Lu, XIAO Chang-run, JIAO Yu-chao
(Naval Engineering Department, Naval University of Engineering, Wuhan 430033, China)
Now submarine main cross rudder and tail control surfaces architectural form X rudder, the rudder have their advantages and disadvantages of these two forms. In the past about the two kinds of steering maneuverability hydrodynamic studies, researchers by submarine six degrees of freedom movement simulation program, the cross rudder and X rudder in the horizontal plane and vertical plane of hydrodynamic performance makes a comprehensive comparison, concludes that X rudder hydrodynamic is better than cross rudder. Based on cross rudder and X rudder Suboff research object, through CFD numerical method to simulate the submarine direct movement and vertical angle of attack motion trajectory of these two kinds of typical significance, and separately calculated the two steering maneuverability hydrodynamic shape submarines, through analysis of calculation results, quantitatively compared under the condition of the rudder area equal cross rudder and the hydrodynamic performance of X rudder submarine, and programing simulation the same conclusion: the hydrodynamic performance of X rudder is better than that of cross rudder.
CFD numerical simulation;cross rudder submarine;X rudder submarine;Suboff;hydrodynamic performance
U661.1
A
1672 – 7649(2017)07 – 0024 – 05
10.3404/j.issn.1672 – 7649.2017.07.005
2016 – 09 – 02;
2016 – 09 – 26
張露(1990 – ),男,碩士研究生,研究方向?yàn)闈撏Р倏v面。