Numerical simulation of fatigue crack propagation in heterogeneous geomaterials under varied loads using displacement discontinuity method

Heterogeneous brittle geomaterials are highly susceptible to cyclic loads.They contain inherent flaws and cracks that grow under fatigue loads and lead to failure.This study presents a numerical model for analyzing fatigue in these materials based on the two-dimensional(2D)boundary element method and linear elastic fracture mechanics.The process is formulated by coupling the displacement discontinuity method with the incorporation technique of dissimilar regions and the governing equations of fatigue.The heterogeneous media are assumed to consist of materials with different properties,and the interfaces are assumed to be completely bonded.In addition,the domains include multiple cracks exposed to constant and variable amplitude cyclic loads.The stress intensity factor is a crucial parameter in fatigue analysis,which is determined using the displacement field around crack tips.An incremental crack growth scheme is applied to calculating the fatigue life.The growth rate values are employed to estimate the length of crack extension when there are multiple cracks.The interaction between cracks is considered,which also includes the coalescence phenomenon.Finally,various structures under different cyclic loads are examined to evaluate the accuracy of this method.The results demonstrate the efficiency of the proposed approach in modeling fatigue crack growth and life estimation.The behavior of life curves for the heterogeneous domain was as expected.These curves illustrate the breakpoints caused by utilizing discrete incremental life equations.At these points,the trend of the curves changed with the material properties and fatigue characteristics of the new material around the crack tips.

Propagation and Coalescence of Blast-Induced Cracks in PMMA Material Containing an Empty Circular Hole Under Delayed Ignition Blasting Load by Using the Dynamic Caustic Method

In this paper,dynamic caustic method is applied to analyze the blast-induced crack propagation and distribution of the dynamic stress field around an empty circular hole in polymethyl methacrylate(PMMA)material under delayed ignition blasting loads.The following experimental results are obtained.(1)In directional-fracture-controlled blasting,the dynamic stress intensity factors(DSIFs)and the propagation paths of the blast-induced cracks are obviously influenced by the delayed ignition.(2) The circular hole situated between the two boreholes poses a strong guiding effect on the coelesence of the cracks,causing them to propagate towards each other when cracks are reaching the circular hole area.(3)Blast-induced cracks are not initiated preferentially because of the superimposed effect from the explosive stress waves on the cracking area.(4) By using the scanning electron microscopy(SEM)method,it is verified that the roughness of crack surfaces changes along the crack propagation paths.

Shear band evolution and acoustic emission characteristics of sandstone containing non-persistent flaws

Direct shear tests were conducted on sandstone specimens under different constant normal stresses to study the coalescence of cracks between non-persistent flaws and the shear sliding characteristics of the shear-formed fault.Digital image correlation and acoustic emission(AE)techniques were used to monitor the evolution of shear bands at the rock bridge area and microcracking behaviors.The experimental results revealed that the shear stresses corresponding to the peak and sub-peak in the stressdisplacement curve are significantly affected by the normal stress.Strain localization bands emerged at both the tip of joints and the rock bridge,and their extension and interaction near the peak stress caused a surge in the AE hit rate and a significant decrease in the AE b value.Short and curvilinear strain bands were detected at low normal stress,while high normal stress generally led to more microcracking events and longer coplanar cracks at the rock bridge area.Furthermore,an increase in normal stress resulted in a higher AE count rate and more energetic AE events during friction sliding along the shearformed fault.It was observed that the elastic energy released during the crack coalescence at the prepeak stage was much greater than that released during friction sliding at the post-peak stage.More than 75%of AE events were located in the low-frequency band(0e100 kHz),and this proportion continued to rise with increasing normal stress.Moreover,more AE events of low AF value and high RA value were observed in specimens subjected to high normal stress,indicating that greater normal stress led to more microcracks of shear nature.

Experimental study on cracking behaviour of moulded gypsum containing two non-parallel overlapping flaws under uniaxial compression

Failure of rock mass that is subjected to compressive loads occurs from initiation, propagation, and linkage of new cracks from preexisting fissures. Our research investigates the cracking behaviour and coalescence process in a brittle material with two non-parallel overlapping flaws using a high-speed camera. The coalescence tensile crack and tensile wing cracks were the first cracks to occur from the preexisting flaws. The initiation stresses of the primary cracks at the two tips of each flaw were simultaneous and decreased with reduced flaw inclination angle. The following types of coalescence cracks were identified between the flaws: primary tensile coalescence crack, tensile crack linkage, shear crack linkage, mixed tensile-shear crack, and indirect crack coalescence. Coalescence through tensile linkage occurred mostly at pre-peak stress. In contrast, coalescence through shear or mixed tensile-shear cracks occurred at higher stress. Overall, this study indicates that the geometry of preexisting flaws affect crack initiation and coalescence behaviour.

