瀝青混凝土心墻壩三維有限元靜動力分析

孔憲京，余 翔，鄒德高1，，周 揚

（1.大連理工大學(xué) 海岸和近海工程國家重點實驗室，遼寧 大連 116024；2.大連理工大學(xué) 水利工程學(xué)院，遼寧 大連 116024）

0 引 言

土石壩是目前最常用的壩型，而以水工瀝青混凝土作為心墻防滲材料的土石壩優(yōu)勢越來越明顯，已成為一種非常有競爭力的壩型［1］.水工瀝青混凝土用于土石壩內(nèi)部防滲已近50年［2］，國際大壩委員會（ICOLD）分別在第42號公報［3］與第84號公報［4］中總結(jié)了這種壩型的設(shè)計與施工經(jīng)驗.在1988年第16屆國際大壩會議上，與會專家一致認為以瀝青作為防滲體的土石壩和混凝土面板堆石壩是未來高壩的適宜壩型.與混凝土面板壩［5］相比，瀝青混凝土心墻壩仍處于發(fā)展階段，相對滯后，尤其在高瀝青混凝土心墻壩的建設(shè)方面缺乏實踐經(jīng)驗.

隨著瀝青混凝土心墻壩在土石壩中逐步得到廣泛應(yīng)用，國內(nèi)外學(xué)者開展了瀝青混凝土心墻壩的數(shù)值計算分析.但目前有關(guān)瀝青混凝土心墻壩的數(shù)值計算以二維［6－8］居多，不能反映三維河谷效應(yīng)；數(shù)值分析均依托實際工程［9－10］開展，主要用于評價具體工程，無法把握瀝青混凝土心墻的響應(yīng)規(guī)律；動力分析主要針對壩體加速度和心墻應(yīng)力反應(yīng)［11－12］，而有關(guān)心墻的動剪應(yīng)變反應(yīng)規(guī)律鮮有涉及.等效線性法［13］是目前土石壩工程中廣泛采用的動力分析方法，動剪應(yīng)變是該法的重要參數(shù)，但目前有關(guān)水工瀝青混凝土的動力試驗只有少數(shù)學(xué)者涉及［14－15］，不能很好地把握動剪應(yīng)變的基本范圍.若根據(jù)已有試驗成果，通過數(shù)值敏感性計算分析，確定心墻的最大動剪應(yīng)變的分布范圍和規(guī)律，對室內(nèi)試驗有一定的指導(dǎo)意義.同時，心墻的動剪應(yīng)變較大區(qū)域也是心墻地震反應(yīng)強烈區(qū)，熟悉其分布規(guī)律也有助于工程實踐.鑒于此，本文開展不同壩高的瀝青混凝土心墻壩三維有限元靜動力數(shù)值計算，針對瀝青混凝土心墻，分析并總結(jié)其應(yīng)力與變形和動剪應(yīng)變規(guī)律，以期為水工瀝青混凝土試驗與瀝青混凝土心墻壩建設(shè)提供依據(jù).

1 計算模型及參數(shù)

本文采用大連理工大學(xué)工程抗震研究所開發(fā)的非線性有限元計算程序——GEODYNA［16］進行計算，分析了靜力狀態(tài)下（竣工和滿蓄期）瀝青混凝土心墻的應(yīng)力與變形以及遭遇地震時的動剪應(yīng)變分布特性及其規(guī)律.非線性分析中采用的本構(gòu)模型、計算參數(shù)等介紹如下.

1.1 本構(gòu)模型

采用鄧肯E－B模型［17］進行靜力計算，主要計算公式及參數(shù)如下：

式中：Et為切線變形模量，需K、n、φ、c和Rf5個試驗常數(shù)；Eur為卸載及再加載的彈性模量，需Kur和n兩個試驗常數(shù)；B為體積模量，需Kb和m兩個試驗常數(shù)；pa為大氣壓，量綱與σ3相同.

周恩來在招待會上說，建設(shè)新民航，人才是主要的。林雨水和其他兩航起義人員都表示要積極參加新中國的建設(shè)，為祖國的航空事業(yè)貢獻力量。在隨后的歲月里，林雨水先后在空軍、海軍航空兵、民航任飛行員和教員，一直珍藏著那張留有周恩來簽名的法幣，并于1979年將法幣捐獻給了中共代表團梅園新村紀念館。

