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

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

Enriched meshless manifold method for two-dimensional crack modeling

The meshless manifold method is based on the partition of unity method and the finite cover approximation theory which provides a unified framework for solving problems dealing with both continuum with and without discontinuities. The meshless manifold method employs two cover systems. The mathematical cover system provides the nodes for forming finite covers of the solution domain and the partition of unity functions. And the physical cover system describes geometry of the domain and the discontinuous surfaces in the domain. The shape functions are derived by the partition of unity and the finite covers approximation theory. In meshless manifold method, the mathematical finite cover approximation theory is used to model cracks that lead to interior discontinuities in the displacement. Therefore, the discontinuity is treated mathematically instead of empirically by the existing methods. However, one cover of a node is divided into two irregular sub-covers when the meshless manifold method is used to model the discontinuity. As a result, the method sometimes causes numerical errors at the tip of a crack. To improve the precision of the meshless manifold method, the enriched methods are introduced in this work for crack problems.     

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

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

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

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

数值流形方法研究及应用进展

基于有限覆盖技术的数值流形方法是一种新的广义的数值方法.该方法的场函数近似原理和有限元、无网格、单位分解等方法相似,但在网格划分、覆盖形式、近似函数等方面有其自身的特点和优势.对该方法近年来在理论研究和应用方面取得的重要进展进行了综述.在理论研究方面, 目前已对不同形式物理覆盖流形单元的性能进行了研究,结果表明流形单元的精度较有限单元高,且提高覆盖函数的阶次能提高单元的精度;同时理论研究已由二维低阶流形方法推广到三维高阶流形方法,由线性流形方法推广到非线性流形方法,由基于能量原理的流形方法推广到基于加权余量的流形方法,非协调流形方法、无网格流形方法等也已开展了研究; 此外,覆盖系统的自动生成、覆盖函数的形式以及边界条件的处理方法等流形方法相关理论的研究也取得了进展.在应用方面,开展了有关岩石破坏和裂纹扩展等非连续变形分析更深入的研究,并已逐步推广到金属塑性变形分析、多孔介质变形分析以及温度场的数值分析等多个领域.针对目前流形方法的研究和应用现状,该文展望了流形方法理论及实现方法的研究方向、及其在计算流体力学、金属成形等大变形问题、多物理场分析等领域的应用前景.      

二维定常不可压缩粘性流动N-S方程的数值流形方法

将流形方法应用于定常不可压缩粘性流动N-S方程的直接数值求解,建立基于Galerkin加权余量法的N-S方程数值流形格式,有限覆盖系统采用混合覆盖形式,即速度分量取1阶和压力取0阶多项式覆盖函数,非线性流形方程组采用直接线性化交替迭代方法和Nowton-Raphson迭代方法进行求解.将混合覆盖的四节点矩形流形单元用于阶梯流和方腔驱动流动的数值算例,以较少单元获得的数值解与经典数值解十分吻合.数值实验证明,流形方法是求解定常不可压缩粘性流动N-S方程有效的高精度数值方法.     

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

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

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

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

基于流形覆盖思想的无网格方法的研究

本语言基于流形思想,利用有限覆盖,单位分解等概念,引入建立在覆盖上的覆盖函数和具有紧支撑特性的单位分解函数,建立场逼近的近似表达,由弱形式的Galerkin变分得到数值分析模型,结合边界条件用于边值问题的求解,由此建立了一类新的无网格数值方法,论文采用这种方法分析了平面弹性问题,分析了体积闭锁现象,h、p型收敛性等,提出了一种选择覆盖大小的方案,且对狭长城采用了椭圆覆盖形式,取得了比较好的效果。     

数值流形单元法数学网格自适应

基于数值流形方法和有限覆盖技术,将有限元法的后验误差估计理论及h型网格自适应技术推广应用到数值流形单元法中,提出了数值流形单元法的后验误差估计方法和数学网格自适应技术,并编制了相应的程序。数值算例表明,经过网格自适应,可以在粗糙的初始网格基础上得到质量比较理想的网格,计算结果可达到用户要求的精度。     

THEORETICAL STUDY OF THREE-DIMENSIONAL NUMERICAL MANIFOLD METHOD

