基于GIS的边坡三维有限元分析前处理.pdfVIP

维普资讯 第 1期 矿 山 测 量 NO．1 2007年 3月 MINESURVEYING Mar．20o7 基于 GIS的边坡三维有限元分析前处理 刘庆元．韦建超 ．戴水财 (中南大学信息物理工程学院，长沙 410083) 摘要：三维有限元数值分析是最新发展起来的边坡稳定性定量分析方法，但其前处理工作非常繁琐， 从而影响了该方法的整体效率。文中对利用GIS技术，实现了从滑坡范围的确定到有限元网格的自 动生成的全部过程 ．提高了边坡玉维有限元计算的前处理工作效率做 了介绍。 关键词：边坡 ；有限元；GIS；前处理 中图分类号：P208 文献标识码 ：B 文章编号 ：1001—358X(2007)01—0029—02 边坡失稳是生态环境和工程建设 中经常遇到 的 数据模型是空间采样点按某种规则连接形成相邻互 广泛且严重的地质灾害．在对边坡稳定性进行分析 不重叠的三角形面。因为采样点可 以不规则分布， 的方法 中．有限单元法是较为常用的定量分析方法， 所以可 以处于最优位置．因此 TIN可以用比其他数 其用于边坡 中应力和位移 的分析 已经 比较完善 ，并 据模型更少的点准确的表示一个表面。由于边坡地 得到了许多学者成功的应用…。有限元用于边坡分 面受风化 ，雨水冲刷 ，人工开挖等影响，地面形状相 析．其基本过程主要包括边坡体的离散化、有限元方 对复杂 ．而地下水位面则 由于水 自重力 的作用下易 程求解 、计算结果 (位移 、应力 、塑性 区)的处理共 3 流动的原因，比较光滑 。地层面，断层只受岩层和地 个部分 。据统计 ，一般手工传统 的有限元法求解的3 质构造 的作用 ．相对光滑。所 以采用基于地表 TIN 个阶段所用的时间分别占总时间的40％一50％、5％ 的三角形网格 向下垂直延伸 ．各层面的三角划分 以 10％及50％ 55％。对于坡 的三维有 限元计算 ， 地面 TIN为基准 ，获取各层面，断面层的TIN网格 ， 主要是因为边坡表面 、滑动面等都十分不规则 ，使得 既考虑了地表起伏 ．又兼顾层面变化 ．满足计算精度 主要时间都用于边坡体的离散及结果的处理。实际 的要求 ．地面采样点不规则分布可 以反映边坡表面 上在边坡的有限元计算与分析 中，所涉及 的对象大 变化关键点的空间特征．以较高精度模拟地面，如图 多都为地理对象 ．所 以利用 GIS对地理对象的可视 1(a)所示 。地面 TIN向下垂直投影 ，TIN顶点的投 化管理和空间分析、以及可二次开发 的强大功能，进 影线与各层 的交点 z值可 以通过相应层面内插得 行基于 GIS的边坡有 限元分析前分析．将会大大提 到 。这样地面与各层面的TIN三角形具有垂直对应 高有限元分析的效率 ．同时也有利于将有 限元分析 关系 ，如 图 1(b)所示 。利用 VBA +ArcEngine对 结果与边坡空间信息的联系起来 ，实现 3DGIS平 台 ArcGIS的3D模块进行扩展可以实现该模型。 下的边坡有限元分析的可视化查询 、显示 。另外 ，基 于GIS的边坡有 限元分析系统 。因集成 了GIS的功 地面 能，更容易与GPS，RS，变形监测机器人等融合，对边 地层 坡变形进行更精确的监测。