• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.12 no.3
• /
• pp.485-494
• /
• 1999

본 연구에서는 요소를 사용하지 않고 절점들만을 이용하여 해석이 가능한 새로운 수치해석기법인 EFG(Element-Free Galerkin)법을 사용하여 임의의 균열의 성장과정을 해석할 수 있는 효율적인 알고리즘을 개발하고, 이를 바탕으로 균열의 성장방향과 경로를 정확히 추정하여 일련의 균열진전해석을 수행할 수 있는 프로그램을 개발하였다. 균열해석에 있어서는 균열선단의 특이성과 균열면의 분연속성을 수치적으로 반영할 수 있는 기법을 도입하여 균열을 모형화하였으며, 선형탄성파괴역학이론에 근거하여 균열해석과정을 정식화하였다. 또한, EFG 형상함수가 kronecker delta 조건을 만족시키지 못함으로써 발생하는 필수경계조건의 처리문제를 penalty법을 이용하여 해결하였다. 개발된 균열진전해석 알고리즘을 정지상태와 성장하는 상태에 있는 모드 Ⅰ, 모드 Ⅱ 및 혼합모드상태의 대표적인 균열문제들에 적용하여 응력확대계수와 균열성장방향 및 균열의 성장경로를 추정하고 이를 이론적·실험적 결과들과 비교함으로써 그 정확성과 효율성을 검증하였다.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.36 no.12
• /
• pp.1521-1528
• /
• 2012

The present paper provides the elastic stress intensity factors (SIFs) and the crack opening displacements (CODs) of a thick-walled pipe with a slanted axial through-wall crack. For estimating these elastic fracture mechanics parameters, systematic three-dimensional elastic finite element (FE) analyses were performed by considering geometric variables, i.e., thickness of pipe, reference crack length, and crack length ratio, affecting the SIFs and CODs. As for loading condition, the internal pressure was considered. Based on the FE results, the SIFs and CODs of slanted axial through-wall cracks in a thickwalled pipe along the crack front and the wall thickness were calculated. In particular, to calculate the SIFs of a thick-walled pipe with a slanted axial through-wall crack from those of a thick-walled pipe with an idealized axial through-wall crack, a slant correction factor representing the effect of the slant crack on the SIFs was proposed.

• Computational Structural Engineering
• /
• v.7 no.3
• /
• pp.167-176
• /
• 1994

A local remeshing algorithm using Delaunay property is developed for the analysis on the phenomenon of quasi-static crack propagation, which is a typical problem of accompanying constantly varying geometry. The algorithm performs both remeshing and refinement. The use of M-integral is demonstrated to simulate crack propagation under mixed mode with the edge spalling problem.

• Journal of Ocean Engineering and Technology
• /
• v.13 no.1 s.31
• /
• pp.39-50
• /
• 1999

특이 가중함수로 표현된 shepard interpolant와 일관조건을 사용하여 무요소법 형성함수를 도출하였다. 따라서 통상의 EFGM(Element Free Galerkin Method)과는 달리 변위로 주어지는 경계조건을 자연스럽게 부과할 수 있다. 수치계산 예로서 외팔보 문제를 다루었는데 보이론과 비교하여 매우 잘 맞는 결과를 보여주고, 유한요소법과의 결합도 자연스럽게 이루어짐을 보인다. 또 penny-shaped 균열을 다루는데, 응력확대계수는 균열 표면의 변위로부처 직접 계산하여 해석해와 비교한다.

• Computational Structural Engineering
• /
• v.8 no.4
• /
• pp.137-148
• /
• 1995

A new path independent contour integral formulus for the distinct calculation of mode I stress intensity factors in two dimensional linear elastic fracture mechanics problems is presented. This method is based on p-convergence concepts and can be easily appended to existing finite element computer codes. In this study, the stress state at crack tip has been investigated and the path independence of J-integral values has been tested with respect to different contours expressed by normalized distance apart from the crack tip. Numerical results by p-convergence for the problems such as centrally cracked panels, single and double edged cracks in rectangular panels have been compared with those by the conventional h-convergence. The comparison demonstrates the accuracy and stability of the proposed method.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.15 no.3
• /
• pp.433-440
• /
• 2002

