• Korean Journal of Fisheries and Aquatic Sciences
• /
• v.31 no.2
• /
• pp.202-209
• /
• 1998

A three-dimensional numerical method based on the Green's integral equation is developed to predict the motion response and drift force in waves for the ocean monitoring facilities. In this method, we use source and doublet distribution, and triangular and rectangular eliments. To eliminate the irregular frequency phenomenon, the method of improved integral equation is applied and the time-mean drift force is calculated by the method of direct pressure integration over the body surface. To conform the validity of the present numerical method, some calculations for the floating sphere are performed and it is shown that the present method provides sufficiently reliable results. As a calculation example for the real facilities, the motion response and the drift force of the vertical cylinder type ocean monitoring buoy with 2.6 m diameter and 3,77 m draft are calculated and discussed. The obtained results of motion response can be used to determine the shape and dimension of the buoy to reduce the motion response, and other data such as the effect of motion reduction due to a damper can be predictable through these motion calculations. Also, the calculation results of drift force can be used in the design procedure of mooring system to predict the maximum wave load acting on the mooring system. The present method has, in principle, no restriction in the application to the arbitrary shape facilities. So, this method can be a robust tool for the design, installation, and operation of various kinds of the floating-type ocean monitoring facilities.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2008.05a
• /
• pp.1634-1638
• /
• 2008

하천의 2차원 흐름 해석, 유사이동 해석, 오염확산 해석을 위한 유체의 수치해석법에는 유한요소법, 유한차분법, 유한차분법의 변형인 유한체적법, 경계적분법 등이 있다. 유체에 대한 수치해석 기법으로 전통적으로 가장 많이 사용되고 있는 방법은 유한차분법이지만, 비구조적 요소망(unstructured mesh)을 이용하여 복잡한 형상을 표현하기가 상대적으로 용이한 유한요소법이 다양한 형태의 하천 해석에는 더욱 적합할 것이다. 본 연구에서는 비구조적 요소망을 advanced front method를 이용하여 생성해 보았다. Advanced front method는 해석하고자 하는 영역에 적절한 절점들을 생성한 후 삼각 요소망을 구성하는 grid based advanced front method와 절점들을 생성하지 않고 해석 영역에 삼각 요소를 바로 구성하는 element based advanced front method로 나눌 수 있다. Grid based advanced front method에서 해석 영역에 적절한 절점을 생성하는 방법으로는 일반적인 격자 구조의 절점 생성 방법을 적용하였으며 경계와의 거리가 가까운 절점은 생성되지 않으며, 삼각 요소를 구성할 때에는 두 개의 인접 절점을 비교하여 최적의 삼각 요소를 구성하게 된다. 단 두 개의 인접 절점만을 비교함으로서 비교적 빠른 시간 안에 최적의 삼각 요소망을 구성할 수 있다. 삼각 요소망을 생성한 후에는 Laplacian smoothing을 이용하여 삼각 요소망의 형질을 개선하였다. Element based advanced front method는 외부 경계에서부터 시작된 Front가 내부 영역으로 확대되어지며 각 Front에서 적절한 절점을 직접 생성하여 바로 삼각 요소를 구성하게 된다. Front에서 생성된 절점은 인접 절점들이 있는지 검색하여 인접 절점이 있다면 생성된 절점은 삭제되어지며 인접 절점이 삼각 요소를 위한 나머지 한 점으로 채택되어진다. Front는 외부 경계와 교차되어지지 않아야 하며 또한 연속된 Front를 효율적으로 관리하기 위해 list 자료 구조를 활용하였다.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.9
• /
• pp.1711-1721
• /
• 1992

A nonlinear dissipative dynamical system can often have multiple attractors. In this case, it is important to study the global behavior of the system by determining the global domain of attraction of each attractor. In this paper we study the global behavior of a forced beam with two mode interaction. The governing equation of motion is reduced to two second-order nonlinear nonautonomous ordinary differential equations. When .omega. /=3.omega.$_{1}$ and .ohm.=.omega $_{1}$, the system can have two asymptotically stable steady-state periodic solutions, where .omega./ sub 1/, .omega.$_{2}$ and .ohm. denote natural frequencies of the first and second modes and the excitation frequency, respectively. Both solutions have the same period as the excitation period. Therefore each of them shows up as a period-1 solution in Poincare map. We show how interpolated mapping method can be used to determine the two four-dimensional domains of attraction of the two solutions in a very effective way. The results are compared with the ones obtained by direct numerical integration.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.27 no.1
• /
• pp.63-72
• /
• 1990

