蔡改貧，郝書灝，余 成，宣律偉

（1. 江西理工大學(xué)機電工程學(xué)院, 江西 贛州 341000；2. 江西省礦冶機電工程技術(shù)研究中心, 江西 贛州 341000）

礦石是由金屬元素礦物和非金屬元素雜質(zhì)經(jīng)過復(fù)雜地質(zhì)作用形成的集合體，具有一定的結(jié)構(gòu)特征，其內(nèi)部存在不同尺度的微結(jié)構(gòu)面，如礦物的節(jié)理、晶粒邊界、粒間空隙、微裂隙等。由于存在自身結(jié)構(gòu)缺陷，且其本身又是由多種礦物顆粒、微小顆粒、層結(jié)構(gòu)、膠結(jié)物等組成的復(fù)雜集合體，因此巖石在變形和破壞過程中不同于連續(xù)介質(zhì)，呈現(xiàn)各向異性特征[1]。

礦石具有多種礦物組分，但有用礦物與大量黏土、石英、碳酸鹽晶體等脈石礦物緊密連結(jié)。解離是指借粉碎作用，使共生組分以單體形式脫離。Gaudin[2]較早對礦物解離問題進(jìn)行了系統(tǒng)性研究，并提出用解離度定量描述礦物解離；王雅蓉等[3]提出了界面形態(tài)特征對礦物單體解離的影響遵循相界面積守恒定律；劉建遠(yuǎn)等[4]從碎后粒度分布、銅金屬粒級分布、銅礦物粒度和顆粒含量分布、銅礦物粒級解離度等方面研究了載荷強度和粉碎方式對銅礦物解離的影響。

不同礦物組分之間的物理力學(xué)性能存在差異，則彼此間黏結(jié)形成的礦物界面的力學(xué)性能亦存在差異，這對礦石單一組分礦物解離具有重要影響。Sun 等[5]通過單變量和多變量回歸方程分析了巖石材料的結(jié)構(gòu)特征受礦物組成的影響，并全面總結(jié)了力學(xué)性能與礦物含量之間的回歸方程。楊小彬等[6]進(jìn)行了花崗巖試件在恒定圍壓下的雙剪模型試驗，分析了花崗巖試件在循環(huán)加卸載條件下的界面位移場、摩擦滑移位移的演化特征。Alneasan 等[7]詳細(xì)研究了壓縮和拉伸作用下兩個不同巖石層之間傾斜內(nèi)部界面裂紋情況，以及裂紋傾角對應(yīng)力強度因子和界面裂紋尖端附近應(yīng)力分布函數(shù)的影響。

礦石是長時間地質(zhì)作用產(chǎn)生的內(nèi)部結(jié)構(gòu)復(fù)雜的聚合物[8]，這使得礦石結(jié)構(gòu)組分及其物理特性非常復(fù)雜[9]。目前，國內(nèi)外學(xué)者普遍通過接觸理論、分形統(tǒng)計理論、細(xì)觀力學(xué)等對顆粒的宏觀或細(xì)觀破碎相關(guān)性質(zhì)進(jìn)行研究[10-11]?；陔x散元理論，相關(guān)學(xué)者將組成礦石的顆粒簡化成球狀顆粒模型，顆粒之間通過“黏結(jié)鍵”連接，從而模擬礦石內(nèi)部顆粒之間的黏結(jié)力。

綜上可知，目前關(guān)于不同組分礦物解離及礦物組分對礦石整體力學(xué)性質(zhì)的影響機制的研究較多，為研究巖石的破碎機理、提高巖石的破碎效率提供了一定的理論指導(dǎo)，但對組分礦物黏結(jié)界面等方面的研究較少，無法實現(xiàn)利用破碎特征對礦物界面進(jìn)行分類，從而無法為有用礦物高效解離提供指導(dǎo)。本研究將基于礦石的巖相分析，建立符合礦石特性的代表性體積單元(representative volume element，RVE)模型，通過模擬礦物界面原位加載實驗，并結(jié)合破碎實驗，對不同組分礦物界面破碎特征進(jìn)行分類，為研究礦石的破碎機理和提高礦石有用礦物的高效快速解離提供一定的理論依據(jù)。