Discrete element modeling on the crack evolution behavior of brittle sandstone containing three fissures under uniaxial compression

Based on experimental restilts of brittle, intact sandstone under uniaxial compression, the micro-parameters were firstly confirmed by adopting particle flow code （PFC2D）. Then, the validation of the simulated models were cross checked with the experimental results of brittle sandstone containing three parallel fissures under uniaxial compression. The simulated results agreed very well with the experimental results, including the peak strength, peak axial strain, and ultimate failure mode. Using the same micro- parameters, the numerical models containing a new geometry of three fissures are constructed to investigate the fissure angle on the fracture mechanical behavior of brittle sandstone under uniaxial compression. The strength and deformation parameters of brittle sandstone containing new three fissures are dependent to the fissure angle. With the increase of the fis- sure angle, the elastic modulus, the crack damage threshold, and the peak strength of brittle sandstone containing three fissures firstly increase and secondly decrease. But the peak axial strain is nonlinearly related to the fissure angle. In the entire process of deformation, the crack initiation and propagation behavior of brittle sandstone containing three fissures under uniaxial compression are investigated with respect to the fissure angle. Six different crack coalescence modes are identified for brittle sandstone containing three fissures under uniaxial compression. The influence of the fissure angle on the length of crack propagation and crack coalescence stress is evaluated. These investigated conclusions are very important for ensuring the stability and safety of rock engineering with intermittent structures.

A state-of-the-art review of mechanical characteristics and cracking processes of pre-cracked rocks under quasi-static compression

The mechanical characteristics and failure behavior of rocks containing flaws or discontinuities have received wide attention in the field of rock mechanics.When external loads are applied to rock materials,stress-induced cracks would initiate and propagate from the flaws,ultimately leading to the irreversible failure of rocks.To investigate the cracking behavior and the effect of flaw geometries on the mechanical properties of rock materials,a series of samples containing one,two and multiple flaws have been widely investigated in the laboratory.In this paper,the experimental results for pre-cracked rocks under quasistatic compression were systematically reviewed.The progressive failure process of intact rocks is briefly described to reveal the background for experiments on samples with flaws.Then,the nondestructive measurement techniques utilized in experiments,such as acoustic emission(AE),X-ray computed tomography(CT),and digital image correlation(DIC),are summarized.The mechanical characteristics of rocks with different flaw geometries and under different loading conditions,including the geometry of pre-existing flaws,flaw filling condition and confining pressure,are discussed.Furthermore,the cracking process is evaluated from the perspective of crack initiation,coalescence,and failure patterns.

Fracture initiation and propagation in intact rock—A review

The initiation and propagation of failure in intact rock are a matter of fundamental importance in rock engineering. At low confining pressures, tensile fracturing initiates in samples at 40%-60% of the uniaxial compressive strength and as loading continues, and these tensile fractures increase in density, ultimately coalescing and leading to strain localization and macro-scale shear failure of the samples. The Griffith theory of brittle failure provides a simplified model and a useful basis for discussion of this process. The Hoek-Brown failure criterion provides an acceptable estimate of the peak strength for shear failure but a cutoff has been added for tensile conditions. However, neither of these criteria adequately explains the progressive coalition of tensile cracks and the final shearing of the specimens at higher confining stresses. Grain-based numerical models, in which the grain size distributions as well as the physical properties of the component grains of the rock are incorporated, have proved to be very useful in studying these more complex fracture processes.

Numerical investigations of the failure mechanism evolution of rock-like disc specimens containing unfilled or filled flaws

The mechanical responses and ultimate failure patterns of rocks are associated with the failure mechanism evolution.In this study,smoothed particle hydrodynamics(SPH)method with the mixed-mode failure model is proposed to probe into failure mechanism evolutions for disc specimens upon loading.The tensile damage model and the Drucker-Prager model are used to calculate the tensile failure and shear failure of the material,respectively.It is concluded that for flaw-unfilled disc specimens,the crack coalescence mechanism in the rock bridge area is affected by the flaw inclination angle and the material property.Considering disc specimens with filled flaws,the incremental rate of tensile damage grows more rapidly when the disc and filling material have a closer ratio of tensile strength to cohesion,which makes the entire specimen response greater brittleness.Furthermore,with the increasing non-uniformity of filling distribution,the incremental rate of tensile-activated damage decreases and the disc specimen performs more ductile.Besides,the influence of the fillings is greater when the flaw inclination angle is approaching 45°.It is proved that the proposed SPH method can be used to simulate the failure mechanism evolution of rocks,which lays a foundation for the study of more complex rock failure.