固定方式對樁柱式圍網(wǎng)網(wǎng)片波浪力學(xué)特性影響研究

桂 福 坤， 陳 天 華， 趙 云 鵬， 董 國 海

( 1.浙江海洋大學(xué) 國家海洋設(shè)施養(yǎng)殖工程技術(shù)研究中心， 浙江 舟山 316022；2.大連理工大學(xué) 海岸和近海工程國家重點實驗室, 遼寧 大連 116024 )

固定方式對樁柱式圍網(wǎng)網(wǎng)片波浪力學(xué)特性影響研究

桂 福 坤*1， 陳 天 華1， 趙 云 鵬2， 董 國 海2

( 1.浙江海洋大學(xué) 國家海洋設(shè)施養(yǎng)殖工程技術(shù)研究中心， 浙江 舟山 316022；2.大連理工大學(xué) 海岸和近海工程國家重點實驗室, 遼寧 大連 116024 )

淺海圍網(wǎng)養(yǎng)殖是新興的一種生態(tài)型養(yǎng)殖模式，其中比較典型的樁柱式圍網(wǎng)，主體由排樁和網(wǎng)片組成．圍網(wǎng)網(wǎng)衣系統(tǒng)安全是決定圍網(wǎng)養(yǎng)殖工程安全的關(guān)鍵所在，因此針對樁柱式圍網(wǎng)中的單元網(wǎng)片，采用數(shù)值方法對其在波浪作用下的力學(xué)特性進行了研究，著重分析了不同網(wǎng)片尺度和固定方式下單元網(wǎng)片的網(wǎng)線張力分布、結(jié)節(jié)偏移和樁柱系縛點受力特性．結(jié)果表明，波浪條件下樁柱式圍網(wǎng)單元網(wǎng)片的網(wǎng)線最大張力主要出現(xiàn)在網(wǎng)片上端兩側(cè)位置，網(wǎng)線最大張力與網(wǎng)片尺寸成正比；單元網(wǎng)片結(jié)節(jié)最大偏移主要出現(xiàn)在網(wǎng)片的中上位置，同樣結(jié)節(jié)最大偏移與網(wǎng)片尺寸成正比．兩者均隨系縛點數(shù)量的增加先減小后在系縛點多于4個時趨于穩(wěn)定；網(wǎng)片與樁柱的系縛點最大受力出現(xiàn)在上端點，其次是下端點，且數(shù)值上均遠(yuǎn)遠(yuǎn)大于樁柱中間的系縛點，并隨網(wǎng)片尺寸的增大而增大，隨系縛點數(shù)量的增加而不斷減?。?/p>

樁柱式圍網(wǎng)養(yǎng)殖；單元網(wǎng)片；數(shù)值模擬；固定方式；波浪

0 引 言

淺海圍網(wǎng)養(yǎng)殖是繼陸基工廠化養(yǎng)殖、岸邊灘涂養(yǎng)殖和圍塘養(yǎng)殖、淺海網(wǎng)箱養(yǎng)殖之后研發(fā)的一種新興的海水生態(tài)養(yǎng)殖模式．樁柱式圍網(wǎng)主要由管樁和網(wǎng)片構(gòu)成．布置方式是圍網(wǎng)養(yǎng)殖工程最為重要的影響因素之一，圍網(wǎng)安全與固定方式緊密相關(guān)，因此，開展圍網(wǎng)網(wǎng)片在不同固定方式下受波浪作用的水力特性研究意義重大．近年來，國內(nèi)外已有許多有關(guān)網(wǎng)衣在波浪作用下的水動力特性研究，宋偉華等[1]采用4種網(wǎng)箱系泊方法，進行了單點系泊網(wǎng)衣構(gòu)件水槽波浪試驗研究，結(jié)果表明單點系泊在波浪條件下構(gòu)件水動力有較大的增加．李玉成等[2]利用模型試驗對單片網(wǎng)衣在波浪作用下的水動力特性進行了研究，其結(jié)果表明網(wǎng)衣密實度對網(wǎng)衣速度力系數(shù)的影響相對較大，網(wǎng)衣所受慣性力項和速度力相比很小，可以忽略．劉莉莉等[3]應(yīng)用有限元方法與集中參數(shù)法建立了張網(wǎng)在波流場中的數(shù)值模型，通過數(shù)值模擬來研究張網(wǎng)在波流作用下的水動力特性，并將數(shù)值計算結(jié)果與水槽模型試驗結(jié)果進行比較，符合良好．Chan等[4]通過數(shù)值模擬對穿過柔性漁網(wǎng)的表面波散射進行了研究分析，清楚地揭示了柔性滲透屏障對波浪傳播的影響規(guī)律．Tsukrov等[5]應(yīng)用有限元法建立了網(wǎng)片在波浪和水流環(huán)境負(fù)荷下的水動力響應(yīng)模型，并將其結(jié)果應(yīng)用到計算張力腿網(wǎng)箱中．葉衛(wèi)富[6]通過模型試驗對浮繩式圍網(wǎng)的水動力特性進行了初步研究，但受試驗條件限制，測試內(nèi)容及結(jié)果的應(yīng)用有一定局限性．