運用等效線性模型［18］進行動力分析.其最大動剪切模量表示為

由于三維計算內(nèi)容較多，為便于對比，主要對河谷心墻中央斷面與心墻壩軸線剖面的應(yīng)力、變形和動剪應(yīng)變進行分析.圖中應(yīng)力以壓為正，順河向位移以向下游為正，豎向位移以向上為正.

1.2 計算模型

假設(shè)大壩坐落于基巖上，計算模型采用的瀝青混凝土心墻均質(zhì)壩的壩高分別為100、150 和200m，上下游壩坡均為1∶2.0，河谷對稱，坡度為1∶1.5，河谷寬80m，心墻厚度沿壩高線性變化，其中心墻頂部厚度取0.5m，底部厚度取壩高的1／100，筑壩材料采用堆石料.圖1為標(biāo)準(zhǔn)計算模型網(wǎng)格圖.

圖1 瀝青混凝土心墻壩三維有限元網(wǎng)格Fig.1 3D－FE mesh of asphalt concrete core dam

目前，有關(guān)堆石料與瀝青混凝土之間的接觸性質(zhì)的研究非常少見.因此，本文計算中并未在堆石體與瀝青混凝土心墻之間設(shè)置接觸面，且當(dāng)不設(shè)置接觸面時，心墻的“懸掛”效應(yīng)更明顯，心墻的應(yīng)力結(jié)果是偏于安全的.

而如果用番茄沙司代替辣椒來烹制川湘味海鮮，效果就大不相同了，既能在視覺上給人紅辣熱烈的感覺，吃起來辣味輕淡（番茄沙司配少量辣椒烹制）鮮味猶存，讓喜歡吃川湘菜的食客也能愛上小海鮮。

1.3 計算參數(shù)

圖9為不同壩高的河谷中央斷面心墻靜位移沿壩高的分布曲線圖.可以看出，由于心墻位于大壩中部，心墻竣工期順河向位移基本為零，心墻滿蓄期順河向位移最大值位于心墻頂部.心墻豎向沉降最大值基本位于心墻中部.心墻位移最大值與壩高并非簡單地呈比例關(guān)系.滿蓄期，上游堆石受浮托力作用，而堆石控制心墻的整體變形，心墻將隨上游堆石一起向壩頂方向移動，心墻豎向沉降減小.

表1 靜力分析參數(shù)Tab.1 Parameters of static analysis

表2 動力分析參數(shù)Tab.2 Parameters of dynamic analysis

圖2 堆石料動剪切模量比與動剪應(yīng)變關(guān)系Fig.2 Relationship between dynamic shear modulus ratio and dynamic shear strain of rockfill materials

圖3 堆石料阻尼比與動剪應(yīng)變關(guān)系Fig.3 Relationship between damping ratio and dynamic shear strain of rockfill materials

圖4 瀝青混凝土動剪切模量比與動剪應(yīng)變關(guān)系Fig.4 Relationship between dynamic shear modulus ratio and dynamic shear strain of asphalt concrete

1.4 地震輸入與計算工況

當(dāng)最大壩高增大（相應(yīng)的工況為1、2 和3）時，心墻的最大動剪應(yīng)變極值變化不大，但最大動剪應(yīng)變極值位置更靠近壩頂，這是因為壩越高壩頂部甩動越明顯.

圖5 規(guī)范譜人工波地震時程曲線Fig.5 Seismic time history curves of artificial waves of standard spectrum

圖6 兩河口地震時程曲線Fig.6 Seismic time history curves of Lianghekou

圖7 糯扎渡地震時程曲線Fig.7 Seismic time history curves of Nuozhadu

本文有關(guān)計算工況列于表3.其中考慮了瀝青混凝土的動剪切模量、壩高、地震波，以及加速度峰值的變化.

圖8 地震波加速度反應(yīng)譜Fig.8 Acceleration amplification spectra of seismic waves

表3 計算工況Tab.3 Cases for calculation

2 計算結(jié)果及分析

式中：K、n由試驗參數(shù)確定，σ′m＝（σ′1＋σ′2＋σ′3）／3，σ′1、σ′2、σ′3為試驗時試樣所受的有效主應(yīng)力.可根據(jù)試驗確定動剪應(yīng)變γ與G／Gmax及阻尼比λ的關(guān)系.