1 礦石巖相分析

1.1 巖相采集設(shè)備和礦石預(yù)處理

巖相采集設(shè)備為配套高性能感光芯片傳感器和測量軟件的數(shù)碼顯微鏡，放大倍數(shù)為50～500倍；其內(nèi)置芯片以靜態(tài)拍攝為主，靜態(tài)分辨率為640×480，最大分辨率為1 600×1 200，其影像分辨率等同靜態(tài)分辨率。

所研究的礦石采集自贛州某黑鎢礦石，其主要有用礦物為金屬鎢，主要脈石礦物為石英和硅質(zhì)巖，如圖1 所示。為獲取礦石內(nèi)部巖相和組分界面區(qū)域巖相，對礦石進(jìn)行隨機切割和界面破壞，切割線如圖1中紅線所示。

圖1 處理前黑鎢礦石試樣Fig. 1 Wolframite ore sample before treatment

1.2 鎢顆粒區(qū)域分布

黑鎢礦石中的鎢礦物以顆粒形態(tài)存在，為查明金屬鎢顆粒在黑鎢礦石中的分布情況，隨機選取兩塊粒徑為50 mm 的黑鎢礦石，觀察礦石表面，發(fā)現(xiàn)石英與硅質(zhì)巖大面積黏結(jié)，但無明顯鎢顆粒存在。將其中一塊進(jìn)行隨機切割，另一塊在石英-硅質(zhì)巖界面處進(jìn)行破碎，通過統(tǒng)計各切割截面和石英-硅質(zhì)巖斷面中的鎢顆粒分布情況，對金屬鎢顆粒進(jìn)行區(qū)域分布描述。利用數(shù)碼顯微鏡對礦石切割截面進(jìn)行觀察，并對金屬鎢顆粒分布區(qū)域進(jìn)行不同倍率巖相采集，如圖2 所示。

圖2 為采集自第1 塊黑鎢礦石切割截面的可見鎢顆粒巖相，受所選黑鎢礦石內(nèi)石英下面硅質(zhì)巖的結(jié)構(gòu)、雜質(zhì)及調(diào)節(jié)焦距等因素的影響，顯微鏡下的石英礦物有明顯顏色差異。觀察巖相發(fā)現(xiàn)，鎢顆粒形狀各異，顆粒邊界曲線既有圓弧、直線等規(guī)則曲線，也有復(fù)雜、隨意的不規(guī)則曲線，且鎢顆粒尺寸差異也較大。分析鎢顆粒礦物周邊組分，對鎢顆粒的空間分布進(jìn)行描述，發(fā)現(xiàn)鎢顆粒主要分布在石英礦物中，石英-硅質(zhì)巖界面處有部分鎢顆粒，硅質(zhì)巖礦物中沒有鎢顆粒。

圖3(a)和圖3(b)分別為第2 塊黑鎢礦石在石英-硅質(zhì)巖界面處切割后得到的硅質(zhì)巖界面和石英界面形貌及局部區(qū)域1.5 倍放大圖?？梢钥闯?，兩塊礦物組分界面處均存在較多的鎢顆粒，且局部放大區(qū)域內(nèi)出現(xiàn)鎢顆粒聚集，但光滑的硅質(zhì)巖界面處未觀察到鎢顆粒分布。此外，石英礦物中能夠觀察到少量鎢顆粒嵌布，但硅質(zhì)巖礦物中未觀察到鎢顆粒。從硅質(zhì)巖和石英礦物的放大圖中可以看出，兩區(qū)域中的鎢顆粒均與石英礦物黏結(jié)在一起。

圖3 石英和硅質(zhì)巖的斷裂界面Fig. 3 Quartz and siliceous rock interface fracture map

1.3 鎢顆粒的特征尺寸分析

