• Journal of Ship and Ocean Technology
• /
• 제10권1호
• /
• pp.10-26
• /
• 2006

The radiation problem for oscillating bodies on the free surface has been formulated by the over-determined Green integral equation, where the boundary condition on the free surface is satisfied by adopting the Kelvin-type Green function and the irregular frequencies are removed by placing additional control points on the free surface surrounded by the body. The B-Spline based higher order panel method is then applied to solve the problem numerically. Because both the body geometry and the potential on the body surface are represented by the B-Splines, that is in polynomials of space parameters, the unknown potential can be determined accurately to the order desired above the constant value. In addition, the potential expressed in B-Spline can be differentiated analytically to get the velocity on the surface without introducing any numerical error. Sample computations are performed for a semispherical body and a rectangular box floating on the free surface for six-degrees of freedom motions. The added mass and damping coefficients are compared with those by the already-validated constant panel method of the same formulation showing strikingly good agreements.

• 대한조선학회지
• /
• 제22권3호
• /
• pp.9-18
• /
• 1985

The numerical calculation for solving boundary-value problem related to potential flows with a free surface is carried out by application of the localized finite element method. Only forced motion of 2-D body in infinitely deep fluid is considered, although this schemes is equally applicable to any first order time-harmonic problems of similar nature. The infinite domain of the fluid is separated into the inner flow field and the outer flow field with common inter-surface boundary. The finite element method is applied to obtain the solution in the inner flow field and the Green functions are utilized to represent the solution in the outer flow field. At the inter-surface boundary, the continuity of the value of potential and the normal derivative of the potential(i.e. matching condition) is conserved. The present method has better computational efficiency than the previous LFEM and the integral equation method of Frank. This enhanced computational efficiency is presumably due to the fact that the present method gives a symmetric coefficient matrix and requires less computational time in calculating the influence coefficient matrix of Green function than the integral equation method. And the irregular frequency desen't exist because the uniqueness of the solution is assured by the such that the exact free surface condition is satisfied on the boundary of the localized finite element region(i.e. inner region). As an example of the above method, the hydrodynamic forces for the circular cylinder and the rectangular cylinders are calculated. In the computed results, the small number of singularity distribution segments($3{\sim}6$) give good result relative to Ursell's and Vugts'.

• 대한수학회논문집
• /
• 제16권1호
• /
• pp.137-152
• /
• 2001

Steady two-dimensional flows due to an applied pressure distribution in water of finite depth are considered. Gravity is included in the dynamic boundary condition. Gravity is included in the dynamic boundary condition. The problem is solved numerically by using the boundary integral equation technique. It is shown that, for both supercritical and subcritical flows, solutions depend on three parameters: (i) the Froude number, (ii) the magnitude of applied pressure distribution, and (iii) the span length of pressure distribution. For supercritical flows, there exist up to two solutions corresponding to the same value of Froude number for positive pressures and a unique solution for negative pressures. For subcritical flows, there are solutions with waves behind the applied pressure distribution. As the Froude number decreases, these waves when the Froude numbers approach the critical values.

• 대한기계학회논문집
• /
• 제16권5호
• /
• pp.838-848
• /
• 1992

본 연구에서는 김성호등의 형상모델인 Fig.1(b)에서 전단하중이 작용하는 경 우에 대한, 복합재료의 탄성층 내부(레진층)에 존재하는 중앙균열의 응력확대계수 산 출을 위하여 균열부위를 제외하고는 섬유층과 레진층이 완전히 접착되었다고 가정한 모델을 다음과 같이 설정 하였다. 접착레진을 주로하는 탄성층(resin rich layer)을 중심으로 상하 각1개의 섬유층(fiber rich layer)과 균질한 특성을 갖는 복합재료의 층으로 단순화하였으며, 복합재료는 레진층이나 섬유층에 비하여 무한히 두꺼우므로 반무한체로 이상화 하였다. 선형탄 이론에 의한 혼합경계조건문제(mixed boundary value problem)로 부터 제2종 Fredholm 적분방정식(fredholm integral equation of a second kind)을 유도하였으며 수치해석적인 방법에 의하여 응력확대계수를 구하였다. 또한, 복합재료의 재료물성 및 균열길이, 섬유두께등이 기하학적 변수에 대하여 응력 확대계수를 산출하였다.

• Structural Engineering and Mechanics
• /
• 제52권1호
• /
• pp.1-17
• /
• 2014

