Delamination between layers occurs during the creasing and subsequent folding of paperboard. Delamination is necessary to provide some stiffness properties, but excessive or uncontrolled delamination can weaken the fold, and therefore needs to be controlled. An understanding of the mechanics of delamination is predicated upon the availability of reliable and properly calibrated simulation tools to predict experimental observations. This paper describes a finite element simulation of paper mechanics applied to the scoring and folding of multi-ply carton board. Our goal was to provide an understanding of the mechanics of these operations and the proper models of elastic and plastic behavior of the material that enable us to simulate the deformation and delamination behavior. Our material model accounted for plasticity and sheet anisotropy in the in-plane and z-direction (ZD) dimensions. We used different ZD stress-strain curves during loading and unloading. Material parameters for in-plane deformation were obtained by fitting uniaxial stress-strain data to Ramberg-Osgood plasticity models and the ZD deformation was modeled using a modified power law. Two-dimensional strain fields resulting from loading board typical of a scoring operation were calculated. The strain field was symmetric in the initial stages, but increasing deformation led to asymmetry and heterogeneity. These regions were precursors to delamination and failure. Delamination of the layers occurred in regions of significant shear strain and resulted primarily from the development of large plastic strains. The model predictions were confirmed by experimental observation of the local strain fields using visual microscopy and linear image strain analysis. The finite element model predicted sheet delamination matching the patterns and effects that were observed in experiments.