董國(guó)慶 何亮

摘要：基于大型靜力三軸壓縮試驗(yàn)、等向壓縮和卸載試驗(yàn)，構(gòu)造了簡(jiǎn)捷的塑性模量表達(dá)式，建立了適用于堆石料的廣義塑性本構(gòu)模型。堆石料的本構(gòu)模型模擬和三軸壓縮試驗(yàn)的對(duì)比結(jié)果表明：廣義塑性模型不僅能反映堆石料在低圍壓下的剪脹性和高圍壓下的剪縮性，而且能反映堆石料因顆粒破碎而引起峰值應(yīng)力比和剪脹應(yīng)力比的非線性變化特性。此外，有限元靜力計(jì)算結(jié)果表明，壩體水平位移和沉降特性均符合高面板堆石壩的一般變形規(guī)律，量值在合理范圍內(nèi)。

關(guān)鍵詞：堆石料;廣義塑性模型;塑性模量;有限元分析

中圖分類(lèi)號(hào)：TV641.4;TU47

文獻(xiàn)標(biāo)志碼：A

doi：10.3969/j.issn.1000- 1379.2019.05 .024

堆石料在壩體內(nèi)高應(yīng)力水平作用下會(huì)產(chǎn)生嚴(yán)重的顆粒破碎，影響其應(yīng)力應(yīng)變特性。因此，高面板堆石壩堆石料本構(gòu)模型的建立必須充分考慮堆石料的顆粒破碎特性。

廣義塑性模型（簡(jiǎn)稱(chēng)P-Z模型）直接定義塑性流動(dòng)方向、加載方向和塑性模量，用一套參數(shù)即可完成土石壩的靜動(dòng)力分析過(guò)程，相比于經(jīng)典彈塑性模型，具有更大的靈活性[1]。該模型雖然在黏土和砂土中得到了廣泛的應(yīng)用，但不能很好地模擬堆石料因顆粒破碎而表現(xiàn)出的應(yīng)力應(yīng)變特性。陳生水等[2]采用峰值摩擦角和特征點(diǎn)摩擦角來(lái)反映堆石料因顆粒破碎而表現(xiàn)出的峰值應(yīng)力比與剪脹應(yīng)力比的非線性變化規(guī)律，構(gòu)建了一個(gè)固相硬度參數(shù)作為堆石料抵抗破碎能力的一種綜合度量。王占軍等[3]建立了能夠反映堆石料顆粒破碎的剪脹方程，構(gòu)造了隨平均應(yīng)力、剪應(yīng)力比和密實(shí)度變化的塑性模量。此外，還有許多學(xué)者探索了考慮堆石料顆粒破碎的廣義塑性模型表達(dá)式[4-10]。目前，上述改進(jìn)模型普遍存在塑性模量表達(dá)式復(fù)雜、參數(shù)多、高圍壓下應(yīng)力應(yīng)變關(guān)系適應(yīng)性較差等問(wèn)題。

筆者基于堆石料的大型靜力三軸試驗(yàn)、等向壓縮和卸載試驗(yàn)結(jié)果，研究考慮顆粒破碎的塑性模量簡(jiǎn)捷表達(dá)式，建立了適用于堆石料的廣義塑性本構(gòu)模型，并開(kāi)展了某高面板堆石壩的有限元分析。

1 堆石料的廣義塑性本構(gòu)模型

高圍壓下顆粒破碎現(xiàn)象對(duì)堆石料的壓縮特性有顯著影響。當(dāng)平均正應(yīng)力p超過(guò)屈服應(yīng)力后，堆石料顆粒開(kāi)始破碎，堆石料的等向壓縮曲線顯著下彎，塑性體積變形不斷累積。為了避免高圍壓下可能出現(xiàn)孔隙比e為負(fù)的不合理現(xiàn)象，采用指數(shù)函數(shù)描述堆石料的壓縮特性[2]：

大量的堆石料三軸壓縮試驗(yàn)資料表明：由于顆粒破碎，因此堆石料的破壞線和剪脹線均呈現(xiàn)出明顯的非線性特性，都隨平均正應(yīng)力p的增大逐漸下彎，且剪脹線的非線性程度低于破壞線的。因此，引入峰值和剪脹摩擦角，構(gòu)建堆石料的峰值和剪脹應(yīng)力比與平均應(yīng)力的關(guān)系公式[2]：