2.1 靜力計算結(jié)果與分析

堆石料［19］與瀝青混凝土靜力分析參數(shù)如表1所示，表2為其動力分析參數(shù)，瀝青混凝土靜動力計算參數(shù)取自西安理工大學(xué)試驗成果［20］.動力計算時，堆石料采用如圖2與3中孔憲京等［18］建議的動剪切模量和阻尼比與動剪應(yīng)變的歸一化關(guān)系的平均值.瀝青混凝土試驗結(jié)果表明，不同固結(jié)比和不同圍壓條件下的阻尼比λ呈帶狀分布，且圍壓對阻尼比影響不大，瀝青混凝土的阻尼比λ在0.10～0.15.本文根據(jù)實驗結(jié)果［20］，將瀝青混凝土的阻尼比λ取為0.121.圖4為溫度（10.5±0.5）℃下瀝青混凝土的動剪切模量比與動剪應(yīng)變歸一化關(guān)系曲線.

圖10是不同壩高的河谷中央斷面心墻豎向應(yīng)力沿壩高的分布曲線圖.由于心墻較薄，心墻的變形受堆石控制，且心墻變形模量小于堆石，因此心墻的變形能力不能充分發(fā)揮，造成心墻豎向應(yīng)力小于自重應(yīng)力.竣工期，豎向應(yīng)力沿壩高遞增分布，在壩基處達到最大.滿蓄期，受浮托力影響，心墻豎向應(yīng)力降低，心墻高度1／5以下區(qū)域，豎向應(yīng)力逐步降低.總的來說，滿蓄期，河谷中部的心墻下部區(qū)域豎向應(yīng)力較低，且高心墻將處于相對更不利的應(yīng)力狀態(tài).

圖9 河谷中央斷面心墻靜位移沿壩高的分布Fig.9 The distribution of static displacement of the core in the valley central section along the dam height

圖10 河谷中央斷面心墻豎向應(yīng)力沿壩高的分布Fig.10 The distribution of vertical stress of the core in the valley central section along the dam height

應(yīng)力水平是判別材料極限平衡區(qū)的一項重要指標(biāo).計算結(jié)果表明，不同壩高的心墻應(yīng)力分布規(guī)律相同.本文僅給出壩高為200m 的壩軸向剖面心墻應(yīng)力水平等值線圖，見圖11.竣工期，最大應(yīng)力水平位于河谷處心墻中下部.滿蓄期，受浮托力影響，心墻豎向應(yīng)力降低，而心墻順河向應(yīng)力變化較小，使河谷處心墻中部應(yīng)力水平大幅降低.壩肩處、頂部及底部心墻應(yīng)力狀態(tài)復(fù)雜，滿蓄期時，這些部位單元的主應(yīng)力方向會發(fā)生偏轉(zhuǎn)，造成應(yīng)力水平增大.表4為不同壩高的心墻應(yīng)力水平極值表.可以看出，隨著壩高的增加，心墻的應(yīng)力水平極值增大，且滿蓄期相對竣工期增幅變大.

2017 年長沙縣完成地區(qū)生產(chǎn)總值1431.1億元，增長11.3%；規(guī)模工業(yè)總產(chǎn)值2662億元，增長12.9%，固定資產(chǎn)投資955.9億元，增長17.2%；財政總收入280.1億元，增長23.2%，率先全省實現(xiàn)全面小康。在全國縣域經(jīng)濟基本競爭力、中國中小城市綜合實力百強排名中均躍居第六位?，F(xiàn)在的長沙縣正乘著改革開放的春風(fēng)，朝著全國縣域經(jīng)濟五強奮進！

圖11 壩軸向剖面心墻應(yīng)力水平分布圖（壩高：200m）Fig.11 The distribution of stress level of the core along the cross section of the dam axis（dam height：200m）

表4 不同壩高的心墻應(yīng)力水平極值Tab.4 Maximum stress level of the core with different dam heights

2.2 動剪應(yīng)變計算結(jié)果與分析

