<p>To investigate the ice load of level ice interacting with a narrow conical structure, a bonded circular elements model is introduced based on discrete element method. In the simulation, the level ice is modelled with bonded circular elements, which can transmit forces and moments to cause the bonded circular element to move as a whole. The bond breaking process is considered the fracture process in the level ice when it interacts with the conical structure. The strength of the bond is determined by the sea ice physical tests. The theoretical foundation of the simulator is explained in detail in this paper. As the main failure mode of the level ice interaction with a conical structure is bending failure, the parameterization of the bonded model is validated with a three-point bending test in the laboratory. To render the numerical fracture size of the level ice consistent with the measured data, single layer and uniform circular elements are used in the simulation. Finally, the ice load of the level ice interacting with the narrow conical structure is simulated, and the characteristics of numerical ice force are compared with the field measured data. It is shown that the numerical ice forces can demonstrate similar behavior observed in the field, including the loading stage, unloading stage and no-load stage.</p>