樁柱式圍網(wǎng)養(yǎng)殖是一種新興的海洋養(yǎng)殖模式，結(jié)構(gòu)上有其特殊性，對圍網(wǎng)網(wǎng)衣系統(tǒng)的水動力研究目前少有報道．淺海圍網(wǎng)工程的網(wǎng)衣系統(tǒng)與深水網(wǎng)箱、捕撈漁具結(jié)構(gòu)上有較多不同，網(wǎng)衣系統(tǒng)安全是決定圍網(wǎng)養(yǎng)殖工程安全的關(guān)鍵所在．本文采用集中質(zhì)量點法對樁柱間單元網(wǎng)片建立波浪作用下網(wǎng)片的運動響應(yīng)數(shù)學(xué)模型，并通過計算機數(shù)值模擬研究不同布置方式下圍網(wǎng)網(wǎng)片的網(wǎng)線張力分布、結(jié)節(jié)偏移和系縛點受力，為樁柱式圍網(wǎng)設(shè)計、制作、海上敷設(shè)和圍網(wǎng)抗風(fēng)浪技術(shù)的研發(fā)等提供理論依據(jù)和技術(shù)支撐．

1 數(shù)值模擬方法

1.1 數(shù)學(xué)模型

圍網(wǎng)網(wǎng)片屬于典型的柔性結(jié)構(gòu)物，可采用集中質(zhì)量法進行模擬．假定網(wǎng)片是由有限的無質(zhì)量彈簧連接的集中質(zhì)量點所構(gòu)成，通過計算集中質(zhì)量點在波浪和邊界條件作用下的位移，來得到網(wǎng)的形狀．模型的集中質(zhì)量點設(shè)于每個網(wǎng)目目腳的兩端，每個集中質(zhì)量點包含網(wǎng)目的一個結(jié)節(jié)和兩個目腳，如圖1(a)所示．把結(jié)節(jié)近似看成圓球，目腳看作圓柱形桿件，所以它具有圓柱的水動力性質(zhì)．

(a) 集中質(zhì)量點示意圖

(b) 單元網(wǎng)片示意圖

圖1 單元網(wǎng)片計算模型示意圖

Fig.1 The schematic diagram of calculation model for unit net panel

網(wǎng)衣集中質(zhì)量點受力主要包括重力、浮力、網(wǎng)線張力、波浪力等，網(wǎng)衣受到的波浪力可根據(jù)莫里森方程[7]來計算．在計算目腳的波浪力時需考慮波浪入射方向與網(wǎng)線夾角關(guān)系，本文在目腳上建立局部坐標(biāo)系(ξ，η，τ)，τ方向為沿目腳方向，ξ軸為在τ和水質(zhì)點相對速度vr組成的平面內(nèi)與τ垂直[8-10]．在整體坐標(biāo)系下，局部坐標(biāo)系的ξ、η、τ軸單位矢量可以通過已知的水質(zhì)點相對速度vr=vw-v0和目腳方向τ=(x1-x0y1-y0z1-z0)向量叉乘得到，為

FD bar=FDξ+FDη+FDτFI bar=FIξ+FIη+FIτ

其中

在整體坐標(biāo)系下將集中質(zhì)量點所包含的結(jié)節(jié)和目腳的受力進行累加：

T=T1+T2+T3+T4FD=FD node+FD barFI=FI node+FI barW=Wnode+WbarB=Bnode+Bbar

