Rare molybdenum resources have been increasingly involved in heavy industries. In this paper, the common unidirectional and cross hot rolling operations, for pure molybdenum plate, are numerically simulated by using MSC. Marc software. An elastic-plastic finite element model is employed together with updated Lagrange method to predict stress and strain fields in the work-piece. The results showed that there was a typical three-dimensional additional compressive stress ( σy > σx > σz ) in deformation zone, while strain could be divided into uniaxial compressive strain and biaxial tensile strain ( Ey > Ex > Ez ). Tensile stress σx increased with the accumulation of reduction and the decrease of friction coefficient at the edge of width spread. More importantly, the interlaced deformation caused by cross commutation was helpful to repair the severe anisotropy created by unidirectional hot rolling. By comparing the theoretical verification of rolling forces and the measured temperatures with the simulated values, eventually, it is demonstrated that the model is aligning well with the actual engineering.