對同一個心墻壩而言，其自振頻率一定，而地震波具有隨機性，即使峰值加速度一定，頻譜特性也不盡相同，因此心墻的反應(yīng)是不同的.為研究地震波變化對心墻動力反應(yīng)的影響，分別采用規(guī)范譜人工波、兩河口，以及糯扎渡場地譜地震波（相應(yīng)的工況為3、9和10）進行動力分析.從圖中可以看出，不同地震波對最大動剪應(yīng)變的分布規(guī)律影響不大，但最大動剪應(yīng)變極值差別較大.

本文分別采用水工抗震設(shè)計規(guī)范譜生成的人工波、兩河口和糯扎渡土石壩工程的場地譜地震波為動力輸入，研究心墻在不同地震波下的動力反應(yīng).圖5～7為3種地震波的時程曲線.圖8為不同地震波的反應(yīng)譜.

試驗條件對瀝青混凝土的動力參數(shù)影響較大.為分析瀝青混凝土動力參數(shù)對心墻動力反應(yīng)的影響，瀝青混凝土的動剪切模量系數(shù)分別取979.4、1 979.4和2 979.4（相應(yīng)的工況為4、3和5）.隨著瀝青混凝土動剪切模量系數(shù)的增大，心墻動剛度增大，心墻的最大動剪應(yīng)變極值減小，同時岸坡處心墻動力反應(yīng)增強，動剪應(yīng)變增大.

圖12 壩軸向剖面心墻最大動剪應(yīng)變分布（單位：%）Fig.12 The distribution of maximum dynamic shear strain of the core along the cross section of the dam axis（unit：%）

雖然被稱為“行業(yè)新兵”，但早在公司成立的三年前，主創(chuàng)團隊就已經(jīng)開始研發(fā)與智能機器人相關(guān)的軟件和算法。黃總認為，相比硬件來說，后臺管理軟件才是機器人系統(tǒng)的靈魂?！拔覀円獙崿F(xiàn)門到門的產(chǎn)品解決方案，就必須不斷堅持研發(fā)，針對客戶倉庫中出現(xiàn)的問題來推出能夠更好為他們解決問題的產(chǎn)品解決方案。比如我們將在未來逐漸推出的自動電拖、輕型車等產(chǎn)品，就是新的嘗試。我相信在明年的CeMAT展會上，觀眾能夠看到牧星帶來的更多新產(chǎn)品，也會看到我們發(fā)布更加重磅的新技術(shù)?！?/p>

為了分析加速度峰值的影響，采用規(guī)范譜地震波，地震加速度峰值分別為0.2g、0.3g、0.4g和0.5g（相應(yīng)的工況為6、7、3 和8）.隨著地震動峰值加速度增大，心墻動力反應(yīng)增強，心墻頂部的最大動剪應(yīng)變極值增大；當(dāng)?shù)卣饎臃逯禐?.5g時，最大動剪應(yīng)變極值與0.4g時相差較小，但心墻岸坡處動剪應(yīng)變較大.

不同工況下的心墻最大動剪應(yīng)變分布規(guī)律基本一致，河谷處心墻頂部區(qū)域的最大動剪應(yīng)變極值均較大.圖12為典型工況下的壩軸向剖面心墻最大動剪應(yīng)變分布圖.

為全面貫徹國家大力發(fā)展特色小鎮(zhèn)建設(shè)精神，《廣西培育特色小鎮(zhèn)意見》提出，到2020年，要培育30個左右全國特色小鎮(zhèn)、100個左右自治區(qū)級特色小鎮(zhèn)、建設(shè)200個左右市級特色小鎮(zhèn)。該意見的出臺為廣西貫徹國家特色小鎮(zhèn)發(fā)展戰(zhàn)略、發(fā)展廣西區(qū)域經(jīng)濟明確了方向。

3 結(jié) 論

（1）瀝青混凝土心墻的變形主要受堆石控制，滿蓄期，壩肩處心墻應(yīng)力水平較高，應(yīng)給予重點關(guān)注.壩增高，心墻頂部向上的豎向位移增大，心墻底部的豎向應(yīng)力降幅增大，壩肩處心墻應(yīng)力水平峰值增大，心墻將處于更不利的應(yīng)力與變形狀態(tài).

