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

    膠結(jié)材料加載破壞的能量積聚與釋放響應(yīng)研究

    2020-04-15 01:59:30王明旭
    關(guān)鍵詞:巖爆分貝張量

    王明旭

    膠結(jié)材料加載破壞的能量積聚與釋放響應(yīng)研究

    王明旭1,2

    (1. 武漢科技大學(xué)資源與環(huán)境工程學(xué)院,湖北,武漢 430081;2. 惠州市國(guó)土資源局,廣東,惠州 516003)

    為研究膠結(jié)材料在0.001 kN/s的速度緩慢加載達(dá)到峰值強(qiáng)度后出現(xiàn)的第二跌落特征,并伴隨發(fā)生在第二跌落區(qū)出現(xiàn)的噪音分貝值的突增現(xiàn)象產(chǎn)生的原因,應(yīng)用能量積聚與釋放的總能量平衡理論,對(duì)膠結(jié)充填體軸向加載試驗(yàn)獲取全應(yīng)力應(yīng)變曲線,并通過電子噪音測(cè)試設(shè)備測(cè)試加載過程中的噪音分貝值的大小,引入可表征試件全應(yīng)力應(yīng)變曲線的即時(shí)彈性模量與割線彈性模量比值,作為柔化正態(tài)張量和可表征試件不同等應(yīng)變變化的應(yīng)變能與總應(yīng)變能比值作為應(yīng)變能張量,通過柔化正態(tài)張量和應(yīng)變能張量的變化大小,表征試件在能量積聚階段和釋放響應(yīng)階段的變化特點(diǎn),特別是在第二跌落階段的能量變化特點(diǎn)。同時(shí)通過噪音分貝測(cè)量方法試驗(yàn),建立噪音分貝與試件能量釋放大小之間的對(duì)應(yīng)關(guān)系。研究結(jié)果表明,只要加載速度較小,達(dá)到0.001 kN/s左右時(shí),膠結(jié)充填體才會(huì)出現(xiàn)明顯的第二跌落區(qū)。通過電子噪音分貝測(cè)試第二跌落區(qū)破壞出現(xiàn)的噪音分貝突增15 dB,相較常規(guī)加載多增加了10 dB以上。同時(shí)將這種緩慢加載方式應(yīng)用在礦巖加載試驗(yàn)中,其加載破壞時(shí)的噪音分貝值突增更加明顯,增加值達(dá)到20~35 dB。表明試件在近似等壓受載環(huán)境下,更易積聚能量,并在破壞時(shí)以突然釋放的形式釋放,說明能量釋放的強(qiáng)度是以前期積聚能量為前提的,這將有助于指導(dǎo)金屬礦山采場(chǎng)高應(yīng)力區(qū)評(píng)價(jià),維護(hù)采場(chǎng)安全。

    膠結(jié)材料;能量釋放;彈性模量;應(yīng)變能;噪音能量

    隨著采深加大,礦山井下巖爆風(fēng)險(xiǎn)也逐漸增大,對(duì)充填空區(qū)未充填前的礦巖受載特征需要有較清晰的認(rèn)識(shí),為此礦巖受載出現(xiàn)的能量積聚與釋放情況值得關(guān)注。膠結(jié)材料充填井下,具備較大的儲(chǔ)能能力,其在與圍巖共同作用時(shí),在各種應(yīng)力作用下,膠結(jié)材料積聚能量。而對(duì)于圍巖,在能量積聚后將伴隨有不同形式的能量釋放。對(duì)巖爆的研究,一些學(xué)者已經(jīng)取得了相關(guān)成果,如何滿潮等[1]設(shè)計(jì)并研發(fā)沖擊巖爆試驗(yàn)系統(tǒng),并進(jìn)行沖擊巖爆試驗(yàn);夏元友等[2]通過自主研發(fā)的巖爆模擬試驗(yàn)裝置對(duì)大尺寸試件進(jìn)行巖爆試驗(yàn),探討試件發(fā)生巖爆的烈度與碎屑分形維數(shù)的關(guān)系;龔劍等[3]通過AHP-TOPSI評(píng)判模型對(duì)巖爆傾向性進(jìn)行預(yù)測(cè);蔡美峰等[4]優(yōu)化應(yīng)力解除測(cè)量技術(shù)并在礦區(qū)深部進(jìn)行現(xiàn)場(chǎng)實(shí)測(cè);賈義鵬等[5]為評(píng)價(jià)巖爆發(fā)生的可能性,提出一種基于粒子群算法和廣義回歸神經(jīng)網(wǎng)絡(luò)模型(PSO-GRNN模型)的巖爆預(yù)測(cè)方法;張文東等[6]、YU Yang等[7]將微震監(jiān)測(cè)技術(shù)用于隧洞工程施工中的巖爆監(jiān)測(cè)預(yù)警;ZHOU Keping等[8]采用熵權(quán)-云模型對(duì)巖爆等級(jí)進(jìn)行預(yù)測(cè);ZHOU Jian等[9]使用隨機(jī)梯度增強(qiáng)(SGB)方法預(yù)測(cè)巖爆破壞;DONG Longjun等[10]認(rèn)為巖爆的主要因素是原位應(yīng)力值、巖石的單軸抗壓強(qiáng)度和抗拉強(qiáng)度值;趙康等[11]利用掃描電鏡(SEM)對(duì)巖爆巖石斷口微觀形貌特征進(jìn)行研究分析,從微觀角度探索巖爆產(chǎn)生的機(jī)理;尚彥軍等[12]在已建立的應(yīng)變型巖爆五因素綜合判據(jù)基礎(chǔ)上表征巖爆趨勢(shì)和強(qiáng)度;馬春馳等[13]采用顆粒流程序PFC3D對(duì)不同圍壓三軸卸荷下的巖爆效應(yīng)進(jìn)行模擬分析;蘇國(guó)韶等[14]利用自主研發(fā)的真三軸巖爆試驗(yàn)系統(tǒng),實(shí)現(xiàn)巖爆碎塊耗能組成的定量化分析。這些都是對(duì)巖體內(nèi)部能量積聚與釋放展開的相關(guān)研究,能夠通過這些方法對(duì)具有脆性破壞特性的硬巖開展巖爆研究。而如何從硬巖過渡到人工膠結(jié)材料,如何借助這些對(duì)硬巖的分析手段和方法研究膠結(jié)材料的能量積聚與釋放,需要?jiǎng)?chuàng)新方法,明晰膠結(jié)材料能量積聚與釋放存在的特點(diǎn)。

    膠結(jié)材料充填后,將形成膠結(jié)材料與礦巖的相互作用,在荷載作用下,不僅礦巖自身在不斷蓄積能量,而且將能量傳遞給膠結(jié)材料。在能量的傳遞過程中,膠結(jié)充填體對(duì)礦巖產(chǎn)生什么樣的影響,是一個(gè)值得深入探索的問題。目前對(duì)于兩者之間的相互作用,比較成熟的研究成果是能量匹配分析法和應(yīng)變非協(xié)調(diào)變形分析法。這兩種方式主要研究膠結(jié)材料和礦巖直接加載破壞前的情況,而對(duì)于礦巖在循環(huán)加卸載以及恒載維持后再加載的情況注意較少,而這正是礦巖或人工澆注體能量積聚與釋放的關(guān)鍵,也是容易產(chǎn)生巖爆或能量突然釋放并伴隨破壞的發(fā)生。在對(duì)膠結(jié)材料進(jìn)行加載試驗(yàn)的過程中,發(fā)現(xiàn)膠結(jié)材料在加載破壞后的應(yīng)力應(yīng)變曲線存在局部凸型特征。通過對(duì)試件的應(yīng)力應(yīng)變曲線進(jìn)行相關(guān)研究,相關(guān)學(xué)者主要是對(duì)巖石試件的應(yīng)力應(yīng)變曲線研究,如姜春露等[15]利用伺服滲透試驗(yàn)結(jié)果,分析脆性巖石全應(yīng)力應(yīng)變過程中滲透性演化和加載變形破壞特征;鄧林等[16]提出基于巖石單軸壓縮應(yīng)力應(yīng)變?nèi)^程曲線的巖爆能量指數(shù)法進(jìn)行巖爆預(yù)測(cè);于德海等[17]對(duì)兩種片巖進(jìn)行力學(xué)試驗(yàn),得到巖石全應(yīng)力應(yīng)變壓縮過程,提出峰前壓密、線彈性、微裂隙穩(wěn)態(tài)擴(kuò)展和微裂隙不穩(wěn)態(tài)擴(kuò)展4個(gè)變形階段以及峰后應(yīng)變軟化、殘余應(yīng)力穩(wěn)定兩個(gè)變形階段;姜永東等[18]認(rèn)為巖石的損傷演化具有分叉和混沌特征。而對(duì)于充填體相關(guān)應(yīng)力應(yīng)變曲線特征的研究相對(duì)較少,如楊天雨等[19]為了探求尾砂膠結(jié)充填體在循環(huán)加卸載下能量的變化特征,利用電液伺服試驗(yàn)機(jī)對(duì)分級(jí)尾砂膠結(jié)充填體進(jìn)行單軸循環(huán)加卸載試驗(yàn);楊逾等[20]認(rèn)為充填材料在壓縮變形的過程中會(huì)吸收能量,對(duì)于煤礦的充填開采,充填材料的能量變化規(guī)律會(huì)對(duì)工作面的能量產(chǎn)生很大的影響;徐文彬等[21]認(rèn)為充填體的變形破壞、能量耗散與圍壓的變化密切相關(guān);劉玉龍等[22]引入充填體變形比能與巖體釋放比能之比的匹配系數(shù)。本文考慮到配比材料加載破壞后的能量積聚釋放,提出了礦巖或人工澆注體能量的釋放大小是與前期能量的積聚多少相關(guān)的設(shè)想,并通過提出的一種蠕變加載能量積聚的試驗(yàn)方法對(duì)這一設(shè)想驗(yàn)證,進(jìn)行了系列試驗(yàn)研究。

    1 微元體破壞的能量積聚與釋放響應(yīng)動(dòng)態(tài)調(diào)整假設(shè)

    1.1 傳統(tǒng)的試件單軸壓縮條件下的巖石變形特征

    巖石試件在單軸壓縮荷載作用下產(chǎn)生變形的過程,存在四個(gè)階段。1)孔隙裂隙壓密階段,即試件中原有張開性結(jié)構(gòu)面或微裂隙逐漸閉合,巖石被壓密,形成早期的非線性變形。在此階段試件橫向膨脹較小,試件體積隨荷載增大而減小。2)彈性變形至微彈性裂隙穩(wěn)定發(fā)展階段,該階段的應(yīng)力-應(yīng)變曲線近似直線型,存在微破裂穩(wěn)定發(fā)展階段。3)非穩(wěn)定破裂發(fā)展階段,此階段微破裂的發(fā)展出現(xiàn)了質(zhì)的變化,破裂不斷發(fā)展,直至試件完全破壞。4)破裂后階段,巖塊承載力達(dá)到峰值強(qiáng)度后,其內(nèi)部結(jié)構(gòu)造成破壞,但試件基本保持整體狀。本階段裂隙快速發(fā)展,交叉且相互聯(lián)合形成宏觀斷裂面。這四個(gè)階段是從試件單軸壓縮加載過程中的變形角度出發(fā),而試件加載快慢對(duì)其產(chǎn)生的影響,特別是每個(gè)較小階段的近恒載過程,其內(nèi)部損傷破壞的各微元體相互作用機(jī)理,是否存在單純從試件變形情況分析無法解釋的情況,從本文提出的蠕變加載方式的能量積聚與釋放響應(yīng)的試驗(yàn)來看,顯然需要進(jìn)一步研究試件在單軸壓縮過程中的內(nèi)部微元體損傷破壞情況。從能量積聚與釋放響應(yīng)的角度考慮,以微元體為分析對(duì)象,為此提出了微元體破壞的統(tǒng)計(jì)學(xué)組合效應(yīng),存在微元體破壞的能量積聚與釋放響應(yīng)動(dòng)態(tài)調(diào)整假設(shè),并認(rèn)為內(nèi)部裂紋擴(kuò)展過程中存在非規(guī)則鋸齒狀破壞。

    1.2 微元體破壞的能量積聚與釋放響應(yīng)動(dòng)態(tài)調(diào)整假設(shè)

    試樣在加載過程中,內(nèi)部損傷不斷積累,當(dāng)內(nèi)部裂紋破壞所需要的載荷達(dá)到之后有一個(gè)微小的能量釋放及泄壓回彈。如果內(nèi)部損傷破壞速度大于外部荷載激發(fā)內(nèi)部破壞速度。加載方式,無論是試驗(yàn)力控制還是位移控制,即加載速率控制在足夠小的范圍,形成試樣在加載過程中的能量積聚環(huán)境。單軸加載速度對(duì)巖石的變形性質(zhì)和強(qiáng)度指標(biāo)有明顯影響,加載速率越大,測(cè)得的彈性模量愈大,獲得的強(qiáng)度指標(biāo)值越高。

    整體表現(xiàn)出的應(yīng)力應(yīng)變曲線,是試件內(nèi)部各個(gè)微元體裂紋擴(kuò)展破壞的統(tǒng)計(jì)學(xué)的組合效應(yīng)。即試件在加載過程中,外部施加荷載導(dǎo)致試件內(nèi)部不斷發(fā)生應(yīng)力重分布。試件的整體穩(wěn)定性是由各微元體的統(tǒng)計(jì)學(xué)穩(wěn)定決定的。在外部荷載作用下,各微元體區(qū)域主要經(jīng)歷如下過程。

    第一過程:各微元體承受荷載逐漸加大;

    第二過程:部分弱微元體開始發(fā)生破壞,產(chǎn)生內(nèi)部損傷;

    第三過程:試件重新形成新的自穩(wěn)體;

    第四過程:以此循環(huán),直到試件最后破壞。

    微元體破壞,形成應(yīng)力的釋放,有助于微元體局部形成增大的Δ抵抗能力,當(dāng)持續(xù)加載的荷載增量沒有超過局部微元體增大的抵抗能力時(shí),試件內(nèi)部進(jìn)行內(nèi)部的自穩(wěn)能力的構(gòu)建,形成更大的自穩(wěn)體。如果加載速度太快,加載速率過大,試件在未形成新的自穩(wěn)體前發(fā)生破壞,不僅試件內(nèi)部能量無法得到積聚的機(jī)會(huì),也會(huì)不斷產(chǎn)生新的破壞,直至試件最后發(fā)生破壞。在加載過程中,試樣內(nèi)部發(fā)生損傷破壞時(shí),當(dāng)試樣內(nèi)部在即將達(dá)到破壞時(shí),存在能量過度積累,當(dāng)破壞發(fā)生時(shí),存在破壞的“滑鐵盧”效應(yīng),致使內(nèi)部積聚的能量過量釋放,致使試樣內(nèi)部有一個(gè)泄壓調(diào)整階段,能夠促發(fā)試樣內(nèi)部向自穩(wěn)能力構(gòu)建的方向發(fā)展,也伴隨能量的積聚。如果在試樣內(nèi)部自穩(wěn)協(xié)調(diào)過程中,加載速度偏快,內(nèi)部自穩(wěn)協(xié)調(diào)未能完成就在新的加載作用下發(fā)生新的破壞,無法形成能量積聚條件,試樣將形成不斷積聚又釋放的動(dòng)態(tài)過程,而如果內(nèi)部自穩(wěn)協(xié)調(diào)能夠較好完成,試樣將積聚大量能量,試樣內(nèi)部的破壞裂紋較少,其在能量積聚的過程中,將會(huì)在最后突然釋放,形成大量的能量釋放。

    1.3 假設(shè)能量積聚與能量釋放響應(yīng)能量守恒

    膠結(jié)材料在加載過程中的應(yīng)力應(yīng)變曲線,在加載速度非常?。?.001 kN/s)的時(shí)候,在膠結(jié)材料達(dá)到最大抗壓強(qiáng)度之后,會(huì)出現(xiàn)第二次跌落階段(見圖1),此次跌落在加載破壞后的應(yīng)力應(yīng)變曲線存在二個(gè)階段:

    1)能量積聚階段:試件在加載達(dá)到最大抗壓強(qiáng)度之前,即應(yīng)力應(yīng)變曲線的即時(shí)彈性模量從緩慢增加到快速增加再到逐漸減小為零的過程中。

    2)能量釋放階段:在能量釋放階段,因試件材料的不同,其全應(yīng)力應(yīng)變曲線呈現(xiàn)兩種情況。①單一的跌落階段。即試件在加載過程中,其全應(yīng)力應(yīng)變曲線的最大抗壓強(qiáng)度右邊部分,隨著應(yīng)變值的增大,其應(yīng)力出現(xiàn)了平滑的減小過程。②雙跌落過程。在部分試件材料中,其全應(yīng)力應(yīng)變曲線的最大抗壓強(qiáng)度的右邊部分出現(xiàn)平滑的跌落與突變跌落過程。

    在膠結(jié)充填體的加載過程中,當(dāng)加載速度非常慢,等同于緩慢恒載時(shí),膠結(jié)充填體出現(xiàn)了雙跌落過程。這種雙跌落現(xiàn)象,也伴隨著試件加載過程中聲音分貝值的增加。膠結(jié)充填體出現(xiàn)了聲音分貝值的突然增加和對(duì)應(yīng)的應(yīng)力應(yīng)變曲線的第二次跌落過程存在的對(duì)應(yīng)關(guān)系,說明了這種跌落過程中伴隨著能量的釋放。

    對(duì)于這種雙跌落現(xiàn)象如何科學(xué)界定?這是由加載速度過慢所導(dǎo)致。因此,可以認(rèn)為膠結(jié)充填體存在的雙跌落過程與其在加載過程中能量積聚與突然釋放有關(guān)。

    圖1 全應(yīng)力應(yīng)變曲線出現(xiàn)的第二跌落區(qū)域

    第二跌落區(qū)是試件達(dá)到最大抗壓強(qiáng)度之后,出現(xiàn)了表壁加載破壞而剝落,而內(nèi)部因?yàn)榫徛虞d導(dǎo)致壓實(shí)積聚能量。隨著試件剝落程度加大,試件內(nèi)部開始發(fā)生破壞,在能量釋放過程中出現(xiàn)。

    應(yīng)力應(yīng)變曲線主要分為兩個(gè)階段:

    第一階段:能量積聚階段,主要是試件達(dá)到最大抗壓強(qiáng)度之前的應(yīng)力應(yīng)變曲線階段。

    第二階段:能量釋放階段,主要是試件達(dá)到最大抗壓強(qiáng)度之后,開始出現(xiàn)較大的損傷破壞,這種損傷破壞從內(nèi)部延伸至表壁,最終導(dǎo)致試件整體破壞。對(duì)于能量釋放階段,脆性試件特別是硬巖試件發(fā)生脆性破壞;而對(duì)于軟巖以及人工材料,在荷載速度加載非常慢時(shí),即等同于階段時(shí)間恒載時(shí),能力釋放階段會(huì)出現(xiàn)雙跌落現(xiàn)象。如劉智高[23]在使用電液伺服巖石三軸試驗(yàn)進(jìn)行粗砂巖的加載時(shí),其應(yīng)力應(yīng)變曲線也出現(xiàn)了明顯的第二跌落現(xiàn)象。

    物體在外界因素影響下產(chǎn)生變形,在微小的時(shí)間間隔內(nèi),物體從一種狀態(tài)過渡到另一種狀態(tài),物體的動(dòng)能和物體自身的應(yīng)變能變化等于物體與外界進(jìn)行的熱交換和外界對(duì)物體做的功,為此根據(jù)熱力學(xué)第一定律,總能量的變化為:

    式中:E—試件的動(dòng)能;—試件的應(yīng)變能;—體力和面力做功;—試件與周圍介質(zhì)的熱量交換。

    兩種加載方式:一種是以0.001 mm/s的速度對(duì)試件(D×H=50 mm×100 mm)進(jìn)行加載,一種是以0.001 kN/s的速度對(duì)試件(D×H=50 mm×100 mm)進(jìn)行加載。后者的全應(yīng)力應(yīng)變曲線在能量釋放階段出現(xiàn)第二跌落區(qū),并且伴隨明顯高于第一種加載方式破壞時(shí)的聲響。

    在緩慢加載過程中,近似等壓過程,試件的動(dòng)能可以忽略,即為0。

    巖石變形前期以彈性應(yīng)變能的方式存儲(chǔ)外界提供的能量,同時(shí)又通過損傷演化等向外界耗散能量;變形的后期以劇烈的能量釋放為主。應(yīng)變能的大小為[24]:

    巖石變形的微觀機(jī)制與能量關(guān)系,巖石的應(yīng)力應(yīng)變曲線,從宏觀上可以分為峰值點(diǎn)前的應(yīng)變硬化階段和峰值點(diǎn)后應(yīng)變軟化階段,對(duì)于峰值點(diǎn)前的應(yīng)變硬化階段對(duì)應(yīng)的能量變化過程:

    巖石將外界傳遞過來的能量轉(zhuǎn)變?yōu)樽陨淼膬?nèi)部應(yīng)變能,所反映的是外界能量和巖石內(nèi)能之間的關(guān)系。

    巖石將內(nèi)部的能量以不同的形式又釋放到外界,即利用率高的能量轉(zhuǎn)變?yōu)槔寐实偷哪芰?,此過程反映巖石內(nèi)能與外界耗散能的變化關(guān)系。

    2 能量積聚與釋放響應(yīng)試驗(yàn)

    按照質(zhì)量比1:8(32.5#普通硅酸鹽水泥與礦山尾礦庫(kù)尾砂)的配比制作膠結(jié)材料(尺寸為D = 50 mm,H=100 mm),濃度控制在68%。養(yǎng)護(hù)14d后進(jìn)行加載試驗(yàn),得到一組5個(gè)試樣的應(yīng)力應(yīng)變曲線(見圖2)。從圖2可以看出,試驗(yàn)加載破壞前的曲線都類同,同時(shí)在加載破壞后的應(yīng)力應(yīng)變曲線中也都存在凸型特征。從試驗(yàn)過程中的相應(yīng)時(shí)刻的加載破壞情況來看,主要是試樣在加載破壞后,由裂紋擴(kuò)展破壞積累成為突然崩塌的整體性破壞,同時(shí)也伴隨著一聲脆響,脆響之時(shí)正是膠結(jié)材料應(yīng)變曲線中出現(xiàn)凸型特征變化之時(shí)。

    圖2 膠結(jié)材料的應(yīng)力應(yīng)變曲線

    2.1 能量積聚階段

    將試樣的應(yīng)力應(yīng)變曲線分為兩部分進(jìn)行分析研究。一是加載破壞前能量積聚階段,二是加載破壞后的殘余應(yīng)力變化和能量釋放階段。從多個(gè)試樣加載破壞前的應(yīng)力應(yīng)變曲線來看,都具有相同的變化趨勢(shì),只是不同的試樣由于內(nèi)部孔隙的差異,存在協(xié)調(diào)變形的次數(shù)不一樣。從圖3可以看出,試樣1在應(yīng)變達(dá)到0.005時(shí)開始接連出現(xiàn)2次協(xié)調(diào)變形階段,而試樣2在應(yīng)變?yōu)?.003就開始出現(xiàn)了應(yīng)變協(xié)調(diào),然后越過應(yīng)變0.005,在應(yīng)變達(dá)到0.0065開始出現(xiàn)接連兩次的應(yīng)變協(xié)調(diào)過程,之后試樣出現(xiàn)穩(wěn)步的應(yīng)力隨應(yīng)變不斷增長(zhǎng)的過程,直到試樣將要破壞時(shí),曲線趨于平緩直至破壞。

    圖3 膠結(jié)材料應(yīng)力應(yīng)變曲線

    為了研究分析的方便,選取編號(hào)為NO.1的試樣進(jìn)行應(yīng)力應(yīng)變曲線特征和相關(guān)的計(jì)算分析。對(duì)NO.1試樣破壞前的應(yīng)力應(yīng)變關(guān)系進(jìn)行曲線擬合,其擬合曲線方程為:

    = 0.06573+170.40464-49890.3036^2+1.43991E 7^3-1.03384E9^4+2.07235E10^5,2=0.9931 (5)

    按照應(yīng)力應(yīng)變曲線,當(dāng)應(yīng)變?yōu)?時(shí),其相應(yīng)的應(yīng)力為0,(5)式出現(xiàn)常數(shù)項(xiàng),允許保留,作為實(shí)際試樣(隨機(jī)選取試樣)加載過程中各種可能的復(fù)雜情況反映。不過,為了與理論解接軌,相應(yīng)的限定條件是應(yīng)變>0。

    對(duì)No.1試樣求積分得到求解累計(jì)應(yīng)變能的解析解:

    Y= c+0.06573+85.20232^2-16630.1012^3+ 0.3599775E7^4-0.206768E9^5+0.34539167E10^6(6)

    依據(jù)邊界條件:|=0=0 代入得到:=0,故

    Y=0.06573+85.20232^2-16630.1012^3+0.3599775E7^4-0.206768E9^5+0.34539167E10^6 (7)

    將No.1試樣的應(yīng)力應(yīng)變曲線中達(dá)到最大破壞強(qiáng)度值時(shí)對(duì)應(yīng)的應(yīng)變值25等分(等分?jǐn)?shù)可進(jìn)行調(diào)整,此處選用了25等分),通過積分得到的應(yīng)變能方程分析等應(yīng)變變化情況下的應(yīng)變能累加值(見圖4)和應(yīng)變能變化值規(guī)律。從圖4可以看出,在相同應(yīng)變變化范圍內(nèi),應(yīng)變能的累加值與應(yīng)力應(yīng)變曲線變化規(guī)律類同,而等應(yīng)變內(nèi)的應(yīng)變能變化呈現(xiàn)先緩慢增加然后逐漸加快的變化趨勢(shì)。

    圖4 應(yīng)變能變化曲線

    同樣對(duì)No.1試樣的曲線擬合方程求一階導(dǎo):

    Y′=170.40464-99780.6072x^+4.31973E7x^2- 4.13536E9x^3+10.36175E10x^4, (8)

    代入應(yīng)變值,即為不同應(yīng)變值對(duì)應(yīng)的即時(shí)彈性模量,計(jì)算出相應(yīng)等應(yīng)變變化的彈性模量變化量。從圖5可以看出,隨著應(yīng)變值的不斷增大,彈性模量的變化量也呈現(xiàn)不斷增大的變化趨勢(shì)。

    圖5 彈性模量的等應(yīng)變變化量

    為了研究比較的方便,提出柔化正態(tài)張量概念,即不同等應(yīng)變下的即時(shí)彈性模量變化量與破壞時(shí)的割線模量的比值。假設(shè)等應(yīng)變的即時(shí)彈性模量變化為△Eεi,割線彈性模量為E,則柔化正態(tài)張量:

    式中:-為試樣破壞時(shí)的最大應(yīng)變值等分。

    通過對(duì)No.1試樣的割線模量計(jì)算,其值E為269.635 MPa,從而得到等應(yīng)變值的柔化正態(tài)張量,建立柔化正態(tài)張量與應(yīng)變能和應(yīng)力之間的相互關(guān)系。將No.1試樣破壞時(shí)的應(yīng)變值均分為25等分,根據(jù)(4)式得出不同等應(yīng)變值對(duì)應(yīng)的應(yīng)變能大小,根據(jù)(1)式得出不同等應(yīng)變值對(duì)應(yīng)的應(yīng)力大小,繪制曲線(見圖6)。

    圖6 應(yīng)變-應(yīng)力-累積應(yīng)變能-柔化正態(tài)張量關(guān)系曲線

    從圖6可以看出,應(yīng)力隨應(yīng)變變化的趨勢(shì),主要存在三個(gè)階段。第一個(gè)階段是增長(zhǎng)緩慢階段,這個(gè)階段主要是應(yīng)變值在0~0.0028之間;第二個(gè)階段是快速增長(zhǎng)階段,這個(gè)階段的應(yīng)變值在0.0028~0.011之間;第三個(gè)階段是增速減緩階段,在最后破壞時(shí)曲線趨于水平。而應(yīng)變能與應(yīng)變的關(guān)系,先是緩慢增長(zhǎng),接著不斷快速增長(zhǎng),應(yīng)變值達(dá)到0.008時(shí)呈現(xiàn)近線性增長(zhǎng)。對(duì)于柔化正態(tài)張量,在加載前期,其隨應(yīng)變?cè)鲩L(zhǎng)緩慢,當(dāng)應(yīng)變達(dá)到0.004時(shí)呈現(xiàn)區(qū)域線性增長(zhǎng)??紤]到應(yīng)變能和柔化正態(tài)張量隨應(yīng)變?cè)诤笃诔霈F(xiàn)的線性變化特性,進(jìn)行相同應(yīng)變下的柔化正態(tài)張量與應(yīng)變能之間的關(guān)系曲線繪制(見圖7)。從圖7可以看出,柔化正態(tài)張量隨應(yīng)變能,先是快速增長(zhǎng),當(dāng)應(yīng)變能達(dá)到0.005 MPa時(shí),其變化呈現(xiàn)線性趨勢(shì)。

    圖7 試樣柔化正態(tài)張量與累積應(yīng)變能關(guān)系曲線

    對(duì)膠結(jié)材料試樣進(jìn)行加載,其應(yīng)力應(yīng)變曲線主要反映了試樣的力學(xué)變形特征及能量大小。關(guān)于反映試樣彈性模量變化特點(diǎn)的柔化正態(tài)張量與反映試樣應(yīng)變能變化的應(yīng)變能張量,兩者之間是否存在關(guān)系?為此提出應(yīng)變能張量,即為等間距應(yīng)變能W與總應(yīng)變能W的比值,則:

    通過NO.1試樣的柔化正態(tài)張量與應(yīng)變能張量進(jìn)行計(jì)算,并擬合曲線(見圖8)。從圖7可以看出,在應(yīng)變能較小時(shí),柔化正態(tài)張量變化曲線幾乎是一條水平線,當(dāng)應(yīng)變能將要達(dá)到最大時(shí),柔化正態(tài)張量變化曲線幾乎是一條豎直線,而中間過程呈現(xiàn)近線性變化特征。

    2.2 能量釋放響應(yīng)階段

    從試樣加載破壞后的應(yīng)力應(yīng)變曲線(見圖9)可以看出,膠結(jié)材料試樣在破壞之后,隨著應(yīng)變的增長(zhǎng),應(yīng)力逐漸降低,呈現(xiàn)光滑的拋物線下降趨勢(shì),但在應(yīng)變后期,出現(xiàn)了應(yīng)力隨應(yīng)變降低非常緩慢,之后出現(xiàn)斷崖式的下降,直至最后的試驗(yàn)結(jié)束。這種變化特點(diǎn)出現(xiàn)在所有膠結(jié)材料的加載試驗(yàn)得出的應(yīng)力應(yīng)變曲線中,說明這是最后破壞時(shí)期存在的共性特點(diǎn)。為了更好厘清產(chǎn)生的原因及所反映的力學(xué)變形特征,對(duì)曲線凸型段進(jìn)行了應(yīng)變能的計(jì)算。

    圖9(a)中No.1試樣的凸型段應(yīng)變能為6.319E~3MPa,(b)圖No.2試樣的凸型段應(yīng)變能為5.524E~3MPa,(c)圖No.1試樣的凸型段(刪除凸型段試驗(yàn)數(shù)據(jù)之后的擬合曲線段)應(yīng)變能為4.446E~3MPa,(d)圖No.2試樣2的凸型段(刪除凸型段試驗(yàn)數(shù)據(jù)之后的擬合曲線段)應(yīng)變能為0.0040407152111206MPa4.041E~3MPa,從而得出No.1試樣非趨勢(shì)變化所形成的應(yīng)變能為1.87E~3MPa,No.2試樣非趨勢(shì)變化所形成的應(yīng)變能為1.48E~3MPa。為此提出了如下設(shè)想,即無論是配比試樣,還是巖石試件,在加載過程中,能量大規(guī)模釋放的時(shí)間存在差異,脆性巖石存在加載破壞時(shí)的能量大量釋放,膠結(jié)材料存在能量釋放的滯后效應(yīng),即這部分能量是試樣在加載破壞后的完全整體破壞時(shí),其能量釋放所引起的體積膨脹導(dǎo)致的應(yīng)變能增加。

    2.3 噪音分貝能量測(cè)量方法

    2.3.1 噪音能量測(cè)量

    對(duì)膠結(jié)充填體進(jìn)行測(cè)試,分貝值在63.4~65.8dB之間變化(見圖10),最后在第二跌落區(qū)域的破壞分貝值達(dá)到79.4 dB。

    圖10 應(yīng)力和分貝隨著應(yīng)變變化曲線

    2.3.2 分貝能量測(cè)量試驗(yàn)

    對(duì)磷礦石進(jìn)行蠕變加載試驗(yàn),當(dāng)達(dá)到最大荷載158.876 kN時(shí),其破壞分貝值達(dá)到101.2 dB,而之前的分貝為76.6~78.9 dB,之后的分貝馬上降為78.4 dB。從破壞所發(fā)生的試件彈射體距離試件最遠(yuǎn)達(dá)到97.8 cm。

    為了驗(yàn)證這種設(shè)想,考慮到膠結(jié)材料自身強(qiáng)度較小,使用鐵礦石進(jìn)行相關(guān)加載,驗(yàn)證能量釋放設(shè)想。取礦山充填現(xiàn)場(chǎng)礦巖,經(jīng)過切割鉆孔打磨之后通過WAW-300微機(jī)控制電液伺服萬能試驗(yàn)機(jī)進(jìn)行軸向加載試驗(yàn)。通過AS804數(shù)字噪音計(jì)進(jìn)行礦巖加載過程中的分貝測(cè)量(見圖11),以期建立能量破壞與分貝之間的對(duì)應(yīng)關(guān)系,化定性分析為定量分析。為了測(cè)量的準(zhǔn)確性,測(cè)試儀器靠近試件,兩者相距15 cm,同時(shí)關(guān)閉門窗和其他試驗(yàn)設(shè)備。試驗(yàn)機(jī)未加載時(shí),實(shí)驗(yàn)室內(nèi)的分貝為30.3~31.1 dB。機(jī)器啟動(dòng)后,相應(yīng)的分貝在63.3~65.4 dB之間變化。每隔10 kN取值一次,相應(yīng)變化較大處,適當(dāng)密集取值。

    圖11 數(shù)字噪音計(jì)監(jiān)測(cè)

    2.3.2 試驗(yàn)結(jié)果及分析

    通過對(duì)礦巖的加載,相應(yīng)的分貝隨試驗(yàn)力的變化關(guān)系曲線見圖12。從圖12可以看出,荷載在50 kN之前,距離礦巖試件15 cm處的數(shù)字噪音計(jì)測(cè)得的分貝值基本沒有變化,但荷載達(dá)到60 kN時(shí),測(cè)量的分貝值出現(xiàn)了小幅增加,之后變化幅度較小,直到最后試樣發(fā)生破壞,分貝值達(dá)到了最高值93.2 dB。

    圖12 分貝試驗(yàn)力關(guān)系曲線

    為了驗(yàn)證礦巖試件能量積聚越多,在最后試件破壞時(shí)產(chǎn)生的分貝數(shù)越大。提出了一種蠕變加載能量積聚法,即將試件加載至某一恒定荷載值后,維持恒載一定時(shí)間,之后再繼續(xù)加載直至破壞。圖13是進(jìn)行蠕變加載能量積聚過程中的應(yīng)變和分貝隨時(shí)間的關(guān)系曲線。從圖13可以看出,在進(jìn)行蠕變能量積聚的試驗(yàn)過程中,軸向(加載方向)的蠕變值有增大的趨勢(shì),中間過程出現(xiàn)了多次應(yīng)變調(diào)整,說明礦巖的蠕變具有間歇式等速蠕變的變化規(guī)律,相應(yīng)的分貝值存在上下起伏變化,在蠕變能量積聚過程中并不具有相應(yīng)的變化規(guī)律。

    圖13 時(shí)間-應(yīng)變-分貝關(guān)系曲線

    考慮到試件在軸向加載破壞時(shí),同批次的最低單軸抗壓強(qiáng)度為114.088 kN,為此,本文以No.3試件為試驗(yàn)試件,將其加載至100 kN后進(jìn)行蠕變式恒載達(dá)到5400 s后再繼續(xù)加載。其破壞時(shí)的分貝值達(dá)到101.2 dB,之后回落到78.4 dB。破壞發(fā)生的巖爆分貝明顯高于未進(jìn)行蠕變時(shí)加載試件,接連的試件加載對(duì)比也驗(yàn)證了這一現(xiàn)象,有力的說明了礦巖試件在加載過程中積聚的能量越多,最后破壞時(shí)釋放的能量越多,相應(yīng)產(chǎn)生的分貝值越大。為了更好分析這種分貝變化產(chǎn)生的能量變化,建立分貝與功率之間的關(guān)系式,即:

    N=10lg(p/0) (11)

    式中:N表示分貝值,p和0為功率值,0是基準(zhǔn)功率量,一般為10~12 W。

    換算成功率關(guān)于分貝的函數(shù)關(guān)系式,即:

    p=10(Ndbi/10)0(12)

    顯然蠕變加載形成的能量值比未進(jìn)行蠕變加載積聚的能量高,相應(yīng)高出值為:

    △=10((101.2-93.2)/10)0=6.610(13)

    這部分高出的能量達(dá)到66.1~79.3 W之間。

    根據(jù)總能量守恒,試件出現(xiàn)的第二跌落區(qū)以及伴隨出現(xiàn)的噪音分貝的徒增所表現(xiàn)出的較大能量釋放,說明能量釋放階段相較常規(guī)加載出現(xiàn)了能量釋放的增加,從而出現(xiàn)試件在能量積聚階段積聚了更多能量,在緩慢加載過程中,主要體現(xiàn)在應(yīng)變能的增加。

    2.3.3 能量積聚與集中釋放條件

    通過試驗(yàn)結(jié)果分析,試樣要能夠積聚釋放,形成集中釋放情況,需要滿足如下三個(gè)條件:

    (一)試樣加載速度非常小,等同于局部加載階段,呈現(xiàn)類靜壓過程。

    (二)試樣是彈塑性體,試樣需要加載至塑性區(qū)。即彈性階段的加載速度不受一個(gè)條件限制,試樣能量積聚并能夠集中釋放的階段不一。對(duì)于強(qiáng)度較低的充填體,其能量積聚后的能量釋放是在充填體靠近表壁部位發(fā)生受拉破壞后,內(nèi)部的積聚能量再次釋放。而對(duì)于強(qiáng)度較高的充填體,則是能量不斷積聚,直到試件達(dá)到最大抗壓強(qiáng)度之時(shí)突然破壞。

    (三)試樣需要有能量集中釋放的臨空面。如果試樣沒有臨空面,試樣在加載過程中的能量積聚會(huì)通過傳導(dǎo)介質(zhì)傳遞出去,而臨空面是隔絕試樣能量多角度傳遞,創(chuàng)造能量集中釋放條件。

    在能量積聚期間的塑形變形階段,試件不斷積累破壞,在破壞的過程中,保持荷載不變,使試件進(jìn)行自我的應(yīng)變協(xié)調(diào),趨向自穩(wěn)方向發(fā)展,從而相對(duì)原先,試件內(nèi)部形成更為穩(wěn)定的致密體,提高了抵御外界加載破壞的能力,有利于將外界的做功以應(yīng)變能的形式儲(chǔ)存在試件內(nèi)部。當(dāng)試件處在能量釋放階段時(shí),由于內(nèi)部形成的更為致密的區(qū)域,抵御破壞發(fā)生的能力更強(qiáng),在荷載超過其承受能力時(shí),將以突然釋放的形式展現(xiàn)。

    3 結(jié)論

    1)膠結(jié)材料加載破壞前存在多次應(yīng)變協(xié)調(diào)階段,加載破壞后應(yīng)力應(yīng)變曲線存在的凸型特征是膠結(jié)材料整體破壞時(shí)的能量釋放。對(duì)膠結(jié)材料的應(yīng)力應(yīng)變曲線進(jìn)行應(yīng)變能和即時(shí)彈性模量的計(jì)算,引入柔化正態(tài)張量和應(yīng)變能張量,建立兩者之間的對(duì)應(yīng)關(guān)系,在應(yīng)變能累積值較小時(shí)其曲線呈現(xiàn)水平變化,累積值將要達(dá)到最大時(shí)呈現(xiàn)近豎直變化,中間階段呈近似的線性變化特征。

    2)提出了一種蠕變加載能量積聚法,為驗(yàn)證試樣加載破壞后的能量釋放大小是依托能量吸收為前提條件,特通過對(duì)礦巖進(jìn)行蠕變恒載保持之后的加載破壞,并首次提出引入數(shù)字噪音計(jì)進(jìn)行試件破壞分貝的測(cè)量。按照預(yù)期設(shè)想,試件進(jìn)行蠕變加載能量積聚的分貝較未蠕變加載高出近20-35 dB,說明這種蠕變加載能量積聚法的可行以及能量釋放的強(qiáng)度是以前期積聚能量為前提的。

    3)試件進(jìn)入塑性區(qū)后,在其荷載是試件平均荷載的85%左右時(shí)開始進(jìn)行蠕變加載,加載時(shí)間為1 h,是否存在在不同塑性荷載階段加載,以及蠕變加載時(shí)間不同,對(duì)其能量釋放造成的影響是否具有對(duì)應(yīng)關(guān)系,值得繼續(xù)研究。

    [1] 何滿潮,劉冬橋,宮偉力,等. 沖擊巖爆試驗(yàn)系統(tǒng)研發(fā)及試驗(yàn)[J].巖石力學(xué)與工程學(xué)報(bào),2014,33(9):1729-1739.

    [2] 夏元友,吝曼卿,廖璐璐,等. 大尺寸試件巖爆試驗(yàn)碎屑分形特征分析[J].巖石力學(xué)與工程學(xué)報(bào),2014,33(7): 1358-1365.

    [3] 龔劍,胡乃聯(lián),崔翔,等. 基于AHP-TOPSIS評(píng)判模型的巖爆傾向性預(yù)測(cè)[J].巖石力學(xué)與工程學(xué)報(bào),2014,33(7): 1442-1448.

    [4] 蔡美峰,冀東,郭奇峰. 基于地應(yīng)力現(xiàn)場(chǎng)實(shí)測(cè)與開采擾動(dòng)能量積聚理論的巖爆預(yù)測(cè)研究[J].巖石力學(xué)與工程學(xué)報(bào),2013,32(10):1973-1980.

    [5] 賈義鵬,呂慶,尚岳全. 基于粒子群算法和廣義回歸神經(jīng)網(wǎng)絡(luò)的巖爆預(yù)測(cè)[J].巖石力學(xué)與工程學(xué)報(bào),2013,32(2): 343-348.

    [6] 張文東,馬天輝,唐春安,等. 錦屏二級(jí)水電站引水隧洞巖爆特征及微震監(jiān)測(cè)規(guī)律研究[J].巖石力學(xué)與工程學(xué)報(bào),2014,33(2):339-348.

    [7] 于洋,馮夏庭,陳天宇,等. 深埋隧洞不同掘進(jìn)方式下即時(shí)型巖爆微震對(duì)比分析[J].東北大學(xué)學(xué)報(bào):自然科學(xué)版, 2014,35(3):429-432.

    [8] Zhou K P, Lin Y, Deng H W, et al. Prediction of rock burst classification using cloud model with entropy weight[J]. Trans. Nonferrous Met. Soc. China,2016(26):1995?2002.

    [9] Zhou J, Shi X Z, Huang R D, et al. Feasibility of stochastic gradient boosting approach for predicting rockburst damage in burst-prone mines[J]. Trans. Nonferrous Met. Soc. China,2016(26):1938?1945.

    [10] Dong L J, Li X B, Peng K, et al. Prediction of rockburst classification using random forest[J]. Trans.Nonferrous Met.Soc.China, 2013(23):472-477.

    [11] 趙康,趙紅宇,賈群燕,等. 巖爆巖石斷裂的微觀結(jié)構(gòu)形貌分析及巖爆機(jī)理[J].爆炸與沖擊,2015,35(6):913-918.

    [12] 尚彥軍,張鏡劍,傅冰駿,等. 應(yīng)變型巖爆三要素分析及巖爆勢(shì)表達(dá)[J].巖石力學(xué)與工程學(xué)報(bào),2013,32(8):1520-1527.

    [13] 馬春馳,李天斌,陳國(guó)慶,等. 硬脆巖石的微觀顆粒模型及其卸荷巖爆效應(yīng)研究[J].巖石力學(xué)與工程學(xué)報(bào), 2015,34(2):217-227.

    [14] 蘇國(guó)韶,陳智勇,蔣劍青,等. 不同加載速率下巖爆碎塊耗能特征試驗(yàn)研究[J].巖土工程學(xué)報(bào), 2016,38(8): 1481-1489.

    [15] 姜春露,孫強(qiáng),姜振泉,等. 脆性巖石全應(yīng)力應(yīng)變過程中滲透性突變研究[J].中南大學(xué)學(xué)報(bào):自然科學(xué)版,2012,43(2): 688-693.

    [16] 鄧林,武君,呂燕,等. 基于巖石應(yīng)力應(yīng)變過程曲線的巖爆能量指數(shù)法[J].鐵道標(biāo)準(zhǔn)設(shè)計(jì),2012(7):108-111.

    [17] 于德海,彭建兵,崔春義,等. 片狀巖石全應(yīng)力應(yīng)變壓縮過程及力學(xué)機(jī)制試驗(yàn)研究[J].應(yīng)用基礎(chǔ)與工程科學(xué)學(xué)報(bào),2010,18(5):792-800.

    [18] 姜永東,鮮學(xué)福,尹光志,等. 巖石應(yīng)力應(yīng)變?nèi)^程的聲發(fā)射及分形與混沌特征[J].巖土力學(xué),2010,31(8):2413-2418.

    [19] 楊天雨,喬登攀. 循環(huán)加卸載下尾砂膠結(jié)充填體能量變化特征的試驗(yàn)研究與分析[J].價(jià)值工程,2017(3):127-130.

    [20] 楊逾,王冰芬,鄭志明. 煤礦矸石充填材料壓實(shí)耗能特性試驗(yàn)研究[J].硅酸鹽通報(bào),2016,35(11):3511-3516.

    [21] 徐文彬,宋衛(wèi)東,王東旭,等. 膠結(jié)充填體三軸壓縮變形破壞及能量耗散特征分析[J].巖土力學(xué),2014,35(12): 3421-3429.

    [22] 劉玉龍,丁德馨,李廣悅,等. 膠結(jié)充填體固化與開挖礦體垂直應(yīng)力匹配特征[J].采礦與安全工程學(xué)報(bào),2013, 30(4): 526-530.

    [23] 劉智高. 基于能量原理的沖擊地壓防治關(guān)鍵層確定實(shí)驗(yàn)研究[D].哈爾濱:黑龍江科技大學(xué),2014:25.

    [24] 趙忠虎,謝和平. 巖石變形破壞過程中的能量傳遞和耗散研究[J].四川大學(xué)學(xué)報(bào):工程科學(xué)版,2008, 40(2):26-31.

    Studyon energy accumulation and release response of cemented materials subjected to loading failure

    WANG Ming-xu1,2

    (1. School of Resource and Environmental Engineering,Wuhan University of Science and Technology, Wuhan, Hubei 430081,China;2. Huizhou Land Resources Bureau Huizhou,,Guangdong 516003,China)

    To study the second drop characteristic that occurs after the peak strength of the cemented material at a rate of 0.001 kN / s, and the occurrence of the sudden increase in the noise decibel value occurred in the second drop zone, the application of energy accumulation. The total stress-strain curve was obtained by the axial loading test of the cemented filling body, and the size of the noise decibel value during the loading process was tested by the electronic noise test equipment. The real elasticity of the stress-strain curve,the ratio of the modulus of elasticity to the modulus of elasticity and the ratio of the total strain energy as the strain tensor could be used to soften the normal tensor and strain the size of the changes in the characterization of the specimen in the energy accumulation stage and the release of the response stage of the characteristics of change, especially in the second drop stage of the energy characteristics. At the same time, the corresponding relationship between the noise decibel and the energy release of the specimen was established by the noise decibel measurement method. The results showed that only the loading speed reaches about 0.001kN / s, the cementing body would appear obvious second drop area. Through the electronic noise decibel test the second drop zone damage noise generated by the sudden increase of 15dB, compared to the conventional load increased by more than 10dB. At the same time, this slow loading method was applied to the rock and rock loading test, and the noise decibel value of the loading was more obvious and the added value is 20-35dB. It was shown that the strength of the energy release was prerequisite for the accumulation of energy in the early period, which would help to guide the mining of metal mining for high stress area evaluation, maintenance stope safety.

    filling; energy release; elastic modulus; strain energy; noise energy

    1674-8085(2020)01-0065-10

    TU 853

    A

    10.3969/j.issn.1674-8085.2020.01.013

    2018-12-26;

    2019-09-21

    國(guó)家自然科學(xué)基金項(xiàng)目(51604195)

    王明旭(1986-),男,湖北監(jiān)利人,博士生,主要從事礦業(yè)工程,充填體與圍名相互作用研究(E-mail: 1765207689@qq.com).

    猜你喜歡
    巖爆分貝張量
    某引水隧洞深埋段花崗巖巖爆演化規(guī)律研究
    使用聲發(fā)射測(cè)試結(jié)果預(yù)判巖爆等級(jí)
    金屬礦山(2022年1期)2022-02-23 11:16:36
    偶數(shù)階張量core逆的性質(zhì)和應(yīng)用
    四元數(shù)張量方程A*NX=B 的通解
    分貝是如何定義的
    分貝是聲波的強(qiáng)度單位嗎
    擴(kuò)散張量成像MRI 在CO中毒后遲發(fā)腦病中的應(yīng)用
    引漢濟(jì)渭工程秦嶺隧洞巖爆數(shù)值模擬與巖爆預(yù)測(cè)研究
    深部開采中巖爆巖塊彈射速度的理論與實(shí)驗(yàn)
    小分貝上街“不講話”了
    一级毛片电影观看| 久久久亚洲精品成人影院| 久久人人97超碰香蕉20202| 黄色毛片三级朝国网站| 亚洲欧洲国产日韩| 寂寞人妻少妇视频99o| 成人亚洲精品一区在线观看| 国产精品国产av在线观看| 欧美av亚洲av综合av国产av | 欧美xxⅹ黑人| 国产日韩欧美亚洲二区| 777久久人妻少妇嫩草av网站| 91精品国产国语对白视频| 一区二区av电影网| 岛国毛片在线播放| 女人高潮潮喷娇喘18禁视频| 妹子高潮喷水视频| 黄色视频在线播放观看不卡| 亚洲情色 制服丝袜| 又粗又硬又长又爽又黄的视频| 日本色播在线视频| 国产欧美亚洲国产| 精品一区二区三卡| 亚洲欧美中文字幕日韩二区| 欧美日韩视频高清一区二区三区二| 中文字幕色久视频| 26uuu在线亚洲综合色| 精品国产露脸久久av麻豆| 在线观看人妻少妇| 久久国内精品自在自线图片| 日韩中文字幕视频在线看片| 亚洲av中文av极速乱| 亚洲精品在线美女| 嫩草影院入口| 啦啦啦在线观看免费高清www| 婷婷色麻豆天堂久久| 90打野战视频偷拍视频| 久久精品aⅴ一区二区三区四区 | 香蕉精品网在线| 97精品久久久久久久久久精品| 99久国产av精品国产电影| 王馨瑶露胸无遮挡在线观看| 国精品久久久久久国模美| 亚洲人成电影观看| 亚洲成人一二三区av| 男人爽女人下面视频在线观看| 日日爽夜夜爽网站| 免费观看性生交大片5| 丰满迷人的少妇在线观看| 欧美中文综合在线视频| 成人手机av| 精品国产露脸久久av麻豆| 久久久欧美国产精品| 一区二区日韩欧美中文字幕| 国产成人欧美| 午夜激情av网站| 亚洲四区av| 一二三四中文在线观看免费高清| 男女高潮啪啪啪动态图| √禁漫天堂资源中文www| 看非洲黑人一级黄片| 1024视频免费在线观看| 满18在线观看网站| 中文字幕另类日韩欧美亚洲嫩草| 午夜影院在线不卡| 蜜桃国产av成人99| 午夜福利影视在线免费观看| 成人影院久久| 色吧在线观看| 亚洲人成电影观看| 香蕉丝袜av| 国产探花极品一区二区| 美女午夜性视频免费| 亚洲精品国产av蜜桃| 欧美人与性动交α欧美软件| 日韩成人av中文字幕在线观看| 国产亚洲最大av| 国产片内射在线| 国产视频首页在线观看| 久久人人爽av亚洲精品天堂| 多毛熟女@视频| 国产精品秋霞免费鲁丝片| 成人毛片60女人毛片免费| 免费观看在线日韩| 久久精品国产亚洲av高清一级| 18禁国产床啪视频网站| 亚洲伊人色综图| 黄色 视频免费看| 亚洲成人一二三区av| 精品少妇一区二区三区视频日本电影 | 观看美女的网站| 爱豆传媒免费全集在线观看| 亚洲国产av影院在线观看| 青草久久国产| 1024视频免费在线观看| 亚洲人成电影观看| 国产黄色免费在线视频| 国产精品偷伦视频观看了| 国产女主播在线喷水免费视频网站| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 极品少妇高潮喷水抽搐| 建设人人有责人人尽责人人享有的| 精品国产国语对白av| 精品第一国产精品| 桃花免费在线播放| 亚洲色图 男人天堂 中文字幕| 中文天堂在线官网| 欧美最新免费一区二区三区| 亚洲久久久国产精品| 青春草视频在线免费观看| 亚洲av在线观看美女高潮| 妹子高潮喷水视频| 欧美精品亚洲一区二区| 一本大道久久a久久精品| av国产精品久久久久影院| 亚洲av福利一区| 国产精品一国产av| 高清视频免费观看一区二区| 国产探花极品一区二区| 男女午夜视频在线观看| 狠狠精品人妻久久久久久综合| 亚洲欧美一区二区三区久久| 午夜日韩欧美国产| 欧美国产精品va在线观看不卡| 午夜激情av网站| av电影中文网址| 日日爽夜夜爽网站| 亚洲色图综合在线观看| 热re99久久国产66热| 婷婷色综合www| 亚洲av电影在线观看一区二区三区| 亚洲精品乱久久久久久| 女人久久www免费人成看片| 国产不卡av网站在线观看| 欧美黄色片欧美黄色片| 国产不卡av网站在线观看| 蜜桃在线观看..| 午夜免费观看性视频| 国产成人午夜福利电影在线观看| 黄色一级大片看看| 精品人妻熟女毛片av久久网站| 老汉色av国产亚洲站长工具| 岛国毛片在线播放| 免费在线观看完整版高清| 日韩,欧美,国产一区二区三区| 国产精品一区二区在线观看99| 最近最新中文字幕大全免费视频 | 日韩中文字幕视频在线看片| 少妇人妻久久综合中文| 国产极品粉嫩免费观看在线| 免费人妻精品一区二区三区视频| 国产精品 国内视频| 麻豆av在线久日| 永久网站在线| 精品福利永久在线观看| 亚洲人成电影观看| 少妇人妻久久综合中文| av网站在线播放免费| 青春草亚洲视频在线观看| 久久青草综合色| 黑人欧美特级aaaaaa片| 我的亚洲天堂| 亚洲色图综合在线观看| 一级毛片我不卡| 美女国产高潮福利片在线看| 免费日韩欧美在线观看| 99九九在线精品视频| 欧美日韩国产mv在线观看视频| 国产免费视频播放在线视频| 五月天丁香电影| 国产黄色视频一区二区在线观看| 日韩中文字幕欧美一区二区 | 伦理电影大哥的女人| 国产免费现黄频在线看| 在线天堂最新版资源| 国产高清不卡午夜福利| 亚洲av欧美aⅴ国产| 国产精品无大码| 亚洲成国产人片在线观看| 久久热在线av| 看免费av毛片| 久久久久久久久免费视频了| 亚洲国产精品一区三区| 国产精品一区二区在线不卡| 黑人猛操日本美女一级片| 日本av手机在线免费观看| 欧美精品国产亚洲| 2018国产大陆天天弄谢| 在线观看国产h片| 亚洲成人手机| 国产黄色视频一区二区在线观看| 久久久久久久大尺度免费视频| 男人添女人高潮全过程视频| 美女中出高潮动态图| 成人午夜精彩视频在线观看| 国产成人一区二区在线| 人体艺术视频欧美日本| 青草久久国产| 女人精品久久久久毛片| 成年女人毛片免费观看观看9 | 国产欧美亚洲国产| 高清视频免费观看一区二区| 老司机影院毛片| 边亲边吃奶的免费视频| 1024视频免费在线观看| 亚洲精品,欧美精品| 亚洲国产av影院在线观看| 久久这里只有精品19| www.精华液| 最近中文字幕高清免费大全6| 色婷婷久久久亚洲欧美| 欧美亚洲日本最大视频资源| 久久精品国产亚洲av天美| 久久久久久久国产电影| 美女视频免费永久观看网站| av不卡在线播放| 男人爽女人下面视频在线观看| 侵犯人妻中文字幕一二三四区| 妹子高潮喷水视频| 亚洲美女搞黄在线观看| av在线app专区| 国产日韩欧美在线精品| 少妇被粗大的猛进出69影院| 亚洲在久久综合| 国产高清不卡午夜福利| 大话2 男鬼变身卡| freevideosex欧美| 男人操女人黄网站| 精品国产一区二区久久| 国产成人精品婷婷| 日韩精品免费视频一区二区三区| 国产一区有黄有色的免费视频| 国产xxxxx性猛交| 欧美变态另类bdsm刘玥| 九草在线视频观看| 国产熟女欧美一区二区| 91午夜精品亚洲一区二区三区| 成人黄色视频免费在线看| 性少妇av在线| 激情五月婷婷亚洲| 2021少妇久久久久久久久久久| 成人二区视频| 国产xxxxx性猛交| 久久久久国产网址| 2022亚洲国产成人精品| av网站在线播放免费| 亚洲婷婷狠狠爱综合网| 婷婷成人精品国产| 国产精品国产av在线观看| 国产女主播在线喷水免费视频网站| 哪个播放器可以免费观看大片| 欧美日韩亚洲国产一区二区在线观看 | 国产精品一二三区在线看| 久久精品国产a三级三级三级| 成人午夜精彩视频在线观看| av有码第一页| 久久免费观看电影| 菩萨蛮人人尽说江南好唐韦庄| 中文精品一卡2卡3卡4更新| 久久精品久久精品一区二区三区| 午夜激情av网站| 国产一级毛片在线| 人妻 亚洲 视频| 一级毛片黄色毛片免费观看视频| 欧美日韩成人在线一区二区| av片东京热男人的天堂| 在线观看国产h片| 国产精品熟女久久久久浪| 国产欧美日韩一区二区三区在线| a级片在线免费高清观看视频| 色94色欧美一区二区| 国产女主播在线喷水免费视频网站| 91国产中文字幕| 18禁裸乳无遮挡动漫免费视频| 黑人欧美特级aaaaaa片| 日韩免费高清中文字幕av| 建设人人有责人人尽责人人享有的| 乱人伦中国视频| 1024香蕉在线观看| 人体艺术视频欧美日本| 久久99蜜桃精品久久| 日韩三级伦理在线观看| 天堂8中文在线网| 中文字幕精品免费在线观看视频| 日韩中文字幕欧美一区二区 | 免费人妻精品一区二区三区视频| 亚洲综合色惰| 美女xxoo啪啪120秒动态图| 国产在视频线精品| 热re99久久精品国产66热6| 美女中出高潮动态图| 日韩,欧美,国产一区二区三区| 欧美xxⅹ黑人| 日韩伦理黄色片| 制服人妻中文乱码| 免费黄色在线免费观看| 日本91视频免费播放| 女性生殖器流出的白浆| 欧美人与性动交α欧美精品济南到 | 在线看a的网站| 欧美日韩一级在线毛片| 亚洲色图综合在线观看| 激情视频va一区二区三区| 成人午夜精彩视频在线观看| 国产白丝娇喘喷水9色精品| 人成视频在线观看免费观看| 日本vs欧美在线观看视频| 人人妻人人添人人爽欧美一区卜| 亚洲五月色婷婷综合| 另类亚洲欧美激情| 欧美日韩成人在线一区二区| 欧美精品av麻豆av| 黄网站色视频无遮挡免费观看| 国产成人一区二区在线| 日本av免费视频播放| 搡老乐熟女国产| 欧美人与性动交α欧美软件| 日日啪夜夜爽| 精品国产一区二区三区久久久樱花| 久久99蜜桃精品久久| 两个人看的免费小视频| 97人妻天天添夜夜摸| 成年动漫av网址| 天天躁夜夜躁狠狠久久av| 欧美精品国产亚洲| 毛片一级片免费看久久久久| 久久久精品区二区三区| 成年女人在线观看亚洲视频| 国产高清不卡午夜福利| 日本色播在线视频| 精品一区在线观看国产| 99久国产av精品国产电影| 三上悠亚av全集在线观看| 欧美老熟妇乱子伦牲交| 欧美xxⅹ黑人| 美女国产视频在线观看| 国产av国产精品国产| 麻豆精品久久久久久蜜桃| 亚洲精品久久成人aⅴ小说| 熟女电影av网| 国产麻豆69| 寂寞人妻少妇视频99o| 26uuu在线亚洲综合色| 欧美少妇被猛烈插入视频| 亚洲综合精品二区| 久久人人97超碰香蕉20202| 亚洲成av片中文字幕在线观看 | 欧美bdsm另类| 综合色丁香网| 久久久久国产精品人妻一区二区| 国产色婷婷99| 黑丝袜美女国产一区| 免费播放大片免费观看视频在线观看| 人人澡人人妻人| 99热全是精品| 午夜福利一区二区在线看| 亚洲一码二码三码区别大吗| 一本色道久久久久久精品综合| 精品卡一卡二卡四卡免费| 在线亚洲精品国产二区图片欧美| 欧美激情 高清一区二区三区| 日韩在线高清观看一区二区三区| 最近最新中文字幕大全免费视频 | 国产在视频线精品| 男女午夜视频在线观看| 美女脱内裤让男人舔精品视频| 国产麻豆69| 亚洲精品日本国产第一区| 18+在线观看网站| 日本午夜av视频| 一区二区三区乱码不卡18| 久久狼人影院| 大香蕉久久成人网| 久久韩国三级中文字幕| 欧美成人午夜免费资源| 久久久国产精品麻豆| 亚洲国产欧美在线一区| av不卡在线播放| 18+在线观看网站| 欧美人与性动交α欧美精品济南到 | 久久久精品国产亚洲av高清涩受| 久久人人爽av亚洲精品天堂| 少妇人妻精品综合一区二区| 国产精品不卡视频一区二区| 国产精品无大码| 午夜福利视频精品| 高清在线视频一区二区三区| 日本-黄色视频高清免费观看| 国产激情久久老熟女| 久久精品久久久久久久性| videosex国产| 国产激情久久老熟女| 亚洲成人一二三区av| 久久久久国产网址| 99热国产这里只有精品6| 亚洲精品国产av成人精品| 日韩一区二区视频免费看| 一级,二级,三级黄色视频| 久久热在线av| 亚洲国产成人一精品久久久| 狠狠婷婷综合久久久久久88av| 2018国产大陆天天弄谢| 亚洲成人一二三区av| 国产在线免费精品| av福利片在线| 美女福利国产在线| 97在线人人人人妻| 欧美激情高清一区二区三区 | 黄片无遮挡物在线观看| 黄色怎么调成土黄色| 美女视频免费永久观看网站| 岛国毛片在线播放| 精品一区二区三区四区五区乱码 | 麻豆乱淫一区二区| 99国产精品免费福利视频| 精品国产乱码久久久久久小说| 国产精品久久久久久av不卡| 国产野战对白在线观看| 制服丝袜香蕉在线| 国产精品 欧美亚洲| 丰满饥渴人妻一区二区三| 欧美日本中文国产一区发布| 最新中文字幕久久久久| 午夜91福利影院| 免费黄频网站在线观看国产| 亚洲一码二码三码区别大吗| 男女国产视频网站| 飞空精品影院首页| 国产成人精品在线电影| 免费不卡的大黄色大毛片视频在线观看| 最近2019中文字幕mv第一页| 国产白丝娇喘喷水9色精品| 国产精品三级大全| 熟妇人妻不卡中文字幕| 美女大奶头黄色视频| 亚洲 欧美一区二区三区| 精品国产一区二区三区四区第35| 91aial.com中文字幕在线观看| 色婷婷久久久亚洲欧美| 国产成人一区二区在线| 午夜av观看不卡| 午夜日韩欧美国产| 国产精品久久久av美女十八| 欧美精品一区二区免费开放| 伦理电影大哥的女人| 精品少妇黑人巨大在线播放| 亚洲国产av影院在线观看| 丝袜脚勾引网站| 亚洲国产日韩一区二区| 中国三级夫妇交换| 国产精品免费视频内射| 人人妻人人爽人人添夜夜欢视频| 日本wwww免费看| 一级a爱视频在线免费观看| 国产免费现黄频在线看| 91久久精品国产一区二区三区| 国产乱人偷精品视频| 精品人妻偷拍中文字幕| 色婷婷av一区二区三区视频| 亚洲国产精品成人久久小说| 久久国产精品大桥未久av| 十分钟在线观看高清视频www| 久久久久久久久免费视频了| 久久国产精品男人的天堂亚洲| 这个男人来自地球电影免费观看 | 人人妻人人添人人爽欧美一区卜| 久久久久久久亚洲中文字幕| 男人操女人黄网站| 亚洲欧美中文字幕日韩二区| 久久国产精品大桥未久av| 亚洲精品久久午夜乱码| 18+在线观看网站| 亚洲精品国产一区二区精华液| 国产日韩欧美亚洲二区| 久久久久精品久久久久真实原创| 天堂俺去俺来也www色官网| 久久久久久久久久久久大奶| 亚洲国产色片| 亚洲国产精品一区二区三区在线| 性色avwww在线观看| 国产免费现黄频在线看| 一区二区三区四区激情视频| 久久精品国产亚洲av涩爱| 搡老乐熟女国产| 亚洲精品一二三| 亚洲欧洲精品一区二区精品久久久 | 亚洲国产av新网站| 国产爽快片一区二区三区| 国产 一区精品| 人妻少妇偷人精品九色| 亚洲欧美精品自产自拍| 美女福利国产在线| 久久久欧美国产精品| 18禁裸乳无遮挡动漫免费视频| 超碰97精品在线观看| 亚洲精品久久久久久婷婷小说| 最近2019中文字幕mv第一页| 日韩一本色道免费dvd| 成年人午夜在线观看视频| 欧美激情高清一区二区三区 | 久久精品国产综合久久久| 考比视频在线观看| 一本色道久久久久久精品综合| 女的被弄到高潮叫床怎么办| 国产精品国产三级专区第一集| 不卡av一区二区三区| 国产男女超爽视频在线观看| 观看美女的网站| 黄片小视频在线播放| 久久久久国产一级毛片高清牌| 久久久久久人妻| 欧美日韩视频精品一区| 在线免费观看不下载黄p国产| 黄频高清免费视频| 热re99久久国产66热| 亚洲在久久综合| 欧美另类一区| 午夜免费观看性视频| 最近最新中文字幕免费大全7| 18禁动态无遮挡网站| 亚洲一码二码三码区别大吗| 久久99热这里只频精品6学生| 两个人免费观看高清视频| av福利片在线| 欧美日韩精品网址| 成人亚洲精品一区在线观看| 亚洲欧洲日产国产| 国产精品久久久av美女十八| 国产黄色免费在线视频| 精品国产国语对白av| 精品少妇黑人巨大在线播放| 精品久久蜜臀av无| 亚洲图色成人| 精品国产乱码久久久久久小说| 亚洲成人一二三区av| 久久99精品国语久久久| 亚洲欧美一区二区三区国产| 一级毛片黄色毛片免费观看视频| 丝袜美足系列| www.av在线官网国产| 天美传媒精品一区二区| 国产精品国产av在线观看| 久久久久国产一级毛片高清牌| 国产亚洲精品第一综合不卡| 菩萨蛮人人尽说江南好唐韦庄| 妹子高潮喷水视频| 2022亚洲国产成人精品| 九色亚洲精品在线播放| 久久99一区二区三区| 免费在线观看视频国产中文字幕亚洲 | 丰满饥渴人妻一区二区三| 人成视频在线观看免费观看| 一二三四在线观看免费中文在| 汤姆久久久久久久影院中文字幕| 一级黄片播放器| 国产精品成人在线| 18禁动态无遮挡网站| 又大又黄又爽视频免费| www.自偷自拍.com| 亚洲av男天堂| 一级爰片在线观看| 成年女人毛片免费观看观看9 | 中文字幕亚洲精品专区| 青草久久国产| 久久人人爽人人片av| a级毛片在线看网站| 亚洲国产精品成人久久小说| 久久人妻熟女aⅴ| 国语对白做爰xxxⅹ性视频网站| 免费大片黄手机在线观看| 天美传媒精品一区二区| 人人妻人人爽人人添夜夜欢视频| 91久久精品国产一区二区三区| 考比视频在线观看| 精品一区在线观看国产| 亚洲欧美成人精品一区二区| 久久午夜综合久久蜜桃| 亚洲男人天堂网一区| av天堂久久9| 午夜福利,免费看| 国产日韩一区二区三区精品不卡| 国产精品不卡视频一区二区| 国产人伦9x9x在线观看 | 国产片内射在线| 宅男免费午夜| 人人妻人人添人人爽欧美一区卜| 日本黄色日本黄色录像| 毛片一级片免费看久久久久| 香蕉丝袜av| 久久精品夜色国产| 香蕉国产在线看| 九九爱精品视频在线观看| 国产一区二区在线观看av| 精品人妻一区二区三区麻豆| 亚洲综合色惰| 91成人精品电影| 又粗又硬又长又爽又黄的视频| 国产色婷婷99| 成人国语在线视频| 欧美97在线视频| 又粗又硬又长又爽又黄的视频| www.自偷自拍.com| 久久99蜜桃精品久久| 婷婷成人精品国产| 欧美日韩一级在线毛片| 日韩制服丝袜自拍偷拍| 国产成人免费观看mmmm| av在线老鸭窝| 大码成人一级视频| 精品国产一区二区三区四区第35| 夫妻午夜视频| 国产成人一区二区在线| 色视频在线一区二区三区| 欧美激情极品国产一区二区三区| 啦啦啦中文免费视频观看日本|