科学计算选讲论文.pdf

科学计算选讲论文 论文题目：无奇异边界积分方程在坝基渗流问题 中的应用 学院: 材料科学与工程 姓名: 秦承群 学号: 1012208017 2013-6-4 无奇异边界积分方程在坝基渗流问题中的应用 引言 物理学、力学和工程技术中的许多问题都可归结为求解偏微分方程的初边值 问题。当只有初始条件而没有边界条件时称为初值问题;反之则称为边值问题。 区域内的偏微分方程称为基本方程或控制方程。常把初始条件和边界条件合称为 边值条件。 工程实际问题中，由于物性不均匀、物体几何形状不规则、边界条件较复杂 等原因，能采用解析法并按照边值条件求解偏微分方程的情况少之又少。一般只 能研究它的近似解法——数值方法，来求解近似解。目前，数值解法己成为人类 探索未知科学领域和进行大型工程设计的一种重要方法和手段。其中边界元法是 一种新型的较精确有效的工程数值分析方法，又称边界积分方程-边界元法。边 界元法是在定义域的边界上划分单元，将满足控制方程的函数去逼近边界条件， 具有单元的未知数少，数据准备简单等优点。 近年来，中外边界元法研究者在奇异积分的处理问题上做了大量的工作，发 展了许多种数值处理技术，如解析和半解析计算CPV 积分，数值计算奇异积分， 域外奇点法等。本文介绍了无奇异边界积分方程在坝基渗流问题中的应用。 坝基渗流问题 地下渗流分析是大坝边坡稳定性分析中一项主要内容，是影响水工建筑物 安全性的主要因素。在大坝的设计和运行管理中，关于渗流的数值计算占有重要 的位置。地下水以空隙水的形式存在于坝基之中，对坝基产生静、动水压力，坝 基面渗透压力是大坝承受的主要载荷之一。目前设计规范规定，在计算闸坝地基 面渗透压力时，均采用简化的直线或折线分布。闸坝地基若不设置帷幕灌浆和排 水孔幕，其渗透压力呈折线分布，在排水孔幕处折减，折减系数与各种孔幕有关。 而用排水改变渗流场来达到改善建筑物的工作状态是工程设计中采用的主要措 施，排水孔和排水管是经常采用的排水结构。实际上，不管在哪种情况下，坝基 面渗透压力分布都与坝基面渗流场紧密相关，不同的坝基渗流场对应不同的坝基 面渗透压力分布。因此，应将二者紧密结合起来。要想知道静水压力或动水压力 对坝基的不利影响，必须准确获取坝基内的水压分布。但地下水渗流分析是一项 比较复杂的工作。地下水在某一方向的渗透性与另一方向相差很大，岩体内的地 下水运动存在渗流的各向异性。前人的方法分析渗流问题主要采用有限差分、有 限单元法、裂隙网络法。目前的研究工作多采用数值分析方法确定渗透压力，再 进一步计算其他的流动特性。常规边界元法分析坝基面及远坝基内点的渗透压力 已经获得了成功，然而精度不够理想，且存在边界奇异积分的问题。本文对平面 正交各向异性单一介质的渗流场进行了无奇异边界元分析。 计算思想 在稳态下，渗流问题可归结为达西方程边值问题的求解，达西定律的一般形式为： 1-1 式中，h 是渗流场中各点的水头值，V 是达西流速，K 是渗透系数。将式(1-1) 代入不可压缩流体在刚体介质中流动连续性方程，可得稳定渗流的微分方程式： 1-2 求解上述微分方程得到坝基渗流场水头分布后，可得坝基面上的渗透压力为： 1-3 式(1-3)中，p 为渗透压力，γ 为水的容重，Z 为位置水头(坐标)。若坝基面水平， 且取坝基面的位置水头: Z=0 ，则p=γh ，即渗透压力与水头成正比。 本文假设水土均不可压缩，各向渗透性为常数，由渗流计算分析知平面正交 各向异性稳定渗流的微分方程式：