• Journal of the Society of Naval Architects of Korea
• /
• v.41 no.4
• /
• pp.38-44
• /
• 2004

The motion response analysis of a submerged elliptic cylinder in waves is presented and the elliptic cylinder is a simplification of the section of submarine in this paper. The method is based on boundary integral method and two-dimensional 3 degree motions are calculated in regular harmonic waves. The fully nonlinear free surface boundary condition is assumed in an numerical domain and this solution is matched along an assumed boundary as a linear solution composed of transient Green function, The large amplitude motions of an elliptic cylinder are directly simulated and effects of wave frequency, wave amplitude and the distance from buoyancy center to gravity center are discussed.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.9
• /
• pp.2138-2144
• /
• 1993

An integral equation representation of cracks was presented which differs from well-known "dislocation layer" representation. In this new representation, the equations are written in terms of the displacement discontinuity across the crack surfaces rather than derivatives of the displacement-discontinuity. It was shown in that the new technique is well-suited to the treatment of kinked cracks. In the present paper, this integral equation representation is coupled to the direct boundary-element method for the treatment of finite bodies containing kinked cracks. The method is demonstrated for two-dimensional finite domains but extension to three-dimensional problems would appear to be possible. The resulting approach is shown to be simple, yet very accurate. accurate.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.11 no.1
• /
• pp.45-53
• /
• 1991

The underground structure, which has infinite or semi-infinite boundary conditions, is subjected by body forces and in-situ stresses. It also has stress concentration, which causes material nonlinear behavior, in the vicinity of the excavated surface. In this paper, some methods which can be used to transform domain integrals into boundary integrals are reviewed in order to analyze the effect of the body forces and the in-situ stresses. First, the domain integral of the body force is transformed into boundary integral by using the Galerkin tensor and divergence theorem. Second, it is transformed by writing the domain integral in cylindrical coordinates and using direct integration. The domain integral of the in-situ stress is transformed into boundary integral applying the direct integral method in cylindrical coordinates. The methodology is verified by comparing the results from the boundary element analysis with those of the finite element analysis. Coupling the above boundary elements with finite elements, the nonlinear behavior that occurs locally in the vicinity of the excavation is analyzed and the results are verified. Thus, it is concluded that the domain integrals of body forces and in-situ stresses could be performed effectively by transforming them into the boundary integrals, and the nonlinear behavior can be reasonably analyzed by coupled nonlinear finite element and boundary element method. The result of this research is expected to he used for the analysis of the underground structures in the effective manner.

• Journal of the Korean Society of Safety
• /
• v.14 no.1
• /
• pp.10-18
• /
• 1999

An integral equation representation of cracks was presented, which differs from well-known "dislocation layer" representation. In this new representation, the integral equation representation of cracks was developed and coupled to the direct boundary-element method for treatment of cracks in finite plane bodies. The method was developed for in-plane(mode I and II) loadings only. In this paper, the method is formulated and applied to various crack problems involving multiple and branch cracks in finite region. The results are compared to exact solutions where available and the method is shown to be very accurate despite of its simplicity.implicity.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.28 no.9
• /
• pp.1284-1296
• /
• 2004

In order to investigate effects of anisotropic fiber packing on stresses in composites, a Volume Integral Equation Method is applied to calculate the elastostatic field in an unbounded isotropic elastic medium containing multiple orthotropic inclusions subject to remote loading, and a Mixed Volume and Boundary Integral Equation Method is introduced for the solution of elastostatic problems in unbounded isotropic materials containing multiple anisotropic inclusions as well as one void under uniform remote loading. A detailed analysis of stress fields at the interface between the isotropic matrix and the central orthotropic inclusion is carried out for square, hexagonal and random packing of orthotropic cylindrical inclusions, respectively. Also, an analysis of stress fields at the interface between the isotropic matrix and the central orthotropic inclusion is carried out, when it is assumed that a void is replaced with one inclusion adjacent to the central inclusion of square, hexagonal and random packing of orthotropic cylindrical inclusions, respectively, due to manufacturing and/or service induced defects. The effects of random orthotropic fiber packing on stresses at the interface between the isotropic matrix and the central orthotropic inclusion are compared with the influences of square and hexagonal orthotropic fiber packing on stresses. Through the analysis of plane elastostatic problems in unbounded isotropic matrix with multiple orthotropic inclusions and one void, it will be established that these new methods are very accurate and effective for investigating effects of general anisotropic fiber packing on stresses in composites.

• Journal of computational fluids engineering
• /
• v.3 no.1
• /
• pp.11-21
• /
• 1998

