本文分别基于有限单元法（Finite Element Method，FEM）和有限离散单元法（Finite Discrete Element Method，FDEM）对冰水池中的典型柱状模型冰开展了单轴压缩数值模拟，并对2种数值方法进行对比研究。FEM数值模拟时利用Drucker Prager塑性模型和延性损伤模型表征了模型冰的非弹性阶段力学行为和损伤。FDEM数值模拟时则利用Cohesive单元模拟模型冰的失效破碎行为。此外，在FDEM数值模型中考虑模型冰晶界和晶间强度差异，仅模拟晶界处的失效。研究发现，基于FDEM的模型冰裂纹扩展形态和物理试验更为相似。最后通过网格映射的方法，改善了数值模型网格的质量，提高了FDEM的计算效率。本文工作可为模型冰力学特性精细模拟研究提供有益参考。
The present study focuses on the uniaxial compression numerical simulation of columnar model ice in ice tank. The simulation is conducted using the Finite Element Method (FEM) and the Finite Discrete Element Method (FDEM), and a comparison is made between the two numerical methods. In the FEM simulation, the Drucker Prager plastic model and ductile damage model are employed to characterize the mechanical behavior and inelastic stage and damage of the model ice. On the other hand, the FDEM simulation utilizes Cohesive element to simulate the failure and fragmentation behavior of model ice. The FDEM numerical model also considers the difference between the ice grain boundary and intergranular strength, focusing solely on the failure at the grain boundary. The simulation results indicate that the crack growth morphology based on the FDEM and physical tests exhibit greater similarity. Additionally, the quality of the numerical model grid is enhanced through the application of the grid mapping method, leading to improved computational efficiency of the FDEM. This research provides valuable insights for accurately simulating the mechanical properties of model ice, and can serve as a useful reference for future endeavors in this field.
2025,47(19): 1-6 收稿日期：2024-11-13
DOI：10.3404/j.issn.1672-7649.2025.19.001
分类号：U661.4
基金项目：工信部高技术船舶科研计划项目（[2021]342）；国家自然科学基金重大项目（52101331，52192690，52192694）；基础科研项目（JCKY2020206B073）
作者简介：余朝歌(1997-),男,博士研究生,研究方向为冰区结构冰载荷
参考文献：
[1] 田于逵, 岳前进, 孔帅, 等. 冰水池中海冰模拟技术现状与关键问题[J]. 船舶力学, 2023, 27(1): 153-162.
[2] KARL-ULRICH EVERS. Model tests with ships and offshore structures in HSVA's ice tanks[C]//Proceedings of the 24th International Conference on Port and Ocean Engineering under Arctic Conditions, Busan, Korea: ISOPE, 2017.
[3] GESA ZIEMER. Hsva model ice – a status report[C]//Proceedings of the ASME 2022 41st International Conference on Ocean, Offshore and Arctic Engineering , Hamburg, Germany: OMAE, 2022.
[4] FRANZ VON BOCK UND POLACH. The mechanical behavior of model-scale ice: experiments, numerical modeling and scalability[D]. Aalto University, 2016.
[5] TIAN Y K. Design and realization of CSSRC small ice model basin for ice related fundamental researches[C]//Proceedings of the 25th International Conference on Port and Ocean Engineering under Arctic Conditions, Delft, The Netherlands: POAC, 2019.
[6] ITTC. Recommended procedures and guidelines: Test methods for model ice properties[S]. 2019.
[7] TIAN Y K. Characterization of uniaxial compression strength for columnar saline model ice in CSSRC small ice model basin[J]. Journal of Ship Mechanics, 2020, 24(12): 1647-1656
[8] JI S Y, et al. DEM simulation of uniaxial compressive and flexural strength of sea ice: parametric study[J]. Journal of Engineering Mechanics, 2017, 143(1): 10-11.
[9] DI S, XUE Y, BAI X, et al. Effects of model size and particle size on the response of sea-ice samples created with a hexagonal-close-packing pattern in discrete-element method simulations[J]. Particuology, 2018, 36: 106-113.
[10] BATEMAN S P, ORZECH M D, CALANTONI J. Simulating the mechanics of sea ice using the discrete element method[J]. Mechanics Research Communications, 2019, 99: 73-78.
[11] 周陈超. 不同破坏准则下极地海冰有限元材料和粘聚单元法的应用[D]. 镇江: 江苏科技大学, 2022.
[12] LU W J. Cohesive zone method based simulations of ice wedge bending: a comparative study of element erosion, CEM, DEM and XFEM[C]//21st IAHR International Symposium on Ice, 2012 Dalian University of Technology Press, Dalian, 2012.
[13] Abaqus Analysis User's Guide, DS Simulia Abaqus, 2016.
[14] ABDELAZIZ A, ZHAO Q, GRASSELLI G. Grain based modelling of rocks using the combined finite-discrete element method[J]. Computers and Geotechnics, 2018, 10(3): 73-81.
[15] TIAN Y K. Experimental investigation on ice loading mechanism of typical offshore structures under the action of level ice[C]//Proceedings of the Thirty-second (2022) International Ocean and Polar Engineering Conference, Shanghai, China: ISOPE, 2022.
