陳昌平,王文,鄭佳成,鄭艷娜,史憲瑩,劉昌鳳、3(.大連海洋大學(xué)海洋與土木工程學(xué)院,遼寧大連603;.馬鞍山市水利局,安徽馬鞍山43000;3.大連理工大學(xué)海岸和近海工程國(guó)家重點(diǎn)實(shí)驗(yàn)室,遼寧大連604)

水流作用下平面金屬網(wǎng)衣水阻力特性數(shù)值模擬

陳昌平1,王文1,鄭佳成2,鄭艷娜1,史憲瑩1,劉昌鳳1、3
(1.大連海洋大學(xué)海洋與土木工程學(xué)院,遼寧大連116023;2.馬鞍山市水利局,安徽馬鞍山243000;3.大連理工大學(xué)海岸和近海工程國(guó)家重點(diǎn)實(shí)驗(yàn)室,遼寧大連116024)

為研究養(yǎng)殖用金屬網(wǎng)衣多孔小直徑網(wǎng)狀結(jié)構(gòu)的水動(dòng)力響應(yīng)特性,采用計(jì)算流體動(dòng)力學(xué)原理,基于Navier-Stokes方程及大渦模擬(LES)湍流模型,并運(yùn)用ABAQUS/CFD模塊中的Gauss-Seidel耦合方法進(jìn)行了水流作用下網(wǎng)衣的流固耦合計(jì)算,在試驗(yàn)驗(yàn)證的基礎(chǔ)上,得到了不同目腳尺寸和網(wǎng)線直徑組合條件下鋅鋁合金網(wǎng)衣的平面受力計(jì)算結(jié)果。結(jié)果表明:當(dāng)網(wǎng)衣目腳尺寸由25 mm增加到35、45 mm時(shí),3種網(wǎng)線直徑對(duì)應(yīng)的網(wǎng)衣受到的水阻力值平均減小幅度分別為29.15%和43.79%;當(dāng)網(wǎng)線直徑由2.5 mm增加到3.2、4.0 mm時(shí),3種目腳尺寸對(duì)應(yīng)的網(wǎng)衣受到的水阻力值平均增加幅度分別為27.01%和55.84%。本研究結(jié)果可為進(jìn)一步探究金屬網(wǎng)箱在波流作用下的整體受力變形提供理論參考。

金屬網(wǎng)衣;流固耦合;水阻力特性;數(shù)值模擬

在近海養(yǎng)殖海域生態(tài)環(huán)境不斷惡化、人類對(duì)食品質(zhì)量安全要求不斷提高的背景下,近岸傳統(tǒng)網(wǎng)箱養(yǎng)殖的弊端逐漸顯現(xiàn)。近年來(lái),深水金屬網(wǎng)箱逐漸進(jìn)入中國(guó)網(wǎng)箱養(yǎng)殖業(yè)的視野,與重力式網(wǎng)箱相比,金屬網(wǎng)箱具有抗浪流、耐腐蝕、輕附著、低污染、養(yǎng)殖空間大和水產(chǎn)品品質(zhì)高等優(yōu)勢(shì),將會(huì)成為中國(guó)極具發(fā)展?jié)摿Φ男滦宛B(yǎng)殖模式。網(wǎng)衣是金屬網(wǎng)箱的重要組成部分,在養(yǎng)殖環(huán)境中承受著復(fù)雜的波浪水流作用,在網(wǎng)箱養(yǎng)殖向外海深水推進(jìn)過(guò)程中其安全性備受關(guān)注,因此,研究金屬網(wǎng)衣在波浪水流作用下的水動(dòng)力特性具有十分重要的意義。

