En este artículo, desarrollamos un modelo elasto-viscoplástico (EVP) para arcilla utilizando la regla de flujo no asociada. Esto se logra mediante el uso de una forma modificada de la teoría EVP sobreeestressada de Perzyna, la mecánica de suelos en estado crítico y la teoría de múltiples superficies. El nuevo modelo incluye seis parámetros, cinco de los cuales son idénticos a los del modelo de mecánica de suelos en estado crítico. El otro parámetro es el índice de compresión secundaria no lineal generalizado. El modelo EVP se implementó en un código consolidado acoplado no lineal utilizando un algoritmo numérico de elementos finitos (AFENA). Luego, probamos el modelo para diferentes arcillas, como la arcilla de Osaka, la arcilla de lodo de la Bahía de San Francisco, la arcilla de caolín y la arcilla de depósito marino de Hong Kong. Los resultados numéricos muestran una buena concordancia con los datos experimentales.