MPE Mathematical Problems in Engineering 1563-5147 1024-123X Hindawi Publishing Corporation 10.1155/2016/9794605 9794605 Research Article The Application of Nonordinary, State-Based Peridynamic Theory on the Damage Process of the Rock-Like Materials http://orcid.org/0000-0002-6080-9069 Gu X. B. 1 Wu Q. H. 2 Ma Ninshu 1 School of Civil Engineering Sichuan University of Science & Engineering Sichuan China scu.edu.cn 2 School of Architecture and Civil Engineering Chengdu University Chengdu China cdu.edu.cn 2016 10112016 2016 30 05 2016 03 08 2016 16 08 2016 2016 Copyright © 2016 X. B. Gu and Q. H. Wu. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Peridynamics has a great advantage over modeling the damage process of rock-like materials, which is assumed to be in a continuum interaction with each other across a finite distance. In the paper, an approach to incorporate classical elastic damage model in the nonordinary, state-based peridynamics is introduced. This method can model the dynamic damage process and stress change of rock-like materials. Then two instances about three-point bend experiment are simulated in the rock-like materials. Finally the conclusions are drawn that numerical results are close to the experimental results. So the method has a great predictable value in the geotechnical engineering.

Talent Introduction Projection in 2016 in Sichuan University of Science & Engineering 2016RCL19 Sichuan Provincial Department of Science 2016GFW0137 Sichuan Province Science and Technology Plan Project 2013JY0119
1. Introduction

The rock-like material is a quasibrittle material; it is widely applied to the geotechnical engineering, for example, tunnels and lots of other underground buildings. Especially for the mountain area, the research on the damage process of rock-like materials is essential, so the prediction for the damage process of rock-like materials becomes a focus problem gradually. But it is difficult to model because the mechanical character of rock-like materials is very complex and the precise damage prediction of rock-like materials is elusive. In the geotechnical engineering, the rock-like mass includes lots of tiny cracks, even before the load is applied. The damage process of the rock-like materials is often caused by tension, comparison, the change of the temperature, and so on, and the stress and displacement fields are influenced by the propagation and coalescence of the crack. Even through the damage process of rock-like materials has been investigated for many years, the damage mechanism and prediction of rock-like material are still not understood. In fact, the damage process of rock-like materials includes the initiation, propagation, and coalescence of the crack, which will lead to a sudden collapse due to brittle damage. How to model the damage character of the rock-like materials has brought great challenge.

Over the past decades, many methods are put forward to model the damage process of rock-like materials. In the finite element-based method, singular crack-tip elements are frequently encountered . Because of the crack-tip stress singularity, an external fracture criterion must be introduced to determine propagation and bifurcation of the cracks, and the nucleation question of the crack is still not solved . In order to overcome the above difficulties, the extended finite element theory  is proposed to simulate the propagation of cracks. Although many crack questions are solved by virtue of the extended finite element theory, external and bifurcation criterion must still be introduced when displacement is discontinuous and when interaction and bifurcation of multiple cracks are involved. Besides, a series of difficulties are encountered for the problem of the three-dimensional cracks by the XFEM. In order to solve the problems of the three-dimensional cracks, such as interactions among cracks and branching phenomenon of multiple cracks, meshless methods are developed . The propagation and coalescence process of cracks can be simulated by Smooth Particle Hydrodynamics (SPH); however, the tensile instability problems are still encountered in method of SPH . In order to avoid the aforementioned lacks, peridynamic theory, which is a numerical method based on the nonlocal thoughts, is introduced to model propagation and bifurcation process of cracks.

The peridynamic theory is a nonlocal meshless method; it is put forward by Silling , at Sandia National Laboratory. It is assumed that particles in a continuum interact with each other across a finite distance, and it formulates problems in terms of intergral equations rather than partial differential equations . Therefore, the peridynamic method can be applied to model the problems of continuous or discontinuous displacements .

