• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.26 no.5A
• /
• pp.861-872
• /
• 2006

The element free Galerkin method is extended by the local partition of unity method to model the cohesive cracks in two dimensional continuum. The shape function of a particle whose domain of influence is completely cut by a crack is enriched by the step enrichment function. If the domain of influence contains a crack tip inside, it is enriched by a branch enrichment function which does not have the LEFM stress singularity. The discrete equations are obtained directly from the standard Galerkin method since the enrichment is only for the displacement field, which satisfies the local partition of unity. Because only particles whose domains of influence are influenced by a crack are enriched, the system matrix is still sparse so that the increase of the computational cost is minimized. The condition for crack growth in dynamic problems is obtained from the material instability; when the acoustic tensor loses the positive definiteness, a cohesive crack is inserted to the point so as to change the continuum to a discontiuum. The crack speed is naturally obtained from the criterion. It is found that this method is more accurate and converges faster than the classical meshless methods which are based on the visibility concept. In this paper, several well-known static and dynamic problems were solved to verify the method.

• Advanced Composite Materials
• /
• v.16 no.3
• /
• pp.223-244
• /
• 2007

This study experimentally and numerically investigated the tensile damage progress in stitched laminates. In particular, it focused on the effects of stitching on the damage progress. First, we experimentally confirmed that ply cracks and delamination appeared under load regardless of stitching. We then performed damage-extension simulation for stitched laminates using a layer-wise finite element model with stitch threads as beam elements, in which the damage (ply cracks and delamination) was represented by cohesive elements. A detailed comparison between observation and the simulated results confirmed that stitching had little effect on the onset and accumulation of ply cracks. Furthermore, we demonstrated that the stitch threads significantly suppressed the extension of the delamination.

• Geomechanics and Engineering
• /
• v.16 no.1
• /
• pp.49-58
• /
• 2018

Open-mode cracks could be commonly observed in layered rocks. A concept model is firstly used to explore the mechanism of the vertical cracks (VCs) in the top layer. Then the crack behaviour of the two-layer model is simulated based on a cohesive zone model (CZM) for layer interfaces and a plastic-damage model for rocks. The model indicates that the tensile stress normal to the VCs changes to compression if the crack spacing to layer thickness ratio is lower than a threshold. The results indicate that there is a threshold for interfacial shear strength that controls the crack patterns of the layered system. If the shear strength is lower than the threshold, the top layer is meshed by the VCs and interfacial cracks (ICs). When the shear strength is higher than the threshold, the top layer is meshed by the VCs and parallel cracks (PCs). If the shear strength is comparative to the threshold, a combining pattern of VCs, PCs and ICs for the top layer can be formed. The evolutions of stress distribution in the crack-bound block indicate that the ICs and PCs can reduce the load transferred for the substrate layer, and thus leads to a crack saturation state.

• 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 Korea Academia-Industrial cooperation Society
• /
• v.9 no.2
• /
• pp.286-290
• /
• 2008

This paper summarizes the components of an explicit aeroelastic solver developed especially for the simulation of dynamic fracture events occurring during the flight of solid propellant rockets. The numerical method combines an explicit Arbitrary Lagrangian Eulerian (ALE) version of the Cohesive Volumetric Finite Element (CVFE) scheme, used to simulate the spontaneous motion of one or more cracks propagating dynamically through a domain with regressing boundaries, and an explicit unstructured finite volume Euler code to follow the flow field during the failure event. A key feature of the algorithm is the ability to adaptively repair and expand the fluid mesh to handle the large geometrical changes associated with grain deformation and crack motion.

• Interaction and multiscale mechanics
• /
• v.1 no.1
• /
• pp.83-101
• /
• 2008

A derivation is given of two-scale models that are able to describe deformation and flow in a fluid-saturated and progressively fracturing porous medium. From the micromechanics of the flow in the cavity, identities are derived that couple the local momentum and the mass balances to the governing equations for a fluid-saturated porous medium, which are assumed to hold on the macroscopic scale. By exploiting the partition-of-unity property of the finite element shape functions, the position and direction of the fractures are independent from the underlying discretization. The finite element equations are derived for this two-scale approach and integrated over time. The resulting discrete equations are nonlinear due to the cohesive crack model and the nonlinearity of the coupling terms. A consistent linearization is given for use within a Newton-Raphson iterative procedure. Finally, examples are given to show the versatility and the efficiency of the approach.

• International Journal of Concrete Structures and Materials
• /
• v.2 no.2
• /
• pp.153-160
• /
• 2008

Fracture of prestressed concrete beams is studied with a novel and robust three-dimensional meshfree method. The meshfree method describes the crack as a set of cohesive crack segments and avoids the representation of the crack surface. It is ideally suited for a large number of cracks. The crack is modeled by splitting particles into two particles on opposite sides of the crack segment and the shape functions of neighboring particles are modified in a way the discontinuous displacement field is captured appropriately. A simple, robust and efficient way to determine, on which side adjacent particles of the corresponding crack segment lies, is proposed. We will show that the method does not show any "mesh" orientation bias and captures complicated failure patterns of experimental data well.

• Journal of the Korea Concrete Institute
• /
• v.17 no.6 s.90
• /
• pp.1053-1064
• /
• 2005

The mechanical damage of concrete is normally attributed to the formation of microcracks and their propagation and coalescence into macroscopic cracks. This physical degradation is caused from progressive and hierarchical damage of the microstructure due to debonding and slip along bimaterial interfaces at the mesoscale. Their growth and coalescence leads to initiation of hairline discrete cracks at the mesoscale. Eventually, single or multiple major discrete cracks develop at the macroscale. In this paper, from this conceptual model of mechanical damage in concrete, the computational efforts were made in order to characterize physical cracks and how to quantify the damage of concrete materials within the laws of thermodynamics with the aid of interface element in traditional finite element methodology. One dimensional effective traction/jump constitutive interface law is introduced in order to accommodate the normal opening and tangential slips on the interfaces between different materials(adhesion) or similar materials(cohesion) in two and three dimensional problems. Mode I failure and mixed mode failure of various geometries and boundary conditions are discussed in the sense of crack propagation and their spent of fracture energy under monotonic displacement control.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2004.10a
• /
• pp.395-402
• /
• 2004

The extended finite element scheme applied to crack problems is reviewed in this paper. As the enrichments of the solution space and the basic formulation are discussed, several examples of the application of the method are given. The examples include a LEFM crack, a cohesive crack, multiple LEFH cracks and dynamic crack propagation problems. It is shown that the extended finite element method is one of the powerful tools to study crack problems.

• Computers and Concrete
• /
• v.27 no.2
• /
• pp.131-141
• /
• 2021

This paper aims at investigating the capability of different FRP/concrete interface models to predict the effect of carbon nanotubes on the flexural behavior of RC beams strengthened with CFRP. Three different interfacial bond models are proposed to simulate the adhesion between CFRP composites and concrete, namely: full bond, nonlinear spring element, and cohesive zone model. 3D Nonlinear finite element model is developed then validated using experimental work conducted by the authors in a previous investigation. Cohesive zone model (CZM) has the best agreement with the experimental results in terms of load-deflection response. CZM is the only bond model that accurately predicted the cracks patterns and failure mode of the strengthened RC beams. The FE model is then expanded to predict the effect of bond strength on the flexural capacity of RC beams strengthened with externally bonded CNTs modified CFRP composites using CZM bond model. The results reveal that the flexural capacity of the strengthened beams increases with increasing the bond strength value. However, only 23% and 22% of the CFRP stress and strain capacity; in the case of full bond; can be utilized before failure.