A computer-implemented press forming simulation method for determining a springback-inducing area of an elastic-plastic material being formed with a forming tool, the press forming simulation method comprising: calculating an element equivalent nodal force vector {f} e from a stress tensor σ representing stress in the elastic-plastic material at a bottom dead point of the forming tool using a finite element method for one or a plurality of finite elements of a target configuration of the elastic-plastic material based on the following equation f e = ∫ Ve L t B t σ dV = : N x i N y i N z i M x i M y i M z i : where [L] represents a coordinate transform matrix and [B] represents a relationship matrix between displacement and strain, the coordinate transform matrix [L] is a transformation matrix for converting the element equivalent nodal force vector in an element coordinate system into a vector in a global coordinate system, the coordinate transform matrix [L] is calculated by using a direction cosine with respect to an element coordinate system (X', Y' and Z') of the X-, Y- and Z axes of the global coordinate system; calculating a total equivalent nodal force vector {f} of areas by integrating a plurality of the element equivalent nodal force vectors {f} e calculated for one or a plurality of finite elements over all of the areas or specified areas of the elastic-plastic material as represented in the following equation f = ∫ V f e dV the total equivalent nodal force vector consisting of cross-sectional force (Nx, Ny, Nz) and bending moment (Mx, My, Mz); and specifying positions of areas with an amount of bending moment (Mx, My, Mz) exceeding a predetermined threshold as springback-inducing areas regarding the bending moment (Mx, My, Mz) of components of the calculated total equivalent nodal force vector. calculating an element equivalent nodal force vector {f} e from a stress tensor σ representing stress in the elastic-plastic material at a bottom dead point of the forming tool using a finite element method for one or a plurality of finite elements of a target configuration of the elastic-plastic material based on the following equation f e = ∫ Ve L t B t σ dV = : N x i N y i N z i M x i M y i M z i : where [L] represents a coordinate transform matrix and [B] represents a relationship matrix between displacement and strain, the coordinate transform matrix [L] is a transformation matrix for converting the element equivalent nodal force vector in an element coordinate system into a vector in a global coordinate system, the coordinate transform matrix [L] is calculated by using a direction cosine with respect to an element coordinate system (X', Y' and Z') of the X-, Y- and Z axes of the global coordinate system; calculating a total equivalent nodal force vector {f} of areas by integrating a plurality of the element equivalent nodal force vectors {f} e calculated for one or a plurality of finite elements over all of the areas or specified areas of the elastic-plastic material as represented in the following equation f = ∫ V f e dV the total equivalent nodal force vector consisting of cross-sectional force (Nx, Ny, Nz) and bending moment (Mx, My, Mz); and specifying positions of areas with an amount of bending moment (Mx, My, Mz) exceeding a predetermined threshold as springback-inducing areas regarding the bending moment (Mx, My, Mz) of components of the calculated total equivalent nodal force vector. <EOS>