然后利用牛頓第二定律建立質(zhì)點運動方程．網(wǎng)衣模型簡化后的質(zhì)點運動方程[11-13]為

(M+ΔM)a=T+FD+FI+W+B

其中

式中：ΔM、M分別表示各集中質(zhì)量點的附加質(zhì)量和質(zhì)量，a表示各集中質(zhì)量點的加速度矢量，T表示集中質(zhì)量點所受到的張力矢量，F(xiàn)D、FI分別表示集中質(zhì)量點的速度力矢量和慣性力矢量，W表示質(zhì)點的重力矢量，B表示質(zhì)點的浮力矢量，ρw表示水體的密度，CD表示速度力項系數(shù)，A表示網(wǎng)線沿波浪方向的投影面積，v、vr分別表示波浪水質(zhì)點的速度矢量和相對速度矢量，CM表示慣性力項系數(shù)，?表示網(wǎng)線構(gòu)件的體積，Cm表示附加的質(zhì)量力項系數(shù)，d為網(wǎng)線直徑，l0為網(wǎng)線原始長度，l為變形后的長度，C1、C2為構(gòu)件材料彈性系數(shù)．

建立的質(zhì)點運動方程為典型的二階偏微分方程，可采用歐拉法、龍格庫塔法等多種方法求解．本文利用Fortran軟件編程計算求解，采用歐拉法即可獲得良好的收斂性．計算中首先根據(jù)網(wǎng)片的初始狀態(tài)，計算網(wǎng)線上的波浪力以及網(wǎng)線變形所產(chǎn)生的拉力，基于質(zhì)點運動方程求解出質(zhì)點的加速度，然后根據(jù)向前歐拉公式[14]：yi+1=yi+h·f(xi，yi)，求出下一時間步長的質(zhì)點位移和速度，從而確定網(wǎng)片形狀，最后以求出的質(zhì)點位置和速度作為網(wǎng)片新的狀態(tài)參數(shù)，重復(fù)上述步驟直至結(jié)束．

1.2 模型驗證

本文引用文獻[15]中的預(yù)加張力放射法系泊網(wǎng)衣波浪試驗情況，對實際漁用網(wǎng)衣J，在波浪周期1.6 s試驗條件下進行了數(shù)值模擬．網(wǎng)衣J尺寸為90.5 cm×80.6 cm(寬×高)，材料為PE，縮結(jié)系數(shù)為0.707，網(wǎng)目大小為4 cm，網(wǎng)線直徑為1.95 mm，網(wǎng)目數(shù)為32×28.5．網(wǎng)衣框架為HDPE空心管材，管材外徑14 mm，框架底部距離池底100 mm，水深0.7 m．圖2給出了各種波高下網(wǎng)衣水平波浪力最大值的模擬結(jié)果．由圖可見，在各種不同的波高條件下，網(wǎng)衣水平波浪力最大值的模擬值和試驗值吻合較好．

圖2 網(wǎng)衣水平波浪力模擬值和試驗值對比

1.3 計算參數(shù)選取

圍網(wǎng)所用網(wǎng)衣有多種材質(zhì)，大多采用PE網(wǎng)和超高相對分子質(zhì)量纖維網(wǎng)，有些圍網(wǎng)也采用金屬網(wǎng)(如銅網(wǎng))．網(wǎng)衣的尺寸也有所不同，單元網(wǎng)片寬度一般為3～5 m，高度一般為6～10 m．本文重點研究波浪對網(wǎng)片不同部位受力分布特性的影響，簡單起見，網(wǎng)片全部采用PE材質(zhì)，網(wǎng)目大小2a=8 cm，網(wǎng)線直徑d=3 mm，水平縮結(jié)系數(shù)0.66，垂直縮結(jié)系數(shù)0.75．采用線性波浪理論，水深10 m．為避免波浪條件下，網(wǎng)片露出水面而導(dǎo)致受力差異，研究時，網(wǎng)片上邊纜水下布置深度為二分之一波高．本文主要關(guān)注網(wǎng)片尺寸和系縛方式對圍網(wǎng)網(wǎng)片的影響，因此波浪要素僅給定一種情況，即取波高H=3 m，波向θ=90°，波周期T=7 s．網(wǎng)片尺寸取寬×高=1 m×6 m，2 m×6 m，3 m×6 m，4 m×6 m和5 m×6 m；系縛方式主要指單元網(wǎng)片在圍網(wǎng)樁柱上的系縛點數(shù)量，共設(shè)置6種情況，即單邊系縛點數(shù)量分別為2、3、4、5、7、13個，如圖3所示．網(wǎng)片四邊用10 mm繩索加勁，如圖1(b)所示．由于計算中網(wǎng)片結(jié)節(jié)數(shù)較多，為減少計算時間、提高計算效率，采用網(wǎng)目群化方法[16-17]將相鄰64個網(wǎng)目合并為一個等效大網(wǎng)目．

