A press-forming process for press-forming a workpiece with a molding tool, the process including a computer implemented method for identifying a shape defect factor, the method comprising: calculating a stress distribution that acts on an entire workpiece before elastic recovery; calculating a first deformation amount based on an elastic recovery of the workpiece in a certain evaluation point before and after the elastic recovery in a case where the stress distribution is applied to the workpiece before elastic recovery; calculating second deformation amounts based on the elastic recovery of the workpiece in the evaluation point before and after the elastic recovery in a case where the workpiece is divided into a plurality of regions set in advance, the stress distribution is divided into a plurality of regions corresponding to the divided regions of the workpiece, and each of the stress distributions divided into the regions is applied to the corresponding regions of the workpiece before elastic recovery, whereby each of the second deformation amounts related to the divided regions is calculated; characterized in that it further comprises the steps of comparing the first deformation amount with the second deformation amounts, and identifying a region having the smallest difference as a main influence region that is a stress distribution region most related to a shape defect of the workpiece in the evaluation point before and after the elastic recovery; calculating a third deformation amount based on the elastic recovery of the workpiece in the evaluation point before and after the elastic recovery in a case where the stress distribution acting on the main influence region is decomposed in each direction component of an orthogonal coordinate system and applied to the workpiece before elastic recovery; and comparing the second deformation amount in a case where the stress distribution of the main influence region is applied with the third deformation amount, and identifying, as the shape defect factor, the stress component in a direction having the smallest difference as a main influence stress component that is the stress component most related to the shape defect of the workpiece in the evaluation point before and after the elastic recovery. calculating a stress distribution that acts on an entire workpiece before elastic recovery; calculating a first deformation amount based on an elastic recovery of the workpiece in a certain evaluation point before and after the elastic recovery in a case where the stress distribution is applied to the workpiece before elastic recovery; calculating second deformation amounts based on the elastic recovery of the workpiece in the evaluation point before and after the elastic recovery in a case where the workpiece is divided into a plurality of regions set in advance, the stress distribution is divided into a plurality of regions corresponding to the divided regions of the workpiece, and each of the stress distributions divided into the regions is applied to the corresponding regions of the workpiece before elastic recovery, whereby each of the second deformation amounts related to the divided regions is calculated; characterized in that it further comprises the steps of comparing the first deformation amount with the second deformation amounts, and identifying a region having the smallest difference as a main influence region that is a stress distribution region most related to a shape defect of the workpiece in the evaluation point before and after the elastic recovery; calculating a third deformation amount based on the elastic recovery of the workpiece in the evaluation point before and after the elastic recovery in a case where the stress distribution acting on the main influence region is decomposed in each direction component of an orthogonal coordinate system and applied to the workpiece before elastic recovery; and comparing the second deformation amount in a case where the stress distribution of the main influence region is applied with the third deformation amount, and identifying, as the shape defect factor, the stress component in a direction having the smallest difference as a main influence stress component that is the stress component most related to the shape defect of the workpiece in the evaluation point before and after the elastic recovery.