基于變形離散元法的重力壩地震開裂分析

楊利福， 常曉林， 周　偉， 程勇剛， 馬　剛

(武漢大學(xué) 水資源與水電工程科學(xué)國家重點實驗室，武漢　430072)

?

基于變形離散元法的重力壩地震開裂分析

楊利福， 常曉林， 周偉， 程勇剛， 馬剛

(武漢大學(xué) 水資源與水電工程科學(xué)國家重點實驗室，武漢430072)

摘要：考慮脆性材料Ⅰ型和Ⅱ型復(fù)合斷裂情況，基于漸進破壞理論建立了混凝土損傷破壞模型，提出一種損傷開裂模型與變形離散單元法耦合的分析法。模擬了Koyna壩在強烈地震作用下的漸進破壞過程，通過與室內(nèi)模型試驗及其他方法所得結(jié)果的對比，驗證了該方法的可靠性和正確性。研究表明：壩體初始裂縫是由地震拉應(yīng)力集中所致，壩體開裂后，裂縫周圍應(yīng)力重新分布，但縫端仍存在拉應(yīng)力集中現(xiàn)象，導(dǎo)致裂縫不斷擴張；裂縫發(fā)展可分4個階段；裂縫貫通后;壩頭與底部劇烈錯動;壩體動能減小;摩擦功急劇增加。

關(guān)鍵詞：重力壩；地震開裂；變形離散元法；分離裂縫模型

近年來在我國西南地震斷裂帶附近地區(qū)已建、在建或設(shè)計了一大批高混凝土重力壩(龍灘、金安橋、觀音巖等)，其中許多樞紐壩高達到了200～300 m級，創(chuàng)造了同類型壩的世界最高紀錄。西南地區(qū)地震烈度多在Ⅷ～Ⅸ度之間(如金安橋水電站設(shè)計地震加速度已達0.399 g)，由于壩高、庫大，一旦遭遇超強震而潰決失事，不僅會造成重大經(jīng)濟損失，而且對下游所形成的次生災(zāi)害將造成難以估量的人民生命財產(chǎn)損失。

目前，研究混凝土重力壩在強震作用下漸進破壞過程的方法主要有模型試驗和數(shù)值計算等方法。朱彤等[1-2]對混凝土仿真材料特性進行了試驗研究，并研究了龍灘水電站碾壓混凝土重力壩在地震載荷作用下的破壞特性；李弋等[3]經(jīng)一系列配合比和力學(xué)性能試驗及動彈性模量超聲波檢測，研制出新型高混凝土壩動力試驗?zāi)Ｐ筒牧希环稌⒌萚4]采用仿真混凝土材料研究了龍開口水電站溢流壩段在不同因素影響下的動力破壞特征和機理；文獻[5]通過室內(nèi)動力模型試驗系統(tǒng)研究了Koyna壩和Pine Flat壩在地震作用下的非線性響應(yīng)；Tinawi等[6]通過一個3.4 m高的混凝土壩模型振動臺試驗和數(shù)值模擬研究了重力壩裂縫擴展和壩體滑移機制。模型試驗在一定程度上能夠反映混凝土重力壩的破壞過程和規(guī)律，但是模型試驗采用不同材料可能導(dǎo)致不同的破壞模式，載荷和邊界條件難以同時滿足相似律要求，試驗結(jié)果也難以與實際情況一致，并且由于試驗觀測手段的限制，模型試驗無法全面揭示重力壩地震開裂破壞機理。

