• International Journal of Naval Architecture and Ocean Engineering
• /
• v.11 no.2
• /
• pp.639-647
• /
• 2019

A realistic numerical simulation technology using a Lagrangian Fluid-Structure Interaction (FSI) model was combined with a fracture algorithm to predict the fluid-ice-structure interaction. The failure of ice was modeled as the tensile fracture of elastic material by applying a novel FSI model based on the Moving Particle Semi-implicit (MPS) method. To verify the developed fracture algorithm, a series of numerical simulations for 3-point bending tests with an ice beam were performed and compared with the experiments carried out in an ice room. For application of the developed FSI model, a dropping water droplet hitting a cantilever ice beam was simulated with and without the fracture algorithm. The simulation showed that the effects of fracture which can occur in the process of a FSI simulation can be studied.

• Structural Engineering and Mechanics
• /
• v.36 no.4
• /
• pp.499-511
• /
• 2010

This paper examines the ice cover effects on the seismic response of concrete gravity dam-reservoir-foundation interaction systems subjected to a horizontal earthquake ground motion. ANSYS program is used for finite element modeling and analyzing the ice-dam-reservoir-foundation interaction system. The ice-dam-reservoir interaction system is considered by using the Lagrangian (displacementbased) fluid and solid-quadrilateral-isoparametric finite elements. The Sariyar concrete gravity dam in Turkey is selected as a numerical application. The east-west component of Erzincan earthquake, which occurred on 13 March 1992 in Erzincan, Turkey, is selected for the earthquake analysis of the dam. Dynamic analyses of the dam-reservoir-foundation interaction system are performed with and without ice cover separately. Parametric studies are done to show the effects of the variation of the length, thickness, elasticity modulus and density of the ice-cover on the seismic response of the dam. It is observed that the variations of the length, thickness, and elasticity modulus of the ice-cover influence the displacements and stresses of the coupled system considerably. Also, the variation of the density of the ice-cover cannot produce important effects on the seismic response of the dam.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.6 no.4
• /
• pp.876-893
• /
• 2014

The resistance performance of an icebreaking cargo vessel according to the variation of waterline angles is investigated numerically and experimentally. A recently developed Finite Element (FE) model is used in our analysis. A resistance test with synthetic ice is performed in the towing tank at Pusan National University (PNU) to compare and validate the computed results. We demonstrate good agreement between the experimental and numerical results. Shipice interaction loads are numerically calculated based on the Fluid Structure Interaction (FSI) method in the commercial FE package LS-DYNA. Test results from model testing with synthetic ice at the PNU towing tank are used to compare and validate the numerical simulations. For each waterline angle, numerical and experimental comparisons were made for three concentrations (90%, 80%, and 60%) of pack ice. Ice was modeled as a rigid body, but the ice density was the same as that used in the experiments. A comparative study according to the variation of stem angles is expected to be conducted in the near future.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.5 no.1
• /
• pp.116-131
• /
• 2013

The resistance performance of an icebreaking cargo vessel in pack ice conditions was investigated numerically and experimentally using a recently developed finite element (FE) model and model tests. A comparison between numerical analysis and experimental results with synthetic ice in a standard towing tank was carried out. The comparison extended to results with refrigerated ice to examine the feasibility of using synthetic ice. Two experiments using two different ice materials gave a reasonable agreement. Ship-ice interaction loads are numerically calculated based on the fluid structure interaction (FSI) method using the commercial FE package LS-DYNA. Test results from model testing with synthetic ice at the Pusan National University towing tank, and with refrigerated ice at the National Research Council's (NRC) ice tank, are used to validate and benchmark the numerical simulations. The designed ice-going cargo vessel is used as a target ship for three concentrations (90%, 80%, and 60%) of pack ice conditions. Ice was modeled as a rigid body but the ice density was the same as that in the experiments. The numerical challenge is to evaluate hydrodynamic loads on the ship's hull; this is difficult because LS-DYNA is an explicit FE solver and the FSI value is calculated using a penalty method. Comparisons between numerical and experimental results are shown, and our main conclusions are given.

• Journal of Advanced Research in Ocean Engineering
• /
• v.2 no.2
• /
• pp.67-80
• /
• 2016

The resistance performance of an icebreaking cargo vessel with varied stem angles is investigated numerically and experimentally. Ship-ice interaction loads are numerically calculated based on the fluid structure interaction (FSI) method using the commercial FE package LS-DYNA. Test results obtained from model testing with synthetic ice at the Pusan National University towing tank and with refrigerated ice at the National Research Council's (NRC) ice tank are used to validate and benchmark the numerical simulations. The designed icebreaking cargo vessel with three stem angles ($20^{\circ}$, $25^{\circ}$, and $30^{\circ}$) is used as the target ship for three concentrations (90%, 80%, and 60%) of pack ice conditions. The comparisons between numerical and experimental results are shown and our main conclusions are given.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.12 no.1
• /
• pp.414-427
• /
• 2020

The finite element method is used to simulate the navigation of an ice-area bulk carrier in broken ice fields. The ice material is defined as elastic, and the simulations are accomplished at four model speeds and three ice concentrations. The movements of ice floes in the simulation are consistent with those in the model test, and the percentage deviation of the numerical ice resistance from the ice resistance in the model test can be controlled to be less than 15 %. The key characteristics of ice loads, including the average ice loads, extreme ice loads, and characteristic frequency, are analyzed thoroughly in a comprehensive manner. Moreover, the effects of sailing speed and ice concentration on the ice loads are analyzed. In particular, the stress distribution of ice floes is presented to help understand how model speed and concentration affect the ice loads. The "ice pressure" phenomenon is observed at 90 % ice concentration, and it is realistically reflected both in the time―and frequency―domain ice force curves.

• Structural Engineering and Mechanics
• /
• v.75 no.4
• /
• pp.487-496
• /
• 2020

Most of the previous works on numerical analysis of galloping of transmission lines are generally based on the quasisteady theory. However, some wind tunnel tests of the rectangular section or hangers of suspension bridges have shown that the galloping phenomenon has a strong unsteady characteristic and the test results are quite different from the quasi-steady calculation results. Therefore, it is necessary to check the applicability of the quasi-static theory in galloping analysis of the ice-covered transmission line. Although some limited unsteady simulation researches have been conducted on the variation of parameters such as aerodynamic damping, aerodynamic coefficients with wind speed or wind attack angle, there is a need to investigate the numerical simulation of unsteady galloping of two-dimensional iced transmission line with comparison to wind tunnel test results. In this paper, it is proposed to conduct a two dimensional (2-D) unsteady numerical analysis of ice-covered transmission line galloping. First, wind tunnel tests of a typical crescent-shapes iced conductor are conducted firstly to check the subsequent quasisteady and unsteady numerical analysis results. Then, a numerical simulation model consistent with the aeroelastic model in the wind tunnel test is established. The weak coupling methodology is used to consider the fluid-structure interaction in investigating a two-dimension numerical simulation of unsteady galloping of the iced conductor. First, the flow field is simulated to obtain the pressure and velocity distribution of the flow field. The fluid action on the iced conduct at the coupling interface is treated as an external load to the conductor. Then, the movement of the conduct is analyzed separately. The software ANSYS FLUENT is employed and redeveloped to numerically analyze the model responses based on fluid-structure interaction theory. The numerical simulation results of unsteady galloping of the iced conduct are compared with the measured responses of wind tunnel tests and the numerical results by the conventional quasi-steady theory, respectively.