A phase field approach enhanced with a cohesive zone model for modeling delamination induced by matrix cracking