• Computational Structural Engineering
• /
• v.10 no.2
• /
• pp.213-223
• /
• 1997

Calculation of elastic wave fields has important applications in a variety of engineering fields including NDE (Non-destructive evaluation). Scattering problems have been investigated by numerous authors with different solution schemes. For simple geometries of the scatterers (e.g., cylinders or spheres), the analysis of steady-state elastic wave scattering has been carried out using analytical techniques. For arbitrary geometries and multiple inclusions, numerical methods have been developed. Special finite element methods, e.g., the infinite element method and a hybrid method called the Global-Local finite element method have also been developed for this purpose. Recently, the boundary integral equation method has been used successfully to solve scattering problems. In this paper, a volume integral equation method (VIEM) is proposed as a new numerical solution scheme for the solution of general elasto-dynamic problems in unbounded solids containing multiple inclusions and voids or cracks. A boundary integral equation method (BIEM) is also presented for elastic wave scattering problems. The relative advantage of the volume and boundary integral equation methods for solving scattering problems is discussed.

• International Journal of Safety
• /
• v.2 no.1
• /
• pp.39-44
• /
• 2003

The present study deals with an approximate integral equation approach to finite deflection of elastic plates with arbitrary plane form. An integral formulation leads to a system of boundary integral equations involving values of deflection, slope, bending moment and transverse shear force along the edge. The basic principles of the development of boundary element technique are reviewed. A computer program for solving for stresses and deflections in a isotropic, homogeneous, linear and elastic bending plate is developed. The fundamental solution of deflection and moment is employed in this program. The deflections and moments are assumed constant within the quadrilateral element. Numerical solutions for sample problems, obtained by the direct boundary element method, are presented and results are compared with known solutions.

• Journal of the Korean Society for Precision Engineering
• /
• v.20 no.8
• /
• pp.175-184
• /
• 2003

A method of the shape design sensitivity analysis based on the boundary integral equation formulation is presented for two-dimensional inhomogeneous thermal conducting solids with multiple domains. Shape variation of the external and interface boundary is considered. A sensitivity formula of a general performance functional is derived by taking the material derivative to the boundary integral identity and by introducing an adjoint system. In numerical analysis, state variables of the primal and adjoint systems are solved by the boundary element method using quadratic elements. Two numerical examples of a compound cylinder and a thermal diffuser are taken to show implementation of the shape design sensitivity analysis. Accuracy of the present method is verified by comparing analyzed sensitivities with those by the finite difference. As application to the shape optimization, an optimal shape of the thermal diffuser is found by incorporating the sensitivity analysis algorithm in an optimization program.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.23 no.2
• /
• pp.243-253
• /
• 1999

The dual reciprocity boundary element method (DRBEM) is used to solve the Graetz problem of laminar flow inside circular duct. In this method the domain integral tenn of boundary integral equation resulting from source term of governing equation is transformed into equivalent boundary-only integrals by using the radial basis interpolation function, and therefore complicate domain discretization procedure Is completely removed. Velocity profile is obtained by solving the momentum equation first and then, using this velocities as Input data, energy equation Is solved to get the temperature profile by advancing from duct entrance through the axial direction marching scheme. DRBEM solution is tested for the uniform temperature and heat flux boundary condition cases. Local Nusselt number, mixed mean temperature and temperature profile inside duct at each dimensionless axial location are obtained and compared with exact solutions for the accuracy test Solutions arc in good agreement at the entry region as well as fully developed region of circular duct, and their accuracy are verified from error analysis.

• Journal of the Chungcheong Mathematical Society
• /
• v.7 no.1
• /
• pp.141-152
• /
• 1994

