蘇超+董安雨+劉崇巍+何志亮+徐匯

摘要：碾壓混凝土重力壩結(jié)構(gòu)簡單，施工方便，應(yīng)用廣泛。在實際工程中，經(jīng)常會遇到地基的各向異性問題。應(yīng)用非線性有限元結(jié)構(gòu)分析軟件ABAQUS，以某碾壓混凝土重力壩的岸坡壩塊為研究對象，建立了碾壓重力壩的三維有限元計算模型。在正常蓄水工況下設(shè)置了六組水平向與鉛直向彈性模量比值不同的計算組，通過對比分析，研究其可不考慮各向異性的材料比值區(qū)間。

關(guān)鍵詞：岸坡壩段；三維有限元；橫觀各向同性地基；應(yīng)力應(yīng)變；彈性模量

中圖分類號：TV642 文獻(xiàn)標(biāo)志碼：A 文章編號：1672-1683（2017）04-0155-06

Abstract：RCC gravity dam has a simple structure，can be constructed conveniently，and is widely used nowadays.In practical engineering，foundation anisotropy problems are often encountered.In this paper，we used the finite element analysis software ABAQUS，took the slope dam blocks of a roller compacted concrete gravity dam as the basic computational structure，and established a 3D finite element calculation model of RCC gravity dam.Under the normal storage condition，six different calculation groups were set that have different ratios of horizontal to vertical elastic modulus.Through comparative analysis，we studied the material ratio range within which anisotropy needs not to be considered.

Key words：bank slope dam section；three-dimensional finite element；transverse isotropic foundation；stress-strain；elastic modulus

碾壓混凝土重力壩結(jié)合了混凝土重力壩的高安全性與土石壩的施工高效性，是一種施工方便、安全度高的壩型[1-3]。目前在設(shè)計碾壓混凝土重力壩時使用的方法與設(shè)計常態(tài)混凝土重力壩時的方法差別不大，使用的都是材料力學(xué)法[4-9]。當(dāng)壩基的材料分布均勻時，材料力學(xué)法的計算結(jié)果能夠較好的反映出壩體的實際受力狀況，但當(dāng)壩基地質(zhì)條件復(fù)雜時就會產(chǎn)生一定誤差，通常采用有限元方法進(jìn)行水工設(shè)計計算[10-13]。在有限元計算中，考慮地基的分區(qū)以及地基的塑性變形較多，但考慮地基的各向異性問題的研究相對較少，在實際問題當(dāng)中應(yīng)該考慮這一問題的影響[14-17]。若壩段處于岸坡位置，建基面通常設(shè)計為臺階狀，應(yīng)力分布將更為復(fù)雜，會出現(xiàn)應(yīng)力集中現(xiàn)象。本文以某碾壓混凝土重力壩的岸坡壩塊為研究對象，在有限元計算中考慮了地基的橫觀各向同性性質(zhì)，設(shè)置了六組水平向與鉛直向不同比值彈性模量的計算組，研究材料橫觀各向同性的影響，通過對比分析，探討需要考慮各向異性的參數(shù)區(qū)間，為設(shè)計人員提供參考依據(jù)。

1 橫觀各向同性計算原理

對于橫觀各向同性問題，通過改變各向同性材料的剛度矩陣便可以解決[18-21]。在有限元計算中材料一般是各向同性材料，即在各個方向上的材料參數(shù)是相同的，這種材料在其彈性矩陣中包含彈性模量E、剪切模量G和泊松比[WTB1X]v[WT]3個參數(shù)，其彈性矩陣表達(dá)式如下：

2 計算模型

本文算例為某碾壓混凝土重力壩，有限元計算壩塊位于右邊岸坡上，為適應(yīng)斜坡壩段的側(cè)滑趨勢增加穩(wěn)定性，建基面設(shè)為臺階型，見圖2。壩塊長22 m，最大壩高94 m，最小壩高84 m，壩底最低高程220 m，最高高程230 m，是擋水壩段。壩身主要由RCC15混凝土和RCC20混凝土澆筑碾壓而成，其表面覆蓋有與壩體同等級的防滲混凝土，厚度為1 m。大壩地基包含微風(fēng)化層和弱風(fēng)化層，其中弱風(fēng)化層地基主要位于壩底前部，其橫觀各向同性性質(zhì)較為明顯。

