水利枢纽有限元分析毕业设计

2 水利枢纽有限元分析毕业设计 第一部分 设 计 说 明 书 1 兴旺沟 水利枢纽工程基本资料 流域概况 乐昌峡 水利枢纽位于韶关市乐昌境内、北江支流武江 乐昌峡 河段内，坝址位于 兴旺沟 火车站附近，下距乐昌市约 14关市 址以上集水面积 4988 武江是北江流域的一级支流。发源于湖南省临武县三峰岭，流经湖南省的临武县、宜章县、郴县、桂阳、汝城等五县，于乐昌市老坪石上游约 3入广东省的乐昌市，经乳源、曲江，与浈江在韶关市沙洲尾汇合注入北江。武江全河长 260域面积 7097 水文 气象 武江流域位于五岭山脉以南，属东亚季风气候区。冬半年受东北季风控制，气候寒冷略干燥，夏半年受西南和东南季风控制，气候炎热多雨。本流域具有山地气候特征。大暴雨发生的频率前汛期（ 4 6月）占 67% 69%，后汛期（ 7 9月）占 31% 33%。 根据乐昌 市 气象站 1959年 2004年资料统计，主要气象特性如下： 气温：多年平均气温 ，极端最高气温 （ 2003年 7月 23日），极端最低气温 1967年 1月 17 日）。 降雨量：多年平均降雨量 1488大日降雨量 220 1973年 6月 28 日）。 3 湿度：多年平均相对湿度 80，最大月平均相对湿度 90（ 1975年 5月）。 风向、风速：多年平均最大风速 s，最大风速 22m/s（ 1970年 8 月 11 日），最多风向为 蒸发量：多年平均蒸发量 1039 水文基本资料 根据工程位置、任务、流域地区洪水组成及已有的水文站点和资料，确定武江的坪石（二）、犁市水文站，浈江的长坝（浈湾）水文站和北江的韶关水位站四站为设计基本站。资料系列均至2006年。 洪水 根据北江上游洪水组成特点， 乐昌峡 水库设计洪水考虑甲、乙两种洪水组成。 甲种洪水：韶关与 乐昌峡 水库同频率， 乐昌峡 韶关湾头区间相应，以武水为主的洪水组成。 乙种洪水：韶关与 乐昌峡 韶关湾头区间同频率 乐昌峡 水库相应，以浈水为主的洪水组成。 乐昌峡 水利枢纽 兴旺沟 坝址采用以武水为主的甲种洪水组成时还同时考虑 乐昌峡 坝址乐昌市小区间的甲、乙两种洪水组成。 典型年选择 ： 甲种洪水： 1994 年 6 月（全流域大水）、 2006 年 7 月（武江大水）。 乙种洪水： 1976 年 6 月（浈江大水）。 泥沙 犁市站多年平均含沙量为 年平均侵蚀模数为 计算坝址多年平均输沙量为 t。 坝基岩石及沙砾石的物理学性质资料 坝基岩石及砂砾石的物理力学性质 1、泊松比 、摩擦系数 、抗剪强度系数 、地基重度 2700N/ 淤沙的力学指标 1、淤沙高度 23m 2、内摩擦角 18 4 3、 浮容重 水库特性表 （ 1）、水库水位 1）正常蓄水 位 ）设计洪水位（ P=1%） ）校核洪水位 (P= ) 汛期运行水位 ）极限死水位 2）、电站下游尾水位 1）设计洪水位 ）校核洪水位 ）正常尾水位 ）最低计算尾水位 3）、溢洪道下游水位 1）设计洪水位 ）校核洪水位 重力坝的工作原理、特点及其荷载 重力坝的工作原理 重力坝是用浆砌石或混凝土材料修筑而成的挡水建筑物。一般做成上游面近似垂直的三角形断面，在坝体和 地基接触面间产生抗剪强度或摩擦力来抵抗水库的水推力，以达到稳定的要求；同时依靠坝体自重产生的压应力来抵消由于水压力所引起的坝体上游侧面的拉应力，以满足坝身强度的要求。为了适应地基变形、温度变化和混凝土的浇筑能力，沿坝轴线用横缝将坝体分隔成若干个独立工作的坝段。 重力坝的特点 （ 1） 结构作用明确，设计方法简单，安全可靠。 （ 2） 重力坝对地形、地质条件适应性强。任何形状的河谷都可以修建重力坝，因为坝体作用于地基面上的压应力不高，所以对地质条件的要求也较拱坝低，甚至在土基上也可以修建高度不大的重力坝。 （ 3）枢纽泄 洪问题容易解决。重力坝可以做成溢流的，也可以在坝身不同高度设置泄水孔， 5 一般不需另设溢洪道或泄水隧洞，枢纽布置紧凑。 （ 4）便于施工导流。在施工期可以利用坝体导流，一般不需要另设导流隧洞。 （ 5）施工方便。大体积混凝土可以采用机械化施工，在放样、立模和混凝土浇筑方面都比较简单，并且补强、修复、维护或扩建也比较方便。 （ 6）重力坝是大体积混凝土建筑物，施工时混凝土的水化发热和散热、硬化收缩将引起坝体内温度和收缩应力，可能使坝体产生裂缝。因此，在浇筑混凝土时，需要有严格的温度控制措施。 作用在重力坝上的 荷载 作用在重力坝上的荷载主要有：坝体及坝基上永久设备的自重，上下游坝面上的静水压力，扬压力、泥沙压力、及地震荷载等。 力坝有限元建模 有限元求解坝体和基岩位移和应力应变响应时，关键是整体刚度矩阵 K。对于一般的水库大坝，其刚度矩阵应该由坝体和基岩等单元刚度矩阵组合而成，即 K=2。其中 受坝体自身刚度的影响，又受基础约束的影响； 示基岩单元的总刚度，主要取决于所考虑的范围大小及边界约束条件。数值模拟计算只能在有限的区域内进行，为了减小计算误差，必须选 取合适的计算范围。 在实际建模过程中，建议地基选取的范围以 2倍的坝高为宜，根据坝体剖面尺寸不一，也可以适当调整范围，由实际情况而定。 重力坝建模参数 1 坝体：混凝土的弹性模量 E=度为 2400Kg/松比为 =拉强度 岩基：弹性模量 E=度为 2700Kg/松比为 = 计的主要内容及内容 1、 确定枢纽布置方案，绘出上游或下游立视图。 2、 主要建筑物设计 （ 1） 当水坝坝段设计：确定当水坝的实用剖面及轮 廓尺寸，并进行稳定、应力分析。绘出最大坝高处当水坝横剖面图。 （ 2） 底孔坝坝段设计：确定底孔坝的基本剖面及轮廓尺寸，进行稳定、应力计算。绘出最大坝高处底孔坝横剖面图。 3、 坝体细部构造设计：包括坝体混凝土分区、分缝廊道系统及坝身排水系统的设计。 6 3 坝型选择与枢纽布置 坝型选择 拟定三种方案：分别选择土石坝、混凝土重力坝、拱坝三种类型进行比较。 表 4种坝型比较表 类型 条件 土石坝 混凝土重力坝 拱坝 筑坝材料 坝址附近要有土石料场 坝址附近要有砂石，砾料场 坝址附近要有砂石，砾料场 地形条件 对地形适应性强 对地形，地质条件适应性强。 拱坝的理想地形是左右两岸对称，岸坡平顺无突变 地质条件 对地基地质条件要求较低 对地基要求较低 坝基处的岩基要求坚硬，完整，质地均匀，稳定 气候条件 适应各种不同气候 对外界气候条件的适应性好 温度变化对坝体应力影响显著 泄洪要求 泄洪建筑物应具有足够的泄洪能力，坝体在洪水位以上要有足够的安全超高，以防洪水满顶，造成坝的失事。 枢纽泄洪问题容易解决。 枢纽泄洪问题容易解决 抗震要求 抗震性好 抗震性好 抗震性不好 根据基本资料显示，坝址下游有足 量的混凝土用砂，石料场，砾石料场储量丰富。武江水量丰富，水质矿化度较低，一般对混凝土不具有侵蚀作用，唯河水含泥沙量较大，需进行沉淀处理，且洪水流量大。故本枢纽要求泄洪能力强 质条件适应性强 ,对地形和地质条件要求较拱坝低 ;枢纽泄洪问题容易解决 ,重力坝可以做成溢流的 ,也可在坝内设不同高程泄水孔 ,不需另设溢洪道或泄水隧洞 ,枢纽布置紧凑 ;便于施工导流 ;重力坝剖面尺寸大 ,因而抵抗洪水漫顶、渗漏、地震和战争破坏的能力都比土石坝强 ;施工方便 ;结构作用明确 ,重力坝沿坝轴线用横缝分成若干坝段 ,各坝段独立工作 ,稳定和应力计算都比其他坝型简单 纽的地形、地质和气象条件等资料都显示 ,适于修混凝土重力坝 . 拱坝 :由于拱坝剖面较薄 ,坝体几何形状复杂 ,对于筑坝材料强度 ,抗渗性和施工质量等要求都比重力坝严格 工程布置的主要因素 与本枢纽地形不符 ,故本枢纽不宜修拱坝 . 土石坝 :在洪水流量较大的河流上 ,土石坝工程的导流 ,泄洪问题比混凝土坝难以解决 坝址地质条件较好 ,建混凝土坝具有较大优越性 ,故本枢纽也不适合修土石坝 . 综上选择该水利枢 纽为混凝土实体重力坝较为合理。 枢纽布置 根据坝址坝型比较选择意见，结合地形、地质条件，针对选定的重力坝坝型拟定下列可能的枢纽组合方案，进行分析比较。 1右岸至左岸依次为挡水坝段、厂房坝段、溢流坝段、挡水坝段 2右岸至左岸依次为溢流坝段、厂房坝段、挡水坝段 3右岸至左岸依次为挡水坝段、溢流坝段、厂房坝段 通过上述枢纽方案各建筑物的布置设计、分析比较，将各枢纽方案特点分述如下： 方案 1：布置合理，溢流坝段泄流的尾水直接进入河道，消能容易解决 。 方案 2：右岸布置溢流坝段尾水不宜处理，需要额外修建尾水渠 段增加投资。 方案 3：左岸直接布置厂房，左岸岸坡需要开挖较多，增加了投资。 由此可见，各枢纽方案建筑工程总投资排列顺序为方案 (1)方案 (2)方案 (3)， 综合考虑地形地质、水文、施工、运行管理等条件，选择方案，即 布置从右岸依次为 75左岸为 4 有限元分析方法与 介 有限元 介 有限元法（ 又称有限单元法或有限元素法，是结构分析的一种数值法 ，自从 1960 年柯劳夫（ 出“有限元法”的名称以来，有限元法的研究工作蓬勃发展，在工程领域中作用日益增强，已成为分析连续体力学问题最新颖和最有效的工具。有限元法的实质就是把具有无限个自由度的连续体，理想化为只有有限个自由度的一组有限个且按一定方式联结在一起的组合域。各单元之间彼此相连接的点 8 称为节点，从研究有限单元的力学特性着手，最后得到一组以节点位移作为未知量的矩阵形式表达的方程，求解出这些未知量，就可以得到整个求解域上的近似解。有限元法采用矩阵表达形式，便于编制计算机程序，以求解大型 线性方程组。利用计算机的便利大大提高了计算速度，使许多复杂的工程问题迎刃而解。 有限元 析过程 限元求解问题的基本过程主要包括：前置处理、有限元求解、后置处理三部分。介绍如下： （ 1） 前置处理： 建立有限元模型所需要输入的资料，如节点、坐标资料、元素内节点排列次序、材料特性等。之后将建好的模型划分网格，是使整个模型系统离散为有限个单元。 （ 2） 有限元求解 通过对模型外部作用力的加载，利用能量最低原理和函数值定理换成一组线性联立方程组，如 元素刚度矩阵计算 K、系统外力 向量的组合 F、线形代数方程 K U=F等，使每组作用力都转换在作用在单元体上的荷载，然后进行计算处理。 （ 3） 后置处理 将解题部分所得的解答如：变位、应力、反力等资料，通过图形接口以各种不同表示方式把等位移图、等应力图等显示出来 5 坝体 的设计及有限元分析计算 挡水坝段设计 轮廓尺寸确定 设计时考虑各影响因素：荷载、地基、运用、材料、施工条件等影响，拟定几个方案比较，选出最优方案。 重力坝承受的主要荷载是水压力和自重，控制剖面尺寸的主要指标是稳定和强 度。因此作用于上游面的水压力及基础扬压力呈三角形分布，所以重力坝的基本剖面是三角形。 坝体高度 为 81 米，坝顶宽度一般取坝高的 8%10%，考虑交通要求 ,经过综合分析 ,坝顶宽度取 8 米。上游坝坡坡率 n=0，坝底宽为 。符合上游边坡 :1:0 1:游边坡 :1:1: 廊道系统的设计 9 帷幕灌浆需要在坝体浇筑到一定高度后进行，以便利用混凝土压重提高灌浆压力，保证灌浆质量。为此，需要在坝踵附近距上游坝面 作用水头、且不小于 4 5m 处设置灌 浆廊道。廊道断面多为城门洞形，宽度和高度应能满足灌浆作业的要求，一般宽为 3m，高位 34m，地面距基岩面不宜小于 廊道宽度。本设计廊道距上游面 坝基 道宽 3m，高 满足要求。 挡水坝段 空间 有限元分析 有限元建模 基本剖面 混凝土重力坝建模的基本参数 （ 1）大坝：弹性模量 E=松比为 =度为 2400kg/拉强度为压强度 0 （ 2）基岩：弹性模量 E=松比为 =度为 2700kg/ 5挡水坝 平面轮廓 本组合坝体荷载图、变形图，应力图 10 图 5基本组合坝体加载图 挡水坝最大荷载为 884559 图 5基本组合坝体变形 图 坝体最大位移为 11 图 5基本组合坝体第一主应力分析图 坝踵处 出现拉应力，拉应力最大值为 应力区深度为 超过廊道中心线，满足规定要求。 图 5基本组合坝体第三主应力分析图 在坝址处出现最大压应力，最大应力值为 于混凝土抗压强度值 足强度要求。 12 特殊组合（一）坝体 荷载图、变形图、应力图 图 5特殊组合（一）加载图 图 5特殊组合（一） 变形 图 坝体最大位移为 13 图 5特殊组 合（一）坝体第一主应力 分析 图 在坝踵处出现拉应力 ，拉应力最大值为 应力区深度为 超过廊道中心线，满足规定要求。 图 5特殊组合（一）坝体第三主应力分析图 在坝址处出现最大压应力，最大应力值为 于混凝土抗压强度值 14 足强度要求。 殊组合（二）坝体 荷载图、变形图、应力图 图 5特殊组合（二） 坝体 加载图 图 5特殊组合（二） 坝体 沉降量图 坝体最大位移为 15 图 5特殊组合（二 ） 坝体第一主应力 分析 图 在坝踵处出现拉应力 ，拉应力最大值为 应力区深度为 超过廊道中心线，满足规定要求。 图 5特殊组合（二）坝体 第三主应力 分析 图 在坝址处出现最大压应力，最大应力值为 于混凝土抗压强度值 足强度要求。 16 坝体抗滑稳定分析： 抗剪强度公式： () 基本组合 满足要求 特殊组合 (一 ) 满足要求 特殊组合 (二 ) 满足要求 按抗剪强度公式： ( ) U c 基本组合 满足要求 特殊组合（一） 满足要求 特殊组合（二） = 满足要求 可见挡水坝段满足抗滑稳定要求 坝踵坝址在三 种工况下符合应力条件，且变形在安全范围内，故设计合格。 厂房坝段设计 轮廓尺寸确定 设计时考虑各影响因素：荷载、地基、运用、材料、施工条件等影响，拟定几个方案比较，选出最优方案。 重力坝承受的主要荷载是水压力和自重，控制剖面尺寸的主要指标是稳定和强度。因此作 用于上游面的水压力及基础扬压力呈三角形分布，所以重力坝的基本剖面是三角形。 坝体高度为 76米，坝顶宽度一般取坝高的 8%10%，考虑交通要求 ,经过综合分析 ,坝顶宽度取 。上游坝坡坡率 n=0，坝底宽为 合上游边坡 :1:0 1:游边坡 :1: 管道的设计 本阶段 初 选 3台 40根据已建工程，设计进口段高程为 127米，为 8米 7米的 矩形，管径为 5米，管道出口高程为 。 房 坝段有限元分析 限元建模基本剖面 17 混凝土重力坝建模的基本参数 （ 1）大坝：弹性模量 E=松比为 =度为 2400kg/拉强度为压强度 0 （ 2）基岩：弹性模量 E=松比为 =度为 2700kg/ 5厂房 坝建模 图 本组合厂房 坝段 稳定分析 18 图 5本组合厂房加载图 图 5本组合厂房 坝段 变形 图 坝体最大位移为 19 图 5本组合厂房坝段 第一 主应力分析图 在坝踵处出现拉应力 ，拉应力最大值为 应力区深度为 超过廊道中心线，满足规定要求。 图 5本组合厂房 坝段 第三主应力分析图 结论： 20 按抗剪强度公式：基本组合 满 足稳定要求 可见挡水坝段满足抗滑稳定要求 。 坝踵坝址在 该 工况下符合应力条件，且指定方向位移变形在安全范围内，故设计合格。 殊组合（一）厂房坝段稳定分析 图 5殊组合（一）加载图 21 图 5殊组合（一）厂房坝段 变形 图 坝体最大位移为 图 5殊组合（一） 厂房坝段 第一 主应力分析图 在坝踵处出现拉应力 ，拉应力最大值为 应力区深度为 超过廊道中心线，满足规定要求。 22 图 5殊组合（一） 厂房坝段 第三 主应力分析图 结论： 按 抗剪强度公式：特殊组合 (一 ) 满足稳定要求 可见挡水坝段满足抗滑稳定要求 。 坝踵坝址在 该 工况下符合应力条件，且变形在安全范围内，故设计合格。 6 结论 这次 兴旺沟 枢纽工程有限元分析，通过建模模拟出挡水坝段和厂房 坝段，对其进行强度、稳 定及抗震分析 ，参照水工建筑物标准，最后得出以下结论： （ 1） 挡水坝段和厂房 坝段应力分析符合规范要求，挡水坝段坝踵的拉应力深度为 于坝底宽度的 7%，满足要求； 溢流坝段的拉应力区深度为 样也小于坝底宽度的 7%，虽然设计遵循材料力学的思想，在坝底不可以出现 拉应力区，但是从有限元方法出发，只要出现的拉应力区不超过坝底宽度的 7%或坝踵到廊道帷幕灌浆中心线的长度就符合要求。 （ 2）挡水坝段和厂房 坝段符合抗滑稳定要求。 挡水坝段在基本组合，特殊组合（一），特殊组合 ( 二 )时的抗滑稳定安全系数 分别为 足要求； 分别为 分别大于规范要求值 23 足要求。厂房 坝段在基本组合，特殊组（一）时的抗滑稳定安全系数 分别为 别大于规范要求 的 满足要求。 （ 3）挡水坝段抗震分析中最大拉应力在规范要求的允许之内，当出现设计情况的地震强度时，坝体不会产生破坏，满足要求。 （ 4）此次设计虽然结果满足水工设计的各种标准，但是其中也存在一些不足，例如，坝体截面尺寸相对有些偏 大，设计有些保守，安全储备有些大，从经济方面考虑还有很大的优化空间，需要在以后的设计中不断进行改进。 总体而言，这次 兴旺沟 水利枢纽有限元分析运用了功能强大的 件，给坝体的分析计算带来了方便 。在不断地运用该软件的过程中对一些操作有了更深入的理解，再加上老师的悉心辅导，终于圆满地完成了此次设计。 第二部分 设 计 计 算 书 24 1 挡水坝段的设计 剖面轮廓及尺寸 坝顶高程的确定 坝 顶高程应高于校核洪水位，坝顶上游防浪墙顶的高程应高于波浪顶高程。防浪强顶至设计洪水位或校核洪水位的高差 h ，可按下式计算 1% h h h 5 / 4 1 / 300 . 0 1 6 6 D 1%h= 0 . 8 2h c 式中： 波浪爬高高度， m ； 25 波浪中心线至静水位的高度， m ； L 波长， m ； H 坝前水深， m ； 0V 计算风速 /指水面以上 10m 处 10风速平均值，水库为正常蓄水位或设计洪水位时，宜采用相应季节 50 年重现期的最大风速，校核洪水位时 ，宜采用相应洪水期最大风速的多年平均值 ； D 库面的波浪吹程 ， 安全超高 m， (1)、设计洪 水位 情况 ： D=2V=22m/s 5 / 4 1 / 300 . 0 1 6 6 D =1m1%h= = . 8 2h c =.0ch m 1% h h h =因此 设计 蓄水位 时 防浪墙 顶高程 =(2)、 校核洪水位： D=2V=s 5 / 4 1 / 300 . 0 1 6 6 D =h= . 8 2h c =26 4.0ch m 1% h h h =此校核洪水位时 防浪墙 顶高 程 = 1 、 2 两种情况的最大值为 坝顶高程 =161m，防浪墙高 用剖面设计 坝底高程为 80m，坝高 H=1611m，上游坝坡坡率取 n=0，下游坝坡坡率取 m=底宽为 顶宽度取 8m。 图 1挡水坝剖面设计 有限元建模 单元类型 首先将坝体简化建立二维模型，所以单元类型采用 划分单元之前应将材料属 27 性附给坝体和地基。以便求解的准确性。 2 维实体结构单元 ， 用于建立 2 维实体结构模型。本单元既可用作平面单元 (平面应力或平面应变 )，也可以用作轴对称单元。本单元有 4 个节点，每个节点有 2 个自由度，分别为 x 和 y 方向的平移。本单元具有塑性、蠕变、辐射膨胀、应力刚度、大变形以及大应变的能力。并有一个选项可以支持额外的位移形状。 元是专为混凝土、岩石等抗压能力远大于抗拉能力的非均匀材料开发的单元。它可以模拟混凝土中的加强钢筋，以及材料的拉裂和压溃现象。它是在三维 8 节点等参元 基础上，增加了针对于混凝土的性能参数和组合式钢筋模型。 元最多可以定义 3 种不同的加固材料，即此单元允许同时拥有四种不同的材料。混凝土 材料具有开裂，压碎，塑性变形和蠕变的能力 ;加强材料则只能承受拉压，不能承受剪切力。 图 1元 图 1何模型图 于平面应变与平面应力问题： 平面应力模型 28 当结构中任何位置点的应力状态（六个应力分量）中某方向上应力分量为零，其他两个方 向上应力不为零，这时三维空间结构处于平面应力状态。在 型需要在总体坐标系的 时它们具有以下行为特征 ： 1） 存在应变 。 2） Z 方向几何尺寸远远小于 X 和 Y 方向尺寸。 3）只承受 面内载荷。 4）只在 面内出现位移， Z 方向没有位移发生。 平面应变模型 当结构中任何位置点的应变状态（六个应力分量）中某方向上应力分量为零，其他两个方向上应变不为零，这时三维空间结构处于平面应变状态。当一个方向的尺寸远远小于其他两个方向的尺寸，并且垂直于大尺寸方向的横截面积 总是保持不变形的结构，如大坝坝体可以近似认为符合该要求。在 是将几何尺寸很大的方向指定为总体坐标系的 型需要在总体坐标系的 时它们具有以下行为特征 ： 1） Z 方向的应变为零，但存在应力。 2） Z 方向几何尺寸远远大于 X 和 Y 方向上尺寸。 3）只承受 面内载荷。 4）位移只在 面内发生。 4 建立模型 1）分析类型 本设计分析类型为结构分析，所以选择分析类型为 2）初始建模 混凝土重力坝建模的基本参数 （ 1）大坝：弹性模量 E=松比为 =度为 2400kg/拉强度为压强度 2 （ 2）基岩：弹性模量 E=松比为 =度为 2700kg/ （ 3）根据水工建筑物抗震设计规范，对于重力坝，反应谱代表值为 ,中表达式为：( 注：设计的基本假定： （ 4）坝体和坝基连续，即坝体和 坝基之间紧密联系在一起； 29 （ 5）坝基和坝体的材料是均匀的； （ 6）基岩模型采用线弹性本构模型。 坝体上下游基础宽度均取坝高的 2倍，基础深取 2倍坝高 建模的方法是先创建关键点，根据关键点创建面，利用布尔运算建立基本的平面模型，然后对该平面模型网格划分，网格划分是坝体于基础的接触面的各个节点是否对应是关键。平面网格划分后对该坝体和基础分别拖拉形成空间模型。 图 1挡水坝 模 网格划分 网格划分基本知识 网格划分是建模中非常重要 的一个环节，它将几何模型转化为由节点和单元构成的有限元模型。网格划分的好坏将直接影响到计算结果的准确性和计算精度，甚至会因为网格划分不合理而导致计算不收敛。 网格划分主要包括以下 3个步骤： 1）定义单元属性（单元类型，实常数，材料属性） ）设定网格尺寸控制 ）执行网格划分 30 格划分质量问题 网格划分的好坏将直接影响到计算结果的准确性和计算精度，甚至会因为网格划分不合理而导致计算不收敛。划分网格之前一般需要对网格密度进行必要的控制。合理的单元网格密度是获得高精度的结果的保证，准确捕捉场量的分布和梯度变化，如应力、应变的分布与梯度变化，应力集中位置需要相对较密的单元网格，为了捕捉振动的波形效应必须在一个波长距离上划分多个单元，为了捕捉塑性区域的产生位置、扩展过程以及最终塑性区域的大小要求在 适当区域划分较密网格等。在实际工程计算中，用户一定要根据自己的分析目的和问题的性质来确定划分网格的密度控制。 很多情况下，例如裂纹分析、应力集中等，都要求某个局部的网格具有比较精细的划分，以便获得精确的结果。网格细化的命令在主菜单中的 菜单中。再者，坝体和基础的接触面的节点是否完全对应也是关键。 格划分具体操作 单元网格的划分也是有限元建模过程的一个重要的部分，网格划分的好坏与否直接影响到计算的精度，在坝趾和坝踵及廊道处的应力较大， 所以这些地方在进行网格划分的时候应该进一步细化，有两种方法可行： 1）在为模型按统一比例划分完后，再在现有单元的基础上选择应力较大的部位再选择合适的比例进一步细化。但是采用进一步细化容易对坝体和地基相接触部位的单元节点对应不上，这对耦合时容易将坝体和地基接触初产生刚性面。对求解的精度带来不利影响。一般采用第二种方法。 2）在划分面之前先将线不等分化，在应力集中处线的分段长度小，反之则划分长度相对较长。然后再将面网格化。 网格划分可以采取先划分线，再划分面。 31 图 1网格划分 载求解 1）自由度约束：岩基作为固定端施加全约束，垂直水流流向的上坝基面施加 X 方向的约束 2）主要荷载有 ： 自重 静水压力和动水压力 扬压力 浪压力 泥沙压力 (6)土压力 (7)地震荷载 总体上分为基本荷载和特殊荷载两种，设计混凝土重力坝时，荷载组合可分为基本组合和特殊组合，基本组合由基本荷载所组成；特殊 除相应的基本组合外，还包括一种或几种特殊荷载。挡水坝段设计分三种情况，分别是： 1）基本组合：设计洪 水位 : 荷载有 :自重 +静水压力 +扬压力 +泥沙压力 2）特殊组合（ 1） :校核洪水位 : 荷载有 :自重 +静水压力 +扬压力 +泥沙压力 3）特殊组合（ 2） :正常蓄水位 : 荷载有 :自重 +静水压力 +扬压力 +泥沙压力 +地震荷载 注 : 在这里只考虑对坝体影响较大的荷载，如上下游水压力、泥沙压力、扬压力、自重。 在 ，加载求解也是很重要的一步， 合适的加载和边界条件将能够更 好的模拟实际 32 情况；而适当的求解过程的控制将直接影响到求解的精度和所花费的机时，甚至收敛与否。 施加荷载的主要内容有： 1）施加荷载 2）求解过程控制 在 序中，载荷包括边界条件和外部或内部作用力。对于本设计载荷主要以下两类： 1）自由度约束：岩基作为固定端施加全约束，垂直水流流向的上坝基面施加 X 方向的约束。 2）主要荷载有：上、下游静水压力；泥沙压力；扬压力。 荷载的梯度计算公式：00( ) ( ) ( )P x P x x x 00( ) ( ) ( )P y P y y y 00( ) ( ) ( )P z P z z z 实体模型加载，选择主菜单 的分类加载子菜单中的命令。 荷载的施加主要考虑三种工况，分别是基本组合、特殊组合（一）和特殊组合（二）。 数据分别如下： 表 1 1 单位： M 上游水位 下游水位 泥沙高程 坝顶高程 坝顶宽 坝底宽 基本组合 殊组合（一） 殊组合（二） 33 图 1基本组合 施加荷载 图 1殊组合（一）施加荷载 34 图 1殊组合（二）施加荷载 应力分析 应力的目的是为了检验大坝在施工期和运用期是否满足强度要求，同时也是为研究解决设计和施工中的某些问题，如混凝土标号分区和某些部位的配筋等提供依据。本设计采用 弹性理论的有限元法 。 1）有限元法计算的应力控制标准： 材料力学法规定坝体内不容许出现主拉应力，而根据均匀弹性体用有限元法计算高重力坝时，上游坝踵都有一个拉应力区 ,拉应力区宽度宜小于坝底宽度的 前，对用有限元计算得的坝踵拉应力如何考虑，尚无统一规定。因坝体和岩基均非完全弹性体，用有限元法计算出上游坝踵处较大的拉应力，可能只是造成局部的塑性屈服，而并不一定威胁大坝的安全。 2）各工况的应力分析： 图形显示 应力：选择菜单 出话框，在 项的左侧列表中选中 侧列表中选择 1st 1 选项，单击 纽观察第一主应力结果。 图形显示 应力：选择菜单 出 35 话框，在 项的左侧列表中选中 侧列表中选择 3st 3 选项，单击 纽观察第三主应力结果。 图 1本组合第一主应力 图 1本组合 第三 主应力 36 图 1本组合 X 方向位移图 图 1本组合 Y 方向位移图 37 图 1本组合 Z 方向位移图 图 1本组合总方向位移图 基本组合 计算量 符号 /单位 数值 第一主应力 坝踵 1（兆帕） 址 38 第三 主应力 坝踵 2（兆帕） 址 方向位移 （） 方向位移 （） 方向位移 （） 方向位移 （） 1特殊组合（ 一 ）第一主应力 39 图 1特殊组合 (一 )第三 主应力 图 1殊组合（一） X 方向位移图 40 图 1殊组合（一） Y 方向位移图 图 1殊组合（一） Z 方向位移图 41 图 1殊组合（一） 总 方向位移图 特殊组合（一） 计算量 符号 /单位 数值 第一主应力 坝踵 1（兆帕） 址 三 主应力 坝踵 2（兆帕） 址 方向位移 （） 方向位移 （） 方向 位移 （） 方向位移 （） 42 图 1特殊组合（二）第一主应力 图 1特殊组合 (二 )第三 主应力 43 图 1殊组合 (二 )X 方向位移图 图 1殊组合 (二 )Y 方向位移图 44 图 1殊组合 (二 )Z 方向位移图 图 1殊组合 (二 )总方向位移图 特殊组合（二） 计算量 符号 /单位 数值 第一主应力 坝踵 1（ 兆帕） 45 坝址 三 主应力 坝踵 2（兆帕） 址 方向位移 （） 方向位移 （） 方向位移 （） 方向位移 （） 滑稳定的计算 抗剪强度公式： () 式中 ： 抗滑稳定安全系数； f 坝体混凝与坝基接触面的抗剪摩擦系数 ，本工程取 U 作用在接触面上的扬压力； W 接触面以上的总铅直力 ； P 接触面以上的总水平力 。 抗滑稳 定安全指标，查教材水工建筑物得基本组合 (1) ， 特殊组合（ 1） =特殊组合（ 2） = 抗剪断强度公式： ( ) U c 式中： 抗滑稳定安全系数； f 坝体混凝与坝基接触面的抗剪 断擦系数，本工程 c 坝体混凝土与坝基接触面的抗剪凝聚力，本工程 取 A 坝基接触面的截面积。 46 抗剪断安全指标，设计规范规定，不分工程级别，基本组合 (1) ， 特殊组合 (1) ， 特殊组合（二）不小于 (1) 有限元抗滑稳定数据的提取 具体步骤如下： 1）利用菜单 T择需要计算的位置。 2） T这里需要提取的结果项是： P 提取结果如下，带入公式即可得到抗剪断强度，抗剪强度。 挡水坝段稳定计算成果表 稳 定 计 算 组合 计算量 符 号 基本组合 特殊组合 (一 ) 特殊组合 (二 ) 垂直力 W(N) 08 09 08 水平力 P (N) 08 09 08 抗剪公式 剪断公式 定指 标分析总结 1）按抗剪强度公式： 基本组合 满足要求 特殊组合 (一 ) 满足稳定要求 特殊组合 (二 ) 满足稳定要求 2）按抗剪断强度公式： 基本组合 满足要求 特殊组合（一） 满足要求 特殊组合（二） =满足要求 可见挡水坝段满足抗滑稳定要求 47 水工建筑物抗震设计 我国受环太平洋地震带的影响，地震活动频繁，历史上多次发生灾害性大地震，全国大部分地区为抗震设防区，近期又处于地震活动上升期，需要重视水工建筑物的抗震设计。 水库蓄水会影响震源活动条件可诱发地震活动，而地震作用是典型的动态作用，在地基随机性运动的影响下，可能使基岩断层活动发生错动 ，沙地层液化，库水对坝产生动水压力，建筑物开裂或倾倒，填土对挡土建筑物产生动土压力、水库库岸崩塌、土石坝坝坡滑动或沉降裂缝等作用效应，而对于这一系列的破坏，对水工建筑物的危害非常大，故需进行抗震设计。在考虑到以上情况和地震作用效应时，降低工程的失效率，增加其可靠度，从而作到设计合理、经济可靠、技术先进的设计要求。 地震的强烈程度，是一次地震活动的规模已其释放的总能量来评价，按里氏震级划分标准，最大震级不超过九级，而我国将地震烈度