Pile penetration in soft ground involves complex mechanisms, including significant alterations in the soil state surrounding the pile, which influence the pile negative skin friction (NSF) over time. However, the pile penetration process is often excluded from finite element analysis. This paper investigates the impact of pile penetration on the generation of NSF and dragload. A stable node-based smoothed particle finite element method (SNS-PFEM) framework is introduced for two-dimensional axisymmetric conditions and coupled consolidation, incorporating the ANICREEP model of soft soil with a modified cutting-plane algorithm. A field case study with penetration process is simulated to verify the numerical model's performance, followed by a parametric analysis on the effect of penetration rate on NSF during consolidation. Results indicate that without the pile penetration process in NSF analysis can result in an unsafely low estimation of NSF and dragload magnitudes. The penetration rate affects dragload only at the initial consolidation stage. As consolidation progresses, dragload converges to nearly the same magnitude across different rates. Additionally, current design methods inadequately predict the beta value (where beta is an empirical factor correlating vertical effective stress of soil with the pile skin friction) and its time dependency, for which a new empirical formula for the time-dependent beta value is proposed and successfully applied to other field cases.
