The nonlinear variation of soil compressibility and permeability with void ratio (i.e., e-log σ′ and e-log k) has been included in the consolidation theory to accurately predict the behavior of soft soil stabilized by vertical drains. However, most current nonlinear consolidation models incorporating the coupled radial-vertical flow are based on some simplified assumptions, while including some features such as the complex implementation of multilayered computations, time-dependent loading and stress distribution with depth. This study hence introduces a novel approach where the spectral method is used to analyze the nonlinear consolidation behavior of multilayered soil associated with coupled vertical-radial drainage. In addition, time- and depth-dependent stress and soil properties at each soil layer are incorporated into the proposed model. Subsequently, the solution is verified against experimental and field data with comparison to previous analytical solutions. The results show greater accuracy of the proposed method in predicting in-situ soil behavior. A parametric study based on the proposed solution indicates that the ratio between the compression and permeability indices (ω = Cc/Ck) has a great impact on the consolidation rate, i.e., the greater the ω, the smaller the consolidation rate. Increasing the load increment ratio and the absolute difference between unity and ω (i.e., |ω − 1|) can exacerbate prediction error if the conventional simplified methods are used.