唐林斌(廣西恒晟建設(shè)集團(tuán)有限公司，廣西桂林 541199)

無(wú)限深透水壩基滲流場(chǎng)與應(yīng)力場(chǎng)耦合研究

唐林斌
(廣西恒晟建設(shè)集團(tuán)有限公司，廣西桂林 541199)

在建造土石壩的工程當(dāng)中，無(wú)限深透水壩的基滲流問(wèn)題顯得非常的重要，它關(guān)乎到土石壩工程土質(zhì)量的好壞。而經(jīng)過(guò)相關(guān)專家的研究后發(fā)現(xiàn)，采用滲流場(chǎng)與應(yīng)力場(chǎng)耦合的方式能夠很好的解決這個(gè)問(wèn)題，那么下面我們就從對(duì)兩場(chǎng)相互影響的分析，以及耦合計(jì)算的方法這兩方面的內(nèi)容，來(lái)具體的討論一下有關(guān)無(wú)限深透水壩基滲流場(chǎng)與應(yīng)力場(chǎng)耦合方面的問(wèn)題。

無(wú)限深透水壩基；滲流場(chǎng)；應(yīng)力場(chǎng)

我國(guó)有些地區(qū)比較容易見(jiàn)到在無(wú)限滲透水地上建設(shè)土石壩?？紤]到壩基的透水性能較強(qiáng)，因此非常有必要對(duì)其進(jìn)行應(yīng)力分析以及滲流計(jì)算。不過(guò)目前所采用的計(jì)算方法是把應(yīng)力場(chǎng)和滲流場(chǎng)單獨(dú)進(jìn)行計(jì)算，如果不關(guān)心應(yīng)力場(chǎng)和滲流場(chǎng)的耦合關(guān)系，造成所計(jì)算出來(lái)的結(jié)果往往也和現(xiàn)實(shí)狀況有很大的出入。而為了解決這一問(wèn)題，相關(guān)專家經(jīng)過(guò)了認(rèn)真的研究，終于找到了將兩者進(jìn)行耦合的方法。

1.兩場(chǎng)相互影響分析

1.1 滲流場(chǎng)對(duì)應(yīng)力場(chǎng)的影響

在無(wú)限深透水地基對(duì)應(yīng)力進(jìn)行計(jì)算，一定要充分的考量它的滲透力。例如在其中的一個(gè)時(shí)段上已經(jīng)獲得了該地基的水頭函數(shù)，那么就能夠計(jì)算出滲流部位的滲水力度。如果用u來(lái)代表滲水力度的話，那么滲流場(chǎng)和應(yīng)力場(chǎng)的關(guān)系則可以用公式u=x(i-z)來(lái)表示。通過(guò)這樣的公式，還能夠計(jì)算出滲透體積力的強(qiáng)度以及偏積分，而在利用有限元算法的情況下，則要把u和f轉(zhuǎn)變成單一的結(jié)點(diǎn)來(lái)計(jì)算。

1.2 應(yīng)力場(chǎng)對(duì)滲透場(chǎng)的影響

無(wú)限深透水地基相對(duì)比較軟弱，如果在外載荷的沖擊下，很有可能對(duì)出現(xiàn)變形的情況，這樣土體孔隙率也會(huì)隨之改變，從而對(duì)壩基的滲水性造成一定的影響。如果只是應(yīng)對(duì)空隙的壓縮能力，那么可以用下面的關(guān)系式來(lái)進(jìn)行表達(dá)：

dn=(1-n)adu。其中n代表的是介質(zhì)孔隙率，u代表的是孔隙水壓力,a則代表的是壓縮系數(shù)。同時(shí)，由于n的函數(shù)是k，而K反應(yīng)的是土體滲透系數(shù)，所以K和n就會(huì)形成如下的關(guān)系式：

K=K0=[]3[]2。其中，n0和n分別代表的是壩基土體的初始以及沉降期間的孔隙度，K0和K分別代表的是壩基土體在初始以及被壓縮以后的滲透系數(shù)。而孔隙率以及滲透系數(shù)之間則是互相依賴的關(guān)系。