We consider a numerical scheme for solving a nonlinear boundary integral equation (BIE) obtained by reformulation the nonlinear boundary value problem (BVP). We give a simple alternative to the standard collocation method for the nonlinear BIE. This method consists of one conventional linear system and another coupled linear system resulting from an auxiliary BIE which is obtained by differentiating both side of the nonlinear interior integral equations. We obtain an analogue BIE through the perturbation of the fundamental solution of Laplace's equation. We procure the super-convergence of approximate solutions.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.26 no.4
• /
• pp.775-786
• /
• 2002

A Mixed Volume and Boundary Integral Equation Method is applied for the effective analysis of elastic wave scattering problems and 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. In the formulation of this method, the continuity condition at each interface is automatically satisfied, and in contrast to finite element methods, where the full domain needs to be discretized, this method requires discretization of the inclusions only. Finally, this method takes full advantage of the pre- and post-processing capabilities developed in FEM and BIEM. Through the analysis of plane elastostatic problems in unbounded isotropic matrix with orthotropic inclusions and voids or isotropic inclusions, and the analysis of plane wave scattering problems in unbounded isotropic matrix with isotropic inclusions and voids, it will be established that this new method is very accurate and effective for solving plane wave scattering problems and plane elastic problems in unbounded solids containing general anisotropic inclusions and voids/cracks or isotropic inclusions.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.2 no.1
• /
• pp.111-129
• /
• 1998

An integral equation representation of cracks was presented, which differs from well-known "dislocation layer" representation. In this new representation, an integral equation representation of cracks was developed and coupled to the direct boundary-element method for treatment of cracks in plane finite bodies. The method was developed for in-plane(modes I and II) loadings only. In this paper, the method is formulated and applied to mode III problems involving smooth or kinked cracks in finite region. The results are compared to exact solutions where available and the method is shown to be very accurate despite of its simplicity.

• Structural Engineering and Mechanics
• /
• v.15 no.5
• /
• pp.535-550
• /
• 2003

In this paper, a boundary-type meshfree method, the boundary radial point interpolation method (BRPIM), is presented for solving boundary value problems of two-dimensional solid mechanics. In the BRPIM, the boundary of a problem domain is represented by a set of properly scattered nodes. A technique is proposed to construct shape functions using radial functions as basis functions. The shape functions so formulated are proven to possess both delta function property and partitions of unity property. Boundary conditions can be easily implemented as in the conventional Boundary Element Method (BEM). The Boundary Integral Equation (BIE) for 2-D elastostatics is discretized using the radial basis point interpolation. Some important parameters on the performance of the BRPIM are investigated thoroughly. Validity and efficiency of the present BRPIM are demonstrated through a number of numerical examples.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2005.11a
• /
• pp.915-918
• /
• 2005

Acoustic holography uses Kirchhoff·Helmholtz integral equation and Green's function which satisfies Dirichlet boundary condition Applications of acoustic holography have been taken to the sound field neglecting the effect of flow. The uniform flow, however, changes sound field and the governing equation, Green's function and so on. Thus the conventional method of acoustic holography should be changed. In this research, one possibility to apply acoustic holography to the sound field with uniform flow is introduced through checking for the plane wave in a duct. Change of Green's function due to uniform flow and one method to derive modified form of Kirchhoff·Heimholtz integral is suggested for 1-dimensional sound field. Derivation results show that using Green's function satisfying Dirichlet boundary condition, we can predict sound pressure in a duct using boundary value.

• Proceedings of the KIEE Conference
• /
• 1999.11b
• /
• pp.101-103
• /
• 1999

This paper presents an analysis of 2-dimensional(2-D) multi-regions problem using boundary integral equation method(BIEM). When compared with finite element method(FEM), there are only a few unknown variables in BIEM because it implements numerical analysis only for the surface or boundary of a model. As a result, a lot of computational memory and time can be saved. Procedure to analyze 2-D multi-regions problem using potentials and its derivatives in a boundary as unknown variables, first, numerical analysis is performed for each of subregions. And then interface continuity condition is applied to the interface between them and Gauss Quadrature Formula are adopted to solve singular integral in a boundary in this paper.