何詠睿 ，朱 晟 ，2，武利強(qiáng) ，2

（1.河海大學(xué)水工結(jié)構(gòu)研究所，江蘇 南京 210098；2.河海大學(xué)水文水資源與水利水電工程科學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室，江蘇 南京 210098）

0 引 言

近年來我國(guó)高土石壩發(fā)展十分迅速，粗粒料作為土石壩的主要材料，其力學(xué)特性對(duì)壩體的變形有著決定性的影響。對(duì)粗粒料進(jìn)行二維數(shù)值試驗(yàn)時(shí)，孔隙率是影響粗粒料力學(xué)性質(zhì)的重要因素，對(duì)粗粒料的強(qiáng)度與變形均有很大的影響。為使二維數(shù)值試驗(yàn)中的二維模型與模擬的三維實(shí)體力學(xué)性質(zhì)等價(jià)，必須實(shí)現(xiàn)二維模型的孔隙率與三維實(shí)體的等價(jià)?，F(xiàn)今常用的二維與三維孔隙率轉(zhuǎn)換公式是由等粒徑顆粒體組成的結(jié)構(gòu)推導(dǎo)而出，僅用密實(shí)度進(jìn)行修正，未考慮粒徑分布對(duì)孔隙率的影響[1-3]，無法實(shí)現(xiàn)二維模型孔隙率和實(shí)體的三維孔隙率完全等價(jià)。鑒于此，本文參考不同土石壩工程的粗粒料級(jí)配，設(shè)計(jì)一系列具有代表性的級(jí)配組，對(duì)每組用二維顆粒流程序（pfc2d）和三維顆粒流程序 （pfc3d）模擬極值孔隙率試驗(yàn)，統(tǒng)計(jì)出不同級(jí)配下二維與三維極值孔隙率，探討粒徑分布及壓實(shí)程度對(duì)二維與三維孔隙率轉(zhuǎn)換關(guān)系的影響，擬合出二維與三維孔隙率的轉(zhuǎn)換公式，并與常用的轉(zhuǎn)換公式進(jìn)行比較。

1 數(shù)值試驗(yàn)的級(jí)配確定

土石壩為當(dāng)?shù)夭牧蠅?，受料?chǎng)原生條件的限制，不同工程所用粗粒料級(jí)配差異較大。為了使所選用的級(jí)配曲線組具有代表性，探尋不同工程粗粒料級(jí)配的共有特征，在采用重型振動(dòng)碾進(jìn)行薄層碾壓施工的工程中，粗粒料的級(jí)配設(shè)計(jì)均以talbot曲線[4]為母線。

參照水布埡和三板溪[5-6]等工程級(jí)配設(shè)計(jì)方法，以talbot曲線為母線設(shè)計(jì)試驗(yàn)級(jí)配。根據(jù)DL/T 5356—2006《水利水電工程粗粒土試驗(yàn)規(guī)程》，選定最大粒徑為60 mm，引入分形維數(shù)定義talbot曲線形狀指數(shù)[7]，設(shè)置分形維數(shù)n分別為2.7、2.6、2.5、2.4、2.3、2.2、2.1、2.0，所得的級(jí)配曲線組見圖1。曲線組構(gòu)成的區(qū)域幾乎涵蓋現(xiàn)代采用連續(xù)級(jí)配理論設(shè)計(jì)的堆石壩進(jìn)行數(shù)值試驗(yàn)的級(jí)配。因此，根據(jù)這組級(jí)配研究得到的規(guī)律適用于現(xiàn)代及今后按連續(xù)級(jí)配理論設(shè)計(jì)的堆石壩粗粒料的數(shù)值試驗(yàn)研究中。

圖1 各分形維數(shù)下的級(jí)配

2 數(shù)值試驗(yàn)

2.1 顆粒流程序模擬粗粒料極值孔隙率試驗(yàn)的可行性分析

由等粒徑顆粒結(jié)構(gòu)推導(dǎo)出的常用二維與三維孔隙率轉(zhuǎn)換關(guān)系為[3]

式中，ε2d為二維孔隙率；ε3d為三維孔隙率；ξ為修正系數(shù)；Dr為相對(duì)密度。

用pfc2d和pfc3d模擬等粒徑顆粒料的極值孔隙率試驗(yàn)，校驗(yàn)得出的二維與三維孔隙率是否符合公式（1），以此決定研究方案的可行性。為加快計(jì)算速度，pfc2d數(shù)值模擬中設(shè)置顆粒直徑為0.005 m、2000個(gè)；pfc3d中設(shè)置顆粒直徑為0.005 m、5000個(gè)。將pfc2d模擬極值孔隙率試驗(yàn)后得出的極值孔隙率按公式（1）轉(zhuǎn)換成三維孔隙率，與pfc2d、pfc3d模擬極值孔隙率試驗(yàn)得到的極值孔隙率統(tǒng)計(jì)于表1。表1中，按公式轉(zhuǎn)換得到的三維孔隙率與pfc3d模擬得到的孔隙率非常近似，證明用顆粒流程序模擬極值孔隙率試驗(yàn)探究二維與三維孔隙率轉(zhuǎn)換關(guān)系的方案可行。