多年來(lái),國(guó)內(nèi)外學(xué)者通過(guò)模型試驗(yàn)、數(shù)值模擬手段針對(duì)纖維類網(wǎng)衣開展了系列研究,主要包括網(wǎng)衣的模型試驗(yàn)相似準(zhǔn)則研究[1-4]、網(wǎng)衣水阻力系數(shù)研究[5-9]、網(wǎng)衣在波流作用下的動(dòng)力特性研究[10-13]等。在此基礎(chǔ)上,相關(guān)學(xué)者開展了一系列由網(wǎng)衣、浮架和錨碇系統(tǒng)構(gòu)成的深水網(wǎng)箱水動(dòng)力特性研究。近年來(lái),隨著金屬網(wǎng)箱的逐漸運(yùn)用,金屬網(wǎng)衣的研究也隨之展開。Tsukrov等[14]以網(wǎng)衣密實(shí)度為參數(shù),對(duì)4種不同密實(shí)度的銅合金平面網(wǎng)衣在水流作用下的受力進(jìn)行了測(cè)試,分析了4種網(wǎng)衣的阻力系數(shù),并與兩種不同密實(shí)度的尼龍網(wǎng)衣測(cè)試結(jié)果進(jìn)行了比較,結(jié)果表明,銅合金網(wǎng)衣的法向阻力系數(shù)明顯小于尼龍網(wǎng)衣,同時(shí)討論了影響網(wǎng)衣法向阻力系數(shù)的相關(guān)因子(網(wǎng)目形狀、網(wǎng)衣密實(shí)度、網(wǎng)衣材質(zhì))。Gansel等[15]采用試驗(yàn)測(cè)試的方法,得到了水流作用下不同孔隙率的圓柱形金屬網(wǎng)衣流阻力值,并與經(jīng)驗(yàn)公式進(jìn)行了比較。Cha等[16]采用試驗(yàn)的方法研究了不同直徑、網(wǎng)目長(zhǎng)度在不同沖角條件下的鏈形銅合金網(wǎng)衣水阻力系數(shù),并和柔性網(wǎng)衣的阻力系數(shù)進(jìn)行了比較。Decew等[17]通過(guò)試驗(yàn)方法獲得了一種鏈形銅合金網(wǎng)衣的等效剛度,并證實(shí)采用傳統(tǒng)的有限元方法可有效模擬該種網(wǎng)衣組成的重力式網(wǎng)箱的體積變形。

考慮到纖維類網(wǎng)衣屬于小直徑多孔柔性結(jié)構(gòu),在網(wǎng)衣的數(shù)值模擬中,通常采用的網(wǎng)衣單元模型主要包括桁架單元、梁?jiǎn)卧?、索單元、彈簧單元等。本研究中針?duì)一種鋅鋁合金網(wǎng)衣,依據(jù)有限元原理,基于ABAQUS軟件,將網(wǎng)衣視為梁?jiǎn)卧Ｐ?并首次運(yùn)用Gauss-Seidel耦合方案,在試驗(yàn)結(jié)果驗(yàn)證的基礎(chǔ)上,計(jì)算了方形金屬網(wǎng)衣在水流作用下的水阻力值,討論了網(wǎng)衣受力與目腳尺寸、網(wǎng)線直徑的關(guān)系,數(shù)值模擬結(jié)果將為進(jìn)一步研究金屬網(wǎng)箱在波流作用下的整體受力變形提供有益參考。

1 數(shù)值模擬方法

1.1 流場(chǎng)模擬

模擬流體為牛頓流體,連續(xù)性方程為

針對(duì)不可壓縮流體,連續(xù)性方程為

其中:t為時(shí)間;v和籽分別為流體的速度和密度。

不可壓縮流體的動(dòng)量守恒N-S方程為

其中:p為流體壓力;滋為流體運(yùn)動(dòng)黏度系數(shù);g為重力加速度。

為便于數(shù)值求解N-S方程,采用Johnston等[18]提出的與不可壓縮流N-S方程等價(jià)的壓力泊松方程(PPE)形式:

式中:f為外力項(xiàng)。

對(duì)公式(4)、(5)通過(guò)下式進(jìn)行離散:

其中:M為質(zhì)量矩陣;A(v)為對(duì)流算子;K為黏性擴(kuò)散算子;G為梯度算子;D為散度算子;L為壓力泊松算子;F為體力。

流場(chǎng)中的湍流模型采用大渦模擬法(LES),大尺度渦通過(guò)濾波后的Navier-Stokes方程直接求解,小尺度的渦則采用亞格子尺度模型(SGS)模擬。

經(jīng)過(guò)濾波函數(shù)的濾波,可得到大尺度渦的不可壓縮Navier-Stokes方程為

亞格子應(yīng)力子ij采用Smagorinsky-Lilly的亞格子尺度模型,具有以下形式:

式中:淄t為亞格子尺度的湍流黏度,且

