Finite-discrete element modelling of sea ice sheet fracture

INTERNATIONAL JOURNAL OF SOLIDS AND STRUCTURES

卷:217

Article

Lilja, Ville-Pekka¹  Polojarvi, Arttu¹  Tuhkuri, Jukka¹  Paavilainen, Jani² 

[1] Aalto Univ, Sch Engn, Dept Mech Engn, POB 14300, FI-00076 Aalto, Finland

[2] Rand Simulat Ltd, Vantaa, Finland

关键词: Dynamic fracture;    Plates;    Numerical algorithms;    Size effect;    Ice and snow;   

DOI  :  10.1016/j.ijsolstr.2020.11.028

来源: Elsevier

【 摘 要 】

A rate-independent, de-cohesive damage model for the fracture modelling of large, cellular, plate-like, quasi-brittle structures is proposed. A hybrid, three-dimensional finite-discrete element method to investigate sea ice sheet fracture is then introduced, followed by three applications. The uniaxial tensile fracture of an ice sheet of varying physical sizes is examined first. The effects of both the size of an ice sheet and the loading rate applied on the effective tensile strength are investigated. The vertical penetration fracture of an ice sheet loaded by a rigid, flat-ended, cylindrical indenter is examined next. The breakthrough loads and strengths of an ice sheet of varying physical sizes are computed, applicable scaling rules as regards to the vertical breakthrough strength searched for. To conclude, the breaking of an ice sheet containing a circular hole by a surfacing, rigid, truncated cone is studied (an axisymmetric contact problem). The loads on the cone are computed and then compared with loads that can be obtained analytically for a case in which a structure is stationary, a sheet moves, and the contact is unilateral. While computing the tensile and the breakthrough strengths, a set of self-similar sheet samples with an inplane size range of 1:16 is examined. The samples are square; have a side length of either L = 10, 20, 40, 80, or 160 m; and a thickness of either h = 0.5, 1.0, or 1.5 m. With the sheets containing holes, only the largest samples (L = 160 m) are investigated. The results indicate that 1) both the tensile and the breakthrough strengths are strong functions of both L and h; ii) the tensile strength is a strong function of the applied loading rate; iii) the failure mode as regards to the vertical penetration fracture changes drastically as a function of L; iv) the model is able to demonstrate both radial and circumferential cracking; and that v) the proposed (in-direct) approach to compute ice loads on a conical offshore structure provides realistic results. (C) 2020 Elsevier Ltd. All rights reserved.

【 授权许可】

Free