水利枢纽有限元分析.doc

1、目录第一部分3设 计 说 明 书31 兴旺沟水利枢纽工程基本资料41.1 流域概况41.2.3 洪水44、地基重度2700N/m351.3.2 淤沙的力学指标51、淤沙高度23m52、内摩擦角1853、浮容重9.5kN/m351.3.3 水库特性表52 重力坝的工作原理、特点及其荷载62.1 重力坝的工作原理62.2 重力坝的特点62.3 作用在重力坝上的荷载62.4 重力坝有限元建模62.5 重力坝建模参数72.6 设计的主要内容及内容73 坝型选择与枢纽布置73.1 坝型选择7拟定三种方案：分别选择土石坝、混凝土重力坝、拱坝三种类型进行比较。7表4-1三种坝型比较表73.2 枢纽布置84

2、有限元分析方法与ANSYS简介95 坝体的设计及有限元分析计算105.1 挡水坝段设计105.2 挡水坝段空间有限元分析105.3 厂房坝段设计185.4 厂房坝段有限元分析186 结论24第二部分25设 计 计 算 书251 挡水坝段的设计26a）进入ANSYS求解器SOLUTION492. 厂房坝段的设计53谢 辞66第一部分设 计 说 明 书1 兴旺沟水利枢纽工程基本资料1.1 流域概况乐昌峡水利枢纽位于韶关市乐昌境内、北江支流武江乐昌峡河段内，坝址位于兴旺沟火车站附近，下距乐昌市约14km，韶关市81.4km，坝址以上集水面积4988km2。武江是北江流域的一级支流。发源于湖南省临武县

3、三峰岭，流经湖南省的临武县、宜章县、郴县、桂阳、汝城等五县，于乐昌市老坪石上游约3km流入广东省的乐昌市，经乳源、曲江，与浈江在韶关市沙洲尾汇合注入北江。武江全河长260km，流域面积7097km2。1.2 水文1.2.1 气象 武江流域位于五岭山脉以南，属东亚季风气候区。冬半年受东北季风控制，气候寒冷略干燥，夏半年受西南和东南季风控制，气候炎热多雨。本流域具有山地气候特征。大暴雨发生的频率前汛期（46月）占67%69%，后汛期（79月）占31%33%。根据乐昌市气象站1959年2004年资料统计，主要气象特性如下：气温：多年平均气温19.5，极端最高气温41.0（2003年7月23日），极端

4、最低气温-4.6（1967年1月17日）。降雨量：多年平均降雨量1488mm，最大日降雨量220 mm（1973年6月28日）。湿度：多年平均相对湿度80，最大月平均相对湿度90（1975年5月）。风向、风速：多年平均最大风速14.8m/s，最大风速22m/s（1970年8月11日），最多风向为NNW。 蒸发量：多年平均蒸发量1039mm。1.2.2 水文基本资料 根据工程位置、任务、流域地区洪水组成及已有的水文站点和资料，确定武江的坪石（二）、犁市水文站，浈江的长坝（浈湾）水文站和北江的韶关水位站四站为设计基本站。资料系列均至2006年。 1.2.3 洪水根据北江上游洪水组成特点，乐昌峡水库

5、设计洪水考虑甲、乙两种洪水组成。 甲种洪水：韶关与乐昌峡水库同频率，乐昌峡韶关湾头区间相应，以武水为主的洪水组成。乙种洪水：韶关与乐昌峡韶关湾头区间同频率乐昌峡水库相应，以浈水为主的洪水组成。乐昌峡水利枢纽兴旺沟坝址采用以武水为主的甲种洪水组成时还同时考虑乐昌峡坝址乐昌市小区间的甲、乙两种洪水组成。典型年选择：甲种洪水：1994年6月（全流域大水）、2006年7月（武江大水）。乙种洪水：1976年6月（浈江大水）。1.2.4 泥沙 犁市站多年平均含沙量为0.173kg/m3，多年平均侵蚀模数为177.47t/km2，经计算坝址多年平均输沙量为 90.34万t。1.3 坝基岩石及沙砾石的物理学性