式中:駐i為沿坐標(biāo)軸i方向的網(wǎng)格尺寸;CS為Sma-gorinsky常數(shù),且

式中:y+為網(wǎng)格到壁面的最近距離;A+為半經(jīng)驗(yàn)常數(shù),取25;CS0為VanDriest常數(shù),取0.1。

本研究中,流場(chǎng)模擬區(qū)域大小為:X方向1200 mm,Y方向1000 mm,Z方向1030 mm。

1.2 網(wǎng)衣模擬

在ABAQUS的Standard模塊里,按照網(wǎng)衣平面尺寸及網(wǎng)線直徑、網(wǎng)目大小,建立了一個(gè)與實(shí)際模型相同的網(wǎng)衣模型,并將建立好的模型導(dǎo)入到Hypermesh分析軟件中。劃分網(wǎng)格時(shí),將導(dǎo)入Hy-permesh軟件的網(wǎng)衣進(jìn)行十字經(jīng)緯交叉。每個(gè)經(jīng)緯十字交叉時(shí),上下左右各取半個(gè)網(wǎng)衣長(zhǎng)度大小的距離,將所取的部分看成一個(gè)獨(dú)立體,整張網(wǎng)衣就是由許多此類獨(dú)立體形成的矩形陣列。網(wǎng)格采用結(jié)構(gòu)化網(wǎng)格,所有單元都由六面體組成。網(wǎng)格的收斂性由ABAQUS分析無(wú)誤后方可完成整個(gè)運(yùn)算過(guò)程。以直徑4 mm、網(wǎng)目45 mm的平面網(wǎng)衣為例,網(wǎng)衣劃分的總單元個(gè)數(shù)為2 280 944個(gè)。圖1為網(wǎng)衣結(jié)構(gòu)單元?jiǎng)澐质疽鈭D。

網(wǎng)衣單元結(jié)構(gòu)運(yùn)動(dòng)方程為

其中:m為質(zhì)量矩陣;c為阻尼矩陣;k為剛度矩陣;x為結(jié)構(gòu)位移;F為作用在單元結(jié)構(gòu)上的外力。

1.3 流固耦合計(jì)算

在ABAQUS中,使用ABAQUS/Standard和ABAQUS/CFD兩個(gè)模塊進(jìn)行流固耦合計(jì)算,流體域和固體域的計(jì)算分開進(jìn)行,在流體域內(nèi)釆用有限體積法和有限單元法結(jié)合,對(duì)線性化的基于壓力的流體動(dòng)力學(xué)控制方程及湍流模型控制方程進(jìn)行迭代求解,得到結(jié)構(gòu)表面流體的速度,通過(guò)伯努利方程,求得結(jié)構(gòu)表面的壓強(qiáng),沿結(jié)構(gòu)表面積分可得到作用于結(jié)構(gòu)的流體水動(dòng)力荷載值;在結(jié)構(gòu)域內(nèi)采用有限單元法計(jì)算結(jié)構(gòu)響應(yīng)。流固之間將每一步計(jì)算的物理量相互反饋給對(duì)方,并得到對(duì)方所得出的物理量進(jìn)行下一步計(jì)算。圖2為Gauss-Seidel流固耦合計(jì)算過(guò)程示意圖。

圖1 網(wǎng)衣單元結(jié)構(gòu)網(wǎng)格示意圖Fig.1M esh of the net structure element

圖2 流固耦合計(jì)算示意圖Fig.2Schematic diagram of Gauss-Seidelcoup ling schem e

2 數(shù)值模擬驗(yàn)證

2.1 水槽試驗(yàn)

銀行規(guī)模與向小企業(yè)貸款力度呈反向關(guān)系，廣西貧困縣縣域經(jīng)濟(jì)主要以中小企業(yè)為主。對(duì)縣域中小企業(yè)提供足夠資金支持，就會(huì)對(duì)縣域經(jīng)濟(jì)第二產(chǎn)業(yè)發(fā)展產(chǎn)生顯著的促進(jìn)作用。發(fā)展本地小商業(yè)銀行，有利于向本地縣域中小企業(yè)貸款，服務(wù)本地經(jīng)濟(jì)。

