We use a centimeter-scale 2D numerical model to investigate the effect of the presence of strong inclusions on strain localization in a visco-elastic ice matrix. The presence of hard particles strengthens the aggregate, but also causes strain localization and the formation of ductile shear zones. Comparing our model to the laboratory uniaxial compression experiments on ice cylinders allows the analysis of transient and steady-state strain distribution in ice matrix, and strain partitioning between ice and second phases through empirical calibration of viscous softening parameters.


The video shows a numerical experiment where a roughly square-centimetre-sized sample of polymineralic material is compressed. The light blue rhombs are strong calcite inclusions, which are embedded in water ice. The reddish areas depict regions of high strain that localise in discrete shear bands within the soft ice matrix.

Cite as:

Cyprych, D., Brune, S., Piazolo, S., Quinteros, J.: Strain localization in polycrystalline material with second phase particles: Numerical modeling with application to ice mixtures: Geochemistry, Geophysics, Geosystems, doi:10.1002/2016GC006471.

Downloads and Resources:

Link to Paper | Link to Supplementary Information | Model animations (mp4): StrainStrain RateViscosity | Model animations (youtube): Strain, Strain Rate, Viscosity