非穩(wěn)定流作用下管涌發(fā)生發(fā)展的細(xì)觀數(shù)值模型試驗(yàn)研究

倪小東，趙帥龍，王媛

?

非穩(wěn)定流作用下管涌發(fā)生發(fā)展的細(xì)觀數(shù)值模型試驗(yàn)研究

倪小東1, 2，趙帥龍1, 2，王媛1, 2

(1. 巖土力學(xué)與堤壩工程教育部重點(diǎn)實(shí)驗(yàn)室，江蘇南京，210098；2. 河海大學(xué)隧道與地下工程研究所，江蘇南京，210098)

為了揭示非穩(wěn)定流作用管涌發(fā)生后的發(fā)展過程及其破壞規(guī)律，進(jìn)行管涌砂槽模型試驗(yàn)和顆粒流數(shù)值模擬，建立6組模型研究水頭抬升速率對(duì)顆粒運(yùn)移過程的影響，分析不同條件下土體內(nèi)部細(xì)顆粒運(yùn)移引起管涌破壞的動(dòng)態(tài)過程。根據(jù)數(shù)值模擬結(jié)果，揭示管涌發(fā)展過程中顆粒細(xì)觀變化規(guī)律及流體的變化規(guī)律，并與砂槽模型試驗(yàn)結(jié)果進(jìn)行比較。研究結(jié)果表明：一次加載方式對(duì)孔口區(qū)域的影響最大，也最不安全，此結(jié)論與工程實(shí)際相吻合，模型試驗(yàn)結(jié)果證實(shí)了該數(shù)值模型的合理性。

非穩(wěn)定流；管涌；離散元；破壞水頭；滲透系數(shù)

管涌是土體在滲流作用下被侵蝕為可動(dòng)細(xì)顆粒，并隨水在孔隙中運(yùn)移流失的過程。在此過程中，土體中的粗顆?？赡鼙患芸?、塌落，最后造成土體的破 壞[1]。近年來，國(guó)內(nèi)外學(xué)者針對(duì)管涌機(jī)理進(jìn)行了大量的研究，提出了不同的模型和計(jì)算方法：RICHARDS等[2]研制了真三軸管涌試驗(yàn)裝置，研究了影響管涌發(fā)展的重要因素；SELLMEIJIER[3?5]給出了單一砂層堤基上臨界水頭的理論解答；毛昶熙等[6]給出了管涌險(xiǎn)情是否會(huì)影響大堤安全的估算公式，但還只能適用于單一粉細(xì)砂的堤基結(jié)構(gòu)；周健[7]從細(xì)觀角度揭示了管涌發(fā)展過程中土體幾何特性和水力特性的復(fù)雜變化，但其簡(jiǎn)化了土體成分，與實(shí)際土體差別較大；STERPI[8]根據(jù)可動(dòng)細(xì)顆粒的連續(xù)方程和滲流方程，建立了研究可動(dòng)細(xì)顆粒流失過程的有限差分模型，雖然考慮了滲流對(duì)土體內(nèi)可動(dòng)細(xì)顆粒流失的影響，但并未考慮可動(dòng)細(xì)顆粒流失引起的土體滲透性變化；SHAMY等[9]通過顆粒流方法研究了砂土液化機(jī)制；倪小東等[10]基于顆粒流方法，考慮顆粒相與液相之間的相互作用力，從細(xì)觀層面研究了管涌發(fā)生的機(jī)理。以上研究大多基于穩(wěn)定流，而非穩(wěn)定滲流作用下土體發(fā)生管涌破壞的研究，對(duì)了解土體的滲透變形特性非常重要。毛昶熙等[11]進(jìn)行了非穩(wěn)定滲流模型試驗(yàn)，給出了管涌險(xiǎn)情發(fā)展的時(shí)間參考值；王媛等[12]研究了巖體水力劈裂非穩(wěn)定滲流的影響機(jī)制；葉祖洋等[13]將濕區(qū)上的非穩(wěn)定滲流問題轉(zhuǎn)化為全域上的一個(gè)新的初邊值問題對(duì)巖體裂隙網(wǎng)絡(luò)非穩(wěn)定滲流加以研究；張健等[14]基于球狀井非穩(wěn)定井流理論對(duì)管涌流場(chǎng)分布進(jìn)行了研究；陳益峰等[15]建立了拋物型變分不等式(PVI)提法，研究了Signorini類型邊界條件下非穩(wěn)定滲流問題；陶同康[16]基于砂模型試驗(yàn)，探討了不透水地基均質(zhì)土壩非穩(wěn)定滲流的物理力學(xué)過程；胡冉等[17]從穩(wěn)定滲流分析拓展到非穩(wěn)定滲流分析；張培文等[18]介紹了飽和?非飽和滲流的計(jì)算程序數(shù)值模擬方法；柳厚祥等[19]根據(jù)尾礦壩形成的特點(diǎn)，進(jìn)行考慮應(yīng)力場(chǎng)與滲流場(chǎng)耦合的非穩(wěn)定滲流分析。然而，大多數(shù)研究只是從表觀破壞的角度切入，并沒有從破壞發(fā)展的時(shí)間角度進(jìn)行研究，并且土體在非穩(wěn)定流作用下的管涌特性還沒有得到系統(tǒng)的描述。因此，有必要對(duì)非穩(wěn)定流作用下土體管涌發(fā)生發(fā)展的規(guī)律特性進(jìn)行研究。本文作者通過數(shù)值模擬，研究不同水頭條件下土體管涌發(fā)生發(fā)展的規(guī)律。采用細(xì)觀離散數(shù)值模擬方法比對(duì)室內(nèi)模型試驗(yàn)共同探究水頭抬升對(duì)管涌的影響。將細(xì)觀數(shù)值模擬獲得的管涌臨界水力梯度與室內(nèi)試驗(yàn)以及相關(guān)理論結(jié)果進(jìn)行比較，為非穩(wěn)定流作用下的管涌發(fā)生發(fā)展提供了一種新的研究手段。