A vorticity-based method for the numerical solution of the two-dimensional incompressible Navier-Stokes equations is presented. The governing equations for vorticity, velocity and pressure variables are expressed in an integro-differential form. The global coupling between the vorticity and the pressure boundary conditions is fully considered in an iterative procedure when numerical schemes are employed. The finite volume method of the second order TVD scheme 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. The 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 present formulation is validated by comparison with data from the literature for the two-dimensional cavity flow driven by shear in a square cavity. We take two types of the cavity now: (ⅰ) driven by non-uniform shear on top lid and body forces for which the exact solution exists, and (ⅱ) driven only by uniform shear (of the classical type).

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.30 no.1
• /
• pp.47-53
• /
• 2017

A boundary layer integral combined with a 1-D isentropic core flow model has been successfully used to determine heat transfer rate on the surface of a supersonic nozzle. However its accuracy is affected by the core flow condition which is used as a boundary condition for the integral calculation. Because flow behavior near a nozzle throat deviates from 1-D isentropic condition due to 2-D flow turning and interaction between core flow and boundary layer, accuracy of heat transfer calculation decreases at a nozzle throat. Therefore, CFD is adopted to deduce improved core flow condition and increase accuracy of boundary layer integral at nozzle throat in this research. Euler model and SST $k-{\omega}$ model is solved by CFD code and used as a boundary condition for boundary layer integral. Developed code is tested in the supersonic nozzle from the previous research and improvement in accuracy is observed, especially at nozzle throat and diverging section of the nozzle. Error between experimental result and calculation result reduced by 16% when a calculation is made based on the SST $k-{\omega}$ model. Method developed in this research is expected to be used in thermal design of the rocket nozzle.

• Advances in nano research
• /
• v.16 no.3
• /
• pp.265-272
• /
• 2024

In this paper, the problem of static bending of axially functionally graded (AFG) nanobeam is formulated with the local stress(Lσ)- and strain-driven(εD) two-phase local/nonlocal integral models (TPNIMs). The novelty of the present study aims to compare the size-effects of nonlocal integral models on bending deflections of AFG Euler-Bernoulli nano-beams. The integral relation between strain and nonlocal stress components based on two types nonlocal integral models is transformed unitedly and equivalently into differential form with constitutive boundary conditions. Purely LσD- and εD-NIMs would lead to ill-posed mathematical formulation, and Purely εD- and LσD-nonlocal differential models (NDM) may result in inconsistent size-dependent bending responses. The general differential quadrature method is applied to obtain the numerical results for bending deflection and moment of AFG nanobeam subjected to different boundary and loading conditions. The influence of AFG index, nonlocal models, and nonlocal parameters on the bending deflections of AFG Euler-Bernoulli nanobeams is investigated numerically. A consistent softening effects can be obtained for both LσD- and εD-TPNIMs. The results from current work may provide useful guidelines for designing and optimizing AFG Euler-Bernoulli beam based nano instruments.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1998.04a
• /
• pp.41-52
• /
• 1998

A volume integral equation method and a mixed volume and boundary integral equation method are presented for the solution of plane elastostatic problems in solids containing orthotropic inclusions and voids. The detailed analysis of the displacement and stress fields are developed for orthotropic cylindrical and elliptic-cylindrical inclusions and voids. The accuracy and effectiveness of the new methods are examined through comparison with results obtained from analytical and boundary integral equation methods. Through the analysis of plane elastostatic problems in unbounded isotropic matrix containing orthotropic inclusions and voids, it is established that these new methods are very accurate and effective for solving plane elastostatic and elastodynamic problems in unbounded solids containing general anisotropic inclusions and voids or cracks.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.11 no.6
• /
• pp.1001-1004
• /
• 1987

In this study mathematical properties of variational solution and solution of the boundary integral equation of the linear elasticity problem are studied. It is first reviewed that a variational solution for the three-dimensional linear elasticity problem exists in the Sobolev space [ $H^{1}$(.OMEGA.)]$^{3}$ and, then, it is shown that a unique solution of the boundary integral equation is identical to the variational solution in [ $H^{1}$(.OMEGA.)]$^{3}$. To represent the boundary integral equation, the Green's formula in the Sobolev space is utilized on the solution domain excluding a ball, with small radius .rho., centered at the point where the point load is applied. By letting .rho. tend to zero, it is shown that, for the linear elasticity problem, boundary integral equation is valid for the variational solution. From this fact, one can obtain a numerical approximatiion of the variational solution by the boundary element method even when the classical solution does not exist.exist.