用适体坐标变换方法求解渗流-中国科技论文在线.PDF

第26卷 第l期 岩 土 工 程 学 报 V01．26 No．1 2004年 1月 Chinese Journal of Geotechnical Er neering Jan．， 2004 用适体坐标变换方法求解渗流 Seepage analysis by boundary—fitted coordinate transformation method 介玉新，揭冠周，李广信 (if／华大学水利水电工程系，北京 looos4) 摘 要：运用传统的有限元方法求解渗流问题时，需要人为进行网格划分等前处理工作。当问题涉及自由变动边界时，求解比较困 难。本文将适体(Boundary-Fitted)s~标变换方法用于岩土工程渗流计算。通过求解Poisson方程自动生成贴合实际物理边界的计算 网格，并将渗流问题变换至规则统一的直角网格系统下进行有限差分离散和数值求解。对不同的几何边界可以实现统一的数值求 解算法，自动化程度高。计算实例表明，这种方法具有简单、快捷的优点，生成的网格质量较高，是一种很有意义的方法。 关键词：适体坐标；有限差分法；渗流；自由面 中图分类号：TV 139．14 文献标识码：A 文章编号：1000—4548(2004)01—0052—05 作者简介：介玉新(1970一)，男，博士，副教授，主要从事土工数值计算、土工合成材料等研究工作。 。 Jm Yu—xin，Jm Guan．zhou，LI Guang．xin (Department of Hydraulic~gineeriag，Tsinghua University，Beijing 100084，Q．na) Abstract：The traditional finite element method(F )involves constructing mesh bv hand，and the computation is c(m~plex when the analysis fe— lates tO seepage with free surface．In this paper，the finite difference method(FDM)based on the boundary-fitted cooldinate(BFC)transforma— tion is presented．The curvilinear grid system，with computational boundary coincident with the physical boundary，is numerically obtained by solving the Poisson equation．The seepage analysis is then done by FDM in a uniform tlansfonned Orthogonal O∞柑i瑚Ite system．and a uniform computer code is compiled which can be consistently applied tO problem in domain of any shape．Case studies showed that the me thod could get satisfactory results and was nlo]~efficient than traditional FEM． Key words_．Boundary-fitted~o~di[1atc；finite difference method(FDM)；seepage；free surface 饱和非饱和法等，但以上各种复杂的算法都具有程序 IJ^ -_■‘_ 刖 -_J 吾 编制困难和耗费机时的缺点…。 为分析堤坝等的稳定性以及选择防渗、排渗设计