圖3 網(wǎng)片系縛方式示意圖

網(wǎng)目群化后，13點系縛為完全系縛，即計算模型中，每個邊側(cè)節(jié)點均被約束．

2 計算結(jié)果與分析

2.1 網(wǎng)片尺寸對圍網(wǎng)網(wǎng)片波浪力學(xué)特性的影響

研究時，給定一種網(wǎng)片系縛方式，即5點系縛，此時系縛點間距與網(wǎng)片橫寬比按網(wǎng)片尺寸的增加依次為1.500、0.750、0.500、0.375、0.300．重點關(guān)注5種不同尺寸的單元網(wǎng)片，在波浪作用下的網(wǎng)線最大張力和結(jié)節(jié)最大偏移，以及網(wǎng)片與樁柱系縛點的受力特性．結(jié)節(jié)最大偏移指的是網(wǎng)衣變形后，結(jié)節(jié)到網(wǎng)片初始平面的最大偏移距離．

2.1.1 網(wǎng)線張力與結(jié)節(jié)偏移 圖4和圖5給出了網(wǎng)線最大張力、結(jié)節(jié)最大偏移與網(wǎng)片尺寸的關(guān)系．由圖可見，兩者均隨網(wǎng)片尺寸增大而增大．網(wǎng)線最大張力出現(xiàn)在網(wǎng)片上端兩側(cè)位置(如圖6實粗線所示)，但尺寸為w0(寬)×l0(長)=1 m×6 m 的網(wǎng)片網(wǎng)線最大張力出現(xiàn)位置有所不同．網(wǎng)片中結(jié)節(jié)最大偏移位置出現(xiàn)在網(wǎng)片中上部(如圖6星號所示)．

圖4 網(wǎng)線最大張力隨網(wǎng)片尺寸的變化

圖5 結(jié)節(jié)最大偏移隨網(wǎng)片尺寸的變化

2.1.2 圍網(wǎng)網(wǎng)片與樁柱系縛點之間的受力特性 圖7給出了網(wǎng)片在圍網(wǎng)樁柱上的5個系縛點(系縛點編號No.為1、4、7、10、13，見圖3)的最大受力值與網(wǎng)片尺寸的關(guān)系．由圖可見，首尾兩端系縛點(No.1和13)的受力遠(yuǎn)大于中間系縛點的受力，且最大受力值隨網(wǎng)片尺寸的增大而迅速增大；中間系縛點(No.4、7、10)的最大受力隨網(wǎng)片尺寸的增大變化相對較?。?/p>

圖6 不同網(wǎng)片尺寸下網(wǎng)線張力最大的部位(實粗線)和結(jié)節(jié)偏移最大的部位(星號)

Fig.6 The position of the maximum twine force (heavy lines) and the maximum offset of the nodes (pentagram) under different net panel′s sizes

圖7 不同尺寸網(wǎng)片各系縛點的最大受力情況

2.2 系縛方式對圍網(wǎng)網(wǎng)片波浪受力特性的影響

在給定網(wǎng)片尺寸(寬×高=3m×6m)條件下，本文研究了系縛方式對圍網(wǎng)網(wǎng)片的網(wǎng)線最大張力和結(jié)節(jié)最大偏移以及網(wǎng)片與樁柱系縛點受力的影響．此時系縛點間距與網(wǎng)片橫寬比按系縛點數(shù)的增加依次為2.000、1.000、0.667、0.500、0.333、0.020．