1 砂槽模型試驗(yàn)

1.1 試驗(yàn)方案

為研究單層堤基管涌破壞過程及機(jī)理，設(shè)計(jì)了長(zhǎng)1.3 m、寬0.5 m、高0.9 m的試驗(yàn)?zāi)Ｐ筒?，其中進(jìn)水室高0.8 m，砂樣長(zhǎng)1 m，寬 0.5 m，高 0.5 m。模型槽進(jìn)水系統(tǒng)、觀測(cè)系統(tǒng)(模型槽的側(cè)壁和底部以及上覆有機(jī)玻璃板中線位置分別開孔安裝測(cè)壓管) 見圖1。為了避免發(fā)生接觸沖刷，在砂層和蓋板之間預(yù)設(shè)5 mm厚填充層，其材質(zhì)為不透水且具有較大可塑性的油泥，油泥外包薄膜并在上下表面涂抹凡士林以阻水。在中心線上距進(jìn)水口0.775，0.875和0.975 m處分別預(yù)留直徑為5 cm的管涌口，試驗(yàn)前將所有孔口封閉，進(jìn)行不同滲徑長(zhǎng)度試驗(yàn)時(shí)分別打開相應(yīng)的孔口，同時(shí)用膠塞封閉其他孔口。試樣采用水下分層拋填的方式裝填，頂面抹平，鋪設(shè)填充層，蓋上有機(jī)玻璃板，用螺桿壓緊，接縫處用玻璃膠密封。為考慮水頭抬升對(duì)流場(chǎng)的影響，首先確定模型發(fā)生破壞的臨界狀態(tài)后，以該水壓力作為臨界水壓力，模擬4種不同加載方式下各工況的情況。具體工況參數(shù)見表1。

單位：m

表1 試驗(yàn)工況

1.2 砂樣選擇與物理參數(shù)

在砂模型制樣過程中，首分別代表骨料顆粒和填充顆粒，粗顆粒由粒徑范圍1~5 mm的粗砂及細(xì)礫混合料組成，細(xì)顆粒由粒徑范圍0.1~0.5 mm的中細(xì)砂粒組成。粗、細(xì)砂物性指標(biāo)見表2，粗、細(xì)砂級(jí)配曲線如圖2所示。

表2 試驗(yàn)用粗、細(xì)砂物性指標(biāo)

1—粗砂級(jí)配曲線；2—細(xì)砂級(jí)配曲線。

1.3 試驗(yàn)方法及結(jié)果

試驗(yàn)水頭邊界模擬采用逐級(jí)抬升水頭方式進(jìn)行，待測(cè)壓管水位穩(wěn)定一段時(shí)間后將水頭抬升至下一級(jí)，直至滲流無法保持穩(wěn)定。滲流穩(wěn)定的判別標(biāo)準(zhǔn)是：滲流量和測(cè)壓管水位基本穩(wěn)定，管涌口水流清澈且沒有砂粒帶出。