目前，混凝土重力壩地震響應(yīng)數(shù)值分析主流方法還是有限單元法。在傳統(tǒng)有限單元法基礎(chǔ)上，文獻[7]采用分離裂縫模型、文獻[8-9]采用彌散裂縫模型、文獻[10-11]采用擴展有限單元法研究了重力壩在地震過程中的非線性響應(yīng)和失效模式。然而，分離裂縫模型需要預(yù)設(shè)裂紋發(fā)展路徑，計算過程中需要重剖分計算網(wǎng)格；彌散裂縫模型通過本構(gòu)模型反映材料開裂后的性能退化，但要求裂縫單元邊界一致，才能得到較精確的結(jié)果；擴展有限元計算結(jié)果對網(wǎng)格存在依賴性，上述不同的裂縫處理方法所計算的結(jié)果也不一致。也有學(xué)者在有限單元法的基礎(chǔ)上引入裂縫帶模型[12]、非線性斷裂力學(xué)模型[13]、同軸旋轉(zhuǎn)開裂模型[14]、界面單元[15]等研究混凝土重力壩開裂，但是由于有限單元法小位移的限制使得系統(tǒng)由連續(xù)介質(zhì)完全變?yōu)榉沁B續(xù)介質(zhì)(壩體裂縫貫通)后的力學(xué)行為難以模擬。此外，Das等[16]采用光滑粒子流體動力學(xué)方法(Smoothed particle hydrodynamics, SPH)研究了地震載荷作用下的大壩失效過程，但SPH法要求粒子規(guī)則排列，對于大變形和沖擊載荷問題，其計算結(jié)果并不穩(wěn)定。劉君等[17]應(yīng)用非連續(xù)變形分析(Discontinuous Deformation Analysis,DDA) 和有限單元法(Finite Element Method,FEM)耦合方法研究了碾壓混凝土重力壩的動力響應(yīng)，但該模型尚未考慮混凝土開裂破壞；侯艷麗等[18]采用分離式裂縫模型與離散元耦合的方法模擬了Koyna壩在地震作用下的破壞過程，該方法需要在可能產(chǎn)生斷裂的結(jié)構(gòu)部位預(yù)設(shè)虛擬裂縫，但其破壞模式受預(yù)設(shè)裂縫位置的影響。

變形離散單元法允許塊體發(fā)生大變形及旋轉(zhuǎn)，能夠很好地模擬塊體的脫開和滑移，為分析強震作用下重力壩開裂過程及裂縫貫通后的力學(xué)響應(yīng)提供了一個新方法。因此，本文考慮脆性材料Ⅰ型和Ⅱ型復(fù)合斷裂情況，基于材料斷裂過程區(qū)理論建立了混凝土開裂模型，將其嵌入到離散單元法中，提出一種基于變形離散元法的混凝土損傷開裂模型，將該模型應(yīng)用于Koyna混凝土重力壩在強震作用下破壞過程分析，研究了其由連續(xù)介質(zhì)向非連續(xù)系統(tǒng)轉(zhuǎn)換過程以及非連續(xù)介質(zhì)系統(tǒng)的地震響應(yīng)，得到Koyna壩動力漸進破壞過程、破壞形態(tài)以及裂縫貫通后非線性響應(yīng)，驗證了該模型在混凝土重力壩地震破壞分析應(yīng)用的適用性。

1混凝土損傷開裂模型

變形離散單元法將許多塊體組裝在一起，引入開裂準則后，可以視其為一個連續(xù)的塊體，即當(dāng)塊體接觸面應(yīng)力狀態(tài)滿足一定條件時，塊體開裂，連續(xù)介質(zhì)向非連續(xù)介質(zhì)轉(zhuǎn)變，這種方法的本質(zhì)也屬于分離裂縫模型，但是離散單元法允許系統(tǒng)發(fā)生大位移、轉(zhuǎn)動甚至分離，只要引入合適的開裂破壞準則，就能很好的反映連續(xù)介質(zhì)向非連續(xù)介質(zhì)轉(zhuǎn)變的過程。Kazerani[19-20]基于變形離散單元法研究了火成巖、沉積巖、砂巖在巴西圓盤試驗和單軸壓縮試驗下的細觀破壞機理，并將破壞形態(tài)與物理試驗比對，發(fā)現(xiàn)數(shù)值試驗與物理試驗結(jié)果非常吻合，驗證了該方法在脆性材料開裂研究中的有效性。

分離式裂縫模型將脆性材料的帶狀微裂區(qū)簡化為一條分離的裂縫，并假定斷裂過程區(qū)外的材料仍為彈性變形區(qū)。當(dāng)斷裂過程區(qū)的尖端應(yīng)力超過抗拉強度后，過程區(qū)向前發(fā)展，當(dāng)外界所做功等于斷裂能時，材料完全開裂破壞。只考慮Ⅰ型和Ⅱ型復(fù)合斷裂情況，定義等效位移后，則斷裂能為：

(1)

(2)

