A phase field model for ductile fracture is coupled with the Modified Mohr Coulomb (MMC) model for plastic damage evolution and subsequent crack growth. An energy-based damage threshold is applied to control degradation due to ductile damage. The model is implemented through a user subroutine. MMC parameters from the literature are utilized and found to be compatible with the model, accurately reproducing material response curves in a variety of loading conditions for 6016-T4 aluminum alloy. The influence of model parameters is demonstrated and additionally the Nakazima test is simulated to demonstrate the capability of the model predicting the formability of the material through a failure locus. The model is found capable of reproducing experimentally observed crack paths and quantitative material behavior.