分別針對(duì)細(xì)料質(zhì)量分?jǐn)?shù)為10%，20%，25%的3組試樣進(jìn)行分析，同時(shí)考慮孔隙率不同情況，將細(xì)料質(zhì)量分?jǐn)?shù)不同的砂樣編號(hào)為Ⅰ，Ⅱ，Ⅲ，孔隙率不同的砂樣編號(hào)為1，2，得到各組砂樣參數(shù)見表3。文中室內(nèi)試驗(yàn)及數(shù)值模型體積分?jǐn)?shù)均選用細(xì)料質(zhì)量分?jǐn)?shù)為25%的模型，其砂樣編號(hào)為Ⅲ-1。

表3 制備砂樣參數(shù)

實(shí)際工程中，管涌型土最優(yōu)細(xì)料質(zhì)量分?jǐn)?shù)一般比較接近于25%，因此，表4給出了細(xì)料質(zhì)量分?jǐn)?shù)為25%時(shí)臨界水力梯度。

表4 臨界水力梯度

水頭在抬升的各個(gè)局部階段可以看作近似屬于穩(wěn)定流的范疇，將各個(gè)水頭抬升階段組合起來，以反映非穩(wěn)定流作用下模型內(nèi)部管涌發(fā)生發(fā)展的規(guī)律。由表4可以看出：滲透破壞時(shí)隨著加載數(shù)的增加各模型中的水力梯度逐漸升高。

2 數(shù)值模型的建立與工況

2.1 細(xì)觀模型與宏觀模型匹配

