• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 1994.10a
• /
• pp.44-54
• /
• 1994

A variational formulation and the associated finite elemet equations have been derived for general three-dimensional deformation of a planar anisotropic rigid-plastic sheet metal which obeys the strain-rate potential proposed by BARLAT et al [13]. By using the natural convected coordinate system, the effect of geometric change and the rotation of planar anisotropic axes are considered efficiently. In order to check the validity of present formulation, a cylindrical cup and a square cup deep drawing test was modeled. good agreement was found between the FE simulation and the experiment. The results have shown that the present formulation for planar anisotropic deformation can be efficiently applied to the analysis of sheet metal working processes for planar anisotropic nonferrous metals.

• Proceedings of the KIEE Conference
• /
• 2000.07b
• /
• pp.903-905
• /
• 2000

For an accurate analysis of three dimensional linear magnetostatic problems, a new boundary integral equation formulation is presented. This formulation adopts difference magnetic field concept and uses single layer magnetic surface charge as unknown. The proposed method is capable of eliminating numerical cancellation errors inside ferromagnetic materials. In additions, computing time and storage memory are reduced by 75% in comparison with the reduced and total scalar potential formulation. Two examples are given to show its efficiency and accuracy.

• Coupled systems mechanics
• /
• v.1 no.1
• /
• pp.59-78
• /
• 2012

The available analytical methods of analysis for laterally loaded piles in level ground cannot be directly applied to such piles in sloping ground. With the commercially available software, the simulation of the appropriate field condition is a challenging task, and the results are subjective. Therefore, it becomes essential to understand the process of development of a user-framed numerical formulation, which may be used easily as per the specific site conditions without depending on other indirect methods of analysis as well as on the software. In the present study, a detailed three-dimensional finite element formulation is presented for the analysis of laterally loaded piles in sloping ground developing the 18 node triangular prism elements. An application of the numerical formulation has been illustrated for the pile located at the crest of the slope and for the pile located at some edge distance from the crest. The specific examples show that at any given depth, the displacement and bending moment increase with an increase in slope of the ground, whereas they decrease with increasing edge distance.

• Interaction and multiscale mechanics
• /
• v.3 no.3
• /
• pp.213-234
• /
• 2010

A multiscale method is presented for analysis of thin slab structures in which the microstructures can not be reduced to two-dimensional plane stress models and thus three dimensional treatment of microstructures is necessary. This method is based on the classical asymptotic expansion multiscale approach but with consideration of the special geometric characteristics of the slab structures. This is achieved via a special form of multiscale asymptotic expansion of displacement field. The expanded three dimensional displacement field only exhibits in-plane periodicity and the thickness dimension is in the global scale. Consequently by employing the multiscale asymptotic expansion approach the global macroscopic structural problem and the local microscopic unit cell problem are rationally set up. It is noted that the unit cell is subjected to the in-plane periodic boundary conditions as well as the traction free conditions on the out of plane surfaces of the unit cell. The variational formulation and finite element implementation of the unit cell problem are discussed in details. Thereafter the in-plane material response is systematically characterized via homogenization analysis of the proposed special unit cell problem for different microstructures and the reasoning of the present method is justified. Moreover the present multiscale analysis procedure is illustrated through a plane stress beam example.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.03b
• /
• pp.2-5
• /
• 2008

This study describes the numerical simulation of three-dimensional droplet formation and the following motion in a cross-junction microchannel by using the Lattice Boltzmann Method (LBM). Our aim is to develop the three-dimensional binary fluids model, consisting of two sets of distribution functions to represent the total fluid density and the density difference, which introduces the repulsive interaction consistent with a free-energy function between two fluids. We validated the LBM code with the velocity profile in a 3-dimensional rectangular channel. Then, we applied our code to the numerical simulation of a binary fluid flow in a cross-junction channel focusing on the investigation of the droplet formulation. Due to the pressure and interfacial-tension effect, one component of the fluids which is injected from one inlet is cut off into many droplets periodically by the other component which is injected from the other inlets. We considered the effect of the boundary conditions for density difference (order parameter) on the wetting of the droplet to the side walls.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2008.03a
• /
• pp.597-601
• /
• 2008