The linear hydrodynamic forces, acting on a forced oscillating cylinder from its mean position on a free surface with a small amplitude, are calculated in the time domain. The integral equation method using a time dependent Green function is employed. The numerical results for the heaving and swaying circular cylinder are shown and give good agreements with others Furthermore it is shown that the use of the Green function, which is expressed by a series expansion or asymptotic expansion according to time range, reduces computing time greatly.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.20 no.6
• /
• pp.789-799
• /
• 2007

In this paper, we propose a stochastic finite element formulation which takes into account the randonmess in the material and geometrical parameters. The formulation is proposed for plate structures, and is based on the weighted integral approach. Contrary to the case of elastic modulus, plate thickness contributes to the stiffness as a third-order function. Furthermore, Poisson's ratio is even more complex since this parameter appears in the constitutive relations in the fraction form. Accordingly, we employ Taylor's expansion to derive decomposed stochastic field functions in ascending order. In order to verify the proposed formulation, the results obtained using the proposed scheme are compared with those in the literature and those of Monte Carlo analysis as well.

• Journal of Korean Society of Steel Construction
• /
• v.15 no.6 s.67
• /
• pp.693-700
• /
• 2003

An evaluation method for flutter derivatives is proposed, using indicial functions of structural members produced by Computational Fluid Dynamics (CFD). Flutter derivatives are obtained by Fourier integration of indicial functions. Instead of direct simulation of oscillating objects, only the calculation of time-dependent lift and moment variations of fixed objects with constant attack angle are necessary.The Finite Element Method (FEM) is developed as a tool for the numerical method. For two rectangular sections having different aspect ratios, the numerical analysis and wind tunnel test are carried out to inspect the adequacy of this study. The results proved to be good, and they could be used for a preliminary design.

• Journal of the Korean Geotechnical Society
• /
• v.35 no.2
• /
• pp.19-27
• /
• 2019

In this study, an analytical model is proposed to estimate the buckling responses of tapered piles. The governing differential equation of the soil-pile system considering the tapering and side friction of the pile and the soil nonhomogeneity is derived, which is numerically integrated by the Runge-Kutta method and then the eigenvalue of bucking load is determined by Regula-Falsi algorithm. For a cylindrical pile, the results obtained from this study are found to compare well with those reported in literature. Illustrative examples for buckling load and stress as well as buckled shape are provided to investigate the effects of dimensionless parameters related to the soil-pile system.

• Journal of the Korean Geosynthetics Society
• /
• v.18 no.3
• /
• pp.69-77
• /
• 2019

This paper presents a new analytical model for estimating the free vibration of tapered friction piles. The governing differential equation for the free vibration of statically axially-loaded piles embedded in non-homogeneous soil is derived. The equation is numerically integrated by the Runge-Kutta method, and then the eigenvalue of natural frequency is determined by the Regula-Falsi scheme. For a cylindrical non-tapered pile, the computed natural frequencies compare well with the available data from literature. Numerical examples are presented to investigate the effects of the tapering, the skin friction resistance, the end condition of the pile, the vertical compressive loading, and the soil non-homogeneity on the natural frequency and mode shape of tapered friction piles.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.18 no.3
• /
• pp.311-320
• /
• 2005

In this paper an enhanced quadratic finite element for static and dynamic analysis of plate structures is presented. The performance of a proposed plate element is improved by the coupled use of non conforming displacement modes, the selective integration scheme, and the assumed shear strain fields. An efficient direct modification method is also applied to this element to solve the problem such as failure of the patch test due to the adoption of non conforming modes. The proposed quadratic finite element does not show any spurious mechanism and does not produce shear locking phenomena even with distorted meshes. It is shown that the results obtained by this element converged to analytical solutions very rapidly tough numerical tests for standard benchmark problems. It is also noted that this element is applicable to transient dynamic analysis of Mindlin plates.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.23 no.5
• /
• pp.525-533
• /
• 2010

Direct time integration method such as Newmark method is numerically performed under the assumption that continuous load function such as constant amplitude load can be treated as a discontinuous load fuction. It is due that the load can be treated as a constant value at the given time period regardless of variation of load at the time increment interval. It means the numerical results should be accompanied by the error due to approximation of load fuction. In contrast, the load function is calculated by convolution integral for the given time interval at finite element equation based on Gurtin's variation equation. Therefore. precise numerical results can be obtained by Gurtin's method because of convolution integral for the continuous load fuction curve even at the variation of load function in the given time interval. In this study, we prove that Gurtin's method can be more suitable than Newmark method in the problem of constant amplitude loading, using the numerical results for the free end of the one-dimensional rod. This study also shows that Gurtin's method is more effective in constant amplitude loading than in constant loading. The accuracy and the validity are verified by comparison between the results of in-house FORTRAN code and ADINA, a commercial software supporting Newmark method.