An FEM model of two-layer steel cylinder with internal hydraulic pressure was established, and the influences of different contact analysis methods on the calculation precision were analyzed systematically by MARC software. Based on the results, it can be concluded that both the discrete and analytical methods in MARC software can be effectively used in contact analysis and the calculation precision can be reached when the mesh density of two contact bodies are in consistent during the numerical simulation. The error of boundary inaccurate description can be eliminated and the calculation precision can be improved by analytical method when subjected to inconsistent mesh density. Based on above analysis, the stress distribution in the spherical backing of a hydraulic cylinder of some forging press was calculated, which proved that analytical method could be used to improve calculation precision. Research work provides an effective method to carry out contact analysis when the mesh density is not inconsistent for the existence of great structure difference between two bodies.