Three dimensional calculations of ductile fracture under mode I plane strain, small scale yielding conditions are carried out using an elastic-viscoplastic constitutive relation for a progressively cavitating solid with two populations of void nucleating second phase particles. Larger inclusions that result in void nucleation at an early stage are modeled discretely while smaller particles that require large strains to nucleate voids are homogeneously distributed. Full field solutions are obtained for eight volume fractions, ranging from 1% to 19%, of randomly distributed larger inclusions. For each volume fraction calculations are carried out for seven random distributions of inclusion centers. Crack growth resistance curves and fracture surface roughness statistics are calculated using standard procedures. The crack growth resistance is characterized in terms of both JIC and the tearing modulus TR. For all volume fractions considered, the computed fracture surfaces are self-affine over a size range of nearly two orders of magnitude with a microstructure independent roughness exponent of 0.53 with a standard error of 0.0023. The cut-off length of the scale invariant regime is found to depend on the inclusion volume fraction. Consideration of the full statistics of the fracture surface roughness revealed other parameters that vary with inclusion volume fraction. For smaller values of the discretely modeled inclusion volume fraction (≤7%), there is a linear correlation between several measures of fracture surface roughness and both JIC and TR. In this regime crack growth is dominated by a void-by-void process. For greater values of the discretely modeled inclusion volume fraction, crack growth mainly involves multiple void interactions and no such correlation is found. © 2013 Elsevier Ltd.