• Journal of the Society of Naval Architects of Korea
• /
• v.29 no.1
• /
• pp.14-28
• /
• 1992

For a planar curve-sided paned with constant or linear density distributions of source or doublet in the singularity methods, Cantaloube and Rehbach(1986) show that the surface integral can be transformed into contour integral by using Stokes' formulas. As an extension of their formulations, this paper deals with a planar polygonal panel for which we derive the closed-forms of the potentials and the velocities induced by the singularity distributions. Test calculations show that the analytical evaluation of the closed-forms is superior to numerical integration(suggested by Cantaloube and Rehbach) of the contour integral. The compact and explicit expressions may produce accurate values of matrix elements of simultaneous linear equations in the singularity methods with much reduced computer tiome.

• Journal of the Society of Naval Architects of Korea
• /
• v.37 no.1
• /
• pp.67-81
• /
• 2000

An accurate and efficient numerical method for two-dimensional nonlinear radiation problem has been developed. The wave motion due to a moving body is described by the assumption of ideal fluid flow, and the governing Laplace equation can be effectively solved by the higher-order boundary element method with the help of the GMRES (Generalized Minimal RESidual) algorithm. The intersection or corner problem is resolved by utilizing the so-called discontinuous elements. The implicit trapezoidal rule is used in updating solutions at new time steps by considering stability and accuracy. Traveling waves caused by the oscillating body are absorbed downstream by the damping zone technique. It is demonstrated that the present method for time marching and radiation condition works efficiently for nonlinear radiation problem. To avoid the numerical instability enhanced by the local gathering of grid points, the regriding technique is employed so that all the grids on the free surface may be distributed with an equal distance. This makes it possible to reduce time interval and improve numerical stability. Special attention is paid to the local flow around the body during time integration. The nonlinear radiation force is calculated by the "acceleration potential technique". Present results show good agreement with other numerical computations and experiments.

• Journal of the Korean Geotechnical Society
• /
• v.22 no.9
• /
• pp.45-59
• /
• 2006

The successful performance of any numerical geotechnical simulation depends on the accuracy and efficiency of the numerical implementation of constitutive model used to simulate the stress-strain (constitutive) response of the soil. The corner stone of the numerical implementation of constitutive models is the numerical integration of the incremental form of soil-plasticity constitutive equations over a discrete sequence of time steps. In this paper a well known two-surface soil plasticity model is implemented using a generalized implicit return mapping algorithm to arbitrary convex yield surfaces referred to as the Closest-Point-Projection method (CPPM). The two-surface model describes the nonlinear behavior of coarse-grained materials by incorporating a bounding surface concept together with isotropic and kinematic hardening as well as fabric formulation to account for the effect of fabric formation on the unloading response. In the course of investigating the performance of the CPPM integration method, it is proven that the algorithm is an accurate, robust, and efficient integration technique useful in finite element contexts. It is also shown that the algorithm produces a consistent tangent operator $\frac{d\sigma}{d\varepsilon}$ during the iterative process with quadratic convergence rate of the global iteration process.

• Journal of the Korean Society of Fisheries and Ocean Technology
• /
• v.39 no.1
• /
• pp.27-32
• /
• 2003

Under the assumption of potential flow, free-surface flow around a hydrofoil is calculated by the high-order spectra1!boundary-integral method, This method is one of the most efficient numerical methods by which the nonlinear interactions between hydrofoil and free-surface can be simulated in time-domain. In this method. the wave potential which represents the nonlinear evolution of free-surface is solved by the high-order spectral method and the body potential which provides the effects of hydrofoil and shed vortex is solved by the boundary-integral method. The calculated free-surface profiles which are generated by a uniformly translating hydrofoil are compared with other experimental results. And they show relatively good agreements each other. As another example, free-surface flow generated by a heaving and translating hydrofoil is calculated and discussed.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.12 no.4
• /
• pp.537-549
• /
• 1999