胃潰瘍屬于臨床多見病，2型糖尿病一旦合并胃潰瘍，若治療不及時，易導(dǎo)致血糖過高所致葡萄糖毒性，損害人體神經(jīng)系統(tǒng)，進而產(chǎn)生一系列嚴重后果，比如上消化道出血、潰瘍穿孔和幽門梗阻等，嚴重影響患者生活質(zhì)量以及身心健康[3]。對于糖尿病患者來說，其尚無根治方法，需終身服藥治療，而長時間服用降糖藥物，易導(dǎo)致人體胃酸大量分泌，而大量分泌胃酸易損害人體胃粘膜，繼而形成胃潰瘍，因此，對于2型糖尿病合并胃潰瘍的治療，應(yīng)以降糖治療和抑制胃酸分泌為基礎(chǔ)。

（2）地震時最大動剪應(yīng)變發(fā)生在河谷中部的心墻頂部區(qū)域，但其極值不大.在加速度峰值為0.5g的規(guī)范譜人工波激勵下，最大動剪應(yīng)變極值約為0.2%；當(dāng)采用加速度峰值為0.4g的兩河口地震波激勵時，最大動剪應(yīng)變極值約為0.3%.因此可以認為，對水工瀝青混凝土進行室內(nèi)試驗時，最大動剪應(yīng)變達到0.5%即可滿足工程要求.

生態(tài)修復(fù)是針對采煤沉陷區(qū)治理的重要方式。根據(jù)沁水縣的自然資源特點和采煤沉陷區(qū)的土地破壞程度，采用地貌重塑、土壤重構(gòu)以及植被重建等方法，對采煤沉陷區(qū)進行治理［8］。對采煤沉陷區(qū)的居民采取采煤沉陷區(qū)搬遷安置，拆除沉陷區(qū)的危舊房屋，并實施土地復(fù)墾。安置時要進行相關(guān)公共服務(wù)設(shè)施的配套建設(shè)，同時創(chuàng)造一定的就業(yè)崗位，促進農(nóng)民就地就近安居和轉(zhuǎn)移產(chǎn)業(yè)。

（3）瀝青混凝土心墻在靜力和動力狀態(tài)下，頂部區(qū)域的變形和動剪應(yīng)變均處于相對不利狀態(tài).因此，設(shè)計時應(yīng)引起注意.

本文研究中，假設(shè)模型均坐落在基巖上，深厚覆蓋層上的瀝青混凝土心墻壩有待于進一步研究.

［1］ WANG Wei－biao，Hoeg K.The asphalt core embankment dam：a very competitive alternative［C］／／ Modern Rockfill Dams — 2009.Beijing：China Water Power Press，2009.

［2］ Hoeg K.Asphaltic Concrete Cores for Embankment Dams［M］.Norway：Stikka Press，1993.

［3］ ICOLD.Bituminous cores for earth and rockfill dam［Z］.Germany：International Commission on Large Dam，1982.

［4］ ICOLD.Bituminous cores for fill dams－state of art［Z］.Paris：International Commission on Large Dams，1992.

［5］ 楊澤艷，周建平，蔣國澄，等.中國混凝土面板堆石壩的發(fā)展［J］.水力發(fā)電，2011，37（2）：18－23.YANG Ze－yan，ZHOU Jian－ping，JIANG Guocheng，etal.Development of CFRD in China［J］.Water Power，2011，37（2）：18－23.（in Chinese）

［6］ 張丙印，李全明，熊 焰，等.三峽茅坪溪瀝青混凝土心墻堆石壩應(yīng)力變形分析［J］.長江科學(xué)院院報，2004，21（2）：18－21.ZHANG Bing－yin，LI Quan－ming，XIONG Yan，etal.Numerical study on Maopingxi rockfill dam with asphalt concrete core in the Three Gorges Project［J］.Journal of Yangtze River Scientific Research Institute，2004，21 （2）：18－21.（in Chinese）

［7］ 李同春，劉曉青，夏頌佑，等.冶勒堆石壩瀝青混凝土心墻型式及尺寸研究［J］.河海大學(xué)學(xué)報，2000，28（2）：109－112.LI Tong－chun，LIU Xiao－qing，XIA Song－you，etal.Study on asphalt concrete core type and size of Yele rockfill dam ［J］.Journal of Hohai University，2000，28（2）：109－112.（in Chinese）

［8］ Feizi－Khankandi S，F(xiàn)halandarzadeh A，Mirghasemi A A，etal.Seismic analysis of the Garmrood embankment dam with asphalt concrete core［J］.Soils and Foundations，2009，49（2）：153－166.

