$hp$-adaptive finite element simulation of a static anti-plane shear crack in a nonlinear strain-limiting elastic solid

TL;DR

This paper develops an $hp$-adaptive finite element method to analyze static anti-plane shear cracks in nonlinear strain-limiting elastic solids, combining residual-based mesh refinement and local polynomial degree adjustment for improved accuracy.

Contribution

It introduces a novel $hp$-adaptive finite element scheme for nonlinear strain-limiting materials with a dual-error estimator, enhancing simulation accuracy for crack problems.

Findings

Method demonstrates high accuracy and convergence in numerical experiments.

Regularized crack-tip fields are characterized for various parameters.

Framework supports extension to dynamic crack propagation problems.

Abstract

An $hp$-adaptive continuous Galerkin finite element method is developed to analyze a static anti-plane shear crack embedded in a nonlinear, strain-limiting elastic body. The geometrically linear material is described by a constitutive law relating stress and strain that is algebraically nonlinear. In this investigation, the constitutive relation utilized is \textit{uniformly bounded}, \textit{monotone}, \textit{coercive}, and \textit{Lipschitz continuous}, ensuring the well-posedness of the mathematical model. The governing equation, derived from the balance of linear momentum coupled with the nonlinear constitutive relationship, is formulated as a second-order quasi-linear elliptic partial differential equation. For a body with an edge crack, this governing equation is augmented with a classical traction-free boundary condition on the crack faces. An $hp$-adaptive finite element scheme is proposed for the numerical approximation of the resulting boundary value problem. The adaptive strategy is driven by a dual-component error estimation scheme: mesh refinement ($h$-adaptivity) is guided by a residual-based a posteriori error indicator of the \textit{Kelly type}, while the local polynomial degree ($p$-adaptivity) is adjusted based on an estimator of the local solution regularity. The performance, accuracy, and convergence characteristics of the proposed method are demonstrated through numerical experiments. The structure of the regularized crack-tip fields is examined for various modeling parameters. Furthermore, the presented framework establishes a robust foundation for extension to more complex and computationally demanding problems, including quasi-static and dynamic crack propagation in brittle materials.