This paper presents the implementation of the Blatz-Ko hyperelastic compressible model in a finite element program to deal with large deformation problems. We show analytically and numerically that the minimum number of increment steps in the Newton-Raphson algorithm depends on material properties and applied loads. We also show that this dependence is related to the orientation preservation principle. So we propose a convergence criteria based on the sign of eigenvalues of the deformation gradient matrix.