顆粒特征尺寸用于描述顆粒的幾何特性，主要包括長度、寬度、周長和面積4 個指標(biāo)，分別通過顆粒最小外接矩形長度、顆粒最小外接矩形寬度、顆粒邊界像素總和以及顆粒內(nèi)部像素總和4 個基本參數(shù)定量描述。長度和寬度通過數(shù)碼顯微鏡的配套測量軟件采用最佳外接矩形法進(jìn)行測量；引入Freeman 鏈碼（freeman chain code，F(xiàn)CC），通過MATLAB 計算鎢顆粒的邊界像素總和以及內(nèi)部像素總和，從而獲得顆粒的周長和面積。計算周長和面積時，需對圖像進(jìn)行中值濾波、掩模濾波以及二值分割預(yù)處理。隨機選取8 粒鎢顆粒，對其特征尺寸進(jìn)行計算，結(jié)果如表1 所示。

分析表1 中數(shù)據(jù)可知，鎢顆粒之間的尺寸差異較大，但鎢顆粒整體尺寸均較小，其中鎢顆粒長度范圍為0.922～2.966 mm，寬度范圍為0.137～2.233 mm，周長范圍為2.230～13.930 mm，面積范圍為0.070 3～3.612 5 mm2。

表1 鎢顆粒的特征尺寸Table 1 Characteristic sizes of different tungsten particles

2 原位加載實驗

2.1 實驗系統(tǒng)搭建

如圖4 所示，原位加載實驗通過無錫歐凱電子研制的萬能材料試驗機實現(xiàn)，主要由控制柜、液壓控制系統(tǒng)、壓力傳感器、加載壓頭等部分組成。通過該實驗裝置能夠?qū)崿F(xiàn)試件的動態(tài)加載和準(zhǔn)靜態(tài)加載。

圖4 原位加載實驗裝置Fig. 4 In-situ loading experimental device

處理后的片狀礦石最大截面尺寸為35 mm×28 mm，厚度為10 mm；試驗機的原配壓頭尺寸為20 cm×20 cm，壓頭尺寸過大，不符合原位加載試驗要求。為成功進(jìn)行加載實驗，考慮礦石各組分礦物中鎢顆粒尺寸以及萬能材料試驗機的結(jié)構(gòu)參數(shù)，自制直徑為2 mm 的原位加載實驗壓頭，如圖5所示。

圖5 自制原位加載壓頭示意圖（單位：mm）Fig. 5 Sketch map of self-made in-situ loading indenter (Unit: mm)

由于自制原位加載壓頭的直徑太小，如果加載速度過快，壓頭在與試樣接觸瞬間會產(chǎn)生極大的應(yīng)力，使礦石瞬間破壞，且壓頭接觸區(qū)域會發(fā)生“過粉碎”現(xiàn)象，導(dǎo)致試樣被加載處產(chǎn)生小粒徑碎石甚至粉末，從而無法統(tǒng)計試樣的破碎形狀和可靠的應(yīng)力、應(yīng)變數(shù)據(jù)。為避免礦石瞬間破壞，加載過程中的加載速率設(shè)置為0.005 和0.010 kN/s，直至試件被加載處出現(xiàn)宏觀破壞。原位加載實驗結(jié)果的可靠性完全取決于對礦石加載處發(fā)生宏觀破壞的瞬間判斷，由于自制原位加載壓頭直徑僅為2 mm，在實驗過程中難以用肉眼判斷加載處是否發(fā)生宏觀破壞，為得到加載處發(fā)生宏觀破壞的準(zhǔn)確時刻，將數(shù)碼顯微鏡安裝到加載裝置上，即可在計算機得到壓頭與加載部位的放大圖像，同時可實時觀察載荷施加過程中礦石加載處的結(jié)構(gòu)變化，從而能夠準(zhǔn)確判斷礦石加載處產(chǎn)生宏觀破碎時刻，該實驗方案原理如圖6 所示。

圖6 原位加載實驗示意圖Fig. 6 Schematic diagram of in-situ loading experiment

2.2 實驗結(jié)果分析

