Local refinement with arbitrary irregular meshes and implementation in numerical manifold method

Abstract The solution accuracy based on finite element interpolation strongly depends on mesh quality. Thus, the numerical manifold method advocates that if a finite element mesh is selected as the mathematical cover, the mesh should be uniform, where all elements are congruent squares or isosceles right triangles for two-dimensional problems. However, the elements close to the points of stress concentrations or singularities must be smaller than those further away from these points. This leads to an irregular mesh in which a few nodes of certain small elements are hung on the edges of adjacent larger elements. This study implements local refinement with arbitrary k-irregular meshes by selecting piecewise polynomials as the local approximations over physical patches and enforcing appropriate constraints on the piecewise polynomials. All the elements are ordinary isoparametric elements, and elements with variable numbers of nodes on edges are not required. No linear dependency issue exists. Furthermore, refinement with any number of mesh layers is developed so that solution error can be distributed as uniformly as possible. Typical examples with different singularities are designed to illustrate the effectiveness of the proposed procedure.

[1]  Kenjiro Terada,et al.  Performance assessment of generalized elements in the finite cover method , 2004 .

[2]  S. Timoshenko,et al.  Theory of elasticity , 1975 .

[3]  Peng Zhang,et al.  Explicit edge-based smoothed numerical manifold method for transient dynamic modeling of two-dimensional stationary cracks , 2021, Engineering Analysis with Boundary Elements.

[4]  S. Im,et al.  Variable-node element families for mesh connection and adaptive mesh computation , 2012 .

[5]  A. Gupta A finite element for transition from a fine to a coarse grid , 1978 .

[6]  Charles E. Augarde,et al.  A new partition of unity finite element free from the linear dependence problem and possessing the delta property , 2010 .

[7]  Thomas-Peter Fries,et al.  Hanging nodes and XFEM , 2011 .

[8]  I. Babuska,et al.  The design and analysis of the Generalized Finite Element Method , 2000 .

[9]  I. Babuska,et al.  Introduction to Finite Element Analysis: Formulation, Verification and Validation , 2011 .

[10]  Zhiye Zhao,et al.  Investigation of linear dependence problem of three-dimensional partition of unity-based finite element methods , 2012 .

[11]  Hyun Gyu Kim Interface element method: Treatment of non-matching nodes at the ends of interfaces between partitioned domains , 2003 .

[12]  J. Fish The s-version of the finite element method , 1992 .

[13]  Massimo Ruzzene,et al.  A geometric multiscale finite element method for the dynamic analysis of heterogeneous solids , 2013 .

[14]  H. H. Zhang,et al.  Extract of stress intensity factors on honeycomb elements by the numerical manifold method , 2012 .

[15]  Jae Hyuk Lim,et al.  Variable‐node elements for non‐matching meshes by means of MLS (moving least‐square) scheme , 2007 .

[16]  Timon Rabczuk,et al.  Effective 2D and 3D crack propagation with local mesh refinement and the screened Poisson equation , 2017 .

[17]  T. Rabczuk,et al.  XLME interpolants, a seamless bridge between XFEM and enriched meshless methods , 2014 .

[18]  Ted Belytschko,et al.  Cracking particles: a simplified meshfree method for arbitrary evolving cracks , 2004 .

[19]  Liu Zhi-jun,et al.  Two-dimensional numerical manifold method with multilayer covers , 2015, Science China Technological Sciences.

[20]  M. Wheeler,et al.  Multigrid on the interface for mortar mixed finite element methods for elliptic problems , 2000 .

[21]  S. Y. Han,et al.  Determination of T-stress for thermal cracks in homogeneous and functionally graded materials with the numerical manifold method , 2021 .

[22]  F. Liu,et al.  Two-dimensional Hermitian numerical manifold method , 2020, Computers & Structures.

[23]  Hong Zheng,et al.  New strategies for some issues of numerical manifold method in simulation of crack propagation , 2014 .

[24]  Shuodao Wang,et al.  A Mixed-Mode Crack Analysis of Isotropic Solids Using Conservation Laws of Elasticity , 1980 .

[25]  Ted Belytschko,et al.  Elastic crack growth in finite elements with minimal remeshing , 1999 .

[26]  Yongtao Yang,et al.  Smoothed numerical manifold method with physical patch‐based smoothing domains for linear elasticity , 2020, International Journal for Numerical Methods in Engineering.

[27]  T. Belytschko,et al.  A three dimensional large deformation meshfree method for arbitrary evolving cracks , 2007 .

[28]  Timon Rabczuk,et al.  Dual‐horizon peridynamics , 2015, 1506.05146.

[29]  Timon Rabczuk,et al.  Element-wise fracture algorithm based on rotation of edges , 2013 .

[30]  H. Zheng,et al.  A domain decomposition based method for two-dimensional linear elastic fractures , 2016 .

[31]  T. Laursen,et al.  A mortar‐finite element formulation for frictional contact problems , 2000 .

[32]  Timon Rabczuk,et al.  Adaptive fourth-order phase field analysis for brittle fracture , 2020 .

[33]  T. Belytschko,et al.  Element‐free Galerkin methods , 1994 .

[34]  H. Nguyen-Xuan,et al.  An extended isogeometric thin shell analysis based on Kirchhoff-Love theory , 2015 .

[35]  Guowei Ma,et al.  Numerical analysis of 2-D crack propagation problems using the numerical manifold method , 2010 .

[36]  Gen-Hua Shi,et al.  Manifold Method of Material Analysis , 1992 .

[37]  Ted Belytschko,et al.  A method for dynamic crack and shear band propagation with phantom nodes , 2006 .

[38]  Timon Rabczuk,et al.  Damage and fracture algorithm using the screened Poisson equation and local remeshing , 2016 .

[39]  Hiroshi Tada,et al.  The stress analysis of cracks handbook , 2000 .

[40]  Luís A.G. Bitencourt,et al.  A coupling technique for non-matching finite element meshes , 2015 .

[41]  Klaus-Jürgen Bathe,et al.  The finite element method enriched by interpolation covers , 2013 .

[42]  Feng Liu,et al.  Complementarity problem arising from static growth of multiple cracks and MLS-based numerical manifold method , 2015 .

[43]  Genki Yagawa,et al.  Non‐matching mesh gluing by meshless interpolation—an alternative to Lagrange multipliers , 2007 .

[44]  S. M. Liu,et al.  Computation of T-stresses for multiple-branched and intersecting cracks with the numerical manifold method , 2019, Engineering Analysis with Boundary Elements.

[45]  Timon Rabczuk,et al.  Dual-horizon peridynamics: A stable solution to varying horizons , 2017, 1703.05910.

[46]  Feng Liu,et al.  The MLS-based numerical manifold method with applications to crack analysis , 2014, International Journal of Fracture.

[47]  Faker Ben Belgacem,et al.  The Mortar finite element method with Lagrange multipliers , 1999, Numerische Mathematik.