Abstract: The purpose of this thesis is to investigate different numerical models evaluating large strain consolidation and creep and to provide assessment of potential impacts of these models to engineering designs. Several models were used to conduct calculations presented in this thesis: CONDES, iConsol.js, PLAXIS and a finite difference code developed in Mathematica. The model developed in Mathematica was used to investigate differences between the classical (linear) and the non-linear consolidation equations with the goal to identify limitations of the classical approach. In addition, a simplified creep model developed in Mathematica was used to evaluate differences between iConsol.js and PLAXIS, two selected models that are available to practicing engineers for the assessment of secondary compression. Results of this study indicate that a model selection should be tailored to specific design requirements after considering results of the available geotechnical investigations and monitoring because all computer codes considered in this study exhibit limitations in terms of the available constitutive models, the ability to account for secondary compression effects, and the flexibility to account for changes in the boundary condition, i.e. model results may not be able to represent actual soil behavior when considering general (pre-defined) material, boundary and loading conditions.