• Proceedings of the Korean Geotechical Society Conference
• /
• 2000.11a
• /
• pp.665-672
• /
• 2000

The problems of discontinuous layer interface are very important in the algorithm and programming for the analysis of multi-layered consolidation using a numerical analysis, finite difference method(F.D,M.). Better results can be obtained by the process for discontinuous layer interface, since it can help consolidation analysis to model the actual ground Explicit method is simple for analysis algorithm and convenient for use except for applying the operator Crank-Nicolson method represents implicit method, which have different analysis method according to weighting factor. This method uses different algorithm according to dimension. And, this paper uses alternative direction implicit method. The purpose of this paper provides an efficient computer algorithm based on numerical analysis using finite difference method which account for multi-layered soils with confined aquifer to determine the degree of consolidation and excess pore pressures relative to time and positions more realistically.

• Geophysics and Geophysical Exploration
• /
• v.3 no.2
• /
• pp.39-47
• /
• 2000

A great deal of computing time and a large computer memory are needed to solve wave equation in a large complex subsurface layers using the finite difference method. The computing time and memory can be reduced by decreasing the number of grid points per minimum wave length. However, the decrease of grids may cause numerical dispersion and poor accuracy. In this study we performed the grid dispersion analysis for several rotated finite difference operators, which was commonly used to reduce grids per wavelength with accuracy in order to determine the solution for the acoustic wave equation in frequency domain. The rotated finite difference operators were to be extended to 81, 121 and 169 difference stars and studied whether the minimum grids could be reduced to 2 or not. To obtain accuracy (numerical errors less than $1\%$) the following was required: more than 13 grids for conventional 5 point difference stars, 9 grids for 9 difference stars, 3 grids for 25 difference stars, and 2.7 grids for 49 difference stars. After grid dispersion analysis for the new rotated finite difference operators, more than 2.5 grids for 81 difference stars, 2.3 grids for 121 difference stars and 2.1 grids for 169 difference stars were needed. However, in the 169 difference stars, there was no solution because of oscillation of the dispersion curves in the group velocity curves. This indicated that the grids couldn't be reduced to 2 in the frequency acoustic wave equation. According to grid dispersion analysis for the determination of grid points, the more rotated finite difference operators, the fewer grid points. However, the more rotated finite difference operators that are used, the more complex the difference equation terms.

• Transactions of Materials Processing
• /
• v.16 no.8
• /
• pp.590-595
• /
• 2007

This paper introduces a new approach to take account of bending history in finite element inverse analysis during sheet metal forming process. A modified membrane element was adopted for finite element inverse analysis so that bending-unbending energy was additionally imposed in the total plastic energy, predicting bending-unbending regions using the geometry of the final shape and tools. An algorithm was applied to a cylindrical cup deep drawing process. The blank shape and the distribution of the thickness strain were compared with those obtained from the incremental finite element analysis in order to evaluate the effect of the bending history. The algorithm reduced the difference between the results of the inverse analysis from those of the incremental analysis due to bending history. The analysis was also carried out with the variation of the thickness of the initial blank to investigate the effect of bending deformation. The results showed that the difference was remarkably reduced as the thickness of the initial blank increased. This indicates that the finite element inverse analysis cooperated with the suggested scheme is useful to obtain more accurate results, especially when bending effects are significant.

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.21 no.2
• /
• pp.211-225
• /
• 2019

Recently, earthquakes have occurred throughout the entire region of Korea and seismic analysis studies have been actively conducted in various fields. SSI analyses studies considering ground have been carried out consistently. However, few comparative analyses have been performed on the dynamic behavior of buildings according to numerical analysis method in the case of the previous dynamic analyses considering grounds. Therefore, in this study, the dynamic analyses were performed on a high-rise building by using both a finite element program MIDAS GTS NX and a finite difference program FLAC 2D. The results were compared and analyzed each other. As a result, both the maximum compressive and tensile bending stresses of above ground and below ground part were estimated to be a little larger by MIDAS GTS NX than by FLAC 2D. However, the maximum horizontal displacement value, the horizontal displacement distribution, and the position of weak part were turned out to be similar in both analysis programs. Therefore, it can be concluded that there is no difference in using either a finite element program or a finite difference program for the convenience of a user for a dynamic analysis.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.26 no.11
• /
• pp.2245-2252
• /
• 2002