6、质资料1.3.1 坝基岩石及砂砾石的物理力学性质1、泊松比 0.22、摩擦系数 0.653、抗剪强度系数 1.04、地基重度2700N/m31.3.2 淤沙的力学指标1、淤沙高度23m2、内摩擦角183、浮容重9.5kN/m3 1.3.3 水库特性表 （1）、水库水位 1）正常蓄水位 154.50m2）设计洪水位（P=1%） 158.42m3）校核洪水位(P=0.1%) 159.18m4) 汛期运行水位 144.50m5）极限死水位 139.50m（2）、电站下游尾水位 1）设计洪水位 116.58m 2）校核洪水位 117.30m 3）正常尾水位 113.87m 4）最低计算尾水位 112.

7、00m（3）、溢洪道下游水位 1）设计洪水位 113.87m2）校核洪水位 116.62m2 重力坝的工作原理、特点及其荷载2.1 重力坝的工作原理 重力坝是用浆砌石或混凝土材料修筑而成的挡水建筑物。一般做成上游面近似垂直的三角形断面，在坝体和地基接触面间产生抗剪强度或摩擦力来抵抗水库的水推力，以达到稳定的要求；同时依靠坝体自重产生的压应力来抵消由于水压力所引起的坝体上游侧面的拉应力，以满足坝身强度的要求。为了适应地基变形、温度变化和混凝土的浇筑能力，沿坝轴线用横缝将坝体分隔成若干个独立工作的坝段。2.2 重力坝的特点（1） 结构作用明确，设计方法简单，安全可靠。（2） 重力坝对地形、地质条件

8、适应性强。任何形状的河谷都可以修建重力坝，因为坝体作用于地基面上的压应力不高，所以对地质条件的要求也较拱坝低，甚至在土基上也可以修建高度不大的重力坝。（3）枢纽泄洪问题容易解决。重力坝可以做成溢流的，也可以在坝身不同高度设置泄水孔，一般不需另设溢洪道或泄水隧洞，枢纽布置紧凑。（4）便于施工导流。在施工期可以利用坝体导流，一般不需要另设导流隧洞。（5）施工方便。大体积混凝土可以采用机械化施工，在放样、立模和混凝土浇筑方面都比较简单，并且补强、修复、维护或扩建也比较方便。（6）重力坝是大体积混凝土建筑物，施工时混凝土的水化发热和散热、硬化收缩将引起坝体内温度和收缩应力，可能使坝体产生裂缝。因此，在

9、浇筑混凝土时，需要有严格的温度控制措施。2.3 作用在重力坝上的荷载 作用在重力坝上的荷载主要有：坝体及坝基上永久设备的自重，上下游坝面上的静水压力，扬压力、泥沙压力、及地震荷载等。2.4 重力坝有限元建模有限元求解坝体和基岩位移和应力应变响应时，关键是整体刚度矩阵K。对于一般的水库大坝，其刚度矩阵应该由坝体和基岩等单元刚度矩阵组合而成，即K=K1+K2。其中K1表示坝体单元的总体刚度矩阵，既受坝体自身刚度的影响，又受基础约束的影响；K2表示基岩单元的总刚度，主要取决于所考虑的范围大小及边界约束条件。数值模拟计算只能在有限的区域内进行，为了减小计算误差，必须选取合适的计算范围。在实际建模过程中

10、，建议地基选取的范围以2倍的坝高为宜，根据坝体剖面尺寸不一，也可以适当调整范围，由实际情况而定。2.5 重力坝建模参数1 坝体：混凝土的弹性模量E=2.4E10Pa，密度为2400Kg/m3，泊松比为=0.167，抗拉强度Fc=1.96E6Pa2 岩基：弹性模量E=2.7E10Pa，密度为2700Kg/m3，泊松比为=0.2。2.6 设计的主要内容及内容1、确定枢纽布置方案，绘出上游或下游立视图。2、主要建筑物设计（1）当水坝坝段设计：确定当水坝的实用剖面及轮廓尺寸，并进行稳定、应力分析。绘出最大坝高处当水坝横剖面图。（2）底孔坝坝段设计：确定底孔坝的基本剖面及轮廓尺寸，进行稳定、应力计算。绘