［9］ 鄧傳貴，沈振中，凌春海.高瀝青混凝土心墻壩動力反應(yīng)特性［J］.南水北調(diào)與水利科技，2012，10（5）：147－152.DENG Chuan－gui，SHEN Zhen－zhong，LING Chun－h(huán)ai.Analysis of dynamic response characteristics for core rockfill dam with high asphalt concrete ［J］.South－to－North Water Transfers and Water Science ＆ Technology，2012，10（5）：147－152.（in Chinese）

［10］ 劉武軍，陳堯隆，司 政，等.瀝青混凝土心墻壩三維有限元動力分析［J］.水資源與水工程學(xué)報，2009，20（6）：152－156.LIU Wu－jun，CHEN Yao－long，SI Zheng，etal.3－D finite element dynamic analysis of asphalt concrete core dam ［J］.Journal of Water Resources ＆Engineering，2009，20（6）：152－156.（in Chinese）

［11］ 朱 晟.瀝青混凝土心墻堆石壩三維地震反應(yīng)分析［J］.巖土力學(xué)，2008，29（11）：2933－2938.ZHU Sheng.3－D seismic response analysis of rockfill dam with asphalt concrete core［J］.Rock and Soil Mechanics，2008，29（11）：2933－2938.（in Chinese）

［12］ Baziar M H，Salemi S，Heidari T.Analysis of earthquake response of an asphalt concrete core embankment dam［J］.International Journal of Civil Engineering，2006，4（3）：192－211.

［13］ 顧淦臣，沈長松，岑威鈞.土石壩地震工程學(xué)［M］.北京：中國水利水電出版社，2009.GU Gan－chen，SHEN Chang－song，CEN Wei－jun.Earthquake Engineering for Earthrock Dam ［M］.Beijing：China Water Power Press，2009.（in Chinese）

［14］ WANG Wei－biao，Hoeg K.Cyclic behavior of asphalt concrete used as impervious core in embankment dams［J］.Journal of Geotechnical and Geoenvironmental Engineering，2011，137（5）：536－544.

［15］ Feizi－Khankandi S，F(xiàn)halandarzadeh A，Mirghasemi A A，etal.Cyclic triaxial tests on asphalt concrete as a water barrier for embankment dams［J］.Soil and Foundations，2008，48（3）：319－332.

［16］ 鄒德高，孔憲京，徐 斌.GEOtechnical DYnamic Nonlinear Analysis—GEODYNA 使用說明［R］.大連：大連理工大學(xué)工程抗震研究所，2008.ZOU De－gao，KONG Xian－jing，XU Bin.GEOtechnical DYnamic Nonlinear Analysis —GEODYNA user manual ［R］.Dalian：The Earthquake Engineering Research Institute，Dalian University of Technology，2008.（in Chinese）

［17］ Duncan J M，Chang Chin－yung.Nonlinear analysis of stress and strain in soils［J］.Journal of the Soil Mechanics and Foundations Division，ASCE，1970，96（5）：1629－1653.

［18］ 孔憲京，婁樹蓮，鄒德高，等.筑壩堆石料的等效動剪切模量與等效阻尼比［J］.水利學(xué)報，2001（8）：20－25.KONG Xian－jing，LOU Shu－lian，ZOU De－gao，etal.The equivalent dynamic shear modulus and equivalent damping ratio of the rockfill for dam［J］.Journal of Hydraulic Engineering，2001（8）：20－25.（in Chinese）

［19］ 酈能惠，孫大偉，李登華，等.300 m 級超高面板堆石壩變形規(guī)律的研究［J］.巖土工程學(xué)報，2009，31（2）：155－160.LI Neng－h(huán)ui，SUN Da－wei，LI Deng－h(huán)ua，etal.Deformation behavior of 300 m high－concrete face rockfill dams［J］.Chinese Journal of Geotechnical Engineering，2009，31（2）：155－160.（in Chinese）

［20］ 西安理工大學(xué).新疆五一水庫試驗結(jié)論［R］.西安：西安理工大學(xué)，2008.Xi′an University of Technology.Test conclusion of Xinjiang Wuyi Reservoir ［R］.Xi′an：Xi′an University of Technology，2008.（in Chinese）