劉圣聰+劉佳茗

摘 要：拖網(wǎng)網(wǎng)板為網(wǎng)衣提供水平張力，是拖網(wǎng)作業(yè)系統(tǒng)非常重要的漁具構(gòu)件之一。該文利用CFD軟件FLUENT研究了展弦比對(duì)矩形曲面網(wǎng)板其水動(dòng)力性能的影響。研究中設(shè)計(jì)了2種不同展弦比λ的網(wǎng)板，在流速1.54m/s，迎流沖角α = 0°～50°時(shí)，建立數(shù)值水槽進(jìn)行數(shù)值模擬，得到網(wǎng)板的阻力系數(shù)Cd、升力系數(shù)Cl、俯仰力矩系數(shù)Cm和升阻比K，對(duì)比不同網(wǎng)板的水動(dòng)力性能差異。結(jié)果表明，1號(hào)網(wǎng)板（λ = 2.0）與2號(hào)網(wǎng)板（λ = 1.5）的最大升力系數(shù)分別為0.85（α =15°）和0.92（α = 42°）；最大升阻比分別為9.59和8.35，俯仰力矩系數(shù)的絕對(duì)值分別為0.013和0.001?？梢?jiàn)，展弦比越大，臨界沖角和Cl值越小，最大升阻比越高，但穩(wěn)定性越差。研究結(jié)果可為拖網(wǎng)網(wǎng)板的結(jié)構(gòu)優(yōu)化設(shè)計(jì)提供參考。

關(guān)鍵詞：拖網(wǎng)網(wǎng)板；展弦比；水動(dòng)力學(xué)性能；數(shù)值模擬

中圖分類號(hào) S141 文獻(xiàn)標(biāo)識(shí)碼 A 文章編號(hào) 1007-7731（2017）18-0092-04

Numerical Simulation on Hydrodynamic Performance of Rectangular Otter Boards with Different Aspect Ratios

Liu Shengcong et al.

（Dalian Tianzheng Industrial Co. Ltd，Dalian 116001，China）

Abstract：The otter board is a vital component for the efficient expansion performance of an otter trawl system. The numerical CFD code ANSYS fluent was conducted for rectangular otter board model to study the effect of aspect ratio on hydrodynamic performance. 2 different kinds of aspect ratio of otter boards were numerical simulated using the FLUENT with the flow velocity at 1.54m/s and the angle of attack measured from 0° to 50°. The cal-culated drag coefficient Cd，lift coefficient Cl，pitching moment coefficient Cm，lift to drag ratio K and the relations curve of these value and angle of attack were compared. The results showed that the maximum lift coefficient of otter boards 1 and 2 was 0.85（α =15°） and 0.92（α =42°），the maximum lift to drag ratio was 9.59 and 8.35，and the absolute value of pitching moment coefficient was 0.013 and 0.001 respectively. The result indicated that the greater the aspect ratio，the smaller the critical angle of attack and lift coefficient，the higher the maximum lift to drag ratio，but the worse of stability. The results would provide reference for the structural optimization design of trawl otter board.

Key words：Otter boards；Aspect ratio；Hydrodynamic performance；CFD analysis

拖網(wǎng)作業(yè)具有靈活、高效、適應(yīng)性強(qiáng)等特點(diǎn)，是最常見(jiàn)的一種捕撈方式，也是目前我國(guó)遠(yuǎn)洋捕撈中最主要的作業(yè)方式之一。網(wǎng)板是單船拖網(wǎng)系統(tǒng)中的重要組成部分，其水動(dòng)力性能的優(yōu)劣，關(guān)系到網(wǎng)具能否得到充分的擴(kuò)張，從而直接影響到拖網(wǎng)漁具的漁獲性能和經(jīng)濟(jì)效益[1-2]。因此，對(duì)網(wǎng)板水動(dòng)力性能的研究以及對(duì)新型網(wǎng)板的開(kāi)發(fā)，一直是國(guó)內(nèi)外漁業(yè)工程研究者的關(guān)注重點(diǎn)。

