The aim of the work is to perform a comparative analysis of the results of analyzing arbitrarily loaded shells of revolution using finite element method in various formulations, namely, in the formulation of the displacement method and in the mixed formulation. Methods. To obtain the stiffness matrix of a finite element a functional based on the equality of the actual work of external and internal forces was applied. To obtain the deformation matrix in the mixed formulation the functional obtained from the previous one by replacing the actual work of internal forces in it with the difference of the total and additional work was used. Results. In the formulation of the displacement method for an eight-node hexahedral solid finite element, displacements and their first derivatives are taken as the nodal unknowns. Approximation of the displacements of the inner point of the finite element was carried out through the nodal unknowns on the basis of the Hermite polynomials of the third degree. For a finite element in the mixed formulation, displacements and stresses were taken as nodal unknowns. Approximation of the target finite element values through their nodal values in the mixed formulation was carried out on the basis of trilinear functions. It is shown on a test example that a finite element in the mixed formulation improves the accuracy of the strength parameters of the shell of revolution stress-strain state.