Geogrid has been extensively used in geotechnical engineering practice due to its effectiveness and economy. Deep insight into the interaction between the backfill soil and the geogrid is of great importance for proper design and construction of geogrid reinforced earth structures. Based on the calibrated model of sand and geogrid, a series of numerical pullout tests are conducted using PFC3D under special considerations of particle angularity and aperture geometry of the geogrid. In this work, interface characteristics regarding the displacement and contact force developed among particles and the deformation and force distribution along the geogrid are all visualized with PFC3D simulations so that new understanding on how geogrid-soil interaction develops under pullout loads can be obtained. Meanwhile, a new variable named fabric anisotropy coefficient is introduced to evaluate the inherent relationship between macroscopic strength and microscopic fabric anisotropy. A correlation analysis is adopted to compare the accuracy between the newly-proposed coefficient and the most commonly used one. Furthermore, additional pullout tests on geogrid with four different joint protrusion heights have been conducted to investigate what extent and how vertical reinforcement elements may result in reinforcement effects from perspectives of bearing resistance contribution, energy dissipation, as well as volumetric response. Numerical results show that both the magnitude and the directional variation of normal contact forces gov the development of macroscopic strength and the reinforcing effects of joint protrusion height can be attributed to the accelerated energy dissipation across the particle assembly and the intensive mobilization of the geogrid.