• <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ī)理研究
    一级,二级,三级黄色视频| 99久久人妻综合| 蜜桃在线观看..| 99国产精品一区二区蜜桃av | 一级a爱视频在线免费观看| 日韩制服丝袜自拍偷拍| 美女高潮到喷水免费观看| 欧美日韩成人在线一区二区| 丰满人妻熟妇乱又伦精品不卡| av电影中文网址| 欧美日韩亚洲高清精品| 老汉色av国产亚洲站长工具| 在线观看免费日韩欧美大片| 在线观看国产h片| 亚洲欧美中文字幕日韩二区| 考比视频在线观看| 日韩大码丰满熟妇| 91老司机精品| 国产精品一二三区在线看| 两个人免费观看高清视频| 亚洲欧美一区二区三区黑人| 国产成人啪精品午夜网站| 热99国产精品久久久久久7| av天堂在线播放| 国产亚洲欧美在线一区二区| 国产一区二区在线观看av| 丁香六月欧美| av在线app专区| 亚洲成色77777| 在线观看免费午夜福利视频| 国产av国产精品国产| 999精品在线视频| 亚洲,欧美精品.| 免费看十八禁软件| 国产人伦9x9x在线观看| 老汉色∧v一级毛片| 久久av网站| 国产成人影院久久av| 男女高潮啪啪啪动态图| 天堂俺去俺来也www色官网| 国产亚洲欧美精品永久| 美女国产高潮福利片在线看| 热99久久久久精品小说推荐| 美女大奶头黄色视频| 一边亲一边摸免费视频| 下体分泌物呈黄色| 欧美激情 高清一区二区三区| 男女边摸边吃奶| 男女边吃奶边做爰视频| 国产精品免费视频内射| 中文字幕人妻丝袜一区二区| 日本午夜av视频| xxxhd国产人妻xxx| 亚洲av在线观看美女高潮| 视频区欧美日本亚洲| 日本一区二区免费在线视频| 亚洲中文日韩欧美视频| 自线自在国产av| 成人免费观看视频高清| 亚洲精品乱久久久久久| 一级a爱视频在线免费观看| 精品一品国产午夜福利视频| 免费高清在线观看视频在线观看| 精品少妇久久久久久888优播| 久久久久久免费高清国产稀缺| 波多野结衣av一区二区av| 日本91视频免费播放| 天天躁日日躁夜夜躁夜夜| 午夜激情av网站| 女人被躁到高潮嗷嗷叫费观| av天堂在线播放| 男女床上黄色一级片免费看| 欧美在线一区亚洲| 美女主播在线视频| 亚洲情色 制服丝袜| 少妇 在线观看| 欧美成人精品欧美一级黄| 国产日韩欧美在线精品| 尾随美女入室| 日本a在线网址| 一区二区三区精品91| 视频区图区小说| 亚洲欧洲日产国产| 黄色怎么调成土黄色| 夜夜骑夜夜射夜夜干| 国产精品国产三级国产专区5o| 伊人久久大香线蕉亚洲五| 国产在线观看jvid| netflix在线观看网站| 久久亚洲国产成人精品v| 女性生殖器流出的白浆| 久久99一区二区三区| 精品免费久久久久久久清纯 | 久久精品熟女亚洲av麻豆精品| 亚洲第一青青草原| 久久久精品免费免费高清| videosex国产| 九草在线视频观看| 少妇人妻 视频| 丰满人妻熟妇乱又伦精品不卡| 亚洲欧美日韩高清在线视频 | 久久99精品国语久久久| 2018国产大陆天天弄谢| 啦啦啦在线免费观看视频4| 黄色视频在线播放观看不卡| 国产视频一区二区在线看| 精品一品国产午夜福利视频| 狠狠婷婷综合久久久久久88av| 国产麻豆69| 欧美另类一区| av国产精品久久久久影院| 欧美性长视频在线观看| 一区二区三区激情视频| 精品亚洲成国产av| 午夜福利,免费看| 肉色欧美久久久久久久蜜桃| 欧美成狂野欧美在线观看| 国产一区亚洲一区在线观看| 人人妻人人添人人爽欧美一区卜| av不卡在线播放| 免费看av在线观看网站| 国产在线观看jvid| av有码第一页| 大香蕉久久成人网| 日本猛色少妇xxxxx猛交久久| e午夜精品久久久久久久| 久久精品成人免费网站| 国产成人一区二区三区免费视频网站 | 亚洲人成网站在线观看播放| 日韩熟女老妇一区二区性免费视频| 伊人亚洲综合成人网| 国产成人一区二区三区免费视频网站 | 一级毛片电影观看| 超碰成人久久| 另类精品久久| 国产真人三级小视频在线观看| 亚洲,欧美精品.| 日本a在线网址| 丰满人妻熟妇乱又伦精品不卡| 久久精品久久久久久久性| 亚洲,欧美,日韩| 精品国产一区二区久久| 日韩 亚洲 欧美在线| 国产精品二区激情视频| 国产一区亚洲一区在线观看| 美女视频免费永久观看网站| 黄频高清免费视频| 精品一区二区三区av网在线观看 | 久久国产亚洲av麻豆专区| av在线app专区| 亚洲黑人精品在线| 黄色毛片三级朝国网站| 午夜免费观看性视频| 中文欧美无线码| 好男人视频免费观看在线| 国产日韩欧美在线精品| 国产一区二区三区av在线| 精品高清国产在线一区| 欧美精品高潮呻吟av久久| 1024视频免费在线观看| 王馨瑶露胸无遮挡在线观看| 亚洲精品一二三| 国产极品粉嫩免费观看在线| 九色亚洲精品在线播放| 亚洲av电影在线进入| 亚洲综合色网址| 亚洲熟女毛片儿| 久久精品人人爽人人爽视色| 伦理电影免费视频| 国产亚洲精品第一综合不卡| 成人手机av| 色播在线永久视频| 久久久欧美国产精品| 99香蕉大伊视频| 欧美 亚洲 国产 日韩一| 午夜福利一区二区在线看| 69精品国产乱码久久久| 国产黄色视频一区二区在线观看| 国产精品一区二区在线观看99| 一级片'在线观看视频| 亚洲精品在线美女| 满18在线观看网站| 黄色毛片三级朝国网站| 国产不卡av网站在线观看| 在线观看免费午夜福利视频| 少妇 在线观看| 女性生殖器流出的白浆| 制服人妻中文乱码| 国产一区有黄有色的免费视频| 亚洲av国产av综合av卡| 91精品伊人久久大香线蕉| 国产精品香港三级国产av潘金莲 | 日本91视频免费播放| 这个男人来自地球电影免费观看| 亚洲欧美清纯卡通| 国产欧美亚洲国产| 日本猛色少妇xxxxx猛交久久| 女人被躁到高潮嗷嗷叫费观| 大型av网站在线播放| 久久99热这里只频精品6学生| 少妇精品久久久久久久| 日本午夜av视频| 欧美精品啪啪一区二区三区 | 99精国产麻豆久久婷婷| 人人妻人人添人人爽欧美一区卜| av欧美777| 操出白浆在线播放| 秋霞在线观看毛片| 欧美xxⅹ黑人| 精品少妇内射三级| 人妻一区二区av| e午夜精品久久久久久久| 精品少妇一区二区三区视频日本电影| 天堂中文最新版在线下载| 久久综合国产亚洲精品| 久久久精品免费免费高清| bbb黄色大片| 亚洲精品一卡2卡三卡4卡5卡 | cao死你这个sao货| 欧美日韩国产mv在线观看视频| 午夜免费男女啪啪视频观看| 自拍欧美九色日韩亚洲蝌蚪91| 飞空精品影院首页| 欧美精品人与动牲交sv欧美| 国产免费视频播放在线视频| 熟女av电影| 国产成人欧美在线观看 | 99热网站在线观看| 在线亚洲精品国产二区图片欧美| 国产精品国产三级专区第一集| 人妻人人澡人人爽人人| 精品人妻熟女毛片av久久网站| 精品久久蜜臀av无| 三上悠亚av全集在线观看| 国产成人精品久久二区二区91| 亚洲精品久久午夜乱码| 国产不卡av网站在线观看| 十分钟在线观看高清视频www| 韩国高清视频一区二区三区| 日韩av不卡免费在线播放| 大片免费播放器 马上看| 久久天堂一区二区三区四区| 国产成人91sexporn| 性色av一级| 亚洲欧美精品自产自拍| 精品国产国语对白av| 亚洲国产精品999| 可以免费在线观看a视频的电影网站| 免费在线观看影片大全网站 | 一本大道久久a久久精品| 亚洲成av片中文字幕在线观看| 又粗又硬又长又爽又黄的视频| 在线亚洲精品国产二区图片欧美| 中文欧美无线码| www日本在线高清视频| 亚洲国产日韩一区二区| 欧美久久黑人一区二区| 午夜免费观看性视频| 亚洲人成电影免费在线| 亚洲欧美日韩高清在线视频 | 麻豆国产av国片精品| 欧美日韩视频高清一区二区三区二| 麻豆乱淫一区二区| 欧美变态另类bdsm刘玥| 母亲3免费完整高清在线观看| 夫妻午夜视频| 狂野欧美激情性bbbbbb| 亚洲男人天堂网一区| 成年人黄色毛片网站| 校园人妻丝袜中文字幕| 蜜桃国产av成人99| 国产伦理片在线播放av一区| 波多野结衣一区麻豆| av在线播放精品| 热99久久久久精品小说推荐| 亚洲av成人精品一二三区| 国产亚洲av高清不卡| 天堂中文最新版在线下载| 午夜福利视频精品| 国产在线视频一区二区| 国产一区二区 视频在线| 国产日韩欧美在线精品| 首页视频小说图片口味搜索 | 亚洲综合色网址| 久久精品国产a三级三级三级| 欧美成狂野欧美在线观看| 777久久人妻少妇嫩草av网站| 亚洲九九香蕉| 午夜91福利影院| 欧美日本中文国产一区发布| 国产一区二区 视频在线| avwww免费| 两人在一起打扑克的视频| 晚上一个人看的免费电影| 汤姆久久久久久久影院中文字幕| 亚洲精品一卡2卡三卡4卡5卡 | 国产黄频视频在线观看| 一边摸一边做爽爽视频免费| 99国产精品一区二区蜜桃av | 亚洲av国产av综合av卡| 亚洲欧美精品综合一区二区三区| 精品福利永久在线观看| 欧美精品啪啪一区二区三区 | 亚洲精品自拍成人| 美女国产高潮福利片在线看| 国产福利在线免费观看视频| 免费在线观看视频国产中文字幕亚洲 | 亚洲国产精品国产精品| 最新的欧美精品一区二区| 脱女人内裤的视频| 午夜福利,免费看| 一级毛片电影观看| 久久毛片免费看一区二区三区| 亚洲综合色网址| 在线观看人妻少妇| 午夜老司机福利片| 99久久人妻综合| 欧美日韩一级在线毛片| 纯流量卡能插随身wifi吗| 久久精品亚洲av国产电影网| 亚洲成人免费电影在线观看 | 国产一级毛片在线| 免费观看a级毛片全部| 亚洲欧美激情在线| 国产黄色免费在线视频| 亚洲欧美一区二区三区国产| 在线看a的网站| 一本久久精品| 久久人妻熟女aⅴ| 天天躁日日躁夜夜躁夜夜| 国产免费现黄频在线看| 午夜福利视频精品| 精品久久久久久久毛片微露脸 | 色播在线永久视频| 亚洲av电影在线进入| 高清不卡的av网站| 91精品三级在线观看| 亚洲精品国产色婷婷电影| 精品视频人人做人人爽| 国产亚洲午夜精品一区二区久久| 久热这里只有精品99| 午夜激情久久久久久久| videos熟女内射| 黄色片一级片一级黄色片| 婷婷成人精品国产| 亚洲熟女精品中文字幕| 欧美成人精品欧美一级黄| 国产精品偷伦视频观看了| 老司机在亚洲福利影院| 亚洲国产欧美一区二区综合| 亚洲五月色婷婷综合| 国产成人精品无人区| 国产免费一区二区三区四区乱码| 大陆偷拍与自拍| 精品国产国语对白av| 90打野战视频偷拍视频| 99九九在线精品视频| 成人亚洲精品一区在线观看| 亚洲国产欧美网| 久久久国产精品麻豆| 久久 成人 亚洲| 国产三级黄色录像| 精品国产国语对白av| 精品久久久久久久毛片微露脸 | a级毛片黄视频| 大型av网站在线播放| 久久久久网色| 久久久亚洲精品成人影院| 欧美精品亚洲一区二区| 亚洲激情五月婷婷啪啪| 亚洲熟女毛片儿| 嫁个100分男人电影在线观看 | 久久这里只有精品19| 成人国产一区最新在线观看 | 一区福利在线观看| 国产亚洲午夜精品一区二区久久| 免费看不卡的av| 国产人伦9x9x在线观看| 99re6热这里在线精品视频| 久久人人爽av亚洲精品天堂| 国产在线视频一区二区| 丝袜美足系列| 国产亚洲精品久久久久5区| 日韩免费高清中文字幕av| 丝袜美腿诱惑在线| 亚洲自偷自拍图片 自拍| 亚洲欧美精品自产自拍| 下体分泌物呈黄色| 汤姆久久久久久久影院中文字幕| 国产精品一二三区在线看| 成人手机av| 免费在线观看影片大全网站 | 建设人人有责人人尽责人人享有的| 免费人妻精品一区二区三区视频| 丝袜在线中文字幕| 国产精品秋霞免费鲁丝片| 日韩制服丝袜自拍偷拍| 亚洲精品乱久久久久久| h视频一区二区三区| 国产精品欧美亚洲77777| 80岁老熟妇乱子伦牲交| 欧美另类一区| 亚洲第一青青草原| 国产av国产精品国产| 国产精品.久久久| 精品少妇内射三级| 天堂中文最新版在线下载| 中文乱码字字幕精品一区二区三区| cao死你这个sao货| 亚洲熟女毛片儿| 制服诱惑二区| 国产爽快片一区二区三区| 日韩制服骚丝袜av| 国产免费又黄又爽又色| 亚洲精品国产色婷婷电影| 亚洲av美国av| 久久 成人 亚洲| 免费高清在线观看日韩| 久久精品亚洲熟妇少妇任你| a级毛片在线看网站| 亚洲中文字幕日韩| 久久精品久久久久久噜噜老黄| 人人妻,人人澡人人爽秒播 | 日本五十路高清| 在线av久久热| 成人午夜精彩视频在线观看| 50天的宝宝边吃奶边哭怎么回事| 女人久久www免费人成看片| 亚洲欧美一区二区三区久久| 久久鲁丝午夜福利片| 在线观看国产h片| 男人添女人高潮全过程视频| 国产熟女午夜一区二区三区| 热99国产精品久久久久久7| 国产精品av久久久久免费| 亚洲欧美清纯卡通| 色精品久久人妻99蜜桃| 免费看av在线观看网站| 菩萨蛮人人尽说江南好唐韦庄| 亚洲成人免费av在线播放| 91麻豆av在线| 熟女av电影| 国精品久久久久久国模美| 成人亚洲欧美一区二区av| 在现免费观看毛片| 丝袜人妻中文字幕| 国产一区二区激情短视频 | 激情视频va一区二区三区| 久9热在线精品视频| 黑人巨大精品欧美一区二区蜜桃| 亚洲国产欧美在线一区| 国产不卡av网站在线观看| 免费在线观看日本一区| 国产熟女欧美一区二区| 久久国产精品影院| 别揉我奶头~嗯~啊~动态视频 | 99精国产麻豆久久婷婷| 老司机在亚洲福利影院| 一区在线观看完整版| 精品人妻1区二区| 亚洲国产欧美在线一区| 免费在线观看完整版高清| 丁香六月欧美| 男女免费视频国产| 亚洲成人免费电影在线观看 | 久久精品熟女亚洲av麻豆精品| 久久影院123| 黄色视频不卡| 久久精品亚洲熟妇少妇任你| 国产精品国产三级专区第一集| 中文字幕亚洲精品专区| 久久久久国产一级毛片高清牌| 久久久欧美国产精品| av在线播放精品| 免费观看人在逋| 亚洲专区国产一区二区| 99精品久久久久人妻精品| 在线观看人妻少妇| 中文字幕最新亚洲高清| 亚洲av日韩精品久久久久久密 | 国产男女超爽视频在线观看| 精品福利观看| 精品熟女少妇八av免费久了| 天天影视国产精品| 两个人看的免费小视频| 久久女婷五月综合色啪小说| 9热在线视频观看99| 香蕉国产在线看| 亚洲国产精品国产精品| av一本久久久久| 国产三级黄色录像| 男人爽女人下面视频在线观看| 久久性视频一级片| 女人高潮潮喷娇喘18禁视频| 丰满人妻熟妇乱又伦精品不卡| 高清视频免费观看一区二区| 天天躁日日躁夜夜躁夜夜| 满18在线观看网站| 国产精品久久久久久精品电影小说| 久久精品国产亚洲av高清一级| 9热在线视频观看99| 后天国语完整版免费观看| 国产高清视频在线播放一区 | 久久国产精品人妻蜜桃| www.av在线官网国产| 99国产精品一区二区三区| 国产女主播在线喷水免费视频网站| 国产亚洲av高清不卡| 日本午夜av视频| 新久久久久国产一级毛片| 欧美国产精品一级二级三级| 亚洲中文av在线| 五月天丁香电影| 超色免费av| 在线亚洲精品国产二区图片欧美| xxx大片免费视频| 免费在线观看完整版高清| 精品国产国语对白av| 老鸭窝网址在线观看| 欧美亚洲日本最大视频资源| 国产精品亚洲av一区麻豆| 91麻豆精品激情在线观看国产 | 99久久人妻综合| 午夜久久久在线观看| 欧美日本中文国产一区发布| 热99国产精品久久久久久7| 欧美日韩成人在线一区二区| 亚洲少妇的诱惑av| 中文字幕最新亚洲高清| 天天添夜夜摸| 亚洲欧美一区二区三区国产| 女警被强在线播放| 欧美黑人精品巨大| 久久99热这里只频精品6学生| 精品国产一区二区三区久久久樱花| 亚洲专区国产一区二区| 久久这里只有精品19| 国产av一区二区精品久久| 中文字幕av电影在线播放| 中文字幕人妻丝袜制服| 一区二区三区激情视频| 欧美日韩一级在线毛片| 国产淫语在线视频| 一区二区三区四区激情视频| 亚洲国产精品一区二区三区在线| 亚洲av综合色区一区| 一本色道久久久久久精品综合| 亚洲精品一二三| 9191精品国产免费久久| 亚洲精品久久久久久婷婷小说| av天堂在线播放| 国产精品一区二区在线不卡| 欧美激情高清一区二区三区| 久久精品久久久久久久性| 狂野欧美激情性xxxx| 亚洲久久久国产精品| 99国产精品一区二区三区| 亚洲美女黄色视频免费看| 午夜福利视频在线观看免费| 青青草视频在线视频观看| 免费观看av网站的网址| 欧美亚洲日本最大视频资源| 色综合欧美亚洲国产小说| 捣出白浆h1v1| 9热在线视频观看99| 国产极品粉嫩免费观看在线| 两个人免费观看高清视频| 午夜福利视频精品| 久久精品成人免费网站| 极品少妇高潮喷水抽搐| 9色porny在线观看| 狠狠婷婷综合久久久久久88av| 性色av一级| 啦啦啦在线观看免费高清www| 捣出白浆h1v1| 午夜激情av网站| 免费不卡黄色视频| 亚洲av美国av| 国产片特级美女逼逼视频| 视频区图区小说| 宅男免费午夜| 老司机亚洲免费影院| 别揉我奶头~嗯~啊~动态视频 | 这个男人来自地球电影免费观看| 男女午夜视频在线观看| 久久久久久久精品精品| 精品福利永久在线观看| 亚洲国产精品国产精品| 黄色a级毛片大全视频| 欧美在线一区亚洲| 亚洲精品一区蜜桃| 国产成人av激情在线播放| 国产女主播在线喷水免费视频网站| 嫩草影视91久久| 久久热在线av| 满18在线观看网站| 丝瓜视频免费看黄片| 亚洲国产精品一区三区| 久久精品国产亚洲av高清一级| 国产爽快片一区二区三区| av天堂久久9| 在线看a的网站| 亚洲图色成人| 国产成人av激情在线播放| 欧美精品人与动牲交sv欧美| 黑人巨大精品欧美一区二区蜜桃| 两个人看的免费小视频| 欧美在线一区亚洲| 亚洲伊人色综图| 久久久久久久精品精品| 老鸭窝网址在线观看| 精品人妻一区二区三区麻豆| 黄色a级毛片大全视频| 国产成人精品无人区| 久久精品国产a三级三级三级|