為驗(yàn)證數(shù)值模擬的準(zhǔn)確性和可靠性,在大連理工大學(xué)船模試驗(yàn)水槽中進(jìn)行了金屬網(wǎng)衣試驗(yàn),水槽規(guī)格為170 m伊7 m伊4 m(長(zhǎng)伊寬伊深)。根據(jù)目前網(wǎng)箱養(yǎng)殖工程實(shí)際情況,確定鋅鋁金屬網(wǎng)衣的3種網(wǎng)目和3種網(wǎng)線直徑,網(wǎng)衣具體參數(shù)如表1所示。

表1 平面網(wǎng)衣參數(shù)Tab.1Properties of a plane net

試驗(yàn)過(guò)程中,網(wǎng)衣固定在定制的鐵質(zhì)框架內(nèi),網(wǎng)衣四邊與框架固定,框架上部為傳感器和數(shù)據(jù)采集系統(tǒng)。圖3為網(wǎng)衣試驗(yàn)布置圖。網(wǎng)衣上邊緣距離水表面為15 cm。試驗(yàn)采用拖曳方法,拖車速度分別為0.2、0.4、0.6、0.8、1.0 m/s。由于采用拖曳方式,框架的存在對(duì)流場(chǎng)的影響較小,因此,可以忽略對(duì)網(wǎng)衣受力的影響。

圖3 網(wǎng)衣試驗(yàn)布置圖Fig.3Experimental set up of net

2.2 結(jié)果比較

根據(jù)網(wǎng)線直徑、目腳尺寸和流速分組情況,網(wǎng)衣試驗(yàn)工況共有30種,考慮不失一般性,圖4給出了鋅鋁合金網(wǎng)衣在5種流速條件、不同網(wǎng)目大小及網(wǎng)線直徑情況下,網(wǎng)衣受力的數(shù)值計(jì)算結(jié)果與試驗(yàn)結(jié)果的對(duì)比情況。計(jì)算公式為

網(wǎng)衣受力試驗(yàn)值=框架與網(wǎng)衣總體受力試驗(yàn)值-單獨(dú)框架拖曳受力試驗(yàn)值。

圖4-A中25 mm伊2.5 mm表示網(wǎng)衣目腳尺寸為25 mm、網(wǎng)線直徑為2.5 mm的網(wǎng)衣;圖4-B中35 mm伊3.2 mm表示目腳尺寸為35 mm、網(wǎng)線直徑為3.2 mm的網(wǎng)衣;圖4-C中45 mm伊4.0 mm表示目腳尺寸為45 mm、網(wǎng)線直徑為4.0 mm的網(wǎng)衣。

由圖4可知,不同工況下網(wǎng)衣受力的數(shù)值模擬計(jì)算結(jié)果與試驗(yàn)結(jié)果吻合較好,3種工況下計(jì)算值與試驗(yàn)值網(wǎng)衣受力值相對(duì)誤差分別為7.59%、4.35%和4.65%,統(tǒng)計(jì)其他各種工況計(jì)算值與試驗(yàn)值相對(duì)誤差,平均值小于10%,說(shuō)明采用該數(shù)值計(jì)算方法具有較好的準(zhǔn)確性。

3 討論

3.1 網(wǎng)衣受力與目腳尺寸的關(guān)系

由表2可知,在3種網(wǎng)線直徑條件下,目腳尺寸為35、45 mm時(shí),各工況下網(wǎng)衣受力平均減小幅度分別為29.15%和43.79%。根據(jù)Morsion等提出的計(jì)算小直徑構(gòu)件流阻力半經(jīng)驗(yàn)半理論方法可知,在純流作用下,平面網(wǎng)衣受力大小與網(wǎng)衣在水流作用方向的投影面積大小成正比關(guān)系,在平面方形編織網(wǎng)外框尺寸及網(wǎng)衣網(wǎng)線直徑相同情況下,目腳尺寸增加,網(wǎng)衣在水流作用方向的投影面積減小,因此,網(wǎng)衣整體受力也隨之減小。

圖4 網(wǎng)衣受力計(jì)算值與試驗(yàn)值比較Fig.4Comparison of the calculation value and meas-ured results of force acting on the net

圖5 不同目腳尺寸條件下網(wǎng)衣受力比較Fig.5Comparison of net forces under different mesh sizes

