章春煒，路軍富，鐘英哲

(成都理工大學(xué), 四川成都 610059)

PFC中應(yīng)力測量圓半徑淺析

章春煒，路軍富，鐘英哲

(成都理工大學(xué), 四川成都 610059)

在PFC軟件中，可用應(yīng)力測量圓來提取模型中的應(yīng)力值，但是不同半徑的應(yīng)力測量圓所提取的應(yīng)力則不同。文章通過數(shù)值模擬與實際計算的方法相結(jié)合，研究不同半徑的應(yīng)力測量圓所提取的應(yīng)力值與實際應(yīng)力值的擬合程度。研究結(jié)果表明：在顆粒半徑最大為30 mm、最小半徑為5 mm情況下，半徑為0.3 m的應(yīng)力測量圓所提取的應(yīng)力值與實際應(yīng)力值擬合度最高。

PFC； 應(yīng)力測量圓； 半徑； 模型； 應(yīng)力

隨著顆粒離散元方法應(yīng)用越來越廣泛，PFC作為一種典型的離散元軟件，在巖石類材料基本特性、顆粒物質(zhì)動力響應(yīng)等問題中具有較多的應(yīng)用。朱煥春[1]結(jié)合闡述了PFC應(yīng)用現(xiàn)狀，并結(jié)合工程實例說明PFC可以針對性的解決一些巖土體問題。周國慶[2]運用PFC 2D研究顆粒流程序中的阻尼參數(shù)的適用性，得出了在PFC 2D中局部阻尼和黏性阻尼的定義及適用范圍。吳順川[3]通過對比分析室內(nèi)卸載巖爆試驗與 數(shù)值模擬研究，得出PGC 3D數(shù)值模擬可以部分代替室內(nèi)試驗。王延可[4]得出用PFC 數(shù)值模擬試驗研究巖爆特性的效果非常好。國內(nèi)學(xué)者對用PFC的實際應(yīng)用得出了較多的研究成果，但是對于PFC軟件內(nèi)置功能的研究較少。

本文通過數(shù)值模擬方法建立砂卵石地層模型，并在重力作用下達(dá)到平衡，研究不同半徑應(yīng)力測量圓所提取的應(yīng)力與實際計算應(yīng)力的擬合程度，得出下最合適的應(yīng)力測量圓的半徑。

1 應(yīng)力測量圓理論

應(yīng)力是一個連續(xù)的量，因此不存在于每個粒子上，因為介質(zhì)是離散的。在離散點的PFC 2D模型中，接觸力和粒子位移都能被計算。這些矢量在研究微觀行為時十分有用，但不能直接轉(zhuǎn)移到一個連續(xù)模型。因此從微觀模型到連續(xù)模型需要通過公式轉(zhuǎn)化，轉(zhuǎn)化公式如：

(1)

由公式(1)可知，平均應(yīng)力主要與測量圓的體積有關(guān)，而測量圓大小則與所取半徑有關(guān)，所以確定所需應(yīng)力測量圓的半徑是重要的一步。通過不同應(yīng)力測量圓的半徑所取得的應(yīng)力與計算所得應(yīng)力相比較，最終確定所需應(yīng)力測量圓的半徑。

2 應(yīng)力測量圓半徑確定

用半徑為0.2 m、0.3 m、0.4 m、0.5 m的應(yīng)力圓獲取中密天然含水率狀態(tài)下PFC 2D離散元模型自重應(yīng)力平衡后的應(yīng)力值，該模型寬50 m，高45 m，土顆粒密度ρ=2.0 g/cm3，顆粒半徑最小為5 mm，最大顆粒半徑為30 mm。

由于該離散元模型中，沒有構(gòu)造應(yīng)力，只有自重應(yīng)力。對于土體自重應(yīng)力計算可通過公式：σ=γz=ρgz得出，并將計算得出模型每一層的自重應(yīng)力與每一層應(yīng)力測量圓測得的自重應(yīng)力進(jìn)行比較，所得對比曲線如圖1所示。

分析圖1可得，在半徑R=0.2 m時應(yīng)力測量圓所得自重應(yīng)力曲線波動較大，并且最大值遠(yuǎn)遠(yuǎn)高于計算值。因此在

(a) R=0.2 m

(b) R=0.3 m

(c) R=0.4 m

(d) R=0.5 m圖2 計算自重應(yīng)力與測量圓自重應(yīng)力對比

該顆粒級配下，用半徑為0.2 m的應(yīng)力測量圓不能準(zhǔn)確測得顆粒之間的應(yīng)力，可以直接排除。對于應(yīng)力測量圓半徑為0.3 m、0.4 m、0.5 m時，其自重應(yīng)力曲線圖波動都比較小。對比應(yīng)力測量圓半徑為0.3 m、0.4 m時，應(yīng)力測量圓半徑為0.5 m時所測得的自重應(yīng)力曲線與計算所得自重應(yīng)力曲線擬合度較差。為了更好地確定各半徑應(yīng)力測量圓所得自重應(yīng)力曲線與其對應(yīng)的計算所得自重應(yīng)力曲線的擬合程度，計算了各半徑應(yīng)力測量圓所得的自重應(yīng)力曲線與目標(biāo)直線的擬合度，擬合度取值見表1。

表1 自重應(yīng)力曲線擬合度取值

分析表1可得，在應(yīng)力測量圓半徑為0.3 m時，其擬合度R2=0.94，遠(yuǎn)遠(yuǎn)大于應(yīng)力測量圓半徑為0.4 m和測量圓半徑為0.5 m時的擬合度，因此在該顆粒級配下，應(yīng)力測量圓半徑應(yīng)選擇R=0.3 m。

3 結(jié)論

(1) 在顆粒半徑最大為30 mm、最小半徑為5 mm的級配下，應(yīng)力測量圓半徑最合適半徑為0.3 m。

(2) 選用不同半徑的應(yīng)力測量圓所提取的應(yīng)力值有較大的差別，因此只有選取合適半徑的應(yīng)力測量圓時，所提取的應(yīng)力才能與實際應(yīng)力比較符合。

(3) 用PFC 2D能較好的模擬出砂卵石地層在自重應(yīng)力平衡后的實際應(yīng)力情況。

[1] 朱煥春．PFC及其在礦山崩落開采研究中的應(yīng)用 [J]. 巖石力學(xué)與工程學(xué)報, 2006,25(9): 1927-1931.

[2] 周國慶, 周杰, 陸勇,等．顆粒流程序(PFC 2D )中阻尼參數(shù)的適用性研究[J]. 中國礦業(yè)大學(xué)學(xué)報，2011，40(5): 667-672.

[3] 吳順川，周喻，高斌．卸載巖爆試驗及 PFC 3D 數(shù)值模擬研究[J]. 中巖石力學(xué)與工程學(xué)報，2010，29(Z2):4082-4088.

[4] 王延可, 李天斌, 陳國慶, 等. 巖爆特性PFC 數(shù)值模擬試驗研究[J]. 現(xiàn)代隧道技術(shù)，2013, 50(4): 98-103.

章春煒(1991～)，男，在讀碩士，研究方向為地質(zhì)工程。

TU478

A

[定稿日期]2017-07-24