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

    非粘接柔性立管抗拉和抗扭剛度計算分析方法對比

    2018-01-09 02:16:58裴曉梅張恩勇李麗瑋周巍偉
    海洋工程裝備與技術(shù) 2017年5期
    關(guān)鍵詞:鎧裝抗拉立管

    裴曉梅, 張恩勇, 李麗瑋, 周巍偉, 曹 靜

    (中海油研究總院,北京 100028)

    非粘接柔性立管抗拉和抗扭剛度計算分析方法對比

    裴曉梅, 張恩勇, 李麗瑋, 周巍偉, 曹 靜

    (中海油研究總院,北京 100028)

    針對非粘接柔性立管抗拉和抗扭剛度計算方法進行對比。以300m水深8英寸(1英寸≈2.54cm)非粘接柔性立管為例,首先介紹了立管截面設(shè)計時計算兩種剛度所使用的理論計算方法,之后在理論計算的基礎(chǔ)上,對非粘接柔性立管進行有限元分析,分別對比了非粘接柔性立管在軸向拉力及扭矩作用下的抗拉和抗扭剛度的計算結(jié)果,并對各層結(jié)構(gòu)的剛度貢獻情況進行了討論。最終得出結(jié)論: 柔性立管在軸向載荷和扭矩作用下,理論計算結(jié)果與有限元計算結(jié)果較為一致,在抗扭剛度計算中理論計算結(jié)果更偏保守。在實際工程設(shè)計中,可以根據(jù)設(shè)計階段選擇合適的計算方法。

    非粘接柔性立管;抗拉剛度;抗扭剛度;接觸

    0 引 言

    為了防止結(jié)構(gòu)失效,立管在設(shè)計階段需要進行嚴格的分析計算和測試。立管結(jié)構(gòu)在海面上主要承受較高的循環(huán)張力載荷,而在海底主要承受較大的外壓、軸向壓縮、彎曲以及扭轉(zhuǎn)組合載荷[1]。同時,非粘接柔性立管(軟管)由不同的層組成,從內(nèi)往外包括骨架層、聚合物材料的內(nèi)部水密層、互鎖結(jié)構(gòu)形式的抗壓鎧裝層、耐磨層、兩層抗拉鎧裝層和外部水密層,各層都有不同的功能作用[2],各層間可以自由地相互滑動。柔性立管的剛度計算是截面設(shè)計中的重要環(huán)節(jié),是關(guān)乎后續(xù)整體分析準確性的重要參數(shù)。根據(jù)理論計算公式和ABAQUS軟件有限元分析,本文針對300m水深8英寸(1英寸≈2.54cm)輸油軟管,分別對柔性立管抗拉和抗扭剛度進行計算和對比,并對各層結(jié)構(gòu)的剛度貢獻情況進行了討論?;诒疚难芯拷Y(jié)果,在實際工程設(shè)計中,可以根據(jù)設(shè)計階段選擇合適的計算方法。

    1 分析對象參數(shù)

    圖1所示為典型的非粘接柔性立管的各層結(jié)構(gòu)。

    圖1 典型的非粘接柔性立管結(jié)構(gòu)Fig.1 Inner structure of a typical unbonded flexible riser

    以300m水深8英寸輸油軟管為例,立管共有16層,各層主要參數(shù)如表1所示。

    表1 8英寸輸油軟管各層參數(shù)Table 1 Layer properties of the 8 inch production riser

    其中骨架層旋轉(zhuǎn)角度為88°,抗壓鎧裝層為89°,抗拉鎧裝層角度分別為35°和-35°。

    2 分析方法介紹

    目前針對非粘接柔性立管剛度的計算分析方法主要有兩類: 理論計算法和有限元計算法。

    理論計算法中所依據(jù)的解析方程主要來自柔性軟管設(shè)計和操作手冊[3],此手冊來源于針對柔性軟管開展的研究項目(FPS 2000)成果報告,手冊中提供了軟管剛度計算公式,比較適用于前期設(shè)計階段中軟管剛度的快速計算。由于計算公式基于眾多假設(shè),在后續(xù)迭代計算中需要與有限元計算相結(jié)合,并采用試驗驗證計算結(jié)果。

    對于多層柔性立管的有限元計算,目前有兩種方法: 傳統(tǒng)的等效梁法和多層法[4]。前者采用一個等效的梁單元模型,通過賦予等效截面屬性和材料屬性來進行分析計算[5]。這種方法的優(yōu)點在于模型簡單,節(jié)省計算時間;缺點是無法對立管內(nèi)部各層的受力及相互作用進行計算和模擬。而柔性立管多層模型的有限元計算分析,則可以彌補等效梁方法的不足,可以模擬層與層之間的接觸以及立管層與層間的摩擦和相對滑移,也可以考慮內(nèi)外壓力的影響。

    3 理論計算方法

    3.1 抗拉剛度

    在柔性立管受拉力作用時,主要由螺旋纏繞的抗拉鎧裝層承擔(dān)抗拉作用。聚合物層(內(nèi)襯層、中間包覆層、外包覆層)在拉力作用下對剛度影響不大,在評估柔性立管抗拉性能時忽略其影響。

    假設(shè)螺旋纏繞線為n層,鎧裝層內(nèi)應(yīng)力與拉力關(guān)系計算公式為

    (1)

    式中:ni為i層內(nèi)鎧裝線數(shù)量;σi為i層內(nèi)鎧裝線應(yīng)力;Ai為i層內(nèi)鎧裝線橫截面積;αi為i層內(nèi)鎧裝線角度。

    由于抗壓鎧裝層旋轉(zhuǎn)角通常近似為90°,從式(1)可以看出,抗壓鎧裝層對平衡軸向載荷影響甚微。因此,軸向載荷主要由抗拉鎧裝層承擔(dān)。

    假定拉力由每一條鋼線平均分擔(dān),由軸向載荷引發(fā)的應(yīng)力可由下式計算:

    (2)

    式中:As為單根抗拉鎧裝層鋼線橫截面積,等于寬度與厚度的乘積。

    當柔性立管受到拉力載荷作用時,拉力會造成軸向變形(延伸),內(nèi)外壓會引起徑向變形(橫截面積變大或變小),而橫截面會有旋轉(zhuǎn)趨勢,旋轉(zhuǎn)變形與其他方向的變形相比通常非常小,忽略不計??估z裝層鋼線應(yīng)變與這些變形之間的關(guān)系可以寫作

    (4)

    ν為直徑減小引起的軸向延伸,可以得到

    T=nEAscosα(cos2α-νsin2α)εxx.

    (5)

    因此,柔性立管軸向剛度可以寫作

    EA=nEAscosα(cos2α-νsin2α).

    (6)

    由于抗拉鎧裝層鋼線為多層,表達式可寫作

    (7)

    上式基于小幾何變形的假設(shè),并且假設(shè)在變形過程中各層始終保持接觸。

    抗拉鎧裝層的內(nèi)徑和厚度確定后,抗拉剛度主要由抗拉鎧裝層鋼線決定。由于抗拉層材料的楊氏模量要比熱塑層大得多,因此熱塑層的影響忽略不計。

    3.2 抗扭剛度

    在計算抗扭剛度時,可以同時考慮聚合物層和金屬螺旋層的貢獻。在考慮聚合物層時,其對扭轉(zhuǎn)剛度的貢獻可采用以下公式:

    (8)

    (9)

    (10)

    式中:Gi為某一層的剪切模量;E為楊氏模量;υ為材料泊松比;doi為聚合物層i層外徑;dii為聚合物層i層內(nèi)徑。

    由于金屬螺旋鎧裝層的楊氏模量與聚合物層相比大得多,扭轉(zhuǎn)力矩主要由螺旋鎧裝層承擔(dān),故忽略聚合物層的影響。而骨架層與抗壓鎧裝層作為互鎖層,纏繞角度太大,對抵抗扭矩貢獻不大,因此主要由抗拉鎧裝層抵抗扭矩。

    扭矩與在螺旋鎧裝層中應(yīng)力的平衡公式可寫作

    (11)

    式中:Mt為扭矩;ni為i層中鎧裝線數(shù)量;σi為i層中應(yīng)力;Asi為i層中鎧裝線橫截面積;ri為i層平均半徑;αi為i層旋轉(zhuǎn)角度。

    假設(shè)抗拉鎧裝層中所有鎧裝鋼線的形狀均相同,公式可寫作

    (12)

    忽略軸向和徑向應(yīng)變,應(yīng)力與扭轉(zhuǎn)應(yīng)變之間的關(guān)系為

    (13)

    結(jié)合以上兩式,可以得到

    (14)

    因此,扭轉(zhuǎn)剛度的近似計算公式為

    (15)

    此公式假設(shè)所有層之間始終保持接觸狀態(tài)。

    4 有限元計算法

    在有限元分析軟件的選擇上,選用ABAQUS/Explicit來進行計算[6]。ABAQUS/Explicit求解器可很容易地分析包括多個獨立物體的相互復(fù)雜接觸問題,無需進行迭代就能確定非線性問題的解答,且需要的磁盤空間和內(nèi)存比較小。

    4.1 各層模型及網(wǎng)格劃分

    在建模階段可借助有限元軟件ANSYS,使用APDL語言將各層模型建好之后進行組裝,導(dǎo)入ABAQUS中進行計算[7]。

    骨架層示意圖如圖2所示;根據(jù)剖面圖建模,如圖3所示。

    圖2 骨架層互鎖截面型式Fig.2 Geometrical profile of carcass

    圖3 骨架層有限元模型Fig.3 Finite element model of the carcass

    抗壓鎧裝層橫剖面示意圖如圖4所示;根據(jù)橫剖面建模,如圖5所示。

    圖4 抗壓鎧裝層截面型式Fig.4 Geometrical profile of the pressure armor layer

    圖5 抗壓鎧裝層有限元模型Fig.5 Finite element model of the pressure armor layer

    模型建成后,總體結(jié)構(gòu)的剖面圖如圖6所示。

    圖6 有限元整體模型剖面圖Fig.6 Cross-section profile of the finite element model

    由于包含接觸問題,模型采用三維六面體8節(jié)點非協(xié)調(diào)單元C3D8I,以獲得更好的應(yīng)力解和避免沙漏問題[8]。

    4.2 邊界條件及載荷施加

    在立管兩端中心建立2個幾何參考點[9],然后將兩端所有節(jié)點與中心參考點之間建立運動耦合約束。立管兩端的邊界條件均施加在立管兩端的參考點上,所有載荷也都只施加在立管頂部的參考點上。摩擦系數(shù)取0.1[10]。

    立管頂部所施加的載荷分為拉伸載荷及扭矩,其中扭矩根據(jù)作用方向分兩種情況進行研究。

    5 計算結(jié)果對比與分析

    5.1 抗拉剛度計算結(jié)果

    5.1.1理論計算結(jié)果

    表2為根據(jù)前文公式計算所得聚合物層的抗拉剛度。

    表2 聚合物層的抗拉剛度Table 2 Axial stiffness of the polymer layers

    耐磨層為螺旋纏繞的塑帶層,厚度很小,在剛度計算過程中不予考慮。柔性立管中的金屬螺旋層包含1層抗壓鎧裝層和2層抗拉鎧裝層,其中抗壓鎧裝層螺旋纏繞角度接近90°,在軸向可自由移動,因此抗壓鎧裝層對抗拉剛度貢獻不明顯,主要由抗拉鎧裝層抵抗軸向載荷。根據(jù)公式計算得到抗拉鎧裝層的抗拉剛度,如表3所示。

    表3 抗拉鎧裝層的抗拉剛度Table 3 Axial stiffness of the tensile layers

    從表3可以看出,抗拉鎧裝層的抗拉剛度為659.62MN;而由表2可知,聚合物層的抗拉剛度值為40.11MN,大約只占抗拉鎧裝層剛度值的6%。

    5.1.2有限元計算結(jié)果

    在軸向載荷作用下,模型整體應(yīng)力云圖如圖7所示。

    圖7 軸向載荷作用下整體模型應(yīng)力云圖Fig.7 Axial stress contour

    從圖7可以看出,在軸向載荷作用下,內(nèi)外抗拉鎧裝層的應(yīng)力相對其他各層都要大。

    通過有限元分析,可以得到軸向載荷與軸向應(yīng)變之間的關(guān)系如圖8所示。

    圖8 軸向載荷-應(yīng)變關(guān)系Fig.8 Axial force and axial strain relationship

    圖8中實線代表立管在軸向載荷作用下的應(yīng)變變化,虛線為趨勢擬合線,其斜率代表有限元計算得到的立管軸向抗拉剛度。

    根據(jù)有限元計算結(jié)果可得,立管抗拉剛度值為666.92MN。

    5.1.3結(jié)果對比

    表4為有限元分析結(jié)果與公式計算結(jié)果的對比情況。

    表4 抗拉剛度計算結(jié)果對比

    從表4可以看出,兩者計算結(jié)果相近,誤差為1%,證明有限元計算結(jié)果和理論計算結(jié)果均正確。

    同時可以看出,立管在軸向載荷作用下,抗拉鎧裝層是影響立管軸向抗拉能力的重要層。

    5.2 抗扭剛度計算結(jié)果

    5.2.1理論計算結(jié)果

    對聚合物層的抗扭剛度進行計算,結(jié)果如表5所示。

    表5 聚合物層的抗扭剛度Table 5 Torsional stiffness of the polymer layers

    在柔性立管三層金屬層中包含一層抗壓鎧裝層和兩層抗拉鎧裝層。由于抗壓鎧裝層纏繞角度接近90°,抗壓鎧裝鋼線在環(huán)向可以自由轉(zhuǎn)動,因此在計算抗扭剛度時不考慮其影響,故而柔性立管的抗扭能力主要由抗拉鎧裝層體現(xiàn)??古偠扔嬎愎角拔囊呀?jīng)給出,最終計算得出抗拉鎧裝層的抗扭剛度如表6所示。

    表6 抗拉鎧裝層的抗扭剛度Table 6 Torsional stiffness of the tensile layers

    從以上結(jié)果可以看出,抗拉鎧裝層的扭轉(zhuǎn)剛度值為7 399.63kN·m2,聚合物層的抗扭剛度值為362.0kN·m2,大約只占抗拉鎧裝層扭轉(zhuǎn)剛度的5%。因此,柔性立管在扭矩作用下,抗拉鎧裝層承受主要扭矩載荷。

    5.2.2有限元計算結(jié)果

    抗拉鎧裝層的軸向纏繞角度為35°。沿Z軸扭矩采用線性增加和線性卸載的方式施加在立管自由端的參考點上,方向分別為順時針和逆時針,大小相同。

    在扭矩作用下,模型整體應(yīng)力云圖如圖9所示。

    從圖9可以看出,在扭矩作用下,內(nèi)外抗拉鎧裝層應(yīng)力較其他各層更大,說明主要由抗拉鎧裝層承受扭矩作用。

    根據(jù)有限元計算,得到扭矩與扭轉(zhuǎn)曲率之間的關(guān)系圖,如圖10所示。

    圖10中曲線斜率代表抗扭剛度,順時針抗扭剛度值為8032kN·m2,逆時針抗扭剛度值為10347kN·m2。逆時針剛度比順時針剛度大,主要

    (a) 逆時針扭矩

    (b) 順時針扭矩圖9 扭矩作用下模型整體應(yīng)力云圖Fig.9 Torsional stress contour

    圖10 扭矩與扭轉(zhuǎn)曲率之間的變化關(guān)系Fig.10 Relationship between torque and twisting curvature

    是因為逆時針扭矩方向與外抗拉鎧裝層旋轉(zhuǎn)方向相同,而外抗拉鎧裝層的半徑要比內(nèi)抗拉鎧裝層大。

    5.2.3結(jié)果對比

    將有限元分析結(jié)果與公式計算結(jié)果進行對比,如表7所示。

    表7 不同方法得到的抗扭剛度值

    從表7可以看出,三者在數(shù)量級上一致,相比較而言理論公式計算值更為保守,而有限元計算結(jié)果偏大,且不同方向的抗扭剛度不同,逆時針抗扭剛度值較大,主要是由于此時扭矩方向與外抗拉鎧裝層旋轉(zhuǎn)方向相同,也說明扭矩作用方向?qū)εまD(zhuǎn)剛度的計算影響比較大,在今后的研究中應(yīng)注意扭矩的作用方向與抗拉鎧裝層中螺旋鋼線的旋轉(zhuǎn)方向之間的關(guān)系。

    6 結(jié) 語

    根據(jù)理論計算公式和ABAQUS軟件有限元分析,針對300m水深8英寸輸油軟管,分別對柔性立管抗拉和抗扭剛度進行計算并對比,最終得出以下結(jié)論:

    (1) 通過理論公式計算結(jié)果與有限元計算結(jié)果的對比,可以看出柔性立管在軸向載荷和扭矩作用下,理論計算結(jié)果與有限元結(jié)算結(jié)果較為一致,在抗扭剛度計算中,理論計算結(jié)果相對有限元計算結(jié)果更偏保守,在今后的柔性立管設(shè)計過程中兩種方法均可使用。

    (2) 由于理論計算法簡便快速,可以在前期設(shè)計階段使用,以便節(jié)省時間;有限元計算法可以在詳細設(shè)計階段用于校核剛度計算結(jié)果,但最終仍需要通過試驗進行驗證。

    (3) 柔性立管抗拉剛度和抗扭剛度主要由抗拉鎧裝層貢獻,聚合物層影響不大。抗拉鎧裝層的寬度、厚度和纏繞角度是主要影響參數(shù)。

    [1] Brack M, Troina L M B, Sousa J R M D. Flexible riser resistance against combined axial compression, bending, and torsion in ultra-deep water depths [C]. OMAE, 2005: 67404.

    [2] Kagoura T, Ishii K I, Abe S, et al. Development of a flexible pipe for pipe-in-pipe technology [J]. Furukawa Review, 2003(24): 69.

    [3] Berge S, Engseth A, Fylling I, et al. Handbook on design and operation of flexible pipes [M]. Trondheim: SINTEF, 1992.

    [4] 陳耕.海洋非粘結(jié)柔性管線的橫截面力學(xué)性能分析[D].哈爾濱: 哈爾濱工程大學(xué),2010.

    Chen Geng. Analysis on cross-section mechanical properties of non-bonded offshore flexible pipe [D]. Harbin: Harbin Engineering University, 2010.

    [5] 欒林昌.非粘合柔性立管力學(xué)性能研究[D].大連: 大連理工大學(xué),2011.

    Luan Lin-chang. Research of mechanical behavior of unbonded flexible risers [D]. Dalian: Dalian University of Technology, 2011.

    [6] Perdrizet T, Leroy J-M, Barbin N, et al. Stresses in armour layers of flexible pipes: comparison of ABAQUS models [C].

    SIMULIA Customer Conference, 2011.

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

    Wang Xin-min. Numerical analysis of engineering structure with ANSYS [M]. Beijing: China Communications Press, 2007.

    [8] 莊茁,張帆,岑松,等.ABAQUS非線性有限元分析實例[M].北京: 科學(xué)出版社,2005.

    Zhuang Zhuo, Zhang Fan, Cen Song, et al. ABAQUS nonlinear finite element analysis cases [M]. Beijing: Science Press, 2005.

    [9] 石亦平,周玉蓉.ABAQUS有限元分析實例詳解[M].北京: 機械工業(yè)出版社,2006.

    Shi Yi-ping, Zhou Yu-rong. Detailed notes to ABAQUS finite element analysis examples [M]. Beijing: China Machine Press, 2006.

    ComparativeAnalysisofCalculationMethodsforAxialandTorsionalStiffnessofUnbondedFlexibleRiser

    PEI Xiao-mei, ZHANG En-yong, LI Li-wei, ZHOU Wei-wei, CAO Jing

    (CNOOCResearchInstitute,Beijing100028,China)

    We present the study of unbonded flexible pipe and provide a case of 8 inch unbonded flexible riser at the water depth of 300m. Firstly, the theoretical calculation formulae are introduced for the analysis of axial and torsional stiffness of the unbonded flexible riser. On the basis of the theoretical calculation, finite element analysis is carried out by modelling all layers separately to reflect the mechanical behavior of each layer. The numerical results are compared with the theoretical results for tensile load and torsion. The different contributions of each layer are discussed. Under the axial load or torque, the theoretical results are consistent with the numerical results, and the theoretical results of torsional stiffness are more conservative. During the actual engineering design, the appropriate calculation method can be chosen according to the design stage.

    unbonded flexible riser; axial stiffness; torsional stiffness; contact

    2017-06-15

    國家科技重大專項(2016ZX05028-005);京直地區(qū)第四屆青年科技與管理創(chuàng)新研究課題(JZTW2017KJ26)

    裴曉梅(1988—),女,碩士,工程師,主要從事海底管道和立管結(jié)構(gòu)方面的研究。

    TE973.1

    A

    2095-7297(2017)05-0307-08

    猜你喜歡
    鎧裝抗拉立管
    重組竹抗拉力學(xué)性能分析
    改性聚酯纖維耐堿性能(抗拉強力保持率)測量不確定度評定
    聯(lián)鎖鎧裝電纜的中外標準解析
    電線電纜(2021年1期)2021-03-29 02:26:42
    淺談鎧裝電纜引入防爆電氣設(shè)備
    防爆電機(2020年6期)2020-12-14 07:17:02
    常見高層建筑物室內(nèi)給水立管材質(zhì)解析
    分析10kV金屬封閉鎧裝式開關(guān)柜內(nèi)部發(fā)熱的原因
    電子制作(2018年14期)2018-08-21 01:38:32
    抗拉柔性鉸鏈的理論建模及有限元分析
    深水鋼懸鏈立管J型鋪設(shè)研究
    一種自制鎧裝電纜裝置及制作方法介紹
    水利信息化(2015年5期)2015-12-21 12:54:40
    The Power of Integration
    Beijing Review(2015年43期)2015-11-25 03:12:04
    久久精品国产清高在天天线| 日日夜夜操网爽| 亚洲性久久影院| 国产高潮美女av| 成人三级黄色视频| 偷拍熟女少妇极品色| 精品欧美国产一区二区三| 国内久久婷婷六月综合欲色啪| 欧美人与善性xxx| 动漫黄色视频在线观看| 日本五十路高清| 婷婷六月久久综合丁香| 色综合站精品国产| 一个人观看的视频www高清免费观看| aaaaa片日本免费| 国产欧美日韩精品一区二区| 久久欧美精品欧美久久欧美| 欧美性猛交╳xxx乱大交人| 18禁黄网站禁片午夜丰满| 欧美一级a爱片免费观看看| 亚洲人成网站在线播放欧美日韩| 国产免费av片在线观看野外av| 欧美激情在线99| 免费看日本二区| 欧美日本视频| 神马国产精品三级电影在线观看| 少妇被粗大猛烈的视频| 国产免费av片在线观看野外av| 久久精品国产亚洲av香蕉五月| 亚洲国产欧美人成| 噜噜噜噜噜久久久久久91| 久久久久久大精品| 我要看日韩黄色一级片| 欧美日韩乱码在线| 丝袜美腿在线中文| 免费大片18禁| 制服丝袜大香蕉在线| 午夜福利在线观看吧| 特大巨黑吊av在线直播| 国产黄片美女视频| 婷婷六月久久综合丁香| 国内精品宾馆在线| 亚洲av中文字字幕乱码综合| 国产男靠女视频免费网站| a级毛片a级免费在线| 国产精品一区www在线观看 | 99久久精品热视频| 国产亚洲av嫩草精品影院| 国产在视频线在精品| 在线a可以看的网站| 一区福利在线观看| 春色校园在线视频观看| 国产精品电影一区二区三区| 日韩欧美国产一区二区入口| 亚洲av一区综合| 欧美色欧美亚洲另类二区| 久久久久久久久中文| 日本免费一区二区三区高清不卡| 舔av片在线| 欧美色欧美亚洲另类二区| 久久久久免费精品人妻一区二区| 少妇被粗大猛烈的视频| 亚洲欧美清纯卡通| 精品久久久噜噜| 国产激情偷乱视频一区二区| 午夜免费成人在线视频| 欧美激情在线99| 蜜桃久久精品国产亚洲av| 欧美不卡视频在线免费观看| 久久亚洲精品不卡| 日本一二三区视频观看| 成人美女网站在线观看视频| 又粗又爽又猛毛片免费看| 露出奶头的视频| 久久久久性生活片| 又黄又爽又免费观看的视频| 亚洲综合色惰| 午夜精品久久久久久毛片777| 联通29元200g的流量卡| av在线观看视频网站免费| 欧美性感艳星| 亚洲图色成人| av在线蜜桃| 在线免费观看的www视频| 禁无遮挡网站| 免费av观看视频| 人人妻,人人澡人人爽秒播| 极品教师在线免费播放| 免费看av在线观看网站| 久久国产精品人妻蜜桃| 国产精品人妻久久久久久| 看十八女毛片水多多多| 日本-黄色视频高清免费观看| 国产高潮美女av| 国产三级中文精品| 美女免费视频网站| 亚洲国产欧洲综合997久久,| 午夜爱爱视频在线播放| 欧美精品啪啪一区二区三区| 亚洲一级一片aⅴ在线观看| 日本一二三区视频观看| 高清在线国产一区| 中文字幕av在线有码专区| 白带黄色成豆腐渣| 成人av一区二区三区在线看| 搡老岳熟女国产| 国产伦精品一区二区三区视频9| 天堂动漫精品| 男人的好看免费观看在线视频| 欧美一级a爱片免费观看看| 国内少妇人妻偷人精品xxx网站| 久久精品综合一区二区三区| 欧美丝袜亚洲另类 | 一个人免费在线观看电影| 欧美极品一区二区三区四区| 日本黄大片高清| 国产一区二区亚洲精品在线观看| 精品日产1卡2卡| 一级黄色大片毛片| a级毛片免费高清观看在线播放| 成年免费大片在线观看| 人妻制服诱惑在线中文字幕| 久久久久久久久久黄片| 中文字幕精品亚洲无线码一区| 日日摸夜夜添夜夜添av毛片 | 天堂影院成人在线观看| 日韩一本色道免费dvd| 极品教师在线免费播放| 国产大屁股一区二区在线视频| 永久网站在线| 制服丝袜大香蕉在线| 精品99又大又爽又粗少妇毛片 | 一进一出抽搐动态| 久久久久国产精品人妻aⅴ院| 亚洲美女视频黄频| 国产淫片久久久久久久久| 精品久久久久久久人妻蜜臀av| 午夜影院日韩av| 免费搜索国产男女视频| 国产精品野战在线观看| 久久精品夜夜夜夜夜久久蜜豆| 春色校园在线视频观看| 午夜日韩欧美国产| 久久99热这里只有精品18| 日本欧美国产在线视频| 免费看光身美女| 免费av毛片视频| av女优亚洲男人天堂| av在线亚洲专区| 久久精品夜夜夜夜夜久久蜜豆| 成年人黄色毛片网站| 小蜜桃在线观看免费完整版高清| 日本 av在线| 午夜精品在线福利| 亚洲国产精品合色在线| 少妇的逼好多水| 日本撒尿小便嘘嘘汇集6| 少妇熟女aⅴ在线视频| 一进一出好大好爽视频| av女优亚洲男人天堂| 欧美+日韩+精品| 日韩欧美一区二区三区在线观看| 国内精品久久久久久久电影| 亚洲成人免费电影在线观看| 美女xxoo啪啪120秒动态图| 亚洲av电影不卡..在线观看| 欧美+亚洲+日韩+国产| 成人午夜高清在线视频| 国产av麻豆久久久久久久| 美女大奶头视频| 我要看日韩黄色一级片| 久9热在线精品视频| 精品久久久久久久久久久久久| 国产探花极品一区二区| 色在线成人网| АⅤ资源中文在线天堂| 日本 欧美在线| 欧美极品一区二区三区四区| 亚洲国产精品合色在线| 22中文网久久字幕| 国产精品久久久久久精品电影| 69人妻影院| 亚洲美女搞黄在线观看 | 十八禁网站免费在线| 很黄的视频免费| 国产乱人伦免费视频| 国产av不卡久久| 一级av片app| 女同久久另类99精品国产91| 国产精品爽爽va在线观看网站| 久久精品国产亚洲av天美| 一进一出抽搐gif免费好疼| 我要搜黄色片| 啦啦啦观看免费观看视频高清| 欧美日本亚洲视频在线播放| 国内精品宾馆在线| 国语自产精品视频在线第100页| 美女高潮喷水抽搐中文字幕| 亚洲av不卡在线观看| 欧美日韩黄片免| 国产精品女同一区二区软件 | 欧美日韩中文字幕国产精品一区二区三区| 国产av一区在线观看免费| 桃色一区二区三区在线观看| av天堂中文字幕网| 一级毛片久久久久久久久女| 欧美最黄视频在线播放免费| 男女做爰动态图高潮gif福利片| 在线免费十八禁| 国内精品一区二区在线观看| 久久久久久大精品| 黄色欧美视频在线观看| 国产伦精品一区二区三区视频9| 国产成人av教育| 国产乱人伦免费视频| 国产一区二区亚洲精品在线观看| 亚洲不卡免费看| 中文字幕av成人在线电影| 欧美在线一区亚洲| 国产伦精品一区二区三区四那| 一区福利在线观看| 成人鲁丝片一二三区免费| 动漫黄色视频在线观看| 超碰av人人做人人爽久久| 老女人水多毛片| 国产午夜精品久久久久久一区二区三区 | 免费高清视频大片| 日本三级黄在线观看| 亚洲av中文av极速乱 | 91午夜精品亚洲一区二区三区 | 国产高潮美女av| 丝袜美腿在线中文| 欧美+日韩+精品| 国产欧美日韩精品亚洲av| 一卡2卡三卡四卡精品乱码亚洲| 悠悠久久av| 91久久精品国产一区二区成人| 国产精品电影一区二区三区| 国产精品av视频在线免费观看| 国产真实乱freesex| 99久久精品一区二区三区| 悠悠久久av| 国产一区二区亚洲精品在线观看| 乱系列少妇在线播放| 欧美潮喷喷水| 国产高潮美女av| 在线免费观看的www视频| 高清毛片免费观看视频网站| 精品99又大又爽又粗少妇毛片 | 久久99热6这里只有精品| av.在线天堂| 一个人免费在线观看电影| 内地一区二区视频在线| 欧美精品啪啪一区二区三区| 欧美日韩精品成人综合77777| 国产精品自产拍在线观看55亚洲| 午夜爱爱视频在线播放| 国产av麻豆久久久久久久| 午夜影院日韩av| 一本一本综合久久| 又紧又爽又黄一区二区| 男女视频在线观看网站免费| 韩国av一区二区三区四区| 久久久久久久久久黄片| 精品久久久久久久久久久久久| 精品人妻视频免费看| 简卡轻食公司| 我的老师免费观看完整版| 91久久精品国产一区二区成人| 亚洲精华国产精华液的使用体验 | 久久国产乱子免费精品| 两人在一起打扑克的视频| 成年女人看的毛片在线观看| 在线观看午夜福利视频| 欧美激情久久久久久爽电影| 欧美一区二区亚洲| 久久久久久久久中文| 亚洲自偷自拍三级| 老女人水多毛片| 久久精品久久久久久噜噜老黄 | 国产高清激情床上av| 国产亚洲欧美98| 成人毛片a级毛片在线播放| 日本 欧美在线| 特级一级黄色大片| av在线天堂中文字幕| 99热精品在线国产| 国产精品久久久久久精品电影| 蜜桃亚洲精品一区二区三区| 国产一区二区三区av在线 | 99久久精品国产国产毛片| 白带黄色成豆腐渣| 69av精品久久久久久| 国产午夜精品论理片| 少妇的逼好多水| 中文在线观看免费www的网站| av国产免费在线观看| 精品无人区乱码1区二区| 亚洲中文字幕一区二区三区有码在线看| 免费人成视频x8x8入口观看| 伊人久久精品亚洲午夜| 在线观看美女被高潮喷水网站| 亚洲在线观看片| 搡女人真爽免费视频火全软件 | 欧美另类亚洲清纯唯美| 999久久久精品免费观看国产| 舔av片在线| 99热精品在线国产| 国产三级在线视频| 禁无遮挡网站| 亚洲第一区二区三区不卡| 欧美成人a在线观看| 成人永久免费在线观看视频| 琪琪午夜伦伦电影理论片6080| 精品日产1卡2卡| 国产午夜精品论理片| 在线播放国产精品三级| 国产探花极品一区二区| 97人妻精品一区二区三区麻豆| 一区二区三区四区激情视频 | 在线观看午夜福利视频| 日本成人三级电影网站| 性色avwww在线观看| 白带黄色成豆腐渣| 日日干狠狠操夜夜爽| 91在线观看av| 久久精品国产自在天天线| 一本久久中文字幕| 国产av麻豆久久久久久久| 五月伊人婷婷丁香| 亚洲va在线va天堂va国产| 少妇的逼水好多| 少妇丰满av| 97热精品久久久久久| 中文字幕久久专区| 国产av在哪里看| 亚洲av第一区精品v没综合| 久久精品国产亚洲av天美| 亚洲中文字幕一区二区三区有码在线看| 亚洲aⅴ乱码一区二区在线播放| 欧美激情国产日韩精品一区| 舔av片在线| 午夜免费男女啪啪视频观看 | 99精品久久久久人妻精品| 亚洲电影在线观看av| 亚洲va在线va天堂va国产| 三级男女做爰猛烈吃奶摸视频| 69av精品久久久久久| 日韩,欧美,国产一区二区三区 | 欧美区成人在线视频| 国产精品久久视频播放| 美女高潮喷水抽搐中文字幕| 99久久成人亚洲精品观看| 男女下面进入的视频免费午夜| 久久精品国产99精品国产亚洲性色| 久久欧美精品欧美久久欧美| 欧美xxxx性猛交bbbb| 三级国产精品欧美在线观看| 熟女人妻精品中文字幕| 男人舔奶头视频| 国产麻豆成人av免费视频| 91久久精品电影网| 非洲黑人性xxxx精品又粗又长| 又粗又爽又猛毛片免费看| 18禁黄网站禁片午夜丰满| 久久热精品热| xxxwww97欧美| 搡老岳熟女国产| 老司机午夜福利在线观看视频| 国产午夜精品久久久久久一区二区三区 | 日韩高清综合在线| 免费观看在线日韩| 久久久国产成人免费| 欧美区成人在线视频| 最后的刺客免费高清国语| 国产亚洲精品久久久久久毛片| 午夜免费激情av| 99久久九九国产精品国产免费| 有码 亚洲区| 免费看av在线观看网站| 成年女人毛片免费观看观看9| 中文字幕久久专区| 久久久久久久久大av| 桃色一区二区三区在线观看| 小说图片视频综合网站| 91久久精品国产一区二区三区| 精品久久久久久久久久久久久| 人人妻人人澡欧美一区二区| 简卡轻食公司| 日本熟妇午夜| 午夜影院日韩av| 亚洲五月天丁香| 日韩在线高清观看一区二区三区 | 免费高清视频大片| 日日夜夜操网爽| 久久精品综合一区二区三区| 色5月婷婷丁香| 精品免费久久久久久久清纯| 亚洲经典国产精华液单| www.www免费av| 亚洲电影在线观看av| 亚洲欧美日韩高清专用| 波多野结衣高清无吗| x7x7x7水蜜桃| 国产精品综合久久久久久久免费| 亚洲五月天丁香| 波多野结衣巨乳人妻| 久久久久久久久大av| 老司机午夜福利在线观看视频| 99久久久亚洲精品蜜臀av| 两个人视频免费观看高清| 日韩一区二区视频免费看| 五月伊人婷婷丁香| 欧美色视频一区免费| 欧美精品啪啪一区二区三区| 国产私拍福利视频在线观看| 麻豆av噜噜一区二区三区| 不卡一级毛片| 一本久久中文字幕| 国产精品综合久久久久久久免费| 麻豆精品久久久久久蜜桃| 91av网一区二区| 亚洲久久久久久中文字幕| 亚洲综合色惰| av黄色大香蕉| 久久精品国产自在天天线| 热99re8久久精品国产| 在线天堂最新版资源| 日本爱情动作片www.在线观看 | 窝窝影院91人妻| 久久亚洲精品不卡| 99热这里只有是精品在线观看| 成人特级黄色片久久久久久久| 乱人视频在线观看| av在线老鸭窝| 少妇丰满av| 桃色一区二区三区在线观看| 亚洲成人久久性| 久久久久国内视频| 亚洲最大成人中文| 久久久久久国产a免费观看| 99热精品在线国产| 免费看美女性在线毛片视频| 亚洲av五月六月丁香网| 桃红色精品国产亚洲av| 嫩草影视91久久| 一级av片app| 国产伦精品一区二区三区视频9| 国产高清三级在线| 国国产精品蜜臀av免费| 国产精品乱码一区二三区的特点| a在线观看视频网站| 动漫黄色视频在线观看| 国产成人影院久久av| 又黄又爽又免费观看的视频| 九九热线精品视视频播放| 国产午夜精品论理片| 亚洲黑人精品在线| 夜夜看夜夜爽夜夜摸| 人妻少妇偷人精品九色| av在线观看视频网站免费| 色av中文字幕| 我的女老师完整版在线观看| 免费观看精品视频网站| 国产精品自产拍在线观看55亚洲| 久久久久国产精品人妻aⅴ院| 18禁裸乳无遮挡免费网站照片| 国产高清视频在线播放一区| 日本黄色片子视频| 国产伦精品一区二区三区视频9| 午夜精品一区二区三区免费看| 狂野欧美激情性xxxx在线观看| 亚洲最大成人中文| a在线观看视频网站| 精品人妻偷拍中文字幕| 非洲黑人性xxxx精品又粗又长| 亚洲精品成人久久久久久| 天堂影院成人在线观看| 波多野结衣高清作品| 赤兔流量卡办理| 两个人视频免费观看高清| 日韩中文字幕欧美一区二区| 国产成人av教育| 日韩中字成人| 中文在线观看免费www的网站| 国产精品日韩av在线免费观看| 内射极品少妇av片p| 午夜老司机福利剧场| 中文字幕av成人在线电影| 亚洲欧美精品综合久久99| 97超视频在线观看视频| 少妇的逼水好多| 亚洲aⅴ乱码一区二区在线播放| 动漫黄色视频在线观看| 午夜激情福利司机影院| 一a级毛片在线观看| 久久久精品大字幕| 丰满乱子伦码专区| 亚洲国产精品合色在线| 久久午夜亚洲精品久久| 欧美日本亚洲视频在线播放| 高清在线国产一区| 色哟哟·www| 熟妇人妻久久中文字幕3abv| 日韩欧美精品免费久久| 欧美+亚洲+日韩+国产| 精品久久久久久成人av| 国产国拍精品亚洲av在线观看| 日日摸夜夜添夜夜添av毛片 | 老熟妇仑乱视频hdxx| 欧美中文日本在线观看视频| 久久精品国产清高在天天线| 人人妻人人看人人澡| 毛片一级片免费看久久久久 | 一本久久中文字幕| 神马国产精品三级电影在线观看| 亚洲不卡免费看| 他把我摸到了高潮在线观看| 天堂影院成人在线观看| 亚洲综合色惰| 国产在视频线在精品| 亚洲av.av天堂| 精品久久久久久久久久久久久| 丰满的人妻完整版| 国产在视频线在精品| 免费av毛片视频| 国产精品久久电影中文字幕| 一区福利在线观看| 日韩欧美在线二视频| 精品国产三级普通话版| 在线播放国产精品三级| 国产精品一及| 十八禁国产超污无遮挡网站| 97超视频在线观看视频| 国内精品久久久久精免费| 亚洲熟妇熟女久久| 午夜视频国产福利| 琪琪午夜伦伦电影理论片6080| 日日摸夜夜添夜夜添小说| 日日摸夜夜添夜夜添小说| 女生性感内裤真人,穿戴方法视频| 国产精品亚洲美女久久久| 春色校园在线视频观看| 国产午夜精品久久久久久一区二区三区 | 99精品久久久久人妻精品| 久久精品国产自在天天线| 欧美潮喷喷水| 男女视频在线观看网站免费| 久久精品人妻少妇| 久久精品国产鲁丝片午夜精品 | 国产乱人伦免费视频| 亚洲人与动物交配视频| 热99在线观看视频| 少妇的逼水好多| 国模一区二区三区四区视频| 日韩欧美 国产精品| 亚洲人成网站在线播| 国产国拍精品亚洲av在线观看| 欧美日韩瑟瑟在线播放| 能在线免费观看的黄片| 成人国产麻豆网| 精品久久久久久久人妻蜜臀av| 久久天躁狠狠躁夜夜2o2o| 12—13女人毛片做爰片一| 国产成人av教育| 自拍偷自拍亚洲精品老妇| 色视频www国产| 亚洲真实伦在线观看| 男人狂女人下面高潮的视频| 欧美一区二区精品小视频在线| 99在线视频只有这里精品首页| 精品人妻偷拍中文字幕| 美女高潮的动态| 久久国产乱子免费精品| 国产精品久久电影中文字幕| 亚洲av.av天堂| 国产欧美日韩精品一区二区| 超碰av人人做人人爽久久| 免费看光身美女| 色精品久久人妻99蜜桃| 日韩欧美 国产精品| 欧洲精品卡2卡3卡4卡5卡区| 亚洲性夜色夜夜综合| 日本黄大片高清| 99热这里只有精品一区| 中文字幕熟女人妻在线| 国产极品精品免费视频能看的| 欧美极品一区二区三区四区| a级毛片免费高清观看在线播放| 九九在线视频观看精品| 亚洲一区高清亚洲精品| 国产亚洲91精品色在线| 国产乱人伦免费视频| 两人在一起打扑克的视频| 亚洲熟妇中文字幕五十中出| 一个人看的www免费观看视频| 精品无人区乱码1区二区| 美女高潮喷水抽搐中文字幕| 婷婷亚洲欧美| 欧美色欧美亚洲另类二区| 亚洲国产欧洲综合997久久,| 国产黄片美女视频| 啦啦啦啦在线视频资源| 亚洲av成人av| 一级黄色大片毛片| 久久草成人影院| 琪琪午夜伦伦电影理论片6080| 亚洲成人免费电影在线观看| 亚洲av不卡在线观看| 中文字幕熟女人妻在线| 精品一区二区三区人妻视频| 欧美国产日韩亚洲一区| 听说在线观看完整版免费高清| 免费在线观看影片大全网站| 别揉我奶头 嗯啊视频|