• Journal of the Korean Institute of Navigation
• /
• v.10 no.1
• /
• pp.129-138
• /
• 1986

A matrix partitioning method is proposed for the 2-D motion analysis of floating bodies. For the numerical solution, the boundary of a floating body is approximated with a series of line segments and the governing integral equation is transformed into a system of linear equations. A new solution procedure of resulting linear equation with complex coefficients is formulated and programmed using a matrix partitioning scheme and the Choleski decomposition. From the case study, it is found that the proposed method is efficient in the motion analysis of floating bodies, especially in the calculation of hydrodynamic coefficients. Also, it requires smaller memory size and less computing time compared with conventional methods.

• Magazine of the Korea Concrete Institute
• /
• v.8 no.6
• /
• pp.205-212
• /
• 1996

This paper presents progressive fracture analysis of concrete using boundary element method and its stabilizing technique. To determine ultimate strength and to predict nonlinear behavior of concrete during progressive crack growth, the modelling of fracture process zone is done based on Dugdale-Barenblatt model with linear tension-softening curve. We regulate displacement and traction boundary integral equation of solids including crack boundary and analyze progressive fracture of concrete beam and compact tension specimen. Also a numerical technique which considers the growth of stress-free crack of concrete during the analysis and removes snapback of postpeak behavior is proposed.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.14 no.5
• /
• pp.1033-1042
• /
• 1994

Dynamic behavior of flexible rectangular liquid containers is analyzed by a two-dimensional coupled boundary element-finite element method. The irrotational motion of inviscid and incompressible ideal fluid is modeled by boundary elements and the motion of structure by finite elements. A singularity free integral formulation is employed for the implementation of boundary element method. Coupling is performed by using compatibility and equilibrium conditions along the interface between the fluid and structure. The fluid-structure interaction effects are reflected into the coupled equation of motion as added fluid mass matrix and sloshing stiffness matrix. By solving the eigen-problem for the coupled equation of motion, natural frequencies and mode shapes of coupled system are obtained. The free surface sloshing motion and hydrodynamic pressure developed in a flexible rectangular container due to horizontal and vertical ground motions are computed in time domain.

• Communications of the Korean Mathematical Society
• /
• v.37 no.3
• /
• pp.939-955
• /
• 2022

In this paper, we study a first-order non-linear singularly perturbed Volterra integro-differential equation (SPVIDE). We discretize the problem by a uniform difference scheme on a Bakhvalov-Shishkin mesh. The scheme is constructed by the method of integral identities with exponential basis functions and integral terms are handled with interpolating quadrature rules with remainder terms. An effective quasi-linearization technique is employed for the algorithm. We establish the error estimates and demonstrate that the scheme on Bakhvalov-Shishkin mesh is O(N^-1) uniformly convergent, where N is the mesh parameter. The numerical results on a couple of examples are also provided to confirm the theoretical analysis.

• Journal of the Society of Naval Architects of Korea
• /
• v.47 no.3
• /
• pp.291-305
• /
• 2010

A three-dimensional time-domain calculation method is of crucial importance in prediction of the motions and wave loads of a ship advancing in a severe irregular sea. The exact solution of the free surface wave-ship interaction problem is very complicated because of the essentially nonlinear boundary conditions. In this paper, an approximate body nonlinear approach based on the three-dimensional time-domain forward-speed free-surface Green function has been presented. The Froude-Krylov force and the hydrostatic restoring force are calculated over the instantaneous wetted surface of the ship while the forces due to the radiation and scattering potentials over the mean wetted surface. The time-domain radiation and scattering potentials have been obtained from a time invariant kernel of integral equations for the potentials which are discretized according to the second-order boundary element method (Hong and Hong 2008). The diffraction impulse-response functions of the Wigley seakeeping model advancing in transient head waves at various Froude numbers have been presented. A simulation of coupled heave-pitch motion of a long rectangular barge advancing in regular head waves of large amplitude has been carried out. Comparisons between the linear and the approximate body nonlinear numerical results of motions and wave loads of the barge at a nonzero Froude number have been made.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.27 no.3
• /
• pp.38-46
• /
• 1990

