• 한국해양공학회지
• /
• 제32권6호
• /
• pp.447-457
• /
• 2018

In this study, a two-dimensional fully nonlinear transient wave numerical tank was developed using a desingularized indirect boundary integral equation method. The desingularized indirect boundary integral equation method is simpler and faster than the conventional boundary element method because special treatment is not required to compute the boundary integral. Numerical simulations were carried out in the time domain using the fourth order Runge-Kutta method. A mixed Eulerian-Lagrangian approach was adapted to reconstruct the free surface at each time step. A numerical damping zone was used to minimize the reflective wave in the downstream region. The interpolating method of a Gaussian radial basis function-type artificial neural network was used to calculate the gradient of the free surface elevation without element connectivity. The desingularized indirect boundary integral equation using an isolated point source and radial basis function has no need for information about the element connectivity and is a meshless method that is numerically more flexible. In order to validate the accuracy of the numerical wave tank based on the desingularized indirect boundary integral equation method and meshless technique, several numerical simulations were carried out. First, a comparison with numerical results according to the type of desingularized source was carried out and confirmed that continuous line sources can be replaced by simply isolated sources. In addition, a propagation simulation of a $2^{nd}$-order Stokes wave was carried out and compared with an analytical solution. Finally, simulations of propagating waves in shallow water and propagating waves over a submerged bar were also carried and compared with published data.

• 한국해양환경ㆍ에너지학회지
• /
• 제6권4호
• /
• pp.45-51
• /
• 2003

본 연구에서는 경계요소법을 이용하여 불투과성 사다리꼴형상 잠제에 의한 파랑의 반사율과 투과율을 이론적으로 수치해석 하였다. 해석기법으로는 유체와 투과성 영역을 동시에 해석할 수 있는 파압함수를 사용하였으며, 이론적으로 간단히 하기 위해, 투과성 영역내의 파랑의 운동은 선형소산계수와 부가질량계수를 도입하여 정식화하였다. 사다리꼴형상 잠제의 소파특성은 잠제의 법면경사의 변화에 대한 의존도가 높은 것을 알 수 있었다.

• Earthquakes and Structures
• /
• 제9권4호
• /
• pp.753-766
• /
• 2015

An investigation has been carried out for the propagation of torsional surface waves in an inhomogeneous prestressed layer over an inhomogeneous half space when the upper boundary plane is assumed to be rigid. The inhomogeneity in density, initial stress (tensile and compressional) and rigidity are taken as an arbitrary function of depth, where as for the elastic half space, the inhomogeneity in density and rigidity is hyperbolic function of depth. In the absence of heterogeneities of medium, the results obtained are in agreement with the same results obtained by other relevant researchers. Numerically, it is observed that the velocity of torsional wave changes remarkably with the presence of inhomogeneity parameter of the layer. Curves are compared with the corresponding curve of standard classical elastic case. The results may be useful to understand the nature of seismic wave propagation in geophysical applications.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2001년도 추계학술대회논문집A
• /
• pp.341-348
• /
• 2001

A Mixed Volume and Boundary Integral Equation Method is applied for the effective analysis of plane elastostatic problems in unbounded solids containing general anisotropic inclusions and voids or isotropic inclusions. It should be noted that this newly developed numerical method does not require the Green's function for anisotropic inclusions to solve this class of problems since only Green's function for the unbounded isotropic matrix is involved in their formulation for the analysis. This new method can also be applied to general two-dimensional elastodynamic and elastostatic problems with arbitrary shapes and number of anisotropic inclusions and voids or isotropic inclusions. Through the analysis of plane elastostatic problems in unbounded isotropic matrix with orthotropic inclusions and voids or isotropic inclusions, it will be established that this new method is very accurate and effective for solving plane elastic problems in unbounded solids containing general anisotropic inclusions and voids or isotropic inclusions.

• 한국소음진동공학회논문집
• /
• 제26권5호
• /
• pp.602-608
• /
• 2016

The paper proposes a practical meshless method for the free vibration analysis of clamped plates having arbitrary shapes by extending the non-dimensional dynamic influence function (NDIF) method, which was developed by the author in 1999. In the proposed method, the domain and boundary of the plate of interest are discretized using only nodes without elements unlike FEM and the system matrices are obtained by making domain nodes and boundary nodes satisfy the governing differential equation and boundary conditions, respectively. However, since the above system matrices are not square ones, the problem of free vibrations of clamped plates is not reduced to an algebraic eigenvalue problem. An additional theoretical treatment is considered to produce an algebraic eigenvalue problem. It is revealed from case studies that the proposed method is valid and accurate.

