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

    靜水壓下熱損傷大理巖動態(tài)拉伸特性研究

    2020-08-31 13:36:46夏開文蔡英鵬潘永慶
    關(guān)鍵詞:靜水壓

    夏開文,蔡英鵬,徐?穎,王?帥,潘永慶

    靜水壓下熱損傷大理巖動態(tài)拉伸特性研究

    夏開文1, 2,蔡英鵬1, 2,徐?穎1, 2,王?帥1, 2,潘永慶3

    (1. 天津大學(xué)水利工程仿真與安全國家重點實驗室,天津 300072;2. 天津大學(xué)建筑工程學(xué)院,天津 300350;3. 中水東北勘測設(shè)計研究有限責(zé)任公司,長春 130021)

    處于高溫高地應(yīng)力等復(fù)雜地質(zhì)環(huán)境下的深部圍巖,可能遭受爆破、地震等動態(tài)荷載影響,在巖石工程開挖區(qū)附近多導(dǎo)致張拉破壞,因此,研究熱損傷巖石在不同地應(yīng)力條件下的動態(tài)拉伸特性,在深部巖石工程中具有重要意義.選取均質(zhì)細(xì)粒房山大理巖,利用自主研發(fā)的多功能分離式霍普金森壓桿(SHPB)系統(tǒng),進(jìn)行4種溫度(25℃、250℃、450℃、700℃)損傷梯度下,4種靜水壓環(huán)境(0MPa、5MPa、10MPa和20MPa)下,房山大理巖巴西圓盤試樣動態(tài)加載試驗.研究結(jié)果表明:①隨著溫度的增加,房山大理巖內(nèi)部微裂紋增多,礦物成分由CaMg(CO3)2向CaCO3和MgO轉(zhuǎn)變,密度和波速均隨著損傷溫度的增加而逐漸降低;②靜水壓條件下試樣拉伸應(yīng)力曲線呈現(xiàn)為雙峰特征,這主要是由于側(cè)向圍壓的存在抑制了試樣的劈裂,使試樣拉伸破壞后繼續(xù)承載造成的;③在固定圍壓環(huán)境下,各溫度梯度下熱損傷大理巖的動態(tài)拉伸強度均具有明顯的率相關(guān)性;同時動態(tài)拉伸強度隨著溫度的升高而明顯降低;④在固定溫度下,圍壓的存在明顯提高了房山大理巖動態(tài)拉伸強度,但拉伸強度的增幅隨著圍壓的增加而減弱.此外,當(dāng)試樣的熱處理溫度超過450℃后,動態(tài)拉伸強度的圍壓效應(yīng)小于25℃和250℃熱處理情況.這可能與高溫處理后的巖石在高圍壓狀態(tài)下發(fā)生了脆性向延性轉(zhuǎn)變有關(guān).

    深部巖石動力學(xué);熱損傷;動態(tài)拉伸強度;分離式霍普金森壓桿;房山大理巖

    巖石是構(gòu)成地球的基本材料,地球上99%的巖石都處于600℃和1GPa以上的高溫高壓狀態(tài)[1].此外,巖石作為一種準(zhǔn)脆性材料,其抗拉強度遠(yuǎn)遠(yuǎn)小于抗壓強度,在深部地下巖石工程中,巷道變形失穩(wěn)、巖爆等一些受動力擾動后地質(zhì)災(zāi)害的發(fā)生往往由巖石的張拉破壞造成.隨著深部巖石開采深度的增加,巖層地應(yīng)力提高,溫度上升,深部巖體表現(xiàn)出與淺部不同的力學(xué)特性[2-3].此外,處于高溫高地應(yīng)力等復(fù)雜地質(zhì)環(huán)境下的深部巖石工程,在受到爆破、地震等動態(tài)荷載作用后,在開挖區(qū)附近多導(dǎo)致張拉破壞.因此,對高溫高壓狀態(tài)下巖石的動態(tài)拉伸特性展開研究具有重要意義.

    目前,研究巖石拉伸特性的方法有兩種:直接拉伸方法和間接拉伸方法.由于直接拉伸相較于巴西圓盤劈裂試驗,存在操作復(fù)雜、費時費力、經(jīng)濟(jì)性差以及試驗過程中容易造成樣品一端應(yīng)力集中和破壞等缺點[4],巴西圓盤劈裂試驗逐漸成為巖石動態(tài)拉伸強度測試標(biāo)準(zhǔn)方法[5-6].鑒于巖石拉伸特性對于地下巖石工程建設(shè)具有重要意義,國內(nèi)外大量學(xué)者不僅高度關(guān)注其靜態(tài)特性,對于巖石的動態(tài)拉伸強度也逐漸引起廣泛重視[5,7-10].在溫度對巖石動態(tài)拉伸特性影響方面,劉石等[11]對高溫?fù)p傷后大理巖進(jìn)行了巴西圓盤劈裂試驗,指出高溫后大理巖的動態(tài)劈裂破壞形式受到?jīng)_擊加載速率和溫度的共同影響;黃彥華等[12]針對高溫后含預(yù)制缺陷巖石的裂紋擴(kuò)展特征進(jìn)行了研究,通過對花崗巖圓盤試樣預(yù)制孔洞并進(jìn)行巴西圓盤劈裂試驗,探討了溫度和孔洞分布對花崗巖的物理、力學(xué)特性及裂紋演化特征的影響;Yin等[13]研究了熱處理對勞倫特花崗巖物理性能和抗拉強度的影響,發(fā)現(xiàn)隨著加熱溫度升高,試樣密度、縱波波速和靜態(tài)拉伸強度均減小,而動態(tài)拉伸強度先升高后降低,并從礦物成分及微裂紋角度對上述現(xiàn)象做出了相應(yīng)解釋;Yao等[14]開展了考慮溫度和加載率等因素影響的巴西圓盤試驗,建立了關(guān)于龍游砂巖動態(tài)拉伸強度隨溫度、加載率變化的本構(gòu)模型,并進(jìn)一步利用X射線、CT技術(shù)及圖像處理技術(shù),對龍游砂巖溫度損傷程度給予跟蹤描述及量化表征.

    深部巖石工程中的圍巖不僅受到溫度的影響,更重要的是受到地應(yīng)力的作用[1,15].夏開文等[3]將深部巖石力學(xué)試驗中試樣的圍壓狀態(tài)分為4種類型:軸向圍壓、側(cè)向圍壓、靜水壓和三軸圍壓狀態(tài).在側(cè)向圍壓對巖石動態(tài)拉伸特性影響方面,袁超[16]、劉世奇[17]進(jìn)行了花崗巖試樣側(cè)向圍壓下動態(tài)直接拉伸試驗,試驗結(jié)果顯示:對于同一側(cè)向壓力,花崗巖的抗拉強度與直接拉伸彈性模量隨應(yīng)變速率增加而呈現(xiàn)增加趨勢;對于同一應(yīng)變速率,其抗拉強度與直接拉伸彈性模量隨側(cè)向壓力增加而呈現(xiàn)遞減趨勢.Wu等[15]對花崗巖巴西圓盤試樣施加靜水壓約束,探討了巖石動態(tài)拉伸強度隨靜水壓(最高20MPa)變化的演化規(guī)律. 夏開文等[18]通過具有預(yù)加載裝置的分離式霍普金森壓桿系統(tǒng)研究了巖石在不同預(yù)拉伸應(yīng)力下的拉伸強度,結(jié)果表明:動態(tài)拉伸強度和總拉伸強度隨著加載率的增加而增加,同時,在相同加載率下,動態(tài)拉伸強度隨著預(yù)拉伸載荷的增加而減小,而總拉伸應(yīng)力與預(yù)拉伸載荷的大小無關(guān).李夕兵等[2]在國內(nèi)較早開展了溫壓耦合作用下巖石的動態(tài)力學(xué)特性研究,對20~300℃和不同軸向靜壓作用下粉砂巖的動態(tài)壓縮強度進(jìn)行研究,探究了溫壓耦合作用下巖石的動態(tài)強度變化規(guī)律.

    綜上所述,眾多學(xué)者對溫度和應(yīng)力作用下巖石的動態(tài)拉伸特性進(jìn)行了大量研究.但是由于受到加載裝置的限制,在模擬地應(yīng)力條件下的動態(tài)拉伸試驗中,僅能實現(xiàn)對樣品的軸向加載[19],而難以實現(xiàn)圍壓加載,因此對于圍壓作用下熱損傷巖石的動態(tài)拉伸特性的研究依然不夠系統(tǒng).而且,當(dāng)前的研究多集中于溫度或軸向應(yīng)力單方面對巖石動態(tài)拉伸特性影響,缺乏綜合考慮二者影響的研究.深部巖石工程中圍巖尤其是處于深部硐室遠(yuǎn)場區(qū)的巖石多處于高溫高地應(yīng)力狀態(tài),因此有必要對圍壓狀態(tài)下熱損傷后巖石的動態(tài)拉伸特性展開系統(tǒng)研究.筆者通過自主研發(fā)的帶圍壓功能的分離式霍普金森壓桿(split Hopkinson pressure bar,SHPB)系統(tǒng),對不同靜水壓條件下熱損傷大理巖的動態(tài)拉伸特性進(jìn)行系統(tǒng)研究,探討不同溫度損傷、不同地應(yīng)力條件下大理巖動態(tài)拉伸破壞演化規(guī)律,為深入開展深部巖石動力學(xué)的研究提供一定理論依據(jù).

    1?試樣制備

    1.1?試樣的選取與制備

    試驗樣品為細(xì)顆粒大理巖,取自北京房山地區(qū),俗稱漢白玉.不僅在深部地下工程如隧道開挖、礦冶爆破中會經(jīng)常遇到大理巖,且由于該種巖石其均質(zhì)性和完整性較好,還廣泛用于建筑物的基礎(chǔ),如歷史建筑和橋梁裝飾.目前,國內(nèi)外學(xué)者對房山大理巖進(jìn)行了大量相關(guān)研究,其主要礦物成分和物理力學(xué)性質(zhì)相對明確[20-22].本試驗選取均勻完整的房山大理巖石塊,參照國際巖石力學(xué)與工程學(xué)會(ISRM)推薦方法[5]以及中國爆破行業(yè)協(xié)會標(biāo)準(zhǔn)《巖石材料巴西圓盤試樣動態(tài)拉伸強度測試方法》(T/CSEB 0001—2018)[6]的相關(guān)要求,將房山大理巖石塊制成48mm×24mm圓盤狀試樣,并對所有試樣端面進(jìn)行剖光.加工后大理巖試樣如圖1所示,試樣的基本物理力學(xué)參數(shù)見表1[20].巖石樣品的水分來源于吸附空氣中游離的水分子.由于巖石從完整均質(zhì)巖塊上取樣,除溫度處理不同外,其他加工、制備、儲藏條件均一致,故認(rèn)為各級溫度處理試樣初始含水率一致.

    圖1?典型大理巖圓盤試樣(450℃)

    表1?大理巖基本物理力學(xué)參數(shù)

    Tab.1?Basic physical and mechanical parameters of marble

    1.2?試樣熱處理

    依據(jù)前人對大理巖溫度損傷研究成果及大理巖經(jīng)歷不同溫度損傷的變化情況,選取了4個級別溫度(25℃、250℃、450℃、700℃)對試樣進(jìn)行熱處理.溫度加載試驗在天津大學(xué)巖石力學(xué)實驗室GWL 1200℃高溫電爐中完成.為避免熱沖擊造成的不均勻損傷,將加熱速度控制在2℃/min.當(dāng)爐溫達(dá)到設(shè)定溫度后,保持恒溫120min,再自然冷卻至室溫.

    大理巖在不同溫度處理階段具有不同的熱損傷過程[20].利用天津大學(xué)巖石力學(xué)實驗室自備的電子天平、游標(biāo)卡尺及聲波測速儀,分別測量各級溫度損傷后試樣的質(zhì)量、體積和縱波波速,得到房山大理巖密度和波速隨溫度損傷變化曲線,如圖2所示.可以看出,隨著處理溫度的提高,大理巖密度和波速均呈下降趨勢.試樣經(jīng)過250℃高溫后,密度有微弱下降,幅度為0.148%,主要是因礦物內(nèi)部的結(jié)晶水析出所致;經(jīng)過450℃溫度損傷后,密度進(jìn)一步降低,降低幅度為0.633%,相較于250℃下降幅度有所增加,分析認(rèn)為,大理巖在450℃經(jīng)歷了內(nèi)部各種礦物的不均勻膨脹,導(dǎo)致巖石內(nèi)部微裂紋增加,體積有所增大,密度降低.當(dāng)熱處理溫度為700℃時,密度大幅降低,幅度高達(dá)3.944%.

    波速隨著溫度的下降趨勢則與密度相反,250℃時,波速大幅降低,幅度為52.6%;當(dāng)處理溫度超過450℃時,波速降低幅度逐漸減小.各級溫度荷載下的波速相較于常溫狀態(tài),分別下降52.6%、69.8%和79.4%.

    由圖2可知,當(dāng)熱損傷溫度小于等于450℃時,巖石的縱波波速下降較密度明顯,分析認(rèn)為上述損傷階段主要為試樣內(nèi)部微孔隙的增加,而縱波波速相較于密度對微孔隙增加反應(yīng)更加敏感,故縱波波速比密度下降更加明顯.而試樣經(jīng)過700℃溫度損傷后,密度急劇降低,而波速降低幅度有限.這是由于試樣內(nèi)部礦物成分發(fā)生化學(xué)反應(yīng),比如CO2的釋放而導(dǎo)致樣品質(zhì)量減小,同時樣品的體積基本保持不變,所以這種改變對密度影響更為顯著.

    圖2?不同溫度下大理巖密度與縱波波速曲線

    通過XRD物相分析不同溫度下試樣的礦物成分變化.常溫條件下該大理巖試樣的主要礦物成分為白云石(97%),含少量石英(3%),檢測結(jié)果如圖3(a)所示.圖3(b)則展示了不同溫度損傷后的大理巖XRD圖譜,可以看出隨著處理溫度的增加,試樣中主要礦物CaMg(CO3)2發(fā)生分解,含量逐漸降低,說明白云石含量逐漸降低;當(dāng)處理溫度達(dá)到700℃時,試樣內(nèi)部主要礦物白云石的特征峰消失,其主要礦物成分為MgO和CaCO3,這是由于CaMg(CO3)2在700℃時發(fā)生分解產(chǎn)生CO2所致,這也是導(dǎo)致試樣密度明顯降低的原因[23-24].

    圖3?大理巖XRD分析結(jié)果

    2?試驗裝置與方案

    2.1?試驗裝置

    作為國際巖石力學(xué)學(xué)會推薦的巖石動力學(xué)標(biāo)準(zhǔn)測試裝置,分離式霍普金森壓桿(SHPB)已成為測試巖石各種動態(tài)力學(xué)參數(shù)的主要工具[5-6].試驗采用天津大學(xué)巖土工程研究所自主研發(fā)的50mm多功能分離式霍普金森壓桿系統(tǒng)進(jìn)行動態(tài)加載試驗[25],該系統(tǒng)由發(fā)射系統(tǒng)、桿件系統(tǒng)、圍壓裝置、軸壓裝置及信號記錄系統(tǒng)組成[26],試驗裝置如圖4所示.該系統(tǒng)所有桿件材質(zhì)均為馬氏體鋼,彈性模量211GPa,縱波波速5280m/s.入射桿、透射桿分別長3.0m和1.8m,吸收桿長0.5m,子彈長0.35m.圍壓裝置由圍壓缸、軸壓缸和液壓油泵組成,圍壓缸最大工作壓力60MPa[26].

    圖4?SHPB系統(tǒng)示意

    在試驗過程中,改變子彈發(fā)射所需的氣壓實現(xiàn)對樣品不同沖擊速度的動態(tài)加載;通過改變軸壓缸和圍壓缸的油壓值,實現(xiàn)對樣品不同靜水壓狀態(tài)的加載.為保證樣品能夠維持在入射桿和透射桿之間,軸壓缸的加載要略高于圍壓缸的加載值.信號記錄系統(tǒng)包括粘貼在透射桿、入射桿以及試樣端面上的應(yīng)變片、應(yīng)變放大器以及可以實時同步記錄儲存信號的示波器.

    為避免空氣中的水分和圍壓缸內(nèi)的液壓油滲入試樣孔隙,影響最終測試結(jié)果,動態(tài)試驗前將所有試樣進(jìn)行蠟封處理,然后安放于圍壓缸中入射桿和透射桿之間.

    2.2?動態(tài)力平衡及數(shù)據(jù)處理

    為忽略動態(tài)測試中的慣性效應(yīng),首先要使試樣兩端實現(xiàn)動態(tài)力平衡.采用直徑9.0mm、厚1.5mm的C11000銅片將入射波形狀由矩形波變?yōu)檎也? Frew等[27]對SHPB壓縮試驗中的脈沖整形技術(shù)有過詳細(xì)討論.圖5為典型桿端力平衡曲線,其中圖5(a)為圍壓為0MPa時桿端力平衡曲線,可見試樣兩端力平衡較好,可忽略其慣性效應(yīng).

    圖5(b)為圍壓為20MPa時桿端力平衡曲線,與圖5(a)不同的是,入射波與透射波起點均高于0kN.這是由于在對試樣進(jìn)行施加靜水壓時,入射桿和透射桿均處于預(yù)壓狀態(tài),且預(yù)加載力值等于設(shè)定的靜水壓力值對應(yīng)的加載力,例如:20MPa對應(yīng)的桿端預(yù)加載力為39.27kN(圖5(b)中虛線框),當(dāng)采用直流耦合(DC)進(jìn)行應(yīng)變數(shù)據(jù)采集時,二者的波形起點則為桿中的預(yù)應(yīng)力值[28].

    基于一維應(yīng)力波理論,試樣與入射桿接觸面的動態(tài)力1和試樣與透射桿接觸面處的動態(tài)力2為

    式中:為桿的彈性模量;為桿截面面積;i()、r()和t()分別為入射波、反射波、透射波的應(yīng)??變[4-6]. Chen 等[28]的理論研究結(jié)果顯示,發(fā)射桿速度和預(yù)加荷載對圍壓狀態(tài)下力平衡的計算公式有影響,在本套試驗裝置條件下,預(yù)加載應(yīng)力小于60MPa時,力平衡條件為

    因此上述靜水壓條件下力平衡的條件為1=2.

    此外,圖5(b)中圍壓條件下透射桿信號是在DC模式下采集獲得,其起始值為透射桿中預(yù)加靜水壓力值.根據(jù)Wu等[15]對巴西圓盤試樣受力分析可知,靜水壓下巴西圓盤試樣動態(tài)拉伸應(yīng)力與預(yù)加靜水壓無關(guān),即在進(jìn)行動態(tài)拉伸強度取值時,應(yīng)減去預(yù)加載荷.值得一提的是,若使用的是交流耦合(AC),透射波信號則直接從0起跳,無需修正.

    因此,當(dāng)試樣兩端達(dá)到力平衡,試樣中慣性效應(yīng)可忽略時,其動態(tài)拉伸應(yīng)力可由靜態(tài)公式求得,即

    式中:為時間;()為試樣動態(tài)拉伸應(yīng)力;pre為預(yù)軸向力,其中當(dāng)無圍壓狀態(tài)時,pre=0;為巴西圓盤樣品的直徑;為樣品的厚度.

    由上述試驗步驟以及式(1)~(4),可得各溫度及靜水壓條件下試樣的動態(tài)拉伸應(yīng)力隨時間的變化曲線.不同靜水壓下典型拉伸應(yīng)力-時間曲線見圖6.加載率由峰值前線性部分的斜率確定(如圖中紅色虛線所示)[5-6].

    圖6?典型拉伸應(yīng)力-時間曲線

    3?試驗結(jié)果分析

    3.1?試樣破壞模式

    在充分考慮地應(yīng)力對應(yīng)的埋深深度以及試驗數(shù)據(jù)有效性,將該試驗中靜水壓的狀態(tài)設(shè)定為4個,分別為0MPa、5MPa、10MPa、20MPa,對應(yīng)的埋深深度約為0m、175m、350m、700m.經(jīng)過溫度預(yù)處理的圓盤試樣,分別進(jìn)行上述4個圍壓條件下的動態(tài)拉伸加載試驗.圖6(a)為無圍壓下動態(tài)拉伸應(yīng)力-時間曲線和破壞試樣.試樣的破壞模式是巖石破壞機理的直觀指標(biāo)[29].該試樣在受軸向動態(tài)沖擊加載后,其圓盤試樣均沿直徑方向劈裂成2個半圓,其劈裂面平行于加載力方向,且無明顯端部壓裂破壞,如圖6(a)所示.同時其動態(tài)拉伸應(yīng)力曲線呈現(xiàn)出單峰特征,與常規(guī)動態(tài)拉伸應(yīng)力-時間曲線形態(tài)一致[30].該峰值應(yīng)力即為動態(tài)巴西圓盤拉伸強度.

    與圖6(a)不同的是,靜水壓條件下動態(tài)拉伸應(yīng)力曲線出現(xiàn)2個峰值.在圍壓狀態(tài)下,當(dāng)試樣受到沿入射桿軸向方向的動態(tài)加載后,試樣內(nèi)部首先產(chǎn)生微裂紋,并沿著桿軸向擴(kuò)張直至破壞,從而形成第1個峰值;由于動態(tài)加載時間極短,試樣周圍的液壓油無法瞬時卸載,對試樣裂紋擴(kuò)張及劈裂存在抑制作??用[15].通過圖6(b)中試樣的破碎形態(tài)可知,試樣除了劈裂成兩部分外,其與桿接觸的兩端有壓裂區(qū)域.這是由于試樣在受拉劈裂后,在靜水壓的作用下繼續(xù)承受荷載,直至發(fā)生壓縮破壞,因此出現(xiàn)第2個波峰.可見,靜水壓條件下第1個峰值對應(yīng)的拉應(yīng)力為動態(tài)拉伸強度.

    3.2?基于溫度-圍壓-加載率的動態(tài)拉伸強度力學(xué)模型

    為合理量化描述圍壓、溫度和加載率對大理巖動態(tài)拉伸強度的影響,客觀分析強度與圍壓、溫度及加載率之間的變化規(guī)律,基于圖7和圖8所示的動態(tài)拉伸強度與圍壓、溫度和加載率之間的關(guān)系,提出了如下公式來擬合試驗結(jié)果,并進(jìn)一步解釋巖石在熱損傷后靜水約束條件下的動態(tài)拉伸行為.

    3.3?動態(tài)拉伸強度熱損傷效應(yīng)

    根據(jù)上述方法獲得的各靜水壓條件下房山大理巖的動態(tài)拉伸強度見圖7.由圖7可知,固定靜水壓條件下,各溫度荷載下試樣的動態(tài)拉伸強度均具有明顯的率相關(guān)性,且隨著溫度的增加,拉伸強度整體呈下降趨勢.例如:由圖7(a)中擬合曲線可知,當(dāng)加載率為800GPa/s、圍壓為0MPa時,25℃、250℃、450℃和700℃溫度荷載下的動態(tài)強度分別為30.43MPa、28.44MPa、26.66MPa和24.50MPa,與常溫狀態(tài)下大理巖的動態(tài)拉伸強度相比,250℃、450℃和700℃下的強度分別下降6.54%、12.39%和19.65%.隨著損傷溫度的增加,巖石動態(tài)拉伸強度的降低幅度也逐漸增加.

    分析認(rèn)為,在250℃溫度損傷條件下,試樣波速出現(xiàn)大幅降低,但密度僅出現(xiàn)微弱下降,表明造成密度降低的原因是由于試樣中缺陷增加,但是缺陷多表現(xiàn)為沒有體積的微裂紋.微裂紋的出現(xiàn)是由于礦物的不均勻膨脹,使得試樣內(nèi)部開始出現(xiàn)微裂紋,微裂紋的出現(xiàn)也導(dǎo)致了拉伸強度的降低.

    當(dāng)溫度荷載升高到450℃后,大理巖的動態(tài)拉伸強度進(jìn)一步下降,這也是由于巖石的不均勻膨脹造成的.由于試樣微裂紋的出現(xiàn)和發(fā)展,使得動態(tài)拉伸強度進(jìn)一步弱化.

    當(dāng)溫度荷載升至700℃時,試樣內(nèi)部主要礦物成分由單一CaMg(CO3)2向兩種主要礦物CaCO3和活性MgO轉(zhuǎn)變[24].恒溫保存2h后,仍會存有少量未分解的白云石,在3種礦物并存下,礦物顆粒連接相比于450℃時更加弱化,巖石總孔隙增加,促使微裂紋進(jìn)一步擴(kuò)展,這由試樣密度明顯下降和波速進(jìn)一步下降可知(圖2).由于礦物成分的改變和巖石內(nèi)部孔隙裂紋的增加,最終導(dǎo)致試樣動態(tài)拉伸強度在700℃時弱化效應(yīng)更為顯著[18,24].

    3.4?動態(tài)拉伸強度靜水壓效應(yīng)

    為進(jìn)一步探究靜水壓條件對不同溫度損傷后試樣的動態(tài)拉伸強度影響規(guī)律,將上述試驗結(jié)果整理成固定溫度下,動態(tài)拉伸強度與靜水壓和加載率的關(guān)系,如圖8所示.圖8(a)~(d)對應(yīng)的固定溫度分別是常溫條件下(25℃)、250℃、450℃和700℃.可見固定溫度下,隨著靜水壓的增加,拉伸強度呈整體增加趨勢.例如:由圖8(a)中擬合曲線可知,在常溫條件、加載率為1000GPa/s時,0MPa、5MPa、10MPa和20MPa靜水壓對應(yīng)的動態(tài)拉伸強度分別為33.12MPa、44.99MPa、49.00MPa和53.28MPa.與0MPa相比,5MPa、10MPa和20MPa下的動態(tài)拉伸強度分別提高35.81%、47.92%和60.87%,可知在固定溫度荷載下,靜水壓越高,試樣的動態(tài)拉伸強度越大,說明圍壓的存在對提高巖石拉伸強度作用明顯;圖8(b)、(c)和(d)均有相似的趨勢.這一趨勢與靜水壓下花崗巖試驗動態(tài)拉伸強度一致[15].這是由于圍壓的存在不僅對巖石初始裂隙產(chǎn)生壓實作用,同時在巖石受動態(tài)加載時對裂紋的擴(kuò)展起到抑制作用,從而提升了試樣的拉伸強度.

    同時,由圖8(a)可以看出,當(dāng)加載率為1000GPa/s時,常溫試樣在靜水壓由0MPa升至5MPa時,相應(yīng)的動態(tài)拉伸強度增幅最大,為35.81%,當(dāng)靜水壓增至10MPa和20MPa時,相應(yīng)增幅逐漸減小,分別為8.92%和8.73%.可見,隨著圍壓的進(jìn)一步增加,圍壓對試樣拉伸強度的弱化效應(yīng)會逐漸減弱.該現(xiàn)象在靜水壓下花崗巖的動態(tài)拉伸強度中也被發(fā)現(xiàn)[15].

    這種高應(yīng)力條件下圍壓效應(yīng)逐漸減弱的現(xiàn)象在高溫處理的大理巖試樣中同樣存在,且較低溫條件下明顯.根據(jù)各溫度下單獨擬合結(jié)果,當(dāng)熱處理溫度為450℃時,靜水壓由5MPa提升到20MPa時,動態(tài)拉伸強度增加幅度(16.8%)小于常溫和250℃熱處理情況(15.7%).這可能是由于高溫處理后的巖石,在高圍壓下巖石材料由脆性逐漸向延性轉(zhuǎn)變造成的[31].由此可以推斷,當(dāng)達(dá)到一定深部時,巖石的拉伸強度不會隨著地應(yīng)力的增加而明顯增加,但仍會受加載率的顯著影響.同時,初始損傷大的巖石,強度受圍壓作用而增強的效應(yīng)會有所減弱.

    4?結(jié)?論

    通過采用改進(jìn)的帶圍壓SHPB裝置對4種溫度(25℃、250℃、450℃、700℃)損傷后的大理巖進(jìn)行了常規(guī)物理力學(xué)性質(zhì)測試、XRD測試及不同靜水壓下動態(tài)拉伸試驗,主要結(jié)論如下.

    (1)房山大理巖在經(jīng)歷各級溫度荷載后,其密度在450℃之后急劇下降,縱波波速在250℃下降明顯,250℃之后下降較平緩.700℃后主要礦物成分由原來的碳酸鈣鎂(CaMg(CO3)2)轉(zhuǎn)換成碳酸鈣(CaCO3)和氧化鎂(MgO).

    (2)在本試驗條件下,圍壓影響作用下的試樣拉伸應(yīng)力曲線表現(xiàn)出雙峰特征,這主要是由于側(cè)向圍壓的存在抑制了試樣的劈裂,使試樣拉伸破壞后繼續(xù)承載造成的,第1個峰值對應(yīng)的拉伸應(yīng)力為動態(tài)拉伸?強度.

    (3)熱損傷房山大理巖動態(tài)拉伸試驗結(jié)果表明,受溫度損傷后的房山大理巖,在不同靜水壓下動態(tài)拉伸強度均具有明顯的率相關(guān)性.溫度對大理巖具有明顯的熱損傷效應(yīng),在同一圍壓和同一加載率情況下,隨著熱處理溫度的增加,動態(tài)拉伸強度出現(xiàn)不同程度的弱化現(xiàn)象,當(dāng)熱處理溫度超過450℃時,拉伸強度弱化現(xiàn)象更為明顯.

    (4)在固定溫度下,靜水壓的存在對房山大理巖的動態(tài)拉伸強度有顯著的提升作用.然而,動態(tài)拉伸強度的增幅隨著圍壓的增加而有所減弱,即隨著地下工程埋深的增加,巖石的動態(tài)拉伸強度將不再顯著增加,但仍存在明顯的率相關(guān)性.當(dāng)試樣的熱處理溫度為450℃和700℃時,動態(tài)拉伸強度的圍壓效應(yīng)小于常溫(25℃)和250℃熱處理情況.這可能與高溫處理后的巖石在高圍壓狀態(tài)下發(fā)生了脆性向延性轉(zhuǎn)變有關(guān).

    [1] 陳?颙,黃庭芳,劉恩儒. 巖石物理學(xué)[M]. 合肥:中國科學(xué)技術(shù)大學(xué)出版社,2009.

    Chen Yong,Huang Tingfang,Liu Enru. Rock Physics[M]. Hefei:University of Science and Technology of China Press,2009(in Chinese).

    [2] 李夕兵,尹土兵,周子龍,等. 溫壓耦合作用下的粉砂巖動態(tài)力學(xué)特性試驗研究[J]. 巖石力學(xué)與工程學(xué)報,2010,29(12):2377-2384.

    Li Xibing,Yin Tubing,Zhou Zilong,et al. Study of dynamic properties of siltstone under coupling effects of temperature and pressure[J]. Chinese Journal of Rock Mechanics and Engineering,2010,29(12):2377-2384(in Chinese).

    [3] 夏開文,徐?穎,陳?榮. 考慮深部賦存條件的巖石動態(tài)破壞試驗研究進(jìn)展[J]. 隧道與地下工程災(zāi)害防治,2019,1(1):58-75.

    Xia Kaiwen,Xu Ying,Chen Rong. Dynamic tests of rocks subjected to simulated deep underground environments[J]. Hazard Control in Tunnelling and Underground Engineering,2019,1(1):58-75(in Chinese).

    [4] 夏開文,周傳波. 巖石動態(tài)力學(xué)參數(shù)測試綜述[J]. 工程爆破,2014,20(2):43-53.

    Xia Kaiwen,Zhou Chuanbo. Review of testing methods for dynamic rock mechanical properties[J]. Engineering Blasting,2014,20(2):43-53(in Chinese).

    [5] Zhou Y X,Xia K,Li X B,et al. Suggested methods for determining the dynamic strength parameters and mode-I fracture toughness of rock materials[J]. International Journal of Rock Mechanics and Mining Sciences,2012,49:105-112.

    [6] 中國爆破行業(yè)協(xié)會. T/CSEB 0001—2018 巖石材料巴西圓盤試樣動態(tài)拉伸強度測試方法[S]. 北京:冶金工業(yè)出版社,2018.

    China Society of Explosives and Blasting. T/CSEB 0001—2018 Testing Method for Determining Dynamic Tensile Strength of Rock Materials by the Brazilian Disc Specimen[S]. Beijing:Metallurgy Industry Press,2018(in Chinese).

    [7] Huang S,Xia K,Yan F,et al. An experimental study of the rate dependence of tensile strength softening of Longyou sandstone[J]. Rock Mechanics and Rock Engineering,2010,43(6):677-683.

    [8] Dai F,Huang S,Xia K,et al. Some fundamental issues in dynamic compression and tension tests of rocks using split Hopkinson pressure bar[J]. Rock Mechanics and Rock Engineering,2010,43(6):657-666.

    [9] Zhang Q B,Zhao J. A review of dynamic experimental techniques and mechanical behaviour of rock materials[J]. Rock Mechanics and Rock Engineering,2014,47(4):1411-1478.

    [10] Xia K,Yao W. Dynamic rock tests using split Hopkinson(Kolsky)bar system—A review[J]. Journal of Rock Mechanics and Geotechnical Engineering,2015,7(1):27-59.

    [11] 劉?石,許金余,白二雷,等. 高溫后大理巖動態(tài)劈裂拉伸試驗研究[J]. 巖土力學(xué),2010,34(12):3500-3504.

    Liu Shi,Xu Jinyu,Bai Erlei,et al. Experimental study of dynamic tensile behaviors of marble after high temperature[J]. Rock and Soil Mechanics,2010,34(12):3500-3504(in Chinese).

    [12] 黃彥華,楊圣奇. 高溫后含孔花崗巖拉伸力學(xué)特性試驗研究[J]. 中國礦業(yè)大學(xué)學(xué)報,2017,46(4):783-791.

    Huang Yanhua,Yang Shengqi. The tensile mechanical behavior of granite containing pre-existing holes after high temperature treatment[J]. Journal of China University of Mining & Technology,2017,46(4):783-791(in Chinese).

    [13] Yin T B,Li X B,Cao W Z,et al. Effects of thermal treatment on tensile strength of Laurentian granite using Brazilian test[J]. Rock Mechanics and Rock Engineering,2015,48(6):2213-2223.

    [14] Yao W,Xu Y,Wang W,et al. Dependence of dynamic tensile strength of Longyou sandstone on heat-treatment temperature and loading rate[J]. Rock Mechanics and Rock Engineering,2015,48(10):3899-3915.

    [15] Wu B,Yao W,Xia K. An experimental study of dynamic tensile failure of rocks subjected to hydrostatic confinement[J]. Rock Mechanics and Rock Engineering,2016,49(10):3855-3864.

    [16] 袁?超. 有側(cè)壓力巖石動態(tài)直接拉伸試驗及動態(tài)斷裂力學(xué)研究[D]. 湘潭:湖南科技大學(xué),2013.

    Yuan Chao. The Experiment of Dynamic Direct Tensile Test for Rock Under Lateral Pressure and Dynamic Fracture Mechanics Study[D]. Xiangtan:Hunan University of Science and Technology,2013(in Chinese).

    [17] 劉世奇. 有側(cè)向壓力的花崗巖動態(tài)直接拉伸力學(xué)特性研究[D]. 武漢:中國科學(xué)院武漢巖土力學(xué)研究所,2007.

    Liu Shiqi. Studies on Dynamic Tensile Properties of Granite Under Confining Pressure[D]. Wuhan:Institute of Rock and Soil Mechanics,Chinese Academy of Sciences,2007(in Chinese).

    [18] 夏開文,姚?偉. 預(yù)加載下巖石的動態(tài)力學(xué)性能研究[J]. 工程爆破,2015(6):7-13.

    Xia Kaiwen,Yao Wei. Dynamic mechanical properties of rock under pre-load[J]. Engineering Blasting,2015(6):7-13(in Chinese).

    [19] Gong F Q,Si X F,Li X B,et al. Dynamic triaxial compression tests on sandstone at high strain rates and low confining pressures with split Hopkinson pressure bar[J]. International Journal of Rock Mechanics and Mining Sciences,2019,113:211-219.

    [20] 高龍山,徐?穎,吳幫標(biāo),等. 溫度損傷大理巖不同含水條件下的動態(tài)壓縮特性研究[J]. 巖石力學(xué)與工程學(xué)報,2018,37(增2):3826-3833.

    Gao Longshan,Xu Ying,Wu Bangbiao,et al. Dynamic compression strength of thermal damaged Fangshan marble on dry and saturated conditions[J]. Chinese Journal of Rock Mechanics and Engineering,2018,37(Suppl 2):3826-3833(in Chinese).

    [21] 宋小林,謝和平,王啟智. 大理巖動態(tài)破裂試樣破壞應(yīng)變[J]. 巖石力學(xué)與工程學(xué)報,2005,24(16):2953-2959.

    Song Xiaolin,Xie Heping,Wang Qizhi. Failure strain of Brazilian disc samples of marble under dynamic split tests[J]. Chinese Journal of Rock Mechanics and Engineering,2005,24(16):2953-2959(in Chinese).

    [22] Yao W,He T,Xia K. Dynamic mechanical behaviors of Fangshan marble[J]. Journal of Rock Mechanics and Geotechnical Engineering,2017,9(5):807-817.

    [23] 陳國飛,楊圣奇. 高溫后大理巖力學(xué)性質(zhì)及其破壞規(guī)律研究[J]. 工程力學(xué),2014,31(8):197-204.

    Chen Guofei,Yang Shengqi. Study on failure mechanical behavior of marble after high temperature[J]. Engineering Mechanics,2014,31(8):197-204(in Chinese).

    [24] 王?沖,劉勇兵,曹占義. 白云石煅燒組織的轉(zhuǎn)變過程[J]. 材料熱處理學(xué)報,2013,34(2):25-28.

    Wang Chong,Liu Yongbing,Cao Zhanyi. Structure transformation of calcining dolomite[J]. Transactions of Materials and Heat Treatment,2013,34(2):25-28(in Chinese).

    [25] 王奇智,夏開文,吳幫標(biāo),等. 預(yù)制平行雙節(jié)理類巖石材料板動態(tài)破壞試驗研究[J]. 天津大學(xué)學(xué)報:自然科學(xué)與工程技術(shù)版,2019,52(10):1099-1108.

    Wang Qizhi,Xia Kaiwen,Wu Bangbiao,et al. Dynamic failure of simulated rock mass plate containing two parallel cracks[J]. Journal of Tianjin University:Science and Technology,2019,52(10):1099-1108 (in Chinese).

    [26] 徐?穎,夏開文,王?帥,等. 一種帶吸能裝置可模擬原位應(yīng)力環(huán)境的分離式霍普金森壓桿:中國,CN110057696A [P]. 2019-07-26.

    Xu Ying,Xia Kaiwen,Wang Shuai,et al. A Split Hopkinson Press Bar with Energy Absorption Devices That Can Be Analogous to In-Situ Stress Environment:CN,CN110057696A[P]. 2019-07-26(in Chinese).

    [27] Frew D J,F(xiàn)orrestal M J,Chen W. Pulse shaping techniques for testing brittle materials with a split Hopkinson pressure bar[J]. Experimental Mechanics,2002,42(1):93-106.

    [28] Chen R,Yao W,Lu F,et al. Evaluation of the stress equilibrium condition in axially constrained triaxial SHPB tests[J]. Experimental Mechanics,2018,58(3):527-531.

    [29] Xu Y,Yao W,Xia K,et al. Experimental study of the dynamic shear response of rocks using a modified punch shear method[J]. Rock Mechanics and Rock Engineering,2019,52:2523-2534.

    [30] Xia K,Yao W,Wu B. Dynamic rock tensile strengths of Laurentian granite:Experimental observation and micromechanical model[J]. Journal of Rock Mechanics and Geotechnical Engineering,2017,9(1):116-124.

    [31] 張強星,劉建鋒,曾?寅,等. 溫度-應(yīng)力耦合條件下膏巖滲透性試驗研究[J]. 巖石力學(xué)與工程學(xué)報,2019,38(9):1-9.

    Zhang Qiangxing,Liu Jianfeng,Zeng Yin,et al. Experimental investigation on permeability of gypsum rock under temperature-stress coupling[J]. Chinese Journal of Rock Mechanics and Engineering,2019,38(9):1-9(in Chinese).

    Experimental Study of Dynamic Tensile Failure of Thermally Damaged Fangshan Marble Subjected to Hydrostatic Confinement

    Xia Kaiwen1, 2,Cai Yingpeng1, 2,Xu Ying1, 2,Wang Shuai1, 2,Pan Yongqing3

    (1. State Key Laboratory of Hydraulic Engineering Simulation and Safety,Tianjin University,Tianjin 300072,China;2. School of Civil Engineering,Tianjin University,Tianjin 300350,China;3. China Water Northeastern Investigation,Design and Research Corporation Limited,Changchun 130021,China)

    The environment of rocks in deep engineering applications is characterized by high temperatures and high in-situ stresses. During excavation,these rocks are affected by blasting,excavation unloading,earthquakes,and other dynamic loads. Therefore,understanding the dynamic tensile characteristics of thermally damaged rocks under in-situ stress conditions is crucial. In this study,Brazilian disk(BD)specimens made of homogeneous fine-grained Fangshan marble(FM) were heat-treated at various temperatures(25℃,250℃,450℃,and 700℃). Subsequently,a modified split Hopkinson pressure bar(SHPB)system with hydrostatic confinement was utilized to conduct dynamic tests on the BD specimens. Four groups of thermally damaged samples under hydrostatic confinements of 0MPa,5MPa,10MPa,and 20MPa were tested at different loading rates. The results showed the following. ①With an increase in the heat-treatment temperature,the internal microcracks in the marble increased and the mineral composition changed from CaMg(CO3)2to CaCO3and MgO. The density and longitudinal wave velocity values of the specimens decreased owing to increases in the heat-treatment temperature. ②The tensile stress-time curves of tests with confining pressure had two peaks,whereas those of tests without confinement had only one peak. This difference is mainly caused by the existence of lateral confining pressure,which suppresses the occurrence of splitting in the specimens such that the specimens continue to bear load even after tensile fracture. ③Under a given confining pressure,the dynamic tensile strength of the thermally damaged FM specimens increased with the loading rate. Further,the dynamic tensile strength decreased with increases in the treatment temperature at a given hydrostatic confining pressure and loading rate. ④At a given treatment temperature,the dynamic tensile strength increased with the confining pressure. However,the increment of the tensile strength decreased with the hydrostatic confinement. In addition,when the heat-treatment temperature of the sample exceeded 450℃,the confining pressure effect of the dynamic tensile strength was clearly lower than those at 25℃ and 250℃. This may be due to the brittle-to-ductile transition of the rock under high confining pressure after high-temperature treatment. These results indicate the dynamic tensile failure characteristics of thermally damaged rocks under hydrostatic confinement and provide theoretical guidance for ensuring the safety of deep-rock engineering construction and disaster prevention.

    deep rock dynamics;thermal damage;dynamic tensile strength;split Hopkinson pressure bar (SHPB);Fangshan marble

    TU45

    A

    0493-2137(2020)11-1136-10

    10.11784/tdxbz202001002

    2020-01-02;

    2020-03-23.

    夏開文(1973—??),男,博士,教授,kaiwen@tju.edu.cn.

    徐?穎,maggie_xu@tju.edu.cn.

    國家自然科學(xué)基金青年基金資助項目(51704211);國家自然科學(xué)基金資助項目(51879184).

    Supported by the National Natural Science Foundation for Young Scientists of China(No. 51704211),the National Natural Science Founda-tion of China(No. 51879184).

    (責(zé)任編輯:金順愛)

    猜你喜歡
    靜水壓
    洗滌對防水織物面料靜水壓性能的影響
    中國纖檢(2022年8期)2022-09-22 07:28:06
    SDL Atlas推出第二代HydroPro耐靜水壓測試儀
    靜水壓測試儀校準(zhǔn)方法探討
    身痛逐瘀湯對人髓核細(xì)胞模型PI3K/Akt信號通路Bad、Caspase-9、GSK-3β表達(dá)的影響
    沖鋒衣面料防水透濕指標(biāo)分析與比較
    深海環(huán)境中靜水壓對2205鋼腐蝕電化學(xué)行為的影響
    海底管道修復(fù)用三通連接器的設(shè)計
    防水性紡織品耐靜水壓性能的探討
    中國纖檢(2018年5期)2018-05-10 03:16:15
    紡織品靜水壓測試標(biāo)準(zhǔn)比較分析
    中國纖檢(2017年8期)2017-12-20 21:18:34
    織物防水透濕性能測試方法及在我國產(chǎn)品標(biāo)準(zhǔn)中的應(yīng)用
    中國纖檢(2017年7期)2017-12-15 13:04:41
    国产私拍福利视频在线观看| 午夜福利成人在线免费观看| 亚洲av中文字字幕乱码综合| 夫妻性生交免费视频一级片| 2021天堂中文幕一二区在线观| 国产成年人精品一区二区| 观看美女的网站| 在线免费十八禁| videos熟女内射| 少妇熟女aⅴ在线视频| 中文天堂在线官网| 亚洲精品乱久久久久久| 一区二区三区四区激情视频| 国产黄色视频一区二区在线观看 | 亚洲国产精品专区欧美| 国模一区二区三区四区视频| 中文字幕免费在线视频6| 中文字幕av在线有码专区| 国产av码专区亚洲av| 成人三级黄色视频| 一级毛片电影观看 | 国产精品99久久久久久久久| 汤姆久久久久久久影院中文字幕 | 我要看日韩黄色一级片| 午夜爱爱视频在线播放| 午夜久久久久精精品| 天天一区二区日本电影三级| av免费观看日本| 欧美一区二区国产精品久久精品| 日本一二三区视频观看| 99热这里只有是精品50| 在线天堂最新版资源| 午夜a级毛片| 久久综合国产亚洲精品| 最近的中文字幕免费完整| 国产精品麻豆人妻色哟哟久久 | 精品久久久久久电影网 | 免费观看在线日韩| 国产欧美另类精品又又久久亚洲欧美| 精品欧美国产一区二区三| 亚洲国产精品国产精品| 九九在线视频观看精品| 亚洲一级一片aⅴ在线观看| 两个人视频免费观看高清| 国产精品一二三区在线看| 亚洲丝袜综合中文字幕| 亚洲精品自拍成人| 欧美性猛交╳xxx乱大交人| 国产黄色视频一区二区在线观看 | 一夜夜www| 国产乱来视频区| 国产激情偷乱视频一区二区| 91久久精品电影网| 97超视频在线观看视频| 久久国内精品自在自线图片| 99久久人妻综合| 搞女人的毛片| 看非洲黑人一级黄片| 青春草亚洲视频在线观看| 成人av在线播放网站| 一级爰片在线观看| 两个人的视频大全免费| 一级毛片电影观看 | 一边摸一边抽搐一进一小说| 久久久a久久爽久久v久久| 国产老妇女一区| 精品国产三级普通话版| videossex国产| 免费黄网站久久成人精品| 久久久久久久久久成人| 久久精品夜夜夜夜夜久久蜜豆| 观看美女的网站| 建设人人有责人人尽责人人享有的 | 国产精品99久久久久久久久| 少妇被粗大猛烈的视频| 校园人妻丝袜中文字幕| 国产精品蜜桃在线观看| 伦理电影大哥的女人| 少妇丰满av| 久久99热这里只频精品6学生 | 日韩欧美 国产精品| 97超视频在线观看视频| 午夜爱爱视频在线播放| 国产精品不卡视频一区二区| 国产真实伦视频高清在线观看| 乱人视频在线观看| 欧美3d第一页| 日韩成人伦理影院| 久久精品久久久久久噜噜老黄 | 日本-黄色视频高清免费观看| 亚洲综合色惰| 亚洲欧美日韩高清专用| 免费av观看视频| 国产亚洲5aaaaa淫片| 国产精品精品国产色婷婷| a级毛色黄片| 尾随美女入室| 最新中文字幕久久久久| 1024手机看黄色片| 亚洲国产日韩欧美精品在线观看| 精华霜和精华液先用哪个| 校园人妻丝袜中文字幕| 又爽又黄a免费视频| 亚洲成人av在线免费| 国产精品久久久久久精品电影| 国语对白做爰xxxⅹ性视频网站| 亚洲天堂国产精品一区在线| 精品国产露脸久久av麻豆 | 国产不卡一卡二| 亚洲国产欧美人成| 狂野欧美白嫩少妇大欣赏| 亚洲av电影不卡..在线观看| 国产熟女欧美一区二区| 久久人人爽人人爽人人片va| 麻豆精品久久久久久蜜桃| 国产成人一区二区在线| 午夜a级毛片| 丝袜美腿在线中文| 如何舔出高潮| 国产精品一区二区在线观看99 | 夜夜看夜夜爽夜夜摸| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国内揄拍国产精品人妻在线| 成人无遮挡网站| 真实男女啪啪啪动态图| 久久久久久九九精品二区国产| 又粗又硬又长又爽又黄的视频| 欧美色视频一区免费| 国产精品电影一区二区三区| 欧美日韩一区二区视频在线观看视频在线 | 国产单亲对白刺激| .国产精品久久| 色吧在线观看| 1024手机看黄色片| av在线亚洲专区| 久久国产乱子免费精品| 一本久久精品| 少妇的逼水好多| 男人舔奶头视频| av在线天堂中文字幕| 久久精品国产99精品国产亚洲性色| 黄色配什么色好看| 搡女人真爽免费视频火全软件| 亚洲国产精品成人综合色| 精品99又大又爽又粗少妇毛片| 亚洲成色77777| 免费电影在线观看免费观看| 国产淫片久久久久久久久| 欧美97在线视频| 久久久亚洲精品成人影院| 国产视频首页在线观看| 亚洲久久久久久中文字幕| 免费av不卡在线播放| 久久久久久久久久黄片| 成年免费大片在线观看| 99视频精品全部免费 在线| 久久99热这里只有精品18| 久久国内精品自在自线图片| 国产伦一二天堂av在线观看| 亚洲自偷自拍三级| 又爽又黄无遮挡网站| 久久久久久久久久久免费av| 国产一区二区在线av高清观看| 少妇熟女欧美另类| 国产黄片美女视频| 91久久精品国产一区二区三区| 女人被狂操c到高潮| 亚洲色图av天堂| 七月丁香在线播放| 免费在线观看成人毛片| 国语自产精品视频在线第100页| 日韩欧美三级三区| 又爽又黄无遮挡网站| 黄色日韩在线| 淫秽高清视频在线观看| 成年av动漫网址| 永久网站在线| 久久久久久久午夜电影| 欧美成人免费av一区二区三区| 在线观看66精品国产| 亚洲av福利一区| 最近的中文字幕免费完整| 日韩国内少妇激情av| 亚洲怡红院男人天堂| 午夜精品在线福利| 国产成人一区二区在线| 亚洲国产欧洲综合997久久,| 中文字幕制服av| 你懂的网址亚洲精品在线观看 | 国产成人freesex在线| 亚洲伊人久久精品综合 | 国产精品国产高清国产av| videossex国产| 91精品国产九色| 夫妻性生交免费视频一级片| 听说在线观看完整版免费高清| 久久久久久久国产电影| 两个人的视频大全免费| 一区二区三区四区激情视频| 联通29元200g的流量卡| 日本午夜av视频| 久久午夜福利片| 久久久精品大字幕| 亚洲乱码一区二区免费版| 少妇猛男粗大的猛烈进出视频 | 免费看日本二区| 亚洲av二区三区四区| 国产精品人妻久久久影院| 久久人人爽人人片av| 亚洲欧美日韩东京热| 国产精品久久久久久精品电影| 最近2019中文字幕mv第一页| 色网站视频免费| АⅤ资源中文在线天堂| 国产精品久久电影中文字幕| 午夜福利在线在线| 最后的刺客免费高清国语| 美女被艹到高潮喷水动态| 中文欧美无线码| 国产成人a∨麻豆精品| 精品国内亚洲2022精品成人| 亚洲成人中文字幕在线播放| 一本一本综合久久| 精品酒店卫生间| 国产国拍精品亚洲av在线观看| 建设人人有责人人尽责人人享有的 | 国产精品无大码| 一区二区三区免费毛片| 色5月婷婷丁香| 色综合亚洲欧美另类图片| 最近的中文字幕免费完整| 成年版毛片免费区| 九九热线精品视视频播放| 老司机福利观看| 中文亚洲av片在线观看爽| 中国国产av一级| 蜜桃亚洲精品一区二区三区| 日本免费在线观看一区| 最近中文字幕高清免费大全6| 一级爰片在线观看| 欧美日韩在线观看h| 18禁裸乳无遮挡免费网站照片| 久久久久久久久中文| 草草在线视频免费看| 水蜜桃什么品种好| 熟女电影av网| av视频在线观看入口| 午夜精品一区二区三区免费看| 人人妻人人澡欧美一区二区| 一级毛片我不卡| 人人妻人人看人人澡| 亚洲av成人精品一二三区| 国产片特级美女逼逼视频| 熟妇人妻久久中文字幕3abv| 中文精品一卡2卡3卡4更新| 美女高潮的动态| 亚洲电影在线观看av| 青青草视频在线视频观看| 久久久久久久亚洲中文字幕| 亚洲美女视频黄频| 欧美丝袜亚洲另类| 22中文网久久字幕| 精品久久久久久久末码| 国产免费男女视频| 国产成人福利小说| 一级爰片在线观看| 草草在线视频免费看| 久久精品夜色国产| 国产精品久久久久久久电影| 日本爱情动作片www.在线观看| 国产免费又黄又爽又色| 亚洲综合色惰| 国产精品久久电影中文字幕| 黄片无遮挡物在线观看| 中文乱码字字幕精品一区二区三区 | 色综合色国产| 激情 狠狠 欧美| 级片在线观看| 欧美日本亚洲视频在线播放| 熟女人妻精品中文字幕| 国产熟女欧美一区二区| 一二三四中文在线观看免费高清| 日日摸夜夜添夜夜添av毛片| 亚洲熟妇中文字幕五十中出| 亚洲精品一区蜜桃| 久久久午夜欧美精品| 欧美成人一区二区免费高清观看| 成人高潮视频无遮挡免费网站| 亚洲av成人av| 亚洲天堂国产精品一区在线| 国内揄拍国产精品人妻在线| 又爽又黄无遮挡网站| 国产免费一级a男人的天堂| 国产免费视频播放在线视频 | 国产成人精品久久久久久| 一个人免费在线观看电影| 午夜精品一区二区三区免费看| 久久精品国产自在天天线| 精品人妻熟女av久视频| 永久网站在线| 亚洲精品国产av成人精品| 伊人久久精品亚洲午夜| 亚洲欧洲日产国产| 九九爱精品视频在线观看| 日韩人妻高清精品专区| 亚洲综合精品二区| 国产又色又爽无遮挡免| 日韩欧美精品v在线| 99九九线精品视频在线观看视频| 中文精品一卡2卡3卡4更新| 国产麻豆成人av免费视频| 国产免费又黄又爽又色| 三级经典国产精品| 国产欧美另类精品又又久久亚洲欧美| 熟女人妻精品中文字幕| 国产成人a区在线观看| 精品熟女少妇av免费看| 高清av免费在线| 色播亚洲综合网| 秋霞伦理黄片| 91aial.com中文字幕在线观看| 色吧在线观看| 秋霞在线观看毛片| 欧美人与善性xxx| 亚洲综合色惰| 精品欧美国产一区二区三| 成人午夜精彩视频在线观看| 亚洲精品日韩av片在线观看| 最近手机中文字幕大全| 成人欧美大片| 欧美日本亚洲视频在线播放| 免费观看性生交大片5| 高清毛片免费看| 国产老妇女一区| 久久婷婷人人爽人人干人人爱| 日韩亚洲欧美综合| 久久这里有精品视频免费| 欧美日本亚洲视频在线播放| 成人午夜高清在线视频| 国产精品不卡视频一区二区| 日产精品乱码卡一卡2卡三| 欧美性猛交╳xxx乱大交人| 老司机福利观看| 亚洲第一区二区三区不卡| 日本猛色少妇xxxxx猛交久久| 国产成人aa在线观看| 18+在线观看网站| 久久国产乱子免费精品| 好男人视频免费观看在线| 精品久久久噜噜| 久久久午夜欧美精品| av视频在线观看入口| 真实男女啪啪啪动态图| 亚洲真实伦在线观看| 狠狠狠狠99中文字幕| 国产 一区 欧美 日韩| 亚洲精品aⅴ在线观看| 日韩成人av中文字幕在线观看| 蜜桃久久精品国产亚洲av| 永久免费av网站大全| 偷拍熟女少妇极品色| 亚洲图色成人| 人妻夜夜爽99麻豆av| 免费看美女性在线毛片视频| 尤物成人国产欧美一区二区三区| 91久久精品电影网| 一区二区三区免费毛片| 国产精品一区二区在线观看99 | 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 久久久久久久午夜电影| 美女被艹到高潮喷水动态| 亚州av有码| 久久久久久久久久久丰满| 97人妻精品一区二区三区麻豆| 免费黄色在线免费观看| 男女啪啪激烈高潮av片| 99久久精品国产国产毛片| 久久鲁丝午夜福利片| 免费观看人在逋| 亚洲欧美成人精品一区二区| 狂野欧美激情性xxxx在线观看| 免费人成在线观看视频色| 国产成年人精品一区二区| 午夜久久久久精精品| 亚洲美女视频黄频| 女人被狂操c到高潮| 日本五十路高清| 免费黄色在线免费观看| 婷婷色麻豆天堂久久 | 国产一区二区在线观看日韩| 亚洲av不卡在线观看| 中文欧美无线码| 国产视频内射| 村上凉子中文字幕在线| 成人二区视频| 2021天堂中文幕一二区在线观| 精品无人区乱码1区二区| 欧美一区二区精品小视频在线| 毛片女人毛片| 成人亚洲欧美一区二区av| 嫩草影院入口| 在线观看66精品国产| 成人无遮挡网站| 九九在线视频观看精品| 亚洲精品国产av成人精品| 99视频精品全部免费 在线| 97热精品久久久久久| 亚洲成色77777| 欧美zozozo另类| 亚洲国产精品合色在线| 丝袜美腿在线中文| 人人妻人人澡人人爽人人夜夜 | 男插女下体视频免费在线播放| 人体艺术视频欧美日本| 免费观看的影片在线观看| 91午夜精品亚洲一区二区三区| 国产亚洲精品久久久com| 亚洲国产高清在线一区二区三| 国产免费视频播放在线视频 | 免费播放大片免费观看视频在线观看 | 午夜久久久久精精品| 婷婷色综合大香蕉| 国产黄色视频一区二区在线观看 | 久久精品久久精品一区二区三区| 赤兔流量卡办理| 天堂√8在线中文| 欧美潮喷喷水| 国产精品爽爽va在线观看网站| 久久久欧美国产精品| 国产av不卡久久| 欧美一级a爱片免费观看看| 国产精品人妻久久久影院| 精品久久久噜噜| 日本猛色少妇xxxxx猛交久久| 麻豆一二三区av精品| 我要看日韩黄色一级片| 久久鲁丝午夜福利片| 男女视频在线观看网站免费| 国产老妇女一区| 亚洲欧美清纯卡通| 国产免费福利视频在线观看| 国产真实乱freesex| 精品一区二区三区人妻视频| 国产精品久久久久久精品电影小说 | 日本免费a在线| 国产成人freesex在线| 女人十人毛片免费观看3o分钟| 亚洲,欧美,日韩| 91在线精品国自产拍蜜月| 久久这里只有精品中国| 亚洲美女搞黄在线观看| 99久久精品热视频| 夜夜爽夜夜爽视频| 成年女人永久免费观看视频| 亚洲精华国产精华液的使用体验| 综合色丁香网| 网址你懂的国产日韩在线| 亚洲五月天丁香| 老司机影院毛片| 午夜精品在线福利| 日韩一本色道免费dvd| 午夜激情欧美在线| 蜜臀久久99精品久久宅男| 日本黄色片子视频| eeuss影院久久| 晚上一个人看的免费电影| 久久精品夜夜夜夜夜久久蜜豆| 色尼玛亚洲综合影院| 五月玫瑰六月丁香| 在线播放无遮挡| 国产高清有码在线观看视频| 熟女电影av网| 狂野欧美白嫩少妇大欣赏| 黄片无遮挡物在线观看| 天堂影院成人在线观看| 国产亚洲午夜精品一区二区久久 | 一区二区三区免费毛片| 搞女人的毛片| 听说在线观看完整版免费高清| 男女那种视频在线观看| av免费观看日本| 建设人人有责人人尽责人人享有的 | 欧美3d第一页| 亚洲av中文av极速乱| 能在线免费观看的黄片| 人妻夜夜爽99麻豆av| 国产亚洲av嫩草精品影院| 国产精品1区2区在线观看.| 色播亚洲综合网| 女人久久www免费人成看片 | 欧美最新免费一区二区三区| av在线天堂中文字幕| 国产亚洲精品av在线| 国产高清国产精品国产三级 | 两个人的视频大全免费| 毛片女人毛片| 日韩精品青青久久久久久| 一卡2卡三卡四卡精品乱码亚洲| 黑人高潮一二区| 亚洲国产日韩欧美精品在线观看| 卡戴珊不雅视频在线播放| 国产免费一级a男人的天堂| 欧美3d第一页| 久久久久久久久久黄片| 亚洲性久久影院| 久久午夜福利片| 久久精品人妻少妇| 高清毛片免费看| 亚洲欧洲日产国产| 女人久久www免费人成看片 | 欧美bdsm另类| 国产精品电影一区二区三区| 亚洲欧美日韩卡通动漫| 一区二区三区高清视频在线| 亚洲国产成人一精品久久久| 久久人人爽人人爽人人片va| 国内精品一区二区在线观看| 亚洲欧美日韩高清专用| 婷婷色综合大香蕉| 成人一区二区视频在线观看| 丝袜美腿在线中文| 精品人妻偷拍中文字幕| 国产91av在线免费观看| 亚洲内射少妇av| 成人高潮视频无遮挡免费网站| 男人和女人高潮做爰伦理| 成人三级黄色视频| 啦啦啦观看免费观看视频高清| 免费大片18禁| 国产午夜福利久久久久久| 久久人人爽人人片av| 菩萨蛮人人尽说江南好唐韦庄 | av又黄又爽大尺度在线免费看 | 狂野欧美激情性xxxx在线观看| 国产精品1区2区在线观看.| 最近的中文字幕免费完整| 亚洲怡红院男人天堂| 99热精品在线国产| 亚洲国产精品成人综合色| 精品人妻偷拍中文字幕| 国产精品一二三区在线看| 神马国产精品三级电影在线观看| 99久久精品热视频| 亚洲av不卡在线观看| 看片在线看免费视频| 日本爱情动作片www.在线观看| 国产精品一区二区性色av| 最近视频中文字幕2019在线8| 国产精品.久久久| 国产黄色小视频在线观看| 熟妇人妻久久中文字幕3abv| 老司机影院成人| av在线蜜桃| 精品99又大又爽又粗少妇毛片| 亚洲av福利一区| 日韩在线高清观看一区二区三区| 亚洲综合精品二区| 九九在线视频观看精品| 又爽又黄无遮挡网站| 夜夜看夜夜爽夜夜摸| 一级av片app| 久久99蜜桃精品久久| 中文在线观看免费www的网站| 久久久久久伊人网av| 欧美日韩国产亚洲二区| 综合色丁香网| 日韩欧美 国产精品| 老司机影院成人| 纵有疾风起免费观看全集完整版 | av卡一久久| 男的添女的下面高潮视频| 建设人人有责人人尽责人人享有的 | 日韩av不卡免费在线播放| 欧美xxxx性猛交bbbb| 日日撸夜夜添| 99热这里只有是精品50| 中文字幕熟女人妻在线| 亚洲三级黄色毛片| 精品一区二区免费观看| 国产日韩欧美在线精品| av在线播放精品| 欧美激情久久久久久爽电影| 午夜激情欧美在线| 婷婷色综合大香蕉| 少妇的逼水好多| 美女cb高潮喷水在线观看| 中文天堂在线官网| 精品久久久久久久久久久久久| 国产国拍精品亚洲av在线观看| 亚洲成色77777| 性插视频无遮挡在线免费观看| 亚洲精品自拍成人| 久久人妻av系列| 九九热线精品视视频播放| 日本免费在线观看一区| 三级毛片av免费| 免费不卡的大黄色大毛片视频在线观看 | 六月丁香七月| 日本五十路高清| 欧美高清性xxxxhd video| 国产精品爽爽va在线观看网站| 熟妇人妻久久中文字幕3abv| 熟女人妻精品中文字幕| av免费在线看不卡| 赤兔流量卡办理| 国产亚洲精品久久久com| 久久久久久九九精品二区国产| 亚洲天堂国产精品一区在线| 亚洲av成人av| 白带黄色成豆腐渣| 少妇熟女欧美另类| 国产真实伦视频高清在线观看| 国产日韩欧美在线精品| 国产黄a三级三级三级人| 国产免费男女视频|