張 強，劉慧娟，李文柱，舒華生

（1.昆明市水利水電勘測設(shè)計研究院，云南昆明650231；2.廣西水利電力職業(yè)技術(shù)學(xué)院，廣西南寧530021）

有限元分析在高烈度地區(qū)土石壩設(shè)計中的應(yīng)用

張 強1，劉慧娟2，李文柱1，舒華生1

（1.昆明市水利水電勘測設(shè)計研究院，云南昆明650231；2.廣西水利電力職業(yè)技術(shù)學(xué)院，廣西南寧530021）

本文采用鄧肯E-B雙曲線非線性模型為本構(gòu)模型，對粘土心墻堆石壩分別進(jìn)行二維和三維有限元分析，計算結(jié)果均較好的反映了地震作用下壩體的位移和應(yīng)力分布情況，計算結(jié)果表明三維有限元分析位移和應(yīng)力分布規(guī)律與二維分析完全吻合，但三維分析位移和應(yīng)力量值明顯小于二維分析，三維分析垂直位移相對二維分析偏小，同時三維計算未能反映濕化產(chǎn)生的沉降位移。由于受地形限制，壩體斷面形狀受山體影響較多，三維計算結(jié)果較為真實地反映了復(fù)雜地形變化特性對壩體地震響應(yīng)的影響，對高烈度地區(qū)、復(fù)雜地形條件下的土石壩設(shè)計具有一定的借鑒意義。

土石壩；有限元分析；位移；應(yīng)力

1 前言

土石壩由于就地取材和施工比較簡單，是應(yīng)用非常普遍的一種壩型，尤其是中等高度以下的壩，數(shù)量尤其多。對于高烈度地區(qū)、復(fù)雜地形條件下的土石壩設(shè)計，進(jìn)行必要的有限元計算，能較好的反映壩體內(nèi)部及不連續(xù)界面的力學(xué)特性，

本文采用有限元方法對海馬箐水庫土石壩在地震作用下的動力響應(yīng)進(jìn)行分析，研究壩體位移及應(yīng)力分布的變化情況，為大壩的設(shè)計提供定性的設(shè)計參考資料。

2 工程概況

海馬箐水庫工程擋水建筑物為粘土心墻堆石壩。工程規(guī)模為中型，工程等別為Ⅲ等。壩址區(qū)域河道走向呈“S”形，受河道地形彎道的限制，大壩平面上呈現(xiàn)出上、下游不規(guī)則的布置形式。壩頂寬8m，壩軸線長158.8m，壩頂高程1786.40m，建基面清基開挖高程1727.20m，最大壩高59.2m。工程區(qū)地震基本烈度≥Ⅸ度。工程設(shè)計按Ⅸ度設(shè)防。

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

3.1 巖土本構(gòu)模型

非線性彈性模型一般適用于安全系數(shù)較大，不發(fā)生屈服的情況，盡管這與實際情況會有一定的差異，但由于避免了使用彈塑性模型的困難，在土石壩的有限元計算中得到了較為廣泛的應(yīng)用。其中，鄧肯E-B雙曲線非線性模型由于能較好地模擬粘土心墻堆石的變形性狀，計算中土體材料的彈性模量和泊松比隨應(yīng)力的變化而變化，計算參數(shù)測定有比較成熟的經(jīng)驗，因而在水利水電工程中得到了較為廣泛的應(yīng)用。

鄧肯E-B雙曲線非線性模型切線模量Et、卸載時切線彈模Eur和切線體積模量B的表達(dá)式為，

式中：Rf為材料的破壞比；K、Kur、Kb、m、n均為試驗常數(shù)；Pa為工程大氣壓；c和φ為材料的黏聚力與內(nèi)摩擦角。

3.2 荷載處理方式

首先完成在壩基自重條件下的應(yīng)力分析，提取對應(yīng)的壩基自重應(yīng)力場（相應(yīng)的應(yīng)變和變形位移清零）作為后續(xù)分析的初始地應(yīng)力場。在初始地應(yīng)力場的基礎(chǔ)上，本工程將壩體模型分為34個填筑層，采用逐層單元依次激活的方式模擬壩體施工的填筑過程。

由于粘土心墻的滲透系數(shù)很小，可認(rèn)為心墻為基本不透水材料，水庫蓄水后，水壓力作用在心墻上游面。穩(wěn)定滲流形成后，浸潤線以下土體處于飽和狀態(tài)，以上土體處于非飽和狀態(tài)。壩體在滲流作用下將承受滲透作用力，其中上游壩殼料由于為完全透水性材料，滲透作用力可簡化為浮托力。浸潤線以下的土體浸水飽和后，有效應(yīng)力發(fā)生了變化，在向下作用的滲透力分量的作用下，土體將表現(xiàn)出壓縮下沉的位移現(xiàn)象，稱為濕化效應(yīng)。二維有限元計算可完整地考慮水的作用效應(yīng)，但三維有限元計算尚不能考慮濕化效應(yīng)。

