SUMMARYNumerical crack propagation schemes were augmented in an elegant manner by the X-FEM method. The use of special tip enrichment functions, as well as a discontinuous function along the sides of the crack allows one to do a complete crack analysis virtually without modifying the underlying mesh, which is of industrial interest, especially when a numerical model for crack propagation is desired. This paper improves the implementation of the X-FEM method for stress analysis around cracks in three ways. First, the enrichment strategy is revisited. The conventional approach uses a 'topological' enrichment (only the elements touching the front are enriched). We suggest a 'geometrical' enrichment in which a given domain size is enriched. The improvements obtained with this enrichment are discussed. Second, the conditioning of the X-FEM both for topological and geometrical enrichments is studied. A preconditioner is introduced so that 'off the shelf' iterative solver packages can be used and perform as well on X-FEM matrices as on standard FEM matrices. The preconditioner uses a local (nodal) Cholesky based decomposition. Third, the numerical integration scheme to build the X-FEM stiffness matrix is dramatically improved for tip enrichment functions by the use of an ad hoc integration scheme. A 2D benchmark problem is designed to show the improvements and the robustness.
SUMMARYThis paper is devoted to the imposition of Dirichlet-type conditions within the extended finite element method (X-FEM). This method allows one to easily model surfaces of discontinuity or domain boundaries on a mesh not necessarily conforming to these surfaces. Imposing Neumann boundary conditions on boundaries running through the elements is straightforward and does preserve the optimal rate of convergence of the background mesh (observed numerically in earlier papers). On the contrary, much less work has been devoted to Dirichlet boundary conditions for the X-FEM (or the limiting case of stiff boundary conditions). In this paper, we introduce a strategy to impose Dirichlet boundary conditions while preserving the optimal rate of convergence. The key aspect is the construction of the correct Lagrange multiplier space on the boundary. As an application, we suggest to use this new approach to impose precisely zero pressure on the moving resin front in resin transfer moulding (RTM) process while avoiding remeshing. The case of inner conditions is also discussed as well as two important practical cases: material interfaces and phase-transformation front capturing.
International audienceThis paper introduces a new algorithm to define a stable Lagrange multiplier space to impose stiff interface conditions within the context of the extended finite element method. In contrast to earlier approaches. we do not work with an interior penalty formulation as, e.g. for Nitsche techniques, but impose the constraints weakly in terms of Lagrange multipliers. Roughly speaking a stable and optimal discrete Lagrange multiplier space has to satisfy two criteria: a best approximation property and a uniform inf-sup condition. Owing to the fact that the interface does not match the edges of the mesh, the choice of a good discrete Lagrange Multiplier space is not trivial. Here we propose a new algorithm for the local construction of the Lagrange Multiplier space and show that a uniform inf-sup condition is satisfied. A counterexample is also presented, i.e. the inf-sup constant depends on the mesh-size and degenerates as it tends to zero. Numerical results in two-dimensional confirm the theoretical ones. Copyrigh
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.