為降低顆粒模型中總的自由度，節(jié)約計(jì)算資源，需降低細(xì)顆粒生成數(shù)目，因此考慮構(gòu)建與室內(nèi)模型試驗(yàn)相似的細(xì)觀顆粒模型。參照離心試驗(yàn)相似原理，結(jié)合多孔介質(zhì)流體運(yùn)動(dòng)相似準(zhǔn)則[20]，將室內(nèi)模型和粒徑縮放為原型的1/10，標(biāo)識(shí)為相似模型Ⅰ，施加10重力場(chǎng)，并將流體黏滯系數(shù)取1/10，則可實(shí)現(xiàn)離心相似(雷諾相似及弗洛德相似)；通過室內(nèi)試驗(yàn)量測(cè)及分析，獲知模型中雷諾數(shù)接近于0.12，采用與室內(nèi)模型試驗(yàn)縮尺相反的技術(shù)手段，將相似模型I中模型尺寸不變，粒徑統(tǒng)一放大10倍，同時(shí)將黏滯系數(shù)變?yōu)樵偷恼{(diào)整為初始值的倍，得出滿足雷諾相似的相似模型Ⅱ，此為數(shù)值顆粒模型建立依據(jù)。

在PFC3D中，基于相似模型Ⅱ方案，為與室內(nèi)模型試驗(yàn)孔口1的結(jié)果進(jìn)行比較，基于PFC3D建立如圖3所示模型，模型長(zhǎng)×寬×高為0.10 m×0.05 m×0.05 m。室內(nèi)模型中骨架顆粒與填充顆粒粒徑差別較大，如按顆分曲線在PFC3D中生成骨架與填充顆粒，需要耗費(fèi)巨大的計(jì)算資源，甚至無法完成顆粒生成工作。根據(jù)反濾層設(shè)計(jì)相關(guān)知識(shí)可知，骨架顆粒與填充顆粒平均粒徑之比大于5時(shí)，填充顆粒在骨架孔隙中能夠自由流動(dòng)，因此本文數(shù)值模型中骨架顆粒與填充顆粒的粒徑比取為5。粗細(xì)顆粒粒徑分別為0.5 mm和0.1 mm，流體計(jì)算單元取邊長(zhǎng)為1 cm的立方體，，，方向單元數(shù)分別為10，5，5。建立與室內(nèi)試驗(yàn)相對(duì)應(yīng)的6組試樣，具體參數(shù)見表5。

表5 數(shù)值試驗(yàn)材料參數(shù)

圖3 數(shù)值模型示意圖

采用PFC3D顆粒隨機(jī)生成器先產(chǎn)生粗顆粒，后在剩余區(qū)域產(chǎn)生細(xì)顆粒。施加水頭邊界，左側(cè)為上游水頭邊界，孔口為自由面邊界，砂槽側(cè)壁和底部均設(shè)為剛性不透水非滑移邊界，模型具體參數(shù)按表6選取[21]，滲透破壞示意圖如圖4所示。將=10 m，=3 m，=5 m處單元設(shè)置為出流孔口，滲透破壞示意圖如圖4所示，模型左側(cè)為水頭邊界，孔口為零壓力邊界，其他邊界均設(shè)為剛性不透水非滑移邊界。

表6 材料特性參數(shù)

圖4 滲透破壞模型示意圖

2.2 臨界狀態(tài)的確定

近似采用滿足達(dá)西流狀態(tài)的雷諾相似準(zhǔn)則建立模型[22?23]，為了研究試樣在非穩(wěn)定流作用下的管涌破壞機(jī)理，首先確定模型發(fā)生破壞的臨界狀態(tài)。模型發(fā)生滲透破壞后，存在一基本特征，如有大量細(xì)顆粒從內(nèi)部流出，并且短時(shí)間內(nèi)，這種現(xiàn)象不會(huì)隨著時(shí)間的延長(zhǎng)而消失，直至模型內(nèi)部形成管涌通道，流量或顆粒流失量穩(wěn)定在某一數(shù)值。

2.2.1 根據(jù)孔口流量判斷

模型在各種承壓水頭下孔口流量隨時(shí)間變化曲線如圖5所示，由圖5可以看出：在水頭較小時(shí)，流量曲線尚能保持為一定值，或者在一定的區(qū)間內(nèi)波動(dòng)，當(dāng)水壓力為38 MPa時(shí)，模型中孔口流量仍能保持穩(wěn)定值，當(dāng)超過該值時(shí)，流量曲線開始呈現(xiàn)上揚(yáng)趨勢(shì)，說明模型發(fā)生滲透破壞。

水壓力/MPa：1—22；2—32；3—38；4—39；5—40。

2.2.2 根據(jù)孔口顆粒流失量判斷

孔口顆粒流失量隨時(shí)間變化關(guān)系曲線如圖6所示。由圖6可以看出：在水壓力較小時(shí)，孔口顆粒流失量隨時(shí)間變化趨勢(shì)逐漸平緩，直至成水平狀，也即不再有顆粒流出，認(rèn)定此時(shí)模型仍能保持穩(wěn)定狀態(tài)，水壓力超過38 MPa后，孔口顆粒流失量隨時(shí)間持續(xù)增加，判定為模型發(fā)生破壞。

水壓力/MPa：1—22；2—32；3—38；4—39；5—40。

綜合圖5和圖6分析結(jié)果可得：模型的臨界水壓力為38 MPa。

2.3 數(shù)值模擬工況

確定了模型發(fā)生破壞的臨界狀態(tài)，以該水頭作為臨界水頭，然后在數(shù)值模擬中對(duì)供水水頭加以控制，分別模擬4種工況，通過不同的方式逐步升高上游水頭使模型發(fā)生管涌破壞，模型共持續(xù)64 s，具體工況參數(shù)見表7。

表7 數(shù)值模擬工況參數(shù)

2.4 室內(nèi)試驗(yàn)、數(shù)值模擬及理論公式比較

數(shù)值模擬獲得的臨界水力梯度與室內(nèi)試驗(yàn)獲得的臨界值比較見圖7。

圖7 室內(nèi)模型與數(shù)值模型臨界水力梯度

由圖7可知：數(shù)值模擬獲得的結(jié)果和室內(nèi)試驗(yàn)獲得的結(jié)果比較接近，遵循加載數(shù)越多，壓力分布調(diào)整的幅度越大、發(fā)生管涌破壞的臨界水頭越小的規(guī)律?，F(xiàn)有研究中，有關(guān)無黏性土中非穩(wěn)定流管涌臨界水力梯度研究尚無相應(yīng)的經(jīng)驗(yàn)公式，因此，將試驗(yàn)結(jié)果與沙金煊公式[24]進(jìn)行比較，可見結(jié)果比較接近，說明采用顆粒流方法對(duì)管涌進(jìn)行研究是可行的，從另一方面也說明所建立的水頭抬升的管涌模型及其計(jì)算結(jié)果是正確的。

3 模擬結(jié)果與影響因素分析

確定模型發(fā)生破壞的臨界狀態(tài)后，以該水壓力作為臨界水壓力，模擬4種不同加載方式下各工況的情況。根據(jù)實(shí)際情況可知，在水頭升高過程中，一次加載可近似看作處于穩(wěn)定流狀態(tài)，多次加載可以看作具有一定非穩(wěn)定程度的非穩(wěn)定流，以反映非穩(wěn)定流作用下模型內(nèi)部管涌發(fā)生發(fā)展的規(guī)律。

3.1 不同加載方式下單元孔隙率變化

不同加載方式下孔口附近及上游單元孔隙率隨時(shí)間的變化曲線如圖8所示。由圖8可以看出：隨著時(shí)間的增加，模型各單元孔隙率呈現(xiàn)明顯的增大趨勢(shì)，孔隙率隨著加載數(shù)增加而減小，單元孔隙率隨著離孔口距離的增加而變大，且多次加載與一次加載條件下單元孔隙率的差值隨著離孔口距離的增加而變大。

(a) x方向第8個(gè)單元；(b) x方向第9個(gè)單元；(c) x方向第10個(gè)單元1—一次加載；2—二次加載；3—三次加載；4—四次加載。

3.2 相同加載方式下不同單元孔隙率變化

相同加載方式下不同孔口附近及上游單元孔隙率隨時(shí)間的變化曲線如圖9所示。由圖9可以看出：模型各單元孔隙率隨著加載數(shù)增加而減小。在圖9(a)中，方向第8個(gè)單元與第9個(gè)、第10個(gè)單元的孔隙率變化趨勢(shì)幾乎一致，而在圖9(b)~(d)中，該單元孔隙率變化趨勢(shì)完全不同，特別是在圖9(d)中，水頭較小時(shí)，孔隙率幾乎沒有變化，水頭增加到一定值時(shí)，孔隙率才發(fā)生明顯的變化。

(a) 一次加載；(b) 二次加載；(c) 三次加載；(d) 四次加載1—x方向第8個(gè)單元；2—x方向第9個(gè)單元；3—x方向第10個(gè)單元。

根據(jù)匯流理論，水力梯度在上游水頭較小時(shí)，尚不足以帶動(dòng)顆粒發(fā)生運(yùn)移，隨著上游水頭增加，水力梯度相應(yīng)增加，顆粒開始啟動(dòng)。一級(jí)加載作用下距離孔口較遠(yuǎn)處的顆粒也可以開始起動(dòng)并逐漸匯向孔口，引起孔口細(xì)料大量聚集，水力梯度激增，直至突破臨界狀態(tài)；而逐級(jí)加載時(shí)細(xì)料起動(dòng)范圍較小，壓力場(chǎng)可以隨著逐級(jí)加載不斷得到調(diào)整，因此加載至臨界水頭時(shí)，一級(jí)加載比多級(jí)加載危險(xiǎn)。

3.3 孔口流量隨時(shí)間變化

分級(jí)加載時(shí)孔口流量隨時(shí)間的變化曲線如圖10所示。由圖10可知：加載分級(jí)數(shù)越多，最終孔口流量越小。相同條件下孔口流量較小主要是由于孔口滲透系數(shù)較小或者水力梯度較小，此處假定是由于孔口周邊滲透系數(shù)減小而引起流量的減小?？梢娝鲿?huì)對(duì)土體的滲透性造成影響，并且這種影響的結(jié)果，偏向有利于土體滲透性減小的方向，與相關(guān)研究[25]所得結(jié)論一致。一定意義上，所測(cè)得的滲透規(guī)律可作為非穩(wěn)定流過程中土體滲透性研究的依據(jù)。