計算模型的坐標(biāo)原點位于建基面上（圖2），坐標(biāo)系的X軸正方向指向下游，Y軸正方向指向左岸，Z軸正方向豎直向上；建模時為了避免地基范圍過小而對壩體應(yīng)力產(chǎn)生影響，所以在壩體上、下游各延伸了兩倍壩高長度，并在建基面向下也延伸兩倍壩高厚度；在邊界約束方面，在地基垂直于X軸的兩個表面約束其X軸方向位移，在垂直于Y軸的兩個表面約束其Y軸方向位移，在底面約束其三個方向上的位移；剖分網(wǎng)格時采用了C3D6單元和C3D8單元，具體剖分結(jié)果見圖3，整個模型離散為12 688個單元和14 642個節(jié)點。

3 橫觀各向同性地基計算參數(shù)

壩基有4種材料分區(qū)，考慮到材料參數(shù)的差異不大，為突出討論研究地基橫觀各向同性的影響，所以在計算模型中將地基設(shè)為勻質(zhì)地基。以地基水平向彈性模量與法向彈性模量的比值不同分為六組進(jìn)行計算（法向減小時水平向增加），最后對比各組的計算結(jié)果。六組中第六組為對照組，各組的地基材料參數(shù)及分組情況見表1。

根據(jù)設(shè)計規(guī)范的規(guī)定，計算中考慮的荷載包括：自重、靜水壓力、揚(yáng)壓力、淤沙壓力和浪壓力。其中混凝土的重度為24 kN/m3，靜水壓力取正常蓄水位時的壓力。大壩在上游設(shè)置了防滲帷幕和排水孔，由于壩體是實體重力壩并且是岸坡壩段，根據(jù)設(shè)計規(guī)范滲透壓力系數(shù)應(yīng)取0.35，計算結(jié)果見表2。淤沙壓力按規(guī)范計算，淤沙容重取為19 kN/m3。

4 結(jié)果對比分析

4.1 X軸方向位移對比

圖4為第一組壩體X軸方向位移云圖，各組壩體在X軸方向上的位移趨勢是大致相同的，隨著高度的增加X軸方向的位移在不斷地增大，且最大水平位移都出現(xiàn)在壩頂?shù)奈恢茫⑶译S著地基法向彈性模量的增大，壩體在X軸方向上的各組最大位移值在不斷地減小。

其最大位移值變化見表3。從表中可以看出，第五組與第六組的最大位移值差0.089 mm，且差值比例為4.34%，小于5%，此時地基橫觀各向同性性質(zhì)對水平方向位移已經(jīng)影響不大了。

4.2 Y軸方向位移對比

壩體在Y軸方向上的位移值較小，特別是壩底位移，其最大位移值幾乎都可以忽略不計，地基橫觀各向同性性質(zhì)對Y軸方向位移影響并不大，所以不以Y軸方向位移大小為評判標(biāo)準(zhǔn)。

4.3 Z軸方向位移對比

圖5為第一組壩體Z軸方向位移云圖，壩體在Z軸方向上的位移趨勢大致相同，且出現(xiàn)最大沉降的位置大致相同。

表4為各組壩體在Z軸方向上的位移值，從表中可以看出隨著地基法向彈性模量的增大，壩體的沉降在快速減小，這說明地基的法向彈性模量對壩體的沉降有很大的影響。對比第五組的數(shù)據(jù)，與對照組的沉降值差為0.042 cm，差值比例為1.89%，此時地基橫觀各向同性性質(zhì)對沉降的影響已經(jīng)很小了。

4.4 拉應(yīng)力結(jié)果對比

各組的拉應(yīng)力分布形式大致相同，并且都在壩體220截面下游拐角處出現(xiàn)了拉應(yīng)力集中區(qū)域（圖6）。壩趾附近上部區(qū)域拉應(yīng)力集中現(xiàn)象是由于附加區(qū)域外伸混凝土的彎曲變形所致。由于該部分結(jié)構(gòu)不在常規(guī)重力壩的基本結(jié)構(gòu)上，不會對重力壩整體穩(wěn)定產(chǎn)生較大不利影響。當(dāng)不考慮應(yīng)力集中區(qū)時，各組壩體的最大拉應(yīng)力都出現(xiàn)在壩體上游左側(cè)底面，見圖7。隨著高度的增加拉應(yīng)力在不斷地減小，并且隨著地基法向彈性模量的增大，壩體的拉應(yīng)力值不斷地減小。

4.5 壓應(yīng)力結(jié)果對比

圖8為第一組的壩體主壓應(yīng)力云圖，各組壩體的主壓應(yīng)力分布形式大致相同，并且都在壩體220截面下游拐角處出現(xiàn)了壓應(yīng)力集中區(qū)域。當(dāng)不考慮應(yīng)力集中區(qū)域時各組壩體的最大壓應(yīng)力都出現(xiàn)在壩體下游面底部。

