# A free, square, point-loaded ice sheet : A finite element-discrete element approach

Research output: Contribution to journal › Article

### Standard

**A free, square, point-loaded ice sheet : A finite element-discrete element approach.** / Lilja, Ville Pekka; Polojärvi, Arttu; Tuhkuri, Jukka; Paavilainen, Jani.

Research output: Contribution to journal › Article

### Harvard

*Marine Structures*, vol. 68, 102644. https://doi.org/10.1016/j.marstruc.2019.102644

### APA

*Marine Structures*,

*68*, [102644]. https://doi.org/10.1016/j.marstruc.2019.102644

### Vancouver

### Author

### Bibtex - Download

}

### RIS - Download

TY - JOUR

T1 - A free, square, point-loaded ice sheet : A finite element-discrete element approach

AU - Lilja, Ville Pekka

AU - Polojärvi, Arttu

AU - Tuhkuri, Jukka

AU - Paavilainen, Jani

PY - 2019/11

Y1 - 2019/11

N2 - Deflection of a free, thin, square, elastic, point-loaded ice sheet is studied. A three-dimensional, combined finite element-discrete element method (FEM-DEM) is applied. Examined is a set of square, self-similar (plan view) ice sheet samples with the side lengths of L=20, 40, 80, and 160 m; thicknesses of h=0.5, 1.0, and 1.5 m; and the average discrete element sizes of l=2h and 3h. The samples are generated and meshed via a centroidal Voronoi tessellation procedure. Each mesh is unstructured and consists of convex polyhedra as the discrete elements and a Delaunay-triangulated, in-plane beam lattice of co-rotational, viscously damped Timoshenko beam finite elements. Two load cases are considered: a sheet sample subjected to i) a vertical, centrally applied point load; and ii) a vertical point load on an edge. In both load cases, a good agreement is found between the deflections computed with the FEM-DEM approach and with FEM. Maximum, non-dimensionalized deflections w¯max (w¯max=q/Fz, where q=q(wmax) is an approximate, resultant reaction force due to the buoyancy, wmax the computed maximum deflection, and Fz the applied point load) scale, in each load case, either by a power-law, w¯max∝hm1(L), or exponentially, w¯max∝10m2(h)L, depending on whether computed over h (for each L) or vice versa. The exponents m1 and m2 are not constants, but functions of L and h, respectively, and approximately equal in either load case. In-plane size estimates are sought for an FEM-DEM sheet sample to represent either a) an infinite, elastic, point-loaded ice sheet; or b) a semi-infinite, elastic, point-loaded (with the load on the free edge) ice sheet with a free edge. A free, square, point-loaded FEM-DEM sheet sample well approximates both a) and b), if L/2>≈5lch, where lch denotes a characteristic length.

AB - Deflection of a free, thin, square, elastic, point-loaded ice sheet is studied. A three-dimensional, combined finite element-discrete element method (FEM-DEM) is applied. Examined is a set of square, self-similar (plan view) ice sheet samples with the side lengths of L=20, 40, 80, and 160 m; thicknesses of h=0.5, 1.0, and 1.5 m; and the average discrete element sizes of l=2h and 3h. The samples are generated and meshed via a centroidal Voronoi tessellation procedure. Each mesh is unstructured and consists of convex polyhedra as the discrete elements and a Delaunay-triangulated, in-plane beam lattice of co-rotational, viscously damped Timoshenko beam finite elements. Two load cases are considered: a sheet sample subjected to i) a vertical, centrally applied point load; and ii) a vertical point load on an edge. In both load cases, a good agreement is found between the deflections computed with the FEM-DEM approach and with FEM. Maximum, non-dimensionalized deflections w¯max (w¯max=q/Fz, where q=q(wmax) is an approximate, resultant reaction force due to the buoyancy, wmax the computed maximum deflection, and Fz the applied point load) scale, in each load case, either by a power-law, w¯max∝hm1(L), or exponentially, w¯max∝10m2(h)L, depending on whether computed over h (for each L) or vice versa. The exponents m1 and m2 are not constants, but functions of L and h, respectively, and approximately equal in either load case. In-plane size estimates are sought for an FEM-DEM sheet sample to represent either a) an infinite, elastic, point-loaded ice sheet; or b) a semi-infinite, elastic, point-loaded (with the load on the free edge) ice sheet with a free edge. A free, square, point-loaded FEM-DEM sheet sample well approximates both a) and b), if L/2>≈5lch, where lch denotes a characteristic length.

KW - Beam lattice network

KW - Centroidal Voronoi tessellation

KW - Finite element-discrete element method

KW - Ice

KW - Plate on a Winkler foundation

KW - Size effect

UR - http://www.scopus.com/inward/record.url?scp=85070686579&partnerID=8YFLogxK

U2 - 10.1016/j.marstruc.2019.102644

DO - 10.1016/j.marstruc.2019.102644

M3 - Article

VL - 68

JO - Marine Structures; Desing, Construction & Safety

JF - Marine Structures; Desing, Construction & Safety

SN - 0951-8339

M1 - 102644

ER -

ID: 36366001