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

    開口管樁高頻振動貫入過程的ALE有限元分析*

    2016-08-19 08:04:53肖勇杰陳福全林良慶
    工程地質(zhì)學(xué)報 2016年3期
    關(guān)鍵詞:管樁徑向阻力

    肖勇杰 陳福全 林良慶

    (福州大學(xué)土木工程學(xué)院 福州 350116)

    ?

    開口管樁高頻振動貫入過程的ALE有限元分析*

    肖勇杰陳福全林良慶

    (福州大學(xué)土木工程學(xué)院福州350116)

    任意拉格朗日-歐拉(ALE)方法吸取了拉格朗日和歐拉法的優(yōu)點,避免了常規(guī)有限元中拉格朗日方法的網(wǎng)格畸變問題,適用于開口管樁高頻振動貫入過程的計算分析。采用ALE有限元方法,建立開口管樁高頻振動貫入過程的數(shù)值模型,對沉樁過程中擠土效應(yīng)、樁側(cè)阻力和土塞效應(yīng)的變化規(guī)律進行了詳細(xì)研究。研究結(jié)果表明:擠土應(yīng)力主要沿徑向傳播,且深層土體受到的擠土應(yīng)力比淺層土體大; 水平擠土位移隨管樁貫入深度的增加而增大,而最大水平擠土位移與管樁貫入深度存在累積效應(yīng); 擠土效應(yīng)的影響范圍約為10倍管徑,因此在施工過程中要給以足夠重視; 樁外側(cè)摩阻力隨貫入深度增加呈近似線性增長,樁內(nèi)側(cè)摩阻力隨貫入深度增加而呈非線性增長,增長速率隨貫入深度增加而逐漸增大; 管內(nèi)土塞處于不完全閉塞狀態(tài),土塞程度由完全非閉塞向部分閉塞過渡。此外,研究了土體模量、樁土界面摩擦系數(shù)、振動頻率和樁徑對土體位移的影響。

    開口管樁高頻振動ALE方法有限元

    0 引 言

    隨著巖土工程技術(shù)的發(fā)展和大型工程建設(shè)的需求,開口管樁在各種建筑基礎(chǔ)中得到越來越廣泛的應(yīng)用(Liu, 2008; Liu et al.,2009; 丁選明等, 2013; 許崧等, 2013)。開口管樁主要沉樁方法包括靜壓沉樁法、錘擊沉樁法和振動沉樁法。當(dāng)開口管樁采用靜壓法或錘擊法貫入時,樁周土體不會產(chǎn)生大變形,常規(guī)有限元能有效模擬開口管樁貫入過程(Liyanapathirana et al.,2000)。然而,開口管樁高頻振動貫入過程是一個非常復(fù)雜的結(jié)構(gòu)與土體相互作用的力學(xué)過程,不僅涉及到高頻振動與擠壓對周圍土層的強烈擾動作用,還涉及到管樁內(nèi)外側(cè)壁與土體間的高頻反復(fù)剪切與滑移,樁周土體產(chǎn)生大變形。采用常規(guī)有限元模擬開口管樁高頻振動貫入過程時,容易出現(xiàn)樁周附近土體大變形、網(wǎng)格高度扭曲和邊界條件改變明顯等數(shù)值計算問題,導(dǎo)致計算結(jié)果不精確,甚至計算不收斂。

    許多學(xué)者針對壓樁大變形問題進行了研究。Van den Berg(1994)較早采用Eulerian有限元法模擬樁貫入土體。Eulerian有限元單元在空間上是固定的,材料在空間流動,單元不會隨著材料的運動而扭曲。但Eulerian單元在處理移動邊界和樁土相互作用的問題上受到較大限制。Feng(1997)利用FLAC有限差分程序建立樁振動貫入的二維模型,分析了土體條件、振動錘參數(shù)和樁體參數(shù)對振動沉樁過程的影響。有限差分單元材料可以采用線性或非線性本構(gòu)模型,當(dāng)材料發(fā)生屈服流動后,網(wǎng)格能夠相應(yīng)發(fā)生變形和移動,從而避免了網(wǎng)格扭曲問題。由于有限差分法無需反復(fù)迭代實現(xiàn)非線性本構(gòu)關(guān)系,其在分析大變形問題時比有限元更具優(yōu)勢。但FLAC有限差分程序的前處理功能較弱,建立復(fù)雜三維模型較為困難。

    近年來,任意拉格朗日-歐拉方法(Arbitrary Lagrangian-Eulerian, ALE)被發(fā)現(xiàn)是處理大變形網(wǎng)格扭曲問題的有效方法。ALE方法最早是由Noh(1964)根據(jù)“混合歐拉-拉格朗日”方法提出的,并用于模擬流體力學(xué)的自由面和邊界條件運動。隨后,該法被分別應(yīng)用于解決二維和三維流體問題(Hirt et al.,1974; Stein et al.,1977)。Belytschko et al.(1978)首次將ALE方法嵌入到有限元程序中,解決了核安全分析中常見的流固耦合問題,主要特點是流體采用ALE方法分析,固體采用拉格朗日方法分析。Hughes et al.(1981)結(jié)合運動學(xué)理論,提出了適用于考慮黏性不可壓縮流體的有限元程序。隨著ALE理論不斷發(fā)展,ALE方法逐漸應(yīng)用于固體力學(xué)的大變形問題。Liu et al.(1986)最早將ALE方法應(yīng)用于固體力學(xué)問題求解,建立了與路徑相關(guān)材料ALE描述的一般格式和非線性ALE有限元分析的顯式計算過程。近來,ALE方法被廣泛應(yīng)用于解決巖土工程問題。Hu et al.(1998)提出了ALE方法的土體大變形數(shù)值程序,該程序結(jié)合了傳統(tǒng)小應(yīng)變有限元方法與自動網(wǎng)格劃分技術(shù),用于解決平面應(yīng)變及軸對稱的圓孔擴張和地基承載力等問題。Walker et al.(2006)采用ALE有限元方法,對靜力觸探試驗(CPT)進行數(shù)值模擬,并與前人研究結(jié)果比較,驗證了數(shù)值模擬結(jié)果的可靠性。結(jié)果表明,ALE方法能有效避免常規(guī)有限元分析中網(wǎng)格畸變及高度扭曲問題,真實地反映圓錐貫入時土體的變形情況,可用于分析巖土貫入問題。Nazem et al.(2006, 2008)基于ALE有限元原理,提出了一種新的網(wǎng)格優(yōu)化技術(shù),并成功地應(yīng)用于求解土體固結(jié)、滲透問題。Henke et al.(2008)利用ALE有限元方法,對開口管樁振動貫入過程中土塞效應(yīng)進行研究。研究發(fā)現(xiàn),管樁振動貫入過程中管內(nèi)土塞處于不完全閉塞狀態(tài),與試驗結(jié)果相一致。Sheng et al.(2009)分別采用ALE法和拉格朗日(Updated Lagrangian, UL)法分析樁的貫入過程。與UL法相比,ALE法對于處理樁貫入過程土體單元網(wǎng)格畸變問題是更有效的,計算成本也相對較低。Ekanayake et al.(2013)采用ALE有限元方法,研究了振動打樁過程中地面振動規(guī)律。

    由于ALE方法吸取了拉格朗日和歐拉法的優(yōu)點,避免了常規(guī)有限元中拉格朗日方法的網(wǎng)格畸變問題,能更好地模擬開口管樁高頻振動貫入過程。因此,采用ALE有限元方法精確分析開口管樁高頻振動貫入過程,有助于更加深入、全面認(rèn)識開口管樁高頻振動貫入的沉樁機理及樁周土體性狀變化過程。本文首先詳細(xì)介紹了ALE方法原理,研究了其求解步驟,及其在高頻打樁全過程分析中的建模過程和方法; 然后采用ALE有限元方法詳細(xì)地研究了開口管樁高頻振動貫入全過程,并將ALE有限元方法的計算結(jié)果與淺層應(yīng)變路徑法(SSPM)、更新的拉格朗日有限元方法(UL)的計算結(jié)果進行比較,探究ALE有限元技術(shù)的優(yōu)點及其分析開口管樁高頻振動貫入過程的適用性; 接著利用ALE有限元方法,深入研究了開口管樁振動沉樁擠土效應(yīng)、樁側(cè)阻力和土塞效應(yīng); 最后分析了土體模量、樁土界面摩擦系數(shù)、振動頻率和樁徑對土體位移的影響。

    1 任意拉格朗日-歐拉法的原理

    任意拉格朗日-歐拉(ALE)法是基于算子分裂技術(shù)建立的(Sheng et al.,2009),并通過更新的拉格朗日步建立基本平衡方程。

    1.1更新的拉格朗日解

    根據(jù)虛功原理,如果虛位移δu滿足位移邊界條件,可以建立平衡方程:

    (1)

    式中,δε為虛位移引起的應(yīng)變張量的變量; σ為Cauchy應(yīng)力張量; b為體力矢量;q是施加在體積Vα邊界Sα上的分布力; tN和tT分別為接觸面Sc上的法向力和切向力; δgN和δgT分別為虛法向間隙和虛切向間隙。

    對于非線性問題,假設(shè)分析從0時刻開始,t時刻所有的狀態(tài)變量都滿足增量形式的平衡方程:

    (2)

    式(2)中的組態(tài)(Vt+Δt和St+Δt)是未知的,需轉(zhuǎn)換為一個已知的組態(tài),通常是采用第二Piola-Kirchhoff應(yīng)力σ和Green-Lagrangian應(yīng)變ε來消除剛體運動對Cauchy應(yīng)力張量的影響。對于大滑移接觸問題,接觸約束是采用當(dāng)前的組態(tài)來描述,更新的拉格朗日方程更適合該類接觸問題。然而,在巖土材料本構(gòu)關(guān)系中,第2 Piola-Kirchhoff應(yīng)力σ和Green-Lagrangian應(yīng)變ε的應(yīng)力-應(yīng)變關(guān)系難以實現(xiàn),則將Cauchy應(yīng)力率分解為一個與應(yīng)變有關(guān)的構(gòu)架獨立的應(yīng)力率(Jaumann應(yīng)力率)和一個與剛體運動有關(guān)的應(yīng)力率,即:

    (3)

    (4)

    式中,Cijkl為Cauchy應(yīng)力和線性應(yīng)變組成的應(yīng)力應(yīng)變張量;Ω為旋轉(zhuǎn)張量。將式(3)代入式(2),則基于Jaumann應(yīng)力率的UL法平衡方程為:

    (5)

    式中,δ(dη)為增量Green-Lagrangian應(yīng)變張量非線性部分的變量; 公式左邊第一項為材料非線性,第2和第3項為幾何非線性。

    1.2ALE求解步驟

    UL步結(jié)束后,因網(wǎng)格與材料一起移動而產(chǎn)生畸變。若將網(wǎng)格和材料的位移分離,使網(wǎng)格獨立于材料移動,可避免這種畸變的發(fā)生。在這種假設(shè)的前提下,可在每一個UL步后加入歐拉步,從而構(gòu)成ALE解。

    在歐拉步中,生成了一個新網(wǎng)格的變形區(qū)域,并根據(jù)材料導(dǎo)數(shù)與網(wǎng)格導(dǎo)數(shù)的關(guān)系,將所有舊網(wǎng)格上的變量轉(zhuǎn)移到新網(wǎng)格上,即:

    (6)

    由于網(wǎng)格位移是材料位移的耦合,故歐拉步中生成的新網(wǎng)格可以是任意的,但需滿足兩個基本條件:(1)新的網(wǎng)格要符合區(qū)域和材料的變形邊界; (2)新舊網(wǎng)格的拓?fù)浣Y(jié)構(gòu)和連通性要保持一致。

    歐拉步的第一個目標(biāo)是將經(jīng)UL步后得到的材料和區(qū)域邊界進行理想分割,把重新分割后所得節(jié)點與初始網(wǎng)格節(jié)點對比,獲得節(jié)點位移,以此對初始網(wǎng)格建立一個彈性分析,來獲得內(nèi)部節(jié)點位移,實現(xiàn)理想分割網(wǎng)格的目的。由于邊界節(jié)點位移是指定的,則分析所用彈性參數(shù)可以是任意的,這種情況下仍可得到與舊網(wǎng)格有相同連通性和拓?fù)浣Y(jié)構(gòu)的新網(wǎng)格,再將狀態(tài)變量從舊網(wǎng)格映射到新網(wǎng)格,則ALE方法完成。綜上所述,ALE方法的主要步驟為:(1)執(zhí)行UL步獲得材料位移; (2)重新分配區(qū)域和材料邊界的節(jié)點; (3)執(zhí)行彈性分析實現(xiàn)網(wǎng)格理想分割; (4)將歷史變量從舊網(wǎng)格映射到新網(wǎng)格。

    2 有限元分析模型

    2.1幾何模型

    應(yīng)用上述ALE計算原理與求解方法,采用Abaqus/Explicit建立開口管樁高頻振動貫入黏土的數(shù)值模型(圖1)。鋼管樁長L=10m,外徑D=80cm,壁厚2cm。管樁與近場域土體采用有限元模擬,有限元土體采用8節(jié)點實體減縮積分單元C3D8R。為了消除振動波反彈的影響,遠(yuǎn)場域土體采用無限元模擬,無限元采用8節(jié)點實體無限單元CIN3D8。打入的管樁采用離散剛體單元R3D4,不考慮管樁變形。有限元分析模型網(wǎng)格(圖2)。

    圖1 有限元模型Fig. 1 Finite element model

    圖2 有限元網(wǎng)格及“拉鏈模型”分析技術(shù)示意圖Fig. 2 Schemes of FEM meshes and zipper-type technique

    由于開口管樁高頻振動貫入過程中,樁周土體被擠壓,且管樁內(nèi)外側(cè)壁與土體間的高頻反復(fù)剪切與滑移,導(dǎo)致樁周土體大變形和網(wǎng)格高度畸變。為解決開口管樁高頻振動貫入過程樁周土體大變形和網(wǎng)格畸變問題,樁周附近1m范圍的土體單元采用ALE網(wǎng)格劃分。

    2.2“拉鏈模型”有限元分析技術(shù)

    近年來,有限元方法廣泛被應(yīng)用于模擬樁的打入過程。Mabsout et al.(1994)最先應(yīng)用軸對稱有限元模型模擬樁打入硬黏土中。但在該模型中,樁預(yù)先打入18m,僅是剩下的錘擊數(shù)目被模擬。隨后,Mabsout et al.將一種“拉鏈模型”有限元分析技術(shù)應(yīng)用于樁的貫入過程,在貫入軸向預(yù)留直徑約為樁徑1%的孔洞,預(yù)留孔洞在貫入過程中擴展,以便建立樁和土的接觸。這種“拉鏈模型”分析技術(shù)也適用于管樁的貫入,即沿管樁壁建一個壁厚t=1mm的剛體管(圖2)。管樁貫入前,管樁與周圍的土設(shè)為無摩擦; 管樁貫入過程中,管樁沿剛體管滑動,土與剛體管分開。采用這種方法,即可建立貫入管樁與周圍土之間的接觸。

    2.3樁土接觸面處理

    管樁壁與土體的接觸采用基于主從原理的Coulomb摩擦接觸模型。由于樁的剛度比土體大得多,則將樁表面設(shè)為主控表面,土表面設(shè)為從屬面。當(dāng)樁土接觸面受拉時,樁與土之間形成裂隙,可以模擬樁與土的脫開過程; 當(dāng)樁土接觸面受壓時,樁土之間可以傳遞剪力,極限剪應(yīng)力為τcrit=μp(μ為接觸面摩擦系數(shù); p為接觸面法向壓力)。

    2.4荷載及邊界條件

    管樁振動貫入的激振力采用力控制模擬,其中靜態(tài)荷載F0=85.5kN,振動頻率f=25Hz,動力荷載幅值Fc=1250kN,激振力Fd=F0+Fcsin(2πft)。

    在有限元區(qū)域的左側(cè)施加水平位移約束,底部施加水平和垂直方向的位移約束。在管樁頂部施加水平向和旋轉(zhuǎn)向的位移約束。無限元區(qū)域無需設(shè)置邊界條件。Abaqus在有限元與無限元的交界面上引入了法向阻尼和切向阻尼,可將反彈回有限元區(qū)域的振動波過濾掉,能夠模擬近場無反射波的情況。

    2.5計算參數(shù)

    土體有限元部分采用Mohr-Coulomb彈塑性模型。由于Abaqus/Explicit采用總應(yīng)力分析方法,土體參數(shù)選用不排水指標(biāo),具體參數(shù)(表1)。

    表1 土體參數(shù)Table1 Parameters of soils

    材料天然重度γ/kN·m-3模量E/MPa泊松比ν黏聚力c/kP內(nèi)摩擦角φ/(°)黏土18100.49280

    對于管樁壁與土的接觸面摩擦角δ的取值,Potyondy(1961)和Acer et al.(1982)的研究表明,對于黏土,取δ/φ′=0.6~0.7是比較合適的。黏土有效內(nèi)摩擦角范圍φ′=20°~32°(李廣信, 2004; 常士驃等, 2007)則δ=12°~22°,即接觸面摩擦系數(shù)μ=tanδ=0.2~0.4,取接觸面摩擦系數(shù)μ=0.25。

    3 數(shù)值模擬結(jié)果及分析

    3.1模型驗證

    3.1.1ALE與SSPM方法對比

    Xu et al.(2006)根據(jù)淺層應(yīng)變路徑法(SSPM),建立了管樁振動貫入過程中地表隆起量的解析式(式7~式8),并將現(xiàn)場試驗結(jié)果與SSPM解析解相比較,證實SSPM法能準(zhǔn)確預(yù)測管樁振動貫入引起的地表隆起量。

    (7)

    (8)

    式中,δz為地表隆起量; Req為管樁等效半徑; r為距管樁軸線的徑向距離; Ls為管樁貫入深度; R為管樁外徑; Ri為管樁內(nèi)徑。

    建立歸一化地表隆起量(δzLs/R2)與歸一化徑向距離(r/Ls)的關(guān)系(圖3)。從圖3中可看出:(1)本文地表隆起量計算值比Xu et al.(2006)研究結(jié)果稍小,因為本文采用的管樁壁厚為2cm, Xu et al.(2006)采用的管樁壁厚為13.8cm,壁厚小的管樁對周圍土體擠土效應(yīng)較小; (2)數(shù)值模擬發(fā)現(xiàn),由于管樁的拖拽效應(yīng),導(dǎo)致管壁周圍出現(xiàn)沉陷區(qū),淺層土體與管壁之間出現(xiàn)脫離,而SSPM法不能反映這一現(xiàn)象; (3)本文計算值更接近SSPM解析解,離散性較低。試驗場地的軟黏土層上存在2m厚的填土,而本文采用均質(zhì)土層進行模擬,因此數(shù)值模擬結(jié)果與SSPM解析解更為吻合。通過圖3 中所描述的地表隆起量與徑向距離的變化關(guān)系,表明本文建立的數(shù)值模型與方法是合理的。

    圖3 地表隆起量隨徑向距離的變化Fig. 3 Variations of ground heave with the radial distance

    3.1.2ALE與UL方法對比

    若改用UL方法,則在上述計算條件下,管樁只能貫入5m; 而采用ALE方法,管樁可以貫入10m。

    圖4分別為采用UL方法和ALE方法開口管樁貫入5m后的網(wǎng)格變形圖。UL方法的網(wǎng)格變形嚴(yán)重,管樁內(nèi)的土體單元穿透管壁,導(dǎo)致計算不能收斂; 采用ALE方法的網(wǎng)格形態(tài)很好,反映出管樁貫入時土體變形的真實情況——管樁內(nèi)側(cè)土受連續(xù)循環(huán)振動,導(dǎo)致土體壓縮變形; 管樁外側(cè)土體向外擠壓。表明ALE方法能夠有效避免常規(guī)有限元中拉格朗日方法的網(wǎng)格畸變問題,保持變形較大處網(wǎng)格形狀仍較規(guī)則,求解結(jié)果與實際變形性狀相符。

    表2為開口管樁貫入5m時,UL與ALE方法求解速度對比。UL方法的增量步和計算時間分別為ALE方法的2倍和3.1倍。UL方法的計算成本高于ALE方法的主要原因是,開口管樁貫入過程中樁周土體產(chǎn)生大變形,UL方法在網(wǎng)格剖分上不能解決大變形問題,只能通過減小增量步長獲得計算收斂,從而導(dǎo)致增量步增加,計算時間也相應(yīng)增加。

    圖4 開口管樁下沉5m后的網(wǎng)格變形圖Fig. 4 Grid deformation map after open-ended pipe pile subsidence of 5ma. UL網(wǎng)格劃分技術(shù);b. ALE網(wǎng)格劃分技術(shù)

    表2 UL與ALE方法計算成本Table2 Computational expense of UL and ALE methods

    增量步計算時間/sUL方法499483174ALE方法246251021

    1. 模型網(wǎng)格單元數(shù): 33480個有限元單元; 810個無限元單元; 2. 計算機主要性能:Intel酷睿四核i5-4570; 內(nèi)存8G

    3.2擠土應(yīng)力

    圖5和圖6 為開口管樁貫入8m后的應(yīng)力等值線云圖。結(jié)果表明,隨管樁貫入,土體徑向應(yīng)力S11和豎向應(yīng)力S33都有所增加。豎向應(yīng)力S33的應(yīng)力泡比徑向應(yīng)力S11小很多,說明振動沉樁過程擠土應(yīng)力主要沿徑向傳播; 所有應(yīng)力最大值均出現(xiàn)在管樁端部,在管樁端部有明顯轉(zhuǎn)折,出現(xiàn)了應(yīng)力集中現(xiàn)象。

    圖5 開口管樁下沉8m后的徑向應(yīng)力S11 等值線云圖Fig. 5 Radial stress S11 nephogram after open-ended pipe pile subsidence of 8m

    圖7為樁周附近土體徑向應(yīng)力隨深度的變化,圖中歸一化深度的0代表地表, 1代表深度為L的土層。研究表明:(1)樁周土受到的最大擠土應(yīng)力出現(xiàn)在樁端到達(dá)相應(yīng)深度的時候; (2)隨著管樁貫入深度的增加,進入管內(nèi)的土體受到的摩擦阻力越來越大,沿徑向擠出的土體越來越多,同時由于上覆土壓力的增加,深層土體受到的擠土應(yīng)力比淺層土體大,與費康(2004)的試驗結(jié)果一致; (3)振動沉樁產(chǎn)生的擠土應(yīng)力豎向影響范圍在樁端下2D左右,在此范圍下的土體徑向應(yīng)力基本為K0狀態(tài)(初始應(yīng)力)下的徑向應(yīng)力。

    圖6 開口管樁下沉8m后的豎向應(yīng)力S33 等值線云圖Fig. 6 Vertical stress S33 nephogram after open-ended pipe pile subsidence of 8m

    圖7 樁周土徑向應(yīng)力隨深度的變化Fig. 7 Variations of radial stresses around the open-ended pipe pile with depth

    圖8 沉樁過程徑向應(yīng)力沿徑向分布(深度為4m)Fig. 8 Distribution of radial stresses with radial distances in open-ended pipe pile penetration(4m in depth)

    圖8為深度4m處土體在沉樁過程中徑向應(yīng)力沿徑向的變化。樁未貫入時,土體僅受自重應(yīng)力作用; 樁貫入2m時,深度4m處土體徑向應(yīng)力略微增加; 樁貫入4m時,樁側(cè)土體徑向應(yīng)力迅速增加; 隨著樁的繼續(xù)貫入,徑向應(yīng)力逐漸減小,這是由于樁貫入過程中樁端處土體徑向應(yīng)力急劇增大導(dǎo)致的。沉樁過程中徑向應(yīng)力沿徑向衰減得很快,在徑向距離10D范圍以外,徑向應(yīng)力接近K0狀態(tài)應(yīng)力??芍駝映翗哆^程產(chǎn)生的擠土應(yīng)力徑向影響范圍為10D。

    3.3擠土位移

    圖9為管樁貫入過程中地表豎向位移隨徑向距離變化,表明:(1)管樁貫入過程中,管樁的拖拽效應(yīng)導(dǎo)致管壁周圍出現(xiàn)沉陷區(qū),淺層土體與管壁之間出現(xiàn)脫離,形成深約9mm的孔縫,降低管樁側(cè)摩阻力,與劉俊偉(2012)的試驗結(jié)果相符合; (2)最大地表豎向位移出現(xiàn)在徑向距離1.5D處,此范圍以外豎向位移逐漸減小,影響范圍為10D; (3)地表豎向位移主要發(fā)生于管樁淺層貫入時,由于上覆土壓力的作用,管樁深層貫入時地表隆起量基本保持不變。

    圖9 地表豎向位移隨徑向距離的變化Fig. 9 Variations of vertical displacement of ground surface with the radial distance

    圖10為管樁貫入過程中徑向距離1.5D處土體的水平位移變化規(guī)律。從圖中可看出,水平位移變化規(guī)律為表層、深層土體位移較小,而中部土體位移較大。分析其原因:表層土體水平位移較小是由于地表無約束,管樁貫入引起的側(cè)向擠土位移由水平變?yōu)榈乇砺∑鸬木壒剩?由于上覆土壓力的作用使中部土體水平位移變大; 管樁高頻振動貫入過程中,振動錘施加在管樁頂部的振動能量使土體產(chǎn)生擠土位移,并以振動波的形式向深層傳遞(Khoubani et al.,2014)因為物質(zhì)阻尼和幾何阻尼的存在,振動波能量隨深度的增大而衰減,導(dǎo)致深層土體位移逐漸減小,且影響范圍在管端下2D左右。

    圖10 水平位移隨深度的變化Fig. 10 Variations of horizontal displacement with depth

    從圖10 還可看出,管端處土體水平位移并未達(dá)到最大值,而是管端以上約1D處土體水平位移達(dá)到最大值,且在深層貫入時基本保持不變。隨著管樁繼續(xù)貫入,土體發(fā)生回彈,上部土體位移稍有減小。由于某一位置的最大擠土位移與管樁的貫入深度存在累積效應(yīng),若周邊有需要保護的管線時,不要誤認(rèn)為只要管樁貫入達(dá)到管線埋置深度時不產(chǎn)生破壞就安全了,要重視隨后超過管線埋置深度的貫入過程對其產(chǎn)生的破壞更大。

    3.4樁側(cè)阻力

    圖11為開口管樁內(nèi)外側(cè)摩阻力隨貫入深度的變化,圖中的內(nèi)、外側(cè)摩阻力是通過樁土接觸面上的剪應(yīng)力沿樁長積分得到(劉漢龍等, 2004)。研究表明:(1)樁外側(cè)摩阻力隨貫入深度增加呈近似線性增長; (2)樁內(nèi)側(cè)摩阻力隨貫入深度增加而呈非線性增長,其增長速率隨貫入深度增加而逐漸增大。

    圖11 內(nèi)外側(cè)摩阻力隨貫入深度的變化Fig. 11 Variations of frictional resistance with penetration depth

    圖12為開口管樁內(nèi)外側(cè)摩阻力分擔(dān)比隨貫入深度的變化,可知內(nèi)外側(cè)摩阻力存在異步發(fā)揮的現(xiàn)象。當(dāng)管樁貫入0~4m時,側(cè)摩阻力主要由外側(cè)摩阻力承擔(dān); 隨著貫入深度增加,外側(cè)摩阻力分擔(dān)比例逐漸減小,內(nèi)側(cè)摩阻力分擔(dān)比例逐漸增加,側(cè)摩阻力由內(nèi)外側(cè)摩阻力共同承擔(dān)。內(nèi)外側(cè)摩阻力的發(fā)揮不同,原因是管樁內(nèi)側(cè)土塞和管樁外側(cè)土的不同邊界條件和位移模式造成的。管樁外側(cè)土可以徑向無限延伸,以剪切變形為主; 管樁內(nèi)側(cè)土塞相當(dāng)于一個一維土柱,往往伴隨著壓縮變形。

    圖12 內(nèi)外側(cè)摩阻力分擔(dān)比隨貫入深度的變化Fig. 12 Variations of ratio of frictional resistance with penetration depth

    3.5土塞效應(yīng)

    開口管樁沉樁過程中大量的土體涌入管內(nèi),形成土塞。土塞會產(chǎn)生完全非閉塞、部分閉塞或完全閉塞3種破壞模式。王磊(2011)考慮土塞與管內(nèi)壁的摩阻力及管樁樁端土的承載力,提出了不完全閉塞(完全非閉塞或部分閉塞)和完全閉塞管樁承載力的表達(dá)式(圖13):

    圖13 開口管樁土塞狀態(tài)示意圖(王磊, 2011)Fig. 13 Soil plug states of an open-ended pipe pile(Wang, 2011)a. 完全非閉塞;b. 部分閉塞; c.完全閉塞

    (9)

    (10)

    式中,Qunplugged、Qplugged分別為開口管樁在土塞不完全閉塞和完全閉塞條件下的單樁承載力; Qsi為管樁內(nèi)側(cè)阻力; Qso為管樁外側(cè)阻力; Qw為管樁樁端阻力; Wp為土塞自重; Qp為土塞端阻力。

    從上式可知,當(dāng)管樁土塞達(dá)到完全閉塞狀態(tài)時管樁內(nèi)側(cè)阻力Qsi為:

    (11)

    管樁內(nèi)土塞處于完全閉塞狀態(tài)時,沒有任何土體繼續(xù)進入管內(nèi),土塞自重不變; 對于同1層土,土塞底部地基極限承載力相同,即土塞端阻力相同。從式(11)可推斷管樁處于完全閉塞狀態(tài)時管樁內(nèi)側(cè)阻力是恒定的。因此,可以利用管樁內(nèi)側(cè)阻力的變化來判斷土塞是否處于閉塞狀態(tài)。

    圖14 開口管樁內(nèi)側(cè)徑向應(yīng)力隨深度的變化Fig. 14 Variations of radial stresses around the inner wall of open-ended pipe pile with depth

    圖14為開口管樁內(nèi)側(cè)徑向應(yīng)力隨深度的變化。表明:(1)隨著貫入深度的增加,管樁內(nèi)側(cè)徑向應(yīng)力逐漸增大,導(dǎo)致管樁內(nèi)側(cè)摩阻力也相應(yīng)增大(圖11),管樁在整個貫入過程中管內(nèi)土塞處于不完全閉塞狀態(tài); (2)當(dāng)管樁貫入0~4m時,管樁內(nèi)側(cè)徑向應(yīng)力小于K0狀態(tài)的初始徑向應(yīng)力,此時管樁內(nèi)土塞處于完全非閉塞狀態(tài); (3)當(dāng)管樁貫入6m時,管樁內(nèi)側(cè)徑向應(yīng)力高于K0狀態(tài)的初始徑向應(yīng)力,此時管樁內(nèi)的土塞處于部分閉塞狀態(tài)。管樁內(nèi)土塞由完全非閉塞狀態(tài)向部分閉塞狀態(tài)過渡的原因是,管樁內(nèi)土體在高頻振動荷載初始作用下,土體受到擾動變松散,管樁容易貫入土體; 隨著貫入深度增加,管樁內(nèi)土體在高頻振動荷載反復(fù)作用下振動壓縮密實,伴隨土塞的壓縮,側(cè)向變形受到鋼管的限制,導(dǎo)致管樁徑向應(yīng)力增加,形成阻止土體進入管內(nèi)的內(nèi)摩阻力。

    4 影響土體位移的參數(shù)分析

    4.1土體模量對土體位移的影響

    圖15為不同土體模量時,管樁貫入6m后徑向距離1.5D處土體水平、豎向位移沿深度變化情況。研究表明:(1)對于水平位移,不同土體模量情況下,其變化規(guī)律趨向一致,即表層、深層土體位移較小,而中部土體位移較大,且隨土體模量增大而減小; (2)對于豎向位移,當(dāng)土體模量較小時,淺層土體產(chǎn)生較大的隆起,而深層土體產(chǎn)生較大的下沉,這是由于土體模量較小時,會產(chǎn)生較大壓縮的緣故。

    圖15 不同土體模量時水平及豎向位移隨深度的變化Fig. 15 Variations of horizontal and vertical displacement with depth at different soil elastic modulus

    圖16 不同樁土摩擦系數(shù)時水平及豎向位移隨深度的變化Fig. 16 Variations of horizontal and vertical displacement with depth at different pile-soil interface friction coefficients

    4.2樁土界面摩擦系數(shù)對土體位移的影響

    圖16為不同樁土界面摩擦系數(shù)時,管樁貫入6m后徑向距離1.5D處土體水平、豎向位移沿深度變化情況。從圖中可看出:(1)隨樁土界面摩擦系數(shù)增大,土體水平位移明顯增大; (2)樁土界面摩擦系數(shù)取值越大,淺層土體的隆起量越小,而深層土體的下沉量越大。這是由于樁土界面摩擦對土體產(chǎn)生拖拽作用,會使向下的豎向位移明顯增大,而徑向表現(xiàn)為較大的水平位移。

    4.3振動頻率對土體位移的影響

    圖17為不同振動頻率時,管樁貫入6m后徑向距離1.5D處土體水平、豎向位移沿深度變化情況。研究表明,振動頻率對土體位移的影響存在一門檻值。當(dāng)振動頻率小于某一值時,土體水平位移和深層土豎向位移隨振動頻率增大而減小; 當(dāng)振動頻率達(dá)到一定值之后,振動頻率對土體位移的影響不明顯。Rodger et al.(1980)通過試驗也發(fā)現(xiàn)振動頻率門檻值的存在。分析認(rèn)為,管樁高頻振動產(chǎn)生的應(yīng)力波會降低土顆粒間接觸壓力,導(dǎo)致土顆粒間摩擦力和樁土之間摩擦力都得到減小; 若振動頻率增大到某一值,土顆粒骨架便產(chǎn)生一種“假液化狀態(tài)”,此時,樁土間的摩阻力將顯著降低甚至喪失,從而導(dǎo)致樁土界面摩擦對土體的拖拽效應(yīng)減弱,土體位移也相應(yīng)減小。

    圖17 不同振動頻率時水平及豎向位移隨深度的變化Fig. 17 Variations of horizontal and vertical displacement with depth at different vibration frequencies

    4.4樁徑對土體位移的影響

    圖18為不同樁徑時,管樁貫入6m后徑向距離1.5D處土體水平、豎向位移沿深度變化情況。從圖中可看出:(1)不同樁徑情況下,土體水平位移仍表現(xiàn)為表層、深層位移較小,而中部位移較大,但隨樁徑的增大,沿樁長深度范圍內(nèi)土體水平位移明顯增大; (2)隨樁徑D的增大,淺層土體的隆起量和深層土體的下沉量均有所增大。

    圖18 不同樁徑時水平及豎向位移隨深度的變化Fig. 18 Variations of horizontal and vertical displacement with depth at different pile diameters

    5 結(jié) 論

    (1)ALE方法吸取了拉格朗日和歐拉法的優(yōu)點,避免了常規(guī)有限元中拉格朗日方法的網(wǎng)格畸變問題,適用于開口管樁高頻振動貫入過程的分析。

    (2)開口管樁高頻振動貫入過程中擠土應(yīng)力主要沿徑向傳播,且深層土體受到的擠土應(yīng)力比淺層土體大; 某一位置的最大擠土位移與管樁的貫入深度存在累積效應(yīng)。

    (3)擠土效應(yīng)(包括擠土應(yīng)力和擠土位移)的影響主要集中在徑向距離為10倍樁徑范圍內(nèi)。

    (4)開口管樁高頻振動貫入過程中樁外側(cè)摩阻力隨貫入深度增加呈近似線性增長; 樁內(nèi)側(cè)摩阻力隨貫入深度增加而呈非線性增長,增長速率隨貫入深度增加而逐漸增大。管樁貫入過程中管內(nèi)土塞處于不完全閉塞狀態(tài); 土塞閉塞程度由完全非閉塞狀態(tài)向部分閉塞狀態(tài)過渡。

    (5)當(dāng)土體模量較小時,淺層土體產(chǎn)生較大的隆起,而深層土體產(chǎn)生較大的下沉; 水平位移隨土體模量增大而減小。樁土界面摩擦系數(shù)取值越大,土體水平位移和深層土體豎向位移均增大,而地表隆起量減小。當(dāng)振動頻率小于某一值時,土體水平位移和深層土豎向位移隨振動頻率增大而減小; 當(dāng)振動頻率達(dá)到一定值之后,振動頻率對土體位移的影響不明顯。隨樁徑的增大,淺層土體的隆起量、深層土體的下沉量和土體水平位移均有所增大。

    Acer Y B,Durgunoglu H T,Tumay M T. 1982. Interface properties of sands[J]. Journal of the Geotechnical Engineering Division,108(4): 648~654.

    Belytschko T B,Kennedy J M. 1978. Computer models for subassembly simulation[J]. Nuclear Engineering and Design,49(1-2): 17~38.

    Chang S P,Zhang S M. 2007. Manual of engineering geology[M]. Beijing: China Architecture & Building Press: 156~162.

    Ding X M,Liu H L. 2013. Time-domain analytical solution of the vibration response of a large-diameter pipe pile subjected to transient concentrated load[J]. Chinese Journal of Geotechnical Engineering,35(6): 1010~1017.

    Ekanayake S D,Liyanapathirana D S,Leo C J. 2013. Influence zone around a closed-ended pile during vibratory driving[J]. Soil Dynamics and Earthquake Engineering,53: 26~36.

    Fei K. 2004. Theory and practice of thin-wall pipe pile using cast-in-situ concrete[Doctoral Thesis][D]. Nanjing: Hohai University.

    Feng Z Y. 1997. Bearing capacity and load-deformation characteristics of vibratory driven piles[Doctoral Thesis][D]. West Lafayette: Purdue University.

    Henke S,Grabe J. 2008. Numerical investigation of soil plugging inside open-ended piles with respect to the installation method[J]. Acta Geotechnica,3(3): 215~223.

    Hirt C W,Amsden A A,Cook J L. 1974. An arbitrary Lagrangian-Eulerian computing method for all flow speeds[J]. Journal of Computational Physics,14(3): 227~253.

    Hu Y,Randolph M F. 1998. A practical numerical approach for large deformation problems in soil[J]. International Journal for Numerical and Analytical Methods in Geomechanics,22(5): 327~350.

    Hughes T J R,Liu W K,Zimmermann T K. 1981. Lagrangian-Eulerian finite element formulation for incompressible viscous flows[J]. Computer Methods in Applied Mechanics and Engineering,29(3): 329~349.

    Khoubani A,Ahmadi M M. 2014. Numerical study of ground vibration due to impact pile driving[J]. Proceedings of the ICE-Geotechnical Engineering,167(1): 28~39.

    Li G X. 2004. Advanced soil mechanics[M]. Beijing: Tsinghua University Press: 143~156.

    Liu E. 2008. Prevention of cracking for large diameter concrete pipe piles[C]∥8th International Conference on the Application of Stress-Wave Theory to Piles[C]. Lisbon:[s.n.].

    Liu H L,F(xiàn)ei K,Zhou Y D,et al. 2004. Numerical simulation of inner frictional resistance of cast-in-situ concrete thin-wall pipe pile[J]. Rock and Soil Mechanics,25(S2): 211~216.

    Liu J W. 2012. Experimental and theoretical studies on the construction effects for jacked open-ended concrete pipe piles[Doctoral Thesis][D]. Hangzhou: Zhejiang University.

    Liu W K,Belytschko T,Chang H. 1986. An arbitrary Lagrangian-Eulerian finite element method for path dependent materials[J]. Computer Method in Applied Mechanics and Engineering,58(2): 227~245.

    Liu R, Yan S W, Li Z H. 2009. Soil plug effect prediction and pile driveability analysis for large-diameter steel piles in ocean engineering[J]. China Ocean Engineering,23(1): 107~118.

    Liyanapathirana D S,Deeks A J,Randolph M F. 2000. Numerical modelling of large deformations associated with driving of open-ended piles[J]. International Journal for Numerical and Analytical Methods in Geomechanics,24(14): 1079~1101.

    Mabsout M E,Tassoulas J L. 1994. A finite element model for the simulation of pile driving[J]. International Journal for Numerical Methods in Engineering,37(2): 257~278.

    Nazem M,Sheng D C,Carter J P. 2006. Stress integration and mesh refinement for large deformation in geomechanics[J]. International Journal for Numerical Methods in Engineering,65(7): 1002~1027.

    Nazem M,Sheng D C,Carter J P,et al. 2008. Arbitrary Lagrangian-Eulerian method for large-strain consolidation problems[J]. International Journal for Numerical and Analytical Methods in Geomechanics,32(9): 1023~1050.

    Noh W. 1964. A time-dependent two-space dimensional coupled Eulerian-Lagrangian code[C]∥Methods in Computational Physics. New York:[s.n.].

    Potyondy J G. 1961. Skin friction between various soils and construction materials[J]. Géotechnique,11(4): 339~353.

    Rodger A A,Littlejohn G S. 1980. A study of vibratory driving in granular soils[J]. Géotechnique, 1980, 30(3): 269~293.

    Sheng D C,Nazem M,Carter J P. 2009. Some computational aspects for solving deep penetration problems in geomechanics[J]. Computational Mechanics,44(4): 549~561.

    Stein L R,Gentry R A,Hirt C W. 1977. Computational simulation of transient blast loading on three dimensional structures[J]. Computer Methods in Applied Mechanics and Engineering,11(1): 57~74.

    Van den Berg P. 1994. Analysis of soil penetration[Doctoral Thesis][D]. Delft: Technische Universiteit Delft.

    Walker J,Yu H S. 2006. Adaptive finite element analysis of cone penetration in clay[J]. Acta Geotechnica,1(1): 43~57.

    Wang L. 2011. The sinking and pulling mechanism of the sleeve of large diameter benoto pile driven by high frequency hammer[D]. Fuzhou: Fuzhou University.

    Xu S,Yan C H,Xu B T. 2013. In-situ test and analysis on mechanical characteristics of piles in sands[J]. Journal of Engineering Geology,21(4): 577~582.

    Xu X T,Liu H L,Lehane B M. 2006. Pipe pile installation effects in soft clay[J]. Geotechnical Engineering,159(4): 285~296.

    常士驃,張?zhí)K民. 2007. 工程地質(zhì)手冊[M]. 北京:中國建筑工業(yè)出版社: 156~162.

    丁選明,劉漢龍. 2013. 大直徑管樁在瞬態(tài)集中荷載作用下的振動響應(yīng)時域解析解[J]. 巖土工程學(xué)報,35(6): 1010~1017.

    費康. 2004. 現(xiàn)澆混凝土薄壁管樁的理論與實踐[博士學(xué)位論文][D]. 南京:河海大學(xué).

    李廣信. 2004. 高等土力學(xué)[M]. 北京:清華大學(xué)出版社: 143~156.

    劉漢龍,費康,周云東,等. 2004. 現(xiàn)澆混凝土薄壁管樁內(nèi)摩阻力的數(shù)值分析[J]. 巖土力學(xué),25(增2): 211~216.

    劉俊偉. 2012. 靜壓開口混凝土管樁施工效應(yīng)試驗及理論研究[博士學(xué)位論文][D]. 杭州:浙江大學(xué).

    王磊. 2001. 大直徑灌注樁套管高頻振動沉拔機理[D]. 福州:福州大學(xué).

    許崧,閻長虹,許寶田. 2013. 砂土地層中樁基受力特征試驗分析[J]. 工程地質(zhì)學(xué)報,21(4): 577~582.

    ALE FINITE ELEMENT ANALYSIS ON PENETRATION PROCESS OF OPEN-ENDED PIPE PILES DRIVEN WITH HIGH FREQUENCY HAMMERS

    XIAO YongjieCHEN FuquanLIN Liangqing

    (School of Civil Engineering of Fuzhou University, Fuzhou350116)

    Arbitrary Lagrangian-Eulerian(ALE)methods couple the advantages of Lagrangian and Eulerian methods, and avoid the mesh distortion problem of Lagrangian method of general finite element. The methods can effectively analyze the penetration process of open-ended pipe piles driven by high frequency hammers. Based on ALE finite element methods, the finite element model of full penetration process of open-ended pipe piles driven by high frequency vibratory hammers is built. The paper studies in detail squeezing effect, frictional resistance and soil plugging effect during pile-sinking. The computational results show that the compacting stress mainly spreads along the horizontal direction, and the compacting stress in deep soil layers is larger than compacting stress in shallow soil layers. The horizontal compacting displacements increase with the increase of penetration depth. But the maximum compacting displacement delays penetration depth. The affecting range of squeezing effect is approximately 10 times pile diameter. So it is necessary to put great emphasis on full penetration process of pipe pile. Outside friction resistance of piles increases linearly with the increase of penetration depth. Inside friction resistance of piles increases nonlinear with the increase of penetration depth. The growth rate increases gradually with the increase of penetration depth. Soil plugs of pipes are incompletely plugged conditions. Degree of soil plugs varies from unplugged conditions to partially plugged conditions. Furthermore, the influence of soil elastic modulus, frictions, vibration frequencies and pile diameters on the compacting displacements are investigated.

    Open-ended pipe pile, High frequency vibration, Arbitrary lagrangian-eulerian method, Finite element method

    10.13544/j.cnki.jeg.2016.03.009

    2015-03-11;

    2015-04-29.

    國家自然科學(xué)基金項目(41272299)資助.

    肖勇杰(1988-),男,博士生,研究方向為土力學(xué)與基礎(chǔ)工程、巖土工程數(shù)值分析. Email: xiao_yongjie@126.com

    簡介: 陳福全(1971-),男,博士,教授,研究方向為土力學(xué)與基礎(chǔ)工程、巖土工程數(shù)值分析、邊坡工程. Email: phdchen@fzu.edu.cn

    TU473.1

    A

    猜你喜歡
    管樁徑向阻力
    基于靜壓預(yù)應(yīng)力混凝土管樁的實施有關(guān)思考
    鼻阻力測定在兒童OSA診療中的臨床作用
    靜壓PHC管樁施工技術(shù)質(zhì)量控制
    淺探徑向連接體的圓周運動
    RN上一類Kirchhoff型方程徑向?qū)ΨQ正解的存在性
    零阻力
    英語文摘(2020年10期)2020-11-26 08:12:12
    基于PID+前饋的3MN徑向鍛造機控制系統(tǒng)的研究
    重型機械(2020年3期)2020-08-24 08:31:40
    濱海軟土地基預(yù)應(yīng)力管樁的應(yīng)用與問題處理
    甘肅科技(2020年20期)2020-04-13 00:30:36
    一類無窮下級整函數(shù)的Julia集的徑向分布
    別讓摩擦成為學(xué)習(xí)的阻力
    久久毛片免费看一区二区三区| 熟女av电影| 美女内射精品一级片tv| 亚洲欧洲精品一区二区精品久久久 | 欧美日韩在线观看h| 久久久久久久久久久久大奶| 一级,二级,三级黄色视频| 国产av码专区亚洲av| 国产白丝娇喘喷水9色精品| 精品久久久久久电影网| 街头女战士在线观看网站| 午夜福利网站1000一区二区三区| 大陆偷拍与自拍| 在线看a的网站| 你懂的网址亚洲精品在线观看| 天堂中文最新版在线下载| 人妻夜夜爽99麻豆av| 肉色欧美久久久久久久蜜桃| 国产熟女欧美一区二区| 肉色欧美久久久久久久蜜桃| 久久人妻熟女aⅴ| 国产一区有黄有色的免费视频| 韩国高清视频一区二区三区| 国产精品99久久99久久久不卡 | 一级毛片我不卡| 女性被躁到高潮视频| 在线免费观看不下载黄p国产| 久久精品夜色国产| 美女视频免费永久观看网站| 99热这里只有精品一区| 亚洲国产精品国产精品| 日本欧美视频一区| 国产在线一区二区三区精| 九九爱精品视频在线观看| 蜜桃久久精品国产亚洲av| 久久精品国产a三级三级三级| 在线观看人妻少妇| 男女无遮挡免费网站观看| 精品一品国产午夜福利视频| 国产av一区二区精品久久| 久久毛片免费看一区二区三区| 亚洲无线观看免费| 久久av网站| 人妻制服诱惑在线中文字幕| 久久久国产精品麻豆| 亚洲久久久国产精品| 精品酒店卫生间| 大话2 男鬼变身卡| 这个男人来自地球电影免费观看 | 中国三级夫妇交换| √禁漫天堂资源中文www| 视频中文字幕在线观看| av在线播放精品| 嫩草影院入口| 国产精品久久久久久精品电影小说| 狂野欧美白嫩少妇大欣赏| 51国产日韩欧美| 亚洲成人一二三区av| 亚洲少妇的诱惑av| 久久婷婷青草| 亚洲国产日韩一区二区| 国产精品一区二区三区四区免费观看| 人妻夜夜爽99麻豆av| 最新中文字幕久久久久| 久久久久精品久久久久真实原创| 亚州av有码| 亚洲国产精品一区三区| 国产精品久久久久久av不卡| 婷婷色综合www| 韩国av在线不卡| 中文字幕制服av| 一区二区日韩欧美中文字幕 | 欧美 日韩 精品 国产| 亚洲精品一区蜜桃| 亚洲三级黄色毛片| 麻豆成人av视频| 日本-黄色视频高清免费观看| 亚洲一区二区三区欧美精品| 26uuu在线亚洲综合色| 欧美另类一区| 亚洲成人手机| 在现免费观看毛片| 在线免费观看不下载黄p国产| 久久久久久久大尺度免费视频| 亚洲少妇的诱惑av| 免费看光身美女| 纵有疾风起免费观看全集完整版| 精品国产乱码久久久久久小说| 少妇的逼好多水| 9色porny在线观看| av卡一久久| 中文字幕精品免费在线观看视频 | 丝袜在线中文字幕| 亚洲图色成人| 在线观看免费视频网站a站| 免费看不卡的av| videosex国产| 国产成人精品久久久久久| 又大又黄又爽视频免费| 久久久国产一区二区| 亚洲色图综合在线观看| 成年人午夜在线观看视频| 飞空精品影院首页| 午夜视频国产福利| 色网站视频免费| 国产精品国产三级专区第一集| 成人国产麻豆网| 久久精品熟女亚洲av麻豆精品| 18禁动态无遮挡网站| 久久这里有精品视频免费| 熟女av电影| 高清午夜精品一区二区三区| 午夜av观看不卡| 国产精品三级大全| 9色porny在线观看| 亚洲av免费高清在线观看| 少妇熟女欧美另类| 亚洲av在线观看美女高潮| 日韩精品免费视频一区二区三区 | 极品少妇高潮喷水抽搐| 久热这里只有精品99| 午夜视频国产福利| 麻豆乱淫一区二区| 飞空精品影院首页| 日本欧美视频一区| 欧美激情极品国产一区二区三区 | 国产成人精品在线电影| 久久久久久久久久久久大奶| 中文天堂在线官网| 少妇人妻久久综合中文| 日韩中字成人| 97在线视频观看| 成人无遮挡网站| 99久久精品国产国产毛片| 亚洲精品中文字幕在线视频| 中文字幕制服av| 国产精品国产av在线观看| 99精国产麻豆久久婷婷| 亚洲美女黄色视频免费看| 亚洲国产欧美在线一区| 精品少妇久久久久久888优播| 一本一本综合久久| 女的被弄到高潮叫床怎么办| 国产亚洲精品第一综合不卡 | 成人亚洲精品一区在线观看| a级毛片在线看网站| 美女国产视频在线观看| 黑人欧美特级aaaaaa片| 下体分泌物呈黄色| 精品亚洲成国产av| 欧美+日韩+精品| 精品熟女少妇av免费看| av免费观看日本| 国产精品国产三级国产av玫瑰| 中文字幕最新亚洲高清| 精品国产国语对白av| 久久精品国产亚洲网站| 亚洲国产成人一精品久久久| 考比视频在线观看| 99国产综合亚洲精品| 欧美日韩成人在线一区二区| 最近最新中文字幕免费大全7| 日韩成人av中文字幕在线观看| 精品一区二区三区视频在线| 久久久久久人妻| 久久久久久久久久成人| 纵有疾风起免费观看全集完整版| 人人妻人人澡人人看| 9色porny在线观看| 久久人人爽人人片av| 国产精品麻豆人妻色哟哟久久| 亚洲欧美成人精品一区二区| 欧美xxⅹ黑人| 免费不卡的大黄色大毛片视频在线观看| 黄色怎么调成土黄色| 韩国av在线不卡| 99久久精品一区二区三区| 中文字幕人妻丝袜制服| 久久久久久久久大av| 99九九线精品视频在线观看视频| 五月天丁香电影| 欧美亚洲日本最大视频资源| 成人综合一区亚洲| 黄色配什么色好看| 亚洲熟女精品中文字幕| 又粗又硬又长又爽又黄的视频| 丝袜美足系列| 在线亚洲精品国产二区图片欧美 | 欧美精品国产亚洲| 精品卡一卡二卡四卡免费| 日韩一区二区三区影片| 大又大粗又爽又黄少妇毛片口| 性高湖久久久久久久久免费观看| 中文欧美无线码| 国产高清国产精品国产三级| 一区二区三区乱码不卡18| 中文字幕免费在线视频6| 熟女av电影| 两个人的视频大全免费| 亚洲天堂av无毛| 久久99精品国语久久久| 久久久久国产精品人妻一区二区| 视频中文字幕在线观看| 丰满少妇做爰视频| 国产精品蜜桃在线观看| 国产免费视频播放在线视频| 亚洲av日韩在线播放| 久久久久久久精品精品| 欧美bdsm另类| 日本av免费视频播放| 黑人猛操日本美女一级片| 久久久久网色| 乱码一卡2卡4卡精品| 久久精品人人爽人人爽视色| 丰满乱子伦码专区| 天堂8中文在线网| 麻豆乱淫一区二区| 成人午夜精彩视频在线观看| 欧美日韩综合久久久久久| 亚洲av电影在线观看一区二区三区| √禁漫天堂资源中文www| 在现免费观看毛片| 青春草视频在线免费观看| 亚洲国产精品一区二区三区在线| 日本vs欧美在线观看视频| 18禁在线无遮挡免费观看视频| 国产69精品久久久久777片| 国产精品国产av在线观看| 亚洲精品亚洲一区二区| a级毛片免费高清观看在线播放| 五月天丁香电影| 午夜影院在线不卡| 老司机影院毛片| 国产精品秋霞免费鲁丝片| av国产精品久久久久影院| 亚洲成人一二三区av| 中国美白少妇内射xxxbb| 久久久久久久久久久久大奶| 精品一区二区免费观看| 国产高清有码在线观看视频| 晚上一个人看的免费电影| 麻豆成人av视频| 最新的欧美精品一区二区| av一本久久久久| 国产精品国产三级国产专区5o| 精品人妻在线不人妻| 男女啪啪激烈高潮av片| 少妇的逼好多水| 曰老女人黄片| 国产av一区二区精品久久| 男女啪啪激烈高潮av片| 亚洲精品国产色婷婷电影| 午夜激情福利司机影院| 久久久精品免费免费高清| 寂寞人妻少妇视频99o| 免费久久久久久久精品成人欧美视频 | 美女中出高潮动态图| 亚洲av日韩在线播放| 日韩av免费高清视频| 夜夜骑夜夜射夜夜干| 国产精品麻豆人妻色哟哟久久| 亚洲国产精品一区二区三区在线| 在线观看三级黄色| 亚洲欧美中文字幕日韩二区| 美女中出高潮动态图| 男女边摸边吃奶| 日韩三级伦理在线观看| 一本久久精品| 久热这里只有精品99| 最新的欧美精品一区二区| 成人漫画全彩无遮挡| 飞空精品影院首页| 在线天堂最新版资源| 久久久久久久久大av| 日韩免费高清中文字幕av| 夜夜爽夜夜爽视频| 91久久精品电影网| 校园人妻丝袜中文字幕| 亚洲精品国产av蜜桃| 亚洲人成77777在线视频| 各种免费的搞黄视频| 欧美 亚洲 国产 日韩一| 秋霞伦理黄片| 成人手机av| 美女福利国产在线| 国产亚洲精品第一综合不卡 | 男女边摸边吃奶| 久久ye,这里只有精品| 国产爽快片一区二区三区| 欧美日韩综合久久久久久| 中国三级夫妇交换| 午夜激情av网站| 十八禁网站网址无遮挡| 美女国产视频在线观看| 国产精品99久久99久久久不卡 | 如日韩欧美国产精品一区二区三区 | 综合色丁香网| 亚洲精品乱码久久久v下载方式| 熟女电影av网| 欧美少妇被猛烈插入视频| 久久久精品区二区三区| 色94色欧美一区二区| 国产国拍精品亚洲av在线观看| 国产精品一区www在线观看| 日日摸夜夜添夜夜添av毛片| 欧美精品一区二区免费开放| 国产成人a∨麻豆精品| 午夜福利视频在线观看免费| 在线看a的网站| 一级毛片电影观看| 飞空精品影院首页| 国产精品一区二区在线观看99| 色吧在线观看| 少妇熟女欧美另类| 色94色欧美一区二区| 五月玫瑰六月丁香| 国产av码专区亚洲av| 日韩一本色道免费dvd| 国产成人a∨麻豆精品| 男人添女人高潮全过程视频| 国产熟女欧美一区二区| 成人国产麻豆网| 亚洲av免费高清在线观看| 国产又色又爽无遮挡免| 国产熟女午夜一区二区三区 | 亚洲中文av在线| 2018国产大陆天天弄谢| 王馨瑶露胸无遮挡在线观看| 天天影视国产精品| 免费高清在线观看视频在线观看| 91在线精品国自产拍蜜月| 精品人妻一区二区三区麻豆| 男的添女的下面高潮视频| 男女边吃奶边做爰视频| 国产无遮挡羞羞视频在线观看| 中国美白少妇内射xxxbb| 草草在线视频免费看| 在现免费观看毛片| 欧美激情国产日韩精品一区| 九九久久精品国产亚洲av麻豆| 日韩亚洲欧美综合| 欧美激情 高清一区二区三区| 在线观看免费视频网站a站| 2022亚洲国产成人精品| 精品少妇久久久久久888优播| 乱码一卡2卡4卡精品| 国产黄色免费在线视频| 黄色毛片三级朝国网站| 精品亚洲成国产av| 国产成人精品无人区| 亚洲精品av麻豆狂野| 精品久久久噜噜| 蜜臀久久99精品久久宅男| 热re99久久国产66热| a级片在线免费高清观看视频| 亚洲av成人精品一二三区| 欧美日韩av久久| 国产一区二区三区综合在线观看 | 亚洲内射少妇av| 纯流量卡能插随身wifi吗| 日本免费在线观看一区| 国产亚洲午夜精品一区二区久久| 女的被弄到高潮叫床怎么办| 制服丝袜香蕉在线| 久久午夜综合久久蜜桃| 中文字幕制服av| 22中文网久久字幕| 黄色毛片三级朝国网站| 成年女人在线观看亚洲视频| 欧美日韩成人在线一区二区| 99久久综合免费| 考比视频在线观看| 亚洲怡红院男人天堂| 国产黄色视频一区二区在线观看| 精品人妻熟女毛片av久久网站| 色婷婷av一区二区三区视频| 国产黄色视频一区二区在线观看| 国产精品.久久久| 有码 亚洲区| 亚洲美女搞黄在线观看| 午夜福利网站1000一区二区三区| 18禁在线播放成人免费| 午夜影院在线不卡| 日韩中字成人| 国产精品99久久99久久久不卡 | 亚洲精品亚洲一区二区| 嫩草影院入口| 热99久久久久精品小说推荐| 国产成人免费观看mmmm| 久久精品国产亚洲网站| 精品国产露脸久久av麻豆| 国产成人av激情在线播放 | 亚洲精品自拍成人| 精品久久久精品久久久| 亚洲性久久影院| 国产精品免费大片| 亚洲国产成人一精品久久久| 91午夜精品亚洲一区二区三区| 久久 成人 亚洲| 新久久久久国产一级毛片| 亚洲天堂av无毛| 日韩 亚洲 欧美在线| 久久精品人人爽人人爽视色| 搡老乐熟女国产| 国产成人精品一,二区| 婷婷色综合大香蕉| 午夜激情av网站| 一级毛片黄色毛片免费观看视频| 成人漫画全彩无遮挡| 国产色婷婷99| 最新中文字幕久久久久| 丝袜脚勾引网站| 夜夜看夜夜爽夜夜摸| 成人影院久久| 午夜影院在线不卡| 国产欧美日韩综合在线一区二区| 秋霞在线观看毛片| 精品熟女少妇av免费看| 久久97久久精品| 国产伦理片在线播放av一区| √禁漫天堂资源中文www| 熟妇人妻不卡中文字幕| 亚洲av国产av综合av卡| 婷婷色综合www| 精品少妇黑人巨大在线播放| 欧美少妇被猛烈插入视频| 午夜福利视频在线观看免费| 3wmmmm亚洲av在线观看| 在线免费观看不下载黄p国产| 免费观看在线日韩| 久久99热这里只频精品6学生| 亚洲av成人精品一区久久| 26uuu在线亚洲综合色| 韩国av在线不卡| 美女大奶头黄色视频| 免费观看的影片在线观看| 女人久久www免费人成看片| 午夜免费鲁丝| 日韩欧美一区视频在线观看| 男人操女人黄网站| 亚洲综合色惰| 午夜久久久在线观看| 建设人人有责人人尽责人人享有的| 99精国产麻豆久久婷婷| 精品少妇久久久久久888优播| 午夜免费男女啪啪视频观看| 国产熟女午夜一区二区三区 | av.在线天堂| 亚洲色图综合在线观看| 成人午夜精彩视频在线观看| 亚洲精品乱久久久久久| 欧美 日韩 精品 国产| 母亲3免费完整高清在线观看 | 国产探花极品一区二区| 国产一区二区三区av在线| 日韩 亚洲 欧美在线| 中文天堂在线官网| 一本大道久久a久久精品| 少妇被粗大猛烈的视频| 热re99久久国产66热| 久久久欧美国产精品| 各种免费的搞黄视频| 三上悠亚av全集在线观看| 国产精品一国产av| 简卡轻食公司| 黑丝袜美女国产一区| 黄色怎么调成土黄色| 97超碰精品成人国产| 黑人高潮一二区| 乱人伦中国视频| 免费观看a级毛片全部| 妹子高潮喷水视频| 亚洲精品亚洲一区二区| 成年人午夜在线观看视频| 人体艺术视频欧美日本| 国精品久久久久久国模美| 18+在线观看网站| 啦啦啦在线观看免费高清www| 极品少妇高潮喷水抽搐| 99热国产这里只有精品6| 免费黄色在线免费观看| 欧美日韩一区二区视频在线观看视频在线| 国产成人av激情在线播放 | 免费观看在线日韩| 大片免费播放器 马上看| 免费看光身美女| 18禁在线播放成人免费| 亚洲欧美成人精品一区二区| 多毛熟女@视频| 男女无遮挡免费网站观看| 日韩三级伦理在线观看| 精品人妻一区二区三区麻豆| 看非洲黑人一级黄片| 热re99久久精品国产66热6| 亚洲精品自拍成人| 亚洲人与动物交配视频| 永久网站在线| 美女内射精品一级片tv| www.av在线官网国产| 在线观看美女被高潮喷水网站| 欧美老熟妇乱子伦牲交| av一本久久久久| 一个人免费看片子| 国产伦精品一区二区三区视频9| 在线亚洲精品国产二区图片欧美 | 青春草国产在线视频| 22中文网久久字幕| 免费日韩欧美在线观看| 欧美三级亚洲精品| 2021少妇久久久久久久久久久| 99国产精品免费福利视频| 精品久久国产蜜桃| 午夜影院在线不卡| 成人国产av品久久久| 少妇人妻久久综合中文| 国产深夜福利视频在线观看| 美女主播在线视频| 91精品国产九色| 这个男人来自地球电影免费观看 | av线在线观看网站| 一区二区av电影网| 久久精品国产a三级三级三级| 成人国产麻豆网| 久久久久久人妻| 国产一区二区在线观看日韩| 如何舔出高潮| 日本午夜av视频| 在线观看www视频免费| 国产一区二区三区av在线| 另类精品久久| 国产熟女午夜一区二区三区 | 国产黄色视频一区二区在线观看| 午夜福利网站1000一区二区三区| 老司机影院成人| 满18在线观看网站| 亚洲精品日本国产第一区| 在线看a的网站| a级毛片黄视频| 寂寞人妻少妇视频99o| 亚洲国产精品成人久久小说| 中文天堂在线官网| 欧美xxxx性猛交bbbb| 亚洲精品乱码久久久久久按摩| 国产精品一区www在线观看| 国产永久视频网站| 91精品伊人久久大香线蕉| 伦理电影大哥的女人| av卡一久久| 亚洲经典国产精华液单| 亚洲高清免费不卡视频| 国产av精品麻豆| 免费观看的影片在线观看| 夜夜看夜夜爽夜夜摸| 国产精品久久久久成人av| 亚洲av日韩在线播放| 欧美日韩av久久| 日韩av在线免费看完整版不卡| 国产日韩欧美视频二区| 18在线观看网站| 国产乱来视频区| 街头女战士在线观看网站| 亚洲av综合色区一区| 欧美日韩在线观看h| 超碰97精品在线观看| 日本色播在线视频| 久热久热在线精品观看| 婷婷色综合大香蕉| 久久久久久人妻| 日日撸夜夜添| 新久久久久国产一级毛片| 97超视频在线观看视频| 黄色欧美视频在线观看| 国产精品免费大片| 少妇的逼水好多| 国产日韩一区二区三区精品不卡 | 制服诱惑二区| 成年av动漫网址| 国产精品国产三级国产av玫瑰| 男人添女人高潮全过程视频| 国产亚洲午夜精品一区二区久久| 80岁老熟妇乱子伦牲交| 久久综合国产亚洲精品| 亚洲中文av在线| 人妻人人澡人人爽人人| 永久免费av网站大全| 亚洲精品成人av观看孕妇| 日本av免费视频播放| 自线自在国产av| 亚洲精品色激情综合| 赤兔流量卡办理| 晚上一个人看的免费电影| 成人黄色视频免费在线看| 中文字幕免费在线视频6| 久久久久久伊人网av| 日本欧美视频一区| 欧美人与善性xxx| 天堂8中文在线网| 国产免费现黄频在线看| 国产成人免费无遮挡视频| 国产熟女欧美一区二区| 欧美日韩av久久| 亚洲色图 男人天堂 中文字幕 | 精品一区二区三区视频在线| 男女啪啪激烈高潮av片| 亚洲精品中文字幕在线视频| 边亲边吃奶的免费视频| 日韩一区二区三区影片| 成人国产麻豆网| 美女中出高潮动态图| 在线观看国产h片| 欧美激情极品国产一区二区三区 | 亚洲色图综合在线观看| 熟女电影av网| 狂野欧美白嫩少妇大欣赏| 欧美人与性动交α欧美精品济南到 |