渗流自由面分析的比例边界有限元法.pdf

[文章编号】1001．246X(2009105．0665．06 渗流自由面分析的比例边界有限元法 李凤志 (南京航空航天大学航空宇航学院，江苏南京210016) [摘要】 为确定渗流自由面，利用比例边界有限元法对二维稳定渗流场进行分析．通过比例坐标与直角坐标 之间的转换，推导渗流问题的比例边界有限元方程．利用比例边界有限元空间降低一个维度，并在降维方向保 持解析的特点，只对求解域边界进行离散．当比例中心选在坝体上游直线边界和坝体下面不可渗透直线边界交 点时，只需离散自由面及其下游边界．并给出一种控制溢出点的控制点法．对二维坝体的稳态渗流场自由面问 题进行分析并与实验结果比较．结论表明，该方法收敛快、结果较精确、数据准备工作量小． [关键词】 渗流；自由面；比例边界有限元 [中图分类号】035；0241 [文献标识码】 A O 引言 土坝及地下水等渗流问题中，研究重点之一是确定渗流自由面．由于渗流自由面的未知性以及边界条件 的复杂性，该问题难于用解析法求解．目前工程上普遍采用有限元nJ、边界元幢]、无单元法¨o和差分法H。等数 值方法进行迭代求解．有限元法求解渗流自由面主要归结为移动网格法和固定网格法如残余流量法、结点虚 流量法和虚单元法等…．前者容易引起自由面附近单元形状畸变，甚至翻转重叠等问题，造成计算错误；后者 由于浸润线穿越单元，处理上存在一定的不便．传统的边界元法，只在边界上划分单元，克服了有限元法的不 但仍需对整个边界划分网格．无单元法与有限元法相比，实现了真正意义上的固定网格，但需花费更多的机 时，基函数的选取及结点布置对计算精度的影响等问题还有待于进一步研究．而差分法与有限元法类似，也 存在求解域全域离散和调整网格节点的问题． 比例边界有限元法是求解偏微分方程的一种新的数值方法．由W0lf和Songt副提出，并首次应用于弹性 力学计算中．该方法充分利用有限元和边界元法的优点，只需在边界上进行离散，降低了一个空间坐标的维 数，而在没有离散的坐标方向利用解析方法求解．相对于边界元法，它不需要基本解，也不存在积分的奇异性 问题，具有较高的计算精度．目前这种方法已被用于求解有限域、无穷域的弹性静力、动力学问题、断裂力学、 绕流场问题等口。“．本文首次将该方法应用于渗流问题的求解．以二维稳态自由面问题为例，验证了比例边 界有限元法求解渗流问题的有效性，计算结果与实验结果相符．当比例中心选在坝体上游直线边界和坝体下 面不可渗透直线边界交点时，只需离散拟定的自由面及其下游边界，从而将数据准备工作量降低到最小． 2 自由面渗流问题的数学模型 符合Darcy定律的二维稳定渗流问题的控制方程和边界条件为 V·(I|}V声)=0，在域Q内， (1) rHl， 在SI上， 声(菇，Y)=尹={也， 在s，上， (2) 【Y， 在S3和S．上， [收稿日期]2008—04—08；[修回日期]2008—10—20 [基金项目】国家自然科学基金及南航引进人才科研启动基金(S0717-014)资助项目 【作者简介]李风志(1975一)，男，黑龙江绥棱县，副教授，博士．从事工程中的热、质传递现象建模与数值仿真工作． 万方数据万方数据 计 算 物 理 第26卷 ， a声 一 ，、 2 一库n丽。qnu， 在S，和S：上， (3) 式中V为梯度算子，k为渗透系数张量，乒为总水头，n 为边界外法线方向．如图1所示，S。和s，分别为上游 V 和下游给定水头边界