Design sensitivity analysis scheme is proposed in an elasto -plastic finite element method with explicit time integration using a direct differentiation method. The direct differentiation is concerned with large deformation, the elasto-plastic constitutive relation, shell elements with reduced integration and the contact scheme. The design sensitivities with respect to the process parameter are calculated with the direct analytical differentiation of the governing equation. The sensitivity results obtained from the present theory are compared with that obtained by the finite difference method in a class of sheet metal forming problems such as hemi-spherical stretching and cylindrical cup deep-drawing. The result shows good agreement with the finite difference method and demonstrates that the preposed sensitivity calculation scheme is a pplicable in the complicated sheet metal forming analysis and design.

• Structural Engineering and Mechanics
• /
• v.46 no.2
• /
• pp.269-294
• /
• 2013

This paper is concerned with the determination of exact buckling loads and vibration frequencies of variable thickness isotropic plates using well known finite difference technique. The plates are subjected to uni, biaxial compression and shear loadings and various combinations of boundary conditions are considered. The buckling load is found out as the in plane load that makes the determinant of the stiffness matrix equal to zero and the natural frequencies are found out by carrying out eigenvalue analysis of stiffness and mass matrices. New and exact results are given for many cases and the results are in close agreement with the published results. In this paper, like finite element method, finite difference method is applied in a very simple manner and the application of boundary conditions is also automatic.

• Proceedings of the IEEK Conference
• /
• 2000.11a
• /
• pp.489-492
• /
• 2000

Film Bulk Acoustic Wave Resonator(FBAR) used in microwave region was analyzed with Finite Difference Time-Domain Methods(FDTD) in this paper. FBAR have been analyzed with one dimensional Mason model analysis or Finite Element methods(FEM), but the first couldn't analyze effect of area variation and spurious characteristics, the second had difficulty in element separation because of thin electrode. So in this paper FBAR was analyzed by Finite Difference Time-Domain Methods and it's results were transformed to frequency domain using Discrete Fourier Transform.

• Journal of the Korean Society for Precision Engineering
• /
• v.18 no.9
• /
• pp.110-121
• /
• 2001

During hot pipe bending using induction heating, the wall of bending outside is thinned by tensile stress. In design requirement, the reduction of wall thickness is not allowed to exceed 12.5%. So in this study, two methods of bending, one is loading of reverse moment and the other is loading of temperature gradient, have been investigated to design pipe bending process that satisfy design requirements. For this purpose, finite element analysis with a bending radius 2Do(outer diameter of pipe) has been performed to calculate proper reverse moment and temperature gradient to be applied. Induction heating process has been analyzed to estimate influence of heating process parameters on heating characteristic by finite difference method. Then pipe bending experiments have been performed for verification of finite element and finite difference analysis results. Experimental results are in good agreement with the results of simulations.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.7 no.4
• /
• pp.320-327
• /
• 1996

This paper presents that the lossy or amplification property of the Finite Difference-Time Domain(FD-TD) method based on the leap-frog scheme is theoretically verified by using a plane wave analysis. The basic algorithm of the FD-TD method is introduced in order to help understanding the analysis procedure. Since our analysis is formulated by the Von Neumann's approach, the stability inequality is also produced as an another outcome.

• Journal of Ocean Engineering and Technology
• /
• v.16 no.5
• /
• pp.15-20
• /
• 2002

In this research, the finite difference lattice Boltzmann method(FDLBM) is used to analyze gravity currents in the lock exchange configuration that occur in many natural and man-made situations. At a lock those are seen when a gate is suddenly opened, and, in the atmosphere, when the thunderstorm outflows make a cold front. At estuaries in the ocean, the phenomenon is found between fresh water from a river and salt water in the sea. Since such interesting phenomena were recognized, pioneers have challenged to make them clear by conducing both experiments and analysis. Most of them were about the currents of liquid or Boussinesq fluids, which are assumed as incompressible. Otherwise, the difference in density of two fluids is small. The finite difference lattice Boltzmann method has been a powerful tool to simulate the flow of compressible fluids. Also, numerical predictions using FDLBM to clarify the gravity currents of compressible fluids exhibit all features, but typically observed in experimental flows near the gravity current head, including the lobe-and-cleft structure at the leading edge.