A branched crack in a semi-infinite plate under uniform tension and bending moment is considered in this study By using the superposition, the stress and moment intensity factors for the branched crack subjected to uniform tension and bending moment we evaluated. The stress intensity factors we obtained by using the finite element method and the J-based mutual integral. The moment intensity factors are calculated by extrapolating the values of the moment new the crack tip. Numerical results lot the normalized stress and moment Intensity factors we shown as functions of the ratio of branched crack length to main crack length and the branching angle.

• Tunnel and Underground Space
• /
• v.25 no.4
• /
• pp.359-372
• /
• 2015

Rock masses include various rock discontinuities such as faults, joints, and bedding planes. These discontinuities appear as complex structures in geometry. In this study, growth patterns of fractures between two stepping shear fracture tips are numerically modeled using PFC2D (Particle Flow Code). The numerical model showed not only incipient growth of fractures at the tips of preexisting fractures but also subsequent growth of the new fractures. It is observed from all of the experiments that the incipient fractures are tensile cracks developed at $30{\sim}57^{\circ}$ to the preexisting fractures and the subsequent growth of these fractures were at low angles to the preexisting fractures this study.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.12 no.4
• /
• pp.691-700
• /
• 1999

본 연구에서는 요소를 사용하지 않는 새로운 해석방법인 EFG(Element-Free Galerkin)법을 사용하여 복수의 초기균열을 지닌 강재가 반복피로하중을 받는 경우 균열들이 점진적으로 성장하여 부재가 파단에 이르는 과정을 해석적으로 규명하였다. 이를 위하여 본 연구에서는 일반적인 피로균열성장법칙을 EFG법을 이용한 균열해석 알고리즘에 적용하여 복수의 균열들이 각각의 응력상태에 따라 차별적으로 성장해 나가는 과정을 해석할 수 있는 알고리즘을 도입하고 이를 바탕으로 다양한 하중상태하에서 복수의 균열들의 성장경로를 추정함과 동시에 이에 따른 잔존수명을 산정할 수 있는 기법을 제시하였다. 본 연구에서 제안된 해석방법을 피로균열 발생빈도가 큰 몇가지의 강부재 형태에 적용해 본 결과 다수균열 함유 부재의 피로균열 성장거동과 균열들의 피로수명을 성공적으로 예측할 수 있었다.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.36 no.5
• /
• pp.295-305
• /
• 2023

This study introduces a smoothed finite-element implementation into the phase-field framework. In recent years, the phase-field method has recieved considerable attention in crack initiation and propagation since the method needs no further treatment to express the crack growth path. In the phase-field method, high strain-energy accuracy is needed to capture the complex crack growth path; thus, it is obtained in the framework of the smoothed finite-element method. The salient feature of the smoothed finite-element method is that the finite element cells are divided into sub-cells and each sub-cell is rebuilt as a smoothing domain where smoothed strain energy is calculated. An adaptive quadtree refinement is also employed in the present framework to avoid the computational burden. Numerical experiments are performed to investigate the performance of the proposed approach, compared with that of the finite-element method and the reference solutions.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.30 no.11 s.254
• /
• pp.1455-1460
• /
• 2006

The relationship between fatigue crack growth behavior and cyclic crack tip opening displacement is studied. An elastic-plastic finite element analysis (FEA) is performed to examine the growth behavior of fatigue crack, where the contact elements are used in the mesh of the crack tip area. We investigate the relationship between the reversed plastic zone size and the changes of the cyclic crack tip opening displacement along the crack growth. We investigate the effect of the element size when predict fatigue crack opening behavior using the cyclic crack tip opening displacement obtained from FEA. The cyclic crack tip opening displacement is related to fatigue crack opening behavior.