This study aims to investigate the tip resistance mechanism of open-ended steel pipe piles under partially plugged conditions by decomposing the load-sharing contribution of the ring zone and the internal soil core. A virtual static loading test was performed using the two-dimensional discrete element method (2D-DEM). Note that the findings of this study were obtained within the range of the 2D-DEM analysis conditions and do not intend to directly reproduce the three-dimensional arching mechanism or to establish equivalence between 2D and 3D responses. Quasi-static conditions were ensured by identifying loading parameters such that the energy residual remained ≤5% during driving, rest, and static loading phases, and the sensitivity criterion |Δqb|/qb ≤ 3% was satisfied when the loading rate was halved or doubled. The primary evaluation range of static loading was set to s/D = 0. 1 (10% D), corresponding to the displacement criterion for confirming the tip resistance in the Japanese design specifications for highway bridges. For reference, the post-peak mechanism was additionally tracked up to s/D = 0. 2 (20% D). Within a fixed evaluation window located immediately beneath the pile tip, high-contact-force (HCF) points were binarized using the threshold τ = μ + σ, and their occupancy ratio φ and normalized force intensity I* were calculated separately for the ring and core regions. A density-based contribution index (“K-density share”) was defined by combining “strength × area” and normalizing by the geometric width. The results suggest that, for the sand conditions and particle-scale ratios examined (D/d₅0 = 25–100), the ring zone tends to carry on the order of 85–90% of the tip resistance within the observed cases up to the ultimate state. Even at high plugging ratios (CRs), the internal soil core gradually increases its occupancy and intensity with settlement; however, high-contact-force struts beneath the ring remain active, and it is suggested that the ring-dominant load-transfer mechanism is generally preserved. In the post-peak plastic regime, the K-density share remains around 60%, indicating that the internal core plays a secondary, confining role rather than becoming dominant. These findings suggest that the conventional plug/unplug classification based on PLR can be supplemented by a combined use of plugging ratio CR (a kinematic indicator) and the ring contribution index (K-density share), potentially enabling a continuous interpretation of plugged and unplugged behaviors and contributing to the establishment of a design backbone for tip resistance evaluation. Calibration of design coefficients, scale regression, and mapping to practical indices such as N-values will be addressed in part II of this study. (Note: “Contribution” in this study refers to the HCF-based density contribution index K-density share, not the reaction–force ratio. )
Katayama et al. (Fri,) studied this question.