Optimizing integration point density for exponential finite element shape functions for phase-field modeling of fracture in functionally graded materials

Abstract

In the realm of fracture mechanics and phase-field modeling, the utilization of exponential finite element (EFE) shape functions has shown promise in accurately predicting fracture responses, particularly in scenarios with intricate crack propagation paths. However, the computational complexity associated with EFE shape functions, including higher integration schemes and orientation requirements, prompts the need for optimization. This study delves into the critical aspect of integration point density and aims to determine the ideal balance between accuracy and computational efficiency in the context of EFE shape functions. The research focuses on functionally graded materials, which inherently exhibit spatial variations in solution quantities due to their graded material system. Through a systematic investigation, the study aims to identify the optimal number of integration points required to maximize the computational advantages offered by EFE shape functions while minimizing the associated computational resources. By conducting a comprehensive analysis across various loading scenarios, such as tension and shear, this study intends to establish a quantitative relationship between integration point density and prediction accuracy for fracture responses.