Simulation of time-dependent crack growth in brittle rocks under constant loading conditions

https://doi.org/10.1016/j.engfracmech.2014.02.008Get rights and content

Highlights

  • A lifetime prediction scheme was developed for rock under constant loadings.

  • Numerical simulations were performed utilizing the developed modeling scheme.

  • The damage process and macroscopic fracture pattern of the models were studied.

  • Typical fracture patterns were observed in the numerical models.

  • Preliminary studies were performed through numerical modeling.

Abstract

Based on the theory of subcritical crack growth, linear elastic fracture mechanics (LEFM), and Charles equation, a lifetime prediction scheme has been developed for rock specimens containing initial microcracks under constant loadings. Numerical simulations were performed utilizing the developed modeling scheme. Lifetimes were obtained through numerical calculation; the damage process and macroscopic fracture pattern of the models were studied. Typical fracture pattern, like tensile cracks and shear bands, were observed. A few preliminary studies were also performed to compare the results with in situ observations. Conclusions were drawn and possible improvements to future research work are proposed.

Introduction

Experiments and theoretical research on the phenomenon of fracture and damage started centuries ago. But no quantitative results were obtained until the great work of Griffith [1], which had lead to the onset of modern fracture mechanics. Griffith analyzed the stresses of a cracked plate, which was studied before by Inglis [2], and develop a theory for crack initiation and propagation. Orowan [3] modified Griffith’s theory by including the influence of plasticity into the energy balance. Irwin [4] introduced the parameter later known as the stress intensity factor K. After the fundamentals of linear elastic fracture mechanics were well established around 1960 [5], researchers had focused more on the influence of materials’ plasticity on the fracture analysis. One trend in recent fracture and damage mechanical research is related to time-dependent studies (life time prediction or time to failure prediction) for rocks (e.g. Kemeny [6], [7], [8], Mishnaevsky [9], Shao et al. [10], Li et al. [11], Rinne [12], Konietzky et al. [13]).

Based on the linear elastic fracture mechanical theory (LEFM), numerical simulations of time-dependent fracture propagation on Westerly Granite have been performed by Konietzky et al. [13]. The innovations of their research work include the simulation of specific initial cracks with certain distributions at the micro-scale, and the simulation of sub-critical and critical fracture growth to describe the time-dependent damage process until failure. It is assumed, that the damage and finally the failure of a rock specimen is the result of the partially parallel growth and coalescence of many initially existing microcracks at the grain size level rather than the growth of one or only a few single cracks. Other studies, like the experimental investigations on dynamic fracture performed by Ravi-Chandar and Knauss [14], [15], also implied this concept. As a further improvement of the work of Konietzky et al. [13], this study includes the influence of orientation distribution of the initial microcracks on the crack growth, and an adopted wing crack propagation scheme to describe the growth of initial microcracks.

Section snippets

Theoretical basis

Irwin [4] introduced stress intensity factor to describe the stress distribution and displacements near the crack tip of brittle materials. A fracture criterion based on the stress intensity factor can be described as follows: failure of the material occurs when the stress intensity factor reaches the critical value, the so-called fracture toughness KC. It was assumed by classical LEFM that the crack will propagate ultrasonically when the stress intensity factor K reaches the fracture toughness

Numerical simulations

The numerical calculations were performed with FLAC in plain strain mode using the internal program language FISH (Itasca [32]). Data from Westerly Granite were used for the numerical model, as is shown in Table 1.

Applications

Pre-cracked specimens under uniaxial load were studied utilizing the proposed modeling scheme. Size and geometry of the numerical models are shown in Fig. 11(a) and (b). The simulation results are compared with those obtained by lab tests on Hwangdeung granite specimens (Lee and Jeon [33]), which are of smaller scale, but the ratio between crack size and specimen size is the same (Fig. 11(b)). The numerical model is divided into 20,000 zones with each zone of the size 0.04 × 0.04 m2. The initial

Conclusions

