El presente trabajo de investigación se centra en el estudio tanto experimental como numérico de la respuesta de suelos no saturados ensayados en un aparato triaxial verdadero con succión controlada. En primer lugar, se realizó una revisión bibliográfica en el estado del conocimiento de los principales equipos utilizados en la actualidad para estudiar el efecto de la tensión intermedia en caminos a succión controlada y se escogió el modelo numérico para interpretar los datos. Los ensayos se modelaron como problemas de contorno mediante el programa Elementos Finitos acoplado (termo-)hidro-mecánico Code_bright, que considera una implementación visco-plástica del modelo constitutivo básico de Barcelona (BBM). Para realizar la parte experimental se ha utilizado un equipo de laboratorio llamado triaxial cúbico desarrollado en de la Universidad de Texas at Arlington, entre cuyas características tiene la particularidad de permitir el control simultáneo de los esfuerzos principales (esfuerzos principales menores, intermedio y mayor), las deformaciones principales y la aplicación y control del estado de succión matricial en el suelo a ensayar. En una primera etapa la muestra cúbica se llevó a un estado de esfuerzos bajo equilibrio hidrostático (HC). A continuación, la muestra se sometió a una trayectoria de esfuerzos de compresión confinada (CTC) hasta lograr un estado crítico incipiente, inmediatamente seguido por una descarga monotónica. Posteriormente, la misma muestra se volvió a equilibrar hidrostáticamente (HC) bajo un mayor esfuerzo isotrópico y, finalmente, sometida a una trayectoria de esfuerzos de compresión con tensión neta media constante (TC). En todo el transcurso del ensayo, la succión matricial se mantuvo constante (s=50 kPa) mediante la técnica de traslación de ejes. Finalizada la etapa experimental se realizó la modelación numérica con la cual se pretendió evaluar la capacidad del modelo numérico para predecir la respuesta de una muestra de suelo bajo diferentes trayectorias de esfuerzos aplicados con succión controlada. Para este propósito se modeló la respuesta obtenida en el laboratorio tanto para la presente investigación como para otros resultados de diferentes autores obtenidos en la literatura. La respuesta general de las comparaciones entre el modelo numérico y los ensayos experimentales, ha llegado a una aproximación razonable. Las discrepancias observadas entre las predicciones del modelo y los resultados de los ensayos podrían deberse a la anisotropía inducida por la compactación estática en condiciones de deformación lateral nula.