黑鎢礦石試件中石英礦物、硅質(zhì)巖礦物和鎢礦物3 種礦物組分均以聚集體形式存在，且各礦物組分相互黏結(jié)，不同礦物組分在黏結(jié)過程中形成了石英-硅質(zhì)巖界面和石英-鎢界面。在兩種礦物界面處進(jìn)行原位加載實驗，獲得礦物界面的最小破碎應(yīng)力，見表2。其中，實驗1～實驗4 的加載速率為0.010 kN/s，實驗5～實驗8 的加載速率為0.005 kN/s。分析表2 中數(shù)據(jù)可知，兩種加載速率下礦物組分界面的抗壓強度較接近，石英-硅質(zhì)巖界面的最小破碎應(yīng)力范圍為1.178 5～1.482 6 GPa，石英-鎢界面的最小破碎應(yīng)力范圍為1.335 5～1.542 0 GPa。

表2 礦物界面的最小破碎應(yīng)力Table 2 Minimum fracture stress of component mineral interface

各組分礦物的最小破碎應(yīng)力雖然反映了礦物破壞的臨界值，但在加載過程中，礦物內(nèi)部的力學(xué)性能變化卻無法通過最小破碎應(yīng)力表現(xiàn)。應(yīng)力-應(yīng)變曲線可以反映在外力作用下礦物發(fā)生的脆性、塑性、屈服、斷裂等各種形變過程。傳統(tǒng)應(yīng)力-應(yīng)變曲線分析方法是將曲線按照峰值應(yīng)力的百分比分為彈性階段、屈服階段、強化階段和局部變形階段，其中彈性階段與屈服階段、屈服階段與強化階段的臨界應(yīng)力分別為峰值應(yīng)力的25%和80%。

圖7 為0.010 和0.005 kN/s 兩種加載速率下石英-硅質(zhì)巖界面原位加載實驗得到的應(yīng)力-應(yīng)變曲線。對比分析發(fā)現(xiàn)，在較高的加載速率作用下，曲線更光滑。分析圖7(a)可知：在彈性階段和屈服階段，4 條曲線的變化趨勢比較相似，在該過程中，加載區(qū)域先產(chǎn)生彈性應(yīng)變，再產(chǎn)生塑性應(yīng)變，使加載區(qū)域內(nèi)部孔隙閉合；在局部變形階段，曲線長度較短，即在較短時間內(nèi)完成界面內(nèi)部裂紋擴展、聚集，并達(dá)到宏觀破壞的效果，說明黏結(jié)界面在較少的微觀裂紋聚集下即可發(fā)生宏觀破壞。分析圖7(b)可知：4 條曲線在彈性階段的形狀或凹或凸，即在載荷作用下，加載區(qū)域產(chǎn)生了較小的彈性應(yīng)變后立即產(chǎn)生塑性應(yīng)變，直至屈服階段，加載區(qū)域在較大應(yīng)力作用后，內(nèi)部微觀缺陷閉合；在強化階段，幾條曲線的應(yīng)力增長緩慢，且應(yīng)力時有下降，但應(yīng)變快速增大，此時內(nèi)部不斷產(chǎn)生微裂紋，并不斷聚集直至形成宏觀破壞。

圖7 不同加載速率下石英-硅質(zhì)巖組分界面的應(yīng)力-應(yīng)變曲線Fig. 7 Stress-strain curves of quartz-siliceous rock component interface under different loading rates

