• 대한조선학회논문집
• /
• 제39권4호
• /
• pp.1-10
• /
• 2002

초기쇄파의 수치모사에는 지금까지 경계적분법이 주로 쓰여왔고, 이 방법은 과도한 계산시간의 문제를 제외하고는 어느 정도 성공적이라고 할 수 있다. 본 논문에서는 쇄파실험을 수치모사하기 위한 새로운 수치기법을 보였다. 이 수치기법은 고차 스펙트럴/경계요소법과 경계적분법을 순차적으로 사용하며, 계산시간을 현저히 줄여준다. 조파 및 파 에너지 집중과정은 고차 스펙트럴/경계요소법에 의해 효율적으로 수치모사되고, 파의 전복과정만이 경계적분법에 의해 계산된다. 계산예에서 높은 입자속도와 가속도 등 쇄파의 두드러진 특성이 보여졌다.

• 한국해양공학회:학술대회논문집
• /
• 한국해양공학회 2006년 창립20주년기념 정기학술대회 및 국제워크샵
• /
• pp.382-385
• /
• 2006

A high-order spectral/boundary-element method is newly adapted as an efficient numerical tool. In this method, the velocity potential is expressed as the sum of surface potential and body potential. Then, surface potential is solved fly using the high-order spectral method and body potential is solved fly using the high-order boundary element method. Through the combination of these two methods, the wave-making problems fly a submerged sphere moving with the large amplitude oscillation are solved in time-domain. With the example calculations, nonlinear effects on free-surface profiles and hydrodynamic forces are shown and discussed.

• 한국해양공학회지
• /
• 제13권1호통권31호
• /
• pp.113-120
• /
• 1999

In this paper, mathematical formulation of the high-order spectral/boundary-element method is shown. This method is one of the most efficient numerical methods by which the nonlinear gravity waves can be simulated in time-domain. Three-dimensional waves generated by a uniformly translating suriace pressure are calculated and discussed. The obtained results are compared with others results, The comparisons show good agreements.

• 한국전산구조공학회논문집
• /
• 제25권5호
• /
• pp.413-420
• /
• 2012

본 논문은 2차원 선형탄성 직접 경계요소법에서 저매개변수 요소를 사용할 때 Kernel의 적분방법에 대하여 논의하였다. 일반적으로 등매개변수 요소의 경우 형상함수로 통칭되는 해의 기저함수와 요소의 적분을 위해 사용되는 사상함수를 동일하게 사용한다. 그러나 본 논문에서는 사상함수의 차수를 낮게 취하여 순수기저절점을 도입하고 그때 직접 경계요소의 Kernel을 적분하기 위한 방법이 모색되었다. 일반적으로 경계요소법의 적분 Kernel의 경우 Log수치적분과 코쉬주치(Cauchy principal value) 등을 통해 해결하는데, 본 논문에서는 대수적 조작을 통해 적분값의 정확도를 높일 수 있도록 새로운 수식을 유도하였다. 본 연구에서 저매개변수 기반의 직접 경계요소에 대한 강건성과 정확도를 검증하기 위해 2차원 타원형 편미분방정식으로 표현되는 평면응력과 평면변형문제에 대해 적용하였다. 적용 예제로는 단순연결영역(simple connected region)의 대표적 문제인 캔틸레버보와 다중연결영역(multiple connected region)의 대표적인 문제인 개구부가 있는 사각평면에 대해 각각 수치해석을 수행한 결과 대폭적인 자유도의 감소에 비해 정확도 측면에는 기존의 방법과 차이가 없음을 볼 수 있었다. 본 논문에서 제시된 방법은 기저함수 고차화 저매개변수 직접 경계요소법(subparametric high order boundary element)과 이에 기초를 둔 저매개변수 고차 이중경계요소법(subparametric high order dual boundary element)의 초석이 될 수 있을 것이다.

• International Journal of Naval Architecture and Ocean Engineering
• /
• 제6권4호
• /
• pp.1000-1013
• /
• 2014

This paper presents some of the efforts by the authors towards numerical prediction of springing of ships. A time-domain Higher Order Boundary Element Method (HOBEM) based on cubic shape function is first presented to solve a complete second-order problem in terms of wave steepness and ship motions in a consistent manner. In order to avoid high order derivatives on the body surfaces, e.g. mj-terms, a new formulation of the Boundary Value Problem in a body-fixed coordinate system has been proposed instead of traditional formulation in inertial coordinate system. The local steady flow effects on the unsteady waves are taken into account. Double-body flow is used as the basis flow which is an appropriate approximation for ships with moderate forward speed. This numerical model was used to estimate the complete second order wave excitation of springing of a displacement ship at constant forward speeds.

• 동력기계공학회지
• /
• 제16권5호
• /
• pp.40-46
• /
• 2012