An interfacial penny-shaped crack between piezoelectric layer and elastic half-space subjected to mechanical and electric loads is investigated. Using Hankel transform technique, the mixed boundary value problem is reduced to a system of singular integral equations. The integral equations are further reduced to a system of algebraic equations with the aid of Jacobi polynomials. The stress intensity factor and energy release rate are determined. Numerical results reveal the effects of electric loadings and material parameters of composite on crack propagation and growth. The results seem useful for design of the piezoelectric composite structures and devices of high performance.

• Structural Engineering and Mechanics
• /
• 제23권2호
• /
• pp.163-175
• /
• 2006

In this paper, the dynamic response of a piezoelectric layer with a penny-shaped crack is investigated. The piezoelectric layer is subjected to an axisymmetrical action of both mechanical and electrical impacts. Two kinds of crack surface conditions, i.e., electrically impermeable and electrically permeable, are adopted. Based upon integral transform technique, the crack boundary value problem is reduced to a system of Fredholm integral equations in the Laplace transform domain. By making use of numerical Laplace inversion the time-dependent dynamic stress and electric displacement intensity factors are obtained, and the dynamic energy release rate is further derived. Numerical results are plotted to show the effects of both the piezoelectric layer thickness and the electrical impact loadings on the dynamic fracture behaviors of the crack tips.

• 대한조선학회논문집
• /
• 제32권3호
• /
• pp.62-71
• /
• 1995

비압축성 가정하에 Navier-Stokes 방정식을 이용하여 비정상 점성유동을 수치해석하기 위해서는 매시간 단계에서 타원형 압력 Poisson 방정식의 해를 구해야 하며, 이에 많은 계산시간이 소요된다. 본 논문에서는 직접법을 이용하여 압력 Poisson 방정식을 수치해석하였으며, 분할수 증가에 따른 소요시간 문제를 다루었다. Green 정리를 압력 Poisson 방정식에 적용하면 주어진 문제는 경계치문제로 변환되고, convolution 형의 영역적분은 F.F.T.를 이용하여 계산시간을 단축할 수 있어, 직접법 이용시 소요시간은 경계치문제의 해를 구하는 데에 좌우된다. 직접법의 검증을 위하여 해석해를 알고 있는 경우에 대하여 수치해석하였고, 물체경계조건과 정합문제에 관하여 수치해석 하였는데. 분할수가 (n.n) 시 O($n^{3}$) 미만의 계산시간으로 수치해석할 수 있었다.

• 한국해양환경ㆍ에너지학회지
• /
• 제13권2호
• /
• pp.83-90
• /
• 2010

해양파 중에서 2개 선체가 가까운 거리에 계류되어있는 경우에 2선체의 운동을 엄밀한 파 중 2 부유체 운동이론에 의하여 계산하였다. 2 선체의 방사 및 산란파 포텐셜은 특이파 수 현상이 없는 주파수영역 Improved Green 적분방 정식으로부터 구하였다. 규칙파 중에서 근접한 2선체가 서로 연결되지 않은 경우와 강체 연결장치로 병렬 계류된 된 경우의 운동 및 강체연결체에 작용하는 6방향 힘과 모오멘트를 해석하는 기법과 수치실험 결과도 보였다.

• 한국지진공학회:학술대회논문집
• /
• 한국지진공학회 1997년도 춘계 학술발표회 논문집 Proceedings of EESK Conference-Fall 1997
• /
• pp.128-135
• /
• 1997

The indirect boundary integral equation is formulated to analyze the behavior of a cavity in a multilayered half-plane subjected to earthquake waves. This formulation uses the fundamental solutions that are numerically calculated by the generalized transmission and reflection coefficient method. The free surface of the cavity without external excitation influences the behavior of the half-plane. Consequently this analysis adds the consideration of scattering-field into the analysis and the total motion field of the cavity is decomposed into the free-field and scattering-field motions. The free-field motion is obtained from the modification of the transmission and reflection coefficient method. The scattering-field motion is calculated is calculated by the indirect boundary value problem which has the ficticious boundaries and sources. In this study, P wave, SV wave, SH wave, and Rayleigh wave are analyzed respectively.

• 대한조선학회논문집
• /
• 제28권1호
• /
• pp.27-37
• /
• 1991

본 논문에서는 축대칭 회전 몰수체가 자유표면하에서 운동하고 있을때 이에 의한 비선형 표면파의 생성과 물체에 작용하는 힘을 계산하였다. 축대칭 포텐시얼 경계치 문제를 해석하기 위하여 경계 적분 방정식을 풀었으며, 시간에 따른 자유표면의 변화를 추적하기 위한 수치적분방법으로 Runge-Kutta 4차 방법을 사용하였다. 이 결과로부터 축대칭 몰수체에 작용하는 힘을 계산하였고, 선형이론과 비선형 이론에 의한 결과를 비교하였다.