表6為各組非應(yīng)力集中區(qū)最大壓應(yīng)力值對比表。對比第五組和第六組數(shù)據(jù)可以看出二者差值不超過0.2MPa，且差值比例為2.3%，小于5%。當(dāng)水平向彈性模量與法向彈性模量之比小于1.15時，地基橫觀各向同性對壩體壓應(yīng)力的影響已經(jīng)很小。

5 結(jié)論

通過以上五組對比可以得到以下幾個結(jié)論，一是斜坡段地基中，當(dāng)?shù)鼗臋M向彈性模量與法向彈性模量之比小于1.15時，可以不考慮地基的各向異性性質(zhì)，以各向同性勻質(zhì)地基替代計算即可。二是壩體應(yīng)力大小與地基的彈性模量有關(guān)，從應(yīng)力對比表中可以看出，當(dāng)?shù)鼗姆ㄏ驈椥阅Ａ吭龃髸r壩體的拉應(yīng)力值和壓應(yīng)力值都會減小。因此在混凝土重力壩的有限元分析中，當(dāng)?shù)鼗臋M向彈性模量和豎向彈性模量差異較大時，在計算中應(yīng)當(dāng)考慮地基的橫觀各向同性，不能再簡單的將其視為各向同性體，否則計算結(jié)果將會出現(xiàn)較大的偏差。

參考文獻(xiàn)（References）：

[1] 程心恕，劉國明，蘇言.高等水工建筑物[M].北京：中國水利水電出版社，2011.（CHENG Xin-nu，LIU Guo-ming，SU Yan.The Higher Hydraulic Structures[M].Beijing：China Water Conservancy and Hydropower Press，2011.（in Chinese））

[2] 李麗.碾壓混凝土筑壩技術(shù)在世界的發(fā)展[J].中國水利.2000（10）：23-24.（LI Li.The development of Roller compacted concrete damming technology in the world[J].China Hydraulic.2000（10）：23-24 （in Chinese））

[3] 王松春，中國碾壓混凝土筑壩技術(shù)及其發(fā)展[J].水利水電科技進(jìn)展，2001（2）：16-18（WANG Song-chun.Technology and Development of Roller-Compacted Concrete Dams in China[J].Advances In Science and Technology of Water Resources.2001（2）：16-18）（in Chinese））DOI：10.3880/j.issn.1006-7647.2001.02.006

[4] 王海波，楊會臣，混凝土重力壩全壩段三維動力分析[J].水力發(fā)電學(xué)報，2011.30（6）：133-137.（WANG Hai-bo，YANG Hui-cheng，Seismic analyses of concrete gravity dam with three-dimensional full-dam model[J].Journal of Hydroelectric Power.2011.30（6）：133-137.（in Chinese））

[5] 董玉文，喻杰，蘇琴.碾壓混凝土重力壩碾壓層厚及層面力學(xué)性能的數(shù)值分析[J].水電能源科學(xué)，2014（5）：56-58.（DONG Yu-wen，YU Jie，SU Qin.Numerical Analysis of Layer Thickness of RCC Gravity Dam and Mechanical Properties of Layer Surface[J].Water Resources and Power，2014（5）：56-58.（in Chinese））

[6] 魏博文，徐鎮(zhèn)凱，徐寶松.碾壓混凝土壩層面影響帶黏彈塑性流變模型[J].水利學(xué)報，2012（9）：1097-1102.（WEI Bo-wen，XU Zhen-kai，XU Bao-song.Viscoelasto-plastic rheological model of effect zone for RCCD[J].Journal of Hydraulic Engineering，2012（9）：1097-1102.（in Chinese））

[7] 商兆濤.碾壓混凝土重力壩靜力及動力特性研究[D].合肥：合肥工業(yè)大學(xué)，2009.（SHANG Zhao-tao.Roller compacted concrete gravity dam static and dynamic research[D].Hefei：Hefei University of Technology，2009.（in Chinese））

[8] 燕喬，宋志誠，張利雷，等.碾壓混凝土壩碾壓層厚對層間面特性影響研究[J].三峽大學(xué)學(xué)報：自然科學(xué)版.2015（01）：11-13.（YAN Qiao，SONG Zhi-cheng，ZHANG Li-lei etc.Influence of RRC layer thickness on interlayer surface features[J].Three Gorges University Journal：Science and Technology.2015（01）：11-13.（in Chinese））DOI：10.13393/j.cnki.issn.1672-948X.2015.01.003

[9] 劉光廷，郝巨濤.碾壓混凝土拱壩壩體應(yīng)力的簡化計算[J].清華大學(xué)學(xué)報：自然科學(xué)版，1996（1）：27-33.（LIU Guang-yan，HAO Ju-tao.Simplified numerical method in stress analysis of RCC arch dam with layered structure[J].Journal of Tsinghua University：Science and Technology，1996（1）：27-33.（in Chinese））

