Generalmente, la evaluación de la falla no drenada por subpresión en lumbreras construidas en suelos se lleva a cabo mediante la determinación de un factor de seguridad utilizando expresiones analíticas. Sin embargo, existen escasos lineamientos en la literatura para la simulación numérica de este mecanismo. En este trabajo se exploran tres metodologías para la evaluación numérica de la falla por subpresión, que consideran diferentes procedimientos para la generación del desequilibrio que deriva en la condición de falla. Los resultados se comparan con el factor de seguridad obtenido a partir de una expresión analítica. Se señala la metodología más adecuada para representar el mecanismo de falla por subpresión, y se realizan observaciones sobre la posible necesidad de una revisión adicional a la estabilidad de la lumbrera a causa de la disminución de la altura del lodo bentonítico.