Thaw consolidation of ice-rich permafrost is a typical problem in cold regions engineering. This paper proposes a three dimensional analysis of large strain thaw consolidation for post-thawed zone of permafrost, which is defined by a moving thawing boundary problem with phase changes. The theory is implemented in a numerical code and the numerical results are compared with thaw consolidation tests. For problems with low water contents, the small and large strain methods provide virtually the same results. For problems with high water contents, however, the large strain theory shows a much better performance.