裂纹扩展分析的无网格流形方法

运用无网格流形方法求解裂纹扩展问题。该方法利用单位分解法和有限覆盖技术建立形函数。形函数的建立不受域内不连续面的影响,可较好地求解裂纹问题。尤其当这种不连续面变得复杂时,更能显示该方法的优点。对于应变局部化问题,该方法的形函数构造较其他方法更为有效,避免了其他方法在建立试函数时不能考虑不连续尖端的缺点。与传统的数值流形方法相比,无网格流形方法的有限覆盖形状更加灵活。它可以用一系列节点的影响域来建立有限覆盖和单位分解函数,具有无网格方法的特性,从而摆脱了传统数值流形方法中网格所带来的困难。与目前的无网格方法相比,由于采用了有限覆盖技术,试函数的构造不受域内不连续面的影响,克服了原有的无网格方法在处理不连续问题时所遇到的困难。在求解裂纹扩展问题时,弹性力学基本方程的弱形式采用加权残数法获得。最后利用无网格流形方法追踪岩体试件在复杂应力状态下裂纹扩展过程。在此给出了数值算例,并将计算结果与实验结果进行对比,说明该方法的正确性和可行性。     

弹性动力学的无网格流形方法

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

考虑裂纹尖端场的数值流形方法

数值流形方法起源于不连续变形分析,主要用于统一求解连续和非连续问题,其核心技术是在分析时采用了双重网格:数学网格提供的节点形成求解域的有限覆盖和权函数;而物理网格为求解的积分域。由于该方法考虑了块体运动学,因此就可以模拟节理岩体裂隙的开裂与闭合过程。但对于裂纹尖端的局部化现象,数值流形方法需要像有限元那样在裂纹尖端设置细密单元。本文在单位分解法的理论基础上,应用裂纹尖端局部函数来扩展原有的数值流形方法的基函数,提出考虑裂纹尖端场的数值流形方法。本文方法扩展了原有数值流形方法对裂纹尖端问题的求解能力,同时对非连续问题也比原有数值流形方法的求解精度高。     

单位分解法、无网格法、数值流形方法之形函数的内在联系

对目前广泛采用的单位分解法、无网格法、数值流形方法的形函数的形成作了分析,据此分析了它们之间的内在联系,为综合应用各种方法提供了必要的思想基础.     

有限覆盖无单元法在裂纹扩展数值分析问题中的应用

有限覆盖无单元法是一种基于有限覆盖技术和无单元法的数值计算方法。有限覆盖技术是数值流形方法的基础,由于数值流形方法具有统一处理连续与非连续问题的能力,而无单元法的前处理比较简单。因此有限覆盖无单元法综合了数值流形方法与无单元方法的优点,能够更有效地处理非连续性问题。本文简要阐述了有限覆盖无单元法的基本理论,着重将这种方法应用于应力强度因子计算和裂纹扩展模拟问题。若干算例数值计算结果表明了这种方法的有效性。     

数值流形方法的网格自动剖分技术及其数值方法

数值流形方法包含数学覆盖与物理覆盖双重网格，数学网格用以构造流形单元的插值函数，物理覆盖确定了流形单元的积分区域。数值流形方法的前处理一直是一个难题。讨论了数值流形方法的网格自动生成技术，解决了数值流形方法的前处理问题。无论是连续性材料还是非连续性材料，数学覆盖都保持不变，因此，借助现有的有限元技术生成数值流形方法的数学网格，利用面向对象的编程方法生成了数值流形方法的物理网格。实例应用表明，这种方法是可行的和有效的。     

四面体有限单元覆盖的三维数值流形方法

构造三维流形单元的覆盖函数和总体位移函数，基于最小势能原理建立四面体有限单元覆盖的三维数值流形方法分析格式，详细推导了三维流形单元的刚度矩阵、初应力矩阵、荷载矩阵、惯性矩阵、位移约束矩阵和接触矩阵，并给出了相应的公式和算例。     