二維有限元計算使用擬靜力法，有關(guān)土條地震作用水平慣性力代表值的計算方法參照水工建筑物抗震規(guī)范進(jìn)行。三維有限元分析由于使用的A B A Q U S軟件無與水工抗震規(guī)范規(guī)相一致的擬靜力法可供選用，本次采用了非線性時程分析法。根據(jù)海馬箐壩址區(qū)域的地址場地條件，選擇了如圖1所示的地震加速度時程（使用前將地震峰值加速度調(diào)幅到0.4g）。

圖1 三維計算使用的地震加速度時程（調(diào)幅前）

3.3 幾何模型建立

由于大壩受地形和河道幾何形狀的影響，大壩斷面沿壩軸線變化較大。為能充分體現(xiàn)不同斷面位置的受力及穩(wěn)定性狀態(tài)，根據(jù)壩體幾何特征沿壩軸線的變化情況，在實際計算值選擇了a、b兩個代表性斷面進(jìn)行計算分析。具體見圖2、圖3。

圖2 斷面a二維模型圖

圖3 斷面b二維模型圖

根據(jù)壩址區(qū)域地形及壩體幾何數(shù)據(jù)構(gòu)造的三維地形及壩體幾何實體形狀如圖4、圖5所示。

圖4 壩址三維地形圖

圖5 壩體三維幾何實體圖

3.4 材料參數(shù)

心墻粘土料及石碴料采用鄧肯-張模型描述材料的應(yīng)力-應(yīng)變關(guān)系，其余部位材料采用摩爾-庫倫模型描述。有限元計算使用材料參數(shù)如表1所示。

表1 材料參數(shù)表

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

4.1 三維計算結(jié)果

在正常蓄水位運行情況下，考慮9度地震作用下（水平地震加速度0.4g）的三維有限元計算應(yīng)力及變形位移計算結(jié)果如圖6～圖9所示。

圖6 三維計算壩體順河流方向位移分布云圖（m）

圖7 三維計算壩體垂直向位移分布云圖（m）

根據(jù)計算結(jié)果，對于水平向位移，上游堆石體的位移表現(xiàn)為向上游方向移動，最大位移量36c m，而下游堆石體的位移表現(xiàn)為向下游方向移動，最大位移量43.8c m。壩體表現(xiàn)出沉陷式的變形特征。對于垂直方向的位移，整個壩體不同部位表現(xiàn)出不同的位移方式，有上抬位移，也有下沉位移。最大上抬位移為38.3c m，最大下沉位移18.1c m，壩體變形位移分布十分復(fù)雜。

圖8 壩體小主應(yīng)力分布云圖（kP a）

圖9 壩體大主應(yīng)力等值線（kP a）

三維有限元計算結(jié)果表明，壩體上、下游堆石區(qū)壩面馬道附近分布有拉性應(yīng)力，量值一般小于50kP a，但在壩頂中部以及兩壩肩建基面附近的較小區(qū)域內(nèi)，最大拉性應(yīng)力超過了300kP a。

4.2 二維和三維計算結(jié)果對比

斷面a二維有限元計算壩體位移和應(yīng)力分布如圖10～圖13所示。壩體最大水平位移56.9c m，發(fā)生在壩頂靠下游側(cè)區(qū)域；最大豎向位移16c m，發(fā)生在壩頂靠下游側(cè)區(qū)域。上游壩體堆石區(qū)較大范圍及心墻靠上游側(cè)的局部區(qū)域出現(xiàn)了拉性應(yīng)力，堆石區(qū)拉性應(yīng)力值為100-400kP a；靠上游側(cè)心墻區(qū)拉性應(yīng)力為100-300kP a。這些拉性應(yīng)力可能導(dǎo)致壩體上游側(cè)堆石區(qū)和靠上游側(cè)心墻局部開裂。

三維有限元計算壩體地震響應(yīng)特性較為復(fù)雜。根據(jù)斷面a的位移和應(yīng)力切片結(jié)果，壩體最大順河向位移下游面28.3c m，最大垂直下沉位移18.7c m。整個壩頂區(qū)域出現(xiàn)較大的拉性應(yīng)力分布區(qū)域，但應(yīng)力值不大，最大拉性應(yīng)力值約90kP a。

圖10 斷面a二維計算順河向位移等值線（m）

圖11 斷面a二維計算豎向位移等值線（m）

圖12 斷面a二維計算大主應(yīng)力等值線（kP a）

圖13 斷面a二維計算9度地震小主應(yīng)力等值線（kP a）