After this theory is put forward, it has been widely applied to model the damage process of different materials. Firstly, the theory is applied to model the fracture process of composite material; for example, the fracture processes in laminated composites subjected to low-velocity impact and in woven composites subject to static indentation are predicted by Askari et al.  and Colavito et al. [10, 11] and the fatigue crack growth analysis is done in layered heterogeneous material system using peridynamic method by Jung and Soek . In addition, the notched laminated composite under biaxial loads is considered by Xu et al. . Then it is applied to model the fracture process of metal material; for example, Foster et al.  use peridynamic viscoplastic theory to model the collapse process of metal and Wu et al. [15, 16] analyze the ductile fracture of metal materials using nonordinary state-based peridynamics. Meanwhile, Sun and Sundararaghavan  model the crystal material using the peridynamic plasticity theory. Immediately, it is extended to model the damage of concrete; for example, Gerstle et al.  model the damage process of concrete by introducing “micropolar peridynamic model.” The damage results of these materials are rather good by using the peridynamics, but for the geotechnical engineering, rock-like materials are an important material, the peridynamic theory is seldom used to model the damage process of rock-like materials, and only Ha et al.  use peridynamics to model the fracturing patterns of rock-like materials in compression, but the stress field description is not considered. In the paper, the application of nonordinary, state-based peridynamics in the damage process of rock-like material will be investigated. Not only the damage process of rock-like materials is described by using the method, but also the change of stress field in the damage process is depicted, so it provides a new idea to model the damage process of rock-like materials by using the peridynamics.

The paper is organized as follows. In Section 2, the state-based peridynamic theory is introduced at first, and the implementation on this method for the specific damage model is discussed. In Section 3, the state-based peridynamic numerical discretization is described. In Section 4, we represent numerical results consisting of (1) the numerical simulation about the three-point bend test of the beam and (2) the numerical simulation about the three-point bend test of Brazilian disk. In Section 5, conclusions are drawn.

2. Model Description

For completeness, nonordinary state-based peridynamics is reviewed briefly; it includes a summary of basic peridynamic equation, the idea of constitutive relation, and its property.

2.1. Basic Theory

Peridynamic theory is put forward by Silling  at Sandia laboratory in the United States; its basic equation of motion is shown as follows:(1)ρxu¨x,t=Htkjuj-uk,xj-xk,t-tjkuk-uj,xk-xj,tdH+bx,t,where ρx is the density of material point x, H is a neighborhood (Figure 1), and bx,t is a prescribed body force density field of material point xk at the instance t, which represents the external force per unit reference volume square. tkjuj-uk,xj-xk,t and tjkuk-uj,xk-xj,t are the force density vector of material points x(k) and x(j), respectively.

The deformation of PD material point.

In the neighborhood of point x, the relative position of any point, x and x, is X_ξ, the deformation state of the bond is defined as Y_x-x, and they are expressed, respectively, as follows:(2)X_ξ=x-xY_x-x=u+x-u+x.

The nonlocal deformation gradient Fx of material point x is shown as the following expression:(3)Fx=HωξY_ξξdVξK-1x,where ωξ is the influence function of the bond. It is the function of the relative origin position ξ between points x and x, and Kx is a nonlocal shape tensor defined by(4)Kx=HωξξξdVξ-1.

Alternatively, this integral can be defined without recourse to peridynamic states as follows:(5)Fx,t=Hωx-xu-ux-xdVξK-1xKx=Hωx-xx-xx-xdVξ-1.

The Green-Lagrange strain sensor can be expressed as(6)C=12FTF-I.

The elastic strain energy W is written as(7)W=12C:ψ:C,where ψ is the fourth-order elastic tensor and is assumed as anisotropic with cubic symmetry.(8)ψijkn=ψijnk=ψknij=ψjikn.

From (7), the second Piola-Kirchhoff stress is written as(9)S=WE=ψ:C.

And Cauchy stress is found:(10)σ=FSdetFFT.

The first Piola-Kirchhoff stress can be written from the Cauchy stress as(11)P¯=detFσF-T.

Equation (9) is substituted into (10); the expression of the first Piola-Kirchhoff stress is rewritten as(12)P¯=FS.

Peridynamic force vector state t_x,tξ can be expressed by the traditional stress state as(13)t_x-x=ωx-xP¯Kxξ.

Substituting (12) into (13), (14) can be obtained as follows:(14)t_x-x=ωx-xFSKxξ.

From equations (1) and (14), the basic equation in the nonordinary, state-based peridynamics can be written as(15)ρxiu¨xi,t=Hωxi-xjFiSiKi-1xj-xi-FjSjKj-1xi-xjdVj+bxi,t.

2.2. Damage Correspondence

