Statistics of A combined incremental-secant mean-field homogenization scheme with per-phase residual strains for elasto-plastic composites