• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.30 no.4
• /
• pp.335-341
• /
• 2017

In this paper, beam finite elements with higher-order derivatives' continuity are formulated and evaluated for various boundary conditions. All the beam elements are based on Euler-Bernoulli beam theory. These higher-order beam elements are often required to analyze structures by using newly developed higher-order beam theories and/or non-classical beam theories based on nonlocal elasticity. It is however rare to assess the performance of such elements in terms of boundary and loading conditions. To this end, two higher-order beam elements are formulated, in which $C^2$ and $C^3$ continuities of the deflection are enforced, respectively. Three different boundary conditions are then applied to solve beam structures, such as cantilever, simply-support and clamped-hinge conditions. In addition to conventional Euler-Bernoulli beam boundary conditions, the effect of higher-order boundary conditions is investigated. Depending on the boundary conditions, the oscillatory behavior of deflections is observed. Especially the geometric boundary conditions are problematic, which trigger unstable solutions when higher-order deflections are prescribed. It is expected that the results obtained herein serve as a guideline for higher-order derivatives' continuous finite elements.

• Journal of the Korea Institute of Military Science and Technology
• /
• v.1 no.1
• /
• pp.128-144
• /
• 1998

The BEM for linear elastic stress analysis is applied to the highly rotating axisymmetric body problem which also involves the thermoelastic effects due to steady-state thermal conduction. The axisymmetric BEM formulation is briefly summarized and an alternative approach for transforming the volume integrals associated with such body force kernels into equivalent boundary integrals is described in a way of using the concept of inner product and vector identity. A discretization scheme for higher order BE is outlined for numerical treatment of the resulting boundary integral equations, and it is consequently illustrated by determining the stress distributions of the turbine rotor disk of the small turbojet engine(ADD 500) for which a FEM stress solution has been furnished by author.

• Journal of the Society of Naval Architects of Korea
• /
• v.39 no.4
• /
• pp.1-10
• /
• 2002

For the time-domain simulation of incipient breaking waves, usually the boundary integral method has been used so far, and it seems to be successful except a problem of too much computation time. The present paper shows a new computation technique for the simulation of breaking wave experiment. This technique uses the high-order spectral/boundary element method and the boundary integral method in sequence, and reduces the computation time remarkably. The wave generation and energy focusing process is efficiently simulated by the high-order spectral/boundary element method. Only the wave over-turning process is simulated by the boundary integral method. In the example calculation result, salient features of breaking waves such as high particle velocities and accelerations are shown.

• Journal of the Society of Naval Architects of Korea
• /
• v.38 no.3
• /
• pp.31-40
• /
• 2001

Low-order panel methods are being used to design marine propellers. Since the potential value over each panel for these methods is assumed to be a constant, the accuracy of prediction is known to be limited. Therefore, a higher order boundary element method(HOBEM) has been studied to enhance the accuracy of prediction. In this paper, a HOBEM representing the body boundary surfaces and physical quantities by a 9-node Lagrangian shape function is employed to analyse the flow around marine propellers in steady potential flow. First, the numerical results for a circular wing with thickness variations are compared with Jordan's linear solution. Then, the computational results of two propellers(DTRC 4119 & DTRC 4842 propeller) are compared with the experimental and numerical results published. The pressure distribution on the surface of the propeller is also compared with experimental data.

• The Journal of Korean Institute of Communications and Information Sciences
• /
• v.17 no.2
• /
• pp.181-186
• /
• 1992

The scattering property of TMz illuminated perfectly conducting and dielectric cylinders of arbitrary cross section are analyzed by the boundary element techniques. The boundary element equations are formulated via Maxwell’s equations, weighted residual or Green’s theorem, and the boundary conditions. The unknown surface fields on the boundaries are then calculated by the boundary element integral equations. Once the surface fields are found, the scattered fields in from a perfectly conducting circular and elliptic cylinders, a dielectric circular and elliptic cylinders are numerically analyzed. A general computer program has been developed using the quadratic elements(higher order boundary elements) and the Gaussian quadrature.

• Journal of Ocean Engineering and Technology
• /
• v.12 no.1
• /
• pp.99-106
• /
• 1998

