A numerical model is presented for one-dimensional large strain consolidation of a layered soil stratum. The algorithm accounts for vertical strain, soil self-weight, conventional constitutive relationships, changing material properties during consolidation, unload/reload, time-dependent loading and boundary conditions, an externally applied hydraulic gradient, and multiple soil layers with different material properties. The model can also accommodate depth-dependent loading and variable preconsolidation stress profiles. Verification checks of the model show excellent agreement with available analytical and numerical solutions for consolidation of soils with multiple layers and depth-dependent loading. An overview of this model is presented and a parametric study is performed to illustrate the effects of large strain and depth-dependent loading on consolidation response of layered soils. The simulations indicate that large errors may occur when stress reduction with depth is not correctly accounted for.