表2 目腳尺寸變化引起的網(wǎng)衣受力變化幅度統(tǒng)計(jì)Tab.2Statistics of relative variation in the net force in differentmesh sizes

3.2 網(wǎng)衣受力與網(wǎng)線直徑的關(guān)系

圖6為網(wǎng)衣目腳尺寸相同、網(wǎng)線直徑不同條件下,鋅鋁合金網(wǎng)衣在各種流速條件下的受力比較,其中l(wèi)=25 mm、l=35 mm、l=45 mm分析表示目腳尺寸為25、35、45 mm。由圖6可知,在目腳尺寸相同條件下,隨網(wǎng)線直徑的增加,網(wǎng)衣受力也隨之增加。以網(wǎng)線直徑2.5 mm為基準(zhǔn),在不同流速作用下,直徑為3.2、4.0 mm時(shí)網(wǎng)衣受力增加幅度平均值見表3。表3中3.2 mm/2.5mm表示網(wǎng)線直徑3.2 mm相對(duì)于網(wǎng)線直徑2.5 mm網(wǎng)衣受力的增加幅度。

由表3可知,在3種網(wǎng)目大小情況下,網(wǎng)線直徑從2.5 mm增加到3.2、4.0 mm時(shí),在不同流速作用下,網(wǎng)衣受力平均增加幅度分別為27.01%和55.84%。由于網(wǎng)線直徑增加,在目腳尺寸不變情況下,網(wǎng)衣在水流方向的投影面積相應(yīng)增加,因此,網(wǎng)衣受到的流阻力也相應(yīng)有所增加。

圖6 不同網(wǎng)線直徑條件下網(wǎng)衣受力比較Fig.6Comparison of net forces under different diame-ters of tw ines

表3 網(wǎng)線直徑變化引起的受力變化幅度統(tǒng)計(jì)Tab.3Statistics of relative variation in the net force indifferent diameters of tw ines

4 結(jié)論

基于Navier-Stokes方程及湍流模型,運(yùn)用CFD流固耦合技術(shù),得到了金屬網(wǎng)衣在水流作用下的阻力值,一種鋅鋁合金網(wǎng)衣在水流作用下,水阻力大小與網(wǎng)目尺寸及網(wǎng)線直徑的關(guān)系為:

(1)當(dāng)目腳尺寸由25 mm增加到35、45 mm時(shí),網(wǎng)線直徑對(duì)應(yīng)的網(wǎng)衣受到的水阻力值平均減小幅度分別為29.15%和43.79%。

(2)當(dāng)網(wǎng)線直徑由2.5 mm增加到3.2、4.0 mm時(shí),目腳尺寸對(duì)應(yīng)的網(wǎng)衣受到的水阻力值平均增加幅度分別為27.01%和55.84%。

[1]Fredriksson D W.Open ocean fish cage and mooring system dy-namics[D].New Hampshire:University of New Hampshire,2001.

[2]Lader P F,Fredheim A,Lien E.Dynamic behaviour of3D nets ex-posed to waves and current[C]//Proceedings of the 20th Interna-tional Conference on Offshore Mechanics and Arctic Engineering. Rio de Janeiro,Brazil:OMAE,2001.

[3]桂福坤,李玉成,張懷慧.網(wǎng)衣受力試驗(yàn)的模型相似條件[J].中國(guó)海洋平臺(tái),2002,17(5):22-25.

[4]李玉成,桂福坤,張懷慧,等.深水養(yǎng)殖網(wǎng)箱試驗(yàn)中網(wǎng)衣相似準(zhǔn)則的應(yīng)用[J].中國(guó)水產(chǎn)科學(xué),2005,12(2):179-187.

[5]Swift M R,Fredriksson DW,Unrein A,et al.Drag force acting on biofouled net panels[J].Aquacultural Engineering,2006,35(3): 292-299.

[6]Zhan JM,Jia X P,LiY S,etal.Analytical and experimental inves-tigation of drag on nets of fish cages[J].Aquacultural Engineer-ing,2006,35(1):91-101.

[7]Zhan Jiemin,Su Wei,Li Yuxiang,et al.Comparison of drags on fish cages of different shapes[J].China Ocean Engineering,2006, 20(3):457-471.

