混凝土重力壩的三維有限元分析

本文格式為Word版，下載可任意編輯——混凝土重力壩的三維有限元分析

摘要：重力壩的主要安好評(píng)價(jià)指標(biāo)是應(yīng)力分析和抗滑穩(wěn)定分析，目前有限元法是舉行重力壩應(yīng)力分析的主要方法。由于其能夠較好的模擬壩體和壩基接觸、可求解大量繁雜的工程問題、計(jì)算精度較高等眾多優(yōu)點(diǎn)，故其在工程計(jì)算中具有廣泛的應(yīng)用前景。利用大型有限元分析軟件Ansys對(duì)某重力壩工程舉行三維有限元分析，并得出有益結(jié)論，夢想對(duì)實(shí)際工程起到參考性作用。

重力壩；Ansys；三維有限元；應(yīng)力分析

前言

重力壩是水利工程創(chuàng)辦中較為常見的壩型，在重力壩的設(shè)計(jì)過程中應(yīng)力分析和穩(wěn)定分析是最為重要的兩個(gè)片面。材料力學(xué)法是重力壩應(yīng)力分析的傳統(tǒng)分析方法，但由于其做了大量假定和簡化，故計(jì)算結(jié)果誤差較大，而有限元法相比材料力學(xué)法具有較高的計(jì)算精度，故有限元法在壩工計(jì)算中具有廣泛的應(yīng)用前景。鑒此，本文采用大型有限元分析軟件ANSYS對(duì)某混凝土重力壩工程舉行應(yīng)力應(yīng)變分析，并分析溢流壩段和非溢流壩段的應(yīng)力和變形規(guī)律，結(jié)果提出相應(yīng)的改善措施。

1工程概況

某水庫位于遼陽市弓長嶺區(qū)南沙村東，地處太子河干流中部，其操縱流域面積為6175平方千米，兼具有防洪、澆灌、發(fā)電、供水多功能的綜合利用的水利樞紐。擋水建筑物為混凝土重力壩，最大壩高45.8m，壩頂寬度為5m，壩頂長度為477m，分成33個(gè)壩段，其中有13個(gè)位溢流壩段，溢流壩段寬度均為13m，重力壩上游面為鉛直，下游面坡比為1：0.75。建基面高程為55m，正常蓄水位為96.60m，對(duì)應(yīng)的下游水位為60m；設(shè)計(jì)洪水位為96.81m，對(duì)應(yīng)的下游水位為65.20m；校核洪水位為100.48m，對(duì)應(yīng)的下游水位為69m。其非溢流壩和溢流壩段剖面圖如圖1、圖2所示。

2重力壩三維有限元數(shù)值模擬

2.1有限元模型

舉行三維有限元數(shù)值模擬需要體驗(yàn)前處理、求解、后處理三個(gè)步驟；其中前處理中的建立重力壩的三維模型是最繁雜的過程。為了簡化建模，此處疏忽了閘墩和工作橋?qū)χ亓螒?yīng)力分析的影響。為了真實(shí)的反映重力壩各點(diǎn)的應(yīng)力狀態(tài)，在建模的過程中需要考慮地基變形對(duì)壩體應(yīng)力和變形的影響。分別沿壩踵和壩趾向上、下游取1倍壩高距離的地基舉行分析，地基的深度也取為1倍的壩高。地基和壩體的計(jì)算參數(shù)如表1所示。

本次重力壩模型是通過AUTOCAD軟件舉行三維建模，建立好物理模型之后將其導(dǎo)入到ANSYS中，通過表1來定義其材料屬性，并采用Solid185單元。定義好材料屬性和單元類型后，采用自由分網(wǎng)對(duì)重力壩模型舉行網(wǎng)格剖分，其三維有限元網(wǎng)格如圖3所示。

2.2載荷施加