11、出最大坝高处底孔坝横剖面图。3、 坝体细部构造设计：包括坝体混凝土分区、分缝廊道系统及坝身排水系统的设计。3 坝型选择与枢纽布置3.1 坝型选择拟定三种方案：分别选择土石坝、混凝土重力坝、拱坝三种类型进行比较。 表4-1三种坝型比较表类型条件土石坝混凝土重力坝拱坝筑坝材料坝址附近要有土石料场坝址附近要有砂石，砾料场坝址附近要有砂石，砾料场地形条件对地形适应性强对地形，地质条件适应性强。拱坝的理想地形是左右两岸对称，岸坡平顺无突变地质条件对地基地质条件要求较低对地基要求较低坝基处的岩基要求坚硬，完整，质地均匀，稳定气候条件适应各种不同气候对外界气候条件的适应性好温度变化对坝体应力影响显著泄洪要求

12、泄洪建筑物应具有足够的泄洪能力，坝体在洪水位以上要有足够的安全超高，以防洪水满顶，造成坝的失事。枢纽泄洪问题容易解决。枢纽泄洪问题容易解决抗震要求抗震性好抗震性好抗震性不好根据基本资料显示，坝址下游有足量的混凝土用砂，石料场，砾石料场储量丰富。武江水量丰富，水质矿化度较低，一般对混凝土不具有侵蚀作用，唯河水含泥沙量较大，需进行沉淀处理，且洪水流量大。故本枢纽要求泄洪能力强.而修混凝土重力坝的优点是对地形、地质条件适应性强,对地形和地质条件要求较拱坝低;枢纽泄洪问题容易解决,重力坝可以做成溢流的,也可在坝内设不同高程泄水孔,不需另设溢洪道或泄水隧洞,枢纽布置紧凑;便于施工导流;重力坝剖面尺寸大,

13、因而抵抗洪水漫顶、渗漏、地震和战争破坏的能力都比土石坝强;施工方便;结构作用明确,重力坝沿坝轴线用横缝分成若干坝段,各坝段独立工作,稳定和应力计算都比其他坝型简单.该枢纽的地形、地质和气象条件等资料都显示,适于修混凝土重力坝.拱坝:由于拱坝剖面较薄,坝体几何形状复杂,对于筑坝材料强度,抗渗性和施工质量等要求都比重力坝严格.地形条件是决定拱坝结构形式,工程布置的主要因素.因其理想地形是狭窄河谷,与本枢纽地形不符,故本枢纽不宜修拱坝.土石坝:在洪水流量较大的河流上,土石坝工程的导流,泄洪问题比混凝土坝难以解决.再根据该枢纽资料内容,坝址地质条件较好,建混凝土坝具有较大优越性,故本枢纽也不适合修土石

14、坝.综上选择该水利枢纽为混凝土实体重力坝较为合理。3.2 枢纽布置根据坝址坝型比较选择意见，结合地形、地质条件，针对选定的重力坝坝型拟定下列可能的枢纽组合方案，进行分析比较。1右岸至左岸依次为挡水坝段、厂房坝段、溢流坝段、挡水坝段2右岸至左岸依次为溢流坝段、厂房坝段、挡水坝段3右岸至左岸依次为挡水坝段、溢流坝段、厂房坝段通过上述枢纽方案各建筑物的布置设计、分析比较，将各枢纽方案特点分述如下：方案1：布置合理，溢流坝段泄流的尾水直接进入河道，消能容易解决。方案2：右岸布置溢流坝段尾水不宜处理，需要额外修建尾水渠段增加投资。方案3：左岸直接布置厂房，左岸岸坡需要开挖较多，增加了投资。由此可见，各枢

15、纽方案建筑工程总投资排列顺序为方案(1)方案(2)方案(3)，综合考虑地形地质、水文、施工、运行管理等条件，选择方案，即布置从右岸依次为83.1m的挡水坝段，52.5m的挡厂房坝段，75m的溢流坝段，最左岸为163.5m的挡水坝段。4 有限元分析方法与ANSYS简介4.1 有限元ANSYS简介有限元法（Finite Element Method/FEM）又称有限单元法或有限元素法，是结构分析的一种数值法，自从1960年柯劳夫（R.W.clough）提出“有限元法”的名称以来，有限元法的研究工作蓬勃发展，在工程领域中作用日益增强，已成为分析连续体力学问题最新颖和最有效的工具。有限元法的实质就是把

