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

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.2 no.1
• /
• pp.11-17
• /
• 1990

Dirichlet boundary value problems originated from unsteady deep water wave propagation are transformed to Boundary Intergral Equation Methods by use of a free surface Green's function and the integral equations are discretized by a cubic spline element method. In order to enhance the stability of the numerical model based on the derived Fredholm integral equation of 1 st kind, the method by Hsiao and MacCamy (1973) is employed. The numerical model is tested against exact solutions for two cases and the model shows very good accuracy.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.05a
• /
• pp.15-28
• /
• 1998

As an alternative for solving the incompressible Navier-Stokes equations, we present a vorticity-based integro-differential formulation for vorticity, velocity and pressure variables. One of the most difficult problems encountered in the vorticity-based methods is the introduction of the proper value-value of vorticity or vorticity flux at the solid surface. A practical computational technique toward solving this problem is presented in connection with the coupling between the vorticity and the pressure boundary conditions. Numerical schemes based on an iterative procedure are employed to solve the governing equations with the boundary conditions for the three variables. A finite volume method is implemented to integrate the vorticity transport equation with the dynamic vorticity boundary condition . The velocity field is obtained by using the Biot-Savart integral derived from the mathematical vector identity. Green's scalar identity is used to solve the total pressure in an integral approach similar to the surface panel methods which have been well-established for potential flow analysis. The calculated results with the present mettled for two test problems are compared with data from the literature in order for its validation. The first test problem is one for the two-dimensional square cavity flow driven by shear on the top lid. Two cases are considered here: (i) one driven both by the specified non-uniform shear on the top lid and by the specified body forces acting through the cavity region, for which we find the exact solution, and (ii) one of the classical type (i.e., driven only by uniform shear). Secondly, the present mettled is applied to deal with the early development of the flow around an impulsively started circular cylinder.

• Kyungpook Mathematical Journal
• /
• v.55 no.4
• /
• pp.997-1030
• /
• 2015

This paper is concerned with integral type boundary value problems of second order singular differential equations with quasi-Laplacian on whole lines. Sufficient conditions to guarantee the existence and non-existence of positive solutions are established. The emphasis is put on the non-linear term $[{\Phi}({\rho}(t)x^{\prime}(t))]^{\prime}$ involved with the nonnegative singular function and the singular nonlinearity term f in differential equations. Two examples are given to illustrate the main results.

• Nonlinear Functional Analysis and Applications
• /
• v.28 no.2
• /
• pp.537-555
• /
• 2023

In this paper, we consider two types of fractional boundary value problems, one of them is an implicit type and the other will be an integro-differential type with nonlocal integral multi-point boundary conditions in the frame of generalized Hilfer fractional derivatives. The existence and uniqueness results are acquired by applying Krasnoselskii's and Banach's fixed point theorems. Some various numerical examples are provided to illustrate and validate our results. Moreover, we get some results in the literature as a special case of our current results.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.22 no.3
• /
• pp.9-18
• /
• 1985

The numerical calculation for solving boundary-value problem related to potential flows with a free surface is carried out by application of the localized finite element method. Only forced motion of 2-D body in infinitely deep fluid is considered, although this schemes is equally applicable to any first order time-harmonic problems of similar nature. The infinite domain of the fluid is separated into the inner flow field and the outer flow field with common inter-surface boundary. The finite element method is applied to obtain the solution in the inner flow field and the Green functions are utilized to represent the solution in the outer flow field. At the inter-surface boundary, the continuity of the value of potential and the normal derivative of the potential(i.e. matching condition) is conserved. The present method has better computational efficiency than the previous LFEM and the integral equation method of Frank. This enhanced computational efficiency is presumably due to the fact that the present method gives a symmetric coefficient matrix and requires less computational time in calculating the influence coefficient matrix of Green function than the integral equation method. And the irregular frequency desen't exist because the uniqueness of the solution is assured by the such that the exact free surface condition is satisfied on the boundary of the localized finite element region(i.e. inner region). As an example of the above method, the hydrodynamic forces for the circular cylinder and the rectangular cylinders are calculated. In the computed results, the small number of singularity distribution segments($3{\sim}6$) give good result relative to Ursell's and Vugts'.

• Journal of Aerospace System Engineering
• /
• v.1 no.1
• /
• pp.25-35
• /
• 2007

In this paper, the interlaminar crack problems of a fiber metal laminate (FML) under generalized plane deformation are studied using the theory of anisotropic elasticity. The crack is considered to be embedded in the matrix interlaminar region (including adhesive zone and resin rich zone) of the FML. Based on Fourier integral transformation and the stress matrix formulation, the current mixed boundary value problem is reduced to solving a system of Cauchy-type singular integral equations of the 1st kind. Within the theory of linear fracture mechanics, the stress intensity factors are defined on terms of the solutions of integral equations and numerical results are obtained for in-plane normal (mode I) crack surface loading. The effects of location and length of crack in the 3/2 and 2/1 ARALL, GLARE or CARE type FML's on the stress intensity factors are illustrated.

• Structural Engineering and Mechanics
• /
• v.35 no.2
• /
• pp.141-173
• /
• 2010

In this paper a boundary element method is developed for the general flexural-torsional buckling analysis of Timoshenko beams of arbitrarily shaped cross section. The beam is subjected to a compressive centrally applied concentrated axial load together with arbitrarily axial, transverse and torsional distributed loading, while its edges are restrained by the most general linear boundary conditions. The resulting boundary value problem, described by three coupled ordinary differential equations, is solved employing a boundary integral equation approach. All basic equations are formulated with respect to the principal shear axes coordinate system, which does not coincide with the principal bending one in a nonsymmetric cross section. To account for shear deformations, the concept of shear deformation coefficients is used. Six coupled boundary value problems are formulated with respect to the transverse displacements, to the angle of twist, to the primary warping function and to two stress functions and solved using the Analog Equation Method, a BEM based method. Several beams are analysed to illustrate the method and demonstrate its efficiency and wherever possible its accuracy. The range of applicability of the thin-walled theory and the significant influence of the boundary conditions and the shear deformation effect on the buckling load are investigated through examples with great practical interest.

• Journal of the Society of Naval Architects of Korea
• /
• v.56 no.1
• /
• pp.11-22
• /
• 2019

In this paper, a Rankine source method is applied and validated to analyze the hydrodynamic response of a three-dimensional floating structure in the frequency domain. The boundary value problems for radiation and diffraction problem are solved by using a desingularized indirect boundary integral equation method (DIBIEM). The DIBIEM is simpler and faster than conventional methods based on the numerical surface integration of Green's function because the singularities of Green's function are located outside of fluid regions. In case of floating structure with complex geometry, it is difficult to desingularize the singularities of Green's function consistently. Therefore a mixed approach is carried out in this study. The mixed approach is partially desingularized except singularities of the body. Wave drift loads are calculated by the middle-field formulation method that is mathematically simple and has fast convergence. In order to validate the accuracy of the developed program, various numerical simulations are carried out and these results are analyzed and compared with previously published calculations and experiments.

• Bulletin of the Korean Mathematical Society
• /
• v.53 no.2
• /
• pp.335-350
• /
• 2016

Based on the notion of $q_k$-derivative introduced by the authors in [17], we prove in this paper existence and uniqueness results for nonlinear second-order impulsive $q_k$-difference equations with anti-periodic boundary conditions. Two results are obtained by applying Banach's contraction mapping principle and Krasnoselskii's fixed point theorem. Some examples are presented to illustrate the results.