• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.16 no.2
• /
• pp.207-216
• /
• 2003

In this paper, the finite element method for the elastodynamic axisymmetric fracture analysis is presented in matrix form through the application of the Galerkin method to the time integral equations of motion with no inertia forces. Isoparametric quadratic quadrilateral element and triangular crack tip singular elements with one-quarter node are used in the mesh division of the finite element model. To show the validity and accuracy of the proposed method, the infinite elastic medium with the penny shaped crack is solved first and compared with the analytical solution and the numerical results by the finite difference method and the boundary element method existing in the published literatures, and then the dynamic stress intensity factors of solid and hollow cylinders of finite dimensions haying penny-shaped cracks and internal and external circumferential tracks are computed in detail.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.37 no.4
• /
• pp.491-496
• /
• 2013

Finite element simulations of tensile tests were performed to determine the equivalent stress - equivalent plastic strain curves, critical equivalent stresses, and critical equivalent plastic strains. Then, the curves were used as inputs to finite element simulations of fracture toughness tests to determine the plane strain fracture toughness. The critical COD was taken as the COD when the equivalent plastic strain at the crack tip reached a critical value, and it was used as a crack growth criterion. The relationship between the critical COD and the critical equivalent plastic strain or the reduction of area was found. The relationship between the plane strain fracture toughness and the product of the critical equivalent stress and the critical equivalent plastic strain was also found.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.13 no.4
• /
• pp.636-642
• /
• 1989

In this paper, 3-dimensional fracture phenomena in the local area near the notch tip between the surface and the center of a composite bar which is consisted of two different materials were investigated by using embedded dyeing grids with the pitch of 50.8.mu. It was confirmed that the maximum strain near the notch tip occurred not at the notch front but at the curvature beginning point of the notch which is in the interior of z/W = 0.3 - 0.35 from the surface and a notch was initiated from that point. And that the strain near the notch front between the center and z/W = 0.65 - 0.7 toward the surface was in a state of compression.

• Journal of the Society of Naval Architects of Korea
• /
• v.35 no.1
• /
• pp.61-71
• /
• 1998

This paper describes a dynamic fracture behaviors of structural elements under elastic or elasto-plastic stress waves in two dimensional space. The governing equation of this problem has the type of hyperbolic partial differential equation, which consists of the equation of motions and incremental elasto-plastic constitutive equations. To solve this problem we introduce Zwas' method which is based on the finite difference method. Additionally, in order to deal with the dynamic behavior of elasto-plastic problems, an elasto-plastic loading path in the stress space is proposed to model the plastic yield phenomenon. Based on the result of this computation, the dynamic stress intensity factor at the crack tip of an elastic material is calculated, and the time history of a plastic zone of a elasto-plastic material is to be shown.

• Journal of Ocean Engineering and Technology
• /
• v.15 no.3
• /
• pp.75-80
• /
• 2001