To describe the damage correspondence within the peridynamic framework, the influence function ωxi-xj is introduced to represent the damage; it can be expressed as follows:(16)ωxi-xj=χt,ξ1+δε,where δ is the horizontal radius of material point xi, ε=x-x is the relative displacement of any point x and point x, χt,ξ is a scalar function, and its expression is given:(17)χt,ξ=1st,ξ<s00st,ξ>s0,where s is the stretch of the bond; it is defined as(18)s=ε+η-εε.s0 is critical stretch; in two dimensions it is defined as (19)s0=Gc6μ/π+16κ-2μ/9π2δ,where κ is volume modulus, μ is the shear modulus, Gc is the critical energy dissipating ratio, and it is related to the fracture toughness KIC.

To model the problem of damage, the concept of local damage value is introduced; it is defined as(20)φx,t=1-Hχx,t,εdVHdVε,where φx,t is local damage value; its limitation is 0φx,t1, where 0 represents original materials and 1 represents complete disconnection of a point from all of the points with which it initially interact.

3. The Discretization

The region is discretized into nodes, each with a known volume in the reference configuration; taken together, the nodes form a grid. So the deformation gradient Fx is discretized into the following form:(21)Fj=i=1NjViωxi-xjyi-yjxi-xjKj-1,where a nonlocal shape sensor Kx is discretized as(22)Kj=i=1NjViωxi-xjxi-xjxi-xj,where Vi is the volume of nodes xi; Nj is the magnitude of node xj in the range of horizontal radius δ of the node xi.

Likewise, the discretization form of basic equation in the nonordinary, state-based peridynamics is shown as follows:(23)ρxiui¨n=j=1NiVjωxj-xiPi¯Ki-1xj-xi-Pj¯Kj-1xi-xj,where Pi¯ and Pj¯ are the first Piola-Kirchhoff stress of nodes i and j, respectively, Ki and Kj are the shape sensor of the nodes i and j.

The acceleration ui¨n can be shown as(24)ui¨n=uin+1-2uin+uin-1Δt2,where the superscript n represents the time step and the subscript represents the number of nodes, so ui¨n represents the acceleration of the node xi at time step n.

4. The Numerical Example

To evaluate the functionality of the proposed state-based peridynamic formulation, the discretized equations were implemented in a Fortran computer code for two-dimensional simulation. The following two examples are investigated, respectively.

4.1. The Three-Point Bend Test of the Beam

Three-point bend test  has been widely used for characterizing the dynamic response of materials as well as for validating numerical methods in the case of the state-based peridynamic formulation with damage. This test consists of a rectangular concrete beam with a preexisting vertical crack; its geometric configuration is shown in Figure 2. Its length is 0.32 m, the height is 0.07 m, the length of the crack in the beam is 0.0233 m, and the mechanic parameters of the concrete beam are listed as follows: the elastic modulus E=32.8 GPa, Poisson’s ratio ν=0.25, and center point of coordinate axis is positioned in the bottom center of concrete beam. x-axis is vertical to the pressure, tension is positive, and compression is negative. The concentrated load is positioned at the upper center of the beam. the magnitude of load L is 0.075 mm/min, the concrete beam model is discretized into 800×175=140000 particles, the distance between the adjacent two particles is Δx=4×10-4 m, time step dt=1.3367×10-8 s, the critical stretch s0=0.002, the density ρ=2650 kg/m3, and δ is adopted as 3Δx. When it arrives at 25000 time steps, the concrete beam is fractured fully. The computed damage contour of concrete beam is shown in Figure 3, the maximum principal stress σ contour in the damage process of concrete beam is shown in Figure 4, the damage result in specific test is given in Figure 5, and the comparison curve between test and numerical solution is shown in Figure 6.

The geometric configuration.

The damage contour of the beam.

1 . 3367 × 10 - 4  s

2.005 × 10 - 4  s

2.6734 × 10 - 4  s

3.34175 × 10 - 4  s

The maximum princess stress σ (Pa) contour of the concrete beam.

1.3367 × 10 - 4  s

2.005 × 10 - 4  s

2.6734 × 10 - 4  s

3.34175 × 10 - 4  s

The experimental result of three-point bend in the beam .

Load versus crack mouth displacement (CMD) curves obtained from the extended nonordinary, state-based peridynamic model and experiment.

