Numerical simulations of parasitic folding in multilayers