陳 忱,湯霄揚,張海瑛,高曉磊,鄭少雄,林昭友
(1.中國船舶及海洋工程研究院,上海 200011;2.廣州海洋地質(zhì)調(diào)查局,廣東 廣州 510075)
海洋地質(zhì)調(diào)查船在開展海洋地質(zhì)調(diào)查等作業(yè)任務(wù)時需要船體保持穩(wěn)定的運動姿態(tài)。相比較于傳統(tǒng)側(cè)推,伸縮型側(cè)推(RDT)裝置,可以提高船舶首部操縱性,滿足船舶動力定位功能,且在尾部螺旋槳有故障時,可作為推進器維持船舶較低的航速前進。近年來,在海洋地質(zhì)作業(yè)類的船舶上應(yīng)用較多。
目前國內(nèi)的相關(guān)研究集中在側(cè)推裝置的振動評估[1]、安裝工藝[2]、減振降噪[3-4]和側(cè)推載荷估算等方面,在RDT 裝置底部封板設(shè)計載荷和結(jié)構(gòu)強度分析方面的研究較少。胡燕萍[5]用經(jīng)驗公式方法估算了海洋科學(xué)考察船側(cè)推筒體封蓋載荷;魏躍峰[6]通過CFD 數(shù)值仿真方法分析了球鼻艏船型船首側(cè)推孔推力減縮的產(chǎn)生原因及其影響因素;李齊垚[7]通過CFD 數(shù)值仿真和模型試驗開展了單側(cè)推裝置在不同流速、流向下的推力減額的研究;韋紅剛[8]通過CFD 數(shù)值仿真方法計算了不同尺度下的側(cè)推孔、軸包套、軸支架的阻力增量。
RDT 裝置底部封板同時承受波浪載荷和靜水載荷的作用,其結(jié)構(gòu)強度將直接影響該類船型的聲學(xué)設(shè)備作業(yè)可靠性和安全性。因此,有必要研究RDT 裝置底部封板的載荷成分和載荷特征,提出RDT 裝置設(shè)計載荷確定方法,并在此基礎(chǔ)上開展底部封板結(jié)構(gòu)強度評估,保證結(jié)構(gòu)安全。
本文的研究對象RDT 裝置安裝于某海洋地質(zhì)調(diào)查船,該船長85 m,寬22 m,型深8.0 m,吃水6.3 m,船體線型設(shè)計為首部采用直行。
船首110#肋位處安裝有動力定位用的RDT 裝置,RDT 裝置底封板作業(yè)承受水流為5 kn,最大作業(yè)工況5 級海況。RDT 裝置采用輪緣推進形式,輪緣內(nèi)部具備槳轂,輪緣導(dǎo)流罩底部有圓形封板即底封板,導(dǎo)流罩和底封板之間有翼板材設(shè)置的加強筋。RDT 裝置結(jié)構(gòu)及裝置與船體安裝配合情況如圖1 所示。
RDT 裝置底封板承受的載荷主要為波浪載荷、靜水載荷和輪緣導(dǎo)流罩的傳遞載荷,波浪載荷包括水流沖擊載荷和船舶在波浪中搖蕩受到的動態(tài)載荷,與RDT 裝置在船舶布置情況、航速、流場和工作海況相關(guān),靜水載荷與裝置的吃水深度有關(guān),輪緣導(dǎo)流罩的傳遞載荷與結(jié)構(gòu)連接方案相關(guān)。
RDT 裝置屬于航行船舶附體,對于RDT 裝置受載分析,需要在整個船體范圍整體考核。對排水型船舶自由面擾流問題的計算分析是將自由面流動作為兩相流來分析,自由面為水和空氣的交界面,使用VOF 方法計算。該數(shù)學(xué)模型的控制方程包括:連續(xù)性方程、動量方程、體積分數(shù)方程,以及湍流模型方程。
連續(xù)性方程和動量方程分別為:
式中,t為時間;ρ為密度;V為控制體;S為控制體的面積;Ud為控制體面積法向向量的速度;U和p分別為速度和壓力;Ui為xi方向上的平均速度分量;τij為黏性應(yīng)力張量,gi為重力矢量;Ii和Ij為方向向量。
本文采用SSTk-ω湍流模型,湍動能k和湍流耗散率 ω的方程如下:
式中:μ為分子黏度;xj為坐標軸;Uj為在xj方向上的平均速度分量;μt為湍流渦黏度;tij為湍流雷諾應(yīng)力張量;Sij為平均應(yīng)變率張量;F1為輔助混合函數(shù);Pω為ω的導(dǎo)出項;β*,σk和σω2為模型常數(shù)。
2.2.1 計算區(qū)域設(shè)置
采用CFD 進行水流沖擊載荷分析,船舶和RDT 裝置配合三維模型如圖2 所示。由于船舶計算水流速度較低,且船舶和RDT 裝置尺度差異較大。為消除尺度效應(yīng)的影響,采用實尺度完整船舶計算。計算區(qū)域的上下邊界為壓力邊界,遠流場邊界為前后截面及遠離船體一側(cè)。
圖2 船-RDT 裝置水動力CFD 計算域及表面網(wǎng)格劃分情況Fig.2 Computational domain of hydrodynamic CFD and surface network division of the RDT device
2.2.2 網(wǎng)格劃分
采用軟件自帶前處理器自動生成全六面體網(wǎng)格,在自由面網(wǎng)格進行加密;對RDT 裝置封板、翼板和輪緣導(dǎo)流罩進行充分加密,以反映RDT 裝置外形尺寸的影響。由于計算流速較低,進行本船CFD 繞流問題計算分析時將壁面邊界層分為了3 個子區(qū)域,即粘性底層、過渡層和慣性子層。
全部網(wǎng)格帶RDT 裝置,計算網(wǎng)格總數(shù)為580 萬。計算區(qū)域設(shè)置及船體-裝置表面網(wǎng)格劃分如圖2 所示。
為充分考慮RDT 裝置所承受的極端載荷,在計算時,以船舶縱向航行,RDT 裝置軸線與船舶中縱向剖面呈90°夾角時受到的載荷最大。計算水流速度設(shè)置為3 kn,4 kn 和5 kn。分別計算船舶及RDT 裝置受到的水流沖擊載荷及水動力載荷。
RDT 裝置在不同航速(水流)受到的水流沖擊載荷如表1 所示。
表1 RDT 裝置在不同流速下的受力計算Tab.1 Force calculation of the RDT device at different flow rates
可以看出,RDT 裝置在水流速度為5 kn 時候,承受的載荷最大,受到的水流沖擊載荷達到5.07 kN。RDT裝置在水動力作用下,承受總波浪載荷如圖3 所示。
圖3 5 kn 航速時的RDT 裝置總波浪載荷圖Fig.3 Diagram of the total wave load of the RDT device at the speed of 5 kn
可以看出,在航速5 kn 的時候,承受到載荷為3 900 Pa。
裝置除了承受水動力載荷之外,還受到靜壓力載荷,兩者合成為RDT 裝置承受的總的載荷。圖4 為RDT裝置在不同航速下所承受的靜水壓力載荷,可以看出,水動力載荷比靜壓載荷小近20 倍,所以裝置以承受靜水壓力載荷為主,正比于吃水深度,最大壓力載荷為0.08 MPa。
圖4 5 kn 航速時的RDT 裝置水動力載荷圖Fig.4 Hydrodynamic load of the RDT device at the speed of 5 kn
采用Ansys 軟件進行有限元仿真,RDT 裝置采用Q235 鋼裝配制造而成,材料參數(shù)如表2 所示。
表2 RDT 裝置材料屬性Tab.2 Material properties of the RDT device
結(jié)構(gòu)模型采用三維有限元軟件進行建模。由于該RDT 裝置的各部分結(jié)構(gòu)分別采用螺栓或焊接等方式進行緊固。邊界條件采用剛固模式設(shè)置,并進行網(wǎng)格劃分,如圖5 所示。
圖5 RDT 裝置三維有限元模型Fig.5 The finite element model of the RDT device
由載荷特征分析可知在船速為5 kn 時,水動力載荷較靜壓載荷大小相差一個數(shù)量級,因此加載時,僅考慮靜水載荷。靜水載荷與水深成正比,將受力加載分為上、中、下3 個區(qū)域簡化加載,如圖6 所示。
考慮到RDT 在實際使用過程中采用導(dǎo)軌方式進行上下升降,本次分析中對其上部節(jié)點進行約束,如圖7所示。
圖7 邊界條件Fig.7 Boundray condition
RDT 裝置底封板在靜水載荷下的整體變形與應(yīng)力分布如圖8 所示。RDT 裝置的最大變形出現(xiàn)封板遠離支架處,為0.068 mm,且變形分布與靜水下的載荷分別趨于一致;RDT 裝置底部封板最大應(yīng)力出現(xiàn)在RDT 導(dǎo)管與支架接觸處,為4.56 MPa,應(yīng)力主要集中于支撐部件的拐角處,遠小于材料的屈服極限235 MPa。
圖8 RDT 變形與應(yīng)力分布圖Fig.8 The deformation and stress distribution diagram of the RDT device
本文通過研究不同作業(yè)工況下RDT 裝置底封板的載荷成分和載荷特征,提出針對RDT 裝置底封板設(shè)計載荷確定方法,并進一步開展了結(jié)構(gòu)強度評估,主要結(jié)論如下:
1)RDT 裝置的波浪載荷隨航速的增加顯著增大,由輪緣導(dǎo)流罩傳遞到底封板的載荷也會相應(yīng)顯著增加。
2)對于本船RDT 裝置的布置方案,5 級海況下的RDT 裝置輪緣導(dǎo)流罩和底封板區(qū)域的波浪載荷成分占比較少,遠小于靜水載荷,可忽略波浪載荷的影響。
3)RDT 裝置底部封板的最大應(yīng)力出現(xiàn)在RDT導(dǎo)管與支架接觸處,設(shè)計時應(yīng)重點關(guān)注該區(qū)域。