Effects of disorder on electron-doped iron pnictides are investigated systematically based on self-consistent Bogoliubov-de Gennes equations. Multiply impurities with same scattering potential (SP) are randomly distributed in a square lattice. Probability distribution functions of normalized order parameters for different impurity concentrations δ imp , different electron doping concentrations δ are investigated for given SPs. Samples are found to be very robust against weak SP, in which order parameters do not have qualitative change even at very large δ imp . While strong SP is able to easily break down the order parameters. For moderate SP, variations of order parameters on and around impurities strongly depend on δ, however the distribution functions of normalized order parameters have similar behavior as δ imp increases. Compared with superconducting (SC) order, the magnetic order is more sensitive to multi-impurity effect. The spatial spin density wave pattern has already been destroyed before the system loses its superconductivity. Dependence of SC order on temperature is similar to that of impurity-free case, with the critical temperature being remarkably suppressed for high δ imp .