A numerical and analytical method for solving nonlinear geomechanics problems has been developed. This method is based on the combination of the Vlasov variational method, the Newton-Kantorovich linearization method and the method of general iterative process. The physical nonlinearity of the rock deformation process is considered in the framework of the main hypotheses of small elastic and plastic deformation Genki-Ilyushin.