IntroductionNumerical manifold method is a new numerical method established on the basis of finitecover of manifold[1,2].By using continuous and non-continuous finite cover system,thenumerical manifold method includes the continuous and non-continuous as …     

瞬态热传导问题的精细积分数值流形方法研究

数值流形方法（NMM）因其特有的双覆盖系统（数学覆盖和物理覆盖）在域离散方面具有独特的优势,而精细时间积分法则具有精度高、无条件稳定、无振荡以及计算结果不依赖于时间步长等特点。发展了用于研究二维瞬态热传导问题的精细积分NMM。结合待求问题的控制方程和边界条件,并基于修正变分原理导出了NMM的总体方程,给出了求解此类时间相依方程的精细时间积分及空间积分策略,选取了两个典型算例对方法的有效性进行了验证,结果表明本文方法可以高效高精度地求解瞬态热传导问题。     

流形元法固定边界约束处理研究

在数值流形方法中,对于材料的固定边界,一般采用罚函数的方法进行处理,即在固定边界上设置刚性弹簧约束其位移来实现固定约束条件的近似满足。罚函数法在理论上不是严格的固定约束处理方法,罚弹簧的布置与弹簧刚度的大小对模拟的效果都会产生影响。基于流形单元上位移函数的组成提出了流形方法固定边界约束处理的新方法,在组成流形单元的物理覆盖上,通过取消相应的覆盖函数在流形单元位移函数中的组成来实现双向固定的约束条件,通过使用只包含单方向位移的覆盖函数使x向固定约束条件和y向固定约束条件得到实现,推导了相应固定约束条件下的流形单元刚度矩阵的数值计算格式。该方法严格满足固定约束的物理意义,简化了固定边界的处理,并经算例证明是有效和准确的,有利于数值流形方法的程序实现和工程应用。     

有限覆盖径向点插值方法理论及其应用

数值流形方法能够统一地处理连续与非连续变形问题,有限覆盖技术是这种方法的核心。无网格方法前处理过程比较简单,径向点插值法是其中的一种计算格式。本文将有限覆盖技术与径向点插值方法相结合发展了有限覆盖径向点插值无网格方法,综合了数值流形方法与点插值方法的各自优点,能够有效地处理连续与非连续性问题,由此所构造的形函数具有Kronecker δ-函数属性,能够有效地处理位移边界条件。本文在阐述了这种方法基本原理的基础上,通过算例分析与数值计算论证了本文所建议方法的可靠性及其有效性。     

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

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

多裂纹扩展的数值流形法

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

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

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

2D numerical manifold method based on quartic uniform B-spline interpolation and its application in thin plate bending

A new numerical manifold （NMM） method is derived on the basis of quartic uniform B-spline interpolation. The analysis shows that the new interpolation function possesses higher-order continuity and polynomial consistency compared with the conven- tional NMM. The stiffness matrix of the new element is well-conditioned. The proposed method is applied for the numerical example of thin plate bending. Based on the prin- ciple of minimum potential energy, the manifold matrices and equilibrium equation are deduced. Numerical results reveal that the NMM has high interpolation accuracy and rapid convergence for the global cover function and its higher-order partial derivatives.     

一个新的八节点非协调曲边四边形平面单元

本文在文[1]和文[2]的基础上,提出构造非协调有限元的新方法。该方法不用一般的变分原理,可适用任意变系数正定和非正定偏微分方程。利用这一方法得到一个新的八节点四边形平面应力单元。与一般有限元相比,位移和应力可提高一阶收敛精度。形成单刚矩阵时,不需要进行数值积分。单元之间的协调条件容易满足,文中给出收敛性证明。文末给出数值算例,表明利用本文的方法,应力和位移均可获得满意的数值精度。     

四节点四边形数值流形方法及其改进

本文对四节点四边形流形元提出了改进措施,将覆盖位移函数用自然坐标表示,使得在一般非规则有限数学覆盖网格下,数值积分变得比较容易,克服了现有四节点四边形流形单元数值积分困难的缺点。数值算例将其应用于复合材料数值模拟,计算结果表明,当覆盖位移函数采用完全一阶等参多项式时,计算精度较传统有限元法有很大改进。在力应集中或应力突变的区域,无需网格加密,只需提高覆盖位移函数的阶次即可。