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

    考慮流固耦合的裝配式管線水擊壓力波速計(jì)算研究

    2023-11-09 06:49:36陳明梁澤山彭小慶許亞能付江華
    機(jī)床與液壓 2023年20期
    關(guān)鍵詞:水擊膠圈密封膠

    陳明,梁澤山,彭小慶,許亞能,付江華

    (1.重慶賽寶工業(yè)技術(shù)研究院有限公司,重慶 401332;2.重慶理工大學(xué)車輛工程學(xué)院,重慶 400054)

    0 前言

    裝配式管線由進(jìn)出口兩端均帶有卡槽的鋼管、管路附件以及各種配套設(shè)備組合裝配而成,具有機(jī)動(dòng)靈活、展收迅速、維修方便等優(yōu)點(diǎn),已經(jīng)在化工工業(yè)、農(nóng)業(yè)灌溉、市政工程以及軍事后勤保障,尤其是油料應(yīng)急輸送保障等領(lǐng)域得到了廣泛的應(yīng)用[1]。

    在管輸過程中,水擊是經(jīng)常出現(xiàn)的一種瞬態(tài)工況,嚴(yán)重時(shí)會(huì)導(dǎo)致超常壓力脈動(dòng),在流固耦合作用(Fluid-Structure Interaction,F(xiàn)SI)下將進(jìn)一步誘發(fā)管道結(jié)構(gòu)振動(dòng)甚至破壞,因此應(yīng)該進(jìn)行全面地分析、預(yù)防和控制[2-4]。水擊壓力波速是水擊分析中非常重要和敏感的參數(shù),受到流體壓力、溫度、含氣量以及管道結(jié)構(gòu)尺寸、管端約束形式、管材等多種因素的綜合影響。準(zhǔn)確計(jì)算管道水擊壓力波速一直是水擊分析領(lǐng)域的熱點(diǎn)。

    周云龍等[5]分析了含氣量對(duì)漿體水擊過程的影響,提出了復(fù)合管道氣液固三相流的漿體水擊波速計(jì)算公式。吳迪等人[6]提出了內(nèi)置軟管的水平管道水擊波速公式,并對(duì)影響水擊波速的因素進(jìn)行了敏感性分析。奚斌等人[7]針對(duì)環(huán)形斷面管道的結(jié)構(gòu)特點(diǎn)和流體運(yùn)動(dòng)特性,推導(dǎo)了環(huán)形斷面管道水擊波速的計(jì)算公式,并分析了液流摩阻項(xiàng)對(duì)公式的影響。WAN和 MAO[8]基于等效彈性模量法,提出了一個(gè)適用于鋼絲骨架復(fù)合型PE管的水擊波速計(jì)算式,并分析了鋼絲骨架對(duì)瞬變響應(yīng)過程的影響。史富全[9]綜合考慮了井筒流體含氣率、氣液兩相流型、固相類型等對(duì)水擊波速的影響,推導(dǎo)了水擊波在氣、液、固三相混合流體中傳播速度計(jì)算公式。郭強(qiáng)等人[10]在薄壁管流固耦合模型的基礎(chǔ)上,建立了適用于厚壁管流固耦合分析的一維數(shù)學(xué)模型,并研究了管道不同厚徑比對(duì)壓力波速與應(yīng)力波速的影響。REZAPOUR和 RIASI[11]以不同濃度的聚丙烯酰胺溶液作為黏彈性流體,在銅管和聚乙烯管中,分別研究了流體的黏彈性對(duì)水擊壓力波的衰減以及波周期的影響。LI等[12]對(duì)液體火箭推進(jìn)劑輸送管道中的低溫冷凝兩相流壓力波傳播問題進(jìn)行了研究,深入分析了流體壓力波的傳播與衰減對(duì)不穩(wěn)定頻率的依賴關(guān)系。KUBRAK和KODURA[13]研究了在管道中分別置入薄壁圓管、厚壁圓管以及實(shí)心圓管時(shí)的水擊壓力波速計(jì)算方法,并通過進(jìn)一步分析認(rèn)為內(nèi)置低體積彈性模量的圓管可對(duì)水擊過程產(chǎn)生阻尼效應(yīng)。胡曉東等[14]搭建了水擊壓力波試驗(yàn)系統(tǒng),對(duì)4種水擊壓力波速計(jì)算方法進(jìn)行了研究與對(duì)比,以平均值、標(biāo)準(zhǔn)差和變異系數(shù)為指標(biāo),分析了4種波速計(jì)算方法的數(shù)值結(jié)果與穩(wěn)定性。

    綜上,目前對(duì)管道水擊壓力波速的研究主要集中在管輸介質(zhì)、管壁材料、管道截面形狀等因素對(duì)波速的影響,而有關(guān)管道連接方式對(duì)波速的影響研究卻鮮有見到。裝配式管線各管段、管件之間均采用連接器加膠圈的連接結(jié)構(gòu),密封膠圈的彈性模量遠(yuǎn)小于鋼材,兩者相差105數(shù)量級(jí)。裝配式管線上數(shù)量眾多的密封膠圈必然會(huì)改變整條管線的彈性模量,從而對(duì)水擊波速產(chǎn)生明顯影響。雖然蒲家寧[15]、蔣明等人[16]研究了裝配式管線中密封膠圈對(duì)水擊響應(yīng)的影響,但均未計(jì)及管線系統(tǒng)中的FSI效應(yīng)。鑒于實(shí)際管線系統(tǒng)中存在流體與管路結(jié)構(gòu)的耦合互動(dòng)作用,考慮FSI效應(yīng)的水擊響應(yīng)更符合實(shí)際情況。因此,本文作者對(duì)裝配式管線中的耦合水擊波速計(jì)算方法進(jìn)行研究,以期進(jìn)一步為準(zhǔn)確地分析裝配式管線系統(tǒng)耦合水擊響應(yīng)奠定基礎(chǔ)。

    1 流固耦合水擊動(dòng)力學(xué)模型

    在管道水擊分析研究領(lǐng)域常采用薄壁管假設(shè),以此可忽略徑向應(yīng)力的影響。但在實(shí)際的管輸系統(tǒng)中,有時(shí)難以界定究竟是薄壁管還是厚壁管。因此,文中不做薄壁管假設(shè)??紤]流體與管道之間的耦合互動(dòng)效應(yīng),耦合水擊動(dòng)力學(xué)模型[17]可表示為

    (1)

    (2)

    (3)

    (4)

    式中:v為流體流速;p為管內(nèi)流體壓力;U為管道平均軸向振速;σ為管壁平均軸向應(yīng)力;K為流體體積彈性模量;Ep為管材的彈性模量;μp為管材的泊松比;D為管道內(nèi)直徑;e為管壁厚度;ρf為流體密度;ρp為管材的密度;g為重力加速度;θ為管道傾角;x為管道軸向;t為時(shí)間。

    上述FSI四方程模型屬于擬線性雙曲型偏微分方程組,其特征方程具有4個(gè)不相等的實(shí)根,即考慮FSI效應(yīng)的水擊波速Cf和軸向應(yīng)力波速Cp:

    (5)

    (6)

    從式(5)(6)可以看出,考慮FSI效應(yīng)時(shí),耦合水擊波速與經(jīng)典水擊波速及軸向應(yīng)力波速均相關(guān),而管材彈性模量Ep則是決定水擊波速和軸向應(yīng)力波速的重要參數(shù)之一。密封膠圈的存在必然會(huì)改變整條裝配式管線的彈性模量,進(jìn)而對(duì)耦合水擊波速造成影響。

    2 考慮密封膠圈影響的耦合水擊波速計(jì)算

    2.1 密封膠圈的形變分析

    裝配式管線各管段、管件之間由卡箍式連接器連接、膠圈密封,其連接結(jié)構(gòu)如圖1所示。當(dāng)受壓液體充滿膠圈內(nèi)的環(huán)形槽時(shí),膠圈膨脹并緊緊地貼在管子卡槽的凸緣上,起到密封作用。

    圖1 連接器(a)、膠圈(b)及連接結(jié)構(gòu)(c)示意

    根據(jù)圖1(c)所示的連接結(jié)構(gòu),密封膠圈微元段的受力情況如圖2所示。

    密封膠圈受到管內(nèi)流體壓力p作用產(chǎn)生形變。在Z方向上(即膠圈圓周方向),膠圈微元段因泊松效應(yīng)而伸長(zhǎng),但變形受到連接器的約束,產(chǎn)生擠壓力pZ。膠圈在三向應(yīng)力的作用下,體積Vr將發(fā)生變化,從體積變化量ΔVr相等的角度,可將膠圈形狀及受力狀態(tài)簡(jiǎn)化,如圖3所示。

    圖3 膠圈當(dāng)量體及受力狀態(tài)

    密封膠圈由丁腈耐油橡膠制成,不是線彈性材料,但在小變形范圍內(nèi)可近似認(rèn)為其滿足虎克定律,應(yīng)力與應(yīng)變成正比。由廣義虎克定律可知,在三向應(yīng)力狀態(tài)下,各方向的線應(yīng)變[17]分別為

    (7)

    (8)

    (9)

    式中:Er為橡膠的彈性模量;μr為橡膠的泊松比;x、y、z分別為膠圈微元段的寬度、厚度、長(zhǎng)度;Δx、Δy、Δz分別為膠圈微元段的寬度增量、厚度增量、長(zhǎng)度增量;σX、σY、σZ分別為X、Y、Z方向上的應(yīng)力。

    值得注意的是,由于膠圈對(duì)連接器的擠壓,連接器也會(huì)發(fā)生一定程度的變形。但連接器厚度較大且外側(cè)有加強(qiáng)筋,故強(qiáng)度和剛度都很大,可近似為剛性體,不產(chǎn)生形變。這樣,膠圈在Z方向無法發(fā)生伸長(zhǎng),即Δz=0。

    在壓力作用下,膠圈的體積應(yīng)變?yōu)?/p>

    (10)

    忽略二階微量,則式(10)變?yōu)?/p>

    (11)

    由式(9)及Δz=0可得:

    σZ=μr(σX+σY)

    (12)

    將式(12)和σX=σY=-p代入式(11)得:

    (13)

    該式即為膠圈受液體壓力后的體積壓縮量。

    2.2 考慮密封膠圈影響的水擊壓力波速

    先考慮一根管子和一個(gè)膠圈的情形,兩者長(zhǎng)度分別為L(zhǎng)p和Lr。管道末端突然關(guān)閥引發(fā)水擊。假定管道內(nèi)初始流速為U0,關(guān)閥后的最終流速為U1,則Δt=(Lp+Lr)/af時(shí)間內(nèi)流入(Lp+Lr)管段內(nèi)的流體體積為

    Vf=(U0-U1)(AfLp+ApLr)/af=

    ΔU(AfLp+ApLr)/af

    (14)

    式中:Af為管道內(nèi)截面積;Ap為管道外截面積。

    鋼管在流體壓力p作用下,其徑向膨脹量ΔD為

    (15)

    則單根鋼管的容積變化量ΔVp為

    (16)

    由于膠圈受到連接器約束,因此膠圈段被壓縮導(dǎo)致的體積減小量就等于其容積的增加量ΔV′r,即:

    (17)

    根據(jù)連續(xù)性原理,流入的流體量等于流體、膠圈的體積壓縮量和管道容積膨脹量之和:

    Vf=ΔVf+ΔVp+ΔV′r

    (18)

    式中:ΔVf為流體的體積壓縮量;ΔVf=p(AfLp+ApLr)/K。

    將式(14)(16)(17)代入式(18)可得:

    (19)

    將慣性水擊壓力表達(dá)式ΔU=p/(ρfaf)代入式(19),即可得計(jì)及密封膠圈影響的經(jīng)典水擊壓力波速計(jì)算式:

    (20)

    上述的推導(dǎo)是基于一根管子和一個(gè)膠圈的情形。若推廣到整條管線,設(shè)全線共有管子np根,密封膠圈nr個(gè),則水擊壓力波速計(jì)算式為

    af=

    (21)

    式中:α=(AfLp+βApLr);β=nr/np。

    將式(21)代入式(5),則可得到考慮密封膠圈影響的裝配式管線耦合水擊波速計(jì)算公式。

    3 驗(yàn)證與分析

    3.1 試驗(yàn)驗(yàn)證

    試驗(yàn)管道系統(tǒng)由室內(nèi)管道、室外環(huán)形管道、電動(dòng)泵、控制閥以及數(shù)據(jù)采集設(shè)備等組成,如圖4所示。

    圖4 試驗(yàn)管道系統(tǒng)示意

    因?yàn)槭軐?shí)驗(yàn)場(chǎng)地限制,試驗(yàn)管道較多地使用了短管和彎頭,單管段平均長(zhǎng)度按4 m計(jì)。各管段、管件之間的連接方式為卡箍式連接,膠圈密封。變頻電動(dòng)離心泵將清水從水池抽吸到試驗(yàn)管道中循環(huán)輸送。當(dāng)管內(nèi)流動(dòng)達(dá)到穩(wěn)定狀態(tài)時(shí),控制電磁閥突然關(guān)閉(或關(guān)閉一定時(shí)間后突然開啟)產(chǎn)生水擊工況。壓力傳感器量程為0~5.0 MPa,精度等級(jí)為0.2級(jí),頻率響應(yīng)為10 kHz。上、下游壓力傳感器之間的管道總長(zhǎng)度L為656 m,分別采用周期法與波前法測(cè)量水擊波速值。試驗(yàn)管道系統(tǒng)的主要參數(shù)見表1。

    表1 試驗(yàn)管道系統(tǒng)主要參數(shù)

    當(dāng)穩(wěn)態(tài)流量達(dá)27.9 m3/h時(shí),發(fā)出關(guān)閥信號(hào)使電磁閥瞬間關(guān)閉,獲得水擊壓力響應(yīng)波形如圖5所示。根據(jù)水擊周期法,若水擊壓力振蕩周期為T,則水擊壓力波速af=4L/T。經(jīng)3次重復(fù)試驗(yàn),測(cè)得試驗(yàn)數(shù)據(jù)如表2所示。

    表2 瞬間關(guān)閥的水擊試驗(yàn)數(shù)據(jù)(周期法)

    圖5 瞬間關(guān)閥的水擊壓力響應(yīng)曲線

    當(dāng)穩(wěn)態(tài)流量為30.1 m3/h時(shí),先將電磁閥關(guān)閉一段時(shí)間后再突然打開,記錄減壓波壓力波前變化曲線,如圖6所示。

    圖6 瞬間開閥的減壓波前曲線

    根據(jù)水擊波前法,若上、下游2個(gè)壓力傳感器檢測(cè)到壓力波前到達(dá)的時(shí)間差為Δt,則水擊壓力波速af=L/Δt。同樣進(jìn)行3次水擊試驗(yàn),試驗(yàn)數(shù)據(jù)如表3所示。

    表3 瞬間開閥的水擊試驗(yàn)數(shù)據(jù)(波前法)

    按表1中的參數(shù)值,分別采用文中提出的方法和Korteweg式進(jìn)行水擊壓力波速計(jì)算,并與試驗(yàn)實(shí)測(cè)值進(jìn)行對(duì)比,如表4所示。

    表4 水擊波速計(jì)算值與試驗(yàn)實(shí)測(cè)值對(duì)比

    由表4可知:水擊波速實(shí)測(cè)平均值在1 040~1 065 m/s之間,明顯小于Korteweg式的經(jīng)典水擊波速計(jì)算值。這說明裝配式管線的連接結(jié)構(gòu)對(duì)水擊響應(yīng)具有一定的阻尼效果,密封膠圈的存在對(duì)水擊波速起到了明顯的衰減作用。文中方法計(jì)算值與試驗(yàn)實(shí)測(cè)值吻合得較好,說明文中方法是正確可行的。相比不考慮FSI效應(yīng)而言,考慮FSI效應(yīng)的耦合水擊波速計(jì)算值更接近試驗(yàn)實(shí)測(cè)值,表明管道結(jié)構(gòu)運(yùn)動(dòng)和變形對(duì)流體壓力波動(dòng)響應(yīng)的影響是不容忽視的。

    3.2 密封膠圈對(duì)不同輸送介質(zhì)中水擊波速的影響

    當(dāng)裝配式管線分別輸送清水、汽油和柴油時(shí),采用文中方法進(jìn)行水擊波速計(jì)算。管道參數(shù)見表1,汽油密度為730 kg/m3,體積模量為916 MPa;柴油密度為840 kg/m3,體積模量為1 410 MPa。根據(jù)上述參數(shù),不同輸送介質(zhì)中的水擊波速計(jì)算值如表5所示。

    表5 不同輸送介質(zhì)中的水擊波速計(jì)算值

    表5中β=1.0是裝配式管線系統(tǒng)中短管件數(shù)量少、可忽略不計(jì)的情況,若不可忽略時(shí),則β>1.0。從表5可知:密封膠圈對(duì)水擊波速的衰減程度隨著其數(shù)量的增多而增大。在不同輸送介質(zhì)中,考慮FSI效應(yīng)的耦合水擊波速值均小于經(jīng)典水擊波速值,這是因?yàn)樵诹鞴恬詈献饔孟?,流體的部分水擊壓能轉(zhuǎn)換成了管道結(jié)構(gòu)的振動(dòng)能量,導(dǎo)致流體水擊波速有所降低。

    圖7所示為密封膠圈對(duì)不同介質(zhì)中水擊波速衰減率的變化趨勢(shì)對(duì)比曲線。圖中的波速衰減率是指不計(jì)密封膠圈影響的波速與計(jì)及密封膠圈影響的波速之差的百分?jǐn)?shù)。

    圖7 密封膠圈對(duì)水擊波速的衰減率曲線

    從圖7可以看出:當(dāng)考慮FSI效應(yīng)時(shí),密封膠圈對(duì)水擊波速的衰減作用均比不考慮FSI效應(yīng)時(shí)小。對(duì)3種不同輸送介質(zhì)而言,密封膠圈對(duì)水中的水擊波速衰減率最大,柴油次之,汽油最小。這主要是因?yàn)樗捏w積模量最大,最不容易被壓縮,此時(shí)密封膠圈的彈性則體現(xiàn)得最明顯,因而對(duì)水中水擊波速的衰減作用也最顯著。

    3.3 管道不同厚徑比對(duì)水擊波速的影響

    分析裝配式管線不同厚徑比e/D對(duì)水擊波速的影響,輸送介質(zhì)為清水。為了顯示文中方法與薄壁管水擊波速模型的區(qū)別,對(duì)式(21)按薄壁管假設(shè)進(jìn)行簡(jiǎn)化,可得計(jì)及膠圈影響的薄壁管經(jīng)典水擊波速表達(dá)式:

    (22)

    將式(22)代入式(5)則得到計(jì)及密封膠圈影響的薄壁管耦合水擊波速計(jì)算式。

    圖8展示了在不同密封膠圈數(shù)量下,耦合水擊波速隨管道厚徑比的變化情況。

    圖8 管道厚徑比對(duì)耦合水擊波速的影響

    從圖8可以看出:在相同厚徑比條件下,采用文中方法的水擊波速計(jì)算值均小于薄壁管模型計(jì)算值,這與文獻(xiàn)[17]中的相關(guān)結(jié)論一致。薄壁管模型假設(shè)管壁軸向應(yīng)力沿徑向?yàn)槌?shù),而文中方法基于厚壁管模型,考慮了管壁軸向應(yīng)力沿徑向變化的實(shí)際情況,采用的管壁截面軸向應(yīng)力平均值小于管內(nèi)壁與流體交界面處的軸向應(yīng)力,因此計(jì)算得出的耦合水擊波速值也相對(duì)較小。當(dāng)e/D<0.05時(shí),文中方法與薄壁管模型的波速計(jì)算值相差小于5%,此時(shí)采用文中方法或薄壁管模型均是可行的;當(dāng)e/D≥0.05時(shí),采用文中方法進(jìn)行計(jì)算則更為可靠。

    (23)

    這是由于裝配式管線采用連接器連接,各管段接頭之間存在軸向間隙。每個(gè)間隙與外側(cè)的密封膠圈之間均形成了環(huán)狀空腔結(jié)構(gòu),可吸收水擊能量,在水擊過程中擔(dān)任著“水擊壓力緩沖環(huán)”的角色,能對(duì)水擊波速起到明顯的衰減作用。當(dāng)e/D小于某值時(shí),管道尚屬薄壁管范疇,雖然e/D增加,但其環(huán)狀空腔體積增量較小,不足以抵消由于管壁增厚、管道徑向膨脹量減小對(duì)水擊波速帶來的增速效果,因此水擊波速呈暫時(shí)增加的趨勢(shì)。當(dāng)e/D超過某值時(shí),隨著e/D增加,環(huán)狀空腔體積增量對(duì)水擊響應(yīng)的衰減作用逐漸占據(jù)主導(dǎo)地位,從而導(dǎo)致水擊波速持續(xù)降低直至極限值。

    4 結(jié)論

    裝配式管線中含有大量密封膠圈,對(duì)全線水擊響應(yīng)的影響顯著。文中在考慮FSI效應(yīng)的基礎(chǔ)上,分析了密封膠圈在水擊壓力下的形變情況,提出了裝配式管線耦合水擊壓力波速的計(jì)算方法。采用水擊試驗(yàn)系統(tǒng)對(duì)文中方法進(jìn)行了試驗(yàn)驗(yàn)證,并分析了密封膠圈和管道厚徑比對(duì)水擊波速的影響,得出如下結(jié)論:

    (1)試驗(yàn)實(shí)測(cè)值與文中方法計(jì)算值吻合較好,表明文中方法適用于裝配式管線水擊波速的預(yù)測(cè)計(jì)算與分析。其中,考慮FSI效應(yīng)的水擊波速計(jì)算值更接近實(shí)測(cè)值,說明對(duì)裝配式管線系統(tǒng)的水擊響應(yīng)分析應(yīng)考慮FSI效應(yīng)的影響。

    (2)密封膠圈對(duì)水擊波速起衰減作用,衰減率隨其數(shù)量的增多而增大。當(dāng)考慮FSI效應(yīng)時(shí),密封膠圈對(duì)水擊波速的衰減作用比不考慮FSI效應(yīng)時(shí)小。管道內(nèi)輸送介質(zhì)的體積模量越大,密封膠圈對(duì)水擊壓力波速的衰減作用越明顯。

    (3)文中方法基于厚壁管模型,當(dāng)e/D<0.05時(shí),文中方法的水擊波速計(jì)算值與薄壁管模型計(jì)算值相對(duì)偏差較?。划?dāng)e/D≥0.05時(shí),采用文中方法進(jìn)行計(jì)算則具有更高的可靠性。裝配式管線各管段之間的連接部位形成了具有吸收水擊能量效果的環(huán)狀空腔結(jié)構(gòu)。隨著e/D的增加,環(huán)狀空腔對(duì)水擊響應(yīng)的衰減作用逐漸凸顯出來,因此管內(nèi)水擊波速呈現(xiàn)短暫增大后又持續(xù)減小的變化趨勢(shì)。

    猜你喜歡
    水擊膠圈密封膠
    溫暖講成語:鵬程萬里
    小讀者之友(2024年5期)2024-06-01 16:10:03
    玻璃生產(chǎn)線膠圈輥加工工藝中的控制要點(diǎn)
    建材世界(2022年6期)2022-12-10 11:53:56
    膠圈控制機(jī)構(gòu)探究
    水擊壓力簡(jiǎn)化計(jì)算在有壓輸水管道設(shè)計(jì)中的應(yīng)用
    陜西水利(2021年1期)2021-04-12 02:25:36
    建筑用硅酮密封膠及其市場(chǎng)觀察
    上海建材(2020年5期)2020-12-15 00:33:26
    新型圓珠筆筆頭密封膠應(yīng)用研究
    WRA-S-25型菱形膠圈的研發(fā)及生產(chǎn)實(shí)踐
    紡織器材(2017年6期)2017-12-14 06:38:02
    海底輸油管道水力瞬變分析模型探究及水擊保護(hù)
    減小硬關(guān)井水擊問題探討
    化工管理(2015年17期)2015-12-22 02:44:58
    “Ω”形膠圈密封機(jī)理研究
    亚洲av不卡在线观看| 亚洲内射少妇av| 国产 一区精品| 九九在线视频观看精品| 国产精品99久久久久久久久| 亚洲成人一二三区av| 亚洲国产最新在线播放| 欧美另类一区| 亚洲内射少妇av| 一区二区av电影网| 成人手机av| 亚洲精品aⅴ在线观看| 国产在线视频一区二区| 97超碰精品成人国产| 午夜福利视频在线观看免费| 国产欧美日韩一区二区三区在线 | 伦理电影大哥的女人| 日韩在线高清观看一区二区三区| 一级爰片在线观看| 18禁观看日本| 欧美 日韩 精品 国产| 丰满饥渴人妻一区二区三| 亚洲欧美中文字幕日韩二区| 亚洲综合精品二区| 中文字幕av电影在线播放| 激情五月婷婷亚洲| 国产精品三级大全| 美女国产视频在线观看| 免费播放大片免费观看视频在线观看| 成人毛片60女人毛片免费| 99视频精品全部免费 在线| 国产在线视频一区二区| 欧美日本中文国产一区发布| 久热这里只有精品99| 亚洲欧美色中文字幕在线| 18禁观看日本| 欧美成人午夜免费资源| 美女视频免费永久观看网站| 一级黄片播放器| 最近2019中文字幕mv第一页| 亚洲精品国产av成人精品| 美女福利国产在线| 国产精品女同一区二区软件| 一级a做视频免费观看| xxx大片免费视频| 亚洲综合色惰| 亚洲av成人精品一二三区| 免费观看a级毛片全部| 久久99蜜桃精品久久| 日韩av不卡免费在线播放| 成年美女黄网站色视频大全免费 | 国产日韩一区二区三区精品不卡 | 亚洲国产最新在线播放| 成人无遮挡网站| 91精品伊人久久大香线蕉| 久久97久久精品| 亚洲精品久久午夜乱码| 日本爱情动作片www.在线观看| 成人18禁高潮啪啪吃奶动态图 | 能在线免费看毛片的网站| 精品卡一卡二卡四卡免费| 精品国产国语对白av| 国产免费一区二区三区四区乱码| 91在线精品国自产拍蜜月| 免费观看性生交大片5| 久久综合国产亚洲精品| 国产高清国产精品国产三级| 亚洲精品色激情综合| 九九久久精品国产亚洲av麻豆| 久久久久精品久久久久真实原创| 亚洲天堂av无毛| 国产精品国产三级专区第一集| 满18在线观看网站| 精品熟女少妇av免费看| 色网站视频免费| 国产成人精品在线电影| 九草在线视频观看| 国产欧美日韩一区二区三区在线 | 欧美+日韩+精品| 欧美一级a爱片免费观看看| 老女人水多毛片| 十八禁高潮呻吟视频| 国产乱人偷精品视频| av线在线观看网站| 国产亚洲精品第一综合不卡 | 如何舔出高潮| 国产精品一二三区在线看| 晚上一个人看的免费电影| 美女视频免费永久观看网站| 如何舔出高潮| 香蕉精品网在线| 蜜桃久久精品国产亚洲av| 熟妇人妻不卡中文字幕| 妹子高潮喷水视频| 国产欧美日韩一区二区三区在线 | 欧美日韩精品成人综合77777| 一区在线观看完整版| 国产 一区精品| 少妇丰满av| 精品久久久久久电影网| 国产欧美另类精品又又久久亚洲欧美| 五月玫瑰六月丁香| 高清av免费在线| 国产视频首页在线观看| 在线观看国产h片| 日韩强制内射视频| 在线观看美女被高潮喷水网站| 天堂8中文在线网| 日本爱情动作片www.在线观看| 91精品伊人久久大香线蕉| 一区二区三区乱码不卡18| 午夜福利影视在线免费观看| 在线 av 中文字幕| 九色亚洲精品在线播放| 久久久精品区二区三区| 亚洲国产欧美日韩在线播放| 女的被弄到高潮叫床怎么办| 十八禁高潮呻吟视频| 国产免费又黄又爽又色| 少妇的逼水好多| 欧美人与善性xxx| 国产精品一区二区三区四区免费观看| 人体艺术视频欧美日本| 成人毛片60女人毛片免费| 日韩av在线免费看完整版不卡| 蜜臀久久99精品久久宅男| 国产午夜精品久久久久久一区二区三区| 亚洲色图 男人天堂 中文字幕 | a级片在线免费高清观看视频| 99国产精品免费福利视频| 人人妻人人澡人人爽人人夜夜| 精品熟女少妇av免费看| 久久免费观看电影| 久久99蜜桃精品久久| 亚洲av成人精品一二三区| 一本—道久久a久久精品蜜桃钙片| 日韩av免费高清视频| 成人影院久久| 下体分泌物呈黄色| 在线观看www视频免费| 欧美日韩亚洲高清精品| 蜜桃久久精品国产亚洲av| 一个人免费看片子| av视频免费观看在线观看| 亚洲性久久影院| 国产成人精品福利久久| 夜夜看夜夜爽夜夜摸| 永久免费av网站大全| 国产黄色免费在线视频| 观看av在线不卡| 精品久久蜜臀av无| 亚洲精品国产色婷婷电影| 久久国产亚洲av麻豆专区| 中文字幕精品免费在线观看视频 | 99久久综合免费| 一本一本综合久久| 人妻系列 视频| 美女视频免费永久观看网站| 亚洲少妇的诱惑av| 少妇精品久久久久久久| 精品视频人人做人人爽| 久久av网站| 你懂的网址亚洲精品在线观看| 91久久精品电影网| 男人操女人黄网站| 菩萨蛮人人尽说江南好唐韦庄| 777米奇影视久久| 久久影院123| 男男h啪啪无遮挡| 国产老妇伦熟女老妇高清| 久热久热在线精品观看| 国产乱人偷精品视频| 亚洲丝袜综合中文字幕| 在线精品无人区一区二区三| 日本黄大片高清| 街头女战士在线观看网站| 成人国语在线视频| 久久综合国产亚洲精品| 亚洲国产最新在线播放| 水蜜桃什么品种好| 女的被弄到高潮叫床怎么办| 成人综合一区亚洲| 性高湖久久久久久久久免费观看| 永久免费av网站大全| 国产精品无大码| 成人国语在线视频| 国产成人av激情在线播放 | 少妇 在线观看| 久久久久久久国产电影| 青春草视频在线免费观看| 人妻一区二区av| 国产高清不卡午夜福利| 午夜日本视频在线| 亚洲欧洲国产日韩| 国产精品免费大片| 最近中文字幕高清免费大全6| 亚洲av免费高清在线观看| 婷婷成人精品国产| 日日啪夜夜爽| 日韩欧美一区视频在线观看| 99国产精品免费福利视频| 一级毛片黄色毛片免费观看视频| 久久热精品热| 欧美人与性动交α欧美精品济南到 | 两个人免费观看高清视频| 一级毛片aaaaaa免费看小| 亚洲成色77777| 这个男人来自地球电影免费观看 | 日韩精品有码人妻一区| 亚洲内射少妇av| 丰满少妇做爰视频| 91精品国产九色| 青春草视频在线免费观看| 少妇高潮的动态图| av网站免费在线观看视频| 国产极品粉嫩免费观看在线 | 最后的刺客免费高清国语| 国产乱来视频区| 在线观看人妻少妇| 99国产精品免费福利视频| 天天躁夜夜躁狠狠久久av| 国产一区二区在线观看av| 一级二级三级毛片免费看| 水蜜桃什么品种好| 亚洲人成网站在线观看播放| 亚洲精品一二三| 免费观看av网站的网址| 91aial.com中文字幕在线观看| 中国三级夫妇交换| 亚洲欧美日韩卡通动漫| 男人爽女人下面视频在线观看| av线在线观看网站| 欧美激情国产日韩精品一区| .国产精品久久| 18在线观看网站| 18禁裸乳无遮挡动漫免费视频| 欧美日韩视频高清一区二区三区二| 国产亚洲精品久久久com| 国产精品蜜桃在线观看| 欧美一级a爱片免费观看看| 波野结衣二区三区在线| 久久99精品国语久久久| 十分钟在线观看高清视频www| 国产在线一区二区三区精| 蜜桃在线观看..| 国产精品久久久久久久电影| 男人操女人黄网站| av一本久久久久| 国模一区二区三区四区视频| 蜜桃国产av成人99| 欧美一级a爱片免费观看看| 免费看光身美女| 国产精品 国内视频| 久久精品国产自在天天线| 大片电影免费在线观看免费| 色视频在线一区二区三区| 美女国产高潮福利片在线看| 久热久热在线精品观看| 日本欧美视频一区| 看十八女毛片水多多多| 如日韩欧美国产精品一区二区三区 | 少妇人妻 视频| 午夜久久久在线观看| 亚洲av电影在线观看一区二区三区| 亚洲国产精品成人久久小说| 中文字幕亚洲精品专区| 久久毛片免费看一区二区三区| .国产精品久久| 搡女人真爽免费视频火全软件| 日韩av免费高清视频| 国模一区二区三区四区视频| 中文字幕制服av| 亚洲无线观看免费| 91久久精品电影网| 国产精品人妻久久久影院| 人妻人人澡人人爽人人| 久久久久久久大尺度免费视频| 免费少妇av软件| 黑丝袜美女国产一区| 又大又黄又爽视频免费| 一区二区三区四区激情视频| 最近2019中文字幕mv第一页| 大又大粗又爽又黄少妇毛片口| 色网站视频免费| 91aial.com中文字幕在线观看| 多毛熟女@视频| 热99国产精品久久久久久7| av线在线观看网站| 黄色配什么色好看| 亚洲av中文av极速乱| 亚洲欧洲日产国产| 久久久久精品性色| 久久久国产欧美日韩av| 一级毛片黄色毛片免费观看视频| 日韩制服骚丝袜av| 免费观看性生交大片5| 欧美丝袜亚洲另类| 国产不卡av网站在线观看| 日韩人妻高清精品专区| 777米奇影视久久| 大香蕉97超碰在线| 成人综合一区亚洲| 亚洲av日韩在线播放| 国产精品欧美亚洲77777| 日韩免费高清中文字幕av| 国产爽快片一区二区三区| 精品少妇黑人巨大在线播放| 日本午夜av视频| 午夜免费鲁丝| 一边摸一边做爽爽视频免费| 国产不卡av网站在线观看| 欧美亚洲 丝袜 人妻 在线| 成人18禁高潮啪啪吃奶动态图 | 久久婷婷青草| 天天操日日干夜夜撸| 多毛熟女@视频| 在线免费观看不下载黄p国产| 久久这里有精品视频免费| 我的老师免费观看完整版| 国产成人免费观看mmmm| 嫩草影院入口| 99久久人妻综合| 亚洲av二区三区四区| 久久精品国产自在天天线| 色视频在线一区二区三区| 亚洲五月色婷婷综合| 免费人成在线观看视频色| 在线看a的网站| 成人国产av品久久久| 少妇 在线观看| 亚洲国产毛片av蜜桃av| 夫妻性生交免费视频一级片| 亚洲天堂av无毛| 国精品久久久久久国模美| 国产一级毛片在线| 伊人久久精品亚洲午夜| 观看av在线不卡| av在线播放精品| 亚洲美女搞黄在线观看| 大香蕉久久网| 日韩不卡一区二区三区视频在线| 国精品久久久久久国模美| 如何舔出高潮| 99国产精品免费福利视频| 精品国产乱码久久久久久小说| 啦啦啦中文免费视频观看日本| a级毛色黄片| 久久精品国产亚洲av涩爱| 欧美日韩精品成人综合77777| 国产av国产精品国产| 欧美3d第一页| 免费看光身美女| 久久99蜜桃精品久久| xxxhd国产人妻xxx| 久久精品国产亚洲av天美| 伦理电影大哥的女人| 精品久久久久久电影网| 亚洲精品日本国产第一区| 五月天丁香电影| 成人午夜精彩视频在线观看| 中文欧美无线码| 亚洲激情五月婷婷啪啪| 在线观看国产h片| 美女视频免费永久观看网站| 国产欧美亚洲国产| 97超视频在线观看视频| 国产精品熟女久久久久浪| 国产av一区二区精品久久| 久久久亚洲精品成人影院| 一级毛片电影观看| 在线观看www视频免费| 最近中文字幕2019免费版| 国产精品国产三级国产av玫瑰| 色婷婷久久久亚洲欧美| 欧美精品高潮呻吟av久久| 精品熟女少妇av免费看| 久久久精品94久久精品| 久久精品国产亚洲网站| 欧美日本中文国产一区发布| 精品国产乱码久久久久久小说| 亚洲欧美清纯卡通| 丝袜喷水一区| a级毛片在线看网站| 国产av国产精品国产| av在线app专区| 天堂中文最新版在线下载| 免费av中文字幕在线| 一边摸一边做爽爽视频免费| 国产毛片在线视频| 观看av在线不卡| 简卡轻食公司| 97超碰精品成人国产| 国产 一区精品| 精品久久久精品久久久| 在线观看免费视频网站a站| 亚洲精品,欧美精品| 日韩成人av中文字幕在线观看| 久久久久久久大尺度免费视频| 人人妻人人澡人人爽人人夜夜| 自拍欧美九色日韩亚洲蝌蚪91| 午夜激情av网站| 亚洲无线观看免费| 亚洲伊人久久精品综合| 最后的刺客免费高清国语| 人妻 亚洲 视频| 久久久久精品久久久久真实原创| 午夜老司机福利剧场| 欧美一级a爱片免费观看看| 亚洲精品自拍成人| 免费大片18禁| 黑人巨大精品欧美一区二区蜜桃 | 少妇高潮的动态图| 国产国拍精品亚洲av在线观看| 亚洲欧美一区二区三区国产| 97超碰精品成人国产| 欧美人与性动交α欧美精品济南到 | 岛国毛片在线播放| 日本vs欧美在线观看视频| 久久久久网色| 久久狼人影院| 亚洲精品亚洲一区二区| 亚洲精品久久午夜乱码| 在线亚洲精品国产二区图片欧美 | 亚洲精品日韩在线中文字幕| 亚洲国产毛片av蜜桃av| 国产av码专区亚洲av| 一边亲一边摸免费视频| 免费观看的影片在线观看| 在线观看人妻少妇| 人成视频在线观看免费观看| 九色成人免费人妻av| 在线观看www视频免费| 久久久久久久久久成人| 日本vs欧美在线观看视频| 最后的刺客免费高清国语| 赤兔流量卡办理| 日本与韩国留学比较| 成人影院久久| 精品人妻一区二区三区麻豆| 夜夜看夜夜爽夜夜摸| 欧美性感艳星| 最近2019中文字幕mv第一页| 免费大片黄手机在线观看| 人人妻人人澡人人看| 男人爽女人下面视频在线观看| 一级毛片aaaaaa免费看小| 女人久久www免费人成看片| 国产精品欧美亚洲77777| 亚洲精品日本国产第一区| 欧美激情 高清一区二区三区| 91aial.com中文字幕在线观看| av电影中文网址| 又黄又爽又刺激的免费视频.| 国产极品天堂在线| 亚洲精品国产色婷婷电影| 91精品国产九色| 久久鲁丝午夜福利片| 中文欧美无线码| 午夜福利影视在线免费观看| 亚洲国产毛片av蜜桃av| 久久国产精品大桥未久av| 欧美少妇被猛烈插入视频| 三级国产精品片| 成人亚洲精品一区在线观看| 亚洲伊人久久精品综合| 日本黄色日本黄色录像| 亚洲精品aⅴ在线观看| av又黄又爽大尺度在线免费看| 成人无遮挡网站| 一区二区日韩欧美中文字幕 | 大香蕉97超碰在线| 美女中出高潮动态图| 只有这里有精品99| 亚洲精品久久午夜乱码| 韩国av在线不卡| 夫妻性生交免费视频一级片| 日日撸夜夜添| 99九九在线精品视频| av黄色大香蕉| 你懂的网址亚洲精品在线观看| √禁漫天堂资源中文www| 韩国高清视频一区二区三区| 夜夜看夜夜爽夜夜摸| 亚洲av.av天堂| 国产亚洲一区二区精品| h视频一区二区三区| 亚洲怡红院男人天堂| 国产一级毛片在线| 国产成人精品在线电影| 日韩大片免费观看网站| 中文乱码字字幕精品一区二区三区| 久久久精品免费免费高清| 26uuu在线亚洲综合色| 国内精品宾馆在线| 成年人免费黄色播放视频| 在线观看三级黄色| 一本一本综合久久| 亚洲成色77777| 久久久国产一区二区| 久久精品国产亚洲av天美| 天堂中文最新版在线下载| 久久青草综合色| 极品少妇高潮喷水抽搐| 久久99精品国语久久久| 国产精品久久久久久av不卡| 亚洲国产毛片av蜜桃av| 亚洲av成人精品一二三区| 国产片特级美女逼逼视频| 少妇高潮的动态图| 大码成人一级视频| 久久久久久久国产电影| av.在线天堂| 自拍欧美九色日韩亚洲蝌蚪91| 在线亚洲精品国产二区图片欧美 | 久久99热这里只频精品6学生| 久久精品国产鲁丝片午夜精品| 免费人妻精品一区二区三区视频| 国模一区二区三区四区视频| 精品一区二区三卡| 搡女人真爽免费视频火全软件| 老司机亚洲免费影院| 久久久久久久久久久免费av| 亚洲精品国产色婷婷电影| 丝袜喷水一区| 69精品国产乱码久久久| 麻豆精品久久久久久蜜桃| 国产精品一区二区三区四区免费观看| 国产精品欧美亚洲77777| 在线亚洲精品国产二区图片欧美 | 久久午夜综合久久蜜桃| 蜜桃久久精品国产亚洲av| 亚洲第一区二区三区不卡| 美女国产视频在线观看| 亚洲国产最新在线播放| 一级爰片在线观看| 亚洲精品自拍成人| 久久国内精品自在自线图片| 蜜桃在线观看..| 97超视频在线观看视频| 国产日韩欧美视频二区| 久久国产精品大桥未久av| 又粗又硬又长又爽又黄的视频| 特大巨黑吊av在线直播| 肉色欧美久久久久久久蜜桃| 国产极品粉嫩免费观看在线 | 熟女人妻精品中文字幕| 一个人看视频在线观看www免费| 免费播放大片免费观看视频在线观看| 美女cb高潮喷水在线观看| 久久女婷五月综合色啪小说| 大片电影免费在线观看免费| 亚洲人成网站在线播| 永久免费av网站大全| 美女内射精品一级片tv| 考比视频在线观看| 九色成人免费人妻av| 中文字幕亚洲精品专区| 边亲边吃奶的免费视频| 国产成人a∨麻豆精品| 亚洲国产精品999| 十分钟在线观看高清视频www| 嘟嘟电影网在线观看| 精品国产露脸久久av麻豆| 亚洲一级一片aⅴ在线观看| 国产精品久久久久久久久免| 国产免费现黄频在线看| 伦理电影免费视频| 久久精品国产a三级三级三级| 欧美亚洲 丝袜 人妻 在线| 丰满少妇做爰视频| 看十八女毛片水多多多| 男女边吃奶边做爰视频| 亚洲第一av免费看| 街头女战士在线观看网站| 久久久久国产精品人妻一区二区| 色94色欧美一区二区| 一本色道久久久久久精品综合| 麻豆精品久久久久久蜜桃| 亚洲一区二区三区欧美精品| 久久韩国三级中文字幕| 在线观看一区二区三区激情| 色吧在线观看| 韩国av在线不卡| 美女脱内裤让男人舔精品视频| 精品午夜福利在线看| 丁香六月天网| 免费高清在线观看视频在线观看| 熟女人妻精品中文字幕| 成人毛片60女人毛片免费| 久久97久久精品| tube8黄色片| 午夜av观看不卡| 久久精品国产亚洲网站| 精品国产露脸久久av麻豆| 大片免费播放器 马上看| 赤兔流量卡办理| 色婷婷av一区二区三区视频| 韩国av在线不卡| 欧美xxⅹ黑人| 午夜福利影视在线免费观看| 国产成人freesex在线| 黄色怎么调成土黄色| 久热这里只有精品99| 日韩av免费高清视频| 亚洲精品乱码久久久v下载方式| 最新中文字幕久久久久| 国产精品国产三级国产专区5o| 日韩av不卡免费在线播放| 午夜91福利影院| 成人亚洲精品一区在线观看| 狂野欧美激情性xxxx在线观看| 精品一品国产午夜福利视频| 成人国语在线视频| 国产高清三级在线| 国产精品 国内视频|