The formation cross section of the nuclide production of proton-induced reactions with intermediate energy is important for a variety of applications. For instance, the mass and charge distribution of residual products produced in the spallation reactions needs to be studied because it can provide useful information for the disposal of nuclear waste and residual radioactivity generated by the spallation neutron target system. With our current work, we have developed the Many Stage Dynamical Model (MSDM) based on the Cascade-Exciton Model (CEM). By introducing Mshnik's recent work on the CEM code, the MSDM code and the Quantum Molecular Dynamics (QMD) plus Statistical Decay Model (SDM) (QMD+SDM) and QMD plus FISSION (QMD+FISSION) models are adopted; we use them to investigate the mass distribution of Nb, Au, and Pb of proton-induced reactions with energies from 100 MeV to 3 GeV. The agreement between the developed MSDM simulations and the measured data as well as the QMD+FISSION model are good in the energy range of 100 MeV to 3 GeV, and deviations mainly show up in the mass range of 90 to 140 in the high energy of protons incident on the Au and Pb target for the MSDM and QMD+FISSION model simulations. The QMD+SDM can reproduce only part of the spallation fragments and cannot reproduce the fission fragments of the measured data.