• Structural Engineering and Mechanics
• /
• v.11 no.6
• /
• pp.591-604
• /
• 2001

In this paper, a finite element with embedded displacement discontinuity which eliminates the need for remeshing of elements in the discrete crack approach is applied for the progressive fracture analysis of concrete structures. A finite element formulation is implemented with the extension of the principle of virtual work to a continuum which contains internal displacement discontinuity. By introducing a discontinuous displacement shape function into the finite element formulation, the displacement discontinuity is obtained within an element. By applying either a nonlinear or an idealized linear softening curve representing the fracture process zone (FPZ) of concrete as a constitutive equation to the displacement discontinuity, progressive fracture analysis of concrete structures is performed. In this analysis, localized progressive fracture simultaneous with crack closure in concrete structures under mixed mode loading is simulated by adopting the unloading path in the softening curve. Several examples demonstrate the capability of the analytical technique for the progressive fracture analysis of concrete structures.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1991.04a
• /
• pp.84-90
• /
• 1991

Underground structures usually consist of rock masses or concretes which can be cracked or have cracks. This study aims to develop an analysis program which can deal with the effect of discontinuous behavior due to those cracks using the block theory. It is assumed that rock masses form blocks along the discontinuity lines, and deformation within the block is relatively small. The behavior of discontinuity plane of the structures is divided into sliding along the discontinuity plane. separation of discontinuity by tensile force, and degradation of asperity angle of discontinuity plane by external force with sliding of rock Basses. These behaviors are implemented using constitutive relation and relevent load-displacement relation defined through normal and shear stiffnesses. Time varying displacements and block velocities are calculated by explicit time stepping algorithm. The effect of rock supports including rockbolts is also considered, and the tending effects which occurs in relatively thin lining is also considered.

• Coupled systems mechanics
• /
• v.7 no.6
• /
• pp.669-690
• /
• 2018

This paper deals with numerical modeling of dynamic failure phenomena in rate-sensitive brittle and/or ductile materials. To this end, a two-dimensional continuum viscodamage-embedded discontinuity model, which is based on our previous work (see Do et al. 2017), is developed. More specifically, the pre-peak nonlinear and rate-sensitive hardening response of the material behavior, representing the fracture-process zone creation, is described by a rate-dependent continuum damage model. Meanwhile, an embedded displacement discontinuity model is used to formulate the post-peak response, involving the macro-crack creation accompanied by exponential softening. The numerical implementation in the context of the finite element method exploiting the second-order mid-point scheme is discussed in detail. In order to show the performance of the model several numerical examples are included.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.9
• /
• pp.2138-2144
• /
• 1993

An integral equation representation of cracks was presented which differs from well-known "dislocation layer" representation. In this new representation, the equations are written in terms of the displacement discontinuity across the crack surfaces rather than derivatives of the displacement-discontinuity. It was shown in that the new technique is well-suited to the treatment of kinked cracks. In the present paper, this integral equation representation is coupled to the direct boundary-element method for the treatment of finite bodies containing kinked cracks. The method is demonstrated for two-dimensional finite domains but extension to three-dimensional problems would appear to be possible. The resulting approach is shown to be simple, yet very accurate. accurate.

• Journal of the Korean Society for Nondestructive Testing
• /
• v.27 no.6
• /
• pp.582-590
• /
• 2007

Imperfectly jointed interface serves as mechanical waveguide for elastic waves and gives rise to two distinct kinds of guided wave propagating along the interface. Contact acoustic nonlinearity (CAN) is known to plays major role in the generation of these interface waves called generalized Rayleigh waves in non-welded interface. Closed crack is modeled as non-welded interface that has nonlinear discontinuity condition in displacement across its boundary. Mathematical analysis of boundary conditions and wave equation is conducted to investigate the dispersive characteristics of the interface waves. Existence of the generalized Rayleigh wave(interface wave) in nonlinear contact interface is verified in theory where the dispersion equation for the interface wave is formulated and analyzed. It reveals that the interface waves have two distinct modes and that the phase velocity of anti-symmetric wave mode is highly dependent on contact conditions represented by linear and nonlinear dimensionless specific stiffness.

• Journal of the Korea Concrete Institute
• /
• v.11 no.3
• /
• pp.181-192
• /
• 1999

