Higher‐order extended FEM for weak discontinuities – level set representation,quadrature and application to magneto‐mechanical problems

In this article, we present the application of bilinear and biquadratic extended FEM (XFEM) formulations to model weak discontinuities in magnetic and coupled magneto‐mechanical boundary value problems. For properly resolving the location of curved interfaces and the discontinuous physical behaviour, the major part of the contribution is devoted to review and develop methods for level set representation of curved interfaces and numerical integration of the weak form in higher‐order XFEM formulations. In order to reduce the complexity of the representation of curved interfaces, an element local approach that allows for an automated computation of the level set values and also improves the compatibility between the level set representation and the integration subdomains is proposed. Integration rules for polygons and strain smoothing are applied in conjunction with biquadratic elements and compared with curved integration subdomains. Eventually, a coupled magneto‐mechanical demonstration problem is modelled and solved by XFEM. For demonstration purposes, a magneto‐mechanical coupling due to magnetic stresses is considered. Errors and convergence rates are analysed for the different level set representations and numerical integration procedures as well as their dependence on the ratio of material parameters at an interface. Copyright © 2012 John Wiley & Sons, Ltd.     

On time integration in the XFEM

The extended finite element method (XFEM) is often used in applications that involve moving interfaces. Examples are the propagation of cracks or the movement of interfaces in two‐phase problems. This work focuses on time integration in the XFEM. The performance of the discontinuous Galerkin method in time (space–time finite elements (FEs)) and time‐stepping schemes are analyzed by convergence studies for different model problems. It is shown that space–time FE achieve optimal convergence rates. Special care is required for time stepping in the XFEM due to the time dependence of the enrichment functions. In each time step, the enrichment functions have to be evaluated at different time levels. This has important consequences in the quadrature used for the integration of the weak form. A time‐stepping scheme that leads to optimal or only slightly sub‐optimal convergence rates is systematically constructed in this work. Copyright © 2009 John Wiley & Sons, Ltd.     

hp‐Generalized FEM and crack surface representation for non‐planar 3‐D cracks

A high‐order generalized finite element method (GFEM) for non‐planar three‐dimensional crack surfaces is presented. Discontinuous p‐hierarchical enrichment functions are applied to strongly graded tetrahedral meshes automatically created around crack fronts. The GFEM is able to model a crack arbitrarily located within a finite element (FE) mesh and thus the proposed method allows fully automated fracture analysis using an existing FE discretization without cracks. We also propose a crack surface representation that is independent of the underlying GFEM discretization and controlled only by the physics of the problem. The representation preserves continuity of the crack surface while being able to represent non‐planar, non‐smooth, crack surfaces inside of elements of any size. The proposed representation also provides support for the implementation of accurate, robust, and computationally efficient numerical integration of the weak form over elements cut by the crack surface. Numerical simulations using the proposed GFEM show high convergence rates of extracted stress intensity factors along non‐planar curved crack fronts and the robustness of the method. Copyright © 2008 John Wiley & Sons, Ltd.     

Extended finite element method coupled with face‐based strain smoothing technique for three‐dimensional fracture problems

In this work, the extended finite element method (XFEM) is for the first time coupled with face‐based strain‐smoothing technique to solve three‐dimensional fracture problems. This proposed method, which is called face‐based smoothed XFEM here, is expected to combine both the advantages of XFEM and strain‐smoothing technique. In XFEM, arbitrary crack geometry can be modeled and crack advance can be simulated without remeshing. Strain‐smoothing technique can eliminate the integration of singular term over the volume around the crack front, thanks to the transformation of volume integration into area integration. Special smoothing scheme is implemented in the crack front smoothing domain. Three examples are presented to test the accuracy, efficiency, and convergence rate of the face‐based smoothed XFEM. From the results, it is clear that smoothing technique can improve the performance of XFEM for three‐dimensional fracture problems. Copyright © 2015 John Wiley & Sons, Ltd.     

The XFEM for high‐gradient solutions in convection‐dominated problems

