新疆EDG水庫瀝青混凝土心墻壩三維有限元分析

楊春雁，吳俊杰，2

（1.新疆水利水電勘測設計研究院，烏魯木齊 830000；2.南京水利科學研究院，南京 210000）

新疆EDG水庫瀝青混凝土心墻壩三維有限元分析

楊春雁1，吳俊杰1，2

（1.新疆水利水電勘測設計研究院，烏魯木齊 830000；2.南京水利科學研究院，南京 210000）

采用鄧肯-張E～μ模型對EDG水庫大壩進行了三維非線性有限元靜力分析，通過對瀝青混凝土心墻壩的應力與應變分析，研究滿蓄期砂礫料壩體的應力、位移分布。結果表明：大壩變形量值和分布規律基本合理，應力應變均在合理范圍，大壩安全穩定。

瀝青混凝土心墻壩；非線性有限元；靜力分析；應力應變

1 工程簡介

EDG水庫具有防洪、供水、灌溉等綜合利用效益。控制灌溉面積10.7萬畝，擔負下游工業園區、礦區和農業灌溉的供水任務，總庫容4 450萬m3，興利庫容3 050萬m3，死庫容850萬m3，為Ⅲ等中型水利工程。水庫樞紐布置：河道中間是瀝青混凝土心墻壩、右岸岸邊為開敞式正槽溢洪道、以及導流兼沖沙放空洞、左岸為輸水隧洞。

2 計算模型選用

瀝青混凝土心墻和砂礫石堆石料在大壩中起到主要受力作用，為了能夠更好地評測大壩的安全合理性，能夠較好地模擬出瀝青混凝土料與堆石料的應力應變關系，對大壩的安全影響評價至關重要[1]。通過大量的試驗，國內外的研究學者發現瀝青混凝土料與堆石料料的應力應變關系，與剪應力水平而變化和加載應力路徑密切相關。結果得到瀝青混凝土料與堆石料的應力應變關系明顯呈現出非線性特性。在國內外，通過大、小三軸試驗儀，測出這種非線性材料的剪應力與軸應變的關系，用這個非線性的應力應變關系反映瀝青混凝土料與堆石料的本構關系。這類非線性材料的本構關系用的較多的是鄧肯-張模型，還有一些是用的南水模型[2]。通過大量的試驗研究，發現南水模型比鄧肯-張模型計算的變形數值要小一些。為了安全起見，本次數值模擬用鄧肯-張的E～μ模型來反應瀝青混凝土料與堆石料的本構關系，分別用切線變形模量Et和切線泊松比μt為計算參數[3]。

上述公式中共有5個常數分別是K、n、φ、c、Rf，通過這5個常數可以很好地反應出材料的本構關系。

根據大量三軸試驗，鄧肯等人發現應變ε1與側向應變ε1之間也存在一個雙曲線關系，之后提出了切線泊松比的數學表達式：[1]

μt是應力的函數，不是泊松比，也不是常量，它是通過3個參數G、F、D來確定的值[4]。

上述兩式（1）和（2）就是鄧肯-張E～μ模型，兩式中一共有8個參數，通過常規的三軸常規試驗來獲取這些參數[5]。這些參數在推導過程中都具有明確的物理意義，所以可以很好的被工程技術人員所理解，經過長時間的使用，已經積累了較為豐富的運用經驗，基本能夠反映出非線性料的本構關系[6]。

3 壩體幾何模型

對大壩進行剖分，本次計算采用六面體實體單元來剖分壩體以及基巖，通過剖分得到整個模型結構的總結點數為11 655個，總單元數為9 888個。瀝青混凝土心墻節點486個，單元192個，心墻幾何尺寸一般為4.7 m，僅在壩頂最上一層個別區域的實體邊界可達12 m。大壩整體剖分網格如圖1所示。

圖1 有限元剖分網格

4 計算成果與分析

圖2為滿蓄期大壩整體變形示意圖（放大200倍），壩體位移分布與國內外同類工程相似，基本符合瀝青混凝土堆石壩體的變形規律。

圖3～6為滿蓄期壩體最大剖面堆石體的上下游方向（順河向）水平位移、豎向位移及應力等值線分布圖。當水庫蓄水后，瀝青混泥土心墻在水荷載的作用下，大壩整體與心墻有向下游移動的趨勢，向上游最大水平位移減小到8.0 cm，位于上游壩體1/3壩高附近；向下游最大水平位移增大為15.7 cm，位于下游壩體2/3壩高附近。最大豎向位移為32.5 cm，大壩最大沉降約為壩高的0.38%，最大值的位置與竣工期相比稍偏向上游，均在已建工程監測值范圍內。

