重力壩抗震性能有限元分析＊

王金龍

(濰坊學院,山東 濰坊 261061)

重力壩抗震性能有限元分析＊

王金龍

(濰坊學院,山東 濰坊 261061)

建立了重力壩結構的有限元分析模型。采用靜力理論和動力理論,對壩體進行了靜力分析和動力分析,得到了壩體變形圖、合位移等值線圖、有效應力等值線圖、壩體振型圖。通過分析可以得到地震載荷對壩體產生的影響,并為壩體的設計提供理論依據。

重力壩;抗震性;有限元;地震

地震是一種常見而又危害極大的自然現象。在臨近地震中心的地方,地面先是上下震動,然后開始水平震動,隨著地震中心距離的增大,兩種震動逐漸減弱,而且上下震動不如水平震動明顯[1]。

1 地震作用理論

1.1 靜力理論

靜力理論創始于意大利,發展于日本。1900年代,日本學者大森房吉、佐野利器等對其發展做出了重要貢獻。靜力理論認為,結構物是剛體,其最大加速度等于地震動最大加速度,所以結構所受地震力為:

式中,W為物體重量;k為水平地震系數,日本稱為工程震度或工程列度。

此方法比較簡單,且用這種方法設計的建筑物大多經受了一般地震的考驗。但該方法完全忽略了結構本身的動力特性的影響,因此只適合于低矮的,剛性較大的建筑,如路基、擋土墻和重力式橋臺等[2]。

1.2 動力理論

動力理論主要包括反應譜理論、隨機振動等非線性振動理論,以及現在數值仿真常用的時程分析法。對于反應譜理論,它考慮了地震時地面的運動特性與結構物自身的動力特性,是以單質點體系在實際地震作用下的反應為基礎來分析結構反應的方法,是當前工程設計應用最為廣泛的抗震設計方法。

圖1 重力壩斷面結構

反應譜僅能給出結構各振型反應的最大值,而丟失了與最大值和振型組合的有關的重要信息,使得難以正確的進行各振型最大值的組合。隨著電子計算機的廣泛應用,將實際地震的加速度輸入結構計算模型,直接分析結構的地震反應已成為可能。可直接獲得地震過程中結構節點各時刻的位移、速度和加速度,從而計算各時刻豎向地震作用和構件的地震內力。這種方法稱為時間歷程分析法,簡稱時程分析法。

2 模型建立

本例選取應用非常廣泛的重力壩斷面結構如圖1所示。壩高120米,壩底寬為76米,壩頂為10米,上游壩面坡度和下游壩面坡度如圖1所示[3]。

因為重力壩結構比較簡單,垂直于長度方向的斷面結構受力分布情況也基本相同,并且大壩的縱向長度遠大于其橫斷面,因此大壩抗震性能分析選用單位斷面進行平面應變分析是可行的。大壩抗震性能分析的計算條件如下:

(1)假設大壩的基礎是嵌入到基巖中,地基是剛性的。

(2)大壩采用的材料參數為:彈性模量E=35GPa,泊松比v=0.2,容重γ=25KN/m3。

(3)計算分析大壩水位為120米。

(4)水的質量密度為1000kg/m3。

(5)大壩設防地震烈度為8,水平方向地震加速度值為0.2g。

因為大壩是縱向很長的實體,故計算模型可以簡化為平面應變問題。選擇PLANE42單元,壩體模型的單元網格如圖2所示。對壩體模型底部施加位移約束,并施加重力加速度和水壓力載荷。

圖2 壩體模型單元網格

圖3 壩體變形圖

3 求解結果

3.1 靜力分析求解結果

壩體變形圖如圖3所示;壩體合位移云圖如圖4所示;壩體有效應力云圖如圖5所示。

圖4 壩體合位移等值線圖

圖5 壩體有效應力等值線圖

3.2 模態分析求解結果

壩體第1階振型如圖6所示;壩體第2階振型如圖7所示。

圖6 壩體第1階振型圖

圖7 壩體第2階振型圖

3.3 模態合并求解后壩體位移云圖

因為響應譜分析是在頻域內進行的,對于結構動力特性依賴于頻率而變化,因此在模態分析后要進行模態合并求解,才能得到壩體結構真實的總體效應。合并模態求解后,得到壩體在各階頻率的真實位移云圖。壩體第1階合位移云圖如圖8所示;壩體第2階合位移云圖如圖9所示;壩體第1階有效應力云圖如圖10所示;壩體第2階有效應力云圖如圖11所示。

圖8 壩體第1階合位移等值線圖

圖9 壩體第2階合位移等值線圖

圖10 壩體第1階有效應力等值線圖

圖11 壩體第2階有效應力等值線圖

4 結束語

對壩體進行了靜力分析和動力分析,通過對重力壩靜力有限元分析,可以知道壩體在靜力載荷作用下的位移場和應力場,從而可以了解壩體的安全性能。通過對重力壩抗震性能有限元分析,可以了解大壩在地震載荷作用下的動力響應特性,從而可以評價大壩在地震載荷作用下的安全性能。

[1]尚曉江,邱峰,趙海峰,等.ANSYS結構有限元高級分析方法與范例應用[M].北京:中國水利水電出版社,2006.

[2]闞前華,譚長建,張娟,等.ANSYS高級工程應用實例分析與二次開發[M].北京:電子工業出版社,2006.

[3]賴永標,胡仁喜,黃書珍.ANSYS11.0土木工程有限元分析典型范例[M].北京:電子工業出版社,2007.

(責任編輯:肖恩忠)

Fin ite-Element Analysis of Earthquake Resistant for Gravity Dam

WANG Jin-long
(Weifang University,Weifang 261061,China)

Finite-element analysismodelof gravity dam isestablished.The dam body was analyzed by using static theories and dynamic theories,by simulation,distortion,disp lacement,effective stress and vibration types of the dam body are obtained.Finally,good results are obtained based on analysis and calculation of a gravity dam.

gravity dam,quake-p roof,finite element,earthquake

2010-09-26

王金龍(1972-),男,山東臨朐人,濰坊學院建筑工程學院副教授,博士。

TU18 文獻標識碼:A 文章編號:1671-4288(2010)06-0113-03