Convection‐dominated problems typically involve solutions with high gradients near the domain boundaries (boundary layers) or inside the domain (shocks). The approximation of such solutions by means of the standard finite element method requires stabilization in order to avoid spurious oscillations. However, accurate results may still require a mesh refinement near the high gradients. Herein, we investigate the extended finite element method (XFEM) with a new enrichment scheme that enables highly accurate results without stabilization or mesh refinement. A set of regularized Heaviside functions is used for the enrichment in the vicinity of the high gradients. Different linear and non‐linear problems in one and two dimensions are considered and show the ability of the proposed enrichment to capture arbitrary high gradients in the solutions. Copyright © 2009 John Wiley & Sons, Ltd.     

A partition of unity‐based ‘FE–Meshfree’ QUAD4 element for geometric non‐linear analysis

The recently published ‘FE–Meshfree’ QUAD4 element is extended to geometrical non‐linear analysis. The shape functions for this element are obtained by combining meshfree and finite element shape functions. The concept of partition of unity (PU) is employed for the purpose. The new shape functions inherit their higher order completeness properties from the meshfree shape functions and the mesh‐distortion tolerant compatibility properties from the finite element (FE) shape functions. Updated Lagrangian formulation is adopted for the non‐linear solution. Several numerical example problems are solved and the performance of the element is compared with that of the well‐known Q4, QM6 and Q8 elements. The results show that, for regular meshes, the performance of the element is comparable to that of QM6 and Q8 elements, and superior to that of Q4 element. For distorted meshes, the present element has better mesh‐distortion tolerance than Q4, QM6 and Q8 elements. Copyright © 2009 John Wiley & Sons, Ltd.     

c‐Type method of unified CAMG and FEA. Part 1: Beam and arch mega‐elements—3D linear and 2D non‐linear

Computer‐aided mesh generation (CAMG) dictated solely by the minimal key set of requirements of geometry, material, loading and support condition can produce ‘mega‐sized’, arbitrary‐shaped distorted elements. However, this may result in substantial cost saving and reduced bookkeeping for the subsequent finite element analysis (FEA) and reduced engineering manpower requirement for final quality assurance. A method, denoted as c‐type, has been proposed by constructively defining a finite element space whereby the above hurdles may be overcome with a minimal number of hyper‐sized elements. Bezier (and de Boor) control vectors are used as the generalized displacements and the Bernstein polynomials (and B‐splines) as the elemental basis functions. A concomitant idea of coerced parametry and inter‐element continuity on demand unifies modelling and finite element method. The c‐type method may introduce additional control, namely, an inter‐element continuity condition to the existing h‐type and p‐type methods. Adaptation of the c‐type method to existing commercial and general‐purpose computer programs based on a conventional displacement‐based finite element method is straightforward. The c‐type method with associated subdivision technique can be easily made into a hierarchic adaptive computer method with a suitable a posteriori error analysis. In this context, a summary of a geometrically exact non‐linear formulation for the two‐dimensional curved beams/arches is presented. Several beam problems ranging from truly three‐dimensional tortuous linear curved beams to geometrically extremely non‐linear two‐dimensional arches are solved to establish numerical efficiency of the method. Incremental Lagrangian curvilinear formulation may be extended to overcome rotational singularity in 3D geometric non‐linearity and to treat general material non‐linearity. Copyright © 2003 John Wiley & Sons, Ltd.     

An adaptive interface‐enriched generalized FEM for the treatment of problems with curved interfaces

An adaptive refinement scheme is presented to reduce the geometry discretization error and provide higher‐order enrichment functions for the interface‐enriched generalized FEM. The proposed method relies on the h‐adaptive and p‐adaptive refinement techniques to reduce the discrepancy between the exact and discretized geometries of curved material interfaces. A thorough discussion is provided on identifying the appropriate level of the refinement for curved interfaces based on the size of the elements of the background mesh. Varied techniques are then studied for selecting the quasi‐optimal location of interface nodes to obtain a more accurate approximation of the interface geometry. We also discuss different approaches for creating the integration sub‐elements and evaluating the corresponding enrichment functions together with their impact on the performance and computational cost of higher‐order enrichments. Several examples are presented to demonstrate the application of the adaptive interface‐enriched generalized FEM for modeling thermo‐mechanical problems with intricate geometries. The accuracy and convergence rate of the method are also studied in these example problems. Copyright © 2015 John Wiley & Sons, Ltd.     

