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

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

HIGHER-ORDER NUMERICAL MANIFOLD METHOD BASED ON ANALYSIS-SUITABLE T-SPLINE 1)

Numerical manifold method (NMM) is a very flexible numerical method which contains and combines finite element method (FEM) and discontinuous deformation analysis (DDA).High-order numerical manifold method can be constructed by increasing the order of the weight function.This method often needs to configure the appropriate edge nodes along the element boundary, the emergence of these nodes increase the complexity of pre-processing, especially for large and complex spatial problems.On the other hand, the level of approximation of NMM can be improved by splitting the elements into smaller ones (known as h-refinement).With regard to the h-refinement, a cover refinement strategy is necessary to overcome the singularity of the stress when simulating crack propagation in NMM.One traditional solution is to refine the entire mesh which can lead to a significant decrease in the computational efficiency.In this paper analysis-suitable T-spline is introduced into NMM and regular rectangular meshes are used as the mathematical cover system.Specifically, analysis-suitable T-spline is linearly independent, forms a partition of unity, and can be locally refined which make it meet the demands of both design and analysis.The basis function of analysis-suitable T-spline is adopted as the weight function in NMM to construct high-order NMM and make the local refinement for feasible adaptive procedure.Two numerical examples are given to demonstrate the accuracy and efficiency of the proposed method and the results show that the higher order AS T-spline based NMM shows higher accuracies when solving both continuous and discontinuous problems.Furthermore, the local mesh refinement using analysis-suitable T-spline reduces the number of degrees of freedom while maintaining calculation accuracy at the same time.