数值模拟在岩体力学中应用.docVIP

数值模拟在岩体力学中应用 　　摘要:岩体由于各种结构面的存在导致其在工程结构和力学性能上呈现非均质、非线性、非连续的各向异性,因此在研究岩体在各种力场作用下变形与破坏规律时非常困难。数值模拟由于其可控性、无破坏性、安全性、允许多次重复等优点而在岩体力学中有着广泛的应用前景。文章从岩体的抗压试验模拟、抗剪试验模拟、岩体渗流几个方面介绍了近年来国内外这方面的研究成果。 　　关键词:岩体力学;数值模拟 　　中图分类号:TD313文献标识码:A文章编号:1006-8937(2010)14-0147-02 　　 　　岩体是在不连续面的切割下形成一定的岩体结构并赋存于一定的地质环境中的地质体。它与其他固体介质的重要区别在于岩体中存在着断层、节理、层面等各种结构面。这些结构面的存在大大改变了岩体的力学性能,导致了岩体在工程结构和力学性能上呈现非均质、非线性、非连续的各向异性等特性。要获得对岩体力学的基本力学参数,通常需要做原位试验。但是由于岩体结构的复杂性,以及试验研究成本昂贵,特别是试验结果可重复性差等因素,仅仅通过现场原位试验,是不能反映岩体工程中的地质信息和岩体力学的真实规律的。数值模拟由于其可控性、无破坏性、安全性、允许多次重复等优点而在岩体力学中有着广泛的应用前景。而且国内外的学者在这方面取得了不少的研究成果,它们主要表现在岩体力学的抗压试验模拟、抗剪试验模拟、岩体渗流试验模拟。 　　1岩体的抗压试验模拟 　　岩石是组成岩体的基本材料,而岩石的损伤和破坏是岩石力学中最复杂、最困难的课题,但是对岩石破坏机理的认识是控制和解决重大岩体工程灾害的关键。学者周辉等采用平面弹塑性细胞自动机模型在细观尺度上模拟应变软化岩石类材料的单轴压缩破坏过程,模拟了不同均质度岩石的单轴压缩破坏过程。得到了岩石单轴压缩条件下的全过程应力――应变曲线,开发了二维岩石弹塑性细胞自动机数值模拟软件。赵吉坤等基于应变空间理论导出了弹塑性损伤细观模型,提出了破坏单元网格消去法来模拟裂纹扩展,实现了岩石三维破裂过程的数值模拟。从而得出考虑岩石细观局部塑性变形时的岩体稳定性要大于不考虑局部塑性时的稳定性。学者赵吉坤指出通过考虑非均匀性特点研究岩石的非线性力学行为,是一种运用连续介质力学方法解决非连续介质力学问题的新型数值分析方法。细观模型并行算法在三维大理石三点弯曲梁无晶体和有晶体情况下破坏试验和三点偏弯曲梁弹塑性损伤破坏力学试验的实现,为此模型的实际工程应用奠定了良好的基础。 　　两位学者均指出岩石材料的非均匀性是影响岩石力学性能的主要因素。学者赵复笑等运用岩石破裂过程分析系统,研究了岩石介质非均匀性对岩石力学性质的影响。他指出:岩样的非均匀程度将会影响岩样的变形非线性行为和强度。{1}随着均质度系数的降低,宏观变形的非线性增强;随着均质度系数的提高,宏观变形的线性规律增强。{2}随着均质度系数的提高,岩样的宏观强度和弹性模量增加。{3}随着均质度系数的提高,主破裂前的前兆信息显著减小。{4}随着均质度系数的提高,岩石的宏观破裂表现为集中破裂模式,反之,表现出弥散性。 　　学者S.Skatulla把电活性聚合物(EPA)引入到岩石材料的非线性分析中,采取应变梯度方法的基础上广义连续框架并扩大到也包括电力耦合行为的EAP,建立了广义连续理论。广义连续理论的基本思路是考虑微观和宏观这两个空间跨度的广义空间。所有的数量包括本构方程是定义在这个广义的空间,而这项工作的常规电力耦合非线性超弹性分析是植根于广义连续。在这种方法下,材料信息的微型空间在这里只是几何规格的微型连续体,自然可以用本构方程表示其变化规律。该方法适应把不同的线性以及非线性超弹性构法、可以归结为电力耦合材料的行为以及规模效应编入处理。模拟数值实例表明了规模效应的重大影响。这些都是小的规模问题展示了像电活化聚合体的结构得到加强这样材料的机电耦合行为。例如:例如纤维等,在某些同向方向其硬度是增加的。但是他同时指出: 在现实的模拟,我们需要从EAP的专门实验中得到材料数据,他的理论并不能直接的应用。在今后的工作重点将是建立一个健全的本质功能函数,为机电耦合问题提供更进一步逼真模型。Yu-Geng Tang 提出了一种基于实地观测对岩土工程问题进行反分析的非线性优化技术(NOT),并结合额外的辅助手段以提高NOT的收敛性和稳定性。他利用有限元软件结合NOT和辅助手段成功地分析了岩土工程中的变形问题。 　　2岩体的抗剪试验模拟 　　岩石节理的抗剪强度是岩土工程设计中重要的力学参数,其影响因素复杂,研究节理抗剪强度与其影响因素之间的关系对于合理确定岩石节理抗剪强度参数具有重要意义。学者徐磊等在对岩石节理进行分形模拟的基础上,建立了岩石分形充填节理直剪试验数值分析模型,采用数值模拟方法,研究了充填度