• Structural Engineering and Mechanics
• /
• v.29 no.5
• /
• pp.469-488
• /
• 2008

Finite element methods have often been used for structural analyses of various mechanical problems. When finite element analyses are utilized to resolve mechanical systems, numerical uncertainties in the initial data such as structural parameters and loading conditions may result in uncertainties in the structural responses. Therefore the initial data have to be as accurate as possible in order to obtain reliable structural analysis results. The typical finite element method may not properly represent discrete systems when using uncertain data, since all input data of material properties and applied loads are defined by nominal values. An interval finite element analysis, which uses the interval arithmetic as introduced by Moore (1966) is proposed as a non-stochastic method in this study and serves a new numerical tool for evaluating the uncertainties of the initial data in structural analyses. According to this method, the element stiffness matrix includes interval terms of the lower and upper bounds of the structural parameters, and interval change functions are devised. Numerical uncertainties in the initial data are described as a tolerance error and tree graphs of uncertain data are constructed by numerical uncertainty combinations of each parameter. The structural responses calculated by all uncertainty cases can be easily estimated so that structural safety can be included in the design. Numerical applications of truss and frame structures demonstrate the efficiency of the present method with respect to numerical analyses of structural uncertainties.

• Proceedings of the Korean Geotechical Society Conference
• /
• 2002.03a
• /
• pp.197-204
• /
• 2002

In this study, multivariate analysis based on domestic data(958 EA) of road tunnel, and suggest the easy prediction equation of Q-system. We generate applicable Q-value to numerical analysis method with using the equation and investigate the behavior as variable Q-value of rock mass induced excavation with discontinuum numerical analysis method, UDEC. In the result of the experiment, we research the application range of Q-value to apply the continuum model to discontinuous rock mass is below 0.7 and we testify the applicability of continuum model as researched Q-value with continuum numerical analysis method, FLAC.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.7 no.6
• /
• pp.1034-1043
• /
• 2015

This is the continuation of a numerical study on vortex shedding from a blunt trailing-edge of a hydrofoil. In our previous work (Lee et al., 2015), numerical schemes for efficient computations were successfully implemented; i.e. multiple domains, the approximation of domain boundary conditions using cubic spline functions, and particle-based domain decomposition for better load balancing. In this study, numerical results through a hybrid particle-mesh method which adopts the Vortex-In-Cell (VIC) method and the Brinkman penalization model are further rigorously validated through comparison to experimental data at the Reynolds number of $2{\times}10^6$. The effects of changes in numerical parameters are also explored herein. We find that the present numerical method enables us to reasonably simulate vortex shedding phenomenon, as well as turbulent wakes of a hydrofoil.

• Tunnel and Underground Space
• /
• v.2 no.1
• /
• pp.152-163
• /
• 1992

The study of rock mass behaviour through a numerical analysis is important for the design, construction and maintenance of large rock caverns. The objectives of the numerical analysis are to design reasonably and construct safely the underground structures, to maintain them soundly after construction and to extend them securely for a desired period of time. Methods of numerical analyses included in this case study are the finite element method, the boundary element method, and the distinct element method. The numerical models are purely elastic, elastoplastic, visco-elastic, visco-plastic, easto-visco-plastic and jointed-discontinuous materials. The results of this case study indicate that the rock mass behaviour could be predicted exactly through continuous comparisons of the numerical results with the in-situ measurements.

• Journal of Ship and Ocean Technology
• /
• v.5 no.4
• /
• pp.29-43
• /
• 2001

In this paper, free surface flows around an advancing two-dimensional rectangular cylinder piercing the free surface are studied using numerical and experimental methods. Especially, wave breaking phenomenon around the cylinder is treated in detail. A series of numerical simulations and experiments were performed for the purpose of comparison. For the numerical simulations, a finite difference method was adopted with a rectangular grid system, and the variation of the free surface was computed by the marker density method. The computational results are compared with the experiments. It is confirmed that the present numerical method is useful for the numerical simulation of nonlinear free surface waves around a piercing body.

• Proceedings of the Korean Institute of Intelligent Systems Conference
• /
• 1993.06a
• /
• pp.1230-1233
• /
• 1993