挪威、日本、法國(guó)等漁業(yè)發(fā)達(dá)國(guó)家對(duì)網(wǎng)板水動(dòng)力性能的改進(jìn)研究較早，相繼開(kāi)發(fā)出大展弦比矩形曲面網(wǎng)板、綜合型網(wǎng)板（橢圓形曲面開(kāi)縫網(wǎng)板）、立式V型曲面網(wǎng)板和矩形V型曲面開(kāi)縫網(wǎng)板等各種高性能網(wǎng)板[3]。A.Sala等[4]通過(guò)模型與原型試驗(yàn)相對(duì)比的方法，研究了2種矩形V型曲面開(kāi)縫網(wǎng)板的水動(dòng)力性能及其對(duì)海底棲息地的破壞影響。Yuki等[5]用計(jì)算流體動(dòng)力學(xué)方法研究了復(fù)翼式曲面網(wǎng)板的水動(dòng)力特性，并通過(guò)軟件可視化將網(wǎng)板周圍的流場(chǎng)形態(tài)展現(xiàn)出來(lái)，與模型試驗(yàn)進(jìn)行對(duì)比，得到較好的吻合效果。Balash等[6]研究了2種柔性蝠翼式網(wǎng)板，分析了網(wǎng)板的膨脹和扭曲對(duì)網(wǎng)板水動(dòng)力特性及穩(wěn)定性的影響，發(fā)現(xiàn)膨脹可以提高升力，扭曲可以減小阻力從而提高升阻比，適當(dāng)?shù)呐で鷮?duì)網(wǎng)板穩(wěn)定性是有利的。

中國(guó)對(duì)拖網(wǎng)網(wǎng)板水動(dòng)力性能的研究起步較晚，網(wǎng)板的結(jié)構(gòu)形式較國(guó)外比較單一。中國(guó)的網(wǎng)板模型試驗(yàn)主要在風(fēng)洞設(shè)備中完成。郭根喜等[7]較系統(tǒng)的對(duì)網(wǎng)板進(jìn)行了研究，包括網(wǎng)板的幾何結(jié)構(gòu)、開(kāi)縫數(shù)量、導(dǎo)流翼相對(duì)高度對(duì)不同類型網(wǎng)板水動(dòng)力性能的影響。張勛等[8]、王錦浩等[9]分別研究了矩形V型曲面網(wǎng)板開(kāi)縫口導(dǎo)流板曲率、展弦比、縫口位置和縫口寬度等參數(shù)對(duì)網(wǎng)板水動(dòng)力性能的影響。王磊等[10]研究了葉板尺度比例變化對(duì)單縫立式曲面網(wǎng)板水動(dòng)力性能的影響，發(fā)現(xiàn)導(dǎo)流板尺度過(guò)小會(huì)增加網(wǎng)板阻力，過(guò)長(zhǎng)則會(huì)降低網(wǎng)板升阻比。劉健等[11]通過(guò)水槽模型試驗(yàn)，研究了展弦比為1.0的矩形曲面網(wǎng)板在不同水平傾角和迎流沖角下的水動(dòng)力性能。endprint

國(guó)內(nèi)目前對(duì)網(wǎng)板水動(dòng)力性能的研究手段多數(shù)是在風(fēng)洞和水槽中進(jìn)行模型試驗(yàn)，對(duì)設(shè)備要求較高，同時(shí)需要大量的人力財(cái)力，且誤差較大。隨著計(jì)算機(jī)模擬仿真技術(shù)的發(fā)展，利用數(shù)值模擬的方法對(duì)網(wǎng)板水動(dòng)力性能進(jìn)行研究，將會(huì)更加省時(shí)省力，研究?jī)?nèi)容也會(huì)更加靈活。本文利用CFD軟件FLUENT對(duì)矩形曲面網(wǎng)板進(jìn)行數(shù)值模擬，研究展弦比變化對(duì)網(wǎng)板水動(dòng)力性能的影響，以期為優(yōu)化拖網(wǎng)網(wǎng)板性能提供參考。

1 數(shù)值模型

1.1 控制方程 本文假設(shè)流體不可壓縮、各向同性，控制方程包括連續(xù)性方程和N-S方程：

[??V=0] （1）

[?V?t+（V??）?V=-1ρ?p+μρ?2V+F] （2）

式中：[V]為速度；[t]為時(shí)間；[p]為壓力；[F]為體積力；[ρ]為密度；[μ]為動(dòng)力粘度。

1.2 參數(shù)定義 描述網(wǎng)板水動(dòng)力特性的參數(shù)主要有雷諾數(shù)Re、阻力系數(shù)Cd、升力系數(shù)Cl、升阻比K、縱向力矩系數(shù)Cm。計(jì)算公式如下：

