耿寶磊，王榮泉，寧德志，張慈珩，趙 鵬

（1.交通運輸部天津水運工程科學研究所港口水工建筑技術國家工程實驗室工程泥沙交通行業(yè)重點實驗室，天津300456；2.大連理工大學海岸與近海工程國家重點實驗室，大連116024）

波浪與直墻前多層透空薄板作用的解析研究

耿寶磊1，王榮泉2，寧德志2，張慈珩1，趙 鵬1

（1.交通運輸部天津水運工程科學研究所港口水工建筑技術國家工程實驗室工程泥沙交通行業(yè)重點實驗室，天津300456；2.大連理工大學海岸與近海工程國家重點實驗室，大連116024）

文章基于勢流理論下均勻透空薄板內部流體連續(xù)方程和忽略對流項的歐拉方程，得到薄板內流體水平速度與壓力差間的關系，進而建立了波浪與直墻前單層和多層透空薄板作用模型，通過匹配薄板兩側的邊界條件，可求得直墻前透空薄板對波浪的反射系數(shù)、透射系數(shù)以及波能的吸收率。算例針對水深2.5 m，入射波浪周期3.5 s的工況，研究了直墻前12 m范圍內設置不同層數(shù)透空薄板時，反射系數(shù)、透射系數(shù)及波能吸收率隨組合參數(shù)G的變化，結果表明當布置四層透空薄板且0.8＜Gr＜4.6、0＜Gi＜1.6時，本算例波浪的吸收率可達90%以上。

波浪；多層透空薄板；吸收率；反射系數(shù)；透射系數(shù)

開孔直立結構具有消能效果好、波壓力小等優(yōu)點，在直立碼頭、防波堤、水槽消波等實際工程中得到了廣泛應用。國內外很多學者對波浪與開孔結構作用進行了很多研究。在試驗研究方面，杜金曼和吳承元［1］研究了多層直立消波裝置的反射系數(shù)和透射系數(shù)的影響因素。戴冠英［2］、張芹和戴冠英［3］開展模型實驗研究了開孔率較小時（0.15＜ε＜0.25）直立全開孔板或全開孔板后帶實體墻兩種結構形式的波浪反射系數(shù)及結構所受的波壓力變化規(guī)律，并給出了反射系數(shù)、透射系數(shù)及波壓力與結構形式、水深、波要素之間的回歸關系式。馬寶聯(lián)等［4］和Tanimoto等［5］分別針對開孔沉箱結構進行了試驗研究，分析了反射系數(shù)與消浪室相對寬度和開孔率等影響因素之間的關系。

在理論研究與數(shù)值分析方面，Yu［6］基于線性勢流理論研究了波浪對透水半無限防波堤的繞射作用。Chwang和Chan［7］基于達西定律分析了波浪與透空結構的作用，分析了透空結構對減小港內波浪和港內共振的影響。Twu和Lin［8］等采用解析的方法研究了多塊透空板不同布置間距對波浪的反射的影響，并發(fā)現(xiàn)前面板的開孔率較大，后面板開孔率較小對波浪有較好的吸收效果。李世森等［9］在波浪數(shù)值水槽模型末端采用沿水槽縱軸方向具有漸變垂向孔隙率的孔隙結構消除反射波，并分析了不同因素對消波性能的影響。陳雪峰等［10］利用數(shù)值水槽對波浪與開孔板作用進行了數(shù)值分析，研究了開孔率等因素對開孔板點壓力分布和反射系數(shù)的影響。

本文基于勢流理論，建立了波浪與直墻前多層透空薄板作用的數(shù)學模型，利用解析方法得到了多層透空薄板對波浪的反射系數(shù)、透射系數(shù)及波能的吸收率，進一步通過具體算例分析了組合參數(shù)G對透空薄板消波性能的影響。

1 透空薄板內流體運動控制方程

圖1為單層透空薄板示意圖，模型采用勢流理論，在均勻薄板內部流體流動滿足連續(xù)方程和忽略對流項的歐拉方程如下［11］：

式中：U為流體的速度矢量；r為流體密度；P為流體壓強，t為時間；f為薄板的阻力系數(shù)；w為入射波浪頻率；Cm為附加質量；e為薄板開孔率。

對于速度U和壓強P，按下式分離出時間因子

可得其復變量u和p的方程為

式中

其實部和虛部分別對應介質的阻力和慣性影響。