It can be seen from Figure 3 that, when time arrives at 1.3367×10-4 s, the crack begins to initiate along the upper tip; when time arrives at 2.005×10-4 s, the crack begins to propagate along the position where the vertical load is applied; when time arrives at 2.6734×10-4 s, the crack continues to propagate and elongate along the position where the vertical load is applied, and two support positions in the bottom of concrete beam begin to be damaged; when time arrives at 3.34175×10-4 s, the crack propagates to the position where the concentrated load is applied, and the whole concrete beam broke down fully. It can be seen from Figure 4 that, when time arrives at 1.3367×10-4 s, the stress concentration begins to occur in the tip of crack, the stress concentration phenomenon is found at the position where the concentrated load is applied, and the magnitude of maximum principal stress in the tip of the crack reaches 1.5×106 Pa; when time arrives at 2.005×10-4 s, the stress concentration in the tip of concrete beam continues to propagate to the position where the concentrated load is applied, and the magnitude of maximum principal stress reaches 3.1×106 Pa; when time arrives at 3.34175×10-4 s, the magnitude of the stress in the tip of concrete beam reaches the maximum, whole concrete beam is destroyed fully, the magnitude of the maximum principal stress reaches 6.02×106 Pa, the numerical result in Figure 3 is close to the experimental observation in Figure 5, and from Figure 6, it can be found that the numerical process in the curve of load versus displacement is consistent with the trend of experimental results. These conclusions demonstrate that nonordinary, state-based peridynamics can provide good prediction for the damage process of rock-like materials. In the example, the difference of the damage result between the numerical prediction and experimental observation exists, and the constitutive relation of material has important influence on the difference, so the improvement for the constitutive relation of rock-like materials is necessary in my future research. It can not only improve the accuracy of numerical prediction, but also model the large deformation problem by altering the constitutive relation of materials, such as impact damage problem and blast. A critical extension is furnished to solve the problems involving severe deformation and damage by using the method.

This test consists of a circular rock sample with a preexisting vertical crack, its geometric configuration is shown in Figure 7, its diameter is 100 mm, there is a crack in the middle of sample, and its inclination angle is 0°. The length of crack is 2b = 30 mm, the mechanic parameters in this rock sample are listed as follows: the elastic modulus E=21 GPa, Poisson’s ratio ν=0.22, and center point of coordinate axis is positioned in the center of circular disk. x-axis is vertical to the pressure, tension is positive, and compression is negative. The compressive load is located at the upper and lower tips of the circular disk and their magnitude is 0.05 m/s; rock sample is discretized into 200×200=40000 particles, the distance between the adjacent two particles is Δx=5×10-4 m, time step dt=1.3367×10-8 s, the critical stretch s0=0.002, the density ρ=2300 kg/m3, and δ is adopted as 3Δx. When it arrives at 700 time steps, the concrete beam is fractured fully. The numerical damage contour of rock sample is shown in Figure 8, the maximum horizontal principal stress σx contour is shown in Figure 9, and the damage result in specific test is given in Figure 10.

The geometric configuration of rock sample.

The numerical contour about the failure of rock sample.

1.3367 × 10 - 4  s

4.01 × 10 - 2  s

6.68 × 10 - 2  s

9.36 × 10 - 2  s

The maximum horizontal princess stress σx (Pa) contour of the rock sample.

3.34 × 10 - 2  s

4.01 × 10 - 2  s

6.68 × 10 - 2  s

9.36 × 10 - 2  s

The experimental result of the rock sample .

It can be seen from Figure 8 that, when time arrives at 4.01×10-2 s, the upper and lower tips of the crack begin to initiate; when time arrives at 6.68×10-2 s, the crack along its tips continues to propagate; when time arrives at 9.36×10-2 s, the crack begins to penetrate along the direction where the concentrated load is applied. The whole rock sample fractures fully. Likewise, it can be seen from Figure 9 that, when time arrives at 3.34×10-2 s, the stress concentration in the tips of crack has taken place before the tips of crack begin to initiate, and the magnitude of the maximum principal stress in the tips of crack reaches 4×106 Pa; when time arrives at 4.01×10-2 s, the crack begins to initiate, and the magnitude of the maximum principal stress reaches 4.5×106 Pa; when time arrives at 6.68×10-2 s, the horizontal stress in the tips of crack continues to increase, and its magnitude reaches 5.5×106 Pa; when time arrives at 9.36×10-2 s, the magnitude of the maximum horizontal stress reaches 7.5×106 Pa; rock sample is fractured fully. The conclusions can be drawn from the comparison between Figures 8 and 10, the numerical results are similar to the experimental ones, and good effect can be obtained for the numerical simulation of the damage process in the rock material by using nonordinary, state-based peridynamic method, so it provides great instructional significance for the prediction of damage process in the future.