Three-dimensional structures of detonation wave propagating in circular tube were investigated. Inviscid fluid dynamics equations coupled with a conservation equation of reaction progress variable were analyzed by a MUSCL-type TVD scheme and four stage Runge-Kutta time integration. Variable-$\gamma$ formulation was used to account for the variable properties between unburned and burned states and the chemical reaction was modeled by using a simplified one-step irreversible kinetics model. The computational code was parallelized based on domain decomposition technique using MPI-II message passing library. The computations were carried out using a home made Windows based PC cluster having 160 AMD AthloxXP and Athlon64 processor. The computational domain consisted of through a roundshaped tube with wall conditions. As an initial condition, analytical ZND solution was distributed over the computational domain with disturbances. The disturbances has circumferential large gradient. The unsteady computational results in three-dimension show the detailed mechanisms of multi-cell mode of detonation wave instabilities resulting diamond shape in smoked-foil record.

• Structural Engineering and Mechanics
• /
• v.17 no.1
• /
• pp.113-140
• /
• 2004

A mixed eight-node hexahedral element formulated via the Hu-Washizu principle as well as the field extrapolation technique is presented. The mixed element with only three translational degrees of freedom at each node can provide extremely accurate and reliable performance for popular benchmark problems such as spacial beams, plates, shells as well as general three-dimensional elasticity problems. Numerical calculations also show that when extremely skewed and coarse meshes and nearly incompressible materials are used, the proposed mixed element can still possess excellent behaviour. The mixed formulation starts with introduction of a parallelepiped domain associated with the given general eight-node hexahedral element. Then, the assumed strain field at the nodal level is constructed via the Hu-Washizu variational principle for that associated parallelepiped domain. Finally, the assumed strain field at the nodal level of the given hexahedral element is established by using the field extrapolation technique, and then by using the trilinear shape functions the assumed strain field of the whole element domain is obtained. All matrices involved in establishing the element stiffness matrix can be evaluated analytically and expressed explicitly; however, a 24 by 24 matrix has to be inverted to construct the displacement extrapolation matrix. The proposed hexahedral element satisfies the patch test as long as the element with a shape of parallelepiped.

• 한국전산유체공학회:학술대회논문집
• /
• 2005.10a
• /
• pp.117-124
• /
• 2005

Splitting algorithms of the incompressible Navier-Stokes equations using P1P1/P2P1 finite element formulation are newly proposed. P1P1 formulation allocates velocity and pressure at the same nodes, while P2P1 formulation allocates pressure only at the vertex nodes and velocity at both the vertex and mid nodes. For comparison of the elapsed time and accuracy of the two methods, they have been applied to the well-known benchmark problems. The three cases chosen are the two-dimensional steady and unsteady flows around a fixed cylinder, decaying vortex, and impinging slot jet. It is shown that the proposed P2P1 semi-splitting method performs better than the conventional P1P1 splitting method in terms of both accuracy and computation time.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.30 no.3 s.246
• /
• pp.262-269
• /
• 2006

Segregation algorithms of the incompressible Wavier-Stokes equations using P1P1/P2P1 finite element formulation are newly proposed. P1P1 formulation allocates velocity and pressure at the same nodes, while P2P1 formulation allocates pressure only at the vertex nodes and velocity at both the vertex and the midpoint nodes. For a comparison of both the elapsed time and the accuracy between the two methods, they have been applied to the well-known benchmark problems. The three cases chosen are the two-dimensional steady and unsteady flows around a fixed cylinder, decaying vortex, and impinging slot jet. It is shown that the proposed P2P1 semi-segregation algorithm performs better than the conventional P1P1 segregation algorithm in terms of both accuracy and computation time.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.22 no.3
• /
• pp.316-325
• /
• 1998

The objective of this study is to investigate the role of large-scale coherent structures in a spatially developing plane mixing layers. To achieve this, we have to look into the mutual interactions between three-dimensional large-scale coherent structures and the mean flow. Our attention will be focused on the energy exchange mechanism between the various modes, and the effects of the nonlinear evolution of the phases of the interacting modes. Linear stability of the three-dimensional viscous shear layer is formulated and solved as the basis for the solution of the nonlinear formulation based on the energy method. The importance of the initial conditions that may affect the evolution of the flow has been examined. It has been numerically calculated the nonlinear effects arising from the interactions among the three-dimensional large-scale coherent structures in a spatially developing plane mixing layers. The results of this study provide useful parametric information for the control of shear layer in practical applications in the mixing and transport augmentation.