论文部分内容阅读
土木建筑、水利水电、铁路公路交通、矿山开采、文物保护、地质灾害治理、能源开发及国防建设等众多的工程实践中,常常会遇到大量的岩土工程问题,尤其是工程岩体稳定性、岩体强度及变形的问题。作为主要工程对象的岩体是在漫长的地质历史发展过程中经受各种地质营力作用而形成的成因和构造复杂,岩性多样且包含了连续介质,裂隙介质和散体介质特征的,兼有弹性、塑性、粘性、流变性的各向异性的非均质地质体。岩体在不同的工程实践中往往表现为加载和卸载两种不同的力学状态,并且由于其固有的裂隙介质特征,又不可避免地赋存于一定形态的液态介质一各种地表水、潜水、承压水等中。由于实际工程岩体其固有的非线性和复杂性,使得求解对应的诸如其稳定性、变形等问题面临较大的困难,而有关的工程结构的设计必须保证该结构在外来因素的作用下具有足够的安全度、经济性和合理性,这必然要求对岩体及其工程结构的受力与变形有一套量化评判体系和优化技术,抓住实际工程问题的本质特征并建立可行的符合几何仿真、本构仿真、受力仿真、过程仿真四原则的求解模型,通过该模型的数值模拟成果来指导岩体工程的设计、施工及运营、管理;而岩体结构面的存在使得基于传统连续介质力学理论的理论分析和数值模拟面临巨大的挑战,物理模拟的试验周期和成本也大大增加,而岩土体工程问题则成为典型的数据有限、了解程度有限类问题,这类问题的解决需要综合应用理论分析、经验判断、物理模拟和数值模拟等方法,数值模拟可以完成目前许多技术手段无法完成的实验,如参数控制,复杂条件下的边界条件的处理,同时数值模拟具有高可重复性,且数值模拟的成本和人力开销等远低于物理模拟,因此研究岩土体工程问题的流形元数值模拟方法是一项具有理论和实际工程应用价值的重要课题。 流形元法是十年在有限覆盖系统基础上发展起来的新的岩体工程数值模拟方法,流形的基本结构是有限覆盖系统和覆盖间的接触函数。与通常的网格不同,覆盖可重叠或折叠产生一个组合空间,因此流形方法是追求一切空间上的所有解法。在流形一般化定义下,流形方法具有两种独立的网格——数学网格和物理网格。数学网格由折叠的覆盖组成,物理网格则由实际的材料边界、裂缝、块体、不同材料的交接面甚至不同的材料相的界面等组成。数学网格所含的数学覆盖定义了近似解的精度,由用户自由选择。常规的有限元分析网格等都可以转换成数学覆盖;而物理网格则定义了积分区域,不能人为选择。流形元法的最大优点可以统一处理连续变形与非连续变形问题, 同时可以统一解决有限元法、非连续变形分析及解析法的计算问题,具有广泛的应用前景。 本文以流形元法的基本理论为基础,主要研究对应的流形元有限覆盖系统的生成算法及流形元法的接触判断、接触转换、流形元法总体控制方程的建立与稀疏存储、动态连通性分析及分片直接求解,并基于泛型程序设计技术及面向对象的系统分析实现了对应的C++流形元法分析工具NMMG,并用具体的计算模型检验了流形元法分析的结果。 具体的研究内容主要为: 1.流形元有限覆盖系统生成法 在计算儿何的实际直线求交法、剪切算法及空间元素拓扑关系确定法等方法和非连续变形分析的裂隙切割块体生成算法的基础上,结合有向图遍历理论,本文构造并实现一个适用于以三角形常应变单元和四节点等参元为数学网格的流形元有限覆盖系统生成算法,并进行了相关算例的检验,结果表明该算法能满足实际问题的容错性及鲁棒性要求;并能较好地扩充到以其它有限元单元类型为数学网格的流形元有限覆盖系统的生成,开发了对应的C++流形元有限覆盖系统的生成类(elass ME_Builder)。2.流形元法中的接触分析 针对二维流形元法实际的角对角、角对边、边对边三种接触类型,在石根华提出的相关接触准则的基础上借鉴离散单元法中对接触的处理方法,利用接触判断的距离准则、角度准则及嵌入准则对流形元法中的接触进行处理,同时研究了接触状态在不同时间步之间的转换方法,开发了对应的e++接触处理类(elass ContaetS)。3.流形元法总体控制方程的建立 在流形元有限覆盖基础上推导了以三角形常应变单元和四节点等参元为数学网格的对应的各种应力边界条件、位移边界条件及各种荷载作用下的流形元单元刚度子矩阵,各种接触类型对应的刚度子矩阵、最终控制方程的集成和无嵌入、无张拉约束条件的处理方法以及对应的基于单纯形积分理论的数值积分法,开发了对应的C++类(class KMatrioes)。4.流形元法总体控制方程求解理论的初步研究 主要研究了基于有向图理论的流形元系统节点编码的优化、对应稀疏矩阵的压缩存储与流形元总体控制方程的动态连通性分析、分片直接求解法及流形元接触的开一合迭代方法,并开发了对应的e++类(elass Solver)。5.基于范型程序、模式设计技术的流形元分析系统的实现 结合最新的范型程序设计技术及模式设计技术,用C++实现了上述算法及必要的辅助类,并最终集成一个实用的二维流形元分析工具一确蒯G(而merieal Manifold Met