Crack propagation is modeled by means of manifold method with first order cover function employed. Quadratic displacement functions are formed on manifold element in this situation. An algorithm to deal with the crack tip’s stopping at any place of an element is presented. Thereby
the length of crack extension in each extension step can be controlled by means of the algorithm
and the shortcoming in crack simulation by manifold method with constant cover function is overcome.