• <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
    多毛熟女@视频| 亚洲一级一片aⅴ在线观看| 精品一区二区三区四区五区乱码 | 高清视频免费观看一区二区| 成人综合一区亚洲| 精品一区二区三卡| 久久精品国产亚洲av涩爱| av线在线观看网站| 精品人妻熟女毛片av久久网站| 美女大奶头黄色视频| 国产精品人妻久久久影院| 在线观看免费高清a一片| 国产精品一二三区在线看| 国产精品久久久av美女十八| 成年美女黄网站色视频大全免费| 日日撸夜夜添| 在线 av 中文字幕| 少妇精品久久久久久久| 亚洲精华国产精华液的使用体验| 视频在线观看一区二区三区| 美女xxoo啪啪120秒动态图| 毛片一级片免费看久久久久| 中文字幕免费在线视频6| 日韩欧美一区视频在线观看| 免费观看在线日韩| 一级毛片 在线播放| 亚洲色图综合在线观看| 建设人人有责人人尽责人人享有的| 久久狼人影院| 九九在线视频观看精品| 99香蕉大伊视频| 午夜免费鲁丝| 成年av动漫网址| 久久鲁丝午夜福利片| 久久99蜜桃精品久久| 久久ye,这里只有精品| 国内精品宾馆在线| 99热网站在线观看| 国产成人av激情在线播放| 9热在线视频观看99| 一级爰片在线观看| 一级毛片黄色毛片免费观看视频| 久久久国产一区二区| 中文乱码字字幕精品一区二区三区| 丝袜喷水一区| 亚洲av欧美aⅴ国产| 2018国产大陆天天弄谢| 99久国产av精品国产电影| 久久久国产欧美日韩av| 久久国产精品大桥未久av| 丝袜人妻中文字幕| 国产激情久久老熟女| www日本在线高清视频| 久久久国产欧美日韩av| 一级毛片黄色毛片免费观看视频| 有码 亚洲区| 成年人午夜在线观看视频| 亚洲伊人色综图| 欧美日本中文国产一区发布| 男人操女人黄网站| 人人妻人人添人人爽欧美一区卜| 妹子高潮喷水视频| 国产伦理片在线播放av一区| av免费观看日本| 久久精品久久精品一区二区三区| 热99国产精品久久久久久7| 久久国产精品男人的天堂亚洲 | 亚洲精品成人av观看孕妇| 午夜免费观看性视频| 最近最新中文字幕大全免费视频 | 久久亚洲国产成人精品v| 一区二区三区精品91| 精品一区在线观看国产| 99热6这里只有精品| 久久久欧美国产精品| 亚洲精品日韩在线中文字幕| 狂野欧美激情性xxxx在线观看| 亚洲欧美一区二区三区国产| 久久久久精品人妻al黑| 国产欧美日韩综合在线一区二区| 一本大道久久a久久精品| 久久久国产欧美日韩av| 大话2 男鬼变身卡| 看非洲黑人一级黄片| 蜜臀久久99精品久久宅男| 久久人人97超碰香蕉20202| 国产精品麻豆人妻色哟哟久久| 精品人妻在线不人妻| 亚洲在久久综合| 制服丝袜香蕉在线| 亚洲一区二区三区欧美精品| 婷婷色综合大香蕉| 国产国拍精品亚洲av在线观看| 91在线精品国自产拍蜜月| 国产精品久久久久成人av| 国产一级毛片在线| 久久国产精品男人的天堂亚洲 | 自拍欧美九色日韩亚洲蝌蚪91| av在线播放精品| 日韩 亚洲 欧美在线| 自线自在国产av| 久久毛片免费看一区二区三区| 欧美xxⅹ黑人| 亚洲欧美成人综合另类久久久| 高清av免费在线| av电影中文网址| a级毛色黄片| 深夜精品福利| 一区在线观看完整版| 久久国产亚洲av麻豆专区| 亚洲成人一二三区av| 日韩,欧美,国产一区二区三区| 波多野结衣一区麻豆| 久久国产亚洲av麻豆专区| 国产亚洲精品第一综合不卡 | 青春草国产在线视频| 日韩,欧美,国产一区二区三区| 草草在线视频免费看| 国产精品久久久久成人av| 涩涩av久久男人的天堂| 黄色一级大片看看| 亚洲精品久久成人aⅴ小说| 国产乱人偷精品视频| av视频免费观看在线观看| 国产日韩欧美在线精品| 久久久久久久久久成人| 最近最新中文字幕大全免费视频 | 一区二区三区精品91| 国产精品国产三级国产av玫瑰| 青春草亚洲视频在线观看| 欧美人与性动交α欧美软件 | 在线观看免费日韩欧美大片| 中文精品一卡2卡3卡4更新| 久久 成人 亚洲| 少妇熟女欧美另类| 国产精品久久久av美女十八| 黄色一级大片看看| 免费人妻精品一区二区三区视频| 亚洲一级一片aⅴ在线观看| 一二三四中文在线观看免费高清| 国产精品国产三级国产专区5o| 久久精品国产亚洲av涩爱| www.色视频.com| 99久国产av精品国产电影| 国产精品国产av在线观看| 国产一区二区三区av在线| 少妇猛男粗大的猛烈进出视频| 一个人免费看片子| 亚洲精品aⅴ在线观看| 亚洲情色 制服丝袜| 纯流量卡能插随身wifi吗| 精品国产露脸久久av麻豆| 亚洲综合色网址| 麻豆乱淫一区二区| 九草在线视频观看| 18禁裸乳无遮挡动漫免费视频| 最近的中文字幕免费完整| 亚洲欧美一区二区三区黑人 | a级毛色黄片| 成人亚洲精品一区在线观看| 国产精品三级大全| 香蕉精品网在线| 国产高清不卡午夜福利| 国产国拍精品亚洲av在线观看| 国产69精品久久久久777片| 亚洲人成77777在线视频| 国产男人的电影天堂91| 欧美3d第一页| 午夜福利在线观看免费完整高清在| √禁漫天堂资源中文www| 亚洲国产欧美日韩在线播放| 22中文网久久字幕| 天堂8中文在线网| 国产在线一区二区三区精| 99视频精品全部免费 在线| 亚洲内射少妇av| 成人二区视频| 国产成人精品在线电影| 国产男女内射视频| 欧美成人精品欧美一级黄| 九九在线视频观看精品| 在线观看国产h片| 韩国精品一区二区三区 | 国产一级毛片在线| 久久人妻熟女aⅴ| 国产日韩一区二区三区精品不卡| 国产高清三级在线| 伊人亚洲综合成人网| 亚洲精品久久午夜乱码| 国产精品久久久av美女十八| 日韩不卡一区二区三区视频在线| a级毛色黄片| 国产av精品麻豆| 男女国产视频网站| 自线自在国产av| 精品人妻偷拍中文字幕| 黄色毛片三级朝国网站| 精品国产一区二区三区四区第35| 国产极品粉嫩免费观看在线| 人妻系列 视频| 国产av精品麻豆| 久久午夜综合久久蜜桃| 91成人精品电影| 久久99精品国语久久久| videossex国产| 亚洲精华国产精华液的使用体验| av又黄又爽大尺度在线免费看| 久久99热6这里只有精品| 久久久久久久久久人人人人人人| 男女国产视频网站| 丝袜在线中文字幕| 捣出白浆h1v1| 免费黄频网站在线观看国产| 亚洲第一av免费看| 国产一区二区激情短视频 | 久久精品人人爽人人爽视色| 国产男人的电影天堂91| 欧美成人精品欧美一级黄| 日本黄大片高清| 久久精品久久精品一区二区三区| 国产又色又爽无遮挡免| 亚洲av男天堂| 中文字幕亚洲精品专区| 国产女主播在线喷水免费视频网站| 人人妻人人添人人爽欧美一区卜| 国产精品国产av在线观看| 国产成人91sexporn| 亚洲精品成人av观看孕妇| 99精国产麻豆久久婷婷| 人人妻人人添人人爽欧美一区卜| 少妇被粗大的猛进出69影院 | 国产亚洲一区二区精品| 肉色欧美久久久久久久蜜桃| 高清不卡的av网站| 亚洲精品一区蜜桃| 欧美日韩成人在线一区二区| a级毛片黄视频| 国产成人精品福利久久| 男女高潮啪啪啪动态图| 国产亚洲av片在线观看秒播厂| 欧美亚洲 丝袜 人妻 在线| 日韩免费高清中文字幕av| 99re6热这里在线精品视频| 爱豆传媒免费全集在线观看| 国产深夜福利视频在线观看| 婷婷色综合大香蕉| 免费黄色在线免费观看| 免费观看av网站的网址| 91精品国产国语对白视频| 一区二区日韩欧美中文字幕 | 黄色 视频免费看| 伊人亚洲综合成人网| 美女内射精品一级片tv| 色5月婷婷丁香| 爱豆传媒免费全集在线观看| 精品人妻一区二区三区麻豆| 少妇人妻 视频| 亚洲五月色婷婷综合| 国产av精品麻豆| 视频中文字幕在线观看| 亚洲丝袜综合中文字幕| 日韩精品有码人妻一区| 97人妻天天添夜夜摸| 亚洲av电影在线进入| 国产精品免费大片| 新久久久久国产一级毛片| 哪个播放器可以免费观看大片| 人妻少妇偷人精品九色| 免费av不卡在线播放| 免费大片18禁| 亚洲欧美成人精品一区二区| www.熟女人妻精品国产 | 国产亚洲最大av| 国产亚洲av片在线观看秒播厂| 69精品国产乱码久久久| 在线免费观看不下载黄p国产| 人体艺术视频欧美日本| av在线播放精品| 国产精品不卡视频一区二区| 国产男人的电影天堂91| 精品人妻一区二区三区麻豆| 精品一区二区三区四区五区乱码 | 蜜臀久久99精品久久宅男| 色婷婷久久久亚洲欧美| 人人妻人人澡人人爽人人夜夜| 国产视频首页在线观看| 日本欧美视频一区| 精品久久久久久电影网| 飞空精品影院首页| 丝袜人妻中文字幕| 欧美日韩av久久| 人体艺术视频欧美日本| 国产免费视频播放在线视频| 日本黄色日本黄色录像| 国产日韩一区二区三区精品不卡| 亚洲欧美成人精品一区二区| 中文字幕最新亚洲高清| 秋霞伦理黄片| 色哟哟·www| 午夜影院在线不卡| 精品一区二区免费观看| 亚洲精品国产av蜜桃| 黄片无遮挡物在线观看| 水蜜桃什么品种好| 久久99一区二区三区| av电影中文网址| 国产白丝娇喘喷水9色精品| 日韩成人av中文字幕在线观看| 国产又色又爽无遮挡免| 国产色婷婷99| 亚洲欧美成人综合另类久久久| 免费看av在线观看网站| 蜜臀久久99精品久久宅男| 国产精品一区www在线观看| 永久网站在线| 春色校园在线视频观看| 国产欧美亚洲国产| 国产成人av激情在线播放| 18+在线观看网站| 2018国产大陆天天弄谢| 天天影视国产精品| 一本色道久久久久久精品综合| 亚洲美女搞黄在线观看| 看非洲黑人一级黄片| 国产精品久久久久久精品古装| 午夜日本视频在线| 久久精品夜色国产| 边亲边吃奶的免费视频| 在线 av 中文字幕| 亚洲国产精品一区二区三区在线| 在现免费观看毛片| 老司机影院毛片| 日本av手机在线免费观看| 午夜免费男女啪啪视频观看| 国产视频首页在线观看| 国产精品国产av在线观看| 一级a做视频免费观看| 欧美精品高潮呻吟av久久| 日本vs欧美在线观看视频| 高清毛片免费看| 婷婷色综合大香蕉| 女人被躁到高潮嗷嗷叫费观| 国产极品粉嫩免费观看在线| 黑人巨大精品欧美一区二区蜜桃 | 精品一区二区三卡| 自拍欧美九色日韩亚洲蝌蚪91| 一个人免费看片子| 爱豆传媒免费全集在线观看| 日本欧美视频一区| 自拍欧美九色日韩亚洲蝌蚪91| 国产精品久久久av美女十八| 大片电影免费在线观看免费| 男人添女人高潮全过程视频| 成人国产av品久久久| 美女内射精品一级片tv| 90打野战视频偷拍视频| 久久韩国三级中文字幕| av有码第一页| 亚洲成人手机| 香蕉丝袜av| 精品久久久精品久久久| 成人亚洲欧美一区二区av| 精品卡一卡二卡四卡免费| 成人免费观看视频高清| 久久国产亚洲av麻豆专区| 国产一区二区激情短视频 | 十八禁高潮呻吟视频| 久久精品久久久久久噜噜老黄| 女的被弄到高潮叫床怎么办| av福利片在线| 99久久人妻综合| 免费观看a级毛片全部| 青春草国产在线视频| 亚洲第一区二区三区不卡| 校园人妻丝袜中文字幕| 国产淫语在线视频| 制服诱惑二区| 一本大道久久a久久精品| 人妻 亚洲 视频| 街头女战士在线观看网站| 国语对白做爰xxxⅹ性视频网站| 午夜免费男女啪啪视频观看| 十分钟在线观看高清视频www| 另类精品久久| 纵有疾风起免费观看全集完整版| 精品少妇久久久久久888优播| 精品一区二区三卡| 日韩精品有码人妻一区| 亚洲精品第二区| 日本-黄色视频高清免费观看| 脱女人内裤的视频| 成年动漫av网址| 99久久精品国产亚洲精品| 在线观看www视频免费| 亚洲欧美一区二区三区黑人| 热re99久久国产66热| 国产精品国产av在线观看| 国产单亲对白刺激| 精品电影一区二区在线| 一边摸一边抽搐一进一小说 | 欧美日韩乱码在线| 中文字幕av电影在线播放| 亚洲,欧美精品.| 黄片播放在线免费| 国产一区二区三区综合在线观看| 男男h啪啪无遮挡| 黄色视频不卡| av电影中文网址| 最近最新中文字幕大全免费视频| 久久国产精品影院| 久久午夜亚洲精品久久| 久久人人97超碰香蕉20202| aaaaa片日本免费| 亚洲国产精品sss在线观看 | 国产1区2区3区精品| 露出奶头的视频| 欧美乱色亚洲激情| 夜夜爽天天搞| 国产亚洲一区二区精品| 操出白浆在线播放| 亚洲美女黄片视频| 亚洲国产中文字幕在线视频| 午夜福利,免费看| 窝窝影院91人妻| 久久国产精品大桥未久av| 中文字幕精品免费在线观看视频| 欧美在线一区亚洲| 黄网站色视频无遮挡免费观看| 中文欧美无线码| 亚洲精品中文字幕在线视频| 老司机午夜十八禁免费视频| 国产又色又爽无遮挡免费看| 丰满迷人的少妇在线观看| 在线天堂中文资源库| 王馨瑶露胸无遮挡在线观看| 狠狠婷婷综合久久久久久88av| 欧美激情久久久久久爽电影 | 亚洲国产精品一区二区三区在线| 黑人巨大精品欧美一区二区mp4| 午夜精品在线福利| 亚洲av成人不卡在线观看播放网| 中文字幕色久视频| 日韩欧美一区二区三区在线观看 | 国产xxxxx性猛交| 免费久久久久久久精品成人欧美视频| 欧美老熟妇乱子伦牲交| 看片在线看免费视频| 国产男女超爽视频在线观看| 人人妻,人人澡人人爽秒播| 视频区图区小说| 国产一区二区三区视频了| 夜夜爽天天搞| 中文亚洲av片在线观看爽 | 999精品在线视频| 欧美亚洲 丝袜 人妻 在线| 青草久久国产| 一本综合久久免费| 天天添夜夜摸| 好看av亚洲va欧美ⅴa在| 成人亚洲精品一区在线观看| 亚洲成人手机| 老司机影院毛片| 一本综合久久免费| 99re6热这里在线精品视频| 中文字幕最新亚洲高清| 久久精品国产亚洲av香蕉五月 | 三上悠亚av全集在线观看| 99riav亚洲国产免费| 在线观看免费视频日本深夜| 久久久久久久久久久久大奶| 曰老女人黄片| cao死你这个sao货| 国产精品国产高清国产av | 欧美人与性动交α欧美软件| 90打野战视频偷拍视频| 亚洲成人手机| av国产精品久久久久影院| 日韩欧美三级三区| 搡老乐熟女国产| 久久狼人影院| 熟女少妇亚洲综合色aaa.| 巨乳人妻的诱惑在线观看| 18禁裸乳无遮挡动漫免费视频| 十八禁高潮呻吟视频| x7x7x7水蜜桃| 嫩草影视91久久| 日本一区二区免费在线视频| 女人被躁到高潮嗷嗷叫费观| 欧美日韩视频精品一区| 999久久久精品免费观看国产| 这个男人来自地球电影免费观看| av免费在线观看网站| av超薄肉色丝袜交足视频| 精品福利观看| 午夜福利在线观看吧| 欧美日韩精品网址| 母亲3免费完整高清在线观看| 久久精品亚洲熟妇少妇任你| 51午夜福利影视在线观看| 国产免费男女视频| 国产精品久久久久久人妻精品电影| 天堂俺去俺来也www色官网| 精品国产一区二区久久| 久久久精品免费免费高清| 亚洲精品av麻豆狂野| 亚洲av第一区精品v没综合| 精品国产亚洲在线| 首页视频小说图片口味搜索| 在线免费观看的www视频| 久久精品国产综合久久久| 中亚洲国语对白在线视频| 国产精品综合久久久久久久免费 | 交换朋友夫妻互换小说| 成人av一区二区三区在线看| 中文字幕色久视频| 日韩欧美一区二区三区在线观看 | 亚洲 国产 在线| 国产精品 国内视频| av福利片在线| 女人被躁到高潮嗷嗷叫费观| 中文字幕人妻熟女乱码| 亚洲专区中文字幕在线| 国产色视频综合| 亚洲欧美精品综合一区二区三区| 免费女性裸体啪啪无遮挡网站| 亚洲欧美一区二区三区久久| 露出奶头的视频| 91在线观看av| 久久精品人人爽人人爽视色| 亚洲色图综合在线观看| 国产精品九九99| 色综合婷婷激情| 老司机影院毛片| av电影中文网址| 狠狠狠狠99中文字幕| 老熟妇乱子伦视频在线观看| 麻豆av在线久日| 亚洲专区国产一区二区| 亚洲精华国产精华精| 亚洲欧美日韩另类电影网站| 国产成人精品无人区| 757午夜福利合集在线观看| а√天堂www在线а√下载 | 国产一区二区激情短视频| 亚洲av片天天在线观看| 丝袜美腿诱惑在线| 国产亚洲精品一区二区www | a级毛片黄视频| 免费看十八禁软件| 国产精品亚洲一级av第二区| 两性夫妻黄色片| 12—13女人毛片做爰片一| 黄色片一级片一级黄色片| 亚洲美女黄片视频| 国产成人影院久久av| 成熟少妇高潮喷水视频| 在线观看一区二区三区激情| 操出白浆在线播放| 国产精品影院久久| 女人精品久久久久毛片| 啦啦啦免费观看视频1| www.精华液| 午夜免费观看网址| 脱女人内裤的视频| 亚洲人成77777在线视频| 欧美人与性动交α欧美软件| 国产精品久久久久久人妻精品电影| 在线观看免费视频日本深夜| 69精品国产乱码久久久| 人人妻,人人澡人人爽秒播| 中出人妻视频一区二区| 国产蜜桃级精品一区二区三区 | 午夜福利视频在线观看免费| 精品午夜福利视频在线观看一区| 国产av一区二区精品久久| 两个人看的免费小视频| 啦啦啦 在线观看视频| 亚洲一区二区三区不卡视频| 亚洲av熟女| 777久久人妻少妇嫩草av网站| 在线观看免费日韩欧美大片| 18禁观看日本| 国产成人欧美在线观看 | 亚洲一区二区三区不卡视频| 制服诱惑二区| 一区二区日韩欧美中文字幕| 韩国av一区二区三区四区| 如日韩欧美国产精品一区二区三区| 国产淫语在线视频| 久久精品国产综合久久久| 制服诱惑二区| 久久精品aⅴ一区二区三区四区| 看黄色毛片网站| 久久性视频一级片| 精品一品国产午夜福利视频| 亚洲成人手机| 午夜精品在线福利| 一区二区日韩欧美中文字幕| 一本综合久久免费| 亚洲人成电影免费在线| 午夜免费观看网址| www日本在线高清视频| 国产成人一区二区三区免费视频网站| 黄片小视频在线播放| 欧美激情久久久久久爽电影 | 久9热在线精品视频| 国产av一区二区精品久久| 一本大道久久a久久精品| 久99久视频精品免费| 他把我摸到了高潮在线观看| 亚洲va日本ⅴa欧美va伊人久久| 99久久国产精品久久久| 午夜福利,免费看| 精品福利观看| 18禁裸乳无遮挡动漫免费视频|