La modélisation proposée, dans cet article, comporte de nombreuses originalités, notamment l'écoulement de l'eau dans un espace poral inter-agrégats avec une loi d'écoulement seuil et le couplage gonflement-consolidation. Le but de ce développement est de montrer la faisabilité d'une telle modélisation qui comportera, dans le cas d'une approximation par éléments finis, cinq variables nodales W1, W2, W3 (déplacements), u (pression interstitielle) et Sr (degré de saturation).