La mecánica de suelos no saturados es una ciencia relativamente joven que se ha venido desarrollado en las últimas décadas. Los temas de investigación en esta área son muy amplios. El desarrollo de modelos constitutivos para suelos se desarrolló primero para suelos saturados. Este tipo de suelos son bifásicos y su comportamiento es más simple en comparación con los suelos no saturados. Sin embargo, estos últimos suelos, son mucho más comunes que los primeros. Para analizar el comportamiento de las grandes obras de ingeniería civil, se requieren modelos constitutivos generales. Algunas características importantes de estos modelos es que son simples y requieren una pequeña cantidad de parámetros con significado físico. En este trabajo, se desarrolla un modelo basado en la teoría del estado crítico y la ecuación de esfuerzos efectivos propuesta por Bishop. El valor del parámetro χ se obtiene mediante un modelo hidráulico propuesto por Zhou et al. (2012) que permite incluir el fenómeno de histéresis de la curva de retención en los suelos. El objetivo de la investigación es desarrollar un modelo simple con una matriz de rigidez simétrica que requiera pocos parámetros que puedan obtenerse mediante pruebas de laboratorio. Se desarrolló un programa de computadora del modelo acoplado en un código para Matlab. Se incluyen comparaciones entre resultados experimentales y numéricos.