In the present study Finite Element Method (FEM) has been applied for modelling of indentation process in polycrystalline metals. Two different models of the material, with linear and non-linear hardening, have been assumed in calculations. The indentation process has been quantitatively described by load-depth curves. The hardness values were obtained from the unloading part of such curves, which are responsible for plastic deformation of the material during indentation process. In addition, the geometry changes of the indent after deformation have been analyzed. The results show that simulation techniques, used in the present study, are useful in estimation of materials hardness. It has been also found that specific combination of materials properties such as hardening exponent, hardening coefficient, tangent modulus and yield stress may result in sink-in or pile-up phenomenon. These processes may influence the indent geometry and hardness measurements consequently.