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

    拋錨貫入深度計(jì)算方法比較研究

    2020-06-12 04:38:34馮士倫朱曉宇唐友剛趙志娟
    海洋技術(shù)學(xué)報(bào) 2020年2期
    關(guān)鍵詞:觸底砂土土質(zhì)

    馮士倫 ,朱曉宇 ,李 焱*,唐友剛 ,趙志娟

    (1.天津大學(xué) 水利工程仿真與安全國(guó)家重點(diǎn)實(shí)驗(yàn)室,天津 300350;2.天津大學(xué) 天津市港口與海洋工程重點(diǎn)實(shí)驗(yàn)室,天津 300350;3.天津大學(xué) 建筑工程學(xué)院,天津 300350)

    隨著世界經(jīng)濟(jì)的發(fā)展,我國(guó)的航運(yùn)業(yè)得到了充分的發(fā)展,商船進(jìn)行拋錨作業(yè)的頻次大幅增加,同時(shí),海底輸油管道等海底設(shè)施也不斷增多,商船拋錨引發(fā)海底管線事故時(shí)有發(fā)生。錨貫入海底過(guò)程中恰好擊中海底油氣管道等海底設(shè)施,可能導(dǎo)致海底管道發(fā)生泄露或設(shè)施損傷,從而對(duì)原油輸送和油田作業(yè)帶來(lái)巨大損失。為避免船舶拋錨造成海底管道等海底設(shè)施的破壞,通常的做法是將海底管線或海底設(shè)施埋入土中一定深度。因此,船舶拋錨后貫入海底土深度計(jì)算,便成為了決定海底管線等海底設(shè)施埋設(shè)深度的關(guān)鍵因素。

    陳峰等[1]基于能量法計(jì)算了黏土與砂土土質(zhì)下大抓力錨的貫入深度,并進(jìn)行了現(xiàn)場(chǎng)試驗(yàn),證實(shí)了能量法可以較好地?cái)M合試驗(yàn)的實(shí)測(cè)落深,揭示了船錨的質(zhì)量和下落高度對(duì)錨的貫入深度的影響規(guī)律。莊元等[2]提出了一種較為簡(jiǎn)單的物理模型,分析并預(yù)測(cè)了船錨在海底土中的貫入深度。李曉松等[3]分析了船舶拋錨時(shí)錨的運(yùn)動(dòng)全程的受力,建立了船錨觸底速度的求解模型,而后使用微分方程法計(jì)算了錨最終的貫入深度。雷震名等[4]將能量法與正交試驗(yàn)法結(jié)合,研究了海底管道對(duì)不同的影響因素的敏感性,得出了管道的最小推薦設(shè)計(jì)埋深。李慶等[5]采用DNV的相關(guān)方法,計(jì)算船錨在下落過(guò)程中的有效沖擊能量和落錨的最終貫穿深度,并用微分方程法計(jì)算了拋錨-拖錨-停滯的最終深度,為穿越航道的海底管道的保護(hù)提供了一定的理論指導(dǎo)。

    Y H Kim等[6]對(duì)魚雷錨的動(dòng)態(tài)貫入深度進(jìn)行了深入研究,研究表明:錨的貫入深度是沖擊速度、剪切強(qiáng)度對(duì)應(yīng)變率的依賴性、剪切過(guò)程中的重塑程度以及土壤敏感性的函數(shù)。Wenkai Wang等[7]在實(shí)驗(yàn)室用9種不同的魚雷錨進(jìn)行了128組試驗(yàn),根據(jù)能量守恒原理和試驗(yàn)測(cè)量結(jié)果,提出了不考慮土壤分離工況下錨桿貫入深度的計(jì)算公式。結(jié)果表明,能量法計(jì)算的侵徹深度與實(shí)驗(yàn)室和現(xiàn)場(chǎng)實(shí)測(cè)的侵徹深度吻合較好。Haixiao Liu等[8]通過(guò)研究粘土和砂土的參數(shù),研究了承載力系數(shù)、粘滯系數(shù)、不排水抗剪強(qiáng)度、內(nèi)摩擦角、側(cè)土壓力系數(shù)等因素的影響,推導(dǎo)了在粘土與砂土中計(jì)算錨極限埋深的公式,結(jié)果表明,在給定合適的參數(shù)值的情況下,可以合理地預(yù)測(cè)粘性土中錨桿的極限承載力。其中,參數(shù)的選擇是影響計(jì)算精度的最重要的因素。

    船舶拋錨貫入海底的深度的影響因素十分復(fù)雜,已知的主要有錨型、錨重、觸底速度、海底土物理特性等。目前,世界上各個(gè)國(guó)家對(duì)錨貫入深度的計(jì)算方法尚未形成統(tǒng)一的認(rèn)識(shí),還沒(méi)有明確提出貫入深度的通用計(jì)算公式。現(xiàn)有的船錨貫入深度研究方法主要有基于DNV規(guī)范的能量法與基于牛頓第二定律的微分方程法。錨的觸底速度大小直接影響錨貫入海底時(shí)的沖擊能量,它也是公式法計(jì)算時(shí)的重要影響因素。因此,無(wú)論哪種計(jì)算方法,觸底速度的精確計(jì)算都十分有必要。本文以不同重量的霍爾錨為例,基于CFD方法計(jì)算拋錨下落過(guò)程的速度及觸底速度,使用正交試驗(yàn)原理分析船錨貫入深度對(duì)不同影響因素的敏感性,為海底管線等海底設(shè)施的埋設(shè)深度提供依據(jù)。

    1 拋錨貫入深度計(jì)算方法

    目前計(jì)算拋錨貫入深度使用較多的是DNV規(guī)范中的能量法與利用牛頓第二定律推導(dǎo)出的微分方程法。

    1.1 能量法

    DNV-RP-F107給出了采用能量法計(jì)算船錨的落深的方法[9],認(rèn)為船錨的沖擊能量完全被土體吸收,基于能量守恒原理,土體對(duì)船錨所做的功等于船錨在下落過(guò)程中的沖擊總能量,由此可以計(jì)算出船錨的貫入深度。

    錨對(duì)海底土的沖擊能量,與結(jié)構(gòu)自身質(zhì)量和附連水質(zhì)量有關(guān),包括附連水在內(nèi)錨體的沖擊總能量為:

    式中:EE為錨觸底時(shí)的沖擊總能量,kJ;Ep為錨觸底時(shí)的動(dòng)能,kJ;EA為附加水動(dòng)力能量,kJ;m為錨的質(zhì)量,kg;ρw為海水的密度,kg/m3;V為附加水的體積,m3;Ca為附加質(zhì)量系數(shù),該值取決于物體的幾何形狀,根據(jù)工程經(jīng)驗(yàn),本次研究中對(duì)于掉落的錨取Ca=1.5。

    海底土質(zhì)吸收船錨的沖擊能量,可以用式(2)計(jì)算[6]:

    式中:Ek為海底土質(zhì)吸收的能量,kJ;h為船錨貫入深度,m;γ'為土的有效單位重量,kN/m3;Sr為錨的形狀系數(shù),與錨冠水平投影面積有關(guān),依據(jù)規(guī)范要求,霍爾錨的形狀系數(shù)取值為0.6,其余錨的形狀系數(shù)可同樣參考該規(guī)范取值;Nr為與內(nèi)摩擦角有關(guān)的土承載力系數(shù)。根據(jù)API RP 2GEO規(guī)范中的推薦方法,如下式所示:

    式中:φ 為土的內(nèi)摩擦角,(°);Nq與Nγ為 φ 的無(wú)量綱函數(shù),Nq=(exp[πtanφ'])(tan2(45°+φ/2));Nγ=1.5(Nq-1)tanφ

    1.2 微分方程法

    錨在貫入海底土質(zhì)的運(yùn)動(dòng)過(guò)程中會(huì)受到重力、土端阻力、土摩擦阻力、拖曳力等的作用,從而使錨做減速運(yùn)動(dòng)并最終停止。貫入深度的影響因素很多,以下方程式主要反映了錨的觸底速度、錨型、錨重、海底土特性等的影響。

    錨貫入海底土的運(yùn)動(dòng)可以使用牛頓第二定律表示:

    式中:M*為錨在土中的有效質(zhì)量,kg;W為錨的水下重,等于錨的重力減去浮力,N;y為錨在土中的貫入深度,m;FBP為土底部承載力,N;FAD為土側(cè)向承載力,N;FH為拖曳力,N。具體計(jì)算見(jiàn)式(5)~式(10):

    式中:Su為土的不排水抗剪強(qiáng)度,N/m2;Nc為土極限承載力系數(shù);Nγ,Nq為土承載力系數(shù),與上節(jié)中相同;δ為土的側(cè)向承載力系數(shù);St為黏土的靈敏度;AF為錨的前端面面積,m2;As為錨的側(cè)面面積,m2;CD為拖曳力系數(shù);ν為錨的速度,m/s;P0為深度為y的有效覆壓,N/m2;γ為土在水下的單位重量,N/m3;C為錨面短邊的長(zhǎng)度,m;α為土和錨之間的摩擦角,(°);q0為平均有效覆壓,N/m2;Kp為被動(dòng)土壓力系數(shù);φ為土的內(nèi)摩擦角,(°)。

    本文開發(fā)了計(jì)算機(jī)程序,基于微分法計(jì)算錨體下落過(guò)程的軌跡和貫入深度。當(dāng)錨的速度減小到0時(shí),對(duì)應(yīng)錨貫入海底的最終深度。

    2 參數(shù)選擇

    2.1 海底土質(zhì)參數(shù)

    當(dāng)船錨接觸到海底土以后,會(huì)對(duì)海底土產(chǎn)生一個(gè)豎直向下的作用力,由于相互作用力的影響,海底土很容易出現(xiàn)變形,并給船錨一個(gè)向上的反作用力。在該力的作用下,船錨會(huì)不斷減速,直到速度減為0。貫入深度與海底土質(zhì)密切相關(guān),海底土一般有淤泥、石子、黏土、砂土等,因此需對(duì)分析海底土質(zhì)的參數(shù)對(duì)于貫入深度的影響。本文以3種典型黏土與1種典型砂土為例,分別計(jì)算船錨的貫入深度,其主要物理參數(shù)如表1所示。

    表1 海底土的主要物理性質(zhì)參數(shù)

    根據(jù)文獻(xiàn)[10],黏土的側(cè)向承載力系數(shù)取值為0.23~0.40,本文取為 0.40。粘土的靈敏度[11]St:一般介于1.5~5,計(jì)算中取St=3;無(wú)量綱Nc為土的極限承載系數(shù)[12],在錨貫入計(jì)算中取為9.0。

    2.2 錨參數(shù)

    使用以上兩種方法計(jì)算船錨貫入量的過(guò)程中,除了錨的觸底速度ν,錨的質(zhì)量M、錨的水平投影面積Af和側(cè)投影面積As同樣是非常重要的參數(shù)。

    霍爾錨是商船拋錨常用的錨型,小型、中型、大型商船拋錨需要分別使用相對(duì)應(yīng)重量的霍爾錨。以下計(jì)算分別選取重量為2 t,5 t,8 t和10 t的霍爾錨,同時(shí)使用其對(duì)應(yīng)的投影面積。

    當(dāng)計(jì)算錨的投影面積時(shí),參照富志禹等[13]提出的船用錨水平投影面積估算公式,計(jì)算得到霍爾錨在不同錨重下錨冠水平投影面積,使用2016年發(fā)布的國(guó)標(biāo)《霍爾錨》(GB/T 546-2016)中霍爾錨的數(shù)據(jù),計(jì)算錨的側(cè)投影面積。

    表2 霍爾錨錨冠投影面積

    根據(jù)DNV-RP F107規(guī)范規(guī)定,物體從球形到更復(fù)雜形式的拖曳阻力系數(shù)范圍為0.6~2.0。根據(jù)文獻(xiàn)[14],霍爾錨的拖曳力系數(shù)Cd取1.0時(shí),其極限速度偏小,受到的阻力偏大。本文保守假定霍爾錨的拖曳力系數(shù)Cd=0.7。

    3 貫入深度計(jì)算

    3.1 拋錨觸底速度分析

    船舶拋錨時(shí),將錨從錨鏈孔處釋放,直至錨接觸海底土的過(guò)程中,錨分別經(jīng)歷了錨穴-水面、水面到泥線兩個(gè)不同的運(yùn)動(dòng)過(guò)程,觸底速度與錨在空氣中運(yùn)動(dòng)和水中運(yùn)動(dòng)過(guò)程有關(guān)。忽略錨鏈的作用,錨在空氣中運(yùn)動(dòng)時(shí)僅受重力作用;錨在水中運(yùn)動(dòng)時(shí)受重力、水阻力的作用。

    采用ICEM軟件對(duì)霍爾錨進(jìn)行建模和流場(chǎng)、網(wǎng)格劃分,編寫錨自由下落的UDF文件,導(dǎo)入FLUENT軟件中,利用CFD原理對(duì)錨下落20 m以及25 m后的下落速度進(jìn)行計(jì)算。使用ICEM軟件建立的簡(jiǎn)化模型如圖1所示。計(jì)算取15 m×15 m×40 m的矩形流體域,并進(jìn)行動(dòng)網(wǎng)格劃分,如圖2所示。

    流域頂部設(shè)置為壓力出口,四周設(shè)置為對(duì)稱邊界,以模擬無(wú)限范圍的流域。使用VOF方法模擬液體和空氣的兩相流動(dòng),采用動(dòng)網(wǎng)格中的自動(dòng)更新網(wǎng)格方法控制結(jié)構(gòu)下落時(shí)的網(wǎng)格質(zhì)量,使用PISO方法進(jìn)行數(shù)值離散,初始時(shí)刻,設(shè)置拋錨的位置處于水面以上5 m。計(jì)算時(shí)間步長(zhǎng)取0.002 s,總計(jì)算時(shí)間為8 s。

    圖1 霍爾錨的簡(jiǎn)化模型

    圖2 計(jì)算域及網(wǎng)格劃分

    根據(jù)上一節(jié)的計(jì)算設(shè)置,得到結(jié)構(gòu)下落8 s后的速度歷程和距離歷程,從而得到結(jié)構(gòu)在水中運(yùn)動(dòng)20 m以及25 m處的下落速度,縱坐標(biāo)0點(diǎn)為水面,水面以上設(shè)為正值,如圖3所示。

    圖3 霍爾錨下落的水深與速度曲線

    從圖3中可以得到如下結(jié)論:(1)霍爾錨在水中運(yùn)動(dòng)20 m和25 m后,觸底速度主要集中在5~7 m/s之間;(2)觸底速度隨水深的增大而增大;(3)結(jié)合文獻(xiàn)[15]中公式4.9,霍爾錨投影面積增加的速度小于質(zhì)量增加的速度,因此,觸底速度隨錨質(zhì)量的增大而增大。

    3.2 錨貫入深度分析

    3.2.1 能量法 分別選取錨重為 2 t,5 t,8 t,10 t的霍爾錨,根據(jù)計(jì)算得到的觸底速度,計(jì)算霍爾錨下落過(guò)程中具有的沖擊總能量,如表3所示。

    表3 不同錨重拋錨觸底速度與沖擊總能量

    針對(duì)4種不同土質(zhì),使用能量法的式(1)~式(2),開展船錨貫入深度計(jì)算,分別得到不同錨重的霍爾錨在水中下落20 m與25 m的相關(guān)計(jì)算結(jié)果,如圖4~圖5所示。

    由圖4~圖5可見(jiàn),不同土質(zhì)對(duì)應(yīng)的霍爾錨貫入深度由大到小依次為淤泥、軟黏土、一般黏土、砂土,這表明霍爾錨貫入海底土的深度與土質(zhì)有關(guān),隨土質(zhì)的變化而發(fā)生變化。同時(shí),不同錨重下的貫入深度不同,10 t重的霍爾錨在水中下落25 m后貫入淤泥的深度最大,為3.89 m;2 t重的霍爾錨在水中下落20 m后貫入砂土的深度最小,為0.86 m,這表明霍爾錨貫入海底土的深度與錨重有關(guān),隨錨重的變化而發(fā)生變化。

    3.2.2 微分方程法 針對(duì)黏土和砂土這兩種海底土質(zhì),基于式(4)~式(10)使用微分方程法,開展船錨貫入深度計(jì)算,分別得到不同錨重下霍爾錨在水中下落20 m與25 m的相關(guān)計(jì)算結(jié)果,如圖6所示。

    由圖6可知,霍爾錨貫入海底土的深度受錨重、海底土質(zhì)、水深等的共同影響?;魻栧^貫入海底土的深度隨錨重的增加而增大;霍爾錨貫入黏土層的深度明顯比砂土層更大;水深增大,霍爾錨的貫入深度也隨之增大,但是對(duì)黏土層的影響比砂土層更大。

    圖4 霍爾錨下落20 m的貫入深度

    圖5 霍爾錨下落25 m的貫入深度

    圖6 霍爾錨不同錨重的貫入深度

    4 正交試驗(yàn)分析

    拋錨貫入深度明顯受到多個(gè)因素的控制和影響。通過(guò)正交試驗(yàn)分析,對(duì)兩種計(jì)算方法進(jìn)行參數(shù)敏感性分析,明確兩種方法影響因素的主次順序,這可為實(shí)際工程計(jì)算提供指導(dǎo)。

    4.1 能量法

    使用能量法計(jì)算拋錨貫入深度過(guò)程中,錨重、觸底速度、土的有效重度、內(nèi)摩擦角是其主要影響因素。根據(jù)文獻(xiàn)[16],黏土的有效重度最大約為17 kN/m3,砂土的內(nèi)摩擦角最大為35°,因此下面取合理有效重度范圍、合理內(nèi)摩擦角范圍的黏土為例,將這些影響因素作為正交試驗(yàn)的試驗(yàn)因子,計(jì)算出實(shí)驗(yàn)所需參數(shù),將其用于下面的正交表中,如表4所示。

    表4 正交試驗(yàn)表

    由表4可知,最大極差值為2.61,最小極差值為0.87,分別對(duì)應(yīng)內(nèi)摩擦角和錨重,這表明對(duì)于貫入深度影響最大的因素為內(nèi)摩擦角。從表4中標(biāo)準(zhǔn)差也可以得到相同的結(jié)論。

    4.2 微分方程法

    對(duì)微分方程法進(jìn)行綜合分析研究,可知,使用微分方程法計(jì)算貫入深度過(guò)程中,需要將海底底質(zhì)分為黏土和砂土兩種不同底質(zhì)狀況。底質(zhì)為黏土?xí)r,錨重、觸底速度、土的有效重度、不排水抗剪強(qiáng)度是影響貫入深度的主要因素;底質(zhì)為砂土?xí)r,錨重、觸底速度、土的有效重度、內(nèi)摩擦角是其主要影響因素。下面以海底底質(zhì)為黏土?xí)r為例,計(jì)算出實(shí)驗(yàn)所需參數(shù),將其用于下面的正交表中,如表5所示。

    表5 正交試驗(yàn)表

    由表5可知,最大極差值為1.39,最小極差值為1.16,分別對(duì)應(yīng)有效重度和錨重,這表明對(duì)于貫入深度影響最大的因素為有效重度。從表5中標(biāo)準(zhǔn)差也可以得到相同的結(jié)論。

    5 結(jié)論

    本文采用能量法與微分方程法計(jì)算了船舶拋錨貫入海底土的深度,得到如下結(jié)論:

    (1)文中采用能量法計(jì)算出的貫入深度值總體來(lái)看比微分方程法要大,這是因?yàn)槭褂媚芰糠ㄓ?jì)算沒(méi)有考慮能量耗散,所以計(jì)算結(jié)果較大。

    (2)使用能量法與微分方程法計(jì)算出的結(jié)果曲線變化趨勢(shì)一致,表明兩種計(jì)算結(jié)果都具有合理性,能量法計(jì)算的貫入深度較大,以此為基準(zhǔn)鋪設(shè)海底管線較為安全。

    (3)通過(guò)正交試驗(yàn)可知,兩種計(jì)算方法均遵循以下規(guī)律:隨著土的有效重度增大,錨受到的阻力增大,其貫入海底的深度隨之減??;土的內(nèi)摩擦角增大,土體之間的摩擦力增大,土更加難以貫入其中,因此錨貫入海底的深度隨之減小。

    (4)船舶拋錨貫入海底的深度大小受到錨的重量、錨的觸底速度、土的有效重度、土的內(nèi)摩擦角、土的不排水抗剪強(qiáng)度等因素的影響,使用能量法計(jì)算過(guò)程中,參數(shù)敏感性依次為:內(nèi)摩擦角>有效重度>觸底速度>錨重;使用微分方程法計(jì)算過(guò)程中,參數(shù)敏感性依次為:土的有效重度>觸底速度>不排水抗剪強(qiáng)度>錨重。

    (5)使用微分方程法計(jì)算過(guò)程中,最大極差值1.39與最小極差值1.16相差0.23,差距不大,說(shuō)明敏感性差距不大。

    因此,對(duì)于使用能量法計(jì)算貫入深度時(shí),內(nèi)摩擦角的選取盡量合理準(zhǔn)確;而在使用微分方程法計(jì)算貫入深度時(shí),4個(gè)主要影響參數(shù)均應(yīng)盡量合理選取,全面考慮。

    猜你喜歡
    觸底砂土土質(zhì)
    飽和砂土地層輸水管道施工降水方案設(shè)計(jì)
    高含鐵大比重土質(zhì)對(duì)泥漿配比的影響
    龍之中華 龍之砂土——《蟠龍壺》創(chuàng)作談
    平安千億回購(gòu) 股指觸底回升
    凍融循環(huán)作用下土質(zhì)河堤的穩(wěn)定性分析
    土質(zhì)文物鹽害中硫酸鈉的研究——從微觀到宏觀
    如何判斷經(jīng)濟(jì)觸底?
    城市淺埋隧道穿越飽和砂土復(fù)合地層時(shí)適宜的施工工法
    市場(chǎng)觸底了嗎?
    溯源侵蝕作用下昔格達(dá)地層土質(zhì)邊坡的穩(wěn)定性
    欧美日韩综合久久久久久| 免费看美女性在线毛片视频| 亚洲欧美成人综合另类久久久| 日韩人妻高清精品专区| 肉色欧美久久久久久久蜜桃 | 免费在线观看成人毛片| 日本欧美国产在线视频| 69人妻影院| 国产黄色小视频在线观看| 久久久精品欧美日韩精品| 偷拍熟女少妇极品色| 综合色丁香网| 精品久久久精品久久久| 国产在视频线精品| 嘟嘟电影网在线观看| 又爽又黄a免费视频| 免费看不卡的av| 欧美最新免费一区二区三区| 菩萨蛮人人尽说江南好唐韦庄| 噜噜噜噜噜久久久久久91| 亚洲四区av| 亚洲不卡免费看| 特级一级黄色大片| 亚洲精品久久久久久婷婷小说| 五月玫瑰六月丁香| 插阴视频在线观看视频| 午夜久久久久精精品| 免费大片黄手机在线观看| 国产极品天堂在线| 丰满乱子伦码专区| 午夜激情福利司机影院| 亚洲av中文字字幕乱码综合| 国产免费一级a男人的天堂| 99久久精品热视频| 国产黄色小视频在线观看| 免费在线观看成人毛片| 久久久欧美国产精品| 亚洲电影在线观看av| 日韩伦理黄色片| 午夜亚洲福利在线播放| 国产精品一二三区在线看| 亚洲丝袜综合中文字幕| 91久久精品电影网| 有码 亚洲区| 97热精品久久久久久| 欧美日韩亚洲高清精品| 国产一区二区三区综合在线观看 | 日韩成人伦理影院| 久久精品人妻少妇| 久久久精品94久久精品| 国产人妻一区二区三区在| 精品人妻一区二区三区麻豆| 九九在线视频观看精品| 国产精品国产三级国产av玫瑰| 两个人视频免费观看高清| 国产黄色小视频在线观看| 91在线精品国自产拍蜜月| 国内精品宾馆在线| 五月天丁香电影| 亚洲国产色片| 久久精品国产亚洲av涩爱| 亚洲色图av天堂| 国产老妇女一区| 大话2 男鬼变身卡| 人体艺术视频欧美日本| 你懂的网址亚洲精品在线观看| 人妻夜夜爽99麻豆av| 婷婷色麻豆天堂久久| 亚洲熟女精品中文字幕| 国产不卡一卡二| 少妇丰满av| 亚洲av.av天堂| 男女国产视频网站| 人妻系列 视频| 精品国内亚洲2022精品成人| 免费观看在线日韩| 床上黄色一级片| 一本一本综合久久| 秋霞伦理黄片| 亚洲精品自拍成人| 搡老乐熟女国产| 午夜激情福利司机影院| 国产精品99久久久久久久久| 黄片wwwwww| 一级毛片久久久久久久久女| 久热久热在线精品观看| 少妇熟女aⅴ在线视频| av在线观看视频网站免费| 男女边吃奶边做爰视频| 国产精品国产三级国产专区5o| 午夜免费男女啪啪视频观看| 少妇人妻一区二区三区视频| 亚洲欧美日韩卡通动漫| 中文字幕制服av| 久久久久免费精品人妻一区二区| 少妇人妻一区二区三区视频| 91久久精品国产一区二区成人| 99热这里只有是精品在线观看| 久久久欧美国产精品| 中文字幕av成人在线电影| 亚洲av成人av| 天堂俺去俺来也www色官网 | 亚洲自偷自拍三级| 亚洲欧美一区二区三区国产| 日韩成人伦理影院| 久久久久精品久久久久真实原创| 久久午夜福利片| 精品亚洲乱码少妇综合久久| av在线蜜桃| 亚洲内射少妇av| 69av精品久久久久久| 69人妻影院| 国产国拍精品亚洲av在线观看| 国产精品人妻久久久影院| 男女边吃奶边做爰视频| 少妇裸体淫交视频免费看高清| 性插视频无遮挡在线免费观看| 男人和女人高潮做爰伦理| 亚洲成人一二三区av| 成人美女网站在线观看视频| 亚洲,欧美,日韩| 成人毛片60女人毛片免费| 热99在线观看视频| 国产av不卡久久| 一区二区三区乱码不卡18| 韩国高清视频一区二区三区| 午夜精品一区二区三区免费看| 亚洲成色77777| 国产极品天堂在线| 又粗又硬又长又爽又黄的视频| 国产免费视频播放在线视频 | 国产精品国产三级专区第一集| 久久久国产一区二区| 丰满少妇做爰视频| 精品人妻一区二区三区麻豆| av在线亚洲专区| 亚洲av国产av综合av卡| 在线免费观看的www视频| 最近视频中文字幕2019在线8| 18禁裸乳无遮挡免费网站照片| 欧美成人精品欧美一级黄| 欧美97在线视频| 国产男女超爽视频在线观看| 人人妻人人澡欧美一区二区| 亚洲18禁久久av| 性插视频无遮挡在线免费观看| 久久99热6这里只有精品| 亚洲丝袜综合中文字幕| 在线天堂最新版资源| 美女主播在线视频| 亚洲精品国产成人久久av| 国产精品嫩草影院av在线观看| 亚洲av不卡在线观看| 国产精品不卡视频一区二区| 永久免费av网站大全| 一级毛片 在线播放| 性插视频无遮挡在线免费观看| 亚洲av不卡在线观看| 高清视频免费观看一区二区 | 欧美日韩国产mv在线观看视频 | 老师上课跳d突然被开到最大视频| 成人高潮视频无遮挡免费网站| www.色视频.com| 嫩草影院新地址| 麻豆乱淫一区二区| 少妇熟女aⅴ在线视频| 日韩人妻高清精品专区| 成人特级av手机在线观看| 麻豆精品久久久久久蜜桃| 欧美成人精品欧美一级黄| 国产有黄有色有爽视频| 亚洲国产av新网站| 欧美成人一区二区免费高清观看| 国产麻豆成人av免费视频| 国产 亚洲一区二区三区 | 久久精品国产鲁丝片午夜精品| 秋霞在线观看毛片| 少妇人妻精品综合一区二区| 国产精品美女特级片免费视频播放器| 国产在视频线在精品| 一级毛片aaaaaa免费看小| 18+在线观看网站| 久久久久免费精品人妻一区二区| 国产大屁股一区二区在线视频| 十八禁国产超污无遮挡网站| 国产免费又黄又爽又色| 中文在线观看免费www的网站| 成人高潮视频无遮挡免费网站| 日韩视频在线欧美| 日韩电影二区| 亚洲在线观看片| 亚洲人成网站高清观看| av天堂中文字幕网| 国产免费一级a男人的天堂| 亚洲精品视频女| 黄片无遮挡物在线观看| 国产亚洲午夜精品一区二区久久 | 永久网站在线| 精品一区二区三区人妻视频| 中文字幕av成人在线电影| 亚洲国产av新网站| 日日撸夜夜添| 精品久久国产蜜桃| 丰满乱子伦码专区| 国精品久久久久久国模美| 国产成人a∨麻豆精品| 天堂俺去俺来也www色官网 | 直男gayav资源| 亚洲国产成人一精品久久久| 国产乱人偷精品视频| 午夜久久久久精精品| 国产精品人妻久久久久久| 全区人妻精品视频| av网站免费在线观看视频 | 一夜夜www| 又粗又硬又长又爽又黄的视频| 亚洲综合色惰| 少妇人妻一区二区三区视频| 中文天堂在线官网| 边亲边吃奶的免费视频| 国产免费又黄又爽又色| 亚洲欧美日韩东京热| 亚洲国产高清在线一区二区三| 亚洲精品乱久久久久久| 亚洲美女视频黄频| 麻豆av噜噜一区二区三区| 欧美高清性xxxxhd video| 人妻制服诱惑在线中文字幕| 亚洲成人精品中文字幕电影| 日韩强制内射视频| 成人av在线播放网站| 午夜精品一区二区三区免费看| 熟女电影av网| 黑人高潮一二区| 国产亚洲av嫩草精品影院| av黄色大香蕉| 免费看光身美女| 九九久久精品国产亚洲av麻豆| 欧美成人一区二区免费高清观看| 亚洲三级黄色毛片| 日韩一本色道免费dvd| 七月丁香在线播放| 男人爽女人下面视频在线观看| 欧美变态另类bdsm刘玥| 午夜福利在线在线| 18禁在线播放成人免费| 国产91av在线免费观看| 国产精品1区2区在线观看.| av免费在线看不卡| 国产在视频线精品| 一级毛片aaaaaa免费看小| 神马国产精品三级电影在线观看| 一夜夜www| 亚洲人成网站高清观看| av在线观看视频网站免费| 国产精品久久视频播放| 成人一区二区视频在线观看| 亚洲精品亚洲一区二区| 日日干狠狠操夜夜爽| 高清在线视频一区二区三区| 肉色欧美久久久久久久蜜桃 | 老师上课跳d突然被开到最大视频| 秋霞伦理黄片| 日韩欧美精品v在线| 久久99热这里只频精品6学生| 成年av动漫网址| 伦理电影大哥的女人| 日韩av在线免费看完整版不卡| 亚洲国产欧美在线一区| 高清欧美精品videossex| 久久久久久久久中文| 高清午夜精品一区二区三区| 久久久久久久久大av| 最近的中文字幕免费完整| 成人亚洲精品av一区二区| 国内少妇人妻偷人精品xxx网站| 亚洲欧美精品自产自拍| 色吧在线观看| 嫩草影院入口| 99久国产av精品国产电影| 久久精品久久久久久久性| 看黄色毛片网站| eeuss影院久久| 日日干狠狠操夜夜爽| 欧美xxxx性猛交bbbb| 80岁老熟妇乱子伦牲交| 日韩一本色道免费dvd| 久久99热这里只频精品6学生| 三级毛片av免费| 成年版毛片免费区| 亚洲伊人久久精品综合| 最近2019中文字幕mv第一页| av专区在线播放| 日本-黄色视频高清免费观看| 狂野欧美激情性xxxx在线观看| 免费黄色在线免费观看| 99久久精品一区二区三区| 国产精品久久久久久久电影| 免费不卡的大黄色大毛片视频在线观看 | 纵有疾风起免费观看全集完整版 | 国产精品嫩草影院av在线观看| 乱码一卡2卡4卡精品| 大陆偷拍与自拍| 中文精品一卡2卡3卡4更新| 2021少妇久久久久久久久久久| 直男gayav资源| 91午夜精品亚洲一区二区三区| 2018国产大陆天天弄谢| 午夜爱爱视频在线播放| h日本视频在线播放| 亚洲伊人久久精品综合| 国语对白做爰xxxⅹ性视频网站| 欧美变态另类bdsm刘玥| 亚洲欧美精品专区久久| 日本爱情动作片www.在线观看| 亚洲在线自拍视频| 国内少妇人妻偷人精品xxx网站| 国产亚洲午夜精品一区二区久久 | 亚洲av电影不卡..在线观看| 欧美不卡视频在线免费观看| 精品一区在线观看国产| 欧美xxxx性猛交bbbb| 大香蕉久久网| 日韩欧美三级三区| 欧美日本视频| 国产av在哪里看| videossex国产| 免费黄频网站在线观看国产| 精品人妻视频免费看| 在线 av 中文字幕| 白带黄色成豆腐渣| 少妇丰满av| 成人午夜高清在线视频| 亚洲精品国产av蜜桃| 天堂av国产一区二区熟女人妻| 久久久久久久久久成人| av免费观看日本| 久久久久久国产a免费观看| 乱系列少妇在线播放| 别揉我奶头 嗯啊视频| 午夜免费男女啪啪视频观看| 人人妻人人澡欧美一区二区| 亚洲av免费在线观看| 两个人视频免费观看高清| 99久国产av精品国产电影| 亚洲不卡免费看| 亚洲成人av在线免费| 日韩欧美精品v在线| 麻豆国产97在线/欧美| 国产美女午夜福利| 中文字幕久久专区| 全区人妻精品视频| 亚洲熟妇中文字幕五十中出| 性色avwww在线观看| 亚洲美女搞黄在线观看| 日韩欧美精品v在线| 一本久久精品| h日本视频在线播放| 久久久精品94久久精品| 国产乱人偷精品视频| 欧美bdsm另类| 久久久久精品性色| 免费看光身美女| 久久综合国产亚洲精品| 天堂网av新在线| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | .国产精品久久| 老司机影院毛片| 免费看a级黄色片| 91久久精品国产一区二区三区| 69人妻影院| 精品久久久精品久久久| 久99久视频精品免费| 欧美日韩亚洲高清精品| 一级毛片aaaaaa免费看小| www.av在线官网国产| 日产精品乱码卡一卡2卡三| 我的女老师完整版在线观看| av在线天堂中文字幕| 毛片女人毛片| 真实男女啪啪啪动态图| 97超视频在线观看视频| 特大巨黑吊av在线直播| 又大又黄又爽视频免费| 熟女人妻精品中文字幕| 亚洲精品视频女| 插阴视频在线观看视频| 色网站视频免费| 舔av片在线| 免费av观看视频| 美女黄网站色视频| 久久99热6这里只有精品| 国产极品天堂在线| 久久国内精品自在自线图片| 夜夜看夜夜爽夜夜摸| 日日干狠狠操夜夜爽| 少妇的逼水好多| 最新中文字幕久久久久| 久久精品夜色国产| 亚洲第一区二区三区不卡| 97在线视频观看| 亚洲av中文字字幕乱码综合| 国产 亚洲一区二区三区 | 性色avwww在线观看| 日日摸夜夜添夜夜添av毛片| a级毛片免费高清观看在线播放| 秋霞在线观看毛片| 不卡视频在线观看欧美| 免费av毛片视频| 亚洲国产最新在线播放| 成人综合一区亚洲| 国产精品国产三级专区第一集| 六月丁香七月| 婷婷色av中文字幕| 久久久久国产网址| 免费播放大片免费观看视频在线观看| 一区二区三区乱码不卡18| 久久精品国产鲁丝片午夜精品| 波野结衣二区三区在线| 国产精品久久久久久精品电影| 能在线免费看毛片的网站| 国产在线一区二区三区精| 美女脱内裤让男人舔精品视频| 日本wwww免费看| 一区二区三区免费毛片| 又黄又爽又刺激的免费视频.| 亚洲国产高清在线一区二区三| 美女黄网站色视频| 日韩 亚洲 欧美在线| 女的被弄到高潮叫床怎么办| 久久久午夜欧美精品| 爱豆传媒免费全集在线观看| 51国产日韩欧美| 97超碰精品成人国产| 日日撸夜夜添| 国产淫片久久久久久久久| 国产欧美日韩精品一区二区| 国产亚洲5aaaaa淫片| 一级二级三级毛片免费看| 免费看日本二区| 人人妻人人澡人人爽人人夜夜 | 丰满少妇做爰视频| 国产精品不卡视频一区二区| 青青草视频在线视频观看| 亚洲在线观看片| 草草在线视频免费看| 男人和女人高潮做爰伦理| 国产一级毛片七仙女欲春2| 日日摸夜夜添夜夜爱| 免费av观看视频| 欧美日韩精品成人综合77777| 69av精品久久久久久| 国产欧美日韩精品一区二区| 国产精品一及| av女优亚洲男人天堂| 男女边摸边吃奶| 男插女下体视频免费在线播放| 欧美激情久久久久久爽电影| 久久午夜福利片| 国产伦在线观看视频一区| 午夜激情久久久久久久| 特大巨黑吊av在线直播| 亚洲一级一片aⅴ在线观看| 嘟嘟电影网在线观看| 天堂中文最新版在线下载 | 一级毛片 在线播放| 国产91av在线免费观看| 亚洲乱码一区二区免费版| 国语对白做爰xxxⅹ性视频网站| 免费观看性生交大片5| 日韩视频在线欧美| or卡值多少钱| 成人午夜精彩视频在线观看| 观看美女的网站| 在线播放无遮挡| 亚洲国产精品成人综合色| 亚洲最大成人中文| 春色校园在线视频观看| 国产精品久久久久久久电影| 水蜜桃什么品种好| 两个人的视频大全免费| 久久久久性生活片| 亚洲精品视频女| 极品教师在线视频| 久久99蜜桃精品久久| 干丝袜人妻中文字幕| 国产精品国产三级专区第一集| 蜜桃亚洲精品一区二区三区| 免费看a级黄色片| 亚洲精品456在线播放app| 亚洲自拍偷在线| 成年av动漫网址| 少妇被粗大猛烈的视频| 99热这里只有是精品在线观看| 成人国产麻豆网| 成人二区视频| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 久久精品久久久久久噜噜老黄| av黄色大香蕉| 国产亚洲午夜精品一区二区久久 | 亚洲精品一二三| 校园人妻丝袜中文字幕| 青春草视频在线免费观看| 99九九线精品视频在线观看视频| 亚洲精华国产精华液的使用体验| 最近中文字幕高清免费大全6| 熟妇人妻久久中文字幕3abv| 亚洲婷婷狠狠爱综合网| 大陆偷拍与自拍| 国产麻豆成人av免费视频| 91久久精品电影网| 国产69精品久久久久777片| 日韩一本色道免费dvd| 色播亚洲综合网| 久久久久久久久久久免费av| 97在线视频观看| 中文字幕亚洲精品专区| 亚洲av中文av极速乱| 女人十人毛片免费观看3o分钟| 日韩欧美一区视频在线观看 | 免费大片18禁| 极品少妇高潮喷水抽搐| 中文欧美无线码| 午夜福利视频1000在线观看| 免费观看在线日韩| 综合色丁香网| 成人综合一区亚洲| 亚洲av成人精品一二三区| 中文字幕免费在线视频6| 男女啪啪激烈高潮av片| 嫩草影院入口| 亚洲在久久综合| 91在线精品国自产拍蜜月| 日韩一区二区视频免费看| 国产亚洲午夜精品一区二区久久 | 欧美另类一区| 在线观看人妻少妇| 精品午夜福利在线看| 国产色婷婷99| 一夜夜www| 一级毛片aaaaaa免费看小| 久久6这里有精品| 国产成人freesex在线| 欧美激情在线99| 国产精品1区2区在线观看.| 一级毛片我不卡| 国产探花极品一区二区| 嫩草影院入口| 午夜福利在线在线| 亚洲一区高清亚洲精品| av国产免费在线观看| videos熟女内射| 久久久久久伊人网av| 18+在线观看网站| 国产高清有码在线观看视频| 亚洲无线观看免费| 99久久精品国产国产毛片| 三级国产精品欧美在线观看| 久久久色成人| 久久精品国产自在天天线| 国产91av在线免费观看| 男女国产视频网站| 久久久久免费精品人妻一区二区| 国产精品精品国产色婷婷| 欧美区成人在线视频| 久久午夜福利片| 成人亚洲欧美一区二区av| 免费在线观看成人毛片| 2021天堂中文幕一二区在线观| 成年免费大片在线观看| 国产欧美另类精品又又久久亚洲欧美| 精品国内亚洲2022精品成人| av又黄又爽大尺度在线免费看| 久久鲁丝午夜福利片| 成人午夜高清在线视频| 成人高潮视频无遮挡免费网站| 国产精品熟女久久久久浪| 寂寞人妻少妇视频99o| 性色avwww在线观看| 欧美日韩一区二区视频在线观看视频在线 | 国产精品.久久久| 好男人视频免费观看在线| 亚洲乱码一区二区免费版| 特大巨黑吊av在线直播| 日本wwww免费看| 老师上课跳d突然被开到最大视频| 国产亚洲午夜精品一区二区久久 | 看十八女毛片水多多多| 欧美性感艳星| 精品一区二区三区人妻视频| 在线免费十八禁| 免费观看a级毛片全部| 亚洲人成网站在线播| 一级毛片 在线播放| 免费看光身美女| 亚洲婷婷狠狠爱综合网| 日韩av在线免费看完整版不卡| 男人狂女人下面高潮的视频| av免费观看日本| 亚洲人成网站在线观看播放| 欧美一区二区亚洲| 一个人看视频在线观看www免费| 人妻少妇偷人精品九色| 国内少妇人妻偷人精品xxx网站| av黄色大香蕉| 日韩欧美三级三区| 精品一区二区三卡| 嫩草影院入口| 午夜福利在线观看免费完整高清在| 男人爽女人下面视频在线观看| 久久精品夜夜夜夜夜久久蜜豆| 伊人久久国产一区二区| 亚洲av男天堂| 一个人看视频在线观看www免费| 亚洲精品影视一区二区三区av|