式中：δeff為等效位移；δs為切向位移,δn為法向位移，以拉為正；E為材料的彈性模量；KΙC、KΙΙC分別為材料Ⅰ型、Ⅱ型斷裂韌度；Gf為斷裂能。

斷裂過程區(qū)理論認為：當(dāng)不考慮塊體內(nèi)部損傷時，界面剛度應(yīng)反映開裂破壞區(qū)域材料的劣化。在峰值強度前，假設(shè)材料處于彈性階段，處于破壞過程區(qū)的接觸面法向和切向剛度應(yīng)該反映局部材料變形特性的劣化，此時，界面上的應(yīng)力與位移關(guān)系是非線性彈性的，即加卸載曲線重合。在峰值強度之后，材料進入損傷階段，界面的承載力降低 ，剛度逐漸下降，當(dāng)剛度?0時，界面完全失效，材料由連續(xù)介質(zhì)轉(zhuǎn)變?yōu)榉沁B續(xù)介質(zhì)(見圖1)。本文采用Mohr-Coulomb準則作為開裂破壞準則，當(dāng)界面的法向應(yīng)力達到抗拉強度后發(fā)生拉伸破壞，當(dāng)界面的切向應(yīng)力超過抗剪強度則相接觸的塊體發(fā)生滑移。采用指數(shù)形式來描述峰值強度前的剛度劣化，則界面的法向和切向應(yīng)力-位移關(guān)系有：

(3)

(4)

式中：kt、ks為法向和切向初始接觸剛度系數(shù)；δct為峰值應(yīng)力對應(yīng)的位移；σn、σs分別為法向和切向應(yīng)力；σt為材料的抗拉強度；δut為材料完全斷裂時的位移；D為損傷因子;kred為殘余剛度系數(shù)；σsmax=knδntan(φ)+c;φ為接觸面內(nèi)摩擦角;c為黏聚力。

圖1　拉應(yīng)力-位移關(guān)系曲線Fig. 1 Stress-displacement curve

因此，斷裂能可以表示為：

(5)

變形離散單元法采用動態(tài)松弛法對控制方程進行顯式時間積分[21]，系統(tǒng)的動能可以表示為：

(6)

當(dāng)界面上的切向應(yīng)力σs大于抗剪強度σsmax時，塊體滑移錯動并產(chǎn)生摩擦功[22]：

(7)

2數(shù)值模型與計算參數(shù)

印度Koyna壩是混凝土壩動力分析的經(jīng)典問題，它是少數(shù)幾個經(jīng)歷強震并具有較完整記錄的受震害的混凝土壩工程，因此本文以Koyna壩作為計算模型，并以同幅實測地震動作為輸入條件。在該地震中，Koyna壩在91.75 m深的水推力和0.474 g與0.312 g的水平與豎向加速度峰值地震作用下，在多個非溢流壩段，于下游壩高66.5 m和上游壩高60.0～66.5 m處發(fā)生了水平裂縫，并在這些裂縫位置發(fā)現(xiàn)了滲漏，說明裂縫已經(jīng)完全貫穿壩體。為了與現(xiàn)有研究成果對比，驗證本文所提方法的可行性，本文模型尺寸與計算參數(shù)與文獻[23]相同。塊體劃分時，在實際壩體開裂區(qū)域周圍(壩高51.5～70.5 m)加密，模型尺寸(見圖2)，計算參數(shù)見表1。

對于普通混凝土，傳遞應(yīng)力=0時，對應(yīng)的位移δut≈0.05 mm，對于高強混凝土δut可達0.08 mm，本文取δut=0.055 mm，混凝土拉應(yīng)力超過材料強度后，其強度損傷因子按照圖3所示規(guī)律演化。本文在計算時，變形塊體內(nèi)部劃分有限差分單元后采用彈性材料模擬，塊體間的接觸特性采用混凝土漸進破壞模型；采用Westergaard附加質(zhì)量考慮動水壓力效應(yīng)；動力計算采用Rayleigh阻尼，阻尼比=0.05，壩體開裂后，壩體振型及固有頻率將改變，為了保證解的穩(wěn)定性，本文采用文獻[21]推薦的方法，即將自振頻率提高2倍而不考慮剛度阻尼；輸入的順河向和豎向的近壩基位置實測加速度地震記錄(見圖4)，動力計算中在壩底輸入經(jīng)積分后的速度波。