2.2.1 最大網(wǎng)線張力與結(jié)節(jié)偏移 圖8和9給出了網(wǎng)線最大張力及結(jié)節(jié)最大偏移與系縛方式的關(guān)系．由圖可見，增加系縛點的數(shù)量N，有助于減小網(wǎng)線的最大張力和結(jié)節(jié)的最大偏移，但當(dāng)系縛點超過一定數(shù)量(本文為4點)后，網(wǎng)線最大張力和結(jié)節(jié)最大偏移不再出現(xiàn)顯著的降低．圖10中，實粗線給出了網(wǎng)線最大張力出現(xiàn)在網(wǎng)片上的位置．由圖可見，兩點系縛時網(wǎng)線最大張力出現(xiàn)在網(wǎng)片上下兩側(cè)位置，其他系縛方式網(wǎng)線最大張力出現(xiàn)在網(wǎng)片上端兩側(cè)位置．系縛點數(shù)量同樣對結(jié)節(jié)最大偏移位置有影響，隨著系縛點數(shù)量增加，最大偏移位置發(fā)生上移，并當(dāng)系縛點達(dá)到一定數(shù)量(本文為4點)后，位置不再有大的變化，如圖10所示．

圖8 網(wǎng)線最大張力隨系縛方式的變化

圖9 結(jié)節(jié)最大偏移隨系縛方式的變化

2.2.2 圍網(wǎng)網(wǎng)片與樁柱系縛點之間的最大受力 圖11給出了網(wǎng)片在不同系縛方式下系縛點的受力情況．由圖可見，網(wǎng)片上端系縛點最大受力均略大于下端系縛點最大受力，且均遠(yuǎn)遠(yuǎn)大于中間系縛點受力．隨著系縛點數(shù)的增加，系縛點最大受力整體上在逐漸減小．

圖10 不同系縛方式下網(wǎng)線張力最大的部位(粗線)和結(jié)節(jié)偏移最大的部位(星號)

Fig.10 The location of the maximum twine force (heavy lines) and the maximum node offset (pentagram) with different fixations

圖11 6種系縛方式下網(wǎng)片各系縛點的最大受力情況

Fig.11 The maximum force at each fixing point with 6 types of fixations

3 討 論

3.1 網(wǎng)片尺寸的影響

樁柱式圍網(wǎng)主體由網(wǎng)片和排樁組成，其中圍網(wǎng)網(wǎng)片起到防止魚類逃脫和敵類侵害等作用，由此可見，保證圍網(wǎng)網(wǎng)衣系統(tǒng)安全十分重要．對圍網(wǎng)網(wǎng)片采用不同的固定方式在波浪中會產(chǎn)生不一樣的效果，其在波浪作用下的網(wǎng)線最大張力、結(jié)節(jié)最大偏移和樁柱系縛點最大受力是圍網(wǎng)規(guī)劃設(shè)計者最為關(guān)心的指標(biāo)，也是保證惡劣海況下圍網(wǎng)養(yǎng)殖安全的首要條件．因此，圍網(wǎng)網(wǎng)片的網(wǎng)線最大張力、結(jié)節(jié)最大偏移和樁柱系縛點最大受力這3個指標(biāo)在樁柱式圍網(wǎng)養(yǎng)殖工程的風(fēng)險評估中都具有十分重要的作用．

本文研究發(fā)現(xiàn)，網(wǎng)線最大張力和結(jié)節(jié)最大偏移都與網(wǎng)片尺寸成正比．尺寸為1m×6m的網(wǎng)片因為寬度小于系縛點間距，同時波浪能量隨水深衰減[18]，導(dǎo)致最大張力發(fā)生在網(wǎng)片兩側(cè)上部與系縛點相連的位置，如圖6所示．而其他網(wǎng)片的寬度都大于系縛點間距，同時伴隨波浪的影響導(dǎo)致最大張力發(fā)生在網(wǎng)片上端兩側(cè)位置．但所有的網(wǎng)線最大張力都發(fā)生在加勁繩上，這也表明對網(wǎng)片四邊進行加勁處理是非常必要的．結(jié)節(jié)最大偏移位置位于網(wǎng)片中上部，跟波浪能量隨水深的變化有關(guān)，在網(wǎng)片此部位的最大偏移以內(nèi)盡量避免有障礙物，否則易發(fā)生網(wǎng)線因磨損而斷裂．研究結(jié)果表明，圍網(wǎng)網(wǎng)片上系縛點的整體受力呈現(xiàn)首尾兩端大中間小的現(xiàn)象，而且隨著網(wǎng)片尺寸的增加，首尾系縛點的受力會不斷增大，表明改變網(wǎng)片尺寸對首尾系縛點產(chǎn)生的影響大，而對中間系縛點的影響不明顯．因此，對網(wǎng)片首尾兩端的系縛點進行加固是非常必要的．