圖8 為0.010 和0.005 kN/s 兩種加載速率下石英-鎢界面原位加載實驗得到的應(yīng)力-應(yīng)變曲線。對比分析發(fā)現(xiàn)，在較高的加載速率作用下，曲線更光滑。分析圖8(a)可知，在彈性階段，曲線近似直線，在屈服階段，曲線呈“內(nèi)凹”形狀，表明彈性階段應(yīng)力可恢復(fù)，屈服階段降低了界面處的孔隙率。彈性階段和屈服階段的應(yīng)力、應(yīng)變區(qū)間不同，表明不同區(qū)域石英-鎢界面的結(jié)構(gòu)性能和力學(xué)性能也存在差異，這是由鎢顆粒的形狀、尺寸差異導(dǎo)致的。在局部變形階段，曲線存在明顯波動，且在相差較小的應(yīng)力范圍內(nèi)產(chǎn)生了較大應(yīng)變，表面石英-鎢界面的破壞需要較多微裂紋聚集，并在峰值應(yīng)力下產(chǎn)生宏觀破壞。分析圖8(b)可知：在彈性階段和屈服階段，實驗1 和實驗3 的應(yīng)力-應(yīng)變曲線趨勢相似，實驗2 和實驗4 的曲線趨勢相似，這是由于顆粒形狀、尺寸差異導(dǎo)致界面力學(xué)性質(zhì)存在差異；在強化階段，4 條曲線的變化趨勢相似，這是由于加載速率較低時，經(jīng)過彈性階段和屈服階段，較小應(yīng)力增長即可產(chǎn)生較大應(yīng)變，并形成足夠的微裂紋，最終微裂紋聚集形成宏觀破壞。

圖8 不同加載速率下石英-鎢組分界面的應(yīng)力-應(yīng)變曲線Fig. 8 Stress-strain curves of quartz-tungsten component interface under different loading rates

綜合分析兩種加載速率下黑鎢礦石礦物組分界面原位加載實驗得到的應(yīng)力-應(yīng)變曲線可知，加載速率對屈服階段和強化階段具有較大的影響。加載速率較大時，強化階段應(yīng)力會突然降低并發(fā)生波動；加載速率較小時，應(yīng)力更穩(wěn)定。

3 RVE 模型構(gòu)建及其破碎模擬

3.1 模型參數(shù)及邊界條件設(shè)定

非線性、多尺度的材料與結(jié)構(gòu)建模平臺DIGIMAT 可生成各種材料微觀結(jié)構(gòu)的實際RVE，用來描述復(fù)合材料的微觀結(jié)構(gòu)及其組分的材料特性。通過以石英礦物為基體相、硅質(zhì)巖礦物和鎢礦物為內(nèi)含物，構(gòu)建礦石RVE 模型，各礦物組分材料的物理屬性見表3。

表3 不同黑鎢礦石礦物組分的材料屬性Table 3 Mineral material properties of different wolframite components

根據(jù)真實組分礦物的形狀選取了代表性硅質(zhì)巖礦物和鎢礦物內(nèi)含物，計算其相關(guān)形態(tài)參數(shù)，使模型中硅質(zhì)巖礦物和鎢礦物的RVE 與石英礦物的RVE 黏結(jié)形成的界面更接近實際形狀，兩種內(nèi)含物的形態(tài)參數(shù)見表4。

表4 內(nèi)含物的形態(tài)參數(shù)Table 4 Inclusion morphological parameters

根據(jù)所定義的材料屬性和內(nèi)含物形態(tài)參數(shù)，由DIGIMAT-FE 模塊可生成礦石RVE 模型，如圖9所示。圖9 為非線性、多尺度的材料與結(jié)構(gòu)建模平臺DIGIMAT 建立的三維黑鎢礦石RVE 模型及其網(wǎng)格模型。圖9 中透明區(qū)域為石英礦物，藍(lán)色區(qū)域為鎢顆粒，紅色區(qū)域為硅質(zhì)巖礦物。

圖9 三維礦石RVE 模型和網(wǎng)格模型Fig. 9 3D ore RVE model and mesh model

ABAQUS 軟件是專業(yè)的有限元分析軟件，通過與建模軟件DIGIMAT 耦合可實現(xiàn)不同條件下礦石破碎模擬。模型構(gòu)建成功后，ABAQUS 軟件能夠讀取DIGIMAT 導(dǎo)出的iges 文件和inp 文件，并在此基礎(chǔ)上結(jié)合原位加載實驗設(shè)置相關(guān)條件，完成相同條件下的原位加載數(shù)值模擬，也可通過更改參數(shù)設(shè)置進(jìn)一步分析礦石的破碎過程。

