Abstract: According to a reasonable hypothesis, the rolling deformation zone of H-beams was divided into two regions and the velocity field in each deformation region was set,lp by the flow function method. A rolling energetic parameter model of H-beams based on the flow function method was established. Parameters of deformation regions were optimized with the Powell method to minimize the value of total power. Finally, the rolling energetic parameters of H-beams were obtained. The Gauss integral method was used to get more accurate results. Calculation results show that errors between the results of the model and the finite element model verified by test do not exceed 1.53% when the elongations of the web and flange are the same. If the hypothesis is not satisfied, the model in this paper also has well calculation accuracy with proper modification. When the ratio of elongation between flange and webA is around 1, the trend of the model's rolling force is the same as the results of the finite element model. When the arm coefficients are reasonable, the rolling force of the model matches with the results of the finite element model well.