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

• Proceedings of the KIEE Conference
• /
• 2000.07b
• /
• pp.828-830
• /
• 2000

This paper describes BIEM(Boundary Integral Equation Method) for computation of three dimensional electric field distribution and numerical method that an equivalent charge density is unknown variable. After computing numerically the surface charge distribution. the distribution of both potential and electric field are obtained. Finally, this numerical method is applied to the concentric sphere and the coaxial cylindrical model and numerical result is compared to the analytic solution.

• Tunnel and Underground Space
• /
• v.4 no.1
• /
• pp.17-23
• /
• 1994

This study consists of two procedures on back analysis and forward analysis which is a basic tool of the former. For a safe and economical construction of underground structures, it is required to identify the structural parameters and analyze the structural behavior as exactly as possible. In this paper, a boundary element method to analyze the behavior of multi-alyered underground structures is studied, in which body forces and initial stresses are considered. That is, each layer is discritized into subregions using infinite fundamental solutions, and terms of body forces and initial stresses are transformed into boundary integral where the applied direct integral method is used. And the system of equations containing body forces and initial stresses are considered. That is, each layer is discritized into subregions using infinite fundamental solutions, and terms of body forces and initial stresses are transformed into boundary integral where the applied direct integral method is used. And the system of equations containing body forces and initial stresses are composed, then the method to solve unknowns is used with applying compatibility and equilibrium conditions between interfaces. As well, the direct search method is applied in back analysis problems. By Powell's method as a technique to search unknown parameters, assuming displacements calculated from boundary element analysis as in-situ displacements, elastic moduli and initial stresses are presumed. As consequences of this study, the results of boundary element analysis of the behavior of multilayered structure considering body forces and initial stresses are agreed with those of finite element analysis. And results of back analysis of elastic moduli and initial stresses in each layers are agreed with exact values with a little difference. Therefore, it is known that this study can be efficiently applied for analyzing the behavior of underground structures including back analysis problems.

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

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.10 no.3
• /
• pp.277-284
• /
• 1986

A numerical procedure (integral method) for calculating the interacting turbulent boundary layer is set up. With this method, some free interactions with various upstream conditions are simulated in order to investigate the influence of upstream state on the interacting turbulent boundary layer. The results obtained by this numerical simulation can be summarized as follows; Free interaction of upstream unstabilized (or separated) turbulent boundary layer is subcritical regardless of its external Mach number, while free interaction of upstream stabilized turbulent boundary layer has two different characteristics (subcritical, supercritical) according to the external Mach number.

• Computational Structural Engineering
• /
• v.9 no.4
• /
• pp.127-134
• /
• 1996

This paper is about a progressive fracture analysis of concrete by boundary element method. From both displacement boundary integral equation and traction boundary integral equation of solids with cracks, a boundary integral equation for crack problem is derived. For the analysis of progressive fracture of concrete, fracture process zone is modelled based on Dugdale-Barenblatt model with linear tension-softening curve. By using the boundary element modeling, the progressive fractures of concrete beam and compact-tension specimens with various loading conditions are analyzed and compared with experiments. The analysis results show that the technique in this paper can predict the maximum strength and the nonlinear behavior of concrete including post-peak behavior.

• The Proceeding of the Korean Institute of Electromagnetic Engineering and Science
• /
• v.4 no.1
• /
• pp.38-43
• /
• 1993

A Hybrid Finite Element Method(HFEM) is applied to solve the electrormagnetic scattering from multi-layered dielectric cylinders. An unbounde region is divided into local boundary regions where a practical differential equation solution is obtained, with the remaining unbounded region represented by a boundary integral equation. If sources, media inhomogeneities, and anisotropies are local, a surgace may be defined to enclose them. Therefore the integral region so defined is bounded, and differential techniques may be used there. Also, in the re- maining unbounded region a boundary integral equation may be formulated using only a simple free - space green's function. Therefore, The local boundary is represented by a boundary - value problem with boundary conditions and solved by the finite element method. The advantage of the proposed method is simple and efficient in the work of electromagnetic scattering. The validity of the results have been verified by comparing results of other method(boundary element method). Examples has been presented to calculate the scattered fields of lossy dielectric cylinders of arbitray cross section.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.20 no.5
• /
• pp.1458-1467
• /
• 1996

A direct differentiationmethod is presented for the shape design sensitivity analysis of axisymmeetric elastic solids. Based on the exisymmetric boundary integralequaiton formulation, a new boundary ntegral equatio for sensitivity analysis is derived by taking meterial derivative to the same integral identity that was used in the adjoint variable melthod. Numerical implementation is performed to show the applicaiton of the theoretical formulation. For a simple example with analytic solution, the sensitivities by present method are compared with analytic sensitivities. As an application to the shape optimization, an optimal shape of a gas turbine disc toinimize the weight under stress constraints is found by incorporating the sensitivity analysis algorithm in an optimizatio program.

• Journal of Korean Association for Spatial Structures
• /
• v.1 no.2 s.2
• /
• pp.101-110
• /
• 2001

A large floating structure is attracting great attention in recent years from the view of ocean space utilization. Its huge scale in the horizontal directions compared with the wavelength and relatively shallow depth make this type of floating structure flexible and its wave-induced motion be characterized by the elastic deformation. In this paper, a boundary integral equation method is proposed to predict the wave-induced dynamic response mat-like floating offshore structure. The structure is modeled as an elastic plate and its elastic deformation is expressed as a superposition of free-vibration modes in air. This makes it straightforward to expand the well-established boundary integral technique for rigid floating bodies to include the hydroelastic effects. In order to validate the theoretical analysis, we compare with the experimental result of reduced model test. Satisfactory agreement is found between theory and experiment.

• Journal of the Korean Society of Industry Convergence
• /
• v.4 no.4
• /
• pp.433-439
• /
• 2001

A large floating structure is attracting great attention in recent years from the view of ocean space utilization. Its huge scale in the horizontal directions compared with the wavelength and relatively shallow depth make this type of floating structure flexible and its wave-induced motion be characterized by the elastic deformation. In this paper, a boundary integral equation method is proposed to predict the wave-induced dynamic response mat-like floating offshore structure. The structure is modeled as an clastic plate and its elastic deformation is expressed as a superposition of free-vibration modes in air. This makes it straightforward to expand the well-established boundary integral technique for rigid floating bodies to include the hydroelastic effects. In order to validate the theoretical analysis, we compare with the experimental result of previous model test. Satisfactory agreement is found between theory and experiment.