何玲麗，田東方，王世梅，陳 勇

(1.水利與環(huán)境國(guó)家級(jí)實(shí)驗(yàn)教學(xué)示范中心，湖北 宜昌 443002；2.三峽大學(xué)水利與環(huán)境學(xué)院，湖北 宜昌 443002；3.三峽庫(kù)區(qū)地質(zhì)災(zāi)害省部共建教育部重點(diǎn)實(shí)驗(yàn)室，湖北 宜昌 443002；4.三峽大學(xué)土木與建筑學(xué)院，湖北 宜昌 443002)

1 研究背景

土石混合體(SRM)是指第四紀(jì)以來(lái)形成的，由具有一定工程尺度、強(qiáng)度較高的巖塊、細(xì)粒土體及孔隙構(gòu)成且具有一定含石量的極端不均勻松散巖土介質(zhì)系統(tǒng)[1]。這類介質(zhì)在我國(guó)各地廣泛分布于邊坡、地基及地下工程中，還經(jīng)常作為工程材料應(yīng)用于大壩、岸堤、路基、機(jī)場(chǎng)等。以SRM的細(xì)觀結(jié)構(gòu)(即其中塊石的形狀、尺寸、分布等信息)為基礎(chǔ)，通過(guò)滲流數(shù)值模擬可以確定水分在土體與塊石中的運(yùn)移過(guò)程。以此為基礎(chǔ)，可開展?jié)B透數(shù)值實(shí)驗(yàn)，揭示細(xì)觀結(jié)構(gòu)對(duì)SRM滲透特征的影響，也可從細(xì)觀尺度研究土石混合體內(nèi)部變形破壞過(guò)程，有助于探索土石混合體破壞機(jī)制。

目前，基于細(xì)觀結(jié)構(gòu)的滲流模擬，多將SRM視為土體和塊石構(gòu)成的二元介質(zhì)，采用達(dá)西定律描述土體和塊石的滲流，利用有限元實(shí)現(xiàn)數(shù)值模擬，如徐文杰等[2]和Yan等[3]的研究。數(shù)值實(shí)驗(yàn)表明，SRM滲透系數(shù)受含石率影響較大，且隨含石率的增大而單調(diào)遞減；同時(shí)，由于土體和塊石的滲流符合達(dá)西定律，故SRM的滲流也符合達(dá)西定律。

然而，實(shí)驗(yàn)表明SRM的滲透系數(shù)隨含石率的增大并非單調(diào)遞減，而是先降后增[4-5]，如圖1a所示；此外，SRM的滲透速度是水力梯度的冪函數(shù)[6-7]，如圖1b所示。從圖1a可見，現(xiàn)有的數(shù)值模擬不能很好地反映SRM的滲透特征。這可能是因?yàn)槲纯紤]土石界面滲透性而造成的。土石界面通常是SRM中的薄弱環(huán)節(jié)和應(yīng)力集中的部位，這將使得土石界面附近土體的孔隙分布與基質(zhì)土體存在一定差異而容易形成集中滲流通道。因此，土石界面的滲透性必須加以考慮。

圖1 土石混合體滲透特性

由于土體和塊石滲流符合達(dá)西定律，那么SRM滲透速度與水力梯度之間的冪函數(shù)關(guān)系，應(yīng)由土石界面的滲流引起。結(jié)合文獻(xiàn)[4，6]的實(shí)驗(yàn)規(guī)律，本文假定界面流量是水力梯度的冪函數(shù)。

綜上，基于SRM細(xì)觀結(jié)構(gòu)進(jìn)行滲流數(shù)值模擬時(shí)，現(xiàn)有方法未考慮土石界面的滲透性，不能正確反映SRM的滲透特性。為此，本文忽略巖塊的滲透性，將SRM視為由土體和土石界面構(gòu)成的二元介質(zhì)，假定土體滲流為二維達(dá)西流，土石界面為一維滲流且流量是水力梯度的冪函數(shù)，基于有限元法實(shí)現(xiàn)考慮土石界面滲透性的SRM滲流數(shù)值模擬。

2 模型的構(gòu)建

2.1 控制方程

SRM滲流模擬區(qū)域如圖2所示，考慮到巖塊的滲透性遠(yuǎn)小于土體，故忽略巖塊的滲透性；則滲流求解域?yàn)橥馏w區(qū)域Ω，以及M條土石界面Γr，r=1，2，…，M。假定土體各向同性且為二維達(dá)西滲流；土石界面為一維滲流。

