基于ansys的楔形体稳定性计算程序的开发.docxVIP

基于ansys的楔形体稳定性计算程序的开发 1 几何模型及分析 刚体横截面平衡法历史悠久，被证明是分析边坡稳定性有用而可靠的方法之一。它具有明确的力学假设和丰富的工程实践经验，可以为工程中最受关注的稳定性系数提供可靠的测试调整措施。它对巩固边坡和规范的应用具有纠正措施。目前, 二维刚体极限平衡分析已在边坡稳定性分析中得到了广泛应用, 积累了丰富的经验并已开发出一系列应用软件。 然而, 在岩质边坡中, 岩体发生失稳破坏本质上是一个空间问题, 其主要形式为由空间几组结构面和临空面切割形成的岩块 (楔形体) 发生失稳, 这就要求根据岩体结构条件, 快速准确地判断出可能失稳块体的形状、规模、滑移边界、稳定性状况并及时提出整治处理措施。针对楔形体稳定性分析, 石根华等人提出了块体理论及DDA方法、, 徐明毅等人提出了危险块石理论及荷载计算方法、, 常用的分析方法包括赤平投影法、矢量分析法等。目前, 国内外已开发出了一些关于楔形体稳定性分析的计算机软件, 但在一些诸如复杂形状块体几何建模、最危险块体搜寻等关键问题上尚未得到有效解决。 笔者结合块体理论和刚体极限平衡理论, 基于通用有限元分析软件ANSYS的三维建模平台, 以工作平面代替“软弱结构面”对边坡进行“切割”, 简便快捷地建立任意复杂块体的三维空间模型;同时, 结合VB语言的可视化编程功能, 对原有楔形块体稳定性分析方法作了改进和完善, 开发出了一套计算分析功能较为完备的楔形体稳定性计算程序, 并以西部某水电站进水口边坡为研究背景, 考虑了块体的自重、锚固力、水的渗流作用及地震作用力等荷载作用, 开展了岩质边坡的稳定性分析, 提出了合理的加固处理措施。 2 3d刚性计算模型和计算参数的三维性质 2.1 斜坡坡度和底板高程 某水电站位于四川省甘孜州雅江县境内的雅砻江干流上, 电站坝址位于雅砻江干流与支流鲜水河的汇合口下游约2 km河段, 水库正常蓄水位高程2 865 m。右岸进水口部位斜坡总体走向为S10°W, 其中进水口底板以下斜坡坡度一般为40°～ 50°, 进水口底板以上斜坡坡度相对较缓, 一般为30°～35°, 底板高程为2 763 m。进水口边坡主要发育有断层f10、f34-1、f34-2、f34-7、f34-8, 边坡岩体中主要发育7组优势裂隙:①N55°～ 90°W/NE∠55°～ 85°、②N55°～ 90°W/SW∠55°～ 85°、③N0°～ 30°E/SE∠70°～ 85°、④N0°～ 30°E/NW∠70°～ 85°、⑤N0°～ 30°E/SE∠5°～ 35°、⑥N0°～ 30°E/SE∠40°～ 60°、⑦N0°～ 35°W/NE∠5°～ 30°, 其中断层f10、f34-3、f34-4及优势裂隙①、③、⑤是控制性软弱结构面。 2.2 工程边坡三维刚体极限平衡模型 为研究进水口工程边坡在空间的稳定性, 根据进水口边坡的地形地貌及开挖后的边坡形态, 建立了进水口工程边坡三维刚体极限平衡模型 (图1) , 以进水口工程边坡与底板交界面中心为界, X方向:向河谷延伸237 m, 向边坡山体延伸309 m;Y方向:向上游延伸447 m, 向下游延伸523 m;Z方向:与X和Y垂直, 铅直向上, 底部取至2 400 m, 上部取至地表, 最大高差770 m。 2.3 抗剪断强度参数确定 在基于三维刚体极限平衡法的岩质边坡稳定性分析中, 岩体物理力学参数主要根据岩体工程分类标准建议的力学参数和设计参数确定 (表1) ;优势裂隙力学参数和连通率由设计参数确定 (表2) 。边坡所发育的f34-1、f34-2等断层属于岩块岩屑型结构面;在抗剪断强度参数中, 摩擦系数f=0.47, 凝聚力c=0.13 MPa。而在正常运行期, 高程2 865 m以下边坡岩体及结构面处于饱水状态, 其容重采取饱和容重, 抗剪断强度的折减系数取0.8。 3 水流口前向的边坡组合的稳定性分析 3.1 进水口工程边坡荷载组合 鉴于该水电站进水口及开关站边坡为Ⅰ级永久性边坡, 根据《水利水电工程边坡设计规范》DL/T5353-2006, 该边坡稳定安全系数持久状况采用1.3～1.25, 短暂状况采用1.2～1.15, 偶然状况采用1.1～1.05。进水口工程边坡在不同工况下所受的荷载组合见表3, 其中暴雨工况采用块体高度折减法假设地下水位, 折减系数K分别取1/2、2/3;地震荷载采用拟静力法, 水平向峰值加速度设计状况采用0.137 g, 校核状况采用0.181 4 g, 地震作用效应折减系数取1/4。 3.2 随机质块体组合模式 结合进水口工程边坡区域地形地质条件, 考虑到边坡结构面众多, 主要为顺坡向且结构面倾角分布范围较大, 易形成空间滑动块体组合, 进水口洞脸边坡坡面产状为N7°W/NE∠5