• Journal of Mechanical Science and Technology
• /
• v.14 no.2
• /
• pp.197-206
• /
• 2000

A boundary integral equation method in the shape design sensitivity analysis is developed for the elasticity problems with axisymmetric non-homogeneous bodies. Functionals involving displacements and tractions at the zonal interface are considered. Sensitivity formula in terms of the interface shape variation is then derived by taking derivative of the boundary integral identity. Adjoint problem is defined such that displacement and traction discontinuity is imposed at the interface. Analytic example for a compound cylinder is taken to show the validity of the derived sensitivity formula. In the numerical implementation, solutions at the interface for the primal and adjoint system are used for the sensitivity. While the BEM is a natural tool for the solution, more generalization should be made since it should handle the jump conditions at the interface. Accuracy of the sensitivity is evaluated numerically by the same compound cylinder problem. The endosseous implant-bone interface problem is considered next as a practical application, in which the stress value is of great importance for successful osseointegration at the interface. As a preliminary step, a simple model with tapered cylinder is considered in this paper. Numerical accuracy is shown to be excellent which promises that the method can be used as an efficient and reliable tool in the optimization procedure for the implant design. Though only the axisymmetric problem is considered here, the method can be applied to general elasticity problems having interface.

• Communications of the Korean Mathematical Society
• /
• v.38 no.1
• /
• pp.281-298
• /
• 2023

A class of systems of Caputo fractional differential equations with integral boundary conditions is considered. A numerical method based on a finite difference scheme on a uniform mesh is proposed. Supremum norm is used to derive an error estimate which is of order κ − 1, 1 < κ < 2. Numerical examples are given which validate our theoretical results.

• Kyungpook Mathematical Journal
• /
• v.47 no.4
• /
• pp.529-548
• /
• 2007

The purpose of this article is to find a relation between the finite difference method and the boundary element method, and propose a new approach deriving a discrete approximation formula as like that of the finite difference method for harmonic functions. We develop a discrete approximation formula on a uniform grid based on the boundary integral formulations. We consider three different boundary integral formulations and derive one discrete approximation formula on the uniform grid for the harmonic function. We show that the proposed discrete approximation formula has the same computational molecules with that of the finite difference formula for the Laplace operator ${\nabla}^2$.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.20 no.7
• /
• pp.2167-2173
• /
• 1996

An approximate weight function technique using the indirect boundary integral equation has been presented for the analysis of stress intensity foactors(SIFs) of a thick pipe. One-term boundary integral was introduced to represent the crack surface displacement field for the displacement based weight function technique. An explicit closed-form SIF solution applicable to symmetric cracked pipes without any modification of the solution including both circumferential and radial cracks has been derived. The necessary information in the analysis is two or three reference SIFs. In most cases the SIF solution were in good agreement with those available in the literature.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 1997.04a
• /
• pp.701-708
• /
• 1997

A direct boundary element method(DBEM) is developed for thin bodies whose surfaces are rigid or compliant. Th eHelmholtz integral equation and its normal derivative integral equation are adopted simultaneously to calculate the pressure on both sides of the thin body, instead of the jump values across it, to account for the different surface conditions of each side. Unlike the usual assumption, the normal velocity is assumed to be discontinuous across the thin body. In this approach, only the neutral surface of the thin body has to be discontinuous across the thin body. In this approach, only the neural surface of the thin body has to be discretized. The method is validated by comparison with analytic and/or numerical results for acoustic scattering and radiation from several surface conditions of the thin body; the surfaces are rigid when stationary or vibrating, and part of the interior surface is lined with a sound-absorbing material.

• Structural Engineering and Mechanics
• /
• v.10 no.4
• /
• pp.393-404
• /
• 2000

Linear stress analysis without body force can be easily solved by means of the boundary element method. Some cases of linear stress analysis with body force can also be solved without a domain integral. However, domain integrals are generally necessary to solve the linear stress problem with arbitrary body forces. This paper shows that the linear stress problem with arbitrary body forces can be solved approximately without a domain integral by the triple-reciprocity boundary element method. In this method, the distribution of arbitrary body forces can be interpolated by the integral equation. A new computer program is developed and applied to several problems.

• Journal of the Society of Naval Architects of Korea
• /
• v.56 no.1
• /
• pp.11-22
• /
• 2019

In this paper, a Rankine source method is applied and validated to analyze the hydrodynamic response of a three-dimensional floating structure in the frequency domain. The boundary value problems for radiation and diffraction problem are solved by using a desingularized indirect boundary integral equation method (DIBIEM). The DIBIEM is simpler and faster than conventional methods based on the numerical surface integration of Green's function because the singularities of Green's function are located outside of fluid regions. In case of floating structure with complex geometry, it is difficult to desingularize the singularities of Green's function consistently. Therefore a mixed approach is carried out in this study. The mixed approach is partially desingularized except singularities of the body. Wave drift loads are calculated by the middle-field formulation method that is mathematically simple and has fast convergence. In order to validate the accuracy of the developed program, various numerical simulations are carried out and these results are analyzed and compared with previously published calculations and experiments.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.14 no.4
• /
• pp.303-309
• /
• 2004

Absorptive material arrangement method for effective interior noise control is proposed. Sound field with arbitrary boundary condition is formulated by Kirchhoff-Helmholtz integral equation. A simple example such as a rectangular cavity will present physical meaning between changing boundary condition and control of sound field. The effect of changing boundary condition is expressed in modal admittance. From this formulation. an admittance map is presented. The admittance map is the figure to represent position where absorptive material is attached for effective interior noise control. The admittance map can be assigned to each resonant frequency. There. however, may be common area of those maps. Then, frequency robust arrangement of absorptive material in noise control will be presented.

• 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.

• Proceedings of the KSME Conference
• /
• 2003.11a
• /
• pp.1702-1707
• /
• 2003

Absorptive material arrangement method for effective interior noise control is proposed. Sound field with arbitrary boundary condition is formulated by Kirchhoff-Helmholtz integral equation. A simple example such as a rectangular cavity will present physical meaning between changing boundary condition and control of sound field. The effect of changing boundary condition is expressed in modal admittance. From this formulation, an admittance map is presented. The admittance map is the figure to represent position where absorptive material is attached. The admittance map can be assigned to each resonant frequency. There, however, may be common area of those maps. Then, frequency robust arrangement of absorptive material in noise control will be presented.