對(duì)有限元網(wǎng)格模型舉行施加載荷包括兩片面：一是邊界處的位移載荷，二是壩體及壩基承受的外力載荷。對(duì)上述模型左右岸垂直的壩基面施加水平方向的位移載荷，令其位移為零；對(duì)根基的地面施加豎直方向的位移載荷，令其位移為零；對(duì)上下游的垂直地基面施加沿水流方向的位移載荷，令其位移為零。

重力壩在正常運(yùn)用期主要承受自重、揚(yáng)壓力、靜水壓力三種荷載，為了簡化計(jì)算，此處只要考靜水壓力和自重荷載。按表1來設(shè)置壩體和壩基材料的物理特性參數(shù)。由于在建壩之前壩基沉降已經(jīng)穩(wěn)定，為了反映壩體的真實(shí)沉降，故此處需做提升，即將地基設(shè)置為無質(zhì)量地基。但需留神采用無質(zhì)量地基那么不能反映地基內(nèi)的真實(shí)應(yīng)力，由于此處只研究壩體的變形和應(yīng)力，故可作此處理。同時(shí)分別在壩體的上游面和下游面以及與水接觸的河谷面施加水荷載，此處計(jì)算采用正常蓄水位、設(shè)計(jì)洪水位、校核洪水位三種工況，其各工況如表2所示。

2.3成果分析

對(duì)上述模型舉行求解便可得到各結(jié)點(diǎn)的位移，從而得到壩體內(nèi)各點(diǎn)的位移、應(yīng)變、應(yīng)力，結(jié)果便可生成大壩的應(yīng)力云圖和位移云圖。此處以校核水位工況下的應(yīng)力云圖和位移云圖為例。在舉行重力壩的有限元分析我們最關(guān)切第一主應(yīng)力、第三主應(yīng)力、豎直方向的應(yīng)力、豎直方向的位移、水平方向的位移。通過計(jì)算其應(yīng)力、位移云圖如圖4-圖9所示。

通過上述應(yīng)力云圖可知豎直方向最大拉應(yīng)力和第一主應(yīng)力均發(fā)生在靠近河谷兩側(cè)的壩踵處，而最大壓應(yīng)力發(fā)生在非溢流壩段的下游壩趾處；其中最大豎向拉應(yīng)力為0.166MPa，最大拉應(yīng)力為0.915MPa，最大壓應(yīng)力1.13MPa；故在校核水位處境下大壩得志強(qiáng)度要求。

通過上述位移云圖可知壩體豎向最大沉降位移發(fā)生在溢流壩段和非溢流壩段靠下游的中上片面區(qū)域，沿水流方向最大位移和總位移均發(fā)生在靠近河谷最近的非溢流壩段頂部區(qū)域；其中最大豎向沉降位移（即下表中的（uy）amx）為1.02mm，沿水流方向最大位移（即下表中的（uz）amx）為1.4mm。

分別計(jì)算三種運(yùn)用工況，求得各工況下的最大應(yīng)力和位移如表3所示。

由計(jì)算結(jié)果可知在三種運(yùn)用工況下重力壩均可得志其強(qiáng)度安好。由分析可知當(dāng)上游水位越高，最大應(yīng)力和位移并不確定越大，由于壩體的最大應(yīng)力和位移與下游水位也有關(guān)系。故在舉行重力壩的應(yīng)力應(yīng)變分析時(shí)不能僅僅只看校核水位處境下的應(yīng)力和變形，而要計(jì)算各種工況并綜合其成果來舉行重力壩的設(shè)計(jì)。

3終止語

（1）重力壩在運(yùn)用期，其局部壩踵處輕易展現(xiàn)拉應(yīng)力集中，為防止壩體混凝土被拉裂，需在拉應(yīng)力較大區(qū)域配置受拉鋼筋。最大壓應(yīng)力通常發(fā)生在非溢流壩段的壩趾處，一般處境下最大壓應(yīng)力遠(yuǎn)小于混凝土的抗壓強(qiáng)度，但為了安好起見，在壓應(yīng)力較大區(qū)域宜采用高強(qiáng)混凝土。（2