Summary
A recently proposed phase‐field model for cohesive fracture is examined. Previous investigations have shown stress oscillations to occur when using unstructured meshes. It is now shown that the use of nonuniform rational B‐splines (NURBS) as basis functions rather than traditional Lagrange polynomials significantly reduces this oscillatory behavior. Moreover, there is no effect on the global structural behavior, as evidenced through load‐displacement curves. The phase‐field model imposes restrictions on the interpolation order of the NURBS used for the three different fields: displacement, phase field, and crack opening. This holds within the Bézier element, but also at the boundaries, where a reduction to 𝒞0‐continuity yields optimal results. Application to a range of cases, including debonding of a hard fiber embedded in a soft matrix, illustrates the potential of the cohesive phase‐field model.