Arbitrary placement of local meshes in a global mesh by the interface‐element method (IEM)

A new method is proposed to place local meshes in a global mesh with the aid of the interface‐element method (IEM). The interface‐elements use moving least‐square (MLS)‐based shape functions to join partitioned finite‐element domains with non‐matching interfaces. The supports of nodes are defined to satisfy the continuity condition on the interfaces by introducing pseudonodes on the boundaries of interface regions. Particularly, the weight functions of nodes on the boundaries of interface regions span only neighbouring nodes, ensuring that the resulting shape functions are identical to those of adjoining finite‐elements. The completeness of the shape functions of the interface‐elements up to the order of basis provides a reasonable transfer of strain fields through the non‐matching interfaces between partitioned domains. Taking these great advantages of the IEM, local meshes can be easily inserted at arbitrary places in a global mesh. Several numerical examples show the effectiveness of this technique for modelling of local regions in a global domain. Copyright © 2003 John Wiley & Sons, Ltd.     

Coupling of non‐conforming meshes in a component mode synthesis method

A common mesh refinement‐based coupling technique is embedded into a component mode synthesis method, Craig–Bampton. More specifically, a common mesh is generated between the non‐conforming interfaces of the coupled structures, and the compatibility constraints are enforced on that mesh via L2‐minimization. This new integrated method is suitable for structural dynamic analysis problems where the substructures may have non‐conforming curvilinear and/or surface interface meshes. That is, coupled substructures may have different element types such as shell, solid, and/or beam elements. The proposed method is implemented into a commercial finite element software, B2000++, and its demonstration is carried out using an academic and industry oriented test problems. Copyright © 2013 John Wiley & Sons, Ltd.     

A semi‐Lagrangean time‐integration approach for extended finite element methods

Many computational problems incorporate discontinuities that evolve in time. The eXtendend Finite Element Method (XFEM) is able to represent discontinuities sharply on fixed arbitrary meshes, but numerical difficulties arise if these discontinuities move in time. We point out that this issue is crucial for interface problems with strongly discontinuous fields on fixed grids. A method using semi‐Lagrangean techniques is proposed to adequately handle time integration based on finite difference schemes in the context of the XFEM. The basic idea is to adapt previous numerical solutions to the current interface position by tracking back virtual Lagrangean particles to their previous positions, where an appropriate solution can be extrapolated from a smooth field. Convergence properties of the proposed method in time and space are thoroughly studied for two one‐dimensional model problems. Finally, the method is applied to the particularly challenging problem of premixed combustion, where the discontinuity appears at the flame front separating the burnt from the unburnt gases. A two‐dimensional and a three‐dimensional expanding flame demonstrates that the method is sufficiently accurate to retain the properties of the overall Nitsche‐type formulation for interface problems with embedded strong discontinuities. Copyright © 2014 John Wiley & Sons, Ltd.     

An adaptive fixed mesh method for modelling and simulating of unsteady two‐phase flows with sharp physical interfaces

We present in this paper a new computational method for simulation of two‐phase flow problems with moving boundaries and sharp physical interfaces. An adaptive interface‐capturing technique (ICT) of the Eulerian type is developed for capturing the motion of the interfaces (free surfaces) in an unsteady flow state. The adaptive method is mainly based on the relative boundary conditions of the zero pressure head, at which the interface is corresponding to a free surface boundary. The definition of the free surface boundary condition is used as a marker for identifying the position of the interface (free surface) in the two‐phase flow problems. An initial‐value‐problem (IVP) partial differential equation (PDE) is derived from the dynamic conditions of the interface, and it is designed to govern the motion of the interface in time. In this adaptive technique, the Navier–Stokes equations written for two incompressible fluids together with the IVP are solved numerically over the flow domain. An adaptive mass conservation algorithm is constructed to govern the continuum of the fluid. The finite element method (FEM) is used for the spatial discretization and a fully coupled implicit time integration method is applied for the advancement in time. FE‐stabilization techniques are added to the standard formulation of the discretization, which possess good stability and accuracy properties for the numerical solution. The adaptive technique is tested in simulation of some numerical examples. With the test problems presented here, we demonstrated that the adaptive technique is a simple tool for modelling and computation of complex motion of sharp physical interfaces in convection–advection‐dominated flow problems. We also demonstrated that the IVP and the evolution of the interface function are coupled explicitly and implicitly to the system of the computed unknowns in the flow domain. Copyright © 2005 John Wiley & Sons, Ltd.     

