<p>In numerical modeling, the progressive failure behavior in the geogrid pullout test has not been effectively reproduced. Therefore, this study proposes a novel method for this modeling. This method proved successful, with the following findings: a) the numerical study exhibited remarkable agreement with the physical modeling; b) the maximum pullout force obtained corresponded to the initial yielding of the zone near the end of the geogrid, whereas the residual pullout force was obtained when all zones adjacent to the soil–geogrid interface reached cohesion at residual; c) a bubble of shear failure zones developed near the end of the geogrid, transitioning from the state of the maximum pullout force to the state of the residual pullout force.</p>