5. Conclusions

An extension of the state-based peridynamic constitutive correspondence framework to incorporate into elastic damage has been proposed. Then three-point bending simulation is performed with damage process. In two examples, the numerical results are compared with the experimental observation. In conclusion, the suggested method provides a new thought for the damage process prediction of the rock-like materials.

Nomenclature ρ ( x ) :

The density of material point x

b ( x , t ) :

The applied body force density of material point x at the time t

u j :

The displacement of material point x(j)

ξ :

The relative origin position between material points x and x

t k j u j - u k , x j - x k , t :

The force density vector of material point x(k)

t j k u k - u j , x k - x j , t :

The force density vector of material point x(j)

X _ ξ :

The relative position of points xand x

Y _ x - x :

The deformation state of the bond between material points x and x

F x :

The deformation gradient of material point x

ω ξ :

The influence function of the bond

V :

The volume of material point

W :

The elastic strain energy

S :

The second Piola-Kirchhoff stress

δ :

The horizontal radius of material point xi

χ t , ξ :

A scalar function of the bond between material points x and x at time t

s :

The stretch of the bond

κ :

The volume modulus

G c :

The critical energy dissipating ratio

φ x , t :

The local damage value

H :

The neighborhood of material points

ν :

Poisson’s ratio

C :

The Green-Lagrange strain sensor

ψ :

The elastic tensor

P ¯ :

The first Piola-Kirchhoff stress

s 0 :

The critical stretch of the bond

μ :

The shear modulus

K IC :

The fracture toughness

u i ¨ n :

The acceleration of the node xi at time step n.

Competing Interests

The authors declare that they have no competing interests.

Acknowledgments

This work is supported by the first batch of Talent Introduction Projection in 2016 in Sichuan University of Science & Engineering (no. 2016RCL19), the projection in the Sichuan Provincial Department of Science (no. 2016GFW0137), and financial contributions from the Sichuan Province Science and Technology Plan Project (2013JY0119).

