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

Tutkimustuotos: Lehtiartikkeli

Standard

Harvard

APA

Vancouver

Author

Bibtex - Lataa

@article{fef02cb0c6564c128ddae738fcd321c2,
title = "A free, square, point-loaded ice sheet : A finite element-discrete element approach",
abstract = "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.",
keywords = "Beam lattice network, Centroidal Voronoi tessellation, Finite element-discrete element method, Ice, Plate on a Winkler foundation, Size effect",
author = "Lilja, {Ville Pekka} and Arttu Poloj{\"a}rvi and Jukka Tuhkuri and Jani Paavilainen",
year = "2019",
month = "11",
doi = "10.1016/j.marstruc.2019.102644",
language = "English",
volume = "68",
journal = "Marine Structures; Desing, Construction & Safety",
issn = "0951-8339",
publisher = "Elsevier Ltd.",

}

RIS - Lataa

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