Three‐dimensional finite elements with embedded strong discontinuities to model material failure in the infinitesimal range

This paper presents new three‐dimensional finite elements with embedded strong discontinuities in the small strain infinitesimal range. The goal is to model localized surfaces of failure in solids, such as cracks at fracture, through enhancements of the finite elements that capture the propagating discontinuities of the displacement field in the element interiors. In this way, such surfaces of discontinuity can be sharply resolved in general meshes not necessarily related to the detailed geometry of the surface, unknown a priori. An important issue is also the consideration of general finite element formulations in the developments (e.g., basic displacement‐based, mixed or enhanced assumed strain finite element formulations), as needed to optimally resolve the continuum problem in the bulk. The actual modeling of the discontinuity effects, including the incorporation of the cohesive law defining the discontinuity constitutive response, is carried out at the element level with the proper enhancement of the discrete strain field of the element. The added elemental degrees of freedom approximate the displacement jumps associated with the discontinuity and are defined independently from element to element, thus allowing their static condensation at the element level without affecting the global mechanical problem in terms of the number and topology of the global degrees of freedom. In fact, this global‐local structure of the finite element methods developed in this work arises naturally from a multi‐scale characterization of these localized solutions, with the discontinuities understood to appear in the small scales, thus leading directly to these computationally efficient numerical methods for their numerical resolution, easily incorporated to an existing finite element code. The focus in this paper is on the development of finite elements incorporating a linear interpolation of the displacement jumps in the general three‐dimensional setting. These interpolations are shown to be necessary for hexahedral elements to avoid the so‐called stress locking that occurs with simpler constant approximations of the jumps (namely, a spurious transfer of stresses across the discontinuity not allowing its full release and, hence, resulting in an overstiff or locked numerical solution). The design of the new finite elements is accomplished in this work by a direct identification of the separation modes to be incorporated in the discrete strain field of the element, rather than from an assumed discontinuous interpolation of the displacements, assuring with this approach their locking‐free response by design. An additional issue addressed in the paper is the geometric characterization and propagation of the discontinuity surfaces in the general three‐dimensional setting of interest here. The paper includes a series of numerical simulations illustrating and evaluating the properties of the new finite elements. Copyright © 2012 John Wiley & Sons, Ltd.