圖2　離散元計算塊體模型(單位：m)Fig.2 Discrete element model(Unit: m)

圖3　損傷因子演化曲線Fig.3 Damage development curve

密度/(kg·m-3)彈性模量/GPa泊松比摩擦系數(shù)抗拉強度/MPa粘聚力/MPa2643.031.030.201.202.903.00

圖4　地震加速度記錄曲線Fig.4 Acceleration of Koyna earthquake time histories

圖5　混凝土單軸拉伸下的應(yīng)力-位移曲線Fig.5 Stress-displacement curve of concrete under uniaxial tension

圖6　不同方法計算結(jié)果破壞形態(tài)比較Fig.6 Failure mode comparison withdifferent methods

3數(shù)值模擬

3.1程序驗證

為檢驗所提混凝土損傷開裂模型與變形離散元法耦合分析方法的可靠性及合理性，本文在商業(yè)離散元軟件UDEC平臺上通過fish語言將式(1)～式(4)編制成相應(yīng)程序，采用混凝土單軸拉伸試驗對程序進行了驗證。由圖5可知,單軸拉伸下的應(yīng)力位移曲線數(shù)值解與理論值基本一致，說明本文所提計算方法能夠模擬混凝土開裂，計算結(jié)果可靠。

文獻[5]給出了Koyna大壩振動臺模型試驗得到的強震破壞模式，文獻[14]采用同軸旋轉(zhuǎn)開裂模型、文獻[23]采用塑性損傷模型給出了Koyna大壩失效模式，見圖6。從圖6可知,本文計算的壩體破壞模式與其他數(shù)值方法計算的壩體失效模式基本一致，與模型試驗結(jié)果很接近, 說明采用本文所提方法模擬混凝土重力壩地震開裂破壞是可行的。本文模擬得到的最終破壞形態(tài)呈“Y”字形：下游折坡處裂縫與坡面基本成90°，裂縫向底部并向上游擴展后出現(xiàn)一個拐點，在此之后裂縫水平向上游擴展，最終形成貫穿性裂縫。

3.2開裂過程分析

圖7給出了Koyna地震波作用下的大壩開裂擴展過程。由圖7可知，大壩從初始裂紋產(chǎn)生到裂縫貫通歷時0.23 s，破壞主要分4個階段。① 下游折坡處裂縫傾斜擴展階段：T>4.31 s 時，在壩體下游折坡處首先出現(xiàn)裂縫，裂縫呈約45°(與坡面呈90°)向下部擴展至水平寬度約1/3壩頸寬后，出現(xiàn)一個拐點；② 中部裂縫向下彎曲擴展階段：在裂縫拐點出現(xiàn)之后，裂縫向下彎曲擴展1/3壩頸寬度；③ 底部裂縫水平擴展階段：T>4.39 s 時，裂縫大致呈水平向向上游面擴展1/3壩頸寬度；④ 中部裂縫水平擴展貫通階段：T>4.42 s 時，底部裂縫停止擴展，在上述拐點(裂縫中部)處出現(xiàn)一條向上游擴展的裂縫，該裂縫先斜向上擴展再呈水平向上游擴展，T=4.54 s時，裂縫擴展至上游面，形成上下游貫通的裂縫。

圖7　Koyna壩開裂擴展圖(單位：m)Fig.7 Crack growth of the Koyna dam(Unit: m)

圖8、圖9為地震過程中壩體主應(yīng)力矢量分布圖(線段長度表示大小，黑色線條表示受壓，灰白色線條表示受拉)從圖可知，在開裂前(T=4.31 s)，壩體上游側(cè)主要受壓，小主應(yīng)力(壓主應(yīng)力)平行于豎直向，下游側(cè)受拉，大主應(yīng)力(拉主應(yīng)力)平行于下游坡面，折坡處拉應(yīng)力應(yīng)力集中現(xiàn)象較為明顯，最大拉主應(yīng)力約3.4 MPa；壩體開裂后(T=4.39 s)，縫端拉應(yīng)力集中較大，且垂直于裂縫擴展方向，開裂區(qū)裂縫周圍雖然還存在拉主應(yīng)力，但是該拉應(yīng)力方向平行于裂縫方向；T=4.42 s，下游側(cè)裂縫周圍的拉應(yīng)力逐漸轉(zhuǎn)化為壓應(yīng)力；當(dāng)裂縫貫通后(T=4.54 s)，下游側(cè)裂縫周圍拉應(yīng)力進一步轉(zhuǎn)換為壓應(yīng)力，此時整個壩體應(yīng)力分布以壓為主，最大壓應(yīng)力絕對值達22.0 MPa。對比不同時刻壩體主應(yīng)力矢量圖可知，壩體下游側(cè)在地震作用下，壩體下游折坡處初始裂縫是由拉應(yīng)力集中導(dǎo)致的，壩體開裂后，裂縫周圍應(yīng)力重分布，但是縫端依然存在拉應(yīng)力集中現(xiàn)象，導(dǎo)致裂縫不斷擴展，最終形成貫穿性裂縫。