16、具有无限个自由度的连续体，理想化为只有有限个自由度的一组有限个且按一定方式联结在一起的组合域。各单元之间彼此相连接的点称为节点，从研究有限单元的力学特性着手，最后得到一组以节点位移作为未知量的矩阵形式表达的方程，求解出这些未知量，就可以得到整个求解域上的近似解。有限元法采用矩阵表达形式，便于编制计算机程序，以求解大型线性方程组。利用计算机的便利大大提高了计算速度，使许多复杂的工程问题迎刃而解。4.2 有限元ANSYS分析过程ANSYS有限元求解问题的基本过程主要包括：前置处理、有限元求解、后置处理三部分。介绍如下：（1） 前置处理：建立有限元模型所需要输入的资料，如节点、坐标资料、元素内节点排

17、列次序、材料特性等。之后将建好的模型划分网格，是使整个模型系统离散为有限个单元。（2） 有限元求解 通过对模型外部作用力的加载，利用能量最低原理和函数值定理换成一组线性联立方程组，如元素刚度矩阵计算K、系统外力向量的组合F、线形代数方程KU=F等，使每组作用力都转换在作用在单元体上的荷载，然后进行计算处理。（3） 后置处理 将解题部分所得的解答如：变位、应力、反力等资料，通过图形接口以各种不同表示方式把等位移图、等应力图等显示出来5 坝体的设计及有限元分析计算5.1 挡水坝段设计5.1.1 轮廓尺寸确定 设计时考虑各影响因素：荷载、地基、运用、材料、施工条件等影响，拟定几个方案比较，选出最优方

18、案。重力坝承受的主要荷载是水压力和自重，控制剖面尺寸的主要指标是稳定和强度。因此作用于上游面的水压力及基础扬压力呈三角形分布，所以重力坝的基本剖面是三角形。 坝体高度为81米，坝顶宽度一般取坝高的8%10%，考虑交通要求,经过综合分析,坝顶宽度取8米。上游坝坡坡率n=0，坝底宽为64.8米。符合上游边坡:1:01:0.2;下游边坡:1:0.61:0.80所允许的范围.5.1.2 廊道系统的设计帷幕灌浆需要在坝体浇筑到一定高度后进行，以便利用混凝土压重提高灌浆压力，保证灌浆质量。为此，需要在坝踵附近距上游坝面0.050.1倍作用水头、且不小于45m处设置灌浆廊道。廊道断面多为城门洞形，宽度和高度

19、应能满足灌浆作业的要求，一般宽为2.53m，高位34m，地面距基岩面不宜小于1.5倍廊道宽度。本设计廊道距上游面5.0m，距坝基4.5m，廊道宽3m，高3.5m均满足要求。5.2 挡水坝段空间有限元分析5.2.1 有限元建模基本剖面混凝土重力坝建模的基本参数（1）大坝：弹性模量E=2.4e10Pa，泊松比为=0.167，密度为2400kg/m3，抗拉强度为ft=1.96e6Pa，抗压强度fc=10e6Pa。（2）基岩：弹性模量E=2.7e10Pa，泊松比为=0.2，密度为2700kg/m3 图5-1 挡水坝平面轮廓5.2.2 基本组合坝体荷载图、变形图，应力图 图5-2 基本组合坝体加载图挡水

20、坝最大荷载为884559pa。图5-3 基本组合坝体变形图坝体最大位移为0.019727m。 图5-4 基本组合坝体第一主应力分析图坝踵处出现拉应力，拉应力最大值为0.689453MPa，拉应力区深度为2.7m，不超过廊道中心线，满足规定要求。图5-5 基本组合坝体第三主应力分析图在坝址处出现最大压应力，最大应力值为4.7 MPa，小于混凝土抗压强度值9.6 MPa，满足强度要求。5.2.3 特殊组合（一）坝体荷载图、变形图、应力图 图5-6 特殊组合（一）加载图 图5-7 特殊组合（一）变形图坝体最大位移为0.019841m。图5-8 特殊组合（一）坝体第一主应力分析图在坝踵处出现拉应力，拉

21、应力最大值为0.716MPa，拉应力区深度为2.46m，不超过廊道中心线，满足规定要求。图5-9 特殊组合（一）坝体第三主应力分析图在坝址处出现最大压应力，最大应力值为4.63MPa，小于混凝土抗压强度值9.6 MPa，满足强度要求。5.2.4特殊组合（二）坝体荷载图、变形图、应力图图5-10 特殊组合（二）坝体加载图图5-11 特殊组合（二）坝体沉降量图坝体最大位移为1.68mm。 图5-12 特殊组合（二）坝体第一主应力分析图在坝踵处出现拉应力，拉应力最大值为0.36MPa，拉应力区深度为1.78m，不超过廊道中心线，满足规定要求。 图5-13 特殊组合（二）坝体第三主应力分析图在坝址处出

