• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.12 no.3
• /
• pp.465-474
• /
• 1999

A recently developed numerical method based on a volume integral formulation is applied to calculate the accurate stress intensity factors at the crack tips in unbounded isotropic solids in the presence of multiple anisotropic inclusions and cracks subject to external loads. In this paper, a detailed analysis of the stress intensity factors are carried out for an unbounded isotropic matrix containing an orthotropic cylindrical inclusion and a crack. The accuracy and effectiveness of the new method are examined through comparison with results obtained from analytical method and finite element method using ANSYS. It is demonstrated that this new method is very accurate and effective for solving plane elastostatic problems in unbounded solids containing anisotropic inclusions and cracks.

• Journal of the Society of Naval Architects of Korea
• /
• v.42 no.4 s.142
• /
• pp.388-395
• /
• 2005

A numerical procedure is proposed as a mesh-size insensitive structural stress definition that gives a stress state at a weld toe with relatively large mesh size. The structural stress values obtained using different finite element types, i.e. shell element and solid element, are examined for typical weld structures. The calculation procedures are performed using the balanced nodal forces and moments obtained from finite element solutions. A consistent formulation based on work equivalent argument has been implemented to transform the balanced nodal forces and moments from shell to line force and line moments at each nodal position. The mesh-insensitivity, the effect of distance $\delta$(where the stress is calculated) and the potential limitations of the structural stress method are examined for various types of weldments. Based on the results from this study, it is expected to develop a more precise stress estimation technique for fatigue strength assessment of welded structures.

• Journal of the Society of Naval Architects of Korea
• /
• v.42 no.4 s.142
• /
• pp.357-367
• /
• 2005

A numerical method is developed for the hydroelastic response of the Very Large Floating Structure considering the sea-bottom topography. The sea-bottom effects on the hydroelastic response of the floating structure is studied. The sea-bottom topography should be considered when the floating structure is constructed near the shore. To investigate the sea-bottom effects, four different sea-bottom topographies are considered in this study. finite-element method based on the variational formulation is used in the fluid domain, The pontoon-type floating structure is modeled as the Kirchhoff plate. The mode superposition method is adopted for the hydroelastic behavior of the floating structure.

• Journal of the Society of Naval Architects of Korea
• /
• v.42 no.4 s.142
• /
• pp.330-340
• /
• 2005

Compressible two-phase flow is analyzed based on the arbitrary Lagrangian-Eulerian (ALE) formulation. For water, Tamman type stiffened equation of state is used. Numerical fluxes are calculated using the ALE two-phase Godunov scheme which assumes only that the speed of sound and pressure can be provided whenever density and internal energy are given. Effects of the approximations of a material interface speed are Investigated h method Is suggested to assign a rigid body boundary condition effectively To validate the developed code, several well-known problems are calculated and the results are compared with analytic or other numerical solutions including a single material Sod shock tube problem and a gas/water shock tube problem The code is applied to analyze the refraction and transmission of shock waves which are impacting on a water-gas interface from gas or water medium.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.6 no.1
• /
• pp.51-60
• /
• 1994

A one-dimensional numerical model has been developed through a Lagragian formulation, which is applicable to advection and diffusion of dissolved matters in storages. The study was conducted to the pollutants released into a canal network in Burnt Store Isles, Punta Gorda, Florida. USA. The hydrodynamic model was developed by using an implicit finite difference scheme. In the computational domain, the network system consists of prismatic channels and storages. The finger canals and small tributaries ale treated as storages in the simplified flow network The numerical results show relatively good agreement with field experiments.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.9 s.240
• /
• pp.1049-1056
• /
• 2005

A conservative pressure-based finite-volume numerical method has been developed for computing flow and heat transfer by using an unstructured grid system. The method admits arbitrary convex polyhedra. Care is taken in the discretization and solution procedures to avoid formulations that are cell-shape-specific. A collocated variable arrangement formulation is developed, i.e. all dependent variables such as pressure and velocity are stored at cell centers. Gradients required for the evaluation of diffusion fluxes and for second-order-accurate convective operators are found by a novel second-order accurate spatial discretization. Momentum interpolation is used to prevent pressure checkerboarding and the SIMPLE algorithm is used for pressure-velocity coupling. The resulting set of coupled nonlinear algebraic equations is solved by employing a segregated approach, leading to a decoupled set of linear algebraic equations fer each dependent variable, with a sparse diagonally dominant coefficient matrix. These equations are solved by an iterative preconditioned conjugate gradient solver which retains the sparsity of the coefficient matrix, thus achieving a very efficient use of computer resources.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.28 no.4
• /
• pp.395-405
• /
• 2004