[Re=Vbρμ] （3）

[Cd=Fd（0.5?ρV2S）] （4）

[Cl=Fl（0.5?ρV2S）] （5）

[K=ClCd] （6）

[Cm=Mz（0.5?ρV2Sb）] （7）

式中：V為來(lái)流速度（m·s-1）；b為網(wǎng)板弦長(zhǎng)（m）；ρ為流體密度（kgm-3）；μ為流體動(dòng)力粘度（kgm-1s-1）；Fd為網(wǎng)板阻力（N）；S為網(wǎng)板面積（m2）；Fl為網(wǎng)板擴(kuò)張力（N）；Mz為俯仰力矩（N·m）。

1.3 網(wǎng)板結(jié)構(gòu)參數(shù) 網(wǎng)板的基本尺寸如圖1所示，L為網(wǎng)板翼展，b為網(wǎng)板翼弦，f為網(wǎng)板的最大厚度，展弦比λ=L/b，迎流沖角α為網(wǎng)板翼弦與來(lái)流方向的夾角，網(wǎng)板模型如圖2所示。因?yàn)榫匦吻婢W(wǎng)板為對(duì)稱結(jié)構(gòu)，因此在計(jì)算時(shí)，可模擬網(wǎng)板的一半，將大大節(jié)省計(jì)算時(shí)間。

圖1 網(wǎng)板基本尺寸

圖2 網(wǎng)板模型

本文研究對(duì)象為2種不同展弦比的矩形曲面網(wǎng)板，迎流沖角α均為0°～50°范圍，網(wǎng)板的具體規(guī)格見(jiàn)表1。

表1 網(wǎng)板規(guī)格

[編號(hào) 翼展

（L/m） 翼弦

（b/m） 展弦比

（λ） 面積

（S/m2） 最大厚度

（f/m） 1 2.14 1.07 2.0 2.3 0.107 2 1.86 1.24 1.5 2.3 0.124 ]

1.4 模型參數(shù)設(shè)置 網(wǎng)板周圍的網(wǎng)格劃分如圖3所示，靠近網(wǎng)板的區(qū)域選擇六面體單元，遠(yuǎn)場(chǎng)區(qū)域選擇四面體單元來(lái)進(jìn)行計(jì)算，第一層網(wǎng)格高度為b/100。計(jì)算區(qū)域如圖4所示，水的密度ρ為998.2kg·m-3，動(dòng)力粘度μ為0.001003kg·m-1·s-1，上游入口處設(shè)為速度入口邊界（velocity inlet），速度為1.54m/s，方向沿x軸正向的均勻流，湍流強(qiáng)度為2.6%，湍流長(zhǎng)度尺度為0.07b。下游出口設(shè)為完全發(fā)展自由出流邊界（outflow）。計(jì)算區(qū)域的底面設(shè)為對(duì)稱邊界（symmetry），網(wǎng)板表面及計(jì)算域其他邊界均設(shè)為不可滑移壁面（no-slip wall）。

圖3 網(wǎng)板周圍網(wǎng)格劃分（α=20°）

圖4 計(jì)算域整體網(wǎng)格劃分

1.5 數(shù)值計(jì)算方法 經(jīng)驗(yàn)證，湍流模型選擇標(biāo)準(zhǔn)[k-ω]模型，控制方程采用有限體積法進(jìn)行離散，對(duì)流項(xiàng)采用二階迎風(fēng)格式，擴(kuò)散項(xiàng)采用二階中心差分格式，計(jì)算步長(zhǎng)為0.004s。

2 結(jié)果與分析

