“…Analytical solutions for the consolidation of a semi-infinite poroelastic body and a single or multi-layered saturated poroelastic material system have been acquired Gibson, 1960a, 1960b;Schiffman and Fungaroli, 1965;Booker and Small, 1982a, 1982bVardoulakis and Harnpattanapanich, 1986;Yue et al, 1994;Yue and Selvadurai, 1995;Senjuntichai and Rajapakse, 1995;Pan, 1999). An efficient analytical algorithm to solve Biot's consolidation of a multi-layered saturated poroelastic material system is the transfer matrix method (Bahar, 1972;Fang, 2001, 2003;Ai et al, 2008aAi et al, , 2008bAi et al, , 2010bAi et al, , 2010cAi and Cheng, 2009a). However, it is difficult to eliminate the positive exponential functions involved in the elements of the transfer matrix, which may lead to unstable numerical results or overflow in the calculation, due to the influence of numerically illconditioned matrices.…”