In vivo aspiration experiments on human livers are analyzed and material parameters for a non-linear-viscoelastic constitutive model are determined. A novel procedure is applied for the inverse analysis that accounts for the initial tissue deformation in the experiment and for the non-homogeneity of liver tissue. A numerical model is used consisting of a surface layer (capsule) and an underlying non-linear-viscoelastic solid (parenchyma). The capsule is modeled as hyperelastic membrane using data from measurements on bovine and human tissue. In a two step optimization procedure the set of constitutive model parameters for the "average" response of liver parenchyma is obtained. The proposed model is in line with literature values of high strain rate elastic modulus obtained from dynamic elastography. The model can be used to predict the nonlinear, time dependent behavior of human liver in computer simulations related to surgery training and planning.