2.1 升力系數(shù)Cl和阻力系數(shù)Cd 網(wǎng)板的升力系數(shù)Cl和阻力系數(shù)Cd隨迎流沖角α的變化曲線如圖5所示。從圖5可以看出，網(wǎng)板的阻力系數(shù)Cd隨迎流沖角α的增大均呈上升趨勢(shì)。1號(hào)網(wǎng)板的升力系數(shù)Cl在α=0°～15°范圍內(nèi)隨著α的增大而增大，在α=15°時(shí)達(dá)到最大值Clmax=0.85，此時(shí)的Cd=0.17。在α=15°～20°范圍時(shí)，Cl有所下降，在α=20°時(shí)Cl=0.73。在α=20°～35°之間，隨著α的增大，升力系數(shù)有略微增加，當(dāng)α>35°后，隨著α的增大，升力系數(shù)呈下降趨勢(shì)。2號(hào)網(wǎng)板的Cl在α=0°～27°范圍內(nèi)先增后降，在α= 15°左右達(dá)到較大值為0.80，在α=27°～42°范圍內(nèi)，Cl持續(xù)上升，在α=42°時(shí)達(dá)到最大值Clmax=0.92，此時(shí)的Cd= 0.82。當(dāng)α>42°后，升力系數(shù)開(kāi)始下降。圖中Cl-exp和Cd-exp分別為2號(hào)網(wǎng)板的物理模型試驗(yàn)所得到的升力系數(shù)和阻力系數(shù)曲線[7]，從圖中可以看出，數(shù)值模擬得到的Cd曲線與物模試驗(yàn)相差不大，吻合較好。試驗(yàn)得到的最大升力系數(shù)Clmax=1.08（α=27°），數(shù)值模擬為Clmax=0.92（α=42°），兩者在數(shù)值上相差不大，但是所對(duì)應(yīng)的沖角有較大差距。產(chǎn)生誤差的原因有以下：（1）數(shù)值模擬沒(méi)有考慮摩擦力；（2）由于物模試驗(yàn)時(shí)間較早，且在風(fēng)洞中進(jìn)行，試驗(yàn)裝置和條件有限制，只在迎流沖角為5°～30°范圍內(nèi)進(jìn)行了試驗(yàn)，數(shù)據(jù)有限。從阻力系數(shù)曲線的擬合效果與升力系數(shù)曲線的變化范圍和趨勢(shì)來(lái)看，數(shù)值模擬的結(jié)果是正確的。

圖5 升力系數(shù)Cl和阻力系數(shù)Cd隨沖角α的變化曲線

2.2 升阻比K和俯仰力矩系數(shù)Cm 如圖6所示，升阻比K隨著迎流沖角α的增加呈先升后降趨勢(shì)。從圖7可以看出，1號(hào)網(wǎng)板的升阻比K較高，最大為K=9.59（α=5°），2號(hào)網(wǎng)板的最大升阻比K=8.35（α=5°），兩種網(wǎng)板最大升阻比時(shí)對(duì)應(yīng)的迎流沖角大致相同。但是，在最大升阻比對(duì)應(yīng)的沖角α=5°時(shí)，網(wǎng)板的升力系數(shù)并不高，無(wú)法充分?jǐn)U張網(wǎng)口，最大升力系數(shù)對(duì)應(yīng)沖角附近升阻比的大小也很重要。1號(hào)網(wǎng)板在α=15°（最大升力系數(shù)對(duì)應(yīng)沖角）時(shí)，升阻比K=5.05，2號(hào)網(wǎng)板在α=42°時(shí)，升阻比K=1.12。2號(hào)網(wǎng)板升阻比K的模擬結(jié)果與試驗(yàn)相比較，兩者在α=5°時(shí)完全吻合，隨著α的增大，K均呈下降趨勢(shì)，數(shù)值模擬結(jié)果與試驗(yàn)結(jié)果的變化趨勢(shì)是一致的。

圖6 升阻比K和俯仰力矩Cm隨迎流沖角α的變化曲線

俯仰力矩也稱為縱向力矩，是指作用在網(wǎng)板的流體動(dòng)力對(duì)其重心所產(chǎn)生的力矩沿橫軸的分量，一般轉(zhuǎn)換為俯仰力矩系數(shù)Cm來(lái)分析比較。俯仰力矩系數(shù)有正負(fù)之分，其絕對(duì)值大小表示俯仰力矩的高低，一般Cm趨于0表示該網(wǎng)板的穩(wěn)定性越優(yōu)。在α=5°左右時(shí)，兩塊網(wǎng)板均達(dá)到最大升阻比K，此時(shí)1號(hào)網(wǎng)板Cm的絕對(duì)值為0.013，2號(hào)網(wǎng)板Cm的絕對(duì)值為0.001。在α=0°～20°范圍內(nèi)，1號(hào)網(wǎng)板Cm的絕對(duì)值始終比2號(hào)網(wǎng)板大，意味著2號(hào)網(wǎng)板的穩(wěn)定性較好。在α>25°后，兩塊網(wǎng)板Cm的絕對(duì)值均小于0.10，比較穩(wěn)定。

