Abstract. We present a new method to simulate deformable volumetric objects interactively using finite elements. With quadratic basis functions and a non-linear strain tensor, we are able to model realistic local compression as well as large global deformation. The construction of the differential equations is described in detail including the Jacobian matrix required to solve the non-linear system. The results show that the bending of solids is reflected more realistically than with the linear refinement previously used in computer graphics. At the same time higher frame rates are achieved as the number of elements can be drastically reduced. Finally, an application to virtual tissue simulation is presented with the objective to improve surgical training. (a) (b) (c) Fig. 1: Interactive liver deformation, (a) 856 linear tetrahedra, (b) 96 quadratic tetrahedra, (c) a high resolution surface mesh attached.