• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1999.04a
• /
• pp.321-326
• /
• 1999

In this study, fatigue crack propagation problem of plate under multi-axial loading is mainly considered To analyze this special problem, recently developed technique called EFGM(Element-Free Galerkin Method), one of the Meshfree Methods, and general fatigue crack growth raw herein Paris law are used Using the Implemented scheme, paths of fatigue cracks by constant-amplitude load fluctuation and multiple-crack growth behavior are examined. The failure mechanism of steel plate due to crack propagation is studied. As a result, an algorithm that treats multiple fatigue crack problems is proposed. A numerical example shows that the prediction of growing paths can be achieved successfully and efficiently by proposed algorithm.

• Interaction and multiscale mechanics
• /
• v.2 no.4
• /
• pp.333-352
• /
• 2009

We introduce a radial basis collocation method to solve axially moving beam problems which involve $2^{nd}$ order differentiation in time and $4^{th}$ order differentiation in space. The discrete equation is constructed based on the strong form of the governing equation. The employment of multiquadrics radial basis function allows approximation of higher order derivatives in the strong form. Unlike the other approximation functions used in the meshfree methods, such as the moving least-squares approximation, $4^{th}$ order derivative of multiquadrics radial basis function is straightforward. We also show that the standard weighted boundary collocation approach for imposition of boundary conditions in static problems yields significant errors in the transient problems. This inaccuracy in dynamic problems can be corrected by a statically condensed semi-discrete equation resulting from an exact imposition of boundary conditions. The effectiveness of this approach is examined in the numerical examples.

• Transactions of the Korean Society of Machine Tool Engineers
• /
• v.12 no.5
• /
• pp.73-79
• /
• 2003

In this paper, RKPM is extended for solving moderately thick and thin beams. General Timoshenko beam theory is used for formulation. Shear locking is the main difficulty in analysis of these kinds of structures. Shear relaxation factor, which is formulated using the difference between bending and shear strain energy, and corrected shear rigidity are introduced to overcome shear locking. Analysis results obtained reveal that RKPM using introduced methods is free of locking and very effectively applicable to deep beams as well as shallow beams.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.9 no.2
• /
• pp.209-218
• /
• 2017

Sea ice is the main factor affecting the safety of the Arctic engineering. However, traditional numerical methods derived from classical continuum mechanics have difficulties in resolving discontinuous problems like ice damage. In this paper, a non-local, meshfree numerical method called "peridynamics", which is based on integral form, was applied to simulate the interaction between level ice and a cylindrical, vertical, rigid structure at different velocities. Ice in the simulation was freshwater ice and simplified as elastic-brittle material with a linear elastic constitutive model and critical equivalent strain criterion for material failure in state-based peridynamics. The ice forces obtained from peridynamic simulation are in the same order as experimental data. Numerical visualization shows advantages of applying peridynamics on ice damage. To study the repetitive nature of ice force, damage zone lengths of crushing failure were computed and conclude that damage zone lengths are 0.15-0.2 times as ice thickness.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.27 no.1
• /
• pp.1-8
• /
• 2014

Using the level set and the meshfree methods, we develop a topological shape optimization method applied to linear elasticity problems. Design gradients are computed using an efficient adjoint design sensitivity analysis(DSA) method. The boundaries are represented by an implicit moving boundary(IMB) embedded in the level set function obtainable from the "Hamilton-Jacobi type" equation with the "Up-wind scheme". Then, using the implicit function, explicit boundaries are generated to obtain the response and sensitivity of the structures. Global nodal shape function derived on a basis of the reproducing kernel(RK) method is employed to discretize the displacement field in the governing continuum equation. Thus, the material points can be located everywhere in the continuum domain, which enables to generate the explicit boundaries and leads to a precise design result. The developed method defines a Lagrangian functional for the constrained optimization. It minimizes the compliance, satisfying the constraint of allowable volume through the variations of boundary. During the optimization, the velocity to integrate the Hamilton-Jacobi equation is obtained from the optimality condition for the Lagrangian functional. Compared with the conventional shape optimization method, the developed one can easily represent the topological shape variations.

• Structural Engineering and Mechanics
• /
• v.65 no.5
• /
• pp.557-572
• /
• 2018

