A methodology for predicting 3D crack propagation in quasi-brittle materials, such as rocks and concrete, is developed to overcome mesh dependency in finite element method (FEM) simulations. Building on the theory of configurational mechanics and extending previous 2D work , the approach is based on a strategy to reorient crack fronts towards a minimum energy configuration. The formulation incorporates fundamental FEM concepts with zero-thickness interface elements using an elastoplastic fracture mechanics model. A MATLAB code has been developed that extends the existing 2D strategy to 3D, incorporating hexahedral and tetrahedral elements, calculation of configurational forces and re-orientation of element edges on the basis of the direction of those forces. Numerical simulations, including beam bending tests with varying notch configurations, confirm the capability of the method to improve fracture predictions by minimizing system energy. The research contributes an initial methodology for 3D crack front reorientation with potential applications in geotechnical and materials engineering.