A numerical model for phase-resolved simulation of nonlinear ocean waves propagating through fragmented sea ice is proposed. This model solves the full time-dependent equations for nonlinear potential flow coupled with a nonlinear thin-plate formulation for the ice cover. The coefficient of flexural rigidity is allowed to vary spatially so that distributions of ice floes can be directly specified in the physical domain. Two-dimensional simulations are performed to examine the attenuation of solitary waves by scattering through an irregular array of ice floes. This is joint work with Emilian Parau (UEA).