• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.34 no.7
• /
• pp.881-889
• /
• 2010

A volume integral equation method (VIEM) is introduced for solving the elastostatic problems related to an unbounded isotropic elastic solid; this solid is subjected to remote uniaxial tension, and it contains multiple interacting isotropic inclusions. The method is applied to two-dimensional problems involving long parallel cylindrical inclusions. A detailed analysis of the stress field at the interface between the matrix and the central inclusion is carried out; square and hexagonal packing of the inclusions are considered. The effects of the number of isotropic inclusions and different fiber volume fractions on the stress field at the interface between the matrix and the central inclusion are also investigated in detail. The accuracy and efficiency of the method are clarified by comparing the results obtained by analytical and finite element methods. The VIEM is shown to be very accurate and effective for investigating the local stresses in composites containing isotropic fibers.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.36 no.1
• /
• pp.59-71
• /
• 2012

A volume integral equation method (VIEM) is introduced for the solution of elastostatic problems in unbounded isotropic elastic solids containing multiple interacting isotropic or anisotropic diamond-shaped inclusions subject to remote uniaxial tension. The method is applied to two-dimensional problems involving long parallel diamond-shaped cylindrical inclusions. A detailed analysis of the stress field at the interface between the matrix and the central inclusion is carried out for square and hexagonal packing of the inclusions. The effects of the number of isotropic or anisotropic diamond-shaped inclusions and of the various fiber volume fractions for the circular inclusions circumscribing its respective diamond-shaped inclusion on the stress field at the interface between the matrix and the central inclusion are also investigated in detail. The accuracy and efficiency of the method are examined through comparison with results obtained using the finite element method.

• Composites Research
• /
• v.23 no.4
• /
• pp.7-13
• /
• 2010

A volume integral equation method (VIEM) is introduced for the solution of elastostatic problems in an unbounded isotropic elastic solids containing interacting multiple anisotropic inclusions subject to remote uniaxial tension. The method is applied to two-dimensional problems involving long parallel cylindrical inclusions. A detailed analysis of stress field at the interface between the matrix and the central inclusion is carried out for square and hexagonal packing of the inclusions. Effects of the number of anisotropic inclusions and various fiber volume fractions on the stress field at the interface between the matrix and the central inclusion are also investigated in detail. The accuracy of the method is validated by solving the single inclusion problem for which solutions are available in the literature.

• Computational Structural Engineering
• /
• v.9 no.4
• /
• pp.127-134
• /
• 1996

This paper is about a progressive fracture analysis of concrete by boundary element method. From both displacement boundary integral equation and traction boundary integral equation of solids with cracks, a boundary integral equation for crack problem is derived. For the analysis of progressive fracture of concrete, fracture process zone is modelled based on Dugdale-Barenblatt model with linear tension-softening curve. By using the boundary element modeling, the progressive fractures of concrete beam and compact-tension specimens with various loading conditions are analyzed and compared with experiments. The analysis results show that the technique in this paper can predict the maximum strength and the nonlinear behavior of concrete including post-peak behavior.

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

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.3 no.4
• /
• pp.223-234
• /
• 1991

In this paper, the harbor tranquility problem is analysed by an improved boundary integral equation method. The effect of the diffracted wave transformation induced by the breakwater and structures located at a harbor mouth is considered. Partial reflection concept is also used to consider energy dissipation effects. The present model is verified by comparing the results of the model for rectangular and semi-circular harbors with the analytic solutions. they show a reasonable agreement. Also the wave height distribution of the HUPO harbor computed using the present model agree well to those from the previous hydraulic model tests. It also shows a good agreement with the results from the time-dependent mild slope equation model.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.21 no.4
• /
• pp.375-382
• /
• 2008

The particular integral formulation for two(2D) and three(3D) dimensional inelastic transient dynamic stress analysis is presented. The elastostatic equation is used for the complementary solution. Using the concept of global shape function, the particular integrals for displacement and traction rates are obtained to approximate acceleration of the inhomogeneous equation. The Houbolt time integration scheme is used for the time-marching process. The Newton-Raphson algorithm for plastic multiplier is used to solve the system equation. Numerical results of four example problems are given to demonstrate the validity and accuracy of the present formulation.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.19 no.5
• /
• pp.502-510
• /
• 2007

The present study is to estimate the effect of diffracted wave fields inside a harbor, around harbor entrance and outer breakwater, when a navigation channel is dredged in the vicinity of the a harbor entrance. The wave field of the problem is considered to be two-dimensional plane and the configuration of the submarine pit on the sea bed is designated by a single rectangular type. The numerical simulation is performed by using the solution of the Greet function based on the boundary integral equation. The results of this study is illustrated by applying the normal incidence and partially reflecting boundaries.

• Geophysics and Geophysical Exploration
• /
• v.3 no.2
• /
• pp.61-69
• /
• 2000

The free-surface boundary conditions are persistent problem in elastic wave equation modeling by finite-difference method, which can be summarized with the degradation of the accuracy of the solution and limited stability range in Poisson's ratio. In this paper, we propose the mid-point averaging scheme as an alternative way of implementing the free-surface boundary conditions, and present the solution to Lamb's problem to verify our approach.

• Journal of Ocean Engineering and Technology
• /
• v.12 no.1
• /
• pp.99-106
• /
• 1998

3차원 파수조에서 완전 비선형파를 시뮬레이션하기 위하여 우선 랜킨 소스를 기저로한 적분방정식을 고차경계요소법을 이용하여 이산화하였다. 그리고 방사경계조건은 파흡수 비치와 포텐셜 스트레칭 기법을 이용하여 모델링하였으며, 비선형 자유표면과 경계조건식은 고차 예측 및 보정 기법을 이용하여 시간 적분하였다. 파흡수 비치는 파의 진행방향 특성에 따라 수조내에 다양하게 배치할 수 있으며 비칭서 흡수가 덜된 파는 수조의 길이 방향 끝단에서 포텐셜 스트레칭 기법에 의하여 반사없이 진행하도록 하였다. 수치실험 결과 일-에너지 보존법칙과 모멘텀-임펄스 보존 법칙이 만족됨으로써 본 수치기법의 효용성이 검증되었다.