假設薄壁只橫向透水，忽略掉流體的垂向分量，可得到薄板內水平速度與壓力差間的關系為

式中b為薄板的厚度，下標0和b分別表示薄板兩側的物理量。

將孔隙中的流速換算成整個板上的流速，并與外部波動速度相匹配，得

通過對不同情況下速度勢與邊界條件的匹配，可以得出直墻前單層和多層透空薄板模型對波浪的反射系數(shù)、透射系數(shù)以及波浪能量損失的表達式。

圖1 透空薄板模型示意圖Fig.1 Sketch of grid plate model

圖2 波浪與直墻前單層透空薄板作用Fig.2 Wave action with one grid plate in front of wall

2 波浪與直墻前單層透空薄板作用

對于圖2所示的波浪與直墻前單層透空薄板的作用問題，在W1和W2區(qū)域內，速度勢f1和f2應滿足以下條件：

自由水面條件

式中g為重力加速度。

水底條件

無窮遠條件

式中φ0為入射勢，k為波數(shù)，i為虛數(shù)單位。

直墻上物面條件

根據(jù)式（8）和式（9），并令

其中，實部Gr主要跟介質引起的阻力系數(shù)有關，虛部Gi主要跟透空介質引起的附加質量系數(shù)有關。式（9）可寫為

忽略透空薄板的厚度，在薄板處的速度連續(xù)條件可寫為

對于匹配條件式（15）和式（16），φ1和φ2中不包括非傳播項［6］，這樣φ1和φ2可寫為

式中A為波幅；h為水深；R和T為待定系數(shù)。將φ1和φ2代入（13）、（14）式可解得

則單層透空薄板的反射系數(shù)Kr和透射系數(shù)Kt分別為

波浪穿過單層透空薄板的能量損失系數(shù)為

將波浪與直墻前透空薄板作用時的能量損失和透射波浪作為入射波浪的消減部分，則透空薄板對波能的吸收率定義為

圖3 波浪與直墻前多層透空薄板作用示意圖Fig.3 Wave action with multi-layer perforated plates in front of the wall

3 波浪與直墻前多層透空薄板作用

按照同樣的方法可以求得波浪與直墻前多層透空薄板作用的情況，圖3是波浪與直墻前兩層、三層和四層透空薄板作用的示意圖。

3.1 波浪與直墻前兩層透空薄板作用

以波浪與直墻前兩層透空薄板作用為例，兩層薄板內水平速度與壓力差間的關系分別為

式中：lmn表示薄板的間距，其中m表示柵格層數(shù)，n表示第n和n+1層的間距；分別將兩層薄板孔隙中的流速換算成整個板上的流速，并與外部波動速度相匹配，得

同樣令G=ε/Ωkb=Gr+iGi，（23）式可寫為

忽略透空薄板的厚度，則在兩層薄板處的速度連續(xù)條件可寫為

速度勢φ1、φ2和φ3的表達式為

式中R1和T1分別表示第一層透空薄板復變量的反射系數(shù)和透射系數(shù)；R2和T2分別表示第二層透空薄板復變量的反射系數(shù)和透射系數(shù)。此外，由于不透水直墻對波浪發(fā)生全反射，反射系數(shù)仍為T2。

將速度勢φ1、φ2和φ3代入上述匹配條件式（24）和（25）可求得

式中

由此可得直墻前兩層透空薄板對波浪的反射系數(shù)Kr和透射系數(shù)Kt分別為

直墻前雙層透空薄板對波能的吸收率同式（21）。

3.2 波浪與直墻前三層透空薄板作用

波浪與直墻前三層透空薄板作用時，各層薄板復變量的反射系數(shù)和透射系數(shù)表達式為

3.3 波浪與直墻前四層透空薄板作用

波浪與直墻前四層透空薄板作用時，各層薄板復變量的反射系數(shù)和透射系數(shù)表達式為

由此可得直墻前四層透空薄板對波浪的反射系數(shù)Kr和透射系數(shù)Kt分別為

4 算例分析

根據(jù)式（7）和（14），可以得到G的表達式為

為研究反射系數(shù)、透射系數(shù)和波能吸收率與G的關系，選取特定的算例進行計算，其中水深h=2.5 m，入射波浪周期為3.5 s，第一層薄板距離直墻的距離為12.0 m。

（1）布置單層透空薄板

