基于流形覆盖的岩体宏细观破裂的颗粒离散元法

本文基于流形覆盖的思想,结合颗粒接触模型,构造一种新的颗粒离散元数值方法即MPD(Manifold Particle Discrete),该方法可以从宏细观角度分析岩体破裂的全过程.通过分析流形覆盖和球颗粒模型,构造三维流形覆盖的四面体单元,得到MPD模型的平衡方程和相容方程.通过数值算例,验证了所提出的数值方法的正确性和分析岩石破裂的可行性.     

基于单位分解法的无网格数值流形方法

在数值流形方法和单位分解法的基础上，提出了无网格数值流形方法. 无网格数值流形 方法在分析时采用了双重覆盖系统，即数学覆盖和物理覆盖. 数学覆盖提供的节点形成求解 域的有限覆盖和单位分解函数；而物理覆盖描述问题的几何区域及其域内不连续性. 与原有 的数值流形方法相比，无网格数值流形方法的数学覆盖形状更加灵活，可以用一系列节点的 影响域来建立数学覆盖和单位分解函数，具有无网格方法的特性，从而摆脱了传统的数值流 形方法中网格所带来的困难. 与无网格方法相比，由于采用了有限覆盖技术，试函数的构造 不受域内不连续的影响，克服了原有的无网格方法在处理不连续问题时所遇到的困难. 详细推导了无网格数值流形方法的试函数和求解方程，最后给出了算例，验证了该方法的正 确性.     

基于FEMLIP的全风化边坡失稳破坏全过程的数值模拟研究

边坡坡角和强度是影响边坡稳定性的重要因素,而边坡失稳往往伴随着大变形的发生,其变形从数十米至数千米不等.目前,传统有限元法在处理大变形问题时常常因网格畸变而导致计算终止.因此,为了实现边坡失稳破坏全过程的模拟,并研究边坡坡角和强度对边坡稳定性的影响,基于Lagrange(拉格朗日)积分点有限元法（FEMLIP）,采用C语言编写了能够模拟边坡失稳滑塌全过程的Ellipsis程序,并通过一个典型案例对该方法的正确性和可行性进行了验证.采用该方法分析了边坡在不同坡角和强度条件下的稳定性和滑坡过程.研究结果表明,Lagrange积分点有限元法可以较准确地模拟边坡的潜在滑移面,并且可以模拟边坡失稳后的滑坡发展过程,为边坡滑坡大变形分析提供了一种新的数值计算方法.     

动态断裂力学的无网格流形方法

运用无网格流形方法求解动态断裂力学问题.该方法利用单位分解法和有限覆盖技术建立形函数，形函数的建立不受域内不连续的影响，可较好地求解裂纹问题.对于局部化问题，该方法的形函数构造较其他方法更为有效，避免了其他方法在建立试函数时没有考虑不连续尖端的缺点.由于采用有限覆盖技术建立试函数，该方法克服了不连续对试函数的影响，尤其当不连续变得复杂时，更能显示该方法在处理不连续方面的优点.在求解动态断裂力学问题时，弹性动力学积分弱形式的推导采用加权残数法，空间离散采用基于单位分解法的无网格流形方法，时间离散主要采用Newmark法.最后给出两个数值算例，将计算结果与解析解对比，说明该方法的正确性和可行性. 关键词： 有限覆盖 无网格流形方法 动态断裂力学 动态应力强度因子     

数值流形方法及其在岩石力学中的应用

数值流形方法是目前岩石力学分析的主要方法之一.该方法起源于不连续变形分析,主要用于统一求解连续和非连续问题,其核心技术是在分析时采用了双重网格:数学网格提供的节点形成求解域的有限覆盖和权函数;而物理网格为求解的积分域.数学网格被用来建立数学覆盖,数学覆盖与物理网格的交集定义为物理覆盖,由物理覆盖的交集形成流形单元.流形方法的优点在于它使用了独立的数学和物理网格,具有和有限元明显不同的定义形式,且数学网格对于同一问题不同的求解精度的需求可以很方便地细化.由于该方法考虑了块体运动学,可以模拟节理岩体裂隙的开裂和闭合过程,因而在岩石力学中得到了广泛应用,近年来许多学者对该方法进行了研究.本文简要叙述了节理岩体的数值方法从连续到非连续的发展过程,详细地介绍了数值流形方法的组成和数值流形方法在岩石力学及其相关领域的研究和发展概况,最后就作者所关心的一些问题,如三维问题的数值流形方法、数值流形方法在物理非线性问题和裂纹扩展问题中的应用、相关的耦合方法等进行了探讨.