Numerical modelling of fracture processes in thermal shock weakened rock

This paper presents some preliminary results of a research project aiming at the simulation of thermal shock assisted percussive drilling. In the present study, a numerical model for transient thermal shock induced damage in rock is presented. This model includes a rock mesostructure description accounting for different mineral properties and a thermo-mechanical constitutive model based on embedded discontinuity finite elements. In the numerical simulations, the thermal shock induced damge process is first simulated. Then the uniaxial compression test on thermally affected numerical rock samples is carried out. The effect of thermal shock is demonstrated by comparison to uniaxial compression test simulation on intact rock. The results show that the thermal-shock assisted rock breakage is a feasible idea to be extended to percussive drilling as well.Peer reviewe

Comparative Numerical Study on the Weakening Effects of Microwave Irradiation and Surface Flux Heating Pretreatments in Comminution of Granite

Thermal pretreatments of rock, such as conventional heating and microwave irradiation, have received considerable attention recently as a viable method of improving the energy efficiency of mining processes that involve rock fracturing. This study presents a numerical analysis of the effects of thermal shock and microwave heating on the mechanical properties of hard, granite-like rock. More specifically, the aim is to numerically assess the reduction of uniaxial compressive strength of thermally pretreated specimens compared to intact ones. We also compare the performance of these two pretreatments (conventional heating and microwave irradiation) in terms of consumed energy and induced damage. Rock fracture is modelled by a damage-viscoplasticity model, with separate damage variables in tension and compression. A global solution strategy is developed for solving the thermo-mechanical problem (conventional heating) and the electromagnetic–thermo-mechanical problem (microwave heating). The electromagnetic part of the microwave heating problem is solved in COMSOL Multiphysics software Version 6.1 first. The electromagnetic solution is used as an input for the thermo-mechanical problem, which is finally solved by means of a staggered explicit solution method. Due to the predominance of the external thermal sources, the thermal and the mechanical parts of the problem in both cases are considered as uncoupled. Three-dimensional finite element simulations are utilized to study the damage-viscoplasticity model. An ore-shaped three-mineral numerical rock specimen is used in uniaxial compression tests.publishedVersionPeer reviewe

Modelling the effect of concrete cement composition on its strength and failure behavior

Typical concrete is a mixture of Portland cement, water, and aggregates. While aggregates have a substantial effect on the concrete strength and fracture behavior, the focus of the present study is on the hardened cement paste which can be further divided into the unreacted core, inner and outer products. In high strength concrete, water-to-cement ratio is low, and thus the distance between cement particles is small. Also, the amount of unreacted (high strength) core is higher, and the porosity is low. When water-to-cement ratio is higher, both the distance between cement particles and the porosity due to capillary pores increases. In the present study, we develop a numerical model based on the embedded discontinuity finite elements to predict the effect of the water-to-cement ratio on the compressive fracture behavior of concrete. Representative 2D plane strain simulations demonstrate that the present method captures the major features of concrete fracture and, particularly, qualitatively predicts the known effects of the water-to-cement ratio on concrete compressive strength.publishedVersionPeer reviewe

Numerical modelling of microwave irradiated rock fracture

Rock fracturing through microwave irradiation has received significant attention recently as a viable pretreatment for improving the energy efficiency of comminution processes. This study presents a numerical analysis on the effects of microwave heating on the mechanical properties of hard rock. In particular, the reduction of the uniaxial compressive and tensile strength of granite-like rock due to microwave irradiation induced damage is numerically assessed. Rock fracture is modelled by a damage-viscoplasticity model, with separate damage variables for tension and compression types of failure. A global solution strategy is developed where first the electromagnetic problem is solved in COMSOL multiphysics software, then its solution is used as an input for the thermomechanical problem, which is finally solved by means of a staggered explicit solution method. Due to the preeminence of the thermal radiation, the thermal and the mechanical parts of the problem are considered as uncoupled. The model behaviour is tested in 3D finite element simulations of three-mineral numerical rock specimens, with mesostructures explicitly defined, pretreated first in a microwave oven and then subjected to uniaxial compression and tension tests. The results show that the compressive and tensile strength of rock can be considerably reduced by the microwave irradiation pretreatment.Peer reviewe