A conformal decomposition finite element method for arbitrary discontinuities on moving interfaces

Interface capturing methods using enriched finite element formulations are well suited for solving multimaterial transport problems that contain weak or strong discontinuities. The conformal decomposition FEM decomposes multimaterial elements of a non‐conforming background mesh into sub‐elements that conform to material interfaces captured using a level set method. As the interface evolves, interfacial nodes move, and background nodes may change material. The present work describes approaches for handling moving interfaces in the context of the conformal decomposition FEM for both weakly and strongly discontinuous fields. Dynamic discretization methods using extrapolation and moving mesh approaches are considered and developed with first‐order and second‐order time integration methods. The moving mesh approach is demonstrated to be a stable method that preserves both weak and strong discontinuities on a variety of one‐dimensional and two‐dimensional test problems, while achieving the expected second‐order error convergence rate in space and time. Copyright © 2014 John Wiley & Sons, Ltd.     

A NURBS‐based interface‐enriched generalized finite element method for problems with complex discontinuous gradient fields

A non‐uniform rational B‐splines (NURBS)‐based interface‐enriched generalized finite element method is introduced to solve problems with complex discontinuous gradient fields observed in the structural and thermal analysis of the heterogeneous materials. The presented method utilizes generalized degrees of freedom and enrichment functions based on NURBS to capture the solution with non‐conforming meshes. A consistent method for the generation and application of the NURBS‐based enrichment functions is introduced. These enrichment functions offer various advantages including simplicity of the integration, possibility of different modes of local solution refinement, and ease of implementation. In addition, we show that these functions well capture weak discontinuities associated with highly curved material interfaces. The convergence, accuracy, and stability of the method in the solution of two‐dimensional elasto‐static problems are compared with the standard finite element scheme, showing improved accuracy. Finally, the performance of the method for solving problems with complex internal geometry is highlighted through a numerical example. Copyright © 2015 John Wiley & Sons, Ltd.     

Lower‐dimensional interface elements with local enrichment: application to coupled hydro‐mechanical problems in discretely fractured porous media

In this study, we develop lower‐dimensional interface elements to represent preexisting fractures in rock material, focusing on finite element analysis of coupled hydro‐mechanical problems in discrete fractures–porous media systems. The method adopts local enrichment approximations for a discontinuous displacement and a fracture relative displacement function. Multiple and intersected fractures can be treated with the new scheme. Moreover, the method requires less mesh dependencies for accurate finiteelement approximations compared with the conventional interface element method. In particular, for coupled problems, the method allows for the use of a single mesh for both mechanical and other related processes such as flow and transport. For verification purposes, several numerical examples are examined in detail. Application to a coupled hydro‐mechanical problem is demonstrated with fluid injection into a single fracture. The numerical examples prove that the proposed method produces results in strong agreement with reference solutions. Copyright © 2012 John Wiley & Sons, Ltd.     

An extended finite element library

