• 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 the Society of Naval Architects of Korea
• /
• v.32 no.1
• /
• pp.42-57
• /
• 1995

An efficient and accurate scheme has been constructed by taking advantages of the hi-quadratic spline scheme and the higher-order boundary element method selectively depending on computation domains. Boundary surfaces are represented by 8-node boundary elements to describe curved surfaces of a ship and its neighboring free surface more accurately. The variation of the velocity potential complies with the characteristics of the 8-node element on the body surface. But on the free surface, it is assumed to follow that of the hi-quadratic spline scheme. By which, the free surface solution is free from numerical damping and has better numerical dispersion property. As numerical examples, steady and unsteady Neumann-Kelvin problems are considered. Numerical results for a submerged spheroid, Series 60($C_B=0.6$) and a modified support the proposed method. Finally, a new upstream radiation condition is derived using a wave equation operator in order to deal with problems for subcritical reduced frequency. The relevance of this operator has been confirmed in the case of unsteady Kelvin source potential.

• 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 Society of Naval Architects of Korea
• /
• v.57 no.5
• /
• pp.271-277
• /
• 2020

Recoverability and vulnerability of navy ships under underwater explosion are critical verification factors in the acquisition phase of navy ships. This paper aims to establish numerical analysis techniques for the underwater explosion of navy ships. Doubly Asymptotic Approach (DAA) Equation of Motion (EOM) of primary shock wave and secondary bubble pulse proposed by Geers-Hunter was introduced. Assuming a non-compressive fluid, reference solution of the DAA EOM of Geers-Hunter using Runge-Kutta method was derived for the secondary bubble pulse phase with an assumed charge conditions. Convergence analyses to determine fluid element size were performed, suggesting that the minimum fluid element size for underwater explosion analysis was 0.1 m. The spherical and cylindrical fluid domains were found to be appropriate for the underwater explosion analyses from the fluid domain shape study. Because the element size of 0.1 m was too small to be applied to the actual navy ships, a very slender beam with the square solid section was selected for the study of fluid domain existence effect. The two underwater explosion models with/without fluid domain provided very similar results in terms of the displacement and stress processes.

• Journal of the Korean Society of Fisheries and Ocean Technology
• /
• v.31 no.4
• /
• pp.423-434
• /
• 1995

It was used boundary element method(BEM) and analysed axisymmetric problem to solve hydraulic cylinder for large vessel acting uniform internal pressure(25N/m super(2)) within elastic limit. This paper was utilized the carbon steel tubes for machine structural purposed model, inner radius was 150mm and outer radius was 250mm, axial length was semi-infinite and the isoparametric element was used. The important results obtained in this study were summarized as follows. Radial, tangential and shearing stress occured the maximum stresses(48, -20 and 34MPa) at the inner radius and the minimum stresses(32, -4 and 18MPa) at the outer radius of the hydraulic cylinder for large vessel. But negative signs have meaning compressive stress and stress diminution ratio was about 0.15MPa/mm. The use of isoparametric element raised accuracy and the increment of input data lessened the error in internal point but computer run-time was increased. The double node was improved the internal solutions to settle discontinuity at corner and the double exponential formula lessened error of stress value at boundary neighborhood. And then coincidence between the analytical and exact results is found to be fairly good, showing that the proposed analytical by BEM is reliable.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.25 no.3
• /
• pp.415-423
• /
• 2001

In is well recognized that the ultrasonic method is one of the most common and reliable nondestructive testing(NDT) methods for the quantitative estimation of defects in welded structures. However, NDT techniques applying for adhesively bonded joints have not been clearly established yet. In this paper, the detection of interface crack by the ultrasonic method was applied for the measurement of interface crack length in the adhesively bonded joints of double-cantilever beam(DCB). The optimum condition of transmission coefficients and experimental accuracy by the ultrasonic method in the adhesively bonded joints have been investigated. The experimental values are in good agreement with the computed results by boundary element method(BEM) and Riplings equation.

• Proceedings of the KSME Conference
• /
• 2000.04a
• /
• pp.97-102
• /
• 2000

It is well recognized that the ultrasonic methods is one of the most common and reliable nondestructive testing(NDT) methods for the quantitative estimation of defects in welded structures. However, NDT techniques applying for adhesively bonded joints have not been clearly established yet. In this paper, the detection of interface crack by the ultrasonic method was applied for the measurement of interfacial crack length in the adhesively bonded joints of double-cantilever beam(DCB). The optimum condition of transmission coefficients in the adhesively bonded joints and it's experimental accuracy by the ultrasonic method have been investigated. The experimental values are in good agreement with the computed results by boundary element method(BEM) and Ripling's equation.

• Journal of the Society of Naval Architects of Korea
• /
• v.36 no.3
• /
• pp.122-133
• /
• 1999

Overlap, excessive wide gaps, and errors between blocks during erection process increases cost and man-hour. In this paper, a system to simulate the assembly process is suggested and a plane block assembly of welding deformation with gravity is simulated in consideration of assembly order, deformations and errors occurring in the cutting and forming processes are not considered and welding deformations are acquired by equivalent stiffness and load method from experiments and hull double bottom plane block is assembled on a assembly order by panel method. It is certified that according to the order of assembly, intermediate product shape affects rigidity which affects welding deformations. Assembly order must be considered in the assembly process. It is certified that the gravity has important role in the assembly process.

• Journal of the Korean Institute of Gas
• /
• v.2 no.1
• /
• pp.23-27
• /
• 1998

To protect the underground pipeline from the mechanical damage and to enhance the cathodic protection effect, the river sand has been backfilled traditionally around the buried pipeline. However, river sand became depleted and expensive. One has to seek for the economic alternative materials. Crashed stone is a good candidate for the backfill material. In this study, how much the particle size and shape of the crashed stone can effect on the gas pipeline coating was examined. A series of laboratory and field test was performed. In the Lab, the increasing loads were applied to the coated pipeline surrounded by the crashed stone, where no significant damage was observed.

• Journal of the Korean Institute of Gas
• /
• v.2 no.1
• /
• pp.40-46
• /
• 1998

Mathematical modeling on the corrosion of the steel casing and main pipe due to the protection current resulting from a cathodic protection system was carried out using boundary element method. The model is consisted of Laplace's equation with non-linear boundary conditions(Tafel equations) and the iterative technique to determine the miexed potential of the steel casing. The model is applied to the normal steel casing section as well as abnormal one with defects such as metal touch and insulation defects. From the modeling procedure, we can calculate the potential distributions and current density distributions of the system. The theoretical results of the qualitatiive corrosion aspect along the steel casing and main pipe agree well with the experimental results within the experimental conditions studied.