선형탄성 파괴해석은 균열을 갖는 변형률 경화재료의 파괴거동을 예측하는데 불충분하기 때문에 최근에는 균열 선단 부에서 대규모 소성 역을 갖는 균열 체에 적용할 수 있는 많은 파괴역학개념이 제안되고 있다. 따라서, 본 연구에서는 대규모항복 조건하의 연성파괴를 보이는 평판을 정확하게 해석할 수 있는 새로운 유한요소모델을 제시하고자 한다. 균열 선단 부의 응력 장을 정의하는데 가장 지배적인 파괴매개변수인 J-적분 값과 소성 역의 크기 및 형상을 J-적분법과 등가영역적분법을 통해 파괴거동을 설명할 수 있도록 증분소성이론에 기초를 둔 p-version 유한요소해석이 채택되었다. 제안된 유한요소모델에 의한 수치해석결과는 이론 해와 h-version 유한요소해석과 비교되었다.

• Computational Structural Engineering
• /
• v.6 no.2
• /
• pp.103-114
• /
• 1993

An advanced time-domain finite element formulation is presented for the displacement and stress analysis of isotropic, linear viscoelastic problems of a hereditary-type constitutive law. The semidiscrete finite element method with linear time-stepping scheme and an elastic-viscoelastic correspondence principle are used in the theoretical development. An efficient treatment of hereditary integral is introduced to improve the numerical accuracy, to reduce the computation time, and to avoid the use of large memory storage. Two-dimensional numerical examples of plane strain and plane stress are solved under the assumption on the material property of being elastic in dilatation and like three-element Voigt model in distorsion, and compared with the analytical solutions and the past numerical results to show the versatility and efficiency of the proposed method.

• Journal of Industrial Technology
• /
• v.20 no.A
• /
• pp.203-209
• /
• 2000

The fluctuation integrals which give useful information in the structure of solution are associated with the mixed direct correlation integral ($C_{12}$) known. Using its weighted arithmetic mean of $C_{11}$ and $C_{22}$ and the activity coefficient model, the fluctuation integrals on solute-solute, solvent-solute, and solvent-solvent can be calculated in the function of mole fraction. In this work, several binary mixtures containing c-hexane were tested and the results on the fluctuation integrals were rather good.

• The Journal of the Acoustical Society of Korea
• /
• v.14 no.5
• /
• pp.58-62
• /
• 1995

Integral equations of self-radiation impedance of a rectangular piston in a rigid infinite baffle are derived using by polar coordinate. The self-radiation impedance is separated by two parts ; self-radiation resistance and self-radiation reactance. Derived integral equations are simulated by numerical method. Based on the numerical results, self-radiation impedance can be obtained in the low and high frequency ranges without any limited conditions.

• Korean Journal of Fisheries and Aquatic Sciences
• /
• v.27 no.6
• /
• pp.782-786
• /
• 1994

Integrated finite difference method (IFDM) is used to solve one dimensional free oscillation problem in the coastal area. To evaluate the solution accuracy of IFDM in free oscillation analysis, two finite difference equations based on area discretization method and point discretization method are derived from the governing equations of free oscillation, respectively. The difference equations are transformed into a generalized eigenvalue problem, respectively. A numerical example is presented, for which the analytical solution is available, for comparing IFDM to conventional finite difference equation (CFDM), qualitatively. The eigenvalue matrices are solved by sub-space iteration method. The numerical results of the two methods are in good agreement with analytical ones, however, IFDM yields better solution than CFDM in lower modes because IFDM only includes first order differential operator in finite difference equation by Green's theorem. From these results, it is concluded that IFDM is useful for the free oscillation analysis in the coastal area.

• Geophysics and Geophysical Exploration
• /
• v.3 no.1
• /
• pp.13-18
• /
• 2000

We describe the concept of the singularity in the integral equation of electromagnetic (EM) modeling in comparison with that in the integral representation of electric fields in EM theory, which would clarify the singular integral problems of the Green's tensor. We have also derived and classified the singular integrals of the Green's tensors in 3-D, 2.5-D and 2-D as well as in the thin sheet integral equations of the EM scattering problem, which have the most important effect on the accuracy of the numerical solution of the problems.