• PDF
• Cite
• Share
Article Contents  Article Contents

# Finite element approximation of nonlocal dynamic fracture models

• * Corresponding author: P. K. Jha

This material is based upon work supported by the U. S. Army Research Laboratory and the U. S. Army Research Office under contract/grant number W911NF1610456

• In this work we estimate the convergence rate for time stepping schemes applied to nonlocal dynamic fracture modeling. Here we use the nonlocal formulation given by the bond based peridynamic equation of motion. We begin by establishing the existence of $H^2$ peridynamic solutions over any finite time interval. For this model the gradients can become large and steep slopes appear and localize when the non-locality of the model tends to zero. In this treatment spatial approximation by finite elements are used. We consider the central-difference scheme for time discretization and linear finite elements for discretization in the spatial variable. The fully discrete scheme is shown to converge to the actual $H^2$ solution in the mean square norm at the rate $C_t\Delta t +C_s h^2/\epsilon^2$. Here $h$ is the mesh size, $\epsilon$ is the length scale of nonlocal interaction and $\Delta t$ is the time step. The constants $C_t$ and $C_s$ are independent of $\Delta t$, and $h$. In the absence of nonlinearity a CFL like condition for the energy stability of the central difference time discretization scheme is developed. As an example we consider Plexiglass and compute constants in the a-priori error bound.

Mathematics Subject Classification: Primary: 34A34, 34B10, 74S05; Secondary: 74H55.

 Citation: • • Figure 1.  Two-point potential $W^\epsilon(S,y - x)$ as a function of strain $S$ for fixed $y - x$

Figure 2.  Nonlocal force $\partial_S W^\epsilon(S,y - x)$ as a function of strain $S$ for fixed $y - x$. Second derivative of $W^\epsilon(S,y-x)$ is zero at $\pm \bar{r}/\sqrt{|y -x|}$

Figure 3.  Geometry

• ## Article Metrics  DownLoad:  Full-Size Img  PowerPoint