布置單層透空薄板時，l11=12.0 m。計算結果見圖4。

圖4 單層透空薄板的反射系數(shù)、透射系數(shù)與波能吸收率變化曲線Fig.4 Reflection coefficient,transmission coefficient and absorption rate for one perforated plate

根據(jù)式（36），當附加質量Cm取為0時，G的實部Gr和虛部Gi與開孔率ε成正比。根據(jù)圖4中Kr的變化曲線，當Gr和Gi較小或較大時，反射系數(shù)都很大，可以解釋為當Gr和Gi較小時，開孔率ε較小，此時透空薄板相當于不透水直墻，故而反射系數(shù)很大；而當Gr和Gi較大時，開孔率ε較大，此時大部分波浪能透過透空板，但是，波浪在后面直墻發(fā)生全反射進而又透過透空板進入入射波場，相當于透空板發(fā)生反射，因此Gr和Gi較大時反射系數(shù)也很大。

同理，對于圖4中Kt的變化曲線，當Gr和Gi較小時，開孔率ε較小，此時透空薄板相當于不透水直墻，故而透射系數(shù)很小；而當Gr和Gi較大時，開孔率ε較大，此時大部分波浪能透過薄板，故透射系數(shù)很大。

從圖4中Ka的變化曲可以看到，當Gr在0.6～1.6之間取值，Gi在0～0.3之間取值時，能使得波浪能量的吸收率達到90%以上，只是Gr和Gi的有效取值范圍較小。

（2）布置雙層透空薄板

布置雙層透空薄板時，l21=l22=6.0 m，計算結果見圖5。從Ka的變化曲線可以看到，當Gr在0.8～2.2之間取值、Gi在0～0.35以內取值時，波浪吸收率能達到90%以上。波浪通過雙層薄板后也只在較小的區(qū)域能達到好的消波效果，但相對單層薄板的情況略有改善。

（3）布置三層透空薄板

布置三層透空薄板時，l31=l32=l33=4.0 m，計算結果見圖6。從Ka的變化曲線可以看到，當Gr在0.7～3.5之間取值、Gi在1.0以內取值時，波浪吸收率能達到90%以上。波浪通過三層薄板后，消波效果相對前兩種布置方式有較為明顯的改善。

（4）布置四層透空薄板

布置四層透空薄板時，l41=l42=l43=l44=3.0 m，計算結果見圖7。從Ka的變化曲線可以看到，當Gr在0.8～4.6之間取值、Gi在1.6以內取值時，波浪吸收率能達到90%以上。波浪通過四層薄板后，消波效果較之前三種情況有明顯增強。

通過圖4～圖7的對比可以發(fā)現(xiàn)，隨著直墻前透空薄板數(shù)量的增加，能使波能吸收率大于90%的Gr和Gi取值范圍在不斷增加，也就是說，通過增加透空薄板的數(shù)量能有效提高系統(tǒng)對波浪的吸收，減小波浪的反射。

圖5 雙層透空薄板的反射系數(shù)、透射系數(shù)與波能吸收率變化曲線Fig.5 Reflection coefficient,transmission coefficient and absorption rate for two perforated plates

圖6 三層透空薄板的反射系數(shù)、透射系數(shù)與波能吸收率變化曲線Fig.6 Reflection coefficient,transmission coefficient and absorption rate for three perforated plates

圖7 四層透空薄板的反射系數(shù)、透射系數(shù)與波能吸收率變化曲線Fig.7 Reflection coefficient,transmission coefficient and absorption rate for four perforated plates

5 結論

本文基于勢流理論建立了波浪與直墻前多層透空薄板作用模型，得到了多層透空薄板裝置的反射系數(shù)、透射系數(shù)以及對波能的吸收率，進一步通過具體算例分析了組合參數(shù)G的變化對反射系數(shù)、透射系數(shù)及吸收率的影響，研究結果表明當0.8＜Gr＜4.6、0＜Gi＜1.6時，本算例波浪的吸收率可達90%以上。

組合參數(shù)G包含了薄板的孔隙率、薄板寬度、薄板的阻力系數(shù)、附加質量以及入射波浪的波數(shù)等，進一步的研究將針對不同參數(shù)的變化給出直墻前多層透空薄板裝置波能吸收率變化的敏感性，研究成果可為新型消波裝置的開發(fā)提供借鑒。