22、现最大压应力，最大应力值为0.32 MPa，小于混凝土抗压强度值9.6 MPa，满足强度要求。坝体抗滑稳定分析：抗剪强度公式： 基本组合 =1.061.05 满足要求特殊组合(一) =1.141.00 满足要求特殊组合(二) =1.176 1.00 满足要求按抗剪强度公式：基本组合 =4.78 3.0 满足要求特殊组合（一）=2.62.5 满足要求特殊组合（二） =5.28 2.3 满足要求可见挡水坝段满足抗滑稳定要求坝踵坝址在三种工况下符合应力条件，且变形在安全范围内，故设计合格。5.3 厂房坝段设计5.3.1 轮廓尺寸确定 设计时考虑各影响因素：荷载、地基、运用、材料、施工条件等影响，拟定

23、几个方案比较，选出最优方案。重力坝承受的主要荷载是水压力和自重，控制剖面尺寸的主要指标是稳定和强度。因此作用于上游面的水压力及基础扬压力呈三角形分布，所以重力坝的基本剖面是三角形。 坝体高度为76米，坝顶宽度一般取坝高的8%10%，考虑交通要求,经过综合分析,坝顶宽度取10.5米。上游坝坡坡率n=0，坝底宽为60.8米。符合上游边坡:1:01:0.2;下游边坡:1:0.61:0.80所允许的范围.5.3.2 管道的设计 本阶段初选3台40MW混流式水轮机机组,根据已建工程，设计进口段高程为127米，为8米7米的矩形，管径为5米，管道出口高程为89.8米。5.4 厂房坝段有限元分析5.4.1 有

24、限元建模基本剖面混凝土重力坝建模的基本参数（1）大坝：弹性模量E=2.4e10Pa，泊松比为=0.167，密度为2400kg/m3，抗拉强度为ft=1.96e6Pa，抗压强度fc=10e6Pa。（2）基岩：弹性模量E=2.7e10Pa，泊松比为=0.2，密度为2700kg/m3 图5-14 厂房坝建模图5.4.4基本组合厂房坝段稳定分析 图5-15基本组合厂房加载图图5-16基本组合厂房坝段变形图坝体最大位移为17.6mm。图5-17基本组合厂房坝段第一主应力分析图在坝踵处出现拉应力，拉应力最大值为3.4MPa，拉应力区深度为3.78m，不超过廊道中心线，满足规定要求。 图5-18基本组合厂房

25、坝段第三主应力分析图结论：按抗剪强度公式：基本组合 =1.052 1.05 满足稳定要求可见挡水坝段满足抗滑稳定要求。坝踵坝址在该工况下符合应力条件，且指定方向位移变形在安全范围内，故设计合格。5.4.5 特殊组合（一）厂房坝段稳定分析图5-19特殊组合（一）加载图 图5-20特殊组合（一）厂房坝段变形图坝体最大位移为17.6mm。图5-21特殊组合（一）厂房坝段第一主应力分析图在坝踵处出现拉应力，拉应力最大值为3.4MPa，拉应力区深度为4.3m，不超过廊道中心线，满足规定要求。图5-22特殊组合（一）厂房坝段第三主应力分析图结论：按抗剪强度公式：特殊组合(一) =1.02 1.00 满足稳

26、定要求可见挡水坝段满足抗滑稳定要求。坝踵坝址在该工况下符合应力条件，且变形在安全范围内，故设计合格。6 结论 这次兴旺沟枢纽工程有限元分析，通过建模模拟出挡水坝段和厂房坝段，对其进行强度、稳定及抗震分析，参照水工建筑物标准，最后得出以下结论： （1）挡水坝段和厂房坝段应力分析符合规范要求，挡水坝段坝踵的拉应力深度为1.3m，小于坝底宽度的7%，满足要求；溢流坝段的拉应力区深度为1.25m，同样也小于坝底宽度的7%，虽然设计遵循材料力学的思想，在坝底不可以出现拉应力区，但是从有限元方法出发，只要出现的拉应力区不超过坝底宽度的7%或坝踵到廊道帷幕灌浆中心线的长度就符合要求。（2）挡水坝段和厂房坝段

