改进的有自由面渗流问题的Bathe算法 ,郑宏戴会超刘德富.pdf

第26卷第4期 岩 土 力 学 Vo1.26 No.4 2005年4月 RockandSoilMechanics Apr.2005 文章编号:100)(一7598一(2005)04-0505-08 改进的有自由面渗流问题的Bathe算法 郑 宏 ,‘2戴会超1，刘德富1 (1.三峡大学，湖北宜昌443002;2.中国科学院武汉岩土力学研究所，武汉430071) 摘 要:建议了一个新的有自由面渗流问题的变分不等式提法，该提法通过将潜在出渗面上的边界条件提为Signorini型条 件，从而从理论上消除了出渗点的奇性，解决了出渗点的定位问题。同时在离散求解时，通过引进依赖于网格参数的连续型 Heaviside函数，克服了Bath。算法中所固有的网格依赖性，提高了这类方法的数值稳定性。 关 键 词:渗流;自由边值问题;Bathe算法;变分不等式 中图分类号:0357.3;TB115 文献标识码:A ImprovedBathesalgorithmforseepageproblemswithfreesurfaces ZHENGHong1,2,DAIHui-chaol,LIUDe-ful (1.ChinaThreeGorgesUniversity,Yichang443002,China;2.InstituteofRockSoilMechanics,ChineseAcademyofSciences,Wuhan430071,China) Abstract:Anewvariationalinequalityformulationisproposedforseepageproblemswithfreesurfaces.Theformulationenforcesa Signorinistypeofboundaryconditiononthepotentialseepageboundaries.Thesingulariytoftheseepagepointsishenceeliminated theoretically;andthelocationoftheseepagepointsisalsodetermined.Meanwhile,throughreplacingthestepHeavisidefunction withacontinuousHeavisidefunctiondependentonameshparameterduringdiscretization,themeshdependencyinherentinBathes algorithmisovercomeandthenumericalstabilityissignificantlyimproved. Keywords:seepage;rfeeboundaryproblems;Bathesalgorithm;variationalinequality 透系数法则类似于切线刚度法。这两类方法在本质 1前 言 上都是力图使干区内的渗流量远远低于湿区内的渗 有自由面的渗流问题的有限元法可分为两种解 流量。 法:调整网格法和固定网格法。 值得指出的是，一些直觉化的固定网格法[[6--9]r 调整网格法方法简单直观，至今仍然被广泛地 虽然在形式上网格保持不变，但却仍然采用的是调 使用。然而，由于在每一个迭代中都需要重新生成 整网格法的思想，其共同之处就是需要在迭代过程 网格，对非均质土体计算量通常都很大，且常常会 中确定自由面的位置，以期能够较准确地计算出被 导致计算结果的不收敛[[1]，特别是当需要与应力进 自由面切割的单元的渗透矩阵。 行祸合分析时，会带来诸多的不便，因而现正被固 有着严格数学基础的变分不等式法12〔-18]总是 定网格法所取代。当然，如果能够结合适体坐标变 试图构造一个定义在固定区域上的新问题，自由面 换技