The phase field method integrates the Griffith theory and damage mechanics approach to predict crack initiation, propagation, and branching within one framework. No crack tracking topology is needed, and complex crack shapes can be captures without user intervention. In this paper, a detailed description of how the phase field method is implemented with explicit dynamics into LS-DYNA is provided. The displacement field and the damage field are solved in a staggered approach and the phase field equation is solved every Nth time step (N is refered to as calculation cycle) to save computational time. An N value smaller than 1/400 of the total time step numbers is suggested. Several simulations are presented to demonstrate the feasibility of this solving scheme.