27、符合抗滑稳定要求。挡水坝段在基本组合，特殊组合（一），特殊组合( 二)时的抗滑稳定安全系数值分别为1.06, 1.14, 1.176,分别大于规范要求值1.05,1.0，1.0，满足要求；值分别为4.78， 2.6， 5.28，分别大于规范要求值 3.0， 2.5 ，2.3，满足要求。厂房坝段在基本组合，特殊组（一）时的抗滑稳定安全系数值分别为1.052 ，1.02，分别大于规范要求的1.05， 1.0，满足要求。（3）挡水坝段抗震分析中最大拉应力在规范要求的允许之内，当出现设计情况的地震强度时，坝体不会产生破坏，满足要求。 （4）此次设计虽然结果满足水工设计的各种标准，但是其中也存在一些不足

28、，例如，坝体截面尺寸相对有些偏大，设计有些保守，安全储备有些大，从经济方面考虑还有很大的优化空间，需要在以后的设计中不断进行改进。 总体而言，这次兴旺沟水利枢纽有限元分析运用了功能强大的ANSYS软件，给坝体的分析计算带来了方便。在不断地运用该软件的过程中对一些操作有了更深入的理解，再加上老师的悉心辅导，终于圆满地完成了此次设计。 第二部分设 计 计 算 书1 挡水坝段的设计1.1 剖面轮廓及尺寸1.1.1 坝顶高程的确定 坝顶高程应高于校核洪水位，坝顶上游防浪墙顶的高程应高于波浪顶高程。防浪强顶至设计洪水位或校核洪水位的高差，可按下式计算 =1.24 式中：波浪爬高高度，； 波浪中心线至静水

29、位的高度，；波长，；坝前水深，； 计算风速，是指水面以上10m处10min的风速平均值，水库为正常蓄水位或设计洪水位时，宜采用相应季节50年重现期的最大风速，校核洪水位时，宜采用相应洪水期最大风速的多年平均值； 库面的波浪吹程，； 安全超高m， (1)、设计洪水位情况： D=2km，V=22m/s=1m =1.24=1.24m =10.4m=0.3m=1.24+0.3+0.5=2.04m因此设计蓄水位时防浪墙顶高程=158.42+2.04=160.46m(2)、校核洪水位：D=2km，V=14.8m/s=0.61m =1.24=0.7564m =7.0m=0.17m=0.7564+0.17+0

30、.4=1.33m因此校核洪水位时防浪墙顶高程=159.18+1.33=160.51m取、两种情况的最大值为160.51m，坝顶高程=161m，防浪墙高0.5m。1.1.2实用剖面设计坝底高程为80m，坝高H=161-80=81m，上游坝坡坡率取n=0，下游坝坡坡率取m=0.8，坝底宽为64.8m，坝顶宽度取8m。 图1-1 挡水坝剖面设计1.2 有限元建模1.2.1 单元类型首先将坝体简化建立二维模型，所以单元类型采用PLANE42。在划分单元之前应将材料属性附给坝体和地基。以便求解的准确性。PLANE42是2 维实体结构单元，用于建立 2 维实体结构模型。本单元既可用作平面单元 (平面应力或

31、平面应变)，也可以用作轴对称单元。本单元有 4 个节点，每个节点有 2 个自由度，分别为 x 和 y 方向的平移。本单元具有塑性、蠕变、辐射膨胀、应力刚度、大变形以及大应变的能力。并有一个选项可以支持额外的位移形状。ANSYS的SOLID65单元是专为混凝土、岩石等抗压能力远大于抗拉能力的非均匀材料开发的单元。它可以模拟混凝土中的加强钢筋，以及材料的拉裂和压溃现象。它是在三维8节点等参元SOLID45的基础上，增加了针对于混凝土的性能参数和组合式钢筋模型。SOLID65单元最多可以定义3种不同的加固材料，即此单元允许同时拥有四种不同的材料。混凝土材料具有开裂，压碎，塑性变形和蠕变的能力;加强材