A numerical solution method of the boundary integral equation for axisymmetric potential flows is presented. Those are represented by ring source and ring vorticity distribution. Strengths of ring source and ring vorticity are approximated by linear functions of a parameter $\zeta$ on a segment. The geometry of the body is represented by a cubic B-spline. Limiting integral expressions as the field point tends to the surface having ring source and ring vorticity distribution are derived upto the order of ${\zeta}ln{\zeta}$. In numerical calculations, the principal value integrals over the adjacent segments cancel each other exactly. Thus the singular part proportional to $\(\frac{1}{\zeta}\)$ can be subtracted off in the calculation of the induced velocity by singularities. And the terms proportional to $ln{\zeta}$ and ${\zeta}ln{\zeta}$ can be integrated analytically. Thus those are subtracted off in the numerical calculations and the numerical value obtained from the analytic integrations for $ln{\zeta}$ and ${\zeta}ln{\zeta}$ are added to the induced velocity. The four point Gaussian Quadrature formula was used to evaluate the higher order terms than ${\zeta}ln{\zeta}$ in the integration over the adjacent segments to the field points and the integral over the segments off the field points. The root mean square errors, $E_2$, are examined as a function of the number of nodes to determine convergence rates. The convergence rate of this method approaches 2.

• Journal of electromagnetic engineering and science
• /
• v.12 no.4
• /
• pp.260-270
• /
• 2012

This work aims to present a theoretical analysis of the electric and magnetic surface current densities of a circular microstrip antenna (CMSA) as a body of revolution. The rigorous analysis of these problems begins with the application of the equivalence principle, which introduces an unknown electric current density on the conducting surface and both unknown equivalent electric and magnetic surface current densities on the dielectric surface. These current densities satisfy the integral equations (IEs) obtained by canceling the tangential components of the electric field on the conducting surface and enforcing the continuity of the tangential components of the fields across the dielectric surface. The formulation of the radiation problems is based on the combined field integral equation. This formulation is coupled with the method of moments (MoMs) as a numerical solution for this equation. The numerical results of the electric and magnetic surface current densities on the outside boundary of a CMSA excited by $TM_{11^-}$ and $TM_{21^-}$ modes are presented. The radiation pattern is calculated numerically in the two principle planes for a CMSA and gives a good results compared with measured results published by other research workers.

• Structural Engineering and Mechanics
• /
• v.54 no.4
• /
• pp.607-622
• /
• 2015

This paper presents a comparative study of analytical method and finite element method (FEM) for analysis of a continuous contact problem. The problem consists of two elastic layers loaded by means of a rigid circular punch and resting on semi-infinite plane. It is assumed that all surfaces are frictionless and only compressive normal tractions can be transmitted through the contact areas. Firstly, analytical solution of the problem is obtained by using theory of elasticity and integral transform techniques. Then, finite element model of the problem is constituted using ANSYS software and the two dimensional analysis of the problem is carried out. The contact stresses under rigid circular punch, the contact areas, normal stresses along the axis of symmetry are obtained for both solutions. The results show that contact stresses and the normal stresses obtained from finite element method (FEM) provide boundary conditions of the problem as well as analytical results. Also, the contact areas obtained from finite element method are very close to results obtained from analytical method; disagree by 0.03-1.61%. Finally, it can be said that there is a good agreement between two methods.

• Journal of the Society of Naval Architects of Korea
• /
• v.31 no.4
• /
• pp.82-98
• /
• 1994

A mathematical model for the hydrodynamic forces acting on the ship manoeuvring in waves is formulated and a numerical method for the problem is developed. The motion of a ship, which manoeuvres in waves, may be thought to have two components; one is a high frequency component due to encounter waves, and the other is a low frequency component due to manoeuvring motion. So the method of two time scale expansion is used to divide linear boundary value problem. For the effects of waves on the manoeuvring motion of a ship, only the second order drift forces are considered. The integral equation for the velocity potential is solved by 3 dimensional panel method and hydrodynamic forces are calculated by direct integral method.

• Proceedings of the Optical Society of Korea Conference
• /
• 1990.02a
• /
• pp.64-68
• /
• 1990

Boundary element method using linear basis function is applied to obtain fields scattered from a 3-D dielectric sphere. Electric field integral equation is used on the surfaces of the dielectric material where its surface is discretized into trilateral cells. For plane wave incidence, scattered fields by a dielectric sphere is calculated and compared with its analytic solution. The total electric fields are calculated on the great circle of the sphere boundary as well as the outside of the sphere in the plane of the wave vector and the polarization vector of the incident electric field.