【摘要】 极地海冰在多种空间尺度上呈现出离散分布特征，本文基于Gilbert-Johnson-Keerthi(GJK)和Expanding Polytope Algorithm(EPA)接触算法发展一种用于极地海冰的离散元方法。该方法允许模拟浮冰的漂移和接触碰撞过程，并基于简化的刚体有限元方法表征海冰的冻结破坏。通过在平整冰场上进行单轴压缩模拟证实程序的稳定性，并研究压缩强度对浮冰场变形的影响。在此基础上，对穿越狭窄通道的冰块集合进行离散元模拟，研究浮冰尺度、初始密集度和圆形度对其流动特性的影响规律。结果表明，在幂律分布的浮冰尺度分布函数中增大幂指数、增加圆形度以及降低初始密集度都会导致浮冰场的流动速度加快。更多还原

【Abstract】 In polar regions, floating ice exhibits distinct characteristics across a range of spatial scales. It is well recognized that the irregular geometry of these ice formations markedly influences their dynamic behavior. This study introduces a polyhedral Discrete Element Method(DEM) tailored for polar ice, incorporating the Gilbert-Johnson-Keerthi(GJK) and Expanding Polytope Algorithm(EPA) for contact detection. This approach facilitates the simulation of the drift and collision processes of floating ice, effectively capturing its freezing and fragmentation. Subsequently, the stability and reliability of this model are validated by uniaxial compression on level ice fields, focusing specifically on the influence of compression strength on deformation resistance. Additionally, clusters of ice floes navigating through narrow channels are simulated. These studies have qualitatively assessed the effects of Floe Size Distribution(FSD), initial concentration, and circularity on their flow dynamics. The higher power-law exponent values in the FSD, increased circularity, and decreased concentration are each associated with accelerated flow in ice floe fields. The simulation results distinctly demonstrate the considerable impact of sea ice geometry on the movement of clusters, offering valuable insights into the complexities of polar ice dynamics.更多还原