表1 等粒徑顆粒組的二維、三維極值孔隙率

2.2 pfc模擬極值孔隙率試驗(yàn)

2.2.1 試驗(yàn)參數(shù)的確定

數(shù)值試驗(yàn)參數(shù)通過一組硅質(zhì)巖堆石料三軸試驗(yàn)曲線反算調(diào)整獲得[8-9]，數(shù)值試驗(yàn)的細(xì)觀參數(shù)見表2。其中，Kn、Ks分別為顆粒的法向及切向接觸剛度；f為顆粒間摩擦系數(shù)。

表2 數(shù)值試驗(yàn)的細(xì)觀參數(shù)

2.2.2 數(shù)值模擬步驟

（1）建立邊界:在pfc3d中建立底半徑為0.15 m，高為0.34 m的圓筒模型。在pfc2d中建立寬為0.3 m、高為0.34 m的矩形模型。

（2）制樣:在邊界中根據(jù)設(shè)計(jì)的級(jí)配曲線制樣，按照水利水電實(shí)驗(yàn)規(guī)程，粗粒料按60～40、40～20、20～10、10～5 mm以及5 mm以下進(jìn)行粒組篩分，但依此生成顆粒后發(fā)現(xiàn)pfc3d中顆粒數(shù)目達(dá)幾十萬個(gè)，遠(yuǎn)遠(yuǎn)超出電腦負(fù)荷。因此，將5 mm以下的顆粒全部用5～10 mm顆粒進(jìn)行等質(zhì)量替換。

（3）模擬松填法求最大孔隙率:顆粒在重力作用下自由下落，通過控制顆粒速度確保顆粒落距2～5 cm，顆粒下落穩(wěn)定后為最大孔隙率狀態(tài)。

（4）模擬振動(dòng)臺(tái)法求最小孔隙率:在顆粒表面置加重底板，使其對(duì)顆粒產(chǎn)生14 kPa的平均壓力。設(shè)定底部墻體運(yùn)動(dòng)方程為正弦曲線，振幅為0.002 m，頻率為251.3 rad/s，模擬振動(dòng)臺(tái)振動(dòng)，追蹤配位數(shù)變化，穩(wěn)定后停止振動(dòng)，顆粒達(dá)到最密實(shí)狀態(tài)。

（5）卸除加重底板，顆粒釋放應(yīng)力自由回彈后為最小孔隙率狀態(tài)。

2.2.3 模擬結(jié)果

根據(jù)公式（1），將pfc2d數(shù)值模擬得到的孔隙率換算成三維孔隙率，與pfc2d、pfc3d模擬得到的極值孔隙率統(tǒng)計(jì)于表3。

3 結(jié)果分析

3.1 數(shù)據(jù)有效性分析

從表3看出，三維孔隙率均大于二維孔隙率，這是由于二維模型中假定所有的顆粒球心均處在一個(gè)平面上，而實(shí)際三維實(shí)體中顆粒球心無規(guī)律分布于空間，更容易形成 “弓形”架構(gòu)。

表3 各級(jí)配的極值孔隙率

由于pfc3d運(yùn)算對(duì)顆粒數(shù)目的限制，制樣中簡(jiǎn)化了0～5 mm的顆粒，不可避免地給極值孔隙率值帶來一定誤差。表3中，分形維數(shù)由2.7減至2.0，極值孔隙率總體上呈先減后增的趨勢(shì)，在分形維數(shù)為2.5左右達(dá)到最小，這與Ochiai M等[10]得到的規(guī)律相吻合。本文旨在探究粗粒料的二維與三維極值孔隙率轉(zhuǎn)換關(guān)系，pfc2d與pfc3d數(shù)值試驗(yàn)中粗粒料級(jí)配的簡(jiǎn)化是同步的，且所得數(shù)據(jù)符合實(shí)際規(guī)律，故得到的轉(zhuǎn)換關(guān)系是可靠的。

3.2 擬合公式

繪制pfc數(shù)值試驗(yàn)得到的三維孔隙率關(guān)于二維孔隙率的散點(diǎn)圖 （見圖2），將二維極值孔隙率按公式 （1）換算得到三維孔隙率，在圖2中繪入公式換算得到的三維孔隙率與二維孔隙率對(duì)應(yīng)點(diǎn)。圖2中，按照公式換算得到的三維孔隙率較pfc3d模擬得到的均偏大，且在Dr=0.0時(shí)相差較大。這是由于公式 （1）由等粒徑顆粒體推導(dǎo)而來，僅用相對(duì)密度進(jìn)行修正，忽略了粒徑分布對(duì)孔隙率的影響，而級(jí)配良好的顆粒堆積體孔隙率明顯小于級(jí)配差的，故按照公式計(jì)算的孔隙率均大于pfc3d模擬得到的。

