Mode I fracture in a nonlinear lattice with viscoelastic forces

We study Mode I fracture in a viscoelastic lattice model with a nonlinear force law, with a focus on the velocity and linear stability of the steady-state propagating solution. This study is a continuation both of the study of the piece-wise linear model in Mode I, and the study of more general nonlinear force laws in Mode III fracture. At small driving, there is a strong dependency of the velocity curve on the dissipation and a strong sensitivity to the smoothness of the force law at large dissipation. At large driving we calculate, via a linear stability analysis, the critical velocity for the onset of instability as a function of the smoothness, the dissipation and the ratio of lattice spacing to critical extension. This critical velocity is seen to be very sensitive to these parameters. We confirm our calculations via direct numerical simulations of the initial value problem.