This paper presents and exercises a general structure for an object‐oriented‐enriched finite element code. The programming environment provides a robust tool for extended finite element (XFEM) computations and a modular and extensible system. The programme structure has been designed to meet all natural requirements for modularity, extensibility, and robustness. To facilitate mesh–geometry interactions with hundreds of enrichment items, a mesh generator and mesh database are included. The salient features of the programme are: flexibility in the integration schemes (subtriangles, subquadrilaterals, independent near‐tip, and discontinuous quadrature rules); domain integral methods for homogeneous and bi‐material interface cracks arbitrarily oriented with respect to the mesh; geometry is described and updated by level sets, vector level sets or a standard method; standard and enriched approximations are independent; enrichment detection schemes: topological, geometrical, narrow‐band, etc.; multi‐material problem with an arbitrary number of interfaces and slip‐interfaces; non‐linear material models such as J2 plasticity with linear, isotropic and kinematic hardening. To illustrate the possible applications of our paradigm, we present 2D linear elastic fracture mechanics for hundreds of cracks with local near‐tip refinement, and crack propagation in two dimensions as well as complex 3D industrial problems. Copyright © 2007 John Wiley & Sons, Ltd.     

Parameter‐free,weak imposition of Dirichlet boundary conditions and coupling of trimmed and non‐conforming patches

We present a parameter‐free domain sewing approach for low‐order as well as high‐order finite elements. Its final form contains only primal unknowns; that is, the approach does not introduce additional unknowns at the interface. Additionally, it does not involve problem‐dependent parameters, which require an estimation. The presented approach is symmetry preserving; that is, the resulting discrete form of an elliptic equation will remain symmetric and positive definite. It preserves the order of the underlying discretization, and we demonstrate high‐order accuracy for problems of non‐matching discretizations concerning the mesh size h as well as the polynomial degree of the order of discretization p. We also demonstrate how the method may be used to model material interfaces, which may be curved and for which the interface does not coincide with the underlying mesh. This novel approach is presented in the context of the p‐version and B‐spline version of the finite cell method, an embedded domain method of high order, and compared with more classical methods such as the penalty method or Nitsche's method. Copyright © 2014 John Wiley & Sons, Ltd.     

hp‐adaptive extended finite element method

This paper discusses higher‐order extended finite element methods (XFEMs) obtained from the combination of the standard XFEM with higher‐order FEMs. Here, the focus is on the embedding of the latter into the partition of unity method, which is the basis of the XFEM. A priori error estimates are discussed, and numerical verification is given for three benchmark problems. Moreover, methodological aspects, which are necessary for hp‐adaptivity in XFEM and allow for exponential convergence rates, are summarized. In particular, the handling of hanging nodes via constrained approximation and an hp‐adaptive strategy are presented. Copyright © 2011 John Wiley & Sons, Ltd.     

On the application of higher‐order elements in the hierarchical interface‐enriched finite element method

This article introduces a new algorithm for evaluating enrichment functions in the higher‐order hierarchical interface‐enriched finite element method (HIFEM), which enables the fully mesh‐independent simulation of multiphase problems with intricate morphologies. The proposed hierarchical enrichment technique can accurately capture gradient discontinuities along materials interfaces that are in close proximity, in contact, and even intersecting with one another using nonconforming finite element meshes for discretizing the problem. We study different approaches for creating higher‐order HIFEM enrichments corresponding to six‐node triangular elements and analyze the advantages and shortcomings of each approach. The preferred method, which yields the lowest computational cost and highest accuracy, relies on a special mapping between the local and global coordinate systems for evaluating enrichment functions. A comprehensive convergence study is presented to show that this method yields similar convergence rate and precision as those of the standard FEM with conforming meshes. Finally, we demonstrate the application of the higher‐order HIFEM for simulating the thermal and deformation responses of several materials systems and engineering problems with complex geometries. Copyright © 2015 John Wiley & Sons, Ltd.     

A high‐order generalized FEM for through‐the‐thickness branched cracks

This paper presents high‐order implementations of a generalized finite element method for through‐the‐thickness three‐dimensional branched cracks. This approach can accurately represent discontinuities such as triple joints in polycrystalline materials and branched cracks, independently of the background finite element mesh. Representative problems are investigated to illustrate the accuracy of the method in combination with various discretizations and refinement strategies. The combination of local refinement at crack fronts and high‐order continuous and discontinuous enrichments proves to be an excellent combination which can deliver convergence rates close to that of problems with smooth solutions. Copyright © 2007 John Wiley & Sons, Ltd.