• Journal of the Korean Society for Precision Engineering
• /
• v.14 no.6
• /
• pp.64-73
• /
• 1997

솔리드 모델러, 자동요소분할 기법, 4면체 특이요소, 응력확대계수의 해석 기능을 통합하여, 3차원 균열의 응력확대계수를 효율적으로 해석할 수 있는 시스템을 개발하였다. 균열을 포함하는 기하모델을 CAD 시스템을 이용하여 정의하고, 경계조건과 재료 물성치 및 절점밀도 분포를 기하모델에 직접 지정함으로써, 퍼지이론 에 의한 절점발생과 데로우니 삼각화법에 의한 요소가 자동으로 생성된다. 특히, 균열 근방에는 4면체 2차 특이요소를 생성시켰으며, 유한요소 해석을 위한 입력 데이터가 자동으로 작성되어 해석코드에 의한 응력 해석이 수행된다. 해석 후, 출력되는 변위를 이용하여 변위외삽법에 의한 응력확대계수가 자동적으로 계산되어 진다. 본 시스템의 효용성을 확인하기 위해, 인장력을 받는 평판내의 표면균열에 대해 해석하여 보았다.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2010.04a
• /
• pp.314-317
• /
• 2010

본 논문에서는 p-수렴 천이요소의 층별이론을 사용하여 균열을 갖는 적층판의 팻취보강 효과에 대해 알아보았다. 이 모델의 면내 변위는 각 구간별로 연속인 선형변화로 가정하였고, 두께방향으로의 면외 변위는 일정한 상수로 가정하여 적용하였다. 변위장의 정의는 적분형 르장드로 다항식을 적용하였다. 또한 에너지 방출률법과 VCCT 방법을 사용하여 응력확대계수를 산출하였다. 수치적분은 별도의 외삽법 없이 각층별의 절점에서 방생하는 적분 값을 바로 얻을 수 있는 가우스-로바토 적분법을 사용하여 계산하였으며, 수치예제를 통해 제안된 모델의 정확도와 기존의 3차원 고체요소를 사용한 것보다 동일한 정확도를 얻기 위해 휠씬 적은 요소 및 자유도가 사용됨을 알 수 있었다.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.5A
• /
• pp.649-656
• /
• 2008

The enhancement of the service life of damaged or cracked structures is a major issue for researchers and engineers. The hierarchic void element based on the integrals of Legendre polynomials is used to characterize the fracture behaviour of unpatched crack as well as repaired crack with bonded composite patches by computing the stress intensity factors and stress contours at the crack tip. Since the equivalent single layer approach is adopted in this study, the proposed element is necessary to represent a discontinuous crack part as a continuum body with zero stiffness. Thus the aspect ratio of this element to represent the crack should be extremely slender. The sensitivity of numerical solution with respect to energy release rate, displacement and stress has been tested to show the robustness of zero stiffness element as the aspect ratio is increased up to 2000. The stiffness derivative method and displacement extrapolation method have been applied to calculate the stress intensity factors of Mode I problem. It is noted that the proposed hierarchical void element can be one of alternatives to analyze the patched crack problems.

• Computational Structural Engineering
• /
• v.9 no.2
• /
• pp.133-142
• /
• 1996

In this paper, two different techniques for mixed-mode type engineering fracture mechanics are investigated to estimate the stress intensity factors by using p-version finite element model. These two techniques are displacement extrapolation with COD and CSD method and J-integral with decomposition method. By decomposing the displacement field obtained from p-version of finite element analysis into symmetric and antisymmetric displacement fields with respect to the crack line, Mode-I and Mode-II stress intensity factors can be determined using aforementioned techniques. The example problems for validating the proposed techniques are centrally and centrally oblique cracked panels under tension. The numerical results associated with the variation of oblique angle and the ratio of crack length and panel width (a /W ratio) are compared with those by theoretical values and empirical solutions in literatures. Very good agreements with the existing solutions are shown.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.23 no.2
• /
• pp.191-198
• /
• 2010

This study investigated that the stress reduction of single edge cracked plates with patch repairs according to different type of patching such as material, size and thickness of patch and adhesive as well as single sided or double sided patches. As a numerical tool, the p-convergent partial layerwise model has been employed. The proposed model is formulated by assuming piecewise linear variation of in-plane displacement and a constant value of out-of-plane displacements across thickness. The integrals of Legendre polynomials are chosen to define displacement fields and Gauss-Lobatto numerical integration is implemented in order to directly obtain maximum values occurred at the nodal points of each layer without other extrapolation techniques. Also, total strain energy release rate method is adopted to obtain stress intensity factors. Numerical examples are presented not only to demonstrate the stress reduction effect in terms of non-dimensional stress intensity factor and deflection with respect to different type of patch repairs, but also the accuracy of proposed model.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.12
• /
• pp.2226-2240
• /
• 1992

This paper studies the applicability of the method of caustics into anisotropic materials under mode I and mixed static loading conditions and introduces the procedure to obtain stress intensity factors(S.I.F.) in anisotropic materials by the method of caustics. The mapping equations for initial and caustic curves in anisotropic materials were introduced and their computer graphical images were compared to the experimental ones to check the validity of the mapping equations proposed in this paper. The agreement between them was found to satisfactory. Two kinds of method to determine S.I.F. in anisotropic materials by the method of caustics were proposed in this paper and applied into the orthotropic materials under various loading conditions. In the case of mode I loading condition, the S.I.F.'s obtained by this paper's methods were found to be quite similar to the results by other method, boundary element method(B.E.M) and in the case of mixed loading condtion, the S.I.F's by this paper and B.E.M. showed a little differences(2.2-24.4%) with respect to the slanted angle of crack.

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.15 no.6
• /
• pp.547-557
• /
• 2013

The mode II fracture toughness and strength due to shear stress are important parameters in the stability of caprock and injection zone with application to geological sequestration of carbon dioxide. In this research, a short beam compression test has been used to determine the shear strength and the mode II fracture toughness for Coconino sandstone. The average value of the shear strength and mode II fracture toughness are estimated to be 23.53 MPa and 1.58 MPa${\surd}$m respectively. The stress intensity factor is suggested by finite element analysis using the displacement extrapolation method. The effect of biaxial stress and water saturation on the fracture toughness has also been investigated. The fracture toughness increases with confining stresses, but decreases by 11.4% in fully saturated condition.