A numerical analysis of weakening of a granitic rock by piezoelectric excitation of quartz

This work presents a numerical model to simulate intergranular damage in a granitic rock by oscillating piezoelectric excitation of quartz dispersed in the structure. The damage evolution at grain boundaries was assumed to be related to fatigue and it was modelled using cohesive elements and a damage evolution model formulated in terms of discipation of mechanical work. An explicit representation of the granular mesostructure was built, and it was subjected to high-voltage alternating-current exitation. The effect of the fatigue damage on the mechanical properties was quantified by simulated tension and compression tests. The numerical results show that the electrical treatment can potentially cause rock weakening due to fatigue, but the model needs to be calibrated with experimental data for a quantitative analysis.Peer reviewe

Weakening of tensile strength of granitic rock by HV-HF-AC actuation of piezoelectric properties of Quartz : a 3D numerical study

High-voltage and high-frequency alternating current (HV-HF-AC) excitation of piezoelectric properties of Quartz is a potential method to induce cracks in granite. This was recently shown in a numerical feasibility study [6], where cracking was induced on cylindrical rock samples made of granite by sinusoidal AC excitation at the frequency of ~100 kHz and the amplitude of ~10 kV. However, this study did not investigate the weakening effect due to this cracking on the tensile strength of the sample. The present study addresses this topic numerically. For this end, a numerical method based on 3D embedded discontinuity finite elements for rock fracture and an explicit time stepping scheme to solve the coupled piezoelectro-mechanical problem is adopted. Rock heterogeneity and anisotropy are accounted for at the mineral mesotructure level. A preliminary numerical simulation demonstrates that the HV-HF-AC treatment reduces the tensile strength of a cylindrical granite sample by 12 %, making it thus a potential non-conventional pre-treatment method in comminution and excavation of Quartz bearing rocks and ores.publishedVersionPeer reviewe

Numerical modelling of dynamic spalling test on rock with an emphasis on the influence of pre-existing cracks

This article deals with numerical modeling of rock fracture under dynamic tensileloading and the related prediction of dynamic tensile strength. A special emphasis is laid on theinfluence of pre-existing natural microcrack populations as well as structural (articial) cracks.For this end, a previously developed 3D continuum viscodamage-embedded discontinuity modelis employed in the explicit dynamic nite element simulations of the spalling test. This modelis capable of modelling the eect of natural microcracks populations always present in rocks aswell as to capture the strain rate hardening eect of quasi-brittle materials. In the numericalsimulations of spalling test on Bohus granite, it is shown that the model can predict the pull-pack velocity of the free end of the intact rock sample and the eect of structural cracks witha good accuracy. According to the simulations, the effect of microcrack populations, modeledhere as pre-embedded discontinuity populations, is weaker than the corresponding eect underquasi-static loadin

Numerical modelling of bit-rock interaction in percussive drilling

publishedVersionNon peer reviewe

Combined Viscoplasticity‑Embedded Discontinuity Model for 3D Description of Rock Failure Under Dynamic Loading

This paper presents a combined viscoplasticity-embedded discontinuity model for 3D analyses of rock failure processes under dynamic loading. Capabilities of a rate-dependent embedded discontinuity model, implemented with the linear tetrahedral element, for mode I (tension) loading induced fractures is extended to compressive (shear) failure description by viscoplastic softening model with the Drucker–Prager yield criterion. The return mapping update formulas are derived for the corner plasticity case exploiting the consistency conditions for both models simultaneously. The model performance is demonstrated in 3D numerical simulations of uniaxial tension and compression test on a heterogeneous rock at various loading rates. These simulations corroborate the conception that the rate sensitivity of rock is a genuine material property in tension while structural (inertia) effects play the major role in compression at high loading rates (up to 1000 s−1). Finally, the model is validated with predicting the experiments of dynamic Brazilian disc test on granite.publishedVersionPeer reviewe