• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    波流復合工況下縫邊式圍欄網(wǎng)衣水動力特性

    2022-02-04 12:10:52馮德軍桂福坤
    農(nóng)業(yè)工程學報 2022年19期
    關(guān)鍵詞:網(wǎng)衣波流網(wǎng)線

    王 順,馮德軍,桂福坤

    波流復合工況下縫邊式圍欄網(wǎng)衣水動力特性

    王順,馮德軍,桂福坤※

    (浙江海洋大學國家海洋設(shè)施養(yǎng)殖工程技術(shù)研究中心,浙江省海洋養(yǎng)殖裝備與工程技術(shù)重點實驗室,舟山 316022)

    圍欄養(yǎng)殖是一種發(fā)展迅速的新型生態(tài)養(yǎng)殖模式。網(wǎng)衣系統(tǒng)是圍欄養(yǎng)殖設(shè)施的核心組成部分,與設(shè)施的安全密切相關(guān)。該研究針對縫邊式圍欄網(wǎng)衣系統(tǒng),采用集中質(zhì)量點法研究其在波流復合工況下的受力特性。研究結(jié)果表明:在波高4~6 m、流速1.0 m/s、入射角度0°時,綱繩最大受力出現(xiàn)在距水面下1~2.5 m處的第二、三根橫向綱繩兩端處;網(wǎng)線最大受力出現(xiàn)在水面下第二根橫向綱繩與中間兩根豎向綱繩連接處;網(wǎng)衣受力呈現(xiàn)對稱分布,整片網(wǎng)衣的載荷主要來自上半部分網(wǎng)衣,約占總荷載的80%;隨著入射角度的增大,綱繩和網(wǎng)線的最大受力以及網(wǎng)衣的最大偏移量逐漸減?。煌炔ǜ邨l件下,入射角度0°~40°時,圍欄網(wǎng)衣受力及最大偏移量下降速率較慢,入射角度40°~80°時,下降速率較快;位于兩側(cè)樁柱上的系縛點的最大受力出現(xiàn)在水面下第二、三根橫向綱繩兩端處。在工程運用過程中,建議將以上最大受力部位進行加固處理,以提升圍欄網(wǎng)衣抗破壞能力。該研究結(jié)果可為圍欄實際工程建設(shè)和維護階段提供技術(shù)支撐。

    水產(chǎn)養(yǎng)殖;模型;力學性能;網(wǎng)衣;波流復合工況

    0 引 言

    圍欄養(yǎng)殖工程,因其有養(yǎng)殖面積大、養(yǎng)殖對象品質(zhì)更近生態(tài)等諸多優(yōu)點備受養(yǎng)殖企業(yè)青睞[1]。起初在池塘、湖泊、沿海灘涂用毛竹插入泥沙底部的方式進行養(yǎng)殖[2],然后由內(nèi)海向優(yōu)質(zhì)水體的外海推進,出現(xiàn)浮繩式、圍欄等結(jié)構(gòu)形式。但是目前所有的圍欄養(yǎng)殖中都存在不同程度的網(wǎng)衣破壞,無論是在網(wǎng)衣強度設(shè)計、結(jié)構(gòu)設(shè)計,還是后期維護階段,圍欄網(wǎng)衣設(shè)施安全性能評估成為當今研究熱點。

    樁柱和網(wǎng)衣是圍欄養(yǎng)殖工程的兩個核心構(gòu)成單元。樁柱一般可參照《港口工程樁基規(guī)范》[3]或《建筑樁基技術(shù)規(guī)范》[4]進行設(shè)計和建造。圍欄養(yǎng)殖工程的建造是在深水抗風浪網(wǎng)箱投入應(yīng)用之后,最近10年迅速發(fā)展的水產(chǎn)養(yǎng)殖產(chǎn)業(yè),抗風浪網(wǎng)箱網(wǎng)衣的研究發(fā)展為圍欄養(yǎng)殖工程網(wǎng)衣提供了方法和手段。最初的研究以試驗[5]為主,后又引入數(shù)值模擬的方法[6-9],建立了通用的網(wǎng)箱的數(shù)值模擬方法和程序[10-14],并被廣泛應(yīng)用于網(wǎng)箱浮架、網(wǎng)衣、錨碇等系統(tǒng)的水動力特性研究。

    圍欄養(yǎng)殖工程和浮繩式圍網(wǎng)都是圈圍一定量的海面,上至海面、下至海底,但是在結(jié)構(gòu)形式上,浮繩式圍網(wǎng)采用浮球支撐網(wǎng)的重力,圍欄則是利用樁柱固定網(wǎng)衣,在支撐網(wǎng)衣上存在一定的不同。葉衛(wèi)富[15]通過模型水槽試驗對浮繩式圍網(wǎng)進行了研究,初步得到圍網(wǎng)設(shè)施水動力變化規(guī)律,但只研究一個入射方向上的水流和波浪工況,也未對圍網(wǎng)形狀變化進行研究,且僅僅針對框架綱和系泊綱進行研究。崔勇等[16]通過物理試驗和數(shù)值模擬手段,對圍網(wǎng)浮子的最大位移和網(wǎng)衣張力進行研究。徐克品等[17]采用正弦波理論和莫里森方法理論對浮繩式圍網(wǎng)網(wǎng)衣進行研究,得到局部網(wǎng)衣在垂直方向的受力情況。以上均為針對浮繩式圍網(wǎng)進行的研究。針對圍欄網(wǎng)衣水動力研究,桂福坤等[18-19]采用集中質(zhì)量點法建立了網(wǎng)衣的運動響應(yīng)數(shù)學模型,得到固定方式和綱繩直徑、網(wǎng)格形狀對圍欄網(wǎng)衣力學性能的影響,但為了研究方便,將網(wǎng)衣全部放置水底,網(wǎng)衣最高處置于波浪波谷處,即波浪作用于網(wǎng)衣時不露出水面。陳天華等[20-21]采用數(shù)值模擬的手段,分析單因素(波高、流速)對網(wǎng)衣的影響,這里需要指出其是針對四邊加綱繩形式的單元網(wǎng)衣進行研究,且也為了研究方便網(wǎng)衣全部置于水底。楊熙等[22]研究了波浪載荷作用下網(wǎng)片出水和完全淹沒條件的水動力特性,分析得出網(wǎng)衣出水與不出水在受力、位移等方面存在較大差異,其研究也是單因素波高對網(wǎng)衣的影響。

    綜上所述,以上研究均存在工況單一等問題,均未考慮波流復合工況及其入射角度對其的影響。該研究針對波流復合工況(波高+流速),結(jié)合不同入射角對圍欄網(wǎng)衣力學性能進行系統(tǒng)的研究。并以實際工程圍欄網(wǎng)衣為研究對象,以確定圍欄網(wǎng)衣的最大受力、最大偏移、系縛點最大受力,為圍欄設(shè)計安裝以及后期維護等提供技術(shù)支撐。

    1 材料和方法

    1.1 圍欄養(yǎng)殖工程原模型

    該研究以實際工程浙江省鹿西島圍欄網(wǎng)衣布置為研究背景(圖1)。該圍欄養(yǎng)殖工程主要由樁柱、網(wǎng)衣、系網(wǎng)桿、抱箍等組成,其特點是長方形網(wǎng)衣固定在系網(wǎng)桿上,系網(wǎng)桿通過抱箍固定在樁柱上。網(wǎng)衣縱向有4根間隔1.2~1.5 m豎綱繩固定,其中左右兩端為了防止養(yǎng)殖魚類的逃逸,通過固定繩固定在系網(wǎng)桿上。網(wǎng)衣橫向根據(jù)水深分布每隔1 m左右布置一根綱繩,上方通過固定繩固定。為了防止養(yǎng)殖魚類從網(wǎng)衣底部逃逸,將鐵鏈加裝在底部綱繩上并嵌入海底。

    1.樁柱 2.抱箍 3.系網(wǎng)桿 4.網(wǎng)線 5.豎向綱繩 6.橫向綱繩 7.系縛點 8.鐵鏈 9.固定繩

    1.2 數(shù)值模擬方法

    圍欄網(wǎng)衣的綱繩和網(wǎng)線屬于柔性結(jié)構(gòu)物,具有小尺度、多自由度、大變形的特點,所以采用集中質(zhì)量點的方法[24]進行模擬。通過計算集中質(zhì)量點在波浪和邊界條件下的運動,來確定網(wǎng)衣的變形,進而確定網(wǎng)衣的運動。網(wǎng)線和綱繩的集中質(zhì)量點分別設(shè)于網(wǎng)目目腳的兩端和綱繩與網(wǎng)線結(jié)節(jié)接觸處,該研究集中質(zhì)量點的計算主要分為4類(圖2a)。以1號結(jié)節(jié)處集中質(zhì)量點為例,網(wǎng)衣集中質(zhì)量點受力主要包括:重力、浮力、張力、波浪力(速度力和慣性力)。1號結(jié)節(jié)處所受重力和浮力由一個結(jié)節(jié)、兩個目腳、一根橫綱和一根豎綱提供(圖2b),其受力分析如圖2c所示。

    注:B為集中質(zhì)量點的浮力,N;W為集中質(zhì)量點的重力,N;FD、FI分別表示集中質(zhì)量點的速度力和慣性力,N;T1、T2、T3、T4為綱繩作用于結(jié)節(jié)的張力,N;T5、T6、T7、T8為網(wǎng)線作用在結(jié)節(jié)上的張力,N;1~4為網(wǎng)衣的4類集中質(zhì)量點。

    集中質(zhì)量點張力其主要計算式[23-24]如下:

    式中為此段網(wǎng)線(或綱繩)直徑,m;0為網(wǎng)線原始長度,m;為變形后的長度,m;(x,y,z)、(x,y,z)分別表示該集中質(zhì)量點的初始位置和移動后的坐標,m;1、2為構(gòu)件材料彈性系數(shù);

    網(wǎng)衣運動主要由波浪力提供,波浪力的計算可以采用莫里森方程[25]進行計算。根據(jù)網(wǎng)衣的結(jié)構(gòu)特點,網(wǎng)線的目腳和綱繩兩端集中質(zhì)量點為圓球,圓球在各個方向的投影面積相等,故水動力系數(shù)取恒定值。網(wǎng)線的目腳和綱繩假設(shè)為圓柱體桿件,桿件在運動過程各個方向投影有差別,故需要考慮波浪入射方向與網(wǎng)線夾角之間的關(guān)系,故需在目腳和綱繩上建立局部坐標系(,,),由于目腳和綱繩的非結(jié)節(jié)處波浪力計算相同,取綱繩的一段,受力分析如圖3所示。

    注:τ軸方向為沿綱繩方向,ξ軸在τ軸和水質(zhì)點相對速度組成的平面內(nèi)與τ軸垂直,η軸與ξ軸和水質(zhì)點組成的平面互相垂直;V為波流共同作用的水質(zhì)點速度,m?s-1;q為波流入射方向與y軸的夾角,(°);其他變量含義與圖2相同。

    該段綱繩在波流復合工況下波浪力主要計算式[25]如下:

    波浪水質(zhì)點速度的主要速度勢計算式[25]如下:

    式中表示水質(zhì)點速度勢;表示波速,m/s;表示波浪波數(shù);表示波浪角頻率,rad/s;表示集中質(zhì)量點深度(z方向)的相對位置,m;表示水深,m;表示集中質(zhì)量點在方向的位置,m;表示水質(zhì)點運動的時間,s;表示由(波高)與?決定的參數(shù);表示波長,m;123表示式(5)累積項F,N。

    在整體坐標系下將各集中質(zhì)量點所包含的結(jié)節(jié)和目腳的受力進行累加,并將其分配到集中質(zhì)量點上,然后利用牛頓第二定律建立集中質(zhì)量點運動方程:

    式中、Δ分別表示各集中質(zhì)量點的質(zhì)量和附加質(zhì)量,kg;表示集中質(zhì)量點的加速度,m/s2。

    利用Fortran軟件進行程序編寫,給定網(wǎng)衣的初始狀態(tài),計算各個Δ時間,網(wǎng)衣質(zhì)點運動的位移、速度和加速度等,最后利用MATLAB軟件對結(jié)果進行可視化分析。

    1.3 模型驗證

    為了驗證網(wǎng)衣數(shù)學模型的正確性和有效性,在浙江海洋大學拖曳力水槽進行驗證性試驗,具體試驗過程同文獻[20]和[23]。拖曳力水槽長度130 m、寬度6 m、深度3.44 m。本試驗中網(wǎng)衣的材料采用聚乙烯(PE)材質(zhì),網(wǎng)的尺寸0.64 m′0.20 m。試驗水深0.7 m,波浪周期在1.6 s,由于造波機最大波高0.35 m,圖4給出波高0.05~0.25 m網(wǎng)衣最大水平波浪力的計算結(jié)果??傻镁W(wǎng)衣水平波浪力最大值的模擬值與試驗值誤差在4.9%,說明該研究數(shù)值模擬計算的結(jié)果可真實反映網(wǎng)衣受力特性。

    1.4 計算參數(shù)及試驗設(shè)計

    實際工程結(jié)構(gòu)中,圍欄網(wǎng)衣有成千上萬個網(wǎng)目,計算機難以處理眾多網(wǎng)目單元形成的計算矩陣,所以采用網(wǎng)目群化的方法。計算原始網(wǎng)目數(shù)為13 824,采用“8×8網(wǎng)目群化”方法[26-27],群化后網(wǎng)目數(shù)為216。圍欄網(wǎng)衣尺寸:3.80 m×5.76 m,由從上而下5根橫向綱繩,從左至右4根豎向綱繩和網(wǎng)線組成,水平縮結(jié)系數(shù)0.66,垂直縮結(jié)系數(shù)0.75,采用材料具體的參數(shù)如表1所示。結(jié)節(jié)水動力參數(shù)[28-29]?。?i>C=1.45,C=0.50;網(wǎng)線水動力參數(shù)取:C=1.45,C=1.45,C=1.45;C=0.20,C=0.20,C=0.20。網(wǎng)衣的運動狀態(tài)是三維的,z方向為網(wǎng)衣集中質(zhì)量點位置信息,運動方向為z方向,初始狀態(tài)下,軸的坐標為0,將網(wǎng)衣露出水面0.5 m,如圖5所示。該研究通過比較波面水質(zhì)點位置()與網(wǎng)衣集中質(zhì)量點的位置()[30]判斷網(wǎng)衣是否出水。當時,水質(zhì)點的速度和加速度都為0;當時,運用式(5)的速度勢公式計算水質(zhì)點速度和加速度。根據(jù)道式波級[31],波高分為十個等級,在中國近海有效波高4.0 m以上就定義為災(zāi)害性波浪,所以設(shè)計波高4、5、6 m,流速1.0 m/s,根據(jù)圍欄養(yǎng)殖工程建設(shè)的位置,將圍欄可能面臨的海浪拍打方向即入射角(式4)分為0°、10°、40°、60°、80°,組合共設(shè)計15組組合工況。波浪和流速同向,其中0°代表波浪垂直拍打圍欄網(wǎng)衣。

    圖4 網(wǎng)衣最大水平波浪力試驗裝置及模擬值和試驗值對比

    表1 網(wǎng)衣參數(shù)

    圖5 網(wǎng)衣運動的初始狀態(tài)

    該研究主要通過數(shù)值模擬的方法,針對圍欄網(wǎng)衣水動力特性的4個方面展開研究:1)研究圍欄網(wǎng)衣綱繩和網(wǎng)線的最大受力;2)研究圍欄網(wǎng)衣的最大偏移量(最大偏移量是指網(wǎng)衣變形后結(jié)節(jié)到網(wǎng)衣初始平面的最大距離);3)研究圍欄網(wǎng)衣網(wǎng)線的受力。4)研究圍欄網(wǎng)衣系縛點的受力。

    2 結(jié)果與分析

    2.1 入射角度對圍欄網(wǎng)衣最大受力的影響

    波高4 m、流速1.0 m/s、入射角0°、10°、40°、60°、80°時,網(wǎng)衣的受力如圖6所示。從圖中可以直接看出隨著入射角的增加,網(wǎng)衣的受力是逐漸減小的。入射角0°時,綱繩和網(wǎng)線受力最大;入射角80°時,綱繩和網(wǎng)線受力最小。波流復合工況下入射角0°時,圍欄網(wǎng)衣受力比較均勻,受力成對稱分布,隨著角度轉(zhuǎn)變圍欄網(wǎng)衣右側(cè)區(qū)域(圖中紅色虛線矩形框)網(wǎng)衣的綱繩和網(wǎng)線受力要明顯大于左側(cè)(圖中黑色虛線矩形框)。為了對比左右兩端網(wǎng)衣受力隨入射角變化情況,研究取由上而下第二根和第三根左右兩端橫向綱繩受力進行對比。研究結(jié)果表明波流復合工況入射角度0°時,網(wǎng)衣左右兩邊綱繩受力相等。隨著入射角度的轉(zhuǎn)變,圍欄網(wǎng)衣受力向右側(cè)傾斜。波流復合工況入射角度10°時,右側(cè)綱繩受力要比左側(cè)綱繩受力大27.04%;波流復合工況入射角度40°時,右側(cè)綱繩受力比左側(cè)綱繩受力大42.44%;波流復合工況入射角度60°時,右側(cè)綱繩受力比左側(cè)綱繩受力大86.01%;波流復合工況下入射角度80°時,右側(cè)綱繩受力比左側(cè)綱繩受力大92.15%。

    注:紅色實線矩形框為綱繩最大受力處,黑色實線橢圓框為網(wǎng)線最大受力處,紅色虛線框為右側(cè)網(wǎng)衣區(qū)域,黑色虛線框為左側(cè)網(wǎng)衣區(qū)域。

    在工程建造和實際運用中,一方面,在網(wǎng)衣設(shè)計制造時,網(wǎng)衣受力最大處需要采用特殊材料或者加固處理;另一方面,在復合惡劣工況結(jié)束后,最大受力處的綱繩和網(wǎng)線需要進行重點檢查,以防養(yǎng)殖魚類的逃逸。圖7為圍欄網(wǎng)衣在不同波高時的綱繩和網(wǎng)線最大受力隨入射角的變化,從最大受力數(shù)值上看,綱繩最大受力在入射角0°~40°時,斜率變化是平穩(wěn)的,受力逐漸減小。但在入射角40°~80°綱繩受力減小速率加快。網(wǎng)線最大受力在入射角10°~60°時,網(wǎng)線最大受力下降速率較快??傮w上看,波流復合工況下正向拍打網(wǎng)衣,網(wǎng)衣承受受力最大。隨著入射角的增大,波流共同作用的于網(wǎng)衣的水質(zhì)點速度在降低,綱繩和網(wǎng)線最大受力減小,角度的變化本質(zhì)上共同作用的水質(zhì)點速度在變化。綜合綱繩和網(wǎng)線受力下降速率,網(wǎng)衣的綱繩和網(wǎng)線受力在入射角40°~60°變化時,受力下降速率較快。

    注:實線為綱繩受力,虛線為網(wǎng)線受力。

    2.2 入射角度對圍欄網(wǎng)衣最大偏移量的影響

    圖8為波流復合工況下,圍欄網(wǎng)衣在不同入射角下的最大偏移圖。從圖中可以看出隨著入射角的增加,網(wǎng)衣運動變形的形狀從大的半橢圓慢慢變小,沿方向運動最大偏移量的值會逐漸減小。最大偏移的位置是跟隨入射角的位置變化。圖9給出網(wǎng)衣最大偏移數(shù)值隨入射角的變化情況。從圖中可以看出入射角在0°~40°時,圍欄網(wǎng)衣最大偏移是等比例減小的。入射角在40°~80°時,圍欄網(wǎng)衣最大偏移量減小速率加快,也間接反映網(wǎng)衣最大受力值較小。入射角度相同的情況下,隨著波高的增加圍欄網(wǎng)衣的最大偏移量也越大。入射角的改變對網(wǎng)衣偏移的影響很大,同等波高條件下,入射角從0°轉(zhuǎn)到40°、80°時最大偏移量分別降低了16.01%、58.40%。從偏移量變化的數(shù)值上看,網(wǎng)衣在入射角40°~80°時,對網(wǎng)衣最大偏移量影響較大。

    圖8 波流不同入射角下的圍欄網(wǎng)衣最大偏移量圖

    圖9 不同波高時圍欄網(wǎng)衣最大偏移值隨入射角的變化

    2.3 入射角度對圍欄網(wǎng)衣網(wǎng)線受力的影響

    研究結(jié)果表明由上而下前3根橫向綱繩之間區(qū)域的網(wǎng)線受力要遠遠大于下而上后3根橫向綱繩之間區(qū)域的網(wǎng)線受力,為了進一步研究該區(qū)域的網(wǎng)線受力,將前3根橫向綱繩和4根豎向綱繩將網(wǎng)衣網(wǎng)線劃分成6個區(qū)域,將網(wǎng)線的目腳按順時針編號目腳1~4,如圖10所示。結(jié)合圖5網(wǎng)衣受力分布,波流入射角0°時,A1、A3和A4、A6區(qū)域網(wǎng)線受力沿中軸線對稱。研究結(jié)果表明A1、A4區(qū)域中,目腳2和3的受力是遠遠大于目腳1和4;在A3、A6區(qū)域中,目腳1和4的受力是遠遠大于目腳2和3,其正好和A1、A3區(qū)域相反。中間部位A2、A5區(qū)域,目腳1、2、3、4的受力幾乎相等。當波流復合工況入射角度10°時,網(wǎng)衣網(wǎng)線受力幾乎和入射角0°相似。但入射角度40°、60°、80°時,A3、A6區(qū)域受力要大于其他區(qū)域網(wǎng)線受力。并且在這兩個區(qū)域網(wǎng)線的目腳1和4承受較大受力。

    注:A1~A6為網(wǎng)衣劃分的6個區(qū)域;1~4為目腳編號。

    2.4 入射角度對圍欄網(wǎng)衣系縛點受力的影響

    波流復合工況下,波高的改變,只會改變系縛點的受力數(shù)值大小。入射角度的改變會改變系縛點的受力部位。在入射角0°時,圍欄網(wǎng)衣系縛點受力成對稱分布,網(wǎng)衣承受的波浪力主要由系縛點P1~P6承擔(圖11),其中系縛點P2、P5的受力最大。當入射角度的增大,系縛點P4、P5受力較大。以P2、P5系縛點為例,在入射角0°時,左右兩邊系縛點受力大小相同;在入射角10°時,P5受力比P2受力大22.62%;在入射角40°時,P5受力比P2受力大36.43%;在入射角60°和80°時,P5受力比P2受力分別大76.28%和97.99%。此外,也進一步說明入射角從40°到60°變化后,網(wǎng)衣受力下降速率快。

    研究結(jié)果表明系縛點受力要大于綱繩和網(wǎng)線受力,需要在建設(shè)過程中仔細觀察系縛點是否固定好;在復合惡劣工況結(jié)束后需對各個系縛點進行重點檢查和加固,需重點關(guān)注系縛點P1、P2、P4、P5。

    a. 圍欄網(wǎng)衣系縛點a. The joint points of net of net enclosureb. 系縛點最大受力b. The maximum force magnitude of joint points

    3 討 論

    實際海況往往是波浪、水流以及不確定其入射方向的復雜環(huán)境,浪大流急是圍欄網(wǎng)衣所要面臨的重要環(huán)境荷載。與以往文章研究不同的是:以往文章僅考慮單因素(波高或流速)工況對圍欄網(wǎng)衣的影響,該研究針對波流復合(波高+流速)工況及其入射角對網(wǎng)衣展開系統(tǒng)的研究。

    對比波流復合工況和單因素(流速)工況[21],以流速1.0 m/s,入射角度0°為例。單因素流速與波流載荷變化顯著不同的是:單因素流速作用下根據(jù)式(4)u為定值、=0,僅存在速度力在各個深度也相同;在靜水面下綱繩和網(wǎng)線都呈現(xiàn)出相同的力學特性,水流作用于網(wǎng)衣具有整體性的特點,不會隨著水深的增加而能量降低。1)從綱繩和網(wǎng)線的最大受力角度看:首先從最大受力的數(shù)值來看,波流復合工況下綱繩和網(wǎng)線最大受力的數(shù)值比單因素工況綱繩和網(wǎng)線受力分別大85.42%和41.30%左右;其次綱繩和網(wǎng)線的最大受力部位發(fā)生改變,單因素工況下綱繩和網(wǎng)線最大部位分別出現(xiàn)在上下緣綱和網(wǎng)片上下端兩側(cè)位置,由于受力均勻,網(wǎng)衣的綱繩還是網(wǎng)線受力幾乎對稱;波流復合工況情況與上同,不在贅述。2)從網(wǎng)衣網(wǎng)線最大偏移量的角度看:首先從數(shù)值來看,在最大位移上波流復合工況位移相對于單因素工況位移幾乎相同;網(wǎng)衣的最大偏移量的位置也發(fā)生改變,單因素工況最大位移出現(xiàn)在網(wǎng)片偏上和偏下兩個位置。3)從系縛點受力的角度來看:首先從受力的數(shù)值來看,波流復合工況系縛點的最大受力在同等條件下大81.20%左右;且最大系縛點受力的部位也發(fā)生改變,單因素工況最大受力出現(xiàn)在上下緣綱左右兩端。

    在圍欄養(yǎng)殖工程中,圍欄網(wǎng)衣在結(jié)構(gòu)、尺寸、大小等有所區(qū)別,但是該研究針對縫邊式圍欄網(wǎng)衣力學特性進行了深入研究,可以為圍欄網(wǎng)衣設(shè)計、加固、以及后期維護等提供可靠依據(jù)。下一步將研究圍欄網(wǎng)衣破壞機理與設(shè)計理論的研究,以期更加科學的指導生產(chǎn)實踐。

    4 結(jié) 論

    該研究針對縫邊式圍欄養(yǎng)殖設(shè)施中的網(wǎng)衣系統(tǒng)在波流復合工況下的水動力特性進行了分析,獲得圍欄網(wǎng)衣在波流復合工況下各力學及偏移參數(shù),得到以下結(jié)論:

    1)圍欄網(wǎng)衣受力是隨波流復合工況入射角度的增加逐漸減小的,入射角度0°時,網(wǎng)衣受力最大,隨入射角度的增大受力越來越小。除入射角度0°外,圍欄網(wǎng)衣受力向右側(cè)傾斜,入射角度10°時,右側(cè)網(wǎng)衣比左側(cè)大27.04%;入射角度80°時,右側(cè)網(wǎng)衣比左側(cè)大92.15%。在實際工程運用中,在圍欄網(wǎng)衣建設(shè)過程中需對綱繩和網(wǎng)線最大受力處進行加固或采取特殊材料。

    2)在圍欄建設(shè)過程中,輔助設(shè)備的安裝應(yīng)避開網(wǎng)衣的運動,以防發(fā)生摩擦破壞。在波流復合工況入射方向0°~40°時,圍欄網(wǎng)衣最大偏移量下降速率相同,入射方向40°~80°時圍欄網(wǎng)衣偏移量下降速率增大。

    3)波流復合工況下,上半部分綱繩(由上而下前三根橫向綱繩)之間的網(wǎng)線受力要大于其他綱繩之間區(qū)域網(wǎng)線受力且左右兩側(cè)區(qū)域網(wǎng)線受力最大。

    4)研究結(jié)果表明波流復合工況下,系縛點最大受力從開始的入射角0°時兩側(cè)受力相同,隨著入射角度的增大,右側(cè)系縛點受力要大于左側(cè)。入射角在40°~80°變化時,左側(cè)系縛點受力顯著小于右側(cè)系縛點受力。在實際工程建設(shè)過程中,需要仔細觀察系縛點是否固定好。

    [1] 宋煒,韓昕辰,謝正麗,等. 我國深遠海大型圍欄養(yǎng)殖發(fā)展現(xiàn)狀與展望[J]. 漁業(yè)科學進展,2022,43(6):111-120.

    Song Wei, Han Xinchen, Xie Zhengli, et al. Development status and the prospect of deep-sea large-scale fence culture in China[J]. Progress in Fishery Sciences, 2022, 43(6): 111-120. (in Chinese with English abstract)

    [2] 江國強. 大黃魚淺海圍欄養(yǎng)殖技術(shù)研究[J]. 現(xiàn)代漁業(yè)信息,2011,26(1):3-4.

    Jiang Guoqiang. A study on enclosure culture techniques for large yellow croaker,() in epeiric sea[J]. Modern Fisheries Information, 2011, 26(1): 3-4. (in Chinese with English abstract)

    [3] 中華人民共和國住房和城鄉(xiāng)建設(shè)部. 建筑樁基技術(shù)規(guī)范:JGJ 94-2008[S].北京:中國建筑工業(yè)出版社,2008.

    [4] 中華人民共和國交通部. 港口工程樁基規(guī)范:JTS 167-4-2012[S].北京:中國標準出版社,2012.

    [5] Fredriksson D W, Swift M, Irish J D, et al. Fish cage and mooring system dynamics using physical and numerical models with field measurements[J]. Aquacultural Engineering, 2003, 27(2): 117-146.

    [6] Lee C W, Kim Y B, Lee G H, et al. Dynamic simulation of a fish cage system subjected to currents and waves[J]. Ocean Engineering, 2008, 35(14/15): 1521-1532.

    [7] Kim T H, Fredriksson D W, DeCew J. Mooring tension and motion characteristics of a submerged fish reef with net in waves and currents using numerical analysis[J]. China Ocean Engineering, 2012, 26(2): 261-276.

    [8] 黃小華,王芳芳,劉海陽,等. 系泊和壓載方式對半潛式漁場平臺動力特性的影響[J]. 農(nóng)業(yè)工程學報,2019,35(15):48-53.

    Huang Xiaohua, Wang Fangfang, Liu Haiyang, et al. Effects of mooring systems and ballast status on dynamic behaviors of semi-submersible offshore fish farm affected by wave[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(15): 48-53. (in Chinese with English abstract)

    [9] 劉海陽,胡昱,黃小華,等. 深水網(wǎng)箱浮架結(jié)構(gòu)的失效及疲勞性能分析[J]. 農(nóng)業(yè)工程學報,2020,36(3):46-54.

    Liu Haiyang, Hu Yu, Huang Xiaohua, et al. Failure and fatigue analysis of offshore structures of offshore cage[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(3): 46-54. (in Chinese with English abstract)

    [10] Song W, Liang Z, Chi H, et al. Experimental study on the effect of horizontal waves on netting panels[J]. Fisheries Science, 2006, 72(5): 967-976.

    [11] Bi C W, Zhao Y P, Dong G H, et al. Experimental investigation of the reduction in flow velocity downstream from a fishing net[J]. Aquacultural Engineering, 2013, 57: 71-81.

    [12] Tang M F. Dong G H. Xu T J, et al. Experimental analysis of the hydrodynamic coefficients of net panels in current[J]. Applied Ocean Research, 2018, 79: 253-261.

    [13] Xu Z, Qin H. Fluid-structure interactions of cage based aquaculture: from structures to organisms[J]. Ocean Engineering, 2020, 217: 107961.

    [14] 黃小華,劉海陽,胡昱,等. 深水養(yǎng)殖網(wǎng)箱浮架變形模擬及結(jié)構(gòu)改進設(shè)計[J]. 農(nóng)業(yè)工程學報,2018,34(15):44-49.

    Huang Xiaohua, Liu Haiyang, Hu Yu, et al. Deformation simulation and structural improvement design for floating collar of deep-water aquaculture net cage[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(15): 44-49. (in Chinese with English abstract)

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

    Ye Weifu. Preliminary Study on the Characteristics of Hydrodynamic on the Aquaculture Seine[D]. Zhoushan: Zhejiang Ocean University, 2012. (in Chinese with English abstract)

    [16] 崔勇,關(guān)長濤,黃濱,等. 浮繩式圍網(wǎng)水動力特性研究[J]. 漁業(yè)現(xiàn)代化,2018,45(5):14-18.

    Cui Yong, Guan Changtao, Huang Bin, et al. Study on hydrodynamic characteristics of float-type net enclosure [J]. Fisheries Modernization, 2018, 45(5):14-18. (in Chinese with English abstract)

    [17] 徐克品,詹建明,胡利永,等. 大型浮繩式圍網(wǎng)養(yǎng)殖設(shè)施的網(wǎng)衣結(jié)構(gòu)受力分析[J]. 海洋通報,2017,36(1):114-120.

    Xu Kepin, Zhan Jianming, Hu Liyong, et al. Net structure force analysis of the large-scale floating rope round net[J]. Ocean Bulletin, 2017, 36(1): 114-120. (in Chinese with English abstract)

    [18] 桂福坤,陳天華,趙云鵬,等. 固定方式對樁柱式圍網(wǎng)網(wǎng)片波浪力學特性影響研究[J]. 大連理工大學學報,2017,57(3):285-292.

    Gui Fukun, Chen Tianhua, Zhao Yunpeng, et al. Study of effect on wave mechanical properties for net panel of pile-column type net enclosure by fixations[J]. Journal of Dalian University of Technology, 2017, 57(3): 285-292. (in Chinese with English abstract)

    [19] 桂福坤,邵振宇,陳天華,等. 綱繩對樁柱式圍網(wǎng)網(wǎng)片波浪力學特性影響研究[J]. 水動力學研究與進展(A輯),2020,35(4):505-514.

    Gui Fukun, Shao Zhenyu, Chen Tianhua, et al. Study on effect of wave mechanical properties for net panel of pile-column type net enclosure[J]. Research and Progress in Hydrodynamics (Series A), 2020, 35(4): 505-514. (in Chinese with English abstract)

    [20] 陳天華,孟昂,桂福坤. 波浪高度及方向?qū)吨絿W(wǎng)養(yǎng)殖系統(tǒng)網(wǎng)片水力特性的影響[J]. 農(nóng)業(yè)工程學報,2017,33(2):245-251.

    Chen Tianhua, Meng Ang, Gui Fukun. Effect of wave height and direction on the hydraulic characteristics of the pile-column type net enclosure aquaculture system[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(2): 245-251. (in Chinese with English abstract)

    [21] 陳天華,潘昀,孟昂,等. 樁柱式圍欄單元網(wǎng)片在水流作用下的水動力特性研究[J]. 水動力學研究與進展(A輯),2017,32(4):511-519.

    Chen Tianhua, Pan Yun, Meng Ang, et al. Study on the hydrodynamic of unit net panel of the column-type net enclosure in currents[J]. Hydrodynamic Research and Progress (Series A), 2017, 32(4): 511-519. (in Chinese with English abstract)

    [22] 楊熙,劉燦,潘昀,等. 網(wǎng)片出水條件下波高和水深對其水動力特性的影響[J]. 浙江海洋大學學報(自然科學版),2020,39(4):347-352.

    Yang Xi, Liu Can, Pan Yun, et al. Effects of wave heights and water depths on hydraulic characteristics of net under the water surface conditions[J]. Journal of Zhejiang Ocean University (Natural Science Edition), 2020, 39(4): 347-352. (in Chinese with English abstract)

    [23] 陳天華. 樁柱式圍網(wǎng)養(yǎng)殖系統(tǒng)水動力特性研究[D]. 舟山:浙江海洋大學,2017.

    Chen Tianhua. Study on Hydrodynamic Characteristics of Pile-column Type Net Enclosure Aquaculture Facility[D]. Zhoushan: Zhejiang Ocean University, 2017. (in Chinese with English abstract)

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

    Zhao Yunpeng. Numerical Investigation on Hydrodynamic Behavior of Deep-water Gravity Cage[D]. Dalian: Dalian University of Technology, 2007. (in Chinese with English abstract)

    [25] 王樹青,梁丙臣. 海洋工程波浪力學[M]. 青島:中國海洋大學出版社,2013:130-174.

    [26] 萬榮,何鑫,王欣欣,等. 一種適用于網(wǎng)箱耐流特性有限元分析的網(wǎng)目群化方法[J]. 中國海洋大學學報(自然科學版),2007(6):885-888.

    Wan Rong, He Xin, Wang Xinxin, et al. A method of mesh group for fish cage anti-current characteristics by finite element analysis[J]. Journal of Ocean University of China (Natural Science Edition), 2007(6): 885-888. (in Chinese with English abstract)

    [27] Li Y C, Zhao Y P, Gui F K, et al. Numerical simulation of the hydrodynamic behavior of submerged plane nets in current[J]. Ocean Engineering, 2006, 33: 2352-2368.

    [28] 李玉成,桂福坤. 平面有結(jié)節(jié)和無結(jié)節(jié)網(wǎng)目試驗及水阻力系數(shù)的選擇[J]. 中國海洋平臺,2005(6):11-17.

    Li Yucheng, Gui Fukun. Expermental study and selection of the drag coefficient of knotted and knot less plane fishing net[J]. China Ocean Platform, 2005(6): 11-17. (in Chinese with English abstract)

    [29] 桂福坤. 深水重力式網(wǎng)箱水動力學特性研究[D]. 大連:大連理工大學,2006.

    Gui Fukui. Hydrodynamic Behavior of Deep-water Gravity Cage[D]. Dalian: Dalian University of Technology, 2006. (in Chinese with English abstract)

    [30] 王順,馮德軍,桂福坤. 樁網(wǎng)連接方式對樁柱式圍網(wǎng)網(wǎng)衣受力的影響[J]. 漁業(yè)科學進展,2022,43(6):29-39.

    Wang Shun, Feng Dejun, Gui Fukun.Study on effect of the connection mode between pile and net on mechanical properties for mesh of pile-column type net enclosure[J].Advances in Fisheries Science, 2022, 43(6): 29-39. (in Chinese with English abstract)

    [31] 郭帥. 深水大型養(yǎng)殖網(wǎng)箱抗臺設(shè)計及水動力性能分析[D]. 鎮(zhèn)江:江蘇科技大學,2020.

    Guo Shuai. Anti-typhoon Design and Hydrodynamic Performance Analysis of the Large Deep-Water Fish Cage[D]. Zhenjiang: Jiangsu University of Science and Technology, 2020. (in Chinese with English abstract)

    Hydrodynamic characteristics of the net of seam-edge type net enclosure under the combined wave-current loads

    Wang Shun, Feng Dejun, Gui Fukun※

    (,,316022,)

    Net enclosure aquacultureis one of the most important ways for farming fish and economic animals.Among them, nets, piles, ropes and iron chains can be used to delineate a certain water column in the shallow seas. More than 10 net enclosure facilities have been built across the countryfrom the Shandong to Guangdong Province in China in recent years.However, only a few studies have been focused on the hydrodynamic performance of net enclosure aquaculture and emerging marine aquaculture facilities. The net enclosure facilities are still in the early stages. It is also lacking under the extreme ocean load resistance of aquaculture engineering structures. In addition, there is no such structural design theory or standard specifications for the net enclosure design and installation, particularly on the safety assessment. In this study, a net enclosure aquaculture facility was numerically investigated using the lumped mass model. A Newton’s second-law-based motion equation was then solved using Euler’s method. Finally, the MATLAB platform was used to visualize the calculations.The results showed that the maximum force of ropes was found in the second horizontal rope and the third horizontal rope under the still water level. Furthermore, the maximum forces of ropes and net twine were get the wave-current directions of 0°, 10°, 40°, 60°, and 80°, respectively, when the wave and current were 4 m, and 1.0 m/s, respectively. The force on the netting was tilted to the right (incident direction), as the incident angle was changed gradually. As such, the force on the right side of the rope was much greater than that left side. The maximum force on the right side of the horizontal rope was 27.04% higher than that on the left side at the wave-current incident direction of 10°. Moreover, the maximum forces on the right side were 42.44%, 86.01% and 92.15% higher than that on the left side in the wave-current incident direction of 40°, 60? and 80°, respectively. The maximum displacementof the net decreased greatly, as the incident angles increased and the maximum positions of the net moved to the left. The force of the net significantly increased with the waves. The load of the entire net was mainly from the top half of the net. The main force of the net twines was distributed in the six areas that were divided by the first three horizontal ropes and the two vertical ropes. The maximum force on the net twines occurs at the position on both sides above the second horizontal rope under the water surface. The force of net twines on the right side of the rope was greater than that on the left side, with the change of the incident direction. The maximum force decreased rapidly at the incident angle of 40°-80° under the combined wave-current conditions.The maximum force of joint points appeared at positions 1, 2, 3, and 4. The numerical simulation was then performed on the net of the seam-edge type net enclosure aquaculture facility in use. The force distribution and deformation were achieved in the seam-edge type net enclosure. The finding can provide the theoretical basis for the design, installation, and reinforcement of aquaculture facilities. In addition, it was necessary to focus on the inspection and reinforcement after the event of extremely bad weather.

    aquaculture; models; mechanics properties; net; wave-current conditions

    10.11975/j.issn.1002-6819.2022.19.007

    S954.1

    A

    1002-6819(2022)-19-0060-09

    王順,馮德軍,桂福坤. 波流復合工況下縫邊式圍欄網(wǎng)衣水動力特性[J]. 農(nóng)業(yè)工程學報,2022,38(19):60-68.doi:10.11975/j.issn.1002-6819.2022.19.007 http://www.tcsae.org

    Wang Shun, Feng Dejun, Gui Fukun. Hydrodynamic characteristics of the net of seam-edge type net enclosure under the combined wave-current loads[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2022, 38(19): 60-68. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2022.19.007 http://www.tcsae.org

    2022-06-17

    2022-09-29

    國家重點研發(fā)計劃項目(2019YFD0900902);國家自然科學基金項目(32002441)

    王順,研究方向為海洋養(yǎng)殖工程裝備技術(shù)。Email:1622426519@qq.com

    桂福坤,博士,教授,博士生導師,研究方向為海洋養(yǎng)殖工程裝備技術(shù)。Email:gui2237@163.com.

    猜你喜歡
    網(wǎng)衣波流網(wǎng)線
    遮蔽效應(yīng)下雙平面網(wǎng)衣和圓形網(wǎng)衣的水動力性能
    基于CFD的剛性養(yǎng)殖網(wǎng)衣流場數(shù)值模擬及不確定度分析
    波流耦合下樁周珊瑚砂沖刷機理研究
    海洋通報(2022年5期)2022-11-30 12:04:02
    波流聯(lián)合作用下海上輸油漂浮軟管動力響應(yīng)分析
    偶感
    關(guān)于如何降低配網(wǎng)線損的研究
    電子測試(2018年14期)2018-09-26 06:04:42
    槽道內(nèi)渦波流場展向渦的分布特征
    遼墓出土網(wǎng)衣編法三例
    遼墓出土銅絲網(wǎng)衣修復與復原報告
    加強配網(wǎng)線損管理 提高企業(yè)經(jīng)濟效益
    亚洲av成人精品一二三区| 999久久久国产精品视频| 日韩 欧美 亚洲 中文字幕| 国产在线视频一区二区| 亚洲美女黄色视频免费看| 日本猛色少妇xxxxx猛交久久| 好男人视频免费观看在线| 亚洲五月色婷婷综合| videos熟女内射| 国产欧美亚洲国产| 国产成人av教育| 啦啦啦中文免费视频观看日本| 久久久久国产精品人妻一区二区| 午夜免费男女啪啪视频观看| 免费看十八禁软件| 男女国产视频网站| 一级黄片播放器| 一区二区日韩欧美中文字幕| 国产黄色免费在线视频| 两个人看的免费小视频| 另类精品久久| 19禁男女啪啪无遮挡网站| 成年人黄色毛片网站| 国产高清国产精品国产三级| 麻豆乱淫一区二区| 午夜福利,免费看| 一级,二级,三级黄色视频| 亚洲综合色网址| 美国免费a级毛片| 老熟女久久久| 久久中文字幕一级| 午夜激情久久久久久久| 狂野欧美激情性bbbbbb| 久久久欧美国产精品| 久久99一区二区三区| 中文字幕色久视频| 成年人午夜在线观看视频| 国产一区有黄有色的免费视频| 国产高清不卡午夜福利| 菩萨蛮人人尽说江南好唐韦庄| 免费一级毛片在线播放高清视频 | 欧美大码av| 下体分泌物呈黄色| 一区二区日韩欧美中文字幕| 国产黄频视频在线观看| 日日夜夜操网爽| 韩国高清视频一区二区三区| 国产不卡av网站在线观看| 日韩一本色道免费dvd| 高清av免费在线| 国产精品熟女久久久久浪| 亚洲少妇的诱惑av| 免费观看av网站的网址| 午夜久久久在线观看| 久久久久网色| 免费看不卡的av| 亚洲一区二区三区欧美精品| 黄色a级毛片大全视频| 激情视频va一区二区三区| 国产伦理片在线播放av一区| 麻豆国产av国片精品| 色视频在线一区二区三区| 亚洲中文字幕日韩| 国产三级黄色录像| 99热网站在线观看| 丝袜脚勾引网站| h视频一区二区三区| 大型av网站在线播放| 国产亚洲欧美精品永久| 麻豆乱淫一区二区| 性高湖久久久久久久久免费观看| 亚洲成国产人片在线观看| 亚洲自偷自拍图片 自拍| 少妇精品久久久久久久| 天天躁夜夜躁狠狠久久av| 国产精品欧美亚洲77777| 97人妻天天添夜夜摸| 人妻 亚洲 视频| 亚洲一码二码三码区别大吗| 青春草亚洲视频在线观看| 一边摸一边做爽爽视频免费| 99re6热这里在线精品视频| 后天国语完整版免费观看| 久久久久久久精品精品| av在线app专区| 中文字幕人妻丝袜制服| 久久ye,这里只有精品| 国产亚洲一区二区精品| 五月开心婷婷网| 极品人妻少妇av视频| av电影中文网址| tube8黄色片| 99国产精品一区二区蜜桃av | 人人妻人人爽人人添夜夜欢视频| 老司机影院毛片| 中文字幕另类日韩欧美亚洲嫩草| 久久人妻熟女aⅴ| 视频区欧美日本亚洲| 亚洲欧美精品自产自拍| 青草久久国产| 18禁裸乳无遮挡动漫免费视频| xxxhd国产人妻xxx| 精品福利观看| 精品人妻1区二区| 中文字幕人妻丝袜一区二区| 亚洲国产精品一区二区三区在线| 欧美日韩福利视频一区二区| 成人影院久久| 国产精品熟女久久久久浪| 国产爽快片一区二区三区| a级片在线免费高清观看视频| 日本a在线网址| 99热网站在线观看| 在线观看免费高清a一片| 丝袜脚勾引网站| 国产不卡av网站在线观看| 午夜av观看不卡| 黄色毛片三级朝国网站| 亚洲国产看品久久| 日本色播在线视频| 不卡av一区二区三区| 成人午夜精彩视频在线观看| 日本色播在线视频| 亚洲自偷自拍图片 自拍| 一区福利在线观看| 婷婷色综合大香蕉| 亚洲精品日本国产第一区| 国产欧美日韩一区二区三 | 国产一区二区在线观看av| 91九色精品人成在线观看| 欧美在线黄色| 午夜av观看不卡| 亚洲av片天天在线观看| 少妇被粗大的猛进出69影院| bbb黄色大片| 最近手机中文字幕大全| 永久免费av网站大全| svipshipincom国产片| 性色av一级| av有码第一页| 亚洲精品国产一区二区精华液| 国产精品一二三区在线看| 久久鲁丝午夜福利片| 亚洲欧洲国产日韩| 久久精品亚洲熟妇少妇任你| 看免费成人av毛片| 大片电影免费在线观看免费| 性色av一级| 人人妻人人澡人人爽人人夜夜| 一本色道久久久久久精品综合| 丝袜脚勾引网站| 热99久久久久精品小说推荐| 色婷婷av一区二区三区视频| 一区二区三区激情视频| 一级毛片黄色毛片免费观看视频| 50天的宝宝边吃奶边哭怎么回事| 国产成人系列免费观看| 99精品久久久久人妻精品| 天天躁日日躁夜夜躁夜夜| 久久精品aⅴ一区二区三区四区| 制服诱惑二区| 老司机影院成人| 宅男免费午夜| 菩萨蛮人人尽说江南好唐韦庄| a级毛片黄视频| 黄片小视频在线播放| 日本wwww免费看| 女性被躁到高潮视频| 亚洲成人国产一区在线观看 | tube8黄色片| 亚洲一区二区三区欧美精品| 一级毛片我不卡| 两个人免费观看高清视频| 午夜av观看不卡| 无限看片的www在线观看| 欧美变态另类bdsm刘玥| 亚洲免费av在线视频| 欧美精品亚洲一区二区| 十八禁网站网址无遮挡| 熟女少妇亚洲综合色aaa.| 亚洲av综合色区一区| 国产av精品麻豆| 黑人猛操日本美女一级片| 午夜精品国产一区二区电影| 日韩av不卡免费在线播放| 久久 成人 亚洲| 美女午夜性视频免费| 制服诱惑二区| 亚洲欧美一区二区三区国产| 黄色视频不卡| 久久午夜综合久久蜜桃| 精品一区二区三卡| 性少妇av在线| 一级黄片播放器| 亚洲人成77777在线视频| 国产成人啪精品午夜网站| 黄色 视频免费看| 天堂中文最新版在线下载| 天堂8中文在线网| 91成人精品电影| 中国美女看黄片| 真人做人爱边吃奶动态| 人妻一区二区av| 性少妇av在线| 久久久久久人人人人人| 麻豆av在线久日| 国产成人a∨麻豆精品| 欧美黑人精品巨大| 18禁裸乳无遮挡动漫免费视频| 建设人人有责人人尽责人人享有的| 久久精品人人爽人人爽视色| 男女边吃奶边做爰视频| 欧美av亚洲av综合av国产av| 黄色片一级片一级黄色片| a 毛片基地| 五月开心婷婷网| 欧美另类一区| 国产精品久久久av美女十八| 操出白浆在线播放| 韩国高清视频一区二区三区| 母亲3免费完整高清在线观看| 啦啦啦在线免费观看视频4| 男人添女人高潮全过程视频| 少妇人妻久久综合中文| 国产无遮挡羞羞视频在线观看| 你懂的网址亚洲精品在线观看| 欧美日韩视频高清一区二区三区二| 涩涩av久久男人的天堂| 日本一区二区免费在线视频| 亚洲熟女毛片儿| 黄频高清免费视频| 男男h啪啪无遮挡| 国产又色又爽无遮挡免| 亚洲综合色网址| 久久久久久久国产电影| 亚洲精品久久久久久婷婷小说| 丰满少妇做爰视频| 国产一区二区激情短视频 | 一级黄色大片毛片| 免费观看av网站的网址| 国产免费一区二区三区四区乱码| 成年人免费黄色播放视频| 2018国产大陆天天弄谢| tube8黄色片| 别揉我奶头~嗯~啊~动态视频 | 国产午夜精品一二区理论片| 国产精品人妻久久久影院| 欧美日韩亚洲高清精品| 国产成人影院久久av| 免费在线观看视频国产中文字幕亚洲 | 在线观看免费视频网站a站| 黄色片一级片一级黄色片| 国产视频一区二区在线看| 一级毛片黄色毛片免费观看视频| 男男h啪啪无遮挡| a级片在线免费高清观看视频| 精品久久蜜臀av无| 精品人妻一区二区三区麻豆| 黄色片一级片一级黄色片| 国产精品国产三级国产专区5o| 国产亚洲精品久久久久5区| 国产一区亚洲一区在线观看| 国产精品熟女久久久久浪| 两个人免费观看高清视频| 在线观看免费日韩欧美大片| 亚洲精品国产色婷婷电影| 精品亚洲乱码少妇综合久久| 精品国产乱码久久久久久小说| 久久久久国产精品人妻一区二区| 纵有疾风起免费观看全集完整版| av不卡在线播放| 成人亚洲精品一区在线观看| 久久99热这里只频精品6学生| 久热这里只有精品99| 在线av久久热| 日本五十路高清| 如日韩欧美国产精品一区二区三区| 国产老妇伦熟女老妇高清| 国产高清不卡午夜福利| 国产av国产精品国产| 国产又色又爽无遮挡免| 国产精品欧美亚洲77777| 波多野结衣av一区二区av| 国产极品粉嫩免费观看在线| 久久久久久免费高清国产稀缺| 老鸭窝网址在线观看| 好男人视频免费观看在线| 精品一品国产午夜福利视频| 十八禁人妻一区二区| 青青草视频在线视频观看| 久久精品久久久久久久性| 欧美激情 高清一区二区三区| 久9热在线精品视频| 午夜日韩欧美国产| 国产熟女欧美一区二区| 七月丁香在线播放| 手机成人av网站| 免费看不卡的av| 欧美日韩福利视频一区二区| 一本大道久久a久久精品| 精品人妻在线不人妻| 久久99精品国语久久久| 好男人视频免费观看在线| 亚洲精品久久成人aⅴ小说| 波多野结衣一区麻豆| 一级毛片我不卡| 欧美精品高潮呻吟av久久| 日韩av在线免费看完整版不卡| 亚洲精品国产av成人精品| 久久久国产精品麻豆| 黄色a级毛片大全视频| 久久毛片免费看一区二区三区| xxx大片免费视频| 精品人妻1区二区| 日本一区二区免费在线视频| 亚洲 国产 在线| 免费观看a级毛片全部| 日韩中文字幕欧美一区二区 | 一边摸一边做爽爽视频免费| 汤姆久久久久久久影院中文字幕| 一级黄片播放器| 一级片免费观看大全| 国产精品久久久久成人av| 女性生殖器流出的白浆| 91精品伊人久久大香线蕉| 日本黄色日本黄色录像| 脱女人内裤的视频| 黄色视频不卡| 在线观看免费高清a一片| 亚洲精品美女久久av网站| 亚洲av电影在线观看一区二区三区| 国产xxxxx性猛交| 国产成人免费观看mmmm| 无遮挡黄片免费观看| 久久这里只有精品19| 狠狠精品人妻久久久久久综合| 日本av免费视频播放| 国产精品一二三区在线看| 丝袜人妻中文字幕| 亚洲精品在线美女| 欧美日韩视频精品一区| 999精品在线视频| 国产精品国产三级专区第一集| 亚洲欧美精品综合一区二区三区| 欧美xxⅹ黑人| 日日夜夜操网爽| 成人午夜精彩视频在线观看| 又大又黄又爽视频免费| 超碰成人久久| 亚洲国产看品久久| 亚洲一卡2卡3卡4卡5卡精品中文| 国语对白做爰xxxⅹ性视频网站| 91麻豆精品激情在线观看国产 | 亚洲av日韩在线播放| 校园人妻丝袜中文字幕| 丝袜喷水一区| 国产一区二区三区av在线| 青春草亚洲视频在线观看| 婷婷色av中文字幕| 免费观看av网站的网址| 亚洲成国产人片在线观看| 侵犯人妻中文字幕一二三四区| 成人亚洲欧美一区二区av| 欧美另类一区| 新久久久久国产一级毛片| 日本wwww免费看| 丝瓜视频免费看黄片| 免费在线观看黄色视频的| 午夜精品国产一区二区电影| 极品人妻少妇av视频| 精品亚洲乱码少妇综合久久| 极品人妻少妇av视频| 亚洲精品一区蜜桃| 可以免费在线观看a视频的电影网站| 波多野结衣av一区二区av| 18禁裸乳无遮挡动漫免费视频| 午夜福利视频在线观看免费| 91麻豆av在线| 亚洲人成网站在线观看播放| 国产精品一区二区免费欧美 | 午夜av观看不卡| 丝袜在线中文字幕| 中文字幕人妻熟女乱码| 人人妻,人人澡人人爽秒播 | av线在线观看网站| 十分钟在线观看高清视频www| 极品少妇高潮喷水抽搐| 国产一区亚洲一区在线观看| 午夜免费观看性视频| 男女边吃奶边做爰视频| 久久av网站| xxxhd国产人妻xxx| 亚洲人成77777在线视频| 在线观看免费高清a一片| 亚洲中文av在线| av在线播放精品| 精品人妻熟女毛片av久久网站| 午夜免费鲁丝| 精品少妇黑人巨大在线播放| 日日摸夜夜添夜夜爱| av有码第一页| 午夜91福利影院| 精品国产一区二区三区久久久樱花| 9191精品国产免费久久| 99精国产麻豆久久婷婷| 老司机亚洲免费影院| 精品免费久久久久久久清纯 | 91麻豆精品激情在线观看国产 | 一级毛片 在线播放| 丁香六月天网| 亚洲国产欧美在线一区| 国产精品香港三级国产av潘金莲 | 亚洲成人国产一区在线观看 | 欧美激情极品国产一区二区三区| 精品少妇内射三级| 精品国产一区二区三区久久久樱花| 欧美日韩国产mv在线观看视频| 亚洲av电影在线进入| 欧美日韩国产mv在线观看视频| 国产av精品麻豆| 日韩免费高清中文字幕av| 少妇裸体淫交视频免费看高清 | 在线av久久热| 亚洲国产看品久久| 久久久精品94久久精品| 尾随美女入室| 中文字幕av电影在线播放| 久久久国产精品麻豆| 99九九在线精品视频| 99国产综合亚洲精品| 国产在视频线精品| 国产免费福利视频在线观看| 国产黄色视频一区二区在线观看| 欧美亚洲 丝袜 人妻 在线| 亚洲成人免费电影在线观看 | 午夜精品国产一区二区电影| 99久久精品国产亚洲精品| 丝袜喷水一区| 高清av免费在线| 少妇被粗大的猛进出69影院| 一级毛片女人18水好多 | 一本综合久久免费| 天堂中文最新版在线下载| 国产老妇伦熟女老妇高清| 亚洲国产欧美一区二区综合| 婷婷色av中文字幕| 国产精品秋霞免费鲁丝片| 乱人伦中国视频| 色婷婷av一区二区三区视频| 美女午夜性视频免费| 亚洲欧洲日产国产| 欧美老熟妇乱子伦牲交| 日日爽夜夜爽网站| 国产有黄有色有爽视频| 97人妻天天添夜夜摸| 欧美亚洲 丝袜 人妻 在线| 欧美亚洲日本最大视频资源| 久久久久久久精品精品| 欧美日本中文国产一区发布| 咕卡用的链子| 男女边摸边吃奶| 飞空精品影院首页| 各种免费的搞黄视频| 赤兔流量卡办理| 国产99久久九九免费精品| 叶爱在线成人免费视频播放| 水蜜桃什么品种好| 99精品久久久久人妻精品| 91成人精品电影| 国产精品欧美亚洲77777| 男女床上黄色一级片免费看| 99九九在线精品视频| av网站免费在线观看视频| 久久精品久久久久久久性| 久久久精品国产亚洲av高清涩受| 男女无遮挡免费网站观看| 国产精品一区二区免费欧美 | e午夜精品久久久久久久| 成年动漫av网址| 亚洲欧美日韩高清在线视频 | 国产在线视频一区二区| 亚洲自偷自拍图片 自拍| 婷婷色麻豆天堂久久| 999久久久国产精品视频| 精品久久蜜臀av无| 深夜精品福利| 亚洲欧美清纯卡通| 最新的欧美精品一区二区| 天天躁夜夜躁狠狠躁躁| 久久久国产一区二区| 亚洲av美国av| av国产久精品久网站免费入址| 母亲3免费完整高清在线观看| 国产成人欧美| 欧美黄色淫秽网站| 宅男免费午夜| 天天躁日日躁夜夜躁夜夜| 亚洲激情五月婷婷啪啪| 欧美+亚洲+日韩+国产| 日韩av免费高清视频| 欧美人与善性xxx| 午夜91福利影院| 日韩av在线免费看完整版不卡| 在线亚洲精品国产二区图片欧美| bbb黄色大片| 乱人伦中国视频| 在线av久久热| 日韩视频在线欧美| 高清av免费在线| 在线观看一区二区三区激情| 亚洲情色 制服丝袜| 午夜福利在线免费观看网站| 久久久久国产精品人妻一区二区| 亚洲av成人精品一二三区| 在线观看免费视频网站a站| 国产av国产精品国产| 高清黄色对白视频在线免费看| 亚洲精品国产区一区二| 精品欧美一区二区三区在线| 亚洲精品美女久久久久99蜜臀 | 亚洲一区二区三区欧美精品| 久久精品国产a三级三级三级| 国产精品 欧美亚洲| 亚洲国产日韩一区二区| 亚洲午夜精品一区,二区,三区| 69精品国产乱码久久久| 免费少妇av软件| 青青草视频在线视频观看| 中文字幕人妻丝袜制服| 妹子高潮喷水视频| 日韩中文字幕视频在线看片| 久久久久精品人妻al黑| 亚洲一区中文字幕在线| 女警被强在线播放| 午夜福利乱码中文字幕| 丝袜美足系列| 亚洲精品久久久久久婷婷小说| 97在线人人人人妻| 精品第一国产精品| 国产亚洲精品久久久久5区| 在线观看www视频免费| 爱豆传媒免费全集在线观看| 国产成人一区二区在线| 老司机亚洲免费影院| 别揉我奶头~嗯~啊~动态视频 | 国产精品 国内视频| 999久久久国产精品视频| 男女边吃奶边做爰视频| 免费高清在线观看日韩| 午夜两性在线视频| 爱豆传媒免费全集在线观看| 日本黄色日本黄色录像| 97人妻天天添夜夜摸| 国产麻豆69| 国产精品久久久人人做人人爽| 91精品国产国语对白视频| 国产精品久久久久久人妻精品电影 | 青春草亚洲视频在线观看| 你懂的网址亚洲精品在线观看| 超碰97精品在线观看| 亚洲国产欧美在线一区| 美女高潮到喷水免费观看| 亚洲专区国产一区二区| 七月丁香在线播放| 十八禁网站网址无遮挡| a 毛片基地| 极品人妻少妇av视频| 久久av网站| 女性被躁到高潮视频| 又大又黄又爽视频免费| 久久久欧美国产精品| 最近手机中文字幕大全| 国产熟女午夜一区二区三区| 精品人妻在线不人妻| 丰满少妇做爰视频| 午夜av观看不卡| 在线观看免费视频网站a站| 亚洲国产欧美网| 国产精品国产三级专区第一集| 少妇精品久久久久久久| 日日爽夜夜爽网站| 这个男人来自地球电影免费观看| 一个人免费看片子| 高清不卡的av网站| 天天操日日干夜夜撸| 国产无遮挡羞羞视频在线观看| 少妇人妻 视频| 黄色a级毛片大全视频| 亚洲视频免费观看视频| 久久99热这里只频精品6学生| 欧美激情极品国产一区二区三区| 少妇猛男粗大的猛烈进出视频| 天天添夜夜摸| 欧美性长视频在线观看| 色精品久久人妻99蜜桃| 亚洲五月色婷婷综合| 一区在线观看完整版| 丝瓜视频免费看黄片| 亚洲av男天堂| 国产一区有黄有色的免费视频| 国产精品三级大全| av欧美777| 在线观看国产h片| 国产一区二区激情短视频 | 香蕉丝袜av| 老司机亚洲免费影院| 男女边摸边吃奶| 蜜桃国产av成人99| 久久这里只有精品19| 青春草视频在线免费观看| 中文字幕最新亚洲高清| 国产欧美日韩精品亚洲av| 啦啦啦视频在线资源免费观看| 最近中文字幕2019免费版| 久久精品亚洲熟妇少妇任你|