Numerical simulations of fracture propagation provide critical insights into mechanical processes, enabling early failure prediction to inform design decisions and optimise material development. Despite progress in damage modelling over 30 years, it remains a complex computational challenge, particularly in ductile fracture, where failure follows substantial elastoplastic deformations. Phase-field models have emerged as the state-of-the-art approach and are widely used in mesh-based schemes such as the Finite Element Method [1]. However, large elastoplastic deformations often lead to mesh distortion errors, compromising the accuracy of simulations with mesh-based numerical schemes. To address this limitation, a Material Point Method (MPM) is coupled with phase-field models to mitigate mesh distortion effects. Unlike conventional deforming mesh approaches, MPM employs a particle-based representation of the deformable body within a fixed background grid, preserving continuum approximation while reducing reliance on high particle densities. This study utilises an Extended B-Spline-Based MPM scheme [2] for the first time, which not only addresses cell-crossing and numerical integration errors associated with MPM but also enhances computational accuracy and facilitates the representation of discrete crack surfaces. The effectiveness of the proposed method is demonstrated through benchmark simulations, illustrating its advantages over traditional mesh-based approaches [1] Yin B., Kaliske M., (2020). A ductile phase-field model based on degrading the fracture toughness: Theory and implementation at small strain, Computer Methods in Applied Mechanics and Engineering, 366, 113068, DOI: https://doi.org/10.1016/j.cma.2020.113068. [2] Yamaguchi Y., Moriguchi S., Terada K., (2001). Extended B-Spline-based implicit material point method, Internatinal Journal for Numerical Methods in Engineering, 122, pp. 1746-1769, DOI: https://doi.org/10.1002/nme.6598.