通常要采用geostudio軟件來(lái)對(duì)壩基的滲水性進(jìn)行計(jì)算，不過(guò)與此同時(shí)，還要做好壩基土體的約束工作，一般情況下，滲透強(qiáng)度是由土體的變形程度來(lái)決定其大小。如果土體的厚度用Mt，沉降量用St來(lái)表示的的話，則會(huì)出現(xiàn)下面的關(guān)系式：n=n0-。而通過(guò)這組關(guān)系式則能夠看出，土層沉降程度主要是由沉降系數(shù)所決定的，因此這就把土層沉降和滲透系數(shù)很好的結(jié)合到了一起，從而得以完成兩場(chǎng)的耦合工作。

2.數(shù)學(xué)模型和計(jì)算方法

2.1 控制方程

如果水體具有大量的空隙，那么在滲流飽和區(qū)就能夠創(chuàng)建相應(yīng)的方程式，而且方程式的要以增量形式為主：{△o}=[D]{ △e}+{m}△u。其中[D]代表四十線彈性矩陣，{m}代表的是等等參向量，而△u則表示孔隙水的壓力。

利用虛功原理，可以把外力轉(zhuǎn)變成內(nèi)力。而內(nèi)力，也就是我們經(jīng)常說(shuō)的變性能，并要把它保存在彈性體里面。加入外力用{F}來(lái)表示的話，那么代表虛功的方程式則為：{e*}{△o}dv={o}{F}dV。

而如果與{△o}=[D]{△e}+{m}△u進(jìn)行結(jié)合的話，那么就可以形成如下的關(guān)系式：[B]T[D][B]{△o}+[B]T{m}[N]{△u}=F。其中[D]代表的是應(yīng)變矩陣，{△o}代表的是位移增量，[N]代表的是形函數(shù)。加入代表剛度系數(shù)的[K]=[B][D][B],而[L]=[B]{m}[N]代表耦合矩陣，那么[D][B]{△o}+[B]T{m}[N]{△u}=F就可以化簡(jiǎn)成[K]{△o}+[L]{△u}={△F}。

而通過(guò)簡(jiǎn)化的關(guān)系式[K]{△o}+[L]{△u}={△F}能夠發(fā)現(xiàn)，滲透場(chǎng)和應(yīng)力場(chǎng)一旦形成耦合關(guān)系，方程中就會(huì)多出{△u}這一項(xiàng)，這種情況的出現(xiàn)則能夠很好的表示出在滲流場(chǎng)里，一旦孔隙水的壓力發(fā)生壓力變化，能夠給應(yīng)力場(chǎng)造成什么程度的影響。

2.2 滲流方程

