Many prototypes and experiments are repeated in order to decide the specification of the optimum corrugated fiberboard box, so if it is possible to predict the compressive strength of the box in advance, the development period can be shortened. Therefore, a method of predicting compressive strength of corrugated fiberboard boxes using nonlinear finite element analysis was studied. Considering geometric nonlinearity and material nonlinearity, phenomena leading to the maximum of the load through bifurcation buckling, initial yielding and material failure were reproduced. As the material properties, the order of the buckling mode and the box compressive strength matched the experimental results by inputting the data of the stressstrain curve obtained by the edgewise crush test of the corrugated fiberboard and the rotational rigidity around the crease line axis. As factors influencing the box compressive strength, it is conceivable that the edgewise crush resistance of the corrugated fiberboard, the bending rigidity of the box corner part, the buckling mode, and the shape of the box. It was proved that estimation of the box compressive strength by numerical simulation is effective as a means to realize optimum design of the box specification and shortening the turnaround time of new product development.