The objective of this course is to familiarize the followers to the geotechnical constitutive laws and the applications of numerical methods to the analysis of boundary value problems in geomechanics. In particular, the course provides discussions of actual material response vs. predictions based on the results of implemented numerical model. It also illustrates the concepts with practical applications in geotechnical engineering. The student will develop their plasticity theory based finite element code, written in Mathcad.