32、料则只能承受拉压，不能承受剪切力。图1-2 PLANE42 单元 图1-3 SOLID65几何模型图1.2.3关于平面应变与平面应力问题：平面应力模型 当结构中任何位置点的应力状态（六个应力分量）中某方向上应力分量为零，其他两个方向上应力不为零，这时三维空间结构处于平面应力状态。在ANSYS中，模型需要在总体坐标系的XOY平面内建立，这时它们具有以下行为特征：1）Z方向上的应力为零，但存在应变。2）Z方向几何尺寸远远小于X和Y方向尺寸。3）只承受XY平面内载荷。4）只在XY平面内出现位移，Z方向没有位移发生。平面应变模型当结构中任何位置点的应变状态（六个应力分量）中某方向上应力分量为零，其他两

33、个方向上应变不为零，这时三维空间结构处于平面应变状态。当一个方向的尺寸远远小于其他两个方向的尺寸，并且垂直于大尺寸方向的横截面积总是保持不变形的结构，如大坝坝体可以近似认为符合该要求。在ANSYS中，总是将几何尺寸很大的方向指定为总体坐标系的Z方向，模型需要在总体坐标系的XOY平面建立模型，这时它们具有以下行为特征：1）Z方向的应变为零，但存在应力。2）Z方向几何尺寸远远大于X和Y方向上尺寸。3）只承受XY平面内载荷。4）位移只在XY平面内发生。1.2. 4 建立模型1）分析类型本设计分析类型为结构分析，所以选择分析类型为Structural。2）初始建模混凝土重力坝建模的基本参数（1）大坝：

34、弹性模量E=2.4e10Pa，泊松比为=0.167，密度为2400kg/m3，抗拉强度为fc=1.96ePa，抗压强度ft=22e6Pa。（2）基岩：弹性模量E=2.7e10Pa，泊松比为=0.2，密度为2700kg/m3。（3）根据水工建筑物抗震设计规范，对于重力坝，反应谱代表值为max=2,Tg=0.2,其中表达式为： 注：设计的基本假定：（4）坝体和坝基连续，即坝体和坝基之间紧密联系在一起；（5）坝基和坝体的材料是均匀的；（6）基岩模型采用线弹性本构模型。坝体上下游基础宽度均取坝高的2倍，基础深取2倍坝高建模的方法是先创建关键点，根据关键点创建面，利用布尔运算建立基本的平面模型，然后对该

35、平面模型网格划分，网格划分是坝体于基础的接触面的各个节点是否对应是关键。平面网格划分后对该坝体和基础分别拖拉形成空间模型。 图1-4 挡水坝ANSYS建模 1.3 网格划分1.3.1 网格划分基本知识网格划分是建模中非常重要的一个环节，它将几何模型转化为由节点和单元构成的有限元模型。网格划分的好坏将直接影响到计算结果的准确性和计算精度，甚至会因为网格划分不合理而导致计算不收敛。网格划分主要包括以下3个步骤：1）定义单元属性（单元类型，实常数，材料属性）Main MeauPreprocessorMaterialMaterial Models2）设定网格尺寸控制 Main MeauPreproce

36、ssorMeshing MeshToolSize ControlsLines set3）执行网格划分Main MeauPreprocessorMeshingMeshMeshTool1.3.2 网格划分质量问题网格划分的好坏将直接影响到计算结果的准确性和计算精度，甚至会因为网格划分不合理而导致计算不收敛。划分网格之前一般需要对网格密度进行必要的控制。合理的单元网格密度是获得高精度的结果的保证，准确捕捉场量的分布和梯度变化，如应力、应变的分布与梯度变化，应力集中位置需要相对较密的单元网格，为了捕捉振动的波形效应必须在一个波长距离上划分多个单元，为了捕捉塑性区域的产生位置、扩展过程以及最终塑性区域的

37、大小要求在适当区域划分较密网格等。在实际工程计算中，用户一定要根据自己的分析目的和问题的性质来确定划分网格的密度控制。很多情况下，例如裂纹分析、应力集中等，都要求某个局部的网格具有比较精细的划分，以便获得精确的结果。网格细化的命令在主菜单中的Preprocessor/meshing/modify mesh子菜单中。再者，坝体和基础的接触面的节点是否完全对应也是关键。1.3.3网格划分具体操作单元网格的划分也是有限元建模过程的一个重要的部分，网格划分的好坏与否直接影响到计算的精度，在坝趾和坝踵及廊道处的应力较大，所以这些地方在进行网格划分的时候应该进一步细化，有两种方法可行：1）在为模型按统一比

