洮河海甸峽堆石面板坝在稳定状态和非稳定状态下渗流分析的研究毕业设计.docVIP

1绪 论 1.1本文研究得意义 大坝作为一种挡水构筑物一般都会涉及到渗流稳定性问题,并且它们的破坏主要都是从渗流开始。如若不及时引起重视进行处理,必将危及大坝的整体稳定性。坝体发生渗流破坏的原因,通常都是由于水力梯度引起的。因此，要研究确定水工构筑物的水力梯度的分布,同时必须通过其进行渗流分析,确定浸润自由面，再对坝体进行渗流稳定分析，以此对大坝的安全运行进行作出合理评估。 ，由于作用下会坝体和坝基渗流现象，运行不利的，可是又是不可避免的。根据统计，我国已建成上万座各种类型水库，经过长年累月的运行，很多坝体在不同程度上都存在严重的病害或者隐患，经过广泛分析研究后，分析得出在可能对坝体造成破坏的所有模式中，渗流都在其中扮演着重要的作用。因此，对大坝渗流的研究是大坝安全检测的一项重要内容，深入的分析和研究坝体渗流状况，目前较为理想的一种计算渗流场的方法是有限单元法。在渗流场的基础上对其进行整体稳定性分析,才能得到较为可靠的结果。因此,用有限单元法计算其渗流场，对挡水构筑物安全、维护、治理等都有巨大的社会效益和经济效益。 有限单元法的基本思路和原理早就有人提出，它是一套求解微分方程的数值计算法，与传统的解题法相比，有物理意义直观明确，解题效能高方法本身的适应性强，解题步骤易规范化，电子计算机配合使用后，有限单元法得到了迅速的应用和发展。 法国工程师达西(H.Dracy)，在1986年通过试验研究提出了线性渗透定律，为渗流理论的发展莫定了基础。达西的试验是通过液体在均质沙土中作匀速流动的情况下进行的。这个研究结果已被后来的学者逐渐推广至整个渗流计算领域中去，达西定律成为渗流分析计算中最基本、最重要的公式,它解决了土体的基本渗透计算理论。达西定律研究得出的公式是对均匀渗流平面上计算平均流速以及渗流区域内任意一点的渗流流速。1857年，法国学者杜比推导出杜比公式来研究非均匀渐变渗流分析的规律。1889年，H.E茹可夫斯基研究发现土体中渗流规律分析计算过程符合拉普拉斯方程，首次推导出了渗流的微分方程。但是，相对于许多数学模型和解析解的研究只能适用于均匀的渗透介质和边界比较简单的条件基础上，因此在实际工程的运用中就要受到一定限制。1914年，P. Forehheiness提出利用绘制流网的方法进行图解渗流试算。流网针对分析整个渗流场这个二维平面问题给出了比较可靠全面的图案。绘制出了流网，整个渗流场问题就会得到解决。直到1922年，为解决边界条件和渗流条件比较复杂的渗流问题，H.H·巴甫洛夫斯基研究出了一个有效的求解渗流场的工具—电拟法。所谓电拟法，就是采用电导液模型进行渗流模拟试验，但仍有一些弊端就是无法对非均质各向异性渗透介质的渗流进行模拟实验，当遇到一些比较复杂的地质和复杂的边界条件时就是去了适用性。后来。 数值方法有限差分法、有限单元法，。 边界元法是在贝蒂互换定理经典力学理论及弗雷德霍姆积分方程基础上建立起来的，限差分法法目前，。1960年，克劳夫为了区别有限差分法，而最早提出并采用“有限单元法”这个名词。的研究和应用在国内外对都引起了高度重视，。 随着现代计算机技术的高度发展和应用，有限元法也得到了迅速发展。近几年来，早期从事渗流分析方面的研究学者，采用有限元方法分析具有自由面的无压渗流问题可是往往由于渗流流域的边界具有的复杂性，渗流区域的范围是未知的，这类渗流问题需要通过反复不断的分析，再反复不断的修改自由面的位置，直至求解、迭代逐步确定。其中变网格法和固定网格法是修改渗流自由面最常用到的分析方法。 在渗流计算过程中由于渗流自由面变网格法在求解渗流问题。为突破这个局限提出了新的研究方法。研究出了（2）固定网格在渗流分析的研究其中Galerkin方法1.4.1常用的渗流计算方法 常用的渗流计算方法分为两大类；如图所示。 图1-1渗流计算方法 Figure 1 The methods of seep computing 理论分析法。水力学法和流体力学法是理论分析法两大类。 水力学法 (2)数值方法。边界元法、有限单元法和差分法是目前最常用的三种数值方法。 边界元法：首先将渗流边界元法相对于有限元法而言具有以下优点：便于对半无限或无限渗透介质的处理、渗流奇异和自由面等问题一般说其计算精度也高于有限单元法。边界元法的缺点：当面对三维非均质各向异性的渗透介质问题时，在应用方面尚存在相当大的困难。 有限单元法：它归纳总结了对古典近似计算方法，既吸收了变分计算中选择试探函数的内核，又继承了有限差分法中离散处理的精髓，并对分析区域进行积分，。法在差分法：控制方程中的导数用节点上的函数值的差商代替进行离散，从而建立以网格节点值为未知的方程组。将微分问题转化为问题。常常由于差分网格的规则受到一定的局限；其次，在对各向异性曲线边界进行模拟分析时