對(duì)不同分形維數(shù)下二維與三維極值孔隙率進(jìn)行相關(guān)性分析可知，二維與三維極值孔隙率呈高度線性相關(guān)，相關(guān)系數(shù)達(dá)0.94，可用直線y=1.4618x+0.08794擬合，并繪制擬合直線于圖2中。顯然，相比于公式 （1），該直線能更好地反應(yīng)粗粒料二維與三維極值孔隙率的關(guān)系，且相對(duì)密度分別為0.0與1.0時(shí)，二維與三維孔隙率的對(duì)應(yīng)點(diǎn)分別密集分布在直線的 x=0.12～0.14、0.18～0.20段。由此推測(cè)，其他相對(duì)密度下粗粒料二維與三維孔隙率對(duì)應(yīng)點(diǎn)會(huì)密集分布在0.14～0.18段，但需進(jìn)一步證實(shí)。

圖2 二維與三維孔隙率對(duì)應(yīng)點(diǎn)

4 結(jié) 論

通過設(shè)計(jì)一系列具有代表性的粗粒料級(jí)配，用pfc2d、pfc3d模擬極值孔隙率試驗(yàn)探究粗粒料二維與三維孔隙率的對(duì)應(yīng)關(guān)系，得出如下結(jié)論:

（1）實(shí)際三維結(jié)構(gòu)比二維模型的孔隙率大。二維模型假定所有的顆粒球心均處在一個(gè)平面上，而三維結(jié)構(gòu)中顆粒球心無規(guī)律分布于空間中，更易形成 “弓形”架構(gòu)，故孔隙更大。

（2）相同壓實(shí)度下，等粒徑顆粒體比級(jí)配良好的顆粒體的孔隙率大。根據(jù)常用二維與三維孔隙率轉(zhuǎn)換公式求得的三維孔隙率會(huì)大于實(shí)際三維孔隙率，且在相對(duì)密度為0時(shí)相差較大。

（3）設(shè)計(jì)級(jí)配下，粗粒料的二維極值孔隙率與三維極值孔隙率呈高度線性相關(guān)，用直線擬合后具有良好的歸一性。相對(duì)密度分別為0.0與1.0時(shí)，二維與三維孔隙率的對(duì)應(yīng)點(diǎn)分別密集分布在直線的兩段，推測(cè)其他密實(shí)度下二維與三維孔隙率對(duì)應(yīng)點(diǎn)也會(huì)密集分布在該直線的某段，但需進(jìn)一步研究證實(shí)。

（4）本次數(shù)值試驗(yàn)探究粗粒料的二維與三維孔隙率對(duì)應(yīng)關(guān)系，沒有考慮粗粒料的顆粒形狀和破碎的影響，考慮顆粒形狀和破碎后粗粒料的二維與三維孔隙率對(duì)應(yīng)關(guān)系尚待進(jìn)一步探究。

［1］MITCHEL J K.巖土工程土性分析原理［M］.高國(guó)瑞，譯.南京:南京工學(xué)院出版社，1988.

［2］日本土質(zhì)工學(xué)會(huì).粗粒料的現(xiàn)場(chǎng)壓實(shí)［M］.北京:水利水電出版社，1999.

［3］張剛.管涌現(xiàn)象細(xì)觀機(jī)理的模型試驗(yàn)與顆粒流數(shù)值模擬研究［D］.上海:同濟(jì)大學(xué)，2007.

［4］王永明，朱晟，任金明，等.筑壩粗粒料力學(xué)特性的縮尺效應(yīng)研究［J］.巖土力學(xué)，2013，34（6）:1799-1807.

［5］李昌彩.水布埡面板堆石壩前期關(guān)鍵技術(shù)研究［M］.北京:水利水電出版社，2005.

［6］趙海斌.超硬巖筑壩材料現(xiàn)場(chǎng)爆破與現(xiàn)場(chǎng)碾壓試驗(yàn)研究［R］.長(zhǎng)沙:中國(guó)水電顧問集團(tuán)中南勘測(cè)設(shè)計(jì)研究院，2006.

［7］朱晟，武利強(qiáng)，魏匡民，等.一種基于分形理論的粗粒料縮尺方法: 中國(guó)，CN103134906A［P］.2013-06-05.

［8］王永明，任金明，武利強(qiáng)，等.pfc模擬筑壩粗粒料力學(xué)特性的有效性分析［J］.水電能源科學(xué)，2012，30（10）:72-75.

［9］周健，池毓蔚，池永，等.砂土雙軸試驗(yàn)的顆粒流模擬［J］.巖土工程學(xué)報(bào)，2000，22（6）:701-704.

［10］OCHIAI M,OZAO R,YAMAZAKI Y,et al.Self-similarity law of particle size distribution and energy law in size reduction of solids［J］.Physica A:Statistical Mechanics and its Applications,1992,191（2）:295-300.