Numerical solution of inverse problem of Takagi-Sugeno fuzzy model is proposed. The method is located on the application of numerical optimization to the fuzzy model. Steepest descent method is used for the numerical optimization. We use the linear approximation of fuzzy model, called pseudo first order approximation, by fixing the membership value on the neighborhood of the corresponding input. It is introduced in order to reduce the difficulty of optimization process. The efficiency of this method is shown by a numerical experiment.

• Journal of computational fluids engineering
• /
• v.7 no.1
• /
• pp.10-19
• /
• 2002

The non-staggered(collocated) grid approach in which all the solution variables are located at the centers of control volumes is very popular for incompressible flow analyses because of its numerical efficiency on the curvilinear or unstructured grids. Rhie and Chow's paper is the first in using non-staggered grid method for SIMPLE algorithm, where pressure weighted interpolation was used to prevent decoupling of pressure and velocity. But it has been known that this non-staggered grid method has stability problems when pressure fields are nonlinear like in natural convection flows. Also Rhie-Chow scheme generates large numerical diffusion near curved walls. The cause of these unwanted problems is too large pressure damping term compared to the magnitude of face velocity. In this study the magnitude of pressure damping term of Rhie-Chow's method is limited to 1∼10% of face velocity to prevent physically unreasonable solutions. The wall pressure extrapolation which is necessary for cell-centered FVM is another source of numerical errors. Some methods are applied in a unstructured FV solver and analyzed in view of numerical accuracy. Here, two natural convection problems are solved to check the effect of the Rhie-Chow's method on numerical stability. And numerical diffusion from Rhie-Chow's method is studied by solving the inviscid flow around a circular cylinder.

• Journal of Korean Society of Steel Construction
• /
• v.10 no.3 s.36
• /
• pp.451-461
• /
• 1998

The solution of anisotropic plate via the classical methods is limited to relatively load and boundary conditions. If these conditions are more complex, the analysis becomes increasingly tedious and even impossible. For many plate problems of considerable practical interest, analytic solutions to the governing differential equations cannot be found. Among the numerical techniques presently available, the finite difference method and the finite element method are powerful numerical methods. The objective of this paper is to compare with each numerical methods for the buckling load and modes of anisotropic composite laminated plates considering shear deformation. In applying numerical methods to solve differential equations of anisotropic plates, this study uses the finite difference method and the finite element method. In determining the eigenvalue by Finite Difference Method, this paper represent good convergence compared with Finite Element Method. Several numerical examples and buckling modes show the effectiveness of various numerical methods and they will give a guides in deciding minimum buckling load and various mode shapes.

• Journal of the Korean Society for Precision Engineering
• /
• v.12 no.3
• /
• pp.32-41
• /
• 1995

In the analysis of metal forming processes by the finite element method, there are many numerical instabilities such as element locking, hourglass mode and shear locking. These instabilities may have a bad effect upon accuracy and convergence. The present work is concerned with improvement of stability and efficiency in two-dimensional rigid-plastic finite element method using various type of elemenmts and numerical intergration schemes. As metal forming examples, upsetting and backward extrusion are taken for comparison among the methods: various element types and numerical integration schemes. Comparison is made in terms of stability and efficiency in element behavior and computational efficiency and a new scheme of adaptive directional reduced integration is introduced. As a result, the finite element computation has been stabilized from the viewpoint of computational time, convergency, and numerical instability.

• Journal of Industrial Technology
• /
• v.10
• /
• pp.69-72
• /
• 1990

In this paper the numerical method for the study of wave reflection from and transmission through submerged permeable breakwaters using the boundary element method is developed. The numerical analysis technique is based on the wave pressure function instead of velocity potential because it is difficult to define the velocity potential in the each region arising the energy dissipation. Also, the non-linear energy dissipation within the submerged porous structure is simulated by introducing the linear dissipation coefficient and the tag mass coefficient equivalent to the non-linear energy dissipation. For the validity of this analysis technique, the numerical results obtained by the present boundary element method are compared with those obtained by the other computation method. Good agreements are obtained and so the validity of the present numerical analysis technique is proved.