• <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
    露出奶头的视频| 午夜福利在线在线| 91在线观看av| 校园春色视频在线观看| 十八禁人妻一区二区| 国产欧美日韩精品亚洲av| 一区二区三区高清视频在线| 无限看片的www在线观看| 久久久精品大字幕| 在线观看日韩欧美| 国产精品爽爽va在线观看网站| 美女免费视频网站| av在线蜜桃| 在线十欧美十亚洲十日本专区| 国产综合懂色| 国产又黄又爽又无遮挡在线| 一边摸一边抽搐一进一小说| 亚洲熟妇熟女久久| 五月伊人婷婷丁香| 欧美三级亚洲精品| 亚洲成av人片在线播放无| 狂野欧美激情性xxxx| 最近最新免费中文字幕在线| 午夜免费观看网址| 午夜福利高清视频| 成人精品一区二区免费| 一级黄片播放器| 嫩草影院入口| 天堂动漫精品| 亚洲av二区三区四区| 法律面前人人平等表现在哪些方面| 欧美黑人巨大hd| 一进一出好大好爽视频| 亚洲中文字幕日韩| 欧美日本视频| 欧美中文综合在线视频| 高清毛片免费观看视频网站| 一级a爱片免费观看的视频| 亚洲av二区三区四区| 中文字幕人妻熟人妻熟丝袜美 | 久久久久久大精品| 人人妻人人澡欧美一区二区| 少妇高潮的动态图| 精品午夜福利视频在线观看一区| 成人无遮挡网站| 给我免费播放毛片高清在线观看| eeuss影院久久| 小说图片视频综合网站| 看黄色毛片网站| 欧美另类亚洲清纯唯美| 99在线视频只有这里精品首页| 极品教师在线免费播放| 亚洲,欧美精品.| 欧美在线黄色| 日韩欧美 国产精品| 色尼玛亚洲综合影院| 日日干狠狠操夜夜爽| 亚洲精品一区av在线观看| 欧美日韩精品网址| 久久久久久久亚洲中文字幕 | 18禁在线播放成人免费| 国产蜜桃级精品一区二区三区| 亚洲av成人不卡在线观看播放网| 国产精品乱码一区二三区的特点| 国产色爽女视频免费观看| 免费在线观看成人毛片| 中亚洲国语对白在线视频| a级毛片a级免费在线| 久久久久久久午夜电影| 啦啦啦免费观看视频1| 91久久精品电影网| 天天一区二区日本电影三级| 美女黄网站色视频| 人妻丰满熟妇av一区二区三区| 日韩 欧美 亚洲 中文字幕| 国产午夜精品论理片| 色av中文字幕| 我要搜黄色片| 国产精品av视频在线免费观看| 国产成人a区在线观看| 真人一进一出gif抽搐免费| 亚洲18禁久久av| 日本与韩国留学比较| 国产精品1区2区在线观看.| 天堂网av新在线| 女人被狂操c到高潮| 欧美成狂野欧美在线观看| 免费大片18禁| 丰满乱子伦码专区| 久久久久久久久久黄片| 精华霜和精华液先用哪个| 精品久久久久久,| 身体一侧抽搐| 亚洲精品一卡2卡三卡4卡5卡| 中文字幕高清在线视频| 亚洲第一电影网av| 亚洲最大成人手机在线| 国产综合懂色| 国产美女午夜福利| tocl精华| 久久久久国内视频| 狠狠狠狠99中文字幕| 琪琪午夜伦伦电影理论片6080| 亚洲精品色激情综合| 日韩欧美一区二区三区在线观看| 国产一区在线观看成人免费| 日韩欧美国产一区二区入口| 悠悠久久av| 日韩中文字幕欧美一区二区| 天天一区二区日本电影三级| 亚洲一区二区三区不卡视频| 午夜视频国产福利| 免费观看精品视频网站| 久久久久久久亚洲中文字幕 | 日本一二三区视频观看| 丝袜美腿在线中文| 香蕉丝袜av| 男插女下体视频免费在线播放| 久久性视频一级片| 桃红色精品国产亚洲av| 久久亚洲真实| 夜夜爽天天搞| 日韩av在线大香蕉| 啪啪无遮挡十八禁网站| 欧美黑人巨大hd| 嫩草影院入口| 久久天躁狠狠躁夜夜2o2o| а√天堂www在线а√下载| 亚洲av电影在线进入| 成人特级av手机在线观看| 成人特级av手机在线观看| 国产一区二区亚洲精品在线观看| 中亚洲国语对白在线视频| 亚洲黑人精品在线| 免费看美女性在线毛片视频| 日韩欧美免费精品| 深爱激情五月婷婷| 国产av不卡久久| 五月玫瑰六月丁香| 精品乱码久久久久久99久播| 岛国视频午夜一区免费看| 97超视频在线观看视频| 宅男免费午夜| 亚洲成人中文字幕在线播放| 日韩欧美 国产精品| av在线蜜桃| 亚洲成人免费电影在线观看| 日韩精品中文字幕看吧| 在线观看舔阴道视频| 窝窝影院91人妻| 首页视频小说图片口味搜索| 可以在线观看的亚洲视频| 在线观看免费午夜福利视频| 在线观看美女被高潮喷水网站 | 国产97色在线日韩免费| 女警被强在线播放| 国产成人aa在线观看| 在线看三级毛片| 成年女人毛片免费观看观看9| 午夜福利成人在线免费观看| 少妇的丰满在线观看| 亚洲精品乱码久久久v下载方式 | 国产野战对白在线观看| 婷婷精品国产亚洲av在线| av中文乱码字幕在线| 人人妻,人人澡人人爽秒播| 2021天堂中文幕一二区在线观| 午夜福利18| 欧美日本亚洲视频在线播放| 婷婷精品国产亚洲av在线| 亚洲欧美日韩无卡精品| 美女大奶头视频| 色在线成人网| eeuss影院久久| 日本精品一区二区三区蜜桃| 国产黄a三级三级三级人| 最近视频中文字幕2019在线8| 中文字幕人妻丝袜一区二区| av中文乱码字幕在线| 国产欧美日韩一区二区精品| 伊人久久精品亚洲午夜| 国内毛片毛片毛片毛片毛片| 久久久久久久久久黄片| 美女高潮喷水抽搐中文字幕| 深夜精品福利| 老汉色av国产亚洲站长工具| 欧美黑人欧美精品刺激| 99精品欧美一区二区三区四区| 国产私拍福利视频在线观看| 天堂动漫精品| 免费高清视频大片| 九色成人免费人妻av| 国产一区二区在线av高清观看| 亚洲精品粉嫩美女一区| 国产一区二区三区在线臀色熟女| 欧美日韩乱码在线| 观看美女的网站| 乱人视频在线观看| 偷拍熟女少妇极品色| 亚洲性夜色夜夜综合| 国产三级黄色录像| 哪里可以看免费的av片| 19禁男女啪啪无遮挡网站| 少妇的逼好多水| 久久久国产精品麻豆| 中文字幕熟女人妻在线| а√天堂www在线а√下载| 亚洲av熟女| 免费人成视频x8x8入口观看| 成人无遮挡网站| 真人一进一出gif抽搐免费| 国产成人av教育| www国产在线视频色| 成熟少妇高潮喷水视频| 亚洲男人的天堂狠狠| 国产高清videossex| 国内精品久久久久精免费| 久久久久免费精品人妻一区二区| 国产毛片a区久久久久| 女同久久另类99精品国产91| 人妻夜夜爽99麻豆av| 午夜免费成人在线视频| 亚洲片人在线观看| 深夜精品福利| 有码 亚洲区| 国产欧美日韩一区二区精品| 欧美另类亚洲清纯唯美| 精品人妻一区二区三区麻豆 | 深夜精品福利| 国产精品日韩av在线免费观看| 淫秽高清视频在线观看| 久久6这里有精品| 国产精华一区二区三区| av女优亚洲男人天堂| 别揉我奶头~嗯~啊~动态视频| 国产亚洲欧美98| 亚洲不卡免费看| 国内精品一区二区在线观看| 亚洲国产精品合色在线| 女同久久另类99精品国产91| 我的老师免费观看完整版| 久久人人精品亚洲av| www.色视频.com| 久久精品综合一区二区三区| 69av精品久久久久久| 久久性视频一级片| 亚洲一区高清亚洲精品| xxx96com| 又黄又粗又硬又大视频| 日韩中文字幕欧美一区二区| 久久久久久人人人人人| 亚洲精品日韩av片在线观看 | 亚洲av第一区精品v没综合| 丝袜美腿在线中文| 别揉我奶头~嗯~啊~动态视频| 色视频www国产| 国产高清三级在线| 99久久99久久久精品蜜桃| 成年版毛片免费区| 好男人电影高清在线观看| 中文字幕人成人乱码亚洲影| 亚洲av免费在线观看| 日本免费a在线| 日韩高清综合在线| 51国产日韩欧美| 亚洲av不卡在线观看| av专区在线播放| 国产乱人伦免费视频| 啦啦啦免费观看视频1| 少妇的逼好多水| 最近最新中文字幕大全免费视频| 国产精华一区二区三区| 欧美黄色淫秽网站| av黄色大香蕉| 看片在线看免费视频| 欧美日韩精品网址| 两个人的视频大全免费| 午夜福利欧美成人| 成年女人永久免费观看视频| 一个人观看的视频www高清免费观看| 一个人观看的视频www高清免费观看| 国产淫片久久久久久久久 | 男人舔奶头视频| 丰满人妻一区二区三区视频av | 99国产精品一区二区蜜桃av| av女优亚洲男人天堂| 欧美色欧美亚洲另类二区| 观看免费一级毛片| 色老头精品视频在线观看| 午夜激情福利司机影院| 久久草成人影院| 免费在线观看日本一区| 岛国在线观看网站| 成人国产一区最新在线观看| 国产欧美日韩一区二区三| 69人妻影院| 日本一二三区视频观看| 又爽又黄无遮挡网站| 啦啦啦韩国在线观看视频| 久久国产精品影院| 亚洲自拍偷在线| 久久久精品欧美日韩精品| 欧美xxxx黑人xx丫x性爽| 美女高潮的动态| 日日摸夜夜添夜夜添小说| 久久精品91蜜桃| 搡女人真爽免费视频火全软件 | 国产精品乱码一区二三区的特点| 成人国产一区最新在线观看| 欧美黑人欧美精品刺激| 日韩高清综合在线| 欧美一级毛片孕妇| 久久久久免费精品人妻一区二区| 在线观看av片永久免费下载| 欧美黄色片欧美黄色片| 成人无遮挡网站| 亚洲一区高清亚洲精品| 日日夜夜操网爽| 在线观看66精品国产| 午夜福利18| 最近视频中文字幕2019在线8| 午夜激情欧美在线| 啦啦啦韩国在线观看视频| 午夜福利在线观看免费完整高清在 | 真人做人爱边吃奶动态| 人妻夜夜爽99麻豆av| 午夜免费观看网址| 免费看十八禁软件| 久久人妻av系列| 又紧又爽又黄一区二区| 婷婷精品国产亚洲av在线| 国产黄a三级三级三级人| 亚洲欧美日韩高清专用| 尤物成人国产欧美一区二区三区| 免费看光身美女| 国产一区在线观看成人免费| 国产精品精品国产色婷婷| 女警被强在线播放| 日日夜夜操网爽| 老汉色∧v一级毛片| 国产主播在线观看一区二区| av在线天堂中文字幕| 亚洲人与动物交配视频| 最近最新中文字幕大全电影3| 成人精品一区二区免费| av福利片在线观看| 九九久久精品国产亚洲av麻豆| 精品免费久久久久久久清纯| 国内精品久久久久久久电影| 国产精品电影一区二区三区| 老汉色∧v一级毛片| 97碰自拍视频| 免费av毛片视频| 亚洲精品久久国产高清桃花| 在线免费观看不下载黄p国产 | 性欧美人与动物交配| 一本综合久久免费| 亚洲真实伦在线观看| 极品教师在线免费播放| 色视频www国产| 国产亚洲精品av在线| 女人高潮潮喷娇喘18禁视频| 欧美性感艳星| 欧美中文日本在线观看视频| 内射极品少妇av片p| 男人舔女人下体高潮全视频| 亚洲七黄色美女视频| 99久久精品国产亚洲精品| 亚洲午夜理论影院| 一区二区三区高清视频在线| 哪里可以看免费的av片| 1024手机看黄色片| 日日夜夜操网爽| 精品人妻偷拍中文字幕| 国产一区二区激情短视频| 国产精品美女特级片免费视频播放器| 欧美最黄视频在线播放免费| 欧美区成人在线视频| 国产精品久久久人人做人人爽| 免费看十八禁软件| 国产精品 欧美亚洲| 麻豆成人av在线观看| 亚洲美女视频黄频| 天天躁日日操中文字幕| 精品99又大又爽又粗少妇毛片 | 黄色女人牲交| 久久久久国内视频| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 午夜精品在线福利| 怎么达到女性高潮| 成人av一区二区三区在线看| 国产一区二区在线av高清观看| bbb黄色大片| 久久久久免费精品人妻一区二区| 亚洲真实伦在线观看| 我要搜黄色片| 亚洲第一欧美日韩一区二区三区| 岛国在线观看网站| 婷婷六月久久综合丁香| 最后的刺客免费高清国语| xxx96com| av中文乱码字幕在线| 精品免费久久久久久久清纯| 国产av在哪里看| 亚洲狠狠婷婷综合久久图片| 久久久久性生活片| 两性午夜刺激爽爽歪歪视频在线观看| 最好的美女福利视频网| 亚洲欧美一区二区三区黑人| 欧美另类亚洲清纯唯美| 国产视频一区二区在线看| 亚洲熟妇中文字幕五十中出| 亚洲国产色片| 国产亚洲精品综合一区在线观看| 久久人人精品亚洲av| 亚洲七黄色美女视频| 在线播放国产精品三级| 久久精品91蜜桃| 亚洲av成人av| 久久久久性生活片| 岛国在线观看网站| 黄片小视频在线播放| 一进一出抽搐动态| 欧美+日韩+精品| 国内揄拍国产精品人妻在线| 性色avwww在线观看| 天堂网av新在线| 亚洲av成人不卡在线观看播放网| 午夜亚洲福利在线播放| 成年女人毛片免费观看观看9| 成人国产综合亚洲| 99久久精品热视频| 亚洲成av人片免费观看| 精品久久久久久久久久免费视频| ponron亚洲| 国产伦人伦偷精品视频| 99热这里只有精品一区| 国产亚洲精品久久久久久毛片| 老汉色av国产亚洲站长工具| 亚洲人成伊人成综合网2020| 国产精品亚洲美女久久久| 少妇丰满av| 日韩高清综合在线| 亚洲精品乱码久久久v下载方式 | 国产高清有码在线观看视频| 99热精品在线国产| 变态另类丝袜制服| 日韩欧美在线乱码| 午夜福利高清视频| 日韩av在线大香蕉| 人人妻人人看人人澡| 亚洲激情在线av| 成人性生交大片免费视频hd| 伊人久久精品亚洲午夜| 99久久精品热视频| tocl精华| 无限看片的www在线观看| 日本 av在线| 一级作爱视频免费观看| 搞女人的毛片| 美女高潮喷水抽搐中文字幕| 午夜激情福利司机影院| 少妇高潮的动态图| 99在线人妻在线中文字幕| 首页视频小说图片口味搜索| 小说图片视频综合网站| 午夜福利免费观看在线| 成人亚洲精品av一区二区| 悠悠久久av| 亚洲av第一区精品v没综合| 黄色丝袜av网址大全| 亚洲av成人不卡在线观看播放网| 九色国产91popny在线| 一级毛片女人18水好多| av天堂在线播放| 日韩欧美精品免费久久 | bbb黄色大片| 99热精品在线国产| 精华霜和精华液先用哪个| 国产精品嫩草影院av在线观看 | 国产精品乱码一区二三区的特点| 久久亚洲精品不卡| 99精品久久久久人妻精品| 日韩国内少妇激情av| 亚洲,欧美精品.| 国产成人欧美在线观看| 亚洲国产高清在线一区二区三| 久久久精品欧美日韩精品| 日本在线视频免费播放| 日韩免费av在线播放| 母亲3免费完整高清在线观看| 国产av麻豆久久久久久久| 午夜福利欧美成人| 最近在线观看免费完整版| 日本熟妇午夜| 久久中文看片网| 亚洲欧美激情综合另类| 欧美区成人在线视频| 欧美日韩乱码在线| 黑人欧美特级aaaaaa片| 色在线成人网| 免费看十八禁软件| 看黄色毛片网站| 国产欧美日韩一区二区三| 无人区码免费观看不卡| 欧美三级亚洲精品| 国产黄色小视频在线观看| 怎么达到女性高潮| 久久精品国产亚洲av香蕉五月| 久久久久性生活片| 国产老妇女一区| 少妇裸体淫交视频免费看高清| 国产成人a区在线观看| 日韩 欧美 亚洲 中文字幕| 别揉我奶头~嗯~啊~动态视频| 人妻夜夜爽99麻豆av| 性色av乱码一区二区三区2| 精品人妻一区二区三区麻豆 | 国产激情欧美一区二区| 草草在线视频免费看| 桃色一区二区三区在线观看| 中文亚洲av片在线观看爽| 午夜亚洲福利在线播放| 99久久精品一区二区三区| 国产精品一及| 蜜桃亚洲精品一区二区三区| 91字幕亚洲| 国产精华一区二区三区| 免费无遮挡裸体视频| 我要搜黄色片| 欧美最黄视频在线播放免费| 国产高清激情床上av| 日本黄色片子视频| 成人国产一区最新在线观看| 久久九九热精品免费| 黄色日韩在线| 91在线精品国自产拍蜜月 | 午夜a级毛片| 欧美丝袜亚洲另类 | 亚洲欧美日韩卡通动漫| 韩国av一区二区三区四区| 久久精品影院6| 高潮久久久久久久久久久不卡| 男女做爰动态图高潮gif福利片| 又爽又黄无遮挡网站| 亚洲av成人不卡在线观看播放网| 十八禁人妻一区二区| 老熟妇乱子伦视频在线观看| 偷拍熟女少妇极品色| a级一级毛片免费在线观看| 美女黄网站色视频| 高清在线国产一区| eeuss影院久久| 久久精品国产99精品国产亚洲性色| 中文字幕av成人在线电影| 一本一本综合久久| 九九在线视频观看精品| 天天躁日日操中文字幕| 日本熟妇午夜| 性色av乱码一区二区三区2| 中文资源天堂在线| 露出奶头的视频| 中文字幕av在线有码专区| 九色成人免费人妻av| 国产成人a区在线观看| 色吧在线观看| 精品人妻一区二区三区麻豆 | 老司机深夜福利视频在线观看| 淫秽高清视频在线观看| 国产精品亚洲美女久久久| 国产三级黄色录像| 波多野结衣巨乳人妻| 亚洲精品美女久久久久99蜜臀| 亚洲天堂国产精品一区在线| 欧美午夜高清在线| www.熟女人妻精品国产| 草草在线视频免费看| 精品99又大又爽又粗少妇毛片 | 美女高潮的动态| 最好的美女福利视频网| 三级毛片av免费| 国产三级中文精品| 熟女人妻精品中文字幕| 91久久精品电影网| 一本一本综合久久| 不卡一级毛片| 国产亚洲精品久久久com| 欧美日韩综合久久久久久 | 成熟少妇高潮喷水视频| 最新中文字幕久久久久| 首页视频小说图片口味搜索| 一区二区三区高清视频在线| 别揉我奶头~嗯~啊~动态视频| ponron亚洲| 免费在线观看成人毛片| 亚洲真实伦在线观看| 色视频www国产| 日韩成人在线观看一区二区三区| 欧美黄色淫秽网站| 一级a爱片免费观看的视频| 欧美极品一区二区三区四区| av欧美777| 99久久成人亚洲精品观看| 亚洲,欧美精品.| 中文字幕久久专区| 琪琪午夜伦伦电影理论片6080| 久久久久久久久中文| av在线天堂中文字幕| 特级一级黄色大片| 欧美zozozo另类| 国产老妇女一区| 成年人黄色毛片网站| 国产高潮美女av| 国产成人av教育| 日韩成人在线观看一区二区三区|