3.2 系縛方式的影響

研究發(fā)現(xiàn)，網(wǎng)片兩側(cè)系縛點受力整體呈現(xiàn)出首尾兩端大中間小的現(xiàn)象，其中首端的系縛點受力最大．因此，對網(wǎng)片首尾兩端的系縛點進行加固是非常必要的．同時，網(wǎng)線最大張力和結(jié)節(jié)最大偏移與系縛點數(shù)量有關(guān)，隨系縛點數(shù)量的增加先減小后趨于穩(wěn)定，表明增加一定數(shù)量的系縛點可以有效降低整個網(wǎng)片中的網(wǎng)線最大張力和結(jié)節(jié)最大偏移，從而降低網(wǎng)線斷裂的風(fēng)險．同時，隨著系縛點數(shù)量的增加，網(wǎng)片首尾系縛點的最大受力也不斷下降，表明系縛點數(shù)量的增加同樣可以分擔(dān)首尾兩端系縛點的受力．從整體受力看，系縛點達(dá)到一定數(shù)量后，網(wǎng)線最大張力趨于穩(wěn)定，此時進一步增加系縛點數(shù)量反而增加了圍網(wǎng)工程的施工難度．因此，圍網(wǎng)網(wǎng)片應(yīng)存在一個臨界系縛點數(shù)量，本文研究該臨界系縛點數(shù)量為4～5個，對于不同圍網(wǎng)網(wǎng)片尺度，需要具體分析．研究同時發(fā)現(xiàn)，當(dāng)采用2點系縛時，由于豎向系縛點之間的間距大于網(wǎng)片的橫向?qū)挾?，整個網(wǎng)片的最大受力位置出現(xiàn)在了左右側(cè)纜上．隨著兩側(cè)系縛點數(shù)量的增加，當(dāng)系縛點的間距小于網(wǎng)片的橫向?qū)挾龋霈F(xiàn)網(wǎng)線張力最大的部位轉(zhuǎn)移至上邊纜兩端，再次證明了對網(wǎng)片四邊做加勁處理的重要性．

4 結(jié) 論

(1)網(wǎng)片尺寸和系縛方式均對網(wǎng)片結(jié)構(gòu)單元的受力和變形特性有著顯著的影響．單元網(wǎng)片中的網(wǎng)線最大張力和結(jié)節(jié)最大偏移均隨著網(wǎng)片尺寸的增大而增大，網(wǎng)線最大張力一般出現(xiàn)于網(wǎng)片頂端兩側(cè)系縛點處，結(jié)節(jié)最大偏移出現(xiàn)于單元網(wǎng)片的中上部．實際應(yīng)用時，建議強化頂部綱繩的設(shè)計考慮．

(2)波浪條件下，網(wǎng)片系縛點最大受力出現(xiàn)于頂部第一個系縛點處，其次為底端系縛點，中間系縛點受力均遠(yuǎn)小于上下兩端系縛點受力．實際應(yīng)用時，建議強化頂部和底部系縛點的設(shè)計考慮．

(3)增加單元網(wǎng)片兩側(cè)系縛點的數(shù)量，有利于降低所有系縛點的受力．但對單元網(wǎng)片中的網(wǎng)線受力而言，當(dāng)系縛點數(shù)量增至4個以上時，網(wǎng)線的最大受力不再顯著下降．因此從這一角度分析，應(yīng)存在一臨界系縛點數(shù)量，本文研究條件下，臨界系縛點數(shù)量為4～5個．實際應(yīng)用應(yīng)保證達(dá)到臨界數(shù)量以上，并保留一定的富余，具體數(shù)量需視具體工程而定．