1—一次加載；2—二次加載；3—三次加載；4—四次加載。

3.4 孔口顆粒流失量隨時(shí)間變化

孔口顆粒流失量隨時(shí)間變化曲線如圖11所示。由圖11可以看出：在不同加載方式的作用下，顆粒流失量最終均能達(dá)到穩(wěn)定狀態(tài)，顆粒流失量達(dá)到穩(wěn)定值所需要的時(shí)間隨著加載分級(jí)數(shù)的增加而變大，最終的顆粒流失量隨著加載分級(jí)數(shù)的增加而減小，說明土樣內(nèi)部的水流速度以及土樣內(nèi)部土顆粒所受到的水流作用力隨著加載分級(jí)數(shù)的增加而減小。主要原因在于，分級(jí)加載過程中平均水力梯度較低，孔口附近區(qū)域由于空間匯流作用，孔口數(shù)倍直徑范圍內(nèi)水力梯度仍能使顆粒發(fā)生運(yùn)移。

1—一次加載；2—二次加載；3—三次加載；4—四次加載。

以四次加載曲線為例，雖然第1次加載僅施加10 MPa壓力，但是孔口仍然流失了一定量的顆粒，因此，整個(gè)模型內(nèi)的孔壓發(fā)生重新分布。由于已經(jīng)流失了一定的顆粒，孔口影響區(qū)域在施加下一級(jí)荷載時(shí)承擔(dān)的水頭較直接加載時(shí)承擔(dān)的水頭小，孔口附近區(qū)域的水力梯度相應(yīng)減小，顆粒受到的滲透力減小，水流能夠帶出孔口的顆粒相應(yīng)減小。

