• Proceedings of the KSR Conference
• /
• 2005.05a
• /
• pp.385-390
• /
• 2005

A specimen with residual stress due to welding was analyzed by three-dimensional cohesive zone model. The residual stress distribution was calculated by simulating welding process, and cohesive elements were located along crack propagation planes. Crack growth is possible since two planes of the cohesive element are separated beyond a maximum load carrying capacity. Stress fields around a crack tip are compared for specimens with and without residual stresses. Load-displacement curves and crack growth behaviors are also examined.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.33 no.6
• /
• pp.552-559
• /
• 2009

The effect of porosity on the crack propagation is studied by using the cohesive zone model. Standard mode I fracture test were done by using compact tension specimens with various porosities. Load-load line displacement curves and ${\delta}_5$-crack resistance curves for various porosities were obtained from experiments. The cohesive zone model proposed by Xu and Needleman was employed to describe the crack propagation in porous media, and the Gurson model is used for constitutive relation of porous materials. These models were implemented into user subroutines of a finite element program ABAQUS. The fracture mode changes from ductile fracture to brittle fracture as the porosity increases. Numerical calculations agree well with experimental results.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.37 no.12
• /
• pp.1453-1463
• /
• 2013

In this study, we investigated the effect of the indenter geometry on the crack characteristics by indentation cracking test and FEA. We conducted various cohesive finite element simulations based on the findings of Lee et al. (2012), who examined the effect of cohesive model parameters on crack size and formulated conditions for crack initiation and propagation. First, we verified the FE model through comparisons with experimental results that were obtained from Berkovich and Vickers indentations. We observed whether nonsymmetrical cracks formed beneath the surface during Berkovich indentation via FEA. Finally, we examined the relation between the crack size and the number of cracks. Based on this relation and the effect of the indenter angle on the crack size, we can predict from the crack size obtained with an indenter of one shape (such as Berkovich or Vickers) the crack size for an indenter of different shape.

• Advances in concrete construction
• /
• v.2 no.3
• /
• pp.229-247
• /
• 2014

A numerical study of the influence of shear-span/depth ratio on the cohesive crack fracture parameters and double - K fracture parameters of concrete is carried out in this paper. For the study the standard bending specimen geometry loaded with four point bending test is used. For four point loading, the shear - span/depth ratio is varied as 0.4, 1 and 1.75 and the ao/D ratio is varied from 0.2, 0.3 and 0.4 for laboratory specimens having size range from 100 - 500 mm. The input parameters for determining the double - K fracture parameters are taken from the developed fictitious crack model. It is found that the cohesive crack fracture parameters are independent of shear-span/depth ratio. Further, the unstable fracture toughness of double-K fracture model is independent of shear-span/depth ratio whereas, the initial cracking toughness of the material is dependent on the shear-span/depth ratio.

• Interaction and multiscale mechanics
• /
• v.2 no.1
• /
• pp.31-44
• /
• 2009

While the classical theory of Griffith is the foundation of modern understanding of brittle fracture, it has a number of significant shortcomings: Griffith theory does not predict crack initiation and path and it suffers from the presence of unphysical stress singularities. In 1998, Francfort and Marigo presented an energy functional minimization method, where the crack (or its absence) as well as its path are part of the problem's solution. The energy functionals act on spaces of functions of bounded variations, where the cracks are related to the discontinuity sets of such functions. The new model presented here uses modified energy functionals to account for molecular interactions in the vicinity of crack tips, resulting in Barenblatt cohesive forces, such that the model becomes free of stress singularities. This is done in a physically consistent way using recently published concepts of Sinclair. Here, for the consistency of the model, it becomes necessary to allow for crack reversibility and to consider local minimizers of the energy functionals. The latter is achieved by introducing different time scales. The model is solved in its global as well as in its local version for a simple one-dimensional example, showing that local minimization is necessary to yield a physically reasonable result.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.22 no.6
• /
• pp.519-525
• /
• 2009

In this study a cohesive zone model was used to simulate the delamination phenomena which occurs by a successive crack initiation and propagation in composite laminates. The cohesive zone model was incorporated to the classical finite element method via cohesive element formulation and then implemented into the user-subroutine UEL of a commercial finite element program Abaqus. To validate the formulation and implementation of the cohesive element the finite element results were compared with the experimental data of double cantilever beam and end notched flexure tests. The numerical results well agree with the experimental load-displacement curves. Also the effect of the elastic stiffness and the size of the cohesive element on the global load-displacement curves were studied numerically. To minimize the mesh-dependency of the crack propagation path and eliminate the zig-zag patterns in the load-displacement curve, cohesive elements should be refined at the crack-tip.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.31 no.4
• /
• pp.199-206
• /
• 2018

In this study, numerical analysis is applied to a two - dimensional model for verifying the general finite element program, Abaqus' s extended finite element method(XFEM). The cohesive element model used in the existing research has a limitation in simulating the actual crack because of the disadvantage that the crack path should be predicted and the element should be inserted. For this reason, the extended finite element method(XFEM), which predicts the path of cracks based on the directionality and specificity of stress, is emerging as a new solution in crack analysis. The validity of the XFEM application was confirmed by comparing the cohesive element analysis with the XFEM analysis by applying the crack path to the self - evident two - dimensional model. Numerical analysis confirms stress distribution and stress specificity immediately before crack initiation and compares it with actual crack initiation path. Based on this study, it is expected that cracks can be simulated by performing actual crack propagation analysis of complex models.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.32 no.5
• /
• pp.297-304
• /
• 2019

In this research, crack propagation in a silicon anode during two-phase lithiation was evaluated using a cohesive zone model. The phase transition from crystalline silicon to lithiated silicon causes compressive yielding due to the high volume expansion rate. Li-ion diffuses from the surface of the silicon to its core, and the complex deformation mechanisms during lithiation cause tensile hoop stress along the surface. The Park-Paulino-Roesler (PPR) potential-based cohesive zone model that guarantees consistent energy dissipation in mixed-mode fracture was adopted to simulate edge crack propagation. It was confirmed that the edge crack propagation characteristics during lithiation from the FEM simulation results coincided with the real experimental results. Crack turning observed from real experiments could also be predicted by evaluating the angles of maximum tensile stress directions.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.42 no.11
• /
• pp.911-918
• /
• 2014

In this paper, a selective activation strategy of cohesive elements using user subroutine UMPC was studied as an efficient solution for the added compliance problem in cohesive zone model crack propagation analyses. The cohesive elements were inserted between every bulk elements in region where cracks were expected to initiate and propagate, but initially not activated by tying the cohesive nodes using multi-point constraints. During analyses, the cohesive elements for which specified criterion was met were selectively activated by releasing the constraints. The effect of initial cohesive stiffness and the release criterion on the crack propagation behavior was carefully investigated.

• Structural Engineering and Mechanics
• /
• v.8 no.5
• /
• pp.439-452
• /
• 1999

A simple numerical scheme suitable for tracing the fracture propagation path for structures idealized by means of Hillerborg's classical cohesive crack model is presented. A direct collocation, multidomain boundary element method is adopted for the required space discretization. The algorithm proposed is necessarily iterative in nature since the crack itinerary is a priori unknown. The fracture process is assumed to be governed by a path-dependent generally nonlinear softening law. The potentialities of the method are illustrated through two examples.