通過(guò)對(duì)虛功原理的縝密分析，能夠很好的推斷出滲流方程：[dv=u*vndA。其中Vn代表邊界通量。

而如果對(duì)滲流方程[dv=u*vndA進(jìn)行由原先離散處理，那么則會(huì)形成如下的方程式：[B]T[Kw][B]{u}dV+〈N〉vndv。其中[Kw]代表的是滲透矩陣，而在經(jīng)過(guò)對(duì)[B]T[Kw][B]{u}dV+〈N〉vndv的轉(zhuǎn)化以后，會(huì)形成〈N〉{m}[B]dv,這個(gè)則代表滲流耦合矩陣。此后，再對(duì)滲流方程[dv=u*vndA從t到t+△t進(jìn)行時(shí)間上的劃分，則會(huì)形成：同時(shí)，還要將時(shí)間步長(zhǎng)滲透到這個(gè)公示當(dāng)中，然后再對(duì)這套公式實(shí)行時(shí)間差分，從而得到-[Kf]{u}+(1-s)[Kf]{u}+[Lf]{e}=△t〈N〉[vn]+(1-s)vn]dA。

而且為了可以列出表示孔隙水壓力增大的公式，則可以將-[Kf]{u}+(1-s)[Kf]{u}+[Lf]{e}=△t〈N〉[vn]+(1-s)vn]dA。轉(zhuǎn)化成為：[L]{△e}-[Kf]{△u}= △t{Q}t+△t+[K]{u}。其中{Q}代表的是邊界節(jié)點(diǎn)的流量，而[Kw]則代表滲透矩陣。

2.3 耦合計(jì)算的步驟

第一步，是要設(shè)立{Ho}，也就是我們經(jīng)常說(shuō)的滲流場(chǎng)的水頭分布的情況，在分布任務(wù)完成了以后，還要計(jì)算出材料的具體參數(shù)，這樣的話就能夠得到滲透系數(shù)的具體數(shù)值，那么此后再通過(guò)表示孔隙水壓力增大的公式[L]{△e}-[Kf]{△u}=△t{Q}t+△t+[K]{u}，最終算出滲流水頭的具體分布情況。

第二步，根據(jù)上面所算出的滲流水頭的具體分布情況，找出水頭的具體變化情況，并根據(jù)水頭的具體變化,來(lái)最終的算出孔隙水壓力。

第三步，根據(jù)計(jì)算出的體力學(xué)參數(shù)，來(lái)得到應(yīng)力場(chǎng)的剛度矩陣，同時(shí)要在土體的上方施以外荷矩陣，然后再根據(jù)控制方程中的[B]T[D][B]{△o}+[B]T{m}[N]{△u}=F，算出壩基土體應(yīng)力場(chǎng)的具體情況。

第四步，通過(guò)對(duì)壩體加以水平約束的方式，來(lái)充分的了解接壩基土體的位移變形情況，并其還有引入方程式[dv=u*vndA，這樣就能夠得出土體孔隙率的具體數(shù)值，然后再引用方程式[B]T[Kw][B]{u}dV+〈N〉vndv，以此來(lái)對(duì)壩基土體的滲透系數(shù)進(jìn)行最終的確認(rèn)。

第五步，通過(guò)引用[L]{△e}-[Kf]{△u}=△t{Q}t+△t+[K]{u}，能夠算出水頭進(jìn)行第二次分布以后的具體情況，這樣就會(huì)對(duì)孔隙水壓力所出現(xiàn)的變化能夠有一個(gè)詳細(xì)的了解。然后再根據(jù)孔隙水壓力所出現(xiàn)的變化，來(lái)對(duì)應(yīng)力場(chǎng)所出現(xiàn)的變化進(jìn)行系統(tǒng)的調(diào)查，并得出空隙率以及壩基土體的滲透系數(shù)的數(shù)值。

第六步，一直要等到迭代工作完結(jié)以后，并且計(jì)算達(dá)到一定精度的情況下，才可以結(jié)束計(jì)算的工作。

3.結(jié)束

通過(guò)以上的內(nèi)容我們能夠了解到，將無(wú)限深透水壩基中的滲流場(chǎng)與應(yīng)力場(chǎng)進(jìn)行耦合是一項(xiàng)非常復(fù)雜的工作，相關(guān)研究人員必須對(duì)兩場(chǎng)相互之間的影響關(guān)系進(jìn)行系統(tǒng)的分析以及研究出耦合計(jì)算的正確方法，這樣才能夠讓滲流場(chǎng)與應(yīng)力場(chǎng)實(shí)現(xiàn)耦合。所以，相關(guān)專家應(yīng)該要刻苦努力，不斷專研，研究出更加完善的耦合方案，從而讓我國(guó)的土石壩的建造技術(shù)邁向一個(gè)新的高度。

[1]樓??；毛海濤；何華祥；邵東國(guó)；王曉菊.無(wú)限深透水壩基滲流場(chǎng)與應(yīng)力場(chǎng)耦合分析[J].水動(dòng)力學(xué)研究與進(jìn)展A輯, 2015，16(2):61-62.

[2]毛海濤；何華祥；邵東國(guó)；王曉菊.無(wú)限深透水壩基上懸掛式防滲墻控滲試驗(yàn)研究[J].西安理工大學(xué), 2015，16(35):1-1.

[3]毛海濤；王曉菊；付亞男；姜海波[J].無(wú)限元法在無(wú)限深透水壩基滲流計(jì)算中的應(yīng)用[J]水利水運(yùn)工程學(xué)報(bào) ,2014,12(11):51-52.

G322

B

1007-6344（2016）04-0264-01