• Journal of Ocean Engineering and Technology
• /
• v.27 no.4
• /
• pp.90-97
• /
• 2013

Oil booms are one of the most widely used types of equipment for the protection of coastal areas against oil spills. In some situations, however, there are several types of oil leaks from the oil boom. Important factors regarding these phenomena include the surrounding ocean environment, such as waves, the density and viscosity of oil, the length of the oil boom skirt, etc. To estimate the performance of the oil boom, it is necessary to predict the behavior of the spilled oil and oil boom. In the present study, the prediction of oil boom performance in waves was carried out using the Pusan-National-University-modified Moving Particle Semi-implicit (PNU-MPS) method, which is an improved version of the original MPS proposed by Koshizuka and Oka (1996). The governing equations, which consist of continuity and Navier-Stokes equations, are solved by Lagrangian moving particles, and all terms expressed by differential operators in the governing equations are replaced by the particle interaction models based on a kernel function. The simulation results were validated through a comparison with the results of Violeau et al. (2007)..

• Transactions of Materials Processing
• /
• v.7 no.3
• /
• pp.233-245
• /
• 1998

A new approach has been proposed for the incremental analysis of the nonsteady state large deformation of planar anisotropic elastic-plastic sheet forming. A mathematical brief review of a constitutive law for the incremental deformation theory has been presented from flow theory using the minimum plastic work path for elastic-plastic material. Since the material embedded coordinate system(Lagrangian quantity) is used in the proposed theory the stress integration procedure is completely objective. A new return mapping algorithm has been also developed from the general midpoint rule so as to achieve numerically large strain increment by successive control of yield function residuals. Some numerical tests for the return mapping algorithm were performed using Barlat's six component anisotropic stress potential. Performance of the proposed algorithm was shown to be good and stable for a large strain increment, For planar anisotropic sheet forming updating algorithm of planar anisotropic axes has been newly proposed. In order to show the effectiveness and validity of the present formulation earing simulation for a cylindrical cup drawing and front fender stamping analysis are performed. From the results it has been shown that the present formulation can provide a good basis for analysis for analysis of elastic-plastic sheet metal forming processes.

• Journal of computational fluids engineering
• /
• v.4 no.1
• /
• pp.19-26
• /
• 1999

This paper describes a numerical method for predicting the incompressible unsteady laminar three-dimensional flows with free-surface. The Navier-Stokes equations governing the flows have been discretized by means of finite-difference approximations, and the resulting equations have been solved via the SIMPLE-C algorithm. The free-surface is defined by the motion of a set of marker particles and the interface behaviour was investigated by means of a "Lagrangian" technique. Using the GALA concept of Spalding, the conventional mass continuity equation is modified to form a volumetric or bulk-continuity equation. The use of this bulk-continuity relation allows the hydrodynamic variables to be computed over the entire flow domain including both liquid and gas regions. Thus, the free-surface boundary conditions are imposed implicitly and the problem formulation is greatly simplified. The numerical procedure is validated by comparing the predicted results of a periodic standing waves problems with analytic solutions. The results show that this numerical method produces accurate and physically realistic predictions of three-dimensional free-surface flows.

• Journal of Ocean Engineering and Technology
• /
• v.31 no.6
• /
• pp.388-396
• /
• 2017

A Lagrangian approach based computational fluid dynamics (CFD) was used to simulate large and/or sharp deformations and fragmentations of interfaces, including free surfaces, through tracing each particle with physical quantities. According to the concept of the particle-based CFD method, it is possible to apply it to both fluid particles and solid particles such as sand, gravel, and rock. However, the presence of more than two different phases in the same domain can make it complicated to calculate the interaction between different phases. In order to solve multiphase problems, particle interaction models for multiphase problems, including surface tension, buoyancy-correction, and interface boundary condition models, were newly adopted into the moving particle semi-implicit (MPS) method. The newly developed MPS method was used to simulate a typical validation problem involving dam breaking. Because the soil and other particles, excluding the water, may have different viscosities, various viscosity coefficients were applied in the simulations for validation. The newly developed and validated MPS method was used to simulate the mobile beds induced by broken dam flows. The effects of the viscosity on soil particles were also investigated.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.8
• /
• pp.2113-2122
• /
• 1994

The 3-D FEM analysis for simulating the stamping operation of planar anisotropic sheet metals with arbitrarily-shaped tools is introduced. An implicit, incremental, updated Lagrangian formulation with a rigid-viscoplastic constitutive equation is employed. Contact and friction are considered through the mesh-normal, which compatibly describes arbitrary tool surfaces and FEM meshes without depending on the explicit spatial derivatives of tool surfaces. The consistent full set of governing relations, comprising equilibrium equation and mesh-normal geometric constraints, is appropriately linearized. The linear triangular elements are used for depicting the formed sheet, based on membrane approximation. Barlat's non-quadratic anisotropic yield criterion(strain-rate potential) is employed, whose in-plane anisotropic properties are taken into account with anisotropic coefficients and non-quadratic function parameter. The planar anisotropic finite element formulation is tested with the numerical simulations of the stamping of an automotive hood inner panel and the drawing of a hemispherical punch. The in-plane anisotropic effects on the formability of both mild steel and aluminum alloy sheet metals are examined.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.15 no.7 s.100
• /
• pp.836-842
• /
• 2005

