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

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

一种建筑材料细观力学数值模拟的新方法

数值流形方法通过数学和物理双重网格,分析连续和非连续问题,已应用于模拟节理岩体裂隙的开裂与闭合问题.但对于裂纹尖端的局部化问题,数值流形方法需要像有限元那样在裂纹尖端设置细密单元.本文利用裂纹尖端解析解将数值流形方法的基函数进行扩展,推导了相应的试函数.从最小势能原理出发提出了断裂力学的数值流形方法,推导了相应的求解方程,将其应用于建筑材料细观力学数值模拟.最后给出两个数值算例,将计算结果与解析解对比,说明该方法的正确性和可行性.     

弹性力学的复变量数值流形方法

数值流形方法通过引入数学和物理双重网格，将插值域和积分域分别定义在两个不同的覆盖 上来完成系统能量泛函积分运算. 当采用高阶函数构造位移函数时，广义节点自由度将大大 增加. 在求解系统的平衡方程中，运算量是与自由度的三次方成正比的，因此数值流形方法 的计算量是较大的. 为此，在复变量理论的基础上，采用一维基函数建立二维问题的逼 近试函数，然后将其应用于弹性力学的数值流形方法，提出了复变量数值流形方法，推导了 弹性力学的复变量数值流形方法的公式. 与传统的数值流形方法相比，复变量数值流形方法 具有计算量小、精度高的优点.     

多裂纹扩展的数值流形法

数值流形法的求解体系建立在两套覆盖(包括数学覆盖和物理覆盖) 和接触环路的基础之上,实现了对连续和非连续问题的统一求解. 在处理裂纹问题时,数学覆盖无需与裂纹重合,方便岩体破坏过程的模拟. 通过在裂纹尖端影响区域内的物理片上增加用于模拟应力奇异性的增强位移函数,发展了扩展的数值流形法. 在此基础上,提出一种多裂纹扩展的控制算法,并给出了裂纹扩展过程中材料体的整体响应. 针对典型的线弹性断裂力学问题, 给出的数值算例表明所建议的方法是正确有效的.     

基于NMM的EFG方法及其裂纹扩展模拟

数值流形方法 (numerucal manifold method,NMM)通过引入数学覆盖和物理覆盖两套系统来统一处理连续和非连续问题.通过用移动最小二乘插值(moving least squares interpolation,MLS)中的节点影响域构造数学覆盖,得到了基于数值流形方法的无网格伽辽金法(element free Galerkin,EFG).该方法在保证前处理简单的同时,又能方便处理如裂纹等不连续问题.建立了适用于小变形和大变形的裂纹扩展计算格式,并通过对曲折裂纹(kinked crack)的处理,在不加密的情况下实现了任意小步长的裂纹扩展,大大提高了在固定网格中模拟裂纹扩展的实用性.大小变形的结果对比表明,按照不考虑构型变化的小变形计算,结果可能偏于危险.     

弹性力学中的一种非协调数值流形方法

通过引入数学和物理双重网格，将插值域与 积分域分别定义在不同的覆盖上，即在数学网格上进行插值函数的构造，物理网格上完成 系统能量泛函积分运算，最后通过覆盖权函数将二者联结在一起. 它的优点是单元网格划 分随意，不受复杂边界形状和二相材料界面的限制，单元可以是任意形状，是较之于有限 元方法更一般的数值模拟方法. 在4节点四边形数值流形方法中，由于单元总体位移函数 包含的完全多项式不完全，使得计算精度不够精确，为此，在单元总体位移函数上附 加非协调位移基本项，使之趋于完全，提出了弹性力学问题的一种改进的数值流形 方法------非协调数值流形方法. 通过内部自由度静力凝聚处理，导出了消除内参后的单元应变矩阵 和单元刚度矩阵，使得在不增加广义节点自由度的前提下，大大提高了数值流形方法的计 算精度和计算效率. 同时对非协调项进行了显式处理，可以对工程实践起到更切实的帮助. 数值试验表明，它们能够保证收敛，有较高的精度，对畸变不敏感，从而证明了该方法的 可行性.     

热传导问题的非协调数值流形方法

数值流形方法通过引入数学与物理双重网格，将插值域与积分域分别定义在两个不同的覆盖上，其优点是网格划分随意，不受复杂边界形状和材料界面的限制，是较之于有限元方法更一般化的数值模拟方法。在计算精度方面，数值流形方法远远高于有限元法。但它的精度还是不够理想。为此本文在单元总体位移场上附加非协调位移基本项，使单元位移函数趋于完全，构造了非协调流形单元来改善流形单元的计算精度和计算效率，并将其应用于热传导问题，推导了势问题的非协调数值流形方法。     

复杂裂纹问题的多边形数值流形方法求解

数值流形方法是一种能统一处理连续和不连续问题的有效数值方法。该方法采用的数学覆盖系统可完全独立于物理域，能很好地求解各类裂纹问题，而n边形单元（n>4）则具有网格划分灵活，求解精度高等优点。本文基于数值流形方法，采用正六边形数学单元求解线弹性复杂裂纹问题。在导出相关方程的基础上对典型裂纹问题进行了分析，通过互能积分法得到了裂尖的应力强度因子，计算结果与参考解吻合得较好。除此之外，文中还对不同单元上的求解精度进行了比较，结果表明采用正六边形单元的求解精度较正四边形单元和正三角形单元上的精度均更高。     

