• Structural Engineering and Mechanics
• /
• 제11권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.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 1991년도 봄 학술발표회 논문집
• /
• 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
• /
• 제7권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.

• 대한기계학회논문집
• /
• 제17권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.

• 비파괴검사학회지
• /
• 제27권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.

• 콘크리트학회논문집
• /
• 제11권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.

• 한국지반공학회논문집
• /
• 제15권4호
• /
• pp.221-234
• /
• 1999

본 논문에서는 불연속암반내 시공되는 터널의 유한요소해석에 있어서 절리요소를 이용한 불연속면의 모델링에 관한 내용을 다루었다. 불연속 암반터널의 모델링이 가능한 유한요소해석 프로그램의 개발을 위해 기존의 유한요소해석 프로그램 GEOFE2D에 불연속면의 모델링이 가능한 절점변위 절리요소를 적용하고 모형실험 및 기존의 상용프로그램과의 비교를 통해 그 타당성을 검증하였다. 또한 검증된 GEOFE2D를 이용하여 불연속면이 터널의 거동에 미치는 영향을 고찰하기 위해 불연속면이 터널을 관통하는 경우에 대한 해석을 수행하고, 그 과정에서 불연속면과 숏크리트 라이닝 교차부에서의 변위 적합조건을 만족시킬 수 있는 불연속면 모델링 기법을 제시하였다. 한편, 해석결과를 분석한 결과 불연속면은 터널 주변의 응력-변형률 상태에 현저한 영향을 미치며, 특히 불연속면이 관통하는 부위에서의 숏크리트 라이닝 축력 및 휨 모멘트가 현저히 증가하는 것으로 나타났다. 따라서, 불연속면을 포함하는 터널의 거동해석시에는 불연속면에 대한 보다 상세한 모델링이 수반되어야 실제 거동에 보다 근접하는 해석결과를 도출시킬 수 있을 것으로 판단된다.

• Computers and Concrete
• /
• 제5권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
• /
• 제7권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
• /
• 제1권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.