<p>This study aims to clarify the influence of the variability in the material properties of Glass Fiber- Reinforced Polymers (GFRP) and the modeling method of GFRP on the slipping behavior of high- strength bolted frictional GFRP joints. Therefore, study conducted a statistical evaluation of the variability in material properties and finite element (FE) analysis. As a result, the initial load-relative displacement relationship in the anisotropic model generally matched the experimental results. Furthermore, it was found that the Poisson's ratio did not have a significant impact on the slipping behavior in any direction. The material property that most influenced the slipping behavior was considered to be the shear modulus in the plate thickness direction.</p>