基于MATLAB的滑坡三维有限元前处理研究-防灾减灾工程及防护工程专业论文.docxVIP

三 三 峡 大 学 硕 士 学 位 论 文 II万方数据 II 万方数据 内 容 摘 要 滑坡是一种常见、多发的地质灾害，其稳定性分析与评价是国内外学者关注的热 点问题，有限单元法是最常用的数值分析方法之一。因地形地貌复杂多变，如何快速 构建滑坡三维数值模型及有限元网格自动生成等前处理工作是该领域研究的难点问 题。目前大型商用有限元软件，如 ANSYS、ABAQUS 等，都有较强大的前处理模块， 但是在几何模型的建立与修改，网格的编辑等方面存在限制，特别是针对滑坡这类地 形地貌复杂的地质体构建有限元模型时需耗费大量的时间和精力，需开发一套专门的 滑坡三维有限元前处理程序。 本文重点研究了滑坡三维有限元前处理中的关键技术问题：三维空间数据的快速 提取及程序实现、基于不规则三角网和基于规则格网的数字高程模型（DEM）表面 建模技术及程序开发，基于滑坡剖面图的滑面 DEM 建模技术，滑坡三维有限元网格 自动生成技术，并基于 MATLAB 平台自主研发了一套滑坡三维有限元前处理软件 （PREFEM），其快速高效的特征使滑坡三维有限元建模、网格剖分等前处理工作效 率显著提高，在该领域具有重要的实践价值与应用前景。主要的研究工作及取得的成 果如下： 1）编写了具有版本兼容性的滑坡三维地形数据提取及可视化的 MATLAB 程序， 实现了滑坡三维地形特征的快速构造：在分析 AutoCAD DXF 格式文件的基本结构， 对比不同版本 DXF 文件的图元类型的基础上，利用 MATLAB 编写了具有版本兼容性 的 DXF 文件数据提取程序，建立 MATLAB 与 CAD 文件的数据接口。该 DXF 文件 数据提取程序可以从 CAD 地形图中直接提取建立 DEM 表面模型所需的三维空间数 据，具有快速高效，能够准确提取地形数据等优点，实现了复杂地形滑坡三维地形地 貌的快速高效的数字可视化。 2）编写了 DEM 表面模型建模的 MATLAB 程序，成功建立了基于不规则三角网 的 TIN 模型和基于规则格网的 Grid 模型。针对地形几何建模中常用的 Delaunay 三角 形开展研究，并基于 Delaunay 三角剖分规则构建了不规则三角网 TIN 模型；详细研 究了 DEM 内插模型中的逐点内插法和整体内插法，针对逐点内插法中插值邻域搜索 算法不收敛的情况，采用了变步长搜索策略，编写了逐点内插法和整体内插法的计算 程序，构建了基于规则格网的 Grid 模型。基于不规则三角网的 TIN 模型虽然能够更 逼真的表示复杂地形表面的起伏状态，但考虑到四面体网格在计算精度、划分数量、 抗畸变程度等方面明显不如六面体网格，所以建立滑坡几何模型时更多的采用基于规 则格网的 Grid 模型。 3）提出了基于滑坡剖面图建立滑面 DEM 模型的方法，该方法是在 DEM 表面模 III万方数据 III 万方数据 型的基础上，首先用三次样条曲线沿顺坡向进行纵向内插，内插出与表面模型格网数 据相匹配的剖面线，然后再用纵向内插的数据作为采样点，结合滑坡边界信息进行横 向内插，计算出滑面区域内的所有网格点的高程，从而建立滑面 DEM 模型。采用该 方法建立的滑面 DEM 表面光滑，边界处过渡自然，精度较高，能够更好的逼近真实 的滑面几何形态。 4）提出了基于多层 DEM 的滑坡三维有限元网格自动生成技术，并利用 MATLAB 编程实现了滑坡三维有限元网格的自动生成。该三维有限元网格自动生成技术是以基 于栅格法的六面体网格生成方法为基础，首先建立各个特征地层的 DEM 模型（DEM 表面模型、滑面 DEM 和滑床底面 DEM），形成多层 DEM，再对多层 DEM 进行缝合 和插值处理，最终形成三维有限元网格模型。在 MATLAB 中采用多层 DEM 的思想 建立滑坡几何实体模型比采用 ANSYS 自底而上或自顶而下的实体建模速度快、自动 化程度高，节约了大量的时间和精力，在建立滑坡这类地形地貌比较复杂的几何模型 时具有明显优势。再者，采用该技术生成的有限元网格相对比较规则，没有产生任何 的畸变单元，网格质量较高。 5）基于 MATLAB 平台采用可视化的界面环境（GUIDE）对滑坡三维有限元前 处理系统进行了 GUI 设计，自主研发了一套滑坡三维有限元前处理软件 PREFEM。 软件 PREFEM 包含数据准备、几何建模、网格生成、数据输出等四个基本模块，其 主要功能是根据滑坡平面图和剖面图快速建立三维有限元计算模型，输出模型离散文 件，同时，该软件能够直接从 AutoCAD 图形中直接提取 CAD 图形的几何实体信息， 还可以作为一种 CAD 数据接口程序使用。采用软件 PREFEM 对花莲树滑坡和生基包 滑坡进行有限元前处理，再把离散模型数据文件导入到 ABAQUS 进行有