The discrete element method (DEM) coupled with the pore-scale finite volume (PFV) method was used to simulate the suffusion and post-suffusion behavior of gap-graded soil with different initial fines content (fc). A series of drained triaxial tests were performed on the non-eroded, eroded and reconstituted samples. The results indicate that the erosion ratio (Er) increases as the fc of the sample increases. The mechanical response of the sample with an initial fc of 15 % is almost unaffected by suffusion. The dilatancy and peak deviatoric stress ratios of sample with initial fc of 25 % and 35 % are significantly lower with increasing erosion ratio. When the Er is greater than 8.2, eroded samples with an initial fc of 35 % collapse during shearing, and the eroded and reconstituted samples behave as dilation and contraction, respectively. As the initial fc increases, the mechanical response of the reconstituted samples differs more and more from that of the corresponding eroded samples. The force chain analysis indicates more force chains in the eroded sample than in the reconstituted sample, resulting in higher dilatancy of the eroded sample with fc of 25 % and 35 % during shearing.