In this paper, a three-dimensional part model is constructed for the finite element analysis of hemming models where hemming defects frequently occur. The roller path is modeled as the boundary condition with the one-dimensional beam element and the revolute joint model. With the constructed part model and the roller movement, a finite element analysis has been pursued in order to identify the hemming load and hemming defects such as wrinkling in the flange region. The analysis result shows that the maximum hemming load occurs in the intake situation while oscillatory behavior of the load is found especially when hemming the curved model because of wrinkle generation. This paper compares the amplitude and the period of wrinkle between the analysis result and the experiment, which shows good agreement with each other.