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

    充填粗糙節(jié)理直剪數(shù)值模擬宏細觀分析

    2018-06-01 05:35:59許萬忠曹日紅
    西南交通大學(xué)學(xué)報 2018年3期
    關(guān)鍵詞:剪切應(yīng)力法向節(jié)理

    許萬忠, 林 杭, 曹日紅

    (1. 昆明理工大學(xué)國土資源工程學(xué)院, 云南 昆明 650093; 2. 中南大學(xué)資源與安全工程學(xué)院, 湖南 長沙 410083; 3. 中國礦業(yè)大學(xué)深部巖土力學(xué)與地下工程國家重點實驗室, 江蘇 徐州 221116)

    經(jīng)歷億萬年的地質(zhì)運動,巖體被各種分界面切割呈破碎狀態(tài),這些軟弱界面包括節(jié)理、裂隙、褶皺等[1].軟弱界面的存在不但使巖體呈各向異性特征,也在很大程度上削弱了巖體的力學(xué)性質(zhì),同時也成為控制各類巖體工程穩(wěn)定性的關(guān)鍵因素.在采礦工程及邊坡工程當中,節(jié)理經(jīng)常受到剪切荷載作用而導(dǎo)致巖體滑移及崩塌的事故.因此,深入研究節(jié)理的的力學(xué)性質(zhì)及破壞行為極為重要.截至目前,已有學(xué)者針對節(jié)理巖體的力學(xué)性質(zhì)及破壞過程進行了相關(guān)研究[2-5]:蔣宇靜等[3]對巖石節(jié)理在不同接觸狀態(tài)下的破壞行為進行了研究;曹平等[6]利用劈裂所得巖石節(jié)理試樣進行直剪實驗并對形貌演化情況進行了分析;Ghazvinian等[7]采用水泥制作鋸齒形節(jié)理試樣并進行了不同法向荷載下的直剪實驗,對剪切峰值強度及殘余強度的演化特征進行了研究.以上研究促進了對節(jié)理剪切強度及形貌參數(shù)演化的理解,但是巖石實驗主要從宏觀角度進行分析,而不同粗糙節(jié)理內(nèi)部顆粒的相互作用及其細觀參數(shù)的影響還需進一步研究.近年來,隨著計算機技術(shù)的不斷發(fā)展,涌現(xiàn)出了多種對裂隙巖石的力學(xué)性質(zhì)和破壞行為進行模擬的數(shù)值計算方法[8-10],其中離散元數(shù)值方法尤為適合裂隙巖體壓剪狀態(tài)下的破裂及損傷演化研究.目前,PFC離散元數(shù)值模擬方法被國內(nèi)外學(xué)者廣泛接受并用來模擬脆性材料的裂紋擴展及損傷演化過程[11-13],其在巖石力學(xué)及損傷的研究成果也尤為突出,主要體現(xiàn)在裂隙巖體加載條件下的裂紋起裂、剪切破壞等[4,14-15].基于PFC二維離散元程序建立不同粗糙程度的巖石充填節(jié)理模型;通過自編FISH程序語言制定巖石節(jié)理間充填層的力學(xué)參數(shù)與性質(zhì),對充填節(jié)理的力學(xué)性質(zhì)與破壞過程進行模擬研究;分析節(jié)理JRC(粗糙度系數(shù))、接觸面強度、充填體強度與抗剪強度間的關(guān)系.

    1 節(jié)理直剪數(shù)值試驗

    1.1 數(shù)值模型的建立

    (1)

    ΔU=viΔt,

    Δθ3=(ω2-ω1)Δt,

    (2)

    (3)

    (4)

    (a) 顆粒接觸示意

    (b) 粘結(jié)力鏈分解圖1 平行粘結(jié)示意Fig.1 Parallel bond model in PFC

    新的力和力矩分別由式(5)、(6)計算得出.

    (5)

    (6)

    作用在粘結(jié)上的最大法向應(yīng)力和切向應(yīng)力為

    (7)

    在粗糙節(jié)理模型建立時,首先生成由顆粒組成的矩形.模型的整體尺寸為:100 mm(長)×100 mm(寬).矩形中的剛性顆粒由程序隨機生成,強度分布符合威布爾函數(shù),其中顆粒最小半徑為0.2 mm,顆粒半徑放大系數(shù)為1.66.模型中其余細觀參數(shù)如表1中所示.PFC模型由顆粒與顆粒粘結(jié)組成.顆粒通過粘結(jié)而形成塊體,顆粒類似于天然巖石中的礦物顆?;蚧炷林械墓攘N镔|(zhì).

    表1 顆粒細觀參數(shù)Tab.1 Microscopic parameters for rock mass

    平行粘結(jié)模型中,顆粒接觸關(guān)系如圖1(a)所示,顆粒是模型的主要框架,屬于剛體,接觸模量、接觸剛度比屬于顆粒固有性質(zhì).而平行粘結(jié)類似于天然巖石中的黏性礦物及混凝土中的水泥一樣充斥于顆粒物質(zhì)的空隙中.其中粘結(jié)半徑系數(shù)、粘結(jié)法向強度、粘結(jié)切向強度均是平行粘結(jié)內(nèi)部參數(shù).粘結(jié)半徑為1表示顆粒兩側(cè)的粘結(jié)鍵長度等于顆粒直徑.

    表1中所示參數(shù)是根據(jù)實際巖石試樣(泥質(zhì)灰?guī)r)的宏觀強度參數(shù)標定后得到的.由于巖石在加載初期內(nèi)部微裂隙受壓閉合而導(dǎo)致應(yīng)力-應(yīng)變曲線出現(xiàn)非線性增長,所以試樣的彈性模量為曲線直線段斜率.從圖2(a)不難看出,數(shù)值計算結(jié)果在峰值強度及彈性模量方面與實驗結(jié)果吻合良好.圖2(b)、2(c)中所示為數(shù)值模型最終破壞模式與實際巖石試樣單軸加載破壞狀態(tài).

    (a) 應(yīng)力-應(yīng)變曲線

    (b) 數(shù)值模擬試件破壞形態(tài)(c) 真實試件破壞形態(tài)圖2 單軸加載下完整試樣應(yīng)力-應(yīng)變曲線對比Fig.2 Comparison of stress-strain curves and failure mode with numerical model and real sample under uniaxial compression

    顯然,數(shù)值計算結(jié)果與實驗結(jié)果一致.PFC作為一款數(shù)值算法,其計算過程中材料及邊界條件均是理想化的,因此和所有數(shù)值算法一樣無法模擬出真實巖石中峰值前出現(xiàn)的軟化現(xiàn)象.在數(shù)值模擬中峰前曲線基本屬于線性或類線性,而真實試樣中存在泥質(zhì)和較軟的伴生礦物,在加載初期會出現(xiàn)一定的軟化現(xiàn)象,從而導(dǎo)致圖2(a)中所示的情況.其次,在標定細觀參數(shù)時,所取的是切線彈性模量,即應(yīng)力-應(yīng)變曲線直線段的斜率,從而在圖2(a)中可以看到,數(shù)值曲線峰前段的斜率和實驗結(jié)果中直線段斜率相同,峰值強度相同,峰值應(yīng)變不同.

    利用Barton JRC標準輪廓線(圖3(a))為參照,生成不同粗糙程度的節(jié)理模型.PFC中需要利用FISH語言自編程序來實現(xiàn)該指令.具體過程為:先將不同粗糙度的巴頓輪廓線轉(zhuǎn)化為坐標,這一過程中為了保證輪廓線的準確度,將每條輪廓線分成100段;通過編程讀取每一段前后坐標值,之后導(dǎo)入PFC中,生成與輪廓線相似的彎曲墻體;利用墻與顆粒接觸的原理,識別接觸點并記錄下接觸點所屬顆粒的ID值;將與輪廓線墻體相接觸的顆粒找出并進行重新標記.從而,生成了不同粗糙度的充填節(jié)理模型.

    如圖3(b)所示為PFC中生成的JRC為14~16的粗糙充填節(jié)理模型.模型中顏色較淺的顆粒代表巖石,中間顏色較深部分的顆粒代表充填物質(zhì).圖3(b)中右邊部分為放大區(qū)域,顆粒球心間的連線即為接觸鍵,而球外側(cè)的線條代表平行粘結(jié).節(jié)理顆粒與圍巖間的虛線輪廓為充填物質(zhì)與圍巖的接觸面.改變節(jié)理處顆粒的強度參數(shù)值與充填物質(zhì)-圍巖接觸面粘結(jié)強度值來對充填節(jié)理剪切荷載下的力學(xué)性質(zhì)進行研究與分析.雖然完整模型中的顆粒是隨機生成的,但是不同粗糙度的節(jié)理面是后期指定生成的,所以前后兩次生成的節(jié)理面一致,不對計算結(jié)果產(chǎn)生影響.

    (a) Barton JRC標準輪廓線(b) JRC為14^16的充填節(jié)理模型圖3 Barton JRC標準輪廓線與充填節(jié)理模型Fig.3 Barton JRC standard profile and the filled joint model

    1.2 直剪試驗?zāi)M

    PFC中通常采用移動控制墻來進行加載,圖4為本文中的直剪示意圖.

    圖4 直剪示意圖Fig.4 Direct shear diagram for this study

    墻體2為法向加載墻,模擬過程中該墻通過伺服機制對模型施加恒定的法向荷載.在離散元計算程序中,由伺服機制獲得的法向荷載并非完全恒定,而是在極小值的誤差范圍內(nèi)跳動,相比于所施加的法向荷載而言,該誤差完全可以忽略不計.之后,程序通過加載墻5進行切向加載,此時墻體3運動方向和速率與墻5保持一致且剪切速率保持為0.03 mm/(106step),而墻1、4及6保持固定.

    2 計算結(jié)果與分析

    2.1 充填節(jié)理剪切力學(xué)性質(zhì)與破壞分析

    本試驗中對節(jié)理模型施加1.0、1.5、2.0、2.5、3.0 MPa 5種不同的法向荷載.不同法向荷載下充填節(jié)理模型將表現(xiàn)出不同的力學(xué)性質(zhì)與破壞程度.圖5為JRC為14~16的充填節(jié)理試樣在法向荷載為1 MPa和3 MPa下,剪應(yīng)力、法向位移以及裂紋累計數(shù)的變化情況.

    從圖5中可看出:兩種不同法向荷載下,充填節(jié)理模型力學(xué)性質(zhì)存在一些差異,整體來看,隨著剪切位移不斷增加,節(jié)理剪切應(yīng)力首先呈現(xiàn)出類線性增長趨勢;當剪切應(yīng)力達到峰值時,繼續(xù)加載,剪切應(yīng)力出現(xiàn)明顯的下降趨勢,并最終在達到一定值后基本保持穩(wěn)定;雖然剪切應(yīng)力保持在一定的范圍,但仍然存在一定的波動,這一現(xiàn)象也與實際剪切曲線相吻合[6];在剪切力未達到峰值之前,隨著剪切位移不斷增加,剪切應(yīng)力增長明顯,而法向位移并未出現(xiàn)明顯的增長趨勢;剪切應(yīng)力達到峰值后,隨著剪切位移的增加,剪切應(yīng)力急劇下降,而法向位移一直處于增長趨勢,同時微裂紋的積累與法向位移呈現(xiàn)出類似的增長趨勢;隨著剪切應(yīng)力的波動,法向位移依然存在一定的增長,而微裂紋增長在剪切殘余階段基本保持穩(wěn)定.

    (1) 剪切應(yīng)力(2) 豎向位移(3) 微裂紋數(shù)量(a) 法向荷載1 MPa(1) 剪切應(yīng)力(2) 豎向位移(3) 微裂紋數(shù)量(b) 法向荷載3 MPa圖5 不同法向荷載下剪切應(yīng)力、法向位移、裂紋數(shù)量與剪切位移的關(guān)系Fig.5 Shear stress, normal displacement, and number of micro-cracks vs. shear displacement under different normal stresses

    圖6為JRC為14~16和18~20的充填節(jié)理破壞情況.圖中節(jié)理面附近聚集了大量的微裂紋,其中黑色裂紋代表拉伸型裂紋,紅色細線代表剪切型裂紋.

    從圖6中可看出:不論JRC為14~16還是18~20的充填節(jié)理,在法向荷載為3 MPa時節(jié)理面損傷程度較大;不論拉伸裂紋還是剪切型裂紋,較高法向荷載下的分布都更為密集.此現(xiàn)象主要是當法向荷載較低時,節(jié)理模型在剪切作用下出現(xiàn)較為明顯的剪脹現(xiàn)象,從而節(jié)理面附近區(qū)域在剪切過程中主要表現(xiàn)為摩擦破壞.充填物在剪切過程中產(chǎn)生變形導(dǎo)致顆粒粘結(jié)斷裂.同時,代表充填物質(zhì)的顆粒在剛度和強度上都遠小于巖石.因此,剪切過程中節(jié)理面間的充填物也將在一定程度上起到“潤滑”的作用,節(jié)理面滑移也更為明顯.隨著法向荷載的增大,節(jié)理模型法向移動受到更大的限制,即使軟弱充填物質(zhì)能增進節(jié)理剪切作用下的滑動程度,但較高法向荷載下節(jié)理面間的吻合系數(shù)較高,從而剪切荷載下節(jié)理面主要表現(xiàn)為剪斷與摩擦的復(fù)合式破壞.

    圖7為JRC為14~16與18~20的充填節(jié)理在不同法向荷載下的細觀粘結(jié)力分布情況.圖中的黑色區(qū)域代表壓應(yīng)力,紅色部分表示剪切應(yīng)力.

    從圖7中可以發(fā)現(xiàn),不同粗糙度的充填節(jié)理在剪切過程中的接觸點位置明顯不同.JRC為14~16的充填節(jié)理接觸點主要集中在模型右側(cè),該處的巖石節(jié)理面吻合度較高,此處最先產(chǎn)生破壞.相比于14~16的節(jié)理面,JRC為18~20的節(jié)理面接觸位置分布處于左側(cè).當法向荷載為3 MPa時,JRC為14~16的節(jié)理接觸位置并未發(fā)生明顯變化,只是在法向荷載較大時,接觸面積和接觸程度明顯增加,此現(xiàn)象從圖7(b)所示的力鏈分布規(guī)律中也能看出.與14~16的節(jié)理類似,法向荷載較大時,JRC為18~20的節(jié)理面接觸程度更大,壓應(yīng)力與剪應(yīng)力分布區(qū)域明顯增加,且接觸位置與分布更為均勻,模型中節(jié)理面上存在多個接觸點.在剪切作用下,接觸點附近的巖石將被剪斷,之后隨著充填物一起產(chǎn)生摩擦作用.

    (a) JRC為14^16(b) JRC為18^20圖6 法向荷載為1MPa和3MPa下的節(jié)理損傷情況Fig.6 Joint morphology damage at normal stresses of 1 MPa and 3 MPa

    (a) JRC為14^16(b) JRC為18^20圖7 法向荷載為1 MPa和3 MPa下的細觀粘結(jié)力分布Fig.7 Distribution of bond stress at normal stresses of 1 MPa and 3 MPa

    2.2 宏細觀參數(shù)分析

    為了研究巖石節(jié)理間充填物性質(zhì)及JRC對充填節(jié)理力學(xué)性質(zhì)的影響,進行3種不同因素下的節(jié)理模型直剪數(shù)值實驗.影響因素包括:JRC、充填節(jié)理顆粒強度、充填物與巖石接觸面強度.其中JRC分為2~4、6~8、10~12、14~16和18~20共5個區(qū)間.在顆粒流中,顆粒的強度參數(shù)包括粘結(jié)強度、顆粒摩擦因數(shù).而本文中的節(jié)理顆粒強度比是代表節(jié)理充填物顆粒的剛度、粘結(jié)強度、摩擦因數(shù)與代表完整巖樣的顆粒所對應(yīng)的強度參數(shù)的比值.本文中5種節(jié)理顆粒強度摩擦因數(shù)比分別為:0.15、0.30、0.45、0.60、0.75.接觸面是指充填物質(zhì)與圍巖接觸面,而接觸面粘結(jié)強度為節(jié)理顆粒與圍巖顆粒間的粘結(jié)強度與完整巖樣顆粒間粘結(jié)強度的比值.不同組內(nèi)的相關(guān)細觀參數(shù)如表2所示.

    表2 各工況組內(nèi)的影響因素及數(shù)值Tab.2 Factors influencing the simulation cases

    節(jié)理模型中原巖的細觀參數(shù)是根據(jù)實際巖石宏觀力學(xué)性質(zhì)進行標定的.充填物質(zhì)的強度則是根據(jù)原巖參數(shù)不斷弱化確定的,進行這一設(shè)置主要是考慮到實際工程當中節(jié)理間的充填物可能是低粘結(jié)的類砂質(zhì)碎屑,也有可能是完整的侵入巖體,只是相比之下強度稍弱.所以進行數(shù)值模擬時,基于原巖強度參數(shù)對節(jié)理處顆粒參數(shù)進行弱化處理,且弱化處理過程中將粘結(jié)及摩擦因數(shù)進行同等弱化.

    根據(jù)表2中A組所示的細觀參數(shù)對充填節(jié)理力學(xué)性質(zhì)開展研究,主要分析節(jié)理粗糙程度對峰值剪切強度及剪切強度參數(shù)的影響.圖8為不同JRC影響下的節(jié)理剪切力學(xué)性質(zhì),由圖8(a)可以看出:隨法向荷載不斷增加,節(jié)理的峰值剪切強度呈線性增長,與室內(nèi)實驗結(jié)果相似[6];在同等法向荷載下,隨著節(jié)理JRC的不斷增加,節(jié)理模型的峰值剪切強度呈現(xiàn)出明顯的增長趨勢;隨著節(jié)理粗糙度的不斷增大其抗剪強度也不斷增大.圖8(b)為節(jié)理抗剪強度參數(shù)隨著JRC的演化情況,隨節(jié)理JRC不斷增大,抗剪強度參數(shù)Cj(粘聚力)和φj(內(nèi)摩擦角)均呈增大的趨勢,與黃飄[17]試驗結(jié)果的規(guī)律相同,相比而言Cj值的增長幅度較小,在JRC值超過12后,φj值出現(xiàn)明顯的增長.

    根據(jù)表2中B組進行數(shù)值計算可以得出節(jié)理間充填物強度對其力學(xué)性質(zhì)的影響,圖9為充填物強度對節(jié)理剪切力學(xué)性質(zhì)的影響情況.與圖8類似,圖9(a)為不同法向荷載下的節(jié)理峰值剪切強度.不同的充填強度比下,隨著法向荷載增大,各組節(jié)理峰值剪切強度呈線性變化.當節(jié)理所處法向荷載一定時,隨充填物強度的增加其峰值剪切強度也出現(xiàn)明顯的增長趨勢.節(jié)理充填強度比為0.15時峰值剪切強度最小,0.75時最大.圖9(b)為不同充填強度對節(jié)理抗剪強度參數(shù)Cj和φj的影響.隨充填物強度比的增加,粘結(jié)力出現(xiàn)了明顯的增長,從0.15時的0.92 MPa增長至0.75時的4.3 MPa.由此可見,充填物強度比對節(jié)理面剪切粘結(jié)力存在較大的影響.隨著充填強度比的增加,節(jié)理的抗剪強度參數(shù)φj呈先下降后增加的趨勢.

    當節(jié)理強度比為0.3時,節(jié)理面的φj處于最低值.隨著充填強度比的不斷增加,與參數(shù)值Cj一樣,內(nèi)摩擦角φj同樣呈明顯的增長趨勢.表2中C組所示即為改變接觸面粘結(jié)強度時各個參數(shù)的分布和變化情況.其中接觸面強度是變量,而JRC固定為2~4,節(jié)理顆粒強度比定為 0.3.

    圖10中為充填物與巖石接觸面粘結(jié)情況對峰值剪切應(yīng)力及抗剪強度參數(shù)的影響.隨著法向荷載的不斷增長,JRC為2~4的節(jié)理面峰值剪切應(yīng)力明顯增長.反觀接觸面粘結(jié)強度對峰值剪切強度的影響,當法向荷載較小時,不同接觸面的粘結(jié)強度比下峰值剪切強度差距較大,可見較低法向荷載下接觸面粘結(jié)強度比對峰值剪切應(yīng)力的影響較大.當法向荷載較大時,接觸面粘結(jié)強度對峰值剪切應(yīng)力的影響程度明顯降低.主要是由于法向荷載較大時上下節(jié)理面被充分壓縮,接觸程度較大,剪切荷載作用下剪切應(yīng)力主要受到粗糙度的控制和影響.

    (a) 法向荷載對峰值強度的影響(b) JRC對粘聚力的影響(c) JRC對內(nèi)摩擦角的影響圖8 JRC對節(jié)理力學(xué)性質(zhì)的影響Fig.8 Influence of JRC value on the mechanical properties of the fill joint

    (a) 法向荷載對峰值強度的影響(b) 充填體強度比對粘聚力的影響(c) 充填體強度比對內(nèi)摩擦角的影響圖9 節(jié)理充填物強度比對節(jié)理力學(xué)性質(zhì)的影響Fig.9 Influence of the filling strength ratio on the mechanical properties of the fill joint

    (a) 法向荷載對峰值強度的影響(b) 接觸面粘結(jié)強度比對粘聚力的影響(c) 接觸面粘結(jié)強度比對內(nèi)摩擦角的影響圖10 接觸面強度比對節(jié)理力學(xué)性質(zhì)的影響Fig.10 Influence of the interface strength ratio on the mechanical properties of the fill joint

    從圖10(b)可以看出,節(jié)理面抗剪強度參數(shù)Cj隨接觸面粘結(jié)的增加同樣呈增長的趨勢,而抗剪強度參數(shù)φj隨接觸面粘結(jié)強度增大而減小,與Cj的變化規(guī)律相反.由此可見,充填節(jié)理接觸面粘結(jié)強度對節(jié)理抗剪強度參數(shù)存在不同程度的影響.

    除前述的幾種影響因素外,充填厚度也對充填節(jié)理剪切力學(xué)性質(zhì)產(chǎn)生極為重要的影響.表2中的D組即為5種不同充填厚度下的直剪模擬方案.此組模擬方案當中,節(jié)理粗糙度選為JRC處于10~12的節(jié)理面.節(jié)理顆粒強度比為0.3,接觸面粘結(jié)強度比為1.0.充填厚度作為唯一變量,分別為 0.4、1.4、3.0、4.0 mm以及5.0 mm.

    圖11中為不同充填厚度下的節(jié)理峰值剪切強度、抗剪強度參數(shù)以及粘結(jié)破壞數(shù)量分布規(guī)律.

    從圖11(a)中看出:隨著法向荷載不斷增加,各充填厚度下的節(jié)理峰值抗剪強度依然呈線性變化;隨著充填厚度的增加,節(jié)理峰值抗剪強度越小;充填厚度為0.4 mm時,各法向荷載下的峰值抗剪強度處于最大值,而厚度為5 mm時均處于最小值;隨著充填厚度的增加,峰值抗剪強度并未線性遞減.

    圖11(b)中所示為充填厚度的變化對抗剪強度參數(shù)的影響.隨著充填厚度的增加,節(jié)理面抗剪強度參數(shù)存在明顯的劣化趨勢.節(jié)理面抗剪強度參數(shù)Cj與φj均隨著充填厚度的增加而減小.節(jié)理面抗剪強度參數(shù)Cj與填厚度呈類線性關(guān)系,而φj呈非線性遞減的關(guān)系.充填節(jié)理的力學(xué)性質(zhì)受充填厚度的影響較大,還要考慮到節(jié)理面起伏度的影響.對于JRC為10~12的模型而言,當充填厚度為5 mm 時,充填度大于100%,此時隨著充填厚度的增加節(jié)理強度主要受充填體本身強度的影響.因此,圖11中不論是峰值強度還是抗剪強度參數(shù)Cj與φj,當充填厚度從4 mm增加至5 mm時,強度參數(shù)變化程度不如前面幾種厚度變化顯著.

    (a) 法向荷載對峰值強度的影響(b) 充填體厚度對粘聚力的影響(c) 充填體厚度對內(nèi)摩擦角的影響圖11 充填體厚度對節(jié)理力學(xué)性質(zhì)的影響Fig.11 Influence of filling thickness on the mechanical properties of the fill joint

    3 結(jié) 論

    利用二維離散元程序建立不同粗糙程度的巖石充填節(jié)理模型,分析JRC、充填強度、接觸面粘結(jié)強度、充填厚度等參數(shù)對節(jié)理剪切強度的影響,主要結(jié)論如下:

    (1) 法向荷載較大時,上下節(jié)理面接觸點分布較為均勻,節(jié)理峰值抗剪應(yīng)力更大,節(jié)理面損傷程度也愈發(fā)明顯.

    (2) 隨著JRC的增加,節(jié)理峰值剪切強度明顯增長,而當法向荷載較大時,這種影響將被放大.節(jié)理面抗剪強度參數(shù)Cj隨著JRC的增加呈類線性增長,而φj隨著JRC的增大呈非線性增長.

    (3) 節(jié)理峰值剪切應(yīng)力隨節(jié)理充填物的強度增加而增大.隨充填物強度比的增加,粘結(jié)力Cj出現(xiàn)明顯的增長趨勢,隨著充填強度比的增加,節(jié)理的抗剪強度參數(shù)φj呈先下降后增加.

    (4) 法向荷載較小時,接觸面粘結(jié)強度比對剪切強度的影響較大;當法向荷載較大時,接觸面粘結(jié)強度對剪切強度的影響程度明顯降低.抗剪強度參數(shù)Cj隨接觸面粘結(jié)強度的增加呈現(xiàn)增長的趨勢,而抗剪強度參數(shù)φj隨著接觸面粘結(jié)強度增大而減小.

    (5) 充填厚度對節(jié)理剪切力學(xué)特性存在明顯的影響,隨著充填厚度的增加,峰值強度、抗剪強度參數(shù)Cj和φj逐漸減小.當充填厚度增加至一定值后,充填節(jié)理剪切力學(xué)性質(zhì)主要受充填體本身的強度影響.

    致謝:中南大學(xué)“創(chuàng)新驅(qū)動計劃”項目(2016CX019);深部巖土力學(xué)與地下工程國家重點實驗室開放基金項目(SKLGDUEK1405).

    [1] 夏才初,孫宗頎. 工程巖體節(jié)理力學(xué)[M]. 同濟大學(xué)出版社,2002:1-2.

    [2] YANG S Q, XU T, H E L, et al. Numerical study on failure behavior of brittle rock specimen containing pre-existing combined flaws under different confining pressure[J]. Archives of Civil and Mechanical Engineering, 2015, 15(4): 1085-1097.

    [3] 蔣宇靜,王剛,李博,等. 巖石節(jié)理剪切滲流耦合試驗及分析[J]. 巖石力學(xué)與工程學(xué)報,2007,26(11): 2253-2259.

    JIANG Yujing, WANG Gang, LI bo, et al. Experimental study and analysis of shear-flow coupling behaviors of rock joints[J]. Chinese Journal of Rock Mechanicss and Engineering, 2007, 26(11): 2253-2259.

    [4] 周喻,MISAR A,吳順川,等. 巖石節(jié)理直剪試驗顆粒流宏細觀分析[J]. 巖石力學(xué)與工程學(xué)報,2012,31(6): 1245-1256.

    ZHOU Yu, MISRA A, WU Shunchuan, et al. Macro-and meso-analysis of rock joint direct shear test using particle flow theory[J]. Chinese Journal of Rock Mechanics and Engineering, 2012, 31(6): 1245-1256.

    [5] 唐雯鈺,林杭. 不同鋸齒高度對軟弱結(jié)構(gòu)面剪切特性的影響[J]. 中南大學(xué)學(xué)報:自然科學(xué)版,2017,48(5): 1300-1307.

    TANG Wenyu, LIN Hang. Influence of dentate discontinuity height on shear properties of soft structure plane[J]. Journal of Central South University: Science and Technology, 2017, 48(5): 1300-1307.

    [6] 曹平,范祥,蒲成志,等. 節(jié)理剪切試驗及其表面形貌特征變化分析[J]. 巖石力學(xué)與工程學(xué)報,2011,30(3): 480-485.

    CAO Ping, FAN Xiang, PU Chengzhi, et al. Shear test of joint and analysis of morohology characteristic evolution of joint surface[J]. Chinese Journal of Rock Mechanics and Engineering, 2011, 30(3): 480-485.

    [7] GHAZVINIAN A, AZINFAR MJ, NOROZI P. Mechanical response of discontinuities of different joint wall contact strengths[J]. Arabian Journal of Geosciences, 2013, 6(11): 4419-4426.

    [8] 陳必光,宋二祥,程曉輝. 二維裂隙巖體滲流傳熱的離散裂隙網(wǎng)絡(luò)模型數(shù)值計算方法[J]. 巖石力學(xué)與工程學(xué)報,2014(1): 43-51.

    CHEN Biguang, SONG Erxiang, CHENG Xiaohui. A numerical method for discrete fracture network model for flow and heat transfer in two-dimensional fractured rocks[J]. Chinese Journal of Rock Mechanics and Engineering, 2014(1): 43-51.

    [9] 杜守繼,職洪濤,周枝華. 巖石節(jié)理剪切過程中應(yīng)力與滲流特性的數(shù)值模擬[J]. 巖石力學(xué)與工程學(xué)報,2008,27(12): 2473-2481.

    DU Shouji, ZHI Hongtao, ZHOU Zhihua. Numerical simulation of stress and seepage properties of rock joint during shearing processes[J]. Chinese Journal of Rock Mechanics and Engineering, 2008, 27(12): 2473-2481.

    [10] BAHAADDINI M, SHARROCK G, HEBBLEWHITE B K. Numerical investigation of the effect of joint geometrical parameters on the mechanical properties of a non-persistent jointed rock mass under uniaxial compression[J]. Computers and Geotechnics, 2013, 49(1): 206-225.

    [11] 吳順川,周喻,高斌. 卸載巖爆試驗及PFC3D數(shù)值模擬研究[J]. 巖石力學(xué)與工程學(xué)報,2010,29(S2): 4082-4088.

    WU Shunchuan, ZHOU Yu, GAO Bin. Study of unloading tests of rock burst and PFC3D numerical simulation[J]. Chinese Journal of Rock Mechanicss and Engineering, 2010, 29(S2): 4082-4088.

    [12] MANOUCHEHRIAN A, SHARIFZADEH M, MARJI M F, et al. A bonded particle model for analysis of the flaw orientation effect on crack propagation mechanism in brittle materials under compression[J]. Archives of Civil and Mechanical Engineering, 2014, 14(1): 40-52.

    [13] 黃達,岑奪豐. 單軸靜-動相繼壓縮下單裂隙巖樣力學(xué)響應(yīng)及能量耗散機制顆粒流模擬[J]. 巖石力學(xué)與工程學(xué)報,2013(9): 1926-1936.

    HUANG Da, CEN Duofeng. Mechanical responses and energy dissipation mechanism of rock specimen with a single fissure under static and dynamic unaxial compression using particle flow code simulations[J]. Chinese Journal of Rock Mechanics and Engineering, 2013(9): 1926-1936.

    [14] PARK J W, SONG J J. Numerical simulation of a direct shear test on a rock joint using a bonded-particle model[J]. International Journal of Rock Mechanics and Mining Sciences, 2009, 46(8): 1315-1328.

    [15] SARFARAZI V, GHAZVINIAN A, SCHUBERT W. Numerical simulation of the process of fracture of echelon rock joints[J]. Rock Mechanics and Rock Engineering, 2014, 47(4): 1355-1371.

    [16] Itasca Consulting Group. PFC2D (partical flow code in 2 dimensions) theory and background[M]. Minnesota: Itasca Consulting Group Inc, 2002: 77-78.

    [17] 黃飄. 不同粗糙度結(jié)構(gòu)面試件剪切試驗特性研究[J]. 建筑監(jiān)督檢測與造價,2015,8(3): 1-5.

    HUANG Piao. Study on shear test properties of different structural plane roughness[J]. Supervision Test and Cost of Construction, 2015, 8(3): 1-5.

    猜你喜歡
    剪切應(yīng)力法向節(jié)理
    落石法向恢復(fù)系數(shù)的多因素聯(lián)合影響研究
    新疆阜康白楊河礦區(qū)古構(gòu)造應(yīng)力場特征
    心瓣瓣膜區(qū)流場中湍流剪切應(yīng)力對瓣膜損害的研究進展
    新疆阜康白楊河礦區(qū)構(gòu)造節(jié)理發(fā)育特征
    中國煤層氣(2018年3期)2018-07-27 11:41:42
    Effect of Magnetic Field on Forced Convection between Two Nanofluid Laminar Flows in a Channel
    低溫狀態(tài)下的材料法向發(fā)射率測量
    剪切應(yīng)力對聚乳酸結(jié)晶性能的影響
    中國塑料(2016年6期)2016-06-27 06:34:24
    落石碰撞法向恢復(fù)系數(shù)的模型試驗研究
    動脈粥樣硬化病變進程中血管細胞自噬的改變及低剪切應(yīng)力對血管內(nèi)皮細胞自噬的影響*
    硫化氫在低剪切應(yīng)力導(dǎo)致內(nèi)皮細胞自噬障礙中的作用
    日韩精品青青久久久久久| 成年av动漫网址| 黄片wwwwww| 色视频www国产| 国产国拍精品亚洲av在线观看| 日韩大尺度精品在线看网址| 国产欧美日韩一区二区精品| 国产精品一及| 级片在线观看| 五月伊人婷婷丁香| 看十八女毛片水多多多| 日本欧美国产在线视频| 18禁黄网站禁片免费观看直播| 精品一区二区三区人妻视频| 国产黄色视频一区二区在线观看 | 男女视频在线观看网站免费| 一a级毛片在线观看| 亚洲人成网站在线播放欧美日韩| 丝袜喷水一区| 乱系列少妇在线播放| 美女cb高潮喷水在线观看| 日日撸夜夜添| 少妇丰满av| 五月伊人婷婷丁香| 国产一区二区在线观看日韩| 日韩av不卡免费在线播放| 日本黄色视频三级网站网址| 色在线成人网| 国产精品久久久久久久久免| 精品人妻熟女av久视频| 97热精品久久久久久| 露出奶头的视频| 菩萨蛮人人尽说江南好唐韦庄 | 91av网一区二区| av.在线天堂| 国产探花极品一区二区| 亚洲人成网站高清观看| 欧美中文日本在线观看视频| 亚洲va在线va天堂va国产| 听说在线观看完整版免费高清| 高清日韩中文字幕在线| 国产精品福利在线免费观看| 白带黄色成豆腐渣| 国产探花极品一区二区| 免费一级毛片在线播放高清视频| 尤物成人国产欧美一区二区三区| 99热全是精品| 禁无遮挡网站| 嫩草影院新地址| 久久精品影院6| 看非洲黑人一级黄片| 永久网站在线| 久久久久国产精品人妻aⅴ院| 日韩一区二区视频免费看| 99热这里只有精品一区| 桃色一区二区三区在线观看| 午夜爱爱视频在线播放| 亚洲高清免费不卡视频| 日韩成人av中文字幕在线观看 | 国产v大片淫在线免费观看| 成人亚洲精品av一区二区| 精品99又大又爽又粗少妇毛片| 日本在线视频免费播放| 国产精品久久电影中文字幕| 直男gayav资源| 亚洲在线自拍视频| 亚洲在线观看片| 少妇的逼水好多| 欧美绝顶高潮抽搐喷水| 亚洲无线观看免费| 亚洲欧美日韩东京热| 国产精品乱码一区二三区的特点| 免费无遮挡裸体视频| 夜夜看夜夜爽夜夜摸| 久久热精品热| 男女啪啪激烈高潮av片| 波多野结衣高清无吗| 大型黄色视频在线免费观看| 91久久精品电影网| 两个人的视频大全免费| 日本免费一区二区三区高清不卡| 中文亚洲av片在线观看爽| 国产精品亚洲一级av第二区| 色哟哟哟哟哟哟| 久久热精品热| 老熟妇仑乱视频hdxx| 22中文网久久字幕| 老熟妇仑乱视频hdxx| 精品人妻偷拍中文字幕| 麻豆久久精品国产亚洲av| av在线蜜桃| 成人av一区二区三区在线看| 午夜福利高清视频| 亚洲精品日韩在线中文字幕 | 日韩国内少妇激情av| 免费av不卡在线播放| 伦理电影大哥的女人| 亚洲经典国产精华液单| 国语自产精品视频在线第100页| 日本成人三级电影网站| 亚洲人成网站在线播放欧美日韩| 日本成人三级电影网站| 国产爱豆传媒在线观看| 嫩草影视91久久| 精华霜和精华液先用哪个| 草草在线视频免费看| 精品国产三级普通话版| 99久久中文字幕三级久久日本| 久久久久久久久久黄片| 精品久久国产蜜桃| 亚洲七黄色美女视频| 在线观看美女被高潮喷水网站| 国产成人a∨麻豆精品| 国产精品久久久久久av不卡| videossex国产| 最后的刺客免费高清国语| 简卡轻食公司| av女优亚洲男人天堂| 夜夜夜夜夜久久久久| 国产精品永久免费网站| 嫩草影视91久久| 永久网站在线| 在线观看av片永久免费下载| 午夜福利在线在线| 国产精品久久久久久精品电影| 五月玫瑰六月丁香| 精品乱码久久久久久99久播| 中文亚洲av片在线观看爽| 国产精品国产高清国产av| 成人欧美大片| 国产亚洲精品av在线| 国产探花极品一区二区| 精品人妻偷拍中文字幕| 久久这里只有精品中国| 美女黄网站色视频| 日韩成人av中文字幕在线观看 | 欧洲精品卡2卡3卡4卡5卡区| 中文字幕久久专区| 亚洲,欧美,日韩| 一级黄色大片毛片| 熟女电影av网| 日韩av不卡免费在线播放| 国产伦精品一区二区三区视频9| 免费不卡的大黄色大毛片视频在线观看 | 秋霞在线观看毛片| 高清毛片免费看| 国产欧美日韩精品一区二区| 欧美不卡视频在线免费观看| 欧美不卡视频在线免费观看| 久久热精品热| 一进一出好大好爽视频| 日本爱情动作片www.在线观看 | 人妻制服诱惑在线中文字幕| 麻豆久久精品国产亚洲av| 91麻豆精品激情在线观看国产| 日本a在线网址| 秋霞在线观看毛片| 老女人水多毛片| 国产精品国产高清国产av| 亚洲精品日韩在线中文字幕 | 久久99热这里只有精品18| 男女下面进入的视频免费午夜| 国产真实伦视频高清在线观看| 日韩人妻高清精品专区| 国产免费一级a男人的天堂| 国产私拍福利视频在线观看| 亚洲第一电影网av| 精品久久久久久成人av| 午夜福利成人在线免费观看| 欧美极品一区二区三区四区| 国产精品野战在线观看| 人人妻人人澡人人爽人人夜夜 | 国产大屁股一区二区在线视频| 麻豆一二三区av精品| 亚洲,欧美,日韩| 精品熟女少妇av免费看| 自拍偷自拍亚洲精品老妇| 色综合亚洲欧美另类图片| 国产亚洲91精品色在线| 日本爱情动作片www.在线观看 | 三级经典国产精品| 亚洲一区二区三区色噜噜| 欧美+亚洲+日韩+国产| 精品久久久久久久久亚洲| 日韩欧美一区二区三区在线观看| 在线免费观看的www视频| 天堂√8在线中文| 中国国产av一级| 男女之事视频高清在线观看| 欧美xxxx黑人xx丫x性爽| 春色校园在线视频观看| 成人av在线播放网站| 99热全是精品| 一级av片app| 久久精品国产鲁丝片午夜精品| 人人妻人人看人人澡| 变态另类成人亚洲欧美熟女| 免费看a级黄色片| 丰满的人妻完整版| 国产精品爽爽va在线观看网站| 中国美女看黄片| 在线免费观看不下载黄p国产| 国产久久久一区二区三区| 国产乱人偷精品视频| 亚洲国产日韩欧美精品在线观看| 三级经典国产精品| 精品一区二区三区人妻视频| 91av网一区二区| 久久人人爽人人爽人人片va| 我要搜黄色片| 亚洲婷婷狠狠爱综合网| 中出人妻视频一区二区| 最近最新中文字幕大全电影3| 午夜免费男女啪啪视频观看 | 国产欧美日韩精品亚洲av| 欧美日韩乱码在线| 毛片一级片免费看久久久久| 日韩欧美精品免费久久| 欧美zozozo另类| 美女内射精品一级片tv| 一本久久中文字幕| 又爽又黄无遮挡网站| 一级毛片aaaaaa免费看小| 变态另类成人亚洲欧美熟女| av在线观看视频网站免费| 97超视频在线观看视频| 一级毛片久久久久久久久女| 在现免费观看毛片| 99久国产av精品| 真人做人爱边吃奶动态| 全区人妻精品视频| 国产高清激情床上av| 免费观看人在逋| 亚洲欧美清纯卡通| 男人狂女人下面高潮的视频| 老女人水多毛片| 99久久九九国产精品国产免费| 亚洲av中文av极速乱| 国产高清激情床上av| 我要看日韩黄色一级片| 久久久色成人| 国内久久婷婷六月综合欲色啪| 18+在线观看网站| 成人亚洲欧美一区二区av| av天堂在线播放| 亚洲18禁久久av| 欧美日韩乱码在线| 别揉我奶头~嗯~啊~动态视频| 免费看av在线观看网站| 国产在视频线在精品| 精品久久久久久久人妻蜜臀av| 国语自产精品视频在线第100页| 亚洲第一区二区三区不卡| 99国产极品粉嫩在线观看| 亚洲熟妇熟女久久| 免费人成视频x8x8入口观看| 亚洲av不卡在线观看| 国产亚洲精品综合一区在线观看| 在线观看免费视频日本深夜| 国产高清激情床上av| 一级毛片我不卡| 久久精品国产亚洲av天美| 久久鲁丝午夜福利片| 日日撸夜夜添| 亚洲aⅴ乱码一区二区在线播放| 国产亚洲欧美98| 赤兔流量卡办理| 久久欧美精品欧美久久欧美| 欧美一区二区亚洲| 午夜精品在线福利| 欧洲精品卡2卡3卡4卡5卡区| 中国美白少妇内射xxxbb| 寂寞人妻少妇视频99o| 久久精品久久久久久噜噜老黄 | 精品无人区乱码1区二区| 十八禁国产超污无遮挡网站| 国产精品亚洲美女久久久| 亚洲欧美精品自产自拍| 欧美日韩在线观看h| 亚洲欧美日韩卡通动漫| 桃色一区二区三区在线观看| 欧美绝顶高潮抽搐喷水| 日韩一本色道免费dvd| 亚洲欧美成人综合另类久久久 | 啦啦啦韩国在线观看视频| 在线观看免费视频日本深夜| 精品一区二区免费观看| 天天躁日日操中文字幕| 精品人妻一区二区三区麻豆 | 午夜福利在线在线| 丰满乱子伦码专区| 国产人妻一区二区三区在| 亚洲国产色片| 精品人妻熟女av久视频| 亚洲三级黄色毛片| av天堂中文字幕网| 免费电影在线观看免费观看| 嫩草影院入口| 国产精品永久免费网站| 亚洲,欧美,日韩| av天堂在线播放| 男女边吃奶边做爰视频| 性插视频无遮挡在线免费观看| 校园春色视频在线观看| 亚洲中文日韩欧美视频| 国产免费一级a男人的天堂| 久久久久久国产a免费观看| 亚洲精品国产av成人精品 | 免费av观看视频| 在线a可以看的网站| 日韩一本色道免费dvd| 久久婷婷人人爽人人干人人爱| 国产精品日韩av在线免费观看| 精品免费久久久久久久清纯| 日本欧美国产在线视频| 欧美高清性xxxxhd video| 国产色婷婷99| 午夜福利在线观看免费完整高清在 | 亚洲图色成人| 亚洲性久久影院| 国产三级中文精品| 精品午夜福利在线看| 国产成人影院久久av| 又爽又黄a免费视频| 又黄又爽又免费观看的视频| 国产精品福利在线免费观看| 美女 人体艺术 gogo| 99久久成人亚洲精品观看| 亚洲av五月六月丁香网| 午夜a级毛片| 69av精品久久久久久| 一卡2卡三卡四卡精品乱码亚洲| 色5月婷婷丁香| 老熟妇仑乱视频hdxx| 亚洲国产欧洲综合997久久,| 国产亚洲av嫩草精品影院| 成人二区视频| 不卡视频在线观看欧美| 日本黄大片高清| avwww免费| 我要看日韩黄色一级片| 99久久精品国产国产毛片| 又爽又黄无遮挡网站| 国产伦精品一区二区三区四那| 亚洲人成网站在线播放欧美日韩| 欧美bdsm另类| 成人漫画全彩无遮挡| 国产成人影院久久av| 真实男女啪啪啪动态图| 欧美三级亚洲精品| 九色成人免费人妻av| 亚洲欧美中文字幕日韩二区| 人人妻人人看人人澡| 国产一区二区三区av在线 | 18禁黄网站禁片免费观看直播| 乱人视频在线观看| 99热这里只有是精品在线观看| a级一级毛片免费在线观看| 亚洲激情五月婷婷啪啪| 亚洲熟妇中文字幕五十中出| 国产精华一区二区三区| 亚洲av不卡在线观看| 成人高潮视频无遮挡免费网站| 精品人妻视频免费看| a级毛片a级免费在线| 一级黄片播放器| 午夜免费男女啪啪视频观看 | 性色avwww在线观看| 天堂网av新在线| 两个人的视频大全免费| 欧美一区二区国产精品久久精品| 在现免费观看毛片| 一级毛片aaaaaa免费看小| 人人妻,人人澡人人爽秒播| 亚洲av成人av| 狂野欧美白嫩少妇大欣赏| 91在线观看av| 最好的美女福利视频网| 超碰av人人做人人爽久久| 欧美色欧美亚洲另类二区| 中文字幕av在线有码专区| 欧美日韩在线观看h| av卡一久久| 国产在视频线在精品| 国产熟女欧美一区二区| 国产一区二区在线av高清观看| 麻豆成人午夜福利视频| 一本精品99久久精品77| 一个人看的www免费观看视频| 搡老岳熟女国产| 久久精品综合一区二区三区| 日韩亚洲欧美综合| 在线播放国产精品三级| 国产精品一区二区三区四区免费观看 | 热99在线观看视频| 99久国产av精品国产电影| 又黄又爽又刺激的免费视频.| 身体一侧抽搐| 久久精品综合一区二区三区| 国产精品伦人一区二区| 亚洲经典国产精华液单| 亚洲熟妇中文字幕五十中出| 久久久久性生活片| 熟妇人妻久久中文字幕3abv| 色在线成人网| 一夜夜www| 又黄又爽又刺激的免费视频.| 大型黄色视频在线免费观看| 色播亚洲综合网| 精品一区二区三区av网在线观看| 蜜臀久久99精品久久宅男| 国产精品女同一区二区软件| 日韩 亚洲 欧美在线| 成人永久免费在线观看视频| 五月伊人婷婷丁香| 男女视频在线观看网站免费| 精品久久久久久久久久免费视频| 国产探花极品一区二区| 久久婷婷人人爽人人干人人爱| 一本一本综合久久| 黄色一级大片看看| 精品久久久久久久久av| 欧美xxxx性猛交bbbb| 国产高清不卡午夜福利| 国产一区二区亚洲精品在线观看| 日本三级黄在线观看| 2021天堂中文幕一二区在线观| 亚洲欧美日韩高清专用| 看免费成人av毛片| 18禁在线播放成人免费| 大又大粗又爽又黄少妇毛片口| 国产真实伦视频高清在线观看| 国产淫片久久久久久久久| 小说图片视频综合网站| 亚洲av一区综合| 黄色欧美视频在线观看| 日韩欧美 国产精品| 少妇裸体淫交视频免费看高清| 91精品国产九色| 中文亚洲av片在线观看爽| 国产国拍精品亚洲av在线观看| 亚洲成人久久爱视频| 欧美激情国产日韩精品一区| 一区福利在线观看| 一级毛片我不卡| 联通29元200g的流量卡| 亚洲av中文av极速乱| 欧美日韩乱码在线| 亚洲美女搞黄在线观看 | 乱人视频在线观看| 亚洲美女搞黄在线观看 | 国内精品宾馆在线| 日韩,欧美,国产一区二区三区 | 亚洲精品国产成人久久av| 在线看三级毛片| av福利片在线观看| 中国国产av一级| 99热精品在线国产| 日本与韩国留学比较| 性色avwww在线观看| 在线国产一区二区在线| 最新中文字幕久久久久| 亚洲va在线va天堂va国产| 99久久精品一区二区三区| 我的女老师完整版在线观看| 一区二区三区四区激情视频 | 欧美色欧美亚洲另类二区| 免费观看在线日韩| 一个人观看的视频www高清免费观看| 99精品在免费线老司机午夜| 国内揄拍国产精品人妻在线| 国产淫片久久久久久久久| 国产午夜精品论理片| 国产精品av视频在线免费观看| 久久午夜亚洲精品久久| 午夜视频国产福利| 3wmmmm亚洲av在线观看| 国产淫片久久久久久久久| 久久亚洲国产成人精品v| 男人舔女人下体高潮全视频| 日本一二三区视频观看| 在线播放无遮挡| 久久精品夜夜夜夜夜久久蜜豆| 国产成人影院久久av| 伦精品一区二区三区| 99热网站在线观看| 校园春色视频在线观看| 国产真实伦视频高清在线观看| 人妻少妇偷人精品九色| 精品熟女少妇av免费看| 国内精品宾馆在线| 国产高清不卡午夜福利| 69av精品久久久久久| 国产伦在线观看视频一区| 搡女人真爽免费视频火全软件 | 日韩大尺度精品在线看网址| 国产蜜桃级精品一区二区三区| 99热这里只有是精品在线观看| 精品不卡国产一区二区三区| 久久久精品欧美日韩精品| 最近最新中文字幕大全电影3| 亚洲精品亚洲一区二区| 欧美在线一区亚洲| 美女被艹到高潮喷水动态| 日韩中字成人| 久久国内精品自在自线图片| 搡老熟女国产l中国老女人| 国产精品一及| 国产精品一二三区在线看| 小说图片视频综合网站| 亚洲最大成人av| 欧美成人免费av一区二区三区| 你懂的网址亚洲精品在线观看 | 在线播放国产精品三级| 在线观看av片永久免费下载| 亚洲av第一区精品v没综合| 久久久精品大字幕| 91av网一区二区| 99热全是精品| 国产高清视频在线播放一区| 久久午夜福利片| 两性午夜刺激爽爽歪歪视频在线观看| 露出奶头的视频| 免费人成在线观看视频色| 丰满人妻一区二区三区视频av| 99热6这里只有精品| 全区人妻精品视频| 最近中文字幕高清免费大全6| 久久久精品大字幕| 国产蜜桃级精品一区二区三区| 国产亚洲欧美98| 日韩欧美 国产精品| 97热精品久久久久久| 亚洲欧美日韩东京热| 噜噜噜噜噜久久久久久91| 精品久久久噜噜| 久久亚洲精品不卡| 99热全是精品| 久久久久久九九精品二区国产| 在线观看午夜福利视频| 亚洲欧美日韩东京热| 久久久久久大精品| 成人美女网站在线观看视频| 日韩精品有码人妻一区| 51国产日韩欧美| 日本黄色视频三级网站网址| 日本爱情动作片www.在线观看 | 亚洲成人精品中文字幕电影| av女优亚洲男人天堂| 波野结衣二区三区在线| 免费人成视频x8x8入口观看| 久久久久国产精品人妻aⅴ院| 久久久精品欧美日韩精品| 成人一区二区视频在线观看| 日韩欧美三级三区| 国产精品久久电影中文字幕| 搡老岳熟女国产| 国产精品一二三区在线看| 寂寞人妻少妇视频99o| 插阴视频在线观看视频| 网址你懂的国产日韩在线| 日日干狠狠操夜夜爽| 国产午夜精品论理片| 长腿黑丝高跟| 老师上课跳d突然被开到最大视频| 国产不卡一卡二| 中文字幕免费在线视频6| 国产一区二区三区av在线 | 午夜精品在线福利| 波野结衣二区三区在线| 免费人成视频x8x8入口观看| 久久鲁丝午夜福利片| 精品免费久久久久久久清纯| 精品少妇黑人巨大在线播放 | 18+在线观看网站| 少妇的逼好多水| 免费观看精品视频网站| 一级毛片电影观看 | 看十八女毛片水多多多| 99九九线精品视频在线观看视频| 成人一区二区视频在线观看| 一进一出好大好爽视频| 成熟少妇高潮喷水视频| 自拍偷自拍亚洲精品老妇| 能在线免费观看的黄片| 日本-黄色视频高清免费观看| 亚洲五月天丁香| 亚洲无线观看免费| 超碰av人人做人人爽久久| 在线国产一区二区在线| 一级毛片电影观看 | 国产av在哪里看| 国产精品国产三级国产av玫瑰| 九九爱精品视频在线观看| 3wmmmm亚洲av在线观看| 国产v大片淫在线免费观看| 亚洲欧美成人综合另类久久久 | 婷婷色综合大香蕉| 亚洲av不卡在线观看| 在线观看一区二区三区| 欧美极品一区二区三区四区| 婷婷色综合大香蕉| 精品日产1卡2卡| 久久国产乱子免费精品| 一本一本综合久久| 黄色日韩在线| 久久午夜福利片| 国产亚洲精品av在线| 久久久久久久久大av| 天堂影院成人在线观看| 久久久久久久久大av| av在线老鸭窝| 免费不卡的大黄色大毛片视频在线观看 | 国产高清视频在线观看网站| 国产色爽女视频免费观看|