For the problem of forming defects in the rounded corner area at the top of automobile energy absorbing box that was prone to excessive material thinning and cracking during the stamping, taking pressing force, stripping force and die clearance as the test factors and the minimization of  maximum thinning rate of energy-absorbing box and the maximization of safe domain proportion of forming limit diagram (FLD) as the quality optimization objective, a variety of approximate models between experimental factors and quality optimization objective were constructed by the Latin Hypercube experiment design method combined with finite element analysis, and the prediction accuracy of the approximate models was analyzed. Furthermore, based on the multi-objective particle swarm algorithm (MOPSO), the multi-objective optimization calculation was carried out within the optimized Kriging approximate model to obtain the Pareto solution set, and based on the entropy weight approximation ideal solution sorting method (TOPSIS), a set of optimal process parameter, combination was determined from Pareto solution set. Finally, the process was simulated and verified by actual stamping production. The experimental results show that the proposed method is effective and can provide useful reference for the stamping production of automobile energy-absorbing box with similar structure.