流形切割及有限元网格覆盖下的三维流形单元生成

 三维流形单元生成和接触搜索算法问题是制约三维数值流形方法发展的瓶颈问题。系统详细地研究三维流形单元的生成方法,在前人工作基础上,采用三维有限元网格生成技术生成数学网格;通过块体数据结构、块体识别算法等方面的改进,将三维块体切割技术发展成流形切割技术,来解决流形块体生成问题;将石根华博士在二维NMM程序中采用的物理覆盖系统编码算法,扩展成三维流形编码算法,进而实现三维流形单元的生成。并在此基础上开发三维流形切割程序3D_MC.f90,可以实现四面体及六面体网格覆盖下任意形状三维流形单元的生成。通过几个例子可以看出,三维流形切割程序生成的流形块体形态、流形单元的节点与单元编码等均满足三维流形单元定义要求,从流形切割角度来看,说明此算法是正确的。     

流形元方法尖模拟裂纹扩展中的应用

介绍了数值流形方法的基本原理及其在模拟裂纹扩展中的实现过程，并就两上算例进行了分析，结果表明流形元方法的有效性。     

流形元方法在模拟裂纹扩展中的应用

介绍了数值流形方法的基本原理及其在模拟裂纹扩展中的实现过程，并就两个算例进行了分析，结果表明流形元方法的有效性。     

高阶扩展数值流形法在裂纹扩展中的应用

 为正确预测工程材料体中裂纹扩展趋势,提出改进的数值流形法。该方法在每个物理片上采用位移函数的一阶泰勒展开形式,使得物理片上的自由度均具有明确的物理意义;且在裂纹尖端附近的物理片上增加扩充位移函数,用于模拟裂纹尖端的应力奇异性,可更为准确地预测裂纹扩展方向。同时,也提出裂纹扩展过程中适用于大小变形的物理覆盖系统的更新算法。针对典型的线弹性断裂问题,给出用该方法求解的数值算例。结果表明,预测的裂纹扩展路径与已有研究结果一致,从而证实方法的有效性及正确性。     

无网格法在断裂力学中的应用

采用罚函数法满足无网格法的本征边界条件,给出了罚参数的选择方案。对单裂纹体的应力场进行了分析,讨论了积分围线对应力强度因子计算结果的影响。数值结果表明,该方法可以有效地计算裂纹尖端的应力集中和应力强度因子,并能准确地模拟裂纹的传播。     

有限元法、无单元法及自然单元法之比较研究

无单元法与有限元法的根本区别在于形函数的构造方法不同，因而得到的形函数形式也不同，但可统一于单位分解。另外，采用Galerkin法对偏微分方程求解时需要对弱形式在求解域内积分，研究表明有限元法和无单元法中对弱形式的积分可统一于单位分解积分方法。无单元法解决了有限元法中前处理复杂以及难以有效处理的诸如裂纹扩展和大变形等问题；然而自身也存在着计算量大、边界条件和不连续的场变量及其导数处理上的困难。自然单元法是新近出现的一种求解PDE的数值方法，它兼有无单元的特性和有限元的优点，可以认为是介于两者之间的一种极具发展前途的数值方法。     

岩石块体数值流形分析网格形成方法之研究

根据数值流形方法的基本理论 ,系统地研究了数值流形元覆盖系统形成的方法。基于平面三角形有限元网格 ,探讨了岩石块体系统数值流形方法分析之数学覆盖和物理覆盖及其流形元分析网格的形成技术 ,研究了覆盖系统两套覆盖编号的自动形成方法。据此研制成计算机分析软件并应用于计算实例 ,取得良好效果。该文所取得的研究成果为数值流形方法应用于大型岩土工程分析与计算奠定了坚实的基础 ,具有广泛的工程应用前景     

四个物理覆盖构成一个单元的流形方法及应用

将由三个物理覆盖构成的流形单元的数值流形方法推广至四个物理覆盖的情形,给出具体的数值实现过程,并将其应用于裂纹扩展的模拟,数值结果令人满意。     

岩体工程数值流形方法的固定边界约束处理方法

在岩体工程数值流形方法分析中,往往存在较多的固定约束边界。目前,数值流形方法一般采用罚函数方法处理固定边界问题,但罚弹簧的布置与大小对数值模拟的效果有一定的影响,且处理方式也比较复杂。基于流形单元上位移函数的组成提出固定约束处理的新方法,将固定边界的约束处理转化为对广义节点的约束处理,改变广义节点上的覆盖函数使固定边界的约束条件得到严格满足,并推导相应的流形单元刚度矩阵。该方法在物理意义上严格满足固定边界的约束条件,同时简化了处理工作,有利于数值流形方法程序的实现和工程应用。