The phase field method for fracture integrates the Griffith theory and damage mechanics approach to predict crack initiation and propagation within one framework. It replaced the discrete representation of crack by diffusive damage and solved it based on a minimization of the global energy storage functional. As a result, no crack tracking topology is needed, and complex crack shapes can be captures without user intervention. However, it is also reported to have an inconsistency between the predicted fracture toughness and the material strength. Recently, Wu's research  has proposed a novel energetic degradation function to handle this issue. This research does some further modifications to the global energy storage functional so that Newton's method can be directly used to solve the energy minimization. With the new energy form, direct implementation of the length-scale independent phase field method into finite element packages like LS-DYNA becomes possible. This paper presents the framework and details of implementing the length-scale independent phase field method into LS-DYNA through a user-defined element and material subroutine. Several numerical examples are presented to compare with the experiment crack shape. Most importantly, this paper is one of the first ones to quantitatively predict accurate force response compared to experiments. These examples verify the accuracy of the new energy form and implementation.