參考原位加載實驗的夾具工作原理和載荷施加方式，在ABAQUS 中對RVE 模型進(jìn)行載荷加載類型選取和邊界條件設(shè)置。載荷種類選擇Surface traction（施加在面上的載荷），根據(jù)原位加載實驗的數(shù)據(jù)，設(shè)置對應(yīng)的載荷峰值，且載荷作用方向為豎直向下；邊界條件種類選取Mechanical-Encastre（力學(xué)-端部固定邊界條件），并選取RVE 模型4 個側(cè)面（與z軸平行端面）進(jìn)行完全固定約束。

3.2 組分界面破碎過程分析

分別在RVE 模型中的石英-硅質(zhì)巖界面和石英-鎢界面選取直徑為2 mm 的圓形區(qū)域作為載荷作用區(qū)域，參考黑鎢礦石原位加載實驗中礦石明顯破碎時的峰值應(yīng)力，依次對石英-硅質(zhì)巖界面和石英-鎢界面施加沿z軸反方向峰值為9.34 和10.31 kN 的面載荷，對RVE 模型進(jìn)行有限元求解。

表5 為在ABAQUS 軟件中應(yīng)用RVE 模型模擬石英-硅質(zhì)巖界面原位加載實驗得到的邊界、內(nèi)部和外部載荷作用區(qū)域結(jié)點的應(yīng)力和應(yīng)變，“*”標(biāo)記的結(jié)點位于石英礦物與硅質(zhì)巖礦物的黏結(jié)界面上。

將表5 中的結(jié)點編號重新按照自然數(shù)編號繪制應(yīng)力、應(yīng)變曲線，如圖10 所示。結(jié)合應(yīng)力分析圖10(a)可知，載荷作用區(qū)域內(nèi)部和邊界的應(yīng)力較大，而載荷作用區(qū)域外部的應(yīng)力較小，載荷作用區(qū)域內(nèi)部、邊界、外部的平均應(yīng)力分別為1.487 18、1.457 83 和0.344 65 GPa。載荷作用區(qū)域邊界和載荷作用區(qū)域內(nèi)部不僅平均應(yīng)力接近，而且兩個區(qū)域的應(yīng)力波動均比較大。此外，因黏結(jié)界面結(jié)點在各區(qū)域均有分布，使得界面不同結(jié)點之間的應(yīng)力差距較大，應(yīng)力范圍為0.230 54～1.621 67 GPa。結(jié)合應(yīng)變分析圖10(b)可知，載荷作用區(qū)域邊界的應(yīng)變遠(yuǎn)大于載荷作用區(qū)域內(nèi)部和外部，且載荷作用區(qū)域內(nèi)部和外部存在部分結(jié)點應(yīng)變重疊，載荷作用區(qū)域邊界、內(nèi)部、外部的平均應(yīng)變分別為1.40×10?4、0.03×10?4和0.60×10?4。石英-硅質(zhì)巖界面不同結(jié)點之間的應(yīng)變差距較大，應(yīng)變范圍為0.38×10?4～1.35×10?4。

表5 石英-硅質(zhì)巖界面結(jié)點的應(yīng)力和應(yīng)變Table 5 Quartz-siliceous rock interface junction stress and strain

圖10 石英-硅質(zhì)巖界面原位加載模擬結(jié)點的應(yīng)力、應(yīng)變分布Fig. 10 Stress and strain distributions of in-situ loading simulation node of the quartz-siliceous rock interface

表6 為ABAQUS 利用RVE 模型模擬石英-鎢界面原位加載實驗得到的邊界、內(nèi)部和外部載荷作用區(qū)域結(jié)點的應(yīng)力和應(yīng)變，其中“*”標(biāo)記結(jié)點位于石英礦物與鎢礦物的黏結(jié)界面上。

