Determination of Energy Release Rate of Penny-shaped Interface Crack on Bimaterial Cylinder

The mixed mode problem (I and II) of a peny-shaped interface cracks in remote tension loading on a bi-material cylinder is studied using finite element method. The energy release rates for the tip of the crack in the interface were calibrated for several different moduli combinations and crack ratios using the modified crack closure integral technique and J-integral method, with numerical results obtained from a commercial finite element program. Numerical results show that non-dimensional value of$\sqrt{G_{II}E^*}/\sqrt[p]{\pi a}$ increases as the crack size or moduli ratio increases. Meanwhile, non-dimensional value of$\sqrt{G_{I}E^*}/\sqrt[p]{\pi a}$ decreases as the moduli ratio increases, but above the moduli ratio of 3 its value decreases then increases again as the crack size increases. Reliability of the numerical analysis in this study was acquired with comparison to an analytical solution for the peny-shaped interface crack in an infinite medium.

동전모양 균열이 존재하는 이상복합체의 에너지해방율 산정

동전모양의 균열이 이상복합 실린더 계면에 존재하는 혼합모드 조건(I, II)에 대해 유한요소법을 사용하여 에너지해방율을 구하였다. 두재료의 탄성비와 노치율을 변화시켜 상업용 FEM 프로그램인 ABAQUS로부터 얻은 결과를 가상 균열법과 J 적분법에 적용하였으며 에너지해방율을 구하여 무차원함수로 표현하였다. 모드 II의 무차원 에너지해방율($\sqrt{G_{II}E^*}/\sqrt[p]{\pi a}$)은 균열길이와 탄성비가 증가되면서 그 값이 증가됨을 알수 있었다. 반면, 모드 I의 무차원 에너지해방율($\sqrt{G_{I}E^*}/\sqrt[p]{\pi a}$)은 탄성비가 증가하면서 그 값이 감소하며, 두재료의 탄성비가 3 이상인 경우에 균열길이가 증가되면서 무차원 에너지해방율이 감소하다가 다시 증가하게 나타났다. 또한 수치해석된 결과치를 무한판 실린더의 응력확대계수에 대한 정해와 비교하여 본 해석의 신뢰성을 확보하였다.