劉航飛, 陳昌平, 鄭艷娜, 梁新宇
(大連海洋大學(xué)海洋與土木工程學(xué)院,遼寧 大連 116023)
水流作用下一種養(yǎng)殖金屬網(wǎng)衣水阻力特性的數(shù)值模擬研究
劉航飛, 陳昌平, 鄭艷娜, 梁新宇
(大連海洋大學(xué)海洋與土木工程學(xué)院,遼寧 大連 116023)
為保障深遠(yuǎn)海金屬網(wǎng)箱養(yǎng)殖系統(tǒng)的結(jié)構(gòu)安全,進行網(wǎng)箱結(jié)構(gòu)的水動力學(xué)研究是一項重要的工作。構(gòu)成金屬網(wǎng)箱主要組成部分的金屬網(wǎng)衣是一種小直徑多孔的網(wǎng)狀結(jié)構(gòu),在波浪和水流作用下的水動力學(xué)特性與一般的海洋工程結(jié)構(gòu)物有顯著不同。本文基于有限元基本原理,采用梁單元模擬金屬網(wǎng)線結(jié)構(gòu),采用連接單元模擬網(wǎng)線接觸部位的相互作用,并運用ABAQUS軟件進行數(shù)值計算,在試驗驗證的基礎(chǔ)上,分析了在水流作用下一種金屬菱形鏈網(wǎng)衣在不同網(wǎng)目尺寸、不同網(wǎng)線直徑情況下的水阻力變化情況。數(shù)值模擬結(jié)果表明在各種工況下,當(dāng)網(wǎng)目尺寸由25 mm增加到35 mm、45 mm時,網(wǎng)衣受力增加幅度的平均值分別為14.71%和38.07%;當(dāng)網(wǎng)線直徑由2.5 mm增加到3.2 mm、4.0 mm時,各工況下網(wǎng)衣受力增加幅度的平均值分別為25.05%和45.06%。研究結(jié)果可以為進一步開展深海養(yǎng)殖金屬網(wǎng)箱水動力特性研究提供積極的基礎(chǔ)。
金屬網(wǎng)衣;水阻力特性;梁單元;連接單元;數(shù)值模擬
近年來近岸網(wǎng)箱養(yǎng)殖業(yè)對海洋環(huán)境的破壞日益嚴(yán)重,且采用的纖維類網(wǎng)箱在波流作用下網(wǎng)衣變形嚴(yán)重,影響了魚類生存空間,導(dǎo)致魚類品質(zhì)降低。深水金屬網(wǎng)箱以其抗浪、抗流、耐腐蝕、輕附著、低污染等優(yōu)勢成為極具發(fā)展?jié)摿Φ男滦宛B(yǎng)殖模式。金屬網(wǎng)衣是網(wǎng)箱的重要組成部分,在工作過程中,承受復(fù)雜的波浪水流作用,不僅是整個網(wǎng)箱系統(tǒng)受力的主要部分,也是整個網(wǎng)箱系統(tǒng)變形較大,容易破壞的部位,網(wǎng)衣的安全性已成為網(wǎng)箱離岸養(yǎng)殖過程中首要關(guān)注的問題。研究金屬網(wǎng)衣的水阻力特性,為網(wǎng)衣的安全設(shè)計提供設(shè)計依據(jù)是一項重要的基礎(chǔ)性工作。
多年來各國學(xué)者采用數(shù)值模擬和模型試驗的方法對纖維類網(wǎng)衣的水動力特性進行了研究:采用有限差分法[1]和有限單元法[2-8]對網(wǎng)衣進行了數(shù)值模擬研究;采用試驗方法對網(wǎng)衣的水阻力系數(shù)、變形和受力等水動力特性進行了研究[9-11],采用彈簧-集中質(zhì)量模型[12-16]和多孔介質(zhì)模型[17]分析了網(wǎng)衣的水流阻力和變形。在金屬網(wǎng)衣研究方面,Drach A等[18-19]基于有限元方法和接觸算法,分析了一種金屬菱形鏈網(wǎng)在靜力加載作用下連接部位的應(yīng)力情況。陳昌平等[20]基于CFD原理,建立了流體與金屬編織網(wǎng)衣之間的流固耦合模型,分析了一種金屬網(wǎng)衣在水流作用下的受力特性。盡管有許多學(xué)者對網(wǎng)衣水動力特性進行了較多研究,但主要針對于平面矩形或菱形網(wǎng)目的纖維類網(wǎng)衣,而對于具有三維空間結(jié)構(gòu)的金屬菱形鏈網(wǎng)的研究相對較少。
深遠(yuǎn)海養(yǎng)殖用金屬菱形鏈網(wǎng)是一種小尺度多孔彈性網(wǎng)狀結(jié)構(gòu),同時網(wǎng)衣在材質(zhì)、網(wǎng)目形狀、編制方式等方面與傳統(tǒng)纖維類柔性網(wǎng)衣相比,具有網(wǎng)衣整體剛度大、網(wǎng)線連接方式特殊等特點,因此在波浪水流作用下的受力和變形特性與纖維類網(wǎng)衣有顯著的不同。文中采用有限單元法,針對網(wǎng)線結(jié)構(gòu)采用梁單元進行模擬,采用連接單元模擬網(wǎng)線交接處的接觸,利用ABAQUS軟件對一種鋅鋁合金菱形鏈網(wǎng)在水流作用下的水阻力特性進行了數(shù)值計算,并在試驗驗證的基礎(chǔ)上分析了金屬網(wǎng)衣水阻力大小與網(wǎng)目尺寸以及網(wǎng)線直徑之間的關(guān)系,研究結(jié)果將為進一步開展整體網(wǎng)箱系統(tǒng)的水動力特性研究打下良好的數(shù)?;A(chǔ)。
金屬菱形鏈網(wǎng)具有較為復(fù)雜的三維空間結(jié)構(gòu)。圖1為網(wǎng)目大小m=45 mm,直徑大小d=4.0 mm的鋅鋁合金菱形鏈網(wǎng)模型圖。在網(wǎng)衣結(jié)構(gòu)建模時,將其投影截面近似為橢圓,長軸a=3.5 cm,短軸b=1.8 cm。然后按照螺距s=6.7 cm進行拉伸,形成一個總長為0.8 m長的鏈條。最后將鏈條進行陣列,組合成一個長寬為0.8 m×0.8 m的金屬網(wǎng)衣。在有限元軟件ABAQUS中建立的模型尺寸采用毫米制。
圖1 金屬菱形鏈網(wǎng)模型圖
金屬網(wǎng)衣采用兩節(jié)點的梁單元進行模擬,局部網(wǎng)格的劃分見圖2,針對圖1所示網(wǎng)衣模型,劃分后的節(jié)點總數(shù)為6 358個,單元總數(shù)為6 336個。在兩根網(wǎng)線交接部位加入連接單元,模擬網(wǎng)線之間的接觸作用,計算時約束接觸處3個位移方向的相對自由度,不約束相對旋轉(zhuǎn)自由度。
圖2 網(wǎng)衣單元劃分示意圖
1.3.1 受力分析
(1)重力與浮力。由于網(wǎng)衣全部浸沒在水面以下,因此將重力和浮力一起考慮。圓形截面桿件的重力和浮力的計算公式[21]為:
(1)
式中:PZ為桿件所受重力和浮力的合力,N;ρf和ρw分別為金屬網(wǎng)衣的密度和水的密度,kg/m3;D為桿件的橫截面直徑,m;l為桿件長度,m;V為桿件體積,m3。
(2)水流阻力。浸沒在水中的網(wǎng)衣受到的水流阻力采用Morison公式[22]計算。為了便于計算,將水流阻力沿桿件的法向和切向分解為法向水流阻力和切向水流阻力,如圖3所示。
(2)
式中:Fn為單位長度的法向阻力,N;Cn為法向阻力系數(shù);Δvn為流體相對于桿的速度在桿法面內(nèi)的分量,m/s。
(3)
式中:Ft為單位長度的切向阻力,N;Ct為切向的摩擦力系數(shù);Δvt為流體相對于桿的速度在桿切面內(nèi)的分量,m/s。
圖3 流體速度沿桿的切向和法向分量
采用Morison公式計算小尺度構(gòu)件在水流作用下的受力時,根據(jù)網(wǎng)衣的結(jié)構(gòu)特點及相關(guān)文獻,本文采用Choo[23]等提出的速度力系數(shù)與雷諾數(shù)的關(guān)系式來計算法向和切向水阻力系數(shù)。
(4)
(5)
式中:s=0.077 215 665+ln(8/Ren);Ren=ρwVnD/μ;μ為水黏性系數(shù),kg/(m·s)。
1.3.2 運動方程的構(gòu)建和求解
網(wǎng)衣采用兩節(jié)點的梁單元進行模擬。在三維空間下,每一個節(jié)點有3個位移自由度,其離散化運動方程[24]:
(6)
(7)
式中:Me、Ce、Ke、Qe分別為單元節(jié)點的質(zhì)量矩陣(kg)、阻尼矩陣(N·s/m)、剛度矩陣(N/m)和載荷向量(N);N為插值函數(shù)矩陣;B為應(yīng)變矩陣;D為彈性矩陣;f為由單元體積力組成的節(jié)點荷載矩陣,N;T為單元表面力組成的節(jié)點荷載矩陣,N;Ve為單元體積(m3);Se為單元面積(m2);ρ為質(zhì)量密度(kg/m3);η為阻尼系數(shù)(N·s/m)。圖4為動力隱式算法流程圖。
圖4 動力隱式算法流程圖
網(wǎng)線交接處采用連接單元模擬,連接單元的兩個節(jié)點分別與交接處兩個梁單元節(jié)點相重合,因此連接單元節(jié)點控制方程與梁單元節(jié)點控制方程相同。離散化運動方程采用ABAQUS/Standard直接積分法中的動力隱式算法進行求解。
本試驗采用原型的鋅鋁合金菱形鏈網(wǎng),依據(jù)目前實際工程采用的網(wǎng)衣情況,結(jié)合研究內(nèi)容,試驗選用方形網(wǎng)衣,平面尺寸為80 cm×80 cm,網(wǎng)線直徑分別為2.5 mm、3.2 mm、4.0 mm,網(wǎng)目尺寸分別為25 mm、35 mm、45 mm,金屬網(wǎng)絲密度為6 300 kg/m3,彈性模量193 Gpa,泊松比為0.3。結(jié)合網(wǎng)線直徑及網(wǎng)目尺寸大小,共進行了9組網(wǎng)衣試驗。試驗過程中,網(wǎng)衣固定在定制的鐵質(zhì)框架上,金屬網(wǎng)衣四邊與框架連接,框架上部為傳感器和數(shù)據(jù)采集系統(tǒng),網(wǎng)衣距離水面15 cm,如圖5所示。試驗采用拖拽方法,拖車速度分別為0.2 m/s、0.4 m/s、0.6 m/s、0.8 m/s、1.0 m/s。網(wǎng)衣受力值為框架與網(wǎng)衣總受力試驗值減去單獨框架拖拽試驗受力值。
圖5 網(wǎng)衣實驗布置圖
根據(jù)網(wǎng)衣試驗的邊界條件及試驗工況,采用ABAQUS中動力隱式算法進行網(wǎng)衣水阻力數(shù)值計算,并將計算結(jié)果與試驗結(jié)果進行比較。圖6為五種流速條件下數(shù)值計算結(jié)果與試驗結(jié)果的比較。
由圖可知,各種工況下數(shù)值計算結(jié)果與試驗結(jié)果隨流速的變化均表現(xiàn)出相同的變化趨勢,計算值與試驗值擬合程度良好。以試驗值為基準(zhǔn),統(tǒng)計5種流速條件下各工況下計算值與試驗值之間的相對誤差見表1。表中每一種網(wǎng)衣相對誤差表示在5種流速條件下得到的網(wǎng)衣受力的計算值與試驗值相對誤差平均值,平均相對誤差表示所有工況下網(wǎng)衣受力的計算值與試驗值相對誤差平均值。由表1可知,各種工況下網(wǎng)衣受力計算值與試驗值的相對誤差最大值為10.69%,最小值為5.23%,相對誤差平均值為8.13%。由圖6及表1可以表明,采用本文建立的網(wǎng)衣水阻力計算模型,可以很好地計算金屬菱形鏈網(wǎng)在水流作用下受到水阻力大小,與試驗結(jié)果誤差較小,可以為進一步進行網(wǎng)箱整體受力數(shù)值計算提供參考。
圖6 網(wǎng)衣受力試驗值與模擬值對比圖
表1 水流阻力相對誤差表
由計算結(jié)果可知,在流速與網(wǎng)線直徑相同,網(wǎng)目尺寸不同條件下,金屬菱形鏈網(wǎng)受到的水阻力有較明顯差異。圖7為在網(wǎng)線直徑分別為2.5 mm、3.2 mm和4.0 mm情況下,對應(yīng)網(wǎng)目尺寸分別為25 mm、35 mm、45 mm的網(wǎng)衣在不同流速條件下所受的水阻力比較。
由圖7可知,在相同網(wǎng)線直徑與相同流速條件下,網(wǎng)衣水阻力均隨著網(wǎng)目尺寸的增加有較明顯的減小,當(dāng)網(wǎng)目為25 mm時,網(wǎng)衣受力最大;網(wǎng)目為35 mm次之;網(wǎng)目為45 mm時受力最小。分別以較小網(wǎng)目尺寸網(wǎng)衣受力值為基準(zhǔn),分析各工況下相對較大網(wǎng)目尺寸的網(wǎng)衣受力減小幅度平均值,統(tǒng)計情況見表2,表中35 mm/25 mm表示網(wǎng)目尺寸為35mm的網(wǎng)衣相對于網(wǎng)目尺寸為25 mm的網(wǎng)衣受力減小幅度,其他表示含義類同。
圖7 網(wǎng)衣水阻力隨網(wǎng)目尺寸的變化
表2 網(wǎng)目尺寸變化引起的網(wǎng)衣水阻力變化幅度統(tǒng)計
由表2可知,在3種網(wǎng)線直徑條件下,各工況下3組對比網(wǎng)衣的受力值平均減小幅度分別為14.71%,27.39%,38.07%。不同網(wǎng)線直徑情況下,隨網(wǎng)目尺寸的增加,網(wǎng)衣受力減小幅度值差異不明顯。由于在水流作用下,平面網(wǎng)衣受力大小與網(wǎng)衣在水流作用方向的投影面積成正比例關(guān)系,因此當(dāng)菱形鏈網(wǎng)平面尺寸及網(wǎng)衣網(wǎng)線直徑相同情況下,網(wǎng)目尺寸增加,網(wǎng)衣在水流作用方向的投影面積減小,因此網(wǎng)衣整體受力也隨之減小。
由計算結(jié)果可知,在流速與網(wǎng)目尺寸相同,網(wǎng)線直徑不同條件下,金屬菱形鏈網(wǎng)受到的水阻力有較明顯差異。圖8為在網(wǎng)目分別為25 mm、35 mm、45 mm情況下,對應(yīng)網(wǎng)線直徑分別為2.5 mm、3.2 mm和4.0 mm的網(wǎng)衣在不同流速條件下網(wǎng)衣所受的水阻力比較。
由圖8可知,在網(wǎng)目尺寸和流速相同條件下,網(wǎng)線直徑增加,其所受的水流阻力隨之有明顯增加。各種工況下3種網(wǎng)目尺寸的網(wǎng)衣在直徑為4.0 mm時受到的水阻力最大,直徑為3.2 mm的網(wǎng)衣次之,網(wǎng)衣直徑為2.5 mm時受到的水阻力最小。分別以較小網(wǎng)線直徑的網(wǎng)衣受力值為基準(zhǔn),分析各工況下相對較大網(wǎng)線直徑的網(wǎng)衣受力增加幅度平均值,統(tǒng)計情況見表3,表中3.2 mm/2.5 mm表示網(wǎng)線直徑為3.2 mm的網(wǎng)衣相對于網(wǎng)線直徑為2.5 mm的網(wǎng)衣受力增加幅度,其他表示含義類同。
圖8 網(wǎng)衣水阻力隨網(wǎng)線直徑的變化
表3 網(wǎng)線直徑變化引起的受力變化幅度統(tǒng)計表
由表3可知,在3種網(wǎng)目尺寸條件下,各工況下3組對比網(wǎng)衣的受力增加平均幅度分別為25.05%、16.0%、45.06%。在水流作用下,平面網(wǎng)衣受力大小與網(wǎng)衣在水流作用方向的投影面積成正比例關(guān)系,因此當(dāng)菱形鏈網(wǎng)平面尺寸及網(wǎng)目尺寸相同情況下,網(wǎng)線直徑增加,網(wǎng)衣在水流作用方向的投影面積增加,因此網(wǎng)衣整體受力也隨之增加。
本文基于有限元基本原理,采用梁單元和連接單元建立了金屬菱形鏈網(wǎng)模型,通過動力隱式求解節(jié)點運動方程。在數(shù)值模擬結(jié)果與試驗結(jié)果擬合良好的基礎(chǔ)上,分析討論了在水流作用下金屬菱形鏈網(wǎng)受力與網(wǎng)目尺寸及網(wǎng)線直徑之間的關(guān)系:在水流作用下,網(wǎng)衣水阻力隨著網(wǎng)目尺寸的增加而逐漸減小,當(dāng)網(wǎng)目尺寸由25 mm增加到35 mm、45 mm時,對應(yīng)的網(wǎng)衣受力的平均減小幅度分別為14.71%和38.07%;網(wǎng)衣水阻力隨著網(wǎng)線直徑的增加而逐漸增大,當(dāng)網(wǎng)線直徑由2.5 mm增加到3.2 mm、4.0 mm時,對應(yīng)的網(wǎng)衣受力的平均增加幅度分別為25.05%和45.06%。
□
[1] BESSONNEAU J S, MARICHAL D. Study of the dynamics of submerged supple nets (applications to trawls)[J]. Ocean Engineering, 1998, 25(7):563-583.
[2] TSUKROV I, EROSHKIN O, FREDRIKSSON D, et al. Finite element modeling of net panels using a consistent net element[J]. Ocean Engineering, 2003, 30(2):251-270.
[3] DECEW J, OSIENSKI M, DRACH A, et al. Effect of the bending stiffness on the volumetric stability of fish cages with copper alloy net[C]// V International Conference on Computational Methods in Marine Engineering. 2013.
[4] LI L, FU S, XU Y, et al. Dynamic responses of floating fish cage in waves and current[J]. Ocean Engineering, 2013, 72(7):297-303.
[5] LI L, FU S, LI R. Dynamic Responses of the Floating Cage System in Current and Waves[C]// ASME 2012, International Conference on Ocean, Offshore and Arctic Engineering. 2012:239-248.
[6] LI L, FU S, XU Y. Nonlinear hydroelastic analysis of an aquaculture fish cage in irregular waves[J]. Marine Structures, 2013, 34(34):56-73.
[7] 崔勇,關(guān)長濤,萬榮,等.基于有限元方法對波流場中養(yǎng)殖網(wǎng)箱的系統(tǒng)動力分析[J].工程力學(xué), 2010,27(5):250-256.
[8] 崔江浩.重力式養(yǎng)殖網(wǎng)箱耐流特性的數(shù)值模擬及仿真[D].青島:中國海洋大學(xué), 2005.
[9] 詹杰民,胡由展,趙陶,等.漁網(wǎng)水動力試驗研究及分析[J].海洋工程, 2002, 20(2):49-53.
[10] 李玉成,桂福坤.平面無結(jié)節(jié)網(wǎng)衣水阻力系數(shù)的試驗研究[J].海洋學(xué)報, 2006, 28(5):145-151.
[11] 李玉成,桂福坤.平面有結(jié)節(jié)和無結(jié)節(jié)網(wǎng)目試驗及水阻力系數(shù)的選擇[J].中國海洋平臺, 2005, 20(6):11-17.
[12] 黃小華,郭根喜,胡昱,等.波流作用下深水網(wǎng)箱受力及運動變形的數(shù)值模擬[J].中國水產(chǎn)科學(xué), 2011,18(2):443-450.
[13] 陳昌平,李玉成,趙云鵬,等.波流共同作用下單體網(wǎng)格式錨碇網(wǎng)箱水動力特性研究[J].水動力學(xué)研究與進展,2009, 24(4):493-502.
[14] CHEN C P, LI Y C, ZHAO Y P, et al. Numerical Analysis on the Effects of Submerged Depth of the Grid and Direction of Incident Wave on Gravity Cage[J]. China Ocean Engineering, 2009, 23(2):233-250.
[15] CHEN C P, ZHAO Y P, LI Y C, et al. Numerical analysis of hydrodynamic behaviors of two net cages with grid mooring system under wave action[J]. China Ocean Engineering, 2012, 26(1):59-76.
[16] ZHAO Y P, WANG X X, DECEW J, et al. Comparative Study of Two Approaches to Model the Offshore Fish Cages[J]. China Ocean Engineering, 2015, 29(3):459-472.
[17] BI C W, ZHAO Y P, DONG G H, et al. A numerical analysis on the hydrodynamic characteristics of net cages using coupled fluid-structure interaction model[J]. Aquacultural Engineering, 2014, 59(2):1-12.
[18] DRACH A, TSUKROV I, DECEW J, et al. Engineering procedures for design and analysis of submersible fish cages with copper net for exposed marine environment[J]. Aquacultural Engineering, 2016, 70:1-14.
[19] DRACH A, TSUKROV I, DECEW J, et al. Design and modeling of submersible fish cages with copper net for open ocean aquaculture[C]// International Conference on Computational Methods in Marine Engineering, 2014.
[20] 陳昌平, 王文, 鄭佳成,等. 水流作用下平面金屬網(wǎng)衣水阻力特性數(shù)值模擬[J]. 大連海洋大學(xué)學(xué)報, 2017, 32(3):363-368.
[21] 蘇煒,詹杰民.水流作用下漁網(wǎng)養(yǎng)殖空間變化的計算方法[J].海洋工程, 2007, 25(1):93-100.
[22] GIGNOUX H, THOMPSON L D, MESSIER R H. Computational Model of Aquaculture Fin-fish Net-pens[C]// Proceedings of an International Conference. University of Hawaii Sea Grant College Program. Maui, Hawaii, 1997:107-129.
[23] CHOO Y I, CASARELLA M J. Hydrodynamic Resistance of Towed Cables[J]. Journal of Hydronautics, 1971, 5(4):126-131.
[24] 王勖成.有限單元法[M].北京:清華大學(xué)出版社, 2003:469-522.
Studyofthenumericalsimulationofcharacteristicsofdifferentflowresistancesonametalchain-linknetforaquacultureundercurrents
LIUHangfei,CHENChangping,ZHENGYanna,LIANGXinyu
(SchoolofOceanandCivilEngineering,DalianOceanUniversity,Dalian,Liaoning116023,China)
It is an important work to study the hydrodynamic characteristics of the net structure to guarantee the safety of the whole net cage system in deep open sea. The most widely used copper alloy net in aquaculture cages is of the “chain-link” type mental net, which is a kind of porous structure with a small diameter, whose hydrodynamic characteristics under the waves and currents is significantly different from general ocean engineering structures. Based on the finite element method, the chain-link nets and the interaction between them are simulated by the beam element and connector element in the ABAQUS software. After verification by the experiment data, the different flow resistances on the net under currents with different net diameters and mesh sizes were analyzed. The numerical simulation results showed that as the mesh size increases from 25mm, 35mm to 45mm, the average value of the flow resistance increases with the amplification of 14.71% and 38.07% respectively under various operating modes. Similarly, the average value of flow resistance increases with the amplification of 25.05% and 45.06% respectively under various operating modes as the net diameter increases from 2.5 mm, 3.2 mm to 4.0 mm. The research provides a positive foundation for further study of the hydrodynamic characteristics of deep-sea aquaculture metal cages.
metal net; flow resistance characteristics; beam element; connector element; numerical simulation
10.3969/j.issn.1007-9580.2017.09.013
2017-11-01
國家自然科學(xué)基金資助項目(31572663)
劉航飛(1995—),男,碩士研究生,研究方向:水工結(jié)構(gòu)工程。E-mail:18340831879@163.com
鄭艷娜(1978—),女,副教授,博士,研究方向:港口、海岸及近海工程。E-mail:zhengyn@dlou.edu.cn
TV131.6
A
1007-9580(2017)06-073-07