A rate-dependent model and its user subroutine for cohesive element method to investigate propagation and branching behavior of dynamic brittle crack