• Journal of the Korean Society for Precision Engineering
• /
• v.19 no.5
• /
• pp.81-88
• /
• 2002

An interface crack of a piezoelectric medium bonded between an elastic layer and a half-space is analyzed using the theory of linear piezoelectricity. Both out-of-plane mechanical and in-plane electrical loads are applied to the piezoelectric laminate. By the use of courier transforms, the mixed boundary value problem is reduced to a singular integral equation which is solved numerically to determine the stress intensity factors. Numerical analyses for various material combinations are performed and the results are discussed.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.30 no.8 s.251
• /
• pp.949-956
• /
• 2006

The method of Greenwood and Williamson is extended to obtain a solution to the coupled non-linear problem of steady-state electrical and thermal conduction across a crack in a conductive layer, for which the electrical resistivity and thermal conductivity are functions of temperature. The problem can be decomposed into the solution of a pair of non-linear algebraic equations involving boundary values and material properties. The new mixed-boundary value problem given from the thermal and electrical boundary conditions for the crack in the conductive layer is reduced in order to solve a singular integral equation of the first kind, the solution of which can be expressed in terms of the product of a series of the Chebyshev polynomials and their weight function. The non-existence of the solution for an infinite conductor in electrical and thermal conduction is shown. Numerical results are given showing the temperature field around the crack.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.25 no.2
• /
• pp.11-18
• /
• 1988

The radiation problem for an oscillating cylinder advancing under the free surface with a constant horizontal velocity is studied using the Green integral equation in the frequency domain. The Green function expressed in terms of the complex exponential, is derived using the damped free surface condition. Special attention is given to the behavior of the numerical solution in the vicinity of the critical Brard number ${\gamma}_c=\omega{\cdot}u/g=0.25$ where $\omega$ is the circular frequency of encounter, u the advancing speed and g the gravitational acceleration. It is shown that the solution is finite in the vicinity of ${\gamma}_c$ although the Green function becomes singular at ${\gamma}_c$. It is also shown that the computed hydrodynamic coefficients agree well with those obtained from the solution of the same problem formulated in the time domain.

• Journal of the Chungcheong Mathematical Society
• /
• v.24 no.2
• /
• pp.247-252
• /
• 2011

We study some properties of the Riemann-Stieltjes integrals with respect to the Riesz-$N\acute{a}gy$-$Tak\acute{a}cs$ distribution $H_{a,p}$ and its inverse $H_{p,a}$ on the unit interval satisfying the equation 1 - a = $a^2$ and p = 1 - a. Using the properties of the dual distributions $H_{a,p}$ and $H_{p,a}$, we compare the Riemann-Stieltjes integrals of $H_{a,p}$ over some essential intervals with that of its inverse $H_{p,a}$ over the related intervals.

• 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 Welding and Joining
• /
• v.29 no.2
• /
• pp.64-71
• /
• 2011

In this paper, crack propagation analyses in the inner diameter (ID) repair weld of the dissimilar metal weldment of a nozzle were performed using a finite element alternating method (FEAM). To calculate the theoretical solution for the crack tip stress intensity factor, a weak type singular integral equation consisted of crack surface traction and dislocation density function was constructed and solved in conjunction with the FEAM. A two-dimensional axisymmetric finite element nozzle model was prepared and ID repair welding was simulated. An initial crack, 10% depth of weld thickness, was assumed and crack propagation trajectory from the initial crack to the 75% depth of thickness was calculated using the FEAM. Crack growth versus time curve was also calculated and compared with the curves obtained from ASME code method. With the method constructed in this paper, crack propagation trajectory and crack growth time were calculated automatically and effectively.

• Journal of the Society of Naval Architects of Korea
• /
• v.28 no.1
• /
• pp.38-52
• /
• 1991

In this paper, a semi-Lagrangian method is used to solve the nonlinear hydrodynamics of a three-dimensional body beneath the free surface in the time domain. The boundary value problem is solved by using the boundary integral method. The geometries of the body and the free surface are represented by the curved panels. The surfaces are discretized into the small surface elements using a bi-cubic B-spline algorithm. The boundary values of $\phi$ and $\frac{\partial{\phi}}{\partial{n}}$ are assumed to be bilinear on the subdivided surface. The singular part proportional to $\frac{1}{R}$ are subtracted off and are integrated analytically in the calculation of the induced potential by singularities. The far field flow away from the body is represented by a dipole at the origin of the coordinate system. The Runge-Kutta 4-th order algorithm is employed in the time stepping scheme. The three-dimensional form of the integral equation and the boundary conditions for the time derivative of the potential Is derived. By using these formulas, the free surface shape and the equations of motion are calculated simultaneously. The free surface shape and fille forces acting on a body oscillating sinusoidally with large amplitude are calculated and compared with published results. Nonlinear effects on a body near the free surface are investigated.