Pile-up around indenter is usually observed during instrumented indentation tests on bulk metallic glass. Neglecting the pile-up effect may lead to errors in evaluating hardness,Young’s modulus,stress-strain response,etc. Finite element analysis was employed to implement numerical simulation of spherical indentation tests on bulk metallic glass. A new model was proposed to describe the pile-up effect. By using this new model,the contact radius and hardness of Zr41.2Ti13.8Cu12.5Ni10Be22.5 bulk metallic glass were obtained under several different indenter loads with pile-up,and the results agree well with the data generated by numerical simulation.