基于MeshPy的3DGIS与三维有限元数值计算无缝耦合方法-岩土力学.PDFVIP

第39 卷第10 期 岩 土 力 学 Vol.39 No. 10 2018 年10 月 Rock and Soil Mechanics Oct. 2018 DOI ：10.16285/j.rsm.2017.0221 基于MeshPy 的3DGIS 与三维有限元 数值计算无缝耦合方法 1 1, 2 3 4 1 1 刘振平 ，杜根明 ，蔡 洁 ，周 凡 ，刘 建 ，卞 康 （1.中国科学院武汉岩土力学研究所 岩土力学与工程国家重点实验室，湖北 武汉 430071 ；2. 中国科学院大学，北京100049； 3. 湖北交投宜昌投资开发有限公司，湖北 宜昌443000 ；4. 湖北公路智能养护科技股份有限公司，湖北 武汉430014 ） 摘 要：3DGIS 与有限元数值计算耦合方法及其系统开发，是目前GIS 与岩土工程交叉学科的重要研究内容。通过对MeshPy、 GRASS GIS 、vtkPython 以及Scipy 等多个开源库的探索与编程实践，以Python 为“黏合剂”，构建了GIS 与三维有限元数值 计算无缝耦合的平台框架。在满足几何形状、Delaunay 特性、角度、面积以及体积等约束条件下，通过MeshPy 包含的Triangle 与TetGen 库获得能用于有限元数值计算的三角形与四面体网格，并阐述了限定Delaunay 三角化（CDT ）节点插入算法与四 面体网格质量控制方法。以排土场为例，针对CDT 网格建立的地表模型，结合vtkPython 封装的模块与算法，提出了根据节 点分布位置综合采用线性与三角形重心坐标的插值方法，有效消除了其他插值手段导致的部分地表高程畸变的现象。借助 Scipy 提供的矩阵运算模块，开发了邓肯-张E- 模型，给出了采用中点增量法进行非线性求解的算法，并通过对室内大型三 轴剪切试验的数值模拟验证了程序的正确性。最后，对准东露天煤矿北部排土场堆排过程开展三维有限元数值计算，结合 GRASS GIS 对排弃物料高度与厚度空间分布规律的分析结果，对排土场边坡位移变化规律进行初步研究，结果表明，排土 场垂直位移与排弃料高度及厚度相关，而最大水平位移则出现在各台阶坡体边缘位置，与实际情况一致。 关 键 词：三维地理信息系统 （3DGIS ）；开源软件；无缝耦合；有限元法 （FEM ）；排土场边坡 中图分类号：TD 73 文献识别码：A 文章编号：1000－7598 (2018) 10－3841 －12 Seamless coupling method of 3DGIS combined with 3DFEM simulation based on MeshPy 1 1, 2 3 4 1 1 LIU Zhen-ping , DU Gen-ming , CAI Jie , ZHOU Fan , LIU Jian , BIAN Kang （1. State Key Laboratory of Geomechanics and Geotechnical Engineering, Institute of Rock and Soil Mechanics, Chinese Academy of Sciences, Wuhan, Hubei