Since this is a problem with an indefinite depth, my suggestion is to simulate it using a 2D model with plane strain elements for the sheet, and rigid bodies for the punch and the die. Then, you can compare the punch's force vs displacement plot with the analytical solution. Of course, you have to rearrange the equation a little to transform theta into a vertical displacement of the punch.