The academic research works about liquid storage tanks are reviewed for the purpose of providing valuable reference to the engineering practice on their aseismic design. A summary of the performance of tanks during past earthquakes is described in this paper. Next, the seismic response of tanks under unidirectional earthquake is reported, supplemented with the dynamic response under multidirectional motions. Then, researches on the influence of soil-structure interaction are brought out to help modify the seismic design approach of tanks in different areas with variable properties of soils. Afterwards, base isolation systems are reported to demonstrate their effectiveness for the earthquake-resistant design of liquid storage tanks. Further, researches about the liquid-structure interaction are reviewed with description of simplified models and numerical analytical methods, some of which consider the elastic effect of tank walls. Moreover, the liquid sloshing phenomenon on the hydrodynamic behaviors of tanks is presented by various algorithms including grid-based and meshfree method. And then the impact of baffles in changing the dynamic characteristics of the liquid-structure system is raised, which shows the energy dissipation by the vortex motion of liquid. In addition, uplifting effect is given to enhance the understanding on the capacity of unanchored tanks and some assessment of their development. At last, the concluding remarks and the aspects of extended research in the field of liquid storage tanks under seismic loads are provided, emphasizing the thermal stress analysis, the replaceable system for base isolation, the liquid-solid interaction and dynamic responses with stochastic excitations.

• Interaction and multiscale mechanics
• /
• v.6 no.2
• /
• pp.211-235
• /
• 2013

An incompressible smoothed particle hydrodynamics (ISPH) method based on the incremental pressure projection method is developed in this study. The Rayleigh-B$\acute{e}$nard convection in a square enclosure is used as a validation case and the results obtained by the proposed ISPH model are compared to the benchmark solutions. The comparison shows that the established ISPH method has a good performance in terms of accuracy. Subsequently, the proposed ISPH method is employed to simulate natural convection from a heated cylinder in a square enclosure. It shows that the predictions obtained by the ISPH method are in good agreements with the results obtained by previous studies using alternative numerical methods. A rotating and heated cylinder is also considered to study the effect of the rotation on the heat transfer process in the enclosure space. The numerical results show that for a square enclosure at, the addition of kinetic energy in the form of rotation does not enhance the heat transfer process. The method is also applied to simulate forced convection from a circular cylinder in an unbounded uniform flow. In terms of results, it turns out that the proposed ISPH model is capable to simulate heat transfer problems with the complex and moving boundaries.

• Advances in environmental research
• /
• v.3 no.1
• /
• pp.45-69
• /
• 2014

The groundwater has been a major source of water supply throughout the ages. Around 50% of the rural as well as urban population in the developing countries like India depends on groundwater for drinking. The groundwater is also an important source in the agriculture and industrial sector. In many parts of the world, groundwater resources are under increasing threat from growing demands, wasteful use and contamination. A good planning and management practices are needed to face this challenge. A key to the management of groundwater is the ability to model the movement of fluids and contaminants in the subsurface environment. It is obvious that the contaminant source activities cannot be completely eliminated and perhaps our water bodies will continue to serve as receptors of vast quantities of waste. In such a scenario, the goal of water quality protection efforts must necessarily be the control and management of these sources to ensure that released pollutants will be sufficiently attenuated within the region of interest and the quality of water at points of withdrawal is not impaired. In order to understand the behaviour of contaminant transport through different types of media, several researchers are carrying out experimental investigations through laboratory and field studies. Many of them are working on the analytical and numerical studies to simulate the movement of contaminants in soil and groundwater of the contaminant transport. With the advent of high power computers especially, a numerical modelling has gained popularity and is indeed of particular relevance in this regard. This paper provides the state of the art of contaminant transport and reviews the allied research works carried out through experimental investigation or using the analytical solution and numerical method. The review involves the investigation in respect of both, saturated and unsaturated, porous media.

• The Sea:JOURNAL OF THE KOREAN SOCIETY OF OCEANOGRAPHY
• /
• v.23 no.4
• /
• pp.204-217
• /
• 2018

This study proposes a cost-effective method to monitor coastal environments using e-Navigation-implemented domestic and international ferries, and to analyze big data of records such as wind, temperature, salinity, waves, and currents that are gathered through e-Navigation system. First, we present the concept and architecture of e-Navigation operation system based on the General Information Center on Maritime Safety and Security. Then, the marine observation system that can be applied to ferries operating in our nation's territory is discussed. Analytical methods, such as spatio-temporal mixed effects model, ensemble method, and meshfree method, in handling real-time big data obtained by the e-Navigation observing system are then explained in detail. This study will support the implementation of the Korean e-Navigation project that focuses on the safety of small vessels such as coasters and fishing vessels.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.27 no.5
• /
• pp.337-343
• /
• 2014

This paper presents a continuum-based shape design sensitivity analysis(DSA) method for crack propagation problems using a reproducing kernel method(RKM), which facilitates the remeshing problem required for finite element analysis(FEA) and provides the higher order shape functions by increasing the continuity of the kernel functions. A linear elasticity is considered to obtain the required stress field around the crack tip for the evaluation of J-integral. The sensitivity of displacement field and stress intensity factor(SIF) with respect to shape design variables are derived using a material derivative approach. For efficient computation of design sensitivity, an adjoint variable method is employed tather than the direct differentiation method. Through numerical examples, The mesh-free and the DSA methods show excellent agreement with finite difference results. The DSA results are further extended to a shape optimization of crack propagation problems to control the propagation path.