A nodal C1C1 approximation procedure combined with partition of unity enrichment has been constructed for fracture. The double-interpolation approximation improves the accuracy of the solution field by up to one order of magnitude with the same degrees of freedom as the constant strain triangular element. The double-interpolation approximation is able to calculate the stress at node directly through interpolation. The enriched form of the double-interpolation approximation enables a way to propagate the crack with minimal remeshing with a smoother crack tip stress field obtained.