圖8　不同時間壩體主應(yīng)力矢量圖Fig.8 Principal stress vector at different time

圖9　不同時間壩體局部主應(yīng)力矢量圖Fig.9 Principal stress vector of the part of dam at different time

3.3開裂后壩體地震響應(yīng)

圖10、圖11為開裂面上A點和B點法向和切向相對位移與接觸力時間歷程圖(法向應(yīng)力以壓為正，法向位移張開為正；切向應(yīng)力與位移符號表示方向)。由圖可知，壩體開裂后，位于裂縫處的應(yīng)力及位移作劇烈變化且向下游滑移，說明壩頭與底部劇烈錯動，壩體由連續(xù)介質(zhì)系統(tǒng)向非連續(xù)介質(zhì)系統(tǒng)轉(zhuǎn)化后，系統(tǒng)固有振動頻率增大[21]。A點最大的張開位移達9.2 mm，B點法向最大張開位移達0.486 mm。由于A點法向張開位移較大，開裂后并未馬上閉合，其法向和切向基本不受力；B點開裂后又閉合，法向和切向接觸力方向變化頻繁，表現(xiàn)出強烈的非線性。圖12為壩體系統(tǒng)摩擦功和動能時間歷程曲線，圖13裂縫貫穿后壩頭滑移圖，從圖12和圖13可知，系統(tǒng)動能達到最大值時，壩體裂縫貫通，壩頭滑移量最大值達4.8 mm；裂縫貫通后，壩體動能急劇減小，摩擦功急劇增加，說明在裂縫貫通后，壩頭與底部劇烈錯動，部分動能以摩擦形式耗散，但整個地震過程中未出現(xiàn)壩頭倒塌現(xiàn)象。

圖10　不同位置裂縫的接觸力Fig.10 Contact force of crack at different location

圖11　不同位置裂縫的相對位移Fig.11 Relative displacement of crack at different location

圖12　摩擦功和動能時間歷程Fig.12 Kinetic energy and friction work time history

圖13　不同時刻壩頭滑移圖Fig.13 Slip map of dam head at different time

4結(jié)論

本文基于脆性材料斷裂過程區(qū)理論, 考慮Ⅰ型和Ⅱ型復(fù)合斷裂情況，建立了混凝土損傷破壞模型，提出一種將損傷開裂模型與變形離散元法耦合的分析方法，對Koyna重力壩強震作用下的破壞過程以及壩體開裂后的動力響應(yīng)進行數(shù)值仿真模擬。主要結(jié)論如下：

(1) 本文采用變形離散元法與混凝土損傷開裂模型耦合分析方法模擬的Koyna重力壩強震失效模式與模型試驗及其他數(shù)值方法所得結(jié)果一致，說明該方法用于模擬混凝土重力壩地震開裂破壞是可行的。

(2) 下游折坡處初始裂縫是由于地震拉應(yīng)力集中所致，壩體開裂后，裂縫周圍應(yīng)力重新分布，但縫端仍存在拉應(yīng)力集中現(xiàn)象，導(dǎo)致裂縫不斷擴張，最終形成貫穿性裂縫。

(3) 在強震載荷作用下，Koyna混凝土重力壩的裂縫發(fā)展可以分4個階段：下游折坡處裂縫傾斜擴展階段，中部裂縫向下彎曲擴展階段，底部裂縫水平擴展階段，中部裂縫水平擴展貫通階段。

