This article presents two methods to determine pressing forces to the M1 copper plate. The first approach employs the upper-bound method, and the second method combines numerical simulation with a central composite design. The upper-bound method uses physical analysis to determine the deformation zone and establish the division of the rigid block model. The findings were identical to the slip-line solution method. The results obtained from the three-dimensional simulation performed using Qform software combined with a central composite design have provided regression equations that help determine the force and punching pressure. The impact of each parameter on the pressing force is similar in both approaches, and the difference between the two methods is insignificant and falls within a dependable range.