3 結(jié)論

本文利用數(shù)值模擬的方法，研究了展弦比對(duì)矩形曲面網(wǎng)板水動(dòng)力學(xué)性能的影響。模擬結(jié)果與物理模型試驗(yàn)結(jié)果吻合較好，驗(yàn)證了模擬的正確性。有以下結(jié)論：

（1）展弦比為2.0的1號(hào)網(wǎng)板的臨界沖角α=15°，Clmax=0.85，Cd=0.17。展弦比為1.5的2號(hào)網(wǎng)板的臨界沖角α=42°，Clmax=0.92，Cd=0.82。因此展弦比越大，臨界沖角越小，對(duì)應(yīng)的Clmax和Cd越小。

（2）1、2號(hào)網(wǎng)板的最大升阻比K分別為9.59和8.35，在迎流沖角α=5°左右時(shí)達(dá)到，此時(shí)的俯仰力矩系數(shù)Cm的絕對(duì)值分別為0.013和0.001。在各自最大升力系數(shù)處，1、2號(hào)網(wǎng)板的升阻比K分別為5.05和1.12，俯仰力矩系數(shù)的絕對(duì)值分別為0.133和0.094。因此，展弦比越大，升阻比K越大，但穩(wěn)定性越差。

參考文獻(xiàn)

[1]張勛，王明彥，徐寶生.拖網(wǎng)網(wǎng)板型式、結(jié)構(gòu)與性能的研究與應(yīng)用進(jìn)展[J].中國(guó)水產(chǎn)科學(xué)，2004，11：107-113.

[2]徐寶生，張勛，王明彥.單船拖網(wǎng)網(wǎng)板的現(xiàn)狀及發(fā)展趨勢(shì)[J].福建水產(chǎn)，2010（1）：86-90.

[3]劉健，黃洪亮，吳越，等.2種立式曲面縫翼式網(wǎng)板水動(dòng)力學(xué)性能的實(shí)驗(yàn)研究[J].南方水產(chǎn)科學(xué)，2015，11（1）：68-74.

[4]Sala A，F(xiàn)arran JDP，Antonijuan J，et al. Performance and impact on the seabed of an existing- and an experimental-otterboard： Comparison between model testing and full-scale sea trials[J]. Fisheries Research，2009，100（2）：156-166.

[5]Takahashi Y，F(xiàn)ujimori Y，Hu F，et al. Design of trawl otter boards using computational fluid dynamics[J]. Fisheries Research，2015，161：400-407.

[6]Balash C，Sterling D，Lustica M，et al. Twist-and-camber effects on the performance of simple hydro-sails for efficiently spreading penaeid trawls[J]. Ocean Engineering，2015，109：161-168.

[7]郭根喜，劉同渝，黃小華.拖網(wǎng)網(wǎng)板動(dòng)力學(xué)理論研究與實(shí)踐[M]. 廣州：廣東省出版集團(tuán)，2008.

[8]張勛，王錦浩，王明彥，等.矩形V型曲面網(wǎng)板水動(dòng)力性能的研究——Ⅰ導(dǎo)流板曲率與網(wǎng)板水動(dòng)力性能的關(guān)系[J]. 中國(guó)水產(chǎn)科學(xué)，2004，11（s1）：5-8.

[9]王錦浩，張勛，王明彥，等.矩形V型曲面網(wǎng)板水動(dòng)力性能的研究——Ⅱ展弦比、縫口位置、縫口寬度對(duì)網(wǎng)板水動(dòng)力性能的影響[J]. 中國(guó)水產(chǎn)科學(xué)，2004，11（z1）：9-13.

[10]王磊，王魯民，馮春雷，等.葉板尺度比例變化對(duì)單縫曲面網(wǎng)板水動(dòng)力性能的影響[J]. 漁業(yè)現(xiàn)代化，2015，42（6）：55-60.

[11]劉健，黃洪亮，陳帥，等.小展弦比立式曲面網(wǎng)板的水動(dòng)力性能[J]. 水產(chǎn)學(xué)報(bào)，2013，37（11）：1742-1749.

（責(zé)編：張宏民）