[10] Nonlinear simulation of arch dam cracking with mixed finite element method[J].Water Science and Engineering，2008，02：88-101.

[11] 莊茁，由小川，廖劍暉，等，基于ABAQUS的有限元分析和應(yīng)用[M].北京：清華大學(xué)出版社，2009.（ZHUANG Chu，YOU Xiao-chuang，LIAO Jian-hui，et al.The finite element analysis and application based on ABAQUS[M].Beijing：Tsinghua University Press，2009.（in Chinese））

[12] 王金昌，陳頁開.ABAQUS在土木工程中的應(yīng)用[M].杭州：浙江大學(xué)出版社，2006.（WANG Jin-chang，CHEN Ye-kai.The application of ABAQUS in civil engineering[M].Hangzhou：Zhejiang University Press，2006.（in Chinese））

[13] 朱向榮，王金昌.ABAQUS軟件中部分土模型簡介及其工程應(yīng)用[J].巖土力學(xué).2004（S2）：144-148.（ZHU Xiang-rong，WANG Jing-chang.Introduction to partly soil models in ABAQUS Software and their application to the geotechnical engineering[J].Rock and Soil Mechanics.2004（S2）：144-148.（in Chinese））DOI：10.3969/j.issn.1000-7598.2004.z2.029

[14] 李美蓉.復(fù)雜地質(zhì)條件下多結(jié)構(gòu)面對重力壩壩基穩(wěn)定性的影響及處理[J].巖土力學(xué)，2014（S1）：328-333.（LI Mei-rong，Impacts of multi-structural planes on gravity dam foundation stability and treatment under complex geological conditions[J].Rock and Soil Mechanics.2014（S1）：328-333（in Chinese））

[15] 徐干成，李成學(xué)，劉平.各向異性和非均質(zhì)地基土上淺基礎(chǔ)的極限承載力[J].巖土工程學(xué)報.2007，29（2）：164-168.（XU Gan-cheng，LI Cheng-xue，LIU Ping.Ultimate bearing capacity of shallow footing on anisotropic and nonhomogeneous clays[J].Journal of Geotechnical Engineering.2007，29（2）（in Chinese））DOI：10.3321/j.issn：1000-4548.2007.02.002

[16] 姚仰平，孔玉俠.橫觀各向同性土強(qiáng)度與破壞準(zhǔn)則的研究[J].水利學(xué)報，2012（1）：43-50.（YAO Yang-ping，KONG Yu-xia.Study on strength and failure criterion of cross-anisotropic soil[J].Journal of Hydraulic Engineering，2012（1）：43-50.（in Chinese））

[17] 靳紹巖，吳嘉丞，張久長，等.各向異性地基條件下橋墩基礎(chǔ)沖刷穩(wěn)定性分析[J].四川建筑，2010，30（6）：91-93.（JIN Shao-yan，WU Jia-ceng，ZHANG Jiu-chang et al.Under the condition of anisotropic foundation pier foundation scour stability analysis[J].Sichuan Architecture，2010，30（6）：91-93.（in Chinese））DOI：10.3969/j.issn.1007-8983.2010.06.037

[18] 阮懷寧，沈家蔭.水工建筑各向異性地基極限平衡分析[J].大連理工大學(xué)學(xué)報1993（S1）：62-67.（YUAN Huai-ning，SHENG Jia-meng.Limit equilibrium analysis of anisotropic foundation in hydreoelectric engineering[J].Journal of Dalian University of Technology 1993（S1）：62-67.（in Chinese））

[19] 李婕.橫觀各向同性地基空間問題的位移函數(shù)解法[J].巖土工程學(xué)報，2007（1）：137-142.（LI Jie，Displacement function method of space problem for transversely isotropic foundation.[J].Journal of Geotechnical Engineering，2007（1）：137-142.（in Chinese））DOI：10.3321/j.issn：1000-4548.2007.01.024

[20] 艾智勇，成怡沖.三維橫觀各向同性成層地基的傳遞矩陣解[J].巖土力學(xué)，2010（S2）：25-30.（AI Zhi-yong，CHENG Yi-chong.Transfer matrix solutions of three-dimensional transversely isotropic multilayered soils[J].Rock and Soil Mechanics，2010（S2）：25-30.（in Chinese））

[21] 艾智勇，李博.橫觀各向同性層狀地基平面應(yīng)變問題的解析層元解[J].巖土工程學(xué)報，2012（10）：1787-1791（AI Zhi-yong，LI Bo.Analytical layer element solutions to plane strain problem of transversely isotropic multilayered soils[J].Journal of Geotechnical Engineering 2012（10）：1787-1791 （in Chinese））