A time-dependent crack propagation scheme is developed to simulate subcritical and critical crack growth of microcracks in brittle rock. The development of macroscopic crack formed by the coalescence of microcracks is simulated utilizing the proposed modeling scheme. The time-related macroscopic failure of a rock sample is studied under different loading conditions. Some preliminary studies are also performed to compare modeling results with typical in situ observations. Some conclusions are

Acknowledgment

The authors thank the anonymous reviewers for their valuable hints and recommendations for the improvement of this paper.

References (34)

  • E. Orowan

    Fracture and strength of solids

    Rep Prog Phys

    (1949)
  • G.R. Irwin

    Analysis of stresses and strains near the end of a crack traversing a plate

    J Appl Mech

    (1957)
  • T.L. Anderson

    Fracture mechanics

    (1995)
  • J. Kemeny

    The time-dependent reduction of sliding cohesion due to rock bridges along discontinuities: a fracture mechanics approach

    Rock Mech Rock Engng

    (2003)
  • L.L. Mishnaevsky

    Determination for the time-to-fracture of solids

    Int J Fracture

    (1996)
  • J.T. Li et al.

    Crack growth time dependence analysis of granite under compressive-shear stresses state

    J Coal Sci Engng (China)

    (2008)
  • Rinne M. Fracture mechanics and subcritical crack growth approach to model time-dependent failure in brittle rock....
  • Cited by (20)

    • Application of equivalent circle and ellipse for pore shape modeling in crack growth problem: A numerical investigation in microscale

      2021, Engineering Fracture Mechanics
      Citation Excerpt :

      Among all parameters, pore shape (or pore geometry) is one of the important parameters controlling the crack propagation and the yield strength [7,8,12,14,15], but has not been fully investigated in microscale. Although, various methods have been proposed to study crack growth and failure in porous media, numerical simulation has attracted many attentions due to its flexibility, unique features and high ability to predict the crack growth path [13,14,16–21]. For example, Jiang et al. [16] used eXtended Finite Element Method (XFEM) to investigate the effects of pores, inclusions and minor cracks on the path of major crack propagation.

    • Weibull linear parallel bond model (WLPBM) for simulating micro-mechanical characteristics of heterogeneous rocks

      2019, Engineering Analysis with Boundary Elements
      Citation Excerpt :

      However, considerable limitations still permeate the quantitative characterization of the mechanical properties of heterogeneous rocks. Alternatively, numerical methods have enjoyed great success in recent years to investigate the effect of the heterogeneity on the microcracking processes of rocks [26–33]. In principle, the developed approaches available for simulating the heterogeneous characteristics of rocks can be classified into two main categories.

    • A unified analytical method calculating brittle rocks deformation induced by crack growth

      2019, International Journal of Rock Mechanics and Mining Sciences
      Citation Excerpt :

      Fracture geometries were kept track of independently of the mesh allowing geometric handling of fracture arrest, closure, and coalescence. Li and Konietzky27 presented a lifetime prediction method for rocks containing initial microcracks under constant stress. Lifetimes and deformations were obtained by numerical calculation.

    • A microcrack growth-based constitutive model for evaluating transient shear properties during brittle creep of rocks

      2018, Engineering Fracture Mechanics
      Citation Excerpt :

      Based on the subcritical crack law and the tensile microcrack mechanism, an analytical micromechanics-based method was proposed to study mechanical behaviors during brittle creep [30], and the relation between shear strength and microcrack size in rocks was analyzed by inserting the Mohr-Coulomb failure criterion [31]. Moreover, embedding the theory of linear elastic fracture mechanics, Chares equation, and subcritical crack law into the FLAC software, a numerical method was proposed to simulate the rock creep deformation caused by crack growth, the spatial effect of rock failure from strain localization was studied, and the lifetime of rock was discussed [27–29]. However, the transient shear properties caused by subcritical crack growth cannot be described during brittle creep in intact rocks.

    View all citing articles on Scopus
    View full text