THE NMM-BASED EFG METHOD AND SIMULATION OF CRACK PROPAGATION

数值流形方法（numerucal manifold method,NMM）通过引入数学覆盖和物理覆盖两套系统来统一处理连续和非连续问题. 通过用移动最小二乘插值（moving least squares interpolation,MLS）中的节点影响域构造数学覆盖,得到了基于数值流形方法的无网格伽辽金法（element free Galerkin,EFG）. 该方法在保证前处理简单的同时,又能方便处理如裂纹等不连续问题. 建立了适用于小变形和大变形的裂纹扩展计算格式,并通过对曲折裂纹（kinked crack）的处理,在不加密的情况下实现了任意小步长的裂纹扩展,大大提高了在固定网格中模拟裂纹扩展的实用性. 大小变形的结果对比表明,按照不考虑构型变化的小变形计算,结果可能偏于危险.     

数值流形方法在流固耦合谐振分析中的应用

数值流形方法(流形法)是石根华博士利用现代数学中流形分析的有限覆盖技术建立起来的新的数值分析方法,统一解决了连续和非连续变形的力学问题,具有广阔的应用前景。本文将流形法应用于交界面耦合的流固振动分析,采用平面矩形数学网格,针对无粘、无旋、不可压缩流体和无阻尼的固体结构,提出分析流固耦合系统简谐振动的高阶流形法公式,其中,采用拉格朗日乘子法引入流场的已知边界条件。本文还初步研究了在特殊的无限远流场中采用解析解覆盖函数的实现技术。文中算例体现了流形法网格划分的方便性和计算的高精度,显示出流形法在数值解和解析解联合运用上的优势。     

EVALUATION ON STRESS INTENSITY FACTOR OF CRACK UNDER DYNAMIC LOAD USING NUMERICAL MANIFOLD METHOD

相较于传统有限元,数值流形方法（numerical manifold method, NMM） 的一个显著优点是在处理裂纹问题时网格无需与裂纹重合,这就方便了岩体破坏过程的模拟. 基于包含裂尖增强函数的NMM,采用Newmark 隐式动力学算法进行时间积分,重点研究了动力载荷条件下裂纹动态应力强度因子（dynamic stress intensity factor,DSIF） 的求解方法. 针对典型的线弹性动力裂纹问题,给出了NMM 的数值算例. 结果表明NMM 能够准确计算动载荷条件下裂纹的DSIF,并且具有较好的收敛性.     

基于适合分析T样条的高阶数值流形方法

数值流形方法是一种非常灵活的数值计算方法,连续体的有限单元方法和块体系统的非连续变形分析方法只是这一数值方法的特例.数值流形方法中高阶位移函数的构造可通过提高权函数的阶次来实现,这种方法往往需要沿单元边界配置适当的边内节点,这些结点的出现增加了前处理的复杂性,特别是对于大型复杂的空间问题.另一方面,在数值流形方法中可通过缩小单元尺寸(h加密)来提高求解精度.当模拟裂纹扩展时,这种细化策略可用来克服裂纹尖端的奇异性.一个传统的解决方案是细化整个网格,但这会导致计算效率的显著降低.将适合分析的T样条(analysis-suitable T-spline,AST)引入数值流形方法中来建立高阶数值流形方法的分析格式,有效的避免了该问题的出现.AST样条基函数具有线性无关,单位分解,局部加密等许多重要性质,使得其非常适合用于工程设计及分析.在引入AST样条后,可通过改变数学覆盖的构造形式建立不同阶次的数值流形方法分析格式;AST样条自身的局部加密性质也使得数值流形方法中的数学网格局部加密更容易实现.算例结果表明:随着AST样条基函数阶次的提高,数值流形方法的计算结果有了明显的改善;基于AST样条基函数的数值流形方法在保持计算精度的前提下降低了自由度的数量.     

冲击载荷作用下岩体破坏规律的数值流形方法模拟研究

为了更好地模拟岩体这种典型的连续与非连续介质共存的复杂结构体在冲击载荷作用下的破坏规律,在近年来新出现的数值流形方法基本理论的基础上,利用其算法程序对岩体中存在的不连续面如节理、裂隙等对岩体破裂效果的影响进行了数值模拟分析,计算结果充分显示了数值流形方法在模拟这类问题中的准确性和有效性。     

基于最佳质量网格的薄板问题的非协调流形方法

对于大部分非协调板单元，使用规则网格能得到很好的效果。但是，当网格不规则时，非协调元的数值特性将变得很差，甚至收敛性得不到保证。为解决网格依赖性问题，许多专家学者提出了改造单元，如拟协调元法和广义协调元法，这些方法能解决收敛性问题，但是数值实践证明没有一种单元能在所有情况下都具有良好的数值特性。考虑到流形方法采用两套完全独立的覆盖系统，可以用规则的数学网格来作为数学覆盖进行插值，取得最佳的插值效果，单元收敛性便能得到保证。再结合适用于流形方法的变分提法，建立起流形方法处理非规则物理边界非协调板单元的一般格式。以ACM薄板单元为例，与ANSYS、拟协调元法和广义协调元法进行了对比，证明本文方法在处理具有曲线边界的薄板弯曲问题时具有收敛快和精度高等优势。