3차원 파수조에서 완전 비선형파를 시뮬레이션하기 위하여 우선 랜킨 소스를 기저로한 적분방정식을 고차경계요소법을 이용하여 이산화하였다. 그리고 방사경계조건은 파흡수 비치와 포텐셜 스트레칭 기법을 이용하여 모델링하였으며, 비선형 자유표면과 경계조건식은 고차 예측 및 보정 기법을 이용하여 시간 적분하였다. 파흡수 비치는 파의 진행방향 특성에 따라 수조내에 다양하게 배치할 수 있으며 비칭서 흡수가 덜된 파는 수조의 길이 방향 끝단에서 포텐셜 스트레칭 기법에 의하여 반사없이 진행하도록 하였다. 수치실험 결과 일-에너지 보존법칙과 모멘텀-임펄스 보존 법칙이 만족됨으로써 본 수치기법의 효용성이 검증되었다.

• Journal of Ocean Engineering and Technology
• /
• v.13 no.4 s.35
• /
• pp.1-9
• /
• 1999

일차원 및 이차원의 단순한 문제에 대하여 계층 요소를 사용한 혼합 차수 유한 요소해의 정확성 및 수렴성을 조사하였다. 이러한 작업은 임의의 차수를 가진 블록들을 조합하여 요소를 구성함으로서 이루어질 수 있다. 블록간의 연결성이 유지될 수 있는 블록의 구성과 요소 생성에 대하여는 코드개발과 관련하여 설명하고 있으며, 서로 다른 차수를 가진 인접 블록간의 해의 연속성에 대하여는 계층 요소의 구성과 관련하여 서술되었다. 수치적 결과는 블록의 차수를 잘 선택함으로서 유한 요소해의 수렴성과 정확성을 증가시킬 수 있음을 보여주고 있으며, 고차 요소 영역을 너무 많이 할당하여 선형 요소의 영역이 너무 적을 경우에는 경계 조건에 따라 오차가 내부로 전파됨을 보여준다. 또한 세분화된 요소에 대한 고차 보간의 경우, 해의 수렴성이 저해될 수 있음이 발견되었다.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.25 no.5
• /
• pp.413-420
• /
• 2012

In this study, the Kernel integration scheme for 2D linear elastic direct boundary element method has been discussed on the basis of subparametric element. Usually, the isoparametric based boundary element uses same polynomial order in the both basis function and mapping function. On the other hand, the order of mapping function is lower than the order of basis function to define displacement field when the subparametric concept is used. While the logarithmic numerical integration is generally used to calculate Kernel integration as well as Cauchy principal value approach, new formulation has been derived to improve the accuracy of numerical solution by algebraic modification. The subparametric based direct boundary element has been applied to 2D elliptical partial differential equation, especially for plane stress/strain problems, to demonstrate whether the proposed algebraic expression for integration of singular Kernel function is robust and accurate. The problems including cantilever beam and square plate with a cutout have been tested since those are typical examples of simple connected and multi connected region cases. It is noted that the number of DOFs has been drastically reduced to keep same degree of accuracy in comparison with the conventional isoparametric based BEM. It is expected that the subparametric based BEM associated with singular Kernel function integration scheme may be extended to not only subparametric high order boundary element but also subparametric high order dual boundary element.

• Journal of Ocean Engineering and Technology
• /
• v.13 no.1 s.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.

• Journal of the Society of Naval Architects of Korea
• /
• v.39 no.1
• /
• pp.8-15
• /
• 2002

A two-dimensional higher order panel method using B-splines has been developed to overcome the disadvantages of the low order panel method and to obtain more accurate solution. The sources and the normal dipoles are distributed on both the body and the free surface. Instead of applying the upwind finite difference schemes to satisfy the linearized free surface and the radiation condition, the derivatives of the basis functions of the B-splines are directly applied to the linearized free surface condition. Numerical damping in the Dawson's method are avoided in the Present computations. In order to validate the present method, numerical computations are carried out for a submerged cylinder and a two-dimensional hydrofoil steadily moving beneath a free surface. The numerical results show that fast convergence and better accuracies have been achieved by the present method.