This paper presents a finite element model for the simulation of ice–structure interaction problems, which are dominated by crushing. The failure mode of ice depends significantly on the strain rate. At low strain rates, the ice behaves ductile, whereas at high strain rates it reacts in brittle mode. This paper focuses on the brittle mode, which is the dominating mode for ship–ice interactions. A multitude of numerical approaches for the simulation of ice can be found in the literature. Nevertheless, the literature approaches do not seem suitable for the simulation of continuous ice–structure interaction processes at low and high confinement ratios in brittle mode. Therefore, this paper seeks to simulate the ice–structure interaction with the finite element method (FEM). The objective of the here introduced Mohr-Coulomb Nodal Split (MCNS) model is to represent the essential material behavior of ice in an efficient formulation. To preserve mass and energy as much as possible, the node splitting technique is applied, instead of the frequently used element erosion technique. The intention of the presented model is not to reproduce individual cracks with high accuracy, because this is not possible with a reasonable element size, due to the large number of crack fronts forming during the ice–structure interaction process. To validate the findings of the model, the simulated maximum ice forces and contact pressures are compared with ice extrusion and double pendulum tests. During validation, the MCNS model shows a very good agreement with these experimental values.