• Proceedings of the KSME Conference
• /
• 2001.11a
• /
• pp.384-388
• /
• 2001

The higher order singularities[1] are systematically examined, and discussed are their complementarity relation with the nonsingular eigenfunctions and their relations to the configurational forces like J-integral and M-integral. By use of the so-called two state conservation laws(Im and Kim[2]) or interaction energy, originally proposed by Eshelby[3] and later treated by Chen and Shield[4], the intensities of the higher order singularities are calculated, and their roles in elasticplastic fracture are investigated. Numerical examples are presented for illustration.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.7 s.94
• /
• pp.1710-1718
• /
• 1993

This paper investigates the problem of a crack approaching two circular holes in an orthotropic infinite plate. The stress intensity factors were obtained by using the modified mapping-collocation method. The present results show excellent agreement with existing solutions for a crack approaching two circular holes in an isotropic infinite plate. In the numerical examples, various types of cross-ply laminated composites were considered. To investigate the effect of orthotropy and geometry(d/R and a/(d-R)) on crack tip singularity, stress intensity factors were considered as functions of the normalized crack length. It is expected that the modified mapping-collocation method can be applied to the analysis of various kinds of cracks existing around the stress-concentration region of composite laminate.

• Journal of the Korean Society of Safety
• /
• v.21 no.1 s.73
• /
• pp.15-20
• /
• 2006

The finite element alternating method is a convenient and efficient method to analyze three-dimensional cracks embedded in an infinite or a finite body because the method has the property that the uncracked body and cracks can be modeled independently. In this paper the method was applied for fatigue crack growth simulation. A surface crack in a cylinder was considered as an initial crack and the crack configurations and stress intensity factors during the crack growth were obtained. In this paper the finite element alternating method proposed by Nikishkov, Park and Atluri was used after modification. In the method, as the required solution for a crack in an infinite body, the symmetric Galerkin boundary element method formulated by Li and Mear was used. And a crack was modeled as distribution of displacement discontinuities, and the governing equation was formulated as singularity-reduced integral equations.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2002.10a
• /
• pp.607-614
• /
• 2002

For the accurate analysis of crack problems, considerable nodal refinement near the crack tip to capture singular stress field with sufficient accuracy to provide a useful computation of stress intensity factor is required. So, in this paper, adaptive nodal refinement scheme is proposed where nodes in restricted cell regions centered at crack tip are arranged in array for enhanced spatial resolution and adaptivity. With only cell-wise adaptive refinement scheme around crack tip fields, singularity of crack tip is sufficiently described to expect a successive crack propagate direction. Through numerical tests, accuracy of the proposed adaptive scheme is investigated and compared with the finite element and experimental results. By this implementation, it is shown that high accuracy is achieved by using iterative cell-wise solution method fur analyzing crack propagation problems.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2001.04a
• /
• pp.335-342
• /
• 2001

In this paper, an improved Element-Free Galerkin (EFG) method is proposed by adding enrichment function to the standard EFG approximation and a discontinuity function is implemented in constructing the shape function across the crack surface. In this method, the singularity and the discontinuity of the crack are efficiently modeled by using initial node distribution to evaluate reliable stress intensity factor, though the standard EFG method requires placing additional nodes near the crack tip. The proposed method enables the initial node distribution to be kept without any additional nodal d.o.f. and expresses the asymptotic stress field near the crack tip successfully. Numerical example verifies the improvement and the effectiveness of the method.

• Journal of the Microelectronics and Packaging Society
• /
• v.8 no.3
• /
• pp.25-30
• /
• 2001

The stress intensity factors for edge cracks located at the bonding interface between the semiconductor chip and the adhesive layer subjected to a uniform transverse tensile strain are investigated. Such cracks might be generated due to a stress singularity in the vicinity of the free surface. The boundary element method (BEM) is employed to investigate the behavior of interface stresses. The amplitude of complex stress intensity factor depends on the crack length, but it has a constant value at large crack lengths. The rapid propagation of interface crack is expected if the transverse tensile strain reaches a critical value.

• Proceedings of the KSME Conference
• /
• 2004.04a
• /
• pp.551-556
• /
• 2004

We propose a 2D 'crack' element for the simulation of propagating crack with minimal remeshing. A regular finite element containing the crack tip is replaced with this novel crack element, while the elements which the crack has passed are split into two transition elements. Singular elements can easily be implemented into this crack element to represent the crack-tip singularity without enrichment. Both crack element and transition element proposed in our formulation are mapped from corresponding master elements which are commonly built using the moving least-square (MLS) approximation only in the natural coordinate. In numerical examples, the accuracy of stress intensity factor $K_I$ is demonstrated and the crack propagation in a plate is simulated.

• Proceedings of the Earthquake Engineering Society of Korea Conference
• /
• 1999.04a
• /
• pp.256-266
• /
• 1999

The angled crack which is the simplest and representative case in the mixed mode crack analysis has stimulated the interests of many investigators during past 20 years. In this study the conventional quadratic isoparametric elements were used in all regions except the crack tip zone where triangular singularity elements with 6 nodes were used. The stress intensity factor of K1 and KII were determined respectively by the displacement correlation method. The finite element analysis program in this paper based on maximum energy release rate criteria and the results obtaiend by this program were compared with those calculated from the maximum circumferential tensile criteria and those by Jenq and Shah's experiments of the same geometry and material properties

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.14 no.3
• /
• pp.309-315
• /
• 2001

The Stress intensity factors for edge cracks located at the bonding interface between the elastic semiconductor chip and the viscoelastic adhesive layer have been investigated. Such cracks might be generated due to stress singularity in the vicinity of the free surface. The domain boundary element method(BEM) has been employed to investigate the behavior of interface stresses. The overall stress intensity factor for the case of a small interfacial edge crack has been computed. The magnitude of stress intensity factors decrease with time due to viscoelastic relaxation.

• Journal of the Korean Society of Safety
• /
• v.19 no.1
• /
• pp.38-43
• /
• 2004

Many analysis methods, including finite element method, have been suggested and used for assessing the integrity of cracked structures. In the paper, in order to analyze arbitrary three dimensional cracks, the finite element alternating method is extended. The crack is modeled by the symmetric Galerkin boundary element method as a distribution of displacement discontinuities, which is formulated as singularity-reduced integral equations. And the finite element method is used to calculate the stress values for the uncracked body only. Applied the proposed method to several example problems for planner cracks in finite bodies, the accuracy and efficiency of the method were demonstrated.