The filling pattern and an adaptive grid refinement based on the finite element method and Eulerian mesh advancement approach have been developed to analyze incompressible transient viscous flow with free surfaces. The governing equation for flow analysis is Navier-Stokes equation including inertia and gravity effects. The mixed FE formulation and predictor-corrector method are used effectively for unsteady numerical simulation. The flow front surface and the volume inflow rate are calculated using the filling pattern technique to select an adequate pattern among four filling patterns at each triangular control volume. By adaptive grid refinement, the new flow field that renders better prediction in flow surface shape is generated and the velocity field at the flow front part is calculated more exactly. In this domain the elements in the surface region are made finer than those in the remaining regions for more efficient computation. Using the proposed numerical technique, the collapse of a water dam has been analyzed to predict flow phenomenon of fluid and the predicted front positions with respect to time have been compared with the reported experimental results.

• KSTLE International Journal
• /
• v.3 no.2
• /
• pp.90-94
• /
• 2002

A press shoe is an element of a machine for squeezing water from wood pulp in the field of manufacturing paper. This is used to compress the pulp enveloped by felt sheet with a large roller. The squeezing farce is made by hydraulic pressure. The press shoe has a mechanism similar to a partial hydrostatic bearing. The pressure profile between press shoe and roller affects their squeezing ability, and partial peak pressure can tear the wet pulp. The curvature of the surface of press shoe varies to reduce the peak pressure and increase the mean pressure simultaneously, Therefore, the prediction of pressure distribution considering partially changed curvature of hydrostatic bearing is very important far designing the press shoe. In this study, the difference formulation of Reynolds equation far partial hydrostatic bearing is derived by direct numerical method and a computer program to calculate the pressure distribution is developed. We investigate the effect of partially changed curvature of bearing surface on the pressure distribution. Other design parameter far hydrostatic bearing such as depth of pocket and relative velocity are also studied.

• Tribology and Lubricants
• /
• v.16 no.1
• /
• pp.39-45
• /
• 2000

In this paper, the effects of shaft and bearing flexibilities are investigated for the accurate modeling of a shaft-bearing system supported by ball bearings. Generally, rolling bearings are modeled by simple rigid pin-joint in the mechanical design. However, they can no longer be modeled by ideal boundary conditions in the advanced applications because the rigid pin-joint model cannot satisfy the current trends of mechanical design decreasing mass and reducing volume. Consequently the flexible support model of ball bearing is investigated using the static analysis module developed by A .B. Jones and T. A. Harris. A simple two-bearing system, supported by two deep groove ball bearings and radially loaded on the shaft midway between the bearings, is utilized to validate the coupled model of shaft-bearing system. Numerical computations using the model indicate that the shaft span length, locating/floating bearing arrangements and applied bearing size are significant factors in determining the mechanical behaviors. The flexible support model of ball bearing can be escaped to over-estimate in the bearing fatigue life. The proposed simple design formulation obtained by numerical simulations can approximately predict a rate of bearing life reduction as a function of shaft span length/shaft diameter (L/d).

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 1996.06a
• /
• pp.97-109
• /
• 1996

In the present work computations are carried out for analysis of complicated sheet metal forming process such as forming of a rear hinge. Finite element formulation using dynamic explicit time integration scheme and step-wise combined Implicit/Explicit scheme are introduced for numerical analysis of sheet metal forming process. The rigid-plastic finite element method based on membrane elements has long been employed as a useful numerical technique for the analysis of sheet metal forming because of its time effectiveness. The explicit scheme in general use is based on the elastic-plastic modelling of material requiring large computation time. In finite element simulation of sheet metal forming processes, the robustness and stability of computation are important requirements since the computation time and convergency become major points of consideration besides the solution accuracy due to the complexity of geometry and boundary conditions. The implicit scheme employs a more reliable and rigorous scheme in considering the equilibrium at each step of deformation, while in the explicit scheme the problem of convergency is eliminated at the cost of solution accuracy. The explicit approach and the implicit approach have merits and demerits, respectively. In order to combine the merits of these two methods a step-wise combined implicit/explicit scheme has been developed.