總結分析了在離散元方法中研究顆粒破碎的兩種方法，黏結模型和顆粒破碎模型。
關鍵詞：顆粒破碎離散元數值模擬
中圖分類號：P61 文獻標識碼：A 文章編號：1672-3791（2015）01（c）-0000-00
Review of particle breakage of rock and soil under
DEM analysistechniques
Zheng Lin， Zhao Zhen-ying
（School of Ocean Science and Engineering， Shanghai Maritime University， Shanghai 201306， China）
Abstract： This paper outlines the significance of particle breakage of rock and soil in theory and engineering. Combined with the presentexperimental study on particle breakageof physics， synthetically expounds the advantages of the numerical simulation technology， and the method of DEM can be more effective to study the particle breakage of rock and soil. Two methods of particle breakage in DEM was summarized and analyzed， the bond modeland the particle failure criterion.
Key words：Particle breakage； DEM； Numerical simulation
巖土材料（包括粘土和砂土）是粒狀土顆粒的集合體。作為一種摩擦材料，它具有壓硬性和剪脹性（含剪脹與剪縮）。除巖土顆粒本身的力學性之外，巖土顆粒間相互作用的力學特性、顆粒間相互接觸的方式等也是影響巖土宏觀力學特性的重要因素。
在很小主應力的條件下，巖土顆粒不容易發生破碎，只有巖土顆粒間的互相作用才能決定土的變形以及強度。而有較大平均主應力的時候，就很容易導致砂土顆粒的破碎，在這種條件下，除了顆粒之間的相互作用是影響砂土變形和強度的主要原因外，還有顆粒破碎對砂土力學特性的影響。相比砂土而言，粘土有很小的顆粒粒徑，因此破碎現象不容易發生，只是在壓力作用極其高的情況下，才會對粘土顆粒發生破碎影響了變形和強度進行考慮。然而在工程實際中會較多地出現砂土顆粒破碎的現象，現在已經有很多學者對巖土顆粒的破碎現象進行了相應的研究。
伴隨著規模不斷擴大的堤壩、隧道等，高壓應力下的巖土工程課題也越來越多，巖土顆粒材料在高應力下的物理力學特性研究也變得越來越重要。顆粒材料在高應力水平下的一種基本現象是顆粒破碎現象，對顆粒材料的力學特性受破碎情況的影響進行分析，不僅有實際工程價值，還具有理論研究意義。
1物理試驗研究現狀
Terzaghi[1]對砂樣進行壓力高達96.5MPa一維壓縮試驗時，發現有十分明顯的顆粒破碎現象。另外，他發現孔隙比與有效應力的對數坐標的斜率在高應力水平下，在應力的增長下增加是很微小的。DeSouza[2]對三種不同的砂進行了壓力高達138MPa的試驗。結果顯示，在壓縮的過程當中，壓縮指數發生變化的地方出現了一個破碎點。分析顆粒結果表明，在應力高于破碎點的時候，具有十分顯著的破碎現象；有較高的相對密度，則對應的破碎點的應力也就越高；增加顆粒棱角度會降低破碎點的應力。Harremoes研究了顆粒礦物成分對其破碎特性的影響，發現在高達138MPa應力的作用下，改變礦物和晶體尺寸不會對破碎點的應力造成影響。不過破碎點以上，具有十分明顯的大量裂隙顆粒的破碎現象。Roberts對相同材料進行了試驗，壓力達到276MPa，與Harremoes和De Souza的結果非常相似。與此同時，DeBeer[3]為了驗證Terzaghi提出的說法，即可以忽略不計在9.8MPa的應力以下產生的破碎；在9.8 MPa以上，隨著應力的增加，極大地降低了破碎。他對均勻砂樣進行了一維壓縮試驗，發現應力為15MPa時，會產生十分明顯的破碎現象，但是在34MPa以上，伴隨著應力的增加會漸次降低破碎。
綜觀目前對顆粒破碎的研究進展，很容易發現人們的注意點主要是在對產生破碎的影響因素的認識上，對破碎本身機理系統性的認識則比較缺乏。顆粒破碎對巖土體力學性質的影響等方面做得還很不夠，尤其是沒有深入研究微觀破碎機理。當前的研究成果還遠遠不能說明巖土體顆粒破碎的力學機理，與實際工程的應用還有很遠的距離。
2數值模擬方法
隨著數值模擬方法以及高性能計算機的發展，巖土顆粒破碎力學特性研究對數值模擬方法的應用越來越廣泛。由于室內試驗難以從細觀結構角度來分析認識顆粒破碎現象的機理和帶來的力學特性的變化，數值模擬的方法可以有效地解決這些問題，為研究巖土體顆粒破碎的力學特性提供了途徑。目前，用于巖土力學的數值模擬方法主要有FEM（有限元）和DEM（離散元）。
2.1有限元方法
解決數學物理和工程問題的數值方法叫做有限元方法，也稱為有限單元法，是矩陣方法在彈性力學和結構力學等領域中的應用和發展。因其有效性和通用性，有限元方法被廣泛地應用在了工程分析中，已成為計算機輔助制造和計算機輔助設計的重要組成部分。
將求解區域離散為一組有限個，且按一定方式互相聯接在一起的單元組合體，是有限元方法的基本思想。由于能按不同的聯接方式組合各單元，且可以具有形狀不同的單元，所以可以幾何形狀復雜的求解域進行模擬。利用在每一個單元內假設近似函數來分片地表達求解域上的未知場函數，是有限元方法作為數值分析方法的一個重要特點。
閆澍旺在特制的壓縮儀中對其壓縮特性進行了試驗，不同粒徑組試樣的相應的變形模量和屈服應力，顆粒級配情況，破碎規律進行了研究，還通過彈塑性有限元分析對石料破碎的規律性進行了分析驗證。米占寬[4]基于Ueng和Chen剪脹方程，通過分析三軸剪切試驗過程中的能量平衡，提出了考慮顆粒破碎的剪脹方程及其參數確定方法。通過對比分析三軸CD剪切試驗成果，和分析工程實例的有限元數值，表明所提模型對材料的剪脹特性有較好的反映。賈宇峰[5]在土體的能量平衡方程中引入顆粒破碎耗能，進而根據能量方程建立了考慮顆粒破碎的粗粒土本構模型。同時結合三軸試驗，根據剪切過程中的能量平衡關系，建立了考慮顆粒破碎影響的摩擦系數公式，對模型參數進行確立，驗證了模型是有效的。
2.2離散元方法
基于分子動力學原理，美國學者Cundall [6-7]提出了離散單元法，是對巖石力學問題進行分析的一種間斷性數值模擬方法。De Josselin de Jong和Verrujit[8]所做的光彈試驗驗證了此種方法的可靠性。因離散單元法在離散物質分析方面具有一定的優勢，已成為近來年解決非連續介質問題的一種有效的數值模擬方法。PFC2D是以離散元方法為基礎，基本單元為剛性圓盤或圓球，通過基本單元的相互作用和運動對實際問題進行模擬的離散元軟件[9]。離散單元法首次提出是用來分析巖石力學，之后發展迅速，應用領域越來越廣泛。Cundall和Strack [10]又將其應用到土力學中，并開發了二維圓盤程序BALL和三維圓球程序TRUBAL，這就是顆粒流程序PFC的雛形。
3 顆粒破碎在離散元數值模擬中的方法
離散元方法在巖土工程領域得到了廣泛應用，現在已經成為巖土工程領域數值模擬的有效方法。目前，基于離散元方法來模擬顆粒破碎的途徑有黏結模型和顆粒破碎模型。
3.1黏結模型
將基于若干離散顆粒通過接觸點粘結形成的顆粒簇視為一個可破碎顆粒，在外力作用下該顆粒簇內某個接觸點處的相互作用力滿足一定強度條件時，則解除該接觸點的粘結（可認為發生一個斷裂事件）。當某部分顆?；蚰硞€顆粒與其余部分顆粒間的接觸粘結都解除后，則從原來的顆粒簇破碎出該部分顆?；蛟擃w粒[11]。
圖1給出了黏結模型發生破碎的示意圖。在加荷開始之前，簇顆粒沒有破碎，顆粒簇的實體面積應包括顆粒簇內部孔隙的面積，見圖1（a）；開始加荷后，顆粒簇會在荷載增加到一定程度時產生破碎，這時，就會釋放出原來歸屬于顆粒簇實體面積的那部分內部孔隙，形成試樣的外部空隙，見圖1（b）所示。在內部空隙不斷釋放的同時，會不斷減小試樣中的實體總面積。
圖1黏結模型破碎示意圖
李凡[12]基于PFC理論，利用顆粒黏結模型，通過記錄試樣黏結破損的數目和空間位置來反映結構性巖土材料破損特性。劉君[13]利用黏結模型對粗粒料進行數值模擬，發現顆粒破碎的細觀演化規律和破碎帶的分布范圍。BuddhimaIndraratna[14]運用黏結模型形成二維尖角顆粒模型，利用細觀力學參數如接觸力和粘結力的分布在循環荷載作用下的發展來解釋顆粒破碎的力學機制。
3.2顆粒破碎模型
基于單個離散顆粒，當其所有接觸點處的接觸力滿足一定條件時，認為該顆粒將發生顆粒破碎模型的破碎[15]。這是由Lobo Guerrero[16]提出的顆粒破碎模型的方法，當顆粒達到破碎條件以后，原顆粒用8個小顆粒進行替代，由此來模擬顆粒破碎的現象。顆粒破碎模型的方法很好的解決了顆粒數量以及發生多級破碎的問題。
圖2給出了顆粒破碎模型發生破碎的示意圖。由于顆粒破碎是沿著顆粒中的孔隙發生破壞，破碎后顆粒中的孔隙會釋放出來，變成試樣的外部孔隙，陰影部分表示顆粒中的孔隙。
圖2顆粒破碎模型破碎示意圖
Lobo Guerrero [16]考慮了配位數對顆粒破碎的影響并基于假定的顆粒破碎模型考察了顆粒破碎演及其對顆粒材料力學行為的影響。Subero與Ghadiri[17]考察了沖擊速度對顆粒破碎模型的影響，并總結了規律。
4 結語
離散元發展到如今相對于有限元有著它的巨大優勢。有限元在研究砂土的力學特征時，采用的是連續介質理論，造成剪切強度和變形的細觀力學機理無法從顆粒層面進行揭示，而對顆粒形狀、孔隙率等試樣要素變化對模擬結果的影響也無法反映出來。不過離散元剛好彌補了這一問題，可以從細觀參數的研究來分析顆粒材料的宏細觀力學特性，在分析巖土體大變形破壞問題及其宏細觀機理研究方面具有獨特的優勢。
在模擬破碎時，黏結模型不用對破碎模式進行考慮，可由計算自動獲得破碎后產物的分布和尺寸，但其與顆粒簇內單個顆粒尺寸分布有關系，并和顆粒間的強度分布關系緊密，再加上顆粒簇內的每個離散顆粒均參與將顆粒簇內的每個離散顆粒都參與接觸力的更新、位置的更新及接觸判斷，故使計算量龐大。在計算量上顆粒破碎模型的方法則具有一定的優勢，不過需要對顆粒破碎模式和破碎準則進行定義，而顆粒的應力分布與微細結構決定了顆粒是否破碎及按照何種模式破碎。但顆粒破碎模型把破碎前的顆粒視為完整的，沒有考慮微細結構，且經典離散單元把顆粒視作剛體，雖對基于單個顆粒的名義應力進行了定義，但沒給出顆粒內名義應力的分布。從變形體角度考慮，顆粒內部應力分布與顆粒表面接觸力的大小及分布密切相關，由此可知顆粒破碎及破碎模式與顆粒配位數有關。
參考文獻：
[1]TERZAGHI K， PECK R B. Soil mechanics inengineering practice[M]. NewYork： John Wiley and Sons Inc， 1948： 65-67.
[2]DE SOUZA J M. Compressibility of sand at high pressure[D].Massachusetts Institute ofTechnology， 1958.
[3]DeBeerEE.Thescaleeffectinthetranspositionofthe resultsofdeepsoundingtestsontheultimatebearing capacityofpilesandcaissonfoundations [J]. Geotechnique， 1963， 13（1）： 39-75.
[4]米占寬， 李國英， 陳生水. 基于破碎能耗的粗顆粒料本構模型[J]. 巖土工程學報， 2012， 34（10）： 1801-1811. （MI Zhan-kuan， LI Guo-ying， CHEN Sheng-shui. Constitutive model for coarse granular materials based on breakage energy[J].ChineseJournalofGeotechnicalEngineering，2012， 34（10）：1801–1811. （inChinese））
[5]賈宇峰. 考慮顆粒破碎的粗粒土本構關系研究[D]. 大連： 大連理工大學， 2008. （JIA Yu-feng. Coarse granular soil constitutive model incorporating particle breakage[D]. Dalian： Dalian University of Technology， 2008. （in Chinese））
[6]Cundall PA. A computer model for simulating progressive large scale movements in blocky system.Proceedings of Symposium of the International Society of Rock Mechanics. Rotterdam： A.A. Balkema， 1971（1）： 8-12.
[7]Cundall PA. The measurement and analysis of acceleration in rock slopes. Ph.D. Dissertation， University of London， Imperial college of Science and Technology， 1971.
[8]Cundall PA， Strack O D L. A discrete numerical method for granular assemblies[J]. Geotechnique， 1979， 29（1）： 47-65.
[9]ItascaConsulting Group， Inc. PFC2D theory and background. Minneapolis： Minnesota， 2004.
[10]楊全文. 離散元法干顆粒接觸模型研究及微機可視化程序設計[D]. 中國農業大學， 2001.
[11]ChengYP，BoltonMD，NakataY.CrushingandplasticdeformationofsoilssimulatedusingDEM[J]. Geotechnique， 2004， 54： 131-141.
[12]李凡. 巖土材料破損特性的顆粒流研究[J]. 土木工程學報， 2007， 40（9）： 78-81. （Li Fan. A study on the breakage properties of geological materials using particle flow simulation[J]， China civil engineering journal， 2007， 40（9）： 78-81. （in Chinese））
[13]劉君， 劉福海， 孔憲京. 考慮破碎的堆石料顆粒流數值模擬[J]. 巖土力學， 2008， 29（增）： 107-112. （Liu Jun，LiuFuhai，KongXianjing.Particleflowcodenumericalsimulationofparticlebreakageofrockfill[J]， Rock and Soil Mechanics， 2008， 29（增）： 107-112. （in Chinese））
[14]J. Lackenby， B. Indraratna， G. McDowell， D. Christie. Effect of confining pressure on ballast degradation and deformation under cyclic triaxial loading[J]. Géotechnique， 2007， 57（6）： 527-536.
[15]McdowellGR，BoltonMD，RobertsonD.Thefractalcrushingofgranularmaterials[J].Journalofthe Mechanics and Solids， 1996， 44： 2079-2102.
[16]Lobo G S. Evaluation of crusing in granular materials using the discrete element method and fractal theory[D]. University of Pittsburgh， 2006.
[17]Subero J， Ghadiri M. Breakage patterns of agglomerates[J]. Powder Technology， 2001， 120： 232-243.