二维Biot固结方程的自然单元法求解.pdfVIP

第 42卷 第 11期 上 海 交 通 大 学 学 报 Vo1．42No．11 2008年 l1月 JOURNALOFSHANGHAIJIAOTONGUNIVERSIy NOv．2008 文章编号：1006—2467(2008)11-1880—04 二 维 Biot固结方程的自然单元法求解 褚衍标， 王建华 (上海交通大学 土木工程系，上海 200030) 摘 要：为了使 自然单元法能够应用于土体等多孔介质的流固耦合计算，通过结合Biot固结理论 及 自然单元法 自身特点，利用经典变分原理推导了固结微分方程的离散形式，并针对二维问题编制 了相应的计算程序．算例结果表明，自然单元法的结果与解析解吻合 良好，其精度高于有限单元法． 从而验证了自然单元法在固结分析中的正确性，拓展了自然单元法的适用范围． 关键词：自然单元法；Biot固结方程；自然相邻插值 ；经典变分原理 中图分类号：TU43；O241 文献标识码：A NaturalElementMethodforBiotPlaneCOnsOljdatiOnAnalysis CHUY 一biao， WANGJian—hua (DepartmentofCivilEngineering，ShanghaiJiaotongUniversity，Shanghai200030，China) Abstract：Thenaturalelementmethod(NEM)isanovelnumericalcomputationalmethodforsolvingpar— tialdifferentialequation．Itisbuiltuponthenotionofthenaturalneighborinterpolation，whichisbased onVoronoidiagram andDelaunaytriangulation．ThispaperfocusedonitsapplicationinsolvingBiotcon— solidationequation．Thediscreteform ofcontrolequationwasobtainedwithclassicalvariationprinciple； thealgorithm routinefor2D conditionwasalsoelaborated．Theresultsofnumericalexamplesshow that theresultsofNEM areinconcordancewiththeanalyticalsolutionandtheprecisionishigherthanthatof FEM ． Keywords：naturalelementmethod(NEM)；Biotconsolidationequation；naturalneighborinterpolation； classicalvariationprinciple Biot固结理论于 1941年建立，可反映孔隙压力 基于给定结点的Voronoi图和Delaunay三角剖分， 消散与土骨架变形之间的相互关系，被称为较完善 是以自然相邻插值为近似函数的无网格方法．作为 的真固结理论．但按 Biot方程求解固结问题的精确 无网格方法，NEM克服了有限单元法依赖于求解 解相当复杂，方程建立 6O多年来，仅在若干特殊情 区域网格划分的缺点，易于处理裂纹扩展、大变形和 况下求得了解析解 ，一般须采用数值方法求解． 移动相边界等方面的问题．同时相对于其他各种无 自然单元法(naturalelementmethod，NEM)是 网格方法，NEM计算量较小，并且 由于形函数满足 近 10多年发展起来的一种求解微分方程的新型数 插值性质，不存在一些采用移动最小二乘插值的无 值方法，由Braun等[1]于 1995年首先提出．该方法