foam-extend4.1-coherent-io/applications/solvers/solidMechanics/viscoElasticSolidFoam/updateCrack.H

{
    forAll (DU.boundaryField(), patchI)
    {
        if
        (
            DU.boundaryField()[patchI].type()
         == cohesiveZoneIncrementalModeIIFvPatchVectorField::typeName
        )
        {
            cohesiveZoneIncrementalModeIIFvPatchVectorField& cohesivePatchDU =
                refCast<cohesiveZoneIncrementalModeIIFvPatchVectorField>
                (
                    DU.boundaryField()[patchI]
                );

            nCrackedFaces = cohesivePatchDU.updateCrack();
            Info << "nCrackedFaces is " << nCrackedFaces << endl;
        }
    }
}