We develop a numerical method to simulate the fracture in kilometer-sized chunks of floating ice in the ocean. Our approach uses a mathematical model that balances deformation energy against the energy required for fracture. We study the strength of ice chunks that contain random impurities due to prior damage or refreezing and what types of fractures are likely to occur. Our model shows that crack direction critically depends on the orientation of impurities relative to surrounding forces.