 the following modifications have been added to the paper:  In section Finite Element Simulation, paragraph 0 "The contact between the probe and the phantom is  modeled with two surface contact pairs covering the palmar region of the probe and the upper layer of the  phantom, respectively. The pressure contact calculation is based on an Augmented Lagrangian Method,  and the friction between the 0 contacts is modeled as Coulomb friction (