• Structural Engineering and Mechanics
• /
• v.30 no.3
• /
• pp.279-296
• /
• 2008

This paper presents an exact integration for the hypersingular boundary integral equation of two-dimensional elastostatics. The boundary is discretized by straight segments and the physical variables are approximated by discontinuous quadratic elements. The integral for the hypersingular boundary integral equation analysis is given in a closed form. It is proven that using the exact integration for discontinuous boundary element, the singular integral in the Cauchy Principal Value and the hypersingular integral in the Hadamard Finite Part can be obtained straightforward without special treatment. Two numerical examples are implemented to verify the correctness of the derived exact integration.

• Structural Engineering and Mechanics
• /
• v.30 no.3
• /
• pp.317-330
• /
• 2008

Stress intensity factors for a planar crack parallel to a bimaterial interface are considered. The formulation leads to a system of hypersingular integral equations whose unknowns are three modes of crack opening displacements. In the numerical analysis, the unknown displacement discontinuities are approximated by the products of the fundamental density functions and polynomials. The numerical results show that the present method yields smooth variations of stress intensity factors along the crack front accurately. The mixed mode stress intensity factors are indicated in tables and figures with varying the shape of crack, distance from the interface, and elastic constants. It is found that the maximum stress intensity factors normalized by root area are always insensitive to the crack aspect ratio. They are given in a form of formula useful for engineering applications.

• Journal of the Korean Mathematical Society
• /
• v.41 no.2
• /
• pp.345-368
• /
• 2004

A preconditioning algorithm is developed in this paper for the iterative solution of the linear system of equations resulting from the p-version boundary element approximation of the three dimensional integral equation with hypersingular operators. The preconditioner is derived by first making the nodal and side basis functions locally orthogonal to the element internal bases, and then by decoupling the nodal and side bases from the internal bases. Its implementation consists of solving a global problem on the wire-basket and a series of local problems defined on a single element. Moreover, the condition number of the preconditioned system is shown to be of order $O((1＋ln/p)^{7})$. This technique can be applied to discretization with triangular elements and with general basis functions.

• Korean Journal of Mathematics
• /
• v.18 no.4
• /
• pp.425-440
• /
• 2010

Using Green's theorem, elliptic boundary value problems can be converted to boundary integral equations. A numerical methods for boundary integral equations are boundary elementary method(BEM). BEM has advantages over finite element method(FEM) whenever the fundamental solutions are known. Helmholtz type equations arise naturally in many physical applications. In a boundary integral formulation for the exterior Neumann there occurs a hypersingular operator which exhibits a strong singularity like $\frac{1}{|x-y|^3}$ and hence is not an integrable function. In this paper we are going to remove this hypersingularity by reducing the regularity of test functions.

• Journal of the Korean Society for Precision Engineering
• /
• v.12 no.10
• /
• pp.89-101
• /
• 1995

선형 탄성 등방성 물체 내에 있는 일반적인 복합모드 크랙 문제들을 해석하기 위한 이중 경계적분방정식의 일반식과 계산해법이 제시되었다. 크랙면이 포함된 물체 해석에 있어서 유일한 해를 얻기 위하여, 한 면상의 점에는 변위 경계적분방정식이 적용되었고 마주하고 있는 상대면 상의 점에는 인력 경계적분방정식이 적용되었다. 인력 및 변위 경계적분방정식의 강특이해 및 초특이해 적분항들은 수치해법을 적용하기 전에 정상화되었다. 정상화과정 중 보정되는 강특이적분항이 상대 크랙면 상의 특이해 요소를 따라 직접 적분되는 것을 격리시키기 위하여, 특이해 적분 경로를 완만한 곡면으로 우회시킨 가상의 비특이해 보조경계로 대치하여 적분값을 계산하였다. 제시된 해법의 정확성과 효율성을 예시하기 위하여, 2차원 및 3차원 크랙 문제의 변형 후 모습과 응력강도계수 계산 결과를 보였다.

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.26 no.3
• /
• pp.297-307
• /
• 2020

In the study, energy flow analysis is performed to predict the performance of silencers. To date, deterministic approaches such as finite element method have been widely used for silencer analysis. However, they have limitations in analyzing large structures and mid-high frequency ranges due to unreasonable computational costs and errors. However, silencers used for ships and off-shore plants are much larger than those used in other engineering fields. Hence, energy governing equation, which is significantly efficient for systems with high modal density, is solved for silencers in ships and off-shore plants. The silencer is divided into two different acoustic media, air and absorption materials. The discontinuity of energy density at interfaces is solved via hypersingular integrals for the 3-D modified Helmholtz equation to analyze multi-domain problems with the energy flow boundary element method. The method is verified by comparing the measurements and analysis results for ship silencers over mid-high frequency ranges. The comparisons confirm good agreement between the measurement and analysis results. We confirm that the applied analysis method is useful for large silencers in mid-high frequency ranges. With the proven procedures, energy flow analysis can be performed for various types of silencer used in ships and off-shore plants in the first stage of the design.

• Journal of electromagnetic engineering and science
• /
• v.17 no.2
• /
• pp.91-97
• /
• 2017

An accurate and efficient evaluation for hypersingular integrals (HIs), strongly singular integrals (SSIs), and weakly singular integrals (WSIs) plays an essential role in the numerical solutions of 3D electromagnetic scattering problems. We derive analytic formulas for WSIs based on Stokes' theorem, which can be expressed in elementary functions. Several numerical examples are presented to validate these analytic formulas. Then, to show the feasibility of the proposed formulations for numerical methods, these formulations are used with the existing analytical expressions of HIs and SSIs to correct the near-field interaction in an iterative physical optics (IPO) scheme. Using IPO, the scattering caused by a dihedral reflector is analyzed and compared with the results of the method of moments and measurement data.