摘 要：隨著計算機(jī)技術(shù)的發(fā)展和廣泛應(yīng)用， 各種數(shù)值分析方法開始應(yīng)用到巖土工程中來， 其中有限單元法因為其對邊界條件的適應(yīng)性以及能很好的模擬各種土的本構(gòu)關(guān)系， 使其在巖土工程計算中得到廣泛的應(yīng)用。本文選用基于 ABAQUS計算軟件環(huán)境下的強(qiáng)度折減有限元法，進(jìn)行土石壩在設(shè)計洪水位作用下的應(yīng)力、位移及穩(wěn)定性分析。

關(guān)鍵詞：土石壩；ABAQUS；強(qiáng)度折減法；穩(wěn)定性

DOI：10.16640/j.cnki.37-1222/t.2017.14.237

0 引言

有限元進(jìn)行邊坡穩(wěn)定性分析受到關(guān)注。與傳統(tǒng)方法相比，有限元法優(yōu)勢明顯。本文選用基于 ABAQUS 計算軟件環(huán)境下的強(qiáng)度折減法，進(jìn)行土石壩在設(shè)計洪水位作用下的應(yīng)力、位移及穩(wěn)定性分析[1]。

1 算例分析

均質(zhì)土石壩，壩頂高程50.0m，底高程0.0m，頂寬度 8.0m，上游壩坡 1：2.5，高程20.0m 處上游壩坡設(shè)有3m 寬馬道，高程 25.0m處下游壩坡設(shè)有 2m 寬馬道，下游壩坡 1：2.0。上游設(shè)計洪水位高程 46.0m，下游水位高程4.0m，下游排水棱體頂寬 3.0，頂高程 5.0m，底高程 0.0m，上游邊坡 1：1.5，下游邊坡1：1.5。

壩體材料如下表所示：

取土石壩部位最大橫斷面為典型剖面， 進(jìn)行二維非線性有限元分析。網(wǎng)格單元類型為 CPE4。左右邊界均采用水平約束，底面采用水平豎向約束。先對整體模型進(jìn)行地應(yīng)力平衡，然后施加靜水壓力分析壩體的應(yīng)力及位移情況，最后分析在設(shè)計洪水作用下的壩體穩(wěn)定性[2]。取河床部位最大橫斷面為計算的典型剖面，進(jìn)行有限元分析。采用的Mohr-Coulomb 模型，摩擦系數(shù)取 0.3。

2 計算結(jié)果應(yīng)力分析

ABAQUS 中以拉為正， ABAQUS 中的小主應(yīng)力對應(yīng)于巖土工程中的大主應(yīng)力。 在設(shè)計洪水作用下土石壩的屈服應(yīng)力、 小主應(yīng)力等值線云圖如下圖 1和圖 2?？梢?，壩體屈服應(yīng)力及小主應(yīng)力有較好的分布規(guī)律，主應(yīng)力值隨壩體斷面高程增加而減小，壩體小主應(yīng)力的最大值為 0.8545Mpa。

圖 3、 圖 4分別為土石壩在自重作用下的水平位移云圖和靜水壓力作用下的水平位移云圖。 位移符號規(guī)定：水平位移向下游為正， 單位為 m。 由圖可知， 大壩在自重作用下， 壩體的水平位移呈大致對稱分布，各自指向坡外方向； 在施加靜水壓力后壩體上游坡腳水平位移發(fā)生明顯向下游傾斜情況， 最大位移為1.182m， 這是由于壩基密度較低， 在設(shè)計洪水位作用下容易發(fā)生大的位移。

3 結(jié)語

依托算例分析，考慮洪水時的約束關(guān)系，建立了相應(yīng)的結(jié)構(gòu)圖仿真模型，并計算模擬典型洪水過程的水庫水位變化曲線。結(jié)果表明：由于洪水歷時短，上、下游壩坡的穩(wěn)定系數(shù)變化趨勢與水庫水位基本一致。本文方法為洪水作用下的土石壩應(yīng)力、位移及穩(wěn)定性分析提供了新的途徑。

參考文獻(xiàn)：

[1]祁慶和.水工建筑物[M].北京：水利水電出版社，1992.

[2]費(fèi)康，張建偉.ABAQUS在巖土工程中的應(yīng)用[M].北京：中國水利水電出版社，2010.

作者簡介：劉強(qiáng)（1991-），男，湖北黃梅縣人，研究生，研究方向：抗拔樁基礎(chǔ)。