The progressive failure following strain localization in concrete can be analyzed effectively using finite element modeling of fracture process zone of concrete with a finite element embedded discontinuity. In this study, a finite element with embedded discontinuous line is utilized for the analysis of progressive failure in concrete. The finite element with embedded discontinuity is a kind of discrete crack element, but the difficulties in discrete crack approach such as remeshing or adding new nodes along with crack growth can be avoided. Using a discontinuous shape function for this element, the displacement discontinuity is embedded within an element and its constitutive equation is modeled from the modeling of fracture process zone. The element stiffness matrix is derived and its dual mapping technique for numerical integration is employed. Then, a finite element analysis program with employed algorithms is developed and failure analysis results using developed finite element program are verified through the comparison with experimental data and other analysis results.

• Journal of the Korean Geotechnical Society
• /
• v.15 no.4
• /
• pp.221-234
• /
• 1999

This paper deals with the application of joint element in the finite element modeling of discontinuities encountered during rock tunneling. A nodal displacement joint element was implemented in a two dimensional finite element program GEOFE2D. The applicability of the joint element for modeling of discontinuities and the numerical stability of the implemented algorithm were examined by comparing the results of reduced small scale model tests as well as commercially available FEM program. The GEOFE2D was then used to analyze a tunnel crossed by a major discontinuity for the purpose of understanding the effect of discontinuity on the tunnel behavior. In addition, a modeling technique for the junction of discontinuity and shotcrete lining was presented. The results of analysis indicated that the stress-strain field around the tunnel is significantly altered by the presence of discontinuity, and that the stresses in the shotcrete lining considerably increase at the junction of the shotcrete lining and the discontinuity. It is therefore concluded that the major discontinuities must be carefully modeled in the finite element analysis of a tunneling problem in order to obtain more reliable results close to actual tunnel behavior.

• Computers and Concrete
• /
• v.5 no.4
• /
• pp.375-388
• /
• 2008

Crack opening governs many transfer properties that play a pivotal role in durability analyses. Instead of trying to combine continuum and discrete models in computational analyses, it would be attractive to derive from the continuum approach an estimate of crack opening, without considering the explicit description of a discontinuous displacement field in the computational model. This is the prime objective of this contribution. The derivation is based on the comparison between two continuous variables: the distribution if the effective non local strain that controls damage and an analytical distribution of the effective non local variable that derives from a strong discontinuity analysis. Close to complete failure, these distributions should be very close to each other. Their comparison provides two quantities: the displacement jump across the crack [U] and the distance between the two profiles. This distance is an error indicator defining how close the damage distribution is from that corresponding to a crack surrounded by a fracture process zone. It may subsequently serve in continuous/discrete models in order to define the threshold below which the continuum approach is close enough to the discrete one in order to switch descriptions. The estimation of the crack opening is illustrated on a one-dimensional example and the error between the profiles issued from discontinuous and FE analyses is found to be of a few percents close to complete failure.

• Computers and Concrete
• /
• v.7 no.6
• /
• pp.483-496
• /
• 2010

This work deals with the determination of crack openings in 2D reinforced concrete structures using the Finite Element Method with a smeared rotating crack model or an embedded crack model. In the smeared crack model, the strong discontinuity associated with the crack is spread throughout the finite element. As is well known, the continuity of the displacement field assumed for these models is incompatible with the actual discontinuity. However, this type of model has been used extensively due to the relative computational simplicity it provides by treating cracks in a continuum framework, as well as the reportedly good predictions of reinforced concrete members' structural behavior. On the other hand, by enriching the displacement field within each finite element crossed by the crack path, the embedded crack model is able to describe the effects of actual discontinuities (cracks). This paper presents a comparative study of the abilities of these 2D models in predicting the mechanical behavior of reinforced concrete structures. Structural responses are compared with experimental results from the literature, including crack patterns, crack openings and rebar stresses predicted by both models.

• Computers and Concrete
• /
• v.1 no.3
• /
• pp.227-248
• /
• 2004

A combined experimental-computational study of a double edge-notched stone specimen subjected to tensile loading is presented. In the experimental part, the load-deformation response and the displacement field around the crack tip are recorded. An Electronic Speckle Pattern Interferometer (ESPI) is used to obtain the local displacement field. The experimental results are used to validate a numerical model for the description of fracture using finite elements. The numerical model uses displacement discontinuities to model cracks. At the discontinuity, a plasticity-based cohesive zone model is applied for monotonic loading and a combined damage-plasticity cohesive zone model is used for cyclic loading. Both local and global results from the numerical simulations are compared with experimental data. It is shown that local measurements add important information for the validation of the numerical model. Consequently, the numerical models are enhanced in order to correctly capture the experimentally observed behaviour.