本文研究結(jié)果可為樁柱式圍網(wǎng)養(yǎng)殖工程設(shè)計和施工提供重要的參考依據(jù)．本文僅針對一種波浪條件進行了研究，有關(guān)波浪要素等其他因素對單元網(wǎng)片波浪力學(xué)特性的影響將于后期進一步研究．

[1] 宋偉華,梁振林,趙芬芳,等. 單點系泊網(wǎng)衣構(gòu)件波浪試驗研究[J]. 海洋與湖沼, 2005, 36(3):199-206.

SONGWeihua,LIANGZhenlin,ZHAOFenfang, et al.Waveexperimentonnettingstructurewithsinglepointmoorings[J].OceanologiaetLimnologiaSinica, 2005, 36(3):199-206. (inChinese)

[2] 李玉成,趙云鵬,桂福坤,等. 波浪條件下漁網(wǎng)水動力特性研究[J]. 中國造船, 2006, 47(s1):265-273.

LIYucheng,ZHAOYunpeng,GUIFukun, et al.Hydrodynamicstudyonasubmergedplanenetsinwave[J].ShipbuildingofChina, 2006, 47(s1):265-273. (inChinese)

[3] 劉莉莉,萬 榮,黃六一,等. 波流場中張網(wǎng)漁具水動力學(xué)特性的數(shù)值模擬[J]. 中國海洋大學(xué)學(xué)報(自然科學(xué)版), 2013, 43(5):24-29.

LIULili,WANRong,HUANGLiuyi, et al.Numericalsimulationofhydrodynamiccharacteristicsofstownetinwaveandcurrent[J].PeriodicalofOceanUniversityofChina, 2013, 43(5):24-29. (inChinese)

[4]CHANAT,LEESWC.Wavecharacteristicspastaflexiblefishnet[J].OceanEngineering, 2001, 28(11):1517-1529.

[5]TSUKROVI,EROSHKINO,FREDRIKSSOND, et al.Finiteelementmodelingofnetpanelsusingaconsistentnetelement[J].OceanEngineering, 2003, 30(2):251-270.

[6] 葉衛(wèi)富. 淺海養(yǎng)殖圍網(wǎng)設(shè)施水動力性能的初步研究[D]. 杭州:浙江海洋學(xué)院, 2012.

YEWeifu.Preliminarystudyonthecharacteristicsofhydrodynamicsontheaquacultureseine[D].Hangzhou:ZhejiangOceanUniversity, 2012. (inChinese)

[7] 王樹青,梁丙臣. 海洋工程波浪力學(xué)[M]. 青島:中國海洋大學(xué)出版社, 2013.

WANGShuqing,LIANGBingchen.WaveMechanicsforOceanEngineering[M].Qingdao:ChinaOceanUniversityPress, 2013. (inChinese)

[8] 黃小華,郭根喜,胡 昱,等. 水流作用網(wǎng)衣過程的數(shù)值模擬[J]. 南方水產(chǎn)科學(xué), 2011, 7(3):56-61.

HUANGXiaohua,GUOGenxi,HUYu, et al.Numericalsimulationofdynamicprocessofcylindernetsincurrent[J].SouthChinaFisheriesScience, 2011, 7(3):56-61. (inChinese)

[9] 黃小華,郭根喜,陶啟友,等. 平面網(wǎng)衣在水流作用下的受力和變形特性數(shù)值模擬研究[J]. 南方水產(chǎn), 2009, 5(3):23-29.

HUANGXiaohua,GUOGenxi,TAOQiyou, et al.Numericalsimulationoftheforceanddeformationofsubmergedplanenetsincurrent[J].SouthChinaFisheriesScience, 2009, 5(3):23-29. (inChinese)

[10] 黃小華,郭根喜,胡 昱,等. 圓形網(wǎng)衣在水流作用下的運動變形特性[J]. 中國水產(chǎn)科學(xué), 2010, 17(2):312-319.

HUANGXiaohua,GUOGenxi,HUYu, et al.Movementanddeformationcharacteristicsofcylindernetsincurrent[J].JournalofFisherySciencesofChina, 2010, 17(2):312-319. (inChinese)

[11] 趙云鵬. 深水重力式網(wǎng)箱水動力特性數(shù)值模擬研究[D]. 大連:大連理工大學(xué), 2007.

