Abstract: We establish the a priori convergence rate for fnite element approximations of a class of nonlocal nonlinear fracture models. We consider state-based peridynamic models where the force at a material point is due to both the strain between two points and the change in volume inside the domain of the nonlocal interaction. The pairwise interactions between points are mediated by a bond potential of multi-well type while multi-point interactions are associated with the volume change mediated by a hydrostatic strain potential. The hydrostatic potential can either be a quadratic function, delivering a linear force-strain relation, or a multi-well type that can be associated with the material degradation and cavitation. We frst show the well-posedness of the peridynamic formulation and that peridynamic evolutions exist in the Sobolev space H². We show that the fnite element approximations converge to the H² solutions uniformly as measured in the mean square norm. For linear continuous fnite elements, the convergence rate is shown to be C_tΔt + C_sh²∕ε², where ε is the size of the horizon, h is the mesh size, and Δ_t is the size of the time step. The constants C_t and C_s are independent of Δt and h and may depend on ε through the norm of the exact solution. We demonstrate the stability of the semi-discrete approximation. The stability of the fully discrete approximation is shown for the linearized peridynamic force. We present numerical simulations with the dynamic crack propagation that support the theoretical convergence rate.

Key words: Nonlocal fracture models, Peridynamic, State-based peridynamic, Numerical analysis, Finite element approximation