圖10和圖11模擬結(jié)果顯示：直接加載至臨界水頭時(shí)對(duì)孔口流量及顆粒流失量的作用較分級(jí)加載時(shí)作用更為明顯。

4 結(jié)論

1) 管涌試驗(yàn)顆粒流模擬結(jié)果和室內(nèi)試驗(yàn)比較接近，說明采用細(xì)觀離散數(shù)值方法可以很好地模擬砂槽模型試驗(yàn)。通過顆粒流模擬可以更直觀地分析室內(nèi)試驗(yàn)不易得到的域內(nèi)顆粒運(yùn)移特征及各單元孔隙率變化特征，可以得到顆粒的運(yùn)動(dòng)狀態(tài)和應(yīng)力狀態(tài)，避免采用宏觀方法研究所面臨的不確定性。

2) 在臨界水頭相同，一次加載至臨界狀態(tài)，模型孔口區(qū)域的孔隙率最大，加載數(shù)越大，孔口附近及上游單元孔隙率越小；單位時(shí)間內(nèi)的水流量越?。活w粒流失量達(dá)到穩(wěn)定值的時(shí)間越久；最終顆粒流失量越??；土樣內(nèi)部的水流速度以及土樣內(nèi)部土顆粒所受到的水流作用力也越小。

3) 在臨界水頭相同，逐級(jí)加載時(shí)細(xì)料起動(dòng)范圍較小，且在逐級(jí)加載過程中，壓力場(chǎng)可以不斷得到調(diào)整，加載至臨界狀態(tài)時(shí)，逐級(jí)加載較一級(jí)加載安全。

4) 顆粒流模擬管涌試驗(yàn)獲得的臨界水力梯度值和室內(nèi)試驗(yàn)以及理論結(jié)果比較接近，表明采用數(shù)值模擬方法可以很好地模擬砂槽模型試驗(yàn)。同時(shí)顆粒流程序由于其本身的特點(diǎn)，可以分析理論較難涉及的非穩(wěn)定流的砂土管涌現(xiàn)象。因此采用顆粒流對(duì)管涌試驗(yàn)進(jìn)行數(shù)值模擬可以克服室內(nèi)試驗(yàn)面臨的填料不均、重復(fù)性不易實(shí)現(xiàn)以及理論研究暫時(shí)無法實(shí)現(xiàn)的困難。

[1] 龔曉南. 土力學(xué)[M]. 北京: 中國(guó)建筑工業(yè)出版社, 2005: 32?59.GONG Xiaonan. Soil mechanics[M]. Beijing: China Building Industry Press, 2005: 32?59.

[2] RICHARDS K S, REDDY K R. True triaxial piping test apparatus for evaluation of piping potential in earth structures[J]. Geotechnical Testing Journal, 2010, 33(1): 1?13.

[3] SELLMEIJER H, LOPEZ DE LA CRUZ J, Van BEEK V M, et al. Fine-tuning of the piping model through small-scale, medium-scale and Ijkdijk experiments[J]. European Journal of Environmental and Civil Engineering, 2011, 15(8): 1139–1154.