The elastic plastic fracture toughness of this material is evaluated by the an unloading compliance method according to the ASTM E813-97 and E1152-97 method on the smooth and side groove 1CT specimens. The effect of smooth and side groove is studied on the material tearing modulus and characterizes the crack tip field under the plane stress and strain. SG-365 steel is observed that J-R curve and Tmat value decrease as 0%, 20%, 30%, and 40%. The 40% side grooved specimen is very useful in estimation of the $J_IC$. Because it is much easier than the smooth specimen to the onset of the ductile tearing by the R curve method. Besides. it improves the accuracy of toughness values, decreases the scattering the them and tunneling and shear lip by the side groove. Applicability of tearing modulus($T_J$ proposed by paris et al as instability panameter for this material is investigated.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.7 no.1
• /
• pp.1-10
• /
• 1983

The photoelastic determination of S.I.F. in Fracture mechanics has been regarded as one of the most effective and practical experimental methods in which stresses are read directly, except a few shortcomings involved in the process of experiment; the difficulties of making a sharp crack tip similar to the practical one and nearly impossibilities of carving an arbitrarily shaped crack on the test plate, etc. To eliminate flaws mentioned above, recently, Kitagawa and Watanabe of Tokyo Univ.developed a new method named"Teflon Insert Method". which has improved experimental accuracy to a considerable extent byt remaining still room for further improvement, that is, the elimination of bonding boundary scars which render photoelastic fringes obscure. In this paper, a newly exploited"Teflon Molding Method" was attempted for the completion of teflon-epoxy experimental method. The experimental results obtained by this method are compared with existent theoretical and experimental values to evaluate its accuracy. As result, 1-6% of margin of errors were appeared in a series of photoelastic experiments which defied any other conventional method in terms of experimental accuracy.perimental accuracy.

• Journal of the Korean Society for Precision Engineering
• /
• v.17 no.1
• /
• pp.153-163
• /
• 2000

The size effect on fracture toughness was investigated by introducing $J-A_2$ theory. For this application,small size specimens were chosen to establish $J-A_2$ assessment curve with FEM analysis. Two-dimensional FEM analysis was conducted with plane strain model using ABAQUS by domain integral method to calculate both crack tip stress and fracture toughness which were used to establish $J-A_2$ curve. The assessment curve predicted the fracture toughness of large specimens very well when compared to the test values. The results showed good prediction for deep crack specimen, though there were acceptable deviations in shallow cracked specimens, presumably caused by constraint effect. When the curve applied to reactor vessel in order to predict end of life fracture toughness with assumption of on-power pressure test condition, it provided the reasonable pressure compared to the existing design value. Better predictions would be possible if more test data were available.

• Proceedings of the Korean Society For Composite Materials Conference
• /
• 2001.05a
• /
• pp.106-109
• /
• 2001

In this work, we investigate the toughening mechanism of the rubber-modified epoxy resin. The fracture toughness($K_{IC}$) is measured using CT specimens for three kinds of rubber-modified epoxy resin with different rubber content. The damage zone and rubber particles around a crack tip of a damaged specimen just before fracture are observed by a polarization microscope and an atomic force microscope(AFM). Both the fracture energy($G_{IC}$) and the size of damage zone increase with the rubber content below l5wt%. The size of the rubber particles can be qualitatively correlated with the $G_{IC}$ and the size of damage zone. The cavitation of the rubber particles inside the damage zone is observed, which is expected to be main toughening mechanism by rubber particles. the stress which causes the cavitation of rubber particles is estimated by the Dugdale model.

• Journal of the Korean Society for Precision Engineering
• /
• v.4 no.4
• /
• pp.43-55
• /
• 1987

본 연구는 Orthotropic 재료에 존재하는 균열선단에서의 응력장의 양상에 대하여 지금까지의 연구들에서 상세하게 고려되지 않았던 여러 탄성계수들의 영향을 평면-변형조건 하에서 고찰하였다. 여러 하중조건과 경계조건 하에서의 일반적인 이방성 재료에 존재하는 균열선단에서의 응력자의 양상도 고찰하였으며, $a_{11}S^4-2a_{16}S^3+(2a_{12}+a_{66})S^2-2a_{26}S+a_{22}=0$인 특성방정식의 4근은 Cardano와 Ferrar의 해법을 쓴 새로운 단순기법에 의해 해석하였고 특성방정식의 4근에 미치는 일반적인 이방성 재료의 물성치 즉 Exx, Eyy, Ezz, Gxy, ${\nu}_{zx}$, ${\nu}_{xy}$,${\nu}_{yz}$, ${\eta}$ ${\eta}$xy,x, ${\eta}$xy, y, ${\eta}$xy,z , 들의 영향들에 대해 상세히 연구되었다. 또한 균열선단부근에서의 응력분포 양상을 재료의 물성치와 경계조건 및 하중조건들로써 체계적으로 연구하였다.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.2
• /
• pp.330-337
• /
• 1995

For a rotating body with cracks, the new energy release rate equation is presented. The derived equation is different from the other researcher's results. It is a path-independent integral which excluded the derivatives of displacements near the crack tip, thereby improving the numerical accuracy of the energy release rate computation. Moreover, as the equation was derived on basis of the energy principle and non-linear elasticity without assumptions, it can applied to the cracked body with arbitrary shape under elastic-plastic deformation. Several examples are treated to demonstrate the efficiency and accuracy of the proposed method compared to existing methods.