• Structural Engineering and Mechanics
• /
• v.50 no.5
• /
• pp.665-677
• /
• 2014

A numerical method, using a special mixed finite element associated with the virtual crack extension technique, has been developed to evaluate the energy release rate for kinking cracks. The element is two dimensional 7-node mixed finite element with 5 displacement nodes and 2 stress nodes. The mixed finite element ensures the continuity of stress and displacement vectors on the coherent part and the free edge effect. This element has been formulated starting from a parent element in a natural plane with the aim to model different types of cracks with various orientations. Example problems with kinking cracks in a homogeneous material and bimaterial are presented to assess the computational accuracies.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.27 no.7
• /
• pp.1159-1169
• /
• 2003

This paper investigates the effects of T-stress and plastic hardening mismatch on the interfacial crack-tip stress field via finite element analyses. Plane strain elastic-plastic crack-tip fields are modeled with both MBL formulation and a full SEC specimen under pure bending. Modified Prandtl slip line fields illustrate the effects of T-stress on crack-tip constraint in homogeneous material. Compressive T-stress substantially reduces the interfacial crack-tip constraint, but increases the J-contribution by lower hardening material, J$\_$L/. For bimaterials with two elastic-plastic materials, increasing plastic hardening mismatch increases both crack-tip stress constraint in the lower hardening material and J$\_$L/. The fracture toughness for bimaterial joints would consequently be much lower than that of lower hardening homogeneous material. The implication of unbalanced J-integral in bimaterials is also discussed.

• Computational Structural Engineering
• /
• v.8 no.1
• /
• pp.173-186
• /
• 1995

It is well known that structures constructed by bonding two or more materials and then subjected to temperature change experience thermal stress. This stress results from thermal expansion mismatch of materials. The present paper derives formulas for the stresses in a bimaterial axisymmetric disk which is subjected to a uniform temperature change. First, an approximate solution following strength-of-materials principles is developed. However, the strength-of-materials solution has difficulty in predicting both the peak value of interfacial stresses and its associated distribution. Next, a solution consistent with the theory of elasticity is developed by way of an eigenfunction expansion approach. The eigenfunction analysis is compared with finite element stress analysis results for a specific numerical example. Finite element analysis results show that the interfacial stresses are adequately predicted by eigenfunction solution. Therefore, the method developed in this paper will be useful in determination of the interfacial stress state.

• Proceedings of the KSME Conference
• /
• 2001.06a
• /
• pp.102-107
• /
• 2001

Interfacial cracks between an isotropic and orthotropic material, subjected to static far field tensile loading are analyzed using the technique of photoelasticity. The fracture parameters are extracted from the full-field isochromatic data and the same are compared with that obtained using boundary collocation method. Dynamic Photoelasticity combined with high-speed digital photography is employed for capturing the isochromatics in the case of propagating interfacial cracks. The normalized stress intensity factors for static crack is greater when $\alpha=90^{\circ}C$ (fibers perpendicular to the interface) than when $\alpha=0^{\circ}C$ (fiber parallel to the interface) and those when $\alpha=90^{\circ}C$ are similar to ones of isotropic material. The dynamic stress intensity factors for interfacial propagating crack are greater when $\alpha=0^{\circ}C$ than $\alpha=90^{\circ}C$. The relationship between complex dynamic stress intensity factor $|K_D|$ and crack speed C is similar to that for isotropic homogeneous materials, the rate of increase of energy release rate G or $|K_D|$ with crack speed is not as drastic as that reported for homogeneous materials.

• Journal of the Korean Society of Safety
• /
• v.11 no.2
• /
• pp.9-15
• /
• 1996

The characteristics on the extension of the CED(Crack energy density) concept to the interface kinked crack problems in a dissimilar material are examined. Each mode contributions of CED are found by symmetric and antisymmetric components and domain independent integrals. Finite element calculation is carried out to simulate the Interface kinked crack growth on bimaterial. The focus is the establishment of fracture criterion with CED and finding the orientation of crack extension. From the results, a prediction about the extension behavior of an interface kinked crack can be done. And we show that CED can be a parameter to indicate fracture criterion at an Interface kinked crack.