(4) 變形離散單元法能夠很好模擬壩體裂縫貫通后非連續(xù)介質(zhì)系統(tǒng)的力學(xué)響應(yīng)，裂縫貫通后，系統(tǒng)固有振動頻率增大，下游側(cè)裂縫基本不受力，上游側(cè)裂縫接觸面在法向上只受壓力作用，切向力方向變化頻繁，壩頭與底部劇烈錯動，壩體動能急劇減小，摩擦功急劇增加，但壩頭未倒塌。

參 考 文 獻

[ 1 ] 朱彤, 林皋, 馬恒春. 混凝土仿真材料特性及其應(yīng)用的試驗研究[J]. 水力發(fā)電學(xué)報, 2004, 23(4): 31-37.

ZHU Tong, LIN Gao, MA Heng-chun. The test research on properties and applications of emulation concrete material[J].Journal of Hydroelectric Engineering,2004,23(4):31-37.

[ 2 ] 朱彤, 林皋,馬恒春，等. 龍灘大壩的動力模型破壞試驗研究[J]. 水電站設(shè)計, 1995, 11(2): 48-54.

ZHU Tong, LIN Gao, MA Heng-chun, et al. Failure test research on dynamic model of long tan dam[J]. Design of Hydroelectric Power Station, 1995, 11(2): 48-54.

[ 3 ] 李弋, 張少杰, 徐進, 等. 混凝土高壩動力模型材料的試驗研究[J]. 巖土力學(xué), 2011, 32(3): 757-760.

LI Yi, ZHANG Shao-jie, XU Jin,et al.Experimental research on dynamic model material of high concrete dam[J].Rock and Soil Mechanics, 2011, 32(3): 757-760.

[ 4 ] 范書立, 陳健云, 周晶, 等. 龍開口水電站溢流壩段動力模型破壞試驗研究[J].水利學(xué)報,2008(增刊1):195-199.

FAN Shu-li, CHEN Jian-yun, ZHOU Jing,et al.Experimental research on over fall sect ion dynamic rupture of longkaikou project[J]. Journal of Hydraulic Engineering, 2008 (Sup1):195-199.

[ 5 ] Panel on earthquake engineering for concrete dams, committee on earthquake engineering. Earthquake engineering for concrete dams: design, performance, and research needs[M]. National Academy Press: 1990.

[ 6 ] Tinawi R, Léger P, Leclerc M, et al. Seismic safety of gravity dams: from shake table experiments to numerical analyses[J].Journal of Structural Engineering,2000,126(4):518-529.

[ 7 ] Mirzabozorg H, Ghaemian M. Non-linear behavior of mass concrete in three-dimensional problems using a smeared crack approach[J]. Earthquake Engineering & Structural Dynamics, 2005, 34(3): 247-269.

[ 8 ] 沈懷至, 周元德, 王進廷. 基于彌散裂縫模型的重力壩簡化地震分析[J]. 水利學(xué)報, 2007, 38(10): 1221-1227.

SHEN Huai-zhi, ZHOU Yuan-de, WANG Jin-ting.Simplified earthquake analysis of concrete gravity dams using smeared crack approach[J].Journal of Hydraulic Engineering, 2007, 38(10): 1221-1227.

[ 9 ] Valamanesh V, Estekanchi H E, Vafai A, et al. Application of the endurance time method in seismic analysis of concrete gravity dams[J]. Scientia Iranica, 2011, 18(3): 326-337.

[10] 方修君, 金峰, 王進廷. 基于擴展有限元法的 Koyna 重力壩地震開裂過程模擬[J]. 清華大學(xué)學(xué)報: 自然科學(xué)版, 2008, 48(12): 2065-2069.

FANG Xiu-jun，JIN Feng，WANG Jin-ting． Seismic fracture simulation of the Koyna gravity dam using an extended finite element method[J].Tsinghua University:Sci & Tech,2008,48(12): 2065-2069．

[11] 張社榮, 王高輝, 龐博慧, 等. 基于 XFEM 的強震區(qū)砼重力壩開裂與配筋抗震措施研究[J]. 振動與沖擊, 2013, 32(6): 137-142.

ZHANG She-rong，WANG Gao-hui，PANG Bo-hui, et al. Seismic cracking and reinforcement analysis of concrete gravity dam based on XFEM[J]. Journal of Vibration and Shock, 2013, 32(6): 137-142.

