当前岩石力学数值计算方法应用探讨.pdfVIP

维普资讯 2007年第 12期 西部探矿工程 当前岩石力学数值计算方法应用探讨 唐广慧 ，刘发祥 ，唐升贵 ，宋 旭 ，伍锡举 (1．昆明理工大学N土资源工程学院，云南 昆明650093；2．贵州省有包地质勘查局，贵州 贵阳550005) 摘 要：对几种常用岩石力学数值计算方法的优越性和局限性进行探讨，分析了岩石力学数值计算方 法存在的瓶颈和当前解决这些问题的途径。 关键词：岩石力学；数值计算；应用探讨 中图分类号：TU458+．4 文献标识码：B 文章编号：1004—5716(2007)12—0025—05 当前岩石力学数值计算方法得到迅猛发展，出现了 将地壳岩石看成粘性流体，对青藏高原中段进行了构 有限差分、有限元、边界元、离散元、块体元、无限元、流 造应力场和速度场的有限元模拟，模拟结果与构造勘 形元及其混合应用等各种数值模拟技术，使复杂岩石力 察、地震活动和GPS测量结果具有很好的一致性，较 学工程问题的设计发生了根本性的变化u】。岩石力学 好地揭示了青藏高原构造活动规律、断层活动强度和东 数值计算不仅已基本上取代了传统线弹性力学实验，而 向位移的驱动机理。 且也在岩石工程非线性实验中显示出极大的优势。在 1．2 边界元法L2] 中国，有限元法已不仅 由线性发展到高度非线性和大变 边界元法在 20世纪70年代中期得到迅速发展，在 形 问题，由二维发展到三维，同时，还可 以考虑粘性流 处理半无限域或无限域问题方面非常方便，适用于解决 变、渗流与应力场耦合、损伤、断裂以及动力效应。国内 岩石力学问题尤其是岩石力学中地下开采的有关问题。 出现 了一些享有声誉的有限元程序，如北京大学 该法只在求解区域的边界上进行离散(剖分单元)，这样 NOIM一83，清华大学的TFINE，西安科技学院NCAP 就把考虑问题的维数降低了一维，这也是边界元法的优 一 2D程序等，这些程序均在许多大型岩石工程中得到 点。另外，边界元法计算精度高，应力和位移具有同样 了良好应用。不同数值计算方法的结合，更能发挥各种 的精度。当仅需要知道物体 内部个别点的解时，有限元 数值方法优势互补的作用，如有限元一边界元的混合、 仍不得不剖分整个物体才能确定个别点的解，而边界元 有限元一离散元的混合、有限元一无限元和有限元一块 则可以在已知边界上的解后根据需要去求物体 内部预 体元的混合采用等。 知点的解，这就可以比有限元大大地节省计算量和费 然而，由于岩体具有非连续、非均质、各向异性、天 用。但边界元法的系数阵是满阵，远比有限元刚度矩阵 然初始地应力影响、地下水影响及复杂边界条件处理等 (带状稀疏阵)的结构复杂。对于面体比比较大的薄壁 诸多复杂性使得当前岩石力学数值计算仍然是一个值 结构等物体，边界元法不如有限元优越。且边界元法对 得探讨的问题。 奇异边界难于处理。另外，边界元法对变系数、非线性 l 常用岩石力学数值计算方法应用分析 等问题较难适应，且它的应用是基于所求解的方程有无 1．1 有限元法L2 基本解 ，因此 ，限制 了边界元法在更广泛领域 的应用。 1966年，布理克(W．Blake)最先应用有限元法解 然而，边界元法和其它数值方法联合使用，可为解决愈 决地下工程岩石力学问题。目前，在岩石力学数值计算 来愈广泛的问题开辟了新的途径。例如，在线弹性区域 方面，有限元法主要用来求解线弹性、弹塑性、粘弹塑 或无限域、半无限域采用边界元法，在非线性的区域采 性、粘塑性等问题，是地下工程岩体应