關(guān)于顆粒流軟件PFC的離散元數(shù)值模擬參數(shù)標(biāo)定

1、關(guān)于顆粒流軟件PFC的離散元數(shù)值模擬參數(shù)標(biāo)定工程背景隨著工程規(guī)模的日益擴(kuò)大，作為建筑材料以及地基介質(zhì)的巖土體一直以來(lái)都是人們研究的熱點(diǎn)。目前，解決巖土類工程問(wèn)題的主要研究手段有三大類，分別是理論分析、室內(nèi)實(shí)驗(yàn)和數(shù)值模擬。其中，數(shù)值模擬又因其便捷、高效以及可控性成為當(dāng)今最為常用的研究方法。常用的數(shù)值模擬方法主要基于有限元分析理論或離散元分析理論，運(yùn)用計(jì)算機(jī)軟件對(duì)實(shí)際的工程問(wèn)題進(jìn)行模擬，從而發(fā)現(xiàn)問(wèn)題、得到內(nèi)在規(guī)律并作出研究結(jié)論，對(duì)于指導(dǎo)工程實(shí)踐具有重大意義。研究方法針對(duì)不連續(xù)介質(zhì)問(wèn)題以及大變形問(wèn)題，近年來(lái)新興的離散元理論逐步代替?zhèn)鹘y(tǒng)的有限元理論被人們采用。其主要理論依據(jù)是牛頓第二定律和力與位移的關(guān)

2、系，求解思路為：將求解空間離散為分離的單元陣，相鄰兩個(gè)單元之間用合理的元件相連接；基本變量是單元間的相對(duì)位移，通過(guò)位移和力的關(guān)系可以分別得到單元之間的法向、切向作用力；求單元上受到的合力及合力矩，由牛頓運(yùn)動(dòng)定律計(jì)算出單元體的加速度；最后，積分得到單元的速度和位移。通過(guò)以上計(jì)算過(guò)程可以得到任意時(shí)刻所有單元體的位移、速度和加速度。在基于離散元分析理論開發(fā)出來(lái)的計(jì)算機(jī)軟件中，最為廣泛使用的是美國(guó)ITASCA公司開發(fā)的UDEC(universaldistinctelementcode)、3DEC(3-dimensionaldistinctelementcode)塊體離散元程序，和PFC2Dpartic

3、leflowcodein2dimensions)、PFC3D(particleflowcodein3dimensions）軟件。前者用于模擬巖石的力學(xué)過(guò)程，后者則可以模擬顆粒流問(wèn)題，如粗粒土、砂土的力學(xué)現(xiàn)象等?，F(xiàn)階段，國(guó)內(nèi)使用PFC軟件進(jìn)行離散元數(shù)值模擬研究的情況并不多見，目前國(guó)內(nèi)已出版的相關(guān)中文書籍也僅有離散元顆粒流軟件（PFC）在道路工程中的應(yīng)用和軟件自帶的全英文版幫助手冊(cè)。PFC軟件從微觀結(jié)構(gòu)角度著手，將土體看成是由土顆粒組成，其宏觀力學(xué)特性取決于顆粒與顆粒之間粘結(jié)（接觸）方式的幾何、物理特性。在PFC中用等厚度的剛性圓盤代表顆粒，顆粒之間允許重疊，同樣遵循牛頓第二定律。顆粒的破壞主要有

4、剪切破壞和張拉破壞兩種方式，當(dāng)顆粒發(fā)生破碎時(shí)土體的宏觀力學(xué)特性也會(huì)發(fā)生相應(yīng)變化，即介質(zhì)內(nèi)顆粒接觸狀態(tài)的變化決定了介質(zhì)的本構(gòu)關(guān)系。因此，在PFC計(jì)算中無(wú)需給定材料的宏觀本構(gòu)關(guān)系和對(duì)應(yīng)的參數(shù)，這些傳統(tǒng)的參數(shù)和力學(xué)特性在程序中可以自動(dòng)得到，需要使用者自行定義的是顆粒的幾何力學(xué)參數(shù)，包括顆粒級(jí)配、剛度、摩擦力、粘結(jié)介質(zhì)強(qiáng)度等細(xì)觀力學(xué)參數(shù)。參數(shù)標(biāo)定方法及常見問(wèn)題需要注意的一點(diǎn)是，運(yùn)用PFC軟件進(jìn)離散元數(shù)值模擬時(shí)，可以自行定義顆粒的物理力學(xué)參數(shù)，也可以根據(jù)實(shí)際情況標(biāo)定出較為真實(shí)的顆粒參數(shù)。而通常情況下，為了使模擬得到的結(jié)果更具有說(shuō)服力，研究者通常會(huì)采用后一種方法獲取顆?；緟?shù)。因此，數(shù)值模擬的首要任務(wù)是

5、進(jìn)行參數(shù)標(biāo)定，而在標(biāo)定過(guò)程中也會(huì)遇到一些共性的問(wèn)題，下面以筆者的實(shí)際模擬過(guò)程為例進(jìn)行簡(jiǎn)要說(shuō)明。筆者的主要研究?jī)?nèi)容是運(yùn)用PFC3.0軟件進(jìn)行砂土的顆粒破碎模擬，選用真實(shí)的室內(nèi)三軸剪切試驗(yàn)，在PFC中對(duì)其進(jìn)行三軸試驗(yàn)的數(shù)值模擬，將模擬得到的宏觀應(yīng)力應(yīng)變曲線與真實(shí)的室內(nèi)實(shí)驗(yàn)得到的應(yīng)力應(yīng)變曲線進(jìn)行比對(duì)，如果二者在線性、斜率和峰值強(qiáng)度上都能夠很好的定量吻合，則說(shuō)明數(shù)值模擬中采用的顆粒參數(shù)是相對(duì)真實(shí)可靠的。在PFC中如何根據(jù)宏觀參數(shù)確定細(xì)觀參數(shù)是一個(gè)困難的問(wèn)題。在標(biāo)定階段，影響宏觀應(yīng)力應(yīng)變曲線的因素主要有顆粒之間的摩擦系數(shù)、顆粒與墻體間的摩擦系數(shù)、顆粒的切向及法向剛度、試樣的初始孔隙率和剪切速率。如何選取

6、適當(dāng)?shù)膮?shù)成為數(shù)值模擬的一個(gè)關(guān)鍵性問(wèn)題，參數(shù)標(biāo)定的具體方案是：先擬定好以上參數(shù)，再通過(guò)控制變量法逐一調(diào)整每個(gè)因素，直至最后的曲線與實(shí)際曲線定量地對(duì)應(yīng)上。值得注意的一點(diǎn)是，先要保證參數(shù)的數(shù)量級(jí)是正確的，在此基礎(chǔ)上再微調(diào)參數(shù)的具體數(shù)值。研究結(jié)論通過(guò)筆者近一個(gè)月的標(biāo)定過(guò)程，最后總結(jié)出一些關(guān)于三軸試驗(yàn)中參數(shù)標(biāo)定的規(guī)律：1)在三軸試驗(yàn)條件下，容易發(fā)生顆粒破碎的砂土的摩擦系數(shù)在0.70.9,二維模擬情況下選取的摩擦系數(shù)盡量控制在0.9以下，否則模擬將失真。(2)切向剛度、法向剛度的數(shù)量級(jí)在e9e10,也可以調(diào)整切向、法向的剛度比；(3)二維情況下的初始孔隙率在0.20.35之間，曲線的線性則主要取決于初始孔隙率；4)應(yīng)力應(yīng)變曲線的峰值強(qiáng)度很大程度上由加載速率