圖2 SRM滲流模擬區(qū)域示意

土體滲流方程為

(1)

式中，H為總水頭；ks為土體飽和滲透系數(shù)；x，y為坐標(biāo)，位于水平面內(nèi)。

在定水頭邊界ΓH上

(2)

沿任意一條土石界面，建立弧坐標(biāo)l；設(shè)在第r條界面Γr上，從土體到界面的流量為qi(l)，則

(3)

式中，l為沿界面的弧長(zhǎng)坐標(biāo)；nx、ny分別為Γr外法線單位向量在x、y軸的分量(土體為內(nèi)部區(qū)域)。

記q(l)為土石界面Γr的沿程流量，假定只考慮從土體流入界面的水量，根據(jù)水量守恒有

?q/?l=qi

(4)

本文假定q為水力梯度J的冪函數(shù)，則

q=aJb

(5)

式中，a、b為土石界面滲透參數(shù)。

記h(l)為土石界面Γr上的沿程水頭，則J=-?h/?l，令kq=aJb-1，結(jié)合式(4)、(5)有

(6)

2.2 方程的Garlerkin弱解形式

(7)

由格林公式可知

(8)

將式(8)代入式(7)，引入邊界條件(3)得

(9)

式(6)的加權(quán)余量格式為

(10)

根據(jù)格林公式，式(10)左端變?yōu)?/p>

(11)

式(11)右端第一項(xiàng)，當(dāng)Γr自行封閉時(shí)(如圖2中的Γ2)，其值為0；當(dāng)Γr不封閉時(shí)(如圖2中的Γ1)，其值根據(jù)邊界條件式(2)確定。這一條件在有限元求解時(shí)，通過(guò)引入定水頭邊界得以滿足。故該項(xiàng)略去不寫。

將式(10)引入式(11)后，可改寫為

(12)

2.3 模型的有限元求解

用有限元求解式(9)和式(12)時(shí)，可借鑒巖體滲流問(wèn)題中的兩域滲流模型，即分別求解兩場(chǎng)(H和h)、迭代確定qi。但土石界面數(shù)量較多，迭代將帶來(lái)較大的計(jì)算開銷以及收斂問(wèn)題。為此，本文假定水頭在土體與界面接觸處是連續(xù)的，即在任意界面Γr上，h=H。將式(9)與式(12)左右對(duì)應(yīng)相加消去qi，可得

(13)

如圖3所示，本文采用適應(yīng)性強(qiáng)的三角形單元離散Ω；土石界面將被離散為兩節(jié)點(diǎn)線單元，節(jié)點(diǎn)整體編號(hào)與三角形單元共用。

圖3 土體和界面單元示意

三角形單元e內(nèi)的He、線單元τ上的Hτ采用單元節(jié)點(diǎn)水頭Hi、Hm的線性插值，即

(14)

([D]+[R]){H}={f}

(15)

本文采用增量法[8]求解式(15)。

3 模型驗(yàn)證

3.2 算例一

算例一與徐文杰等[2]的滲透數(shù)值實(shí)驗(yàn)進(jìn)行對(duì)比。試樣細(xì)觀結(jié)構(gòu)見文獻(xiàn)[2]中圖6；本文計(jì)算時(shí)，不考慮塊石滲透性，土石界面滲透參數(shù)取a=b=0，其他條件同文獻(xiàn)[2]。對(duì)比結(jié)果見表1，其中，ku為試樣等效滲透系數(shù)與土體滲透系數(shù)之比。

表1 不同含石率相對(duì)滲透系數(shù)ku的對(duì)比

由表1可知，本文與文獻(xiàn)[2]結(jié)果相差均小于2%，表明本文方法和程序可以模擬不考慮界面滲透性的情形。由于塊石透水性相比土體而言通常小幾個(gè)數(shù)量級(jí)，而塊石數(shù)量將隨含石率增大而增多從而降低試樣滲透性。因此，當(dāng)只考慮土體和塊石的滲透性時(shí)，SRM的滲透性必然隨塊石的增多而減小。

3.3 算例二

