This study is concerned with plasticity and ductile fracture modelling aspects of press hardened steels, where multiple phases of different yield strength, work hardening and ductility are present. Macroscopic constitutive properties are determined using a semi-analytical approach termed Mean Field Homogenization (MFH), based on the properties of individual phases and micro-topology. This is combined with a phenomenological ductile fracture initiation criterion formulated in stress space at the microscopic scale. Localization enhanced kinematics are introduced at the macro-scale when instability is signaled by a local bifurcation analysis of the homogenized material. Discontinuities in the strain field are introduced into a quadrilateral shell based on the weak discontinuity approach. Fracture properties and flow curves up to large strains are determined from single phase tensile experiments using digital image correlation. The predictions of the integrated MFH and fracture model is compared with experimental results on steel sheets with different thermal histories, in terms of fracture prediction and and overall mechanical properties