[8]SongWeihua,Liang Zhenlin,ChiHongfu,et al.Experimental study on the effect of horizontal waves on netting panels[J].Fisheries Science,2006,72(5):967-976.

[9]李玉成,桂福坤.平面無(wú)結(jié)節(jié)網(wǎng)衣水阻力系數(shù)的試驗(yàn)研究[J].海洋學(xué)報(bào),2006,28(5):145-151.

[10]Shimizu T,Takagi T,Suzuki K,et al.Refined calculation model for NaLA,a fishing net shape simulator,applicable to gill nets [J].Fisheries Science,2004,70(3):401-411.

[11]Tsukrov I,Eroshkin O,Fredriksson D,etal.Finite elementmodel-ing of net panels using a consistent net element[J].Ocean Engi-neering,2003,30(2):251-270.

[12]Takagi T,Suzuki K,Hiraishi T.Development of the numerical simulation method of dynamic fishing net shape[J].Nippon Su-isan Gakkaishi,2002,68(3):320-326.

[13]Lader P F,Fredheim A.Dynamic properties of a flexible net sheet in waves and current-a numerical approach[J].Aquacultural Engineering,2006,35(3):228-238.

[14]Tsukrov I,Drach A,DeCew J,et al.Characterization of geometry and normal drag coefficients of copper nets[J].Ocean Engineer-ing,2011,38(17-18):1979-1988.

[15]Gansel L C,McClimans T A,Myrhaug D.The effects of fish cages on ambient currents[J].Journal of Offshore Mechanics and Arc-tic Engineering,2012,134(1):011303.

[16]Cha B J,Kim H Y,Bae JH,et al.Analysis of the hydrodynamic characteristics of chain-link woven copper alloy nets for fish ca-ges[J].Aquacultural Engineering,2013,56:79-85.

[17]Decew JC,OsienskiM,Drach A,etal.Effectof the bending stiff-ness on the volumetric stability of fish cages with copper alloy netting[C]//Proceedings of the 5th International Conference on Computational Methods in Marine Engineering.Spain:MARINE, 2013:1-10.

[18]Johnston H,Liu Jianguo.Accurate,stable and efficient Navier-Stokes solvers based on explicit treatment of the pressure term [J].Journal of Computational Physics,2004,199(1):221-259.

Numerical simulation on hydrodynam ic characteristics of the planemetal net under current

CHEN Chang-ping1,WANGWen1,ZHENG Jia-cheng2,ZHENG Yan-na1, SHIXian-ying1,LIU Chang-feng1,3
(1.College of Marine and Civil Engineering,Dalian Ocean University,Dalian 116023,China;2.Water Conservancy Bureau of Ma蒺anshan City, Ma蒺anshan 243000,China;3.State Key laboratory of Coastal and Offshore Engineering,Dalian University of Technology,Dalian 116024,China)

In order to study the hydrodynamic characteristics of the metal net in aquaculture,which is a special kind of porous net structurewith small diameter,the computational fluid dynamicswas used in the force calculation of the plane net.Based on the Navier-Stokes Equation and the turbulencemodel of Large Eddy Simulation(LES), the fluid-structure interaction was simulated by the Gauss-Seidelmethod in the ABAQUS/CFDmodel.Verified by the experiment,the force on the plane netwith different net size and different thread diameter was computed and analyzed.The simulation results showed thatwhen the net size increased from 25 mm to 35 mm,and 45 mm,the drag forces of the net with three different diameter thread were decreased by 29.15%and 43.79%respectively. When the thread diameter increases from 2.5 mm to 3.2 mm,and 4.0 mm,the average increase of drag forces of the netwith three different net size were increased by 27.01%and 55.84%respectively.The results provide a theoretical reference for the further study of themetal cage under the action of wave current.

metal net;fluid-solid coupling;hydrodynamic characteristics;mathematical simulation

TV131.6

A

10.16535/j.cnki.dlhyxb.2017.03.017

2095-1388(2017)03-0363-06

2016-12-30

國(guó)家自然科學(xué)基金資助項(xiàng)目(31572663,51239002,51409039);遼寧省教育廳科研項(xiàng)目(L2014280)

陳昌平(1971—),男,博士,教授。E-mail:ccp@dlou.edu.cn

鄭艷娜(1978—),女,博士,副教授。E-mail:zhengyn@dlou.edu.cn