38、例划分完后，再在现有单元的基础上选择应力较大的部位再选择合适的比例进一步细化。但是采用进一步细化容易对坝体和地基相接触部位的单元节点对应不上，这对耦合时容易将坝体和地基接触初产生刚性面。对求解的精度带来不利影响。一般采用第二种方法。2）在划分面之前先将线不等分化，在应力集中处线的分段长度小，反之则划分长度相对较长。然后再将面网格化。网格划分可以采取先划分线，再划分面。 图1-5 网格划分1.4 加载求解1）自由度约束：岩基作为固定端施加全约束，垂直水流流向的上坝基面施加X方向的约束2）主要荷载有：自重 静水压力和动水压力 扬压力 浪压力 泥沙压力 (6)土压力 (7)地震荷载 总体上分为基本荷

39、载和特殊荷载两种，设计混凝土重力坝时，荷载组合可分为基本组合和特殊组合，基本组合由基本荷载所组成；特殊除相应的基本组合外，还包括一种或几种特殊荷载。挡水坝段设计分三种情况，分别是：1）基本组合：设计洪水位 158.42米:荷载有:自重+静水压力+扬压力+泥沙压力2）特殊组合（1）:校核洪水位159.18米:荷载有:自重+静水压力+扬压力+泥沙压力3）特殊组合（2）:正常蓄水位 154.5米:荷载有:自重+静水压力+扬压力+泥沙压力+地震荷载注: 在这里只考虑对坝体影响较大的荷载，如上下游水压力、泥沙压力、扬压力、自重。在ANSYS中，加载求解也是很重要的一步，合适的加载和边界条件将能够更好的模

40、拟实际情况；而适当的求解过程的控制将直接影响到求解的精度和所花费的机时，甚至收敛与否。施加荷载的主要内容有：1）施加荷载2）求解过程控制在ANSYS程序中，载荷包括边界条件和外部或内部作用力。对于本设计载荷主要以下两类：1）自由度约束：岩基作为固定端施加全约束，垂直水流流向的上坝基面施加X方向的约束。2）主要荷载有：上、下游静水压力；泥沙压力；扬压力。荷载的梯度计算公式： 实体模型加载，选择主菜单Solution/Define loads/Apply中的分类加载子菜单中的命令。荷载的施加主要考虑三种工况，分别是基本组合、特殊组合（一）和特殊组合（二）。数据分别如下：表11单位：M上游水位下游水

41、位泥沙高程坝顶高程坝顶宽坝底宽基本组合158.42116.58103.00161.00864.80特殊组合（一）159.18117.30103.00161.00864.80特殊组合（二）154.50113.87103.00161.00864.80 图1-6 基本组合施加荷载图1-7特殊组合（一）施加荷载图1-8特殊组合（二）施加荷载1.5 应力分析应力的目的是为了检验大坝在施工期和运用期是否满足强度要求，同时也是为研究解决设计和施工中的某些问题，如混凝土标号分区和某些部位的配筋等提供依据。本设计采用弹性理论的有限元法 。1）有限元法计算的应力控制标准： 材料力学法规定坝体内不容许出现主拉应力，

42、而根据均匀弹性体用有限元法计算高重力坝时，上游坝踵都有一个拉应力区,拉应力区宽度宜小于坝底宽度的0.07倍或小于坝踵至帷幕中心线的距离。目前，对用有限元计算得的坝踵拉应力如何考虑，尚无统一规定。因坝体和岩基均非完全弹性体，用有限元法计算出上游坝踵处较大的拉应力，可能只是造成局部的塑性屈服，而并不一定威胁大坝的安全。2）各工况的应力分析：图形显示S1主应力：选择菜单Main MenuGeneral PostprocPlot ResultsNodal Sol弹出Contour Nodal Solution Data对话框，在 Contour Nodal Solution Data选项的左侧列表中选中Stress选项，右侧列表中选择1st principal S1选项，单击OK按纽观察第一主应力结果。图形显示S3主应力：选择菜单Main MenuGeneral PostprocPlot ResultsNodal Sol弹出Contour Nodal Solution Data对话框，在 Contour Nodal Solution Data选项的左侧列表中选中Stress选项，右侧列表中选择3st principal S3选项，单击OK按纽观察第三主应力结果。 图1-9基本组合第一主应力 图1-10基本组合第三主应力 图1-1