This paper presents numerical studies of stiff fiber pullout behaviors of fiber reinforced cementitious composites based on a progressive damage model. The ongoing debonding process is simulated. Interfacial stress distribution for different load levels is analyzed. A parametric study, including bond strength and the homogeneity index on the pullout behaviors is carried out. The numerical results indicate that the bond stress decreases gradually from loaded end to embedded end along fiber-cement interface. The debonding initially starts from loaded end and propagates to embedded end as load increasing. The embedded length and bond strength affect the load-loaded end displacement curves significantly. The numerical results have a general agreement with the experimental investigation.