A numerical method for crack growth is described in which the crack is not regarded as a single discontinuity that propagates continuously. Instead, the crack is represented by a set of overlapping cohesive segments. These cohesive segments are inserted into finite elements as discontinuities in the displacement field by exploiting the partition-of-unity property of shape functions. The cohesive segments can be incorporated at arbitrary locations and orientations and are not tied to any particular mesh direction. The evolution of decohesion of the segments is governed by a cohesive law. The independent specification of bulk and cohesive constitutive relations leads to a characteristic length being introduced into the formulation. The formul...