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

In this paper, we further generalize the work of Lin and Abel to the case of the first and the second order derivatives of energy release rates for two-dimensional, multiply cracked systems. The direct integral expressions are presented for the energy release rates and their first and second order derivatives. The salient feature of this numerical method is that the energy release rates and their first and second order derivatives can be computed in a single analysis. It is demonstrated through a set of examples that the proposed method gives expectedly decreasing, but acceptably accurate results for the energy release rates and their first and second order derivatives. The computed errors were approximately 0.5% for the energy release rates, $3\sim5%$ for their first order derivatives and $10\sim20%$ for their second order derivatives for the mesh densities used in the examples. Potential applications of the present method include a universal size effect model and a probabilistic fracture analysis of cracked structures.

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 2005.06a
• /
• pp.337-340
• /
• 2005

The effect of residual stress on fatigue crack growth was investigated in terms of finite element analysis. Simulations were performed on a CT specimen in plane strain. An interface-cohesive element that accounts for damage accumulation due to fatigue along the notch direction has been used. Numerical results show that fatigue crack growth rate slows down when compressive residual stress field exists in front of the crack tip.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.13 no.2
• /
• pp.221-228
• /
• 1989

The objective of this paper is to develop a numerical technique for analyzing crack driving forces in multi-phase materials. The analysis was based on finite element method coupled with a virtual crack extension technique which is known as the most efficient tool in computational fracture mechanics analysis. The modified J-integral method, proposed by Miyamoto and Kikuchi for the analysis of dual-phase material was carried out by subtracting the J-values for contours surrounding each phase boundary from the J-values for overall contour. It was shown that the proposed numerical procedure, based on the modified J-integral coupled with a virtual crack extension technique, can be used as an effective numerical tool for determining crack driving forces in multi-phase materials.

• Tunnel and Underground Space
• /
• v.21 no.1
• /
• pp.66-81
• /
• 2011

In this study, it was aimed to develop a thermal-hydraulic-mechanical coupled fracture mechanics code that models a fracture initiation, propagation and failure of underground rock mass due to thermal and hydraulic loadings. The development was based on a 2D FRACOD (Shen & Stephasson, 1993), and newly developed T-M and H-M coupled analysis modules were implemented into it. T-M coupling in FRACOD employed a fictitious heat source and time-marching method, and explicit iteration method was used in H-M coupling. The validity of developed coupled modules was verified by the comparison with the analytical result, and its applicability to the fracture initiation and propagation behavior due to temperature changes and hydraulic fracturing was confirmed by test simulations.

• Journal of the Korean Society of Safety
• /
• v.11 no.4
• /
• pp.143-150
• /
• 1996

This is an explicit-Implicit, finite element analysis for linear as well as nonlinear hygrothermal stress problems. Additional features, such as moisture diffusion equation, crack element and virtual crack extension(VCE ) method for evaluating J-integral are implemented in this program. The Linear Elastic Fracture Mechanics(LEFM) Theory is employed to estimate the crack driving force under the transient condition for and existing crack. Pores in materials are assumed to be saturated with moisture in the liquid form at the room temperature, which may vaporize as the temperature increases. The vaporization effects on the crack driving force are also studied. The Ideal gas equation is employed to estimate the thermodynamic pressure due to vaporization at each time step after solving basic nodal values. A set of field equations governing the time dependent response of porous media are derived from balance laws based on the mixture theory Darcy's law Is assumed for the fluid flow through the porous media. Perzyna's viscoplastic model incorporating the Von-Mises yield criterion are implemented. The Green-Naghdi stress rate is used for the invariant of stress tensor under superposed rigid body motion. Isotropic elements are used for the spatial discretization and an iterative scheme based on the full newton-Raphson method is used for solving the nonlinear governing equations.

• Journal of the Korean Professional Engineers Association
• /
• v.20 no.3
• /
• pp.5-14
• /
• 1987

