Many material and structure failures are initiated at the interface between two distinctive materials, such as the debonding of interfaces of matrix/fiber of composites or bulk matrix/inclusion in metals at micro-scales. Accurate and robust modeling of these interface failures is critical to the design of novel material and structure with improved functionality. This paper presents a robust implementation of a discontinuous Galerkin finite element method in MOOSE based on the Incomplete Interior Penalty Galerkin scheme for modeling material interface behavior. This DG-based interface method naturally accommodates important extrinsic cohesive zone material models and is capable of modeling the failures of material interfaces across which material properties are highly distinctive. We demonstrate a good performance of the proposed method using single fiber and multiple fiber debonding problems.
