Failure Surface Morphology and Uplift Bearing Capacity Calculation of Soil Surrounding Glass Fiber-Reinforced Plastic Screw Anchors

[Objective] As a new type of support anchor, the Glass Fiber-Reinforced Plastic (GFRP) screw anchor is increasingly used in engineering practice due to its simple construction, light weight, high strength, corrosion resistance, and cost-effectiveness. Previous calculation formulas for uplift bearing capacity of screw anchors mostly based on artificially defined boundaries between shallow and deep embedment, resulting in large differences in critical embedment depth ratio H/D (where H is the embedment depth and D is the anchor plate diameter) given by different scholars. This study aims to propose a calculation formula for uplift bearing capacity based on the generalized unified failure surface morphology, providing theoretical support for the design of GFRP screw anchors in engineering practice. [Methods] Using finite element numerical simulation software ABAQUS, numerical simulations of the uplift performance of GFRP screw anchors under different embedment ratios were conducted. The maximum uplift bearing capacity of vertical GFRP screw anchors in soil and the evolution of the failure surface morphology of anchored soil were investigated. [Results and Conclusion] (1) The load-displacement curves of GFRP screw anchor during uplift in anchored soil could be divided into three stages: the elastic stage, the local plastic stage, and the penetration failure stage. In the elastic stage, the load and displacement showed linear relationship, with the slope of the curve representing the equivalent stiffness of the anchor plate-soil system. In the local plastic stage, plastic deformation occurred in the local soil zone, deformation of the anchored soil gradually increased, system stiffness decreased, load-displacement relationship became non-linear (indicated by a gradually decreasing curve slope), and the interaction between the GFRP screw anchor and the anchored soil began to exhibit non-linear characteristics. In the penetration failure stage, the resistance of the anchored soil reached its ultimate value, cracks in the anchored soil became interconnected, the soil failure surface formed, and the curve entered a stable phase where the load no longer changed with displacement increase. The uplift bearing capacity factor Nγ showed a trend of initial increase followed by gradual stabilization with embedment ratio, which matched well with experimental results from other scholars, thereby verifying the simulation’s validity. The turning point occurred near H/D=9, suggesting that an embedment ratio H/D=9 could be considered as the critical embedment ratio distinguishing shallow from deep embedment. (2) Numerical simulation showed that during uplift, the anchored soil could be divided into three zones based on stress state: conical active zone formed by compaction above the GFRP screw anchor, passive zone extending outward along screw anchor edge influenced by extrusion from soil within failure surface, and transition zone between active and passive zones in plastic state. The soil within the failure surface consistently extruded the soil outside failure surface. Except for the outermost region, which might be in a state of at-rest earth pressure, the soil on the failure surface was generally in a state of passive earth pressure. With increasing embedment ratio, the failure surface morphology of the anchored soil underwent a dynamic evolution from “trumpet” shape➝“goblet” shape➝“light bulb” shape. No distinct “critical embedment ratio” between shallow and deep embedment for GFRP screw anchors existed. (3) By integrating the maximum uplift bearing capacities obtained from numerical simulations and the evolution patterns of the failure morphology in the anchored soil, a generalized unified failure surface model that did not require distinguishing between shallow and deep embedment for the soil anchored by GFRP screw anchors was proposed. This model was derived by introducing a cubic term based on the “trumpet” shape (with inclined linear boundaries) failure surface. On this basis, a calculation formula for the uplift bearing capacity under different embedment depths of GFRP screw anchors was derived, and the formula was compared and validated with numerous experimental results from domestic and international studies. The results demonstrated that the proposed calculation formula could effectively predict the maximum uplift bearing capacity under different embedment depth ratios.