2 高面板堆石壩有限元分析

2.1 計(jì)算模型和計(jì)算參數(shù)

以某高面板堆石壩為例進(jìn)行靜力有限元分析，該壩高256 m，壩頂寬12 m，上游壩坡坡比為1：1.6，下游壩坡坡比為1：1.4或1：1.3。有限元計(jì)算網(wǎng)格見(jiàn)圖1。

由靜力大三軸壓縮試驗(yàn)、等向壓縮和卸載試驗(yàn)結(jié)果獲得的計(jì)算參數(shù)見(jiàn)表1，其中參數(shù)h_s、n和K通過(guò)等向壓縮和卸載試驗(yàn)確定，其他參數(shù)通過(guò)三軸壓縮試驗(yàn)確定。

對(duì)比主堆石料的本構(gòu)模型模擬曲線與三軸壓縮試驗(yàn)曲線（見(jiàn)圖2.σ1 -σ3為偏應(yīng)力，εa為軸向應(yīng)變，εv為體積應(yīng)變），可得出以下結(jié)論：試驗(yàn)曲線和模型模擬曲線吻合較好：廣義塑性模型不僅可以較好地模擬堆石料在低圍壓下的剪脹性和高圍壓下的剪縮性，而且能反映堆石料因顆粒破碎而引起峰值應(yīng)力比和剪脹應(yīng)力比的非線性特性。

2.2 有限元分析結(jié)果

圖3為竣工期大壩最大剖面的順河向水平位移和沉降分布圖。竣工期壩體的水平位移大致以壩軸線為界分為上游水平位移和下游水平位移，各自指向坡外方向。下游壩體的水平位移比上游的大，上、下游水平位移最大值分別為12 cm和38 cm;竣工期壩體的沉降呈現(xiàn)出規(guī)律性環(huán)狀分布，從壩體內(nèi)部向四周遞減，最大沉降量為114 cm，發(fā)生在約1/2壩高處。

圖4為蓄水期大壩最大剖面的沉降和水平位移分布圖。與竣工期相比，在水壓力作用下，壩軸線上游側(cè)指向上游的水平位移被壓回，大壩的上游水平位移減小，下游水平位移增大，沉降增大。上、下游水平位移最大值分別為4 cm和43 cm，沉降量為121 cm。與竣工期相比，上游水平位移最大值減小了8 cm，下游水平位移最大值增大了5 cm，沉降量增大了7 cm。

3 結(jié)論

（1）相比于其他文獻(xiàn)中的塑性模量表達(dá)式，本文構(gòu)造的塑形模量表達(dá)式更為簡(jiǎn)捷，且能較好地反映顆粒破碎堆石料的塑性變形規(guī)律。

（2）本文建立的廣義塑性本構(gòu)模型可以較好地模擬堆石料在低圍壓下表現(xiàn)出的剪脹性和高圍壓下的剪縮性。

（3）有限元計(jì)算結(jié)果顯示，壩體水平位移和沉降均符合高面板堆石壩的一般變形規(guī)律，量值在合理范圍內(nèi)，表明該模型應(yīng)用于壩體的靜力有限元計(jì)算是可行的。

參考文獻(xiàn)：

[1] ZIENKIEWICZ O C，CHAN A H C，PASTOR M， et al.Computational Geomechanics with Special Reference toEarthquake Engineering[M].New York： John Wiley& Sons，1999：55-71.

[2] 陳生水，傅中志，韓華強(qiáng)，一個(gè)考慮顆粒破碎的堆石料彈塑性本構(gòu)模型[J].巖土工程學(xué)報(bào)，2011，33（10）：1489-1495.

[3] 王占軍，陳生水，傅中志，堆石料剪脹特性與廣義塑性模型[J].巖土力學(xué)，2015，36（7）：1931-1938.

[4]朱晟，魏匡民，林道通，筑壩土石料的統(tǒng)一廣義塑性本構(gòu)模型[J].巖土工程學(xué)報(bào)，2014，36（8）：1394-1399.