Geotechnical and geological problems involve the description of the behavior of materials such as soil and rock, and their eventual interaction with fluids and structures. In general, the evolution of these problems is characterized by large deformations and displacements, discontinuities, heterogeneities and complex constitutive behavior. Addressing these problems requires numerical techniques that take these characteristics into account, without numerical drawbacks associated with element distortion as occurs in the finite element method (FEM). In this thesis is developed a computational algorithm based on the material point method (MPM) to approximate the solution of the governing equations to the mentioned phenomena. The algorithm is based on a three-dimensional dynamic formulation of the continuum considering large deformations. Rayleigh damping and non-viscous local damping are incorporated to model dynamic and quasi-static problems. The dynamic generation of pore pressures is formulated assuming the saturated porous medium and a single material point to discretize the mixture. Different techniques are evaluated to mitigate spurious pressure in impact problems on saturated media. Different constitutive models are implemented to model the failure surface and the soil mass flow process during slope instabilities, as well as the genesis, evolution and failure zone quantification in geological processes. To address the discretization of large-scale geological problems using MPM, a methodology is proposed and validated with the discretization of the Daguangbao landslide, in China. In order to decrease the computational time, the algorithm is implemented according to the parallel programming paradigm.