基于有限元法的土石坝渗流及稳定研究分析-research and analysis on seepage and stability of earth-rock dam based on finite element method.docxVIP

目录1绪论 11.1研究背景 11.2国内外渗流计算发展概况 31.2.1国外渗流计算发展概况 31.2.2国内渗流计算发展概况 41.3土石坝渗流计算的任务及内容 41.4本文研究的主要内容 51.5本文的创新之处 61.6本章小结 62渗流有限元分析的理论基础及基本方法 72.1概述 72.2渗流的几个基本概念 72.2.1渗流作用力 72.2.2渗透系数 72.3渗流基本定律——达西定律 82.4渗流的基本方程 92.4.1渗流的连续性方程 92.4.2渗流的基本微分方程 102.4.3渗流基本微分方程定解条件 112.5渗流计算方法简介 122.6土石坝渗流的有限元法 132.6.1有限元法概述 132.6.2有限元法的分析过程 132.7本章小结 143基于 ANSYS 的土石坝二维稳定渗流的数值模拟 15ANSYS 软件介绍153.1.1概述 15ANSYS 软件的功能15APDL 语言简介163.2ANSYS 渗流计算原理分析 163.2.1理论基础相似 163.2.2微分方程相似 163.2.3定解条件相似 173.3ANSYS 热分析简介183.3.1ANSYS 热分析介绍183.3.2稳态热分析过程 183.4基于 ANSYS 的土石坝的浸润线求解193.5ANSYS 计算二维稳定渗流算例203.5.1算例概况 203.5.2模型网格划分 203.5.3ANSYS 求解及结果分析213.6本章小结 224边坡稳定性分析 234.1稳定分析的几种方法 234.1.1极限平衡方法 234.1.2数值分析法 244.2有限元强度折减法 244.2.1安全系数 244.2.2破坏标准 254.2.3本构模型 254.3边坡稳定计算 254.3.1典型考核题 254.3.2ANSYS 软件计算及结果分析264.4本章小结 295工程实例分析 305.1工程概况 305.2渗流场数值计算 315.2.1建模并划分单元网格 315.2.2ANSYS 求解及结果分析315.3坝坡稳定分析 335.3.1不考虑渗透力作用下的上游坝坡稳定分析 345.3.2不考虑渗透力作用下的下游坝坡稳定分析 355.4本章小结 376结论与展望 386.1结论 386.2展望 38致谢 39参考文献 40作者简介 43插图和附录清单1. 图 1微分单元体各面上进出流量示意图 92. 图 2渗流计算方法 133. 图 3三角形单元 134. 图 4ANSYS 模块结构示意图 155. 图 5心墙坝横断面图 206. 图 6模型网格划分图 217. 图 7 压力水头云图 218. 图 8 流速分布矢量图 229. 图 9剖面尺寸 2610. 图 10圆弧滑裂面 2611. 图 11模型网格划分图 2812. 图 12求解迭代不收敛过程图 2813. 图 13折减系数为 1.11 时的塑性云图 2914. 图 14心墙坝横断面图 3015. 图 15模型网格划分图 3116. 图 16部分压力水头云图 3217. 图 17压力水头等值线图 3218. 图 18心墙内压力水头等值线图 3319. 图 19下游坝脚处浸润线出口图 3320. 图 20折减系数为 1 时坝体上游水平位移及塑性区云图 3421. 图 21折减系数为 1.65 时坝体上游水平位移及塑性区云图 3422. 图 22折减系数为 1.655 时坝体上游水平位移及塑性区云图 3523. 图 23折减系数为 1.656 时坝体上游塑性区云图 3524. 图 24折减系数为 1 时坝体下游水平位移及塑性区云图 3725. 图 25折减系数为 1.8 时坝体下游水平位移及塑性区云图 3626. 图 26折减系数为 1.808 时坝体下游水平位移及塑性区云图 3627. 图 27折减系数为 1.809 时坝体下游塑性区云图 3728. 表 1世界大坝按坝型分类统计 129. 表 2世界土石坝数量按国家或地区统计表 . 130. 表 3不同种土体的渗透系数参考数值 831. 表 4浸润线位置结果比较 2232. 表 5材料性质 2533. 表 6材料性质 341绪论1.1研究背景土石坝又称当地材料坝，它是主要利用当地的土石材料填筑而成的一种挡水坝 [1]。 土石坝作为历史最悠久的一种坝型，在大多数国家建设数量居于首位。据国际大坝 会议统计，截止到 2003 年底，全世界共建 15m 高度以上的大坝 41413 座，其中土 石坝 33958 座，占 82.7%，见表 1。我国于公元 400 年东晋时期，开始修建荆江大 堤，堤高