Abstract

A new enriched finite element technique, named the Discontinuity-Enriched Finite Element Method (DE-FEM), was introduced recently for solving problems with both weak and strong discontinuities in 2-D. In this mesh-independent procedure, enriched degrees of freedom are added to new nodes collocated at the intersections between discontinuities and the sides of finite elements of the background mesh. In this work we extend DE-FEM to 3-D and describe in detail the implementation of a geometric engine capable of handling interactions between discontinuities and the background mesh. Several numerical examples in linear elastic fracture mechanics demonstrate the capability and performance of DE-FEM in handling discontinuities in a fully mesh-independent manner. We compare convergence properties and the ability to extract stress intensity factors with standard FEM. Most importantly, we show DE-FEM provides a stable formulation with regard to the condition number of the resulting system stiffness matrix.

Highlights

  • Enriched finite element methods have fundamentally changed the modeling of problems containing discontinuities

  • We show Discontinuity-Enriched Finite Element Method (DE-FEM) provides a stable formulation with regard to the condition number of the resulting system stiffness matrix. ⃝c 2019 The Authors

  • In this paper we extend DE-FEM to 3-D, and we discuss in detail its computer implementation in a displacement-based finite element code

Read more

Summary

Introduction

Enriched finite element methods have fundamentally changed the modeling of problems containing discontinuities. S.J. van den Boom, F. van Keulen et al / Computer Methods in Applied Mechanics and Engineering 355 (2019) 1097–1123 a single formulation to model problems containing both weak and strong discontinuities (C−1-continuous field) by placing enriched degrees of freedom (DOFs) only to nodes created along discontinuities. Creating matching meshes in general can be quite challenging depending on the morphology of the problem, because of the strict conditions on mesh quality that are required—elements with bad aspect ratios could reduce the approximation accuracy [4]; robustness is still questionable, especially for 3-D FE meshes, as issues remain in generating meshes that correctly match boundaries [5] Numerical techniques such as Universal Meshes [6] and the Conforming to Interface Structured Adaptive Mesh Refinement (CISAMR) [7] have been proposed to modify meshes locally to discontinuities. We showcase DE-FEM in the challenging problem of intersecting weak and strong discontinuities, which is handled effortlessly by means of the hierarchical implementation ala HIFEM

Formulation
Geometric engine
Hierarchical bookkeeping
Enrichment functions
Single discontinuity
Multiple discontinuities
Convergence
Stress intensity factors
Stability
Summary and conclusions
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call