• Journal of applied mathematics & informatics
• /
• v.7 no.3
• /
• pp.895-904
• /
• 2000

The contact problem of two elastic bodies of arbitrary shape with a general kernel form, investigated from Hertz problem, is reduced to an integral equation of the second kind with Cauchy kernel. A numerical method is adapted to determine the unknown potential function between the two surfaces under certain conditions. Many cases are derived and discussed from the work.

• The Transactions of The Korean Institute of Electrical Engineers
• /
• v.60 no.2
• /
• pp.357-361
• /
• 2011

The terminal current in voltage driven systems is an essential role for characterizing the pattern of electric discharge such as corona, breakdown, etc. Until now, to evaluate this terminal current, Sato's equation has been widely used in areas of high voltage and plasma discharge. Basically Sato's equation was derived by using the energy balance equation and its final form described physical meaning explicitly. To give more general abilities in Sato's equation, we present a generalized approach by directly using the Poynting's theorem incorporating the finite element method. When the magnetic field effect or the time-dependent voltage source is considered, this generalized energy method can be easily applicable to those problems with any dielectric media such as gas, fluid, and solid. As an alternative approach, the integral Ohm's law resulting in small numerical errors has an ability to be applied to multi-port systems. To test the generalized energy method and integral Ohm's law, first, the results from two prosed methods were compared to those from Sato's approach and an analytic solution in parallel plane electrodes. After verification, the generalized method was applied to the tip-sphere electrodes for evaluating the terminal current with three carriers and the Fowler-Nordheim field emission condition. From these results, we concluded that the generalized energy method can be a consistent technique for evaluating the discharge current with various dielectric materials or large magnetic field.

• Structural Engineering and Mechanics
• /
• v.43 no.5
• /
• pp.561-582
• /
• 2012

Despite popularity of FEM in analysis of static and dynamic structural problems and the routine applicability of FE softwares, analytical methods based on simple mathematical relations is still largely sought by many researchers and practicing engineers around the world. Development of such analytical methods for analysis of free vibration of non-prismatic beams is also of primary concern. In this paper a new and simple method is proposed for determination of vibration frequencies of non-prismatic beams under variable axial forces. The governing differential equation is first obtained and, according to a harmonic vibration, is converted into a single variable equation in terms of location. Through repetitive integrations, integral equation for the weak form of governing equation is derived. The integration constants are determined using the boundary conditions applied to the problem. The mode shape functions are approximated by a power series. Substitution of the power series into the integral equation transforms it into a system of linear algebraic equations. Natural frequencies are determined using a non-trivial solution for system of equations. Presented method is formulated for beams having various end conditions and is extended for determination of the buckling load of non-prismatic beams. The efficiency and convergence rate of the current approach are investigated through comparison of the numerical results obtained to those obtained using available finite element software.

• Journal of Advanced Navigation Technology
• /
• v.7 no.2
• /
• pp.180-190
• /
• 2003

By using accurate closed-form Green's functions obtained from real-axis integration method, the full-wave analysis of CPW discontinuities are performed in space domain. In solving MPIE(Mixed Potential Integral Equation), Galerkin's scheme is employed with the linear basis functions on the triangular elements in air-dielectric boundary. In the singular integral arising when the observation point and source point coincides, the surface integral is transformed into the line integral and the integral is evaluated by regular integration. By using the Green's function from the real-axis integration method, the discontinuities are characterized accurately.

• Proceedings of the Korea Electromagnetic Engineering Society Conference
• /
• 2003.11a
• /
• pp.54-58
• /
• 2003

In this study, the solution of Gel'fand-Levitan-Marchenko integral equation in the inverse scattering problem is efficiently solved for arbitrarily specified spectra pattern which are reflected from the restricted potential. The procedure is based on the successive approach without iterations. This method lessens the truncation errors which plague conventional design schemes using specific windows for reflection coefficients. It is shown that the method is adequate for the synthesis of the continuously varying one-dimensional potential of the nonuniformly distributed dielectric constants.

• Journal of Electrical Engineering and Technology
• /
• v.2 no.1
• /
• pp.129-135
• /
• 2007

In this paper, we present a set of numerical schemes to solve the Muller integral equation for the analysis of electromagnetic scattering from arbitrarily shaped three-dimensional (3-D) dielectric bodies by applying the method of moments (MoM). The piecewise homogeneous dielectric structure is approximated by planar triangular patches. A set of the RWG (Rao, Wilton, Glisson) functions is used for expansion of the equivalent electric and magnetic current densities and a combination of the RWG function and its orthogonal component is used for testing. The objective of this paper is to illustrate that only some testing procedures for the Muller integral equation yield a valid solution even at a frequency corresponding to an internal resonance of the structure. Numerical results for a dielectric sphere are presented and compared with solutions obtained using other formulations.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.12
• /
• pp.3219-3226
• /
• 1994

A model is constructed to evaluate the stress intensity factors for a center crack subjected to polynomial anti-symmetric loading in a layered material. A Fredholm integral equation is derived by Fourier integral transform method. The integral equation is numerically analyzed to evaluate the effects of the ratios of shear modulus, Poisson's ratio and crack length to layer thickness as well as the number of layers on the stress intensity factor. The stress intensity factors are approached to constant values as the number of layers increase and decrease as the polynomial power of the loading increase. In case of the E-glass/Epoxy composite, dimensionless stress intensity factor is affected by cracked-resin layer thickness.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.6
• /
• pp.1382-1387
• /
• 1994

A model is constructed to evaluate the stress intensity factors for a center crack subjected to anti-symmetric loading in a layered material. A Fredholm integral equation is derived using the Fourier integral transform method. The integral equation is numerically analyzed to evaluate the effects of stress intensity factor on the shear modulus, Poisson's ratio and crack length to layer thickness. In case of the isotropic homogeneous material, the values of stress intensity factor derived in the present study agree with the previous solutions.

• The Journal of the Acoustical Society of Korea
• /
• v.27 no.8
• /
• pp.411-417
• /
• 2008

An alternative formulation of the Helmholtz integral equation derived to express the pressure field explicitly in terms of the velocity vector of a radiating surface is used to solve acoustic radiation and fluid/structure interaction problems. This formulation, derived for arbitrary sources, is similar in form to the Rayleigh's formula for planar sources. Because the surface pressure field is expressed explicitly as a surface integral of the surface velocity, which can be implemented numerically using standard Gaussian quadratures, there is no need to use BEM to solve a set of simultaneous equations for the surface pressure at the discretized nodes. Furthermore the non-uniqueness problem inherent in methods based on Helmholtz integral equation is avoided. Validation of this formulation is demonstrated for some simple geometries.

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