Nonlinear dynamic characteristics of active piezolaminated plates are investigated with respect to the thermopiezoelastic behaviors. For largely deformed structures with small strain, the incremental total Lagrangian formulation is presented based on the virtual work principles. A multi-field layer-wise finite shell element is proposed for assuring high accuracy and non-linearity of displacement, electric and thermal fields. For dynamic consideration of thermopiezoelastic snap-through phenomena, the implicit Newmark's scheme with the Newton-Raphson iteration is implemented for the transient response of various piezolaminated models with symmetric or eccentric active layers. The bifurcate thermal buckling of symmetric structural models is first investigated and the characteristics of piezoelectric active responses are studied for finding snap-through piezoelectric potentials and the load-path tracking map. The thermoelastic stable and unstable postbuckling, thermopiezoelastic snap-through phenomena with several attractors are proved using the nonlinear time responses for various initial conditions and damping loss factors. Present results show that thermopiezoelastic snap-through phenomena can result in the difficulty of buckling and postbuckling control of intelligent structures.

• Journal of Mechanical Science and Technology
• /
• v.19 no.spc1
• /
• pp.388-394
• /
• 2005

Recently, the authors have developed a method for real-time dynamics of multibody systems, which combines a semi-recursive formulation to derive the equations of motion in dependent relative coordinates, along with an augmented Lagrangian technique to impose the loop closure conditions. The following numerical integration procedures, which can be grouped into the so-called structural integrators, were tested : trapezoidal rule, Newmark dissipative schemes, HHT rule, and the Generalized-${\alpha}$ family. It was shown that, for large multi body systems, Newmark dissipative was the best election since, provided that the adequate parameters were chosen, excellent behavior was achieved in terms of efficiency and robustness with acceptable levels of accuracy. In the present paper, the performance of the described method in combination with another group of integrators, the Implicit Runge-Kutta family (IRK), is analyzed. The purpose is to clarify which kind of IRK algorithms can be more suitable for real-time applications, and to see whether they can be competitive with the already tested structural family of integrators. The final objective of the work is to provide some practical criteria for those interested in achieving real-time performance for large and complex multibody systems.

• Korea-Australia Rheology Journal
• /
• v.21 no.3
• /
• pp.193-200
• /
• 2009

We present a direct simulation technique for two-dimensional mold-filling simulations of fluids filled with a large number of circular disk-like rigid particles. It is a direct simulation in that the hydrodynamic interaction between particles and fluid is fully considered. We employ a pseudo-concentration method for the evolution of the flow front and the DLM (distributed Lagrangian multipliers)-like fictitious domain method for the implicit treatment of the hydrodynamic interaction. Both methods allow the use of a fixed regular discretization during the entire computation. The discontinuous Galerkin method has been used to solve the concentration evolution equation and the rigid-ring description has been introduced for freely suspended particles. A buffer zone, the gate region of a finite area subject to the uniform velocity profile, has been introduced to put discrete particles into the computational domain avoiding any artificial discontinuity. From example problems of 450 particles, we investigated the particle motion and effects of particles on the flow for both Newtonian and shear-thinning fluid media. We report the prolonged particle movement toward the wall in case of a shear-thinning fluid, which has been interpreted with the shear rate distribution.

• Journal of the Korean Society for Precision Engineering
• /
• v.15 no.9
• /
• pp.117-126
• /
• 1998

For investigating rubber pad sheet metal forming process, the rubber pad deformation characteristics as well as the contact problem of rubber pad-sheet metal has been analyzed. In this paper, the behavior of the rubber deformation is represented by hyper-elastic constitutive relations based on a generalized Mooney-Rivlin model. Finite element procedures for the two-dimensional responses, employing total Lagrangian formulations are implemented in an implicit form. The volumetric incompressibility condition of the rubber deformation is included in the formulation by using penalty method. The sheet metal is characterized by elasto-plastic material with strain hardening effect and analyzed by a commercial code. The contact procedure and interface program between rubber pad and sheet metal are implemented. Inflation experiment of circular rubber pad identifies the behaviour of the rubber pad deformation during the process. The various form dies and scaled down apparatus of the rubber-pad forming process are fabricated for simulating realistic forming process. The obtaining experimental data and FEM solutions were compared. The numerical solutions illustrate fair agreement with experimental results. The forming pressure distribution according to the dimensions of sheet metal and rubber pads, various rubber models and rubber material are also compared and discussed.

• Steel and Composite Structures
• /
• v.29 no.6
• /
• pp.749-758
• /
• 2018

This article derived a hybrid coupling technique using the higher-order displacement polynomial and three soft computing techniques (teaching learning-based optimization, particle swarm optimization, and artificial bee colony) to predict the optimal stacking sequence of the layered structure and the corresponding frequency values. The higher-order displacement kinematics is adopted for the mathematical model derivation considering the necessary stress and stain continuity and the elimination of shear correction factor. A nine noded isoparametric Lagrangian element (eighty-one degrees of freedom at each node) is engaged for the discretisation and the desired model equation derived via the classical Hamilton's principle. Subsequently, three soft computing techniques are employed to predict the maximum natural frequency values corresponding to their optimum layer sequences via a suitable home-made computer code. The finite element convergence rate including the optimal solution stability is established through the iterative solutions. Further, the predicted optimal stacking sequence including the accuracy of the frequency values are verified with adequate comparison studies. Lastly, the derived hybrid models are explored further to by solving different numerical examples for the combined structural parameters (length to width ratio, length to thickness ratio and orthotropicity on frequency and layer-sequence) and the implicit behavior discuss in details.