ZHAOYunpeng.Numericalinvestigationonhydrodynamicbehaviorofdeep-watergravitycage[D].Dalian:DalianUniversityofTechnology, 2007. (inChinese)

[12] 趙云鵬,李玉成,董國海,等. 水流作用下重力式網(wǎng)箱網(wǎng)衣張力分布[J]. 漁業(yè)現(xiàn)代化, 2008, 35(6):5-8, 63.

ZHAOYunpeng,LIYucheng,DONGGuohai, et al.Tensiondistributionofgravitycagenetinsteadycurrent[J].FisheryModernization, 2008, 35(6):5-8, 63. (inChinese)

[13] 陳小芳. 鲆鰈類方形網(wǎng)箱水動力特性數(shù)值模擬[D]. 大連:大連理工大學(xué), 2012.

CHENXiaofang.Numericalsimulationofthehydrodynamicbehaviorofthesquarenetcageforflatfish[D].Dalian:DalianUniversityofTechnology, 2012. (inChinese)

[14] 劉師少. 計算方法[M]. 北京:科學(xué)出版社, 2005:151-172.

LIUShishao.ComputingMethod[M].Beijing:SciencePress, 2005:151-172. (inChinese)

[15] 宋偉華. 網(wǎng)衣波浪水動力學(xué)研究[D]. 青島:中國海洋大學(xué), 2006.

SONGWeihua.Studyonthemechanicsofwaveforcesonthenettingpanel[D].Qingdao:OceanUniversityofChina, 2006. (inChinese)

[16]BESSONNEAUJS,MARICHALD.Studyofthedynamicsofsubmergedsupplenets(applicationstotrawls) [J].OceanEngineering, 1998, 25(7):563-583.

[17] 蘇 煒,詹杰民. 等效網(wǎng)面法在模擬網(wǎng)的水動力特性中的應(yīng)用[J]. 水動力學(xué)研究與進展A輯, 2007, 22(3):267-272.

SUWei,ZHANJiemin.Applicationofanequivalentmethodtomodeldynamicpropertiesoffishnet[J].JournalofHydrodynamics,Ser.A, 2007, 22(3):267-272. (inChinese)

[18] 竺艷蓉. 海洋工程波浪力學(xué)[M]. 天津:天津大學(xué)出版社, 1991.

ZHUYanrong.WaveMechanicsforOceanEngineering[M].Tianjin:TianjinUniversityPress, 1991. (inChinese)

Study of effect on wave mechanical properties for net panel of pile-column type net enclosure by fixations

GUI Fukun*1, CHEN Tianhua1, ZHAO Yunpeng2, DONG Guohai2

( 1.National Engineering Research Center for Marine Aquaculture, Zhejiang Ocean University, Zhoushan 316022, China；2.State Key Laboratory of Coastal and Offshore Engineering, Dalian University of Technology, Dalian 116024, China )

The net enclosure aquaculture (NEA) in shallow water is a new developing ecological aquaculture mode. The pile-column type net enclosure is one of the most typical structures, which is mainly composed of piles and net panels. The net system is the most important part which affects the security of the whole NEA system, therefore the mechanical performance of unit net panels under different wave conditions is investigated using numerical method. The effects of net panel size and fixations on the tension force in twines, offset of nodes and loads at fixing points are studied. The experimental results show that the maximum twine force occurs mainly at the upper two sides of the unit net panel and increases as the size of net panel increases. The maximum node offset appears at the mid-upper part of net panel and also increases with the increase of the size of net panel. Both of the maximum twine force and the maximum node offset decrease as the fixing points number increases. However, they show little change when the number of the fixing points is over four. Regarding of the force at the fixing points, the largest force appears at the upper, then the bottom fixing points, and they are much higher than those at the other fixing points. The force increases with the increase of the size of net panel, and decreases with the increase of the fixing points number.

pile-column type net enclosure aquaculture; unit net panel; numerical simulation; fixations; wave

1000-8608(2017)03-0285-08

2016-07-05；

2017-03-28．

國家自然科學(xué)基金資助項目(51239002)；國家海洋局公益性行業(yè)科研專項經(jīng)費資助項目(201505025-2)；舟山市海洋專項(2015C41001).

桂福坤*(1976-)，男，教授，E-mail:gui2237@163.com.

S955.9

A

10.7511/dllgxb201703011