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

    水流作用下重力式網(wǎng)箱浮架結(jié)構(gòu)的動(dòng)力響應(yīng)研究

    2016-11-25 05:25:11李坤鵬許條建董國(guó)海趙云鵬
    漁業(yè)現(xiàn)代化 2016年1期
    關(guān)鍵詞:架結(jié)構(gòu)網(wǎng)箱圓環(huán)

    李坤鵬, 許條建,2, 董國(guó)海, 趙云鵬

    (1大連理工大學(xué)海岸和近海工程國(guó)家重點(diǎn)實(shí)驗(yàn)室,遼寧 大連,116024 2大連理工大學(xué)船舶工程學(xué)院,遼寧 大連,116024

    ?

    水流作用下重力式網(wǎng)箱浮架結(jié)構(gòu)的動(dòng)力響應(yīng)研究

    李坤鵬1, 許條建1,2, 董國(guó)海1, 趙云鵬1

    (1大連理工大學(xué)海岸和近海工程國(guó)家重點(diǎn)實(shí)驗(yàn)室,遼寧 大連,116024 2大連理工大學(xué)船舶工程學(xué)院,遼寧 大連,116024

    近年來(lái)深水網(wǎng)箱養(yǎng)殖設(shè)施得到了廣泛的應(yīng)用,浮架作為網(wǎng)箱結(jié)構(gòu)的重要組成部分,其安全性對(duì)網(wǎng)箱的設(shè)計(jì)至關(guān)重要?;谟邢迒卧ú捎肧HELL單元建立了浮架結(jié)構(gòu)的荷載-變形數(shù)值模型,對(duì)水流作用下浮架的應(yīng)力和變形進(jìn)行了數(shù)值模擬,并開(kāi)展物理模型試驗(yàn)對(duì)該數(shù)值模型進(jìn)行驗(yàn)證。結(jié)果表明,該數(shù)值模型可以準(zhǔn)確模擬浮架的變形,采用該數(shù)值模型,分析了不同流速條件下浮架的變形和應(yīng)力,數(shù)值模擬結(jié)果顯示:隨著流速的增加,浮架的變形和應(yīng)力逐漸增加。相同流速條件下,注水下潛后,浮架的應(yīng)力和變形能夠顯著減?。灰虼?,在強(qiáng)流條件下,可以采用浮架注水的方式,使浮架處于下潛狀態(tài),以改善浮架結(jié)構(gòu)的變形和應(yīng)力分布。

    養(yǎng)殖設(shè)施;海洋環(huán)境;有限單元法;動(dòng)力響應(yīng)

    近年來(lái)深海網(wǎng)箱養(yǎng)殖技術(shù)發(fā)展迅速,國(guó)內(nèi)常用的是圓形重力式網(wǎng)箱,其具有成本低、易于操作和便于管理等優(yōu)點(diǎn)。重力式網(wǎng)箱由網(wǎng)衣系統(tǒng)、沉子系統(tǒng)、浮架系統(tǒng)和錨碇系統(tǒng)組成。浮架系統(tǒng)是重力式網(wǎng)箱的重要組成部分,維護(hù)養(yǎng)殖空間的重要構(gòu)件,決定了整個(gè)網(wǎng)箱養(yǎng)殖設(shè)施的安全。在水流作用下浮架有較大的剛體位移,浮架結(jié)構(gòu)也將產(chǎn)生復(fù)雜的三維變形,并產(chǎn)生嚴(yán)重的應(yīng)力集中現(xiàn)象[1],因此,非常有必要對(duì)浮架的應(yīng)力和變形進(jìn)行研究。

    國(guó)內(nèi)外專家采用模型實(shí)驗(yàn)和數(shù)值模擬的方法對(duì)網(wǎng)箱養(yǎng)殖設(shè)施的浮架結(jié)構(gòu)進(jìn)行了廣泛的研究。在浮架水動(dòng)力特性研究方面,將浮架簡(jiǎn)化為水平浮管,分析浮管在波浪作用下的受力和運(yùn)動(dòng)特性[2-3];基于三維水彈性理論分析浮架在波浪作用下的動(dòng)力特性[4-5];采用統(tǒng)計(jì)學(xué)法發(fā)現(xiàn)浮架縱蕩和升沉運(yùn)動(dòng)響應(yīng)隨著波浪頻率的增加而減小[6-7]。在浮架變形和應(yīng)力研究方面,通過(guò)建立有限元模型分析了浮架在荷載作用下的應(yīng)力和變形[8];運(yùn)用曲梁理論建立了圓環(huán)平面內(nèi)和平面外變形的控制方程[9];通過(guò)有限元軟件建模分析了簡(jiǎn)化的浮架雙浮管模型在不同計(jì)算工況準(zhǔn)靜態(tài)荷載施加下的變形-應(yīng)力特性[10]。

    對(duì)于浮架結(jié)構(gòu)的變形和應(yīng)力,前人主要是基于靜態(tài)分析,本文采用有限單元法建立了水流作用下浮架結(jié)構(gòu)的數(shù)值模型,分析動(dòng)荷載條件下浮架結(jié)構(gòu)的變形和受力,研究浮架下潛對(duì)于其變形和受力的影響。

    1 數(shù)值模型

    1.1 浮架結(jié)構(gòu)

    浮架通常處于漂浮狀態(tài),雙浮管是主要提供浮力并承受風(fēng)、浪、流荷載的部件。本文將浮架簡(jiǎn)化成雙浮管結(jié)構(gòu),將上部立管、扶手、豎向支撐和斜向支撐的質(zhì)量都集中到雙浮管上,雙浮管的結(jié)構(gòu)尺寸見(jiàn)表1。雙浮管結(jié)構(gòu)采用HDPE材料,材料的力學(xué)性能參數(shù)采用Fredriksson[11]的試驗(yàn)數(shù)據(jù),彈性模量和屈服強(qiáng)度分別為667 MPa和24.1 MPa,泊松比0.42,密度950 kg/m3。采用自下而上(由線網(wǎng)格創(chuàng)建面網(wǎng)格)的建模方式,使用自由四邊形剖分網(wǎng)格,共包含10 302個(gè)單元和30 870個(gè)節(jié)點(diǎn)。雙浮管結(jié)構(gòu)有限元模型如圖1所示,浮管采用單點(diǎn)錨碇形式固定,水深20 m。

    表1 雙浮管結(jié)構(gòu)幾何參數(shù)

    圖1 雙浮管結(jié)構(gòu)有限元模型Fig.1 Finite element model of double floating pipes structure

    1.2 控制方程

    用有限元?jiǎng)恿W(xué)方法分析浮架結(jié)構(gòu)在瞬態(tài)荷載下的應(yīng)力和變形,基本運(yùn)動(dòng)方程[12-13]如下:

    (1)

    雙浮管單元采用8節(jié)點(diǎn)、6個(gè)自由度的殼單元[14],可考慮拉壓、扭轉(zhuǎn)、彎曲變形,在考慮彎曲效應(yīng)時(shí),與軸向位移自由度有關(guān)的單元質(zhì)量系數(shù)項(xiàng)應(yīng)乘以系數(shù)(Ma/Mt)。Ma與Mt按下式計(jì)算:

    Mt=(mw+mint+madd)L

    (2)

    Ma=(mw+mint)L

    (3)

    (4)

    (5)

    式中:ρ—浮管密度,kg/m3;εin—初應(yīng)變;mint—浮管內(nèi)部流體與附件單位長(zhǎng)度的質(zhì)量,kg/m;CI—外部流體附加質(zhì)量系數(shù);ρw—流體密度,kg/m3。

    采用Newmark方法求解控制方程(1),在每個(gè)時(shí)間步內(nèi)采用修正的Newton-Raphson法進(jìn)行迭代[15]。

    1.3 環(huán)境荷載

    在水流作用下,浮架除了受到水流作用力,還受到垂直方向的重力和浮力,由于水流對(duì)浮架的立柱、扶手、豎向支撐和斜支撐的作用較小,可考慮將這部分重量集成到主浮管,將浮架系統(tǒng)簡(jiǎn)化為一個(gè)雙浮管結(jié)構(gòu)。

    (1)浮架的水流力F采用Morison方程計(jì)算[16]:

    (6)

    (2)浮架單元除受到水流力外,還受到重力G和浮力Ff,其中微元的重力Gi為:

    Gi=G/N

    (7)

    式中:G—浮架總重,N;N—浮架微元總個(gè)數(shù)。

    (3)浮架在運(yùn)動(dòng)過(guò)程中處于漂浮狀態(tài),入水體積隨時(shí)間變化,浮力也隨時(shí)間變化,各時(shí)刻浮力Ffi為:

    (8)

    式中:Si—微元入水橫截面積,m2;li—微元長(zhǎng)度,m;g—重力加速度,m/s2;r—浮管半徑,m;φi—微元入水部分對(duì)應(yīng)的中心角,rad。

    2 模型實(shí)驗(yàn)

    2.1 模型結(jié)構(gòu)布置

    為了驗(yàn)證該數(shù)值模型,將數(shù)值模擬的結(jié)果與模型試驗(yàn)進(jìn)行對(duì)比。圓環(huán)拉伸的模型試驗(yàn)在大連理工大學(xué)海岸和近海工程國(guó)家重點(diǎn)實(shí)驗(yàn)室進(jìn)行,模型放置于水平桌面,圓環(huán)135°和225°處與錨繩相連,錨繩長(zhǎng)22 cm,錨繩另一端施加固定約束,在圓環(huán)0°位置施加水平拉力。共做了12組實(shí)驗(yàn),拉力依次為0.98 N、1.96 N、2.94 N、3.92 N、4.90 N、5.88 N、6.86 N、7.84 N、8.82 N、9.80 N、10.78 N和11.76 N。圓環(huán)幾何參數(shù):環(huán)直徑84.6 cm,浮管外徑和內(nèi)徑分別為15.32 mm和12.66 mm,泊松比0.42,彈性模量6.67×108Pa。

    2.2 模型驗(yàn)證

    數(shù)值模型采用殼單元模擬圓環(huán),桿單元模擬錨繩,采用自由四邊形剖分網(wǎng)格,共包含1 280個(gè)單元和3 840個(gè)節(jié)點(diǎn),在圓環(huán)135°和225°纜繩處施加固定約束,在0°處施加水平荷載。圖2為圓環(huán)在3種拉力(0 N、4.9 N、9.8 N)作用下的物理模型試驗(yàn)和數(shù)值模擬得到的圓環(huán)變形。結(jié)果表明,數(shù)值模擬與物理模型試驗(yàn)獲得的圓環(huán)變形吻合較好,隨著拉力的增加,圓環(huán)沿拉伸方向的變形逐漸增加并呈紡錘形。

    圖2 拉力作用下圓環(huán)的變形Fig.2 The deformation of circle ring under the action of point loads

    給出12組拉力作用下的圓環(huán)的直徑伸長(zhǎng)量,并將數(shù)值模擬結(jié)果與物理模型實(shí)驗(yàn)進(jìn)行定量比較(圖3)。結(jié)果顯示,數(shù)值模擬得到的圓環(huán)徑向變形與物理模型實(shí)驗(yàn)得到的實(shí)測(cè)結(jié)果誤差在10%以內(nèi),數(shù)值模擬能夠較好的模擬圓環(huán)在拉力作用下的變形。

    圖3 拉力作用下圓環(huán)徑向變形Fig.3 The radial deformation of the circular ring under the action of point loads

    3 結(jié)果與討論

    3.1 流速的影響

    采用殼體有限元模型模擬雙浮管在水流作用下的變形和應(yīng)力,分析流速對(duì)浮管的變形和應(yīng)力的影響。共采用4組流速,分別為0.5、1.0、1.5和2.0 m/s,流向?yàn)閄軸正向。圖4表示不同流速下雙浮管結(jié)構(gòu)的變形立面和平面圖,結(jié)果表明:當(dāng)浮架為彈性變形時(shí),雙浮管以平面內(nèi)變形為主,平面外存在復(fù)雜的小變形。流速為0.5 m/s時(shí),浮架平面外幾乎沒(méi)有變形;隨著流速的增加,在未連接錨繩的背流側(cè)發(fā)生了嚴(yán)重的平面外翹曲。

    圖5以浮架在X方向的徑向伸長(zhǎng)量來(lái)定量分析浮架的變形,當(dāng)流速由0.5 m/s增加到2.0 m/s時(shí),浮架徑向伸長(zhǎng)量由23.5 cm增加到30.9 cm。

    (a) v=0.5 m/s (b) v=1.0 m/s (c) v=1.5 m/s (d) v=2.0 m/s圖4 不同流速條件下浮架的變形Fig.4 The deformation of float collar under different flow velocity

    圖5 不同流速條件下浮架X方向徑向伸長(zhǎng)量Fig.5 The X elongation of float collar under different flow velocity

    水流作用下浮架的最大應(yīng)力出現(xiàn)在內(nèi)外浮管的連接構(gòu)件上。圖6為不同流速條件下浮架結(jié)構(gòu)的最大應(yīng)力和浮架節(jié)點(diǎn)應(yīng)力大于2 MPa的節(jié)點(diǎn)個(gè)數(shù)。流速較小時(shí),隨著流速的增加,浮架結(jié)構(gòu)的最大應(yīng)力增加。當(dāng)流速大于1.5 m/s時(shí),浮管的最大應(yīng)力不再增加,但應(yīng)力值大于2 MPa的節(jié)點(diǎn)個(gè)數(shù)增加。流速為2 m/s時(shí),浮架自動(dòng)下潛至水面以下,使得最大應(yīng)力減小。

    3.2 浮架注水下潛對(duì)其變形的影響

    下潛至水面以下的浮架在受到水流作用時(shí),整體的變形和應(yīng)力會(huì)有所不同。采用內(nèi)外浮管注水改變浮架線密度的方法研究浮架下潛對(duì)其變形和應(yīng)力的影響。考慮兩種情況,一種是將3 941.9 kg的水全部注入到內(nèi)外浮管中;另一種是將一半質(zhì)量(1 971.0 kg)的水注入到外浮管中。兩種情況下浮架的線密度分別為24.3 kg/m和13.1 kg/m,未注水的浮架線密度為9.2 kg/m。浮架在纜繩力、浮架重力和浮力的作用下達(dá)到一種平衡狀態(tài),當(dāng)流速增加時(shí),浮架會(huì)自動(dòng)下潛至水面以下,使得纜繩力、浮架重力和浮力達(dá)到新的平衡。本節(jié)采用的流速為2 m/s。

    圖6 不同流速條件下浮架的應(yīng)力Fig.6 The stress of float collar under different flow velocity

    圖7顯示了不同注水情況下浮架在水流作用下整個(gè)下潛的動(dòng)態(tài)過(guò)程,當(dāng)浮架未注水時(shí),在水流的作用下,浮架下潛至30 cm左右,在注一半的水和注滿水的情況下,下潛深度分別達(dá)到72.6 cm和190.1 cm左右。

    圖7 不同注水情況下浮架的下潛深度Fig.7 The submergence depth of float collar under different water injection condition

    圖8給出了不同注水情況下浮架的變形,總體而言浮架處于彈性變形,下潛狀態(tài)的浮架比漂浮狀態(tài)的浮架變形小,浮架平面外也沒(méi)有產(chǎn)生嚴(yán)重的扭轉(zhuǎn)現(xiàn)象。因?yàn)楦〖茏⑺聺撝螅〖芪⒍蔚母×椭亓χg的差值減小,導(dǎo)致作用在浮架上彎曲荷載減小,使得浮架的平面外變形減小;

    此外,流速相同時(shí),注水不同,浮架下潛的深度不同,水面處水質(zhì)點(diǎn)的速度較大,隨著水深的增加而減小,浮架下潛后受到的水流作用力較小。

    圖9給出了不同注水條件下浮架的應(yīng)力值,結(jié)果表明:注水增加,浮架的線密度增加,浮架的最大應(yīng)力和節(jié)點(diǎn)應(yīng)力大于2 MPa的節(jié)點(diǎn)個(gè)數(shù)隨著線密度的增加顯著減??;浮架線密度增加時(shí),浮架最大應(yīng)力減小,從2.91 MPa降到2.77 MPa,同時(shí)節(jié)點(diǎn)應(yīng)力大于2 MPa的個(gè)數(shù)也相應(yīng)減少,應(yīng)力分布更加均勻。

    (a) 未注水 (b) 注一半水 (c) 注滿水圖8 不同注水情況下浮架的變形Fig.8 The deformation of float collar under different water injection condition

    圖9 不同注水情況的應(yīng)力值Fig.9 The stress value of float collar under different water injection condition

    4 結(jié)論

    采用殼體有限元方法建立浮架結(jié)構(gòu)的數(shù)值模型,并采用物理模型試驗(yàn)對(duì)該數(shù)值模型進(jìn)行了驗(yàn)證,結(jié)果表明:該數(shù)值模型能夠較好的模擬浮架結(jié)構(gòu)的變形。采用該數(shù)值模型分析了不同流速條件下浮架變形和受力,研究了浮架下潛對(duì)其變形和應(yīng)力的影響,結(jié)果表明:當(dāng)浮架為彈性變形時(shí),雙浮管以平面內(nèi)變形為主,平面外存在復(fù)雜的小變形;流速較小時(shí),隨著流速的增加,浮架結(jié)構(gòu)的最大應(yīng)力增加。當(dāng)流速大于1.5m/s時(shí),浮管的最大應(yīng)力不再增加,但應(yīng)力值大于2MPa的節(jié)點(diǎn)個(gè)數(shù)增加。浮架的下潛能夠有效改善浮架的變形和應(yīng)力分布,采用浮架注水的方式使浮架下潛是一種較為合理的提高深水網(wǎng)箱養(yǎng)殖設(shè)施安全性的措施。

    [1] 桂福坤, 王煒霞, 張懷慧. 網(wǎng)箱工程發(fā)展現(xiàn)狀及展望[J]. 大連水產(chǎn)學(xué)院學(xué)報(bào), 2002(1):70-78.

    [2] NEWMAN J N. Marine hydrodynamics [M]. Massachusetts:

    The MIT Press, 1977:357-362.

    [3] FALTINSEN O M. Sea loads on ships and offshore structure [M]. Cambridge:Cambridge University Press, 1990:174-186.

    [4] FU S X,MOAN T. Dynamic analyses of floating fish cage collars in waves [J]. Aquacultural Engineering, 2012, 47:7-15

    [5] FU S X,MOAN T,CHEN X J,etal. Hydroelastic analysis of flexible floating interconnected structures[J]. Ocean Engineering, 2007, 34:1516-1531.

    [6] XU T J,DONG G H,ZHAO Y P,etal. Hydrodynamic characteristics of floating pipes in random waves [C]// Proceedings of the 32nd International Conference on Ocean, Offshore and Arctic Engineering, 2013.

    [7] 許條建. 離岸組合式網(wǎng)箱水動(dòng)力特性研究[D]. 大連:大連理工大學(xué), 2014.

    [8] 郭根喜. 圓形網(wǎng)箱浮架結(jié)構(gòu)有限元分析[J]. 湛江:廣東海洋大學(xué)學(xué)報(bào), 2011, 31(4):69-74.

    [9] 郝雙戶. 深水重力式網(wǎng)箱浮架結(jié)構(gòu)的流固耦合研究[D]. 大連:大連理工大學(xué), 2008.

    [10]王笛清. 重力式網(wǎng)箱浮架結(jié)構(gòu)變形-應(yīng)力特性有限元分析[D]. 大連:大連理工大學(xué), 2014.

    [11]FREDRIKSSON D W, TSUKROV I. Development of structural modeling techniques for evaluating HDPE plastic net pens used in marine aquaculture [J]. Ocean Engineering, 2007, 34(16):2124-2137.

    [12]崔勇,關(guān)長(zhǎng)濤,萬(wàn)榮,等.基于有限元方法對(duì)波流場(chǎng)中養(yǎng)殖網(wǎng)箱的系統(tǒng)動(dòng)力分析[J].北京:工程力學(xué), 2010(5):250-256.

    [13]CUI Y,GUAN C T,WAN R,etal. Numerical simulation of a flatfish cage system in waves and currents[J]. Aquacultural Engineering, 2013, 56(1):26-33.

    [14]王新敏. ANSYS工程結(jié)構(gòu)數(shù)值分析[M]. 北京:人民交通出版社, 2007:11-12.

    [15]SUHEY D J, KIM H N, NIEZRECKI C.Numercial modeling and design of inflatable structures-Application to open-ocean-aquaculture cages [J]. Aquacultural Engineering, 2005, 33:285-303.

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

    Study on the dynamic response of float collar of gravity cages subject to water currents

    LI Kunpeng1, XU Tiaojian1,2, DONG Guohai1, ZHAO Yunpeng1

    (1 State Key Laboratory of Coastal and Offshore Engineering, Dalian University of Technology, Dalian 116024, China 2 School of Naval Architecture and Ocean Engineering, Dalian University of Technology, Dalian 116024, China)

    The offshore culture cage has been widely applied in the open sea recently. Floating collar is an important component for the fish cage structure, and the safety of floating collar is very important for the design of fish cage. The numerical model of floating collar was built by using the SHELL element based on the finite element method, the stress and deformation of the floating collar in current were simulated, and the numerical model were verified through the physical model experiments. The results showed that the numerical model could accurately simulate the deformation of the floating collar. Based on the numerical model of the floating collar, the stress and deformation of the floating collar under different water velocity were analyzed. The numerical simulation results indicated that the stress and deformation of floating collar were increased gradually with the increase of the velocity. The maximum deformation and stress of the floating collar in the floating condition was significantly larger than that for the floating collar in the submerge condition. In conclusion, the deformation and stress distribution of floating collar under high water current could be improved obviously when the floating collar was injected with water to make it be in submersible condition.

    culture facilities; marine environment; finite element method; dynamic response

    2015-10-10

    2016-01-19

    國(guó)家自然科學(xué)基金項(xiàng)目“網(wǎng)箱養(yǎng)殖設(shè)施非線性動(dòng)力響應(yīng)及其破壞機(jī)理(No.51409037、No.51239002)”;博士后基金項(xiàng)目“復(fù)雜荷載作用下網(wǎng)箱的力學(xué)特性及破壞機(jī)理(No.2014M560211、No.2015T80254)”

    李坤鵬(1990—),女,碩士研究生,研究方向:深水網(wǎng)箱結(jié)構(gòu)的動(dòng)力響應(yīng)。E-mail:likunpeng2623@163.com

    許條建(1983—),男,碩導(dǎo),博士,研究方向:海洋牧場(chǎng)結(jié)構(gòu)水動(dòng)力特性。E-mail:tjxu@dlut.edu.cn

    10.3969/j.issn.1007-9580.2016.01.006

    S953.4

    A

    1007-9580(2016)01-030-06

    猜你喜歡
    架結(jié)構(gòu)網(wǎng)箱圓環(huán)
    海洋平臺(tái)板架結(jié)構(gòu)簡(jiǎn)化評(píng)估對(duì)比研究
    某輕卡車(chē)架結(jié)構(gòu)強(qiáng)度性能優(yōu)化
    加權(quán)全能量最小的圓環(huán)形變
    豬圓環(huán)病毒病的發(fā)生、診斷和防治
    一例鴨圓環(huán)病毒病的診斷
    叉車(chē)門(mén)架結(jié)構(gòu)分析和定期檢驗(yàn)探討
    圓環(huán)上的覆蓋曲面不等式及其應(yīng)用
    10kV配網(wǎng)箱變局放檢測(cè)案例分析
    電子制作(2017年24期)2017-02-02 07:14:45
    全國(guó)首例網(wǎng)箱豢養(yǎng)江豚繁育成功
    池塘小網(wǎng)箱養(yǎng)殖黃鱔
    久久久久国内视频| 日韩免费av在线播放| 国内久久婷婷六月综合欲色啪| 精品午夜福利视频在线观看一区| 美女免费视频网站| av免费在线观看网站| 亚洲第一电影网av| 国产午夜精品久久久久久| 国产免费av片在线观看野外av| 成人特级黄色片久久久久久久| 亚洲av中文字字幕乱码综合 | 神马国产精品三级电影在线观看 | 91字幕亚洲| 18禁国产床啪视频网站| 久久精品成人免费网站| 国产在线精品亚洲第一网站| 亚洲中文字幕一区二区三区有码在线看 | 一本久久中文字幕| 最近最新中文字幕大全电影3 | 51午夜福利影视在线观看| 国内精品久久久久久久电影| 日日摸夜夜添夜夜添小说| 日韩中文字幕欧美一区二区| 国产精品,欧美在线| 人成视频在线观看免费观看| 久久伊人香网站| 又黄又爽又免费观看的视频| 久久精品aⅴ一区二区三区四区| 精品免费久久久久久久清纯| 午夜老司机福利片| 国产伦在线观看视频一区| 日韩精品青青久久久久久| 国产精品二区激情视频| 看片在线看免费视频| 在线国产一区二区在线| 中文字幕高清在线视频| 色播在线永久视频| 成年版毛片免费区| 国产精品乱码一区二三区的特点| 亚洲熟妇熟女久久| 亚洲精品美女久久av网站| 中文字幕人妻丝袜一区二区| 日本精品一区二区三区蜜桃| 亚洲精品色激情综合| 精品国产美女av久久久久小说| 久久中文看片网| avwww免费| 国产精品自产拍在线观看55亚洲| 丝袜人妻中文字幕| 免费在线观看成人毛片| 中文字幕精品免费在线观看视频| 人人妻人人澡欧美一区二区| 欧美成人性av电影在线观看| 校园春色视频在线观看| 变态另类成人亚洲欧美熟女| 丰满的人妻完整版| 麻豆国产av国片精品| 国产成人av激情在线播放| 欧美亚洲日本最大视频资源| 国产av不卡久久| av福利片在线| 9191精品国产免费久久| 波多野结衣av一区二区av| 欧美日韩亚洲综合一区二区三区_| 两人在一起打扑克的视频| 国产伦在线观看视频一区| 久久国产乱子伦精品免费另类| av在线天堂中文字幕| 成人国语在线视频| 在线十欧美十亚洲十日本专区| 啦啦啦免费观看视频1| а√天堂www在线а√下载| 久久草成人影院| 欧美性长视频在线观看| 美女免费视频网站| 在线观看免费日韩欧美大片| 国产成人av激情在线播放| 悠悠久久av| 欧美激情极品国产一区二区三区| bbb黄色大片| 一区二区三区国产精品乱码| 国产又爽黄色视频| 在线观看免费日韩欧美大片| 9191精品国产免费久久| 久久国产精品男人的天堂亚洲| 亚洲aⅴ乱码一区二区在线播放 | 国产精品香港三级国产av潘金莲| 国产精品电影一区二区三区| 中文字幕久久专区| 又黄又粗又硬又大视频| 国产亚洲av高清不卡| 亚洲av成人不卡在线观看播放网| 在线观看一区二区三区| 美女高潮喷水抽搐中文字幕| 黄频高清免费视频| 国产成人精品无人区| 亚洲aⅴ乱码一区二区在线播放 | 欧美黄色片欧美黄色片| 国产精品爽爽va在线观看网站 | 国产真实乱freesex| 最新在线观看一区二区三区| 亚洲avbb在线观看| 最近最新中文字幕大全电影3 | 亚洲成人免费电影在线观看| 在线观看66精品国产| 亚洲国产欧美日韩在线播放| 村上凉子中文字幕在线| 亚洲第一青青草原| 国产精品精品国产色婷婷| 欧美色欧美亚洲另类二区| 中文字幕精品免费在线观看视频| 婷婷六月久久综合丁香| 亚洲电影在线观看av| 色综合亚洲欧美另类图片| 日韩欧美免费精品| 最近最新中文字幕大全免费视频| 精品国产乱子伦一区二区三区| 国产一区二区三区在线臀色熟女| 真人一进一出gif抽搐免费| 国产在线精品亚洲第一网站| 又黄又爽又免费观看的视频| 一区二区三区高清视频在线| 午夜福利在线在线| 老司机福利观看| 中文在线观看免费www的网站 | 久久香蕉激情| √禁漫天堂资源中文www| 看片在线看免费视频| 免费无遮挡裸体视频| 丝袜在线中文字幕| 亚洲成人精品中文字幕电影| 免费在线观看完整版高清| 亚洲一码二码三码区别大吗| av片东京热男人的天堂| 人妻丰满熟妇av一区二区三区| 正在播放国产对白刺激| 日日爽夜夜爽网站| 亚洲av中文字字幕乱码综合 | 99re在线观看精品视频| avwww免费| or卡值多少钱| 一本久久中文字幕| 男人操女人黄网站| 欧美丝袜亚洲另类 | 中出人妻视频一区二区| 国产成人啪精品午夜网站| 免费在线观看黄色视频的| 成年免费大片在线观看| 十分钟在线观看高清视频www| 最新在线观看一区二区三区| 成人一区二区视频在线观看| 亚洲专区国产一区二区| 亚洲av成人一区二区三| 一本综合久久免费| 黄网站色视频无遮挡免费观看| 村上凉子中文字幕在线| 午夜福利欧美成人| 精品国产美女av久久久久小说| 女生性感内裤真人,穿戴方法视频| 99国产综合亚洲精品| 黄片大片在线免费观看| 午夜激情福利司机影院| 制服人妻中文乱码| 99久久99久久久精品蜜桃| 国产精品免费视频内射| 午夜福利欧美成人| 老鸭窝网址在线观看| 50天的宝宝边吃奶边哭怎么回事| 国产不卡一卡二| 午夜福利视频1000在线观看| 搡老熟女国产l中国老女人| 99re在线观看精品视频| 亚洲成人久久爱视频| 久久久久国产精品人妻aⅴ院| 国产精品,欧美在线| 午夜激情福利司机影院| 动漫黄色视频在线观看| 俄罗斯特黄特色一大片| 正在播放国产对白刺激| 亚洲国产毛片av蜜桃av| 在线观看日韩欧美| 欧美成人午夜精品| 久久精品国产亚洲av高清一级| 黑人欧美特级aaaaaa片| 搡老岳熟女国产| 两性夫妻黄色片| 女人爽到高潮嗷嗷叫在线视频| 婷婷六月久久综合丁香| 可以在线观看的亚洲视频| 韩国av一区二区三区四区| 国产亚洲精品第一综合不卡| 亚洲熟妇熟女久久| 女警被强在线播放| 亚洲av中文字字幕乱码综合 | 香蕉丝袜av| 亚洲成av人片免费观看| 精品午夜福利视频在线观看一区| 中文字幕高清在线视频| 嫩草影视91久久| 一级毛片精品| 51午夜福利影视在线观看| 在线观看免费午夜福利视频| 成人18禁在线播放| www日本在线高清视频| 亚洲国产欧美日韩在线播放| 在线观看一区二区三区| 在线天堂中文资源库| 又大又爽又粗| svipshipincom国产片| 一区二区三区国产精品乱码| 亚洲va日本ⅴa欧美va伊人久久| 身体一侧抽搐| 18禁黄网站禁片免费观看直播| 在线免费观看的www视频| 18美女黄网站色大片免费观看| 欧美色视频一区免费| 欧美日韩福利视频一区二区| 中文在线观看免费www的网站 | 日韩国内少妇激情av| 精品一区二区三区视频在线观看免费| 美女免费视频网站| 日日干狠狠操夜夜爽| e午夜精品久久久久久久| 精品第一国产精品| 久久久久久久久中文| 久热爱精品视频在线9| 午夜福利在线在线| 亚洲精品美女久久久久99蜜臀| 后天国语完整版免费观看| 免费看日本二区| 国产激情欧美一区二区| 日本熟妇午夜| 男人舔奶头视频| 一个人免费在线观看的高清视频| 麻豆成人午夜福利视频| 国产精品电影一区二区三区| 国产在线精品亚洲第一网站| 色婷婷久久久亚洲欧美| 在线视频色国产色| 香蕉丝袜av| 久久这里只有精品19| 亚洲中文日韩欧美视频| 欧美又色又爽又黄视频| 国产又色又爽无遮挡免费看| 久久亚洲真实| 午夜久久久在线观看| 成人18禁高潮啪啪吃奶动态图| 亚洲熟女毛片儿| 久久久久久国产a免费观看| 丝袜人妻中文字幕| 欧美国产日韩亚洲一区| 亚洲 欧美 日韩 在线 免费| 国产精品一区二区三区四区久久 | 人人妻人人澡人人看| 天堂√8在线中文| 色播在线永久视频| 精品不卡国产一区二区三区| av在线播放免费不卡| 在线天堂中文资源库| 欧美成狂野欧美在线观看| 99久久综合精品五月天人人| 淫秽高清视频在线观看| 99国产精品一区二区蜜桃av| 亚洲国产中文字幕在线视频| 人妻丰满熟妇av一区二区三区| 此物有八面人人有两片| 久久香蕉激情| 精品不卡国产一区二区三区| 性色av乱码一区二区三区2| 丁香六月欧美| 哪里可以看免费的av片| 在线播放国产精品三级| 男女下面进入的视频免费午夜 | 久久久久久久久中文| 美女高潮喷水抽搐中文字幕| 制服人妻中文乱码| 国产蜜桃级精品一区二区三区| 99久久精品国产亚洲精品| 亚洲精品国产区一区二| 久久九九热精品免费| 国产精品99久久99久久久不卡| 国产精品精品国产色婷婷| 一级毛片女人18水好多| 亚洲熟妇中文字幕五十中出| 欧美久久黑人一区二区| 欧美激情 高清一区二区三区| 亚洲人成伊人成综合网2020| www日本在线高清视频| 亚洲国产精品sss在线观看| 日本免费a在线| 国产国语露脸激情在线看| 99精品在免费线老司机午夜| ponron亚洲| 亚洲人成77777在线视频| 久久久久久亚洲精品国产蜜桃av| 亚洲 欧美一区二区三区| 动漫黄色视频在线观看| 国产成+人综合+亚洲专区| 亚洲成人久久性| 免费女性裸体啪啪无遮挡网站| 精品第一国产精品| 老汉色∧v一级毛片| 亚洲精品色激情综合| 国产av又大| 中文字幕高清在线视频| 免费在线观看完整版高清| www.精华液| 亚洲国产精品sss在线观看| 精品一区二区三区四区五区乱码| 日本免费一区二区三区高清不卡| 亚洲人成电影免费在线| 成人av一区二区三区在线看| 男人操女人黄网站| 老司机靠b影院| 高清毛片免费观看视频网站| 黑人巨大精品欧美一区二区mp4| 日韩av在线大香蕉| 亚洲专区字幕在线| 香蕉丝袜av| 91在线观看av| 国内毛片毛片毛片毛片毛片| 又黄又粗又硬又大视频| 亚洲av第一区精品v没综合| 99久久99久久久精品蜜桃| av有码第一页| 精品高清国产在线一区| 90打野战视频偷拍视频| 最近最新中文字幕大全电影3 | 熟女少妇亚洲综合色aaa.| 国产激情欧美一区二区| 99久久久亚洲精品蜜臀av| 亚洲男人天堂网一区| 99在线视频只有这里精品首页| av电影中文网址| 精品久久久久久,| 男女视频在线观看网站免费 | 两人在一起打扑克的视频| 欧美日本视频| 天天躁狠狠躁夜夜躁狠狠躁| 少妇裸体淫交视频免费看高清 | 搡老妇女老女人老熟妇| 男男h啪啪无遮挡| 女生性感内裤真人,穿戴方法视频| 啦啦啦观看免费观看视频高清| 日日爽夜夜爽网站| 久久99热这里只有精品18| 91麻豆av在线| 亚洲片人在线观看| 十分钟在线观看高清视频www| 精品人妻1区二区| 中文字幕精品亚洲无线码一区 | 18禁裸乳无遮挡免费网站照片 | 国产成人精品无人区| 99国产精品一区二区三区| 成人亚洲精品一区在线观看| 中文资源天堂在线| 淫秽高清视频在线观看| 久久久久免费精品人妻一区二区 | 最新美女视频免费是黄的| 国产欧美日韩一区二区三| svipshipincom国产片| 亚洲精品久久成人aⅴ小说| www.精华液| 久热这里只有精品99| √禁漫天堂资源中文www| 久久精品aⅴ一区二区三区四区| 美女扒开内裤让男人捅视频| 免费av毛片视频| 国产人伦9x9x在线观看| 国产精品久久久久久精品电影 | 老汉色av国产亚洲站长工具| 人妻久久中文字幕网| 久久精品人妻少妇| 最近最新免费中文字幕在线| 精品第一国产精品| 久久久久久九九精品二区国产 | 岛国视频午夜一区免费看| 操出白浆在线播放| 一区二区三区国产精品乱码| 国产精品香港三级国产av潘金莲| 黄片大片在线免费观看| 国产99久久九九免费精品| 一进一出抽搐gif免费好疼| 亚洲精品色激情综合| 久99久视频精品免费| 9191精品国产免费久久| 一区福利在线观看| 久久性视频一级片| 热re99久久国产66热| 韩国精品一区二区三区| 色精品久久人妻99蜜桃| 久久人妻福利社区极品人妻图片| 在线播放国产精品三级| 99re在线观看精品视频| 狂野欧美激情性xxxx| 黄色 视频免费看| cao死你这个sao货| 国产亚洲精品av在线| 日本精品一区二区三区蜜桃| 欧美成人免费av一区二区三区| 欧美精品亚洲一区二区| 久久久国产精品麻豆| 国产亚洲欧美98| 国产爱豆传媒在线观看 | 国产精品九九99| 黄色丝袜av网址大全| av天堂在线播放| 青草久久国产| 少妇熟女aⅴ在线视频| 久久久国产精品麻豆| 最新美女视频免费是黄的| 亚洲第一av免费看| 在线观看日韩欧美| 久久久久亚洲av毛片大全| 最新在线观看一区二区三区| 国产精品香港三级国产av潘金莲| 最近最新中文字幕大全电影3 | 法律面前人人平等表现在哪些方面| 国产免费av片在线观看野外av| 最新美女视频免费是黄的| 亚洲av成人av| 欧美日本视频| 色播亚洲综合网| 国产精品亚洲一级av第二区| 99精品久久久久人妻精品| 亚洲全国av大片| 午夜精品久久久久久毛片777| 亚洲人成网站在线播放欧美日韩| 香蕉丝袜av| 我的亚洲天堂| 精品国产一区二区三区四区第35| 高清毛片免费观看视频网站| 国产片内射在线| 亚洲精品一区av在线观看| 日韩av在线大香蕉| 一区二区三区国产精品乱码| 亚洲精品美女久久久久99蜜臀| 精品久久久久久久久久久久久 | 成人一区二区视频在线观看| 免费电影在线观看免费观看| 久久 成人 亚洲| 淫秽高清视频在线观看| 久久久久国产精品人妻aⅴ院| 中国美女看黄片| 身体一侧抽搐| 桃色一区二区三区在线观看| 啦啦啦免费观看视频1| 国产av一区在线观看免费| 波多野结衣高清作品| 日日摸夜夜添夜夜添小说| 欧美激情极品国产一区二区三区| 国内精品久久久久精免费| 国产精品久久久人人做人人爽| 欧美日韩亚洲国产一区二区在线观看| 久久久久亚洲av毛片大全| 亚洲色图 男人天堂 中文字幕| 婷婷亚洲欧美| 免费看美女性在线毛片视频| 成人一区二区视频在线观看| 国产又爽黄色视频| 很黄的视频免费| 91成年电影在线观看| 亚洲久久久国产精品| 国内久久婷婷六月综合欲色啪| 亚洲国产精品999在线| 国产av又大| 亚洲熟女毛片儿| 亚洲av第一区精品v没综合| 亚洲精品中文字幕一二三四区| 桃色一区二区三区在线观看| 亚洲全国av大片| 一进一出抽搐动态| 精品一区二区三区视频在线观看免费| a级毛片在线看网站| 婷婷精品国产亚洲av| 亚洲aⅴ乱码一区二区在线播放 | 人人妻人人看人人澡| 亚洲人成77777在线视频| 国产精品亚洲美女久久久| 搡老妇女老女人老熟妇| 亚洲国产精品999在线| aaaaa片日本免费| 神马国产精品三级电影在线观看 | 亚洲精品国产一区二区精华液| 国产精品香港三级国产av潘金莲| 午夜福利成人在线免费观看| 久久精品夜夜夜夜夜久久蜜豆 | 久久久久久亚洲精品国产蜜桃av| 久久欧美精品欧美久久欧美| 午夜老司机福利片| 97碰自拍视频| 美女 人体艺术 gogo| 伊人久久大香线蕉亚洲五| 一级作爱视频免费观看| 亚洲全国av大片| 很黄的视频免费| 神马国产精品三级电影在线观看 | www.熟女人妻精品国产| 久久亚洲精品不卡| 成人三级做爰电影| 国产在线观看jvid| 夜夜夜夜夜久久久久| 无遮挡黄片免费观看| 国产日本99.免费观看| av天堂在线播放| 亚洲人成伊人成综合网2020| 日韩欧美国产一区二区入口| 国内精品久久久久久久电影| 国产精品久久久久久亚洲av鲁大| 在线观看免费午夜福利视频| 亚洲全国av大片| 99热6这里只有精品| 人人妻人人澡欧美一区二区| 国产日本99.免费观看| 色在线成人网| 亚洲男人的天堂狠狠| 欧美黄色片欧美黄色片| 成人永久免费在线观看视频| 禁无遮挡网站| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲在线自拍视频| 18禁美女被吸乳视频| 久久久久久免费高清国产稀缺| 日本撒尿小便嘘嘘汇集6| 一进一出好大好爽视频| 久久久久久久久中文| 最近最新免费中文字幕在线| 神马国产精品三级电影在线观看 | 午夜两性在线视频| 最新在线观看一区二区三区| 岛国视频午夜一区免费看| 国产成+人综合+亚洲专区| 大香蕉久久成人网| 亚洲国产精品sss在线观看| 国产成人啪精品午夜网站| 丰满的人妻完整版| 精品第一国产精品| АⅤ资源中文在线天堂| 色综合婷婷激情| 免费女性裸体啪啪无遮挡网站| 日韩精品中文字幕看吧| 久久精品成人免费网站| 亚洲一区二区三区色噜噜| 亚洲一卡2卡3卡4卡5卡精品中文| 一进一出抽搐gif免费好疼| 在线观看免费午夜福利视频| 18禁观看日本| x7x7x7水蜜桃| 日韩中文字幕欧美一区二区| 99re在线观看精品视频| 国产v大片淫在线免费观看| 欧美在线黄色| 成人三级做爰电影| 在线永久观看黄色视频| 精品久久久久久久久久免费视频| 天堂√8在线中文| 色综合婷婷激情| 中文字幕另类日韩欧美亚洲嫩草| 99国产精品一区二区三区| 亚洲黑人精品在线| 中文亚洲av片在线观看爽| 真人一进一出gif抽搐免费| 老熟妇仑乱视频hdxx| 午夜激情福利司机影院| 亚洲国产精品sss在线观看| 亚洲国产中文字幕在线视频| 亚洲激情在线av| 精品久久久久久久久久免费视频| 国产av不卡久久| 色综合婷婷激情| 亚洲欧美一区二区三区黑人| 亚洲欧美精品综合一区二区三区| 国产精华一区二区三区| 亚洲成人久久性| x7x7x7水蜜桃| 日韩 欧美 亚洲 中文字幕| 黄色丝袜av网址大全| 长腿黑丝高跟| 久久 成人 亚洲| 亚洲黑人精品在线| 99国产极品粉嫩在线观看| 日本一本二区三区精品| 久久中文字幕人妻熟女| 久久久久国产精品人妻aⅴ院| 高清毛片免费观看视频网站| 身体一侧抽搐| 中文字幕精品亚洲无线码一区 | 国产91精品成人一区二区三区| 人妻丰满熟妇av一区二区三区| 午夜福利免费观看在线| 无限看片的www在线观看| 免费女性裸体啪啪无遮挡网站| 亚洲最大成人中文| 最近最新免费中文字幕在线| 国产一级毛片七仙女欲春2 | 亚洲专区中文字幕在线| 精品一区二区三区四区五区乱码| 亚洲美女黄片视频| 99热这里只有精品一区 | 麻豆国产av国片精品| 美女午夜性视频免费| 真人做人爱边吃奶动态| а√天堂www在线а√下载| 国产亚洲av高清不卡| 操出白浆在线播放| 在线视频色国产色| 亚洲自拍偷在线| 国产精品野战在线观看| 精品久久久久久久末码| 亚洲国产精品999在线| 精品久久久久久久人妻蜜臀av| 久久精品成人免费网站| 97人妻精品一区二区三区麻豆 | 国产av不卡久久| 啦啦啦观看免费观看视频高清|