The nano-scaled structural heterogeneity is closely associated with the magnetic softness and mechanical properties of Fe-based amorphous alloys. However, limited by the shape and size of ribbons, the details of the atomic-scale structural evolution in ribbons under external stimulation are still missing. In this work, Maxwell-Voigt model with 2 K units and 1 Mx unit was employed to investigate the structural response of FeBCCu amorphous ribbons by substitution of P for B during creeping process under nanoindentation. The P-added alloy exhibited a more sensitive mechanical response at various loading rates, especially for the reduction of hardness and the increase of maximum indentation depth and creep displacement. With the help of creep analysis using above Maxwell-Voigt model, both types of defects were inclined to be activated at a quasi-static loading mode in P-added alloy, and the small size defects seemed to be easily activated and finally coupled into large defects in the fast loading mode, revealing the atomic diffusion and reconstruction process of activated defects. These types of defects evolved into the response of structural heterogeneity to external stimulation, manifesting a large deformation behavior during creeping. This work enlightens the comprehensive insights into the evolution of structural defects in Fe-based amorphous ribbons and has certain guiding significance for the structural design of high-performance alloys.