• Journal of Welding and Joining
• /
• v.21 no.6
• /
• pp.71-76
• /
• 2003

For the purpose of clarifying the influence of welding residual stress to the fatigue crack propagations behaviour, an analytical investigation based on finite element method is performed to examine the opening behaviour of tip-closed crack in the compressive residual stress. A finite element model comprised of contact elements for the crack plane and plane stress elements for the base material is used to evaluate crack opening stress of the crack existing in the residual stress field. Also an analytical method based on the superposition principle to estimate the length of opened part of tip closed crack and the stress distribution adjacent to the crack during uploading is applied to the finite element model. The software for the analysis is ABAQUS, which is a general purpose finite element package. The results show that stresses distributed on the crack surfaces are reduced and approached to zero as the applied stresses are increased up to crack tip opening stress and no mechanical discontinuity is found at the boundary of contact elements and plane stress elements. It is verified that the opening behavior of the fatigue crack in the residual stress can be predicted by finite element method with the proposed analytical method.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.27 no.10
• /
• pp.1703-1711
• /
• 2003

An elastic-plastic finite element analysis is performed to investigate detailed closure behaviour of fatigue cracks in residual stress fields and the numerical results are compared with experimental results. The finite element analysis performed under plane stress using contact elements can predict fatigue crack closure behaviour. The mesh of constant element size along crack surface can not predict the opening level of fatigue crack. Specially, the mesh of element sizes depending upon the reversed plastic zone size included the effect of crack opening point can precisely predict the opening level. By using the concept of the mesh of element sizes depending upon the reversed plastic zone size included the effect of crack opening point, the opening level of fatigue crack can be determined very well.

• 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.

• Journal of Welding and Joining
• /
• v.21 no.6
• /
• pp.77-83
• /
• 2003

The finite element analysis was performed for the cracks existing in residual stress fields in order to investigate the effects of configuration of residual stress distribution to the fatigue crack opening behaviour. And the variation of stress distributions adjacent to the crack caused by uploading was examined. The finite element model with contact elements for the crack plane and plane stress elements for the base material and the analytical method based on the superposition principle to estimate crack opening behaviour and the stress distribution adjacent to the crack subjected to uploading were used. The results of the analysis showed that crack opening behaviors and variations of stress distribution caused by uploading were changed depending on the configuration of residual stress distribution. When the crack existed in the region of compressive residual stress and the configuration of compressive residual stress distribution were inclined, a partial crack opening just behind of a crack tip occurred during uploading. Based on the above results, it was clarified that the crack opening behaviour in the residual stress field could be predicted accurately by the finite element analysis using these analytical method and model.

• 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.

• Journal of Welding and Joining
• /
• v.21 no.7
• /
• pp.71-77
• /
• 2003

In this study, an investigation based on the superposition principle to predict residual stress redistribution caused by crack propagation itself initially through residual tensile stress field was performed by finite element method. The tendency in residual stress redistribution caused by crack propagation recognized both from the analytical results and experimental result was the residual stress concentration consecutively occurred in the vicinity of crack tip even the situation that the crack propagated to the region initially residual compressive stress existed. The software for the analysis is ABAQUS, which is a general purpose finite element package. The analytical method that attempt to take the plastic deformation at the crack tip due to tensile residual stress into the consideration of residual stress redistribution caused by crack propagation was proposed. The plastic zone size at the tip of fatigue crack and redistributed residual stresses were calculated by finite element method on the bases of the concept of Dugdale model. Comparing these analytical results with experimental results, it is verified that the residual stress redistribution caused by crack propagation can be predicted by finite element method with the proposed analytical method.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.31 no.10
• /
• pp.1009-1016
• /
• 2007

Finite element alternating method has been suggested and used effectively to obtain the fracture parameters in assessing the integrity of cracked structures. The method obtains the solution from alternating independently between the FEM solution for an uncracked body and the crack solution in an infinite body. In the paper, the finite element alternating method is extended in order to obtain the elastic-plastic stress fields of a three dimensional inner crack. The three dimensional crack solutions for an infinite body were obtained using symmetric Galerkin boundary element method. As an example of a three dimensional inner crack, a penny-shaped crack in a finite body was analyzed and the obtained elastc-plastic stress fields were compared with the solution obtained from the finite element analysis with fine mesh. It is noted that in the region ahead of the crack front the stress values from FEAM are close to the values from FEM. But large discrepancy between two values is observed near the crack surfaces.

• Journal of Mechanical Science and Technology
• /
• v.14 no.4
• /
• pp.401-407
• /
• 2000

An elastic-plastic finite element analysis is performed to investigate detailed closure behavior of fatigue cracks and the numerical results are compared with experimental results. The finite element analysis performed under plane stress using 4-node isoparametric elements can predict fatigue crack closure behavior. The mesh of constant element size along crack surface can not predict the opening level of fatigue crack. The crack opening level for the constant mesh size increases linearly from initial crack growth. The crack opening level for variable mesh size, is almost flat after crack tip has passed the monotonic plastic zone. The prediction of crack opening level using the variable mesh size proportioning the reversed plastic zone size with the opening stress intensity factors presents a good agreement with the experimental data regardless of stress ratios.

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 2000.11a
• /
• pp.517-522
• /
• 2000

When a body including a crack inside is subjected to the compressive forces, the crack is closed and sliding occurs on the crack surfaces. In this work, a subsurface crack subjected to a static or moving compressive load is analyzed with the finite element method considering friction on the crack surface. The friction on the crack surface is assumed to follow the Coulomb friction law. A numerical method based on the finite element method and iterative method is applied in this work. And the result is compared with the frictional contact of crack by ANSYS using contact 12 element. The numerical results of two methods are compared with the wellknown analytical solutions, and the accuracy of iterative method is checked..

• Journal of Mechanical Science and Technology
• /
• v.18 no.11
• /
• pp.1989-1995
• /
• 2004

Finite element analysis(FEA) is the most popular numerical method to simulate plasticity-induced fatigue crack closure and can predict fatigue crack closure behavior. Finite element analysis under plane stress state using 4-node isoparametric elements is performed to investigate the detailed closure behavior of fatigue cracks and the numerical results are compared with experimental results. The mesh of constant size elements on the crack surface can not correctly predict the opening level for fatigue crack as shown in the previous works. The crack opening behavior for the size mesh with a linear change shows almost flat stress level after a crack tip has passed by the monotonic plastic zone. The prediction of crack opening level presents a good agreement with published experimental data regardless of stress ratios, which are using the mesh of the elements that are in proportion to the reversed plastic zone size considering the opening stress intensity factors. Numerical interpolation results of finite element analysis can precisely predict the crack opening level. This method shows a good agreement with the experimental data regardless of the stress ratios and kinds of materials.