연속체의 해석에 있어서, 특별한 경우를 제외하고는, 구조물의 개략적인 거동을 파악해야 될 경우가 종종 있다. 이러한 요구에 부응하기 위해서 강체요소법(Rigid Element Method)이라 불리우는 새로운 해석법이 개발되었다. 강체요소법은 원래 평정연구실에서 벽식프리캐스트 철근콘크리트 구조물의 탄소성해석을 하기 위해서 개발된 해석법에 착안하여, 내수벽과 같은 연속체에 적용함으로서 시작된 수치해석법이다. 그 후 저자들은 도통쉘, 구형쉘 혹은 이들이 조합된 쉘구조물에 적용할 수 있도록 개발 확장하였다. 강체요소법의 기본개념은 연속체의 분해된 각 요소를 강체(rigid body)라고 가정하고, 각 요소들은 요소의 강성으로 치환된 가상스프링으로 서로 연결되어 있다고 가정하여, 이 가상스프링의 거동을 평가함으로서 전체구조물의 거동을 파악하는 해석법이다. 이때 요소의 주변에 취해진 스프링은 해석을 단순화하기 위해서 축력, 면내전단력 및 면외전단력만을 전달한다고 가정하고, 요소의 강체변위(자유도)는 요소내의 임의의 한 점에서 취하며, 이 점에서의 강체변위(rigid displacements)는 요소의 주변에 취해진 스프링을 통하여 다른 요소로 전달된다. 상기와 같은 강체요소법의 개념을 연속체의 탄성 및 탄소성해석에 적용하면, 해석적 개념이 단순할 뿐만 아니라 구조물 전체의 자유도수를 대폭 줄여 컴퓨터 계산시간을 절약할 수 있는 잇점이 있고, 거시적인 모델(macroscopic modeling)과 미시적인 모델 (microscopic modeling)의 중간적인 성격을 가지기 때문에 구조물의 파괴상황에 대해서도 그 개략을 파악할 수 있다. 본 논문에서는 강체요소법을 보다 일반화된 해석법으로 개발, 확장하기 위해서 종전에 단층스프링시스템(single-layer spring system)으로 해석이 어려웠던 문제점들을 보완한 복층프링시스템(double-layer spring system)을 사용함으로서 휨, 비틀림의 효과를 파악할 수 있는 이론적 개념을 적용한 새로운 구요소, 원통요소 및 평면요소를 개발하고, 이러한 강체요소들의 적합매트릭스의 유도 및 해석저긴 방법을 정식화하였다. 또 휨, 비틀림 및 전단력의 효과를 고려한 사각형원통요소 및 능형원 통요소를 이용하여 원통쉘의 탄성 및 탄소성해석할 수 있는 프로그램을 개발하고, 이 프로그램으로 캔틸레버로된 연속형철근콘크리트 원통쉘의 탄성 및 탄소성해석에 적용하여 구조물의 거동에 관한 수치해석의 결과, 즉 내력의 분포, 균열의 진전, 파괴의 상황 및 변형의 상태 등을 파악해 보았다.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.38 no.7
• /
• pp.727-734
• /
• 2014

This paper describes multi-scale based ductile fracture simulation using finite element (FE) damage analysis. The maximum and crack initiation loads of cracked components were predicted using proposed virtual testing method. To apply the local approach criteria for ductile fracture, stress-modified fracture strain model was adopted as the damage criteria with modified calibration technique that only requires tensile and fracture toughness test data. Element-size-dependent critical damage model is also introduced to apply the proposed ductile fracture simulation to large-scale components. The results of the simulation were compared with those of the tests on SA333 Gr. 6 full-scale pipes at $288^{\circ}C$, performed by the Battelle Memorial Institute.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.40 no.8
• /
• pp.693-701
• /
• 2016

In this study, notch defects are evaluated using fracture mechanics. To understand the effects of notch defects, FE analysis is conducted to predict the limit load and J-integral for middle-cracked and single-edge cracked plates with various sizes of notch under tension and bending. As the radius of the notch increases, the energy release rate also increases, although the limit load remains constant. The values of fracture toughness($J_{IC}$) of SM490A are determined for various notch radii through FE simulation instead of conducting an experiment. As the radius of the notch increases, the energy release rate also increases, together with a more significant increase in fracture toughness. To conclude, as the notch radius increases, the resistance to crack propagation also increases.