[4] KOENDERS M A, SELLMEIJER J B. Liquefaction of unstable slopes[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 2005, 16(8): 571?579.

[5] Van BEEK V M, BEZUIJEN A, SELLMEIJERJ B, et al. Initiation of backward erosion piping in uniform sands[J]. Géotechnique, 2014, 64(12): 927?941.

[6] 毛昶熙, 段祥寶, 蔡金傍, 等. 懸掛式防滲墻控制管涌發(fā)展的試驗(yàn)研究[J]. 水利學(xué)報(bào), 2005, 36(1): 42?50. MAO Changxi, DUAN Xiangbao, CAI Jinbang, et al. Experimental study on piping development control by means of suspended cut-off wall[J]. Journal of Hydraulic Engineering, 2005, 36(1): 42?50.

[7] 周健, 姚志雄, 張剛, 等. 砂土滲流過程的細(xì)觀數(shù)值模擬[J]. 巖土工程學(xué)報(bào), 2007, 29(7): 977?981. ZHOU Jian, YAO Zhixiong, ZHANG Gang. Mesomechanical simulation of seepage flow in sandy soil[J]. Chinese Journal of Geotechnical Engineering, 2007, 29(7): 977?981.

[8] STERPI D. Effects of the erosion and transport of fine particles due to seepage flow[J]. International journal of Geomechanics, 2003, 3(1): 111?122.

[9] SHAMY U E, AYDIN F. Multiscale modeling of flood-induced piping in river levees[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2008, 134(9): 1385?1398.

[10] 倪小東, 王媛, 王飛, 等. 管涌的砂槽試驗(yàn)研究及顆粒流模擬[J]. 四川大學(xué)學(xué)報(bào)(工程科學(xué)版), 2009, 41(6): 51?57. NI Xiaodong, WANG Yuan, WANG Fei, et al. Study on piping by sand-bank model and simulation by PFC3D[J]. Journal of Sichuan University (Engineering Science Edition), 2009, 41(8): 51?57.

[11] 毛昶熙, 段祥寶, 蔡金傍, 等. 洪峰過程非穩(wěn)定滲流管涌試驗(yàn)研究與理論分析[J]. 水利學(xué)報(bào), 2005, 36(9): 1105?1120. MAO Changxi, DUAN Xiangbao, CAI Jinbang, et al. Piping experimental study and theoretical analysis of unsteady seepage flow during flood peak[J]. Journal of Hydraulic Engineering, 2005, 36(9): 1105?1120.

[12] 王媛, 顏青青. 巖體水力劈裂非穩(wěn)定滲流影響機(jī)制初探[J]. 巖石力學(xué)與工程學(xué)報(bào), 2012, 31(10): 2016?2021. WANG Yuan, YAN Qingqing. Primary research of influence of unsteady seepage processes on hydraulic fracturing in rock mass[J]. Chinese Journal of Rock Mechanics and Engineering, 2012, 31(10): 2016?2021.

[13] 葉祖洋, 姜清輝, 姚池, 等. 巖體裂隙網(wǎng)絡(luò)非穩(wěn)定滲流分析與數(shù)值模擬[J]. 巖土力學(xué), 2013, 34(4): 1171?1190. YE Zuyang, JIANG Qinghui, YAO Chi, et al. Formulation and simulation of non-steady seepage flow through fracture network in rock masses[J]. Rock and Soil Mechanics, 2013, 34(4): 1171?1190.

[14] 張健, 翟劍峰, 王仙美, 等. 基于球狀井非穩(wěn)定井流理論的管涌流場(chǎng)分布研究[J]. 水利水運(yùn)工程學(xué)報(bào), 2014(1): 49?55. ZHANG Jian, ZHAI Jianfeng, WANG Xianmei, et al. Analysis of piping field distributions based on unsteady spherical well flow theory[J]. Hydro-Science and Engineering, 2014(1): 49?55.

[15] CHEN Yifeng, HU Ran, ZHOU Chuangbing, et al. A new parabolic variational inequality formulation of signorini’s condition for non-steady seepage problems with complex seepage control systems[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 2011, 35: 1034?1058.

[16] 陶同康. 堤壩非穩(wěn)定滲流試驗(yàn)研究與計(jì)算[J]. 巖土工程學(xué)報(bào), 1981, 3(1): 81?93. TAO Tongkang. Experiment study and computation of unsteady seepage in embankment dams[J]. Chinese Journal of Geotechnical Engineering, 1981, 3(1): 81?93.

[17] 胡冉, 陳益峰, 周創(chuàng)兵, 等. 非穩(wěn)定滲流問題的變分不等式方法及工程應(yīng)用[J]. 水動(dòng)力學(xué)研究與進(jìn)展, 2011, 26(2): 239?251. HU Ran, CHEN Yifeng, ZHOU Chuangbing, et al. Avariational inequality approach for non-steady seepage problems and its application in engineering practices[J]. Chinese Journal of Hydrodynamics, 2011, 26(3): 239?251.

[18] 張培文, 劉德富, 黃達(dá)海, 等. 飽和?非飽和非穩(wěn)定滲流數(shù)值模擬[J]. 巖土力學(xué), 2003, 24(6): 927?930. ZHANG Peiwen, LIU Defu, HUANG Dahai, et al. Saturated-unsaturated unsteady seepage flow numerical simulation[J]. Rock and Soil Mechanics, 2003, 24(6): 927?930.

[19] 柳厚祥, 李寧, 廖雪, 等. 考慮應(yīng)力場(chǎng)與滲流場(chǎng)耦合的尾礦壩非穩(wěn)定滲流分析[J]. 巖石力學(xué)與工程學(xué)報(bào), 2004, 23(17): 2870?2875. LIU Houxiang, LI Ning, LIAO Xue, et al. Unsteady seepage analysis of tailing dams considering coupling of stress and seepage fields[J] . Chinese Journal of Rock Mechanics and Engineering, 2004, 23(17): 2870?2875.

[20] WANG Yuan, NI Xiaodong. Hydro-mechanical analysis of piping erosion based on similarity criterion at micro-level by PFC3D[J]. European Journal of Environmental and Civil Engineering, 2013, 17(S1): 187?204.

[21] 倪小東. 土體細(xì)觀滲透破壞顆粒流理論與試驗(yàn)研究[D]. 南京: 河海大學(xué)土木工程學(xué)院, 2010: 86?87. NI Xiaodong. Study on seepage failure in meso-level by particle flow code (PFC) and model test[D]. Nanjing: Hohai University. School of Civil Engineering, 2010: 86?87.

[22] DEAN E T R. Discussion on scale-modeling of fluid flow in geotechnical centrifuges[J]. Soils and Foundations, 2001, 41(4): 108?110.

[23] THUSYANTHAN I, MADABHUSHI S. Scaling of seepage flow velocity in centrifuge models[J]. Acta Gastroenterologica Latinoamericana, 2008, 38(2): 105?115.

[24] 沙金煊. 多孔介質(zhì)中的管涌研究[J]. 水利水運(yùn)科學(xué)研究, 1981(3): 89?93.SHA Jinxuan. Piping research on porous medium[J]. Hydro- science and Engineering, 1981(3): 89?93.

[25] 陳亮, 雷文, 張紅宇, 等.非穩(wěn)定流作用下管涌發(fā)生發(fā)展的室內(nèi)試驗(yàn)及理論分析[J]. 巖土工程學(xué)報(bào), 2013, 35(4): 665?662.CHEN Liang, LEI Wen, ZHANG Hongyu, et al. Laboratory simulation and theoretical analysis of piping mechanism under unsteady flows[J]. Chinese Journal of Geotechnical Engineering, 2013, 35(4): 655?662

Numerical analysis of generation and evolution of piping mechanism in meso-level under unsteady flow

NI Xiaodong1, 2, ZHAO Shuailong1, 2, WANG Yuan1, 2

(1. Key Laboratory of Ministry of Education for Geomechanics and Embankment Engineering, Nanjing 210098, China; 2. Institute of Tunnel and Underground Engineering, HoHai University, Nanjing 210098, China)

In order to show the meso-mechanical of piping development under unsteady flow, the dynamic process of piping was studied by sand through model test and numerical simulation of PFC in meso-level. Six models were built to study the variation of water level on particle migration process and the piping failure caused by movement of soil particles. Based on the results of the numerical simulation, the migration and movement of soil and the changes of flowing fields were studied and compared with the results of sand through model test. The results show that one-time loading on the orifice area has the greatest impact, and is the most insecure. The results are identical with those from the engineering experiences. The simulation results agree well with the model test results.

unsteady flow; piping; distinct element method; damage water head; permeability coefficient

10.11817/j.issn.1672-7207.2016.09.033

TV139.1

A

1672?7207(2016)09?3154?08

2015?06?12；

2015?09?18

國(guó)家自然科學(xué)基金資助項(xiàng)目(51309086)；中央高?；究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)資金資助項(xiàng)目(2014B04914)；教育部博士點(diǎn)新教師基金資助項(xiàng)目(20110094120002) (Project(51309086)supported by the National Natural Science Foundation of China; Project(2014B04914) supported by the Fundamental Research Funds for the Central Universities; Project(20110094120002) supported by the Doctoral Fund of Youth Scholars of Ministry of Education of China)

倪小東，博士，副研究員，從事滲透變形方面的研究；E-mail: lulingnxd@126.com

(編輯 趙俊)