Abstract: Slope-stabilizing piles typically are commonly permitted significant horizontal displacement while satisfying the requirement of bearing capacity. The in-situ pile tests of such piles indicate that with minimal thrust and pile top displacements less than 10mm, cracking occurs near the sliding surface of the pile body, leading to non-elastic flexural deformations. Despite this, slope-stabilizing piles are still conventionally treated as elastic in the current design calculations and numerical simulations, resulting in significant discrepancies between expected and actual behaviors. To address such issues, using tested piles as a case study, a concrete pile model with actual reinforcement was established using the Diana finite element program. This model incorporates material nonlinearities such as the total strain crack model, the Von-Mises model, and the hardening soil model to realistically simulate the behaviour of the concrete pile, steel bars, and soils, respectively. The analysis further considered the boundary and geometric nonlinearities inherent in pile-soil and soil-rock interactions. Results show high agreement with measured data for displacements at the pile top or along the pile body; the bending moments in sections of the pile near the top and bottom, where no cracks appeared, aligned well with measurements. The corresponding load that led to the advent of cracking, and the positions of cracking also highly consistent with experimental observations. For the first time, numerical simulations revealed that as thrust increases, a double semi- "inverted cone" wedge-shaped shear failure zone forms in front of the anti-slide pile, consistent with experiment descriptions of triangular wedge deformations and soil shear failures leading to localized reductions in soil resistance. The phenomenon of the decrease of the soil resistances in some places induced by the shear failure of the soil mass was also coincident with the actual situation. The above results demonstrate that the proposed method can evidently enhance the level of the design calculation and analysis of slope-stabilizing piles and demonstrate potential for wider application.