斷面b二維有限元計算壩體位移和應(yīng)力分布如圖14～圖17所示。最大水平位移52c m，分別發(fā)生在上、下游壩坡中部；最大豎向位移16c m，發(fā)生在壩頂部靠上游側(cè)。上游壩體堆石區(qū)一定范圍出現(xiàn)了拉性應(yīng)力，最大應(yīng)力值約為50kP a；上游側(cè)建基面附近靠上游側(cè)心墻出現(xiàn)了約50kP a的拉性應(yīng)力。這些拉性應(yīng)力區(qū)可能導(dǎo)致壩體和局部心墻開裂，但心墻未出現(xiàn)貫穿性的裂縫，裂縫僅發(fā)生在靠上游的局部區(qū)域。

根據(jù)斷面b的位移和應(yīng)力切片結(jié)果，三維有限元壩體最大順河向位移下游面28.3c m，上游面13.6c m；最大垂直下沉位移3.2c m，上抬位移18c m。上游側(cè)建基面附近心墻與過渡層接觸區(qū)域分布有少量拉性應(yīng)力區(qū)，最大拉性應(yīng)力分260kP a。

圖14 斷面b二維計算順河向位移等值線（m）

圖15 斷面b二維計算豎向位移等值線（m）

圖16 斷面b二維計算大主應(yīng)力等值線（kP a）

圖17 斷面b二維計算小主應(yīng)力等值線（kP a）

二維和三維計算均表明，9度地震作用下，壩體出現(xiàn)較大的水平變形位移，壩體不同程度出現(xiàn)拉性應(yīng)力區(qū)，局部區(qū)域拉性應(yīng)力超過了200kP a，這些拉性應(yīng)力可能會導(dǎo)致上游堆石區(qū)、下游堆石區(qū)、部分心墻區(qū)域以及壩頂區(qū)域局部開裂。地震作用下粘土心墻大部分區(qū)域仍保持為壓性應(yīng)力狀態(tài)，未發(fā)現(xiàn)粘土心墻中有貫穿性的拉性應(yīng)力區(qū)，即心墻未出現(xiàn)貫穿性的裂縫。并且，由于受地形限制，壩體斷面形狀受山體影響較多，尤其是凸入壩體內(nèi)部的山體部分或壩體附著在山體的部分受到的影響較大，使得這些區(qū)域附近出現(xiàn)了較高局部拉性應(yīng)力區(qū)。

同時可以看出，三維有限元分析位移和應(yīng)力分布規(guī)律與二維分析基本吻合，但三維分析位移和應(yīng)力量值明顯小于二維分析，堆石體馬道附近區(qū)域在浮托力作用下出現(xiàn)一定量的上抬位移。三維分析垂直位移相對二維分析偏小的主要原因是三維模型約束剛度明顯大于二維，同時三維計算未能反映濕化產(chǎn)生的沉降位移。

5 結(jié)論

本文采用鄧肯E-B雙曲線非線性模型為本構(gòu)模的管理，抓好校審制的崗位職責(zé)，避免校審記錄流于形式，并能真正起到追溯作用，對出院的設(shè)計成果減少會簽評審，嚴(yán)格會議評審。

（3）加強技術(shù)標(biāo)準(zhǔn)、強條的培訓(xùn)及管理工作。在計算書、報告、圖紙的流程卡中明確填寫所采用的主要標(biāo)準(zhǔn)、涉及的強條的內(nèi)容，使校審各環(huán)節(jié)均能對照檢查，從根本上解決技術(shù)標(biāo)準(zhǔn)有效性及強條執(zhí)行的問題。

（4）進(jìn)一步推廣勘測設(shè)計質(zhì)量管理計算機(jī)網(wǎng)絡(luò)信息化系統(tǒng)，促進(jìn)勘測設(shè)計單位質(zhì)量管理的規(guī)范化、程序化，確保勘測設(shè)計過程有章可循，有據(jù)可查，并為全行業(yè)信息化勘測設(shè)計檔案管理系統(tǒng)打好基礎(chǔ)。

4.2 建議

（1）加強水行政主管部門的勘測設(shè)計市場監(jiān)管。主管部門要嚴(yán)格勘測設(shè)計成果審查審批，認(rèn)真開展質(zhì)量評價工作，充分利用市場監(jiān)管、審查審批、質(zhì)量評價等信息，建議水行政主管單位嚴(yán)格設(shè)計成果審查專家的資格管理。

（2）進(jìn)一步加強項目法人代建制。對于一些臨時組建項目法人，技術(shù)管理人員數(shù)量及管理水平不足的項目，建議推行建設(shè)管理“代建制”，以實現(xiàn)建設(shè)項目管理“建設(shè)、監(jiān)管、使用”多位一體，相應(yīng)制約腐敗，提高建設(shè)管理的專業(yè)化、規(guī)范化，促進(jìn)勘測設(shè)計市場健康發(fā)展。

T V 223

：A

：1672-2469（2015）09-0081-05

10.3969/j.i s s n.1672-2469.2015.09.026

張 強（1984年—），男，工程師。