Crack propagation is a typical discontinuous problem, which is difficult to be simulated with traditional finite element method. In order to solve the defects of traditional methods, extended finite element method introduced level set function method to describe crack shape so as to realize initiation and development of crack in the inner element. The problem of crack propagation simulation in traditional method replying on the life and death of element was then solved. This paper used Python langrage to develop a crack automatically propagation program package in ABAQUS platform. I mode edge crack of a finite flat and the crack propagation in three point bending beam were simulated. Compared to the exact results, this method can effectively simulate the process of crack propagation. The propagation of multiple cracks was simulated at the same time. The results showed that the method in this paper not only took advantage of ABAQUS software to simulate complex loading and boundary, but also reflected program package’s superiority on crack analysis. This method provided a convenient and effective way to solve real complex problems.