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

    穿越河流管道環(huán)焊縫裂紋缺陷安全性研究*

    2023-07-12 08:27:42段宇航董紹華魏昊天楊娜
    石油機械 2023年6期
    關鍵詞:母材屈服河流

    段宇航 董紹華 魏昊天 楊娜

    (1.國家管網(wǎng)集團工程技術創(chuàng)新有限公司 2.中國石油大學(北京)管道技術與安全研究中心)

    0 引 言

    目前,世界各國廣泛采用管道輸送石油天然氣。2021年,我國僅天然氣的表觀消費量就高達3 726×109m3,同比增長12.7%,增速遠超石油和煤炭[1]。為了滿足市場需要,2018至2030年天然氣管道建設將迎來新的大發(fā)展時期,預計此期間將新建油氣管道10萬km。2025年,我國油氣管網(wǎng)總長度將達到24萬km,覆蓋率進一步擴大,結構優(yōu)化升級,儲運能力明顯提升[2-3]。管道工程具有施工距離長、跨越度大和穿越地理情況復雜等特點。每年的管道施工要穿越近百條河流,管道穿越河流時,由于河床表層的土壤不均勻或河水長期沖擊作用而使管線附近的泥沙沖蝕、河床凹凸不平,造成露出管道設施致使其懸空,且受力情況較為復雜。長輸管道通過焊接連接,環(huán)焊縫區(qū)域作為管道的主要失效點,在建設和服役過程中難免會存在未焊透、夾渣、裂紋等焊接問題,從而造成安全隱患[4-5]。因此,穿越河流管道環(huán)焊縫處極可能發(fā)生事故,威脅管道的安全運營。

    在穿越河流管道方面,由小川等[6]采用ABAQUS有限元軟件,研究了地質條件發(fā)生變化時對天然氣輸氣管線強度和剛度的影響,首次在國內利用有限元法進行應力分析,對管道的懸空跨度和埋深進行了分析計算;符煒東等[7]采用ANSYS有限元軟件中的Workbench平臺,對涪陵-王場輸氣管道某穿越河流段懸空管道進行流固耦合計算,得到懸空管段應力特征云圖和管道破壞模式。在含裂紋缺陷管道環(huán)焊縫安全性方面,ExxonMobil公司采用靜裂紋計算模型,分析了管道裂紋位置、深度及長度、管道壁厚、錯邊量等因素對管道環(huán)焊縫的影響,并利用切線法提出了極限拉伸應變計算模型;ZHAO H.S.等[8]將等效應力應變關系方法(ESSRM)和傳統(tǒng)裂紋驅動力J積分的估計法結合,建立了塑性極限載荷計算模型;吳鍇等[9-10]考慮地質災害常引起的永久地表變形,基于參數(shù)化建模技術,研究了裂紋長度、裂紋深度和強度匹配系數(shù)對裂紋尖端張開位移演化的影響。

    國內外眾多學者一直致力于對穿跨越段管道或管道環(huán)焊縫事故的研究,但是很少將兩者結合。為此,筆者以穿越某河流的X70管道環(huán)焊縫為主要研究對象,基于管道局部懸跨段及管土相互作用的簡化力學模型,考慮裂紋尖端處應力的奇異性,建立穿越河流的含裂紋缺陷管道環(huán)焊縫ABAQUS非線性有限元分析模型,研究屈服強度、焊縫匹配系數(shù)、內壓及土壤承載力參數(shù)對管道斷裂韌性值J積分的影響規(guī)律,得到影響管道環(huán)焊縫裂紋起裂的關鍵因素。以期為穿越河流管道安全運行提供理論基礎和參考。

    1 某穿越河流管道實際描述

    1.1 管道工況分析

    對所研究的天然氣長輸管道定期進行外檢測和三軸高清漏磁內檢測。檢測結果表明:穿越河流的管道采用直埋方式敷設,在水流長時間沖刷作用下,土體松軟液化,導致局部完全裸露而懸空,且其懸空段存在環(huán)焊縫裂紋缺陷??紤]無損檢測的精度以及焊接接頭各區(qū)域幾何結構特點等因素,定義臨界裂紋缺陷尺寸為50.0 mm×3.0 mm,即取裂紋長度為50.0 mm,裂紋深度為3.0 mm。該穿越管道的管材為X70,設計壓力為10 MPa,規(guī)格參數(shù)為?1 016 mm×25 mm;外防腐涂層采用環(huán)氧粉末,輸送介質為天然氣;穿越河流河水平均密度取1 020 kg/m3,流速為0~1.5 m/s。

    1.2 管道受力分析

    穿越河流管道懸跨段受力情況復雜,受到的靜力載荷包括:管道自重和內部介質重力以及水流沖刷管道產生的浮力,計算公式為:

    (1)

    (2)

    式中:G為單位長度管道及內部介質重量,N/m;Ff為單位長度管道同水流產生的浮力,N/m;D為管道的外徑,m;d為管道內徑,m;ρo為介質的密度,kg/m3;ρp為管道材質的密度,kg/m3;g為重力加速度,m/s2;ρw為管道所處河水的密度,kg/m3。

    根據(jù)Morison理論可知,管道表面流速是動態(tài)變化的[11],故懸跨管道主要受到平行于河流方向上的拖曳力及慣性力、垂直方向上的漩渦升力作用。單位長度管道受到的水動力的計算式為:

    (3)

    (4)

    (5)

    式中:FL為作用在管道上的旋渦升力,N/m;CL為升力系數(shù),0.31;u為河水流速,m/s;fs為渦激振蕩的泄放頻率,s-1,fs=Stu/D,其中St=0.2;FD為作用在管道上的拖曳力,N/m;CD為升力系數(shù);FⅠ為作用在管道上的慣性力,N/m;CⅠ為慣性力系數(shù)。

    由于河流在垂直方向上的速度、加速度近似為0,故管道垂直于河流方向上受到的慣性力及拖曳力忽略不計[12]。圖1為懸跨管道幾何結構及受載示意圖。圖2中,B為管道壁厚

    圖1 懸跨管道幾何結構及受載剖面示意圖Fig.1 Geometric structure and load section of the suspended section

    1.3 管道材料性能表征

    取實際載荷工況下的X70管道母材及環(huán)焊縫進行常溫下的單軸拉伸試驗(參照GB/T 228.1—2010)和準靜態(tài)斷裂韌性試驗。在X70無縫鋼管上截取管材,加工為標準單軸拉伸試樣,試樣的長度方向為管道軸向。樣品母材及焊縫區(qū)各2根,共計4根,母材以B1、B2標識,焊縫區(qū)材料以W1、W2標識。經(jīng)試驗,得到母材區(qū)及焊縫區(qū)材料強度性能、臨界斷裂韌性值等,如表1所示。

    表1 管道材料參數(shù)Table 1 Pipe material parameters

    根據(jù)試驗得到的工程應力-應變曲線,依據(jù)下式 繪出其真實應力-應變曲線,如圖2和圖3所示[13]。

    σtrue=σnom(1+εnom)

    (6)

    εtrue=ln(1+εnom)

    (7)

    式中:σtrue為真實應力,MPa;σnom為試驗得到的工程應力,MPa;εtrue為真實應變,%;εnom為試驗所得的工程應變,%。

    圖2 試驗所得管材工程應力-應變曲線Fig.2 Measured engineering stress-strain curve of pipe

    圖3 管材真實應力-應變曲線Fig.3 True stress-strain curve of pipe

    2 含表面裂紋的管道有限元模型

    本文主要根據(jù)ABAQUS軟件模擬管道實際工況,對X70管材進行有限元分析。選取校正跨長為16 m的管道,其跨肩長度選取懸跨段的,即為8 m;環(huán)焊縫位于管道最中央,寬度為10 mm,兩端設置土彈簧單元??紤]裂紋和管道承受載荷以及幾何結構上的對稱性,同時為了提高計算效率,通過非線性數(shù)值仿真軟件ABAQUS建立含局部外表面裂紋管道環(huán)焊縫的對稱模型。

    2.1 管土相互作用模型

    依據(jù)ASCE準則、DNV-RP-F114及SY/T 7056—2016等相關的準則確定土單元本構關系。通常情況下,土彈簧的模型本構關系將埋設管道的土壤視為均勻的,通過計算土壤的相對位移與極限承載力間的關系,來確定土彈簧的剛度。土彈簧模型如圖4所示,KU、Kh、Kv分別表示3個方向上的土剛度系數(shù)[14]。

    圖4 土彈簧模型Fig.4 Spring model of soil

    選取的土質為松黏土,土壤的抗剪強度梯度值取0.5 kPa/m,河床土壤的黏聚力取10 kPa,內摩擦角取20°,容量為21 000 N/m3。穿越河流管道懸跨段埋深1.3 m。松黏土彈簧剛度值如表2所示,管土的本構關系如圖5所示。

    表2 松黏土彈簧剛度計算值Table 2 Calculated spring stiffness for loose clay soils

    圖5 土單元的本構關系Fig.5 Constitutive relation of soil elements

    2.2 有限元模型

    由圖5可知,土體并不具有嚴格的線性性質。此處充分考慮土壤的非線性作用,并未使用ABAQUS自帶傳統(tǒng)線性彈簧模型,而是創(chuàng)新性地用連接器建立管土非線性耦合單元,如圖6所示。采取實體切割劃分方法在懸跨中建立管道環(huán)焊縫模型。管體網(wǎng)格劃分采用八節(jié)點完全積分單元(C3D8),對環(huán)焊縫附近區(qū)域網(wǎng)格加密,以保證缺陷處計算的準確性?;赑ython參數(shù)化編程方法,采用“Canoe”型靜裂紋法模擬管道環(huán)焊縫位置表面型裂紋缺陷,并可基于圍線積分自動輸出J積分斷裂參數(shù)值。為更好地模擬裂紋在張開過程中的鈍化行為,裂紋尖端采用半徑為0.02 mm的楔形退化單元表征裂紋尖端應力奇異。模型采用蜘蛛網(wǎng)型的局部網(wǎng)格加密方式對裂紋面附近的網(wǎng)格進行細化,以消除變形過程中的應力奇異效應,改善計算收斂性,以及提高模型的計算精度。裂紋處有限元模型如圖7所示。

    圖6 管土非線性作用模型Fig.6 Nonlinear pipe-soil interaction model

    圖7 管道環(huán)焊縫裂紋模型Fig.7 Pipeline girth weld crack model

    模型共設置4個邊界條件,并對管道模型對稱面YZ面和環(huán)焊縫XY面采用對稱約束,土體3個方向上的位移轉角均設置為0,埋入土壤中管道最右端完全固定。設2個加載步:第一步,對管道內壁施加內壓載荷;第二步,保持內壓載荷不變對管道懸跨段管體施加表面應力,模擬管道受到的外力合力。

    2.3 模型驗證

    為了驗證穿越河流管道有限元數(shù)值模型的準確性及可靠性,需將施加載荷后有限元計算結果與解析法計算出的結果進行對比。在不設置管道環(huán)焊縫時,驗證穿越河流管道模型的準確性,將模型計算出的管道最大變形量與基于Winkler彈性地基梁理論[15]計算結果進行比較,可知2種方法的誤差為8.39%,在可接受范圍內。

    基于BS 7910—2019標準驗證裂紋模型的準確性,將有限元計算得到的J積分通過式(8)轉換為應用拉應力強度因子(KI),將其與BS 7910標準中關于管道外表面周向裂紋的應力強度因子解析法(式(9))計算出的結果進行對比,來驗證所建裂紋有限元模型的準確性[16-17]。

    (8)

    (9)

    Yσ=Mfw{ktmMkmMmpm+ktbMkbMb×

    [pb+(km-1)pm]}

    (10)

    式中:J為從一個裂紋面到另一個裂紋面包圍裂紋前緣的線或面積分,用來表征裂紋尖端周圍的局部應力-應變場,N/mm;K1為應用拉應力強度因子,N/mm;E為彈性模量,MPa;σ為應力,N/mm2;Y為應力強度修正系數(shù),無量綱;a為原缺陷或再特種缺陷,mm;fw為有限寬度修正因子;ktm為結構不邊疆性產生的對膜應力集中系數(shù);ktb和km分別為結構不連續(xù)性及結構錯位產生的應力集中系數(shù);M為膨脹修正因子;Mkm、Mm、Mkb、Mb均為應力強度放大系數(shù),無量綱(根據(jù)BS7910—2019規(guī)范查找使用);pb為一次彎曲應力,MPa;pm為薄膜應力,MPa。

    圖8 裂紋有限元模型的驗證結果Fig.8 Validation results of the crack finite element model

    由于BS7910—2019規(guī)范中提供的方法適用于計算線彈性情況,所以在有限元仿真時輸入純彈性材料參數(shù)(彈性模量E為210 GPa;泊松比υ為0.3)。圖8給出了?420 mm×18 mm管道3種環(huán)向裂紋深度在不同情況下應力強度因子值的對比,圖中a為裂紋深度,t為壁厚。有限元計算結果與BS 7910—2019計算所得結果較為一致,誤差均未超過10%。因此,建立的管道表面裂紋有限元仿真模型較為準確。

    2.4 結果分析

    根據(jù)圖9可視化云圖,可以看出其應力在裂紋缺陷處最大,遠離管道裂紋區(qū)域的數(shù)值緩慢減小,紅色部分位于管道環(huán)焊縫裂紋區(qū)域。在復雜組合應力的作用下,管道環(huán)焊縫缺陷處的等效應力最大值遠超過X70管線鋼材料屈服強度,而其裂紋J積分值達到材料臨界斷裂韌性值。因此,該穿越河流管道發(fā)生裂紋擴展,最終可能發(fā)展為環(huán)焊縫處斷裂失效。

    圖9 穿越河流管道環(huán)焊縫等效應力云圖Fig.9 Equivalent stress nephogram of girth weld of river-crossing pipeline

    3 基于J 積分的敏感性分析

    采用前文所述有限元模型,系統(tǒng)研究管道屈服強度、焊縫匹配系數(shù)、內壓及土壤承載力變化,及其對穿越河流管道環(huán)焊縫裂紋處的影響規(guī)律,以識別影響管道環(huán)焊縫裂紋起裂的主要因素。

    3.1 材料特性影響因素

    裂紋位于環(huán)焊縫中心線位置,因此環(huán)焊接頭材料主要考慮母材以及焊縫金屬。基于材料應力-應變關系的本構模型,認為焊縫以及母材的應力-應變均滿足冪次強化規(guī)律;在彈性階段,其應力-應變關系滿足下式:

    σ=Eεσ<σy

    (11)

    式中:σy為材料的屈服強度,MPa;ε為材料應變,%。對于焊縫與母材彈性模量E取210 GPa。

    在塑性部分,管道母材與環(huán)焊縫的應力與塑性應變關系滿足下式:

    (12)

    式中:εp為塑性應變;εy=σy/E,為屈服應力σy對應的應變;n為材料的應變硬化系數(shù),焊縫與母材的硬化系數(shù)相同取n=0.06[18]。

    定義環(huán)焊縫的屈服強度與母材的屈服強度之比稱作強度匹配系數(shù)m:

    (13)

    式中:σSW為環(huán)焊縫的屈服強度,MPa;σSB為母材的屈服強度,MPa。

    對于母材屈服強度的影響考慮分別為600、550和500 MPa 3種情況;以及母材屈服強度σy為600 MPa時強度匹配系數(shù)為0.8、0.9、1.0和1.1的4種情況。圖10分別展示了不同屈服強度母材及強度匹配系數(shù)條件下焊縫的真實應力—應變關系曲線。

    基于應力的方法以力或應力作為評價和設計的參數(shù)條件,此時材料強度對結果的影響較大。橫坐標河流流速是已知管道所受合力,其大小通過1.2章節(jié)中水流橫向作用載荷公式(3)推算而來。圖11和圖12分別為不同屈服強度、不同焊縫強度匹配系數(shù)對應的J積分隨河流流速變化的曲線。

    圖11 管道屈服強度對J積分的影響 Fig.11 Influence of pipe yield strength on J-integral

    圖12 焊縫強度匹配系數(shù)對J積分的影響Fig.12 Influence of weld strength matching coefficient on J-integral

    從圖11中可以看出,隨著河流流速的逐漸增加,穿越河流管道環(huán)焊縫處J積分逐漸增大,且同一流速下,管材屈服強度越低,對應的J積分值越大。分析認為,懸跨管道所受合力為重力方向,河流流速越大,管道承受的載荷越大,J積分也隨之增大;管材屈服強度高,則管道整體承載能力強,則J積分隨流速增長較緩。由圖12可知,焊縫強度匹配系數(shù)越高裂紋越不易進入屈服,因此相同載荷條件下強度匹配系數(shù)越高的工況,對應的J積分值就越小,且匹配系數(shù)越低,J積分與河流流速的關系曲線也更加陡峭。

    3.2 內壓影響

    穿越河流管道受內壓引起的環(huán)向應力與外界載荷作用力的共同作用而處于雙軸應力狀態(tài)。通常來說,在僅有內壓作用的情況下,增加內壓會導致J積分增加,但影響較小。僅施加內壓對管道J積分結果影響如圖13所示。圖14為模擬實際工況下含表面初始裂紋X70管道環(huán)焊縫,在內壓變化為7、8、9、10和11 MPa下對應的J積分隨河流流速的變化曲線。

    從圖14可以看出,在內壓及外界載荷共同作用下,裂紋處J積分隨著河流流速的增大而增大,且低內壓條件下J積分的增加幅度明顯低于高內壓條件下的結果。由1.3結論可知,試驗所得環(huán)焊縫臨界斷裂韌性值JQ0.2BL為192.094 kJ/m2,該值可作為焊縫區(qū)的臨界斷裂韌性值JIc,可基于裂紋起裂的失效準則來確定不同內壓條件下的安全河流流速,結果如圖14所示。

    圖13 僅內壓作用下對J積分的影響Fig.13 Influence of pure internal pressure on J-integral

    圖14 雙軸應力作用下內壓對J積分的影響Fig.14 Influence of internal pressure on J-integral under biaxial stress

    3.3 土壤承載力影響

    穿越河流管道入土端土壤性質對管道的變形會產生一定程度的影響。管土間的相互作用可用土彈簧進行描述,埋深、管道類型及土壤特征均會影響土彈簧剛度取值,導致有限元計算結果的差異。如圖15所示為等比例改變土壤的極限承載力,即土壤剛度(土彈簧剛度),此條件下分析J積分隨河流流速的變化規(guī)律。

    圖15 土壤承載力對J積分的影響Fig.15 Influence of soil bearing capacity on J-integral

    從圖15中可以看出,管道周圍土壤承載力對環(huán)焊縫J積分有著一定的影響作用。即管道在相同載荷作用下,J積分隨流速增大而逐漸增大,且在相同流速下,土彈簧剛度越小的情況對應的J積分值越大。這是由于管道周圍土壤的彈性支撐作用,土彈簧剛度的增加會削弱管道整體的彎曲變形程度,減緩了復雜載荷作用下管道的變形程度。所以土彈簧剛度越大,J積分越小。

    4 結 論

    (1)穿越河流管道受到復雜應力的作用,運行條件惡劣。管道環(huán)焊縫處是管道薄弱環(huán)節(jié),最易發(fā)生斷裂失效,而其結構完整性是管道完整性評估的關鍵。以穿越某河流的X70天然氣長輸管道為主要研究對象進行應力分析,包括靜載荷分析和由Morison理論計算的荷漩渦載荷及波浪流載。

    (2)基于ABAQUS軟件,建立了含裂紋穿越河流懸跨管道環(huán)焊縫非線性有限元分析模型。模型創(chuàng)新性地使用連接器建立管土非線性耦合單元,避免使用ABAQUS自帶傳統(tǒng)彈簧模型僅考慮土壤線性性質的弊端?;赑ython參數(shù)化編程方法,建立了可導出管道裂紋驅動力的有限元模型,能夠準確描述裂紋尖端鈍化、起裂等斷裂行為。開發(fā)了以J積分為表征的裂紋驅動力智能后處理提取方法,且經(jīng)驗證后發(fā)現(xiàn)模型輸出結果準確。

    (3)采用單因素敏感性分析法研究不同影響因素對J積分的影響規(guī)律。針對含相同裂紋型缺陷的管道,可以通過提高環(huán)焊縫材料屈服強度、強度匹配系數(shù)、降低管道內壓來提高管道的抵抗性。由于環(huán)焊縫斷裂影響參數(shù)較多,后續(xù)需進一步研究管道類型、裂紋幾何尺寸等因素的影響,最終考慮各因素之間相互耦合作用下對含裂紋管道環(huán)焊縫承載能力的影響規(guī)律。

    猜你喜歡
    母材屈服河流
    牙被拔光也不屈服的史良大律師秘書
    紅巖春秋(2022年1期)2022-04-12 00:37:34
    銅與不銹鋼激光焊接頭金相組織
    河流
    小太陽畫報(2019年4期)2019-06-11 10:29:48
    The Classic Lines of A Love so Beautiful
    流放自己的河流
    散文詩(2018年20期)2018-05-06 08:03:44
    Sn/Al真空釬焊壓力下母材氧化膜的破碎效果
    母材及焊絲對轉向架用耐候鋼焊接接頭性能的影響
    電焊機(2016年8期)2016-12-06 08:28:44
    勇敢
    當河流遇見海
    雙聯(lián)法冶煉鐵基非晶合金母材的工業(yè)化技術開發(fā)
    山東冶金(2015年5期)2015-12-10 03:27:40
    色播亚洲综合网| 首页视频小说图片口味搜索| 两个人看的免费小视频| h日本视频在线播放| 国产男靠女视频免费网站| 国产主播在线观看一区二区| 亚洲欧美一区二区三区黑人| 美女扒开内裤让男人捅视频| 成年版毛片免费区| 国产蜜桃级精品一区二区三区| 久久中文字幕一级| 九九热线精品视视频播放| 俄罗斯特黄特色一大片| 久久草成人影院| 亚洲国产欧美网| 欧美日韩黄片免| 宅男免费午夜| 这个男人来自地球电影免费观看| 黑人欧美特级aaaaaa片| 黑人巨大精品欧美一区二区mp4| 国产精品香港三级国产av潘金莲| a级毛片在线看网站| 99热这里只有精品一区 | 人人妻人人澡欧美一区二区| 国产三级黄色录像| 男女视频在线观看网站免费| 69av精品久久久久久| 国产乱人视频| 91麻豆精品激情在线观看国产| 听说在线观看完整版免费高清| 亚洲七黄色美女视频| 99精品久久久久人妻精品| 欧美国产日韩亚洲一区| 久久精品国产综合久久久| 亚洲熟女毛片儿| 搡老岳熟女国产| 国产伦精品一区二区三区视频9 | 丁香六月欧美| 国产精品野战在线观看| 又黄又粗又硬又大视频| 午夜两性在线视频| 国产午夜福利久久久久久| 好男人电影高清在线观看| 美女午夜性视频免费| 美女黄网站色视频| 视频区欧美日本亚洲| а√天堂www在线а√下载| 18禁裸乳无遮挡免费网站照片| 黄色片一级片一级黄色片| 国产亚洲欧美98| 中文在线观看免费www的网站| 精品国产乱子伦一区二区三区| 无限看片的www在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 国产黄片美女视频| 婷婷精品国产亚洲av在线| 视频区欧美日本亚洲| 无人区码免费观看不卡| 特级一级黄色大片| 久久精品国产亚洲av香蕉五月| 制服丝袜大香蕉在线| 免费看十八禁软件| 日韩av在线大香蕉| 国产精品98久久久久久宅男小说| 亚洲av五月六月丁香网| av国产免费在线观看| 精品久久久久久,| 亚洲狠狠婷婷综合久久图片| 女人被狂操c到高潮| 老司机深夜福利视频在线观看| 国产三级在线视频| 99久久无色码亚洲精品果冻| 最近最新中文字幕大全电影3| 成人特级黄色片久久久久久久| 丰满人妻一区二区三区视频av | 国产成人精品无人区| 性欧美人与动物交配| 亚洲自偷自拍图片 自拍| 成人三级黄色视频| 亚洲黑人精品在线| 亚洲在线观看片| 每晚都被弄得嗷嗷叫到高潮| 淫秽高清视频在线观看| 18禁黄网站禁片午夜丰满| 黄色丝袜av网址大全| 成年女人毛片免费观看观看9| 日本 欧美在线| 午夜激情福利司机影院| 久久国产精品人妻蜜桃| 亚洲精品一卡2卡三卡4卡5卡| 午夜福利在线在线| 老汉色av国产亚洲站长工具| 天天躁日日操中文字幕| 无人区码免费观看不卡| 免费在线观看日本一区| 亚洲精品在线美女| 亚洲激情在线av| 精品免费久久久久久久清纯| 欧美中文日本在线观看视频| 99国产精品99久久久久| 亚洲va日本ⅴa欧美va伊人久久| av天堂在线播放| 色综合站精品国产| а√天堂www在线а√下载| 99精品久久久久人妻精品| 桃色一区二区三区在线观看| 亚洲av电影不卡..在线观看| 成人精品一区二区免费| 麻豆久久精品国产亚洲av| 男女做爰动态图高潮gif福利片| 亚洲性夜色夜夜综合| 精品久久久久久久毛片微露脸| a级毛片在线看网站| 特大巨黑吊av在线直播| 在线观看舔阴道视频| 日韩免费av在线播放| 97超视频在线观看视频| 亚洲五月天丁香| 又爽又黄无遮挡网站| 一级毛片女人18水好多| 亚洲欧美精品综合一区二区三区| 久久精品影院6| 国内精品久久久久久久电影| 变态另类成人亚洲欧美熟女| 欧美黄色片欧美黄色片| 美女高潮的动态| 国产毛片a区久久久久| 日韩欧美一区二区三区在线观看| 亚洲人与动物交配视频| 手机成人av网站| 18禁国产床啪视频网站| 99久久精品国产亚洲精品| 国产伦在线观看视频一区| 国产一区二区在线观看日韩 | 欧美性猛交╳xxx乱大交人| 久久草成人影院| av天堂中文字幕网| 国产三级黄色录像| 成人三级做爰电影| 一进一出抽搐gif免费好疼| 听说在线观看完整版免费高清| 成人亚洲精品av一区二区| 亚洲人成网站在线播放欧美日韩| 一个人免费在线观看的高清视频| 亚洲第一电影网av| 亚洲精品在线观看二区| 俺也久久电影网| 一级毛片精品| 日本 av在线| 久久性视频一级片| 熟妇人妻久久中文字幕3abv| 久久这里只有精品中国| 亚洲国产色片| 老熟妇乱子伦视频在线观看| 欧美日韩中文字幕国产精品一区二区三区| 国产久久久一区二区三区| 精品国产超薄肉色丝袜足j| 精品久久久久久久毛片微露脸| 色av中文字幕| 伊人久久大香线蕉亚洲五| 一个人免费在线观看电影 | 久久久水蜜桃国产精品网| 免费在线观看亚洲国产| 全区人妻精品视频| 亚洲精品国产精品久久久不卡| 色视频www国产| 网址你懂的国产日韩在线| www.999成人在线观看| 97超级碰碰碰精品色视频在线观看| 欧美不卡视频在线免费观看| 精品久久久久久久毛片微露脸| 欧美极品一区二区三区四区| 亚洲专区中文字幕在线| 国模一区二区三区四区视频 | 99国产精品一区二区蜜桃av| 色综合婷婷激情| av欧美777| 别揉我奶头~嗯~啊~动态视频| 一个人看视频在线观看www免费 | 操出白浆在线播放| 99视频精品全部免费 在线 | 亚洲 欧美 日韩 在线 免费| 两个人视频免费观看高清| 国产精品一区二区免费欧美| 18禁黄网站禁片免费观看直播| 久久九九热精品免费| 亚洲无线观看免费| 日韩三级视频一区二区三区| 色尼玛亚洲综合影院| 欧美一区二区精品小视频在线| 欧美成人性av电影在线观看| 国产1区2区3区精品| 国产野战对白在线观看| 亚洲国产欧美人成| 在线观看美女被高潮喷水网站 | 无遮挡黄片免费观看| 色哟哟哟哟哟哟| 九九久久精品国产亚洲av麻豆 | 精品久久久久久成人av| 少妇熟女aⅴ在线视频| 久久久国产欧美日韩av| 精品久久久久久久毛片微露脸| 在线免费观看不下载黄p国产 | 一区二区三区高清视频在线| 老鸭窝网址在线观看| 波多野结衣巨乳人妻| 亚洲片人在线观看| 欧美丝袜亚洲另类 | 久久久国产成人免费| 亚洲欧美日韩东京热| 欧美黑人欧美精品刺激| 亚洲avbb在线观看| 久久久久久人人人人人| 桃红色精品国产亚洲av| 国产免费男女视频| 国产精品日韩av在线免费观看| 熟女少妇亚洲综合色aaa.| 99久久久亚洲精品蜜臀av| 亚洲av第一区精品v没综合| 91老司机精品| 国产成人系列免费观看| 久久这里只有精品19| 国产精品亚洲一级av第二区| 俄罗斯特黄特色一大片| 精品久久久久久久久久免费视频| av欧美777| 国产成人一区二区三区免费视频网站| 成人特级黄色片久久久久久久| 免费大片18禁| 看黄色毛片网站| 国产免费av片在线观看野外av| 美女 人体艺术 gogo| 黑人巨大精品欧美一区二区mp4| 叶爱在线成人免费视频播放| 窝窝影院91人妻| 国产不卡一卡二| 国产成人av教育| 国产亚洲av高清不卡| 国产精品女同一区二区软件 | 偷拍熟女少妇极品色| 精品久久久久久久久久免费视频| 国产免费av片在线观看野外av| 亚洲人与动物交配视频| 日本黄大片高清| 亚洲熟妇熟女久久| 成年免费大片在线观看| 黄色 视频免费看| 长腿黑丝高跟| 亚洲精品456在线播放app | 99久国产av精品| 亚洲成人久久性| 69av精品久久久久久| 人人妻人人看人人澡| 久久中文看片网| 中出人妻视频一区二区| 婷婷精品国产亚洲av在线| 毛片女人毛片| 成年版毛片免费区| 亚洲av熟女| 少妇的逼水好多| 91麻豆精品激情在线观看国产| 国产97色在线日韩免费| 好男人电影高清在线观看| www.www免费av| 免费在线观看视频国产中文字幕亚洲| or卡值多少钱| 国产伦精品一区二区三区四那| 国产精品av视频在线免费观看| av片东京热男人的天堂| 老汉色∧v一级毛片| 亚洲精华国产精华精| 91九色精品人成在线观看| 国产精品一区二区三区四区久久| 国产激情久久老熟女| 精品人妻1区二区| 欧美色视频一区免费| av天堂中文字幕网| 欧美黑人巨大hd| 熟妇人妻久久中文字幕3abv| 久久香蕉国产精品| 看免费av毛片| 日本精品一区二区三区蜜桃| 最新美女视频免费是黄的| 亚洲国产看品久久| 久久草成人影院| 国产精品久久久久久亚洲av鲁大| 国产美女午夜福利| 午夜视频精品福利| 久久久国产成人免费| 精品不卡国产一区二区三区| a级毛片在线看网站| 欧美绝顶高潮抽搐喷水| 免费看美女性在线毛片视频| 99久久99久久久精品蜜桃| 免费高清视频大片| 淫妇啪啪啪对白视频| 久久精品国产亚洲av香蕉五月| 欧美中文日本在线观看视频| 亚洲成人精品中文字幕电影| 免费在线观看亚洲国产| 99国产极品粉嫩在线观看| 欧美色视频一区免费| 九九在线视频观看精品| 精品久久久久久久末码| 国产精品1区2区在线观看.| 亚洲精品中文字幕一二三四区| netflix在线观看网站| 波多野结衣高清作品| 中文字幕最新亚洲高清| 在线永久观看黄色视频| 精品久久久久久久久久久久久| 五月伊人婷婷丁香| 免费看日本二区| 午夜福利在线在线| 欧美色欧美亚洲另类二区| 日韩免费av在线播放| 男人和女人高潮做爰伦理| 三级男女做爰猛烈吃奶摸视频| 精品免费久久久久久久清纯| 亚洲最大成人中文| 99国产精品一区二区蜜桃av| 欧美不卡视频在线免费观看| 亚洲黑人精品在线| 国产精品98久久久久久宅男小说| 精品久久久久久久末码| 麻豆成人av在线观看| 日韩av在线大香蕉| 免费无遮挡裸体视频| 一个人免费在线观看电影 | 成人无遮挡网站| 岛国在线观看网站| 亚洲aⅴ乱码一区二区在线播放| 亚洲一区二区三区不卡视频| 免费看美女性在线毛片视频| 一个人看视频在线观看www免费 | 蜜桃久久精品国产亚洲av| 搡老岳熟女国产| 女人被狂操c到高潮| 国产男靠女视频免费网站| 热99re8久久精品国产| 两个人的视频大全免费| 国产美女午夜福利| 亚洲精品在线观看二区| 夜夜夜夜夜久久久久| 国产精品久久久久久亚洲av鲁大| 国产视频内射| 久久久国产欧美日韩av| 亚洲第一欧美日韩一区二区三区| 午夜福利欧美成人| 女警被强在线播放| 亚洲专区国产一区二区| 波多野结衣高清作品| 日本一二三区视频观看| 看免费av毛片| 国产精品久久久av美女十八| 热99re8久久精品国产| 男人舔女人的私密视频| 免费在线观看亚洲国产| 桃色一区二区三区在线观看| 午夜免费成人在线视频| 日本黄色片子视频| 巨乳人妻的诱惑在线观看| 欧美丝袜亚洲另类 | svipshipincom国产片| 久久久国产欧美日韩av| 人妻久久中文字幕网| 欧美在线黄色| 欧美国产日韩亚洲一区| 在线免费观看的www视频| 中文在线观看免费www的网站| 九九久久精品国产亚洲av麻豆 | 97人妻精品一区二区三区麻豆| 欧美3d第一页| 中文字幕人成人乱码亚洲影| www.精华液| 99久久精品热视频| 黑人欧美特级aaaaaa片| 久久久水蜜桃国产精品网| 可以在线观看毛片的网站| 国产三级黄色录像| 欧美乱码精品一区二区三区| 搡老妇女老女人老熟妇| 国语自产精品视频在线第100页| 国产精品久久久久久人妻精品电影| 日韩精品中文字幕看吧| 麻豆av在线久日| 久久香蕉精品热| 免费观看的影片在线观看| 日本黄色视频三级网站网址| 一级作爱视频免费观看| www.自偷自拍.com| 国产69精品久久久久777片 | 看免费av毛片| 噜噜噜噜噜久久久久久91| 欧美乱码精品一区二区三区| 岛国在线免费视频观看| 变态另类成人亚洲欧美熟女| 午夜精品在线福利| 国产乱人视频| 久久久国产成人精品二区| www日本黄色视频网| 人妻久久中文字幕网| 亚洲美女视频黄频| 精品久久久久久久久久免费视频| 亚洲国产欧美人成| 成人永久免费在线观看视频| 国产精品久久久久久人妻精品电影| 免费av不卡在线播放| 午夜免费激情av| av福利片在线观看| 精品欧美国产一区二区三| 免费在线观看影片大全网站| 久久久久性生活片| 色噜噜av男人的天堂激情| 蜜桃久久精品国产亚洲av| 一进一出好大好爽视频| 怎么达到女性高潮| 舔av片在线| 99国产精品一区二区三区| 亚洲国产精品成人综合色| 一二三四在线观看免费中文在| 桃红色精品国产亚洲av| 亚洲专区中文字幕在线| 少妇的丰满在线观看| 亚洲国产看品久久| 国产高清三级在线| 久久久久久国产a免费观看| 欧美日本视频| 亚洲七黄色美女视频| 国产精品99久久99久久久不卡| 性色avwww在线观看| 身体一侧抽搐| 国产成人av教育| 日韩成人在线观看一区二区三区| 法律面前人人平等表现在哪些方面| 69av精品久久久久久| 亚洲国产日韩欧美精品在线观看 | 两人在一起打扑克的视频| 91av网一区二区| 色尼玛亚洲综合影院| 欧美黑人欧美精品刺激| 亚洲av成人一区二区三| 亚洲欧美日韩卡通动漫| 亚洲熟女毛片儿| 国产真实乱freesex| 免费看a级黄色片| 一卡2卡三卡四卡精品乱码亚洲| 黄色女人牲交| 99热精品在线国产| 91av网站免费观看| 欧美日韩综合久久久久久 | 亚洲熟女毛片儿| 日韩有码中文字幕| 18禁黄网站禁片免费观看直播| 亚洲第一电影网av| 麻豆国产97在线/欧美| 香蕉久久夜色| 18美女黄网站色大片免费观看| 色哟哟哟哟哟哟| 欧美色欧美亚洲另类二区| 午夜福利欧美成人| 九色国产91popny在线| 琪琪午夜伦伦电影理论片6080| 又黄又粗又硬又大视频| 性色avwww在线观看| 精品国产乱码久久久久久男人| 国产探花在线观看一区二区| 深夜精品福利| 日本a在线网址| 黑人欧美特级aaaaaa片| 国产精品久久久人人做人人爽| 久久久久精品国产欧美久久久| 久久久成人免费电影| 高清在线国产一区| 黄片小视频在线播放| 夜夜躁狠狠躁天天躁| 久久中文看片网| 日韩欧美国产在线观看| 色噜噜av男人的天堂激情| 熟女少妇亚洲综合色aaa.| avwww免费| 99热这里只有精品一区 | 九九久久精品国产亚洲av麻豆 | 国产精品亚洲一级av第二区| 中国美女看黄片| 国产毛片a区久久久久| 亚洲aⅴ乱码一区二区在线播放| 手机成人av网站| 久久天躁狠狠躁夜夜2o2o| h日本视频在线播放| 嫩草影院入口| 国产精品99久久99久久久不卡| a级毛片在线看网站| h日本视频在线播放| 免费看美女性在线毛片视频| 亚洲av片天天在线观看| 床上黄色一级片| 久久人妻av系列| 老汉色av国产亚洲站长工具| 国内精品久久久久精免费| 亚洲,欧美精品.| xxx96com| 国产精品自产拍在线观看55亚洲| 一二三四在线观看免费中文在| 日韩有码中文字幕| 日韩 欧美 亚洲 中文字幕| 精品一区二区三区av网在线观看| 听说在线观看完整版免费高清| 亚洲欧美精品综合一区二区三区| bbb黄色大片| 99re在线观看精品视频| 日韩大尺度精品在线看网址| 午夜亚洲福利在线播放| 国产99白浆流出| 久久久国产成人精品二区| 精品乱码久久久久久99久播| 国产伦一二天堂av在线观看| 美女cb高潮喷水在线观看 | www.熟女人妻精品国产| 久久久国产精品麻豆| 国语自产精品视频在线第100页| 久久久色成人| 日韩成人在线观看一区二区三区| 亚洲人成网站高清观看| 搡老妇女老女人老熟妇| 1000部很黄的大片| 欧美极品一区二区三区四区| 麻豆国产97在线/欧美| 国产精品国产高清国产av| 真人一进一出gif抽搐免费| 在线视频色国产色| 男女那种视频在线观看| 香蕉丝袜av| 精品国产超薄肉色丝袜足j| 桃色一区二区三区在线观看| 在线免费观看的www视频| 中文字幕高清在线视频| 久久中文字幕人妻熟女| 国产精品99久久久久久久久| 亚洲国产精品合色在线| 国产av不卡久久| 国产又黄又爽又无遮挡在线| 久久这里只有精品中国| 欧美日韩亚洲国产一区二区在线观看| 国产一区二区在线观看日韩 | 成年人黄色毛片网站| 午夜福利在线观看免费完整高清在 | 一个人观看的视频www高清免费观看 | 小说图片视频综合网站| 中国美女看黄片| 亚洲一区二区三区色噜噜| 久久久国产成人免费| 国产不卡一卡二| 国产激情偷乱视频一区二区| 色精品久久人妻99蜜桃| 国产男靠女视频免费网站| 两个人看的免费小视频| 欧美日韩乱码在线| 12—13女人毛片做爰片一| 精品久久久久久成人av| 一个人免费在线观看的高清视频| 日本成人三级电影网站| 欧美黄色淫秽网站| 天堂√8在线中文| 老司机在亚洲福利影院| aaaaa片日本免费| 亚洲av美国av| 在线十欧美十亚洲十日本专区| 真人一进一出gif抽搐免费| 男女下面进入的视频免费午夜| 免费一级毛片在线播放高清视频| 麻豆av在线久日| 日本与韩国留学比较| 午夜精品一区二区三区免费看| 国产成人精品无人区| 成人国产一区最新在线观看| 亚洲欧洲精品一区二区精品久久久| 精品国产乱码久久久久久男人| av在线蜜桃| 欧美另类亚洲清纯唯美| 国产一区二区三区在线臀色熟女| 精品国产亚洲在线| 亚洲av电影不卡..在线观看| a级毛片a级免费在线| 国产亚洲精品久久久久久毛片| 免费看日本二区| 国产精品日韩av在线免费观看| 女同久久另类99精品国产91| 免费看日本二区| 琪琪午夜伦伦电影理论片6080| 国产亚洲av高清不卡| 久久精品国产99精品国产亚洲性色| 天堂√8在线中文| 一二三四在线观看免费中文在| 美女高潮的动态| 19禁男女啪啪无遮挡网站| 三级男女做爰猛烈吃奶摸视频| 国产精品亚洲av一区麻豆| 国产午夜精品久久久久久| 三级男女做爰猛烈吃奶摸视频| 又爽又黄无遮挡网站| 亚洲中文av在线| 夜夜看夜夜爽夜夜摸| 人妻久久中文字幕网| 噜噜噜噜噜久久久久久91| 亚洲一区二区三区色噜噜| 日本在线视频免费播放| 国产精品女同一区二区软件 | 婷婷丁香在线五月| 淫秽高清视频在线观看| 亚洲中文字幕一区二区三区有码在线看 | 日韩人妻高清精品专区| 中文字幕最新亚洲高清| 成人午夜高清在线视频| 日本五十路高清|