圖7～9為瀝青混凝土心墻的大、小主應力及剪應力水平等值線圖，根據圖7、8可以得到竣工期的大、小主應力均為壓應力，沒有出現主拉應力區。大主應力最大值為1.833 MPa，小主應力最大值為1.100 MPa，蓄水期的大主應力的最大值為1.885 MPa，小主應力最大值為1.238 MPa。在心墻下部1/3的壩高處，個別單元的第三主應力有較小的拉應力出現，其最大極值為-0.045 MPa。這不會因出現小范圍的拉應力而引起拉裂，影響瀝青混凝土心墻的防滲性能。圖9為瀝青心墻的的剪應力水平等值線圖，從圖中可以得到心墻最大值剪應力值出現在左岸與心墻接觸面附近為0.36，但均小于1。

圖2 滿蓄期大壩整體變形圖（放大200倍）

圖3 滿蓄期最大剖面水平位移等值線圖（cm）

圖4 滿蓄期最大剖面垂直位移等值線圖（cm）

圖5 滿蓄期最大剖面大主應力等值線圖（MPa）

圖6 滿蓄期最大剖面小主應力等值線圖（MPa）

圖7 滿蓄期瀝青混凝土心墻軸線剖面大主應力等值線圖（MPa）

5 結論

根據室內大、小三軸試驗資料，整理出鄧肯-張模型的力學參數，對EDG瀝青混凝土心墻堆石壩進行三維非線性有限元靜力計算，結果表明：靜力計算得到的大壩變形量值和分布規律基本合理，最大沉降發生在蓄水期為為壩高的0.38%。壩體位移分布在國內外工程的經驗范圍內。瀝青混凝土心墻的大小主應力均為壓應力，有很少的單元出現拉應力區域，但是拉應力極值都小于瀝青材料的拉應力值，可知瀝青心墻在竣工期與蓄水期時不會出現拉伸破壞。同時瀝青心墻的剪應力水平值均小于1，基本沒有單元在竣工期與蓄水期時產生剪切破壞，這樣也降低了瀝青心墻水力劈裂的可能性，研究成果基本可以作為大壩設計的依據。

圖8 滿蓄期瀝青混凝土心墻軸線剖面小主應力等值線圖（MPa）

圖9 滿蓄期瀝青混凝土心墻軸線剖面應力水平等值線圖

[1] 謝定義，姚仰平，黨發寧.高等土力學[M].北京：高等教育出版社，2008.

[2] 顧淦臣，黃金明.混凝土面板堆石壩的堆石本構模型與應力變形分析[J].水力發電學報，1991（1）：12-14.

[3] 吳俊杰，鳳 煒，劉 亮，等.山前傾斜平原水庫全庫盤防滲復合土工膜應力應變分析[J].水電能源科學，2013（12）：22.

[4] 吳 強.應力應變有限元法在面板堆石壩設計中的應用[J].廣西水利水電，2009（1）：23-27.

[5] 唐 崗.斧子口水利樞紐溢流壩三維有限元分析[J].廣西水利水電，2013（3）：28-32.

[6] 紀 亮，張建海，何昌榮，等.雙江口水電站心墻堆石壩應力應變數值分析[J].廣西水利水電，2008（4）：4-7.

（責任編輯：周 群）

3D finite element analysis for asphalt concrete core wall dam of EDG Reservoir in Xinjiang

YANG Qiu-yan1,WU Jun-jie1,2
（1.Xinjiang Water Conservancy and Hydropower Investigation and Design Institute,Urumqi 830000,China;2.Nanjing Hydraulic Research Institute,Liaoning,Nanjing 210000,China）

3D finite element static analysis was conducted with the Duncan-Chang model for the main dam of EDG Reservoir.Based on stress and strain analysis of the asphalt concrete core wall dam,the stress and displacement dis?tribution of sand gravel dam body under full impoundment was studied.The results of study demonstrate that main dam is safe and stable with deformation value and distribution rule being rational,stress and strain within accept?able range.

Asphalt concrete core wall dam;non-linear finite element;static analysis;stress and strain

TV641.41

B

1003-1510（2016）02-0051-03

2016-01-06

楊春雁（1980-），女，新疆烏魯木齊人，新疆水利水電勘測設計研究院工程師，學士，主要從事水利工程勘察、設計工作。