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

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

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

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

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

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

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

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

Biot平面固结分析的流形单元法

基于数值流形方法和有限覆盖技术，提出了适用于Biot固结分析的三节点平面协调流形元。由于土骨架位移和孔隙水压力的节点覆盖函数(Lagrange插值函数)阶次可分别任意选择，该单元是一组满足位移和孔压插值阶次不同且所有节点具有相同自由度数的新型u-p混合模式单元，并且更加方便编程。数值分析表明，位移和孔压的节点覆盖函数阶次分别取一次和零次的流形单元(T1-0)是该组单元中最为有效的。与等价四边形等参元相比，T1-0流形元能给出精度更高的初期孔压和位移。     

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

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

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

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

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

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

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

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

混凝土断裂的连续-非连续方法

采用有限元形函数作为单位分解函数,位移间断用富集节点的附加自由度表示,建立了允许在单元内部位移非连续的局部富集公式以表征混凝土的开裂区域．富集基函数由节点形函数和节点形函数与间断函数的乘积的并集构成．非连续位移的扩展路径完全与网格结构无关．不同于以非协调应变为基础的嵌入非连续模型,对单元的类型没有限制而且间断位移可以贯穿单元边界．局部富集思想与扩展有限元类似,但富集点自由度保持节点位移的物理意义不变,使相邻单元无需进行富集运算．在变分公式中引入混凝土粘结本构定律,推导了考虑断裂过程区非线性影响的基本方程．对混凝土粘结裂纹扩展的数值模拟说明了该计算方法的有效性．     

Incompatible numerical manifold method for fracture problems

The incompatible numerical manifold method （INMM） is based on the finite cover approximation theory, which provides a unified framework for problems dealing with continuum and discontinuities. The incompatible numerical manifold method employs two cover systems as follows. The mathematical cover system provides the nodes for forming finite covers of the solution domain and the weighted functions, and the physical cover system describes geometry of the domain and the discontinuous surfaces therein. In INMM, the mathematical finite cover approximation theory is used to model cracks that lead to interior discontinuities in the process of 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 INMM 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 INMM, the analytical solution is used at the tip of a crack, and thus the cover displacement functions are extended with higher precision and computational efficiency. Some numerical examples are given.     

An effective approach for modeling fluid flow in heterogeneous media using numerical manifold method

A major challenge of modeling fluid flow in heterogeneous media is to model the material interfaces, which may be arbitrarily oriented or intersected with Dirichlet, Neumann, or other boundaries, making it difficult to mesh and accurately satisfy the boundary constraints. In order to solve these problems, we derived a new continuous approach in the numerical manifold method (NMM). NMM is an ideal method to handle boundaries, considering its flexibility and efficiency with fixed mathematical mesh and its integration precision. With the two‐cover‐meshing system, we construct physical covers containing gradient jump terms defined as extended degrees of freedom to realize the refraction law across material interfaces. In the global equilibrium equations, the jump terms are naturally considered with the energy‐work seepage model. In this approach, high accuracy is expected from the newly constructed jump function together with simplex integration. Moreover, high mesh efficiency is realized by fixed triangular mathematical mesh with algorithms fully considering interfaces intersecting with Dirichlet, Neumann, or other boundaries and simplex integration on elements in arbitrary shapes. The new approach was coded into our NMM fluid flow model. We calculated examples involving fluid flow through a domain including (1) a single interface, (2) an idealized fault represented by multiple material interfaces, (3) intersected interfaces, and (4) an octagonal inclusion. We compared the simulated results to analytical solutions or results with denser mesh to test precision and efficiency and thereby proved that the new approach is accurate, efficient, and flexible, especially when considering intense geometric change or intersections. Copyright © 2015 John Wiley & Sons, Ltd.     

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

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