Moës N. Dolbow J. Belytschko T. A finite element method for crack growth without remeshing International Journal for Numerical Methods in Engineering 1999 46 1 131 150 10.1002/(sici)1097-0207(19990910)46:1<131::aid-nme726>3.0.co;2-j 2-s2.0-0033349534 Belytschko T. Black T. Elastic crack growth in finite elements with minimal remeshing International Journal for Numerical Methods in Engineering 1999 45 5 601 620 10.1002/(sici)1097-0207(19990620)45:5<601::aid-nme598>3.0.co;2-s 2-s2.0-0032649861 Wagner G. J. Moës N. Liu W. K. Belytschko T. The extended finite element method for rigid particles in Stokes flow International Journal for Numerical Methods in Engineering 2001 51 3 293 313 10.1002/nme.169 MR1829318 Belytschko T. Krongauz Y. Organ D. Meshless methods: an overview and recent developments Computer Methods in Applied Mechanics and Engineering 1996 139 1–4 3 47 10.1016/s0045-7825(96)01078-x Libersky L. D. Petsehek A. G. Smooth particle hydrodynamics with strength of materials Advances in the Free-Lagrange Method Including Contributions on Adaptive Gridding and the Smooth Particle Hydrodynamics Method: Proceedings of the Next Free-Lagrange Conference Held at Jackson Lake Lodge, Moran, WY, USA 3–7 June 1990 1990 395 Berlin, Germany Springer 26 248 Lecture Notes in Physics 10.1007/3-540-54960-9_58 Silling S. A. Reformulation of elasticity theory for discontinuities and long-range forces Journal of the Mechanics and Physics of Solids 2000 48 1 175 209 10.1016/s0022-5096(99)00029-0 MR1727557 2-s2.0-0346055020 Silling S. A. Bobaru F. Peridynamic modeling of membranes and fibers International Journal of Non-Linear Mechanics 2005 40 2-3 395 409 10.1016/j.ijnonlinmec.2004.08.004 2-s2.0-6344294007 Silling S. A. Bathe K. J. Dynamic fracture modeling with a meshfree peridynamic code Computational Fluid and Solid Mechanics 2003 Amsterdam, The Netherlands Elsevier 641 644 10.1016/B978-008044046-0.50157-3 Askari E. Xu J. Silling S. Peridynamic analysis of damage and failure in composites Proceedings of the 44th AIAA Aerospace Sciences Meeting and Exhibit (AIAA '06) January 2006 Reno, Nev, USA 1123 1134 2-s2.0-34250771050 Colavito K. W. Kilic B. Celik E. Madenci E. Askari E. Silling S. A. Effect of void content on stiffness and strength of composites by a peridynamic analysis and static indentation test Proceedings of the 48th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics and Materials Conference April 2007 Waikiki, Hawaii, USA 23 26 Paper 2007–2251 Colavito K. W. Kilic B. Celik E. Madenci E. Askari E. Silling S. Effect of nano particles on stiffness and impact strength of composites Proceedings of the 48th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics, and Materials Conference April 2007 Honolulu, Hawaii, USA Paper 2007-2001 2-s2.0-34547514536 Jung J. Soek J. Fatigue crack growth analysis in layered heterogeneous material systems using peridynamic approach Composite Structures 2016 152 403 407 10.1016/j.compstruct.2016.05.077 Xu J. Askari A. Weckner O. Razi H. Silling S. Damage and failure analysis of composite laminates under biaxial loads Proceedings of the 48th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics, and Materials Conference April 2007 Honolulu, Hawaii, USA AIAA 7450 7458 10.2514/6.2007-2315 2-s2.0-34547524848 Foster J. T. Silling S. A. Chen W. W. Viscoplasticity using peridynamics International Journal for Numerical Methods in Engineering 2010 81 10 1242 1258 10.1002/nme.2725 2-s2.0-74949139199 Wu C. T. Ren B. A stabilized non-ordinary state-based peridynamics for the nonlocal ductile material failure analysis in metal machining process Computer Methods in Applied Mechanics and Engineering 2015 291 197 215 10.1016/j.cma.2015.03.003 MR3345247 2-s2.0-84928333104 Wu C. T. Ma N. Takada K. Okada H. A meshfree continuous-discontinuous approach for the ductile fracture modeling in explicit dynamics analysis Computational Mechanics 2016 58 3 391 409 10.1007/s00466-016-1299-6 Sun S. Sundararaghavan V. A peridynamic implementation of crystal plasticity International Journal of Solids and Structures 2014 51 19-20 3350 3360 10.1016/j.ijsolstr.2014.05.027 2-s2.0-84908099034 Gerstle W. Sau N. Silling S. Peridynamic modeling of concrete structures Nuclear Engineering and Design 2007 237 12-13 1250 1258 10.1016/j.nucengdes.2006.10.002 2-s2.0-34249999485 Ha Y. D. Lee J. Hong J.-W. Fracturing patterns of rock-like materials in compression captured with peridynamics Engineering Fracture Mechanics 2015 144 176 193 10.1016/j.engfracmech.2015.06.064 2-s2.0-84936871810 Weckner O. Abeyaratne R. The effect of long-range forces on the dynamics of a bar Journal of the Mechanics and Physics of Solids 2005 53 3 705 728 10.1016/j.jmps.2004.08.006 MR2116266 2-s2.0-12344302299 Zhou X.-P. Gu X.-B. Wang Y.-T. Numerical simulations of propagation, bifurcation and coalescence of cracks in rocks International Journal of Rock Mechanics & Mining Sciences 2015 80 241 254 10.1016/j.ijrmms.2015.09.006 2-s2.0-84943550301 Jenq Y. S. Shah S. P. Mixed-mode fracture of concrete International Journal of Fracture 1988 38 2 123 142 10.1007/BF00033002 2-s2.0-0024091584 Haeri H. Shahriar K. Marji M. F. Moarefvand P. Experimental and numerical study of crack propagation and coalescence in pre-cracked rock-like disks International Journal of Rock Mechanics and Mining Sciences 2014 67 20 28 10.1016/j.ijrmms.2014.01.008 2-s2.0-84894071009