將表6 中的結(jié)點編號重新按照自然數(shù)編號繪制應(yīng)力、應(yīng)變曲線，如圖11 所示。分析圖11(a)可知，載荷作用區(qū)域內(nèi)部和邊界的應(yīng)力較大，而載荷作用區(qū)域外部應(yīng)力較小，載荷作用區(qū)域內(nèi)部、邊界、外部的平均應(yīng)力分別為1.675 82、1.690 05、0.561 44 GPa。此外，因界面結(jié)點在各區(qū)域均有分布，使得不同界面結(jié)點之間的應(yīng)力差距較大，其應(yīng)力范圍為0.513 55～1.614 53 GPa。分析圖11(b)可知，載荷作用區(qū)域邊界和內(nèi)部的應(yīng)變略大于載荷作用區(qū)域外部，但載荷作用區(qū)域邊界和內(nèi)部的最小應(yīng)變小于載荷作用區(qū)域外部，且3 個區(qū)域的應(yīng)變波動范圍均比較大，其中載荷作用區(qū)域邊界、內(nèi)部、外部的平均應(yīng)變分別為1.01×10?4、1.44×10?4和0.84×10?4。石英-鎢界面不同結(jié)點之間的應(yīng)變差距較大，應(yīng)變范圍為0.57×10?4～1.90×10?4。

表6 石英-鎢界面結(jié)點的應(yīng)力和應(yīng)變Table 6 Quartz-tungsten interface junction stress and strain

圖11 石英-鎢界面原位加載模擬結(jié)點的應(yīng)力、應(yīng)變分布Fig. 11 Stress and strain distributions of in-situ loading simulation node of the quartz-tungsten interface

綜合分析石英-硅質(zhì)巖界面和石英-鎢界面在原位載荷下的破碎過程，發(fā)現(xiàn)載荷作用區(qū)域邊界和礦物界面重合結(jié)點的應(yīng)力、應(yīng)變均比較大，即宏觀破壞發(fā)生在載荷作用區(qū)域邊界，且礦物界面處破壞特征更明顯。但石英-鎢界面3 個載荷作用區(qū)域的應(yīng)變波動范圍均較大，部分結(jié)點的應(yīng)變重疊，平均應(yīng)變差距小，即石英-鎢界面在原位載荷下不僅在載荷作用區(qū)域邊界和內(nèi)部發(fā)生明顯破壞，在載荷作用區(qū)域外部也存在破壞。

4 結(jié) 論

通過對鎢礦石進(jìn)行巖相分析，獲取了有用礦物尺寸特征及其空間分布，并在此基礎(chǔ)上構(gòu)建符合組分礦物微觀結(jié)構(gòu)屬性的RVE 模型，利用礦石RVE 模型模擬礦物界面原位加載，結(jié)合原位加載實驗，研究了礦石不同組分礦物黏結(jié)界面的破碎力學(xué)特征。

(1) 所研究的黑鎢礦石中，有用礦物以顆粒狀分布在礦石內(nèi)部，主要分布在石英礦物內(nèi)及其與硅質(zhì)巖礦物黏結(jié)界面處，純石英礦物區(qū)域存在少量鎢顆粒分布，硅質(zhì)巖礦物區(qū)域無鎢顆粒分布，且鎢顆粒的面積、周長等尺寸均較小，顆粒之間尺寸差異較大。

(2)對鎢礦石試樣進(jìn)行加載速率分別為0.005、0.010 kN/s 的原位加載實驗，分析實驗結(jié)果可知，石英-硅質(zhì)巖界面的最小破碎應(yīng)力范圍為1.178 5～1.482 6 GPa，石英-鎢界面的最小破碎應(yīng)力范圍為1.335 5～1.542 0 GPa。加載速率對礦物界面形變的影響較大，且加載速率較大時，強化階段時應(yīng)力易突然降低并不斷波動。

(3) 通過對礦石RVE 模型原位破碎模擬結(jié)果分析，發(fā)現(xiàn)破碎主要發(fā)生在載荷作用區(qū)域邊界，且石英-鎢界面產(chǎn)生的應(yīng)力、應(yīng)變均大于石英-硅質(zhì)巖界面，但兩組界面中石英礦物的破碎力學(xué)特征參數(shù)大于鎢礦物和硅質(zhì)巖礦物，即界面組成礦物中石英礦物優(yōu)先形成破壞。