參考文獻：

［1］杜金曼,吳承元.多層透空直立式消波器［J］.水利水運科學研究,1984（4）：79-88. DU J M,WU C Y.Multi-layer perforated vertical wave absorber［J］.Hydro-Science and Engineering,1984（4）：79-88.

［2］戴冠英.波浪作用下開孔直立結構的反射與透射性能［J］.水利水運科學研究,1993（3）：291-300. DAI G Y.Reflection and transmission performances of vertical perforated structure under wave action［J］.Hydro-Science and Engineering,1993（3）：291-300.

［3］張芹,戴冠英.波浪對開孔直立結構作用力的試驗研究［J］.水利水運科學研究,1994（4）：367-374. ZHANG Q,DAI G Y.Experimental study on interaction between waves and vertical perforated structures［J］.Hydro-Science and Engineering,1994（4）：367-374.

［4］馬寶聯(lián).波浪與開孔直墻式防波堤的相互作用［D］.大連：大連理工大學,2004.

［5］Tanimoto K,Yoshimoto Y.Theoretical and experimental study of reflection coefficient for wave dissipating caisson with a permeable front wall［J］.Report of the Port and Harbour Research Institute,1982,21（3）：44-77.

［6］Yu X.Diffraction of water waves by porous breakwaters［J］.Journal of waterway,port,coastal,and ocean engineering,1995,121（6）：275-282.

［7］Chwang A T,Chan A T.Interaction between porous media and wave motion［J］.Annual Review of Fluid Mechanics,1998,30（1）：53-84.

［8］Twu S W,Lin D T.Wave reflection by a number of thin porous plates fixed in a semi-infinitely long flume［J］.Coastal Engineering Proceedings,1990,1（22）：21-23.

［9］李世森,江鳴.數(shù)值波浪水槽末段孔隙結構消波性能的研究［J］.港工技術,2012（49）：1-4. LI S S,JIANG M.Study on Wave-absorbing Performance of Porosity Solid in End-piece of Numerical Wave Flume［J］.Port Engineering Technology,2012（49）：1-4.

［10］陳雪峰,張梅,李玉成.開孔率對開孔板消浪效果影響的數(shù)值模擬研究［J］.水道港口,2013（34）：285-292. CHEN X F,ZHANG M,LI Y C.Numerical simulation study on effect of porosity on performance of perforated structures［J］.Journal of Waterway and Harbor,2013（34）：285-292.

［11］李玉成,滕斌.波浪對海上建筑物的作用［M］.北京.海洋出版社出版社，2015.

Analytical study on wave interaction with multi-layer perforated plates in front of vertical wall

GENG Bao-lei1,WANG Rong-quan2,NING De-zhi2,ZHANG Ci-heng1,ZHAO Peng1
（1.Tianjin Research Institute for Water Transport Engineering，National Engineering Laboratory for Port Hydraulic Construction Technology,Key Laboratory of Engineering Sediment,Ministry of Transport,Tianjin300456,China;2. State Key Laboratory of Coastal and Offshore Engineering,Dalian University of Technology,Dalian116024,China）

Based on potential flow theory,an analytical study on wave action with one grid plate or multi-grid plates was carried out.Continuity equation and Euler equation without convective term were used to get the relationship between horizontal velocity and pressure difference of the fluid in the thin plate.By matching the boundary conditions at two sides of the plate,reflection coefficient,transmission coefficient and absorption rate were developed. Examples with water depth 2.5 m,incident wave period 3.5 s and wave absorbing distance 12 m were carried out to study the effects of different plates numbers with an assembly parameterG,and the results show that the absorption rate of four-layer perforated plates is bigger than 90%when 0.8＜Gr＜4.6 and 0＜Gi＜1.6.

wave;multi-grid plates;absorption rate;reflection coefficient;transmission coefficient

U 656.2；O 242.1

A

1005-8443（2017）01-0008-08

2016-06-17；

2016-07-27

交通運輸部應用基礎研究計劃項目（2014329224380）；國家自然科學基金（51409135）；天津市應用基礎與前沿技術研究計劃（15JCQNJC07300）；中央級科研院所基本科研業(yè)務費專項（TKS150203）

耿寶磊（1980-），男，河北省衡水人，副研究員，主要從事波浪理論及波浪與結構物作用研究。

Biography：GENG Bao-lei（1980-），male，associate professor.