• 전자공학회논문지SC
• /
• 제41권2호
• /
• pp.45-52
• /
• 2004

본 논문에서는 슬라이딩 라인 함수의 각 요소를 간편 간접추론 퍼지 논리의 입력변수로 사용하여 슬라이딩 모드 제어기의 경계층을 모사하는 간편 간접추론 퍼지논리 경계층을 갖는 슬라이딩 모드 제어기를 설계하였다. 제안된 제어시스템은 별도의 비선형 함수를 구할 필요가 없으며, 네 개의 퍼지규칙으로 이루어지므로 간단하고 안정도 증명이 쉬운 장점을 가진다. 본 논문의 유효성을 비선형 시변 시스템에 적용하여 고찰하였다.

• 수산해양기술연구
• /
• 제27권2호
• /
• pp.120-124
• /
• 1991

이상의 결과를 정리하면 다음과 같이 요약할 수 있다. (1)Karman의 적분 방정식에 미소 거칠기 영향을 고려함으로서 디퓨져 닥트 표면의 경계층 계산에 응용한 결과 모멘트 적분법 및 실험치와의 비교에서 개선된 결과임이 확인되었다. (2) 국소 거칠기의 효과를 주는 방법으로는 Cole의 벽 및 와 법칙에 Clauser의 거칠기 함수와, 압력 기울기를 고려한 부가 형상계수 값으로 적분 방정식에 응용할 수 있다. (3) 국소 거칠기 분포에 의하여 경계층 특성을 교란시켜 표면 마찰력 계수를 줄일 수 있어 마찰력 손실을 줄일 수 있는 방안이 제시되었다.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 1997년도 봄 학술발표회 논문집
• /
• pp.3-10
• /
• 1997

Meshless particle method is a numerical technique which does not use the concept of element. This method can easily handle special engineering problems which cause difficulty in the use of finite element method, however it has a drawback that essential boundary condition is not satisfied. In this paper, several studies for satisfying essential boundary conditions and enhancing the accuracy of solutions are discussed. Particular emphasis is placed on a new numerical technique in which finite elements are used on the boundaries to satisfy the essential boundary conditions and meshless particle method is used in the interior domain. For coupling of the two methods interface elements are introduced into the zone between the subdomains using meshless particle method and finite element method. The shape functions and the approximated displacement functions of the interface element are derived with the ramp function based on the shape function of finite elements. The whole numerical procedures are formulated by Galerkin method. Several numerical examples for enhancing the accuracy of solution in the meshless particle method and a new coupling method are presented.

• 한국안전학회지
• /
• 제24권1호
• /
• pp.37-42
• /
• 2009

The present study is to estimate the effect of wave height affecting at the front face of breakwater systems, when a submarine trench is dredged in the distant offshore from outer breakwater. The wave diffraction field, which is important hydraulic factor in the ocean, is considered to be two dimensional(2D) plane and the configuration of the submarine trench on the sea bed designated by single horizontal long-rectangular system. The numerical simulation is performed by using Green function based on the boundary integral equation and meshed at moving boundary conditions. The results of present numerical simulations are illustrated by applying the normal incidence. It is shown that the ratios of wave height reduction at the front face of breakwater systems are approximately 20% by the effect of placing long trench on the sea bed. This study can effectively be utilized for safety assessment to various breakwater systems in the ocean field.

• 대한기계학회논문집B
• /
• 제30권3호
• /
• pp.270-277
• /
• 2006

Inverse radiation problems were solved for estimating boundary temperature distribution in a way of function estimation approach in an axisymmetric absorbing, emitting and scattering medium, given the measured radiative data. In order to reduce the computational time fur the calculation of sensitivity matrix, automatic differentiation and Broyden combined method were adopted, and their computational precision and efficiency were compared with the result obtained by finite difference approximation.. In inverse analysis, the effects of the precision of sensitivity matrix, the number of measurement points and measurement error on the estimation accuracy had been inspected using quasi-Newton method as an inverse method. Inverse solutions were validated with the result acquired by additional inverse methods of conjugate-gradient method or Levenberg-Marquardt method.