. In this case, out-of-plane deformations
are much smaller than in-plane deformation. To avoid
dealing with highly deformed mesh, an Eulerian approach was
used for the in-plane direction. The dotted region in Fig. 1b corresponds
to the model domain that is meshed. The bottom half of the
pre-consolidated sample covers half of the domain and the rest is
initially filled with air.