[12] Guanglun W, Pekau O A, Chuhan Z, et al. Seismic fracture analysis of concrete gravity dams based on nonlinear fracture mechanics[J]. Engineering Fracture Mechanics, 2000, 65(1): 67-87.

[13] El-Aidi B, Hall J F. Non-linear earthquake response of concrete gravity dams part 1: Modelling[J]. Earthquake Engineering & Structural Dynamics,1989,18(6):837-851.

[14] Calayir Y, Karaton M. Seismic fracture analysis of concrete gravity dams including dam-reservoir interaction[J]. Computers & Structures, 2005, 83(19): 1595-1606.

[15] 徐海濱, 杜修力, 楊貞軍. 基于預(yù)插黏性界面單元的 Koyna 重力壩強震破壞過程分析[J]. 振動與沖擊, 2014, 33(17): 74-79.

XU Hai-bin，DU Xiu-li，YANG Zhen-jun.Seismic failure analysis of Koyna gravity dam using cohesive interface elements[J]. Journal of Vibration and Shock, 2014, 33(17): 74-79.

[16] Das R, Cleary P W. A mesh-free approach for fracture modelling of gravity dams under earthquake[J].International Journal of Fracture, 2013, 179(1/2): 9-33.

[17] 劉君, 陳健云, 孔憲京, 等. 基于 DDA 和 FEM 耦合方法的碾壓混凝土壩抗震安全性分析[J]. 大連理工大學(xué)學(xué)報, 2004, 43(6): 793-798.

LIU Jun,CHEN Jian-yun, KONG Xian-jing, et al. Seismic stability analysis of rol led compacted concrete dam based on approach of DDA coupled with FEM[J].Journal of Dalian University of Technology, 2004, 43(6): 793-798.

[18] 侯艷麗.砼壩-地基破壞的離散元方法與斷裂力學(xué)的耦合模型研究[D]. 北京:清華大學(xué), 2005.

[19] Kazerani T. Effect of micromechanical parameters of microstructure on compressive and tensile failure process of rock[J].International Journal of Rock Mechanics and Mining Sciences, 2013, 64: 44-55.

[20] Kazerani T, Yang Z Y, Zhao J. A discrete element model forpredicting shear strength and degradation of rock joint by using compressive and tensile test data[J]. Rock Mechanics and Rock Engineering, 2012, 45(5): 695-709.

[21] 王泳嘉, 邢紀波. 離散單元法及其在巖土力學(xué)中的應(yīng)用[M]. 沈陽: 東北工學(xué)院出版社, 1991.

[22] UDEC (Universal Distinct Element Code) user’s manual[M]. Version 4.01, Itasca Consulting Group. Inc., 2008.

[23] Dassult Systems Inc. ABAQUS Standard User’s manual[M]. Version 6.11.1, 2011.

Seismic cracking analysis of a gravity dam based on deformable distinct element method

YANGLi-fu,CHANGXiao-lin,ZHOUWei,CHENGYong-gang,MAGang(State Key Laboratory of Water Resources and Hydropower Engineering Science, Wuhan University, Wuhan 430072, China)

Abstract:Considering brittle materials’ mode-I and mode-II mixed fracture, a concrete damage failure model was established based on the progressive failure theory. An analysis method considering the coupling between the model and the deformable distinct element method was proposed. The progressive failure process of Koyna dam under the action of strong earthquake was simulated. The reliability and correctness of the proposed method was verified by comparing the simulated results with those of the model tests and other methods. The study results indicated that the initial crack of the dam is caused due to tensile stress concentration induced by earthquake; the dam cracking causes the redistribution of stress but there is a tension stress concentration phenomenon around the crack tip to cause crack expanding; the crack growth can be divided into four phases; after fracture penetration, severe diastrophism happens between the head and bottom of the dam, the dam’s kinetic energy decreases but the friction work dramatically increases.

Key words:gravity dam; seismic fracture; deformable distinct element method; discrete crack model

中圖分類號:TU457

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

DOI:10.13465/j.cnki.jvs.2016.07.008

通信作者周偉 男，教授，1975年生

收稿日期：2015-01-04修改稿收到日期：2015-03-06

基金項目：國家自然科學(xué)基金資助項目(51179139)

第一作者 楊利福 男，博士生，1987年生