本算例與王宇等[4,6]的SRM滲透性實(shí)驗(yàn)進(jìn)行對(duì)比。實(shí)驗(yàn)試樣尺寸、邊界條件如圖4所示。試樣中塊石粒徑2～5 mm。文獻(xiàn)[6]中試樣的含石率分別為20%、30%、40%、50%、60%、70%；由于未提供試樣細(xì)觀結(jié)構(gòu)，本文利用文獻(xiàn)[9]方法生成相應(yīng)試樣的細(xì)觀結(jié)構(gòu)；每個(gè)試樣中塊石粒徑分別為2、3、4、5 mm，各粒徑塊石的體積分?jǐn)?shù)均為0.25。限于篇幅，圖4只展示了含石率為20%和70%時(shí)試樣的細(xì)觀結(jié)構(gòu)。

圖4 算例二的試樣尺寸、邊界及細(xì)觀結(jié)構(gòu)(灰色區(qū)域代表巖塊)

圖5 平均流速-水力梯度關(guān)系對(duì)比

數(shù)值模擬時(shí)，土體ks=6.17×10-7m·s-1，土石界面滲透參數(shù)a、b通過(guò)試算確定，具體取值見表2。

表2 算例二計(jì)算參數(shù)及實(shí)驗(yàn)結(jié)果的偏差

數(shù)值模擬所得試樣中的平均流速與文獻(xiàn)[6]的對(duì)比見圖5。本文定義err描述某一含石率時(shí)的n組模擬和實(shí)驗(yàn)數(shù)據(jù)之間的偏差，按式(16)計(jì)算，即

(16)

從圖5可以看出，本文方法可再現(xiàn)SRM的非達(dá)西流特性；從表2所列誤差err來(lái)看，含石率較高(如60%、70%)時(shí)，模擬結(jié)果與實(shí)驗(yàn)結(jié)果較接近，含石率為30%～50%時(shí)偏差較大。

文獻(xiàn)[6]給出了水力梯度等于50時(shí)，滲透系數(shù)隨含石率的變化規(guī)律。數(shù)值模擬所得滲透系數(shù)與文獻(xiàn)[6]對(duì)比見表3。

由表3可知，盡管有時(shí)數(shù)值模擬和實(shí)驗(yàn)結(jié)果相差較大(最大約35%)，但本文方法能再現(xiàn)SRM試樣的滲透性隨含石率的增大而先降后增的現(xiàn)象。

圖6 試樣平均流速與水力梯度的關(guān)系

表3 滲透系數(shù)結(jié)果對(duì)比

4 土石界面滲透參數(shù)影響分析

圖6為a、b取不同數(shù)值時(shí)，試樣平均流速與水力梯度的關(guān)系。由圖6可知，b值越大，非線性越強(qiáng)；a值越大，滲透速度越大；含石率越高，滲透特性對(duì)參數(shù)變化越敏感，即隨a值的增加滲透速度增大得更快，隨b值的增加非線性更強(qiáng)。這是因?yàn)楹试酱?，土石界面?shù)量越多，對(duì)試樣滲透特性影響越大。

本節(jié)初步研究了土石界面滲透參數(shù)(a、b)對(duì)SRM滲透性的影響。其他因素的影響將在后續(xù)研究中進(jìn)一步考慮，例如土石界面流量與水力梯度的其他函數(shù)形式、塊石形狀及分布的隨機(jī)性等。

5 結(jié) 論

(1)將SRM視為由土體和土石界面構(gòu)成的二元介質(zhì)，假定土體滲流為二維達(dá)西流，土石界面為一維滲流且流量是水力梯度的冪函數(shù)，水頭在土體與土石界面處連續(xù)；從而避免采用迭代方法確定土體和界面間的流量交換，構(gòu)建了土體和土石界面滲流場(chǎng)同步求解有限元模型，實(shí)現(xiàn)了考慮土石界面滲透性的土石混合體滲流模擬。

(2)算例驗(yàn)證表明，所建方法可模擬不考慮或考慮土石界面滲透性的情形。當(dāng)考慮土石界面滲透性時(shí)，本文方法可較好地再現(xiàn)土石混合體滲透性隨含石率增大先降后增以及滲流速度與水力梯度的非線性關(guān)系等實(shí)驗(yàn)現(xiàn)象。

(3)滲透數(shù)值實(shí)驗(yàn)表明，土石界面滲透參數(shù)對(duì)SRM滲透特性影響為：a值越大，滲透速度越大；b值越大，非線性越強(qiáng)；含石率越高，SRM滲透特性對(duì)參數(shù)的變化越敏感。