• Journal of Mechanical Science and Technology
• /
• v.16 no.2
• /
• pp.165-174
• /
• 2002

In the present study, interfacial cracks between an isotropic and orthotropic material, subjected to static far field tensile loading are analyzed using the technique of photoelasticity. The fracture parameters are extracted from the full-field isochromatic data and the same are compared with that obtained using boundary collocation method. Dynamic photoelasticity combined with high-speed digital photography is employed for capturing the isochromatics in the case of propagating interfacial cracks. The normalized stress intensity factors for static cracks are greate. when ${\alpha}$: 90$^{\circ}$(fibers perpendicular to the interface) than when ${\alpha}$=0$^{\circ}$(fibers parallel to the interface), and those when ${\alpha}$=90$^{\circ}$are similar to ones of isotropic material. The dynamic stress intensity factors for interfacial propagating cracks are greater when ${\alpha}$=0$^{\circ}$ than ${\alpha}$=90$^{\circ}$. For the velocity ranges (0.1 < C/C$\sub$s1/<0.7) observed in this study, the complex dynamic stress intensity factor │K$\sub$D/│increases with crack speed c, however, the rate of increase of │K$\sub$D/│with crack speed is not as drastic as that reported for homogeneous materials.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.29 no.2 s.233
• /
• pp.293-302
• /
• 2005

In this paper, the thick-walled laminated, orthotropic as well as bimaterial, composite hollow spheres under impact pressure are analyzed in detail by using the semi-discrete finite element method with the Houbolt time-integration scheme which results in unconditionally stable transient numerical results. Numerical results are obtained by using the self-constructed spherically symmetric (one-dimensional) and axially symmetric (two-dimensional) finite element programs, and compared with the previous solutions by other researchers, being shown some of which are incorrect. The finite element package Nastran is also adopted for numerical comparison.

• Structural Engineering and Mechanics
• /
• v.66 no.5
• /
• pp.649-663
• /
• 2018

The scaled boundary finite element method is extended to evaluate the dynamic stress intensity factors and T-stress with a numerical procedure based on the improved continued-fraction. The improved continued-fraction approach for the dynamic stiffness matrix is introduced to represent the inertial effect at high frequencies, which leads to numerically better conditioned matrices. After separating the singular stress term from other high order terms, the internal displacements can be obtained by numerical integration and no mesh refinement is needed around the crack tip. The condition numbers of coefficient matrix of the improved method are much smaller than that of the original method, which shows that the improved algorithm can obtain well-conditioned coefficient matrices, and the efficiency of the solution process and its stability can be significantly improved. Several numerical examples are presented to demonstrate the increased robustness and efficiency of the proposed method in both homogeneous and bimaterial crack problems.

• Journal of the Korean Society for Precision Engineering
• /
• v.9 no.2
• /
• pp.108-115
• /
• 1992

이종재료의 접합면에 수직으로 존재하는 크랙에 대하여 경계요소 해석을 수행하여, 그 결과 실용가능한 수치 근사해을 얻을 수 있었다. 크랙을 정확히 모델링하기 위하여 크랙표면을 분리영역으로 하는 영역분할법을 채택하였으며, 해의 정확성을 향상시키기 위하여 등매개 2차요소로의 경계분할과 함께 크랙선단에서 표면력의 특이성을 나타내도록 하였다. 응력확대계수는 크랙표면상 절점의 상대변위를 이용하여 결정하였다. 또한 이종 재료내 크랙에 대하여 응력확대계수를 간단히 구할 수 있는 간편해석법을 제안하고 이의 적용 가능한 범위를 제 시하였다.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.21 no.3
• /
• pp.414-422
• /
• 1997

The characteristics on the extension of the CED(Crack Energy Density) concept to the interface kinked crack problems in a dissimilar are examined. Each mode contributions of CED are found by symmetric and antisymmetric conponents and domain independent integrals. Finite element calculation is carried out to simulate the interface kinked crack growth on a bimaterial. The focus is the establishment of fracture criterion with CED and finding the orientation of crack extension. From the results, a prediction about the extension behavior of an interface kinked crack can be done. And we show that CED can be a parameter to indicate fracture criterion at an interface kinked crack.