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

    土體率效應(yīng)對(duì)動(dòng)力錨沉貫深度影響

    2017-02-09 05:23:20君,治,
    關(guān)鍵詞:板形剪應(yīng)變魚雷

    劉 君, 李 明 治, 韓 聰 聰

    ( 1.大連理工大學(xué) 海岸和近海工程國(guó)家重點(diǎn)實(shí)驗(yàn)室, 遼寧 大連 116024;2.上??睖y(cè)設(shè)計(jì)研究院有限公司, 上海 200434 )

    土體率效應(yīng)對(duì)動(dòng)力錨沉貫深度影響

    劉 君*1, 李 明 治1,2, 韓 聰 聰1

    ( 1.大連理工大學(xué) 海岸和近海工程國(guó)家重點(diǎn)實(shí)驗(yàn)室, 遼寧 大連 116024;2.上??睖y(cè)設(shè)計(jì)研究院有限公司, 上海 200434 )

    動(dòng)力安裝錨(簡(jiǎn)稱動(dòng)力錨)是一種依靠自重進(jìn)行安裝、適用于深海錨泊系統(tǒng)的新型錨固基礎(chǔ).在動(dòng)力沉貫過(guò)程中,土體經(jīng)受很高的剪應(yīng)變率(約25 s-1),因此必須考慮土體率效應(yīng)影響.總結(jié)了以往關(guān)于土體率效應(yīng)的研究成果,整理了常用3種率效應(yīng)模型中參數(shù)的取值范圍,采用運(yùn)動(dòng)微分方程研究了不同率效應(yīng)模型參數(shù)對(duì)動(dòng)力錨沉貫深度的影響.結(jié)果表明:土體率效應(yīng)會(huì)造成動(dòng)力錨的沉貫深度顯著減小,可偏小30%~40%;當(dāng)率效應(yīng)參數(shù)較大時(shí),錨的沉貫深度對(duì)參考應(yīng)變率的變化比較敏感;摩擦項(xiàng)與端承項(xiàng)的率效應(yīng)比值對(duì)動(dòng)力錨的沉貫深度影響較大,需要進(jìn)一步研究該比值的確切取值方法;給出了可以滿足任意形狀和深度的矩形基礎(chǔ)承載力計(jì)算公式.

    海洋基礎(chǔ);動(dòng)力錨;率效應(yīng);沉貫深度;承載力系數(shù)

    0 引 言

    動(dòng)力安裝錨(dynamically installed anchors,DIAs)簡(jiǎn)稱動(dòng)力錨,是一種依靠自重進(jìn)行安裝、適用于深海錨泊系統(tǒng)的新型錨固基礎(chǔ)[1].在安裝時(shí),首先將錨釋放至距離海床表面一定高度處(50~150 m),然后松開安裝錨鏈,錨在自重作用下貫入海床中.動(dòng)力錨需要較高的貫入速率(15~25 m/s)以達(dá)到較深的沉貫深度,提供足夠的抗拔承載力.動(dòng)力錨的沉貫過(guò)程是一個(gè)土體經(jīng)受高剪應(yīng)變率的過(guò)程,因此必須考慮土體率效應(yīng)影響.

    土體率效應(yīng)是指土的抗剪強(qiáng)度隨剪應(yīng)變率提高而增大的現(xiàn)象[2].土體率效應(yīng)大小取決于貫入物體的形狀和尺寸、貫入速度和土體的密度、滲透性及超固結(jié)比(OCR)等特性.大部分海床土為低滲透性黏土或鈣質(zhì)土,所以對(duì)于動(dòng)力錨的沉貫過(guò)程,必須考慮率效應(yīng)導(dǎo)致的土體強(qiáng)度增加效應(yīng).土體率效應(yīng)對(duì)動(dòng)力錨的沉貫深度影響程度如何,不同的率效應(yīng)模型又有怎樣的影響目前還不得而知.

    本文首先總結(jié)以往關(guān)于土體率效應(yīng)的研究結(jié)果,整理常用3種率效應(yīng)模型中的參數(shù)取值范圍.在此基礎(chǔ)上,采用基于牛頓第二運(yùn)動(dòng)微分方程研究率效應(yīng)參數(shù)對(duì)動(dòng)力錨沉貫深度的影響,旨在揭示動(dòng)力錨貫入過(guò)程中率效應(yīng)對(duì)沉貫深度的影響范圍.

    1 率效應(yīng)模型

    表征率效應(yīng)的主要模型有3種:冪函數(shù)型、半對(duì)數(shù)型和反雙曲正弦型,如下式所示:

    (1a)

    (1b)

    (1c)

    (2)

    參考剪應(yīng)變率的取值應(yīng)考慮兩方面的因素:首先,參考剪應(yīng)變率應(yīng)該足夠大以確保土體處于不排水狀態(tài);其次在保證土體處于不排水狀態(tài)的前提下,參考剪應(yīng)變率應(yīng)足夠小,盡可能使剪應(yīng)變率與參考剪應(yīng)變率的比值有2~3個(gè)量級(jí)的變化量,以便對(duì)比剪應(yīng)變率對(duì)率效應(yīng)的影響.

    2 已有的率效應(yīng)研究成果

    確定土體強(qiáng)度的試驗(yàn)方法主要有三軸試驗(yàn)、錐形貫入儀試驗(yàn)(CPT)以及十字板剪切試驗(yàn).近年來(lái),T形和球形貫入儀及動(dòng)態(tài)錐形貫入儀(FFP)也越來(lái)越多地用于測(cè)量土體強(qiáng)度.不同類型試驗(yàn)用到的貫入儀形狀、剛度、粗糙度等也不同,得到的率效應(yīng)大小有所差別[2-4],但差別不大.海床土大多為黏性土[5],所以本文主要總結(jié)了黏土和鈣質(zhì)土土體率效應(yīng)試驗(yàn)結(jié)果.

    2.1 三軸試驗(yàn)

    因三軸試驗(yàn)可嚴(yán)格控制土樣的排水狀態(tài)和外加荷載,所以其能精確測(cè)量土的不排水抗剪強(qiáng)度.表1總結(jié)了以往的三軸試驗(yàn)結(jié)果.從表1可以看出,剪應(yīng)變率每增長(zhǎng)一個(gè)量級(jí),土強(qiáng)度的增長(zhǎng)范圍為8%~20%.根據(jù)率效應(yīng)的定義,可反算得出率效應(yīng)參數(shù)β、λ和λ′的范圍分別為0.034~0.079、0.080~0.200和0.027~0.067.

    2.2 十字板剪切試驗(yàn)

    十字板剪切試驗(yàn)耗時(shí)短、操作簡(jiǎn)單,在確定原位黏土不排水抗剪強(qiáng)度試驗(yàn)中得到了廣泛應(yīng)用[11].表2列出了幾種有代表性的試驗(yàn)結(jié)果.

    從表2可以看出,除了Schlue等[14]的結(jié)果,剪應(yīng)變率每增長(zhǎng)一個(gè)量級(jí)對(duì)應(yīng)的土強(qiáng)度增長(zhǎng)范圍為10%~30%.有一點(diǎn)需要說(shuō)明的是,土體的含水量越高,其峰值強(qiáng)度的率效應(yīng)會(huì)越明顯.Schlue等[14]的試驗(yàn)探測(cè)了3種含水量的土樣,含水量分別是71%、112%、177%.對(duì)于含水量為177%的土樣,其流動(dòng)機(jī)制近乎于黏滯流體,因此得到的率效應(yīng)參數(shù)比其他結(jié)果要大出很多.另一點(diǎn)需要說(shuō)明的是,恒定的率效應(yīng)參數(shù)(β、λ、λ′)可能會(huì)低估高剪應(yīng)變率情況下土體的率效應(yīng)[15].這一點(diǎn)在O′Loughlin等[17]關(guān)于魚雷錨動(dòng)力沉貫的離心模型試驗(yàn)中也得到了印證,率效應(yīng)參數(shù)會(huì)隨沉貫速度的增大而增大.因此,當(dāng)剪應(yīng)變率的變化范圍超過(guò)4個(gè)量級(jí)時(shí),恒定的率效應(yīng)參數(shù)可能不再適合描述土體率效應(yīng)[15-16].

    表1 三軸試驗(yàn)得到的率效應(yīng)研究結(jié)果

    表2 十字板剪切試驗(yàn)得到的率效應(yīng)研究結(jié)果

    2.3 T形貫入儀和球形貫入儀試驗(yàn)結(jié)果

    T形和球形貫入儀作為新型測(cè)量土體不排水抗剪強(qiáng)度的貫入儀,在海洋工程領(lǐng)域應(yīng)用前景廣泛.當(dāng)沉貫深度足夠大時(shí),貫入儀周圍土體呈全流機(jī)制,土體塑性區(qū)只在貫入儀附近的區(qū)域發(fā)展,此時(shí)承載力系數(shù)為常數(shù).關(guān)于T形和球形貫入儀確定土體率效應(yīng)的試驗(yàn)如表3所示.

    如前文所述,同一種土樣的率效應(yīng)會(huì)隨著試驗(yàn)方法的不同而出現(xiàn)偏差.表3中同一種土樣(Burswood黏土)通過(guò)T形和球形貫入儀測(cè)出的率效應(yīng)參數(shù)有所不同[20].從表3中可知,剪應(yīng)變率每增加一個(gè)量級(jí)對(duì)應(yīng)的土強(qiáng)度增長(zhǎng)量范圍是10%~38%.

    2.4 其他試驗(yàn)

    除了上述試驗(yàn)方法,還有一些其他試驗(yàn)方法可以確定土體率效應(yīng),如錐形貫入試驗(yàn)、單剪試驗(yàn)以及FFP試驗(yàn),相關(guān)結(jié)果如表4所示.錐形貫入儀側(cè)壁摩阻項(xiàng)對(duì)應(yīng)的率效應(yīng)參數(shù)明顯高于端承阻力項(xiàng)的率效應(yīng)參數(shù),因?yàn)樨炄雰x側(cè)壁附近的剪切帶要比錐尖處更窄,對(duì)應(yīng)的剪應(yīng)變率更大[24].根據(jù)Einav等[24]的研究成果,表征率效應(yīng)的冪函數(shù)公式(1a)可改寫為

    (3)

    式中:μ是一個(gè)量綱為一的量,端承阻力項(xiàng)和側(cè)壁摩阻項(xiàng)的率效應(yīng)對(duì)應(yīng)的μ分別取1和(1/β-1).

    從表4中可以得出,除了Dayal等[16]的結(jié)果外,剪應(yīng)變率每增加一個(gè)量級(jí)對(duì)應(yīng)的土強(qiáng)度提高12%~35%(端承阻力項(xiàng)).在Dayal等[16]的試驗(yàn)中,土強(qiáng)度范圍為3~81 kPa,結(jié)果表明強(qiáng)度越低,率效應(yīng)越明顯.強(qiáng)度為3 kPa的土樣對(duì)應(yīng)端承阻力項(xiàng)的Δsu=150%;而強(qiáng)度在9~81 kPa的土樣對(duì)應(yīng)的Δsu=25%~3%.這一低強(qiáng)度土體率效應(yīng)顯著增長(zhǎng)的現(xiàn)象與Schlue等[14]的試驗(yàn)結(jié)果相似.這可能是對(duì)于某一種黏土,存在一個(gè)強(qiáng)度臨界值,當(dāng)土體強(qiáng)度低于該臨界值時(shí),土體會(huì)呈現(xiàn)明顯的黏滯性.

    表3 T形和球形貫入儀試驗(yàn)得到的率效應(yīng)研究結(jié)果

    表4 土體率效應(yīng)的其他類型試驗(yàn)研究結(jié)果

    2.5 率效應(yīng)參數(shù)取值范圍匯總

    綜上所述,剪應(yīng)變率每增加一個(gè)量級(jí)對(duì)應(yīng)的土強(qiáng)度增長(zhǎng)范圍是8%~38%.影響率效應(yīng)的因素可歸為以下幾類:

    (1)不同的試驗(yàn)方法由于貫入儀形狀、尺寸、粗糙度等不同得到的率效應(yīng)參數(shù)會(huì)有所不同;

    (2)低強(qiáng)度和低超固結(jié)比的土體可能會(huì)表現(xiàn)出更明顯的率效應(yīng)[6,12-13,15];

    (3)當(dāng)剪應(yīng)變率改變量超過(guò)4個(gè)量級(jí)時(shí),恒定的率效應(yīng)參數(shù)可能會(huì)低估土體率效應(yīng)[15-16];

    (4)可能存在一個(gè)臨界強(qiáng)度值,當(dāng)土強(qiáng)度低于該臨界值時(shí),土體率效應(yīng)會(huì)顯著提高.

    通常海洋工程中采用的T形貫入儀尺寸為L(zhǎng)T×DT=250 mm×40 mm(LT為T形桿的長(zhǎng)度,DT為T形桿的直徑),球形貫入儀直徑為113 mm,錐形貫入儀直徑為36~44 mm[25].如果3種貫入儀的貫入速度均為20 mm/s,則貫入儀的應(yīng)變率范圍在0.18 s-1到0.5 s-1.在動(dòng)力錨沉貫問(wèn)題數(shù)值模擬中,推薦的參考應(yīng)變率為0.1 s-1[26].綜合表1~4的參考應(yīng)變率,實(shí)際工程中貫入儀的應(yīng)變率以及本文所要研究動(dòng)力錨的速度及等效直徑等因素,本文選取的參考應(yīng)變率分別為0.01、0.1、0.5 s-1,盡可能使動(dòng)力錨的剪應(yīng)變率與參考剪應(yīng)變率的比值有2~3個(gè)量級(jí)的變化量,但不超過(guò)4個(gè)量級(jí).根據(jù)上述總結(jié),不同率效應(yīng)模型對(duì)應(yīng)的率效應(yīng)參數(shù)的取值范圍如表5所示.

    表5 應(yīng)變率參數(shù)及參考剪應(yīng)變率的取值范圍

    3 率效應(yīng)對(duì)動(dòng)力錨沉貫深度的影響

    3.1 動(dòng)力錨沉貫深度預(yù)測(cè)程序的驗(yàn)證

    本文僅僅研究土體率效應(yīng)對(duì)錨沉貫深度的影響,故在忽略錨鏈的附加質(zhì)量及錨鏈與水相互作用的情況下,根據(jù)True[2]的分析,動(dòng)力錨在沉貫時(shí)的受力如圖1所示.

    (a) 魚雷錨

    (b) 板形重力安裝錨

    圖1 動(dòng)力錨沉貫時(shí)受到的作用力

    Fig.1 Forces acting on dynamically installed anchors during penetrating in clay

    依據(jù)牛頓第二定律,動(dòng)力錨沉貫深度的微分方程如下式所示:

    (4)

    式中:m是錨的質(zhì)量,t是時(shí)間,z是錨在t時(shí)刻的沉貫深度,Ws是錨在水中的重量,F(xiàn)b是錨的垂直于錨軸線方向的橫截面受到的端承阻力,F(xiàn)f是錨側(cè)壁受到的摩阻力,F(xiàn)γ是土體對(duì)錨的浮力,F(xiàn)d是錨受到的拖曳阻力.

    Fb=Ncsu0Ab

    Ff=αsu0Af

    Fγ=γ′V

    Fd=0.5CdρApv2

    其中Nc是錨的橫截面對(duì)應(yīng)形狀的基礎(chǔ)承載力系數(shù),對(duì)于魚雷錨和板形重力安裝錨的扁平狀翼板,可取Skempton[27]給出的深埋條形基礎(chǔ)的承載力系數(shù)Nc=7.5,對(duì)于魚雷錨的錨軸,可按照Ma等[28]給出的CPT的承載力系數(shù)公式,取為12.0;su0是土的靜態(tài)不排水剪切強(qiáng)度;Ab是t時(shí)刻錨沒(méi)入土中部分的橫截面積.

    Af是錨的表面積;α是錨土界面摩擦因數(shù),通常取為土體靈敏度的倒數(shù).

    γ′是土的浮重度,V是錨擠開土的體積.

    Cd是拖曳力系數(shù),通常錨的長(zhǎng)細(xì)比越大,形狀越復(fù)雜,Cd值越大,根據(jù)?ye[29]給出的研究結(jié)果,魚雷錨的Cd取0.65,根據(jù)Kim等[26]給出的結(jié)果,板形重力安裝錨的Cd取0.63;ρ是土體的浮密度;Ap是錨垂直于軸線方向的投影面積;v是錨在t時(shí)刻的速度.

    Rf1和Rf2分別表征土體率效應(yīng)在錨的端承阻力項(xiàng)和側(cè)壁摩阻項(xiàng)中的增強(qiáng)效應(yīng).Rf1的計(jì)算如式(1)、(2)所示.在計(jì)算魚雷錨和板形重力安裝錨的沉貫深度時(shí),式(2)中的d分別代表魚雷錨的錨軸直徑和板形重力安裝錨的等效直徑(與橫截面積等效的圓的直徑).雖然式(3)給出了冪函數(shù)形式的Rf1和Rf2的不同,且近幾年一些學(xué)者[22-23]通過(guò)FFP試驗(yàn)得出Rf2/Rf1的取值范圍在1.0~2.0,但具體的比值尚無(wú)定論.本文的主旨并非探討動(dòng)力錨的絕對(duì)沉貫深度,因此Rf2/Rf1的具體數(shù)值就顯得不是特別重要.所以如無(wú)特殊說(shuō)明,以下算例的Rf2/Rf1均根據(jù)O′Loughlin等[17]的取值方法取為1.0.

    許多學(xué)者在式(4)基礎(chǔ)上對(duì)動(dòng)力錨的沉貫深度進(jìn)行了研究,其中O′Loughlin等[17]依據(jù)離心機(jī)模型試驗(yàn)并結(jié)合式(4)做了較多關(guān)于魚雷錨沉貫深度的研究工作.本文首先根據(jù)式(4)編制了計(jì)算動(dòng)力錨沉貫深度的Matlab程序,模擬了O′Loughlin 等[17]的兩種魚雷錨的離心機(jī)模型試驗(yàn),并進(jìn)行了對(duì)比,以此驗(yàn)證本文所編制程序的正確性.從圖2所示對(duì)比結(jié)果可以看到兩者的吻合程度很好,表明本文所編制程序是準(zhǔn)確的.然后根據(jù)前文總結(jié)的土體率效應(yīng)參數(shù)的取值范圍,使用該程序分析了土體率效應(yīng)參數(shù)對(duì)魚雷錨和板形重力安裝錨沉貫深度的影響.

    (a) 無(wú)尾翼魚雷錨

    (b) 有尾翼魚雷錨

    圖2 魚雷錨沉貫深度

    Fig.2 Penetration depth of torpedo anchors

    3.2 率效應(yīng)參數(shù)對(duì)魚雷錨沉貫深度的影響

    魚雷錨的形狀及幾何參數(shù)如圖3(a)所示,錨的質(zhì)量為98 t,錨的初始沉貫速度是20 m/s,土的不排水抗剪強(qiáng)度su0=(5+2z) kPa,土的浮容重為6 kN/m3.

    (a) 魚雷錨 (b) 板形重力安裝錨

    從表6可以看出,土體率效應(yīng)參數(shù)(β、λ、λ′)對(duì)魚雷錨的沉貫深度影響很大,甚至可達(dá)30%以上.錨土之間的摩擦因數(shù)α對(duì)魚雷錨沉貫深度影響也較大,從其變化的幅值來(lái)看可達(dá)27.8%.因此魚雷錨的沉貫過(guò)程分析必須細(xì)致考慮土體率效應(yīng)和錨土界面摩擦效應(yīng)對(duì)錨的沉貫深度的影響.

    表6 魚雷錨沉貫深度變化范圍

    (a) ref=0.1 s-1 (b) β=0.05 (c) α=0.33

    前文述及,由于魚雷錨的翼板厚寬比接近0.1,很多學(xué)者將翼板受到的端承阻力按照深埋條形基礎(chǔ)的承載力進(jìn)行計(jì)算,并使用Skempton[27]給出的承載力系數(shù)Nc=7.5計(jì)算動(dòng)力錨沉貫過(guò)程中翼板受到的端承阻力.

    然而通過(guò)梳理近20年來(lái)其他學(xué)者[30-32]的研究成果發(fā)現(xiàn):深埋條形基礎(chǔ)的承載力系數(shù)并非7.5.如圖5所示,Skempton[27]認(rèn)為條形基礎(chǔ)埋深超過(guò)4倍基礎(chǔ)寬度時(shí),即為深埋情況,承載力系數(shù)達(dá)到7.5后不再增長(zhǎng).而Salgado等[30]、Edwards等[31]以及Quoc[32]給出的結(jié)果是:條形基礎(chǔ)的埋深達(dá)到了5倍基礎(chǔ)寬度時(shí),其承載力系數(shù)仍然有增長(zhǎng)的趨勢(shì),且高出7.5較多.因此本文采用ABAQUS有限元軟件計(jì)算了長(zhǎng)寬比為0.05、0.10 的矩形基礎(chǔ)以及條形基礎(chǔ)在較大埋深(最大埋深達(dá)到了150倍基礎(chǔ)寬度)情況下的承載力系數(shù).

    圖5 不同埋深情況的條形基礎(chǔ)承載力系數(shù)

    由于對(duì)稱性,采用1/4模型計(jì)算,有限元網(wǎng)格如圖6所示.土體采用理想彈塑性材料,滿足Tresca屈服準(zhǔn)則;基礎(chǔ)的剛度通常遠(yuǎn)大于土體,因此將基礎(chǔ)視為剛體.通過(guò)網(wǎng)格收斂性分析,采用B/50(B為基礎(chǔ)寬度)作為土體最小網(wǎng)格尺寸,既能滿足精度要求又能有效控制計(jì)算時(shí)間.基礎(chǔ)棱角處的土體的塑性應(yīng)變較大,因此棱角處土體網(wǎng)格密度較大,而遠(yuǎn)離該區(qū)域的土體網(wǎng)格密度較?。?jīng)過(guò)比較分析,土體單元類型選擇一階完全積分的六面體雜交單元(C3D8H).土體與基礎(chǔ)之間接觸光滑.土體底面為鉸支約束,側(cè)壁為滾軸約束.加載方式為給基礎(chǔ)一個(gè)向下的位移荷載,直至土體的抗力達(dá)到一個(gè)屈服平臺(tái)為止.圖7所示為B/L= 0.05,D/B= 10情況下承載力系數(shù)隨著基礎(chǔ)量綱一位移的變化情況(L為基礎(chǔ)長(zhǎng)度,D為基礎(chǔ)上表面的埋深).圖中的橫坐標(biāo)為基礎(chǔ)的量綱一位移,E為土體的彈性模量,δ為基礎(chǔ)的位移.

    所有工況的承載力系數(shù)結(jié)果如圖8中的散點(diǎn)所示.依據(jù)計(jì)算結(jié)果采用Origin擬合了一個(gè)計(jì)算薄板基礎(chǔ)在較大埋深情況下的承載力系數(shù)計(jì)算公式,如下式所示:

    Nc=5.14×{1.022+0.063B/L+c1ln(1+c2(D/B))}

    (5)

    其中

    圖8 矩形基礎(chǔ)的極限承載力系數(shù)

    3.3 率效應(yīng)參數(shù)對(duì)板形重力安裝錨沉貫深度的影響

    板形重力安裝錨形狀及幾何參數(shù)如圖3(b)所示.錨質(zhì)量為39.8 t,錨的初始沉貫速度是20 m/s.土體不排水抗剪強(qiáng)度su0=(2.4+1.1)zkPa,土體浮容重為6 kN/m3.率效應(yīng)參數(shù)的取值范圍及對(duì)沉貫深度的影響如表7所示,表中的沉貫深度變化范圍的計(jì)算方法與魚雷錨的相同.

    從表7可以看出,土體率效應(yīng)參數(shù)(β、λ、λ′)對(duì)板形重力安裝錨的沉貫深度影響更大,甚至可達(dá)40%以上.錨土之間的摩擦因數(shù)α對(duì)板形重力安裝錨沉貫深度影響也較大,從其變化的幅值來(lái)看可達(dá)22.0%.因此板形重力安裝錨的沉貫過(guò)程分析必須細(xì)致考慮土體率效應(yīng)和錨土界面摩擦效應(yīng)對(duì)錨的沉貫深度的影響.

    表7 板形重力安裝錨沉貫深度變化范圍

    (a) ref=0.10 s-1

    (b)λ=0.14

    (c)α=0.33

    圖9 率效應(yīng)參數(shù)的相互影響(板形重力安裝錨)

    Fig.9 The mutual effect of strain-rate effect parameters (gravity installed plate anchor)

    4 結(jié) 論

    (1)不論是魚雷錨,還是板形重力安裝錨,土體率效應(yīng)均會(huì)造成錨的沉貫深度偏小31%~42.1%,所以在評(píng)估魚雷錨和板形重力安裝錨的沉貫深度時(shí),必須考慮率效應(yīng).

    (2)當(dāng)率效應(yīng)參數(shù)(β、λ、λ′)較大時(shí),錨的沉貫深度對(duì)參考應(yīng)變率的變化比較敏感.因此當(dāng)率效應(yīng)參數(shù)較大時(shí),需謹(jǐn)慎選擇參考應(yīng)變率.

    (3)摩擦項(xiàng)和端承項(xiàng)的率效應(yīng)比值Rf2/Rf1對(duì)魚雷錨和板形重力安裝錨的沉貫深度影響較大,但目前關(guān)于Rf2/Rf1的取值尚無(wú)定論,因此需要進(jìn)一步研究Rf2/Rf1的確切取值方法.

    (4)目前常用的率效應(yīng)模型均可以很好地描述土體率效應(yīng),在本文所研究的參數(shù)范圍內(nèi),3種模型給出基本一致的影響.因此選擇任何一個(gè)模型來(lái)描述土體率效應(yīng)對(duì)沉貫深度的影響均可.但計(jì)算翼板的端承阻力時(shí),采用以往的條形基礎(chǔ)的承載力系數(shù)是偏于危險(xiǎn)的.

    (5)給出了可以滿足任意形狀和深度的矩形基礎(chǔ)承載力計(jì)算公式.

    [1] RICHARDSON M D, O′LOUGHLIN C D, RANDOLPH M F,etal. Setup following installation of dynamic anchors in normally consolidated clay [J]. Journal of Geotechnical and Geoenvironmental Engineering, 2009, 135(4):487-496.

    [2] TRUE D G. Undrained vertical penetration into ocean bottom soils [D]. Berkeley: University of California, 1976.

    [3] GRAHAM J, CROOKS J H, BELL A L. Time effects on the stress-strain behaviour of natural soft clays [J]. Géotechnique, 1983, 33(3):327-340.

    [4] LEHANE B M, O′LOUGHLIN C D, RANDOLPH M F,etal. Rate effects on penetrometer resistance in kaolin [J]. Géotechnique, 2009, 59(1):41-52.

    [5] RANDOLPH M, CASSIDY M, GOURVENEC S,etal. Challenges of offshore geotechnical engineering [C] // Proceedings of the 16th International Conference on Soil Mechanics and Geotechnical Engineering: Geotechnology in Harmony with the Global Environment. Rotterdam: Millpress Science Publishers, 2005:123-176.

    [6] SHEAHAN T C, LADD C C, GERMAINE J T. Rate-dependent undrained shear behavior of saturated clay [J]. Journal of Geotechnical Engineering-ASCE, 1996, 122(2):99-108.

    [7] YIN Z Y, HICHER P Y, RIOU Y,etal. An elasto-viscoplastic model for soft clay [J]. Geotechnical Special Publication, 2006(150):312-319.

    [8] YIN Zhenyu, KARSTUNEN M, HICHER P Y. Evaluation of the influence of elasto-viscoplastic scaling functions on modelling time-dependent behaviour of natural clays [J]. Soils and Foundations, 2010, 50(2):203-214.

    [9] TORISU S S, PEREIRA J M, DE GENNARO V,etal. Strain-rate effects in deep marine clays from the gulf of Guinea [J]. Géotechnique, 2012, 62(9):767-775.

    [10] PALIX E, WU H, CHAN N,etal. Liwan 3-1:How Deepwater Sediments from South China Sea Compare with Gulf of Guinea Sediments [C]. Houston: Offshore Technology Conference, 2013:OTC-24010-MS.

    [11] BISCONTIN G, PESTANA J M. Influence of peripheral velocity on vane shear strength of an artificial clay [J]. Geotechnical Testing Journal, 2001, 24(4):423-429.

    [12] AUBENY C P, SHI Han. Effect of rate-dependent soil strength on cylinders penetrating into soft clay [J]. IEEE Journal of Oceanic Engineering, 2007, 32(1):49-56.

    [13] BOUKPETI N, WHITE D, RANDOLPH M,etal. Characterization of the solid-fluid transition of fine-grained sediments [C] // Proceedings of the 28th International Conference on Ocean, Offshore and Arctic Engineering 2009, OMAE2009. New York: ASME, 2009:293-303.

    [14] SCHLUE B F, MOERZ T, KREITER S. Influence of shear rate on undrained vane shear strength of organic harbor mud [J]. Journal of Geotechnical and Geoenvironmental Engineering, 2010, 136(10):1437-1447.

    [15] ABELEV A, VALENT P. Strain-rate dependence of strength of the gulf of Mexico soft sediments [J]. IEEE Journal of Oceanic Engineering, 2013, 38(1):25-31.

    [16] DAYAL U, ALLEN J H. The effect of penetration rate on the strength of remolded clay and sand samples [J]. Canadian Geotechnical Journal, 1975, 12(3):336-348.

    [17] O′LOUGHLIN C D, RICHARDSON M D, RANDOLPH M F,etal. Penetration of dynamically installed anchors in clay [J]. Géotechnique, 2013, 63(11):909-919.

    [18] CHUNG S F, RANDOLPH M F, SCHNEIDER J A. Effect of penetration rate on penetrometer resistance in clay [J]. Journal of Geotechnical and Geoenvironmental Engineering, 2006, 132(9):1188-1196.

    [19] YAFRATE N J, DEJONG J T. Influence of penetration rate on measured resistance with full-flow penetrometers in soft clay [J]. Geotechnical Special Publication, 2007(173):1-10.

    [20] LOW H E, MAYNARD M L, RANDOLPH M F,etal. Geotechnical characterisation and engineering properties of Burswood clay [J]. Géotechnique, 2011, 61(7):575-591.

    [21] GANESAN S A, BOLTON M D. Characterisation of a high plasticity marine clay using a T-bar penetrometer [J]. Underwater Technology, 2013, 31(4):179-185.

    [22] STEINER A, KOPF A J, L′HEUREUX J,etal. In situ dynamic piezocone penetrometer tests in natural clayey soils - a reappraisal of strain-rate corrections [J]. Canadian Geotechnical Journal, 2014, 51(3):272-288.

    [23] CHOW S H, O′LOUGHLIN C D, RANDOLPH M F. Soil strength estimation and pore pressure dissipation for free-fall piezocone in soft clay [J]. Géotechnique, 2014, 64(10):817-827.

    [24] EINAV I, RANDOLPH M. Effect of strain rate on mobilised strength and thickness of curved shear bands [J]. Géotechnique, 2006, 56(7):501-504.

    [25] LOW H E, LUNNE T, ANDERSEN K H,etal. Estimation of intact and remoulded undrained shear strengths from penetration tests in soft clays [J]. Géotechnique, 2010, 60(11):843-859.

    [26] KIM Y H, HOSSAIN M S. Dynamic installation of OMNI-Max anchors in clay:numerical analysis [J]. Géotechnique, 2015, 65(12):1029-1037.

    [27] SKEMPTON A W. The bearing capacity of clays [C] // Proceedings of Building Research Congress. London: [s. n.], 1951:180-189.

    [28] MA Hongliang, ZHOU Mi, HU Yuxia,etal. Interpretation of layer boundaries and shear strengths for soft-stiff-soft clays using CPT data:LDFE analyses [J]. Journal of Geotechnical and Geoenvironmental Engineering, 2016, 142(1):04015055.

    [29] ?YE I. Simulation of trajectories for a deep penetrating anchor:Report No. 250:2000 [R]. Trondheim:CFD Norway, 2000.

    [30] SALGADO R, LYAMIN A V, SLOAN S W,etal. Two- and three-dimensional bearing capacity of foundations in clay [J]. Géotechnique, 2004, 54(5):297-306.

    [31] EDWARDS D H, ZDRAVKOVIC L, POTTS D M. Depth factors for undrained bearing capacity [J]. Géotechnique, 2005, 55(10):755-758.

    [32] QUOC V N. Numerical modelling of the undrained vertical bearing capacity of shallow foundations [D]. Toowoomba: University of Southern Queensland, 2008.

    Influence of soil strain-rate effect on embedment depth of dynamically installed anchors

    LIU Jun*1, LI Mingzhi1,2, HAN Congcong1

    ( 1.State Key Laboratory of Coastal and Offshore Engineering, Dalian University of Technology, Dalian 116024, China; 2.Shanghai Investigation, Design & Research Institute Company, Limited, Shanghai 200434, China )

    Dynamically installed anchors, which are installed through self-weight, are newly developed anchors applied to deep-water mooring systems. During the dynamical embedment procedure, the soil is subjected to high shear strain-rate (about 25 s-1). Hence, the soil strain-rate effect should be necessarily taken into consideration. The ranges of soil strain-rate parameters of three commonly used models are summarized based on the published literatures referred to soil strain-rate effect, and the influence of strain-rate effect on the embedment depth of dynamically installed anchors based on the motion differential equation is discussed in the following. The research results indicate that the embedment depth of dynamically installed anchors would drastically decrease (30%-40%) due to the soil strain-rate effect. The embedment depth is sensitive to the reference strain-rate when strain-rate parameters are in high values. The shear strain-rate effect ratio of the frictional resistance to the end bearing resistance has significant effect on the embedment depth and further investigation is necessitated to define this ratio. A new formula is proposed to calculate the bearing capacity of rectangular footing with any shape and embedment depth.

    offshore foundation; dynamically installed anchors; strain-rate effect; embedment depth; bearing capacity factor

    1000-8608(2017)01-0068-10

    2016-07-09;

    2016-11-20.

    國(guó)家自然科學(xué)基金資助項(xiàng)目(51479027,51539008).

    劉 君*(1972-),男,博士,教授,E-mail:junliu@dlut.edu.cn;李明治(1990-),男,碩士生,E-mail:731426627@qq.com;韓聰聰(1990-),女,博士生,E-mail:1059281811@qq.com.

    TU47

    A

    10.7511/dllgxb201701010

    猜你喜歡
    板形剪應(yīng)變魚雷
    魚雷造反
    考慮后續(xù)工序要求的軋機(jī)板形標(biāo)準(zhǔn)曲線研究
    基于Legendre多項(xiàng)式的板形模式識(shí)別優(yōu)化仿真及應(yīng)用
    軍事島 魚雷人
    改良隔震橡膠支座的試驗(yàn)及模擬分析*
    魚雷也瘋狂
    水泥改良黃土路基動(dòng)力穩(wěn)定性評(píng)價(jià)參數(shù)試驗(yàn)研究
    鄢家橋水庫(kù)土壩剪應(yīng)變及穩(wěn)定分析
    以板形為基礎(chǔ)的建筑物凍害防治措施
    利用GPS分析新疆地區(qū)形變特征
    科技視界(2015年27期)2015-10-08 10:43:12
    久久精品综合一区二区三区| 日韩欧美 国产精品| 亚洲国产成人一精品久久久| 免费大片18禁| 日本猛色少妇xxxxx猛交久久| 黑人高潮一二区| 久久久久久伊人网av| 男女下面进入的视频免费午夜| 免费播放大片免费观看视频在线观看| 日日啪夜夜撸| 国产欧美日韩精品一区二区| 精品亚洲乱码少妇综合久久| 最近最新中文字幕大全电影3| 性色av一级| 最近的中文字幕免费完整| 五月开心婷婷网| 中文资源天堂在线| av在线亚洲专区| 欧美高清成人免费视频www| 亚洲精品国产av蜜桃| 婷婷色综合www| 男人爽女人下面视频在线观看| av免费观看日本| 亚洲av欧美aⅴ国产| 中文字幕av成人在线电影| 少妇被粗大猛烈的视频| 亚洲欧美精品专区久久| 亚洲国产日韩一区二区| 久久热精品热| 中国三级夫妇交换| 观看美女的网站| av国产精品久久久久影院| 免费观看的影片在线观看| 韩国av在线不卡| 97在线人人人人妻| 亚洲天堂国产精品一区在线| 国产伦精品一区二区三区四那| 国产毛片在线视频| 美女脱内裤让男人舔精品视频| 白带黄色成豆腐渣| 亚洲欧美精品自产自拍| 国产精品无大码| 777米奇影视久久| 免费观看性生交大片5| 亚洲欧美一区二区三区国产| 亚洲av不卡在线观看| 亚洲av在线观看美女高潮| 国产精品国产av在线观看| 全区人妻精品视频| 国产真实伦视频高清在线观看| 高清av免费在线| 亚洲婷婷狠狠爱综合网| 国产精品一二三区在线看| 成年女人在线观看亚洲视频 | 日本免费在线观看一区| 黄片wwwwww| 亚洲国产精品专区欧美| 91在线精品国自产拍蜜月| 在线 av 中文字幕| 天天躁日日操中文字幕| 一本一本综合久久| 亚洲成人一二三区av| 亚洲三级黄色毛片| 亚洲精品一二三| 肉色欧美久久久久久久蜜桃 | 久久久色成人| 香蕉精品网在线| 国产伦精品一区二区三区四那| 男女无遮挡免费网站观看| 国产毛片在线视频| 亚洲经典国产精华液单| 国产免费福利视频在线观看| 日韩中字成人| 九九久久精品国产亚洲av麻豆| 久久影院123| 亚洲精品aⅴ在线观看| 秋霞在线观看毛片| 伊人久久国产一区二区| 免费人成在线观看视频色| 一个人看的www免费观看视频| 街头女战士在线观看网站| 久久精品久久久久久噜噜老黄| 成人漫画全彩无遮挡| 热99国产精品久久久久久7| 香蕉精品网在线| 亚洲在线观看片| 国产一区二区三区综合在线观看 | 国产精品国产三级国产专区5o| 国产爽快片一区二区三区| 国产成人a区在线观看| 成人高潮视频无遮挡免费网站| 色5月婷婷丁香| 一区二区三区精品91| 欧美少妇被猛烈插入视频| 国产成人a∨麻豆精品| videossex国产| 国产精品伦人一区二区| 欧美一级a爱片免费观看看| 亚洲精华国产精华液的使用体验| 国产高清不卡午夜福利| 国产黄频视频在线观看| 亚洲,欧美,日韩| 欧美区成人在线视频| 国产精品成人在线| 亚洲成色77777| 国产精品福利在线免费观看| 日韩大片免费观看网站| 观看美女的网站| 最近手机中文字幕大全| 美女视频免费永久观看网站| 亚洲三级黄色毛片| 大陆偷拍与自拍| 亚洲精品自拍成人| 少妇人妻久久综合中文| 91在线精品国自产拍蜜月| 免费大片18禁| 晚上一个人看的免费电影| 另类亚洲欧美激情| 免费观看av网站的网址| 黄片wwwwww| 毛片一级片免费看久久久久| 99热这里只有是精品在线观看| 亚洲在线观看片| 国产成人精品福利久久| 国产爱豆传媒在线观看| 最近最新中文字幕免费大全7| 18禁在线无遮挡免费观看视频| 最近2019中文字幕mv第一页| 熟妇人妻不卡中文字幕| 国产欧美亚洲国产| 免费高清在线观看视频在线观看| 欧美性猛交╳xxx乱大交人| 最近的中文字幕免费完整| 国产免费又黄又爽又色| 女人久久www免费人成看片| 3wmmmm亚洲av在线观看| 超碰av人人做人人爽久久| 青春草国产在线视频| 成年女人在线观看亚洲视频 | 高清欧美精品videossex| 极品教师在线视频| 亚洲丝袜综合中文字幕| av在线蜜桃| 免费电影在线观看免费观看| 一区二区三区乱码不卡18| 国产在线一区二区三区精| 一级a做视频免费观看| 亚洲av免费高清在线观看| 老师上课跳d突然被开到最大视频| 高清午夜精品一区二区三区| 色网站视频免费| 身体一侧抽搐| 97超视频在线观看视频| 日韩精品有码人妻一区| 18禁裸乳无遮挡动漫免费视频 | 国产日韩欧美在线精品| 久久久久久久久久久丰满| 久久精品国产a三级三级三级| 欧美激情国产日韩精品一区| 久久精品国产亚洲av天美| av专区在线播放| 日日摸夜夜添夜夜添av毛片| 别揉我奶头 嗯啊视频| 日韩在线高清观看一区二区三区| 男女那种视频在线观看| 国产高清三级在线| 免费黄频网站在线观看国产| 中文字幕制服av| 少妇猛男粗大的猛烈进出视频 | 观看免费一级毛片| 中文字幕制服av| 久久精品熟女亚洲av麻豆精品| 精品国产一区二区三区久久久樱花 | 成人亚洲精品一区在线观看 | 插逼视频在线观看| 一区二区av电影网| 亚洲av福利一区| 一区二区三区乱码不卡18| 久久6这里有精品| 少妇熟女欧美另类| 特级一级黄色大片| 欧美bdsm另类| 精品久久久久久电影网| 女人十人毛片免费观看3o分钟| 丝袜脚勾引网站| 亚洲人成网站在线观看播放| 91久久精品电影网| 国内揄拍国产精品人妻在线| 亚洲国产欧美在线一区| 一个人观看的视频www高清免费观看| 少妇 在线观看| 蜜桃久久精品国产亚洲av| 麻豆成人午夜福利视频| 五月天丁香电影| 精品99又大又爽又粗少妇毛片| 伦精品一区二区三区| 国产真实伦视频高清在线观看| 国产欧美另类精品又又久久亚洲欧美| 深爱激情五月婷婷| 天天躁日日操中文字幕| 免费黄色在线免费观看| 性色av一级| videos熟女内射| 成人免费观看视频高清| 亚洲av国产av综合av卡| 精品国产乱码久久久久久小说| 午夜激情久久久久久久| 精品久久久精品久久久| 国产v大片淫在线免费观看| 久久久欧美国产精品| 成人毛片60女人毛片免费| 99久久精品热视频| 亚洲精品久久午夜乱码| 韩国av在线不卡| 另类亚洲欧美激情| 午夜福利高清视频| 91aial.com中文字幕在线观看| 97在线视频观看| 男插女下体视频免费在线播放| 亚洲国产精品国产精品| 男女那种视频在线观看| 国产精品一区二区在线观看99| 精品人妻熟女av久视频| 尤物成人国产欧美一区二区三区| 欧美日韩综合久久久久久| 国产精品人妻久久久影院| 国产片特级美女逼逼视频| 99热这里只有是精品50| 成人黄色视频免费在线看| 中文精品一卡2卡3卡4更新| 久久影院123| 欧美成人精品欧美一级黄| 国产成人freesex在线| 22中文网久久字幕| 国产成年人精品一区二区| 日韩欧美精品免费久久| 久久久a久久爽久久v久久| 肉色欧美久久久久久久蜜桃 | av播播在线观看一区| 草草在线视频免费看| 久久97久久精品| 国产69精品久久久久777片| 欧美变态另类bdsm刘玥| 午夜日本视频在线| 亚洲av电影在线观看一区二区三区 | 干丝袜人妻中文字幕| 特级一级黄色大片| eeuss影院久久| 亚洲精品影视一区二区三区av| av线在线观看网站| 99视频精品全部免费 在线| 亚州av有码| 在线免费十八禁| 国产在线男女| av在线观看视频网站免费| 嫩草影院新地址| 熟妇人妻不卡中文字幕| av.在线天堂| 国产精品国产三级国产av玫瑰| 2021少妇久久久久久久久久久| 欧美成人午夜免费资源| 日韩成人伦理影院| 欧美最新免费一区二区三区| freevideosex欧美| 一区二区av电影网| 国国产精品蜜臀av免费| 舔av片在线| 在线 av 中文字幕| 丝瓜视频免费看黄片| 九草在线视频观看| 午夜福利在线在线| 最近中文字幕高清免费大全6| 欧美高清性xxxxhd video| 丰满乱子伦码专区| 黄色视频在线播放观看不卡| 久久久久久久亚洲中文字幕| 亚洲欧美日韩东京热| 26uuu在线亚洲综合色| 天天躁日日操中文字幕| 久久久久网色| 赤兔流量卡办理| h日本视频在线播放| 亚洲国产欧美在线一区| 一本一本综合久久| 秋霞伦理黄片| 国产女主播在线喷水免费视频网站| 少妇被粗大猛烈的视频| 午夜亚洲福利在线播放| 51国产日韩欧美| 波多野结衣巨乳人妻| 日本-黄色视频高清免费观看| 天堂网av新在线| 人人妻人人看人人澡| 人妻夜夜爽99麻豆av| av在线天堂中文字幕| 欧美 日韩 精品 国产| 国产探花极品一区二区| 国产爱豆传媒在线观看| 99九九线精品视频在线观看视频| 乱系列少妇在线播放| 亚洲,一卡二卡三卡| 中文字幕亚洲精品专区| 人妻 亚洲 视频| 新久久久久国产一级毛片| 国产免费福利视频在线观看| 综合色av麻豆| 日本三级黄在线观看| 欧美三级亚洲精品| 夫妻午夜视频| 亚洲人成网站在线播| 国产成人精品久久久久久| 国产一区二区亚洲精品在线观看| 久热久热在线精品观看| 亚洲第一区二区三区不卡| 听说在线观看完整版免费高清| 精品亚洲乱码少妇综合久久| 欧美精品人与动牲交sv欧美| 天美传媒精品一区二区| 最近最新中文字幕大全电影3| 亚洲精品日韩在线中文字幕| 卡戴珊不雅视频在线播放| 免费人成在线观看视频色| 久久久久精品久久久久真实原创| 简卡轻食公司| 99精国产麻豆久久婷婷| 中国国产av一级| 国产欧美日韩一区二区三区在线 | 人妻 亚洲 视频| 国产午夜精品久久久久久一区二区三区| 免费不卡的大黄色大毛片视频在线观看| 热99国产精品久久久久久7| 大又大粗又爽又黄少妇毛片口| 亚洲精品aⅴ在线观看| 黑人高潮一二区| 99久久精品国产国产毛片| 国产高清三级在线| 制服丝袜香蕉在线| 我的女老师完整版在线观看| 性插视频无遮挡在线免费观看| 久热久热在线精品观看| 免费av观看视频| 国产成人精品久久久久久| 日韩av免费高清视频| 最近的中文字幕免费完整| .国产精品久久| 免费高清在线观看视频在线观看| 久久综合国产亚洲精品| 国产一级毛片在线| 久久久久久久久久人人人人人人| 亚洲欧美成人综合另类久久久| 嫩草影院新地址| 在线亚洲精品国产二区图片欧美 | 久久久久久久大尺度免费视频| 久久这里有精品视频免费| 22中文网久久字幕| 激情 狠狠 欧美| 亚洲精品乱码久久久v下载方式| 亚洲精品aⅴ在线观看| 青春草亚洲视频在线观看| 下体分泌物呈黄色| 人人妻人人看人人澡| 亚洲精品视频女| 啦啦啦中文免费视频观看日本| 欧美日韩亚洲高清精品| 久久女婷五月综合色啪小说 | av在线播放精品| 真实男女啪啪啪动态图| 久久久久久久久久久免费av| 我要看日韩黄色一级片| 亚洲精品乱码久久久v下载方式| 纵有疾风起免费观看全集完整版| 听说在线观看完整版免费高清| 蜜桃久久精品国产亚洲av| 在线观看av片永久免费下载| 免费看光身美女| 免费不卡的大黄色大毛片视频在线观看| 日本爱情动作片www.在线观看| 亚洲国产精品国产精品| 亚洲精品视频女| 精品久久久久久久末码| 国产 一区 欧美 日韩| 国产精品一区二区性色av| 国产亚洲午夜精品一区二区久久 | av在线蜜桃| 亚洲国产欧美人成| 又爽又黄无遮挡网站| 午夜福利在线观看免费完整高清在| 91久久精品国产一区二区三区| 狠狠精品人妻久久久久久综合| 日韩av不卡免费在线播放| 国产亚洲最大av| 久久99热6这里只有精品| 特级一级黄色大片| 亚洲精品乱久久久久久| 国产中年淑女户外野战色| 人妻制服诱惑在线中文字幕| 自拍欧美九色日韩亚洲蝌蚪91 | 亚洲真实伦在线观看| 男人和女人高潮做爰伦理| 日本熟妇午夜| 大片电影免费在线观看免费| 熟妇人妻不卡中文字幕| 熟女av电影| av在线蜜桃| 能在线免费看毛片的网站| 国产精品.久久久| 成人漫画全彩无遮挡| 男女那种视频在线观看| 在线观看一区二区三区| 久久人人爽人人爽人人片va| 国产淫片久久久久久久久| 黄色欧美视频在线观看| 成人国产麻豆网| 日韩一区二区三区影片| 十八禁网站网址无遮挡 | 亚洲国产色片| 成年女人在线观看亚洲视频 | 免费看av在线观看网站| 午夜福利网站1000一区二区三区| 一二三四中文在线观看免费高清| 免费av观看视频| 国产成年人精品一区二区| 精品人妻一区二区三区麻豆| 日本与韩国留学比较| 国模一区二区三区四区视频| 97超碰精品成人国产| 亚洲国产欧美人成| 国产精品麻豆人妻色哟哟久久| 一区二区av电影网| 国产精品国产三级国产专区5o| av在线播放精品| 美女被艹到高潮喷水动态| 亚洲av男天堂| 色视频www国产| 人妻夜夜爽99麻豆av| 午夜福利视频1000在线观看| 丝瓜视频免费看黄片| 日本色播在线视频| 国产 精品1| 精品熟女少妇av免费看| 久久精品久久精品一区二区三区| 亚洲欧美日韩东京热| 观看美女的网站| 亚洲精品一二三| 激情 狠狠 欧美| 国产精品国产av在线观看| 国产亚洲午夜精品一区二区久久 | 午夜福利高清视频| 免费观看无遮挡的男女| 国产欧美日韩精品一区二区| 菩萨蛮人人尽说江南好唐韦庄| 亚洲成人久久爱视频| 国产黄片美女视频| 欧美少妇被猛烈插入视频| 中文在线观看免费www的网站| 亚洲天堂av无毛| 国产熟女欧美一区二区| 国产熟女欧美一区二区| 丰满少妇做爰视频| 免费看光身美女| 国产精品三级大全| 成人一区二区视频在线观看| 自拍欧美九色日韩亚洲蝌蚪91 | 国产亚洲5aaaaa淫片| 精品久久久精品久久久| 成人亚洲精品一区在线观看 | 精品一区在线观看国产| 国产在线一区二区三区精| 菩萨蛮人人尽说江南好唐韦庄| 亚洲久久久久久中文字幕| 自拍偷自拍亚洲精品老妇| 中文乱码字字幕精品一区二区三区| 全区人妻精品视频| 嫩草影院入口| 九草在线视频观看| 搡老乐熟女国产| 欧美日本视频| 激情五月婷婷亚洲| 日韩在线高清观看一区二区三区| 日韩强制内射视频| 在线观看人妻少妇| 我的女老师完整版在线观看| 老女人水多毛片| 亚洲欧美成人综合另类久久久| 内地一区二区视频在线| 3wmmmm亚洲av在线观看| 免费看日本二区| 国产一区二区三区av在线| 狂野欧美白嫩少妇大欣赏| 永久网站在线| 大片免费播放器 马上看| 内射极品少妇av片p| 久久久精品欧美日韩精品| 日韩欧美一区视频在线观看 | 亚洲av二区三区四区| 一二三四中文在线观看免费高清| 亚洲成色77777| 卡戴珊不雅视频在线播放| 人妻制服诱惑在线中文字幕| 精品久久久久久久末码| 熟女电影av网| 丰满少妇做爰视频| 免费观看av网站的网址| 日韩免费高清中文字幕av| 国产高清不卡午夜福利| 18禁在线无遮挡免费观看视频| 国产老妇女一区| 亚洲va在线va天堂va国产| av线在线观看网站| 久久久久久伊人网av| 观看免费一级毛片| 日韩三级伦理在线观看| 中文字幕久久专区| 九色成人免费人妻av| 小蜜桃在线观看免费完整版高清| 丰满乱子伦码专区| 久久久欧美国产精品| 欧美性猛交╳xxx乱大交人| 亚洲精华国产精华液的使用体验| 亚洲av福利一区| 午夜免费观看性视频| 国产男女内射视频| 看黄色毛片网站| 国产高清国产精品国产三级 | 午夜免费男女啪啪视频观看| 中国国产av一级| 久久鲁丝午夜福利片| 视频中文字幕在线观看| 十八禁网站网址无遮挡 | 国产乱人视频| 一级黄片播放器| 蜜臀久久99精品久久宅男| 欧美3d第一页| 日韩伦理黄色片| 亚洲天堂国产精品一区在线| 亚洲av福利一区| 99久久精品一区二区三区| 综合色av麻豆| 十八禁网站网址无遮挡 | 插逼视频在线观看| 国产精品嫩草影院av在线观看| 国产高潮美女av| 18禁在线无遮挡免费观看视频| 九九在线视频观看精品| 男插女下体视频免费在线播放| av在线蜜桃| 欧美xxxx性猛交bbbb| 欧美日韩在线观看h| 国产欧美日韩一区二区三区在线 | 日韩欧美精品免费久久| 熟女电影av网| 人妻少妇偷人精品九色| 婷婷色麻豆天堂久久| 啦啦啦啦在线视频资源| 哪个播放器可以免费观看大片| 欧美成人午夜免费资源| 中国三级夫妇交换| 中文字幕久久专区| 国产精品精品国产色婷婷| 精品酒店卫生间| 久久国内精品自在自线图片| 日韩一本色道免费dvd| 久久精品国产亚洲网站| 亚洲天堂av无毛| 亚洲精品一区蜜桃| 亚洲精品乱久久久久久| 久热这里只有精品99| 校园人妻丝袜中文字幕| 三级男女做爰猛烈吃奶摸视频| 久久国内精品自在自线图片| 天天躁夜夜躁狠狠久久av| 美女高潮的动态| eeuss影院久久| 最近最新中文字幕大全电影3| 男女边吃奶边做爰视频| 欧美老熟妇乱子伦牲交| 免费观看无遮挡的男女| 天堂中文最新版在线下载 | 高清视频免费观看一区二区| av卡一久久| 制服丝袜香蕉在线| 韩国av在线不卡| 精品国产三级普通话版| 午夜精品一区二区三区免费看| 国产视频首页在线观看| 精品视频人人做人人爽| 久久精品综合一区二区三区| 久久久欧美国产精品| 日韩av免费高清视频| 亚洲在线观看片| 国产一区亚洲一区在线观看| 人妻 亚洲 视频| 亚洲色图综合在线观看| 亚洲国产高清在线一区二区三| 午夜爱爱视频在线播放| 亚洲自拍偷在线| 内地一区二区视频在线| 国产成人精品婷婷| 大又大粗又爽又黄少妇毛片口| 久久综合国产亚洲精品| 在线a可以看的网站| 黄色视频在线播放观看不卡| 亚洲天堂av无毛| 老女人水多毛片| 搡老乐熟女国产| 亚洲av二区三区四区| 国产 精品1| 美女主播在线视频| 国产 一区 欧美 日韩| 国产亚洲最大av| 亚洲成色77777| 人妻少妇偷人精品九色| 亚洲成人中文字幕在线播放| 亚洲天堂国产精品一区在线| 国产亚洲最大av|