As camadas argilosas saturadas submetidas a um carregamento vertical sofrem deformações verticais ao longo do tempo, isto é, adensamento. Essas deformações quando acentuadas provocam danos estruturais, funcionais ou estéticos nas edificações. Nesse sentido, um modelo matemático é desenvolvido para a previsão do adensamento das argilas saturadas levando-se em conta os conceitos da mecânica dos solos. A formulação matemática considera a existência de uma parcela viscosa da tensão efetiva vertical e incorpora na equação diferencial parcial (EDP) modelada uma lei de viscosidade não linear como uma propriedade fundamental. A EDP do adensamento é obtida em termos de dissipação do excesso de poropressão normalizada. Em vista da complexidade do modelo, o método das diferenças finitas (MDF) é utilizado para a obtenção da solução teórica considerando variação no tempo e no espaço. A solução numérica mostra que a parcela viscosa da tensão efetiva vertical atua retardando o fenômeno do adensamento. Além disso, o ajuste dos valores teóricos aos experimentais mostra que a modelagem para o fator viscosidade igual a um é promissora.