본 연구에서는 자유표면 유동문제를 효율적으로 계산하기 위한 방법으로서, 수치 파수조를 구현하여 잠수체에 의한 조파현상을 시간영역에서 다룰 수 있는 수치해법에 대하여 소개하였다. 그리고 이를 이용하여 수행된 연구내용들을 검토하고 양력물체의 경우를 포함하기 위하여 개선된 수학적 정식화 및 수치해법의 개요를 다루었다. 임의의 운동을 하는 양력물체에 의한 조파현상을 전산기로 구현하는 수치 파수조는 과중한 계산시간이 문제가 되는데, 이는 수치 Kutta 조건의 구현과 양력표면 후방의 wake 영역을 계산에서 고려해야 하기 때문이다. 따라서 한층 더 수치계산의 효율성이 중요하다고 판단되므로, 본 연구에서 소개된 3차원 고차 스펙트럴/경계요소법(High-Order Spectral/Boundary Element Method)은 자유표면 요소수를 N이라 할 때 그 계산량이 NlogN에 비례(N이 클 때는 거의 선형적으로 비례)하여 증가하므로 기존의 방법들 보다 매우 효율적인 수치해법이라 할 수 있다. 향후, 본 연구의 타당성 검증을 위한 수치코드의 개선과 여러 가지 수치계산결과 비교 등의 노력이 더 필요하다고 생각된다.

• 한국해양공학회지
• /
• 제20권6호
• /
• pp.67-74
• /
• 2006

A high-order spectral/boundary-element method is newly adapted as an efficient numerical tool. This method is one of the most efficient numerical methods by which the nonlinear gravity waves can be simulated and hydrodynamic forces also can be calculated in time-domain. In this method, the velocity potential is expressed as the sum of surface potential and body potential. Then, surface potential is solved by using the high-order spectral method and body potential is solved by using the high-order boundary element method. By the combination of these two methods, the wave-making problems by a submerged sphere oscillating with large amplitude under the free~surface are solved in time-domain. Through the example calculations, nonlinear effects on free-surface profiles and hydrodynamic forces are shown and discussed.

• 동력기계공학회지
• /
• 제14권6호
• /
• pp.75-82
• /
• 2010

In this study, the topics for free-surface wave simulation, nonlinear hydrodynamic force, and the critical resonance frequency of so-called ${\tau}=U{\omega}/g$=1/4 are discussed. A high-order spectral/boundary element method is newly adapted as an efficient numerical tool. This method is one of the most efficient numerical methods by which the nonlinear gravity waves can be simulated and hydrodynamic forces also can be calculated in time domain. In this method, the velocity potential is expressed as the sum of surface potential and body potential. Then, surface potential is solved by using the high-order spectral method and body potential is solved by using the high-order boundary element method. By the combination of these two methods, the wave-making problems by a submerged sphere oscillating with forward speed under the free-surface are solved in time domain.

• Coupled systems mechanics
• /
• 제12권4호
• /
• pp.365-389
• /
• 2023

Dynamic analysis of a typical concrete gravity dam-reservoir system is formulated by FE-(FE-TE) approach (i.e., Finite Element-(Finite Element-Truncation Element)). In this technique, dam and reservoir are discretized by plane solid and fluid finite elements. Moreover, the H-W (i.e., Hagstrom-Warburton) high-order condition imposed at the reservoir truncation boundary. This task is formulated by employing a truncation element at that boundary. It is emphasized that reservoir far-field is excluded from the discretized model. The formulation is initially reviewed which was originally proposed in a previous study. Thereafter, the response of Pine Flat dam-reservoir system is studied due to horizontal and vertical ground motions for two types of reservoir bottom conditions of full reflective and absorptive. It should be emphasized that study is carried out under high order of H-W condition applied on the truncation boundary. The initial part of study is focused on the time harmonic analysis. In this part, it is possible to compare the transfer functions against corresponding responses obtained by FE-(FE-HE) approach (referred to as exact method). Subsequently, the transient analysis is carried out. In that part, it is only possible to compare the results for low and high normalized reservoir length cases. Therefore, the sensitivity of results is controlled due to normalized reservoir length values.

• Structural Engineering and Mechanics
• /
• 제73권1호
• /
• pp.97-108
• /
• 2020

This study presents a 4 node, 11 DOF/node plate element based on higher order shear deformation theory for lamina composite plates. The theory accounts for parabolic distribution of the transverse shear strain through the thickness of the plate. Differential field equations of composite plates are obtained from energy methods using virtual work principle. Differential field equations of composite plates are obtained from energy methods using virtual work principle. These equations were transformed into the operator form and then transformed into functions with geometric and dynamic boundary conditions with the help of the Gâteaux differential method, after determining that they provide the potential condition. Boundary conditions were determined by performing variational operations. By using the mixed finite element method, plate element named HOPLT44 was developed. After coding in FORTRAN computer program, finite element matrices were transformed into system matrices and various analyzes were performed. The current results are verified with those results obtained in the previous work and the new results are presented in tables and graphs.