›› 2013, Vol. 34 ›› Issue (7): 1867-1873.

• 基础理论与实验研究 • 上一篇    下一篇

岩土体渗流自由面问题的重心插值无网格方法

李树忱1,王兆清2,袁 超1   

  1. 1. 山东大学 岩土与结构工程研究中心,济南 250061;2. 山东建筑大学 工程力学研究所,济南 250101
  • 收稿日期:2012-06-29 出版日期:2013-07-10 发布日期:2013-07-15
  • 通讯作者: 王兆清,男,1965年生,博士,副教授,主要从事工程数值分析方法研究。E-mail: sdjzuwang@126.com E-mail:shuchenli@sdu.edu.cn
  • 作者简介:李树忱,男,1973年生,博士,教授,博士生导师,主要从事深部岩体力学特性和隧道灾害防治方面的研究与教学工作。
  • 基金资助:

    国家重点基础研究发展计划(973)项目(No. 2010CB732002);国家自然科学基金面上项目(No. 51179098);教育部博士点基金博导类 (No. 20120131110031);新世纪优秀人才支持计划资助。

Barycentric interpolation meshless method for solving rock and soil seepage free surface problems

LI Shu-chen1,WANG Zhao-qing2,YUAN Chao1   

  1. 1. Research Center of Geotechnical and Structural Engineering, Shandong University, Jinan 250061, China; 2. Institute of Engineering Mechanics, Shandong Jianzhu University, Jinan 250101, China
  • Received:2012-06-29 Online:2013-07-10 Published:2013-07-15

摘要: 岩土体的渗透破坏、地下工程的防渗设计等无不与渗流计算有关。针对渗流自由面问题,提出一种重心拉格朗日插值的配点型无网格方法。由于渗流自由面问题的求解区域是不规则区域,该方法通过将不规则求解区域嵌入一个正则矩形区域,在正则区域上采用重心拉格朗日插值近似未知函数,利用配点法离散渗流问题的控制方程,将重心拉格朗日插值的微分矩阵离散成代数方程表达的矩阵形式。将自由面上的边界条件通过重心拉格朗日插值离散,通过置换方程法和附加方程法施加边界条件,利用正则区域上的重心插值配点法,通过迭代确定最终自由面的位置。数值算例表明所提出的无网格方法对于求解渗流自由面问题的正确性和高精度。

关键词: 渗流自由面问题, 无网格方法, 重心拉格朗日插值, 微分矩阵, 配点法, 正则区域法

Abstract: The seepage failure of rock and soil and the seepage design for underground engineering are closely related to the seepage calculation. A meshless collocation method based on barycentric Lagrange interpolation for solving seepage free surface problems in 2D is presented. Embedded the irregular seepage domain into a regular rectangular domain, the seepage governing equation is extended into the regular rectangular domain. The unknown function of seepage in seepage domain is approximated by barycentric Lagrange interpolation in regular rectangular domain. Then, the barycentric Lagrange interpolation collocation method (BLICM) can be used to solve seepage problems in the regular rectangular domain. The governing equation of seepage problem is discretized by using BLICM into algebraic equations. Using notations of differentiation matrix and Kronecker product of matrices, the system of algebraic equations can be rewritten as matrix form. The boundary conditions are discretized by barycentric Lagrange interpolation and applied by using replacement equation method or/and additional equation method. The final location of free surface is obtained by iteration computation of BLICM in regular rectangular domain. The numerical results compared with others numerical methods indicate the correctness and high accuracy of the proposed method.

Key words: seepage free surface problem, meshless method, barycentric Lagrange interpolation, differentiation matrix, collocation method, regular domain method

中图分类号: 

  • O 302
[1] 张翠莲. 节理岩体力学性能不确定性分析方法研究[J]. , 2016, 37(9): 2721-2727.
[2] 冯 春 ,李世海 ,孙厚广 ,李志刚,. 颗粒接触无网格法及其在边坡成灾范围模拟中的应用[J]. , 2016, 37(12): 3608-3617.
[3] 申林方 ,王志良 ,谢建斌 . 基于边界配点法浅埋隧道开挖引起 周围土体变形的分析 [J]. , 2012, 33(S2): 297-301.
[4] 马文涛 ,李 宁 ,师俊平 . 基于单位分解的扩展径向点插值无网格法[J]. , 2012, 33(12): 3795-3800.
[5] 栾茂田 ,樊 成 ,黎 勇 ,杨 庆 , . 有限覆盖径向点插值法及其在土工问题中的应用[J]. , 2006, 27(12): 2143-2148.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
[1] 亓 乐,施建勇,曹 权. 刚性桩复合地基垫层合理厚度确定方法[J]. , 2009, 30(11): 3423 -3428 .
[2] 喻 军,童立元,刘松玉,汤劲松. 基于优势面理论的隧道控水与模拟分析[J]. , 2009, 30(12): 3825 -3830 .
[3] 杨 慧,曹 平,江学良. 水-岩化学作用等效裂纹扩展细观力学模型[J]. , 2010, 31(7): 2104 -2110 .
[4] 刘 琦,卢耀如,张凤娥,熊康宁. 动水压力作用下碳酸盐岩溶蚀作用模拟实验研究[J]. , 2010, 31(S1): 96 -101 .
[5] 张 鹏,陈剑平,邱道宏. 基于粗糙集的隧道围岩质量可拓学评价[J]. , 2009, 30(1): 246 -250 .
[6] 李 响,贾明涛,王李管,白云飞. 基于蒙特卡罗随机模拟的矿岩崩落块度预测研究[J]. , 2009, 30(4): 1186 -1190 .
[7] 尹盛斌,丁红岩. 深基坑施工时地表沉降预测的时序-投影寻踪回归模型[J]. , 2011, 32(2): 369 -374 .
[8] 霍 明,陈建兵,章金钊. 东北岛状多年冻土区公路路基清基试验研究[J]. , 2009, 30(S2): 263 -268 .
[9] 徐文杰. 大型土石混合体滑坡空间效应与稳定性研究[J]. , 2009, 30(S2): 328 -333 .
[10] 曹佳文 ,彭振斌 ,彭文祥 ,何忠明 ,尹 泉. 充气锚杆在砂土中的模型试验研究[J]. , 2011, 32(7): 1957 -1962 .